火山巖發育區地震成像技術探討

摘要：隨著油氣勘探程度的不斷提高，已陸續發現一些以火山巖和火山碎屑巖為儲集層的油氣田。但是，由于火山巖對地震信號屏蔽作用明顯，地震波在不同層系、不同類型火山巖中的速度差異大，導致火山巖發育區的速度建模和成像極其困難。EP地區處于陽江—一統斷裂帶，資源潛力大，且發育大量構造圈閉，同時，該地區發育古近系及上覆地層的多期火山巖，以往處理成果中火山巖內幕成像不清晰，難以刻畫火山巖相和準確評估火山巖對油氣生、儲、蓋條件的影響。為此，針對火山巖多次波壓制、火山巖下弱信號恢復、火山巖速度精細刻畫等方面的關鍵問題開展新的處理技術探索。首先采用MWD（Model based Water-layer related Demultiple）技術預測火山巖多次波；然后利用AVO剩余振幅補償和低頻算子外推隨機噪聲衰減技術提高火山巖下弱信號質量；最后，聯合應用初至波層析和反射波網格層析技術構建精確的疊前深度偏移速度模型。重處理后的疊前深度偏移成果顯示，火山巖及下伏地層成像質量顯著提高，表明所提重處理技術方案對提高火山巖發育區地震成像質量有較好的借鑒意義。

關鍵詞：火山巖，多次波，SRME，MWD，初至波層析反演，網格層析反演

中圖分類號：P631文獻標識碼：A DOI：10.13810/j.cnki.issn.1000-7210.2024.04.017

Seismic imaging technology in volcanic rock areas

ZHANG Weiwei1，YUAN Quanshe2，XU Leyi1，ZHANG Zhizhong1，GUO Jia1，GUO Jianqing3

（1.Nanhai East Petroleum Research Institute，CNOOC Limited Shenzhen Branch，Shenzhen，Guangdong 518067，China；

2.China National Offshore Oil Corporation，Beijing 100005，China；3.GRIDagang Branch，BGP，CNPC，Tianjin 300280，China）

Abstract：With the continuous improvement of oil and gas exploration，some oil and gas fields with volcanic rocks and volcaniclastic rocks as reservoirs have been discovered.However，due to the significant shielding ef?fect of volcanic rocks on seismic signals，the significant velocity variations of seismic waves in different layers and types of volcanic rocks，velocity modeling and imaging in volcanic rock development areas are extremely dif?ficult.The EP region is located in the Yangjiang?Yitong Fault Zone，with great resource potential and a large number of structural traps.Meanwhile，the Paleogene and its overlying strata have developed multiple periods and layers of volcanic rocks.Due to unclear imaging of the volcanic rocks and their underlying strata，it is diffi?cult to give a good definition of volcanic facies and to accurately evaluate the impact of volcanic rocks on the con?ditions of hydrocarbon generation，storage，and capping.This paper explores some new methods and new tech?nology attempted for key issues such as multiple wave suppression，weak signal recovery under volcanic rocks，and precise velocity modeling of volcanic rocks.Firstly，MW D（Model based Water-layer related Demultiple）technology is used to predict volcanic multiples.Then，AVO residual amplitude compensation and low?fre?quency operator extrapolation random noise attenuation techniques are used to improve the quality of weak sig?nals under volcanic rocks.Finally，first break tomography and grid?based tomography are used for accurate PSDM（Prestack Deep Migration）velocity modeling.The reprocessed PSDM results demonstrate substantial improvement in imaging of volcanic rocks and their underlying strata，indicating that the reprocessing techno?logy scheme can be a good reference for seismic imaging in volcanic rock areas.

Keywords：volcanic rocks，multiple waves，SRME，MW D，first break tomography，grid?based tomography

張衛衛，袁全社，徐樂意，等.火山巖發育區地震成像技術探討［J］.石油地球物理勘探，2024，59（4）：800-810.

ZHANG Weiwei，YUAN Quanshe，XU Leyi，et al.Seismic imaging technology in volcanic rock areas［J］.Oil Geophysical Prospecting，2024，59（4）：800-810.

0引言

火山巖指高溫熔融的巖漿在地下或噴出地表后冷凝而成的巖石[1]。火山巖和火山碎屑巖在各時代地層中都有廣泛分布，已陸續發現一些以其為儲集層的油氣田[2]。

長期以來，火山巖及下伏地層的成像問題一直是火山巖發育區地震資料處理的巨大挑戰。火山巖主要在兩方面影響地震資料處理：一是受火山巖屏蔽的影響，火山巖下伏地層反射能量弱，噪聲淹沒了有效信號，難以有效識別地層真實的反射特征；二是火山巖發育區速度的劇烈變化導致地震波場復雜、時間域構造成像畸變，使速度建模和成像極為困難[3]。地球物理工作者針對上述兩方面挑戰進行了長期探索。黎書琴等[4]通過高精度三維觀測系統優化設計和低頻檢波器的應用，提高了火山巖地區地震數據品質。彭浩天等[5]通過低頻補償和Q偏移技術，以提高火山巖成像質量。韓站一等[6]從Zoep?pritz方程出發，分析了火成巖地層中地震波的透射損失，通過模型正演揭示了火成巖非均質性對透射能量的衰減作用，認為火成巖對下伏地層的影響因素除了高速屏蔽外還有介質非均質性導致的散射問題。姜巖等[7]根據廣義O'Doherty-Anstey公式，探索了一種利用聲波曲線求取透射系數的方法，以補償強波阻抗界面對下伏地層反射振幅的影響。張濤等[8-9]根據火山巖特殊的地質特征，利用VSP測井信息為約束條件，解決三維網格層析反演速度建模的多解性問題。郭樹祥[10]根據火山巖的波場特征，正演模擬了火山巖地震波場頻散以及低通濾波效應。王興等[11]利用全區的聲波測井速度平滑插值得到初始速度模型，以減小火山巖速度誤差。梁海龍等[12]基于高密度地震數據，嘗試利用Q疊前深度偏移技術提高火山巖地層的成像精度。這些研究成果大多依賴原始資料和鉆井資料的品質，在工區鉆井資料不足或火山巖空間分布較復雜時，仍然難以有效解決火山巖及下伏地層的成像問題。

本文依托EP地區地震采集資料研究火山巖發育區的成像處理技術。針對長期困擾該地區的火山巖及下伏構造成像瓶頸問題，在火山巖多次波壓制、火山巖下弱信號恢復、火山巖速度精細刻畫三方面進行了新的探索，有效改善了火山巖及下伏地層成像質量。

1研究區概況

EP地區處于陽江—一統斷裂帶，資源潛力大，且發育有大量構造圈閉，是油氣勘探有利的戰略接替區。該地區古近系及上覆地層發育多期火山巖，且形態復雜，明顯屏蔽地震信號；火山巖多次波對火山巖下伏地層反射的干擾作用強；同時，不同層系、不同類型火山巖速度差異大。而且，工區位于水深50～90 m的淺海地區，反射波在水層中反復振蕩產生極強能量的多次波，掩蓋了火山巖下伏地層真實的反射特征。另外，工區地震數據為常規的窄方位拖纜數據，目前僅有的一口探井也未在火山巖區域。

火山巖對地震數據的影響導致該地區地震成像品質差，嚴重制約了油氣勘探、開發工作的進程。圖1為工區以往疊前深度偏移成果的時間域剖面。在750～3000 ms的時間范圍內，可見明顯的不同層系火山巖強振幅特征，但在藍色圓圈區域，復雜火山巖體的邊界不清楚；在火山巖下，箭頭所示位置地層和控洼斷裂成像不清楚，導致地層產狀難以落實，深部疑似半地塹的凹陷結構難以識別。

2關鍵技術及應用效果

本文在針對拖纜地震數據特點完成涌浪噪聲壓制等基本的前置處理后，根據EP火山巖發育區的地震資料特點，針對制約地震成像的關鍵問題，制定了圖2所示的地震成像處理技術方案。首先采用三維MWD（Model based Water-layer related Demul?tiple）技術提高火山巖微屈多次波預測效果；然后，通過AVO剩余振幅補償及低頻算子外推隨機噪聲衰減技術合理補償地震波的透射損失，并提高火山巖下伏地層信噪比；最后，綜合應用回折波層析反演和反射波網格層析，精細刻畫火山巖導致的速度縱、橫向變化，為疊前深度偏移成像奠定了基礎。

2.1多次波壓制

2.1.1主要多次波類型

EP地區水深為50～90 m，自由表面（海水面）、海底、火山巖頂是工區最主要的三個強反射界面。與這三個強反射界面對應的三種多次波（圖3）：第一種傳播路徑如紅色箭頭所示，海底反射波在水層中反復振蕩形成海底的水層多次波；第二種為黃色箭頭所示的火山巖全程多次波；第三種為綠色箭頭所示，火山巖一次反射在海水層反復振蕩，形成與火山巖和水層相關的火山巖微屈多次波。其中，火山巖全程多次波由于反射周期長，多次波旅行時大于目的層的一次波旅行時，反射能量大多位于主要目的層以下，對地震成像的影響較小。而且，由于多次波僅與火山巖和海水面相關，與海底反射無關，SRME（Surface Related Multiple Elimination）等自由表面多次波預測技術可以較好地構建多次波模型。而海底的水層多次波和火山巖微屈多次波是一次波在海底與自由表面之間的水層反復振蕩產生的多次反射，當海底反射缺失時，SRME等依賴地震數據褶積預測多次波的方法難以發揮作用，需要采用新的方法預測多次波模型。

三種多次波全部和自由表面相關，理論上可以采用SRME等數據驅動的自由表面多次波預測方法預測多次波模型。然而，在水深非常淺的海域，SRME的應用效果往往不理想。首先，海洋拖纜采集的最小炮檢距一般在200 m左右，因為缺少近炮檢距數據，大部分水底反射信號缺失，接收到的信息主要來自臨界角之外的水底反射能量，對較小炮檢距的水層多次波沒有貢獻；其次，由于觀測到的一次水底反射被直達波干擾，因而不能正確重構水底一次波[13]。

正演記錄能夠揭示近道水底反射信號缺失對SRME多次波預測的影響。選擇一個二維正演模擬觀測系統建立多次波正演模型（圖4）。從圖5a的正演記錄可以發現，在火山巖一次反射（紅色箭頭）下，存在明顯的火山巖反射在水層反復振蕩形成的短周期多次波（藍色箭頭），而海底近道反射信息的缺失嚴重制約了SRME對火山巖多次波的預測效果（圖5b）。為彌補SRME技術在淺水區多次波預測的不足，本文采用模型驅動的三維MWD多次波預測方法預測火山巖微屈多次波。

2.1.2 MWD多次波預測技術原理

MWD是一種模型驅動的水層相關多次波預測技術，首先利用已知的海底曲面模型構建海底的格林函數，然后將接收地震記錄與格林函數褶積，恢復淺水環境的多次波信息。

理論上，如圖3所示的火山巖微屈多次波S→R可以通過地震數據與水層一次反射的褶積得到[14]

M（S，R，ω）=∫D（S，Ri，ω）P（R，ω）dR i（1）

式中：M是多次波模型；D為地震數據；P為海底一次反射，具有多次波預測算子的作用；ω為角頻率。

為解決淺水區因海底近道反射缺失導致的多次波預測問題，可以利用海底測深信息計算水層格林函數，用格林函數代替水層的一次反射預測多次波

M（S，R，ω）=D r（S，Ri，ω）G（R，Ri，ω）（2）

式中：G為水層的格林函數；Ω為多次波預測孔徑。

由于利用格林函數代替海底一次反射和地震數據褶積，MWD擺脫了對海底反射的依賴，從而解決了淺水區海底近道反射信息缺失導致的多次波預測問題。圖6的正演記錄顯示了MWD在淺水區預測多次波的優勢。對于最小炮檢距為200 m的正演記錄（圖6a），MWD獲得的火山巖多次波模型（圖6b）和正演記錄的多次波（圖6a）反射特征非常相似。

圖7為EP地區SRME與MWD多次波預測獲得的多次波模型對比。在原始疊加剖面上（圖7a），紅色箭頭指示火山巖反射層，藍色箭頭處可見明顯的火山巖微屈多次波。MWD預測的火山巖多次波模型（圖7b）較SRME（圖7c）更接近真實的多次波反射特征。

MWD的技術原理決定了其難以預測與水層無關的火山巖全程多次波。由于火山巖全程多次波可以由火山巖一次反射褶積得到，不需要海底反射信息，因此，SRME能夠得到較好的多次波模型。根據EP火山巖發育區多次波的特點，本文分別采用MWD和SRME進行多次波預測，然后利用多模型同時自適應相減技術[15]，在最小平方意義下完成MWD和SRME多次波模型與地震數據的最佳匹配，并從地震數據中減去多次波數據。圖8為自適應相減前、后的疊加剖面及差剖面。由圖可見，自適應相減后，不僅海底及正常反射層的多次波得到有效壓制，火山巖的強多次波能量也得到很好的衰減（圖8b）。

2.2火山巖下弱信號恢復

2.2.1火山巖地震波透射損失

當地震縱波從一種介質進入相互密接的另一種介質時，會發生地震波的反射和透射。反射和透射系數可以用Zoeoppritz方程表示

式中：α為縱波入射角或反射角；β為P-SV波反射角；α'為縱波透射角；β'為P-SV波透射角。ρ1、ρ2分別為上覆介質和下伏介質的密度；vP1、vP2分別為上覆介質和下伏介質縱波速度；vS1、vS2分別為上覆介質和下伏介質橫波速度；RPP、RPS分別為反射縱波、反射橫波反射系數；TPP、TPS分別為透射縱波、透射橫波透射系數。

當地震波垂直入射時，即α為零時，地震數據相當于自激自收，因此，α、α'、β、β'都為零，根據式（3）可得

從式（4）可知，上覆地層與下伏地層的波阻抗差異越大，反射波能量越強，而透射波能量越小。由于火山巖和圍巖在速度、密度上存在明顯差異，導致地震波穿過火山巖時透射損失增大，即使前期進行了球面擴散補償等一系列振幅恢復處理，火山巖下伏地層仍然存在振幅補償不足的問題。在多次波等噪聲得到有效壓制后，通過剩余振幅補償進一步增強火山巖下弱信號能量，改善地震數據縱、橫向能量的一致性是提高地震資料品質的重要一步。

2.2.2 AVO剩余振幅補償

在地震勘探中，地震波是非垂直入射的，理論上求取疊前地震數據的透射損失需要求解式（3）。但由于式（3）涉及的地球物理參數較多，地震資料處理中又難以獲得準確的參數信息，因此，透射系數的求取非常困難。剩余振幅補償是基于統計地震數據補償地震波能量損失的一種方法，不需要地層反射系數等先驗信息。在完成球面擴散補償等的基礎上，通過剩余振幅補償可以較好補償透射損失等因素導致的前期振幅恢復不足部分。

AVO剩余振幅補償是在分時窗統計振幅隨炮檢距變化的基礎上補償振幅的一種方式，該方法在地震資料處理中分兩步實施。

第一步為剩余振幅補償函數的計算。首先按等間距劃分成若干個炮檢距組，分別對每個炮檢距組沿時間方向滑動時窗，統計全部輸入CMP的地震道振幅絕對值，計算疊前數據的AVO剩余振幅補償函數。

第n個炮檢距組、第k個時窗的剩余振幅補償函數為

式中：M為輸入的總CMP數；I為總地震道數；J為第k個時窗內的總樣點數；xn為第n個剩余振幅補償函數的參考炮檢距；An（m，i，j）為第m個CMP、第n個炮檢距組、第i道、第k個時窗的第j個樣點的地震道振幅值；P為常數。利用式（5）可以計算不同時窗、不同炮檢距分組的剩余補償函數。

第二步為剩余振幅補償函數的平滑應用。先對計算出的補償函數沿炮檢距組和時間方向進行多點平滑，消除局部異常；然后將補償函數應用于對應的地震道。

由于AVO剩余振幅補償在橫向上按指定的炮檢距分組，縱向上按指定的時窗長度進行振幅統計，最終計算出不同時窗、不同炮檢距分組的補償函數，因此在增加振幅橫向一致性的同時，對透射引起的縱向能量損失有一定的補償作用。從圖9所示剩余振幅補償前（圖9a）、后（圖9b）疊加剖面可以發現，后者火山巖下的反射信號能量得到進一步補償，地震數據能量的縱、橫向一致性得到進一步改善。

2.2.3低頻算子外推隨機噪聲衰減

AVO剩余振幅補償在恢復弱信號振幅的同時，會導致隨機噪聲能量的同步放大。具有低頻算子外推功能的隨機噪聲衰減技術能夠利用低頻優勢頻段預測地震信號，有利于提高火山巖下地震數據的噪聲壓制效果。

本文采用的低頻算子外推隨機噪聲衰減方法將常規RNA（Randon Noise Attenuation）方法擴展到四維空間，即將三維疊前地震數據視為一個擬四維數據體，四個維度分別為線號、CMP號、CMP內的炮檢距、記錄時間，沿時間方向進行傅里葉變換后，利用頻率—空間域預測理論求取每一個頻率成分的預測算子，把預測算子應用于三維疊前地震數據，達到衰減三維疊前隨機噪聲的目的。另外，該方法具有低頻算子外推的功能，其基本技術原理表述為：一組道間時差為Δt1的地震波W1（t），經過傅氏變換后，各道的地震數據可以用一個復數序列表示

w 1（f），w 1（f）e-i2πfΔt1?w 1（f）e-i2πf2Δt1，…（6）

預測算子e-i2πfΔt1既可視作頻率為f、道間時差為Δt1的預測算子，也可以視作頻率為f/2、道間時差為2Δt1的預測算子，同時可以視為頻率為2f、道間時差為Δt1/2的地震數據的預測算子。在地震資料處理中，可以用火山巖高頻散射噪聲范圍之外的低頻成分預測地震數據，以提高隨機噪聲衰減效果。從圖9所示剩余振幅補償前、后及隨機噪聲衰減后的對比剖面可以看出，通過低頻算子外推隨機噪聲衰減，提高了火山巖下地震數據信噪比。

2.3火山巖速度精細刻畫及成像

2.3.1火山巖地震波場特點及建模流程

準確的火山巖速度模型是提高火山巖及下伏構造成像質量的必要前提。EP工區位于水深不足100 m的淺水區，地震資料采集的最小炮檢距為200 m。圖10所示的正演記錄展示了相似地震采集條件下地震數據的波場特征，在最小炮檢距為0時，圖中紅色箭頭所示近道位置可見海底雙曲線反射；當最小炮檢距增大到200 m時，由于出射角超過臨界角，海底一次反射無法被檢波器所接收，僅在綠線以外區域記錄了明顯的回折波信息。從圖11所示實際炮集記錄上同樣可以發現，紅色箭頭所示淺層火山巖及上覆地層的反射波信息相對不足，而在綠線以外的遠炮檢距區域，淺部地層的回折波信息豐富，綠色箭頭所示的初至波位置可見明顯的火山巖高速、強振幅特征。

考慮EP火山巖發育區的地震波場特點，采用如圖12所示的技術流程開展疊前深度偏移速度建模。

技術流程劃分為三個階段。

（1）建立初始PSDM速度模型。用多次波壓制、火山巖弱信號恢復處理后的地震數據開展PSTM速度分析，建立PSTM速度模型并通過約束速度反演將其轉化為深度域層速度模型；利用海底測深文件將海水層速度嵌入PSTM深度域速度模型中，形成如圖13a所示的初始PSDM速度模型。

（2）淺層速度模型更新。拾取單炮初至后，通過菲涅耳帶初至波層析反演，利用中淺層火山巖及上覆地層的走時信息更新淺層速度，形成如圖13b所示的速度模型。

（3）深層速度模型更新。由于初至波層析反演的有效深度受排列長度限制，而且理論上初至波層析反演也難以準確反演高速層下低速地質體的速度，因此，需要利用反射波層析反演進一步優化淺層火山巖以下深部地層的速度。首先通過速度掃描確定深層火山巖相對準確的速度范圍，更新火山巖速度，減小火山巖的速度誤差；然后利用反射波網格層析，精確刻畫深部火山巖及火山巖下伏地層的速度變化，形成如圖13c所示的最終疊前深度偏移速度模型。

2.3.2初至波層析速度建模

和傳統的基于最小走時的初至波層析相比，本文采用的初至波層析技術是一種基于菲涅耳帶初至波走時運算的網格層析技術。

由于地震波是有限帶寬的，地震波主要能量不只沿著最小走時路徑傳播，也在包含最小走時路徑的第一菲涅耳帶內傳播。檢波器接收到的初至波是第一菲涅耳帶內介質的綜合效應。

對于EP地區而言，由于一條菲涅耳帶覆蓋的網格數遠大于一條最小走時射線，菲涅耳帶射線形成的射線密度遠大于最小走時射線。因此，一些由于火山巖等速度異常體引起的射線陰影區，在菲涅耳射線下可能不復存在，有助于提高淺層速度模型的精度。

2.3.3反射波網格層析速度建模

反射波網格層析技術的有效反演深度不受排列長度限制，能夠有效提高深層速度模型的精度。同時，可以自動拾取地震反射傾角和反射同相軸，對反演過程進行約束，減少網格層析速度的多解性[6]。層析方程可簡單表示為

DΔS=Δt（7）

式中：D為給定反演網格的射線追蹤矩陣；ΔS為慢度誤差矩陣；Δt為剩余殘差。實際地震資料處理中首先拾取全局剩余延遲，建立剩余旅行時殘差矩陣，然后通過最小平方算法求解最終剩余速度。

初至波和反射波層析聯合應用的速度建模方案充分利用了火山巖發育區的地震波場信息，提高了速度精度和速度模型準確性。圖14為速度模型切片對比，在最終速度模型上，速度的空間變化與構造的空間展布吻合較好，火山巖的高速異常也被精確刻畫（圖14c）。

圖15是最終速度模型與偏移剖面的疊合顯示，速度模型由淺至深與構造及火山巖分布吻合較好。高精度的速度模型為改善火山巖及下伏構造成像質量奠定了基礎。

2.4應用效果

本文圍繞EP地區火山巖及下伏地層的成像問題，在工區以往疊前深度偏移處理方案的基礎上針對性地優化了成像處理技術。首先通過MWD技術有效提高了火山巖多次波的壓制效果，克服了以往處理方案單純依賴SRME存在的問題；其次，在流程中增加了AVO剩余振幅補償和低頻算子外推隨機噪聲壓制處理，有效補償了火山巖下伏地層反射信號的能量，提高了地震數據信噪比；最后，通過初至波和反射波層析技術的聯合應用，克服了以往單純依賴反射波層析進行速度建模的不足，綜合利用回折波和反射波信息有效提高了火山巖發育區速度模型精度，疊前深度偏移成像質量有了顯著改善。圖16是本次PSDM重處理成果的時間域剖面，和圖1相比，在火山巖反射層以下，控洼斷裂成像更清晰，斷面干脆，復式半地塹結構更容易被識別，洼陷內波組特征清楚，小、微斷裂成像有一定程度改善；受火山巖屏蔽的地層成像得到較好恢復；基于重處理資料可以相對更準確地落實火山巖下構造形態及規模。

圖17是重處理前、后偏移成果的局部放大顯示，在藍色橢圓所示位置，火山巖體的邊界成像有了明顯改善（圓圈所示），有利于準確落實火山巖及正常沉積地層的分布范圍，降低了勘探、開發風險。

3結論

本文采用多次波壓制、信號恢復等一系列技術重新處理了EP地區火山巖發育區的地震資料，得出以下結論。

（1）在淺水火山巖發育區，由于海底反射的缺失，SRME等數據驅動的多次波預測技術難以預測火山巖微屈多次波。MWD是一種以海底測深為先驗信息的模型驅動的多次波預測技術，在淺水火山巖微屈多次波以及其他與水層相關多次波的預測方面更具優勢。

（2）在測井資料不足的情況下，基于振幅統計的AVO剩余振幅補償技術是補償火山巖下伏地層反射能量不足的一種適用方法；在恢復弱振幅的基礎上，應用具有低頻算子外推功能的隨機噪聲衰減技術，可以進一步壓制火山巖下伏地層的隨機噪聲，提高資料信噪比。

（3）淺水火山巖發育區回折波發育，初至波層析可以較準確地刻畫淺層火山巖速度變化；以準確的淺層速度模型為基礎，進一步開展反射波網格層析迭代，能夠明顯提高整體速度模型的準確性和模型精度，為獲得高品質的疊前深度偏移成果奠定基礎。

參考文獻

［1］徐夕生，邱檢生.火成巖巖石學[M].北京：科學出版社，2010.

XU Xisheng，QIU Jiansheng.Igneous Petrology[M].Science Press，Beijing，2010.

［2］趙澄林，朱筱敏.沉積巖石學[M].北京：石油工業出版社，2001.

ZHAO Chenglin，ZHU Xiaomin.Sedimentary Petro?logy[M].Petroleum Iindustry Press，Beijing，2001.

［3］熊翥.復雜地區地震數據處理思路[M].北京：石油工業出版社，2002.

XIONG Zhu.Processing Ideas for Seismic Data in Complex Areas[M].Petroleum Industry Press，Bei?jing，2002.

［4］黎書琴，王曉陽，張孟，等.川西地區二疊系火山巖地震采集關鍵技術應用及效果[J].石油地球物理勘探，2022，57（1）：62?73.

LI Shuqin，WANG Xiaoyang，ZHANG Meng，et al.Application and effect of key seismic acquisition tech?nologies for permian volcanic rocks in western Sichuan basin[J].Oil Geophysical Prospecting，2022，57（1）：62?73.

［5］彭浩天，冉崎，陳康，等.GeoEast系統在四川盆地JY地區二疊系火山巖地震資料處理中的應用[J].石油地球物理勘探，2022，57（增刊1）：6?15.

PENG Haotian，RAN Qi，CHEN Kang，et al.Appli?cation of GeoEast system in seismic data processing of Permian volcanic rocks in JY area of Sichuan Basin[J].Oil Geophysical Prospecting，2022，57（S1）：6?15.

［6］韓站一，尚新民，孫成禹，等.火成巖下伏地層地震反射波場特征研究[J].石油物探，2009，48（3）：226?231.

HAN Zhanyi，SHANG Xinmin，SUN Chengyu，et al.Study on seismic reflection wave field characteris?tics of formations below igneous rock[J].Geophysical Prospecting for Petroleum，2009，48（3）：226?231.

［7］姜巖，李文艷，徐仲達，等.地震傳播的透射損失補償方法及正演模擬分析[J].大慶石油地質與開發，2007（3）：132?134.

JIANG Yan，LI Wenyan，XU Zhongda，et al.Trans?mission dissipation method and forward simulation analysis of seismic propagation[J].Petroleum Geologyamp;Oilfield Development in Daqing，2007（3）：132?134.

［8］張濤，王小衛，田彥燦，等.火成巖速度—深度模型建立方法研究[J].石油地球物理勘探，2015，50（3）：431?435.

ZHANG Tao，WANG Xiaowei，TIAN Yancan，et al.Igneous velocity?depth model building[J].Oil Geo?physical Prospecting，2015，50（3）：431?435.

［9］張濤，王小衛，王孝，等.特殊巖性體的速度—深度建模方法[J].石油地球物理勘探，2020，55（增刊1）：49?55.

ZHANG Tao，WANG Xiaowei，WANG Xiao，et al.The velocity?depth modeling method of special litho?logic bodies[J].Oil Geophysical Prospecting，2020，55（S1）：49?55.

［10］郭樹祥.火山巖發育區地震處理技術及需求分析——以勝利油田濟陽拗陷探區為例[J].石油地球物理勘探，2013，48（2）：228?238.

GUO Shuxiang.Processing technology of seismic data acquired in igneous rock area：example of Jiyang depression in Shengli oilfiled[J].Oil Geophysical Pros?pecting，2013，48（2）：228?238.

［11］王興，曾同生，宋雅瑩，等.滴西火山巖區地震成像關鍵技術[J].石油地球物理勘探，2017，52（增刊1）：47?54.

WANG Xing，ZENG Tongsheng，SONG Yaying，et al.Key techniques for volcanic reservoir imaging in the Dixi area[J].Oil Geophysical Prospecting，2017，52（S1）：47?54.

［12］梁海龍，陳志德.利用高密度地震技術精細刻畫疊置火山巖[J].石油地球物理勘探，2021，56（4）：833?844.

LIANG Hailong，CHEN Zhide.Depiction of superim?posed volcanic edifices in detail by high?density seismic technique[J].Oil Geophysical Prospecting，2021，56（4）：833?844.

［13］弗斯丘爾著;陳浩林，張保慶，等譯.地震多次波去除技術的過去、現在和未來[M].北京：石油工業出版社，2010.

［14］SUN W，WANG H.Model?based water?layer?relateddemultiple with sparse constraints[C].SEG Technical Program Expanded Abstracts，2014，33：4152?4156.

［15］BARNES S R，HEGGE R F.Improved shallow wa?ter demultiple with 3D multi?model subtraction[C].SEG Technical Program Expanded Abstracts，2015，34：4460?4464.

（本文編輯：趙君）

作者簡介

張衛衛高級工程師，碩士，1979年生；2002年獲江漢石油學院地球物理系勘探技術與工程專業學士學位，2005年獲長江大學地球探測與信息技術碩士學位；現就職于中海石油（中國）有限公司深圳分公司，主要從事珠江口盆地油氣勘探、開發相關技術研究及管理工作。