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

    基于DIMINE-MIDAS/GTS的采空區(qū)穩(wěn)定性分析

    2020-10-16 01:19:10任青云
    有色金屬(礦山部分) 2020年5期
    關(guān)鍵詞:礦柱中段力學(xué)

    任青云

    (湖南有色金屬研究院,長沙 410100)

    地下礦山開采的重點安全問題是開采形成的采空區(qū)對后續(xù)開采產(chǎn)生的巖層移動問題,開采過程中伴隨的地層沉降及地壓活動[1-2]。礦山開采面臨的重要安全問題主要表現(xiàn)在,開采區(qū)域的不斷擴大,地層隨之進(jìn)行持續(xù)的運動,由于開采現(xiàn)場會出現(xiàn)不同的地質(zhì)條件,例如溶洞、斷層、空區(qū)等可能會產(chǎn)生移動變形的大范圍空間,此時會造成地下礦山的地壓顯現(xiàn),導(dǎo)致頂板出現(xiàn)大變形,繼而影響圍巖發(fā)生冒落和片幫等現(xiàn)象。因此地下礦山長期開采導(dǎo)致的采空區(qū)穩(wěn)定性是影響礦山生命財產(chǎn)安全及后續(xù)開采的重要因素[3-4]。

    影響采空區(qū)的穩(wěn)定性及礦山地表位移沉降的主要原因有開采礦體的埋深、采場的高度,礦體的厚度、傾角、巖石力學(xué)特性及礦山采用的采礦工藝。采用大型有限元數(shù)值模擬[5-6]前處理模式可以直觀地表示出地下采空區(qū)的分布情況及整改礦山的輪廓,其后處理模式能夠計算出礦山的地表位移沉降、采空區(qū)是否發(fā)生拉應(yīng)力破壞、礦柱是否產(chǎn)生主應(yīng)力集中及塑性區(qū)變形破壞情況。本文根據(jù)南川河石灰?guī)r礦礦區(qū)地質(zhì)資料,結(jié)合礦體情況及開采設(shè)計,對南川河石灰?guī)r礦進(jìn)行取樣試驗,采用數(shù)值模擬方法分析礦山開采后巖層移動變形情況,并結(jié)合傳統(tǒng)的理論計算方法對礦柱、頂板穩(wěn)定性進(jìn)行分析,以此確定無支護(hù)狀態(tài)下采空區(qū)跨度,從而評判采空區(qū)的穩(wěn)定性及地表位移沉降的影響范圍和程度。

    1 工程概況

    南川河石灰?guī)r礦位于瀏陽市澄潭江鎮(zhèn)南川河畔,距瀏陽市城區(qū)直距 21 km,公司始建于1972年,原采用露天開采,1996年開始地下開采,主要采用斜坡道開拓,對角式通風(fēng)。礦山長期使用的采礦方法是房柱采礦法。礦山總的回采順序是沿礦體傾斜由上而下循階段進(jìn)行,按邊探邊采的原則進(jìn)行開采。采場回采順序是后退式開采,礦房內(nèi)采用分層開采。礦山2014年之前開采130 m中段,2014年至今開采100 m中段,目前礦山正在+100 m中段北部進(jìn)行生產(chǎn)。

    由于石灰?guī)r的巖性較軟,并且礦山開采歷史悠久,地下空區(qū)存在極大的不穩(wěn)定性及不確定性,為了掌握空區(qū)對下部礦體回采的影響,必須針對采空區(qū)進(jìn)行穩(wěn)定性研究和分析。根據(jù)礦山的實際生產(chǎn)現(xiàn)狀,針對空區(qū)進(jìn)行了現(xiàn)場調(diào)查,并通過統(tǒng)計分析獲得礦山不同開采水平的采空區(qū)規(guī)模、尺寸。此次現(xiàn)場調(diào)查針對+130 m中段和+100 m中段共統(tǒng)計出130個采空區(qū),+100 m中段采空區(qū)部分調(diào)查數(shù)據(jù)見表1?,F(xiàn)場調(diào)查發(fā)現(xiàn),礦山部分采空區(qū)已形成采空區(qū)群;大部分空區(qū)寬度約12~14 m,采空區(qū)未發(fā)生冒頂片幫事故,側(cè)壁處理較好,有極少量的浮石[7]。溶洞構(gòu)造溶蝕形成,并有軟塑性、硬塑溶洞泥質(zhì)充填,與地表水源連通的構(gòu)造少。

    表1 +100 m中段采空區(qū)特征統(tǒng)計表(部分)

    2 礦巖力學(xué)參數(shù)室內(nèi)試驗

    2.1 礦巖力學(xué)參數(shù)試驗

    為了合理評價南川河石灰?guī)r礦工程地質(zhì)條件,為后續(xù)數(shù)值模擬提供依據(jù),開展了礦巖力學(xué)參數(shù)室內(nèi)試驗,測試的內(nèi)容包括礦巖的容重、單軸抗壓強度、抗拉強度、彈性模量、泊松比等基本物理力學(xué)參數(shù)。針對礦山實際情況,選取了石灰?guī)r巖樣進(jìn)行加工,并對各巖樣進(jìn)行了相關(guān)參數(shù)測定。巖樣總計22個,其中5個用于開展單軸壓縮試驗、5個用于巴西劈裂試驗、12個用于變角剪試驗。

    根據(jù)數(shù)值模擬計算分析需要的基本數(shù)據(jù),主要進(jìn)行了如下四項內(nèi)容的實驗室試驗:1)礦、巖的天然容重;2)礦、巖的單軸抗壓力學(xué)特性;3)巖體的劈裂抗拉強度試驗;4)巖體的泊松比與彈性模量試驗。

    試驗開始前,先對巖樣加工。首先,按照標(biāo)準(zhǔn)巖石力學(xué)試驗所要求的規(guī)格,制成試樣尺寸為Φ50 mm×100 mm的標(biāo)準(zhǔn)巖石試樣;然后擦干試樣,在自然狀態(tài)下風(fēng)干;最后,試驗在250 t全數(shù)字型液壓伺服剛性巖石力學(xué)試驗系統(tǒng)(RMT150,圖1)及30 t萬能材料試驗機上完成。

    圖1 RMT150型試驗設(shè)備Fig.1 RMT150 test equipment

    2.1.1 單軸抗壓試驗

    試樣尺寸為直徑Ф50 mm,高徑比約為2∶1,巖石單軸抗壓試驗結(jié)果如表2所示。

    表2 巖石單軸抗壓強度及靜力受壓彈性模量和泊松比測試

    2.1.2 抗拉強度試驗

    巖石的抗拉強度試驗如圖2所示,各種試驗結(jié)果見表3和表4。

    圖2 巴西劈裂試驗Fig.2 Brazilian splitting experiment

    表3 各試樣劈裂抗拉試驗參數(shù)測定值

    表4 變角剪試驗結(jié)果

    根據(jù)單軸抗壓試驗和單軸抗拉試驗,將室內(nèi)試驗結(jié)果匯總?cè)绫?所示。

    表5 礦巖力學(xué)參數(shù)試驗結(jié)果匯總表

    2.2 巖體強度參數(shù)分析

    由于自然巖體內(nèi)部有很多缺陷,比如節(jié)理、裂隙和孔洞等,因此自然巖體與小試樣巖石的力學(xué)性質(zhì)是有差異的,導(dǎo)致實驗室測得小試樣巖石力學(xué)參數(shù)不能直接作為數(shù)值模擬的力學(xué)參數(shù)。因此,巖體力學(xué)工程界,大多數(shù)依據(jù)室內(nèi)巖石力學(xué)實驗,采用理論力學(xué)的方法,估算出巖體的力學(xué)參數(shù),用于有限元數(shù)值模擬計算。

    根據(jù)HOEK,CARRANZA-TORRES提出的巖體破壞經(jīng)驗準(zhǔn)則:

    (1)

    巖體的mb、S和a由GSI和上式計算的mi確定。

    (2)

    (3)

    (4)

    根據(jù)確定的mb、S和a值,可以計算出巖體的單軸抗壓強度和抗拉強度:

    σc=σci·sa

    (5)

    (6)

    巖體的內(nèi)摩擦角φ和黏聚力c可以由以下公式進(jìn)行計算:

    (7)

    (8)

    針對南川河石灰?guī)r礦工程地質(zhì)條件,取石灰?guī)rGSI=70,根據(jù)室內(nèi)巖石力學(xué)參數(shù)試驗結(jié)果,利用HOEK和CARRANZA-TORRES提出的巖體破壞經(jīng)驗準(zhǔn)則及巖體分類RMR值,采用Roclab軟件進(jìn)行計算。選取巖體的力學(xué)參數(shù)如表6所示。

    表6 巖體物理力學(xué)參數(shù)匯總表

    3 采空區(qū)穩(wěn)定性數(shù)值分析

    3.1 三維礦體模型的建立

    制約礦山地表位移沉降及采空區(qū)穩(wěn)定性的重要因素,包括礦區(qū)地表的地形地質(zhì)條件及井下采空區(qū)的分布狀態(tài),因此在三維有限元數(shù)值模擬計算過程中,缺少完整的地形模塊,不能真實地反映礦山的實際情況。本研究通過基于DIMINE強度的前處理功能,利用二維地質(zhì)剖面圖或者M(jìn)APGIS地形數(shù)據(jù)庫生成三維地質(zhì)模型,對勘探線剖面進(jìn)行初步的處理,僅保留勘探線與礦體解譯線[8]。將礦山剖面圖經(jīng)過CAD處理之后導(dǎo)入到DIMINE中,進(jìn)行線串運算使勘探線“立起來”,對立起來的三維線文件進(jìn)行坐標(biāo)運算,使各點坐標(biāo)為其真實坐標(biāo)值。通過上述基本工作的完成,可以聯(lián)接各礦體解譯線剖面,對于比較復(fù)雜的礦體解譯線,可以構(gòu)造中線,通過中線與剖面線來生成實體,礦體實體模型見圖3。在DIMINE中生成的礦山實體模型,無法直接用于有限元數(shù)值模擬的運算,需要將此模型導(dǎo)入到MIDAS/GTS中,利用MIDAS/GTS的網(wǎng)格劃分功能,結(jié)合采礦工藝,進(jìn)一步建立采區(qū)模型,劃分網(wǎng)格單元,進(jìn)行數(shù)值模擬分析。

    圖3 南川河石灰?guī)r礦三維實體模型Fig.3 3D solid model of Nanchuanhe limestone mine

    3.2 理論基礎(chǔ)及相關(guān)參數(shù)設(shè)定

    根據(jù)礦山的地層參數(shù),確定有限元數(shù)值模擬的巖體力學(xué)參數(shù)見表6。計算過程礦巖及廢石均采用摩爾-庫倫((Mohr-Coulomb)屈服準(zhǔn)則,該屈服準(zhǔn)則的控制方程為:

    (9)

    最大拉應(yīng)力屈服準(zhǔn)則函數(shù)為:

    ft=σ3-σt

    (10)

    式中:c—黏聚力,MPa;φ—內(nèi)摩擦角,(°);σ1—最大主應(yīng)力,MPa;σ3—最小主應(yīng)力,MPa;σt—抗拉強度,MPa。MIDAS/GTS有限元計算程序,默認(rèn)當(dāng)fs大于0時,表明巖體已發(fā)生了剪切破壞;當(dāng)ft大于0時,表明巖體已發(fā)生了拉伸破壞。

    為了模型建立與分析結(jié)果的準(zhǔn)確性,本數(shù)值模擬進(jìn)行了以下假設(shè):模型的礦巖體為均質(zhì)、各向同性材料,模擬的三維網(wǎng)格模型底部為固定約束,四周為限制水平約束,只在豎直方向有位移變形。初始階段考慮地應(yīng)力與自重的影響,位移清零,表明礦山在原始狀態(tài)已處于穩(wěn)定的固結(jié)沉降。計算模擬不考慮地震、爆破等動荷載的影響,不分析地下水的影響;礦山地質(zhì)條件較復(fù)雜,礦巖穩(wěn)固性在空間分布上具有較大的隨機性,礦體為斷層分割,且礦巖和圍巖中均有穿插,這些地質(zhì)現(xiàn)象的影響在巖體參數(shù)折減時已做考慮,在模擬過程中不予另行考慮。

    根據(jù)地下采空區(qū)對礦山圍巖的影響范圍,設(shè)置有限元數(shù)模網(wǎng)格模型的邊界固定條件如下:1)模型X方向的邊界為X方向固定邊界,取u=0,v≠0,w≠0(u為X軸方向的位移,v為Y軸方向的位移,w為Z軸方向的位移);2)模型Y方向邊界為Y方向固定邊界,取u≠0,v=0,w≠0;3)模型底部為全固定邊界,取u=0,v=0,w=0;4)模型地表部分為自由面,不限制位移約束。礦體邊界約束情況及三維網(wǎng)絡(luò)質(zhì)量情況見圖4和圖5。

    圖4 三維網(wǎng)格模型+邊界條件Fig.4 3D grid model + boundary conditions

    圖5 網(wǎng)格質(zhì)量檢查Fig.5 Grid quality inspection

    3.3 數(shù)值模擬計算結(jié)果分析

    3.3.1 位移場分析

    采空區(qū)的穩(wěn)定性分析,位移是評判礦山穩(wěn)定性的一個重要指標(biāo),因為大量的地表變形,會使礦山發(fā)生地表塌陷,且位移變形量超過巖體的變形限值,會發(fā)生采空區(qū)頂板冒頂片幫事故。從圖6~9的位移云圖得出,該礦山產(chǎn)生的位移變化情況,在礦體的頂板位移值小于0,說明由于產(chǎn)生了采空區(qū),發(fā)生了頂板下沉變形情況,但隨著礦體的開挖過程進(jìn)行,在采空區(qū)的底部及巷道的底板,豎直方向位移大于0,是因為礦體開挖深度加大,采空區(qū)及巷道周圍發(fā)生了主應(yīng)力集中現(xiàn)象,水平方向的應(yīng)力大于垂直方向的應(yīng)力,使巷道底板產(chǎn)生上拱的位移變形。

    從圖6豎向位移云圖可以看出,+130 m中段礦體開采時,采空區(qū)頂板豎向位移最大值為-1.74 mm,礦體開采底板巷道發(fā)生向上最大位移為+2.39 mm。

    圖6 +130 m中段頂板位移云圖Fig.6 Displacement cloud image of the +130 m middle section roof

    圖7為礦柱的位移云圖,+130 m中段采空區(qū)礦柱位移變形為受壓沉降,最大位移值為-1.8 mm,+100 m中段礦柱由于+130 m中段礦體已全部回采,+100 m礦柱所受側(cè)向壓力大于豎向壓力,變形為部分礦柱產(chǎn)生上拱位移,最大位移值為+2.5 mm。

    圖7 礦柱位移云圖Fig.7 Cloud diagram of pillar displacement

    圖8為+100 m中段礦體開采后采空區(qū)頂板產(chǎn)生的位移沉降變形云圖,目前在+100 m中段回采時采空區(qū)頂板最大位移沉降為-2.4 mm,采空區(qū)底板巷道上拱位移值為+2.94 mm。

    圖8 +100 m中段空區(qū)頂板位移云圖Fig.8 Displacement cloud image of the +100 m middle section roof

    采空區(qū)形成的地表沉降有微弱的影響,地表沉降值不足1 mm,該位移沉降整體上可以忽略,采空區(qū)地表沉降等值線范圍見圖9,位移變形值數(shù)據(jù)見表7。

    圖9 地表沉降位移云圖與沉降等值線圖Fig.9 Surface displacement cloud map and subsidence contour map

    表7 采空區(qū)位移變化值計算結(jié)果

    3.3.2 應(yīng)力場分析

    礦體開采后對圍巖造成的影響是一個非線性的不可逆的加載過程,它使得處于初始地應(yīng)力狀態(tài)下的圍巖進(jìn)行應(yīng)力重新分布,最后達(dá)到新的平衡,礦體開挖引起徑向應(yīng)力釋放、切向應(yīng)力增加,導(dǎo)致礦體開采頂板和底板產(chǎn)生拉應(yīng)力;從應(yīng)力云圖可以看出礦體開采區(qū)域與礦柱接觸處發(fā)生較大的應(yīng)力集中,周邊最大。采空區(qū)影響范圍一定深度后,逐步變化為恢復(fù)到原巖應(yīng)力狀態(tài)。采空區(qū)形成后的主應(yīng)力值都小于0,表明主應(yīng)力(壓應(yīng)力)為受壓狀態(tài)。

    圖10 最小主應(yīng)力云圖Fig.10 The minimum principal stress cloud diagram

    圖11 最大主應(yīng)力云圖Fig.11 The maximum principal stress cloud diagram

    根據(jù)計算結(jié)果圖10、圖11得出結(jié)論:在采空區(qū)的頂板和底板產(chǎn)生拉應(yīng)力破壞,其中+130 m中段采空區(qū)頂板產(chǎn)生的最大拉應(yīng)力約為2.27 MPa;+100 m中段采空區(qū)頂板產(chǎn)生的最大拉應(yīng)力為1.85 MPa,已接近巖體的抗拉強度,表明在采空區(qū)的頂板部分已發(fā)生拉應(yīng)力破壞,容易產(chǎn)生冒頂片幫事故。從最大應(yīng)力云圖可以得出:在礦柱與空區(qū)的邊界圍巖出現(xiàn)壓應(yīng)力集中現(xiàn)象,最大壓應(yīng)力集中發(fā)生在礦山+100 m中段東北方向的采空區(qū)與13號礦柱交界處,產(chǎn)生的最大壓應(yīng)力值為12.66 MPa,+130 m中段的(老采空區(qū))礦柱發(fā)生最大壓應(yīng)力集中值為8.88 MPa,表明礦柱發(fā)生受壓破壞。應(yīng)力結(jié)果見表8。

    表8 最大、最小主應(yīng)力值計算結(jié)果

    3.3.3 塑性區(qū)分析

    圖12為分析得出的采空區(qū)塑性變形情況,從計算結(jié)果可以得出,在采空區(qū)的附近發(fā)生了塑性區(qū)變形破壞,塑性區(qū)面積占總面積約44.5%,最大塑性變形為2.06×10-3,礦體開采采空區(qū)周邊礦柱僅有少量單元發(fā)生塑性破壞,塑性區(qū)變化主要出現(xiàn)在采空區(qū)頂板圍巖。

    圖12 塑性區(qū)云圖Fig.12 Cloud image of plastic zone

    4 結(jié)論

    利用DIMINE強大的前處理功能,建立了南川河石灰?guī)r礦的三維地質(zhì)模型,導(dǎo)入到MIDAS/GTS,形成了包含采空區(qū)、礦體、地表地形的三維有限元數(shù)值模型。利用基于DIMINE-MIDAS/GTS的耦合技術(shù),分析了礦體開采后礦山的位移變形、應(yīng)力集中現(xiàn)象及塑性區(qū)破壞情況,分析論證該礦山采空區(qū)的穩(wěn)定性,結(jié)果表明:

    1)+100 m中段采空區(qū)頂板產(chǎn)生的最大位移沉降為-2.4 mm,采空區(qū)的底板與巷道底板發(fā)生最大正向位移+2.94 mm。采空區(qū)形成的地表沉降有微弱的影響,沉降值不足1 mm,表明該礦山開采對地表影響較小。

    2)采空區(qū)頂板發(fā)生了拉應(yīng)力破壞。+130 m中段老采空區(qū)最大拉應(yīng)力為2.27 MPa,接近圍巖的抗拉強度值,采空區(qū)容易發(fā)生冒頂片幫事故。+100 m中段采空區(qū)與13號礦柱接觸處產(chǎn)生了最大壓應(yīng)力值集中12.66 MPa,礦柱處于壓應(yīng)力破壞。

    3)采空區(qū)開采后產(chǎn)生的塑性區(qū)變形主要分布在采空區(qū)的四周邊界,產(chǎn)生的塑性區(qū)變形范圍占總面積的44.5%,其中最大塑性區(qū)變形量為2.06×10-3,采空區(qū)周邊礦柱僅有少量單元發(fā)生塑性破壞,塑性區(qū)主要存在于采空區(qū)頂板。

    4)綜合考慮巖體位移、應(yīng)力、塑性形變等計算結(jié)果,礦體開采采空區(qū)頂板產(chǎn)生少量的豎直位移,地表發(fā)生微量的沉降,部分采空區(qū)頂板發(fā)生塑性變形,只在+130 m中段老采空區(qū)發(fā)生拉應(yīng)力破壞,表明礦山采空區(qū)基本處于穩(wěn)定狀態(tài)。

    猜你喜歡
    礦柱中段力學(xué)
    力學(xué)
    尿檢時如何取中段尿
    某銅鈷礦柱浸試驗研究
    云南化工(2021年9期)2021-12-21 07:43:40
    弟子規(guī)·余力學(xué)文(十)
    快樂語文(2021年11期)2021-07-20 07:41:32
    弟子規(guī)·余力學(xué)文(四)
    快樂語文(2020年30期)2021-01-14 01:05:28
    東天山中段晚古生代剪切帶疊加特征及構(gòu)造控礦作用
    邁好從低段到中段的“坎”
    考試周刊(2018年15期)2018-01-21 17:24:26
    基于動態(tài)強度折減法礦柱漸進(jìn)失穩(wěn)演變規(guī)律
    力學(xué) 等
    鎖骨中段骨折的處理
    天堂√8在线中文| 国产一区二区激情短视频| 色综合婷婷激情| 大码成人一级视频| av网站免费在线观看视频| 久久久久久大精品| 在线观看一区二区三区激情| 精品电影一区二区在线| 热re99久久精品国产66热6| 在线国产一区二区在线| 国产高清国产精品国产三级| 老司机午夜十八禁免费视频| 香蕉国产在线看| av网站免费在线观看视频| 视频区欧美日本亚洲| 在线天堂中文资源库| 99热只有精品国产| 亚洲精品一二三| 亚洲精品久久成人aⅴ小说| 国产av在哪里看| 12—13女人毛片做爰片一| 精品国产亚洲在线| 老司机亚洲免费影院| netflix在线观看网站| 亚洲av成人一区二区三| 91老司机精品| 亚洲欧洲精品一区二区精品久久久| 女人被狂操c到高潮| 欧美日韩福利视频一区二区| 1024视频免费在线观看| 激情视频va一区二区三区| 一边摸一边做爽爽视频免费| 午夜两性在线视频| 少妇 在线观看| 国产精品一区二区精品视频观看| 女性被躁到高潮视频| 黑人操中国人逼视频| 夜夜躁狠狠躁天天躁| 我的亚洲天堂| 人人澡人人妻人| 99精品在免费线老司机午夜| 91老司机精品| 亚洲情色 制服丝袜| 国产成人av激情在线播放| 50天的宝宝边吃奶边哭怎么回事| 一级片'在线观看视频| 最新在线观看一区二区三区| 女人精品久久久久毛片| 桃色一区二区三区在线观看| 色婷婷av一区二区三区视频| 日本vs欧美在线观看视频| 性少妇av在线| 别揉我奶头~嗯~啊~动态视频| 久久久久国内视频| 99在线视频只有这里精品首页| 久久久久国内视频| 亚洲av熟女| 免费看十八禁软件| 日韩大码丰满熟妇| 这个男人来自地球电影免费观看| 久久久久久久久久久久大奶| 一级片免费观看大全| 国产精品成人在线| 国产色视频综合| 久久午夜亚洲精品久久| 成人手机av| 后天国语完整版免费观看| 久久性视频一级片| 欧美精品一区二区免费开放| 在线观看免费午夜福利视频| 国产成人一区二区三区免费视频网站| 国产av又大| 中文字幕人妻熟女乱码| 热99re8久久精品国产| 国产成人欧美在线观看| 久久草成人影院| 正在播放国产对白刺激| 久久影院123| 精品午夜福利视频在线观看一区| 97超级碰碰碰精品色视频在线观看| 老司机午夜十八禁免费视频| 桃色一区二区三区在线观看| 桃色一区二区三区在线观看| 人妻久久中文字幕网| 色尼玛亚洲综合影院| 一级毛片精品| 9191精品国产免费久久| 日日摸夜夜添夜夜添小说| 窝窝影院91人妻| 国产成人精品久久二区二区免费| 日日摸夜夜添夜夜添小说| 9191精品国产免费久久| 高清欧美精品videossex| 激情在线观看视频在线高清| 国产乱人伦免费视频| 中文亚洲av片在线观看爽| 在线观看免费视频日本深夜| 国产亚洲精品久久久久久毛片| 又紧又爽又黄一区二区| 免费在线观看影片大全网站| 午夜精品国产一区二区电影| 老熟妇仑乱视频hdxx| 久久人人97超碰香蕉20202| 日本免费一区二区三区高清不卡 | 伦理电影免费视频| 国产aⅴ精品一区二区三区波| 两个人看的免费小视频| 国产成+人综合+亚洲专区| 中文亚洲av片在线观看爽| 午夜两性在线视频| 成人av一区二区三区在线看| 悠悠久久av| 黑人巨大精品欧美一区二区mp4| 亚洲五月色婷婷综合| 超碰成人久久| 日韩精品中文字幕看吧| 亚洲一区中文字幕在线| 在线观看免费视频日本深夜| 国产无遮挡羞羞视频在线观看| 一进一出抽搐gif免费好疼 | 免费人成视频x8x8入口观看| 国产伦人伦偷精品视频| 国产精品 国内视频| 久久久久九九精品影院| 午夜精品在线福利| 99精品久久久久人妻精品| 在线观看免费午夜福利视频| 久久久国产成人精品二区 | 丝袜在线中文字幕| 亚洲三区欧美一区| 亚洲精品国产区一区二| 乱人伦中国视频| 国产精品影院久久| 国产在线观看jvid| 一区二区日韩欧美中文字幕| 一a级毛片在线观看| 亚洲熟妇中文字幕五十中出 | 国产精品 欧美亚洲| 午夜成年电影在线免费观看| 黑人猛操日本美女一级片| 久久影院123| xxx96com| 欧美乱码精品一区二区三区| 欧美日韩瑟瑟在线播放| 在线观看66精品国产| 精品一区二区三区视频在线观看免费 | 免费不卡黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 成人三级做爰电影| 麻豆一二三区av精品| 亚洲午夜精品一区,二区,三区| 999久久久国产精品视频| 成人三级黄色视频| 美女高潮到喷水免费观看| av在线天堂中文字幕 | 成人手机av| 久久久久久亚洲精品国产蜜桃av| 亚洲avbb在线观看| 麻豆一二三区av精品| 亚洲五月色婷婷综合| 自线自在国产av| 在线观看一区二区三区| 18禁国产床啪视频网站| 国产精品日韩av在线免费观看 | 亚洲中文av在线| 亚洲七黄色美女视频| 日本wwww免费看| 国产精品一区二区免费欧美| 一边摸一边抽搐一进一出视频| 亚洲精品中文字幕在线视频| 国产精品香港三级国产av潘金莲| 国产亚洲欧美在线一区二区| 精品福利永久在线观看| 别揉我奶头~嗯~啊~动态视频| 久久久国产一区二区| 色播在线永久视频| 午夜两性在线视频| 波多野结衣高清无吗| 99热只有精品国产| 最近最新免费中文字幕在线| 欧美日韩视频精品一区| 在线av久久热| 久久久久久久久免费视频了| 国产99久久九九免费精品| 九色亚洲精品在线播放| 丝袜在线中文字幕| 日韩成人在线观看一区二区三区| 亚洲激情在线av| 亚洲aⅴ乱码一区二区在线播放 | 最新在线观看一区二区三区| 淫妇啪啪啪对白视频| 亚洲精品一卡2卡三卡4卡5卡| 日本免费一区二区三区高清不卡 | 国产主播在线观看一区二区| 国产欧美日韩精品亚洲av| 色精品久久人妻99蜜桃| 国产一区二区在线av高清观看| 久久精品91蜜桃| 免费在线观看亚洲国产| 91老司机精品| av在线播放免费不卡| 亚洲自偷自拍图片 自拍| 亚洲精品av麻豆狂野| 老熟妇乱子伦视频在线观看| 999久久久国产精品视频| 日韩欧美一区二区三区在线观看| 成在线人永久免费视频| 69精品国产乱码久久久| 亚洲成国产人片在线观看| 丁香六月欧美| 麻豆成人av在线观看| 亚洲国产精品sss在线观看 | 午夜免费激情av| 国产麻豆69| 黄片播放在线免费| 色播在线永久视频| 成人特级黄色片久久久久久久| 波多野结衣av一区二区av| 欧美黑人欧美精品刺激| 日日干狠狠操夜夜爽| 久99久视频精品免费| 制服人妻中文乱码| 色婷婷久久久亚洲欧美| 涩涩av久久男人的天堂| a级片在线免费高清观看视频| 国产有黄有色有爽视频| 18禁黄网站禁片午夜丰满| 美女扒开内裤让男人捅视频| 怎么达到女性高潮| 欧美黄色片欧美黄色片| 琪琪午夜伦伦电影理论片6080| 精品福利永久在线观看| e午夜精品久久久久久久| 亚洲欧美日韩高清在线视频| 亚洲av五月六月丁香网| 日韩欧美国产一区二区入口| 人人妻人人爽人人添夜夜欢视频| 69av精品久久久久久| 免费不卡黄色视频| 欧美日韩国产mv在线观看视频| 老司机靠b影院| 男女做爰动态图高潮gif福利片 | 天天躁夜夜躁狠狠躁躁| 亚洲av成人av| www.www免费av| 久久伊人香网站| 日韩国内少妇激情av| 日韩精品免费视频一区二区三区| 国产精品av久久久久免费| 99香蕉大伊视频| 久久影院123| 欧美不卡视频在线免费观看 | 久热这里只有精品99| 叶爱在线成人免费视频播放| 亚洲第一欧美日韩一区二区三区| 91九色精品人成在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 免费搜索国产男女视频| 欧美日韩亚洲国产一区二区在线观看| 日日摸夜夜添夜夜添小说| 女同久久另类99精品国产91| 亚洲成人免费av在线播放| 国产一区二区激情短视频| 久久九九热精品免费| 日韩欧美在线二视频| 欧美日韩精品网址| 国产成人系列免费观看| 香蕉丝袜av| e午夜精品久久久久久久| netflix在线观看网站| 悠悠久久av| 国产精品二区激情视频| 成人黄色视频免费在线看| 亚洲色图综合在线观看| 亚洲精品国产区一区二| 亚洲精品一二三| 免费看十八禁软件| 成人18禁高潮啪啪吃奶动态图| 欧美国产精品va在线观看不卡| 波多野结衣av一区二区av| 久久中文字幕一级| 免费不卡黄色视频| 亚洲午夜理论影院| 18禁观看日本| 亚洲中文av在线| 午夜精品在线福利| 国产成人av激情在线播放| 欧美成人午夜精品| 天堂动漫精品| 亚洲黑人精品在线| 国产精品亚洲av一区麻豆| 9191精品国产免费久久| 欧美日本亚洲视频在线播放| 免费观看人在逋| 亚洲国产中文字幕在线视频| 久热爱精品视频在线9| 日韩欧美三级三区| 欧美乱妇无乱码| 高潮久久久久久久久久久不卡| 亚洲精品成人av观看孕妇| 少妇被粗大的猛进出69影院| 狠狠狠狠99中文字幕| 欧美一级毛片孕妇| 色综合婷婷激情| 国产野战对白在线观看| av在线天堂中文字幕 | www.熟女人妻精品国产| av免费在线观看网站| 日韩欧美一区二区三区在线观看| 中亚洲国语对白在线视频| 黑人猛操日本美女一级片| 丰满饥渴人妻一区二区三| 国产黄a三级三级三级人| 欧美日韩一级在线毛片| 中文字幕av电影在线播放| 黑人巨大精品欧美一区二区蜜桃| 色综合站精品国产| 99久久国产精品久久久| 中文字幕高清在线视频| 午夜免费成人在线视频| 丝袜美腿诱惑在线| 精品午夜福利视频在线观看一区| 国产伦一二天堂av在线观看| 日韩 欧美 亚洲 中文字幕| 国产黄色免费在线视频| 欧美久久黑人一区二区| 国产高清视频在线播放一区| 午夜福利免费观看在线| 成人国语在线视频| 好看av亚洲va欧美ⅴa在| 亚洲精品久久成人aⅴ小说| 超色免费av| 人妻丰满熟妇av一区二区三区| 午夜精品国产一区二区电影| 国产成人精品在线电影| 亚洲免费av在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美不卡视频在线免费观看 | 久久精品影院6| 久久这里只有精品19| 国产熟女xx| 久热这里只有精品99| 男人操女人黄网站| 成人精品一区二区免费| 日韩成人在线观看一区二区三区| 国产伦人伦偷精品视频| 黄色丝袜av网址大全| 一进一出抽搐gif免费好疼 | 中文字幕高清在线视频| 欧美日本亚洲视频在线播放| 在线av久久热| 久久国产精品影院| 久久精品影院6| 色婷婷av一区二区三区视频| 国产精品久久久久成人av| 精品人妻在线不人妻| 国产精华一区二区三区| 美女福利国产在线| 亚洲国产毛片av蜜桃av| 交换朋友夫妻互换小说| 免费av中文字幕在线| 久久婷婷成人综合色麻豆| 国产免费现黄频在线看| 成年女人毛片免费观看观看9| 欧美精品亚洲一区二区| www.999成人在线观看| 很黄的视频免费| 亚洲午夜精品一区,二区,三区| 亚洲男人天堂网一区| 精品国产一区二区久久| 欧美激情久久久久久爽电影 | 亚洲第一青青草原| 久99久视频精品免费| 国产免费男女视频| 精品国产美女av久久久久小说| 久久香蕉激情| 又黄又粗又硬又大视频| 午夜福利,免费看| 亚洲欧美激情综合另类| 国产精品99久久99久久久不卡| 亚洲精品美女久久av网站| 操美女的视频在线观看| 亚洲午夜精品一区,二区,三区| 岛国视频午夜一区免费看| 国产国语露脸激情在线看| 亚洲色图av天堂| 一a级毛片在线观看| 中亚洲国语对白在线视频| 久久午夜综合久久蜜桃| 亚洲精品在线观看二区| 中亚洲国语对白在线视频| 国产亚洲精品综合一区在线观看 | 人人妻人人爽人人添夜夜欢视频| av福利片在线| 国产亚洲精品综合一区在线观看 | 啦啦啦免费观看视频1| 高清毛片免费观看视频网站 | ponron亚洲| 国产伦人伦偷精品视频| 国产av一区在线观看免费| 欧美+亚洲+日韩+国产| 啪啪无遮挡十八禁网站| 国产成人免费无遮挡视频| 欧美在线黄色| 国产色视频综合| 欧美乱色亚洲激情| 国产精品免费一区二区三区在线| 国产97色在线日韩免费| 国产视频一区二区在线看| 美女高潮到喷水免费观看| av网站在线播放免费| 每晚都被弄得嗷嗷叫到高潮| 免费搜索国产男女视频| 久久精品亚洲熟妇少妇任你| 久久精品影院6| 搡老乐熟女国产| 无遮挡黄片免费观看| 夜夜爽天天搞| xxxhd国产人妻xxx| 两个人免费观看高清视频| 午夜精品久久久久久毛片777| av电影中文网址| 91麻豆精品激情在线观看国产 | 妹子高潮喷水视频| 在线看a的网站| 精品国产乱子伦一区二区三区| 精品久久久久久,| 久久久久国内视频| 久久国产亚洲av麻豆专区| 在线观看一区二区三区| 国产av又大| 欧美人与性动交α欧美软件| 久久久水蜜桃国产精品网| svipshipincom国产片| 久久欧美精品欧美久久欧美| 一边摸一边抽搐一进一小说| 9191精品国产免费久久| 人人妻人人爽人人添夜夜欢视频| 国产精品亚洲av一区麻豆| 欧美亚洲日本最大视频资源| 91av网站免费观看| 两个人看的免费小视频| 日韩成人在线观看一区二区三区| 国产黄色免费在线视频| 一级毛片精品| 日韩欧美一区视频在线观看| 久久久国产成人免费| 亚洲欧美一区二区三区久久| 长腿黑丝高跟| 成人18禁高潮啪啪吃奶动态图| 纯流量卡能插随身wifi吗| 久久精品91无色码中文字幕| 可以在线观看毛片的网站| 午夜福利欧美成人| 首页视频小说图片口味搜索| 少妇粗大呻吟视频| 亚洲五月天丁香| 成人精品一区二区免费| 国产精品久久久人人做人人爽| 日本vs欧美在线观看视频| 在线观看日韩欧美| 亚洲 欧美一区二区三区| 免费在线观看完整版高清| 级片在线观看| av在线播放免费不卡| 成熟少妇高潮喷水视频| av在线播放免费不卡| 十分钟在线观看高清视频www| 悠悠久久av| 亚洲五月天丁香| 免费在线观看视频国产中文字幕亚洲| 长腿黑丝高跟| 国产高清videossex| 色在线成人网| 视频区图区小说| 妹子高潮喷水视频| 精品福利永久在线观看| 高清欧美精品videossex| 国产亚洲欧美在线一区二区| 正在播放国产对白刺激| 亚洲成人精品中文字幕电影 | 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩av久久| 黄色视频,在线免费观看| 国产在线精品亚洲第一网站| 男人舔女人下体高潮全视频| 国产精品美女特级片免费视频播放器 | 免费观看精品视频网站| 亚洲欧美日韩无卡精品| 一进一出抽搐动态| 在线观看www视频免费| 麻豆av在线久日| 国产乱人伦免费视频| 交换朋友夫妻互换小说| 无限看片的www在线观看| 99国产综合亚洲精品| 91精品国产国语对白视频| 99国产精品免费福利视频| 久久热在线av| 免费在线观看完整版高清| xxx96com| 国产精品久久久人人做人人爽| 欧美日韩亚洲国产一区二区在线观看| 精品久久久久久,| 成人国语在线视频| 国产成人系列免费观看| 十分钟在线观看高清视频www| 国产精品影院久久| 99精品久久久久人妻精品| 热re99久久精品国产66热6| 久久九九热精品免费| 国产亚洲精品久久久久久毛片| 亚洲黑人精品在线| 亚洲国产精品合色在线| 久久中文字幕一级| 女人爽到高潮嗷嗷叫在线视频| 国产成人精品在线电影| 99国产精品99久久久久| 亚洲视频免费观看视频| 国产精品久久电影中文字幕| 精品国产乱码久久久久久男人| 国产亚洲精品综合一区在线观看 | 亚洲精品一区av在线观看| 免费搜索国产男女视频| 无限看片的www在线观看| 天天添夜夜摸| 亚洲午夜理论影院| 久久精品人人爽人人爽视色| 中文字幕av电影在线播放| 日韩三级视频一区二区三区| 免费在线观看日本一区| 人人妻人人爽人人添夜夜欢视频| 久久热在线av| 欧美黄色片欧美黄色片| www.熟女人妻精品国产| 国产精品国产av在线观看| √禁漫天堂资源中文www| 日本欧美视频一区| 亚洲欧美精品综合久久99| 亚洲精品一区av在线观看| 亚洲欧美激情综合另类| 涩涩av久久男人的天堂| 国产区一区二久久| 午夜影院日韩av| 日本wwww免费看| av免费在线观看网站| 午夜福利在线观看吧| 在线观看免费高清a一片| 亚洲精品国产精品久久久不卡| 久久国产亚洲av麻豆专区| 久久久久久久午夜电影 | 老汉色av国产亚洲站长工具| 高潮久久久久久久久久久不卡| 亚洲久久久国产精品| 亚洲情色 制服丝袜| 99精品欧美一区二区三区四区| 久久午夜综合久久蜜桃| 亚洲激情在线av| 久久午夜综合久久蜜桃| 亚洲视频免费观看视频| 国产片内射在线| 久热这里只有精品99| √禁漫天堂资源中文www| av中文乱码字幕在线| 首页视频小说图片口味搜索| 亚洲国产精品sss在线观看 | 又黄又爽又免费观看的视频| 国产成人精品久久二区二区91| www日本在线高清视频| 日日干狠狠操夜夜爽| 午夜激情av网站| 国产成人影院久久av| 久久精品91蜜桃| 久久香蕉激情| 大型黄色视频在线免费观看| 免费在线观看亚洲国产| 午夜影院日韩av| 免费不卡黄色视频| 亚洲自拍偷在线| 国产成人av教育| 亚洲精品在线观看二区| 精品一区二区三区四区五区乱码| 热re99久久精品国产66热6| 大香蕉久久成人网| 夜夜看夜夜爽夜夜摸 | 正在播放国产对白刺激| 他把我摸到了高潮在线观看| 我的亚洲天堂| 美女高潮到喷水免费观看| 国产精品一区二区免费欧美| 久久精品影院6| 中文字幕高清在线视频| 久久久精品欧美日韩精品| 91在线观看av| 亚洲狠狠婷婷综合久久图片| 欧美日韩国产mv在线观看视频| 黄频高清免费视频| 国产精品99久久99久久久不卡| 一边摸一边抽搐一进一小说| 亚洲精品国产区一区二| av福利片在线| 婷婷精品国产亚洲av在线| 成年女人毛片免费观看观看9| 久久亚洲精品不卡| xxx96com| 日韩三级视频一区二区三区| 国产成人精品久久二区二区免费| 国内毛片毛片毛片毛片毛片| 超碰成人久久| 视频在线观看一区二区三区| 久久精品91蜜桃| 无人区码免费观看不卡| 国产成人系列免费观看| 黑人猛操日本美女一级片| 亚洲成av片中文字幕在线观看|