簡(jiǎn)介: 本文基于洪水風(fēng)險(xiǎn)的定義,從洪災(zāi)發(fā)生的概率,洪災(zāi)后果評(píng)估,風(fēng)險(xiǎn)的評(píng)價(jià)指標(biāo),防洪標(biāo)準(zhǔn)經(jīng)濟(jì)優(yōu)化決策原則等方面闡述了中荷兩國(guó)在洪水風(fēng)險(xiǎn)分析研究方面的異同,初步分析了存在差別的原因以及兩種方法的特點(diǎn)。在分別應(yīng)用中國(guó)和荷蘭現(xiàn)有風(fēng)險(xiǎn)評(píng)估方法分析安慶市堤防圈的洪水風(fēng)險(xiǎn)的基礎(chǔ)上,探討了將兩國(guó)的風(fēng)險(xiǎn)分析方法進(jìn)行整合的思路。算例表明,將堤防的工程風(fēng)險(xiǎn)與洪水本身出現(xiàn)的水文風(fēng)險(xiǎn)相結(jié)合,并考慮損失在不同洪水位下的差異的方法是與洪災(zāi)發(fā)生機(jī)理相符的,綜合考慮了各種不確定因素的,切合我國(guó)洪災(zāi)特點(diǎn)的一種值得進(jìn)一步深入研究的新方法。
關(guān)鍵字:洪水風(fēng)險(xiǎn)、概率、洪災(zāi)損失、經(jīng)濟(jì)優(yōu)化
1 前言
2003年9月2日至12月1日在荷蘭交通水運(yùn)部進(jìn)行了為期三個(gè)月的學(xué)習(xí),作為水利部防洪減災(zāi)中心與荷蘭交通水運(yùn)部(DWW)合作研究項(xiàng)目的一部分內(nèi)容,雙方交換洪水風(fēng)險(xiǎn)研究領(lǐng)域的相關(guān)知識(shí),在查閱荷方提供的大量有關(guān)洪水風(fēng)險(xiǎn)文獻(xiàn)的基礎(chǔ)上[1-7],結(jié)合自己在洪災(zāi)風(fēng)險(xiǎn)分析方面的研究實(shí)踐和一些體會(huì),對(duì)中荷雙方洪水風(fēng)險(xiǎn)分析技術(shù)進(jìn)行了比較研究,考察荷蘭的工程可靠性分析與國(guó)內(nèi)常規(guī)洪水風(fēng)險(xiǎn)分析相結(jié)合的可能性。選擇長(zhǎng)江流域下游的安慶堤防作為研究實(shí)例,運(yùn)用中國(guó)、荷蘭的常規(guī)洪水分析方法分別計(jì)算了安慶堤防保護(hù)區(qū)的洪水風(fēng)險(xiǎn),并初步探討了將兩國(guó)方法整合起來,適宜我國(guó)洪災(zāi)特點(diǎn)的洪水風(fēng)險(xiǎn)評(píng)價(jià)方法。
2 洪水風(fēng)險(xiǎn)的定義
不同的人對(duì)風(fēng)險(xiǎn)具有不同的理解,有的人認(rèn)為風(fēng)險(xiǎn)是事件發(fā)生的概率,有的人認(rèn)為風(fēng)險(xiǎn)是事件發(fā)生的結(jié)果,還有的人認(rèn)為風(fēng)險(xiǎn)是危害性及其后果的產(chǎn)物。在荷蘭洪水風(fēng)險(xiǎn)被明確的定義為洪災(zāi)發(fā)生的概率與其所導(dǎo)致后果的乘積。因此洪災(zāi)發(fā)生概率的確定以及洪災(zāi)發(fā)生后果的評(píng)估是荷蘭洪水風(fēng)險(xiǎn)分析的重點(diǎn)。
我國(guó)對(duì)洪水風(fēng)險(xiǎn)定義尚不明確,洪水發(fā)生頻率、工程安全風(fēng)險(xiǎn),洪水水力特征(水深、淹沒歷時(shí)、流速、洪水到達(dá)時(shí)間等)、洪災(zāi)期望損失等都曾被用來描述洪水風(fēng)險(xiǎn),其中洪災(zāi)期望損失是現(xiàn)階段刻畫洪水風(fēng)險(xiǎn)的最高層次指標(biāo)。該指標(biāo)定義與荷蘭的洪水風(fēng)險(xiǎn)定義相符,為了便于比較,在本文后續(xù)的敘述和討論中,洪水風(fēng)險(xiǎn)都定義為洪災(zāi)發(fā)生概率與洪災(zāi)后果的乘積,也即洪災(zāi)期望損失。
3 洪災(zāi)發(fā)生概率
荷蘭已成海堤與河堤的設(shè)計(jì)主要基于超過某一特定頻率的水位,即所謂的超頻率法(Probability of excceedance),采用水文頻率分析方法得到超過某一特定頻率的水位(設(shè)計(jì)水位),再加上一定的安全超高確定堤防的高度。認(rèn)為這一特定的頻率近似表明了被保護(hù)區(qū)域的受淹可能性。20世紀(jì)90年代初,荷蘭防洪設(shè)施技術(shù)咨詢委員會(huì)(TAW)開始了洪水風(fēng)險(xiǎn)的研究,提出了洪災(zāi)發(fā)生概率的概念(Probability of flooding),一個(gè)區(qū)域發(fā)生洪災(zāi)的概率是指由于防洪保護(hù)設(shè)施破壞而使整個(gè)被保護(hù)區(qū)域遭受淹沒的概率.
它與超頻率洪水法的區(qū)別在于:
。1)將整個(gè)堤防圈作為研究對(duì)象,而不只是超頻率洪水法中所注重的單個(gè)堤防斷面;
。2)考慮了堤防的各種失穩(wěn)模式,這與超頻率洪水法中主要考慮溢流和漫頂兩種失穩(wěn)模式顯著不同;
(3)預(yù)先系統(tǒng)地考慮了洪災(zāi)頻率計(jì)算過程中的不確定因素(洪水發(fā)生的隨機(jī)性,模型的可靠性以及數(shù)據(jù)的真實(shí)性等)。在超頻率洪水法中,往往通過在設(shè)計(jì)洪水位的基礎(chǔ)上另加一定的安全超高來考慮這些不確定因素[1]。TAW計(jì)算了荷蘭中部4個(gè)堤防圈的洪災(zāi)發(fā)生頻率,結(jié)果證明運(yùn)用洪災(zāi)概率法能夠較容易地辨識(shí)堤防圈強(qiáng)度的薄弱點(diǎn)(安全制約瓶頸),并進(jìn)行排序,幫助更好地了解防洪系統(tǒng)(或其部分)改善的重點(diǎn),為逐步提高堤防圈保護(hù)標(biāo)準(zhǔn)提供科學(xué)依據(jù),有利于提高防洪建設(shè)的投資效率。洪災(zāi)頻率計(jì)算是荷蘭洪水風(fēng)險(xiǎn)分析研究中至關(guān)重要的一部分。目前,荷蘭已經(jīng)在全國(guó)展開了計(jì)算53個(gè)堤防圈保護(hù)區(qū)洪災(zāi)發(fā)生概率的工作,并不斷改進(jìn)模型,使計(jì)算概率盡可能接近保護(hù)區(qū)實(shí)際受淹的概率。
我國(guó)堤防防洪標(biāo)準(zhǔn)的確定類似于荷蘭的超頻率洪水法,堤防的設(shè)防水位與一定的頻率相對(duì)應(yīng)。以設(shè)計(jì)水位或控泄流量的重現(xiàn)期來判別堤防保護(hù)區(qū)的受淹安全度。認(rèn)為一旦洪水超過了設(shè)計(jì)水位,即洪水的重現(xiàn)期大于堤防的防洪標(biāo)準(zhǔn)時(shí),保護(hù)區(qū)就會(huì)受淹,且其淹沒概率與超標(biāo)準(zhǔn)洪水的出現(xiàn)概率(水文頻率)相等;洪水位低于設(shè)計(jì)水位,洪水的重現(xiàn)期小于堤防的防洪標(biāo)準(zhǔn)時(shí),堤防則不發(fā)生破壞,保護(hù)區(qū)是不會(huì)被淹的。這種方法主要考慮了洪水的水文風(fēng)險(xiǎn)因子和水力風(fēng)險(xiǎn)因子,而忽視了工程風(fēng)險(xiǎn)因子。實(shí)際上,設(shè)計(jì)水位被超出的概率與淹沒概率并不相同。由于工程風(fēng)險(xiǎn)的存在,即使在洪水位低于設(shè)計(jì)水位時(shí),保護(hù)區(qū)也有被淹沒的可能,相反堤壩也經(jīng)常能夠很好地?fù)踝】刂扑徊⒕哂凶銐虻陌踩珒?chǔ)備。更接近于保護(hù)區(qū)實(shí)際受淹概率的算法是應(yīng)該全面考慮水文風(fēng)險(xiǎn)、水力風(fēng)險(xiǎn)以及防洪工程風(fēng)險(xiǎn)。
4 洪災(zāi)后果評(píng)估方法
洪災(zāi)的后果包括洪災(zāi)所引起的經(jīng)濟(jì)損失,人員傷亡,對(duì)環(huán)境、社會(huì)、文化等帶來的損害等。荷蘭開發(fā)了較為成熟的洪災(zāi)經(jīng)濟(jì)損失模型,并計(jì)劃出臺(tái)有關(guān)洪災(zāi)后果評(píng)估的操作規(guī)范[2]。
荷蘭經(jīng)濟(jì)損失評(píng)估公式為:
其中,——第類資產(chǎn)的洪災(zāi)損失系數(shù),與淹沒深度有關(guān)
——第類資產(chǎn)的數(shù)量
——單位類資產(chǎn)的估計(jì)最大損失值
荷蘭在綜合考慮國(guó)家統(tǒng)計(jì)局所采用的分類方法和財(cái)產(chǎn)承災(zāi)特性的基礎(chǔ)上,將受淹資產(chǎn)劃分為不同的類別,并盡可能建立資產(chǎn)損失類型與國(guó)家、地方或者部委的統(tǒng)計(jì)類型之間的一一對(duì)應(yīng)關(guān)系。每類資產(chǎn)的數(shù)量和單位最大可能損失值在相關(guān)的文件中可以方便地查閱到。損失系數(shù)相當(dāng)于我國(guó)的洪災(zāi)損失率,它揭示了材料的破壞程度,該系數(shù)取決于相應(yīng)的淹沒特征。根據(jù)歷史洪災(zāi)記錄(1953年歐洲大洪水)建立資產(chǎn)的損失系數(shù)與淹沒特征(淹沒水深)的函數(shù)關(guān)系,對(duì)于資料不足而難以建立損失系數(shù)函數(shù)的類型,則通過征詢建筑工程師、企業(yè)管理者和經(jīng)營(yíng)者的方式來近似確定。全國(guó)范圍內(nèi)采用統(tǒng)一的洪災(zāi)損失系數(shù)函數(shù)關(guān)系,也就是說,在同等淹沒條件下,分別位于荷蘭南部和荷蘭北部同等類型的房屋通過模型計(jì)算得到的洪水經(jīng)濟(jì)損失是相同的。另外,荷蘭還將被評(píng)估區(qū)域分為低風(fēng)險(xiǎn)區(qū)和高風(fēng)險(xiǎn)區(qū),洪水發(fā)生重現(xiàn)期大于25年一遇的區(qū)域?qū)儆诘惋L(fēng)險(xiǎn)區(qū),小于25年一遇的劃歸高風(fēng)險(xiǎn)區(qū)。高風(fēng)險(xiǎn)區(qū)域由于財(cái)產(chǎn)的密集度低,各種承災(zāi)體的承災(zāi)能力強(qiáng),防災(zāi)救災(zāi)措施完備,它在遭受同等大小的洪水時(shí)損失要比低風(fēng)險(xiǎn)區(qū)域的損失小,因此在實(shí)際計(jì)算中,按高風(fēng)險(xiǎn)區(qū)域的最大可能損失比低風(fēng)險(xiǎn)區(qū)的最大可能損失小25%來考慮這種差別。
荷蘭也將經(jīng)濟(jì)損失分為直接經(jīng)濟(jì)損失和間接損失兩大類。財(cái)產(chǎn)損失是直接經(jīng)濟(jì)損失,而農(nóng)業(yè)、公司企業(yè)的停產(chǎn)停業(yè)損失以及防洪救災(zāi)費(fèi)用屬間接損失。直接損失和間接損失的損失系數(shù)函數(shù)關(guān)系相同,但直接損失和間接損失的最大可能損失值在不同行業(yè)中有著較大差別(如表1)
表1 荷蘭洪災(zāi)最大可能損失值一覽表
Table 1 Maximum damage amounts in the Netherlands
損失類型 | 單位 | 最大損失數(shù)(荷蘭盾) | ||
直接損失 | 間接損失 | |||
土地利用 | 農(nóng)業(yè)直接損失 | M2 | 2.5 | 1 |
花房直接損失 | M2 | 50 | 20 | |
企業(yè)公司 | 工礦業(yè)直接損失 | Jobs | 394000 | 315000 |
建筑業(yè)直接損失 | Jobs | 96000 | 171000 | |
商業(yè)直接損失 | Jobs | 97000 | 331000 | |
金融保險(xiǎn)業(yè)直接損失 | Jobs | 156000 | 180000 | |
交通、通訊業(yè)直接損失 | Jobs | 450000 | 151000 |
荷蘭不僅采用經(jīng)濟(jì)損失描述洪災(zāi)的后果,而且將由洪災(zāi)引起的人員傷亡作為衡量洪災(zāi)后果的一個(gè)重要方面。認(rèn)為因?yàn)?zāi)死亡率與洪水淹沒深度和洪水流速有直接的關(guān)系,初步建立了人口洪災(zāi)死亡率與這兩個(gè)參數(shù)的函數(shù)關(guān)系,但沒有考慮居住區(qū)域的建筑物因素以及避難的可能性對(duì)洪水死亡率的影響。估計(jì)由洪水造成的人員死亡數(shù)字是一個(gè)較為復(fù)雜的問題,因?yàn)樗婕霸S多因素,其它區(qū)域的數(shù)據(jù)或歷史數(shù)據(jù)不能直接用于現(xiàn)在的荷蘭。此外,目前已知數(shù)據(jù)及其存在的死難者數(shù)目和相關(guān)因素的關(guān)系尚不能定量地確定下次洪水可能導(dǎo)致的死亡數(shù)目。因此,洪災(zāi)傷亡定量評(píng)估模型的研究只是初步的。
總體上來講,我國(guó)的洪災(zāi)損失評(píng)估方法與荷蘭的方法基本相同,但由于我國(guó)洪澇災(zāi)害的復(fù)雜性和受災(zāi)區(qū)域承災(zāi)體的多樣性,無論是受淹資產(chǎn)的統(tǒng)計(jì)還是洪災(zāi)損失率函數(shù)關(guān)系的確定都是基于某一具體的受災(zāi)區(qū)域的,也就是說在不同的區(qū)域間,同一類受淹資產(chǎn)的單價(jià)和損失率函數(shù)關(guān)系都可能是不同的。這種方法雖然具有針對(duì)性,但調(diào)查和計(jì)算的工作量都會(huì)成倍加大,不可避免地存在大量重復(fù)工作,況且事無俱細(xì)反而難以保證評(píng)估精度?梢钥紤]將全國(guó)劃分為不同的洪水災(zāi)害典型區(qū)域(流域),在每一類區(qū)域試采用統(tǒng)一的損失率函數(shù)關(guān)系和受淹資產(chǎn)統(tǒng)計(jì)標(biāo)準(zhǔn),具體區(qū)域都可參照此標(biāo)準(zhǔn)評(píng)估洪災(zāi)損失,這樣既簡(jiǎn)化了評(píng)估過程,又兼顧了我國(guó)洪水災(zāi)害的復(fù)雜性。在洪水損失率關(guān)系的確定方式上也有可借鑒之處,我國(guó)一般都是根據(jù)歷史典型洪水統(tǒng)計(jì)調(diào)查資料建立損失率關(guān)系,一方面這種資料非常有限,不可能收集到所有財(cái)產(chǎn)損失類型在各種不同水深和淹沒等級(jí)下的損失資料。另一方面,因?yàn)椴粩嗟挠懈鞣N減災(zāi)措施的實(shí)施,實(shí)際的洪災(zāi)損失并不完全反應(yīng)預(yù)測(cè)意義上的洪災(zāi)損失。應(yīng)用綜合的損失率確定方法,在進(jìn)行歷史洪水調(diào)查的同時(shí)結(jié)合工程人員、設(shè)計(jì)人員及相關(guān)專家的意見,既可以克服資料短缺的困境,又有利于提高洪水損失率的精確性。隨著社會(huì)經(jīng)濟(jì)的發(fā)展,人們?cè)絹碓疥P(guān)注生命的價(jià)值,為了全面反應(yīng)洪水災(zāi)害的后果,有必要對(duì)洪水災(zāi)害人員傷亡進(jìn)行定量評(píng)估。
5 洪水風(fēng)險(xiǎn)指標(biāo)
荷蘭的洪水風(fēng)險(xiǎn)指標(biāo)大體分為兩類:從經(jīng)濟(jì)損失來講的經(jīng)濟(jì)風(fēng)險(xiǎn)指標(biāo),從人員傷亡角度來講的個(gè)人風(fēng)險(xiǎn)指標(biāo)和社會(huì)風(fēng)險(xiǎn)指標(biāo)[3]。
用洪水所造成的經(jīng)濟(jì)損失來表征洪災(zāi)后果的風(fēng)險(xiǎn)指標(biāo)稱為經(jīng)濟(jì)風(fēng)險(xiǎn)指標(biāo),類似于我國(guó)的洪災(zāi)期望損失。不同點(diǎn)在于我國(guó)認(rèn)為淹沒的概率等于超標(biāo)準(zhǔn)洪水發(fā)生的概率,未考慮標(biāo)準(zhǔn)內(nèi)洪水的淹沒可能性。計(jì)算超標(biāo)準(zhǔn)洪水的洪災(zāi)期望損失時(shí),考慮了洪水頻率大小與損失的關(guān)系,洪水出現(xiàn)的頻率越小,對(duì)應(yīng)的洪量則越大,直接經(jīng)濟(jì)損失亦越大,這種方法從概率統(tǒng)計(jì)理論上講是可取的,并也在一定程度上反映了洪水損失隨機(jī)的特性,所以期望損失值的求取是一個(gè)積分的過程。荷蘭在計(jì)算經(jīng)濟(jì)風(fēng)險(xiǎn)時(shí),認(rèn)為一個(gè)堤防圈如果在最薄弱點(diǎn)處(或者其它某處)破壞,即使在來水量級(jí)并不相同的情況下,也假設(shè)所造成的洪災(zāi)損失是相等的。(如圖1)。采用這種假設(shè)的原因是,荷蘭地處低洼地區(qū),地勢(shì)相對(duì)平坦,堤防圈的保護(hù)范圍比較小,在來水很小的時(shí)候,整個(gè)保護(hù)圈就有可能全部被淹沒。尤其在比較洪水風(fēng)險(xiǎn)的相對(duì)大小而確定堤防圈的防護(hù)薄弱點(diǎn)時(shí),區(qū)分不同級(jí)別洪水造成洪水損失的現(xiàn)實(shí)意義不大。
如果將洪災(zāi)后果用人員死亡數(shù)來表達(dá),荷蘭引入個(gè)人風(fēng)險(xiǎn)和社會(huì)風(fēng)險(xiǎn)的概念來衡量洪水風(fēng)險(xiǎn)。并在理論上研究了個(gè)人可接受風(fēng)險(xiǎn)程度和社會(huì)可接受風(fēng)險(xiǎn)程度。將可接受風(fēng)險(xiǎn)的理論用于防洪工程設(shè)施標(biāo)準(zhǔn)的選擇等問題上,便于與其它社會(huì)問題進(jìn)行風(fēng)險(xiǎn)比較。
個(gè)人風(fēng)險(xiǎn)是基于事故統(tǒng)計(jì)資料的一個(gè)指標(biāo),個(gè)人可接受風(fēng)險(xiǎn)是社會(huì)對(duì)風(fēng)險(xiǎn)接受程度的最小部分,個(gè)人根據(jù)所要承擔(dān)的風(fēng)險(xiǎn)和可能帶來的利益評(píng)判風(fēng)險(xiǎn)是否可接受,其中個(gè)人承擔(dān)風(fēng)險(xiǎn)的自愿程度是一個(gè)很重要的方面。某項(xiàng)活動(dòng)發(fā)生事故的允許概率取決于該項(xiàng)活動(dòng)的年事故死亡率與表征參加活動(dòng)自愿程度的綜合系數(shù)的乘積。
社會(huì)風(fēng)險(xiǎn)描述了災(zāi)難發(fā)生的可能性與特定人群在該災(zāi)難中傷亡數(shù)的關(guān)系。年期望人員傷亡值是表述社會(huì)風(fēng)險(xiǎn)的最簡(jiǎn)捷的指標(biāo)。為了區(qū)分公眾對(duì)大概率小傷亡與小概率大傷亡事件接受程度的不同,也常采用年期望人員傷亡數(shù)與標(biāo)準(zhǔn)方差和表示風(fēng)險(xiǎn)好惡指數(shù)之積的和來表述社會(huì)風(fēng)險(xiǎn)。
6 堤防安全標(biāo)準(zhǔn)的選擇(經(jīng)濟(jì)優(yōu)化決策方法)
基于人員傷亡的個(gè)人和社會(huì)可接受風(fēng)險(xiǎn)程度對(duì)防洪標(biāo)準(zhǔn)選擇的影響目前仍在理論研究階段。荷蘭大部分堤防標(biāo)準(zhǔn)(主要指堤防高度)的選擇仍基于經(jīng)濟(jì)風(fēng)險(xiǎn)指標(biāo),即認(rèn)為最優(yōu)經(jīng)濟(jì)標(biāo)準(zhǔn)是當(dāng)整個(gè)系統(tǒng)的費(fèi)用達(dá)到最低。費(fèi)用包括堤防加高費(fèi)用與預(yù)期經(jīng)濟(jì)損失,這兩項(xiàng)都是洪水出現(xiàn)概率(水文頻率)的函數(shù),在總費(fèi)用與洪水出現(xiàn)概率(水文頻率)建立了函數(shù)關(guān)系之后,可以進(jìn)一步求解出最優(yōu)經(jīng)濟(jì)堤防保護(hù)標(biāo)準(zhǔn)以及相應(yīng)的堤防高度[4]。如果把人員傷亡按經(jīng)濟(jì)折價(jià)考慮到洪災(zāi)損失當(dāng)中去,會(huì)引起優(yōu)化破壞概率(防護(hù)標(biāo)準(zhǔn))的下降和堤壩最優(yōu)高度的提高,人的生命損失對(duì)堤壩高度的影響程度決定于人的生命價(jià)值的衡量標(biāo)準(zhǔn)。如果假設(shè)人的生命價(jià)值等于每位居民的凈國(guó)民生產(chǎn)總值的現(xiàn)金價(jià)格,其計(jì)算實(shí)例表明,其影響程度是相當(dāng)小的,即使增加十倍,這一影響仍是非常有限的。在進(jìn)行洪災(zāi)可接受風(fēng)險(xiǎn)程度評(píng)價(jià)時(shí),如何用經(jīng)濟(jì)術(shù)語來考慮人的生命損失仍在研討之中[3]。
我國(guó)在防洪工程方案優(yōu)選時(shí)常用到效益-費(fèi)用分析方法,防洪工程的效益指在有無工程條件下洪災(zāi)損失的差值。最優(yōu)的方案是凈效益(=效益-費(fèi)用)最大的方案。經(jīng)過公式推導(dǎo),凈效益最大的方案等同于工程投入與工程建成后期望損失值之和最小的方案。盡管表述方式有差別,荷蘭與我國(guó)防洪工程方案優(yōu)選的實(shí)質(zhì)是相同的。
7 安慶堤防保護(hù)范圍風(fēng)險(xiǎn)計(jì)算
選擇長(zhǎng)江下游的安慶市堤防作為典型研究區(qū)域,前期研究計(jì)算了基于不同洪水位的包含堤防圈各類失險(xiǎn)模式在內(nèi)的綜合工程風(fēng)險(xiǎn),確定了最易失險(xiǎn)的堤防斷面,假設(shè)堤防在此潰口,根據(jù)潰口處的流量過程以及安慶市的地形條件,近似框定洪水淹沒范圍和淹沒水深[8]。結(jié)合安慶市現(xiàn)狀社會(huì)經(jīng)濟(jì)發(fā)展?fàn)顩r調(diào)查資料(2001年),并對(duì)1954年和1998年安慶市洪水損失統(tǒng)計(jì)資料進(jìn)行分析的基礎(chǔ)上,運(yùn)用已有的較為成熟的損失評(píng)估模型計(jì)算了基于不同洪水位的洪水淹沒經(jīng)濟(jì)損失。分別運(yùn)用荷蘭的洪水風(fēng)險(xiǎn)計(jì)算方法、國(guó)內(nèi)傳統(tǒng)的洪水風(fēng)險(xiǎn)計(jì)算方法以及綜合工程失險(xiǎn)概率和洪水水文概率的方法對(duì)安慶市的洪澇風(fēng)險(xiǎn)進(jìn)行評(píng)估。
。1)在確定淹沒區(qū)的洪水風(fēng)險(xiǎn)時(shí),進(jìn)行了以下假設(shè):
l 洪水風(fēng)險(xiǎn)區(qū)在不同洪水位下的淹沒概率等于超洪水位的概率與潰口處該洪水位下堤防工程失險(xiǎn)概率之積。
l 在運(yùn)用荷蘭的方法時(shí),假設(shè)堤防保護(hù)范圍的總的淹沒概率等于不同洪水位下的淹沒概率之和,與總淹沒概率相對(duì)應(yīng)的經(jīng)濟(jì)損失采用洪水在設(shè)計(jì)水位時(shí)堤防潰破所造成的經(jīng)濟(jì)損失。
l 運(yùn)用國(guó)內(nèi)傳統(tǒng)方法計(jì)算時(shí),假設(shè)在設(shè)計(jì)水位(安慶堤防的設(shè)計(jì)水位是:19.29m)和設(shè)計(jì)水位以下的淹沒概率等于零,當(dāng)水位超過設(shè)計(jì)水位后,堤防保護(hù)范圍的淹沒概率等于洪水水位出現(xiàn)的水文概率。也就是說,只要洪水水位超過設(shè)計(jì)水位,堤防100%會(huì)失險(xiǎn)。
l 綜合考慮洪水位的出現(xiàn)概率和工程失險(xiǎn)概率,并考慮到不同洪水位下?lián)p失的差異時(shí),假設(shè)在洪水水位較低時(shí)(相當(dāng)于5年一遇),即使發(fā)生洪水,損失非常小,在進(jìn)行風(fēng)險(xiǎn)計(jì)算時(shí)可以忽略不計(jì)。
(2)計(jì)算結(jié)果
對(duì)洪水位進(jìn)行概率分析,得出超不同洪水位的頻率,經(jīng)工程可靠性分析后得出在不同等級(jí)洪水位下的堤防失險(xiǎn)概率[8]。依據(jù)假設(shè),求取超水位頻率與堤防失險(xiǎn)概率之積確定洪水風(fēng)險(xiǎn)區(qū)在不同洪水位下的淹沒概率。并運(yùn)用洪災(zāi)損失評(píng)估模型計(jì)算了洪水風(fēng)險(xiǎn)區(qū)對(duì)應(yīng)于不同洪水位的淹沒損失,如表2所示。
表2 不同洪水位下淹沒概率及洪災(zāi)損失計(jì)算結(jié)果
Table 2 Flooding probabilities and damage estimates for the given flood events
洪水位(m) | (1) | 17.29 | 18.29 | 18.79 | 19.29 | 20.79 | 21.19 |
超洪水位頻率 | (2) | 0.100 | 0.036 | 0.014 | 0.010 | 0.003 | 0.001 |
工程失險(xiǎn)概率 | (3) | 0.1578 | 0.3450 | 0.4382 | 0.5017 | 0.9474 | 0.9832 |
淹沒概率 | (4)=(2)×(3) | 0.0158 | 0.0124 | 0.0061 | 0.0044 | 0.0028 | 0.001 |
經(jīng)濟(jì)損失(×104 Yuan) | (5) | 3175.7 | 3969.14 | 5152.23 | 6707.11 | 13420 | 21949.6 |
運(yùn)用三種不同的方法計(jì)算洪災(zāi)風(fēng)險(xiǎn)時(shí),所采用的公式分別是:
荷蘭方法:
國(guó)內(nèi)傳統(tǒng)方法:
新方法:
式中,,洪水位
,設(shè)計(jì)洪水位
,洪水位概率分布函數(shù)
,風(fēng)險(xiǎn)區(qū)淹沒概率函數(shù)
,洪災(zāi)損失值,是洪水位的函數(shù)
,與設(shè)計(jì)洪水位相對(duì)應(yīng)的洪災(zāi)損失
可以看出,新方法嘗試將荷蘭方法與國(guó)內(nèi)傳統(tǒng)方法相結(jié)合,與荷蘭方法的不同點(diǎn)在于認(rèn)為洪災(zāi)損失是洪水位的函數(shù)又考慮洪水量級(jí)不同所造成損失的差異,而與國(guó)內(nèi)傳統(tǒng)方法的不同在于考慮了當(dāng)洪水位低于設(shè)計(jì)洪水位時(shí)由于工程風(fēng)險(xiǎn)可能造成的洪災(zāi)損失,也即考慮標(biāo)準(zhǔn)內(nèi)的洪水風(fēng)險(xiǎn)。洪災(zāi)風(fēng)險(xiǎn)計(jì)算結(jié)果如表3所示:
表3 洪災(zāi)風(fēng)險(xiǎn)計(jì)算結(jié)果表
Table 3 The results of flood risk assessment
評(píng)價(jià)方法 | 荷 蘭 | 中 國(guó) | 新 方 法 |
洪災(zāi)風(fēng)險(xiǎn)(104 Yuan /y) | 99 | 95 | 154 |
。3)討論
算例研究的主要目的是實(shí)際應(yīng)用中荷兩國(guó)的洪水風(fēng)險(xiǎn)分析方法,計(jì)算過程中采用了簡(jiǎn)化的方法進(jìn)行洪水水位的頻率分析,洪水淹沒特征只是近似的框算,洪水損失計(jì)算結(jié)果不可避免地存在一定的誤差。用超洪水位頻率與該水位下工程失險(xiǎn)的概率之積來表示洪水淹沒概率的方法尚有待進(jìn)一步研究。運(yùn)用荷蘭的風(fēng)險(xiǎn)分析方法時(shí),采用設(shè)計(jì)洪水位下的洪災(zāi)損失與保護(hù)區(qū)的總淹沒概率相對(duì)應(yīng)也沒有進(jìn)行充分的論證。但可以明確的是:綜合中荷兩國(guó)現(xiàn)有風(fēng)險(xiǎn)評(píng)估方法的特點(diǎn),將堤防的工程風(fēng)險(xiǎn)與洪水本身出現(xiàn)的水文風(fēng)險(xiǎn)相結(jié)合,并考慮損失在不同洪水位下的差異的方法不失為一種更科學(xué)、更符合我國(guó)洪災(zāi)特點(diǎn)的方法,一方面,當(dāng)洪水位未達(dá)到設(shè)計(jì)水位時(shí),就存在淹沒的風(fēng)險(xiǎn),而且占有相當(dāng)比例,算例中,發(fā)生在標(biāo)準(zhǔn)之內(nèi)洪水風(fēng)險(xiǎn)約占總風(fēng)險(xiǎn)的1/3,應(yīng)該把這部分風(fēng)險(xiǎn)考慮在內(nèi)。我國(guó)實(shí)際發(fā)生的洪水災(zāi)害事件也證明了這一實(shí)事,1998年長(zhǎng)江大水影響較大的六處潰口中,有四處是在洪水水位并未達(dá)到設(shè)計(jì)水位的情況下由于堤基管涌而出險(xiǎn)的[9]。另外,我國(guó)與荷蘭洪災(zāi)的特點(diǎn)不同,防洪工程保護(hù)范圍較大,其內(nèi)地勢(shì)起伏較大,洪水水位的高低(進(jìn)入保護(hù)區(qū)洪量的大。┲苯雨P(guān)系到洪災(zāi)損失的大小,因此還有必要發(fā)展荷蘭的分析方法,考慮洪災(zāi)損失的這種隨機(jī)性。
8 結(jié)論
由于中荷兩國(guó)的洪水災(zāi)害特點(diǎn)存在著較大差異,形成了各具特色的洪水防御體系,在洪水風(fēng)險(xiǎn)分析技術(shù)上的側(cè)重點(diǎn)也有所不同。在分析我國(guó)洪災(zāi)具體特點(diǎn)的基礎(chǔ)上(洪災(zāi)損失的隨機(jī)性),借鑒荷蘭在洪災(zāi)分析方法上的先進(jìn)經(jīng)驗(yàn)(風(fēng)險(xiǎn)定義明確,淹沒概率分析全面),改進(jìn)常規(guī)分析方法的不足,是進(jìn)一步提高我國(guó)洪水風(fēng)險(xiǎn)分析技術(shù)的有效途徑。
參考文獻(xiàn)
[1] Technical Advisory Committee on Water Defence. From probability of exceedance to probability of flooding [M]. Delft, The Netherlands: Ministry of Transport, public Works and water Management, 2000. 3-5.
[2] N. Vrisou van Eck and M. Kok. Standard method for Predicting Damage and Casualties as a Result of Floods [M]. Delft, The Netherlands: Ministry of Transport, public Works and water Management, 2001.22-41.
[3] J.K. Vrijling, H.G. Voorman. Reliability-based design of flood defences [A]. Flood defence’2002 [C]. Beijing: Science Press, 2002, 123-133.
[4] S.N.Jonkman, P.H.A.J.M. van Gelder, J.K. Vrijling. Loss of life models for sea and river floods [A]. Flood defence’2002 [C]. Beijing: Science Press, 2002, 196-206.
[5]Dushmanta Dutta, Srikantha Herath, Katumi Musiake. A mathematical model flood loss estimation[J]. Journal of hydrology, 2003,277(2003): 24-49.
[6]J.K. Vrijling. Probabilistic design of water defense system in the Netherlands [J]. Reliability engineering and system safety, 2001,74 (2001): 337-344.
[7]S.E.van Manen, M.Brinkhuis. Pilot case flood risk [A]. Applications of statistics and probability in civil engineering [C]. Rotterdam: Millpress, 2003.881-887.
[8]Wu Xingzheng. Probabilistic risk analysis and safety evaluation of dikes [D]. Research report performed in DWW, the Netherlands: 2003. 97-105.
[9]洪慶余、羅鐘毓,等,長(zhǎng)江防洪與98大洪水[M].北京,中國(guó)水利水電出版社,1999.171-175