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

    采用流管元測量冰期流量及流速分布的方法

    2021-01-28 05:39:36潘佳佳郭新蕾李甲振郭永鑫
    水利學(xué)報 2020年12期
    關(guān)鍵詞:冰封冰蓋測流

    潘佳佳,郭新蕾,王 濤,付 輝,李甲振,郭永鑫

    (流域水循環(huán)模擬與調(diào)控國家重點實驗室,中國水利水電科學(xué)研究院,北京 100038)

    1 研究背景

    北半球較高緯度的地區(qū)約有60%以上河流在冬季會遭遇河冰問題,我國北緯30度以北占國土3/4的流域均可能出現(xiàn)流凌和冰封。不僅天然河流,北方人工大型輸水明渠冬季也常常會形成冰蓋。冬季冰期流量是寒區(qū)河渠水文、水力學(xué)測量和計算的重要水力參數(shù)之一,更是凌汛洪水預(yù)測預(yù)報的主要因子。它對于估算河渠輸水輸冰能力、冰塞冰壩預(yù)報、河道及其水利工程的安全設(shè)計起著關(guān)鍵作用。

    北方河流傳統(tǒng)冰期測流需要沿河道橫斷面在冰層上鉆孔10個以上,通過各孔測量垂線上不同水深處的流速。然后利用斷面上的流速積分得到相應(yīng)的總流量。但是,這種測量方法不能測量初封期或冰蓋不穩(wěn)定斷面的流量。而且寒冷惡劣氣候條件下耗費的人力、財力大,測量時間長,不利于大面積長歷時的流量測量[1-4]。長時間的冰蓋測流還置測量人員于落水高風(fēng)險中,時有冰蓋破裂人員墜水的事故[5-6]。這嚴重制約了冰期流量采集,影響北方河流冬季防凌減災(zāi)[7-8]。

    受季節(jié)性氣溫變化影響,高緯度寒區(qū)河流冬季流量觀測均需考慮冰蓋的影響[9-11]。Walker[12]針對美國愛荷華州的三條河流在1987—1988年冬季開展了6周的現(xiàn)場測量,分別采用冰期和非冰期的流量比例法、同流域的水文類比法、考慮河冰糙率的修正水位流量曲線法、單點指示流速法和均勻流法觀測了冰蓋下的斷面流量,并對不同方法的準確性進行了對比分析。結(jié)果表明不同測量方法均具有較大的偏差和不確定性,更為準確的流量測量還是依靠斷面流速分布測量再積分。美國地質(zhì)勘查局和加拿大水務(wù)局在1988—1991開展了冰期流量聯(lián)合觀測研究,在美國共測量了115個冰蓋下的流量及2300個垂線流速分布,在加拿大測量了83個流量和約1800個垂線流速分布[13]。研究表明0.5倍水深處的流速與水深平均流速存在系統(tǒng)偏差,修正系數(shù)約為0.88,而0.2與0.8倍水深處流速的平均值能很好反映水深平均流速[14]。我國《河流流量測驗規(guī)范》(GB50179-2015)詳細介紹了明流條件下河渠流量的測驗方法,如流速儀法、浮標法、比降面積法、聲學(xué)多普勒法、電波流速儀法、水工建筑物法等[15]。但在冬季結(jié)冰期,河流表面被冰雪覆蓋,這些明流測驗方法如浮標法、比降面積法、水工建筑物法等均無法開展。因此,該國標也推薦在冰蓋上鉆孔,再通過流速儀法或聲學(xué)多普勒法測量垂線流速分布,其中過水?dāng)嗝婷娣e應(yīng)扣除水浸冰和冰花面積[15]。如按以上方法測量冰層下的流量,需要測量沿水面布線方向所有垂線上的流速分布,進而利用各點垂線平均流速(十一點法、五點法、三點法等)和各部分小斷面面積,積分估算得到整個冰蓋斷面下的總流量。寒區(qū)天然河流水面寬度較大,布設(shè)垂向測點多時,如垂向測線大于50條,則測線位置流速分布的測量工作量很大。

    為了提高冰期直接測流的效率,大量研究推導(dǎo)了理論公式以簡化冰蓋下的流量觀測?;趷垡蛩固沟谋诿婧痛裁娣謪^(qū)理論[16],Calkins等劃分了河道斷面冰蓋影響區(qū)及河床影響區(qū)[17]。在斷面平均流速和水力坡度恒定的條件下,他們得到河道綜合糙率系數(shù),考慮了床面和冰面的綜合影響。該方法能有效估算冰蓋影響下的斷面流量,但需要給定冰蓋及河床糙率[18-20]?;诰C合糙率系數(shù),Shen和Ackermann 在忽略橫向環(huán)流和界面水流內(nèi)摩擦力的條件下,推導(dǎo)了冰蓋下單寬流量橫向分布的公式。該公式能較好地估算部分冰封和完全冰封河流的流量,模擬結(jié)果與日本和美國河流實測資料相符[21]。受暢流期Sharifi等提出的SKM方法啟發(fā)[22-23],楊開林從冰蓋下的雷諾平均圣維南方程出發(fā),推導(dǎo)了冰蓋下的水深平均流速分布,可基于已有河道地形計算冰下流量[20]。該方法隨后被槐文信等推廣到部分冰蓋和完全冰蓋下的流量、水深平均流速及水位流量關(guān)系分析,在率定阻力系數(shù)的條件下具有較好準確性[24-26]。這些研究有效推動了冰下測流的理論分析,但缺少與現(xiàn)場觀測技術(shù)的結(jié)合。

    針對以上研究的不足和冰蓋下流量觀測的難題,本文提出將流管元法和愛因斯坦的斷面分區(qū)方法相結(jié)合的冰期測流方法。該方法能準確計算實驗室復(fù)式斷面冰蓋下的流速分布,也能模擬部分冰封和完全冰封條件下黃河典型斷面的流量及單寬流量分布。新方法能有效減少觀測工作量,提高冬季測流效率,能廣泛用于冰期和非冰期的流速及流量觀測。

    2 理論原理

    2.1 理論方程基于愛因斯坦的斷面分區(qū)理論[16],北方河流冰封斷面可以劃分為河床影響區(qū)及冰蓋影響區(qū),具體見圖1。假設(shè)河道斷面平均流速與不同分區(qū)的平均流速一致,河道斷面總體水力坡度與分區(qū)的水力坡度一致,則可以通過以下方程分別計算不同分區(qū)的流量:

    式中:Q為流量;f為阻力系數(shù);R為水力半徑;A為面積;J為水力坡度;g為重力加速度;下標i、b分別指冰蓋影響區(qū)和河床影響區(qū),具體示意見圖1。結(jié)合式(1)和(2)可得冰蓋下的總流量為:

    若考慮部分河道斷面,如圖2垂線C-C左側(cè)斷面過流量,類似式(1)至(3)可以得到部分過流斷面的流量為:

    式中:下標p示意部分斷面水力參數(shù),具體示意見圖2。如果忽略C-C斷面水流間的內(nèi)摩擦力,在恒定均勻流條件下水流重力沿底坡的分量由床面和冰面的摩擦力平衡,則存在:

    式中:P為相應(yīng)區(qū)域的濕周;τ為壁面切應(yīng)力,下標與前面定義一致。圖1和2中冰面影響區(qū)與河床影響區(qū)的分界線上水流切應(yīng)力為零,因此式(7)中床面阻力與床面影響區(qū)的水流重力分量平衡,冰面阻力與冰面影響區(qū)的水流重力分量平衡,則:

    定義冰面與床面阻力系數(shù)之比為:

    則有:

    式中U為斷面平均流速。將式(8)除以(9)得:

    圖1 自然河流冰期測流示意圖

    圖2 冰封期過流斷面示意圖

    將式(6)除以式(3)可得:

    將式(10)(12)(13)代入式(14)可得:

    方程(15)顯示部分河道斷面流量與總流量之比是斷面面積、濕周、水力半徑及冰面與床面阻力系數(shù)比的函數(shù),這與Shen和Ackermann的結(jié)論一致。針對自然河流寬淺河道特征,Shen和Ackermann得到如下公式:

    式中:α為待率定的斷面特征系數(shù),為部分過流斷面與總過流斷面間濕周、水力半徑之比及冰面與床面阻力系數(shù)之比的函數(shù)。在實驗條件下,α一般取0.25;在自然河流中,α取值在1/2至2/3之間。

    2.2 測量步驟根據(jù)以上理論分析,冬季冰蓋下的流量觀測可采用如下步驟:

    (1)利用雙頻雷達[27-28]量測冰蓋下水深沿橫斷面的分布,或者在封河以前量測斷面地形高程,依據(jù)冰封期的水面高程計算橫斷面的水深分布。將斷面劃分為m-1個條形單元,記錄m個測點處的水深,具體參考圖2。

    (2)從左側(cè)端點計算起點距,利用數(shù)值差分計算每個條形單元的面積,計算公式如下:

    式中:Ai為第k個條形單元的面積;y為距離左岸的起點距;為y點處的水深,其中下標k,k+1分別為第k個條形單元的左右端點標號。

    (3)分別計算第i單元左側(cè)累計的斷面面積和累計水力半徑:

    式中:Ayk為1到k單元累計的斷面面積;Ryk為相應(yīng)的水力半徑;而χyk為對應(yīng)的濕周長度。

    (5)利用方程(23)計算左側(cè)相對單寬流量在各個測點的值:

    (6)從右側(cè)端點計算起點距,重復(fù)步驟(2)—(5),得到各個測點右側(cè)相對單寬流量的值,最后各個測點的單寬流量為:

    (7)在冰蓋任一測點k處鉆孔,測量0.2 倍水深處流速Uk0.2和0.8 倍水深處流速Uk0.8,k點處實際的單寬流量為:

    該斷面總的流量為:

    式中:U為某一垂線水深平均流速;下標0.2示意0.2倍水深處的流速;下標0.8示意0.8倍水深處的流速。k測點的水深平均流速為:

    任意一點的單寬流量為:

    3 應(yīng)用案例

    3.1 實驗驗證為了驗證本文新提出的冰期測流方法,選取文獻[25]實驗條件下復(fù)式斷面冰蓋下的水深平均流速驗證式(26)至式(28)的準確性[25]。實驗在一復(fù)式斷面順直水槽開展,長20 m,寬1 m。因為水槽沿中軸線對稱,選取一半斷面以簡化計算,具體斷面形態(tài)見圖3。他們采用塑料泡沫板模擬完全冰封和部分冰封工況,在恒定均勻流條件下,實驗測量了中間斷面水深平均流速分布,具體見文獻[25]。在α率定為0.25 的條件下,圖4 對比了方程(27)計算的斷面流速分布與實驗資料和文獻[25]的SKM方法計算值。由于實驗斷面沿中軸線對稱,圖4給出了一半斷面的水深平均流速分布。結(jié)果顯示本文計算結(jié)果與實驗觀測資料吻合良好。新方法計算準確性與文獻[25]的SKM方法類似,地形突變處新方法計算精度更高。新的冰期測量方法能應(yīng)用于暢流河段、部分冰封及完全冰封河道水深平均流速、單寬流量及總流量觀測。

    圖3 實驗條件復(fù)式斷面地形及水位分布

    圖4 新方法計算的水深平均流速分布與實測值及文獻[25]計算結(jié)果的對比

    3.2 黃河驗證采用實驗資料驗證本文新提出的冰期測量方法后,將該方法進一步應(yīng)用到天然河道部分冰蓋及完整冰蓋下的流量計算。王春青等在2014—2015的冬季開展了黃河內(nèi)蒙古河段的冰情觀測,詳細觀測了溫布壕、打不素太和三湖河口斷面的冰蓋厚度、水下地形及單寬流量分布,該測量結(jié)果為新方法的應(yīng)用提供了寶貴的原型資料[29]。在斷面特征系數(shù)α率定為0.5的基礎(chǔ)上,圖5顯示了溫布壕斷面2015年3月13日實測與計算的單寬流量分布以及觀測的冰厚及水下地形分布。其中,垂線高程坐標z顯示了冰蓋底部距離表面的厚度及河床距離表面的距離,圖6和圖7的坐標與此一致。結(jié)果顯示新方法計算的單寬流量與原型觀測值符合良好,在給定冰厚和河床地形條件下,式(28)能準確計算溫布壕斷面的單寬流量分布。圖6顯示了打不素太斷面2015年2月9日實測的冰蓋厚度、水下地形高程及采用鉆孔法測量的單寬流量分布。結(jié)合該圖給定的水下地形高程及冰蓋位置,采用2.2節(jié)的步驟1到7計算了相應(yīng)單寬流量分布。計算值與實測值的對比結(jié)果見圖6,兩者吻合良好。新方法能準確模擬實際復(fù)雜地形下的單寬流量分布,計算的深泓線上流量最大,與現(xiàn)場觀測一致??傮w而言,本文新提出的理論公式能準確模擬不同斷面、不同冰厚和復(fù)雜地形下的單寬流量分布和斷面總流量。新方法能極大減少冰蓋鉆孔的數(shù)量,進而提高冰期測流效率,具有良好的應(yīng)用前景。作為我國第二大河流,黃河每年都會經(jīng)歷100多天的冰期,準確觀測黃河冬季流量對防凌減災(zāi)具有重大意義。下一步研究可將冰水情一體化雷達和本文理論公式結(jié)合,以高效準確地獲取冰期流量,支撐黃河冬季原型觀測。

    圖5 黃河內(nèi)蒙古河段溫布壕斷面2015年3月13日實測冰蓋厚度、地形分布及單寬流量的驗證

    圖6 黃河內(nèi)蒙古河段打不素太斷面2015年2月9日實測冰蓋厚度、地形分布及單寬流量的驗證

    圖7 黃河內(nèi)蒙古河段三湖河口斷面2015年1月19日實測冰蓋厚度、地形分布及單寬流量的驗證

    圖7進一步顯示了三湖河口2015年1月19日部分冰封下實測與計算的單寬流量以及相應(yīng)觀測的表面冰厚及河床高程分布。結(jié)果顯示新方法能準確計算部分冰蓋下暢流區(qū)與冰蓋區(qū)的單寬流量,具有較高計算精度。此外,該方法能避免脆弱冰蓋上的鉆孔,減少寒冷天氣下的測量時間,進而保障測量人員安全,比傳統(tǒng)鉆孔測流方法適用范圍更廣。新的冰期測流方法能進一步推廣應(yīng)用到其它河流暢流、部分冰封及完全冰封條件下的流量、單寬流量及水深平均流速觀測。該方法僅需測量一個測線上兩點的流速,結(jié)合冰水情一體化雷達獲取的水深及冰厚,就能給出斷面總流量、單寬流量及水深平均流速分布,能有效提高冰期測流效率。需要指出的是,本文提出的方法忽略了環(huán)流和斷面間動量交換的影響,會影響公式在彎道的應(yīng)用。但通過斷面特征系數(shù)α的率定,該公式能普遍用于一般河流流量的觀測,具體彎道斷面的應(yīng)用尚待進一步研究。

    4 結(jié)論

    本文將愛因斯坦的床面和壁面分區(qū)方法推廣應(yīng)用到冰封斷面,將斷面過流區(qū)分為冰蓋影響區(qū)及河床影響區(qū);基于流管元法提出一種冰期測量水深平均流速、單寬流量和斷面總流量的新方法。新方法只需測量一條垂線上0.2和0.8倍水深兩點處的流速,能極大提高冰期測流的工作效率,降低工作量。新方法應(yīng)用于實驗室復(fù)式斷面水槽完全冰封和部分冰封的水深平均流速分布計算及黃河內(nèi)蒙古河段三個不同斷面完全冰封和部分冰封的單寬流量計算。結(jié)果顯示該方法具有較高精度和可靠性,能進一步推廣應(yīng)用到其它河流暢流、部分冰蓋及冰蓋下的流量觀測。

    猜你喜歡
    冰封冰蓋測流
    軍事文摘(2024年6期)2024-04-30 03:13:59
    神秘的冰封之灣
    渠道斷面自動測流系統(tǒng)在位山灌區(qū)測水量水中的應(yīng)用
    水文測流技術(shù)方法與進展分析
    石河子科技(2022年4期)2022-03-24 05:45:28
    格陵蘭島的冰蓋悄悄融化
    參花(下)(2022年1期)2022-01-15 00:45:01
    感受極地冰封之美
    冰封奇觀
    奇聞怪事(2021年2期)2021-04-20 12:12:01
    冰封在“空中”的湖魚
    曹店灌區(qū)渠首測流存在的問題及對策
    山東水利(2018年6期)2018-03-24 13:00:35
    長距離輸水工程的冰期冰蓋數(shù)值模擬研究
    国产伦精品一区二区三区四那| 黄色一级大片看看| 国内精品宾馆在线| 久久狼人影院| 精品国产国语对白av| 午夜91福利影院| 高清欧美精品videossex| 国产 精品1| 在线免费观看不下载黄p国产| 热99国产精品久久久久久7| 国产亚洲精品久久久com| 国产精品久久久久成人av| 夜夜骑夜夜射夜夜干| 偷拍熟女少妇极品色| 欧美一级a爱片免费观看看| 少妇精品久久久久久久| 免费少妇av软件| 亚洲国产毛片av蜜桃av| 日本欧美视频一区| 秋霞在线观看毛片| 欧美精品亚洲一区二区| 亚洲国产精品一区三区| 91精品伊人久久大香线蕉| 热99国产精品久久久久久7| 欧美丝袜亚洲另类| 熟妇人妻不卡中文字幕| 中文字幕免费在线视频6| 在线免费观看不下载黄p国产| 国产淫片久久久久久久久| 热99国产精品久久久久久7| 久久精品国产亚洲网站| av线在线观看网站| 国产深夜福利视频在线观看| 国产乱人偷精品视频| 成人综合一区亚洲| 日韩,欧美,国产一区二区三区| 国产精品久久久久久久电影| 只有这里有精品99| 日本vs欧美在线观看视频 | 一本—道久久a久久精品蜜桃钙片| 国国产精品蜜臀av免费| 青青草视频在线视频观看| 久久6这里有精品| 中文字幕精品免费在线观看视频 | 成人亚洲欧美一区二区av| 亚洲欧美日韩东京热| 少妇高潮的动态图| 夫妻性生交免费视频一级片| 各种免费的搞黄视频| 99九九在线精品视频 | 中文精品一卡2卡3卡4更新| 插阴视频在线观看视频| 亚洲国产精品国产精品| 性高湖久久久久久久久免费观看| 日日撸夜夜添| 黄色怎么调成土黄色| 插阴视频在线观看视频| 高清欧美精品videossex| 欧美日韩av久久| 97精品久久久久久久久久精品| 国产亚洲一区二区精品| 蜜桃久久精品国产亚洲av| 国产在线一区二区三区精| 在线观看一区二区三区激情| 91精品一卡2卡3卡4卡| 插逼视频在线观看| 免费观看a级毛片全部| 一级毛片电影观看| 在线观看免费高清a一片| 国产精品久久久久成人av| 少妇 在线观看| 极品教师在线视频| 在线亚洲精品国产二区图片欧美 | 麻豆乱淫一区二区| 人妻系列 视频| 极品教师在线视频| 国产精品麻豆人妻色哟哟久久| 国产又色又爽无遮挡免| 国产黄片美女视频| 国产免费福利视频在线观看| 国产精品三级大全| 嘟嘟电影网在线观看| 国产成人freesex在线| 欧美日韩在线观看h| 人人澡人人妻人| 亚洲精品第二区| 少妇被粗大的猛进出69影院 | 少妇的逼水好多| 91午夜精品亚洲一区二区三区| 国产一区亚洲一区在线观看| 免费黄色在线免费观看| 99国产精品免费福利视频| 亚洲真实伦在线观看| 中文精品一卡2卡3卡4更新| a级毛片免费高清观看在线播放| 国产黄片美女视频| 中文天堂在线官网| 各种免费的搞黄视频| 国产成人a∨麻豆精品| 在线免费观看不下载黄p国产| 新久久久久国产一级毛片| 久久av网站| 精品国产露脸久久av麻豆| 日韩人妻高清精品专区| 五月伊人婷婷丁香| 亚洲内射少妇av| 卡戴珊不雅视频在线播放| 久久人人爽人人片av| 熟女av电影| 亚洲精品,欧美精品| 有码 亚洲区| 亚洲情色 制服丝袜| 国产女主播在线喷水免费视频网站| 99精国产麻豆久久婷婷| 伦理电影大哥的女人| 99精国产麻豆久久婷婷| 成人国产麻豆网| 欧美日韩视频高清一区二区三区二| 国产高清有码在线观看视频| 成人影院久久| 久久精品久久精品一区二区三区| 少妇 在线观看| 久久狼人影院| 男女边摸边吃奶| 久久人人爽人人爽人人片va| 能在线免费看毛片的网站| 亚洲精品第二区| 国产乱来视频区| 天堂俺去俺来也www色官网| 极品少妇高潮喷水抽搐| 91久久精品电影网| 一级毛片久久久久久久久女| 亚洲av在线观看美女高潮| 全区人妻精品视频| 亚洲欧洲日产国产| 另类精品久久| 日本免费在线观看一区| 成人影院久久| 久久精品国产a三级三级三级| 你懂的网址亚洲精品在线观看| 一级a做视频免费观看| 国产视频内射| 人体艺术视频欧美日本| 欧美另类一区| 九九爱精品视频在线观看| 黑丝袜美女国产一区| 777米奇影视久久| 久久久久国产精品人妻一区二区| 日日摸夜夜添夜夜添av毛片| 久久久久久伊人网av| 在线 av 中文字幕| 精品人妻熟女毛片av久久网站| 国产 一区精品| 少妇猛男粗大的猛烈进出视频| 乱码一卡2卡4卡精品| 交换朋友夫妻互换小说| 国产高清有码在线观看视频| h视频一区二区三区| 午夜av观看不卡| 亚洲高清免费不卡视频| 日韩视频在线欧美| 街头女战士在线观看网站| 久久韩国三级中文字幕| a级一级毛片免费在线观看| 精华霜和精华液先用哪个| 插阴视频在线观看视频| 两个人免费观看高清视频 | 亚洲欧美清纯卡通| 亚洲,欧美,日韩| 九九久久精品国产亚洲av麻豆| 亚洲av免费高清在线观看| kizo精华| 亚洲欧洲国产日韩| 欧美xxxx性猛交bbbb| 欧美激情极品国产一区二区三区 | 久久精品久久久久久噜噜老黄| 一区二区三区乱码不卡18| 日韩精品免费视频一区二区三区 | 日本黄色日本黄色录像| 黑人高潮一二区| 国产中年淑女户外野战色| 国产精品一区二区三区四区免费观看| 国产精品欧美亚洲77777| 亚洲精品456在线播放app| 色网站视频免费| 九色成人免费人妻av| 另类精品久久| 国内少妇人妻偷人精品xxx网站| 嫩草影院新地址| videos熟女内射| 22中文网久久字幕| 91精品国产国语对白视频| av黄色大香蕉| 国产成人精品福利久久| 91久久精品电影网| 天天躁夜夜躁狠狠久久av| 看非洲黑人一级黄片| 欧美国产精品一级二级三级 | 插逼视频在线观看| 亚洲国产精品专区欧美| 啦啦啦视频在线资源免费观看| 亚洲精品中文字幕在线视频 | 亚洲国产欧美在线一区| 久久久久久久大尺度免费视频| 夜夜看夜夜爽夜夜摸| 成人影院久久| 又黄又爽又刺激的免费视频.| 中文欧美无线码| 国产成人精品无人区| 免费久久久久久久精品成人欧美视频 | 伊人久久国产一区二区| 亚洲av在线观看美女高潮| 人妻人人澡人人爽人人| 人妻少妇偷人精品九色| 最近2019中文字幕mv第一页| 有码 亚洲区| 日韩亚洲欧美综合| 涩涩av久久男人的天堂| 成年人免费黄色播放视频 | 少妇被粗大猛烈的视频| 欧美少妇被猛烈插入视频| 美女cb高潮喷水在线观看| 妹子高潮喷水视频| 国产精品国产av在线观看| 中文天堂在线官网| 中文精品一卡2卡3卡4更新| 一级二级三级毛片免费看| 春色校园在线视频观看| 国产亚洲一区二区精品| kizo精华| 在线观看免费日韩欧美大片 | 久久午夜综合久久蜜桃| 99热全是精品| 国产成人精品福利久久| 亚洲欧洲日产国产| 夜夜爽夜夜爽视频| 中文欧美无线码| av在线播放精品| 下体分泌物呈黄色| 香蕉精品网在线| tube8黄色片| 黄色配什么色好看| 中文欧美无线码| 日日摸夜夜添夜夜爱| 国产91av在线免费观看| 亚洲经典国产精华液单| 国产一区二区在线观看av| 热re99久久国产66热| 久久久久久人妻| 免费看光身美女| 精华霜和精华液先用哪个| 欧美性感艳星| 国产亚洲精品久久久com| 亚洲av成人精品一区久久| 国产av精品麻豆| 日产精品乱码卡一卡2卡三| 99九九线精品视频在线观看视频| 久久99热6这里只有精品| 日韩,欧美,国产一区二区三区| 国产一区二区在线观看av| 三级国产精品片| 你懂的网址亚洲精品在线观看| 久久青草综合色| 亚洲丝袜综合中文字幕| 欧美精品国产亚洲| av黄色大香蕉| 国产永久视频网站| 亚洲国产精品一区三区| 欧美日韩国产mv在线观看视频| 国产91av在线免费观看| 人妻一区二区av| 国产精品无大码| 少妇裸体淫交视频免费看高清| 欧美日韩亚洲高清精品| 国产av精品麻豆| 日日撸夜夜添| 国产 精品1| 亚洲欧美成人精品一区二区| 最新的欧美精品一区二区| 亚洲情色 制服丝袜| 国产精品伦人一区二区| 国产精品99久久久久久久久| 久久99精品国语久久久| 国产日韩一区二区三区精品不卡 | 一区二区av电影网| 免费观看的影片在线观看| 久久婷婷青草| 免费高清在线观看视频在线观看| 91在线精品国自产拍蜜月| 99热这里只有精品一区| 国产av码专区亚洲av| 国产精品女同一区二区软件| 乱人伦中国视频| 国产一区二区三区av在线| 简卡轻食公司| 国产av精品麻豆| 午夜久久久在线观看| 成年人午夜在线观看视频| 精品一区二区三卡| 午夜福利网站1000一区二区三区| 国产成人精品婷婷| 国产精品99久久99久久久不卡 | 精品一区二区三区视频在线| 久热这里只有精品99| 高清黄色对白视频在线免费看 | 两个人免费观看高清视频 | 免费在线观看成人毛片| 不卡视频在线观看欧美| 欧美bdsm另类| 91精品伊人久久大香线蕉| 久久久久久久久久久丰满| 天堂中文最新版在线下载| a级毛片免费高清观看在线播放| 亚洲成人一二三区av| 亚洲中文av在线| 日本vs欧美在线观看视频 | 国产亚洲5aaaaa淫片| 中文乱码字字幕精品一区二区三区| 大陆偷拍与自拍| 91精品伊人久久大香线蕉| 国产成人freesex在线| 成人漫画全彩无遮挡| 国产亚洲av片在线观看秒播厂| 亚洲电影在线观看av| 99热网站在线观看| 一本一本综合久久| 国产精品国产av在线观看| 蜜桃久久精品国产亚洲av| 欧美三级亚洲精品| 夫妻午夜视频| 日韩一区二区三区影片| 久久久久久久久久人人人人人人| 午夜久久久在线观看| 一级av片app| av.在线天堂| 激情五月婷婷亚洲| 中文字幕免费在线视频6| av.在线天堂| 丰满迷人的少妇在线观看| 亚洲av不卡在线观看| .国产精品久久| 99热网站在线观看| 亚洲欧美一区二区三区国产| 成人美女网站在线观看视频| 免费高清在线观看视频在线观看| 午夜老司机福利剧场| 精品视频人人做人人爽| 亚洲精品久久久久久婷婷小说| 六月丁香七月| h视频一区二区三区| 日本午夜av视频| 国产淫片久久久久久久久| 欧美日本中文国产一区发布| 日本av手机在线免费观看| 人人妻人人看人人澡| 国产亚洲欧美精品永久| 国产高清三级在线| 黑人高潮一二区| 国产色爽女视频免费观看| 成人国产av品久久久| 人妻系列 视频| freevideosex欧美| 韩国av在线不卡| 亚洲精品久久久久久婷婷小说| 国产黄色视频一区二区在线观看| 日产精品乱码卡一卡2卡三| 国产日韩欧美视频二区| 国产真实伦视频高清在线观看| 天堂俺去俺来也www色官网| 一区二区三区精品91| 三级国产精品片| 一区二区三区精品91| 国产亚洲一区二区精品| 五月伊人婷婷丁香| 久久精品国产亚洲av天美| 最新中文字幕久久久久| 亚洲图色成人| 婷婷色综合大香蕉| 亚洲国产成人一精品久久久| 精品一区二区免费观看| 国产精品麻豆人妻色哟哟久久| 黑人高潮一二区| 国产免费又黄又爽又色| 免费av中文字幕在线| 午夜免费男女啪啪视频观看| 亚洲av电影在线观看一区二区三区| 中文资源天堂在线| 久久6这里有精品| 亚洲成色77777| 大片免费播放器 马上看| 国产高清不卡午夜福利| 日本欧美国产在线视频| 一级,二级,三级黄色视频| 少妇熟女欧美另类| 国产精品无大码| 特大巨黑吊av在线直播| 精品久久久久久电影网| 日本欧美视频一区| 91久久精品国产一区二区三区| 91精品一卡2卡3卡4卡| 噜噜噜噜噜久久久久久91| 亚洲精华国产精华液的使用体验| 青青草视频在线视频观看| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品第二区| 欧美精品人与动牲交sv欧美| 国产中年淑女户外野战色| 久久99精品国语久久久| 18禁动态无遮挡网站| 国产精品成人在线| 亚洲国产色片| 亚洲av欧美aⅴ国产| 97精品久久久久久久久久精品| 国产男女超爽视频在线观看| 久久久久久久精品精品| 欧美激情极品国产一区二区三区 | 综合色丁香网| av有码第一页| 免费观看的影片在线观看| 亚洲国产精品一区二区三区在线| 国产中年淑女户外野战色| 91aial.com中文字幕在线观看| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 丝袜喷水一区| 丰满乱子伦码专区| 色5月婷婷丁香| 日本爱情动作片www.在线观看| 女性生殖器流出的白浆| 一区二区三区乱码不卡18| 三级国产精品欧美在线观看| 日韩一区二区三区影片| av福利片在线| 涩涩av久久男人的天堂| 又粗又硬又长又爽又黄的视频| 赤兔流量卡办理| 在线观看美女被高潮喷水网站| 国产精品一区二区三区四区免费观看| 亚洲欧美清纯卡通| 天天操日日干夜夜撸| 国产乱来视频区| 3wmmmm亚洲av在线观看| 又爽又黄a免费视频| 国产一区亚洲一区在线观看| 国产av一区二区精品久久| 一级毛片久久久久久久久女| 色94色欧美一区二区| 国产亚洲精品久久久com| 久久久久久久久久人人人人人人| 国产高清有码在线观看视频| 一级毛片黄色毛片免费观看视频| 99久国产av精品国产电影| 久久久亚洲精品成人影院| 国产视频内射| 亚洲国产av新网站| 丝袜喷水一区| 搡老乐熟女国产| 国产伦理片在线播放av一区| 美女福利国产在线| 成年美女黄网站色视频大全免费 | 在线观看免费视频网站a站| 五月开心婷婷网| 一二三四中文在线观看免费高清| 中文字幕人妻丝袜制服| 亚洲精品成人av观看孕妇| 婷婷色综合大香蕉| 国产精品麻豆人妻色哟哟久久| 丰满乱子伦码专区| 精品酒店卫生间| 久久国内精品自在自线图片| 日韩av在线免费看完整版不卡| 嘟嘟电影网在线观看| 日韩制服骚丝袜av| 午夜免费观看性视频| 日韩成人av中文字幕在线观看| av天堂中文字幕网| 国产黄色视频一区二区在线观看| 3wmmmm亚洲av在线观看| 日本欧美视频一区| 视频中文字幕在线观看| 高清午夜精品一区二区三区| 韩国av在线不卡| 又大又黄又爽视频免费| 免费少妇av软件| 青春草亚洲视频在线观看| 久久久国产精品麻豆| 亚洲精品日本国产第一区| 久久99蜜桃精品久久| 成年av动漫网址| 亚洲自偷自拍三级| 人妻一区二区av| 99九九在线精品视频 | 中文字幕制服av| 丰满人妻一区二区三区视频av| 九九久久精品国产亚洲av麻豆| 青春草视频在线免费观看| 在线观看美女被高潮喷水网站| 99精国产麻豆久久婷婷| 99热6这里只有精品| 成人国产麻豆网| 久久狼人影院| 亚洲,一卡二卡三卡| 亚洲欧美精品自产自拍| 韩国高清视频一区二区三区| 精品国产一区二区久久| 免费人成在线观看视频色| 久久久久久久精品精品| 制服丝袜香蕉在线| 久久精品久久久久久久性| 精品国产一区二区久久| 国产av一区二区精品久久| 啦啦啦视频在线资源免费观看| 在线看a的网站| 国产无遮挡羞羞视频在线观看| av免费在线看不卡| 午夜免费男女啪啪视频观看| 王馨瑶露胸无遮挡在线观看| 亚洲精品亚洲一区二区| 蜜桃久久精品国产亚洲av| 国产精品免费大片| 国产精品一区二区三区四区免费观看| 精品一区二区三区视频在线| 最黄视频免费看| 亚洲av成人精品一区久久| 人人澡人人妻人| av.在线天堂| 日本黄色片子视频| 久久久久视频综合| 十八禁网站网址无遮挡 | 亚洲欧美精品专区久久| 一个人免费看片子| 一边亲一边摸免费视频| 赤兔流量卡办理| 少妇的逼好多水| 国产精品欧美亚洲77777| 成年av动漫网址| 一区在线观看完整版| 国产精品蜜桃在线观看| 日本免费在线观看一区| 国产免费一区二区三区四区乱码| 色94色欧美一区二区| 美女内射精品一级片tv| 看免费成人av毛片| 午夜福利视频精品| 日韩成人伦理影院| 色婷婷av一区二区三区视频| 国产精品久久久久久久久免| 性高湖久久久久久久久免费观看| 高清av免费在线| 99久久精品一区二区三区| 亚洲欧洲日产国产| 日日啪夜夜爽| 日韩大片免费观看网站| 少妇人妻精品综合一区二区| 这个男人来自地球电影免费观看 | 亚洲欧洲精品一区二区精品久久久 | 一区二区av电影网| 少妇精品久久久久久久| 国产精品欧美亚洲77777| 欧美人与善性xxx| 18禁在线无遮挡免费观看视频| 成年美女黄网站色视频大全免费 | 熟女人妻精品中文字幕| 天天躁夜夜躁狠狠久久av| 亚洲精品成人av观看孕妇| 亚洲av国产av综合av卡| 欧美3d第一页| 久久精品久久久久久噜噜老黄| 久久久久久久国产电影| 久久久久久伊人网av| 精品一区二区三卡| 噜噜噜噜噜久久久久久91| 一个人看视频在线观看www免费| 久久6这里有精品| 99精国产麻豆久久婷婷| 国产精品熟女久久久久浪| 欧美精品一区二区免费开放| 国产成人精品一,二区| 中文字幕av电影在线播放| 欧美精品人与动牲交sv欧美| 久久精品国产鲁丝片午夜精品| 成人18禁高潮啪啪吃奶动态图 | 自拍偷自拍亚洲精品老妇| 国产日韩欧美亚洲二区| 黄色怎么调成土黄色| 男女免费视频国产| 九九爱精品视频在线观看| 久久久久人妻精品一区果冻| 永久网站在线| 自线自在国产av| 成人漫画全彩无遮挡| 久久99一区二区三区| 久久久久久久国产电影| 国产精品一区二区在线不卡| 中文字幕人妻丝袜制服| 亚洲精品国产av蜜桃| 噜噜噜噜噜久久久久久91| 制服丝袜香蕉在线| 内地一区二区视频在线| 丰满迷人的少妇在线观看| 久久99一区二区三区| 久久精品国产a三级三级三级| 国产精品女同一区二区软件| 亚洲精品aⅴ在线观看| 亚洲精品乱码久久久v下载方式| 一级二级三级毛片免费看| 国产女主播在线喷水免费视频网站| 国产在线一区二区三区精| 18禁在线播放成人免费| 交换朋友夫妻互换小说| 涩涩av久久男人的天堂| 一区二区三区四区激情视频| 日韩人妻高清精品专区| 国产高清三级在线| 国产成人a∨麻豆精品|