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

    編碼激勵在巖石衰減(1/Q)檢測中的應(yīng)用

    2022-01-25 07:06:40朱維王尚旭吳何珍翟鴻宇和泰名李長征
    地球物理學(xué)報 2022年1期
    關(guān)鍵詞:旁瓣信噪比超聲波

    朱維,王尚旭*,吳何珍,翟鴻宇,和泰名,李長征

    1 中國石油大學(xué)(北京)油氣資源與探測國家重點實驗室,北京 102249 2 中國地震局地球物理研究所,北京 100081 3 黃河水利科學(xué)研究院,鄭州 450003

    0 引言

    巖石超聲檢測技術(shù)是研究地震波傳播特征的重要物理模擬手段,實驗室超聲檢測結(jié)果能夠為人工地震和天然地震探測地下介質(zhì)變化提供參考和理論依據(jù),是實驗地震學(xué)的重要內(nèi)容.利用超聲波檢測巖石的聲學(xué)參數(shù),如彈性模量、聲速和衰減系數(shù)等,能夠獲得巖石內(nèi)部結(jié)構(gòu)和巖石強度.

    巖石聲學(xué)特性的廣泛深入研究開始于20世紀(jì)60年代,主要分為巖石超聲波速度和衰減的測量及其應(yīng)用研究.一方面,利用巖石超聲波速不僅可以計算巖石的彈性參數(shù)(Winkler and Liu,1996),還可監(jiān)測巖石所受應(yīng)力的變化(Gladwin,1982).另一方面,由于超聲波衰減對介質(zhì)物性變化的敏感性高于超聲波速度,因而能更精準(zhǔn)地反映巖石的內(nèi)部結(jié)構(gòu)及應(yīng)力狀態(tài).然而,衰減測量十分容易受到外部因素的影響,因而一直是實驗地震學(xué)研究的難點.Auberger和Rinehart(1961)利用250~1000 kHz的超聲波檢測了八種不同巖石的衰減系數(shù),發(fā)現(xiàn)以上頻段內(nèi)巖石的衰減與頻率并不滿足線性關(guān)系.Toks?z等(1979)通過對比干燥和飽和巖石的超聲波衰減說明流體會增加衰減.何國梁等(2007)利用超聲波衰減系數(shù)研究了砂巖高溫前后內(nèi)部裂隙的發(fā)展情況,認(rèn)為衰減可以用來評價巖石強度.梁利喜等(2015)通過數(shù)值實驗研究了超聲波衰減特性與孔洞結(jié)構(gòu)的關(guān)系,發(fā)現(xiàn)孔洞形狀、尺寸、分布和密度都對衰減有影響.然而,巖石超聲波衰減的測量十分依賴于檢測信號的信噪比,信噪比越高檢測結(jié)果就越可靠.

    通過提高信號發(fā)射功率可以增強信號能量并提高信噪比.但是,由于受到儀器參數(shù)和技術(shù)規(guī)范的限制,信號發(fā)射功率仍然有限.增加信號長度也能提高信號能量,但會降低信號的分辨率.因此,編碼激勵技術(shù)應(yīng)用而生,它通過發(fā)射一長串調(diào)制編碼信號提高平均峰值功率和增強信號能量;在接收端,它利用匹配濾波器或失配濾波器將時域上的長串信號壓縮為一段窄脈沖,來達(dá)到提高接收信號信噪比的同時增加系統(tǒng)分辨率的目的.目前編碼激勵技術(shù)已被應(yīng)用于醫(yī)學(xué)超聲成像 (O′Donnell,1992;Misaridis et al.,2000;Zhao et al.,2007;Fu et al.,2014;Zhou et al.,2014)、無損檢測(Wei et al.,2010;Li et al.,2013;Arora and Mulaveesala,2017;Lasaygues et al.,2018)以及地震探測領(lǐng)域(Bernhardt and Peacock,1978;Chapman et al.,1981;Barbier,1982;Park et al.,1996;Wong and Stewart,2006;葛洪魁等,2006).常用編碼信號可以被分為連續(xù)編碼,如頻率調(diào)制信號,和相位編碼,如巴克序列和正交Golay互補序列等(Chiao and Hao,2003;Zhao et al.,2007;Fu et al.,2014).其中,線性調(diào)頻信號(LFM)是最簡單且最重要的編碼,已被用來檢測人體內(nèi)部器官病變(O′Donnell,1992;Misaridis et al.,2000;Misaridis and Jensen,2005)、材料內(nèi)部缺陷(Wei et al.,2010;周正干等,2014)以及巖石的速度(李長征等,2013).LFM信號優(yōu)點突出,但對硬件要求較高,需要多級的發(fā)射器和用于匹配濾波的多級解碼濾波器.二進(jìn)制編碼信號只需要二相脈沖發(fā)生器和解碼器,因而成本更低(Zhao et al.,2007),最近也逐漸受到人們重視(Hu et al.,2006;Leavens et al.,2009;Mulaveesala and Ghali,2011;Zhou et al.,2014;Fu et al.,2014;Lu et al.,2015;朱維和吳何珍,2017;Zhang et al.,2018;Wu et al.,2020).由于Barker編碼具有完美的自相關(guān)性能,它已被用于醫(yī)學(xué)高頻超聲成像(劉凱和高上凱,2007)、混凝土缺陷檢測(李長征和王銳,2013)以及非接觸超聲無損檢測等方面.最近,它被引入到巖石的超聲波速度的精確檢測中(Wu et al.,2020;李長征等,2013;朱維和吳何珍,2017).

    然而,要將該項技術(shù)引入到巖石品質(zhì)因子Q值檢測中仍存在較大困難.朱維和吳何珍(2017)指出衰減引起的波形畸變會降低編碼信號的脈沖壓縮性能,可能會影響衰減的檢測.因此為了獲取準(zhǔn)確的Q值,必須選擇脈沖壓縮性能受衰減影響更小的編碼信號.為了評估衰減對編碼信號脈沖壓縮性能的影響,我們利用數(shù)值實驗對編碼激勵技術(shù)進(jìn)行了測試.通過使用編碼激勵方法的實驗室實例,給出了檢測巖石品質(zhì)因子Q值的流程,并介紹了如何利用衰減監(jiān)測三軸壓縮時巖石的變形過程.

    1 基本原理

    1.1 衰減檢測

    一般來說,穿過衰減介質(zhì)后,超聲波的頻域表達(dá)式可以寫為:

    W(ω)=S(ω)*Gtr(ω)*Hm(ω)+N(ω),

    (1)

    式中,*代表乘積,ω是角頻率,S(ω)是激勵信號頻譜,N(ω)是噪聲頻譜,Gtr(ω)和Hm(ω)分別是換能器和衰減介質(zhì)的傳遞函數(shù):

    (2)

    式中,c和Q分別是衰減介質(zhì)的速度和品質(zhì)因子,z為超聲波在介質(zhì)中的傳播距離.實驗室中,巖石樣品的衰減往往是相對于已知速度和衰減的參考樣品進(jìn)行測定的,這就是被人熟知的譜比法.如果選定的參考樣品近似為無衰減介質(zhì)(如鋁樣),則兩個樣品的振幅譜之比可以寫為:

    (3)

    式中,τ為透射系數(shù),下標(biāo)“r”代表參考樣品,SNR是信號與噪聲的振幅之比:

    SNR=|W(ω)|/|N(ω)|,

    (4)

    如果SNR很高,即1/SNR~0,那么方程(3)可以表示為:

    (5)

    式中,α為衰減系數(shù):

    (6)

    因此,衰減系數(shù)可以通過方程(5)的線性擬合之斜率k來確定.于是品質(zhì)因子Q就可以利用方程(6)計算得到.然而,當(dāng)SNR很低時噪聲的影響無法被忽略,斜率k不能代表衰減系數(shù).因此,低信噪比信號不能通過譜比法準(zhǔn)確檢測衰減.為了提高衰減檢測精度,提高接收信號信噪比成為關(guān)鍵.本研究將編碼激勵方法引入到巖石衰減檢測中,通過脈沖壓縮來提高接收信號信噪比.

    1.2 編碼激勵

    如圖1所示,編碼信號s(t)激勵換能器向巖石樣品發(fā)射超聲波,超聲波在穿透巖石后在另一端被接收換能器接收并轉(zhuǎn)換成電信號:

    圖1 編碼激勵超聲測試系統(tǒng)與脈沖壓縮原理Fig.1 Ultrasonic testing system for coded excitation and the principle of the pulse compression

    y(t)=(s(t)?hTtr(t)?hm(t))?hRtr(t)+n(t),

    (7)

    式中,?表示卷積,hm(t)為巖石樣品傳輸特性函數(shù),n(t)為隨機噪聲,hTtr(t)為發(fā)射換能器響應(yīng)函數(shù),hRtr(t)為接收換能器響應(yīng)函數(shù),通常二者是相等的.在得到接收信號y(t)后,需通過匹配濾波器對其進(jìn)行脈沖壓縮(或稱解碼).設(shè)匹配濾波器為spc(t),并令gtr=hTtr(t)?hRtr(t),則脈沖壓縮結(jié)果可表示為:

    ypc(t)=(s(t)?gtr(t)?hm(t)+n(t))?spc(t)

    =s(t)?spc(t)?gtr(t)?hm(t)+n(t)?spc(t),

    (8)

    理想情況下,

    s(t)?spc(t)=κδ(t),

    (9)

    式中κ為大于1的常數(shù),則(8)式可寫為:

    ypc(t)=k·gtr(t)?hm(t)+n(t)?spc(t),

    (10)

    式中n(t)?spc(t)即噪聲與接收信號的互相關(guān).由于它們是不相關(guān)的,因此經(jīng)過脈沖壓縮后噪聲得到了明顯的抑制.

    常用的編碼信號主要可以分為頻率編碼和相位編碼.每個編碼信號會使用不同的調(diào)制函數(shù).對于線性調(diào)頻信號(LFM),其調(diào)制函數(shù)可以寫為:

    s(t)=sin(2π(f0-B/2)t+πμt2),0≤t≤T

    (11)

    式中f0是換能器的中心頻率,B是LFM的帶寬,μ=B/T是掃頻速率,T是LFM的持續(xù)時間.圖2a展示了LFM信號的波形.對于相位編碼,其調(diào)制函數(shù)由(12)式給定:

    (12)

    式中,p(t)是決定了相位編碼信號的頻譜的載頻函數(shù),c是長度為N的編碼,Tp是每片編碼的持續(xù)時間.圖2e展示了正弦載頻調(diào)制的巴克編碼信號(BS).

    圖2 編碼信號及相應(yīng)脈沖壓縮信號(a)LFM信號;(b)匹配濾波后(a)的脈沖壓縮信號;(c)幅度削減的LFM信號;(d)匹配濾波后(c)的脈沖壓縮信號;(e)BS信號;(f)傳統(tǒng)匹配濾波后(e)的脈沖壓縮信號;(g)BS信號的逆濾波器;(h)利用逆濾波器匹配濾波后(e)的脈沖壓縮信號.Fig.2 Coded signals and their corresponding pulse compressed signals(a)LFM signal;(b)Pulse compressed signal of (a)after matched filtered;(c)Tapered LFM signal;(d)Pulse compressed signal of (c)after matched filtered;(e)BS signal;(f)Pulse compressed signal of (e)after traditionally matched filtered;and (g)Inverse filter of BS;(h)Pulse compressed signal of (e)after matched filtered by inverse filter.

    然而,不同編碼信號的脈沖壓縮性能是不同的,通??梢杂尚旁氡仍鲆?GSNR)、距離旁瓣水平(RSL)以及主瓣寬度(MLW)進(jìn)行評價.信噪比增益通常是指信號匹配濾波前后信噪比的增加值,理論上GSNR由編碼信號的時間帶寬積(TBP)決定:

    GSNR=10log10(TBP)(dB).

    (13)

    一般而言,編碼信號通過匹配濾波后,在主瓣附近還會出現(xiàn)比較明顯的旁瓣,為了衡量主瓣與旁瓣的相對大小可以用峰值旁瓣水平(PSL)進(jìn)行評價.PSL越低說明編碼信號的脈沖壓縮性能越好.

    (14)

    但是旁瓣往往是無法避免的,一般采用數(shù)學(xué)方法抑制這些旁瓣對速度或Q值檢測的影響.對頻率編碼信號(如LFM信號),需要采用諸如振幅削減技術(shù)的方法對發(fā)射源進(jìn)行處理,從而達(dá)到抑制旁瓣的目的,但會損失一定的軸向分辨率.根據(jù)圖2b可以看到其脈沖壓縮波形存在嚴(yán)重旁瓣.圖2c為振幅削減后的線性調(diào)頻信號(TLFM),從圖2d可以看到其脈沖壓縮波形旁瓣得到了明顯抑制.對于相位編碼信號(如Barker編碼),可以通過設(shè)計新的匹配濾波器(如偽逆濾波器)抑制旁瓣;圖2f展示了BS信號通過傳統(tǒng)匹配濾波后的脈沖壓縮波形,而圖2h則展示了利用偽逆濾波器濾波后的脈沖壓縮波形,可以看到旁瓣得到了明顯抑制.最后一個評價指標(biāo)——主瓣寬度一般是指編碼信號自相關(guān)函數(shù)主瓣的寬度,它是衡量信號軸向分辨率的重要指標(biāo),主瓣寬度越窄說明系統(tǒng)的分辨率越高,因此編碼信號經(jīng)過脈沖壓縮后必須具有較窄主瓣寬度.

    選擇合適的編碼信號以后,我們就可以進(jìn)行衰減檢測了.利用編碼信號檢測超聲波衰減的基本流程如圖3所示.

    圖3 利用編碼激勵技術(shù)檢測巖石衰減的流程圖Fig.3 Flow chart of the measurement of the attenuation in rock using coded excitation technique

    2 數(shù)值實驗

    由于衰減會嚴(yán)重影響編碼信號的脈沖壓縮性能,對不同編碼信號檢測衰減的能力評估是十分必要的.我們選擇正弦信號(Sine)、BS信號以及TLFM信號作為激勵信號生成了40組接收波形,其Q值范圍為5~200,合成信號都加入了相同功率的噪聲.圖4展示了Q=5時利用不同激發(fā)源生成的接收波形及相應(yīng)的脈沖壓縮波形.從圖中可以看到,脈沖壓縮后噪聲被明顯抑制.

    圖4 Q=5時,不同激勵源激發(fā)得到的理論接收信號及脈沖壓縮信號(a)、(b)和(d)分別是Sine信號、TLFM和BS信號作為激勵源的接收信號;(c)和(e)分別是匹配濾波后(b)的脈沖壓縮信號和利用逆濾波器匹配濾波后(d)的脈沖壓縮信號.Fig.4 The synthetic received signals and pulse compressed signals obtained by different excitation sources when Q=5Received signals excited with (a)Sine signal,(b)TLFM signal and (d)BS signal;(c)Pulse compressed signal of (b)after matched filtered;(e)Pulse compressed signal of (d)after matched filtered by inverse filter.

    圖5是衰減對編碼信號脈沖壓縮性能影響的評估結(jié)果.可以看到,衰減會嚴(yán)重影響編碼信號的脈沖壓縮性能.一方面,TLFM信號的GSNR會隨衰減的增大而減小,而BS信號卻并不受衰減的影響.由于不同頻率的能量衰減不同,高頻信息容易缺少或丟失,從而導(dǎo)致接收信號的中心頻率向低頻偏移,因而導(dǎo)致頻率編碼接收波形畸變十分嚴(yán)重;而相位編碼信號的頻率信息包含在載波信號中,因此相位編碼激勵源接收信號失配較小.另一方面,當(dāng)Q值>5時,BS信號的GSNR和MLW都大于TLFM信號;當(dāng)Q=5時,則正好相反.然而,由于TLFM信號的信噪比增益明顯低于BS信號,利用編碼信號求取Q值時,TLFM信號不如BS信號適合Q值的檢測.

    為了進(jìn)一步檢測編碼信號對衰減的檢測能力,我們對這些波形采用譜比法求取了Q值,如圖5d所示.根據(jù)計算結(jié)果,BS與TLFM激勵源檢測的Q值誤差都低于正弦激勵源的結(jié)果,這表明編碼激勵源能夠提高Q值的檢測精度.同時,從圖5d可以看到,Q<20時,TLFM激勵源計算的Q值誤差大于BS激勵源.這是因為衰減較大時,TLFM信號的脈沖壓縮性能受衰減影響較大,接收波形失配嚴(yán)重,從而導(dǎo)致了更大的Q值的計算誤差,因此TLFM信號并不適合應(yīng)用于較大衰減介質(zhì)的Q值檢測.當(dāng)Q>20時,TLFM激勵源計算的Q值誤差小于BS激勵源,這是因為TLFM激勵源的接收信號經(jīng)過脈沖壓縮后保存了更多高頻能量的緣故.綜合以上結(jié)果可知,編碼激勵源能提高Q值的檢測精度;對于衰減較大的介質(zhì),BS編碼激勵源檢測的Q值更準(zhǔn)確,而對衰減較小的介質(zhì),TLFM激勵源檢測的Q值更準(zhǔn)確.

    圖5 利用TLFM和BS激勵源合成的接收波形計算得到的脈沖壓縮波形的(a)信噪比增益,(b)峰值旁瓣水平,(c)主瓣寬度以及(d)Q值檢測相對誤差隨品質(zhì)因子的變化Fig.5 (a)The gain in SNR (GSNR),(b)peak sidelobes level (PSL)and (c)main-lobe width (MLW)of the waveforms compressed from the synthetic signals excited by TLFM and BS,and (d)the relative errors of Q varied with quality factor

    3 實際應(yīng)用

    我們利用編碼激勵技術(shù)測試了不同樣品的Q 值并監(jiān)測了巖石三軸加載過程時的衰減變化.

    3.1 測試1:巖石Q值檢測

    本測試采用中心頻率為1 MHz的超聲波換能器測試了有機玻璃、人工砂巖(環(huán)砂比為3∶8)及天然砂巖(采自重慶武隆)的Q值(圖6).利用單脈沖信號(這里采用Sine信號)和編碼信號檢測的Q值列于表1中.表中還列出了利用標(biāo)準(zhǔn)測試方法獲得的Q值.從表中可以看到對于有機玻璃,TLFM編碼激勵源檢測的Q值最接近參考值,而對于人工砂巖和天然砂巖BS激勵源檢測的Q值最接近參考值.這與數(shù)值實驗的結(jié)論是一致的.

    3.2 測試2:巖石破裂過程中的衰減監(jiān)測

    巖石破裂過程中內(nèi)部會經(jīng)歷舊微裂紋閉合、新微裂紋產(chǎn)生、集結(jié)并最終發(fā)生宏觀破裂等過程.由于巖石微結(jié)構(gòu)的變化會引起巖石衰減的變化,因此通過監(jiān)測衰減就能反映巖石受壓過程中內(nèi)部微結(jié)構(gòu)的變化.本測試對測試1中天然砂巖開展了三軸壓縮實驗.加載設(shè)備為SANS電液伺服壓力機,其最大載荷為300 kN.巖石樣品置于圍壓容器中保持5 MPa圍壓不變,軸壓采用應(yīng)力控制模式(加載速率300 N·s-1).軸向位移隨應(yīng)力的變化如圖7中黑線所示.為了監(jiān)測巖石變形過程中超聲波衰減的變化,我們將換能器置于巖石樣品上下兩端用于發(fā)射和接收超聲波形.整個加載過程共記錄23組不同激勵源的檢測波形.利用這些波形我們計算了不同加載階段巖石的衰減(圖7).

    圖7 品質(zhì)因子Q值及軸向位移隨應(yīng)力的變化QSine,QBS和QTLFM分別代表Sine信號、BS信號和TLFM信號激勵源檢測波形得到的Q值;Da代表軸向位移.Fig.7 Quality factor Q and axial displacement varied with the stressQSine,QBS and QTLFM represent the quality factor measured with the signals excited by Sine,BS and TLFM,respectively.Da represents the axial displacement.

    根據(jù)前人研究,巖石變形過程中,衰減會經(jīng)歷先變小達(dá)到一定值后保持穩(wěn)定最后由于新裂紋的產(chǎn)生又逐漸減小的過程(Lockner et al.,1977).從圖7中衰減隨應(yīng)力的變化曲線來看,三種信號都能很好地反映這一規(guī)律.這說明編碼信號可以用于監(jiān)測受載巖石的變形和破裂過程.同時,實驗結(jié)果表明應(yīng)力小于25 MPa時,巖石的Q值小于20且TLFM激勵源檢測的Q值明顯小于單脈沖和BS激勵源的檢測結(jié)果.這再次驗證了數(shù)值實驗得出的結(jié)論——衰減較大時TLFM激勵源的檢測結(jié)果誤差更大.由于巖石受壓變形過程中Q值變化較大,要在整個過程中都獲得更小誤差的Q值,這就要求編碼信號具有更穩(wěn)定的脈沖壓縮性能.因此Barker編碼信號比TLFM信號更適合監(jiān)測巖石的變形過程.

    4 結(jié)論

    利用編碼激勵技術(shù)檢測巖石的Q值是可行的,結(jié)果是可靠的.通過對編碼激勵與脈沖壓縮的基本原理、方法、實驗和應(yīng)用的系統(tǒng)分析,得出編碼激勵技術(shù)可以通過脈沖壓縮提高信噪比,但是巖石的衰減會削弱編碼信號的脈沖壓縮性能的結(jié)論.雖然幅度削減的線性調(diào)頻信號的脈沖壓縮性能比正弦調(diào)制的Barker編碼更好,但是當(dāng)Q值小于20時即巖石衰減較大時,TLFM的信噪比增益下降更為嚴(yán)重,TLFM激勵源檢測的Q值誤差更大.我們利用正弦信號,正弦調(diào)制的Barker編碼信號和幅度削減的線性調(diào)頻信號測量了有機玻璃、人工砂巖以及天然砂巖Q值.實驗結(jié)果表明編碼激勵技術(shù)能夠通過脈沖壓縮來抑制噪聲,提高Q值檢測精度.另外,我們還利用編碼激勵源檢測了砂巖在三軸壓縮過程中Q值的變化并監(jiān)測了巖石的變形和破裂過程.實驗結(jié)果表明,與TLFM激勵源相比BS激勵源在衰減較大時檢測的Q值誤差更小,因而更適合監(jiān)測衰減不斷變化的受壓巖石的變形過程.

    盡管我們已經(jīng)論證了編碼激勵在超聲頻段測試巖石Q值時的可靠性與適用性,但是當(dāng)巖石中黏土礦物含量較高或包含流體時,巖石的彈性波速度和Q值都可能出現(xiàn)頻散效應(yīng).當(dāng)然,對于超聲頻段的彈性波而言,這種頻散可能只有在巖石衰減很強時才會顯著.但是,由于頻散效應(yīng)的存在,要將編碼激勵在超聲頻段的應(yīng)用推廣至地震頻段仍需進(jìn)一步的研究.

    猜你喜歡
    旁瓣信噪比超聲波
    基于圓柱陣通信系統(tǒng)的廣義旁瓣對消算法
    一種基于線性規(guī)劃的頻率編碼旁瓣抑制方法
    基于深度學(xué)習(xí)的無人機數(shù)據(jù)鏈信噪比估計算法
    基于Niosll高精度超聲波流量計的研究
    電子制作(2018年17期)2018-09-28 01:56:38
    低信噪比下LFMCW信號調(diào)頻參數(shù)估計
    電子測試(2018年11期)2018-06-26 05:56:02
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    基于加權(quán)積分旁瓣最小化的隨機多相碼設(shè)計
    蝙蝠的超聲波
    超聲波流量計的研究
    電子制作(2016年21期)2016-05-17 03:52:47
    超聲波流量計在蘭干渠上的應(yīng)用
    国产又黄又爽又无遮挡在线| 国产高清三级在线| 我的女老师完整版在线观看| 亚洲专区国产一区二区| a在线观看视频网站| 国产伦人伦偷精品视频| 亚洲国产欧美人成| 日韩免费av在线播放| 在线观看66精品国产| 亚洲久久久久久中文字幕| 免费搜索国产男女视频| 噜噜噜噜噜久久久久久91| 日本黄色视频三级网站网址| 亚洲第一欧美日韩一区二区三区| 男女床上黄色一级片免费看| 最近中文字幕高清免费大全6 | 久久精品影院6| 久久精品夜夜夜夜夜久久蜜豆| 日韩精品青青久久久久久| 男人舔奶头视频| 精品人妻偷拍中文字幕| 伦理电影大哥的女人| 免费看光身美女| 老司机深夜福利视频在线观看| av专区在线播放| 啦啦啦观看免费观看视频高清| 午夜激情欧美在线| 亚洲av免费高清在线观看| 国产精品久久久久久精品电影| 好男人电影高清在线观看| 亚洲欧美日韩无卡精品| 欧美黑人欧美精品刺激| 99国产精品一区二区蜜桃av| 午夜激情福利司机影院| 国产 一区 欧美 日韩| 日本在线视频免费播放| 午夜视频国产福利| 在线天堂最新版资源| 国产三级在线视频| 精品乱码久久久久久99久播| 日本 av在线| 久久精品综合一区二区三区| 噜噜噜噜噜久久久久久91| 一区二区三区四区激情视频 | 欧美乱妇无乱码| 无人区码免费观看不卡| 一本精品99久久精品77| 99久久精品一区二区三区| 色综合婷婷激情| .国产精品久久| 久久精品人妻少妇| 他把我摸到了高潮在线观看| 最近中文字幕高清免费大全6 | www日本黄色视频网| 亚洲精品粉嫩美女一区| 男人的好看免费观看在线视频| 非洲黑人性xxxx精品又粗又长| 国产精品久久电影中文字幕| 欧美+日韩+精品| av在线蜜桃| www日本黄色视频网| 天堂√8在线中文| 女人十人毛片免费观看3o分钟| 最近在线观看免费完整版| 国产亚洲欧美98| 国产亚洲欧美98| 两个人视频免费观看高清| 观看免费一级毛片| 女人被狂操c到高潮| 欧美性感艳星| 久久草成人影院| 俺也久久电影网| 日本免费a在线| 成人毛片a级毛片在线播放| 制服丝袜大香蕉在线| 免费电影在线观看免费观看| 国产高清激情床上av| 日本与韩国留学比较| 久久热精品热| 国内少妇人妻偷人精品xxx网站| 国产精品一区二区三区四区免费观看 | 亚洲成人中文字幕在线播放| 亚洲自拍偷在线| 精品一区二区三区人妻视频| 国产精品一及| 少妇的逼好多水| 一夜夜www| 青草久久国产| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产成+人综合+亚洲专区| 三级国产精品欧美在线观看| 男女做爰动态图高潮gif福利片| 国产一区二区在线av高清观看| 一区二区三区激情视频| 亚洲五月婷婷丁香| 欧美在线黄色| 亚洲美女黄片视频| av在线观看视频网站免费| 午夜福利在线在线| 97碰自拍视频| 又粗又爽又猛毛片免费看| 国产精品一区二区三区四区免费观看 | 桃红色精品国产亚洲av| 美女cb高潮喷水在线观看| 亚洲av不卡在线观看| 亚洲欧美日韩高清专用| 夜夜夜夜夜久久久久| 国产精品亚洲美女久久久| 久久人人精品亚洲av| 亚洲七黄色美女视频| 亚洲人成网站在线播| 在线观看66精品国产| 亚洲自偷自拍三级| 又爽又黄无遮挡网站| 欧美激情久久久久久爽电影| 亚洲熟妇中文字幕五十中出| 久久久国产成人精品二区| 成人高潮视频无遮挡免费网站| 久久精品夜夜夜夜夜久久蜜豆| 久久久久九九精品影院| 中文字幕av在线有码专区| 日本与韩国留学比较| 69av精品久久久久久| 国产大屁股一区二区在线视频| 精品一区二区三区人妻视频| 天美传媒精品一区二区| 淫妇啪啪啪对白视频| 久久久久久九九精品二区国产| 舔av片在线| 久久草成人影院| av福利片在线观看| 国产精品爽爽va在线观看网站| 日本黄色片子视频| 亚洲av五月六月丁香网| 欧美精品国产亚洲| 久久草成人影院| 国产精品1区2区在线观看.| 深夜精品福利| 首页视频小说图片口味搜索| h日本视频在线播放| 97超视频在线观看视频| 最新在线观看一区二区三区| 在线观看免费视频日本深夜| 久久精品综合一区二区三区| 色尼玛亚洲综合影院| 身体一侧抽搐| 精品免费久久久久久久清纯| 欧美一级a爱片免费观看看| 亚洲精品亚洲一区二区| 国产精品久久视频播放| 亚洲成av人片在线播放无| 亚洲精品成人久久久久久| 日韩欧美精品免费久久 | 国产一区二区在线av高清观看| 日本黄色片子视频| 亚洲精品在线美女| 18禁在线播放成人免费| 色尼玛亚洲综合影院| 悠悠久久av| 在线观看美女被高潮喷水网站 | 亚洲成人精品中文字幕电影| 午夜福利视频1000在线观看| 亚洲,欧美,日韩| 他把我摸到了高潮在线观看| 成年免费大片在线观看| 不卡一级毛片| 乱人视频在线观看| 亚洲在线自拍视频| 免费看a级黄色片| 哪里可以看免费的av片| 美女cb高潮喷水在线观看| 天堂影院成人在线观看| 一级黄色大片毛片| www.色视频.com| 免费人成视频x8x8入口观看| 亚洲av电影不卡..在线观看| 国产精品不卡视频一区二区 | 两人在一起打扑克的视频| 真人做人爱边吃奶动态| 亚洲专区中文字幕在线| 五月玫瑰六月丁香| 桃红色精品国产亚洲av| 特大巨黑吊av在线直播| 国产熟女xx| 国产av在哪里看| 乱人视频在线观看| 亚洲人成电影免费在线| av在线老鸭窝| 欧美不卡视频在线免费观看| 男人舔奶头视频| 亚洲国产欧美人成| 国产精品嫩草影院av在线观看 | av专区在线播放| 国产熟女xx| 每晚都被弄得嗷嗷叫到高潮| 亚州av有码| 亚洲五月天丁香| 99久国产av精品| 亚洲一区二区三区不卡视频| av天堂在线播放| 国产精品自产拍在线观看55亚洲| 无人区码免费观看不卡| 久久午夜福利片| 精品午夜福利视频在线观看一区| 久久精品国产亚洲av涩爱 | 欧美黄色片欧美黄色片| 亚洲精品久久国产高清桃花| 人人妻人人看人人澡| 亚洲片人在线观看| 久久久久性生活片| 成人高潮视频无遮挡免费网站| 国内少妇人妻偷人精品xxx网站| 国产av一区在线观看免费| 伊人久久精品亚洲午夜| 久久婷婷人人爽人人干人人爱| 中亚洲国语对白在线视频| 天堂影院成人在线观看| 99久久九九国产精品国产免费| 久久精品久久久久久噜噜老黄 | 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩精品亚洲av| bbb黄色大片| 精品国内亚洲2022精品成人| 亚洲中文字幕日韩| 国产精品一区二区三区四区久久| 简卡轻食公司| 欧美日本视频| 午夜福利欧美成人| 国产精品久久久久久人妻精品电影| 精品午夜福利视频在线观看一区| 99精品在免费线老司机午夜| 免费av不卡在线播放| 日韩成人在线观看一区二区三区| 国产精品日韩av在线免费观看| 亚洲欧美日韩高清专用| 亚洲三级黄色毛片| 免费在线观看影片大全网站| 床上黄色一级片| 日韩 亚洲 欧美在线| 国产av麻豆久久久久久久| 国产高潮美女av| 亚洲av日韩精品久久久久久密| 欧美日韩福利视频一区二区| 成人精品一区二区免费| 男人舔奶头视频| 俄罗斯特黄特色一大片| 久久这里只有精品中国| 亚洲人与动物交配视频| 真实男女啪啪啪动态图| 在线天堂最新版资源| 国产一区二区在线av高清观看| 成年免费大片在线观看| 国产亚洲精品久久久久久毛片| 亚洲狠狠婷婷综合久久图片| 免费看日本二区| 国产精品亚洲美女久久久| 国产视频内射| 久久久色成人| 老女人水多毛片| 日本黄色视频三级网站网址| 免费高清视频大片| 午夜激情福利司机影院| 亚洲,欧美,日韩| 搡女人真爽免费视频火全软件 | 丰满人妻熟妇乱又伦精品不卡| 熟妇人妻久久中文字幕3abv| 两个人视频免费观看高清| 午夜福利成人在线免费观看| 国产伦精品一区二区三区四那| 波多野结衣高清无吗| 欧美成人性av电影在线观看| 高清在线国产一区| 在线a可以看的网站| 天堂动漫精品| 男人狂女人下面高潮的视频| 757午夜福利合集在线观看| 草草在线视频免费看| 91久久精品国产一区二区成人| 中文字幕免费在线视频6| 国产成+人综合+亚洲专区| 婷婷丁香在线五月| 欧美绝顶高潮抽搐喷水| 在线观看一区二区三区| 18禁黄网站禁片午夜丰满| 日韩精品中文字幕看吧| 亚洲精品粉嫩美女一区| 亚洲第一欧美日韩一区二区三区| 国产精品1区2区在线观看.| www.999成人在线观看| 人人妻,人人澡人人爽秒播| 一个人观看的视频www高清免费观看| 啦啦啦韩国在线观看视频| 深夜a级毛片| 少妇的逼好多水| 国产69精品久久久久777片| 熟女电影av网| 听说在线观看完整版免费高清| 国产在视频线在精品| 午夜影院日韩av| 亚洲最大成人手机在线| 国产探花极品一区二区| 高清日韩中文字幕在线| 美女高潮的动态| 999久久久精品免费观看国产| 久久这里只有精品中国| 好男人在线观看高清免费视频| 中文字幕av在线有码专区| 日本 av在线| 日本撒尿小便嘘嘘汇集6| 国产老妇女一区| 免费看美女性在线毛片视频| 国产在线精品亚洲第一网站| 日韩成人在线观看一区二区三区| 好男人在线观看高清免费视频| 欧美高清成人免费视频www| 久久国产乱子免费精品| 国产精品嫩草影院av在线观看 | 欧美3d第一页| 国产精品久久久久久久电影| 美女大奶头视频| 给我免费播放毛片高清在线观看| 淫妇啪啪啪对白视频| av天堂在线播放| 亚洲av一区综合| 国产精品久久电影中文字幕| 国产一级毛片七仙女欲春2| 欧美性猛交黑人性爽| 欧美成人免费av一区二区三区| 日本黄大片高清| 国内毛片毛片毛片毛片毛片| 中文在线观看免费www的网站| 91久久精品电影网| 黄色女人牲交| av视频在线观看入口| 日韩欧美在线乱码| 久久精品国产亚洲av香蕉五月| 久久九九热精品免费| 免费在线观看日本一区| 欧美日韩黄片免| 亚洲av.av天堂| 最近在线观看免费完整版| 久久久久久久久久成人| 日韩欧美三级三区| 欧美色欧美亚洲另类二区| av天堂中文字幕网| 国产黄a三级三级三级人| 三级毛片av免费| 欧美最新免费一区二区三区 | 午夜激情福利司机影院| 国产黄片美女视频| 国产一区二区亚洲精品在线观看| 极品教师在线视频| 国产精品一区二区性色av| 一级毛片久久久久久久久女| 真实男女啪啪啪动态图| 丁香欧美五月| 午夜福利在线观看吧| 日本熟妇午夜| 丝袜美腿在线中文| 赤兔流量卡办理| 亚洲综合色惰| 日本 欧美在线| 午夜a级毛片| 男女之事视频高清在线观看| 国产精品伦人一区二区| 成人国产综合亚洲| 国产黄a三级三级三级人| 国产亚洲欧美98| 日本黄色视频三级网站网址| 欧美一区二区亚洲| 国产在线精品亚洲第一网站| 亚洲真实伦在线观看| 国产黄色小视频在线观看| 欧美成人免费av一区二区三区| 成人国产综合亚洲| 亚洲精品色激情综合| 麻豆av噜噜一区二区三区| 国产精品1区2区在线观看.| 自拍偷自拍亚洲精品老妇| 超碰av人人做人人爽久久| 丰满人妻熟妇乱又伦精品不卡| 此物有八面人人有两片| 日韩高清综合在线| 国产精品亚洲美女久久久| 男女床上黄色一级片免费看| 精品午夜福利视频在线观看一区| 日韩免费av在线播放| 麻豆成人午夜福利视频| 日韩欧美国产一区二区入口| 美女高潮的动态| 亚洲国产精品999在线| 国产精品伦人一区二区| 国产精品av视频在线免费观看| 日本在线视频免费播放| 欧美黄色片欧美黄色片| av专区在线播放| 午夜福利免费观看在线| 老司机深夜福利视频在线观看| 91在线精品国自产拍蜜月| 国内精品一区二区在线观看| 一边摸一边抽搐一进一小说| 国产一级毛片七仙女欲春2| 国产淫片久久久久久久久 | 噜噜噜噜噜久久久久久91| 久久久久免费精品人妻一区二区| 一个人免费在线观看电影| 少妇熟女aⅴ在线视频| 亚洲黑人精品在线| 99视频精品全部免费 在线| 久久九九热精品免费| 国产高潮美女av| 少妇人妻精品综合一区二区 | 久久久精品欧美日韩精品| 婷婷色综合大香蕉| 日韩av在线大香蕉| 老女人水多毛片| 成人国产综合亚洲| 午夜福利在线在线| 午夜老司机福利剧场| 欧美激情国产日韩精品一区| 免费高清视频大片| 精品久久久久久久久久久久久| 午夜免费成人在线视频| 国产淫片久久久久久久久 | 国产午夜精品久久久久久一区二区三区 | 午夜久久久久精精品| 国产真实伦视频高清在线观看 | 亚洲avbb在线观看| 国产av在哪里看| 久99久视频精品免费| 国产精品久久久久久久电影| 91久久精品国产一区二区成人| 亚洲无线观看免费| 国产黄a三级三级三级人| 国产色爽女视频免费观看| 亚洲美女视频黄频| 日韩欧美免费精品| 免费在线观看日本一区| 91麻豆精品激情在线观看国产| 亚洲av电影不卡..在线观看| 国产野战对白在线观看| 精品欧美国产一区二区三| av天堂中文字幕网| 波多野结衣高清作品| 91麻豆av在线| 国产男靠女视频免费网站| 性插视频无遮挡在线免费观看| www.熟女人妻精品国产| 一区二区三区高清视频在线| 欧美xxxx黑人xx丫x性爽| 少妇熟女aⅴ在线视频| 哪里可以看免费的av片| 97超级碰碰碰精品色视频在线观看| 免费av不卡在线播放| 美女xxoo啪啪120秒动态图 | 99国产极品粉嫩在线观看| 国产真实乱freesex| 国内精品美女久久久久久| 老鸭窝网址在线观看| 欧美+日韩+精品| 色综合站精品国产| 国产高清三级在线| 午夜福利在线在线| 亚洲最大成人av| 免费人成在线观看视频色| 久久亚洲精品不卡| 首页视频小说图片口味搜索| 国产精品av视频在线免费观看| 国产成人a区在线观看| a级毛片免费高清观看在线播放| 亚洲无线观看免费| 在线观看午夜福利视频| ponron亚洲| 又爽又黄无遮挡网站| 色综合亚洲欧美另类图片| 99热精品在线国产| 久久久精品欧美日韩精品| 人妻制服诱惑在线中文字幕| 天天躁日日操中文字幕| av在线蜜桃| 色视频www国产| 亚洲国产精品sss在线观看| 别揉我奶头 嗯啊视频| 久久人人爽人人爽人人片va | 亚洲自偷自拍三级| 好看av亚洲va欧美ⅴa在| 免费搜索国产男女视频| 国产高潮美女av| 国产单亲对白刺激| 久久中文看片网| 国产激情偷乱视频一区二区| 日日干狠狠操夜夜爽| 动漫黄色视频在线观看| 亚洲成a人片在线一区二区| 国产精品av视频在线免费观看| 人妻丰满熟妇av一区二区三区| 我要看日韩黄色一级片| 三级男女做爰猛烈吃奶摸视频| 日韩亚洲欧美综合| 美女大奶头视频| 久久国产精品人妻蜜桃| 国产一区二区三区在线臀色熟女| 精品福利观看| 国产精品美女特级片免费视频播放器| 欧美在线黄色| 国产色婷婷99| 亚洲七黄色美女视频| 精品一区二区三区视频在线观看免费| 99久久精品国产亚洲精品| 国产主播在线观看一区二区| 窝窝影院91人妻| 两个人视频免费观看高清| 国产一级毛片七仙女欲春2| 最新在线观看一区二区三区| 91麻豆精品激情在线观看国产| 少妇的逼水好多| 99久久无色码亚洲精品果冻| 国产免费一级a男人的天堂| 桃色一区二区三区在线观看| 精品一区二区免费观看| 小蜜桃在线观看免费完整版高清| 久久久久久久精品吃奶| 国产精品综合久久久久久久免费| 国产v大片淫在线免费观看| 黄片小视频在线播放| 久久久久免费精品人妻一区二区| 日韩欧美精品免费久久 | 欧美日韩亚洲国产一区二区在线观看| 欧美精品国产亚洲| 亚洲国产精品合色在线| 国产欧美日韩一区二区精品| 一级a爱片免费观看的视频| 成人美女网站在线观看视频| 淫秽高清视频在线观看| 久久国产精品影院| 在线看三级毛片| 国产精品一及| 国产成+人综合+亚洲专区| 美女cb高潮喷水在线观看| 午夜福利18| 亚洲av成人精品一区久久| 观看美女的网站| 日韩中文字幕欧美一区二区| 欧美丝袜亚洲另类 | 亚洲欧美日韩卡通动漫| 一区二区三区高清视频在线| 亚洲第一欧美日韩一区二区三区| 久久精品国产清高在天天线| 91在线观看av| 欧美黑人巨大hd| 中文字幕av成人在线电影| 亚洲不卡免费看| 琪琪午夜伦伦电影理论片6080| h日本视频在线播放| 在线免费观看不下载黄p国产 | 中国美女看黄片| 午夜免费成人在线视频| 最近最新中文字幕大全电影3| 国内毛片毛片毛片毛片毛片| 亚洲av五月六月丁香网| 午夜免费激情av| 老熟妇仑乱视频hdxx| 色综合亚洲欧美另类图片| 久久精品国产99精品国产亚洲性色| 久久欧美精品欧美久久欧美| 美女高潮的动态| 3wmmmm亚洲av在线观看| 免费黄网站久久成人精品 | 国产亚洲精品久久久com| 免费在线观看影片大全网站| 欧美xxxx性猛交bbbb| 最好的美女福利视频网| 亚洲欧美精品综合久久99| avwww免费| 亚洲人成网站高清观看| 欧美成人性av电影在线观看| 国产乱人视频| 乱码一卡2卡4卡精品| 少妇裸体淫交视频免费看高清| 人妻夜夜爽99麻豆av| 亚洲美女搞黄在线观看 | 波多野结衣高清作品| 少妇高潮的动态图| 亚洲最大成人手机在线| 成人国产一区最新在线观看| 国产精品伦人一区二区| 免费av毛片视频| 丁香六月欧美| av视频在线观看入口| 三级国产精品欧美在线观看| 欧美丝袜亚洲另类 | 国产欧美日韩一区二区三| 天堂网av新在线| 一个人免费在线观看的高清视频| 18+在线观看网站| av视频在线观看入口| 国产精品久久久久久久电影| 一进一出抽搐动态| 婷婷精品国产亚洲av在线| 蜜桃亚洲精品一区二区三区| 色综合亚洲欧美另类图片| 免费人成在线观看视频色| 欧美一区二区亚洲| 国产黄片美女视频| 三级毛片av免费| 久久久久久久久久黄片| 国产精品国产高清国产av| 亚洲人成伊人成综合网2020| 日本精品一区二区三区蜜桃| 国产乱人伦免费视频| 国产69精品久久久久777片| 精品久久久久久久久av| 亚洲人成网站在线播放欧美日韩| 亚洲精品456在线播放app | 很黄的视频免费|