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

    基于Rhino精準曲面建模的常村煤礦原巖地應(yīng)力場反演方法

    2021-11-16 03:25:36祝凌甫伊丙鼎曲秋揚
    礦業(yè)安全與環(huán)保 2021年5期
    關(guān)鍵詞:常村應(yīng)力場主應(yīng)力

    祝凌甫,伊丙鼎,曲秋揚,胡 濱

    (1.中煤科工開采研究院有限公司,北京 100013; 2.煤炭科學(xué)研究總院開采研究分院,北京 100013;3.煤炭資源高效開采與潔凈利用國家重點實驗室,北京 100013)

    煤礦井下各種地質(zhì)構(gòu)造復(fù)雜,原巖地應(yīng)力場分布極不均勻,并且受限于煤礦井下有限的操作環(huán)境和高昂的測試費用,很難通過實測的方式獲取整個煤礦區(qū)的地應(yīng)力場。而采用數(shù)值模擬的方法可以對煤礦區(qū)地應(yīng)力場進行地應(yīng)力反演,進而獲取整個煤礦區(qū)的地應(yīng)力場分布,該方法可以為煤礦井下開采和支護等現(xiàn)場實踐提供更加全面和可靠的依據(jù)。

    國內(nèi)外學(xué)者在地應(yīng)力反演和分析方面做了大量研究[1-11]??导t普等[12]針對華亭礦區(qū)的大型向斜構(gòu)造附近的地應(yīng)力場進行了反演;蘇生瑞[13]通過離散元數(shù)值模擬反演了斷層附近的地應(yīng)力場;張勇慧等[14]采用數(shù)值模擬的方法對水電站地下廠房地質(zhì)模型進行了計算,反演出地下廠房及周邊應(yīng)力場的分布情況;孫猛[15]采用有限元軟件ANSYS和有限差分軟件FLAC3D反演了平頂山礦區(qū)原巖地應(yīng)力場;張寧博[16]通過數(shù)值模擬反演了東榮礦的地應(yīng)力場,同時分析了斷層對應(yīng)力場和沖擊地壓發(fā)生的影響。

    地應(yīng)力反演結(jié)果的精確性取決于很多因素,包括地質(zhì)模型與現(xiàn)場實際的相似程度、數(shù)值模擬過程中各層巖石力學(xué)參數(shù)賦值的準確性,以及模型計算過程中載荷施加情況等。其中,在三維地質(zhì)建模方面,通過煤礦區(qū)大量的鉆孔柱狀圖對不同巖層進行標識和擬合為目前公認準確度較高的方法,常用的三維建模軟件有Rhino、Surfer和ANSYS等。一般煤巖層均為起伏不平的曲面,而部分建模軟件在建立曲面方面效果一般,無法準確反映實際地層延伸情況,因此,選擇一個能夠精準建立曲面模型的建模軟件對于地應(yīng)力反演的準確度至關(guān)重要。

    Rhino建模軟件作為專業(yè)的建模軟件具有以下幾個特點:①Rhino建模方式多樣化,效果好,其還帶有分析功能和渲染功能,基本滿足了3D建模的所有需求;②自動快速形成曲面,曲面精度高,與現(xiàn)場實際貼合性強;③兼容性佳,其默認3D模型保存格式是3dm,也可以多種格式保存,Rhino軟件經(jīng)過轉(zhuǎn)換,可以把三維文件轉(zhuǎn)換成二維圖形和線條文件;④可操作性強,界面繪圖指令豐富,既可以采用編程的方式繪圖,還可以采用直接手繪的方式作圖;并且Rhino建模軟件建立的三維地質(zhì)模型與FLAC3D、UDEC及3DEC等數(shù)值模擬軟件具有非常好的兼容性,可以直接將三維模型導(dǎo)入這些軟件進行數(shù)值計算。

    選取潞安礦區(qū)常村煤礦+470 m水平煤巖層為研究對象,采用Rhino軟件構(gòu)建其三維地質(zhì)模型,采用有限差分數(shù)值模擬軟件FLAC3D進行地應(yīng)力場的計算和反演,最后將地應(yīng)力反演結(jié)果與現(xiàn)場實測結(jié)果進行對比,評估地應(yīng)力反演的準確性,以及Rhino建模的適用性。

    1 Rhino三維地質(zhì)模型及數(shù)值模型的建立

    1.1 工程概況

    常村煤礦主采煤層平均厚度為5.98 m。礦井分為2個水平,筆者只針對+470 m水平開展建模和計算。該水平埋深在470~600 m內(nèi)。煤層直接頂為粉砂巖,基本頂為中粒砂巖,老底為泥巖。井田內(nèi)地質(zhì)構(gòu)造以褶曲為主,小斷層相對比較多,大斷層主要是aF9斷層和FB-1斷層,礦井三維地質(zhì)模型的建立只考慮了這2個斷層的存在。

    1.2 Rhino三維地質(zhì)模型的建立思路

    由于煤礦井下地層結(jié)構(gòu)的不可見性,無法了解從地表到煤層這部分地層的準確分布情況,因此采用一種簡化的辦法模擬地表到煤層的地層分布結(jié)構(gòu)??偣卜譃?個步驟:首先對煤礦鉆孔柱狀圖進行簡化,選取不同標志層作為關(guān)鍵地層分界面;然后,在地質(zhì)鉆孔資料中將各標志層的標高進行標識,按照標高將這些標識點連接并擬合成曲面,該曲面與地層的延伸方向相似度比較高,可以表示數(shù)值模擬模型中地層的分布。

    分析常村煤礦的鉆孔柱狀圖等地質(zhì)資料,從鉆孔資料的精確性和分布的均勻性兩個方面考慮,選取了15個地質(zhì)鉆孔,各鉆孔標識點分布情況如圖1所示。按照上述的簡化思路和曲面擬合的原則,將選取的地質(zhì)鉆孔的標志層進行標識,連接擬合出精度最好的曲面來表示巖層的分布。

    圖1 各鉆孔標識點分布圖

    將擬合出的曲面地層合并,通過Rhinoceros軟件建立常村煤礦+470 m水平煤巖層的三維地質(zhì)模型,如圖2所示,其長×寬×高=7 300 m×3 800 m×750 m。常村煤礦煤層賦存最大深度為750 m,即褶曲構(gòu)造的核心部位,為了更加準確模擬整個常村煤礦的煤層內(nèi)部地應(yīng)力分布情況,選擇750 m作為所建立的地質(zhì)模型的高度。

    圖2 Rhinoceros軟件構(gòu)建的常村煤礦三維地質(zhì)模型

    考慮反演的精度和計算機的能力,將地質(zhì)模型進行單元劃分,劃分精度選取為40 m左右,總共生成了281 733個單元,180 026個節(jié)點。將劃分好網(wǎng)格的模型導(dǎo)入FLAC3D數(shù)值模擬軟件,生成FLAC3D數(shù)值計算模型,如圖3所示。本構(gòu)關(guān)系采用Mohr-Coulomb屈服準則。由于斷層帶附近巖石承載能力基本為零,斷層帶的模擬采取巖石力學(xué)參數(shù)弱化的方式實現(xiàn)。

    圖3 常村煤礦FLAC3D數(shù)值計算模型

    2 數(shù)值模型的計算

    2.1 計算原理

    煤巖體水平應(yīng)力一般由兩部分組成:一是由自重應(yīng)力引起的水平應(yīng)力分量;二是由構(gòu)造運動引起的水平應(yīng)力分量[1]。最大水平主應(yīng)力計算公式如下:

    σH=σHg+σHt

    (1)

    (2)

    (3)

    式中:σH為最大水平主應(yīng)力;σHg為由自重應(yīng)力產(chǎn)生的水平應(yīng)力分量;σHt為由構(gòu)造應(yīng)力場產(chǎn)生的水平應(yīng)力分量;σV為垂直應(yīng)力;μ為泊松比;εHt為構(gòu)造應(yīng)力引起的水平應(yīng)變;E為彈性模量。

    選取某個地應(yīng)力測點為初始點進行應(yīng)力反演(即S3膠帶下山測點,其最大、垂直和最小主應(yīng)力分別為12.87、9.33、6.71 MPa,最大主應(yīng)力為水平構(gòu)造應(yīng)力,為走滑型應(yīng)力狀態(tài))。由式(2)可知,由重力側(cè)壓效應(yīng)所產(chǎn)生的水平應(yīng)力分量為4.00 MPa;由式(1)計算可得,由構(gòu)造應(yīng)力場產(chǎn)生的水平應(yīng)力分量為8.87 MPa,構(gòu)造應(yīng)力引起的水平應(yīng)變εHt為1.66×10-3。因此,基于以上計算過程可知,為了還原地應(yīng)力場原始賦存狀況,需要施加給模型側(cè)邊的水平應(yīng)變?yōu)?.66×10-3。

    2.2 計算過程

    常村煤礦+470 m水平三維模型總共分為5個邊界,在應(yīng)力場反演過程中需要分別對其施加不同的邊界載荷條件或約束條件,這5個邊界在計算模型中的數(shù)學(xué)坐標,以及表示方法如圖4 所示。

    圖4 模型邊界數(shù)學(xué)坐標表示方法

    模擬過程分2個階段:①自重應(yīng)力場演化。固定模型的側(cè)面和底部位移,施加重力加速度,模型在自重作用下達到平衡。②構(gòu)造應(yīng)力場演化。解除模型側(cè)面的法向水平位移固定邊界條件,施加速度邊界條件,計算直至平衡。應(yīng)力場演化模擬的邊界條件如圖5所示。

    圖5 應(yīng)力場演化模擬的邊界條件

    3 地應(yīng)力反演結(jié)果及分析

    基于FLAC3D數(shù)值模擬計算結(jié)果,為了對比同一埋深下的地應(yīng)力變化情況,采用固定埋深這一變量,輸出不同坐標的地應(yīng)力值,繪制三向主應(yīng)力折線圖,分析煤層內(nèi)部和基本頂內(nèi)應(yīng)力場的分布規(guī)律。依據(jù)模型數(shù)學(xué)坐標,設(shè)置了4條勘探線來監(jiān)測礦區(qū)內(nèi)部應(yīng)力場的變化。具體為:在沿x方向的煤層和基本頂內(nèi)部,每隔300 m的距離設(shè)置程序輸出三向主應(yīng)力量值,繪制煤層和基本頂內(nèi)部y=600 m、y=2 500 m、y=3 500 m三向主應(yīng)力沿x軸的變化折線圖,分析主應(yīng)力的變化趨勢;在沿y方向的煤層和基本頂內(nèi)部,每隔100 m的距離設(shè)置程序輸出三向主應(yīng)力量值,繪制煤層和基本頂內(nèi)部x=7 000 m三向主應(yīng)力沿y軸的變化折線圖,分析主應(yīng)力的變化趨勢。各主應(yīng)力變化情況如圖6~9所示。通過該方法,定量分析了常村煤礦+470 m水平井下煤巖應(yīng)力場的分布及演化。同時,由于本文構(gòu)建的三維地質(zhì)模型為曲面模型,更接近地層真實賦存情況,地層起伏等情況導(dǎo)致在遠離大型褶曲和斷層的地理區(qū)域內(nèi),地應(yīng)力反演的數(shù)據(jù)也會隨地層起伏在一定范圍內(nèi)波動。

    圖6 y=600 m、z=470 m勘探線基本頂主應(yīng)力隨x坐標變化

    圖7 x=7 000 m、z=470 m勘探線主應(yīng)力隨y坐標變化

    圖8 y=2 500 m、z=470 m勘探線主應(yīng)力隨x坐標變化

    圖9 y=3 500 m、z=470 m勘探線主應(yīng)力隨x坐標變化

    通過分析可以得出,常村煤礦+470 m水平基本頂最大主應(yīng)力量值(構(gòu)造應(yīng)力)穩(wěn)定在10~18 MPa內(nèi)。遠離斷層區(qū)域內(nèi),在煤巖層褶曲構(gòu)造附近,最大主應(yīng)力達到10 MPa和18 MPa兩個極值,其中,向斜構(gòu)造軸部的應(yīng)力值最高,達到18 MPa;其他范圍內(nèi)的最大主應(yīng)力值基本穩(wěn)定在12 MPa,最小主應(yīng)力值穩(wěn)定在6 MPa。因此,在向斜構(gòu)造區(qū)域內(nèi)進行開采活動時,一定要密切監(jiān)測頂板應(yīng)力值變化情況,當應(yīng)力值產(chǎn)生突變或者遠高于其他區(qū)域時,應(yīng)采用加強支護或者卸壓的方法維護采場和巷道的圍巖穩(wěn)定。在斷層附近區(qū)域,地應(yīng)力產(chǎn)生突變,斷層端部會產(chǎn)生應(yīng)力集中,這也符合斷裂力學(xué)中裂紋擴展力學(xué)機制。但在斷層核部則進入完全松弛區(qū),地應(yīng)力值降低為 6 MPa 左右或者進入完全卸壓區(qū),巖石承載力基本為0,地應(yīng)力值接近于0。因此,在斷層區(qū)域內(nèi)進行開采活動時,一定要密切監(jiān)測應(yīng)力變化和巖石強度,一般情況下,應(yīng)對斷層附近區(qū)域采用加強支護或者預(yù)注漿的方法保證工作面圍巖的強度和穩(wěn)定性。

    煤炭科學(xué)研究總院開采研究分院經(jīng)過10余年的現(xiàn)場地應(yīng)力測量工作,在潞安礦區(qū)常村煤礦進行了17個不同地點的地應(yīng)力測量,獲取了17組地應(yīng)力測量數(shù)據(jù),實測數(shù)據(jù)見表1。

    表1 常村煤礦地應(yīng)力實測數(shù)據(jù)

    備注:①由于地應(yīng)力現(xiàn)場實測過程中未列出地應(yīng)力測點的地理坐標,無法在地應(yīng)力反演結(jié)果中準確定位,未給出每個地應(yīng)力測點處的地應(yīng)力反演結(jié)果;② 9號測點的異常性有可能是由不可避免的測量誤差、處于地質(zhì)異常區(qū)等因素導(dǎo)致。

    經(jīng)過對比分析實測數(shù)據(jù)與通過數(shù)值模擬獲取的地應(yīng)力反演結(jié)果,可以發(fā)現(xiàn),潞安礦區(qū)常村煤礦現(xiàn)場實測地應(yīng)力最大主應(yīng)力分布在10.17~13.57 MPa內(nèi),平均值為12.10 MPa;最小主應(yīng)力分布在5.30~7.30 MPa內(nèi),平均值為6.40 MPa;垂直主應(yīng)力分布在8.03~11.54 MPa內(nèi),平均值為9.20 MPa。數(shù)值模擬計算結(jié)果顯示,在遠離地質(zhì)構(gòu)造(斷層和褶曲等)區(qū)域內(nèi),最大主應(yīng)力基本在12 MPa左右波動,最小主應(yīng)力基本在6 MPa左右波動,中間主應(yīng)力變化范圍相對比較大,但也基本在9 MPa左右上下浮動,與地應(yīng)力實測的平均值和變化范圍一致性非常高。綜上分析,兩者均顯示,常村煤礦+470 m水平最大水平主應(yīng)力穩(wěn)定在12 MPa左右,最小水平主應(yīng)力穩(wěn)定在6 MPa左右,垂直主應(yīng)力穩(wěn)定在9 MPa左右,屬于走滑型應(yīng)力狀態(tài)。因此,在地應(yīng)力量值方面,地應(yīng)力實測數(shù)據(jù)和地應(yīng)力反演數(shù)據(jù)具有比較好的一致性。在三維地質(zhì)模型建立與現(xiàn)場實際貼合度比較高的條件下,采用數(shù)值模擬的方法進行地應(yīng)力反演具有可行性,能夠真實反映煤礦井下原巖地應(yīng)力場的分布,并且可以彌補在地質(zhì)構(gòu)造附近很難獲取準確的地應(yīng)力實測數(shù)據(jù)的不足。

    4 結(jié)論

    1)利用Rhino建模軟件可以比較精確地對煤礦區(qū)的地質(zhì)模型進行構(gòu)建,以曲面代替巖層更加真實地反映了煤礦井下地層的實際分布情況,褶曲等地質(zhì)構(gòu)造可以被真實體現(xiàn),原巖地應(yīng)力場的數(shù)值模擬反演結(jié)果更加準確可靠。

    2)大型地質(zhì)構(gòu)造區(qū)域原巖地應(yīng)力場與其他區(qū)域有較大的差異,向斜軸部原巖地應(yīng)力場集中程度最大,從斷層中部到斷層端部再到遠離斷層區(qū)域地應(yīng)力場呈現(xiàn)“小—大—原”的分布規(guī)律。

    3)通過對比常村煤礦現(xiàn)場實測的17個地應(yīng)力測點量值與地應(yīng)力反演量值,常村煤礦+470 m水平最大主應(yīng)力穩(wěn)定在12 MPa左右,最小主應(yīng)力穩(wěn)定在6 MPa左右,垂直主應(yīng)力穩(wěn)定在9 MPa左右,屬于走滑型應(yīng)力狀態(tài),地應(yīng)力實測結(jié)果和地應(yīng)力反演結(jié)果具有高度的一致性。

    猜你喜歡
    常村應(yīng)力場主應(yīng)力
    潞安環(huán)保能源開發(fā)股份有限公司常村煤礦企業(yè)簡介
    煤(2024年4期)2024-04-11 08:58:20
    節(jié)能減潞安環(huán)保能源開發(fā)股份有限公司常村煤礦企業(yè)簡介耗宣傳
    煤(2024年3期)2024-03-13 09:40:30
    復(fù)合斷層對地應(yīng)力的影響研究
    地音監(jiān)測技術(shù)在常村煤礦沖擊地壓預(yù)警中的應(yīng)用
    中國煤炭(2016年1期)2016-05-17 06:11:32
    鋁合金多層多道窄間隙TIG焊接頭應(yīng)力場研究
    焊接(2016年9期)2016-02-27 13:05:22
    深部沿空巷道圍巖主應(yīng)力差演化規(guī)律與控制
    考慮中主應(yīng)力后對隧道圍巖穩(wěn)定性的影響
    考慮斷裂破碎帶的丹江口庫區(qū)地應(yīng)力場與水壓應(yīng)力場耦合反演及地震預(yù)測
    基于位移相關(guān)法的重復(fù)壓裂裂縫尖端應(yīng)力場研究
    斷塊油氣田(2014年5期)2014-03-11 15:33:49
    岸坡應(yīng)力場及卸荷帶劃分量化指標研究
    国产黄色视频一区二区在线观看| 亚洲va在线va天堂va国产| 18禁裸乳无遮挡动漫免费视频 | 在线观看三级黄色| a级毛片免费高清观看在线播放| 91在线精品国自产拍蜜月| 国产成人免费无遮挡视频| 欧美性猛交╳xxx乱大交人| 毛片一级片免费看久久久久| 18禁裸乳无遮挡免费网站照片| 国语对白做爰xxxⅹ性视频网站| 一区二区三区乱码不卡18| 成人漫画全彩无遮挡| 国产成人免费无遮挡视频| 少妇裸体淫交视频免费看高清| 亚洲成色77777| 白带黄色成豆腐渣| 天天一区二区日本电影三级| 国产精品.久久久| 国产熟女欧美一区二区| 国产乱人视频| 2021天堂中文幕一二区在线观| 3wmmmm亚洲av在线观看| 丰满人妻一区二区三区视频av| 国产男人的电影天堂91| 久久久精品欧美日韩精品| 国产成人免费无遮挡视频| 水蜜桃什么品种好| 免费看日本二区| 欧美日本视频| 亚洲成人精品中文字幕电影| 啦啦啦中文免费视频观看日本| 三级国产精品片| 小蜜桃在线观看免费完整版高清| 神马国产精品三级电影在线观看| 男女国产视频网站| 国产69精品久久久久777片| 看黄色毛片网站| 边亲边吃奶的免费视频| 丰满乱子伦码专区| 韩国av在线不卡| 少妇的逼水好多| 午夜福利高清视频| 我要看日韩黄色一级片| 亚洲成人一二三区av| 国产精品人妻久久久久久| 男女边摸边吃奶| 久久人人爽人人爽人人片va| 最后的刺客免费高清国语| 国产爽快片一区二区三区| 国产一区二区三区综合在线观看 | 亚洲国产欧美在线一区| 成年人午夜在线观看视频| 国产成人a区在线观看| 全区人妻精品视频| 2018国产大陆天天弄谢| 免费在线观看成人毛片| 在线精品无人区一区二区三 | 美女xxoo啪啪120秒动态图| 亚洲精品成人av观看孕妇| 国产淫语在线视频| freevideosex欧美| 最近中文字幕高清免费大全6| 超碰97精品在线观看| 日韩在线高清观看一区二区三区| 三级男女做爰猛烈吃奶摸视频| 看免费成人av毛片| 亚洲丝袜综合中文字幕| 18禁动态无遮挡网站| 狠狠精品人妻久久久久久综合| 少妇 在线观看| 秋霞在线观看毛片| 久久精品国产鲁丝片午夜精品| 一级毛片aaaaaa免费看小| 97超碰精品成人国产| 黄色日韩在线| 亚洲精品乱码久久久久久按摩| 精品国产乱码久久久久久小说| 亚洲av免费在线观看| 国产综合精华液| 亚洲精品日本国产第一区| 如何舔出高潮| 国产在线男女| 波野结衣二区三区在线| 免费观看av网站的网址| kizo精华| 男女无遮挡免费网站观看| 可以在线观看毛片的网站| 新久久久久国产一级毛片| 国产精品麻豆人妻色哟哟久久| 一级av片app| 成年人午夜在线观看视频| 国产91av在线免费观看| 日本一本二区三区精品| 亚洲国产精品999| 亚洲久久久久久中文字幕| 成人综合一区亚洲| 色综合色国产| 汤姆久久久久久久影院中文字幕| 免费电影在线观看免费观看| 少妇裸体淫交视频免费看高清| 五月天丁香电影| 亚洲成人一二三区av| 99久久精品一区二区三区| 国产精品国产av在线观看| 亚洲国产高清在线一区二区三| 在线观看三级黄色| 免费看日本二区| h日本视频在线播放| 免费av毛片视频| 亚洲综合精品二区| 天堂网av新在线| 午夜福利高清视频| 69av精品久久久久久| 日本三级黄在线观看| 99视频精品全部免费 在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产 一区 欧美 日韩| 在线a可以看的网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久久久久免费av| 各种免费的搞黄视频| 好男人在线观看高清免费视频| 美女被艹到高潮喷水动态| 亚洲精品亚洲一区二区| 少妇猛男粗大的猛烈进出视频 | 一边亲一边摸免费视频| 国产精品伦人一区二区| 国产av码专区亚洲av| 欧美精品国产亚洲| 久久国内精品自在自线图片| 超碰97精品在线观看| 热99国产精品久久久久久7| 最后的刺客免费高清国语| 午夜福利视频1000在线观看| 最近手机中文字幕大全| 亚洲人与动物交配视频| 涩涩av久久男人的天堂| 伊人久久精品亚洲午夜| 日韩一本色道免费dvd| 亚洲av免费在线观看| 精品久久久久久久末码| 免费看日本二区| 欧美老熟妇乱子伦牲交| 国产女主播在线喷水免费视频网站| 国产一区二区在线观看日韩| 综合色av麻豆| 国产精品嫩草影院av在线观看| 国产大屁股一区二区在线视频| 色综合色国产| 国产一区二区亚洲精品在线观看| 狠狠精品人妻久久久久久综合| 高清午夜精品一区二区三区| 久久影院123| 中文字幕亚洲精品专区| 欧美日韩一区二区视频在线观看视频在线 | freevideosex欧美| 色播亚洲综合网| 尤物成人国产欧美一区二区三区| 亚洲激情五月婷婷啪啪| 国模一区二区三区四区视频| 日韩在线高清观看一区二区三区| 最近最新中文字幕免费大全7| 下体分泌物呈黄色| 黄色日韩在线| 舔av片在线| 中文字幕制服av| 精品熟女少妇av免费看| 午夜免费鲁丝| 欧美97在线视频| 神马国产精品三级电影在线观看| 久久久久九九精品影院| 欧美+日韩+精品| 久久99热这里只有精品18| 国产精品嫩草影院av在线观看| 中国美白少妇内射xxxbb| 久久久久精品久久久久真实原创| 在线免费观看不下载黄p国产| 不卡视频在线观看欧美| 高清av免费在线| 日本黄色片子视频| 成人毛片60女人毛片免费| 亚洲精品国产色婷婷电影| 午夜精品国产一区二区电影 | av在线天堂中文字幕| 精品一区二区三区视频在线| 午夜福利网站1000一区二区三区| av网站免费在线观看视频| 国产精品一及| 国产视频首页在线观看| 97热精品久久久久久| 成人鲁丝片一二三区免费| 成人综合一区亚洲| 免费观看av网站的网址| 国产成人91sexporn| 啦啦啦中文免费视频观看日本| 亚洲精品第二区| av专区在线播放| 寂寞人妻少妇视频99o| 日韩制服骚丝袜av| 一区二区三区乱码不卡18| 亚洲美女搞黄在线观看| 成人亚洲精品一区在线观看 | 久久亚洲国产成人精品v| 国产精品国产三级专区第一集| 亚洲自偷自拍三级| 91精品国产九色| 天堂中文最新版在线下载 | 尤物成人国产欧美一区二区三区| 人妻夜夜爽99麻豆av| 视频中文字幕在线观看| 一级毛片aaaaaa免费看小| 少妇丰满av| 99久国产av精品国产电影| 尾随美女入室| 亚洲国产av新网站| 下体分泌物呈黄色| 美女脱内裤让男人舔精品视频| 国产精品精品国产色婷婷| 久久久久性生活片| 麻豆久久精品国产亚洲av| 午夜免费鲁丝| 看黄色毛片网站| 亚洲精品日本国产第一区| 人人妻人人看人人澡| 99久久精品热视频| 晚上一个人看的免费电影| 一级毛片黄色毛片免费观看视频| 丝袜美腿在线中文| 小蜜桃在线观看免费完整版高清| 大码成人一级视频| 亚洲av不卡在线观看| 国产免费一级a男人的天堂| 欧美日韩视频高清一区二区三区二| 亚洲人成网站高清观看| 久久精品国产自在天天线| 亚洲欧美中文字幕日韩二区| 色5月婷婷丁香| freevideosex欧美| 精品一区二区三区视频在线| 一个人看视频在线观看www免费| 国产一区二区三区av在线| 久久久色成人| 99视频精品全部免费 在线| 视频中文字幕在线观看| kizo精华| 久久久久久国产a免费观看| 高清毛片免费看| 久久女婷五月综合色啪小说 | 国产精品一及| 国产成人a∨麻豆精品| 男人添女人高潮全过程视频| 国产午夜福利久久久久久| 一级毛片 在线播放| 国产精品一区www在线观看| 波多野结衣巨乳人妻| 久久99精品国语久久久| 少妇猛男粗大的猛烈进出视频 | 97超视频在线观看视频| 国产黄色视频一区二区在线观看| 色5月婷婷丁香| 精品国产三级普通话版| 中文乱码字字幕精品一区二区三区| 欧美日本视频| 精品一区二区免费观看| 精品一区二区三区视频在线| 伦精品一区二区三区| www.色视频.com| 99久国产av精品国产电影| 一个人看视频在线观看www免费| 免费av毛片视频| 在线免费十八禁| kizo精华| 99久久精品热视频| 中文乱码字字幕精品一区二区三区| 亚洲av中文字字幕乱码综合| 婷婷色麻豆天堂久久| 亚洲精品影视一区二区三区av| 国产日韩欧美在线精品| 精品人妻视频免费看| 久久久久九九精品影院| 免费看不卡的av| 精品久久久久久久久亚洲| 欧美成人a在线观看| 人人妻人人看人人澡| 免费电影在线观看免费观看| 国产精品一二三区在线看| 国产熟女欧美一区二区| 视频中文字幕在线观看| 日日啪夜夜爽| 男人舔奶头视频| 中国国产av一级| 国产亚洲av片在线观看秒播厂| 欧美日韩一区二区视频在线观看视频在线 | 一区二区av电影网| 欧美成人一区二区免费高清观看| 国产午夜精品一二区理论片| 99久久精品热视频| 麻豆乱淫一区二区| 国产精品.久久久| 亚州av有码| 黄色视频在线播放观看不卡| 一个人看视频在线观看www免费| 五月玫瑰六月丁香| 少妇熟女欧美另类| 精品国产露脸久久av麻豆| 亚洲最大成人av| 亚洲精品乱码久久久v下载方式| 国产精品一区二区性色av| 天美传媒精品一区二区| 又爽又黄a免费视频| 汤姆久久久久久久影院中文字幕| 啦啦啦中文免费视频观看日本| 晚上一个人看的免费电影| 日韩国内少妇激情av| 最近最新中文字幕免费大全7| 免费黄色在线免费观看| 高清视频免费观看一区二区| 日韩伦理黄色片| 午夜福利视频1000在线观看| 秋霞伦理黄片| 精品久久久久久久久亚洲| 中文乱码字字幕精品一区二区三区| 日韩成人伦理影院| 日韩制服骚丝袜av| 国产视频内射| 日本-黄色视频高清免费观看| 日韩成人av中文字幕在线观看| 国产高清三级在线| 午夜日本视频在线| 久久精品久久精品一区二区三区| 妹子高潮喷水视频| 肉色欧美久久久久久久蜜桃| 观看美女的网站| 人人妻,人人澡人人爽秒播 | 日韩 欧美 亚洲 中文字幕| 老鸭窝网址在线观看| 国产亚洲最大av| 久久久精品区二区三区| 人体艺术视频欧美日本| 国产高清国产精品国产三级| 少妇 在线观看| 久久久久精品人妻al黑| 操美女的视频在线观看| 2021少妇久久久久久久久久久| 人妻 亚洲 视频| 久久精品久久久久久噜噜老黄| av不卡在线播放| 亚洲第一青青草原| 亚洲免费av在线视频| 天堂俺去俺来也www色官网| 日韩熟女老妇一区二区性免费视频| 日韩中文字幕视频在线看片| 999精品在线视频| 免费观看av网站的网址| 精品第一国产精品| 亚洲少妇的诱惑av| 成人手机av| 婷婷色综合www| 欧美人与善性xxx| 欧美日韩一级在线毛片| 亚洲三区欧美一区| 亚洲天堂av无毛| 999精品在线视频| 午夜av观看不卡| 韩国高清视频一区二区三区| 日本色播在线视频| 久久久久久久大尺度免费视频| 国产精品 国内视频| 曰老女人黄片| 丝袜人妻中文字幕| 韩国高清视频一区二区三区| 中文欧美无线码| 天堂俺去俺来也www色官网| 一区二区三区精品91| 一边摸一边抽搐一进一出视频| 久久国产精品大桥未久av| 欧美日韩国产mv在线观看视频| 黑丝袜美女国产一区| 欧美乱码精品一区二区三区| 国产av精品麻豆| 午夜91福利影院| 99精国产麻豆久久婷婷| 久久精品亚洲熟妇少妇任你| 大话2 男鬼变身卡| 一二三四在线观看免费中文在| 一级片免费观看大全| 水蜜桃什么品种好| 汤姆久久久久久久影院中文字幕| 性色av一级| 久久性视频一级片| 免费看av在线观看网站| 最黄视频免费看| 在线亚洲精品国产二区图片欧美| 日本vs欧美在线观看视频| 亚洲精品一区蜜桃| 国产成人免费观看mmmm| 日韩电影二区| av一本久久久久| 亚洲精华国产精华液的使用体验| 中国国产av一级| 尾随美女入室| 啦啦啦视频在线资源免费观看| 国产一区二区激情短视频 | 亚洲精品国产色婷婷电影| 国产又爽黄色视频| 成人黄色视频免费在线看| 国产精品久久久久久精品电影小说| 欧美黑人欧美精品刺激| bbb黄色大片| 飞空精品影院首页| 中文天堂在线官网| 黄色怎么调成土黄色| avwww免费| 国产精品一国产av| 亚洲四区av| 久久久国产精品麻豆| 黄片播放在线免费| 欧美97在线视频| 大码成人一级视频| 亚洲伊人色综图| 欧美日韩视频精品一区| 一二三四在线观看免费中文在| www.av在线官网国产| 欧美国产精品一级二级三级| 国产野战对白在线观看| 亚洲色图 男人天堂 中文字幕| 下体分泌物呈黄色| 男女之事视频高清在线观看 | 中文字幕人妻丝袜制服| 久久热在线av| 日韩电影二区| 精品久久蜜臀av无| 亚洲图色成人| 免费在线观看黄色视频的| 国产无遮挡羞羞视频在线观看| 久久久亚洲精品成人影院| 大片免费播放器 马上看| 大香蕉久久网| 国产一区二区 视频在线| 精品视频人人做人人爽| 日韩一卡2卡3卡4卡2021年| 18禁动态无遮挡网站| 男女之事视频高清在线观看 | 亚洲欧美一区二区三区黑人| 日韩av免费高清视频| 丰满饥渴人妻一区二区三| 性高湖久久久久久久久免费观看| 天堂8中文在线网| 在线精品无人区一区二区三| 制服人妻中文乱码| 成人午夜精彩视频在线观看| 亚洲欧洲国产日韩| 亚洲精品成人av观看孕妇| 亚洲人成网站在线观看播放| 亚洲国产精品一区二区三区在线| 国产成人免费观看mmmm| 男男h啪啪无遮挡| 18禁动态无遮挡网站| 欧美黑人欧美精品刺激| 婷婷色麻豆天堂久久| 午夜久久久在线观看| 五月天丁香电影| videos熟女内射| 久久精品人人爽人人爽视色| 制服诱惑二区| 在线观看国产h片| 国产精品久久久久久精品电影小说| 日本91视频免费播放| videos熟女内射| 国产xxxxx性猛交| 久久国产精品大桥未久av| 国产福利在线免费观看视频| 亚洲色图 男人天堂 中文字幕| √禁漫天堂资源中文www| 女人爽到高潮嗷嗷叫在线视频| 另类精品久久| 国产精品香港三级国产av潘金莲 | 亚洲四区av| 国产乱人偷精品视频| 另类亚洲欧美激情| 丝袜美腿诱惑在线| 日韩不卡一区二区三区视频在线| 黑人欧美特级aaaaaa片| 欧美成人午夜精品| 18禁国产床啪视频网站| 麻豆av在线久日| 国产高清国产精品国产三级| 亚洲精品成人av观看孕妇| 国产成人午夜福利电影在线观看| 一级,二级,三级黄色视频| 国产男女内射视频| 国产片特级美女逼逼视频| tube8黄色片| 国产伦人伦偷精品视频| 色吧在线观看| 亚洲国产最新在线播放| 熟妇人妻不卡中文字幕| 九草在线视频观看| 国产国语露脸激情在线看| 在线观看人妻少妇| 精品国产露脸久久av麻豆| 女人精品久久久久毛片| av电影中文网址| 欧美精品高潮呻吟av久久| 狠狠婷婷综合久久久久久88av| 色吧在线观看| 亚洲自偷自拍图片 自拍| tube8黄色片| 满18在线观看网站| 午夜激情久久久久久久| 老司机在亚洲福利影院| 久久久精品免费免费高清| 国产精品 国内视频| 一级毛片电影观看| 国产精品免费视频内射| 亚洲视频免费观看视频| 香蕉国产在线看| av在线播放精品| 97在线人人人人妻| 狠狠婷婷综合久久久久久88av| 亚洲精品乱久久久久久| 欧美人与性动交α欧美软件| 国产成人系列免费观看| 满18在线观看网站| 夫妻午夜视频| 97精品久久久久久久久久精品| 亚洲第一青青草原| 亚洲天堂av无毛| 国产免费一区二区三区四区乱码| 高清在线视频一区二区三区| 成年美女黄网站色视频大全免费| av网站在线播放免费| 久热爱精品视频在线9| 久久av网站| 欧美日韩一级在线毛片| 亚洲婷婷狠狠爱综合网| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品 国内视频| 人人澡人人妻人| 亚洲国产成人一精品久久久| 久久久国产精品麻豆| 777米奇影视久久| 一二三四在线观看免费中文在| 高清视频免费观看一区二区| 在线观看免费午夜福利视频| 亚洲精品国产av成人精品| 日韩一区二区视频免费看| 日本av免费视频播放| 色播在线永久视频| 精品国产一区二区三区四区第35| 麻豆av在线久日| 久久精品久久精品一区二区三区| 国产野战对白在线观看| 久久免费观看电影| 啦啦啦在线免费观看视频4| 国产在线一区二区三区精| 欧美精品人与动牲交sv欧美| 欧美日韩视频高清一区二区三区二| 韩国精品一区二区三区| 人成视频在线观看免费观看| 国产视频首页在线观看| 熟妇人妻不卡中文字幕| 国产福利在线免费观看视频| 国产日韩一区二区三区精品不卡| 人人妻,人人澡人人爽秒播 | 国产老妇伦熟女老妇高清| 免费在线观看完整版高清| 一本—道久久a久久精品蜜桃钙片| 水蜜桃什么品种好| 性高湖久久久久久久久免费观看| 大码成人一级视频| 如何舔出高潮| 国产免费福利视频在线观看| 久久99热这里只频精品6学生| 在线观看一区二区三区激情| 观看av在线不卡| av不卡在线播放| 不卡视频在线观看欧美| 日韩成人av中文字幕在线观看| 捣出白浆h1v1| 亚洲成人手机| 午夜福利在线免费观看网站| 精品少妇内射三级| 久久鲁丝午夜福利片| 黄片播放在线免费| 亚洲成人国产一区在线观看 | 美女福利国产在线| 蜜桃在线观看..| 久久这里只有精品19| 亚洲五月色婷婷综合| 波野结衣二区三区在线| 两个人看的免费小视频| 免费少妇av软件| 日韩伦理黄色片| 久久人人爽人人片av| 丝袜美腿诱惑在线| av在线播放精品| 热99久久久久精品小说推荐| 极品人妻少妇av视频| 蜜桃国产av成人99| 欧美在线一区亚洲| 午夜日韩欧美国产| 久久精品国产亚洲av涩爱| 日韩av免费高清视频| 午夜福利视频在线观看免费| 美女主播在线视频| 一边摸一边做爽爽视频免费| 最近最新中文字幕大全免费视频 | 97人妻天天添夜夜摸| 久久久亚洲精品成人影院| 亚洲精品第二区| 亚洲精品日韩在线中文字幕| 亚洲三区欧美一区| 香蕉国产在线看|