科研之友
微信
新浪微博
Facebook
分享链接
摘要
As a highly efficient near-surface geophysical technology, cross-hole GPR can provide high-resolution images between two boreholes. Waveform inversion, of which the resolution reaches sub-wavelength, can reconstruct accurate subsurface dielectric parameters such as permittivity and conductivity. We can get accurate and high-resolution subsurface dielectric information by combining waveform inversion and cross-hole GPR data. In the waveform inversion of field data, the source wavelet is unknown. Normally, the source wavelet can be estimated by using a deconvolution method. In this method, the source wavelet is a new unknown parameter added in the inversion and updated with iteration. When the results of inversion are same as the true models, the estimated source wavelet is same as the true source wavelet. This method is useful in the inversion of synthetic data, but it does not perform well in the field data. Lots of intervention are needed to choose the best source wavelet after the estimation. The reason may be that the subsurface medium of field data is highly complex and the signal-to-noise ratio of the field data is low. @@@ In this paper, we realize a source-independent time-domain waveform inversion. Firstly, the observed wavefields are convolved with a reference trace of the modeled wavefield, then the modeled wavefields are convolved with a reference trace of observed wavefield. A new object function is based on the convolved wavefields. In theory, no matter which source wavelet is used in the forward modeling. The source wavelet of the observed and the modeled wavefields are equally convolved with both terms in the object function, so that the effect of the source wavelet is removed. Another important feature of this object function is that the modeled wavefields act as a low-pass filter of the wavefield based on the frequency range of the source used for modeling. Though it is very easy to employ a frequency-selection strategy, we don't discuss the strategy here. For each shot gather, we always choose the trace which is nearest to the antenna as the reference trace in this study. @@@ To check the effect of our algorithm, we simulate three different synthetic data: 1) Single cylindrical body of permittivity. 2) Single cylindrical body of conductivity. And 3) layered media with multiple embedded cylindrical inclusions of permittivity and conductivity. Finally, we apply the source-independent waveform inversion to the field data of Guizhou and Xiuyan. @@@ The first two results of synthetic data show that our algorithm can provide high-resolution images. When permittivity and conductivity are simultaneously updated, the results of synthetic data are good. The layers are well reconstructed and the two pipes can be clearly distinguished. But there are many false features which are caused by the artifacts of sources. We have to do convolution and crosscorrelation to compute the gradients. These convolution and crosscorrelation operations increase the nonlinearity of the inversion. The artifacts of the sources become stronger after the convolution and crosscorrelation is hard to remove when the complex models are of a small scale. Then we apply the source-independent waveform inversion to the field data of Guizhou and Xiuyan. The signal-to-noise ratio of Guizhou field data is high, so both the source independent and source-estimated waveform inversion can obtain accurate results. We think the results of source-independent waveform inversion are better. The reason may be the difference between the estimated source wavelet and the true source wavelet, though they are very similar. However, it is difficult to get the ideal estimated source wavelet of Xiuyan field data because its signal-to-noise ratio is low. So the results of the Xiuyan field data can better reflect the value of our source-independent waveform inversion.
- 出版日期2016-12
- 单位吉林大学
