Non-Linear Parameter Estimation Method Based on T-Matrix

Wang Ben-Feng,Wu Ru-Shan,Chen Xiao-Hong,Lu Wen-Kai

DOI: https://doi.org/10.6038/cjg20160628

2016-01-01

Abstract:With the development of seismic exploration and exploitation, it is necessary to study the accurate parameter distribution finely describing the subsurface. Full Waveform Inversion (FWI) can help us achieve this goal, however, traditional FWI theory is based on the Born approximation and the performance is dependent on the initial model. This is one of the main factors which prevent FWI being used widely in practice. In order to study the accurate parameter distribution, we studied nonlinear estimation method based on T-matrix which overcomes the initial model dependence.The Inverse Thin -Slab Propagator (ITSP) is studied in detail both theoretically and numerically based on Transmission Matrix (T-matrix) and De Wolf approximation using the 2D acoustic wave equation with constant density, which can overcome the limitations of Born approximation. The ITSP method considers all orders of scattering effects and has no initial model dependence, besides, it is a regularized method in each correction step and has no divergence effects. One correction sweep (half from upper part, half from lower part) is involved and is computationally efficient. Based on the ITSP method, the velocity perturbation can be achieved and the velocity distribution can be obtained finally with the background velocity.We designed three numerical models including ellipse ball with positive velocity perturbations, Gaussian ball with negative velocity perturbations and ellipse ball with both positive and negative velocity perturbations. With the ITSP method, we achieve the velocity distribution based on the known T-matrix and the T-matrix estimation method will be left for future research. The reconstructed results are consistent with the true velocity distributions and the relative errors are minor enough for these three models, which demonstrate the validity of the proposed method. Using the ITSP method, only one correction sweep is involved to achieve the velocity distribution, and it is time efficient. Therefore, this nonlinear estimation method may have wide applications in the future.The accurate velocity distribution is achieved based on the known T-matrix using the ITSP method. 2D numerical examples of the designed models including different perturbation types, demonstrate the validity of the proposed method. The ITSP method considers all scattering effects and overcomes the initial model dependence of Born approximation, besides, it is a regularized method in each correction step and overcomes the divergence effect of Born series for strong velocity perturbation or large perturbation volume. It is also time efficient because only one correction sweep is involved (half from upper part, half from lower part). This method is based on the known T-matrix, suitable for smooth media, then T-matrix estimation method and the improved ITSP method for complex media will be developed in future work.