馬勇驥+曾凡永+李君

【摘 要】本文介紹了三種水文干旱指標(biāo)，闡述了水文干旱重現(xiàn)期計算的原理與方法，利用樟樹水文站1970-2009年旬徑流系列進(jìn)行實(shí)例研究，其結(jié)果表明：（1）各指標(biāo)對水文干旱的識別效果較為一致；（2）三個指標(biāo)得出干旱歷時進(jìn)行擬合均可采用指數(shù)分布，干旱歷時擬合函數(shù)的選擇存在不同；（3）徑流距平指標(biāo)與SDI指數(shù)應(yīng)用于基于Copula函數(shù)的干旱重現(xiàn)期計算效果相對較好，Z指數(shù)相對較差。

【關(guān)鍵詞】水文干旱；干旱指標(biāo)；診斷分析；干旱頻率

中圖分類號： P333.3 文獻(xiàn)標(biāo)識碼： A 文章編號： 2095-2457（2017）23-0009-004

【Abstract】In this paper， three kinds of hydrological and drought indices are introduced， and the principle and method of hydrological drought recurrence period are expounded. Using the runoff series of Camphor Hydrological Station from 1970 to 2009 for example， the results show that： （1） （2） The exponential distribution can be obtained by fitting the drought duration with the three indexes， and the choice of fitting function of drought duration varies； （3） The runoff anomaly index and SDI index are applied to the study based on Copula The function of the drought return period calculation effect is relatively good， Z index is relatively poor.

【Key words】Hydrological drought； Drought index； Diagnostic analysis； Drought frequency

干旱是一種偏離正常狀態(tài)的自然現(xiàn)象，其表現(xiàn)為水文循環(huán)的一個或多個要素長時間缺水，在人與自然密切相連的今天，干旱通常會給人類社會帶來一定的災(zāi)害即旱災(zāi)，其主要危害表現(xiàn)在：農(nóng)作物減產(chǎn)、城市供水危機(jī)、土地荒漠化、河湖水體污染加重等方面。水文干旱是氣象干旱延續(xù)和發(fā)展，與農(nóng)業(yè)干旱緊密相關(guān)，還是社會經(jīng)濟(jì)干旱的起因之一，水文干旱研究不僅有著重要的理論意義，而且有助于對其他類型干旱的理解和研究。

目前，相關(guān)學(xué)者針對極限水文干旱頻率與重現(xiàn)期的研究也取得了不少成果。如馮國章將最大正游程長概率密度函數(shù)轉(zhuǎn)換為最大負(fù)游程長概率密度函數(shù)，導(dǎo)出了最大負(fù)游程長概率密度函數(shù)，并將其分別作為極限水文干旱歷時的概率密度函數(shù)和概率分布函數(shù)，提出了水文干旱的劃分標(biāo)準(zhǔn)以及確定水資源抗旱年限的方法[1]。周玉良等利用徑流量距平百分率和徑流量累積頻率識別水文干旱特征變量，以適線法確定單個干旱特征變量的分布曲線之后，利用Copula函數(shù)構(gòu)建干旱歷時與干旱烈度間的聯(lián)合分布[2]。陸桂華等將基于帕爾默干旱指數(shù)（Palmer Drought Severity Index，PDSI）和標(biāo)準(zhǔn)化降水指數(shù)（Standardized Precipitation Index，SPI）的綜合氣象干旱指數(shù)應(yīng)用于游程理提取出干旱特征變量，并應(yīng)用Copula方法計算干旱歷時和干旱烈度的聯(lián)合分布概率，進(jìn)而推算干旱極限事件重現(xiàn)水平[3]。

水文干旱頻率與重現(xiàn)期的研究，對于指導(dǎo)抗旱救災(zāi)有著非常重要的實(shí)際應(yīng)用價值，但是目前應(yīng)用于水文干旱特征變量識別的指標(biāo)多種多樣，基于不同指標(biāo)識別出的水文干旱特征變量應(yīng)用于干旱頻率與重現(xiàn)期的計算效果的差異不得而知。

1 水文干旱指標(biāo)

1.1 徑流距平指標(biāo)

1.3 SDI指數(shù)

標(biāo)準(zhǔn)化徑流指數(shù)（SDI）與標(biāo)準(zhǔn)化降水指數(shù)（SPI）類似，兩者計算原理相同，由實(shí)測徑流資料就可以求得，標(biāo)準(zhǔn)化徑流指數(shù)能夠用于不同的區(qū)域旱澇狀況的比較，也能夠靈活的給出不同空間尺度的值。SDI基于長系列徑流資料，采用gamma函數(shù)、對數(shù)正態(tài)分布函數(shù)或其他函數(shù)分布擬合徑流系列，再經(jīng)正態(tài)標(biāo)準(zhǔn)化求得SDI的值，不同的取值表征不同的旱澇程度，如用對數(shù)正態(tài)分布對徑流進(jìn)行擬合，則SDI指數(shù)的計算公式如下：

2 水文干旱重現(xiàn)期分析方法

2.1 水文干旱特征變量

水文干旱具有持續(xù)性與水分短缺性特征，因而同其他各種類型干旱一樣皆存在干旱歷時、干旱強(qiáng)度、干旱烈度三個主要特征變量。圖1給出了根據(jù)游程理論所描述的水文干旱的特征變量示意圖。

圖1 水文干旱特征變量示意圖

其中X（t）為時段t時水文干旱指標(biāo)，X0為干旱指標(biāo)對應(yīng)的截取水平，若時段干旱指標(biāo)連續(xù)處于截取水平以下則可判斷發(fā)生干旱事件。以基于水文變量絕對值單因素水文干旱指標(biāo)為例，圖中兩個陰影部分分別為一次干旱事件，其中D為干旱歷時，M為干旱強(qiáng)度，S為干旱烈度，三個干旱特征的含義分述如下：

2.1.1 干旱歷時D

干旱歷時為干旱事件的持續(xù)時間，即為一次水文干旱事件從開始到結(jié)束的時間間隔。水文干旱歷時從時間尺度上來描述水文干旱事件，反應(yīng)的是水文干旱的在時間上的連續(xù)性。干旱歷時D=Ts-Te+1，其中Ts為水文干旱開始時段，Te為水文干旱結(jié)束時段，則

2.1.2 干旱強(qiáng)度M

干旱強(qiáng)度是指水文干旱事件歷時內(nèi)時段平均缺水程度，具體為干旱歷時內(nèi)，干旱指標(biāo)與截取水平的平均差值，如圖1中M所示。endprint

2.1.3 干旱烈度S

干旱烈度是指水文干旱事件的干旱歷時內(nèi)累計干旱缺水程度，干旱烈度是從干旱累積缺水程度上對干旱事件進(jìn)行描述。反映了干旱缺水嚴(yán)重程度，具體如圖2中陰影部分面積S，即為水文干旱事件的干旱烈度。

2.2 重現(xiàn)期計算

以上兩式中EL為干旱間隔的期望值（以年為單位），數(shù)值為干旱歷時期望與非干旱歷時期望之和，通過對比聯(lián)合重現(xiàn)期與同現(xiàn)重現(xiàn)期的公式可知對給定場次的水文干旱同現(xiàn)重現(xiàn)期的總是大于前者， 從有利于防旱抗旱的角度考慮，本文僅研究同現(xiàn)重現(xiàn)期。

3 實(shí)例研究

贛江是江西省第一大河流，為長江八大支流之一。流域內(nèi)總?cè)丝?250萬（2005年統(tǒng)計），耕地面積98.3萬公頃，其中水田88萬公頃，旱地9.4萬公頃。雖然贛江流域雨量充沛、水資源豐富，但自中華人名共和國成立上溯到宋朝咸平三年的1000年間，贛江流域共發(fā)生大小旱災(zāi)307年次。1943流域年發(fā)生特大旱災(zāi)，中華人民共和國成立后，發(fā)生特大干旱旱災(zāi)的年份有：1963年、1978年、1986年。1963年春旱、夏旱連秋旱。夏季降水量比歷時同期減少35%到56%，秋季降水量比歷史同期減少25%至33%，流域受災(zāi)面積達(dá)61.2萬公頃，糧食減產(chǎn)95萬噸。本文采用了樟樹水文站1970-2009年共計40年的徑流量系列進(jìn)行水文干旱頻率研究，考慮到干旱發(fā)生過程緩慢的特點(diǎn)，在此次應(yīng)用中，以旬作為分析時段。

3.1 干旱事件識別

3.1.1 識別規(guī)則

根據(jù)干旱的持續(xù)性，結(jié)合干旱指標(biāo)的特點(diǎn)規(guī)定水文干旱事件的識別規(guī)則如下：

（1） 連續(xù)三個旬以上（包括三個旬）水文干旱指標(biāo)處于截取水平以下判定為發(fā)生水文干旱事件；

（2） 干旱歷時為水文干旱指標(biāo)連續(xù)處于截取水平以下旬?dāng)?shù)（大于等于三個旬）；

（3） 干旱強(qiáng)度為干旱歷時所對應(yīng)水文干旱指標(biāo)均值的絕對值；

（4） 干旱烈度為干旱歷時所對應(yīng)的水文干旱指標(biāo)之和的絕對值。

3.1.2 截取水平分析

徑流距平指標(biāo)截取水平參照水利部公益性行業(yè)專項《土壤墑情多源信息綜合與預(yù)測預(yù)警示范研究》的研究成果加以確定。Z指數(shù)與SDI指數(shù)的截取水平，依據(jù)徑流距平指標(biāo)截取水平所相應(yīng)的頻率，再據(jù)此反查Z指數(shù)與SDI指數(shù)的頻率曲線來求得。徑流距平指標(biāo)、Z指數(shù)與SDI指數(shù)的截取水平如表1所示。

3.1.3 識別結(jié)果

在已確定的識別規(guī)則和截取水平下，應(yīng)用三個水文干旱指標(biāo)對樟樹水文站1970-2009年旬徑流系列進(jìn)行水文干旱識別，并對識別的出的干旱事件的干旱歷時、干旱強(qiáng)度以及干旱烈度進(jìn)行統(tǒng)計，其結(jié)果如表2所示。

3.2 干旱重現(xiàn)期計算

基于三個水文干旱指標(biāo)識別出的水文干旱場次及對應(yīng)的水文干旱特征變量系列，利用Copula函數(shù)建立干旱歷時與干旱強(qiáng)度聯(lián)合概率分布函數(shù)，并分三個指標(biāo)對水文干旱歷時與干旱強(qiáng)度的聯(lián)合頻率，條件概率以及重現(xiàn)期進(jìn)行分析。

3.2.1 邊際分布擬合

應(yīng)用適線法對三個指標(biāo)得出的干旱歷時與干旱強(qiáng)度特征變量所采用的擬合函數(shù)進(jìn)行選擇，并應(yīng)用K-S檢驗(yàn)驗(yàn)證合理性，函數(shù)選擇結(jié)果及K-S檢驗(yàn)結(jié)果如表3.3所示：

可以看出K-S統(tǒng)計量皆小于不同顯著性水平下的，表明可以接受為其選定的擬合函數(shù)。

3.2.2 Copula函數(shù)的選擇

利用相關(guān)性指標(biāo)法計算計算布Copula函數(shù)的參數(shù)，分別應(yīng)用適線法、RMSE法、AIC法與BI對三個指標(biāo)得出的干旱歷時與干旱強(qiáng)度特征變量聯(lián)合分布所采用的Coupla函數(shù)進(jìn)行選擇，并應(yīng)用K-S檢驗(yàn)驗(yàn)證合理性。選擇結(jié)果及檢驗(yàn)結(jié)果如表4所示

3.2.2 計算結(jié)果

利用水文于旱歷時D和干旱強(qiáng)度M的邊際分布FD（d）與FM（m），以及選定的G-H Copula函數(shù)所建立D與M的聯(lián)合概率分布，分別計算出各場次水文干旱事件歷時與強(qiáng)度的同現(xiàn)重現(xiàn)期（以年為單位）。應(yīng)用各水文干旱指標(biāo)所得指標(biāo)重現(xiàn)期的等值線分別如圖2至4所示。

由于人們往往關(guān)注的是旱情嚴(yán)重、重現(xiàn)期大的干旱事件，故將采用三個指標(biāo)所得重現(xiàn)期大于20年的干旱事件摘錄如下表5，并進(jìn)行分析。

對比三個指標(biāo)所得重現(xiàn)期大于20年的干旱事件可以發(fā)現(xiàn)，（1）徑流距平與SDI指數(shù)所得出重現(xiàn)期最大的事件基本一致，且重現(xiàn)期較為接近，而Z指數(shù)相應(yīng)的結(jié)果則與前兩個指標(biāo)所得相差較大；（2）表中徑流距平指標(biāo)和Z指數(shù)所得干旱事件中有3場基本一致，但3場對應(yīng)重現(xiàn)期相差較大，Z指數(shù)與SDI指數(shù)所得干旱事件中有3場基本一致，對應(yīng)重現(xiàn)期同樣相差較大，而徑流距平與SDI指數(shù)所得干旱事件中有4場基本一致，且重現(xiàn)期相差不大。

4 結(jié)論

本文應(yīng)用三種水文干旱指標(biāo)，對樟樹水文站1970-2009年旬徑流系列進(jìn)行水文干旱識別提取特征變量，利用Copula函數(shù)構(gòu)建水文干旱歷時與強(qiáng)度特征聯(lián)合分布，進(jìn)行重現(xiàn)期計算，得出以下結(jié)論：

（1）各指標(biāo)對水文干旱的識別效果較為一致。利用徑流距平指標(biāo)、Z指數(shù)和SDI指數(shù)三個水文干旱指標(biāo)，分別識別出大小共計77、75、79次水文干旱事件。

（2）對三個指標(biāo)得出干旱歷時進(jìn)行擬合均可采用指數(shù)分布并用極大似然法進(jìn)行參數(shù)估計，對干旱強(qiáng)度的擬合不同的指標(biāo)所采用的分布則不盡相同，由徑流距平指標(biāo)、Z指數(shù)以及SDI指數(shù)得出的干旱強(qiáng)度系列分別采用正態(tài)分布、對數(shù)正態(tài)分布與gamma分布較為合適。干旱歷時與干旱強(qiáng)度聯(lián)合分布則均采用G-H Copula。

（3）徑流距平指標(biāo)與SDI指數(shù)在最大重現(xiàn)期以及20年以上重現(xiàn)期計算中結(jié)果較為一致，Z指數(shù)與其他兩個指數(shù)相應(yīng)的重現(xiàn)期計算結(jié)果則相差較大。說明徑流距平指標(biāo)與SDI指數(shù)應(yīng)用于基于Copula函數(shù)的干旱重現(xiàn)期計算效果相對較好。

【參考文獻(xiàn)】

[1]馮國章.極限水文干旱歷時概率分布的解析與模擬研究[J].地理學(xué)報.1994，49（5）：457-466.

[2]周玉良，袁瀟晨，金菊良，等.基于Copula的區(qū)域水文干旱頻率分析.地理科學(xué)，2011，31（11）：1383-1388.

[3]陸桂華，閆桂霞，吳志勇，等.基于copula函數(shù)的區(qū)域干旱分析方法[J].水科學(xué)進(jìn)展，2010，21（2）：188-193.

[4]Hisdal H，Tallaksen L M.Estimation of regional meteorological and hydrological drought characteristics：a case study for Denmark[J].Journal of Hydrology，2003，281：230-247.

[5]Konstantinos M A，Elizabeth A C，Andrew W，et al Twentieth-Century Drought in the Conterminous United States[J].Journal of Hydrometeor，2005，6 （6）：985-1001.endprint