• <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)用
    e午夜精品久久久久久久| 欧美激情极品国产一区二区三区| 99国产精品一区二区蜜桃av| e午夜精品久久久久久久| 99国产精品99久久久久| 亚洲全国av大片| 国产一区二区三区视频了| 嫩草影视91久久| 亚洲久久久国产精品| 久久久久国内视频| 亚洲成a人片在线一区二区| 日日夜夜操网爽| 老司机在亚洲福利影院| 人人妻,人人澡人人爽秒播| 日韩欧美一区视频在线观看| 美女高潮喷水抽搐中文字幕| 午夜亚洲福利在线播放| cao死你这个sao货| 一级作爱视频免费观看| 欧美一级毛片孕妇| 免费无遮挡裸体视频| 日本五十路高清| 国产av一区二区精品久久| 午夜福利免费观看在线| 非洲黑人性xxxx精品又粗又长| 最新美女视频免费是黄的| 在线国产一区二区在线| 久久人妻av系列| 精品一区二区三区av网在线观看| 欧洲精品卡2卡3卡4卡5卡区| 91麻豆av在线| 亚洲精品国产一区二区精华液| 亚洲五月婷婷丁香| 亚洲精品粉嫩美女一区| 亚洲av成人av| 一区二区日韩欧美中文字幕| 国产在线观看jvid| 此物有八面人人有两片| 不卡一级毛片| 九色亚洲精品在线播放| 在线十欧美十亚洲十日本专区| 国产成人精品久久二区二区免费| 久久青草综合色| netflix在线观看网站| 老熟妇乱子伦视频在线观看| 少妇熟女aⅴ在线视频| 91精品国产国语对白视频| 法律面前人人平等表现在哪些方面| 在线av久久热| 别揉我奶头~嗯~啊~动态视频| 国产成人精品无人区| 久久久精品欧美日韩精品| 国产一区二区激情短视频| 国产精品久久久人人做人人爽| 精品国产一区二区三区四区第35| 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩一区二区三区在线| 在线观看www视频免费| 日韩免费av在线播放| 91字幕亚洲| 日韩欧美在线二视频| 欧美另类亚洲清纯唯美| x7x7x7水蜜桃| 欧美成狂野欧美在线观看| 一区二区三区精品91| 日韩高清综合在线| 99久久精品国产亚洲精品| 大码成人一级视频| 9色porny在线观看| 香蕉久久夜色| 99国产精品一区二区三区| 女同久久另类99精品国产91| 性少妇av在线| 99在线视频只有这里精品首页| 中文字幕色久视频| 亚洲男人的天堂狠狠| 久久青草综合色| 精品国产国语对白av| 99久久99久久久精品蜜桃| 天天一区二区日本电影三级 | 曰老女人黄片| 怎么达到女性高潮| 中文字幕精品免费在线观看视频| 精品不卡国产一区二区三区| 欧美性长视频在线观看| 久久久久久久久免费视频了| 伦理电影免费视频| 桃色一区二区三区在线观看| 精品熟女少妇八av免费久了| 桃红色精品国产亚洲av| 这个男人来自地球电影免费观看| 亚洲av熟女| 99国产精品一区二区三区| 身体一侧抽搐| 国产欧美日韩精品亚洲av| 在线观看免费午夜福利视频| 在线观看舔阴道视频| 亚洲色图综合在线观看| 午夜福利影视在线免费观看| 久久人人精品亚洲av| 国产亚洲精品久久久久5区| 狂野欧美激情性xxxx| 欧美精品啪啪一区二区三区| 久久精品国产清高在天天线| 大陆偷拍与自拍| av视频在线观看入口| 亚洲色图 男人天堂 中文字幕| 精品日产1卡2卡| 亚洲av五月六月丁香网| 一a级毛片在线观看| 国产三级黄色录像| 操美女的视频在线观看| 久久精品亚洲精品国产色婷小说| 男女做爰动态图高潮gif福利片 | 色综合欧美亚洲国产小说| www.熟女人妻精品国产| 久久久久国产精品人妻aⅴ院| 久久国产精品男人的天堂亚洲| 一夜夜www| 亚洲欧美精品综合久久99| 亚洲一区中文字幕在线| 日本精品一区二区三区蜜桃| 色精品久久人妻99蜜桃| 桃色一区二区三区在线观看| 欧美不卡视频在线免费观看 | 九色亚洲精品在线播放| 国产精品一区二区三区四区久久 | 日本免费a在线| 操美女的视频在线观看| 真人做人爱边吃奶动态| avwww免费| 国产精品久久视频播放| 窝窝影院91人妻| 最好的美女福利视频网| a在线观看视频网站| av天堂久久9| 纯流量卡能插随身wifi吗| 天天添夜夜摸| 国产欧美日韩一区二区精品| 看黄色毛片网站| 涩涩av久久男人的天堂| 91国产中文字幕| 精品久久久久久久人妻蜜臀av | 波多野结衣一区麻豆| 男男h啪啪无遮挡| 亚洲国产高清在线一区二区三 | 少妇裸体淫交视频免费看高清 | 男女下面插进去视频免费观看| 熟女少妇亚洲综合色aaa.| 亚洲欧美激情在线| 夜夜夜夜夜久久久久| 亚洲自偷自拍图片 自拍| 一a级毛片在线观看| 久久久久久国产a免费观看| 后天国语完整版免费观看| 俄罗斯特黄特色一大片| 午夜福利视频1000在线观看 | 美女国产高潮福利片在线看| 12—13女人毛片做爰片一| 韩国精品一区二区三区| 精品人妻在线不人妻| 狠狠狠狠99中文字幕| 这个男人来自地球电影免费观看| 国产精品久久电影中文字幕| 手机成人av网站| 巨乳人妻的诱惑在线观看| 黄频高清免费视频| 国产极品粉嫩免费观看在线| 制服人妻中文乱码| 国产1区2区3区精品| 亚洲一区二区三区色噜噜| 亚洲精品久久成人aⅴ小说| 日韩欧美国产一区二区入口| 国产私拍福利视频在线观看| or卡值多少钱| 亚洲五月婷婷丁香| 免费不卡黄色视频| 欧美大码av| 欧美最黄视频在线播放免费| 咕卡用的链子| 麻豆成人av在线观看| 乱人伦中国视频| 国产主播在线观看一区二区| 免费高清在线观看日韩| 久久久久九九精品影院| 亚洲欧美精品综合久久99| 看免费av毛片| 精品不卡国产一区二区三区| 丰满的人妻完整版| 日本 欧美在线| 成人18禁在线播放| 精品人妻1区二区| 69av精品久久久久久| 国产在线精品亚洲第一网站| 亚洲国产精品999在线| 十八禁人妻一区二区| 国产精品野战在线观看| 国产精品免费一区二区三区在线| 亚洲免费av在线视频| 一级a爱视频在线免费观看| 国产一区二区激情短视频| 少妇的丰满在线观看| 国产99久久九九免费精品| 每晚都被弄得嗷嗷叫到高潮| 亚洲av熟女| 窝窝影院91人妻| 欧美av亚洲av综合av国产av| 琪琪午夜伦伦电影理论片6080| aaaaa片日本免费| 国产精品国产高清国产av| 亚洲熟妇中文字幕五十中出| 欧洲精品卡2卡3卡4卡5卡区| 美国免费a级毛片| 日韩 欧美 亚洲 中文字幕| 久久午夜亚洲精品久久| 中文字幕高清在线视频| 最近最新免费中文字幕在线| 国产精品一区二区免费欧美| 亚洲第一av免费看| 国产熟女午夜一区二区三区| 欧美不卡视频在线免费观看 | 久久精品国产亚洲av香蕉五月| 操出白浆在线播放| 麻豆成人av在线观看| 天天添夜夜摸| 久久久久久久久免费视频了| 精品国产一区二区三区四区第35| 女人精品久久久久毛片| 欧美日韩黄片免| 搡老妇女老女人老熟妇| 看黄色毛片网站| 国产91精品成人一区二区三区| 亚洲一区高清亚洲精品| 亚洲av成人不卡在线观看播放网| 午夜福利视频1000在线观看 | 视频在线观看一区二区三区| 国产伦一二天堂av在线观看| 国产精品九九99| 欧美国产日韩亚洲一区| 国产精品综合久久久久久久免费 | 亚洲第一电影网av| 97碰自拍视频| 国产伦人伦偷精品视频| 精品欧美一区二区三区在线| 又紧又爽又黄一区二区| 国产精品九九99| 正在播放国产对白刺激| 午夜免费成人在线视频| 久久性视频一级片| 老司机深夜福利视频在线观看| 午夜福利一区二区在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 韩国av一区二区三区四区| 国产av一区二区精品久久| 桃色一区二区三区在线观看| 日韩精品青青久久久久久| www.www免费av| 国产成+人综合+亚洲专区| 亚洲熟妇中文字幕五十中出| 久久人妻av系列| 精品国产一区二区三区四区第35| 欧美中文日本在线观看视频| 久久精品91蜜桃| 国产精品自产拍在线观看55亚洲| 精品国内亚洲2022精品成人| 国产精品二区激情视频| 日本在线视频免费播放| 中文字幕人成人乱码亚洲影| 在线av久久热| 丝袜在线中文字幕| av在线天堂中文字幕| 性欧美人与动物交配| 免费在线观看视频国产中文字幕亚洲| 无遮挡黄片免费观看| 国产成年人精品一区二区| www.www免费av| svipshipincom国产片| 亚洲精品一卡2卡三卡4卡5卡| 国产精品香港三级国产av潘金莲| 亚洲一码二码三码区别大吗| 亚洲免费av在线视频| 丝袜美足系列| 午夜免费成人在线视频| 国产午夜精品久久久久久| 精品熟女少妇八av免费久了| 久久人妻av系列| 精品少妇一区二区三区视频日本电影| 在线天堂中文资源库| 两人在一起打扑克的视频| 老熟妇乱子伦视频在线观看| 亚洲av熟女| 精品久久久精品久久久| 午夜福利,免费看| 久久久精品国产亚洲av高清涩受| 国产成+人综合+亚洲专区| 亚洲av美国av| 精品一区二区三区四区五区乱码| 777久久人妻少妇嫩草av网站| 变态另类丝袜制服| 国产欧美日韩综合在线一区二区| 国产人伦9x9x在线观看| 成人特级黄色片久久久久久久| 久久婷婷人人爽人人干人人爱 | 久久 成人 亚洲| 老熟妇仑乱视频hdxx| 亚洲人成77777在线视频| 如日韩欧美国产精品一区二区三区| 久久影院123| 亚洲 国产 在线| 国产三级黄色录像| 国内精品久久久久久久电影| 在线国产一区二区在线| 亚洲少妇的诱惑av| 免费看美女性在线毛片视频| 国产精品久久电影中文字幕| 免费在线观看影片大全网站| 国产精品自产拍在线观看55亚洲| 成人18禁在线播放| ponron亚洲| 国产精品久久久久久人妻精品电影| 91老司机精品| 免费在线观看影片大全网站| 长腿黑丝高跟| 免费在线观看亚洲国产| 亚洲精品国产色婷婷电影| 亚洲国产日韩欧美精品在线观看 | 欧美一级毛片孕妇| 亚洲欧美日韩另类电影网站| 亚洲av成人av| 日本撒尿小便嘘嘘汇集6| 亚洲自偷自拍图片 自拍| 久久国产精品男人的天堂亚洲| 18美女黄网站色大片免费观看| 制服人妻中文乱码| 丝袜在线中文字幕| 亚洲性夜色夜夜综合| av片东京热男人的天堂| 一本大道久久a久久精品| 国产蜜桃级精品一区二区三区| 好男人在线观看高清免费视频 | 91在线观看av| 国产免费av片在线观看野外av| 热99re8久久精品国产| 一边摸一边做爽爽视频免费| 亚洲性夜色夜夜综合| 久久午夜综合久久蜜桃| 99久久久亚洲精品蜜臀av| 免费搜索国产男女视频| 国产成人免费无遮挡视频| 男女做爰动态图高潮gif福利片 | 丝袜美足系列| 久久精品影院6| 亚洲在线自拍视频| 免费高清视频大片| 一级毛片高清免费大全| 久久久久久大精品| 国产精品98久久久久久宅男小说| 在线天堂中文资源库| 精品国内亚洲2022精品成人| av有码第一页| 麻豆一二三区av精品| 亚洲中文av在线| 欧美最黄视频在线播放免费| 女性生殖器流出的白浆| 一级毛片高清免费大全| 国产精品自产拍在线观看55亚洲| 一本大道久久a久久精品| 少妇裸体淫交视频免费看高清 | 大型av网站在线播放| 久久久久国产精品人妻aⅴ院| 两个人视频免费观看高清| 中文字幕高清在线视频| 亚洲久久久国产精品| 亚洲男人的天堂狠狠| 精品国产一区二区三区四区第35| 多毛熟女@视频| 国产av一区二区精品久久| 在线播放国产精品三级| 男女午夜视频在线观看| 午夜免费激情av| 一边摸一边抽搐一进一出视频| 91av网站免费观看| 亚洲午夜精品一区,二区,三区| 国产亚洲精品久久久久5区| 国产99久久九九免费精品| 国产成年人精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 好男人电影高清在线观看| 亚洲成av人片免费观看| 悠悠久久av| 麻豆国产av国片精品| 久热爱精品视频在线9| 亚洲精品中文字幕一二三四区| 真人做人爱边吃奶动态| 深夜精品福利| 欧美中文综合在线视频| 午夜a级毛片| 黄色丝袜av网址大全| 欧美激情 高清一区二区三区| 欧美久久黑人一区二区| 99在线人妻在线中文字幕| 麻豆一二三区av精品| 国产精品自产拍在线观看55亚洲| 在线十欧美十亚洲十日本专区| 欧美成人性av电影在线观看| 精品福利观看| 99riav亚洲国产免费| 国产一级毛片七仙女欲春2 | 看片在线看免费视频| 免费高清在线观看日韩| 欧美日韩瑟瑟在线播放| 麻豆av在线久日| 动漫黄色视频在线观看| 精品熟女少妇八av免费久了| 国产亚洲精品第一综合不卡| 色av中文字幕| 自线自在国产av| 99香蕉大伊视频| 亚洲专区中文字幕在线| 一本久久中文字幕| 国产精品免费视频内射| 国产亚洲av嫩草精品影院| 女人被躁到高潮嗷嗷叫费观| 极品人妻少妇av视频| 村上凉子中文字幕在线| svipshipincom国产片| 日韩欧美三级三区| 国产激情欧美一区二区| 亚洲九九香蕉| 欧美不卡视频在线免费观看 | 亚洲av熟女| 变态另类成人亚洲欧美熟女 | 女人被躁到高潮嗷嗷叫费观| 给我免费播放毛片高清在线观看| 国产精品综合久久久久久久免费 | 午夜视频精品福利| 在线av久久热| 在线永久观看黄色视频| 日日夜夜操网爽| 黄色视频,在线免费观看| www日本在线高清视频| 1024视频免费在线观看| 一进一出抽搐gif免费好疼| 人人妻人人澡人人看| 91成人精品电影| 国产成人啪精品午夜网站| 亚洲片人在线观看| 女人精品久久久久毛片| 999久久久精品免费观看国产| 两个人视频免费观看高清| 久久亚洲精品不卡| 亚洲五月婷婷丁香| 757午夜福利合集在线观看| 黄色视频不卡| 国产av一区在线观看免费| 亚洲国产毛片av蜜桃av| 一级作爱视频免费观看| 亚洲熟女毛片儿| 夜夜爽天天搞| 国产伦一二天堂av在线观看| 国产精品久久久久久人妻精品电影| 久热爱精品视频在线9| 久久午夜亚洲精品久久| 久久久久亚洲av毛片大全| 国产麻豆成人av免费视频| 亚洲精品在线观看二区| 成熟少妇高潮喷水视频| 人人妻人人爽人人添夜夜欢视频| 久久青草综合色| 90打野战视频偷拍视频| 亚洲自偷自拍图片 自拍| 中亚洲国语对白在线视频| 亚洲人成网站在线播放欧美日韩| 别揉我奶头~嗯~啊~动态视频| 久久欧美精品欧美久久欧美| 精品日产1卡2卡| av天堂在线播放| 久99久视频精品免费| 丝袜在线中文字幕| 久久精品国产99精品国产亚洲性色 | 少妇粗大呻吟视频| 久久性视频一级片| 不卡一级毛片| 两人在一起打扑克的视频| av在线播放免费不卡| 人人澡人人妻人| 午夜福利免费观看在线| 亚洲午夜精品一区,二区,三区| 成年人黄色毛片网站| 最好的美女福利视频网| 亚洲欧洲精品一区二区精品久久久| 亚洲熟妇熟女久久| 搞女人的毛片| 久久九九热精品免费| 久久人人97超碰香蕉20202| 婷婷丁香在线五月| 男女床上黄色一级片免费看| 国产亚洲精品一区二区www| cao死你这个sao货| 国语自产精品视频在线第100页| 免费搜索国产男女视频| 久久精品综合一区二区三区| 国产探花极品一区二区| 国产亚洲av嫩草精品影院| av中文乱码字幕在线| 国产精品永久免费网站| 国产精品三级大全| 亚洲av五月六月丁香网| 久99久视频精品免费| 日韩中字成人| 午夜精品一区二区三区免费看| 超碰av人人做人人爽久久| 国产精品亚洲一级av第二区| aaaaa片日本免费| 丰满的人妻完整版| netflix在线观看网站| videossex国产| 男人舔奶头视频| 成人二区视频| 免费不卡的大黄色大毛片视频在线观看 | 小蜜桃在线观看免费完整版高清| 日韩av在线大香蕉| 22中文网久久字幕| 日本在线视频免费播放| 丰满乱子伦码专区| 一级黄色大片毛片| 国产日本99.免费观看| 九九久久精品国产亚洲av麻豆| 国产在线男女| 欧美成人一区二区免费高清观看| 人人妻,人人澡人人爽秒播| 尤物成人国产欧美一区二区三区| 舔av片在线| 欧美区成人在线视频| 午夜激情福利司机影院| 久久久久久久久大av| 在线观看舔阴道视频| 高清毛片免费观看视频网站| 日本与韩国留学比较| 国产高清三级在线| 精品免费久久久久久久清纯| 精品一区二区三区人妻视频| 成人亚洲精品av一区二区| 麻豆国产97在线/欧美| 亚洲色图av天堂| 国产成人av教育| 男女做爰动态图高潮gif福利片| 赤兔流量卡办理| 天天躁日日操中文字幕| 成人国产一区最新在线观看| 久久久久久久久久久丰满 | 中文字幕久久专区| 欧美在线一区亚洲| 国内少妇人妻偷人精品xxx网站| 老司机福利观看| 国产一区二区亚洲精品在线观看| av福利片在线观看| 麻豆av噜噜一区二区三区| 美女黄网站色视频| 不卡视频在线观看欧美| 精品久久久久久久末码| 伊人久久精品亚洲午夜| 午夜爱爱视频在线播放| 九九久久精品国产亚洲av麻豆| 亚洲精华国产精华液的使用体验 | 永久网站在线| 啪啪无遮挡十八禁网站| 国产精品99久久久久久久久| 男女之事视频高清在线观看| 看免费成人av毛片| 日韩高清综合在线| 日韩欧美 国产精品| 在线看三级毛片| 五月伊人婷婷丁香| 日本在线视频免费播放| 国产精品一区二区三区四区久久| 大又大粗又爽又黄少妇毛片口| 天堂√8在线中文| 久久久久久久亚洲中文字幕| 欧美精品国产亚洲| 国产伦精品一区二区三区四那| 乱人视频在线观看| 一个人看的www免费观看视频| 亚洲最大成人手机在线| 欧美又色又爽又黄视频| 91麻豆av在线| 草草在线视频免费看| 国产aⅴ精品一区二区三区波| 美女高潮喷水抽搐中文字幕| 日韩大尺度精品在线看网址| 少妇的逼好多水| 小蜜桃在线观看免费完整版高清| 久久久精品欧美日韩精品| 欧美+日韩+精品| 日日干狠狠操夜夜爽| 白带黄色成豆腐渣| 日本一二三区视频观看| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久久久免| 国产私拍福利视频在线观看| 波多野结衣巨乳人妻| 一级毛片久久久久久久久女| 欧美日韩瑟瑟在线播放| 久久午夜福利片| 动漫黄色视频在线观看| 久99久视频精品免费| 日本与韩国留学比较| 九九爱精品视频在线观看| 99久久中文字幕三级久久日本| 又黄又爽又刺激的免费视频.| 亚洲精品粉嫩美女一区| 欧美+日韩+精品| 亚洲精华国产精华精| 欧美三级亚洲精品| 国产日本99.免费观看|