% Encoding: UTF-8 @Conference{louboutin2018segeow, abstract = {In order to obtain accurate images of the subsurface, anisotropic modeling and imaging is necessary. However, the twenty-one parameter complete wave-equation is too computationally expensive to be of use in this case. The transverse tilted isotropic wave-equation is then the next best feasible representation of the physics to use for imaging. The main complexity arising from transverse tilted isotropic imaging is to model the receiver wavefield (back propagation of the data or data residual) for the imaging condition. Unlike the isotropic or the full physics wave-equations, the transverse tilted isotropic wave-equation is not not self-adjoint. This difference means that time-reversal will not model the correct receiver wavefield and this can lead to incorrect subsurface images. In this work, we derive and implement the adjoint wave-equation to demonstrate the necessity of exact adjoint modeling for anisotropic modeling and compare our result with adjoint-free time-reversed imaging.}, author = {Mathias Louboutin and Philipp A. Witte and Felix J. Herrmann}, booktitle = {SEG Technical Program Expanded Abstracts}, doi = {10.1190/segam2018-2996274.1}, keywords = {Adjoint, anisotropy, finite-differences, Imaging, RTM, SEG, TTI}, month = {10}, note = {(SEG, Anaheim)}, pages = {331-335}, title = {Effects of wrong adjoints for RTM in TTI media}, url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/louboutin2018SEGeow/louboutin2018SEGeow.html}, year = {2018}, bdsk-url-1 = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/louboutin2018SEGeow/louboutin2018SEGeow.html}, bdsk-url-2 = {http://dx.doi.org/10.1190/segam2018-2996274.1} } @Article{virieuxFWI, Title = {An overview of full-waveform inversion in exploration geophysics}, Author = {Virieux,~J. and Operto,~S.}, Journal = {Geophysics}, Year = {2009}, Number = {6}, Pages = {1--26}, Volume = {74}, Doi = {10.1190/1.3238367} } @Article{symesMVAFWI, Title = {Migration velocity analysis and waveform inversion}, Author = {Symes,~W.~W.}, Journal = {Geophysical Prospecting}, Year = {2008}, Pages = {765--790}, Volume = {56}, Doi = {10.1111/j.1365-2478.2008.00698.x} } @Article{vanleeuwenLocMin, Title = {Mitigating local minima in full-waveform inversion by expanding the search space}, Author = {van~Leeuwen,~T. and Herrmann,~F.~J.}, Journal = {Geophys. J. Int.}, Year = {2013}, Pages = {661--667}, Volume = {195}, Doi = {10.1093/gji/ggt258} } @article{vtiwri, author = "Aghamiry, H. and Gholami, A. and Operto, S.", title = "Multi-Parameter ADMM-Based Wavefield Reconstruction Inversion in VTI Acoustic Media", journal= "", year = "2019", volume = "2019", number = "1", pages = "1-5", doi = "https://doi.org/10.3997/2214-4609.201900871", url = "https://www.earthdoc.org/content/papers/10.3997/2214-4609.201900871", publisher = "European Association of Geoscientists & Engineers", issn = "2214-4609", type = "", } @incollection{wang2016full, title={Full-waveform inversion with the reconstructed wavefield method}, author={Wang, Chao and Yingst, David and Farmer, Paul and Leveille, Jacques}, booktitle={SEG Technical Program Expanded Abstracts 2016}, pages={1237--1241}, year={2016}, publisher={Society of Exploration Geophysicists} } @article{huang2018volume, title={Volume source-based extended waveform inversion}, author={Huang, Guanghui and Nammour, Rami and Symes, William W}, journal={Geophysics}, volume={83}, number={5}, pages={R369--R387}, year={2018}, publisher={Society of Exploration Geophysicists} } @article{devito-api, author = {Louboutin, M. and Lange, M. and Luporini, F. and Kukreja, N. and Witte, P. A. and Herrmann, F. J. and Velesko, P. and Gorman, G. J.}, title = {Devito (v3.1.0): an embedded domain-specific language for finite differences and geophysical exploration}, journal = {{Geoscientific Model Development}}, volume = {12}, year = {2019}, number = {3}, pages = {1165--1187}, url = {https://www.geosci-model-dev.net/12/1165/2019/}, doi = {10.5194/gmd-12-1165-2019} } @Article{devito-compiler, author = { {Luporini}, F. and {Lange}, M. and {Louboutin}, M. and {Kukreja}, N. and {H{\"u}ckelheim}, J. and {Yount}, C. and {Witte}, P. and {Kelly}, P.~H.~J. and {Gorman}, G.~J. and {Herrmann}, F.~J. }, title = { Architecture and performance of Devito, a system for automated stencil computation }, journal = { CoRR }, volume = { abs/1807.03032 }, month = { jul }, year = { 2018 }, url = { http://arxiv.org/abs/1807.03032 }, archiveprefix = { arXiv }, eprint = { 1807.03032 } } @article {witte2018alf, title = {A large-scale framework for symbolic implementations of seismic inversion algorithms in Julia}, journal = {Geophysics}, volume = {84}, number = {3}, year = {2019}, note = {(Geophysics)}, month = {03}, pages = {F57-F71}, keywords = {FWI, inversion, LSRTM, Modeling, software}, doi = {10.1190/geo2018-0174.1}, url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2019/witte2018alf/witte2018alf.pdf}, author = {Philipp A. Witte and Mathias Louboutin and Navjot Kukreja and Fabio Luporini and Michael Lange and Gerard J. Gorman and Felix J. Herrmann} } @article{petersWRI, author = "Peters, B. and Herrmann, F.J. and van Leeuwen, T.", title = "Wave-equation Based Inversion with the Penalty Method - Adjoint-state Versus Wavefield-reconstruction Inversion", journal= "", year = "2014", volume = "2014", number = "1", pages = "1-5", doi = "https://doi.org/10.3997/2214-4609.20140704", url = "https://www.earthdoc.org/content/papers/10.3997/2214-4609.20140704", publisher = "European Association of Geoscientists & Engineers", issn = "2214-4609", type = "", } @article{zhang-tti, author = {Yu Zhang and Houzhu Zhang and Guanquan Zhang}, title = {A stable TTI reverse time migration and its implementation}, journal = {GEOPHYSICS}, volume = {76}, number = {3}, pages = {WA3-WA11}, year = {2011}, doi = {10.1190/1.3554411}, url = {https://doi.org/10.1190/1.3554411}, eprint = {https://doi.org/10.1190/1.3554411} } @InBook{bubesatti2016, abstract = { Second-order pseudo-acoustic systems are commonly used in reverse-time migration in tranversely isotropic media to model quasi-P-wave propagation, and have also been used in full-waveform inversion. Systems with the shear wavespeeds set to zero can su er from instabilities. Systems with positive shear wavespeeds behave better, but they can occasionally have di culties when material parameters are strongly variable, particularly when applying the adjoint system in full-waveform inversion. We present here a new family of 2 × 2 second-order pseudo-acoustic systems with positive shear wavespeeds for tranversely isotropic media that are time-reversible, self-adjoint, and conserve a complete energy. They are thus stable for both forward and adjoint computations. In addition, they require only two-thirds of the number of spatial derivatives compared with previous systems. Presentation Date: Tuesday, October 18, 2016 Start Time: 9:15:00 AM Location: 146 Presentation Type: ORAL }, author = {Kenneth Bube and John Washbourne and Raymond Ergas and Tamas Nemeth}, booktitle = {SEG Technical Program Expanded Abstracts 2016}, doi = {10.1190/segam2016-13878451.1}, eprint = {https://library.seg.org/doi/pdf/10.1190/segam2016-13878451.1}, pages = {1110-1114}, title = {Self-adjoint, energy-conserving second-order pseudoacoustic systems for VTI and TTI media for reverse time migration and full-waveform inversion}, url = {e}, year = {2016}, bdsk-url-1 = {https://library.seg.org/doi/abs/10.1190/segam2016-13878451.1}, bdsk-url-2 = {http://dx.doi.org/10.1190/segam2016-13878451.1} } @Article{symes2015iwave, author = {Symes, William W}, journal = {THE RICE INVERSION PROJECT}, pages = {85}, title = {IWAVE structure and basic use cases}, url = {http://www.trip.caam.rice.edu/reports/2014/book.pdf#page=89}, year = {2015}, bdsk-url-1 = {http://www.trip.caam.rice.edu/reports/2014/book.pdf#page=89} } @conference {rizzuti2019SEGadf, title = {A dual formulation for time-domain wavefield reconstruction inversion}, booktitle = {SEG Technical Program Expanded Abstracts}, year = {2019}, note = {(SEG, San Antonio)}, month = {09}, pages = {1480-1485}, abstract = {We illustrate a dual formulation for full-waveform inversion potentially apt to large 3-D problems. It is based on the optimization of the wave equation compliance, under the constraint of data misfit not exceeding a prescribed noise level. In the Lagrangian formulation, model and wavefield state variables are complemented with multipliers having the same dimension of data ("dual data" variables). Analogously to classical wavefield reconstruction inversion, the wavefield unknowns can be projected out in closed form, by solving a version of the augmented wave equation. This leads to a saddle-point problem whose variables are only model and dual data. As such, this formulation represents a model extension, and is potentially robust against local minima. The classical downsides of model extension methods and wavefield reconstruction inversion are here effectively mitigated: storage of the dual variables is affordable, the augmented wave equation is amenable to time-marching finite-difference schemes, and no continuation strategy for penalty parameters is needed, with the prospect of 3-D applications.}, keywords = {3D, Full-waveform inversion, SEG, Time-domain}, doi = {10.1190/segam2019-3216760.1}, url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/rizzuti2019SEGadf/rizzuti2019SEGadf.html}, presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/rizzuti2019SEGadf/rizzuti2019SEGadf_pres.pdf}, software = {https://github.com/slimgroup/Software.rizzuti2019SEGadf}, author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Emmanouil Daskalakis and Felix J. Herrmann} } @conference {peters2019SEGans, title = {A numerical solver for least-squares sub-problems in 3D wavefield reconstruction inversion and related problem formulations}, booktitle = {SEG Technical Program Expanded Abstracts}, year = {2019}, note = {(SEG, San Antonio)}, month = {09}, pages = {1536-1540}, abstract = {Recent years saw a surge of interest in seismic waveform inversion approaches based on quadratic-penalty or augmented-Lagrangian methods, including Wavefield Reconstruction Inversion. These methods typically need to solve a least-squares sub-problem that contains a discretization of the Helmholtz equation. Memory requirements for direct solvers are often prohibitively large in three dimensions, and this limited the examples in the literature to two dimensions. We present an algorithm that uses iterative Helmholtz solvers as a black-box to solve the least-squares problem corresponding to 3D grids. This algorithm enables Wavefield Reconstruction Inversion and related formulations, in three dimensions. Our new algorithm also includes a root-finding method to convert a penalty into a constraint on the data-misfit without additional computational cost, by reusing precomputed quantities. Numerical experiments show that the cost of parallel communication and other computations are small compared to the main cost of solving one Helmholtz problem per source and one per receiver.}, keywords = {Full-waveform inversion, least-squares, numerical linear algebra, penalty method, quadratic constraints, root finding, SEG}, doi = {10.1190/segam2019-3216638.1}, url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/peters2019SEGans/peters2019SEGans.pdf}, presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/peters2019SEGans/peters2019SEGans_pres.pdf}, author = {Bas Peters and Felix J. Herrmann} } @article {vanLeeuwen2013GJImlm, title = {Mitigating local minima in full-waveform inversion by expanding the search space}, journal = {Geophysical Journal International}, volume = {195}, year = {2013}, month = {10}, pages = {661-667}, abstract = {Wave equation based inversions, such as full-waveform inversion and reverse-time migration, are challenging because of their computational costs, memory requirements and reliance on accurate initial models. To confront these issues, we propose a novel formulation of wave equation based inversion based on a penalty method. In this formulation, the objective function consists of a data-misfit term and a penalty term, which measures how accurately the wavefields satisfy the wave equation. This new approach is a major departure from current formulations where forward and adjoint wavefields, which both satisfy the wave equation, are correlated to compute updates for the unknown model parameters. Instead, we carry out the inversions over two alternating steps during which we first estimate the wavefield everywhere, given the current model parameters, source and observed data, followed by a second step during which we update the model parameters, given the estimate for the wavefield everywhere and the source. Because the inversion involves both the synthetic wavefields and the medium parameters, its search space is enlarged so that it suffers less from local minima. Compared to other formulations that extend the search space of wave equation based inversion, our method differs in several aspects, namely (i) it avoids storage and updates of the synthetic wavefields because we calculate these explicitly by finding solutions that obey the wave equation and fit the observed data and (ii) no adjoint wavefields are required to update the model, instead our updates are calculated from these solutions directly, which leads to significant computational savings. We demonstrate the validity of our approach by carefully selected examples and discuss possible extensions and future research.}, doi = {10.1093/gji/ggt258}, url = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalJournalInternational/2013/vanLeeuwen2013GJImlm/vanLeeuwen2013mlm.pdf}, author = {Tristan van Leeuwen and Felix J. Herrmann} } @Article{duveneck, author = {Eric Duveneck and Peter M. Bakker}, title = {Stable P-wave modeling for reverse-time migration in tilted TI media}, journal = {GEOPHYSICS}, volume = {76}, number = {2}, pages = {S65-S75}, year = {2011}, doi = {10.1190/1.3533964}, url = {https://doi.org/10.1190/1.3533964}, eprint = {https://doi.org/10.1190/1.3533964} } @conference {rizzuti2020EAGEtwri, title = {Time-domain wavefield reconstruction inversion for large-scale seismics}, year = {2020}, note = {Accepted in EAGE}, month = {1}, abstract = {Wavefield reconstruction inversion is an imaging technique akin to full-waveform inversion, albeit based on a relaxed version of the wave equation. This relaxation aims to beat the multimodality typical of full-waveform inversion. However it prevents the use of time-marching solvers for the augmented equation and, as a consequence, cannot be straightforwardly employed to large 3D problems. In this work, we formulate a dual version of wavefield reconstruction inversion amenable to explicit time-domain solvers, yielding a robust and scalable inversion technique.}, keywords = {3D, EAGE, Full-waveform inversion, Time-domain}, url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2020/rizzuti2020EAGEtwri/rizzuti2020EAGEtwri.html}, author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Felix J. Herrmann} } @book{rockafellar-1970a, added-at = {2008-03-02T02:12:02.000+0100}, address = {Princeton, N. J.}, author = {Rockafellar, R. Tyrrell}, biburl = {https://www.bibsonomy.org/bibtex/223aa07ea525f6dd11585fc2037a0daf1/dmartins}, callnumber = {UniM Maths 516.08 R59}, description = {robotica-bib}, interhash = {30830becb0a2c5ebca5946b895d9740a}, intrahash = {23aa07ea525f6dd11585fc2037a0daf1}, keywords = {imported}, notes = {A SRL reference.}, publisher = {Princeton University Press}, series = {Princeton Mathematical Series}, timestamp = {2008-03-02T02:14:11.000+0100}, title = {Convex analysis}, year = 1970 } @Article{ thomsen1986, author = {Leon Thomsen}, title = "{Weak elastic anisotropy}", journal = {Geophysics}, year = 1986, month = october, volume = 51, number = 10, pages = {1964-1966} } @article{witteJUDI2019, author = {Philipp A. Witte and Mathias Louboutin and Navjot Kukreja and Fabio Luporini and Michael Lange and Gerard J. Gorman and Felix J. Herrmann}, title = {A large-scale framework for symbolic implementations of seismic inversion algorithms in Julia}, journal = {GEOPHYSICS}, volume = {84}, number = {3}, pages = {F57-F71}, year = {2019}, doi = {10.1190/geo2018-0174.1}, URL = {https://doi.org/10.1190/geo2018-0174.1}, eprint = {https://doi.org/10.1190/geo2018-0174.1} }