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.展开更多
Strong fluctuation of seabed,abrupt variation in depth and dip of seabed bring seismic imaging problems,such as irregular reflection waves,obvious multiple waves,serious lateral wave development,poor imaging on base s...Strong fluctuation of seabed,abrupt variation in depth and dip of seabed bring seismic imaging problems,such as irregular reflection waves,obvious multiple waves,serious lateral wave development,poor imaging on base surface and depression structure,low signal-to-noise ratio of middle and deep layers.In this paper,Gaussian beam migration imaging method is used to analyze the imaging effect of rugged seabed in deep water area,and the ray tracing method of wavefront construction method is used to analyze the kinematic characteristics of seismic waves.By improving the design of seismic data acquisition and observation system,imaging quality of fine structures is improved.展开更多
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.展开更多
Gaussian beam migration （GBM） is an effec- tive and robust depth seismic imaging method, which overcomes the disadvantage of Kirchhoff migration in imaging multiple arrivals and has no steep-dip limitation of one-wa...Gaussian beam migration （GBM） is an effec- tive and robust depth seismic imaging method, which overcomes the disadvantage of Kirchhoff migration in imaging multiple arrivals and has no steep-dip limitation of one-way wave equation migration. However, its imaging quality depends on the initial beam parameters, which can make the beam width increase and wave-front spread with the propagation of the central ray, resulting in poor migration accuracy at depth, especially for exploration areas with complex geological structures. To address this problem, we present an adaptive focused beam method for shot-domain prestack depth migration. Using the infor- mation of the input smooth velocity field, we first derive an adaptive focused parameter, which makes a seismic beam focused along the whole central ray to enhance the wave- field construction accuracy in both the shallow and deep regions. Then we introduce this parameter into the GBM, which not only improves imaging quality of deep reflectors but also makes the shallow small-scale geological struc- tures well-defined. As well, using the amplitude-preserved extrapolation operator and deconvolution imaging condi- tion, the concept of amplitude-preserved imaging has been included in our method. Typical numerical examples and the field data processing results demonstrate the validity and adaptability of our method.展开更多
Since the wave equation of magnetoteiluric (MT)field is similar to the one of seismic , the migration techniques used in seismic can be applied to MT data . In this paper we make use of the principle of reflector mapp...Since the wave equation of magnetoteiluric (MT)field is similar to the one of seismic , the migration techniques used in seismic can be applied to MT data . In this paper we make use of the principle of reflector mapping (i. e. U/D imaging principle ) to image MT data . That is, the MT wavefield observed on the surface of the earth can be resolved into upgoing and downgoing waves , the waves are extrapolated downward by the phase - shift method or the phase - shift plus interpolation (PSPI )method . Conductivity interfaces of the medium could be found by using the time coincidence of the upgoing and downgoing waves . Theoretical calculations show that the migration technique of MT data presented here is very effective . It can not only enhance the lateral resolution of MT data , but also obtain the visual image of subsurface interfaces . As compared with the conventional 2 - D inversion , this procedure is more simple in calculation and can be easily put into practice on a personal computer and is able to obtain the MT depth section , which is similar to seismic section .展开更多
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 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）.
基金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.
基金Supported by projects of National Natural Science Foundation of China (No. 42074150)National Key R&D Program of China (No. 2017YFC0601305)。
文摘Strong fluctuation of seabed,abrupt variation in depth and dip of seabed bring seismic imaging problems,such as irregular reflection waves,obvious multiple waves,serious lateral wave development,poor imaging on base surface and depression structure,low signal-to-noise ratio of middle and deep layers.In this paper,Gaussian beam migration imaging method is used to analyze the imaging effect of rugged seabed in deep water area,and the ray tracing method of wavefront construction method is used to analyze the kinematic characteristics of seismic waves.By improving the design of seismic data acquisition and observation system,imaging quality of fine structures is improved.
文摘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.
文摘Gaussian beam migration （GBM） is an effec- tive and robust depth seismic imaging method, which overcomes the disadvantage of Kirchhoff migration in imaging multiple arrivals and has no steep-dip limitation of one-way wave equation migration. However, its imaging quality depends on the initial beam parameters, which can make the beam width increase and wave-front spread with the propagation of the central ray, resulting in poor migration accuracy at depth, especially for exploration areas with complex geological structures. To address this problem, we present an adaptive focused beam method for shot-domain prestack depth migration. Using the infor- mation of the input smooth velocity field, we first derive an adaptive focused parameter, which makes a seismic beam focused along the whole central ray to enhance the wave- field construction accuracy in both the shallow and deep regions. Then we introduce this parameter into the GBM, which not only improves imaging quality of deep reflectors but also makes the shallow small-scale geological struc- tures well-defined. As well, using the amplitude-preserved extrapolation operator and deconvolution imaging condi- tion, the concept of amplitude-preserved imaging has been included in our method. Typical numerical examples and the field data processing results demonstrate the validity and adaptability of our method.
文摘Since the wave equation of magnetoteiluric (MT)field is similar to the one of seismic , the migration techniques used in seismic can be applied to MT data . In this paper we make use of the principle of reflector mapping (i. e. U/D imaging principle ) to image MT data . That is, the MT wavefield observed on the surface of the earth can be resolved into upgoing and downgoing waves , the waves are extrapolated downward by the phase - shift method or the phase - shift plus interpolation (PSPI )method . Conductivity interfaces of the medium could be found by using the time coincidence of the upgoing and downgoing waves . Theoretical calculations show that the migration technique of MT data presented here is very effective . It can not only enhance the lateral resolution of MT data , but also obtain the visual image of subsurface interfaces . As compared with the conventional 2 - D inversion , this procedure is more simple in calculation and can be easily put into practice on a personal computer and is able to obtain the MT depth section , which is similar to seismic section .
基金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.
