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

    基于特征分析的設計降雨過程分析與計算

    2022-09-20 01:26:26龐婷婷柳金杰呂樹峰滿子豪
    河南科學 2022年7期
    關鍵詞:雨型場次歷時

    龐婷婷, 柳金杰, 呂樹峰, 滿子豪, 董 森, 何 帥

    (1.河北工程大學河北省智慧水利重點實驗室,邯鄲 056038; 2.河北省邯鄲市水文勘測研究中心,邯鄲 056001;3.邯鄲市漳滏河灌溉供水管理處,邯鄲 056001; 4.上海交通大學設計學院,上海 200204)

    設計降雨過程是工程水文計算的重要部分,降雨過程的客觀與否會直接影響區(qū)域工程水文計算的準確性,進而會影響工程的運行管理.2021年7月20日,河南鄭州發(fā)生罕見特大暴雨災害,造成了重大的人員傷亡和財產損失,引發(fā)了公眾對城市暴雨洪澇應急管理的熱議,也凸顯了解析降雨時程特征及準確描繪暴雨洪水過程對于城市洪澇運行管理的重要性.

    國內外有關設計降雨的研究大多采用區(qū)域經驗性暴雨公式結合芝加哥雨型法[1-2]、Pilgrim&Cordery 雨型法[3],但它們僅適用于城市洪澇模擬或風險評估研究[4-6],得出的結果與實際降雨過程的差異性較大. 為改進套用固定雨型計算的不足,有學者嘗試用其他方法進行設計降雨. 張建平等[7]引入量化指標暴雨衰減指數(shù)對雨型進行設計,結果表明采用暴雨衰減指數(shù)設計的雨型具有區(qū)域代表性;唐明等[8]分析了不同典型暴雨放大方法的優(yōu)缺點,最終得出“復合雨型”可以兼顧城市排水與流域排澇系統(tǒng)的能力復核,雨型生成過程相對簡單、便于應用的結論;朱秀迪等[9]采用Circular統(tǒng)計法分析了北京市不同類型降水的精細化特征,對于改進區(qū)域數(shù)值模式和提高氣象預測水平有較為重要的現(xiàn)實意義;程昶等[10]采用牛頓-高斯法計算了各重現(xiàn)期下的暴雨強度分公式和暴雨強度總公式的參數(shù),并通過誤差檢驗發(fā)現(xiàn)暴雨強度分公式的精度符合要求,適用于計算相應重現(xiàn)期下各歷時的暴雨強度. 盡管國內外學者在設計暴雨雨型、推求暴雨強度公式上做了大量研究,但這些研究在時間尺度上都集中于長短歷時暴雨的設計,鮮有研究對汛期和非汛期的典型降雨過程進行設計和計算.

    本研究選取河北省永年站1980—2018年4月1日至10月31日的逐小時降雨數(shù)據(jù)為研究數(shù)據(jù),在分析其數(shù)百場降雨年際和年內分布特征的基礎上,采用以降雨量和降雨歷時為主的排頻法分別對永年站汛期、非汛期的典型降雨過程進行設計和分析,以期為城市洪澇計算、灌區(qū)干旱應急管理提供關鍵支撐. 本研究也可為其他區(qū)域典型設計降雨的計算提供參考.

    1 氣象站概況與數(shù)據(jù)處理

    1.1 氣象站概況

    本研究選擇的典型氣象站——永年站是國家級氣象站,位于河北省邯鄲市永年區(qū)境內,永年區(qū)地理坐標為東經114°20′~114°52′,北緯36°35′~36°56′,如圖1所示. 永年站位于半濕潤半干旱地區(qū),屬暖溫帶大陸性季風氣候,年平均氣溫14 ℃,最冷月份1月的平均氣溫-2.5 ℃,極端最低氣溫-20 ℃,最熱月份7月的平均氣溫27 ℃,極端最高氣溫42.5 ℃,全年無霜期200 d,年日照2557 h. 該氣象站冬季寒冷干燥、春季風多雨少、秋季天高氣爽、夏季炎熱多雨的氣象特征與華北地區(qū)的氣象特征較為一致[11].

    圖1 研究區(qū)位置圖Fig.1 Study area location

    1.2 數(shù)據(jù)來源與數(shù)據(jù)處理

    本研究只考慮降雨,不考慮雪、霜等其他降水. 結合區(qū)域調研的結果,選取永年站1980—2018年4月1日至10月31日的逐小時降雨數(shù)據(jù)為研究數(shù)據(jù). 所有數(shù)據(jù)來源于“中國國家級地面氣象站逐小時降水數(shù)據(jù)集(V1.0)”,所有數(shù)據(jù)都經過質量控制,且每組數(shù)據(jù)都設有質量控制碼. 其中,每組降雨數(shù)據(jù)由20—21 時開始記錄,至次日19—20 時結束,單位為0.1 mm. 當涉及連續(xù)降雨時,將發(fā)生降雨的前后6 h 未檢測到雨情的降雨過程劃分為一場降雨事件[12-13]. 由于數(shù)據(jù)較為龐大,以區(qū)域是否產流作為初步篩選條件,依據(jù)華北平原產流特性[14],將降雨量小于10 mm 的降雨場次剔除,僅將降雨總量≥10 mm 的392 場降雨過程用于典型降雨過程的設計.

    2 降雨分布特征分析

    2.1 計算方法

    以篩選后的永年站長序列逐小時降雨觀測數(shù)據(jù)為基礎,統(tǒng)計該氣象站各年總降雨量,并采用算數(shù)平均法對各年總降雨量均值進行計算. 降雨年際特征分布分析采用5 a滑動平均法[15],其表達式如式(1)所示.

    式中:Yt為t點滑動平均值. 分析降雨年內分布特征時,對每年4—10月各月分配的降雨量求平均值處理.

    2.2 降雨年際分布特征分析

    通過算數(shù)平均法可計算得出永年站1980—2018年的多年平均降雨量為405.2 mm,年際降雨量變化過程如圖2 所示. 由圖2 可知,永年站各年總降雨量在多年平均降雨量附近上下波動,且波動范圍較大,其中2016年的總降雨量最大,為765.7 mm,1997年的總降雨量最小,為183.6 mm,極值比為4.17,年降雨量最大值和最小值與多年平均降雨量之比分別為1.89、0.45,說明永年站年際降雨量變化較大.

    圖2 1980—2018年永年站年際降雨量變化Fig.2 Inter-annual rainfall variation of Yongnian Station from 1980 to 2018

    分析1980—2018 年的5 a 滑動平均值可知,降雨量最大值為501.52 mm,發(fā)生在2012—2016 年;降雨量最小值為301.98 mm,發(fā)生在1985—1989年,兩者相差199.54 mm. 對比多年5 a滑動平均值可知:2006年以前永年站降雨量相對偏少,其中1980—1982年、1991年、1999—2000年的降雨量略多,2006年后永年站的降雨量整體偏多,僅在2010—2011年略有下降. 總體上看,1980—2018年永年站的降雨量在多年平均降雨量附近上下波動,具有明顯的豐平枯交替變化特征.

    2.3 降雨年內分布特征分析

    依據(jù)2.1小節(jié)所述方法分析永年站年內降雨量變化,結果如圖3所示. 整體上看,永年站年內降雨量分布相對集中,降雨量最多的三個月集中于6—8月,分別占統(tǒng)計降雨總量的11.4%、31.6%、26.9%,且7月份所占比例最大;9月份降雨量也相對較多,為40.4 mm,約占統(tǒng)計降雨總量的10.0%;降雨量最少的三個月份分別為4 月、5 月、10 月,這三個月份的降雨量分別占統(tǒng)計降雨總量的5.4%、9.3%、5.6%. 其中,6—9月的總降雨量達到323.2 mm,占統(tǒng)計降雨總量的79.8%,4月、5月和10 月的總降雨量為82 mm,占統(tǒng)計降雨總量的20.2%.

    圖3 1980—2018年永年站年內降雨量變化Fig.3 Intra-annual rainfall variation of Yongnian Station from 1980 to 2018

    此外,分別計算了永年站1980—2018 年中每年6—9 月的降雨總量及其占該年降雨總量的比例. 結果顯示,1980—2018 年永年站每年6—9 月的降雨總量占該年降雨總量的比例在57%~97%之間. 其中,1981年6—9月的降雨總量占該年降雨總量的比例最大,達96.2%;1991年6—9月的降雨總量占該年降雨總量的比例最小,為57.9%. 永年站的年內降雨分布特征和其他學者[16-17]對華北地區(qū)年內降雨分布特征的分析結論一致.

    2.4 汛期和非汛期的降雨統(tǒng)計分析

    分別對永年站1980—2018年的汛期及非汛期降雨量數(shù)據(jù)進行統(tǒng)計,結果如圖4所示. 結果顯示,39年間永年站汛期降雨量的最大值是660.2 mm,發(fā)生在2016年,最小值是134.8 mm,發(fā)生在1986年;非汛期降雨量的最大值是183.1 mm,發(fā)生在1983年,最小值是22.7 mm,發(fā)生在1997年. 2016年永年站汛期和非汛期的降雨量差值最大,為554.7 mm.

    圖4 1980—2018年永年站汛期及非汛期各年降雨量Fig.4 Annual rainfall in flood season and non-flood season of Yongnian Station from 1980 to 2018

    盡管永年站1983年的非汛期降雨量最多,但因為其汛期降雨量僅為288 mm,所以會導致其年內降雨量偏少,年際降雨量也偏少. 雖然永年站2016年的非汛期降雨量較少,僅為105.5 mm,但其汛期降雨量卻高達660.2 mm,最終使得其年內降雨量及年際降雨量偏多. 以上結果進一步說明,汛期降雨量可顯著影響區(qū)域年內、年際降雨分布特征.

    綜上,考慮到永年站汛期與非汛期降雨量的差異性以及汛期降雨量對年內、年際降雨分布特征影響的顯著性,本研究將降雨數(shù)據(jù)按照汛期、非汛期進行劃分,并分別推求不同時期的典型降雨過程,以保證最大限度地反映永年站所在區(qū)域的降雨特征,為產流規(guī)律研究提供準確、可靠的設計降雨過程. 經統(tǒng)計最終得出,1980—2018年永年站汛期共308場降雨,非汛期共84場降雨.

    3 典型設計降雨計算與分析

    3.1 計算方法

    降雨量、降雨歷時、降雨峰值、峰值比是用于表述降雨特征的重要變量[18]. 降雨量是最能直觀體現(xiàn)降雨總體特征的最常用的指標. 降雨歷時是指從降雨開始時刻至降雨結束時刻的時長,在此以h為單位,其長短也可直觀地反映一場降雨的總體特征. 降雨峰值是指一場降雨中某一時刻降雨量的最大值,單位mm,它是反映降雨極值特征的一個重要指標. 峰值比可用來反映降雨峰值出現(xiàn)的時間及位置,是較好反映降雨過程的一個定性指標,其計算公式如下:

    式中:ti表示降雨峰值出現(xiàn)的時間;T表示降雨總歷時;R表示峰值比,該峰值比可以具體地反映出降雨峰值出現(xiàn)的時間.

    在參考文獻[19-20]的基礎上,結合永年站降雨分布特征,以峰值比作為判定要素,可將典型降雨過程分為6 種,如圖5 所示. 圖5(a)、(b)、(c)展示的三種雨型均為單峰型雨型,其降雨峰值分別位于前部、中部、后部,圖5(d)、(e)、(f)展示的三種雨型均為雙峰型雨型,其降雨峰值特點分別為前大后小、前后相等、前小后大.本研究首先根據(jù)降雨量和降雨歷時這兩個能直觀反映降雨總體特征的指標對永年站1980—2018年的降雨場次進行分級和排頻計算,然后依據(jù)降雨峰值、峰值比這兩個指標再次對篩選出的降雨場次進行分級和排頻計算,最后通過對比分析求得永年站各時期典型降雨過程的特征值,具體過程如下.

    圖5 不同峰值特征的典型降雨過程示意圖Fig.5 Typical rainfall processes with different peak characteristics

    3.1.1 以降雨量等級為主的排頻計算

    首先依據(jù)降水量等級劃分標準[21],結合永年站實測大暴雨數(shù)據(jù),將1980—2018 年永年站所有降雨場次的降雨量分為6 個等級,如表1 所示. 然后以不同等級的降雨量為變量進行排頻,將發(fā)生頻次最高的降雨量等級對應的所有降雨場次篩選出來,并以篩選后的這些降雨數(shù)據(jù)作為計算樣本(記為樣本1)來進行典型降雨過程的設計. 最后計算該等級降雨量的平均值,并將其作為典型降雨過程的降雨總量.

    表1 降雨量等級劃分表Tab.1 Classification table of rainfall grades

    將樣本1 中的降雨場次分別以降雨峰值和降雨歷時為變量進行排頻,其中降雨峰值每1 mm為一個等級,降雨歷時每1 h為一個等級,對樣本1中發(fā)生頻次最高的降雨峰值及降雨歷時等級對應的降雨場次篩選出來,并對篩選出來的所有降雨場次的降雨峰值和降雨歷時分別求平均值,即為典型降雨過程的降雨峰值和降雨歷時.

    從樣本1中將發(fā)生頻次最高的降雨歷時等級對應的降雨場次篩選出來,記為樣本2. 參考公式(2),再對樣本2的降雨場次以峰值比為變量進行排頻,將發(fā)生頻次最高的峰值比對應的降雨場次篩選出來,并對篩選出來的降雨場次的峰值比求平均值,即為典型降雨過程的峰值比.

    通過以上方法可得到一套典型降雨過程的降雨量、降雨峰值、降雨歷時、峰值比這四種基本變量的初值. 為方便理解,在此記為A過程特征值.

    3.1.2 以降雨歷時等級為主的排頻計算

    首先以所有降雨場次的降雨歷時為變量進行等級劃分,并進行排頻,每1 h為一個等級,將發(fā)生頻次最高的降雨歷時等級對應的所有降雨場次篩選出來,用于下一步典型降雨過程的設計分析計算,在此記為樣本3. 然后計算該等級降雨歷時的平均值,即為典型降雨過程的降雨歷時.

    將樣本3中的所有降雨場次分別以降雨量和峰值比為變量進行排頻,對發(fā)生頻次最高的降雨量及峰值比等級對應的降雨場次篩選出來,再對篩選出來的所有降雨場次的降雨量和峰值比分別求平均值,即為典型降雨過程的降雨量及峰值比.

    從樣本3中將發(fā)生頻次最高的降雨量等級所對應的降雨場次篩選出來,記為樣本4,再將這些降雨場次(樣本4)以降雨峰值為變量進行排頻,對發(fā)生頻次最高的降雨峰值等級對應的降雨場次篩選出來,并對這些篩選出來的降雨場次的降雨峰值求平均值,即為典型降雨過程的降雨峰值.

    通過以上方法可得到一套典型降雨過程的降雨量、降雨峰值、降雨歷時、峰值比這四種基本變量的初值. 為方便理解,在此記為B過程特征值.

    3.1.3 典型降雨過程推求方法

    對以上兩種特征排頻法分別推求出的典型降雨過程的特征值進行比較分析:

    1)A、B兩種過程所得出的特征值基本一致. 若A、B兩種過程得出的降雨量和降雨峰值的等級一致,且A、B兩種過程的均值相差不大,則對兩種過程的結果再次求均值,即為設計典型降雨過程的降雨量及降雨峰值;若A、B兩種過程得出的峰值比及降雨歷時等級均一致,則B過程得出的峰值比和降雨歷時即為設計典型降雨過程的峰值比和降雨歷時.

    2)A、B兩種過程所得出的特征值有差異. 若A、B兩種過程得出的降雨量等級不一致,不管A、B兩種過程的均值相差多少,都以A過程求得的降雨量作為設計典型降雨過程的降雨量;若A、B兩種過程得出的降雨峰值等級不一致,則以不同的降雨峰值等級為變量對所有降雨場次進行排頻求均值,選擇與其最接近的特征值初值作為設計典型降雨過程的降雨峰值;若A、B兩種過程得出的降雨歷時等級不一致,不管A、B兩種過程的結果相差多少,均以B過程求得的降雨歷時作為設計典型降雨過程的降雨歷時;由于A、B兩種過程得出的峰值比都是在降雨歷時的基礎上進行排頻的,因此當A、B兩種過程得出的峰值比不一致時,同樣選擇B過程求得的峰值比作為設計典型降雨過程的峰值比,即對樣本3進行排頻計算后得出的峰值比.

    依據(jù)上述推求方法,即可推求出典型降雨過程的降雨量、降雨歷時、降雨峰值、峰值比這四個變量最終的特征值. 然后,同樣以出現(xiàn)頻次為依據(jù),以峰值比的峰型及位置為節(jié)點,以降雨峰值為控制點,將不同典型降雨過程(樣本3)前后的增加和減少過程進行逐小時繪制,分別選取不同時段中出現(xiàn)頻次最高的雨強(mm/h),即可得出不同時段雨強之間的比值,之后再以降雨量和降雨歷時為控制點,進行不同時段的同倍比放大,即可設計出典型降雨過程. 典型降雨過程的具體推求過程如圖6所示.

    圖6 典型降雨過程的推求過程Fig.6 The designed workflow of typical rainfall process

    3.2 汛期典型降雨統(tǒng)計分析

    先依據(jù)3.1.1 小節(jié)方法,以降雨量為變量對永年站1980—2018 年汛期的308 場降雨進行排頻計算. 經統(tǒng)計,1980—2018年永年站汛期發(fā)生頻次最高的降雨量等級為10.0~24.9 mm的中雨,其對應的降雨場次共179 場,占汛期總降雨場次的58.1%,其降雨量均值為16.2 mm. 這179場降雨中,降雨峰值為6~7 mm等級的降雨發(fā)生的頻次最多,達24場,占所選降雨量等級場次的13.4%,且這24場降雨的降雨峰值均值為6.2 mm;降雨歷時等級為2 h的降雨發(fā)生的頻次最多,達26場,占所選降雨量等級場次的14.5%;峰值比的峰值特征為單峰靠前的降雨發(fā)生的頻次最多,達18場,占所選降雨歷時等級場次的69.2%.

    再根據(jù)3.1.2小節(jié)方法,以降雨歷時為變量對永年站1980—2018年汛期的308場降雨進行排頻計算. 經統(tǒng)計,1980—2018年永年站汛期發(fā)生頻次最高的降雨歷時等級為3 h,其對應的降雨場次的共40場,占汛期總降雨場次的13.0%. 這40場降雨中,發(fā)生頻次最高的降雨量等級為10.0~24.9 mm的中雨,其對應的降雨場次共23場,占所選降雨歷時等級場次的57.5%,且這23場降雨的降雨量均值為16.7 mm;峰值比的峰值特征為單峰靠后的降雨發(fā)生的頻次最多,共20場,占所選降雨歷時等級場次的50%;降雨峰值為11~12 mm等級的降雨發(fā)生的頻次最多,共11場,占所選降雨量等級場次的47.8%,且這11場降雨的降雨峰值均值為11.6 mm.

    在此基礎上,按照3.1.3小節(jié)的方法即可推求出永年站汛期典型降雨過程最終的特征值,如表2所示. 由表2可知,兩種方法計算出的發(fā)生頻次最高的降雨量等級均為10.0~24.9 mm的中雨,且得出的降雨量均值相差不大,降雨峰值、降雨歷時和峰值比的結果都存在較大差異. 就降雨量而言,對兩種方法的結果求平均值,即可得到永年站汛期典型降雨過程的降雨量,為16.5 mm;就降雨峰值而言,對永年站汛期所有降雨場次的降雨峰值進行排頻,發(fā)現(xiàn)降雨峰值等級為5~6 mm 的降雨發(fā)生的頻次最高,共24 場,占汛期總降雨場次的7.8%,且這24 場降雨的降雨峰值均值為5.8 mm,再對比兩種方法得出的降雨峰值均值,最后確定永年站汛期典型降雨過程的降雨峰值為6.2 mm;因峰值比和降雨歷時出現(xiàn)差異時,應以降雨歷時為變量排頻得出的結果為主,故最終確定永年站汛期典型降雨過程的降雨歷時為3 h,峰值比的峰值特征為單峰靠后.

    表2 基于不同排頻方法得出的永年站汛期典型降雨過程特征值Tab.2 Characteristic values of typical rainfall process in flood season of Yongnian Station obtained by different frequency analysis methods

    依據(jù)表2結果,結合3.1.3小節(jié)所述方法可得出永年站汛期典型降雨過程(圖7). 由圖7 可知,永年站汛期典型降雨過程的雨型屬于單峰型,降雨量逐漸增多,降雨歷時短,降雨峰值位于降雨過程的后部,降雨峰值時段降雨量占汛期降雨過程總降雨量的37.6%,降雨峰值時段的雨強大,降雨量在整個降雨過程后期更為集中. 根據(jù)實地調研并參考其他文獻[17]可知,華北地區(qū)汛期降雨多為對流雨,對流雨具有強度大、雨量多、歷時短的特點,通常表現(xiàn)為雷暴雨、短時強降雨、陣性降雨等. 由此可知,本研究結果符合區(qū)域汛期降雨特征.

    圖7 永年站汛期典型降雨過程Fig.7 Typical rainfall process in flood season of Yongnian Station

    3.3 非汛期典型降雨統(tǒng)計分析

    先根據(jù)3.1.1小節(jié)計算方法,以降雨量為變量對永年站1980—2018年非汛期的84場降雨進行排頻計算.經統(tǒng)計,1980—2018年永年站非汛期發(fā)生頻次最高的降雨量等級為10.0~24.9 mm 的中雨,其對應的降雨場次共59 場,占非汛期總降雨場次的70.2%,其降雨量均值為15.8 mm. 這59 場降雨中,降雨峰值為2~3 mm等級的降雨發(fā)生的頻次最多,達10 場,占所選降雨量等級場次的16.9%,且這10 場降雨的降雨峰值均值為2.6 mm;降雨歷時等級為10 h的降雨發(fā)生的頻次最多,達7場,占所選降雨量等級場次的11.9%;峰值比的峰值特征為單峰居中的降雨發(fā)生的頻次最多,達4場,占所選降雨歷時等級場次的57.1%.

    再根據(jù)3.1.2 小節(jié)方法,以降雨歷時為變量對永年站1980—2018 年非汛期的84 場降雨進行排頻計算.經統(tǒng)計,1980—2018年永年站非汛期發(fā)生頻次最高的降雨歷時等級為10 h,其對應的降雨場次共9場,占非汛期總降雨場次的10.7%. 這9場降雨中,發(fā)生頻次最高的降雨量等級為10.0~24.9 mm的中雨,其對應的降雨場次共7場,占所選降雨歷時等級場次的77.8%,且這7場降雨的降雨量均值為15.9 mm;峰值比的峰值特征為單峰居中的降雨發(fā)生的頻次最多,共3場,占所選降雨歷時等級場次的33.3%;降雨峰值為3~4 mm等級的降雨發(fā)生的頻次最多,共2場,占所選降雨量等級場次的28.6%,且這2場降雨的降雨峰值均值為3.2 mm.

    在此基礎上,按照3.1.3 小節(jié)方法即可推求出永年站非汛期典型降雨過程最終的特征值,如表3 所示.由表3可知,兩種方法計算出的發(fā)生頻次最高的降雨量等級均為10.0~24.9 mm的中雨,且得出的降雨量均值相差不大,故對兩種方法得出的降雨量再次求平均值即可得到永年站非汛期典型降雨過程的降雨量,為15.9 mm. 因為兩種方法得出的降雨歷時和峰值比完全一致,所以永年站非汛期典型降雨過程的降雨歷時為10 h,峰值比的峰值特征為單峰居中. 由于兩種方法計算得出的降雨峰值等級不一致,因此需對永年站非汛期所有降雨場次的降雨峰值進行排頻,結果顯示降雨峰值等級為2~3 mm的降雨發(fā)生的頻次最多,共15場,占非汛期總降雨場次的17.9%,且這15場降雨的降雨峰值均值為2.8 mm,與兩種方法得出的降雨峰值均值進行對比后,最終確定的永年站非汛期典型降雨過程的降雨峰值為2.6 mm.

    表3 基于不同排頻方法得出的永年站非汛期典型降雨過程特征值Tab.3 Characteristic values of typical rainfall process in non-flood season of Yongnian Station obtained by different frequency analysis methods

    同樣,可得永年站非汛期典型降雨過程(圖8). 由圖8可知,永年站非汛期典型降雨過程的雨型也屬于單峰型,降雨峰值位于中部,降雨峰值時段的降雨量僅占非汛期降雨過程總降雨量的16.4%,前期持續(xù)降雨且降雨量呈階梯狀逐步增加,后期降雨量呈緩慢下降趨勢,且降雨峰值前的降雨量要低于降雨峰值后的降雨量. 根據(jù)調研可知,華北地區(qū)非汛期降雨大多以鋒面雨為主,具有持續(xù)時間長、雨強小等降雨特征. 由以上分析可知,本研究結果符合區(qū)域非汛期降雨特征.

    圖8 永年站非汛期典型降雨過程Fig.8 Typical rainfall process in non-flood season of Yongnian Station

    綜合對比永年站汛期和非汛期的典型降雨過程可知,汛期典型降雨過程的雨強較大、降雨過程短而急,非汛期典型降雨過程的雨強相對較小,降雨過程長而緩;就雨型而言,汛期及非汛期典型降雨過程的降雨量均集中在中后部.

    華北地區(qū)屬于典型的暖溫帶半濕潤大陸性季風氣候,其汛期降雨量的變化是影響其年降雨量變化的主要因素. 由于華北地區(qū)汛期降雨具有時程短、變化急、雨量大的特點,因此該類型降雨的發(fā)生往往具有突發(fā)性. 若此類型降雨發(fā)生在西部山丘地區(qū),則極易產生山洪且極易伴有滑坡、泥石流等次生災害發(fā)生;若此類型降雨發(fā)生在城市區(qū)域,則極易形成城中“看?!爆F(xiàn)象,且會嚴重影響居民交通出行,甚至造成人員傷亡. 因此,應針對區(qū)域汛期的降雨特征開展典型降雨過程設計并用于相應的工程防洪排澇計算,以進一步減少區(qū)域財產損失. 由于非汛期降雨的降雨時程長且往往發(fā)生在4—5月冬小麥關鍵生長期,因此其是影響區(qū)域糧食產量的重要因素,故應針對區(qū)域非汛期的降雨特征開展典型降雨過程設計并用于相應的灌區(qū)設計及其灌溉制度的科學制定,以優(yōu)化農業(yè)水資源的可持續(xù)利用,進一步增加灌區(qū)農作物產量及其用水效益和效率.

    此外,針對不同時期的降雨特征開展人工降雨設計及其產流試驗和規(guī)律分析,可為深入了解區(qū)域產流機理及水資源特性等研究提供關鍵支撐.

    4 結論

    基于永年站1980—2018年4月1日至10月31日的逐小時降雨數(shù)據(jù),對其年際和年內降雨分布特征進行了分析,然后以降雨量、降雨歷時、降雨峰值和峰值比為關鍵要素,通過排頻計算分別推求出永年站汛期及非汛期的典型降雨過程,得出以下結論:

    1)永年站近39 年的降雨過程具有明顯的豐平枯交替變化特征,且大部分屬于枯水年,只有1980—1982年、1991年、1999—2000年以及2007—2018年屬于豐水年;每年6—9月的降雨量占比較大,在57%~97%之間. 汛期降雨量對區(qū)域降雨年內、年際分布特征影響顯著.

    2)永年站汛期典型降雨過程具有短而急的特點,降雨量主要集中在后部,設計典型降雨過程的降雨量為16.5 mm、降雨峰值為6.2 mm、降雨歷時為3 h、雨型為單峰靠后;非汛期典型降雨過程具有長而緩的特點,降雨量主要集中在中后部,設計典型降雨過程的降雨量為15.9 mm、降雨峰值為2.6 mm、降雨歷時為10 h、雨型為單峰居中.

    猜你喜歡
    雨型場次歷時
    長江上游高洪水期泥沙輸移特性
    量詞“只”的形成及其歷時演變
    常用詞“怠”“惰”“懶”的歷時演變
    概化的累計暴雨量百分數(shù)法在太湖區(qū)域設計暴雨雨型研究的應用
    江蘇水利(2020年9期)2020-10-09 02:53:52
    天津市設計暴雨雨型的演變
    水資源保護(2020年1期)2020-02-05 02:52:26
    基于運行場次用時誤差的載人設備故障預警可視化研究
    排考場次分配方法及其SQL實現(xiàn)
    對《紅樓夢》中“不好死了”與“……好的”的歷時考察
    紅樓夢學刊(2019年5期)2019-04-13 00:42:36
    深圳市流域暴雨雨型及變化趨勢分析
    古今字“兌”“說”“悅”“?!睔v時考察
    国产精品一区二区三区四区久久| 国产精品女同一区二区软件 | 狠狠狠狠99中文字幕| 男人舔奶头视频| 午夜福利在线在线| 日本与韩国留学比较| 欧美在线黄色| 美女黄网站色视频| 亚洲熟妇中文字幕五十中出| 丰满的人妻完整版| 丰满人妻一区二区三区视频av | 久久久成人免费电影| 88av欧美| 色播亚洲综合网| 中亚洲国语对白在线视频| 亚洲狠狠婷婷综合久久图片| 午夜亚洲福利在线播放| 桃红色精品国产亚洲av| x7x7x7水蜜桃| 国产精品久久电影中文字幕| 操出白浆在线播放| 少妇人妻精品综合一区二区 | 日韩欧美 国产精品| 99国产极品粉嫩在线观看| 男人舔女人下体高潮全视频| 伊人久久大香线蕉亚洲五| 精品一区二区三区视频在线观看免费| 久久精品人妻少妇| 99久久精品热视频| 亚洲av美国av| 国产一区二区亚洲精品在线观看| 日本a在线网址| 男插女下体视频免费在线播放| 亚洲人成电影免费在线| 中文字幕精品亚洲无线码一区| 三级毛片av免费| 亚洲欧美激情综合另类| 国产单亲对白刺激| 99热这里只有是精品50| 91九色精品人成在线观看| 久久国产精品影院| 亚洲专区中文字幕在线| 日韩免费av在线播放| 欧美av亚洲av综合av国产av| www.999成人在线观看| 亚洲性夜色夜夜综合| 两个人视频免费观看高清| 国产成人影院久久av| 99久久精品国产亚洲精品| 亚洲久久久久久中文字幕| 在线播放国产精品三级| 欧美日韩亚洲国产一区二区在线观看| 日本撒尿小便嘘嘘汇集6| 国产成人欧美在线观看| 精品午夜福利视频在线观看一区| 亚洲第一欧美日韩一区二区三区| 狂野欧美激情性xxxx| 亚洲av不卡在线观看| 精品久久久久久久久久久久久| 黄色丝袜av网址大全| 久久久久久久久中文| 国产毛片a区久久久久| 我的老师免费观看完整版| 好男人电影高清在线观看| 国产男靠女视频免费网站| 亚洲国产精品sss在线观看| 国产三级黄色录像| 一个人观看的视频www高清免费观看| 免费看十八禁软件| 少妇人妻精品综合一区二区 | 欧美在线黄色| 国产亚洲欧美在线一区二区| 欧美黑人巨大hd| 久久天躁狠狠躁夜夜2o2o| 精品欧美国产一区二区三| 久久精品国产99精品国产亚洲性色| 亚洲国产精品久久男人天堂| 看黄色毛片网站| 黄色日韩在线| 伊人久久大香线蕉亚洲五| 国产高清激情床上av| 九九热线精品视视频播放| a在线观看视频网站| x7x7x7水蜜桃| 午夜久久久久精精品| 99国产精品一区二区三区| 免费在线观看影片大全网站| 一级a爱片免费观看的视频| 国产一区二区在线观看日韩 | 一本久久中文字幕| 久久午夜亚洲精品久久| 91麻豆av在线| 成人一区二区视频在线观看| 在线观看舔阴道视频| 黄色日韩在线| 国产精品综合久久久久久久免费| 国产蜜桃级精品一区二区三区| 亚洲国产欧美人成| 欧美一区二区亚洲| 18禁裸乳无遮挡免费网站照片| 一卡2卡三卡四卡精品乱码亚洲| 在线观看舔阴道视频| 91久久精品电影网| av欧美777| 亚洲熟妇熟女久久| 两个人视频免费观看高清| 久久性视频一级片| 午夜精品久久久久久毛片777| 亚洲人与动物交配视频| 亚洲av不卡在线观看| 国产色婷婷99| 麻豆一二三区av精品| 内射极品少妇av片p| 国产乱人视频| 国产精品日韩av在线免费观看| 国产精华一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品一区二区www| ponron亚洲| 18禁国产床啪视频网站| 我要搜黄色片| 国产麻豆成人av免费视频| 级片在线观看| 给我免费播放毛片高清在线观看| 女警被强在线播放| 内地一区二区视频在线| 窝窝影院91人妻| 99视频精品全部免费 在线| 国产爱豆传媒在线观看| 亚洲精品成人久久久久久| 中文字幕人妻丝袜一区二区| 日本免费a在线| 国产淫片久久久久久久久 | 久久久久久久久大av| 亚洲精品粉嫩美女一区| 欧美性感艳星| www.www免费av| 一级黄片播放器| 日本撒尿小便嘘嘘汇集6| 51国产日韩欧美| 亚洲人与动物交配视频| 国产亚洲欧美98| 免费观看精品视频网站| 亚洲最大成人中文| 999久久久精品免费观看国产| 一a级毛片在线观看| 国产一区二区在线av高清观看| 香蕉丝袜av| 亚洲人成电影免费在线| 国产aⅴ精品一区二区三区波| 少妇高潮的动态图| 欧美在线黄色| 欧美性猛交黑人性爽| e午夜精品久久久久久久| 婷婷亚洲欧美| 中文字幕人妻丝袜一区二区| 特大巨黑吊av在线直播| 在线观看av片永久免费下载| 一个人免费在线观看电影| 高清在线国产一区| 亚洲专区中文字幕在线| 18美女黄网站色大片免费观看| 国内精品久久久久精免费| av天堂在线播放| 色播亚洲综合网| 高潮久久久久久久久久久不卡| 欧美又色又爽又黄视频| 亚洲精华国产精华精| 国产成人av激情在线播放| 久久久久亚洲av毛片大全| 久久久久久久精品吃奶| 午夜精品一区二区三区免费看| 女人高潮潮喷娇喘18禁视频| 桃红色精品国产亚洲av| 亚洲最大成人中文| 亚洲国产欧美人成| 精品午夜福利视频在线观看一区| 亚洲 欧美 日韩 在线 免费| 欧美成人性av电影在线观看| 午夜精品一区二区三区免费看| 2021天堂中文幕一二区在线观| 国产精品 国内视频| 看免费av毛片| 亚洲自拍偷在线| a在线观看视频网站| 在线国产一区二区在线| 亚洲,欧美精品.| 欧美色欧美亚洲另类二区| h日本视频在线播放| 九色成人免费人妻av| 免费在线观看亚洲国产| 美女被艹到高潮喷水动态| 不卡一级毛片| 精品国产美女av久久久久小说| 一区二区三区免费毛片| 看片在线看免费视频| 18禁黄网站禁片免费观看直播| 黄片大片在线免费观看| 欧美精品啪啪一区二区三区| 有码 亚洲区| 亚洲av不卡在线观看| 国产一区二区三区视频了| 制服丝袜大香蕉在线| 最新美女视频免费是黄的| 精品人妻1区二区| 嫩草影院入口| 欧美+日韩+精品| 国产精品久久久久久精品电影| 最近最新中文字幕大全免费视频| 亚洲avbb在线观看| 色哟哟哟哟哟哟| 欧美日韩乱码在线| 国产精品一区二区免费欧美| 亚洲国产精品sss在线观看| 一边摸一边抽搐一进一小说| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产单亲对白刺激| 99久久久亚洲精品蜜臀av| 最近最新中文字幕大全免费视频| 脱女人内裤的视频| 日韩欧美一区二区三区在线观看| 亚洲一区高清亚洲精品| 99久久久亚洲精品蜜臀av| 女人十人毛片免费观看3o分钟| 欧美一级a爱片免费观看看| 亚洲精品成人久久久久久| 久久久国产成人免费| 国产精品国产高清国产av| 天堂√8在线中文| 人妻久久中文字幕网| ponron亚洲| 久久天躁狠狠躁夜夜2o2o| 久久精品国产综合久久久| 欧美午夜高清在线| 欧美极品一区二区三区四区| 婷婷精品国产亚洲av在线| 高清毛片免费观看视频网站| 日本熟妇午夜| 久久精品影院6| 色av中文字幕| 最好的美女福利视频网| 亚洲五月婷婷丁香| 一个人看视频在线观看www免费 | netflix在线观看网站| 黄片小视频在线播放| 国产97色在线日韩免费| 国产午夜福利久久久久久| 搡老妇女老女人老熟妇| 午夜免费男女啪啪视频观看 | 亚洲成av人片在线播放无| 国产伦精品一区二区三区视频9 | 久久久久久九九精品二区国产| 丝袜美腿在线中文| 特级一级黄色大片| 亚洲熟妇熟女久久| 国产aⅴ精品一区二区三区波| 无遮挡黄片免费观看| 久久久久九九精品影院| 精品福利观看| 一二三四社区在线视频社区8| 国产精品久久电影中文字幕| 精品一区二区三区视频在线 | 91久久精品国产一区二区成人 | 在线免费观看的www视频| 青草久久国产| 日本成人三级电影网站| 夜夜爽天天搞| 少妇熟女aⅴ在线视频| 久久久久久久午夜电影| 夜夜躁狠狠躁天天躁| 在线观看一区二区三区| 国产成人欧美在线观看| 男女午夜视频在线观看| 欧美成人a在线观看| 国产精品 欧美亚洲| 国产午夜福利久久久久久| 久久久久久久精品吃奶| 久久精品亚洲精品国产色婷小说| 一个人免费在线观看的高清视频| 国产精品电影一区二区三区| 搡女人真爽免费视频火全软件 | 老汉色∧v一级毛片| 色播亚洲综合网| 人人妻人人看人人澡| 亚洲精品在线观看二区| av国产免费在线观看| 亚洲无线观看免费| 在线观看av片永久免费下载| 国产精品 国内视频| 成人一区二区视频在线观看| 人妻夜夜爽99麻豆av| 偷拍熟女少妇极品色| 亚洲男人的天堂狠狠| av在线蜜桃| 听说在线观看完整版免费高清| 欧美另类亚洲清纯唯美| 久久久久免费精品人妻一区二区| 在线观看av片永久免费下载| а√天堂www在线а√下载| 亚洲不卡免费看| 日本免费一区二区三区高清不卡| 国产精品女同一区二区软件 | 午夜福利在线观看免费完整高清在 | 在线免费观看的www视频| 麻豆国产av国片精品| 亚洲无线观看免费| 免费电影在线观看免费观看| 午夜精品久久久久久毛片777| 特级一级黄色大片| 99精品欧美一区二区三区四区| 少妇的丰满在线观看| 一级毛片女人18水好多| 九九在线视频观看精品| 高清日韩中文字幕在线| 在线观看日韩欧美| 看片在线看免费视频| 中文亚洲av片在线观看爽| 国产高清视频在线观看网站| 男人和女人高潮做爰伦理| 色尼玛亚洲综合影院| 老汉色∧v一级毛片| 国产精品女同一区二区软件 | 欧美性猛交╳xxx乱大交人| 亚洲国产中文字幕在线视频| 国产一区二区三区视频了| 18禁在线播放成人免费| 黄色视频,在线免费观看| 中文字幕人妻熟人妻熟丝袜美 | 99热这里只有是精品50| 亚洲精品色激情综合| 老汉色∧v一级毛片| 免费看十八禁软件| 久久国产精品影院| 午夜老司机福利剧场| 亚洲色图av天堂| 高清在线国产一区| 免费av毛片视频| 精品久久久久久成人av| 最新中文字幕久久久久| 亚洲avbb在线观看| 国产主播在线观看一区二区| 亚洲国产中文字幕在线视频| 日本 av在线| 99久久综合精品五月天人人| 亚洲成人中文字幕在线播放| 成熟少妇高潮喷水视频| 男人舔奶头视频| 高清毛片免费观看视频网站| 亚洲va日本ⅴa欧美va伊人久久| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美日韩一级在线毛片| 两性午夜刺激爽爽歪歪视频在线观看| 波多野结衣高清作品| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品在线观看二区| 中国美女看黄片| 国产激情偷乱视频一区二区| ponron亚洲| 亚洲欧美一区二区三区黑人| 亚洲第一电影网av| 国产精品99久久99久久久不卡| 高清日韩中文字幕在线| 亚洲成人久久爱视频| 精品福利观看| 18禁在线播放成人免费| 一级毛片女人18水好多| 一区福利在线观看| 狂野欧美白嫩少妇大欣赏| 欧美日韩精品网址| av在线天堂中文字幕| 日本一本二区三区精品| 亚洲精品成人久久久久久| 欧美日韩一级在线毛片| 欧美最新免费一区二区三区 | 韩国av一区二区三区四区| 给我免费播放毛片高清在线观看| 91九色精品人成在线观看| 九色成人免费人妻av| 国产日本99.免费观看| 欧美日韩精品网址| 亚洲av电影在线进入| 最近最新中文字幕大全电影3| 少妇人妻一区二区三区视频| 99在线视频只有这里精品首页| 黄色日韩在线| 搡老妇女老女人老熟妇| 国产熟女xx| 亚洲精品美女久久久久99蜜臀| 校园春色视频在线观看| 97碰自拍视频| av女优亚洲男人天堂| 变态另类成人亚洲欧美熟女| 18禁在线播放成人免费| 欧美一级a爱片免费观看看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲一区高清亚洲精品| 日韩欧美国产一区二区入口| 午夜精品在线福利| 欧美精品啪啪一区二区三区| 亚洲性夜色夜夜综合| 人人妻人人看人人澡| www.999成人在线观看| 国产av在哪里看| 69人妻影院| 国产精品香港三级国产av潘金莲| av在线天堂中文字幕| 最新中文字幕久久久久| 热99在线观看视频| 2021天堂中文幕一二区在线观| 国内揄拍国产精品人妻在线| 国产精品久久久人人做人人爽| 99久久无色码亚洲精品果冻| 午夜视频国产福利| 此物有八面人人有两片| 母亲3免费完整高清在线观看| 亚洲av电影在线进入| 国产视频内射| av天堂中文字幕网| 亚洲狠狠婷婷综合久久图片| 亚洲第一电影网av| 精品久久久久久,| 国产精品嫩草影院av在线观看 | 亚洲avbb在线观看| 欧美成人一区二区免费高清观看| 少妇人妻精品综合一区二区 | 又爽又黄无遮挡网站| 亚洲国产精品合色在线| 久久久久久久久大av| 最近最新免费中文字幕在线| 欧美成人免费av一区二区三区| 亚洲午夜理论影院| 欧美乱妇无乱码| 国产高清videossex| 一个人看的www免费观看视频| 搡老岳熟女国产| 制服丝袜大香蕉在线| 在线看三级毛片| 在线a可以看的网站| 国产一区二区亚洲精品在线观看| 一级作爱视频免费观看| 久久伊人香网站| 国产精品一区二区免费欧美| 中文字幕人成人乱码亚洲影| 少妇人妻一区二区三区视频| 久久九九热精品免费| 精品国内亚洲2022精品成人| av中文乱码字幕在线| 亚洲av电影在线进入| 国产97色在线日韩免费| 亚洲无线观看免费| 国产成人福利小说| 黄色丝袜av网址大全| 日韩欧美精品v在线| 午夜福利在线观看吧| 国产亚洲精品av在线| 欧美午夜高清在线| 亚洲欧美激情综合另类| 免费无遮挡裸体视频| 欧美一区二区国产精品久久精品| 国产精品一区二区免费欧美| 精品人妻1区二区| 在线观看美女被高潮喷水网站 | 日韩欧美三级三区| 日韩欧美在线二视频| 国产午夜精品久久久久久一区二区三区 | 中文字幕高清在线视频| 3wmmmm亚洲av在线观看| 亚洲国产色片| 国语自产精品视频在线第100页| 亚洲熟妇熟女久久| 99久国产av精品| 身体一侧抽搐| 免费人成视频x8x8入口观看| 免费一级毛片在线播放高清视频| 久久亚洲精品不卡| 国产视频内射| 久久久久久九九精品二区国产| 99国产精品一区二区蜜桃av| 国产激情偷乱视频一区二区| 97超视频在线观看视频| 欧美成人免费av一区二区三区| 日本与韩国留学比较| 毛片女人毛片| 激情在线观看视频在线高清| www.999成人在线观看| e午夜精品久久久久久久| 麻豆成人av在线观看| 久久精品影院6| 亚洲欧美精品综合久久99| 我的老师免费观看完整版| 最近最新中文字幕大全免费视频| 国产伦一二天堂av在线观看| 亚洲激情在线av| 日本黄色片子视频| 好男人电影高清在线观看| 给我免费播放毛片高清在线观看| 亚洲精品一卡2卡三卡4卡5卡| 搡老妇女老女人老熟妇| 国产熟女xx| 午夜福利免费观看在线| 国产综合懂色| 美女黄网站色视频| 欧美成人免费av一区二区三区| 国产乱人视频| 91在线精品国自产拍蜜月 | 国内精品美女久久久久久| а√天堂www在线а√下载| 亚洲第一电影网av| ponron亚洲| 欧美在线一区亚洲| 精品一区二区三区av网在线观看| 女人高潮潮喷娇喘18禁视频| 91在线精品国自产拍蜜月 | 最新中文字幕久久久久| 色av中文字幕| 精品国产三级普通话版| 亚洲不卡免费看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品乱码久久久v下载方式 | 亚洲一区高清亚洲精品| 色在线成人网| 一级黄片播放器| 欧美区成人在线视频| 少妇裸体淫交视频免费看高清| 成人18禁在线播放| 免费在线观看影片大全网站| a级毛片a级免费在线| 欧美成人免费av一区二区三区| 国内少妇人妻偷人精品xxx网站| 在线观看美女被高潮喷水网站 | 国产高清激情床上av| 欧美日韩福利视频一区二区| 岛国在线观看网站| 最近最新中文字幕大全免费视频| 日韩国内少妇激情av| 国内毛片毛片毛片毛片毛片| 久久久久久久久中文| 亚洲熟妇熟女久久| 真实男女啪啪啪动态图| 国产私拍福利视频在线观看| 黄色丝袜av网址大全| 丰满的人妻完整版| 最近最新中文字幕大全电影3| 美女黄网站色视频| 国产在视频线在精品| eeuss影院久久| 成年免费大片在线观看| 午夜免费成人在线视频| 两人在一起打扑克的视频| 色在线成人网| 免费看日本二区| 男人舔女人下体高潮全视频| 内地一区二区视频在线| 亚洲片人在线观看| 亚洲男人的天堂狠狠| 久久亚洲精品不卡| 特大巨黑吊av在线直播| 日韩 欧美 亚洲 中文字幕| 99久久精品国产亚洲精品| 精品一区二区三区人妻视频| 天天躁日日操中文字幕| 欧美激情久久久久久爽电影| 亚洲欧美一区二区三区黑人| 国产成+人综合+亚洲专区| a级毛片a级免费在线| 国产激情偷乱视频一区二区| 久久久久久久午夜电影| 一区二区三区国产精品乱码| 亚洲性夜色夜夜综合| 国产免费男女视频| 99国产综合亚洲精品| 成人一区二区视频在线观看| 国产精品一及| 欧美大码av| 久久精品国产亚洲av香蕉五月| 内地一区二区视频在线| www日本黄色视频网| 亚洲一区高清亚洲精品| 婷婷精品国产亚洲av在线| 99精品久久久久人妻精品| 在线观看免费视频日本深夜| 精品免费久久久久久久清纯| 好看av亚洲va欧美ⅴa在| 国产精品爽爽va在线观看网站| 亚洲av中文字字幕乱码综合| 国产三级黄色录像| 啦啦啦观看免费观看视频高清| 欧美性猛交黑人性爽| 亚洲无线观看免费| 美女高潮喷水抽搐中文字幕| 特大巨黑吊av在线直播| 人人妻,人人澡人人爽秒播| 免费观看的影片在线观看| 成人精品一区二区免费| 五月玫瑰六月丁香| 久久久久久久亚洲中文字幕 | 成年女人看的毛片在线观看| 免费看十八禁软件| 亚洲人成网站在线播放欧美日韩| 99riav亚洲国产免费| 久久精品91蜜桃| 亚洲欧美日韩高清专用| ponron亚洲| 国产中年淑女户外野战色| 村上凉子中文字幕在线| 蜜桃亚洲精品一区二区三区| 色噜噜av男人的天堂激情| 老司机福利观看| 久久这里只有精品中国| 国产久久久一区二区三区| 国产欧美日韩一区二区三| 亚洲av成人精品一区久久| 变态另类丝袜制服| 99久久99久久久精品蜜桃| 国产黄色小视频在线观看|