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

    氣候變化條件下降雨-融水型泥石流流量預測方法

    2021-01-16 02:49:09楊宗佶付校龍游勇
    人民長江 2021年12期
    關鍵詞:融水冰川泥石流

    楊宗佶 付校龍 游勇

    摘要:全球氣候變化背景下,目前唯一的冰川泥石流流量經(jīng)驗模型無法定量反映雨熱耦合變化的影響。根據(jù)度日模型冰雪融水當量和泥石流流量雨洪計算方法,提出了考慮雨熱耦合條件下降雨-融水型泥石流流量的計算模型。通過統(tǒng)計51 a降雨和溫度觀測資料的變化趨勢,建立了帕隆藏布流域降雨-融水型泥石流流量計算重現(xiàn)期標準,并通過對典型降雨-融水型泥石流流量的實測資料進行驗證,表明該方法準確性相比經(jīng)驗模型明顯增加。進一步以帕隆藏布流域15條典型降雨-融水型泥石流為例,開展氣候變化條件下泥石流流量的計算和預測。結(jié)果顯示:該方法可量化雨熱耦合條件對泥石流流量的影響,可為氣候變化條件下冰川泥石流的防治關鍵參數(shù)計算和預測提供參考依據(jù)。

    關 鍵 詞:降雨-融水型泥石流; 泥石流流量; 度日模型; 雨熱耦合變化; 氣候變化; 帕隆藏布流域

    中圖法分類號: P642.23

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2021.12.006

    0 引 言

    冰川泥石流廣泛分布于中國青藏高原[1-2]及世界各地的高山地區(qū)[3-5]。隨著全球氣候變暖,冰川消融加劇[6],將導致冰川泥石流頻繁發(fā)生。國內(nèi)外對冰川泥石流展開的相關研究主要是成災條件和機理方面。對于冰川泥石流的形成過程,部分學者認為冰磧物在其中發(fā)揮了重要作用[7-8],還有學者認為冰川泥石流的形成是復合過程[9],并將其啟動過程總結(jié)為侵蝕-滑移型[10]。對于冰川泥石流的形成條件,前人認為其暴發(fā)與水熱組合關系密切,且多發(fā)生在濕熱環(huán)境下[1,11];還有部分學者研究認為其起動和起動方式與氣溫上升下冰川融水[12-13]、黏粒含量[14]以及初始含水率[15]有關。但是目前有關冰川泥石流流量的計算方法的研究較少,只有一套中科院成都山地所提出的經(jīng)驗公式[16],無法反映溫度和降雨兩個關鍵因素對冰川泥石流流量的影響。因此本研究以帕隆藏布流域為研究區(qū)域,以降雨與氣溫[17]為關鍵參數(shù),構(gòu)建了疊加小流域暴雨洪水計算公式和改進度日模型[18-20]的降雨-融水型泥石流流量的計算模型和預測方法,為氣候變化背景下冰川泥石流的工程防治關鍵參數(shù)的選取和預測預報提供了新的方法。

    1 研究方法

    雨熱耦合條件下降雨-融水型泥石流流量的計算方法分別考慮全流域的降雨量和冰雪流域的冰雪融水當量。

    1.1 改進冰雪融水當量計算方法

    全流域24 h內(nèi)冰雪的消融水當量M1根據(jù)冰雪消融與氣溫之間關系建立的度日模型[18-20]進行計算:

    M1=DDF·T24(S1F)(1)

    式中:T24為一天的平均氣溫,℃,考慮冰雪區(qū)域不同海拔位置的氣溫梯度,T24=T-(ELE-2 736)×0.47[21];S1為冰雪融水面積,km2;DDF為冰川或雪的度日因子,mm/(d·℃)。

    改進的度日模型是在劉金平等[22]擬合的公式基礎上進行了改進,考慮了坡度的影響,由于冰雪的消融和其與空氣的接觸面積密切相關[23],不同坡度投影面積下冰雪區(qū)域與空氣的接觸面積不同,因而改進公式如下:

    DDF=(0.009×ELE-0.934×LAT-8.1)×cosθ(2)

    式中:ELE代表海拔高度,m;LAT代表緯度,(°);θ為冰雪融水區(qū)域以上范圍的平均坡度,(°)。

    1.2 降雨-融水型泥石流流量計算

    (1) 改進算法。

    改進《四川省中小流域暴雨洪水計算手冊》中的計算公式,疊加融水和降雨量,求得全流域在降雨-融水條件下的洪水洪峰流量。

    Q=0.278(H24+M1)24n-1τn-μF(3)

    式中:Q為洪水洪峰流量,m3/s;H24是全流域降雨量,采用24 h降雨量的統(tǒng)計數(shù)值,mm;M1為24 h融水量,mm;n為暴雨參數(shù);τ為流域匯流時間,h,τ=0.278LmJ1/3Q1/4[24];m為匯流參數(shù);J為流域平均縱比降;μ為平均下滲強度,mm/h;F為全流域面積,km2。

    本文采用配方法計算泥石流的峰值流量,其計算公式為

    Qc=(1+φc)QDu(4)

    φc=γc-γwγs-γc(5)

    式中:Qc為泥石流洪峰流量,m3/s;φc為泥石流洪峰流量修正系數(shù);Du為堵塞系數(shù);γc為泥石流容重,t/m3;γw為清水比重,t/m3;γs為泥石流中固體物質(zhì)比重,t/m3。

    (2) 傳統(tǒng)算法。

    目前唯一能考慮冰川泥石流流量計算的經(jīng)驗模型是由中科院成都山地所提出的經(jīng)驗公式[16]:

    Qc=(Q2+Q0)(1+φc)d(6)

    式中:Q0為流域內(nèi)非冰川區(qū)的洪峰流量,m3/s;Q2為降雨型冰川消融流量,m3/s,Q2=F1(0.05H+2.1),H為降雨量,mm;d為冰川消融洪峰系數(shù),d=1+7.6(F1/F)+0.05θ0;F1為冰川面積,km2;θ0為冰川坡度,(°)。

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

    帕隆藏布流域是中國最典型的海洋性冰川區(qū)和雅魯藏布江降水中心且具有雨熱同期的特點,升溫和暴雨成為誘發(fā)大型泥石流的有利條件。本文對波密縣1961~2011年51 a的降雨和氣溫數(shù)據(jù)通過皮爾遜-Ⅲ型曲線進行概率分析統(tǒng)計計算(見圖1~2),在皮爾遜-Ⅲ型曲線的分析中得到在不同頻率下年最高日平均氣溫以及日降雨量如表1所列。

    由于波密縣的氣溫呈現(xiàn)整體上升的趨勢,因此將這51 a的日最高平均氣溫擬合如圖2所示,在對2011年以前的泥石流實例驗證中的溫度指標減去了這個升溫趨勢,公式為TPa=TP-0.03612×(2011-a);TPa為第a年頻率為P的年最高日平均氣溫,℃;TP為波密縣統(tǒng)計不同頻率年最高日平均氣溫,℃;a為年份。對以后泥石流的預測則是在不同情景條件下疊加升溫趨勢,公式為TPn=TP+0.03612×n,TPn為n年一遇的年最高日平均氣溫,℃。

    3 模型驗證

    本文選取了帕隆藏布流域前期研究和實測數(shù)據(jù)相對較豐富的古鄉(xiāng)溝和培龍溝兩條典型冰川泥石流開展實例驗證(見圖3)。歷史統(tǒng)計數(shù)據(jù)如表2所列。

    通過皮爾遜-Ⅲ型曲線計算不同頻率下的溫度TP以及TPa得到1983,1984年和2005年的不同頻率下的最高日平均氣溫如表3所列。

    通過改進算法和傳統(tǒng)算法計算古鄉(xiāng)溝與培龍溝不同頻率下的泥石流數(shù)據(jù)如表4~5所列。

    改進算法和傳統(tǒng)算法在不同頻率下計算得到的古鄉(xiāng)溝和培龍溝泥石流流量與實測數(shù)據(jù)對比如圖4~5所示。

    如表2~3所列,以溫度為衡量指標,古鄉(xiāng)溝2005年約為27 a一遇的泥石流,由改進算法求得的流量為831.33 m3/s,與前人[25]計算的33 a一遇的泥石流規(guī)模849.65 m3/s比較符合。以溫度為衡量指標,培龍溝1983年和1984年的泥石流頻率分別為36 a一遇和154 a一遇,由改進算法求得的流量分別為2 808 m3/s與5 259 m3/s,與歷史記載流量2 950 m3/s與5 245 m3/s相似度較高。而由圖4~5和表4~5對比分析可知,傳統(tǒng)算法的計算結(jié)果普遍偏大,改進算法對古鄉(xiāng)溝和培龍溝的計算準確性相比傳統(tǒng)方法分別提高47.26%與21.12%~32.5%。

    根據(jù)歷史資料揭示,古鄉(xiāng)溝泥石流的規(guī)模在1953年以后在逐漸減小,且20世紀70年代以后泥石流流量都普遍小于1 000 m3/s[28],而傳統(tǒng)經(jīng)驗算法計算古鄉(xiāng)溝20 a一遇流量達1 156 m3/s,嚴重偏大,相比之下改進算法精度更高。

    4 氣候變化條件下泥石流預測

    4.1 波密地區(qū)泥石流溝基本數(shù)據(jù)

    本文選擇了15條典型降雨-融水型泥石流溝進行氣候變化條件下泥石流流量預測(見圖3),基礎數(shù)據(jù)如表6所列。

    4.2 波密地區(qū)泥石流冰雪消融量

    波密地區(qū)的泥石流冰雪融水當量采用度日模型公式(1) 進行計算。以表1的不同頻率下的年最高日平均氣溫為標準,疊加最近51 a升溫趨勢下的日平均氣溫與降雨量,如表7所列。計算得到的全流域24 h冰雪融水當量如表8所列。

    4.3 波密地區(qū)泥石流溝計算結(jié)果

    結(jié)合全流域24 h冰雪融水量和降雨量,通過傳統(tǒng)算法(式(6))和改進算法(式(3)、(4))計算不同頻率下泥石流溝的泥石流流量的結(jié)果如表9和圖6所示。

    將上述15條泥石流溝按冰雪融水面積占全流域面積大于1/4和小于1/4分為兩類(見圖7和圖8)。

    由上述圖表可知,傳統(tǒng)經(jīng)驗算法的計算結(jié)果很大程度上取決于冰川區(qū)域,在冰川面積占比大于1/4的區(qū)域(見圖7),改進算法的計算結(jié)果稍低于偏向保守的傳統(tǒng)算法。而在冰川面積占比小于1/4的區(qū)域(見圖8),傳統(tǒng)算法計算的結(jié)果則遠小于改進算法,因而在冰川占比較小的區(qū)域,傳統(tǒng)算法的結(jié)果則會過于偏小。而現(xiàn)如今冰川的面積在逐年減少[29],所以相比傳統(tǒng)依賴于冰川面積的經(jīng)驗算法,改進算法在今后預測氣候變化條件下泥石流流量的實用性和準確性更強。

    5 結(jié)論和討論

    (1) 建立了降雨-融水型泥石流溫度頻率確定標準,構(gòu)建了波密地區(qū)氣候變化條件下的升溫趨勢,提出了基于改進的度日模型和雨洪法的改進降雨-融水型泥石流流量算法。

    (2) 以古鄉(xiāng)溝和培龍溝為例,改進算法的結(jié)果與歷史數(shù)據(jù)對比,誤差分別為2.24%與0.27%~4.85%,準確性比傳統(tǒng)經(jīng)驗模型提高了47.26%與21.12%~32.5%。

    (3) 對于冰川占比大于1/4的4條溝,傳統(tǒng)算法與改進算法誤差在3%~20%之間,由于傳統(tǒng)算法較為保守,因此傳統(tǒng)算法的結(jié)果均大于本文的改進算法,所以傳統(tǒng)的經(jīng)驗算法對于冰川面積占比大于1/4的泥石流流域仍然具有不錯的參考價值。但是對于冰川面積占比較小的泥石流溝,計算結(jié)果則會偏低甚至無法作為參考,因此改進算法適用性更廣。

    (4) 本研究在用度日模型計算冰雪融水的時候疊加了波密縣近51 a來的升溫趨勢,使預測結(jié)果更加符合氣候變化的趨勢。

    參考文獻:

    [1] 呂儒仁,李德基.西藏波密冬茹弄巴的冰雪融水泥石流[J].冰川凍土,1989,11(2):148-160.

    [2] GE Y G,CUI P,SU F H,et al.Case history of the disastrous debris flows of Tianmo Watershed in Bomi County,Tibet,China:some mitigation suggestions[J].Journal of Mountain Science,2014,11(5):1253-1265.

    [3] JACKSON L E,HUNGR O,GARDNER J S,et al.Cathedral Mountain debris flows,Canada[J].Bulletin of Engineering Geology and the Environment,1989,40(1):35-54.

    [4] WILKERSON F D,SCHMID G L.Distribution of debris flows in Glacier National Park,Montana,U.S.A.[J].Journal of Mountain Science,2008,5(4):318-326.

    [5] HUGGEL C,HAEBERLI W,KAAB A,et al.An assessment procedure for glacial hazards in the Swiss Alps[J].Canadian Geotechnical Journal,2004,41(6):1068-1083.

    [6] 何秋樂,匡星星,梁四海,等.1966~2015年長江源冰川融水變化及其對徑流的影響:以冬克瑪?shù)缀恿饔驗槔齕J].人民長江,2020,51(2):77-85,130.

    [7] 施雅風,楊宗輝,謝自楚,等.西藏古鄉(xiāng)地區(qū)的冰川泥石流[J].科學通報,1964,9(6):542-544.

    [8] 杜榕桓,章書成.西藏高原東南部冰川泥石流的特征[J].冰川凍土,1981,3(3):10-16.

    [9] FLEISCHMAN C M.Debris Flow[M].Beijing:Science Press,1985.

    [10] RUOFU S.Basic characteristics and measuring method of debris flow[M].Chongqing:Science and Technology Literature Press Chongqing Branch,1986.

    [11] 齊云龍,鄧明楓.川藏公路波密段九絨溝泥石流形成機制研究[J].災害學,2019,34(3):123-127.

    [12] 陳曉清,崔鵬,楊忠,等.聶拉木縣沖堆普2002年泥石流成因分析及防治對策[J].冰川凍土,2006,28(5):776-781.

    [13] ZIMMERMANN M,HAEBERLI W.Climatic change and debris flow activity in high-mountain areas:a case study in the Swiss Alps[J].Catena Supplement,1992,22(1):59-72.

    [14] 陳寧生,王政,田樹峰,等.冰磧土體起動泥石流的特征研究[J].第四紀研究,2019,39(5):1235-1245.

    [15] 潘蕾,魏學利,張遠芳,等.初始含水率對冰川泥石流的起動影響分析[J].水土保持學報,2017,31(6):116-122.

    [16] 羅德富,毛濟周.川藏公路南線(西藏境內(nèi))山地災害及防治對策[M].北京:科學出版社,1995.

    [17] 屈永平,肖進,潘義為.藏東南地區(qū)冰川泥石流形成條件初步分析:以天摩溝冰川泥石流為例[J].水利水電技術(shù),2018,49(12):180-187.

    [18] SINGH P,KUMAR N,ARORA M.Degree-day factors for snow and ice for Dokriani Glacier,Garhwal Himalayas[J].Journal of Hydrology,2000,235(1):1-11.

    [19] HOCK R.A distributed temperature-index ice-and snowmelt model including potential direct solar radiation[J].Journal of Glaciology,1999,45(149):101-111.

    [20] HOCK R.Temperature index melt modelling in mountain areas[J].Journal of Hydrology,2003,282(1/4):104-115.

    [21] 鄧明楓.帕隆藏布流域冰川降雨型泥石流頻發(fā)成因分析[D].北京:中國科學院大學,2019.

    [22] 劉金平,張萬昌.雅魯藏布江流域度日因子空間變化[J].中國科學院大學學報,2018,35(5):131-138.

    [23] HOCK R.Glacier melt:a review of processes and their modelling[J].Progress in Physical Geography,2005,29(3):162-391.

    [24] 熊文華,孫立霞,李誠,等.瀾滄江與怒江流域洪峰流量公式推求[J].人民長江,2020,51(11):82-86.

    [25] 魯安新,鄧曉峰,趙尚學,等.2005年西藏波密古鄉(xiāng)溝泥石流暴發(fā)成因分析[J].冰川凍土,2006,28(6):956-960.

    [26] 程尊蘭,吳積善.西藏東南部培龍溝泥石流堵塞壩的形成機理[A].第八屆海峽兩岸山地災害與環(huán)境保育學術(shù)研討會[C]∥烏魯木齊:中國水土保持學會,2011:39-45.

    [27] 國家防汛總指揮部辦公室,中國科學院水利部成都山地災害與環(huán)境研究所.山洪泥石流滑坡災害及防治[M].北京:科學出版社,1994.

    [28] 朱平一,羅德富.西藏古鄉(xiāng)溝泥石流發(fā)展趨勢[J].山地學報,1997,15(4):296-299.

    [29] 劉時銀,姚曉軍,郭萬欽,等.基于第二次冰川編目的中國冰川現(xiàn)狀[J].地理學報,2015,70(1):3-16.

    (編輯:劉 媛)

    Discharge prediction method for rainfall-melting water triggered glacial debris flows under climate change scenarios

    YANG Zongji1,F(xiàn)U Xiaolong1,2,YOU Yong1

    (1.Institute of Mountain Hazards and Environment,CAS,Chengdu 610041,China; 2.University of Chinese Academy of Sciences,Beijing 100049,China)

    Abstract:

    Under the background of global warming,current empirical model for the glacial debris flow discharge calculation cannot quantitatively consider the coupling effect of rainfall water and the snow/ice melting water.We proposed a discharge calculation model that considers coupling effect of both rain and melting water for glacial debris flows by studying ice-snow melt equivalent in the improved degree-day model as well as the current storm flood calculation method for rainfall triggered debris flows.The criteria of return period for the glacial debris flows discharge prediction was established based on the rainfall and temperature change trends fitted by the 51 years long period observation data in the Palong Zangbo watershed.Then this proposed model was verified by the observation data of 2 typical glacial debris flows along the Palong Zangbo River,and the accuracy of the proposed method was significantly improved by comparing to the current empirical model.Afterwards,the discharge estimation of 15 glacial debris flow gullies in the Palong Zangbo watershed were analyzed under various climate change scenarios.The results show that the proposed method can quantify the coupling effect of rain-thermal on debris flow and can provide a reference for the calculation and prediction of key parameters for the prevention and control of glacial debris flows.

    Key words:

    rainfall-melting water debris flows;debris flows discharge;degree-day model;rain-thermal coupling changes;climate change;Palong Zangbo watershed

    猜你喜歡
    融水冰川泥石流
    我在融水唱貝江
    融水糯米柚種植栽培管理與技術(shù)探究
    為什么冰川會到處走?
    冰川會發(fā)出聲音嗎?
    泥石流
    雜文月刊(2018年21期)2019-01-05 05:55:28
    長途跋涉到冰川
    例舉高中數(shù)學解題切入點的找尋
    “民謠泥石流”花粥:唱出自己
    海峽姐妹(2017年6期)2017-06-24 09:37:36
    泥石流
    機械班長
    国产在线精品亚洲第一网站| 女性生殖器流出的白浆| www.999成人在线观看| 深夜精品福利| 十八禁网站免费在线| 国产精品.久久久| 色精品久久人妻99蜜桃| 中亚洲国语对白在线视频| 午夜精品在线福利| 成人黄色视频免费在线看| 欧美乱妇无乱码| 黄色视频不卡| 久久国产乱子伦精品免费另类| 亚洲欧美一区二区三区久久| 99精品久久久久人妻精品| 丰满的人妻完整版| 国产精品乱码一区二三区的特点 | 久久久久视频综合| 叶爱在线成人免费视频播放| 9热在线视频观看99| 99香蕉大伊视频| 国产精品一区二区在线观看99| 亚洲欧美精品综合一区二区三区| 精品国产乱子伦一区二区三区| av网站免费在线观看视频| a级片在线免费高清观看视频| 欧美一级毛片孕妇| 国产成人精品无人区| 麻豆国产av国片精品| 一进一出抽搐动态| 国产区一区二久久| 最新在线观看一区二区三区| 成人国语在线视频| 成年女人毛片免费观看观看9 | 搡老熟女国产l中国老女人| 国产乱人伦免费视频| 久久久久久久国产电影| 国产真人三级小视频在线观看| 少妇被粗大的猛进出69影院| 老汉色av国产亚洲站长工具| 日韩欧美免费精品| 亚洲国产精品一区二区三区在线| 亚洲国产欧美一区二区综合| 在线观看免费日韩欧美大片| tocl精华| 18禁黄网站禁片午夜丰满| 亚洲成a人片在线一区二区| 久久久久国产一级毛片高清牌| 国产激情久久老熟女| av一本久久久久| 久久久久久久久免费视频了| 日日爽夜夜爽网站| 国产在视频线精品| 一边摸一边抽搐一进一出视频| 黄色成人免费大全| 妹子高潮喷水视频| 欧美日韩亚洲高清精品| 韩国精品一区二区三区| 自线自在国产av| 国产深夜福利视频在线观看| 精品亚洲成a人片在线观看| 一级a爱片免费观看的视频| 欧美av亚洲av综合av国产av| 久久精品熟女亚洲av麻豆精品| 久久天堂一区二区三区四区| 成人亚洲精品一区在线观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产毛片av蜜桃av| 国产一区在线观看成人免费| 国产精品久久久人人做人人爽| 中文字幕人妻丝袜一区二区| 国产亚洲欧美在线一区二区| 黄频高清免费视频| 人成视频在线观看免费观看| 亚洲精品自拍成人| 69精品国产乱码久久久| 午夜精品久久久久久毛片777| 18禁黄网站禁片午夜丰满| 两人在一起打扑克的视频| 校园春色视频在线观看| 中文字幕最新亚洲高清| av片东京热男人的天堂| 国产国语露脸激情在线看| 高清视频免费观看一区二区| 人妻一区二区av| 99国产综合亚洲精品| 欧美成狂野欧美在线观看| 中文欧美无线码| 国产精品 国内视频| 日韩熟女老妇一区二区性免费视频| 国产淫语在线视频| 人人妻,人人澡人人爽秒播| 精品久久久精品久久久| 亚洲精品乱久久久久久| 欧美日韩福利视频一区二区| 人人妻人人澡人人爽人人夜夜| 午夜福利视频在线观看免费| 欧美人与性动交α欧美软件| 亚洲avbb在线观看| a级毛片在线看网站| 99热只有精品国产| 免费在线观看影片大全网站| 免费一级毛片在线播放高清视频 | 成人手机av| 精品国产乱码久久久久久男人| 国产亚洲欧美在线一区二区| av有码第一页| 免费看a级黄色片| 久久久久国产精品人妻aⅴ院 | 亚洲七黄色美女视频| 国产精品成人在线| 亚洲aⅴ乱码一区二区在线播放 | 老熟妇仑乱视频hdxx| 丰满饥渴人妻一区二区三| 亚洲国产欧美日韩在线播放| 国产免费av片在线观看野外av| 啦啦啦 在线观看视频| 高清黄色对白视频在线免费看| 国产欧美日韩综合在线一区二区| 91麻豆av在线| 精品久久久久久电影网| 老司机午夜十八禁免费视频| 亚洲人成电影免费在线| 国产精品 国内视频| 色94色欧美一区二区| 久久久精品免费免费高清| 久久久水蜜桃国产精品网| 亚洲中文av在线| 香蕉丝袜av| 国产成人影院久久av| 欧美色视频一区免费| 女人精品久久久久毛片| 大码成人一级视频| 国产亚洲欧美精品永久| 久久草成人影院| 淫妇啪啪啪对白视频| 免费观看精品视频网站| 国产成人精品久久二区二区免费| 国产成人精品久久二区二区免费| 国产精品电影一区二区三区 | 老熟妇仑乱视频hdxx| av欧美777| 99在线人妻在线中文字幕 | 看免费av毛片| 老司机深夜福利视频在线观看| 老司机深夜福利视频在线观看| av免费在线观看网站| 丰满饥渴人妻一区二区三| 热re99久久精品国产66热6| 老司机靠b影院| 国产99白浆流出| 欧美日韩亚洲高清精品| 人人妻,人人澡人人爽秒播| 久久久精品免费免费高清| 日本a在线网址| 精品一品国产午夜福利视频| 国产亚洲精品第一综合不卡| 日韩欧美国产一区二区入口| 999久久久精品免费观看国产| 超碰97精品在线观看| 极品少妇高潮喷水抽搐| a级毛片黄视频| 丝瓜视频免费看黄片| 久热这里只有精品99| 欧美黑人欧美精品刺激| 免费观看人在逋| 91麻豆精品激情在线观看国产 | 天天影视国产精品| 成年人免费黄色播放视频| 99国产精品一区二区蜜桃av | 亚洲av成人不卡在线观看播放网| 黑人巨大精品欧美一区二区蜜桃| 日韩欧美三级三区| 国产精品一区二区在线观看99| 国产精品久久久久久人妻精品电影| 1024视频免费在线观看| 亚洲成a人片在线一区二区| 亚洲在线自拍视频| 热99国产精品久久久久久7| 中文亚洲av片在线观看爽 | 国产免费av片在线观看野外av| 日本精品一区二区三区蜜桃| 亚洲七黄色美女视频| 久久性视频一级片| xxx96com| 在线观看66精品国产| 丰满的人妻完整版| 99riav亚洲国产免费| 精品久久蜜臀av无| 热99国产精品久久久久久7| 亚洲av成人不卡在线观看播放网| 男女床上黄色一级片免费看| 18禁美女被吸乳视频| 99国产综合亚洲精品| 在线观看日韩欧美| 亚洲精品国产一区二区精华液| 丁香六月欧美| 每晚都被弄得嗷嗷叫到高潮| 99久久人妻综合| 99精品欧美一区二区三区四区| 性色av乱码一区二区三区2| 成人av一区二区三区在线看| 黑人巨大精品欧美一区二区蜜桃| 精品人妻熟女毛片av久久网站| 18在线观看网站| 性少妇av在线| 一区二区三区国产精品乱码| 久久久久久免费高清国产稀缺| 国产高清视频在线播放一区| 黄网站色视频无遮挡免费观看| 成人av一区二区三区在线看| 在线天堂中文资源库| 亚洲一区二区三区欧美精品| 三级毛片av免费| 我的亚洲天堂| 欧美日韩av久久| 日韩精品免费视频一区二区三区| 狠狠婷婷综合久久久久久88av| av片东京热男人的天堂| 国产成+人综合+亚洲专区| 精品久久久久久,| 久久人妻福利社区极品人妻图片| 亚洲人成伊人成综合网2020| 看黄色毛片网站| 欧美亚洲日本最大视频资源| av视频免费观看在线观看| 国产淫语在线视频| 国产有黄有色有爽视频| 国产单亲对白刺激| 国产在视频线精品| 久久久久久免费高清国产稀缺| 精品人妻1区二区| 五月开心婷婷网| 久久精品国产亚洲av高清一级| bbb黄色大片| 怎么达到女性高潮| 亚洲欧美日韩另类电影网站| 欧美黄色片欧美黄色片| 午夜亚洲福利在线播放| 人妻久久中文字幕网| 黄色 视频免费看| 久久中文字幕人妻熟女| 人人妻人人添人人爽欧美一区卜| tube8黄色片| 亚洲五月婷婷丁香| 777久久人妻少妇嫩草av网站| 亚洲精品国产精品久久久不卡| 少妇的丰满在线观看| 久久国产精品影院| 欧美日韩亚洲综合一区二区三区_| 18禁裸乳无遮挡动漫免费视频| 国产亚洲一区二区精品| av超薄肉色丝袜交足视频| 丁香六月欧美| 日韩欧美三级三区| 高清在线国产一区| 亚洲欧美色中文字幕在线| 日本a在线网址| 久久 成人 亚洲| 一a级毛片在线观看| 成人手机av| 热99国产精品久久久久久7| 久久香蕉国产精品| xxx96com| 黄色女人牲交| 亚洲国产中文字幕在线视频| 午夜影院日韩av| 午夜精品在线福利| 日韩中文字幕欧美一区二区| 母亲3免费完整高清在线观看| 日日夜夜操网爽| 一进一出好大好爽视频| 精品一区二区三区四区五区乱码| 天堂俺去俺来也www色官网| 国产精品av久久久久免费| 成人影院久久| 久久狼人影院| 欧美黄色淫秽网站| 午夜精品国产一区二区电影| 后天国语完整版免费观看| 国产精品久久久久成人av| av在线播放免费不卡| 精品久久久久久电影网| 国产精品一区二区精品视频观看| 人人妻人人爽人人添夜夜欢视频| 亚洲欧洲精品一区二区精品久久久| 村上凉子中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 侵犯人妻中文字幕一二三四区| 亚洲av日韩精品久久久久久密| www.熟女人妻精品国产| 国产又色又爽无遮挡免费看| 久久久国产成人精品二区 | 国产高清视频在线播放一区| 中亚洲国语对白在线视频| 国产亚洲精品一区二区www | 欧美日韩精品网址| 午夜福利在线免费观看网站| 高潮久久久久久久久久久不卡| av中文乱码字幕在线| 亚洲熟女精品中文字幕| 9热在线视频观看99| 少妇裸体淫交视频免费看高清 | 欧美性长视频在线观看| 国产精品一区二区在线观看99| 久久人人爽av亚洲精品天堂| 黄片小视频在线播放| 国产区一区二久久| 亚洲伊人色综图| 国产欧美日韩一区二区精品| 黄网站色视频无遮挡免费观看| 久久久久久久精品吃奶| 欧美精品啪啪一区二区三区| 极品教师在线免费播放| 一级毛片高清免费大全| 色94色欧美一区二区| 50天的宝宝边吃奶边哭怎么回事| 在线免费观看的www视频| 国产精品 欧美亚洲| tube8黄色片| 精品国产一区二区久久| 国产精品影院久久| 热99re8久久精品国产| 国产无遮挡羞羞视频在线观看| av网站免费在线观看视频| 不卡av一区二区三区| 精品久久久久久久毛片微露脸| 精品一区二区三区av网在线观看| 亚洲av美国av| 丰满人妻熟妇乱又伦精品不卡| 国产一区有黄有色的免费视频| e午夜精品久久久久久久| 老熟妇仑乱视频hdxx| 国产av又大| 香蕉久久夜色| 免费高清在线观看日韩| 丝袜美足系列| 美女高潮到喷水免费观看| 欧美大码av| 啦啦啦免费观看视频1| 国产成人欧美在线观看 | 午夜免费成人在线视频| 欧美日韩乱码在线| 日本撒尿小便嘘嘘汇集6| 一级a爱片免费观看的视频| 成人免费观看视频高清| 国产亚洲欧美98| 亚洲av成人一区二区三| 日韩欧美一区视频在线观看| 视频在线观看一区二区三区| 国产成人av激情在线播放| 男人舔女人的私密视频| 亚洲男人天堂网一区| 久久精品91无色码中文字幕| 少妇猛男粗大的猛烈进出视频| 国产一区在线观看成人免费| 啪啪无遮挡十八禁网站| 淫妇啪啪啪对白视频| 亚洲美女黄片视频| 老司机影院毛片| 夜夜夜夜夜久久久久| 亚洲中文字幕日韩| 亚洲精品国产精品久久久不卡| 久久久国产一区二区| 身体一侧抽搐| 亚洲性夜色夜夜综合| 国产片内射在线| 9热在线视频观看99| 黄色视频,在线免费观看| 欧美日韩亚洲国产一区二区在线观看 | www.熟女人妻精品国产| 高潮久久久久久久久久久不卡| 精品久久久久久,| 欧美日韩视频精品一区| 国产成+人综合+亚洲专区| 99精品欧美一区二区三区四区| 亚洲综合色网址| 久久久精品区二区三区| 少妇裸体淫交视频免费看高清 | 国内毛片毛片毛片毛片毛片| 欧美丝袜亚洲另类 | 亚洲成人国产一区在线观看| 亚洲aⅴ乱码一区二区在线播放 | svipshipincom国产片| 老司机靠b影院| 中文字幕人妻丝袜一区二区| 久久人人97超碰香蕉20202| 天堂中文最新版在线下载| 亚洲欧美日韩高清在线视频| 极品少妇高潮喷水抽搐| 亚洲精品自拍成人| 我的亚洲天堂| 亚洲国产精品sss在线观看 | 日韩免费高清中文字幕av| 看片在线看免费视频| 国产欧美日韩一区二区三区在线| 免费人成视频x8x8入口观看| 人人妻人人添人人爽欧美一区卜| 老熟妇仑乱视频hdxx| а√天堂www在线а√下载 | videos熟女内射| 99热国产这里只有精品6| 一级黄色大片毛片| 99国产综合亚洲精品| 脱女人内裤的视频| 在线永久观看黄色视频| 99久久精品国产亚洲精品| 亚洲欧美一区二区三区久久| 成年人免费黄色播放视频| 精品久久久久久久久久免费视频 | 国产精品国产高清国产av | 久久ye,这里只有精品| 午夜福利免费观看在线| 丁香欧美五月| 超色免费av| 丝袜美足系列| 国产精品国产av在线观看| 免费在线观看亚洲国产| 精品久久久精品久久久| 伦理电影免费视频| 亚洲精品一卡2卡三卡4卡5卡| 一区二区日韩欧美中文字幕| 日本vs欧美在线观看视频| 久久国产精品影院| 成在线人永久免费视频| 成人18禁在线播放| 丝瓜视频免费看黄片| 精品卡一卡二卡四卡免费| 啦啦啦 在线观看视频| 亚洲av熟女| 一级毛片高清免费大全| 国产亚洲欧美在线一区二区| 精品久久久久久久久久免费视频 | 中文字幕精品免费在线观看视频| 国产精品二区激情视频| 国内久久婷婷六月综合欲色啪| 国产亚洲精品第一综合不卡| 99香蕉大伊视频| 午夜两性在线视频| 亚洲五月天丁香| 99国产精品一区二区三区| 大陆偷拍与自拍| 国产激情欧美一区二区| 午夜免费观看网址| 精品久久久久久电影网| 欧美日韩瑟瑟在线播放| 女性生殖器流出的白浆| 久久久久久亚洲精品国产蜜桃av| svipshipincom国产片| 精品久久蜜臀av无| 国产精品亚洲一级av第二区| 欧美激情极品国产一区二区三区| 国产片内射在线| 丁香欧美五月| 免费在线观看黄色视频的| 国产亚洲精品久久久久久毛片 | 精品第一国产精品| 精品国产一区二区三区四区第35| av福利片在线| 成人特级黄色片久久久久久久| 叶爱在线成人免费视频播放| av天堂久久9| 精品视频人人做人人爽| 久热这里只有精品99| 亚洲一区高清亚洲精品| 国产高清视频在线播放一区| 18在线观看网站| 国产成人av教育| 女人精品久久久久毛片| 欧美成人午夜精品| 俄罗斯特黄特色一大片| 在线国产一区二区在线| 国产免费现黄频在线看| 脱女人内裤的视频| 亚洲国产精品sss在线观看 | 亚洲午夜精品一区,二区,三区| 免费高清在线观看日韩| 精品久久久久久久久久免费视频 | 999久久久精品免费观看国产| 在线观看午夜福利视频| 97人妻天天添夜夜摸| 国产伦人伦偷精品视频| 亚洲中文日韩欧美视频| 日本黄色日本黄色录像| 精品国内亚洲2022精品成人 | 亚洲色图 男人天堂 中文字幕| 欧美日韩精品网址| 少妇裸体淫交视频免费看高清 | 午夜老司机福利片| 精品卡一卡二卡四卡免费| 色94色欧美一区二区| 亚洲精品久久午夜乱码| 国产午夜精品久久久久久| 黄色毛片三级朝国网站| 黄色女人牲交| aaaaa片日本免费| 国产欧美日韩一区二区三区在线| 一级片免费观看大全| 一级毛片高清免费大全| 午夜福利,免费看| 9色porny在线观看| 黄色女人牲交| 国产成人免费观看mmmm| 老熟女久久久| 欧美日韩视频精品一区| 欧美亚洲日本最大视频资源| 1024视频免费在线观看| 一二三四在线观看免费中文在| 人成视频在线观看免费观看| 成人精品一区二区免费| 大香蕉久久网| 成人精品一区二区免费| 久久精品人人爽人人爽视色| 亚洲中文日韩欧美视频| 黄色视频,在线免费观看| 欧美日韩国产mv在线观看视频| 欧美激情 高清一区二区三区| 日本wwww免费看| 中文字幕制服av| 日本一区二区免费在线视频| 亚洲精品成人av观看孕妇| 国产精品二区激情视频| 女性生殖器流出的白浆| 黑人猛操日本美女一级片| 中文字幕最新亚洲高清| 欧美亚洲 丝袜 人妻 在线| 男男h啪啪无遮挡| 99在线人妻在线中文字幕 | av线在线观看网站| 一区二区三区精品91| 欧美久久黑人一区二区| 成年动漫av网址| 搡老乐熟女国产| 在线观看免费视频网站a站| 老司机午夜福利在线观看视频| 精品一区二区三卡| 国产片内射在线| 久久精品国产亚洲av香蕉五月 | 搡老岳熟女国产| 欧洲精品卡2卡3卡4卡5卡区| 久热这里只有精品99| 欧美精品av麻豆av| 亚洲av成人av| 成人影院久久| 欧美激情 高清一区二区三区| 搡老乐熟女国产| 99热只有精品国产| 黑人猛操日本美女一级片| 人妻 亚洲 视频| 夜夜夜夜夜久久久久| 成人黄色视频免费在线看| 欧美国产精品va在线观看不卡| 性色av乱码一区二区三区2| 美女国产高潮福利片在线看| 一级片免费观看大全| 麻豆乱淫一区二区| 超色免费av| 每晚都被弄得嗷嗷叫到高潮| 99久久精品国产亚洲精品| 久久久国产欧美日韩av| 亚洲成a人片在线一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人欧美在线观看 | 久9热在线精品视频| videos熟女内射| 国产激情久久老熟女| 欧美黑人精品巨大| 国产高清videossex| 在线观看免费视频网站a站| 精品第一国产精品| 欧美日本中文国产一区发布| 精品国产乱子伦一区二区三区| 日韩熟女老妇一区二区性免费视频| 看片在线看免费视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美国产精品va在线观看不卡| 亚洲一码二码三码区别大吗| 国产亚洲精品久久久久久毛片 | 动漫黄色视频在线观看| 国产欧美日韩一区二区精品| 国产高清视频在线播放一区| 757午夜福利合集在线观看| av超薄肉色丝袜交足视频| 美女国产高潮福利片在线看| 亚洲av第一区精品v没综合| 亚洲精品一卡2卡三卡4卡5卡| 免费人成视频x8x8入口观看| 国产黄色免费在线视频| 男男h啪啪无遮挡| 亚洲一码二码三码区别大吗| 亚洲中文字幕日韩| 国产精华一区二区三区| 97人妻天天添夜夜摸| 一级,二级,三级黄色视频| 热99国产精品久久久久久7| 国产有黄有色有爽视频| 亚洲自偷自拍图片 自拍| 国产精品九九99| 757午夜福利合集在线观看| 中亚洲国语对白在线视频| 欧美国产精品va在线观看不卡| 亚洲五月婷婷丁香| 免费在线观看影片大全网站| 色婷婷久久久亚洲欧美| 成人国语在线视频| 国产在线一区二区三区精| av网站在线播放免费| 国产欧美日韩一区二区三| 亚洲 欧美一区二区三区| avwww免费| 夜夜爽天天搞| 国产亚洲av高清不卡| 日韩有码中文字幕| 侵犯人妻中文字幕一二三四区| 中文字幕人妻丝袜制服|