黃顯峰, 黃雪晴, 方國華, 王宗志, 朱麗向

(1.河海大學(xué) 水利水電學(xué)院,江蘇 南京 210098; 2.南京水利科學(xué)研究院 水文水資源研究所,江蘇 南京 210029; 3.連云港市通榆河北延送水工程管理處,江蘇 連云港 222004)

洪水資源利用風(fēng)險(xiǎn)效益量化研究

黃顯峰1, 黃雪晴1, 方國華1, 王宗志2, 朱麗向3

(1.河海大學(xué) 水利水電學(xué)院,江蘇 南京 210098; 2.南京水利科學(xué)研究院 水文水資源研究所,江蘇 南京 210029; 3.連云港市通榆河北延送水工程管理處,江蘇 連云港 222004)

結(jié)合洪水資源利用的風(fēng)險(xiǎn)因素,對(duì)水庫汛限水位調(diào)整后的洪水資源利用風(fēng)險(xiǎn)效益進(jìn)行量化研究。汛限水位是涉及水庫防洪風(fēng)險(xiǎn)和興利效益的關(guān)鍵水位。汛限水位調(diào)整后,可得到每個(gè)汛限水位對(duì)應(yīng)的風(fēng)險(xiǎn)率,并將水庫因水位調(diào)整而產(chǎn)生的增蓄水量按比例分配,綜合考慮工業(yè)、農(nóng)業(yè)、生活、生態(tài)4個(gè)方面的風(fēng)險(xiǎn)效益,采用C-D生產(chǎn)函數(shù)、分?jǐn)傁禂?shù)法、能值水價(jià)法和能值理論,分別對(duì)這4種效益進(jìn)行量化。以連云港市石梁河水庫為例,計(jì)算了各汛限水位對(duì)應(yīng)的風(fēng)險(xiǎn)率和效益值,并根據(jù)汛限水位、風(fēng)險(xiǎn)率、風(fēng)險(xiǎn)效益之間的關(guān)系,在考慮風(fēng)險(xiǎn)最小且效益最大的前提下,得到石梁河水庫汛限水位合理的控制范圍為23.9～24.9 m。

洪水利用;風(fēng)險(xiǎn);效益;C-D生產(chǎn)函數(shù);能值法

我國洪澇災(zāi)害頻發(fā),威脅著人民的生命財(cái)產(chǎn)安全。因此,越來越多學(xué)者開始了洪水資源利用研究,并在洪水資源可利用量和效益計(jì)算方面取得了一定成果[1-2]。但洪水利用存在風(fēng)險(xiǎn),如何控制風(fēng)險(xiǎn)、提高效益,在洪水資源化的過程中是一個(gè)尤為重要的課題。

國外學(xué)者比較側(cè)重于洪水風(fēng)險(xiǎn)的防范研究,對(duì)洪水利用效益提及較少。Roy Brouwer等[3]對(duì)荷蘭的洪水管理替代政策進(jìn)行了效益評(píng)估研究,采用CBA(Cost-Benefit Analysis)和MCA(Multi-Criteria Analysis)兩種方法對(duì)其產(chǎn)生的效益進(jìn)行了計(jì)算和對(duì)比。由于荷蘭更偏重于防洪減災(zāi),故此研究針對(duì)的是洪水管理政策產(chǎn)生的效益,其先進(jìn)的效益計(jì)算方法可作為參考。我國許多地區(qū)已開展了洪水資源化利用工作,在不同地區(qū)、不同條件下都有了相應(yīng)的利用方法。殷峻逞等[4]針對(duì)大連市碧流河水庫,建立了洪水資源效益分析模型,對(duì)該水庫洪水資源利用效益做出了定量分析。胥衛(wèi)平等[5]建立了雨洪資源利用的經(jīng)濟(jì)效益評(píng)價(jià)指標(biāo)體系及其經(jīng)濟(jì)價(jià)值評(píng)價(jià)模型,并將其應(yīng)用于西安市。馮峰等[6]構(gòu)建了社會(huì)效益分度測評(píng)模型,將洪水資源利用產(chǎn)生的社會(huì)效益進(jìn)行了定量評(píng)價(jià)。

洪水資源化的實(shí)施有多種途徑,目前，研究較多的是水庫汛限水位的調(diào)整。汛限水位是涉及水庫防洪風(fēng)險(xiǎn)和興利效益的一個(gè)關(guān)鍵水位,必須正確、全面地處理風(fēng)險(xiǎn)與效益的關(guān)系[7]。若只顧增大效益,勢必會(huì)減少防洪庫容,增大洪水風(fēng)險(xiǎn);若過度控制風(fēng)險(xiǎn),則預(yù)留的防洪庫容過大,造成不必要的棄水,降低興利效益。因此,本文在考慮汛限水位調(diào)整的基礎(chǔ)上,對(duì)風(fēng)險(xiǎn)率和效益值進(jìn)行計(jì)算,并提出工業(yè)、農(nóng)業(yè)、生活、生態(tài)4個(gè)方面效益量化的方法。

1 研究方法

1.1 洪水資源利用風(fēng)險(xiǎn)識(shí)別與估計(jì)

洪水資源作為一種非常規(guī)水資源,要在保證水利工程安全、生態(tài)環(huán)境不受影響的前提下，開展洪水資源利用。為滿足這些條件,必須進(jìn)行洪水資源利用風(fēng)險(xiǎn)識(shí)別。通過識(shí)別,洪水資源利用的風(fēng)險(xiǎn)因素主要有抬高水庫汛限水位的風(fēng)險(xiǎn)、水庫上游來水不確定性風(fēng)險(xiǎn)、不利生態(tài)環(huán)境影響風(fēng)險(xiǎn)、水質(zhì)不確定性風(fēng)險(xiǎn)、洪水調(diào)度管理風(fēng)險(xiǎn)等[8]。本文僅對(duì)汛限水位調(diào)整而導(dǎo)致的風(fēng)險(xiǎn)進(jìn)行估計(jì),主要計(jì)算的風(fēng)險(xiǎn)因子包括大壩安全風(fēng)險(xiǎn)率和下泄流量超下游河道安全泄量風(fēng)險(xiǎn)率。

運(yùn)用數(shù)理統(tǒng)計(jì)方法,以水庫上游來水的歷時(shí)統(tǒng)計(jì)資料為基礎(chǔ)，進(jìn)行調(diào)洪演算,找出汛限水位對(duì)應(yīng)的使水庫達(dá)到設(shè)計(jì)標(biāo)準(zhǔn)水位的洪水,以該洪水發(fā)生的概率作為汛限水位對(duì)應(yīng)的大壩安全風(fēng)險(xiǎn)率,計(jì)算公式如下:

(1)

式中:f1(Hi)為汛限水位為Hi時(shí)出現(xiàn)超設(shè)計(jì)標(biāo)準(zhǔn)水位的風(fēng)險(xiǎn)率;N為資料系列中歷次洪水的總數(shù)；n為給定的汛限水位Hi調(diào)洪演算發(fā)生的庫水位超設(shè)計(jì)標(biāo)準(zhǔn)水位的洪水次數(shù)。

以水庫下泄流量不超過河道安全泄量為準(zhǔn),用模糊數(shù)學(xué)法，計(jì)算下泄流量超過下游河道安全泄量的風(fēng)險(xiǎn)率,計(jì)算公式如下:

f2(Hi)=∑Vi∈AμA(Vi)·p(Vi) ×100%。

(2)

式中:f2(Hi)為汛限水位為Hi時(shí)下泄流量超下游河道安全泄量風(fēng)險(xiǎn)率;p(Vi)為汛限水位為Hi時(shí)下泄流量隸屬于區(qū)間Vi的概率;μA(Vi)為不同區(qū)間Vi的隸屬度。

水庫洪水資源利用的綜合風(fēng)險(xiǎn)率為:

f(H)=w1f1(Hi)+w2f2(Hi)。

(3)

式中:w1為大壩安全風(fēng)險(xiǎn)率的權(quán)重值;w2為下泄流量超下游河道安全泄量風(fēng)險(xiǎn)率的權(quán)重值。

1.2 洪水資源利用效益量化方法

洪水資源產(chǎn)生的效益即采取各種措施蓄存洪水,增供的水資源量產(chǎn)生的直接、間接效益,如促進(jìn)工、農(nóng)業(yè)生產(chǎn),提高人民生活水平,改善生態(tài)環(huán)境等而產(chǎn)生的效益。以調(diào)整后汛限水位為基礎(chǔ),庫容增加值作為洪水資源利用量,計(jì)算相應(yīng)的蓄水效益。主要研究利用洪水資源緩解水資源供需緊張的壓力,所以，在效益方面未考慮發(fā)電效益,只考慮工業(yè)供水、農(nóng)業(yè)灌溉、生活供水和生態(tài)供水4個(gè)方面效益。

1.2.1 C-D生產(chǎn)函數(shù)法量化工業(yè)供水效益

洪水資源利用產(chǎn)生的工業(yè)供水效益一般參考水資源效益來進(jìn)行量化,常用的量化方法包括殘值法、實(shí)際單位經(jīng)濟(jì)貢獻(xiàn)法、分?jǐn)傁禂?shù)法及缺水損失法等。在實(shí)際運(yùn)用中,這些方法都存在一定不足,如殘值法的影子價(jià)格較難測算,其實(shí)際操作性不強(qiáng);實(shí)際單位經(jīng)濟(jì)貢獻(xiàn)法是殘值法的改進(jìn),增加了可操作性,但該方法所計(jì)算的工業(yè)用水效益偏大。綜合考慮各方法的優(yōu)缺點(diǎn),本文選用經(jīng)濟(jì)學(xué)理論中的柯布-道格拉斯(C-D)生產(chǎn)函數(shù)法對(duì)工業(yè)供水效益進(jìn)行量化計(jì)算,該方法考慮了工業(yè)生產(chǎn)中非水要素的投入回報(bào),同時(shí)通過折算用水彈性消除了科技技術(shù)進(jìn)步、政策調(diào)整等對(duì)工業(yè)總產(chǎn)值的影響,使工業(yè)用水效益計(jì)算結(jié)果更為準(zhǔn)確[9]。

1)C-D生產(chǎn)函數(shù)構(gòu)建。在計(jì)算工業(yè)供水效益時(shí),只考慮工業(yè)供水對(duì)工業(yè)生產(chǎn)的影響,其他投入不做研究,將水作為第三種投入列入生產(chǎn)函數(shù)投入組合之中。由此,參考C-D函數(shù)形式,考慮工業(yè)供水的C-D生產(chǎn)函數(shù)可表示為:

Q=A·Kα·Lβ·Wλ。

(4)

式中:Q為工業(yè)總產(chǎn)值;A為效率系數(shù);K為工業(yè)固定資產(chǎn);L為工業(yè)勞動(dòng)力;W為工業(yè)用水量;α為固定資產(chǎn)彈性;β為勞動(dòng)力彈性;λ為用水彈性。

2)用水彈性計(jì)算。對(duì)式(4)構(gòu)建的生產(chǎn)函數(shù)求偏導(dǎo)數(shù),得到工業(yè)用水彈性,為了消除規(guī)模彈性的影響,引入折算用水彈性λ′:

(5)

(6)

3)效益計(jì)算。對(duì)工業(yè)用水求偏導(dǎo)數(shù)得出單位工業(yè)用水效益BW,根據(jù)以下公式求得洪水資源工業(yè)用水效益:

(7)

E工=BW·W工。

(8)

式中:E工為洪水資源工業(yè)用水效益;BW為單位工業(yè)用水效益;W工為洪水資源分配給工業(yè)用水的水量。

1.2.2 分?jǐn)傁禂?shù)法量化農(nóng)業(yè)灌溉效益

《水利建設(shè)項(xiàng)目經(jīng)濟(jì)評(píng)價(jià)規(guī)范》(SL 72—2013)中,推薦采用分?jǐn)傁禂?shù)法、影子水價(jià)法、缺水損失法計(jì)算農(nóng)業(yè)灌溉效益。其中,影子水價(jià)法理論清晰,但需要采用研究區(qū)的大量社會(huì)、經(jīng)濟(jì)、環(huán)境等資料進(jìn)行測算,具有一定的地域局限性;缺水損失法符合經(jīng)濟(jì)理論,但需要采用大量長系列資料,且涉及優(yōu)化分配問題。本文選用運(yùn)用最為廣泛的分?jǐn)傁禂?shù)法來量化農(nóng)業(yè)灌溉效益,其關(guān)鍵在于效益分?jǐn)傁禂?shù)的確定[10],可通過統(tǒng)計(jì)法與試驗(yàn)法確定分?jǐn)傁禂?shù),因試驗(yàn)法存在所得結(jié)果局限性較大、試驗(yàn)造價(jià)高等問題,采用統(tǒng)計(jì)法更為方便。農(nóng)業(yè)灌溉效益量化公式為:

E農(nóng)=Vqk=(I/W)fqk。

(9)

式中:E農(nóng)為洪水資源農(nóng)業(yè)灌溉效益;V為單方水價(jià)值;q為增供水量;k為水資源利用率;I為產(chǎn)業(yè)增加值;W為總用水量;f為分?jǐn)傁禂?shù)。

1.2.3 能值水價(jià)法量化生活供水效益

目前,尚無生活用水效益計(jì)算方法的具體規(guī)范,生活用水效益常用水價(jià)法來粗略估算[11]。常用的水價(jià)計(jì)算方法沒有體現(xiàn)出生活用水的資源價(jià)值、工程價(jià)值和環(huán)境價(jià)值,運(yùn)用能值價(jià)值理論計(jì)算水價(jià),可以計(jì)算出水資源作為一種自然環(huán)境資源進(jìn)入經(jīng)濟(jì)社會(huì)的真實(shí)、完整的價(jià)格,用以量化的洪水資源對(duì)生活供水效益更為全面完整。

生活用水的能值投入包括自然環(huán)境資源投入及社會(huì)經(jīng)濟(jì)資源投入,通過能值轉(zhuǎn)換率將各項(xiàng)投入指標(biāo)轉(zhuǎn)換為統(tǒng)一的能值,從而可得到生活供水全部的能值投入,根據(jù)能值貨幣比率將總能值以貨幣形式表達(dá),即可得到能值計(jì)算方法下的水價(jià)。水價(jià)的計(jì)算公式為:

P=Pz+Pg+Ph，

(10)

(11)

式中:P為居民生活用水水價(jià);Pz為資源水價(jià);Pg為工程水價(jià);Ph為環(huán)境水價(jià);Cyzr為水資源形成過程中投入的物質(zhì)的能值;Cgzr、Cgsj分別為供水過程中所投入自然環(huán)境資源和社會(huì)經(jīng)濟(jì)資源的能值;Cpzr、Cpsj分別為污水處理過程中所投入自然環(huán)境資源和社會(huì)經(jīng)濟(jì)資源的能值;Qg為城市居民生活供水量;Qp為城市居民生活污水排放量;EDR為該市當(dāng)年能值與貨幣比率。

在能值理論下的水價(jià)計(jì)算的基礎(chǔ)上,采用水價(jià)法計(jì)算生活用水效益,即城市居民用水價(jià)格與洪水資源生活用水可利用量之積視為生活用水效益,其計(jì)算公式如下:

E生=P·W生。

(12)

式中:E生為城市居民生活用水效益;W生為洪水資源分配給生活用水的水量。

1.2.4 能值法量化生態(tài)供水效益

生態(tài)用水效益是洪水資源利用效益的重要組成部分,因?yàn)槠鋸?fù)雜性,且涉及范圍較廣,量化較為困難。目前,研究者們提出的常用量化方法主要有影子價(jià)格法、替代工程法、費(fèi)用分析法和旅行費(fèi)用法等。這些方法一般從經(jīng)濟(jì)學(xué)的角度出發(fā),只針對(duì)生態(tài)環(huán)境某一方面的效益進(jìn)行量化,沒有從生態(tài)經(jīng)濟(jì)系統(tǒng)整體出發(fā)。而能值法可根據(jù)生態(tài)環(huán)境用水的特點(diǎn)和種類,在定性分析生態(tài)環(huán)境用水產(chǎn)生的效益的基礎(chǔ)上,對(duì)其各效益因子分別進(jìn)行量化。

洪水資源產(chǎn)生的生態(tài)環(huán)境用水效益，可由洪水資源分配給生態(tài)環(huán)境的用水量占總生態(tài)環(huán)境用水量的比值乘以總生態(tài)環(huán)境用水效益得到[12],即

(13)

式中:EMf為洪水資源生態(tài)用水效益;W生為洪水資源分配給生態(tài)環(huán)境的用水量；τs為物種的能值轉(zhuǎn)換率;N為研究區(qū)域內(nèi)水生生物種數(shù);R為水生生物活動(dòng)面積占全球面積的比例；mi為水體對(duì)第i種污染物的納污能力;τi為第i種污染物的能值轉(zhuǎn)換率；E為水體蒸發(fā)能量;τq為蒸汽能值轉(zhuǎn)化率；Wt為湖泊或濕地的年蓄水量;τt為相應(yīng)水體的能值轉(zhuǎn)化率。

2 實(shí)例研究

2.1 連云港市概況

連云港市位于江蘇省東北部,地處淮河流域沂沭泗河水系的最下游,流域主要洪水入海通道新沂河、新沭河經(jīng)連云港市入海,承擔(dān)上游近80 000 km2流域面積的泄洪任務(wù)。有新沂河、新沭河、龍王河3條過境河道,當(dāng)?shù)囟嗄昶骄^境洪水量為60.48億m3,為洪水資源利用提供了水量保證。該市的水庫湖泊眾多,總數(shù)為江蘇省內(nèi)首位。大型水庫有3座,分別為石梁河水庫、塔山水庫和安峰山水庫。其中,石梁河水庫擔(dān)負(fù)著調(diào)蓄洪水的重要職能,水庫集水面積約為15 365 km2,總庫容5.31億m3,設(shè)計(jì)洪水位26.81 m,校核洪水位27.95 m,最大下泄流量10 131 m3/s。

2.2 洪水資源利用風(fēng)險(xiǎn)效益估算

2.2.1 工業(yè)供水效益

根據(jù)《連云港市統(tǒng)計(jì)年鑒》《連云港市水資源公報(bào)》及Wind數(shù)據(jù)庫[13-14],統(tǒng)計(jì)2000—2012年間連云港市的工業(yè)總產(chǎn)值、工業(yè)固定資產(chǎn)投資凈值、工業(yè)勞動(dòng)力投入及工業(yè)用水總量。將工業(yè)固定資產(chǎn)投資、工業(yè)勞動(dòng)力和工業(yè)用水量作為自變量,以工業(yè)總產(chǎn)值為因變量,分別取自然對(duì)數(shù),運(yùn)用EViews軟件對(duì)數(shù)據(jù)進(jìn)行回歸分析,可得,連云港市工業(yè)用水效益的生產(chǎn)函數(shù)和折算用水彈性為:

lnQ=-2.456+0.717lnK+0.653lnL+1.126lnW，

根據(jù)式(5)可以得到工業(yè)總產(chǎn)值的效益增量,再由式(7)得到每年的估算單方水效益和供水總效益,連云港市多年平均工業(yè)用水單方水效益約為2.4元/m3。

2.2.2 農(nóng)業(yè)灌溉效益

根據(jù)連云港市2014年統(tǒng)計(jì)年鑒與水資源公報(bào)[13-14],連云港市2014年的農(nóng)業(yè)產(chǎn)業(yè)增加值為12.4億元,農(nóng)業(yè)灌溉用水量約22.1億m3;根據(jù)相關(guān)研究資料,該市的農(nóng)業(yè)用水效益分?jǐn)傁禂?shù)取值0.45;水資源利用率取值0.9。因此,農(nóng)業(yè)單方水效益為:

0.227(元/m3)。

2.2.3 生活供水效益

運(yùn)用能值理論,對(duì)連云港市生活用水水價(jià)進(jìn)行計(jì)算,能值轉(zhuǎn)化率由Buenfil[15]提供的計(jì)算方法得到,以連云港市供水系統(tǒng)能值投入計(jì)算為例進(jìn)行計(jì)算。由相關(guān)資料可知,2013年連云港市生活用水量為2.01億m3,生活污水排放量為1.09億m3,根據(jù)式(11)可得到資源水價(jià)、工程水價(jià)、環(huán)境水價(jià)分別為:0.22、4.98、0.85元/m3,則水價(jià)為:

P=Pz+Pg+Ph=6.05(元/m3)。

2.2.4 生態(tài)供水效益

結(jié)合連云港市牛山站、青口站及石梁河測站的實(shí)測蒸發(fā)數(shù)據(jù)、部分骨干河道資料、《連云港市水資源綜合規(guī)劃》等相關(guān)資料[13-14],計(jì)算連云港生態(tài)用水量及能值貨幣比率，得到連云港生態(tài)環(huán)境單方用水效益為:

E生單= 3.21×1021/(7.44×1011×1.085×108)=

39.77(元/m3)。

根據(jù)2014年連云港市統(tǒng)計(jì)年鑒與水資源公報(bào)數(shù)據(jù),該市水資源在工業(yè)、農(nóng)業(yè)、生活及生態(tài)供水的分配比例分別為16.5%、69.0%、12.1%、2.4%;洪水資源利用分配比例也取以上比例。由不同汛限水位得到增蓄的水量,水量按比例分配后與各類單方水效益的乘積,即洪水資源利用該類風(fēng)險(xiǎn)效益值。石梁河水庫原汛限水位為23.5 m,按水位每抬高0.2 m計(jì)算效益值。工業(yè)、農(nóng)業(yè)灌溉、生活、生態(tài)用水效益計(jì)算見表1,各類效益值隨水位抬高的變化曲線如圖1所示。

表1 洪水資源利用風(fēng)險(xiǎn)效益

圖1 汛限水位抬高時(shí)各類用水效益變化曲線

由表1和圖1可知,隨著水庫汛限水位的抬高,水庫蓄存水量增多,使得各類用水效益值隨之增大。其中,生態(tài)供水分配比例最小,用水量最少,但是效益值增長速率最快,效益與水量投入比最大;而農(nóng)業(yè)用水量最多,但效益值最小,隨汛限水位抬高的增長速率也最小。

不同的汛限水位對(duì)應(yīng)的風(fēng)險(xiǎn)率不同。運(yùn)用數(shù)理統(tǒng)計(jì)方法,以石梁河水庫上游來水的歷時(shí)統(tǒng)計(jì)資料為基礎(chǔ)進(jìn)行調(diào)洪演算,找出汛限水位對(duì)應(yīng)的使水庫達(dá)到設(shè)計(jì)標(biāo)準(zhǔn)水位的洪水,以該洪水發(fā)生概率作為汛限水位對(duì)應(yīng)的大壩安全風(fēng)險(xiǎn)率;以水庫下泄流量不超過新沭河安全泄量為準(zhǔn),計(jì)算下泄流量超過下游河道安全泄量的風(fēng)險(xiǎn)率。在計(jì)算中,為了保證不與大壩安全風(fēng)險(xiǎn)率計(jì)算重合,使調(diào)洪演算的最高水位不高于水庫設(shè)計(jì)標(biāo)準(zhǔn)值26.81 m。按式(1)和式(2)分別計(jì)算大壩安全風(fēng)險(xiǎn)率和下泄流量超下游河道安全泄量風(fēng)險(xiǎn)率,按式(3)計(jì)算綜合風(fēng)險(xiǎn)率,根據(jù)相關(guān)資料,大壩安全風(fēng)險(xiǎn)造成的風(fēng)險(xiǎn)損失占總風(fēng)險(xiǎn)損失的比重較大,因此,大壩安全風(fēng)險(xiǎn)率權(quán)重值取0.6,下泄流量超下游河道安全泄量風(fēng)險(xiǎn)率權(quán)重值取0.4。得到風(fēng)險(xiǎn)率計(jì)算結(jié)果見表2,風(fēng)險(xiǎn)隨水位抬高的變化曲線如圖2所示。綜合風(fēng)險(xiǎn)率與綜合效益的關(guān)系如圖3所示。

表2 綜合風(fēng)險(xiǎn)效益與風(fēng)險(xiǎn)率關(guān)系

圖2 汛限水位抬高時(shí)各類風(fēng)險(xiǎn)率變化曲線

由表2和圖2可知,隨著水庫汛限水位的抬高,大壩安全風(fēng)險(xiǎn)率和下泄流量超下游安全泄量風(fēng)險(xiǎn)率逐漸增加,增長速率也逐漸加快。當(dāng)汛限水位在24.9 m以下變動(dòng)時(shí),大壩安全風(fēng)險(xiǎn)率變化平緩,風(fēng)險(xiǎn)率值不超過5%,認(rèn)為取24.9 m汛限水位值，遇20年一遇洪水,能夠保證水庫大壩安全風(fēng)險(xiǎn)處于可控范圍。經(jīng)過調(diào)洪演算分析,水庫下泄流量對(duì)不同汛限水位對(duì)應(yīng)的風(fēng)險(xiǎn)率值影響明顯,汛限水位達(dá)到24.5 m以上時(shí),下泄流量每增加200 m3/s,可以使風(fēng)險(xiǎn)率減少5%左右。如果在保證大壩防洪安全風(fēng)險(xiǎn)率不變的情況下，繼續(xù)抬高汛限水位，必須增大水庫下泄流量,會(huì)增加下游河道的防洪風(fēng)險(xiǎn),使洪水利用風(fēng)險(xiǎn)從大壩向下游河道堤防發(fā)生轉(zhuǎn)移。

由圖3可知,總效益值隨風(fēng)險(xiǎn)率的增加而增長,當(dāng)風(fēng)險(xiǎn)率控制在5%以內(nèi)(即水位不超過23.9 m)時(shí),效益增長速率很快；超過5%后,效益增長速率逐漸變慢。綜上所述,在考慮風(fēng)險(xiǎn)最小且效益最大的情況下,汛限水位應(yīng)控制在23.9～24.9 m為宜。

圖3 綜合風(fēng)險(xiǎn)效益與風(fēng)險(xiǎn)率變化曲線

3 結(jié)語

在考慮洪水風(fēng)險(xiǎn)因素前提下,對(duì)洪水資源利用效益進(jìn)行了量化研究。主要考慮了大壩安全風(fēng)險(xiǎn)因素和下泄流量超下游安全泄量兩種風(fēng)險(xiǎn)因素,其中大壩安全風(fēng)險(xiǎn)因素是主要風(fēng)險(xiǎn)因素。將汛限水位調(diào)整后的增蓄水量按工業(yè)、農(nóng)業(yè)、生活、生態(tài)分配,分別進(jìn)行效益量化,工業(yè)用水效益用C-D生產(chǎn)函數(shù)法量化,農(nóng)業(yè)灌溉效益用分?jǐn)傁禂?shù)法量化,生活和生態(tài)用水效益用能值法量化。結(jié)合連云港市的石梁河水庫進(jìn)行分析,隨著汛限水位的抬高,風(fēng)險(xiǎn)效益和風(fēng)險(xiǎn)率都呈現(xiàn)增長趨勢,風(fēng)險(xiǎn)效益增長速率比較穩(wěn)定,而風(fēng)險(xiǎn)率增長速率則逐漸加快。根據(jù)水位、效益、風(fēng)險(xiǎn)三者的關(guān)系曲線圖,考慮風(fēng)險(xiǎn)最小且效益最大的情況下,給出了是石梁河水庫汛限水位合理的控制范圍。對(duì)于如何確定具體的汛限水位值,可以以風(fēng)險(xiǎn)最小和效益最大為目標(biāo),建立多目標(biāo)決策模型,有待進(jìn)一步研究。

[1]王宗志,程亮,劉友春,等．流域洪水資源利用的現(xiàn)狀與潛力評(píng)估方法[J].水利學(xué)報(bào),2004,45(4):474-481.

[2]王銀堂,胡慶芳,張書函,等.流域雨洪資源利用評(píng)價(jià)及利用模式研究[J].中國水利,2009(15):13-16.

[3]BROUWER R,EK R V.Integrated ecological,economic and social impact assessment of alternative flood control policies in the Netherlands[J].Ecological Economies,2004,50(1/2):1-21.

[4]殷峻逞,侯召成.一種基于供水風(fēng)險(xiǎn)的洪水資源效益分析方法[C]∥中國自然資源學(xué)會(huì)2004年學(xué)術(shù)年會(huì)論文集(上冊),2004.

[5]胥衛(wèi)平,李括.城市雨洪資源利用經(jīng)濟(jì)價(jià)值評(píng)價(jià)研究[J].環(huán)境科學(xué)與技術(shù),2009,32(7)：190-194.

[6]馮峰,何宏謀.河流洪水資源利用的社會(huì)效益分度測評(píng)定量評(píng)估研究[J].水力發(fā)電學(xué)報(bào),2013,32(4):25-31.

[7]FENG P,HAN S.Flood risk analysis of heightening limited water level of reservoir in flood season[J].Journal of Tianjin University Science and Technology,2007,40(5):525-529.

[8]黃顯峰,陳穎欽,方國華,等.洪水資源利用多目標(biāo)風(fēng)險(xiǎn)決策模型研究[J].珠江現(xiàn)代建設(shè),2014(6):1-3.

[9]趙寶璋,王曉妍.水利工程供水經(jīng)濟(jì)效益計(jì)算方法研究[J].水利經(jīng)濟(jì),2000(6):31-38.

[10]王震.基于分?jǐn)傁禂?shù)法估算治太工程的灌溉經(jīng)濟(jì)效益[J].灌溉排水學(xué)報(bào),2011,30(5):139-142.

[11]朱麗向,黃顯峰.連云港市雨洪資源利用模式及效益分析[J].水利科技與經(jīng)濟(jì),2014,20(10):21-24.

[12]樊艷麗,聶華.山西省水資源產(chǎn)業(yè)貢獻(xiàn)度及生態(tài)用水效益分析[J].中國水土保持科學(xué),2014(2):99-104.

[13]連云港市統(tǒng)計(jì)局.連云港統(tǒng)計(jì)年鑒2014[M].連云港:連云港市統(tǒng)計(jì)局,2014.

[14]連云港市水利局.2014年連云港水資源公報(bào)[Z].連云港:連云港市水利局,2014.

[15]BUENFIL A A.Emergy evaluation of water[D].Gainesville,U.S.:Department of Environmental Engineering Sciences,University of Florida,2001.

(責(zé)任編輯：宰松梅)

Study on Risk-benefit Quantification of Floodwater Resources Utilization

HUANG Xianfeng1, HUANG Xueqing1, FANG Guohua1, WANG Zongzhi2, ZHU Lixiang3

(1.College of Water Resources and Hydropower Engineering, Hohai University, Nanjing 210098, China; 2.Nanjing Hydraulic Research Institute, Nanjing 210029, China; 3.Tongyu River North-extended Water Supply Project Management Division of Lianyungang, Lianyungang 222004, China)

This paper combines the risk factors of floodwater resources utilization to quantify the risk-benefit of floodwater resources after adjusting limit water level. The limit water level is a key level involved the flood risk and benefits of the reservoir. After an adjustment of limited water level, each risk rate correspond to each water level is obtained. The increasing water storage capacity due to the adjustment of the flood control level is divided in portion. Taking into account the industrial, agricultural, living, ecological risk benefits of four areas, the 4 kinds of benefits are quantified by C-D production function, the sharing coefficient method, water price based on energy and the energy theory. Taking Shilianghe reservoir of Lianyungang City as the research case, the risk rate and benefit value corresponding to each flood limit water level are calculated. And according to the relationship between flood control level, risk ratio and risk benefit, in considering the minimum risk and maximum benefit under the premise, the results show the reasonable control range of flood control level of Shilianghe reservoir is from 23.9 m to 24.9 m.

flood utilization; risk; benefit; C-D production function; energy analysis method

2016-09-06

國家重點(diǎn)研發(fā)計(jì)劃課題(2016YFC0400909);江蘇省自然科學(xué)基金項(xiàng)目(BK20130849);江蘇省水利科技項(xiàng)目(2014064)。

黃顯峰(1980—),男,湖北黃岡人,副教授,碩導(dǎo),主要從事水資源規(guī)劃、水庫優(yōu)化調(diào)度、洪水資源利用方面的研究。E-mail:hxfhuang2005@163.com。

10.3969/j.issn.1002-5634.2016.06.009

TV213

A

1002-5634(2016)06-0049-06