采用敏感性分析和最差情況分析進(jìn)行地下水污染風(fēng)險(xiǎn)-wenkub
2023-07-14 17:25:57 本頁面
【正文】 （Calow，1998）。有幾種解決參數(shù)不確定性的方法。盡管蒙特卡羅模擬方法的功能強(qiáng)大，而且漸近收斂，但是計(jì)算效率卻很差。通過敏感性分析，可以確定由于輸入變量和參數(shù)造成的輸出變化，這是一種檢驗(yàn)?zāi)Ｐ椭休敵鲎兞繉?duì)輸入變量的靈敏性的分析方法。zalka市擁有25400人口，位于匈牙利東部，靠近羅馬尼亞和烏克蘭邊界。M225。另一個(gè)地下水污染源是在1971～1997年間利用的污水氧化池（Nauner，2000）。工業(yè)活動(dòng)是第五個(gè)地下水污染源。zalka位于匈牙利平原，是帕諾尼亞（Pannonian）山間盆地的一部分。上更新統(tǒng)由中細(xì)砂和粉砂組成，也可以作為含水層，但滲透性比上更新統(tǒng)要差。表1 更新統(tǒng)的6個(gè)水文地質(zhì)單元描述水文地質(zhì)單元平均厚度（m）描述地層第一層65薄砂、粉砂和粘土層上更新統(tǒng)第二層25粘土或粘質(zhì)砂層中更新統(tǒng)第三層7砂層第四層40粘土或粘質(zhì)砂層第五層20粗砂礫層和粘土互層的各向異性單元下更新統(tǒng)第六層100厚砂礫和粘土互層三、地下水流動(dòng)模型根據(jù)MODFLOW來求解微分方程。按照每月的抽水資料，確定模型中從井場(chǎng)抽取的地下水量。為了進(jìn)行計(jì)算，每個(gè)單元的長(zhǎng)寬比不超過10。采用Beyer公式和Zamarin公式對(duì)第六層的6個(gè)土樣進(jìn)行了粒度分析，這兩個(gè)公式是分析粒度與傳導(dǎo)系數(shù)之間關(guān)系的經(jīng)驗(yàn)公式。第五層和第六層是滲透性最強(qiáng)的更新統(tǒng)地層；第二層和第四層是滲透性最差的更新統(tǒng)地層，它們形成了地下水流動(dòng)的天然屏障。表3 采用“試錯(cuò)法”校正前后的水力傳導(dǎo)系數(shù)初始值校正值Kh110－5 m/s區(qū)1：10－4 m/s區(qū)2：10－5 m/s區(qū)3：10－5 m/s區(qū)4：10－5 m/sKv110－6 m/s區(qū)1：10－5 m/s區(qū)2：10－6 m/s區(qū)3：10－6 m/s區(qū)4：10－6 m/sKh210－7 m/s10－8 m/sKv210－8 m/s10－9 m/sKh310－5 m/s10－5 m/sKv310－6 m/s10－6 m/sKh410－7 m/s10－8 m/sKv410－8 m/s10－9 m/sKh510－4 m/s10－4 m/sKv510－5 m/s10－5 m/sKh610－4 m/s10－5 m/sKv610－5 m/s10－6 m/s有效入滲率30 mm/年25 mm/年采用PEST進(jìn)行自動(dòng)校正，這是一個(gè)參數(shù)估計(jì)程序。在本研究過程中，通過自動(dòng)參數(shù)估計(jì)與試錯(cuò)校正獲得的平均絕對(duì)誤差相同。應(yīng)當(dāng)注意到，出于計(jì)算的原因，模型中Kraszna河的深度要大于實(shí)際觀測(cè)值。采用MODFLOW的區(qū)域預(yù)算模塊計(jì)算穩(wěn)定態(tài)的水均衡，可以了解該區(qū)不同地層、河流、入滲、抽水和邊界之間能量的相互作用。第三層是細(xì)砂層，從第二層獲得的水量幾乎全部流入第四層，沿側(cè)邊界輸入的水量和輸出的水量相差不大。這說明通過抽水井獲取的水量中，有一部分是來自垂直補(bǔ)給，僅通過側(cè)邊界補(bǔ)給，不能滿足6181m3/天的抽水量。（一）邊界條件在運(yùn)移模型的邊界，假定濃度梯度。確定彌散性時(shí)要復(fù)雜一些，彌散值與測(cè)試或觀測(cè)范圍有關(guān)（Zheng and Bennett，1995）。由于資料不足，選擇橫向彌散度比縱向彌散度小一
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1