各向异性渗流固结方程
- 格式：pdf
- 大小：421.03 KB
- 文档页数：5
文档推荐
-
各向异性渗流固结方程页数:5
-
062三维各向异性稳态渗流页数:7
-
渗流的基本定律(达西定律)页数:38
-
061二维各向异性稳态渗流页数:8
下载文档原格式
合集下载
地下水渗流基本方程及数学模型总结
p减少 有效应力增大会引起固体颗粒的压缩变形。固 体颗粒的压缩变形比多孔介质中空隙的变形小得多，通
常可忽略。
（二）含水层的状态方程
含水层弹性存储的概念： 弹性储存：当地下水水头（水压）降低（或升高）时， 含水层、弱透水层释放（或储存）地下水的性质。 含水层弹性存储的物理意义：
（承压含水层）弹性储存与（潜水）重力储存不同；
第一步：化简方程左端项： 当渗流满足达西定律，且取坐标与各向异性主轴方向一致，有：
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 ( )
常可忽略。
（二）含水层的状态方程
含水层弹性存储的概念： 弹性储存：当地下水水头（水压）降低（或升高）时， 含水层、弱透水层释放（或储存）地下水的性质。 含水层弹性存储的物理意义：
（承压含水层）弹性储存与（潜水）重力储存不同；
第一步：化简方程左端项： 当渗流满足达西定律，且取坐标与各向异性主轴方向一致，有：
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 ( )
第一章_渗流相关方程(2)
g H H p g 1 p t t t
（1-70）
令 s ( n ) g 则 根据连续性原理有：
M H s x y z t t
p
1 dV v Vv d
2 贮水率和贮水系数 考虑承压含水层受力情况，取一水平横截面AB，按 Terzaghi（1883~1963）观点： 式中 ——上覆荷重引起的总应力（total stress）； s ——作用在固体颗粒上的粒间应力 (intergranular stress)； ——横截面面积中颗粒与颗粒接触面积所占的水平面 积比； p——水的压强。 Terzaghi令 s = ' 称为有效应力（effective stress）。 很小，(1- )p ≈ p，因此有：
Vv＝nVb；Vs＝(1-n)Vb n dV v 1 n dV s 1 dV s 1 dV v Vb d Vb d Vs d Vv d
(1-58)
Hale Waihona Puke (1 n) s n p
1 dVs
(1-59)
式中 s Vs d ——多孔介质固体颗粒压缩系数，表示多 孔介质中固体颗粒本身的压缩性的指标，s<<p； ——多孔介质中孔隙压缩系数 （Compressibility of the pores of a porous medium），表 示多孔介质中孔隙的压缩性的指标。 n——多孔介质的孔隙度。 因 (1 n)s ，故 n p 。
dn 因Vs=constant，故dV b Vb 1 n
dV b d (z ) Vb dz
只在垂直方向上有压缩，
故 d ( z ) z dp
（1-62） （1-63）
（1-70）
令 s ( n ) g 则 根据连续性原理有：
M H s x y z t t
p
1 dV v Vv d
2 贮水率和贮水系数 考虑承压含水层受力情况，取一水平横截面AB，按 Terzaghi（1883~1963）观点： 式中 ——上覆荷重引起的总应力（total stress）； s ——作用在固体颗粒上的粒间应力 (intergranular stress)； ——横截面面积中颗粒与颗粒接触面积所占的水平面 积比； p——水的压强。 Terzaghi令 s = ' 称为有效应力（effective stress）。 很小，(1- )p ≈ p，因此有：
Vv＝nVb；Vs＝(1-n)Vb n dV v 1 n dV s 1 dV s 1 dV v Vb d Vb d Vs d Vv d
(1-58)
Hale Waihona Puke (1 n) s n p
1 dVs
(1-59)
式中 s Vs d ——多孔介质固体颗粒压缩系数，表示多 孔介质中固体颗粒本身的压缩性的指标，s<<p； ——多孔介质中孔隙压缩系数 （Compressibility of the pores of a porous medium），表 示多孔介质中孔隙的压缩性的指标。 n——多孔介质的孔隙度。 因 (1 n)s ，故 n p 。
dn 因Vs=constant，故dV b Vb 1 n
dV b d (z ) Vb dz
只在垂直方向上有压缩，
故 d ( z ) z dp
（1-62） （1-63）
水利堤坝工程中渗透参数的选取及渗流计算方法评价
水利堤坝工程中渗透参数的选取及渗流计算方法评价水利堤坝工程中渗透参数的选取及渗流计算方法评价摘要：渗流是引起涉水工程破坏的重要原因，因此渗流计算是水利水电工程涉水工程设计中不可或缺的步骤。
渗透参数的选取与渗流方法的选择，直接影响对工程渗流稳定性的评价。
本文结合笔者多年工作经验，就水利水电工程设计中渗透参数的选取与渗流计算的几种方法进行了初步的分析，并总结出渗流计算注意的一些问题，提高了计算结果准确性，对进一步采取防渗措施提供参考。
关键词：水利工程渗流计算堤坝设计引言堤防工程的设计与施工准则要求保证堤防建筑物能抵御洪水的威胁。
由于堤防大多沿天然河岸修建，因此，堤防基础的渗透稳定问题普遍存在。
本文主要针对堤防渗流参数的选用并对渗流计算方法进行了评价。
1、渗流计算目的（1）坝体（堤身）浸润线的位置。
（2）渗透压力、水力坡降和流速。
（3）通过坝体（堤身）或堤基的渗流量。
（4）坝体（堤身）整体和局部渗流稳定性分析。
2、计算工况及渗透系数的选用岩土工程参数的选用需要根据满足给定保证率时，通过实验方法选用。
不同工况需要选用不同的参数，否则就无法满足工程设计所需要的保证率。
2.1常规堤防工程常规的堤防工程计算提出了三种水位组合，此三种水位组合的渗流计算目的及相应土体的渗透系数选取原则主要为：（1）临水侧为高水位，背水坡为相应水位。
本组合的计算目的：①计算背水坡可能最高的逸出点位置、背水坡逸出段及背水坡基础表面出逸比降，用于背水坡渗流安全复核、反滤层及排水设施设计；②背水坡面可能最高的浸润线，用于背水边坡稳定计算；③当堤身、堤基土的渗透系数大于10-3cm∕s时，计算渗流量，用于分析防渗措施对本工程运行要求的可行性和背水坡排水设施设计（对于大坝均要求进行渗流量计算）。
对上述第①、②种计算目的工况，堤身、堤基的渗透系数则取小值平均值，对第③种计算目的工况则取大值平均值。
（2）临水侧为高水位，背水坡为低水位或无水。
本组合的计算目的：①背水坡面可能最高的浸润线，用于背水坡边坡稳定计算，相应各土体的渗透系数取小值平均值；②复核局部渗流稳定及进行反滤层设计，则进行局部渗流稳定性复核土体的渗透系数取小值，其上、下部位土体的渗透系数取大值平均值。
2地下水渗流基本方程及数学模型
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
*范围值：n×10-3~ n×10-5； 范围值：0.05~ 0.30。实际测出的值往往小于理论值。
上述两参数之间的不同，还在于潜水含水层存在滞后疏干现象。 弹性释水与重力给水: 对于含水层而言，由于受埋藏条件的限制，抽水时，水的给 出存在着不同。 潜水含水层在抽水过程中，大部分水在重力作用下排出，疏干作用于水位变动带（
为反映含水层地下水运动的普遍规律，我们选定在各向异性多孔介 质中建立地下三维不稳定流动连续性方程。
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
由于渗流场中各点的渗流速度大小、方向都不同，为了反映液体运动的 质量守恒关系，需要在三维空间中建立微分方程形式表达的连续性方程。
则有：
即：
将
代入整理得：
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
所以有
上式为三维流微分方程，也可写成：
物理意义：渗流空间内任一单位体积含水层在单位时间内流入与流出该体 积含水层中的弹性水量的变化量，即单位体积含水层的水量均衡方程。
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
= =
由含水层状态方程，
地下水动力学
安徽理工大学 地球与环境学院 水资源与规划系
Ch2 地下水渗流微分方程及数学模型
因为 则可得到： 所以有 ，Z为定值，则
于是连续性方程变为：
FLAC3D渗流说明书(中文版)
FLAC3D流固耦合（手册翻译）1.1简介FLAC3D通过具有渗流性的实体（比如土）来模拟流体的流动。
流动模型的建立可以独立于力学计算而自动完成，或者说可以与力学模型同时建立，这样就可以考虑流体与土体之间的相互作用。
流固耦合的一种类型是“固结”，即:空隙水压力逐渐消散而导致土体的沉降。
这个过程包括两种力学反映：一，空隙水压的改变导致有效应力的变化，这将影响到土体的力学反映（如：有效应力的减小可能导致塑性区的产生）；二，力学实体中某一区域的流动会随着空隙水压的改变而改变。
该程序可以计算完全饱和情况下的流动，也可以模拟具有自由水面的流动。
模拟具有自由水面的流动时，自由水面以上的部分空隙水压等于0，气相将不参与计算。
对于不考虑毛细水压力颗粒较粗的材料可以采用这种模拟方法。
流体计算就有以下特点：1 根据各项同性和各项异性的渗流计算，相应采用两种流体运动定律。
流动中的null材料用来模拟流动范围内的非渗流材料。
2 不同区域可以拥有不同的流动模型（isotropic, anisotropic or null）和模型参数。
3 可以事先指定流体的压力、流量、非渗流区边界条件。
4 流体源可以以电源，也可以以体源的形式插入到材料中，这些源对应于流体的流入或流出，可以随着时间而变化。
5 对于完全饱和流动，可以采用显式和隐式两种算法，但对于非饱和流动则只能采用显示计算。
6 任何力学和温度计算模型都可以与流体模型一起使用，在耦合计算中，可以考虑饱和体的压缩性和热膨胀性。
7.流体与力学计算的耦合通过提供比奥系数来实现。
和不排水温度系数β8．与温度的耦合计算可以通过提供线性热膨胀系数αt（undrained thermal coefficient,可能翻译的不对）来实现。
9.热-流体计算以线性理论为基础，假定材料参数为常数，不考虑对流。
流体与实体的温度保持局部平衡。
非线性行为可以采用FISH语言改变孔隙压力、材料特性来实现。
土力学地基基础课件第三章渗流固结理论
渗流固结理论的重要性
渗流固结理论在土木工程、水利工程 、地质工程等领域具有广泛的应用价 值。
它对于理解土体的力学行为、预测土 体的变形和稳定性、优化工程设计和 施工具有重要意义。
渗流固结理论的应用领域
01
02
03
水利工程
水库、堤防、水电站等水 利设施的设计和安全评估。
土木工程
高层建筑、高速公路、桥 梁等基础设施的建设和安 全评估。
渗透试验
通过测量土体的渗透系数、 渗透速度等参数，研究土 体的渗透特性。
现场试验方法
现场观测
通过在土体中埋设传感器和监测 仪器，实时监测土体的渗流和固
结过程。
触探试验
通过触探设备对土体进行触探，测 量土体的物理性质和强度特性。
旁压试验
通过旁压设备对土体施加压力，测 量土体的变形和强度特性。
数值模拟方法
三维固结理论通过求解偏微分方程组， 得到土体在固结过程中任意时刻的孔隙
水压力分布、土层沉降和位移场。
04
渗流固结理论的实验研究
室内试验方法
室内模型试验
通过模拟实际土体中的渗 流和固结过程，研究土体 的变形和强度特性。
土工离心机试验
利用离心加速度模拟土体 应力状态，研究土体在复 杂应力状态下的渗流和固 结行为。
06
结论
渗流固结理论的发展趋势
数值模拟与实验研究的结 合
随着计算机技术的进步，数值 模拟方法在渗流固结理论的研 究中越来越受到重视。通过与 实验研究相结合，可以更准确 地模拟复杂条件下的土体渗流 和固结过程。
多场耦合分析
考虑土体的应力、应变、渗流 和温度等多场耦合效应，对土 体的复杂行为进行更全面的分 析。
渗流固结理论可以用于分析地 下水的流动规律和土体的渗透 性能，为地下水控制提供理论 支持。
第一章 渗流理论基础-2-专
2Z , 2 x 2Z 2Z 2Z 2 0 2 2 y 则方程 x y 为Laplace方程。
H H vx K ; vy K x y
由达西定律知：
v ; v y 流函数满足的条件： x y x H H 有 K ; K x y y x
(2)平行层面的等效渗透系数总是大于垂直 层面的等效渗透系数。 证明：（以二层为例）
K1M 1 K 2 M 2 M 1 M 2 K p Kv M1 M 2 M1 M 2 K1 K 2 M 1M 2 0 K1M 1 K 2 M 2 M 1 M 2
K1 K 2
n
取等效渗透系数Kv，那么单宽流量为：
二式相等得：
因此，
Mq q n M i K vb b i 1 K i
M Kv n Mi i 1 K i
Kv
M
i 1 n
n
i
Mi i 1 K i
此式为层状岩层垂直于层面的等效渗透系数。 说明：(1)当某一层的Ki较小时，Mi/Ki较大，Kv变小； 当Ki→0时， Mi/Ki→∞，Kv→0，也就是说，垂直于层 面的等效渗透系数主要取决于渗透系数最小的分层。
1.2.3非线性运动方程 Re小于1—10时，地下水流为线性流，用 Darcy定律描述； Re大于1—10时，地下水 流为非线性流，用下列定律描述： Forchheimer公式： 1901年福希海默提出Re>10时： J=av+bv2 Chezy公式 1912年克拉斯诺波里斯基提出紊流公式： 1
vK J2
dq=vxac+vybc 因，ac=dy，bc=-dx 所以，dq= vxdy-vydx 把 v ; v
x
H H vx K ; vy K x y
由达西定律知：
v ; v y 流函数满足的条件： x y x H H 有 K ; K x y y x
(2)平行层面的等效渗透系数总是大于垂直 层面的等效渗透系数。 证明：（以二层为例）
K1M 1 K 2 M 2 M 1 M 2 K p Kv M1 M 2 M1 M 2 K1 K 2 M 1M 2 0 K1M 1 K 2 M 2 M 1 M 2
K1 K 2
n
取等效渗透系数Kv，那么单宽流量为：
二式相等得：
因此，
Mq q n M i K vb b i 1 K i
M Kv n Mi i 1 K i
Kv
M
i 1 n
n
i
Mi i 1 K i
此式为层状岩层垂直于层面的等效渗透系数。 说明：(1)当某一层的Ki较小时，Mi/Ki较大，Kv变小； 当Ki→0时， Mi/Ki→∞，Kv→0，也就是说，垂直于层 面的等效渗透系数主要取决于渗透系数最小的分层。
1.2.3非线性运动方程 Re小于1—10时，地下水流为线性流，用 Darcy定律描述； Re大于1—10时，地下水 流为非线性流，用下列定律描述： Forchheimer公式： 1901年福希海默提出Re>10时： J=av+bv2 Chezy公式 1912年克拉斯诺波里斯基提出紊流公式： 1
vK J2
dq=vxac+vybc 因，ac=dy，bc=-dx 所以，dq= vxdy-vydx 把 v ; v
x
刘月田各向异性油田
刘月田各向异性油藏渗流理论与开发方法1各向异性油藏特点？注水开发时如何设计调整？渗透率具有方向性的油藏叫做各向异性油藏。
各向异性油藏有两大类。
一类是裂缝作用造成的，称为裂缝各向异性油藏：另一类是沉积作用形成的，称为沉积各向异性油藏。
（1）特点：渗流速度方向一般情况下不与压力梯度保持一致：渗流速度大小随压力梯度大小和方向两者发生变化。
渗流速度的大小都会改变，两者方向一般情况下不平行；当且仅当压力梯度方向在渗透率主轴上时，渗流速度的方向和压力梯度方向平行。
各向异性汕藏中，渗流速度的大小和方向由位势梯度的大小和方向共同决泄，无论位势梯度的大小还是方向发生变化，渗流速度的大小都会改变；当且仅岂位势梯度方向在渗透率主轴上时，渗流速度的方向和位势梯度方向平行。
由于各向异性渗透率对井网具有破坏与重组作用，会明显改变原来的井网形式，所以当井网中同一注采单元内任意两口井的连线与各向异性渗透率主方向平行，各向异性油藏变换为等价各向同性油藏时，井网注采单元不会被破坏，只是形状发生变化。
一般情况下，各向异性油藏布井方法如下：1.井排方向与渗透率主方向平行或垂直。
井排方向指同一注采单元内任意两井连线；渗透率主方向指裂缝方向或沉积过程中的古水流方向。
2.各向异性油藏井网设计的计算公式为：苴中，a和d分别为各向异性油藏设计井网的井距和排距，a，和d份别为等价各向同性油藏井网的井距和排距，Kx和Ky分别是各向异性渗透率的最大和最小主值。
各向异性油藏在各方向的总体导流能力等价于渗透率为SK冬的各向同性油藏。
并向异性油藏井网的开发效果可以用等价各向同性油藏井网来表示，而各向同性汕藏井网设计及其开发效果分析技术早已成熟。
2各向异性油藏水平井特点及设计方法？跟直井井网相比复杂程度成倍增加：（1）水平井网需考虑渗透率主方向、井排方向和水平井段方向三者之间的两两匹配关系，直井只有渗透率主方向和井排方向的关系；（2）水平井网需考虑井距、排距和水平井段长度的两两匹配关系，直井只有井距和排距的关系。
材料的各向异性对渗流场的影响分析
材料的各向异性对渗流场的影响分析胡瑞;曾胜【摘要】在土石坝设计中,材料渗透参数的选取是坝体防渗设计的重要部分,它对坝的渗透稳定起着至关重要的作用.考虑非饱和区的影响,在材料各向异性的情况下,分别对其进行渗流计算分析,计算结果表明:在非饱和区的存在和材料渗透参数的各向异性的条件下,单宽流量明显增大;非稳定渗流下,自由面变化显著.【期刊名称】《三峡大学学报（自然科学版）》【年(卷),期】2010(032)001【总页数】3页(P17-19)【关键词】土石坝;自由面;渗流场;非饱和区【作者】胡瑞;曾胜【作者单位】湖南五凌电力工程有限公司,长沙,415000;湖南五凌电力工程有限公司,长沙,415000【正文语种】中文【中图分类】TV223.4在土石坝的设计中,材料渗透参数的选取是坝体的防渗设计的重要部分,它对坝的渗透稳定起着至关重要的作用[1].同时,在水利工程建设中,各向异性渗流问题愈来愈受到关注,如在各向异性材料的土坝内会带来坝体内水头压力、浸润线抬高、在坝下游的出渗坡降增大等现象,相对于各向同性渗流场而言其危害性更大[2-4].利用数值模拟方法计算堤坝工程的非饱和区渗流,并在材料各向异性的基础上进一步研究非饱和区渗流作用下堤坝自由面变化情况,为堤坝非饱和区渗流在工程上的应用提供参考.1 非饱和土的渗流方程在饱和土壤中,引起水分转移的力是重力和水的压力,在非饱和土中,支配着土壤水在液态下整体转移的是重力和水的表面张力,Richards等曾在1931年就证明非饱和土中的渗流与饱和土一样符合达西定律和连续方程.若将达西定律代入连续方程(忽略渗透过程中总应力的改变和土颗粒骨架的变形)并以总水头h作为未知量,当渗透的主方向与坐标轴一致时,非饱和土渗流的二维微分方程就可表示为式中,kx,ky分别为x和y方向上的渗透系数;θw为体积含水率;h为总水头;t为时间. 令y为位置水头,则h=uw/γw+y.若mw为土水特征曲线的斜率,则∂θw=mw∂uw =mwγw∂(h-y),式(1)就可以写为由于y是常数,式(2)简化为当计算稳定渗流时,则从式(4)也可以看出非饱和土中的渗透系数k不再是常数,而是含水量的一个函数,此函数被称为非饱和土的渗透系数函数.因此,要进行渗流分析“需具备以下基本条件: (1)材料的渗透系数函数;(2)确定边界条件.当上述条件已知时,就可以对以上方程进行求解,从而得到非饱和土渗流场.2 土坝饱和-非饱和渗流分析2.1 计算模型[5]某均质土坝,坝高50 m,顶宽10 m,上游坝坡为1∶2.5,下游坝坡为1∶2,上游水位为45 m,下游水位为5m,渗透系数为0.29m/d,孔隙率为0.464.大坝的模型如图1所示,分2091个节点,2000个四边形单元.图1 土坝计算模型图2.2 计算工况的选取考虑非饱和区材料渗透系数各向异性,分别以不同的水库水位下降速度来进行渗流计算分析,其工况组合见表1～2.表1 工况计算表工况编号材料各向异性下降速度/(m◦d-1) 1 K x=Ky V=0.5 2 K x=Ky V=0.1 3 K x=2Ky V=0.5 4 K x=2Ky V=0.1 5 K x=3Ky V=0.5 6 K x=3Ky V=0.1表2 非稳定渗流计算时步表时步时间增量/d V=0.5 V=0.1消逝时间/d V=0.5V=0.1 1 8 40 8 40 2 8 40 16 80 3 8 40 24 120 4 8 40 32 160 5 8 40 40 200 6 8 40 48 240 7 8 40 56 280 8 40 64 320 9 8 40 72 360 10 8 40 80 400 11 8 40 8 440 12 8 80 16 480 13 8 120 24 520 14 8 160 32 560 15 8 200 40 600 16 8 240 48 640 17 8 280 56 680 18 8 320 64 720 19 8 360 72 760 20 8 400 80 800 82.2.1 稳定渗流计算结果分析由于非饱和区的存在,浸润线附近的饱和区的渗流水会穿过浸润线而进入非饱和区;增加水平方向的渗透系数,更多的水流向下游,负压区逐渐减小,等势线向下游推移,如图2～4所示.从表3可以看出,随着水平渗透系数的增加,更多的水向下游渗流,逸出点略有抬升,意味着浸润线抬升,渗流量明显增大,说明材料的各向异性对其影响较大.表3 稳定渗流计算结果对照表工况单宽渗流量Q /(m3◦d-1)逸出点Y坐标/m逸出点坡降K x=Ky 2.0834 21.25 0.365 K x=2Ky 3.9877 22.5 0.364 K x=3Ky 5.8283 23.75 0.3552.2.2 非稳定渗流计算结果图5～10可以看出(图中的数字为时步数),库水位下降浸润线也下降,当不考虑材料的各向异性时,水位迅速下降,坝内的孔隙水压力来不急消散,自由面的形状较陡,同时在坝的上游面形成严重的“倒流”现象,这将影响上游坝坡稳定(如图5所示).在同一下降速度时,随着水平渗透系数的增加,自由面的形状由陡逐渐趋于平缓,变化较为显著,说明坝内的孔隙水压力消散迅速,下降形成的水“倒流”现象不严重.尤其是工况5和工况6时,由于水平渗透系数较大,在水库下降速度较快(V=0.5m/d)时,其坝内自由面变价也较为平缓,对上游坝坡的稳定性影响不大,说明材料的各向异性对渗流场的影响较为显著.3 结语由于非饱和区的存在,流水会穿过浸润线而进入非饱和区,坝内非饱和区域也存在渗流.在渗透参数的各向异性的条件下,随着水平渗透系数的增加,渗透坡降变化不大,而单宽流量增大较明显;在非稳定渗流中,即使水库水位快速下降,坝内的孔隙水压力也能迅速消散,自由面变化也较显著,其对上游坝坡稳定产生的影响比不考虑材料各向异性时的要小,说明渗透参数的各向异性对其影响较大.单宽流量的增大,在保证渗透稳定的条件下可根据兴利的要求来选取排水设施,但它也会给下游坝坡带来影响;水库降落时,孔隙水压力的消散迅速,它同时也给上游坝坡带来安全隐患,对坝体的渗透稳定也同样存在着影响,故都应做好防渗排水措施.而实际工程中,存在着土体非饱和性、不均匀性、材料各向异性的特点,所以对材料渗透参数的各向异性在非饱和条件下来进行渗流分析较为合理、更为贴近工程实际.参考文献:[1] 韩瑜,李松德.土石坝心墙内浸润线的逸出高度[J].西安理工大学学报,2000(2):204-206.[2] 毛昶熙.渗流计算分析与计算[M].北京:水利电力出版社,1990:62-68.[3] 张家发.土坝饱和与非饱和稳定渗流场的有限元分析[J].长江科学院院报,1994(9):41-45.[4] 刘晓庆.高土石坝稳定-非稳定渗流有限元分析[J].三峡大学学报:自然科学版,2008,30(3):30-33.[5] 王晓章.非稳定渗流及其作用下的边坡稳定数值模拟[D].西安:西安理工大学,2004:62-68.。
渗流方程的原理和应用
渗流方程的原理和应用1. 渗流方程的基本概念渗流方程（Darcy’s law）是描述岩石或土壤中流体渗透运动的基本方程。
它是流体力学中的一种基本方程，由法国工程师亨利·达西（Henry Darcy）在1856年提出。
渗流方程表示了渗透流量与渗透率之间的关系。
渗透率是描述岩石或土壤中孔隙的互连程度的物理量，它决定了流体在岩石或土壤中的移动能力。
达西理论适用于岩石、土壤、沙砾等多孔介质中的流体渗透问题。
2. 渗流方程的数学表达根据渗流方程，渗透流量（Q）等于渗透率（K）乘以梯度（∇h）。
数学表达式如下：Q = -K * ∇h其中，Q表示渗透流量，K表示渗透率，∇h表示压力梯度。
3. 渗流方程的应用3.1 地下水资源评估渗流方程在地下水资源评估中起着重要作用。
通过对地下水流动的模拟和预测，可以评估地下水资源的分布、储量和可利用性。
利用渗流方程可以计算地下水的流量和流速，并研究不同参数对地下水流动的影响。
3.2 污染物迁移研究渗流方程在研究污染物在地下水中的迁移、扩散和传输过程中也得到了广泛应用。
通过模拟污染物在地下水中的迁移行为，可以评估污染物对地下水质量的影响，指导环境保护和水资源管理。
3.3 石油开采渗流方程在石油开采领域的应用也非常重要。
通过研究岩石的渗透性和岩石中原油、天然气等流体的运移规律，可以指导石油开采工程的设计和操作。
渗流方程在石油开采领域的应用可以优化采油方案，提高油田开采效率。
3.4 地下工程渗流方程在地下工程中的应用也很广泛。
地下工程包括地下建筑、隧道、地下储气库等，渗流方程可以用于模拟地下水的流动，评估地下工程的稳定性和可行性。
4. 渗流方程的局限性渗流方程是基于一些假设和简化条件推导出来的，因此在某些情况下可能存在局限性。
例如，渗流方程假设介质是均质、各向同性的，但实际介质往往是非均质和各向异性的。
在研究介质非均质性和各向异性时，需要引入更复杂的模型和方法。
此外，渗流方程还假设流体是层流流动，不考虑湍流效应。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
with the symmetry property aij= a+ The partial components of this tensor do not have the conventional significance. If we consider a cube of unit size of the bulk material, u represents the total normal tension force applied to the fluid part of the faces of the cube. Denoting by p the hydrostatic pressure of the fluid in the pores we may write
M. A. BIOT* Shell Development Company, New York City, I\lew York (Received May 5, 1954)
The author’s previous theory of elasticity and consolidation for isotropic materials m. Appl. Phys. 12, 15.5-164 (1941)] is extended to the general case of anisotropy. The method of derivation is also different and more direct. The particular cases of transverse isotropy and complete isotropy are discussed. 1. INTRODUCTION
EQUATIONS FOR ANISOTROPIC CASE THE
T
sample of bulk volume vb. It is understood that the term “porosity” refers as is customary to the effective porosity, namely, that encompassing only the intercommunicating void spaces as opposed to those pores which are sealed off. In the following, the word “pore” will refer to the effective pores while the sealed pores will be considered as part of the solid. It will be noted that a property of the porosity f is that it represents also a ratio of areas
12, 578-581 (1942).
J. Appl. Phys. 13,35-40
The remaining components uZZ,uZy, etc., of the tensor are the forces applied to that portion of the cube faces occupied by the solid. We shall now call our attention to this system of fluid and solid as a general elastic system with conservation properties. The solid skeleton is considered’to have compressibility and shearing rigidity, and the*.fluid may be compressible. The deformation of a unit cube is assumed to be completely reversible. By deformation is meant here that determined by both strain tensors in the solid and the fluid which will now be defined. The average displacement components of the solid is designated by uZ, uy, u,, and that of the fluid by U,, U,, U,.
u=-fp.
(2.4)
Let us consider an elastic skeleton with a statistical distribution of interconnected pores. This porosity is usually denoted by f= v,/ T/b, where V, * Consultant. 1 M. A.~Biot, 2 M. A. Biot, 8 M. A. Biot (1941). 4 M. A. Biot (2.1)
HE theory of consolidation deals with the settlement under loading of a porous deformable solid containing a viscous fluid. In a previous publication’ a consolidation theory was developed for isotropic materials. The purpose of the present paper is to extend the theory to the most general case of anisotropy. The method by which the theory is derived is also more general and direct. The same physical assumption is introduced, that the skeleton is purely elastic and contains a compressible viscous fluid. The theory may therefore also be considered as a generalization of the theory of elasticity to porous materials. It is applicable to the prediction of the time history of stress and strain in a porous solid in which fluid seepage occurs. The general equations derived in Sec. 2 are applied to the case of transverse isotropy in Sec. 3. This is a case of particular interest in the application of the theory to soils and natural rock formations, since transverse isotropic is the type of symmetry usually acquired by rock under the influence of gravity. For an isotropic material the equations reduce to a simple form given in Sec. 4. They are shown to coincide with the equations derived in reference 1. Application of the theory to specific cases was made previously,2-4 and it was shown that the operational calculus offers a very powerful tool for the solution of consolidation problems in which a load is applied to the material at a given instant and the time history of the settlement is to be calculated. These methods are directly applicable to the more general nonisotropic case. More general solutions of the equations have been developed and will be presented in a forthcoming publication. 2. GENERAL
Theory of Elasticity and Consolidation for a Porous Anisotropic Solid
M. A. BIOT
Reprinted from JOURNAL OF APPLIEDPHYSICS,Vol. 26, No. 2, pp. 182-185, February, 19.55
is the volume of the pores contained in a
J. Appl. Phys. 12, 155-164 (1941). J. Appl. Phys. 12,4X-430 (1941). and F. M. Clingan, J. Appl. Phys. and F. M. Clingan,
f =WSa,
(2.2)
i.e., the fraction S, occupied by the pores in any crosssectional area Sb of the bulk material. It must beassumed, of course, that the pores are randomly distributed in location but not necessarily in direction. That this relation holds may be ascertained by integrating S$Sb over a length of unity in a direction normal to the cross section Sb, The value of this integral then represents the fraction f of the volume occupied by the pores. It is seen that the ratio SP/Sb is also independent of the direction of the cross section. The stress tensor in the porous material is UZ,-tU C,Y ҦI2 auv+u i UZZ UZY UZS UYZ 7 o,,+u I (2.3)