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

    無資料地區(qū)設(shè)計(jì)洪水計(jì)算方法應(yīng)用研究

    2021-03-30 09:08:22陸泳舟王世昭閆旭輝
    關(guān)鍵詞:匯流洪峰流量計(jì)算方法

    朱 勇,陸泳舟,王世昭,閆旭輝

    (上??睖y設(shè)計(jì)研究院有限公司,上海 200434)

    設(shè)計(jì)洪水是確定水利水電工程建設(shè)規(guī)模及制定運(yùn)行管理規(guī)劃的重要依據(jù),設(shè)計(jì)洪水的可靠性估計(jì)是工程設(shè)計(jì)安全性的基礎(chǔ)[1]。對(duì)于有長序列流量、水位等實(shí)測資料的工程,可以在開展資料系列三性審查的基礎(chǔ)上,利用水文頻率計(jì)算方法推求設(shè)計(jì)洪水。然而,國內(nèi)外很多地區(qū)由于水情測報(bào)系統(tǒng)建設(shè)較晚,水文站的覆蓋密度不大,觀測時(shí)限也不足30年,因此在水利水電工程設(shè)計(jì)時(shí)沒有直接可用的流量資料,這類地區(qū)本文稱為無資料地區(qū)[2]。對(duì)于無資料地區(qū),需要通過其他水文要素推求設(shè)計(jì)洪水[3],主要方法有兩種,一是由設(shè)計(jì)暴雨通過產(chǎn)匯流計(jì)算得到設(shè)計(jì)洪水,具體的計(jì)算方法包括推理公式法、瞬時(shí)單位線法、SCC法等[4];二是由附近多個(gè)水文站的實(shí)測流量探尋地區(qū)洪水規(guī)律,通過本地區(qū)洪水和降雨、地形等相關(guān)因子的關(guān)系計(jì)算設(shè)計(jì)洪水,具體的計(jì)算方法包括地區(qū)洪峰模數(shù)法、區(qū)域回歸法等[5]。

    由于不同地區(qū)的流域特性不同,水情、雨情測報(bào)系統(tǒng)完善程度也不同,加之水文要素本身具有隨機(jī)性,在國內(nèi)外水利水電工程設(shè)計(jì)領(lǐng)域,尚沒有通用的設(shè)計(jì)洪水計(jì)算方法。實(shí)際工作中,往往從有利于水安全保障的角度考慮,在多種方法的計(jì)算值中選最大者,或者排除明顯偏大或者偏小值,根據(jù)人為經(jīng)驗(yàn)確定。本文重點(diǎn)研究、分析無資料地區(qū)設(shè)計(jì)洪水計(jì)算中常用的幾種方法在我國不同地區(qū)的使用情況,以實(shí)際生產(chǎn)中常遇的問題為關(guān)注點(diǎn),為水利規(guī)劃設(shè)計(jì)中相關(guān)問題的解決提供參考。

    1 國內(nèi)常用計(jì)算方法概述

    1.1 推理公式法

    (1)計(jì)算方法與參數(shù)取值

    推理公式法又稱“合理化”法,是最早根據(jù)暴雨資料推求設(shè)計(jì)洪水流量的方法之一,至今已有一百多年歷史[6]。我國各省市水文手冊在處理降雨和產(chǎn)匯流規(guī)律的方法基本相同,計(jì)算公式如下:

    對(duì)于全面匯流,即tc>τ時(shí),

    Qm=0.278(Sp/τn-μ)F

    (1)

    對(duì)于部分匯流,即tc<τ時(shí),

    (2)

    計(jì)算參數(shù)中除了集水面積、匯流河長、平均比降等流域地形參數(shù)外,還包括產(chǎn)流參數(shù)μ和匯流參數(shù)m。在我國各省市水文手冊中,產(chǎn)流參數(shù)μ一般由下墊面類型給出相應(yīng)取值[7],匯流參數(shù)m對(duì)計(jì)算成果的敏感性較高,計(jì)算方法有兩種,一是采用水文手冊的計(jì)算公式,通常和流域地形特征參數(shù)表現(xiàn)為指數(shù)關(guān)系,二是根據(jù)SL44—2006《水利水電工程設(shè)計(jì)洪水計(jì)算規(guī)范》中表B.2.2的取值條件,綜合考慮地區(qū)、植被和流域形狀等因素取值。

    (2)實(shí)例

    福建漳州市某流域入??诩娣e88.5km2,匯流河長16km,平均比降5.29‰,設(shè)計(jì)暴雨由地區(qū)實(shí)測雨量確定。依據(jù)福建省水文手冊,匯流參數(shù)m計(jì)算值為0.129,10年一遇設(shè)計(jì)洪峰流量667m3/s;依據(jù)華東地區(qū)小流域洪水計(jì)算手冊,匯流參數(shù)m計(jì)算值為0.773,10年一遇洪峰流量510m3/s;參考SL44—2006中“南方地區(qū)、植被一般、耕地較多植被地區(qū)”的要求,匯流參數(shù)m計(jì)算值為0.9,10年一遇設(shè)計(jì)洪峰流量560m3/s,三者差異很大。在缺少水文資料的情況下,經(jīng)過實(shí)地訪問,根據(jù)該河道近幾十年來的水位情況,初步判斷SL44—2006的計(jì)算值較合理,可見除流域地形特征外,下墊面情況也是重要的影響因素。

    1.2 瞬時(shí)單位線法

    (1)計(jì)算方法與參數(shù)取值

    綜合瞬時(shí)單位線法是在納希(Nash)單位線的基礎(chǔ)上,借助有水文資料流域的單位線要素或瞬時(shí)單位線的參數(shù)與其流域自然地理因子之間的相關(guān)關(guān)系,推求無資料地區(qū)設(shè)計(jì)洪水的方法。瞬時(shí)單位線是指無窮小時(shí)段內(nèi)流域上均勻單位凈雨在出口斷面所形成的地面徑流過程,其數(shù)學(xué)表達(dá)式為:

    (3)

    式中,μ(o,t)—t時(shí)刻瞬時(shí)單位線的縱高;Γ(n)—n階的伽瑪函數(shù);n—流域調(diào)蓄能力,相當(dāng)于線性水庫的個(gè)數(shù);K—調(diào)蓄系數(shù),相當(dāng)于流域匯流時(shí)間的參數(shù);e—自然對(duì)數(shù)的底。

    上式中計(jì)算參數(shù)較多,計(jì)算方法有兩種,一是在我國各省水文手冊中,根據(jù)實(shí)測資料,建立了不同地區(qū)瞬時(shí)單位線參數(shù)和地形參數(shù)的關(guān)系式,以此為計(jì)算依據(jù);二是根據(jù)流域數(shù)字地面高程模型直接計(jì)算[8- 9],如河海大學(xué)石朋、芮孝芳等經(jīng)過長期深入研究[10]提出的綜合運(yùn)用坡面匯流和河道匯流的計(jì)算方法。

    (2)實(shí)例

    鎮(zhèn)江市通勝地區(qū)某河道流域面積150.26km2,匯流河長27.7km,平均比降14.2‰,設(shè)計(jì)暴雨由地區(qū)相關(guān)規(guī)劃確定。根據(jù)江蘇省水文手冊中的瞬時(shí)單位線計(jì)算方法,50年一遇洪峰流量1467m3/s,20年一遇洪峰流量1186m3/s。根據(jù)數(shù)字地形高程推算的瞬時(shí)單位線進(jìn)行計(jì)算[11],50年一遇洪峰流量2655m3/s,20年一遇洪峰流量2192m3/s。兩種方法的結(jié)果差異很大,瞬時(shí)單位線對(duì)比如圖1所示。從該地區(qū)水文規(guī)律來看,后者更合理,原因在于江蘇省水文手冊編制時(shí)間距今已有近40年,受當(dāng)時(shí)的計(jì)算技術(shù)和資料儲(chǔ)備限制,其對(duì)瞬時(shí)單位線的參數(shù)的取值僅以蘇南山丘、混合區(qū)和蘇北山丘、混合區(qū)區(qū)分計(jì)算,顯得過于籠統(tǒng),且不考慮匯流河長對(duì)匯流參數(shù)的影響,而根據(jù)數(shù)字地形高程推算的瞬時(shí)單位線,充分考慮了現(xiàn)狀的地形地貌以及河網(wǎng)形態(tài),對(duì)匯流時(shí)間的模擬更真實(shí)。

    圖1 兩種方法瞬時(shí)單位線對(duì)比圖

    1.3 區(qū)域回歸法

    (1)計(jì)算方法與參數(shù)取值

    對(duì)于無資料地區(qū),設(shè)計(jì)洪水的推求是通過建立同一水文子區(qū)域內(nèi)鄰近站點(diǎn)洪峰流量與流域特征值的回歸方程來實(shí)現(xiàn)的[12]。

    首先是劃分水文區(qū)域。確定不同流域站點(diǎn)的洪水系列是否具有“水文相似性”是區(qū)域洪水頻率分析中一個(gè)關(guān)鍵的環(huán)節(jié)。通過選擇聚類因子進(jìn)行分類。先根據(jù)流域地貌和氣候特征的相似情況,初步選定幾個(gè)子區(qū)域,再根據(jù)過往的區(qū)域化研究成果和影響徑流的物理特征值綜合確定候選特征變量。然后通過一個(gè)循環(huán)過程來劃定區(qū)域邊界和選定特征變量。特征變量采用最小二乘法(OLS)來初選,所選的特征變量間應(yīng)避免相關(guān)性,在取舍具有相關(guān)性的特征變量時(shí)應(yīng)將其對(duì)回歸方程的有效性影響作為主要標(biāo)準(zhǔn)。對(duì)初選的每個(gè)水文分區(qū),建立多元回歸方程。在做回歸分析前,將所有的流量和流域特征值均取以10為底的對(duì)數(shù)。假定因變量為一個(gè)或多個(gè)自變量的線性函數(shù),基于所有變量對(duì)數(shù)轉(zhuǎn)換后的回歸關(guān)系可以表示為:

    lgQT=lgK+algA+blgB+…+nlgN

    (4)

    式中,QT—T年一遇的洪峰流量估值;K—回歸常數(shù);A、B、N—流域特征值;a、b、n—回歸系數(shù)。

    經(jīng)多次重復(fù)劃定子區(qū)域邊界,做OLS/GLS回歸,評(píng)價(jià)其殘差,最后可以得到合理的水文子區(qū)域。GLS回歸方程的特征變量一般為流域面積、平均流域坡度、降雨強(qiáng)度指數(shù)和主河道坡度等。

    通用的洪峰流量與流域特征變量的關(guān)系式如下:

    Y=xB+ε

    (5)

    式中,Y—T年一遇峰量事件組成的(n×1)向量;x—流域和氣候特征值組成的(n×p)矩陣;B—未知的回歸系數(shù)組成的(p×1)向量;ε—(n×1)隨機(jī)誤差向量。

    在確定了回歸方程后,需要計(jì)算方程的回歸系數(shù),常用的計(jì)算方法有普通最小二乘法(OLS)和加權(quán)最小二乘法(WLS)。在進(jìn)行回歸分析時(shí),需要根據(jù)特征變量的顯著性逐步剔除顯著性較小的特征變量,逐步得到最終回歸方程的系數(shù),并依據(jù)預(yù)測誤差的平方和最小或預(yù)測標(biāo)準(zhǔn)差最小的原則確定GLS回歸方程。

    (2)實(shí)例

    福建省萩蘆溪流域某水庫工程位于萩蘆溪干流上游,庫區(qū)上游河道以峽谷為主。工程壩址控制流域面積為75.8km2,主河長14.8km,比降43.61‰。萩蘆溪流域內(nèi)沒有水文測站,鄰近的水文測站分布位置如圖2所示。鄰近流域有渡里、瀨溪、東圳、太平口等水文站,均有30年以上流量測驗(yàn)序列和歷史洪水調(diào)查記錄,即各站均有合理可靠的設(shè)計(jì)洪水成果。

    圖2 工程和水文站位置示意圖

    本次研究以反映流域物理特性的洪峰統(tǒng)計(jì)參數(shù)(Cv、Cs、一階樣本矩L)相近為原則,確定瀨溪、東圳、太平口、東張、造水和赤橋站與壩址處可作為一個(gè)水文子區(qū)域。選取流域面積、河道縱坡、平均高程、年均降雨量、年均徑流量等5項(xiàng)指標(biāo)作為評(píng)價(jià)指標(biāo),對(duì)水文子區(qū)域建立區(qū)域回歸方程,然后對(duì)數(shù)線性化后,采用OLS多元回歸分析各指標(biāo)的顯著性。通過逐步OLS多元回歸分析,剔除顯著性較小的指標(biāo),最終確定流域面積(F)、河道比降(J)、平均高程(H)、年平均降雨量(P)為流域特征變量。區(qū)域回歸方程如下:

    lgQj,T=α0+α1lgFj+α2lgJj+α3lgHj+α4lgPj

    (6)

    式中,Qj,T—第j個(gè)流域T年一遇洪峰流量;α0—回歸常數(shù);α1、α2、α3、α4—回歸系數(shù)。

    經(jīng)最小二乘法確定回歸常數(shù)α0和回歸系數(shù)α1、α2、α3、α4,計(jì)算得壩址處100年一遇設(shè)計(jì)洪峰流量1186m3/s,50年一遇設(shè)計(jì)洪峰流量968m3/s。

    作為研究對(duì)比,如果僅以洪峰模數(shù)估算設(shè)計(jì)洪水,建立地區(qū)各站集水面積~洪峰流量對(duì)數(shù)關(guān)系,則計(jì)算得壩址處100年一遇設(shè)計(jì)洪峰流量1026m3/s,50年一遇設(shè)計(jì)洪峰流量849m3/s。相對(duì)而言,區(qū)域回歸法計(jì)算值比洪峰模數(shù)法大10%~15%,結(jié)果更合理,原因是該方法綜合考慮了集水面積、河道比降、降雨量等諸多因素對(duì)洪水的影響,尤其在浙南、閩南等山區(qū),相近地區(qū)不同流域往往河道比降、暴雨特性有較大差異,因此僅僅考慮集水面積,用周邊地區(qū)其他站點(diǎn)的洪峰模數(shù)估算設(shè)計(jì)洪水存在片面性。

    2 計(jì)算方法的適用條件分析

    2.1 推理公式法

    推理公式法的結(jié)構(gòu)形式簡單,但理論假定是流域內(nèi)產(chǎn)流強(qiáng)度在時(shí)間和空間上保持恒定不變,即降雨量和地面下滲量在空間分布上具備一致性,應(yīng)用條件比較局限。在我國多數(shù)省份的水文手冊中,推理公式的適用范圍都限制在匯水面積小于100km2以下的流域,同時(shí),在實(shí)際使用中匯流參數(shù)m敏感性高,對(duì)洪水值影響很大。因此,推理公式在使用中,在依據(jù)水文手冊計(jì)算的基礎(chǔ)上,對(duì)于下墊面組成復(fù)雜、或者下墊面條件明顯有變化的流域,應(yīng)當(dāng)實(shí)地調(diào)查,分析計(jì)算結(jié)果的合理性。尤其是很多山丘區(qū)村莊經(jīng)改造后,坡面匯流變成管網(wǎng)匯流[13],這種情況要格外重視。

    2.2 瞬時(shí)單位線法

    瞬時(shí)單位線法不考慮凈雨和下墊面的不均勻性,其適用范圍遠(yuǎn)大于推理公式法,一般適用于800~1000km2以下的流域,在設(shè)計(jì)工作中被廣泛使用。由于我國多數(shù)省份的水文手冊編制于二十世紀(jì)八九十年代,而隨著用地條件的變化和河網(wǎng)水系的調(diào)整,有些地區(qū)產(chǎn)匯流特性較之當(dāng)初已經(jīng)有了明顯變化,手冊中匯流參數(shù)的計(jì)算方法已不滿足現(xiàn)狀需求[14]。因此,在使用瞬時(shí)單位線法時(shí),參數(shù)取值不應(yīng)拘泥于地區(qū)水文手冊。用數(shù)字地形高程推算瞬時(shí)單位線的方法需要多次雨洪數(shù)據(jù)率定參數(shù),無資料地區(qū)需要臨時(shí)建站測驗(yàn),工作量比較大,暫無法普及運(yùn)用。對(duì)于防洪任務(wù)重、投資大的工程,應(yīng)充分重視瞬時(shí)單位線成果的合理性分析。

    2.3 區(qū)域回歸法

    區(qū)域回歸法是根據(jù)鄰近地區(qū)洪水和高程、比降、降雨等物理因素的規(guī)律計(jì)算設(shè)計(jì)洪水,對(duì)數(shù)據(jù)量的要求很高,但該方法避免了雨洪模型計(jì)算中諸如雨量點(diǎn)面折算系數(shù)、下滲率、匯流時(shí)間等諸多參數(shù)計(jì)算中存在的失真風(fēng)險(xiǎn),并充分考慮了歷史洪水的作用,而歷史暴雨調(diào)查的難度遠(yuǎn)大于歷史洪水,因此其應(yīng)用不局限于對(duì)降雨空間均勻性較高的中小流域。在設(shè)計(jì)過程中,應(yīng)認(rèn)識(shí)到區(qū)域回歸法的優(yōu)勢,搜集工程鄰近水文站的洪水?dāng)?shù)據(jù),分析各水文站洪水序列的一致性和歷史洪水重現(xiàn)期,在合理確定鄰近各站設(shè)計(jì)洪水和各因素回歸系數(shù)的基礎(chǔ)上計(jì)算設(shè)計(jì)洪水。

    3 結(jié)語

    本文以設(shè)計(jì)工作中最常用的三種計(jì)算方法為研究對(duì)象,闡述了實(shí)際應(yīng)用中面臨的問題。通過多個(gè)工程實(shí)例的洪水計(jì)算分析,采用推理公式、瞬時(shí)單位線法等雨洪模型計(jì)算設(shè)計(jì)洪水時(shí),應(yīng)重視流域?qū)嶋H地貌條件和水文手冊、相關(guān)規(guī)范中取值要求的契合度,采用地區(qū)洪水規(guī)律計(jì)算設(shè)計(jì)洪水時(shí),應(yīng)綜合考慮氣候、地形等多因素的影響。必要時(shí),盡可能采取多種方法計(jì)算,并結(jié)合實(shí)地訪問,確定設(shè)計(jì)洪水的合理值。本文所列實(shí)例均為中小流域天然產(chǎn)匯流情況,對(duì)于來水情況復(fù)雜的無資料地區(qū)[15],目前研究尚淺,將在之后的工作中繼續(xù)推進(jìn)。

    猜你喜歡
    匯流洪峰流量計(jì)算方法
    浮力計(jì)算方法匯集
    退耕還林工程對(duì)渭河洪峰流量的影響
    流域匯流計(jì)算方法在風(fēng)峪沙河段的應(yīng)用探討
    佛岡縣潖江流域年洪峰流量P-Ⅲ分布參數(shù)估算
    隨機(jī)振動(dòng)試驗(yàn)包絡(luò)計(jì)算方法
    大南川流域設(shè)計(jì)洪峰流量計(jì)算分析
    某特小流域設(shè)計(jì)洪峰流量計(jì)算分析
    不同應(yīng)變率比值計(jì)算方法在甲狀腺惡性腫瘤診斷中的應(yīng)用
    基于ARM CortexM0的智能光伏匯流采集裝置設(shè)計(jì)及應(yīng)用
    一種球載雷達(dá)匯流環(huán)設(shè)計(jì)
    久久久久精品久久久久真实原创| 男男h啪啪无遮挡| 午夜日本视频在线| 午夜视频国产福利| 日韩一区二区三区影片| 久久97久久精品| 欧美国产精品一级二级三级| av免费观看日本| 熟女电影av网| 欧美激情国产日韩精品一区| 新久久久久国产一级毛片| 国产成人aa在线观看| 国产视频首页在线观看| 少妇的丰满在线观看| 欧美日韩精品成人综合77777| 母亲3免费完整高清在线观看 | av视频免费观看在线观看| 在线观看www视频免费| 亚洲欧美一区二区三区国产| 国产成人aa在线观看| 欧美 亚洲 国产 日韩一| 亚洲内射少妇av| 街头女战士在线观看网站| 亚洲一区二区三区欧美精品| 狠狠婷婷综合久久久久久88av| 亚洲成av片中文字幕在线观看 | 五月玫瑰六月丁香| 街头女战士在线观看网站| 人妻少妇偷人精品九色| 国产成人精品无人区| 成年人免费黄色播放视频| 欧美精品亚洲一区二区| xxxhd国产人妻xxx| 亚洲精品久久久久久婷婷小说| 成年美女黄网站色视频大全免费| 在线免费观看不下载黄p国产| 一级片'在线观看视频| av免费观看日本| 亚洲精品久久久久久婷婷小说| 精品国产露脸久久av麻豆| 宅男免费午夜| 精品亚洲成a人片在线观看| 久久人人97超碰香蕉20202| 三级国产精品片| 五月天丁香电影| 欧美精品亚洲一区二区| 国产精品嫩草影院av在线观看| 午夜福利视频精品| 99久久综合免费| 在线观看免费高清a一片| 伦理电影免费视频| 男女国产视频网站| 色哟哟·www| 国产成人精品婷婷| 国产国拍精品亚洲av在线观看| 九色成人免费人妻av| 国产精品免费大片| 丰满乱子伦码专区| 亚洲av电影在线观看一区二区三区| 制服人妻中文乱码| 国产日韩一区二区三区精品不卡| 男女国产视频网站| 97在线视频观看| 欧美精品av麻豆av| 制服人妻中文乱码| 日本av手机在线免费观看| 亚洲色图综合在线观看| videossex国产| 欧美 亚洲 国产 日韩一| 下体分泌物呈黄色| 久久国内精品自在自线图片| 色婷婷久久久亚洲欧美| 97人妻天天添夜夜摸| 国产精品熟女久久久久浪| 一本久久精品| 热re99久久国产66热| 黄片播放在线免费| 亚洲av.av天堂| 国产麻豆69| 国产一区亚洲一区在线观看| 色婷婷av一区二区三区视频| 亚洲伊人久久精品综合| www.色视频.com| 久久精品久久久久久噜噜老黄| 成人手机av| 国产 一区精品| 一区二区av电影网| 日日撸夜夜添| 桃花免费在线播放| 99久久中文字幕三级久久日本| 免费播放大片免费观看视频在线观看| 日韩精品有码人妻一区| 精品国产一区二区久久| 国产 一区精品| 亚洲天堂av无毛| 国产精品国产三级专区第一集| 亚洲av日韩在线播放| 久久人妻熟女aⅴ| av网站免费在线观看视频| 精品福利永久在线观看| 日韩av不卡免费在线播放| 韩国av在线不卡| 自拍欧美九色日韩亚洲蝌蚪91| 日韩制服骚丝袜av| 18禁裸乳无遮挡动漫免费视频| 亚洲精品色激情综合| 久久ye,这里只有精品| 久久热在线av| 黑人欧美特级aaaaaa片| 色视频在线一区二区三区| 国精品久久久久久国模美| 看免费av毛片| 亚洲色图综合在线观看| 国产成人免费无遮挡视频| 国产成人a∨麻豆精品| 午夜91福利影院| 中文乱码字字幕精品一区二区三区| 亚洲国产色片| 51国产日韩欧美| 国产无遮挡羞羞视频在线观看| 超碰97精品在线观看| 女人被躁到高潮嗷嗷叫费观| 极品少妇高潮喷水抽搐| 激情视频va一区二区三区| 美女视频免费永久观看网站| 免费黄网站久久成人精品| 狠狠精品人妻久久久久久综合| 亚洲av在线观看美女高潮| 国产色爽女视频免费观看| 在线观看一区二区三区激情| 91成人精品电影| 久久精品国产a三级三级三级| 免费黄频网站在线观看国产| 日韩欧美一区视频在线观看| 国产熟女午夜一区二区三区| 国产高清不卡午夜福利| 亚洲人成77777在线视频| 老女人水多毛片| av在线老鸭窝| 青春草国产在线视频| 另类精品久久| 黄色 视频免费看| 亚洲av在线观看美女高潮| 亚洲人成网站在线观看播放| 99re6热这里在线精品视频| 亚洲欧美成人精品一区二区| 91精品三级在线观看| 国产一区二区激情短视频 | 日韩制服骚丝袜av| 精品久久久久久电影网| 精品久久国产蜜桃| 精品第一国产精品| 9191精品国产免费久久| 51国产日韩欧美| 97在线人人人人妻| 80岁老熟妇乱子伦牲交| kizo精华| 看免费成人av毛片| 精品国产一区二区三区久久久樱花| 国产av码专区亚洲av| 亚洲一码二码三码区别大吗| 欧美精品人与动牲交sv欧美| 日韩伦理黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品一区二区在线不卡| 国产精品久久久久久av不卡| 水蜜桃什么品种好| 久久久久久人人人人人| 久久久久久久亚洲中文字幕| 欧美+日韩+精品| 中文字幕人妻丝袜制服| 2022亚洲国产成人精品| 99热国产这里只有精品6| 国产国拍精品亚洲av在线观看| 免费观看性生交大片5| 一本久久精品| 亚洲伊人色综图| 人人妻人人爽人人添夜夜欢视频| 亚洲一区二区三区欧美精品| 久久精品久久久久久噜噜老黄| 黄色一级大片看看| 亚洲天堂av无毛| 91精品国产国语对白视频| 九九爱精品视频在线观看| 岛国毛片在线播放| 少妇被粗大猛烈的视频| 亚洲国产欧美在线一区| av不卡在线播放| 免费黄频网站在线观看国产| av网站免费在线观看视频| 有码 亚洲区| 亚洲欧美一区二区三区国产| 伦精品一区二区三区| 五月伊人婷婷丁香| 久久午夜福利片| 久久久a久久爽久久v久久| 欧美xxⅹ黑人| www.熟女人妻精品国产 | 成人免费观看视频高清| 99久久中文字幕三级久久日本| 中文乱码字字幕精品一区二区三区| 亚洲国产毛片av蜜桃av| 欧美xxⅹ黑人| 色吧在线观看| 欧美国产精品va在线观看不卡| 激情视频va一区二区三区| 少妇的逼水好多| 22中文网久久字幕| 最近最新中文字幕免费大全7| 日本av手机在线免费观看| a 毛片基地| 亚洲欧美色中文字幕在线| 国产熟女欧美一区二区| 精品久久久久久电影网| 亚洲欧美日韩卡通动漫| 中文字幕亚洲精品专区| 日韩电影二区| 色网站视频免费| 桃花免费在线播放| 好男人视频免费观看在线| 乱人伦中国视频| 看十八女毛片水多多多| 国产一区二区在线观看av| 9191精品国产免费久久| 国产日韩欧美亚洲二区| 国产高清不卡午夜福利| 一区二区av电影网| 亚洲国产最新在线播放| tube8黄色片| 国产精品国产三级专区第一集| 一级毛片我不卡| 18+在线观看网站| 亚洲国产精品专区欧美| 黄色 视频免费看| 国产女主播在线喷水免费视频网站| 亚洲精品久久午夜乱码| 午夜福利乱码中文字幕| 成人国语在线视频| 9热在线视频观看99| 天天影视国产精品| 日日撸夜夜添| 自线自在国产av| 这个男人来自地球电影免费观看 | 黑人欧美特级aaaaaa片| 午夜视频国产福利| 99久久中文字幕三级久久日本| 亚洲成人一二三区av| 国产福利在线免费观看视频| 亚洲综合精品二区| 国产免费现黄频在线看| 一区二区三区四区激情视频| 亚洲,欧美,日韩| 国产亚洲欧美精品永久| 久久精品国产a三级三级三级| 一区二区三区精品91| av不卡在线播放| 伦理电影大哥的女人| 狂野欧美激情性xxxx在线观看| 久久久国产欧美日韩av| 国产国语露脸激情在线看| 国产精品国产三级国产专区5o| 制服丝袜香蕉在线| 夜夜骑夜夜射夜夜干| 欧美日本中文国产一区发布| 国产日韩一区二区三区精品不卡| 在线观看一区二区三区激情| 99热网站在线观看| 在线亚洲精品国产二区图片欧美| 如日韩欧美国产精品一区二区三区| 欧美日韩av久久| 国产av精品麻豆| av国产精品久久久久影院| 两个人免费观看高清视频| 亚洲综合色惰| 久久精品国产鲁丝片午夜精品| 亚洲色图 男人天堂 中文字幕 | 国产av码专区亚洲av| 久久精品久久精品一区二区三区| 日韩成人伦理影院| 丰满乱子伦码专区| 男人添女人高潮全过程视频| 亚洲av男天堂| av线在线观看网站| 久久精品国产自在天天线| 在线免费观看不下载黄p国产| 国产成人精品久久久久久| 精品少妇内射三级| 国产欧美亚洲国产| 久久青草综合色| 日韩在线高清观看一区二区三区| 久久久精品94久久精品| 亚洲精品久久成人aⅴ小说| 国产成人免费观看mmmm| 精品人妻一区二区三区麻豆| 极品人妻少妇av视频| 国产精品久久久久久精品古装| 亚洲激情五月婷婷啪啪| 狂野欧美激情性bbbbbb| 日韩人妻精品一区2区三区| 又黄又粗又硬又大视频| av.在线天堂| 在线观看免费高清a一片| 欧美+日韩+精品| 欧美精品av麻豆av| 蜜桃在线观看..| 免费久久久久久久精品成人欧美视频 | 日韩制服丝袜自拍偷拍| 日韩精品有码人妻一区| 内地一区二区视频在线| 国产国语露脸激情在线看| 婷婷成人精品国产| 成人黄色视频免费在线看| a级毛色黄片| 精品人妻偷拍中文字幕| 国产黄色视频一区二区在线观看| 日韩一本色道免费dvd| www日本在线高清视频| 全区人妻精品视频| 春色校园在线视频观看| 欧美xxⅹ黑人| 国产精品嫩草影院av在线观看| 精品国产露脸久久av麻豆| 黄色毛片三级朝国网站| 伦理电影大哥的女人| 久久青草综合色| 欧美人与善性xxx| 制服人妻中文乱码| 国产毛片在线视频| 人人澡人人妻人| 久久亚洲国产成人精品v| 亚洲欧洲日产国产| 欧美变态另类bdsm刘玥| av天堂久久9| 久久久精品免费免费高清| 国产亚洲欧美精品永久| 久热这里只有精品99| 高清av免费在线| 精品一区二区三卡| 色哟哟·www| 男女无遮挡免费网站观看| av又黄又爽大尺度在线免费看| 男女午夜视频在线观看 | 成年女人在线观看亚洲视频| 午夜久久久在线观看| 天天影视国产精品| 少妇的逼好多水| 一级爰片在线观看| av又黄又爽大尺度在线免费看| 久久久a久久爽久久v久久| 一级毛片 在线播放| av黄色大香蕉| 国产一区亚洲一区在线观看| 久久狼人影院| 激情五月婷婷亚洲| 蜜桃在线观看..| 成人影院久久| 亚洲性久久影院| 国产精品女同一区二区软件| 亚洲成av片中文字幕在线观看 | 国语对白做爰xxxⅹ性视频网站| 国产男女超爽视频在线观看| av电影中文网址| 美女脱内裤让男人舔精品视频| 丰满少妇做爰视频| 午夜福利在线观看免费完整高清在| 咕卡用的链子| 日本黄大片高清| 国内精品宾馆在线| 精品人妻偷拍中文字幕| 三级国产精品片| 大香蕉久久成人网| 欧美日韩亚洲高清精品| 日韩精品有码人妻一区| 国产日韩欧美亚洲二区| 777米奇影视久久| 99热网站在线观看| 伊人久久国产一区二区| 免费不卡的大黄色大毛片视频在线观看| 在线观看人妻少妇| 久久久国产一区二区| 日韩精品免费视频一区二区三区 | 久久久国产欧美日韩av| 午夜激情久久久久久久| 国产精品国产av在线观看| 欧美少妇被猛烈插入视频| 热99久久久久精品小说推荐| 色视频在线一区二区三区| 国产在线视频一区二区| 午夜福利视频精品| 全区人妻精品视频| 9热在线视频观看99| 天天操日日干夜夜撸| 水蜜桃什么品种好| 久久99蜜桃精品久久| 亚洲av综合色区一区| av免费观看日本| 中文字幕免费在线视频6| 亚洲美女视频黄频| 国产一区有黄有色的免费视频| 亚洲国产欧美在线一区| 午夜视频国产福利| 视频中文字幕在线观看| 看十八女毛片水多多多| 欧美人与性动交α欧美软件 | 久久精品久久久久久久性| 亚洲熟女精品中文字幕| av线在线观看网站| 亚洲av福利一区| 一级爰片在线观看| 国产男女内射视频| 欧美激情 高清一区二区三区| 观看美女的网站| 亚洲精品美女久久av网站| 美女福利国产在线| 永久网站在线| 99热6这里只有精品| 久久99蜜桃精品久久| 欧美 日韩 精品 国产| 亚洲av成人精品一二三区| 国产亚洲av片在线观看秒播厂| 免费不卡的大黄色大毛片视频在线观看| 成人综合一区亚洲| 91久久精品国产一区二区三区| 国产毛片在线视频| 99久久综合免费| 免费观看av网站的网址| 亚洲人与动物交配视频| 成人黄色视频免费在线看| 日韩制服丝袜自拍偷拍| 在线观看三级黄色| a级片在线免费高清观看视频| 午夜免费鲁丝| 免费高清在线观看日韩| 中文乱码字字幕精品一区二区三区| 亚洲激情五月婷婷啪啪| 晚上一个人看的免费电影| 亚洲国产av影院在线观看| 午夜影院在线不卡| 亚洲欧洲精品一区二区精品久久久 | 99久久综合免费| 狠狠婷婷综合久久久久久88av| 性色avwww在线观看| 男女下面插进去视频免费观看 | 搡女人真爽免费视频火全软件| 男女午夜视频在线观看 | 日韩,欧美,国产一区二区三区| 黄色毛片三级朝国网站| 永久网站在线| 黄网站色视频无遮挡免费观看| 成人黄色视频免费在线看| 乱人伦中国视频| 欧美精品国产亚洲| 亚洲精品一二三| 国精品久久久久久国模美| 国产av码专区亚洲av| 水蜜桃什么品种好| 哪个播放器可以免费观看大片| 久久精品国产鲁丝片午夜精品| 亚洲综合色网址| 日韩成人伦理影院| 亚洲婷婷狠狠爱综合网| av网站免费在线观看视频| 热re99久久精品国产66热6| 26uuu在线亚洲综合色| av在线老鸭窝| 国产精品一区二区在线不卡| 少妇高潮的动态图| 免费日韩欧美在线观看| 国产精品国产三级专区第一集| 久热久热在线精品观看| 国产精品成人在线| 亚洲国产成人一精品久久久| av又黄又爽大尺度在线免费看| 久久久欧美国产精品| 亚洲欧美日韩卡通动漫| 日本黄大片高清| 日韩一区二区三区影片| 精品久久久久久电影网| 精品人妻在线不人妻| 最近中文字幕高清免费大全6| 大陆偷拍与自拍| 曰老女人黄片| 欧美日韩视频精品一区| 新久久久久国产一级毛片| 国国产精品蜜臀av免费| 黄色配什么色好看| 五月天丁香电影| 中文欧美无线码| 亚洲少妇的诱惑av| 狂野欧美激情性xxxx在线观看| 女性生殖器流出的白浆| 久久久久精品性色| 国产熟女午夜一区二区三区| 午夜福利视频在线观看免费| 国产高清国产精品国产三级| 亚洲经典国产精华液单| 亚洲欧美一区二区三区国产| 免费在线观看黄色视频的| 夜夜爽夜夜爽视频| 久久久久久久精品精品| av.在线天堂| 三级国产精品片| 日韩制服骚丝袜av| 国产午夜精品一二区理论片| 久久 成人 亚洲| 成人毛片60女人毛片免费| 欧美3d第一页| 日韩中文字幕视频在线看片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩免费高清中文字幕av| a级毛片在线看网站| 校园人妻丝袜中文字幕| 成年美女黄网站色视频大全免费| 99热这里只有是精品在线观看| 精品久久久久久电影网| 免费看不卡的av| 好男人视频免费观看在线| 日韩av不卡免费在线播放| 免费黄网站久久成人精品| 国产激情久久老熟女| 九草在线视频观看| 久久久久精品性色| 又大又黄又爽视频免费| 国产av一区二区精品久久| 午夜福利视频精品| 日韩一本色道免费dvd| 国产爽快片一区二区三区| 国产精品免费大片| 日韩熟女老妇一区二区性免费视频| 丝袜脚勾引网站| 婷婷色综合大香蕉| 亚洲丝袜综合中文字幕| 国产成人免费无遮挡视频| 纵有疾风起免费观看全集完整版| 三上悠亚av全集在线观看| 日本猛色少妇xxxxx猛交久久| 欧美激情极品国产一区二区三区 | 国产熟女午夜一区二区三区| 美女中出高潮动态图| 亚洲欧美日韩卡通动漫| 22中文网久久字幕| 久久精品国产a三级三级三级| 国产精品三级大全| 国产精品一区二区在线观看99| 狂野欧美激情性bbbbbb| 亚洲国产精品专区欧美| 亚洲国产精品999| 又黄又粗又硬又大视频| 一区在线观看完整版| 青春草视频在线免费观看| 中文字幕av电影在线播放| 欧美精品一区二区大全| 亚洲精品国产av蜜桃| 免费av中文字幕在线| 免费高清在线观看日韩| 国产亚洲一区二区精品| 五月玫瑰六月丁香| 国产欧美亚洲国产| 狂野欧美激情性bbbbbb| 99热网站在线观看| 亚洲,欧美,日韩| www日本在线高清视频| 看非洲黑人一级黄片| 成人午夜精彩视频在线观看| 午夜91福利影院| 亚洲精品,欧美精品| 亚洲第一av免费看| 国产1区2区3区精品| 国产精品 国内视频| 亚洲图色成人| 搡老乐熟女国产| 全区人妻精品视频| 如何舔出高潮| 一级毛片电影观看| 国产精品久久久久久久电影| 丝袜在线中文字幕| 免费av中文字幕在线| 最近中文字幕2019免费版| 黄网站色视频无遮挡免费观看| 国产一区亚洲一区在线观看| 黄色一级大片看看| 99热6这里只有精品| 18禁裸乳无遮挡动漫免费视频| 日本黄色日本黄色录像| 免费黄色在线免费观看| 亚洲成国产人片在线观看| 亚洲国产日韩一区二区| 国产 精品1| 国产成人91sexporn| 欧美激情国产日韩精品一区| 国产精品久久久久久精品电影小说| 五月伊人婷婷丁香| 少妇的逼水好多| 又黄又粗又硬又大视频| 精品国产乱码久久久久久小说| 久久久欧美国产精品| 嫩草影院入口| 国产黄频视频在线观看| 国产xxxxx性猛交| 热re99久久国产66热| 超色免费av| 亚洲国产看品久久| 亚洲婷婷狠狠爱综合网| 香蕉国产在线看| 亚洲国产看品久久| 精品国产国语对白av| 如日韩欧美国产精品一区二区三区| 午夜福利乱码中文字幕| 亚洲综合精品二区| 免费在线观看黄色视频的| 欧美老熟妇乱子伦牲交| 五月天丁香电影| 色婷婷av一区二区三区视频| 亚洲精品美女久久av网站| 亚洲五月色婷婷综合|