水資源短缺風險評價

頁關系矩陣可表示為式中：表示因素的評價對等級的隸屬度,因而矩陣中第行為對第個因素的單因素評判結果.在評價計算中代表了各個因素對綜合評判重要性的權系數(shù),因此滿足；同時,模糊變換也即退化為普通矩陣計算,即取合成運算,即用模型計算,可得綜合評判.通過模糊綜合評判模型,我們又對水資源短缺風險進行了分析,建立起相應的風險指標體系.該指標體系分為3個層次,共由16個指標組成.如表6所示,上述權系數(shù)的確定可用層次分析法(AHP)得到.由上述分析可以看出,評價因素集對應評語集,而評判矩陣中即為某因素對應等級的隸屬度,其值可根據(jù)各評價因素的實際數(shù)值對照各因素的分級指標推求.六、模型的求解北京市水資源短缺風險的模糊綜合評判模型求解：（1）因素集,其中指短缺性,指危險性,指易損性,指承險性.（2）評判集,其中：低；：較低；：中；：較高；：高.（3）單因素評判.依據(jù)我們的調查問卷的數(shù)據(jù),利用層次分析法,我們計算出了短缺性對水資源短缺風險的影響程度=(0.150.150.30.350.15),危險性對水資源短缺風險的影響程度=(0.150.30.30.150.1),易損性對水資源短缺風險的影響程度=(0.20.350.350.050.05)承險性對水資源短缺風險的影響程度=(0.20.40.20.150.05),便得到→(0.150.150.30.350.15)→(0.150.30.30.150.1)→(0.20.350.350.050.05)→(0.20.40.20.150.05)即得到1個到得模糊映射由此單因素評判可誘導出模糊關系,即得單因素評判矩陣（4）綜合評判.同樣利用層次分析法的到短缺性、危險性、易損性、承險性關于水資源短缺風險的權重分配.如下圖2：圖2取合成用算,即用模型Ⅰ：（主因素決定型）,計算可求得綜合評判為這表明水資源短缺危險程度較高,需要政府相關部門及全人類的高度注視.下面再用模糊綜合評價的另1種方法即最大隸屬原則,對北京市水資源短缺風險進行評價.我們將評語級分為5個級別,各評價因素分級指標見下表：表7：水資源短缺風險評價分級指標水資源短缺風險(短缺性）(危險性）（易損性）(承險性）(低）≤0.200≤0.200≤0.200≥0.800（較低）0.200～0.4000.200～0.4000.200～0.4000.601～0.800(中）0.401～0.6000.401～0.6000.401～0.6000.401～0.600(較高）0.601～0.8000.601～0.8000.601～0.8000.200～0.400(高）≥0.800≥0.800≥0.800≤0.200我們在整理、分析調查問卷中用1表示水資源短缺低風險,2表示水資源短缺較低風險,3表示水資源短缺中風險,4表示水資源短缺較高風險,5表示水資源短缺高風險.通過求每個風險因子的風險等級的平均值,就得到短缺性、危險性、易損性、承險性的等級劃分圖如下表：表8：水資源短缺風險指標等級劃分類別低風險較低風險中風險較高風險高風險短缺性0.60.244.62危險性1.83.63.61.81.2易損性23.53.50.50.5承險性1.63.21.61.30.4依據(jù)上表可構造短缺性、危險性、易損性、承險性隸屬函數(shù)分別為：將0.6帶入,于是類似地,可算出其他指標的隸屬度,得到單因素評價矩陣為用（主因素突出型）計算可求得綜合評判為對進行歸1化,得按最大隸屬度原則北京市水資源風險處于較高風險等級,可見水資源供需狀況極度危險.七、未來五年水資源狀況的預測與分析根據(jù)對《北京統(tǒng)計年鑒》[4]中有關水資源情況的分析,本文采取趨勢預測法:基于歷史統(tǒng)計數(shù)據(jù)的分析,選取1定長度的、具有可靠性、1致性和代表性的統(tǒng)計數(shù)據(jù)作為樣本,進行回歸分析,并以相關性顯著的回歸方程進行趨勢外延.為了使數(shù)據(jù)更加精確,采取了Excel軟件進行數(shù)據(jù)的描點作圖（圖3—圖7）,從圖表中可以看出北京市水資源在各個階段的總體變化趨勢,為了對未來的水資源數(shù)據(jù)進行預測,又用MATLAB軟件進行了數(shù)據(jù)擬合,得到了擬合曲線的函數(shù)表達式.圖3：總用水變化趨勢圖4：農(nóng)業(yè)用水變化趨勢圖5：工業(yè)用水變化趨勢圖6：第三產(chǎn)業(yè)及生活等其他用水變化趨勢圖7：水資源總量變化趨勢以上是用Excel軟件對1979—2008年從總用水量、農(nóng)業(yè)用水量、工業(yè)用水量、第三產(chǎn)業(yè)及生活等其他用水量和水資源總量來描點,對它們總體的變化趨勢進行分析,進而用MATLAB軟件對它們未來五年的水資源情況進行預測.在MATLAB程序中,為了使擬合函數(shù)的表達式的誤差最小,避免大數(shù)運算帶來的截斷誤差,我們用1—30分別代表1979—2008（年）,所用的程序如下（以工業(yè)用水數(shù)據(jù)為例）：Formatlong;x=1:1:30y=[38.23,26,24,36.6,34.7,39.31,38,27.03,38.66,39.18,21.55,35.86,42.29,22.44,19.67,45.42,30.34,45.87,22.25,37.7,14.22,16.86,19.2,16.1,18.4,21.4,23.2,24.5,23.8,34.2]plot(x,y,'k.','markersize',25)p4=polyfit(x,y,4)t=1:1:30s=polyval(p4,t)holdonplot(t,s,'r-','linewidth',2)plot(t,s,'b--','linewidth',2)grid;a=polyfit(t,y,4)總用水量趨勢預測所得四次多項式擬合曲線的函數(shù)表達式（這里用x表示年份,y表示水量,下同）為：=0.00044588652621-0.0304696800778+0.67862309754160-5.67157995768602+55.71015915119385(1)農(nóng)業(yè)用水量趨勢預測所得四次多項式擬合曲線的函數(shù)表達式為：=0.00017119427989-0.01317566138966+0.33904269110566-3.71699485945023+34.97676294331482（2）（3）工業(yè)用水趨勢預測所得四次多項式擬合曲線的函數(shù)表達式為：=0.00013215070897-0.00832866640558+0.15155240236764-0.93093484246700+14.92019746536991(3)（4）第三產(chǎn)業(yè)及生活等其他用水趨勢預測所得四次多項式擬合曲線的函數(shù)表達式為：=0.00014035806804-0.00881793151789+0.18477810891512-0.99897078802151+5.78234428725811（4）（5）水資源總量趨勢預測所得三次多項式擬合曲線的函數(shù)表達式為：=0.00567170437536-0.27452014650985+3.19607030989301+25.33436708629838（5）由以上各函數(shù)表達式,將未來五年的年數(shù)對應在函數(shù)中,由于在前面我們用1—30來代表1979—2008年來減少誤差,因此,這里用31—37來代表2009—2015年,將其帶入函數(shù)中,可以算得到未來五年內各水資源的數(shù)據(jù),如下表所示：表9未來五年北京市水資源狀況預測單位：億立方米年份總用水量農(nóng)業(yè)用水工業(yè)用水第三產(chǎn)業(yè)及生活等其他用水水資源總量200936.111311.1555.627719.31463829.564430201038.245010.9825.976221.258184132.350395201141.365911.0626.652923.603235.676287201245.6386511.4997.710826.4063939.57613201351.2390912.20119.206329.72756244.08397201458.353513.381211.198633.6300649.2338201567.19715.056213.750638.1805855.05972通過對以上數(shù)據(jù)的分析可以得到北京的用水量、用水結構、水資源存量的相關信息.7.1.用水量變化分析1980年到1990年,北京市用水總量呈明顯下降趨勢,年均減少總用水量0.087億立方米.進入1990年以來,年用水總量間的變幅則急劇縮小,介于穩(wěn)定的40.01億立方米到46.43立方米之間,今后五年內將穩(wěn)步上升,具體數(shù)據(jù)如表8所示7.2.用水結構變化分析北京市用水結構及其變化大體可按工、農(nóng)業(yè)和第三產(chǎn)業(yè)及生活等其他用水等3個方面進行分析.（1）農(nóng)業(yè)用水比重縮小,呈繼續(xù)縮減態(tài)勢自1980年以來,農(nóng)業(yè)用水作為北京市的用水大戶,其用水量的減少趨勢最為明顯,由1980年的31.83億立方米降至2000年的16.49億立方米,1980年—1990年、1990年—2000年和1996年—2000年年均減少量分別為0.247億立方米、0.477億立方米和0.638億立方米.農(nóng)業(yè)用水占全市總用水量的比重也呈下降趨勢由1980年的58.13%降為2000年的40.82%,近五年則平均以0.946%的份額下降,其下降趨勢仍無停止跡象.（2）工業(yè)用水呈減少趨勢,近年趨于穩(wěn)定工業(yè)用水亦呈負增長態(tài)勢,1980年用水量和占總用水量比重中分別為13.5億立方米及32.08%,到2000年下降為10.52億立方米和26.04%,年均遞減0.142億立方米,但近年來這種下降趨勢已明顯減緩.1997年—2000年工業(yè)用水總量介于10.5億立方米—11.0億立方米之間,變幅為4.5%.未來五年用水量將逐步上升,但幅度會越來越小.（3）第三產(chǎn)業(yè)及生活等其他用水持續(xù)增加,與工、農(nóng)業(yè)用水情況相反,城市及生活用水量從1980年的4.94億立方米迅速遞增為2000年的13.39億立方米,所占比重從9.79%增長到33.14%,而且不同階段的年均增加量呈逐步上升趨勢.1980年—1990年、1990年—2000年和1996年—2000年城市及生活用水量年均增加量分別為0.265億立方米、0.577億立方米和0.5億立方米,相應占總用水量比重的年均增加值分別為0.666%、1.456%和1.588%.未來五年持續(xù)增加的狀況不變.預計未來北京市用水結構總體趨勢為:總用水量不會發(fā)生大的變化,工業(yè)用水基本保持穩(wěn)定或略有增加,農(nóng)業(yè)用水量和占總用水量的比重仍將呈下降趨勢,生活用水量與比重將持續(xù)遞增.7.3．水資源存量變化分析水資源總量先呈減少趨勢然后逐步上升,水資源總量從1980年和1990年分別為26億立方米和35.86億立方米,呈上升狀況,1990年—2000水資源總量從35.86億立方米減少到16.86億立方米,進入2000年以后,由數(shù)據(jù)可知,水資源總量又穩(wěn)步上升,在今后五年里,如果沒有其他因素的干擾,水資源總量會逐步上升,但上升的幅度會越來越小.八、模型的評價與改進雖然,算子有很好的代數(shù)性質,但也存在著缺陷,它常常出現(xiàn)綜合評判的結果不易分辨的情況,因此,模型Ⅰ需要改進,下面介紹改進數(shù)學模型的方法,即將原模型中的算子改用其他算子.模型Ⅰ：綜合評判的著眼點是考慮主要因素,其他因素對結果影響不大,為了避免出現(xiàn)決策結果不易分辨的情況,以下對模糊綜合評判決策模型進行改進.模型Ⅱ：（主因素突出型）計算可求得綜合評判為對進行歸1化,得通過觀察,表明水資源短缺危險程度較高,這與模型Ⅰ：的結果是1致的,即表明水資源短缺風險突出,其水資源開采、利用、再生等治理迫在眉睫.但為了避免權重與主要因素有關而忽略次要因素,我們還可以對其進行改進.模型Ⅲ：（主因素突出型）這里的為有界和,即計算可求得綜合評判為我們能明顯的看出水資源短缺危險程度較高,這與模型Ⅰ：和模型Ⅱ：的結果是1致的,即表明水資源短缺風險突出,是不能被忽略的.模型Ⅲ：在實際應用中,主因素（權重最大的因素）在綜合評價中起決定作用,為了避免其帶來的負面影響,我們用另外1種模型即模型Ⅳ進行檢驗.模型Ⅳ：（加權平均模型）計算可求得綜合評判為模型Ⅳ對所有因素以權重大小均衡兼顧,適用于考慮各因素起作用的情況.從模型Ⅳ可以檢驗出前三種模型在數(shù)據(jù)分析及用算過程中是適合的,并沒有帶來較大的偏差.九、結論（1）本文基于模糊綜合評判模型建立了水資源短缺風險評價模型,同時考慮到水資源系統(tǒng)的模糊不確定性,可對水資源短缺風險的影響程度給予綜合評價.社會經(jīng)濟因素是北京市水資源短缺的主要致險因子.（2）由模糊綜合評價模型可以得出北京市水資源短缺風險處于較高水平,根據(jù)我們建立的北京市水資源系統(tǒng)風險評價指標體系及評價模型,對北京市進行水資源系統(tǒng)風險評價,得出北京市水資源短缺風險級別為較高風險,從而為北京市水資源系統(tǒng)管理管控及水資源系統(tǒng)風險控制提供依據(jù).對水資源采取有效的風險管理管控措施已刻不容緩.（3）北京市用水結構變化總的來說呈現(xiàn)以下趨勢:總用水量趨于平穩(wěn),工業(yè)用水和農(nóng)業(yè)用水從量上和占總用水量的比重上都有所下降,而生活用水卻迅速遞增.下面是給北京市水資源管理管控相關部門的研究報告.北京市水資源短缺研究報告北京市水行政主管部門：由于北京市水資源短缺已經(jīng)成為影響和制約首都社會和經(jīng)濟發(fā)展的主要因素.我們對北京市水資源資料的分析,對水資源短缺的風險因子進行了重要性分析.概括來說,造成水資源短缺的成因主要源于以下四方面:1.環(huán)境因素;2.工業(yè)因素;3.農(nóng)業(yè)因素,4.社會經(jīng)濟因素.在這些因素中,社會經(jīng)濟因素對水資源造成短缺的作用越來越大,隨著人口的增長,居民生活水平的提高帶來的居民生活用水的迅速增長,城市建設、環(huán)境質量的提高以及服務業(yè)的蓬勃發(fā)展造成了公共用水的增加,共同推動了北京市生活用水迅速增長.經(jīng)過我們對北京市水資源狀況分析,北京市用水結構變化總的來說呈現(xiàn)以下趨勢:總用水量趨于平穩(wěn),工業(yè)用水和農(nóng)業(yè)用水從量上和占總用水量的比重上都有所下降,而生活用水卻迅速遞增;預計今后境內自產(chǎn)水量變化不大.地表水資源可隨調蓄和聯(lián)合調度能力增加而有所增加,但地下水資源已多年超采,能夠保持現(xiàn)有資源量不使地下水位繼續(xù)下降已是最好的結果.密云、官廳兩大水庫目前和今后很長1段時間內都將是北京主要的水源,但來水減少之勢不可逆轉.現(xiàn)在北京的出境水大部分是污水,深度開發(fā)利用這部分水資源,使之成為可能利用的新水源,是緩解北京市嚴重缺水的重要措施.因此,北京市未來的水資源狀況不容樂觀,在沒有外來水源引入的條件下,供水量不會發(fā)生明顯的變化.由以上原因北京市水行政主管部門應該從以下幾個方面做出改善,來保護水資源以及降低水資源的短缺風險.（1）降低水資源短缺風險1、對北京市的水資源短缺風險必須進行調控,這些調控措施主要有需水管理管控和供水管理管控.需水管理管控的核心是抑制水資源需求的過度膨脹,促進水資源的可持續(xù)利用,節(jié)水防污型社會建設是需水管理管控中最重要的系統(tǒng)工程之1；供水管理管控措施主要有提供污水處理率和污水利用率、對當?shù)厮Y源進行挖潛、增加雨洪利用、增加海水利用等等.2、南水北調工程是解決北京市水資源短缺風險的根本措施.如果將南水北調工程加入到風險分析模型中,2010年北京市的水資源短缺風險率、脆弱性等將得到大幅度的改善,所以當前最重要的工作是保證南水北調東、中線工程按規(guī)劃準時建成,保證首都圈的經(jīng)濟社會可持續(xù)發(fā)展.密云、官廳兩大水庫目前和今后很長1段時間內都將是北京主要的水源,但來水減少之勢不可逆轉.現(xiàn)在北京的出境水大部分是污水和雨洪水,深度開發(fā)利用這部分水資源,使之成為可能利用的新水源,是緩解北京市嚴重缺水的重要措施.(2)提高污水資源化程度.污水是被污染、使用價值不高的水資源,污水資源化是指污水處理后變成可供用水部門使用的合格水源.實現(xiàn)污水資源化,提高水資源重復利用率,是緩解北京市水資源緊缺、防治水污染的重要途徑.目前,北京市大部分污水未經(jīng)處理直接外排,不僅浪費水資源,而且造成環(huán)境污染.因此,建議北京市北郊、西郊污水處理廠盡快上馬.污水處理廠的建設,是實現(xiàn)污水資源化的核心,處理過的污水,除農(nóng)灌外,還可回用于部分工業(yè)企業(yè)、市政及城市景觀用水.建筑中水回用是污水資源化的另1途徑.建筑中水回用指建筑物的各種排水經(jīng)處理后回用于建筑物和建筑小區(qū)雜用的供水系統(tǒng),它可以作為沖廁、綠化、洗車、澆灑道路、做冷卻水等用.(3)發(fā)展節(jié)水型產(chǎn)業(yè).根據(jù)北京市的水資源狀況和環(huán)保要求,21世紀北京市的產(chǎn)業(yè)結構應以節(jié)水型、清潔型產(chǎn)業(yè)為主.推廣高新技術,生產(chǎn)和使用節(jié)水材料和節(jié)水設備,采用節(jié)水工藝,努力降低萬元產(chǎn)值耗水量.做到科學用水,合理配置水資源,這是建立資源節(jié)約型城市發(fā)展體系的必然.三、為了實現(xiàn)節(jié)水目標,必須積極采取各類有效的節(jié)水措施,既包括工程措施,也包括行之有效的非工程措施.對于工、農(nóng)業(yè)和生活用水的節(jié)水措施如下:1、農(nóng)業(yè)節(jié)水措施（1）改造配套灌區(qū)工程:農(nóng)業(yè)節(jié)水重點是結合大中型灌區(qū)改造配套工程規(guī)劃,續(xù)建、配套和技術改造現(xiàn)有灌區(qū)節(jié)水工程,減少輸水過程中的蒸發(fā)滲漏損失,提高渠系水的利用效率.改進作物灌溉制度:傳統(tǒng)的灌溉制度是建立在充分灌溉的基礎上,每次灌水都滿足“及時足量”的灌溉要求,造成大量蒸發(fā)和滲漏損失,因此必須根據(jù)農(nóng)作物灌溉的實際需求,改變灌溉制度和灌溉方式,采用經(jīng)濟灌溉定額制度和實施有限灌溉,減少田間灌溉水的無效損失量.（2）加強用水管理管控,科學調配水量:由于灌區(qū)大多數(shù)渠系龐大、用戶眾多,管理管控的水平對用戶效率起很大的作用,因此加強用水管理管控,實行科學調配水量,維護工程良好狀態(tài),將會取得明顯的節(jié)水效果.全面實施灌溉用水定額管理管控制度,實行用水總量控制和定額管理管控,可有效減少灌溉水分的損失,提高用水效率.（3）強化農(nóng)業(yè)用水定量和價格管理管控:①使用地表水的渠灌區(qū),強化定量收費,累進加價制度;②地下水井灌區(qū)實行計量收費,配額管理管控的辦法.③實行合理的水價政策.要調整水費政策,以成本定價水價,按引水量計征水費以及對超用水加價收費,從而促進節(jié)水.2、工業(yè)節(jié)水措施工業(yè)節(jié)水主要為產(chǎn)業(yè)結構的調整、用水裝置的更新?lián)Q代、生產(chǎn)工藝的改進、節(jié)水器具的推廣、管理管控水平的提高等方面.工業(yè)節(jié)水要抓四個方面:①全面調整工業(yè)結構,有相關計劃、有重點的限制和壓縮高耗水工業(yè)相關項目的建設;②通過各種行政手段加強用水管理管控,相關計劃用水和嚴格控制廢污水的排放,可獲得明顯的節(jié)水效果.3)生活節(jié)水措施目前城市供水管網(wǎng)和用水器具水的漏失率較高,要采取有效措施減少跑、冒、滴、漏現(xiàn)象,要進1步普及節(jié)水型器具.總之,要采取1切節(jié)水措施,挖掘潛力,提高水的利用率,使有限的水資源適應人民生活水平日益提高的要求.(4)建立合理的水資源價格體系.合理的水資源價格是實現(xiàn)水資源可持續(xù)利用的關鍵.美國研究表明:當生活用水實行按水量計費時,節(jié)約用水10%～30%,城市用水峰值系數(shù)下降30%～50%;當水價上漲5倍時,工業(yè)萬元耗水量可減少1/2[11,12].由此可見,合理的水資源價格體系對緩解水資源供需矛盾非常有益.21世紀北京市的水資源管理管控應充分利用這1經(jīng)濟杠桿的作用,合理調節(jié)不同類型水資源的價格,制定出切合實際的工業(yè)用水、生活用水及其它用水的水費標準.(5)加強宣傳教育,提高市民的節(jié)水意識.水是1種有價值的有限資源,要解決水資源浪費問題,必須加強節(jié)水宣傳,開展科普教育,提高市民用水的道德觀,形成社會節(jié)水的共識.節(jié)水意識是多因素相互作用的共同結果,它包括多方面的內涵,如經(jīng)濟的、政治的、道德倫理及法律方面等,但經(jīng)濟是主要的,法律是保證.水價過低,會給人們以錯誤的導向,不利于節(jié)水.因此,提高水價是提高市民節(jié)水意識的關鍵,并促使市民盡快建立起水資源價值觀.參考文獻[1]姜啟源.數(shù)學模型[M].高等教育出版社,1993[2]徐全智,楊晉浩.數(shù)學建模入門[M].成都:電子科技大學出版社,1996[3]陳汝棟,于延榮.數(shù)學模型與數(shù)學建模[M].國防工業(yè)出版社,2009[4]衷平,沈珍瑤,楊志峰,等.石羊河流域水資源短缺風險敏感因子的確定[J].干旱區(qū)資源與環(huán)境,2005,19(2):81-86.[4]北京市統(tǒng)計局.北京統(tǒng)計年鑒[M].北京:中國統(tǒng)計出版社,2008-2009[5]北京市人民政府,水利部.21世紀初期首都水資源可持續(xù)利用規(guī)劃[M].北京:中國水利水電出版社,2001.[6]張井躍.模糊數(shù)學方法及其應用[M].北京:煤炭工業(yè)出版社,1992,55-100.[7]劉濤,邵東國.水資源系統(tǒng)風險評估方法研究[J].武漢大學學報(工學版),2005,38(6):66-71.[8]劉濤,邵東國,顧文權.基于層次分析法的供水風險綜合評價模型[J].武漢大學學報(工學版),2006,39(4):25-28.[9]WANGHong-rui,QIANLong-xia,XUXin-yi,etal.Modelforevaluatingwatershortageriskbasedonfuzzyprobabilityanditsapplication.JournalofHydraulicEngineering,2009,40(7):813-821.[10]HydrologicCentralStationofBeijing.BeijingWaterResourcesBulletin.2008-2009.[11]HANYu-ping.ManagementTheoryandPracticeonWaterShortageRisk.Zhengzhou:YellowRiverConservancyPress,2008:36-40.附錄1：為了計算的準確性和有效性,我們借助VB程序來使我們得到的結果更加具有可靠性,下面是1段程序,以供參考.1．算權重Dima(10,10)AsDouble,b(10)AsDoublePrivateSubForm_Load()ShowDimmAsInteger,nAsInteger,iAsInteger,jAsInteger,piAsDouble,sumAsDoublem=InputBox("輸入階數(shù)")i=1:j=1Fori=1TomForj=1Toma(i,j)=InputBox("輸入矩陣")NextjNextii=1:j=1:pi=1Fori=1TomForj=1Tompi=pi*a(i,j)Nextjb(i)=pi:pi=1Nextii=1Fori=1Tomb(i)=b(i)^(1/m)NextiFori=1TomPrintb(i)NextiEndSub2．算評價指標值Dima(10)AsDouble,b(10)AsDouble,c(10)AsDouble,d(10)AsDoublePrivateSubForm_Load()ShowDimiAsInteger,nAsInteger,sumAsDoublen=InputBox("輸入階數(shù)")sum=0Fori=1Tona(i)=InputBox("輸入第1列數(shù)")NextiFori=1Tonb(i)=InputBox("輸入第二列數(shù)")NextiFori=1Tonc(i)=b(i)*a(i)sum=sum+c(i)NextiFori=1TonPrintc(i)NextiPrintsumEndSub附錄二調查問卷親愛的同學們：你們好,本調查問卷旨在對北京市水資源短缺風險進行綜合評價