江苏盐城市地下含水系统真三维模型.docxVIP

江苏盐城市地下含水系统真三维模型 0 含水层机理评价 在第三纪和第四纪层中，有较强的含水量砂层和相对较弱的粘土层组成了一个完整的含水层系统。以往在对这种系统进行地下水资源评价时,往往将粘性土弱含水层概化为越流层来处理,采用准三维模型对各含水层逐一进行评价计算。这种概化方式不仅忽略了粘性土弱含水层的释水机理,而且大大简化了各含水层之间的水文地质机制,导致了计算结果的重复或缺失。要正确评价此类含水系统的地下水资源量,唯有将各含水层有机地结合起来,将其视为一个统一的复合含水层系统,采用真三维模型来评价。随着地下水动态监测网络的不断完善、计算机科学的飞速发展,为地下水系统的三维模拟计算提供了可能。本文以盐城市为例,详细论述了复合含水层系统地下水资源评价三维数值模型的理论与方法,在建立盐城市复合含水层系统三维数值模型的基础上,结合盐城市对各含水层地下水位的控制要求,评价出了各含水层的最优开采量。 1 地下水开采量划分 盐城市是江苏省苏北沿海的新兴城市,东临黄海,西接泰兴、淮阴,南接南通,北接连云港,面积约1.5×104km2。区内堆积了一套巨厚的新第三系、第四系松散沉积层,构成了巨厚的松散孔隙含水系统,并根据岩性和时代,从上到下分别划分为全新世孔隙潜水含水层,晚更新世第Ⅰ承压含水层,中更新世第Ⅱ承压含水层,早更新世第Ⅲ承压含水层和中、上新世第Ⅳ承压含水层,其间以弱含水的粘性土层相分隔(图1)。由于潜水及第Ⅰ承压含水层中的地下水主要以咸水和半咸水为主,因此,地下水的开采主要以第Ⅱ承压含水层、第Ⅲ承压含水层和第Ⅳ承压含水层为主。自20世纪90年代以来,地下水开采量剧增,地下水位大幅度下降,各含水层均已形成了区域性的地下水位降落漏斗。局部地段由于地下水开采过于集中,已引发了地面沉降和上部咸水入侵等地质环境问题,且严重程度日益增强。因此,正确认识地下水系统的结构与功能,选取合适的数学模型,合理评价出地下水系统的可开采量,对地下水系统的开采进行优化控制,已成为当务之急。 2 地下水流态概化 研究区内有供水意义的含水层为第Ⅱ、第Ⅲ、第Ⅳ承压含水层,各含水层通过其间的粘性土弱含水层发生强烈的水力联系,故将第Ⅱ、第Ⅲ、第Ⅳ承压含水层和其间的粘性土弱含水层作为一个统一的复合含水层系统,并将其概化为三度空间上的非均质各向异性。东部以黄海为边界,南部、北部和西部均以行政区域为界,有水位实测点的地段按第一类边界处理,其余均按第二类边界处理,系统的顶、底部均为隔水边界。整个地下水流态概化为三维非稳定流。区内的抽水井均以乡镇为单位,按大井处理。 3 kxnx+kz.z,t 根据研究区的水文地质概念模型,建立下列与之相适应的数学模型: ??x(kxx?h?x)+??y(kyy?h?y)+??z(kzz?h?z)?W=SS?h?t,(x,y,z)∈Ω,h(x,y,z,t)|t=0=h0(x,y,z),(x,y,z)∈Ω,h(x,y,z,t)|Γ1=h1(x,y,z,t)?(x,y,z)∈Γ1,kxx?h?nx+kyy?h?ny+kzz?h?nz|Γ2=?q(x,y,z,t)?(x,y,z)∈Γ2。??x(kxx?h?x)+??y(kyy?h?y)+??z(kzz?h?z)-W=SS?h?t,(x,y,z)∈Ω,h(x,y,z,t)|t=0=h0(x,y,z),(x,y,z)∈Ω,h(x,y,z,t)|Γ1=h1(x,y,z,t)?(x,y,z)∈Γ1,kxx?h?nx+kyy?h?ny+kzz?h?nz|Γ2=-q(x,y,z,t)?(x,y,z)∈Γ2。 式中:kxx,kyy,kzz分别为各向异性主方向渗透系数(m/d);h为点(x,y,z)在t时刻的水头值(m);W为源汇项(1/d);h0为计算域初始水头值(m);h1为第一类边界的水头值(m);SS为弹性储水率(1/m);t为时间(d);Ω为计算域;Γ1为第一类边界;Γ2为第二类边界:nx、ny、nz为第二类边界x、y、z方向上的外法线方向;q(x,y,z,t)为第二类边界上单位面积的补给量(m/d)。 4 地下水渗透系数及流场分析 上述模型采用有限差分法求解,并采用强隐式法(SIP)联立迭代求解代数方程组。将整个研究区在平面上剖分成100×60的矩形网格单元,在垂向上将其剖分成第Ⅱ、Ⅲ、Ⅳ等3个承压含水层和2个粘性土弱含水层,共5层,每层的有效计算单元为3 076个,共计15 380个(图2)。采用间接法识别、验证模型中的参数(渗透系数k和弹性储水率SS),以2000年12月31日至2002年6月30日作为模型识别的时段,以2002年6月30日至2002年12月31日作为模型验证的时段,将每个季度作为一个地下水开采制度期,共分8个制度期。各含水层均有一定数量的观测井用来进行水位拟合,其中第Ⅱ承压含水层有34个