地下水渗流基本方程及数学模型
地下水动力学(第一章 渗流理论基础-2-专)
∂2H ∂2ψ ∂2H ∂2ψ −K = ; −K =− 2 ∂x∂y ∂y2 ∂y∂x ∂x
二、流网及其性质
流网：在渗流场内，取一组流线和一组等势线 组成的网格。 流网的性质： 流网的性质： 1. 在各向同性介质中，流线与等势线处处垂直， 故流网为正交网格。 证明：等水头线和流线的梯度为：
gradH = ∇H = ∂H ∂H i+ j ∂x ∂y
一般地下水流都为Darcy流。 思考题
§1—3 岩层透水特征分类和渗透系数张量 一、岩层透水特征分类 据岩层透水性随空间坐标的变化情况，将岩层 分为均质的和非均质的两类。 均质岩层：在渗流场中，所有点都具有相同的 渗透系数。 非均质岩层：在渗流场中，不同点具有不同的 渗透系数。 非均质岩层有两种类型：一类透水性是渐变的， 另一类透水性是突变的。 均质、非均质:指 与空间坐标的关系 与空间坐标的关系， 均质、非均质 指K与空间坐标的关系，即不同位 是否相同； 置K是否相同； 是否相同
K1M1 + K2M2 M1 + M2 Kp − Kv = − M1 M2 M1 + M2 + K1 K2 M1M2 = >0 (K1M1 + K2M2 )(M1 + M2 )
(K1 − K2 )
2
渗流模型知识点总结图
渗流模型知识点总结图渗流模型是描述地下水流动和传输的数学模型，它可以帮助我们理解和预测水在地下的流动情况。
渗流模型可以应用于地下水资源管理、地下水污染治理、水文地质等领域，具有重要的实用价值。
下面是关于渗流模型的一些重要知识点总结。
1. 渗流方程渗流模型的数学描述基于渗流方程，它描述了地下水在多孔介质中的流动规律。
渗流方程通常采用达西定律和杜安-卡丁方程进行描述，它们可以用来描述地下水的渗流速度、渗透率、孔隙度等参数之间的关系。
2. 边界条件在渗流模型中，边界条件是描述模型边界上的地下水流动情况的重要参数。
常见的边界条件包括：Dirichlet边界条件、Neumann边界条件和混合边界条件。
这些边界条件可以帮助我们对地下水流动的边界条件进行准确描述，是渗流模型计算的基础。
3. 初始条件渗流模型中的初始条件是指模型开始计算时的地下水流动情况。
初始条件通常是指地下水位和地下水流动速度的初始数值，它们是模型计算的起点。
在模型计算中，初始条件的准确性对计算结果具有重要影响。
4. 离散化方法为了解决渗流方程，通常需要将其离散化。
常见的离散化方法包括有限差分法、有限元法和边界元法等。
这些方法可以将连续的渗流方程转化为离散的问题，通过计算机进行数值计算，得到地下水流动的数值解。
5. 模型验证渗流模型的验证是指利用现场观测数据来验证模型的准确性和可靠性。
验证通常包括比对模型计算结果和现场观测数据，评估模型的拟合程度，以及对模型参数的敏感性分析等。
模型验证可以帮助我们了解模型的适用范围和局限性，提高模型的预测准确性。
6. 模型应用渗流模型在地下水资源管理、地下水污染治理、水文地质和地下水开采等领域有着广泛的应用。
通过渗流模型，我们可以模拟地下水流动过程，预测地下水位和地下水流向，并为地下水资源的合理开发和保护提供科学依据。
此外，渗流模型也可以帮助我们理解地下水污染的传播规律，优化地下水治理方案。
总的来说，渗流模型是描述地下水流动和传输的重要工具，它可以帮助我们理解地下水资源的分布和变化规律，为地下水资源管理和保护提供科学依据。
渗流立方定律
渗流立方定律渗流立方定律是渗流理论中的一项基本定理，也称为达西定律。
它描述了流体通过孔隙介质的速率与孔隙直径的关系。
渗流立方定律的名称源于其方程的形式，即渗流速率与孔隙直径的立方成正比。
在此文中，我们将详细介绍渗流立方定律及其应用。
1. 渗流立方定律的原理和表达式渗流立方定律反映了渗透流动的速度与介质孔隙结构的特征有关，其数学表达式为：Q=kH^3ΔP/μL其中，Q表示单位时间内通过介质的液体（气体）体积，k表示孔隙介质渗透系数，H表示介质厚度，ΔP表示单位长度介质压力差，μ表示介质的动力黏度，L表示介质中液体（气体）通过的距离。
渗流立方定律反映了孔隙介质中液体渗透速率与渗透孔隙的物理结构性质之间的关系。
具体来说，渗透速率随着孔隙直径的增加而增加，并呈现出孔隙直径的立方次幂关系，即Q∝d^3。
因此，该定理也被称为“孔隙方肆立定律”。
渗流立方定律是地下渗透流动理论的基础，广泛应用于水文地质、土壤力学、石油勘探等领域。
下面列举几个应用：(1) 水文地质学。
渗流立方定律可以用于描述地下水的渗透速率。
在地下水资源开发中，可以根据渗流立方定律确定不同孔隙介质的渗透系数，以评估地下水资源的开采潜力和水文地质条件。
(2) 土壤力学。
渗流立方定律可以用于研究土壤中水分的输运规律和渗透特性，对土壤侵蚀、滑坡和沉降等问题有重要意义。
(3) 石油勘探。
渗流立方定律可以用于预测油气藏中的渗透能力和产能。
通过测量油气藏中不同孔隙介质的孔隙直径和自然渗透试验，可以计算得到渗透系数，从而预测油田的产量和石油资源的分布。
地下水渗流基本方程及数学模型总结
常可忽略。
（二）含水层的状态方程
含水层弹性存储的概念： 弹性储存：当地下水水头（水压）降低（或升高）时， 含水层、弱透水层释放（或储存）地下水的性质。 含水层弹性存储的物理意义：
（承压含水层）弹性储存与（潜水）重力储存不同；
第一步：化简方程左端项： 当渗流满足达西定律，且取坐标与各向异性主轴方向一致，有：
H v x K xx x
H v y K yy y
H v z K zz z
( v x ) H H H ( K xx ) [ K xx (K xx )] x x x x x x x
§5 描述地下水运动的数学模型及解算方法
一、各向异性含水层中地下水三维流的基本微 分方程的推导 二、地下水运动微分方程的各种形式 三、地下水运动数学模型的建立及求解
§5 描述地下水运动的数学模型及解算方法
一、各向异性含水层中地下水三维流基本微分方程的推导 为反映含水层地下水运动的普遍规律，研究选定在各向 异性多孔介质中建立地下三维不稳定流动连续性方程。 水均衡的基本思想，对某一研究对象：
描述地下水运动的数学模型及解算方法二地下水运动微分方程的各种形式zzyyxxzzyyxx使潜水面边界处理的简单化直接近似地在微分方程中处理dsdh此时1潜水面比较平缓等水头面呈铅直水流基本水平可忽略渗流速度的垂直分量v2隔水底板水平铅垂剖面上各点的水头都相等各点的水力坡度和渗流速度都相等sin可以近似地用tg代替此即著名的dupuit假设
m d( )
m
1 d d ( )
地下水运动基本定律、基本微分方程和数学模型
第二节 数学模型
导水系数T
当水力坡度为1时，通过整个含水层上 的单位宽度流量。即：
T=K·M
第二节 数学模型
水的状态方程 对于给定质量的水体积，增加一个压力
dPw，水体积产生一定的压缩，根据质量守 恒定律：
ρVw=常数 取全微分有：
ρdVw+Vwdρ=0
由于dPw=dH
第二节 数学模型
达西定律的实质是水流在流动过程中消耗的 能量与流速和渗流长度成正比，与含水层的 渗透系数成反比。
HH1H2
vL K
达西定律的适用范围
Re ud
第一节 达西定律
当雷诺数Re<100时，适用； 当雷诺数Re>100时，不适用； 在天然情况下，绝大多数地下水运动服从达西定律。
第二节 数学模型
地下水渗流连续性方程 表示：在渗流场中的 任何局部，都必须满足质量守恒和能量守恒。
第二节 数学模型
对于稳定渗流，且假定n、ρ不变，则为地
下水稳定流的连续性方程：
x(K x H x) y(K y H y) z(K z H z) 0
第二节 数学模型
形式相似，意义有所差别
x(Tx H x) y(Ty H y)* H t
承压水二维流的微分方程：
第二节 数学模型
当水头变化很小时，即ΔH<0.1h时，对均质 各向同性的潜水有
2xH2 2yH2 T THt
T=Kh h为潜水含水层平均厚度
第二节 数学模型
H(x,y),t0 z,H t0(x,y),z
第二节 数学模型
注：对于稳定流来说，定解条件中没有初始条 件，因为地下水作稳定流时其运动要素是不随 时间而变化的。
2地下水渗流基本方程及数学模型
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
*范围值：n×10-3~ n×10-5； 范围值：0.05~ 0.30。实际测出的值往往小于理论值。
上述两参数之间的不同，还在于潜水含水层存在滞后疏干现象。 弹性释水与重力给水: 对于含水层而言，由于受埋藏条件的限制，抽水时，水的给 出存在着不同。 潜水含水层在抽水过程中，大部分水在重力作用下排出，疏干作用于水位变动带（
为反映含水层地下水运动的普遍规律，我们选定在各向异性多孔介 质中建立地下三维不稳定流动连续性方程。
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
由于渗流场中各点的渗流速度大小、方向都不同，为了反映液体运动的 质量守恒关系，需要在三维空间中建立微分方程形式表达的连续性方程。
则有：
即：
将
代入整理得：
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
所以有
上式为三维流微分方程，也可写成：
物理意义：渗流空间内任一单位体积含水层在单位时间内流入与流出该体 积含水层中的弹性水量的变化量，即单位体积含水层的水量均衡方程。
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
= =
由含水层状态方程，
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
因为 则可得到： 所以有 ，Z为定值，则
于是连续性方程变为：
1渗流基本理论-7
§6 渗流基本微分方程
§6 渗流基本微分方程
2、假设 除与承压含水层基本微分方程有相同假设条件外： （1）当弱透水层的渗透系数K1比主含水层的渗透系数 K小很多时，近似认为水基本上是垂直地通过弱透水 层，折射90º 后在主含水层中基本上是水平流动的。 （如K1与K相差较小时，用等效渗透系数，非越流） （2）主含水层中水头看作是整个含水层厚度上水头的 平均值，即： 1 M H H ( x, y , t ) H ( x, y, z , t )dz M 0 （3）和主含水层释放的水及相邻含水层的越流量相比 ，弱透水层本身释放的水量小到可以忽略不计。
§6 渗流基本微分方程
（2）渗流场中任何一个局部，都必须满足质量守恒和 能量守恒。 4、数学意义 表示渗流空间内任一点任一时刻的渗流规律。 5、讨论 （1）各向同性介质
∂ ∂H ∂ ∂H ∂ ∂H ∂H ( K )+ ( K )+ ( K ) =μ s ∂x ∂x ∂y ∂y ∂z ∂z ∂t
§6 渗流基本微分方程
上次课复习
1、渗流连续性方程—地下水质量守恒定律 （1）表达式
ρ v y ) ∂( ∂( ρ v x ) ∂( ρ vz ) ∂ [ + + ] Δ xΔ yΔ z = ( ρ nΔ xΔ yΔ z ) ∂x ∂y ∂z ∂t
（2）物理含义 某一渗流场中，流入流出单元体的质量差等于单元 体内液体质量的变化。 （3）实质（机理） 水头变化引起含水层弹性释水（贮水）
§6 渗流基本微分方程
3、微分方程的物理意义
∂ ∂H ∂ ∂H ∂ ∂H ∂H ( K xx )+ ( K yy )+ ( K zz ) =μ s ∂x ∂x ∂y ∂y ∂z ∂z ∂t
渗流模型知识点总结
渗流模型知识点总结一、渗流模型概述渗流模型是研究地下水运动及地下水资源管理的一种数学工具。
地下水是地球上水资源的重要组成部分，渗流模型的研究对于有效管理和可持续利用地下水资源具有重要意义。
渗流模型通过数学方法描述地下水在多孔介质中的流动过程，可以预测地下水位、地下水流速、地下水补给和排泄等重要参数，为地下水资源的管理和保护提供科学依据。
二、渗流模型的分类根据模型所涉及的方程和假设的不同，渗流模型可以分为不同的类型。
常见的渗流模型包括：1. 饱和渗流模型：描述地下水在孔隙中完全饱和的情况下的流动规律。
2. 非饱和渗流模型：描述地下水在孔隙中部分饱和或完全不饱和的情况下的流动规律。
3. 二维渗流模型：描述地下水在平面内的流动规律。
4. 三维渗流模型：描述地下水在空间内的流动规律。
根据模型的时间跨度，渗流模型又可以分为：1. 静态渗流模型：描述地下水在静态条件下的分布情况。
2. 动态渗流模型：描述地下水在时间上的变化规律。
渗流模型还可以根据所使用的计算方法不同来进行分类，主要包括有限元模型、有限差分模型、边界元模型等。
三、渗流模型的基本方程1. 边界条件：渗流模型中通常需要给定一定的边界条件，常见的包括恒定头水边界条件、恒定流量边界条件等。
2. 连续方程：描述地下水流线和水位分布的方程，通常为黎曼－莱布尼茨方程。
3. 速度场方程：描述地下水在多孔介质中的流速分布，通常为达西定律或理想渗流方程。
4. 保温方程：描述地下水的运动过程中能量守恒的方程。
5. 变渗透率方程：描述多孔介质中渗透率随深度和位置变化的方程。
以上方程是渗流模型中最基本的方程，通过这些方程可以描述地下水在多孔介质中的流动规律。
四、渗流模型的建立和求解建立一个合适的渗流模型是研究地下水运动的关键。
渗流模型的建立通常需要以下几个步骤：1. 收集地下水数据：包括地下水位、渗透率、孔隙度等信息。
2. 建立地下水模型：通过建立连续方程、速度场方程和边界条件等方程，构建地下水的数学模型。
渗流的基本方程
渗流的基本方程渗流是指在多孔介质中流动的现象，是水文地质学中的重要研究内容之一。
多孔介质是由许多微小的孔隙组成的，例如岩石、土壤、砂土等。
渗流的基本方程描述了多孔介质中流动的物理过程，是渗流理论的核心。
渗流的基本方程可以通过守恒原理来推导，主要包括质量守恒方程和达西-里查德森方程。
质量守恒方程是描述渗流速度分布的方程，它表达了单位时间内通过单位面积的流体质量与孔隙介质中流体储量的变化率之间的关系。
在水平地层中，质量守恒方程可以简化为二维平面问题。
其数学表示为：div（φηρv）= ∂(φηρ)/∂t + div(φηρvq)其中，div表示散度，φ表示孔隙度，η表示介质有效渗透率，ρ表示流体密度，v表示流体速度矢量，q表示流体产生或消失速率。
达西-里查德森方程则是描述渗透压梯度与渗流速度之间的关系。
达西-里查德森方程是根据流体密度不变、黏性流体和渗透性线性增大的假设下推导出来的，经过实验验证，在一定渗透条件下仍然适用。
其数学表示为：v = -K∇h其中，v表示流体速度，K表示渗透性系数，∇h表示渗透压梯度。
通过质量守恒方程和达西-里查德森方程，可以进一步推导得到渗流方程，用于描述多孔介质中任意截面内渗流速度和渗透压梯度之间的关系。
渗流方程可以用一维形式表示为：q = -K∇h其中，q表示单位面积内的流量，K表示有效渗透率，∇h表示渗透压梯度。
渗流方程是多孔介质中流动现象的数学表达式，通过解这个方程，可以求解出多孔介质中的流动速度分布、渗透压梯度分布等有关渗流过程的重要参数。
在实际应用中，渗流方程可以用来预测地下水位变化、估算地下水资源、探测地下水污染传播等。
渗流方程的求解通常依赖于一些边界条件和初值条件。
边界条件是指在孔隙介质的边界上给定的约束条件，如给定流速、压力等。
初值条件是指在求解过程中给定的初始条件。
总之，渗流的基本方程是描述多孔介质中流动现象的数学表达式，包括质量守恒方程和达西-里查德森方程。
渗流方程的原理及应用
渗流方程的原理及应用1. 渗流方程的概述渗流方程是描述岩石或土壤中流体流动行为的基本方程。
它通过描述孔隙介质中流体的质量守恒和动量守恒来推导。
渗流方程可以用于分析地下水运动、油气开采、含水层污染控制等领域。
本文将介绍渗流方程的基本原理以及其在实际应用中的一些典型案例。
2. 渗流方程的基本原理2.1. 渗流方程的基本形式渗流方程的基本形式可以表示为：∂(φS) / ∂t + ∇·(φv) = Q其中，φ是孔隙度，S是饱和度，t是时间，v是流体速度，Q是源项。
渗流方程左侧的第一项表示孔隙度和饱和度的变化率，右侧的第二项表示流体速度的散度，两者相加等于源项的贡献。
2.2. 渗流方程的边界条件渗流方程的边界条件包括初值条件和边界条件。
初值条件是表示渗流方程在初始时刻的状态，边界条件是表示渗流方程在边界上的行为。
例如，对于地下水运动的问题，初值条件通常是地下水位，边界条件可以是水位或流量。
3. 渗流方程的应用案例3.1. 地下水资源管理渗流方程在地下水资源管理中有着广泛的应用。
通过对渗流方程的求解，可以模拟地下水的运动和分布，预测地下水位变化趋势，评估地下水资源的可持续利用能力。
这对于地下水资源的合理规划和管理至关重要。
例如，某地区地下水严重超采，可以通过建立地下水数值模型，解决哪些井应该停用、如何合理分配地下水资源等问题。
3.2. 油气田开发渗流方程在油气田开发中也发挥着重要的作用。
通过对渗流方程的求解，可以模拟油气在储层中的运动和分布，预测油气井的产量和优化开发方案。
这对于提高油气田的开发效率、降低开发成本具有重要意义。
例如，在水驱油开发中，可以通过对渗流方程的模拟，确定注水井的位置和注水量，以实现最优的驱油效果。
3.3. 土壤污染控制渗流方程在土壤污染控制中的应用越发重要。
通过对渗流方程的求解，可以模拟污染物在土壤中的迁移和扩散过程，预测污染物的浓度分布和扩散范围，为土壤污染源治理提供科学依据。
第一章渗流理论基础
绪 论地下水动力学：是研究地下水在孔隙岩石、裂隙岩石和岩溶岩石中运动规律的科学。
它是模拟地下水流基本状态和地下水中溶质运移过程，对地下水从数量上和质量上进行定量评价和合理开发利用，以及兴利防害的理论基础。
第一章 渗流理论基础§1—1 渗流的基本概念一、地下水在含水岩石中的运动1 多孔介质：具有孔隙的岩石。
含水介质一般分为三类：孔隙介质：含有孔隙水的岩层。
裂隙介质：含裂隙水的岩层。
岩溶（Karst ）介质：含岩溶水的岩层。
二、地下水和多孔介质的性质1 地下水的状态方程地下水的状态方程：实际上是地下水的体积和密度随压力变化的方程。
等温条件下，水的压缩系数为：设初始压强p 0时，水的体积为V 0，当压强变到p 时，体积变为V ，由上式得：用Taylor 级数展开，舍去高次项，得到如下的状态方程：V = V 0[1-β（p-p 0）]ρ=ρ0[1-β（p-p 0）]2 多孔介质的某些性质（1）多孔介质的孔隙性孔隙度：指孔隙体积和多孔介质总体积之比。
有效孔隙：互相连通的、不为结合水所占据的那一部分孔隙。
有效孔隙度：指有效孔隙体积和多孔介质总体积之比。
死端孔隙：一端与其它孔隙连通，另一端是封闭的，其中的地下水是相对停滞的。
（2）多孔介质的压缩性天然条件下，一定深度处的多孔介质，要受到上覆岩层荷重的压力。
荷重增加，将引起多孔介质的压缩。
多孔介质的压缩系数：dp dV V 1-=βdpd ρρβ1-=()()000000p p p p p p VV e V V e V V dp VdV ----==-=⎰⎰βββ多孔介质的压缩包括固体颗粒的压缩和孔隙的压缩。
即：V b =V s +V v 上式令上式变为：α=（1-n ）αs +n αp固体骨架的压缩性比孔隙的压缩性小的多，上式变为：α=n αp 三、贮水率和贮水系数1. 水位变化对含水层厚度的影响有效应力地下水位下降，水压力减小，有效应力增大，多孔介质被压缩。
1渗流基本理论-7
§6 渗流基本微分方程
3、 Boussinesq方程—潜水基本微分方程 （1）假设条件 ①符合Dupuit假设； ②忽略水的压缩和骨架的压缩—不符合弹性释（贮）水 规律。原因：潜水面是个自由面，相对压强为0； ③潜水含水层隔水底板水平； ④潜水面存在水量的垂向交换W（ W为潜水面处单位水 平面积、单位时间的入渗量, W> 0 ，入渗；W< 0 ， 蒸发） 。
（6）各向同性柱坐标系（x = rcosθ、y = rsin θ） 1 H 1 2 H 2 H s H (r ) 2 2 2 r r r r K t z 2 H 1 H 1 2 H 2 H s H 或 2 2 2 2 K t r r r r z
导压系数（a）—压力传导系数 描述含水层水头变化的传导速度的参数，其数值等 于含水层的导水系数与贮水系数之比或渗透系数与贮 水率之比。
a=
T
μ
*
=
K
μs
（2）均质各向同性介质
∂2H ∂2H ∂2H μ s ∂H 1 ∂H + 2 + 2 = = 2 K ∂t a ∂t ∂x ∂y ∂z
§6 渗流基本微分方程
§6 渗流基本微分方程
（7）有源（流入）汇（流出）项W或 一般指垂向补给或排泄。 和W分别为三维流和平面二维流的源汇。分别定义 为单位体积含水层和单位水平面积含水层柱体中，单 位时间内产生（为正值）或消耗（为负值）的水量。
∂ ∂H ∂ ∂H ∂ ∂H ∂H ( xx K )+ ( yy K )+ ( zz K )+W =μ s ∂x ∂x ∂y ∂y ∂z ∂z ∂t
流体力学讲义 第十二章 渗流
流体力学讲义第十二章渗流第十二章渗流概述一、概念1.渗流（Seepage Flow）：是指流体在孔隙介质中的流动。
2.地下水流动：在土建工程中，渗流主要是指水在地表以下的土壤和岩石层中的流动，简称为地下水流动。
判断：地下水的流动与明渠流都是具有自由液面的流动。
错二、渗流理论的应用1.生产建设部门；如水利、化工、地质、采掘等部门。
2.土建方面的应用给水方面排灌工程方面水工建筑物建筑施工方面三、渗流问题确定渗流量：如确定通过闸坝地基或井等的渗流流量。
确定渗流浸润线的位置：如确定土坝坝体内的浸润线以及从井中抽水所形成的地下水面线的位置。
确定渗流压力：如确定渗流作用于闸坝底面上的压力。
估计渗流对土壤的破坏作用：计算渗流流速，估计发生渗流破坏的可能性，以便采取防止渗流破坏的措施。
四、土壤的水力特性不均匀系数：（12-1）式中：d60，d10——土壤颗粒经过筛分时分别有60%，10%重的颗粒能通过筛孔直径。
孔隙率n：是指单位总体积中孔隙所占的体积，。
沙质土：n=0.35~0.45；天然粘土、淤泥：n=0.4-0.6。
1.透水性透水性（hydraulic permeability）：是指土或岩石允许水透过本身的性能。
通常用渗透系数k来衡量，k值越大，表示透水性能越强。
均质土壤（homogeneous soil）：是指渗流中在同一方向上各处透水性能都一样的土壤。
非均质土壤（heterogeneous soil）：是指渗流中在同一方向上各处透水性能不一样的土壤。
1各向同性土壤（isotropic soil）：是指各个方向透水性都一样的土壤。
各向异性土壤（anisotropic soil）：是指各个方向透水性不一样的土壤。
2.容水度容水度（storativity）：是指土壤能容纳的最大水体积与土壤总体积之比，数值与土壤孔隙率相等。
3.持水度持水度（retention capacity）：是指在重力作用下仍能保持的水体积与土的总体积之比。
地下水数值模拟
j
i
~ fij
j
j1
Gi
d
,
gij
j1 Gi d , j n
g~ij
j1 Gi d , j n
Ei
D
W T
Gi dxdy
iHi
n j 1
1 j1 j
j1gij
g~ij
Hj
g~ij j gij
H
j 1
j1
fij
~ fij
H n
j
~ fij
j
fij
iHi
n j 1
M j M j1
H
Gi n
Gi
H n
ds
W DT
Gi dxdy
n j 1
j1 j
j1 j1
j
H
j j j1 j
H
Gi
j1
n
j1 j1
j
H n
j j j1 j
H n
j 1
Gi
d
D
W T
Gi dxdy
令：fij
G d , j1
地下水数值模拟
一、基本原理
• 基本思想
——将微分方程得基本解化为边界积分方程, 将边界剖分为有限个单元,在离散得区域边 界上将边界积分方程化为代数方程求解。
• 边界元 ——区域内满足控制方程,边界上近似满足边界条件
• 有限元、有限差 ——区内近似满足控制方程,边界上满足边界条件
一、基本原理
• 特点
u x
v x
u y
v y
dxdy
v
2u x 2
2u y 2
dxdy
v
u n
ds
渗流力学基本方程
渗流力学基本方程渗流力学是研究岩石、土壤等多孔介质中流体运动和物理现象的学科。
渗流力学基本方程是描述多孔介质中流体运动的数学方程组，它是研究渗流问题的基础。
渗流力学基本方程由质量守恒方程和达西定律组成。
质量守恒方程是指在多孔介质中，流体的质量在空间和时间上保持不变。
达西定律是指渗流速度与渗流力的关系，它描述了多孔介质中流体运动的规律。
质量守恒方程是渗流力学中最基本的方程之一。
它可以表示为：∂(ϕρ)/∂t + ∇·(qρ) = S，其中ϕ是多孔介质的孔隙度，ρ是流体的密度，t是时间，q是流体的渗流速度，S是源项。
这个方程描述了多孔介质中流体的密度随时间和空间的变化情况，以及流体的流动和质量的变化。
达西定律是描述多孔介质中渗流速度与渗流力的关系的方程。
它可以表示为：q = -K∇h，其中q是流体的渗流速度，K是多孔介质的渗透性，h是流体的流动势。
这个方程表明，渗流速度与渗流力负梯度成正比，且与多孔介质的渗透性有关。
渗透性越大，流体的渗流速度越大。
渗流力学基本方程还可以通过引入渗透率、渗透率张量等概念，进一步描述多孔介质中流体运动的规律。
渗透率是描述多孔介质对流体渗流的阻力的参数，它与多孔介质的孔隙度、孔隙结构、流体粘度等因素有关。
渗透率张量是描述多孔介质中渗透率随不同方向的变化的张量。
渗流力学基本方程在地下水资源开发、地下水污染治理、石油开采、岩土工程等领域具有重要的应用价值。
通过建立和求解渗流力学基本方程，可以预测多孔介质中流体的运动规律，指导工程设计和实际操作。
渗流力学基本方程是描述多孔介质中流体运动的基本数学方程组。
通过研究和解析这些方程，可以深入理解多孔介质中流体运动的规律，为工程实践提供理论依据和技术支持。
渗流力学基本方程在地下水资源开发、地下水污染治理、石油开采、岩土工程等领域具有广泛的应用前景。
岩土工程渗流：第5章 地下水井流理论
5.1.3 其他井流情 况与相关的约定
向井中注水的情况实际上和抽水的情况一样，在计算与分 析中可以同样利用抽水的相关理论。
本章除特别提到的条件外，一般都用了以下的假定： (1)含水层中水流服从Darcy定律 (2)在水头下降的瞬间水就释放出来（专指潜水浸润面） (3)含水层均质、各向同性、无限延伸 (4)含水层底部水平，承压含水层等厚 (5)抽水前的地下水面是水平 (6)忽略弱透水层的贮水性
17
计算影响半径R的经验公式和经验值
18
19
2）两个观测孔的井流公式
Q 2 KM H0 hw
ln R rw
sw
H0
hw
Q
2 KM
ln R
rw
边界条件并不一定在r=R和r=rw，可以是任何已知点的 水头值（ri，Hi）或降深值（ri，si）。
承压井流量 一般形式为：
Q
2
KM H1
计算机广泛应用 前的时代产物
14
5.2.2 承压井的齐姆（Thiem）公式
1）引用影响半径
影响半径R的取值，理论上没有严格证明，实际上 是困扰人们的一个问题。 该问题分为两类：
无限大区域中R取值； 有界区域R取值
2）两个观测孔的井流公式
由于井壁影响，应用时有误差，工程中常用两个 观测孔的井流公式。
s(x, y,t) H0 (x, y, 0) H (x, y,t)
井中心处降深最大，离井越远降深越小，整个水头下降 区呈漏斗状，称为降落漏斗。
潜水井的降落漏斗在含水层内部扩展，有自由面。
承压水井的降落漏斗没有自由面，是水头的降低区。
如无其他来源，潜水井抽出的水量相当于降落漏斗的含 水层体积的重力疏干，而承压水井的水来自于因降深漏 斗处的水头降低造成含水层的弹性释水。（非稳定井）
地下水渗流耦合力学数值模型
地下水渗流耦合力学数值模型
在地下水渗流耦合力学数值模型中，地下水渗流方程描述了地
下水在多孔介质中的流动过程。
该方程基于达西定律和连续介质力
学原理，考虑了渗透性、孔隙度和渗透率等参数，通过计算流体的
速度和压力分布来描述地下水的运动。
与此同时，围岩力学方程描述了围岩的应力和变形行为。
这些
方程基于弹性力学理论或塑性力学理论，考虑了围岩的弹性模量、
泊松比、强度和变形特性等参数。
通过计算围岩的应力和变形分布，可以了解围岩的稳定性和变形情况。
地下水渗流耦合力学数值模型的基本原理是将地下水渗流方程
和围岩力学方程耦合在一起，形成一个联立的数学模型。
模型通过
离散化方法，如有限元法或有限差分法，将复杂的连续问题转化为
离散的代数方程组。
然后，通过迭代计算的方式，求解这个方程组，得到地下水渗流和围岩的应力和变形场。
地下水渗流耦合力学数值模型在工程领域有广泛的应用。
例如，在地下水资源开发中，可以用于模拟地下水开采对周围围岩的影响，评估地下水资源的可持续利用性。
在地下工程中，可以用于分析地
下水渗流对围岩稳定性的影响，评估工程的安全性。
在地下储气库或储水库设计中，可以用于模拟地下水渗流和围岩变形的过程，优化工程设计。
总之，地下水渗流耦合力学数值模型是一种重要的数值模拟方法，可以帮助我们理解地下水和围岩之间的相互作用，为地下工程和地下水资源管理提供科学依据。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
*范围值：n×10-3~ n×10-5； 范围值：0.05~ 0.30。实际测出的值往往小于理论值。
上述两参数之间的不同，还在于潜水含水层存在滞后疏干现象。 弹性释水与重力给水: 对于含水层而言，由于受埋藏条件的限制，抽水时，水的给 出存在着不同。 潜水含水层在抽水过程中，大部分水在重力作用下排出，疏干作用于水位变动带（
在水位下降为Δ H时，有
即作用于固体骨架上的力增加了H。
。
作用于骨架上力的增加会引起含水层的压缩，而水压力的减少将导 致水的膨胀。
含水层本来就充满了水，骨架的压缩和水的膨胀都会引起水从含水 层中释出，前者就象用手挤压充满了水的海绵会挤出水—样。
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Vv＝nVb；Vs＝(1-n)Vb
式中 ——多孔介质固体颗粒压缩系数，表示多孔介质中固 体颗粒本身的压缩性的指标，s<<p； ——多孔介质中孔隙压缩系数 （Compressibility of the pores of a porous medium），表示多孔介质中孔 隙的压缩性的指标。 n——多孔介质的孔隙度。 因 ，故 。
Ch2 地下水渗流微分方程及数学模型
§2-1 渗流连续方程
一、含水层的状态方程 含水层的状态方程主要包括地下水的状态方程和多孔介质的状态方程。 1、地下水的状态方程 Hooke定律：
式中：E——体积弹性系数（体积弹性模量），20℃时，
E=2.1×105N/cm2。其倒数为压缩系数。 等温条件下，水的压缩系数(coef. of compressibility)为
多孔介质压缩系数的表达式为：
式中，Vb＝Vs+Vv——多孔介质中所取单元体的总体积； Vs——单元体中固体骨架（solid matrix）体积；
Vv——为其中的孔隙（voids）体积。
——介质表面压强；
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
积分（p→p0，V→V0）改写得：
体积：
密度：
按Taylor级数展开，得到近似方程：
和
因
（质量守恒），故有
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
2 、多孔介质的状态方程 多孔介质压缩系数（Coefficient of compressibility）表示多孔介质在 压强变化时的压缩性的指标，用表示。
用于承压水也适用于潜水。对于平面二维非稳定流地下水运动，当研究整
个含水层厚度上的释水情况时，用贮水系数来体现。
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
地下水动力学
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
3、 贮水率和贮水系数 考虑承压含水层受力情况，取一水平横截面AB，按Terzaghi（1883~1963 ）观点：
式中
——上覆荷重引起的总应力（total stress）； ——作用在固体颗粒上的粒间应力 (intergranular stress)； ——横截面面积中颗粒与颗粒接触面积所占的水平面积比； p ——水的压强。 很小，(1- )p ≈ p
贮水系数又称释水系数或储水系数，为含水层水头变化一个单位时，从底
面积为一个单位，高度等于含水层厚度的柱体中所释放（或贮存）的水量 ；指面积为一个单位、厚度为含水层全厚度M的含水层柱体中，当水头改
变一个单位时弹性释放或贮存的水量，无量纲。既适用于承压含水层，也
适用于潜水含水层。 贮水率是描述地下水三维非稳定流或剖面二维流中的水文地质参数，既适
Terzaghi令 ，因此有：
=
称为有效应力（effective stress）。
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
图1—1 一个可压缩的承压含水层（J. Bear）
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
因Vs=constant，故 故 只在垂直方向上有压缩，
上两式表示垂直厚度变化、孔隙度变化与水的压强变化的关 系。 水头降低时含水层释出水的特征，取面积为1m2、厚度为l m (即体积为l m3)的 含水层，考察当水头下降1m时释放的水量。此时，有效应力增加了H＝ g×1=g。
介质压缩体积减少所释放出的水量（dVb）为
Ch2 地下水渗流微分方程及数学模型
思考题 1、在松散层覆盖的岩溶地区，建立高楼群时，出现地表塌陷，如何解释这一现象？
2、一列火车经过一个松散含水层附近的观测孔，该孔水位如何变化？
3、我国不同地区因开地下水引起的地面沉降如何理解？
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
*——贮水系数（storativity）。
*=sM
贮水系数*和贮水率s都是表示含水层弹性释水能力的参数，在地下水 动力学计算中具有重要的意义。
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
贮水率 表示含水层水头变化一个单位时，从单位体积含
水层中，因水体积膨胀（压缩）以及骨架的压缩（或伸长）而释放（或储 存）的弹性水量。单位1/L。
与水体积膨胀所释放出的水量（dV）之和
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
上述二者之和所释放出的水量为
或 式中 s ——贮水率[释水率]（specific storativity），量纲 [L-1]，为弹性 释水[贮水] ；
式中 M——含水层厚度(m)；
饱水带）和包气带两部分，由于包气带的存在，使得饱水带中水的释放存在延滞和滞后 现象。