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

    完整井抽降水引起的側(cè)向有界越流承壓含水層變形解析研究

    2020-08-14 07:13:06蔡兵華李忠超馮慶高
    安全與環(huán)境工程 2020年4期
    關(guān)鍵詞:無(wú)界有界水井

    蔡兵華,李忠超 ,馮慶高,袁 祥

    (1.武漢市市政建設(shè)集團(tuán)有限公司,湖北 武漢 430023;2.中國(guó)地質(zhì)大學(xué)(武漢)工程學(xué)院,湖北 武漢 430074)

    基坑工程中常采用降水井抽取承壓水的方式來(lái)降低地下水水位以滿足施工要求,但抽取承壓水不可避免地會(huì)引起地層沉降[1-5]。為了避免抽降承壓水對(duì)周邊建筑物或地下構(gòu)筑物造成的危害,須將基坑降水引起的地層變形控制在合理的范圍內(nèi)。因此,分析抽降承壓水引起的地層變形顯得尤為重要。

    針對(duì)基坑降水引起的地層變形問(wèn)題,許多學(xué)者進(jìn)行了大量研究。如Terzaghi[6]提出土的有效應(yīng)力原理,解釋了降水會(huì)引起地層中孔隙水壓力的減小致使其有效應(yīng)力增加,從而導(dǎo)致土層產(chǎn)生壓縮變形;Shen等[7]通過(guò)結(jié)合三維滲流和一維土體固結(jié)的方法,提出了關(guān)于地下水開(kāi)采引起地面沉降的計(jì)算方法;許燁霜等[8]在考慮了地層壓縮性的情況下,分析計(jì)算了降水引起的地面沉降問(wèn)題;Shen等[9]在考慮土體固結(jié)過(guò)程中壓縮系數(shù)和滲透系數(shù)變化的條件下,建立了三維滲流和一維固結(jié)的部分耦合模型,分析了開(kāi)采地下水對(duì)地面沉降的影響。但是更貼合實(shí)際的是應(yīng)考慮地下水滲流和土體變形的相互作用,使?jié)B流和應(yīng)力實(shí)現(xiàn)真正意義上的流固耦合。為此,有不少學(xué)者基于Biot固結(jié)理論來(lái)計(jì)算抽水所引起的土體變形。如Derski[10]最早利用Biot固結(jié)理論分析了抽水所引起的土體變形問(wèn)題;Booker等[11]利用Fourier積分變換,求解得到了彈性半空間體內(nèi)點(diǎn)源穩(wěn)態(tài)抽水作用下土體的Boit固結(jié)解析解;艾智勇等[12-13]利用傳遞矩陣法解決了飽和半空間或飽和成層地基中點(diǎn)源非穩(wěn)定抽水問(wèn)題;Selvadurai等[14]應(yīng)用Laplace-Hankel積分變換求解得到彈性半空間無(wú)限體內(nèi)某一圓盤抽水區(qū)域條件下,非穩(wěn)態(tài)抽水引起的土體固結(jié)變形的解析解。此外,還有諸多學(xué)者[15-19]對(duì)此問(wèn)題進(jìn)行了研究。

    現(xiàn)有的計(jì)算研究普遍建立在含水層是側(cè)向無(wú)限延伸(側(cè)向無(wú)界)假定的基礎(chǔ)上,缺少以在距離抽水井r=R(本文將R稱為有界半徑)處的地下水水位沒(méi)有變化為邊界條件的水位降深解答。而在工程實(shí)踐中,特別是基坑工程中,往往在距離抽水井r=R處存在補(bǔ)給源(如回灌源)使得該處的地下水水位不發(fā)生變化或在r=R處的地下水水位不受抽水的影響,這樣的含水系統(tǒng)被認(rèn)為是側(cè)向有限延伸的,被稱為側(cè)向有限含水系統(tǒng)。此外,降水通常是通過(guò)豎向抽水井來(lái)實(shí)現(xiàn)的,應(yīng)考慮井的完整性這一客觀事實(shí)。因此,有必要對(duì)不同側(cè)向邊界條件下完整井抽降承壓水引起的含水層變形進(jìn)行分析。

    在抽水過(guò)程中通常認(rèn)為完整井附近的地下水壓力和流速僅隨含水層的平面方向x、y或r、θ變化,而不隨垂直方向z變化,同時(shí)還具有軸對(duì)稱的特點(diǎn)?;谏鲜鐾暾樗奶攸c(diǎn),本文在平面應(yīng)力假定的基礎(chǔ)上[20-21],從Biot固結(jié)方程出發(fā),應(yīng)用Laplace變換求解了越流承壓含水層在側(cè)向有界的情況下,完整井抽降水引起的承壓含水層豎向變形和徑向變形的二維解析表達(dá),并重點(diǎn)分析了不同側(cè)向邊界對(duì)完整井抽降承壓水引起的含水層變形的影響。

    1 模型建立及求解

    1.1 基本方程

    側(cè)向有界抽水承壓含水層二維變形計(jì)算示意圖見(jiàn)圖1。側(cè)向有界越流承壓含水層中存在一個(gè)定流量Q抽水的完整井,抽水承壓含水層均質(zhì)、等厚,含水層內(nèi)的滲流僅沿徑向上發(fā)生。其中,抽水承壓含水層厚度為B、滲透系數(shù)為K、彈性模量為E、泊松比為ν,且抽水承壓含水層上覆一個(gè)滲透系數(shù)為K′、厚度為B′的弱透水層,在距離抽水井r=R處存在補(bǔ)給源使得該處地下水水位沒(méi)有變化。

    圖1 側(cè)向有界抽水承壓含水層二維變形計(jì)算示意圖Fig.1 Schematic diagram of two-dimensional deformation calculation for pumped confined aquifer of laterally finite extent in radial direction

    (1)

    (2)

    (3)

    式中:ur和uz分別為徑向位移和豎向位移。

    (4)

    (5)

    再將公式(5)對(duì)r積分,可得:

    (6)

    式中:f稱為無(wú)量綱積分常數(shù),由通過(guò)在距離抽水井r=R處的邊界條件來(lái)確定。

    結(jié)合方程(3)~(6)并化簡(jiǎn),可得到平面應(yīng)力狀態(tài)下土體的基本方程如下:

    (7)

    (8)

    (9)

    1.2 固結(jié)控制方程的求解

    考慮抽水承壓含水層上覆土層的越流補(bǔ)給,則土體質(zhì)量守恒方程為

    (10)

    式中:K和K′分別為抽水承壓含水層和上覆弱透水層的滲透系數(shù)(cm/s);B和B′分別為抽水承壓含水層和上覆弱透水層的厚度;γw為水的容重(kg/m3)。

    將方程(7)代入方程(10),可得:

    (11)

    完整井以定流量Q抽水時(shí)的井壁邊界條件為

    (12)

    且在r=0處徑向變形為零,即:

    ur(r→0,t)=0

    (13)

    同時(shí),在r=R處應(yīng)滿足如下條件:

    p(r=R,t)=0

    (14)

    σrr(r=R,t)=0

    (15)

    將方程(11)對(duì)時(shí)間t進(jìn)行Laplace變換,可得:

    (16)

    方程(16)的通解為

    (17)

    式中:I0為第一類0階修正貝塞爾函數(shù);K0為第二類0階修正貝塞爾函數(shù)。

    根據(jù)方程(12)和(14)可以確定方程(17)中的待定常數(shù)A1和A2,因此可以得到在Laplace空間下,孔隙水壓力的解析解為

    (18)

    隨后對(duì)方程(7)~(9)進(jìn)行Laplace變換,并將方程(18)代入,可得:

    (19)

    (20)

    (21)

    將公式(20)對(duì)r積分,并考慮徑向變形滿足方程(13)這一條件,可得:

    (22)

    式中:I1為第一類1階修正貝塞爾函數(shù);K1為第二類1階修正貝塞爾函數(shù)。

    假定抽水承壓含水層底部無(wú)豎向變形,則含水層頂部豎向變形可表達(dá)為

    (23)

    (24)

    在r=R時(shí),由方程(19)、(21)和(22),可得:

    (25)

    (26)

    (27)

    (28)

    其中:

    (29)

    (30)

    (31)

    (32)

    (33)

    至此,得到了拉氏空間下的半解析解答,由孔隙水壓力和位移的具體表達(dá)式(31)~(33)可知,很難對(duì)其直接進(jìn)行拉氏逆變換而得到時(shí)間域的解答,為此本文采用常用的Stehfest數(shù)值逆變換法[22]來(lái)求解實(shí)域解。此外,當(dāng)不考慮越流影響或R→∞時(shí),方程(31)~(33)可以簡(jiǎn)化為完全承壓含水層或側(cè)向無(wú)限含水層情況下孔隙水壓力和變形的表達(dá)式,具體在此不做詳述。

    2 分析與討論

    本文在上述得到的解析解的基礎(chǔ)上,首先通過(guò)有限元數(shù)值計(jì)算的方法對(duì)解析解進(jìn)行對(duì)比驗(yàn)證,然后重點(diǎn)分析越流承壓含水層側(cè)向有限延伸或無(wú)限延伸條件下,完整井抽降承壓水引起的含水層豎向變形和徑向變形的特征。在無(wú)特別說(shuō)明的情況下,抽水引起的含水層豎向變形和徑向變形均表示為正值。

    2.1 有限元計(jì)算的對(duì)比驗(yàn)證

    根據(jù)上述得到的完整井抽降承壓水引起的地層變形的解析解答,將本文計(jì)算結(jié)果與COMSOL Multiphysics的有限元計(jì)算結(jié)果進(jìn)行了對(duì)比。依據(jù)圖1建立了完整井抽降承壓水的簡(jiǎn)化二維有限元計(jì)算模型。模型的水平方向尺寸設(shè)為1 000 m,由于抽水承壓含水層上覆弱透水層的越流補(bǔ)給是通過(guò)設(shè)置越流項(xiàng)來(lái)實(shí)現(xiàn)的,故僅需在豎直方向建立厚度為20 m的承壓含水層即可;降水井作為線源位于模型中心位置,在線源上即進(jìn)水濾管處設(shè)置流量邊界以實(shí)現(xiàn)抽水過(guò)程,其中抽水流量Q=10 m3/h、抽水時(shí)間t=300 h;同時(shí),上覆弱透水層的越流補(bǔ)給量在COMSOL Multiphysics的滲流方程中設(shè)置添加越流項(xiàng)即可,其最終考慮越流補(bǔ)給的土體質(zhì)量平衡方程正如方程(10)一樣;此外,模型兩側(cè)設(shè)置水平方向的位移約束和定水頭補(bǔ)給邊界,模型底部設(shè)置水平和豎直兩個(gè)方向的位移約束,為不透水邊界,并在距離降水井兩側(cè)300 m位置處設(shè)置定水頭邊界,以保證抽水過(guò)程中該處地下水水位沒(méi)有變化。整個(gè)有限元計(jì)算模型共劃分為1 146個(gè)三角形單元,土體本構(gòu)模型選用線彈性模型。在有限元計(jì)算和解析計(jì)算中都采用統(tǒng)一的土體參數(shù)值,具體土體參數(shù)值見(jiàn)表1。

    表1 土體參數(shù)值Table 1 Values of soil parameters

    表1中的土體參數(shù)取值可以如實(shí)地反映出含水砂層和黏土弱透水層的基本物理力學(xué)特性,以下在對(duì)抽水承壓含水層變形的影響分析中土層的基本參數(shù)取值與表1中一致。本文計(jì)算結(jié)果與有限元計(jì)算結(jié)果的對(duì)比,見(jiàn)圖2。

    圖2 本文計(jì)算結(jié)果與有限元計(jì)算結(jié)果的對(duì)比Fig.2 Comparison of the calculation results obtained by the proposed solution and the finite element numerical solution

    由圖2可見(jiàn),利用本文計(jì)算方法得到的完整井抽降承壓水引起的含水層豎向位移(Uz)和徑向位移(Ur)與有限元方法的計(jì)算結(jié)果具有較好的一致性,驗(yàn)證了本文計(jì)算方法的正確可靠性。

    2.2 抽水承壓含水層豎向變形的影響分析

    在距抽水井距離r=25 m的位置處,側(cè)向無(wú)界和有界(有界半徑R=300 m)條件下越流承壓含水層中豎向位移Uz與抽水時(shí)間t的關(guān)系曲線,見(jiàn)圖3。

    圖3 側(cè)向無(wú)界和有界(R=300 m)條件下越流承壓 含水層Uz-t的關(guān)系曲線圖Fig.3 Relation curves of the vertical deformation (Uz) versus the pumping time (t) in leaky confined aquifer of infinite and finite extent in radial direction (R=300 m)

    由圖3可見(jiàn),在整個(gè)抽水階段中,無(wú)論承壓含水層側(cè)向有界還是無(wú)界,其豎向位移Uz都隨著抽水時(shí)間t的增長(zhǎng)而增大,并在后期達(dá)到穩(wěn)定;當(dāng)抽水時(shí)間t小于10 h時(shí),尚未達(dá)到滲流穩(wěn)定階段,抽水引起的側(cè)向有界和無(wú)界含水層的豎向位移基本是一致的,而在抽水后期,側(cè)向有界含水層相對(duì)來(lái)說(shuō)更易達(dá)到穩(wěn)定滲流階段,因此側(cè)向無(wú)界含水層的豎向變形要大于側(cè)向有界含水層。

    當(dāng)抽水時(shí)間t=300 h時(shí),側(cè)向無(wú)界和有界(R=200 m)條件下越流承壓含水層的豎向位移Uz與距抽水井距離r的關(guān)系曲線,見(jiàn)圖4。

    圖4 側(cè)向無(wú)界和有界(R=200 m)條件下越流承壓 含水層中Uz-r的關(guān)系曲線圖Fig.4 Relation curves of the vertical deformation (Uz) versus the distance (r) in leaky confined aquifer of infinite and finite extent in radial direction (R=200 m)

    由圖4可見(jiàn),從含水層土體變形的整體規(guī)律上看,無(wú)論越流承壓含水層是側(cè)向有界還是無(wú)界,抽水引起的含水層豎向位移Uz均隨著與抽水井距離r的增大而減小;在抽水滲流達(dá)到穩(wěn)定階段時(shí),側(cè)向無(wú)界含水層的豎向變形量要大于側(cè)向有界含水層,且這種變形量的差異在距離抽水井越遠(yuǎn)的位置越發(fā)明顯,這是由于越接近側(cè)向邊界(R=200 m),受到側(cè)向邊界補(bǔ)給作用的影響更加明顯,從而減小了抽水井對(duì)豎向變形的影響。

    當(dāng)有界半徑R取不同值時(shí),越流承壓含水層的豎向位移Uz見(jiàn)圖5,為了比較分析,側(cè)向無(wú)界含水層中抽水的情況也考慮其中。

    圖5 不同有界半徑R對(duì)越流承壓含水層豎向位移 Uz的影響Fig.5 Impact of finite radius (R) on the vertical deformation of leaky confined aquifer

    由圖5可見(jiàn),在較長(zhǎng)的抽水時(shí)間(t=300 h)內(nèi),隨著有界半徑R的增大,側(cè)向有界含水層的豎向位移Uz越接近于側(cè)向無(wú)界含水層;在抽水達(dá)到穩(wěn)定滲流階段后,隨著有界半徑R的增大,側(cè)向有界含水層的豎向變形范圍和具體量值均隨之增大,這是由于側(cè)向邊界越接近于抽水井,其補(bǔ)給作用對(duì)抽水的影響則越明顯,從而減小了抽水井所能引起的含水層豎向變形范圍和豎向變形量。

    2.3 抽水含水層徑向變形的影響分析

    在距抽水井距離r=25 m位置處,抽水引起的側(cè)向無(wú)界和有界(R=200 m)越流承壓含水層的徑向變形Ur與抽水時(shí)間t的關(guān)系曲線,見(jiàn)圖6。

    圖6 側(cè)向無(wú)界和有界(R=200 m)條件下越流承壓 含水層Ur-t的關(guān)系曲線Fig.6 Relation curves of the vertical deformation (Uz) versus the pumping time (t) in leaky confined aquifer of infinite and finite extent in radial direction (R=200 m)

    由圖6可見(jiàn),在整個(gè)抽水階段中,側(cè)向無(wú)界和有界越流承壓含水層的徑向變形Ur均隨著抽水時(shí)間t的增長(zhǎng)而增大,并在達(dá)到穩(wěn)定滲流時(shí)保持不變,且此時(shí)含水層徑向位移達(dá)到最大值;此外,在相同的抽水時(shí)間,距離抽水井同一位置處側(cè)向無(wú)界承壓含水層的徑向位移要小于側(cè)向有界承壓含水層,這是由于在側(cè)向有界承壓含水層中,側(cè)向邊界位置(有界半徑R)處的地下水水位始終保持不變,抽水過(guò)程中會(huì)形成一個(gè)相比于側(cè)向無(wú)界承壓含水層中更大的水頭差和滲流力作用,從而引起含水層更大的徑向變形。

    在不同抽水時(shí)間t=50 h和t=300 h時(shí),側(cè)向無(wú)界和有界(R=200 m)條件下越流承壓含水層的徑向位移Ur與距抽水井距離r的關(guān)系曲線,見(jiàn)圖7。

    圖7 不同抽水時(shí)間側(cè)向無(wú)界和有界條件下越流承壓 含水層Ur-r的關(guān)系曲線Fig.7 Relation curves of the radial deformation (Ur) versus the distance (r) in leaky confined aquifer of infinite and finite extent in radial direction under different pumping time

    由圖7可見(jiàn),在同一抽水時(shí)間,側(cè)向無(wú)界和有界承壓含水層的徑向位移Ur沿徑向的分布規(guī)律基本一致,即隨著距抽水井距離r越遠(yuǎn),側(cè)向有界和無(wú)界承壓含水層徑向位移Ur呈現(xiàn)先增大后減小的趨勢(shì),這是由于靠近抽水井位置附近的側(cè)向有界和無(wú)界承壓含水層處于壓縮變形區(qū),距離抽水井較遠(yuǎn)位置的側(cè)向有界和無(wú)界承壓含水層處于拉張變形區(qū),而側(cè)向有界和無(wú)界承壓含水層的最大徑向變形會(huì)出現(xiàn)在這兩個(gè)區(qū)域的交界位置處,因此越流承壓含水層的最大徑向變形出現(xiàn)在距抽水井一定距離處。

    當(dāng)有界半徑R取不同值時(shí),越流承壓含水層的徑向位移Ur見(jiàn)圖8,為了比較分析,側(cè)向無(wú)界含水層中抽水的情況也考慮其中。

    圖8 不同有界半徑R對(duì)越流承壓含水層徑向位移 Ur的影響Fig.8 Impact of finite radius (R) on the radial deformation of leaky confined aquifer(Ur)

    由圖8可見(jiàn),在較長(zhǎng)的抽水時(shí)間(t=300 h)內(nèi),隨著有界半徑R的增大,抽水井引起的側(cè)向有界承壓含水層的徑向變形量和變形范圍均隨之增大,這是由于有界半徑R越大,側(cè)向邊界對(duì)抽水的影響越小,抽水所能引起的含水層徑向變形量和變形范圍也就越大;同時(shí),隨著側(cè)向邊界與抽水井相距越遠(yuǎn),含水層的拉張變形區(qū)也相對(duì)距離抽水井越遠(yuǎn),進(jìn)而導(dǎo)致含水層壓縮變形區(qū)和拉張變形區(qū)的交界位置也越遠(yuǎn)離抽水井,最終呈現(xiàn)如圖8所示含水層最大徑向變形位置隨著有界半徑R的增大而遠(yuǎn)離抽水井這一特征;此外,側(cè)向無(wú)界條件下含水層的徑向變形規(guī)律與側(cè)向有界條件有所差異,這是由于側(cè)向邊界的補(bǔ)給作用減小了抽水井對(duì)含水層徑向變形的影響范圍,從而導(dǎo)致側(cè)向無(wú)界承壓含水層的徑向變形范圍要大于側(cè)向有界承壓含水層。

    3 結(jié) 論

    本文建立并得到了完整井抽降水引起的側(cè)向有界越流承壓含水層豎向位移和徑向位移的解析解,并重點(diǎn)分析了完整井抽水條件下側(cè)向邊界對(duì)越流承壓含水層豎向變形和徑向變形的影響,得到主要結(jié)論如下:

    (1) 無(wú)論抽水承壓含水層側(cè)向有界還是無(wú)界,在整個(gè)抽水階段,承壓含水層的豎向變形在抽水前期均隨著抽水時(shí)間的增長(zhǎng)而增大,并在后期達(dá)到穩(wěn)定保持不變,抽水引起的承壓含水層豎向變形也隨著距抽水井距離的增大而減??;在達(dá)到滲流穩(wěn)定階段后,側(cè)向無(wú)界承壓含水層的豎向變形量要大于側(cè)向有界承壓含水層,且這種變形量的差異在距離抽水井越遠(yuǎn)的位置越發(fā)明顯。

    (2) 在抽水達(dá)到穩(wěn)定滲流階段時(shí),隨著有界半徑R的增大,抽水引起的側(cè)向有界承壓含水層的豎向變形范圍和具體變形量值均隨之增大,且其豎向變形的分布規(guī)律也越接近于側(cè)向無(wú)界承壓含水層的情況。

    (3) 在整個(gè)抽水階段,側(cè)向無(wú)界和有界越流承壓含水層的徑向變形均隨著抽水時(shí)間的增長(zhǎng)而增大,并在達(dá)到穩(wěn)定滲流時(shí)保持不變且承壓含水層達(dá)到最大徑向位移值;在同一抽水時(shí)間,隨著距抽水井距離的增大,側(cè)向無(wú)界和有界承壓含水層的徑向位移呈現(xiàn)出先增大后減小的變化規(guī)律,即越流承壓含水層的徑向位移的最大值出現(xiàn)在距抽水井一定距離處,且側(cè)向無(wú)界承壓含水層的徑向變形量要小于側(cè)向有界承壓含水層。

    (4) 在抽水達(dá)到穩(wěn)定滲流時(shí),隨著有界半徑R的增大,抽水引起的側(cè)向有界承壓含水層的徑向變形量和變形范圍也隨之增大,其最大徑向變形位置也隨之遠(yuǎn)離抽水井,且抽水引起的側(cè)向無(wú)界承壓含水層的徑向變形范圍要大于側(cè)向有界承壓含水層。

    猜你喜歡
    無(wú)界有界水井
    復(fù)Banach空間的單位球上Bloch-型空間之間的有界的加權(quán)復(fù)合算子
    愛(ài)的大禮物 智能小怪獸 無(wú)界Pro
    車主之友(2022年6期)2023-01-30 07:58:18
    山西發(fā)現(xiàn)一口2000余年前的大型木構(gòu)水井
    朗智無(wú)界 盛享未來(lái)——與朗盛聚合物添加劑業(yè)務(wù)部的深入研討
    上海建材(2020年5期)2020-12-15 00:33:32
    水井的自述
    一類具低階項(xiàng)和退化強(qiáng)制的橢圓方程的有界弱解
    凡水井處皆聽(tīng)單田芳
    淺談?wù)?xiàng)有界周期數(shù)列的一些性質(zhì)
    烏龜與水井
    半無(wú)界區(qū)域上半線性薛定諤方程初邊值問(wèn)題解的破裂及其生命跨度的估計(jì)
    美女内射精品一级片tv| 欧美xxⅹ黑人| 午夜影院在线不卡| 在线观看www视频免费| 一区二区三区四区激情视频| 国产欧美日韩综合在线一区二区| 亚洲一码二码三码区别大吗| 精品人妻在线不人妻| 亚洲五月色婷婷综合| 美女内射精品一级片tv| 男女高潮啪啪啪动态图| 婷婷色综合大香蕉| 寂寞人妻少妇视频99o| 捣出白浆h1v1| 又黄又爽又刺激的免费视频.| 老司机影院成人| 天堂8中文在线网| 日日撸夜夜添| 免费观看av网站的网址| 欧美人与善性xxx| 国产高清三级在线| 热re99久久精品国产66热6| 精品国产国语对白av| 国产精品一区二区在线不卡| 波多野结衣一区麻豆| a级片在线免费高清观看视频| 九色亚洲精品在线播放| 久久韩国三级中文字幕| 男女啪啪激烈高潮av片| 欧美老熟妇乱子伦牲交| 中文字幕精品免费在线观看视频 | 黄色怎么调成土黄色| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品,欧美精品| 香蕉丝袜av| 欧美日本中文国产一区发布| 一级爰片在线观看| 亚洲av综合色区一区| av不卡在线播放| 在线观看美女被高潮喷水网站| 99久久中文字幕三级久久日本| 午夜福利视频在线观看免费| 少妇熟女欧美另类| 又黄又粗又硬又大视频| 亚洲情色 制服丝袜| 成人国产av品久久久| 亚洲成av片中文字幕在线观看 | 卡戴珊不雅视频在线播放| 日韩中文字幕视频在线看片| 五月伊人婷婷丁香| 永久免费av网站大全| 久久久国产精品麻豆| 午夜福利,免费看| 巨乳人妻的诱惑在线观看| 黄色 视频免费看| 9热在线视频观看99| 观看美女的网站| √禁漫天堂资源中文www| 搡女人真爽免费视频火全软件| 日韩一本色道免费dvd| 最近中文字幕2019免费版| 亚洲精品一二三| 久久久a久久爽久久v久久| 国产一级毛片在线| 国产爽快片一区二区三区| 天堂俺去俺来也www色官网| 大香蕉97超碰在线| 亚洲 欧美一区二区三区| 一级a做视频免费观看| 国产一区二区在线观看av| 成人手机av| 天天操日日干夜夜撸| 色94色欧美一区二区| 欧美激情 高清一区二区三区| 久久久久久久亚洲中文字幕| 亚洲国产成人一精品久久久| 久久精品夜色国产| 2021少妇久久久久久久久久久| 亚洲欧美日韩卡通动漫| 久久久精品94久久精品| 一本色道久久久久久精品综合| 五月伊人婷婷丁香| 亚洲精品日韩在线中文字幕| 国产av精品麻豆| 久久久久精品人妻al黑| 欧美变态另类bdsm刘玥| 国产精品麻豆人妻色哟哟久久| 中文字幕免费在线视频6| 麻豆精品久久久久久蜜桃| www日本在线高清视频| 国产成人a∨麻豆精品| 国产日韩一区二区三区精品不卡| 又粗又硬又长又爽又黄的视频| 亚洲av男天堂| 国产一区二区在线观看av| 日本午夜av视频| 免费高清在线观看视频在线观看| 国产成人精品无人区| 九九在线视频观看精品| videosex国产| 久久鲁丝午夜福利片| 国产探花极品一区二区| 国产精品国产av在线观看| 精品国产一区二区三区四区第35| 国产爽快片一区二区三区| 婷婷成人精品国产| 日韩 亚洲 欧美在线| 99久国产av精品国产电影| 咕卡用的链子| 久久韩国三级中文字幕| 久久青草综合色| 免费观看性生交大片5| 黄片播放在线免费| 国产片内射在线| 国产精品欧美亚洲77777| 国产男女内射视频| 中文精品一卡2卡3卡4更新| 七月丁香在线播放| 亚洲国产精品专区欧美| 草草在线视频免费看| 国产国语露脸激情在线看| 欧美人与性动交α欧美精品济南到 | 欧美日韩一区二区视频在线观看视频在线| videosex国产| 看十八女毛片水多多多| 国产欧美亚洲国产| 国产在线视频一区二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲色图 男人天堂 中文字幕 | www.熟女人妻精品国产 | 久久精品熟女亚洲av麻豆精品| 在线观看一区二区三区激情| 国产精品麻豆人妻色哟哟久久| 午夜影院在线不卡| 老女人水多毛片| 男女午夜视频在线观看 | 欧美日韩视频精品一区| 国产 一区精品| 中文欧美无线码| 欧美日韩视频精品一区| 制服诱惑二区| 中文字幕人妻熟女乱码| 免费播放大片免费观看视频在线观看| av福利片在线| 国产极品粉嫩免费观看在线| 搡女人真爽免费视频火全软件| 日本wwww免费看| 啦啦啦在线观看免费高清www| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品人妻久久久影院| 大码成人一级视频| 国产精品免费大片| 亚洲欧洲日产国产| 一级a做视频免费观看| 制服人妻中文乱码| 在线观看国产h片| 一本久久精品| 韩国av在线不卡| 97精品久久久久久久久久精品| 最近2019中文字幕mv第一页| 少妇熟女欧美另类| 妹子高潮喷水视频| 中文字幕制服av| 亚洲精品一区蜜桃| 欧美精品高潮呻吟av久久| 青春草视频在线免费观看| av视频免费观看在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一本久久精品| 国产老妇伦熟女老妇高清| 中文乱码字字幕精品一区二区三区| 少妇 在线观看| 啦啦啦啦在线视频资源| 狠狠精品人妻久久久久久综合| 伊人久久国产一区二区| 亚洲第一av免费看| av视频免费观看在线观看| 亚洲久久久国产精品| 亚洲av电影在线进入| 欧美精品一区二区大全| 欧美人与性动交α欧美精品济南到 | 香蕉精品网在线| 男女午夜视频在线观看 | 国产精品人妻久久久影院| 黑人猛操日本美女一级片| 看免费av毛片| 亚洲国产最新在线播放| 一区二区三区乱码不卡18| 人妻人人澡人人爽人人| 欧美日韩av久久| 99热这里只有是精品在线观看| 天天躁夜夜躁狠狠躁躁| 国产精品麻豆人妻色哟哟久久| 成人午夜精彩视频在线观看| 中国美白少妇内射xxxbb| 国产精品欧美亚洲77777| 久久久久久久大尺度免费视频| 精品熟女少妇av免费看| 22中文网久久字幕| 亚洲精品自拍成人| 欧美人与性动交α欧美精品济南到 | 精品卡一卡二卡四卡免费| 五月天丁香电影| 婷婷成人精品国产| 丁香六月天网| 日韩欧美一区视频在线观看| 深夜精品福利| 国产有黄有色有爽视频| 久久免费观看电影| 亚洲天堂av无毛| 亚洲国产av新网站| www.熟女人妻精品国产 | 另类亚洲欧美激情| 国产欧美亚洲国产| 寂寞人妻少妇视频99o| 国产精品成人在线| 亚洲综合色惰| 男女高潮啪啪啪动态图| 久久国产亚洲av麻豆专区| 曰老女人黄片| 国产精品久久久久久精品古装| 丁香六月天网| 欧美人与善性xxx| 欧美日韩一区二区视频在线观看视频在线| 男女啪啪激烈高潮av片| 美女主播在线视频| 一本色道久久久久久精品综合| 成人18禁高潮啪啪吃奶动态图| 蜜桃国产av成人99| 国产午夜精品一二区理论片| 国产成人精品福利久久| 成年美女黄网站色视频大全免费| 日韩一本色道免费dvd| 26uuu在线亚洲综合色| 国产精品人妻久久久影院| 国产精品一区二区在线不卡| 在线观看美女被高潮喷水网站| 午夜福利视频精品| 99久久精品国产国产毛片| 熟女人妻精品中文字幕| 99热这里只有是精品在线观看| 18禁动态无遮挡网站| 久久精品国产a三级三级三级| 一级爰片在线观看| 日韩成人伦理影院| 欧美激情极品国产一区二区三区 | 亚洲av电影在线观看一区二区三区| 久久久久久久精品精品| 婷婷色综合www| a级毛色黄片| 老熟女久久久| 免费人成在线观看视频色| 亚洲第一av免费看| 亚洲,欧美精品.| 亚洲人成77777在线视频| a级毛片在线看网站| 日韩精品有码人妻一区| 午夜久久久在线观看| 黄色怎么调成土黄色| 亚洲综合精品二区| 97精品久久久久久久久久精品| 中文字幕亚洲精品专区| 日韩视频在线欧美| 伦理电影大哥的女人| 国产免费一级a男人的天堂| 中文字幕人妻熟女乱码| 国产精品免费大片| 亚洲精品av麻豆狂野| 成年人午夜在线观看视频| 久久久久久人人人人人| 男人添女人高潮全过程视频| 少妇高潮的动态图| 天天操日日干夜夜撸| 一级毛片 在线播放| 丰满迷人的少妇在线观看| 18禁在线无遮挡免费观看视频| 国产成人aa在线观看| 亚洲av在线观看美女高潮| 国产成人午夜福利电影在线观看| 欧美最新免费一区二区三区| 大话2 男鬼变身卡| 亚洲欧美成人综合另类久久久| 制服人妻中文乱码| www日本在线高清视频| 欧美日韩亚洲高清精品| 老女人水多毛片| 国产探花极品一区二区| 纵有疾风起免费观看全集完整版| 69精品国产乱码久久久| 九色成人免费人妻av| 精品卡一卡二卡四卡免费| 中文字幕制服av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲欧美中文字幕日韩二区| 精品99又大又爽又粗少妇毛片| 人妻 亚洲 视频| 精品国产露脸久久av麻豆| 国产精品国产三级专区第一集| 最近最新中文字幕免费大全7| 人人妻人人添人人爽欧美一区卜| 夜夜骑夜夜射夜夜干| 欧美日韩亚洲高清精品| 午夜免费男女啪啪视频观看| 久久人妻熟女aⅴ| 人人妻人人澡人人看| av天堂久久9| 中文字幕精品免费在线观看视频 | 国产精品偷伦视频观看了| 亚洲综合色网址| 9191精品国产免费久久| 婷婷色av中文字幕| 日韩一区二区视频免费看| 久久99一区二区三区| 十分钟在线观看高清视频www| 亚洲色图综合在线观看| a级毛色黄片| 成人午夜精彩视频在线观看| 精品熟女少妇av免费看| 69精品国产乱码久久久| 久久久久人妻精品一区果冻| 欧美精品亚洲一区二区| 欧美精品高潮呻吟av久久| 少妇高潮的动态图| 蜜臀久久99精品久久宅男| 国产精品欧美亚洲77777| 春色校园在线视频观看| 纯流量卡能插随身wifi吗| 国产又色又爽无遮挡免| 制服丝袜香蕉在线| 成人国产麻豆网| 少妇猛男粗大的猛烈进出视频| 建设人人有责人人尽责人人享有的| 免费人成在线观看视频色| 亚洲欧美精品自产自拍| 国产精品国产三级国产av玫瑰| 好男人视频免费观看在线| 建设人人有责人人尽责人人享有的| 久久久精品免费免费高清| 人成视频在线观看免费观看| 最近最新中文字幕免费大全7| 卡戴珊不雅视频在线播放| 晚上一个人看的免费电影| 亚洲,欧美精品.| 免费人妻精品一区二区三区视频| 日韩av免费高清视频| 2021少妇久久久久久久久久久| 免费黄频网站在线观看国产| 亚洲av电影在线进入| 欧美bdsm另类| 校园人妻丝袜中文字幕| 青春草国产在线视频| 午夜视频国产福利| 成人毛片a级毛片在线播放| 成人国产av品久久久| 国产免费福利视频在线观看| 成人国语在线视频| 亚洲av免费高清在线观看| tube8黄色片| 中文字幕人妻丝袜制服| 精品久久蜜臀av无| 欧美人与性动交α欧美软件 | 久久久久久伊人网av| 王馨瑶露胸无遮挡在线观看| 亚洲国产精品国产精品| 一级毛片黄色毛片免费观看视频| 日韩大片免费观看网站| 精品久久久精品久久久| 国产一区有黄有色的免费视频| 中文字幕人妻丝袜制服| 国产精品人妻久久久影院| 欧美国产精品va在线观看不卡| 久久精品国产亚洲av涩爱| 日韩一本色道免费dvd| 欧美国产精品一级二级三级| 最近2019中文字幕mv第一页| 亚洲人成网站在线观看播放| 日韩人妻精品一区2区三区| 色哟哟·www| 一二三四在线观看免费中文在 | 人成视频在线观看免费观看| 精品久久国产蜜桃| 少妇 在线观看| 亚洲,欧美精品.| 哪个播放器可以免费观看大片| 精品少妇久久久久久888优播| 香蕉精品网在线| 伦理电影大哥的女人| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色综合www| 一区二区日韩欧美中文字幕 | 2018国产大陆天天弄谢| 亚洲国产精品成人久久小说| 午夜视频国产福利| 久久99精品国语久久久| 国产成人av激情在线播放| 国产爽快片一区二区三区| 免费高清在线观看日韩| av又黄又爽大尺度在线免费看| a 毛片基地| 中国美白少妇内射xxxbb| 国产 一区精品| 美女内射精品一级片tv| 少妇精品久久久久久久| 天天操日日干夜夜撸| 免费av中文字幕在线| 久久99热6这里只有精品| 亚洲国产日韩一区二区| www.色视频.com| 国产老妇伦熟女老妇高清| 免费看不卡的av| 99热网站在线观看| 五月天丁香电影| 久久精品国产亚洲av涩爱| 亚洲av男天堂| 人妻少妇偷人精品九色| 亚洲精品一二三| 欧美日韩精品成人综合77777| 最近手机中文字幕大全| 欧美变态另类bdsm刘玥| 成人黄色视频免费在线看| 在线观看免费高清a一片| 一个人免费看片子| 母亲3免费完整高清在线观看 | 国产成人精品无人区| 久久久久人妻精品一区果冻| 国产欧美日韩一区二区三区在线| 人体艺术视频欧美日本| 国产激情久久老熟女| 在线观看www视频免费| 国产白丝娇喘喷水9色精品| 色吧在线观看| 少妇 在线观看| 成人免费观看视频高清| 黑人高潮一二区| 少妇熟女欧美另类| 亚洲人成77777在线视频| 又粗又硬又长又爽又黄的视频| 青春草国产在线视频| 90打野战视频偷拍视频| 亚洲内射少妇av| 美女视频免费永久观看网站| 久久久久久久久久人人人人人人| 免费av中文字幕在线| 自线自在国产av| 22中文网久久字幕| 丁香六月天网| 一级毛片黄色毛片免费观看视频| 国产视频首页在线观看| 高清在线视频一区二区三区| 母亲3免费完整高清在线观看 | 欧美日韩成人在线一区二区| 国精品久久久久久国模美| 国产精品麻豆人妻色哟哟久久| 亚洲精品久久久久久婷婷小说| 国产成人a∨麻豆精品| 秋霞伦理黄片| 亚洲美女视频黄频| 亚洲成人一二三区av| 国产黄色免费在线视频| a 毛片基地| 1024视频免费在线观看| 老司机影院成人| 午夜福利在线观看免费完整高清在| 精品少妇内射三级| 亚洲精品一二三| 中文精品一卡2卡3卡4更新| 欧美日韩亚洲高清精品| 丝袜人妻中文字幕| 咕卡用的链子| 亚洲成av片中文字幕在线观看 | 大话2 男鬼变身卡| 啦啦啦视频在线资源免费观看| 日韩视频在线欧美| 日韩伦理黄色片| 免费黄色在线免费观看| 国产av精品麻豆| 中国国产av一级| 国产xxxxx性猛交| 免费在线观看黄色视频的| 亚洲国产精品成人久久小说| 日韩熟女老妇一区二区性免费视频| 九色成人免费人妻av| 国产爽快片一区二区三区| 人成视频在线观看免费观看| 日韩成人av中文字幕在线观看| 午夜福利视频在线观看免费| 中文欧美无线码| 啦啦啦啦在线视频资源| 国产黄频视频在线观看| 天堂俺去俺来也www色官网| 亚洲精品乱码久久久久久按摩| 欧美 亚洲 国产 日韩一| 成人免费观看视频高清| 狂野欧美激情性bbbbbb| 国产成人精品福利久久| 中国国产av一级| 久久久精品免费免费高清| 国产一区亚洲一区在线观看| 免费黄网站久久成人精品| 欧美日韩视频精品一区| 国产 精品1| 精品熟女少妇av免费看| 一级毛片我不卡| 女人久久www免费人成看片| 成人漫画全彩无遮挡| 精品人妻在线不人妻| 爱豆传媒免费全集在线观看| 久久精品国产a三级三级三级| 午夜av观看不卡| 国产乱来视频区| 精品人妻偷拍中文字幕| 亚洲综合精品二区| 精品酒店卫生间| 日本爱情动作片www.在线观看| 欧美+日韩+精品| 狂野欧美激情性xxxx在线观看| 午夜免费观看性视频| 亚洲国产欧美在线一区| 男人舔女人的私密视频| 亚洲熟女精品中文字幕| 久久99热6这里只有精品| 一级片免费观看大全| 国语对白做爰xxxⅹ性视频网站| 国产欧美亚洲国产| 精品少妇久久久久久888优播| 国产成人一区二区在线| 国产视频首页在线观看| 两个人免费观看高清视频| 精品午夜福利在线看| 久久久国产欧美日韩av| 亚洲av.av天堂| 亚洲图色成人| 日韩伦理黄色片| 香蕉精品网在线| 免费高清在线观看视频在线观看| 国产 精品1| 一级片免费观看大全| 97人妻天天添夜夜摸| 久久99一区二区三区| 久久国内精品自在自线图片| 人体艺术视频欧美日本| 好男人视频免费观看在线| 成人亚洲精品一区在线观看| 久久综合国产亚洲精品| 人人妻人人爽人人添夜夜欢视频| 制服人妻中文乱码| 日韩在线高清观看一区二区三区| 欧美少妇被猛烈插入视频| 亚洲精华国产精华液的使用体验| 精品久久蜜臀av无| 亚洲欧美色中文字幕在线| 亚洲一级一片aⅴ在线观看| 女人被躁到高潮嗷嗷叫费观| 日本色播在线视频| 亚洲av综合色区一区| 免费人妻精品一区二区三区视频| 久久人人爽人人片av| 日产精品乱码卡一卡2卡三| 精品一区二区三区视频在线| 亚洲欧美一区二区三区黑人 | 一区二区三区四区激情视频| 国产探花极品一区二区| 一级片'在线观看视频| 国产又爽黄色视频| 精品午夜福利在线看| 久久久久久久亚洲中文字幕| 十分钟在线观看高清视频www| 国产亚洲av片在线观看秒播厂| 少妇人妻 视频| 国产一级毛片在线| av福利片在线| 欧美成人精品欧美一级黄| 一区在线观看完整版| 亚洲国产成人一精品久久久| 亚洲美女黄色视频免费看| 久久久精品区二区三区| 免费人妻精品一区二区三区视频| 免费黄色在线免费观看| 亚洲精品乱码久久久久久按摩| 侵犯人妻中文字幕一二三四区| av网站免费在线观看视频| 精品少妇内射三级| 国产欧美亚洲国产| 亚洲精品av麻豆狂野| 国产成人精品在线电影| 热re99久久国产66热| 亚洲经典国产精华液单| 亚洲成人手机| 亚洲人成77777在线视频| 好男人视频免费观看在线| 久久国产精品大桥未久av| 国产精品久久久久久av不卡| 国产免费现黄频在线看| 国产1区2区3区精品| 久久久久久久久久久免费av| 在线观看免费高清a一片| 久久久国产精品麻豆| 亚洲婷婷狠狠爱综合网| 久久毛片免费看一区二区三区| 国产国拍精品亚洲av在线观看| 精品99又大又爽又粗少妇毛片| 免费高清在线观看视频在线观看| 久久久精品区二区三区| 母亲3免费完整高清在线观看 | 国产欧美日韩综合在线一区二区| 久久久久久久国产电影| 日韩不卡一区二区三区视频在线| 欧美成人午夜免费资源| 成人黄色视频免费在线看| 亚洲国产成人一精品久久久| 日韩av不卡免费在线播放| 亚洲欧美一区二区三区国产| 国产黄色免费在线视频| 亚洲精品成人av观看孕妇| 成人手机av|