新疆博州地區(qū)水資源時(shí)間分布規(guī)律與綜合評(píng)價(jià)研究

烏仁娜

(新疆博爾塔拉水文勘測(cè)局，新疆 博爾塔拉 834000)

水資源安全問題已經(jīng)成為了國(guó)內(nèi)相關(guān)學(xué)科研究的熱點(diǎn)[1- 2]。是能否持續(xù)保障國(guó)家可持續(xù)發(fā)展的前提，水資源安全水平是一個(gè)地區(qū)水資源優(yōu)化配置前提下，城市水資源整體水平的體現(xiàn)，對(duì)其分布規(guī)律和整體水平的評(píng)價(jià)是實(shí)施水資源安全發(fā)展戰(zhàn)略的前提[4]。研究水資源安全問題最主要的就是研究水資源的分布規(guī)律與水資源綜合評(píng)價(jià)的度量標(biāo)準(zhǔn)[3]。

目前針對(duì)水資源分布規(guī)律，已有了部分研究。李文清[5]研究了沂州市河川徑流分布規(guī)律，指出掌握各流域河川徑流量在時(shí)空上的分布規(guī)律，可為緩解地區(qū)水資源壓力提供科學(xué)依據(jù)；劉斌濤等[6]研究了云南省水資源時(shí)空分布的格局，指出云南省水資源分布極不均衡，且春季水資源較少；呂立新[7]研究了吉林省地表水資源的分布規(guī)律，指出年徑流量與下墊面的不同有關(guān)。

目前針對(duì)水資源安全評(píng)價(jià)，同樣有了部分研究。趙西寧等[8]在黃土高原評(píng)價(jià)了雨水資源化綜合效益，基于AHP法建立了水資源安全指標(biāo)綜合評(píng)價(jià)體系；趙洪杰和唐善德[9]、周維博和李佩成[10]分別在黑河中游灌區(qū)和干旱半干旱地區(qū)灌區(qū)，利用層次分析法，建立了水資源效益綜合評(píng)價(jià)指標(biāo)體系，并對(duì)灌區(qū)水資源綜合效益進(jìn)行了評(píng)價(jià)，結(jié)果與實(shí)際情況較為一致。

新疆位于我國(guó)西北部，氣候干燥，降水量少，多年平均降水量?jī)H為150mm左右，其水資源含量較低，目前對(duì)新疆地區(qū)水資源分布規(guī)律與綜合評(píng)價(jià)的研究較少。本文以新疆博州地區(qū)為研究區(qū)域，基于水文站實(shí)測(cè)資料，研究當(dāng)?shù)厮Y源年際及年內(nèi)變化規(guī)律，并依據(jù)TOPSIS模型構(gòu)建水資源綜合評(píng)價(jià)體系。

1 研究方法與研究區(qū)域概況

1.1 研究區(qū)域概況

本文選擇新疆博州地區(qū)精河山口站、溫泉站和沙爾托海站為代表站點(diǎn)，基于當(dāng)?shù)?988—2016年近29年的實(shí)測(cè)數(shù)據(jù)，對(duì)當(dāng)?shù)厮Y源時(shí)間分布規(guī)律進(jìn)行分析，站點(diǎn)經(jīng)度最大相差1°53′，最大落差1072.6m，具體研究區(qū)域及站點(diǎn)分布圖，如圖1所示及文獻(xiàn)[11]。

圖1 研究區(qū)域位置圖

圖2 水資源綜合效益評(píng)價(jià)指標(biāo)體系構(gòu)建

1.2 研究方法

1.2.1 評(píng)價(jià)指標(biāo)體系選擇

本文從工程效益、管理效益和社會(huì)效益3個(gè)方面對(duì)水資源效益進(jìn)行評(píng)價(jià)?；贕1法和熵權(quán)法計(jì)算指標(biāo)權(quán)重，依據(jù)TOPSIS模型建立新疆博州地區(qū)水資源效益綜合評(píng)價(jià)模型，綜合判定區(qū)域水資源綜合效益。其指標(biāo)體系構(gòu)建如圖2所示。

1.2.2 綜合權(quán)重確定

為得出最能反映每個(gè)指標(biāo)體系的權(quán)重，選擇主觀權(quán)重賦值法的G1法和客觀權(quán)重賦值法的熵權(quán)法對(duì)不同指標(biāo)進(jìn)行權(quán)重計(jì)算，最終得出綜合權(quán)重，以下為綜合權(quán)重的計(jì)算步驟。

(1)G1法權(quán)重賦值

假設(shè)評(píng)價(jià)指標(biāo)體系中有n個(gè)評(píng)價(jià)指標(biāo)α1，α2，…，αn，由專家經(jīng)驗(yàn)判斷不同指標(biāo)的重要性，設(shè)置不同指標(biāo)的重要性程度為γn，則不同權(quán)重賦值為：

(1)

ωk+1=ωk×γkk=1，2，3……n

(2)

式中，ωk、ωk+1—第k、k+1個(gè)指標(biāo)的權(quán)重；γk—該指標(biāo)的重要性程度。

(2)熵權(quán)法權(quán)重賦值

假設(shè)有m個(gè)對(duì)象，n個(gè)被評(píng)價(jià)指標(biāo)，熵定義為：

(3)

(4)

式中，Hj—計(jì)算的信息熵；m—評(píng)價(jià)對(duì)象個(gè)數(shù)；n—待評(píng)價(jià)指標(biāo)個(gè)數(shù)；bij—計(jì)算的相對(duì)隸屬度；fij—相對(duì)隸屬度所占比例。

對(duì)各項(xiàng)指標(biāo)進(jìn)行無(wú)量綱化處理，得出相對(duì)隸屬度：

(5)

式中，bij—所求的相對(duì)隸屬度；aij—每個(gè)計(jì)算指標(biāo)；ajmax—某項(xiàng)指標(biāo)最大值；ajmin—某項(xiàng)指標(biāo)最小值。

計(jì)算第j項(xiàng)指標(biāo)的變異度：

Dj=1-Hj

(6)

式中，Dj—所求的變異度；Hj—指標(biāo)的信息熵。

計(jì)算該指標(biāo)所對(duì)應(yīng)的熵權(quán)重(客觀權(quán)重)：

(7)

式中，ωbj—所求的客觀權(quán)重；Dj—某項(xiàng)指標(biāo)的變異度。

(3)綜合權(quán)重賦值

根據(jù)專家經(jīng)驗(yàn)得出的主觀權(quán)重與計(jì)算得出的客觀權(quán)重，得出綜合權(quán)重：

(8)

式中，ωj—所求的綜合權(quán)重；ωaj—通過G1法所得的主觀權(quán)重；ωbj—熵權(quán)法所得的客觀權(quán)重。

1.2.3 TOPSIS水資源效益綜合評(píng)價(jià)模型構(gòu)建

(1)形成無(wú)量綱決策矩陣

假設(shè)目標(biāo)體系的方案層矩陣為N=(N1，N2……，Nm)，指標(biāo)層矩陣為A=(A1，A2……，An)，矩陣Z表示N對(duì)于C的比值，無(wú)量綱化矩陣用M表示可用下式表示。

(9)

(10)

(2)計(jì)算最優(yōu)解，找出各方案與最優(yōu)解的相對(duì)貼近程度

將無(wú)量綱化矩陣M，與計(jì)算得出的權(quán)重相乘，得出加權(quán)決策矩陣R：

R=ωj×mij

(11)

找出不同指標(biāo)的最優(yōu)解，采用文獻(xiàn)[12]的方法找出最優(yōu)解。通過下式計(jì)算不同指標(biāo)與最優(yōu)解的正負(fù)理想解距離，通過公式(14)計(jì)算相對(duì)貼近度矩陣。

(12)

(13)

(14)

根據(jù)ε值的大小，對(duì)不同方案進(jìn)行排序，其值越大，則方案更優(yōu)。

2 結(jié)果與分析

2.1 年際變化趨勢(shì)分析

博州地區(qū)3個(gè)水文站點(diǎn)實(shí)測(cè)徑流值的年際變化及累積距平值如圖3所示。圖3顯示，3個(gè)水文站1988—2016年多年徑流量基本變化較平穩(wěn)，呈緩慢上升的趨勢(shì)，精河山口站、溫泉站和沙爾托海站徑流量多年變化趨勢(shì)分別為0.391、0.287、

0.209mm/10a，上升趨勢(shì)均不明顯。3個(gè)站點(diǎn)在1988—1997年呈下降趨勢(shì)，在1997—2007年呈上升趨勢(shì)，在2007—2016年呈降低趨勢(shì)，3個(gè)站點(diǎn)徑流變化趨勢(shì)基本一致。

2.2 年內(nèi)變化趨勢(shì)分析

博州地區(qū)3個(gè)水文站點(diǎn)實(shí)測(cè)徑流值的年內(nèi)變化及累積距平值如圖4所示。圖4顯示，3個(gè)水文站1—12月年內(nèi)基本呈先升高后降低的二次拋物線形式，3個(gè)站點(diǎn)徑流值的年內(nèi)分布基本在4—10月變化較劇烈，在其余月份基本不變，同時(shí)4—10月的徑流值較大，占了全年徑流了的65.4%～87.4%。各站點(diǎn)累積距平曲線變化趨勢(shì)基本一致，呈“減少-增加-減少”趨勢(shì)，3個(gè)時(shí)段分別為1—5月、5—9月和9—12月。

2.3 水資源評(píng)價(jià)結(jié)果分析

根據(jù)標(biāo)準(zhǔn)，將水資源評(píng)價(jià)指標(biāo)體系進(jìn)行分級(jí)，采用標(biāo)準(zhǔn)差分級(jí)法進(jìn)行，該方法適用于呈正態(tài)分布的數(shù)據(jù)集，具體指標(biāo)分級(jí)情況見表1。

通過調(diào)查3個(gè)水文站附近區(qū)域的各項(xiàng)資料，得出多年平均資料用于進(jìn)行3個(gè)站的水資源效益評(píng)價(jià)，通過式(1)～(8)計(jì)算不同指標(biāo)的客觀權(quán)重和主觀權(quán)重，最終得到不同指標(biāo)的最終權(quán)重，結(jié)果見表2。通過表1中的數(shù)據(jù)，結(jié)合表2中的權(quán)重和式(11)～(14)，可得出水資源效益具體等級(jí)評(píng)價(jià)標(biāo)準(zhǔn)，結(jié)果見表3。通過實(shí)測(cè)數(shù)據(jù)，可得出不同站點(diǎn)多年評(píng)價(jià)的水資源效益相對(duì)貼近度值，結(jié)果見表4。表4顯示，3個(gè)站點(diǎn)水資源效益由高到低排名依次為精河山口站、溫泉站、沙爾托海站，精河山口站水資源評(píng)價(jià)等級(jí)為一般，其余2個(gè)站點(diǎn)為較差，該評(píng)價(jià)結(jié)果與當(dāng)?shù)貙?shí)際情況基本一致。

圖3 當(dāng)?shù)貜搅髂觌H變化及累積距平值

圖4 當(dāng)?shù)貜搅髂陜?nèi)變化及累積距平值

指標(biāo)好較好一般較差差ⅠⅡⅢⅣⅤ人均地區(qū)生產(chǎn)總值C11/元≥80006000～80004000～60002000～4000≤2000當(dāng)?shù)赜盟緾12/(106m3/億)≥7560～7545～6030～45≤30已建水利工程覆蓋率C13/%≥7555～7535～5515～35≤15工程總投資C14/萬(wàn)元≥15001200～1500900～1200600～900≤600新型作物推廣面積比C21/%≥5040～5030～4020～30≤20農(nóng)民年收入C22/元≥60005000～60004000～50003000～4000≤3000作物水分生產(chǎn)率C23≥1.801.60～1.801.40～1.601.20～1.40≤1.20政府支持程度C24≥8.07.0～8.06.0～7.05.0～6.0≤5.0計(jì)劃用水水平C31≥6.56.0～6.55.5～6.05.0～5.5≤5.0節(jié)水機(jī)制健全水平C32≥5.455.05～5.454.65～5.054.25～4.65≤4.25新型節(jié)水工程覆蓋率C33/%≥5540～5525～4010～25≤10技術(shù)推廣體系水平C34≥0.850.70～0.850.55～0.700.40～0.55≤0.40

表2 不同指標(biāo)權(quán)重確定

表3 不同等級(jí)劃分標(biāo)準(zhǔn)

表4 不同站點(diǎn)水資源綜合效益排名

3 結(jié)語(yǔ)

通過分析實(shí)測(cè)數(shù)據(jù)，掌握了新疆博州地區(qū)近年來(lái)的水資源時(shí)間分布規(guī)律，指出近年來(lái)當(dāng)?shù)厮Y源呈緩慢上升的趨勢(shì)，且年內(nèi)呈二次拋物線趨勢(shì)。水資源效益評(píng)價(jià)是一個(gè)多指標(biāo)、多屬性的過程，文章從工程效益、社會(huì)效益和管理效益3個(gè)層面建立的水資源效益綜合評(píng)價(jià)模型體系，可較好地反映當(dāng)?shù)厮Y源的發(fā)展情況，同時(shí)本文權(quán)重確定的方法，可從主觀和客觀多角度反映評(píng)價(jià)指標(biāo)的現(xiàn)狀，具有一定的現(xiàn)實(shí)意義。本文僅建立了12項(xiàng)評(píng)價(jià)指標(biāo)體系，在今后的研究中，可圍繞建立更多指標(biāo)體系角度進(jìn)行。