• <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ī)律
    閩東北鷲峰山不同跡地與不同植被恢復模式對徑流的影響
    陷落柱的成因淺析
    免费在线观看亚洲国产| 精品欧美一区二区三区在线| 淫妇啪啪啪对白视频| 99国产精品一区二区三区| 成人国产一区最新在线观看| 99精品欧美一区二区三区四区| 久久久久久久精品吃奶| 多毛熟女@视频| 久久精品国产综合久久久| 亚洲精品中文字幕在线视频| 欧美黑人欧美精品刺激| 波多野结衣一区麻豆| 国产高清videossex| 女性被躁到高潮视频| 男女午夜视频在线观看| 免费观看人在逋| 亚洲人成伊人成综合网2020| 黄色毛片三级朝国网站| 日韩欧美国产一区二区入口| 一个人免费在线观看的高清视频| 两人在一起打扑克的视频| 日日夜夜操网爽| 看免费av毛片| 久久久久久久午夜电影| 亚洲成人国产一区在线观看| 日韩免费av在线播放| 好男人电影高清在线观看| 他把我摸到了高潮在线观看| 精品免费久久久久久久清纯| 青草久久国产| 操美女的视频在线观看| 久久精品国产亚洲av高清一级| 久久人妻福利社区极品人妻图片| a在线观看视频网站| 黄色片一级片一级黄色片| 免费久久久久久久精品成人欧美视频| 久久午夜综合久久蜜桃| 日日爽夜夜爽网站| 国产精品一区二区精品视频观看| 最好的美女福利视频网| 国产三级黄色录像| 亚洲人成电影观看| 精品国产一区二区久久| 亚洲在线自拍视频| 亚洲精品一卡2卡三卡4卡5卡| 97碰自拍视频| 亚洲片人在线观看| 搡老妇女老女人老熟妇| 国产精品秋霞免费鲁丝片| 可以免费在线观看a视频的电影网站| 看片在线看免费视频| 亚洲中文字幕日韩| 久久国产乱子伦精品免费另类| 美女高潮喷水抽搐中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 91大片在线观看| 国产伦一二天堂av在线观看| 在线观看www视频免费| 巨乳人妻的诱惑在线观看| 亚洲无线在线观看| 国产av又大| 色综合欧美亚洲国产小说| 免费在线观看日本一区| 岛国在线观看网站| 亚洲一区中文字幕在线| 国产免费男女视频| 纯流量卡能插随身wifi吗| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区综合在线观看| 国产精品久久视频播放| 成年女人毛片免费观看观看9| 亚洲情色 制服丝袜| 激情在线观看视频在线高清| 欧美激情久久久久久爽电影 | 91av网站免费观看| 老鸭窝网址在线观看| 9191精品国产免费久久| 欧美精品啪啪一区二区三区| 手机成人av网站| 欧美亚洲日本最大视频资源| 久热这里只有精品99| 久久久水蜜桃国产精品网| 久久久久久久久久久久大奶| 免费看美女性在线毛片视频| 日本a在线网址| 国产麻豆69| 看免费av毛片| 九色亚洲精品在线播放| 亚洲专区中文字幕在线| 国产真人三级小视频在线观看| 99热只有精品国产| 国产精品久久视频播放| 亚洲免费av在线视频| 亚洲精品av麻豆狂野| 手机成人av网站| 大型黄色视频在线免费观看| 国产男靠女视频免费网站| 国产蜜桃级精品一区二区三区| 午夜福利免费观看在线| 中国美女看黄片| 97超级碰碰碰精品色视频在线观看| 久久午夜综合久久蜜桃| 亚洲自偷自拍图片 自拍| 麻豆久久精品国产亚洲av| 久久狼人影院| 久久久久久大精品| 禁无遮挡网站| 法律面前人人平等表现在哪些方面| 精品久久久久久久人妻蜜臀av | 亚洲国产看品久久| 黑丝袜美女国产一区| 亚洲av第一区精品v没综合| 啦啦啦免费观看视频1| 国产精品日韩av在线免费观看 | 一本综合久久免费| 韩国精品一区二区三区| 国产91精品成人一区二区三区| 美女高潮到喷水免费观看| 亚洲av熟女| 18禁裸乳无遮挡免费网站照片 | 老鸭窝网址在线观看| aaaaa片日本免费| 黄片播放在线免费| 女性被躁到高潮视频| 日本五十路高清| 欧美日本中文国产一区发布| 欧美在线一区亚洲| 欧美黄色淫秽网站| 99国产精品一区二区三区| 久久婷婷人人爽人人干人人爱 | 韩国精品一区二区三区| av超薄肉色丝袜交足视频| 欧美大码av| 成人av一区二区三区在线看| 99国产精品一区二区蜜桃av| 两个人看的免费小视频| 国产三级在线视频| www.www免费av| 高潮久久久久久久久久久不卡| ponron亚洲| 成在线人永久免费视频| 9色porny在线观看| 色在线成人网| 给我免费播放毛片高清在线观看| netflix在线观看网站| 国产在线观看jvid| 一区在线观看完整版| 看片在线看免费视频| 免费看a级黄色片| 亚洲aⅴ乱码一区二区在线播放 | 国产精品一区二区免费欧美| 午夜影院日韩av| 69av精品久久久久久| 亚洲天堂国产精品一区在线| 男女午夜视频在线观看| 久久午夜亚洲精品久久| 两性午夜刺激爽爽歪歪视频在线观看 | 日本三级黄在线观看| 成人特级黄色片久久久久久久| videosex国产| 久久久久精品国产欧美久久久| 一进一出抽搐gif免费好疼| 国产精品精品国产色婷婷| √禁漫天堂资源中文www| 免费在线观看影片大全网站| 人成视频在线观看免费观看| 亚洲情色 制服丝袜| 波多野结衣高清无吗| 亚洲一码二码三码区别大吗| 国产1区2区3区精品| 亚洲欧美日韩无卡精品| 一区二区三区精品91| 丝袜美足系列| 久久中文看片网| 国内精品久久久久久久电影| 亚洲欧美激情在线| 日韩欧美免费精品| 50天的宝宝边吃奶边哭怎么回事| av在线播放免费不卡| 后天国语完整版免费观看| 91九色精品人成在线观看| 精品人妻在线不人妻| 国产熟女午夜一区二区三区| 欧美精品啪啪一区二区三区| 久久精品人人爽人人爽视色| 亚洲国产精品久久男人天堂| 成人av一区二区三区在线看| 国产精品久久久久久人妻精品电影| 91老司机精品| 美女大奶头视频| 久久九九热精品免费| 亚洲一区二区三区不卡视频| 变态另类丝袜制服| 日本在线视频免费播放| 少妇粗大呻吟视频| 一级片免费观看大全| 真人做人爱边吃奶动态| 久热这里只有精品99| 18禁黄网站禁片午夜丰满| 亚洲一区中文字幕在线| 国产精品自产拍在线观看55亚洲| 50天的宝宝边吃奶边哭怎么回事| 久久婷婷人人爽人人干人人爱 | 精品一品国产午夜福利视频| 窝窝影院91人妻| 国产真人三级小视频在线观看| 久久午夜亚洲精品久久| 九色亚洲精品在线播放| 国产亚洲精品久久久久5区| 久久久久久久久免费视频了| 午夜a级毛片| 国产极品粉嫩免费观看在线| 日本黄色视频三级网站网址| 国产成人精品久久二区二区91| 亚洲 欧美一区二区三区| 国产精品香港三级国产av潘金莲| 国产成人av教育| 国产精品一区二区精品视频观看| 制服丝袜大香蕉在线| 国产片内射在线| 少妇粗大呻吟视频| 午夜福利在线观看吧| 男男h啪啪无遮挡| 精品一品国产午夜福利视频| 热99re8久久精品国产| 一级黄色大片毛片| 18禁美女被吸乳视频| 黄片播放在线免费| 丰满的人妻完整版| 久久精品国产99精品国产亚洲性色 | 亚洲第一青青草原| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美日韩一区二区三区在线| 欧美日韩乱码在线| 在线国产一区二区在线| 亚洲精品美女久久久久99蜜臀| 成熟少妇高潮喷水视频| 纯流量卡能插随身wifi吗| 亚洲国产欧美一区二区综合| 看黄色毛片网站| 国产精品免费一区二区三区在线| 婷婷丁香在线五月| 色在线成人网| tocl精华| 黑人巨大精品欧美一区二区mp4| 精品久久久久久久毛片微露脸| 欧美大码av| 国产精品影院久久| 久久香蕉激情| 精品少妇一区二区三区视频日本电影| 欧美在线一区亚洲| 看黄色毛片网站| 亚洲av美国av| av天堂久久9| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品秋霞免费鲁丝片| 九色亚洲精品在线播放| 亚洲精品美女久久av网站| 老熟妇乱子伦视频在线观看| 久久久久国内视频| 午夜福利成人在线免费观看| 在线十欧美十亚洲十日本专区| 人人妻人人澡欧美一区二区 | 男人操女人黄网站| 一级毛片女人18水好多| 国产精品免费视频内射| 国产欧美日韩一区二区三| 久久人妻熟女aⅴ| 国产欧美日韩综合在线一区二区| √禁漫天堂资源中文www| 欧美乱码精品一区二区三区| 欧美日韩黄片免| 午夜免费观看网址| a在线观看视频网站| 国产精品一区二区在线不卡| 精品国产超薄肉色丝袜足j| 岛国视频午夜一区免费看| 色精品久久人妻99蜜桃| 我的亚洲天堂| 成人国产一区最新在线观看| 亚洲九九香蕉| 国产精品久久电影中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 国产精品,欧美在线| 两性夫妻黄色片| 一区二区三区国产精品乱码| 人人妻人人澡欧美一区二区 | 欧美乱色亚洲激情| 一区二区日韩欧美中文字幕| 亚洲色图 男人天堂 中文字幕| 波多野结衣巨乳人妻| 乱人伦中国视频| 两性夫妻黄色片| 亚洲专区国产一区二区| 电影成人av| 久久精品亚洲精品国产色婷小说| 国产精品av久久久久免费| 免费观看精品视频网站| 亚洲久久久国产精品| 色尼玛亚洲综合影院| 国产成人系列免费观看| 欧美中文日本在线观看视频| 久久国产乱子伦精品免费另类| 国产成人av教育| 国产成人欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产欧美日韩在线播放| 久久性视频一级片| 国产99白浆流出| 1024视频免费在线观看| 国产乱人伦免费视频| 搡老妇女老女人老熟妇| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品粉嫩美女一区| 亚洲精品av麻豆狂野| 国产精品二区激情视频| 黄色丝袜av网址大全| av在线播放免费不卡| 亚洲成国产人片在线观看| 波多野结衣高清无吗| 中文字幕人妻熟女乱码| 久久久久久人人人人人| 国产极品粉嫩免费观看在线| 伊人久久大香线蕉亚洲五| 国产av又大| 日本黄色视频三级网站网址| 欧美乱妇无乱码| 一二三四在线观看免费中文在| 久久久久久久精品吃奶| 成年版毛片免费区| 国产精品久久视频播放| 国产视频一区二区在线看| 人成视频在线观看免费观看| 欧美精品亚洲一区二区| 中国美女看黄片| 老司机在亚洲福利影院| a级毛片在线看网站| 色婷婷久久久亚洲欧美| 亚洲av熟女| 少妇熟女aⅴ在线视频| 久久人妻福利社区极品人妻图片| 熟妇人妻久久中文字幕3abv| 精品免费久久久久久久清纯| 日韩精品中文字幕看吧| 国产精品久久久久久人妻精品电影| 欧美国产精品va在线观看不卡| videosex国产| 咕卡用的链子| 国产精品精品国产色婷婷| 亚洲午夜理论影院| 免费高清在线观看日韩| 成人永久免费在线观看视频| 一夜夜www| 中文字幕精品免费在线观看视频| 狠狠狠狠99中文字幕| 搞女人的毛片| 日韩av在线大香蕉| 一级a爱视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放 | 在线观看舔阴道视频| 亚洲狠狠婷婷综合久久图片| 亚洲一区二区三区不卡视频| 日日干狠狠操夜夜爽| 欧美精品啪啪一区二区三区| 麻豆成人av在线观看| 18美女黄网站色大片免费观看| 久热爱精品视频在线9| 国产精品久久视频播放| 亚洲三区欧美一区| 亚洲色图 男人天堂 中文字幕| 国产一级毛片七仙女欲春2 | 国产1区2区3区精品| 国产一区二区三区综合在线观看| 男男h啪啪无遮挡| 露出奶头的视频| 久久久久久久久中文| videosex国产| 可以在线观看毛片的网站| 国产又爽黄色视频| a在线观看视频网站| 国产精品一区二区在线不卡| 亚洲精品中文字幕在线视频| 久久精品国产清高在天天线| 欧美人与性动交α欧美精品济南到| 在线观看舔阴道视频| 精品人妻1区二区| 久久人妻熟女aⅴ| 国产成人啪精品午夜网站| 成人国产综合亚洲| 国产一卡二卡三卡精品| 久久精品成人免费网站| 国产激情久久老熟女| 妹子高潮喷水视频| 91成年电影在线观看| 很黄的视频免费| 成人国产综合亚洲| 日本精品一区二区三区蜜桃| 在线av久久热| 亚洲成人免费电影在线观看| 51午夜福利影视在线观看| 丁香六月欧美| 在线永久观看黄色视频| av片东京热男人的天堂| 亚洲精品在线美女| 国产精品久久久久久人妻精品电影| 丝袜美足系列| 在线免费观看的www视频| 色综合欧美亚洲国产小说| 波多野结衣高清无吗| 国产精华一区二区三区| 亚洲一区中文字幕在线| 午夜福利免费观看在线| 精品国产乱子伦一区二区三区| 国产主播在线观看一区二区| 亚洲成国产人片在线观看| 亚洲 国产 在线| 日韩一卡2卡3卡4卡2021年| 最好的美女福利视频网| 欧美成人免费av一区二区三区| 亚洲片人在线观看| 人人澡人人妻人| 久久久水蜜桃国产精品网| 久久久久久久久中文| 日韩欧美免费精品| 亚洲av五月六月丁香网| 给我免费播放毛片高清在线观看| 男人舔女人的私密视频| 成年女人毛片免费观看观看9| 日韩精品免费视频一区二区三区| 两个人看的免费小视频| 久久久久久大精品| www日本在线高清视频| av超薄肉色丝袜交足视频| 一区二区三区激情视频| 国产精品美女特级片免费视频播放器 | 大陆偷拍与自拍| www日本在线高清视频| 久久天躁狠狠躁夜夜2o2o| 亚洲精品在线美女| 99久久国产精品久久久| 巨乳人妻的诱惑在线观看| 亚洲在线自拍视频| 精品乱码久久久久久99久播| 看片在线看免费视频| 免费看美女性在线毛片视频| 一二三四在线观看免费中文在| 青草久久国产| 午夜免费成人在线视频| 午夜久久久在线观看| 亚洲成人国产一区在线观看| 麻豆一二三区av精品| 国产av精品麻豆| 欧美绝顶高潮抽搐喷水| 久久精品91无色码中文字幕| 国产高清视频在线播放一区| 国产成人欧美在线观看| 国产私拍福利视频在线观看| 亚洲色图综合在线观看| 丝袜美腿诱惑在线| 国产精品99久久99久久久不卡| 青草久久国产| 人人妻人人澡欧美一区二区 | 宅男免费午夜| 国产99久久九九免费精品| av在线天堂中文字幕| 午夜福利影视在线免费观看| 亚洲七黄色美女视频| 国产视频一区二区在线看| 操出白浆在线播放| 中文字幕人成人乱码亚洲影| 国产日韩一区二区三区精品不卡| 亚洲av电影在线进入| 亚洲午夜精品一区,二区,三区| 禁无遮挡网站| av中文乱码字幕在线| 精品欧美一区二区三区在线| 亚洲av美国av| 搡老熟女国产l中国老女人| 老鸭窝网址在线观看| 色哟哟哟哟哟哟| 精品国产一区二区久久| 国产一区在线观看成人免费| 久久久久久久久免费视频了| 国产成人精品在线电影| 中文字幕精品免费在线观看视频| 一级片免费观看大全| 12—13女人毛片做爰片一| 十分钟在线观看高清视频www| av在线播放免费不卡| 我的亚洲天堂| 国产伦人伦偷精品视频| 国产亚洲精品久久久久久毛片| av福利片在线| 国产成人系列免费观看| 精品久久久久久久久久免费视频| 一个人免费在线观看的高清视频| 777久久人妻少妇嫩草av网站| 午夜精品久久久久久毛片777| 免费看a级黄色片| 国内久久婷婷六月综合欲色啪| АⅤ资源中文在线天堂| 搡老熟女国产l中国老女人| 18美女黄网站色大片免费观看| 午夜福利高清视频| 大码成人一级视频| 亚洲精品美女久久久久99蜜臀| 无限看片的www在线观看| 一本久久中文字幕| 欧美日本亚洲视频在线播放| 国产精品一区二区在线不卡| 成在线人永久免费视频| 在线观看免费视频日本深夜| 日本vs欧美在线观看视频| 国产主播在线观看一区二区| 免费一级毛片在线播放高清视频 | 动漫黄色视频在线观看| 首页视频小说图片口味搜索| 欧美激情 高清一区二区三区| 国产亚洲精品综合一区在线观看 | 日日干狠狠操夜夜爽| 欧美日韩黄片免| 女人被狂操c到高潮| 精品乱码久久久久久99久播| 国产精品亚洲美女久久久| 嫁个100分男人电影在线观看| 久久欧美精品欧美久久欧美| 伦理电影免费视频| av在线天堂中文字幕| 亚洲人成网站在线播放欧美日韩| av片东京热男人的天堂| 在线免费观看的www视频| 美国免费a级毛片| 男人的好看免费观看在线视频 | 欧美一区二区精品小视频在线| 成年女人毛片免费观看观看9| 亚洲成av片中文字幕在线观看| 老司机深夜福利视频在线观看| 如日韩欧美国产精品一区二区三区| 亚洲国产毛片av蜜桃av| 人妻丰满熟妇av一区二区三区| а√天堂www在线а√下载| 身体一侧抽搐| 一卡2卡三卡四卡精品乱码亚洲| 国产一区二区在线av高清观看| 亚洲色图av天堂| 欧美日韩一级在线毛片| 日韩精品免费视频一区二区三区| 国产成人免费无遮挡视频| 日本 欧美在线| 97碰自拍视频| av天堂久久9| 91av网站免费观看| 身体一侧抽搐| 免费久久久久久久精品成人欧美视频| 在线天堂中文资源库| 成年版毛片免费区| 国产麻豆69| 久久中文看片网| 欧美 亚洲 国产 日韩一| 欧美成狂野欧美在线观看| 国产又爽黄色视频| 久久 成人 亚洲| 国产精品久久视频播放| 国产av在哪里看| 丰满人妻熟妇乱又伦精品不卡| 嫩草影院精品99| 精品国产亚洲在线| 可以免费在线观看a视频的电影网站| 免费高清视频大片| 国产乱人伦免费视频| x7x7x7水蜜桃| 18禁国产床啪视频网站| 亚洲av片天天在线观看| 久久性视频一级片| 国产不卡一卡二| 午夜日韩欧美国产| 久久性视频一级片| 亚洲国产看品久久| 午夜日韩欧美国产| 午夜免费激情av| 丁香欧美五月| 真人做人爱边吃奶动态| a在线观看视频网站| 搡老岳熟女国产| 真人做人爱边吃奶动态| 高清毛片免费观看视频网站| 美女免费视频网站| 国产精品野战在线观看| 日韩高清综合在线| 免费少妇av软件| 免费一级毛片在线播放高清视频 | www.自偷自拍.com| 正在播放国产对白刺激| 精品免费久久久久久久清纯| 久久婷婷成人综合色麻豆| √禁漫天堂资源中文www| av在线天堂中文字幕| 麻豆久久精品国产亚洲av| 日本三级黄在线观看| 日韩欧美免费精品| 欧美一区二区精品小视频在线| 亚洲欧洲精品一区二区精品久久久| 黄频高清免费视频| 91字幕亚洲| xxx96com| 国产成人精品久久二区二区91| 波多野结衣高清无吗| 高潮久久久久久久久久久不卡| av超薄肉色丝袜交足视频| 久久人人97超碰香蕉20202| 九色国产91popny在线| cao死你这个sao货| 国产野战对白在线观看| 亚洲 欧美一区二区三区| 久久久久国内视频|