Quantifying uncertainty in inverse scattering problems set in layered environments
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2025
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866909571412393984 |
|---|---|
| author | Abugattas, Carolina Carpio, Ana Cebrián, Elena Oleaga, Gerardo |
| author_facet | Abugattas, Carolina Carpio, Ana Cebrián, Elena Oleaga, Gerardo |
| contents | The attempt to solve inverse scattering problems often leads to optimization and sampling problems that require handling moderate to large amounts of partial differential equations acting as constraints. We focus here on determining inclusions in a layered medium from the measurement of wave fields on the surface, while quantifying uncertainty and addressing the effect of wave solver quality. Inclusions are characterized by a few parameters describing their material properties and shapes. We devise algorithms to estimate the most likely configurations by optimizing cost functionals with Bayesian regularizations and wave constraints. In particular, we design an automatic Levenberg-Marquardt-Fletcher type scheme based on the use of algorithmic differentiation and adaptive finite element meshes for time dependent wave equation constraints with changing inclusions. In synthetic tests with a single frequency, this scheme converges in few iterations for increasing noise levels. To attain a global view of other possible high probability configurations and asymmetry effects we resort to parallelizable affine invariant Markov Chain Monte Carlo methods, at the cost of solving a few million wave problems. This forces the use of prefixed meshes. While the optimal configurations remain similar, we encounter additional high probability inclusions influenced by the prior information, the noise level and the layered structure, effect that can be reduced by considering more frequencies. We analyze the effect on the calculations of working with adaptive and fixed meshes, under a simple choice of non-reflecting boundary conditions in truncated layered domains for which we establish wellposedness and convergence results. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2504_05776 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Quantifying uncertainty in inverse scattering problems set in layered environments Abugattas, Carolina Carpio, Ana Cebrián, Elena Oleaga, Gerardo Numerical Analysis Optimization and Control Computational Physics Data Analysis, Statistics and Probability Geophysics The attempt to solve inverse scattering problems often leads to optimization and sampling problems that require handling moderate to large amounts of partial differential equations acting as constraints. We focus here on determining inclusions in a layered medium from the measurement of wave fields on the surface, while quantifying uncertainty and addressing the effect of wave solver quality. Inclusions are characterized by a few parameters describing their material properties and shapes. We devise algorithms to estimate the most likely configurations by optimizing cost functionals with Bayesian regularizations and wave constraints. In particular, we design an automatic Levenberg-Marquardt-Fletcher type scheme based on the use of algorithmic differentiation and adaptive finite element meshes for time dependent wave equation constraints with changing inclusions. In synthetic tests with a single frequency, this scheme converges in few iterations for increasing noise levels. To attain a global view of other possible high probability configurations and asymmetry effects we resort to parallelizable affine invariant Markov Chain Monte Carlo methods, at the cost of solving a few million wave problems. This forces the use of prefixed meshes. While the optimal configurations remain similar, we encounter additional high probability inclusions influenced by the prior information, the noise level and the layered structure, effect that can be reduced by considering more frequencies. We analyze the effect on the calculations of working with adaptive and fixed meshes, under a simple choice of non-reflecting boundary conditions in truncated layered domains for which we establish wellposedness and convergence results. |
| title | Quantifying uncertainty in inverse scattering problems set in layered environments |
| topic | Numerical Analysis Optimization and Control Computational Physics Data Analysis, Statistics and Probability Geophysics |
| url | https://arxiv.org/abs/2504.05776 |