• <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
    泥石流
    機械班長
    超碰97精品在线观看| 男女床上黄色一级片免费看| 丝袜美腿诱惑在线| 只有这里有精品99| 少妇人妻久久综合中文| 交换朋友夫妻互换小说| 啦啦啦在线观看免费高清www| 久久免费观看电影| 一区二区三区精品91| 国产97色在线日韩免费| 国产精品人妻久久久影院| 欧美黄色淫秽网站| 美女主播在线视频| 国产在线免费精品| 99国产精品免费福利视频| 黄片小视频在线播放| 人体艺术视频欧美日本| 宅男免费午夜| 捣出白浆h1v1| 久久鲁丝午夜福利片| 成年人免费黄色播放视频| 免费高清在线观看日韩| 久久精品国产亚洲av高清一级| 亚洲欧美精品自产自拍| 免费女性裸体啪啪无遮挡网站| 亚洲欧美精品自产自拍| 免费看十八禁软件| 日韩大码丰满熟妇| √禁漫天堂资源中文www| 国产精品久久久久成人av| 亚洲精品一区蜜桃| 丝袜美腿诱惑在线| 爱豆传媒免费全集在线观看| 亚洲中文日韩欧美视频| 好男人电影高清在线观看| 国产高清视频在线播放一区 | 麻豆乱淫一区二区| avwww免费| 国产有黄有色有爽视频| 日韩大片免费观看网站| 国产精品久久久久成人av| videos熟女内射| 人人妻人人澡人人爽人人夜夜| 国产欧美日韩综合在线一区二区| 久久国产亚洲av麻豆专区| 亚洲精品国产av蜜桃| 欧美黄色片欧美黄色片| 国产色视频综合| 国产高清videossex| 水蜜桃什么品种好| 久久久欧美国产精品| 十八禁高潮呻吟视频| 男女床上黄色一级片免费看| 久久99精品国语久久久| 国产精品一区二区在线不卡| 国产日韩欧美视频二区| 亚洲成人国产一区在线观看 | 亚洲av综合色区一区| 免费在线观看完整版高清| 亚洲精品久久久久久婷婷小说| 精品一区在线观看国产| 亚洲三区欧美一区| 精品少妇黑人巨大在线播放| 欧美人与性动交α欧美精品济南到| 日本av手机在线免费观看| 久久精品久久久久久久性| 另类亚洲欧美激情| 国产一区二区三区av在线| 我的亚洲天堂| 在线亚洲精品国产二区图片欧美| 波野结衣二区三区在线| 免费不卡黄色视频| 午夜福利影视在线免费观看| 少妇精品久久久久久久| 欧美xxⅹ黑人| www.精华液| 日韩 亚洲 欧美在线| 考比视频在线观看| 久久九九热精品免费| xxxhd国产人妻xxx| 国产在线一区二区三区精| 岛国毛片在线播放| 日韩一卡2卡3卡4卡2021年| 国语对白做爰xxxⅹ性视频网站| av在线app专区| 99国产综合亚洲精品| 两人在一起打扑克的视频| 一本综合久久免费| 男女无遮挡免费网站观看| 亚洲av欧美aⅴ国产| www.精华液| 少妇 在线观看| 免费看不卡的av| 亚洲精品久久久久久婷婷小说| 制服人妻中文乱码| 蜜桃国产av成人99| 久久久国产精品麻豆| 精品久久久久久电影网| 亚洲精品国产av成人精品| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一码二码三码区别大吗| 男女边吃奶边做爰视频| 一级黄片播放器| 国产一卡二卡三卡精品| 欧美精品一区二区免费开放| 国产免费视频播放在线视频| 亚洲av在线观看美女高潮| 99国产综合亚洲精品| 啦啦啦在线观看免费高清www| 亚洲欧美一区二区三区国产| 欧美黄色片欧美黄色片| av网站免费在线观看视频| 久久99精品国语久久久| 日韩制服骚丝袜av| 中文字幕av电影在线播放| 久久人人97超碰香蕉20202| 午夜福利在线免费观看网站| 各种免费的搞黄视频| 男人爽女人下面视频在线观看| av视频免费观看在线观看| 欧美亚洲 丝袜 人妻 在线| 黄色怎么调成土黄色| 中文字幕制服av| 国产男女超爽视频在线观看| 成人三级做爰电影| 亚洲成人免费电影在线观看 | 久久精品国产综合久久久| 在线观看一区二区三区激情| 日韩大码丰满熟妇| 日日爽夜夜爽网站| tube8黄色片| 国产成人精品无人区| 成人影院久久| 亚洲国产成人一精品久久久| 精品国产乱码久久久久久男人| 高清av免费在线| 日韩中文字幕欧美一区二区 | 亚洲精品日本国产第一区| 人妻 亚洲 视频| 国产精品国产三级专区第一集| 50天的宝宝边吃奶边哭怎么回事| 精品人妻1区二区| 久久精品国产亚洲av高清一级| 欧美国产精品一级二级三级| 国产精品久久久久久精品古装| 国产欧美亚洲国产| 可以免费在线观看a视频的电影网站| 老司机在亚洲福利影院| 亚洲久久久国产精品| 人人妻人人澡人人看| 国产在视频线精品| 欧美激情 高清一区二区三区| 国产精品一区二区在线观看99| 人人妻人人澡人人看| 欧美 亚洲 国产 日韩一| 精品一区二区三卡| 免费人妻精品一区二区三区视频| 亚洲中文日韩欧美视频| 久热这里只有精品99| 久久性视频一级片| 精品久久久精品久久久| 久久久欧美国产精品| 日韩伦理黄色片| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产成人一精品久久久| 满18在线观看网站| 我的亚洲天堂| 亚洲欧美色中文字幕在线| 色精品久久人妻99蜜桃| 99国产精品99久久久久| 巨乳人妻的诱惑在线观看| 亚洲精品日韩在线中文字幕| 男男h啪啪无遮挡| 午夜影院在线不卡| 如日韩欧美国产精品一区二区三区| 亚洲中文日韩欧美视频| 久久久精品免费免费高清| 婷婷色综合大香蕉| 不卡av一区二区三区| 18禁国产床啪视频网站| 菩萨蛮人人尽说江南好唐韦庄| 亚洲午夜精品一区,二区,三区| 欧美黄色片欧美黄色片| 久久久久视频综合| 亚洲精品国产区一区二| 最近最新中文字幕大全免费视频 | 美女高潮到喷水免费观看| 国产精品人妻久久久影院| 建设人人有责人人尽责人人享有的| av有码第一页| 亚洲欧洲日产国产| 只有这里有精品99| 99国产精品一区二区三区| 免费在线观看完整版高清| 99热网站在线观看| 秋霞在线观看毛片| 看免费成人av毛片| 中文精品一卡2卡3卡4更新| 日韩中文字幕欧美一区二区 | 一级a爱视频在线免费观看| 欧美人与善性xxx| 99久久精品国产亚洲精品| 侵犯人妻中文字幕一二三四区| 国产欧美亚洲国产| www.精华液| 中文字幕亚洲精品专区| 久久免费观看电影| 人妻人人澡人人爽人人| 91精品三级在线观看| 欧美黄色淫秽网站| 国产高清视频在线播放一区 | 啦啦啦在线观看免费高清www| 婷婷成人精品国产| cao死你这个sao货| 少妇猛男粗大的猛烈进出视频| 美女午夜性视频免费| 后天国语完整版免费观看| 我要看黄色一级片免费的| 亚洲,欧美精品.| 国产成人免费观看mmmm| 婷婷成人精品国产| 色精品久久人妻99蜜桃| 一级毛片我不卡| 国产一卡二卡三卡精品| 91麻豆av在线| 天天躁夜夜躁狠狠躁躁| 丝袜人妻中文字幕| 老司机深夜福利视频在线观看 | 国产一区二区三区av在线| 欧美性长视频在线观看| 国产在线免费精品| av天堂久久9| 在线观看人妻少妇| 日本色播在线视频| 一级a爱视频在线免费观看| 夫妻性生交免费视频一级片| 久久综合国产亚洲精品| 欧美精品高潮呻吟av久久| 国产一区二区三区av在线| 欧美国产精品va在线观看不卡| 精品少妇黑人巨大在线播放| 欧美另类一区| 日本91视频免费播放| 欧美日韩视频精品一区| 亚洲精品日本国产第一区| av有码第一页| 久久综合国产亚洲精品| 大陆偷拍与自拍| 亚洲精品乱久久久久久| 免费在线观看影片大全网站 | av电影中文网址| 性高湖久久久久久久久免费观看| av片东京热男人的天堂| 国产欧美日韩一区二区三区在线| 欧美日韩成人在线一区二区| 欧美人与性动交α欧美软件| 亚洲av电影在线进入| 久久ye,这里只有精品| av在线老鸭窝| 亚洲国产av新网站| 在线观看www视频免费| 亚洲一码二码三码区别大吗| 国语对白做爰xxxⅹ性视频网站| 99国产精品免费福利视频| 国产男女内射视频| 免费黄频网站在线观看国产| 少妇人妻 视频| 国产日韩欧美视频二区| 国产精品 国内视频| 亚洲av片天天在线观看| 黄色一级大片看看| 久久热在线av| 国产午夜精品一二区理论片| xxx大片免费视频| 国产欧美日韩一区二区三区在线| 欧美精品一区二区免费开放| 秋霞在线观看毛片| 亚洲国产欧美一区二区综合| 99国产精品一区二区三区| 看十八女毛片水多多多| 精品少妇久久久久久888优播| 精品视频人人做人人爽| 久久久久国产一级毛片高清牌| 欧美+亚洲+日韩+国产| 大香蕉久久成人网| 亚洲一区中文字幕在线| 日本a在线网址| 国产成人a∨麻豆精品| 深夜精品福利| 中文字幕最新亚洲高清| h视频一区二区三区| 亚洲伊人久久精品综合| 亚洲精品国产av蜜桃| 在线观看免费日韩欧美大片| 亚洲专区国产一区二区| 亚洲国产欧美一区二区综合| 手机成人av网站| 99精国产麻豆久久婷婷| 丝袜美足系列| 美女高潮到喷水免费观看| 黑人欧美特级aaaaaa片| 欧美老熟妇乱子伦牲交| 欧美国产精品一级二级三级| 精品一品国产午夜福利视频| 大码成人一级视频| 在线观看人妻少妇| 久久久久久亚洲精品国产蜜桃av| av在线app专区| 女警被强在线播放| 国产爽快片一区二区三区| 亚洲五月婷婷丁香| 国产亚洲一区二区精品| 嫁个100分男人电影在线观看 | 69精品国产乱码久久久| 精品国产超薄肉色丝袜足j| 精品福利永久在线观看| 9热在线视频观看99| 色网站视频免费| 亚洲 欧美一区二区三区| 18禁国产床啪视频网站| 黑人猛操日本美女一级片| 国产真人三级小视频在线观看| 国产一区二区在线观看av| 国产高清视频在线播放一区 | 免费在线观看日本一区| 一二三四社区在线视频社区8| 精品一区在线观看国产| 赤兔流量卡办理| 国产精品免费大片| av欧美777| 一本色道久久久久久精品综合| 自线自在国产av| 亚洲天堂av无毛| 亚洲中文字幕日韩| 国产黄频视频在线观看| 黄网站色视频无遮挡免费观看| av网站免费在线观看视频| 欧美精品啪啪一区二区三区 | 嫁个100分男人电影在线观看 | 亚洲欧洲国产日韩| 乱人伦中国视频| 国产国语露脸激情在线看| 少妇粗大呻吟视频| 别揉我奶头~嗯~啊~动态视频 | 99re6热这里在线精品视频| 在线观看国产h片| 久久精品aⅴ一区二区三区四区| 七月丁香在线播放| 午夜福利视频在线观看免费| 中文字幕色久视频| 国产真人三级小视频在线观看| 高清av免费在线| 国产亚洲av片在线观看秒播厂| 日韩一卡2卡3卡4卡2021年| 丝瓜视频免费看黄片| 99精品久久久久人妻精品| 免费看av在线观看网站| 99香蕉大伊视频| 飞空精品影院首页| 亚洲欧洲国产日韩| 男女午夜视频在线观看| 精品久久蜜臀av无| 搡老岳熟女国产| 最新在线观看一区二区三区 | 婷婷色综合大香蕉| 伊人久久大香线蕉亚洲五| 欧美日韩一级在线毛片| 18禁观看日本| 国产精品秋霞免费鲁丝片| 三上悠亚av全集在线观看| 最黄视频免费看| 在线亚洲精品国产二区图片欧美| 久久这里只有精品19| 99久久人妻综合| 亚洲欧美成人综合另类久久久| 少妇粗大呻吟视频| 日韩制服骚丝袜av| 欧美av亚洲av综合av国产av| 亚洲精品国产一区二区精华液| 美女中出高潮动态图| 自拍欧美九色日韩亚洲蝌蚪91| 欧美av亚洲av综合av国产av| 日韩av免费高清视频| 国产主播在线观看一区二区 | 国产三级黄色录像| 亚洲欧美一区二区三区久久| 欧美日韩亚洲综合一区二区三区_| netflix在线观看网站| 国产成人精品在线电影| 又黄又粗又硬又大视频| 男女免费视频国产| 一二三四在线观看免费中文在| 亚洲精品日韩在线中文字幕| 亚洲精品国产区一区二| 视频区图区小说| 欧美亚洲日本最大视频资源| 国产成人一区二区三区免费视频网站 | 国产成人一区二区在线| 脱女人内裤的视频| 老司机影院毛片| 一本久久精品| 人成视频在线观看免费观看| 国产精品免费视频内射| 成年人黄色毛片网站| 欧美 日韩 精品 国产| 国产极品粉嫩免费观看在线| 精品人妻1区二区| 一本色道久久久久久精品综合| 巨乳人妻的诱惑在线观看| 美女扒开内裤让男人捅视频| 亚洲成人国产一区在线观看 | 久久中文字幕一级| 在现免费观看毛片| 日韩一本色道免费dvd| 亚洲精品成人av观看孕妇| 久久性视频一级片| 久久久久久久大尺度免费视频| 欧美日韩福利视频一区二区| 精品人妻一区二区三区麻豆| 欧美日韩亚洲国产一区二区在线观看 | 精品免费久久久久久久清纯 | 亚洲人成电影观看| av不卡在线播放| 国产一区二区三区综合在线观看| 久久人人爽av亚洲精品天堂| 女人久久www免费人成看片| 99九九在线精品视频| 精品免费久久久久久久清纯 | 69精品国产乱码久久久| 极品少妇高潮喷水抽搐| 国产男女内射视频| 深夜精品福利| 精品卡一卡二卡四卡免费| 成人亚洲精品一区在线观看| 亚洲自偷自拍图片 自拍| 精品福利永久在线观看| 久久久久久免费高清国产稀缺| 精品国产乱码久久久久久小说| 午夜免费成人在线视频| 国产精品免费大片| 亚洲,欧美,日韩| 一区二区三区四区激情视频| 香蕉丝袜av| a 毛片基地| 亚洲五月婷婷丁香| 老熟女久久久| 精品欧美一区二区三区在线| 午夜老司机福利片| 欧美 日韩 精品 国产| 亚洲av成人精品一二三区| 欧美av亚洲av综合av国产av| 亚洲欧洲精品一区二区精品久久久| 国产精品二区激情视频| 亚洲欧美一区二区三区久久| 伊人亚洲综合成人网| 日韩中文字幕欧美一区二区 | 男人添女人高潮全过程视频| 国产高清国产精品国产三级| 性高湖久久久久久久久免费观看| 在现免费观看毛片| 91国产中文字幕| 久久人人97超碰香蕉20202| 一本—道久久a久久精品蜜桃钙片| 亚洲精品一区蜜桃| 午夜老司机福利片| 亚洲图色成人| 视频区图区小说| 晚上一个人看的免费电影| 操出白浆在线播放| 天天影视国产精品| 成人18禁高潮啪啪吃奶动态图| 国产一区二区激情短视频 | 一个人免费看片子| 午夜久久久在线观看| 人妻一区二区av| 91精品伊人久久大香线蕉| 国产免费福利视频在线观看| 激情视频va一区二区三区| 在线观看免费午夜福利视频| 午夜影院在线不卡| 亚洲情色 制服丝袜| 久久女婷五月综合色啪小说| 成年人午夜在线观看视频| 久久99精品国语久久久| 精品亚洲乱码少妇综合久久| www.精华液| 捣出白浆h1v1| 欧美性长视频在线观看| 两性夫妻黄色片| 极品少妇高潮喷水抽搐| 伦理电影免费视频| 人人妻人人澡人人爽人人夜夜| 男人爽女人下面视频在线观看| 亚洲美女黄色视频免费看| 十八禁高潮呻吟视频| 欧美国产精品一级二级三级| 女人久久www免费人成看片| www日本在线高清视频| 中文欧美无线码| 韩国精品一区二区三区| √禁漫天堂资源中文www| 成年人免费黄色播放视频| 欧美 亚洲 国产 日韩一| 只有这里有精品99| 精品亚洲成a人片在线观看| 国产亚洲午夜精品一区二区久久| 老汉色av国产亚洲站长工具| 日本欧美国产在线视频| 午夜av观看不卡| tube8黄色片| 久久国产精品人妻蜜桃| 美女脱内裤让男人舔精品视频| 亚洲欧美日韩另类电影网站| 黄色片一级片一级黄色片| 国产淫语在线视频| 校园人妻丝袜中文字幕| 国产免费福利视频在线观看| 999精品在线视频| 首页视频小说图片口味搜索 | 丰满少妇做爰视频| 1024视频免费在线观看| 久久人人97超碰香蕉20202| 亚洲国产中文字幕在线视频| 国产亚洲欧美在线一区二区| 水蜜桃什么品种好| 国产极品粉嫩免费观看在线| 欧美黑人精品巨大| 精品高清国产在线一区| 美女大奶头黄色视频| 18禁观看日本| 国产亚洲午夜精品一区二区久久| 性高湖久久久久久久久免费观看| 国产精品一区二区精品视频观看| 久久精品久久久久久噜噜老黄| 一区福利在线观看| 男男h啪啪无遮挡| 天天躁夜夜躁狠狠久久av| 国产欧美日韩一区二区三区在线| 午夜久久久在线观看| 国产三级黄色录像| 一二三四在线观看免费中文在| 国产麻豆69| 国产欧美日韩一区二区三 | 免费高清在线观看视频在线观看| 男女午夜视频在线观看| 9色porny在线观看| 亚洲欧美精品自产自拍| 少妇粗大呻吟视频| 多毛熟女@视频| 女人久久www免费人成看片| 欧美 亚洲 国产 日韩一| 亚洲熟女精品中文字幕| 日韩视频在线欧美| 啦啦啦在线免费观看视频4| 国产精品久久久久久人妻精品电影 | 亚洲国产精品一区三区| 欧美日韩亚洲综合一区二区三区_| 建设人人有责人人尽责人人享有的| 亚洲免费av在线视频| 亚洲,欧美,日韩| www日本在线高清视频| 在线观看免费午夜福利视频| 考比视频在线观看| 一级毛片我不卡| 大香蕉久久成人网| 巨乳人妻的诱惑在线观看| 啦啦啦在线免费观看视频4| 电影成人av| 亚洲精品一区蜜桃| 亚洲情色 制服丝袜| 大话2 男鬼变身卡| 少妇 在线观看| 欧美人与性动交α欧美精品济南到| 亚洲国产精品国产精品| 色综合欧美亚洲国产小说| 久久99热这里只频精品6学生| 免费在线观看视频国产中文字幕亚洲 | 日韩精品免费视频一区二区三区| 国产精品一区二区免费欧美 | 亚洲激情五月婷婷啪啪| 亚洲精品在线美女| 日本wwww免费看| 一区福利在线观看| 校园人妻丝袜中文字幕| 亚洲av电影在线进入| 国产精品久久久久成人av| 久久久久久久久久久久大奶| 亚洲精品在线美女| 一区二区av电影网| 在线观看免费视频网站a站| 成人18禁高潮啪啪吃奶动态图| 亚洲av日韩精品久久久久久密 | 国产精品久久久人人做人人爽| 精品一区二区三区av网在线观看 | 热re99久久国产66热| 精品国产乱码久久久久久男人| 侵犯人妻中文字幕一二三四区| 首页视频小说图片口味搜索 | 色婷婷久久久亚洲欧美| 人人妻人人澡人人看| 亚洲国产精品999| 一边摸一边抽搐一进一出视频| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕人妻丝袜一区二区| 黄色视频在线播放观看不卡| 色婷婷av一区二区三区视频| 中文字幕人妻丝袜一区二区| 国产精品一区二区精品视频观看| 日韩大片免费观看网站| 两个人免费观看高清视频| 人人澡人人妻人| 欧美激情极品国产一区二区三区| 欧美国产精品va在线观看不卡| 亚洲欧美精品自产自拍| 一二三四在线观看免费中文在|