0 引言

2013 年 11 月在中國共產黨十八屆三中全會的綱領性文件 《中共中央關于全面深化改革若干重大問題的決定》 中提出了“探索編制自然資源資產負債表，對領導干部實行自然資源資產離任審計，建立生態環境損害責任終身追究制”的重大決策[1],探索編制自然資源資產負債表是新時期加強生態文明體制機制、健全資源節約利用、建設美麗中國的根本戰略需求所在[2]. 水資源資產負債表是自然資源資產負債表在水資源主體上的具體應用，是今后落實最嚴格水資源管理制度、對領導干部實行自然資源資產離任審計和終身追究生態環境損害責任等的重要依據[3]. 水資源價格是編制水資產負債表的前提條件，如何對水資源賦予經濟價值是其核心與難點所在。 目前采用直接出售價格法、水權交易法、供水成本法、完全成本法、影子價格法、CGE 模型法等方法，但不同估價方法對同一資產估算結果存在不盡相同、可靠性以及可比性不強等缺陷，水資源價格核算仍然是一個充滿探索、實驗的研究領域[4-5].

水資源是生態環境可持續發展的決定性因素，在水上如何合理、高效地利用有限的水資源是張掖市發展需要解決的根本問題，其中合理制定水價是基本手段之一[6-7]. 本文綜合考慮水資源、水環境、經濟社會等影響因素，將熵權法和層次分析法相結合，采用模糊數學綜合評價模型，確定張掖市2010-2014年5 a水資源價值及價格動態變化趨勢，為水資產負債表的編制提供方法數據參考和支撐。

1 水資源資產價格核算指標體系

1. 1 研究區域

甘肅省張掖市位于祁連山北麓黑河流域中游地段，屬于典型的干旱綠洲農業區，張掖市集中了全黑河流域95%的耕地、90%以上的人口和89%的國內生產總值，是黑河流域中游社會經濟發展和水資源利用最具影響的區域[8],該區域具有干旱區內陸河流域的典型特征，平均多年人均水資源量為 1 350 m3,地均水量為 7 950 m3·km-2,分別占全國平均值的 60%和 30%,水資源短缺，是中國首批建設節水型社會試點城市之一[9].

1. 2 數據獲取

根據 2010-2014 年的 《張掖市統計年鑒》 和2010-2014 《張掖市水資源公報》，依據采用《國家地面水環境質量標準》（GB3838-2002）的污染程度分類制定水質等級標準以及 2010 年張掖市環境質量報告書，借鑒文獻有關研究成果[10-12],指標選取代表性好、針對性強，易于量化，選擇水量、水質、人口、社會經濟、生態環境5個方面的準則層和 18 個指標作為張掖市水資源價格評價指標體系（表 1）。 各指標的分級情況可根據各評價因素的實際數值，對照各因素的分級指標推求，分為高、較高、中等、較低、低5類。

2 水資源價格模糊綜合評價

模糊綜合評價的基本思想是應用模糊關系合成的原理，根據被評價對象本身存在的形態或類屬上的亦此亦彼性，從數量上對其所屬給以刻畫和描述[13]. 鑒于水資源系統的復雜性和模糊性，用模糊數學的概念和方法對水資源資產價值進行動態定量評價，由水資源價值評價模型和水資源價格計算模型兩部分組成。

2. 1 水資源價值模糊綜合評價模型

（1） 建立評價對象的指標集 V = f {X1,X2,X3,…，Xn},水資源價值要素集 V = {水量，水質，人口，社會經濟，生態環境};

（2） 建立評價集 U = {μ1,μ2,μ3,…，μn},水資源價值的評價向量 W={高，較高，中等，較低，低};

（3） 單因素評價矩陣單因素評價矩陣Ri是分別以水量、水質、人口、社會經濟、生態環境各評價指標的評價向量。 Ri定義如下：

Ri=ωо μ （1）

式中： ω為單因素中各評價因子的權重；μ 為通過各自的隸屬度函數確定的各評價因子的隸屬度矩陣。

（4） 模糊綜合評價

V=Aо R （2）

式中： A =（a1, a2,…， an）為各因素對水資源價值的指標權重，且滿足。R 為以水量、水質、人口、社會經濟、生態環境各評價指標評價向量 Ri為基礎構造的綜合矩陣 .

2. 2 相對隸屬度的確定

選用降半梯形分布，建立一元線性隸屬函數：

當j=1時，



當j=2、3、4時，



當j=5時，



式中： x 為評價因子的實際值； xi,j- 1、xij為評價因子相鄰兩等級的設定標準值； i 為評價因子標號；j=2,3…n,n=5; μn（x）為評價因子 i 的隸屬度。

2. 3 確定權重

2. 3. 1 層次分析法步驟

（1）對 每 一 層 次 各 因 素 的 相 對 重 要 性 1,2,…，9 及它們的倒數作為指標間相對重要性標度，構建指標兩兩對比判斷矩陣；（2） 計算判斷矩陣的特征根和特征向量，得到指標層對于要素層的重要性權值進行層次單排序；（3） 對判斷矩陣進行一致性檢驗；（4） 進行層次總排序，得出各個評價指標對目標層的重要性權值。

2. 3. 2 熵值法應用步驟[14-15]

（1） 對 R 矩陣中的 Xij進行級差標準化處理得到新的矩陣R':

正向指標：



式中： xij為指標初始值；xij‘為標準化處理后的值；max（xij）、min（xij）為相應指標初始最大值、最小值。