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

    非飽和流對邊坡內(nèi)水文過程及穩(wěn)定性的影響分析

    2021-07-01 02:52:20王銳孫杰趙紅軍孔俊
    河北工業(yè)科技 2021年3期
    關鍵詞:邊坡穩(wěn)定性滑坡

    王銳 孫杰 趙紅軍 孔俊

    摘要:為了探究滲流水文過程對非飽和土質(zhì)邊坡穩(wěn)定性的影響機制,建立了能夠模擬地下水流和分析邊坡穩(wěn)定性的耦合模型?;?種表述邊坡水文過程的HSB模型開展建模對比分析,揭示邊坡內(nèi)非飽和滲流對水文過程的影響機制;通過耦合水文模型與邊坡穩(wěn)定極限平衡模型,比較了淺層邊坡在不同水力滲流條件下,因鄰近水位升降所造成穩(wěn)定性系數(shù)變化的差異,探討非飽和側向滲流過程對淺層邊坡穩(wěn)定性的影響機制。分析結果表明,土體非飽和區(qū)的儲水、釋水能力和滲透能力影響著非飽和區(qū)域流的滲流流速大小;在淺層邊坡中,非飽和側向滲流作用可以促進水體在邊坡內(nèi)部的遷移而引起邊坡內(nèi)部水位的快速抬高,致使邊坡穩(wěn)定性降低,并導致邊坡失穩(wěn)時間提前。研究結果可以為臨近地域水位變動區(qū)邊坡工程實施提供相關技術參考。

    關鍵詞:水文地質(zhì)學;滑坡;非飽和流;邊坡穩(wěn)定性;流域水文模型

    中圖分類號:P641文獻標識碼:ADOI: 10.7535/hbgykj.2021yx03012

    Abstract: In order to explore the influence mechanism of seepage hydrological process on the stability of unsaturated soil slope, a coupling model was established to simulate groundwater flow and analyze slope stability. Based on three Hillslope-storage Boussinesq models describing the hydrological process of the slope, the modeling comparative analysis was carried out to reveal the influence mechanism of unsaturated seepage in the slope on the hydrological process. By coupling the hydrological model with the slope stability limit equilibrium model, the variations of stability coefficients of shallow slope caused by the rise and fall of adjacent water levels under different hydraulic seepage conditions were compared, and the influence mechanism of unsaturated lateral flow process on the stability of shallow slope was discussed. The results show that the water storage and release capacity and the permeability of soil in the unsaturated area affect the seepage velocity in unsaturated area. In the shallow slope, unsaturated lateral seepage action can promote the water movement in the slope, which will lead to the rapid rise of the internal water level of the slope and the reduction of the slope stability, resulting in the slope instability time. The research results can provide relevant technical reference for the implementation of slope engineering in the adjacent area with water level fluctuation.

    Keywords:hydrogeology; landslide; unsaturated flow; slope stability; watershed hydrological model

    中國幅員遼闊、人口眾多,是世界上受到地質(zhì)災害影響較為嚴重的國家之一。隨著人類生產(chǎn)生活需求的不斷增長,開采資源的程度逐年遞增,改造自然的力度不斷加大,致使地質(zhì)災害的發(fā)生頻率逐年增加?;率且环N受多因素影響的地貌改造過程[1],也是最常見的一種地質(zhì)災害類型[2]。據(jù)統(tǒng)計,全球每年有數(shù)千人因滑坡事故而死亡或受傷[3],給生產(chǎn)生活帶來不可估量的損失。

    國內(nèi)外研究成果表明邊坡內(nèi)部水文過程對邊坡穩(wěn)定性具有不可忽視的影響[4-5],若在邊坡失穩(wěn)預警工作中忽略非飽和流的作用,邊坡內(nèi)部水位的計算將出現(xiàn)低估問題,進而影響邊坡穩(wěn)定安全系數(shù),使得邊坡預警系統(tǒng)難以準確預測滑坡發(fā)生時間和滑塌面積的大小。選擇合適的流域水文模型有利于更加高效地描述淺層邊坡內(nèi)部水文過程。相較于二維/三維的變飽和模型(如Richards方程),即使面對小規(guī)模問題也需要進行大量的計算[6],Hillslope-storage Boussinesq(簡稱HSB)模型由于擁有更高的計算效率和魯棒性而受到較多關注。目前滑坡的相關研究都側重于分析降雨情況下的邊坡穩(wěn)定性[7-10],對以側向流為主要誘因的淺層滑坡研究較少。

    本文創(chuàng)新性地引入HSB水文模型來模擬邊坡內(nèi)部的水位變化過程,引入3種對非飽和流處理方式不同的水文模型:Troch-HSB模型[11]、Hilberts-HSB模型[12]和Kong-HSB模型[13]。其中Troch-HSB模型為考慮飽和流影響的模型;Hilberts-HSB模型進一步考慮了非飽和區(qū)的垂向流作用;Kong-HSB模型則在Hilberts-HSB模型的基礎上又考慮了非飽和區(qū)的側向流。通過比對上述3種HSB模型由于對非飽和流處理方式不同而在模型的理論參數(shù)上產(chǎn)生的差異,分析非飽和流在邊坡內(nèi)部水文過程中的影響。再將3種模型與邊坡穩(wěn)定性模型耦合到一起,模擬淺層邊坡中的水文過程和邊坡穩(wěn)定性變化過程,進而分析非飽和側向滲流對邊坡穩(wěn)定計算結果的影響,從而為類似工程提供有價值的參考依據(jù)。

    第3期王銳,等:非飽和流對邊坡內(nèi)水文過程及穩(wěn)定性的影響分析河北工業(yè)科技第38卷

    1模型的比較與分析

    1.1變飽和邊坡水文模型的比較

    為將Boussinesq地下水方程用于描述邊坡內(nèi)部的水文過程,TROCH等[11]引入有效含水量S(x,t)來描述邊坡的幾何特征和內(nèi)部介質(zhì)的孔隙度,并假設垂直坡面方向上水壓滿足靜壓分布且垂直于x的水平方向無流動,得到了一維的水動力方程。

    1.2特征參數(shù)分析

    由式(4)可知,非飽和流的作用可分為2類:一是非飽和區(qū)垂向流,非飽和區(qū)中垂向流對邊坡儲水和排水都起著重要作用,在水位變化過程中非飽和區(qū)也會發(fā)生吸水或排水,從而影響坡內(nèi)部的水位變化和含水量的空間分布,這主要反映在參數(shù)有效孔隙度f中;二是非飽和區(qū)側向流,會促進邊坡內(nèi)的水分運移,這主要反映在參數(shù)滲透流量q中。通過比較式(4)可知:

    圖2所示為Kong-HSB模型和Troch-HSB模型中有效孔隙度的比值隨水位的變化。從圖中可以看出,對于砂土,在坡內(nèi)水位較低時,Kong-HSB模型和Troch-HSB模型的有效孔隙度幾乎相等,在坡內(nèi)較高水位才會出現(xiàn)明顯差異;而對于壤土,Kong-HSB模型和Troch-HSB模型有效孔隙度的差異隨著水位的抬高不斷拉大,即使當坡內(nèi)水位趨近于0時,Kong-HSB模型的有效孔隙度也僅相當于Troch-HSB模型的0.9倍。

    可見,在相同含水量變化時考慮非飽和流的影響下,水位的波動更大,同時在相同水位變化下,考慮非飽和流情況下的側向滲透流量更大。

    2非飽和區(qū)對水文變化的影響機制

    為揭示模型有效孔隙度參數(shù)f對邊坡內(nèi)部水位波動的影響機制,設定初始水位分別為0.4,0.8,1.0 m,含水量變化取值為0.05 m3,分別利用Kong-HSB模型和Troch-HSB模型模擬邊坡內(nèi)部含水量的空間分布。2種模型在各工況的水位變化量如表1所列,含水量的空間分布如圖3所示,虛線引向坐標軸表示的是邊坡內(nèi)部水位高度。從表1和圖3可以發(fā)現(xiàn),在保證初始水分分布和含水量變化量相同的情況下,Kong-HSB與Troch-HSB模型水位變化差異會受到土壤屬性和邊坡初始水位的影響,這與圖2中反映的有效孔隙度的相對關系一致。

    根據(jù)圖3與式(6)可知:在砂土中非飽和區(qū)持水性弱,非飽和區(qū)的含水量相對較小,Kong-HSB模型與Troch-HSB模型的有效孔隙度幾乎相等,所以2種模型計算得到的水位變化幾乎相同,只有當初始水位足夠高時,有效孔隙度差異所導致的水位波動差異才會開始顯現(xiàn);在壤土中,由于非飽和區(qū)厚度大且含水量也大,即使在初始水位很低的情況下2個模型的有效孔隙度也會有顯著差異,這種差異造成了如表1所示的2種模型水位波動的不同,且當初始水位升高后,2種模型的有效孔隙度的差異增大導致水位波動差異隨之增大。

    3非飽和流對邊坡穩(wěn)定性的影響

    3.1邊坡穩(wěn)定性分析模型

    相比于邊坡的長度和寬度,所討論的淺層邊坡厚度足夠小,可以采用極限平衡法結合半無限邊坡假定來分析邊坡穩(wěn)定性[16-17]。假定土條兩側側向力的影響可以忽略,即可以認為圖4中條塊的側面法向力E、剪切力X的增量ΔE,ΔX均為0。

    3.2非飽和流對邊坡穩(wěn)定系數(shù)的影響

    將3種HSB水文模型與邊坡穩(wěn)定性模型耦合,形成了Troch-HSB-SM模型、Hilberts-HSB-SM模型和Kong-HSB-SM模型。將非飽和流作用所造成邊坡內(nèi)部水位和含水量分布差異的相關參數(shù)引入到邊坡穩(wěn)定性計算,分析邊坡穩(wěn)定安全系數(shù)變化和失穩(wěn)時間差異。

    為分析河流筑壩、流域匯流情況下邊坡鄰近地域水位升降所引起的側向流動對邊坡穩(wěn)定性的影響,以壤土為例設計了2種坡長不同的工況:考慮坡長分別為10,50 m時,邊坡失穩(wěn)的擾動因素均為坡頂邊界水位抬升,詳細信息見表2,土壤參數(shù)設定如表3,其中的土壤參數(shù)α,β,n1,n2,Ks來自TROCH的研究[14]。

    模擬結果顯示,隨著時間的推移,水分滲入邊坡引起水位抬高,模型得到的邊坡穩(wěn)定系數(shù)逐漸減小。在初始階段3種模型穩(wěn)定系數(shù)變化非常接近,這是因為此時水位抬升所造成的水壓力占主導地位。之后對于工況1中較短的10 m壤土邊坡,Kong-HSB-SM模型預測的邊坡穩(wěn)定性的下降最快,也最早出現(xiàn)邊坡失穩(wěn)現(xiàn)象,其次是Hilberts-HSB-SM模型和Troch-HSB-SM模型。工況2中壤土邊坡坡長為50 m時,Kong-HSB-SM模型和Hilberts-HSB-SM模型預測的邊坡穩(wěn)定系數(shù)F幾乎同時以相似的速度下降并接近邊坡失穩(wěn),說明非飽和側向流對邊坡穩(wěn)定性的影響也會受到邊坡長度的影響,邊坡長度越短,非飽和流誘發(fā)的邊坡失穩(wěn)越明顯,失穩(wěn)時間也越提前。

    3.3非飽和流對邊坡穩(wěn)定影響的動力學機制

    根據(jù)圖7邊坡邊界流量-時間過程曲線來分析坡頂與坡腳的滲水過程,發(fā)現(xiàn)在入滲初期時,由于坡頂邊界處的水位不斷抬高產(chǎn)生的水力梯度使得此時的坡頂入滲量較大;而當坡頂邊界處水位抬高至恒定時,坡頂一側地下水位依舊升高,水力梯度減弱使得坡頂入滲量減少;當坡頂一側地下水位升至與其邊界水位近乎相同時,此階段水力梯度僅受坡度影響,此時的坡頂入滲量保持恒定。同樣在滲透初期,由于上邊界入滲水未到達坡腳,坡腳一側水位和滲出量保持恒定;隨著上游入滲水到達坡腳,坡腳水位抬高產(chǎn)生的水力梯度加大了坡腳處的滲出量;當坡腳滲出量增大到一定值時,會使得坡頂入滲量與坡腳滲出量達到平衡。

    根據(jù)圖7與上述分析,非飽和區(qū)側向流的入滲過程可分為3個階段,如圖8所示:第1階段為滲透初期,與坡腳非飽和區(qū)的滲出量相比,坡頂一側較大的水力梯度使更多的水分從坡頂非飽和區(qū)滲入邊坡,這個階段非飽和區(qū)側向流促進水分滲入邊坡,增加邊坡含水量,降低邊坡穩(wěn)定性;第2階段,升高的坡頂?shù)叵滤皇蛊渌μ荻葴p小,此時坡頂非飽和區(qū)入滲量減小,而坡腳非飽和區(qū)滲出量保持不變或略有增加,此階段非飽和區(qū)側向流作用微弱;第3階段,坡腳一側水力梯度增大使坡腳非飽和區(qū)滲出量增加,此時坡腳非飽和區(qū)滲出量大于坡頂非飽和區(qū)滲入量,非飽和區(qū)側向流協(xié)助邊坡的排水,降低邊坡內(nèi)的水位,緩解邊坡的失穩(wěn)性。

    而對滑坡失穩(wěn)問題而言,非飽和區(qū)側向流入滲過程一般處于第1或第2階段。根據(jù)上述理論可解釋圖6:對于10 m邊坡工況,在失穩(wěn)前非飽和區(qū)側向流長時間處于第1階段,非飽和側向流作用明顯,所以Kong-HSB-SM模型預測的邊坡穩(wěn)定性下降最快;而對于長度50 m邊坡工況,在失穩(wěn)前非飽和區(qū)側向流的作用長時間處于第2階段,此階段非飽和側向流作用微弱,所以Kong-HSB-SM模型和Hilberts-HSB-SM模型預測結果接近。

    4結語

    通過分析計算,比對了3種表述邊坡水文過程的HSB模型,分析非飽和滲流對邊坡內(nèi)部水文過程的影響和作用機制,并耦合HSB模型與邊坡穩(wěn)定極限平衡模型,模擬了淺層邊坡在不同水動力條件下的安全系數(shù)變化情況,探討非飽和側向滲流對淺層邊坡穩(wěn)定性的影響。研究主要結論如下。

    1)對于非飽和流顯著的壤土,在相同含水量變化量的條件下,考慮非飽和流模型的水位波動更大;在相同水位的條件下,考慮非飽和流的模型水分運移更多。

    2)在以側向流為主要滑坡誘發(fā)因素的淺層滑坡中,壤土邊坡中非飽和流作用可以促進水分在邊坡內(nèi)部的運移,導致邊坡穩(wěn)定性的降低,忽視非飽和流的邊坡失穩(wěn)時間預測會滯后。

    3)非飽和區(qū)側向流的作用可劃分成3個階段,第1階段非飽和區(qū)側向流促進水分滲入邊坡,降低邊坡穩(wěn)定性;第2階段非飽和區(qū)側向流作用減弱;第3階段非飽和區(qū)側向流能協(xié)助邊坡的排水,降低邊坡內(nèi)的水位,緩解邊坡的失穩(wěn)性。

    本文假定研究對象為均質(zhì)各向同性邊坡,但在自然環(huán)境中非均質(zhì)性邊坡更為常見,且本文未能考慮更復雜的邊坡幾何形態(tài)(如起伏、發(fā)散、收斂等特征),這些可能會改變邊坡內(nèi)部側向流特征,對研究結果有一定的影響;今后可以通過引入更多的參數(shù)來描述邊坡地質(zhì)和形態(tài),探究其非飽和流的作用。另外,本文主要分析了不受降雨入滲影響的淺層滑坡,在后續(xù)研究中可考慮降雨和側向流的共同作用促成的淺層滑坡。

    參考文獻/References:

    [1]宋玉延. 滑坡治理雙排抗滑樁受力特征及設計計算方法研究[D]. 武漢: 武漢理工大學, 2013.

    SONG Yuyan. Study on Force Characteristics and Design and Calculation[D]. Wuhan: Wuhan University of Technology, 2013.

    [2]SINGH T N, VERMA A K, SARKAR K. Static and dynamic analysis of a landslide[J]. Geomatics Natural Hazards and Risk, 2010, 1(4): 323-338.

    [3]PETLEY D. Global patterns of loss of life fromlandslides[J]. Geology, 2012, 40(10): 927-930.

    [4]CHEN X Z, CUI Y F. The formation ofthe Wulipo landslide and the resulting debris flow in Dujiangyan City, China[J]. Journal of Mountain Science, 2017, 14(6): 1100-1112.

    [5]CONTE E, DONATO A, TRONCONE A. A simplified method for predicting rainfall-induced mobility of activelandslides[J]. Landslides, 2017, 14(1): 35-45.

    [6]BRESCIANI E, DAVY P, DREUZY J R. Is the Dupuit assumption suitable for predicting the groundwater seepage area in hillslopes?[J]. Water Resources Research, 2014, 50(3): 2394-2406.

    [7]陳靈芝, 古振榮, 李揚. 降雨入滲條件下非飽和土高邊坡穩(wěn)定性有限元分析[J]. 齊齊哈爾大學學報(自然科學版), 2020, 36(6): 77-82.

    CHEN Lingzhi, GU Zhenrong, LI Yang. Finite element analysis of high slope stability of unsaturated soil under rainfall infiltration[J]. Journal of Qiqihar University(Natural Science Edition), 2020, 36(6): 77-82.

    [8]劉鋒, 胡江洋, 劉軍勇, 等. 連續(xù)降雨入滲作用下非飽和土路基穩(wěn)定性研究[J]. 路基工程, 2020(1): 55-60.

    LIU Feng, HU Jiangyang, LIU Junyong, et al. Study on the stability of unsaturated soil subgrade under continuous rainfall infiltration[J]. Subgrade Engineering, 2020(1): 55-60.

    [9]毋遠召, 馬文禮, 魏占璽, 等. 降雨及滲流條件下層狀砂質(zhì)板巖邊坡變形破壞模式研究[J]. 河北工業(yè)科技, 2020, 37(4): 230-236.

    WU Yuanzhao, MA Wenli, WEI Zhanxi, et al. Research on deformation and failure modes of layered sandy slate slopes under rainfall and seepage conditions[J]. Hebei Journal of Industrial Science and Technology, 2020, 37(4): 230-236.

    [10]周娟, 朱峰. 中西部地區(qū)黏性土質(zhì)邊坡降雨動態(tài)滲流的非飽和特性分析[J]. 河北工業(yè)科技, 2020, 37(4): 273-279.

    ZHOU Juan, ZHU Feng. Analysis on unsaturated characteristics of dynamic rainfall seepage of clay soil slope in Central and Western China[J]. Hebei Journal of Industrial Science and Technology, 2020, 37(4): 273-279.

    [11]TROCH P A, PANICONI C, van LOON E E. Hillslope-storage Boussinesq model for subsurface flow and variable source areas along complex hillslopes: 1. Formulation and characteristic response[J]. Water Resources Research, 2003, 39(11): 1316-1328.

    [12]HILBERTS A G J, TROCH P A, PANICONI C, et al. Low-dimensional modeling ofhillslope subsurface flow: Relationship between rainfall, recharge, and unsaturated storage dynamics[J]. Water Resources Research, 2007, 43(3): 3445-3459.

    [13]KONG J, SHEN C J, LUO Z Y, et al. Improvement of thehillslope-storage Boussinesq model by considering lateral flow in the unsaturated zone[J]. Water Resources Research, 2016, 52(4): 2965-2984.

    [14]TROCH P A. Conceptual Basin-Scale Runoff Process Models for Humid Catchments: Analysis, Synthesis and Applications[D]. Ghent: Ghent University, 1993.

    [15]van GENUCHTEN M T. A closed-form equation for predicting the hydraulic conductivity of unsaturatedsoils[J]. Soil Science Society of America Journal, 1980, 44(5): 892-898.

    [16]李寧, 許建聰, 欽亞洲.降雨誘發(fā)淺層滑坡穩(wěn)定性的計算模型研究[J]. 巖土力學, 2012, 33(5): 209-214.

    LI Ning, XU Jiancong, QIN Yazhou,et al. Research on calculation model for stability evaluation of rainfall-induced shallow landslides[J]. Rock and Soil Mechanics, 2012,33(5): 209-214.

    [17]彭子茂, 黃震. 模擬非飽和流的新型相對滲透系數(shù)模型[J]. 長江科學院院報, 2020, 37(7): 115-119.

    PENG Zimao, HUANG Zhen. A relative permeability model for unsaturated flow in consideration of hysteresis effect[J]. Journal of Yangtze River Scientific Research Institute, 2020, 37(7): 115-119.

    猜你喜歡
    邊坡穩(wěn)定性滑坡
    水庫滑坡地下水動態(tài)響應規(guī)律及浸潤線計算模型——以石榴樹包滑坡為例
    滑坡推力隱式解與顯式解對比分析——以河北某膨脹土滑坡為例
    淺談公路滑坡治理
    北方交通(2016年12期)2017-01-15 13:52:59
    某水電站混凝土生產(chǎn)系統(tǒng)邊坡穩(wěn)定性復核分析
    某邊坡地質(zhì)災害隱患點治理工程勘查
    綠色科技(2016年20期)2016-12-27 16:23:03
    基于GCSF算法的邊坡穩(wěn)定性分析軟件設計
    公路邊坡穩(wěn)定性及防護加固技術
    組合式邊坡支護技術在深大基坑工程中的實踐分析
    基于Fluent的滑坡入水過程數(shù)值模擬
    淺談邊坡穩(wěn)定性地質(zhì)問題的解決措施
    欧美激情高清一区二区三区 | www日本在线高清视频| 韩国高清视频一区二区三区| 国产伦理片在线播放av一区| 亚洲av日韩在线播放| 久久精品久久久久久噜噜老黄| 中文字幕最新亚洲高清| 精品久久久精品久久久| 97在线人人人人妻| 免费观看人在逋| 久久久精品国产亚洲av高清涩受| 日本av手机在线免费观看| 国产黄色视频一区二区在线观看| 成年动漫av网址| 精品免费久久久久久久清纯 | 咕卡用的链子| 九九爱精品视频在线观看| 精品人妻熟女毛片av久久网站| 汤姆久久久久久久影院中文字幕| 黑人欧美特级aaaaaa片| 999精品在线视频| 久久亚洲国产成人精品v| 久久天躁狠狠躁夜夜2o2o | 黄色视频在线播放观看不卡| 国产精品国产av在线观看| 亚洲精品aⅴ在线观看| 日韩中文字幕视频在线看片| 美女国产高潮福利片在线看| 日韩精品免费视频一区二区三区| 天天操日日干夜夜撸| 亚洲成国产人片在线观看| 亚洲国产最新在线播放| 欧美精品高潮呻吟av久久| 色网站视频免费| 亚洲av成人精品一二三区| 免费观看性生交大片5| 婷婷色av中文字幕| 九草在线视频观看| 久久久久精品久久久久真实原创| 欧美日韩国产mv在线观看视频| 久久人人97超碰香蕉20202| 看免费成人av毛片| 亚洲国产av新网站| 综合色丁香网| www.自偷自拍.com| 香蕉丝袜av| 99热网站在线观看| 日本欧美国产在线视频| 午夜福利一区二区在线看| 国产精品 国内视频| 亚洲成国产人片在线观看| 中文字幕精品免费在线观看视频| 久久婷婷青草| 国产女主播在线喷水免费视频网站| 亚洲国产欧美在线一区| 一二三四在线观看免费中文在| 日韩精品免费视频一区二区三区| 丝袜人妻中文字幕| 免费观看性生交大片5| 久久久久久久大尺度免费视频| 精品少妇黑人巨大在线播放| 亚洲av成人精品一二三区| 男男h啪啪无遮挡| 天天躁夜夜躁狠狠躁躁| 一级毛片黄色毛片免费观看视频| av网站在线播放免费| 在线 av 中文字幕| 1024香蕉在线观看| 只有这里有精品99| av在线播放精品| 天天躁夜夜躁狠狠躁躁| 国产有黄有色有爽视频| 十分钟在线观看高清视频www| 亚洲欧美一区二区三区久久| 国产亚洲精品第一综合不卡| 国产成人精品无人区| 中文精品一卡2卡3卡4更新| 久久精品亚洲熟妇少妇任你| 午夜免费男女啪啪视频观看| 日本色播在线视频| 欧美乱码精品一区二区三区| avwww免费| 母亲3免费完整高清在线观看| 日本猛色少妇xxxxx猛交久久| 一级a爱视频在线免费观看| 久热这里只有精品99| 日韩一本色道免费dvd| 夫妻午夜视频| 超碰97精品在线观看| 久久鲁丝午夜福利片| 日本一区二区免费在线视频| 超碰成人久久| 欧美97在线视频| 伊人亚洲综合成人网| 久久久久精品国产欧美久久久 | 黄片播放在线免费| 三上悠亚av全集在线观看| 看十八女毛片水多多多| 啦啦啦在线观看免费高清www| 五月天丁香电影| 亚洲精品久久成人aⅴ小说| 天天躁夜夜躁狠狠久久av| 激情五月婷婷亚洲| 国产欧美日韩一区二区三区在线| 99精国产麻豆久久婷婷| 欧美另类一区| 女的被弄到高潮叫床怎么办| 精品国产乱码久久久久久男人| 久久久久视频综合| 国产片特级美女逼逼视频| 国产日韩欧美在线精品| av网站免费在线观看视频| 国产亚洲最大av| 久久精品久久久久久久性| 老汉色av国产亚洲站长工具| 国产一区二区在线观看av| 又黄又粗又硬又大视频| 国产精品久久久人人做人人爽| 日韩大码丰满熟妇| 美女福利国产在线| 精品一品国产午夜福利视频| 成人手机av| 久久97久久精品| 亚洲少妇的诱惑av| 久久久亚洲精品成人影院| 欧美少妇被猛烈插入视频| 九色亚洲精品在线播放| 五月开心婷婷网| 99热网站在线观看| 亚洲熟女毛片儿| 国产在线视频一区二区| 成人三级做爰电影| 久久久久久人人人人人| 大话2 男鬼变身卡| 青春草亚洲视频在线观看| 黄色 视频免费看| 国产成人精品在线电影| 激情五月婷婷亚洲| 日日啪夜夜爽| 国产精品国产三级专区第一集| 啦啦啦中文免费视频观看日本| 女人高潮潮喷娇喘18禁视频| 亚洲欧美精品自产自拍| 人人妻人人添人人爽欧美一区卜| 免费黄网站久久成人精品| 黄频高清免费视频| 久久综合国产亚洲精品| 午夜福利视频在线观看免费| 国产黄色免费在线视频| 精品久久久精品久久久| 高清视频免费观看一区二区| 日韩精品有码人妻一区| 少妇 在线观看| 亚洲成人国产一区在线观看 | 91精品国产国语对白视频| 国产爽快片一区二区三区| 你懂的网址亚洲精品在线观看| 国产激情久久老熟女| 国产亚洲一区二区精品| 久久人妻熟女aⅴ| 男女高潮啪啪啪动态图| 极品少妇高潮喷水抽搐| 十八禁人妻一区二区| 一二三四中文在线观看免费高清| 久久婷婷青草| 日韩一卡2卡3卡4卡2021年| 女性被躁到高潮视频| a级片在线免费高清观看视频| av福利片在线| 国产黄色视频一区二区在线观看| 国产av精品麻豆| 男女边吃奶边做爰视频| 黑人巨大精品欧美一区二区蜜桃| 婷婷色麻豆天堂久久| 国产免费视频播放在线视频| 亚洲国产精品国产精品| 久久久久久久久久久久大奶| 国产免费一区二区三区四区乱码| 另类亚洲欧美激情| 七月丁香在线播放| 国产精品久久久久久精品电影小说| 免费在线观看完整版高清| 日韩一卡2卡3卡4卡2021年| 午夜激情av网站| 尾随美女入室| 少妇人妻久久综合中文| 丝袜脚勾引网站| 亚洲精品美女久久av网站| 最近最新中文字幕免费大全7| 97精品久久久久久久久久精品| 亚洲国产日韩一区二区| 国产97色在线日韩免费| 欧美亚洲日本最大视频资源| 欧美精品人与动牲交sv欧美| 美国免费a级毛片| 亚洲视频免费观看视频| 狂野欧美激情性xxxx| 热re99久久精品国产66热6| 大陆偷拍与自拍| 国产精品亚洲av一区麻豆 | 性少妇av在线| 免费看av在线观看网站| a级毛片在线看网站| 卡戴珊不雅视频在线播放| 国产免费又黄又爽又色| 日韩不卡一区二区三区视频在线| 国产精品一二三区在线看| 在线天堂最新版资源| 国产一区二区 视频在线| 国产日韩欧美亚洲二区| 在线观看免费日韩欧美大片| 咕卡用的链子| 国产成人精品福利久久| 久久影院123| 日日啪夜夜爽| 中文字幕人妻熟女乱码| 国产乱来视频区| 老司机影院毛片| 久久久久精品国产欧美久久久 | 婷婷色综合www| 久久亚洲国产成人精品v| 桃花免费在线播放| 男人舔女人的私密视频| 一区二区三区精品91| 天天躁日日躁夜夜躁夜夜| 欧美黄色片欧美黄色片| 国产激情久久老熟女| 日韩制服丝袜自拍偷拍| 色精品久久人妻99蜜桃| av不卡在线播放| 99re6热这里在线精品视频| 成年女人毛片免费观看观看9 | 欧美在线黄色| 亚洲国产av新网站| 久久久久久久大尺度免费视频| 伊人亚洲综合成人网| 丝袜美足系列| 九草在线视频观看| 国产精品 欧美亚洲| 久久人人97超碰香蕉20202| 视频在线观看一区二区三区| 日本欧美国产在线视频| 99国产综合亚洲精品| 精品第一国产精品| 午夜影院在线不卡| 亚洲av男天堂| 丝瓜视频免费看黄片| 满18在线观看网站| 最近手机中文字幕大全| 国产精品.久久久| 色综合欧美亚洲国产小说| 999久久久国产精品视频| 午夜日本视频在线| 老司机影院成人| 性高湖久久久久久久久免费观看| 天美传媒精品一区二区| 久久久欧美国产精品| 国产97色在线日韩免费| 亚洲av在线观看美女高潮| 久久精品亚洲av国产电影网| 侵犯人妻中文字幕一二三四区| 国产av精品麻豆| 日本爱情动作片www.在线观看| 国产淫语在线视频| 欧美日韩视频高清一区二区三区二| 91成人精品电影| 两性夫妻黄色片| 丝袜脚勾引网站| 99国产综合亚洲精品| 最近中文字幕高清免费大全6| 欧美97在线视频| 国产一区二区在线观看av| 超碰97精品在线观看| 男的添女的下面高潮视频| 王馨瑶露胸无遮挡在线观看| 精品酒店卫生间| 成人18禁高潮啪啪吃奶动态图| 国产色婷婷99| 18禁裸乳无遮挡动漫免费视频| 999久久久国产精品视频| av在线观看视频网站免费| 一区二区三区乱码不卡18| 十八禁网站网址无遮挡| 精品一区二区三区av网在线观看 | 深夜精品福利| 丝袜人妻中文字幕| a 毛片基地| 香蕉丝袜av| 美女福利国产在线| 女性生殖器流出的白浆| 一区二区三区四区激情视频| 国产探花极品一区二区| 国产日韩一区二区三区精品不卡| 蜜桃在线观看..| 午夜激情久久久久久久| 国产野战对白在线观看| av片东京热男人的天堂| 老熟女久久久| 国产一区二区在线观看av| 欧美 日韩 精品 国产| 免费高清在线观看视频在线观看| 少妇人妻精品综合一区二区| 亚洲国产欧美一区二区综合| 欧美日韩精品网址| 男的添女的下面高潮视频| 国产精品久久久久成人av| 女人久久www免费人成看片| 9色porny在线观看| 宅男免费午夜| 国产高清国产精品国产三级| 久久天堂一区二区三区四区| av天堂久久9| 日韩欧美一区视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | e午夜精品久久久久久久| 如何舔出高潮| 99精品久久久久人妻精品| 热re99久久精品国产66热6| 午夜日本视频在线| 麻豆乱淫一区二区| 又大又爽又粗| 日韩中文字幕欧美一区二区 | 久久久久人妻精品一区果冻| 1024香蕉在线观看| 亚洲国产精品一区三区| 一本一本久久a久久精品综合妖精| 国产精品欧美亚洲77777| 日韩 欧美 亚洲 中文字幕| 国产成人欧美在线观看 | 黄频高清免费视频| 在现免费观看毛片| 亚洲国产精品一区二区三区在线| 色婷婷久久久亚洲欧美| 大话2 男鬼变身卡| 曰老女人黄片| 国产女主播在线喷水免费视频网站| 欧美人与性动交α欧美软件| 国产精品一区二区在线观看99| 看免费av毛片| 欧美激情极品国产一区二区三区| 91老司机精品| 成年人免费黄色播放视频| 黄色怎么调成土黄色| 欧美黑人欧美精品刺激| 观看美女的网站| 亚洲,欧美精品.| av天堂久久9| 亚洲视频免费观看视频| 精品第一国产精品| 91精品三级在线观看| 91老司机精品| 啦啦啦 在线观看视频| 黄片播放在线免费| 人人妻人人爽人人添夜夜欢视频| 天美传媒精品一区二区| 大香蕉久久网| 嫩草影院入口| 中文字幕另类日韩欧美亚洲嫩草| 老司机亚洲免费影院| 一级毛片黄色毛片免费观看视频| 国产精品熟女久久久久浪| av.在线天堂| 欧美人与性动交α欧美软件| 高清黄色对白视频在线免费看| 又黄又粗又硬又大视频| 亚洲一码二码三码区别大吗| 亚洲国产精品成人久久小说| 人成视频在线观看免费观看| 国产精品99久久99久久久不卡 | 在线亚洲精品国产二区图片欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 在线亚洲精品国产二区图片欧美| 日韩制服骚丝袜av| 午夜福利视频在线观看免费| 午夜激情久久久久久久| 午夜福利一区二区在线看| 一级毛片电影观看| 亚洲男人天堂网一区| 青青草视频在线视频观看| 中文字幕人妻丝袜一区二区 | 国产伦人伦偷精品视频| 国产精品久久久久久精品电影小说| 人成视频在线观看免费观看| 一级黄片播放器| 欧美日韩成人在线一区二区| 性色av一级| 如日韩欧美国产精品一区二区三区| 丰满饥渴人妻一区二区三| 国产黄频视频在线观看| 晚上一个人看的免费电影| 欧美日韩福利视频一区二区| 亚洲少妇的诱惑av| 自拍欧美九色日韩亚洲蝌蚪91| 日韩 欧美 亚洲 中文字幕| 国产男女内射视频| 黄片小视频在线播放| 久久久久久久久免费视频了| av有码第一页| 亚洲欧美清纯卡通| 天天操日日干夜夜撸| 最近2019中文字幕mv第一页| 欧美日韩成人在线一区二区| 中文字幕亚洲精品专区| 日日撸夜夜添| 男人添女人高潮全过程视频| 国产亚洲最大av| 亚洲男人天堂网一区| 黄色怎么调成土黄色| 九草在线视频观看| 一级毛片我不卡| 亚洲欧美精品自产自拍| 亚洲五月色婷婷综合| 91国产中文字幕| 如何舔出高潮| 亚洲av男天堂| 在线观看免费视频网站a站| 久久久久精品国产欧美久久久 | 久久久久久人人人人人| 日日撸夜夜添| 成年动漫av网址| 欧美精品av麻豆av| 精品少妇黑人巨大在线播放| 王馨瑶露胸无遮挡在线观看| 在线看a的网站| 亚洲欧美一区二区三区国产| 欧美日韩视频精品一区| 十八禁高潮呻吟视频| 可以免费在线观看a视频的电影网站 | 日本vs欧美在线观看视频| 欧美av亚洲av综合av国产av | 成年人免费黄色播放视频| 99久久综合免费| 亚洲第一区二区三区不卡| av女优亚洲男人天堂| 丝袜人妻中文字幕| 亚洲专区中文字幕在线 | 成人午夜精彩视频在线观看| 欧美精品av麻豆av| 精品酒店卫生间| 亚洲综合精品二区| 久久天躁狠狠躁夜夜2o2o | 国产一区二区三区综合在线观看| 亚洲国产最新在线播放| 国产男人的电影天堂91| 多毛熟女@视频| 人人妻人人爽人人添夜夜欢视频| 精品少妇内射三级| 男女下面插进去视频免费观看| 久久久久国产精品人妻一区二区| 别揉我奶头~嗯~啊~动态视频 | 久久久久久人人人人人| 中国三级夫妇交换| 少妇被粗大的猛进出69影院| 日本av手机在线免费观看| 波野结衣二区三区在线| 国产精品久久久久久精品电影小说| 汤姆久久久久久久影院中文字幕| 久久久久久人妻| 狂野欧美激情性xxxx| 另类亚洲欧美激情| 久久国产亚洲av麻豆专区| 欧美日韩综合久久久久久| 在现免费观看毛片| 青草久久国产| 国产日韩欧美亚洲二区| 男女国产视频网站| 国产一区二区激情短视频 | 大片免费播放器 马上看| 免费观看av网站的网址| 在线免费观看不下载黄p国产| 两性夫妻黄色片| 亚洲av在线观看美女高潮| 国产亚洲欧美精品永久| 91精品三级在线观看| av国产久精品久网站免费入址| 亚洲精品aⅴ在线观看| 妹子高潮喷水视频| 女人精品久久久久毛片| 午夜免费男女啪啪视频观看| 久久韩国三级中文字幕| 欧美精品一区二区大全| 欧美成人精品欧美一级黄| 一级毛片 在线播放| 在线观看免费视频网站a站| 观看美女的网站| 精品国产一区二区三区久久久樱花| a 毛片基地| av免费观看日本| 丝袜喷水一区| 国产片内射在线| 女性生殖器流出的白浆| 下体分泌物呈黄色| 精品国产乱码久久久久久小说| 少妇精品久久久久久久| 两个人看的免费小视频| 免费高清在线观看日韩| 国产黄色视频一区二区在线观看| 久久精品久久久久久噜噜老黄| 日本欧美视频一区| 黄色怎么调成土黄色| 精品国产乱码久久久久久男人| 精品少妇一区二区三区视频日本电影 | 久久99热这里只频精品6学生| 精品国产露脸久久av麻豆| 亚洲精品第二区| 少妇被粗大猛烈的视频| 男女无遮挡免费网站观看| 观看美女的网站| 肉色欧美久久久久久久蜜桃| 黄网站色视频无遮挡免费观看| 极品人妻少妇av视频| 人妻 亚洲 视频| 亚洲色图综合在线观看| 国产女主播在线喷水免费视频网站| 婷婷成人精品国产| 亚洲国产欧美在线一区| 久久精品国产亚洲av涩爱| 欧美 亚洲 国产 日韩一| 汤姆久久久久久久影院中文字幕| 日日啪夜夜爽| 中国国产av一级| 大陆偷拍与自拍| 国产亚洲av高清不卡| 日韩制服骚丝袜av| 国产不卡av网站在线观看| 精品少妇内射三级| 亚洲精品自拍成人| 嫩草影视91久久| 欧美少妇被猛烈插入视频| videosex国产| 97精品久久久久久久久久精品| 亚洲欧美成人综合另类久久久| 一区福利在线观看| 中文字幕高清在线视频| 亚洲国产欧美一区二区综合| 自线自在国产av| √禁漫天堂资源中文www| 性高湖久久久久久久久免费观看| 国产精品一区二区精品视频观看| 日本wwww免费看| 成人午夜精彩视频在线观看| 亚洲国产欧美网| 18在线观看网站| videos熟女内射| 久久人人爽av亚洲精品天堂| 精品免费久久久久久久清纯 | 女人久久www免费人成看片| 水蜜桃什么品种好| 一区二区av电影网| 国产精品久久久人人做人人爽| 好男人视频免费观看在线| 这个男人来自地球电影免费观看 | 激情五月婷婷亚洲| 国产av精品麻豆| 久久精品久久久久久久性| 狂野欧美激情性xxxx| 不卡视频在线观看欧美| 美女国产高潮福利片在线看| bbb黄色大片| 久久毛片免费看一区二区三区| 久久精品国产综合久久久| 大香蕉久久成人网| 欧美精品高潮呻吟av久久| 久久精品久久精品一区二区三区| 国产成人精品久久二区二区91 | 两性夫妻黄色片| 日韩视频在线欧美| 美女福利国产在线| 免费高清在线观看视频在线观看| 亚洲国产毛片av蜜桃av| xxxhd国产人妻xxx| 男的添女的下面高潮视频| 国产深夜福利视频在线观看| 青春草亚洲视频在线观看| 精品免费久久久久久久清纯 | 色网站视频免费| 久久久精品区二区三区| 中文精品一卡2卡3卡4更新| 亚洲中文av在线| 国产亚洲av高清不卡| 久久久久视频综合| 国产成人精品福利久久| 在线免费观看不下载黄p国产| 成人国产av品久久久| 色综合欧美亚洲国产小说| 亚洲中文av在线| 欧美成人午夜精品| 国产成人欧美在线观看 | 精品一区二区免费观看| 一本一本久久a久久精品综合妖精| 色婷婷久久久亚洲欧美| 久久久精品国产亚洲av高清涩受| 国产伦理片在线播放av一区| 国产老妇伦熟女老妇高清| 精品一区二区三区四区五区乱码 | 中文字幕另类日韩欧美亚洲嫩草| 欧美 日韩 精品 国产| www.熟女人妻精品国产| 99国产综合亚洲精品| 如何舔出高潮| 亚洲欧洲精品一区二区精品久久久 | 亚洲免费av在线视频| 男女床上黄色一级片免费看| 国产老妇伦熟女老妇高清| 日韩制服骚丝袜av| 99re6热这里在线精品视频| 国产精品99久久99久久久不卡 | 蜜桃在线观看..| 热99国产精品久久久久久7| 视频区图区小说| 熟女av电影| 下体分泌物呈黄色| 亚洲免费av在线视频| 黄色一级大片看看| 日韩一区二区三区影片|