Frequency domain waveform inversion of cross-hole GPR data based on a logarithmic objective function
-
摘要: 波形反演在探地雷达领域的应用已有十余年历史,但绝大部分算例属于时间域波形反演.频率域波形反演由于能够灵活地选择迭代频率并可以使用不同类型的目标函数,因而更加多样化.本文的频率域波形反演基于时间域有限差分(FDTD)法,采用对数目标函数,可在每一次迭代过程中同时或者单独反演介电常数和电导率.文中详细推导了频率域波形反演的理论公式,给出对数目标函数下的梯度表达式,并使用离散傅氏变换(DFT)实现数据的时频变换,能够有效地减少大模型反演的内存需求.在后向残场源的时频域转换过程中,提出仅使用以当前频点为中心的一个窄带数据,可以消除高频无用信号的干扰,获得可靠的反演结果.为加速收敛,采用每迭代十次则反演频率跳跃一定频带宽度的反演策略.实验证明适当的频率跳跃能够在不降低分辨率的基础上有效地提高反演效率.通过两组不同情形下合成数据反演的分析对比,证明基于对数目标函数的波形反演结果准确可靠.最后,将该方法应用到一组实际数据,得到较好的反演结果.Abstract: Ground penetrating radar (GPR) tomography plays an important role in geology, hydrogeology and engineering investigations. Traditional tomography (i.e., the first-arrival times and maximum first-cycle amplitudes) based on ray theory cannot provide high-resolution images because only a fraction of the information contained in the radar data is used in the inversion. In recent years, waveform inversion is one of the biggest concerns because it can provide sub-wavelength resolution images. Waveform inversion has been applied in GPR over ten years, but most of the results are computed in the time domain. In the frequency domain, the choice of inverted frequencies is flexible and different types of objective functions can be used, therefore the results of frequency-domain waveform inversion are more diversified than that of the time domain.#br#In this paper, waveform inversion is implemented by means of a finite-difference time-domain solution of Maxwell's equations and a logarithmic objective function is applied. Permittivity and conductivity can be updated simultaneously or separately at each iterative step. The derivation process of the formulas is described in detail and we show the specific expression of gradient under the logarithmic objective function. It is important to note that the gradient formula of the frequency domain waveform inversion is different from the gradient formula of the time domain waveform inversion. The reason is that the cost function is essentially different. Discrete Fourier transform (DFT) is applied to transform time domain data into the frequency domain, which only increases a few calculations in the inversion. The method can greatly reduce memory requirement when the inverted model is on a large scale. When transforming the back-residual source, we present that only a narrow-band data whose center is the current frequency is used. The method can effectively reduce the interference of high frequency information, so reliable inversion results can be obtained.#br#In order to accelerate the convergence, a special frequency strategy is applied. The inversion frequency skips a number of bandwidth after every ten times iterative step. Results turn out that the strategy can efficiently improve the inversion efficiency and does not influence the resolution. To verify the effect of our algorithm, we test it on two different synthetic data. Then we compare with the results of conventional waveform inversion: (1) Two small cylindrical bodies of permittivity in a layered medium. (2) A layered medium with multiple embedded cylindrical inclusions of permittivity and conductivity. Finally, we apply the logarithmic algorithm to field data.#br#The results on synthetic data show that the logarithmic waveform inversion is better than the conventional waveform inversion. When permittivity and conductivity are inverted at the same time, the results of logarithmic waveform inversion can accurately reconstruct both shape and location of the inclusions. This is because that the amplitude of the residual is normalized by the amplitude of the modeled wavefield. The inversion results of the field data also prove the practical value of logarithmic waveform inversion.
-
Key words:
- Waveform inversion /
- Logarithmic objective function /
- Permittivity /
- Conductivity /
- Cross-hole radar
-
[1] Cai W Y, Qin F H, Schuster G T. 1996. Electromagnetic velocity inversion using 2-D Maxwell's equations. Geophysics, 61(4): 1007-1021.
[2] Chung W, Ha T, Ha W, et al. 2007. Robust seismic waveform inversion using back-propagation algorithm.//77th Annual International Meeting, Expanded Abstracts: 1780-1784.
[3] Ernst J R, Maurer H, Green A G, et al. 2007a. Full-waveform inversion of crosshole radar data based on 2-D finite-difference time-domain solutions of Maxwell's equations. IEEE Transactions on Geoscience and Remote Sensing, 45(9): 2807-2828.
[4] Ernst J R, Green A G, Maurer H, et al. 2007b. Application of a new 2D time-domain full-waveform inversion scheme to crosshole radar data. Geophysics, 72(5): J53-J64.
[5] Fu L, Liu S X, Liu L B, et al. 2014. Airborne ground penetrating radar numerical simulation and reverse time migration. Chinese J. Geophys. (in Chinese), 57(5): 1636-1646, doi: 10.6038/cjg20140526.
[6] Kuroda S, Takeuchi M, Kim H J. 2005. Full waveform inversion algorithm for interpreting cross-borehole GPR data.//75th Annual International Meeting, SEG, Expanded Abstracts: 1176-1179.
[7] Liu L B, Lane J W, Quan Y L. 1998. Radar attenuation tomography using the centroid frequency downshift method. Journal of Applied Geophysics, 40(1-3): 105-116.
[8] Liu S X, Sato M. 2002. Electromagnetic logging technique based on borehole radar. IEEE Transactions on Geoscience and Remote Sensing, 40(9): 2083-2092.
[9] Liu S X, Wu J J, Zhou J F, et al. 2011. Numerical simulations of borehole radar detection for metal ore. IEEE Geoscience and Remote Sensing Letters, 8(2): 308-312.
[10] Liu S X, Lei L L, Fu L, et al. 2014. Application of pre-stack reverse time migration based on FWI velocity estimation to ground penetrating radar data. Journal of Applied Geophysics, 107: 1-7.
[11] Meles G A, van der Kruk J, Greenhalgh S A, et al. 2010. A new vector waveform inversion algorithm for simultaneous updating of conductivity and permittivity parameters from combination crosshole/borehole-to-surface GPR data. IEEE Transactions on Geoscience and Remote Sensing, 48(9): 3391-3407.
[12] Pettinelli E, di Matteo A, Mattei E, et al. 2009. GPR response from buried pipes: Measurement on field site and tomographic reconstructions. IEEE Transactions on Geoscience and Remote Sensing, 47(8): 2639-2645.
[13] Pratt R G. 1990. Inverse theory applied to multi-source cross-hole tomography, Part 2: Elastic wave-equation method. Geophysical Prospecting, 38(3): 311-329.
[14] Pratt R G, Worthington M H. 1990. Inverse theory applied to multi-source cross-hole tomography, Part 1: Acoustic wave-equation method. Geophysical Prospecting, 38(3): 287-310.
[15] Pratt G, Shin C, Hicks G J. 1998. Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion. Geophysical Journal International, 133(2): 341-362.
[16] Pratt R G. 1999. Seismic waveform inversion in the frequency domain, Part 1: Theory and verification in a physical scale model. Geophysics, 64(3): 888-901.
[17] Pratt R G, Shipp R M. 1999. Seismic waveform inversion in the frequency domain, Part 2: Fault delineation in sediments using crosshole data. Geophysics, 64(3): 902-914.
[18] Shin C, Min D J. 2006. Waveform inversion using a logarithmic wavefield. Geophysics, 71(3): R31-R42.
[19] Sirgue L, Pratt R G. 2004. Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies. Geophysics, 69(1): 231-248.
[20] Song Z M, Willianson P R, Pratt R G. 1995. Frequency-domain acoustic-wave modeling and inversion of crosshole data: Part Ⅱ: Inversion method, synthetic experiments and real-data results. Geophysics, 60(3): 796-809.
[21] Tarantola A. 1984. Inversion of seismic reflection data in the acoustic approximation. Geophysics, 49(8): 1259-1266.
[22] Tarantola A. 1986. A strategy for nonlinear elastic inversion of seismic reflection data. Geophysics, 51(10): 1893-1903.
[23] Wang F, Liu S X, Qu X X, et al. 2013. Crosshole radar traveltime tomography without ray tracing using the high accuracy fast marching method. Chinese J. Geophys. (in Chinese), 56(11): 3896-3907, doi: 10.6038/cjg20131131.
[24] Wang Z L, Zhou H, Li G F. 2007. Inversion of ground-penetrating radar data for 2D electric parameters. Chinese J. Geophys. (in Chinese), 50(3): 897-904.
[25] Williamson P R, Worthington M H. 1993. Resolution limits in ray tomography due to wave behavior: Numerical experiments. Geophysics, 58(5): 727-735.
[26] Wu J J, Liu S X, Li Y P, et al. 2014. Study of cross-hole radar tomography using full-waveform inversion. Chinese J. Geophys. (in Chinese), 57(5): 1623-1635, doi: 10.6038/cjg20140525.
[27] Zhou H, Chen H M, Li Q Q, et al. 2014. Waveform inversion of ground-penetrating radar data without needing to extract exciting pulse. Chinese J. Geophys. (in Chinese), 57(6): 1968-1976, doi: 10.6038/cjg20140627.
[28] Zhou H, Sato M, Takenaka T, et al. 2007. Reconstruction from antenna-transformed radar data using a time-domain reconstruction method. IEEE Transactions on Geoscience and Remote Sensing, 45(3): 689-696.
计量
- 文章访问数: 1849
- PDF下载数: 1580
- 施引文献: 0