To reduce the spatial simulation error generated by the finite difference method,previous researchers compute the optimal finite-difference weights always by minimizing the error of spatial dispersion relation.However...To reduce the spatial simulation error generated by the finite difference method,previous researchers compute the optimal finite-difference weights always by minimizing the error of spatial dispersion relation.However,we prove that the spatial simulation error of the finite difference method is associated with the dot product of the spatial dispersion relation of the finite-difference weights and the spectrum of the seismic wavefield.Based on the dot product relation,we construct a L_(2) norm cost function to minimize spatial simulation error.For solving this optimization problem,the seismic wavefield infor-mation in wavenumber region is necessary.Nevertheless,the seismic wavefield is generally obtained by costly forward modeling techniques.To reduce the computational cost,we substitute the spectrum of the seismic wavelet for the spectrum of the seismic wavefield,as the seismic wavelet plays a key role in determining the seismic wavefield.In solving the optimization problem,we design an exhaustive search method to obtain the solution of the L_(2) norm optimization problem.After solving the optimization problem,we are able to achieve the finite-difference weights that minimize spatial simulation error.In theoretical error analyses,the finite-difference weights from the proposed method can output more accurate simulation results compared to those from previous optimization algorithms.Furthermore,we validate our method through numerical tests with synthetic models,which encompass homogenous/inhomogeneous media as well as isotropic and anisotropic media.展开更多
The vertical cable method for acquiring and processing pre-stack 3-D marine seismic data is based on the technology developed by the US Navy for antisubmarine warfare. In order to achieve the maximum utili- zation of ...The vertical cable method for acquiring and processing pre-stack 3-D marine seismic data is based on the technology developed by the US Navy for antisubmarine warfare. In order to achieve the maximum utili- zation of vertical cable field data, a new separation method of the up-going and down-going wave fields of the vertical cable data processing was developed in this paper, which is different from the separation of the down-going and up-going wave fields of normal VSP data processing. In tests with synthetic modeling data and actual field data, this newly developed method performs well and is also computationally simpler without pre-assumption conditions.展开更多
We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this te...We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete （NAD） operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction （LWC） method and the eighth-order staggered-grid （SG）method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.展开更多
Stability is the key to inverse Q-filtering. In this paper we present a stable approach to inverse Q-filtering, based on the theory of wavefield downward continuation. It is implemented in a layered manner, assuming a...Stability is the key to inverse Q-filtering. In this paper we present a stable approach to inverse Q-filtering, based on the theory of wavefield downward continuation. It is implemented in a layered manner, assuming a layered-earth Q model. For each individual constant Q layer, the seismic wavefield recorded at the surface is first extrapolated down to the top of the current layer and a constant Q inverse filter is then applied to the current layer. When extrapolating within the overburden, a stable wavefield continuation algorithm in combination with a stabilization factor is applied. This avoids accumulating inverse Q-filter errors within the overburden. Within the current constant Q layer, we use Gabor spectral analysis on the signals to pick time-variant gain-constrained frequencies and then deduce the corresponding gain-constrained amplitudes to stabilize the inverse Q-filtering algorithm. The algorithm is tested and verified application to field data.展开更多
Prestack reverse time migration（PSTM） is a common imaging method; however low-frequency noises reduce the structural imaging precision. Thus, the suppression of migration noises must be considered. The generation me...Prestack reverse time migration（PSTM） is a common imaging method; however low-frequency noises reduce the structural imaging precision. Thus, the suppression of migration noises must be considered. The generation mechanism of low-frequency noises is analyzed and the up-, down-, left-, and right-going waves are separated using the Poynting vector of the acoustic wave equation. The computational complexity and memory capacitance of the proposed method are far smaller than that required when using the conventional separation algorithm of 2D Fourier transform. The normalized wavefield separation crosscorrelation imaging condition is used to suppress low-frequency noises in reverse time migration and improve the imaging precision. Numerical experiments using the Marmousi model are performed and the results show that the up-, down-, left-, and right-going waves are well separated in the continuation of the wavefield using the Poynting vector. We compared the imaging results with the conventional method, Laplacian filtering, and wavefield separation with the 2D Fourier transform. The comparison shows that the migration noises are well suppressed using the normalized wavefield separation cross-correlation imaging condition and higher precision imaging results are obtained.展开更多
Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton...Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation（capillary pressure）.So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.展开更多
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group(No.SNKJ2022A06-R23)+1 种基金the Funds of Creative Research Groups of China(No.41821002)the Major Scientific and Technological Projects of CNPC(No.ZD2019-183-003).
文摘To reduce the spatial simulation error generated by the finite difference method,previous researchers compute the optimal finite-difference weights always by minimizing the error of spatial dispersion relation.However,we prove that the spatial simulation error of the finite difference method is associated with the dot product of the spatial dispersion relation of the finite-difference weights and the spectrum of the seismic wavefield.Based on the dot product relation,we construct a L_(2) norm cost function to minimize spatial simulation error.For solving this optimization problem,the seismic wavefield infor-mation in wavenumber region is necessary.Nevertheless,the seismic wavefield is generally obtained by costly forward modeling techniques.To reduce the computational cost,we substitute the spectrum of the seismic wavelet for the spectrum of the seismic wavefield,as the seismic wavelet plays a key role in determining the seismic wavefield.In solving the optimization problem,we design an exhaustive search method to obtain the solution of the L_(2) norm optimization problem.After solving the optimization problem,we are able to achieve the finite-difference weights that minimize spatial simulation error.In theoretical error analyses,the finite-difference weights from the proposed method can output more accurate simulation results compared to those from previous optimization algorithms.Furthermore,we validate our method through numerical tests with synthetic models,which encompass homogenous/inhomogeneous media as well as isotropic and anisotropic media.
文摘The vertical cable method for acquiring and processing pre-stack 3-D marine seismic data is based on the technology developed by the US Navy for antisubmarine warfare. In order to achieve the maximum utili- zation of vertical cable field data, a new separation method of the up-going and down-going wave fields of the vertical cable data processing was developed in this paper, which is different from the separation of the down-going and up-going wave fields of normal VSP data processing. In tests with synthetic modeling data and actual field data, this newly developed method performs well and is also computationally simpler without pre-assumption conditions.
基金This research was supported by the National Natural Science Foundation of China （Nos. 41230210 and 41204074）, the Science Foundation of the Education Department of Yunnan Province （No. 2013Z152）, and Statoil Company （Contract No. 4502502663）.
文摘We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete （NAD） operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction （LWC） method and the eighth-order staggered-grid （SG）method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.
基金This research is sponsored by the National"973"Project(No.2007CB209603)the"863"Project(No.2006AA06Z108)
文摘Stability is the key to inverse Q-filtering. In this paper we present a stable approach to inverse Q-filtering, based on the theory of wavefield downward continuation. It is implemented in a layered manner, assuming a layered-earth Q model. For each individual constant Q layer, the seismic wavefield recorded at the surface is first extrapolated down to the top of the current layer and a constant Q inverse filter is then applied to the current layer. When extrapolating within the overburden, a stable wavefield continuation algorithm in combination with a stabilization factor is applied. This avoids accumulating inverse Q-filter errors within the overburden. Within the current constant Q layer, we use Gabor spectral analysis on the signals to pick time-variant gain-constrained frequencies and then deduce the corresponding gain-constrained amplitudes to stabilize the inverse Q-filtering algorithm. The algorithm is tested and verified application to field data.
基金supported by the National Natural Science Foundation of China(No.41174087,41204089)the National Oil and Gas Major Project(No.2011ZX05005-005)
文摘Prestack reverse time migration（PSTM） is a common imaging method; however low-frequency noises reduce the structural imaging precision. Thus, the suppression of migration noises must be considered. The generation mechanism of low-frequency noises is analyzed and the up-, down-, left-, and right-going waves are separated using the Poynting vector of the acoustic wave equation. The computational complexity and memory capacitance of the proposed method are far smaller than that required when using the conventional separation algorithm of 2D Fourier transform. The normalized wavefield separation crosscorrelation imaging condition is used to suppress low-frequency noises in reverse time migration and improve the imaging precision. Numerical experiments using the Marmousi model are performed and the results show that the up-, down-, left-, and right-going waves are well separated in the continuation of the wavefield using the Poynting vector. We compared the imaging results with the conventional method, Laplacian filtering, and wavefield separation with the 2D Fourier transform. The comparison shows that the migration noises are well suppressed using the normalized wavefield separation cross-correlation imaging condition and higher precision imaging results are obtained.
基金supported by the 973 Program (Grant No.2007CB209505)the National Natural Science Foundation of China (Grant No.40674061,40704019)
文摘Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation（capillary pressure）.So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.
