Abstract:Full Waveform Inversion (FWI) is an advanced geophysical inversion technique. In fields such as oil exploration and geology, FWI is used for providing images of subsurface structures with higher resolution. The conventional algorithm minimises the misfit error by calculating the least squares of the wavefield solutions between observed data and simulated data, followed by gradient direction and model update increment. Since the gradient is calculated by forward and backward wavefields, the high-accuracy model update relies on accurate forward and backward wavefield modelling. However, the quality of wavefield solutions obtained in practical situations could be poor and does not meet the requirements of high-resolution FWI. Specifically, the low-frequency wavefield is easily affected by noise and downsampling, which influences data quality, while the high-frequency wavefield is susceptible to spatial aliasing effects that produce imaging artefacts. Therefore, we propose using an algorithm called Sparse Relaxation Regularised Regression to optimise the wavefield solution in frequency domain FWI, which is the forward and backward wavefield obtained from the Helmholtz equation, and thus improve the accuracy of the FWI. The sparse relaxed regularised regression algorithm combines sparsity and regularisation, allowing the broadband FWI to reduce the effects of noise and outliers, which can provide data supplementation in the low-frequency band and anti-aliasing in the high-frequency band. Our numerical examples demonstrate the wavefield optimization effect of the sparse relaxed regularised regression-based algorithm in various cases. The accuracy and stability of the improved algorithm are verified in comparison to the Tikhonov regularisation algorithm.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to improve the quality of low - frequency data and the spatial aliasing effect of high - frequency data in full waveform inversion (FWI) in order to achieve higher - resolution imaging of underground structures. Specifically, the paper focuses on the following aspects:
1. **Quality problems of low - frequency data**:
- In practical applications, low - frequency data is easily affected by noise and undersampling, resulting in a decline in data quality. This will affect the accuracy of FWI because low - frequency information is crucial for overcoming nonlinear problems and promoting the convergence of the inversion process to the global optimal solution.
2. **Spatial aliasing effect of high - frequency data**:
- High - frequency data is easily affected by the spatial aliasing effect, producing imaging artifacts. This will lead to the loss of high - frequency information, thereby affecting the accurate recovery of the underground velocity structure.
3. **Limitations of traditional methods**:
- Traditional FWI algorithms rely on the least - squares method to reduce the error between observed data and simulated data. However, this method is not effective in dealing with low - density sampling data and is prone to getting trapped in local minima.
To solve the above problems, the paper proposes an algorithm based on sparse relaxation regularized regression (SR3) to optimize the forward and backward wavefields in frequency - domain FWI, thereby improving the accuracy and resolution of FWI. By introducing sparsity and regularization constraints, the SR3 algorithm can reduce the influence of noise and outliers, provide data supplementation in the low - frequency band, and achieve anti - aliasing effects in the high - frequency band.
### Specific solutions
- **Sparse relaxation regularized regression (SR3)**:
- It combines sparsity and regularization, optimizes the wavefield solution through an auxiliary matrix, compresses the nonlinear wavefield loss space, and narrows the singular value range.
- By selecting appropriate regularization parameters, a balance can be found between sparse constraints and the smoothness of the solution, thereby providing more stable denoising and data completion effects.
- **Numerical experiments**:
- The superiority of the SR3 algorithm was verified using two homogeneous medium models and two synthetic models (the 2004 BP model and the SEG/EAGE Overthrust model).
- The experimental results show that the SR3 algorithm performs excellently in dealing with low - density sampling data, resisting high - frequency spatial aliasing effects, and improving the resolution of high - background - noise data sets.
Through these methods, the paper aims to improve the accuracy and stability of FWI in practical applications, especially in terms of low - frequency data quality and high - frequency data spatial aliasing effects.