第七章 地下水運(yùn)動(dòng)中的若干專門(mén)問(wèn)題.doc

第七章 地下水運(yùn)動(dòng)中的若干專門(mén)問(wèn)題1 非飽和帶的地下水運(yùn)動(dòng)一、關(guān)于非飽和帶水分的基本知識(shí)1. 含水率，飽和度和田間持水量包氣帶中的空隙，一部分被水充填，另一部分被空氣充填。含水率()：表示單位積中水所占的體積，式中：（Vw）0典型單元體中水的體積；V0典型單元體的體積飽和度：巖石的空隙空間中被水占據(jù)部分所占的比例。式中：（V0）0典型單元體中的空隙體積含水率與飽和度的關(guān)系：=nSw式中：n孔隙度。田間持水量：在長(zhǎng)時(shí)間重力排水后仍然保留在土中的水量。2. 毛管壓力毛管壓強(qiáng)：在多孔介質(zhì)的孔隙中，液體和氣體接觸是，二者存在壓力差，這個(gè)壓力差稱毛管壓強(qiáng)。用pc表示pc=pa-pw式中：pa空氣的壓強(qiáng)；pw水的壓強(qiáng)毛管壓強(qiáng)取決于界面的曲率，曲率愈大（液面愈彎曲，毛管壓強(qiáng)愈大。以上毛管壓強(qiáng)是以絕對(duì)壓強(qiáng)為基準(zhǔn)，如果以相對(duì)壓強(qiáng)為基準(zhǔn)，這時(shí)有：pc=pa-pw pa pc=-pw 毛管壓強(qiáng)相對(duì)大氣壓強(qiáng)為負(fù)值。即，非飽和帶孔隙中的水處于小于大氣壓強(qiáng)的情況下。非飽和帶水流中任何點(diǎn)的水頭式中：z位置水頭；hc=pc/r 毛管壓力水頭； H=z-hc3土壤水分特征的曲線水分特征曲線：反映毛管壓力水頭（或毛管壓強(qiáng)）和土壤含水率或飽和度關(guān)系的曲線。如圖：隨著含水率的減少，毛管壓力增加，當(dāng)含水率減小到某一值時(shí)，壓強(qiáng)繼續(xù)增大時(shí)，含水率不在減小。相應(yīng)的飽和度為：影響特征曲線的因素：（1）不同質(zhì)地的土壤，其水分特征曲線不同。一般說(shuō)，土壤的粘粒含量愈高。同一負(fù)壓條件下土壤的含水率愈大，或者同一含水率下其負(fù)壓愈高。這是因?yàn)?，粘粒含量增多。使土壤中?xì)小孔隙發(fā)育的緣故。（2）土壤結(jié)構(gòu)。如圖，為一砂壤土不同干容重的水分特征曲線，在同一負(fù)壓下，土壤愈密實(shí)，（大），相應(yīng)的含水率一般也大。原因，土壤愈密實(shí)，大孔隙數(shù)量減少，中孔隙增多。（3）溫度的影響。溫度升高，水的粘滯性下降，所以表面張力降低，在同樣的負(fù)壓下，含水率要低一些。（4）土壤水分變化過(guò)程的影響。對(duì)于同一土壤，土壤脫濕（由濕變干）過(guò)程測(cè)得的水分特征曲線不同，如圖，在相同的負(fù)壓下，排水（脫濕）時(shí)的含水率要大于吸濕時(shí)的含水率。這種現(xiàn)象稱為滯后現(xiàn)象。（5）容水度：毛管壓力水頭變化一個(gè)單位時(shí)，從單位體積土中釋放出的水體積。數(shù)值上等于，水分特征曲線的斜率的負(fù)倒數(shù)。它是含水率和毛管壓強(qiáng)的函數(shù)，可用或表示。二、非飽和帶水運(yùn)動(dòng)的基本方程非飽和帶地下水的運(yùn)動(dòng)，也可以用達(dá)西定律描述，這時(shí)的滲透系數(shù)是變化的，與含水率和毛細(xì)壓力水頭有關(guān)，是和的函數(shù)，其關(guān)系如圖，隨著含水率的增大，滲透系數(shù)增大，隨毛管壓力水頭的減小，滲透系數(shù)增大。在非飽和帶中，定律的表達(dá)式為：在三個(gè)坐標(biāo)軸的分量為：2 基本微分方程第一章推的滲流連續(xù)性方程，如下：在飽水帶中，全部孔隙被水充滿，等式右端用孔隙度，在非飽和帶中，部分孔隙被水充滿，所以用含水率取代，并兩邊除（近似為常數(shù)），得：將vx、vy、vz代入上式，得：二式為非飽和流的基本方程3 基本方程的幾種形式（1）以含水率為因變量的表達(dá)式將H換成，將H=z-hc代入上（1）式，得：上式進(jìn)一步變換定義K（）/C（）=D（）為擴(kuò)散系數(shù)，得：對(duì)于垂向一維流動(dòng)，去掉前兩項(xiàng)，得：軸向上取正值，軸向下取負(fù)值。（2）以毛管壓力水頭為因變量得表達(dá)式：將H=z-hc代入（2）式，得：代入上式：對(duì)于垂向一維流動(dòng)，去掉前兩項(xiàng)，得：軸向上取正值，軸向下取負(fù)值。2 水動(dòng)力彌散理論用來(lái)模擬地下水中污染物和化學(xué)成分得運(yùn)移過(guò)程，預(yù)測(cè)地下水污染得發(fā)展趨勢(shì)。一 水動(dòng)力彌散現(xiàn)象及其機(jī)理例1. 在一口井中注入一種示蹤劑，示蹤劑在隨地下水向前流動(dòng)得過(guò)程中，向外圍擴(kuò)散，形成一個(gè)以中心點(diǎn)濃度最大，向四周濃度逐漸減小的過(guò)渡帶，并且隨示蹤劑遷移的距離增大，過(guò)渡帶也越來(lái)越寬。如圖（書(shū)中）例2. 在均勻流的砂柱中，用含有示蹤劑濃度為的水去替代，在砂柱另一端測(cè)量示蹤劑濃度，得曲線如圖（書(shū)）。水在流動(dòng)過(guò)程中并非一個(gè)突變界面，而是一個(gè)過(guò)液帶。這種現(xiàn)象稱水動(dòng)力彌散。水動(dòng)力彌散是機(jī)械彌散和分子擴(kuò)散所引起的。1 機(jī)械彌散液體在多孔介質(zhì)中運(yùn)動(dòng)的三種情況：（1）由于液體粘性的作用和結(jié)合水的摩擦阻力，使得靠近孔隙壁的水流速度趨于零?？紫吨行牟课涣魉僮畲蟆？紫洞笮〔灰唬斐刹煌紫吨g沿軸部的最大流速有差異；（1）由于空隙的彎彎曲曲，水流方向也隨之不斷地改變。由于上述三種情況，造成了地下水質(zhì)點(diǎn)運(yùn)動(dòng)速度，在大小和方向上的不均一，造成了示蹤劑有的運(yùn)動(dòng)快，有的運(yùn)動(dòng)慢，從而形成了上述過(guò)渡帶。這種由于速度不均一所造成的這種物質(zhì)運(yùn)移現(xiàn)象稱為機(jī)械彌散。2 分子擴(kuò)散一般溶質(zhì)都有由濃度高向濃度低的地方運(yùn)移的性質(zhì)，以求濃度趨于均一。這種由于液體中所含溶質(zhì)的濃度不均一而引起的物質(zhì)運(yùn)移現(xiàn)象叫分子擴(kuò)散。分子擴(kuò)散服從定律式中：Is單位時(shí)間內(nèi)通過(guò)單位面積的溶質(zhì)的質(zhì)量；dc/ds溶質(zhì)在溶液中的濃度c沿s方向變化的濃度梯度；Dd擴(kuò)散系數(shù)。機(jī)械彌散和分子擴(kuò)散是同時(shí)出現(xiàn)的，當(dāng)流速較大時(shí)，機(jī)械彌散是主要的；當(dāng)流速甚小時(shí)，分子擴(kuò)散的作用就變得明顯。水動(dòng)力彌散還分為沿水流方向和垂直與水流方向的彌散，沿水流方向的彌散稱縱向彌散，垂直水流方向的彌散稱橫向彌散。二 水動(dòng)力彌散系數(shù)分子擴(kuò)散服從定律：式中：D為分子擴(kuò)散系數(shù)；I為由于分子擴(kuò)散在單位時(shí)間內(nèi)通過(guò)單位面積的溶質(zhì)質(zhì)量。機(jī)械彌散也服從定律：式中：D為機(jī)械擴(kuò)散系數(shù)；I為由于機(jī)械擴(kuò)散在單位時(shí)間內(nèi)通過(guò)單位面積的溶質(zhì)質(zhì)量。由于水動(dòng)力彌散是分子彌散和機(jī)械擴(kuò)散共同作用的，定義水動(dòng)力彌散系數(shù)：水動(dòng)力彌散定律如下：式中：I單位時(shí)間內(nèi)通過(guò)面積的溶質(zhì)質(zhì)量；D水動(dòng)力彌散系數(shù)；Dc/ds濃度梯度。如果我們?nèi)》较蚺c流速方向一致，軸和軸與流速方向垂直，上式可用下式表示：三 對(duì)流彌散方程及其定解條件如圖，以滲流區(qū)內(nèi)任一點(diǎn)為中心，取一無(wú)限小的六面體單元，各邊長(zhǎng)為x、y、z，選擇x軸與P點(diǎn)處的平均流速方向一致。（即縱向彌散方向?yàn)檩S方向）在對(duì)流彌散問(wèn)題中，包括兩個(gè)子問(wèn)題：其一，溶質(zhì)隨地下水的流動(dòng)或流出單元體；其二，溶質(zhì)通過(guò)自身的彌散流入或流出單元體。水動(dòng)力彌散引起的物質(zhì)運(yùn)移：設(shè)，沿軸方向溶質(zhì)的質(zhì)量變化率為dIx/dx，如果假設(shè)Ix為在abcd面上，單位時(shí)間內(nèi)通過(guò)單位面積溶質(zhì)的質(zhì)量，那么，時(shí)間內(nèi)通過(guò)面流入單元體的溶質(zhì)質(zhì)量為：Ixnyzt因?yàn)檠剌S方向溶質(zhì)的質(zhì)量變化率為dIx/dx，所經(jīng)距離x后，變化了（dIx/dx）x，所以，t時(shí)間內(nèi)，通過(guò)abcd面流出單元體的溶質(zhì)質(zhì)量為：所以，沿x軸方向流入與流出單元體的溶質(zhì)質(zhì)量差為：同理，沿y軸方向和z軸方向流入與流出單元體的溶質(zhì)質(zhì)量差為：所以，通過(guò)彌散單元體內(nèi)溶質(zhì)質(zhì)量的變化為：隨地下水流的物質(zhì)運(yùn)移：設(shè)沿x軸方向，在abcd面地下水的流速為vx，則單位時(shí)間通過(guò)abcd面單位面積流入單元體的水量為vx11，流入的溶質(zhì)質(zhì)量為vxc：（為溶質(zhì)的濃度），那么，在t時(shí)間內(nèi)流入abcd面溶質(zhì)量為：vxcyzt設(shè)沿方向，通過(guò)單位面積溶質(zhì)質(zhì)量的變化率為：經(jīng)x距離的變化量為：在abcd面，單位時(shí)間單位面積流出單元體的溶質(zhì)質(zhì)量為：在t時(shí)間內(nèi)流出面的溶質(zhì)質(zhì)量為：所以，沿軸方向流入與流出單元體的溶質(zhì)的質(zhì)量差為：同理，沿軸和軸方向流入與流出單元體的溶質(zhì)質(zhì)量差為：所以，隨地下水流流入與流出單元體的溶質(zhì)質(zhì)量差為：t 時(shí)間內(nèi)，流入和流出單元體總的溶質(zhì)質(zhì)量差為：另外，設(shè)單元體內(nèi)溶質(zhì)濃度隨時(shí)間的變化率為dc/dt，那么，t時(shí)間內(nèi)單元體內(nèi)體積溶質(zhì)濃度變化量為：所以，t時(shí)間內(nèi)單元體內(nèi)溶質(zhì)質(zhì)量變化量為：上述二量應(yīng)相等，并消去xyzt，得：代入上式，并兩邊同除以，則得：式中，u實(shí)際流速。上式為對(duì)流彌散方程。如果，有其它源匯項(xiàng)時(shí)，并設(shè)單位時(shí)間單位體積含水層內(nèi)由源匯引起的溶質(zhì)質(zhì)量的變化量為f。在上式的左邊加一項(xiàng)f即可。關(guān)于溶質(zhì)運(yùn)移的數(shù)學(xué)模型除微分方程外，還應(yīng)有定解條件：初始條件：初始時(shí)刻的濃度分布。表達(dá)式如下：C(x，y，z，0)= C0(x，y，z)邊界條件：有兩類。一類是已知濃度的邊界條件，表示如下：式中：1表示一類邊界。另一類是通量邊界，即單位時(shí)間內(nèi)通過(guò)單位邊界面積的溶質(zhì)質(zhì)量已知。如：隔水邊界：補(bǔ)給邊界：所以，要確定一個(gè)水動(dòng)力彌散問(wèn)題的解，即求得濃度的分布，要給出下列信息：（1）微分方程；（2）研究空間區(qū)域和時(shí)間區(qū)域；（3）研究區(qū)域水頭場(chǎng)的分布；（4）有關(guān)參數(shù)，如彌散度和等；（4）定解條件。四 一維彌散問(wèn)題解設(shè)投放示蹤劑前，含含水層中示蹤劑的濃度為0，然后在河渠中連續(xù)注入濃度為C0的示蹤劑。在均勻流情況下，ux=u為常數(shù)（1）數(shù)學(xué)模型：定解問(wèn)題的解：其中當(dāng)時(shí)，上解可近似為：（2）利用實(shí)驗(yàn)資料求縱向彌散系數(shù)有一個(gè)觀測(cè)孔時(shí)，觀測(cè)孔距河渠距高已知，在觀孔中可測(cè)得不同時(shí)刻的，從而可求得。求參步驟： 據(jù)實(shí)驗(yàn)資料作Ci/C0t關(guān)系曲線，如圖 在圖上找出Ci/C0為0.84和0.16二點(diǎn)，并讀出其橫坐標(biāo)t0.84和t0.16。 代入下式求縱向彌散系數(shù)式中：u地下水實(shí)際流速3 海岸帶含水層中的咸淡水界面天然條件下，在海岸帶含水層中的地下水一般是流向海的，由于海水比淡水的比重大，海水體將位于淡水體的下方，是楔型，如圖。處于平衡狀態(tài)。當(dāng)在海岸邊抽取淡水時(shí)，這時(shí)淡水的水位下降，打破了原來(lái)的平衡，引起海水向內(nèi)陸的入侵，以達(dá)到新的平衡，這時(shí)界面向陸地推進(jìn)。該現(xiàn)象為海水入侵。海水與淡水是可以溶混的，由于水動(dòng)力彌散，在海水與淡水之間形成了一個(gè)過(guò)渡帶，在過(guò)渡帶中地下水的礦化度由小變大，直到海水礦化度。過(guò)渡帶的寬度在不同的地區(qū)，其寬度不同，當(dāng)其寬度較小，且與含水層的厚度相比較小時(shí)，可以以為海水與淡水之間是一突變界面；否則寬度較寬時(shí)，則作為水動(dòng)力彌散問(wèn)題加以研究。一 作突變界面處理靜止界面的近似解當(dāng)?shù)秃Ｋ幱谝环N平衡狀態(tài)時(shí)，界面是靜止的，如圖，假設(shè)淡水的容重為rf，海水的容重為rs。在界面上的點(diǎn)，受海水的壓力為：shs在界面上的點(diǎn)，受到淡水的壓力為：f（hs+hf）此二壓力應(yīng)相等：f（hs+hf）= shs解得：令則：hs=hf一般海水密度為1.025g/cm3，容重s=10045N/m3，淡水密度為1.000g/cm3，容重s=9800N/m3。代入可求得：=40 hs=40hf說(shuō)明：在離海岸任一距離上，穩(wěn)定界面在海面以下得深度為該處 淡水高出海面得倍。以上僅是一種近似解法。二 確定界面的形狀及海水入侵的范圍（一） 厚度固定的水平承壓含水層中的界面問(wèn)題水流是穩(wěn)定流，如圖，設(shè)，原點(diǎn)位于坡腳（點(diǎn)），x軸的正向指向海。含水層厚度為M，承壓水頭為H，假設(shè)地下水為水平流，忽略垂向分速度。由Darcy定律，有：K=Kf（Kf含水層中淡水的滲透系數(shù)）H=Hf（Hf含水層中淡水的水頭）由前面知，hs=d+h(x) hs=H d+h (x)= H兩邊對(duì)求導(dǎo)，得：（與無(wú)關(guān)）代入上式得：對(duì)上式整理得：兩邊積分得：當(dāng)x=0時(shí)，h(x)=M，代入上式得：代入上式得：上式表明界面得形狀是一條拋物線。利用此式可確定X處的h（x）。另外，由d+h (x)=H得：h (x)=H-d代入（1）式得：兩邊積分，得：當(dāng)x=0時(shí)，H=H0，代入上式，得：代入，得：此式表明承壓水含水層淡水的水頭面的形狀也是一條拋物線。當(dāng)x=L時(shí)，h=0，這時(shí)d+h=d=H H=d /代入上式，得：此式表示出海水入侵深度與流向海的淡水流量q0和界面坡腳以上測(cè)壓水頭H0之