李中原， 王國重， 張繼宇， 程煥玲， 趙蓓蓓， 李永麗

(1.河南省水文水資源局， 河南 鄭州 450003； 2.黃河水文水資源科學研究院， 河南 鄭州 450004； 3.黃河養(yǎng)護集團有限公司， 河南 鄭州 450003； 4.河南省水土保持監(jiān)督監(jiān)測總站，河南 鄭州 450008； 5.河南省水利勘測設計研究有限公司， 河南 鄭州 450008)

1 研究背景

彰武水庫位于洹河大峽谷出口處，控制流域面積120 km2，總庫容7 063×104m3，興利庫容1 703×104m3，其主要作用是保障安鋼、電廠用水及農(nóng)業(yè)引蓄灌溉，兼營發(fā)電、養(yǎng)殖。水庫設計灌溉面積3.3×104hm2，可養(yǎng)殖水面200 hm2[1]。

作為重要的工農(nóng)業(yè)用水和地下水補給源，該水庫水環(huán)境的安全與否，對安陽市經(jīng)濟社會的可持續(xù)發(fā)展意義重大[2-3]。2011年苗利芳等[4]首次對“十一五”以來安陽市地表水水質(zhì)進行評估，認為彰武水庫處于輕度富營養(yǎng)化狀態(tài)。2017年李世華等[5]進一步指出安陽市“十二五”期間地表水系水質(zhì)總體污染程度為重度污染。2017年唐敏等[6]根據(jù)安陽市環(huán)境保護監(jiān)測中心2006-2015年的實測數(shù)據(jù)，對彰武水庫的營養(yǎng)化狀況進行了深入研究，認為十年間，水庫營養(yǎng)狀態(tài)已由中營養(yǎng)發(fā)展為輕度富營養(yǎng)，還有進一步惡化的趨勢。

水庫水環(huán)境安全與否直接關系著流域的供水和飲用水安全[7]，水庫水環(huán)境安全評估是保障飲水安全和水庫可持續(xù)利用、實現(xiàn)人與自然和諧的基礎，對于水環(huán)境安全進行評估的基本思路是：首先構(gòu)建評估指標體系，然后選用評估方法、確定各指標的權(quán)重進行評估，最后分析評估結(jié)果的合理性[8]。構(gòu)建評估指標體系國際上常用的模型主要有國際經(jīng)濟合作與發(fā)展組織提出的“壓力-狀態(tài)-響應(PSR)”模型、聯(lián)合國可持續(xù)發(fā)展委員會推薦的“驅(qū)動力-狀態(tài)-響應(DSR)”模型、歐洲環(huán)境署建議的“驅(qū)動力-壓力-狀態(tài)-影響-響應(DPSIR)”模式[9]。其中，PSR模式所構(gòu)建的指標體系結(jié)構(gòu)簡單、符合邏輯，常常被人們選用[10]；DSR模型中驅(qū)動力指標與狀態(tài)指標之間的邏輯關系模糊，應用相對較少[11]；DPSIR模式是對PSR體系的改進，添加了“驅(qū)動力”和“影響”兩類指標，使指標體系在邏輯上更為完善[12]。常用的評估方法主要有層次分析法、主成分分析法、綜合指數(shù)法等[13]。王國重等[14]采用PSR模型構(gòu)建了河南省大型水庫生態(tài)安全評估指標體系，以熵值法計算了各指標的權(quán)重，通過綜合指數(shù)法對水庫2005、2015年的生態(tài)安全狀況進行了評估；喻立等[15]通過DPSIR模型建立了沙湖濕地健康評價指標體系，根據(jù)層次分析法確定指標權(quán)重，并用模糊綜合評價方法對該濕地進行了評估；李永麗[16]基于DPSIR模型創(chuàng)建了陸渾水庫水環(huán)境安全評估指標體系，根據(jù)綜合營養(yǎng)狀態(tài)指數(shù)認為該水庫長期處于中營養(yǎng)狀態(tài)。

上述評估方法都存在著不同程度的缺陷：層次分析法確定權(quán)重的方式過于主觀，會造成極大的偏差；主成分分析法能夠計算指標的權(quán)重，但會導致信息的缺失；綜合評估方法難以克服評估過程中的隨機性和模糊性問題[17]。上述評估方法的共性就是沒有克服評估過程中主客觀因素的影響，如對同一個評估指標而言，不同方法設置界限的方式或者所設定的界限不盡相同[18]。云模型是1995年李德毅教授在研究概率論和模糊數(shù)學的基礎上提出來的，能夠建立起定性概念與定量表示之間的轉(zhuǎn)換關系，可有效處理評估過程中的不確定性[19-20]。本文以彰武水庫為典型，以DPSIR模型來建立水庫環(huán)境安全的評估指標體系，嘗試采用云模型對該水庫的環(huán)境安全狀況進行評估。

2 基于DPSIR模型的彰武水庫環(huán)境安全評估指標體系

2.1 DPSIR指標體系

借鑒環(huán)境安全和生態(tài)安全指標體系[21-22]，根據(jù)水庫水環(huán)境安全概念和水庫功能，確定了彰武水庫環(huán)境安全評估指標體系：驅(qū)動力(降雨量)、壓力(污染物入庫量)、狀態(tài)(污染物濃度)、影響(水質(zhì)狀況)、響應(透明度、溶解氧)，共11項指標，分屬5個方案層，各層間的關系為：污染物(TN、TP、高錳酸鹽指數(shù))隨降雨產(chǎn)生的地表徑流進入水庫，增加了入庫污染物數(shù)量和濃度，影響水體水質(zhì)，使得指標透明度和溶解氧發(fā)生相應的變化。水庫的指標體系詳見圖1。

圖1 彰武水庫水環(huán)境安全評估指標體系

根據(jù)彰武水庫2015年水質(zhì)監(jiān)測數(shù)據(jù)，對其水環(huán)境安全狀況進行評估。其中雨量數(shù)據(jù)來自庫區(qū)雨量站的監(jiān)測值；污染物的入庫量則通過入庫河流的流量、入庫排污口監(jiān)測數(shù)據(jù)和水庫水質(zhì)監(jiān)測數(shù)據(jù)計算得到；水質(zhì)狀態(tài)指標TN、TP、Chl-a、CODMn和響應指標透明度、溶解氧的數(shù)值由河南省環(huán)境監(jiān)測中心提供；依據(jù)《地表水環(huán)境質(zhì)量標準( GB3838-2002)》中規(guī)定的II類水水質(zhì)目標，判斷各月水質(zhì)是否達標。2015年各月水質(zhì)指標的監(jiān)測值詳見表1。

表1 2015年彰武水庫水質(zhì)監(jiān)測數(shù)據(jù)

2.2 指標權(quán)重的確定

設n個評估指標P1,P2,…,Pn，其中的m個樣本為S1,S2,…,Sm，Qij為樣本Si在指標Pj時的實測值，這樣就產(chǎn)生了m個樣本n個指標的矩陣(Qij)m×n。該矩陣中各指標的數(shù)量級、量綱等差異很大，需要對該矩陣中的數(shù)據(jù)進行規(guī)范化處理，處理公式見公式(1)、(2)所示[23]。

(1)

(2)

式中：Qij為第i個樣本中第j個指標的實測值；yij為經(jīng)過標準化處理后的第i個樣本中第j個指標的值； maxQij和minQij分別為實測值Qij的最大值和最小值。

熵值法因為受主觀因素影響較小，常被研究者用于確定指標的權(quán)重[24-26]。這里也采用熵值法，第j項指標的熵值Ej可由斯蒂林公式得到，詳見公式(3)。而某項指標的效用價值則取決于該指標的熵值Ej與1的差值。

(3)

式中：常數(shù)k=(lnm)-1。

hj=1-Ej

(4)

則第j項指標的權(quán)重為：

(5)

各指標經(jīng)規(guī)范化處理后，再由公式(3)～(5)計算得到各個指標的權(quán)重，見表2所示。

3 云模型水庫水環(huán)境安全綜合評估

3.1 云模型的生成

云模型可構(gòu)建水庫環(huán)境安全評估中定性概念與定量之間的映射。設U為用精確數(shù)值表示的定量論域，C為U上水庫水環(huán)境安全的一個定性概念，當x∈U且x為定性概念C的一次隨機實現(xiàn)，則x在U上的分布顯示為云，每個x稱為一個云滴，多個云滴形成一朵云，可表征定性概念與數(shù)域空間之間的映射關系[27]。云模型以數(shù)字特征值Ex、En、He和云滴數(shù)目N作為輸入，以N個云滴數(shù)的定量值及每個云滴的隸屬度作為輸出。生成水庫水環(huán)境安全評估的云模型具體步驟如下[28]：

(3)設云滴x為定性概念的一次具體量化值，計算其確定度μ；

(6)

(4)重復步驟(1)～(3)，直到產(chǎn)生N個云滴。

3.2 水庫水環(huán)境安全的云模型評估

采用云模型建立水庫水環(huán)境安全評估模型的步驟如下：

(1)按照《地表水環(huán)境質(zhì)量標準( GB3838-2002)》，結(jié)合水庫功能，確定了各水質(zhì)指標的分級標準，如表3所示。

表3 水庫水質(zhì)評價指標分級標準

(2)按照各水質(zhì)指標不同安全等級所在的閾值范圍，確定各級別云模型相應的各數(shù)字特征Ex、En、He。

云模型數(shù)字特征的計算公式如下[29]：

(7)

式中：umax、umin分別表示變量的最大和最小邊界；k為常數(shù)，令k=0.1。若評估指標的評估等級區(qū)間為：Ⅰ(0，a)、Ⅱ(a，b)、Ⅲ(b，c)、Ⅳ(c，+∞)，則各等級云模型數(shù)字特征的計算如表4所示，參數(shù)a、b、c的取值是根據(jù)表3中設置的各評估指標的分級標準而定，具體詳見表5。

表4 云模型數(shù)字特征的計算規(guī)則

根據(jù)表4，當評估因子的取值范圍為有界區(qū)域時，可通過正態(tài)云發(fā)生器生成正態(tài)云模型；當評估因子的取值無限制范圍時，可通過正態(tài)云、半升云、半降云發(fā)生器形成組合云。

(3)計算各指標不同評估等級的確定度，并乘以相應的指標權(quán)重，再將各指標同等級別的確定度與其權(quán)重的乘積相加，所在級別中的最大值就是水庫水環(huán)境對應的安全等級。

4 評估結(jié)果與分析

根據(jù)確定的評估指標和安全等級，分別計算了相應的3個數(shù)字特征值，生成各指標在不同等級下的云模型。

結(jié)合2015年彰武水庫水質(zhì)實測平均數(shù)據(jù)，得到了不同等級下各指標的確定度，再根據(jù)指標權(quán)重，加權(quán)求和得到了不同安全等級的綜合確定度，其最大值即為水庫2015年所處的安全等級。

彰武水庫各評估指標相應的云模型評估等級下的確定度如表6所示， 2015年不同安全等級的綜合確定度詳見表7。

表5 各評估指標中參數(shù)a、b、c的取值

表6 各評估指標相應等級下的確定度

表7 不同安全等級的綜合確定度

由表6、7可知：2015年彰武水庫水環(huán)境安全評估等級為較安全，但有3個指標為不安全級別，分別為TN濃度、TP濃度、水質(zhì)達標率，降雨量指標處于較不安全級別,這正是其水質(zhì)為劣V類的原因所在。彰武水庫周邊廢棄的煤礦礦山較多，屬于計劃經(jīng)濟時代的歷史遺留問題，這些工業(yè)廢水、廢渣沿著溝道排放、堆積，再加上農(nóng)業(yè)化肥和農(nóng)藥在降雨和灌溉退水作用下排入水體，使水中N、P含量增加[30]；水產(chǎn)養(yǎng)殖尤其是大量的網(wǎng)箱養(yǎng)殖[31]，加大了水體中自源性污染，養(yǎng)殖飼料、魚類排泄物及人為無序開發(fā)，增加了水體中TN含量，使藻類等浮游生物大量繁殖，溶解氧含量降低。上述原因?qū)е抡梦渌畮焖|(zhì)難以達標。

5 結(jié)論與建議

通過DPSIR模型建立了彰武水庫水環(huán)境安全評估的指標體系，采用熵值法確定了各個指標的權(quán)重，應用云模型對其2015年的水環(huán)境安全進行了評估，主要結(jié)論如下：

(1)該水庫水環(huán)境處于較安全等級，但存在安全隱患；由于廢棄的工礦企業(yè)較多，其工業(yè)廢水和周邊農(nóng)田中化肥農(nóng)藥在降雨徑流和灌溉退水作用下進入水體，再加上傳統(tǒng)的水產(chǎn)養(yǎng)殖模式，造成水庫N、P濃度增加，水質(zhì)全年為劣V類，水庫管理部門應引起重視。

(2)云模型能夠揭示事物的隨機性與模糊性，對評估指標進行分級，利用熵值法確定各指標的權(quán)重，根據(jù)各指標不同等級的綜合確定度，實現(xiàn)了定性概念與定量數(shù)值表達之間的轉(zhuǎn)換，在一定程度上克服了評估過程中的不確定性。

針對水庫周邊廢棄的工礦企業(yè)，一方面需對其治理整頓、拍賣或者租售，使其能發(fā)揮效益，對其廢污水實施總量控制、達標排放，另一方面采取水土保持措施治理水土流失和面源污染，同時倡導生態(tài)健康養(yǎng)殖，取締網(wǎng)箱養(yǎng)殖，確保水庫水質(zhì)安全。