MODFLOW中文原理说明

___Y1

已于 2024-11-09 17:11:05 修改

阅读量618

点赞数 5

文章标签：算法人工智能机器学习

于 2024-11-09 14:41:47 首次发布

本文链接：https://blog.csdn.net/weixin_44259522/article/details/143644742

版权

译者的话

摘要

第一章概述

第一节目的

第二节本书内容的安排

第三节致谢

第二章多孔介质中地下水流动的有限差分公式数学模型

第一节离散化的概念

第二节有限差分公式

第三节迭代求解

第四节差分方程的求解形式

第五节计算单元的类型和边界条件的处理

第六节有关垂向离散的讨论

第三章 MODFLOW程序设计

第一节总体结构

第二节模型边界和含水层边界

第三节水均衡计算

第四节内存分配

第五节三维数组下标的定义

第六节数据输人

第七节结果输出控制

第八节主程序

第四章基本子程序包(BAS)

第一节概念及程序化

第二节子程序包的选择及其输人文件

第三节边界数组(IBOUND)

第四节初始条件

第五节时间离散

第六节输出

第七节水均衡计算

第八节基本子程序包输入数据及格式

第九节输入数据说明

第十节 BAS子程序包输入样单

第十一节输出控制部分的输入数据

第十二节输人数据说明

第十三节输人控制选择输入样单

第五章计算单元间渗流子程序包(BCF)

第一节概念及程序化

第二节水力传导系数的计算

第三节承压含水层中水平方向上水力传导系数的计算

第四节潜水含水层中水平方向上的水力传导系数的计算

第五节垂向水力传导密的计算

第六节非饱和带垂向渗流的计算

第七节贮水系数的处理

第八节贮水系数的转换

第九节特殊处理方式的使用及其局限性

第十节输入数据

第十一节 BCF子程序包输入数据及格式

第十二节输入数据说明

第十三节 BCF子程序包输人样单

第六章河流子程序包(RIV)

第一节概念及程序化

第二节河流子程序包输入数据及格式

第三节输人数据说明

第四节河流子程序包输入样单

第七童补给子程序包(RCH)

第一节概念及程序化

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

其中， $Q_{Ri,j}$ 是水平面上某一计算单元(i,j)的补给率，用单位时间水流体积量表示。 $I_{i,j}$ 是施加在该计算单元面积 $DELR_j*DELC_i$ 上的补给通量(用单位时间长度单位表示)。

补给率 $Q_{Ri,j}$ 施加到位于(i,j)的柱体上某一单个计算单元。没必要让补给发生在同一垂向柱体的不同深度上，这是因为自然补给多是从顶部进人地下水系统的。

最简单的情形是地下水系统的上界面出现在第一层内；可是，当潜水面升高或降低时，地下水系统的上界面在各点的位置将随时间而变化。稍后将要叙述如何通过MODFLOW规定的三种选择，对接受补给的每个垂向柱体指定接受补给的计算单元。补给(RCH)子程序包可用于模拟降水以外的
其它水源补给，如人工补给。若需要在同一垂向柱体上不止一个计算单元接受补给，那么可用井流子程序包。

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

补给子程序包将每个应力期(除非选择了使用前一应力期的补给通量)的补给通量
1.读入到一个二维数组RECH,中，然后将补给通量的值乘以计算单元在水平面上的面积
DELR*DELC,，求得Q的值，并将其储存在RECH数组中。如何选择垂向柱体上接受补
给的计算单元，可通过补给方式选择符NRCHOP和补给层号数组IRCH来指定。补给方式
选择项包括:(1)补给指定到模拟层第一层;(2)补给指定到垂向柱体上任一计算单元，并由
二维数组IRCH,中的层号确定;(3)假若计算单元上方没有定水头计算单元，补给指定到位
于最上面的有效计算单元。对于第一种以及第二种选择，如果指定接受补给的计算单元是
不透水的，实际上没有接受补给。至于第三种选择，因为已经假定补给将被定水头水源截
获，如果垂向柱体上有定水头计算单元以及其上方有无效计算单元，那么该垂向柱体实际
上没有采纳补给。模型读取的补给通量，必须采用与其它模型参数相一致的长度和时间单
位，如米/天。

给定水平面上计算单元位置(i，j)以及由补给选择项所确定的垂向位置k，相应的补给量Q.，将会在矩阵方程求解时，从RHS*(公式(26)或(29))中减去。每次迭代都对所
有接受补给的计算单元进行此项计算。按原有定义，补给是不受含水层水头影响的，所以
水头系数HCOF.*不添加任何值。
在实际工作中，应仔细考虑所研究的问题及上述三种选择，再决定选用哪一种选择。
例如，图38所示情形，用一个横剖面模型来模拟补给，河流向地下水渗流，以及地下水渗
流进入河流(图38-a)。按照第五章描述的潜水水平方向水力传导系数计算的规定，型
中有效计算单元在模拟过程中被逐渐调整以近似表达地下水面，结果每个垂向柱体中最上
方的有效计算单元大致落在潜水面上。这个过程给出了有效、定水头、以及无效计算单元
的最终分布情况，见图38-b.
图38-c说明采用上述第一种选择时，模型的补给分布情况。该选择只允许对模型最顶
层补给。当潜水面低于第一层之底面标高时，含水层将无法获得应有的补给。这样做显然
会低估整个潜水面所获得的补给量。
图38-d说明采用第二种选择时的补给分布。这里假定用户在模拟之前已根据预先估计
的潜水面的位置而指定了接受补给的计算单元。但由于预先估计的潜水面与模拟所得的最
终结果有些出人，用户指定为接受补给的计算单元中，有四个已转换成无效计算单元，因
而无法接受补给。
图38-e是采用第三种选择时的模拟结果。这种选择最适合这种特定情形。在这项选择
下，除定水头计算单元用来表示河流外，补给进人每个垂向柱体最上面的有效计算单元。
从而模拟了补给条件下的潜水面形状
就典型的降水补给来说，第三种选择最为常用。用户不需判定哪一个是垂向柱体中最
上面的有效计算单元。因为程序具备了在模拟过程中对此作自动判断的功能。不过，在第
一层中有不透水计算单元的地方，补给不应该穿透至下伏层位时，也用得着第一种选择，
否则MODFLOW会将补给加至不透水层之下的含水层之中。当然，第三种选择仍然可用于
这种情形，只要将不可透水的计算单元的补给率指定为零即可。用户应选用最方便指定输
人数据的选择项。与此类似，当第一层以外的层位露出地表，同时指定层位的补给不应该
穿透不透水计算单元进人下伏含水层，或许用得着第二种选择。在衡量所作选择是否最佳
时，应考虑的其它因素是，比起第一、三种选择来，第二种选择占用更多的计算机存储空间，而第三种选择比第一、二项选择需用较多计算时间。