① 模型建立的方法和步骤
一、模型建立的方法
GMS软件有三种建立确定性模型的方法,包括概念模型法、网格法和Solids法。本书中所选择的方法为Solids法。不管是利用网格法或者概念模型法建模,对含水层结构进行合理的概化是其中一个重要环节,所建模型的准确性很大程度上取决于对实际水文地质条件的正确判断。若轻视对具体水文地质条件的研究,过多依赖模拟技术建立的模型,通常与实际问题相差甚远,也没有使用价值(魏加华等,2003)。当地层出现尖灭、垂向上具有多元结构、水文地质条件比较复杂时,前两种方法不能准确描述此类地层结构,也不能验证基于地质统计学插值求得的含水层顶底板高程是否与实际的钻孔资料相符。GMS中的实体模块Solids利用钻孔资料可以建立地层的三维结构可视化模型,Solids模型定义了地层结构的空间分布,可以切割生成三维显示任意方向的地层剖面(王丽霞等,2011)。
二、模型建立的步骤
利用Solids建模的步骤:
(1)在钻孔模块(borehole)中定义钻孔的坐标位置及垂向上的层位(horizon)。层位即不同地层的交线或岩性分界线。由于地层沉积通常是连续的,因此层位按照一定的次序排列。然而实际地层一般比较复杂,钻孔资料常出现地层缺失现象,遇到此种情况,将缺失的层位空出,使Solids得到的剖面和实际地层剖面相符合。
(2)根据实际的钻孔资料将相应的层位用弧线连接,同时注意地层尖灭的标示。层位连接后生成不同多边形,每个多边形表示相应的地层或岩性。
(3)在地图模块Maps中定义不规则三角网格TIN,来表示地层单元插值的表面边界。
(4)在实体模块Solids选择恰当的插值方法,由horizons生成其相应地层的Solids。如果有N个horizons则有N-1个Solids,Solids生成后即可以在模型上切割任意剖面来检验模型的三维空间结构。
(5)根据Solids数来确定所需网格的最小层数,生成三维网格并进行MODFLOW的初始化。将Solids记录的地层空间信息转成MODFLOW中含水层的顶底板标高,至此地下水三维空间结构模型建立完成。
三、建模过程中可能遇到的问题及解决方法
地下水三维可视化模型建立,首先要基本查明灌区的水文地质条件。了解灌区的地貌、地质条件、构造发育、各地层厚度等信息,需要收集和整理地下水的相关资料,包括灌区水文地质报告、构造图、地质地貌图、水文地质剖面图、电子版地理底图、等高线图、含水层顶底板高程等值线图以及钻孔数据资料等。再结合水文地质条件对含水层资料进行整理和概化。利用GMS建立地下水三维可视化模型时,尤其是在大区域建模中,可能出现3类问题(张永波等,2007;孙红梅等,2008)。
1.由于钻孔分布不均匀而导致的地层缺失
在大区域建模中,由于研究区范围较大,各部分研究程度不同,一般会引起钻孔分布的不均匀。通过不均匀分布的钻孔资料建立水文地质结构模型,可能致使部分地层产生缺失,导致结构模型失真。另外,钻孔分布均匀程度是一个相对概念,对于地形平缓、地层结构相对简单的地区,少量钻孔基本可以比较清楚地反映地层结构;对于地形起伏较大、地层结构比较复杂、构造比较发育的地区,需要较多的有效钻孔,才可能准确揭示地层分布及构造发育状况,然而实际工作中完全实现是不可能的。对于此种问题,根据研究区的地质地貌图、构造分布图及前人绘制的剖面图,对已有的钻孔数据资料进行分析和整理,在具有控制点作用的位置可以适当虚拟部分钻孔数据或者各层面的高程数据,以准确反映该区域地层结构和构造。采用扩充后的钻孔数据资料建立水文地质结构模型,可以弥补由于钻孔资料缺乏而导致的部分地层的缺失。
2.由于钻孔不够深而引起的下伏地层抬升
在钻探工作中,往往有些钻孔深度不够,不能完整地揭露地层。根据这样的钻孔数据建立水文地质结构模型时,系统默认将钻孔底部的标高作为上一层的底部界面。这样就造成下伏地层的抬升。对于这种情况,根据前人绘制的地层等厚度线及剖面图,结合四周钻孔数据对该钻孔资料进行修正,修正后的钻孔资料可以比较准确地反映地层结构。采用修正后的数据资料建立水文地质结构模型,可以有效地控制下伏地层的抬升。
3.由于钻孔资料过细而引起的地层混杂
在野外纪录的钻孔资料中,局部有透镜体形成的地层,透镜体分布的连续性相对较差。采用过细的资料建模,计算机不能分辨透镜体及连续地层,容易出现地层混杂,即将某个钻孔的透镜体地层和另一个或其他几个钻孔的连续地层分界面相连接,导致生成错误的地层结构。对于这种情况,根据该区域剖面图整理资料时,将透镜体区分出来,忽略较小的透镜体,针对较大的透镜体则另外生成地层结构。
此外,在插值计算中,由于计算方法的不同,产生的结果也许会有很大差异,这需要在进行插值计算时,根据不同的具体条件选择适当的插值方法。