孫少游，馬靜晨，劉全國(guó)，李世君,楊 楠，潘 懋

（1. 北京市地質(zhì)工程勘察院，北京 100048；2. 北京大學(xué)，北京 100087；3. University of British Columbia, Vancouver, Canada V6P4Y7。）

中國(guó)北方巖溶水脆弱性評(píng)價(jià)模型及改進(jìn)

孫少游1，馬靜晨1，劉全國(guó)1，李世君1,楊 楠2,3，潘 懋2

（1. 北京市地質(zhì)工程勘察院，北京 100048；2. 北京大學(xué)，北京 100087；3. University of British Columbia, Vancouver, Canada V6P4Y7。）

根據(jù)北方地區(qū)巖溶含水層的特征，基于改進(jìn)歐洲方法，利用ArcGIS構(gòu)建了包括指標(biāo)體系、權(quán)重、評(píng)分體系的北方地區(qū)巖溶含水層脆弱性評(píng)價(jià)模型。并運(yùn)用該模型，對(duì)北京房山巖溶含水層脆弱性進(jìn)行了評(píng)價(jià)。結(jié)果表明:巖溶水脆弱性較高的區(qū)域位于巖溶發(fā)育的地區(qū)尤其是溶洞周圍，而在范圍較小的研究區(qū)，降雨及植被覆蓋因素的影響較小。模擬結(jié)果顯示，北京房山巖溶區(qū)90%處于脆弱性較低的范圍，只有研究區(qū)西北部巖溶發(fā)育強(qiáng)烈的小部分區(qū)域脆弱性較高，建議加強(qiáng)進(jìn)一步調(diào)研或應(yīng)加強(qiáng)保護(hù)。

歐洲方法; ArcGIS; 巖溶水; 脆弱性; 房山

0 前言

巖溶水是地下水的主要賦存形式之一，肩負(fù)著世界大約25%人口的飲水任務(wù)（Derek et al，2007）。中國(guó)巖溶地下水分布面積約82.83×104km2，可開(kāi)采資源量約為870.02×108km3/a，占全國(guó)地下水可開(kāi)采資源總量的6.7%，也是某些地區(qū)主要的供水水源。賦存巖溶水的巖溶含水層地質(zhì)成因獨(dú)特，其一旦遭受污染，恢復(fù)治理將極其困難，特別在以巖溶水作為唯一水源的地區(qū)，其影響更是致命的。因此，對(duì)巖溶含水層進(jìn)行脆弱性評(píng)價(jià)并及時(shí)采取保護(hù)措施是十分必要的（趙玉國(guó)，2011）。

地下水脆弱性指的是污染物質(zhì)到達(dá)含水層系統(tǒng)某一特定位置的傾向性與可能性（鐘佐燊，2005；Derek et al，2007）。國(guó)內(nèi)外專家通過(guò)理論研究與實(shí)例驗(yàn)證已經(jīng)提出了多種地下水脆弱性評(píng)價(jià)的方法（彭穩(wěn)等，2010；趙玉國(guó)，2011；Goldseheider，2003），其中使用最廣泛的是DRASTIC，但它主要適用于孔隙含水層，對(duì)巖溶含水層脆弱性評(píng)價(jià)的適用性較差（鐘佐燊，2005；Daly et al，2002）。為此，又先后出現(xiàn)了針對(duì)巖溶含水層脆弱性評(píng)價(jià)的不同方法，如EPIK、越南模型、Slovenia模型、二元法、COP方法等（邢立亭等，2009；徐慧珍等，2006；章程等，2007），其中以歐洲模型考慮的指標(biāo)最為全面,評(píng)價(jià)體系也最為完善，是目前對(duì)巖溶含水層脆弱性指標(biāo)闡述最詳細(xì)的方法（彭穩(wěn)等，2010；張強(qiáng)等，2009），中國(guó)的蔣勇軍等人在對(duì)我國(guó)西南地區(qū)巖溶含水層脆弱性評(píng)價(jià)時(shí)，以歐洲方法為藍(lán)本進(jìn)行了改進(jìn)（張強(qiáng)，2009；楊桂芳，2002）。

近30年來(lái)，盡管地下水脆弱性研究方面的發(fā)展迅速，關(guān)于巖溶區(qū)地下水脆弱性研究卻進(jìn)展緩慢。本文改進(jìn)的巖溶水脆弱性評(píng)價(jià)模型基于較完善的歐洲方法，選取中國(guó)北方典型巖溶區(qū)北京市房山區(qū)為試點(diǎn)，結(jié)合《中華人民共和國(guó)環(huán)境保護(hù)行業(yè)標(biāo)準(zhǔn)——地下水脆弱性評(píng)價(jià)技術(shù)指南建議稿》，運(yùn)用ArcGIS作為工作環(huán)境，對(duì)評(píng)價(jià)因子進(jìn)行優(yōu)化，形成了一套適用于中國(guó)北方巖溶水脆弱性評(píng)價(jià)的程序化方法，提高了巖溶水脆弱性評(píng)價(jià)的準(zhǔn)確率。

1 歐洲方法

歐洲方法是基于“起源—路徑—目標(biāo)”模型，可應(yīng)用在地下水資源和水源地保護(hù)中，在此模型中將潛在污染物質(zhì)的釋放區(qū)視為“起源”；“目標(biāo)”是地下水面或水井、泉；“路徑”包括從起源到目標(biāo)之間的全部路徑（Goldseheider，2003）。歐洲方法的地下水脆弱性評(píng)價(jià)模型考慮四個(gè)因素：覆蓋層O（Overlying layers）、徑流特征C（Concentration of flow）、降水條件P（Precipitation regime）、巖溶系統(tǒng)發(fā)育K（Karstic network development）。

2 評(píng)價(jià)模型及改進(jìn)

本文的巖溶水脆弱性評(píng)價(jià)模型的指標(biāo)及原理如圖1所示。

以秦嶺淮河為界，中國(guó)南、北方巖溶在巖溶發(fā)育程度、巖溶水含水介質(zhì)、巖溶水動(dòng)態(tài)、賦存條件等方面都存在一系列差異，而國(guó)內(nèi)的《中華人民共和國(guó)環(huán)境保護(hù)行業(yè)標(biāo)準(zhǔn)——地下水脆弱性評(píng)價(jià)技術(shù)指南建議稿》并沒(méi)有對(duì)中國(guó)南、北方巖溶發(fā)育的特點(diǎn)做出詳細(xì)區(qū)分。另外，歐洲方法目前的主要試點(diǎn)在歐洲發(fā)達(dá)地區(qū)，這些地區(qū)基礎(chǔ)地質(zhì)資料豐富，巖溶勘察詳實(shí)，而在我國(guó)大部分巖溶區(qū)，勘察資料獲取存在一定的困難。基于上述問(wèn)題，評(píng)價(jià)模型個(gè)別因子的選取應(yīng)根據(jù)我國(guó)北方與歐洲方法的特點(diǎn)上有所區(qū)分。本文對(duì)歐洲模型進(jìn)行了如下改進(jìn)，詳見(jiàn)表1：

圖1 巖溶水脆弱性評(píng)價(jià)模型原理圖Fig.1 The Theory of Evaluation M odel of Vulnerability of Karst Water

（1）中國(guó)北方的年均降雨量在50～800mm左右，極端降雨（暴雨＞50mm/d，大雨＞25mm/d）次數(shù)極少，故本文將降雨條件P的次級(jí)因素概化為年均降雨量Ar（我國(guó)北方各地的降雨主要集中在夏季，故取值選取近30年來(lái)6—8月降雨量的均值）與年均極端降雨次數(shù)Se。P值越大，代表該地區(qū)脆弱性越低。

（2）考慮到中國(guó)北方巖溶發(fā)育的落水洞與消水河發(fā)育極少，但有少量溶洞發(fā)育，故將歐洲方法中的落水洞、消水河的影響合并為溶洞因素，且溶洞是否在降雨情況下處于活躍狀態(tài)仍需要鑒定。因此，用Tv因子來(lái)鑒定溶洞在全年中所起的作用。其中，距離溶洞越近的區(qū)域，其脆弱性越高。

（3）巖溶區(qū)與非巖溶區(qū)植被覆蓋情況及坡度對(duì)脆弱性的影響是相反的，因此，需以分水嶺為界區(qū)分快速補(bǔ)給區(qū)和緩慢補(bǔ)給區(qū)?？焖傺a(bǔ)給區(qū)是溶洞周邊區(qū)域，地形坡度越大、植被覆蓋越少，其脆弱性越高，緩慢補(bǔ)給區(qū)與此相反。

（4）因?yàn)镃、O、P、K各因素對(duì)脆弱性指數(shù)的貢獻(xiàn)程度是不同的，故利用專家打分法對(duì)各因素賦權(quán)重后進(jìn)行疊加計(jì)算。

3 改進(jìn)模型的應(yīng)用

房山位于北京市西南，其地理坐標(biāo)為東經(jīng)115°25′—116°15′，北緯39°30′—39°55′，屬于大清河水系，境內(nèi)地勢(shì)西北高，東南低。以鳳凰山—上寺嶺—霞云嶺分水嶺為界，北部為大石河，南部為拒馬河。北部大石河切割較淺，由霞云嶺向北經(jīng)薊縣系地層進(jìn)入寒武系、奧陶系灰?guī)r地層之中，沿岸現(xiàn)巖溶溶蝕作用形成的洼地，并受構(gòu)造作用影響形成佛子莊泉、河北泉等多處泉。南部低山丘陵區(qū)廣泛分布薊縣系霧迷山組碳酸鹽巖，受拒馬河切割影響，形成兩岸陡峭的峰叢、峰林等巖溶地貌，且?guī)r溶洞穴較發(fā)育，形成石花洞、云水洞等重點(diǎn)景觀巖溶洞穴。山區(qū)巖溶裂隙發(fā)育，地下水入滲補(bǔ)給條件好。

根據(jù)房山區(qū)的水文地質(zhì)資料，結(jié)合以上改進(jìn)模型進(jìn)行的巖溶水脆弱性評(píng)價(jià)結(jié)果如圖2所示。

表1 模型因子改進(jìn)前后對(duì)比表Tab.1 M odel factor comparison before and after improvement

圖2 房山區(qū)巖溶水脆弱性評(píng)價(jià)結(jié)果Fig.2 Evaluation of the vulnerability of karst water in fangshan

4 結(jié)果分析

根據(jù)此次收集的房山區(qū)基礎(chǔ)地質(zhì)、水位、降雨以及植被覆蓋于地形坡度等資料，利用本文基于ArcGIS改進(jìn)的北方巖溶水脆弱性評(píng)價(jià)模型得到的房山區(qū)巖溶水脆弱性評(píng)價(jià)結(jié)果顯示：

（1）降雨因素對(duì)整個(gè)研究區(qū)的影響完全相同，這是因?yàn)檠芯繀^(qū)面積較小，整個(gè)研究區(qū)的降雨情況基本相同。

（2）覆蓋層因素于土壤類型以及土地利用類型有關(guān)，在整個(gè)研究區(qū)的中部，覆蓋層產(chǎn)生的脆弱性較低；在西北和東南的部分區(qū)域，脆弱性評(píng)價(jià)最低，因?yàn)橥寥李愋图靶再|(zhì)較好，對(duì)地下水產(chǎn)生了較強(qiáng)的保護(hù)能力。

（3）落水洞于消水河對(duì)研究區(qū)的影響較大，因?yàn)槿芏粗苓厱?huì)產(chǎn)生坡度較大的地帶，這一地帶使得地表水可以快速流入溶洞，從而將地表污染物以及與地表水混合的其他污染物一并帶入。研究區(qū)中有小面積脆弱性高的地區(qū)主要位于溶洞附近，因污染物通過(guò)點(diǎn)狀快速入滲進(jìn)入溶洞，通過(guò)溶洞迅速進(jìn)入巖溶含水層。

（4）研究區(qū)4/5的面積處于脆弱性中等或低的地區(qū)，一次無(wú)需特別的調(diào)研與保護(hù)，但是在大型溶洞周邊，模擬產(chǎn)生的紅色區(qū)域指示脆弱性較高，因此需要特別注意或者進(jìn)行進(jìn)一步的調(diào)研。

5 結(jié)論

本文構(gòu)建了北方地區(qū)巖溶含水層脆弱性評(píng)價(jià)模型，以北京市房山區(qū)為例進(jìn)行了應(yīng)用及驗(yàn)證，表明構(gòu)建的模型基本符合實(shí)際情況，可用于北方地區(qū)巖溶含水層脆弱性評(píng)價(jià)。巖溶含水層脆弱性級(jí)別高的區(qū)域主要位于灰?guī)r裸露區(qū)、滲漏河道和泉源區(qū)，因此應(yīng)該進(jìn)行進(jìn)一步的調(diào)研或采取相應(yīng)的措施加強(qiáng)保護(hù)。

Derek C, W illiam s F，2007. Karst Hyd rogeology and Geomorphology [M]. Wiley Revised Edition: 551.

Goldseheider N, 2003. Karst groundwater vulnerability mapping: application of a new method in the Swabian Alb, Germany [J]. Hydrogeology Journal, (13): 555-564.

彭穩(wěn)，裴建國(guó)，2010. 巖溶含水層脆弱性評(píng)價(jià)方法探討[J]. 水資源保護(hù)，6(6)：9-12.

徐慧珍，高贊東，2006. 巖溶地區(qū)地下水防污性能評(píng)價(jià)——PI方法[J]. 新疆地質(zhì)，24(3)：318-321.

楊桂芳，姚長(zhǎng)宏，2002. 我國(guó)西南巖溶區(qū)地下水敏感性評(píng)價(jià)模型研究[J]. 科技進(jìn)展，25(2)：83-85.

章程，蔣勇軍，Michele，等，2007. 巖溶地下水脆弱性評(píng)價(jià)“二元法”及其在重慶金佛山的應(yīng)用[J]. 中國(guó)巖溶，26(4)：334-340.

鐘佐燊，2005. 地下水防污性能評(píng)價(jià)方法探討[J]. 地學(xué)前緣. 12 (S1)：3-11.

張強(qiáng)，蔣勇軍，林玉石，等，2009. 基于歐洲模型的巖溶地下水脆弱性風(fēng)險(xiǎn)性評(píng)價(jià)[J]. 人民長(zhǎng)江，40(13)：51-54.

趙玉國(guó)，2011. 基于GIS的巖溶地區(qū)地下水脆弱性評(píng)價(jià)——以重慶市老龍洞地下河流域?yàn)槔齕D]. 西南大學(xué).

Evaluation M odel and Its Im provement of Vulnerability of Karst Water in North China

SUN Shaoyou1, MA Jingchen1, LIU Quanguo1, LI Shijun1, YANG Nan2,3, PAN Mao2

(1. Beijing Institute of Geological & Prospecting Engineering, Beijing100048; 2. Peking University, Beijing 100871; 3. University of British Columbia, Vancouver, Canada, V6P4Y7)

An improved model was proposed based on the characteristics of karst aquifers in North China. This improved model, including index system, weight and scoring system, comes from the European Model and take advantage of ArcGIS as the eff cient tool to simulate the vulnerability of Fangshan karst aquifers in Beijing. The result shows that the regions w ith high vulnerability to karst water are in the areas where karst landform is relative well-developed especially around the karst caves. Rainfall and vegetation factors are less affected in smaller study areas. The simulation result show that 90% of Fangshan karst area in Beijing is in a low vulnerability area, only a small part of the region w ith strong karst development in the northwest of the study area has higher vulnerability. It is suggested to strengthen further research or strengthen protection in high vulnerable areas.

European Model; ArcGIS; Karst water; Vulnerability; Fangshan District

A

1007-1903（2017）02-0087-04

10.3969/j.issn.1007-1903.2017.02.017

北京巖溶水資源勘查評(píng)價(jià)工程（京發(fā)改[2011]1215號(hào)）

孫少游（1983- ），男，碩士，工程師，主要研究方向：地下水資源。E-mail：418171619@qq.com