基于水資源配置網絡圖的水系連通評價體系研究
——以廊坊市為例

梁彩欣，張守平

(重慶交通大學，國家內河航道整治工程技術研究中心，重慶 400074)

1 研究背景

國外早在20世紀中便開展了一系列關于水系連通方面的研究，在提出河湖水系連接概念的基礎上，同時也對其指標、體系、方法和應用進行了研究探討[1-6]。我國在水系連通研究方面起步較晚，但發展迅速。2005年首次將水系連通性作為一個重要指標列入健康長江指標體系中[7]，在2010年陳雷提出“河湖連通是提高水資源配置能力的重要途徑”[8]后，關于河湖水系連通的相關研究快速發展，先后經歷了定性分析與定量分析兩個階段。學者們[9-12]首先對于“河湖水系連通是什么”進行了一系列研究，通過探討河湖水系連通歷史演變過程，分析水系連通概念框架、功能機理與驅動因素，構建水系連通評價體系。除定性分析之外，學者們[13-17]基于不同地區需求，改進圖論與水文模擬方法建立HEC_HMS水文模型、思索水系盒維數與連通度相關性、探討平原河網區水系連通對改善水環境的效果評估、構建水系連通形態指標與水系連通功能指標定量關系式。但是以上研究多聚焦于衛星地圖上真實存在的河流湖泊水系，缺乏為滿足人類生產生活用水和社會持續穩定發展需求所建設水利工程對水系連通變化的具體研究。基于以上認識，本文從二元水循環水資源配置網絡圖出發，結合生產生活用水水平評價指標以及自然水系連通形態指標，構建一套用于描述在不同時間節點下的調水工程對自然—人工二元結合水系連通度影響的綜合評價指標體系，并以廊坊市為例開展應用研究，力圖為在人類活動影響下的水系連通度的預測評估提供依據。

2 綜合評價方法

現階段對于天然水系連通度的指標選取已初具雛形，往往通過結構形態指標和連通形態指標進行評價，但在指標選取中忽略人類活動所產生的取水、排水行為對于水系連通的影響。基于二元水循環水資源配置網絡圖中對計算單元水傳輸系統(線)的定義，結合生產生活用水水平評價指標以及自然水系連通形態指標，運用主成分分析法對所選取的指標進行歸一化處理，以確定各個指標的貢獻率與權重系數，然后，通過集對分析法對評價指標的同一度、對立度、差異度進行分析，得出不同時段水系連通度所對應等級。

2.1 網絡圖原理與概化

水資源優化配置模型采用規劃與優化相結合的方法，基于線性規劃構建供需平衡分析系統和耗水平衡分析系統，以水量平衡原理為基礎，對網絡圖中各個控制節點、水庫、計算單元建立平衡方程和約束方程，在目標函數最大(或最小)情況下，通過反復迭代求解決策向量的最優值，以此作為社會經濟系統和生態環境系統對水資源需求的參考結果[18]。系統網絡圖是構建水資源優化配置模型的基礎，主要由節點(點)、計算單元水傳輸系統(線)、水資源分區水傳輸系統(面)3類要素構成：

(1)節點(點)可分為需水節點、供水節點和輸水節點3類，如圖1中區縣鎮級計算單元、湖泊濕地、入海匯點、水閘、引水節點、匯水節點均為網絡圖中點元素。

(2)計算單元水傳輸系統(線)分為天然循環與人工側支循環，其中河流系統屬于天然水循環系統，而人工側支水循環系統包括調水系統、排水系統、供水系統。在任意兩個節點間連接若干條有向連線表示的水傳輸關系，并構成的網狀圖形。

(3)水資源分區水傳輸系統(面)由一個或多個計算單元組成，且流域內所有的面按照水傳輸關系連接起來，構成了水資源分區之間的水傳輸系統。

對元素進行抽象簡化后，依據水量平衡原理，以及各節點按照其上、下游連線的實際對應關系和邏輯關系，對水資源優化配置系統網絡圖進行繪制(圖1)，以節點、水傳輸系統構成的網狀圖形反映廊坊市水資源系統、經濟系統和生態環境系統間內在的邏輯關系，而節點與節點之間上下游連線則展現出廊坊市內天然和人工側支兩個循環系統的實際對應關系，為建立評價自然—人工二元結合水系連通度的指標體系提供研究前提。

2.2 KMO檢驗與主成分分析

KMO(Kaiser-Meyer-Olkin)檢驗統計量是用于比較變量間簡單相關系數和偏相關系數的指標，主要應用于多元統計的因子分析。KMO統計量是取值在0和1之間，當所有變量間的簡單相關系數平方和遠遠大于偏相關系數平方和時，KMO值越接近于1，意味著變量間的相關性越強，原有變量越適合作因子分析；當所有變量間的簡單相關系數平方和接近0時，KMO值越接近于0，意味著變量間的相關性越弱，原有變量越不適合做因子分析。

主成分分析法通過降維處理可消除評價指標間的相關影響，減少指標選擇的工作量。在綜合評價函數中，各主成分的權數為其貢獻率，它反映了該主成分包含原始數據的信息量占全部信息量的比重，這樣確定權數是客觀的、合理的，克服了主觀確定權數的不足。

2.3 集對分析法

集對分析(set pair analysis，SPA)是趙克勤在1989年提出的一種新的系統分析方法[19]，它研究兩個事物對共同屬性的同異反程度，即聯系度。集對分析方法解決了不同評價方法對同一個評價對象給出的不同排序結果的協調性問題，其研究結果更加貼近實際，能更好地解決多屬性、不確定系統的評價問題[20]。

設有N個評價指標，其中有S個評價指標優于Ⅰ級標準，有P個指標劣于Ⅲ級標準，有F個指標界于Ⅰ和Ⅲ之間，則可得到聯系度μ的表達式：

(1)

其中，令a=S/N，b=F/N，c=P/N，則a、b、c依次稱為同一度、差異度、對立度。且可簡寫為：

μ=a+bi+cj

(2)

具體計算步驟如下：

(1)對各年份連通度綜合評價值按照指標等級進行篩選，得出各個年份聯系度μ；

(2)判斷指標值類型為成本型還是效益性，成本型指標是指實測值越小評價等級越高的指標，效益型是指實測值越大評價等級越高，計算各個評價指標對各個年份的聯系度及聯系度均值。

成本型指標：

(3)

效益型指標：

(4)

式(3)、(4)中:SⅠ、SⅡ、SⅢ分別為評價標準的Ⅰ、Ⅱ、Ⅲ級的門限值；μ為各個年份對應的水系連通度的指標值；m為第m個年份；n為第n個評價指標。

(3)計算聯系主值數，取i=0，j=-1，并代入步驟(2)中計算的各個年份所對應的a、b、c，得到聯系主值數。

(4)將[-1,1]區間均分為3等份，從右至左分別對應I級、Ⅱ級、Ⅲ級，將得到的聯系主值數與評價等級進行對比，聯系主值數越大則水系的連通度越好，得到該年份水系連通度對應等級。

3 應用實例

廊坊市地處河北平原中北部、海河流域中下游，地形地貌98%為平原。全市流域面積50 km2以上河流110條，長度2 021 km；流域面積50 km2以下渠系約4 252 條，長度5 096 km。全市河網密度為0.32 km/km2，屬于海河流域水系較發達區域。廊坊市是河北省唯一沒有地表水庫的設區市，境內全部為過境河流，受上游城市水量調控影響巨大，絕大多數河道常年無水，水資源分配不平衡問題嚴重。2017年，人均水資源占有量為114 m3，僅為河北省人均水資源量的61.8%、全國人均水資源量的5.5%，缺水形勢嚴峻。同年廊坊市整體水資源開發利用率為51.5%，其中，地表水資源開發利用率為14.6%，地下水開發利用率為155%，地下水超采嚴重。總體而言，廊坊市水資源極度匱乏。

廊坊市被京津二市的通州、武清分割成南北互不相連的兩部分，總體為2個水資源分區，北部為三河市、大廠縣、香河縣(俗稱北三縣)；南部為廣陽區、安次區、固安縣、永清縣、霸州市、文安縣、大城縣(俗稱南七縣)，共計10個行政分區。將城鎮生活生產以及農村生活生產按照鎮級行政區劃分為獨立的計算單元，共計91個。廊坊市內無水庫，因此將廊坊市內水閘均作為獨立水庫節點，共計39個。為了便于在水資源分區進行耗水平衡分析，各水資源分區進出口斷面均設置水資源分區斷面控制節點，根據實際供水和受水單元狀況確定供水節點，共確定控制節點48個。根據以上分析和概化，確定系統網絡圖包括2個水資源分區、10個行政分區、91個計算單元、39個水庫節點、48個控制節點等，具體見圖2。

3.1 評價指標選取與篩選

為了更加科學、系統地進行綜合評價，滿足指標體系的先進性、可測性、層次性原則，在自然水系連通形態指標的基礎上，即水系環度α、節點連接率β、水系連通度γ3個評價指標的基礎上，從水資源優化配置系統網絡圖中對計算單元水傳輸系統(線)的定義出發，根據人工側支水循環系統所包括的水傳輸關系，結合生產生活用水水平評價指標，選取城市人均綜合生活用水量[L/(人·d)]、單位GDP耗水(m3/萬元)、管網漏損率(%)、城市污水處理率(%)、城市污水再生利用率(%)這5個指標。其中水系環度α、節點連接率β、水系連通度γ是從水系之間連通形態的角度出發，運用景觀生態學中描述廊道之間連通性的指標，來反映自然水系網絡圖中每個節點的連通性能強弱以及物質能量交換能力大小。而生產生活用水水平評價指標則可以從用水、取水、排水等方面描述人工側支水循環水系連通度情況。這8個指標結合成一套用于評價自然—人工二元水系連通度的綜合評價指標體系，其具體指標及說明見表1。

表1 評價體系指標Tab.1 Index of evaluation system

至于評價標準，為使評價指標具有可比性，本文參考鄧紅霞[21]對于評價指標的選取規則：①對于“百分率”和“單位均值”類型指標，可參考國內不同城市以及全國平均數值，基于發展速度的社會、經濟發展水平，制訂出適合于不同發展時期的Ⅰ級(優)、Ⅱ級(中)和Ⅲ級(差)3級標準；②對于無法用上述類型表示的指標，可結合當地發展規劃或全國平均水平，制訂Ⅰ級(該地區遠期目標值)、Ⅱ級(近期目標值)和Ⅲ級(現狀值)3級標準。指標對應等級與來源見表2。

表2 指標等級與來源Tab.2 Indicator levels and sources

為了確定評價指標相關性，對樣本數據進行KMO檢驗與巴特利特檢驗，結果見表3。

表3 KMO和巴特利特檢驗Tab.3 KMO and Bartlett test

通過軟件分析得到KMO的值為0.768，大于0.7為可接受范圍，巴特利特檢驗的F值為402.967，P=0<0.05，說明樣本數據適合采用主成分分析。通過主成分分析，相關系數矩陣的特征值及方差貢獻率如表3所示，相關系數矩陣特征向量如表4所示。通過表4我們可以看到，提取所有成分的解釋總方差為86.878%，大于85%，說明解釋能力很好。因此，直接計算相應的特征向量值與指標權重，見表5。

表4 相關系數矩陣的特征值及方差貢獻率Tab.4 Eigenvalue and variance contribution rate of correlation coefficient matrix

通過表5指標權重可以看到排名靠前的指標有管網漏損率X3、城市污水回用率X5、城市污水處理率X4，這3個指標均用于描述人工側支水循環水系連通發育程度，說明在自然-人工二元水系中，人類活動對二元水系連通的影響較大，對于水系連通形態與功能的研究應不單只專注于自然水系。

表5 相關系數矩陣特征向量Tab.5 Eigenvector of correlation coefficient matrix

3.2 實證分析

為了驗證水系連通度評價指標體系和相關權重的有效性和可行性，結合調水工程的建設啟用時間，選取過去(2006年、2010年)、現狀(2015年、2020年)、未來(2025年、2030年、2035年)共7個年份，不同年份對應的綜合評價值見表6，數據來源于《廊坊市水資源公報》、《水資源三調廊坊匯總》、《廊坊經濟統計年鑒2017》、《廊坊市節水型社會建設“十二五”規劃》、《廊坊市水利綜合規劃(2013年版)》、《廊坊水文水資源實用手冊》、《廊坊市城市總體規劃(2012-2030)》、《廊坊市生態綜合規劃報告》、《廊坊水文水資源實用指南》。

表6 水系連通度指標綜合評價值Tab.6 Comprehensive evaluation value of water system connectivity index

根據公式(2)求得8個指標在7個年份的同一度、對立度和差異度值，結果如表7。根據式(3)、(4)計算得到各個評價指標對各個年份的聯系度及其聯系度均值，計算結果見表8。

表7 各個年份的聯系度、同一度、差異度、對立度Tab.7 The degree of connection,identity,difference and opposition of each year

表8 各個評價指標對各個年份的聯系度及聯系度均值Tab.8 Connection degree and mean value of each evaluation index to each year

取i=0，j=-1以及表8中所得的聯系度均值，代入公式(2)，得到各個年份對應的聯系主值數，將得到的聯系主值數與評價等級進行對比，得到該年份水系連通度對應等級，見表9。最終結果為在2006年、2010年、2015年這3個時間節點上，廊坊市水系連通度等級為Ⅲ級，處于較差狀態；在南水北調中線啟動之后，廊坊市2020年、2025年水系連通度等級升為Ⅱ級，處于良好狀態；遠期南水北調東線啟動將使得廊坊市水系連通度進一步提升，水系連通度在2030、2035年處于Ⅰ級。按照本文構建的評價指標體系和指標權重得到的廊坊市過去、現狀、未來的水系連通度與《廊坊市十三五規劃》、《廊坊市生態綜合規劃》對廊坊市現狀調查與未來規劃所評價結果基本一致，說明本文構建的評價體系具有可靠性和有效性。

表9 聯系主值數與等級評價結果Tab.9 Relation between the number of main values and the result of grade evaluation

4 結 語

(1)以原有的自然水系連通評價體系為基礎，基于二元水循環水資源配置網絡圖中對計算單元水傳輸系統(線)的定義，依據科學性、目的性、全面性和代表性，選出了8個指標，從天然循環與人工側支循環兩方面，結合自然水系連通形態指標與生產生活用水水平評價指標，構建自然-人工二元結合水系連通度指標評價體系，用于評價隨著經濟社會發展與水利工程的加入其水系連通度變化。

(2)為了選取和建立自然-人工二元結合水系連通度指標評價體系，本文選用主成分分析法和集對分析法進行分析。主成分分析法可對所選取的指標之間相關性進行分析篩選，使得較少指標包含多數指標的信息，提高分析效率，并且可以克服在指標權重確定方面的主觀判斷。集對分析可以分析水系連通度評價指標間的同異反度，得出不同時段水系連通度所對應等級，結果顯示廊坊市過去、現狀、未來的水系連通度與現狀調查及未來規劃所評價結果基本一致，說明本文構建的評價體系具有可靠性和有效性，可以為之后的相關研究提供參考借鑒。

(3)但研究工作也存在一些不足，由于社會水循環渠系管網數量龐大且難以提取，本文所采用的指標可能無法全面描述人工側支循環水系網絡的具體情況，今后應從增加樣本數據長度、細化評價體系指標等方面來優化該評價指標體系，提高研究成果的精度和實用價值。

□