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

    輸水渡槽對河道行洪風險的影響分析

    2023-10-21 01:58:58劉開顏龔來紅謝亨旺劉博懿
    中國農(nóng)村水利水電 2023年10期
    關鍵詞:山溪行洪渡槽

    劉開顏,付 湘,龔來紅,謝亨旺,劉博懿

    (1. 武漢大學 水資源工程與調(diào)度全國重點實驗室,湖北 武漢 430072; 2. 江西省灌溉試驗中心站,江西 南昌 330201)

    0 引 言

    河流涉水工程可以滿足防洪、灌溉、發(fā)電、供水等方面效益,在地區(qū)經(jīng)濟社會中有重要地位。但是涉水工程的建設會占用河道原有的過水斷面,壅高水位,對河道行洪能力造成不利影響,可能導致洪水漫溢。而洪水沖擊下對工程本身造成破壞,可能導致工程潰決。兩種災害形式都將對當?shù)氐纳钌a(chǎn)和周邊生態(tài)環(huán)境造成不可估計的損失。因此,我們不僅需要對涉水工程對河道行洪能力的影響進行深入研究,還應該關注洪水作用下的工程安全風險。

    在涉水工程影響行洪能力的研究中,張大茹等[1]利用MIKE21 建立平面二維水流數(shù)學模型,考慮建筑物設計標準、橋梁尺寸、結(jié)構等影響因素,分析跨河工程對河道行洪能力的影響。何貞俊等[2]采用物理模型試驗不同水流條件、不同橋梁間距及碼頭群布置形式下,涉水工程疊加所引起的水位壅高值變化規(guī)律。郭立兵等[3]對涉水建筑物概化處理,依據(jù)圣維南方程基本原理,采用有限差分法離散求解方程組,建立了河道一維數(shù)值仿真模型,得出涉水建筑物在不同運行工況下對河道行洪產(chǎn)生不利影響。劉麗紅等[4]構建非結(jié)構化水動力學模型,將進洪退洪過程作為整體對濛洼蓄洪區(qū)洪水演進過程進行實時模擬,分析橋梁工程建設對蓄洪區(qū)行洪能力的影響。WANG 等[5]采用MIKE21 建立二維數(shù)值模型,從最大水位和最大流速兩方面定量分析橋墩對中國嘉陵江行洪的影響。HAN 等[6]針對瓶頸河段河道寬度小、斷面窄深的問題,采用水平二維數(shù)值模型計算河流平面形狀和斷面,模擬不同洪水條件下洪水水力因子和施工后階段的變化橋梁建設對河流防洪的影響。TOMASZ等[7]使用程序模擬洪水水面線,進而用淹沒范圍和淹沒深度評估橋梁對波蘭瓦爾塔河行洪的影響。PERUZZI 等[8]建立了包含復雜形狀水工結(jié)構的三維流體動力學模型,從水位、流量系數(shù)、流線、能量損失等方面探究水工結(jié)構對水流的影響。

    在洪水作用下的工程安全風險評估中,王仲玨等[9]基于水文、氣象、地質(zhì)、工程結(jié)構的隨機特征,采用蒙特卡羅法及一階二次矩法分析洪水作用下南水北調(diào)中線工程渡槽滑移失效的風險。殷丹等[10]基于糙率的不確定性,采用情景模擬法評估河道防洪能力以及堤段設計洪水標準下的河道行洪風險。琚烈紅等[11]針對海堤防洪問題,采用故障樹法進行風險識別,并運用直接積分法和蒙特卡洛方法評估海堤洪水災害風險。

    鑒于各類涉水工程具有不同的結(jié)構特點和功能特性,對于河道行洪能力的影響需要做針對性研究,并且綜合考慮河道行洪能力和工程安全,可以更加全面地評估涉水建筑物修建的合理性以及管理的必要性。本文提出針對修建于河道內(nèi)的跨水渡槽,考慮工程安全及阻水效果,評估其對河道行洪風險的影響。以江西省贛撫平原崗前渡槽為研究對象,從渡槽本身的結(jié)構穩(wěn)定性以及河道水面高度是否超越堤防高度兩方面,探究其對清豐山溪行洪風險的影響。以清豐山溪不同設計頻率的洪水為輸入,通過概化渡槽結(jié)構,依據(jù)洪水特性分析并計算渡槽所受到的洪水水平?jīng)_擊力、洪水豎向上托力、浮力、渡槽自身重力、渡槽內(nèi)水體重力以及風力,進而求解維護渡槽抗滑抗傾穩(wěn)定性所允許清豐山溪行洪的流量閾值。其次,依托一維水動力模型和基于拉丁超立方抽樣改進的蒙特卡羅法隨機模擬法,計算在有無渡槽的兩種情況下清豐山溪行洪風險,以期為崗前渡槽和清豐山溪的后續(xù)管理提供科學參考。

    1 研究區(qū)域概況和數(shù)據(jù)來源

    1.1 研究區(qū)域概況

    崗前渡槽位于江西省贛撫平原灌區(qū),是江西省最大的渡槽,槽身及槽身支撐段均為鋼筋混凝土結(jié)構,其運行六十余年后,逐漸出現(xiàn)混凝土老化剝蝕、槽身開裂漏水等安全問題,被定性為四類渡槽。

    清豐山溪位于江西省中部,發(fā)源于升華山,區(qū)域內(nèi)水系復雜,流經(jīng)崗前渡槽下接撫河故道,最終流入鄱陽湖。

    崗前渡槽跨越清豐山溪排洪道,二者具體的地理位置見圖1。清豐山溪存在河道淤堵嚴重,行洪不暢的問題,所以研究崗前渡槽對清豐山溪行洪能力的影響,進而權衡河道治理與工程安全問題具有重大意義,可以為灌區(qū)管理局的決策提供參考。

    圖1 清豐山溪與崗前渡槽地理位置示意圖Fig.1 Sketch of Gangqian aqueduct

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

    工程及地質(zhì)數(shù)據(jù)來自《江西省贛撫平原灌區(qū)“十四五”續(xù)建配套與現(xiàn)代化改造工程二期可研報告》、《江西省贛撫平原灌區(qū)崗前渡槽安全鑒定報告》和《贛撫平原初設報告》(CS3-2002.06),工程規(guī)范采用《灌溉與排水渠系建筑物設計規(guī)范》(SL482-2011)、《水利水電工程等級劃分及洪水標準(SL252-2017)》。水文數(shù)據(jù)來自清豐山溪的兩個控制站——吳石水文站(115°58'E, 28°21'N)和崗前渡槽下水位站(115°57'E,28°21'N)。

    2 研究方法

    河道行洪風險與河道水面線緊密聯(lián)系,而水面線的推求過程中存在諸多不確定性因子。將主要的不確定性因子納入考慮,構建河道行洪風險模型,通過蒙特卡羅方法隨機模擬水面線,抽取大量樣本后得到較為準確的河道行洪風險。

    2.1 影響河道行洪風險的不確定性因子的隨機模擬

    2.1.1 河流洪峰流量的不確定性

    超標準洪水是影響河道行洪風險的主要因素。對于某種工況,不同標準的洪水對河道行洪風險有不同的影響。根據(jù)《報告》建議,采用P-Ⅲ分布模擬清豐山溪的洪峰流量。

    P-Ⅲ分布表達式為:

    式中:Q為洪峰流量;m3/s,α、β、a0為皮爾遜Ⅲ型分布的形狀、尺度和位置參數(shù);Γ(α)為α的伽馬函數(shù)。

    3個參數(shù)α、β、a0的計算方法如下:

    式中:EX為總體均值;Cs為偏態(tài)系數(shù);Cv為變差系數(shù)。

    2.1.2 河道初值水位的不確定性

    根據(jù)上述河流洪峰流量的不確定性,運用清豐山溪的水位—流量關系曲線,可得到河流初值水位的不確定性。

    2.1.3 河道糙率的不確定性

    河道糙率作為反映河流阻力的綜合性系數(shù),是河道水力計算中重要而敏感的參數(shù)[12]。若缺乏河道糙率具體數(shù)值,可以根據(jù)天然河道糙率表,假設其服從正態(tài)分布[13]。

    2.1.4 河道兩岸堤防高程的不確定性

    堤防高程由于設計誤差及人類活動存在高低不平的情況,造成其不確定性。一般可以假設其符合正態(tài)分布,根據(jù)實測情況,計算均值和方差。

    2.2 渡槽抗滑抗傾穩(wěn)定性安全系數(shù)的確定

    根據(jù)《灌溉與排水渠系建筑物設計規(guī)范》(SL482-2011),渡槽抗滑穩(wěn)定性安全系數(shù)K1需滿足下式要求:

    式中:N表示作用于槽身的豎向力,N;P表示作用于槽身的水平力,N;fb表示支座的摩擦系數(shù)。

    渡槽抗傾覆穩(wěn)定性安全系數(shù)公式K2應滿足:

    式中:K2表示抗傾覆穩(wěn)定安全系數(shù);∑MV表示所有垂直力繞背水面支點轉(zhuǎn)動的傾覆力矩總和,kN·m;∑MP表示所有水平力繞背水面支點轉(zhuǎn)動的力矩總和,kN·m。

    渡槽受力圖如圖2所示。

    圖2 渡槽受力示意圖Fig.2 Sketch of the force on aqueduct

    其中橫向風壓力、渡槽自重、渡槽內(nèi)水重、靜水壓力、靜水浮托力、洪水橫向沖擊力、洪水豎向上托的計算公式如下:

    (1)作用于槽身的橫向風壓力。

    式中:WZ表示風荷載標準值,KN/m2;W0表示基本風壓,KN/m2;μS表示風載體型系數(shù);μZ表示風壓高度變化系數(shù);βZ表示風振系數(shù);S表示迎風面面積,m2。

    (2)渡槽自重。

    式中:m表示渡槽自身重量,kg。

    (3)渡槽內(nèi)水重。

    式中:Gw表示水重,N;B表示渡槽寬度,m;b表示渡槽槽身厚度,m;h表示渡槽內(nèi)水深,m。

    (4)靜水壓力。

    式中:Fs表示靜水壓力,N;;H表示清風山溪淹沒崗前渡槽水深,m;L表示渡槽長度,m。

    (5)靜水浮托力。

    式中:Ff表示靜水浮托力,N;VP表示渡槽淹沒部分的體積,m3。

    (6)洪水橫向沖擊力。MALAVASI 和GUADAGNIN[14]根據(jù)流水具有的動能平衡原理提出水流對橋梁的水平作用力計算公式,以此為參考計算渡槽承受的洪水沖擊力:

    式中:δ為渡槽阻水系數(shù);A為渡槽在垂直于水流方向的投影面積,m2;v為斷面平均流速,m/s。

    (7)洪水豎向托力。洪水不是穩(wěn)定的層流,存在一定的波動。洪水中任一點的瞬時流速均可視為由時均流速-v和脈動流速u'組成,時均流速為穩(wěn)定的層流,其方向平行于渡槽底面,不對渡槽產(chǎn)生豎向上托力。而由于波動引起的脈動流速u'方向是不確定的,將其分解為水平分量ux'和豎向分量uy',渡槽底部受到由uy'產(chǎn)生的上托力[15]。

    研究表明,梁體底面豎向流速水壓的分布呈直線衰減,在迎水面底部水壓為ρgΔh,在背水面底部水壓為0,具體見圖2。故渡槽底部中間的豎向上托力為:

    式中:Δh為渡槽迎水面的壅水高度,m。

    D’ Aubuisson 公式是目前工程中常用的壅水高度經(jīng)驗計算公式[16]:

    式中:μ為測收縮系數(shù);hb表示渡槽后水深,m;Δβ表示渡槽排架總寬度,m。

    2.3 河道行洪風險計算的蒙特卡洛方法

    本文研究洪水使渡槽失穩(wěn)以及洪水超過堤防高程的河道行洪風險,風險率可用下式表示:

    式中:Z(x)表示行洪水位高程,m;M(x)表示堤防高程,m;K1由式(3)求得;K2由式(4)求得。

    Z(x)由圣維南方程組求解,不考慮恒定性的影響,忽略時間因素的作用,進而得到如下形式的恒定非均勻流水面線基本方程[17]。

    式中:α為流速水頭系數(shù);ξ為局部損失系數(shù);ω為流量模數(shù);n為糙率;R為水力半徑。

    河道行洪風險采用常規(guī)蒙特卡羅方法與基于拉丁超立方體的蒙特卡羅方法進行計算。常規(guī)蒙特卡羅方法通過計算機模擬產(chǎn)生處于(0,1)區(qū)間的隨機數(shù),在已知概率分布的情況下,通過反變換法得到各變量的隨機數(shù),代入到相應的風險率計算函數(shù)中,從而求得風險率的無偏估計量。該方法不受極限狀態(tài)函數(shù)復雜程度的影響,也不需要將狀態(tài)函數(shù)線性化和隨機變量等效正態(tài)化,當抽樣數(shù)量足夠多時,可以達到非常高的計算精度[18],是預測和估算可靠度常用的方法之一,但也存在計算次數(shù)多耗時長的缺點[19]。

    基于拉丁超立方體的蒙特卡羅方法對常規(guī)蒙特卡羅法的隨機抽樣技術進行改,具有更高的運算效率和計算精度。拉丁超立方抽樣是一種多維分層抽樣法,由MCKAY 等人[20]于1979年提出,其本質(zhì)是通過產(chǎn)生分布更加均勻的樣本,在抽取更少樣本的同時獲得更高的精度[21],該方法下每個子區(qū)間僅產(chǎn)生一個隨機數(shù),避免重復的隨機數(shù)降低計算效率。公式如下:

    式中:i= 1,2,…,N;U為[0,1]區(qū)間內(nèi)均勻分布的隨機數(shù);Ui是第i個子區(qū)間中的隨機數(shù)。

    3 結(jié)果與分析

    3.1 清豐山溪洪峰流量隨機模擬

    根據(jù)吳石站1961-2018 年共58 年流量資料系列,采用公式(1)~(2)計算清豐山溪洪峰流量的P-Ⅲ型分布,統(tǒng)計參數(shù)為:EX=490 m3/s,CV=0.63,CS=3.0CV,累積密度函數(shù)(CDF)見圖3。

    圖3 清豐山溪洪峰流量概率分布圖Fig.3 PDF of peak discharge of Qingfeng Stream

    圖3表明清豐山溪洪峰流量理論分布概率值和經(jīng)驗分布概率值相近,進一步通過SPSS軟件分析得到二者R2為0.994,可知P-Ⅲ分布對清豐山溪洪峰流量的擬合較好。

    3.2 渡槽抗滑抗傾穩(wěn)定安全系數(shù)計算

    根據(jù)《報告》中“現(xiàn)狀崗前渡槽為四類渡槽,當其運行水深不超過2.0 m時,渡槽結(jié)構滿足安全運行要求”,設定渡槽的4種運行工況,即控制渡槽內(nèi)水位h分別為0.5、1.0、1.5、2.0 m。運用Matlab R2020b 編程得到偽隨機數(shù)后采用反變換法得到符合P-Ⅲ分布的清豐山溪洪峰流量,根據(jù)公式(3)~(12)計算渡槽的抗滑穩(wěn)定系數(shù)K1和抗傾穩(wěn)定系數(shù)K2,計算結(jié)果見圖4。

    圖4 不同槽內(nèi)水位下渡槽抗滑抗傾穩(wěn)定系數(shù)Fig.4 Coefficient of anti-sliding and anti-overturn stability of aqueduct under different h

    由圖4(a)可知,當h分別為0.5、1.0、1.5、2.0 m 時,符合渡槽的抗滑穩(wěn)定性系數(shù)規(guī)范值[K1]=1.05 的流量閾值依次為759.9、865.6、972.1以及1 079.4 m3/s。

    由圖4(b)可知,當h分別為0.5、1.0、1.5、2.0 m 時,符合渡槽的抗傾穩(wěn)定性系規(guī)范值[K2]=1.1 的流量閾值依次為868.5、1 024.3、1 192.8以及1 361.7 m3/s。

    相同清豐山溪洪水流量的情況下,隨著渡槽內(nèi)水位升高,抗滑穩(wěn)定系數(shù)和抗傾穩(wěn)定系數(shù)均變大,渡槽具有更好的穩(wěn)定性,這說明在洪水來臨時,應適當抬高槽內(nèi)水位保證渡槽結(jié)構的穩(wěn)定性。

    綜合抗滑抗傾穩(wěn)定性系數(shù)對應的流量閾值,若上游發(fā)生了20 年一遇的洪水,洪峰流量達到1 200 m3/s,在渡槽內(nèi)水位分別為0.5,1.0,1.5,2.0 m 的4 種工況下,清豐山溪(崗前渡槽段)的行洪能力分別下降36.68%、27.87%、19.00%、10.05%;

    3.3 河道行洪風險計算

    根據(jù)清豐山溪兩岸堤防高程實測值得到均值和方差,將其看作服從N(26.1,12)。河床地質(zhì)成分主要有粉質(zhì)黏土、細砂及淤泥質(zhì)粉質(zhì)黏土,因為缺乏從河道糙率的實測值,故依據(jù)《天然河道糙率表》,以及附近水域相關文獻所使用的糙率值[22,23],假定其服從N(0.022,0.0012)。運用Matlab R2020b 編程得到堤防高程和糙率的隨機模擬值。

    根據(jù)公式(13)~(16),基于水面線計算結(jié)果,采用常規(guī)MC法和LHS-MC 法,計算在沒有渡槽和有渡槽且槽內(nèi)水位控制在0.5 m 兩種情況下,不同模擬次數(shù)所對應的河道行洪風險,具體見圖5。

    圖5 常規(guī)MC法和LHS-MC法下不同模擬次數(shù)對應的河道行洪風險Fig.5 Flood risk corresponding to different simulation times under MC and LHS-MC

    結(jié)合圖5和模擬數(shù)值,MC法在50 000次左右時可以達到誤差小于0.001,LHS-MC 法在38 000 次左右可以達到誤差小于0.001,因此,LHS-MC 法比常規(guī)MC 法所需的模擬次數(shù)更少,風險收斂更快。

    無渡槽時,河道行洪風險收斂后的穩(wěn)定值約為0.007;有渡槽時,河道行洪風險收斂后的穩(wěn)定值約為0.34。說明渡槽對河道行洪能力有負面影響,致使河道防洪標準低于設計值。

    運用相同方法計算有渡槽情況下,渡槽內(nèi)水位分別為1.0,1.5及2.0 m 時的河道行洪風險,結(jié)果依次為0.045、0.016、0.012,這說明提高渡槽內(nèi)水位可以大幅降低河道行洪風險。

    控制渡槽內(nèi)水位為2.0 m,計算不同清豐山溪流量下的河道行洪風險,具體見圖6。

    圖6 當h=2.0 m時,清豐山溪流量和行洪風險關系Fig.6 h=2.0 m, relationship between discharge and flood risk of Qingfeng Stream

    由圖6 可知,盡管控制渡槽內(nèi)水位為其所能承擔的最高水位2.0 m,當清豐山溪流量大于1 600 m3/s時,河道行洪風險仍會顯著上升,此時可以考慮報廢渡槽以保證清豐山溪的行洪能力,盡量避免河道兩旁的農(nóng)田和房屋承受洪澇災害損失。

    綜合研究所得結(jié)果,渡槽本身結(jié)構穩(wěn)定性的要求會降低清豐山溪的可行洪流量,適當抬升渡槽內(nèi)水位可以減少這種不利影響。此外,渡槽使得河道過水斷面減小,水面線升高,造成河道行洪風險上升。所以當清豐山溪上游發(fā)生超設計流量1 200 m3/s 的洪水時,應盡可能抬高渡槽內(nèi)水位,降低河道行洪風險;當洪水流量超過1 600 m3/s時應考慮報廢渡槽。

    5 結(jié) 論

    本文提出了輸水渡槽對河道行洪能力的影響分析方法,以江西省贛撫平原灌區(qū)崗前渡槽和清豐山溪為例,考慮洪峰流量、初值水位、河道糙率以及堤防高程的不確定性,計算河道行洪風險,結(jié)果如下。

    (1)基于洪水力學特性,從維護渡槽抗滑抗傾穩(wěn)定性角度出發(fā),反推清豐山溪行洪流量閾值。當渡槽內(nèi)水位分別為0.5、1.0、1.5、2.0 m 時,行洪流量閾值依次為759.9、865.6、972.1 和1 079.4 m3/s。相比于清豐山溪設計流量1 200 m3/s,河道行洪能力分別下降36.68%、27.87%、19.00%、10.05%。

    (2)沒有渡槽的情況下,河道行洪風險約0.007;有渡槽的情況下,考慮渡槽壅水高度影響,在(1)所述4種工況下,河道行洪風險依次為0.34、0.045、0.016、0.012,說明渡槽對河道行洪能力有負面影響,致使河道防洪標準低于設計值,但提高渡槽內(nèi)水位可以有效削減負面影響。

    (3)當清豐山溪流量超過1 200 m3/s 時,應適當提升渡槽內(nèi)水位降低行洪風險;當清豐山溪流量超過1 600 m3/s 時,即使渡槽內(nèi)水位處于安全范圍內(nèi)可承受的最高值,行洪風險仍隨著河道流量增大顯著增加,此時應考慮報廢渡槽。

    (4)基于拉丁超立方體的蒙特卡羅隨機模擬法比常規(guī)蒙特卡羅隨機模擬法具有更好的收斂性,但兩種方法達到同一精度的模擬次數(shù)差別不大,對于渡槽對河道行洪風險影響評估,兩種隨機模擬方法均適用。

    猜你喜歡
    山溪行洪渡槽
    山溪
    幼兒100(2023年42期)2023-11-21 09:36:38
    淺談水利工程渡槽基礎承臺施工
    傍山道路泥石流渡槽設計
    中華建設(2020年5期)2020-07-24 08:55:46
    滹沱河(南水北調(diào)—機場路)南岸高標準行洪區(qū)行洪通道布局存在問題及對策建議
    河北水利(2020年6期)2020-07-21 01:58:52
    山溪
    淮河正陽關至渦河口段行洪區(qū)單獨運用時對淮河干流水位的影響研究
    治淮(2020年1期)2020-04-11 01:32:54
    袁曉峰攝影作品欣賞
    河北水利(2017年6期)2017-09-05 11:04:13
    山溪
    GPS控制網(wǎng)在大型河流行洪能力分析中的應用
    基于ANSYS的渡槽三維有限元分析
    国产男女超爽视频在线观看| 国产国语露脸激情在线看| 色婷婷av一区二区三区视频| 啦啦啦视频在线资源免费观看| 国产乱来视频区| 飞空精品影院首页| 免费少妇av软件| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | videossex国产| 久久人人爽人人片av| 精品亚洲成国产av| 99视频精品全部免费 在线| 99热国产这里只有精品6| 久久综合国产亚洲精品| 999精品在线视频| 日韩在线高清观看一区二区三区| 母亲3免费完整高清在线观看 | 满18在线观看网站| 日韩av免费高清视频| 国语对白做爰xxxⅹ性视频网站| 欧美日韩亚洲高清精品| 黄色一级大片看看| 一区二区av电影网| 成人无遮挡网站| 99热这里只有精品一区| 久久鲁丝午夜福利片| 大香蕉久久成人网| 欧美日韩视频高清一区二区三区二| 久久久久国产精品人妻一区二区| 国产不卡av网站在线观看| 日韩一本色道免费dvd| 午夜激情久久久久久久| 欧美精品亚洲一区二区| 欧美97在线视频| 欧美最新免费一区二区三区| 日本vs欧美在线观看视频| 美女内射精品一级片tv| 在线观看人妻少妇| 在线看a的网站| 各种免费的搞黄视频| 成人无遮挡网站| 精品亚洲乱码少妇综合久久| 日本-黄色视频高清免费观看| 国产精品秋霞免费鲁丝片| 亚洲精品av麻豆狂野| 国产精品久久久久久精品古装| 在线播放无遮挡| 男女啪啪激烈高潮av片| 亚洲精品国产色婷婷电影| 国产精品一区www在线观看| 国产色婷婷99| 男女边摸边吃奶| 色94色欧美一区二区| 亚洲av成人精品一二三区| 久久97久久精品| 三上悠亚av全集在线观看| 日韩,欧美,国产一区二区三区| 欧美精品一区二区免费开放| 日日爽夜夜爽网站| 高清在线视频一区二区三区| 久久人妻熟女aⅴ| 国产精品国产三级国产av玫瑰| 一区二区三区四区激情视频| 18禁动态无遮挡网站| 在线精品无人区一区二区三| 亚洲欧洲国产日韩| 国产在视频线精品| 成年女人在线观看亚洲视频| 亚洲精品亚洲一区二区| 欧美日韩视频高清一区二区三区二| 热re99久久精品国产66热6| 国产成人一区二区在线| 成人黄色视频免费在线看| 成人18禁高潮啪啪吃奶动态图 | 男女免费视频国产| 欧美另类一区| 久热久热在线精品观看| .国产精品久久| 国产精品秋霞免费鲁丝片| 欧美日韩国产mv在线观看视频| 久久久久久久久久久丰满| 性高湖久久久久久久久免费观看| 亚洲精品乱久久久久久| 免费黄色在线免费观看| 国产欧美亚洲国产| 中国三级夫妇交换| 日本欧美国产在线视频| 午夜免费男女啪啪视频观看| 日本午夜av视频| 亚洲精品av麻豆狂野| 蜜桃久久精品国产亚洲av| 黑丝袜美女国产一区| 一级毛片 在线播放| 亚洲综合精品二区| 国产不卡av网站在线观看| 国产午夜精品久久久久久一区二区三区| 一级二级三级毛片免费看| 成人二区视频| 国产精品.久久久| 亚洲国产av影院在线观看| 最后的刺客免费高清国语| 91成人精品电影| 日本-黄色视频高清免费观看| 日韩av免费高清视频| 亚洲五月色婷婷综合| 91精品一卡2卡3卡4卡| av在线观看视频网站免费| 国产精品女同一区二区软件| 久久精品国产a三级三级三级| 你懂的网址亚洲精品在线观看| 国产欧美另类精品又又久久亚洲欧美| 午夜激情福利司机影院| .国产精品久久| 国产欧美另类精品又又久久亚洲欧美| 最近2019中文字幕mv第一页| 久久久久久久久久成人| 久久久久久伊人网av| 亚洲欧美清纯卡通| 午夜视频国产福利| 99热国产这里只有精品6| 一级毛片我不卡| 97在线人人人人妻| 亚洲精品日韩在线中文字幕| 少妇猛男粗大的猛烈进出视频| 精品人妻一区二区三区麻豆| 久久热精品热| 黄片播放在线免费| 久久精品国产亚洲av天美| 你懂的网址亚洲精品在线观看| 少妇的逼好多水| a级毛色黄片| 国产在线视频一区二区| 亚洲精品色激情综合| 狂野欧美激情性bbbbbb| videosex国产| 日韩av免费高清视频| 建设人人有责人人尽责人人享有的| 国产在线视频一区二区| 在线观看三级黄色| 亚洲精品中文字幕在线视频| 久久精品人人爽人人爽视色| 国产男人的电影天堂91| 蜜桃国产av成人99| 成年美女黄网站色视频大全免费 | 亚洲中文av在线| 久久久久久久久大av| 国内精品宾馆在线| 精品少妇黑人巨大在线播放| 亚洲精品乱久久久久久| 成人二区视频| 九九在线视频观看精品| 制服人妻中文乱码| 欧美成人午夜免费资源| 亚洲精品乱码久久久久久按摩| 日本爱情动作片www.在线观看| 精品亚洲乱码少妇综合久久| 天天操日日干夜夜撸| av天堂久久9| 丝袜在线中文字幕| 美女大奶头黄色视频| 国产综合精华液| 亚洲图色成人| 只有这里有精品99| 国产精品熟女久久久久浪| 蜜桃在线观看..| 少妇丰满av| 免费观看的影片在线观看| 日韩一区二区视频免费看| 日本与韩国留学比较| 久久韩国三级中文字幕| 97超碰精品成人国产| 最新中文字幕久久久久| 亚洲精品日本国产第一区| 91精品一卡2卡3卡4卡| 国产伦精品一区二区三区视频9| 曰老女人黄片| 国产成人精品福利久久| 一边亲一边摸免费视频| 九九爱精品视频在线观看| 大话2 男鬼变身卡| 男男h啪啪无遮挡| 在线观看一区二区三区激情| 69精品国产乱码久久久| 2021少妇久久久久久久久久久| 赤兔流量卡办理| 久久精品国产自在天天线| 在线 av 中文字幕| 国产精品.久久久| 日韩av免费高清视频| 欧美 日韩 精品 国产| 在线天堂最新版资源| 久久久精品94久久精品| 国产日韩欧美亚洲二区| 日本av免费视频播放| 一本久久精品| 考比视频在线观看| 亚洲一级一片aⅴ在线观看| 久久久久国产精品人妻一区二区| 肉色欧美久久久久久久蜜桃| 欧美丝袜亚洲另类| 人成视频在线观看免费观看| 日本av手机在线免费观看| 丝袜喷水一区| 国产精品偷伦视频观看了| 99久久综合免费| 国产成人精品一,二区| 亚洲美女黄色视频免费看| 一本久久精品| 亚洲中文av在线| av又黄又爽大尺度在线免费看| 春色校园在线视频观看| 欧美亚洲日本最大视频资源| 亚洲欧洲日产国产| 国语对白做爰xxxⅹ性视频网站| 精品亚洲成a人片在线观看| 国产男女超爽视频在线观看| 日本wwww免费看| 高清不卡的av网站| 香蕉精品网在线| 国产免费一级a男人的天堂| 一区二区三区四区激情视频| 91精品一卡2卡3卡4卡| 色吧在线观看| 午夜精品国产一区二区电影| 这个男人来自地球电影免费观看 | 人人妻人人澡人人爽人人夜夜| 日本猛色少妇xxxxx猛交久久| 国产精品蜜桃在线观看| 一区二区三区精品91| 高清欧美精品videossex| 亚洲在久久综合| 国产亚洲欧美精品永久| 新久久久久国产一级毛片| 久久99蜜桃精品久久| 日韩免费高清中文字幕av| 国产毛片在线视频| av一本久久久久| 国产在线一区二区三区精| 超色免费av| 免费观看的影片在线观看| 久久精品国产亚洲av涩爱| 在线观看三级黄色| 国产亚洲av片在线观看秒播厂| 超色免费av| 一本一本综合久久| 久久国产精品大桥未久av| 中文字幕最新亚洲高清| 99久久精品一区二区三区| 国产精品秋霞免费鲁丝片| 国产日韩欧美视频二区| 人成视频在线观看免费观看| 久久国内精品自在自线图片| 99精国产麻豆久久婷婷| 少妇猛男粗大的猛烈进出视频| 免费看不卡的av| 爱豆传媒免费全集在线观看| 精品人妻一区二区三区麻豆| 搡女人真爽免费视频火全软件| 97精品久久久久久久久久精品| 99视频精品全部免费 在线| 精品少妇久久久久久888优播| 久久这里有精品视频免费| 日日撸夜夜添| 久久精品国产a三级三级三级| 91国产中文字幕| 曰老女人黄片| 熟女电影av网| 午夜91福利影院| 国产不卡av网站在线观看| 黑人高潮一二区| 高清毛片免费看| 91aial.com中文字幕在线观看| a级毛片黄视频| 免费av不卡在线播放| 观看美女的网站| 亚洲精品久久成人aⅴ小说 | 久久精品国产鲁丝片午夜精品| av在线观看视频网站免费| 亚洲精品日本国产第一区| 夫妻性生交免费视频一级片| 国产精品女同一区二区软件| 欧美激情 高清一区二区三区| 大码成人一级视频| 国产亚洲精品第一综合不卡 | 亚洲高清免费不卡视频| 人妻制服诱惑在线中文字幕| 久久人人爽av亚洲精品天堂| av线在线观看网站| 日本免费在线观看一区| 国产高清国产精品国产三级| 国产精品久久久久久精品古装| 黄色一级大片看看| 精品久久国产蜜桃| 嘟嘟电影网在线观看| 夜夜看夜夜爽夜夜摸| 黄片播放在线免费| 久久久a久久爽久久v久久| 母亲3免费完整高清在线观看 | 少妇精品久久久久久久| 亚洲欧洲国产日韩| 欧美成人精品欧美一级黄| 久久免费观看电影| 国产老妇伦熟女老妇高清| 亚洲av二区三区四区| 免费大片18禁| 亚洲,欧美,日韩| 欧美另类一区| 少妇的逼水好多| 大片电影免费在线观看免费| 黄色欧美视频在线观看| 久久久久视频综合| 老女人水多毛片| 国产极品粉嫩免费观看在线 | freevideosex欧美| 久久精品夜色国产| 视频在线观看一区二区三区| 人妻少妇偷人精品九色| 国产免费一级a男人的天堂| 日韩免费高清中文字幕av| 亚洲av二区三区四区| 一区二区三区四区激情视频| 欧美激情极品国产一区二区三区 | 国产精品蜜桃在线观看| 国产一区二区在线观看日韩| 亚洲av不卡在线观看| 免费人妻精品一区二区三区视频| 亚洲精品成人av观看孕妇| 边亲边吃奶的免费视频| a 毛片基地| 天堂中文最新版在线下载| 精品卡一卡二卡四卡免费| 伦理电影大哥的女人| 综合色丁香网| 精品人妻熟女毛片av久久网站| 欧美性感艳星| 欧美+日韩+精品| 亚洲av电影在线观看一区二区三区| 久久婷婷青草| 日韩三级伦理在线观看| 日韩视频在线欧美| 精品久久国产蜜桃| 少妇的逼水好多| 多毛熟女@视频| 一级毛片aaaaaa免费看小| 精品午夜福利在线看| 新久久久久国产一级毛片| 亚洲国产av新网站| av在线app专区| 国产成人91sexporn| 亚洲人成网站在线播| 制服诱惑二区| 丝袜在线中文字幕| 久久久久精品久久久久真实原创| 丝瓜视频免费看黄片| 免费播放大片免费观看视频在线观看| 久久精品国产亚洲av天美| 视频区图区小说| 99热这里只有精品一区| 日本vs欧美在线观看视频| 丝瓜视频免费看黄片| 少妇人妻 视频| 久久久久久人妻| 国产 一区精品| 久久久精品区二区三区| 国产一区二区三区av在线| 三级国产精品片| 国模一区二区三区四区视频| 久久女婷五月综合色啪小说| 一二三四中文在线观看免费高清| 黄色配什么色好看| 一区二区三区四区激情视频| 黄色一级大片看看| 久久国产亚洲av麻豆专区| 欧美激情极品国产一区二区三区 | 免费黄网站久久成人精品| 最新中文字幕久久久久| 久热久热在线精品观看| 亚洲美女视频黄频| 亚洲少妇的诱惑av| 97超视频在线观看视频| 黄色欧美视频在线观看| a级毛片免费高清观看在线播放| av专区在线播放| 黑人猛操日本美女一级片| 少妇被粗大的猛进出69影院 | 王馨瑶露胸无遮挡在线观看| 看十八女毛片水多多多| 精品国产国语对白av| 波野结衣二区三区在线| av电影中文网址| 大陆偷拍与自拍| 春色校园在线视频观看| 亚洲av欧美aⅴ国产| 午夜激情av网站| 精品人妻偷拍中文字幕| 秋霞在线观看毛片| 免费观看a级毛片全部| 特大巨黑吊av在线直播| 色婷婷av一区二区三区视频| 国产色婷婷99| av天堂久久9| 久久久久久久久久成人| 国产精品久久久久久av不卡| 国产熟女欧美一区二区| 黄色视频在线播放观看不卡| 久久久午夜欧美精品| 久久精品国产亚洲网站| 美女主播在线视频| 十八禁高潮呻吟视频| 99热网站在线观看| 精品久久久噜噜| 国产精品三级大全| 亚洲,欧美,日韩| 亚洲欧美日韩卡通动漫| 熟女电影av网| 水蜜桃什么品种好| 国产精品一国产av| 亚洲精品一区蜜桃| 色吧在线观看| 久热这里只有精品99| 青春草亚洲视频在线观看| 成人午夜精彩视频在线观看| 人人妻人人添人人爽欧美一区卜| 国产成人精品福利久久| 夜夜骑夜夜射夜夜干| 交换朋友夫妻互换小说| 日本黄色片子视频| 日韩伦理黄色片| av在线app专区| 欧美精品国产亚洲| 热99久久久久精品小说推荐| 91在线精品国自产拍蜜月| 伦理电影大哥的女人| 人妻人人澡人人爽人人| 久久99一区二区三区| 水蜜桃什么品种好| 成人黄色视频免费在线看| 丝瓜视频免费看黄片| 爱豆传媒免费全集在线观看| 亚洲人成网站在线观看播放| 一区二区三区四区激情视频| 国产精品久久久久久久电影| 人妻少妇偷人精品九色| 午夜av观看不卡| 18禁观看日本| 亚洲精品乱久久久久久| 国产精品人妻久久久久久| 久久青草综合色| 夫妻性生交免费视频一级片| 亚洲欧美中文字幕日韩二区| 国产极品天堂在线| 自线自在国产av| 久久久久视频综合| 在线播放无遮挡| 啦啦啦视频在线资源免费观看| 99视频精品全部免费 在线| 美女cb高潮喷水在线观看| 免费大片黄手机在线观看| 大片免费播放器 马上看| 全区人妻精品视频| 亚洲欧美一区二区三区国产| 国产精品一国产av| 国产亚洲av片在线观看秒播厂| 日本猛色少妇xxxxx猛交久久| 国产毛片在线视频| 午夜91福利影院| 91久久精品电影网| 国产精品久久久久久久久免| 久久久久久伊人网av| 精品一品国产午夜福利视频| 免费大片18禁| 99视频精品全部免费 在线| 黄色视频在线播放观看不卡| 婷婷色综合大香蕉| 精品亚洲成国产av| 国产精品一区二区在线不卡| .国产精品久久| 国产精品久久久久久精品电影小说| 国产极品天堂在线| 乱人伦中国视频| 一边亲一边摸免费视频| 亚洲性久久影院| 一区二区三区免费毛片| 久久久亚洲精品成人影院| 一级黄片播放器| 天天躁夜夜躁狠狠久久av| 精品久久久精品久久久| 久久热精品热| 亚洲精品456在线播放app| 大话2 男鬼变身卡| 一本一本综合久久| 如日韩欧美国产精品一区二区三区 | 久热久热在线精品观看| 亚洲精品色激情综合| 国产精品不卡视频一区二区| 热re99久久国产66热| 人成视频在线观看免费观看| 欧美日韩亚洲高清精品| 啦啦啦在线观看免费高清www| 在线播放无遮挡| 王馨瑶露胸无遮挡在线观看| 一级,二级,三级黄色视频| 久久久久久久国产电影| 自拍欧美九色日韩亚洲蝌蚪91| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 最新中文字幕久久久久| 成年女人在线观看亚洲视频| 久久精品久久精品一区二区三区| 肉色欧美久久久久久久蜜桃| 日本黄色片子视频| 一本—道久久a久久精品蜜桃钙片| 成人午夜精彩视频在线观看| 国产成人a∨麻豆精品| av网站免费在线观看视频| 久久久久久久大尺度免费视频| 亚洲av福利一区| 午夜福利,免费看| 国语对白做爰xxxⅹ性视频网站| 18在线观看网站| 纯流量卡能插随身wifi吗| 在线观看三级黄色| 成年人免费黄色播放视频| 免费观看的影片在线观看| 在线观看免费日韩欧美大片 | 中文字幕最新亚洲高清| 亚洲成人av在线免费| 韩国av在线不卡| 国产成人91sexporn| 国产成人午夜福利电影在线观看| 99九九在线精品视频| 18在线观看网站| av在线app专区| 亚洲av电影在线观看一区二区三区| 99精国产麻豆久久婷婷| 国产精品一二三区在线看| 最近最新中文字幕免费大全7| 日韩在线高清观看一区二区三区| 一区二区三区乱码不卡18| 亚洲激情五月婷婷啪啪| av黄色大香蕉| 久久ye,这里只有精品| 最近中文字幕2019免费版| 国产视频内射| 亚洲精品色激情综合| 国产探花极品一区二区| a级毛片黄视频| 国产极品粉嫩免费观看在线 | 五月玫瑰六月丁香| 人妻制服诱惑在线中文字幕| 人妻少妇偷人精品九色| 亚洲国产av影院在线观看| 蜜桃在线观看..| 日韩免费高清中文字幕av| 91aial.com中文字幕在线观看| 久久久久久久亚洲中文字幕| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品一区二区三区在线| 亚洲精品456在线播放app| 一级片'在线观看视频| 岛国毛片在线播放| 久久久久国产精品人妻一区二区| 99热网站在线观看| 高清视频免费观看一区二区| 国产亚洲精品久久久com| 哪个播放器可以免费观看大片| 日韩电影二区| 亚洲高清免费不卡视频| av在线观看视频网站免费| 又黄又爽又刺激的免费视频.| 国产精品女同一区二区软件| 亚洲av日韩在线播放| 久久久久久久国产电影| 男女国产视频网站| 欧美日韩一区二区视频在线观看视频在线| 亚洲成人一二三区av| 少妇的逼好多水| 蜜桃在线观看..| 亚洲综合精品二区| 国产熟女午夜一区二区三区 | 亚洲欧美日韩卡通动漫| 男女啪啪激烈高潮av片| 国产亚洲av片在线观看秒播厂| 人人妻人人添人人爽欧美一区卜| 交换朋友夫妻互换小说| 国产乱人偷精品视频| 成人18禁高潮啪啪吃奶动态图 | 午夜91福利影院| 国产精品无大码| 日韩精品免费视频一区二区三区 | 老司机亚洲免费影院| 又大又黄又爽视频免费| a级片在线免费高清观看视频| av电影中文网址| 欧美日韩精品成人综合77777| 少妇被粗大的猛进出69影院 | 啦啦啦啦在线视频资源| 少妇的逼好多水| 最近中文字幕2019免费版| 美女视频免费永久观看网站| 一级毛片 在线播放| 亚洲精品美女久久av网站| 最近的中文字幕免费完整| 成人国语在线视频| 99re6热这里在线精品视频| 精品人妻熟女av久视频| 亚洲国产精品专区欧美| 麻豆成人av视频| 水蜜桃什么品种好| 国产精品国产三级国产av玫瑰| 国产在线视频一区二区| 亚洲欧美日韩另类电影网站| 日韩中文字幕视频在线看片| 久久久国产精品麻豆| 另类精品久久| a级毛片免费高清观看在线播放|