• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于最大熵原理的水文干旱指標(biāo)計算方法研究

    2018-07-05 08:54:00洪興駿郭生練王樂
    南水北調(diào)與水利科技 2018年2期

    洪興駿 郭生練 王樂

    摘要:確定水文變量的概率分布,是計算各類標(biāo)準(zhǔn)化干旱指標(biāo)的關(guān)鍵。提出了基于最大熵原理(POME)的月徑流分布構(gòu)建方法,采用多階矩作為求解最大熵的約束條件,以拉格朗日乘子法估計分布參數(shù),計算了不同時間尺度的標(biāo)準(zhǔn)化徑流干旱指數(shù)SDI,評估了漢江30個子流域歷史水文干旱情勢。結(jié)果表明:與Normal、Gamma、Weibull、Pearson Type Ⅲ等常用概率分布相比,POME分布模型可以最大程度地利用實測水文數(shù)據(jù)中的信息,有效擬合不同時間尺度的累積月徑流量,表現(xiàn)出良好的適用性;隨著時間尺度的增大,采用不同分布擬合同一尺度的月徑流量,差別逐漸減小,月徑流量概率分布特征趨于正態(tài)化。成果可為推求干旱指標(biāo),研究干旱特征的統(tǒng)計規(guī)律,進行干旱頻率分析等提供新的手段。

    關(guān)鍵詞:最大熵原理;水文干旱;標(biāo)準(zhǔn)化徑流干旱指數(shù);漢江流域

    中圖分類號:TV121文獻標(biāo)志碼:A文章編號:

    16721683(2018)02009307

    Abstract:

    Selecting the appropriate probability distribution function (PDF) for hydrological variables is of significant importance to calculating standardized drought indices.In this study,we used the Principle of Maximum Entropy (POME) method to model the PDFs of aggregated monthly streamflow on varying time scales for 30 subbasins of the Hanjiang River Basin.The first three original moments of the cumulative monthly streamflow data were chosen as the constraint functions for maximizing the entropy by the Lagrange Multiplier.The Streamflow Drought Index (SDI) was computed based on monthly streamflow records derived from several theoretical probability distributions such as POME, Normal,Gamma, Weibull,and Pearson Type Ⅲ.Results showed that the POMEbased PDFs could make the best use of the information from observed records while avoiding mistakenly introducing redundant information.They showed satisfying applicability.We found that the PDFs of cumulative monthly streamflow would trend towards normalization as the time scale increased.The proposed method can be a practical tool for calculating hydrological drought indices, analyzing drought characteristics,and performing drought frequency analysis.

    Key words:POME;hydrological drought;Streamflow Drought Index (SDI);Hanjiang River basin

    干旱是一種由水循環(huán)異常引起的水分持續(xù)性短缺現(xiàn)象,不僅會對生態(tài)環(huán)境造成嚴(yán)重破壞,同時對工農(nóng)業(yè)生產(chǎn)和城市供水等也會產(chǎn)生負(fù)面影響,長期制約著經(jīng)濟社會和人類文明發(fā)展。干旱指標(biāo)是研究干旱現(xiàn)象的基礎(chǔ),也是衡量干旱發(fā)生與否與量級大小最為直觀的定量表達。天然狀況下,翔實的實測水文資料是計算干旱指標(biāo),評估區(qū)域干旱演變特征的基礎(chǔ)。將實測水文要素時間序列看作服從某種理論概率分布的隨機變量,通過等概率轉(zhuǎn)換原理將其標(biāo)準(zhǔn)化,以評估其偏離正常情形的程度,是最常見的干旱指標(biāo)構(gòu)建方法。

    Mckee等采用Gamma分布來描述降水量的統(tǒng)計規(guī)律,提出了著名的標(biāo)準(zhǔn)化降水指數(shù)(Standardized Precipitation Index,SPI)[1]。此后,各類仿效SPI的標(biāo)準(zhǔn)化干旱指數(shù)如雨后春筍般被相繼提出,用來描述水文循環(huán)各個環(huán)節(jié)要素的盈缺情勢??紤]到生態(tài)系統(tǒng)與社會經(jīng)濟活動和地表-地下徑流的直接聯(lián)系和高度依賴性,眾多學(xué)者依托實測或模擬的徑流和地下水位,將SPI的計算框架移植到水文干旱研究之中,提出了許多標(biāo)準(zhǔn)化水文干旱指數(shù),如標(biāo)準(zhǔn)化徑流指數(shù)(Standardized Runoff Index,SRI)[2],徑流干旱指數(shù)(Streamflow Drought Index,SDI)[3],標(biāo)準(zhǔn)化地下水位指數(shù)(Standardized Groundwater level Index,SGI)[4]等。

    然而,由于水文氣象變量的時空異質(zhì)性,水文分布線型的選擇對于標(biāo)準(zhǔn)化干旱指數(shù)的計算精度有著重要的影響。特別是受樣本數(shù)量的限制,如果線型選擇不當(dāng),可能給干旱評估帶來較大的誤差[5]。受降水時空分布不確定性和流域調(diào)蓄的影響,時段河川徑流的概率分布相較降水往往表現(xiàn)出更為復(fù)雜的特征。VincenteSerrano等(2012)就指出由于分布線型各自的局限性,沒有任何一個概率分布可以保證對各月份的實測徑流量均有較好的擬合效果,因此他采用不同的分布擬合不同月份的徑流量,擇優(yōu)確定理論分布,雖然在一定程度上提高了擬合的精度,但計算過程頗為繁雜[6]。因而,如果能以一個相對統(tǒng)一的分布形式描述徑流的概率特征,將對快速且穩(wěn)健的水文干旱評估產(chǎn)生有益影響。

    Jaynes(1957)提出了著名的最大熵原理(Principle of Maximum Entropy,POME)[7,8],用于構(gòu)建隨機變量的概率分布。自其被引入到水文頻率分析研究后[9],POME近幾十年間在該領(lǐng)域已經(jīng)得到了迅速的發(fā)展[1013],在干旱頻率分析方面也得到了廣泛的應(yīng)用[1416],如吳孝情等(2014)[17]和Hong等(2015)[18]就通過POME構(gòu)建了降水量的概率分布,計算了SPI指標(biāo),指出其在分布擬合方面具有令人滿意的適用性。過往應(yīng)用研究主要集中在以降水量為對象的氣象干旱上,對表征河流生境及地表與地下水體間轉(zhuǎn)化關(guān)系的地表-地下徑流關(guān)注較少。本文擬以河川徑流為對象,采用POME在不預(yù)先假設(shè)徑流序列概率分布線型的前提下,根據(jù)實測徑流資料的統(tǒng)計特征信息,推求其概率分布函數(shù),用以計算標(biāo)準(zhǔn)化水文干旱指標(biāo)。

    1研究方法

    1.1標(biāo)準(zhǔn)化徑流干旱指數(shù)(SDI)

    本研究采用標(biāo)準(zhǔn)化徑流干旱指數(shù)(Streamflow Drought Index,SDI)表征一段時間內(nèi)的流域水文干旱情況[2]。SDI的計算公式如下:

    SDI=Φ-1[FX(x)][JY](1)

    式中:Φ-1(·)為求高斯分布累積概率的逆函數(shù)。X為某一時間尺度[WTB1X]k(k[WTBX]=1,3,6,12個月等)下的累積月徑流量;f(x)和FX(x)=∫x0f(x)dx分別為其概率密度和概率函數(shù)。當(dāng)SDI>0時,表示流域處于較濕潤的狀態(tài),反之當(dāng)SDI<0時,表示流域處于相對干旱狀態(tài)。

    1.2基于最大熵原理的月徑流分布計算方法

    熵理論認(rèn)為隨機變量不確定性的大小可以用概率分布函數(shù)來描述,而信息量的大小則可以用被消除不確定性的多少來表示。Shannon于1948年將玻爾茲曼熵的概念引入到信息論中,作為量度一個隨機變量不確定性或信息量的定量指標(biāo)[19]。假設(shè)隨機變量(本文中指累積月徑流量)X的概率密度函數(shù)為f(x),則其信息熵H(x)的定義如下:

    H(x)=-∫baf(x)lnf(x)dx[JY](2)

    式中:a、b分別為隨機變量的積分上、下界。

    Jaynes進一步提出了POME的概念,認(rèn)為在所有相容的分布中,如果挑選滿足一定約束條件下能使信息熵達到極大值的分布作為隨機變量的分布,就意味著對數(shù)據(jù)不足的人為添加信息最少,從而求得的分布有最小的偏差,最合乎自然。因此,可將尋找最適宜概率分布的過程轉(zhuǎn)化為求解最大熵的過程,使得導(dǎo)出的理論分布與系統(tǒng)的已知信息相一致。在水文頻率分析中,常選用樣本各階矩,作為求解最大熵的約束信息[20]。本文中,不假設(shè)各尺度累積月徑流的頻率分布線型,以保持理論分布與樣本的前三階原點矩的期望一致作為月徑流分布的約束條件,可以將推求統(tǒng)一形式的月徑流概率分布的過程歸結(jié)為求解如下數(shù)學(xué)規(guī)劃問題。

    實際上,作為一種形式靈活的分布構(gòu)建形式,通過設(shè)置不同的約束方程組,許多常用分布線型都可以通過POME導(dǎo)出。一些常用分布參數(shù)的POME估計可見文獻[23]。

    2研究區(qū)域與資料

    漢江是長江中游最大的支流,位于東經(jīng)106°15′-114°20′E,北緯30°10′-34°20′N,發(fā)源于秦嶺南麓,干流流經(jīng)陜西、湖北兩省,于武漢市注入長江,干流全長1 577 km,流域面積約159 000 km2(見圖1)。漢江流域地處北亞熱帶季風(fēng)區(qū),幅員遼闊,光、熱和水資源空間差異大。加上區(qū)內(nèi)垂直地帶性十分顯著,是我國降水變率較大、旱澇災(zāi)害多發(fā)的地區(qū)之一[24]。漢江流域在國家水資源安全戰(zhàn)略布局中的地位舉足輕重,其上游丹江口水庫是南水北調(diào)中線工程的水源地,中下游的江漢平原則是湖北省重要的經(jīng)濟走廊,因而,研究漢江流域水文干旱演變特征,對于保障長江流域水安全具有重要的現(xiàn)實意義。

    本文選取文獻[25]所采用的漢江流域內(nèi)30個子流域的長系列徑流作為研究對象。我國自20世紀(jì)90年代以來進行了大規(guī)模的流域開發(fā),可能在一定程度上破壞了天然水文序列的一致性[26],為保證資料一致性,本文徑流資料序列基本上僅覆蓋20世紀(jì)50年代至20世紀(jì)80年代,各子流域位置見圖1。徑流量數(shù)據(jù)資料由出口斷面逐日流量資料經(jīng)時段累積和單位換算而來。分別采用MannKendall[27,28]和Pettitt[29]方法對各子流域不同尺度逐月實測徑流進行趨勢和變點檢驗,采用自相關(guān)系數(shù)檢驗月徑流序列的獨立性,置信水平均取95%。經(jīng)檢驗,所采用各月徑流量數(shù)據(jù)均滿足一致性假設(shè)。

    3研究結(jié)果與分析

    3.1月徑流分布擬合效果

    采用前述的POME方法,對漢江30個子流域不同時間尺度逐月累積徑流量進行頻率分析,擬合月徑流分布。為對比所提POME分布的適用性,選取若干常用的統(tǒng)計分布函數(shù),包括兩參數(shù)Normal(NO)、Gamma(GM)、Weibull(WB)和三參數(shù)Pearson Type Ⅲ(PⅢ)[HJ1.91mm]分布,對漢江各子流域不同時間尺度月徑流進行擬合,參數(shù)估計方法也采用文獻[22]中的基于最大熵原理的拉格朗日乘子法,并與基于三階矩POME分布計算的月徑流分布進行比較,選用的尺度包括[WTB1X]k[WTBZ]=1,3,6,12個月。采用AndersonDarling(AD)檢驗方法[30]判斷所選用分布是否可以作為不同尺度累積月徑流的理論分布,置信度水平為95%。不同時間尺度下各分布通過AD檢驗的次數(shù)見表1。不同分布下每一個時間尺度最佳的擬合次數(shù)為:30個流域×12個月=360次,由表可知,在四個時間尺度中,POME分布通過檢驗的次數(shù)均為最高,與另一常用的三參數(shù)PⅢ分布效果相當(dāng)。其它分布中,NO分布在較小時間尺度上擬合效果欠佳,而指數(shù)分布族的GM分布和極值分布族的WB分布效果尚可。

    采用各流域由不同備選分布擬合的月徑流量理論分位數(shù)與實測徑流數(shù)據(jù)的經(jīng)驗分位數(shù)的確定性系數(shù)R2QQQ衡量擬合效果,統(tǒng)計箱形圖見圖2。從圖中可以看出,POME分布和PⅢ分布對于不同尺度的月徑流量均能取得較好的擬合效果,確定性系數(shù)基本都在0.97以上。由其是POME分布(即紅色的“箱子”),有最高的平均表現(xiàn)和最小的分布范圍,表現(xiàn)穩(wěn)定。正態(tài)分布(NO)對于偏態(tài)的小尺度累積徑流量擬合效果欠佳,隨著時間尺度的延長,正態(tài)分布的適線效果逐步增強。而同為兩參數(shù)的兩個偏態(tài)分布效果[CM(22]居中。綜上,對于各個累積時段,所提出的POME[CM)]

    方法構(gòu)建的理論分布均能較好地擬合相應(yīng)時段的徑流量。

    為進一步在微觀序列層面上展示所提分布對實測徑流量樣本的擬合能力,以漢江流域的挽魚溝子流域為例,繪制各時間尺度代表月份累積徑流量理論與經(jīng)驗分布擬合效果圖如圖3所示。由圖3可以看出,各時間尺度月累積徑流量的經(jīng)驗點據(jù)與理論曲線擬合良好,令人滿意。

    3.2干旱指標(biāo)計算及干旱特征分析

    由上一小節(jié)的分析結(jié)果,基于POME良好的擬合效果,可將它作為各尺度累積月徑流量的理論分

    布,計算相應(yīng)的SDI指標(biāo)。挽魚溝流域各累積尺度SDI長序列見圖4。由圖可以看出,在各個時間尺度下,挽魚溝流域在研究期內(nèi)旱澇事件交替發(fā)生。從1個月尺度的SDI序列可以看出,該流域各月徑流的年際變化較大,旱、澇等極端水文事件頻發(fā)。隨著累積時間尺度的增大,干旱事件數(shù)量減少。這主要是由于隨著納入累積的回溯月份增加,一些夾雜在連續(xù)干旱時期內(nèi)略高于零值的SDI(可能由于突然的一場降雨)被抹去,而長期水量豐沛時期中間略小于零值的SDI則被剔除,使得原本被割裂的并不獨立的兩場干旱事件得以連續(xù)。

    進一步地,按照通用的游程理論[31],以SDI=0為閾值,識別流域歷史水文干旱事件,統(tǒng)計干旱特征如表2所示。從表中可以看出,隨著累積尺度的增大,可以觀測到的旱澇事件呈現(xiàn)頻次(數(shù)量)、最大單[CM(22]月徑流虧缺量(即最小SDI值的絕對值)減小,但平

    3.3不同分布線型下水文干旱指標(biāo)計算對比分析

    從以上分析可以注意到,隨著時間尺度的增大,累積月徑流量的分布有逐漸趨于正態(tài)化的趨勢。為驗證這一假設(shè),將各流域基于正態(tài)分布與其他備選分布擬合計算的各尺度SDI序列的確定性系數(shù)R2,SDIQQ的統(tǒng)計箱形圖繪于圖5中。可以看出,對于所有理論分布,隨著徑流量累積時段的增加,由不同理論分布計算得到的SDI指標(biāo)之間的確定性系數(shù)逐漸增大,差異逐漸減小。這表明累積徑流量的分布線型隨著累積時段的增加,逐漸趨于正態(tài)化。

    同樣,以漢江流域挽魚溝子流域出口斷面1,3,6與12個月尺度POME、GM、WB與PⅢ分布與NO分布的SDI為例,繪制相關(guān)關(guān)系見圖6。從圖中不難看出,在時間尺度選得較小的情況下,基于不同分布計算得到的干旱指標(biāo)相關(guān)點據(jù)較為分散,特別是在尾部兩端,這表明此時由不同分布計算得到的SDI值差異較大,這對于穩(wěn)健的干旱評估顯然是不利的。隨著徑流累積時間尺度的增大,不同分布擬合結(jié)果的差異逐步縮小,基本上可以互相替代,這也表明徑流量的理論分布趨于正態(tài)化。這與Angelidis等[32]的研究結(jié)果一致。這一性質(zhì)可在隨機變量的獨立同分布假設(shè)適當(dāng)放寬的條件下,由大數(shù)定律解釋。

    4結(jié)論與展望

    為計算標(biāo)準(zhǔn)化徑流干旱指數(shù)SDI,需要首先確定各時間尺度累積月徑流量的理論概率分布。為使待求的理論分布在充分利用樣本信息,反映變量統(tǒng)計特征的同時,人為假定最小,提出了基于最大熵原理(POME)的累積月徑流量分布,并采用漢江流域30個子流域的歷史資料進行驗證,得到以下結(jié)論。

    (1)基于三階矩約束的最大熵原理求得的理論分布,可以有效擬合漢江各個子流域不同時間尺度的累積月徑流量,表現(xiàn)出良好的時空適應(yīng)性,為水文隨機變量的頻率分析提供了新的工具。

    (2)隨著時間尺度增大,采用不同分布擬合同一尺度的月徑流量的差距有所減少,月徑流量概率分布特征趨于正態(tài)化。這一性質(zhì)在選用累積尺度較大的徑流計算相應(yīng)的SDI值時具有一定的應(yīng)用前景,直接采用正態(tài)分布,能夠大大簡化計算過程。

    (3)作為一種通用的概率分布參數(shù)估計方法,POME可以為推求干旱指標(biāo)和干旱特征的統(tǒng)計規(guī)律,進行干旱頻率分析提供擬合手段。[HJ2.2mm]

    參考文獻(References):

    [1]MCKEE T B,DOESKEN N J,KLEIST J.The relationship of drought frequency and duration to time scales[A].In:Proceedings of the 8th Conference on Applied Climatology[C].Boston,MA:American Meteorological Society,1993:179183.

    [2]SHUKLA S,WOOD A W.Use of a standardized runoff index for characterizing hydrologic drought[J].Geophysical Research Letters,2008,35(2):L02405.DOI:10.1029/2007GL032487.

    [3]NALBANTIS I,TSAKIRIS G.Assessment of hydrological drought revisited[J].Water Resources Management,2009,23(5):881897.DOI:10.1007/s1126900893051.

    [4]BLOOMFIELD J P,MARCHANT B P.Analysis of groundwater drought building on the standardised precipitation index approach[J].Hydrology and Earth System Sciences,2013,17:47694787.DOI:10.5194/hess1747692013.

    [5]MISHRA A K,SINGH V P.A review of drought concepts[J].Journal of Hydrology,2010,391 (12):202216.DOI:10.1016/j.jhydrol.2010.07.012.

    [6]VICENTESERRANO S M,L?PEZMORENO J I,BEGUER?A S,et al.Accurate computation of a streamflow drought index[J].Journal of Hydrologic Engineering,2012,17(2):318332.DOI:10.1061/(ASCE)HE.19435584.0000433.

    [7]JAYNES E T.Information theory and statistical mechanics,Ⅰ[J].Physical Review,1957a,106(4):620630.DOI:10.1103/PhysRev.106.620.

    [8]JAYNES E T.Information theory and statistical mechanics,Ⅱ[J].Physical Review,1957b,108(2):171190.DOI:10.1103/PhysRev.108.171.

    [9]SONUGA J O.Principle of maximum entropy in hydrologic frequency analysis[J].Journal of Hydrology,1972,17(3):177191.DOI:10.1016/00221694(72)900030.

    [10][ZK(#]SINGH V P.On application of the Weibull distribution in hydrology[J].Water Resources Management,1987,1(1):3343.DOI:10.1007/BF00421796.

    [11]PAPALEXIOU S M,KOUTSOYIANNIS D.Entropy based derivation of probability distributions:A case study to daily rainfall[J].Advances in Water Resources,2012,45:5157.DOI:10.1016/j.advwatres.2011.11.007.

    [12]劉登峰,王棟,王遠(yuǎn)坤.考慮矩階數(shù)的最大熵方法及其在水文隨機模擬中的應(yīng)用[J].水力發(fā)電學(xué)報,2015,34(9):2028.(LIU D F,WANG D,WANG Y K.POME with constraints of moments and its application to hydrological stochastic simulations[J].Journal of Hydroelectric Engineering,2015,34(9):2028.(in Chinese))DOI:10.11660/slfdxb.20150903.

    [13]孔祥銘,董艷艷,李薇,等.最大熵方法在香溪河流域徑流分析中的應(yīng)用[J].水電能源科學(xué),2016,34(2):1720.(KONG X M,DONG [JP+1]Y Y,LI W,et al.Application of maximum entropy method in analysis of Xiangxihe Watershed streamflow[J].Water Resources and Power,2016,34(2):1720.(in Chinese))

    [14]張明,金菊良,王國慶,等.基于最大熵分布模擬的干旱頻率分析[J].水力發(fā)電學(xué)報,2013,32(1):101106.(ZHANG M,JIN JL,WANG G Q,et al.Drought frequency analysis using stochastic simulation with maximum entropy model[J].Journal of Hydroelectric Engineering,2013,32(1):101106.(in Chinese))

    [15]陳海濤,黃鑫,邱林,等.基于最大熵原理的區(qū)域農(nóng)業(yè)干旱度概率分布模型[J].水利學(xué)報,2013,44(2):221226.(CHEN H T,HUANG X,QIU L,et al.Probability distribution model of regional agricultural drought degree based on the maximum entropy principle[J].Journal of Hydraulic Engineering,2013,44(2):221226.(in Chinese))DOI:10.13243/j.cnki.slxb.2013.02.009.

    [16]HAO Z,SINGH V P.Entropybased method for bivariate drought analysis[J].Journal of Hydrologic Engineering,2013a,18:780786.DOI:10.1061/(ASCE)HE.19435584.0000621.

    [17]吳孝情,陳曉宏,何艷虎,等.基于最大熵理論的降雨SPI研究[J].水文,2014,34(6):1420.(WU X Q,CHEN X H,HE Y H,et al.Study on SPI index based on maximum entropy theory[J].Journal of China Hydrology,2014,34(6):1420.(in Chinese))DOI:10.3969/j.issn.10000852.2014.06.003.

    [18]HONG X,GUO S,XIONG L,et al.Spatial and temporal analysis of drought using entropybased standardized precipitation index:a case study in Poyang Lake basin,China[J].Theoretical and Applied Climatology,2015,122(34):543556.DOI:10.1007/s007040141312y.

    [19]SHANNON C E.A mathematical theory of communication[J].The Bell System Technical Journal,1948,27(3):379423,623656.DOI:10.1002/j.15387305.1948.tb01338.x.

    [20]Hao Z,Singh V P.Entropybased method for extreme rainfall [JP+1]analysis in Texas[J].Journal of Geophysical Research:Atmosphere,2013b,118:263273.DOI:10.1029/2011JD017394.

    [21]KAPUR J N.Maximum entropy models in science and engineering[M].New York:Wiley,1989.

    [22]MEAD L R,Papanicolaou N.Maximum entropy in the problem of moments[J].Journal of Mathematical Physics,1984,25(8):24042417.DOI:10.1063/1.526446.

    [23]肖可以,宋松柏.最大熵原理在水文頻率參數(shù)估計中的應(yīng)用[J].西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版),2010,38(2):197205.(XIAO K Y,SONG S B.Application of maximum entropy principle in the estimation of hydrological frequency parameter[J].Journal of Northwest A & F University (Natural Science Edition),2010,38(2):197205.(in Chinese))DOI:10.13207/j.cnki.jnwafu.2010.02.035.

    [24]陳華,郭生練,柴曉玲,等.漢江丹江口以上流域降水特征及變化趨勢分析[J].人民長江,2005,36(11):2931.(CHEN H,GUO S L,CHAI X L,et al.Characteristics and trends of precipitation in the upper Hanjiang basin[J].Yangtze River,2005,36(11):2931.(in Chinese))DOI:10.16232/j.cnki.10014179.2005.11.011.

    [25]XIONG L,GUO S.A twoparameter monthly water balance model and its application[J].Journal of Hydrology,1999,216(1):111123.DOI:10.1016/S00221694(98)002972.

    [26]XIONG L,YU K,GOTTSCHALK L.Estimation of the distribution of annual runoff from climatic variables using copulas[J].Water Resources Research,2014,50(9):71347152.DOI:10.1002/2013WR015159.

    [27]MANN H B.Nonparametric tests against trend[J].Econometrica:Journal of the Econometric Society,1945,13:245259.DOI:10.2307/1907187.

    [28]KENDALL M G.Rank correlation methods[M].London:Charles Griffin,1975.

    [29]PETTITT A N.A nonparametric approach to the changepoint problem[J].Applied Statistics,1979,28:126135.DOI:10.2307/2346729.

    [30]ANDERSON T W,DARLING D A.A test of goodness of fit[J].Journal of the American Statistical Association,1954,49(268):765769.DOI:10.2307/2281537.

    [31]YEVJEVICH V.Hydrology paper[M].USA:Colorado State University Press,1967.

    [32]ANGELIDIS P,MARIS F,KOTSOVINOS N,et al.Computation of drought index SPI with alternative distribution functions[J].Water Resources Management,2012,26(9):24532473.DOI:10.1007/s1126901200260.

    十分钟在线观看高清视频www | 免费大片黄手机在线观看| 在线观看www视频免费| 日本黄色片子视频| 日本黄色片子视频| 久久精品国产亚洲av天美| 极品少妇高潮喷水抽搐| 亚洲av在线观看美女高潮| 亚洲国产成人一精品久久久| 亚洲精品第二区| 久久 成人 亚洲| 国产精品人妻久久久影院| 免费人妻精品一区二区三区视频| 欧美日韩综合久久久久久| 赤兔流量卡办理| 老熟女久久久| 99久国产av精品国产电影| 在线观看美女被高潮喷水网站| 日本与韩国留学比较| 精品国产露脸久久av麻豆| 国产又色又爽无遮挡免| 国产在线一区二区三区精| 精品一区二区三区视频在线| 久久久久久久精品精品| 亚洲自偷自拍三级| 少妇被粗大猛烈的视频| 男的添女的下面高潮视频| 看十八女毛片水多多多| 成人午夜精彩视频在线观看| 亚洲国产成人一精品久久久| 国产精品99久久99久久久不卡 | 国产毛片在线视频| 热re99久久国产66热| 亚洲av.av天堂| 欧美亚洲 丝袜 人妻 在线| 国产高清不卡午夜福利| 蜜桃在线观看..| 欧美一级a爱片免费观看看| 特大巨黑吊av在线直播| 久久久久网色| 精品少妇内射三级| 成人影院久久| 自拍欧美九色日韩亚洲蝌蚪91 | 少妇熟女欧美另类| 一本久久精品| 色视频在线一区二区三区| 九色成人免费人妻av| 亚洲综合色惰| 91aial.com中文字幕在线观看| 久久99一区二区三区| 色网站视频免费| 欧美激情国产日韩精品一区| 久久精品熟女亚洲av麻豆精品| 高清毛片免费看| 热re99久久精品国产66热6| 国产成人精品福利久久| 欧美日本中文国产一区发布| 精品亚洲成国产av| 看免费成人av毛片| 蜜桃久久精品国产亚洲av| 免费在线观看成人毛片| 热99国产精品久久久久久7| 亚洲精品中文字幕在线视频 | 午夜福利影视在线免费观看| 夫妻午夜视频| 精品久久久久久久久亚洲| 你懂的网址亚洲精品在线观看| av专区在线播放| 亚洲国产欧美日韩在线播放 | 尾随美女入室| 日韩欧美一区视频在线观看 | 自拍偷自拍亚洲精品老妇| 99国产精品免费福利视频| 国产在线一区二区三区精| 精品一区二区三区视频在线| 亚洲精品久久午夜乱码| 一级毛片 在线播放| 亚洲精品aⅴ在线观看| 国产精品99久久久久久久久| 美女中出高潮动态图| 91精品一卡2卡3卡4卡| 人妻 亚洲 视频| 另类亚洲欧美激情| 在线观看免费视频网站a站| 精品视频人人做人人爽| 成人综合一区亚洲| 成人美女网站在线观看视频| 丁香六月天网| 久久久久久久久久久丰满| 亚洲av在线观看美女高潮| av线在线观看网站| 十八禁高潮呻吟视频 | 亚洲精品自拍成人| 男女边摸边吃奶| av网站免费在线观看视频| 欧美精品一区二区大全| 高清黄色对白视频在线免费看 | 亚洲欧美日韩卡通动漫| 黄色毛片三级朝国网站 | 熟女av电影| 国产伦在线观看视频一区| 人妻一区二区av| 精品少妇黑人巨大在线播放| 久久精品国产亚洲av涩爱| 中文资源天堂在线| 亚洲,欧美,日韩| 人妻一区二区av| 日韩,欧美,国产一区二区三区| 国产高清不卡午夜福利| 日韩一本色道免费dvd| 少妇被粗大猛烈的视频| 伦理电影免费视频| 婷婷色综合www| 国产无遮挡羞羞视频在线观看| 性色av一级| 建设人人有责人人尽责人人享有的| 亚洲欧美日韩另类电影网站| 国产精品一区www在线观看| 欧美 亚洲 国产 日韩一| 一本大道久久a久久精品| 久久国内精品自在自线图片| 亚洲精品国产av成人精品| 国产免费一区二区三区四区乱码| 亚洲欧美日韩卡通动漫| 日日爽夜夜爽网站| 午夜福利视频精品| 大片免费播放器 马上看| 美女福利国产在线| 亚洲国产最新在线播放| videossex国产| 亚洲综合精品二区| 我要看黄色一级片免费的| 一本大道久久a久久精品| 免费观看的影片在线观看| 久久国内精品自在自线图片| 亚洲av免费高清在线观看| 高清欧美精品videossex| 99久国产av精品国产电影| 国产老妇伦熟女老妇高清| 男人爽女人下面视频在线观看| 五月天丁香电影| 国产精品久久久久久精品电影小说| 亚洲欧美清纯卡通| 观看美女的网站| 久久精品久久精品一区二区三区| 亚洲,一卡二卡三卡| 久久6这里有精品| 毛片一级片免费看久久久久| 精品少妇内射三级| 国产精品久久久久成人av| 丝袜脚勾引网站| 国产中年淑女户外野战色| 女的被弄到高潮叫床怎么办| 麻豆成人午夜福利视频| 国产精品久久久久成人av| 美女内射精品一级片tv| 观看美女的网站| 亚洲av电影在线观看一区二区三区| 国产国拍精品亚洲av在线观看| 草草在线视频免费看| 国产 一区精品| 久久久午夜欧美精品| 久久免费观看电影| 三级国产精品欧美在线观看| 日日啪夜夜撸| 亚洲精品,欧美精品| 国产成人免费观看mmmm| 亚洲欧洲精品一区二区精品久久久 | 乱人伦中国视频| 亚洲第一av免费看| 另类亚洲欧美激情| 性高湖久久久久久久久免费观看| 永久网站在线| 免费人妻精品一区二区三区视频| 国内揄拍国产精品人妻在线| 国产一区亚洲一区在线观看| 搡老乐熟女国产| 91精品伊人久久大香线蕉| 久久国内精品自在自线图片| 国产高清有码在线观看视频| 人妻人人澡人人爽人人| 曰老女人黄片| 婷婷色av中文字幕| 乱码一卡2卡4卡精品| 两个人免费观看高清视频 | 在线观看三级黄色| 国产高清有码在线观看视频| 午夜激情福利司机影院| 日本欧美视频一区| 你懂的网址亚洲精品在线观看| 又大又黄又爽视频免费| 亚洲精品久久午夜乱码| 国产日韩欧美视频二区| av在线老鸭窝| 99久久精品一区二区三区| 日本黄色日本黄色录像| 一级毛片黄色毛片免费观看视频| 天天躁夜夜躁狠狠久久av| 国产 一区精品| 色吧在线观看| 我要看日韩黄色一级片| 久久韩国三级中文字幕| 一本一本综合久久| 国产69精品久久久久777片| 国产深夜福利视频在线观看| 日韩精品免费视频一区二区三区 | 亚洲精品视频女| 日本色播在线视频| 日韩三级伦理在线观看| 国产高清国产精品国产三级| 91精品国产国语对白视频| 大香蕉97超碰在线| 久久99热6这里只有精品| 亚洲av在线观看美女高潮| 久久这里有精品视频免费| 国产一区二区三区综合在线观看 | 一区二区三区乱码不卡18| 少妇人妻精品综合一区二区| 自线自在国产av| 国产av一区二区精品久久| 女性生殖器流出的白浆| 日韩欧美 国产精品| 哪个播放器可以免费观看大片| 成人国产av品久久久| 久久久久久久久大av| 最近最新中文字幕免费大全7| 最近最新中文字幕免费大全7| 亚洲国产欧美在线一区| 日本猛色少妇xxxxx猛交久久| 成人影院久久| 欧美性感艳星| 日韩成人av中文字幕在线观看| 一二三四中文在线观看免费高清| 嫩草影院入口| 老司机亚洲免费影院| 免费黄色在线免费观看| 性色av一级| 国产精品人妻久久久久久| 免费观看在线日韩| 精品久久久久久电影网| 精品少妇久久久久久888优播| 老司机影院成人| 午夜精品国产一区二区电影| 少妇被粗大的猛进出69影院 | 日本色播在线视频| 日韩av免费高清视频| 日韩大片免费观看网站| 天美传媒精品一区二区| 我要看黄色一级片免费的| 欧美日韩亚洲高清精品| 欧美变态另类bdsm刘玥| 天天躁夜夜躁狠狠久久av| 亚洲欧洲日产国产| 少妇 在线观看| av黄色大香蕉| 久久久久国产网址| 蜜臀久久99精品久久宅男| 18禁裸乳无遮挡动漫免费视频| 国产av国产精品国产| 免费在线观看成人毛片| 日本欧美视频一区| 欧美区成人在线视频| 一区在线观看完整版| 黄色怎么调成土黄色| 日韩 亚洲 欧美在线| 一个人免费看片子| 久久久国产精品麻豆| 日韩伦理黄色片| 91精品国产国语对白视频| 乱人伦中国视频| 人妻夜夜爽99麻豆av| 精品一区在线观看国产| 热re99久久精品国产66热6| 人妻制服诱惑在线中文字幕| 久久精品久久久久久噜噜老黄| 日日爽夜夜爽网站| 99热这里只有精品一区| 欧美区成人在线视频| xxx大片免费视频| 曰老女人黄片| 日韩电影二区| 亚洲精品视频女| 久久99热这里只频精品6学生| 黑人巨大精品欧美一区二区蜜桃 | 精品国产国语对白av| 精品亚洲成a人片在线观看| av.在线天堂| 九色成人免费人妻av| 在线免费观看不下载黄p国产| 国产一区有黄有色的免费视频| www.av在线官网国产| 国产亚洲5aaaaa淫片| 精品久久久久久久久av| 亚洲精品国产av成人精品| 熟妇人妻不卡中文字幕| 国产伦在线观看视频一区| 22中文网久久字幕| 亚洲av成人精品一二三区| 美女xxoo啪啪120秒动态图| 熟女人妻精品中文字幕| 亚洲国产毛片av蜜桃av| 女性生殖器流出的白浆| 国产 一区精品| 青春草亚洲视频在线观看| 熟女人妻精品中文字幕| 免费少妇av软件| 你懂的网址亚洲精品在线观看| 日韩精品有码人妻一区| 偷拍熟女少妇极品色| 国产中年淑女户外野战色| 精品一品国产午夜福利视频| 日本-黄色视频高清免费观看| 桃花免费在线播放| 大片免费播放器 马上看| 久久久久久久大尺度免费视频| 午夜福利,免费看| 男女啪啪激烈高潮av片| 啦啦啦在线观看免费高清www| 亚洲国产最新在线播放| 亚洲美女搞黄在线观看| 人妻制服诱惑在线中文字幕| av播播在线观看一区| 综合色丁香网| 欧美丝袜亚洲另类| a级片在线免费高清观看视频| 三级经典国产精品| av免费观看日本| 久久综合国产亚洲精品| 免费在线观看成人毛片| 国产又色又爽无遮挡免| 日韩av不卡免费在线播放| 大片电影免费在线观看免费| 国产精品嫩草影院av在线观看| 中文在线观看免费www的网站| 久久久久久久久久人人人人人人| 国内揄拍国产精品人妻在线| 卡戴珊不雅视频在线播放| 久久精品国产鲁丝片午夜精品| 欧美日本中文国产一区发布| 如日韩欧美国产精品一区二区三区 | 狠狠精品人妻久久久久久综合| 99国产精品免费福利视频| 国产精品偷伦视频观看了| 久久国内精品自在自线图片| 日日啪夜夜爽| 久久女婷五月综合色啪小说| 久久精品久久久久久噜噜老黄| 中文字幕av电影在线播放| 久久久久久久国产电影| 日韩,欧美,国产一区二区三区| 看非洲黑人一级黄片| 91精品国产国语对白视频| 99热这里只有精品一区| 色94色欧美一区二区| 精品一区二区免费观看| 性色avwww在线观看| 91成人精品电影| 寂寞人妻少妇视频99o| 高清av免费在线| 少妇 在线观看| 欧美区成人在线视频| 人人澡人人妻人| 国产淫语在线视频| 色婷婷av一区二区三区视频| 天堂俺去俺来也www色官网| 99视频精品全部免费 在线| 日韩成人av中文字幕在线观看| 日本欧美视频一区| 69精品国产乱码久久久| 高清不卡的av网站| 丝瓜视频免费看黄片| 人体艺术视频欧美日本| 亚洲av男天堂| 国产精品熟女久久久久浪| av专区在线播放| 国产精品国产三级专区第一集| 观看美女的网站| 国产免费又黄又爽又色| 久久免费观看电影| 综合色丁香网| 国产伦精品一区二区三区四那| 欧美bdsm另类| 女性被躁到高潮视频| 国产色婷婷99| 国产精品偷伦视频观看了| 日本爱情动作片www.在线观看| 国产爽快片一区二区三区| 国产av码专区亚洲av| 日韩强制内射视频| 国产无遮挡羞羞视频在线观看| 国产成人精品福利久久| 亚洲综合色惰| 免费黄色在线免费观看| 国产探花极品一区二区| 又大又黄又爽视频免费| av国产久精品久网站免费入址| 精品一区二区三卡| freevideosex欧美| 成人特级av手机在线观看| 免费看日本二区| 国产精品久久久久久久电影| 亚洲欧洲日产国产| 99热这里只有是精品在线观看| 久久鲁丝午夜福利片| 亚洲国产av新网站| 日韩 亚洲 欧美在线| 一个人看视频在线观看www免费| 草草在线视频免费看| 日本欧美国产在线视频| 国产毛片在线视频| 亚洲国产av新网站| av女优亚洲男人天堂| 午夜福利,免费看| 亚洲av男天堂| 丝袜脚勾引网站| 国产在线免费精品| 中国三级夫妇交换| 国产69精品久久久久777片| 夜夜看夜夜爽夜夜摸| 中国国产av一级| 亚洲人成网站在线播| 在线观看av片永久免费下载| 最近中文字幕高清免费大全6| 在线观看人妻少妇| 亚洲av中文av极速乱| 午夜福利影视在线免费观看| 亚洲精华国产精华液的使用体验| 亚洲一级一片aⅴ在线观看| 少妇熟女欧美另类| av天堂中文字幕网| 欧美另类一区| 纯流量卡能插随身wifi吗| 成年美女黄网站色视频大全免费 | 免费不卡的大黄色大毛片视频在线观看| 亚洲欧美一区二区三区国产| 亚洲中文av在线| 2018国产大陆天天弄谢| 国产无遮挡羞羞视频在线观看| 欧美区成人在线视频| 亚洲精品国产成人久久av| 熟妇人妻不卡中文字幕| 日韩欧美 国产精品| 丰满少妇做爰视频| 国产精品99久久久久久久久| 丰满人妻一区二区三区视频av| 久久人人爽人人片av| 色94色欧美一区二区| 亚洲欧洲精品一区二区精品久久久 | 涩涩av久久男人的天堂| 久久久午夜欧美精品| 亚洲精品国产av成人精品| 性高湖久久久久久久久免费观看| 国产精品久久久久久久久免| 在线精品无人区一区二区三| 一级毛片 在线播放| 精品99又大又爽又粗少妇毛片| 国产一区二区三区av在线| 午夜激情福利司机影院| 国产白丝娇喘喷水9色精品| 亚洲精品乱久久久久久| 最近2019中文字幕mv第一页| 精品少妇黑人巨大在线播放| 国产成人91sexporn| 日本av手机在线免费观看| 欧美变态另类bdsm刘玥| 极品人妻少妇av视频| 男人爽女人下面视频在线观看| 亚洲无线观看免费| 国产精品秋霞免费鲁丝片| 精品一区二区三区视频在线| 人妻一区二区av| 色吧在线观看| 国产日韩欧美视频二区| 久久久国产精品麻豆| 黑人高潮一二区| a级毛色黄片| 国产高清不卡午夜福利| 国国产精品蜜臀av免费| 晚上一个人看的免费电影| 欧美三级亚洲精品| 久久午夜综合久久蜜桃| 国产在线一区二区三区精| 久久久久久久久久久久大奶| 亚洲精品久久午夜乱码| 国产高清三级在线| 男人舔奶头视频| 男男h啪啪无遮挡| a级一级毛片免费在线观看| 国产极品粉嫩免费观看在线 | 观看av在线不卡| 毛片一级片免费看久久久久| 精品一品国产午夜福利视频| av福利片在线| a 毛片基地| 亚洲经典国产精华液单| 欧美国产精品一级二级三级 | 国产精品偷伦视频观看了| 最后的刺客免费高清国语| 97超碰精品成人国产| 亚洲伊人久久精品综合| 亚洲精品成人av观看孕妇| 黄色日韩在线| 国产精品三级大全| 精品一品国产午夜福利视频| 欧美日韩精品成人综合77777| 在线观看免费高清a一片| 久热这里只有精品99| 国产精品久久久久久久电影| 亚洲内射少妇av| 欧美成人午夜免费资源| 亚洲精品自拍成人| 人妻人人澡人人爽人人| 久久久久久久久久成人| 国产亚洲欧美精品永久| 纯流量卡能插随身wifi吗| 国产女主播在线喷水免费视频网站| 91在线精品国自产拍蜜月| 亚洲丝袜综合中文字幕| 亚洲怡红院男人天堂| 青春草国产在线视频| 成年美女黄网站色视频大全免费 | 欧美精品国产亚洲| 亚洲人成网站在线播| 亚洲av国产av综合av卡| 亚洲av成人精品一二三区| 国产乱来视频区| 国产白丝娇喘喷水9色精品| 国产精品女同一区二区软件| 日韩视频在线欧美| 美女国产视频在线观看| 欧美日韩亚洲高清精品| 一区二区三区精品91| 国产精品一区www在线观看| 亚洲电影在线观看av| 国产精品久久久久久av不卡| 精品少妇黑人巨大在线播放| 成人影院久久| 日韩一区二区视频免费看| 爱豆传媒免费全集在线观看| av线在线观看网站| av福利片在线| 国产精品伦人一区二区| 汤姆久久久久久久影院中文字幕| 下体分泌物呈黄色| 亚洲欧美精品自产自拍| 国模一区二区三区四区视频| 少妇高潮的动态图| 最近中文字幕高清免费大全6| 欧美高清成人免费视频www| 在线观看三级黄色| 丝袜喷水一区| 极品人妻少妇av视频| 亚洲人成网站在线播| 日本色播在线视频| 青春草国产在线视频| 日韩一本色道免费dvd| 成年人午夜在线观看视频| 成人亚洲精品一区在线观看| 亚洲色图综合在线观看| 午夜福利,免费看| 亚洲欧洲精品一区二区精品久久久 | 老司机亚洲免费影院| 日本黄色日本黄色录像| 天堂8中文在线网| 人妻一区二区av| 在线观看国产h片| 国产在线视频一区二区| 91在线精品国自产拍蜜月| 一个人看视频在线观看www免费| 成人漫画全彩无遮挡| 蜜桃在线观看..| 久久人人爽av亚洲精品天堂| 欧美激情极品国产一区二区三区 | 日韩人妻高清精品专区| 在线观看www视频免费| 成人国产麻豆网| 涩涩av久久男人的天堂| 国产成人精品婷婷| 午夜免费观看性视频| 亚洲精品色激情综合| 秋霞伦理黄片| 午夜久久久在线观看| 尾随美女入室| 高清不卡的av网站| 人体艺术视频欧美日本| 国产精品一区www在线观看| 欧美亚洲 丝袜 人妻 在线| 免费看日本二区| 99热国产这里只有精品6| 狠狠精品人妻久久久久久综合| 26uuu在线亚洲综合色| 欧美 日韩 精品 国产| 亚洲国产精品专区欧美| videos熟女内射| 日韩大片免费观看网站| 又黄又爽又刺激的免费视频.| av又黄又爽大尺度在线免费看| 国产免费又黄又爽又色| 亚洲内射少妇av| 亚洲人成网站在线观看播放| 日韩强制内射视频| 亚洲性久久影院| 久久6这里有精品| av天堂中文字幕网| 一本久久精品| 国产免费一级a男人的天堂| 乱系列少妇在线播放| freevideosex欧美| 狂野欧美激情性bbbbbb| 亚洲真实伦在线观看| 国产精品久久久久久精品古装| 国产 一区精品| 国产爽快片一区二区三区| 七月丁香在线播放| 成人综合一区亚洲|