李東慶， 魏建新， 狄?guī)妥專?高 峰， 蘭曉光

( 1. 中國石油大學(xué)(北京) 地球物理與信息工程學(xué)院，北京 102249； 2. 中國石油大學(xué)(北京) 油氣資源與探測國家重點實驗室，北京 102249； 3. 中國石油大學(xué)(北京) CNPC物探重點實驗室，北京 102249 )

?

蘑菇狀火山巖模型的二維地震正演模擬分析

李東慶1,2,3， 魏建新1,2,3， 狄?guī)妥?,2,3， 高 峰1,2,3， 蘭曉光1,2,3

( 1. 中國石油大學(xué)(北京) 地球物理與信息工程學(xué)院，北京 102249； 2. 中國石油大學(xué)(北京) 油氣資源與探測國家重點實驗室，北京 102249； 3. 中國石油大學(xué)(北京) CNPC物探重點實驗室，北京 102249 )

為了研究蘑菇狀火山巖儲層的地震響應(yīng)特征，根據(jù)地質(zhì)資料構(gòu)建蘑菇狀火山巖地質(zhì)模型，分別采用物理模擬和數(shù)值模擬方法對蘑菇狀火山巖模型進(jìn)行正演模擬。結(jié)果表明，蘑菇狀火山巖的地震響應(yīng)特征總體趨勢一致，即火山巖上部外形為蘑菇頭狀，與其結(jié)構(gòu)形態(tài)類似；火山巖下部產(chǎn)生“蘑菇無根”現(xiàn)象，在細(xì)節(jié)上存在差異。物理模擬方法得到的地震波信息更加豐富，且在模擬過程中地震波的傳播方式更符合地下實際情況；數(shù)值模擬方法沒有噪音干擾，資料品質(zhì)好，且可以得到波場快照信息。兩種正演模擬方法具有互補性，為研究復(fù)雜地質(zhì)情況下的火山巖地震響應(yīng)特征提供更好的支持。

蘑菇狀火山巖； 物理模擬； 數(shù)值模擬； 地震響應(yīng)

0 引言

我國東部的松遼盆地和西部的準(zhǔn)噶爾盆地等深層油氣藏勘探成果表明，火山巖具有良好的油氣勘探前景，逐漸成為全球油氣勘探的重要領(lǐng)域，越來越受到油氣勘探界的重視[1-6]。火山巖儲層具有埋藏深、與地層接觸關(guān)系復(fù)雜多變、非均質(zhì)性強及速度高等特點，造成火山巖地震響應(yīng)特征復(fù)雜化，進(jìn)而給火山巖儲層的識別、預(yù)測及地質(zhì)規(guī)律的認(rèn)識帶來困難[7-9]。需要通過綜合地質(zhì)、鉆井、測井和地震等資料[10-12]，運用地震反演等技術(shù)對火山巖的儲層信息進(jìn)行更準(zhǔn)確的描述，而地震反演的基礎(chǔ)是地震正演模擬技術(shù)。采用地震正演模擬技術(shù)，研究火山巖的地震響應(yīng)特征可以減少地震反演的多解性，增加地震解釋的準(zhǔn)確度，從而實現(xiàn)對火山巖內(nèi)部機構(gòu)的準(zhǔn)確劃分與刻畫。因此，地震正演模擬技術(shù)在火山巖地震響應(yīng)及儲層預(yù)測研究中的重要性日益凸顯[13-14]。

王建花等[15]認(rèn)為，由于火山巖往往作為高速異常體侵入地層，因此火山巖的存在可以引起能量屏蔽。毛小平等[16]、李素華等[17]采用波動方程數(shù)值模擬方法，分別研究層狀和不同波阻抗組合條件下的火山巖理論模型的地震響應(yīng)特征，認(rèn)為火山巖的存在導(dǎo)致下伏地層的上翹，且對下伏地層成像具有屏蔽作用，進(jìn)而對地震解釋的準(zhǔn)確性產(chǎn)生影響。蘇巍等引入非平穩(wěn)隨機介質(zhì)模型構(gòu)建方法，建立丘狀火山巖地質(zhì)模型，并采用交錯網(wǎng)格有限差分法進(jìn)行正演模擬，分析由淺至深各反射層及沉積相帶的地震反射特征，認(rèn)為火山體的頂、底界面振幅強弱及反射形態(tài)，與火山巖體的厚度和上、下圍巖的巖性組合有關(guān)[18]。黃林軍等利用Tesseral 2.5D正演軟件，對兩個不同密度和速度的層狀火山巖地質(zhì)模型進(jìn)行數(shù)值模擬分析，論證火山巖儲層地震響應(yīng)特征為連續(xù)低速強振幅反射，同時預(yù)測風(fēng)城組下部另一套火山巖儲層的存在[19]。

目前，有關(guān)火山巖儲層的正演模擬研究存在兩點不足：一是有關(guān)蘑菇狀火山巖的正演模擬研究較少，由于火山巖的噴發(fā)模式及噴發(fā)時期不同，使得火山巖的外部形態(tài)復(fù)雜多變，常見的外部形態(tài)有菱形、橢球形、丘形、蘑菇狀、層狀及楔狀等，其中蘑菇狀的火山巖按照地震相分類的劃分屬于蘑菇狀地震相。通常蘑菇狀地震相的內(nèi)部及其上部的披覆構(gòu)造中易形成氣藏和油藏，是一種典型的火山巖油氣藏[14]，因此有必要對蘑菇狀火山巖進(jìn)行正演模擬，分析地震響應(yīng)特征，增強對蘑菇狀火山巖儲層的認(rèn)識。二是有關(guān)火山巖的正演模擬是基于地震數(shù)值模擬展開的，缺少物理模擬研究。數(shù)值模擬常常受計算方法和假設(shè)條件的限制，使得正演模擬達(dá)不到理想效果，而地震物理模擬能較真實地模擬野外實際地質(zhì)情況，在地震勘探領(lǐng)域得到廣泛的應(yīng)用[20-27]。

筆者利用地質(zhì)資料，建立蘑菇狀火山巖地質(zhì)模型，分別進(jìn)行地震物理模擬和數(shù)值模擬；兩種正演模擬方法除了具有各自的優(yōu)勢外，還具有一定的互補性，將兩種技術(shù)結(jié)合可以更好地認(rèn)識蘑菇狀火山巖的地震響應(yīng)特征，為火山巖地震響應(yīng)特征的研究提供思路。

1 火山巖地質(zhì)模型

根據(jù)某火山巖地區(qū)地質(zhì)解釋、地震資料和測井資料等信息，構(gòu)建蘑菇狀火山巖地質(zhì)模型(見圖1)。該模型中火山巖的外部形態(tài)為蘑菇狀，火山通道與其圍巖之間夾角較大，近似垂直，蘑菇頭頂部上覆地層為白云巖，下伏地層為沉積砂礫巖。模型長度為20.00 km，深度為9.08 km，每層的地震波速度和巖石密度參數(shù)見表1。由于物理模型的采集是在水槽中進(jìn)行的，將該模型的頂?shù)自O(shè)計為相應(yīng)厚度的水層。

圖1 蘑菇狀火山巖地質(zhì)模型Fig.1 Mushroom-shaped volcanic geological model

根據(jù)火山巖地質(zhì)模型，制作相應(yīng)的火山巖物理模型。由于物理模型在制作過程中受到制作材料特性的限制，使得制作的地層密度和速度參數(shù)不能完全與火山巖地質(zhì)模型中的密度和速度參數(shù)相同，需要通過一定的比例轉(zhuǎn)換后才能與地質(zhì)模型接近，最終確定物理模型與火山巖地質(zhì)模型的密度比為1∶1.61，速度比為1∶2，轉(zhuǎn)換前后的密度和速度參數(shù)見表1。蘑菇狀火山巖實物圖及制作完成的火山巖模型實物圖見圖2，物理模型的長度和寬度為1 000.0 mm，高度為334.4 mm。根據(jù)物理模型參數(shù)，制作與其參數(shù)相同的數(shù)值模擬模型。

表1 火山巖地質(zhì)模型參數(shù)

圖2 蘑菇狀火山巖實物和物理模型Fig.2 The material object and physical model of mushroom-shaped volcanic rock

2 數(shù)據(jù)采集

2.1 物理模擬

物理模型數(shù)據(jù)采集由高精度自動定位采集系統(tǒng)完成。設(shè)計二維地震數(shù)據(jù)采集的觀測系統(tǒng)：單邊放炮，共194炮，炮間距為80 m，每炮216道，道間距為40 m，最小偏移距為400 m，采樣間隔為2 ms，采用的震源子波是中心頻率約為230 kHz的多相位混合子波。

第80炮單炮的物理和數(shù)值模擬記錄見圖3(a)，可以看出物理模擬記錄中地震波信息豐富，除了直達(dá)波和一次反射波外含有大量的多次波，記錄上標(biāo)出各反射界面的同相軸，因為存在高速火山巖體對下伏地層的屏蔽作用及多次波的干擾，反射界面很難在記錄中被準(zhǔn)確地分辨出來。單炮記錄產(chǎn)生多次波的原因是物理模擬數(shù)據(jù)采集是在水槽中完成的，模型頂部水層使得地震波從地下界面反射到水層后，有一部分波再次發(fā)生反射轉(zhuǎn)而向地下傳播，最終產(chǎn)生多次波。

圖3 第80炮單炮記錄的物理模擬和數(shù)值模擬記錄Fig.3 The 80th shot record of physical modeling data and numerical simulation data

多次波是一種干擾波，它與一次反射波互相干涉迭加，破壞有效波的識別和追蹤，進(jìn)而影響目的層的成像，因此需要對多次波進(jìn)行壓制。分別測試聚束、拉冬、F-K濾波等壓制多次波方法，結(jié)合物理模擬數(shù)據(jù)特點，選擇拉冬濾波和內(nèi)切結(jié)合的方法對多次波進(jìn)行壓制，多次波衰減前、后速度譜見圖4。由圖4可以看出，多次波衰減后的速度譜質(zhì)量明顯變好，能夠提高速度分析的精度。

圖4 多次波衰減速度譜及NMO道集Fig.4 The velocity spectrum and NMO trace set before and after multiple wave attenuation

2.2 數(shù)值模擬

采用WaveMOD二維正演模擬軟件對火山巖模型進(jìn)行數(shù)值模擬[28]。該軟件具有計算效率高、占用內(nèi)存小等特點，適用于非均勻介質(zhì)、復(fù)雜構(gòu)造及復(fù)雜巖性介質(zhì)的地震波場模擬。鑒于物理模型數(shù)據(jù)采集以縱波為主，數(shù)值模擬采用聲波方程模擬，其觀測系統(tǒng)與物理模擬的相同。綜合考慮模型的尺寸及計算效率，使用20 m網(wǎng)格對數(shù)值模型進(jìn)行網(wǎng)格剖分，激發(fā)震源采用主頻為20 Hz的雷克子波。

第80炮單炮記錄見圖3(b)，在單炮記錄中標(biāo)出各反射界面的同相軸。與圖3(a)相比，各反射面上的同相軸位置基本相同，但在能量上有差別，直達(dá)波能量較強，在遠(yuǎn)偏移距處，物理模擬數(shù)據(jù)能量的衰減明顯比數(shù)值模擬的快。另外，數(shù)值模擬的地震波信息沒有物理模擬的豐富。

3 正演模擬結(jié)果

對物理模擬和數(shù)值模擬的地震記錄進(jìn)行常規(guī)資料處理，其處理流程基本一致，不同的是物理模擬數(shù)據(jù)需要進(jìn)行多次波壓制處理；然后對兩種數(shù)據(jù)進(jìn)行疊前時間偏移和疊前深度偏移成像，在成像過程中兩種數(shù)據(jù)的偏移參數(shù)也基本一致，成像剖面見圖5和圖6。在物理模擬和數(shù)值模擬偏移剖面上，蘑菇狀火山巖的反射特征整體上是一致的，火山巖上部外形為蘑菇頭狀，與其結(jié)構(gòu)形態(tài)類似；火山巖下部通道在地震剖面上沒有得到相應(yīng)的反射特征，產(chǎn)生“蘑菇無根”的現(xiàn)象，但兩者反射特征在細(xì)節(jié)上存在一定的差異。

3.1 蘑菇頭邊界及內(nèi)部反射特征

蘑菇頭頂部邊界(見圖5和圖6中①)同相軸成像不清晰，并且同相軸連續(xù)性較差，原因是上覆地層與火山巖的波阻抗差異較小，反射因數(shù)為0.035 6。由第121炮(t=0.65 s)的波場快照(見圖7(a)中①)可以看出，蘑菇頭頂部的反射很弱，導(dǎo)致頂界面同相軸成像不清晰。相比蘑菇頭頂部的弱反射特征，蘑菇頭的兩個下側(cè)邊界(見圖5和圖6中②)同相軸能量很強，表現(xiàn)為連續(xù)性好、反射特征強，主要原因是蘑菇頭下伏地層與高速火山巖的波阻抗差異較大，使得同相軸振幅強度增強。由第121炮(t=0.75 s)的波場快照(見圖7(b)中①)可以看出，下側(cè)邊界的反射很強，是蘑菇狀火山巖外部形態(tài)識別的重要標(biāo)志，為地震層位標(biāo)定提供較好的依據(jù)。

圖5 物理模擬的疊前時間和疊前深度偏移剖面Fig.5 The prestack time and prestack depth migration section of physical modeling

蘑菇頭兩翼的上側(cè)邊界(見圖5和圖6中③)同相軸基本上未成像，原因是T5層和T6層的波阻抗差大而產(chǎn)生強反射，對下伏地層有屏蔽作用，影響地震波的向下傳播。由第121炮(t=0.65 s)的波場快照(見圖7(a)中②)可以看出，地震波沒有反射而影響邊界成像。當(dāng)?shù)卣鸩ɡ^續(xù)傳播，由第121炮(t=0.75 s)的波場快照(見圖7(b)中②)可以看出，地震波不再受上覆地層的屏蔽作用影響而產(chǎn)生反射。

火山巖上部即蘑菇頭部分，其內(nèi)部反射特征在偏移剖面上表現(xiàn)為明顯的空白反射(見圖5和圖6)，與正演模擬方法和成像方法無關(guān)，因此利用蘑菇頭頂部邊界的弱反射、下側(cè)邊界的強反射，以及蘑菇頭內(nèi)部的空白反射特征可以在地震剖面上確定火山巖上部的外部形態(tài)和空間展布。

3.2 下伏地層反射特征

在物理模擬和數(shù)值模擬的疊前時間偏移剖面上，火山巖底部的同相軸(見圖5(a)和圖6(a)中④)出現(xiàn)嚴(yán)重的上凸現(xiàn)象。原因是火山巖具有高速特性，使得下伏地層的反射波穿過火山巖體時反射旅行時變短，導(dǎo)致火山巖下伏地層的反射旅行時也變短。在兩個疊前深度偏移剖面(見圖5(b)和圖6(b))中未見到上凸現(xiàn)象，說明上凸現(xiàn)象是由疊前時間偏移方法的固有缺點引起的。當(dāng)速度存在橫向劇烈變化或者高陡構(gòu)造時，疊前時間偏移繞射能量和反射能量不能準(zhǔn)確地聚焦和歸位，進(jìn)而產(chǎn)生上凸現(xiàn)象，因此在利用時間偏移結(jié)果對火成巖進(jìn)行解釋時，應(yīng)該注意剩余繞射對火山巖反射特征的影響。

3.3 蘑菇根部反射特征

在物理模擬和數(shù)值模擬的疊前深度偏移剖面上，火山巖的下部通道在地震剖面上產(chǎn)生“蘑菇無根”現(xiàn)象，但“蘑菇無根”現(xiàn)象的程度不同。在數(shù)值模擬的疊前時間和深度偏移剖面上，蘑菇根部表現(xiàn)為空白反射樣；在物理模擬的疊前時間偏移剖面上，蘑菇根的上部沒有反射，蘑菇根的下部出現(xiàn)的反射也無法與實際的結(jié)構(gòu)對應(yīng)；在物理模擬的疊前深度偏移剖面上，雖然沒有產(chǎn)生明顯的蘑菇根的反射特征，但可以看到與蘑菇根部的結(jié)構(gòu)特征對應(yīng)的異常反射(見圖5(b))。

對“蘑菇無根”現(xiàn)象進(jìn)行分析，選取第121炮(t=0.85、1.10 s)的波場快照(見圖7(c)、(d))。由圖7(c)、(d)可以看出，由于整個火山通道的構(gòu)造陡峭而近似垂直，從地震波開始進(jìn)入火山通道到傳播到火山通道的底部時，火山通道內(nèi)接收到的地震波的能量很弱。另外，根據(jù)模型設(shè)計參數(shù)，火山通道與兩側(cè)的地層T8、T9的波阻抗差異很小，接觸面的反射因數(shù)分別為0.109 0和0.046 0；因此即使在火山通道內(nèi)可以接收到地震波傳來的能量，在火山通道的兩側(cè)邊界也很難有充足的反射。雖然通過疊前深度偏移成像可以有效地改善高陡構(gòu)造的成像問題，但是火山通道與兩側(cè)地層的波阻抗差異很小，在數(shù)值模擬的剖面上經(jīng)過深度偏移處理也很難對蘑菇根部的構(gòu)造進(jìn)行成像。

此外，在模型制作和信號采集過程中，物理模擬可能給地震記錄帶來一些噪音，在進(jìn)行疊前時間偏移處理時噪音對蘑菇根部的成像產(chǎn)生影響。在物理模擬的疊前深度偏移剖面上，在蘑菇根部有明顯的異常反射特征，與數(shù)值模擬的疊前深度偏移剖面相比，兩個剖面在蘑菇根部有同樣的難以成像原因，但兩者反射特征有很大差異。因為物理模擬得到的二維地震記錄是從三維物理模型中采集到的，數(shù)值模擬得到的二維地震記錄是根據(jù)二維數(shù)學(xué)模型經(jīng)過波動方程正演而得到的，即兩者地震波在地層接觸面處的傳播規(guī)律存在差異，前者更接近于野外的傳播方式，后者一般按照嚴(yán)格的數(shù)學(xué)條件進(jìn)行傳播，所以即使在構(gòu)造陡峭且波阻抗差異很小的情況下，在物理模擬的疊前深度偏移剖面上，也可以看到與實際火山巖結(jié)構(gòu)特征接近的異常反射。對于存在構(gòu)造陡峭且構(gòu)造與兩側(cè)地層之間的波阻抗差異比較小的地質(zhì)構(gòu)造，使用二維的數(shù)值模擬時得到的地震反射特征可能與實際的反射特征不符，即二維的數(shù)值模擬存在一定的局限性。這種情況下應(yīng)該使用三維正演模擬，或者根據(jù)三維模型進(jìn)行二維正演模擬，得到的反射特征才更真實可靠。

3.4 精度分析

為了對比物理模擬和數(shù)值模擬的精度，分別抽取火山巖體兩翼尖滅點處的CDP道集，火山巖成像的寬度分別為4.92 km和5.24 km，實際設(shè)計寬度為5.16 km，數(shù)值模擬的成像寬度更接近于實際寬度，還要考慮物理模型在制作過程中存在1 mm左右的誤差。因此，物理模擬的精度雖然沒有數(shù)值模擬的高，但也可以達(dá)到理想的精度。

4 結(jié)論

(1)對蘑菇狀的火山巖模型分別進(jìn)行地震物理模擬和數(shù)值模擬，其疊前時間和疊前深度偏移剖面上蘑菇頭邊界及內(nèi)部的反射特征基本一致，利用蘑菇頭頂部邊界的弱反射、下側(cè)邊界的強反射，以及蘑菇頭內(nèi)部的空白反射特征可以在地震剖面上初步確定蘑菇狀火山巖蘑菇頭的外部形態(tài)和空間展布。

(2)在疊前時間偏移剖面上，火山通道的下伏地層出現(xiàn)嚴(yán)重的上凸現(xiàn)象，與偏移成像方法有關(guān)，可以通過疊前深度偏移避免。在利用疊前時間偏移結(jié)果對火山巖進(jìn)行解釋時，需要注意剩余繞射對火山巖反射特征的影響。兩種模擬結(jié)果的偏移剖面上產(chǎn)生“蘑菇無根”的現(xiàn)象，但是“蘑菇無根”的程度及產(chǎn)生的原因并不全完相同。

(3)在一定情況下，二維的物理模擬相比二維的數(shù)值模擬能更真實地模擬地下情況，地震波信息也更加豐富。數(shù)值模擬的資料品質(zhì)好，并且可以利用波場信息進(jìn)一步解釋正演得到的地層反射特征。兩種模擬方法可以相互補充，在進(jìn)行火山巖的地震正演模擬時可以綜合運用物理模擬和數(shù)值模擬。

[1] 馮子輝,朱映康,張元高,等.松遼盆地營城組火山機構(gòu)—巖相帶的地震響應(yīng)[J].地球物理學(xué)報,2011,54(2):556-562.

Feng Zihui, Zhu Yingkang, Zhang Yuangao, et al. The sseismic response of the volcanic edifice-facies zone of Yingcheng formation in the Songliao basin [J]. Chinese Journal of Geophysics, 2011,54(2):556-562.

[2] 王力圓,唐華風(fēng),王璞珺,等.松遼盆地長嶺斷陷火山巖相地震特征及其響應(yīng)關(guān)系[J].地球物理學(xué)進(jìn)展,2012,27(1):271-278.

Wang Liyuan, Tang Huafeng, Wang Pujun, et al. The seismic reflection characteristics and response relationship of volcanic rocks lithofacies in Changling fault depression Songliao basin [J]. Progress in Geophysics, 2012,27(1):271-278.

[3] 王志萍,呂丁友,王保全.渤中地區(qū)火山巖特征及其與斷裂關(guān)系[J].特種油氣藏,2015,22(4):20-24.

Wang Zhiping, Lv Dingyou, Wang Baoquan. Volcanic rock characteristics and its relationship with the rupture in the Bozhong area [J]. Special Oil and Gas Reservoirs, 2015,22(4): 20-24.

[4] 陳永波,潘建國,許多年,等.準(zhǔn)噶爾盆地西北緣火山巖儲層的綜合地球物理預(yù)測[J].石油物探,2010,49(4):364-372.

Chen Yongbo, Pan Jianguo, Xu Duonian, et al. Integrated geophysical prediction of volcanic reservoir in northwestern margin area of Junggar basin [J]. Geophysical Prospecting for Petroleum, 2010,49(4):364-372.

[5] 王文廣,鄭民,盧雙舫,等.徐家圍子斷陷火山巖氣藏主控因素及有利區(qū)[J]. 東北石油大學(xué)學(xué)報,2015,39(4):45-53.

Wang Wenguang, Zheng Min, Lu Shuangfang, et al. Main controlling factors and favorable area prediction of gas reservoir of deep volcanic rock in Xujiaweizi fault depression [J]. Journal of Northeast Petroleum University, 2015,39(4):45-53.

[6] 商豐凱,陳林,王林,等.車排子凸起火山巖油藏成藏主控因素和模式[J].東北石油大學(xué)學(xué)報,2015,39(5):13-22.

Shang Fengkai, Chen Lin, Wang Lin, et al. Main controlling factors and hydrocarbon accumulation modes of volcanic rock oil reservoirs in Chepaizi uplift [J]. Journal of Northeast Petroleum University, 2015,39(5):13-22.

[7] 張子樞,吳邦輝.國內(nèi)外火山巖油氣藏研究現(xiàn)狀及勘探技術(shù)[J].天然氣勘探與開發(fā),1994,16(1):1-26.

Zhang Zishu, Wu Banghui. Research status and exploration technique investigation of volcanic oil & gas reservoirs at home and abroad [J]. Natural Gas Exploration and Development, 1994,16(1):1-26.

[8] 董亦思,黃文彪,于興河,等.英臺斷陷火山巖氣藏特征及有利區(qū)預(yù)測[J].東北石油大學(xué)學(xué)報,2015,39(6):30-37.

Dong Yisi, Huang Wenbiao, Yu Xinghe, et al. Characteristics and main control factors of volcanic gas reservoirs in Yingtai fault depression [J]. Journal of Northeast Petroleum University, 2015,39(6):30-37.

[9] 石磊,李書兵,黃亮,等.火山巖儲層研究現(xiàn)狀與存在的問題[J].西南石油大學(xué)學(xué)報:自然科學(xué)版,2009,31(5):68-72.

Si Lei, Li Shubing, Huang Liang, et al. Current research situation and the problems of volcanic rock reservoirs [J]. Journal of Southwest Petroleum University: Science & Technology Edition, 2009,31(5):68-72.

[10] 王璞珺,張功成,蒙啟安,等.地震火山地層學(xué)及其在我國火山巖盆地中的應(yīng)用[J].地球物理學(xué)報,2011,54(2):597-610.

Wang Pujun, Zhang Gongcheng, Meng Qian, et al. Applications of seismic volcano stratigraphy to the volcanic rifted basins of China [J]. Chinese Journal of Geophysics, 2011,54(2):597-610.

[11] 謝春臨,陳樹民,姜傳金,等.趨勢面分析與相干體技術(shù)在火山巖預(yù)測中的應(yīng)用[J].地球物理學(xué)報,2011,54(2):368-373.

Xie Chunlin, Chen Shumin, Jiang Chuanjin, et al. Application of a method combining trend surface analysis and coherence cube technology to volcanic prediction [J]. Chinese Journal of Geophysics, 2011,54(2):368-373.

[12] 王春燕,張雷,金雪英.蝕變火山巖儲層特征及氣水層識別方法[J].東北石油大學(xué)學(xué)報,2013,37(3):18-23.

Wang Chunyan, Zhang Lei, Jin Xueying. Characteristics and fluid identification of altered volcanic reservoirs [J]. Journal of Northeast Petroleum University, 2013,37(3):18-23.

[13] 王偉鋒,高斌,衛(wèi)平生,等.火山巖儲層特征與油氣成藏模式研究[J].地球物理學(xué)進(jìn)展,2012,27(6):2478-2491.

Wang Weifeng, Gao Bin, Wei Pingsheng, et al. Research of volcanic reservoir characters and hydrocarbon accumulation models [J]. Progress in Geophysics, 2012,27(6):2478-2491.

[14] 孫淑艷,李艷菊,彭莉,等.火成巖地震識別及構(gòu)造描述方法研究[J].特種油氣藏,2003,10(1):47-50.

Sun Shuyan, Li Yanju, Peng Li, et al. Study on seismic identification and structure description of igneous rock [J]. Secial oil and Gas Reservoirs, 2003,10(1):47-50.

[15] 王建花,李慶忠,邱睿.淺層強反射界面的能量屏蔽作用[J].石油地球物理勘探,2003,38(6):589-596.

Wang Jianhua, Li Qingzhong, Qiu Rui. Energy shielding action of shallow strong reflector [J]. Oil Geophysical Prospecting, 2003,38(6):589-596.

[16] 毛小平,何大偉,辛廣柱,等.厚層火山巖地震響應(yīng)研究[J].石油地球物理勘探,2005,40(6):670-676.

Mao Xiaoping, He Dawei, Xin Guanzhu, et al. Study of seismic response for thick-layered volcanic rock [J]. Oil Geophysical Prospecting, 2005,40(6):670-676.

[17] 李素華,王云專,盧齊軍,等.火山巖波動方程正演模擬研究[J].石油物探,2008,47(4):361-366.

Li Suhua, Wang Yunzhuan, Lu Qijun, et al. Wave equation forward modeling for volcanic rock [J]. Geophysical Prospecting for Petroleum, 2008,47(4):361-366.

[18] 蘇巍,張春芬.火山巖地震響應(yīng)特征的正演模擬[J].世界地質(zhì),2011,30(3):429-434.

Su Wei, Zhang Chunfen. Forward modeling of volcanic seismic response characteristics [J]. Global Geology, 2011,30(3):429-434.

[19] 黃林軍,楊巍,王彥軍,等.模型正演技術(shù)在火山巖儲層識別中的應(yīng)用——以準(zhǔn)噶爾盆地烏夏地區(qū)二疊系火山巖儲層為例[J].天然氣地球科學(xué),2011,22(3):539-542.

Huang Linjun, Yang Wei, Wang Yanjun, et al. Application of forward solution model in identifying volcanic reservoir: A case of permian volcanic reservoir in Wuxia area of Junggar basin [J]. Natural Gas Geoscience, 2011,22(3):539-542.

[20] Pant D R, Greenhalgh S A, Watson S. Swismic reflection scale model facility [J]. Exploration Geophysics, 1989,19:499-512.

[21] Yike Mo, Stewart A G, Johan O A, et al. The development and testing of a 2D laboratory seismic modeling system for heterogeneous structure investigations [J]. Journal of Applied Geophysics, 2015,116:224-235.

[22] 孫進(jìn)忠,郭鐵栓,唐文榜,等.我國超聲地震模型試驗的理論研究與實踐[J].地球物理學(xué)報,1997,40(S1):266-274.

Sun Jinzhong, Guo Tieshuan, Tang Wenbang, et al. Theory practice and application of ulltrasonic seismic modeling experiment in China [J]. Chinese Journal of Geophysics, 1997,40(S1):266-274.

[23] 尹志恒,李向陽,魏建新,等.利用物理模型三維縱波數(shù)據(jù)分析HTI介質(zhì)的方位各向異性[J].地球物理學(xué)報,2012,55(11):3805-3812.

Yin Zhiheng, Li Xiangyang, Wei Jianxin, et al. A physical modeling study on the 3D P-wave azimuthal anisotropy in HTI media [J]. Chinese Journal of Geophysics, 2012,55(11):3805-3812.

[24] 司文朋,狄?guī)妥?魏建新.地震物理模型實驗數(shù)據(jù)中多次轉(zhuǎn)換干擾波的分析[J].石油物探,2015,54(1):17-23.

Si Wenpeng, Di Bangrang, Wei Jianxin. Analysis of multiple-converted interference wave in seismic physical modeling data [J]. Geophysical Prospecting for Petroleum, 2015,54(1):17-23.

[25] Wu M S, Di B R, Wei J X, et al. Lark-scale complex physical modeling and precision analysis [J]. Applied Geophysics, 2014,11(2):245-251.

[26] Wang S X, Li X Y, Di B R. Reservoir fluid substitution effects on seismic profile interpretation: A physical modeling experiment [J]. Geophysical Research Letters, 2010,37(10):10306.

[27] 李凡異,魏建新,狄?guī)妥?碳酸鹽巖溶洞橫向尺度變化的地震響應(yīng)正演模擬[J].石油物探,2009,48(6):557-562.

Li Fanyi, Wei Jianxin, Di Bangrang. Forward simulation of seismic response in carbonate caverns with lateral scale [J]. Geophysical Prospecting for Petroleum, 2009,48(6):557-562.

[28] 裴正林.三維各向同性介質(zhì)彈性波方程交錯網(wǎng)格高階有限差分法模擬[J].石油物探,2005,44(4):308-315.

Pei Zhenglin. Numerical simulation of elastic wave equation in 3D isotropic media with staggered-grid high-order difference method [J]. Geophysical Prospecting for Petroleum, 2005,44(4):308-315.

2016-08-03；編輯：任志平

國家科技重大專項(2016ZX05007-006)

李東慶(1988-)，男，博士研究生，主要從事地震物理模擬方面的研究。

TE19;P631.4

A

2095-4107(2016)05-0094-08

DOI 10.3969/j.issn.2095-4107.2016.05.011