劉慧艷

(凌源市水利建筑勘測設計院，遼寧 凌源 122500)

0 引 言

青龍河流域降雨量較大，支流縱橫交錯，河道坡度大，洪水總量多，又無大流域調(diào)蓄，形成峰高量大的特點，其極具破壞性。大雨成大災，小雨成小災，歷年成災，特別是近些年，氣候急劇變化，局部降雨頗多，成災的機率增加，給當?shù)剞r(nóng)業(yè)的高效發(fā)展以及沿河居民財產(chǎn)安全帶來極大威脅。為合理開發(fā)青龍河水資源，減輕其支流星干河、五道梁河、都源河、湯道河、大北河以及下游的防洪負擔，解決超標準洪水威脅和流域洪澇災害的問題，先后在中上游修建了中小型水庫、分洪閘、攔河閘等水利樞紐，立項實施了筑堤加固和疏浚治理工程。

在城市供水、農(nóng)業(yè)灌溉、流域防洪等方面水庫閘壩工程起著重要作用，同時也帶來了許多不利影響。近年來，隨著上游用水量增加河流接受的生產(chǎn)生活廢水量不斷增多，青龍河流域逐漸呈現(xiàn)生態(tài)環(huán)境惡化、斷流、干旱等問題。河流閘壩控制工程大幅度削減了青龍河主河道水體自凈能力，其水質(zhì)標準已超過Ⅴ類地表水。河流的生態(tài)系統(tǒng)結(jié)構(gòu)受嚴重污染的河水已出現(xiàn)破壞，青龍河的生態(tài)功能不斷退化，嚴重威脅著地下水質(zhì)及當?shù)鼐用竦娘嬎踩?，對區(qū)域經(jīng)濟可持續(xù)發(fā)展產(chǎn)生明顯的制約。依據(jù)區(qū)域規(guī)劃，青龍河主要承擔著下游農(nóng)業(yè)灌溉、排除農(nóng)田與城區(qū)澇水、中下游景觀娛樂、安全下泄上游洪水等功能，故定量評價青龍河生態(tài)影響非常重要[1-5]。

國內(nèi)外研究評價河流生態(tài)影響的成果豐碩，如Karr、White等給出了河流生態(tài)健康狀況評價的生物完整性指數(shù)法和表征河流水文特征的參數(shù)指標；Brierley等提出了河流形態(tài)結(jié)構(gòu)的框架及其評價控制流程。對于河流健康評價我國相關(guān)研究相對較晚，許多學者對國內(nèi)河流生態(tài)健康借鑒國外理論與經(jīng)驗開展創(chuàng)新性探究，形成了豐碩的研究理論和方法，如均方差法、專家調(diào)查法、層次分析法、主成分分析和熵值法等，并取得理想的實踐應用效果。宋剛福等對河流生態(tài)健康利用改進的方法實行全面、逐層的評價[6]。河流生態(tài)評價覆蓋范圍廣、涉及因素多，因素集通常存在屬性復雜性、模糊性、隨機性和多層次等特點，不同參評因素間具有復雜關(guān)系，甚至可能出現(xiàn)不相容的單因素評判結(jié)果。鑒于此，本研究對河流生態(tài)影響評價中引入定量與定性相互轉(zhuǎn)換的云理論，有機地結(jié)合評價等級的隨機性與模糊性，從而構(gòu)建生態(tài)評價云模型，以實現(xiàn)不確定性評價河流生態(tài)狀況。

1 云模型

1.1 云概念

定義1：設U、C是用精確數(shù)值表征的定量論域以及論域U上的定性概念，若C的一次隨機實現(xiàn)為定量值x且x∈U，則有穩(wěn)定傾向的隨機數(shù)μ(x)∈[0,1]就是因素x對C的隸屬函數(shù)，其數(shù)學形式為：

μ(x)：U→[0，1]，?x∈U，x→μ(x)

(1)

論域U上x的分布以及每一個x稱為云和云滴，可見隸屬函數(shù)與任意因素x均存在一對多的轉(zhuǎn)化關(guān)系，兩者并非傳統(tǒng)的模糊隸屬關(guān)系。對此，云模型利用三個數(shù)字特征期望Ex、熵En、超熵He表征定量特性和定性概念，并記為(Ex，En，He)。

1.2 正態(tài)云模型

定義2：依據(jù)以上定義，若存在條件x-N(Ex,En'2)，且En'-N(En,He2)，則利用下式計算C的確定度：

(2)

則論域U上x的分布稱為正態(tài)云，正態(tài)云能夠客觀地反映各指標的特征以及各類數(shù)據(jù)分布情況，所以其普適性較強。

定性概念的產(chǎn)生利用正態(tài)云發(fā)生器實現(xiàn)，其相關(guān)運算有：①以(Eni'2)與He2為方差、以Ex與En為期望值，生成隨機數(shù)xi=NORM(Ex，Eni'2)和Eni'2=NORM(En，He2)，其中NORM( )為正態(tài)分布函數(shù)；②然后借助Matlab軟件和公式μi=e-[(xi-Ex)2/2(Eni')2]生成正態(tài)分布隨機數(shù)。

通過運算可以在各評價等級邊界內(nèi)生成一個云滴(xi，μi)，重復以上步驟直至達到預先設定的n次要求，從而生成擁有n個云滴數(shù)的云，并實現(xiàn)定性概念的表征。

1.3 云評價模型

步驟1：設U={u1,u2,…,un}、V={v1,v2,…,vm}、W={w1,w2,…,wn}為生態(tài)影響評價的有限因素集、評語集和權(quán)重集，m、n為等級標準數(shù)與參評因子數(shù)。定性概念的數(shù)字特征利用云模型及確定的單因素評價標準反映，記為(Ex，En，He)。設U={u1，u2，…，un}、V={v1,v2,…,vm}、W={w1,w2,…,wn}為生態(tài)影響參評因素i對應等級j的上、下邊界值，則可以用云模型表征等級j下因素i的定性概念，即：

Exi,j=(xi,j1+xi,j2)/2

(3)

兩個相鄰等級的界限值通常為邊界值，該條件下隸屬于兩種級別，對此可視為隸屬于兩個相鄰級別的程度相等，即：

(4)

超熵Hei,j代表云滴的凝聚程度，其數(shù)值越大則正態(tài)云越厚，熵的不確定性越高，結(jié)合試驗或?qū)＜医?jīng)驗綜合確定Eni,j的大小。

步驟2：根據(jù)每個等級下各參評因子的云數(shù)字特征(Exi,j,Eni,j,Hei,j)，采用所有待評樣本的特征值和正向云發(fā)生器生成相應的隸屬度矩陣，其表達式為：

(5)

隨機隸屬矩陣R反映了評語集V與因素集U之間的關(guān)系，因素集U所對應的權(quán)重模糊子集為W，并利用模糊數(shù)學法轉(zhuǎn)換處理隸屬度矩陣R與權(quán)重集W生成模糊子集B，其表達式為：

B=W?R

(6)

式中：⊕為論域V的一種運算方式；B=(b1，b2，…，bm)，代表參評樣本隸屬于某等級的程度。

步驟3：依據(jù)模糊子集B與最大隸屬準則決策過程中，考慮到隸屬度準則在隸屬等級相等情況下可能發(fā)生失效的問題，為保證評價結(jié)果的科學性與準確性，將評價等級利用模糊可辨識別模型確定，其數(shù)學公式為：

(7)

(8)

2 青龍河生態(tài)影響評價

2.1 構(gòu)建評價指標體系

文章結(jié)合青龍河流域?qū)嶋H情況以及河流生態(tài)健康相關(guān)研究成果，遵循動態(tài)性與連續(xù)性、靈活性與可比性、可操作性與可行性、層次性與系統(tǒng)性等原則，綜合利用專家咨詢法、理論分析法和頻度統(tǒng)計法，從反映河流生態(tài)變化趨勢、生態(tài)影響現(xiàn)狀、生物狀況、棲息地質(zhì)量、水環(huán)境情況、水文條件等角度合理選擇參評因子，并構(gòu)建梯階層次評價體系，如表1所示。

表1 青龍河生態(tài)影響評價體系及其現(xiàn)狀值

采用百分制的專家調(diào)研打分法量化評價體系中的定性指標，如水生生物存活狀況、水域濕地狀況、垂向透水性、橫向連通性、河岸與河床穩(wěn)定性指標；然后利用下述量化公式計算其它定量指標：

2.2 生態(tài)影響評價計算

將青龍河流域按照防洪規(guī)劃和歷史資料劃分為上游、中游、下游段，并利用專家賦分法和定量因子計算公式，確定各河段生態(tài)健康參評因子值，如表1。設定流域處于基本自然狀態(tài)的基點為流域開發(fā)程度低于5%，并以分值閾0.8-1.0、0.6-0.8、0.4-0.6、0.2-0.4、0.0-0.2代表健康、基本健康、亞健康、病態(tài)、瀕于崩潰5個級別，所對應的生態(tài)影響程度依次為很小(Ⅰ級)、一般(Ⅱ級)、較嚴重(Ⅲ級)、嚴重(Ⅳ級)、非常嚴重(Ⅴ級)，各等級的生態(tài)學意義如表2所示。

表2 生態(tài)影響評價等級

根據(jù)國家適用標準、現(xiàn)有研究成果和青龍河相關(guān)歷史資料，通過邀請領(lǐng)域內(nèi)的專家參與評判以及綜合對比分析多個區(qū)域，經(jīng)適當調(diào)整明確不同等級的評分標準，并對所有參評因子權(quán)重按照文獻中的層次分析法及其流程確定，如表3所示。

表3 不同等級的指標標準值

2.3 結(jié)果與分析

根據(jù)各參評指標的現(xiàn)狀值、權(quán)重值和標準特征值，對青龍河流域上、中、下游段的生態(tài)影響狀況利用云模型評價，如表4所示。為進一步驗證云模型的可行性與合理性，對比分析了改進拉開檔次法、拉開檔次法和層次分析法的評價結(jié)果。

表4 青龍河生態(tài)影響評價值

對青龍河流域上、中、下游的生態(tài)影響狀況利用云模型進行評價，其綜合指數(shù)為0.546、0.531、0.560，依據(jù)等級劃分標準青龍河總體處于亞健康狀態(tài)，人類活動對河流系統(tǒng)已產(chǎn)生顯著干擾，水生生物的繁衍和生存受到較大影響。其中，生態(tài)狀況最差段為中游，而下游段較好，其原因為中上游用水過度導致下游河段經(jīng)常出現(xiàn)斷流現(xiàn)象，下游徑流量年內(nèi)存在較大變化。由各子系統(tǒng)權(quán)值可知，河流生態(tài)系統(tǒng)受水環(huán)境狀況的影響最大，青龍河整體上達到中污染程度，部分河段處于嚴重污染現(xiàn)狀。2015-2018年實測資料顯示，青龍河水質(zhì)狀況除汛期外均達到劣Ⅴ類；此外，從河岸帶狀況和河流形態(tài)結(jié)構(gòu)上，青龍河流域也面臨著嚴峻形勢，由于不合理的開發(fā)和濱河濕地侵占等，現(xiàn)階段亟待加強生態(tài)治理及環(huán)境保護。評價結(jié)果能夠客觀、真實的反映青龍河實際情況，可見對于河流生態(tài)影響評價云模型具有較強的適用性、可靠性。

3 結(jié) 論

文章以青龍河流域為例，對其上、中、下游的生態(tài)影響狀況利用能夠?qū)崿F(xiàn)定量與定性轉(zhuǎn)化的云模型進行評價。研究表明，青龍河流域總體達到亞健康狀態(tài)，其中生態(tài)狀況最差段為中游，而下游段較好，為恢復自然水文情勢有必要實施閘壩生態(tài)調(diào)度，通過非工程與工程措施實現(xiàn)生態(tài)補水和修復。

云模型具有考慮評判的隨機性與參評因素概念模糊性的特點，通過對比改進拉開檔次法、拉開檔次法和層次分析法評價結(jié)果，進一步驗證了云模型的實用性、可靠性，對于多因素復雜問題的評價該模型具有較高的適用性。