• <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)眉一顧,妖嬈橫波處
    www.熟女人妻精品国产| videos熟女内射| 欧美在线黄色| 精品国产露脸久久av麻豆| 一级片免费观看大全| 精品久久久久久电影网| 午夜免费鲁丝| 国产成人欧美| 一本一本久久a久久精品综合妖精| 夫妻性生交免费视频一级片| 一区二区av电影网| 久久影院123| 欧美国产精品va在线观看不卡| 三上悠亚av全集在线观看| 久久天堂一区二区三区四区| 美女大奶头黄色视频| 成人国语在线视频| 丁香六月欧美| 成年av动漫网址| 人妻一区二区av| 电影成人av| 少妇的丰满在线观看| 中文字幕高清在线视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲av日韩精品久久久久久密 | 亚洲久久久国产精品| 中文天堂在线官网| 久久久久国产一级毛片高清牌| 91精品伊人久久大香线蕉| 妹子高潮喷水视频| 欧美黄色片欧美黄色片| 十八禁人妻一区二区| 欧美另类一区| 韩国av在线不卡| 亚洲精品久久成人aⅴ小说| 精品国产乱码久久久久久男人| 国产国语露脸激情在线看| 天堂8中文在线网| 久热爱精品视频在线9| 精品视频人人做人人爽| 国产成人精品无人区| 18禁裸乳无遮挡动漫免费视频| 亚洲精品av麻豆狂野| 久久精品久久久久久久性| 国产精品久久久人人做人人爽| 夫妻性生交免费视频一级片| 日日撸夜夜添| 中文字幕高清在线视频| 美女国产高潮福利片在线看| 久久毛片免费看一区二区三区| 国产成人精品无人区| 精品午夜福利在线看| 2021少妇久久久久久久久久久| 亚洲国产精品国产精品| 成人亚洲欧美一区二区av| 精品国产乱码久久久久久小说| 黑人欧美特级aaaaaa片| 亚洲欧美激情在线| 亚洲成人国产一区在线观看 | 黑人猛操日本美女一级片| 老司机亚洲免费影院| 亚洲伊人色综图| 久久婷婷青草| 成年人免费黄色播放视频| videosex国产| 大话2 男鬼变身卡| 天天躁夜夜躁狠狠久久av| 亚洲精品视频女| 狂野欧美激情性bbbbbb| 久久精品久久久久久久性| 97在线人人人人妻| 久久久精品免费免费高清| 99国产精品免费福利视频| 精品国产乱码久久久久久男人| 在线观看免费日韩欧美大片| 成年女人毛片免费观看观看9 | 日韩av在线免费看完整版不卡| 亚洲欧洲日产国产| 欧美成人午夜精品| 捣出白浆h1v1| 精品一区二区三区四区五区乱码 | 日日撸夜夜添| 男女国产视频网站| 国产免费一区二区三区四区乱码| 久久久久久免费高清国产稀缺| 一级爰片在线观看| 伦理电影大哥的女人| 人人妻人人澡人人爽人人夜夜| 亚洲精品日韩在线中文字幕| 一级毛片 在线播放| 成人漫画全彩无遮挡| 我要看黄色一级片免费的| 在线免费观看不下载黄p国产| www.自偷自拍.com| 狠狠精品人妻久久久久久综合| 日韩 欧美 亚洲 中文字幕| 精品人妻在线不人妻| 欧美日韩亚洲综合一区二区三区_| 十八禁人妻一区二区| 成人国产麻豆网| 亚洲av成人精品一二三区| 欧美日本中文国产一区发布| 女人被躁到高潮嗷嗷叫费观| 精品人妻一区二区三区麻豆| 国产色婷婷99| 亚洲精品美女久久久久99蜜臀 | 国产一区二区三区综合在线观看| 男女免费视频国产| 操出白浆在线播放| 亚洲国产中文字幕在线视频| 国产1区2区3区精品| 亚洲精品成人av观看孕妇| 亚洲av在线观看美女高潮| 国产精品国产av在线观看| 午夜福利免费观看在线| 麻豆av在线久日| 国产成人欧美在线观看 | 麻豆av在线久日| 午夜老司机福利片| 制服人妻中文乱码| 亚洲欧美清纯卡通| 亚洲国产欧美日韩在线播放| 欧美黑人精品巨大| 亚洲av日韩精品久久久久久密 | 色网站视频免费| a级毛片在线看网站| 国产麻豆69| 国产日韩一区二区三区精品不卡| 日本欧美国产在线视频| 亚洲精品一区蜜桃| 亚洲成av片中文字幕在线观看| 伊人亚洲综合成人网| 久久ye,这里只有精品| 激情视频va一区二区三区| 成年人免费黄色播放视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美 亚洲 国产 日韩一| 人妻 亚洲 视频| 亚洲激情五月婷婷啪啪| 精品第一国产精品| 人人妻,人人澡人人爽秒播 | 别揉我奶头~嗯~啊~动态视频 | 免费观看性生交大片5| 亚洲av在线观看美女高潮| 51午夜福利影视在线观看| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩一级在线毛片| 毛片一级片免费看久久久久| 在线观看一区二区三区激情| 一区在线观看完整版| 99九九在线精品视频| 巨乳人妻的诱惑在线观看| 麻豆av在线久日| 最近中文字幕高清免费大全6| 女人高潮潮喷娇喘18禁视频| 色精品久久人妻99蜜桃| av卡一久久| 在线观看一区二区三区激情| 观看av在线不卡| av电影中文网址| 在线观看免费视频网站a站| 汤姆久久久久久久影院中文字幕| 一本久久精品| 久久天堂一区二区三区四区| 国产97色在线日韩免费| 亚洲欧美一区二区三区国产| 性高湖久久久久久久久免费观看| 一级毛片电影观看| 亚洲av在线观看美女高潮| 国产成人精品久久二区二区91 | 久久精品熟女亚洲av麻豆精品| 最新的欧美精品一区二区| 免费观看性生交大片5| 在线观看人妻少妇| 制服诱惑二区| 免费日韩欧美在线观看| 在线天堂中文资源库| 女人精品久久久久毛片| 99热全是精品| 色吧在线观看| 另类精品久久| 最近最新中文字幕大全免费视频 | 性色av一级| 夫妻性生交免费视频一级片| 久久热在线av| 国产精品久久久av美女十八| 成人黄色视频免费在线看| 中文字幕av电影在线播放| 亚洲第一av免费看| 日本午夜av视频| 中文字幕精品免费在线观看视频| 亚洲精品美女久久久久99蜜臀 | 欧美人与性动交α欧美软件| 久久久久精品人妻al黑| 久久久国产一区二区| 婷婷色av中文字幕| 国产女主播在线喷水免费视频网站| 久久热在线av| 美女午夜性视频免费| 国产精品偷伦视频观看了| 咕卡用的链子| 国产av一区二区精品久久| 美女主播在线视频| av女优亚洲男人天堂| 91精品三级在线观看| 操美女的视频在线观看| 少妇人妻精品综合一区二区| 夫妻性生交免费视频一级片| 人人妻人人添人人爽欧美一区卜| 在线观看一区二区三区激情| 亚洲欧美成人综合另类久久久| a 毛片基地| 亚洲精品美女久久av网站| 国产精品人妻久久久影院| 亚洲熟女毛片儿| 夫妻性生交免费视频一级片| 国产精品 欧美亚洲| 久久久久国产精品人妻一区二区| 成人午夜精彩视频在线观看| 国产熟女欧美一区二区| av福利片在线| 青春草亚洲视频在线观看| 97精品久久久久久久久久精品| 人妻 亚洲 视频| 国产女主播在线喷水免费视频网站| 日韩中文字幕视频在线看片| 精品久久久精品久久久| 国产乱人偷精品视频| 国产男女超爽视频在线观看| 精品一区二区免费观看| 精品国产国语对白av| 另类精品久久| 蜜桃在线观看..| 亚洲国产欧美一区二区综合| av一本久久久久| 少妇 在线观看| 欧美黄色片欧美黄色片| 别揉我奶头~嗯~啊~动态视频 | 十八禁人妻一区二区| 国产日韩一区二区三区精品不卡| 一级a爱视频在线免费观看| 国产一卡二卡三卡精品 | 国产1区2区3区精品| 亚洲国产欧美日韩在线播放| 亚洲激情五月婷婷啪啪| 超碰成人久久| 国产高清不卡午夜福利| 观看美女的网站| 精品国产一区二区三区四区第35| 日韩不卡一区二区三区视频在线| 亚洲精品中文字幕在线视频| 午夜日本视频在线| 亚洲欧美清纯卡通| 国产成人精品在线电影| 啦啦啦中文免费视频观看日本| h视频一区二区三区| 久久久精品国产亚洲av高清涩受| 妹子高潮喷水视频| 亚洲精品aⅴ在线观看| 热re99久久国产66热| 亚洲精品一二三| 亚洲精品乱久久久久久| 狠狠精品人妻久久久久久综合| 一级,二级,三级黄色视频| 99九九在线精品视频| 高清视频免费观看一区二区| 国产人伦9x9x在线观看| 欧美97在线视频| 成人国产麻豆网| 黄网站色视频无遮挡免费观看| 亚洲欧美精品自产自拍| 大片免费播放器 马上看| 久久天堂一区二区三区四区| 亚洲欧洲日产国产| 国产色婷婷99| 又大又黄又爽视频免费| 亚洲综合精品二区| 天天影视国产精品| 国产精品久久久久久人妻精品电影 | 波多野结衣一区麻豆| 久久99热这里只频精品6学生| 亚洲熟女毛片儿| 中文字幕人妻丝袜一区二区 | 国产精品久久久久久久久免| 久久久久久久久久久免费av| 人人妻人人澡人人看| 精品久久久精品久久久| 免费少妇av软件| 欧美日韩成人在线一区二区| 亚洲 欧美一区二区三区| 女的被弄到高潮叫床怎么办| 国产成人精品福利久久| 一区二区三区精品91| 一区二区三区四区激情视频| 久久久久久人妻| 80岁老熟妇乱子伦牲交| 午夜福利网站1000一区二区三区| 搡老乐熟女国产| 18禁观看日本| 中文字幕另类日韩欧美亚洲嫩草| 色综合欧美亚洲国产小说| 精品少妇一区二区三区视频日本电影 | 蜜桃在线观看..| 欧美日韩一级在线毛片| 欧美国产精品va在线观看不卡| 亚洲精华国产精华液的使用体验| 欧美日韩一级在线毛片| 高清欧美精品videossex| 19禁男女啪啪无遮挡网站| 丰满迷人的少妇在线观看| 成年动漫av网址| 最近最新中文字幕免费大全7| 欧美日韩亚洲高清精品| 制服诱惑二区| 老汉色∧v一级毛片| 1024视频免费在线观看| 精品少妇一区二区三区视频日本电影 | 亚洲av男天堂| 纯流量卡能插随身wifi吗| 亚洲精品第二区| 免费女性裸体啪啪无遮挡网站| 精品少妇内射三级| 久久女婷五月综合色啪小说| 亚洲,一卡二卡三卡| 欧美激情 高清一区二区三区| 免费av中文字幕在线| 亚洲国产毛片av蜜桃av| 成人国语在线视频| 国产精品久久久av美女十八| 精品午夜福利在线看| 精品亚洲成国产av| 各种免费的搞黄视频| 伊人久久国产一区二区| 操美女的视频在线观看| 97在线人人人人妻| 十分钟在线观看高清视频www| 久久精品国产a三级三级三级| 波多野结衣一区麻豆| 激情视频va一区二区三区| 亚洲av国产av综合av卡| 亚洲精品第二区| a级毛片在线看网站| 久久久久精品性色| a级毛片在线看网站| av女优亚洲男人天堂| 亚洲欧美一区二区三区黑人| 中文乱码字字幕精品一区二区三区| 香蕉丝袜av| 咕卡用的链子| 国产av精品麻豆| 嫩草影院入口| 99国产综合亚洲精品| 99久久综合免费| 性高湖久久久久久久久免费观看| 亚洲婷婷狠狠爱综合网| 欧美日韩综合久久久久久| 丝袜美足系列| 丝袜美腿诱惑在线| 欧美国产精品va在线观看不卡| 欧美日韩一区二区视频在线观看视频在线| 欧美人与性动交α欧美软件| 最近最新中文字幕大全免费视频 | 日韩视频在线欧美| 久久久国产一区二区| 精品国产一区二区三区四区第35| 亚洲国产中文字幕在线视频| 少妇被粗大的猛进出69影院| 亚洲欧洲精品一区二区精品久久久 | 国产成人精品久久二区二区91 | 天天添夜夜摸| 高清欧美精品videossex| 一级毛片黄色毛片免费观看视频| 我要看黄色一级片免费的| 欧美激情高清一区二区三区 | 亚洲av国产av综合av卡| 天堂俺去俺来也www色官网| 超色免费av| 亚洲精品成人av观看孕妇| 婷婷色麻豆天堂久久| 欧美成人午夜精品| 男人操女人黄网站| 中文字幕亚洲精品专区| 亚洲国产av影院在线观看| 国产深夜福利视频在线观看| 色吧在线观看| 久久狼人影院| 成人影院久久| 高清视频免费观看一区二区| 日日爽夜夜爽网站| 亚洲成人av在线免费| 女人爽到高潮嗷嗷叫在线视频| 天天添夜夜摸| 在现免费观看毛片| 女性生殖器流出的白浆| 熟女少妇亚洲综合色aaa.| 久久久久久久大尺度免费视频| 日本av免费视频播放| 老司机靠b影院| 操出白浆在线播放| 涩涩av久久男人的天堂| 又粗又硬又长又爽又黄的视频| avwww免费| 丰满乱子伦码专区| 中文字幕人妻丝袜一区二区 | 亚洲第一区二区三区不卡| 国产精品熟女久久久久浪| 欧美精品一区二区大全| 国产精品国产三级国产专区5o| 国产精品国产三级专区第一集| 久久久精品区二区三区| 久久精品久久久久久噜噜老黄| 天天躁夜夜躁狠狠躁躁| 日本一区二区免费在线视频| 国产精品一二三区在线看| 男女下面插进去视频免费观看| 日本爱情动作片www.在线观看| 精品一区二区三卡| 久久精品aⅴ一区二区三区四区| 国产麻豆69| 自线自在国产av| 香蕉国产在线看| 一级黄片播放器| 欧美精品人与动牲交sv欧美| 色播在线永久视频| 精品少妇黑人巨大在线播放| 一级片'在线观看视频| 亚洲欧美色中文字幕在线| 最近最新中文字幕免费大全7| 亚洲图色成人| 香蕉国产在线看| 在线免费观看不下载黄p国产| 成人毛片60女人毛片免费| 男女下面插进去视频免费观看| 亚洲综合精品二区| 免费在线观看完整版高清| 男女边吃奶边做爰视频| 精品国产一区二区三区四区第35| 成人国产麻豆网| 美女视频免费永久观看网站| 一级a爱视频在线免费观看| 国产免费现黄频在线看| 国产精品秋霞免费鲁丝片| 久久久久久人妻| 人人妻人人澡人人看| 老司机靠b影院| 中文欧美无线码| 建设人人有责人人尽责人人享有的| 久久久久久久久久久久大奶| 久久久久精品性色| 久久久久久久久久久久大奶| 一边摸一边抽搐一进一出视频| 国产成人精品久久二区二区91 | 亚洲美女搞黄在线观看| 国产极品粉嫩免费观看在线| 久久久久久人妻| 日本午夜av视频| 女的被弄到高潮叫床怎么办| 亚洲天堂av无毛| 国产无遮挡羞羞视频在线观看| 这个男人来自地球电影免费观看 | 十八禁人妻一区二区| 青春草国产在线视频| 又大又爽又粗| 亚洲精品成人av观看孕妇| 国产精品国产三级专区第一集| 日韩一区二区三区影片| 久久久亚洲精品成人影院| 曰老女人黄片| 男的添女的下面高潮视频| 国产精品熟女久久久久浪| 亚洲国产精品一区二区三区在线| 亚洲综合精品二区| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美网| 国产高清不卡午夜福利| 精品国产国语对白av| 最近中文字幕2019免费版| 国产免费一区二区三区四区乱码| 久久国产精品大桥未久av| 精品卡一卡二卡四卡免费| 亚洲,欧美,日韩| 校园人妻丝袜中文字幕| 国产伦理片在线播放av一区| 久久热在线av| 国产精品二区激情视频| 看免费av毛片| 不卡视频在线观看欧美| 黄色毛片三级朝国网站| 免费高清在线观看日韩| 亚洲精品国产av成人精品| 你懂的网址亚洲精品在线观看| netflix在线观看网站| 午夜福利免费观看在线| 久久精品国产综合久久久| 精品国产乱码久久久久久男人| 亚洲成人手机| 色网站视频免费| 高清欧美精品videossex| 亚洲国产精品一区二区三区在线| 老司机亚洲免费影院| 亚洲欧美一区二区三区黑人| 夜夜骑夜夜射夜夜干| 欧美日韩亚洲综合一区二区三区_| 人妻人人澡人人爽人人| 精品国产乱码久久久久久男人| 国语对白做爰xxxⅹ性视频网站| 精品国产一区二区三区久久久樱花| 9191精品国产免费久久| 99精国产麻豆久久婷婷| 精品久久久精品久久久| 国产有黄有色有爽视频| 中文字幕精品免费在线观看视频| 亚洲欧美精品自产自拍| 中文乱码字字幕精品一区二区三区| 国产成人精品无人区| 国产黄色视频一区二区在线观看| 男的添女的下面高潮视频| 国产精品免费视频内射| 久久精品熟女亚洲av麻豆精品| 欧美日韩成人在线一区二区| 少妇人妻精品综合一区二区| 久久天躁狠狠躁夜夜2o2o | av电影中文网址| 成人亚洲欧美一区二区av| 一级毛片我不卡| 欧美日韩一区二区视频在线观看视频在线| 一本久久精品| 亚洲七黄色美女视频| 国产精品久久久久久精品电影小说| 欧美在线一区亚洲| av女优亚洲男人天堂| 日韩av不卡免费在线播放| 少妇人妻精品综合一区二区| 汤姆久久久久久久影院中文字幕| 日韩中文字幕欧美一区二区 | 亚洲国产av影院在线观看| 高清在线视频一区二区三区| 国产精品亚洲av一区麻豆 | 精品久久蜜臀av无| 午夜影院在线不卡| 国产成人av激情在线播放| 国产伦人伦偷精品视频| 蜜桃国产av成人99| 又粗又硬又长又爽又黄的视频| 国产在线视频一区二区| 亚洲精品日本国产第一区| 中文字幕人妻熟女乱码| 午夜福利乱码中文字幕| 亚洲成人手机| 亚洲人成电影观看| 少妇精品久久久久久久| 国产成人av激情在线播放| 亚洲欧美精品自产自拍| 晚上一个人看的免费电影| 一边亲一边摸免费视频| 日韩熟女老妇一区二区性免费视频| 九色亚洲精品在线播放| 亚洲伊人久久精品综合| 黄片小视频在线播放| 夜夜骑夜夜射夜夜干| 国产精品一国产av| 日韩视频在线欧美| 亚洲 欧美一区二区三区| 欧美精品人与动牲交sv欧美| 亚洲图色成人| 久久ye,这里只有精品| 欧美日韩亚洲综合一区二区三区_| 精品国产超薄肉色丝袜足j| 精品午夜福利在线看| www.av在线官网国产| 亚洲欧美一区二区三区国产| 日本vs欧美在线观看视频| 久久av网站| 五月天丁香电影| 人人妻人人澡人人看| 久久天堂一区二区三区四区| 夜夜骑夜夜射夜夜干| 日本色播在线视频| 久久精品aⅴ一区二区三区四区| 精品卡一卡二卡四卡免费| 一区二区三区精品91| 免费少妇av软件| 国产成人免费无遮挡视频| 美女大奶头黄色视频| 高清在线视频一区二区三区| 欧美日韩视频高清一区二区三区二| 久久精品国产亚洲av高清一级| 久久久久久久国产电影| 中文字幕最新亚洲高清| 久久人妻熟女aⅴ| 中文字幕制服av| 丰满乱子伦码专区| 你懂的网址亚洲精品在线观看| 飞空精品影院首页| 欧美少妇被猛烈插入视频| 欧美日韩亚洲综合一区二区三区_| 婷婷成人精品国产| 男女边吃奶边做爰视频| 欧美日韩亚洲综合一区二区三区_| 丝袜在线中文字幕| www.熟女人妻精品国产| 免费女性裸体啪啪无遮挡网站| 国产欧美日韩一区二区三区在线| 亚洲精品国产色婷婷电影| 免费观看av网站的网址| 久久久久久久国产电影| 青春草亚洲视频在线观看| 满18在线观看网站| 国产在线一区二区三区精| 观看美女的网站| 色精品久久人妻99蜜桃| 亚洲国产毛片av蜜桃av| 久久韩国三级中文字幕| 久久性视频一级片| 青春草视频在线免费观看|