摘要：以山東省鄒城市某巖溶水源地為例,研究數(shù)值模擬方法在巖溶地下水水源地保護(hù)區(qū)劃分中的應(yīng)用,以及水文地質(zhì)參數(shù)對(duì)于水源地保護(hù)區(qū)劃分結(jié)果的影響。基于GMS軟件構(gòu)建了研究區(qū)三維非穩(wěn)定流地下水?dāng)?shù)值模型,采用質(zhì)點(diǎn)追蹤技術(shù),通過(guò)反向追蹤示蹤粒子在運(yùn)移100d、1000d以及25a后的位置及遷移軌跡,對(duì)該巖溶水源地進(jìn)行保護(hù)區(qū)劃分;改變數(shù)值模型中的主要水文地質(zhì)參數(shù),以變化后的保護(hù)區(qū)面積作為衡量參數(shù)敏感度的指標(biāo),并計(jì)算各參數(shù)的敏感度系數(shù),研究水文地質(zhì)參數(shù)對(duì)于保護(hù)區(qū)劃分范圍的影響。結(jié)果表明:當(dāng)參數(shù)改變幅度在20%范圍內(nèi)時(shí),一層垂向滲透系數(shù)(VK1)的敏感度系數(shù)最大可達(dá)到2.63×10^-3,二層垂向滲透系數(shù)(VK2)的敏感度系數(shù)最大可達(dá)到3.64×10^-3;垂向滲透系數(shù)的敏感度明顯高于其他參數(shù),說(shuō)明垂向滲透系數(shù)對(duì)保護(hù)區(qū)劃分范圍的影響最大。因此,在應(yīng)用數(shù)值模擬法對(duì)巖溶地下水源地進(jìn)行保護(hù)區(qū)劃分時(shí),應(yīng)尤其注意模型中各含水層垂向滲透系數(shù)取值的準(zhǔn)確性和合理性。