We propose a novel method for seismic noise attenuation by applying nonstationary polynomial fitting （NPF）, which can estimate coherent components with amplitude variation along the event. The NPF with time-varying ...We propose a novel method for seismic noise attenuation by applying nonstationary polynomial fitting （NPF）, which can estimate coherent components with amplitude variation along the event. The NPF with time-varying coefficients can adaptively estimate the coherent components. The smoothness of the polynomial coefficients is controlled by shaping regularization. The signal is coherent along the offset axis in a common midpoint （CMP） gather after normal moveout （NMO）. We use NPF to estimate the effective signal and thereby to attenuate the random noise. For radial events-like noise such as ground roll, we first employ a radial trace （RT） transform to transform the data to the time-velocity domain. Then the NPF is used to estimate coherent noise in the RT domain. Finally, the coherent noise is adaptively subtracted from the noisy dataset. The proposed method can effectively estimate coherent noise with amplitude variations along the event and there is no need to propose that noise amplitude is constant. Results of synthetic and field data examples show that, compared with conventional methods such as stationary polynomial fitting and low cut filters, the proposed method can effectively suppress seismic noise and preserve the signals.展开更多
Through solving the Zoeppritz's partial derivative equations, we have obtained accurate partial derivatives of reflected coefficients of seismic wave with respect to Pand S-wave velocities.With those partial deriv...Through solving the Zoeppritz's partial derivative equations, we have obtained accurate partial derivatives of reflected coefficients of seismic wave with respect to Pand S-wave velocities.With those partial derivatives, a multi-angle inversion is developed for seismic wave velocities.Numerical examples of different formation models show that if the number of iterations goes over 10, the relative error of inversion results is less than 1%, whether or not there is interference among the reflection waves.When we only have the reflected seismograms of P-wave, and only invert for velocities of P-wave, the multi-angle inversion is able to obtain a high computation precision.When we have the reflected seismograms of both P-wave and VS-wave, and simultaneously invert for the velocities of P-wave and VS-wave, the computation precisions of VS-wave velocities improves gradually with the increase of the number of angles, but the computation precision of P-wave velocities becomes worse.No matter whether the reflected seismic waves from the different reflection interface are coherent or non-coherent, this method is able to achieve a higher computation precision.Because it is based on the accurate solution of the gradient of SWRCs without any additional restriction, the multi-angle inversion method can be applied to seismic inversion of total angles.By removing the difficulties caused by simplified Zoeppritz formulas that the conventional AVO technology struggles with, the multiangle inversion method extended the application range of AVO technology and improved the computation precision and speed of inversion of seismic wave velocities.展开更多
Elastic impedance （EI） inversion has been widely used in industry to estimate kinds of elastic parameters to distinguish lithology or even fluid. However, it is found that conventional three-term elastic impedance f...Elastic impedance （EI） inversion has been widely used in industry to estimate kinds of elastic parameters to distinguish lithology or even fluid. However, it is found that conventional three-term elastic impedance formula is unstable even with slight random noise in seismic data, due to the m-conditioned co- efficient matrix of elastic parameters. We presented two-term Fatti elastic impedance inversion method, which is more robust and accurate than conventional three-term elastic impedance inversion. In our method, density is ignored to increase the robustness of inversion matrix. Besides, P-impedance and S-impedance, which are less sensitive to random noise, are inverted instead of Vp and Vs in conventional three-term elastic impedance. To make the inversion more stable, we defined the range of K value as a con- straint. Synthetic tests claim that this method can obtain promising results with low SNR （signal noise ratio） seismic data. With the application of the method in a 2D line data, we achieved λp, μp and Vp/Vs sections, which matched the drilled well perfectly, indicating the potential of the method in reservoir prediction.展开更多
Energy loss in porous media containing fluids is typically caused by a variety of dynamic mechanisms.In the Biot theory,energy loss only includes the frictional dissipation between the solid phase and the fluid phase,...Energy loss in porous media containing fluids is typically caused by a variety of dynamic mechanisms.In the Biot theory,energy loss only includes the frictional dissipation between the solid phase and the fluid phase,resulting in underestimation of the dispersion and attenuation of the waves in the low frequency range.To develop a dynamic model that can predict the high dispersion and strong attenuation of waves at the seismic band,we introduce viscoelasticity into the Biot model and use fractional derivatives to describe the viscoelastic mechanism,and finally propose a new wave propagation model.Unlike the Biot model,the proposed model includes the intrinsic dissipation of the solid frame.We investigate the effects of the fractional order parameters on the dispersion and attenuation of the P-and S-waves using several numerical experiments.Furthermore,we use several groups of experimental data from different fluid-saturated rocks to testify the validity of the new model.The results demonstrate that the new model provides more accurate predictions of high dispersion and strong attenuation of different waves in the low frequency range.展开更多
基金supported by the National Basic Research Program of China (973 program, grant 2007CB209606) the National High Technology Research and Development Program of China (863 program, grant 2006AA09A102-09)
文摘We propose a novel method for seismic noise attenuation by applying nonstationary polynomial fitting （NPF）, which can estimate coherent components with amplitude variation along the event. The NPF with time-varying coefficients can adaptively estimate the coherent components. The smoothness of the polynomial coefficients is controlled by shaping regularization. The signal is coherent along the offset axis in a common midpoint （CMP） gather after normal moveout （NMO）. We use NPF to estimate the effective signal and thereby to attenuate the random noise. For radial events-like noise such as ground roll, we first employ a radial trace （RT） transform to transform the data to the time-velocity domain. Then the NPF is used to estimate coherent noise in the RT domain. Finally, the coherent noise is adaptively subtracted from the noisy dataset. The proposed method can effectively estimate coherent noise with amplitude variations along the event and there is no need to propose that noise amplitude is constant. Results of synthetic and field data examples show that, compared with conventional methods such as stationary polynomial fitting and low cut filters, the proposed method can effectively suppress seismic noise and preserve the signals.
基金supported by Funding Project for Academic Human Resources Development in Institutions of Higher Learning Under the Jurisdiction of Beijing Municipality(PHR(IHLB))(Grant No.PHR201107145)
文摘Through solving the Zoeppritz's partial derivative equations, we have obtained accurate partial derivatives of reflected coefficients of seismic wave with respect to Pand S-wave velocities.With those partial derivatives, a multi-angle inversion is developed for seismic wave velocities.Numerical examples of different formation models show that if the number of iterations goes over 10, the relative error of inversion results is less than 1%, whether or not there is interference among the reflection waves.When we only have the reflected seismograms of P-wave, and only invert for velocities of P-wave, the multi-angle inversion is able to obtain a high computation precision.When we have the reflected seismograms of both P-wave and VS-wave, and simultaneously invert for the velocities of P-wave and VS-wave, the computation precisions of VS-wave velocities improves gradually with the increase of the number of angles, but the computation precision of P-wave velocities becomes worse.No matter whether the reflected seismic waves from the different reflection interface are coherent or non-coherent, this method is able to achieve a higher computation precision.Because it is based on the accurate solution of the gradient of SWRCs without any additional restriction, the multi-angle inversion method can be applied to seismic inversion of total angles.By removing the difficulties caused by simplified Zoeppritz formulas that the conventional AVO technology struggles with, the multiangle inversion method extended the application range of AVO technology and improved the computation precision and speed of inversion of seismic wave velocities.
基金the sponsorship of the National Natural Science Foundation of China (Nos.41004096 and 41230318) for funding this research
文摘Elastic impedance （EI） inversion has been widely used in industry to estimate kinds of elastic parameters to distinguish lithology or even fluid. However, it is found that conventional three-term elastic impedance formula is unstable even with slight random noise in seismic data, due to the m-conditioned co- efficient matrix of elastic parameters. We presented two-term Fatti elastic impedance inversion method, which is more robust and accurate than conventional three-term elastic impedance inversion. In our method, density is ignored to increase the robustness of inversion matrix. Besides, P-impedance and S-impedance, which are less sensitive to random noise, are inverted instead of Vp and Vs in conventional three-term elastic impedance. To make the inversion more stable, we defined the range of K value as a con- straint. Synthetic tests claim that this method can obtain promising results with low SNR （signal noise ratio） seismic data. With the application of the method in a 2D line data, we achieved λp, μp and Vp/Vs sections, which matched the drilled well perfectly, indicating the potential of the method in reservoir prediction.
基金the National Natural Science Foundation of China(Grant Nos.91730306 and 41390452)the Shengli Oilfield Company(Grant No.30200020-18ZC0613-0030)。
文摘Energy loss in porous media containing fluids is typically caused by a variety of dynamic mechanisms.In the Biot theory,energy loss only includes the frictional dissipation between the solid phase and the fluid phase,resulting in underestimation of the dispersion and attenuation of the waves in the low frequency range.To develop a dynamic model that can predict the high dispersion and strong attenuation of waves at the seismic band,we introduce viscoelasticity into the Biot model and use fractional derivatives to describe the viscoelastic mechanism,and finally propose a new wave propagation model.Unlike the Biot model,the proposed model includes the intrinsic dissipation of the solid frame.We investigate the effects of the fractional order parameters on the dispersion and attenuation of the P-and S-waves using several numerical experiments.Furthermore,we use several groups of experimental data from different fluid-saturated rocks to testify the validity of the new model.The results demonstrate that the new model provides more accurate predictions of high dispersion and strong attenuation of different waves in the low frequency range.
