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

    長江下游多維徑流豐枯遭遇及豐枯演變分析

    2021-12-09 02:43卞佳琪于慧李強王炎良李明生林夢然
    人民長江 2021年10期
    關鍵詞:徑流

    卞佳琪 于慧 李強 王炎良 李明生 林夢然

    摘要:九江站、湖口站和大通站是長江下游的3個重要水文站,識別三者之間的徑流關系、豐枯演變和豐枯遭遇特征,對識別長江下游多維徑流關系具有重要的意義?;诨瑒酉嚓P系數(shù)法分析了長江下游多維徑流相關性的時變特征,基于Copula理論構建了長江下游多維徑流聯(lián)合分布和徑流演變聯(lián)合分布,分析了豐枯遭遇及豐枯演變特性。結(jié)果表明:① 徑流相關性時變強度由高至低為九江-湖口站、湖口-大通站、九江-大通站;各站徑流具有較高的豐枯同步性,由高至低為九江-大通站、湖口-大通站、九江-湖口站,豐枯異步性特征反之。② 九江-湖口站、九江-大通站、湖口-大通站的豐枯同步率最高時期分別為非汛期、汛期、非汛期;各站徑流豐枯狀態(tài)的自轉(zhuǎn)和互轉(zhuǎn)概率均較穩(wěn)定,目前無極端情況出現(xiàn)。③ 九江、大通水文站的年和汛期徑流最易豐枯互轉(zhuǎn),非汛期徑流最易連豐、連枯。④ 湖口水文站的年和汛期徑流最易連豐、連枯,非汛期徑流最易豐枯互轉(zhuǎn)。研究成果可為長江下游的多維徑流關系分析提供理論參考。

    關鍵詞:徑流; 豐枯演變; 豐枯遭遇; Copula函數(shù); 多維聯(lián)合分布; 長江下游

    中圖法分類號: TV121

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2021.10.018

    0引 言

    長江流域水系繁多、水資源豐富,年際、年內(nèi)主要暴雨區(qū)位置不同,使其不同時期的徑流來源與組成有著較大區(qū)別,研究長江下游多時空維度徑流的豐枯演變與豐枯遭遇規(guī)律,對于識別其徑流關系有一定的參考價值。分析長江下游3個重要徑流控制站九江站、湖口站和大通站之間不同時間尺度徑流的水量關系、徑流演變特征和徑流豐枯遭遇特征,有利于識別長江下游多時空維度徑流的組成和演變特征,且對了解復雜水力聯(lián)系、水量交換特征下的江湖關系有著重要的理論支撐作用。在進行多維徑流的豐枯演變和豐枯遭遇分析時,傳統(tǒng)的方法為使用有限的實測徑流數(shù)據(jù),通過單純的頻率分析計算各情景的發(fā)生頻率。由于實測徑流系列的長度有限,傳統(tǒng)的統(tǒng)計方法無法較完整地描述豐枯組合類型的內(nèi)部規(guī)律性,而Copula函數(shù)具有“能將任意兩個序列的邊緣分布連接起來構建新的聯(lián)合序列,且在新的聯(lián)合序列的邊緣分布中包含了原序列的所有信息”的特點,能夠靈活地構造邊緣分布為任意分布的水文變量聯(lián)合分布,且邊緣分布與聯(lián)合相關結(jié)構能分開構建[1]。因此采用Copula函數(shù)構建多變量聯(lián)合分布從而進行頻率計算,越來越被廣泛用于頻率分析和風險分析中。

    已有相關學者采用Copula函數(shù)對長江流域不同區(qū)域的徑流豐枯遭遇及豐枯演變規(guī)律展開了相關研究。張忠波等[2]提出了三峽工程供水期相鄰時段豐枯遭遇聯(lián)系的Copula聯(lián)合分布函數(shù);邴建平[3]基于Copula函數(shù)研究了長江干流與鄱陽湖發(fā)生不同量級洪水的遭遇概率,定量評估了人類活動影響下多因素對洪水遭遇的影響;吳海鷗等[4]運用Copula函數(shù)構建了鄱陽湖水系多維徑流聯(lián)合分布模型,采用特豐、偏豐、平、偏枯、特枯的徑流劃分,定量分析了鄱陽湖5大水系的豐枯遭遇同步聯(lián)合概率的變化特征;石衛(wèi)等[5]構建了三峽工程壩址控制站與洞庭湖、鄱陽湖入湖控制站徑流序列的Copula聯(lián)合分布模型,分析了三峽工程與兩湖天然徑流量豐枯遭遇概率及其變化規(guī)律。當前,采用Copula函數(shù)進行水文分析的研究較多[6-8],且方法逐漸成熟,但關于長江下游多時間維度(全年期、汛期、非汛期)和多空間維度(不同控制站)的徑流豐枯演變和豐枯遭遇研究較少?;诖?,本文以長江下游九江站、湖口站、大通站為研究對象,基于滑動相關系數(shù)檢驗法分析了多維徑流相關關系的年際變化特征;基于Copula理論構建了長江下游多維徑流聯(lián)合分布,從時間尺度和空間尺度分析了多維徑流的豐枯遭遇概率特征;基于Copula理論構建了長江下游各控制站的徑流演變聯(lián)合分布,從時間尺度和空間尺度分析了各控制站的徑流豐枯演變概率特征,為長江下游的多維徑流關系分析提供了理論參考。

    1研究區(qū)域概況

    研究對象為長江下游區(qū)域的九江水文站、湖口水文站和大通水文站。對于長江下游徑流控制站大通水文站來說,主要以上中游來水為主(九江水文站為控制站),鄱陽湖水系來水也是其重要組成部分(湖口水文站為控制站),其中,九江水文站位于長江中游干流尾閭,湖口水文站位于鄱陽湖水系湖口水道。大通水文站位于安徽省池州市,設立于1922年10月,觀測資料包括水位、流量等,集水面積約170萬km2。該站以下主要有淮河、滁河、青弋江、水陽江、秦淮河等小支流匯入,干流區(qū)間入江流量約占大通站流量的2%~3%左右,故大通站的流量特征基本代表長江下游河段的整體來水特征。研究區(qū)域與地理位置如圖1所示。

    3結(jié)果及分析

    3.1長江下游多維徑流年際變化特征

    基于長江下游九江、湖口和大通水文站1955~2019年徑流量數(shù)據(jù)(九江站1955~1987年徑流量由大通站與湖口站擬合外延所得),采用滑動相關系數(shù)法分析九江、湖口和大通水文站之間的全年期、汛期、非汛期徑流量相關關系演變規(guī)律。分別取步長為12,14,16,18,20,22,獲得九江與湖口、九江與大通、湖口與大通水文站的各步長滑動徑流相關系數(shù)序列與平均值序列,如圖2~4所示。結(jié)果表明:

    (1) 九江與湖口水文站的全年期徑流量、汛期徑流量、非汛期徑流量滑動相關系數(shù)變化范圍分別為[0.003,0.882],[0.016,0.842],[0.154,0.804],各時間尺度徑流量的相關性年際變化均較劇烈,有較強的時變特征;九江站與湖口站的各時間尺度徑流滑動相關系數(shù)均于20世紀80年代前呈上升趨勢,80年代至90年代呈下降趨勢,90年代至今呈周期性上升趨勢,且周期較前期較短。

    (2) 九江與大通水文站的全年期徑流量、汛期徑流量、非汛期徑流量滑動相關系數(shù)變化范圍分別為[0.912,0.994],[0.827,0.983],[0.771,0.981],各時間尺度徑流量自20世紀50年代以來穩(wěn)定保持著較高的相關性,年際變化不明顯;九江站與大通站的全年期徑流、汛期徑流滑動相關系數(shù)均于80年代前呈現(xiàn)上升趨勢,80年代至90年代呈現(xiàn)下降趨勢,90年代至今呈周期性上升趨勢,且周期較前期較短。非汛期徑流滑動相關系數(shù)于1985年前呈上升趨勢,1985年至1995年呈下降趨勢,1995年至今較為穩(wěn)定。

    (3) 湖口與大通水文站的全年期徑流量、汛期徑流量、非汛期徑流量滑動相關系數(shù)變化范圍分別為[0.362,0.945],[0.282,0.923],[0.446,0.951],各時間尺度徑流量的相關系數(shù)的變化區(qū)間與強度均介于九江-湖口站和九江-大通站之間;湖口站與大通站的全年期徑流、汛期徑流、非汛期徑流滑動相關系數(shù)均于20世紀80年代前呈現(xiàn)上升趨勢,80年代至90年代呈現(xiàn)下降趨勢,90年代至今呈現(xiàn)周期性上升趨勢。

    分析可知,長江下游九江、湖口、大通水文站的年徑流在不同時期,其相關性、年際變化規(guī)律不盡相同,尤其九江與湖口水文站之間,因此,對長江下游來水系統(tǒng)的豐枯遭遇研究是必要的。

    3.2邊緣分布擬合

    以長江下游九江、湖口、大通水文站的全年期、汛期、非汛期徑流量資料為基礎進行頻率計算,擬合出的P-Ⅲ分布函數(shù)參數(shù)、擬合優(yōu)度指標(RMSE)、χ2檢驗統(tǒng)計量如表3所列。RMSE均較小,擬合效果較好,且在顯著水平0.01下,P-Ⅲ分布的χ2值均小于χ2檢驗臨界值,即通過了χ2檢驗,因此擬合的九江、湖口、大通水文站的全年期、汛期、非汛期徑流量分布函數(shù)均合理。

    3.3長江下游徑流豐枯遭遇分析

    3.3.1聯(lián)合分布模型

    采用Archimedean Copula函數(shù)分別構建長江下游多維徑流的全年期徑流量、汛期徑流量、非汛期徑流量二維聯(lián)合分布函數(shù),用RMSE最小和AIC信息準則值最小來選擇Copula函數(shù)[15],擬合出的Copula函數(shù)參數(shù)、擬合優(yōu)度指標如表4所列。結(jié)果表明:用G-H Copula函數(shù)描述九江-湖口、九江-大通、湖口-大通站的全年期、汛期徑流量聯(lián)合分布均最優(yōu),用ClaytonCopula函數(shù)描述九江-湖口、九江-大通、湖口-大通站的非汛期徑流量聯(lián)合分布均最優(yōu)。

    為檢驗Copula函數(shù)的擬合精度,繪制累積頻率與經(jīng)驗點據(jù)擬合圖進行模型檢驗,如圖5所示,可以看出,采用本文方法能較好地擬合長江下游多維徑流的聯(lián)合分布函數(shù)。當已知多個邊緣分布的隨機變量時,Copula函數(shù)是個非常好的聯(lián)合分布建模工具。

    3.3.2時間尺度豐枯遭遇特征

    長江下游多維徑流豐枯遭遇的概率如表5所列,長江下游多時空尺度徑流量遭遇的聯(lián)合分布等值線圖如圖6所示。結(jié)果表明,從時間尺度來看,長江下游全年期、汛期、非汛期徑流的豐枯遭遇特征如下:

    (1) 九江與湖口水文站在全年期同時遭遇豐、平、枯的概率分別為23.55%,7.82%,22.14%,在汛期同時遭遇豐、平、枯的概率分別為21.51%,7.26%,20.24%,在非汛期同時遭遇豐、平、枯的概率分別為22.18%,7.89%,25.17%。全年期、汛期、非汛期豐枯同步概率分別為53.51%,49.01%,55.24%,表明九江與湖口站的徑流豐枯同步率最高為非汛期,全年期次之,最低為汛期。

    (2) 九江與大通水文站在全年同時遭遇豐、平、枯的概率分別為33.46%,16.37%,32.62%,在汛期同時遭遇豐、平、枯的概率分別為33.78%,16.98%,33.00%,在非汛期同時遭遇豐、平、枯的概率分別為31.98%,16.33%,34.12%。全年期、汛期、非汛期豐枯同步概率分別為82.45%,83.76%,82.43%,表明九江與大通水文站的徑流豐枯同步率極高,最高為汛期,全年期次之,最低為非汛期。

    (3) 湖口與大通水文站在全年同時遭遇豐、平、枯的概率分別為28.14%,10.06%,26.73%,在汛期同時遭遇豐、平、枯的概率分別為25.80%,8.70%,24.34%,在非汛期同時遭遇豐、平、枯的概率分別為27.50%,11.10%,30.66%。全年期、汛期、非汛期豐枯同步概率分別為64.93%,58.84%,69.26%,表明湖口與大通水文站的徑流豐枯同步率最高為非汛期,全年期次之,最低為汛期。

    3.3.3空間尺度豐枯遭遇特征

    結(jié)合長江下游多維徑流豐枯遭遇的概率表(見表5),對長江下游多時空尺度徑流量遭遇的聯(lián)合分布等值線圖(見圖6)進行分析。結(jié)果表明,從空間尺度來看,長江下游九江、湖口、大通水文站的徑流豐枯遭遇特征如下:

    (1) 在全年期、汛期、非汛期時間尺度下,任意2個控制站同時遭遇豐/枯的概率均分別在22%,20%,22%以上,均為所有狀態(tài)概率最大的兩個值,表明長江下游各控制站的年徑流量豐枯狀態(tài)具有較高的同步性。

    (2) 在任意時間尺度下,九江與湖口、九江與大通、湖口與大通水文站徑流豐枯同步的概率分別在49%,82%,58%以上,且對于長江下游3個重要水文站兩兩之間的徑流豐枯同步情況而言,九江與大通水文站的徑流豐枯同步性>湖口與大通水文站的徑流豐枯同步性>九江與湖口水文站的徑流豐枯同步性。

    (3) 從豐枯異步的角度來分析,全年期、汛期、非汛期九江與湖口水文站豐枯異步的概率分別為12.12%,15.52%,10.54%,全年期、汛期、非汛期九江與大通水文站豐枯異步的概率分別為0.30%,0.20%,0.24%,全年期、汛期、非汛期湖口與大通水文站豐枯異步的概率分別為5.18%,8.56%,2.94%。表明在全年期、汛期、非汛期時間尺度下,九江與湖口水文站徑流豐枯異步的概率>湖口與大通水文站徑流豐枯異步的概率>九江與大通水文站徑流豐枯異步的概率,且?guī)缀醪粫霈F(xiàn)九江和大通水文站一豐一枯的情況。

    3.4長江下游徑流演變規(guī)律分析

    3.4.1聯(lián)合分布模型

    采用Archimedean Copula函數(shù)分別構建長江下游各水文控制站前一年與后一年的全年期徑流量、汛期徑流量、非汛期徑流量的聯(lián)合分布函數(shù),用RMSE最小和AIC信息準則值最小來選擇Copula函數(shù),擬合出的Copula函數(shù)參數(shù)、擬合優(yōu)度指標如表6所列。結(jié)果表明:采用AMH Copula擬合九江、湖口、大通水文站的年徑流演變聯(lián)合分布最優(yōu),擬合九江、大通水文站的汛期徑流演變聯(lián)合分布最優(yōu);采用Clayton Copula函數(shù)擬合湖口水文站的汛期和大通水文站的非汛期徑流演變聯(lián)合分布均最優(yōu);采用G-H Copula函數(shù)擬合九江水文站的非汛期徑流演變聯(lián)合分布最優(yōu);采用Frank Copula函數(shù)擬合湖口水文站的非汛期徑流演變聯(lián)合分布最優(yōu)。

    3.4.2時間尺度豐枯演變特征

    長江下游多維徑流豐枯演變的概率如表7所列。結(jié)果表明,從時間尺度來看,長江下游全年期、汛期、非汛期徑流的豐枯演變特征如下:

    (1) 就徑流豐枯狀態(tài)自轉(zhuǎn)移而言,九江站全年期、汛期、非汛期徑流豐枯狀態(tài)保持的概率分別為25.83%,26.01%,32.92%;湖口站全年期、汛期、非汛期徑流豐枯狀態(tài)保持的概率分別為29.75%,30.61%,26.86%;大通站全年期、汛期、非汛期徑流豐枯狀態(tài)保持的概率分別為25.59%,23.99%,31.99%。結(jié)果表明:九江站非汛期徑流豐枯狀態(tài)保持概率最大,汛期次之,全年期最小;湖口站汛期徑流豐枯狀態(tài)保持概率最大,全年期次之,非汛期最小;大通站非汛期徑流豐枯狀態(tài)保持概率最大,全年期次之,汛期最小。

    (2) 就徑流豐枯狀態(tài)互轉(zhuǎn)移而言,九江站全年期、汛期、非汛期徑流豐枯狀態(tài)互轉(zhuǎn)的概率分別為30.44%,30.24%,23.56%,湖口站全年期、汛期、非汛期徑流出現(xiàn)豐枯狀態(tài)互轉(zhuǎn)的概率分別為26.52%,25.72%,29.40%,大通站全年期、汛期、非汛期徑流出現(xiàn)豐枯狀態(tài)互轉(zhuǎn)的概率分別為30.68%,32.30%,24.40%。結(jié)果表明:九江站全年期徑流豐枯狀態(tài)互轉(zhuǎn)概率最大,汛期次之,非汛期最小;湖口站非汛期徑流豐枯狀態(tài)互轉(zhuǎn)概率最大,全年期次之,汛期最小;大通站汛期徑流豐枯狀態(tài)互轉(zhuǎn)概率最大,全年期次之,非汛期最小。

    3.4.3空間尺度豐枯演變特征

    由長江下游多維徑流豐枯演變的概率表(見表7)分析可知,從空間尺度來看,長江下游九江站、湖口站、大通站的徑流豐枯演變特征如下:

    (1) 在全年期、汛期、非汛期時間尺度下,九江站與大通站的豐枯演變概率一致性較高,表明九江站與大通站的徑流豐枯狀態(tài)同步性較九江站與湖口站、大通站與湖口站的同步性略高。由經(jīng)驗可知,大通站位于九江站與湖口站匯流區(qū)域的下游,其徑流近似由九江站和湖口站的徑流組成,且九江站的流量級較湖口站高,大通站的徑流演變規(guī)律受九江站徑流的影響顯著高于湖口站。故本文的計算結(jié)果滿足經(jīng)驗分析的規(guī)律。

    (2) 總體來說,長江下游各控制站不同時間尺度徑流的豐、平、枯狀態(tài)之間的轉(zhuǎn)移較為穩(wěn)定,概率在6.26%~16.76%,表明長江下游徑流的豐枯演變目前無極端異常的情況出現(xiàn)。其中,九江站全年期和汛期徑流最易出現(xiàn)豐枯互轉(zhuǎn),非汛期徑流最易出現(xiàn)連豐/連枯;湖口站全年期和汛期徑流最易出現(xiàn)連豐/連枯,非汛期徑流最易出現(xiàn)豐枯互轉(zhuǎn);大通站情況與九江站類似,全年期和汛期徑流最易出現(xiàn)豐枯互轉(zhuǎn),非汛期徑流最易出現(xiàn)連豐/連枯。

    4結(jié) 論

    本文從時空尺度定量分析了長江下游多維徑流豐枯遭遇和豐枯演變特征、識別了各情景的發(fā)生概率,主要結(jié)論如下:

    (1) 分析了九江、湖口、大通水文站的年、汛期、非汛期徑流的相關性時變特征。研究表明:九江-湖口站各尺度徑流量的相關性均有較強的時變特征,九江-大通站各尺度徑流量自20世紀50年代起穩(wěn)定保持著較高的相關性,湖口-大通站各尺度徑流量的相關性變化強度均介于九江-湖口站和九江-大通站之間。

    (2) 分析了長江下游多維徑流的豐枯遭遇特征。研究表明:從時間尺度看,九江與湖口水文站的徑流豐枯同步率為非汛期>年>汛期,九江與大通水文站為汛期>年>非汛期,湖口與大通水文站為非汛期>年>汛期;從空間尺度看,九江與大通水文站的徑流豐枯同步性>湖口與大通水文站>九江與湖口水文站,九江與湖口水文站徑流豐枯異步性>湖口與大通水文站>九江與大通水文站。

    (3) 分析了長江下游多維徑流的豐枯演變特征。研究表明:從時間尺度看,九江站徑流豐枯保持率為非汛期>汛期>年,湖口站為汛期>年>非汛期,大通站為非汛期>年>汛期,各站豐枯互轉(zhuǎn)率則相反;從空間尺度看,九江與大通站的豐枯同步性較高,各站的豐枯演變目前無極端情況出現(xiàn),概率在6.26%~16.76%。九江、大通水文站的年和汛期最易豐枯互轉(zhuǎn),非汛期最易連豐、連枯;湖口水文站年和汛期最易連豐、連枯,非汛期最易豐枯互轉(zhuǎn)。

    以上研究結(jié)果可為長江下游的多維徑流關系分析、防汛抗旱風險分析等提供理論參考。

    參考文獻:

    [1]陳璐.Copula函數(shù)理論在多變量水文分析計算中的應用研究[M].武漢:武漢大學出版社,2013.

    [2]張忠波,張雙虎,王浩.基于Copula函數(shù)的三峽工程供水期豐枯遭遇分析[J].人民長江,2012,43(3):5-8.

    [3]邴建平.長江-鄱陽湖江湖關系演變趨勢與調(diào)控效應研究[D].武漢:武漢大學,2018.

    [4]吳海鷗,涂新軍,杜奕良,等.基于Copula函數(shù)的鄱陽湖水系徑流豐枯遭遇多維分析[J].湖泊科學,2019,31(3):801-813.

    [5]石衛(wèi),雷靜,李書飛,等.三峽工程與兩湖河川徑流豐枯遭遇研究[J].人民長江,2019,50(8):91-96.

    [6]謝華,羅強,黃介生.基于三維copula函數(shù)的多水文區(qū)豐枯遭遇分析[J].水科學進展,2012,23(2):186-193.

    [7]TU X,SINGH V P,CHEN X,et al.Uncertainty and variability in bivariate modeling of hydrological droughts[J].Stochastic Environmental Research and Risk Assessment,2016,30(5):1317-1334.

    [8]TU X,DU Y L,SINGH V P,et al.Joint distribution of design precipitation and tide and impact of sampling in a coastal area[J].International Journal of Climatology,2017,38(S1):290-302.

    [9]陳廣圣,王義民,劉登峰,等.基于滑動相關系數(shù)法的流域徑流變異診斷[J].自然災害學報,2016,25(1):11-18.

    [10]SKLAR A.Fonctions de répartition àn dimensions et leursmarges[J].Publication de lInstitut de Statistique de lUniversitede Paris,1959,8:229-231.

    [11]張堯庭.連接函數(shù)(Copula)技術與金融風險分化[J].統(tǒng)計研究,2002(4):48-51.

    [12]SALVADORI G,DE MICHELE C.On the use of copulas in hydrology:theory and practice[J].Journal of Hydrologic Engineering,2007,12(4):369-380.

    [13]SCHWEIZER B.Introduction to Copulas[J].Journal ofHydrologic Engineering,2007,12(4):346.

    [14]AGHAKOUCHAK A.A multivariate approach for persistence-based drought prediction:application to the 2010-2011 east africadrought[J].Journal of Hydrology,2015,526:127-135.

    [15]張強,李劍鋒,陳曉宏,等.基于Copula 函數(shù)的新疆極端降水概率時空變化特征[J].地理學報,2011,66(1):3-12.

    (編輯:謝玲嫻)

    Abstract:Jiujiang Station,Hukou Station and Datong Station are three crucial hydrological stations in the lower reaches of the Changjiang River.Analyzing their runoff relationships,wetness-dryness evolutionand encounter characteristics is of great significance for identifying multi-dimensional runoff relationships in the lower reaches of the Changjiang River.Based on the sliding correlation coefficient method,we analyzed the time-varying characteristic of multi-dimensional runoff correlations in the lower reaches of the Changjiang River.Then we constructed a joint distribution of multi-dimensional runoff and runoff evolution in the lower reaches of the Changjiang River based on the Copula theory,and analyzed the wetness-dryness evolutionand encounter characteristics.The results show that the time-varying intensity of runoff correlation is the highest in Jiujiang-HukouStation,followed by Hukou-Datong Station,and the minimum is Jiujiang-Datong Station.The runoff at each station has a high synchronization of wetness-dryness state,the highest isJiujiang-Datong Station,followed by Hukou-Datong Station,the lowest isJiujiang-Hukou Station,while theasynchronous feature of wetness-dryness state is reverse.The periods with the highest synchronization rates of wetness-dryness state for Jiujiang-Hukou Station,Jiujiang-Datong Station and Hukou-Datong Station are non-flood period,flood period,and non-flood period,respectively.The self-change and mutual change probabilities of runoff wetness-dryness state at each station are relatively stable,and there are no extreme situations at present.The annual runoff and flood period runoff at Jiujiang and Datong Stationare most likely to be in the wetness-dryness mutual changestate,and the runoff in non-flood period is most likely to be in the continuous wetness and continuous dryness state.While the annual runoff and the flood periodrunoff of Hukou Station is most likely to be in the continuous wetness and continuous dryness state,and the runoff in non-flood period is most likely to be in the mutual wetness-drynesschange state.The research results can provide theoretical references for the analysis of multi-dimensional runoff relationships in the lower reaches of the Changjiang River.

    Key words:runoff;wetness-dryness evolution;wetness-dryness encounter;Copula function;multi-dimensional joint distribution;lower reaches of the Changjiang River

    猜你喜歡
    徑流
    西南岔河徑流特性實例分析
    汾河水庫控制流域徑流突變及其驅(qū)動因素
    圍場縣不澄河流域水資源狀況分析
    寬城縣長河流域水資源狀況分析
    1961年—2013年渭河流域降水與徑流變化特征
    雙柏縣降雨徑流特性分析
    西安市主城區(qū)土地利用變化及其對地表徑流的影響
    石灰?guī)r坡面花生不同覆蓋度水土及養(yǎng)分流失規(guī)律
    閩東北鷲峰山不同跡地與不同植被恢復模式對徑流的影響
    陷落柱的成因淺析
    91久久精品国产一区二区成人| 亚洲第一电影网av| 国产免费男女视频| 国产高潮美女av| 在线免费观看的www视频| 天天躁日日操中文字幕| 日本黄色片子视频| 午夜日韩欧美国产| 91字幕亚洲| 五月玫瑰六月丁香| 给我免费播放毛片高清在线观看| 国产aⅴ精品一区二区三区波| 白带黄色成豆腐渣| 美女免费视频网站| 搡老妇女老女人老熟妇| 男人狂女人下面高潮的视频| 在线国产一区二区在线| 国产色爽女视频免费观看| 精品国产亚洲在线| 99久国产av精品| 精品一区二区三区视频在线| 国产成人影院久久av| 免费观看的影片在线观看| 一边摸一边抽搐一进一小说| 日本 欧美在线| 在线观看午夜福利视频| 国产精华一区二区三区| 亚洲欧美日韩高清专用| 热99re8久久精品国产| 美女高潮的动态| 女人被狂操c到高潮| 国产美女午夜福利| 免费一级毛片在线播放高清视频| 婷婷精品国产亚洲av在线| 国产一区二区亚洲精品在线观看| 男人和女人高潮做爰伦理| 男女那种视频在线观看| 日本黄色视频三级网站网址| 亚洲av美国av| 简卡轻食公司| 久久久久久久久中文| 中文字幕av成人在线电影| 国产精品一区二区三区四区久久| 国产精品爽爽va在线观看网站| 一区福利在线观看| 亚洲av熟女| 99国产精品一区二区蜜桃av| 久久久久久久久久黄片| 草草在线视频免费看| 亚洲综合色惰| 欧美精品国产亚洲| a级毛片a级免费在线| 日韩欧美免费精品| 亚洲av熟女| 1024手机看黄色片| 99久久精品国产亚洲精品| 蜜桃亚洲精品一区二区三区| 级片在线观看| 亚洲av中文字字幕乱码综合| 国产精品人妻久久久久久| 91在线观看av| 亚洲精品日韩av片在线观看| 欧美zozozo另类| 99热这里只有是精品在线观看 | www日本黄色视频网| 韩国av一区二区三区四区| 9191精品国产免费久久| 亚洲成人久久性| 日本精品一区二区三区蜜桃| 欧美xxxx性猛交bbbb| 波多野结衣巨乳人妻| 亚洲专区中文字幕在线| 免费观看人在逋| 国产单亲对白刺激| 久久国产精品影院| 久久人人爽人人爽人人片va | 男女那种视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 精品欧美国产一区二区三| 欧美xxxx黑人xx丫x性爽| 嫁个100分男人电影在线观看| 国产人妻一区二区三区在| 夜夜夜夜夜久久久久| 亚洲精品乱码久久久v下载方式| 国产aⅴ精品一区二区三区波| 波多野结衣巨乳人妻| 亚洲av五月六月丁香网| 不卡一级毛片| 美女免费视频网站| 久久99热这里只有精品18| 九色成人免费人妻av| 一夜夜www| 一本一本综合久久| 真人做人爱边吃奶动态| 丁香欧美五月| 久久久久九九精品影院| a级一级毛片免费在线观看| 婷婷精品国产亚洲av在线| 18禁在线播放成人免费| 丰满人妻熟妇乱又伦精品不卡| 小说图片视频综合网站| 亚洲av美国av| xxxwww97欧美| 国语自产精品视频在线第100页| 九九久久精品国产亚洲av麻豆| 欧美成人一区二区免费高清观看| 国产精品野战在线观看| www日本黄色视频网| 亚洲五月天丁香| 男女那种视频在线观看| 中文字幕av成人在线电影| 99国产综合亚洲精品| 一级av片app| 婷婷亚洲欧美| 有码 亚洲区| 很黄的视频免费| 99国产精品一区二区三区| 少妇丰满av| 久久久久国内视频| 免费av观看视频| av国产免费在线观看| 欧美黑人巨大hd| 亚洲人成网站高清观看| 日韩欧美 国产精品| 我要看日韩黄色一级片| 亚洲av不卡在线观看| 黄色日韩在线| 亚洲最大成人av| 成年女人毛片免费观看观看9| h日本视频在线播放| 久久中文看片网| 欧美日韩福利视频一区二区| 国内精品一区二区在线观看| 99热只有精品国产| 91字幕亚洲| 丁香欧美五月| 色综合站精品国产| 色综合欧美亚洲国产小说| 看十八女毛片水多多多| 日本精品一区二区三区蜜桃| 精品人妻熟女av久视频| 亚洲色图av天堂| 99热只有精品国产| 免费看日本二区| 婷婷精品国产亚洲av| 欧美精品国产亚洲| 成人亚洲精品av一区二区| 小说图片视频综合网站| 麻豆国产97在线/欧美| 国产亚洲欧美98| 国产毛片a区久久久久| 丰满人妻一区二区三区视频av| 久99久视频精品免费| 久久精品国产亚洲av涩爱 | 欧美色欧美亚洲另类二区| 免费看光身美女| 国产91精品成人一区二区三区| 久99久视频精品免费| 老司机深夜福利视频在线观看| 久久久成人免费电影| 亚洲国产色片| 亚洲av五月六月丁香网| 九九热线精品视视频播放| 一级av片app| 亚洲人成网站在线播放欧美日韩| 国产精品日韩av在线免费观看| 久久草成人影院| 国语自产精品视频在线第100页| 国产高清视频在线观看网站| 男人狂女人下面高潮的视频| 久久香蕉精品热| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 永久网站在线| 免费人成视频x8x8入口观看| 中文字幕高清在线视频| 美女黄网站色视频| 国产v大片淫在线免费观看| 中文亚洲av片在线观看爽| 女生性感内裤真人,穿戴方法视频| 欧美中文日本在线观看视频| 亚洲成人久久性| 天美传媒精品一区二区| 成人美女网站在线观看视频| 最新中文字幕久久久久| 久久久久性生活片| 亚洲第一区二区三区不卡| 国产一区二区三区在线臀色熟女| 亚洲精品456在线播放app | 午夜福利18| 国产精品久久视频播放| 亚洲无线观看免费| 一二三四社区在线视频社区8| 又爽又黄无遮挡网站| 日韩中文字幕欧美一区二区| 最近视频中文字幕2019在线8| 久久精品国产亚洲av天美| 成人av一区二区三区在线看| 久久久久国内视频| 成人永久免费在线观看视频| 国产一区二区在线av高清观看| 亚洲av不卡在线观看| 一区二区三区高清视频在线| 99riav亚洲国产免费| 欧美高清成人免费视频www| 国产高清激情床上av| 国产美女午夜福利| 欧美三级亚洲精品| 国产精品自产拍在线观看55亚洲| 国产成人av教育| 国产日本99.免费观看| 国产一区二区激情短视频| 亚洲av成人av| 毛片一级片免费看久久久久 | 又爽又黄无遮挡网站| 精品不卡国产一区二区三区| 男女之事视频高清在线观看| 亚洲中文日韩欧美视频| 国产精品一区二区三区四区久久| 丰满人妻熟妇乱又伦精品不卡| 好男人电影高清在线观看| 啪啪无遮挡十八禁网站| 精品久久久久久成人av| 欧美激情国产日韩精品一区| 亚洲国产精品合色在线| 九九在线视频观看精品| 午夜福利欧美成人| 简卡轻食公司| 婷婷精品国产亚洲av| or卡值多少钱| 中文字幕久久专区| 热99在线观看视频| xxxwww97欧美| 亚洲专区中文字幕在线| 亚洲欧美清纯卡通| 国产伦一二天堂av在线观看| 男女床上黄色一级片免费看| 老司机午夜福利在线观看视频| 精品一区二区三区人妻视频| 一区二区三区免费毛片| 中文资源天堂在线| 成人特级av手机在线观看| av黄色大香蕉| 在线观看一区二区三区| 国产亚洲欧美98| 色哟哟哟哟哟哟| 男女视频在线观看网站免费| 天天一区二区日本电影三级| 噜噜噜噜噜久久久久久91| 中文字幕人妻熟人妻熟丝袜美| 国产三级在线视频| 久久久国产成人免费| 少妇高潮的动态图| 极品教师在线免费播放| 欧美潮喷喷水| 欧美黄色片欧美黄色片| 精品久久久久久成人av| 少妇人妻一区二区三区视频| x7x7x7水蜜桃| 国产高潮美女av| 色精品久久人妻99蜜桃| 天天一区二区日本电影三级| 18禁在线播放成人免费| 看黄色毛片网站| 91麻豆精品激情在线观看国产| 男人狂女人下面高潮的视频| 精品人妻1区二区| 桃红色精品国产亚洲av| 国内精品一区二区在线观看| 我要搜黄色片| 国产精品不卡视频一区二区 | 国产成人a区在线观看| 一区福利在线观看| 精品日产1卡2卡| 99久久久亚洲精品蜜臀av| 欧美黑人巨大hd| 成人永久免费在线观看视频| 国产精品精品国产色婷婷| 精品久久久久久久久久免费视频| 丰满的人妻完整版| 亚洲精品乱码久久久v下载方式| 国产三级在线视频| 99热只有精品国产| 日本成人三级电影网站| 午夜福利高清视频| 久久久久久久午夜电影| 亚洲精品乱码久久久v下载方式| 日本免费一区二区三区高清不卡| 亚洲欧美日韩高清专用| 一级黄色大片毛片| 精品久久久久久久末码| 黄色日韩在线| 男女视频在线观看网站免费| 757午夜福利合集在线观看| 成人特级av手机在线观看| 精品一区二区三区av网在线观看| 久久亚洲精品不卡| 亚洲美女搞黄在线观看 | 欧美激情久久久久久爽电影| 九色成人免费人妻av| 国产色爽女视频免费观看| 搡女人真爽免费视频火全软件 | 天美传媒精品一区二区| 久久久久久久久久成人| 噜噜噜噜噜久久久久久91| 欧美潮喷喷水| 国产一区二区激情短视频| 亚洲av成人精品一区久久| 婷婷六月久久综合丁香| 久久伊人香网站| 婷婷色综合大香蕉| 精品日产1卡2卡| 午夜福利18| 亚洲人成伊人成综合网2020| 午夜福利在线在线| 99久久精品国产亚洲精品| 国产爱豆传媒在线观看| 一级作爱视频免费观看| 午夜福利高清视频| 国产精品嫩草影院av在线观看 | 欧美又色又爽又黄视频| 麻豆国产av国片精品| 动漫黄色视频在线观看| 成年人黄色毛片网站| 国产精品亚洲一级av第二区| 成人无遮挡网站| 又黄又爽又免费观看的视频| 一个人免费在线观看电影| 色综合婷婷激情| 成年女人永久免费观看视频| 少妇的逼水好多| 变态另类丝袜制服| 国产精品电影一区二区三区| 国产成人福利小说| 真人做人爱边吃奶动态| or卡值多少钱| 一本综合久久免费| 综合色av麻豆| 国产熟女xx| 露出奶头的视频| 国产在线男女| 综合色av麻豆| 好看av亚洲va欧美ⅴa在| 午夜福利免费观看在线| 韩国av一区二区三区四区| 国产精品爽爽va在线观看网站| 国产激情偷乱视频一区二区| 九九热线精品视视频播放| 看片在线看免费视频| 精品久久国产蜜桃| 赤兔流量卡办理| 精品久久久久久久久av| 免费在线观看日本一区| 国产精品三级大全| 在线国产一区二区在线| 深爱激情五月婷婷| 亚洲欧美日韩东京热| 国产黄色小视频在线观看| 韩国av一区二区三区四区| 最近视频中文字幕2019在线8| 免费av观看视频| 露出奶头的视频| 在线国产一区二区在线| 日韩欧美免费精品| 久久久久国产精品人妻aⅴ院| 国产伦人伦偷精品视频| 成人高潮视频无遮挡免费网站| 日本a在线网址| 99国产精品一区二区三区| 嫩草影院新地址| 亚洲最大成人中文| 嫩草影院入口| 九九在线视频观看精品| 真人一进一出gif抽搐免费| 亚洲精品在线美女| 精品人妻熟女av久视频| av在线观看视频网站免费| 国产亚洲欧美98| 精品久久久久久,| 成年免费大片在线观看| eeuss影院久久| 精品99又大又爽又粗少妇毛片 | 午夜老司机福利剧场| 欧美日韩综合久久久久久 | 18禁黄网站禁片午夜丰满| 成人国产一区最新在线观看| 国产69精品久久久久777片| 亚洲自拍偷在线| 欧美3d第一页| 一a级毛片在线观看| 国产久久久一区二区三区| 他把我摸到了高潮在线观看| 夜夜躁狠狠躁天天躁| 床上黄色一级片| 成人精品一区二区免费| 18禁在线播放成人免费| 成人国产综合亚洲| 高清毛片免费观看视频网站| 精品一区二区三区人妻视频| 少妇人妻精品综合一区二区 | 亚洲成av人片免费观看| 日韩人妻高清精品专区| aaaaa片日本免费| 神马国产精品三级电影在线观看| 美女 人体艺术 gogo| 亚洲精品在线观看二区| 特大巨黑吊av在线直播| 男女做爰动态图高潮gif福利片| 女人被狂操c到高潮| 色哟哟哟哟哟哟| 女生性感内裤真人,穿戴方法视频| 美女高潮喷水抽搐中文字幕| 国产成人aa在线观看| 免费av观看视频| 久久精品国产自在天天线| 午夜免费激情av| 黄色一级大片看看| 中国美女看黄片| 欧美日韩综合久久久久久 | 最近中文字幕高清免费大全6 | 国产探花极品一区二区| 我要看日韩黄色一级片| 亚洲va日本ⅴa欧美va伊人久久| av视频在线观看入口| 看十八女毛片水多多多| 嫁个100分男人电影在线观看| 国产午夜福利久久久久久| 1024手机看黄色片| 久久欧美精品欧美久久欧美| 国产精品亚洲av一区麻豆| 亚洲精品日韩av片在线观看| 亚洲不卡免费看| 精品乱码久久久久久99久播| 精品久久久久久久人妻蜜臀av| 18禁在线播放成人免费| 国产黄片美女视频| 午夜精品一区二区三区免费看| 亚洲国产欧美人成| 国产探花极品一区二区| 日韩欧美精品免费久久 | 三级毛片av免费| 亚洲精品在线观看二区| 精品午夜福利视频在线观看一区| 国产主播在线观看一区二区| 真实男女啪啪啪动态图| 亚洲内射少妇av| 免费av观看视频| 天美传媒精品一区二区| 免费av不卡在线播放| 日本一二三区视频观看| 成人特级黄色片久久久久久久| 久9热在线精品视频| 每晚都被弄得嗷嗷叫到高潮| 一级黄片播放器| 大型黄色视频在线免费观看| 成年女人毛片免费观看观看9| 热99在线观看视频| 亚洲天堂国产精品一区在线| 美女高潮喷水抽搐中文字幕| 国产激情偷乱视频一区二区| 午夜a级毛片| 在线免费观看不下载黄p国产 | 在线观看一区二区三区| 午夜免费激情av| 欧美不卡视频在线免费观看| 麻豆成人午夜福利视频| 欧美+日韩+精品| 免费高清视频大片| 内地一区二区视频在线| 欧美日本亚洲视频在线播放| 能在线免费观看的黄片| 不卡一级毛片| 免费在线观看成人毛片| 少妇丰满av| 观看免费一级毛片| 色综合站精品国产| 免费看光身美女| 在线观看午夜福利视频| 成人毛片a级毛片在线播放| 两人在一起打扑克的视频| 欧美日本亚洲视频在线播放| 国产单亲对白刺激| 国内揄拍国产精品人妻在线| 久久国产精品人妻蜜桃| 九九热线精品视视频播放| 欧美日韩瑟瑟在线播放| 国产精品嫩草影院av在线观看 | 亚洲精品日韩av片在线观看| 国产又黄又爽又无遮挡在线| 久久久精品欧美日韩精品| 婷婷亚洲欧美| 午夜精品久久久久久毛片777| 日韩欧美精品v在线| 一本精品99久久精品77| 色综合站精品国产| 欧美不卡视频在线免费观看| 99在线视频只有这里精品首页| 村上凉子中文字幕在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 老司机福利观看| 国产成年人精品一区二区| 91九色精品人成在线观看| 久久午夜福利片| 大型黄色视频在线免费观看| 亚洲久久久久久中文字幕| 国产国拍精品亚洲av在线观看| 国产免费男女视频| 欧美成狂野欧美在线观看| 欧美日本亚洲视频在线播放| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 蜜桃久久精品国产亚洲av| 久久精品影院6| 日本熟妇午夜| 麻豆国产av国片精品| 91av网一区二区| 精品久久久久久久久久免费视频| 免费观看人在逋| 久久精品91蜜桃| 麻豆av噜噜一区二区三区| 国产私拍福利视频在线观看| 日日摸夜夜添夜夜添av毛片 | 国产精品一区二区免费欧美| 俄罗斯特黄特色一大片| 亚洲中文字幕日韩| 成人鲁丝片一二三区免费| 黄片小视频在线播放| 国产探花在线观看一区二区| 一夜夜www| 色哟哟哟哟哟哟| 如何舔出高潮| 亚洲内射少妇av| 欧美不卡视频在线免费观看| 一进一出抽搐动态| 久久99热这里只有精品18| av在线老鸭窝| 美女xxoo啪啪120秒动态图 | 99久久精品国产亚洲精品| 久久精品国产清高在天天线| 中文资源天堂在线| 亚洲欧美清纯卡通| 国产真实伦视频高清在线观看 | 欧美日韩综合久久久久久 | 亚洲天堂国产精品一区在线| 18禁在线播放成人免费| 成年女人看的毛片在线观看| 亚洲第一欧美日韩一区二区三区| 国产国拍精品亚洲av在线观看| 国产精品国产高清国产av| 熟女人妻精品中文字幕| 黄色女人牲交| 日日摸夜夜添夜夜添av毛片 | 日韩大尺度精品在线看网址| 欧美zozozo另类| 大型黄色视频在线免费观看| 精品久久久久久久久久久久久| 淫秽高清视频在线观看| 又紧又爽又黄一区二区| 日本黄色片子视频| 亚洲国产欧洲综合997久久,| 亚洲国产精品合色在线| 欧美xxxx性猛交bbbb| 免费av毛片视频| 欧美日韩国产亚洲二区| 亚洲欧美日韩高清专用| 国内揄拍国产精品人妻在线| 免费av观看视频| 午夜福利视频1000在线观看| av国产免费在线观看| av专区在线播放| 亚洲精华国产精华精| 免费看a级黄色片| 国产在线精品亚洲第一网站| 91在线精品国自产拍蜜月| 午夜福利高清视频| 99久久成人亚洲精品观看| 欧美一区二区亚洲| 不卡一级毛片| 特大巨黑吊av在线直播| 精品乱码久久久久久99久播| a级一级毛片免费在线观看| 美女xxoo啪啪120秒动态图 | 真人做人爱边吃奶动态| 欧美一区二区国产精品久久精品| 9191精品国产免费久久| 日本免费a在线| 国产精品国产高清国产av| www.www免费av| 国产三级在线视频| 免费在线观看日本一区| 最近最新免费中文字幕在线| av在线观看视频网站免费| 成人美女网站在线观看视频| 免费搜索国产男女视频| 亚洲真实伦在线观看| 高潮久久久久久久久久久不卡| 国产一区二区激情短视频| 人人妻人人看人人澡| 国产精品精品国产色婷婷| 国产乱人伦免费视频| 久久久久九九精品影院| 午夜福利高清视频| 又爽又黄无遮挡网站| 日韩欧美一区二区三区在线观看| 窝窝影院91人妻| 性色av乱码一区二区三区2| 黄色日韩在线| 别揉我奶头~嗯~啊~动态视频| 国产男靠女视频免费网站| 内地一区二区视频在线| 高清日韩中文字幕在线| 国产av麻豆久久久久久久| 国产伦人伦偷精品视频| av福利片在线观看| 久9热在线精品视频| 免费看a级黄色片|