边界条件 源汇 等特殊问题的模拟与处理-陈崇希.pptVIP

对于沙岸滨海地区，一般地说,只有潜水层 (当海水位高于潜水层底面时) 才能以海岸线 近似作为该模拟层的边界，其下诸承压含水层的海边界位置要视其向海底的延伸等条件而定（图4-14-2）。 Nawang等［1990］在马来西亚某滨海地区建立地下水模型时,他们认识到陆区承压含水层应向海底延伸，但又不知其延伸的距离，不得已，只好假定承压含水层的顶板延伸5km和15km两种情况的海底边界来建模。美国夏威夷某滨海地区建立地下水模型，也无依据地向海底延伸８英里。 如何确定含水系统的海边界问题成为滨海地区地下水模型的关键。 4.14.2 承压含水层的海边界确定方法 采用海上勘探的方法探明海底含水系统的岩性结构与海底地形,在技术上并无困难，但要耗费大量资金和很长的工期。能否探索一种理论上站得住脚而又便于应用的方法呢? 陈崇希等（1986~1988）在北海市滨海区工作时认为承压含水层顶板应向海底延伸，提出等效（排泄）边界的概念，并利用潮汐波动效应确定边界和含水层参数的方法(陈崇希等,1989a，1990a，Chen等,1990,1992)，我们将其概括为两种方法。 方法一：承压含水潮汐效应的模拟 人们早已了解到，滨海区承压含水层观测孔中的水位随海潮而有观律地波动（可称为潮汐效应 ），这是普遍存在且非常重要的信息。我们的思路是：不同的海底边界条件应该有不同的波动观律，是否可以利用这些信息来近似确定承压含水层的海底边界条件及水文地质参数呢? 滨海区承压含水层观测孔水位的潮汐效应，除海潮外还受什么因素控制呢？显然，另一个重要因素是承压含水层顶板向海向延伸的长度，在其末端出露海水处（可视为承压水的“露头”） 直接与承压含水层发生联系；再者是海水通过顶板的越流与地下水的联系；其三是海潮荷载作用的潮汐效应(Todd,1959;陈崇希,1975,1983)。因此，就准三维模型而言，考虑上述三个因素，地下水流控制方程可表述为： 图4-14-1 潮汐效应机理示意图 该定解问题与无海边界问题的差别是：①水流的控制方程（4-14-1a）左端多了刻画海水越流及海潮荷载作用的第3、4项。其中：Hsea为海水位；TE为潮汐效率(陈崇希,1981,21页;1983,18页;2011,38页］。这两项仅在海域承压含水层顶板分布范围Dsea内起作用。②方程（4-14-1d）是承压含水层海底露头的第一类边界条件。 上述定解问题表述的是，在潮汐作用和边界作用下（若模拟区内无抽/注水、无降雨及无蒸发作用时）承压水的水头动态。为了刻画潮汐单独作用引起的水头动态，可将陆区Dland的边界取至足够远处（通常海岸往里3~4km内外即可），以至对承压水的波动区无影响。或者仅模拟几个海潮波动周期te（一个海潮波动周期约25小时)，这时陆区Dland边界条件的变化对承压水的波动区的影响与海潮相比，通常可忽略不计。另外，承压水水头波动区的初始水头分布是很难测得的，对此可以置H0=0，只要模拟几个海潮波动周期就可获得稳定的波动水头分布。 如果以波动值s取代水头值H（包括海潮）并以平均水头为零.，则海潮波动引起承压含水层水头波动的问题可改写为： 这就是滨海区承压含水层潮汐效应定解问题。如果该含水系统还包括潜水层，只要在陆区加上弱透水层与潜水层即可；其它情况类推。该方程组可用来模型识别求取有关参数，包括承压含水层顶板的海底分布。同样，该定解问题稍作改变，可用于考虑蜜度变化的海水入侵的溶质运移模型(成建梅等,1999) 方法二：“等效（排泄）边界” 1986~1988年，陈崇希等在北海滨海水源地进行地下水资源评价中，提出“等效（排泄）边界”概念及确定方法(陈崇希等,1989a,1990a;Chen等,1990,1992)。 上述影响滨海区承压水波动的三要素中，只有承压含水层海底露头以及海水通过顶板的越流两者与地下水发生水力联系，而海潮荷载作用的潮汐效应并不发生水量交换，它是“假水力联系”（陈崇希,1975,22页; 1981,21页;1983,18页; 陈崇希等2011,38页），且对承压含水层水头仅起振荡作用，其平均值为零。另外，就数值模型而言，全区的模型识别和预测，通常要数十年或更长的时间，若去拟合承压水波动动态，则要求以波动的海水位作为其边界条件，这种情况的时间步长（通常不大于半小时）不能放大，为此大大地增加了运算量，以致一般难以承担，实际上通常这么做也是没有必要的。