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

    基于經(jīng)驗小波變換的瞬時譜分析技術及其在致密砂巖氣藏檢測中的應用

    2023-09-23 08:15:42陳偉楊浪陳名德連辰浩江源
    長江大學學報(自科版) 2023年4期
    關鍵詞:子波譜分析時頻

    陳偉,楊浪,陳名德,連辰浩,江源

    1.油氣資源與勘探技術教育部重點實驗室(長江大學),湖北 武漢 430100 2.長江大學地球物理與石油資源學院,湖北 武漢 430100

    隨著全球原油需求日益增長,油田勘探程度不斷提升,促使勘探重點轉向隱蔽油氣藏,該類油氣藏受埋藏深度影響,多由致密砂巖組成,具有超高壓、超深、低孔、低滲等特點,這對勘探技術帶來新的挑戰(zhàn),而瞬時譜分析技術給油氣勘探向更深處發(fā)展帶來希望。地震信號是非線性的、非平穩(wěn)的信號,它們的能量分布可以通過時頻分析獲得,因此可以在時域和頻域,即聯(lián)合時頻分布中完全建立信號能量或強度的分布。時頻分析方法可以觀察時域和頻域信號的演化,并顯示信號的局部時頻特征。因此,時頻分析方法已被廣泛應用于地震波分析中。瞬時譜分析技術應用到地球物理領域中檢測油氣藏的存在現(xiàn)已成為熱點。瞬時譜分析技術主要依賴于時頻分析方法的優(yōu)劣,如今廣泛應用的時頻分析方法主要分為兩類,一類為線性時頻分析方法,如Gabor變換,小波變換,S變換等;另一類為非線性時頻分析方法,如Hilbert-Huang變換等。

    Fourier變換只能將時間域的信號整體變換到頻率域,將信號的頻率成分和分布表示出來,雖然使信號擁有了頻率分辨率,卻丟失了時間分辨率,不能反映出信號的頻率特征隨時間的變化情況。為了解決該問題,GABOR[1]于1946年首次提出在時頻域內(nèi)分解一維信號的方法,即短時Fourier變換。短時Fourier變換中最具代表的為Gabor變換。Gabor變換選取高斯函數(shù)作為窗函數(shù),在信號的每一小段時間間隔[2]內(nèi)進行Fourier變換,該窗口在時間軸上進行平移,對每一段信號進行局部Fourier分析后,再將變換得到的頻率信息放到時間-頻率的二維平面上,從而得到時頻譜。Gabor變換選取高斯函數(shù)作為窗函數(shù)的優(yōu)點在于高斯函數(shù)的Fourier變換還是高斯函數(shù),不會對變換結果產(chǎn)生不良影響,可以同時提供時域和頻域局部化的信息。另外,根據(jù)Heisenberg測不準原理,變換得到的結果不能同時擁有極高的時間分辨率和頻率分辨率,極大地影響了分析結果的好壞。而高斯函數(shù)窗口面積已達到Heisenberg測不準原理的下界,能兼顧時間軸和頻率軸的分辨率。Gabor變換的基函數(shù)不能成為正交基[3],為保證不丟失信息,只能采用非正交基參與信號分析,使得計算復雜度更高。Gabor變換的窗口固定不變,無法針對高低頻信號給出可以調(diào)節(jié)的時頻窗,割斷了頻率與窗口長度的內(nèi)在聯(lián)系[4];對短時突發(fā)信號進行分析時,可能產(chǎn)生時間或頻率上的模糊現(xiàn)象[5]。MORLET[6]等人從短時Fourier變換中推出小波變換,能對信號同時在時間域和頻率域內(nèi)進行局部化分析。與Fourier變換不同,小波變換有兩個變量:尺度和平移量。尺度控制小波函數(shù)的伸縮,平移量控制小波函數(shù)的平移。尺度就對應于頻率(反比),平移量就對應于時間。尺度增大,則時窗伸展,頻寬收縮,帶寬變窄,中心頻率降低,時間分辨率降低而頻率分辨率增高;尺度減小,則相反[7]。這符合實際問題中高頻信號持續(xù)時間短、低頻信號持續(xù)時間長的自然規(guī)律。小波變換能準確反映待測信號的幅頻特性。在小波變換中,時頻窗口的長度和寬度隨著參數(shù)尺度的變化而變化[8]。雖然小波變換克服了短時Fourier變換單一分辨率分析的不足,但是由于引入的尺度因子與頻率沒有直接的聯(lián)系,且在小波變換中沒有明顯表現(xiàn)出來,因此小波變換的結果不是一種真正的時頻譜。小波分析的另一個問題是其具有的自適應特點,一旦基本小波被選定,就必須用它來分析所有待分析的數(shù)據(jù)。STOCKWELL等[9]結合短時Fourier變換和小波變換提出S變換,S變換相對于小波變換的不同在于其公式中加入了相位項,用于進行相位校正。它與小波變換一樣存在基本小波,其基本小波由簡諧波函數(shù)與高斯函數(shù)相乘組成,其中,簡諧波函數(shù)在時間域內(nèi)作伸縮變換,而高斯函數(shù)則作伸縮和平移變換[10]。S變換既克服了短時Fourier變換和小波變換的不足,也繼承了它們的優(yōu)點[11-12]。但是S變換的窗函數(shù)隨頻率以固定的趨勢改變,基本小波一經(jīng)選定則無法更改。PINNEGAR等[13]、高靜懷等[14]在S變換的基礎上提出廣義S變換,通過各種可調(diào)的屬性參數(shù)控制窗函數(shù),使其具有更高的適應性。

    上述方法都是在Fourier變換的基礎上發(fā)展起來的,但是Fourier變換存在時間域和頻率域的相互制約,無法兼顧時間軸和頻率軸的分辨率,并且在處理非線性非平穩(wěn)信號時效果較差。1998年,HUANG等[15]首次提出希爾伯特-黃變換,這是一種用于處理非線性、非平穩(wěn)信號的時頻分析新方法,該方法由經(jīng)驗模態(tài)分解[16]和Hilbert譜分析兩部分組成。經(jīng)驗模態(tài)分解就是將非平穩(wěn)信號按頻率由高到低分解為一系列固有模態(tài)函數(shù),這一系列固有模態(tài)分量能反映出原始信號在不同頻率上的局部特征。對每個固有模態(tài)分量進行Hilbert變換,得到每一個固有模態(tài)分量隨時間變換的瞬時頻率和瞬時幅值,求得振幅-頻率-時間的三維譜分布[17]。把所有固有模態(tài)分量的Hilbert譜進行綜合之后即可求出原始信號的Hilbert譜,以此來得到信號的時頻屬性。希爾伯特-黃變換很適合處理非線性非平穩(wěn)信號,在地震資料處理和瞬時譜分析中有著廣泛應用。然而經(jīng)驗模態(tài)分解存在模態(tài)混疊現(xiàn)象。模態(tài)混疊主要是指同一個固有模態(tài)分量當中出現(xiàn)了不同尺度或頻率的信號,或者同一尺度或頻率的信號被分解到多個不同的固有模態(tài)分量當中[18],嚴重影響了數(shù)據(jù)分析的結果。為了克服該問題,學者提出了總體經(jīng)驗模態(tài)分解算法[19]。總體經(jīng)驗模態(tài)分解方法是一種噪音分析法,即在每次信號分解的過程中對原始信號添加高斯白噪聲,該方法能夠很好地解決經(jīng)驗模態(tài)分解算法在分解過程中出現(xiàn)的端點效應[20]和模態(tài)混疊現(xiàn)象??傮w經(jīng)驗模態(tài)分解方法在原始信號中加入了高斯白噪聲,所以分解后的結果很難重構得到原始信號,分解不具備很好的完備性。其分解依賴添加白噪聲幅值和集成次數(shù),如果參數(shù)選擇不合適,不僅不能抑制模態(tài)混淆,而且會出現(xiàn)偽分量;且也無法保證分解得到的分量滿足固有模態(tài)分量的定義條件[21]。總體經(jīng)驗模態(tài)分解方法需要的計算量非常大,計算效率較低,不利于實際地震資料的處理。TORRES等[22]在2010年提出完備總體經(jīng)驗模態(tài)分解,與總體經(jīng)驗模態(tài)分解不同的是,該方法向原始信號中加入了正負成對形式的白噪音,這樣重構信號中的殘余輔助信號被有效消除,從而簡化計算量。完備總體經(jīng)驗模態(tài)分解既繼承了經(jīng)驗模態(tài)分解的正交性、自適應性以及能夠處理非平穩(wěn)信號的特點,又在總體經(jīng)驗模態(tài)分解的基礎上簡化了計算量,幾乎可以重構完整的原始信號。

    經(jīng)驗小波變換作為一種全新的自適應性算法在信號處理中得到應用。經(jīng)驗小波變換能夠更好地分解出原始信號中固有的本征信號,因此其較經(jīng)驗模態(tài)分解和總體經(jīng)驗模態(tài)分解算法有更高的自適應性。經(jīng)驗小波變換以成熟的小波理論為基礎,其數(shù)學理論基礎非常充分,自身具有較高的計算效率??紤]到經(jīng)驗小波變換是基于經(jīng)驗模態(tài)分解的最新自適應算法,本文先介紹了經(jīng)驗模態(tài)分解和經(jīng)驗小波變換的基本原理,然后提出基于經(jīng)驗小波變換的瞬時譜分析技術,最后將其應用到致密砂巖氣藏的檢測中。本文首次將經(jīng)驗小波變換應用到瞬時譜分析技術中檢測致密砂巖氣藏,得到了較好的效果。

    1 基本原理

    1.1 經(jīng)驗模態(tài)分解

    HUANG等[15]人認為,只有固有模態(tài)分量的瞬時頻率具有物理意義。經(jīng)驗模態(tài)分解基于簡單的假設,即任何數(shù)據(jù)都包含不同的簡單固有模式。這些振蕩模式中的每一種都由固有模態(tài)分量表示,具有以下定義:①在整個信號段中,極值的數(shù)量和過零的數(shù)量必須相等或相差最多一個;②在任何時候,由局部最大值定義的包絡的平均值和由局部最小值定義的包絡為零。固有模態(tài)函數(shù)定義的第一個條件類似于靜態(tài)高斯過程的窄帶要求。第二個條件使用局部要求而非全局要求,這使得由對稱波形引起的不需要的波動不會出現(xiàn)在瞬時頻率中。

    固有模態(tài)函數(shù)的定義表明,數(shù)據(jù)的內(nèi)部振動模式由固有模態(tài)分量表征。每個固有模態(tài)分量只涉及一種振蕩模式。振蕩也將相對于“局部均值”對稱。線性或非線性的固有模態(tài)分量可以在簡諧分量中具有恒定的幅度和頻率,也可以在時間函數(shù)中具有可變幅度和頻率。經(jīng)驗模態(tài)分解通過將多組分信號簡化為單組分函數(shù)的集合[23],解決了從多組分信號中計算有意義瞬時頻率的難題。最終的復雜信號是由固有模態(tài)分量的重疊形成的。經(jīng)驗模態(tài)分解的目的是獲得固有模態(tài)分量。通過經(jīng)驗模態(tài)分解方法分解的固有模態(tài)分量,可以通過篩選過程根據(jù)信號本身的相鄰極值點之間的延遲來定義和區(qū)分。

    HUANG等[24]表明,為了獲得關于經(jīng)驗模態(tài)分解詳細而經(jīng)驗性的統(tǒng)計知識,經(jīng)驗模態(tài)分解基本上充當了與小波分解相關的二元濾波器組[25],經(jīng)驗模態(tài)分解方法將信號分解為固有模態(tài)分量的步驟如下:

    1)識別時間序列信號X(t)的所有極值點,使用X(t)的所有極大點確定上包絡u(t),并用所有的極小點確定下包絡v(t),X(t)滿足:

    v(t)≤X(t)≤u(t)

    (1)

    然后,上包絡線和下包絡線的平均曲線m(t)是:

    (2)

    設h1(t)=X(t)-m(t),那么h1(t)就是固有模態(tài)分量。

    2)由于過沖和俯沖,包絡樣條將近似產(chǎn)生一個新的極值,它會影響原始極值點的位置和大小。因此,h1(t)不完全滿足固有模態(tài)分量條件。為了獲得所需的h1(t),讓h1(t)代替X(t)。對應于h1(t),上包絡線為u1(t),下包絡線為v1(t),然后重復這個過程:

    (3)

    h2(t)=h1(t)-m1(t)

    (4)

    ?

    (5)

    hk(t)=hk-1-mk-1(t)

    (6)

    重復該過程直到每個hk(t)滿足固有模態(tài)分量條件。然后,獲得第一個固有模態(tài)分量C1(t)和信號r1(t)的剩余部分:

    C1(t)=hk(t)

    (7)

    r1(t)=X(t)-C1(t)

    (8)

    在信號的剩余部分繼續(xù)使用經(jīng)驗模態(tài)分解。分解繼續(xù)進行,直到信號的剩余部分單調(diào)或其值小于預定值。分解后得到的固有模態(tài)分量和余量為:

    r2(t)=r1(t)-C2(t)

    ?

    rn(t)=rn-1(t)-Cn(t)

    (9)

    X(t)可以表示為固有模態(tài)分量和余量之和:

    X(t)=C1(t)+C2(t)+…+Cn(t)+rn(t)

    (10)

    1.2 經(jīng)驗小波變換

    經(jīng)驗小波變換的目標是通過構建自適應小波來提取不同的單分量。該方法的步驟如下:

    1)將快速Fourier變換應用于信號f(t),從而獲得頻譜X(w),其中f(t)是離散信號,t={ti}i=1,2,…,M(M表示樣本數(shù))。在Fourier頻譜中找出最大值M={Mi}1,2,…,N(N表示最大值的數(shù)量),并推導它們的相應頻率W={Wi}i=1,2,…,N。

    2)獲得Fourier頻譜和邊界集的適當分段。將每個段的邊界Ωi定義為兩個連續(xù)最大值的中心:

    (11)

    式中:Wi和Wi+1為兩個頻率;邊界集是Ω={Ωi}i=1,2,…,N-1。

    3)定義一組由一個低通濾波器和基于邊界的N-1個帶通濾波器組成N個小波濾波器。尺度函數(shù)φ1(W)和經(jīng)驗小波ψi(W)的Fourier變換的表達式如下:

    (12)

    (13)

    (14)

    求取尺度函數(shù)和小波函數(shù)以提取不同單分量。近似系數(shù)由分析信號f(t)的內(nèi)積與經(jīng)驗尺度函數(shù)表示:

    (15)

    類似地,細節(jié)系數(shù)是通過經(jīng)驗小波的分析信號f(t)的內(nèi)積得到的:

    (16)

    式中:Wf(i,t)表示第t個時間點的第i個濾波器組的細節(jié)系數(shù)。

    怎樣將Fourier譜進行分段在經(jīng)驗小波變換中至關重要,其直接關系到對原始信號進行分解后的自適應程度。經(jīng)驗小波變換將原始信號進行不同的分割,比如對某個中心頻率的緊支撐部分進行分割。假設斷點數(shù)目為N,意味著需要N+1個邊界。除了起點0和終點σ以外,還需要N-1個邊界。為了找到這些邊界,首先對信號頻譜的局部極大值點進行降序排列(起點0和終點σ包括在內(nèi))。假設找到了M個極大值點,下面兩種情況將會出現(xiàn):①M≥N:算法發(fā)現(xiàn)了足夠的極值點以便于分割原始信號,但只取前N-1個極大值點;②M≤N:信號沒有預期的那么多模態(tài),將這M個極值點保留,并添加一些近似值直到極值點達到N個。

    (17)

    (18)

    這樣原始信號可以通過下面的式子來重構:

    (19)

    為簡單起見,經(jīng)驗小波變換所蘊含的經(jīng)驗模態(tài)函數(shù)可以定義為:

    (20)

    (21)

    1.3 基于經(jīng)驗小波變換的瞬時譜分析技術

    瞬時譜分析(instantaneous spectral analysis)能夠對地震信號進行連續(xù)時頻分析,觀察任一時間點處振幅變化情況及頻率分布。

    基于經(jīng)驗小波變換的瞬時譜分析技術實現(xiàn)步驟為:①利用經(jīng)驗小波變換將某一單道地震信號分解成按頻率高低分布的固有模態(tài)分量;②求取各個固有模態(tài)分量中的瞬時頻率分布;③將各個固有模態(tài)分量的頻率分布依據(jù)時間疊加得到該地震道的時頻剖面;④依次選取二維地震記錄中每個單道記錄,重復進行步驟①~③;⑤從各時頻剖面中提取特定頻率分量,根據(jù)對應的地震道位置依次疊加,最終得到特定頻率的瞬時譜。

    評判瞬時譜分析剖面時遵循以下原則:①時頻分布沿著頻率的振幅疊加值與信號的瞬時振幅值近似相等;②時頻分布沿著時間的振幅疊加值與信號的瞬時頻率值近似相等;③地震剖面上明顯的構造,在瞬時譜剖面上也必須表現(xiàn)出來;④地震剖面上明顯的構造異常旁瓣,在瞬時譜剖面上不能單獨出現(xiàn);⑤一個單獨的同相軸在瞬時譜剖面上是連續(xù)的。

    2 合成數(shù)據(jù)

    為驗證本文所提方法的有效性,現(xiàn)構造一個典型的非平穩(wěn)信號,分別利用經(jīng)驗小波變換和其他時頻分析方法對其進行對比,如短時Fourier變換,S變換,經(jīng)驗模態(tài)分解。

    該信號具體形態(tài)(見圖1)為20 Hz背景余弦波,在300 ms處疊加100 Hz Morlet子波,在1 070 ms和1 100 ms處疊加了兩個30 Hz Richer子波。在1 300~1 700 ms之間,信號具有3種不同的頻率組分:7,30和40 Hz,其中7 Hz的頻率是不連續(xù)的,并且包含少于一個周期的部分,出現(xiàn)在1.37,1.51和1.65 s。該信號具有顯著特征,為典型的非平穩(wěn)信號,適合用于檢驗時頻分析方法的優(yōu)劣。

    圖1 測試信號Fig.1 Test signal

    2.1 信號分解對比

    經(jīng)驗模態(tài)分解將原始信號分解成了7個固有模態(tài)分量(見圖2(a)),在300 ms處信號頻率大幅突變,導致固有模態(tài)分量1(IMF1)中雖提取出了Morlet子波的高頻成分,但混雜了20 Hz背景余弦波、1 100 ms處的Ricker子波等低頻成分。同樣,固有模態(tài)分量2(IMF2)、固有模態(tài)分量3(IMF3)、固有模態(tài)分量4(IMF4)中均是各種組分的高頻部分和低頻部分混雜在同一固有模態(tài)分量中,難以識別出信號包含的具體頻率組分,使信號分析復雜化,這屬于典型的模態(tài)混疊現(xiàn)象。圖2(b)是經(jīng)驗小波變換分解原始信號的結果,IMF1表示1 300~1 700 ms之間的低頻組分(7 Hz)。IMF2中包含20 Hz背景余弦波、1 300~1 700 ms處低頻組分及1 070 ms和1 100 ms處兩個Ricker子波。而IMF4僅提取了300 ms處的Morlet子波,沒有其他組分混雜。IMF3則為少量Morlet子波剩余組分和Ricker子波剩余組分。經(jīng)驗小波變換分解效果明顯優(yōu)于經(jīng)驗模態(tài)分解分解效果。

    注:IMF為固有模態(tài)分量。圖2 不同方法的信號分解結果對比Fig.2 Comparison of signal decomposition results of different methods

    2.2 時頻譜對比

    圖3(a)是經(jīng)過128ms時間窗的短時Fourier變換得到的時頻譜,30 Hz背景余弦波、300 ms處的100 Hz Morlet子波、1 300~1 700 ms處的低頻部分的輪廓已凸顯出來,但時頻分辨率過低,能量顯示較為分散,只能觀察到各頻率組分的大致分布情況,受Heisenberg測不準原理制約,缺乏實際研究價值。

    圖3 不同方法處理原始信號得到的時頻譜對比Fig.3 Comparison of time spectrum obtained by processing raw signals with different methods

    圖3(b)是通過S變換得到的時頻譜,在300 ms處的頻率分辨率很低,無法識別出100 Hz Morlet子波,1 300~1 700 ms處的時頻分辨率較短時Fourier變換相比有明顯提升,能量顯示更為集中。

    圖3(c)顯示通過經(jīng)驗模態(tài)分解得到的瞬時譜,圖中30 Hz背景余弦波,1 300~1 700 ms處低頻組分以及1 100 ms附近的兩個Ricker子波均清晰可見,時頻分辨率較前兩種方法有質的飛躍。在信號兩端可見明顯的發(fā)散現(xiàn)象,屬于典型的端點效應,導致200~1 800 ms間有大量低頻成分混雜,特別是在300 ms處100 Hz Morlet子波部分,其上下均有能量泄露,信號失真十分嚴重。

    經(jīng)驗小波變換處理后的瞬時譜(見圖3(d))首次將1 300~1 700 ms處的7 Hz頻率特征準確地刻畫出來。30 Hz背景余弦波、Ricker子波以及1 300~1 700 ms處其他頻率組分均以最高時頻分辨率呈現(xiàn)出來,且300 ms處高頻突變部分標識準確,能量非常集中,信號兩端未出現(xiàn)端點效應。對比上述4種時頻分析方法處理結果,可以發(fā)現(xiàn)經(jīng)驗小波變換明顯優(yōu)于其他3種方法。

    3 實際地震資料處理

    為了驗證本文方法在含氣檢測中的有效性,針對中國某地區(qū)的實際資料進行瞬時譜分析。該區(qū)儲集層物性總體較差,屬致密-超致密砂巖區(qū)塊。該區(qū)天然氣的富集、產(chǎn)出和高產(chǎn)與裂縫的發(fā)育程度關系極為密切??紤]到裂縫的存在以及天然氣的影響使地震波速度明顯降低,造成物性參數(shù)轉換上具有不同程度的突變特征,因此該地區(qū)的地震資料具有非線性非平穩(wěn)的特性,適合利用本文的經(jīng)驗小波變換對該區(qū)儲集層的含氣特征進行研究。

    取實際地震資料的部分疊后剖面(1 550至2 000道中1 500~2 000 ms,見圖4)進行應用。地震資料經(jīng)過瞬時譜分析技術處理后常出現(xiàn)含氣區(qū)下方振幅異常的低頻陰影現(xiàn)象,低頻陰影現(xiàn)象是含氣檢測的主要依據(jù)。用基于經(jīng)驗小波變換的瞬時譜分析技術處理地震資料,提取3~7 Hz低頻、13~17 Hz中頻以及23~27 Hz高頻的瞬時譜剖面,結果如圖5所示。圖5(a)3~7 Hz低頻瞬時譜剖面,在含氣儲層中信號能量迅速衰減,在其下方(紅色橢圓內(nèi))存在強大的能量分布,出現(xiàn)低頻陰影現(xiàn)象;圖5(c)23~27 Hz高頻瞬時譜剖面,當頻率增大時,含氣儲層中能量增大,其下方的低頻陰影現(xiàn)象消失。結合低頻和高頻瞬時譜剖面,對比能量分布差異可以很好地確定含氣區(qū)域。

    圖4 中國某地區(qū)實際資料疊后地震剖面Fig.4 Post-stack seismic profile of actual data in a region of China

    注:黃色橢圓表示儲層位置,紅色橢圓表示低頻陰影處,紅色豎線表示A井所在位置。圖5 基于經(jīng)驗小波變換的瞬時譜分析技術處理得到的瞬時譜剖面Fig.5 Instantaneous spectrum section processed by instantaneous spectrum analysis technique based on empirical wavelet transform

    4 結束語

    經(jīng)驗小波變換是一種全新的自適應分析算法,具有很高的自適應性,其時頻分析解決了模態(tài)混疊問題,能夠有效地去除線性噪音,也不易丟失數(shù)據(jù)。本文將經(jīng)驗小波變換應用到瞬時譜分析技術中,得到了一種新的瞬時譜分析技術,與基于傳統(tǒng)時頻分析方法的瞬時譜分析技術相比解決了Heisenberg測不準原理的制約等難題,具有更高的時間和頻率分辨率。本文將基于經(jīng)驗小波變換的瞬時譜分析技術應用于實際致密砂巖氣藏資料的含氣檢測中,一維地震道分析效果優(yōu)于其他方法,二維地震道分析中低頻陰影更加明顯,低頻與高頻的瞬時譜剖面結合能夠更好地反映地層細節(jié),更有效地識別致密砂巖氣藏的含氣巖層。

    猜你喜歡
    子波譜分析時頻
    納譜分析技術(蘇州)有限公司
    色譜(2022年5期)2022-04-28 02:49:10
    一類非線性動力系統(tǒng)的孤立子波解
    Cr12MoV冷作模具鋼滲鉻層界面能譜分析
    Rotenberg模型中一類遷移算子的譜分析
    沉香GC-MS指紋圖譜分析
    中成藥(2016年8期)2016-05-17 06:08:26
    地震反演子波選擇策略研究
    基于時頻分析的逆合成孔徑雷達成像技術
    對采樣數(shù)據(jù)序列進行時頻分解法的改進
    雙線性時頻分布交叉項提取及損傷識別應用
    基于倒雙譜的地震子波估計方法
    久久婷婷成人综合色麻豆| 神马国产精品三级电影在线观看 | 亚洲av日韩精品久久久久久密| 舔av片在线| 最近在线观看免费完整版| www.精华液| 又黄又粗又硬又大视频| 五月伊人婷婷丁香| 亚洲性夜色夜夜综合| 99精品久久久久人妻精品| 欧美大码av| 国内精品一区二区在线观看| 女人高潮潮喷娇喘18禁视频| 最近视频中文字幕2019在线8| 欧美日韩福利视频一区二区| 麻豆一二三区av精品| 久久香蕉精品热| 日韩 欧美 亚洲 中文字幕| 俄罗斯特黄特色一大片| 亚洲成a人片在线一区二区| 黄色毛片三级朝国网站| www.自偷自拍.com| 叶爱在线成人免费视频播放| 久久精品国产99精品国产亚洲性色| 日日夜夜操网爽| 亚洲熟女毛片儿| 国产免费av片在线观看野外av| 狠狠狠狠99中文字幕| av片东京热男人的天堂| 亚洲中文日韩欧美视频| 午夜福利在线在线| 日本免费a在线| xxxwww97欧美| 亚洲最大成人中文| 国产高清激情床上av| 国产成人精品久久二区二区91| 国产成人欧美在线观看| 色老头精品视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| cao死你这个sao货| 天堂影院成人在线观看| 亚洲精品在线美女| 国产高清视频在线观看网站| 黄色片一级片一级黄色片| 日日夜夜操网爽| 12—13女人毛片做爰片一| 波多野结衣高清作品| 91国产中文字幕| 精品第一国产精品| 色哟哟哟哟哟哟| 88av欧美| 午夜精品久久久久久毛片777| 久久久精品大字幕| 久久婷婷成人综合色麻豆| 长腿黑丝高跟| 亚洲人成伊人成综合网2020| 久久久久免费精品人妻一区二区| 亚洲精品在线观看二区| 九九热线精品视视频播放| 日本一区二区免费在线视频| 好看av亚洲va欧美ⅴa在| 日本五十路高清| 亚洲狠狠婷婷综合久久图片| 亚洲欧美精品综合久久99| 亚洲最大成人中文| 香蕉av资源在线| 久久久久国内视频| 脱女人内裤的视频| 一区二区三区高清视频在线| 亚洲第一电影网av| 亚洲av片天天在线观看| 18禁黄网站禁片免费观看直播| 少妇裸体淫交视频免费看高清 | 不卡一级毛片| 亚洲精品美女久久久久99蜜臀| 欧美成狂野欧美在线观看| 亚洲男人的天堂狠狠| 岛国在线观看网站| 亚洲男人天堂网一区| 91国产中文字幕| 怎么达到女性高潮| 国产av在哪里看| 亚洲九九香蕉| 蜜桃久久精品国产亚洲av| 国产乱人伦免费视频| 亚洲美女视频黄频| 人成视频在线观看免费观看| 欧美大码av| 日本黄色视频三级网站网址| 视频区欧美日本亚洲| 99热6这里只有精品| 久久精品国产综合久久久| 男女午夜视频在线观看| 亚洲国产欧美一区二区综合| 91老司机精品| 亚洲成人久久性| 久久久久久国产a免费观看| 99久久国产精品久久久| 又黄又爽又免费观看的视频| 国产私拍福利视频在线观看| 欧美色视频一区免费| 日日摸夜夜添夜夜添小说| 日韩精品青青久久久久久| 嫩草影院精品99| 欧美中文日本在线观看视频| 国产黄片美女视频| 色综合亚洲欧美另类图片| 日本一二三区视频观看| 亚洲午夜理论影院| 日本 av在线| 1024视频免费在线观看| 国产精品一区二区免费欧美| 婷婷精品国产亚洲av在线| 久9热在线精品视频| 亚洲第一电影网av| 欧美 亚洲 国产 日韩一| 999久久久精品免费观看国产| 国产精品久久久久久亚洲av鲁大| 99re在线观看精品视频| 亚洲18禁久久av| 欧美成人性av电影在线观看| 麻豆成人午夜福利视频| 欧美极品一区二区三区四区| 亚洲人成网站高清观看| 熟妇人妻久久中文字幕3abv| 老司机靠b影院| 美女午夜性视频免费| 岛国视频午夜一区免费看| 嫩草影院精品99| 99在线视频只有这里精品首页| 可以在线观看毛片的网站| 99国产极品粉嫩在线观看| 18禁美女被吸乳视频| 亚洲第一电影网av| 999精品在线视频| 亚洲国产欧美人成| 欧美中文综合在线视频| 999久久久国产精品视频| 制服丝袜大香蕉在线| 欧美中文综合在线视频| 久久精品91蜜桃| 久久久久久国产a免费观看| 性欧美人与动物交配| 免费无遮挡裸体视频| 亚洲男人天堂网一区| 无人区码免费观看不卡| 一二三四在线观看免费中文在| 中文字幕人妻丝袜一区二区| 非洲黑人性xxxx精品又粗又长| 亚洲va日本ⅴa欧美va伊人久久| 国产精品av视频在线免费观看| 两性夫妻黄色片| 日韩欧美在线乱码| 99久久精品热视频| 无人区码免费观看不卡| 国产亚洲av高清不卡| 床上黄色一级片| 老司机午夜十八禁免费视频| 夜夜看夜夜爽夜夜摸| 国产黄a三级三级三级人| 欧美3d第一页| 精品一区二区三区视频在线观看免费| 国产黄色小视频在线观看| 国产精品 欧美亚洲| 国产亚洲精品久久久久久毛片| 亚洲国产欧美人成| 老司机福利观看| 动漫黄色视频在线观看| 性色av乱码一区二区三区2| 亚洲欧美精品综合久久99| 少妇熟女aⅴ在线视频| 宅男免费午夜| 成年免费大片在线观看| 三级男女做爰猛烈吃奶摸视频| 老汉色av国产亚洲站长工具| 久久草成人影院| 亚洲av五月六月丁香网| 久久久久久大精品| 欧美大码av| 久久欧美精品欧美久久欧美| 久久中文字幕人妻熟女| 日韩欧美国产在线观看| 亚洲国产精品sss在线观看| 高清毛片免费观看视频网站| 亚洲成av人片在线播放无| 国产主播在线观看一区二区| 亚洲免费av在线视频| 亚洲专区国产一区二区| 91成年电影在线观看| 欧美另类亚洲清纯唯美| 日韩精品青青久久久久久| 男女视频在线观看网站免费 | 亚洲性夜色夜夜综合| 高清在线国产一区| 国产成人系列免费观看| 日本一区二区免费在线视频| 白带黄色成豆腐渣| 露出奶头的视频| 国产久久久一区二区三区| 精品一区二区三区视频在线观看免费| 色播亚洲综合网| 国产激情欧美一区二区| 国产野战对白在线观看| 亚洲七黄色美女视频| 欧美久久黑人一区二区| 午夜福利高清视频| 制服丝袜大香蕉在线| 国产在线精品亚洲第一网站| 国产精品香港三级国产av潘金莲| 国产欧美日韩一区二区精品| 久热爱精品视频在线9| 亚洲无线在线观看| 国内精品久久久久精免费| 99re在线观看精品视频| 亚洲成a人片在线一区二区| 熟女电影av网| 国产成人精品久久二区二区91| 亚洲 欧美一区二区三区| 五月伊人婷婷丁香| 午夜精品久久久久久毛片777| 妹子高潮喷水视频| 日本 欧美在线| 精品高清国产在线一区| 欧美黑人欧美精品刺激| 成人18禁在线播放| 少妇人妻一区二区三区视频| bbb黄色大片| www日本在线高清视频| 一二三四社区在线视频社区8| 精品欧美一区二区三区在线| 神马国产精品三级电影在线观看 | 伊人久久大香线蕉亚洲五| 欧美黑人巨大hd| 中文字幕av在线有码专区| 少妇被粗大的猛进出69影院| 亚洲欧美激情综合另类| 老司机靠b影院| 一个人观看的视频www高清免费观看 | 中文字幕精品亚洲无线码一区| 在线观看66精品国产| 国产精品自产拍在线观看55亚洲| 亚洲 欧美一区二区三区| 亚洲美女视频黄频| av片东京热男人的天堂| 精品一区二区三区视频在线观看免费| 精品久久久久久久末码| 在线十欧美十亚洲十日本专区| 床上黄色一级片| 亚洲第一电影网av| 午夜a级毛片| 日本撒尿小便嘘嘘汇集6| 久久久久久免费高清国产稀缺| 国产伦人伦偷精品视频| 欧美日本视频| 亚洲国产欧美人成| 国内精品久久久久精免费| 日韩精品免费视频一区二区三区| 两人在一起打扑克的视频| 在线观看www视频免费| 18禁美女被吸乳视频| 亚洲av中文字字幕乱码综合| 一二三四社区在线视频社区8| 欧美中文综合在线视频| 中文字幕人成人乱码亚洲影| 夜夜夜夜夜久久久久| 动漫黄色视频在线观看| 高清毛片免费观看视频网站| 亚洲自偷自拍图片 自拍| 一进一出好大好爽视频| 两性夫妻黄色片| 99热这里只有是精品50| 国产亚洲av高清不卡| 岛国视频午夜一区免费看| 两性夫妻黄色片| 国产精品电影一区二区三区| 少妇粗大呻吟视频| 亚洲狠狠婷婷综合久久图片| 亚洲第一电影网av| 国产99白浆流出| 亚洲av成人av| 久久久精品国产亚洲av高清涩受| 亚洲午夜理论影院| 99re在线观看精品视频| 高清毛片免费观看视频网站| 欧美一级a爱片免费观看看 | 亚洲欧美日韩高清专用| 亚洲国产高清在线一区二区三| cao死你这个sao货| 岛国在线观看网站| 国产高清有码在线观看视频 | 国产黄a三级三级三级人| 正在播放国产对白刺激| 久久婷婷人人爽人人干人人爱| 夜夜夜夜夜久久久久| 国产爱豆传媒在线观看 | 欧美+亚洲+日韩+国产| 嫁个100分男人电影在线观看| 日韩精品中文字幕看吧| 亚洲欧美一区二区三区黑人| 欧美日韩国产亚洲二区| 亚洲人成77777在线视频| 精品国产乱子伦一区二区三区| 国产高清激情床上av| 欧美中文日本在线观看视频| 日韩欧美精品v在线| bbb黄色大片| 三级国产精品欧美在线观看 | 制服人妻中文乱码| 亚洲激情在线av| 亚洲avbb在线观看| 一a级毛片在线观看| 免费在线观看亚洲国产| 神马国产精品三级电影在线观看 | 五月伊人婷婷丁香| 欧美成人免费av一区二区三区| 最好的美女福利视频网| 国产视频内射| 免费看a级黄色片| 18禁黄网站禁片午夜丰满| 不卡一级毛片| 欧美人与性动交α欧美精品济南到| 亚洲色图av天堂| 听说在线观看完整版免费高清| 可以在线观看毛片的网站| 中文字幕人妻丝袜一区二区| 精品国产乱子伦一区二区三区| 久久中文看片网| 伦理电影免费视频| 亚洲成人国产一区在线观看| 精品第一国产精品| 中文字幕人妻丝袜一区二区| 看免费av毛片| 哪里可以看免费的av片| 日日夜夜操网爽| www.999成人在线观看| 制服人妻中文乱码| 深夜精品福利| 精品国产乱子伦一区二区三区| 丰满人妻一区二区三区视频av | 后天国语完整版免费观看| 观看免费一级毛片| 黄色 视频免费看| 欧美性猛交╳xxx乱大交人| 精品久久久久久久久久久久久| 国产精品综合久久久久久久免费| 欧美性猛交╳xxx乱大交人| svipshipincom国产片| 国产亚洲精品久久久久久毛片| 琪琪午夜伦伦电影理论片6080| 99热这里只有精品一区 | 国产精品亚洲美女久久久| 88av欧美| 91麻豆av在线| 久久久久久免费高清国产稀缺| 制服丝袜大香蕉在线| 久久久久久久久免费视频了| 国产亚洲精品第一综合不卡| 久久久久九九精品影院| netflix在线观看网站| 日韩欧美一区二区三区在线观看| 亚洲人与动物交配视频| 亚洲一区中文字幕在线| 久久精品影院6| 久久久国产成人精品二区| 伊人久久大香线蕉亚洲五| 国产成人av教育| 男人舔女人的私密视频| 国产视频内射| 日韩欧美精品v在线| 又爽又黄无遮挡网站| 97人妻精品一区二区三区麻豆| 性欧美人与动物交配| 欧美av亚洲av综合av国产av| 又爽又黄无遮挡网站| 国产成人系列免费观看| 蜜桃久久精品国产亚洲av| 国产精品av视频在线免费观看| 亚洲第一欧美日韩一区二区三区| 一进一出抽搐gif免费好疼| 色播亚洲综合网| 一级毛片精品| 国产一区二区在线观看日韩 | aaaaa片日本免费| 两性午夜刺激爽爽歪歪视频在线观看 | 琪琪午夜伦伦电影理论片6080| 国产精品永久免费网站| 亚洲精品在线观看二区| 18禁裸乳无遮挡免费网站照片| 久久精品国产综合久久久| 美女黄网站色视频| 国产69精品久久久久777片 | 精品人妻1区二区| 午夜免费激情av| 欧美绝顶高潮抽搐喷水| 久久精品影院6| 国产激情久久老熟女| 欧美又色又爽又黄视频| 久久香蕉激情| 又粗又爽又猛毛片免费看| 国内精品久久久久久久电影| 午夜影院日韩av| 欧美一级毛片孕妇| 免费av毛片视频| 国产精品一区二区免费欧美| 国产一区二区激情短视频| 亚洲av中文字字幕乱码综合| 91av网站免费观看| 国产精品98久久久久久宅男小说| 午夜免费激情av| 我的老师免费观看完整版| 成人av在线播放网站| 此物有八面人人有两片| 久久久久国产一级毛片高清牌| 日韩 欧美 亚洲 中文字幕| 国产精品 欧美亚洲| 国产私拍福利视频在线观看| 成人欧美大片| 免费观看人在逋| 757午夜福利合集在线观看| 亚洲一码二码三码区别大吗| 久久香蕉精品热| 麻豆国产av国片精品| 国产爱豆传媒在线观看 | 两个人免费观看高清视频| 国产aⅴ精品一区二区三区波| 亚洲精品色激情综合| 窝窝影院91人妻| 国产精品爽爽va在线观看网站| 9191精品国产免费久久| 亚洲18禁久久av| 国产aⅴ精品一区二区三区波| 亚洲av成人精品一区久久| 男人舔女人的私密视频| 色精品久久人妻99蜜桃| 一级毛片高清免费大全| 在线永久观看黄色视频| 成人国语在线视频| 日韩精品青青久久久久久| 一本一本综合久久| 在线a可以看的网站| 两个人视频免费观看高清| 中文资源天堂在线| 国产精品九九99| 久久人妻av系列| 男女做爰动态图高潮gif福利片| 久久久久国产一级毛片高清牌| 日韩欧美国产一区二区入口| 欧美人与性动交α欧美精品济南到| 中文亚洲av片在线观看爽| 夜夜爽天天搞| 欧美 亚洲 国产 日韩一| 天天添夜夜摸| 中文字幕高清在线视频| 黄片小视频在线播放| 免费观看精品视频网站| 亚洲性夜色夜夜综合| 欧美不卡视频在线免费观看 | 国产精品1区2区在线观看.| 欧美乱妇无乱码| 国产精品久久久久久亚洲av鲁大| 亚洲av电影不卡..在线观看| 国产黄色小视频在线观看| 欧美精品啪啪一区二区三区| 午夜免费激情av| 大型黄色视频在线免费观看| 成年女人毛片免费观看观看9| 久久99热这里只有精品18| 亚洲在线自拍视频| 亚洲九九香蕉| 黑人巨大精品欧美一区二区mp4| av福利片在线| 亚洲国产精品999在线| 国产片内射在线| 日本精品一区二区三区蜜桃| 久久久久免费精品人妻一区二区| 1024视频免费在线观看| 男人舔女人的私密视频| 一卡2卡三卡四卡精品乱码亚洲| 国产精品免费一区二区三区在线| 真人做人爱边吃奶动态| 国产高清有码在线观看视频 | 国产精品久久久久久人妻精品电影| 久久香蕉精品热| 午夜老司机福利片| 久久这里只有精品中国| 国产激情欧美一区二区| 欧美成人性av电影在线观看| 999精品在线视频| 欧美丝袜亚洲另类 | 十八禁人妻一区二区| 一区二区三区高清视频在线| 婷婷丁香在线五月| 亚洲精品美女久久久久99蜜臀| 女人爽到高潮嗷嗷叫在线视频| 99久久无色码亚洲精品果冻| 我的老师免费观看完整版| 久久久国产成人精品二区| 1024香蕉在线观看| 国产亚洲精品av在线| 国产精品久久视频播放| 此物有八面人人有两片| 欧美国产日韩亚洲一区| 亚洲国产欧美人成| 老司机在亚洲福利影院| 一级作爱视频免费观看| 亚洲国产精品sss在线观看| 久久国产精品人妻蜜桃| 国产精品自产拍在线观看55亚洲| 又紧又爽又黄一区二区| 国内精品久久久久久久电影| 一进一出好大好爽视频| 亚洲午夜理论影院| 又大又爽又粗| 午夜福利在线在线| 天堂影院成人在线观看| 听说在线观看完整版免费高清| 在线a可以看的网站| 中文字幕精品亚洲无线码一区| 亚洲美女视频黄频| 看黄色毛片网站| 久久亚洲真实| 好男人电影高清在线观看| 欧美极品一区二区三区四区| 亚洲精品粉嫩美女一区| 欧美激情久久久久久爽电影| 99精品久久久久人妻精品| 午夜福利欧美成人| 啦啦啦韩国在线观看视频| 搡老妇女老女人老熟妇| 免费在线观看影片大全网站| 校园春色视频在线观看| 国产精品九九99| 午夜激情av网站| 三级毛片av免费| 男女做爰动态图高潮gif福利片| 一级毛片高清免费大全| 麻豆av在线久日| 婷婷精品国产亚洲av在线| 国产高清videossex| 久久久久久大精品| 国产亚洲av高清不卡| 两人在一起打扑克的视频| 久久中文看片网| 又黄又粗又硬又大视频| 最近在线观看免费完整版| 亚洲国产看品久久| 可以在线观看毛片的网站| 国产精品久久久久久人妻精品电影| 亚洲专区中文字幕在线| 国产片内射在线| 久久中文字幕人妻熟女| 国产99白浆流出| 老司机福利观看| av有码第一页| 欧美精品亚洲一区二区| 99久久国产精品久久久| 香蕉av资源在线| 国产伦一二天堂av在线观看| 毛片女人毛片| 午夜福利欧美成人| 看黄色毛片网站| 久久久久九九精品影院| 国产成人系列免费观看| 窝窝影院91人妻| 精品久久久久久久末码| 国产亚洲精品第一综合不卡| 夜夜躁狠狠躁天天躁| 大型av网站在线播放| 欧美+亚洲+日韩+国产| 国产人伦9x9x在线观看| 国产亚洲精品久久久久5区| 老司机靠b影院| 欧美中文综合在线视频| 最近视频中文字幕2019在线8| 精品国产亚洲在线| 成人亚洲精品av一区二区| 亚洲无线在线观看| 亚洲国产高清在线一区二区三| 国产v大片淫在线免费观看| 变态另类丝袜制服| 又紧又爽又黄一区二区| 久久精品国产99精品国产亚洲性色| 亚洲熟妇熟女久久| 亚洲专区国产一区二区| 国产精品亚洲av一区麻豆| 波多野结衣高清作品| 国产黄片美女视频| 18禁裸乳无遮挡免费网站照片| 午夜免费成人在线视频| 欧美一级a爱片免费观看看 | 日韩欧美 国产精品| 国产久久久一区二区三区| 国产精品av久久久久免费| 亚洲精品中文字幕在线视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品av久久久久免费| 两个人的视频大全免费| 久久伊人香网站| 在线观看美女被高潮喷水网站 | 午夜免费激情av| 国产一区二区在线观看日韩 | 久久久国产成人精品二区| 变态另类丝袜制服| 精品不卡国产一区二区三区| 久久久久国产一级毛片高清牌| 琪琪午夜伦伦电影理论片6080| 国产精品久久久久久久电影 | 国产69精品久久久久777片 | x7x7x7水蜜桃| 国产欧美日韩一区二区三| 欧美日韩国产亚洲二区| av福利片在线| 久久中文看片网| 欧美不卡视频在线免费观看 | 1024香蕉在线观看|