第七章 补给子程序包概念及程序化

使用补给子程序包(RCH）的目的是模拟地下水系统的面状补给。面状补给的通常是由降水入渗补给地下水系统形成。本模型中的补给定义为:

其中，是水平面上某一计算单元（i，j）的补给率，用单位时间水流体积量表示。

是施加在该计算单元面积上的补给通量（用单位时间长度单位表示)。

补给率施加到位于（i，j）的柱体上某一单个计算单元。

没必要让补给发生在同一垂向柱体的不同深度上，这是因为自然补给多是从顶部进入地下水系统的。

最简单的情形是地下水系统的上界面出现在第一层内；可是，当潜水面升高或降低时，地下水系统的上界面在各点的位置将随时间而变化。

稍后将要叙述如何通过MODFLOW规定的三种选择，对接受补给的每个垂向柱体指定接受补给的计算单元。

补给(RCH)子程序包可用于模拟降水以外的其它水源补给，如人工补给。

若需要在同一垂向柱体上不止一个计算单元接受补给，那么可用井流子程序包。

井流子程序包允许对模型中任一计算单元进行补给或排泄。

补给子程序包将每个应力期（除非选择了使用前一应力期的补给通量）的补给通量读入到一个二维数组中，然后将补给通量的值乘以计算单元在水平面上的面积，求得的值，并将其储存在RECH数组中。

如何选择垂向柱体上接受补给的计算单元，可通过补给方式选择符和补给层号数组来指定。

补给方式选择项包括:

⑴补给指定到模拟层第一层;

⑵补给指定到垂向柱体上任一计算单元，并由二维数组,中的层号确定;

⑶假若计算单元上方没有定水头计算单元，补给指定到位于最上面的有效计算单元。

对于第一种以及第二种选择，如果指定接受补给的计算单元是不透水的，实际上没有接受补给。

至于第三种选择，因为已经假定补给将被定水头水源截获，如果垂向柱体上有定水头计算单元以及其上方有无效计算单元，那么该垂向柱体实际上没有采纳补给。

模型读取的补给通量，必须采用与其它模型参数相一致的长度和时间单位，如米/天。
 

给定水平面上计算单元位置(i，j）以及由补给选择项所确定的垂向位置k，相应的补给量，将会在矩阵方程求解时，从(公式(26)或(29))中减去。每次迭代都对所有接受补给的计算单元进行此项计算。按原有定义，补给是不受含水层水头影响的，所以水头系数不添加任何值。

在实际工作中，应仔细考虑所研究的问题及上述三种选择，再决定选用哪一种选择。例如，图38所示情形，用一个横剖面模型来模拟补给，河流向地下水渗流，以及地下水渗流进入河流（图38-a)。按照第五章描述的潜水水平方向水力传导系数计算的规定，模型中有效计算单元在模拟过程中被逐渐调整以近似表达地下水面，结果每个垂向柱体中最上方的有效计算单元大致落在潜水面上。这个过程给出了有效、定水头、以及无效计算单元的最终分布情况，见图38-b。
 

图38-c说明采用上述第一种选择时，模型的补给分布情况。该选择只允许对模型最顶层补给。当潜水面低于第一层之底面标高时，含水层将无法获得应有的补给。这样做显然会低估整个潜水面所获得的补给量。
 

图38-d说明采用第二种选择时的补给分布。这里假定用户在模拟之前已根据预先估计的潜水面的位置而指定了接受补给的计算单元。但由于预先估计的潜水面与模拟所得的最终结果有些出入，用户指定为接受补给的计算单元中，有四个已转换成无效计算单元，因而无法接受补给。
图38-e是采用第三种选择时的模拟结果。这种选择最适合这种特定情形。在这项选择下，除定水头计算单元用来表示河流外，补给进入每个垂向柱体最上面的有效计算单元。从而模拟了补给条件下的潜水面形状。
 

就典型的降水补给来说，第三种选择最为常用。用户不需判定哪一个是垂向柱体中最上面的有效计算单元。因为程序具备了在模拟过程中对此作自动判断的功能。不过，在第一层中有不透水计算单元的地方，补给不应该穿透至下伏层位时，也用得着第一种选择，否则MODFLOW会将补给加至不透水层之下的含水层之中。当然，第三种选择仍然可用于这种情形，只要将不可透水的计算单元的补给率指定为零即可。用户应选用最方便指定输人数据的选择项。与此类似，当第一层以外的层位露出地表，同时指定层位的补给不应该穿透不透水计算单元进人下伏含水层，或许用得着第二种选择。在衡量所作选择是否最佳时，应考虑的其它因素是，比起第一、三种选择来，第二种选择占用更多的计算机存储空间,而第三种选择比第一、二项选择需用较多计算时间。
 

图38.补给子程序包中三种选择下接受补给的计算单元。