The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward c...The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations.展开更多
Prestack depth migration for seismic reflection data is commonly used tool for imaging complex geological structures such as salt domes, faults, thrust belts, and stratigraphic structures. Phase shift plus interpolati...Prestack depth migration for seismic reflection data is commonly used tool for imaging complex geological structures such as salt domes, faults, thrust belts, and stratigraphic structures. Phase shift plus interpolation (PSPI) algorithm is a useful tool to directly solve a wave equation and the results have natural properties of the wave equation. Amplitude and phase characteristics, in particular, are better preserved. The PSPI algorithm is widely used in hydrocarbon exploration because of its simplicity, efficiency, and reduced efforts for computation. However, meaningful depth image of 3D subsurface requires parallel computing to handle heavy computing time and great amount of input data. We implemented a parallelized version of 3D PSPI for prestack depth migration using Open-Multi-Processing (Open MP) library. We verified its performance through applications to 3D SEG/EAGE salt model with a small scale Linux cluster. Phase-shift was performed in the vertical and horizontal directions, respectively, and then interpolated at each node. This gave a single image gather according to shot gather. After summation of each single image gather, we got a 3D stacked image in the depth domain. The numerical model example shows good agree- ment with the original geological model.展开更多
With increased computational power, reverse-time prestack depth migration(RT-PSDM) has become a preferred imaging tool in seismic exploration, yet its use has remained relatively limited in ground-penetrating radar...With increased computational power, reverse-time prestack depth migration(RT-PSDM) has become a preferred imaging tool in seismic exploration, yet its use has remained relatively limited in ground-penetrating radar(GPR) applications. Complex topography alters the wavefield kinematics making for a challenging imaging problem. Model simulations show that topographic variation can substantially distort reflection amplitudes due to irregular wavefield spreading, attenuation anomalies due to irregular path lengths, and focusing and defocusing effects at the surface. The effects are magnified when the topographic variations are on the same order as the depth of investigation––a situation that is often encountered in GPR investigations. Here, I use a full wave-equation RT-PSDM algorithm to image GPR data in the presence of large topographic variability relative to the depth of investigation. The source and receiver wavefields are propagated directly from the topographic surface and this approach inherently corrects for irregular kinematics, spreading and attenuation. The results show that when GPR data are acquired in areas of extreme topography, RT-PSDM can accurately reconstruct reflector geometry as well as reflection amplitude.展开更多
In this paper the methods of wave theory based prestack depth migration and their implementation are studied. Using the splitting of wave operator, the wavefield extrapolation equations are deduced and the numerical s...In this paper the methods of wave theory based prestack depth migration and their implementation are studied. Using the splitting of wave operator, the wavefield extrapolation equations are deduced and the numerical schemes are presented. The numerical tests for SEG/EAEG model with MPI are performed on the PC-cluster. The numerical results show that the methods of single-shot (common-shot) migration and synthesized-shot migration are of practical values and can be applied to field data processing of 3D prestack depth migration.展开更多
Kirchhoff beam migration is a beam migration method, which focuses on rapid imaging of geological structures. Although this imaging method ignores the amplitude information in the calculation process, it can calculate...Kirchhoff beam migration is a beam migration method, which focuses on rapid imaging of geological structures. Although this imaging method ignores the amplitude information in the calculation process, it can calculate multi-arrival traveltime. This migration method takes into account both imaging accuracy and computational efficiency. Kirchhoff beam migration employs coarse grid techniques in several key steps such as traveltime calculation, weight function calculation, and imaging calculation. The selection of the coarse mesh size has an important influence on the computational efficiency and imaging accuracy of the migration imaging method. This paper will analyze this influence and illustrate the analysis results by the Marmousi data sets.展开更多
基金This paper is supported by the National Natural Science Foundation of China (No. 40474047)State Key Laboratory of Geological Processes and Mineral Resources (No. GPMR200654)the Focused Subject Program of Beijing (No. XK104910598).
文摘The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations.
基金This project is sponsored by the National Natural Science Foundation (40474041), CNPC Young Innovation Fund (04E7040), the Post-doctoral Research Station of Zhongyuan 0ilfield, Jiangsu 0ilfield, and CNPC Geophysical Laboratories at the China University of Petroleum (East China).
文摘Prestack depth migration for seismic reflection data is commonly used tool for imaging complex geological structures such as salt domes, faults, thrust belts, and stratigraphic structures. Phase shift plus interpolation (PSPI) algorithm is a useful tool to directly solve a wave equation and the results have natural properties of the wave equation. Amplitude and phase characteristics, in particular, are better preserved. The PSPI algorithm is widely used in hydrocarbon exploration because of its simplicity, efficiency, and reduced efforts for computation. However, meaningful depth image of 3D subsurface requires parallel computing to handle heavy computing time and great amount of input data. We implemented a parallelized version of 3D PSPI for prestack depth migration using Open-Multi-Processing (Open MP) library. We verified its performance through applications to 3D SEG/EAGE salt model with a small scale Linux cluster. Phase-shift was performed in the vertical and horizontal directions, respectively, and then interpolated at each node. This gave a single image gather according to shot gather. After summation of each single image gather, we got a 3D stacked image in the depth domain. The numerical model example shows good agree- ment with the original geological model.
基金The Herbette Fondation at the University of Lausanne, Switzerland
文摘With increased computational power, reverse-time prestack depth migration(RT-PSDM) has become a preferred imaging tool in seismic exploration, yet its use has remained relatively limited in ground-penetrating radar(GPR) applications. Complex topography alters the wavefield kinematics making for a challenging imaging problem. Model simulations show that topographic variation can substantially distort reflection amplitudes due to irregular wavefield spreading, attenuation anomalies due to irregular path lengths, and focusing and defocusing effects at the surface. The effects are magnified when the topographic variations are on the same order as the depth of investigation––a situation that is often encountered in GPR investigations. Here, I use a full wave-equation RT-PSDM algorithm to image GPR data in the presence of large topographic variability relative to the depth of investigation. The source and receiver wavefields are propagated directly from the topographic surface and this approach inherently corrects for irregular kinematics, spreading and attenuation. The results show that when GPR data are acquired in areas of extreme topography, RT-PSDM can accurately reconstruct reflector geometry as well as reflection amplitude.
基金This work was supported by Major State Basic Research Program of Peoples's Republic of China(No.G1999032800)Major Project(No.49894190)the National Natural Science Foundation of China(Grant No.40004003).All numerical experiments were completed on the PC-cluster in the State Key Lab of Scientific/Engineering Computing.
文摘In this paper the methods of wave theory based prestack depth migration and their implementation are studied. Using the splitting of wave operator, the wavefield extrapolation equations are deduced and the numerical schemes are presented. The numerical tests for SEG/EAEG model with MPI are performed on the PC-cluster. The numerical results show that the methods of single-shot (common-shot) migration and synthesized-shot migration are of practical values and can be applied to field data processing of 3D prestack depth migration.
基金Supported by projects of the Natural Science Foundation of China(No.41804100)the China Postdoctoral Science Foundation(No.2018M640910)the Fundamental Research Funds for the Central Universities(No.2682018CX36)
文摘Kirchhoff beam migration is a beam migration method, which focuses on rapid imaging of geological structures. Although this imaging method ignores the amplitude information in the calculation process, it can calculate multi-arrival traveltime. This migration method takes into account both imaging accuracy and computational efficiency. Kirchhoff beam migration employs coarse grid techniques in several key steps such as traveltime calculation, weight function calculation, and imaging calculation. The selection of the coarse mesh size has an important influence on the computational efficiency and imaging accuracy of the migration imaging method. This paper will analyze this influence and illustrate the analysis results by the Marmousi data sets.
基金supported by the National Natural Science Foundation of China(No.62175100)Spark Program of Earthquake Sciences of CEA(No.XH22015A)+1 种基金Henan Province Seismic Structure Exploration project(YCZC-2020-950)Special Fund of Chinese Central Government for Basic Scientific Research Operations in Commonweal Research Institutes(No.IGCEA1902)。