胡鐵松+張兵堂+程曉峰+熊威

摘要：提出了水量分配方案后評價指標(biāo)體系，建立了基于非負矩陣分解的水量分配方案后評價模型，并以1983年國務(wù)院批復(fù)的引灤水量分配方案為研究對象進行了水量分配方案診斷式后評價。該方法利用非負矩陣分解算法將指標(biāo)矩陣最大程度分解為基向量和權(quán)向量，根據(jù)權(quán)向量元素的數(shù)值對水量分配方案進行評級，根據(jù)基向量元素的數(shù)值判斷評價指標(biāo)的實施狀況。通過對引灤水量分配方案質(zhì)量、實施過程、實施效果和可持續(xù)性的調(diào)查分析，得出方案綜合后評價結(jié)果為優(yōu)，但在來用水可持續(xù)性、對水環(huán)境影響及水質(zhì)目標(biāo)保證率等方面需加以改進。評價結(jié)果表明，基于非負矩陣分解原理的后評價方法評價效果較好，可為水量分配方案等項目的后評價提供科學(xué)依據(jù)。

關(guān)鍵詞：水量分配；后評價；指標(biāo)體系；非負矩陣分解；引灤水量分配

中圖分類號：TV213 文獻標(biāo)識碼：A 文章編號：1672-1683（2017）03-0013-07

Abstract：This paper established an evaluation index system for water allocation and developed a new post-evaluation model based on non-negative matrix factorization （NMF） for diagnostic evaluation.The model used NMF to decompose the data matrix into base vectors and weight vectors to the greatest extent.Then the base vector was used to assess the evaluation indexes，and the weight vector was used for grading the evaluation scheme.This model was applied to evaluate the water allocation scheme for Luanhe River in terms of scheme quality，implementation process，implementation effect，and sustainability.The scheme was graded as excellent，but needs improvement regarding water sustainability，impact on the water environment，and guarantee rate of water quality targets.Results showed that the post-evaluation model based on NMF can produce satisfactory outcomes and can provide a scientific basis for post-evaluation of water allocation schemes.

Key words：water allocation；post-evaluation；index system；non-negative matrix factorization；Luanhe River water diversion

水量分配是指對水資源可利用總量或可分配水量向行政區(qū)域進行的逐級分配，其方案制定與實施關(guān)系到區(qū)域經(jīng)濟、社會和環(huán)境的綜合效益[1]。但由于水文情勢和社會經(jīng)濟需求的不斷變化，水量分配方案的制定并非是一蹴而就的。如黃河、科羅拉多河和墨累-達令河等流域的水量分配方案在制定與實施過程中都存在著諸多問題，需要進行多次評價與完善[2-4]。開展水量分配方案后評價目的是總結(jié)方案制定實施過程中的經(jīng)驗和教訓(xùn)，診斷存在的問題及原因，促進水量分配方案的進一步落實、修正和完善。

后評價牽涉的指標(biāo)數(shù)量多，評價信息龐雜，評價過程中常遇到三大難題：其一，評價過程常涉及高維數(shù)據(jù)的分析處理；其二，難以從數(shù)量較多的評價指標(biāo)中準(zhǔn)確分析出優(yōu)勢和劣勢指標(biāo)，從而進行診斷性后評價；其三，難以確定合適的指標(biāo)權(quán)重。為此，許多學(xué)者對綜合評價方法的理論和實踐進行了研究，評價方法從最初的前后對比法、成功度法、邏輯框架法發(fā)展到層次分析法、理想點法、模糊綜合評價法、灰色系統(tǒng)評價法和神經(jīng)網(wǎng)絡(luò)法等[5-10]，其中被廣泛應(yīng)用的是理想點法和層次分析法，但在應(yīng)用中均受一定的限制。其中，理想點法在進行后評價時，需給定理想點和權(quán)重，且只能得出最終的綜合評價結(jié)果，無法分析某項指標(biāo)的優(yōu)劣，因而只知評價結(jié)果，無法進行診斷性后評價；層次分析法需通過兩兩對比分析得出指標(biāo)權(quán)重，在應(yīng)對數(shù)目較大的后評價指標(biāo)體系時計算復(fù)雜，缺乏效率。

非負矩陣分解是一種高效處理矩陣數(shù)據(jù)的方法，能夠大幅消減原始矩陣的維數(shù)，從而提取出矩陣中的關(guān)鍵信息[11]。目前非負矩陣分解在很多領(lǐng)域都得到了廣泛的應(yīng)用，Lee和Seung首次提出NMF理論時便將其用于人臉識別[12]，其后M Arngren[13]、F Kaiser[14]、Y Li[15]等分別將其運用到圖像處理、聲音識別、數(shù)據(jù)挖掘等領(lǐng)域。水量分配方案后評價指標(biāo)數(shù)量較多，評價矩陣信息龐雜，因此本文基于非負矩陣分解原理建立一種新的后評價模型，使評價過程更為簡便，同時能夠?qū)崿F(xiàn)診斷式后評價。

1 基于NMF原理的后評價方法

非負矩陣分解原理（Non-negative matrix factorization，NMF）是由Lee和Seung1999年在《Nature》上提出的一種新的矩陣分解算法，與其他特征抽取算法相比，該算法對分解結(jié)果非負的限制具有實現(xiàn)上的簡便性、分解形式和分解結(jié)果上的可解釋性等優(yōu)點[16-17]。

NMF算法用數(shù)學(xué)語言可描述為[18]：對于一個n×m階的非負矩陣V，可將其分解為一個n×r階的非負矩陣W和一個r×m階的非負矩陣H的乘積，即V=W×H。

假設(shè)水量分配方案后評價中有n個評價指標(biāo)和m個評價等級，請k名專家對n項指標(biāo)進行評分，統(tǒng)計各評價等級專家人數(shù)占比，由此可建立評價數(shù)據(jù)矩陣V：

根據(jù)非負矩陣的分解原理，當(dāng)r=1時矩陣V可近似分解為非負列向量W與非負行向量H的乘積。為最大限度地分解矩陣V，Lee和Seung[12]提出了一種基于歐氏距離的目標(biāo)函數(shù)：

當(dāng)且僅當(dāng)，W和H是穩(wěn)定點時，歐氏距離最小且不再變化。該迭代公式具有收斂性和解的存在性。

分析式（3）可知，當(dāng)歐氏距離E達到最小值時，W和H均不具備唯一性。當(dāng)H中的所有元素同時放大k倍，W中所有元素會變成原來的1/k。因此，為保證歐式距離最小時非負分解的唯一性，增加如下獨立條件：

矩陣分解偏度δ*表示歐氏距離最小時矩陣分解后的信息丟失程度。偏度δ*越小，矩陣分解后的信息丟失越少，分解結(jié)果越真實可信。當(dāng)分解偏度δ*小于01時，認為分解結(jié)果是可以接受的。

在求得列向量W*=（w*1，w*2，…，w*n）T與行向量H*=（h*1，h*2，…，h*m）后，由于基向量集合所有專家對各項指標(biāo)的評價等級信息，權(quán)向量綜合了各位專家對所有指標(biāo)評價等級信息，因此可以用基向量元素的大小衡量各項指標(biāo)的優(yōu)劣，用權(quán)向量元素的大小判斷水量分配方案的綜合評價等級。

2 水量分配方案后評價指標(biāo)體系的構(gòu)建與評價標(biāo)準(zhǔn)2.1 水量分配方案后評價指標(biāo)體系的構(gòu)建

評價指標(biāo)體系的建立是后評價工作的基礎(chǔ)。目前國內(nèi)尚缺乏用于水量分配方案后評價的指標(biāo)體系，因此需借鑒其他項目后評價指標(biāo)體系的建立思路形成初選指標(biāo)體系，并采用德爾菲法對其進行篩選，其建立流程見圖1。

圖1中立法后評價是針對已通過并實施的法律進行評估，需評估法律實施過程、效果和立法質(zhì)量[19-20]；而建設(shè)項目后評價的指標(biāo)體系依據(jù)行業(yè)特性的不同而不同，國家計委要求重點建設(shè)項目的后評價內(nèi)容需包含過程評價、經(jīng)濟效益評價、影響評價和持續(xù)性評價[21-23]。水量分配方案和立法均是通過條文發(fā)揮效用，因此需評估方案本身的質(zhì)量、實施過程和效果；同時和立法不同是水量分配方案的既定目標(biāo)是否能長期持續(xù)會受到氣候變化、人類活動等因素的影響，因此需評價方案的可持續(xù)性[24]。

根據(jù)以上分析初步制定了水量分配方案后評估指標(biāo)體系并采用德爾菲法篩選確定，最終形成了以方案質(zhì)量、實施過程、實施效果及其可持續(xù)性為準(zhǔn)則的水量分配方案后評估指標(biāo)體系，該指標(biāo)體系由1個目標(biāo)、4個準(zhǔn)則、10個內(nèi)涵和24項指標(biāo)構(gòu)成（表1）。指標(biāo)體系的設(shè)計是基于目標(biāo)導(dǎo)向的構(gòu)造原則，即要求水量分配能更好地保護生態(tài)、節(jié)約水資源和促進水資源可持續(xù)利用。

2.2 單向指標(biāo)評價標(biāo)準(zhǔn)

指標(biāo)體系建立后，為統(tǒng)一標(biāo)準(zhǔn)，在具體應(yīng)用時需制定單項指標(biāo)評價標(biāo)準(zhǔn)。標(biāo)準(zhǔn)依據(jù)指標(biāo)的評估目標(biāo)和相關(guān)的文件制定。限于篇幅對單項指標(biāo)評價標(biāo)準(zhǔn)的制定不做一一贅述，僅選用一項定性指標(biāo)和一項定量指標(biāo)為例加以說明。其中定性指標(biāo)選用“分配要素合理性”為例。按水利部頒布的《水量分配暫行辦法》及《水量分配方案制定技術(shù)大綱》 對分配要素的要求，以是否考慮水質(zhì)、是否考慮地下水、是否考慮生態(tài)用水、分配對象是耗損水量還是用水量為考核因素，推薦水量分配要素合理性評價標(biāo)準(zhǔn)如下。

（1）方案給出地表水耗損量及相應(yīng)的用水量分配成果，進行水質(zhì)水量雙控制，并考慮河道內(nèi)生態(tài)用水，則分配要素合理性為優(yōu)。

（2）方案給出地表水耗損量分配成果，并考慮河道內(nèi)生態(tài)用水，則分配要素合理性為良。

（3）方案給出地表水耗損量分配成果，則分配要素合理性為中。

（4）方案僅給出地表水取用水量分配成果，則分配要素合理性為差。

定量指標(biāo)選用“年分配水量偏離系數(shù)”為例。當(dāng)申請用水小于分水指標(biāo)時，偏離系數(shù)是指實際供水和申請用水的偏差值，以檢驗是否滿足需求或造成資源浪費；當(dāng)申請用水大于分水指標(biāo)時，偏離系數(shù)是指實際供水和分水指標(biāo)的偏差值，以衡量分配方案的約束效用。即分配水量偏差系數(shù)計算公式如下：

式中：i為年份；Si為第i年實際供水；yi為第i年申請用水；Fi為第i年分水指標(biāo)；pi為第i年偏差。

在具體評價時，推薦以下評價標(biāo)準(zhǔn)：若偏差系數(shù)小于10%，評價級別為優(yōu)；偏差系數(shù)為10%～20%，評價級別為良；偏差系數(shù)為21%～50%，評價級別為中；偏差系數(shù)大于50%，評價級別為差。

3 引灤水量分配方案后評價

3.1 引灤工程水量分配方案

自1983年國務(wù)院批復(fù)引灤水量分配方案以來，該方案已實施30余年，為天津市和河北省創(chuàng)造了巨大的經(jīng)濟、社會和環(huán)境效益[25]。但與此同時該方案還存在諸多問題，如在分配水量時未考慮下游生態(tài)用水等，同時受氣候變化與人類活動影響，相同保證率條件下潘家口和大黑汀水庫的入庫水量和可供水量已遠低于1983年方案批復(fù)時[26]。因此，有必要進行引灤水量分配方案診斷式后評價，為水量分配方案的修訂提供依據(jù)。

引灤水量分配對象是潘家口水庫水量和潘家口至大黑汀水庫的區(qū)間來水，其中潘家口水庫水量分配見表2。

潘大區(qū)間水量分配依據(jù)1983年印發(fā)的《引灤水量分配與供水調(diào)度管理辦法》：當(dāng)潘大區(qū)間年來水量等于或大于2.2億m3時，大黑汀水庫年供水量為2.1億m3；當(dāng)區(qū)間年來水量小于2.2億m3時，年供水量為區(qū)間年來水量減去庫損0.12億m3。

3.2 方案后評價計算

評價標(biāo)準(zhǔn)制定后，選用長期從事水量分配方案制定的3位專家和4位從事引灤水量分配工作的地方工作者組成專家組，根據(jù)方案的實施情況和單項指標(biāo)評價標(biāo)準(zhǔn)對引灤水量分配方案進行評價，評價等級為優(yōu)、良、中、差。統(tǒng)計結(jié)果中各等級的專家人數(shù)占比，可得數(shù)據(jù)矩陣V。對矩陣V采用維度為1的非負矩陣分解，可確定各評價指標(biāo)對應(yīng)的基元素數(shù)值，見表3；同時可得出各評價級別對應(yīng)的權(quán)向量元素數(shù)值，見表4；計算指標(biāo)體系中準(zhǔn)則層所對應(yīng)評價指標(biāo)的基元素均值，見表5。 3.3 方案后評價結(jié)果分析

由于指標(biāo)矩陣的分解偏度小于0.1，因此評價結(jié)果是可以接受的。

為分析求得的基向量和權(quán)向量的含義，分別作基元素數(shù)值與單項指標(biāo)專家打分均值相關(guān)關(guān)系圖和權(quán)向量元素數(shù)值與評價等級專家人數(shù)占比之和的相關(guān)關(guān)系圖見圖2、圖3。

由圖2可知，基元素數(shù)值與專家打分均值的相關(guān)系數(shù)為0978 3，因此基向量中元素數(shù)值大小可以代表評價指標(biāo)的優(yōu)劣；同時由圖3可知，權(quán)向量中元素數(shù)值與級別專家占比之和相關(guān)系數(shù)為0980 6，因此權(quán)向量元素大小可衡量水量分配方案的綜合評價等級。由表4可知，后評價級別中“優(yōu)”對應(yīng)的權(quán)向量元素值最大，因此引灤水量分配方案的評價等級為“優(yōu)”。

基于上述計算結(jié)果，分析引灤水量分配各單項指標(biāo)執(zhí)行情況和四個準(zhǔn)則的評價結(jié)果，可得出以下結(jié)論。

（1）方案質(zhì)量評價。

評價結(jié)果為良。

由表5可以看出，該準(zhǔn)則層對應(yīng)的指標(biāo)層元素均值為0046 3，在4個準(zhǔn)則中排名第二。分析表3中對應(yīng)數(shù)據(jù)可知，這首先得益于方案的編制和審批均合理合法，符合《水量分配暫行辦法》的相關(guān)要求；其次，方案與灤河流域早期規(guī)劃的銜接性良好，此后《灤河流域綜合規(guī)劃報告》和《海河流域水資源綜合規(guī)劃》等的制定也是基于1983年的引灤水量分配方案；最后，確定的區(qū)域間分水比例符合天津市和河北省的實際情況，用水公平性良好。

不足之處有兩點，一是方案未考慮河道內(nèi)生態(tài)用水，因此分配要素合理性評分略低；二是與三條紅線控制指標(biāo)銜接性較差，這是由于方案編制之初三條紅線控制指標(biāo)還未設(shè)立，因此需對水量分配方案進行相關(guān)修正。

（2）方案實施過程評價。

評價結(jié)果為優(yōu)。

由表5可以看出，其對應(yīng)的指標(biāo)層元素均值為0048 6，是4項準(zhǔn)則層中最高的。分析表3中對應(yīng)數(shù)據(jù)可知，這首先得益于引灤局制定了詳細的供水計劃和調(diào)水方案，保障了方案實施。其次，來水量預(yù)測精度較高，這與實際調(diào)查分析的結(jié)果是一致的，以調(diào)度年2012年-2013年為例，潘家口水庫來水量預(yù)測誤差為58%，潘大區(qū)間來水量預(yù)測誤差為104%，均滿足《水文預(yù)報規(guī)范》對預(yù)測精度的要求。同時，確定與下達的分水指標(biāo)較為合理。2003年-2015年分水指標(biāo)與方案的平均偏離系數(shù)僅為623%。最后，水量分配監(jiān)控體系完整，控制斷面設(shè)置合理且完善。

（3）方案總體實施效果評價。

評價結(jié)果為中。

由表5可以看出，其對應(yīng)的指標(biāo)層均值為0035 8，低于方案質(zhì)量和實施過程的評價結(jié)果。實施效果中較好的方面包括：一是規(guī)范了區(qū)域的水資源管理制度，明確了灤河流域的水權(quán)歸屬，對區(qū)域的水資源管理有良好的示范作用。二是促進了區(qū)域的經(jīng)濟社會良好發(fā)展，如引灤供水緩解了天津、唐山兩地區(qū)的用水緊張局面，取得了較大的經(jīng)濟效益和社會效益。第三，天津市多年平均分配水量偏差系數(shù)為1207%，河北省為2160%，可以評級為良好，這與計算結(jié)果是一致的。

引灤水量分配實施的不利效果是其對灤河下游水生態(tài)及環(huán)境的不利影響。引灤工程中的潘家口、大黑汀水庫建成蓄水后，一直采用滴水不漏的運行方式，導(dǎo)致下游河道干涸，地下水位大幅下降，對沿河各縣的經(jīng)濟社會和生態(tài)環(huán)境有較嚴(yán)重的負面影響。

（4）方案可持續(xù)性評價。

評價結(jié)果為差。

由表5可以看出，其對應(yīng)的指標(biāo)層均值為0034 7，是四個準(zhǔn)則層元素中最低的。這主要是因為隨著水文情勢和用水需求的不斷變化，入庫水量減少、用水量增加，可分配水量不具有可持續(xù)性。這與實際計算分析的結(jié)果是一致的，潘家口水庫2000年前后來水量情況如圖4所示，2000年-2013年系列均值比1980年-1999年系列減少了555%。

1980年-2013年天津、河北兩地區(qū)的歷年引灤水量見圖5。由圖5可知天津市引灤水量呈先增大后減小的趨勢，但變化幅度不大；河北省年引灤水量呈顯著下降趨勢?？傮w上引灤工程的年供水量呈下降趨勢。因此引灤工程用水不具有可持續(xù)性。

相比之下，引灤水量分配政策制度的可持續(xù)性則較好。引灤水量分配分工明確，實施機構(gòu)完善，組織管理體系清晰，且有應(yīng)對特殊情況的協(xié)商制度。同時，管理制度可持續(xù)性也較好，其實施過程中的管理制度包括計劃用水制度、審批制度、監(jiān)督考核制度、應(yīng)急響應(yīng)制度和合同管理制度等。

4 結(jié)論與展望

4.1 結(jié)論

（1）本文提出了基于非負矩陣分解的水量分配方案后評價方法，該方法采用非負矩陣法對評價矩陣進行一維非負分解，并根據(jù)得到的權(quán)向量元素大小確定后評價等級，依據(jù)基向量元素大小衡量各評價指標(biāo)的優(yōu)劣，避免了常見方法中需給定理想點和權(quán)重等的缺陷，計算過程簡單，意義明確。

（2）從方案質(zhì)量、實施過程、實施效果和可持續(xù)性四個方面對引灤水量分配方案進行了系統(tǒng)分析，得出引灤水量分配方案綜合評價結(jié)果為“優(yōu)”的結(jié)論，這主要得益于方案具有合理性和可操作性，實施過程嚴(yán)謹(jǐn)且促進了區(qū)域水資源管理。同時指出了引灤水量分配方案及其實施過程中存在的不足，如用水和來水可持續(xù)性差，水質(zhì)目標(biāo)達標(biāo)率差和對區(qū)域的水環(huán)境有不利影響等。

4.2 展望

根據(jù)本文研究，建議今后在引灤水量分配方案修訂及其實施過程中注意以下問題：首先應(yīng)根據(jù)近些年的來水情況和用水情況重新計算可分配水量和分水比例；其次是加強對水環(huán)境和水生態(tài)的保護，在實施過程中需考慮水質(zhì)要求；三是要修正方案使其符合三條紅線的要求。

同時，本文的研究盡管有其理論和實踐意義，但不可避免地存在一些不足之處。首先，本文提出的指標(biāo)體系主要是基于前人的研究成果，缺乏指標(biāo)必要性和有效性的論證；其次，在制定單項指標(biāo)評價標(biāo)準(zhǔn)時較為主觀，不夠嚴(yán)謹(jǐn)；最后，在應(yīng)用非負矩陣分解法進行后評價時，缺乏對該方法優(yōu)劣性的評價。在未來的研究中，還需要進一步完善。

參考文獻（References）：

[1] 呂蘭軍.《水量分配暫行辦法》給水文帶來的機遇與挑戰(zhàn)[J].水利發(fā)展研究，2008（8）：44-45.（LV Lan-jun.Opportunities and challenges for hydrology by Interim Measures For Water Distribution[J].Water Resources Development Research，2008（8）：44-45.（in Chinese））DOI：10.3969/j.issn.1671-1408.2008.08.005

[2] 楊彥明.科羅拉多河水權(quán)分配的啟示[J].水利發(fā)展研究，2004，4（9）：46-49.（YANG Yan-ming.Enlightenment of water right allocation in Colorado River[J].Water Resources Development Research，2004，4（9）：46-49.（in Chinese）） DOI：10.3969/j.issn.1671-1408.2004.09.014

[3] 馬建琴，劉杰，夏軍，等.黃河流域與澳大利亞墨累一達令流域水管理對比分析[J].河南農(nóng)業(yè)科學(xué)，2009（7）：69-73.（MA Jian-qin，LIU Jie，XIA Jun，et al.The comparative research on water management institution between the Yellow River Basin and Australian Murray-Darling Basin[J].Journal of Henan Agricultural Sciences，2009（7）：69-73.（in Chinese）） DOI：10.3969/j.issn.1004-3268.2009.07.019

[4] 彭祥，胡和平.黃河可供水量分配方案的制度評價[J].人民黃河，2006，28（4）：41-43.（PENG Xiang，HU He-ping.The institution evaluation on the distribution scheme of water supply in the Yellow River[J].Yellow River，2006，28（4）：41-43.（in Chinese））DOI：10.3969/j.issn.1000-1379.2006.04.019

[5] 劉麗艷，吳玉芹，郭振苗.淺談節(jié)水灌溉項目后評價[J].中國水利，2004（9）：39-41.（LIU Li-yan，WU Yu-qin，GUO Zhen-miao，Evaluation on water-saving irrigation[J].China Water Resources，2004（9）：39-41.（in Chinese）） DOI：10.3969/j.issn.1000-1123.2004.09.012

[6] 葉珍.基于AHP的模糊綜合評價方法研究及應(yīng)用[D].廣州：華南理工大學(xué)，2010.（YE Zhen.Study and application of fuzzy comprehensive evaluation based on AHP.Guangzhou：South China University of Technology，2010.（in Chinese））

[7] 段云龍，周靜斌，申曉靜.基于熵權(quán)TOPSIS法的房地產(chǎn)項目后評價模型研究[J].項目管理技術(shù)，2011（9）：40-44.（DUAN Yun-long，ZHOU Jing-bin，SHEN Xiao-jing.Evaluation of estate project based on entropy weight TOPSIS method[J].Project Management Technology，2011（9）：40-43.（in Chinese）） DOI：10.3969/j.issn.1672-4313.2011.09.008

[8] 徐健，吳瑋，黃天寅，等.改進的模糊綜合評價法在同里古鎮(zhèn)水質(zhì)評價中的應(yīng)用[J].河海大學(xué)學(xué)報：自然科學(xué)版，2014，42（2）：143-149.（XU Jian，WU Wei，HUANG Tian-yin et al.Application of improved fuzzy comprehensive evaluation to water quality evaluation in Tongli Town[J].Journal of Hohai University：Natural Sciences，2014，42（2）：143-149.（in Chinese）） DOI：10.3876/j.issn.1000-1980.2014.02.009

[9] 華瑤，周雨.邏輯框架法在電網(wǎng)建設(shè)項目后評價中的應(yīng)用[J].工業(yè)技術(shù)經(jīng)濟，2011，30（1）：97-102.（HUA Yao，ZHOU Yu.The application of logical framework approach in post evaluation of grid item[J].Journal of Industrial Technological Economics，2011，30（1）：97-102.（in Chinese）） DOI：10.3969/j.issn.1004-910X.2011.01.017

[10] 李元松，陳文峰，李新平，等.基于模糊神經(jīng)網(wǎng)絡(luò)的邊坡穩(wěn)定性評價方法[J].武漢理工大學(xué)學(xué)報，2013，35（1）：113-118.（LI Yuan-song，CHEN Wen-feng，LI Xin-ping et al.Stability assessment of rock slope based on fuzzy neural network[J].Journal of Wuhan University of Technology，2013，35（1）：113-118.（in Chinese）） DOI：10.3963/j.issn.1671-4431.2013.01.022

[11] 徐泰燕，郝玉龍.非負矩陣分解及其應(yīng)用現(xiàn)狀分析[J].武漢工業(yè)學(xué)院學(xué)報，2010，29（1）：109-114.（XU Taiyan，HAO Yu-long.The current research situation analysis of non-negative matrix factorization and applications[J].Journal of Wuhan Polytechnic University，2010，29（1）：109-114.（in Chinese）） DOI：10.3969/j.issn.1009-4881.2010.01.030

[12] Daniel D.Lee，H.Sebastian Seung.Learning the parts of objects by nonnegative matrix factorization[J].Nature，1999，401（6755）：788-791.DOI：10.1038/44565

[13] M Arngren，MN Schmidt，J Larsen.Unmixing of hyperspectral images using Bayesian non-negative matrix factorization with volume prior[J].Journal of Signal Processing Systems，2011，65（3）：479-496.DOI：10.1007/s11265-010-0533-2

[14] F.Kaiser，T.Sikora.Music structure discovery in popular music using non-negative matrix factorization[C]，Proceedings of International Society for Music Information Retrieval，2010：429-434.

[15] Y Li，A Ngom.The non-negative matrix factorization toolbox for biological data mining[J].Source Code for Biology and Medicine，2013，8（1）：10-24.DOI：10.1186/1751-0473-8-10

[16] 黃帥棟，衛(wèi)志農(nóng)，高宗和等.基于非負矩陣分解的相關(guān)向量機短期負荷預(yù)測模型[J].電力系統(tǒng)自動化.2012，36（11）：62-66.（HUANG Shuai-dong，WEI Zhi-nong，GAO Zong-he et al.A short-term load forecasting model based on relevance vector machine with nonnegative matrix factorization[J].Automation of Electric Power Systems，2012，36（11）：62-66.（in Chinese）） DOI：10.3969/j.issn.1000-1026.2012.11.011

[17] 林濱，張彥斌，陳善工.基于非負矩陣分解算法的工程陶瓷磨削表面損傷檢測[J].光學(xué)精密工程，2012，20（11）：2459-2464.（LIN Bin，ZHANG Yan-bin，CHEN Shan-gong.Damage detection of en neering ceramics ground surface based on NMF[J].Optics and Precision Engineering，2012，20（11）：2459-2464.（in Chinese）） DOI：10.3788/OPE.20122011.2459

[18] 王文俊，張軍英.一種非負矩陣分解的快速方法[J].計算機工程與應(yīng)用，2009，45（25）：1-2.（WANG Wen-jun，ZHANG Jun-ying.Fast algorithm for nonnegative matrix factorization[J].Computer Engineering and Applications，2009，45（25）：1-2.（in Chinese）） DOI：10.3778/j.issn.1002-8331.2009.25.001

[19] 孫樹曼.我國立法后評估指標(biāo)體系的構(gòu)建[D].濟南：山東大學(xué)，2009.（SUN Shu-man.The construction of the post-assessment index system of China[D].Jinan：Shandong University，2009.（in Chinese））

[20] 王保民，崔東曉.歐盟立法評估制度研究[J].行政與法，2011，0（5）：101-105.（WANG Bao-min，CUI Dong-xiao.A research on the legislative evaluation institution of European Community[J].Ublic Administration & Law，2011，0（5）：101-105.（in Chinese）） DOI：10.3969/j.issn.1007-8207.2011.05.025

[21] 安中仁，張文潔.水利建設(shè)項目后評價[J].水利建設(shè)與管理，2010，30（9）：81-84.（AN Zhong-ren，ZHANG Wen-jie.Post-evaluation of water conservancy construction project[J].Water Resources Development & Management，2010，30（9）：81-84.（in Chinese）） DOI：10.3969/j.issn.1005-4774.2010.09.035

[22] 薛敏，陸惠民.完善我國政府投資項目后評價體系研究[J].工程管理學(xué)報，2011，25（3）：318-322.（XUE Min，LU Hui-min.Research on improving the post-evaluation system of government investment projects in China[J].Journal of Engineering Management，2011，25（3）：318-322.（in Chinese）） DOI：10.3969/j.issn.1674-8859.2011.03.016

[23] 蔡文婷.政府投資項目后評價指標(biāo)體系及評價方法研究[D].沈陽：沈陽建筑學(xué)院，2011.（CAI Wen-ting.The research of the government investment project post evaluation index system and methods[D].Shenyang：Shenyang Jianzhu University，2011.（in Chinese））

[24] 許新宜.關(guān)于水量分配原則的幾點看法[J].南水北調(diào)與水利科技，2011，9（2）：15-17.（XU Xin-yi.Several Perspectives on the Principles for River Basin Water All ocation in China[J].South-to-North Water Diversion and Water Science & Technology，2011，9（2）：15-17.（in Chinese）） DOI：10.3742/SP.J.1201.2011.02015

[25] 王立林，王鴻艷，齊德海，等.引灤入津工程水資源現(xiàn)狀分析及解決措施[J].海河水利，2008（5）：9-11.（WANG Li-lin，WANG Hong-yan，QI De-hai et al.Water resources present situation analysis and solutions of water diversion from Luan River to Tianjin[J].Haihe Water Resources，2008（5）：9-11.（in Chinese）） DOI：10.3969/j.issn.1004-7328.2008.05.004

[26] 張利平，李凌程，夏軍，等.氣候波動和人類活動對灤河流域徑流變化的定量影響分析[J].自然資源學(xué)報，2015，30（4）：664-672.（ZHANG Li-ping，LI Ling-cheng，XIA Jun et al.Quantitative assessment of the impact of climate variability and human activities on runoff change in the Luanhe River Catchment[J].Journal of Natural Resources，2015，30（4）：664-672.（in Chinese）） DOI：10.11849/zrzyxb.2015.04.012