1、吉林大學環(huán)境工程系吉林大學環(huán)境工程系20142014年年5 5月月水文地質(zhì)學系統(tǒng)結(jié)構(gòu)水文地質(zhì)學系統(tǒng)結(jié)構(gòu)含水層性質(zhì)水動力性質(zhì)孔隙介質(zhì)裂隙介質(zhì)巖溶介質(zhì)雙重介質(zhì)水文地球化學巖性構(gòu)造地貌地層補給排泄徑流水壓力重 力應 力環(huán)境水文地球化學機械過濾吸附交替作用溶解、沉淀作用生物作用對流、彌散作用生物地球化學作用污染質(zhì)在地下環(huán)境中的主要衰減作用第一講第一講 污染物的生物地球化學降解污染物的生物地球化學降解Biogeochemical Degradation of Pollutants氧化還原 反應方程 KJ/molCH2O 氧還原NO3-還原錳還原鐵還原SO42-還原產(chǎn)甲烷CH2O + O2 CO2 + H

2、2O5 CH2O+4NO3- 2N2+4HCO3- +CO2+3H2OCH2O+3CO2+H2O+2MnO2Mn2+4HCO3-CH2O+7CO2+4Fe(OH)3Fe2+8HCO3-+3H2OCH2O+SO42-H2S+2HCO3- CH2OCO2+CH4 -475-448-349-114-77-58 O2NO3- Mn4+Fe3+SO42-（1）地下環(huán)境中的主要氧化還原物質(zhì) 地下環(huán)境氧化還原系統(tǒng)中主要包括氣體（O2、N2、CH4、CO2）、溶解性組分（NO3-、NH4+、CH2O、Fe2+、Mn2+、SO42-、HS-、H+）、固體（FeOOH、MnO2、FeCO3、MnCO3）和存在于固

3、體中的離子交換態(tài)組分（Fe2+、Mn2+、NH4+）。1. 1. 污染物生物地球化學降解原理污染物生物地球化學降解原理1. 1. 污染物生物地球化學降解原理污染物生物地球化學降解原理1. 1. 污染物生物地球化學降解原理污染物生物地球化學降解原理包氣帶含水層地下水流向地下水位SO42-還原帶產(chǎn)甲烷帶NO3-還原帶氧還原帶鐵錳還原帶垃圾填埋場1. 1. 污染物生物地球化學降解原理污染物生物地球化學降解原理 （2）氧化還原帶形成的原理v 生物地球化學原理 微生物在利用O2、NO3-、鐵錳氧化物和氫氧化物、SO42-等最終電子受體并通過好氧呼吸、硝酸鹽還原、鐵錳還原、硫酸鹽還原和產(chǎn)甲烷等代謝途徑對有

4、機污染物進行降解，在此過程中，由于最終電子受體爭奪最終電子能力大小的不同，導致氧化還原反應發(fā)生的難易程度也不同，所以在污染暈中形成了一系列不同的氧化還原環(huán)境，或者氧化還原梯度，也就是順序氧化還原帶，即產(chǎn)甲烷帶、 硫酸鹽還原帶、鐵還原帶、錳還原帶、硝酸鹽還原帶、氧還原帶等。其實質(zhì)也是由于地下環(huán)境中的各種最終電子受體爭奪電子能力的不同所致。v 熱力學原理 O2、NO3- 、Mn4+、Fe3+和SO42-發(fā)生氧化還原反應所需要能量不同，造成了在不同的氧化還原環(huán)境中發(fā)生不同的氧化還原反應的難易程度和次序的不同，從而導致污染暈在時間和空間上依次出現(xiàn)不同的氧化還原環(huán)境，也就是順序氧化還原帶。1. 1. 污

5、染物生物地球化學降解原理污染物生物地球化學降解原理2. 2. 氧化還原帶中物質(zhì)變化特征氧化還原帶中物質(zhì)變化特征圖7.10 氧化還原帶動態(tài)發(fā)展理想模型abcedtabcdeMGZ/SRZIRZNRZ/ORZ低活性帶鐵氧化物ORZ低鐵帶低鐵還原性S豐鐵帶好氧污染源泄漏還原性物質(zhì)c氧化還原帶發(fā)展演化的理論模型3. 3. 氧化還原帶動態(tài)發(fā)展演化過程氧化還原帶動態(tài)發(fā)展演化過程第二講第二講 多相流多相流Multiphase Flow2 2. .多相流多相流（1）基本概念v非水相液體（Nonaqueous Phase Liquids, NAPL） 不與水發(fā)生混合的液體稱為NAPL;v重非水相液體（Dense

6、 Nonaqueous Phase Liquids，DNAPL） 比水重的非水相液體；v輕非水相液體（Light Nonaqueous Phase Liquids，LNAPLs） 比水輕的非水相液體v自由相（Free Phase） 在水動力條件下可流動的連續(xù)NAPL。v殘余相（Residual Phase）“困”在孔隙中，在水動力條件下不易流動的NAPL 。Mobile Phase NAPL or (Free-Phase NAPL) is a continuous mass of NAPL that can flow under a hydraulic gradientNAPL包氣帶和含水層中

7、相分布示意圖Residual Phase NAPL is trapped in the pore spaces between the soil particles, and cannot be easily moved hydraulically2 2. .多相流多相流WaterContinuousSlug ofNAPLVolume of NAPL VNAPLVolume of Open Pore Space (Effective Porosity) VPoreSolidIndividualNAPL BlobsPoreNAPLNVVSCan be estimated in the lab v

8、ia the following:Where:rb = soil bulk density g/cm3rn = NAPL bulk densityg/cm3n = porosityTPH = Total Petroleum Hydrocarbons mgNAPL/kgdry soil610nTPHSnbNrr2 2. .多相流多相流v 飽和度：指流體充填體積占總孔隙體積的比例。v 剩余飽和度：因毛細作用使其固定不變時的飽和度Partitioning Tracer Test for SNv Chemicals are introduced consisting of conservative a

9、nd partitioning tracers (PTT)v NAPL - Water Partition Coefficient KN.v KN = CN/Cw where CN is the conc. in NAPL and Cw is water conc.v At steady state, the effect is to delay or retard the PT transport rate according to a retardation factorv R = 1 + (KN SN)/(1 - SN) = tp/tn2 2. .多相流多相流Tp is travel t

10、ime for the Partitioning tracerTn is the travel time for the conservative tracer (Br-)Sorption is assumed insignificantPT compounds often used include alcohols such as ethanol， n-pentanol (1.4)，n-hexanol ，2,2 dimethyl-3-pentanol ，n-heptanolv Measured BTC for PTT with bromidev Time difference between

11、 different chemicals allows calculation of SN Partitioning Tracer Test for SNC/C0Effluent Volume (L)500100015002000Bromiden-pentanol 2,2 DMP2 2. .多相流多相流v 表面張力：液體與其他物質(zhì)接觸（如固體、不想混合的液體或氣體）將會產(chǎn)生界面能。這種能量是由于界面與內(nèi)部分子吸引力不同而引起的，這種現(xiàn)象稱為表面張力，定義為分子離開單位面積界面需要做的功，單位為達因/cm。2 2. .多相流多相流v 濕潤性：一種或另一種液體會優(yōu)先蔓延或濕潤整個固體表面。界面角和

12、表面張力之間的關系如下：應在密度大的液體中測出。如果小于90，液體L將優(yōu)先濕潤固體；如果大于90，液體G將優(yōu)先濕潤固體。也就是兩種液體爭奪一固體表面，其中一種占主導地位，覆蓋固體表面。WaterTesting fluid: water Background fluid: NAPL 110AIRSystem air:water air: NAPL water:NAPL air:NAPL:waterWetting Fluid water NAPL water waterorganicairNon-Wetting Fluid air air NAPL airGeneral Wetting Relat

13、ionships for Air, Water, NAPL:= Contact AngleNAPLSOLIDConclusion: Water isthe wetting fluidTesting fluid: NAPL Background fluid: airConclusion: NAPL isthe wetting fluidNAPLSOLIDNAPL is Non wetting fluid2 2. .多相流多相流v 毛細壓力（毛管壓力） 當兩種不相混合的流體接觸時，接觸面發(fā)生彎曲。每相靠近界面的孔隙壓力不一樣，這種壓力差就叫毛細管壓力。毛細管壓力在包氣帶為負值，當視為張力時為正值。

14、 球形空氣-水界面的曲率半徑r、毛細壓力Pc、表面張力之間的關系為：v 毛細管壓力與表面張力成正比，和曲率半徑成反比。曲率半徑取決于孔徑以及所填充溶液的量。v 毛細管壓力是兩種不想混溶的流體特性函數(shù)，即使在同一種介質(zhì)中，毛細管壓力也會由于水和NAPL比例的不同而不同。2 2. .多相流多相流v 相對滲透率 當兩種不相混溶的流體同時流動時，共同占有了有效孔隙。在給定飽和度的情況下，流體的固有滲透率與巖石總的固有滲透率之比。2 2. .多相流多相流100%0Water Saturation100%NAPL Saturation00100%NAPLWaterResidualSaturation of

15、 NAPLRelativePermeability80%60%40%20%SorIrreversibleSaturation of Waterv 束縛飽和度 在沒有任何附加水流動時的含水量。 在油-水兩相流系統(tǒng)，若含水量小于束縛飽和度，油能夠流動，水被毛細作用力控制；若含油量小于油的殘余飽和度時，水能夠流動，而油不能。II MixedFlowIDNAPL FlowIII WaterFlowIncreasing Water SaturationKr = relativepermeability1.000.800.600.400.200.100.080.060.040.020.011.000.80

16、0.600.400.200.100.080.060.040.020.01Increasing DNAPL SaturationKr = relativepermeability2 2. .多相流多相流Vadose ZoneResidualSaturation ofNAPL inVadose ZoneInfiltration, LeachingG round W ater FlowDissolved ContaminantPlume From NAPL ResidualSaturationVadoseZoneGroundwaterFlowInfiltration andLeachingFree-

17、Phase Product Layerin Saturated ZoneDissolvedContaminant PlumeResidualSaturationof LNAPLin soil LNAPL Release to Vadose ZoneLNAPL at the Water Table（2）NAPL遷移的概念模型2 2. .多相流多相流2 2. .多相流多相流LNAPL在含水層中的遷移2 2. .多相流多相流DNAPL在含水層中的遷移After Waterloo Centre for Groundwater Research, 1989.ResidualDNAPLDNAPL Pool

18、sSandLow PermeableStratigraphic UnitClayGroundwater FlowDissolvedContaminantPlume2 2. .多相流多相流透鏡體對DNAPL遷移的影響DNAPL在裂隙中的遷移SandVadose ZoneResidual DNAPLFractured Rock or Fractured ClayAfter Waterloo Centre for Ground Water Research, 1989.2 2. .多相流多相流Composite DNAPL SiteResidualDNAPLFractured ClayResidua

19、l DNAPLDNAPL PoolSandSandDissolved Contaminant PlumesImpermeable BoundaryAfter Waterloo Centre for Ground Water Research, 1989.2 2. .多相流多相流2 2. .多相流多相流NAPL migration in primary vs. secondary porosity(A) NAPL migration through primary porosity. (B) NAPL migration through secondary porosity features (

20、fractures, root holes, etc.).ABSourceSourceDNAPLDNAPL2 2. .多相流多相流LNAPL Moving UpgradientConfiningLayerGroundwater FlowLNAPL Plume Moving Upgradient to Stratigraphic TrapDissolved Contaminant PlumeSource Zone Forcing LNAPL into Confined Aquifer2 2. .多相流多相流Dissolving NAPLFresh WaterResidual NAPL Conta

21、ining Soluble HydrocarbonsContaminatedGroundwater FormingDissolved HydrocarbonPlume2 2. .多相流多相流Apparent LNAPL Thickness True LNAPL ThicknessMonitoring WellVadose ZoneCapillary FringeSaturatedZoneApparent LNAPL ThicknesshfhwrwrLNAPLrLNAPLHf：LNAPL漂浮厚度Hw：井中LNAPL厚度w：水的密度LNAPL：LNAPL的密度v Thickness in scre

22、ened well casing only indicator of actual thickness in the unitv Thickness in the well may be 2 to 10 times larger than actual thickness of mobile LNAPL in the surrounding aquiferv No methods exist to improve these estimates2 2. .多相流多相流DNAPL Detection ProblemsScreen Too High Above Confining Unit No

23、Accumulation Screen Too Deep Into Confining Unit Too Much Accumulation Screen Penetrates Clay LensAccumulation in Wrong Place SaturatedZoneVadoseZoneSaturatedZoneClayClayClayClayDNAPLDNAPLDNAPLDNAPL2 2. .多相流多相流第三講第三講 飽和介質(zhì)中的溶質(zhì)運移飽和介質(zhì)中的溶質(zhì)運移Mass Transport in Saturated Media（1）濃度梯度作用下的遷移 水中的溶質(zhì)從濃度高的地方向濃度低

24、的地方遷移，這個過程叫做分子擴散，或擴散。流體擴散的量與濃度梯度成正比，稱為Fick第一定律。F=-Dd（dC/dx） F：單位時間通過單位面積的溶質(zhì)質(zhì)量； Dd：擴散系數(shù)（L2/T） C：溶質(zhì)濃度（M/L3)負號表示運動沿濃度減小的方向進行。不同物質(zhì)的Dd值可以查到。對于濃度隨時間的系統(tǒng)，適用于Fick第二定律，其一維表達形式為：1 1. .濃度梯度作用下的遷移濃度梯度作用下的遷移 ：為濃度隨時間的變化。 在多孔介質(zhì)中，由于離子在礦物顆粒周圍擴散時需要經(jīng)過更長的路徑，因此，引入有效擴散系數(shù)D*。D*=Dd Ci是t時刻距離污染源x處的濃度；C0是初始濃度，erfc是余誤差函數(shù)。是一個與彎曲度

25、有關的系數(shù)，通常小于1，通過擴散試驗獲得。濃度分布服從正態(tài)分布或高斯分布，還可用統(tǒng)計特征值來描述。Fick第二定律一維表達式的解：1 1. .濃度梯度作用下的遷移濃度梯度作用下的遷移 地下水運動中攜帶溶解性固體使之發(fā)生遷移，這個過程叫做對流遷移。Vx為平均線速度；K為滲透系數(shù)；ne為有效孔隙度；dh/dl為水力梯度。 由于對流作用，一維質(zhì)量通量Fx為：Fx=vxneC2 2. .對流遷移對流遷移 地下水以大于和小于平均線速度的速率運動著。原因有：v孔隙中心處的流速比孔隙壁處大；v在直線距離相同時，多孔介質(zhì)中有的流體質(zhì)點會比其他質(zhì)點沿著更長的路徑運動；v孔隙的大小不一，大孔流速快。 若含有溶質(zhì)的

26、地下水以相同的速度運動，那么他們就會取代不含溶質(zhì)的水，并在兩種水體之間產(chǎn)生一個突變界面。然而，由于侵入的含有溶質(zhì)的水不是以相同速度運動的，所以這兩種水會沿著流動路徑發(fā)生混合，這種混合就稱為機械彌散，他導致在水流的推進邊緣被稀釋。沿著流動路徑方向發(fā)生的混合稱為縱向彌散；溶質(zhì)鋒面沿著垂直于流動方向發(fā)生的混合，稱為橫向彌散。3 3. .機械彌散機械彌散 在地下水流動過程中，分子擴散過程不能從機械彌散中分離出來，將二者結(jié)合起來定義一個參數(shù)稱為水動力彌散系數(shù)D，即DL=Lvi+D*DT=Tvi+D*DL：平行于水流運動方向的水動力彌散系數(shù)（縱向）；DT：垂直于水流運動主方向的水動力彌散系數(shù)（橫向）；L：

27、縱向彌散度；T：橫向彌散度。4 4. .水動力彌散水動力彌散FzFzz dzFyFyydyFxFxx dx（4）對流-彌散方程推導T：橫向彌散度。2 2. .飽和介質(zhì)中的物質(zhì)遷移飽和介質(zhì)中的物質(zhì)遷移FxFyFz （1）吸附過程 吸附過程包括吸附、化學吸附、吸收和離子交換。v吸附：包含溶質(zhì)附著到固體表面的過程。v陽離子交換：陽離子被吸附到帶負電的粘土礦物附近并通過靜電力附著在那里；v化學吸附：溶質(zhì)因化學反應而與沉積層、土壤或巖石表面結(jié)合。v吸收：若含水層顆粒是多孔的，且使溶質(zhì)能夠通過擴散進入顆粒并被吸附到其表面上。3 3. .溶質(zhì)轉(zhuǎn)化、阻滯和衰減溶質(zhì)轉(zhuǎn)化、阻滯和衰減（2）表面反應平衡1）線性等溫

28、吸附 溶質(zhì)的濃度C和固體吸附的溶質(zhì)的量C*之間存在一個直接的線性關系，即：C*=KdCC*：固體單位容重所吸附溶質(zhì)質(zhì)量（mg/kg);C：吸附平衡時溶液中溶質(zhì)的濃度（mg/L）；Kd為分配系數(shù)（L/kg），等于等溫吸附線的斜率。3 3. .溶質(zhì)轉(zhuǎn)化、阻滯和衰減溶質(zhì)轉(zhuǎn)化、阻滯和衰減 整理后，得： 阻滯因子Rf由下式給出：3 3. .溶質(zhì)轉(zhuǎn)化、阻滯和衰減溶質(zhì)轉(zhuǎn)化、阻滯和衰減 如果地下水平均流速為vx，則濃度為初始濃度一半的溶質(zhì)鋒面的平均速度vc為： 線性等溫吸附模型有兩個缺陷：v沒有限制固體吸附溶質(zhì)的量；v點少時，易將C* - C曲線形式誤解為直線。2）Freundlich等溫吸附3 3. .溶質(zhì)

29、轉(zhuǎn)化、阻滯和衰減溶質(zhì)轉(zhuǎn)化、阻滯和衰減 對此式進行線性化： logC*=logK+NlogCC*=KCN 整理后得：阻滯因子Rff由下式給出： 若N1，鋒面不斷擴展；若N1，鋒面會自銳；若N=1，就是線性等溫吸附3 3. .溶質(zhì)轉(zhuǎn)化、阻滯和衰減溶質(zhì)轉(zhuǎn)化、阻滯和衰減 是與結(jié)合能有關的吸附常數(shù)（L/mg）是固體所能吸附的溶質(zhì)的最大量（mg/kg）Langmuir等溫吸附也可表示為： 整理后得： 阻滯因子Rfl由下式給出： 3）Langmuir等溫吸附 第四講第四講 地下水污染數(shù)值模擬地下水污染數(shù)值模擬Numerical Simulation of Groundwater Pollution 數(shù)值模擬

30、的目的：v測試一種假設，或提高對某一給定含水層系統(tǒng)的認識；v了解物理、化學或生物過程；v設計地下水污染修復系統(tǒng)；v預測未來狀況或模擬壓力對地下水的影響；v資源管理4 4. .地下水污染數(shù)值模擬地下水污染數(shù)值模擬一、建模的基本信息（一）模型范圍及邊界 一般以相對獨立的水文地質(zhì)單元為模擬區(qū)范圍，便于準確刻畫自然邊界，提高精度。實際中，由于工作范圍、資料精度以及模擬所限，模擬去可能不是一個完整的地下水系統(tǒng)，這種情況下就要人為劃定模擬邊界。1、對以河流河流或其它地表水體為邊界的情況，可根據(jù)河流與地下水的關系，設置為一類（有直接水力聯(lián)系，已知河流水位）、二類（河流與含水層脫節(jié)，并已知河流滲入補給量）或三

31、類邊界條件（河流下有弱透水層，已知河流水位）。2、以地下水分水嶺、不透水斷層為邊界，可以定義為隔水邊界，即二類邊界，邊界量設為零。3、以山區(qū)-平原交界線為邊界的情況，一般可定義為二類邊界，山區(qū)向平原的流入量可根據(jù)地下水流場或山區(qū)匯水面積等方法確定。4、對于人為邊界，可設置為通用水頭邊界5、泉、排水溝、內(nèi)河河流邊界等可直接處理成含水層的補給量。4 4. .地下水污染數(shù)值模擬地下水污染數(shù)值模擬（二）水文地質(zhì)結(jié)構(gòu) 根據(jù)鉆孔巖性、水文地質(zhì)剖面圖、地下水開采情況（深度）以及研究目的，劃分模擬區(qū)含水巖組、弱透水層、和隔水層，再根據(jù)地下水流動特征，確定水文地質(zhì)結(jié)構(gòu)是單層結(jié)構(gòu)還是多層結(jié)構(gòu)。各層的頂板可用其標高等值線來表示，各層厚度可用其厚度等值線表示。（三）水文地質(zhì)參數(shù) 水文地質(zhì)參數(shù)主要包括兩類：一類是為求含水層補排量所需的參數(shù)，如降水入滲系數(shù)、河流滲漏系數(shù)、蒸發(fā)系數(shù)等；另一類表明含水層介質(zhì)性質(zhì)參數(shù)，如水動力參數(shù)、有關的溶質(zhì)運移參數(shù)、熱運動參數(shù)等。 常用的水文地質(zhì)參數(shù)為含水層的滲透系數(shù)、承壓含水層的導水系數(shù)、潛水含水層的給水度、承壓含水層的儲水系數(shù)和釋水系數(shù)。（四）初始流場 對于模擬非穩(wěn)定流，需要提供模擬初始時刻的地下水流場，地下水水位觀測點太少或分布不均，無法控制模擬區(qū)地下水流場。4 4. .地下水污染數(shù)值模擬地下水污染