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

    基于最大熵準則的Wigner-Ville分布與微型古河道刻畫方法及應用

    2022-01-08 08:48:22徐天吉程冰潔牛雙晨秦正曄王貞貞
    石油勘探與開發(fā) 2021年6期
    關鍵詞:時頻振幅頻譜

    徐天吉 ,程冰潔 ,牛雙晨,秦正曄,王貞貞

    (1. 電子科技大學資源與環(huán)境學院,成都 611731;2. 電子科技大學長三角研究院(湖州),浙江湖州 313000;3. 成都理工大學“油氣藏地質及開發(fā)工程”國家重點實驗室,成都 610059;4. 成都理工大學“地球探測與信息技術”教育部重點實驗室,成都 610059;5. 四川中成煤田物探工程院有限公司,成都 610072)

    0 引言

    平直、蛇曲、辮狀等類型多樣的河流相,多年前就已經成為油氣領域的研究熱點。在俄羅斯西西伯利亞、美國懷俄明州及中國的鄂爾多斯盆地、四川盆地等區(qū)域,已經發(fā)現(xiàn)了許多河流相大型和中型油氣田。但是隨著中國油氣勘探開發(fā)程度的不斷深化,目前河流相資源目標已經由大型、中型逐漸轉向微型油氣藏。微型古河道廣泛分布在陸相與海陸過渡相的淺、中、深等各類地層,油氣儲量豐富,勘探潛力巨大,但具有河道寬度窄、砂體薄、油氣聚集復雜等不利特征[1-2]。此外,隨著埋藏深度的增加,地震反射表現(xiàn)出能量弱、低頻為主、頻寬窄等不利特征,垂向與橫向分辨率降低,制約了窄古河道、薄儲集層及其他油氣信息的精確提取。微型古河道的識別,已經成為增加井位部署、水平井鉆采等勘探開發(fā)風險的突出難題。

    針對微型古河道的識別難題,目前主要從地震采集、處理和解釋等角度著手,利用高分辨地震數據獲取古河道的橫向寬度、垂向厚度等關鍵參數和空間展布特征。就高分辨地震采集而言,盡管在高密度、“兩寬一高”(寬方位、寬頻帶、高密度)等先進觀測系統(tǒng)下采集到的地震資料,提升了淺中層油氣藏預測與描述精度,但是面對深層地質目標,至今未能突破震源能量衰減快、檢波器難以接收到高頻有效信號等激發(fā)與接收方面的“瓶頸”問題。就高分辨地震資料處理與解釋而言,盡管形成了反Q補償[3]、譜整形[4]、頻寬拓展[5]、薄層反射系數反演[6]、低頻擴展深度學習[7]、高分辨阻抗反演[8-9]等方法,可以使地震數據“無中生有”地出現(xiàn)更多的反射軸、波阻抗界面增加等現(xiàn)象,但存在信噪比降低、波組有效性評估難、甚至出現(xiàn)假反射界面等特征,必然降低地質預測的可靠性。同時,埋藏越深、規(guī)模越小的微型古河道的識別難度越大;尤其是極窄、流向復雜的微型古河道,需要從地震波的傳播能量、反射頻率等角度出發(fā),不斷探索突破地震分辨率極限的有效方法,才能逐步解決微型古河道的識別難題。

    識別微型古河道的關鍵,是基于高分辨時頻分析方法提取地震信號的高分辨時頻響應特征,利用地震波的能量、頻率、相位等信息,有效刻畫微型古河道的寬度、流向等特征。目前,反Q補償、頻寬拓展等許多高分辨處理與解釋方法,均以時頻分析為基礎,卻未實現(xiàn)高分辨時頻響應特征提取,致使類似微型古河道的小規(guī)模地質目標難以實現(xiàn)高清識別。其實,自1946年Gabor首次提出了聲波信號的頻譜圖后[10],迅速形成了線性、雙線性和非線性 3大類核函數不同的時頻分析方法[11]。其中,Gabor變換、短時傅里葉變換(STFT)、連續(xù)小波變換(CWT)、S變換(ST)和廣義 S變換(GST)等屬于線性時頻分析方法,雖具有計算效率較高的優(yōu)勢,但多數具有窗口依賴性,且受海森堡(Heisenberg W)不確定性原理控制[12],時間和頻率的分辨能力不能同時提升。自適應 STFT[13]、模態(tài)分解[14-15]、同步擠壓小波變換[16]、匹配追蹤(MP)[17]等屬于非線性時頻分析方法,雖具有分辨率較高的優(yōu)勢,但計算量較大,魯棒性不足,難以適應低信噪比信號時頻分析。Wigner-Ville分布(WVD)[18]、偽Wigner-Ville分布(PWVD)、平滑偽Wigner-Ville分布(SPWVD)和Choi-Williams分布(CWD)[19]等屬于雙線性(又稱二次型)時頻分析方法,不受海森堡不確定性原理約束[12],優(yōu)良的數學性質為高分辨時頻分析奠定了良好的理論基礎。

    本文基于信號最大熵準則(Maximum Entropy Principle),在克服 WVD交叉項干擾的前提下,研究增強 WVD頻譜聚焦特性和提升地震信號時頻分辨率的方法,以突出窄細古河道響應特征。并采用仿真信號、正演數值模擬等理論數據開展試驗與論證分析,以形成微型古河道識別新方法,為井位部署、水平鉆井等提供支撐。

    1 方法原理

    1.1 最大熵準則

    熵(Entropy),是指離散隨機事件或特定信息的出現(xiàn)概率,由Shannon C E在1948年首次定義,主要用于表征隨機變量或特征信息的不確定性[20]。1957年,Jaynes E T 提出了最大熵準則,即在先驗信息準確的前提下,最能代表當前信息狀態(tài)的概率分布,是熵最大的概率分布,換言之,當符合先驗信息的概率分布具有不確定性時,熵最大的信息分布是合理的[21]。

    1967年,Burg將最大熵準則引入信號頻譜分析中,形成了最大熵時頻分析Burg算法[22]。Burg認為,熵是信號X的分布函數,僅與X的取值概率有關,而與實際取值無關。當信號X的取值為x時,分布密度函數為p(x),則其熵H(x)為:

    此時,Burg定義信號X的功率譜熵為:

    當利用(2)式求取功率譜時,需要在H[P(f)]最大的條件下,遞推求解信號X的自相關函數,使功率譜誤差最小,分辨率獲得提升。其實,Burg算法等價于線性自回歸AR(Auto Regressive)預測模型,最大熵功率譜P(f)的計算公式[22]為:

    1.2 離散Wigner-Ville分布與功率譜

    1992年Boashash定義[23]連續(xù)時序信號x(t)的解析信號為z(t)時,Wigner-Ville分布(WVD)為:

    然而,實測的信號并非連續(xù)信號,實際應用時,需要對W(t,f)進行離散化。設x(t)、z(t)的離散表達式分別為x(n)和z(n)(0 ≤n≤N-1),對x(n)作希爾伯特(Hilbert)變換,在計算出z(n)后即可獲得W(t,f)的離散表述:

    據WVD的性質可知,離散WVD是對瞬時自相關函數進行離散傅立葉變換。離散信號的瞬時自相關函數為:

    對于每一個樣點n,kn(l)將組成一個卷積核序列kz(n),如(7)式。當n=0和n=N- 1 時,卷積核序列kz(n)=kz(0)和kz(n) =kz(N- 1)最短;當時,最長。以kn(0)為中心,對kz(n)左右補充0值,可將kz(n)擴充為長度為N的卷積核序列。

    對第n個采樣點的卷積核序列kz(n),進行離散傅立葉變換,即可獲得WVD瞬時功率譜,如下:

    其中,Wz(n,f)的子項可表述為:

    通過計算每一個采樣點n對應的瞬時功率譜,最終可以得到離散信號x(n)的Wigner-Ville功率譜,即:

    當然,受 WVD雙線性的影響,在利用瞬時自相關函數kn(l)前后數據相乘時,將產生交叉項干擾,不利于高精度功率譜特征分析。

    1.3 最大熵準則約束下的Wigner-Ville功率譜計算方法

    當瞬時自相關函數kn(l)的一階導數為0時,求得的Wigner-Ville分布即為最大熵譜。由于一維最大熵譜與自回歸模型AR譜是等效的,故只需要求得AR模型參數、計算出 AR模型譜,就能獲得最大熵譜[24]。這樣,可以避免利用(10)式計算功率譜,繞開瞬時自相關函數前后數據相乘時產生的交叉項干擾。

    基于 Burg算法求取最大熵功率譜P(f)時,(3)式與(10)式是等效的,只需基于Levinson-Durbin遞推規(guī)則[25],在求取預測誤差E、自回歸系數a等 AR模型的參數后,代入(3)式就能獲得最大熵準則約束下的WVD功率譜。Burg算法通過前向預測誤差和后向預測誤差,遞推 AR模型的自回歸系數[26-27]。將設為初始條件,遞推公式如下:

    通過離散信號x(n)及其解析信號z(n),可以計算預測誤差的初始功率,即

    利用m階前向預測誤差和后向預測誤差,可以計算預測誤差的平均功率,即

    據Levinson-Durbin遞推約束規(guī)則,設am(j)為M階AR模型在階次序號為m時的第j個自回歸系數,Em為m階次時的最小預測誤差,則其遞推關系為

    這樣,將am( 1),am( 2),…,am(m- 1)和Em代入(3)式,就能計算出 WVD最大熵功率譜。WVD最大熵功率譜在時間域和頻率域的分布,即為Wigner-Ville最大熵時頻譜(簡稱MEWVD)。

    概括起來,Burg算法的WVD最大熵功率譜的計算步驟[28]主要包括:①設m=1,計算P0(f)、EF,0(n)與EB,0(n);②計算am(m);③計算階次為m時的第j個自回歸系數④計算Pm(f);⑤計算階次為m時的EF,m(n)和EB,m(n);⑥設m=m+1,重復②—⑤,直到Pm(f)不再明顯減小時,利用Levinson-Durbin遞推關系,可以獲得AR模型的參數,最終計算出WVD最大熵功率譜。

    在時間域,MEWVD并未將瞬時自相關函數前后數據相乘,而僅使用當前n點附近的數據,避免了交叉項干擾;在頻率域,受最大熵準則約束,使每個n點對應的功率譜更精確,能量更聚焦。同時,在最大熵的前提下,采用Levinson-Durbin約束遞推未知瞬時自相關函數,使計算的功率譜誤差減小,提升了信號時頻分辨能力。

    需要指出的是,最大熵準則在數學、經濟、信息等領域被廣泛應用,主要優(yōu)勢在于其不僅保留了信息分布的不確性時,同時將預測風險降到了最低。在油氣領域,常常面對橫向較窄、縱向較薄的古河道、砂壩、點礁等微型地質體,地震波衰減、調諧、噪聲等因素增加了微型地質體的預測難度,計算高精度、高聚焦性的地震功率譜,是識別微型地質體的重要基礎。MEWVD在熵最大時計算最能代表微型地質體的地震頻譜響應特征,不僅能實現(xiàn)時頻特征的分類取優(yōu),而且可以在噪聲環(huán)境下提高信號分析的魯棒性[29],有利于提升微型地質體識別的可靠性。

    2 方法實驗

    2.1 仿真地震信號的MEWVD實驗及定量分析

    地震波在地下傳播時,受地層埋深、孔隙度、填充流體等因素影響,將產生能量衰減,引起動力學、運動學和幾何學等特征發(fā)生變化。采用振幅、吸收系數、頻率、相位和走時等參數,可以仿真描述地震波傳播過程,如下:

    上式中,仿真地震信號xF(t)分別由衰減分量x1(t)和平穩(wěn)分量x2(t)疊加而成。其中,x1(t)包含了吸收系數、振幅、頻率、相位和走時等參數,x2(t)不包含吸收系數,傳播能量將不會產生衰減。

    為了對比分析地震波傳播過程中的時頻特征,利用xF(t)的參數,構造出了關聯(lián)特征較強的高、低頻能量強弱不一樣的仿真地震信號。如表1,信號xF1和xF2的起始能量A1和傳播時間t相同,xF1中A2、f2、j2為0,表示xF1僅由衰減分量構成,且其A1、f1和Q均較大,類似強能量、高頻和強衰減的地震信號。xF2包含了xF1,僅通過改變f2,在xF2的基礎上增加了低頻、弱振幅的平穩(wěn)分量,意味著xF2在t時間內的傳播能量、振幅等將隨頻率變化而改變。如圖 1a和圖 1b中第 1列,分別顯示了xF1和xF2的波形特征??梢姡瑇F1隨著t的變化,振幅減弱,能量衰減很快;xF2隨著t的變化,振幅減弱,能量也在衰減,但衰減分量的控制作用越來越弱,在信號傳播約0.25 s之后,信號主要受平穩(wěn)分量控制,且由于f2差異,使xF1振蕩周期短、波形“瘦”,xF2振蕩周期長、波形“胖”。

    表1 仿真地震信號的參數描述

    同時,對比不同的時頻分析方法,發(fā)現(xiàn) MEWVD具有非常優(yōu)良的時頻聚焦特性和分辨率。如圖1a和圖1b,第 2列顯示了仿真地震信號的快速傅里葉變換(FFT)頻譜,揭示信號xF1的能量峰值聚集在80 Hz,xF2的能量峰值聚集在10 Hz和80 Hz。第3至10列,分別顯示了 STFT、GST、CWT、WVD、PWVD、SPWVD、CWD和 MEWVD等方法的時頻效果。對比分析,盡管xF1和xF2的最強能量均聚集在 80 Hz高頻位置,但是,不同的時頻分析方法顯示的能量聚焦程度與頻率分辨能力存在差異,且MEWVD聚焦效果與分辨率最優(yōu)、WVD次之、STFT最弱。GST和CWT的時頻分辨能力相似,對80 Hz高頻強振幅與10 Hz低頻弱振幅均能有效識別,但能量聚焦程度、時間分辨能力等較MEWVD差。此外,WVD時頻特征還存在交叉項干擾,PWVD雖在一定程度上壓制了交叉項干擾,但能量聚焦性減弱,降低了xF2的低頻、弱振幅信息的分辨能力;SPWVD和CWD基本壓制了交叉項干擾,但也幾乎無法有效識別xF2中低頻、弱振幅響應,較 PWVD的分辨率損失更嚴重。

    圖1 仿真地震信號及不同時頻分析方法呈現(xiàn)出的頻譜響應差異

    2.2 窄薄數值模型MEWVD實驗及效果分析

    為了進一步證實 Wigner-Ville最大熵時頻譜(MEWVD)的頻譜聚焦性和時頻高分辨能力,采用表2所示參數,設計出二維數值模型,包含了埋深、寬度和厚度各異的窄薄地質體,其“透鏡狀”的地質特征與微型古河道的橫截面類似。基于二維地震波動方程正演數值模擬方法,以35 Hz的Rick子波作為點震源激發(fā),以10 m采樣間距、1 ms采樣率接收,獲得了正演模擬記錄。模型④—⑨的砂體厚度逐漸變薄、寬度逐漸變窄;模型④—⑥較寬、厚度大于 1/4波長(約30 m),頂底界面的反射同相軸明顯,而⑦—⑨較窄、厚度小于1/4波長,地震記錄表現(xiàn)為諧波反射,且⑦和⑧表現(xiàn)出強調諧反射,極薄、極窄的⑨則表現(xiàn)為弱調諧反射(見圖2)。

    表2 數值模型參數描述

    圖2 窄薄各異模型的單道模擬地震信號及MEWVD等頻譜特征

    分別抽取代表模型④—⑨的單道地震信號,計算MEWVD后,可分析模型的時頻差異。如圖2,顯示了模型④—⑨的單道地震信號及MEWVD頻譜特征,可見,窄薄各異的數值模型的反射信號,MEWVD均表現(xiàn)出了較強的頻譜聚焦性和高分辨特征。其中,模型④—⑥較寬、較厚,單道地震信號表現(xiàn)出明確的頂底界反射,MEWVD在頻率軸方向表現(xiàn)出較窄的頻譜分布。模型⑦—⑨較窄、較薄,單道地震信號僅顯示出了頂界面的反射,底界面受波場調諧作用而難以識別。MEWVD在頻率軸方向表現(xiàn)出較寬的頻譜分布,且厚度越薄、頻譜分布范圍越寬。

    在單頻剖面上,MEWVD的頻譜聚焦性非常顯著,不同頻率的功率譜展現(xiàn)出了不同的時頻分辨率。如圖3,顯示了正演地震信號的瞬時振幅及25,35,45 Hz的MEWVD頻譜特征。其中,圖3a顯示瞬時振幅受時頻分辨能力的局限,僅較厚的模型④表現(xiàn)出了頂底界面清晰的“雙軸”強瞬時振幅特征,而厚度逐漸減薄的模型⑤—⑨卻表現(xiàn)為“單軸”強瞬時振幅分布,連厚度大于1/4波長的模型⑤和⑥的頂、底界面也難以有效識別。圖3b—圖3c分別顯示了3種頻率的MEWVD頻譜分布,且均比瞬時振幅具有更高的分辨率,但不同頻率的 MEWVD頻譜響應各異。其中,45 Hz的MEWVD頻譜分辨率最高,可以準確識別窄薄模型⑤和⑥的頂底界面及橫向寬度。35 Hz的MEWVD頻譜分辨率次之,也能夠有效識別⑤和⑥的頂底界面及橫向寬度,對極窄、極薄模型⑧和⑨的頻譜聚焦響應優(yōu)良,且受不同寬度的影響,二者的頻譜分布厚度相近、寬度差異顯著。進一步分析極窄、極薄模型⑧和⑨,二者厚度相同、寬度差 5倍,瞬時振幅橫向差異不明顯,MEWVD頻譜異常顯著。尤其在圖3d中,45 Hz高頻聚焦性極好,顯示出了厚度一樣、寬度差 5倍的MEWVD頻譜異常特征。同時,雖然25 Hz的MEWVD頻譜分辨率總體不如35 Hz和45 Hz的頻譜,但是對極窄、極薄模型⑧和⑨的頻譜聚焦響應卻有明顯差異。其中,35 Hz和45 Hz的頻譜表現(xiàn)為“團狀”現(xiàn)象,而25 Hz的頻譜卻表現(xiàn)出“層狀”頻譜聚焦現(xiàn)象,表明Wigner-Ville最大熵譜與地震信號的強弱無關,當調諧作用使模型頂底界面的功率譜分布不確定時,熵最大的功率譜即為最合理的分布。

    圖3 正演地震信號的瞬時振幅及25,35,45 Hz的MEWVD頻譜特征

    總之,仿真地震信號和窄薄正演數值模型實驗表明,MEWVD具有良好的頻譜聚焦性和時頻高分辨能力。MEWVD能夠精確表征信號的時頻分布,且不同頻率的MEWVD能夠準確獲取不同尺度窄薄模型的頻譜響應;同時,由于Wigner-Ville最大熵譜僅與頻譜分布概率相關,當地質體較窄或薄度小于1/4波長而發(fā)生調諧作用時,利用地震反射同相軸不能直接識別,熵最大的功率譜即為地質體的最合理的分布空間。因此,綜合不同頻率的MEWVD頻譜特征,可以識別類似窄薄模型的微型古河道的空間分布。

    3 應用案例

    四川盆地是一個多旋回富氣沉積盆地,截至2020年,在震旦系、寒武系、志留系、石炭系、三疊系、侏羅系等27套層系中,已經發(fā)現(xiàn)了構造、構造-巖性、巖性-構造和巖性4類氣藏,形成了四川盆地東北部、西部、南部和北部 4大工業(yè)氣區(qū)[30-31]。其中,川西氣區(qū)的侏羅系沙溪廟組氣藏,主要分布在西鄰龍門山推覆構造帶、東接川中隆起區(qū)的川西坳陷。在川西東斜坡地區(qū),發(fā)現(xiàn)了埋深小于3 100 m的大型巖性-構造氣田——中江氣田。

    中江氣田面積為2 350 km2,區(qū)內侏羅系沙溪廟組以淺水三角洲沉積為主,在三角洲平原亞相中發(fā)育分流河道、分流間灣、河口壩、決口扇等沉積微相。其中,分流河道主要發(fā)育細—中粒巖屑長石砂巖、長石巖屑砂巖、巖屑砂巖和巖屑石英砂巖等砂巖儲集層,平均孔隙度達 8.94%,平均滲透率達0.5×10-3μm2。目前,針對資源豐富的沙溪廟組氣藏,已投產工業(yè)氣井約150口,日產氣超過 220′104m3。然而,隨著開發(fā)程度的深化,微型分流河道氣藏更加隱蔽,儲集層以條帶狀疊置薄層河道砂體為主,橫向寬度窄,縱向厚度薄,精細刻畫難度極大。需要探索窄細微型古河道的精確識別方法,為氣藏描述、水平鉆井、水力壓裂等提供支撐。

    中江氣田侏羅系中統(tǒng)沙溪廟組劃分為上(J2s1)、中(J2s2)、下(J2s3)共 3段;其中,J2s3段由 J2s31、J2s32和 J2s33組成。依據古河道和砂巖沉積特征,J2s33層被細分為J2s33-1和J2s33-2小層。目前,J2s33-2小層是主力產氣層,層內網狀和條帶狀古河道主要沿北東流向南西,延伸長度18~58 km,寬度0.2~3.5 km,砂厚6.5~45.2 m。儲集層呈席狀、箱狀和鐘型分布,表現(xiàn)為低頻、“亮點”、“空白”等地震反射特征,沿古河道方向,地震反射總體較連續(xù)、穩(wěn)定,垂直河道方向,呈“透鏡狀”或弱反射異常。

    近幾年,隨著中江氣田 J2s33-2小層氣藏的深入開發(fā),寬度小于500 m、厚度小于35 m的窄薄古河道砂巖成為重點目標,微型古河道的精確識別也成為氣藏開發(fā)的關鍵環(huán)節(jié)。然而,受埋藏深度、寬度、厚度、巖性、物性等多種因素的影響,不同的微型古河道呈現(xiàn)出明顯的地震反射差異。如圖4a圖顯示J2s33-2小層的地震主頻約30 Hz,頻帶8~70 Hz;圖4b顯示E井J2s33-2小層古河道出現(xiàn)砂體GR低值異常,地震反射為局部雙波峰“透鏡狀”異常,預測河道寬約354 m、厚約35 m;圖4c顯示D井J2s33-2小層古河道出現(xiàn)GR低值異常,預測河道砂巖寬約310 m、厚約30 m,較E井河道更窄、更薄,“透鏡狀”反射更弱;圖4d顯示C井J2s33-2小層古河道出現(xiàn)GR低值異常,預測河道砂巖寬約248 m、厚約30 m,與D、E井相比較,“透鏡狀”異常消失,呈弱反射響應。顯然,通過對比圖4中的C、D、E井微型古河道地震反射特征,隨著寬度變窄、厚度減薄,河道隱藏性增強,地震反射越弱,識別難度越大。

    圖4 目標地層地震振幅譜及典型井微型分流河道地震反射剖面

    采用Wigner-Ville最大熵時頻譜計算方法,提取了中江氣田沙溪廟組的MEWVD頻譜特征,發(fā)現(xiàn)不同寬度和厚度的微型古河道呈現(xiàn)出了不同的MEWVD異常,揭示了窄細古河道的空間展布。如圖5a顯示的瞬時振幅剖面上,E井J2s33-2小層古河道出現(xiàn)瞬時振幅強異常,異常范圍較大,聚焦性較差,分辨率較低;圖5b—圖5d顯示的MEWVD頻譜剖面上,E井J2s33-2小層古河道出現(xiàn)強異常,異常范圍隨著頻率的增加而減小,聚焦能力增強,分辨率提高。比較瞬時振幅與25 Hz、35 Hz的MEWVD頻譜特征,在45 Hz的MEWVD剖面上,分辨能力更高,縱向頻譜被分開,橫向頻譜收窄,頻譜異常由“團狀”變化為“雙層狀”,與砂體GR低值異常更加吻合,有效地揭示出了E井J2s33-2小層35 m薄古河道砂巖的頂底界面??梢?,利用MEWVD頻譜異常,能夠有效刻畫窄細古河道的縱向分布。

    圖5 E井瞬時振幅及25,35,45 Hz的MEWVD頻譜特征

    同時,利用MEWVD頻譜異常,還能刻畫微型古河道的橫向展布。如圖6a瞬時振幅高值異??坍嫵隽讼鄬^寬的Ⅰ—Ⅳ、Ⅵ、Ⅷ號主河流的沿層展布;但在黑色框示區(qū)域,河道很窄、河道砂巖很薄,地震反射能量較強,分辨率低,利用瞬時振幅無法識別窄細古河道的空間位置、寬度和流向。在圖 6b中,利用MEWVD方法獲得的25 Hz頻譜出現(xiàn)高值異常,不僅刻畫出了較寬的Ⅰ—Ⅳ、Ⅵ、Ⅷ號主河流的沿層展布,而且中低值異??坍嫵隽烁蘑鹾廷魈柡拥赖目臻g展布;尤其在左側黑色框內,Ⅴ和Ⅶ號河流的寬度和流向非常清晰。與Ⅰ—Ⅷ微型河道比,Ⅸ號河道更窄,且流向彎曲變化較大,隱蔽性更強,利用瞬時振幅和25 Hz頻譜難以識別,而在35 Hz和45 Hz的MEWVD頻譜特征中,Ⅸ號曲流河道的沿層特征卻非常清晰。由圖6可見,MEWVD較瞬時振幅分辨更高,且不同頻率的MEWVD頻譜特征,可以識別不同寬度的窄細河道。當然,由于任何物質都有相應的敏感頻段,河道也不例外,Ⅴ和Ⅶ號河道對低頻更敏感,故在35 Hz和45 Hz的MEWVD頻譜特征中反而不如25 Hz刻畫更清晰(見圖6黑色方框),但整體規(guī)律仍然是高頻頻譜刻畫效果更佳。

    圖6 瞬時振幅及25,35,45 Hz的MEWVD分別刻畫的J2s33-2小層微型古河道的展布

    采用RGB融合(紅、綠、藍混合生成其他色彩)顯示方式,將不同頻率的MEWVD頻譜進行融合,能更清晰的展示不同寬度和厚度河道的空間分布。如圖7a顯示了采用 RGB融合方式,將 25,35,45 Hz的MEWVD頻譜進行融合,有效刻畫出J2s33-2小層Ⅰ—Ⅸ號古河道的展布;尤其是具有極窄特征的Ⅴ、Ⅶ和Ⅸ號河道,以及Ⅷ號河道以東、Ⅵ號河道以西的更多微型古河道(紅色箭頭所指),也精確地揭示出了其寬度和流向等空間信息。當然,結合區(qū)內沉積構造背景,利用巖心分析與測井解釋資料,可以證實采用MEWVD

    頻譜異常識別微型古河道的可靠性。如圖7b所示,在縱向上,D井主要發(fā)育分流河道、分流間彎、河口壩等三角洲前緣沉積微相。其中,J2s33-2小層以分流河道沉積微相為主,聲波時差和伽馬測井曲線表現(xiàn)出低值響應特征,深側向和淺側向電阻率曲線表現(xiàn)出高值異常。在分流河道中部2 869.57~2 871.49 m井段的鉆井取心顯示,古河道水平層理發(fā)育,在較強的水動力作用下,沉積形成了厚度約30 m的致密砂巖儲集層。受巖性、物性、厚度等因素影響,儲集層與圍巖的阻抗差異不大,在圖4c所示的“透鏡狀”河道反射較弱,但在圖6和圖7卻具有較突出的MEWVD頻譜異常??梢?,綜合地震反射、鉆井取心、測井響應等信息,利用MEWVD頻譜異??坍嬑⑿凸藕拥谰哂休^高可靠性。

    圖7 中江氣田J2s33-2小層RGB融合MEWVD刻畫河道分布及D井沉積相測井解釋

    此外,MEWVD頻譜異常還可以刻畫斷層分布特征,但斷層與古河道具有顯著的MEWVD差異。如圖6d和圖7a顯示,中江氣田沉積構造環(huán)境整體比較穩(wěn)定,在試驗工區(qū)南部僅發(fā)育少量近東西向的斷層,裂縫不發(fā)育。斷層的MEWVD頻譜異常表現(xiàn)為上盤與下盤邊界頻譜較強、中間的頻譜較弱。斷層與古河道雖皆呈細長條帶狀,但古河道邊界頻譜較弱、中間頻譜較強??梢?,斷層與古河道的MEWVD頻譜異常雖有相似之處,卻也存在顯著差異。本質上,二者的共性與差異皆由地質環(huán)境決定,這些地震響應被MEWVD異常精確的刻畫了出來。

    總之,MEWVD頻譜聚焦性較強,時頻分辨率較高,且不同頻率的MEWVD對窄細不同的河道具有不同程度的敏感性,有效的揭示了中江氣田沙溪廟組微型古河道的寬度、砂巖厚度和流向等空間信息。

    4 結論

    熵最大時的頻譜特征,最能代表地震信號的頻譜分布狀態(tài),MEWVD不僅避免了 WVD產生的交叉項干擾,而且最大程度的增強了頻譜聚焦性,有效提升了地震信號的時頻分辨率。

    仿真地震信號和窄薄模型正演模擬信號試驗揭示,較STFT、GST、CWT、WVD、PWVD、SPWVD、CWD等方法,MEWVD能更加精確的提取信號頻譜特征;尤其是針對窄薄各異的數值模型,由于最大熵的差異而表現(xiàn)出了特殊的頻譜響應,采用不同頻率的MEWVD能夠準確的識別尺度各異的窄薄地質體。

    基于最大熵準則與 Wigner-Ville分布的微型古河道識別方法,充分發(fā)揮了MEWVD頻譜聚焦性較強、時頻分辨率較高等優(yōu)勢,有效的識別出了中江氣田不同尺度的微型古河道,且不同頻率的MEWVD對窄細各異的河道具有不同程度的敏感性,準確地刻畫了微型古河道的寬度、砂巖厚度和流向等空間信息。

    符號注釋:

    a——自回歸系數,又稱為預測誤差因子,無因次;am——序號為m的自回歸系數,無因次;1A——衰減分量的振幅,dB;A2——平穩(wěn)分量的振幅,dB;E——預測誤差,W/Hz;EB,m——序號為m的后向預測誤差,W/Hz;EF,m——序號為m的前向預測誤差,W/Hz;f——信號的頻率,Hz;f1——衰減分量的頻率,Hz;f2——平穩(wěn)分量的頻率,Hz;GR——自然伽馬,API;H——熵函數,無因次;i——虛數單位;j——自回歸系數am的序號,無因次;kn——自相關函數,無因次;kz——kn的卷積核,無因次;l——離散信號擴展采樣點序號,無因次;m——am、EF,m、EB,m、wn等函數的序號,無因次;M——AR模型的階數,無因次;n——離散信號采樣點序號,無因次;N——離散信號長度,無因次;p——分布密度函數,無因次;P——功率譜,W/Hz;P0——初始功率譜,W/Hz;Pm——平均功率譜,W/Hz;Q——吸收系數,無因次;Rlld——深側向電阻率,W·m;Rlls——淺側向電阻率,W·m;t——時間,s;W——連續(xù)WVD函數,W/Hz;Wz——離散WVD函數,W/Hz;wn——Wz的子項,W/Hz;X——信號,無因次;x——某時段信號,dB;x1——衰減分量,dB;x2——平穩(wěn)分量,dB;Dt——時間采樣率,s;xF——仿真信號,dB;j1——衰減分量的相位,rad;j2——平穩(wěn)分量的相位,rad;t——時延,s;z——離散信號,dB。

    猜你喜歡
    時頻振幅頻譜
    一種用于深空探測的Chirp變換頻譜分析儀設計與實現(xiàn)
    一種基于稀疏度估計的自適應壓縮頻譜感知算法
    測控技術(2018年7期)2018-12-09 08:58:22
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    認知無線電頻譜感知技術綜述
    基于時頻分析的逆合成孔徑雷達成像技術
    對采樣數據序列進行時頻分解法的改進
    雙線性時頻分布交叉項提取及損傷識別應用
    高清欧美精品videossex| 99在线视频只有这里精品首页| 91成人精品电影| 少妇 在线观看| 久久久久亚洲av毛片大全| 色在线成人网| 波多野结衣av一区二区av| 两个人免费观看高清视频| 精品久久久久久,| 日韩精品免费视频一区二区三区| 国产一区二区三区综合在线观看| 国产有黄有色有爽视频| 嫁个100分男人电影在线观看| 少妇裸体淫交视频免费看高清 | 99香蕉大伊视频| 免费不卡黄色视频| 一级毛片精品| 色哟哟哟哟哟哟| 久久国产精品男人的天堂亚洲| 久久精品国产清高在天天线| 亚洲国产精品一区二区三区在线| 免费人成视频x8x8入口观看| 色尼玛亚洲综合影院| 国产成人免费无遮挡视频| 欧美 亚洲 国产 日韩一| 亚洲中文字幕日韩| 国产蜜桃级精品一区二区三区| 久久天堂一区二区三区四区| 日韩视频一区二区在线观看| 日韩有码中文字幕| 可以在线观看毛片的网站| 亚洲专区国产一区二区| 久久精品aⅴ一区二区三区四区| 国产欧美日韩精品亚洲av| 久久久久久大精品| 欧美日韩亚洲高清精品| 久9热在线精品视频| 18禁裸乳无遮挡免费网站照片 | 99久久精品国产亚洲精品| 午夜精品在线福利| 高清毛片免费观看视频网站 | 搡老乐熟女国产| 免费在线观看日本一区| 欧美性长视频在线观看| 又紧又爽又黄一区二区| 三上悠亚av全集在线观看| 欧美日韩黄片免| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品999在线| 国产一区二区三区综合在线观看| 欧美成狂野欧美在线观看| 久久国产精品人妻蜜桃| 欧美一区二区精品小视频在线| 久久精品国产清高在天天线| 真人一进一出gif抽搐免费| 久久午夜综合久久蜜桃| 亚洲三区欧美一区| 亚洲第一青青草原| 国产精品美女特级片免费视频播放器 | 男女午夜视频在线观看| 国产精品一区二区三区四区久久 | 日本黄色视频三级网站网址| 欧美久久黑人一区二区| 人妻久久中文字幕网| 免费久久久久久久精品成人欧美视频| 亚洲成国产人片在线观看| 久久精品国产99精品国产亚洲性色 | 午夜两性在线视频| 超碰97精品在线观看| 久久人人97超碰香蕉20202| 午夜免费成人在线视频| 午夜影院日韩av| 国产亚洲精品第一综合不卡| 国产精品爽爽va在线观看网站 | 老司机午夜十八禁免费视频| 国产精品电影一区二区三区| 超碰97精品在线观看| av在线播放免费不卡| 午夜福利在线免费观看网站| 国产国语露脸激情在线看| svipshipincom国产片| 无限看片的www在线观看| 免费在线观看影片大全网站| 欧美黑人欧美精品刺激| 中亚洲国语对白在线视频| 国产精品亚洲av一区麻豆| 啦啦啦在线免费观看视频4| 一级a爱视频在线免费观看| 久久国产乱子伦精品免费另类| 欧美日韩中文字幕国产精品一区二区三区 | 精品人妻1区二区| 一二三四社区在线视频社区8| 亚洲黑人精品在线| 久久精品国产99精品国产亚洲性色 | 制服诱惑二区| 成人三级做爰电影| 色老头精品视频在线观看| 国产三级在线视频| 国产亚洲欧美在线一区二区| 国产成人av教育| 久久青草综合色| 老熟妇乱子伦视频在线观看| 一区二区日韩欧美中文字幕| 高清毛片免费观看视频网站 | 精品国产亚洲在线| 天堂√8在线中文| 大香蕉久久成人网| 亚洲情色 制服丝袜| 欧美激情久久久久久爽电影 | 村上凉子中文字幕在线| 久久久久久久精品吃奶| 伦理电影免费视频| 精品少妇一区二区三区视频日本电影| 国产欧美日韩一区二区三| 高清黄色对白视频在线免费看| 女警被强在线播放| 精品高清国产在线一区| 久久国产乱子伦精品免费另类| 国产亚洲欧美在线一区二区| 免费不卡黄色视频| 久久精品亚洲精品国产色婷小说| 村上凉子中文字幕在线| 久久人妻福利社区极品人妻图片| 亚洲成a人片在线一区二区| 在线国产一区二区在线| 啦啦啦免费观看视频1| 欧美日韩乱码在线| 亚洲一区高清亚洲精品| 99久久99久久久精品蜜桃| 国产野战对白在线观看| 人成视频在线观看免费观看| 美女扒开内裤让男人捅视频| 亚洲在线自拍视频| 日日摸夜夜添夜夜添小说| 99久久国产精品久久久| 热99re8久久精品国产| 人妻丰满熟妇av一区二区三区| 动漫黄色视频在线观看| 99精品欧美一区二区三区四区| 视频在线观看一区二区三区| 免费久久久久久久精品成人欧美视频| 免费日韩欧美在线观看| 巨乳人妻的诱惑在线观看| 久久人妻福利社区极品人妻图片| 99久久精品国产亚洲精品| 美国免费a级毛片| 久久国产乱子伦精品免费另类| 午夜福利在线免费观看网站| 亚洲avbb在线观看| 97人妻天天添夜夜摸| 国产主播在线观看一区二区| 两性夫妻黄色片| 日本vs欧美在线观看视频| 黄色毛片三级朝国网站| 天天躁狠狠躁夜夜躁狠狠躁| 久久亚洲真实| 成人精品一区二区免费| 亚洲国产精品sss在线观看 | 久久久精品国产亚洲av高清涩受| 男男h啪啪无遮挡| 国产成人欧美| 大陆偷拍与自拍| 亚洲欧美精品综合久久99| 日韩精品中文字幕看吧| 午夜福利,免费看| 日韩人妻精品一区2区三区| 免费在线观看亚洲国产| 亚洲色图综合在线观看| 亚洲午夜精品一区,二区,三区| 亚洲伊人色综图| 国产精品日韩av在线免费观看 | 女人被狂操c到高潮| 无限看片的www在线观看| 亚洲国产精品sss在线观看 | 国产单亲对白刺激| 国产亚洲精品综合一区在线观看 | 久久精品人人爽人人爽视色| 久久久精品国产亚洲av高清涩受| 亚洲一区二区三区不卡视频| 国产精品影院久久| 大型av网站在线播放| 日本三级黄在线观看| 国产精品1区2区在线观看.| 亚洲视频免费观看视频| 热99re8久久精品国产| 美国免费a级毛片| 国产在线观看jvid| 亚洲aⅴ乱码一区二区在线播放 | 久热这里只有精品99| 精品电影一区二区在线| 欧美+亚洲+日韩+国产| 亚洲精品久久成人aⅴ小说| 亚洲中文av在线| 丝袜美腿诱惑在线| 视频在线观看一区二区三区| 一级作爱视频免费观看| 成人18禁在线播放| 99国产精品99久久久久| 亚洲一区中文字幕在线| 午夜福利一区二区在线看| 免费搜索国产男女视频| 精品一区二区三区av网在线观看| 亚洲一区二区三区色噜噜 | 午夜精品国产一区二区电影| 国产无遮挡羞羞视频在线观看| 久久国产精品影院| 亚洲成人精品中文字幕电影 | 女人被躁到高潮嗷嗷叫费观| 婷婷精品国产亚洲av在线| 两个人免费观看高清视频| 天天影视国产精品| 免费av毛片视频| 午夜福利,免费看| 午夜福利免费观看在线| а√天堂www在线а√下载| 操出白浆在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲情色 制服丝袜| 国产成+人综合+亚洲专区| 91字幕亚洲| 叶爱在线成人免费视频播放| 亚洲在线自拍视频| av视频免费观看在线观看| 久久精品91无色码中文字幕| 国产成人精品久久二区二区免费| 日本欧美视频一区| www.www免费av| 久久国产亚洲av麻豆专区| 999久久久精品免费观看国产| 一二三四社区在线视频社区8| av福利片在线| 香蕉久久夜色| 人人妻人人爽人人添夜夜欢视频| 热99re8久久精品国产| 亚洲一区中文字幕在线| 人人妻人人爽人人添夜夜欢视频| 性少妇av在线| 亚洲人成网站在线播放欧美日韩| 美女午夜性视频免费| 亚洲五月婷婷丁香| 9色porny在线观看| 国产亚洲精品一区二区www| 成人特级黄色片久久久久久久| 亚洲av成人不卡在线观看播放网| 日韩精品免费视频一区二区三区| 久久久久国内视频| 中文字幕av电影在线播放| av天堂在线播放| 夜夜爽天天搞| 可以在线观看毛片的网站| 国产成人系列免费观看| 男人操女人黄网站| 国产亚洲精品一区二区www| 丝袜在线中文字幕| 一a级毛片在线观看| 长腿黑丝高跟| 国产成人精品在线电影| 免费观看人在逋| 69精品国产乱码久久久| 大码成人一级视频| 精品国产乱子伦一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 久久香蕉精品热| 日韩欧美一区视频在线观看| 丰满饥渴人妻一区二区三| 亚洲国产精品合色在线| 婷婷精品国产亚洲av在线| 国产成人免费无遮挡视频| 久久久久久亚洲精品国产蜜桃av| 国产成人免费无遮挡视频| 亚洲av成人一区二区三| 精品福利观看| 午夜a级毛片| 老鸭窝网址在线观看| 校园春色视频在线观看| 热99国产精品久久久久久7| 成人18禁高潮啪啪吃奶动态图| 人成视频在线观看免费观看| 91av网站免费观看| 久久久久九九精品影院| 电影成人av| 亚洲狠狠婷婷综合久久图片| 国产熟女午夜一区二区三区| 伦理电影免费视频| 一区二区三区国产精品乱码| 夜夜看夜夜爽夜夜摸 | 自线自在国产av| 国产一区在线观看成人免费| 久久人妻福利社区极品人妻图片| 十八禁人妻一区二区| 久久影院123| 日韩国内少妇激情av| 国产精品久久久av美女十八| 精品久久久久久久毛片微露脸| 99国产综合亚洲精品| 国产精品av久久久久免费| x7x7x7水蜜桃| 日本a在线网址| 国产极品粉嫩免费观看在线| xxx96com| 9热在线视频观看99| 黑人巨大精品欧美一区二区蜜桃| 91成人精品电影| 亚洲第一青青草原| 国产精品久久视频播放| 正在播放国产对白刺激| 亚洲专区字幕在线| 自线自在国产av| 国产精品久久久人人做人人爽| av有码第一页| 啪啪无遮挡十八禁网站| av片东京热男人的天堂| 亚洲 欧美一区二区三区| 热re99久久精品国产66热6| 久热这里只有精品99| 一级毛片女人18水好多| 精品国产国语对白av| 国产免费男女视频| 国产亚洲精品第一综合不卡| 亚洲七黄色美女视频| 久久精品亚洲av国产电影网| 欧洲精品卡2卡3卡4卡5卡区| 国产精品乱码一区二三区的特点 | www.www免费av| 国产高清国产精品国产三级| 欧美乱色亚洲激情| 欧美亚洲日本最大视频资源| 亚洲国产毛片av蜜桃av| 日本五十路高清| 午夜亚洲福利在线播放| 国产精品秋霞免费鲁丝片| 美女 人体艺术 gogo| 免费在线观看日本一区| 欧美精品亚洲一区二区| 99久久久亚洲精品蜜臀av| 黄片小视频在线播放| 久久中文字幕一级| 午夜福利在线观看吧| 97超级碰碰碰精品色视频在线观看| 99香蕉大伊视频| 亚洲av成人一区二区三| 精品国产超薄肉色丝袜足j| 久久狼人影院| 亚洲全国av大片| 新久久久久国产一级毛片| 他把我摸到了高潮在线观看| 免费人成视频x8x8入口观看| 久久久久精品国产欧美久久久| 激情在线观看视频在线高清| 欧美av亚洲av综合av国产av| 欧美成人性av电影在线观看| 老司机福利观看| 老熟妇乱子伦视频在线观看| 两个人免费观看高清视频| 亚洲一码二码三码区别大吗| 搡老乐熟女国产| 天堂俺去俺来也www色官网| 国产欧美日韩一区二区三| 久9热在线精品视频| netflix在线观看网站| 又紧又爽又黄一区二区| 午夜免费成人在线视频| 亚洲精品美女久久久久99蜜臀| 久久中文字幕人妻熟女| 久久九九热精品免费| 国产精品 欧美亚洲| 亚洲专区国产一区二区| 日韩精品免费视频一区二区三区| 国产黄a三级三级三级人| 操出白浆在线播放| 免费观看精品视频网站| 韩国精品一区二区三区| 在线天堂中文资源库| 成年人黄色毛片网站| 免费在线观看亚洲国产| 国产成人精品久久二区二区免费| www日本在线高清视频| 亚洲精品中文字幕在线视频| 999久久久国产精品视频| 十八禁网站免费在线| 99国产精品一区二区蜜桃av| 丰满人妻熟妇乱又伦精品不卡| 韩国av一区二区三区四区| 日韩免费av在线播放| 淫秽高清视频在线观看| 成人国语在线视频| 国产精品久久久久成人av| 好男人电影高清在线观看| 一二三四社区在线视频社区8| 老鸭窝网址在线观看| 十分钟在线观看高清视频www| 国产在线精品亚洲第一网站| 亚洲成人免费电影在线观看| 大型黄色视频在线免费观看| 国产免费现黄频在线看| 欧美日韩亚洲国产一区二区在线观看| ponron亚洲| 午夜精品在线福利| av免费在线观看网站| 男女之事视频高清在线观看| 国产成+人综合+亚洲专区| 精品电影一区二区在线| 黄色视频不卡| 在线观看免费日韩欧美大片| 久久久久国产一级毛片高清牌| 天天躁狠狠躁夜夜躁狠狠躁| 精品卡一卡二卡四卡免费| 精品第一国产精品| 久久婷婷成人综合色麻豆| 国产欧美日韩一区二区三| 老鸭窝网址在线观看| 50天的宝宝边吃奶边哭怎么回事| 女同久久另类99精品国产91| av在线播放免费不卡| 叶爱在线成人免费视频播放| 国产成人一区二区三区免费视频网站| 91麻豆av在线| 中文字幕精品免费在线观看视频| 国产深夜福利视频在线观看| 99国产综合亚洲精品| 婷婷丁香在线五月| 黄色成人免费大全| 99re在线观看精品视频| 国产精品二区激情视频| 丝袜美腿诱惑在线| 桃色一区二区三区在线观看| 欧美国产精品va在线观看不卡| 9色porny在线观看| 侵犯人妻中文字幕一二三四区| 美女高潮到喷水免费观看| 亚洲男人的天堂狠狠| 老熟妇仑乱视频hdxx| a在线观看视频网站| 日韩国内少妇激情av| 国产午夜精品久久久久久| 国产精品一区二区在线不卡| av天堂在线播放| 欧美人与性动交α欧美精品济南到| 热99国产精品久久久久久7| 国产精品日韩av在线免费观看 | 脱女人内裤的视频| 精品乱码久久久久久99久播| 高清av免费在线| 亚洲av成人不卡在线观看播放网| 亚洲黑人精品在线| 麻豆一二三区av精品| 亚洲少妇的诱惑av| 国产精品久久久久久人妻精品电影| 黑人欧美特级aaaaaa片| 久久久国产欧美日韩av| 69精品国产乱码久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 免费日韩欧美在线观看| 在线天堂中文资源库| 9色porny在线观看| 一区福利在线观看| 精品日产1卡2卡| 亚洲美女黄片视频| 久久久精品国产亚洲av高清涩受| 亚洲人成网站在线播放欧美日韩| 精品国产乱码久久久久久男人| www.熟女人妻精品国产| 老司机靠b影院| 俄罗斯特黄特色一大片| 成人18禁在线播放| 亚洲一区二区三区色噜噜 | 一边摸一边抽搐一进一小说| 日韩免费高清中文字幕av| 99热只有精品国产| 婷婷丁香在线五月| 午夜精品国产一区二区电影| 久久人妻福利社区极品人妻图片| 久久久久精品国产欧美久久久| 丝袜在线中文字幕| 啪啪无遮挡十八禁网站| 国产激情欧美一区二区| 少妇被粗大的猛进出69影院| 午夜老司机福利片| 自线自在国产av| 高清毛片免费观看视频网站 | 99久久久亚洲精品蜜臀av| 亚洲国产欧美一区二区综合| 亚洲 欧美一区二区三区| 久久精品aⅴ一区二区三区四区| 女警被强在线播放| 妹子高潮喷水视频| 日日干狠狠操夜夜爽| 午夜福利影视在线免费观看| 国产一区二区三区综合在线观看| 久久伊人香网站| 国产真人三级小视频在线观看| 国产精品秋霞免费鲁丝片| 亚洲精品在线美女| 久久久久国内视频| 啦啦啦在线免费观看视频4| 欧美最黄视频在线播放免费 | 99re在线观看精品视频| 国产精品秋霞免费鲁丝片| 亚洲欧美日韩另类电影网站| 欧美日韩亚洲综合一区二区三区_| 欧美大码av| av视频免费观看在线观看| 亚洲一区高清亚洲精品| 国产野战对白在线观看| 日本一区二区免费在线视频| 亚洲中文字幕日韩| 露出奶头的视频| 看免费av毛片| 黄片小视频在线播放| 伦理电影免费视频| 亚洲成人精品中文字幕电影 | 国产精品 欧美亚洲| 亚洲片人在线观看| 91大片在线观看| 嫩草影院精品99| 在线观看66精品国产| bbb黄色大片| 男女下面进入的视频免费午夜 | 精品国产一区二区三区四区第35| 欧美精品一区二区免费开放| www日本在线高清视频| 88av欧美| www.精华液| 亚洲自拍偷在线| 精品一区二区三区四区五区乱码| www日本在线高清视频| 国内毛片毛片毛片毛片毛片| 咕卡用的链子| 最近最新免费中文字幕在线| 亚洲 欧美 日韩 在线 免费| 一进一出抽搐动态| 欧美日韩亚洲高清精品| 久99久视频精品免费| 黄网站色视频无遮挡免费观看| 露出奶头的视频| 成人18禁高潮啪啪吃奶动态图| 精品一区二区三卡| 国产97色在线日韩免费| 精品人妻1区二区| 精品国产一区二区久久| 国产99白浆流出| 一级,二级,三级黄色视频| 黄网站色视频无遮挡免费观看| tocl精华| 天天躁夜夜躁狠狠躁躁| 97超级碰碰碰精品色视频在线观看| 天堂俺去俺来也www色官网| 免费搜索国产男女视频| 美国免费a级毛片| 国产一区二区三区在线臀色熟女 | 自线自在国产av| 国产精品成人在线| 亚洲熟女毛片儿| 两性夫妻黄色片| 好男人电影高清在线观看| 黑人猛操日本美女一级片| 欧美老熟妇乱子伦牲交| 成人影院久久| 美女大奶头视频| www国产在线视频色| 国产男靠女视频免费网站| 99久久国产精品久久久| 亚洲全国av大片| 欧美日韩福利视频一区二区| 国产无遮挡羞羞视频在线观看| 久久久久久久精品吃奶| 美女国产高潮福利片在线看| 午夜亚洲福利在线播放| 欧美不卡视频在线免费观看 | 成年人免费黄色播放视频| 亚洲成a人片在线一区二区| 侵犯人妻中文字幕一二三四区| 18禁美女被吸乳视频| 91在线观看av| 视频区欧美日本亚洲| 国产亚洲精品综合一区在线观看 | 在线天堂中文资源库| 亚洲国产欧美网| 国产亚洲精品久久久久5区| 高清欧美精品videossex| 国产伦人伦偷精品视频| a级毛片黄视频| 69av精品久久久久久| 亚洲av五月六月丁香网| 亚洲国产欧美一区二区综合| 国产精品99久久99久久久不卡| 午夜亚洲福利在线播放| 国产免费男女视频| cao死你这个sao货| 女性被躁到高潮视频| 亚洲人成77777在线视频| 麻豆一二三区av精品| 午夜福利在线免费观看网站| 99国产精品一区二区三区| 1024视频免费在线观看| 不卡一级毛片| 久久久国产成人免费| 亚洲精品国产区一区二| 欧美精品一区二区免费开放| 亚洲色图av天堂| 亚洲欧美一区二区三区黑人| 757午夜福利合集在线观看| 亚洲一码二码三码区别大吗| 亚洲午夜精品一区,二区,三区| 亚洲九九香蕉| 97人妻天天添夜夜摸| 国产精品偷伦视频观看了| 身体一侧抽搐| 女人精品久久久久毛片| 丁香欧美五月| 99久久综合精品五月天人人| 婷婷六月久久综合丁香| 久久影院123| 久久天堂一区二区三区四区| 淫妇啪啪啪对白视频|