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

    基于復(fù)合雨量指標(biāo)的山洪預(yù)警方法及其應(yīng)用

    2017-11-03 09:13:16文川
    水利水電快報(bào) 2017年10期
    關(guān)鍵詞:山洪雨量洪水

    文川

    (華北水利水電大學(xué) 水利學(xué)院,河南 鄭州 450046)

    基于復(fù)合雨量指標(biāo)的山洪預(yù)警方法及其應(yīng)用

    王路王文川李浩田丹石佳勇

    (華北水利水電大學(xué)水利學(xué)院,河南鄭州450046)

    為滿足山洪災(zāi)害預(yù)警需要,提出了一種由前期影響雨量、本場降雨量和雨強(qiáng)組成的復(fù)合雨量指標(biāo)。選擇江西省彭坊小流域?yàn)檠芯繀^(qū),根據(jù)該流域的實(shí)測雨洪資料及預(yù)警流量,建立了出口斷面流量與復(fù)合雨量指標(biāo)之間的關(guān)系,確定了復(fù)合雨量指標(biāo)閾值。結(jié)果表明,流量與復(fù)合雨量指標(biāo)值關(guān)系密切, 復(fù)合雨量指標(biāo)值越大,越容易發(fā)生山洪災(zāi)害;通過連續(xù)計(jì)算復(fù)合雨量指標(biāo)并與閾值比較,可快速預(yù)測山洪災(zāi)害,效果較好。復(fù)合雨量指標(biāo)只需要雨量信息,計(jì)算簡便,可為研究區(qū)山洪預(yù)警提供參考。

    山洪預(yù)警;預(yù)警流量;復(fù)合雨量指標(biāo);前期影響雨量;閾值

    據(jù)統(tǒng)計(jì),我國1950~2000年洪澇災(zāi)害死亡人數(shù)為26.3萬人,其中山丘區(qū)的死亡人數(shù)占總死亡人數(shù)的68.4%[1]。2003~2013年,我國因山洪災(zāi)害死亡人數(shù)的比例上升到79.5%[2],山洪災(zāi)害已成為我國防洪減災(zāi)的薄弱環(huán)節(jié)和造成人員傷亡的突出災(zāi)種。準(zhǔn)確及時(shí)的山洪預(yù)報(bào)預(yù)警,可為處在山洪危險(xiǎn)區(qū)的人口和資產(chǎn)贏得寶貴的避險(xiǎn)時(shí)間,因此迫切需要加強(qiáng)突發(fā)性山洪災(zāi)害預(yù)報(bào)預(yù)警技術(shù)的研究。

    目前,山洪預(yù)警方法大致可概括為兩類:一是順推法,二是逆推法。順推法是由降雨預(yù)報(bào)山洪防護(hù)斷面的流量,一旦流量達(dá)到閾值,便發(fā)布預(yù)警信息。如:基于TOPKAPI模型的中小河流洪水預(yù)報(bào)系統(tǒng)[3];基于HEC-HMS模型的山洪預(yù)報(bào)系統(tǒng)[4];水利部水文局開發(fā)的中小河流山洪預(yù)警預(yù)報(bào)系統(tǒng)[5];河南省水利廳組織開發(fā)的山洪災(zāi)害預(yù)警預(yù)報(bào)系統(tǒng)[6]。此類方法的優(yōu)點(diǎn)是理論較完善,不足在于對資料缺乏地區(qū)模型參數(shù)率定難[7],對技術(shù)支撐平臺要求高。逆推法是由山洪防護(hù)斷面的預(yù)警流量反推出產(chǎn)生該流量的臨界雨量,一旦降雨量達(dá)到閾值,便發(fā)布山洪預(yù)警信息。如:美國水文研究中心研發(fā)的基于動(dòng)態(tài)臨界雨量的山洪預(yù)警系統(tǒng)FFGS[8,9];劉志雨等探討的以動(dòng)態(tài)臨界雨量為指標(biāo)的山洪預(yù)警方法[10,11];江錦紅等提出的雙指標(biāo)暴雨臨界曲線法[12];劉嬡媛等討論的資料匱乏地區(qū)山洪災(zāi)害臨界雨量確定方法[13];陳瑜彬等提出的基于馬氏距離識別的流域臨界雨量擬定方法[14]。因山洪源短流急、突發(fā)性強(qiáng),基于逆推思路的臨界雨量對預(yù)警更為實(shí)用。

    目前,臨界雨量大多采取事先設(shè)定若干種初始土壤含水量(如:流域蓄水容量的0.2,0.5或0.8等)逆推出固定歷時(shí)(如:1 h,2 h或3 h等)臨界值的形式。但這種形式的臨界雨量存在一定局限性:①南方濕潤地區(qū)的山洪成災(zāi)流量都是由連續(xù)多時(shí)段降雨在出口斷面產(chǎn)生的流量疊加形成的,很少純粹由1 h或2 h降雨引發(fā);②土壤含水量是一個(gè)連續(xù)變量,當(dāng)初始土壤含水量剛好不等于事先設(shè)定的有限種狀態(tài)值時(shí)則不便使用。

    為此,本文結(jié)合典型流域,從分析出口斷面流量與各影響因素的相關(guān)性入手,提出一種復(fù)合雨量指標(biāo)及其閾值確定方法,并檢驗(yàn)其應(yīng)用效果。

    1 研究區(qū)概況

    彭坊水文站位于東經(jīng)114°20′、北緯27°14′,屬于贛江水系,在江西省吉安市安??h境內(nèi)西南部。流域年均降水量1 808.7 mm,年均徑流深1 050.5 mm,平均徑流系數(shù)0.58;4~8月的降水量約占年降水量的72.2%,徑流深約占年徑流深的71.7%,屬于典型的濕潤區(qū)。流域絕大部分面積為森林覆蓋的山區(qū),僅在沿河道兩側(cè)分布有少量耕地和草地,大部分面積上的土壤質(zhì)地為砂粘土。流域蓄水容量為120 mm。流域內(nèi)有龍下、坪江頭、官田、炎里、老洲、由路、深坳、彭坊8個(gè)國家基本雨量站。彭坊水文站以上流域(簡稱“彭坊流域”)的地理位置及雨量站分布如圖1所示。

    圖1 彭坊流域地理位置及雨量站點(diǎn)分布

    2 研究思路與方法

    2.1 構(gòu)建復(fù)合雨量指標(biāo)的思路

    流域初始土壤含水量、本次降雨的雨量和雨強(qiáng)是影響流域出口斷面流量的主要因素。為避開對土壤含水量的依賴,本文以前期影響雨量代替初始土壤含水量,從而將影響流量的因素全部歸結(jié)為降雨量。從前期影響雨量、本次雨量和雨強(qiáng)三要素出發(fā)來構(gòu)建復(fù)合雨量指標(biāo),并建立出口斷面流量與復(fù)合雨量指標(biāo)的定量關(guān)系。

    2.2 確定復(fù)合雨量指標(biāo)及其閾值的步驟

    (1) 根據(jù)山洪防治斷面的大斷面資料確定上灘水位,由該斷面的水位流量關(guān)系確定山洪預(yù)警流量。

    (2) 采用地區(qū)經(jīng)驗(yàn)公式推求斷面以上小流域的瞬時(shí)單位線的參數(shù)n和K。根據(jù)流域瞬時(shí)單位線的概念,n和K的乘積就是單位線滯時(shí),相當(dāng)于平均流域匯流時(shí)間。最大流域匯流時(shí)間則根據(jù)經(jīng)驗(yàn)近似取為2。

    (3) 選取代表性洪水,分兩種情況判斷每場洪水需計(jì)算復(fù)合雨量指標(biāo)的時(shí)刻t′:一是洪峰未達(dá)到預(yù)警流量的場次洪水,選洪峰流量的出現(xiàn)時(shí)刻;二是洪峰超過預(yù)警流量的場次洪水,選首個(gè)超預(yù)警流量的出現(xiàn)時(shí)刻。采用泰森多邊形法計(jì)算流域面平均雨量。根據(jù)匯流時(shí)間,統(tǒng)計(jì)每場降雨過程中t′時(shí)刻之前0~2 h之間雨強(qiáng)、2 h~之間雨強(qiáng)、2 h~2之間雨量,以及t′-2時(shí)刻的前期影響雨量。

    分析出口斷面流量與上述各降雨統(tǒng)計(jì)量的單相關(guān)關(guān)系,從中選出對出口斷面流量有較大貢獻(xiàn)的時(shí)段雨強(qiáng)(記為PI)和時(shí)段雨量(記為Pc)。參考文獻(xiàn)[15]中“降雨驅(qū)動(dòng)指標(biāo)”的概念,按式(1)計(jì)算與t′時(shí)刻流量對應(yīng)的復(fù)合雨量指標(biāo):

    CRI=PI×P=PI×(Pa+Pc)

    (1)

    式中,CRI為復(fù)合雨量指標(biāo),mm2/h;PI代表本次降雨中對t′時(shí)刻流量起主要作用的時(shí)段雨強(qiáng),mm/h;Pa為前期影響雨量,mm;Pc代表本次降雨中對t′時(shí)刻流量起主要作用的時(shí)段雨量,mm。t′~2時(shí)刻的前期影響雨量由此前連續(xù)15 d的雨量的遞推計(jì)算得到。因研究區(qū)屬于濕潤地區(qū),取遞推Pa的日折減系數(shù)為0.85,并在計(jì)算過程中控制Pa≤Wm,mm;P為有效累積雨量,等于Pa與Pc之和,mm。各降雨要素計(jì)算起止時(shí)間如圖2所示。

    圖2 各降雨要素計(jì)算起止時(shí)間示意

    (4) 根據(jù)所選代表性暴雨洪水的復(fù)合雨量指標(biāo)及其相應(yīng)流量,建立出口斷面流量與復(fù)合雨量指標(biāo)的定量關(guān)系。再結(jié)合警戒流量,確定復(fù)合雨量指標(biāo)閾值。

    3 研究成果

    3.1 預(yù)警流量

    根據(jù)彭坊水文站大斷面資料及水位流量關(guān)系曲線,確定上灘水位為45.3 m,預(yù)警流量為181 m3/s。

    表1 典型洪水的洪峰及影響要素

    3.2 流域匯流時(shí)間

    基于DEM數(shù)據(jù)提取得到彭坊流域的集水面積為122 km2、主河道長度為34.1 km、主河道平均比降為0.007 6。參考《江西省暴雨洪水查算手冊》[16],彭坊水文站所在的贛江水系洲湖水屬于江西省第Ⅲ產(chǎn)匯流分區(qū)。根據(jù)流域面積,取瞬時(shí)單位線的參數(shù)n=2,并采用式(2)計(jì)算m1。

    m1=n·K=5.624(F/J)0.143(I/10)0.134 8Lg(F/J)-0.500

    (2)

    式中,F(xiàn)為流域面積,km2;J為主河道平均比降,‰;I為時(shí)段平均雨強(qiáng),mm/h。

    將流域面積F、主河道平均比降J、時(shí)段平均雨強(qiáng)I(取I=10 mm/h)代入式(2)得m1=8,則K=4。因此,彭坊流域的平均匯流時(shí)間為8 h,最大匯流時(shí)間近似取16 h。

    3.3 復(fù)合雨量指標(biāo)及其閾值

    彭坊水文站從1978年開始觀測流量,雨量站除了龍下站1980年設(shè)站、彭坊站1963年設(shè)站外,其余6站均為1978年設(shè)站。在1983~2013年間,選擇彭坊小流域15場暴雨洪水(其中6場達(dá)到警戒流量)時(shí)段長為0.5 h的配套水文資料,采用泰森多邊形法計(jì)算流域面平均雨量。確定每場洪水的t′時(shí)刻、計(jì)算t′時(shí)刻之前0~2 h之間雨強(qiáng)、2~8 h之間雨強(qiáng)、2~16 h之間雨量,以及t′提前16 h時(shí)刻的前期影響雨量Pa,結(jié)果見表1。

    根據(jù)表1樣本,分別進(jìn)行流量與前2 h雨強(qiáng)、前3~8 h雨強(qiáng)、前3~16 h雨量、前3~16 h雨量與16 h前Pa之和的單相關(guān)分析,以判別哪些因子對流量起主要作用,單相關(guān)分析結(jié)果見表2。

    表2 流量與其影響因子的單相關(guān)分析結(jié)果

    可見,相比于前3~8 h雨強(qiáng),前2 h雨強(qiáng)與流量的相關(guān)關(guān)系不太密切,這是由于流域平均匯流時(shí)間為8 h,前2 h降雨的產(chǎn)流量大部分還沒來得及到達(dá)出口斷面。因此,選前3~8 h雨強(qiáng)作為PI。按式(1)計(jì)算的復(fù)合雨量指標(biāo)結(jié)果見表3。

    表3 場次暴雨的PI、P及CPI指標(biāo)值

    各場暴雨洪水中,復(fù)合雨量指標(biāo)CRI與對應(yīng)流量的相關(guān)關(guān)系見圖3。由表3和圖2可見,達(dá)到警戒流量的場次暴雨的CRI值均大于未達(dá)到警戒流量的場次暴雨的CRI值,且流量與復(fù)合雨量指標(biāo)值之間相關(guān)關(guān)系密切,流量整體上隨著復(fù)合雨量指標(biāo)值的增大而增大,表明所構(gòu)建的復(fù)合雨量指標(biāo)能有效表征其所形成的流量。

    圖3 復(fù)合雨量指標(biāo)與相應(yīng)流量的關(guān)系

    將彭坊流域出口斷面警戒流量181 m3/s代入復(fù)合雨量指標(biāo)與對應(yīng)流量的回歸方程,得到復(fù)合雨量指標(biāo)閾值為2 046 mm2/h。但從表3中可見20100618洪水的復(fù)合雨量指標(biāo)為1 704 mm2/h,產(chǎn)生的流量為189 m3/s,20110613洪水的復(fù)合雨量指標(biāo)為1 981 mm2/h,產(chǎn)生的流量為201 m3/s。由于回歸方程是按照最小二乘法確定的,代表平均情況,而由于實(shí)際受各種不確定因素的影響,比2 046 mm2/h略小的復(fù)合雨量指標(biāo)值也可能產(chǎn)生超警戒流量。從偏安全角度考慮,最終推薦彭坊流域的山洪復(fù)合雨量指標(biāo)閾值為1 700 mm2/h。

    3.4 山洪預(yù)警實(shí)例

    以1995年6月29~30日洪水、2010年6月24日洪水為例,說明復(fù)合雨量指標(biāo)如何用于山洪預(yù)警,并檢驗(yàn)其應(yīng)用效果。

    圖4(a)是1995年6月29~30日洪水期間每隔0.5 h的雨量及復(fù)合雨量指標(biāo)變化過程;圖4(b)是1995年6月29~30日洪水期間每隔0.5 h的流量及復(fù)合雨量指標(biāo)變化過程。

    圖4 復(fù)合雨量指標(biāo)在1995年6月29~30日山洪預(yù)警中的應(yīng)用

    圖5 復(fù)合雨量指標(biāo)在2010年6月24日山洪預(yù)警中的應(yīng)用

    按照本文提出的方法,根據(jù)截至1995年6月30日1:30的實(shí)測雨量資料,計(jì)算出6月30日3:30的復(fù)合雨量為1 417 mm2/h,根據(jù)截至6月30日2:00的實(shí)測雨量資料,計(jì)算出6月30日4:00的復(fù)合雨量為1 864 mm2/h,根據(jù)山洪復(fù)合雨量指標(biāo)閾值1 700 mm2/h,判斷4:00將發(fā)生山洪災(zāi)害,從而做出預(yù)警。實(shí)測流量表明6月30日4:00的流量為113 m3/s,5:00流量就漲到236 m3/s。而且從圖4(b)可見,流量與復(fù)合雨量指標(biāo)的走勢吻合度高。

    圖5(a)是2010年6月24日洪水期間每隔0.5 h的雨量及復(fù)合雨量指標(biāo)變化過程;圖5(b)是2010年6月24日洪水期間每隔0.5 h的流量及復(fù)合雨量指標(biāo)變化過程。

    同樣,根據(jù)截至2010年6月24日11: 00掌握的雨量資料,計(jì)算出6月24日13:00的復(fù)合雨量為1 653 mm2/h,根據(jù)截至6月24日11:30掌握的雨量資料,計(jì)算出6月24日13:30 的復(fù)合雨量為1 721 mm2/h,根據(jù)山洪復(fù)合雨量指標(biāo)閾值為1 700 mm2/h,判斷13:00將發(fā)生山洪災(zāi)害,從而做出預(yù)警。實(shí)測流量表明6月24日13:30的流量為160 m3/s,14:30流量就漲到180 m3/s,達(dá)到警戒流量。

    可見,復(fù)合雨量指標(biāo)對兩場洪水的預(yù)警是有效的。兩場洪水中預(yù)報(bào)的山洪發(fā)生時(shí)間均比實(shí)際發(fā)生時(shí)間提前了1 h,這與選取了略偏安全的指標(biāo)閾值有關(guān)。

    4 結(jié) 論

    本文研究了復(fù)合雨量指標(biāo)構(gòu)建及其在彭坊小流域的應(yīng)用,主要結(jié)論如下:

    (1) 山洪防護(hù)斷面的流量與相應(yīng)的復(fù)合雨量指標(biāo)之間存在密切的相關(guān)關(guān)系,復(fù)合雨量指標(biāo)值越大,在流域出口斷面形成的流量就越大,越容易發(fā)生山洪災(zāi)害。

    (2) 經(jīng)綜合分析,從偏安全角度確定彭坊流域的復(fù)合雨量指標(biāo)閾值為1 700 mm2/h。由于流域內(nèi)最近2 h的降雨大部分還未及到達(dá)出口斷面,可利用實(shí)時(shí)雨量觀測資料計(jì)算2 h后的復(fù)合雨量指標(biāo),與指標(biāo)閾值作對比,就能快速判斷2 h后是否會(huì)發(fā)生山洪災(zāi)害。該方法在兩場山洪中的應(yīng)用表明預(yù)警效果較好。

    (3) 復(fù)合雨量指標(biāo)計(jì)算只需要雨量資料,不需要土壤含水量信息,不受事先設(shè)定的若干種初始土壤含水量的限制,簡便易行,可為研究區(qū)山洪預(yù)警提供依據(jù)。

    [1] 劉志雨.山洪預(yù)警預(yù)報(bào)技術(shù)研究與應(yīng)用[J].中國防汛抗旱,2012,22(2):41-45.

    [2] 國家防汛抗旱總指揮部,中華人民共和國水利水利部.中國水旱災(zāi)害公報(bào)2013 [R].北京:國家防汛抗旱總指揮部,中華人民共和國水利水利部,2014.

    [3] 劉志雨.基于GIS 的分布式托普卡匹水文模型在洪水預(yù)報(bào)中的應(yīng)用[J].水利學(xué)報(bào),2004(5):70-75.

    [4] 徐宗學(xué),程磊.分布式水文模型研究與應(yīng)用進(jìn)展[J].水利學(xué)報(bào),2010,41(9):1009-10 17.

    [5] 郭良,唐學(xué)哲,孔凡哲.基于分布式水文模型的山洪災(zāi)害預(yù)警預(yù)報(bào)系統(tǒng)研究及其應(yīng)用[J].中國水利, 2007(14) : 38-41.

    [6] 胡健偉,劉志雨. 中小河流山洪預(yù)警預(yù)報(bào)系統(tǒng)開發(fā)設(shè)計(jì)及應(yīng)用[J].水文, 2011, 31(3) : 18-21.

    [7] 葉金印,李致家,吳勇拓. 一種用于缺資料地區(qū)山洪預(yù)警方法研究與應(yīng)用[J].水力發(fā)電學(xué)報(bào), 2013, 32(3): 15-19.

    [8] Georgakakos KP. Analytical results for operational flash flood guidance[J].Journal of Hydrology,2006,317: 81-103.

    [9] Hapuarachchi H,Wang Q J,Pagano TC.A review of advances inflash flood forecasting[J].Hydrological processes,2011,25(18):2771-2784.

    [10]劉志雨,楊大文,胡健偉. 基于動(dòng)態(tài)臨界雨量的中小河流山洪預(yù)警方法及其應(yīng)用[J].北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,46(3): 317-321.

    [11]梁家志,劉志雨.中小河流山洪監(jiān)測與預(yù)警預(yù)測技術(shù)研究[M].北京:科學(xué)出版社,2010.

    [12]江錦紅,邵利萍.基于降雨觀測資料的山洪預(yù)警標(biāo)準(zhǔn)[J].水利學(xué)報(bào),2014,41(4):458-463.

    [13]劉嬡媛,胡昌偉,張紅萍,等.資料匱乏地區(qū)山洪災(zāi)害臨界雨量確定方法分析[J].水利水電技術(shù),2014,45(8):15-17.

    [14]陳瑜彬,楊文發(fā),陳力,等.一種基于馬氏距離識別法的流域臨界雨量擬定方法應(yīng)用探討[J].水利水電技術(shù),2016,47(1):21-25.

    [15]李明熹.臺灣土石流監(jiān)測及警戒系統(tǒng)之綜述[J].水土保持研究.2009,16(6):239-242.

    [16]江西省水文局.江西省暴雨洪水查算手冊[R].南昌:江西省水文局,2010.

    2017-05-15

    國家自然科學(xué)基金項(xiàng)目(51509088);華北水利水電大學(xué)大學(xué)生創(chuàng)新訓(xùn)練項(xiàng)目(2016XB079)

    王路,女,華北水利水電大學(xué)水利學(xué)院.

    王文川,男,華北水利水電大學(xué)水利學(xué)院,博士, 教授,博士生導(dǎo)師.

    1006-0081(2017)10-0032-05

    TV122

    A

    (編輯李慧)

    猜你喜歡
    山洪雨量洪水
    寧夏紅柳溝流域水沙變化及產(chǎn)沙分析
    優(yōu)雅地表達(dá)
    洪水時(shí)遇到電線低垂或折斷該怎么辦
    基于小波去噪的稱重雨量數(shù)據(jù)分析
    又見洪水(外二首)
    天津詩人(2017年2期)2017-11-29 01:24:25
    SL—1 型雨量傳感器故障分析排除和維護(hù)
    西藏科技(2016年5期)2016-09-26 12:16:40
    遭遇暴雨山洪如何避險(xiǎn)自救
    洪水來了
    幼兒畫刊(2016年8期)2016-02-28 21:00:52
    湖北省山洪溝治理思路淺析
    中國水利(2015年9期)2015-02-28 15:13:20
    “慧圖杯”山洪災(zāi)害防治知識競賽
    中國水利(2015年1期)2015-02-28 15:11:55
    久久狼人影院| 99riav亚洲国产免费| 真人做人爱边吃奶动态| 国产av不卡久久| 亚洲第一青青草原| 日韩成人在线观看一区二区三区| 亚洲 国产 在线| 国产爱豆传媒在线观看 | 国产精品自产拍在线观看55亚洲| 中文字幕人妻丝袜一区二区| 久久婷婷人人爽人人干人人爱| www日本黄色视频网| 久久久久久九九精品二区国产 | 日韩欧美在线二视频| 1024视频免费在线观看| 久久久水蜜桃国产精品网| 日韩欧美一区二区三区在线观看| 国产精品久久久久久亚洲av鲁大| 欧美成人午夜精品| 亚洲国产日韩欧美精品在线观看 | 成人18禁在线播放| 国产日本99.免费观看| 一进一出抽搐动态| 母亲3免费完整高清在线观看| 90打野战视频偷拍视频| 久久天堂一区二区三区四区| 美女午夜性视频免费| 日本五十路高清| 亚洲国产欧美日韩在线播放| 免费无遮挡裸体视频| 午夜日韩欧美国产| xxx96com| 正在播放国产对白刺激| av片东京热男人的天堂| 好男人在线观看高清免费视频 | 成人欧美大片| 国产高清videossex| 精品欧美一区二区三区在线| 国产av不卡久久| 一区福利在线观看| 99久久久亚洲精品蜜臀av| 热re99久久国产66热| 丰满人妻熟妇乱又伦精品不卡| 国产日本99.免费观看| 最近最新免费中文字幕在线| 亚洲,欧美精品.| 国产一区在线观看成人免费| 丝袜在线中文字幕| 婷婷亚洲欧美| 色老头精品视频在线观看| 中文字幕最新亚洲高清| 欧美日韩亚洲国产一区二区在线观看| 国产日本99.免费观看| 国产视频内射| 中文字幕精品亚洲无线码一区 | 国产国语露脸激情在线看| av在线天堂中文字幕| 亚洲自拍偷在线| 午夜激情av网站| 侵犯人妻中文字幕一二三四区| 人人妻人人看人人澡| 成人三级做爰电影| 久久九九热精品免费| 叶爱在线成人免费视频播放| 国产伦人伦偷精品视频| 欧美在线一区亚洲| 在线观看免费日韩欧美大片| 免费在线观看黄色视频的| 亚洲人成伊人成综合网2020| 少妇的丰满在线观看| 99热这里只有精品一区 | 国产亚洲av高清不卡| 精品不卡国产一区二区三区| 精品国产国语对白av| 最新美女视频免费是黄的| 人妻久久中文字幕网| 亚洲第一欧美日韩一区二区三区| 免费在线观看黄色视频的| 在线看三级毛片| 在线观看免费午夜福利视频| 久久久久久久精品吃奶| 在线免费观看的www视频| 99国产极品粉嫩在线观看| av电影中文网址| 在线观看舔阴道视频| 国产aⅴ精品一区二区三区波| 久久久久久久久中文| 色婷婷久久久亚洲欧美| 在线观看免费视频日本深夜| 国产精品综合久久久久久久免费| 久久久久久九九精品二区国产 | 亚洲精品一区av在线观看| 国产成年人精品一区二区| 哪里可以看免费的av片| 欧美性猛交黑人性爽| 亚洲五月婷婷丁香| 午夜精品久久久久久毛片777| 日韩有码中文字幕| 首页视频小说图片口味搜索| 制服人妻中文乱码| 好看av亚洲va欧美ⅴa在| 黑人巨大精品欧美一区二区mp4| 国产爱豆传媒在线观看 | 最好的美女福利视频网| 久久中文看片网| 亚洲av熟女| 久久精品91无色码中文字幕| av福利片在线| 91麻豆av在线| 男女那种视频在线观看| 欧美日本视频| 精品福利观看| 日本免费一区二区三区高清不卡| 午夜福利18| 宅男免费午夜| 免费看美女性在线毛片视频| 午夜福利在线观看吧| 两性夫妻黄色片| 国产亚洲av嫩草精品影院| 禁无遮挡网站| 青草久久国产| 国产亚洲av嫩草精品影院| 香蕉国产在线看| 国产麻豆成人av免费视频| 国产熟女午夜一区二区三区| 午夜福利视频1000在线观看| 久9热在线精品视频| 男人操女人黄网站| 给我免费播放毛片高清在线观看| 麻豆av在线久日| av超薄肉色丝袜交足视频| 啪啪无遮挡十八禁网站| 大型av网站在线播放| 精品免费久久久久久久清纯| a级毛片在线看网站| 亚洲精品久久成人aⅴ小说| 久久精品国产亚洲av香蕉五月| 免费电影在线观看免费观看| 亚洲精品色激情综合| 男女床上黄色一级片免费看| 99在线人妻在线中文字幕| 91成年电影在线观看| 男女之事视频高清在线观看| 国产黄色小视频在线观看| 国产三级黄色录像| 国产一区二区三区在线臀色熟女| or卡值多少钱| 欧美日韩亚洲综合一区二区三区_| 久久午夜亚洲精品久久| 成人亚洲精品av一区二区| 99久久综合精品五月天人人| 夜夜夜夜夜久久久久| 国产成人一区二区三区免费视频网站| 免费在线观看黄色视频的| 亚洲人成电影免费在线| 十八禁人妻一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| av片东京热男人的天堂| 国内精品久久久久精免费| 可以在线观看毛片的网站| 黄色毛片三级朝国网站| 看黄色毛片网站| 亚洲午夜理论影院| 成人18禁高潮啪啪吃奶动态图| 久久国产精品人妻蜜桃| 午夜免费激情av| 成人亚洲精品一区在线观看| 波多野结衣av一区二区av| 精品不卡国产一区二区三区| 免费在线观看亚洲国产| 长腿黑丝高跟| 少妇粗大呻吟视频| 国产精品乱码一区二三区的特点| 天堂影院成人在线观看| 国产av又大| 国产一区二区三区在线臀色熟女| 亚洲成人国产一区在线观看| 麻豆久久精品国产亚洲av| 天天一区二区日本电影三级| √禁漫天堂资源中文www| cao死你这个sao货| 91成年电影在线观看| 久久婷婷人人爽人人干人人爱| 18禁观看日本| 观看免费一级毛片| 精品一区二区三区四区五区乱码| av超薄肉色丝袜交足视频| 日韩av在线大香蕉| 长腿黑丝高跟| 丰满的人妻完整版| 黄色 视频免费看| 757午夜福利合集在线观看| 精品福利观看| 91老司机精品| 欧美大码av| 精品久久久久久久久久免费视频| 人人妻人人澡人人看| 老司机福利观看| 欧美绝顶高潮抽搐喷水| 妹子高潮喷水视频| 村上凉子中文字幕在线| 午夜福利欧美成人| 日韩欧美一区二区三区在线观看| 女性生殖器流出的白浆| 一级作爱视频免费观看| 黄色a级毛片大全视频| 精品少妇一区二区三区视频日本电影| 一级a爱片免费观看的视频| 久久久久久大精品| 日韩有码中文字幕| 精品国产乱子伦一区二区三区| 色播在线永久视频| 欧美亚洲日本最大视频资源| 久久精品aⅴ一区二区三区四区| 两个人视频免费观看高清| e午夜精品久久久久久久| 亚洲精品av麻豆狂野| 999久久久精品免费观看国产| 久久精品国产亚洲av香蕉五月| 搡老妇女老女人老熟妇| 日本一区二区免费在线视频| √禁漫天堂资源中文www| 欧美性长视频在线观看| 精品一区二区三区视频在线观看免费| 男女午夜视频在线观看| 香蕉av资源在线| 亚洲精品国产精品久久久不卡| 免费看美女性在线毛片视频| 嫁个100分男人电影在线观看| 在线观看日韩欧美| 国产精品亚洲av一区麻豆| www.熟女人妻精品国产| 女性被躁到高潮视频| 夜夜看夜夜爽夜夜摸| 自线自在国产av| 成人三级黄色视频| 国产精品自产拍在线观看55亚洲| 亚洲国产精品成人综合色| 中文字幕人妻熟女乱码| 亚洲国产精品久久男人天堂| 国产亚洲精品av在线| 午夜视频精品福利| 在线免费观看的www视频| 18禁裸乳无遮挡免费网站照片 | 色综合欧美亚洲国产小说| 看免费av毛片| 一级a爱视频在线免费观看| 999久久久国产精品视频| 亚洲国产毛片av蜜桃av| 成在线人永久免费视频| 免费在线观看黄色视频的| 久久国产亚洲av麻豆专区| 丝袜在线中文字幕| 在线天堂中文资源库| 免费人成视频x8x8入口观看| 亚洲一卡2卡3卡4卡5卡精品中文| 成人国语在线视频| 国产在线精品亚洲第一网站| 亚洲av片天天在线观看| 国产又黄又爽又无遮挡在线| 国产精品野战在线观看| 熟女电影av网| 丰满的人妻完整版| 亚洲国产欧美一区二区综合| 欧美色欧美亚洲另类二区| 人人澡人人妻人| 亚洲av五月六月丁香网| 天天躁狠狠躁夜夜躁狠狠躁| 国产欧美日韩精品亚洲av| 黄频高清免费视频| 日韩一卡2卡3卡4卡2021年| 成人18禁在线播放| 国产熟女xx| 麻豆久久精品国产亚洲av| 亚洲国产精品sss在线观看| 啦啦啦 在线观看视频| 麻豆国产av国片精品| 中文字幕最新亚洲高清| 国产精品99久久99久久久不卡| 国内揄拍国产精品人妻在线 | 国产精品日韩av在线免费观看| 国产成人精品无人区| 最新在线观看一区二区三区| 国产成+人综合+亚洲专区| 人人妻人人澡欧美一区二区| 一夜夜www| 久久久久久亚洲精品国产蜜桃av| 国产成+人综合+亚洲专区| 亚洲第一av免费看| 亚洲国产精品成人综合色| 亚洲国产精品久久男人天堂| 在线观看免费视频日本深夜| 88av欧美| 精品欧美国产一区二区三| 窝窝影院91人妻| 亚洲中文字幕日韩| 波多野结衣高清作品| 99精品久久久久人妻精品| 成人精品一区二区免费| 婷婷精品国产亚洲av在线| 日韩视频一区二区在线观看| 女人被狂操c到高潮| 欧美日韩瑟瑟在线播放| 久热这里只有精品99| 日韩大码丰满熟妇| 怎么达到女性高潮| 男人舔女人的私密视频| 成在线人永久免费视频| 特大巨黑吊av在线直播 | 亚洲国产精品sss在线观看| 婷婷亚洲欧美| 国产欧美日韩一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 香蕉国产在线看| 丰满人妻熟妇乱又伦精品不卡| 手机成人av网站| 国产男靠女视频免费网站| 国产av不卡久久| 给我免费播放毛片高清在线观看| 亚洲av电影不卡..在线观看| 国产视频一区二区在线看| 国产不卡一卡二| 成人免费观看视频高清| 欧美乱妇无乱码| 制服人妻中文乱码| 色哟哟哟哟哟哟| 深夜精品福利| 中国美女看黄片| 亚洲男人的天堂狠狠| 国产精品亚洲美女久久久| 国内揄拍国产精品人妻在线 | 日韩大码丰满熟妇| 国产成人啪精品午夜网站| 亚洲中文字幕一区二区三区有码在线看 | 免费av毛片视频| 欧美成人性av电影在线观看| 正在播放国产对白刺激| 日韩精品中文字幕看吧| 女人高潮潮喷娇喘18禁视频| 亚洲精品粉嫩美女一区| 制服诱惑二区| 国产99久久九九免费精品| 夜夜看夜夜爽夜夜摸| 国产av一区二区精品久久| 日本a在线网址| 国产激情久久老熟女| 亚洲中文字幕一区二区三区有码在线看 | 老鸭窝网址在线观看| 成人三级做爰电影| 欧美日韩精品网址| 午夜福利欧美成人| 亚洲七黄色美女视频| 可以免费在线观看a视频的电影网站| 99国产精品一区二区三区| 国产在线观看jvid| 亚洲五月天丁香| 怎么达到女性高潮| 成人国产一区最新在线观看| 国产精品av久久久久免费| 99久久国产精品久久久| 黑人操中国人逼视频| 搡老岳熟女国产| a级毛片在线看网站| 久久午夜综合久久蜜桃| 韩国精品一区二区三区| 亚洲九九香蕉| 精品欧美一区二区三区在线| 色婷婷久久久亚洲欧美| 国产久久久一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 久久狼人影院| 麻豆av在线久日| 看黄色毛片网站| 国产精品影院久久| 色精品久久人妻99蜜桃| 99国产精品99久久久久| 91成人精品电影| 国产成人影院久久av| 最近最新中文字幕大全电影3 | 丁香欧美五月| 亚洲av熟女| 制服人妻中文乱码| 日韩欧美国产一区二区入口| 欧美日韩亚洲综合一区二区三区_| 夜夜爽天天搞| 精品久久蜜臀av无| 免费在线观看日本一区| 中文亚洲av片在线观看爽| 亚洲熟妇熟女久久| 中文字幕精品亚洲无线码一区 | 此物有八面人人有两片| 欧美日韩亚洲综合一区二区三区_| 欧美黄色淫秽网站| 精品国产超薄肉色丝袜足j| 亚洲精品美女久久av网站| 国产麻豆成人av免费视频| 国产熟女xx| 亚洲精品色激情综合| 中出人妻视频一区二区| 久久婷婷成人综合色麻豆| 老鸭窝网址在线观看| 日日爽夜夜爽网站| 国产精品香港三级国产av潘金莲| 午夜福利在线观看吧| 国产又黄又爽又无遮挡在线| 97碰自拍视频| 中亚洲国语对白在线视频| 黑人操中国人逼视频| 亚洲第一青青草原| 免费在线观看黄色视频的| 中文字幕高清在线视频| 丰满的人妻完整版| 免费在线观看亚洲国产| 欧美日本视频| 日韩欧美一区视频在线观看| 久久香蕉激情| 色婷婷久久久亚洲欧美| 日本精品一区二区三区蜜桃| 国产国语露脸激情在线看| 欧美精品啪啪一区二区三区| 18禁国产床啪视频网站| 一进一出抽搐动态| 一本综合久久免费| 久久精品国产清高在天天线| 久久香蕉激情| 午夜激情av网站| 99国产综合亚洲精品| 啦啦啦观看免费观看视频高清| 美女免费视频网站| 国产色视频综合| 可以在线观看毛片的网站| 亚洲国产精品sss在线观看| 国产欧美日韩一区二区精品| 久久精品亚洲精品国产色婷小说| 免费电影在线观看免费观看| 国产成人一区二区三区免费视频网站| 又黄又爽又免费观看的视频| 久久久久久久久久黄片| 免费在线观看黄色视频的| 日韩国内少妇激情av| 亚洲精品美女久久久久99蜜臀| 狂野欧美激情性xxxx| 麻豆久久精品国产亚洲av| 国产伦人伦偷精品视频| 亚洲人成网站高清观看| 国内揄拍国产精品人妻在线 | 99在线视频只有这里精品首页| 中文字幕精品免费在线观看视频| 国产午夜精品久久久久久| 亚洲色图 男人天堂 中文字幕| 99国产精品一区二区三区| 国产一区二区在线av高清观看| 欧美一级毛片孕妇| 午夜福利高清视频| 免费看十八禁软件| 欧美绝顶高潮抽搐喷水| 国产黄a三级三级三级人| 在线播放国产精品三级| 在线视频色国产色| 国产熟女午夜一区二区三区| 国产精品98久久久久久宅男小说| 国产精品精品国产色婷婷| 啦啦啦韩国在线观看视频| 他把我摸到了高潮在线观看| 国产v大片淫在线免费观看| 亚洲av成人不卡在线观看播放网| 99riav亚洲国产免费| 曰老女人黄片| 免费观看人在逋| 成人一区二区视频在线观看| 午夜a级毛片| 午夜老司机福利片| 国语自产精品视频在线第100页| 一个人免费在线观看的高清视频| xxx96com| avwww免费| 亚洲成a人片在线一区二区| 国产精品久久久久久亚洲av鲁大| 亚洲五月婷婷丁香| 一二三四在线观看免费中文在| 国产黄色小视频在线观看| 精品第一国产精品| 国产成人系列免费观看| 欧美人与性动交α欧美精品济南到| 少妇熟女aⅴ在线视频| 俄罗斯特黄特色一大片| 看黄色毛片网站| 久久久水蜜桃国产精品网| 男女那种视频在线观看| 日韩欧美一区视频在线观看| 少妇的丰满在线观看| 男人舔奶头视频| 欧美色视频一区免费| 久久香蕉精品热| 国产精品 国内视频| 天天一区二区日本电影三级| 黄色视频,在线免费观看| 波多野结衣巨乳人妻| 亚洲国产欧洲综合997久久, | 精品久久久久久,| 国产麻豆成人av免费视频| 色哟哟哟哟哟哟| 男人的好看免费观看在线视频 | 一级a爱视频在线免费观看| 日本a在线网址| av电影中文网址| 一个人免费在线观看的高清视频| 两人在一起打扑克的视频| 日韩欧美 国产精品| 午夜精品在线福利| 久久香蕉精品热| 又大又爽又粗| av免费在线观看网站| 69av精品久久久久久| 天天一区二区日本电影三级| 少妇裸体淫交视频免费看高清 | 国产区一区二久久| 精品电影一区二区在线| 别揉我奶头~嗯~啊~动态视频| 丝袜在线中文字幕| 国产色视频综合| 精品国产超薄肉色丝袜足j| 久久性视频一级片| 曰老女人黄片| 大型av网站在线播放| 欧美乱色亚洲激情| 天天添夜夜摸| 看片在线看免费视频| 男女床上黄色一级片免费看| 日韩欧美在线二视频| 美女免费视频网站| 日本撒尿小便嘘嘘汇集6| 午夜福利在线在线| 欧美人与性动交α欧美精品济南到| 亚洲五月色婷婷综合| 一个人观看的视频www高清免费观看 | 人妻丰满熟妇av一区二区三区| 日韩欧美一区视频在线观看| 十分钟在线观看高清视频www| 男人舔女人的私密视频| 青草久久国产| 一边摸一边做爽爽视频免费| 亚洲国产精品成人综合色| 日韩精品中文字幕看吧| 香蕉丝袜av| 热re99久久国产66热| 国产av不卡久久| 午夜久久久久精精品| 婷婷六月久久综合丁香| 熟女少妇亚洲综合色aaa.| 一区福利在线观看| 亚洲精品国产一区二区精华液| 天天躁狠狠躁夜夜躁狠狠躁| 国产日本99.免费观看| 天天躁夜夜躁狠狠躁躁| 一本一本综合久久| 亚洲人成网站在线播放欧美日韩| 亚洲国产精品成人综合色| 岛国在线观看网站| 性色av乱码一区二区三区2| 久久久久免费精品人妻一区二区 | 18美女黄网站色大片免费观看| www日本在线高清视频| 黑丝袜美女国产一区| 好看av亚洲va欧美ⅴa在| 午夜福利免费观看在线| 国产99久久九九免费精品| 亚洲,欧美精品.| 正在播放国产对白刺激| 久热这里只有精品99| 欧美另类亚洲清纯唯美| 欧美黑人精品巨大| 亚洲av电影不卡..在线观看| 亚洲欧美日韩无卡精品| 最新在线观看一区二区三区| 波多野结衣av一区二区av| 成人亚洲精品av一区二区| 午夜福利一区二区在线看| 制服人妻中文乱码| 精品国产美女av久久久久小说| 亚洲男人的天堂狠狠| 精品国产一区二区三区四区第35| 色综合欧美亚洲国产小说| 欧美在线黄色| 精品欧美一区二区三区在线| 亚洲欧美日韩无卡精品| 成人国语在线视频| 91成人精品电影| 最新在线观看一区二区三区| 日韩成人在线观看一区二区三区| e午夜精品久久久久久久| 午夜免费成人在线视频| 在线看三级毛片| av天堂在线播放| 在线看三级毛片| 巨乳人妻的诱惑在线观看| 狂野欧美激情性xxxx| 99re在线观看精品视频| 校园春色视频在线观看| avwww免费| 日韩欧美在线二视频| 国产av一区二区精品久久| 高清在线国产一区| 午夜免费观看网址| 观看免费一级毛片| 久久香蕉国产精品| netflix在线观看网站| 欧美性猛交╳xxx乱大交人| 中文资源天堂在线| 国产伦人伦偷精品视频| 少妇熟女aⅴ在线视频| 69av精品久久久久久| 亚洲色图av天堂| 国产精华一区二区三区| 老司机靠b影院|