• <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ì)問題的解決措施
    国产探花极品一区二区| 日日摸夜夜添夜夜添av毛片| 中文字幕亚洲精品专区| 中文字幕制服av| 亚洲欧美精品综合久久99| 一级av片app| 老师上课跳d突然被开到最大视频| 午夜精品一区二区三区免费看| 看黄色毛片网站| 国产精品乱码一区二三区的特点| 小说图片视频综合网站| 人妻夜夜爽99麻豆av| 神马国产精品三级电影在线观看| 亚洲成人中文字幕在线播放| 国产一级毛片七仙女欲春2| 亚洲精品自拍成人| 激情 狠狠 欧美| 国产真实乱freesex| 99热这里只有是精品50| 一个人看的www免费观看视频| 亚洲国产精品sss在线观看| 亚洲最大成人av| 亚洲欧美中文字幕日韩二区| 午夜精品一区二区三区免费看| 国产精品久久久久久久电影| 国产精品久久久久久精品电影| 99国产精品一区二区蜜桃av| 蜜桃亚洲精品一区二区三区| 黄色一级大片看看| 日韩制服骚丝袜av| 男女视频在线观看网站免费| 免费黄网站久久成人精品| av在线亚洲专区| 亚洲精品日韩在线中文字幕| 一级二级三级毛片免费看| 51国产日韩欧美| 美女被艹到高潮喷水动态| 永久网站在线| 麻豆国产97在线/欧美| 2021少妇久久久久久久久久久| 亚洲精品乱码久久久v下载方式| 内射极品少妇av片p| 免费看a级黄色片| 有码 亚洲区| 男人和女人高潮做爰伦理| 亚洲天堂国产精品一区在线| 亚洲av二区三区四区| 少妇丰满av| 深爱激情五月婷婷| 18禁动态无遮挡网站| 人妻制服诱惑在线中文字幕| 日韩中字成人| 日日干狠狠操夜夜爽| 成人二区视频| 亚洲无线观看免费| 国产 一区精品| 18禁在线播放成人免费| 又爽又黄无遮挡网站| 久久久久免费精品人妻一区二区| 婷婷色综合大香蕉| 丝袜美腿在线中文| 99热这里只有精品一区| 人体艺术视频欧美日本| 国产成人91sexporn| 国产精品久久视频播放| 国产又黄又爽又无遮挡在线| 精品不卡国产一区二区三区| 又粗又硬又长又爽又黄的视频| 亚洲四区av| 观看免费一级毛片| 日韩精品青青久久久久久| 91久久精品国产一区二区三区| 国产美女午夜福利| 欧美日韩国产亚洲二区| 国内揄拍国产精品人妻在线| 亚洲欧美清纯卡通| 免费一级毛片在线播放高清视频| 精品国内亚洲2022精品成人| 日本黄色视频三级网站网址| 美女黄网站色视频| 最后的刺客免费高清国语| 午夜久久久久精精品| 人体艺术视频欧美日本| 丰满乱子伦码专区| 99热6这里只有精品| 日韩在线高清观看一区二区三区| 蜜臀久久99精品久久宅男| 国产国拍精品亚洲av在线观看| 免费观看的影片在线观看| 99热这里只有是精品50| 边亲边吃奶的免费视频| 18禁在线播放成人免费| 国产精品久久电影中文字幕| 国产精品1区2区在线观看.| 九九久久精品国产亚洲av麻豆| 日日摸夜夜添夜夜爱| 蜜桃久久精品国产亚洲av| 成人亚洲欧美一区二区av| av免费在线看不卡| 久久精品夜色国产| 久久久久国产网址| 日韩欧美三级三区| 国产片特级美女逼逼视频| 一级毛片aaaaaa免费看小| 在线观看av片永久免费下载| 午夜亚洲福利在线播放| 日韩欧美国产在线观看| 亚洲精品aⅴ在线观看| 狂野欧美激情性xxxx在线观看| 国产视频内射| 一个人看视频在线观看www免费| 国产一区二区三区av在线| 一级二级三级毛片免费看| 建设人人有责人人尽责人人享有的 | 乱码一卡2卡4卡精品| 网址你懂的国产日韩在线| 美女国产视频在线观看| 国产片特级美女逼逼视频| 国产午夜精品一二区理论片| 亚洲精品成人久久久久久| 亚洲精品影视一区二区三区av| 国产亚洲av片在线观看秒播厂 | 久久久亚洲精品成人影院| 97在线视频观看| 建设人人有责人人尽责人人享有的 | 神马国产精品三级电影在线观看| 在线观看美女被高潮喷水网站| 成人鲁丝片一二三区免费| 人妻制服诱惑在线中文字幕| 黑人高潮一二区| 国产人妻一区二区三区在| 一级毛片aaaaaa免费看小| 久久99精品国语久久久| 亚洲国产精品成人综合色| 级片在线观看| 久久欧美精品欧美久久欧美| 亚洲五月天丁香| 国产又黄又爽又无遮挡在线| 亚洲精品乱码久久久久久按摩| 偷拍熟女少妇极品色| 18禁动态无遮挡网站| 插逼视频在线观看| 嫩草影院精品99| 久久草成人影院| 国产精品一区二区在线观看99 | 久久综合国产亚洲精品| 日本三级黄在线观看| 午夜福利在线观看免费完整高清在| 国产一区二区在线av高清观看| 精品无人区乱码1区二区| 国产免费男女视频| 黑人高潮一二区| 免费在线观看成人毛片| 亚洲国产精品合色在线| 男女边吃奶边做爰视频| 国产av一区在线观看免费| 大话2 男鬼变身卡| 最近视频中文字幕2019在线8| 国产精品蜜桃在线观看| 在线a可以看的网站| 天堂√8在线中文| 青春草国产在线视频| 免费看光身美女| 亚洲国产精品专区欧美| 非洲黑人性xxxx精品又粗又长| 不卡视频在线观看欧美| 欧美精品国产亚洲| 一级av片app| 三级男女做爰猛烈吃奶摸视频| 一级毛片久久久久久久久女| 亚洲精品456在线播放app| 亚洲性久久影院| 男的添女的下面高潮视频| 欧美激情久久久久久爽电影| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久久久久免费av| 在线a可以看的网站| 国产 一区精品| 亚州av有码| 久久99热这里只频精品6学生 | 国产在线一区二区三区精 | 草草在线视频免费看| 亚洲aⅴ乱码一区二区在线播放| 中文字幕精品亚洲无线码一区| 天堂网av新在线| 人体艺术视频欧美日本| 久久99热这里只有精品18| 免费观看性生交大片5| 日韩一区二区三区影片| 麻豆一二三区av精品| 久久久亚洲精品成人影院| 国产av不卡久久| 免费搜索国产男女视频| 最近的中文字幕免费完整| 亚洲在线观看片| 欧美激情久久久久久爽电影| 国产乱人视频| 国产成人免费观看mmmm| 久久久久久久久久久丰满| 夫妻性生交免费视频一级片| 亚洲不卡免费看| 成人欧美大片| 亚洲欧美日韩东京热| av国产免费在线观看| 嫩草影院精品99| 亚洲欧美一区二区三区国产| 精品久久久噜噜| 欧美精品一区二区大全| 最近最新中文字幕免费大全7| 国产真实伦视频高清在线观看| 一个人免费在线观看电影| 精品无人区乱码1区二区| 欧美97在线视频| 天堂中文最新版在线下载 | 激情 狠狠 欧美| 韩国av在线不卡| 午夜免费男女啪啪视频观看| 亚洲精品亚洲一区二区| 久久精品综合一区二区三区| 国产人妻一区二区三区在| 免费一级毛片在线播放高清视频| 22中文网久久字幕| 麻豆国产97在线/欧美| 亚洲人与动物交配视频| 午夜免费男女啪啪视频观看| av在线观看视频网站免费| 九九热线精品视视频播放| 在线免费观看不下载黄p国产| 亚洲国产日韩欧美精品在线观看| 久久99蜜桃精品久久| 禁无遮挡网站| 白带黄色成豆腐渣| 深爱激情五月婷婷| 日韩一区二区视频免费看| 日本午夜av视频| 1024手机看黄色片| 美女脱内裤让男人舔精品视频| 久久综合国产亚洲精品| 麻豆一二三区av精品| 99热网站在线观看| 国产精品蜜桃在线观看| 久久热精品热| 亚洲精品影视一区二区三区av| 亚洲欧美成人精品一区二区| 亚洲内射少妇av| 亚洲av电影不卡..在线观看| 亚洲精品,欧美精品| 一本一本综合久久| 色综合色国产| 亚洲av一区综合| av在线老鸭窝| 草草在线视频免费看| 午夜激情福利司机影院| 中文天堂在线官网| 午夜福利视频1000在线观看| 国内精品一区二区在线观看| 国产精品女同一区二区软件| 联通29元200g的流量卡| 成人国产麻豆网| a级毛色黄片| 午夜福利在线观看免费完整高清在| 91在线精品国自产拍蜜月| av专区在线播放| 亚洲av男天堂| 国产单亲对白刺激| 欧美又色又爽又黄视频| 国产真实乱freesex| 国产老妇女一区| 欧美性猛交╳xxx乱大交人| 日韩视频在线欧美| 日韩三级伦理在线观看| 直男gayav资源| 国产 一区精品| 国产麻豆成人av免费视频| 人体艺术视频欧美日本| videossex国产| 国产精品乱码一区二三区的特点| 在线观看美女被高潮喷水网站| 久久这里有精品视频免费| 国产精品,欧美在线| 国产精品蜜桃在线观看| 自拍偷自拍亚洲精品老妇| 国产淫片久久久久久久久| 亚洲成色77777| 嫩草影院新地址| kizo精华| 成年女人永久免费观看视频| 亚洲三级黄色毛片| 亚洲av成人av| 十八禁国产超污无遮挡网站| 美女黄网站色视频| 国产 一区 欧美 日韩| 亚洲精品成人久久久久久| 99热全是精品| 久久精品人妻少妇| 在线免费十八禁| 色网站视频免费| 69av精品久久久久久| 国产精品人妻久久久久久| 长腿黑丝高跟| 一级爰片在线观看| 亚州av有码| 高清在线视频一区二区三区 | 久久精品国产亚洲网站| 久久欧美精品欧美久久欧美| av女优亚洲男人天堂| 精品国产一区二区三区久久久樱花 | 超碰av人人做人人爽久久| 一级黄片播放器| 国产精品99久久久久久久久| 婷婷六月久久综合丁香| 熟女人妻精品中文字幕| 色5月婷婷丁香| 精品一区二区免费观看| 简卡轻食公司| 精品少妇黑人巨大在线播放 | 久久精品国产亚洲av天美| 日本黄色视频三级网站网址| 日韩强制内射视频| 国产91av在线免费观看| 身体一侧抽搐| 激情 狠狠 欧美| 亚洲精品一区蜜桃| 麻豆久久精品国产亚洲av| 日韩成人av中文字幕在线观看| 国产在线一区二区三区精 | 麻豆乱淫一区二区| 日本五十路高清| 99在线视频只有这里精品首页| 欧美日韩综合久久久久久| av视频在线观看入口| 亚洲无线观看免费| 午夜亚洲福利在线播放| 国产精品乱码一区二三区的特点| 国产黄色小视频在线观看| 国产欧美日韩精品一区二区| 亚洲av成人精品一二三区| 一级av片app| 亚洲成人精品中文字幕电影| 最近的中文字幕免费完整| 久久久久久九九精品二区国产| 欧美xxxx黑人xx丫x性爽| 黄色配什么色好看| 91久久精品国产一区二区成人| 久久久久国产网址| 免费搜索国产男女视频| 老司机影院毛片| 国产精品久久久久久av不卡| 小说图片视频综合网站| 禁无遮挡网站| 午夜福利网站1000一区二区三区| 亚洲av免费高清在线观看| 免费无遮挡裸体视频| 日韩av在线大香蕉| 久久精品久久精品一区二区三区| 日韩视频在线欧美| 亚洲国产欧美人成| 91精品一卡2卡3卡4卡| 91久久精品国产一区二区三区| 看片在线看免费视频| 国产一区二区亚洲精品在线观看| 久久久久久久久中文| 色综合站精品国产| 日日摸夜夜添夜夜添av毛片| 天美传媒精品一区二区| 免费av观看视频| 男女国产视频网站| 午夜免费男女啪啪视频观看| 免费大片18禁| 欧美日本视频| 久久久久久大精品| 国产亚洲av片在线观看秒播厂 | 麻豆国产97在线/欧美| 1024手机看黄色片| 亚洲在线自拍视频| 亚洲精品,欧美精品| 视频中文字幕在线观看| videos熟女内射| 久久精品人妻少妇| 亚洲高清免费不卡视频| 精品国产露脸久久av麻豆 | 国产美女午夜福利| 亚洲真实伦在线观看| 少妇熟女欧美另类| 亚洲成人精品中文字幕电影| 欧美极品一区二区三区四区| 久久久久久伊人网av| 最近最新中文字幕大全电影3| 日日撸夜夜添| 99热这里只有是精品50| 成年女人永久免费观看视频| 日韩在线高清观看一区二区三区| 国产黄片视频在线免费观看| 不卡视频在线观看欧美| 美女内射精品一级片tv| 九九在线视频观看精品| 亚洲国产色片| 汤姆久久久久久久影院中文字幕 | 精品人妻一区二区三区麻豆| 亚洲av福利一区| 成人欧美大片| 欧美一区二区亚洲| 国产精品电影一区二区三区| 看非洲黑人一级黄片| 日韩 亚洲 欧美在线| 全区人妻精品视频| 欧美成人一区二区免费高清观看| 少妇熟女aⅴ在线视频| 精品午夜福利在线看| 老司机福利观看| 黄片无遮挡物在线观看| 亚洲国产色片| 欧美日韩综合久久久久久| 午夜亚洲福利在线播放| 亚洲精华国产精华液的使用体验| 中文精品一卡2卡3卡4更新| 国产69精品久久久久777片| 美女脱内裤让男人舔精品视频| 国产视频首页在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产 一区精品| 午夜激情福利司机影院| 国产极品精品免费视频能看的| 日本av手机在线免费观看| 久久久精品欧美日韩精品| 国产成人aa在线观看| 亚洲国产精品合色在线| 欧美zozozo另类| 亚洲成人av在线免费| 国产精品熟女久久久久浪| 日韩人妻高清精品专区| 色视频www国产| 91av网一区二区| av线在线观看网站| АⅤ资源中文在线天堂| 免费av观看视频| 日日干狠狠操夜夜爽| 在线观看66精品国产| 日本一二三区视频观看| 亚洲激情五月婷婷啪啪| 亚洲欧美精品综合久久99| 亚洲伊人久久精品综合 | 天天躁夜夜躁狠狠久久av| 少妇人妻精品综合一区二区| 国产午夜精品久久久久久一区二区三区| 嘟嘟电影网在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久久a久久爽久久v久久| 又粗又硬又长又爽又黄的视频| 18+在线观看网站| 一级毛片久久久久久久久女| 一级黄色大片毛片| 久久久色成人| 亚洲自拍偷在线| 久久精品国产亚洲av天美| 国产精品一区二区三区四区久久| 中国美白少妇内射xxxbb| 亚洲国产精品成人综合色| 亚洲电影在线观看av| 亚洲国产精品久久男人天堂| 亚洲一级一片aⅴ在线观看| 国产精品野战在线观看| 日日干狠狠操夜夜爽| 中文亚洲av片在线观看爽| 能在线免费观看的黄片| 日本午夜av视频| 国产又黄又爽又无遮挡在线| 国产亚洲最大av| 免费看光身美女| 亚洲va在线va天堂va国产| 日韩一本色道免费dvd| 日本黄色片子视频| 久久久久久久久大av| 久久精品熟女亚洲av麻豆精品 | 一级毛片我不卡| 午夜福利网站1000一区二区三区| 男人狂女人下面高潮的视频| 亚洲一级一片aⅴ在线观看| 小蜜桃在线观看免费完整版高清| 99国产精品一区二区蜜桃av| 国产精品一区二区性色av| 桃色一区二区三区在线观看| 99热网站在线观看| 小蜜桃在线观看免费完整版高清| 91久久精品国产一区二区成人| 永久网站在线| 成人鲁丝片一二三区免费| 亚洲综合色惰| 中文字幕av在线有码专区| 久久99热这里只频精品6学生 | 欧美成人a在线观看| 成人毛片60女人毛片免费| 干丝袜人妻中文字幕| 高清日韩中文字幕在线| 波多野结衣巨乳人妻| 91久久精品国产一区二区成人| 大话2 男鬼变身卡| 久久精品国产99精品国产亚洲性色| 婷婷六月久久综合丁香| 看片在线看免费视频| 欧美zozozo另类| 人妻夜夜爽99麻豆av| 成人二区视频| 国产高清视频在线观看网站| 欧美日韩国产亚洲二区| 最近手机中文字幕大全| 男人舔奶头视频| 亚洲人成网站在线播| 一边亲一边摸免费视频| 伊人久久精品亚洲午夜| 国产精品嫩草影院av在线观看| 亚洲真实伦在线观看| 又爽又黄a免费视频| 欧美性猛交黑人性爽| 黄色配什么色好看| 亚洲欧美中文字幕日韩二区| 色噜噜av男人的天堂激情| 青春草视频在线免费观看| 亚洲av一区综合| 日本黄色视频三级网站网址| 干丝袜人妻中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 欧美人与善性xxx| 午夜免费激情av| 日产精品乱码卡一卡2卡三| 国产av一区在线观看免费| 九九热线精品视视频播放| h日本视频在线播放| 日韩 亚洲 欧美在线| 亚洲av成人精品一二三区| 天堂av国产一区二区熟女人妻| 亚洲国产精品成人综合色| 午夜福利高清视频| 亚洲图色成人| 久久久精品欧美日韩精品| 美女脱内裤让男人舔精品视频| 国产乱来视频区| 国产精品,欧美在线| 国产成人一区二区在线| 男插女下体视频免费在线播放| 蜜臀久久99精品久久宅男| videos熟女内射| 一级毛片久久久久久久久女| 日本黄大片高清| 成人美女网站在线观看视频| 深爱激情五月婷婷| 亚洲中文字幕一区二区三区有码在线看| 桃色一区二区三区在线观看| 亚洲怡红院男人天堂| 午夜免费男女啪啪视频观看| 国语对白做爰xxxⅹ性视频网站| 午夜福利在线在线| 真实男女啪啪啪动态图| 国产精品女同一区二区软件| 身体一侧抽搐| 国产91av在线免费观看| 男女国产视频网站| 中文亚洲av片在线观看爽| 又爽又黄无遮挡网站| 长腿黑丝高跟| 久久久久久久久久久免费av| 亚洲国产精品成人久久小说| 热99在线观看视频| 波多野结衣巨乳人妻| 丰满少妇做爰视频| 人人妻人人澡欧美一区二区| 观看免费一级毛片| 91在线精品国自产拍蜜月| 欧美一区二区精品小视频在线| 日本三级黄在线观看| 亚洲av电影在线观看一区二区三区 | 最近最新中文字幕大全电影3| 少妇裸体淫交视频免费看高清| 成年女人看的毛片在线观看| 免费看日本二区| 床上黄色一级片| 在线观看66精品国产| 久久精品久久久久久久性| 国产亚洲91精品色在线| 国语自产精品视频在线第100页| 久久精品国产亚洲网站| 国产精品人妻久久久影院| 亚洲最大成人中文| 久久精品夜色国产| 国产精品国产三级国产av玫瑰| 亚洲最大成人中文| 免费黄色在线免费观看| 免费不卡的大黄色大毛片视频在线观看 | 国产午夜精品一二区理论片| 国产黄色小视频在线观看| 在线a可以看的网站| 久久久亚洲精品成人影院| 非洲黑人性xxxx精品又粗又长| 亚洲av中文av极速乱| 99热6这里只有精品| 国产亚洲精品av在线| 亚洲一级一片aⅴ在线观看| 天堂av国产一区二区熟女人妻| 少妇裸体淫交视频免费看高清| 免费av不卡在线播放| 性插视频无遮挡在线免费观看| 久久精品91蜜桃| 热99re8久久精品国产| 少妇丰满av| av在线老鸭窝| 色综合站精品国产| 人妻制服诱惑在线中文字幕| 禁无遮挡网站| 亚洲精品,欧美精品| 国产av一区在线观看免费| 有码 亚洲区| 日韩,欧美,国产一区二区三区 | 欧美+日韩+精品| 国语对白做爰xxxⅹ性视频网站| 亚洲精品乱码久久久久久按摩|