@conference {gahlot2023IMAGEtsm,
	title = {Time-lapse seismic monitoring of geological carbon storage with the nonlinear joint recovery model},
	booktitle = {International Meeting for Applied Geoscience and Energy},
	year = {2023},
	note = {(IMAGE, Houston)},
	month = {08},
	abstract = {During time-lapse seismic monitoring, weak 4D signal below the level of inversion or migration artifacts poses challenges. To address these, low-cost randomized non-replicated acquisitions and linear joint recovery model (JRM) have been introduced to take advantage of the shared information between different vintages in the time-lapse seismic data and subsurface structure undergoing localized changes. Since the relationship between seismic data and subsurface properties is seldomly linear, we propose a more versatile nonlinear JRM, which extends linear JRM to nonlinear forward modeling.},
	keywords = {CCS, JRM, monitoring, SEG},
	url = {https://slimgroup.github.io/IMAGE2023/NonLinear-JRM/abstract.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2023/gahlot2023IMAGEtsm},
	author = {Abhinav Prakash Gahlot and Mathias Louboutin and Felix J. Herrmann}
}
@conference {bisbas2020IPDPStbf,
	title = {Temporal blocking of finite-difference stencil operators with sparse "off-the-grid" sources},
	booktitle = {IEEE International Parallel and Distributed Processing Symposium},
	year = {2021},
	note = {(IPDPS, virtual)},
	month = {5},
	abstract = {Stencil kernels dominate a range of scientific applications including seismic and medical imaging, image processing, and neural networks. Temporal blocking is a performance optimisation that aims to reduce the required memory bandwidth of stencil computations by re-using data from the cache for multiple time steps. It has already been shown to be beneficial for this class of algorithms. However, optimising stencils for practical applications remains challenging. These computations often include sparsely located operators, not aligned with the computational grid ("off-the-grid"). For example, our work is motivated by sources that inject a wavefield and measurements interpolating grid values. The resulting data dependencies make the adoption of temporal blocking much more challenging. We propose a methodology to inspect these data dependencies and reorder the computation, leading to performance gains in stencil codes where temporal blocking has not been applicable. We implement this novel scheme in the Devito domain-specific compiler toolchain. Devito implements a domain-specific language embedded in Python to generate optimised partial differential equation solvers using the finite-difference method from high-level symbolic problem definitions. We evaluate our scheme using isotropic acoustic, anisotropic acoustic and isotropic elastic wave propagators of industrial significance. Performance evaluation, after auto-tuning, shows that this enables substantial performance improvement through temporal blocking, over highly-optimised vectorized spatially-blocked code of up to 1.6x.},
	keywords = {finite-difference, HPC, performance, wave-equation},
	doi = {10.1109/IPDPS49936.2021.00058},
	url = {https://arxiv.org/pdf/2010.10248.pdf},
	author = {George Bisbas and Fabio Luporini and Mathias Louboutin and Rhodri Nelson and Gerard Gorman and Paul H. J. Kelly}
}
@conference {rizzuti2021SIAMGSwrid,
	title = {Time-domain Wavefield Reconstruction Inversion for large-scale seismic inversion},
	booktitle = {SIAM Conference on Mathematical and Computational Issues in the Geosciences},
	year = {2021},
	note = {(SIAM GS, virtual)},
	month = {06},
	abstract = {Data fitting based on a wave equation model has driven the modern development in seismic imaging, and it is now relevant for large scale inversion. A wider adoption is currently being limited by available computational resources and the classical limitations of non-linear optimization, which often produces suboptimal local minima. A class of methods based on the inexact solution of the wave equation, known as wavefield reconstruction inversion, have shown remarkable robustness towards local minima, however the solvers required for the inexact wave equation scale badly to 3D problems. We propose a reformulation of wavefield reconstruction inversion which involves the solution of the classical time-domain wave equation, and it is practical for large-scale problems. Applications to tilted-transverse isotropy synthetic datasets clearly demonstrate the advantage over traditional full-waveform inversion.},
	keywords = {FWI, SIAM, time domain, WRI},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMGS/2021/rizzuti2021SIAMGSwrid/rizzuti2021SIAMGSwrid.pdf},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Felix J. Herrmann}
}
@conference {rizzuti2020EAGEtwri,
	title = {Time-domain wavefield reconstruction inversion for large-scale seismics},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2020},
	note = {Accepted in EAGE, withdrawn due to logistical challenge},
	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}
}
@conference {louboutin2020SIAMTXLAtdw,
	title = {Time-Domain Wavefield Reconstruction Inversion In Tilted Transverse Isostropic Media},
	booktitle = {SIAM Texas-Louisiana},
	year = {2020},
	note = {SIAMTXLA},
	month = {10},
	abstract = {We introduce a generalization of time-domain wavefield reconstruction inversion to anisotropic acoustic modeling. Wavefield reconstruction inversion has been extensively researched in recent years for its ability to mitigate cycle skipping. The original method was formulated in the frequency domain with acoustic isotropic physics. However, frequency-domain modeling requires sophisticated iterative solvers that are difficult to scale to industrial-size problems and more realistic physical assumptions, such as tilted transverse isotropy, object of this study. The work presented here is based on a recently proposed dual formulation of wavefield reconstruction inversion, which allows time-domain propagator that are suitable to both large scales and more accurate physics.},
	keywords = {FWI, Time-domain, WRI},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMTXLA/2020/louboutin2020SIAMTXLAtdw/louboutin2020SIAMTXLAtdw_pres.pdf},
	author = {Mathias Louboutin and Gabrio Rizzuti and Felix J. Herrmann}
}
@conference {zhang2020SEGtli,
	title = {Transfer learning in large-scale ocean bottom seismic wavefield reconstruction},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2020},
	note = {(SEG, virtual)},
	month = {09},
	pages = {1666-1670},
	abstract = {Achieving desirable receiver sampling in ocean bottom acquisition is often not possible because of cost considerations. Assuming adequate source sampling is available, which is achievable by virtue of reciprocity and the use of modern randomized (simultaneous-source) marine acquisition technology, we are in a position to train convolutional neural networks (CNNs) to bring the receiver sampling to the same spatial grid as the dense source sampling. To accomplish this task, we form training pairs consisting of densely sampled data and artificially subsampled data using a reciprocity argument and the assumption that the source-site sampling is dense. While this approach has successfully been used on the recovery monochromatic frequency slices, its application in practice calls for wavefield reconstruction of time-domain data. Despite having the option to parallelize, the overall costs of this approach can become prohibitive if we decide to carry out the training and recovery independently for each frequency. Because different frequency slices share information, we propose the use the method of transfer training to make our approach computationally more efficient by warm starting the training with CNN weights obtained from a neighboring frequency slices. If the two neighboring frequency slices share information, we would expect the training to improve and converge faster. Our aim is to prove this principle by carrying a series of carefully selected experiments on a relatively large-scale five-dimensional data synthetic data volume associated with wide-azimuth 3D ocean bottom node acquisition. From these experiments, we observe that by transfer training we are able t significantly speedup in the training, specially at relatively higher frequencies where consecutive frequency slices are more correlated.},
	keywords = {SEG, transfer learning, wavefield reconstruction},
	doi = {10.1190/segam2020-3427882.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2020/zhang2020SEGtli/zhang2020SEGtli.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2020/zhang2020SEGtli/zhang2020SEGtli_pres.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2020/zhang2020SEGtli/zhang2020SEGtli_pres.mp4},
	software = {https://github.com/slimgroup/Software.SEG2020},
	author = {Mi Zhang and Ali Siahkoohi and Felix J. Herrmann}
}
@conference {yang2017SEGWStds,
	title = {Time domain sparsity promoting LSRTM with surface-related multiples in shallow-water case},
	booktitle = {SEG Workshop on Multiples: Status, Progress, Challenges and Open issues; Houston},
	year = {2017},
	note = {(SEG Workshop, Houston)},
	month = {09},
	abstract = {In the SRME relation, multiples are expressed as the multi- dimensional convolution between the the vertical derivative of the surface-free Green{\textquoteright}s function and the down-going receiver wavefield. This relation leads to the methods that separate the surface-related multiples and primaries. Instead of imaging separated multiples trivially, we introduced the SRME relation into wave equation by areal source injection, which costs nothing extra to involve multiples into the forward wavefields. The related image contains the phantoms components from cross-correlation between different-orders of events. Especially in shallow water, due to the strong magnitudes of multiples and their overlap with primaries, the artifacts contain strong phantoms not only of ocean bottom but also subsurface layer interfaces. We use sparsity-promoting inversions where the nicely curved structure of the subsurface models are detected to help cleaning up the phantoms. We reduce the costs by combin- ing randomized source subsampling with linearized Bregman method and get cleaned image with one data pass.},
	keywords = {multiple, SEG, workshop},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/yang2017SEGWStds/yang2017SEGWStds.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/yang2017SEGWStds/yang2017SEGWStds_post.pdf},
	author = {Mengmeng Yang and Felix J. Herrmann}
}
@conference {yang2016SEGtds,
	title = {Time-domain sparsity-promoting least-squares migration with source estimation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {4225-4229},
	abstract = {Traditional reverse-time migration (RTM) gives images with wrong amplitudes and low resolution. Least-squares RTM (LS-RTM) on the other hand, is capable of obtaining true-amplitude images as solutions of $\ell_2$-norm norm minimization problems by fitting the synthetic and observed reflection data. The shortcoming of this approach is that solutions of these $\ell_2$ problems are typically smoothed, tend to be overfitted, and computationally too expensive because it requires compared to standard RTM too many iterations. By working with randomized subsets of data only, the computational costs of LS-RTM can be brought down to an acceptable level while producing artifact-free high-resolution images without overfitting the data. While initial results of these "compressive imaging" methods were encouraging various open issues remain including guaranteed convergence, algorithmic complexity of the solver, and lack of on-the-fly source estimation for LS-RTMs with wave-equation solvers based on time-stepping. By including on-the-fly source-time function estimation into the method of Linearized Bregman (LB), on which we reported before, we tackle all these issues resulting in a easy-to-implement algorithm that offers flexibility in the trade-off between the number of iterations and the number of wave-equation solves per iteration for a fixed total number of wave-equation solves. Application of our algorithm on a 2D synthetic shows that we are able to obtain high-resolution images, including accurate estimates of the wavelet, for a single pass through the data. The produced image, which is by virtue of the inversion deconvolved with respect to the wavelet, is roughly of the same quality as the image obtained given the correct source function.},
	keywords = {least squares, RTM, SEG, source estimation, time domain},
	doi = {10.1190/segam2016-13850609.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/yang2016SEGtds/yang2016SEGtds.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/yang2016SEGtds/yang2016SEGtds_pres.pdf},
	author = {Mengmeng Yang and Philipp A. Witte and Zhilong Fang and Felix J. Herrmann}
}
@conference {kumar2016SEGtjm,
	title = {Time-jittered marine acquisition{\textendash}-a rank-minimization approach for 5D source separation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {119-123},
	abstract = {Simultaneous source marine acquisition has been recognized as an economic way of improving spatial sampling and speedup acquisition time, where a single- (or multiple-) source vessel fires at jittered source locations and time instances. Consequently, the acquired simultaneous data volume is processed to separate the overlapping shot records resulting in densely sampled data volume. It has been shown in the past that the simultaneous source acquisition design and source separation process can be setup as a compressed sensing problem, where conventional seismic data is reconstructed from simultaneous data via a sparsity-promoting optimization formulation. While the recovery quality of separated data is reasonably well, the recovery process can be computationally expensive due to transform-domain redundancy. In this paper, we present a computationally tractable rank-minimization algorithm to separate simultaneous data volumes. The proposed algorithm is suitable for large-scale seismic data, since it avoids singular-value decompositions and uses a low-rank based factorized formulation instead. Results are illustrated for simulations of simultaneous time-jittered continuous recording for a 3D ocean-bottom cable survey.},
	keywords = {5D, marine, SEG, source separation, time-jittered acquisition},
	doi = {10.1190/segam2016-13878249.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/kumar2016SEGtjm/kumar2016SEGtjm.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/kumar2016SEGtjm/kumar2016SEGtjm_pres.pdf},
	author = {Rajiv Kumar and Shashin Sharan and Haneet Wason and Felix J. Herrmann}
}
@conference {louboutin2015SEGtcs,
	title = {Time compressively sampled full-waveform inversion with stochastic optimization},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2015},
	note = {(SEG, New Orleans)},
	month = {10},
	pages = {5153-5157},
	abstract = {Time-domain Full-Waveform Inversion (FWI) aims to image the subsurface of the earth accurately from field recorded data and can be solved via the reduced adjoint-state method. However, this method requires access to the forward and adjoint wavefields that are meet when computing gradient updates. The challenge here is that the adjoint wavefield is computed in reverse order during time stepping and therefore requires storage or other type of mitigation because storing the full time history of the forward wavefield is too expensive in realistic 3D settings. To overcome this challenge, we propose an approximate adjoint-state method where the wavefields are subsampled randomly, which drastically the amount of storage needed. By using techniques from stochastic optimization, we control the errors induced by the subsampling. Examples of the proposed technique on a synthetic but realistic 2D model show that the subsampling-related artifacts can be reduced significantly by changing the sampling for each source after each model update. Combination of this gradient approximation with a quasi-Newton method shows virtually artifact free inversion results requiring only 5\% of storage compared to saving the history at Nyquist. In addition, we avoid having to recompute the wavefields as is required by checkpointing.},
	keywords = {acoustic, Full-waveform inversion, inversion, SEG, Stochastic optimization, Subsampling, Time-domain},
	doi = {10.1190/segam2015-5924937.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2015/louboutin2015SEGtcs/louboutin2015SEGtcs.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2015/louboutin2015SEGtcs/louboutin2015SEGtcs_poster.pdf},
	author = {Mathias Louboutin and Felix J. Herrmann}
}
@conference {witte2015IIPFWItdf,
	title = {Time-domain FWI in TTI media},
	booktitle = {Inaugural Full-Waveform Inversion Workshop},
	year = {2015},
	note = {(Natal, Brazil)},
	month = {08-09},
	abstract = {We develop an inversion workflow for tilted transverse isotropic (TTI) media using a purely acoustic formulation of the wave equation. The anisotropic modeling kernel is used for the forward modeling operator, as well as for the adjoint Jacobian to back propagate the data residual, thus providing the true gradient of the FWI objective function.},
	keywords = {FWI, time domain, TTI},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IIPFWI/witte2015IIPFWItdf/witte2015IIPFWItdf_pres.pdf},
	author = {Philipp A. Witte and Mathias Louboutin and Felix J. Herrmann}
}
@conference {oghenekohwo2015IIPFWItlf,
	title = {Time-lapse FWI with distributed compressed sensing},
	booktitle = {Inaugural Full-Waveform Inversion Workshop},
	year = {2015},
	note = {(Natal, Brazil)},
	month = {08-09},
	abstract = {In this talk, I will discuss our most recent application of full-waveform inversion (FWI) to time-lapse seismic data inversion. Specifically, I will illustrate how a joint recovery model (JRM) can be used in a joint inversion framework, leveraging ideas from distributed compressed sensing (DCS). The key idea involves exploiting the common information in the baseline and monitor gradients during the inversion.},
	keywords = {distributed compressed sensing, FWI, time lapse},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IIPFWI/oghenekohwo2015IIPFWItlf/oghenekohwo2015IIPFWItlf_pres.pdf},
	author = {Felix Oghenekohwo and Rajiv Kumar and Ernie Esser and Felix J. Herrmann}
}
@conference {oghenekohwo2014EAGEtls,
	title = {Time-lapse seismic without repetition: reaping the benefits from randomized sampling and joint recovery},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2014},
	month = {06},
	abstract = {In the current paradigm of 4-D seismic, guaranteeing repeatability in acquisition and processing of the baseline and monitor surveys ranks highest amongst the technical challenges one faces in detecting time-lapse signals. By using recent insights from the field of compressive sensing, we show that the condition of survey repeatability can be relaxed as long as we carry out a sparsitypromoting program that exploits shared information between the baseline and monitor surveys. By inverting for the baseline and monitor survey as the common "background", we are able to compute high-fidelity 4-D differences from carefully selected synthetic surveys that have different sets of source/receivers missing. This synthetic example is proof of concept of an exciting new approach to randomized 4-D acquisition where time-lapse signal can be computed as long as the survey details, such as source/receiver locations are known afterwards.},
	keywords = {4-D seismic, EAGE, joint recovery, time-lapse},
	doi = {10.3997/2214-4609.20141478},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/oghenekohwo2014EAGEtls/oghenekohwo2014EAGEtls.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2014/oghenekohwo2014EAGEtls/oghenekohwo2014EAGEtls_pres.pdf},
	author = {Felix Oghenekohwo and Ernie Esser and Felix J. Herrmann}
}
@conference {wason2013SEGtjo,
	title = {Time-jittered ocean bottom seismic acquisition},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {32},
	year = {2013},
	month = {09},
	pages = {1-6},
	abstract = {Leveraging ideas from the field of compressed sensing, we show how simultaneous or blended acquisition can be setup as a {\textendash} compressed sensing problem. This helps us to design a pragmatic time-jittered marine acquisition scheme where multiple source vessels sail across an ocean-bottom array firing airguns at {\textendash} jittered source locations and instances in time, resulting in better spatial sampling, and speedup acquisition. Furthermore, we can significantly impact the reconstruction quality of conventional seismic data (from jittered data) and demonstrate successful recovery by sparsity promotion. In contrast to random (under)sampling, acquisition via jittered (under)sampling helps in controlling the maximum gap size, which is a practical requirement of wavefield reconstruction with localized sparsifying transforms. Results are illustrated with simulations of time-jittered marine acquisition, which translates to jittered source locations for a given speed of the source vessel, for two source vessels.},
	keywords = {Acquisition, blending, deblending, interpolation, jittered sampling, marine, OBC, SEG},
	doi = {10.1190/segam2013-1391.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/wason2013SEGtjo/wason2013SEGtjo.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2013/wason2013SEGtjo/wason2013SEGtjo_pres.pdf},
	author = {Haneet Wason and Felix J. Herrmann}
}
@conference {hennenfent2004SEGtta,
	title = {Three-term amplitude-versus-offset (AVO) inversion revisited by curvelet and wavelet transforms},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2004},
	pages = {211-214},
	publisher = {SEG},
	organization = {SEG},
	abstract = {We present a new method to stabilize the three-term AVO inversion using Curvelet and Wavelet transforms. Curvelets are basis functions that effectively represent otherwise smooth objects having discontinuities along smooth curves. The applied formalism explores them to make the most of the continuity along reflectors in seismic images. Combined with Wavelets, Curvelets are used to denoise the data by penalizing high frequencies and small contributions in the AVO-cube. This approach is based on the idea that rapid amplitude changes along the ray-parameter axis are most likely due to noise. The AVO-inverse problem is linearized, formulated and solved for all (x, z) at once. Using densities and velocities of the Marmousi model to define the fluctuations in the elastic properties, the performance of the proposed method is studied and compared with the smoothing along the ray-parameter direction only. We show that our method better approximates the true data after the denoising step, especially when noise level increases. \copyright2004 Society of Exploration Geophysicists},
	keywords = {Presentation, SEG, SLIM},
	doi = {10.1190/1.1851201},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2004/Hennenfent04SEGtta/Hennenfent04SEGtta.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2004/Hennenfent04SEGtta/Hennenfent04SEGtta_pres.pdf},
	author = {Gilles Hennenfent and Felix J. Herrmann}
}
