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

    基于SRI指數(shù)的嘉陵江流域1982-2020年旱澇急轉(zhuǎn)特征分析研究

    2025-04-03 00:00:00徐靖雯姬廣興張亞麗

    摘 要:以嘉陵江流域?yàn)檠芯繉ο?,基于北碚水文?982-2020年月尺度徑流數(shù)據(jù),計(jì)算不同時(shí)間尺度的標(biāo)準(zhǔn)化徑流指數(shù)(SRI),采用游程理論分析嘉陵江流域的旱澇急轉(zhuǎn)事件特征,并結(jié)合Mann-Kendall趨勢檢驗(yàn)法分析標(biāo)準(zhǔn)化徑流指數(shù)變化趨勢.研究結(jié)果表明:嘉陵江流域1982-2009年洪澇化趨勢先減弱后增強(qiáng),2010-2020年呈現(xiàn)短期增強(qiáng)趨勢后又減弱.嘉陵江流域1982-2020年的1、5、6、8、11和12月份呈現(xiàn)不明顯的干旱化減弱趨勢;7、9和10月份干旱化趨勢增強(qiáng),2、3和4月份干旱化趨勢顯著減弱.北碚水文站旱澇急轉(zhuǎn)事件共發(fā)生39次,旱轉(zhuǎn)澇事件急轉(zhuǎn)點(diǎn)集中在4月份,對應(yīng)旱澇事件的強(qiáng)度也較高.旱澇急轉(zhuǎn)事件近40年發(fā)生頻率逐步升高,并呈現(xiàn)逐步加重的趨勢,應(yīng)引起足夠的重視.

    關(guān)鍵詞:標(biāo)準(zhǔn)化徑流指數(shù);嘉陵江流域;旱澇急轉(zhuǎn);Mann-Kendall檢驗(yàn)法

    中圖分類號(hào):P333;TV121 文獻(xiàn)標(biāo)志碼:A文章編號(hào):1000-2367(2025)03-0042-08

    近年來,由于全球氣候變化和人類活動(dòng)的影響,極端天氣增多,干旱與洪澇等重大災(zāi)害頻發(fā),造成了極其嚴(yán)重的損失[1-2.數(shù)據(jù)顯示,1950-2000年的5項(xiàng)氣候?yàn)?zāi)害(干旱、洪澇、臺(tái)風(fēng)、凍害、干熱風(fēng))中,干旱災(zāi)害發(fā)生頻次為各項(xiàng)災(zāi)害之首,洪澇災(zāi)害僅居其次[3.我國每年干旱和洪澇造成農(nóng)作物受災(zāi)面積約占農(nóng)作物總受災(zāi)面積的77%,平均2~3年我國就會(huì)遭遇一次嚴(yán)重的干旱災(zāi)害[4.研究表明,旱澇急轉(zhuǎn)(drought-flood abrupt alternation,DFAA)事件[5對自然環(huán)境及人類社會(huì)造成的影響遠(yuǎn)大于單一的干旱或者洪水事件6-7.

    從20世紀(jì)90年代開始,旱澇異常的研究開始成為國內(nèi)外廣泛關(guān)注的熱點(diǎn)問題,而旱澇急轉(zhuǎn)現(xiàn)象是旱澇異常的典型代表[8,在中國華南、長江流域及西南等地區(qū)時(shí)常發(fā)生[9-10.楊家偉等[6采用標(biāo)準(zhǔn)化加權(quán)平均降雨(standardized weighted average precipitation,SWAP)指數(shù)對長江流域旱澇急轉(zhuǎn)時(shí)空規(guī)律進(jìn)行研究,結(jié)果表明長江流域旱澇急轉(zhuǎn)事件發(fā)生頻率呈不同程度的上升趨勢;閃麗潔等[11通過長期徑流旱澇急轉(zhuǎn)指數(shù)(long term runoff drought-flood abrupt alternation index,LRDFAI)和短期徑流旱澇急轉(zhuǎn)指數(shù)(short term runoff drought-flood abrupt alternation index,SRDFI)分析長江中下游旱澇急轉(zhuǎn)時(shí)空演變特征,發(fā)現(xiàn)長江中下游旱澇急轉(zhuǎn)主要表現(xiàn)為以澇轉(zhuǎn)旱事件為主;文獻(xiàn)[12-14]分別對漢江流域、洞庭湖流域和鄱陽湖流域開展了旱澇急轉(zhuǎn)現(xiàn)象的趨勢變化特征研究.然而,對于長江流域旱澇演變研究還主要集中在長江中下游地區(qū),關(guān)于長江上游流域旱澇演變特征的研究較少.嘉陵江是長江水系中流域面積最大的支流,流域內(nèi)水資源與氣象要素的時(shí)空變化對下游三峽庫區(qū)水資源配置和持續(xù)利用具有重要的作用.此外,嘉陵江水量表現(xiàn)出夏季水量大冬季水量少的特點(diǎn),由于降水集中,經(jīng)常發(fā)生洪澇和干旱災(zāi)害,歷史上嘉陵江流域曾發(fā)生過多次洪水和斷流事件[15.旱澇事件的發(fā)生不僅會(huì)使得嘉陵江農(nóng)業(yè)供水緊張,甚至工業(yè)及城市用水都會(huì)受到重大影響[16.因此,開展嘉陵江流域旱澇演變特征研究具有非常重要的意義.

    水文干旱指數(shù)是研究旱澇變化特征的基礎(chǔ),也是科學(xué)評估風(fēng)險(xiǎn)的重要依據(jù)[17.標(biāo)準(zhǔn)化徑流指數(shù)(standardized runoff index,SRI)[9適用于多時(shí)間尺度分析(1、3、6、12月等),以及資料缺乏、地形復(fù)雜的區(qū)域,能夠綜合反映水文和氣象過程.目前在國內(nèi)外水文旱澇識(shí)別研究中應(yīng)用較多[18-20,而從SRI的角度來識(shí)別長江流域旱澇急轉(zhuǎn)特征的研究相對比較缺乏.鑒于此,本文選取嘉陵江為研究對象,基于SRI深入分析嘉陵江流域旱澇急轉(zhuǎn)時(shí)空演變規(guī)律,從而有效地為流域水資源管理和防洪抗旱提供參考依據(jù).

    1 研究區(qū)概況

    嘉陵江是長江水系中流域面積最大的支流,達(dá)16萬km2,于重慶渝中區(qū)朝天門匯入長江主河道,被譽(yù)為重慶的“母親河”.嘉陵江的徑流變化對流域內(nèi)的自然生態(tài)環(huán)境以及沿江地區(qū)的經(jīng)濟(jì)發(fā)展有著深刻的影響.嘉陵江流域四季分明,夏季溫度在33 ℃左右,降水較為充沛,全年平均降雨量400 mm左右.嘉陵江流域受季風(fēng)氣候影響較為顯著,降水集中.此外,嘉陵江流經(jīng)四川,是川渝經(jīng)濟(jì)的重要組成部分,該流域內(nèi)主要以農(nóng)業(yè)生產(chǎn)為主.

    2 數(shù)據(jù)來源

    本文研究資料采用嘉陵江流域北碚水文站(29°50′N,106°26′E)1982-2020年月尺度徑流數(shù)據(jù),數(shù)據(jù)來源于長江水文網(wǎng)(http://www.cjh.com.cn/hlns_zyyjcg.html )的長江水文公報(bào).

    3 研究方法

    3.1 標(biāo)準(zhǔn)化徑流指數(shù)

    本文通過計(jì)算標(biāo)準(zhǔn)化徑流指數(shù)值進(jìn)行嘉陵江流域旱澇識(shí)別,MCKEE等[21提出了標(biāo)準(zhǔn)化降水指數(shù)(standardized precipitation index,SPI)的定義,即利用伽馬概率密度函數(shù)求一定時(shí)間內(nèi)降水量累計(jì)概率,SHUKLA等[22于2008年參照這一概念提出了標(biāo)準(zhǔn)化徑流指數(shù)SRI,兩個(gè)指標(biāo)的計(jì)算方法類似,反映徑流量的變化.

    SRI計(jì)算方法如下:假設(shè)某一確定時(shí)間段河流的徑流量x,則滿足伽馬分布概率密度函數(shù)為:

    f(x)=∫x0(xα-1αΓ(α))e-(x/β)dx,

    式中,α、β為參數(shù),α>0、β>0,用極大似然法對α、β進(jìn)行估算.

    一定時(shí)間尺度的徑流量x的累計(jì)率為:

    F(x)=∫x0f(x)dx.

    對各項(xiàng)的累計(jì)概率進(jìn)行正標(biāo)準(zhǔn)化得相應(yīng)的SRI:

    當(dāng)0<F(x)≤0.5時(shí),令k=ln[(1/F(x)2)],則SRI=-(k-(c0-c1k+c2k2/1+d1k+d2k2+d3k3));

    當(dāng)0.5<F(x)≤1時(shí),令k=ln(1/[1-F(x)]2),則SRI=k-(c0+c1k+c2k2/1+d1k+d2k2+d3k3),其中,c0=2.515 517;c1=0.802 853;c2=0.010 328;d1=1.432 788;d2=0.189 269;d3=0.001 308.

    計(jì)算SRI分析嘉陵江流域的旱澇演變特征,1月尺度的SRI反映流域短時(shí)間的旱澇狀況,12月尺度的SRI反映流域內(nèi)長時(shí)間旱澇的演變特征.此外,本文的旱澇等級確定參考國家標(biāo)準(zhǔn)《氣象干旱等級》(GB/T 20481-2006)[23,見表1.

    3.2 游程理論

    游程理論[24是一種分析時(shí)間序列的方法,是量化旱澇變量變化的理論規(guī)律的基礎(chǔ),根據(jù)游程理論三閾值方法25識(shí)別描述嘉陵江流域干旱和洪澇事件.

    在進(jìn)行旱澇分析時(shí),基于我國氣象干旱等級劃分標(biāo)準(zhǔn),并參考WU等[26有關(guān)研究,確定旱澇事件的劃分閾值.本研究取干旱事件的截取水平R0為-0.5(輕旱),當(dāng)SRI在某月開始小于-0.5,定義為一次干旱事件發(fā)生,到SRI值大于0.5,定義為此次干旱事件結(jié)束.從干旱開始到干旱結(jié)束的這段時(shí)間間隔定義為干旱事件的持續(xù)時(shí)間.干旱持續(xù)時(shí)間內(nèi)如果兩次干旱事件之間有1月的SRI值大于-0.5而小于0,則這兩次干旱可合并為一次干旱.

    干旱烈度S′定義為一次干旱事件中各指標(biāo)值與游程理論確定后的干旱閾值(R0)差的累計(jì)和,由于其為負(fù)數(shù),故一般取其絕對值.公式如下:

    S′=∑Tt=1|SRIt| (SRIt≤-0.5),

    式中,t為某次干旱(即SRI≤-0.5)開始的時(shí)刻,t=1;T為本次干旱結(jié)束的時(shí)刻.

    干旱強(qiáng)度指的是一次干旱事件中的平均缺水量,計(jì)算方式為某次干旱事件中干旱烈度與干旱歷時(shí)的比值.

    類似地,取洪澇事件的截取水平R0為0.5,洪澇過程識(shí)別與干旱過程識(shí)別同理.

    當(dāng)先后存在1個(gè)干旱和洪澇事件,且兩者之間的時(shí)間間隔小于或等于1月,定義為一次水文旱澇急轉(zhuǎn)事件,其中,一次旱澇急轉(zhuǎn)事件的起始時(shí)間是這次干旱事件的發(fā)生時(shí)間,結(jié)束時(shí)間是洪澇事件的結(jié)束時(shí)間,時(shí)段T內(nèi)發(fā)生旱澇急轉(zhuǎn)事件強(qiáng)度記為K,表2為旱澇急轉(zhuǎn)強(qiáng)度等級劃分標(biāo)準(zhǔn)[12.

    3.3 Mann-Kendall趨勢檢驗(yàn)法

    Mann-Kendall趨勢檢驗(yàn)法[27,目前得到廣泛應(yīng)用的趨勢性非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法.具體計(jì)算過程如下:假定x為時(shí)間序列變量,n為時(shí)間序列長度.

    S=∑ni=2i-1j=1sign(x0),

    式中,S為符號(hào)函數(shù),當(dāng)x0小于、等于或大于0時(shí),符號(hào)S等于-1,0和1.

    var(S)=(n(n-1)(2n-5)/18),

    式中,n為數(shù)據(jù)點(diǎn)個(gè)數(shù),統(tǒng)計(jì)檢驗(yàn)值Z計(jì)算式為:Z=((S-1)/var(S)),S>0;Z=0,S=0;Z=((S+1)/var(S)),S<0,若計(jì)算的絕對Z統(tǒng)計(jì)量大于正態(tài)分布表中Z統(tǒng)計(jì)量(1.65、1.96)的臨界值,則表示趨勢通過了90%和95%的顯著和強(qiáng)顯著水平.

    衡量趨勢大小的指標(biāo)為β=M((xi-xj/i-j)),式中,1<j<i<n(數(shù)據(jù)點(diǎn)個(gè)數(shù)),M代表求序列的中位數(shù).正的β值表示“上升趨勢”,負(fù)的β值表示“下降趨勢”.

    4 結(jié)果分析

    4.1 流域旱澇年際變化趨勢特征分析

    為更好地分析水文干旱趨勢和突變特征,選取1982-2020年嘉陵江流域北碚水文站的徑流深數(shù)據(jù),時(shí)間尺度分別設(shè)定為月尺度(1月)、年尺度(12月)兩種時(shí)間尺度,可以較好地識(shí)別流域內(nèi)短期和長期的旱澇狀況.如圖2(a),2(b)別為月尺度和年尺度的SRI序列變化趨勢圖.

    不同時(shí)間尺度下的SRI波動(dòng)幅度不同,但干旱趨勢大致相同.月尺度的SRI波動(dòng)較強(qiáng),表現(xiàn)出對水文干旱的強(qiáng)烈反應(yīng).在月尺度的SRI序列中,達(dá)到特旱水平的是1987年3月、1997年8至9月、1998年1月、2002年10月、2003年1月、2006年7月.分析年尺度SRI可得:嘉陵江流域在1982-2020年呈現(xiàn)旱澇交替現(xiàn)象,但干旱年份居多;1982-2009年洪澇化趨勢先減弱后增強(qiáng),2010年后洪澇化趨勢呈短期增強(qiáng)趨勢后又減弱;此外,1992-2008年、2014-2017年嘉陵江流域干旱化趨勢明顯,旱災(zāi)較為嚴(yán)重;2009-2013年、2018-2020年SRI為正值,澇災(zāi)嚴(yán)重.

    通過Mann-Kendall趨勢檢驗(yàn)對1982-2020年嘉陵江流域北碚水文站月尺度和年尺度的SRI序列進(jìn)行趨勢檢驗(yàn).結(jié)果如表3所示.

    從年尺度的SRI序列來看,統(tǒng)計(jì)量β值小于0,且Z值絕對值0.169 4小于1.65,表明干旱化呈不顯著增強(qiáng)的趨勢;2、3和4月份北碚水文站Mann-Kendall統(tǒng)計(jì)值通過了顯著性檢驗(yàn),發(fā)生的顯著變化趨勢是上升趨勢,表明干旱化趨勢均減弱;7、9和10月份呈現(xiàn)顯著下降趨勢,表明洪澇化趨勢減弱,干旱化趨勢增強(qiáng);1、5、6、8、11和12月份呈現(xiàn)不明顯的干旱化減弱趨勢.

    4.2 流域不同年代旱澇年際變化特征分析

    根據(jù)上文計(jì)算得出的SRI,進(jìn)一步識(shí)別嘉陵江北碚水文站的旱澇事件,將1982-2020年按年代劃分并分別統(tǒng)計(jì)北碚水文站發(fā)生旱澇事件的次數(shù)、歷時(shí)和烈度以及發(fā)生旱澇事件最大歷時(shí)和最大烈度,見表4.

    由表4 可知,20世紀(jì)80年代澇災(zāi)發(fā)生程度和歷時(shí)總體比旱災(zāi)高,整體偏澇;20世紀(jì)90年代,北碚水文站發(fā)生干旱事件的次數(shù)增加,平均歷時(shí)、平均烈度、最長歷時(shí)、最大烈度增大;21世紀(jì)00年代干旱事件發(fā)生頻次進(jìn)一步增多,但平均歷時(shí)、平均烈度、最長歷時(shí)和最大烈度均減少;21世紀(jì)10年代干旱平均歷時(shí)、平均烈度、最長歷時(shí)和最大烈度減小.

    20世紀(jì)90年代和20世紀(jì)80年代相比,北碚水文站洪澇事件發(fā)生次數(shù)不變,洪澇程度減輕,總體上洪澇化趨勢減弱;21世紀(jì)00年代到21世紀(jì)10年代洪澇事件的發(fā)生次數(shù)和平均歷時(shí)、平均烈度、最長洪澇歷時(shí)和最大烈度增加,表明洪澇事件頻繁,且洪澇事件程度加重.20世紀(jì)80年代和21世紀(jì)10年代洪澇事件發(fā)生的頻率和程度高于干旱事件,而20世紀(jì)90年代到21世紀(jì)00年代干旱事件的平均歷時(shí)、平均烈度和最長歷時(shí)都高于洪澇事件.

    4.3 流域旱澇急轉(zhuǎn)事件演變特征分析

    表5統(tǒng)計(jì)了嘉陵江流域北碚水文站旱澇急轉(zhuǎn)事件特征,包括急轉(zhuǎn)事件的起始和終止時(shí)間,強(qiáng)度和急轉(zhuǎn)點(diǎn)月份.北碚水文站旱澇急轉(zhuǎn)事件共發(fā)生39次.旱澇急轉(zhuǎn)次數(shù)最多的是21世紀(jì)10年代,旱澇急轉(zhuǎn)事件發(fā)生14次,20世紀(jì)90年代旱澇急轉(zhuǎn)事件發(fā)生次數(shù)最少,為6次.21世紀(jì)以來,旱澇急轉(zhuǎn)事件發(fā)生次數(shù)明顯增多.

    圖3給出了嘉陵江流域各月旱澇急轉(zhuǎn)頻次和強(qiáng)度的變化,由圖3可知,嘉陵江流域旱澇急轉(zhuǎn)多發(fā)生在5月份,這是因?yàn)榧瘟杲饔虼杭据^為干旱,夏季突發(fā)降水容易造成旱澇急轉(zhuǎn).旱澇急轉(zhuǎn)強(qiáng)度最大值出現(xiàn)在4月份,該時(shí)段對應(yīng)的旱澇急轉(zhuǎn)頻次也較多,因此發(fā)生洪澇災(zāi)害的風(fēng)險(xiǎn)比較大.嘉陵江流域的旱澇急轉(zhuǎn)事件主要發(fā)生在5、6、8和9月份,且時(shí)段內(nèi)旱澇急轉(zhuǎn)事件的強(qiáng)度也較高.同時(shí),旱轉(zhuǎn)澇近40年發(fā)生頻率依次升高,并且旱澇急轉(zhuǎn)事件的強(qiáng)度也呈升高的趨勢,應(yīng)予以足夠重視.

    5 結(jié)論與討論

    本文采用嘉陵江流域北碚水文站1982-2020年的徑流深數(shù)據(jù),計(jì)算不同尺度的標(biāo)準(zhǔn)化徑流指數(shù),采用游程理論定義旱澇事件,并根據(jù)旱澇急轉(zhuǎn)條件篩選出旱澇急轉(zhuǎn)事件,計(jì)算發(fā)生頻次和烈度,得出主要結(jié)論如下:

    (1)嘉陵江流域1982-2009年洪澇化趨勢先減弱后增強(qiáng),2010-2020年洪澇化趨勢呈短期的增強(qiáng)趨勢后又減弱.近20年來流域內(nèi)洪澇事件發(fā)生愈加頻繁,且場次洪澇程度依次加重,并且呈現(xiàn)出洪澇極端化趨勢.

    (2)嘉陵江流域1982-2020年的1、5、6、8、11和12月份呈現(xiàn)不明顯的干旱化減弱趨勢;7、9和10月份干旱化趨勢增強(qiáng),2、3和4月份干旱化趨勢顯著減弱.

    (3)嘉陵江北碚水文站旱澇急轉(zhuǎn)事件共發(fā)生39次,旱轉(zhuǎn)澇事件急轉(zhuǎn)點(diǎn)均集中在4月份,對應(yīng)旱澇事件的強(qiáng)度也較高.旱澇急轉(zhuǎn)事件近40年發(fā)生頻率依次升高,并伴有逐步加重的趨勢,應(yīng)予以足夠重視.

    雖然本研究采用SRI指數(shù)和游程理論識(shí)別了嘉陵江流域1982-2020年旱澇事件特征以及旱澇急轉(zhuǎn)事件演變規(guī)律,但仍有所不足,有待繼續(xù)研究和完善.

    (1)本文計(jì)算了嘉陵江流域北碚水文站1982-2020年不同尺度的標(biāo)準(zhǔn)化徑流指數(shù)值,分析嘉陵江流域旱澇事件演變特征以及變化規(guī)律.但僅選用標(biāo)準(zhǔn)化徑流指數(shù)指標(biāo)對研究區(qū)旱澇狀況進(jìn)行討論,不能對流域存在多種旱澇類型的情況進(jìn)行分析,因此,未來的研究工作可以利用更多指標(biāo)參數(shù)如標(biāo)準(zhǔn)化降水指數(shù)(SPI)、標(biāo)準(zhǔn)化加權(quán)平均降雨(SWAP)指數(shù)等分別對氣象旱澇和水文旱澇進(jìn)行研究.

    (2)本研究僅采用位于嘉陵江下游北碚水文站1982-2020年的徑流數(shù)據(jù)計(jì)算多時(shí)間尺度SRI值,分析其水文旱澇狀況.但隨著氣候變化加劇,洪澇、干旱等自然災(zāi)害頻發(fā),嘉陵江流域發(fā)生旱澇急轉(zhuǎn)事件的無常性大大增加,同一條河流的不同站點(diǎn)旱澇急轉(zhuǎn)事件特征會(huì)有很大的差異,因此未來的研究工作可以著重分析嘉陵江多個(gè)水文站點(diǎn)如亭子口和武勝站的旱澇急轉(zhuǎn)事件的空間演變規(guī)律.

    (3)氣象旱澇和水文旱澇往往有密切的聯(lián)系:降水的虧缺和盈余會(huì)導(dǎo)致氣象旱澇,而地表水與地下水的虧缺和盈余會(huì)引起水文旱澇,氣象旱澇往往會(huì)造成水文旱澇.后續(xù)研究中,可以嘗試準(zhǔn)確描述嘉陵江流域氣象旱澇和水文旱澇的關(guān)系,以便在氣象旱澇事件發(fā)生初期,及時(shí)采取應(yīng)對措施,減輕水文旱澇造成的災(zāi)害等級.

    參 考 文 獻(xiàn)

    [1] 巢清塵.氣候變化背景下我國水旱災(zāi)害演變趨勢及其應(yīng)對思考[J].中國減災(zāi),2024(15):8-11.

    CHAO Q C.Evolution trend of flood and drought disasters in China under the background of climate change and its countermeasures[J].Disaster Reduction in China,2024(15):8-11.

    [2]李娜.干旱氣象災(zāi)害研究現(xiàn)狀和進(jìn)展[J].農(nóng)業(yè)科學(xué),2021,11(5):413-416.

    LI N.Research Status and Progress of Drought Meteorological Disasters[J].Hans Journal of Agricultural Sciences,2021,11(5):413-416.

    [3]宋艷玲.全球干旱指數(shù)研究進(jìn)展[J].應(yīng)用氣象學(xué)報(bào),2022,33(5):513-526.

    SONG Y L.Global research progress of drought indices[J].Journal of Applied Meteorological Science,2022,33(5):513-526.

    [4]宋琳琳,張強(qiáng),任余龍,等.PDSI及sc_PDSI干旱指數(shù)在中國西南地區(qū)適用性分析[J].中國沙漠,2021,41(2):242-251.

    SONG L L,ZHANG Q,REN Y L,et al.The applicable analysis of PDSI and self_calibrating PDSI drought indices in Southwest China[J].Journal of Desert Research,2021,41(2):242-251.

    [5]ZSCHEISCHLER J,MARTIUS O,WESTRA S,et al.A typology of compound weather and climate events[J].Nature Reviews Earth amp; Environment,2020,1:333-347.

    [6]楊家偉,陳華,侯雨坤,等.基于氣象旱澇指數(shù)的旱澇急轉(zhuǎn)事件識(shí)別方法[J].地理學(xué)報(bào),2019,74(11):2358-2370.

    YANG J W,CHEN H,HOU Y K,et al.A method to identify the drought-flood transition based on the meteorological drought index[J].Acta Geographica Sinica,2019,74(11):2358-2370.

    [7]趙東升,張家誠,鄧思琪,等.1960-2018年中國西南地區(qū)旱澇急轉(zhuǎn)的時(shí)空變化特征[J].地理科學(xué),2021,41(12):2222-2231.

    ZHAO D S,ZHANG J C,DENG S Q,et al.Spatio-temporal characteristics of drought-flood abrupt alternation in the Southwest China from 1960 to 2018[J].Scientia Geographica Sinica,2021,41(12):2222-2231.

    [8]吳志偉,李建平,何金海,等.大尺度大氣環(huán)流異常與長江中下游夏季長周期旱澇急轉(zhuǎn)[J].科學(xué)通報(bào),2006,51(14):1717-1724.

    WU Z W,LI J P,HE J H,et al.Large-scale atmospheric circulation anomaly and long-term drought and flood sharp turn in summer in the middle and lower reaches of the Yangtze River[J].Chinese Science Bulletin,2006,51(14):1717-1724.

    [9]張?jiān)品?,翟麗妮,林沛榕,?長江中下游典型流域旱澇與旱澇/澇旱急轉(zhuǎn)演變規(guī)律及其驅(qū)動(dòng)因子研究[J].武漢大學(xué)學(xué)報(bào)(工學(xué)版),2021,54(10):887-897.

    ZHANG Y F,ZHAI L N,LIN P R,et al.Variation characteristics and driving factors of drought and flood and their abrupt alternations in a typical basin in the middle reaches of Yangtze River[J].Engineering Journal of Wuhan University,2021,54(10):887-897.

    [10]張玉琴,李棟梁.華南汛期旱澇急轉(zhuǎn)及其大氣環(huán)流特征[J].氣候與環(huán)境研究,2019,24(4):430-444.

    ZHANG Y Q,LI D L.Drought-flood abrupt alternation and its atmospheric circulation characteristics during flood season in Southern China[J].Climatic and Environmental Research,2019,24(4):430-444.

    [11]閃麗潔,張利平,陳心池,等.長江中下游流域旱澇急轉(zhuǎn)時(shí)空演變特征分析[J].長江流域資源與環(huán)境,2015,24(12):2100-2107.

    SHAN L J,ZHANG L P,CHEN X C,et al.Spatio-temporal evolution characteristics of droughtflood abrupt alternation in the middle and lower reaches of the Yangtze River basin[J].Resources and Environment in the Yangtze Basin,2015,24(12):2100-2107.

    [12]趙英,陳華,楊家偉,等.基于SWAP和SRI的漢江流域旱澇急轉(zhuǎn)時(shí)空特征分析[J].人民長江,2020,51(4):94-99.

    ZHAO Y,CHEN H,YANG J W,et al.Spatio-temporal variations of drought-flood abrupt alternation events in Hanjiang River Basin based on SWAP and SRI[J].Yangtze River,2020,51(4):94-99.

    [13]胡毅鴻,李景保.1951-2015年洞庭湖區(qū)旱澇演變及典型年份旱澇急轉(zhuǎn)特征分析[J].農(nóng)業(yè)工程學(xué)報(bào),2017,33(7):107-115.

    HU Y H,LI J B.Analysis on evolution of drought-flood and its abrupt alternation in typical year from 1951 to 2015 in Dongting Lake area[J].Transactions of the Chinese Society of Agricultural Engineering,2017,33(7):107-115.

    [14]王容,李相虎,薛晨陽,等.1960—2012年鄱陽湖流域旱澇急轉(zhuǎn)事件時(shí)空演變特征[J].湖泊科學(xué),2020,32(1):207-222.

    WANG R,LI X H,XUE C Y,et al.Spatio-temporal variations of the drought-flood abrupt alternation events in the Lake Poyang Basin from 1960 to 2012[J].Journal of Lake Sciences,2020,32(1):207-222.

    [15]葉磊,周建中,曾小凡,等.氣候變化下SPEI指數(shù)在嘉陵江流域的干旱評估應(yīng)用[J].長江流域資源與環(huán)境,2015,24(6):943-948.

    YE L,ZHOU J Z,ZENG X F,et al.Application of spei for the changes of drought in Jialing river basin under climate change[J].Resources and Environment in the Yangtze Basin,2015,24(6):943-948.

    [16]鄧慧平,李愛貞,劉厚風(fēng),等.氣候波動(dòng)對萊州灣地區(qū)水資源及極端旱澇事件的影響[J].地理科學(xué),2000,20(1):56-60.

    DENG H P,LI A Z,LIU H F,et al.Impacts of climate fluctuations on water resources and extremes in the Laizhou Bay area[J].Scientia Geographica Sinica,2000,20(1):56-60.

    [17]孫浚凱.基于SRI的渾河流域水文干旱特征分析[J].人民珠江,2021,42(2):25-30.

    SUN J K.Analysis on drought characteristics in the Hunhe River Basin based on SRI[J].Pearl River,2021,42(2):25-30.

    [18]趙子萌,邵薇薇,曹永強(qiáng),等.遼寧省旱澇急轉(zhuǎn)特征與趨勢研究[J].水利水電技術(shù)(中英文),2024,55(6):32-43.

    ZHAO Z M,SHAO W W,CAO Y Q,et al.Study on the characteristics and trend of drought-flood abrupt alternation in Liaoning Province[J].Water Resources and Hydropower Engineering,2024,55(6):32-43.

    [19]劉玉冰,劉波,王文鵬,等.基于日尺度旱澇急轉(zhuǎn)指數(shù)的重慶市旱澇急轉(zhuǎn)事件演變特征分析[J].水資源與水工程學(xué)報(bào),2024,35(4):20-28.

    LIU Y B,LIU B,WANG W P,et al.Evolution of drought-flood abrupt alternation events in Chongqing based on the daily dry-wet abrupt alternation index[J].Journal of Water Resources and Water Engineering,2024,35(4):20-28.

    [20]ZHANG B W,CHEN Y,CHEN X W,et al.Spatial-temporal variations of drought-flood abrupt alternation events in Southeast China[J].Water,2024,16(3):498.

    [21]MCKEE T B,J DOESKEN N,J KLEIST.The relationship of drought frequency and duration to time scales[J].Proceedings of the 8th Conference on Applied Climatology,1993:179-183.

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

    [23]中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局.機(jī)械制圖、圖樣畫法、圖線[M].北京:中國標(biāo)準(zhǔn)出版社,2004.

    [24]ENDT P M,M PATTERN D,M BUECHNER W,et al.An objective approach to definitions and investigations of continental hydrologic droughts vujica yevjevich:fort collins,Colorado state university,1967,19 p.(hydrology paper no.23)[J].Journal of Hydrology,1969,7(3):353.

    [25]CHANG J X,LI Y Y,WANG Y M,et al.Copula-based drought risk assessment combined with an integrated index in the Wei River Basin,China[J].Journal of Hydrology,2016,540:824-834.

    [26]WU J F,CHEN X W,YAO H X,et al.Non-linear relationship of hydrological drought responding to meteorological drought and impact of a large reservoir[J].Journal of Hydrology,2017,551:495-507.

    [27]姜瑤,徐宗學(xué),王靜.基于年徑流序列的五種趨勢檢測方法性能對比[J].水利學(xué)報(bào),2020,51(7):845-857.

    JIANG Y,XU Z X,WANG J.Comparison among five methods of trend detection for annual runoff series[J].Journal of Hydraulic Engineering,2020,51(7):845-857.

    Analysis on the characteristics of drought-flood abrupt alternation

    in the Jialing River Basin from 1982 to 2020 based on SRIXu Jingwen1,2, Ji Guangxing1, Zhang Yali1

    (1. College of Resources and Environment, Henan Agricultural University, Zhengzhou 450002, China;

    2. Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China)

    Abstract: Basing on the monthly runoff data of Beibei hydrological station in the Jialing River Basin from 1982 to 2020, we calculate the different time scales of standardized runoff index(SRI). Using the run-time theory to analyze the drought-flood abrupt alternation characteristic and utilizing Mann-Kendall test to analyze the trend of SRI change in the Jialing River Basin. The results show that:During 1982 to 2009, the trend of flooding in the basin weakened first and then increased, from 2010 to 2020, the trend of flooding increased for a short period and then decreased again. From 1982 to 2020, there was an insignificant weakening trend of drought in January, May, June, August, November and December. The trend of drying increased in July, September and October, and decreased significantly in February, March and April. A total of 39 times of drought-flood abrupt alternation events were identified at the Beibei hydrological station of the Jialing River Basin. The abrupt alternation points of drought-flood abrupt alternation events were concentrated in April and the intensity of drought-flood events was also high. The frequency of drought-flood abrupt alternation events in" near 40 years has gradually increased, and shows a trend of gradual aggravation, which should be given sufficient attention.

    Keywords: standardized runoff index; Jialing River Basin; drought-flood abrupt alternation; Mann-Kendall test

    [責(zé)任編校 趙曉華 劉洋]

    男女之事视频高清在线观看| 亚洲不卡免费看| 一级毛片aaaaaa免费看小| 伊人久久精品亚洲午夜| 国产一区二区三区av在线 | av女优亚洲男人天堂| 一本精品99久久精品77| 九九在线视频观看精品| 高清日韩中文字幕在线| 一级毛片我不卡| 午夜精品国产一区二区电影 | 你懂的网址亚洲精品在线观看 | av在线观看视频网站免费| 日韩三级伦理在线观看| 好男人在线观看高清免费视频| 久久精品国产鲁丝片午夜精品| 黄色配什么色好看| 最新中文字幕久久久久| 少妇猛男粗大的猛烈进出视频 | АⅤ资源中文在线天堂| 中国美白少妇内射xxxbb| 免费无遮挡裸体视频| 欧美又色又爽又黄视频| 婷婷六月久久综合丁香| 成人毛片a级毛片在线播放| 午夜a级毛片| 中国国产av一级| 少妇猛男粗大的猛烈进出视频 | 97人妻精品一区二区三区麻豆| 亚洲av.av天堂| 国产一区二区亚洲精品在线观看| 午夜亚洲福利在线播放| 毛片女人毛片| 男人狂女人下面高潮的视频| 干丝袜人妻中文字幕| 久久久久久久久久久丰满| 亚洲欧美成人精品一区二区| 婷婷精品国产亚洲av| 久久久午夜欧美精品| 精品一区二区免费观看| 99在线人妻在线中文字幕| 97在线视频观看| 欧美日韩乱码在线| 午夜老司机福利剧场| 少妇熟女欧美另类| 欧美国产日韩亚洲一区| 国产探花在线观看一区二区| 日日摸夜夜添夜夜添av毛片| 亚洲成人av在线免费| 在线观看av片永久免费下载| av国产免费在线观看| 国产精品一区二区三区四区久久| 嫩草影院精品99| 成人无遮挡网站| 麻豆精品久久久久久蜜桃| 99久久成人亚洲精品观看| 在线播放无遮挡| 日韩av不卡免费在线播放| 欧美性猛交黑人性爽| 51国产日韩欧美| 国产激情偷乱视频一区二区| 国产午夜精品久久久久久一区二区三区 | 国产一区二区亚洲精品在线观看| 欧美丝袜亚洲另类| 成人永久免费在线观看视频| 色在线成人网| 国产精品女同一区二区软件| 亚洲人成网站在线播| 欧美激情久久久久久爽电影| 香蕉av资源在线| 国产成人影院久久av| 观看免费一级毛片| 午夜福利高清视频| 免费观看人在逋| a级毛色黄片| 亚洲人成网站在线观看播放| 91av网一区二区| 精品免费久久久久久久清纯| 精品久久久久久久久久久久久| 99视频精品全部免费 在线| 亚洲一级一片aⅴ在线观看| 给我免费播放毛片高清在线观看| 成人精品一区二区免费| 日韩精品青青久久久久久| 亚洲真实伦在线观看| 亚洲人成网站高清观看| 99热这里只有精品一区| 真实男女啪啪啪动态图| 久久久久久久久久成人| 亚洲成人久久性| 久久这里只有精品中国| 免费大片18禁| 国产精品一区www在线观看| 亚洲美女搞黄在线观看 | 久久中文看片网| 久久6这里有精品| 99热6这里只有精品| 国产乱人偷精品视频| 黄片wwwwww| av国产免费在线观看| 国产精品免费一区二区三区在线| 九九热线精品视视频播放| 一夜夜www| 最后的刺客免费高清国语| 国产亚洲欧美98| 国产aⅴ精品一区二区三区波| 亚洲成av人片在线播放无| 女同久久另类99精品国产91| 国产三级中文精品| 亚洲国产欧洲综合997久久,| av卡一久久| 乱系列少妇在线播放| 日韩欧美国产在线观看| 搡老熟女国产l中国老女人| 亚洲人成网站在线播放欧美日韩| 又黄又爽又免费观看的视频| av天堂中文字幕网| av天堂在线播放| 日韩三级伦理在线观看| 欧美在线一区亚洲| 热99re8久久精品国产| 老师上课跳d突然被开到最大视频| 内射极品少妇av片p| 国内少妇人妻偷人精品xxx网站| 久久婷婷人人爽人人干人人爱| 九九热线精品视视频播放| 久久国内精品自在自线图片| 亚洲欧美精品综合久久99| 真人做人爱边吃奶动态| 久久精品91蜜桃| 91麻豆精品激情在线观看国产| 两个人视频免费观看高清| 成人亚洲精品av一区二区| 熟女电影av网| 久久综合国产亚洲精品| 亚洲第一电影网av| 此物有八面人人有两片| 一区二区三区免费毛片| 人人妻人人看人人澡| 无遮挡黄片免费观看| 国产成人a∨麻豆精品| 国产高清三级在线| 免费看a级黄色片| 亚洲aⅴ乱码一区二区在线播放| 老师上课跳d突然被开到最大视频| 悠悠久久av| 18禁在线播放成人免费| 99热这里只有是精品50| 亚洲自拍偷在线| 久久午夜福利片| 国产探花在线观看一区二区| 深夜a级毛片| 久久午夜福利片| 99精品在免费线老司机午夜| 少妇丰满av| 男插女下体视频免费在线播放| 天堂av国产一区二区熟女人妻| 国产一区二区三区av在线 | 欧美中文日本在线观看视频| 韩国av在线不卡| 国产精品久久视频播放| 国产v大片淫在线免费观看| 久久久久久久久久成人| 22中文网久久字幕| 国产精品一区二区三区四区久久| 久久国内精品自在自线图片| 青春草视频在线免费观看| 你懂的网址亚洲精品在线观看 | 中文字幕精品亚洲无线码一区| 色尼玛亚洲综合影院| 伊人久久精品亚洲午夜| 国产男靠女视频免费网站| 久久久精品94久久精品| 亚洲av免费高清在线观看| 中文字幕人妻熟人妻熟丝袜美| 大香蕉久久网| 91av网一区二区| 18禁在线播放成人免费| 一级黄片播放器| 国产欧美日韩精品一区二区| 亚洲中文字幕一区二区三区有码在线看| 国产精品一及| 亚洲精品影视一区二区三区av| 成人亚洲精品av一区二区| 成人高潮视频无遮挡免费网站| 亚洲成人中文字幕在线播放| 亚洲四区av| av在线亚洲专区| 99热网站在线观看| 亚洲一区高清亚洲精品| 老熟妇仑乱视频hdxx| 日本a在线网址| 国产在线精品亚洲第一网站| 国产一区二区在线观看日韩| 免费在线观看成人毛片| 国产高清不卡午夜福利| 嫩草影院精品99| av黄色大香蕉| 3wmmmm亚洲av在线观看| 热99re8久久精品国产| 三级经典国产精品| 亚洲av一区综合| 亚洲天堂国产精品一区在线| 日产精品乱码卡一卡2卡三| 在线观看午夜福利视频| 成人永久免费在线观看视频| 久久人人爽人人爽人人片va| 99久国产av精品| 精品一区二区三区人妻视频| 一夜夜www| 中文字幕久久专区| 国产一区二区在线观看日韩| 久久久久性生活片| 亚洲成人久久性| 日日摸夜夜添夜夜爱| 国产欧美日韩精品一区二区| 日韩欧美三级三区| 91精品国产九色| 亚洲成人久久性| 亚洲国产精品成人综合色| 99精品在免费线老司机午夜| 高清日韩中文字幕在线| 午夜精品国产一区二区电影 | 日韩精品有码人妻一区| 午夜激情福利司机影院| 在线观看66精品国产| 色哟哟·www| 国产久久久一区二区三区| 国产精品亚洲一级av第二区| 国产精品女同一区二区软件| 午夜精品一区二区三区免费看| 国产 一区精品| 国产欧美日韩一区二区精品| 国产精品国产高清国产av| 97在线视频观看| 久久精品国产鲁丝片午夜精品| 国产三级中文精品| 国产精华一区二区三区| 国产黄片美女视频| 国产高清有码在线观看视频| 午夜免费男女啪啪视频观看 | 成人精品一区二区免费| 欧美精品国产亚洲| 中文字幕人妻熟人妻熟丝袜美| 亚洲美女黄片视频| 三级毛片av免费| 国产精品亚洲一级av第二区| 卡戴珊不雅视频在线播放| av在线老鸭窝| 22中文网久久字幕| 国产午夜福利久久久久久| 午夜福利在线观看免费完整高清在 | 免费观看精品视频网站| 不卡一级毛片| 天堂动漫精品| 美女cb高潮喷水在线观看| 久久久a久久爽久久v久久| av.在线天堂| 久久久久久国产a免费观看| 可以在线观看的亚洲视频| 人妻丰满熟妇av一区二区三区| av.在线天堂| 亚洲av免费在线观看| 国产精华一区二区三区| 久久精品国产鲁丝片午夜精品| 色综合色国产| 菩萨蛮人人尽说江南好唐韦庄 | 国产在线精品亚洲第一网站| 12—13女人毛片做爰片一| 国国产精品蜜臀av免费| 桃色一区二区三区在线观看| 不卡视频在线观看欧美| 久久这里只有精品中国| 给我免费播放毛片高清在线观看| a级一级毛片免费在线观看| 久久人人精品亚洲av| 欧美在线一区亚洲| 老女人水多毛片| 欧美一区二区精品小视频在线| 91久久精品国产一区二区三区| 国产v大片淫在线免费观看| 欧美区成人在线视频| 亚洲精品国产成人久久av| 女人被狂操c到高潮| 三级经典国产精品| 女人十人毛片免费观看3o分钟| 国产亚洲精品av在线| 欧美精品国产亚洲| 国产aⅴ精品一区二区三区波| av国产免费在线观看| 少妇熟女欧美另类| 亚洲av中文av极速乱| 国产精品嫩草影院av在线观看| 国产精品一区二区三区四区免费观看 | 51国产日韩欧美| 免费av不卡在线播放| 亚洲性夜色夜夜综合| 国产在线精品亚洲第一网站| 我的女老师完整版在线观看| 蜜桃久久精品国产亚洲av| 久久欧美精品欧美久久欧美| 最近手机中文字幕大全| 在线观看66精品国产| 综合色丁香网| 看黄色毛片网站| 免费不卡的大黄色大毛片视频在线观看 | 欧美最黄视频在线播放免费| 丰满人妻一区二区三区视频av| 亚洲精品456在线播放app| 日韩欧美 国产精品| 亚洲精品影视一区二区三区av| 日韩在线高清观看一区二区三区| 亚洲熟妇中文字幕五十中出| 国产 一区 欧美 日韩| www日本黄色视频网| 韩国av在线不卡| 人妻夜夜爽99麻豆av| 卡戴珊不雅视频在线播放| 久久亚洲国产成人精品v| 精品午夜福利在线看| 色尼玛亚洲综合影院| av国产免费在线观看| 久久久久性生活片| 此物有八面人人有两片| 热99re8久久精品国产| 亚洲一区二区三区色噜噜| 亚洲高清免费不卡视频| 又黄又爽又免费观看的视频| av视频在线观看入口| 国产亚洲精品综合一区在线观看| 精品久久久久久久久av| av在线老鸭窝| 国产精品国产三级国产av玫瑰| 亚洲人成网站在线播| 无遮挡黄片免费观看| 少妇人妻精品综合一区二区 | 欧美最新免费一区二区三区| 精品免费久久久久久久清纯| 内射极品少妇av片p| or卡值多少钱| 国产精品人妻久久久久久| 亚洲av成人av| 午夜精品一区二区三区免费看| 久久久久久伊人网av| 精品午夜福利视频在线观看一区| 亚洲乱码一区二区免费版| 国产麻豆成人av免费视频| 国产一级毛片七仙女欲春2| 俺也久久电影网| 一个人看的www免费观看视频| 午夜亚洲福利在线播放| 成人永久免费在线观看视频| 国产精品女同一区二区软件| 熟妇人妻久久中文字幕3abv| 欧美日韩乱码在线| 搞女人的毛片| 嫩草影院新地址| 搡老熟女国产l中国老女人| 一级毛片电影观看 | 乱人视频在线观看| 中文资源天堂在线| 亚洲性夜色夜夜综合| 亚洲天堂国产精品一区在线| 在线观看av片永久免费下载| 亚洲aⅴ乱码一区二区在线播放| 久久韩国三级中文字幕| 日本撒尿小便嘘嘘汇集6| 精品不卡国产一区二区三区| 日本免费a在线| 狂野欧美白嫩少妇大欣赏| 高清午夜精品一区二区三区 | avwww免费| 99热6这里只有精品| 不卡一级毛片| 亚洲高清免费不卡视频| 亚洲精品日韩av片在线观看| 99久久无色码亚洲精品果冻| 国产av麻豆久久久久久久| 丰满的人妻完整版| 老司机福利观看| 深夜a级毛片| 秋霞在线观看毛片| 国产亚洲精品久久久久久毛片| 九九热线精品视视频播放| 亚洲久久久久久中文字幕| 久久久午夜欧美精品| 久久午夜福利片| 日韩在线高清观看一区二区三区| 一区福利在线观看| 欧美色视频一区免费| 亚洲欧美日韩东京热| 波多野结衣高清作品| 蜜桃亚洲精品一区二区三区| 国产免费一级a男人的天堂| 丝袜美腿在线中文| a级一级毛片免费在线观看| 婷婷六月久久综合丁香| 人妻少妇偷人精品九色| 人人妻人人澡人人爽人人夜夜 | 性欧美人与动物交配| av在线观看视频网站免费| 干丝袜人妻中文字幕| 91av网一区二区| 校园春色视频在线观看| 色吧在线观看| 亚洲在线自拍视频| 国产精品日韩av在线免费观看| 看黄色毛片网站| 偷拍熟女少妇极品色| 伊人久久精品亚洲午夜| 中国美女看黄片| 亚洲电影在线观看av| 精品一区二区三区视频在线| av天堂在线播放| 中出人妻视频一区二区| 久久99热这里只有精品18| 丝袜喷水一区| 日本精品一区二区三区蜜桃| 看黄色毛片网站| 在线观看美女被高潮喷水网站| 亚洲最大成人中文| 一级毛片我不卡| 色哟哟哟哟哟哟| 亚洲国产精品成人久久小说 | 嫩草影院精品99| 搞女人的毛片| 午夜老司机福利剧场| 欧美三级亚洲精品| 亚洲人与动物交配视频| 日本色播在线视频| 国内久久婷婷六月综合欲色啪| 在线免费观看的www视频| 亚洲丝袜综合中文字幕| 免费看美女性在线毛片视频| 99久久九九国产精品国产免费| 亚洲精品久久国产高清桃花| 成人毛片a级毛片在线播放| 在线看三级毛片| 亚洲成人中文字幕在线播放| 老师上课跳d突然被开到最大视频| 欧美精品国产亚洲| 久久久久免费精品人妻一区二区| 少妇猛男粗大的猛烈进出视频 | 三级国产精品欧美在线观看| 国产黄色视频一区二区在线观看 | 免费不卡的大黄色大毛片视频在线观看 | 亚洲成人精品中文字幕电影| 欧美色欧美亚洲另类二区| 亚洲精品在线观看二区| 人妻久久中文字幕网| 久久99热这里只有精品18| 亚洲人成网站高清观看| 激情 狠狠 欧美| 人人妻人人看人人澡| 99热全是精品| 亚洲精品日韩av片在线观看| 中文在线观看免费www的网站| 久久中文看片网| 大又大粗又爽又黄少妇毛片口| 综合色av麻豆| 久久久久国产精品人妻aⅴ院| 最后的刺客免费高清国语| 国产高潮美女av| a级毛色黄片| 日韩中字成人| 国产aⅴ精品一区二区三区波| 长腿黑丝高跟| 99热精品在线国产| 五月玫瑰六月丁香| 国产综合懂色| 国产精品乱码一区二三区的特点| 国产av在哪里看| 久久国内精品自在自线图片| 国产白丝娇喘喷水9色精品| 免费看美女性在线毛片视频| 国产在线精品亚洲第一网站| 午夜影院日韩av| 久久国内精品自在自线图片| 精品福利观看| 69人妻影院| 免费一级毛片在线播放高清视频| 91麻豆精品激情在线观看国产| 直男gayav资源| 久久亚洲精品不卡| 极品教师在线视频| 美女免费视频网站| 成人鲁丝片一二三区免费| 亚洲精品一区av在线观看| 一进一出抽搐动态| 悠悠久久av| 国内精品一区二区在线观看| 日韩三级伦理在线观看| 观看免费一级毛片| 五月伊人婷婷丁香| 亚洲成a人片在线一区二区| 日日啪夜夜撸| 熟女人妻精品中文字幕| 九九在线视频观看精品| 性色avwww在线观看| 国产成人精品久久久久久| 免费在线观看成人毛片| 久久久成人免费电影| 国产av麻豆久久久久久久| 日日干狠狠操夜夜爽| 久久精品91蜜桃| 国产视频内射| 国产成人91sexporn| 少妇的逼水好多| 免费搜索国产男女视频| 国产高清视频在线播放一区| 乱系列少妇在线播放| 大型黄色视频在线免费观看| 亚洲欧美中文字幕日韩二区| 人人妻人人看人人澡| 国产精品不卡视频一区二区| 久久精品综合一区二区三区| 国产精品美女特级片免费视频播放器| 国产精品久久久久久亚洲av鲁大| 美女高潮的动态| 婷婷色综合大香蕉| 欧美成人精品欧美一级黄| 成人国产麻豆网| 亚洲国产精品国产精品| 99热只有精品国产| 国产色婷婷99| 色哟哟·www| 亚洲成av人片在线播放无| 免费看美女性在线毛片视频| 伊人久久精品亚洲午夜| 六月丁香七月| 国产熟女欧美一区二区| 禁无遮挡网站| 一区二区三区高清视频在线| 可以在线观看的亚洲视频| 两个人的视频大全免费| 国产淫片久久久久久久久| 成人毛片a级毛片在线播放| 一进一出抽搐动态| av天堂在线播放| 一级av片app| 成人国产麻豆网| av黄色大香蕉| 国产一区二区在线av高清观看| 亚洲成av人片在线播放无| 我要搜黄色片| 欧美激情国产日韩精品一区| 亚洲欧美日韩卡通动漫| 婷婷精品国产亚洲av| 高清毛片免费观看视频网站| 亚洲中文字幕一区二区三区有码在线看| av女优亚洲男人天堂| 亚洲av中文av极速乱| 一级av片app| 日本 av在线| 一本精品99久久精品77| 国产av不卡久久| 欧美高清成人免费视频www| 国产精品爽爽va在线观看网站| 国产免费男女视频| 日日干狠狠操夜夜爽| 国产高潮美女av| 亚洲久久久久久中文字幕| 亚洲精品国产成人久久av| 九九在线视频观看精品| 韩国av在线不卡| 在线免费观看的www视频| 国产探花极品一区二区| 毛片女人毛片| 亚洲成人中文字幕在线播放| 国产精品野战在线观看| 久久久久国内视频| 亚洲真实伦在线观看| 亚洲图色成人| 精品久久久久久久久av| 亚洲国产精品合色在线| 乱系列少妇在线播放| 日韩欧美一区二区三区在线观看| 国产午夜精品论理片| 男插女下体视频免费在线播放| 美女大奶头视频| 免费观看精品视频网站| 免费人成视频x8x8入口观看| 午夜福利在线在线| 中文字幕人妻熟人妻熟丝袜美| 亚洲最大成人av| 亚洲第一区二区三区不卡| 亚洲av电影不卡..在线观看| 精品久久久久久久久av| 亚洲成人av在线免费| 男人狂女人下面高潮的视频| 九色成人免费人妻av| 国产毛片a区久久久久| 精品无人区乱码1区二区| av免费在线看不卡| 少妇熟女aⅴ在线视频| 三级毛片av免费| 国产精品国产三级国产av玫瑰| 久久鲁丝午夜福利片| 国产麻豆成人av免费视频| 国产精品免费一区二区三区在线| 久久鲁丝午夜福利片| 非洲黑人性xxxx精品又粗又长| 国产精品国产三级国产av玫瑰| 日韩精品中文字幕看吧| 最近的中文字幕免费完整| 一进一出好大好爽视频| 特级一级黄色大片| 亚洲国产精品成人久久小说 | 久久精品综合一区二区三区| 黄色一级大片看看| 久久人人爽人人片av| 波多野结衣高清作品| 看十八女毛片水多多多| 久久99热6这里只有精品| 长腿黑丝高跟|