• <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超视频在线观看视频| 国产三级在线视频| 亚洲欧美精品自产自拍| 国语自产精品视频在线第100页| av在线亚洲专区| 国产亚洲精品av在线| 国产色爽女视频免费观看| 国产一区有黄有色的免费视频 | 久久久久精品久久久久真实原创| 最新中文字幕久久久久| 精品久久久久久电影网 | 超碰97精品在线观看| 亚洲人成网站在线播| 色尼玛亚洲综合影院| 国产在视频线在精品| 国产免费又黄又爽又色| 精品国内亚洲2022精品成人| 啦啦啦观看免费观看视频高清| 波多野结衣巨乳人妻| 色5月婷婷丁香| 91精品一卡2卡3卡4卡| 人人妻人人看人人澡| 亚洲成av人片在线播放无| 日韩,欧美,国产一区二区三区 | 99国产精品一区二区蜜桃av| 卡戴珊不雅视频在线播放| 亚洲aⅴ乱码一区二区在线播放| videossex国产| 好男人在线观看高清免费视频| 国产黄色视频一区二区在线观看 | 欧美一区二区精品小视频在线| 日产精品乱码卡一卡2卡三| 日本色播在线视频| 国产成人免费观看mmmm| 久久99热这里只有精品18| 亚洲av免费在线观看| 99热网站在线观看| 一边亲一边摸免费视频| 色综合站精品国产| 久久99热这里只频精品6学生 | 久久久久性生活片| 18禁动态无遮挡网站| 国产精品永久免费网站| 啦啦啦观看免费观看视频高清| 午夜免费男女啪啪视频观看| 国产免费男女视频| 亚洲成色77777| 国产熟女欧美一区二区| 天天躁夜夜躁狠狠久久av| 国产 一区精品| 搡女人真爽免费视频火全软件| 日日撸夜夜添| 丰满乱子伦码专区| 一夜夜www| 丰满少妇做爰视频| 99在线视频只有这里精品首页| 麻豆乱淫一区二区| 一卡2卡三卡四卡精品乱码亚洲| 男女那种视频在线观看| 国产又色又爽无遮挡免| 免费看光身美女| 国产精品久久久久久久电影| 麻豆久久精品国产亚洲av| 亚洲精品久久久久久婷婷小说 | 性色avwww在线观看| 亚洲最大成人手机在线| 国产美女午夜福利| 一级av片app| 久久精品国产亚洲av天美| 国产色爽女视频免费观看| 成人一区二区视频在线观看| 久久精品国产亚洲av涩爱| 网址你懂的国产日韩在线| 国产女主播在线喷水免费视频网站 | 热99re8久久精品国产| 嘟嘟电影网在线观看| 蜜桃亚洲精品一区二区三区| 久久久a久久爽久久v久久| 赤兔流量卡办理| 精品少妇黑人巨大在线播放 | 最近中文字幕高清免费大全6| 性色avwww在线观看| 99热全是精品| 久久精品国产亚洲av天美| 一级毛片aaaaaa免费看小| 特级一级黄色大片| 亚洲电影在线观看av| 日韩av在线免费看完整版不卡| 久久国产乱子免费精品| 丰满人妻一区二区三区视频av| 日本猛色少妇xxxxx猛交久久| 男女边吃奶边做爰视频| av在线播放精品| 精品欧美国产一区二区三| 国产精品国产三级国产专区5o | 91aial.com中文字幕在线观看| 搡女人真爽免费视频火全软件| 亚洲内射少妇av| 欧美激情在线99| 国产精品久久久久久精品电影小说 | 一级毛片我不卡| 人妻系列 视频| 国产在线一区二区三区精 | 成人毛片a级毛片在线播放| 少妇人妻一区二区三区视频| 最新中文字幕久久久久| 有码 亚洲区| 欧美zozozo另类| 99在线视频只有这里精品首页| 特大巨黑吊av在线直播| 亚洲三级黄色毛片| 国产亚洲精品av在线| 九九在线视频观看精品| 久久久亚洲精品成人影院| 色综合亚洲欧美另类图片| 男女国产视频网站| av国产久精品久网站免费入址| 国产黄色视频一区二区在线观看 | 免费看a级黄色片| 插逼视频在线观看| 午夜福利在线在线| 在线观看一区二区三区| 亚洲av成人av| av视频在线观看入口| 老司机影院成人| 久久精品国产亚洲av天美| 国产精品不卡视频一区二区| 你懂的网址亚洲精品在线观看 | 国产视频内射| 国产精品麻豆人妻色哟哟久久 | 久久久久国产网址| 亚洲第一区二区三区不卡| 亚洲国产精品sss在线观看| АⅤ资源中文在线天堂| 国产免费视频播放在线视频 | 少妇人妻精品综合一区二区| 女的被弄到高潮叫床怎么办| 国产精品麻豆人妻色哟哟久久 | 黄色一级大片看看| 三级国产精品欧美在线观看| av播播在线观看一区| 久久久久久久久久黄片| 午夜亚洲福利在线播放| 18禁裸乳无遮挡免费网站照片| 女人被狂操c到高潮| 99热这里只有是精品50| 久久久久国产网址| 精品一区二区三区人妻视频| 日韩精品有码人妻一区| 亚洲最大成人av| 春色校园在线视频观看| 国产精品久久久久久久久免| 久久亚洲国产成人精品v| 亚洲第一区二区三区不卡| 最新中文字幕久久久久| 欧美性猛交╳xxx乱大交人| 男人的好看免费观看在线视频| 免费看av在线观看网站| 日韩视频在线欧美| 女人久久www免费人成看片 | 久久久国产成人免费| 别揉我奶头 嗯啊视频| 亚洲丝袜综合中文字幕| 国产精品一区二区三区四区免费观看| 男的添女的下面高潮视频| 国产伦精品一区二区三区视频9| 久久久精品欧美日韩精品| av国产久精品久网站免费入址| 国产亚洲一区二区精品| 亚洲人成网站高清观看| 国产成人freesex在线| 国产女主播在线喷水免费视频网站 | 欧美变态另类bdsm刘玥| 欧美日韩综合久久久久久| 国产一区二区在线观看日韩| 床上黄色一级片| av在线播放精品| 国产视频内射| 亚洲人成网站在线观看播放| 国产免费男女视频| 日韩欧美 国产精品| 99久久精品热视频| 青青草视频在线视频观看| 99在线视频只有这里精品首页| 久久精品国产亚洲av天美| 亚洲精华国产精华液的使用体验| 国产老妇女一区| 哪个播放器可以免费观看大片| 国产一区二区在线av高清观看| 日韩欧美精品免费久久| 免费看av在线观看网站| 日韩一区二区视频免费看| 小说图片视频综合网站| 国产午夜精品久久久久久一区二区三区| 国产亚洲午夜精品一区二区久久 | eeuss影院久久| 精品久久久久久久久av| 日韩高清综合在线| 亚洲人成网站在线播| 少妇熟女欧美另类| 国产成人午夜福利电影在线观看| 日本黄大片高清| 午夜a级毛片| 深爱激情五月婷婷| 午夜激情欧美在线| 中文字幕久久专区| 亚洲电影在线观看av| 热99在线观看视频| 亚洲av一区综合| 午夜激情欧美在线| 欧美潮喷喷水| 日本免费在线观看一区| 亚洲一级一片aⅴ在线观看| 一边摸一边抽搐一进一小说| 99久国产av精品| 一个人看视频在线观看www免费| 日产精品乱码卡一卡2卡三| 中文字幕熟女人妻在线| 亚洲成人久久爱视频| 天堂网av新在线| 欧美最新免费一区二区三区| 久久久久精品久久久久真实原创| 天美传媒精品一区二区| 三级国产精品片| 女人被狂操c到高潮| 成人一区二区视频在线观看| 免费av观看视频| 精品国产一区二区三区久久久樱花 | 国产精品久久久久久av不卡| 国产精品一区www在线观看| 亚洲在线自拍视频| 色网站视频免费| 亚洲怡红院男人天堂| 国产精品久久视频播放| 亚洲av电影不卡..在线观看| 久久久久久久久久成人| 亚洲精品日韩在线中文字幕| 精品午夜福利在线看| 亚洲精品国产av成人精品| 国产精品永久免费网站| 久久久久国产网址| 成人三级黄色视频| 偷拍熟女少妇极品色| 99热这里只有是精品在线观看| 一级毛片aaaaaa免费看小| 日韩欧美在线乱码| 最近最新中文字幕大全电影3| 日韩欧美 国产精品| 99久国产av精品| 一级爰片在线观看| 大香蕉久久网| 黄片wwwwww| 国产免费视频播放在线视频 | 久久草成人影院| 国产精品国产高清国产av| 亚洲精品成人久久久久久| 我要看日韩黄色一级片| 免费在线观看成人毛片| 精品久久久久久久人妻蜜臀av| 日韩一本色道免费dvd| 久久久久久久午夜电影| 亚洲欧美中文字幕日韩二区| 午夜精品一区二区三区免费看| 一级黄色大片毛片| 草草在线视频免费看| 99久久无色码亚洲精品果冻| 日本三级黄在线观看| 色综合色国产| 日韩中字成人| av卡一久久| 精品人妻一区二区三区麻豆| 国产高清不卡午夜福利| 亚洲国产精品专区欧美| 天天躁日日操中文字幕| 国产在线一区二区三区精 | 最近最新中文字幕免费大全7| 国语对白做爰xxxⅹ性视频网站| 亚洲精品一区蜜桃| 日本免费在线观看一区| 黑人高潮一二区| 99久国产av精品国产电影| 一级毛片aaaaaa免费看小| 久久热精品热| 国产男人的电影天堂91| 日韩强制内射视频| 韩国av在线不卡| 国产又色又爽无遮挡免| 日韩人妻高清精品专区| 91aial.com中文字幕在线观看| 91在线精品国自产拍蜜月| 六月丁香七月| 亚洲欧美精品综合久久99| 黄色一级大片看看| 亚洲美女视频黄频| 丰满少妇做爰视频| 丝袜喷水一区| 边亲边吃奶的免费视频| 久久久久久大精品| 久久精品久久久久久久性| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品,欧美精品| 日本av手机在线免费观看| 国产单亲对白刺激| 女人久久www免费人成看片 | 桃色一区二区三区在线观看| 国产在线男女| 蜜臀久久99精品久久宅男| 99国产精品一区二区蜜桃av| 蜜桃久久精品国产亚洲av| 中文亚洲av片在线观看爽| 国产成年人精品一区二区| 国产老妇伦熟女老妇高清| 看黄色毛片网站| 日产精品乱码卡一卡2卡三| 日本av手机在线免费观看| 久久久久久久久久久丰满| 97超碰精品成人国产| 纵有疾风起免费观看全集完整版 | 只有这里有精品99| 国产精品精品国产色婷婷| 又爽又黄无遮挡网站| 国产精品一区二区三区四区免费观看| 久久精品久久久久久噜噜老黄 | 白带黄色成豆腐渣| 菩萨蛮人人尽说江南好唐韦庄 | 国产免费又黄又爽又色| 午夜福利网站1000一区二区三区| 国产精品爽爽va在线观看网站| 久久草成人影院| 小说图片视频综合网站| 国产一区二区在线av高清观看| 亚洲一区高清亚洲精品| 熟妇人妻久久中文字幕3abv| 黄片wwwwww| 男人的好看免费观看在线视频| 免费人成在线观看视频色| 午夜福利在线在线| 久久6这里有精品| 久久人妻av系列| 久久热精品热| 日本爱情动作片www.在线观看| 亚洲天堂国产精品一区在线| 我的老师免费观看完整版| 91久久精品国产一区二区三区| 国产一级毛片在线| 亚洲av二区三区四区| 亚洲欧洲日产国产| 久久精品国产99精品国产亚洲性色| 亚洲国产高清在线一区二区三| 亚洲av成人av| 在线免费观看不下载黄p国产| 欧美激情国产日韩精品一区| 色综合亚洲欧美另类图片| 亚洲欧美一区二区三区国产| 日本熟妇午夜| 99热网站在线观看| 国产精品一区www在线观看| 精品少妇黑人巨大在线播放 | 中文字幕免费在线视频6| 亚洲欧美精品自产自拍| 国产精品精品国产色婷婷| 成年av动漫网址| 纵有疾风起免费观看全集完整版 | 免费观看a级毛片全部| 国产亚洲精品久久久com| 欧美激情国产日韩精品一区| 99久久中文字幕三级久久日本| av女优亚洲男人天堂| 国产探花极品一区二区| 日日撸夜夜添| 久久精品国产亚洲网站| 亚洲成人av在线免费| 麻豆国产97在线/欧美| 免费看日本二区| 精品不卡国产一区二区三区| 天堂√8在线中文| 国产精品蜜桃在线观看| 日本黄色视频三级网站网址| or卡值多少钱| 午夜福利成人在线免费观看| 寂寞人妻少妇视频99o| 超碰av人人做人人爽久久| 亚洲国产精品国产精品| 在线观看av片永久免费下载| 嫩草影院新地址| 亚洲人成网站高清观看| 2021天堂中文幕一二区在线观| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站在线观看播放| 少妇丰满av| 午夜福利视频1000在线观看| av在线亚洲专区| 亚洲av.av天堂| 精品人妻熟女av久视频| 高清视频免费观看一区二区 | 日韩精品有码人妻一区| 亚洲av免费高清在线观看| 国产午夜精品论理片| 在线免费观看的www视频| 日韩一区二区三区影片| 伦理电影大哥的女人| 久久久精品94久久精品| 直男gayav资源| 亚洲av中文av极速乱| 一本久久精品| 人人妻人人看人人澡| 美女国产视频在线观看| 两个人视频免费观看高清| 高清毛片免费看| 中文字幕av成人在线电影| 欧美激情在线99| 国产精品野战在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产老妇伦熟女老妇高清| www.av在线官网国产| 男女边吃奶边做爰视频| 免费黄网站久久成人精品| 中文字幕av成人在线电影| 午夜老司机福利剧场| 国产伦理片在线播放av一区| 国产精品一区二区三区四区久久| 又爽又黄无遮挡网站| 天堂中文最新版在线下载 | 国产精品蜜桃在线观看| 国产高潮美女av| 亚洲高清免费不卡视频| 波多野结衣巨乳人妻| 丰满人妻一区二区三区视频av| 国产精品伦人一区二区| 纵有疾风起免费观看全集完整版 | 日日啪夜夜撸| 亚洲五月天丁香| 亚洲性久久影院| 亚洲国产欧美人成| 在线播放国产精品三级| 免费看a级黄色片| 日本一二三区视频观看| av在线亚洲专区| АⅤ资源中文在线天堂| 欧美成人免费av一区二区三区| 国内少妇人妻偷人精品xxx网站| 亚洲五月天丁香| 中国美白少妇内射xxxbb| 床上黄色一级片| 国产亚洲av片在线观看秒播厂 | 欧美日本视频| 中文精品一卡2卡3卡4更新| 亚洲图色成人| 国产亚洲5aaaaa淫片| 99久国产av精品国产电影| 亚洲国产精品合色在线| 久久久久久大精品| 国产一区有黄有色的免费视频 | 国产伦理片在线播放av一区| 国产精品国产三级专区第一集| 亚洲最大成人手机在线| 成年版毛片免费区| 综合色丁香网| 免费人成在线观看视频色| 亚洲一级一片aⅴ在线观看| 内地一区二区视频在线| 日本色播在线视频| 成年免费大片在线观看| 久久精品人妻少妇| 1024手机看黄色片| 日本黄色视频三级网站网址| 精品久久久久久久末码| 黄色一级大片看看| 国产黄色小视频在线观看| 伊人久久精品亚洲午夜| 狠狠狠狠99中文字幕| 成人亚洲欧美一区二区av| 简卡轻食公司| www日本黄色视频网| 久久人妻av系列| 久久久精品94久久精品| 女人久久www免费人成看片 | 国产午夜精品久久久久久一区二区三区| 欧美另类亚洲清纯唯美| 久久人人爽人人爽人人片va| 搞女人的毛片| 99久久精品国产国产毛片| 亚洲精品自拍成人| 中文字幕亚洲精品专区| 国产亚洲91精品色在线| 搞女人的毛片| 大香蕉97超碰在线| 亚洲精品自拍成人| 久久精品夜色国产| 日韩精品有码人妻一区| 国产精品电影一区二区三区| 久久精品久久精品一区二区三区| 国产午夜福利久久久久久| 国产精品国产三级国产专区5o | 欧美潮喷喷水| 国产高清有码在线观看视频| 日本黄大片高清| 久久精品久久久久久久性| 国产成人91sexporn| 成人午夜精彩视频在线观看| 国产免费视频播放在线视频 | 色哟哟·www| 久久鲁丝午夜福利片| 2021少妇久久久久久久久久久| 久久99精品国语久久久| 中文字幕人妻熟人妻熟丝袜美| 国产精品日韩av在线免费观看| 黄色配什么色好看| 99在线视频只有这里精品首页| 国产一区二区三区av在线| 一卡2卡三卡四卡精品乱码亚洲| 中文精品一卡2卡3卡4更新| 欧美人与善性xxx| 亚洲va在线va天堂va国产| 2022亚洲国产成人精品| 级片在线观看| 一级二级三级毛片免费看| 成人亚洲欧美一区二区av| 非洲黑人性xxxx精品又粗又长| 久久亚洲精品不卡| 国产伦精品一区二区三区四那| 日日摸夜夜添夜夜爱| 午夜福利网站1000一区二区三区| 老司机福利观看| 又粗又硬又长又爽又黄的视频| 伊人久久精品亚洲午夜| 男人和女人高潮做爰伦理| 国产色爽女视频免费观看| 成人国产麻豆网| 国产欧美日韩精品一区二区| 国产免费视频播放在线视频 | 国产国拍精品亚洲av在线观看| 日韩国内少妇激情av| 久久久色成人| eeuss影院久久| 国产淫语在线视频| 卡戴珊不雅视频在线播放| 精品人妻偷拍中文字幕| 久久精品国产亚洲av天美| 九九热线精品视视频播放| 在线观看66精品国产| 免费在线观看成人毛片| 国产精品一及| 亚洲无线观看免费| 99久国产av精品国产电影| 97人妻精品一区二区三区麻豆| 国产免费又黄又爽又色| 女人十人毛片免费观看3o分钟| 日韩欧美 国产精品| 亚洲久久久久久中文字幕| 一级毛片电影观看 | 天美传媒精品一区二区| 麻豆久久精品国产亚洲av| 亚洲av不卡在线观看| 国产一区二区在线av高清观看| 看十八女毛片水多多多| 99热6这里只有精品| 高清日韩中文字幕在线| 麻豆成人av视频| 国产精品嫩草影院av在线观看| 亚洲精品成人久久久久久| 可以在线观看毛片的网站| 欧美成人a在线观看| 欧美又色又爽又黄视频| 久久久a久久爽久久v久久| 校园人妻丝袜中文字幕| 欧美精品一区二区大全| 别揉我奶头 嗯啊视频| 乱人视频在线观看| 久久精品91蜜桃| 18禁动态无遮挡网站| 国产精品日韩av在线免费观看| 老司机影院毛片| 桃色一区二区三区在线观看| 五月玫瑰六月丁香| 欧美不卡视频在线免费观看| 少妇人妻精品综合一区二区| 亚洲一级一片aⅴ在线观看| 国产成人freesex在线| 久99久视频精品免费| 国产精品久久久久久精品电影| 91狼人影院| 国产片特级美女逼逼视频| 午夜a级毛片| 国产精品不卡视频一区二区| 美女cb高潮喷水在线观看| 3wmmmm亚洲av在线观看| 麻豆一二三区av精品| 18禁裸乳无遮挡免费网站照片| 少妇人妻一区二区三区视频| 成人欧美大片| 一级毛片我不卡| 五月伊人婷婷丁香| 午夜免费男女啪啪视频观看| 免费观看的影片在线观看| 日韩中字成人| kizo精华| 欧美色视频一区免费| 日韩制服骚丝袜av| 亚洲在久久综合| 狂野欧美白嫩少妇大欣赏| 国产精品国产三级国产av玫瑰| 九草在线视频观看| 免费在线观看成人毛片| 亚洲av福利一区| 免费看光身美女| 日本欧美国产在线视频| 日产精品乱码卡一卡2卡三| 亚洲人成网站在线播| 国产av不卡久久| 亚洲最大成人中文| 精品欧美国产一区二区三| 少妇人妻一区二区三区视频| 久久久a久久爽久久v久久|