黏性介質(zhì)薄砂層時(shí)頻特征響應(yīng)分析

范宇婷

( 1. 中國(guó)科學(xué)院 地質(zhì)與地球物理研究所，北京 100029； 2. 中國(guó)海洋石油總公司 信息化部，北京 100010 )

?

黏性介質(zhì)薄砂層時(shí)頻特征響應(yīng)分析

范宇婷1,2

( 1. 中國(guó)科學(xué)院 地質(zhì)與地球物理研究所，北京 100029； 2. 中國(guó)海洋石油總公司 信息化部，北京 100010 )

地下介質(zhì)并不是理想的彈性介質(zhì)，地層存在黏性吸收作用，對(duì)地震子波的頻帶及峰值頻率等參數(shù)產(chǎn)生影響.基于波場(chǎng)延拓理論，采用波動(dòng)方程正演模擬方法，對(duì)3種不同黏性吸收程度及厚度的薄砂層進(jìn)行正演模擬；采用廣義S變換，分析反射復(fù)合波的瞬時(shí)頻譜.結(jié)果表明：薄砂層瞬時(shí)頻譜的峰值頻率與無(wú)黏性吸收作用下相比大幅度變化，導(dǎo)致由峰值頻率預(yù)測(cè)薄砂層厚度出現(xiàn)較大誤差.陷頻頻率受地層黏性吸收作用比較小，相對(duì)比較穩(wěn)定.利用陷頻頻率出現(xiàn)的周期性與薄砂層厚度之間的關(guān)系預(yù)測(cè)薄砂層厚度，可以為復(fù)雜地質(zhì)環(huán)境下有利薄儲(chǔ)層的定性識(shí)別和厚度定量預(yù)測(cè)提供支持.

黏性吸收； 薄砂層； 正演模擬； 瞬時(shí)頻譜； 廣義S變換

0 引言

地下介質(zhì)既有彈性又存在黏滯性[1]，地震波在傳播過程中存在黏性吸收作用，導(dǎo)致地震波高頻能量衰減，地震波的峰值頻率降低，譜寬變窄[2-3].Yuan Chunfang、Richer N H等以Kelvin-Voigt模型[4]和Maxwell模型[5]為基礎(chǔ)，提出基本黏彈性介質(zhì)理論[6-7].楊仁虎等指出地震波在地下介質(zhì)傳播過程中，存在高頻能量的吸收與衰減[8].黏性吸收作用是地層固有存在的特性，因此在地震波場(chǎng)模擬及偏移成像處理中，必須考慮地層的黏性吸收作用，才能更好地研究地震波在地下介質(zhì)中各類傳播現(xiàn)象[9].王小杰等研究地震波在黏彈性介質(zhì)中的傳播與反射特征[10]，得出在黏彈性介質(zhì)中需對(duì)地層加入品質(zhì)因子，再計(jì)算反射系數(shù)，才更接近真實(shí)地下介質(zhì).

目前，人們[11-12]應(yīng)用褶積模型，考慮震源子波頻譜特征和反射系數(shù)譜變化規(guī)律兩種方式，研究薄砂層地震響應(yīng)特性，分析震源子波類型、頻率、相位和譜寬等特性對(duì)薄砂層地震響應(yīng)特征的影響，以及不同厚度薄砂層反射系數(shù)譜的變化規(guī)律，給出峰值頻率或陷頻頻率與薄砂層厚度之間的定量關(guān)系[13-18].薄砂層正演模擬方法由時(shí)域褶積方法向波動(dòng)方程方面轉(zhuǎn)變[19]，使模擬后的波場(chǎng)信息更加豐富、更加接近采集實(shí)際.薄砂層分析由傅里葉變換向廣義S變換[20]和匹配追蹤[21-23]等先進(jìn)時(shí)頻分析方法方向轉(zhuǎn)變，注重反射波頻譜的瞬時(shí)特性，預(yù)測(cè)方法也由定性轉(zhuǎn)向定量[24].

薄砂層的正演模擬和解釋識(shí)別是建立在地層為彈性介質(zhì)假設(shè)基礎(chǔ)上的，忽略地層的黏性吸收作用，現(xiàn)有的方法無(wú)法適用于實(shí)際的地震資料解釋，主要存在問題：首先，黏性吸收使反射波的振幅隨著傳播時(shí)間的增加而逐漸降低，使得利用振幅定量預(yù)測(cè)薄砂層厚度的一類方法出現(xiàn)較大誤差.其次，黏性吸收作用的頻散作用使地震子波產(chǎn)生時(shí)變扭曲，導(dǎo)致地震波的走時(shí)出現(xiàn)誤差.由于地震子波高頻成分存在損失，導(dǎo)致地震子波的波形變胖，造成薄砂層的視厚度的定性判斷出現(xiàn)較大誤差.再次，黏性吸收作用導(dǎo)致地震子波的頻帶寬度變窄，峰值頻率向低頻移動(dòng)，使得以頻率方法為主的薄砂層厚度定量預(yù)測(cè)出現(xiàn)較大誤差.因此，有必要研究在黏性吸收作用下薄砂層時(shí)頻響應(yīng)特征變化規(guī)律，尋找對(duì)黏性吸收參數(shù)變化不敏感的參數(shù)，以更好地指導(dǎo)薄砂層的定性解釋和厚度的定量預(yù)測(cè).

根據(jù)波場(chǎng)延拓理論，采用波動(dòng)方程正演模擬方法，利用3種不同黏性吸收程度的黏性介質(zhì)，對(duì)不同厚度的薄砂層進(jìn)行波動(dòng)方程正演模擬，從獲取炮記錄中抽取零偏移距地震道；利用時(shí)頻分析方法得到反射復(fù)合波的瞬時(shí)頻譜，考察不同厚度、不同黏性吸收程度薄砂層的時(shí)頻特征變化規(guī)律；從頻域振幅譜、峰值頻率、陷頻頻率，以及時(shí)域波形、振幅、視厚度中優(yōu)選對(duì)黏性吸收參數(shù)敏感的屬性特征參數(shù)，為黏性吸收介質(zhì)薄儲(chǔ)層的定性識(shí)別和厚度的定量預(yù)測(cè)提供指導(dǎo).

1 正演模擬方法

采用波動(dòng)方程相移正演方法對(duì)薄砂層進(jìn)行正演模擬.設(shè)震源子波的頻譜為f(ω)，在黏性吸收介質(zhì)條件下，各水平界面處的地震波正傳波場(chǎng)為

P′(kx,zi,ω)=f(ω)e-jkzizi,

(1)

式中：P′(kx,zi,ω)為地震波的正傳波場(chǎng)；j為虛數(shù)單位；zi為第i個(gè)反射界面所在的深度；kzi為垂直方向上的波數(shù)，可表示為

(2)

式中：ci為第i層的地震波傳播速度；Qi為第i層的品質(zhì)因子，表征地層的黏性吸收特性的強(qiáng)弱；kx為水平波數(shù).

當(dāng)?shù)卣鸩▊鞑サ阶钌罱缑嫣帟r(shí)，經(jīng)過界面反射后從z=zi傳播到z=zi-1處的反傳波場(chǎng)，與z=zi-1處產(chǎn)生的反射波場(chǎng)迭加得到z=zi-1處的總反射波場(chǎng)：

P″(kx,zi-1,ω)=P″(kx,zi,ω)e-j(zi-zi-1)kzi+P′(kx,zi-1,ω)Rzi-1.

(3)

將迭加后的波場(chǎng)P″(kx,zi-1,ω)作為下一界面處延拓的初始值，依次迭代逐步向上延拓，直到延拓到地面被各接收點(diǎn)所接收，形成地面反射波場(chǎng)，記為P″(kx,0,ω)，將地面波場(chǎng)作反傅里葉變換：

(4)

式中：P(x,t)為薄砂層的地震正演數(shù)據(jù).

2 地質(zhì)模型及數(shù)值模擬

建立黏性吸收介質(zhì)薄砂層地質(zhì)模型(見圖1).該地質(zhì)模型由上下均勻背景泥巖和中間薄砂層組成，其中薄砂層頂界的反射因數(shù)設(shè)為-0.5，薄砂層底界的反射因數(shù)設(shè)為0.5.整個(gè)介質(zhì)模型的速度為2 000 m/s；設(shè)定上下背景泥巖的厚度，使薄砂層的中心位于2 ms位置；為了研究方便，無(wú)論是背景泥巖還是薄砂層砂巖，設(shè)定相同的品質(zhì)因子.按照黏性吸收程度的不同，分別考慮強(qiáng)吸收(Q=66)、中等吸收(Q=100)和弱吸收(Q=200)條件，考察不同黏性吸收程度對(duì)薄砂層時(shí)頻特性的影響；地震子波采用峰值頻率為30 Hz零相位雷克子波；采用波動(dòng)方程正演模擬獲取炮集數(shù)據(jù)；采取中間放炮、兩邊接收的采集方式，檢波器之間距離設(shè)為10 m；每個(gè)炮集共包含1 000道數(shù)據(jù)；地震道采樣點(diǎn)個(gè)數(shù)為2 000，采樣間隔為2 ms.薄砂層厚度以2 m為間隔從1 m變化到29 m，3種黏性吸收條件共模擬45套地震數(shù)據(jù)體.

圖1 黏性吸收薄砂層地層模型Fig.1 The thin layer model in viscous absorption media

從炮集中抽取零偏移距地震道，分別讀取時(shí)域最大振幅和薄砂層的視厚度(當(dāng)頂?shù)追瓷湟驍?shù)相反時(shí)，分別讀取薄砂層界面處反射子波的峰、谷對(duì)應(yīng)的時(shí)刻).采用廣義S變換對(duì)該地震道進(jìn)行時(shí)頻變換，在時(shí)頻分析圖上薄砂層對(duì)應(yīng)2 ms位置處，獲取薄砂層反射波的瞬時(shí)頻譜，并計(jì)算瞬時(shí)振幅譜；從瞬時(shí)振幅譜分別讀出峰值頻率、陷頻頻率及峰值頻率的振幅幅度，統(tǒng)計(jì)分析45套時(shí)頻分析數(shù)據(jù)；考察不同黏性吸收條件下不同厚度薄砂層的時(shí)頻域特征的變化規(guī)律.

3 數(shù)據(jù)分析

3.1 峰值頻率

圖2 峰值頻率隨薄砂層厚度變化關(guān)系Fig.2 The crossplot for peak frequency versus thickness of thin layer

在不同黏性吸收條件下，薄砂層峰值頻率總體變化規(guī)律：在Q不變條件下，隨著薄砂層厚度的減薄(大于5 m)，峰值頻率逐漸增大(見圖2)；當(dāng)薄砂層厚度小于5 m時(shí)，峰值頻率隨著薄砂層厚度的減薄而逐漸減小；總體變化規(guī)律與無(wú)黏性吸收地層的類似.當(dāng)Q逐漸減小時(shí)，峰值頻率逐漸降低，并且地層越薄，降低的幅度越大.這說(shuō)明在黏性吸收條件下，瞬時(shí)頻譜的峰值頻率向低頻方向移動(dòng)；在無(wú)黏性吸收條件下，當(dāng)薄砂層厚度等于1/2波長(zhǎng)時(shí)，峰值頻率發(fā)生躍變，由低頻峰轉(zhuǎn)變?yōu)楦哳l峰.在黏性吸收條件下，高頻能量被相對(duì)大幅度衰減，使得瞬時(shí)頻譜中并不存在高頻峰，峰值頻率并不發(fā)生躍變，隨著薄砂層厚度的增加而逐漸降低.在給出Q條件下，地層厚度越大時(shí)，不同Q之間與無(wú)黏性地層相比，最小頻率差為6 Hz(地層厚度為25 m)，最大峰值頻率差為18 Hz(地層厚度為1 m)，說(shuō)明在薄砂層(小于1/4波長(zhǎng))條件下黏性吸收大幅改變峰值頻率，從而使由峰值頻率預(yù)測(cè)薄砂層厚度出現(xiàn)較大的誤差.

3.2 瞬時(shí)頻譜最大幅度

在不同黏性吸收條件下，峰值頻率幅度的總體變化規(guī)律：隨著薄砂層厚度的減薄，瞬時(shí)頻譜的峰值頻率幅度先增大后減小，總體趨勢(shì)與無(wú)黏性吸收地層的一致(見圖3).在黏性吸收條件下，并不是像無(wú)黏性吸收地層一樣在1/4波長(zhǎng)處取得最大值，在1/2波長(zhǎng)處取得極小值，而是極大值點(diǎn)向薄砂層增厚方向移動(dòng).峰值頻率幅度隨著Q減小而逐漸降低，其極大值所對(duì)應(yīng)薄砂層厚度逐漸增大(當(dāng)Q=200時(shí)，為15 m；當(dāng)Q=100時(shí)，為17 m；當(dāng)Q=66時(shí)，為21 m；當(dāng)無(wú)黏性吸收時(shí)，為13 m，在1/4波長(zhǎng)處)，到達(dá)幅度極大值后，幅度隨薄砂層厚度變化逐漸趨于平緩.

3.3 陷頻頻率

在不同黏性吸收條件下，陷頻頻率隨薄砂層厚度變化的總體規(guī)律：在存在與不存在黏性吸收條件下薄砂層陷頻頻率近似相等，總體上隨著薄砂層厚度減薄，陷頻頻率有逐漸增大的趨勢(shì)(見圖4).一般薄砂層瞬時(shí)頻譜只存在單一陷頻點(diǎn)，但是在不同Q的瞬時(shí)頻譜中可以看到只有少數(shù)薄砂層厚度存在陷頻點(diǎn).當(dāng)Q為200時(shí)，19 m以上地層存在陷頻點(diǎn)；當(dāng)Q為100時(shí)，21 m以上地層存在陷頻點(diǎn)；當(dāng)Q為66時(shí)，23 m以上地層存在陷頻點(diǎn).在無(wú)黏性吸收條件下，15 m以上地層存在陷頻點(diǎn).隨著Q的減小，具有陷頻點(diǎn)的薄砂層數(shù)量減少，主要原因是位于陷頻點(diǎn)右側(cè)的高頻峰易于受到強(qiáng)烈的高頻衰減，當(dāng)衰減到一定程度時(shí)，高頻峰的能量全部損失殆盡，在瞬時(shí)頻譜上表現(xiàn)為無(wú)明顯陷頻頻率點(diǎn).

薄砂層頻譜屬于周期性頻譜，陷頻頻率位置呈周期性分布，并且相鄰兩個(gè)陷頻頻率差是薄層時(shí)間厚度的倒數(shù)[25]，即

Pf=fn-fn-1=1/τ,

(5)

式中：Pf為陷頻周期；fn為第n個(gè)陷頻頻率；τ為薄層的時(shí)間厚度.

由圖4可知，地層的黏性吸收作用對(duì)薄砂層的陷頻頻率影響較小，且當(dāng)薄砂層厚度大于15 m時(shí)，陷頻頻率與地層厚度呈反比關(guān)系，因此，當(dāng)?shù)貙哟嬖陴ば晕諘r(shí)，可以利用陷頻頻率預(yù)測(cè)地層厚度，其精確程度好于峰值頻率的.

3.4 瞬時(shí)譜寬

定義振幅譜最大幅度的10%對(duì)應(yīng)的頻率差為薄砂層局部頻譜的瞬時(shí)譜寬.當(dāng)?shù)貙訜o(wú)黏性吸收(Q=∞)時(shí)，瞬時(shí)譜寬的變化規(guī)律：當(dāng)薄砂層厚度大于1/2波長(zhǎng)時(shí)，瞬時(shí)譜寬隨著薄砂層厚度減薄逐漸增大，瞬時(shí)譜寬為高頻峰的瞬時(shí)譜寬；當(dāng)薄砂層厚度不大于1/2波長(zhǎng)時(shí)，隨著薄砂層厚度的減薄，瞬時(shí)譜寬逐漸增寬；在地層厚度為1 m時(shí)略有減小，瞬時(shí)譜寬為低頻峰的譜寬.

圖3 瞬時(shí)頻譜最大幅度隨薄砂層厚度變化關(guān)系Fig.3 The crossplot for maximum amplitude of frequency spectrum versus thickness of thin layer

圖4 瞬時(shí)頻譜陷頻頻率隨薄砂層厚度變化關(guān)系Fig.4 The notched frequency of instantaneous amplitude spectrum versus thickness of thin layer

圖5 瞬時(shí)譜寬隨薄砂層厚度變化關(guān)系Fig.5 The bandwidth of spectrum vary with single layer thickness

當(dāng)?shù)貙哟嬖陴ば晕諘r(shí)，由于高頻峰具有強(qiáng)烈衰減作用，瞬時(shí)頻譜只存在低頻峰，因此瞬時(shí)譜寬基本上為低頻峰的譜寬(見圖5)，并且瞬時(shí)譜寬隨著薄砂層厚度的減薄而逐漸增寬.高Q值增寬的幅度較迅速，低Q值增寬的幅度較緩.當(dāng)?shù)貙雍穸刃∮?5 m時(shí)，相同地層厚度的譜寬在Q為無(wú)窮大時(shí)最寬，在Q為66時(shí)譜寬最小；在地層厚度為15～23 m條件下，在Q為200時(shí)，頻譜最寬，無(wú)黏性吸收時(shí)譜寬最小；其他Q值介于兩者之間.當(dāng)?shù)貙雍穸却笥?3 m時(shí)，在Q為66時(shí)，瞬時(shí)譜寬最寬；在Q為200時(shí)，瞬時(shí)譜寬最小.

3.5 時(shí)域最大振幅及視厚度

不同Q對(duì)地震波的走時(shí)產(chǎn)生一定的影響，進(jìn)而影響由峰谷走時(shí)計(jì)算的視厚度.以1/4波長(zhǎng)為界，當(dāng)薄砂層厚度小于1/4波長(zhǎng)時(shí)，由峰谷對(duì)應(yīng)走時(shí)計(jì)算的視厚度大于薄砂層的真實(shí)厚度，且Q越小，偏差越大；地層厚度為15～19 m時(shí)，由Q為66計(jì)算的走時(shí)大于理論值，其他的小于理論值；當(dāng)?shù)貙雍穸却笥?9 m時(shí)，所有由Q計(jì)算的視厚度小于理論值.由Q為66與200計(jì)算的視厚度更接近于理論值(見圖6).

反射復(fù)合波時(shí)域最大振幅隨薄砂層厚度變化的基本趨勢(shì)是一致的(見圖7)，即最大振幅先隨薄砂層厚度的減薄，振幅先增大；在達(dá)到最大值后，隨薄砂層厚度的減小而減小.時(shí)域最大振幅受Q影響較大，Q越小，時(shí)域振幅越小，且最大值并不在1/4波長(zhǎng)處取得，而是向薄砂層厚度增大方向偏移.

圖6 薄砂層的視厚度隨實(shí)際厚度變化關(guān)系Fig.6 The comparison between apparent thickness and the true thickness of thin layer

圖7 時(shí)域最大振幅隨薄砂層厚度變化關(guān)系Fig.7 The maximum amplitude in time-domain versus thickness of thin layer

4 討論

譜寬與瞬時(shí)頻譜的峰值頻率幅度密切相關(guān).當(dāng)薄砂層厚度小于1/4波長(zhǎng)時(shí)，薄砂層具有升頻降幅的作用[26]，瞬時(shí)頻譜的峰值幅度是逐漸降低的，且頻譜的最大幅度變化率是比較平緩的.因此，瞬時(shí)譜寬是測(cè)量瞬時(shí)頻譜的基部，幅度越小，測(cè)量的譜寬越接近基部，譜寬越大.當(dāng)?shù)貙哟嬖陴ば晕諘r(shí)，隨黏性吸收作用的增強(qiáng)，高頻能量全部損失殆盡，因此頻譜寬度逐漸減小.當(dāng)薄砂層厚度大于1/4波長(zhǎng)時(shí)，瞬時(shí)頻譜的幅度隨著地層厚度的增加先減小后趨于穩(wěn)定，因此頻譜寬度變化不大.當(dāng)薄砂層厚度大于1/2波長(zhǎng)時(shí)，由于陷頻頻率存在，導(dǎo)致一部分高頻能量得以增強(qiáng)，使得譜寬隨著Q的減小進(jìn)一步增強(qiáng).

根據(jù)正演模擬結(jié)果，薄砂層的頻率特征參數(shù)受地層黏性吸收作用影響較大，基于頻率參數(shù)的薄砂層厚度預(yù)測(cè)方法將受到不同程度的影響.因此，在進(jìn)行薄砂層厚度預(yù)測(cè)前，應(yīng)對(duì)采集到的地面地震資料進(jìn)行黏性吸收補(bǔ)償疊前時(shí)間偏移成像處理[27]，才能正確補(bǔ)償?shù)貙拥酿ば晕兆饔茫瑥亩档陀勺兓蟮姆逯殿l率預(yù)測(cè)薄砂層厚度所帶來(lái)的誤差.

此外，在黏性介質(zhì)傳播過程中，地震波將根據(jù)地層吸收特性的強(qiáng)弱，出現(xiàn)不同程度的頻散，使得不同頻率的地震波以不同的速度進(jìn)行傳播，導(dǎo)致地震波的走時(shí)發(fā)生變化，使得由峰谷走時(shí)方法預(yù)測(cè)薄砂層視厚度出現(xiàn)較大的偏差.

5 結(jié)論

(1)將品質(zhì)因子引入到波動(dòng)方程正演模擬中，在波場(chǎng)延拓中考慮地層的黏性吸收特性的影響，使正演模擬得到的薄砂層時(shí)頻響應(yīng)規(guī)律更加接近地下實(shí)際.

(2)地層存在黏性吸收作用，薄砂層瞬時(shí)頻譜的峰值頻率與無(wú)黏性吸收作用下相比大幅變化，導(dǎo)致由峰值頻率預(yù)測(cè)薄砂層厚度出現(xiàn)較大誤差.在進(jìn)行薄砂層厚度預(yù)測(cè)前，必須進(jìn)行地震資料的高分辨疊前成像處理，正確補(bǔ)償?shù)貙拥酿ば晕兆饔茫拍芙档陀煞逯殿l率預(yù)測(cè)薄砂層厚度的誤差.

(3)由于黏性吸收作用造成頻散，導(dǎo)致地震波的走時(shí)發(fā)生變化.由峰谷走時(shí)預(yù)測(cè)的視厚度受黏性吸收作用影響，使得測(cè)量結(jié)果與理論值出現(xiàn)較大的偏差.

(4)陷頻頻率受地層黏性吸收作用比較小，在黏性吸收條件下，可以利用陷頻頻率出現(xiàn)的周期性與薄砂層厚度的定量關(guān)系預(yù)測(cè)薄砂層厚度.

[1] 杜啟振,楊慧珠.線性黏彈性各向異性介質(zhì)速度頻散和衰減特征研究[J].物理學(xué)報(bào),2002,51(9):2101-2018. Du Qizhen, Yang Huizhu. Velocity dispersion and attenuation in anisotropic linear viscoelastic media [J]. Acta Physica Sinica, 2002,51(9):2101-2018.

[2] 鄭鴻明,薛為平,蔣琳,等.黏彈性介質(zhì)中地震波的吸收衰減補(bǔ)償方法[J].新疆石油地質(zhì),2006,27(5):556-558. Zheng Hongming, Xue Weiping, Jiang Lin, et al. Compensation method of seismic wave absorption and attenuation in viscoelastic media [J]. Xinjiang Petroleum Geology, 2006,27(5):556-558.

[3] 孫成禹.地震波理論與方法[M].東營(yíng):中國(guó)石油大學(xué)出版社,2007. Sun Chengyu. Seismic wave theory and method [M]. Dongying: China University of Petroleum Press, 2007.

[4] Yuan Chunfang. Seismic wave propagating in Kelvin-Voigt homogeneous visco-elastic media [J]. Science in China: Series D Earth Sciences, 2006,49(2):147-153.

[5] Richer N H. Transient waves in visco-elastic media [M]. Amsterdam: Elsevier Scientific Pub. Co., 1977:1-278.

[6] Aki K, Richards P G. Quantitative seismology [M]. 2nd ed. Sansalito: University Science Books, 2002.

[7] 張中杰.地震各向異性研究進(jìn)展[J].地球物理學(xué)進(jìn)展,2002,17(2):281-293. Zhang Zhongjie. A review of the seismic anisotropy and its applications[J]. Progress in Geophysics, 2002,17(2):281-293.

[8] 楊仁虎,常旭,劉伊克,等.黏彈性介質(zhì)地震波場(chǎng)正演數(shù)值模擬[C]∥中國(guó)地球物理.合肥:中國(guó)科學(xué)技術(shù)大學(xué)出版社,2009:428-429. Yang Renhu, Chang Xu, Liu Yike, et al. Viscoelastic wave inhomogeneous isotropic medium is speech based on numerical simulation [C]∥Geophysical Chinese. Hefei: University of Science and Technology of China Press, 2009:428-429.

[9] 范家參.黏彈性介質(zhì)中的地震波[J].地震研究,2001,24(4):358-359. Fan Jiashen. Seismic waves is viscoelastic medium [J]. Journal of Seismological Reseearch, 2001,24(4):358-359.

[10] 王小杰,印興耀,吳國(guó)忱,等.黏彈性介質(zhì)地震波傳播特征及反射特征研究[J].物探化探計(jì)算技術(shù),2012,34(3):258-266. Wang Xiaojie, Yin Xingyao, Wu Guochen, et al. Study on propagation characteristics and reflection characteristics of seismic wave in viscoelastic medium [J]. Computing Techniques for Geophysical and Geochemical Exploration, 2012,34(3):258-266.

[11] 崔鳳林,管葉君.時(shí)頻分析——薄互層結(jié)構(gòu)研究的新途徑[J].石油物探,1992,31(2):1-15. Cui Fenglin, Guan Yejun. Time-Frequency analysis: A new way for tin interbeds examination [J]. Geophysical Prospecting for Petroleum, 1992,31(2):1-15.

[12] Knapp R W. Energy distribution in wavelets and implications on resolving power [J]. Geophysics, 1993,58(1):39-46.

[13] De Voogd N, den Rooijen N. Thin-layer response and spectral bandwidth [J]. Geophysics, 1983,48(1):12-18.

[14] 凌云,西曉宇,孫德勝,等.薄儲(chǔ)層疊后反演影響因素分析與地震屬性解釋研究[J].石油物探,2008,47(6):531-558. Ling Yun, Xi Xiaoyu, Sun Desheng, et al. Analysis on affecting factors of post-stack inversion and seismic attribute interpretation of thin reservoir [J]. Geophysical Prospecting for Petroleum, 2008,47(6):531-558.

[15] 孫魯平,鄭曉東,首皓,等.薄層地震峰值頻率與厚度關(guān)系研究[J].石油地球物理勘探,2010,45(2):254-259. Sun Luping, Zheng Xiaodong, Shou Hao, et al. The studies on relationship between thin layer seismic peak frequency and its thickness [J]. Geophysical Prospecting for Petroleum, 2010,45(2):254-259.

[16] Okaya D A. Spectral properties of the earth's contribution to seismic resolution [J]. Geophysics, 1995,60(1):244-251.

[17] Marfurt K J, Kirlin R L. Narrow-band spectral analysis and thin-bed tuning [J]. Geophysics, 2001,66(4):1274-1283.

[18] Partyka G A, Gridley J M, Lopez J. Interpretational applications of spectral decomposition in reservoir characterization [J]. The Leading Edge, 1999,18(3):353-360.

[19] 周麗,顧漢明,馬靈偉,等.基于波動(dòng)方程正演模擬分析薄砂層含不同流體的AVO特征[J].地質(zhì)科技情報(bào),2013,32(2):169-173. Zhou Li, Gu Hanming, Ma Lingwei, et al. Analysis of AVO characterisitics on thin sand reservoir with various fluid based seismic modeling [J]. Geological Science and Technology Information, 2013,32(2):169-173.

[20] 熊曉軍,賀振華,黃德濟(jì),等.廣義S變換在地震高分辨處理中的應(yīng)用[J].勘探地球物理進(jìn)展,2006,6(29):415-418. Xiong Xiaojun, He Zhenhua, Huang Deji, et al. Application of generalized S transform in seismic high resolution processing [J]. Progress in Exploration Geophysics, 2006,6(29):415-418.

[21] Mallat S, Zhang Z. Matching pursuits with time-frequency dictionaries [J]. IEEE Transactions on signal processing, 1993,41(12):3397-3415.

[22] 楊昊,鄭曉東,李勁松,等.基于匹配追蹤的薄層自動(dòng)解釋方法[J].石油地球物理勘探,2013,48(3):429-435. Yang Hao, Zheng Xiaodong, Li Jinsong, et al. Thin-bed automatic interpretation based on matching pursuit [J]. Geophysical Prospecting for Petroleum, 2013,48(3):429-435.

[23] 宋維琪,朱衛(wèi)星,孫英杰.復(fù)數(shù)子波匹配追蹤算法識(shí)別薄層砂體[J].地球物理學(xué)進(jìn)展,2007,22(6):1796-1801. Song Weiqi, Zhu Weixing, Sun Yingjie. Identify bed layer sandbody by complex wavelet matching algorithm [J]. Progress in Geophysics, 2007,22(6):1796-1801.

[24] 劉震,張萬(wàn)選,張厚福.儲(chǔ)層厚度定量解釋方法研究[J].石油地球物理勘探,1991,26(6):777-784. Liu Zhen, Zhang Wanxuan, Zhang Houfu. A research into quantitative interpretation of reservoir thickness [J]. Geophysical Prospecting for Petroleum, 1991,26(6):777-784.

[25] 黃緒德.薄層陷頻法[J].勘探地球物理進(jìn)展,2002,25(5):1-6. Huang Xude. Discussion on notches in thin-bed [J]. Progress in Exploration Geophysics, 2002,25(5):1-6.

[26] 李雪英,陳樹民,王建民,等.薄層時(shí)頻特征的正演模擬[J].地球物理學(xué)報(bào),2012,55(10):3410-3419. Li Xueying, Chen Shumin, Wang Jianmin, et al. Forward modeling studies on the time-frequency characteristics of thin layers [J]. Chinese Journal of Geophysics, 2012,55(10):3410-3419.

[27] 李雪英,文慧儉,張桐,等.黏性介質(zhì)疊前時(shí)間偏移方法[J].東北石油大學(xué)學(xué)報(bào),2013,37(4):98-104. Li Xueying, Wen Huijian, Zhang Tong, et al. The method of prestack time migration in anelastic media [J]. Journal of Northeast Petroleum University, 2013,37(4):98-104.

2015-04-17；編輯：任志平

國(guó)家自然科學(xué)基金創(chuàng)新研究群體項(xiàng)目(41021063)

范宇婷(1983-)，女，博士研究生，工程師，主要從事固體地球物理方面的研究.

P631

A

2095-4107(2015)03-0001-06

DOI 10.3969/j.issn.2095-4107.2015.03.001