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

    黏彈介質(zhì)槽波超前探測(cè)數(shù)值模擬與波場(chǎng)分析

    2023-11-06 15:55:36何帝秀姬廣忠焦文杰張亞偉
    煤田地質(zhì)與勘探 2023年10期
    關(guān)鍵詞:波組橫波夾角

    何帝秀,姬廣忠,焦文杰,張亞偉,于 坤

    (安徽理工大學(xué) 地球與環(huán)境學(xué)院,安徽 淮南 232001)

    巷道掘進(jìn)作為煤礦生產(chǎn)工作的重要環(huán)節(jié)之一,若煤巷掘進(jìn)工作面(以下簡(jiǎn)稱工作面)前方存在異常地質(zhì)構(gòu)造如斷層、陷落柱、采空區(qū)等,很容易影響地下煤礦開采工作的安全性,嚴(yán)重時(shí)會(huì)引發(fā)礦井水害、頂板垮落、煤與瓦斯突出等地質(zhì)災(zāi)害,給人員帶來(lái)安全隱患,對(duì)煤礦生產(chǎn)造成經(jīng)濟(jì)損失。因此,采用物探方法進(jìn)行工作面前方異常體超前探測(cè)工作尤其重要,常用于巷道超前探測(cè)的物探技術(shù)主要有反射地震波法、瞬變電磁法、地質(zhì)雷達(dá)法等,而地震超前探測(cè)法中的槽波超前探測(cè)是利用僅在煤層中傳播的槽波對(duì)工作面前方地質(zhì)體進(jìn)行超前探測(cè),在探測(cè)斷層、陷落柱、采空區(qū)等地質(zhì)構(gòu)造方面有優(yōu)勢(shì)。由于煤礦的全空間工作環(huán)境與工作面處的強(qiáng)反射波和繞射波干擾,造成了地震波的復(fù)雜性,因此,利用數(shù)值模擬方法研究含巷道煤層中反射槽波的傳播規(guī)律,可以為實(shí)際巷道槽波超前探測(cè)提供一定的理論指導(dǎo)[1-5]。槽波數(shù)值模擬中通常將煤巖介質(zhì)假設(shè)為完全彈性介質(zhì),而實(shí)際煤層更接近于黏彈性介質(zhì)[6-7]。

    近年來(lái),諸多學(xué)者進(jìn)行了槽波超前探測(cè)數(shù)值模擬研究,如叢皖平等[8]提出將Rayleigh 波運(yùn)用于煤礦井下獨(dú)頭巷道超前探測(cè);楊思通等[9]進(jìn)行了Rayleigh 型槽波數(shù)值模擬并分析了工作面前方斷層的波場(chǎng)特征,提出Rayleigh 槽波能量在垂直分量能量比水平分量更強(qiáng);Yang Sitong 等[10]提出當(dāng)震源位于巷道前方煤層中,且檢波器沿著巷道中的煤壁布測(cè)時(shí),使用反射Rayleigh 槽波超前探測(cè)煤礦巷道地質(zhì)結(jié)構(gòu)。蔣錦朋等[11]提出震源為垂向集中力源激發(fā)時(shí),Rayleigh 槽波可應(yīng)用于槽波超前探測(cè);孫華超等[12]將Love 型槽波運(yùn)用于小斷層超前探測(cè)并分析其波場(chǎng)特征,并認(rèn)為x分量的Love 槽波信噪比好;同年呼邦兵等[13]提出利用x和z分量的Rayleigh 型槽波對(duì)采空區(qū)進(jìn)行超前探測(cè);之后Sun Huachao 等[14]提出將x分量反射Love 槽波應(yīng)用于陷落柱的超前探測(cè);梁紅波等[15]構(gòu)建了4 個(gè)斷層走向與巷道夾角不同的模型進(jìn)行槽波超前探測(cè)模擬,提出夾角小于40°時(shí)偏移處理后的角度變化不大。

    較多學(xué)者利用Kelvin-Voigt 模型研究黏彈性介質(zhì)的地震波傳播特征,例如苑春方等[16]研究了地震波在Kelvin-Voigt 均勻黏彈性介質(zhì)中的傳播速度和衰減系數(shù);劉瑞珣等[17]發(fā)現(xiàn)采用Kelvin-Voigt 模型比其他黏彈模型更加符合地下介質(zhì)是黏彈性質(zhì)的實(shí)際情況;楊思通等[18]研究了Kelvin-Voigt 黏彈介質(zhì)地震波衰減的影響因素,提出地震縱波振幅隨波阻抗降低和震源子波頻率增大衰減加快;嚴(yán)紅勇等[19]研究了二維非均勻各向同性Kelvin-Voigt 黏彈性介質(zhì)模擬,提出介質(zhì)的黏滯性對(duì)轉(zhuǎn)換波影響較大,且黏滯系數(shù)越大地震波能量衰減越大;姬廣忠等[20]進(jìn)行了Kelvin-Voigt 黏彈性TI 煤層介質(zhì)Love 槽波三維模擬,同年張壹等[21]對(duì)Kelvin-Voigt 等典型黏彈性模型介質(zhì)中的地震波衰減特征進(jìn)行了歸納分析;焦文杰等[22]進(jìn)行了Kelvin-Voigt黏彈性各向同性煤層介質(zhì)槽波三維數(shù)值模擬。

    目前黏彈介質(zhì)下槽波超前探測(cè)的理論研究很少,為更好地模擬實(shí)際煤礦工作的槽波超前探測(cè),筆者基于Kelvin-Voigt 黏彈性理論和黏彈三維一階速度-應(yīng)力方程,采用高階有限差分方法進(jìn)行獨(dú)頭巷道槽波超前探測(cè)的三維數(shù)值模擬,分析了煤層黏彈性及其變化、斷層的落差和斷層面與巷道的夾角等因素對(duì)反射槽波傳播的影響,為實(shí)際槽波超前探測(cè)提供理論支持。

    1 原理

    1.1 槽波超前探測(cè)原理

    槽波作為一種沿煤層傳播的導(dǎo)波,傳播距離遠(yuǎn)、能量強(qiáng)、具有明顯的頻散特性。根據(jù)不同的振動(dòng)形式,通常將槽波分為兩類,其一是Rayleigh 型槽波,由P 波與SV 波相互干涉形成,質(zhì)點(diǎn)在與煤層面相互垂直、與傳播方向平行的平面內(nèi)振動(dòng);其二是Love 型槽波,由SH 波相互干涉形成的,質(zhì)點(diǎn)在平行煤層面的平面內(nèi)垂直于波傳播方向上振動(dòng)[23-24]。槽波大部分能量不向圍巖輻射,攜帶大量地質(zhì)構(gòu)造信息,加之煤層的波導(dǎo)性,槽波被應(yīng)用于超前探測(cè)[2]。槽波超前探測(cè)法為槽波反射法的一種特殊情況,兩者原理類似。槽波超前探測(cè)原理為,激發(fā)震源后,槽波沿煤巷向工作面前方傳播,若在超前方向遇斷層面等反射界面,槽波會(huì)完全或不完全被阻斷,產(chǎn)生反射槽波,由布置在巷道側(cè)幫的檢波器接收反射槽波信號(hào),通過(guò)分析處理反射槽波記錄,實(shí)現(xiàn)對(duì)工作面前方構(gòu)造的探測(cè)[25-27],如圖1 所示。

    圖1 槽波超前探測(cè)原理Fig.1 Principle of in-seam wave advance detection

    1.2 Kelvin-Voigt 一階速度-應(yīng)力方程

    實(shí)際煤層具有黏彈特性,目前常用于地震正演模擬的黏彈性介質(zhì)模型主要有Kelvin-Voigt 模型、Maxwell 模型、標(biāo)準(zhǔn)線性模型等。其中Kelvin-Voigt模型由一個(gè)彈簧元件和一個(gè)阻尼器元件并聯(lián)構(gòu)成[17],如圖2 所示。

    圖2 Kelvin-Voigt 模型Fig.2 Kelvin-Voigt model

    Kelvin-Voigt 黏彈性介質(zhì)本構(gòu)關(guān)系為:

    式中:σ為應(yīng)力向量;σ1為 彈性單元體應(yīng)力向量;σ2為黏彈性單元體應(yīng)力向量;ε為應(yīng)變向量;C為彈性矩陣;η為阻尼器黏滯矩陣;t為時(shí)間,s;Q為品質(zhì)因子;ω為角頻率。

    用質(zhì)點(diǎn)振動(dòng)速度表示位移的一階導(dǎo)數(shù),來(lái)減少位移二階導(dǎo)數(shù)的使用,簡(jiǎn)化整個(gè)計(jì)算過(guò)程。在無(wú)外力影響或外力影響消失后,根據(jù)Kelvin-Voigt 模型,獲得三維黏彈一階速度-應(yīng)力方程[22]:

    式中:σxx、σyy、σzz為正應(yīng)力;τxy、τxz、τyz為剪應(yīng)力;vx、vy、vz為質(zhì)點(diǎn)振動(dòng)速度;x、y和z分別為三維空間中的平面水平方向、平面豎直方向和縱向。

    1.3 交錯(cuò)網(wǎng)格高階有限差分法

    應(yīng)用交錯(cuò)網(wǎng)格高階有限差分方法,將三維Kelvin-Voigt 一階速度-應(yīng)力方程進(jìn)行離散化。高階交錯(cuò)網(wǎng)格有限差分方法具有占用內(nèi)存小、通用性強(qiáng)、高精度和高穩(wěn)定性的特點(diǎn)。該方法是將一般的矩形差分網(wǎng)格替換為交錯(cuò)網(wǎng)格,借助交錯(cuò)網(wǎng)格的半程計(jì)算,來(lái)得到足夠高階的空間和時(shí)間精度的差分格式,以此來(lái)計(jì)算空間導(dǎo)數(shù)和時(shí)間導(dǎo)數(shù)。即在交錯(cuò)網(wǎng)格技術(shù)中,不僅要求空間網(wǎng)格的交錯(cuò),在時(shí)間網(wǎng)格上也要求交錯(cuò),變量的導(dǎo)數(shù)是在相應(yīng)的變量網(wǎng)格點(diǎn)間的中點(diǎn)上計(jì)算[28-29]。

    2 模型與參數(shù)

    模型在x、y、z方向的尺寸分別為300 m×100 m×25 m,模型網(wǎng)格dx、dy、dz分別為1 m×1 m×0.25 m,煤層厚度為5 m。巷道在x、y、z方向的大小分別為150 m×4 m×4 m,位于x方向10~160 m、y方向48~52 m、z方向10.5~14.5 m,工作面位于x=160 m 處,斷層位于x=240 m 處,如圖3 所示。煤層頂?shù)装鍘r性相同,巷道賦空值,模型各介質(zhì)參數(shù)見表1,邊界采用PML 吸收邊界。

    表1 模型介質(zhì)參數(shù)Table 1 Medium parameters of the model

    圖3 含斷層三維模型Fig.3 Schematic of the 3D model with faults

    圖4 為模型的觀測(cè)系統(tǒng),震源坐標(biāo)為(110,54,12.5),采用主頻150 Hz 的雷克子波爆炸震源,延遲時(shí)間20 ms。檢波器與震源同一水平線布置,過(guò)斷層布測(cè)于x方向80~280 m,當(dāng)落差大于1/2 倍煤厚時(shí),斷層面后的測(cè)線位于煤層頂板,道間距1 m,共201 個(gè)檢波點(diǎn)。數(shù)值模擬采樣間隔為0.05 ms。

    圖4 含斷層xoy 切面圖(z=12.5 m)Fig.4 xoy section with faults (z=12.5 m)

    3 黏彈煤層介質(zhì)對(duì)反射槽波的影響

    3.1 煤層介質(zhì)性質(zhì)對(duì)反射槽波的影響

    為研究煤層不同介質(zhì)對(duì)反射槽波的影響,設(shè)計(jì)了2 個(gè)三維煤系含斷層模型,分別是完全彈性介質(zhì)模型和黏彈性介質(zhì)模型(煤層QP=100,QS=50),落差d=2.5 m(1/2 煤厚),斷面與巷道垂直(夾角α=90°)。

    圖5a、圖5c 和圖5e 分別是黏彈介質(zhì)模型正演模擬合成地震記錄的x分量、y分量和z分量,圖5b、圖5d和圖5f 分別是完全彈性介質(zhì)模型正演模擬合成地震記錄的x分量、y分量和z分量,圖中橫軸為檢波器編號(hào),縱軸為傳播時(shí)間。結(jié)合模型與時(shí)距關(guān)系分析可知,激發(fā)震源產(chǎn)生直達(dá)波,直達(dá)波沿巷道向前傳播,在工作面處率先發(fā)生反射,形成工作面反射波,直達(dá)波繼續(xù)向前傳播遇斷層產(chǎn)生斷層反射波;直達(dá)波沿巷道向后傳播,在巷道頭處亦會(huì)發(fā)生反射,形成巷道頭反射波。

    圖5 三分量合成地震記錄Fig.5 Three components of synthetic seismograms

    整體上看圖5,各分量上相同波列在不同介質(zhì)模型地震記錄中的走時(shí)與分布相似,波振幅強(qiáng)度不同。圖5a 和圖5b 為分別是黏彈介質(zhì)和完全彈性介質(zhì)合成地震記錄的x分量,依據(jù)時(shí)距曲線將易識(shí)別的波列按照被接收時(shí)間先后排序:直達(dá)縱波、直達(dá)橫波、直達(dá)槽波、工作面反射波、反射橫波、反射槽波和巷道頭反射波。圖5a、圖5b 中,波組①由震源激發(fā)形成,出現(xiàn)時(shí)間最早,速度約為3 000 m/s,判斷為直達(dá)縱波。波組②與直達(dá)縱波同相,振幅能量較強(qiáng),速度約為1 700 m/s,判斷為直達(dá)橫波。波組③能量大,有頻散現(xiàn)象,速度約為900 m/s,判斷為直達(dá)槽波。波組④波列與直達(dá)波反相,出現(xiàn)較早,對(duì)其轉(zhuǎn)換波組不進(jìn)行分析,結(jié)合模型判斷為工作面反射波。波組⑥由直達(dá)橫波在斷面處轉(zhuǎn)換形成,速度約為1 700 m/s,判斷為斷層反射橫波,簡(jiǎn)稱為反射橫波。波組⑦能量大,同相軸斜率比反射橫波小,速度約為960 m/s,判斷為斷層反射槽波,為反射Rayleigh 槽波。在傳播100 ms 后,出現(xiàn)與直達(dá)波同相的系列干擾波,結(jié)合模型判斷其為巷道頭反射波,對(duì)其轉(zhuǎn)換波組不進(jìn)行分析,統(tǒng)一標(biāo)記為波組⑧。直達(dá)波遇斷層面后,斷面對(duì)不同地震波的阻斷作用不同,透射縱波與透射橫波能量衰減較小,槽波穿過(guò)斷層后能量明顯衰減;反射橫波能量大,有清晰地震記錄,而反射Rayleigh 槽波能量絕大部分衰減。完全彈性介質(zhì)中各波列傳播完整,斷層反射波能量衰減小,且黏彈介質(zhì)的反射橫波波列延續(xù)時(shí)間比完全彈性的短。

    圖5c 和圖5d 為分別是黏彈介質(zhì)和完全彈性介質(zhì)的y分量,依據(jù)時(shí)距曲線可以識(shí)別出:直達(dá)縱波、直達(dá)橫波、直達(dá)槽波、工作面反射波、反射橫波、反射槽波和巷道頭反射波。波組⑥由直達(dá)橫波在斷面處轉(zhuǎn)換形成,速度約為1 700 m/s,判斷為反射橫波。波組⑦波列呈現(xiàn)頻散特征,波振幅大,速度約為900 m/s,判斷為反射槽波,為反射Love 槽波。結(jié)合模型與時(shí)距關(guān)系,判斷工作面反射波附近的雙曲線干擾波為吸收邊界反射波。直達(dá)波遇斷層面后,透射橫波與透射槽波能量衰減較小,透射縱波能量明顯衰減;反射Love 槽波能量衰減較小,其他斷層反射波能量相對(duì)衰減較大。完全彈性介質(zhì)中各波列能量大于黏彈介質(zhì)。

    圖5e 和圖5f 為分別是黏彈介質(zhì)和完全彈性介質(zhì)的z分量,依據(jù)時(shí)距曲線可以識(shí)別出:直達(dá)縱波、直達(dá)橫波、直達(dá)槽波、工作面反射波、P-S 波、反射橫波、反射槽波和巷道頭反射波。波組⑤由直達(dá)縱波在斷面處轉(zhuǎn)換形成,波列較寬,速度約為1 800 m/s,判斷為直達(dá)縱波的轉(zhuǎn)換橫波,記為P-S 波。波組⑥與波組⑤相比斜率無(wú)明顯變化,由直達(dá)橫波在斷面處轉(zhuǎn)換形成,速度約為1 700 m/s,判斷為反射橫波。波組⑦波列呈現(xiàn)頻散特征,能量大,速度約為900 m/s,判斷為反射槽波,為反射Rayleigh 槽波。直達(dá)縱波遇斷層產(chǎn)生反射時(shí),另一部分能量轉(zhuǎn)換為反射縱波,可能由于縱波波速較快導(dǎo)致反射縱波與工作面反射波混雜,未見明顯同相軸。直達(dá)波遇斷層面后,透射縱波能量增大,橫波能量無(wú)明顯變化,槽波能量衰減較大;在斷面處產(chǎn)生的P-S 波、反射橫波和反射Rayleigh 槽波能量一定程度上都減小。黏彈介質(zhì)中反射Rayleigh 槽波能量明顯隨傳播距離增大而衰減,完全彈性介質(zhì)中無(wú)明顯能量變化。

    x分量上主要以反射橫波為主,反射Rayleigh 槽波能量較小;y分量上反射橫波波列范圍大,反射Love 槽波能量強(qiáng);z分量上P-S 波、反射橫波和反射Rayleigh 槽波波形特征皆明顯,信噪比好。與完全彈性介質(zhì)模型相比,黏彈介質(zhì)模型中地震波能量較小,地震波在黏彈介質(zhì)中傳播能量衰減速度比完全彈性介質(zhì)快,導(dǎo)致部分波列在傳播過(guò)程中衰減消失,沒有完全彈性介質(zhì)的波列完整。由此可知,黏彈介質(zhì)模型更加符合實(shí)際煤層中地震波傳播與衰減特征,更加貼近實(shí)際超前探測(cè)記錄。

    3.2 煤層品質(zhì)因子變化對(duì)反射槽波的影響

    為研究煤層不同品質(zhì)因子對(duì)反射槽波的影響,設(shè)計(jì)了3 個(gè)不同煤層Q值的三維黏彈性煤系含斷層模型,d=2.5 m,夾 角α=90°,QP、QS分別為40 與20、70 與35 和100 與50。

    取測(cè)線第1 道至第81 道進(jìn)行后續(xù)含斷層模型的反射槽波對(duì)比分析。圖6a、圖6b 和圖6c 分別為QP=40、QS=20 與QP=70、QS=35 及QP=100、QS=50 模型的x分量。圖7、圖8 分別為不同Q值模型正演模擬的y和z分量。根據(jù)時(shí)距關(guān)系3 個(gè)分量上皆可分辨出工作面反射波、反射橫波、反射槽波和巷道頭反射波,以及z分量的P-S 波。各分量上相同波列在不同Q值模型地震記錄中的波形相似,波振幅強(qiáng)度不同。x分量上,反射Rayleigh 槽波和巷道頭反射波能量隨Q值增大明顯增大,工作面反射波與反射橫波能量隨Q值增大而增大。當(dāng)QP=40、QS=20 時(shí),反射Rayleigh 槽波無(wú)明顯同相軸。y分量上,各波列能量受Q值影響與x分量類似。隨Q值增大,反射Love 槽波、工作面反射波、反射橫波和巷道頭反射波能量增大。當(dāng)QP=40、QS=20 時(shí),反射Love 槽波與反射橫波能量小,波形模糊。z分量上,各波列能量受Q值影響與x分量類似。隨Q值增大,反射Rayleigh 槽波和巷道頭反射波能量明顯增大;工作面反射波、反射橫波能量增大;P-S 波能量變化較小。當(dāng)QP=40、QS=20 時(shí),反射Rayleigh 槽波與反射橫波能量小。

    圖6 不同品質(zhì)因子模型地震記錄(x 分量,測(cè)線第1 道到第81 道)Fig.6 Seismograms of models with different quality factors (x component,the 1st to the 81th traces of the survey lines)

    圖7 不同品質(zhì)因子模型地震記錄(y 分量)Fig.7 Seismograms of models with different quality factors (y component)

    圖8 不同品質(zhì)因子模型地震記錄(z 分量)Fig.8 Seismograms of models with different quality factors (z component)

    各分量上煤層Q值變化對(duì)不同地震波能量衰減影響不同,各波列的能量衰減都隨Q值增大而減小。隨Q值增大不同分量的工作面反射波和巷道頭反射波能量增大;反射槽波在三分量上隨Q值增大能量都明顯增大;z分量的P-S 波受Q值影響??;反射橫波在x分量受Q值影響小,在z分量受影響較大。QP=40、QS=20 時(shí)三分量上的反射槽波能量皆衰減大,傳播時(shí)間越久能量越小,y分量上反射Love 槽波與反射橫波能量都較小,因此,可利用x分量反射橫波與z分量P-S 波進(jìn)行超前探測(cè)。

    4 斷層變化對(duì)反射槽波的影響

    4.1 斷層落差變化對(duì)反射槽波的影響

    為研究斷層落差對(duì)反射槽波的影響,設(shè)計(jì)了3 個(gè)三維黏彈性煤系含斷層模型,QP為100,QS為50,夾角α=90°,斷層落差d分別為3.75、5 和7.5 m。

    圖9 為斷層落差3.75、5 和7.5 m 模型的x分量,圖10、圖11 分別為不同落差模型正演模擬的y和z分量,根據(jù)時(shí)距關(guān)系3 個(gè)分量上皆可分辨出工作面反射波、反射橫波、反射槽波和巷道頭反射波,以及z分量的P-S 波。各分量上相同波列在不同落差模型地震記錄中的走時(shí)相似,波振幅強(qiáng)度不同。x分量上,工作面反射波和巷道頭反射波同相軸振幅隨落差變化不明顯,能量變化較小。隨落差增大,斷面對(duì)波的阻斷作用增大,反射橫波和反射Rayleigh 槽波能量增大。y分量上,各波列能量受落差影響與x分量類似。隨落差增大,反射橫波和反射Love 槽波能量增大。z分量上,各波列能量受落差影響與x、y分量不同。斷層落差由2.5 m 增大至3.75 m 時(shí),地震記錄中的P-S 波和反射Rayleigh 槽波波形愈加明顯,能量逐漸增大;落差由3.75 m 增大至7.5 m 時(shí),直達(dá)波遇斷面后轉(zhuǎn)換波組能量迅速衰減,P-S 波和反射Rayleigh 槽波能量減小。落差2.5 m 增大至5 m,反射橫波能量增大;落差5 m 增大至7.5 m,反射橫波振幅變化較小,波列延續(xù)時(shí)間變短。

    圖9 不同落差模型地震記錄(x 分量)Fig.9 Seismograms of models with different fault throws (x component)

    根據(jù)上述分析可知,不同分量上地震波受落差變化的影響不同。工作面反射波和巷道頭反射波在三分量上隨落差增大能量變化都較小,說(shuō)明斷層落差對(duì)其影響??;在x和y分量上,斷層落差越大斷面對(duì)波的阻斷作用越大,反射橫波和反射槽波能量隨落差增大而增大;z分量上,落差小于3/4 煤厚時(shí),反射槽波、P-S 波和反射橫波隨落差增大能量增大;落差超過(guò)煤厚時(shí),地震波的轉(zhuǎn)換作用急劇減弱,由直達(dá)波轉(zhuǎn)換形成的P-S 波和反射Rayleigh 槽波能量衰減大,反射橫波能量變化較小。在斷層上下盤完全錯(cuò)開后,z分量上反射Rayleigh 槽波能量小,可利用y分量的反射Love 槽波以及x和z分量的反射橫波進(jìn)行超前探測(cè)。

    4.2 斷層夾角變化對(duì)反射槽波的影響

    為研究不同斷層夾角對(duì)反射槽波的影響,設(shè)計(jì)了3 個(gè)三維黏彈性煤系含斷層模型,QP均為100,QS均為50,d=2.5 m,夾角α分別為75°、60°和45°。

    圖12 為夾角75°、60°和45°模型的x分量,圖13、圖14 分別為不同夾角模型正演模擬的y和z分量,根據(jù)時(shí)距關(guān)系3 個(gè)分量上皆可分辨出工作面反射波、反射橫波、反射槽波和巷道頭反射波,以及z分量的P-S 波。各分量上相同波列在不同夾角模型地震記錄中的走時(shí)不同,波振幅強(qiáng)度不同。x分量上,工作面反射波和巷道頭反射波同相軸無(wú)明顯變化,走時(shí)與能量變化較??;反射橫波和反射Rayleigh 槽波同相軸斜率增大,波列間距減小,反射Rayleigh 槽波能量增大。y分量上,工作面反射波和巷道頭反射波同相軸無(wú)明顯變化,反射橫波與反射Love 槽波被接收時(shí)間提前。斷層夾角由90°減小至75°時(shí),反射橫波與反射Love 槽波能量增大,同相軸斜率增大;斷層夾角由75°減小至45°時(shí),反射橫波與反射Love 槽波能量減小,波列間距明顯減小。z分量上,各波列受夾角影響與y分量類似。P-S 波、反射橫波和反射Rayleigh 槽波被接收到的時(shí)間提前,同相軸斜率隨之增大,波能量呈現(xiàn)先增后減趨勢(shì)。

    圖13 不同夾角模型地震記錄(y 分量)Fig.13 Seismograms of models with different angles between a fault plane and a tunnel (y component)

    圖14 不同夾角模型地震記錄(z 分量)Fig.14 Seismograms of models with different angles between a fault plane and a tunnel (z component)

    圖15 為夾角90°、75°、60°和45°模型的z分量波場(chǎng)快照在z=12.5 m 的切片圖,結(jié)合模型與各波列速度可知,傳播時(shí)間200 ms 切面上波列主要為巷道頭干擾波、P-S 波、反射橫波、反射Rayleigh 槽波、工作面繞射波和透射波。其中波組⑤的寬度最大,傳播速度較快,判斷為P-S 波;波組⑦波形最清晰,傳播速度較慢,判斷為反射Rayleigh 槽波;波組⑥傳播速度介于波組⑤和波組⑦之間,波形較清晰,判斷為反射橫波;波組⑧與其他波列傳播方向相反,產(chǎn)生于巷道頭,判斷為巷道頭干擾波;波組⑨產(chǎn)生于工作面,波形混雜于斷層反射波之中,判斷為工作面繞射波。隨夾角減小巷道頭反射波波形無(wú)明顯變化;P-S 波、反射橫波和反射Rayleigh 槽波的反射角度隨夾角減小而增大,測(cè)線位置的各波列清晰度在夾角小于75°后明顯降低;工作面繞射波能量小于斷層反射波,對(duì)波場(chǎng)影響較小,地震記錄上難識(shí)別,對(duì)超前探測(cè)影響亦小。

    圖15 不同夾角模型200 ms 波場(chǎng)快照Fig.15 200 ms wave field snapshots of models with different angles between a fault plane and a tunnel

    結(jié)合圖14 地震記錄與圖15 波場(chǎng)快照分析可知,夾角角度越小,直達(dá)波離斷層的相對(duì)距離越小,斷層反射波產(chǎn)生的時(shí)間越早,反射角度越大。夾角由90°減小至75°時(shí),檢波器與斷層面間距較大,反射槽波有足夠的發(fā)育空間,各波列間隔較大,波形特征明顯,且近似垂直反射,反射Rayleigh 槽波能量大;夾角由75°減小至45°時(shí),測(cè)線與斷層間距變小,各斷層反射波出現(xiàn)時(shí)間明顯提前,反射Rayleigh 槽波的反射角度增大,能量隨之減小。P-S 波與反射橫波同理。隨夾角減小各斷層反射波波列間隔減小,各波列難分離。

    根據(jù)上述分析可知,不同分量上夾角對(duì)不同地震波的影響不同,且通過(guò)反射槽波的同相軸斜率可初步判斷斷層夾角大小。工作面反射波和巷道頭反射波能量在3 個(gè)分量上變化較小,說(shuō)明受斷層夾角影響較??;隨夾角角度減小,P-S 波、反射橫波和反射槽波同相軸斜率增大。在x分量上,隨斷層夾角角度減小,反射橫波與反射槽波被接收時(shí)間提前,能量增大;在y和z分量上,P-S 波、反射橫波和反射槽波在夾角60°至90°模型中能量大,波形特征明顯、易識(shí)別,槽波超前探測(cè)效果好;夾角小于60°時(shí),反射槽波與其他波列難分離,探測(cè)分辨率隨之降低,可利用x分量的反射Rayleigh 槽波進(jìn)行超前探測(cè)。

    5 結(jié)論

    a.槽波在黏彈介質(zhì)中能量衰減大,各分量上的反射槽波能量遠(yuǎn)小于直達(dá)槽波能量,z分量上反射Rayleigh 槽波隨傳播距離增大能量衰減速度明顯大于完全彈性介質(zhì)。說(shuō)明黏彈介質(zhì)中槽波的傳播更加符合實(shí)際煤層槽波傳播過(guò)程中能量的衰減特征。

    b.在x、y和z分量上接收到的信號(hào)存在差異,x分量上主要以反射橫波為主,反射Rayleigh 槽波能量較弱,不易識(shí)別;y分量上反射Love 槽波信噪比好;z分量上波列連續(xù)性較好,且P-S 波、反射橫波和反射Rayleigh 槽波波形特征明顯。

    c.煤層Q值大小對(duì)反射槽波傳播有較大的影響。隨煤層Q值增大反射槽波能量增大,超前探測(cè)效果越好。QP=40、QS=20 時(shí)3 個(gè)分量的反射槽波能量皆小,同相軸都不明顯,可利用x分量的反射橫波與z分量的P-S 波進(jìn)行超前探測(cè)。

    d.斷層落差與夾角在三分量上對(duì)反射槽波影響不同。斷層落差小于煤厚時(shí),y和z分量上反射槽波能量大,槽波超前探測(cè)效果都較好;落差大于煤厚時(shí),z分量的反射Rayleigh 槽波能量小,可利用y分量的Love 槽波以及x和z分量的反射橫波進(jìn)行超前探測(cè)。斷層夾角大于60°時(shí),y和z分量上反射槽波能量大,槽波超前探測(cè)效果都較好;夾角小于60°時(shí),x和y分量上反射槽波能量大,槽波超前探測(cè)都較好。

    e.很多實(shí)際煤層Q值條件與QP=40、QS=20 斷層模型相似,三分量上反射槽波能量都較小,槽波超前探測(cè)效果差,因此,可根據(jù)其他Q值斷層模型數(shù)值模擬的分析結(jié)果,綜合應(yīng)用三分量上的反射槽波、反射橫波和P-S 波,相互驗(yàn)證補(bǔ)充,提高煤巷槽波超前探測(cè)的精確度和可靠性。

    猜你喜歡
    波組橫波夾角
    橫波技術(shù)在工程物探中的應(yīng)用分析
    心電圖碎裂QRS波對(duì)于急性心肌梗死診斷的臨床應(yīng)用價(jià)值
    探究鐘表上的夾角
    求解異面直線夾角問題的兩個(gè)路徑
    心電圖異常J波評(píng)估急性ST段抬高型心梗患者短期預(yù)后的價(jià)值分析
    心電圖碎裂QRS波與急性心肌梗死患者室性心律失常及左心室收縮功能的相關(guān)性
    血糖波動(dòng)對(duì)2 型糖尿病合并肺癌患者化療效果的影響研究
    任意夾角交叉封閉邊界內(nèi)平面流線計(jì)算及應(yīng)用
    直線轉(zhuǎn)角塔L形絕緣子串夾角取值分析
    廣西電力(2016年5期)2016-07-10 09:16:44
    揚(yáng)眉一顧,妖嬈橫波處
    欧美bdsm另类| 国产亚洲精品av在线| 2021少妇久久久久久久久久久| 亚洲精品中文字幕在线视频 | 亚洲av福利一区| 午夜福利高清视频| 亚洲不卡免费看| 亚洲国产最新在线播放| 日本猛色少妇xxxxx猛交久久| 亚洲精品日韩av片在线观看| 精品熟女少妇av免费看| 麻豆成人av视频| 又黄又爽又刺激的免费视频.| 一区二区三区四区激情视频| 国产成人91sexporn| 欧美三级亚洲精品| 欧美高清成人免费视频www| 一级黄片播放器| av卡一久久| 国产欧美另类精品又又久久亚洲欧美| 亚洲av国产av综合av卡| 国产白丝娇喘喷水9色精品| 十八禁国产超污无遮挡网站| 少妇人妻精品综合一区二区| 国精品久久久久久国模美| av免费在线看不卡| 免费高清在线观看视频在线观看| 成年女人在线观看亚洲视频 | 能在线免费看毛片的网站| 菩萨蛮人人尽说江南好唐韦庄| 国内精品美女久久久久久| 日韩国内少妇激情av| 日韩在线高清观看一区二区三区| 成年女人在线观看亚洲视频 | 久久久久九九精品影院| 亚洲最大成人av| 久久久久久久久中文| 亚洲精品日本国产第一区| 激情 狠狠 欧美| 亚洲精品日韩在线中文字幕| 久久久久免费精品人妻一区二区| 成年av动漫网址| 成人亚洲精品av一区二区| 99热6这里只有精品| 青青草视频在线视频观看| 国产av不卡久久| 成人国产麻豆网| 午夜福利高清视频| 成人毛片60女人毛片免费| 国产老妇女一区| 久久久欧美国产精品| 成人一区二区视频在线观看| 久久人人爽人人片av| 日韩欧美精品免费久久| 99九九线精品视频在线观看视频| 久久人人爽人人爽人人片va| 日韩欧美精品v在线| h日本视频在线播放| 精品久久久久久久人妻蜜臀av| 日韩一区二区视频免费看| 成人亚洲精品av一区二区| 亚洲四区av| 免费观看无遮挡的男女| 欧美另类一区| 国产欧美另类精品又又久久亚洲欧美| 在线观看一区二区三区| 欧美日韩在线观看h| 22中文网久久字幕| 国产精品.久久久| 精品少妇黑人巨大在线播放| 久久久久久久久久久丰满| 最近的中文字幕免费完整| 性插视频无遮挡在线免费观看| 嫩草影院精品99| 日日撸夜夜添| 婷婷色综合大香蕉| 在线观看av片永久免费下载| 欧美区成人在线视频| 亚洲,欧美,日韩| 精品国内亚洲2022精品成人| 亚洲久久久久久中文字幕| 欧美xxxx性猛交bbbb| 97热精品久久久久久| 免费观看性生交大片5| 亚洲欧洲日产国产| 国产 亚洲一区二区三区 | 欧美不卡视频在线免费观看| kizo精华| 亚洲无线观看免费| 久久久久久久久久成人| 直男gayav资源| 麻豆精品久久久久久蜜桃| 久久久久网色| 人妻少妇偷人精品九色| 99久久精品国产国产毛片| 最后的刺客免费高清国语| 最近手机中文字幕大全| 免费观看的影片在线观看| 91在线精品国自产拍蜜月| 亚洲av成人av| 精华霜和精华液先用哪个| 大话2 男鬼变身卡| 男的添女的下面高潮视频| 国产三级在线视频| 国产免费福利视频在线观看| 国产毛片a区久久久久| 免费看光身美女| 国产精品人妻久久久久久| 80岁老熟妇乱子伦牲交| 欧美人与善性xxx| 国产一区二区三区综合在线观看 | 乱码一卡2卡4卡精品| 久久精品国产亚洲av天美| 高清在线视频一区二区三区| 极品教师在线视频| 亚洲精品国产av成人精品| 亚洲国产精品成人综合色| 国产黄片视频在线免费观看| 午夜福利在线观看吧| 嫩草影院精品99| 欧美性猛交╳xxx乱大交人| 一级毛片电影观看| 男女边吃奶边做爰视频| 97在线视频观看| 欧美另类一区| 久久久成人免费电影| 美女主播在线视频| 亚洲精品456在线播放app| 婷婷色av中文字幕| 中文字幕免费在线视频6| 亚洲不卡免费看| 老师上课跳d突然被开到最大视频| 亚洲av电影不卡..在线观看| 乱系列少妇在线播放| 亚洲欧美中文字幕日韩二区| 高清欧美精品videossex| av国产久精品久网站免费入址| 天天躁日日操中文字幕| 一夜夜www| 国产av国产精品国产| 偷拍熟女少妇极品色| 亚洲aⅴ乱码一区二区在线播放| 一级毛片我不卡| 在线 av 中文字幕| 简卡轻食公司| 午夜福利视频1000在线观看| 2021少妇久久久久久久久久久| 好男人在线观看高清免费视频| 午夜激情福利司机影院| 天堂俺去俺来也www色官网 | 国产v大片淫在线免费观看| 熟妇人妻不卡中文字幕| 亚洲久久久久久中文字幕| 日本三级黄在线观看| 蜜桃亚洲精品一区二区三区| 18禁动态无遮挡网站| 亚洲精品一区蜜桃| 国产伦精品一区二区三区四那| 久久亚洲国产成人精品v| 亚洲av中文字字幕乱码综合| av在线亚洲专区| 成年版毛片免费区| 欧美日韩综合久久久久久| 亚洲精品aⅴ在线观看| 深夜a级毛片| 少妇丰满av| 日本wwww免费看| 亚洲一级一片aⅴ在线观看| 国产精品日韩av在线免费观看| 亚洲av电影在线观看一区二区三区 | 久久久久久久久久久丰满| 成年av动漫网址| 国产亚洲一区二区精品| 能在线免费观看的黄片| 国产一级毛片在线| 午夜精品国产一区二区电影 | 我要看日韩黄色一级片| 校园人妻丝袜中文字幕| 六月丁香七月| 亚洲av免费高清在线观看| 亚洲人与动物交配视频| 国产v大片淫在线免费观看| 99热这里只有精品一区| 又黄又爽又刺激的免费视频.| 黄色欧美视频在线观看| 午夜免费观看性视频| 中文字幕av在线有码专区| 91精品一卡2卡3卡4卡| 国产亚洲精品av在线| 日本熟妇午夜| 免费av不卡在线播放| 国产精品一区二区在线观看99 | 日韩欧美精品v在线| 久久久久久久久大av| 高清日韩中文字幕在线| 夫妻午夜视频| 一夜夜www| 麻豆成人午夜福利视频| 麻豆成人午夜福利视频| 最近2019中文字幕mv第一页| 国产精品美女特级片免费视频播放器| 女人被狂操c到高潮| 国产高清三级在线| a级一级毛片免费在线观看| 国产成人freesex在线| 国产免费福利视频在线观看| 久久国产乱子免费精品| 国产精品久久久久久精品电影| 国产成人福利小说| 国产精品女同一区二区软件| 嫩草影院入口| 成人性生交大片免费视频hd| 亚洲不卡免费看| 18+在线观看网站| 国产精品一区二区三区四区免费观看| 亚洲色图av天堂| 成人毛片60女人毛片免费| 三级国产精品欧美在线观看| 国产真实伦视频高清在线观看| 97超视频在线观看视频| 最近中文字幕高清免费大全6| 五月天丁香电影| 又爽又黄无遮挡网站| 十八禁国产超污无遮挡网站| 赤兔流量卡办理| 日本免费在线观看一区| 欧美三级亚洲精品| 国产黄频视频在线观看| 亚洲精品自拍成人| 国产真实伦视频高清在线观看| 天天躁夜夜躁狠狠久久av| 亚洲成人精品中文字幕电影| 肉色欧美久久久久久久蜜桃 | 尾随美女入室| 麻豆久久精品国产亚洲av| 国产精品无大码| 一级毛片 在线播放| 欧美一区二区亚洲| 性色avwww在线观看| 国产精品国产三级国产专区5o| 日韩欧美国产在线观看| 久久久久免费精品人妻一区二区| 一级毛片黄色毛片免费观看视频| 日本与韩国留学比较| 最近最新中文字幕免费大全7| 欧美一区二区亚洲| 国产淫片久久久久久久久| 一级av片app| 久久久久久久亚洲中文字幕| 色5月婷婷丁香| 最近中文字幕2019免费版| 亚洲精品自拍成人| 日本wwww免费看| 亚洲精品国产成人久久av| 欧美成人a在线观看| 简卡轻食公司| 亚洲精品国产av蜜桃| av专区在线播放| 国产精品一区二区三区四区免费观看| 国产永久视频网站| 一级毛片 在线播放| 日本av手机在线免费观看| 精品一区二区三区视频在线| 日日啪夜夜爽| 十八禁网站网址无遮挡 | 伊人久久精品亚洲午夜| 麻豆久久精品国产亚洲av| 美女国产视频在线观看| 亚洲av电影在线观看一区二区三区 | 亚洲自拍偷在线| a级毛片免费高清观看在线播放| 18禁在线无遮挡免费观看视频| 免费观看的影片在线观看| 三级国产精品片| 日韩视频在线欧美| 欧美激情久久久久久爽电影| 国产精品福利在线免费观看| 一级毛片电影观看| 久久久久网色| 成年av动漫网址| 熟妇人妻久久中文字幕3abv| 亚洲精品aⅴ在线观看| 亚洲国产精品专区欧美| 在线观看美女被高潮喷水网站| 精品人妻熟女av久视频| 中文字幕久久专区| 成人亚洲欧美一区二区av| 国产极品天堂在线| 在线免费观看的www视频| 日韩av在线大香蕉| 干丝袜人妻中文字幕| 欧美高清成人免费视频www| 国产免费一级a男人的天堂| 99久久中文字幕三级久久日本| 欧美激情在线99| 99久国产av精品| 国产色婷婷99| 精品不卡国产一区二区三区| 五月伊人婷婷丁香| 两个人视频免费观看高清| 久久精品夜色国产| 免费少妇av软件| 精品不卡国产一区二区三区| 免费观看a级毛片全部| 国产一区二区亚洲精品在线观看| 日本一二三区视频观看| 日本三级黄在线观看| 国产成人一区二区在线| 日韩一本色道免费dvd| 国产精品人妻久久久影院| 97热精品久久久久久| 中文字幕亚洲精品专区| 特级一级黄色大片| av在线观看视频网站免费| 男女国产视频网站| 国产老妇伦熟女老妇高清| av黄色大香蕉| 白带黄色成豆腐渣| 日日啪夜夜撸| 午夜激情久久久久久久| 18禁裸乳无遮挡免费网站照片| 99久久人妻综合| 大香蕉久久网| 欧美bdsm另类| 69av精品久久久久久| 特大巨黑吊av在线直播| 97超碰精品成人国产| 日韩一本色道免费dvd| 高清av免费在线| 97热精品久久久久久| 又爽又黄a免费视频| 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品专区欧美| 日韩成人av中文字幕在线观看| 人妻系列 视频| 亚洲无线观看免费| 欧美日韩在线观看h| 午夜亚洲福利在线播放| 久久久久性生活片| 中文资源天堂在线| 亚洲成人精品中文字幕电影| 舔av片在线| 亚洲欧美一区二区三区国产| 夜夜爽夜夜爽视频| 久久99热这里只有精品18| 一个人观看的视频www高清免费观看| 欧美高清性xxxxhd video| .国产精品久久| 丰满人妻一区二区三区视频av| 老司机影院成人| 成人综合一区亚洲| 日韩欧美 国产精品| 国产精品日韩av在线免费观看| 免费观看性生交大片5| 国产淫片久久久久久久久| 精品亚洲乱码少妇综合久久| 看黄色毛片网站| 国产三级在线视频| 一夜夜www| 尤物成人国产欧美一区二区三区| 国产有黄有色有爽视频| 国产精品人妻久久久影院| 久99久视频精品免费| 中文字幕av在线有码专区| 蜜臀久久99精品久久宅男| 18禁在线播放成人免费| 亚洲国产成人一精品久久久| 久久鲁丝午夜福利片| 美女主播在线视频| 国产成人a∨麻豆精品| av在线老鸭窝| 精品一区二区三区人妻视频| 久久精品综合一区二区三区| 国产69精品久久久久777片| 91狼人影院| 婷婷色综合大香蕉| 91精品伊人久久大香线蕉| 国产乱人偷精品视频| 插逼视频在线观看| 亚洲自偷自拍三级| 国产黄色视频一区二区在线观看| 日韩在线高清观看一区二区三区| 自拍偷自拍亚洲精品老妇| 日韩成人av中文字幕在线观看| 99热这里只有精品一区| 黄色日韩在线| 高清日韩中文字幕在线| 午夜福利视频1000在线观看| 丝袜喷水一区| 国产免费一级a男人的天堂| 精品久久国产蜜桃| 韩国av在线不卡| 高清视频免费观看一区二区 | 精品久久久精品久久久| 国产精品麻豆人妻色哟哟久久 | 成人无遮挡网站| 国产成人精品一,二区| 22中文网久久字幕| 亚洲av电影在线观看一区二区三区 | 亚洲精品色激情综合| 精品一区二区三区视频在线| 午夜视频国产福利| 中文欧美无线码| 午夜福利在线观看吧| 欧美一级a爱片免费观看看| 日韩精品青青久久久久久| 国产免费视频播放在线视频 | 日韩制服骚丝袜av| 亚洲精品日本国产第一区| 18+在线观看网站| 黄片wwwwww| 日本欧美国产在线视频| 亚洲国产日韩欧美精品在线观看| 少妇被粗大猛烈的视频| 亚洲最大成人手机在线| 国产精品麻豆人妻色哟哟久久 | 免费观看无遮挡的男女| 精品久久久久久久久久久久久| av网站免费在线观看视频 | 青春草国产在线视频| 亚洲经典国产精华液单| 中文精品一卡2卡3卡4更新| 最近的中文字幕免费完整| 免费大片黄手机在线观看| 99久久中文字幕三级久久日本| 2021少妇久久久久久久久久久| av在线亚洲专区| 亚洲国产精品专区欧美| 国产精品美女特级片免费视频播放器| 欧美激情国产日韩精品一区| 亚洲国产av新网站| 97精品久久久久久久久久精品| 久久鲁丝午夜福利片| av又黄又爽大尺度在线免费看| 久久精品国产亚洲网站| 大陆偷拍与自拍| 国产一区二区亚洲精品在线观看| 极品教师在线视频| 亚洲av日韩在线播放| 日韩成人av中文字幕在线观看| 在线观看人妻少妇| 亚洲精品日韩在线中文字幕| 国产免费视频播放在线视频 | 婷婷色综合大香蕉| 国精品久久久久久国模美| 一级爰片在线观看| 亚洲精品成人久久久久久| 亚洲不卡免费看| 人妻制服诱惑在线中文字幕| 免费电影在线观看免费观看| 非洲黑人性xxxx精品又粗又长| 成人高潮视频无遮挡免费网站| 国产精品国产三级国产av玫瑰| 乱系列少妇在线播放| 日本熟妇午夜| 国产欧美另类精品又又久久亚洲欧美| 午夜激情福利司机影院| 女人久久www免费人成看片| 精品午夜福利在线看| 日本黄色片子视频| 色网站视频免费| 日韩强制内射视频| 成年av动漫网址| 男的添女的下面高潮视频| 啦啦啦中文免费视频观看日本| 真实男女啪啪啪动态图| 国产精品国产三级专区第一集| 六月丁香七月| 特级一级黄色大片| 色哟哟·www| 亚洲av在线观看美女高潮| 精品久久久久久久久亚洲| 中文乱码字字幕精品一区二区三区 | 老师上课跳d突然被开到最大视频| 91午夜精品亚洲一区二区三区| 毛片一级片免费看久久久久| 一级黄片播放器| 女人久久www免费人成看片| 精品一区在线观看国产| 乱码一卡2卡4卡精品| 免费观看a级毛片全部| 亚洲国产色片| 免费无遮挡裸体视频| 狂野欧美激情性xxxx在线观看| 国产综合懂色| 日韩中字成人| 国产中年淑女户外野战色| 夜夜看夜夜爽夜夜摸| 日本爱情动作片www.在线观看| 精品国内亚洲2022精品成人| 久久精品夜夜夜夜夜久久蜜豆| 午夜激情福利司机影院| 国产av码专区亚洲av| 美女xxoo啪啪120秒动态图| 国产精品麻豆人妻色哟哟久久 | 成人av在线播放网站| .国产精品久久| 成人漫画全彩无遮挡| 午夜精品一区二区三区免费看| 久久久久精品久久久久真实原创| 99热全是精品| 国产av码专区亚洲av| 久久97久久精品| 赤兔流量卡办理| 精品一区二区免费观看| 免费看av在线观看网站| 中文字幕av成人在线电影| 18禁裸乳无遮挡免费网站照片| 一个人看视频在线观看www免费| 国产黄a三级三级三级人| 伊人久久国产一区二区| 只有这里有精品99| 成人性生交大片免费视频hd| 黄色日韩在线| 国产黄片视频在线免费观看| 国产又色又爽无遮挡免| 久久久久久久久久黄片| 久久久久久久国产电影| 国产精品福利在线免费观看| 能在线免费观看的黄片| 全区人妻精品视频| 性插视频无遮挡在线免费观看| 日韩成人伦理影院| 亚洲内射少妇av| 又爽又黄a免费视频| 看非洲黑人一级黄片| 啦啦啦啦在线视频资源| 日韩人妻高清精品专区| 午夜精品在线福利| 免费av不卡在线播放| 一个人免费在线观看电影| kizo精华| 一级a做视频免费观看| 在线观看人妻少妇| 国产伦精品一区二区三区视频9| 国产一级毛片在线| 亚洲真实伦在线观看| 成人av在线播放网站| 亚洲欧美一区二区三区黑人 | 国产色婷婷99| 97人妻精品一区二区三区麻豆| 亚洲人与动物交配视频| 亚洲欧美精品自产自拍| videossex国产| 国产毛片a区久久久久| 久久综合国产亚洲精品| 国产伦精品一区二区三区四那| 日韩欧美三级三区| 国产在线一区二区三区精| 美女主播在线视频| 久久精品久久久久久久性| 91午夜精品亚洲一区二区三区| 99久国产av精品国产电影| 久久久久久久久久人人人人人人| 亚洲精品久久午夜乱码| 日日摸夜夜添夜夜添av毛片| 男女边吃奶边做爰视频| 精品久久久噜噜| 日韩欧美三级三区| 日日啪夜夜爽| 亚洲美女视频黄频| 亚洲精品国产av蜜桃| 亚洲成人中文字幕在线播放| 国产精品国产三级国产专区5o| 男女啪啪激烈高潮av片| 在线免费十八禁| 两个人的视频大全免费| 久久久色成人| 国产免费一级a男人的天堂| 日本一二三区视频观看| 欧美日韩综合久久久久久| 欧美成人a在线观看| 女人久久www免费人成看片| 两个人视频免费观看高清| 国产探花极品一区二区| 国产黄色小视频在线观看| 国产伦在线观看视频一区| 青青草视频在线视频观看| 麻豆成人av视频| 天美传媒精品一区二区| 天天一区二区日本电影三级| 校园人妻丝袜中文字幕| 国产探花极品一区二区| 日本av手机在线免费观看| 内射极品少妇av片p| 69人妻影院| 女人久久www免费人成看片| 日日摸夜夜添夜夜爱| 伦理电影大哥的女人| 国精品久久久久久国模美| 不卡视频在线观看欧美| 一级毛片久久久久久久久女| 亚洲国产日韩欧美精品在线观看| or卡值多少钱| 色5月婷婷丁香| 69人妻影院| 白带黄色成豆腐渣| 欧美xxxx黑人xx丫x性爽| 人体艺术视频欧美日本| 男插女下体视频免费在线播放| 熟女电影av网| 老司机影院成人| 国产乱人视频| 国产69精品久久久久777片| 欧美日韩精品成人综合77777| 欧美高清性xxxxhd video| 精品一区二区免费观看| 国产大屁股一区二区在线视频| 久久久久久久久久黄片| 久久精品熟女亚洲av麻豆精品 | 午夜福利在线观看吧| av.在线天堂| 亚洲欧美中文字幕日韩二区| 久久久亚洲精品成人影院| 99热6这里只有精品|