基于X形擴散濾波的逆時偏移低頻噪聲壓制方法

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

Abstract： Thereverse time migration（RTM） imaging method is highly favored for itsability to conduct imaging on complex geological structures without being constrained by the imaging angle. However，this method has its inherent low-frequency noise，thus significantly degradingthe migration imagingquality. Therefore，accurately suppressing the low-frequency noise of RTM is crucial for high-quality migration imaging results. Based on the conventional RTM approaches，this study conducts a detailed analysis of the mechanisms that generate low-frequency noise in RTM. Specifically，the Laplacian filtering method iscurrntly the most commonly adopted suppression method，but it still retains a portion of the low-frequency noise and alters the amplitude information of the seismic data.To this end，alow-frequency noise suppression method based on X-shaped diffusion filtering is proposed.This method is aglobal nonlinear iterative filtering approach that introduces a diffusion function derived from thermodynamics into the RTM noise suppression process，with the RTM data containing low-frequency noise taken as the input.Then，byadjusting the diffusion coeficient and the number of iterations，an X-shaped complementary denoising operator for complex dipping structures is introduced，and a damping function for difusion is set to flexibly control the difusion extent. This allows for the separation of effective signals and low-frequency noise in relatively amplitude-preserving conditions.Numerical simulations on both simple and complex models demonstrate that the X-shaped diffusion filtering method can quickly and ac curately suppress low-frequency noise in RTM and produce high-quality migration results. Keywords： reverse time migration，low-frequency noise，noise suppression，X-shaped diffusion filtering 姜大建，李振春，徐凱，等．基于X形擴散濾波的逆時偏移低頻噪聲壓制方法[J].石油地球物理勘探，2025， 60（4） ：912-922. JIANG Dajian，LI Zhenchun，XU Kai，et al. Low-frequency noise suppression for reverse time migration based on X-shaped diffusion filtering[J]. Oil Geophysical Prospecting，2025，60（4） ：912-922.

0 引言

在地震勘探中，疊前偏移成像方法起著至關重要的作用。它主要包括Kirchhoff積分偏移法1、基于單程波方程的偏移法以及基于雙程波方程的逆時偏移（Reverse TimeMigration，RTM）法[3]。其中，基于雙程波動方程的逆時偏移技術具有諸多優點，如可以正確成像不同傳播類型波場、適應速度橫向變化劇烈的區塊等。隨著計算機技術的快速發展，逆時偏移成像方法的應用從簡單構造轉向復雜構造，從最初的疊后資料逐步擴展到疊前地震資料，并進一步推廣至彈性波領域。研究者們不斷提出新的理論和方法以適應不同的地質條件[4-6]。

成像質量的優化是逆時偏移中的一個重要研究領域。由于逆時偏移中存在著固有的低頻噪聲，而如何準確、有效地壓制低頻噪聲去獲取高質量的偏移剖面，一直以來都是地球物理學者們努力的方向[7-11]。學者們主要從最小二乘偏移成像振幅補償、使用無反射雙程波動方程進行波場構建、角度域道集疊加和濾波方法等四個方面進行了大量研究。曲英銘等[12-13]在起伏海底界面和黏質介質偏移成像中分別引入最小二乘逆時偏移方法，不僅能對起伏海底界面和黏質介質界面進行準確成像，還能壓制成像過程中的低頻噪聲。郭振波等[4將最小二乘逆時偏移方法引人黏聲介質中，不僅對振幅實現了補償，還壓制了低頻噪聲。雖然最小二乘逆時偏移成像方法能壓制偏移成像過程中的低頻噪聲，但是該方法計算量比較大。

還有一些學者通過改進波場延拓方法、引入方向阻尼因子、優化成像條件等手段，不斷提高成像精度并減少對偏移速度場精確估計的依賴。Baysal等[15通過采用無反射波動方程進行波場延拓，該方程能顯著降低邊界處的反射系數，從而提高了成像的清晰度。但是該方法只能應用于疊前資料，而不能應用于疊后資料。隨后，Fletcher等[16在雙程波動方程中引入方向阻尼因子衰減模型內部界面處產生的非必要反射，提升了成像質量。柯璇等17將一步法波場延拓法引人RTM，實現了高精度成像效果。Zhang等[8提出了偽空間域聲波方程，進一步提高了偏移成像質量。

在成像過程中，低頻噪聲往往聚集在大角度成像區域，研究者們提出了角度域濾波、波場分解、直接濾波等方法對其進行壓制。Fomel[19提出了平面波分解法并用于獲取角度域道集，通過限制成像角度來抑制低頻噪聲。Yoon等[20]將Poynting矢量引入波場分解方法。王保利等21提出了提取逆時偏移的角度域共成像點道集方法以壓制低頻噪聲。Yang等22將傳統的Poynting矢量與通過Hilbert變換波場派生的輔助向量相結合，進一步降低了算法的不穩定性，并獲得了更高質量的角度集資料。Zhang等[23對常規Laplace濾波進行了振幅補償，提升了保幅性。Fei等[24提出了基于Hilbert變換的全波場分解法。杜啟振等25發現逆時偏移中低頻噪聲的壓制方法主要分為兩類，一類是逆時偏移在波場重構的過程中減少背向反射；另一類是在成像過程中選擇小角度的波進行成像，并且在成像后對其做濾波處理，即有目標、有選擇地進行成像。陳可洋[26]提出了一種相對保幅的低頻噪聲壓制方法，對偏移噪聲進行了準確的壓制。但是擴散濾波方法在實現的過程中，只在相鄰四個方位進行濾波處理，因而導致擴散濾波在傾斜、復雜構造區域對邊緣信息的保留效果一般。

為此，本文從常規逆時偏移出發，對逆時偏移的低頻噪聲產生機制進行詳細分析。然后，介紹三類常用的低頻噪聲壓制方法，并在此基礎上提出針對復雜傾斜構造的基于X形擴散濾波的低頻噪聲壓制方法。該方法是一種全局非線性迭代濾波方法，將熱力學領域的擴散函數引人逆時偏移噪聲壓制中，通過調節控制擴散系數和迭代次數，實現低頻噪聲和有效信號的分離。最后，通過對簡單模型和復雜模型進行試算以驗證方法效果。將基于X形擴散濾波的方法引入逆時偏移中，能快速準確地壓制逆時偏移中的低頻噪聲，獲取高質量的偏移成像結果，可為后續的地震資料處理和解釋提供幫助。

1方法原理

1.1偏移噪聲產生機理

采用雙程波動方程[27-31]進行逆時偏移成像，相比單程波更能改善復雜區域的成像效果，但往往會帶來由背向散射項引起的低頻噪聲[32-33]。低頻噪聲是逆時偏移中最主要、最常見的問題之一，它會降低成像剖面質量，對后續成像資料的處理和解釋造成干擾。

根據時間一致性原理，成像點位于同時刻震源正向延拓波場與接收逆向延拓波場上[34]。基于雙程波偏移的互相關成像條件表示為

式中： u_s，u_r 分別表示震源波場和檢波器接收波場；t為時間。

為解釋低頻噪聲的物理意義，將震源波場和接

收波場分別根據成像點處波傳播的方向分為上行波場和下行波場，用得到的四個分量將基于全波場的互相關成像條件分為四個分量，即

式中： u_s^u，u_s^d 分別表示震源波場的上行波場和下行波場； u_r^u，u_r^d 分別表示接收波場的上行波場和下行波場。

成像公式可改寫為

I₁，I₂ 的成像分別如圖1a和圖1b所示。由圖可見，震源波場的下行波場與檢波點波場的上行波場進行互相關運算后，其結果類似于單程波偏移的成像條件，主要呈現為一次反射波的成像效果；震源波場的上行波場與檢波點波場的下行波場的互相關結果呈現出能量較弱的一次波成像結果。這兩個互相關項都在真實反射點上成像，因此稱 I₁，I₂ 為前向散射項。

圖1前向散射項 （a，b） 和背向散射項 （c，d） 波傳播示意圖

如圖1c、圖1d所示 ?，I₃ 與 I₄ 分別代表震源波場的上行波場與檢波點波場的上行波場之間、震源波場的下行波場與檢波點波場的下行波場之間的互相關，這些互相關項在偏移剖面中造成假成像點，因此稱為背向散射項。

由圖1可見，此時反射角相對較大， I₂ 在頻率域表現出低頻特征。但當偏移速度變化較劇烈或反射界面很復雜時，回轉波的出現會使成像結果出現假像，如圖2所示。此時， ?_J3 在成像剖面中表現為假成像點，且其能量可能強于多次波的成像結果。

圖2回轉波假成像點

由上述可得，用雙程波動方程模擬的震源波場在遇到波阻抗界面時，由于正傳波場和反傳波場都會產生背向散射，成像結果中不僅包含真實成像點，還包含了假成像點和低頻噪聲。此外，由于地下淺層部分背向散射能量較強，低頻噪聲主要存在于深度較淺的區域。

1.2偏移噪聲壓制方法

根據對假成像的不同壓制策略，低頻噪聲的壓制方法主要分為三類：修改波動方程、優化成像條件和成像后直接進行濾波處理[35]。

1.2.1 修改波動方程

修改波動方程的噪聲壓制法的核心是通過對波動方程進行改進，使低頻噪聲能量在傳播的過程被削弱。Baysal等3通過恒定波阻抗，將入射時的界面反射完全消除；Fletcher等5在選定地下界面后，通過在反射界面上引入阻尼系數消除背向散射。這一類方法在實際操作中流程簡單，實現方便，但由于其同時消除了回轉波、棱柱波等具有多個傳播方向的波，使雙程波逆時偏移的優勢喪失，因此不適用于高陡、傾角、巖丘等結構。

1.2.2優化成像條件

優化成像條件法的核心思想是對成像條件進行改造，開展選擇性成像，使最終的偏移結果中只包含有效波能量。常使用的修改條件主要有波場分解類和角度域衰減類。波場分解類方法是通過分離震源波場和檢波場的上、下行波場，根據前述噪聲產生機理，只對形成真實反射面的散射項進行互相關。角度域衰減類方法是通過提取角度域共成像點道集，只保留其中較小角度部分進行成像，以對低頻噪聲進行壓制，該類方法在實際操作中通常耗時長。

1.2.3 成像后濾波

成像后濾波法是一種直接針對成像結果進行濾波處理的簡便方法。這種方法無需復雜的預處理步驟，因此操作起來相對簡單、直接。其中常用的濾波方法包含拉普拉斯濾波法、導數濾波法、高通濾波法以及最小二乘濾波法等[36]。每種方法都有其獨特的應用場景和優勢。通過這些濾波方法，可以有效地改善成像質量，提高成像的清晰度和準確性。其中，最常用的是拉普拉斯濾波法。

1.3基于X形擴散濾波低頻噪聲壓制方法

成像后濾波法直接對成像結果進行低頻噪聲濾除，不需要在波場延拓和成像時做處理，操作起來簡單、快捷。在地震資料成像領域，拉普拉斯濾波法是目前最常用的濾波方法，但它依舊保留一部分低頻噪聲，且改變地震數據的振幅信息。針對這些問題，本文提出了X形擴散濾波法。下面將介紹這兩種濾波方法。

1.3.1 拉普拉斯濾波去噪原理

為直觀介紹拉普拉斯濾波法對地震偏移剖面的噪聲壓制原理，本文在波數域對拉普拉斯算子的頻率域表達式進行推導。Zhang等23提出了針對逆時偏移結果中低頻噪聲的拉普拉斯濾波法。在二維空間中，拉普拉斯算子表示為

將式（4）進行傅里葉變換，得到波數域的表達式

FFT（?²）-（k_x²+k_z²）=-∣k₁∣²

式中：FFT為快速傅里葉正變換； k_x，k_z 分別表示水平、垂直方向上的波數； k_l 為成像域波數矢量。

k₁ 可以表示為

k_l=k_R-k_S

式中： k_s 表示震源波場的波數矢量； k_R 表示檢波點波場的波數矢量。

k_l，k_s，k_R 的矢量關系如圖3所示，根據余弦定理可得

式中： θ 為反射角； v 為速度； ω=2πf， 其中 f 為頻率。

將式（7）代入式（5），得到拉普拉斯算子的頻率域表達式

由此可見，拉普拉斯濾波對成像結果的作用實際上等同于對逆時偏移的角道集進行了 cos²θ 的衰減處理，并同時給振幅帶來了 -4ω²/γ² 的影響。

圖3成像波數域矢量計算示意圖

1.3.2 X形擴散濾波去噪原理

擴散濾波法是一種全局非線性迭代濾波方法，主要應用于熱力學領域，它可以在相對保幅的條件下將有效信號與低瀕噪聲分開。各向同性非線性擴散濾波方程的微分形式記為

u_i，j-1^t-4u_i，j^t）

式中 i，j 分別為水平、垂直方向的網格點。

為了簡單起見，用 N 代替 t 來表示迭代次數，則

Δu_i，j-1^N-4u_i，j^N）

在實際程序實現的過程中，由于只在相鄰四個方位進行濾波處理，導致擴散濾波在傾斜、復雜構造區域對邊緣信息的保留效果一般。為此，Zhang等[引入針對復雜傾斜構造的X形補充去噪算子，該算子的 3×3 矩陣表示為

引入X形補充去噪算子后的X形擴散濾波的差分公式表示為

由于X形擴散濾波的結果受擴散系數 α 、迭代次數 N 的影響，因此在面對更加復雜多樣的地質構造時，可以靈活地選取合適的參數以獲得高質量的成像結果，這相當于人為設置擴散阻尼函數來調控擴散的程度。

2模型試算

式中： α 表示擴散系數（一般來說 0lt;αlt;0.3 ，擴散系數越大，平滑程度越大）； t 表示擴散時間，在計算過程中用迭代次數 N 表示； u₀ 表示初始的偏移圖像； u 表示被濾除的低頻噪聲。可用 u₀-u 表示逆時偏移結果中的有效信號。

空間階數選為三階時，用有限差分法對擴散濾波方程（式（9））進行數值離散化，得到

為了對比拉普拉斯濾波法和X形擴散濾波法對逆時偏移低頻噪聲的壓制效果，首先設計一個簡單二層平層模型進行測試。模型的上層速度為2000m/s ，下層速度為 3000m， s。設置20炮震源，均勻分布在地表，炮間距為 100m 。同時，在地表間隔 5m 均勻設置400個檢波器。時間采樣間隔為0.8ms ，采樣時間長度為 1.6s 。

原始逆時偏移成像結果如圖4a所示，由于照明不均及背向散射，地層淺層部位存在能量較強的低頻噪聲。分別對逆時偏移結果做拉普拉斯濾波、擴散濾波 （α=0.05，N=50） X形擴散濾波 （α=0，05 ，N=50. 處理，效果分別如圖 4b～ 圖4d所示。由圖可見，三種方法都可以壓制低頻噪聲，但與拉普拉斯濾波、擴散濾波方法相比，X形擴散濾波方法對逆時偏移成像結果中的低頻噪聲的壓制更加徹底。

圖4常規逆時偏移及不同濾波方法后的偏移成像效果對比

為驗證X形擴散濾波法對復雜模型的適用性，選用鹽丘模型的疊后逆時偏移成像結果，分別對其進行拉普拉斯濾波處理和X形擴散濾波處理。鹽丘速度模型如圖5所示，深度為 1.5km ，水平距離為 6.5km ，兩個方向的網格間距均為 10m 。采用主頻為 30Hz 的雷克子波作為震源，設置100個炮點，均勻分布于地表，炮間距為 65m ，每隔 10m 均勻設置650個檢波器。時間采樣間隔為 0.8ms ，時間采樣長度為 4.0s 。該模型的逆時偏移結果以及兩種濾波法處理后的成像剖面如圖6所示。

由于鹽丘模型含有強反差速度結構，過強的一次反射波能量在成像過程中會掩蓋掉其他波形的能量。由圖6a可見，在常規逆時偏移結果中，水平構造的成像效果相對較好，但鹽丘頂部相對復雜的構造幾乎難以辨認。將圖6b、圖6c選取一部分區域（紅框）進行局部放大顯示（圖7）。觀察鹽丘頂部構造可以看出，X形擴散濾波法相較于拉普拉斯濾波法，在相對復雜的地質構造處，能在保證成像軸連續性的同時消除更多成像假象（如圖7中紅色箭頭所示）。

圖5鹽丘速度模型

圖6鹽丘模型不同方法逆時偏移成像效果對比

由于擴散系數 α 是一個非常關鍵的參數，它決定了擴散的速率，也就決定了對成像結果的平滑程度。為觀察不同擴散系數對濾波效果的影響，對Marmousi模型（圖8）疊后逆時偏移結果進行測試。分別進行拉普拉斯濾波和擴散系數為0.001、0.01、0.1、0.25的X形擴散濾波。為便于觀察，截取偏移結果濾波處理后的局部（圖8紅框）放大顯示（圖9）進行對比。

圖7圖6b、圖6c紅框區域放大顯示

圖8Marmousi速度模型和逆時偏移成像

圖9原始逆時偏移（a）及其不同方法處理的逆時偏移（b～f）剖面效果對比

由圖9b、圖9e、圖9f可見，拉普拉斯濾波法和擴散系數較大的X形擴散濾波法都會保留一部分由于強反射界面造成的低噪聲。由圖 9c～ 圖9f可知，對于X形擴散濾波，擴散系數 α 越大，成像結果的平滑程度也越強。

通過改變擴散系數的大小，可以平衡平滑與邊緣保護。如果擴散系數過大，則會導致圖像模糊；如果擴散系數過小，則會保留過多的噪聲和細節，所以擴散系數的選取至關重要。經試驗，對于常規的成像結果，擴散系數選取 0.01～0.05 時，可以實現在保留有效信息的同時盡可能地壓制噪聲。

由圖9b、圖9d可見，相較于拉普拉斯濾波法，選擇合適擴散系數 （α=0.01） 的X形擴散濾波處理后的偏移結果中，同向軸更加清晰、連續性更好，這在一定程度上提升了資料的信噪比。

3 實際資料測試

利用中國東部M探區實際數據測試本文方法。M區長 x 寬為 1.3km×3.5km ，二維線總炮數共550炮，炮檢距為 4km ，道間距為 12.5m ，震源主頻為 30Hz 。

圖10a為逆時偏移的成像結果，從淺層可以明顯看出低頻噪聲。經過拉普拉斯濾波后的逆時偏移成像（圖10b），有效能量有明顯的損失，特別是紅色橢圓區域。而本文提出的X形擴散濾波能有效壓制低頻噪聲，且更加有效地保證了成像結果的質量（圖10c）。

圖10M區常規逆時偏移及不同濾波方法后逆時偏移的成像結果對比

4結束語

本文主要闡述了逆時偏移的基本原理、成像條件、低頻噪聲的產生機理以及三類針對低頻噪聲的有效壓制方法。首先，探討了逆時偏移中常用的成像條件，特別是互相關成像條件，因其實現簡便且效率高，目前在實踐中廣泛應用。然后，對RTM算法中低頻噪聲的生成機理以及當前主流的噪聲壓制方法進行了概述并總結。詳細介紹了兩種低頻噪聲壓制方法，即拉普拉斯濾波方法和X形擴散濾波法。通過對簡單模型和復雜模型的數值模擬，驗證了又形擴散濾波法在低頻噪聲壓制方面的有效性和優勢。并且，對擴散系數的選取進行了測試，給出了擴散系數選取的參考范圍，驗證了所提方法的有效性。該方法可為地震資料的后續處理與解釋提供幫助。

參考文獻

[1]SCHNEIDER W A. Integral formulation for migra tionin two and three dimensions[J]. Geophysics， 1978，43（1）：49-76.

[2] CLAERBOUTJF.Toward aunified theoryof reflec tor mapping[J].Geophysics，1971，36（3）：467-481.

[3] BAYSALE，KOSLOFFDD，SHERWOODJW C.Reverse time migration[J]. Geophysics，1983， 48（11）： 1514-1524.

[4] ZHAN G，PESTANA R C，STOFFA PL. An acoustic wave equation for pure P wave in 2D TTI media[C]//SEG Technical Program Expanded Abstracts， 2011，30：168-173.

[5] FLETCHERRP，NICHOLSD，CAVALCA M. Wavepath-consistent effective Q estimation forQ-com pensated reverse-time migration[C]//Extended Abstracts of 74th EAGE Conference and Exhibition， 2012：A020-A025.

[6] SHRAGGE J. Reverse time migration from topogra phy[J].Geophysics，2014，79（4）： S141-S152.

[7] FORNBERG B. Classroom note： Calculationof weights in finite difference formulas[J]. SIAM Review，1998，40（3）：685-691.

[8] LIUY，SEN M K.Time-space domain dispersionrelation-based finite-difference method with arbitrary even-order accuracy for the 2D acoustic wave equation [J].Journal of Computational Physics，2Ol3，232（1）： 327-345.

[9] ALTERMAN Z，KARAL F C. Propagation of elasticwavesinlayered mediabyfinitedifferencemethods [J].Bulletin of the Seismological Society of America， 1968，58（1）：367-398.

[10] ZHANGJH，YAOZX.Optimized finite-difference operator for broadband seismic wave modeling[J]. Geophysics，2013，78（1）：A13-A18.

［I」胡目多，劉威，宋冢文，等．速度一應力聲波萬程混 合交錯網格有限差分數值模擬及逆時偏移[J].石油 地球物理勘探，2023，58（5）：1084-1100. HU Ziduo，LIU Wei，SONG Jiawen，et al.Mixed staggered grid finite difference numerical simulation and reverse time migration of velocity-stress acoustic equation[J]. OilGeophysicalProspecting，2023，58（5）： 1084-1100.

[12］曲英銘，黃崇鵬，劉暢，等．基于海底電纜數據的聲 彈耦合介質多參數最小二乘逆時偏移[J].應用地球 物理，2019，16（3）：327-337，396. QUYingming，HUANG Chongpeng，LIU Chang，et al. Multiparameter least-squares reverse time migration for acoustic-elastic coupling based on ocean bottom cable data[J].Applied Geophysics，20l9，16（3）： 327-337，396.

[13］曲英銘，李金麗，王云超，等．最小平方逆時偏移中 黏彈性和各向異性的校正：以渤海灣地區地震數據為 例[J].地球物理學報，2019，62（6）：2203-2216. QU Yingming，LI Jinli，WANG Yunchao，et al.Cor rection of viscoelasticity and anisotropy in least-squares reverse time migration： A Bohai Bay seismic case study [J].Chinese Journal of Geophysics，2019，62（6）： 2203-2216.

[14］郭振波，李振春．最小平方逆時偏移真振幅成像[J]. 石油地球物理勘探，2014，49（1）：113-120. GUO Zhenbo，LI Zhenchun. True-amplitude imaging based on least-squares reverse time migration[J]. Oil Geophysical Prospecting，2014，49（1）：113-120.

[15] BAYSAL E，KOSLOFF D D，SHERWOOD JW C.A two-way nonreflecting wave equation[J]. Geophysics，1984，49（2）：132-141.

[16] FLETCHERRP，FOWLERPJ，KITCHENSIDE P，et al. Suppressing unwanted internal reflections in prestackreverse-timemigration[J]. Geophysics， 2006，71（6）：E79-E82.

[17］柯璇，石穎．基于一步法波場延拓的正演模擬和逆 時偏移成像[J].地球物理學報，2017，60（11）：4468- 4479. KE Xuan，SHI Ying. Forward simulation and reverse time migration imaging based on one step wavefield extrapolation[J]. Chinese Journal of Geophysics，2017， 60（11） ： 4468-4479.

[18]ZHANG X，WANG X，LIU B，et al. Reverse time migration based on the pseudo-space-domain first-ordervelocity-stress acoustic wave equation[J]. Frontiers in Earth Science，2021，9： 690513.

[19]FOMEL S. Applications of plane-wave destruction filters[J].Geophysics，2002，67（6）： 1946-1960.

[20]YOON K，MARFURT K J. Reverse-time migration using the poynting vector[J]. Exploration Geophysics， 2006，37：102-107.

[21］王保利，高靜懷，陳文超，等．逆時偏移中用 Poynting矢量高效地提取角道集[J].地球物理學報， 2013，56（1）：262-268. WANG Baoli，GAO Jinghuai，CHEN Wenchao，et al.Extracting efficiently angle gathers using Poynting vector during reverse time migration[J]. Chinese Journal of Geophysics，2013，56（1）：262-268.

[22]YANG K，DONG X，WANG X，et al. Deprimary reverse time migration angle gathers with a stabilized Poynting vector[J].Geophysics，2022，87（5）：S267- S279.

[23] ZHANG Y，SUN J. Practical issues of reverse time migration： True amplitude gathers noise removal and harmonic-source encoding[C]// Beijing 2Oo9 International Geophysical Conferenceand Exposition， 2009：204.

[24] FEITW，LUOY，YANGJ，etal.Removing 1 images in reverse time migration： The concept of deprimary[J].Geophysics，2015，80（6）： S237-S244.

[25] 杜啟振，朱釔同，張明強，等．疊前逆時深度偏移低 頻噪聲壓制策略研究[J].地球物理學報，2013，56（7）： 2391-2401. DU Qizhen，ZHU Yitong，ZHANG Mingqiang，et al.A study on the strategy of low wavenumber noise suppression for prestack reverse time depth migration [J].Chinese Journal of Geophysics，2013，56（7）： 2391-2401.

[26] 陳可洋．一種相對保幅的低頻逆時噪聲壓制方法及 其應用[J].油氣藏評價與開發，2014，4（3）：50-54. CHEN Keyang. A relative amplitude preservation lowfrequency reverse time noise suppression scheme and its application[J]. Reservoir Evaluation and Development，2014，4（3）：50-54.

[27]LIU Y，SEN M K. A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation[J].Geophysics，2010，75（2）： A1-A6.

[28] 劉洋，李承楚，牟永光．任意偶數階精度有限差分法 數值模擬[J]．石油地球物理勘探，1998，33（1）：1- 10. LIU Yang，LI Chengchu，MOU Yongguang. Finite difference numerical modeling of any even-order accuracy[J].Oil Geophysical Prospecting，1998，33（1） ： 1- 10.

[29]LIU Y，SEN MK. Finite-difference modeling with adaptive variable -length spatial operators[J]. Geophysics，2011，76（4）：T79-T89.

[30] 董良國，馬在田，曹景忠，等．一階彈性波方程交錯 網格高階差分解法[J].地球物理學報，200O，43（3）： 411-419. DONG Liangguo，MA Zaitian，CAO Jingzhong，et al.A staggered-grid high-order difference method of one-order elastic wave equation[J]. Chinese Journal of Geophysics，2000，43（3） ： 411-419.

[31］馬雄，申天賜，田永凈，等．頻率域穩定化的黏聲波 逆時偏移方法及其應用[J].石油地球物理勘探， 2024，59（6）：1315-1329. MA Xiong，SHEN Tianci，TIAN Yongjing，et al. Stabilized viscoacoustic reverse time migration in the frequency domain and its application[J]. Oil Geophysi calProspecting，2024，59（6）：1315-1329.

[32] YOONK，MARFURTKJ，STARRW.Challenges inreverse time migration[C]//SEG Technical Pro gramExpanded Abstracts，2004，23：1057-1060.

[33] 張浩，馮興強，付昌，等．基于卷積神經網絡的傾角 域彈性波逆時偏移噪聲壓制方法[J].石油物探， 2021，60（3）：376-384. ZHANGHao，FENG Xingqiang，FUChang，etal. Noise suppression during elastic reverse time migration in the dip-angle domain using a convolution neural network[J]. Geophysical Prospecting for Petroleum， 2021，60（3）：376-384.

[34] 王文化．高精度有限差分法波動方程正演及逆時偏 移研究[D].成都：成都理工大學，2021. WANG Wenhua.Research on High-precision Finite Difference Method for Wave Equation Forward Modeling and Reverse Time Migration[D]. Chengdu： Chengdu University of Technology，2021.

[35] 黃杰．TI介質擬聲波方程正演模擬與逆時偏移[D]. 青島：中國石油大學（華東），2018. HUANGJie.Forward ModelingandReverse Time Migration of qP-wave Equation in TI Media[D]. Qingdao：China University of petroleum（East China）， 2018.

[36] 吳天麒．逆時偏移角道集優化處理方法[D].大慶：東 北石油大學，2021. WU Tianlin. Optimization Method of Angle-domain gather Based on Reverse-time Migration[D]. Daqing： Northeast PetroleumUniversity，2021.

[37] ZHANGL，REN Y，CHENL，etal.A newdenoisingstrategy and the X-shaped supplement denoising operator targetingtime-reversed mirror imaging technique[J]. Journal of GeophysicsandEngineering， 2021，18（4）：503-514.

（本文編輯：謝結來）

作者簡介

姜大建高級工程師，1979年生；2002年獲中國石油大學（北京）地質工程專業學士學位；現就職于中國石化石油勘探開發研究院地球物理技術研究所，主要從事地震信號處理、成像等方面的研究。