焦恒，王若帆，江峰，易瑞

貴州省開(kāi)陽(yáng)縣光生壩地下河系統(tǒng)地下水脆弱性評(píng)價(jià)

焦恒，王若帆，江峰，易瑞

（貴州省地質(zhì)礦產(chǎn)勘查開(kāi)發(fā)局114地質(zhì)大隊(duì)，貴州 遵義 563000）

本文結(jié)合開(kāi)陽(yáng)縣地下水系統(tǒng)情況，以地下河為單元，采用層次分析法，取固有脆弱性以巖石地層（R）、含水介質(zhì)類(lèi)型（W）、巖溶類(lèi)型（K）、巖溶發(fā)育程度（D）、地形地貌（T）、土層厚度（S）、植被覆蓋率（V）因子進(jìn)行分值量化;在固有脆弱性評(píng)價(jià)的基礎(chǔ)上增加人為工程活動(dòng)程度（E）、污染源（P）、人口密集程度（C）作為評(píng)價(jià)因子，對(duì)該區(qū)地下水的固有脆弱性和特殊脆弱性進(jìn)行了綜合評(píng)價(jià)。評(píng)價(jià)結(jié)果為該區(qū)域地下水開(kāi)發(fā)利用和保護(hù)提供有價(jià)值的參考依據(jù)。

地下水脆弱性；層次分析法；評(píng)價(jià)因子；開(kāi)陽(yáng)縣

開(kāi)陽(yáng)縣位于黔中腹地，總面積2026km2。具有“喀斯特生態(tài)世界公園”的美譽(yù)，是全國(guó)首個(gè)循環(huán)經(jīng)濟(jì)磷煤化工生態(tài)工業(yè)示范基地縣和貴州省經(jīng)濟(jì)強(qiáng)縣。2019年開(kāi)陽(yáng)縣政府堅(jiān)持生態(tài)優(yōu)先，努力實(shí)現(xiàn)百姓富、生態(tài)美的統(tǒng)一目標(biāo)。

表1 各評(píng)價(jià)因子評(píng)分表

近年來(lái)地下水問(wèn)題嚴(yán)重制約該區(qū)域經(jīng)濟(jì)發(fā)展，開(kāi)展地下水脆弱性評(píng)價(jià)，進(jìn)行地下水脆弱性分區(qū)，可供決策管理者制訂地下水管理戰(zhàn)略參考。

1 研究區(qū)的地質(zhì)背景

光生壩地下河系統(tǒng)（王若帆等，2020）位于貴州省開(kāi)陽(yáng)縣，地下河系統(tǒng)西側(cè)、南側(cè)及北側(cè)均以地表分水嶺為界，南東側(cè)以寒武系第二統(tǒng)金頂山組（∈2）頂界為界，北東側(cè)以地表分水嶺為界，構(gòu)成相對(duì)獨(dú)立地下河系統(tǒng)（圖1）。

地質(zhì)構(gòu)造主要為近南北向石頭寨斷層及受構(gòu)造作用發(fā)育的橫張、縱張裂隙發(fā)育。

碳酸鹽巖含水巖組包括：Pt33b∈1、∈23、∈2、∈3、∈3、∈3-4O1、P2。其中：Pt33b∈1、∈3、∈3、∈3-4O1巖性以白云巖為主，含水介質(zhì)發(fā)育溶孔、溶隙，含水較均勻，富水性中等；∈23、∈2、P2巖性以石灰?guī)r為主，含水介質(zhì)發(fā)育裂隙、溶洞、管道為主，含水極不均勻，富水性強(qiáng)?；鶐r裂隙含水巖組包括：∈1-2、∈21、∈22、∈2、P2，含水介質(zhì)以風(fēng)化裂隙、構(gòu)造裂隙為主，含水較均勻，富水性弱－貧乏，分布在補(bǔ)給區(qū)。

系統(tǒng)內(nèi)總體地勢(shì)南西高北東低，地貌類(lèi)型為峰叢溝谷，地形起伏大，相對(duì)高差100～300m，溝谷多呈南西-北東展布。該系統(tǒng)為裸露型開(kāi)放式地下河系統(tǒng)，系統(tǒng)內(nèi)地下水主要接受南西、西側(cè)一帶大氣降水匯流成地表溪溝注入伏流入口（K245、K246）、落水洞（K233、K236）補(bǔ)給地下水，受地形展布方向及近南北向橫張節(jié)理控制，地下水總體由西向東沿地下河管道徑流，受北西-南東向溝谷切割于光生壩（S1192）一帶以地下河出口的形式出露于地表。

圖1 光生壩地下河系統(tǒng)水文地質(zhì)略圖

1.地下河管道及出口（左：編號(hào)，右上：流量L/s：右下：高程m）；2.受污染的地下河管道及出口；3.下降泉；4.表層巖溶泉；5.地下水開(kāi)采井；6.熱響應(yīng)試驗(yàn)孔；7.伏流入口；8.有水落水洞；9.落水洞；10.巖溶洼地；11.污水處理廠；12.斷層；13.地質(zhì)界線及地層代號(hào)；14.富水塊段及編號(hào)；15.地下水流向

2 研究區(qū)環(huán)境水文地質(zhì)條件

2.1 地球化學(xué)背景

地下水主要賦存于Pt33b∈1、∈23、∈2、∈3、∈3、∈3-4O1、P2地層中，巖性以灰?guī)r（CaCO3）、白云巖（MgCO3）為主，SiO2、Al2O3等成分含量極少；地下水化學(xué)組份在圖內(nèi)相對(duì)單一，陽(yáng)離子以Ca2+、Mg2+，陰離子以HCO3-為主；其他重金屬元素及有毒、有害組分含量極少（吉勤克補(bǔ)子和李強(qiáng)，2020）。

2.2 環(huán)境地質(zhì)背景

調(diào)查發(fā)現(xiàn)南西側(cè)開(kāi)陽(yáng)縣城城區(qū)生活污水通過(guò)排污管道送至H139號(hào)及H144號(hào)污水處理廠處理后排放。訪問(wèn)H139號(hào)污水處理廠處理能力有限，暴雨季節(jié)，雨、污水混合通過(guò)K245、K246號(hào)伏流入口污染該區(qū)地下水。

3 地下水脆弱性評(píng)價(jià)

按中國(guó)地質(zhì)調(diào)查局發(fā)布《地下水脆弱性評(píng)價(jià)技術(shù)要求》（GWI-D3），地下水脆弱性是指地下水系統(tǒng)對(duì)人類(lèi)（或）自然的敏感性。一般分為固有脆弱性和特殊脆弱性（孫才志和潘俊，1999；孫才志和林山杉，2001；雷靜，2002；范琦等，2007；郭曉靜等，2010）。

固有脆弱性是指在天然狀態(tài)下含水層對(duì)污染所表現(xiàn)內(nèi)部固有的敏感性。不考慮污染源或污染物的性質(zhì)和類(lèi)型，是靜態(tài)、不可變和人為不可控的。

表2 層次分析法（AHP）評(píng)價(jià)表

特殊脆弱性是指含水層對(duì)特定的污染物或人類(lèi)活動(dòng)所表現(xiàn)的敏感屬性。它與污染源和人類(lèi)活動(dòng)有段，是動(dòng)態(tài)、可變和人為可控制的。

3.1 評(píng)價(jià)方法

1）固有脆弱性根據(jù)系統(tǒng)內(nèi)巖性分布情況及水文地質(zhì)特征對(duì)系統(tǒng)進(jìn)行分區(qū)，以巖石地層（R）、含水介質(zhì)類(lèi)型（W）、巖溶類(lèi)型（K）、巖溶發(fā)育程度（D）、地形地貌（T）、土層厚度（S）、植被覆蓋率（V）進(jìn)行分值量化。特殊脆弱性評(píng)價(jià)是在固有脆弱性評(píng)價(jià)的基礎(chǔ)上增加人為工程活動(dòng)程度（E）、污染源（P）、人口密集程度（C）作為評(píng)價(jià)因子（江峰等，2021）。

2）采用層次分析法（AHP）確定評(píng)價(jià)因子的權(quán)重。

3）建立符合工作區(qū)實(shí)際的地下水環(huán)境脆弱性模型。模型采用下式：

表3 固有脆弱性層次分析法（AHP）各評(píng)價(jià)因子權(quán)重分配表

表4 特殊脆弱性層次分析法（AHP）各評(píng)價(jià)因子權(quán)重分配表

評(píng)價(jià)因子RWKDTSVEPC 權(quán)重0.100.090.080.120.070.100.070.120.130.11

DI為脆弱性綜合指數(shù)；R、W、K等為上述評(píng)價(jià)因子；角標(biāo)為各評(píng)價(jià)因子權(quán)重；角標(biāo)為各評(píng)價(jià)因子量化分值。

3.2 評(píng)價(jià)因子量化

將巖石地層（R）、含水介質(zhì)類(lèi)型（W）、巖溶類(lèi)型（K）、巖溶發(fā)育程度（D）、地形地貌（T）、土層厚度（S）、植被覆蓋率（V）、人為工程活動(dòng)程度（E）、污染源（P）、人口密集程度（C）評(píng)價(jià)因子，按照[1,5]之間的整數(shù)取分值。詳見(jiàn)表1。

3.3 評(píng)價(jià)指標(biāo)權(quán)重計(jì)算

本次評(píng)價(jià)采用層次分析法（AHP）來(lái)確定評(píng)價(jià)指標(biāo)的權(quán)重（表2）。

先將矩陣A按列求和并進(jìn)行歸一化處理，得到矩陣B；再將矩陣B按行求和得到矩陣V，將矩陣V進(jìn)行歸一化處理得到矩陣W，矩陣W既為各影響因子權(quán)重，見(jiàn)表3、表4。

判斷矩陣的一致性檢驗(yàn)：所謂一致性是指判斷思維的邏輯一致性。采用一致性比率（式2）進(jìn)行驗(yàn)證。

表5 地下水環(huán)境脆弱性等級(jí)劃分一覽表

表6 地下水特殊脆弱性等級(jí)劃分一覽表

CR為一致性比率；CI為一致性指標(biāo)（采用式（3）計(jì)算）；RI為隨機(jī)一次性指標(biāo)。

CI為一致性指標(biāo)；λ為最大特征根的特征向量；為階數(shù)。

根據(jù)計(jì)算得：當(dāng)CR＜0.1時(shí)，A的不一致程度在容許范圍內(nèi)，稀釋可用A的特征向量作為權(quán)向量。

3.4 地下水環(huán)境脆弱性等級(jí)劃分

本次參照《地下水脆弱性評(píng)價(jià)技術(shù)要求》（GWI-D3）對(duì)地下水的脆弱性劃分結(jié)果結(jié)合工作區(qū)實(shí)際情況，將工作區(qū)地下水環(huán)境脆弱性劃分為五個(gè)等級(jí)，即：脆弱性很低（很難污染）、脆弱性較低（較難污染）、脆弱性較高低（較易污染）、脆弱性很高（很易污染）、脆弱性極高（極易污染），見(jiàn)表5、表6。

3.5 地下水固有脆弱性評(píng)價(jià)

1）系統(tǒng)內(nèi)亞區(qū)分區(qū)

根據(jù)系統(tǒng)內(nèi)的巖性特征、水文地質(zhì)條件等，將光生壩地下河系統(tǒng)分為3個(gè)亞區(qū)，即：石灰?guī)r分布區(qū)（A1）、白云巖分布區(qū)（A2）、碎屑巖分布區(qū)（A3），分區(qū)依據(jù)見(jiàn)表7。

表7 光生壩地下河系統(tǒng)地下水環(huán)境脆弱性分區(qū)一覽表

2）脆弱性評(píng)價(jià)

按上述分區(qū)結(jié)果及評(píng)價(jià)方法，對(duì)各個(gè)分區(qū)評(píng)價(jià)因子賦分值，再按照上述模型求取DI值，最后按表6進(jìn)行評(píng)價(jià)，具體評(píng)價(jià)結(jié)果見(jiàn)表8、圖2。

圖2 光生壩地下河系統(tǒng)固有脆弱性評(píng)價(jià)分區(qū)圖

1.脆弱性極高(極易污染)；2.脆弱性很高(很易污染)；3.脆弱性較高(較易污染)；4.脆弱性很低(很難污染)；5.脆弱性分區(qū)界線及代碼；6.地下河管道及出口；7.受污染的地下河管道及出口；8.下降泉；9.表層巖溶泉；10.地下水開(kāi)采井；11.伏流入口；12.落水洞；13.污水處理廠

3.6 系統(tǒng)內(nèi)特殊脆弱性評(píng)價(jià)

1）系統(tǒng)內(nèi)分區(qū)

根據(jù)系統(tǒng)內(nèi)的巖性特征、水文地質(zhì)條件等，將光生壩地下河系統(tǒng)分為3個(gè)亞區(qū)，在根據(jù)污染源分布情況進(jìn)一步分亞區(qū)，分區(qū)依據(jù)見(jiàn)表9。

2）脆弱性評(píng)價(jià)

按上述分區(qū)結(jié)果及評(píng)價(jià)方法，對(duì)各個(gè)分區(qū)評(píng)價(jià)因子賦分值，再按照上述模型求取DI值，最后按表7進(jìn)行評(píng)價(jià)，具體評(píng)價(jià)結(jié)果見(jiàn)表10、圖3。

4 結(jié)論

根據(jù)開(kāi)陽(yáng)縣光生壩地下河系統(tǒng)情況，選擇10個(gè)指標(biāo)因子，采用層次分析法對(duì)地下水脆弱性進(jìn)行評(píng)價(jià)，地下水脆弱性分布圖較好反映了研究區(qū)實(shí)際情況，評(píng)價(jià)脆弱性為很高或極高的地區(qū)是污染相對(duì)嚴(yán)重地區(qū)。具有一定的相關(guān)性。

結(jié)果中，光生壩地下河系統(tǒng)劃分為3個(gè)亞區(qū)：石灰?guī)r分布區(qū)（A1）、白云巖分布區(qū)（A2）、碎屑巖分布區(qū)（A3），僅與巖性相關(guān)，未反映“地下河系統(tǒng)”徑流特征；未結(jié)合隨地下河流動(dòng)污染物彌散擴(kuò)散機(jī)理進(jìn)行評(píng)價(jià)。是今后值得深入研究方面。

圖3 光生壩地下河系統(tǒng)特殊脆弱性評(píng)價(jià)分區(qū)圖

1.脆弱性極高(極易污染) 2.脆弱性很高(很易污染) 3.脆弱性較高(較易污染) 4.脆弱性很低(很難污染) 5.脆弱性分區(qū)界線及代碼 6.地下河管道及出口 7.受污染的地下河管道及出口 8.下降泉 9.表層巖溶泉 10.地下水開(kāi)采井 11.伏流入口 12.落水洞 13.污水處理廠

表8 光生壩地下河系統(tǒng)地下水環(huán)境脆弱性評(píng)價(jià)分區(qū)一覽表

表9 光生壩地下河系統(tǒng)特殊脆弱性分區(qū)一覽表

表10 光生壩地下河系統(tǒng)特殊脆弱性評(píng)價(jià)分區(qū)一覽表

王若帆，趙翠，焦恒，江峰．2020．貴州省1：5萬(wàn)水文地質(zhì)調(diào)查成果報(bào)告養(yǎng)龍司幅（G48E005020）、開(kāi)陽(yáng)縣幅（G48E006020）[R]．

吉勤克補(bǔ)子，李強(qiáng)．2020．貴州息烽溫泉水文地球化學(xué)特征及地質(zhì)成因研究[J]．四川地質(zhì)學(xué)報(bào)，40(3)：434-438.

孫才志，潘?。?999．地下水脆弱性的概念、評(píng)價(jià)方法與研究前景[J]．水科學(xué)進(jìn)展， 10(4)： 444-449.

孫才志，林山杉．2001．地下水脆弱性概念的發(fā)展過(guò)程與評(píng)價(jià)現(xiàn)狀及研究前景[J] ．吉林地質(zhì)，19(1)：30-36.

雷靜．2002．地下水環(huán)境脆弱性的研究[D]．北京：清華大學(xué)．

范琦，王貴玲，藺文靜，陳浩．2007．地下水脆弱性評(píng)價(jià)方法的探討及實(shí)例[J]．水利學(xué)報(bào)， 38(5)：601-605.

郭曉靜，周金龍，靳孟貴，卜新峰． 2010．地下水脆弱性研究綜述[J]．地下水， 32(3)：1-5.

江峰，劉漢武， 吉勤克補(bǔ)子，王若帆，焦恒．2021．單因子水質(zhì)標(biāo)識(shí)指數(shù)法在貴州省洋水河流域地下水水質(zhì)評(píng)價(jià)中的應(yīng)用[J]．四川地質(zhì)學(xué)報(bào)， 41(1)：151-153+176.

Groundwater Vulnerability Assessment of the Guangshengba Underground River System in Kaiyang, Guizhou

JIAO Heng WANG Ruo-fan JIANG Feng YI Rui

(The 114th Geological Party, Guizhou Bureau of Geology and Exploration and Exploitation of Mineral Resources, Zunyi, Guizhou 563000)

This paper evaluates comprehensively the inherent vulnerability and special vulnerability of groundwater of the Guangshengba underground river system in Kaiyang, Guizhou by means of analytic hierarchy process and evaluation factors such as lithostratigraphy (R), aqueous media type (W), karst type (K), degree of karstification development (D),topography and geomorphology (T), soil thickness (S) vegetation coverage (V),degree of engineering activity (E) , pollution source (P) and population density (C).

vulnerability of groundwater; analytic hierarchy process; evaluation factor; Kaiyang, Guizhou

2020-08-20

貴州省地勘基金項(xiàng)目“貴州省1：5萬(wàn)水文地質(zhì)調(diào)查”（2019-1號(hào)）,貴州省地質(zhì)礦產(chǎn)勘查開(kāi)發(fā)局地質(zhì)科研項(xiàng)目“黔中“三磷”分布區(qū)與地下水有關(guān)的地質(zhì)環(huán)境問(wèn)題精準(zhǔn)識(shí)別研究—以開(kāi)陽(yáng)縣洋水背斜分散排泄系統(tǒng)為例”黔地礦科合[2020]23號(hào)資助。

焦恒（1992—），男，湖北荊州人，工程師，從事水工環(huán)地質(zhì)研究工作

P641.8

A

1006-0995（2021）02-0283-06

10.3969/j.issn.1006-0995.2021.02.020