《地下水动力学第一讲v讲义教材》由会员分享，可在线阅读，更多相关《地下水动力学第一讲v讲义教材（41页珍藏版）》请在金锄头文库上搜索。
1、1市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用教材一、教材1、地下水动力学薛禹群 主编， 地质出版社，2003年（第二版）2、GMS模型文档文件（数值方法、应用说明文件） （参考书）：3、地下水非稳定流计算和地下水资源评价张蔚榛 主编， 科学出版社，1983年4、地下水运动与资源评价陈雨孙 著，建筑工业出版社， 1986年5、Dynamics of Fluids in Porous MediaBear J. 1972 （图书馆有中译本）6、Introduction to Ground-Water HydraulicsBennett , G. D. 197
2、62市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用内容二、教学内容1）地下水动力学的基本内容（基础知识与理论）以地下水动力学薛禹群主编内容为主2）地下水动力学的数值计算方法（有限差分法应用）以MODFOLLOW-88版说明书为参考教材3）GMS软件应用介绍以GMS应用教学文档（Tutor I， II）为参考教材三、考核形式：平时作业，笔试（开卷）内容：1）基本知识；2）软件应用报告3市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础第一讲 理论基础知识一、地下水动力学研究对象、内容与方法（一）研究对象 （二）研究内
3、容（三）研究方法 （四）地下水动力学研究的里程碑事件（五）地下水动力学的应用领域 （六）有待进行深入研究问题二、渗流理论基础（一）基本参数定义（二）渗流、流网（三）渗流基本方程1、连续性方程2、运动方程（1）承压含水层运动方程（2）越流含水层运动方程（3）潜水含水层运动方程（四）定解条件与数学模型1、定解条件2、数学模型45市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础3、水储量更新时间（数据摘选自“水资源科学与实验研究”，沈振荣等）地球上参加水文循环的水量约5757.7万km3。储水体一次更新年限大气水约8天湖泊水约17年深层地下水约1400年山
4、地冰川约1600年海洋水约2500年极地冰川需万年以上6市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础4、地下水的分类（1）Subsurface Water 与 Groundwater 的区别？（2）按地下水的存在形式分类 包气带（Zone of areation）土壤水Soil Water过渡区Gravitational (Intermediate vadose) Water毛细水Capillary Water饱和带（Zone of saturation）重力水（地下水）Groundwater熔岩含水层（Zone of rock flowage）
5、与其它物质以化学结合形式的水Internal Water潜水面多孔介质含水层7市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础（3）按含水层构造将饱和带分为（如图所示）1）潜水含水层（Unconfined Aquifer）2）承压含水层（Confined Aquifer）3）越流含水层（Leaky Aquifer）5、本门课程的研究对象地下水动力学是研究在重力作用下饱和的多孔介质（的潜水含水层、承压含水层和越流含水层）中水的运动规律的科学。多孔介质：孔隙介质、裂隙介质和溶岩（喀斯特Karst）。多孔介质中水的形态：汽、固、液三态； 其中液态水：吸着水
6、、薄膜水、重力水。作业：地下水需研究的问题包含什么内容？我国水资源开发利用存在的问题与特点？8市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础（二）研究内容针对：1、潜水含水层； 2、承压水含水层； 3、越流含水层在不同的地质特性参数（如均质、非均质，各向同性、各向异性），在不同的流场描述（一维、二维、三维流动），各运动参数随时间的变化特性（恒定、非恒定）等流动问题的研究，得出地下水流动规律、各参数的时空分布，以达到对地下水的开发与管理提供科学的定量化数据。（三）研究方法1、解析法：特点（常微分方程、偏微分方程，时间变量）2、实验法（模型试验，现场观
7、测试验，电模拟）：特点3、数值模拟：特点9市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础（四）地下水动力学研究的里程碑事件1、Darcy 定律（Henry Darcy 1856年）2、潜水井流问题的裘布依方程（Dupuit 1863年）3、承压井的非稳定流方程及其解析解（1）O. E. Meinzer，1928年开始观注地下水的非恒定流与承压含水层的储水性质；（2）C. V. Theis（泰斯）1935年给出承压含水层非恒定流动的泰斯公式4、1960年左右，利用计算机进行地下水方程的数值分析具有达标性的分析软件系统（1）Feflow原东德的WASY
8、公司，1978发布第一版，数值方法：有限元；（2）MODFLOW美国地质勘察局，1988年发布第一版，数值方法：有限差分法；（3）GMS（Groundwater Model System）以Modflow以及其他地下水问题的分析软件为内核，添加上前处理与后处理模块构成地下水问题的综合软件分析系统（本教学所用为：GMS6.0版）10市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础（五）地下水动力学的应用领域1、工程建设方面：水利、地质、石油、建筑等2、水资源评价、开发利用与管理方面3、水环境分析（六）有待进行深入研究问题1、介质方面：多相流、裂隙介质、
9、溶岩流动；2、机制方面：溶质与热量在地下水中的运动与运移机制；3、方法研究：数值解法，随机理论，并行计算；4、优化管理：地下水开采许可（井群分布、抽水量、价格导向等）。作业：上网查找有关地下水动力学的发展与工程应用资料，写出简要报告（500字以内）。11市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础二、渗流理论基础（一）基本参数定义1、与地下水贮存相关的参数（1）孔隙度（porosity) n有效孔隙度（effective porosity）ne有效孔隙度：相通的，不为结合水占据的孔隙度。死端孔隙（见教材p.4，图1-1）： 所贮水量不参与地下水运
10、动，但参与抽水或灌水运动。12市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础（2）给水度（Specific yield）y也称有效孔隙度，在面积为积为 A的柱体中，当柱体潜水面下降一个h时时所释释放出的水的体积积V，或对对三相图图，指给给定体积积Vb的饱饱和含水体所能释释放（贮贮存）的水的体积积，即：给给水度反映了土壤的给给水性，即饱饱和土壤在重力作用下能自由排出水量的性能。给给水度的影响因素：对对面积为积为 A的柱体，设设其潜水面下降h，排出水体积积量是t 的函数V（t），可定义义 1）完全给给水度y：2）瞬时给时给 水度t：3）平均给给水度t：
11、13市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础（3）含水率与田间持水量0（specific retention, field capacity)1）含水率：在非饱和带中的土壤在重力和表面张力作用下典型单元体中所保持的水量Vw和土壤总体积Vb之比2）饱和度（Saturation）Sw：3）田间持水量0 ：是由于附着力与内聚力作用土壤所吸附的水分。4）给水度与田间持水量及孔隙度的关系（见图） 14市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础（4）贮水率与贮水系数1）水体压缩系数（Water compr
12、essibility）2）多孔介质质中的固体颗颗粒压压缩缩系数（solid compressibility）s3）多孔介质质中孔隙压缩压缩 系数率（Porous compressibility）p4）多孔介质压缩质压缩 系数（Aquifer compressibility）15市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础5）贮贮水率（Specific storage）s单单位：m3/(m3m)=1/m ；量纲纲：L-1；含义义：体积为积为 1立方米的承压压含水层层在水头头下降1米时时所释释放的水体体积积量。其中 在水头头不降到承压压含水层层隔水顶
13、顶板以下时时，含水层层只能引起含水层层的弹弹性释释水，一般认为弹认为弹 性释释水在整个含水层层内瞬时时完成。对对潜水含水层层，当水头头下降时时可引起两部分排水。含水层层上部为为重力排水，用给给水度y表示重力排水的能力；深层饱层饱 和水部分引起弹弹性释释水，用贮贮水率（释释水率）s表示。 s = 10-3 10-5 ； 而 y = 0.05 0.3； 在潜水层层，一般忽略不计弹计弹 性释释水量。4）贮贮水系数（Storage Coefficient）* 对对二维维承压压含水层层，有贮贮水系数反映整个含水层层厚度的释释水能力i） 表示水头下降1米由水体膨胀所释放的水体体积；ii） 表示水头下降1米
14、由含水层压缩所挤出的水体体积。16市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础2、与地下水运动有关的参数（1）Darcy 定律与渗透系数K（Hydraulic Conductivity）1856年法国工程师Henry Darcy 在装满砂的圆筒中进行实验，得出如下结果： 其中 K：渗透系数；单位：m/s， 一般用：m/d。 定义水力坡度： Darcy定律推广为三维情况：17市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础1）对渗透系数的讨论K：综合反映了多孔介质的渗透特性，其与多孔介质的颗粒结构（大小组
15、 成、空间分布）和渗透流体的性质（粘性、密度）有关。按颗粒结构的不同，可将K分为：均质与非均质；各向同性与各向异性。非均质：K（x，y，z）；各向异性：在同一特征单元体，沿不同轴向，K值不同；均匀性与各向异性概念之区别。为表示沿任意面上的渗透系数，用二阶张量表示K： Darcy定律的推广： 对各向同性含水层： 渗透系数矩阵为：18市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础2）层状岩层的等效渗透系数i）岩层与水流流向平行分布的等效Kpii）岩层与水流方向垂直分布的等效Kv总有： Kp Kv。19市政系水资源与水工研究所马长明 地下水动力学讲稿地下
16、水动力学理论、数值技术与软件应用理论基础3）Darcy 定律的适用范围雷诺数的定义：其中 d：土壤的平均颗粒直径。地下水流态层流：紊流：Darcy 定律适用范围：见教材p. 15说明及图1-13。20市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础（2）渗透率（intrinsic permeability）k在渗透系数K中消除流体性质不同的影响，仅考虑颗粒结构特性，即：其中 k称为渗透率，仅与多孔介质的粒径有关。系数c与颗粒形状有关。k的量纲：L2；单位常用：cm2及 darcy。1 darcy = 9.87 X 10 -9 cm2 。（3）导水系数（Transmissivity）T对二维均质含水层，其厚度为M，则定义导水系数：T = KM 单位：cm2/d； 量纲：L2T-1导水系数T反映了通过整个含水层的输水能力，表示水头降为1米下整个含水层的单位宽度的流量值。21市政系水资源与水工研究所马长明 地下水动力学讲稿地下水动力学理论、数值技术与软件应用理论基础（4）越流系数（coefficient of leakage）与 越流因素（leak
