专利名称：一种含水介质模型的多尺度随机耦合建模方法
技术领域：
本发明涉及水文地质、地质领域，更具体地说是一种建立具有多尺度性质的非均质含水介质或者岩性介质模型的方法。
背景技术：
建立含水介质模型对于分析研究区水文地质条件、建立水文地质概念模型和获取水文地质参数等具有重要的意义。转移概率地统计方法(TPROGS, TransitionProbabilityGeostatistics),使用基于马尔科夫链性质的含水介质材质类型转移概率,替代指示变异函数或协方差实现指示地统计，从而根据钻孔岩性数据建立含水介质三维空间模型。含水介质及其参数具有非均质性和尺度特性，其往往表现为大尺度上的结构性和小尺度上的随机性特征。传统上，使用转移概率地统计方法，是基于研究区全部钻孔含水介质材质类型分布数据基础上，统计转移概率从而进行模拟的，其模拟结果往往只能体现单一尺度上的含水介质分布信息，难以体现水文地质的结构性和随机性耦合特征；而且，在自然界中不同区域的含水介质材质类型分布规律并不相同，即含水介质材质类型的空间变异性是一种非平稳的随机过程，而传统方法使用所有钻孔数据计算转移概率，其无法描述不同水文地质结构体之间的含水介质材质类型分布特征的差异性。而且这种问题不仅在转移概率地统计建模时存在，对于序贯指示模拟(Sequentialindicatorsimulation)和指示克里金方法(IndicatorKriging)等，同样存在类似问题。因此，迫切需要提出一种新的处理方法以解决含水介质建模时存在的上述问题。
发明内容
本发明是为避免上述现有技术所存在的不足之处，提供一种利用钻孔岩性数据的含水介质模型的多尺度随机耦合建模方法，以便使模拟结果即可以体现含水介质的结构性，又可以体现含水介质的随机性耦合特征。本发明为解决技术问题采用如下技术方案:本发明含水介质模型的多尺度随机耦合建模方法，其特点是按如下步骤进行:(I)建立父尺度模型:以被测区域为中心划定扩展区域，所述扩展区域包含所述被测区域；将所述扩展区域中的钻孔岩性数据按岩性的含水性质分类为m类，分别标记为父尺度材质类别i，i=l，2，…，m，利用所述父尺度材质类别i使用转移概率地统计方法进行模拟，建立初始尺度的含水介质三维网格模型，以所述初始尺度的含水介质三维网格模型作为父尺度模型；(2)建立父尺度模型中所有材质类别分别对应的子尺度模型:(a)以被测区域为子尺度模型的子尺度模拟范围，将所述子尺度模拟范围内的钻孔岩性数据按岩性的含水性质分为n类，标记为子尺度材质类别j，j=l,2, -, n, n ^ m ；
(b)统计所述子尺度模拟范围内的父尺度材质类别i所对应的子尺度材质类别j，使用转移概率地统计方法进行模拟，建立子尺度模拟范围内的与父尺度材质类别i所对应的各子尺度模型i ;所述子尺度模型i在空间直角坐标系中X、Y、Z三个方向的网格大小分别为父尺度模型相同方向的网格大小的整数分之一倍；所述子尺度模型i的位置和体积按如下方法设置:若所述子尺度模拟范围的边界位于父尺度模型的网格线上，则所述子尺度模型i的位置与所述子尺度模拟范围的位置相同，所述子尺度模型i的体积等于所述子尺度模拟范围的体积；若所述子尺度模拟范围的边界不位于父尺度模型的网格线上，则将子尺度模拟范围向外扩大到子尺度模拟范围的边界位于父尺度模型的网格线上，得外扩子尺度模拟范围，所述子尺度模型i的位置与所述外扩子尺度模拟范围的位置相同，所述子尺度模型i的体积等于所述外扩子尺度模拟范围的体积；(3)父尺度模型与子尺度模型的融合:将子尺度模拟范围内的父尺度材质类别i所在区域用子尺度模型i的相同区域替换，得到被测区域含水介质模型。与已有技术相比，本发明有益效果体现在:1、本发明使用具有结构性特征的父尺度模型作为先验知识，约束提取各材质类别对应的子尺度建模的数据，然后建立子尺度的随机模型，最后再耦合两个尺度的模型，使得建立的含水介质模型兼具有父尺度的结构性和子尺度的随机性；2、本发明在建立各子尺度模型时，所使用的建模数据是父尺度模型对应材质类型控制的子尺度的钻孔材质类型数 据，因此，子尺度模型计算的转移概率更符合含水介质材质类型的空间分布规律；3、本发明在父尺度和子尺度两个尺度上分别对钻孔岩性数据进行概化，使得不同尺度使用不同详细程度的含水介质分类数据，从而使建立的含水介质模型具有相应尺度的概化特性。
图1为本发明实施例1中m=2时建立的父尺度模型示意图；图2为本发明实施例1中n=4、子尺度模拟范围为整个父尺度模型的区域时，父尺度模型与子尺度模型的融合示意图。
具体实施例方式实施例1:本实例以美国密歇根州渥太华县辖区内主体部分的矩形区域为被测区域。本实施例中含水介质模型的多尺度随机耦合建模方法，按如下步骤进行:(I)建立父尺度模型:以所述被测区域为中心划定扩展区域，扩展区域包含被测区域。本实施例中被测区域与扩展区域相同。将扩展区域中的钻孔岩性数据按照表I的父尺度钻孔岩性与含水介质材质类别对照关系分类为I含水材质、II隔水材质两类，分别标记为父尺度材质类别i，i=l，2，所述父尺度材质类别i数据格式为:ID、X、Y、Z、H、材质类别，同一钻孔垂向上不同材质段ID相同，X和Y表示材质段的平面坐标，Z表示材质段顶部高程，H表示材质段长度，材质类别用I或II表示。利用所述父尺度材质类别i使用转移概率地统计方法进行模拟，建立初始尺度的含水介质三维网格模型，以所述初始尺度的含水介质三维网格模型作为父尺度模型；模拟结果如图1所示，模型的参数如下:X方向:最小值:477599(m)，网格大小:682(m)，网格个数:60(个)；Y方向:最小值:245725(m)，网格大小:682(m)，网格个数:59(个)；Z方向:最小值:130.37(m)，网格大小:3.29(m)，网格个数:40(个)；转移概率地统计是基于马尔科夫链性质进行材质转移概率统计，然后使用序贯指示模拟，并使用退火模拟优化的一套方法，该方法具体原理参照(Carle，S.F.(1999).T-PROGS:TransitionProbabiIityGeostatisticalSoftware(Version2.1), HydrologicSciencesGraduateGroup, UniversityofCalifornia, Davis.)。表1父尺度钻孔岩性与含水介质材质类别对照关系
权利要求
1.一种含水介质模型的多尺度随机耦合建模方法，其特征是按如下步骤进行: (1)建立父尺度模型: 以被测区域为中心划定扩展区域，所述扩展区域包含所述被测区域； 将所述扩展区域中的钻孔岩性数据按岩性的含水性质分类为m类，分别标记为父尺度材质类别i，i=l，2，…，m，利用所述父尺度材质类别i使用转移概率地统计方法进行模拟，建立初始尺度的含水介质三维网格模型，以所述初始尺度的含水介质三维网格模型作为父尺度模型； (2)建立父尺度模型中所有材质类别分别对应的子尺度模型: Ca)以被测区域为子尺度模型的子尺度模拟范围，将所述子尺度模拟范围内的钻孔岩性数据按岩性的含水性质分为n类，标记为子尺度材质类别j，j=l,2, -, n, n ^ m ； (b)统计所述子尺度模拟范围内的父尺度材质类别i所对应的子尺度材质类别j，使用转移概率地统计方法进行模拟，建立子尺度模拟范围内的与父尺度材质类别i所对应的各子尺度模型i ;所述子尺度模型i在空间直角坐标系中X、Y、Z三个方向的网格大小分别为父尺度模型相同方向的网格大小的整数分之一倍；所述子尺度模型i的位置和体积按如下方法设置: 若所述子尺度模拟范围的边界位于父尺度模型的网格线上，则所述子尺度模型i的位置与所述子尺度模拟范围的位置相同，所述子尺度模型i的体积等于所述子尺度模拟范围的体积； 若所述子尺度模拟范围的边界不位于父尺度模型的网格线上，则将子尺度模拟范围向外扩大到子尺度模拟范围的边界位于父尺度模型的网格线上，得外扩子尺度模拟范围，所述子尺度模型i的位置与所述外扩子尺度模拟范围的位置相同，所述子尺度模型i的体积等于所述外扩子尺度模拟范围的体积； (3)父尺度模型与子尺度模型的融合: 将子尺度模拟范围内的父尺度材质类别i所在区域用子尺度模型i的相同区域替换，得到被测区域含水介质模型。
全文摘要
本发明公开了一种含水介质模型的多尺度随机耦合建模方法，其特征是以被测区域为中心划定扩展区域，扩展区域包含被测区域；将扩展区域内的钻孔岩性数据按岩性的含水性质进行分类，使用转移概率地统计方法建立父尺度模型，然后建立父尺度模型中被测区域内所有材质类别分别对应的子尺度模型，最后将子尺度模型与父尺度模型融合。本发明使用父尺度模型的材质空间分布作为先验知识，约束选取父尺度模型中各材质类别分别对应的子尺度的建模数据，使得建立子尺度模型的数据是使用父尺度同一材质控制的子尺度的数据，建模数据统计特性更符合该材质的空间分布规律，且使得建立的耦合含水介质模型兼具有父尺度的结构性和子尺度的随机性优点。
文档编号G06F17/50GK103246783SQ20131019523
公开日2013年8月14日 申请日期2013年5月23日 优先权日2013年5月23日
发明者马雷, 钱家忠, 赵卫东 申请人:合肥工业大学