We present a finite difference （FD） method for the simulation of seismic wave fields in fractured medium with an irregular （non-fiat） free surface which is beneficial for interpreting exploration data acquired in ...We present a finite difference （FD） method for the simulation of seismic wave fields in fractured medium with an irregular （non-fiat） free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations （in a curvilinear coordinate system） in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry （HTI）. Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures.展开更多
The mapping method is a forward-modeling method that transforms the irregular surface to horizontal by mapping the rectangular grid as curved; moreover, the wave field calculations move from the physical domain to the...The mapping method is a forward-modeling method that transforms the irregular surface to horizontal by mapping the rectangular grid as curved; moreover, the wave field calculations move from the physical domain to the calculation domain. The mapping method deals with the irregular surface and the low-velocity layer underneath it using a fine grid. For the deeper high-velocity layers, the use of a fine grid causes local oversampling. In addition, when the irregular surface is transformed to horizontal, the flattened interface below the surface is transformed to curved, which produces inaccurate modeling results because of the presence of ladder-like burrs in the simulated seismic wave. Thus, we propose the mapping method based on the dual-variable finite-difference staggered grid. The proposed method uses different size grid spacings in different regions and locally variable time steps to match the size variability of grid spacings. Numerical examples suggest that the proposed method requires less memory storage capacity and improves the computational efficiency compared with forward modeling methods based on the conventional grid.展开更多
Presently the research based on the accurate seismic imaging methods for surface relief, complex structure, and complicated velocity distributions is of great significance. Reverse-time migration is considered to be o...Presently the research based on the accurate seismic imaging methods for surface relief, complex structure, and complicated velocity distributions is of great significance. Reverse-time migration is considered to be one of highly accurate methods. In this paper, we propose a new non-reflecting recursive algorithm for reverse-time migration by introducing the wave impedance function into the acoustic wave equation and the algorithm for the surface relief case is derived from the coordinate transformation principle. Using the exploding reflector principle and the zero-time imaging condition of poststack reverse- time migration, poststack numerical simulation and reverse-time migration with complex conditions can be realized. The results of synthetic and real data calculations show that the method effectively suppresses unwanted internal reflections and also deals with the seismic imaging problems resulting from surface relief. So, we prove that this method has strong adaptability and practicality.展开更多
On the basis of analyzing and researching the current algorithms of cycle-slip detection and correction, a new method of cycle-slip detection and correction is put forward in this paper, that is, a reasonable cycle-sl...On the basis of analyzing and researching the current algorithms of cycle-slip detection and correction, a new method of cycle-slip detection and correction is put forward in this paper, that is, a reasonable cycle-slip detection condition and algorithm with corresponding program COMPRE (COMpass PRE-processing) to detect and correct cycle-slip automatically, compared with GIPSY and GAMIT software, for example, it is proved that this method is effective and credible to cycle-slip detection and correction in GPS data pre-processing.展开更多
The IMU（inertial measurement unit） error equations in the earth fixed coordinates are introduced firstly. A fading Kalman filtering is simply introduced and its shortcomings are analyzed, then an adaptive filtering ...The IMU（inertial measurement unit） error equations in the earth fixed coordinates are introduced firstly. A fading Kalman filtering is simply introduced and its shortcomings are analyzed, then an adaptive filtering is applied in IMU/GPS integrated navigation system, in which the adaptive factor is replaced by the fading factor. A practical example is given. The resuits prove that the adaptive filter combined with the fading factor is valid and reliable when applied in IMU/GPS integrated navigation system.展开更多
The methodology of catchment extraction especially from regular grid digital elevation models (DEMs) is briefly reviewed. Then an efficient algorithm, which combines vector process and traditional neighbourhood raster...The methodology of catchment extraction especially from regular grid digital elevation models (DEMs) is briefly reviewed. Then an efficient algorithm, which combines vector process and traditional neighbourhood raster process, is designed for extracting the catchments and subcatchments from depressionless DEMs. The catchment area of each river in the grid DEM data is identified and delineated, then is divided into subcatchments as required. Compared to traditional processes, this method for identifying catchments focuses on the boundaries instead of the area inside the catchments and avoids the boundary intersection phenomena. Last, the algorithm is tested with a set of DEMs of different sizes, and the result proves that the computation efficiency and accuracy are better than existent methods.展开更多
A terrain-following coordinate （a-coordinate） in which the computational form of pressure gradient force （PGF） is two-term （the so-called classic method） has significant PGF errors near steep terrain. Using the ...A terrain-following coordinate （a-coordinate） in which the computational form of pressure gradient force （PGF） is two-term （the so-called classic method） has significant PGF errors near steep terrain. Using the covariant equations of the a-coordinate to create a one-term PGF （the covariant method） can reduce the PGF errors. This study investigates the factors inducing the PGF errors of these two methods, through geometric analysis and idealized experiments. The geometric analysis first demonstrates that the terrain slope and the vertical pressure gradient can induce the PGF errors of the classic method, and then generalize the effect of the terrain slope to the effect of the slope of each vertical layer （φ）. More importantly, a new factor, the direction of PGF （a）, is proposed by the geometric analysis, and the effects of φ and a are quantified by tan φ.tan a. When tan φ.tan a is greater than 1/9 or smaller than -10/9, the two terms of PGF of the classic method are of the same order but opposite in sign, and then the PGF errors of the classic method are large. Finally, the effects of three factors on inducing the PGF errors of the classic method are validated by a series of idealized experiments using various terrain types and pressure fields. The experimental results also demonstrate that the PGF errors of the covariant method are affected little by the three factors.展开更多
A new approach is proposed to use the covariant scalar equations of the a-coordinate （the covariant method）, in which the pressure gradient force （PGF） has only one term in each horizontal momentum equation, and t...A new approach is proposed to use the covariant scalar equations of the a-coordinate （the covariant method）, in which the pressure gradient force （PGF） has only one term in each horizontal momentum equation, and the PGF errors are much reduced in the computational space. In addition, the validity of reducing the PGF errors by this covariant method in the computational and physical space over steep terrain is investigated. First, the authors implement a set of idealized experiments of increasing terrain slope to compare the PGF errors of the covariant method and those of the classic method in the computational space. The results demonstrate that the PGF errors of the covariant method are consistently much-reduced, compared to those of the classic method. More importantly, the steeper the terrain, the greater the reduction in the ratio of the PGF errors via the covariant method. Next, the authors use geometric analysis to further investigate the PGF errors in the physical space, and the results illustrates that the PGF of the covariant method equals that of the classic method in the physical space; namely, the covariant method based on the non-orthogonal a-coordinate cannot reduce the PGF errors in the physical space. However, an orthogonal method can reduce the PGF errors in the physical space. Finally, a set of idealized experiments are carried out to validate the results obtained by the geometric analysis. These results indicate that the covariant method may improve the simulation of variables relevant to pressure, in addition to pressure itself, near steep terrain.展开更多
Today, the GNSS （global navigation satellite system） is used for more complicate and accurate applications such as monitoring or stake out works. The truth lies in the fact that in the most of the times not enough a...Today, the GNSS （global navigation satellite system） is used for more complicate and accurate applications such as monitoring or stake out works. The truth lies in the fact that in the most of the times not enough attention is paid to the antenna＇s setup. Usually, gross errors are found in the antenna＇s centering, leveling and in the measurement of its height, which are significant. In this paper, a thoroughly analysis of the above mentioned errors is carried out. The influence of these errors in the calculation of the X, Y, Z Cartesian geocentric coordinates and the ~, 2, h ellipsoid geodetic coordinates of a point P on the earth＇s surface, is analyzed and is presented in several diagrams. Also a new convenient method for the accurate measurement of the antenna＇s height is presented and it is strongly proposed. The conclusions outline the magnitude of these errors and prove the significance of the antenna＇s proper setup at the accurate GNSS applications.展开更多
基金sponsored by the Knowledge Innovation Program of the Chinese Academy of Sciences No.KZCX2-YW-132)the Important National Science and Technology Specific Projects(No.2008ZX05008-006)the National Natural Science Foundation of China Nos.41074033,40721003,40830315,and 40874041)
文摘We present a finite difference （FD） method for the simulation of seismic wave fields in fractured medium with an irregular （non-fiat） free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations （in a curvilinear coordinate system） in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry （HTI）. Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures.
基金financially supported by the National Natural Science Foundation of China(Nos.41104069 and 41274124)the National 973 Project(Nos.2014CB239006 and 2011CB202402)+1 种基金the Shandong Natural Science Foundation of China(No.ZR2011DQ016)Fundamental Research Funds for Central Universities(No.R1401005A)
文摘The mapping method is a forward-modeling method that transforms the irregular surface to horizontal by mapping the rectangular grid as curved; moreover, the wave field calculations move from the physical domain to the calculation domain. The mapping method deals with the irregular surface and the low-velocity layer underneath it using a fine grid. For the deeper high-velocity layers, the use of a fine grid causes local oversampling. In addition, when the irregular surface is transformed to horizontal, the flattened interface below the surface is transformed to curved, which produces inaccurate modeling results because of the presence of ladder-like burrs in the simulated seismic wave. Thus, we propose the mapping method based on the dual-variable finite-difference staggered grid. The proposed method uses different size grid spacings in different regions and locally variable time steps to match the size variability of grid spacings. Numerical examples suggest that the proposed method requires less memory storage capacity and improves the computational efficiency compared with forward modeling methods based on the conventional grid.
基金supported by the National Natural Science Foundation of China (Grant No. 40974073)the National 863 Program (Grant No.2007AA060504)the National 973 Program (Grant No. 2007CB209605) and CNPC Geophysical Laboratories
文摘Presently the research based on the accurate seismic imaging methods for surface relief, complex structure, and complicated velocity distributions is of great significance. Reverse-time migration is considered to be one of highly accurate methods. In this paper, we propose a new non-reflecting recursive algorithm for reverse-time migration by introducing the wave impedance function into the acoustic wave equation and the algorithm for the surface relief case is derived from the coordinate transformation principle. Using the exploding reflector principle and the zero-time imaging condition of poststack reverse- time migration, poststack numerical simulation and reverse-time migration with complex conditions can be realized. The results of synthetic and real data calculations show that the method effectively suppresses unwanted internal reflections and also deals with the seismic imaging problems resulting from surface relief. So, we prove that this method has strong adaptability and practicality.
文摘On the basis of analyzing and researching the current algorithms of cycle-slip detection and correction, a new method of cycle-slip detection and correction is put forward in this paper, that is, a reasonable cycle-slip detection condition and algorithm with corresponding program COMPRE (COMpass PRE-processing) to detect and correct cycle-slip automatically, compared with GIPSY and GAMIT software, for example, it is proved that this method is effective and credible to cycle-slip detection and correction in GPS data pre-processing.
基金Supported by the National Natural Science Foundation of China (No.40274002 No.40474001).
文摘The IMU（inertial measurement unit） error equations in the earth fixed coordinates are introduced firstly. A fading Kalman filtering is simply introduced and its shortcomings are analyzed, then an adaptive filtering is applied in IMU/GPS integrated navigation system, in which the adaptive factor is replaced by the fading factor. A practical example is given. The resuits prove that the adaptive filter combined with the fading factor is valid and reliable when applied in IMU/GPS integrated navigation system.
文摘The methodology of catchment extraction especially from regular grid digital elevation models (DEMs) is briefly reviewed. Then an efficient algorithm, which combines vector process and traditional neighbourhood raster process, is designed for extracting the catchments and subcatchments from depressionless DEMs. The catchment area of each river in the grid DEM data is identified and delineated, then is divided into subcatchments as required. Compared to traditional processes, this method for identifying catchments focuses on the boundaries instead of the area inside the catchments and avoids the boundary intersection phenomena. Last, the algorithm is tested with a set of DEMs of different sizes, and the result proves that the computation efficiency and accuracy are better than existent methods.
基金jointly supported by the National Basic Research Program of China[973 Program,grant number 2015CB954102]National Natural Science Foundation of China[grant numbers41305095 and 41175064]
文摘A terrain-following coordinate （a-coordinate） in which the computational form of pressure gradient force （PGF） is two-term （the so-called classic method） has significant PGF errors near steep terrain. Using the covariant equations of the a-coordinate to create a one-term PGF （the covariant method） can reduce the PGF errors. This study investigates the factors inducing the PGF errors of these two methods, through geometric analysis and idealized experiments. The geometric analysis first demonstrates that the terrain slope and the vertical pressure gradient can induce the PGF errors of the classic method, and then generalize the effect of the terrain slope to the effect of the slope of each vertical layer （φ）. More importantly, a new factor, the direction of PGF （a）, is proposed by the geometric analysis, and the effects of φ and a are quantified by tan φ.tan a. When tan φ.tan a is greater than 1/9 or smaller than -10/9, the two terms of PGF of the classic method are of the same order but opposite in sign, and then the PGF errors of the classic method are large. Finally, the effects of three factors on inducing the PGF errors of the classic method are validated by a series of idealized experiments using various terrain types and pressure fields. The experimental results also demonstrate that the PGF errors of the covariant method are affected little by the three factors.
基金supported by the National Basic Research Program of China(973 Program)[grant number 2015CB954102]the National Natural Science Foundation of China[grant number41305095],[grant number 41175064]
文摘A new approach is proposed to use the covariant scalar equations of the a-coordinate （the covariant method）, in which the pressure gradient force （PGF） has only one term in each horizontal momentum equation, and the PGF errors are much reduced in the computational space. In addition, the validity of reducing the PGF errors by this covariant method in the computational and physical space over steep terrain is investigated. First, the authors implement a set of idealized experiments of increasing terrain slope to compare the PGF errors of the covariant method and those of the classic method in the computational space. The results demonstrate that the PGF errors of the covariant method are consistently much-reduced, compared to those of the classic method. More importantly, the steeper the terrain, the greater the reduction in the ratio of the PGF errors via the covariant method. Next, the authors use geometric analysis to further investigate the PGF errors in the physical space, and the results illustrates that the PGF of the covariant method equals that of the classic method in the physical space; namely, the covariant method based on the non-orthogonal a-coordinate cannot reduce the PGF errors in the physical space. However, an orthogonal method can reduce the PGF errors in the physical space. Finally, a set of idealized experiments are carried out to validate the results obtained by the geometric analysis. These results indicate that the covariant method may improve the simulation of variables relevant to pressure, in addition to pressure itself, near steep terrain.
文摘Today, the GNSS （global navigation satellite system） is used for more complicate and accurate applications such as monitoring or stake out works. The truth lies in the fact that in the most of the times not enough attention is paid to the antenna＇s setup. Usually, gross errors are found in the antenna＇s centering, leveling and in the measurement of its height, which are significant. In this paper, a thoroughly analysis of the above mentioned errors is carried out. The influence of these errors in the calculation of the X, Y, Z Cartesian geocentric coordinates and the ~, 2, h ellipsoid geodetic coordinates of a point P on the earth＇s surface, is analyzed and is presented in several diagrams. Also a new convenient method for the accurate measurement of the antenna＇s height is presented and it is strongly proposed. The conclusions outline the magnitude of these errors and prove the significance of the antenna＇s proper setup at the accurate GNSS applications.
