陳 月
(遼寧省丹東水文局,遼寧 丹東 118000)
最嚴格水資源管理對于促進經(jīng)濟社會可持續(xù)發(fā)展和解決復雜的水環(huán)境、水資源問題具有重要意義,并且可為水資源管理格局與規(guī)劃提供重要的參考依據(jù)。水體污染、用水效率低和水資源過度開發(fā)為目前最嚴格水資源管理面臨的主要問題,不僅明確了水資源保護、節(jié)約與優(yōu)化配置的任務要求和工作中心,而且體現(xiàn)了水資源利用的各個環(huán)節(jié)[1- 2]。ANN智能法為評價非線性、高緯度、多指標問題較為常用的方法,可對綜合復雜問題進行系統(tǒng)的處理。BP神經(jīng)網(wǎng)絡為ANN智能法應用較為普遍的模型之一,因有較好的非線性、高維度復雜問題的處理能力,被廣泛應用于水資源綜合評價方面。然而,BP網(wǎng)絡在實際應用過程中通常存在運算效率低、陷入局部最優(yōu)以及結構復雜等問題,SVM支持向量機是一種用于解決非線性映射和模式識別的前饋神經(jīng)網(wǎng)絡,可根據(jù)最小化風險旋球原理和學習統(tǒng)計中的VC維方法提升系統(tǒng)的泛化能力,由于具有嚴謹?shù)臄?shù)學基礎已成為計算機學習領域新的研究方向,SVM模型在回歸預測與模式識別方面具有廣泛的應用前景。RBF徑向基神經(jīng)網(wǎng)絡是一種能夠?qū)崿F(xiàn)高維空間線性可分的前饋型模型,其基本原理是對隱含層空間利用隱單元的“基”進行構造,從而實現(xiàn)低維模型的輸入數(shù)據(jù)與高維隱含層數(shù)據(jù)之間的轉(zhuǎn)換,在高維空間內(nèi)完成線性可分。RBF由于具有良好的逼近能力,并且可對各類非線性函數(shù)進行不同精度的逼近,從而有效避免了傳統(tǒng)網(wǎng)絡易陷入局部極值的問題,在類型分析、數(shù)據(jù)分類、時間序列分析以及非線性問題處理等領域具有廣泛的應用前景和價值[3- 11]。
本文在詳細分析了最嚴格水資源管理內(nèi)涵的基礎上,基于AHP層次分析建立了包含20項典型指標的評價體系,然后根據(jù)納污能力、用水效率、水資源開發(fā)以及用水總量控制紅線,將最嚴格水資源管理類型劃分為Ⅰ~Ⅴ級5個評價標準。通過深入探討SVM與RBF基本原理構建了評價模型,然后以遼寧省現(xiàn)狀年、規(guī)劃水平年為例驗證了模型的可信度與可行性,為區(qū)域水資源的開發(fā)利用和管理規(guī)劃控制提供一定理論支撐。
為保障水資源和經(jīng)濟社會的可持續(xù)健康發(fā)展我國必須實行最嚴格水資源管理制度,它是由水資源管理存在的突出問題、日趨嚴峻的水資源形勢、水環(huán)境持續(xù)惡化、用水矛盾突出以及水資源時空分布不均等基本國情所決定的。針對當前水體污染嚴重、粗放型水資源利用以及過度開發(fā)等突出問題提出的最嚴格水資源管理制度,其基本內(nèi)涵是基于“四項制度”與“三條紅線”監(jiān)理管理體系,為嚴格控制入河排污總量、遏制用水浪費以及用水總量控制建立相應的控制紅線。
我國于2009年就提出了最嚴格水資源管理的理念,然而在真正實行最嚴格水資源管理與各省市“制度出臺”之間還存在較遠的距離。系統(tǒng)、全面的指標體系是最嚴格水資源管理分析的重要步驟和關鍵內(nèi)容,影響可評價結果的因素較多,且因素之間通常具有模糊的、不確定性復雜關系。評價指標體系是由多個子系統(tǒng)構成的多層級復雜的整體,各層級之間相互作用、互相影響,具有十分復雜的關系。因此,指標的選取不僅要反映區(qū)域水資源特征、內(nèi)涵等基本內(nèi)容,并且能夠覆蓋與水資源管理的各個方面,同時考慮系統(tǒng)本身的發(fā)展狀況以及各子系統(tǒng)之間的協(xié)調(diào)作用。評價指標應盡可能的內(nèi)涵豐富、度量簡便、靈敏性強,并且能夠最大限度地反映系統(tǒng)特征,本文參考最嚴格水資源管理相關文獻和已有考核標準,在遵循代表性、系統(tǒng)性、可量化性、獨立性、層次性等原則的基礎上,分別從4個方面利用AHP層次分析法選擇20項典型指標。然后根據(jù)用水總量、用水效率、納污能力、水資源開發(fā)控制紅線,將最嚴格水資源管理類型劃分為Ⅰ~Ⅴ級5個評價標準見表1。
在沒有增加復雜度的條件下可利用SVM實現(xiàn)非線性變化的線性回歸,對非線性復雜問題進行線性空間的求解,從而較好地解決了傳統(tǒng)神經(jīng)網(wǎng)絡易出現(xiàn)的局部極值、高緯度非線性的小樣本復雜問題,其基本流程如下:
f(x)=wΦ(x)+b
(1)
式中,Φ(x)—低維空間中非線性映射函數(shù)。
表1 最嚴格水資源管理評價指標與分級標準
(2)
(3)
式中,f(x)、y—回歸函數(shù)返回的預測值及其對應的真實值;C—懲罰因子,訓練誤差大于不敏感線性損失ε的樣本懲罰程度隨C值的增大而增大,由于不敏感線性損失ε規(guī)定了誤差要求,因此回歸函數(shù)誤差隨ε值的降低而減少;對上述公式進行對偶形式轉(zhuǎn)換時可引入Lagrange函數(shù),具體的轉(zhuǎn)換過程詳見文獻。
(3)計算對偶函數(shù)的最優(yōu)解a與a*,然后根據(jù)計算結果可有:
(4)
(5)
式中,Nnsv—支持向量機個數(shù)。
(4)根據(jù)公式(1)并輸入w*、b*計算結果可有回歸函數(shù),其表達式為:
(6)
RBF是由輸出、隱含與輸入層組成的具有收斂速度快、訓練精度高、結構簡單、能夠根據(jù)不同精度逼近各類連續(xù)函數(shù)等優(yōu)點的模型,相對于BP網(wǎng)絡RBF模型具有較高的學習速度、較強的模式識別與逼近能力。RBF網(wǎng)絡選擇徑向基Gaussian函數(shù)時的激活函數(shù)如下:
(7)
可利用下式作為RBF網(wǎng)絡的輸出結果:
(8)
式中,yj—第j個節(jié)點的網(wǎng)絡實際輸出;wij,φi(x)—分別為連接權值以及第i個節(jié)點的輸出。
引入樣本的期望輸出為d,則可采用下式作為基函數(shù)方差σ:
(9)
RBF網(wǎng)絡的學習方法主要有正交最小二乘法、監(jiān)督選取中法、自組織以及隨機選取中心法等,本研究對RBF網(wǎng)絡學習選擇自組織選擇中心法,受文章篇幅限制不再詳細介紹,具體的計算流程詳見文獻。
(1)歸一化處理。各項評價指標的屬性與量綱存在差異,對評價結果的影響也不盡相同,因此在最嚴格水資源管理評價之前應先對各指標初始數(shù)據(jù)作歸一化處理。數(shù)值越大則評價結果越優(yōu)型指標的歸一化處理方法如公式(10),而對于數(shù)值越小則評價結果越優(yōu)型指標取歸一化處理后的倒數(shù)并乘以100,處理公式如下:
(10)
(2)設計輸入與輸出樣本。相對于傳統(tǒng)BP,網(wǎng)絡SVM是基于最小化風險原則解決小樣本回歸問題,因此在應用過程中具有明顯的優(yōu)質(zhì)。在各個評價等級閥值之間本研究利用隨機生成的方式構成10組樣本,然后隨機選擇5組訓練樣本,并將其他樣本作為檢驗樣本。因此,隨機抽取得到50組樣本,檢驗與訓練樣本各20組,不同模型的輸出模式見表2。
表2 各模型期望輸出與評價樣本
(3)統(tǒng)計性能評價指標。對傳統(tǒng)BP網(wǎng)絡、RBF模型以及LM-BP的性能選擇運行時間、決定系數(shù)R2、RMAPE均方根相對誤差、RMSE均方根絕對誤差、MAE平均絕對誤差、maxRE最大相對誤差以及MRE平均相對誤差7個評價指標進行分析,從而對網(wǎng)絡訓練過程中可能存在的“過擬合”與“欠擬合”問題進行控制,有助于提高模型的泛化能力。決定系數(shù)R2值在0~1區(qū)間,模型的性能隨趨近于1的程度的增大而增大;模型性能與其他指標值呈負相關性,即其他指標值越小則模型性能越優(yōu)。
(4)網(wǎng)絡訓練與評價。本研究利用libsvm工具箱以及Matlab語言分別建立傳統(tǒng)的BP網(wǎng)絡、RBF與SVR模型,然后對表2中各評價樣本與期望輸出進行檢驗與訓練,為獲取較好的評價精度進行多次反復調(diào)試,最終確定的各模型參數(shù)分別如下:
①SVR模型。對于給定的核函數(shù),影響評價精度的主要參數(shù)為模型核函數(shù)參數(shù)g及其懲罰因子C。結合相關文獻選擇核函數(shù)為徑向基函數(shù),設置的g及C參數(shù)的搜索空間相同,均為2-2~26,對參數(shù)g與C利用交叉驗證法確定,不敏感系數(shù)ε為0.001,進步程度均為0.5。
②傳統(tǒng)BP網(wǎng)絡。由于目前還未形成較為理想的BP網(wǎng)絡傳遞函數(shù)、訓練次數(shù)、期望誤差、隱層節(jié)點數(shù)等參數(shù)確定方法,因此主要還是憑借以往經(jīng)驗選取參數(shù)值。根據(jù)現(xiàn)有文獻最終確定為最嚴格水資源管理評級的BP網(wǎng)絡結構為20- 25- 1,訓練函數(shù)為traingdx,輸出與隱含層傳遞函數(shù)為logsig和tansig,最大迭代次數(shù)設定為1000次、期望誤差為0.0001時模型達到最優(yōu)。
③RBF模型。在徑向基函數(shù)中RBF網(wǎng)絡的期望誤差與擴散速度spread利用循環(huán)訓練算法程序確定,模型達到最優(yōu)性能時分別為0.0001和2。
選擇7項指標運行5次的平均值作為模型優(yōu)劣評判的主要依據(jù),根據(jù)上述各模型參數(shù)設定值最終的運算輸出,見表3。
表3 各模型隨機運行5次的優(yōu)劣性能指標輸出值
根據(jù)表3各指標輸出值可知,RBF和SVR模型的7項性能指標在隨機運行5次、隨機選取與生成樣本條件下的均值優(yōu)于傳統(tǒng)的BP網(wǎng)絡,由此表明在預測評價過程中RBF與SVR模型具有較強的泛化能力與較高的預測精度。本文從運行時間、泛化能力以及評價精度3個方面選擇RBF與SVR模型,并以遼寧省為例評價分析了該區(qū)域最嚴格水資源管理。
遼寧省總面積為14.8萬km2,主要包括沈陽、大連、撫順、鞍山、本溪、丹東等14個省轄市;屬于大陸性季風氣候,年降水量在600~1100mm之間,由于其獨特的地理特征和環(huán)境氣候,年際水資源量呈現(xiàn)出典型的周期性豐枯變化特征,其中70%以上的降水量集中在汛期6—9月,并多以暴雨和強降雨的形式出現(xiàn)。年均水資源量為341.79億m3,境內(nèi)主要包括大凌河、遼河、渾河、太子河和繞陽河等河流,東南部與渤海相鄰水源面積較大,降水量較為充足,而西北區(qū)域降水較少并且風沙干旱現(xiàn)象較為頻繁。該區(qū)域人均水資源占有量僅為全國人均占有量的8.6%,水資源儲蓄量嚴重不足,遼寧省屬于我國嚴重缺水區(qū)域,并且水資源供需問題隨著近年來人口的不斷增長、經(jīng)濟的快速發(fā)展以及需水量的大幅增加而日趨突出,再加上工農(nóng)業(yè)用水效率低、水體污染嚴重等問題,進一步降低了區(qū)域水資源承載能力。該區(qū)域最嚴格水資源管理主要有水資源緊缺、降水量時空分布不均、用水效率低、水資源利用方式粗放以及局部地區(qū)水環(huán)境惡化等問題。由于區(qū)域經(jīng)濟發(fā)展不平衡,其用水水平存在一定差異,因此研究該區(qū)域最嚴格水資源管理對于區(qū)域經(jīng)濟發(fā)展和水資源優(yōu)化配置具有重要意義[12- 14]。
根據(jù)遼寧省水資源公報、統(tǒng)計年鑒以及政府工作報告等資料提取現(xiàn)狀年各指標初始數(shù)據(jù),然后結合當?shù)亟?jīng)濟社會與水資源開發(fā)利用現(xiàn)狀分別確定規(guī)劃水平年2025、2030、2035年見表4。
對遼寧省2018、2025、2030、2035年最嚴格水資源管理分別利用已訓練好的RBF、SVR模型進行計算分析,然后將各評價等級閥值進行模擬運算,并以模型輸出結果作為等級評判依據(jù)見表5。
表4 遼寧省最嚴格水資源管理評價指標初始數(shù)據(jù)
表5 遼寧省最嚴格水資源管理輸出
根據(jù)上表計算結果可知,在2018、2025、2030、2035年基于RBF與SVR模型的遼寧省最嚴格水資源管理分別處于Ⅱ級、Ⅲ級、Ⅳ級與V級狀態(tài),SVR和RBF模型評價結果保持良好的一致性,評價結果較好地反映了遼寧省不同年份的水資源管理實際水平,可為該區(qū)域水資源管理規(guī)劃以及最嚴格水資源管理制度的實施提供一定的參考依據(jù);由于遼寧省2018年還處于“制度出臺”階段,根據(jù)模型輸出值可以看出,其管理考核制度在最嚴格水資源評價中處于空白狀態(tài),這也是影響評價結果的掛件因素;雖然在2025年評價結果為Ⅲ級,處于理想水平,但是模型輸出值接近理想狀態(tài),由此表明該區(qū)域?qū)嵤┑淖顕栏袼Y源管理已取得了明顯的成效但仍具有較大的提升空間;在2035年模型輸出為V級最理想狀態(tài),輸出結果處于Ⅳ級的上限水平,由此表明遼寧省水資源管理隨著最嚴格制度的實施逐漸趨于最理想水平。然而由于受到萬元GDP、農(nóng)業(yè)用水效率等因素影響,該區(qū)域仍具有較大的最嚴格水資源管理空間。
在我國水資源管理發(fā)展進程中,基于“思想制度”與“三條紅線”的最嚴格水資源管理制度的實施必將對水資源的永續(xù)發(fā)展與我國經(jīng)濟社會的可持續(xù)發(fā)展產(chǎn)生深遠而重要的影響。通過科學的方法評價分析最嚴格水資源管理的實施效果對于促進區(qū)域水資源利用良性循環(huán)、推動制度的實行具有非常重要的現(xiàn)實意義。本文在詳細分析了最嚴格水資源管理基本內(nèi)涵的基礎上,分別從4個不同方面選擇了20項典型指標并建立了評價體系與分級標準;然后對模型的評價性能與相關參數(shù)利用隨機選取法與隨機內(nèi)插法形成的數(shù)據(jù)樣本進行評價,在此基礎上對比分析了RBF、SVR以及傳統(tǒng)BP網(wǎng)絡模型的泛化能力和評價精度,得出的主要結論如下:
(1)RBF和SVR模型的7項性能指標在隨機運行5次、隨機選取與生成樣本條件下的均值優(yōu)于傳統(tǒng)的BP網(wǎng)絡。這兩個模型具有較強的泛化能力與較高的預測精度,可作為遼寧省最嚴格水資源管理評價模型。
(2)在研究期間基于RBF與SVR模型的遼寧省最嚴格水資源管理分別處于Ⅱ級、Ⅲ級、Ⅳ級與V級狀態(tài),遼寧省水資源管理隨著最嚴格制度的實施逐漸趨于最理想水平。然而由于受到萬元GDP、農(nóng)業(yè)用水效率等因素影響,該區(qū)域仍具有較大的最嚴格水資源管理空間。
(3)SVR和RBF模型評價結果保持良好的一致性,評價結果較好地反映了遼寧省不同年份的水資源管理實際水平,可為該區(qū)域水資源管理規(guī)劃以及最嚴格水資源管理制度的實施提供一定的參考依據(jù)。