水動力學基本微分方程課件

水動力學基本微分方程水動力學基本微分方程2-1

地下水動力學基本微分方程Date2Date2一、含水層的彈性理論1.含水層的彈性釋水（以承壓水為例）

從承壓含水層中抽出的水，由兩部分組成含水層所貯存水的彈性釋放側向補給（來自遠方）現(xiàn)取一處于平衡狀態(tài)的承壓含水層柱體。設含水層上覆巖層對含水層中的固體顆粒和地下水產(chǎn)生的應力為σ，骨架上的反壓強為，水的頂托力為P。當水處于平衡狀態(tài)時，當從含水層中抽水、水頭下降△h時，圖Date3一、含水層的彈性理論1.含水層的彈性釋水（以承壓水為例）圖D（1）H下降△h，水體積膨脹，從而釋放出一定體積的水；（2）σ保持不變，骨架所受壓力增加，因為固體顆粒接近于剛性體不可壓縮，所以壓力增加引起含水層壓縮，使含水層空隙中的部分地下水被擠出。這兩點就是彈性釋水的機理。

Date4（1）H下降△h，水體積膨脹，從而釋放出一定體積的水；Dat上述分析表明：H降低，承壓含水層釋放部分地下水；H增大，承壓含水層貯存部分地下水，這部分水量稱為彈性貯存量。

彈性貯水量的大小與含水層的巖性和結構有關，為了表征含水層彈性釋水（儲水）的能力，下面將給出彈性貯水率和貯水系數(shù)的概念。Date5上述分析表明：H降低，承壓含水層釋放部分地下水；H2.含水層的貯水率和貯水系數(shù)１.貯水率（Specificstorativity）用表示

當含水層中的水頭降低（或升高）一個單位時，單位時間內(nèi)在單位體積含水層中，由于水的彈性膨脹（或壓縮）及含水層的彈性壓縮（或膨脹）釋放（或貯存）的水量，稱為貯水率,也稱單位貯存量，量綱為［L-1］。Date62.含水層的貯水率和貯水系數(shù)１.貯水率（Specifics

２.貯水系數(shù)（storagecoefficient）

m為含水層厚度，用于二維流計算。無量綱，大部分含水層介于10-5～10-3之間物理意義：在單位面積、厚度為m的含水層柱體中，當水頭降低（或升高）一個單位時，單位時間內(nèi)從含水層中釋放（或貯存）的水量。Date7２.貯水系數(shù)（storagecoefficient）

3.給水度

對潛水含水層而言，當水頭下降時，引起兩部分排水：①含水層下部飽水部分的彈性釋水，其釋水能力用表示；②上部潛水面下降部分引起的重力疏干排水，這部分給水能力用給水度

（Specificyield）表示；給水度的物理意義：當含水層中水頭下降一個單位時，在單位體積含水層中，由重力疏干所排出的水量。Date83.給水度Date84.貯水率與給水度的區(qū)別①彈性釋水由減壓引起，為壓力變化所給出的水量，為重力疏干排出的水量；②貯水率與整個含水層厚度上的巖性、液體性質(zhì)有關，給水度僅與水位波動帶的巖性、液體性質(zhì)有關；③彈性釋放瞬時完成；重力疏干具明顯的滯后效應；④數(shù)量級：約10-5～10-3；約0.1～0.3；∵＞＞∴在潛水含水層中，通常只考慮重力疏干，忽略彈性釋水。Date94.貯水率與給水度的區(qū)別Date9二、承壓含水層中滲流基本微分方程（一）方程的建立設：①地下水、含水層均為彈性體；②彈性釋水瞬時完成；③為常數(shù)（constant）；在滲流場中取一無限小的平行六面體，作為均衡單元，如圖示，六面體邊長分別為dx、dy、dz，下面分析dt時段內(nèi)，均衡單元中的質(zhì)量守恒問題。

均衡區(qū)均衡期Date10二、承壓含水層中滲流基本微分方程均衡區(qū)均衡期Date10

依據(jù)質(zhì)量守恒和能量守恒定律，建立承壓含水層中滲流基本微分方程。

Ｑx為單位時間內(nèi)通過abcd斷面流入的水量。在dt內(nèi)，沿x方向通過abcd斷面流入均衡單元的水量,通過a＇b＇c＇d＇斷面從均衡單元流出的水量為Date11依據(jù)質(zhì)量守恒和能量守恒定律，建立承壓含水層中滲流基dt時段內(nèi)分別沿oy、oz進入單元的水量為而流出的水量為

在dt內(nèi)，均衡單元貯存量的變化量為:

據(jù)水均衡原理得：

Date12dt時段內(nèi)分別沿oy、在dt內(nèi)，均衡單元貯存量的變化化簡為：根據(jù)達西定律：Date13化簡為：根據(jù)達西定律：Date13

上式就是非均質(zhì)各向異性承壓含水層中地下水三維非穩(wěn)定運動的基本微分方程。對各向異性介質(zhì)，取坐標軸方向與主滲透方向一致。（二）方程的化簡和討論

１.對于均質(zhì)各向同性含水層，Ｋ為常數(shù)，這時簡化為：Date14上式就是非均質(zhì)各向異性承壓含水層中地下水２.對于二維的情況，常用和T表示(各項均乘以m）

３.當含水層有垂直水量交換時，其量常用Ｗ表示，稱為源（匯）項，含水層的源（匯）項可是t和位置的函數(shù)W=W(x,y,z,t)。

當從含水層中抽水或從垂直方向有水流出含水層時，Ｗ為負，稱為匯；當給含水層中注水或從垂直方向有水流入含水層時，Ｗ為正，稱為源；Date15２.對于二維的情況，常用和T表示(各項均乘以m）３

二維：W為單位時間從單位體積含水層中流入或流出的水量

三維：W為單位時間從垂直方向單位面積含水層中流入或流出的水量(補給強度或蒸發(fā)強度)

Date16二維：W為單位時間從單位體積含水層中流入或流出的水量三維Date17Date17

可得到三維、二維相應的穩(wěn)定狀態(tài)下的滲流基本微分方程。對于均質(zhì)各向同性的三維流來說：

5.若化為柱坐標（三維各向同性介質(zhì)）Date18可得到三維、二維相應的穩(wěn)定狀態(tài)下的滲流基本三、越流含水層中滲流基本微分方程（一）什么是越流？

1.半承壓含水層：一個主含水層的上層和（或）下層為弱透水層，主含水層通過弱透水層與相鄰含水層發(fā)生水力聯(lián)系，但它本身具有承壓性，主含水層稱為半承壓含水層。2.越流：當半承壓含水層與相鄰含水層間存在水頭差時，地下水便通過弱透水層從高水頭含水層向低水頭含水層產(chǎn)生垂向流動，這種現(xiàn)象稱為越流。簡言之，相鄰含水層在水頭差的作用下，通過弱透水層與主含水層發(fā)生水力聯(lián)系的現(xiàn)象稱為越流。

越流的方向：由兩相鄰含水層的水位決定。Date19三、越流含水層中滲流基本微分方程（一）什么是越流？Date1Date20Date20（二）越流含水層中滲流基本微分方程

1.假定a.忽略弱透水層的彈性釋水；b.水流在弱透水層中是垂向運動，而在主含水層中折射為水平運動；

2.方程的建立

在主含水層中取一微分柱體（其長寬分別為dx、dy，高為含水層厚度m）作為均衡單元。下面分析在dt時段內(nèi)，微分柱體的水均衡問題。Date21（二）越流含水層中滲流基本微分方程Date21P(x，y）Date22P(x，y）Date22沿x方向流入單元體的水量：

流出：

沿y方向流入單元體的水量：

流出：

沿z方向流入單元體的水量：

流出：

Date23沿x方向流入單元體的水量：沿y方向流入單元體的水量：沿z方向流入量－流出量＝：單元體內(nèi)貯存量的變化為：根據(jù)水均衡原理得：Date24流入量－流出量＝：單元體內(nèi)貯存量的變化為：Date24代入上式得：這就是不考慮弱透水層彈性釋水條件下，非均質(zhì)各向異性越流含水層系統(tǒng)地下水非穩(wěn)定運動的基本微分方程。Date25代入上式得：這就是不考慮弱透水層彈性釋水條件下，非均質(zhì)（三）越流因素和越流系數(shù)１.越流因素（Leakagefactor）

上式中，若介質(zhì)為均質(zhì)各向同性介質(zhì)，T＝Constant式中：Ｔ為主含水層的導水系數(shù)Ｋi，mi分別為弱透水層的滲透系數(shù)和厚度Date26（三）越流因素和越流系數(shù)式中：Ｔ為主含水層的導水系數(shù)DatＢ1、Ｂ2分別稱為上、下兩弱透水層的越流因素

越流因素Ｂ是反映弱透水層隔水性能的參數(shù)。Ｂ越大，越流量Ｑ越小，對于隔水層Ｂ無窮大，越流量為零。

Date27Ｂ1、Ｂ2分別稱為上、下兩弱透水層的越流因素越流因素Ｂ２.越流系數(shù)（Coefficientofleakage）

因為物理意義：當主含水層和相鄰含水層間的水頭差等于一個長度單位時，通過單位面積含水層上的越流量。越流系數(shù)反映越流量的大小，越大，相同水頭下的越流量也越大。

Date28２.越流系數(shù)（Coefficientofleakage）四、潛水含水層中滲流基本微分方程（一）Dupuit假定

潛水面通常不是水平面，潛水含水層中存在著流速的垂直分量，潛水面本身又是滲流區(qū)邊界，隨時間而變化。為了建立潛水含水層中滲流基本微分方程，引出了Dupuit假定：假設潛水面比較平緩，潛水面上任意一點P有：Date29四、潛水含水層中滲流基本微分方程（一）Dupuit假定假設潛

相當于忽略了滲透速度的垂直分量，代替，在鉛垂面上各點的水頭都是相等的；或者說，水頭不隨深度而變化，同一鉛直面上各點的水力坡度和滲透速度都相等，滲透速度可表示為：

在此條件下，通過寬度為B的鉛直平面沿x方向的流量Ｑx為：

h-含水層厚度，當隔水底板水平時，h=HDate30相當于忽略了滲透速度的垂直分量（二）潛水含水層中滲流基本微分方程1.方程的建立

在Dupuit假定下，考慮一維問題，取平行于xoz平面的單位寬度進行研究。

首先取一微分柱體（其長度為?x，寬為１，高為整個含水層厚度）作為均衡單元，下面分析在dt時段內(nèi)，微分柱體的水均衡問題。Date31（二）潛水含水層中滲流基本微分方程首先取一微分柱從上游斷面流入：從下游斷面流出：

在?t時間內(nèi)，垂直方向的補給量為：

由于潛水面的上升而引起的均衡區(qū)內(nèi)的水的增量為：

Date32從上游斷面流入：在?t時間內(nèi)，垂直方向的由于潛根據(jù)連續(xù)性原理，上面兩個增量應相等，即代入上式得到

上式為非均質(zhì)各向異性潛水含水層中一維地下水非穩(wěn)定運動的基本微分方程，也稱為Ｊ.Boussinesq方程。Date33根據(jù)連續(xù)性原理，上面兩個增量應相等，即代入上式得到上對于均質(zhì)含水層，Ｋ＝Constant

式中：：水位下降時稱為給水度，水位上升時稱為飽和差；W：降雨入滲強度（+）或蒸發(fā)強度（—）；h：含水層厚度；H：含水層的水位（平均值）；K：含水層滲透系數(shù)；注意：

a.Ｈ為整個含水層厚度上的平均值；b.Ｈ、h均為未知，所以該方程為二階非線性偏微分方程；c.該方程不適于水力梯度較大地段；不能計算任一點的Ｈ。Date34對于均質(zhì)含水層，Ｋ＝Constant式中：：水位下降時稱

2.方程討論

a.二維情況：b.均質(zhì)含水層，Ｋ＝Constantc.當隔水層底板水平時，H＝h；Date352.方程討論a.