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

    1959—2018年大伙房水庫流域旱澇急轉(zhuǎn)特征分析

    2019-06-28 07:16:16岳樹朋
    人民珠江 2019年6期
    關(guān)鍵詞:大伙房旱澇降雨量

    岳樹朋

    (遼寧省大伙房水庫管理局有限責(zé)任公司,遼寧 撫順 113000)

    IPCC 第五次評估報(bào)告中指出20世紀(jì)中葉以來,極端事件的強(qiáng)度和頻率發(fā)生明顯變化[1]。1961年以來,中國區(qū)域性的高溫事件、氣象干旱事件和強(qiáng)降雨事件增多。其中,干旱和洪澇災(zāi)害是最為突出的2類極端氣候事件。近年來許多學(xué)者開展了很多具有區(qū)域特色的旱澇事件研究。吳志偉等[2]通過對比5—6月和7—8月降水的差異定義了一個長周期旱澇急轉(zhuǎn)指數(shù),對長江中下游地區(qū)夏季旱澇急轉(zhuǎn)現(xiàn)象進(jìn)行研究;孫鵬等[3]分析了東江流域汛期長周期和短周期降雨旱澇急轉(zhuǎn)現(xiàn)象的趨勢變化和時空分布特征;方永旭等[4]通過旱澇急轉(zhuǎn)指數(shù)研究了王英水庫灌區(qū)水旱交替特征,并探討了旱澇急轉(zhuǎn)對農(nóng)業(yè)的影響;何慧等[5]采用長周期旱澇急轉(zhuǎn)指數(shù)研究了1961—2014年華南地區(qū)夏季旱澇急轉(zhuǎn)的演變規(guī)律,結(jié)果表明,華南地區(qū)旱澇急轉(zhuǎn)強(qiáng)度增強(qiáng),發(fā)生頻率增大;張屏等[6]通過計(jì)算各旬降水量的距平及距平百分率,在確定旱澇等級的基礎(chǔ)上對旱澇急轉(zhuǎn)事件進(jìn)行定義,并對淮北市旱澇急轉(zhuǎn)進(jìn)行規(guī)律分析;羅蔚等[7]對鄱陽湖入湖總水量變化分析,發(fā)現(xiàn)鄱陽湖旱澇和澇旱交替出現(xiàn),進(jìn)入21世紀(jì),旱澇和澇旱轉(zhuǎn)換的時間間隔變短、震蕩加??;閃麗潔等[8]分析長江中下游流域旱澇交替,發(fā)現(xiàn)流域內(nèi)旱澇急轉(zhuǎn)以澇轉(zhuǎn)旱為主,旱轉(zhuǎn)澇主要發(fā)生在長江北岸,澇轉(zhuǎn)旱主要發(fā)生在長江南岸;吳志偉等[9]揭示了華南夏季旱澇急轉(zhuǎn)事件與季節(jié)性嚴(yán)重旱澇事件的顯著差異,并對華南旱澇急轉(zhuǎn)夏季大氣環(huán)流異常特征和前期海氣異常信號進(jìn)行了分析。

    旱澇急轉(zhuǎn)指的是在研究時期內(nèi),一段時間發(fā)生干旱,一段時間發(fā)生洪澇,旱澇交替出現(xiàn)的情形[10]。上述研究表明不同流域旱澇急轉(zhuǎn)特征不盡相同,以往針對大伙房水庫流域內(nèi)旱澇異常的研究大多立足于年總降水量,對季節(jié)內(nèi)旱澇變化情況研究較少,而季節(jié)內(nèi)旱澇變化對水資源配置和工農(nóng)業(yè)生產(chǎn)十分重要。目前在水文序列分析中,對于趨勢變化及預(yù)測的方法很多,比較成熟的方法有R/S分析法、Mann-Kendall檢驗(yàn)法等。R/S分析法[11]是在20世紀(jì)中葉由英國水文學(xué)家Havold Edwin Hurst提出的一種判斷時間序列是否具有標(biāo)度不變性的統(tǒng)計(jì)分析方法。王孝禮等[12]在國內(nèi)首次將R/S分析法應(yīng)用于水文時間序列分析中,該方法能揭示水文時間序列變異特征。Mann-Kendall方法[13]廣泛應(yīng)用于檢驗(yàn)氣候與徑流的時間序列變化中,能很好的體現(xiàn)時間序列的趨勢特征?;诖耍疚母鶕?jù)大伙房水庫流域 1959—2018 年6—9月逐月降水量資料,利用長、短周期旱澇急轉(zhuǎn)指數(shù)和滑動平均法、線性趨勢法、Mann-Kendall趨勢分析法、R/S分析法從典型旱澇急轉(zhuǎn)年、旱澇急轉(zhuǎn)趨勢和旱澇急轉(zhuǎn)強(qiáng)度變化及未來趨勢等方面著手,研究大伙房水庫流域內(nèi)旱澇急轉(zhuǎn)特征,為流域內(nèi)水資源科學(xué)配置、城市防災(zāi)、減災(zāi)工作提供理論參考和數(shù)據(jù)支持。

    1 資料與方法

    1.1 研究區(qū)概況與數(shù)據(jù)資料

    大伙房水庫位于遼河的大支流渾河中上游,坐落在遼寧省撫順市境內(nèi),控制流域面積5 437 km2,總庫容21.87億m3。壩址以上流域多年平均徑流量15.13億m3,洪水主要發(fā)生在6—9月,洪水成因主要是此期間的大暴雨,壩址以上流域多年平均降水量792.1 mm。

    本文選取大伙房水庫流域(圖1)內(nèi)13個雨量站1959—2018年6—9月的實(shí)測降雨量數(shù)據(jù)進(jìn)行分析研究。

    圖1 大伙房水庫流域

    1.2 研究方法

    本文以6—9月這4個月的降水量為研究對象,在消除量綱對實(shí)驗(yàn)結(jié)果影響的基礎(chǔ)上,采用長、短周期旱澇急轉(zhuǎn)指數(shù),來分析大伙房水庫旱澇急轉(zhuǎn)的特征。長周期的變化時間尺度為 2個月,短周期則以月為單位。長周期旱澇急轉(zhuǎn)指的是,6—7 月旱,8—9 月澇稱為“旱轉(zhuǎn)澇”,反之,則為“澇轉(zhuǎn)旱”。短周期旱澇急轉(zhuǎn)指的是相鄰月份,一個月份為旱,另一個月份為澇,稱其為“旱轉(zhuǎn)澇”,反之,則為“澇轉(zhuǎn)旱”。并運(yùn)用滑動平均法、線性趨勢法、Mann-Kendall趨勢分析法和R/S分析法對旱澇急轉(zhuǎn)強(qiáng)度變化趨勢進(jìn)行分析。

    a)降雨量標(biāo)準(zhǔn)化處理

    (1)

    b)長周期旱澇急轉(zhuǎn)指標(biāo)

    LDFAL=(R89-R67)·(|R67|+|R89|)·1.8-|R67+R78|

    (2)

    式中R67——夏季6—7月標(biāo)準(zhǔn)化降雨量;R89——夏季8—9月標(biāo)準(zhǔn)化降雨量;|R67|+|R89|——旱澇強(qiáng)度項(xiàng);1.8-|R67+R89|——權(quán)重系數(shù),作用是增加長周期旱澇急轉(zhuǎn)事件所占權(quán)重,降低全旱或全澇事件權(quán)重。降雨量距平標(biāo)準(zhǔn)差小于-0.5為偏旱,小于-1為顯著偏旱,大于0.5為偏澇,大于1為顯著偏澇。

    c)短周期旱澇急轉(zhuǎn)指標(biāo)

    SDFAL=(Rj-Ri)·(|Ri+Rj|)·3.2-|Ri+Rj|,j=i+1(i=6,7,8)

    (3)

    選用3.2-|Ri+Rj|作為權(quán)重系數(shù),分別計(jì)算6—7、7—8、8—9月的月間旱澇急轉(zhuǎn)指數(shù)。

    d)Mann-Kendall趨勢分析法。Mann-Kendall法是非參數(shù)統(tǒng)計(jì)檢驗(yàn)的方法 ,目前,Mann-Kendall趨勢及突變檢驗(yàn)法對水文站實(shí)測徑流時間序列趨勢變化及突變點(diǎn)的分析已被廣泛應(yīng)用于檢驗(yàn)氣候與降雨的時間序列變化中。Mann-Kendall法能很好地揭示時間序列的趨勢變化。

    e)R/S分析法。R/S (Re-scaled Range )是一種時間序列統(tǒng)計(jì)法。利用R/S分析法可以計(jì)算得Hurst指數(shù)H,該指數(shù)可以反映具有統(tǒng)計(jì)特性的非線性數(shù)據(jù)系列的持續(xù)性,適用于降雨序列。通過Hurst 指數(shù),可判定降雨時序的狀態(tài)持續(xù)性及其記憶長度。

    2 旱澇急轉(zhuǎn)特征

    2.1 旱澇急轉(zhuǎn)典型年

    根據(jù)式(2)統(tǒng)計(jì)的長周期旱澇急轉(zhuǎn)指數(shù),選取其中最高(低)前6 a(表1),高LDFAL年中,6—7月標(biāo)準(zhǔn)化降雨量中有4 a小于-0.5(2 a小于-1.0),說明6個高LDFAL里4 a偏旱,2 a顯著偏旱,而8—9月標(biāo)準(zhǔn)化降雨量中有4 a大于0.5(1 a大于1.0,2 a大于2.0),有4 a偏澇,1 a顯著偏澇;低LDFAL年中,6—7月標(biāo)準(zhǔn)化降雨量中6 a均大于0.5(2 a大于2.0),說明6個低LDFAL里有6 a偏澇,2 a顯著偏澇,而8—9月標(biāo)準(zhǔn)化降雨量中有5 a小于-0.5(3 a小于-1.0),有5 a偏旱,3 a顯著偏旱。

    在高LDFAL年,8—9月的標(biāo)準(zhǔn)降雨量均大于6—7月,即降雨量6—7月到8—9月有增加趨勢;在低LDFAL年,8—9月的標(biāo)準(zhǔn)降雨量均小于6—7月,即降雨量6—7月到8—9月有減少趨勢,由此可知,高LDFAL年對應(yīng)“旱轉(zhuǎn)澇”年,低LDFAL年對應(yīng)“澇轉(zhuǎn)旱”年。

    表1 1959—2018年最高(低)LDFAL年及其標(biāo)準(zhǔn)化降雨量

    根據(jù)上述分析結(jié)果,在LDFAL中選出指數(shù)大于 1 的年份為典型“旱轉(zhuǎn)澇”年,小于-1的年份為典型“澇轉(zhuǎn)旱”年(表2)。需要說明的是,在低LDFAL年中,1995年8—9月的標(biāo)準(zhǔn)化降雨量為0.15,沒有達(dá)到-0.5個標(biāo)準(zhǔn)差,但6—7月標(biāo)準(zhǔn)化降雨量為3.87,大于3倍標(biāo)準(zhǔn)差,差異異常明顯,同樣的情況還有1967年;2002年6—7月的標(biāo)準(zhǔn)化降雨量為0.31,沒有達(dá)到0.5個標(biāo)準(zhǔn)差,但8—9月標(biāo)準(zhǔn)化降雨量為-1.11,大于1倍標(biāo)準(zhǔn)差,差異明顯,同樣的情況還有2003年,故將這4 a選為澇轉(zhuǎn)旱典型年;在高LDFAL年中,1964年6—7月的標(biāo)準(zhǔn)化降雨量為0.18,沒有達(dá)到-0.5個標(biāo)準(zhǔn)差,但8—9月標(biāo)準(zhǔn)化降雨量為2.02,大于2倍標(biāo)準(zhǔn)差,差異明顯;1965年8—9月的標(biāo)準(zhǔn)化降雨量為0.26,沒有達(dá)到0.5個標(biāo)準(zhǔn)差,但6—7月標(biāo)準(zhǔn)化降雨量為-1.48,大于1倍標(biāo)準(zhǔn)差,差異明顯,同樣的情況還有1972年,故將這3 a選為旱轉(zhuǎn)澇典型年。大伙房水庫流域在1959—2018年間共發(fā)生14次典型旱澇急轉(zhuǎn)現(xiàn)象。

    表2 1959—2018年旱澇急轉(zhuǎn)典型年及其標(biāo)準(zhǔn)化降雨量

    根據(jù)式(3)統(tǒng)計(jì)的短周期旱澇急轉(zhuǎn)指數(shù),再選取最高(低)的SDFAL年,得到低SDFAL年的旱澇情況。與長周期典型旱澇年選擇類似,得出典型的“旱轉(zhuǎn)澇”年(表3—5);大伙房水庫流域在1959—2018年間6—7月共發(fā)生9次典型旱澇急轉(zhuǎn)現(xiàn)象,7—8月共發(fā)生6次典型旱澇急轉(zhuǎn)現(xiàn)象,8—9月共發(fā)生11次典型旱澇急轉(zhuǎn)現(xiàn)象。

    表3 1959—2018年6—7月旱澇急轉(zhuǎn)典型年及其標(biāo)準(zhǔn)化降雨量

    表4 1959—2018年7—8月旱澇急轉(zhuǎn)典型年及其標(biāo)準(zhǔn)化降雨量

    表5 1959—2018年8、9月旱澇急轉(zhuǎn)典型年及其標(biāo)準(zhǔn)化降雨量

    2.2 旱澇急轉(zhuǎn)典型年降雨特征

    旱澇急轉(zhuǎn)年分為“旱轉(zhuǎn)澇”年和“澇轉(zhuǎn)旱”年,根據(jù)長周期旱澇急轉(zhuǎn)指數(shù)的計(jì)算結(jié)果,1987年LDFAL指數(shù)為3.15,屬于典型“旱轉(zhuǎn)澇”年,1975年LDFAL指數(shù)為-7.29,屬于典型“澇轉(zhuǎn)旱”年,繪制1975—1987年6—9月逐日降雨量(圖2)。

    圖2 大伙房流域長周期典型旱澇急轉(zhuǎn)年逐日降雨量

    在1975年,6—7月降雨量為553.4 mm,8—9月降雨量為131.9 mm,6—7月降雨次數(shù)多,其中7月發(fā)生6次大降雨,最大值達(dá)92.7 mm;8—9月降雨量與降雨次數(shù)均少,超過10 mm降雨過程僅5次,最大值為20.6 mm,因此認(rèn)為7月以澇為主,8—9月為持續(xù)性干旱,對應(yīng)低LDFAL年的澇轉(zhuǎn)旱過程。

    在1987年,6—7月降雨量為240.9 mm,8—9月降雨量為344.1 mm,6月降雨量少,且次降雨量較小,大多未能形成有效降雨,7月降雨較6月有所增加,但總體雨量較少;8—9月降雨量與降雨次數(shù)增多,8月發(fā)生5次大的降雨過程,9月中后期降雨量減少,總體上,6—9月降雨過程較為溫和,無極端降雨過程,因此認(rèn)為6月為持續(xù)干旱,7月為間歇性干旱,8—9月以澇為主,對應(yīng)高LDFAL年的旱轉(zhuǎn)澇過程。

    根據(jù)短周期旱澇急轉(zhuǎn)指數(shù)的計(jì)算結(jié)果,繪制6—7、7—8、8—9月典型旱澇急轉(zhuǎn)年逐日降雨量(圖3)。

    a)6—7月

    b)7—8月

    c)8—9月圖3 大伙房流域短周期典型旱澇急轉(zhuǎn)年逐日降雨量

    2010年,6月平均降雨67.2 mm,降雨稀少,僅發(fā)生2次大于10 mm降雨過程,且持續(xù)多日無雨,表現(xiàn)為持續(xù)干旱;7月降雨增多,平均降雨386.5 mm,在7月下旬發(fā)生4次降雨過程,最大1 d降雨量105.1 mm,7月以澇為主且雨量集中,對應(yīng)高SDFAL6—7月的旱轉(zhuǎn)澇過程。

    2014年,6月平均降雨187.8 mm,上旬降雨稀少,中、下旬雨量增多,總體以澇為主,7月平均降雨78.1 mm,上、下旬幾乎無降雨過程,中旬降雨有所增加,但雨量小,除21日降雨28.9 mm外,其余降雨均不超過10 mm,表現(xiàn)為持續(xù)干旱,對應(yīng)低SDFAL6—7月的澇轉(zhuǎn)旱過程。

    2011年,7月平均降雨83.4 mm,降雨稀少,僅發(fā)生2次大于10 mm降雨過程,且持續(xù)多日無雨,表現(xiàn)為持續(xù)干旱;8月降雨增多,平均降雨218.4 mm,在8月下旬發(fā)生3次大的降雨過程,其余時間降雨稀少,8月以澇為主,雨量相對集中,對應(yīng)高SDFAL7—8月的旱轉(zhuǎn)澇過程。

    1991年,7月平均降雨292.8 mm,降雨次數(shù)多,雨量大,以澇為主,8月平均降雨87.9 mm,僅有3次大于10 mm的降雨,且最大值僅為20.3 mm,表現(xiàn)為間歇性干旱,對應(yīng)低SDFAL7—8月的澇轉(zhuǎn)旱過程。

    1992年,8月平均降雨56.6mm,降雨稀少,僅發(fā)生2次大于10 mm降雨過程,且持續(xù)多日無雨,表現(xiàn)為持續(xù)干旱;9月降雨增多,但總體雨量較少,平均降雨104.2 mm,在9月發(fā)生3次大的降雨過程,其余時間降雨稀少,9月相對較澇,對應(yīng)高SDFAL8—9月的旱轉(zhuǎn)澇過程。

    2005年,8月平均降雨338.2 mm,上、中旬降雨次數(shù)多,雨量大且降雨集中,以澇為主,9月平均降雨24.0 mm,降雨次數(shù)少且無大于10 mm的降雨過程,表現(xiàn)為持續(xù)干旱,對應(yīng)低SDFAL8—9月的澇轉(zhuǎn)旱過程。

    綜上所述,長、短周期旱澇急轉(zhuǎn)指數(shù)可以很好地反映大伙房水庫流域 6—9 月旱澇急轉(zhuǎn)變化的特征。

    3 旱澇急轉(zhuǎn)變化分析

    3.1 旱澇急轉(zhuǎn)變化趨勢分析

    LDFAL可以反映研究區(qū)域是“旱轉(zhuǎn)澇”,還是“澇轉(zhuǎn)旱”的過程,繪制出 1959—2018 年大伙房流域LDFAL指數(shù)時間序列(圖4)。大伙房水庫流域在1959—2018年期間,發(fā)生典型“旱轉(zhuǎn)澇”5 a,平均12 a/次,發(fā)生典型“澇轉(zhuǎn)旱”9 a,平均6.67 a/次,旱澇急轉(zhuǎn)發(fā)生頻率為4.28 a/次;且在1987年后沒有“旱轉(zhuǎn)澇”發(fā)生,2005年以后無典型旱澇急轉(zhuǎn)現(xiàn)象發(fā)生,曲線起伏平緩,LDFAL變化趨勢呈不明顯的微弱上升。由上所述,隨著時間推移,LDFAL逐漸趨于平緩,長周期旱澇急轉(zhuǎn)現(xiàn)象發(fā)生概率降低。

    圖4 1959—2018年大伙房流域LDFAL指數(shù)時間序列

    根據(jù)SDFAL計(jì)算結(jié)果繪制出 1959—2018 年大伙房流域6—7、7—8、8—9月的SDFAL指數(shù)時間序列(圖5)。大伙房水庫流域在1959—2018年期間,6—7月發(fā)生典型“旱轉(zhuǎn)澇”6 a,平均10 a/次,發(fā)生典型“澇轉(zhuǎn)旱”3 a,平均20 a/次,旱澇急轉(zhuǎn)發(fā)生頻率為6.67 a/次,且在1965—2009年間SDFAL相對平緩,發(fā)生“旱轉(zhuǎn)澇”和“澇轉(zhuǎn)旱”各一次,隨著時間推移,SDFAL變化呈下降趨勢,即由“旱轉(zhuǎn)澇”主導(dǎo)轉(zhuǎn)向以“澇轉(zhuǎn)旱”為主導(dǎo);7—8月發(fā)生典型“旱轉(zhuǎn)澇”2 a,平均30 a/次,發(fā)生典型“澇轉(zhuǎn)旱”4 a,平均15 a/次,旱澇急轉(zhuǎn)發(fā)生頻率為10 a/次;且在1993年以后SDFAL變化相對平緩,隨著時間推移,SDFAL變化呈上升趨勢,即由“澇轉(zhuǎn)旱”主導(dǎo)轉(zhuǎn)向以“旱轉(zhuǎn)澇”為主導(dǎo);8—9月發(fā)生典型“旱轉(zhuǎn)澇”7 a,平均8.57 a/次,8—9月發(fā)生典型“澇轉(zhuǎn)旱”4 a,平均15 a/次,旱澇急轉(zhuǎn)發(fā)生頻率為5.45 a/次;且在1993年之前未發(fā)生“澇轉(zhuǎn)旱”,1993年以后未發(fā)生“旱轉(zhuǎn)澇”,隨著時間推移,SDFAL變化呈下降趨勢,即由“旱轉(zhuǎn)澇”主導(dǎo)轉(zhuǎn)向以“澇轉(zhuǎn)旱”為主導(dǎo)。

    a)6—7月

    b)7—8月

    c)8—9月圖5 1959—2018年大伙房流域SDFAL指數(shù)時間序列

    3.2 旱澇急轉(zhuǎn)強(qiáng)度變化趨勢分析

    LDFAL可以反映研究區(qū)域的旱澇急轉(zhuǎn)過程,但是,在一定程度上未能反映旱澇急轉(zhuǎn)的強(qiáng)度,因此,為了更好地反映旱澇急轉(zhuǎn)的強(qiáng)度變化,對LDFAL和SDFAL取絕對值,并繪制出 1959—2018 年旱澇急轉(zhuǎn)強(qiáng)度的年際變化(圖6、7)。

    圖6 1959—2018年大伙房流域LDFAL強(qiáng)度時間序列

    滑動平均法可以使序列高頻震蕩對變化趨勢分析的影響得以弱化[14],本文旱澇急轉(zhuǎn)強(qiáng)度以5 a進(jìn)行滑動平均,觀察圖6,LDFAL強(qiáng)度序列呈下降趨勢,表示旱澇轉(zhuǎn)換強(qiáng)度減弱。在60 a內(nèi),LDFAL強(qiáng)度序列呈現(xiàn)階段性,1972—1982 年、1985—1994年為相對偏強(qiáng)期;其中,1975、1979 年旱澇急轉(zhuǎn)情況比較嚴(yán)重。

    觀察圖7,6—7月SDFAL強(qiáng)度序列呈上升趨勢,旱澇轉(zhuǎn)換強(qiáng)度增強(qiáng),7—8月SDFAL強(qiáng)度序列呈下降趨勢,旱澇轉(zhuǎn)換強(qiáng)度減弱,8—9月SDFAL強(qiáng)度序列呈不明顯下降趨勢,旱澇轉(zhuǎn)換強(qiáng)度減弱。

    a)6—7月

    b)7—8月圖7 1959—2018年大伙房流域SDFAL強(qiáng)度時間序列

    c)8—9月續(xù)圖7 1959—2018年大伙房流域SDFAL強(qiáng)度時間序列

    利用Mann-Kendall趨勢分析法,計(jì)算統(tǒng)計(jì)量U值,如果U>0,表明有上升趨勢,如果U<0,表明有下降趨勢。當(dāng)U大于置信區(qū)間時,即U>U0.05/2=1.96,表示序列趨勢變化顯著。

    對大伙房水庫流域的旱澇急轉(zhuǎn)強(qiáng)度計(jì)算結(jié)果見表6,長周期旱澇急轉(zhuǎn)強(qiáng)度呈下降趨勢,短周期6—7月旱澇急轉(zhuǎn)強(qiáng)度呈上升趨勢,7—8月旱澇急轉(zhuǎn)強(qiáng)度呈下降趨勢,8—9月旱澇急轉(zhuǎn)強(qiáng)度呈下降趨勢。

    表6 旱澇急轉(zhuǎn)強(qiáng)度Mann-Kendall統(tǒng)計(jì)值成果

    利用Mann-Kendall法揭示序列的趨勢特征,R/S分析法揭示序列的持續(xù)性[15],結(jié)合2種方法綜合得出未來的趨勢特征[16]。Hurst系數(shù)H的取值范圍為(0,1),當(dāng)H=0.5時,旱澇急轉(zhuǎn)強(qiáng)度序列為隨機(jī)序列,目前的強(qiáng)度不會影響未來趨勢;當(dāng)H>0.5時,強(qiáng)度序列具有持續(xù)性,存在長期記憶性,即下一個狀態(tài)將持續(xù)上一個狀態(tài)的態(tài)勢,且其記憶性不隨時間標(biāo)度而變化,當(dāng)H<0.5時,系統(tǒng)是逆持續(xù)性的,即下一個狀態(tài)與上一個狀態(tài)的態(tài)勢相反,這種逆持續(xù)性行為的強(qiáng)度取決于距離0的遠(yuǎn)近,H距離0值越近,負(fù)相關(guān)性越顯著。

    對大伙房水庫流域旱澇急轉(zhuǎn)強(qiáng)度未來趨勢預(yù)測,計(jì)算結(jié)果見表7。在2018年后的一定時間內(nèi),大伙房水庫流域長周期旱澇急轉(zhuǎn)強(qiáng)度呈現(xiàn)不明顯的下降趨勢,短周期6—7、7—8月的旱澇急轉(zhuǎn)強(qiáng)度呈現(xiàn)不明顯的上升趨勢,8—9月的旱澇急轉(zhuǎn)強(qiáng)度呈現(xiàn)不明顯的下降趨勢。

    表7 旱澇急轉(zhuǎn)強(qiáng)度未來趨勢變化特征

    4 結(jié)論

    大伙房水庫流域降水量年內(nèi)分配差異明顯,主要集中在6—9月份。流域處于盛行西風(fēng)帶,降雨特征受歐亞大氣環(huán)流影響較大,環(huán)流形勢和太平洋副熱帶高壓相互作用,常形成大暴雨;其次,臺風(fēng)、華北氣旋等也是流域降雨的重要成因。流域的降雨系統(tǒng)主要有西南系統(tǒng)(華北氣旋、渤海氣旋、江淮氣旋、臺風(fēng)、倒槽)、西北系統(tǒng)(冷鋒、蒙古氣旋、蒙古低壓)、東北系統(tǒng)(東北氣旋、東北低壓)。隨著全球氣溫升高,使流域大氣環(huán)流、降雨系統(tǒng)發(fā)生變化,極端氣候現(xiàn)象增多,導(dǎo)致流域降雨特征發(fā)生變化,進(jìn)而影響流域旱澇急轉(zhuǎn)規(guī)律。

    本文基于對大伙房水庫流域內(nèi)13個雨量站降雨資料的分析,得到流域旱澇急轉(zhuǎn)變化的特點(diǎn),有以下幾個重要結(jié)論。

    a)長、短周期旱澇急轉(zhuǎn)指數(shù)可以很好地反映大伙房水庫流域 6—9 月旱澇急轉(zhuǎn)的特征。

    b)大伙房水庫流域長周期旱澇急轉(zhuǎn)發(fā)生頻率為4.28 a/次,LDFAL指數(shù)1992年后變化平緩,旱澇急轉(zhuǎn)情況不突出,短周期6—7、7—8、8—9月旱澇急轉(zhuǎn)發(fā)生頻率依次為6.67、10.00、5.45 a/次。

    c)長周期LDFAL變化逐漸趨于平緩,旱澇急轉(zhuǎn)現(xiàn)象發(fā)生概率降低;短周期6—7、8—9月SDFAL變化呈下降趨勢,由“旱轉(zhuǎn)澇”主導(dǎo)轉(zhuǎn)向以“澇轉(zhuǎn)旱”主導(dǎo);7—8月SDFAL變化呈上升趨勢,由“澇轉(zhuǎn)旱”主導(dǎo)轉(zhuǎn)向以“旱轉(zhuǎn)澇”為主導(dǎo)。

    d)2018年后的一定時間內(nèi),大伙房水庫流域長周期旱澇急轉(zhuǎn)強(qiáng)度呈現(xiàn)不明顯的下降趨勢,短周期6—7月、7—8月的旱澇急轉(zhuǎn)強(qiáng)度呈現(xiàn)不明顯的上升趨勢,8—9月的旱澇急轉(zhuǎn)強(qiáng)度呈現(xiàn)不明顯的下降趨勢。

    本文僅從降水的角度對大伙房水庫流域內(nèi)旱澇急轉(zhuǎn)現(xiàn)象進(jìn)行研究,總結(jié)歸納了流域內(nèi)旱澇急轉(zhuǎn)變化特征,分析了旱澇急轉(zhuǎn)變化規(guī)律及未來趨勢,具有一定的研究意義,但未考慮溫度、大氣環(huán)流等因素的約束,后續(xù)將結(jié)合多因素對研究區(qū)旱澇急轉(zhuǎn)規(guī)律做進(jìn)一步的探討、分析。

    猜你喜歡
    大伙房旱澇降雨量
    降雨量與面積的關(guān)系
    四川省1960-2017年旱澇時空變化特征分析
    大伙房采風(fēng)作品選登
    中華詩詞(2018年12期)2018-03-25 13:46:36
    遼寧省大伙房輸水工程TBM預(yù)制管片襯砌方案及施工
    大伙房輸水隧洞TBM法施工風(fēng)險(xiǎn)識別及對策
    氣候變化背景下西北干旱區(qū)旱澇的變化規(guī)律
    洞庭湖區(qū)降雨特性分析
    北方農(nóng)牧交錯帶氣候變化與旱澇響應(yīng)特征
    寧夏旱澇災(zāi)害對氣候變化的響應(yīng)
    羅甸縣各鄉(xiāng)鎮(zhèn)實(shí)測降雨量分析及應(yīng)用研究
    女人久久www免费人成看片| 欧美精品一区二区免费开放| 亚洲美女黄色视频免费看| 亚洲天堂av无毛| 十分钟在线观看高清视频www | 日本欧美国产在线视频| 国产一区有黄有色的免费视频| 亚洲国产日韩一区二区| 少妇丰满av| 狂野欧美激情性xxxx在线观看| 日韩欧美精品免费久久| 中文字幕人妻熟人妻熟丝袜美| 麻豆成人av视频| 亚洲一区二区三区欧美精品| 青春草亚洲视频在线观看| 在线 av 中文字幕| 精品一品国产午夜福利视频| 久久久亚洲精品成人影院| 国产成人精品婷婷| 精品人妻熟女av久视频| av专区在线播放| 三级国产精品欧美在线观看| 久久人人爽人人爽人人片va| 久久国内精品自在自线图片| 成人综合一区亚洲| 精品人妻偷拍中文字幕| 国产亚洲午夜精品一区二区久久| 亚洲精品日本国产第一区| 少妇精品久久久久久久| 夜夜看夜夜爽夜夜摸| 欧美精品国产亚洲| 曰老女人黄片| 久久久欧美国产精品| 国产永久视频网站| 我要看日韩黄色一级片| 国产精品久久久久久精品古装| 熟女av电影| 日日啪夜夜撸| 国产黄片视频在线免费观看| 欧美成人午夜免费资源| 欧美区成人在线视频| 少妇人妻 视频| 久久6这里有精品| 免费看日本二区| 亚洲在久久综合| av播播在线观看一区| 亚洲情色 制服丝袜| 久久精品国产亚洲网站| 国产毛片在线视频| 26uuu在线亚洲综合色| 久久99热这里只频精品6学生| 国产69精品久久久久777片| 亚洲精华国产精华液的使用体验| 欧美精品一区二区免费开放| 亚洲色图综合在线观看| 日日啪夜夜爽| 纵有疾风起免费观看全集完整版| 成人国产麻豆网| 一级毛片电影观看| 国产精品99久久久久久久久| 女人久久www免费人成看片| 天堂中文最新版在线下载| 91aial.com中文字幕在线观看| 欧美日韩在线观看h| 国内揄拍国产精品人妻在线| 亚洲综合色惰| 一级毛片电影观看| 久久精品久久精品一区二区三区| 中文字幕av电影在线播放| 成人无遮挡网站| 女性被躁到高潮视频| 色视频www国产| 久久精品久久久久久久性| 久久精品国产亚洲av天美| 日本欧美视频一区| 欧美97在线视频| 亚洲,欧美,日韩| 欧美xxⅹ黑人| 亚洲欧美成人综合另类久久久| 午夜久久久在线观看| 黄色毛片三级朝国网站 | 人妻制服诱惑在线中文字幕| 日韩不卡一区二区三区视频在线| 老司机亚洲免费影院| 久久久国产欧美日韩av| 国产精品一区二区性色av| 亚洲国产欧美在线一区| 国产精品国产三级国产专区5o| a级片在线免费高清观看视频| 日本av免费视频播放| 黄色怎么调成土黄色| 香蕉精品网在线| 五月开心婷婷网| 搡女人真爽免费视频火全软件| 久久午夜福利片| 自线自在国产av| 日韩熟女老妇一区二区性免费视频| 免费高清在线观看视频在线观看| 亚洲av成人精品一区久久| 国产爽快片一区二区三区| 成人美女网站在线观看视频| 亚洲精华国产精华液的使用体验| 99re6热这里在线精品视频| 国产黄片美女视频| 日韩欧美 国产精品| 久久久久久久久久人人人人人人| 国产永久视频网站| 免费观看a级毛片全部| 十分钟在线观看高清视频www | 亚洲av男天堂| 午夜激情久久久久久久| av一本久久久久| 免费观看无遮挡的男女| 简卡轻食公司| 如何舔出高潮| 日韩视频在线欧美| 看免费成人av毛片| 日韩人妻高清精品专区| 看非洲黑人一级黄片| 有码 亚洲区| 日韩免费高清中文字幕av| 免费黄色在线免费观看| 亚洲不卡免费看| 亚洲精华国产精华液的使用体验| 国产女主播在线喷水免费视频网站| 人人妻人人澡人人看| 国产亚洲91精品色在线| 高清av免费在线| 少妇人妻 视频| 精品久久久久久电影网| av一本久久久久| 久久国内精品自在自线图片| 国产精品一二三区在线看| 精品亚洲成国产av| a级片在线免费高清观看视频| 99视频精品全部免费 在线| 久久免费观看电影| 少妇人妻 视频| av天堂久久9| 水蜜桃什么品种好| 大又大粗又爽又黄少妇毛片口| av天堂中文字幕网| 久久久亚洲精品成人影院| 日本av免费视频播放| 免费观看a级毛片全部| 日本午夜av视频| 久久久久久久久大av| 国产一区亚洲一区在线观看| 亚洲精品国产成人久久av| 亚洲婷婷狠狠爱综合网| 欧美日韩av久久| 内地一区二区视频在线| 成年女人在线观看亚洲视频| .国产精品久久| 哪个播放器可以免费观看大片| 久久久久久久久大av| 色视频在线一区二区三区| 欧美精品亚洲一区二区| 国产一区有黄有色的免费视频| 久久国产精品大桥未久av | 免费观看的影片在线观看| 纯流量卡能插随身wifi吗| 免费看光身美女| 国产亚洲5aaaaa淫片| 在线观看三级黄色| 六月丁香七月| 国产91av在线免费观看| 久久鲁丝午夜福利片| 久久国产乱子免费精品| 一二三四中文在线观看免费高清| 亚洲av综合色区一区| 国产又色又爽无遮挡免| av在线观看视频网站免费| 国产精品秋霞免费鲁丝片| 黑人巨大精品欧美一区二区蜜桃 | 欧美bdsm另类| 亚洲,一卡二卡三卡| 丝瓜视频免费看黄片| 久久精品国产亚洲av涩爱| 亚洲精品aⅴ在线观看| 日韩中字成人| 成人影院久久| 麻豆乱淫一区二区| 秋霞伦理黄片| 亚洲av国产av综合av卡| 亚洲内射少妇av| 天美传媒精品一区二区| 亚洲精品视频女| 中文精品一卡2卡3卡4更新| 亚洲欧美日韩卡通动漫| 秋霞在线观看毛片| 久久精品国产亚洲av天美| av有码第一页| 26uuu在线亚洲综合色| 国产视频首页在线观看| 一区二区三区乱码不卡18| 国产乱人偷精品视频| 大香蕉久久网| 免费观看a级毛片全部| 偷拍熟女少妇极品色| 人妻少妇偷人精品九色| 国产免费又黄又爽又色| 欧美区成人在线视频| 99久久精品热视频| 91精品伊人久久大香线蕉| 男人狂女人下面高潮的视频| 99国产精品免费福利视频| 国产 一区精品| 青春草亚洲视频在线观看| 国产 精品1| 亚洲图色成人| 在线播放无遮挡| 亚洲精品aⅴ在线观看| 亚洲伊人久久精品综合| 午夜老司机福利剧场| 中文字幕制服av| 亚洲无线观看免费| 交换朋友夫妻互换小说| 免费不卡的大黄色大毛片视频在线观看| 男女国产视频网站| 国产 精品1| 久久精品熟女亚洲av麻豆精品| 成人毛片60女人毛片免费| 国产精品国产三级国产av玫瑰| 嫩草影院入口| 乱人伦中国视频| 国产高清不卡午夜福利| 一级毛片 在线播放| 人妻人人澡人人爽人人| av线在线观看网站| 国产男女超爽视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人手机| 在线 av 中文字幕| 99久久人妻综合| 在线精品无人区一区二区三| 九九在线视频观看精品| 亚洲av.av天堂| 一区二区三区免费毛片| 乱码一卡2卡4卡精品| 成人毛片a级毛片在线播放| 丰满乱子伦码专区| 亚洲精品成人av观看孕妇| 伊人久久精品亚洲午夜| 久久久久久久久大av| 曰老女人黄片| 人妻制服诱惑在线中文字幕| 日韩三级伦理在线观看| 国产中年淑女户外野战色| 国产探花极品一区二区| 99久久精品热视频| 国产精品熟女久久久久浪| 免费黄网站久久成人精品| 成人免费观看视频高清| 中文欧美无线码| 亚洲精品国产色婷婷电影| 狂野欧美激情性bbbbbb| 一级毛片aaaaaa免费看小| 寂寞人妻少妇视频99o| 亚洲精品色激情综合| 国内精品宾馆在线| 中文字幕久久专区| 伊人亚洲综合成人网| 噜噜噜噜噜久久久久久91| 建设人人有责人人尽责人人享有的| 精品人妻熟女毛片av久久网站| 少妇丰满av| 国产精品成人在线| 亚洲国产精品一区三区| 热re99久久精品国产66热6| 亚洲中文av在线| 国产成人免费观看mmmm| 夜夜骑夜夜射夜夜干| 久久婷婷青草| 久久99蜜桃精品久久| 一二三四中文在线观看免费高清| 又黄又爽又刺激的免费视频.| 在线看a的网站| 伦理电影免费视频| 国产视频内射| 免费大片黄手机在线观看| 九九久久精品国产亚洲av麻豆| 国产av国产精品国产| 特大巨黑吊av在线直播| 免费大片黄手机在线观看| 美女主播在线视频| 高清不卡的av网站| 纯流量卡能插随身wifi吗| 波野结衣二区三区在线| 日日撸夜夜添| 久久久久网色| 少妇的逼水好多| 免费观看a级毛片全部| 丰满少妇做爰视频| 日韩精品免费视频一区二区三区 | 国精品久久久久久国模美| 久久久国产欧美日韩av| 日韩av不卡免费在线播放| 午夜免费鲁丝| 国产日韩欧美亚洲二区| 久久久久久久久久久久大奶| 久久精品久久精品一区二区三区| 99热这里只有是精品50| 这个男人来自地球电影免费观看 | 桃花免费在线播放| 曰老女人黄片| a级片在线免费高清观看视频| 深夜a级毛片| 精品国产国语对白av| 十八禁高潮呻吟视频 | 久久久久视频综合| 亚洲精品日本国产第一区| 久久毛片免费看一区二区三区| 中文天堂在线官网| 一级二级三级毛片免费看| 国产91av在线免费观看| 大话2 男鬼变身卡| 丰满迷人的少妇在线观看| 免费黄网站久久成人精品| 亚洲国产精品一区二区三区在线| 美女xxoo啪啪120秒动态图| 国产精品久久久久久久久免| 久久久亚洲精品成人影院| 老司机影院成人| 夜夜看夜夜爽夜夜摸| 午夜91福利影院| 少妇高潮的动态图| 热re99久久国产66热| 日本黄色日本黄色录像| 国产精品女同一区二区软件| 久久久欧美国产精品| √禁漫天堂资源中文www| 青春草视频在线免费观看| 国产精品国产三级国产av玫瑰| 久久精品国产鲁丝片午夜精品| 精品国产乱码久久久久久小说| 国产精品久久久久久精品古装| 在线观看免费日韩欧美大片 | 两个人免费观看高清视频 | 一个人看视频在线观看www免费| 精品一区在线观看国产| 91精品一卡2卡3卡4卡| 丰满乱子伦码专区| 欧美老熟妇乱子伦牲交| 在线观看美女被高潮喷水网站| 日韩成人伦理影院| 久久久精品94久久精品| 在线观看免费视频网站a站| 又粗又硬又长又爽又黄的视频| 中文字幕亚洲精品专区| 我的老师免费观看完整版| 中文字幕av电影在线播放| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲5aaaaa淫片| 啦啦啦视频在线资源免费观看| 中文字幕免费在线视频6| .国产精品久久| 国产亚洲91精品色在线| 色视频www国产| 亚洲中文av在线| 国产国拍精品亚洲av在线观看| 亚洲精品第二区| 久久久国产精品麻豆| 自拍欧美九色日韩亚洲蝌蚪91 | 人妻夜夜爽99麻豆av| 亚洲欧美日韩卡通动漫| 成年女人在线观看亚洲视频| 亚洲欧美成人综合另类久久久| 少妇人妻精品综合一区二区| 插逼视频在线观看| 国产精品女同一区二区软件| 国产成人精品一,二区| 美女主播在线视频| 人妻系列 视频| 国产有黄有色有爽视频| 午夜91福利影院| 黄色日韩在线| 天美传媒精品一区二区| 国产av一区二区精品久久| 成人特级av手机在线观看| 久久久欧美国产精品| av天堂久久9| 国产毛片在线视频| 国精品久久久久久国模美| 黄色一级大片看看| 大片免费播放器 马上看| 久久国产精品大桥未久av | 午夜日本视频在线| 女性被躁到高潮视频| 亚洲电影在线观看av| av在线播放精品| 亚洲人成网站在线观看播放| 18禁裸乳无遮挡动漫免费视频| 在线亚洲精品国产二区图片欧美 | 国产精品一区二区三区四区免费观看| 久久久久视频综合| 亚洲国产欧美日韩在线播放 | 丰满迷人的少妇在线观看| 91成人精品电影| 国产一区二区在线观看日韩| 国产精品久久久久久精品电影小说| 涩涩av久久男人的天堂| 免费高清在线观看视频在线观看| 久久午夜福利片| 麻豆成人av视频| 日本猛色少妇xxxxx猛交久久| 黑人巨大精品欧美一区二区蜜桃 | 久久婷婷青草| 亚洲av福利一区| 99久久精品一区二区三区| 亚洲精品国产色婷婷电影| 日韩成人伦理影院| 久久 成人 亚洲| 我的老师免费观看完整版| 国精品久久久久久国模美| 久久精品熟女亚洲av麻豆精品| 国产精品偷伦视频观看了| 少妇猛男粗大的猛烈进出视频| 成人美女网站在线观看视频| 成人黄色视频免费在线看| 乱系列少妇在线播放| 精品国产露脸久久av麻豆| 高清不卡的av网站| 日韩一区二区视频免费看| 欧美日韩精品成人综合77777| 下体分泌物呈黄色| 亚洲成人av在线免费| 欧美日韩视频高清一区二区三区二| 免费人成在线观看视频色| 桃花免费在线播放| 在线亚洲精品国产二区图片欧美 | 黄片无遮挡物在线观看| 亚洲精品日本国产第一区| 一级毛片黄色毛片免费观看视频| 国产高清有码在线观看视频| 午夜免费男女啪啪视频观看| 欧美精品国产亚洲| 亚洲不卡免费看| 边亲边吃奶的免费视频| 91在线精品国自产拍蜜月| 黄色欧美视频在线观看| 国产亚洲5aaaaa淫片| 性色avwww在线观看| 人人妻人人爽人人添夜夜欢视频 | 能在线免费看毛片的网站| 80岁老熟妇乱子伦牲交| 日韩制服骚丝袜av| 成人二区视频| 2022亚洲国产成人精品| 麻豆成人午夜福利视频| 你懂的网址亚洲精品在线观看| 久久韩国三级中文字幕| 伦精品一区二区三区| 我要看日韩黄色一级片| 成人综合一区亚洲| 我要看日韩黄色一级片| 色94色欧美一区二区| 国产精品熟女久久久久浪| 一区二区三区乱码不卡18| 我要看黄色一级片免费的| 人人妻人人澡人人爽人人夜夜| 综合色丁香网| 国产毛片在线视频| 中国国产av一级| 蜜桃在线观看..| 国产亚洲午夜精品一区二区久久| 亚洲精品国产av蜜桃| 亚洲av中文av极速乱| 精品卡一卡二卡四卡免费| 97超碰精品成人国产| 免费播放大片免费观看视频在线观看| 下体分泌物呈黄色| 高清午夜精品一区二区三区| 午夜免费鲁丝| 久久婷婷青草| 日本vs欧美在线观看视频 | 边亲边吃奶的免费视频| 亚洲四区av| 最近中文字幕2019免费版| 91精品伊人久久大香线蕉| 性高湖久久久久久久久免费观看| 一级二级三级毛片免费看| 日日摸夜夜添夜夜添av毛片| 男女国产视频网站| 在线观看免费视频网站a站| av又黄又爽大尺度在线免费看| 国产精品久久久久久久久免| 成人毛片60女人毛片免费| 成人特级av手机在线观看| 亚洲人成网站在线播| 日韩av在线免费看完整版不卡| 精品久久久精品久久久| 欧美激情极品国产一区二区三区 | 九九在线视频观看精品| videossex国产| 久久久久久久久久久免费av| 91久久精品国产一区二区成人| 看免费成人av毛片| 国产黄片视频在线免费观看| 久久精品国产亚洲av涩爱| 嫩草影院新地址| 免费大片黄手机在线观看| 久久6这里有精品| 欧美日本中文国产一区发布| 黄色怎么调成土黄色| 中国美白少妇内射xxxbb| 久久久精品94久久精品| 亚洲不卡免费看| 国内少妇人妻偷人精品xxx网站| 六月丁香七月| 99热这里只有是精品50| 免费高清在线观看视频在线观看| 亚洲va在线va天堂va国产| 美女脱内裤让男人舔精品视频| 精品久久久久久久久av| 免费少妇av软件| 女人精品久久久久毛片| 在线观看人妻少妇| 日日撸夜夜添| 国产色爽女视频免费观看| 国产欧美另类精品又又久久亚洲欧美| 又黄又爽又刺激的免费视频.| 国产欧美另类精品又又久久亚洲欧美| 老女人水多毛片| 我的女老师完整版在线观看| 高清在线视频一区二区三区| 日日撸夜夜添| 伊人久久国产一区二区| 国产黄频视频在线观看| 久久精品国产亚洲网站| 在线看a的网站| 欧美日韩国产mv在线观看视频| 国产成人freesex在线| 成人漫画全彩无遮挡| 亚洲性久久影院| 久久人人爽人人爽人人片va| 午夜影院在线不卡| 国产精品人妻久久久影院| 观看免费一级毛片| 一本—道久久a久久精品蜜桃钙片| 在线观看人妻少妇| 国产男人的电影天堂91| 日韩一区二区视频免费看| 精品一区在线观看国产| 熟女av电影| 国产黄片视频在线免费观看| 国产成人精品福利久久| 亚洲欧美日韩卡通动漫| 九草在线视频观看| 久久国产亚洲av麻豆专区| 另类亚洲欧美激情| 国产亚洲午夜精品一区二区久久| 老熟女久久久| 最黄视频免费看| 嘟嘟电影网在线观看| 国产午夜精品久久久久久一区二区三区| 久热久热在线精品观看| 69精品国产乱码久久久| 在线观看www视频免费| 欧美精品一区二区免费开放| 91精品国产九色| 久久精品久久久久久噜噜老黄| 国产乱人偷精品视频| 国产精品人妻久久久影院| 日韩一区二区视频免费看| 99九九线精品视频在线观看视频| 伊人久久精品亚洲午夜| 精品少妇内射三级| 亚洲电影在线观看av| 欧美日韩av久久| 男女边摸边吃奶| 妹子高潮喷水视频| 热re99久久国产66热| 国产成人一区二区在线| 中文字幕av电影在线播放| 国产在线免费精品| 高清不卡的av网站| 成人亚洲精品一区在线观看| 嫩草影院入口| 久久久久久久久久久免费av| 日本爱情动作片www.在线观看| 97超碰精品成人国产| 久久久久久久精品精品| 国产成人freesex在线| 国产免费视频播放在线视频| 性高湖久久久久久久久免费观看| 六月丁香七月| 欧美xxⅹ黑人| 午夜激情福利司机影院| av福利片在线| 人人妻人人看人人澡| 国模一区二区三区四区视频| 永久免费av网站大全| 天堂俺去俺来也www色官网| 久久久久精品久久久久真实原创| 在线观看人妻少妇| 国产精品免费大片| 80岁老熟妇乱子伦牲交| 国产精品一区二区在线不卡| 国产毛片在线视频| 91在线精品国自产拍蜜月| 国产高清不卡午夜福利| 美女cb高潮喷水在线观看| 亚洲av欧美aⅴ国产| 99热这里只有是精品在线观看| 少妇丰满av| 国产亚洲av片在线观看秒播厂| 午夜激情久久久久久久| 日韩在线高清观看一区二区三区| 91在线精品国自产拍蜜月| 国产高清不卡午夜福利| 丝袜脚勾引网站| 亚洲精品自拍成人| 丝袜喷水一区| 91精品伊人久久大香线蕉|