The numerical dispersion phenomenon in the finite-difference forward modeling simulations of the wave equation significantly affects the imaging accuracy in acoustic reflection logging.This issue is particularly prono...The numerical dispersion phenomenon in the finite-difference forward modeling simulations of the wave equation significantly affects the imaging accuracy in acoustic reflection logging.This issue is particularly pronounced in the reverse time migration(RTM)method used for shear-wave(S-wave)logging imaging.This not only affects imaging accuracy but also introduces ambiguities in the interpretation of logging results.To address this challenge,this study proposes the use of a least-squares difference coefficient optimization algorithm aiming to suppress the numerical dispersion phenomenon in the RTM of S-wave reflection imaging logging.By optimizing the difference coefficients,the high-precision finite-difference algorithm serves as an effective operator for both forward and backward RTM processes.This approach is instrumental in eliminating migration illusions,which are often caused by numerical dispersion.The effectiveness of this optimized algorithm is demonstrated through numerical results,which indicate that it can achieve more accurate forward imaging results across various conditions,including high-and low-velocity strata,and is effective in both large and small spatial grids.The results of processing real data demonstrate that numerical dispersion optimization effectively reduces migration artifacts and diminishes ambiguities in logging interpretations.This optimization offers crucial technical support to the RTM method,enhancing its capability for accurately modeling and imaging S-wave reflections.展开更多
Pre-stack depth migration velocity analysis is one of the keys to influencing the imaging quality of pre-stack migration.In this paper we cover a residual curvature velocity analysis method on angle-domain common imag...Pre-stack depth migration velocity analysis is one of the keys to influencing the imaging quality of pre-stack migration.In this paper we cover a residual curvature velocity analysis method on angle-domain common image gathers（ADCIGs） which can depict the relationship between incident angle and migration depth at imaging points and update the migration velocity.Differing from offset-domain common image gathers（ODCIGs）,ADCIGs are not disturbed by the multi-path problem which contributes to imaging artifacts,thus influencing the velocity analysis.On the basis of horizontal layers,we derive the residual depth equation and also propose a velocity analysis workflow for velocity scanning.The tests to synthetic and field data prove the velocity analysis methods adopted in this paper are robust and valid.展开更多
The imaging of offset VSP data in local phase space can improve the image of the subsurface structure near the well.In this paper,we present a migration scheme for imaging VSP data in a local phase space,which uses th...The imaging of offset VSP data in local phase space can improve the image of the subsurface structure near the well.In this paper,we present a migration scheme for imaging VSP data in a local phase space,which uses the Gabor-Daubechies tight framebased extrapolator（G-D extrapolator） and its high-frequency asymptotic expansion to extrapolate wavefields and also delineates an improved correlation imaging condition in the local angle domain.The results for migrating synthetic and real VSP data demonstrate that the application of the high-frequency G-D extrapolator asymptotic expansion can effectively decrease computational complexity.The local angle domain correlation imaging condition can be used to weaken migration artifacts without increasing computation.展开更多
In wave-equation migration and demigration,the cross-correlation imaging/forwarding step implicitly injects an additional copy of the source wavelet,so that the amplitude spectrum of the wavelet is applied redundantly...In wave-equation migration and demigration,the cross-correlation imaging/forwarding step implicitly injects an additional copy of the source wavelet,so that the amplitude spectrum of the wavelet is applied redundantly(effectively imposing a wavelet-spectrum weighting,often akin to an amplitude-squared bias).This redundancy degrades structural fidelity and amplitude balance yet is frequently overlooked.We(i)formalize the mechanism by which cross-correlation duplicates the source-wavelet amplitude effect in both migration and demigration,and(ii)introduce a source-equalized operator that removes the redundancy by deconvolving(or dividing by)the wavelet amplitude spectrum in the imaging condition and its demigration counterpart,while leaving phase/kinematics intact.Using a band-limited Ricker wavelet on a two-layer model and on Marmousi,we show that,if unmanaged,the redundant wavelet spectrum broadens main lobes,introduces ringing,and suppresses vertical resolution in migrated images,and inflates spectrum mismatches between demigrated and observed data even when peak times agree.With our correction,images recover observed-data-consistent bandwidth and sharpened interfaces,and demigrated data also exhibit improved spectrum conformity and reduced amplitude misfit.The results clarify when source amplitudes matter,why cross-correlation makes them redundantly matter,and how a lightweight spectral correction restores physically meaningful amplitude behavior in wave-equation migration/demigration.展开更多
Mesenchymal stromal cell transplantation is an effective and promising approach for treating various systemic and diffuse diseases.However,the biological characteristics of transplanted mesenchymal stromal cells in hu...Mesenchymal stromal cell transplantation is an effective and promising approach for treating various systemic and diffuse diseases.However,the biological characteristics of transplanted mesenchymal stromal cells in humans remain unclear,including cell viability,distribution,migration,and fate.Conventional cell tracing methods cannot be used in the clinic.The use of superparamagnetic iron oxide nanoparticles as contrast agents allows for the observation of transplanted cells using magnetic resonance imaging.In 2016,the National Medical Products Administration of China approved a new superparamagnetic iron oxide nanoparticle,Ruicun,for use as a contrast agent in clinical trials.In the present study,an acute hemi-transection spinal cord injury model was established in beagle dogs.The injury was then treated by transplantation of Ruicun-labeled mesenchymal stromal cells.The results indicated that Ruicunlabeled mesenchymal stromal cells repaired damaged spinal cord fibers and partially restored neurological function in animals with acute spinal cord injury.T2*-weighted imaging revealed low signal areas on both sides of the injured spinal cord.The results of quantitative susceptibility mapping with ultrashort echo time sequences indicated that Ruicun-labeled mesenchymal stromal cells persisted stably within the injured spinal cord for over 4 weeks.These findings suggest that magnetic resonance imaging has the potential to effectively track the migration of Ruicun-labeled mesenchymal stromal cells and assess their ability to repair spinal cord injury.展开更多
We build a high-resolution early aftershock catalog for the 2023 SE Türkiye seismic sequence with PALM,a seamless workflow that sequentially performs phase picking,association,location,and matched filter for cont...We build a high-resolution early aftershock catalog for the 2023 SE Türkiye seismic sequence with PALM,a seamless workflow that sequentially performs phase picking,association,location,and matched filter for continuous data.The catalog contains 29,519 well-located events in the two mainshocks rupture region during 2023-02-01–2023-02-28,which significantly improves the detection completeness and relocation precision compared to the public routine catalog.Employing the new PALM catalog,we analyze the structure of the seismogenic fault system.We find that the Eastern Anatolian Fault(EAF)that generated the first M_(W)7.9 mainshock is overall near-vertical,whereas complexities are revealed in a small-scale,such as subparallel subfaults,unmapped branches,and stepovers.The seismicity on EAF is shallow(<15 km)and concentrated in depth distribution,indicating a clear lock-creep transition.In contrast,the SürgüFault(SF)that is responsible for the second M_(W)7.8 mainshock is shovel-shaped for the nucleation segment and has overall low dip angles(~40°–80°).Aftershocks on the SF distribute in a broad range of depth,extending down to~35 km.We also analyze the temporal behavior of seismicity,discovering no immediate foreshocks within~5 days preceding the first mainshock,and no seismic activity on the SF before the second mainshock.展开更多
Correctly locating the tunnel lining cavity is extremely important tunnel quality inspection.High-accuracy imaging results are hard to obtain because conventional one-way wave migration is greatly aff ected by lateral...Correctly locating the tunnel lining cavity is extremely important tunnel quality inspection.High-accuracy imaging results are hard to obtain because conventional one-way wave migration is greatly aff ected by lateral velocity change and inclination limitation and because the diff racted wave cannot be accurately returned to the real spatial position of the lining cavity.This paper presents a tunnel lining cavity imaging method based on the groundpenetrating radar(GPR)reverse-time migration(RTM)algorithm.The principle of GPR RTM is described in detail using the electromagnetic wave equation.The finite-difference timedomain method is employed to calculate the backward extrapolation electromagnetic fi elds,and the zero-time imaging condition based on the exploding-reflector concept is used to obtain the RTM results.On this basis,the GPR RTM program is compiled and applied to the simulated and observed GPR data of a typical tunnel lining cavity GPR model and a physical lining cavity model.Comparison of RTM and Kirchhoff migration results reveals that the RTM can better converge the diff racted waves of steel bar and cavity to their true position and have higher resolution and better suppress the eff ect of multiple interference and clutter scattering waves.In addition,comparison of RTM results of diff erent degrees of noise shows that RTM has strong anti-interference ability and can be used for the accurate interpretation of radar profi le in a strong interference environment.展开更多
One of the problems experienced in marine geophysical exploration is that the layered features in the migration imaging profile are blurred and the seismic energy reflected is weaker in the middle or lower parts. In t...One of the problems experienced in marine geophysical exploration is that the layered features in the migration imaging profile are blurred and the seismic energy reflected is weaker in the middle or lower parts. In this study we model the seismic wavefield records in the undulating seafloor when there is both a slight change and significant change in seafloor topography to analyze its influence on the seismic reflection data and migration imaging profiles. We compare and analyze the wave field records collected at the same point on the original and modified velocity models, and the cross-bonding resulting migration imaging profiles. The results show that whether the seismic reflection data collection is performed along the direction of the survey line or against the direction of the survey line, slight changes in the seafloor topography have little effect on the wave field records and the migration profile, while significant changes in the seafloor topography have great effect on both the wave field records and migration profile.展开更多
With the development of seismic exploration,passive-source seismic data has attracted increasing attention.Ambient noise passive seismic sources exists widely in nature and industrial production.Passive seismic data i...With the development of seismic exploration,passive-source seismic data has attracted increasing attention.Ambient noise passive seismic sources exists widely in nature and industrial production.Passive seismic data is important in logging while drilling(LWD),large-scale structural exploration,etc.In this paper,we proposed a passive multiple reverse time migration imaging(PMRTMI)method based on wavefield decomposition and normalized imaging conditions method.This method differs from seismic interferometry in that it can use raw passive seismic data directly in RTM imaging without reconstruction of virtual active gather,and we use the wavefield decomposition method to eliminate the low frequency noise in RTM.Further,the energy normalized imaging condition is used in full wavefield decomposition,which can not only enhance the image quality of both edge and deep information but also overcome the wrong energy problem caused by uneven distribution of passive sources;furthermore,this method exhibits high efficiency.Finally,numerical examples with the Marmousi model show the effectiveness of the method.展开更多
The high frame rate(HFR)imaging technique requires only one emission event for imaging.Therefore,it can achieve ultrafast imaging with frame rates up to the kHz regime,which satisfies the frame rate requirements for i...The high frame rate(HFR)imaging technique requires only one emission event for imaging.Therefore,it can achieve ultrafast imaging with frame rates up to the kHz regime,which satisfies the frame rate requirements for imaging moving tissues in scientific research and clinics.Lu’s Fourier migration method is based on a non-diffraction beam to obtain HFR images and can improve computational speed and efficiency.However,in order to obtain high-quality images,Fourier migration needs to make full use of the spectrum of echo signals for imaging,which requires a large number of Fast Fourier Transform(FFT)points and increases the complexity of the hardware when the echo frequency is high.Here,an efficient algorithm using the spectrum migration technique based on the spectrum’s distribution characteristics is proposed to improve the imaging efficiency in HFR imaging.Since the actual echo signal spectrum is of limited bandwidth,low-frequency and high-frequency parts with low-energy have little contribution to the imaging spectrum.We transform the effective part that provides the main energy in the signal spectrum to the imaging spectrum while the ineffective spectrum components are not utilized for imaging.This can significantly reduce the number of Fourier transform points,improve Fourier imaging efficiency,and ensure the imaging quality.The proposed method is evaluated on simulated and experimental datasets.Results demonstrated that the proposed method could achieve equivalent image quality with a reduced point number for FFT compared to the complete spectrum migration.In this paper,it only requires a quarter of the FFT points used in the complete spectrum migration,which can improve the computational efficiency;thus,it is more suitable for real-time data processing.The proposed spectrum migration method has a specific significance for the study and clinical application of HFR imaging.展开更多
In view of the seismic exploration problem of thin sand reservoirs in the Songliao Basin, this paper puts forward a migration imaging method using CGP （common geophone point） stacked cylindrical waves. By this means...In view of the seismic exploration problem of thin sand reservoirs in the Songliao Basin, this paper puts forward a migration imaging method using CGP （common geophone point） stacked cylindrical waves. By this means, seismic data should be acquired from a midpoint shooting layout system with small shot-point spacing and small geophone interval. Using such seismic data, CGP gathers are first stacked to compose a cylindrical wave section. The cylindrical wave section is migrated and imaged by means of the ray path downward continuation of the down-going wave and the wave equation downward continuation of the upgoing wave. The results from the modeling analysis and the data processing of the TK8157 seismic line in the Songliao Basin shows that the proposed migration imaging method has higher seismic resolution and fidelity. Furthermore, the proposed method is proven to be more effective for discovering small sand bodies, small faults, stratigraphic pinch-outs, and so on.展开更多
With the increasing complexity of prospecting objectives,reverse time migration( RTM) has attracted more and more attention due to its outstanding imaging quality. RTM is based on two-way wave equation,so it can avoid...With the increasing complexity of prospecting objectives,reverse time migration( RTM) has attracted more and more attention due to its outstanding imaging quality. RTM is based on two-way wave equation,so it can avoid the limits of angle in traditional one-way wave equation migration,image reverse branch,prism waves and multi-reflected wave precisely and obtain accurate dynamic information. However,the huge demands for storage and computation as well as low frequency noises restrict its wide application. The normalized cross-correlation imaging conditions based on wave field decomposition are derived from traditional cross-correlation imaging condition,and it can eliminate the low-frequency noises effectively and improve the imaging resolution. The practical procedure includes separating source and receiver wave field into one-way components respectively,and conducting cross-correlation imaging condition to the post-separated wave field. In this way,the resolution and precision of the imaging result will be promoted greatly.展开更多
Based on the study on electromagnetic field migration by Zhdanov, we have proposed an improved method for the weak points in the research. Firstly, the initial background resistivity should be determined by using 1-D ...Based on the study on electromagnetic field migration by Zhdanov, we have proposed an improved method for the weak points in the research. Firstly, the initial background resistivity should be determined by using 1-D inversion results. Then in the process of continuation, the results are corrected and calculated layer by layer by the iteration method, so that more exact resistivity can be obtained. Secondly, an improved algorithm for finite-difference equation is studied. According to the property of electromagnetic migration field, the algorithm is designed by means of grids varying with geometric progression in the longitudinal direction. Being improved by the techniques mentioned above, better results are obtained by the new method, which has been verified by both the theory model and practical data.展开更多
In seismic exploration,it is a critical task to image and interpret different seismic signatures over complex geology due to strong lateral velocity contrast,steep reflectors,overburden strata and dipping flanks.To un...In seismic exploration,it is a critical task to image and interpret different seismic signatures over complex geology due to strong lateral velocity contrast,steep reflectors,overburden strata and dipping flanks.To understand the behavior of these seismic signatures,nowadays Reverse Time Migration(RTM)technique is used extensively by the oil&gas industries.During the extrapolation phase of RTM,the source wavefield needs to be saved,which needs high storage memory and large computing time.These two are the main obstacles of RTM for production use.In order to overcome these disadvantages,in this study,a second-generation improved RTM technique is proposed.In this improved form,a shift operator is introduced at the time of imaging condition of RTM algorithm which is performed automatically both in space and time domain.This effort is made to produce a better-quality image by minimizing the computational time as well as numerical artefacts.The proposed method is applied over various benchmark models and validated by implementing over one field data set from the Jaisalmer Basin,India.From the analysis,it is observed that the method consumes a minimum of 45%less storage space and reduce the execution time by 20%,as compared to conventional RTM.The proposed RTM is found to work efficiently in comparison to the conventional RTM both in terms of imaging quality and minimization of numerical artefacts for all the benchmark models as well as field data.展开更多
The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal ...The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal and produce a high precision subsurface image when the object is buried in a complicated subsurface environment.To better handle the multi-off set GPR data,we propose a three-dimensional(3D)prestack RTM algorithm.The high-order fi nite diff erence time domian(FDTD)method,with the accuracy of eighth-order in space and second-order in time,is applied to simulate the forward and backward extrapolation electromagnetic fi elds.In addition,we use the normalized correlation imaging condition to obtain pre-stack RTM result and the Laplace fi lter to suppress the low frequency noise generated during the correlation process.The numerical test of 3D simulated GPR data demonstrated that 3D RTM image shows excellent coincidence with the true model.Compared with 2D RTM image,the 3D RTM image can more clearly and accurately refl ect the 3D spatial distribution of the target,and the resolution of the imaging results is far better.Furthermore,the application of observed GPR data further validates the eff ectiveness of the proposed 3D GPR RTM algorithm,and its fi nal image can more reliably guide the subsequent interpretation.展开更多
Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by whi...Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by which static correction is completely abandoned before migration and surely the imaging quality is remarkably improved. To obtain some artificial shot gathers, high-order staggered-grid finite-difference （FD） method is adapted to model acoustic wave propagation. Since the shot gathers are always disturbed by regular interferences, the statics still must be applied to supporting the interference elimination by apparent velocity filtering method. Then all the shot gathers should be removed back to their original positions by reverse statics. Finally, they are migrated by pre-stack reverse-time depth migration and imaged. The numerical experiments show that the DPDM can ideally avoid the mistakes caused by statics and increase imaging precision.展开更多
Anisotropy correction is necessary during the processing of converted PS- wave seismic data to achieve accurate structural imaging, reservoir prediction, and fracture detection. To effectively eliminate the adverse ef...Anisotropy correction is necessary during the processing of converted PS- wave seismic data to achieve accurate structural imaging, reservoir prediction, and fracture detection. To effectively eliminate the adverse effects of S-wave splitting and to improve PS- wave imaging quality, we tested methods for pre-stack migration imaging and anisotropic correction of PS-wave data. We based this on the propagation rules of seismic waves in a horizontal transverse isotropy medium, which is a fractured medium model that reflects likely subsurface conditions in the field. We used the radial （R） and transverse （T） components of PS-wave data to separate the fast and slow S-wave components, after which their propagation moveout was effectively extracted. Meanwhile, corrections for the energies and propagation moveouts of the R and T components were implemented using mathematical rotation. The PS-wave imaging quality was distinctly improved, and we demonstrated the reliability of our methods through numerical simulations. Applying our methods to three-dimensional and three-component seismic field data from the Xinchang-Hexingchang region of the Western Sichuan Depression in China, we obtained high-quality seismic imaging with continuous reflection wave groups, distinct structural features, and specific stratigraphic contact relationships. This study provides an effective and reliable approach for data processing that will improve the exploration of complex, hidden lithologic gas reservoirs.展开更多
Kirchhoff beam migration is a simplified Gaussian beam migration,which omits the dynamic information and can calculate multi-arrival traveltime,so it is a high-precision and fast seismic imaging method.In the imaging ...Kirchhoff beam migration is a simplified Gaussian beam migration,which omits the dynamic information and can calculate multi-arrival traveltime,so it is a high-precision and fast seismic imaging method.In the imaging process,extracting common image gathers can be used for velocity analysis,improving the accuracy of modeling and imaging quality.Compared with the conventional common image gathers extracting methods,the angle-domain common image gathers extracting method can avoid the artifacts caused by multi-arrival seismic waves.The authors present a new method of extracting common image gathers in angle-domain from Kirchhoff beam migration and verify the method by numerical calculations.展开更多
Reasonable field acquisition geometry can not only guide seismic exploration to obtain sufficient geological information of target body,but also reduce acquisition cost to the maximum.In this study,building on convent...Reasonable field acquisition geometry can not only guide seismic exploration to obtain sufficient geological information of target body,but also reduce acquisition cost to the maximum.In this study,building on conventional ray-based geometry design methods,we incorporate imaging results as a constraint to optimize the geometry design and evaluate its effectiveness.Firstly,the geological model of the target layer is established based on the geological data of the study area and surface seismic data combined with exploration tasks.Then,the ray-tracing method is employed to simulate and assess the proposed geometry design,verifying whether its parameters meet the exploration requirements.Finally,the imaging effect of the designed geometry on the target layer is tested by the cross-well seismic reverse time migration method.This methodology was applied to design the cross-well seismic acquisition geometry for offshore deviated wells in the X Oilfield.The simulation results demonstrate that the imaging-driven geometry design approach effectively guides field operations,enhances the imaging quality of the target layer,and reduces acquisition costs.展开更多
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.展开更多
基金supported by Scientific Research and Technology Development Project of CNPC(2021DJ4002,2022DJ3908).
文摘The numerical dispersion phenomenon in the finite-difference forward modeling simulations of the wave equation significantly affects the imaging accuracy in acoustic reflection logging.This issue is particularly pronounced in the reverse time migration(RTM)method used for shear-wave(S-wave)logging imaging.This not only affects imaging accuracy but also introduces ambiguities in the interpretation of logging results.To address this challenge,this study proposes the use of a least-squares difference coefficient optimization algorithm aiming to suppress the numerical dispersion phenomenon in the RTM of S-wave reflection imaging logging.By optimizing the difference coefficients,the high-precision finite-difference algorithm serves as an effective operator for both forward and backward RTM processes.This approach is instrumental in eliminating migration illusions,which are often caused by numerical dispersion.The effectiveness of this optimized algorithm is demonstrated through numerical results,which indicate that it can achieve more accurate forward imaging results across various conditions,including high-and low-velocity strata,and is effective in both large and small spatial grids.The results of processing real data demonstrate that numerical dispersion optimization effectively reduces migration artifacts and diminishes ambiguities in logging interpretations.This optimization offers crucial technical support to the RTM method,enhancing its capability for accurately modeling and imaging S-wave reflections.
基金supported by the National 863 Program (Grant No.2006AA06Z206,Sustained supported)the National Science and Technology Major Project (Grant No.2008ZX05006-004)SinoPec Group Marine Facies Research (Grant No.08370502000410)
文摘Pre-stack depth migration velocity analysis is one of the keys to influencing the imaging quality of pre-stack migration.In this paper we cover a residual curvature velocity analysis method on angle-domain common image gathers（ADCIGs） which can depict the relationship between incident angle and migration depth at imaging points and update the migration velocity.Differing from offset-domain common image gathers（ODCIGs）,ADCIGs are not disturbed by the multi-path problem which contributes to imaging artifacts,thus influencing the velocity analysis.On the basis of horizontal layers,we derive the residual depth equation and also propose a velocity analysis workflow for velocity scanning.The tests to synthetic and field data prove the velocity analysis methods adopted in this paper are robust and valid.
基金supported by the National Hi-Tech Research and Development Program of China (Grant No.2006AA09A102-11)the National Natural Science Fund of China (Grant No.40730424 and 40674064)
文摘The imaging of offset VSP data in local phase space can improve the image of the subsurface structure near the well.In this paper,we present a migration scheme for imaging VSP data in a local phase space,which uses the Gabor-Daubechies tight framebased extrapolator（G-D extrapolator） and its high-frequency asymptotic expansion to extrapolate wavefields and also delineates an improved correlation imaging condition in the local angle domain.The results for migrating synthetic and real VSP data demonstrate that the application of the high-frequency G-D extrapolator asymptotic expansion can effectively decrease computational complexity.The local angle domain correlation imaging condition can be used to weaken migration artifacts without increasing computation.
基金supported by the National Natural Science Foundation of China(42430303)Strategy Priority Research Program(Category B)of the Chinese Academy of Sciences(XDB0710000)+2 种基金National Natural Science Foundation of China(42288201)the National Key R&D Program of China(2023YFF0803203)the IGGCAS start-up funding(Grant No.E251510101).
文摘In wave-equation migration and demigration,the cross-correlation imaging/forwarding step implicitly injects an additional copy of the source wavelet,so that the amplitude spectrum of the wavelet is applied redundantly(effectively imposing a wavelet-spectrum weighting,often akin to an amplitude-squared bias).This redundancy degrades structural fidelity and amplitude balance yet is frequently overlooked.We(i)formalize the mechanism by which cross-correlation duplicates the source-wavelet amplitude effect in both migration and demigration,and(ii)introduce a source-equalized operator that removes the redundancy by deconvolving(or dividing by)the wavelet amplitude spectrum in the imaging condition and its demigration counterpart,while leaving phase/kinematics intact.Using a band-limited Ricker wavelet on a two-layer model and on Marmousi,we show that,if unmanaged,the redundant wavelet spectrum broadens main lobes,introduces ringing,and suppresses vertical resolution in migrated images,and inflates spectrum mismatches between demigrated and observed data even when peak times agree.With our correction,images recover observed-data-consistent bandwidth and sharpened interfaces,and demigrated data also exhibit improved spectrum conformity and reduced amplitude misfit.The results clarify when source amplitudes matter,why cross-correlation makes them redundantly matter,and how a lightweight spectral correction restores physically meaningful amplitude behavior in wave-equation migration/demigration.
基金supported by the National Key R&D Program of China,Nos.2017YFA0104302(to NG and XM)and 2017YFA0104304(to BW and ZZ)
文摘Mesenchymal stromal cell transplantation is an effective and promising approach for treating various systemic and diffuse diseases.However,the biological characteristics of transplanted mesenchymal stromal cells in humans remain unclear,including cell viability,distribution,migration,and fate.Conventional cell tracing methods cannot be used in the clinic.The use of superparamagnetic iron oxide nanoparticles as contrast agents allows for the observation of transplanted cells using magnetic resonance imaging.In 2016,the National Medical Products Administration of China approved a new superparamagnetic iron oxide nanoparticle,Ruicun,for use as a contrast agent in clinical trials.In the present study,an acute hemi-transection spinal cord injury model was established in beagle dogs.The injury was then treated by transplantation of Ruicun-labeled mesenchymal stromal cells.The results indicated that Ruicunlabeled mesenchymal stromal cells repaired damaged spinal cord fibers and partially restored neurological function in animals with acute spinal cord injury.T2*-weighted imaging revealed low signal areas on both sides of the injured spinal cord.The results of quantitative susceptibility mapping with ultrashort echo time sequences indicated that Ruicun-labeled mesenchymal stromal cells persisted stably within the injured spinal cord for over 4 weeks.These findings suggest that magnetic resonance imaging has the potential to effectively track the migration of Ruicun-labeled mesenchymal stromal cells and assess their ability to repair spinal cord injury.
基金jointly supported by the National Key R&D Program (No.2022YFF0800601)the Istanbul Technical University Research Fund (ITU-BAP)+1 种基金the Alexander von Humboldt Foundation Research Fellowship Award for providing computing facilities through the Humboldt-Stiftung Follow-Up Programthe University of California,Riverside。
文摘We build a high-resolution early aftershock catalog for the 2023 SE Türkiye seismic sequence with PALM,a seamless workflow that sequentially performs phase picking,association,location,and matched filter for continuous data.The catalog contains 29,519 well-located events in the two mainshocks rupture region during 2023-02-01–2023-02-28,which significantly improves the detection completeness and relocation precision compared to the public routine catalog.Employing the new PALM catalog,we analyze the structure of the seismogenic fault system.We find that the Eastern Anatolian Fault(EAF)that generated the first M_(W)7.9 mainshock is overall near-vertical,whereas complexities are revealed in a small-scale,such as subparallel subfaults,unmapped branches,and stepovers.The seismicity on EAF is shallow(<15 km)and concentrated in depth distribution,indicating a clear lock-creep transition.In contrast,the SürgüFault(SF)that is responsible for the second M_(W)7.8 mainshock is shovel-shaped for the nucleation segment and has overall low dip angles(~40°–80°).Aftershocks on the SF distribute in a broad range of depth,extending down to~35 km.We also analyze the temporal behavior of seismicity,discovering no immediate foreshocks within~5 days preceding the first mainshock,and no seismic activity on the SF before the second mainshock.
基金supported by the National Natural Science Foundation of China (Nos. 41764005, 41604039, 41604102, and 41574078)Guangxi Natural Science Foundation of China (Nos. 2016GXNSFBA380082 and 2016GXNSFBA380215)+2 种基金Guangxi Young and Middle-aged Teacher Basic Ability Improvement Project (No. KY2016YB199)Guangxi Collaborative Innovation Center for Exploration of Hidden Nonferrous Metal Deposits and Development of New Materials Project (No. GXYSXTZX2017-II-5)Guangxi Scholarship Fund of Guangxi Education Department。
文摘Correctly locating the tunnel lining cavity is extremely important tunnel quality inspection.High-accuracy imaging results are hard to obtain because conventional one-way wave migration is greatly aff ected by lateral velocity change and inclination limitation and because the diff racted wave cannot be accurately returned to the real spatial position of the lining cavity.This paper presents a tunnel lining cavity imaging method based on the groundpenetrating radar(GPR)reverse-time migration(RTM)algorithm.The principle of GPR RTM is described in detail using the electromagnetic wave equation.The finite-difference timedomain method is employed to calculate the backward extrapolation electromagnetic fi elds,and the zero-time imaging condition based on the exploding-reflector concept is used to obtain the RTM results.On this basis,the GPR RTM program is compiled and applied to the simulated and observed GPR data of a typical tunnel lining cavity GPR model and a physical lining cavity model.Comparison of RTM and Kirchhoff migration results reveals that the RTM can better converge the diff racted waves of steel bar and cavity to their true position and have higher resolution and better suppress the eff ect of multiple interference and clutter scattering waves.In addition,comparison of RTM results of diff erent degrees of noise shows that RTM has strong anti-interference ability and can be used for the accurate interpretation of radar profi le in a strong interference environment.
基金The National Natural Science Foundation of China under contract Nos 41504084 and 41274120
文摘One of the problems experienced in marine geophysical exploration is that the layered features in the migration imaging profile are blurred and the seismic energy reflected is weaker in the middle or lower parts. In this study we model the seismic wavefield records in the undulating seafloor when there is both a slight change and significant change in seafloor topography to analyze its influence on the seismic reflection data and migration imaging profiles. We compare and analyze the wave field records collected at the same point on the original and modified velocity models, and the cross-bonding resulting migration imaging profiles. The results show that whether the seismic reflection data collection is performed along the direction of the survey line or against the direction of the survey line, slight changes in the seafloor topography have little effect on the wave field records and the migration profile, while significant changes in the seafloor topography have great effect on both the wave field records and migration profile.
基金sponsored by the Natural Science Foundation of China(No.41874139)the Natural Science Foundation of China(No.41674124)Jilin Province Foundation for Excellent Youths(No.20190103139JH)
文摘With the development of seismic exploration,passive-source seismic data has attracted increasing attention.Ambient noise passive seismic sources exists widely in nature and industrial production.Passive seismic data is important in logging while drilling(LWD),large-scale structural exploration,etc.In this paper,we proposed a passive multiple reverse time migration imaging(PMRTMI)method based on wavefield decomposition and normalized imaging conditions method.This method differs from seismic interferometry in that it can use raw passive seismic data directly in RTM imaging without reconstruction of virtual active gather,and we use the wavefield decomposition method to eliminate the low frequency noise in RTM.Further,the energy normalized imaging condition is used in full wavefield decomposition,which can not only enhance the image quality of both edge and deep information but also overcome the wrong energy problem caused by uneven distribution of passive sources;furthermore,this method exhibits high efficiency.Finally,numerical examples with the Marmousi model show the effectiveness of the method.
基金supported by National Natural Science Foundation of China,http://www.nsfc.gov.cn/.Peng H.received the project No.62071165.
文摘The high frame rate(HFR)imaging technique requires only one emission event for imaging.Therefore,it can achieve ultrafast imaging with frame rates up to the kHz regime,which satisfies the frame rate requirements for imaging moving tissues in scientific research and clinics.Lu’s Fourier migration method is based on a non-diffraction beam to obtain HFR images and can improve computational speed and efficiency.However,in order to obtain high-quality images,Fourier migration needs to make full use of the spectrum of echo signals for imaging,which requires a large number of Fast Fourier Transform(FFT)points and increases the complexity of the hardware when the echo frequency is high.Here,an efficient algorithm using the spectrum migration technique based on the spectrum’s distribution characteristics is proposed to improve the imaging efficiency in HFR imaging.Since the actual echo signal spectrum is of limited bandwidth,low-frequency and high-frequency parts with low-energy have little contribution to the imaging spectrum.We transform the effective part that provides the main energy in the signal spectrum to the imaging spectrum while the ineffective spectrum components are not utilized for imaging.This can significantly reduce the number of Fourier transform points,improve Fourier imaging efficiency,and ensure the imaging quality.The proposed method is evaluated on simulated and experimental datasets.Results demonstrated that the proposed method could achieve equivalent image quality with a reduced point number for FFT compared to the complete spectrum migration.In this paper,it only requires a quarter of the FFT points used in the complete spectrum migration,which can improve the computational efficiency;thus,it is more suitable for real-time data processing.The proposed spectrum migration method has a specific significance for the study and clinical application of HFR imaging.
文摘In view of the seismic exploration problem of thin sand reservoirs in the Songliao Basin, this paper puts forward a migration imaging method using CGP （common geophone point） stacked cylindrical waves. By this means, seismic data should be acquired from a midpoint shooting layout system with small shot-point spacing and small geophone interval. Using such seismic data, CGP gathers are first stacked to compose a cylindrical wave section. The cylindrical wave section is migrated and imaged by means of the ray path downward continuation of the down-going wave and the wave equation downward continuation of the upgoing wave. The results from the modeling analysis and the data processing of the TK8157 seismic line in the Songliao Basin shows that the proposed migration imaging method has higher seismic resolution and fidelity. Furthermore, the proposed method is proven to be more effective for discovering small sand bodies, small faults, stratigraphic pinch-outs, and so on.
文摘With the increasing complexity of prospecting objectives,reverse time migration( RTM) has attracted more and more attention due to its outstanding imaging quality. RTM is based on two-way wave equation,so it can avoid the limits of angle in traditional one-way wave equation migration,image reverse branch,prism waves and multi-reflected wave precisely and obtain accurate dynamic information. However,the huge demands for storage and computation as well as low frequency noises restrict its wide application. The normalized cross-correlation imaging conditions based on wave field decomposition are derived from traditional cross-correlation imaging condition,and it can eliminate the low-frequency noises effectively and improve the imaging resolution. The practical procedure includes separating source and receiver wave field into one-way components respectively,and conducting cross-correlation imaging condition to the post-separated wave field. In this way,the resolution and precision of the imaging result will be promoted greatly.
文摘Based on the study on electromagnetic field migration by Zhdanov, we have proposed an improved method for the weak points in the research. Firstly, the initial background resistivity should be determined by using 1-D inversion results. Then in the process of continuation, the results are corrected and calculated layer by layer by the iteration method, so that more exact resistivity can be obtained. Secondly, an improved algorithm for finite-difference equation is studied. According to the property of electromagnetic migration field, the algorithm is designed by means of grids varying with geometric progression in the longitudinal direction. Being improved by the techniques mentioned above, better results are obtained by the new method, which has been verified by both the theory model and practical data.
文摘In seismic exploration,it is a critical task to image and interpret different seismic signatures over complex geology due to strong lateral velocity contrast,steep reflectors,overburden strata and dipping flanks.To understand the behavior of these seismic signatures,nowadays Reverse Time Migration(RTM)technique is used extensively by the oil&gas industries.During the extrapolation phase of RTM,the source wavefield needs to be saved,which needs high storage memory and large computing time.These two are the main obstacles of RTM for production use.In order to overcome these disadvantages,in this study,a second-generation improved RTM technique is proposed.In this improved form,a shift operator is introduced at the time of imaging condition of RTM algorithm which is performed automatically both in space and time domain.This effort is made to produce a better-quality image by minimizing the computational time as well as numerical artefacts.The proposed method is applied over various benchmark models and validated by implementing over one field data set from the Jaisalmer Basin,India.From the analysis,it is observed that the method consumes a minimum of 45%less storage space and reduce the execution time by 20%,as compared to conventional RTM.The proposed RTM is found to work efficiently in comparison to the conventional RTM both in terms of imaging quality and minimization of numerical artefacts for all the benchmark models as well as field data.
基金This work is supported by the National Natural Science Foundation of China(No.41604039,41604102,41764005,41574078)Guangxi Natural Science Foundation project(No.2020GXNSFAA159121,2016GXNSFBA380215).
文摘The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal and produce a high precision subsurface image when the object is buried in a complicated subsurface environment.To better handle the multi-off set GPR data,we propose a three-dimensional(3D)prestack RTM algorithm.The high-order fi nite diff erence time domian(FDTD)method,with the accuracy of eighth-order in space and second-order in time,is applied to simulate the forward and backward extrapolation electromagnetic fi elds.In addition,we use the normalized correlation imaging condition to obtain pre-stack RTM result and the Laplace fi lter to suppress the low frequency noise generated during the correlation process.The numerical test of 3D simulated GPR data demonstrated that 3D RTM image shows excellent coincidence with the true model.Compared with 2D RTM image,the 3D RTM image can more clearly and accurately refl ect the 3D spatial distribution of the target,and the resolution of the imaging results is far better.Furthermore,the application of observed GPR data further validates the eff ectiveness of the proposed 3D GPR RTM algorithm,and its fi nal image can more reliably guide the subsequent interpretation.
文摘Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by which static correction is completely abandoned before migration and surely the imaging quality is remarkably improved. To obtain some artificial shot gathers, high-order staggered-grid finite-difference （FD） method is adapted to model acoustic wave propagation. Since the shot gathers are always disturbed by regular interferences, the statics still must be applied to supporting the interference elimination by apparent velocity filtering method. Then all the shot gathers should be removed back to their original positions by reverse statics. Finally, they are migrated by pre-stack reverse-time depth migration and imaged. The numerical experiments show that the DPDM can ideally avoid the mistakes caused by statics and increase imaging precision.
基金supported by the National Natural Science Foundation of China(Grant No.41574099)the National Key Science and Technology Special Projects(grant No.2016ZX05002004-005)
文摘Anisotropy correction is necessary during the processing of converted PS- wave seismic data to achieve accurate structural imaging, reservoir prediction, and fracture detection. To effectively eliminate the adverse effects of S-wave splitting and to improve PS- wave imaging quality, we tested methods for pre-stack migration imaging and anisotropic correction of PS-wave data. We based this on the propagation rules of seismic waves in a horizontal transverse isotropy medium, which is a fractured medium model that reflects likely subsurface conditions in the field. We used the radial （R） and transverse （T） components of PS-wave data to separate the fast and slow S-wave components, after which their propagation moveout was effectively extracted. Meanwhile, corrections for the energies and propagation moveouts of the R and T components were implemented using mathematical rotation. The PS-wave imaging quality was distinctly improved, and we demonstrated the reliability of our methods through numerical simulations. Applying our methods to three-dimensional and three-component seismic field data from the Xinchang-Hexingchang region of the Western Sichuan Depression in China, we obtained high-quality seismic imaging with continuous reflection wave groups, distinct structural features, and specific stratigraphic contact relationships. This study provides an effective and reliable approach for data processing that will improve the exploration of complex, hidden lithologic gas reservoirs.
基金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 simplified Gaussian beam migration,which omits the dynamic information and can calculate multi-arrival traveltime,so it is a high-precision and fast seismic imaging method.In the imaging process,extracting common image gathers can be used for velocity analysis,improving the accuracy of modeling and imaging quality.Compared with the conventional common image gathers extracting methods,the angle-domain common image gathers extracting method can avoid the artifacts caused by multi-arrival seismic waves.The authors present a new method of extracting common image gathers in angle-domain from Kirchhoff beam migration and verify the method by numerical calculations.
基金funded by the Young Scientists Fund of the National Natural Science Foundation of China(42304135)the scientific research project of Gansu Coal Geology Bureau(2023-07).
文摘Reasonable field acquisition geometry can not only guide seismic exploration to obtain sufficient geological information of target body,but also reduce acquisition cost to the maximum.In this study,building on conventional ray-based geometry design methods,we incorporate imaging results as a constraint to optimize the geometry design and evaluate its effectiveness.Firstly,the geological model of the target layer is established based on the geological data of the study area and surface seismic data combined with exploration tasks.Then,the ray-tracing method is employed to simulate and assess the proposed geometry design,verifying whether its parameters meet the exploration requirements.Finally,the imaging effect of the designed geometry on the target layer is tested by the cross-well seismic reverse time migration method.This methodology was applied to design the cross-well seismic acquisition geometry for offshore deviated wells in the X Oilfield.The simulation results demonstrate that the imaging-driven geometry design approach effectively guides field operations,enhances the imaging quality of the target layer,and reduces acquisition costs.
基金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.
