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

    基于螺旋纏繞光纖的分布式聲波傳感三分量信號(hào)采集機(jī)理研究

    2022-01-28 06:27:26曹丹平律軍軍孫上饒馬國(guó)旗尹教建
    石油物探 2022年1期
    關(guān)鍵詞:軸向坐標(biāo)系分量

    曹丹平,律軍軍,孫上饒,馬國(guó)旗,尹教建

    (1.中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島266580;2.中國(guó)石油大學(xué)(華東)山東省深層油氣重點(diǎn)實(shí)驗(yàn)室,山東青島266580;3.哈利法大學(xué)文理學(xué)院,阿布扎比127788;4.中國(guó)石油大學(xué)(華東)理學(xué)院,山東青島266580)

    分布式光纖聲波傳感(DAS)技術(shù)是一種以光纖本身作為傳感器,基于光的瑞利散射實(shí)現(xiàn)振動(dòng)信號(hào)采集的新型光纖傳感技術(shù),具有成本低、測(cè)量精度高、不易受電磁干擾和安裝簡(jiǎn)便等優(yōu)點(diǎn)[1]。相比于常規(guī)單點(diǎn)式和準(zhǔn)分布式傳感器,DAS更適合用于長(zhǎng)距離或時(shí)間及空間分辨率要求較高的應(yīng)用場(chǎng)景,并在石油勘探、管道泄漏監(jiān)測(cè)以及邊防安全監(jiān)測(cè)等方面被廣泛應(yīng)用。

    近年來(lái),研究人員致力于用DAS替代常規(guī)檢波器進(jìn)行微地震監(jiān)測(cè)、淺地表地震勘探以及井中地震勘探等,并在基礎(chǔ)理論、采集、處理和解釋等領(lǐng)域取得了長(zhǎng)足的進(jìn)步[2-8]。但分布式聲波傳感的原理決定了DAS對(duì)光纖軸向應(yīng)變敏感,而對(duì)徑向應(yīng)變不敏感,即直光纖DAS只能采集到沿光纖軸向方向的單分量振動(dòng)信號(hào)[9-10]。隨著地震勘探技術(shù)的進(jìn)步,油氣工業(yè)對(duì)三分量地震信號(hào)的采集需求越來(lái)越高,直光纖DAS只能采集單分量信號(hào)的特點(diǎn),無(wú)法滿足油氣地球物理行業(yè)的地震資料采集需求。根據(jù)DAS光纖對(duì)軸向振動(dòng)敏感的特點(diǎn),如果將光纖的軸向方向分別調(diào)整到振動(dòng)所在位置處的觀測(cè)坐標(biāo)系X,Y,Z3個(gè)正交方向,即可獲得對(duì)應(yīng)位置處的三分量振動(dòng)信號(hào),也就是光纖檢波器的原理[11],但光纖檢波器屬于傳統(tǒng)的單點(diǎn)探測(cè),缺乏分布式傳感技術(shù)能夠開(kāi)展連續(xù)采集的優(yōu)勢(shì)。受光纖檢波器和分布式聲傳感技術(shù)的啟發(fā),在實(shí)際應(yīng)用中,可以通過(guò)對(duì)光纖進(jìn)行纏繞布置獲得更加豐富的光纖軸向方向組合,即螺旋纏繞分布式光纖聲波傳感(HWC-DAS),從而在確保獲得X,Y,Z3個(gè)方向振動(dòng)信號(hào)的同時(shí)實(shí)現(xiàn)信息的連續(xù)采集。

    HWC-DAS技術(shù)增強(qiáng)了DAS對(duì)不同方向分量振動(dòng)信號(hào)的敏感性,可以采集到多分量振動(dòng)信號(hào),實(shí)現(xiàn)了對(duì)直光纖DAS技術(shù)的拓展,并推動(dòng)了DAS技術(shù)在多分量地震勘探中的應(yīng)用。在螺旋纏繞光纖多分量信號(hào)采集的基礎(chǔ)理論研究方面,KUVSHINOV[12]對(duì)比了HWC-DAS中的光纖纏繞形狀對(duì)DAS響應(yīng)的影響,并分析了光纖與井壁在完全耦合或不耦合情況下體波、面波對(duì)DAS的角度敏感特征。NING等[13-14]先后提出利用沿單根和多根光纖測(cè)量的6個(gè)應(yīng)變投影重建全應(yīng)變張量的多分量DAS數(shù)據(jù)重構(gòu)方法,討論了不同標(biāo)距長(zhǎng)度對(duì)多分量應(yīng)變信息恢復(fù)的影響,并通過(guò)數(shù)值模擬驗(yàn)證了該方法的有效性。INNANEN[15-16]在應(yīng)變張量估計(jì)的背景下,提出一種沿任意彎曲的纏繞軸纏繞的螺旋光纖模型,并討論了在多分量應(yīng)變張量重建過(guò)程中標(biāo)距長(zhǎng)度和重建間隔對(duì)P波和S波信號(hào)重建的影響,該模型也被用于地震彈性波的正、反演[17-20]。WUESTEFELD等[21]在復(fù)雜速度模型中使用精確的射線追蹤方法確定入射角,并進(jìn)一步討論了介質(zhì)物性和DAS光纜包裹材料與最佳纏繞角之間的關(guān)系。

    同時(shí),在利用螺旋光纖開(kāi)展多分量振動(dòng)信號(hào)的采集與應(yīng)用方面也得到了快速發(fā)展。KARRENBACH等[22]進(jìn)行了井中三分量DAS數(shù)據(jù)的采集,并與常規(guī)三分量采集的數(shù)據(jù)進(jìn)行了對(duì)比,表明了多分量DAS技術(shù)的應(yīng)用潛力。HORNMAN等[23-24]利用螺旋纏繞光纖進(jìn)行地震波數(shù)據(jù)采集,研究了HWC-DAS的角度敏感性,同時(shí),對(duì)比了垂直分量檢波器、HWC-DAS以及直光纖DAS采集的實(shí)際數(shù)據(jù)之間的差異,表明HWC-DAS對(duì)縱橫波和面波均有較高的敏感性。在實(shí)際應(yīng)用中,DAS采集的是應(yīng)變數(shù)據(jù),與常規(guī)動(dòng)圈式或壓電式檢波器采集的數(shù)據(jù)既有聯(lián)系也有差異,在地震資料處理過(guò)程中,為了實(shí)現(xiàn)同類型地震數(shù)據(jù)對(duì)比分析,可以將DAS采集到的應(yīng)變信號(hào)或應(yīng)變率信號(hào)轉(zhuǎn)變?yōu)樗俣刃盘?hào)[7,19]或加速度信號(hào)。YURIKOV等[25]利用HWC-DAS對(duì)地震頻帶內(nèi)的巖石物理屬性進(jìn)行測(cè)量,實(shí)驗(yàn)結(jié)果與傳統(tǒng)儀器的測(cè)量結(jié)果吻合度較高,表明了DAS在巖石物理實(shí)驗(yàn)中的巨大潛力。陳亞峰等[26]采用螺旋纏繞光纖技術(shù)研究了熱力管道泄漏監(jiān)測(cè)系統(tǒng),表明螺旋纏繞光纖技術(shù)有助于提高管線泄漏定位精度,具有良好的檢測(cè)效果。

    光纖的軸向敏感性導(dǎo)致光纖對(duì)不同類型的地震波具有不同的響應(yīng)機(jī)理,從地震波的傳播特征出發(fā)研究地震波入射角、速度特征和介質(zhì)物性對(duì)光纖軸向應(yīng)變的影響有助于更好地理解DAS在地震資料采集中的作用[12,27-29]。由于縱波和橫波耦合在一起,在常規(guī)三維地震勘探中,采集到的是X,Y,Z3個(gè)相互正交的振動(dòng)信號(hào),需要通過(guò)對(duì)3個(gè)分量信號(hào)進(jìn)行處理間接獲得縱橫波信息,因此,討論螺旋光纖軸向應(yīng)變與X,Y,Z三分量振動(dòng)信號(hào)之間的關(guān)系有助于更好地與常規(guī)檢波器采集到的資料進(jìn)行對(duì)比。螺旋光纖采集到的振動(dòng)信號(hào)相當(dāng)于實(shí)際X,Y,Z3個(gè)方向的振動(dòng)信號(hào)在光纖軸向方向上的投影,厘清光纖采集信號(hào)與實(shí)際三分量振動(dòng)信號(hào)之間的理論關(guān)系對(duì)分布式光纖三分量地震信號(hào)采集具有重要的理論意義和指導(dǎo)作用。

    本文聚焦光纖纏繞方式與三分量振動(dòng)信號(hào)之間的關(guān)系開(kāi)展基礎(chǔ)研究,在假設(shè)螺旋光纖與纏繞軸完全耦合的前提下,首先,推導(dǎo)了螺旋纏繞光纖軸向應(yīng)變與三分量振動(dòng)信號(hào)之間的映射關(guān)系。然后,分析了不同條件下三分量振動(dòng)信號(hào)對(duì)光纖軸向應(yīng)變的影響。最后,采用波動(dòng)方程模擬的三分量數(shù)字信號(hào)對(duì)螺旋光纖軸向應(yīng)變所接收到的等效信號(hào)進(jìn)行驗(yàn)證,并進(jìn)行敏感性分析。

    1 螺旋纏繞光纖三分量信號(hào)采集理論模型

    1.1 螺旋纏繞光纖軸向應(yīng)變的影響因素

    圖1 DAS螺旋光纖纏繞幾何模型

    根據(jù)圖1b所示幾何關(guān)系,可以得到:

    L2=C2+H2

    (1)

    對(duì)公式兩邊同時(shí)進(jìn)行微分得到:

    2L·dL=2C·dC+2H·dH

    (2)

    L·ΔL=C·ΔC+H·ΔH

    (3)

    等式兩邊同時(shí)除以L2,整理可得:

    (4)

    由于纏繞軸圓周周長(zhǎng)C=2πr,ΔC/C=Δr/r,所以公式(4)可以表示為:

    (5)

    為了討論方便,在完全各向同性介質(zhì)中,可以采用從A點(diǎn)纏繞至A′點(diǎn)的整體應(yīng)變來(lái)代表P點(diǎn)處的應(yīng)變,即:

    (6)

    此時(shí),螺旋光纖的軸向應(yīng)變可表示為:

    (7)

    (8)

    由公式(8)可定性地獲得螺旋光纖軸向應(yīng)變與真實(shí)地震波振動(dòng)引起的實(shí)際三分量應(yīng)變之間的關(guān)系,即當(dāng)光纖纏繞角度α為固定值時(shí),ezsin2α為定值,也就是說(shuō)當(dāng)DAS光纖以固定角度均勻繞軸纏繞時(shí),螺旋光纖采集到的軸向應(yīng)變中ez分量的絕對(duì)數(shù)值始終保持不變,而該分量在最終采集到的軸向應(yīng)變總數(shù)值中所占比重則隨著光纖纏繞角度的變化而變化;當(dāng)光纖纏繞角度α增大時(shí),cos2α的值逐漸減小,sin2α的值逐漸增大,螺旋纏繞光纖采集到的軸向應(yīng)變信號(hào)中ez分量的占比增大,而ex和ey分量的占比則減小;即纏繞角越大,ez分量對(duì)螺旋纏繞光纖最終采集到的軸向應(yīng)變信號(hào)影響越大,纏繞角越小時(shí),則ex和ey分量對(duì)螺旋纏繞光纖最終采集到的軸向應(yīng)變信號(hào)影響越大。

    從圖1中還可以看出,在一個(gè)纏繞周期內(nèi),光纖接收到的軸向應(yīng)變還與光纖所在的空間位置有關(guān)。由于纏繞角α可以由h,l和c中的任意兩個(gè)參數(shù)確定,也可以將螺旋纏繞光纖看作三維空間坐標(biāo)系中由h,r和θ3個(gè)參數(shù)決定的螺旋曲線X(h,r,α),表示為:

    X(h,r,α)=(h,rcosθ,rsinθ)

    (9)

    式中:θ為旋轉(zhuǎn)角,該參數(shù)與h,r和α有關(guān)。

    光纖在繞軸纏繞前行的過(guò)程中導(dǎo)致旋轉(zhuǎn)角產(chǎn)生周期性的規(guī)律變化,為了探究θ與h,r和α之間的具體關(guān)系,由圖1中一個(gè)纏繞周期內(nèi)光纖上任意點(diǎn)P的位置可知:

    (10)

    (11)

    由于圖1b為圖1a的平面展開(kāi),因此有:

    (12)

    根據(jù)圖1b中幾何關(guān)系有:

    |AM|=|AP|cosα

    |PM|=|AP|sinα

    (13)

    整理得:

    rθ=l·cosα

    h=l·sinα

    (14)

    公式(14)中的旋轉(zhuǎn)角θ可以表示為:

    (15)

    將公式(15)帶入公式(9),可以得到由h,r和α表示的螺旋纏繞光纖上任意點(diǎn)的三維空間位置關(guān)系公式:

    (16)

    式中:α∈(0,90°)。

    由公式(16)可知,螺旋光纖纏繞后任意點(diǎn)P的位置可以由纏繞軸半徑r、光纖纏繞角α和該點(diǎn)對(duì)應(yīng)的纏繞軸高度h唯一確定,當(dāng)然,纏繞軸高度h也可以用該點(diǎn)的光纖纏繞長(zhǎng)度l代表,兩個(gè)參數(shù)可以通過(guò)纏繞角直接轉(zhuǎn)換,由于在每一個(gè)纏繞周期內(nèi)纏繞角都是從0逐步增加到2π,因此后續(xù)的討論將在一個(gè)周期范圍內(nèi)開(kāi)展。

    1.2 螺旋光纖軸向應(yīng)變與三分量振動(dòng)的映射關(guān)系

    為了探究螺旋光纖軸向應(yīng)變與三分量振動(dòng)之間的映射關(guān)系,根據(jù)DAS螺旋纏繞光纖示意(圖1),建立了如圖2a所示的DAS螺旋光纖局部坐標(biāo)系相對(duì)關(guān)系示意,圖2b為光纖從A點(diǎn)纏繞到任意點(diǎn)P時(shí)所對(duì)應(yīng)旋轉(zhuǎn)角θ的示意。為了方便對(duì)比,在圖中標(biāo)注了A,A1,A2,A3,A′等5個(gè)特殊點(diǎn),這5個(gè)點(diǎn)分別表示旋轉(zhuǎn)角θ為0,π/2,π,3π/2和2π時(shí)所對(duì)應(yīng)光纖上的點(diǎn),通過(guò)這5個(gè)點(diǎn)有助于更加清晰地呈現(xiàn)光纖軸向應(yīng)變?cè)谝粋€(gè)纏繞周期內(nèi)的變化規(guī)律。首先定義地震波振動(dòng)引起的實(shí)際三分量應(yīng)變ex,ey,ez所在方向分別代表全局坐標(biāo)系A(chǔ)-xyz中x軸、y軸和z軸的坐標(biāo)方向??紤]到光纖對(duì)軸向應(yīng)變敏感,為了便于討論,沿光纖的軸向方向來(lái)定義局部坐標(biāo)系,A-lAmAnA表示螺旋光纖上A點(diǎn)的局部坐標(biāo)系,P-lmn表示螺旋光纖上P點(diǎn)的局部坐標(biāo)系。其中,nA和n方向始終沿著光纖的軸向方向。

    圖2 DAS螺旋光纖坐標(biāo)系相對(duì)關(guān)系示意(a)和P點(diǎn)所對(duì)應(yīng)的纏繞角θ示意(b)

    光纖上P點(diǎn)處的局部坐標(biāo)系P-lmn可以看作是全局坐標(biāo)系A(chǔ)-xyz通過(guò)坐標(biāo)系旋轉(zhuǎn)得到的。整個(gè)坐標(biāo)旋轉(zhuǎn)的過(guò)程可以分解為兩步:①先沿著y軸正方向,將全局坐標(biāo)系A(chǔ)-xyz的A-xz平面繞y軸順時(shí)針旋轉(zhuǎn)π/2-α,獲得纏繞光纖在A點(diǎn)處的局部坐標(biāo)系A(chǔ)-lAmAnA;②再沿著z軸正方向,將坐標(biāo)系A(chǔ)-lAmAnA繞z軸順時(shí)針旋轉(zhuǎn)θ,即可獲得P點(diǎn)處的局部坐標(biāo)系P-lmn。

    由于光纖對(duì)軸向應(yīng)變敏感,因此,螺旋光纖上任意點(diǎn)處所測(cè)得的軸向應(yīng)變實(shí)際上是地震波振動(dòng)引起的X,Y,Z三分量應(yīng)變?cè)诼菪饫w所在點(diǎn)的軸向方向上的投影。根據(jù)本文定義,由地震波振動(dòng)引起的實(shí)際三分量應(yīng)變的方向和光纖軸向應(yīng)變的方向分別為全局坐標(biāo)系中3個(gè)坐標(biāo)軸的方向和局部坐標(biāo)系中n軸的方向,因此,求解該投影關(guān)系的關(guān)鍵在于準(zhǔn)確確定局部坐標(biāo)系中的n軸與全局坐標(biāo)系中x軸,y軸和z軸之間的位置關(guān)系(圖3)。

    圖3a展示了圖2a中光纖軌跡上任意點(diǎn)所處局部坐標(biāo)系中n軸繞z軸旋轉(zhuǎn)的軌跡,其中A,A1,A2和A3分別表示旋轉(zhuǎn)角θ為0,π/2,π和3π/2時(shí),n軸所處位置點(diǎn),根據(jù)圖3中標(biāo)注的關(guān)系可知,AA2平行于x軸,A1A3平行于y軸,nA表示螺旋光纖纏繞起點(diǎn)A的軸向方向,n表示螺旋光纖上任意點(diǎn)P的軸向方向,紅實(shí)線箭頭表示局部坐標(biāo)系中n軸繞z軸旋轉(zhuǎn)的運(yùn)動(dòng)方向;圖3b中θnx,θny和θnz分別表示局部坐標(biāo)系中的n軸分別與全局坐標(biāo)系中的x軸、y軸和z軸之間的夾角。

    圖3 局部坐標(biāo)系中n軸繞z軸旋轉(zhuǎn)軌跡示意(a)和n軸與全局坐標(biāo)系坐標(biāo)軸之間的夾角關(guān)系(b)

    由圖3a中局部坐標(biāo)系中n軸繞z軸的旋轉(zhuǎn)軌跡可知,局部坐標(biāo)系中的n軸與全局坐標(biāo)系z(mì)軸的夾角θnz不隨旋轉(zhuǎn)角θ的變化而變化,且僅與光纖的纏繞角有關(guān),即

    (17)

    在局部坐標(biāo)系n軸繞z軸旋轉(zhuǎn)的同時(shí),n軸分別與x軸和y軸的夾角θnx和θny隨著旋轉(zhuǎn)角θ以正弦(余弦)規(guī)律呈現(xiàn)出連續(xù)的均勻變化趨勢(shì),為了建立該變化規(guī)律,表1列出了光纖中任意點(diǎn)P纏繞至A,A1,A2,A3和A′點(diǎn)時(shí),旋轉(zhuǎn)角θ和夾角θnx之間的對(duì)應(yīng)關(guān)系。

    根據(jù)表1呈現(xiàn)出的規(guī)律性角度關(guān)系,可以采用三

    表1 螺旋光纖旋轉(zhuǎn)θ后旋轉(zhuǎn)角與θnx之間的角度對(duì)應(yīng)關(guān)系

    角函數(shù)重建出如下的θnx隨纏繞角α和旋轉(zhuǎn)角θ之間的變化關(guān)系表達(dá)式:

    (18)

    同理可得:

    (19)

    將公式(15)帶入公式(18)和公式(19)可得:

    (20)

    在獲得局部坐標(biāo)系n軸與全局坐標(biāo)系坐標(biāo)軸之間的夾角之后,可以構(gòu)建全局坐標(biāo)系A(chǔ)-xyz與局部坐標(biāo)系P-lmn之間的旋轉(zhuǎn)映射關(guān)系,采用上述所得的n軸與x軸,y軸和z軸之間的方位角關(guān)系,可以將坐標(biāo)旋轉(zhuǎn)矩陣表示為[30]:

    (21)

    式中:Rlx,Rly,Rlz,Rmx,Rmy,Rmz,Rnx,Rny和Rnz分別表示l軸,m軸,n軸與x軸,y軸和z軸之間的方向余弦。其中,影響光纖軸向應(yīng)變的3個(gè)系數(shù)為:

    Rnx=cosθnx

    Rny=cosθny

    Rnz=cosθnz

    (22)

    通過(guò)矩陣乘法,可以將全局坐標(biāo)系A(chǔ)-xyz與局部坐標(biāo)系P-lmn之間的映射關(guān)系表示為:

    (23)

    光纖上任意點(diǎn)的軸向應(yīng)變與光纖所在位置處由地震波振動(dòng)引起的實(shí)際三分量應(yīng)變之間的關(guān)系可以表示為[30-31]:

    e=Re0RT

    (24)

    式中:e表示局部坐標(biāo)系中對(duì)應(yīng)的應(yīng)變張量;e0表示全局坐標(biāo)系中對(duì)應(yīng)的應(yīng)變張量;R表示坐標(biāo)旋轉(zhuǎn)矩陣;RT表示坐標(biāo)旋轉(zhuǎn)矩陣R的轉(zhuǎn)置。

    (25)

    整理公式(25)可以得到螺旋光纖軸向應(yīng)變與實(shí)際三分量應(yīng)變之間的映射關(guān)系公式:

    (26)

    公式(26)也可以簡(jiǎn)寫為:

    m=Gm0

    (27)

    式中:m表示螺旋光纖軸向應(yīng)變;G代表由旋轉(zhuǎn)矩陣R得到的擴(kuò)展矩陣;m0表示光纖所在位置處的實(shí)際三分量應(yīng)變。

    公式(27)代表了從實(shí)際三分量信號(hào)到螺旋光纖軸向應(yīng)變之間的投影過(guò)程,反過(guò)來(lái),在實(shí)際應(yīng)用中就可以使用采集到的光纖軸向應(yīng)變m重構(gòu)光纖所在位置處的真實(shí)三分量應(yīng)變m0,且理論上可以通過(guò)求取最小二乘解得到:

    m0=(GTG)-1GTm

    (28)

    因此,在實(shí)際應(yīng)用中,可以通過(guò)公式(28)來(lái)實(shí)現(xiàn)由螺旋光纖采集到的軸向應(yīng)變重構(gòu)出光纖所在位置處的實(shí)際地震波振動(dòng)的三分量應(yīng)變,此即利用螺旋纏繞光纖開(kāi)展三分量信號(hào)恢復(fù)的理論依據(jù)。

    2 螺旋纏繞光纖軸向應(yīng)變對(duì)三分量應(yīng)變的系數(shù)敏感性分析

    由圖4可知,在纏繞角度α一定的條件下,螺旋光纖采集到的軸向應(yīng)變中ez分量的值始終保持不變,ex和ey分量隨著旋轉(zhuǎn)角度的變化分別呈余弦和正弦的周期性變化,且無(wú)論在多大的纏繞角度下,兩個(gè)分量隨旋轉(zhuǎn)角的正余弦的變化規(guī)律不會(huì)改變,但幅度會(huì)隨著纏繞角度的增大而減小。隨著纏繞角度α增大,螺旋光纖采集到的軸向應(yīng)變中ez分量的值增大,但在一個(gè)旋轉(zhuǎn)周期內(nèi)的值均相同,而ex和ey分量的值逐漸減小,其值隨旋轉(zhuǎn)角的周期性變化規(guī)律保持不變。根據(jù)該變化規(guī)律可知,當(dāng)纏繞角度較小時(shí),光纖軸向應(yīng)變中ez分量的影響較小,而ex和ey分量的影響較大,特別是當(dāng)纏繞角度為0時(shí),光纖軸向應(yīng)變中采集到的ez分量為0,光纖軸向應(yīng)變?nèi)坑蒭x和ey分量來(lái)確定,表明此時(shí)DAS光纖只能采集到徑向應(yīng)變信息;當(dāng)纏繞角度較大時(shí),光纖軸向應(yīng)變中ez分量的影響增大,而ex和ey分量的影響減小,特別是當(dāng)纏繞角度為90°時(shí),此時(shí)螺旋纏繞光纖為直光纖,光纖軸向應(yīng)變中只存在ez分量,無(wú)ex和ey分量。由此可見(jiàn),隨著纏繞角度α的增大,DAS螺旋光纖對(duì)地震波振動(dòng)應(yīng)變的ez分量敏感性逐漸增強(qiáng),對(duì)ex和ey分量的敏感性逐漸減弱。

    圖4 纏繞角α分別為5°(a),20°(b),35°(c),50°(d),65°(e)和80°(f)時(shí),螺旋光纖采集到的軸向應(yīng)變中ex,ey和ez 3個(gè)分量的比重變化趨勢(shì)

    3 基于三分量振動(dòng)信號(hào)的螺旋光纖軸向應(yīng)變敏感性數(shù)值分析

    為了進(jìn)一步驗(yàn)證DAS螺旋光纖對(duì)地震波振動(dòng)產(chǎn)生的三分量應(yīng)變的敏感性變化規(guī)律,采用各向同性介質(zhì)三維波動(dòng)方程數(shù)值模擬得到的振動(dòng)信號(hào)開(kāi)展數(shù)值測(cè)試。為了便于對(duì)比和分析,選擇了一個(gè)相位一致但振幅存在差異的三分量信號(hào)進(jìn)行對(duì)比。其中,X,Y,Z3個(gè)分量的振動(dòng)信號(hào)如圖5a所示,3個(gè)分量所對(duì)應(yīng)的振幅譜如圖5b所示。

    通過(guò)將X,Y,Z3個(gè)分量的振動(dòng)信號(hào)投影到螺旋光纖的軸向應(yīng)變方向上,即可獲得螺旋光纖軸向方向上的應(yīng)變特征,圖6a對(duì)比了纏繞角度α分別為5°,35°和80°時(shí)螺旋光纖上采集到的軸向應(yīng)變,圖6b對(duì)比了3個(gè)不同纏繞角度下所采集到的信號(hào)振幅譜。

    由圖6a可知,當(dāng)纏繞角度為5°時(shí),螺旋纏繞光纖軸向應(yīng)變振幅最大值明顯低于ez分量應(yīng)變的振幅最大值,但與ex和ey分量振幅最大值接近(類似圖5a),且在圖6b的振幅譜也具有相同的特點(diǎn),即纏繞角較小時(shí),螺旋纏繞光纖軸向方向上采集到的信號(hào)主要以X和Y分量成分為主,Z分量的影響較小。當(dāng)纏繞角度增大到35°時(shí),螺旋光纖軸向方向上獲得的信號(hào)幅值和振幅譜特征開(kāi)始逐漸接近ez分量的信號(hào)和振幅譜特征。而當(dāng)纏繞角度增大到80°時(shí),螺旋光纖軸向方向上獲得的信號(hào)幅值十分接近ez分量的幅值,同時(shí)在振幅譜上也具有非常明顯的接近趨勢(shì),即當(dāng)纏繞角較大時(shí),螺旋光纖軸向方向采集到的信號(hào)以Z分量信號(hào)為主。由此可見(jiàn),隨著纏繞角度α增大,螺旋光纖軸向方向采集到的信號(hào)中,Z分量的成分逐漸增多,而X和Y分量的成分逐漸減少。

    圖5 DAS光纖所在位置處由地震波振動(dòng)產(chǎn)生的三分量應(yīng)變信號(hào)(a)及對(duì)應(yīng)的振幅譜(b)

    圖6 不同纏繞角度時(shí)螺旋光纖軸向方向上采集到的信號(hào)(a)及對(duì)應(yīng)振幅譜(b)

    由此可見(jiàn),X和Y分量在纏繞角較小的時(shí)候?qū)β菪饫w軸向方向的應(yīng)變影響較大,由于X和Y兩個(gè)分量隨旋轉(zhuǎn)角具有正余弦變化規(guī)律,為了進(jìn)一步厘清X分量和Y分量對(duì)螺旋光纖軸向方向應(yīng)變特征的影響,在小纏繞角條件下分別對(duì)兩個(gè)特殊旋轉(zhuǎn)角和π的位置開(kāi)展數(shù)值測(cè)試。圖7a展示了纏繞角度α為5°,旋轉(zhuǎn)角度θ分別為π/2和π時(shí)螺旋光纖軸向應(yīng)變上所采集到的應(yīng)變,圖7b展示了對(duì)應(yīng)的振幅譜。

    圖7 相同纏繞角度不同旋轉(zhuǎn)角條件下螺旋光纖軸向方向上采集到的信號(hào)(a)及對(duì)應(yīng)的振幅譜(b)

    由圖7a可知,當(dāng)旋轉(zhuǎn)角度θ從π/2變化到π時(shí),螺旋光纖軸向方向上采集到的應(yīng)變幅值減小,同時(shí)振幅譜幅值也相應(yīng)減小。依據(jù)圖5可知,本例中由地震波振動(dòng)產(chǎn)生的Y分量信號(hào)的幅值和振幅譜幅值均大于X分量信號(hào),表明當(dāng)旋轉(zhuǎn)角度θ從π/2變化到π時(shí),螺旋光纖軸向方向上獲得的信號(hào)中X分量信號(hào)的影響增大,而Y分量信號(hào)的影響減小。同理當(dāng)旋轉(zhuǎn)角度θ從π增大到3π/2時(shí),螺旋光纖軸向方向上獲得的信號(hào)中X分量信號(hào)的影響減小,而Y分量信號(hào)的影響增大。由此可見(jiàn),當(dāng)纏繞角固定時(shí),在一個(gè)纏繞周期內(nèi)螺旋光纖軸向方向上采集的信號(hào)中的X分量信號(hào)的影響增大時(shí),Y分量信號(hào)的影響將減小;而X分量信號(hào)的影響減小時(shí),Y分量信號(hào)的影響將增大,兩個(gè)分量信號(hào)對(duì)螺旋光纖軸向方向上采集到的信號(hào)特征具有一定的周期性和互補(bǔ)性。

    4 結(jié)論

    本文通過(guò)建立單根光纖螺旋纏繞的數(shù)學(xué)模型,推導(dǎo)了螺旋光纖軸向應(yīng)變的表達(dá)式和螺旋光纖上任意點(diǎn)的三維空間位置關(guān)系式,明確了螺旋光纖軸向應(yīng)變與光纖所在位置處3個(gè)方向真實(shí)應(yīng)變之間的數(shù)學(xué)映射關(guān)系,并在螺旋光纖的一個(gè)纏繞周期內(nèi),定量分析了不同纏繞角度、不同旋轉(zhuǎn)角度條件下三分量振動(dòng)信號(hào)對(duì)光纖軸向應(yīng)變的影響,并通過(guò)模擬的三分量數(shù)字信號(hào)對(duì)螺旋光纖軸向方向所接收到的有效信號(hào)進(jìn)行了驗(yàn)證,得出如下結(jié)論:

    1)螺旋纏繞光纖通過(guò)軸向方向的改變能夠采集到光纖所在位置處的三分量振動(dòng)特征,在光纖軸向方向上采集到的信號(hào)強(qiáng)度與實(shí)際三分量信號(hào)強(qiáng)度、光纖纏繞角度和光纖纏繞過(guò)程中各點(diǎn)所在的相對(duì)空間位置有關(guān)。

    2)當(dāng)光纖纏繞軸方向與Z分量信號(hào)方向一致時(shí),隨著光纖纏繞角度α增大,DAS螺旋光纖對(duì)沿纏繞軸方向的地震波振動(dòng)(即Z分量)信號(hào)的敏感性增強(qiáng),而對(duì)垂直于纏繞軸方向的(即X和Y分量)信號(hào)的敏感性減弱。當(dāng)螺旋光纖的纏繞角度固定不變時(shí),螺旋光纖獲得的信號(hào)中Z分量信號(hào)不受旋轉(zhuǎn)角度θ的影響,而X和Y分量信號(hào)則分別以余弦和正弦關(guān)系隨旋轉(zhuǎn)角度θ呈現(xiàn)出周期性的變化規(guī)律,兩個(gè)分量在數(shù)值上具有互補(bǔ)性。

    3)螺旋光纖采集到的信號(hào)取決于光纖上任意點(diǎn)的軸向方向與X,Y,Z3個(gè)分量之間的投影關(guān)系,上述結(jié)論在討論時(shí)假設(shè)光纖纏繞軸方向與Z方向一致,可以直接適用于在直井中開(kāi)展井中地震觀測(cè)時(shí)的情況,而在地面地震勘探時(shí)通??梢詫⒐饫w纏繞軸方向調(diào)整至與地面平行的方向進(jìn)行觀測(cè)。在實(shí)際資料觀測(cè)時(shí),光纖纏繞軸方向與X,Y,Z3個(gè)分量可能均存在一定的夾角,但上述結(jié)論依然成立,不過(guò)規(guī)律性變化特征會(huì)減弱。

    本文重點(diǎn)在均勻介質(zhì)假設(shè)條件下討論了X,Y,Z3個(gè)正交分量振動(dòng)信號(hào)對(duì)螺旋纏繞光纖軸向應(yīng)變的影響,對(duì)不同類型地震波在三維空間3個(gè)方向的振動(dòng)信號(hào)具有普適性。實(shí)際地震勘探中廣泛發(fā)育縱波、橫波和面波等多種類型的波,不同類型的波在不同的傳播方向上的振動(dòng)分量差異較大甚至相位相反,容易導(dǎo)致螺旋光纖軸向應(yīng)變上接收到的信號(hào)發(fā)生更加劇烈的變化,今后需要有針對(duì)性的單獨(dú)討論不同介質(zhì)條件下不同類型地震波的傳播特征和入射角等差異對(duì)光纖應(yīng)變特征的影響,還應(yīng)該充分考慮螺旋纏繞光纖的物性、標(biāo)距長(zhǎng)度以及光纖與纏繞軸之間的耦合程度等因素對(duì)螺旋光纖軸向應(yīng)變響應(yīng)的影響,從而有效提高螺旋纏繞光纖對(duì)不同類型地震波振動(dòng)信號(hào)的采集能力,為真實(shí)三分量振動(dòng)信號(hào)的有效恢復(fù)提供更加合理的依據(jù)。

    猜你喜歡
    軸向坐標(biāo)系分量
    大型立式單級(jí)引黃離心泵軸向力平衡的研究
    帽子的分量
    一物千斤
    智族GQ(2019年9期)2019-10-28 08:16:21
    荒銑加工軸向切深識(shí)別方法
    論《哈姆雷特》中良心的分量
    解密坐標(biāo)系中的平移變換
    坐標(biāo)系背后的故事
    分量
    基于重心坐標(biāo)系的平面幾何證明的探討
    微小型薄底零件的軸向車銑實(shí)驗(yàn)研究
    日本黄色片子视频| 精品人妻熟女av久视频| 亚洲18禁久久av| 国产精品1区2区在线观看.| 久久精品久久久久久久性| 婷婷六月久久综合丁香| 亚洲,欧美,日韩| 亚洲欧美精品综合久久99| 欧美性感艳星| 亚洲欧美精品专区久久| 少妇熟女aⅴ在线视频| 日本黄色片子视频| 激情 狠狠 欧美| 日本色播在线视频| 国产精品99久久久久久久久| 成人av在线播放网站| 天堂√8在线中文| 国内少妇人妻偷人精品xxx网站| 国产色婷婷99| 亚洲欧美一区二区三区国产| 精品久久国产蜜桃| 中文字幕av成人在线电影| 国产午夜精品论理片| 亚洲国产成人一精品久久久| 亚洲欧美日韩卡通动漫| 精品国产三级普通话版| 免费观看在线日韩| 午夜视频国产福利| 成年免费大片在线观看| 日日干狠狠操夜夜爽| 欧美精品国产亚洲| 只有这里有精品99| 亚洲色图av天堂| 可以在线观看毛片的网站| 久久精品国产亚洲网站| 女人十人毛片免费观看3o分钟| 亚洲伊人久久精品综合 | 成人毛片a级毛片在线播放| 亚洲人成网站在线观看播放| 熟女电影av网| 精品久久久久久久久久久久久| 久久精品影院6| 一个人看视频在线观看www免费| 成人av在线播放网站| 18禁在线播放成人免费| 精品一区二区免费观看| 夜夜看夜夜爽夜夜摸| 久久人妻av系列| 亚洲精品国产av成人精品| 国产av不卡久久| 男人舔奶头视频| 亚洲最大成人中文| 老女人水多毛片| 99热全是精品| 91在线精品国自产拍蜜月| 精品一区二区免费观看| 国产av在哪里看| 国语自产精品视频在线第100页| 一级爰片在线观看| 亚洲欧美中文字幕日韩二区| 日本猛色少妇xxxxx猛交久久| 亚洲精品乱码久久久久久按摩| 高清视频免费观看一区二区 | 69av精品久久久久久| 亚洲一级一片aⅴ在线观看| 91av网一区二区| 全区人妻精品视频| 两性午夜刺激爽爽歪歪视频在线观看| 男人狂女人下面高潮的视频| 最近手机中文字幕大全| 中文字幕熟女人妻在线| 黑人高潮一二区| 热99在线观看视频| 免费看a级黄色片| 黄色日韩在线| 97超碰精品成人国产| 一级爰片在线观看| 男人舔女人下体高潮全视频| 亚洲成av人片在线播放无| 精品人妻一区二区三区麻豆| 国产乱来视频区| 日韩欧美精品免费久久| 最近视频中文字幕2019在线8| 久久国产乱子免费精品| 国产午夜精品一二区理论片| 日日摸夜夜添夜夜爱| 亚洲中文字幕一区二区三区有码在线看| 成年免费大片在线观看| 精华霜和精华液先用哪个| 国产成人a∨麻豆精品| 欧美zozozo另类| 亚洲精品乱久久久久久| 亚洲乱码一区二区免费版| 99久久精品热视频| 国产av一区在线观看免费| 国产成年人精品一区二区| 纵有疾风起免费观看全集完整版 | 伦理电影大哥的女人| 禁无遮挡网站| 噜噜噜噜噜久久久久久91| 亚洲av男天堂| 久久久午夜欧美精品| 色综合色国产| 青春草视频在线免费观看| 亚洲国产色片| 91在线精品国自产拍蜜月| 亚洲自偷自拍三级| 亚洲成人中文字幕在线播放| 国产av码专区亚洲av| 久久久久免费精品人妻一区二区| 精品99又大又爽又粗少妇毛片| 一本久久精品| av在线天堂中文字幕| 日本免费a在线| 久久99热这里只有精品18| 亚洲精品aⅴ在线观看| 中文字幕av在线有码专区| 久久久国产成人精品二区| 99热全是精品| 亚洲欧美中文字幕日韩二区| 亚洲va在线va天堂va国产| 国产免费视频播放在线视频 | 亚洲av免费在线观看| www.色视频.com| 久久久久精品久久久久真实原创| 91aial.com中文字幕在线观看| 国产精品综合久久久久久久免费| 国产免费福利视频在线观看| 免费无遮挡裸体视频| 少妇熟女aⅴ在线视频| 国产免费男女视频| 少妇的逼水好多| 欧美97在线视频| 99热这里只有是精品50| 午夜久久久久精精品| 日本三级黄在线观看| 中文字幕制服av| 如何舔出高潮| 免费看av在线观看网站| 久久人人爽人人片av| 久久精品夜色国产| 乱系列少妇在线播放| 国产精品熟女久久久久浪| 最近手机中文字幕大全| 一区二区三区乱码不卡18| 日日摸夜夜添夜夜爱| 青春草亚洲视频在线观看| 亚洲熟妇中文字幕五十中出| 国产探花极品一区二区| 亚洲欧美一区二区三区国产| 菩萨蛮人人尽说江南好唐韦庄 | av黄色大香蕉| 日本av手机在线免费观看| 91久久精品电影网| 国产色爽女视频免费观看| 丰满少妇做爰视频| 免费播放大片免费观看视频在线观看 | 欧美人与善性xxx| 97超碰精品成人国产| 18+在线观看网站| 日韩在线高清观看一区二区三区| 日韩成人伦理影院| av在线老鸭窝| videos熟女内射| 亚洲国产精品久久男人天堂| 99久国产av精品| 亚洲在线自拍视频| 国产白丝娇喘喷水9色精品| 国产精品无大码| 搡女人真爽免费视频火全软件| 日本色播在线视频| av国产免费在线观看| 3wmmmm亚洲av在线观看| 91av网一区二区| 最新中文字幕久久久久| 国产精品av视频在线免费观看| 免费播放大片免费观看视频在线观看 | 国产精品乱码一区二三区的特点| 亚洲性久久影院| 熟妇人妻久久中文字幕3abv| 99久久九九国产精品国产免费| 久久久久久久久久久免费av| 国产黄片美女视频| 伦理电影大哥的女人| 在线a可以看的网站| 九九在线视频观看精品| 你懂的网址亚洲精品在线观看 | 欧美一区二区国产精品久久精品| 国产精品日韩av在线免费观看| av又黄又爽大尺度在线免费看 | 久久久精品欧美日韩精品| 日韩强制内射视频| 深爱激情五月婷婷| 男女国产视频网站| 亚洲欧洲日产国产| 成年女人永久免费观看视频| 欧美一区二区国产精品久久精品| 欧美一区二区国产精品久久精品| 精品久久久久久久久av| 综合色av麻豆| 久久久久久大精品| 国产精品精品国产色婷婷| 国产亚洲av嫩草精品影院| 日韩高清综合在线| 日日干狠狠操夜夜爽| 搡老妇女老女人老熟妇| 一二三四中文在线观看免费高清| 免费大片18禁| 日日干狠狠操夜夜爽| 一级毛片我不卡| 99久国产av精品国产电影| 久久久久久久国产电影| 欧美激情久久久久久爽电影| 久久久久精品久久久久真实原创| 少妇猛男粗大的猛烈进出视频 | 亚洲人成网站高清观看| 成人美女网站在线观看视频| 欧美一级a爱片免费观看看| 小蜜桃在线观看免费完整版高清| 波多野结衣巨乳人妻| 免费一级毛片在线播放高清视频| a级一级毛片免费在线观看| 国产精品99久久久久久久久| 看免费成人av毛片| 欧美最新免费一区二区三区| 夫妻性生交免费视频一级片| 国产免费福利视频在线观看| av国产免费在线观看| 婷婷六月久久综合丁香| 国产一区亚洲一区在线观看| 日韩欧美三级三区| 亚洲av免费在线观看| АⅤ资源中文在线天堂| 嫩草影院精品99| 久久久精品欧美日韩精品| 嫩草影院新地址| 性色avwww在线观看| 在线免费十八禁| 欧美3d第一页| 亚洲精品456在线播放app| 中文字幕亚洲精品专区| 天堂影院成人在线观看| 国产精品电影一区二区三区| 尾随美女入室| 色5月婷婷丁香| 国产伦理片在线播放av一区| 免费无遮挡裸体视频| 成年女人看的毛片在线观看| 国内少妇人妻偷人精品xxx网站| 久久欧美精品欧美久久欧美| 欧美bdsm另类| 久久亚洲精品不卡| 国产男人的电影天堂91| 精品午夜福利在线看| 日韩av在线免费看完整版不卡| 亚洲中文字幕日韩| 少妇人妻一区二区三区视频| 久久婷婷人人爽人人干人人爱| 国产精品日韩av在线免费观看| 少妇猛男粗大的猛烈进出视频 | 小说图片视频综合网站| 欧美又色又爽又黄视频| 国内揄拍国产精品人妻在线| av又黄又爽大尺度在线免费看 | 一区二区三区乱码不卡18| 亚洲欧美中文字幕日韩二区| 特大巨黑吊av在线直播| 久久精品熟女亚洲av麻豆精品 | 在现免费观看毛片| 亚洲一区高清亚洲精品| 级片在线观看| 亚洲精华国产精华液的使用体验| 国产精品一区二区在线观看99 | 亚洲美女搞黄在线观看| 亚洲丝袜综合中文字幕| av天堂中文字幕网| 日韩人妻高清精品专区| 26uuu在线亚洲综合色| 久久精品综合一区二区三区| 亚洲欧美日韩卡通动漫| 99热全是精品| 国产大屁股一区二区在线视频| 国产伦在线观看视频一区| 国产女主播在线喷水免费视频网站 | 国产又色又爽无遮挡免| 亚洲精品乱码久久久久久按摩| 久久99热这里只频精品6学生 | 联通29元200g的流量卡| 少妇高潮的动态图| 亚洲精品自拍成人| 一边亲一边摸免费视频| av卡一久久| 亚洲一区高清亚洲精品| 天堂网av新在线| 伦理电影大哥的女人| 婷婷色麻豆天堂久久 | 午夜精品国产一区二区电影 | 精品不卡国产一区二区三区| 国产 一区精品| 国产视频首页在线观看| 欧美3d第一页| 亚洲国产日韩欧美精品在线观看| 黄色日韩在线| 国产女主播在线喷水免费视频网站 | 国产探花在线观看一区二区| 国产精品一区二区三区四区免费观看| 亚洲精品乱码久久久久久按摩| 我的女老师完整版在线观看| 综合色丁香网| 国产精品美女特级片免费视频播放器| 观看免费一级毛片| 国产伦在线观看视频一区| 免费无遮挡裸体视频| 超碰97精品在线观看| 国产三级在线视频| 看十八女毛片水多多多| 亚洲丝袜综合中文字幕| 麻豆乱淫一区二区| 国产日韩欧美在线精品| 国产av码专区亚洲av| 久久婷婷人人爽人人干人人爱| 久久99热这里只频精品6学生 | 夜夜爽夜夜爽视频| 国产一级毛片七仙女欲春2| 精品久久久久久电影网 | 久久久久久久久中文| 欧美三级亚洲精品| 色综合亚洲欧美另类图片| 69av精品久久久久久| 国产精品久久久久久久久免| 少妇的逼水好多| 22中文网久久字幕| 久久人人爽人人爽人人片va| 国模一区二区三区四区视频| 内地一区二区视频在线| 国产成人freesex在线| 中文欧美无线码| 我要搜黄色片| 三级国产精品片| 天堂av国产一区二区熟女人妻| 小说图片视频综合网站| 国产精品一及| 91在线精品国自产拍蜜月| 99热全是精品| 一级毛片aaaaaa免费看小| 欧美最新免费一区二区三区| 五月伊人婷婷丁香| 舔av片在线| 最近手机中文字幕大全| 最近2019中文字幕mv第一页| 国产探花在线观看一区二区| 亚洲电影在线观看av| 亚州av有码| 精品久久久久久久末码| 中文在线观看免费www的网站| 最近中文字幕2019免费版| 久久久久久伊人网av| 午夜福利在线观看吧| 国产黄色视频一区二区在线观看 | 婷婷色综合大香蕉| 午夜亚洲福利在线播放| 国产精品一二三区在线看| 日韩国内少妇激情av| 国产精品无大码| 狠狠狠狠99中文字幕| 久久精品久久久久久久性| 97超碰精品成人国产| 国产三级在线视频| 综合色丁香网| 成年av动漫网址| 久久久久久久久久久免费av| 久久这里只有精品中国| 一级毛片我不卡| 亚洲av成人精品一区久久| 精品一区二区三区人妻视频| 精品久久国产蜜桃| 狂野欧美白嫩少妇大欣赏| 国产精品不卡视频一区二区| 少妇人妻一区二区三区视频| 久久6这里有精品| 内地一区二区视频在线| 色网站视频免费| 国产高清视频在线观看网站| 男人和女人高潮做爰伦理| 啦啦啦韩国在线观看视频| 国内少妇人妻偷人精品xxx网站| 国产视频首页在线观看| 国产单亲对白刺激| 亚洲国产精品合色在线| 久久婷婷人人爽人人干人人爱| 国产精品熟女久久久久浪| 成人美女网站在线观看视频| 精品人妻偷拍中文字幕| 久久久久久大精品| 国产69精品久久久久777片| 久久久久国产网址| 亚洲欧洲国产日韩| videossex国产| 十八禁国产超污无遮挡网站| 国产高潮美女av| 色综合站精品国产| 干丝袜人妻中文字幕| 少妇高潮的动态图| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲四区av| 国产免费又黄又爽又色| 男的添女的下面高潮视频| 噜噜噜噜噜久久久久久91| 又爽又黄无遮挡网站| 一个人看的www免费观看视频| 免费不卡的大黄色大毛片视频在线观看 | 又粗又爽又猛毛片免费看| 一级毛片久久久久久久久女| 亚洲人成网站在线播| 国产极品精品免费视频能看的| 麻豆国产97在线/欧美| 美女高潮的动态| 少妇高潮的动态图| 国产精品一区二区三区四区久久| 亚洲人与动物交配视频| 伦精品一区二区三区| 国产一区二区亚洲精品在线观看| 韩国av在线不卡| 亚洲色图av天堂| 国产高清视频在线观看网站| 免费黄色在线免费观看| 九草在线视频观看| 国产欧美另类精品又又久久亚洲欧美| 成人性生交大片免费视频hd| 免费av观看视频| 久久精品国产自在天天线| 亚洲最大成人手机在线| 99久久精品一区二区三区| 国产亚洲5aaaaa淫片| 欧美zozozo另类| 中文字幕人妻熟人妻熟丝袜美| 国内少妇人妻偷人精品xxx网站| 一级毛片久久久久久久久女| 久久精品国产鲁丝片午夜精品| 亚洲精品aⅴ在线观看| 我要搜黄色片| 亚洲国产精品久久男人天堂| 国产午夜福利久久久久久| 久久久精品欧美日韩精品| 欧美另类亚洲清纯唯美| 1000部很黄的大片| 校园人妻丝袜中文字幕| 天堂影院成人在线观看| 亚洲精品亚洲一区二区| 亚洲欧美日韩卡通动漫| 黄色欧美视频在线观看| 国产精品无大码| 毛片女人毛片| 日韩成人伦理影院| 国产成人免费观看mmmm| 婷婷色av中文字幕| 国产免费又黄又爽又色| 久久99精品国语久久久| 国模一区二区三区四区视频| 小蜜桃在线观看免费完整版高清| 午夜a级毛片| 亚洲成人久久爱视频| 老司机影院毛片| 噜噜噜噜噜久久久久久91| 色尼玛亚洲综合影院| 男女啪啪激烈高潮av片| 在线观看美女被高潮喷水网站| 久久人妻av系列| 国内少妇人妻偷人精品xxx网站| 97超视频在线观看视频| 国产麻豆成人av免费视频| 国产精品野战在线观看| 天堂中文最新版在线下载 | 99在线视频只有这里精品首页| 波多野结衣巨乳人妻| 国产精品无大码| 午夜福利视频1000在线观看| 国产精品不卡视频一区二区| 精品国产三级普通话版| 国产美女午夜福利| 一级二级三级毛片免费看| 熟女人妻精品中文字幕| 亚洲,欧美,日韩| 色吧在线观看| 国产一级毛片在线| 午夜福利网站1000一区二区三区| 一级二级三级毛片免费看| 尤物成人国产欧美一区二区三区| 日韩欧美 国产精品| 日韩成人av中文字幕在线观看| 久久久久久久久久久丰满| 麻豆久久精品国产亚洲av| 嫩草影院新地址| 国产日韩欧美在线精品| 日本五十路高清| 毛片一级片免费看久久久久| 啦啦啦韩国在线观看视频| 深夜a级毛片| 国产在视频线精品| 联通29元200g的流量卡| 少妇丰满av| 狂野欧美白嫩少妇大欣赏| 99视频精品全部免费 在线| 午夜亚洲福利在线播放| 欧美潮喷喷水| 国产精品一及| 欧美极品一区二区三区四区| 日韩人妻高清精品专区| 亚洲综合精品二区| 不卡视频在线观看欧美| 天堂√8在线中文| 国产在视频线在精品| 亚洲欧美精品专区久久| 99在线视频只有这里精品首页| .国产精品久久| 岛国在线免费视频观看| 九草在线视频观看| 狂野欧美激情性xxxx在线观看| 日韩一区二区视频免费看| 日韩欧美三级三区| 一卡2卡三卡四卡精品乱码亚洲| 神马国产精品三级电影在线观看| 99久久精品国产国产毛片| 国产精品一区二区三区四区免费观看| 久久久久久大精品| 一级爰片在线观看| 久久久成人免费电影| 99久国产av精品| 日产精品乱码卡一卡2卡三| 欧美色视频一区免费| av天堂中文字幕网| 91久久精品国产一区二区成人| 亚洲精华国产精华液的使用体验| 久久人人爽人人爽人人片va| 又爽又黄a免费视频| 男女下面进入的视频免费午夜| 中文字幕精品亚洲无线码一区| 欧美日韩一区二区视频在线观看视频在线 | 美女黄网站色视频| 亚洲av电影不卡..在线观看| 日本与韩国留学比较| 精品无人区乱码1区二区| 永久免费av网站大全| 久久国内精品自在自线图片| 国产精品一区二区在线观看99 | 久久国产乱子免费精品| 久久精品久久久久久久性| 久久久亚洲精品成人影院| 好男人在线观看高清免费视频| 国产熟女欧美一区二区| 亚洲五月天丁香| 99久久九九国产精品国产免费| 校园人妻丝袜中文字幕| 中文资源天堂在线| 乱人视频在线观看| av专区在线播放| 国产国拍精品亚洲av在线观看| 国产精品麻豆人妻色哟哟久久 | www.色视频.com| 99久久成人亚洲精品观看| 少妇被粗大猛烈的视频| 又爽又黄无遮挡网站| 中文欧美无线码| 少妇熟女欧美另类| 亚洲精品亚洲一区二区| 老司机影院成人| 成年女人看的毛片在线观看| 男插女下体视频免费在线播放| 在线a可以看的网站| 97热精品久久久久久| 三级经典国产精品| 亚洲av熟女| 欧美日韩国产亚洲二区| 淫秽高清视频在线观看| 97超碰精品成人国产| 国产精品爽爽va在线观看网站| 精品久久久久久久人妻蜜臀av| 国产精品三级大全| 免费一级毛片在线播放高清视频| 一边亲一边摸免费视频| 久久精品熟女亚洲av麻豆精品 | 日本-黄色视频高清免费观看| 美女cb高潮喷水在线观看| 亚洲中文字幕日韩| 丰满少妇做爰视频| 亚洲国产高清在线一区二区三| 久久久久九九精品影院| 成年女人看的毛片在线观看| 男人的好看免费观看在线视频| 在线a可以看的网站| 一区二区三区乱码不卡18| 99视频精品全部免费 在线| 性色avwww在线观看| av.在线天堂| 久久久久网色| 又粗又爽又猛毛片免费看| 久久精品夜色国产| 欧美又色又爽又黄视频| 久久99热6这里只有精品| 亚洲欧美精品专区久久| 成人鲁丝片一二三区免费| 91av网一区二区| 成人亚洲精品av一区二区| 日本wwww免费看| 高清在线视频一区二区三区 | 国产高清有码在线观看视频| 久久鲁丝午夜福利片| 嫩草影院入口| 国产av码专区亚洲av| 国产精品熟女久久久久浪| 国内少妇人妻偷人精品xxx网站| 国产黄色视频一区二区在线观看 | 国产片特级美女逼逼视频| 日本免费一区二区三区高清不卡| 国产成人aa在线观看|