地下水動(dòng)力學(xué)地下水流基本微分方程及定解條件(1)

1、第二章第二章 地下水流基本微分方程及定解條件地下水流基本微分方程及定解條件基本理論：連續(xù)性假設(shè)基本理論：連續(xù)性假設(shè)+ +達(dá)西定律達(dá)西定律+ +水均衡原理水均衡原理對(duì)各種水流問題建立基本微分方程及數(shù)學(xué)模型對(duì)各種水流問題建立基本微分方程及數(shù)學(xué)模型： 按空間維數(shù)：一維、二維（平面二維、剖面二維）、三維按空間維數(shù)：一維、二維（平面二維、剖面二維）、三維 按含水層類型：承壓水流、潛水流、多層（越流聯(lián)系）等按含水層類型：承壓水流、潛水流、多層（越流聯(lián)系）等求解數(shù)學(xué)模型（利用解析法），得到一些典型解析解：求解數(shù)學(xué)模型（利用解析法），得到一些典型解析解： 裘布依穩(wěn)定井流模型裘布依穩(wěn)定井流模型 無越流承壓含水層

2、中的完整井流（泰斯模型）無越流承壓含水層中的完整井流（泰斯模型） 無越流潛水含水層中的完整井流（博爾頓模型無越流潛水含水層中的完整井流（博爾頓模型-考慮滯后給水、考慮滯后給水、 紐曼模型紐曼模型-考慮流速垂直分量和彈性儲(chǔ)量）考慮流速垂直分量和彈性儲(chǔ)量） 越流系統(tǒng)中的承壓完整井流模型越流系統(tǒng)中的承壓完整井流模型應(yīng)用：應(yīng)用： 預(yù)測(cè)抽水條件下的水頭變化；預(yù)測(cè)抽水條件下的水頭變化； 利用抽水試驗(yàn)資料求含水層參數(shù)。利用抽水試驗(yàn)資料求含水層參數(shù)。第二章第二章 地下水流基本微分方程及定解條件地下水流基本微分方程及定解條件教學(xué)目標(biāo)：教學(xué)目標(biāo)：準(zhǔn)確理解滲流連續(xù)性概念準(zhǔn)確理解滲流連續(xù)性概念掌握達(dá)西定律和質(zhì)量守恒原

3、理的應(yīng)用掌握達(dá)西定律和質(zhì)量守恒原理的應(yīng)用掌握建立地下水基本微分方程的思想方法掌握建立地下水基本微分方程的思想方法幾種典型的地下水流方程的推導(dǎo)幾種典型的地下水流方程的推導(dǎo) 潛水剖面二維流、平面二維流潛水剖面二維流、平面二維流 承壓水二維流承壓水二維流 三維流三維流邊界條件概化，初始條件確定方法與原則邊界條件概化，初始條件確定方法與原則能夠用數(shù)學(xué)模型描述實(shí)際問題能夠用數(shù)學(xué)模型描述實(shí)際問題第二章第二章 地下水流基本微分方程及定解條地下水流基本微分方程及定解條件件主要內(nèi)容：主要內(nèi)容：建立連續(xù)性方程建立連續(xù)性方程分析含水層與巖石、流體壓縮性關(guān)系分析含水層與巖石、流體壓縮性關(guān)系建立不同含水層地下水流微分方

4、程建立不同含水層地下水流微分方程討論邊界條件及初始條件討論邊界條件及初始條件 用數(shù)學(xué)模型描述實(shí)際問題用數(shù)學(xué)模型描述實(shí)際問題 2.1 2.1 水和多孔介質(zhì)的壓縮性水和多孔介質(zhì)的壓縮性水的壓縮方程（地下水的狀態(tài)方程水的壓縮方程（地下水的狀態(tài)方程)0)(00000ln)( 1VVeVVppdVVdpppVVppdpdVVVdVdp11假定水近似地符合彈性變形，依虎克定律，有 p 為水壓； V 為水的體積；為水的體積彈性壓縮(或膨脹)系數(shù)E為體積彈性模量 。V隨p增大而減小，即dV/dp0 積分1E壓縮系數(shù)：壓縮系數(shù)：?jiǎn)挝粔毫ψ兓瘯r(shí)引起的液體單位體積的單位壓力變化時(shí)引起的液體單位體積的變化量，單位為平

5、方米每牛。其倒數(shù)為體積模量，變化量，單位為平方米每牛。其倒數(shù)為體積模量，單位為帕斯卡。水體壓縮系數(shù)與壓力和溫度有關(guān)。單位為帕斯卡。水體壓縮系數(shù)與壓力和溫度有關(guān)。水的壓縮方程水的壓縮方程0( )22323()23000(0)(0)(0)( )( )1!2!1.2!3!1()()().2!3!nnxppffff xf xxxxnxxexepppppp )()(1 )(1)(100000000)(0ppVVppVVppVVppepp0)(0VVepp按Taylor級(jí)數(shù)展開由于很小，且p變化不大，故)()()(00000000ppVVppVVVppVVV水的壓縮方程水的壓縮方程pVV0pVVVdVdp

6、11 (1 5)dddpdp 由于VV0變化不大，故由于ddmmdVdVmV)1()(水的壓縮方程水的壓縮方程多孔介質(zhì)的壓縮方程多孔介質(zhì)的壓縮方程vbvsbbbdVdVVVVVdVdpdpd1 bbVdVd1假定多孔介質(zhì)近似地符合彈性變形，依虎克定律，有 為巖土的體積彈性壓縮系數(shù)。如果上部荷載不變，則由于骨架部分體積不變VvnVb；Vs(1-n)Vb111svsvbbsvdVdVdVdVnnV dV dVdV d VvnVb；Vs(1-n)Vb式中 多孔介質(zhì)固體顆粒壓縮系數(shù)，表示多孔介質(zhì)中固體顆粒本身的壓縮性的指標(biāo)，sp； 多孔介質(zhì)中孔隙壓縮系數(shù)（Compressibility of the

7、pores of a porous medium），表示多孔介質(zhì)中孔隙的壓縮性的指標(biāo)。n多孔介質(zhì)的孔隙度。 1-8因 ，故 。ddVVddVVvbsb11ddVVnddVVnvvss1psnn)1 (ddVVsss1ddVVvvp1(1)spnpn1-9說明本節(jié)假設(shè)：假定多孔介質(zhì)變形符合彈性規(guī)律，對(duì)研究含水層釋水時(shí)可用；但對(duì)研究地面沉降問題時(shí)，應(yīng)用有所差異。水的壓縮方程水的壓縮方程pVVVdVdp1dpd多孔介質(zhì)的壓縮方程多孔介質(zhì)的壓縮方程)1 (epedHdpbbVdVd1psnn)1 (pn測(cè)壓水頭圖2-2-1 飽和含水介質(zhì)中受力情況2.2 2.2 水和多孔介質(zhì)的壓縮性水和多孔介質(zhì)的壓縮性

8、spphp地下水彈性儲(chǔ)存概念地下水彈性儲(chǔ)存概念取一典型處于平衡狀態(tài)的飽和地層柱體來研究，這里只考慮垂直一維壓密，忽略側(cè)面上粒間力(包括內(nèi)聚力和摩擦力)的作用。含水層上覆巖土體、地表建筑物和大氣壓力等荷載形成的總壓應(yīng)力總壓應(yīng)力 由粒間應(yīng)力的垂向分量垂向分量 s s和孔隙水應(yīng)力孔隙水應(yīng)力p p兩者來平衡. ps)1 (為單位水平面積中顆粒間接觸面積的水平投影.由于 0， 入滲vW0， 蒸發(fā)W定義為潛水面處單位水平面積、單位時(shí)間的入滲量定義為潛水面處單位水平面積、單位時(shí)間的入滲量布西涅斯克方程推導(dǎo)),(txxh),(txxq),(txq),(txhtthth圖2-4-2 潛水剖面二維流均衡要素圖 x

9、研究剖面二維流（x-z)均衡單元體：長(zhǎng)度為x，寬度為1個(gè)單位的含水層柱體均衡 V=流入流出均衡時(shí)段 ttqtqtQV1txWtqtqVtxxtx,x方向流入流出差tqtqtxxtx,z方向流入 上邊界潛水面邊界txWtxWtWV1下邊界隔水邊界 V=0dtxttxxhhV1,x,z方向流入流出差單元體內(nèi)水體積的增量表現(xiàn)為水位在該時(shí)段內(nèi)的上升或下降dtxttxtxxtxxhhVtxWtqtqV1,),(txxh),(txxq),(txq),(txhtthth圖2-4-2 潛水剖面二維流均衡要素圖 xxhKhqthWxqtxxd0, 0水均衡方程:dtxttxtxxtxdtxttxtxxtxthh

10、WxqqxhhtxWqq,thWxhKhxd布西涅斯克方程推導(dǎo)thWxhKhxd注意:1. d為重力給水度，是在重力作用下單位水平面積的潛水含水層柱體在其潛水面下降或上升單位水頭時(shí)釋放或儲(chǔ)存的水量。2. 上述推導(dǎo)運(yùn)用了裘布依假定，忽略垂向分流速，因此可近似用x代替s.3. 布西涅斯克方程將刻畫潛水面邊界問題簡(jiǎn)單化為用W直接表示。布西涅斯克微分方程線性化布西涅斯克方程線性化布西涅斯克方程tWhxKxhhtWhxKxhhhthWxhKxhhhdmmddmm即端的再以平均值代替方程左，有并令方程兩端均乘方法代替并視為常量用平均值中的將方程左端第一項(xiàng)括號(hào)方法,21,: )2(,: ) 1 (2tWyK

11、yhxKxhthWyhKyhxhKxhtHWyhKhyxhKhxdmmdmmd:,性化后的方程分別為相應(yīng)第一種和第二種線所得布西涅斯克方程為推廣到三維流問題線性化布西涅斯克方程線性化布西涅斯克方程22222222222222,(/)(/),11mdmdmdKhxytTTKhaTLTLTaxytKhrrrtrrrt對(duì)于均質(zhì)含水層，且無垂向補(bǔ)排時(shí)，可得引入和 分別為潛水含水層的導(dǎo)水系數(shù)和水力傳播系數(shù)對(duì)于軸對(duì)稱潛水流問題 采用極坐標(biāo)系和線性化布西涅斯克方程線性化布西涅斯克方程2.5 2.5 地下水流動(dòng)定解條件及數(shù)學(xué)模型地下水流動(dòng)定解條件及數(shù)學(xué)模型地下水流動(dòng)控制微分方程n潛水二維不穩(wěn)定流動(dòng)控制方程n承

12、壓水二維不穩(wěn)定流動(dòng)控制方程定解條件n邊界條件n初始條件tHyHZHKyxHZHKxdyyxx)()(tHyHMKyxHMKxeyyxx)()(2.5 2.5 地下水流動(dòng)定解條件及數(shù)學(xué)模型地下水流動(dòng)定解條件及數(shù)學(xué)模型1 數(shù)學(xué)模型的有關(guān)概念數(shù)學(xué)模型的有關(guān)概念 同一形式的偏微分方程代表了整個(gè)一大類的地下水流的運(yùn)動(dòng)規(guī)律，而對(duì)于不同邊界性質(zhì)、不同邊界形狀的含水層，水頭的分布是不同的。而且對(duì)于偏微分方程而言，方程本身并不包含反映特定滲流區(qū)條件的全部信息，方程可能存在無數(shù)個(gè)解，如需要從大量的可能解中求得與特定區(qū)域條件相對(duì)應(yīng)的唯一特解，就必須提供反映特定區(qū)域特征的信息。這些信息包括： （1）微分方程中的有關(guān)參

13、數(shù)， 當(dāng)這些參數(shù)確定后，微分方程才能被確定下來。WT,（2）滲流區(qū)范圍和形狀，當(dāng)微分方程所對(duì)應(yīng)的區(qū)域被確定之后才能對(duì)方程求解。（3）邊界條件邊界條件(boundary conditions)：表示滲流區(qū)邊界所處的條件，用以表示水頭H（或滲流量q）在滲流區(qū)邊界上所應(yīng)滿足的條件，也就是滲流區(qū)內(nèi)水流與其周圍環(huán)境相互制約的關(guān)系。（4）初始條件初始條件(initial conditions)：表示滲流區(qū)的初始狀態(tài)，某一選定的初始時(shí)刻(t=0)滲流區(qū)內(nèi)水頭H的分布情況。將邊界條件邊界條件和初始條件初始條件并稱為定解條件定解條件, 指水頭、流量等滲流運(yùn)動(dòng)要素在流場(chǎng)邊界上的已知變化規(guī)律，這種變化規(guī)律是由流場(chǎng)外

14、部條件引起的，但它不斷地影響流場(chǎng)內(nèi)部的滲流過程并在整個(gè)期間一直起作用。微分方程和定解條件一起構(gòu)成滲流場(chǎng)的數(shù)學(xué)模型。數(shù)學(xué)模型。 數(shù)學(xué)模型數(shù)學(xué)模型：描述某一研究區(qū)地下水流運(yùn)動(dòng)的數(shù)學(xué)方程與其定解條件共同構(gòu)成的表示某一實(shí)際問題的數(shù)學(xué)結(jié)構(gòu)。亦即從物理模型出發(fā)，用簡(jiǎn)潔的數(shù)學(xué)語言，即一組數(shù)學(xué)關(guān)系式來刻畫它的數(shù)量關(guān)系和空間形式，從而反映所研究地質(zhì)體的地質(zhì)、水文地質(zhì)條件和地下水運(yùn)動(dòng)的基本特征，達(dá)到復(fù)制或再現(xiàn)一個(gè)實(shí)際水流系統(tǒng)基本狀態(tài)的目的的一種數(shù)學(xué)結(jié)構(gòu)。其中微分方程表示地下水的流動(dòng)規(guī)律，定解條件表明研究對(duì)象所處的特定環(huán)境條件，即所研究的地下水流的真實(shí)狀態(tài)。定解問題是給定了方程（或方程組）和相應(yīng)定解條件的數(shù)學(xué)物理問

15、題。建立模型是指建立數(shù)學(xué)模型的過程。二、定解條件二、定解條件(1) 第一類邊界條件第一類邊界條件(Dirichlet條件)：如果在某一部分邊界(設(shè)為S1或1)上，各點(diǎn)在每一時(shí)刻的水頭都是已知的，則這部分邊界就稱為第一類邊界或給定水頭的邊界，表示為:或 給定水頭邊界不一定就是定水頭邊界。可以作為第一類邊界條件來處理的情況： 河流或湖泊切割含水層，兩者有直接水力聯(lián)系時(shí)，這部分邊界就可以作為第一類邊界處理。此時(shí)，水頭是一個(gè)由河湖水位的統(tǒng)計(jì)資料得到的關(guān)于t的函數(shù)。但要注意，某些河、湖底部及兩側(cè)沉積有一些粉砂、亞粘土和粘土，使地下水和地表水的直接水力聯(lián)系受阻，就不能作為第一類邊界條件來處理。11( ,

16、, , )|( , , , ),( , , )sH x y z tx y z tx y zs11( , , , )|( , , ),( , )H x y z tx y tx y 區(qū)域內(nèi)部的抽水井、注水井或疏干巷道也可以作為給定水頭的內(nèi)邊界來處理。此時(shí)，水頭通常是按某種要求事先給定，例如給定抽水井的允許降深等。上面介紹的都只是給定水頭的邊界。注意，給定水頭邊界不一定是定水頭邊界。 排泄地下水的溢出帶、沖溝或排水渠的邊界也可近似看作給定水頭邊界。(2)第二類邊界條件(Neumam條件)：當(dāng)知道某一部分邊界(設(shè)為S2或2)單位面積(二維空間為單位寬度)上流入(流出時(shí)用負(fù)值)的流量q時(shí)，稱為第二類邊界

17、或給定流量的邊界。相應(yīng)的邊界條件表示為:式中，n為邊界S2或2的外法線方向。q1和q2則為已知函數(shù)，分別表示S2上單位面積和2上單位寬度的側(cè)向補(bǔ)給量。常見的這類邊界條件： 隔水邊界（流線、分水嶺）：0nH2221),(),(),(),(yxtyxqnHTStzyxtzyxqnHK或 抽水井或注水井： 補(bǔ)給或排泄地下水的河渠邊界上，如已知補(bǔ)給量。（3）第三類邊界條件：某邊界上H和 的線性組合是已知的，即有： 又稱混合邊界條件，為已知函數(shù)。邊界為弱透水層（滲透系數(shù)為K1，厚度或?qū)挾葹閙1），wwrQnHT2（1-108）nH（1-109）nHHnH11Km),()(11tzyxqHHmKnHKn在

18、s3上,在 3上,浸潤(rùn)曲線的邊界條件：當(dāng)浸潤(rùn)曲線下降時(shí)，從浸潤(rùn)曲線邊界流入滲流區(qū)的單位面積流量q為：式中，為給水度，為浸潤(rùn)曲線外法線與鉛垂線間的夾角。（1-110）（1-111）qnHKc2（1-112）（1-113）03HHnHKnS03HHMnHTnqnHKc2cos*tHq 一、定解條件潛水面邊界 tHWzHKxHKztzyxHpdzzxx2)(),(0tHWzHKyHKxHKdzzyyxx22)()(三維條件下 3) 初始條件初始條件:某一選定的初始時(shí)刻(t=0)滲流區(qū)內(nèi)水頭H的分布情況?；?其中，H0為D上的已知函數(shù)。DzyxzyxHtzyxHt),(),(| ),(00（1-114

19、）DzyxyxHtyxHt),(),(| ),(00（1-115） (1)線性、非線性模型模型由線性方程所組成，稱為線性模型，如均質(zhì)各向同性承壓二維流方程。模型由非線性方程所組成，稱為非線性模型，如潛水模型方程。 (2)靜態(tài)、動(dòng)態(tài)模型 根據(jù)模型中未知變量與時(shí)間的關(guān)系進(jìn)行劃分，若未知變量與時(shí)間無關(guān)，如穩(wěn)定流模型，稱為靜態(tài)模型，反之，則為動(dòng)態(tài)模型。 (3)集中、分布參數(shù)模型模型中不含有空間坐標(biāo)變量的模型，稱為集中參數(shù)模型，如抽水井流量與降深之間的經(jīng)驗(yàn)公式。模型中含有空間坐標(biāo)變量的模型，稱為分布參數(shù)模型。 (4)確定性與隨機(jī)性模型確定性模型確定性模型：數(shù)學(xué)模型中各變量之間有嚴(yán)格的數(shù)學(xué)關(guān)系的模型。 隨

20、機(jī)性模型隨機(jī)性模型：數(shù)學(xué)關(guān)系式中含有一個(gè)或多個(gè)隨機(jī)變量的模型。三、三、 滲流數(shù)學(xué)模型的分類滲流數(shù)學(xué)模型的分類 用確定性模型來描述實(shí)際地下水流時(shí)，如前述，必須具備下列條件： 有一個(gè)(或一組)能描述這類地下水運(yùn)動(dòng)規(guī)律的偏微分方程；同時(shí)，確定了相應(yīng)滲流區(qū)的范圍、形狀和方程中出現(xiàn)的各種參數(shù)值。 給出相應(yīng)的定解條件。對(duì)所建立的模型進(jìn)行檢驗(yàn)，即把模型預(yù)測(cè)的結(jié)果與通過抽水試驗(yàn)或其它試驗(yàn)對(duì)含水層施加某種影響后所得到的實(shí)際觀測(cè)結(jié)果或一個(gè)地區(qū)地下水動(dòng)態(tài)長(zhǎng)期觀測(cè)資料進(jìn)行比較，看兩者是否一致。若不一致，就要對(duì)模型進(jìn)行校正，即修正條件(1)和(2)直至滿意擬合為止。這一步驟稱為識(shí)別模型或校正模型。 經(jīng)過校正后的模型，能

21、代表所研究的地質(zhì)體，或者說是實(shí)際水流系統(tǒng)的復(fù)制品了，因而可以根據(jù)需要，用這個(gè)模型進(jìn)行計(jì)算或預(yù)測(cè)，例如預(yù)測(cè)礦床疏干時(shí)的涌水量及地下水污染情況預(yù)測(cè)等。 解(即滿足條件和的解)是存在的(存在性)； 解是唯一的(唯一性)；要求所提問題的解存在和唯一是不言而喻的。 解對(duì)原始數(shù)據(jù)是連續(xù)依賴的(穩(wěn)定性)。即穩(wěn)定性的要求，意味著當(dāng)參數(shù)或定解條件發(fā)生微小變化時(shí)，所引起的解的變化也是很微小的。只有有了這條保證，當(dāng)參數(shù)和定解條件的數(shù)據(jù)有某些誤差時(shí)，所求得的解才能仍然接近于真解；否則，解是不可信的，并應(yīng)該認(rèn)為此時(shí)的數(shù)學(xué)模型是有毛病的。在實(shí)際工作中，原始數(shù)據(jù)有某種誤差，在所難免，所以這個(gè)條件很重要。 適定問題（Well

22、 posed problem ）是指數(shù)學(xué)模型滿足（1）解是存在的（存在性），（2）解是唯一的（唯一性），（3）解對(duì)原始數(shù)據(jù)是連續(xù)依賴的（穩(wěn)定性）這三個(gè)條件的問題。只要有一條不滿足就是不適定問題。正問題是根據(jù)數(shù)學(xué)模型、給定的含水層水文地質(zhì)參數(shù)和定解條件求解水頭的問題，又稱水頭預(yù)報(bào)問題。逆問題（inverse problem）是根據(jù)數(shù)學(xué)模型、動(dòng)態(tài)觀測(cè)資料或抽水試驗(yàn)資料反過來確定含水層水文地質(zhì)參數(shù)的問題。4. 建立數(shù)學(xué)模型的基本要點(diǎn)建立數(shù)學(xué)模型的基本要點(diǎn)(1)確定研究區(qū)的范圍及滲流區(qū)的邊界；(2)確定滲流區(qū)的水力特征（包括埋藏條件、滲流狀態(tài)、介質(zhì)特征）；(3)確定滲流區(qū)的邊界條件；(4)確定滲流區(qū)的源匯項(xiàng)；(5)選擇微分方程；(6)確定滲流區(qū)的初始條件。EDBAABqC線河流#1流#2#3均質(zhì)、各向同性的潛水含水層，地下水流為平面非穩(wěn)定流，且與河水有直接的水力聯(lián)系。