Crosshole ground-penetrating radar full-waveform inversion by combining optimal-transport and least-squares distances

被引:0
作者
Hunziker, Jurg [1 ]
Meles, Giovanni [1 ]
Linde, Niklas [1 ]
机构
[1] Univ Lausanne, Inst Earth Sci, Lausanne, Switzerland
关键词
Ground penetrating radar; Full-waveform inversion; Optimal-transport distance; Open-source software; 2-D; GRADIENT;
D O I
10.1016/j.jappgeo.2025.105655
中图分类号
P [天文学、地球科学];
学科分类号
07 ;
摘要
Gradient-based full-waveform inversion using the least-squares distance converges to a local minimum if the starting model is not close enough to the global minimum. In this article, we propose a full-waveform inversion algorithm for crosshole ground penetrating radar tomography that is independent of the starting model. This is achieved by using an optimal-transport distance, which features one broad global minimum and almost no local minima for common ground penetrating radar full-waveform inversion problems. As calculating an optimaltransport distance is computationally expensive and as the broad global minimum of the optimal-transport distance is impacted by numerical approximations in the vicinity of the optimal solution, we propose an algorithm that uses the optimal-transport distance only in the first few iterations of the inversion and switches to the least-squares distance once a solution close enough to the true solution has been found. An additional feature of the proposed algorithm is the sparse and explicit calculation of the gradient using random sets of master points and subsequent interpolation. This avoids the high sensitivity close to the antenna locations, which prohibits meaningful model updates in these regions. Furthermore, this approach smooths the gradient and, thus, the model without an explicit model-regularization term in the objective function. Finally, calculating the gradient explicitly allows to easily implement other distance measures as no update of the gradient calculation becomes necessary. As the explicit calculation of the gradient can be done in parallel and only a limited number of master points is necessary, the additional computational cost is limited. We demonstrate the capabilities of the algorithm on a simple synthetic dataset showing that the proposed algorithm works much better on our example dataset than if only least-squares or optimal-transport distances are used.
引用
收藏
页数:14
相关论文
共 35 条
[1]   Pilot points method incorporating prior information for solving the groundwater flow inverse problem [J].
Alcolea, Andres ;
Carrera, Jesus ;
Medina, Agustin .
ADVANCES IN WATER RESOURCES, 2006, 29 (11) :1678-1689
[2]   Image appraisal for 2-D and 3-D electromagnetic inversion [J].
Alumbaugh, DL ;
Newman, GA .
GEOPHYSICS, 2000, 65 (05) :1455-1467
[3]   Laplace-domain waveform inversion versus refraction-traveltime tomography [J].
Bae, Ho Seuk ;
Pyun, Sukjoon ;
Shin, Changsoo ;
Marfurt, Kurt J. ;
Chung, Wookeen .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2012, 190 (01) :595-606
[4]   Electromagnetic Full Waveform Inversion Based on Quadratic Wasserstein Metric [J].
Deng, Jian ;
Zhu, Peimin ;
Kofman, Wlodek ;
Jiang, Jinpeng ;
Yuan, Yuefeng ;
Herique, Alain .
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, 2022, 70 (12) :11934-11945
[5]   APPLICATION OF THE WASSERSTEIN METRIC TO SEISMIC SIGNALS [J].
Engquist, Bjoern ;
Froese, Brittany D. .
COMMUNICATIONS IN MATHEMATICAL SCIENCES, 2014, 12 (05) :979-988
[6]   Application of a new 2D time-domain full-waveform inversion scheme to crosshole radar data [J].
Ernst, Jacques R. ;
Green, Alan G. ;
Maurer, Hansruedi ;
Holliger, Klaus .
GEOPHYSICS, 2007, 72 (05) :J53-J64
[7]  
Ha W, 2012, GEOPHYSICS, V77, pR199, DOI [10.1190/geo2011-0411.1, 10.1190/GEO2011-0411.1]
[8]   Bayesian full-waveform tomography with application to crosshole ground penetrating radar data [J].
Hunziker, Jurg ;
Laloy, Eric ;
Linde, Niklas .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2019, 218 (02) :913-931
[9]   Numerical modeling of ground-penetrating radar in 2-D using MATLAB [J].
Irving, James ;
Knight, Rosemary .
COMPUTERS & GEOSCIENCES, 2006, 32 (09) :1247-1258
[10]  
Kingsbury D, 2015, P1, DOI [DOI 10.1021/bk-2015-1214.ch001, 10.48550/arXiv.1312.6114]