3D Reverse Time Migration Research Articles

Million-core scalable 3D anisotropic reverse time migration on the Sugon exascale supercomputer

Reverse time migration (RTM) plays a crucial role in high-resolution seismic imaging of the Earth’s interior. However, scaling it across millions of cores in parallel to process large-scale seismic datasets poses significant computational challenges, because the conventional storage solutions are insufficient to deal with the I/O and memory bottlenecks. To address this issue, we present a highly scalable 3D RTM algorithm for vertically transverse isotropic (VTI) media, optimized for the Sugon exascale supercomputer, utilizing over 1,024,000 cores with optimal weak-scaling efficiency. Through cache optimizations tailored for the new DCU accelerator architecture, our approach achieves a 4-6 × speedup compared to conventional methods on a single accelerator. Moreover, based on the lossy compression and boundary-saving techniques, we reduce storage requirements by 266 times, which allows for the effective utilization of million-core computing resources and ensures scalability efficiency when handling large-scale datasets for complex geophysical tasks. Finally, when applied to a industrial dataset, the method demonstrates robust scalability and high efficiency, making it well-suited for large-scale geophysical exploration.

Compensating Acquisition Footprint for Amplitude-Preserving Angle Domain Common Image Gathers Based on 3D Reverse Time Migration

Angle domain common image gathers (ADCIGs) play a crucial role in seismic exploration, offering prestack underground illumination information that aids in validating migration velocity and conducting prestack amplitude versus angle (AVA) analysis for reservoir characterization. This paper introduces an innovative approach for compensating amplitude errors caused by irregular seismic acquisition geometries in ADCIGs. By incorporating an angle domain illumination compensation factor, the proposed method effectively modifies these errors, preserving the amplitude of seismic reflectivity in the prestack angle domain. The effectiveness of the proposed approach is validated through comprehensive tests conducted on synthetic and field data examples. The results demonstrate the capability of the method to enhance the quality of ADCIGs derived from 3D reverse time migration (RTM), yielding accurate and reliable amplitude preservation. While the illumination compensation factor assumes a vertically linear velocity model, the method holds promise for extension to more complex media and diverse migration techniques. This suggests its applicability and adaptability beyond the specific assumptions considered in this study. In conclusion, this paper presents an innovative angle domain illumination compensation factor that significantly improves the quality of ADCIGs by addressing amplitude errors arising from irregular seismic acquisition geometries. The experimental validation using synthetic and field data confirms the effectiveness of the proposed method within the context of 3D RTM. Furthermore, the technique holds potential for broader application in more complex subsurface scenarios and various migration methodologies.

3D reverse-time migration for pure P-wave in orthorhombic media

Compared with the transverse isotropic (TI) medium, the orthorhombic anisotropic medium has both horizontal and vertical symmetry axes and it can be approximated as a set of vertical fissures developed in a group of horizontal strata. Although the full-elastic orthorhombic anisotropic wave equation can accurately simulate seismic wave propagation in the underground media, a huge computational cost is required in seismic modeling, migration, and inversion. The conventional coupled pseudo-acoustic wave equations based on acoustic approximation can be used to significantly reduce the cost of calculation. However, these equations usually suffer from unwanted shear wave artifacts during wave propagation, and the presence of these artifacts can significantly degrade the imaging quality. To solve these problems, we derived a new pure P-wave equation for orthorhombic media that eliminates shear wave artifacts while compromising computational efficiency and accuracy. In addition, the derived equation involves pseudo-differential operators and it must be solved by 3D FFT algorithms. In order to reduce the number of 3D FFT, we utilized the finite difference and pseudo-spectral methods to conduct 3D forward modeling. Furthermore, we simplified the equation by using elliptic approximation and implemented 3D reverse-time migration (RTM). Forward modeling tests on several homogeneous and heterogeneous models confirm that the accuracy of the new equation is better than that of conventional methods. 3D RTM imaging tests on three-layer and SEG/EAGE 3D salt models confirm that the ORT media have better imaging quality.

3D Reverse-Time Migration Imaging for Multiple Cross-Hole Research and Multiple Sensor Settings of Cross-Hole Seismic Exploration

The two-dimensional (2D) cross-hole seismic computed tomography (CT) imaging acquisition method has the potential to characterize the target zone optimally compared to surface seismic surveys. It has wide applications in oil and gas exploration, engineering geology, etc. Limited to 2D hole velocity profiling, this method cannot acquire three-dimensional (3D) information on lateral geological structures outside the profile. Additionally, the sensor data received by cross-hole seismic exploration constitute responses from geological bodies in 3D space and are potentially affected by objects outside the well profiles, distorting the imaging results and geological interpretation. This paper proposes a 3D cross-hole acoustic wave reverse-time migration imaging method to capture 3D cross-hole geological structures using sensor settings in multi-cross-hole seismic research. Based on the analysis of resulting 3D cross-hole images under varying sensor settings, optimizing the observation system can aid in the cost-efficient obtainment of the 3D underground structure distribution. To verify this method’s effectiveness on 3D cross-hole structure imaging, numerical simulations were conducted on four typical geological models regarding layers, local high-velocity zones, large dip angles, and faults. The results verify the model’s superiority in providing more reliable and accurate 3D geological information for cross-hole seismic exploration, presenting a theoretical basis for processing and interpreting cross-hole data.

Research on the 3D Reverse Time Migration Technique for Internal Defects Imaging and Sensor Settings of Pressure Pipelines.

Although pressure pipelines serve as a secure and energy-efficient means of transporting oil, gas, and chemicals, they are susceptible to fatigue cracks over extended periods of cyclic loading due to the challenging operational conditions. Their quality and efficiency directly affect the safe operation of the project. Therefore, a thorough and precise characterization approach towards pressure pipelines can proactively mitigate safety risks and yield substantial economic and societal benefits. At present, the current mainstream 2D ultrasound imaging technology faces challenges in fully visualizing the internal defects and topography of pressure pipelines. Reverse time migration (RTM), widely employed in geophysical exploration, has the capability to visualize intricate geological structures. In this paper, we introduced the RTM into the realm of ultrasonic non-destructive testing, and proposed a 3D ultrasonic RTM imaging method for internal defects and sensor settings of pressure pipelines. To accurately simulate the extrapolation of wave field in 3D pressure pipelines, we set the absorbing boundary and double free boundary in cylindrical coordinates. Subsequently, using the 3D ultrasonic RTM approach, we attained higher-precision 3D imaging of internal defects in the pressure pipelines through suppressing imaging artifacts. By comparing and analyzing the imaging results of different sensor settings, the design of the observation system is optimized to provide a basis for the imaging and interpretation of actual data. Both simulations and actual field data demonstrate that our approach delivers top-notch 3D imaging of pipeline defects (with an imaging range accuracy up to 97.85%). This method takes into consideration the complexities of multiple scattering and mode conversions occurring at the base of the defects as well as the optimal sensor settings.

Three-Dimensional Ultrasonic Reverse-Time Migration Imaging of Submarine Pipeline Nondestructive Testing in Cylindrical Coordinates

Submarine pipelines are a safe and energy-efficient mode of gas transport. However, due to the complex manufacturing process and harsh operating environment, submarine pipelines are subject to fatigue cracks under long-term cyclic loading. A comprehensive and high-precision characterization strategy for submarine pipelines can effectively prevent potential safety hazards and have significant economic and social repercussions. As a matter of fact, pipeline defects cannot be reliably detected with current traditional 2D methods. On the other hand, in ultrasonic testing, cylindrical geometry increases the complexity of the 3D wave field in the submarine pipeline space and significantly influences the accuracy of the detection results. In this paper, we put forward a novel method for 3D ultrasonic image testing that is suitable for cylindrical coordinates. In order to accurately simulate the ultrasonic signal received from pipelines, we generalize the 3D staggered-grid finite-difference method from Cartesian coordinates to cylindrical ones and simulate the full wave field in the 3D pipeline space. Then, signal processing is performed on the ultrasound simulation records, and 3D reverse-time migration imaging of submarine pipeline defects can be effectively achieved using the reverse-time migration method and cross-correlation imaging conditions. The results obtained from simulations and real field data show that the proposed method provides high-quality 3D imaging of defects in pipelines, taking into account multiple scattering and mode conversion information at the bottom of the defects.

Atlantis — 20 years of seismic innovation finally removes the shroud of mystery

The Atlantis Field has gone through more than two decades of continuous seismic imaging efforts, during which time many innovative technologies were incubated, the most recent one being the successful application of full-waveform inversion (FWI) in salt environments. This technique led to a significant improvement in the subsalt image. However, imaging challenges remain for the Atlantis reservoirs, primarily due to the complex overburden salt geometries and the highly compartmentalized reservoir. Even with an improved velocity model from FWI, the conventional reverse time migration (RTM) images still suffer from illumination issues and contain strong migration swings that hinder the subsalt imaging and subsequent interpretations. Furthermore, early versions of FWI employed an acoustic assumption, leading to visible salt halos at the salt boundaries in the velocity model, which adversely impacted the reservoir imaging. In the last 12 months, elastic time-lag FWI (TLFWI) and FWI-derived reflectivity (FDR) imaging using long-offset ocean-bottom node data have minimized these imaging issues at Atlantis, providing another step change in subsalt understanding. Although the 3D RTM images using the elastic FWI velocity model are similar overall to their acoustic counterparts, the 4D time-lapse RTM images at Atlantis show noticeable improvements. Furthermore, FDR images derived from elastic FWI velocities show obvious benefits over the acoustic ones. With a more accurate modeling engine that allows for better match between synthetic and real data, FDR imaging shows improved illumination, higher signal-to-noise ratio, and better reservoir details over acoustic FDR imaging. This recent advancement in using elastic TLFWI has had immediate positive effects in facilitating the Atlantis Field's current and future development.

High-Precision Depth Domain Migration Method in Imaging of 3D Seismic Data in Coalfield

Fault structures developed in coal seams, which are often associated with roof collapse, water inrush, gas outburst, and other accidents, are common geological hazards in coal exploration and development. The accurate detection of micro-structures such as small faults has always been a research focus to ensure safety in coalfields. Three-dimensional (3D) seismic research is one of the most efficient methods for obtaining the structural characteristics of coal areas and identify small faults in coal seams, but it is difficult for traditional seismic data imaging technologies to meet the high-precision demand of current coal exploration. Aiming at the characteristics of 3D seismic data in coalfields, we calculated the difference coefficients based on the optimized inversion algorithm and proposed a variable-density acoustic equation optimized with a temporal–spatial staggered-grid finite difference forward algorithm. On this basis, by combining normalized cross-correlation imaging conditions and GPU/CPU collaborative parallel processing technology, we developed an efficient and high-precision 3D reverse time migration method suitable for 3D seismic data in coalfields. Numerical tests verified the accuracy and efficiency of the proposed migration method for the imaging of coal-measure strata with small fault structures and could effectively identify 5 m small faults in the coal seam. The migration test of 3D seismic data in real coalfields showed that our 3D reverse time migration method has good practicability for high-precision imaging of 3D seismic data in coalfields and is an effective method for the precise imaging of small faults in coal measures.

Memory-efficient source wavefield reconstruction and its application to 3D reverse time migration

Reverse time migration (RTM) generally uses the zero-lag crosscorrelation imaging condition, requiring the source and receiver wavefields to be known at the same time step. However, the receiver wavefield is calculated in time-reversed order, opposite to the order of the forward-propagated source wavefield. The inconvenience can be resolved by storing the source wavefield on a computer memory/disk or by reconstructing the source wavefield on the fly for multiplication with the receiver wavefield. The storage requirements for the former approach can be very large. Hence, we have followed the latter route and developed an efficient source wavefield reconstruction method. During forward propagation, the boundary wavefields at N layers of the spatial grid points and a linear combination of wavefields at [Formula: see text] layers of the spatial grid points are stored. During backward propagation, it reconstructs the source wavefield using the saved wavefields based on a new finite-difference stencil ( M is the operator length parameter, and [Formula: see text]). Unlike existing methods, our method allows a trade-off between accuracy and storage by adjusting N. A maximum-norm-based objective function is constructed to optimize the reconstruction coefficients based on the minimax approximation using the Remez exchange algorithm. Dispersion and stability analyses reveal that our method is more accurate and marginally less stable than the method that requires storage of a combination of boundary wavefields. Our method has been applied to 3D RTM on synthetic and field data. Numerical examples indicate that our method with [Formula: see text] can produce images that are close to those obtained using a conventional method of storing M layers of boundary wavefields. The memory usage of our method is [Formula: see text] times that of the conventional method.

3D reverse time migration using a wavefield domain dynamic approach

3D Reverse-time migration (RTM) is a powerful technique for imaging complex geologic structures. This approach requires a significant computational effort, demanding a high amount of memory for storing the source's wavefields, consequently leading to a high cost to perform the imaging condition. Thus, this work aims to reduce these problems by introducing a dynamic approach (DA) that considers the sparsity of the wavefield in the first periods of propagation. The RTM combined with the DA (RTM-DA) approximates the computational domain to the propagation domain, which is the region delimited by the wavefront. In practical terms, the computational domain expands together with the wavefront, reflecting in a very significant economy of memory and reduction of processing time when compared to the conventional RTM, which we denominate as a static approach (RTM-SA). To reduce the sparsity of the wavefields in the first periods of propagation, we have built an empirical 3D filter that maps each timestep of the wavefield and gives the coordinates to approximate the computational domain to the propagation domain. We compare both approaches using the 3D SEG/EAGE Salt model and demonstrate that the RTM-DA is more efficient than the RTM-SA in terms of memory consumption and computational time, preserving the quality of the seismic image.

Three dimensional reverse time migration on common offset GPR data for void imaging

Three dimensional (3D) Reverse time migration (RTM), as an indispensable imaging technique, has been widely applied to the observed seismic data imaging. However, the RTM of Ground Penetrating radar (GPR) is usually implemented in its two-dimensional (2D) from. This paper proposes a 3D RTM algorithm on common offset GPR data for void imaging. In this algorithm, the 3D finite difference time domain (FDTD) method is applied to calculate the backward extrapolation of electromagnetic fields. And the time zero imaging condition is applied to obtain the final RTM results. Comparisons of the imaging results of 3D and 2D RTM demonstrated that the 3D RTM imaging results were superior to those of 2D RTM. And 3D spatial distribution of Void can be imaged by the proposed 3D RTM algorithm with high accuracy, which can help the subsequent interpretation of GPR profile.

Auto-Tuning of Dynamic Scheduling Applied to 3D Reverse Time Migration on Multicore Systems

Reverse time migration (RTM) is an algorithm widely used in the oil and gas industry to process seismic data. It is a computationally intensive task that suits well in parallel computers. Methods such as RTM can be parallelized in shared memory systems through scheduling iterations of parallel loops to threads. However, several aspects, such as memory size and hierarchy, number of cores, and input size, make optimal scheduling very challenging. In this paper, we introduce a run-time strategy to automatically tune the dynamic scheduling of parallel loops iterations in iterative applications, such as the RTM, in multicore systems. The proposed method aims to reduce the execution time of such applications. To find the optimal granularity, we propose a coupled simulated annealing (CSA) based auto-tuning strategy that adjusts the chunk size of work that OpenMP parallel loops assign dynamically to worker threads during the initialization of a 3D RTM application. Experiments performed with different computational systems and input sizes show that the proposed method is consistently better than the default OpenMP schedulers, static, auto, and guided, causing the application to be up to 33% faster. We show that the possible reason for this performance is the reduction of cache misses, mainly level L3, and low overhead, inferior to 2%. Having shown to be robust and scalable for the 3D RTM, the proposed method could also improve the performance of similar wave-based algorithms, such as full-waveform inversion (FWI) and other iterative applications.

Three-Dimensional Reverse Time Migration of Ground-Penetrating Radar Signals

Three-dimensional (3D) ground-penetrating radar (GPR) systems and 3D seismic imaging techniques have been developing fast and evolving rapidly in the last decade. Ray-based migration methods have been successfully applied to processing 3D GPR signals based on the similarity between electromagnetic and seismic waves. However, reverse time migration (RTM) of 3D GPR signals has not been well studied in the past. In this paper, we present a 3D RTM based on Maxwell’s equations for 3D GPR surveys. Migration recovers the true subsurface structure from a distorted and unfocused time profile and suppresses common electromagnetic clutter and noise. RTM based on Maxwell’s equations can consider conductivity directly and compensate for the attenuation within a high-conductivity zone. Compared with 2D RTM, 3D RTM back-propagates both in-line and cross-line signals simultaneously and can include complex 3D permittivity and conductivity models. We have integrated a parallel finite-difference time-domain (FDTD) algorithm based on a hybrid MPI and OpenMP scheme to reduce the computational cost of 3D problems. The 3D RTM experiments on an anomaly of “EM” shape and a realistic sand dune model demonstrate the effective recovery of 3D subsurface structures.

Distributed-Memory Load Balancing With Cyclic Token-Based Work-Stealing Applied to Reverse Time Migration

Reverse time migration (RTM) is a prominent technique in seismic imaging. Its resulting subsurface images are used in the industry to investigate with higher confidence the existence and the conditions of oil and gas reservoirs. Because of its high computational cost, RTM must make use of parallel computers. Balancing the workload distribution of an RTM is a growing challenge in distributed computing systems. The competition for shared resources and the differently-sized tasks of the RTM are some of the possible sources of load imbalance. Although many load balancing techniques exist, scaling up for large problems and large systems remains a challenge because synchronization overhead also scales. This paper proposes a cyclic token-based work-stealing (CTWS) algorithm for distributed memory systems applied to RTM. The novel cyclic token approach reduces the number of failed steals, avoids communication overhead, and simplifies the victim selection and the termination strategy. The proposed method is implemented as a C library using the one-sided communication feature of the message passing interface (MPI) standard. Results obtained by applying the proposed technique to balance the workload of a 3D RTM system present a factor of 14.1% speedup and reductions of the load imbalance of 78.4% when compared to the conventional static distribution.

Tilted transverse isotropic reverse time migration with angle gathers: Implementation and efficiency

We have developed a 3D reverse time migration (RTM) implementation, with an extended imaging condition, in the presence of tilted transverse isotropy (TTI), using finite-difference time-domain solvers on the rotated staggered grid (RSG) and the Lebedev grid. We have evaluated dispersion analysis of both schemes for a sample TTI medium. Using synthetic and real seismic data in realistic execution configurations, we found surprising inconsistency with quantitative cost estimates in the literature. For fixed accuracy, the RSG scheme proved most efficient in our tests, and not two to three times less efficient, as had been posited. Of course, variability may arise with particular medium properties, computer architectures, or software implementations. Having analyzed wavefield modeling performance in detail, we then quantified the realistic costs associated with the imaging condition in RTM, especially extended imaging conditions. For large extensions of the imaging condition, the associated cost dominates the total execution time and adds very significantly to the total cost of the migration. Our implementation targets both shared and distributed memory parallelism, making it suitable for general CPU clusters, nonuniform memory access architectures, and Intel Xeon-Phi. We have also discussed other aspects of our implementation in detail. In particular, we have explored an alternative method to achieve concurrent communication and computation in RTM.

The AWWE-based hybrid absorbing boundary condition for finite-difference modeling and its application in reverse-time migration

The reverse-time migration (RTM) crosscorrelation imaging condition requires that the forward-propagated source wavefield and the backward-propagated receiver wavefield must be obtained at the same time. The easiest way to get the source wavefield is to save the entire time history of the full wavefield into computer memory. However, this strategy requires huge amount of data storage. It is impossible for large-scale 3D RTM. To reduce the computer memory cost, the back-propagated source wavefield is reconstructed by using the stored boundary wavefield. Its computer memory is proportional to the saved boundary grid points. For high order of spatial finite-difference (FD) schemes, more boundary grid points are needed to be stored, which consumes a large amount of the computer memory required for RTM. To further reduce the computer memory cost, we adopt the hybrid absorbing boundary condition (ABC) combined with the arbitrarily wide-angle wave equations (AWWEs). In our method, three boundary grid points can obtain good absorption. The source wavefield can be accurately reconstructed by using these points and the mirror-image symmetry method. Numerical experiments demonstrate the correctness and effectiveness of the proposed method. We compared our method with the conventional hybrid ABC method based on the 15°one way wave equations (OWWEs). Comparisons show that our method with three boundary grid points can achieve the same absorption as the conventional method with ten boundary grid points. For twentieth order of accuracy in space, our method uses only about 30% of memory requirement and about 59% of computation time required by the conventional method.

3D weak-dispersion reverse time migration using a stereo-modeling operator

Reliable 3D imaging is a required tool for developing models of complex geologic structures. Reverse time migration (RTM), as the most powerful depth imaging method, has become the preferred imaging tool because of its ability to handle complex velocity models including steeply dipping interfaces and large velocity contrasts. Finite-difference methods are among the most popular numerical approaches used for RTM. However, these methods often encounter a serious issue of numerical dispersion, which is typically suppressed by reducing the grid interval of the propagation model, resulting in large computation and memory requirements. In addition, even with small grid spacing, numerical anisotropy may degrade images or, worse, provide images that appear to be focused but position events incorrectly. Recently, stereo-operators have been developed to approximate the partial differential operator in space. These operators have been used to develop several weak-dispersion and efficient stereo-modeling methods that have been found to be superior to conventional algorithms in suppressing numerical dispersion and numerical anisotropy. We generalized one stereo-modeling method, fourth-order nearly analytic central difference (NACD), from 2D to 3D and applied it to 3D RTM. The RTM results for the 3D SEG/EAGE phase A classic data set 1 and the SEG Advanced Modeling project model demonstrated that, even when using a large grid size, the NACD method can handle very complex velocity models and produced better images than can be obtained using the fourth-order and eighth-order Lax-Wendroff correction (LWC) schemes. We also applied 3D NACD and fourth-order LWC to a field data set and illustrated significant improvements in terms of structure imaging, horizon/layer continuity and positioning. We also investigated numerical dispersion and found that not only does the NACD method have superior dispersion characteristics but also that the angular variation of dispersion is significantly less than for LWC.

Reducing the computer memory requirement for 3D reverse-time migration with a boundary-wavefield extrapolation method

Reverse-time migration (RTM) using the crosscorrelation imaging condition requires that the forward-propagated source wavefield and the backward-propagated receiver wavefield be accessible within the imaging domain at the same time step. There are two categories of methods to balance the computer memory requirement and the computational complexity of RTM: checkpointing methods and source-wavefield reconstruction methods. We have developed a new source-wavefield reconstruction method to improve the balance between the computer memory requirement and the computational complexity of RTM. During the forward simulation of the source wavefield, we stored boundary wavefields only at one or two layers of spatial grid points and reconstructed the back-propagated source wavefield at the same time step as that of the back-propagated receiver wavefield, using a high-order wave-equation extrapolation scheme. One conventional RTM method uses boundary wavefields stored at multiple layers of spatial grid points and a high-order finite-difference (FD) scheme to reconstruct the back-propagated source wavefield. For an FD scheme with the eighth or sixteenth order of accuracy in space, our new method used only 37.5% of the computer memory required by this conventional method to store the boundary wavefields. This reduction of computer memory usage is significant because storing the boundary wavefields consumes most of the computer memory required for 3D migration using reconstructed source wavefields. Moreover, our method maintained the spatial order of accuracy of the FD scheme for the entire imaging domain, whereas some conventional methods reduce the spatial-order accuracy of the FD scheme near the boundaries to back-propagate the source wavefield to decrease the computer memory requirement. We validated our method using synthetic seismic data. Our method produced 2D and 3D migration images of complex subsurface structures as accurate as those yielded using an RTM method without reducing the spatial order of accuracy near the boundaries.

An efficient hybrid pseudospectral/finite-difference scheme for solving the TTI pure P-wave equation

The pure P-wave equation for modelling and migration in tilted transversely isotropic (TTI) media has attracted more and more attention in imaging seismic data with anisotropy. The desirable feature is that it is absolutely free of shear-wave artefacts and the consequent alleviation of numerical instabilities generally suffered by some systems of coupled equations. However, due to several forward–backward Fourier transforms in wavefield updating at each time step, the computational cost is significant, and thereby hampers its prevalence. We propose to use a hybrid pseudospectral (PS) and finite-difference (FD) scheme to solve the pure P-wave equation. In the hybrid solution, most of the cost-consuming wavenumber terms in the equation are replaced by inexpensive FD operators, which in turn accelerates the computation and reduces the computational cost. To demonstrate the benefit in cost saving of the new scheme, 2D and 3D reverse-time migration (RTM) examples using the hybrid solution to the pure P-wave equation are carried out, and respective runtimes are listed and compared. Numerical results show that the hybrid strategy demands less computation time and is faster than using the PS method alone. Furthermore, this new TTI RTM algorithm with the hybrid method is computationally less expensive than that with the FD solution to conventional TTI coupled equations.

2.5D reverse-time migration

Reverse-time migration (RTM) in 2.5D offers an alternative to improve resolution and amplitude when imaging 2D seismic data. Wave propagation in 2.5D assumes translational invariance of the velocity model. Under this assumption, we implement a finite-difference (FD) modeling algorithm in the mixed time-space/wavenumber domain to simulate the velocity and pressure fields for acoustic wave propagation and apply it in RTM. The 2.5D FD algorithm is truly parallel, allowing an efficient implementation in clusters. Storage and computing time requirements are strongly reduced compared to a full 3D FD simulation of the wave propagation. This feature makes 2.5D RTM much more efficient than 3D RTM, while achieving improved modeling of 3D geometrical spreading and phase properties of the seismic waveform in comparison to 2D RTM. Together with an imaging condition that compensates for uneven illumination and/or the obliquity factor, this allows recover of amplitudes proportional to the earth’s reflectivity. Numerical experiments using synthetic data demonstrate the better resolution and improved amplitude recovery of 2.5D RTM relative to 2D RTM.