【正文】 入第三層。根據(jù)“試錯(cuò)法”獲得的水力傳導(dǎo)系數(shù)見表3。模型的東部邊界是一條河流，這條河流的坡度為13cm/km，水力傳導(dǎo)系數(shù)約為10－6m/s。與污染處理系統(tǒng)沒(méi)有聯(lián)系的住宅區(qū)是第四個(gè)地下水污染源，這些住宅區(qū)的污染坑沒(méi)有用混凝土加固，因此，污水容易到達(dá)地下，特別是在地下水位較高時(shí)。解決不確定性問(wèn)題最直接的方法是敏感性分析和（或）最差情況分析。在該研究區(qū)存在幾個(gè)地下水污染源，而且位于飲用水井附近。盡管蒙特卡羅模擬方法的功能強(qiáng)大，而且漸近收斂，但是計(jì)算效率卻很差。另一個(gè)地下水污染源是在1971～1997年間利用的污水氧化池（Nauner，2000）。表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來(lái)求解微分方程。第五層和第六層是滲透性最強(qiáng)的更新統(tǒng)地層；第二層和第四層是滲透性最差的更新統(tǒng)地層，它們形成了地下水流動(dòng)的天然屏障。采用MODFLOW的區(qū)域預(yù)算模塊計(jì)算穩(wěn)定態(tài)的水均衡，可以了解該區(qū)不同地層、河流、入滲、抽水和邊界之間能量的相互作用。確定彌散性時(shí)要復(fù)雜一些，彌散值與測(cè)試或觀測(cè)范圍有關(guān)（Zheng and Bennett，1995）。首先，對(duì)第一層到第二層和第五層到第六層的邊界條件、水力傳導(dǎo)系數(shù)、河流參數(shù)以及沿垂直方向的水流量進(jìn)行分析。如前所述，第一層邊界的測(cè)壓水頭增加了2m，第五層和第六層的邊界條件減小了2m，粘土層的水力傳導(dǎo)系數(shù)增加了10倍，這就意味著水平和垂直水力傳導(dǎo)系數(shù)分別為110－6m/s和110－7m/s。通過(guò)敏感性分析，研究了邊界條件和參數(shù)值對(duì)污染物向下游遷移的影響；通過(guò)最差情況分析，確定污染物向下所能達(dá)到的最大深度。有效孔隙度對(duì)濃度分布情況影響不大，而彌散度增加會(huì)使污染源以下污染物的濃度增加，而且會(huì)造成污染物沿與水流垂直的方向運(yùn)移，使污染物到達(dá)較深的地方。計(jì)算主要污染源下125和35m深度的污染物濃度，結(jié)果表明，在主要污染源以下15m，污染物濃度低于污染源以下5m處污染物濃度的幾倍，深度為25m和35m，污染物的濃度相當(dāng)?shù)?。沒(méi)有考慮由于擴(kuò)散造成的運(yùn)移，因?yàn)樵诟邼B透性的環(huán)境中，擴(kuò)散引起的運(yùn)移量可以忽略不計(jì)（Carges and Baehr，1998）。（五）結(jié) 果計(jì)算東西向剖面的測(cè)壓水位，結(jié)果表明：第一層在河流西邊的位置上，地下水補(bǔ)給河流。在第一層進(jìn)行了12次試驗(yàn)，在第五層進(jìn)行了3次試驗(yàn)，在第六層進(jìn)行了8次試驗(yàn)。下更新統(tǒng)的沉積物主要由砂礫組成，是該區(qū)的含水層系統(tǒng)，是重要的飲用水源。t233。由于資料的缺乏和模型參數(shù)的異質(zhì)性，往往會(huì)造成模擬結(jié)果的不確定性。計(jì)算機(jī)模型是常用的進(jìn)行地下水流動(dòng)和污染預(yù)測(cè)的工具。M225。研究區(qū)內(nèi)更新統(tǒng)沉積物厚度約為260m，下更新統(tǒng)厚度為110m，中更新統(tǒng)厚度為90m，上更新統(tǒng)厚度為60m。（三）水力傳導(dǎo)系數(shù)利用抽水試驗(yàn)、排泄和水位下降資料以及粒度分布資料，確定研究區(qū)不同地層的水力傳導(dǎo)系數(shù)。自動(dòng)參數(shù)估計(jì)并沒(méi)有減小誤差，因此，在進(jìn)一步分析工作中，采用的是試錯(cuò)校正。四、運(yùn)移模型采用兩種不同的方法，對(duì)運(yùn)移情況進(jìn)行模擬，這兩種方法分別是采用MODPATH進(jìn)行顆粒示蹤（Pollo