TY - JOUR
T1 - Parsimonious seismic tomography with Poisson voronoi projections
T2 - Methodology and validation
AU - Fang, Hongjian
AU - Van Der Hilst, Robert D.
AU - De Hoop, Maarten V.
AU - Kothari, Konik
AU - Gupta, Sidharth
AU - Dokmanić, Ivan
N1 - Publisher Copyright:
© Seismological Society of America
PY - 2019
Y1 - 2019
N2 - Ill-posed seismic inverse problems are often solved using Tikhonov-type regularization, that is, incorporation of damping and smoothing to obtain stable results. This typically results in overly smooth models, poor amplitude resolution, and a difficult choice between plausible models. Recognizing that the average of parameters can be better constrained than individual parameters, we propose a seismic tomography method that stabilizes the inverse problem by projecting the original high-dimension model space onto random low-dimension subspaces and then infers the high-dimensional solution from combinations of such subspaces. The subspaces are formed by functions constant in Poisson Voronoi cells, which can be viewed as the mean of parameters near a certain location. The low-dimensional problems are better constrained, and image reconstruction of the subspaces does not require explicit regularization. Moreover, the low-dimension subspaces can be recovered by subsets of the whole dataset, which increases efficiency and offers opportunities to mitigate uneven sampling of the model space. The final (high-dimension) model is then obtained from the low-dimension images in different subspaces either by solving another normal equation or simply by averaging the low-dimension images. Importantly, model uncertainty can be obtained directly from images in different subspaces. Synthetic tests show that our method outperforms conventional methods both in terms of geometry and amplitude recovery. The application to southern California plate boundary region also validates the robustness of our method by imaging geologically consistent features as well as strong along-strike variations of San Jacinto fault that are not clearly seen using conventional methods.
AB - Ill-posed seismic inverse problems are often solved using Tikhonov-type regularization, that is, incorporation of damping and smoothing to obtain stable results. This typically results in overly smooth models, poor amplitude resolution, and a difficult choice between plausible models. Recognizing that the average of parameters can be better constrained than individual parameters, we propose a seismic tomography method that stabilizes the inverse problem by projecting the original high-dimension model space onto random low-dimension subspaces and then infers the high-dimensional solution from combinations of such subspaces. The subspaces are formed by functions constant in Poisson Voronoi cells, which can be viewed as the mean of parameters near a certain location. The low-dimensional problems are better constrained, and image reconstruction of the subspaces does not require explicit regularization. Moreover, the low-dimension subspaces can be recovered by subsets of the whole dataset, which increases efficiency and offers opportunities to mitigate uneven sampling of the model space. The final (high-dimension) model is then obtained from the low-dimension images in different subspaces either by solving another normal equation or simply by averaging the low-dimension images. Importantly, model uncertainty can be obtained directly from images in different subspaces. Synthetic tests show that our method outperforms conventional methods both in terms of geometry and amplitude recovery. The application to southern California plate boundary region also validates the robustness of our method by imaging geologically consistent features as well as strong along-strike variations of San Jacinto fault that are not clearly seen using conventional methods.
UR - http://www.scopus.com/inward/record.url?scp=85077110744&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85077110744&partnerID=8YFLogxK
U2 - 10.1785/0220190141
DO - 10.1785/0220190141
M3 - Article
AN - SCOPUS:85077110744
SN - 0895-0695
VL - 91
SP - 343
EP - 355
JO - Seismological Research Letters
JF - Seismological Research Letters
IS - 1
ER -