The main factors affecting seismic exploration is the propagation velocity of seismic waves in the medium. In the past, during marine seismic data processing, the propagation velocity of sea water was generally taken ...The main factors affecting seismic exploration is the propagation velocity of seismic waves in the medium. In the past, during marine seismic data processing, the propagation velocity of sea water was generally taken as a constant 1500 m/s. However, for deep water exploration, the sound velocity varies with the season, time, location, water depth, ocean currents, and etc.. It also results in a layered velocity distribution, so there is a difference of seismic traveltime, ray paths, and amplitude, which affect the migration imaging results if sea water propagation velocity is still taken as constant for the propagation wavefield. In this paper, we will start from an empirical equation of seismic wave velocity in seawater with changes of temperature, salinity, and depth, consider the variation of their values, build a seawater velocity model, and quantitatively analyze the impact of seawater velocity variation on seismic traveltime, ray paths, and amplitude in the seawater velocity model.展开更多
To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, t...To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.展开更多
Irregular surface flattening,which is based on a boundary conforming grid and the transformation between curvilinear and Cartesian coordinate systems,is a mathematical method that can elegantly handle irregular surfac...Irregular surface flattening,which is based on a boundary conforming grid and the transformation between curvilinear and Cartesian coordinate systems,is a mathematical method that can elegantly handle irregular surfaces,but has been limited to obtaining first arrivals only.By combining a multistage scheme with the fast-sweeping method(FSM,the method to obtain first-arrival traveltime in curvilinear coordinates),the reflected waves from a crustal interface can be traced in a topographic model,in which the reflected wavefront is obtained by reinitializing traveltimes in the interface for upwind branches.A local triangulation is applied to make a connection between velocity and interface nodes.Then a joint inversion of first-arrival and reflection traveltimes for imaging seismic velocity structures in complex terrains is presented.Numerical examples all perform well with different seismic velocity models.The increasing topographic complexity and even use of a high curvature reflector in these models demonstrate the reliability,accuracy and robustness of the new working scheme;checkerboard testing illustrates the method's high resolution.Noise tolerance testing indicates the method's ability to yield practical traveltime tomography.Further development of the multistage scheme will allow other later arrivals to be traced and used in the traveltime inversion.展开更多
A new method based on Huygens' principle and Fermat's principle is presented to calculate seismic traveltime and ray-paths in this paper. The algorithm can be easily understood and programmed, and can be appli...A new method based on Huygens' principle and Fermat's principle is presented to calculate seismic traveltime and ray-paths in this paper. The algorithm can be easily understood and programmed, and can be applied to heteroge neous media in order to eliminate its disadvantage of slowness, we have improved the basic algorithm to speed its calculation to practical stage without the loss of its accuracy.展开更多
Identifying potential hazards is crucial for maintaining the structural stability of opencast mining area.To address the limitations of irregular structure and sparse microseismic events in opencast mining monitoring,...Identifying potential hazards is crucial for maintaining the structural stability of opencast mining area.To address the limitations of irregular structure and sparse microseismic events in opencast mining monitoring,this paper proposes an active-source imaging method for identifying potential hazards precisely based on velocity structure.This method innovatively divides the irregular structure into unstructured grids and introduces a damping and smoothing regularization operator into the inversion process,mitigating the ill-posedness caused by the sparse distribution of events and rays.Numerical and laboratory experiments were conducted to verify the reliability and effectiveness of the proposed method.The results demonstrate the competitive performance of the method in identifying hazard areas of varying sizes and numbers.The proposed method shows potential for meeting hazard identification requirements in the complex opencast mining structure.Furthermore,field experiments were conducted on an rare earth mine slope.It confirms that the proposed method provides a more concrete and intuitive scheme for stability monitoring for the microseismic monitoring system.This paper not only demonstrates the application of acoustic structure velocity imaging technology in detecting unstructured potential hazard regions but also provides valuable insights into the construction and maintenance of stable opencast mining area.展开更多
3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic m...3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.展开更多
The picking efficiency of seismic first breaks(FBs)has been greatly accelerated by deep learning(DL)technology.However,the picking accuracy and efficiency of DL methods still face huge challenges in low signal-to-nois...The picking efficiency of seismic first breaks(FBs)has been greatly accelerated by deep learning(DL)technology.However,the picking accuracy and efficiency of DL methods still face huge challenges in low signal-to-noise ratio(SNR)situations.To address this issue,we propose a regression approach to pick FBs based on bidirectional long short-term memory(Bi LSTM)neural network by learning the implicit Eikonal equation of 3D inhomogeneous media with rugged topography in the target region.We employ a regressive model that represents the relationships among the elevation of shots,offset and the elevation of receivers with their seismic traveltime to predict the unknown FBs,from common-shot gathers with sparsely distributed traces.Different from image segmentation methods which automatically extract image features and classify FBs from seismic data,the proposed method can learn the inner relationship between field geometry and FBs.In addition,the predicted results by the regressive model are continuous values of FBs rather than the discrete ones of the binary distribution.The picking results of synthetic data shows that the proposed method has low dependence on label data,and can obtain reliable and similar predicted results using two types of label data with large differences.The picking results of9380 shots for 3D seismic data generated by vibroseis indicate that the proposed method can still accurately predict FBs in low SNR data.The subsequent stacked profiles further illustrate the reliability and effectiveness of the proposed method.The results of model data and field seismic data demonstrate that the proposed regression method is a robust first-break picker with high potential for field application.展开更多
3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be effi...3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.展开更多
The intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity mo...The intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity model used for the conventional intersection method is based on homogeneous or laterally homogeneous media, which is too simple. In order to improve the accuracy, we have modified the existing intersection method. In the modified approach, the earthquake loci are not assumed to be circular or hyperbolic and calculation accuracy is improved using a minimum traveltime tree algorithm for tracing rays. The numerical model shows that the modified method can locate earthquakes in complex velocity models.展开更多
Traditional AVO forward modeling only considers the impact of reflection coefficients at the interface on seismic wave field amplitude and ignores various propagation effects. Introducing wave propagation effects incl...Traditional AVO forward modeling only considers the impact of reflection coefficients at the interface on seismic wave field amplitude and ignores various propagation effects. Introducing wave propagation effects including geometric spreading, transmission loss, attenuation into seismic wave propagation, multi-wave amplitude-preserved AVO forward modeling for horizontally layered media based on ray theory is proposed in this paper. We derived the multi-wave geometric spreading correction formulas for horizontally layered media in order to describe the geometric spreading effect of multi-wave propagation. Introducing the complex traveltime directly, we built the relationship between complex traveltime and quality factor without the help of complex velocity to describe the attenuation of viscoelastic media. Multi-wave transmission coefficients, obtained by solving the Zoeppritz equations directly, is used to describe the transmission loss. Numerical results show that the effects of geometric spreading, attenuation, and transmission loss on multi-wave amplitude varies with offset and multi-wave amplitude-preserved AVO forward modeling should consider the reconstructive effect of wave propagation on reflection amplitude.展开更多
The conventional long-offset nonhyperbolic moveout equation is derived for the transverse isotropic media with a vertical symmetric axis（VTI）.It cannot be extended to the transverse isotropic media with an arbitrary...The conventional long-offset nonhyperbolic moveout equation is derived for the transverse isotropic media with a vertical symmetric axis（VTI）.It cannot be extended to the transverse isotropic media with an arbitrary spatial orientation of symmetry axis（ATI）.In this paper,we optimize a modified long-offset nonhyperbolic moveout equation for ATI media based on the conventional nonhyperbolic moveout equation and the exact analytical solution of the quartic moveout coefficient（A_4） and NMO velocity for ATI media that were derived in our previous work.Compared with the exact traveltimes of the ray-tracing algorithm for anisotropic media,this optimized equation can be used to calculate the traveltime varying with survey line azimuth in arbitrary strong or weak ATI media.It can replace the time-consuming, multi-offset,and multi-azimuth ray tracing method for forward modeling of long-offset reflection traveltimes in ATI media,which is useful to further anisotropic parameter inversion using long-offset nonhyperbolic moveout.展开更多
Reflection-based inversion that aims to reconstruct the low-to-intermediate wavenumbers of the subsurface model, can be a complementary to refraction-data-driven full-waveform inversion(FWI), especially for the deep t...Reflection-based inversion that aims to reconstruct the low-to-intermediate wavenumbers of the subsurface model, can be a complementary to refraction-data-driven full-waveform inversion(FWI), especially for the deep target area where diving waves cannot be acquired at the surface. Nevertheless, as a typical nonlinear inverse problem, reflection waveform inversion may easily suffer from the cycleskipping issue and have a slow convergence rate, if gradient-based first-order optimization methods are used. To improve the accuracy and convergence rate, we introduce the Hessian operator into reflection traveltime inversion(RTI) and reflection waveform inversion(RWI) in the framework of second-order optimization. A practical two-stage workflow is proposed to build the velocity model, in which Gauss-Newton RTI is first applied to mitigate the cycle-skipping problem and then Gauss-Newton RWI is employed to enhance the model resolution. To make the Gauss-Newton iterations more efficiently and robustly for large-scale applications, we introduce proper preconditioning for the Hessian matrix and design appropriate strategies to reduce the computational costs. The example of a real dataset from East China Sea demonstrates that the cascaded Hessian-based RTI and RWI have good potential to improve velocity model building and seismic imaging, especially for the deep targets.展开更多
Estimation of an accurate macro velocity model plays an important role in seismic imag- ing and model parameter inversion. Full waveform inversion （FWI） is the classical data-domain inver- sion method. However, the ...Estimation of an accurate macro velocity model plays an important role in seismic imag- ing and model parameter inversion. Full waveform inversion （FWI） is the classical data-domain inver- sion method. However, the misfit function of FWI is highly nonlinear, and the local optimization cannot prevent convergence of the misfit function toward local minima. To converge to the global minimum, FWI needs a good initial model or reliable low frequency component and long offset data. In this article, we present a wave-equation-based reflection traveltime tomography （WERTT） method, which can pro- vide a good background model （initial model） for FWI and （least-square） pre-stack depth migration （LS-PSDM）. First, the velocity model is decomposed into a low-wavenumber component （background velocity） and a high-wavenumber component （reflectivity）. Second, the primary reflection wave is pre- dicted by wave-equation demigration, and the reflection traveltime is calculated by an automatic picking method. Finally, the misfit function of the 12-norm of the reflection traveltime residuals is mini- mized by a gradient-based method. Numerical tests show that the proposed method can invert a good background model, which can be used as an initial model for LS-PSDM or FWI.展开更多
Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well descri...Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.展开更多
The Tilted tilted transversely isotropic(TTI)media,a kind of anisotropic medium,widely exists within the earth.For faster calculation of travel times in the TTI anisotropic media,we modifi ed a minimum traveltime tree...The Tilted tilted transversely isotropic(TTI)media,a kind of anisotropic medium,widely exists within the earth.For faster calculation of travel times in the TTI anisotropic media,we modifi ed a minimum traveltime tree algorithm with high effi ciency by dynamical modifi cation of the secondary wave propagation region during the spread of seismic waves.To manage the wavefront points in the modified version,we used a novel minimum heap sorting technique to reduce the time spent on selecting secondary waves points.In this study,seismic group velocities were obtained from analytical solutions in terms of phase angle,and the corresponding phase angles were determined by binary search rather than approximate equations for weakly anisotropic media.For the most time-consuming part of the secondary wave traveltime calculation,the parallel computation was initially performed using multiple cores and threads.Numerical examples showed that the improved method can calculate seismic travel times and ray paths faster and accurately in a 3D TTI medium.For four cores and eight threads,the computing speed increased by six times when compared to the conventional method.展开更多
2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcani...2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcanic region. It is used in this study that seismic traveltime inversion for simultaneous determination of 2-D velocity and interface structure of the crust and upper mantle. The result shows that, under Changbaishan-Tianchi crater, there exists a low-velocity body in the shape of an inverted triangle, and the crustal reflecting boundaries and Moho all become lower by a varying margin of 2-6 km, forming a crustal root which is assumed to be the Changbaishan-Tianchi volcanic system. Finally, we make a comparison between our 2-D velocity model and the result from the studies by using trial-and-error forward modeling with SEIS83.展开更多
The objective function of full waveform inversion is a strong nonlinear function,the inversion process is not unique,and it is easy to fall into local minimum.Firstly,in the process of wavefield reconstruction,the wav...The objective function of full waveform inversion is a strong nonlinear function,the inversion process is not unique,and it is easy to fall into local minimum.Firstly,in the process of wavefield reconstruction,the wave equation is introduced into the construction of objective function as a penalty term to broaden the search space of solution and reduce the risk of falling into local minimum.In addition,there is no need to calculate the adjoint wavefield in the inversion process,which can significantly improve the calculation efficiency;Secondly,considering that the total variation constraint can effectively reconstruct the discontinuous interface in the velocity model,this paper introduces the weak total variation constraint to avoid the excessive smooth estimation of the model under the strong total variation constraint.The disadvantage of this strategy is that it is highly dependent on the initial model.In view of this,this paper takes the long wavelength initial model obtained by first arrival traveltime tomography as a prior model constraint,and proposes a weak total variation constrained wavefield reconstruction inversion method based on first arrival traveltime tomography.Numerical experimental results show that the new method reduces the dependence on the initial model,the interface description is more accurate,the error is reduced,and the iterative convergence efficiency is significantly improved.展开更多
The Northeastern Tibetan plateau records Caledonian Qilian orogeny and Cenozoic reactivation by continental collision between the Indian and Asian plates. In order to provide the constraint on the Qilian orogenic mech...The Northeastern Tibetan plateau records Caledonian Qilian orogeny and Cenozoic reactivation by continental collision between the Indian and Asian plates. In order to provide the constraint on the Qilian orogenic mechanism and the expansion of the plateau,wide-angle seismic data was acquired along a 430 km-long profile between Jingtai and Hezuo. There is strong height variation along the profile,which is dealt by topography flattening scheme in our crustal velocity structure reconstruction. We herein present the upper crustal P-wave velocity structure model resulting from the interpretation of first arrival dataset from topography-dependent eikonal traveltime tomography. With topography flattening scheme to process real topography along the profile,the evenness of ray coverage times of the image area（upper crust）is improved,which provides upper crustal velocity model comparable to the classic traveltime tomography（with model expansion scheme to process irregular surface）. The upper crustal velocity model shows zoning character which matcheswith the tectonic division of the Qaidam-Kunlun-West Qinling belt,the Central and Northern Qilian,and the Alax blocks along the profile. The resultant upper crustal P-wave velocity model is expected to provide important base for linkage between the mapped surface geology and deep structure or geodynamics in Northeastern Tibet.展开更多
Marchenko imaging obtains the subsurface reflectors using one-way Green’s functions,which are retrieved by solving the Marchenko equation.This method generates an image that is free of spurious artifacts due to inter...Marchenko imaging obtains the subsurface reflectors using one-way Green’s functions,which are retrieved by solving the Marchenko equation.This method generates an image that is free of spurious artifacts due to internal multiples.The Marchenko imaging method is a target-oriented technique;thus,it can image a user specified area.In the traditional Marchenko method,an accurate velocity model is critical for estimating direct waves from imaging points to the surface.An error in the velocity model results in the inaccurate estimation of direct waves.In turn,this leads to errors in computation of one-way Green’s functions,which then affects the final Marchenko images.To solve this problem,in this paper,we propose a self-adaptive traveltime updating technique based on the principle of equal traveltime to improve the Marchenko imaging method.The proposed method calculates the time shift of direct waves caused by the error in the velocity model,and corrects the wrong direct wave according to the time shift and reconstructs the correct Green’s functions.The proposed method improves the results of imaging using an inaccurate velocity model.By comparing the results from traditional Marchenko and the new method using synthetic data experiments,we demonstrated that the adaptive traveltime updating Marchenko imaging method could restore the image of geological structures to their true positions.展开更多
基金supported by the Major Projects of National Science and Technology Sub-topics(2011ZX05025-001-05)
文摘The main factors affecting seismic exploration is the propagation velocity of seismic waves in the medium. In the past, during marine seismic data processing, the propagation velocity of sea water was generally taken as a constant 1500 m/s. However, for deep water exploration, the sound velocity varies with the season, time, location, water depth, ocean currents, and etc.. It also results in a layered velocity distribution, so there is a difference of seismic traveltime, ray paths, and amplitude, which affect the migration imaging results if sea water propagation velocity is still taken as constant for the propagation wavefield. In this paper, we will start from an empirical equation of seismic wave velocity in seawater with changes of temperature, salinity, and depth, consider the variation of their values, build a seawater velocity model, and quantitatively analyze the impact of seawater velocity variation on seismic traveltime, ray paths, and amplitude in the seawater velocity model.
基金National Natural Science Foundation of China (49894190-024) and Geophysical Prospecting Key Laboratory Foundation of China National Petroleum Corporation.
文摘To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.
基金financial support for this work contributed by the National Key Research and Development Program of China(grant nos.2016YFC0600302,2016YFC0600101 and 2016YFC0600201)the National Natural Science Foundation of China(grants 41604075,41430213,41574092 and 41474068)
文摘Irregular surface flattening,which is based on a boundary conforming grid and the transformation between curvilinear and Cartesian coordinate systems,is a mathematical method that can elegantly handle irregular surfaces,but has been limited to obtaining first arrivals only.By combining a multistage scheme with the fast-sweeping method(FSM,the method to obtain first-arrival traveltime in curvilinear coordinates),the reflected waves from a crustal interface can be traced in a topographic model,in which the reflected wavefront is obtained by reinitializing traveltimes in the interface for upwind branches.A local triangulation is applied to make a connection between velocity and interface nodes.Then a joint inversion of first-arrival and reflection traveltimes for imaging seismic velocity structures in complex terrains is presented.Numerical examples all perform well with different seismic velocity models.The increasing topographic complexity and even use of a high curvature reflector in these models demonstrate the reliability,accuracy and robustness of the new working scheme;checkerboard testing illustrates the method's high resolution.Noise tolerance testing indicates the method's ability to yield practical traveltime tomography.Further development of the multistage scheme will allow other later arrivals to be traced and used in the traveltime inversion.
文摘A new method based on Huygens' principle and Fermat's principle is presented to calculate seismic traveltime and ray-paths in this paper. The algorithm can be easily understood and programmed, and can be applied to heteroge neous media in order to eliminate its disadvantage of slowness, we have improved the basic algorithm to speed its calculation to practical stage without the loss of its accuracy.
基金Project(2021YFC2900500)supported by the National Key Research and Development Program of China。
文摘Identifying potential hazards is crucial for maintaining the structural stability of opencast mining area.To address the limitations of irregular structure and sparse microseismic events in opencast mining monitoring,this paper proposes an active-source imaging method for identifying potential hazards precisely based on velocity structure.This method innovatively divides the irregular structure into unstructured grids and introduces a damping and smoothing regularization operator into the inversion process,mitigating the ill-posedness caused by the sparse distribution of events and rays.Numerical and laboratory experiments were conducted to verify the reliability and effectiveness of the proposed method.The results demonstrate the competitive performance of the method in identifying hazard areas of varying sizes and numbers.The proposed method shows potential for meeting hazard identification requirements in the complex opencast mining structure.Furthermore,field experiments were conducted on an rare earth mine slope.It confirms that the proposed method provides a more concrete and intuitive scheme for stability monitoring for the microseismic monitoring system.This paper not only demonstrates the application of acoustic structure velocity imaging technology in detecting unstructured potential hazard regions but also provides valuable insights into the construction and maintenance of stable opencast mining area.
基金The authors thank the funds supported by the China National Nuclear Corporation under Grants Nos.WUQNYC2101 and WUHTLM2101-04National Natural Science Foundation of China(42074132,42274154).
文摘3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.
基金financially supported by the National Key R&D Program of China(2018YFA0702504)the National Natural Science Foundation of China(42174152)+1 种基金the Strategic Cooperation Technology Projects of China National Petroleum Corporation(CNPC)and China University of Petroleum-Beijing(CUPB)(ZLZX2020-03)the R&D Department of China National Petroleum Corporation(2022DQ0604-01)。
文摘The picking efficiency of seismic first breaks(FBs)has been greatly accelerated by deep learning(DL)technology.However,the picking accuracy and efficiency of DL methods still face huge challenges in low signal-to-noise ratio(SNR)situations.To address this issue,we propose a regression approach to pick FBs based on bidirectional long short-term memory(Bi LSTM)neural network by learning the implicit Eikonal equation of 3D inhomogeneous media with rugged topography in the target region.We employ a regressive model that represents the relationships among the elevation of shots,offset and the elevation of receivers with their seismic traveltime to predict the unknown FBs,from common-shot gathers with sparsely distributed traces.Different from image segmentation methods which automatically extract image features and classify FBs from seismic data,the proposed method can learn the inner relationship between field geometry and FBs.In addition,the predicted results by the regressive model are continuous values of FBs rather than the discrete ones of the binary distribution.The picking results of synthetic data shows that the proposed method has low dependence on label data,and can obtain reliable and similar predicted results using two types of label data with large differences.The picking results of9380 shots for 3D seismic data generated by vibroseis indicate that the proposed method can still accurately predict FBs in low SNR data.The subsequent stacked profiles further illustrate the reliability and effectiveness of the proposed method.The results of model data and field seismic data demonstrate that the proposed regression method is a robust first-break picker with high potential for field application.
基金supported by NSFC(Nos.41274120,41404085,and 41504084)
文摘3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.
基金This work is supported by the National Natural Science Foundation of China(40674044)the Special Foundation for Basic Professional Scientific Research (DQJB06A02)
文摘The intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity model used for the conventional intersection method is based on homogeneous or laterally homogeneous media, which is too simple. In order to improve the accuracy, we have modified the existing intersection method. In the modified approach, the earthquake loci are not assumed to be circular or hyperbolic and calculation accuracy is improved using a minimum traveltime tree algorithm for tracing rays. The numerical model shows that the modified method can locate earthquakes in complex velocity models.
基金sponsored by the National Natural Science Foundation of China (Grant No. 41074098)the National Basic Research Program of China (973 Program) (Grant No. 2007CB209606)
文摘Traditional AVO forward modeling only considers the impact of reflection coefficients at the interface on seismic wave field amplitude and ignores various propagation effects. Introducing wave propagation effects including geometric spreading, transmission loss, attenuation into seismic wave propagation, multi-wave amplitude-preserved AVO forward modeling for horizontally layered media based on ray theory is proposed in this paper. We derived the multi-wave geometric spreading correction formulas for horizontally layered media in order to describe the geometric spreading effect of multi-wave propagation. Introducing the complex traveltime directly, we built the relationship between complex traveltime and quality factor without the help of complex velocity to describe the attenuation of viscoelastic media. Multi-wave transmission coefficients, obtained by solving the Zoeppritz equations directly, is used to describe the transmission loss. Numerical results show that the effects of geometric spreading, attenuation, and transmission loss on multi-wave amplitude varies with offset and multi-wave amplitude-preserved AVO forward modeling should consider the reconstructive effect of wave propagation on reflection amplitude.
基金the National Natural Science Foundation of China(Grant No.40874028)the Special Fund (Grant No.2008ZX05008-006-004).
文摘The conventional long-offset nonhyperbolic moveout equation is derived for the transverse isotropic media with a vertical symmetric axis（VTI）.It cannot be extended to the transverse isotropic media with an arbitrary spatial orientation of symmetry axis（ATI）.In this paper,we optimize a modified long-offset nonhyperbolic moveout equation for ATI media based on the conventional nonhyperbolic moveout equation and the exact analytical solution of the quartic moveout coefficient（A_4） and NMO velocity for ATI media that were derived in our previous work.Compared with the exact traveltimes of the ray-tracing algorithm for anisotropic media,this optimized equation can be used to calculate the traveltime varying with survey line azimuth in arbitrary strong or weak ATI media.It can replace the time-consuming, multi-offset,and multi-azimuth ray tracing method for forward modeling of long-offset reflection traveltimes in ATI media,which is useful to further anisotropic parameter inversion using long-offset nonhyperbolic moveout.
基金supported by National Natural Science Foundation of China (42074157)the National Key Research and Development Program of China (2018YFC0310104)the Strategic Priority Research Program of the Chinese Academy of Science(XDA14010203)。
文摘Reflection-based inversion that aims to reconstruct the low-to-intermediate wavenumbers of the subsurface model, can be a complementary to refraction-data-driven full-waveform inversion(FWI), especially for the deep target area where diving waves cannot be acquired at the surface. Nevertheless, as a typical nonlinear inverse problem, reflection waveform inversion may easily suffer from the cycleskipping issue and have a slow convergence rate, if gradient-based first-order optimization methods are used. To improve the accuracy and convergence rate, we introduce the Hessian operator into reflection traveltime inversion(RTI) and reflection waveform inversion(RWI) in the framework of second-order optimization. A practical two-stage workflow is proposed to build the velocity model, in which Gauss-Newton RTI is first applied to mitigate the cycle-skipping problem and then Gauss-Newton RWI is employed to enhance the model resolution. To make the Gauss-Newton iterations more efficiently and robustly for large-scale applications, we introduce proper preconditioning for the Hessian matrix and design appropriate strategies to reduce the computational costs. The example of a real dataset from East China Sea demonstrates that the cascaded Hessian-based RTI and RWI have good potential to improve velocity model building and seismic imaging, especially for the deep targets.
基金the financial support from National Key Research and Development Program of China(2021YFC2900500)Funds for International Cooperation and Exchange of the National Natural Science Foundation of China(52161135301).
基金financially supported by the National Natural Science Foundation of China(No.41374117)the‘973’Project of China(No.2011 CB201002)the Great and Special Projects of China(Nos.2011ZX05003-003,2011ZX05005-005-008HZ,and 2011ZX05006-002)
文摘Estimation of an accurate macro velocity model plays an important role in seismic imag- ing and model parameter inversion. Full waveform inversion （FWI） is the classical data-domain inver- sion method. However, the misfit function of FWI is highly nonlinear, and the local optimization cannot prevent convergence of the misfit function toward local minima. To converge to the global minimum, FWI needs a good initial model or reliable low frequency component and long offset data. In this article, we present a wave-equation-based reflection traveltime tomography （WERTT） method, which can pro- vide a good background model （initial model） for FWI and （least-square） pre-stack depth migration （LS-PSDM）. First, the velocity model is decomposed into a low-wavenumber component （background velocity） and a high-wavenumber component （reflectivity）. Second, the primary reflection wave is pre- dicted by wave-equation demigration, and the reflection traveltime is calculated by an automatic picking method. Finally, the misfit function of the 12-norm of the reflection traveltime residuals is mini- mized by a gradient-based method. Numerical tests show that the proposed method can invert a good background model, which can be used as an initial model for LS-PSDM or FWI.
基金supported financially by the Ministry of Science and Technology of China(2011CB808904)the National Natural Science Foundation of China(Nos.41021063,41174075,41004034,41174043,and 41274090)
文摘Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.
基金funded by the National Key R&D Program of China (No. 2020YFA0710600)National Science Foundation of China (No. 41374098)the Special Fund of the Institute of Geophysics,China Earthquake Administration (No. DQJB19B40)
文摘The Tilted tilted transversely isotropic(TTI)media,a kind of anisotropic medium,widely exists within the earth.For faster calculation of travel times in the TTI anisotropic media,we modifi ed a minimum traveltime tree algorithm with high effi ciency by dynamical modifi cation of the secondary wave propagation region during the spread of seismic waves.To manage the wavefront points in the modified version,we used a novel minimum heap sorting technique to reduce the time spent on selecting secondary waves points.In this study,seismic group velocities were obtained from analytical solutions in terms of phase angle,and the corresponding phase angles were determined by binary search rather than approximate equations for weakly anisotropic media.For the most time-consuming part of the secondary wave traveltime calculation,the parallel computation was initially performed using multiple cores and threads.Numerical examples showed that the improved method can calculate seismic travel times and ray paths faster and accurately in a 3D TTI medium.For four cores and eight threads,the computing speed increased by six times when compared to the conventional method.
基金Key Project (95-11-02-01) from China Seismological Bureau.Contribution No. RCEG200129, Research Center of Exploration Geophysi
文摘2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcanic region. It is used in this study that seismic traveltime inversion for simultaneous determination of 2-D velocity and interface structure of the crust and upper mantle. The result shows that, under Changbaishan-Tianchi crater, there exists a low-velocity body in the shape of an inverted triangle, and the crustal reflecting boundaries and Moho all become lower by a varying margin of 2-6 km, forming a crustal root which is assumed to be the Changbaishan-Tianchi volcanic system. Finally, we make a comparison between our 2-D velocity model and the result from the studies by using trial-and-error forward modeling with SEIS83.
基金supported by National Key R&D Program of China under contract number 2019YFC0605503CThe Major projects of CNPC under contract number(ZD2019-183-003)+2 种基金the Major projects during the 14th Five-year Plan period under contract number 2021QNLM020001the National Outstanding Youth Science Foundation under contract number 41922028the Funds for Creative Research Groups of China under contract number 41821002.
文摘The objective function of full waveform inversion is a strong nonlinear function,the inversion process is not unique,and it is easy to fall into local minimum.Firstly,in the process of wavefield reconstruction,the wave equation is introduced into the construction of objective function as a penalty term to broaden the search space of solution and reduce the risk of falling into local minimum.In addition,there is no need to calculate the adjoint wavefield in the inversion process,which can significantly improve the calculation efficiency;Secondly,considering that the total variation constraint can effectively reconstruct the discontinuous interface in the velocity model,this paper introduces the weak total variation constraint to avoid the excessive smooth estimation of the model under the strong total variation constraint.The disadvantage of this strategy is that it is highly dependent on the initial model.In view of this,this paper takes the long wavelength initial model obtained by first arrival traveltime tomography as a prior model constraint,and proposes a weak total variation constrained wavefield reconstruction inversion method based on first arrival traveltime tomography.Numerical experimental results show that the new method reduces the dependence on the initial model,the interface description is more accurate,the error is reduced,and the iterative convergence efficiency is significantly improved.
基金financial support for this work by the Ministry of Science and Technology of China (2011CB808904)the Ministry of Land and Resources of China (SinoProbe-02-02 or 201011041,SinoProbe-03-02 or 201011047)the National Nature Science Foundation of China (41174075,41021063,41274090 and 41174043)
文摘The Northeastern Tibetan plateau records Caledonian Qilian orogeny and Cenozoic reactivation by continental collision between the Indian and Asian plates. In order to provide the constraint on the Qilian orogenic mechanism and the expansion of the plateau,wide-angle seismic data was acquired along a 430 km-long profile between Jingtai and Hezuo. There is strong height variation along the profile,which is dealt by topography flattening scheme in our crustal velocity structure reconstruction. We herein present the upper crustal P-wave velocity structure model resulting from the interpretation of first arrival dataset from topography-dependent eikonal traveltime tomography. With topography flattening scheme to process real topography along the profile,the evenness of ray coverage times of the image area（upper crust）is improved,which provides upper crustal velocity model comparable to the classic traveltime tomography（with model expansion scheme to process irregular surface）. The upper crustal velocity model shows zoning character which matcheswith the tectonic division of the Qaidam-Kunlun-West Qinling belt,the Central and Northern Qilian,and the Alax blocks along the profile. The resultant upper crustal P-wave velocity model is expected to provide important base for linkage between the mapped surface geology and deep structure or geodynamics in Northeastern Tibet.
基金supported by the National Natural Science Foundation of China(No.41874167)the National Science and Technology Major Project of China(No.2017YFB0202904)the National Natural Science Foundation of China(No.41904130)。
文摘Marchenko imaging obtains the subsurface reflectors using one-way Green’s functions,which are retrieved by solving the Marchenko equation.This method generates an image that is free of spurious artifacts due to internal multiples.The Marchenko imaging method is a target-oriented technique;thus,it can image a user specified area.In the traditional Marchenko method,an accurate velocity model is critical for estimating direct waves from imaging points to the surface.An error in the velocity model results in the inaccurate estimation of direct waves.In turn,this leads to errors in computation of one-way Green’s functions,which then affects the final Marchenko images.To solve this problem,in this paper,we propose a self-adaptive traveltime updating technique based on the principle of equal traveltime to improve the Marchenko imaging method.The proposed method calculates the time shift of direct waves caused by the error in the velocity model,and corrects the wrong direct wave according to the time shift and reconstructs the correct Green’s functions.The proposed method improves the results of imaging using an inaccurate velocity model.By comparing the results from traditional Marchenko and the new method using synthetic data experiments,we demonstrated that the adaptive traveltime updating Marchenko imaging method could restore the image of geological structures to their true positions.
