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

    西安典型段地面沉降分層標(biāo)觀測(cè)及數(shù)值模擬

    2014-01-20 14:20:54王啟耀彭建兵蔣臻蔚滕宏泉
    巖土力學(xué) 2014年11期
    關(guān)鍵詞:承壓水含水層西安

    王啟耀 ,彭建兵 ,蔣臻蔚 ,滕宏泉

    (1.長(zhǎng)安大學(xué) 建筑工程學(xué)院,西安 710061;2.長(zhǎng)安大學(xué) 地質(zhì)工程與測(cè)繪學(xué)院,西安 710054 3.長(zhǎng)安大學(xué) 西部礦產(chǎn)資源與地質(zhì)工程教育部重點(diǎn)實(shí)驗(yàn)室,西安 710054;4.陜西省地質(zhì)環(huán)境監(jiān)測(cè)總站,西安 710054)

    1 引 言

    20 世紀(jì),隨著工業(yè)化、城市化的快速發(fā)展,對(duì)地下水需求的急劇增加導(dǎo)致了全球多個(gè)地區(qū)大范圍的地面沉降,誘發(fā)了大量的環(huán)境問題。由于特殊的地質(zhì)環(huán)境,不少地區(qū)的地面沉降還伴隨有不同程度的地裂縫,產(chǎn)生了比一般單純的地面沉降更為嚴(yán)重的災(zāi)害問題[1-3]。西安是我國最為典型的因過度抽水作用產(chǎn)生較大地面沉降和地裂縫的地區(qū),災(zāi)害尤為嚴(yán)重[4-5]。

    自1924年日本研究確定地面沉降與抽水具有明確的相關(guān)性之后,地面沉降的研究得到極大的關(guān)注,1969~2010年間共召開了8 屆地面沉降國際討論會(huì),進(jìn)行了專門的交流和探討。在抽水地面沉降的機(jī)制、數(shù)值計(jì)算、災(zāi)害防治等方面開展了廣泛的研究,取得了很大的成果,并就抽水地裂縫的機(jī)制提出各種模式[6-8]。近年來,較多地關(guān)注逐漸轉(zhuǎn)移到先存斷裂對(duì)地面沉降和地裂縫的影響,Bell 等[9]、Sheng 等[10]、Helm 分析了美國亞利桑那州以及拉斯維加斯發(fā)育的地裂縫,在考慮水動(dòng)力作用的基礎(chǔ)上,結(jié)合已經(jīng)存在的地質(zhì)構(gòu)造建立了不同的成因模式[10]。Martin 等[11-12]利用Abaqus 軟件附帶的擴(kuò)展有限元法(XFEM)模塊模擬了斷層帶附近地裂縫的產(chǎn)生。蔣臻蔚等[13]也對(duì)這個(gè)問題開展了一些初步的研究。

    由于地裂縫涉及到斷裂的起裂、擴(kuò)展或不連續(xù)面的張開、滑移問題,目前大范圍的地面沉降三維模擬計(jì)算還難以解決這類復(fù)雜問題[14-15]。本文基于西安地下水的監(jiān)測(cè)和較為詳細(xì)的地層剖面資料,采用比奧固結(jié)理論和不連續(xù)面的接觸分析,建立了包含2 條地裂縫的典型地段的二維抽水沉降計(jì)算模型,對(duì)抽水作用下先存斷裂存在下的地面沉降進(jìn)行了盡可能精細(xì)的模擬計(jì)算,結(jié)合該計(jì)算段的監(jiān)測(cè)資料,分析了抽水作用下地層壓縮量垂向上的分布特點(diǎn)及地面沉降水平位置上存在差異的原因。

    2 西安地面沉降分層標(biāo)監(jiān)測(cè)

    西安的地面沉降觀測(cè)可追溯到1959年,但地裂縫的定點(diǎn)監(jiān)測(cè)和短水準(zhǔn)剖面監(jiān)測(cè)在1989年之后才展開。1989年陜西環(huán)境監(jiān)測(cè)總站在西安地質(zhì)技校內(nèi)的f7地裂縫兩側(cè),建立了西安市惟一的地面沉降分層標(biāo)觀測(cè)點(diǎn)(見圖1)。地面沉降分層標(biāo)埋深為97.31~367.29 m,每季度末定期監(jiān)測(cè),對(duì)比分析土層變形量與地裂縫儀器監(jiān)測(cè)的活動(dòng)量,掌握地裂縫南、北兩側(cè)367 m 以上土層的變形特征及地裂縫南、北兩側(cè)最深標(biāo)底所測(cè)不受抽汲地下水影響土層之間的高差變化量及其隨時(shí)間的變化規(guī)律,至今繼續(xù)保留使用[16]。

    各分層標(biāo)的年度和累積活動(dòng)量如圖2、3 所示。1991年以來的地面沉降分層標(biāo)監(jiān)測(cè)結(jié)果表明,西安地面沉降強(qiáng)度有逐年減弱的總體趨勢(shì),2002年左右進(jìn)入基本穩(wěn)定期。各組分層標(biāo)沉降時(shí)間曲線總體形態(tài)極為相似,但沉降量變化較大,表明不同深度范圍地層的壓密量不同。淺層地層0~104.90 m(S2標(biāo))沉降量變化幅度很小,說明淺層的潛水段地層受深部承壓水水位降低影響很小。埋深104.90~259.91 m(S4標(biāo))段是西安市承壓水主要開采段。從1997年來,由于承壓水開采量逐年減小,水頭降幅減小,部分地段水頭還有所回升,使得分層標(biāo)近年來沉降量總體回落,埋深259.91 m 的S4標(biāo)和367.28 m 的S1標(biāo)由于受地表各種因素影響較小,沉降量逐年減小趨勢(shì)尤其明顯。

    圖1 西安地面沉降分層標(biāo)設(shè)置示意圖Fig.1 Layerwise marks of land subsidence in Xi’an

    圖2 分層標(biāo)年度活動(dòng)量Fig.2 Annual movement of the Layerwise marks

    圖3 分層標(biāo)各標(biāo)累計(jì)活動(dòng)量Fig.3 Accumulated movement of Layerwise marks

    截至目前,地裂縫南側(cè)洼地的259.91~367.28 m段(第3 承壓含水層)地層累計(jì)壓密量為98.58 mm;187.28~259.91 m 段(第2 承壓含水層)地層累計(jì)壓密量為248.50 mm;104.90~187.28 m 段(第1承壓含水層)地層累計(jì)壓密量為525.49 mm;0~104.90 m段(潛水層0地層累計(jì)壓密量為148.85 mm。潛水層、第1~3 承壓含水層的壓密量占總壓密量百分比分別為14.9%、51.4%、24.3%和9.6%。由以上分析可以得出,104.90~187.28 m(第1 承壓含水層)地層為主要壓密層,187.28~259.91 m 段(第2 承壓含水層)為次壓密層。

    地裂縫北側(cè)黃土梁區(qū)的地層沉降量(N 標(biāo),0~360 m)與南側(cè)洼地S1標(biāo)(0~367 m)相比要小得多,截止到目前,北標(biāo)N 累計(jì)沉降量為415.68 mm,南標(biāo)S1累計(jì)沉降量為1 021.42 mm,兩者相差了605.74 mm。由于南北兩側(cè)水位變化不大,按照以往的解釋,這個(gè)差異完全是由南北兩側(cè)地層上的差異造成的,但根據(jù)設(shè)置分層標(biāo)時(shí)的鉆孔資料,兩側(cè)的地層差異并不顯著,那么造成這種巨大差異沉降的原因只能是斷層,但斷層的存在對(duì)地面沉降的影響真有如此強(qiáng)烈嗎?

    3 基于庫侖接觸單元的沉降數(shù)值模擬

    3.1 數(shù)值模擬基本原理和技術(shù)

    進(jìn)行先存斷裂面時(shí)的抽水致裂模擬計(jì)算必須解決2個(gè)關(guān)鍵問題,一個(gè)是滲流固結(jié)耦合問題,另一個(gè)是地裂縫模擬問題。

    地裂縫的出現(xiàn)與地層的水平運(yùn)動(dòng)關(guān)系較大,因此滲流固結(jié)計(jì)算必須采用可以考慮水平變形的固結(jié)理論,本文采用Biot 固結(jié)理論。

    從探槽揭露及試驗(yàn)的情況看,西安地裂縫主要是一個(gè)接觸面的問題,即兩盤土體沿地裂縫發(fā)生張開和滑移的問題,因此,本文采用FLAC2D中的接觸面單元來模擬地裂縫。

    FLAC2D中的接觸面采用的是基于A、B 面接觸判斷的庫侖滑動(dòng)和張裂模型。該接觸面模型的好處是與實(shí)際情況較為吻合,每個(gè)參數(shù)都具有明確的物理意義,能夠較好地模擬地裂縫的開啟及垂直錯(cuò)動(dòng)。

    3.2 計(jì)算模型與參數(shù)

    本次計(jì)算在垂直地裂縫的方向截取包括f6、f72 條地裂縫的祭臺(tái)村-地質(zhì)技校段進(jìn)行計(jì)算。計(jì)算段裂縫間距為1 200 m,兩側(cè)各取800 m,總計(jì)算長(zhǎng)度為2 800 m。計(jì)算深度為420~440 m 不等,頂部地形的刻畫與實(shí)際一致,保持梁洼結(jié)構(gòu)的特點(diǎn)。地層從上到下包括Q3、Q2黃土、Q2沖湖積粉質(zhì)黏土夾砂土、Q1沖洪積粉質(zhì)黏土夾砂土,具體情況基本按照實(shí)際剖面確定。

    計(jì)算模型為平面應(yīng)變模型,邊界條件:地面為自由位移,自由排水邊界;底部為固定位移,不排水邊界;兩側(cè)邊為水平固定,豎向自由,不排水邊界。由于實(shí)際抽水情況的復(fù)雜性,加上本文并不是真正的三維計(jì)算,很難模擬真正的抽水情況,所以計(jì)算中以水位控制作為各個(gè)時(shí)間段的抽水強(qiáng)度指標(biāo)。為了得到與監(jiān)測(cè)數(shù)據(jù)較為一致的水位,模型中設(shè)置了多個(gè)抽水井,通過調(diào)整每個(gè)井的抽水強(qiáng)度,使得不同時(shí)間段的地下水位基本與實(shí)際吻合。最后建立的計(jì)算模型見圖4(為了顯示剖面上的變化,圖形在橫向上進(jìn)行了壓縮)。

    西安400 m 深度范圍內(nèi)具有3 層承壓水,但根據(jù)監(jiān)測(cè),3 層承壓水聯(lián)通較好,具有一致的承壓水位,因此計(jì)算中將承壓水統(tǒng)一考慮。為了模擬承壓水,把潛水層底部的黏土層設(shè)置為不透水層,其上的初始水壓根據(jù)潛水位計(jì)算賦值,其下的承壓水水位根據(jù)1960年初始監(jiān)測(cè)得到的數(shù)據(jù)賦值。

    圖4 西安地面沉降地裂縫計(jì)算模型圖Fig.4 Calculating model of subsidence and ground fissures in Xi’an

    由于西安的地面沉降從60年代就開始出現(xiàn),水位資料從1960年有了較詳細(xì)的記錄[16]。本次計(jì)算從1960年開始,滲流計(jì)算時(shí)的初始水頭根據(jù)1960年的地下水承壓水位圖確定,從左到右為393~405 m,后面的變化根據(jù)相應(yīng)的水位動(dòng)態(tài)曲線確定,如圖5 所示。從水位變化圖上看,計(jì)算段承壓水頭東南高、西北低,水位歷史最大降深為80 m。

    模型計(jì)算參數(shù)主要包括土層物理力學(xué)性質(zhì)參數(shù)和地裂縫的參數(shù),土層參數(shù)主要根據(jù)陜西省環(huán)境監(jiān)測(cè)總站在西安地質(zhì)技工學(xué)校地裂縫監(jiān)測(cè)工作時(shí)的鉆孔資料確定,見表1。

    圖5 計(jì)算初始的1960年承壓水位(單位:m)Fig.5 The water level of confined groundwater on 1960(unit:m)

    表1 模型計(jì)算土層參數(shù)Table 1 Soils parameters of calculation model

    圖6 計(jì)算剖面上水位動(dòng)態(tài)曲線Fig.6 Dynamic curves of confined groundwater level within calculating section

    地裂縫參數(shù)根據(jù)資料分析和經(jīng)驗(yàn)綜合確定,剪切剛度 Ks=1×107MPa/m,法向剛度Kn=1×108MPa/m,黏聚力c=0 MPa,內(nèi)摩擦角φ=10°,抗拉強(qiáng)度T=0 MPa。

    3.3 數(shù)值計(jì)算結(jié)果與分析

    根據(jù)實(shí)際土層、地貌和水位動(dòng)態(tài)情況計(jì)算得到的地面沉降和地裂縫的發(fā)展結(jié)果如圖7、8 所示。

    由計(jì)算結(jié)果可以發(fā)現(xiàn),地面沉降和地裂縫的發(fā)展與地下水的開采具有很好的對(duì)應(yīng)關(guān)系。圖7、8顯示,地裂縫在1960~1971年間由于地下水開采規(guī)模較小,地裂縫發(fā)展緩慢,垂直位錯(cuò)f6地裂縫僅為2.5 cm,f7地裂縫為1.2 cm,地裂縫活動(dòng)不明顯;但隨著地下水開采強(qiáng)度的提高,1971~1976年間地下水位急劇下降,地面沉降和地裂縫也迅速發(fā)展,地裂縫垂直錯(cuò)距達(dá)16.0 cm(f6地裂縫)和18.1 cm(f7地裂縫),地裂縫災(zāi)害非常明顯;隨后的1977~1997年間,地下水開采稍有減緩,但地下水位仍然穩(wěn)定地下降,直到部分位置水位降落到第1 層隔水底板之下,地裂縫和地面沉降繼續(xù)發(fā)展;1998年之后,由于地下水開采減少,水位逐漸回升,但到2002年左右,地面沉降和地裂縫年度活動(dòng)量才顯著減小,說明地面沉降與水位變化相比具有一定的滯后性。

    圖7 地面沉降地裂縫隨時(shí)間發(fā)展曲線Fig.7 Vertical displacement duration curves of land subsidence and ground fissures

    圖8 各年地表形態(tài)變化圖(放大50 倍)Fig.8 Patterns of surface configuration over the years(magnified 50 times)

    計(jì)算得到的0~100、100~300、300~400 m深度的壓縮量如圖9 所示。由圖可見,西安抽水引起的地層壓縮主要集中在100~300 m,截止2008年,該段累積壓縮量為140.8 cm,占了總壓縮量的82.2%,而上部的潛水層及300 m 以下的深部含水層壓縮量都較小,其中100 m 以上地層累積為19.9 cm,占總壓縮量的11.6%,300 m 以下地層累積為10.5 cm,占總壓縮量的6.2%,與分層標(biāo)實(shí)測(cè)資料得到的結(jié)果較為接近。

    地面沉降初期,地裂縫活動(dòng)不明顯時(shí)地裂縫南北兩側(cè)的沉降量差異不大,隨著沉降的進(jìn)一步發(fā)展,地裂縫活動(dòng)越來越強(qiáng)烈,其兩側(cè)梁洼兩地的沉降差異越來越大,如果僅比較1989年以后的沉降,計(jì)算得到的南側(cè)洼地沉降為72 cm,而北側(cè)梁部的沉降僅為34 cm,遠(yuǎn)小于南側(cè)洼地,與監(jiān)測(cè)情況較為一致。

    如果不改變計(jì)算模型的其他情況,只改變斷層的強(qiáng)度參數(shù)內(nèi)摩擦角的大小,計(jì)算得到的北側(cè)梁部、南側(cè)洼地和地裂縫的垂直活動(dòng)量如圖10 所示。

    圖9 分層累積壓縮量圖Fig.9 Layered compression over the years

    圖10 地面沉降和地裂縫活動(dòng)與斷層面強(qiáng)度關(guān)系Fig.10 Relationships between land subsidence and activity of ground fissure

    由圖10 可見,斷裂對(duì)地面沉降具有顯著的影響,由于斷裂的存在,上盤(梁部)地面沉降被減小了,而下盤(洼區(qū))的沉降得到放大,并且這種影響隨著斷層強(qiáng)度的降低而更為明顯。如果斷裂的強(qiáng)度參數(shù)與周圍圍巖接近,則類似于沒有斷裂,地裂縫也不出現(xiàn),兩側(cè)的沉降也會(huì)趨于接近。

    4 結(jié) 論

    (1)西安市的潛水層、第1~3 承壓含水層的壓密量占總壓密量百分比分別為14.9%、51.4%、24.3%和9.6%。104.90~187.28 m(第1 承壓含水層)地層為最主要壓密層。

    (2)西安市的地面沉降和地裂縫與抽水有直接的對(duì)應(yīng)關(guān)系,其發(fā)展變化分為1959~1971年的緩慢發(fā)展階段,1971~1976年的加劇發(fā)展階段,1976~2002年的穩(wěn)定發(fā)展階段,2002 以后的減速階段。從其發(fā)展過程來看,在沒有大的地質(zhì)構(gòu)造活動(dòng)的情況下,只要今后控制好地下水的開采,地面沉降和地裂縫就不會(huì)有大的發(fā)展變化。

    (3)地裂縫南北梁洼部位的沉降量差異非常明顯,已經(jīng)不能簡(jiǎn)單地按照傳統(tǒng)的地層差異沉降來解釋,先存斷裂在這里起了很大的作用,由于斷裂的存在以及其正斷層的特點(diǎn),沉降過程中上盤(南盤)的沉降得到了放大,而下盤(北盤)沉降縮小。正反的結(jié)果導(dǎo)致地面沉降在斷裂位置出現(xiàn)差異放大現(xiàn)象,即產(chǎn)生地裂縫。數(shù)值計(jì)算的結(jié)果與監(jiān)測(cè)情況較為吻合,很好地解釋了這個(gè)現(xiàn)象。

    需要注意的是,盡管本次計(jì)算模型盡量與實(shí)際情況相吻合,但由于地質(zhì)情況的復(fù)雜性,抽水情況的多變性,該問題本質(zhì)上是個(gè)三維問題,二維計(jì)算的結(jié)果與實(shí)際監(jiān)測(cè)在數(shù)值上并不完全一致,但基本的變化過程和規(guī)律較為一致,說明二維計(jì)算在對(duì)復(fù)雜情況下抽水地面沉降問題的機(jī)制研究方面還是非常合適的。

    [1]HOLZER T L.Preconsolidation stress of aquifer systems in areas of induced land subsidence[J].Water Resources Research,1981,17(3):693-704.

    [2]SCHUMANN H H.,CRIPE L S.Land subsidence and earth fissures caused by groundwater depletion in southern Arizona,U.S.A.[C]//Proceedings of the 3rd International Symposium on Land Subsidence,Venice.Wallingford,UK:International Association of Hydrological Science Press,1986:841-851.

    [3]HERMAN B.Land subsidence and cracking due to groundwater depletion[J].Ground Water,1977,15(5):358-364

    [4]張家明.西安地裂縫研究[M].西安:西北大學(xué)出版社,1990.

    [5]王景明.地裂縫及其災(zāi)害的理論與應(yīng)用[M].西安:陜西科學(xué)技術(shù)出版社,2000.

    [6]CARPENTER M C.South-Central Arizona:Earth fissures and subsidence complicate development of desert water resources[C]//Land Subsidence in the United States.U.S:Geological Survey Circular,1999:65-78.

    [7]HOLZER T L.Ground failure in areas of subsidence due to ground-water decline in the United States[C]//Proc.2nd Internat Symp.On Land Subsidence,Internat.Assoc.Sci.Hydrol.Anaheim,Calif.:[s.n.],1984:13-17.

    [8]HOLZER T L,DAVIS S N,LOFGREN B E.Faulting caused by groundwater extraction in South central Arizona[J].Journal of Geophysical Research,1979,84(2):603-612.

    [9]BELL J W,PRICE J G,MIFFLIN M D.Subsidence induced fissuring along preexisting faults in Las Vegas Valley,Nevada[C]//Proceedings,Association of Engineering Geologists,35th Annual Meeting.Los Angeles:[s.n.],1992:66-75.

    [10]SHENG Z P,DONALD C H,JIANG Li.Mechanisms of earth fissuring caused by groundwater withdrawal[J].Environmental &Engineering Geoscience,2003,Ⅸ(4):351-362.

    [11]HERNANDEZ-MARIN M,THOMAS J B.Faultcontrolled deformation and stress from pumping-induced groundwater flow[J].Journal of Hydrology,2012,428-429:80-93

    [12]HERNANDEZ-MARIN M,THOMAS J B.Controls on initiation and propagation of pumping-induced earth fissures:Insights from numerical simulations[J].Hydro-Geology Journal,2010,18(6):1773-1785

    [13]王啟耀,蔣臻蔚,彭建兵.抽水作用下先存斷裂活化滑移機(jī)制研究[J].水文地質(zhì)工程地質(zhì),2013,40(2):108-112.WANG Qi-yao,JIANG Zhen-wei,PENG Jian-bing.Mechanism of reactivation and slip of the preexisting fault under pumping[J].Hydrogeology &Engineering Geology,2013,40(2):108-112.

    [14]薛禹群.我國地面沉降模擬現(xiàn)狀及需要解決的問題[J].水文地質(zhì)工程地質(zhì),2003,30(5):1-5.XUE Yu-qun.Present status of modeling land subsidence in China and problems to be solved[J].Hydrogeology and Engineering Geology,2003,30(5):1-5.

    [15]張?jiān)?薛禹群.抽水地面沉降數(shù)學(xué)模型的研究現(xiàn)狀與展望[J].中國地質(zhì)災(zāi)害與防治學(xué)報(bào),2002,13(2):1-7.ZHANG Yun,XUE Yu-qun.Present situation and prospect on the mathematical model of land subsidence due to pumping[J].The Chinese Journal of Geological Hazard and Control,2002,13(2):1-7.

    [16]彭建兵,張勤,黃強(qiáng)兵,等.汾渭盆地典型地區(qū)地裂縫地面沉降監(jiān)測(cè)與防治[R].西安:長(zhǎng)安大學(xué),2009.

    猜你喜歡
    承壓水含水層西安
    全球多個(gè)含水層里的水正快速流失
    地鐵深基坑承壓水控制研究
    西安2021
    Oh 西安
    深層承壓水污染途徑及防治研究
    《西安人的歌》突如其來?
    追根溯源 回到西安
    美國西部奧加拉拉含水層水位下降原因初探
    全球地下含水層下降驚人:要被抽干了
    巖溶含水層水流模型研究進(jìn)展
    97人妻精品一区二区三区麻豆| 欧美极品一区二区三区四区| 99热这里只有精品一区| 国产精品电影一区二区三区| 国产老妇女一区| 久久精品91蜜桃| 国产精品无大码| 国产成年人精品一区二区| 偷拍熟女少妇极品色| 一区二区三区免费毛片| 国产一级毛片在线| 国产午夜精品论理片| 人妻夜夜爽99麻豆av| 99热网站在线观看| 嫩草影院入口| 国产久久久一区二区三区| 欧美激情国产日韩精品一区| 中文亚洲av片在线观看爽| 91午夜精品亚洲一区二区三区| 国产成年人精品一区二区| 亚洲婷婷狠狠爱综合网| av黄色大香蕉| 在线观看免费视频日本深夜| 色综合亚洲欧美另类图片| 国产高清有码在线观看视频| 如何舔出高潮| 人妻制服诱惑在线中文字幕| 精华霜和精华液先用哪个| 综合色av麻豆| 国产淫片久久久久久久久| 老师上课跳d突然被开到最大视频| 伊人久久精品亚洲午夜| 哪个播放器可以免费观看大片| 久久久精品94久久精品| 日韩亚洲欧美综合| 亚洲丝袜综合中文字幕| 男女啪啪激烈高潮av片| 变态另类丝袜制服| 亚洲人与动物交配视频| 少妇的逼水好多| 国内精品美女久久久久久| 久久婷婷人人爽人人干人人爱| 我要搜黄色片| 91精品国产九色| 97超碰精品成人国产| 一进一出抽搐gif免费好疼| 高清日韩中文字幕在线| 亚洲精品久久国产高清桃花| 女人被狂操c到高潮| 国产v大片淫在线免费观看| 国产单亲对白刺激| 欧美3d第一页| 男插女下体视频免费在线播放| 国产精品爽爽va在线观看网站| 欧美日韩综合久久久久久| 人妻夜夜爽99麻豆av| 啦啦啦观看免费观看视频高清| 欧美色视频一区免费| 亚洲18禁久久av| 自拍偷自拍亚洲精品老妇| 男人舔奶头视频| 国产乱人视频| 亚洲乱码一区二区免费版| 男人舔女人下体高潮全视频| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久久久久精品电影| 国产精品,欧美在线| 亚洲欧美成人综合另类久久久 | 网址你懂的国产日韩在线| 婷婷亚洲欧美| 亚洲在久久综合| 久久鲁丝午夜福利片| 亚洲最大成人中文| 日日摸夜夜添夜夜添av毛片| 亚洲人成网站高清观看| 一个人看视频在线观看www免费| 插逼视频在线观看| 精品午夜福利在线看| 能在线免费看毛片的网站| 亚洲美女搞黄在线观看| 永久网站在线| a级毛片免费高清观看在线播放| 成人性生交大片免费视频hd| 91精品一卡2卡3卡4卡| 日韩av不卡免费在线播放| 欧美一区二区国产精品久久精品| 久久精品久久久久久久性| 亚洲精华国产精华液的使用体验 | 中出人妻视频一区二区| 看十八女毛片水多多多| 国产精品国产三级国产av玫瑰| 亚洲精品国产av成人精品| av黄色大香蕉| 亚洲精品自拍成人| 中文字幕熟女人妻在线| 国产成人午夜福利电影在线观看| 免费av毛片视频| 色哟哟哟哟哟哟| 久久亚洲国产成人精品v| 成人鲁丝片一二三区免费| 性插视频无遮挡在线免费观看| 国内揄拍国产精品人妻在线| 寂寞人妻少妇视频99o| 国产色婷婷99| 一进一出抽搐动态| 91在线精品国自产拍蜜月| 欧美在线一区亚洲| 亚洲欧美中文字幕日韩二区| 99久国产av精品| 久久人妻av系列| 插逼视频在线观看| 成年女人永久免费观看视频| 干丝袜人妻中文字幕| 亚洲欧美日韩东京热| 国产一区二区三区av在线 | 国产伦在线观看视频一区| 日韩一本色道免费dvd| 久久精品综合一区二区三区| 日本黄色片子视频| 亚洲av二区三区四区| 国产蜜桃级精品一区二区三区| 亚洲婷婷狠狠爱综合网| 美女被艹到高潮喷水动态| 成人毛片a级毛片在线播放| 哪里可以看免费的av片| 国产单亲对白刺激| 国产黄片美女视频| 中文精品一卡2卡3卡4更新| 精品久久久久久久末码| 乱码一卡2卡4卡精品| 久久精品91蜜桃| 99热只有精品国产| 美女 人体艺术 gogo| 亚洲欧美成人精品一区二区| 国内揄拍国产精品人妻在线| 男女视频在线观看网站免费| 亚洲精品乱码久久久v下载方式| 久久久久久久久久久免费av| 国产成人福利小说| 欧美成人免费av一区二区三区| 一个人看的www免费观看视频| 亚洲美女视频黄频| 简卡轻食公司| 成人高潮视频无遮挡免费网站| 国产高清激情床上av| 97热精品久久久久久| 久久久a久久爽久久v久久| 两个人的视频大全免费| 欧美日韩乱码在线| 亚洲人成网站在线播| 夫妻性生交免费视频一级片| 国产一区二区在线av高清观看| 少妇熟女欧美另类| 久久亚洲国产成人精品v| 国产一区二区在线观看日韩| 在现免费观看毛片| 日韩一区二区视频免费看| 精品人妻偷拍中文字幕| 九色成人免费人妻av| 99热6这里只有精品| 国产伦在线观看视频一区| 人人妻人人澡欧美一区二区| 色播亚洲综合网| 久久久久久久久久久丰满| 国产亚洲精品久久久久久毛片| 国产精品免费一区二区三区在线| 麻豆乱淫一区二区| 久久鲁丝午夜福利片| 久久99热这里只有精品18| 亚洲欧美清纯卡通| 精品久久久久久久人妻蜜臀av| 免费看a级黄色片| 黄色日韩在线| 看非洲黑人一级黄片| 最近手机中文字幕大全| 国产在线男女| 中文欧美无线码| 99久久无色码亚洲精品果冻| 欧美变态另类bdsm刘玥| 一区福利在线观看| 成人性生交大片免费视频hd| 久久久久久久久久久免费av| 青春草亚洲视频在线观看| 欧美极品一区二区三区四区| 国内精品美女久久久久久| 国产视频首页在线观看| 国产av一区在线观看免费| 两个人的视频大全免费| 69av精品久久久久久| 日韩欧美国产在线观看| 2022亚洲国产成人精品| 亚洲精品久久国产高清桃花| 久久久久免费精品人妻一区二区| 我的女老师完整版在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区三区四区免费观看| АⅤ资源中文在线天堂| 欧美成人免费av一区二区三区| 亚洲第一区二区三区不卡| 97在线视频观看| 嘟嘟电影网在线观看| 小说图片视频综合网站| 噜噜噜噜噜久久久久久91| 亚洲av电影不卡..在线观看| 黄色一级大片看看| 亚洲av男天堂| 日韩高清综合在线| 一本精品99久久精品77| 免费无遮挡裸体视频| 久久婷婷人人爽人人干人人爱| 日韩三级伦理在线观看| 国产不卡一卡二| 狂野欧美激情性xxxx在线观看| 1024手机看黄色片| 中文字幕人妻熟人妻熟丝袜美| av在线亚洲专区| 中文字幕av成人在线电影| 日韩欧美三级三区| 给我免费播放毛片高清在线观看| 久久精品国产99精品国产亚洲性色| 国产午夜福利久久久久久| 国产久久久一区二区三区| av在线老鸭窝| 变态另类成人亚洲欧美熟女| 久99久视频精品免费| 蜜桃久久精品国产亚洲av| 在线观看66精品国产| 男人舔奶头视频| 久久综合国产亚洲精品| 性欧美人与动物交配| 国产精品一区二区三区四区免费观看| 99热网站在线观看| 男人的好看免费观看在线视频| 日韩欧美国产在线观看| 免费av不卡在线播放| 国产极品天堂在线| 一夜夜www| 精品久久久久久久久久久久久| 九九在线视频观看精品| 久久久精品94久久精品| 欧美成人免费av一区二区三区| 亚洲国产精品合色在线| 婷婷六月久久综合丁香| 成人无遮挡网站| 午夜精品国产一区二区电影 | 国产精品久久久久久精品电影小说 | 国产在线精品亚洲第一网站| 精品熟女少妇av免费看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 中国美女看黄片| 国产视频内射| 两性午夜刺激爽爽歪歪视频在线观看| 日韩精品有码人妻一区| videossex国产| 久久精品国产亚洲av天美| 欧美一级a爱片免费观看看| av.在线天堂| 丰满乱子伦码专区| 国内揄拍国产精品人妻在线| 毛片一级片免费看久久久久| 最近最新中文字幕大全电影3| 综合色丁香网| 午夜激情欧美在线| 亚洲欧美成人精品一区二区| 久久精品人妻少妇| 亚洲av第一区精品v没综合| 国产久久久一区二区三区| 亚洲中文字幕日韩| 久久精品综合一区二区三区| 婷婷色av中文字幕| 三级男女做爰猛烈吃奶摸视频| 亚洲18禁久久av| 精品无人区乱码1区二区| 成人高潮视频无遮挡免费网站| 亚洲av.av天堂| 成人av在线播放网站| 看黄色毛片网站| 亚洲精品日韩在线中文字幕 | 人妻久久中文字幕网| 中文资源天堂在线| 波多野结衣高清无吗| 99九九线精品视频在线观看视频| 在线免费十八禁| 国产精品野战在线观看| 亚洲av中文av极速乱| 久久久久久久午夜电影| 日日干狠狠操夜夜爽| 校园春色视频在线观看| 欧美日韩国产亚洲二区| 男人舔奶头视频| 激情 狠狠 欧美| 成年免费大片在线观看| 在线播放无遮挡| 久久人人精品亚洲av| av视频在线观看入口| 校园春色视频在线观看| 亚洲电影在线观看av| 亚洲精品粉嫩美女一区| 国产亚洲av片在线观看秒播厂 | 欧美人与善性xxx| 床上黄色一级片| 精品人妻一区二区三区麻豆| 亚洲欧美成人综合另类久久久 | 国产成人freesex在线| 哪里可以看免费的av片| 国产日韩欧美在线精品| 国产私拍福利视频在线观看| 久久久久久久久久久丰满| 免费看a级黄色片| 国产成人aa在线观看| 亚洲国产高清在线一区二区三| 免费搜索国产男女视频| 国产精品国产高清国产av| or卡值多少钱| 日韩在线高清观看一区二区三区| 国产一区二区亚洲精品在线观看| 免费观看精品视频网站| 国产精品一区二区三区四区久久| 国产高清三级在线| 亚洲av男天堂| 国产av一区在线观看免费| 26uuu在线亚洲综合色| 麻豆精品久久久久久蜜桃| 精品一区二区三区人妻视频| 在线观看午夜福利视频| 变态另类成人亚洲欧美熟女| 97热精品久久久久久| 九草在线视频观看| 中文字幕熟女人妻在线| 精品99又大又爽又粗少妇毛片| 久久国内精品自在自线图片| 亚洲久久久久久中文字幕| 白带黄色成豆腐渣| av免费在线看不卡| 老司机影院成人| 欧美变态另类bdsm刘玥| 日韩精品青青久久久久久| 免费搜索国产男女视频| 黄色欧美视频在线观看| 波多野结衣高清无吗| 99国产精品一区二区蜜桃av| 最近最新中文字幕大全电影3| 日韩精品青青久久久久久| 久久6这里有精品| 蜜桃亚洲精品一区二区三区| 国产熟女欧美一区二区| 国产一区二区在线av高清观看| 一级av片app| 熟女电影av网| 美女xxoo啪啪120秒动态图| 欧洲精品卡2卡3卡4卡5卡区| 不卡一级毛片| 我要看日韩黄色一级片| 麻豆一二三区av精品| av又黄又爽大尺度在线免费看 | 国产精品一二三区在线看| 国产精品综合久久久久久久免费| 久久久久九九精品影院| 在线a可以看的网站| 99九九线精品视频在线观看视频| av在线亚洲专区| 亚洲最大成人av| 亚洲丝袜综合中文字幕| 床上黄色一级片| 日韩视频在线欧美| 亚洲国产欧美人成| 久久久欧美国产精品| 国产男人的电影天堂91| av在线老鸭窝| 国产成人福利小说| 午夜视频国产福利| 成人特级av手机在线观看| 久久这里有精品视频免费| 丰满人妻一区二区三区视频av| 尤物成人国产欧美一区二区三区| 99久久中文字幕三级久久日本| 亚洲不卡免费看| 日韩精品有码人妻一区| 日韩欧美国产在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品国产成人久久av| 国产精品麻豆人妻色哟哟久久 | 午夜a级毛片| 国产一级毛片在线| 欧美日本亚洲视频在线播放| 日本黄色片子视频| 97人妻精品一区二区三区麻豆| 日本黄色片子视频| 欧美日韩精品成人综合77777| 嘟嘟电影网在线观看| 精品人妻视频免费看| 一进一出抽搐动态| 亚洲自拍偷在线| 桃色一区二区三区在线观看| 国产精品无大码| 亚洲国产色片| 久久欧美精品欧美久久欧美| 三级毛片av免费| 久久欧美精品欧美久久欧美| 色视频www国产| 两个人的视频大全免费| 亚洲高清免费不卡视频| 在线观看66精品国产| 亚洲国产精品成人综合色| 赤兔流量卡办理| 26uuu在线亚洲综合色| 国产精品蜜桃在线观看 | 欧美性感艳星| av在线亚洲专区| 成人漫画全彩无遮挡| 久久这里有精品视频免费| 国产精品日韩av在线免费观看| 麻豆一二三区av精品| 直男gayav资源| 午夜a级毛片| 亚洲人成网站在线播| 草草在线视频免费看| 亚洲精品乱码久久久v下载方式| 美女国产视频在线观看| 国产av一区在线观看免费| 高清日韩中文字幕在线| 国产av麻豆久久久久久久| 久久久久久久午夜电影| 亚洲欧美成人精品一区二区| 日本av手机在线免费观看| 岛国毛片在线播放| 精品不卡国产一区二区三区| 天天躁日日操中文字幕| 99视频精品全部免费 在线| 最后的刺客免费高清国语| 九九在线视频观看精品| 亚洲国产精品成人久久小说 | 中文字幕免费在线视频6| 少妇丰满av| 亚洲图色成人| 亚洲国产欧美在线一区| 精品久久久久久久久久久久久| 日本五十路高清| 国产精品一区www在线观看| 国产中年淑女户外野战色| 欧美激情国产日韩精品一区| 床上黄色一级片| 亚洲天堂国产精品一区在线| 黄片无遮挡物在线观看| 最近的中文字幕免费完整| 日韩欧美国产在线观看| 国产在线男女| 久久久久久久久大av| 蜜桃亚洲精品一区二区三区| 91午夜精品亚洲一区二区三区| 18禁裸乳无遮挡免费网站照片| 国产精品乱码一区二三区的特点| 一级毛片电影观看 | 日韩成人av中文字幕在线观看| 日韩欧美国产在线观看| 3wmmmm亚洲av在线观看| 精品一区二区三区视频在线| 免费av不卡在线播放| 精品免费久久久久久久清纯| 少妇的逼水好多| 丰满人妻一区二区三区视频av| 中出人妻视频一区二区| 亚洲精品乱码久久久久久按摩| 国产亚洲av嫩草精品影院| 国产白丝娇喘喷水9色精品| 欧美性猛交╳xxx乱大交人| 国产色婷婷99| 国产精品一及| 亚洲精品乱码久久久v下载方式| 国产午夜精品久久久久久一区二区三区| 麻豆国产av国片精品| 夫妻性生交免费视频一级片| 最近中文字幕高清免费大全6| kizo精华| 亚洲内射少妇av| 国产在线精品亚洲第一网站| www.色视频.com| 老司机福利观看| 国产精品无大码| 国产日本99.免费观看| 在线观看av片永久免费下载| 精品人妻熟女av久视频| 日韩人妻高清精品专区| 国产精品嫩草影院av在线观看| 小说图片视频综合网站| 国产高清激情床上av| 一级毛片电影观看 | 亚洲熟妇中文字幕五十中出| 免费看日本二区| 中文字幕av成人在线电影| 高清毛片免费观看视频网站| 成人永久免费在线观看视频| 国产成人福利小说| 观看免费一级毛片| 九九久久精品国产亚洲av麻豆| 中出人妻视频一区二区| 国产极品天堂在线| 精品人妻熟女av久视频| 小说图片视频综合网站| 亚洲国产精品成人综合色| 精品少妇黑人巨大在线播放 | 国产高清有码在线观看视频| 欧美变态另类bdsm刘玥| 天堂√8在线中文| 精品人妻一区二区三区麻豆| 99久久精品热视频| 校园春色视频在线观看| av在线观看视频网站免费| 干丝袜人妻中文字幕| 国产一区二区在线观看日韩| 观看免费一级毛片| 九九热线精品视视频播放| 天天一区二区日本电影三级| 亚洲va在线va天堂va国产| 亚洲精品乱码久久久久久按摩| 一级av片app| av免费在线看不卡| 精品久久久久久久人妻蜜臀av| 少妇裸体淫交视频免费看高清| 一级黄色大片毛片| 午夜激情欧美在线| 国产不卡一卡二| 精品久久久久久久久久免费视频| 一级毛片电影观看 | 亚洲色图av天堂| 国产亚洲av嫩草精品影院| 亚洲精品亚洲一区二区| 99久国产av精品| 国产精品,欧美在线| 欧美高清成人免费视频www| 亚洲av熟女| 欧美三级亚洲精品| 日本黄色片子视频| 丰满的人妻完整版| 精品久久久久久久人妻蜜臀av| 青青草视频在线视频观看| 婷婷精品国产亚洲av| 毛片女人毛片| 日本与韩国留学比较| 久久精品国产亚洲av天美| 国产一区二区激情短视频| 日本熟妇午夜| 国产精品嫩草影院av在线观看| 人人妻人人澡欧美一区二区| 直男gayav资源| 精品久久国产蜜桃| 久久久久国产网址| 好男人在线观看高清免费视频| 此物有八面人人有两片| 在线播放国产精品三级| 激情 狠狠 欧美| 在线a可以看的网站| 99热6这里只有精品| 菩萨蛮人人尽说江南好唐韦庄 | 两性午夜刺激爽爽歪歪视频在线观看| 欧美xxxx黑人xx丫x性爽| 久久热精品热| 国产久久久一区二区三区| 久久久精品94久久精品| 国产黄片美女视频| 三级毛片av免费| 日韩欧美一区二区三区在线观看| a级一级毛片免费在线观看| av视频在线观看入口| 亚洲一区二区三区色噜噜| 欧美3d第一页| 亚洲av免费在线观看| 一个人看视频在线观看www免费| 一本久久中文字幕| 国产亚洲5aaaaa淫片| 中文在线观看免费www的网站| 国产精品精品国产色婷婷| 全区人妻精品视频| 青春草国产在线视频 | 啦啦啦啦在线视频资源| 亚洲一级一片aⅴ在线观看| 中文资源天堂在线| 美女内射精品一级片tv| 国产亚洲av片在线观看秒播厂 | 久久精品国产自在天天线| 麻豆一二三区av精品| 一级毛片我不卡| 日韩欧美精品v在线| 啦啦啦韩国在线观看视频| 91久久精品国产一区二区成人| 男人舔奶头视频| 日本一本二区三区精品| 欧美zozozo另类| 久久久午夜欧美精品| 亚洲av成人精品一区久久| av视频在线观看入口| 此物有八面人人有两片| 国产在线精品亚洲第一网站| 亚洲aⅴ乱码一区二区在线播放| 在线观看av片永久免费下载| 国产在视频线在精品| 亚洲av第一区精品v没综合| 国产午夜福利久久久久久| 精品久久久久久成人av| 国产精品精品国产色婷婷| 成人毛片60女人毛片免费| 国产精品久久电影中文字幕| 激情 狠狠 欧美| 人妻夜夜爽99麻豆av| 国产精品永久免费网站| 一个人免费在线观看电影| 日韩av不卡免费在线播放| 国产高清三级在线| 久久韩国三级中文字幕| 亚洲欧美精品专区久久| 国产av麻豆久久久久久久| 高清毛片免费观看视频网站| 看非洲黑人一级黄片| 日日摸夜夜添夜夜添av毛片| 插阴视频在线观看视频|