@article {rizzuti2020dfw,
	title = {A dual formulation of wavefield reconstruction inversion for large-scale seismic inversion},
	journal = {Geophysics},
	volume = {86},
	number = {6},
	year = {2021},
	note = {(Geophysics)},
	month = {10},
	pages = {1ND-Z3},
	abstract = {Many of the seismic inversion techniques currently proposed that focus on robustness with respect to the background model choice are not apt to large-scale 3D applications, and the methods that are computationally feasible for industrial problems, such as full waveform inversion, are notoriously limited by convergence stagnation and require adequate starting models. We propose a novel solution that is both scalable and less sensitive to starting models or inaccurate parameters (such as anisotropy) that are typically kept fixed during inversion. It is based on a dual reformulation of the classical wavefield reconstruction inversion, whose empirical robustness with respect to these issues is well documented in the literature. While the classical version is not suited to 3D, as it leverages expensive frequency-domain solvers for the wave equation, our proposal allows the deployment of state-of-the-art time-domain finite-difference methods, and is potentially mature for industrial-scale problems.},
	keywords = {3D, finite-difference, Full-waveform inversion, wave equation},
	doi = {10.1190/geo2020-0743.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2021/rizzuti2020dfw/rizzuti2020dfw.html},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Felix J. Herrmann}
}
@presentation {rizzuti2021ML4SEISMICdfw,
	title = {A dual formulation of wavefield reconstruction inversion for large-scale seismic inversion},
	year = {2021},
	month = {11},
	abstract = {Many of the seismic inversion techniques currently proposed that focus on robustness with respect to the background model choice are not apt to large-scale 3D applications, and the methods that are computationally feasible for industrial problems, such as full waveform inversion, are notoriously limited by convergence stagnation and require adequate starting models. We propose a novel solution that is both scalable and less sensitive to starting models or inaccurate parameters (such as anisotropy) that are typically kept fixed during inversion. It is based on a dual reformulation of the classical wavefield reconstruction inversion, whose empirical robustness with respect to these issues is well documented in the literature. While the classical version is not suited to 3D, as it leverages expensive frequency-domain solvers for the wave equation, our proposal allows the deployment of state-of-the-art time-domain finite-difference methods, and is potentially mature for industrial-scale problems.},
	keywords = {FWI, ML4SEISMIC, SLIM, wave-equation, WRI},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ML4SEISMIC/2021/rizzuti2021ML4SEISMICdfw/Tue-10-55-Rizzuti.pdf},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Felix J. Herrmann}
}
@techreport {louboutin2021NIPSmcte,
	title = {Low-memory stochastic backpropagation with multi-channel randomized trace estimation},
	number = {TR-CSE-2021-1},
	year = {2021},
	month = {06},
	abstract = {Thanks to the combination of state-of-the-art accelerators and highly optimized open software frameworks, there has been tremendous progress in the performance of deep neural networks. While these developments have been responsible for many breakthroughs, progress towards solving large-scale problems, such as video encoding and semantic segmentation in 3D, is hampered because access to on-premise memory is often limited. Instead of relying on (optimal) checkpointing or invertibility of the network layers{\textendash}-to recover the activations during backpropagation{\textendash}-we propose to approximate the gradient of convolutional layers in neural networks with a multi-channel randomized trace estimation technique. Compared to other methods, this approach is simple, amenable to analyses, and leads to a greatly reduced memory footprint. Even though the randomized trace estimation introduces stochasticity during training, we argue that this is of little consequence as long as the induced errors are of the same order as errors in the gradient due to the use of stochastic gradient descent. We discuss the performance of networks trained with stochastic backpropagation and how the error can be controlled while maximizing memory usage and minimizing computational overhead.},
	keywords = {Convolutions, HPC, Low Memory, machine learning, randomized linear algebra},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2021/louboutin2021NIPSmcte/louboutin2021NIPSmcte.pdf},
	software = {https://github.com/slimgroup/XConv},
	author = {Mathias Louboutin and Ali Siahkoohi and Rongrong Wang and Felix J. Herrmann}
}
@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 {rizzuti2019SIAMGStwri,
	title = {A dual formulation for time-domain wavefield reconstruction inversion},
	booktitle = {SIAM Conference on Mathematical and Computational Issues in the Geosciences},
	year = {2019},
	note = {(SIAM GS)},
	month = {03},
	abstract = {Wavefield reconstruction inversion (WRI) is based on the combination of data misfit and PDE-residual penalty terms, and it aims at the reconstruction of the optimal property model and field pair. A full-space optimization would require the storage of field variables for every frequency or time sample (and every source). Remedies via variable projection are well studied, but they entail the solution of a modified wave equation, which is costly both in frequency and time domain. Here, we explore an efficient alternative to WRI, where properties are complemented, instead, by elements of the dual data space. The starting point of our proposal is the denoising reformulation of WRI: the PDE compliance is optimized under the constraint of a given data fit. Since the field variable can be eliminated from the associated saddle-point problem, we obtain a new objective functional of only properties and Lagrange multipliers. Its evaluation, and gradient computation, involves the solution of the original wave equation, which is amenable to time-domain discretization. The optimization scheme of choice can now leverage on affordable storage of the iterates.},
	keywords = {WRI},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SIAMGS/2019/rizzuti2019SIAMGStwri/rizzuti2019SIAMGStwri.pdf},
	author = {Gabrio Rizzuti and Mathias Louboutin and Rongrong Wang and Emmanouil Daskalakis and Felix J. Herrmann}
}
@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}
}
@article {sharan2018fsp,
	title = {Fast sparsity-promoting microseismic source estimation},
	journal = {Geophysical Journal International},
	volume = {216},
	number = {1},
	year = {2019},
	note = {(Geophysical Journal International)},
	month = {01},
	pages = {164-181},
	abstract = {Microseismic events are generated during hydraulic fracturing of unconventional reservoirs and carry information on fracture locations and the origin times associated with these microseismic events. For drilling purposes and to prevent hazardous situations, we need to have accurate knowledge on the fracture locations as well as on their size, and their density. Because microseismic waves can travel far distances, microseismic data collected at the surface and or in boreholes can help us to monitor hydraulic fracturing. While so-called back propagation or time-reversal methods are able to focus recorded energy back onto the sources when a reasonable velocity model is available, these methods suffer from blurring especially in situations where the data acquisition suffers from lack of aperture, sparse sampling, and noise. As a result, these methods typically cannot resolve sources in close proximity, a desired feature since we need this information if we want to follow the fracture evolution in space and time. In that situation, we need to estimate the locations and the associated source-time functions for closely spaced microseismic sources along the active fractures. To overcome the limitations of time-reversal methods, we propose a wave-equation based inversion approach where we invert for the complete source wavefield in both space and time. By promoting sparsity on the source wavefield in space, we negate the effects of non-radiating sources during the inversion and obtain high-resolution intensity plots and high-fidelity estimates for the source-time functions. We obtain these results relatively fast by accelerating the linearized Bregman method with a dual formulation. Through experiments, we demonstrate that our method is computationally feasible, robust to noise, and works for closely spaced sources with overlapping source-time functions in complex geological settings.},
	keywords = {induced seismicity, Joint inversion, waveform inversion},
	doi = {10.1093/gji/ggy415},
	url = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalJournalInternational/2018/sharan2018fsp/sharan2018fsp.pdf},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@conference {sharan2018SEGada,
	title = {A debiasing approach to microseismic},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2018},
	note = {(SEG, Anaheim)},
	month = {10},
	pages = {2942-2946},
	abstract = {Microseismic data is often used to locate fracture locations and their origin in time created by fracking. Although surface microseismic data can have large apertures and is easier to acquire than the borehole data, it often suffers from poor signal to noise ratio (S/R). Poor S/R poses a challenge in terms of estimating the correct location and source-time function of a microseismic source. In this work, we propose a denoising step in combination with a computationally cheap debiasing based approach to locate microseismic sources and to estimate their source-time functions with correct amplitude from extremely noisy data. Through numerical experiments, we demonstrate that our method can work with closely spaced microseismic sources with source-time functions of different peak amplitudes and frequencies. We have also shown the ability of our method with the smooth velocity model.},
	keywords = {Full-waveform inversion, Imaging, microseismic, noise, SEG, wave equation},
	doi = {10.1190/segam2018-2997252.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/sharan2018SEGada/sharan2018SEGada.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/sharan2018SEGada/sharan2018SEGada_pres.pdf},
	author = {Shashin Sharan and Rajiv Kumar and Rongrong Wang and Felix J. Herrmann}
}
@conference {yang2018SEGrde,
	title = {Removing density effects in LS-RTM with a low-rank matched filter},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2018},
	note = {(SEG, Anaheim)},
	month = {10},
	pages = {4251-4255},
	abstract = {Least-squares reverse-time migration faces difficulties when it inverts the data containing strong components related to density variation with velocity-only Born modeling operator. The strong density perturbation will be inverted as strong dummy velocity perturbations, which influence the amplitudes and phase of the velocity perturbations in the inverted model. The traditional method is to invert the additional density variations by developing Born operator with respect to both density and velocity or modify the image condition. In this work, we develop a matched-filter based LS-RTM for velocity-only Born modeling operator, which removes the artifacts in the imaging created by the strong density variation. This method doesn{\textquoteright}t call for extra work of finite difference stencil and is more general. In the experiment part, we use a complex discontinuous layered medium with strong density variations at the ocean bottom, and show the efficacy of the propose formulation.},
	keywords = {density effect, Least square, low rank, SEG},
	doi = {10.1190/segam2018-2997814.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/yang2018SEGrde/yang2018SEGrde.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/yang2018SEGrde/yang2018SEGrde_pres.pdf},
	author = {Mengmeng Yang and Rajiv Kumar and Rongrong Wang and Felix J. Herrmann}
}
@article {fang2017sewri,
	title = {Source estimation for wavefield-reconstruction inversion},
	journal = {Geophysics},
	volume = {83},
	number = {4},
	year = {2018},
	note = {(Geophysics)},
	pages = {R345-R359},
	abstract = {Source estimation is essential to all the wave-equation-based seismic inversions, including full-waveform inversion and the recently proposed wavefield-reconstruction inversion. When the source estimation is inaccurate, errors will propagate into the predicted data and introduce additional data misfit. As a consequence, inversion results that minimize this data misfit may become erroneous. To mitigate the errors introduced by the incorrect and pre-estimated sources, an embedded procedure that updates sources along with medium parameters is necessary for the inversion. So far, such a procedure is still missing in the context of wavefield-reconstruction inversion, a method that is, in many situations, less prone to local minima related to the so-called cycle skipping, compared to full-waveform inversion through exact data-fitting. While wavefield-reconstruction inversion indeed helps to mitigate issues related to cycle skipping by extending the search space with wavefields as auxiliary variables, it relies on having access to the correct source functions. In this paper, we remove the requirement of having the accurate source functions by proposing a source estimation technique specifically designed for wavefield-reconstruction inversion. To achieve this task, we consider the source functions as unknown variables and arrive at an objective function that depends on the medium parameters, wavefields, and source functions. During each iteration, we apply the so-called variable projection method to simultaneously project out the source functions and wavefields. After the projection, we obtain a reduced objective function that only depends on the medium parameters and invert for the unknown medium parameters by minimizing this reduced objective. Numerical experiments illustrate that this approach can produce accurate estimates of the unknown medium parameters without any prior information of the source functions.},
	keywords = {estimation, FWI, source, variable projection, WRI},
	doi = {10.1190/geo2017-0700.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2018/fang2017sewri/fang2017sewri.html},
	author = {Zhilong Fang and Rongrong Wang and Felix J. Herrmann}
}
@conference {sharan2018IEEEIUSspoa,
	title = {Sparsity-promoting photoacoustic imaging with source estimation},
	booktitle = {2018 IEEE International Ultrasonics Symposium (IUS)},
	year = {2018},
	note = {(IEEE IUS, Kobe)},
	month = {10},
	pages = {206-212},
	abstract = {Photoacoustics has emerged as a high-contrast imaging modality that provides optical absorption maps inside of tissues, therefore complementing morphological information of conventional ultrasound. The laser-generated photoacoustic waves are usually envelope-detected, thus disregarding the specific waveforms generated by each photoabsorber. Here we propose a sparsity-promoting image reconstruction method that allows the estimation of each photoabsorber{\textquoteright}s source-time function. Preliminary studies showed the ability to reconstruct the optical absorption map of an in silico vessel phantom. By using a sparsity-promoting imaging method, absorption maps and source-time functions can still be recovered even in situations where the number of transducers is decreased by a factor of six. Moreover, the recovery is able to attain higher resolution than conventional beamforming methods. Because our method recovers the source-time function of the absorbers, it could potentially also be used to distinguish different types of photoabsorbers, or the degree of aggregation of exogenous agents, under the assumption that these would generate different source-time functions at the moment of laser irradiation.},
	keywords = {Acceleration, estimation, IEEE, Image reconstruction, Phantoms, photoacoustic imaging, source-time function, sparsity promoting inversion, Transducers, wave equation},
	doi = {10.1109/ULTSYM.2018.8580037},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/IEEEIUS/2018/sharan2018IEEEIUSspoa/sharan2018IEEEIUSspoa.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IEEEIUS/2018/sharan2018IEEEIUSspoa/sharan2018IEEEIUSspoa_pres.pdf},
	author = {Shashin Sharan and Rajiv Kumar and Diego S. Dumani and Mathias Louboutin and Rongrong Wang and Stanislav Emelianov and Felix J. Herrmann}
}
@conference {wang2017SEGdff,
	title = {A denoising formulation of full-waveform inversion},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2017},
	note = {(SEG, Houston)},
	month = {09},
	pages = {1594-1598},
	abstract = {We propose a wave-equation-based subsurface inversion method that in many cases is more robust than the conventional Full-Waveform Inversion. The new formulation is written in a denoising form that allows the synthetic data to match the observed ones up to a small error. Compared to the Full-Waveform Inversion, our method treats the noise arising from the data measuring/recording process and that from the synthetic modelling process separately. Comparing to the Wavefields Reconstruction Inversion, the new formulation mitigates the difficulty of choosing the penalty parameter λ. To solve the proposed optimization problem, we develop an efficient frequency domain algorithm that alternatively updates the model and the data. Numerical experiments confirm strong stability of the proposed method by comparisons between the results of our algorithm with that from both plain FWI and a weighted formulation of the FWI.},
	keywords = {denoising, FWI, SEG},
	doi = {10.1190/segam2017-17794690.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/wang2017SEGdff/wang2017SEGdff.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/wang2017SEGdff/wang2017SEGdff_pres.pdf},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@conference {sharan2017SEGhrf,
	title = {High resolution fast microseismic source collocation and source time function estimation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2017},
	note = {(SEG, Houston)},
	month = {09},
	pages = {2778-2783},
	abstract = {Sparsity promotion based joint microseismic source collocation and source time function estimation, using Linearized Bregman algorithm, although simple to implement, suffers from slow convergence. This is due to the fact that Linearized Bregman algorithm has only first order of convergence. In this work, we propose to accelerate the existing Linearized Bregman algorithm using the L-BFGS algorithm. Without any initial guess for the source location or source time function, our method is able to estimate the source location and source time function for kinematically correct smooth velocity model. We demonstrate the effectiveness of our method for multiple sources spaced within half a wavelength. We also compare our results with Linearized Bregman based method in {\textquoteleft}{\textquoteleft}2.5{\textquoteleft}{\textquoteleft}D instead of {\textquoteleft}{\textquoteleft}2{\textquoteleft}{\textquoteleft}D.},
	keywords = {algorithm, Full-waveform inversion, microseismic, passive imaging, SEG, wave equation},
	doi = {10.1190/segam2017-17787749.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/sharan2017SEGhrf/sharan2017SEGhrf.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/sharan2017SEGhrf/sharan2017SEGhrf_pres.pdf},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@presentation {wang2017SINBADFnra,
	title = {Noise robust and time-domain formulations of Wavefield Reconstruction Inversion},
	year = {2017},
	publisher = {SINBAD},
	abstract = {We propose a wave-equation-based subsurface inversion method that in many cases is more robust than conventional Full-Waveform Inversion. The new formulation is written in a denoising form that allows the synthetic data to match the observed ones up to a small error. Compared to regular Full-Waveform Inversion, our method treats the noise arising from the data meassuring/recording process and that from the synthetic modelling process separately. Compared to Wavefields Reconstruction Inversion, the new formulation mitigates the difficulty of choosing the penalty parameter λ. To solve the proposed optimization problem, we develop an efficient frequency domain algorithm that alternatively updates the model and the data. Numerical experiments confirm strong stability of the proposed method by comparisons between the results of our algorithm with that from both plain FWI and a weighted formulation of the FWI. We also discuss a new memory efficient time-domain formulation for Wavefield Reconstruction Inversion based on duality.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/wang2017SINBADFnra/wang2017SINBADFnra.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/wang2017SINBADFnra/wang2017SINBADFnra.mov},
	author = {Rongrong Wang and Mathias Louboutin and Bas Peters and Emmanouil Daskalakis and Felix J. Herrmann}
}
@presentation {sharan2017SINBADFtts,
	title = {Tracking the spatial-temporal evolution of fractures by microseismic source collocation},
	year = {2017},
	publisher = {SINBAD},
	abstract = {Unlike conventional reservoirs, unconventional plays are not naturally viable for economical production of oil and gas. They require stimulation by injecting high-pressure fluid causing fractures in the rocks. These fractures make the medium more permeable, hence, the extraction of oil and gas becomes feasible. For drilling purposes and to prevent potentially hazardous situations, we need to have good knowledge of the location of these fractures. Also, we need to have good knowledge about how these fractures originated in time. Hydraulic fracturing changes stress in rocks, which results in the emission of microseismic waves. The opening of cracks due to high pressure fluid injection during hydraulic fracturing mainly causes this change in stress in the rocks. Therefore, microseismic events are mostly localized along these fractures and have finite energy along time. To accurately track the evolution of fractures in both space and time, we need to locate closely spaced microseismic events along these fractures activating at very small time intervals. A naive approach can be the back propagation of the observed data to find out a point in space and time where maximum focusing of back propagating energy occurs. This point corresponds to the location and origin time of a microseismic source. This approach, although simpler, suffers from low resolution and requires scanning of complete 4D volume (3D in space and 1D in time). Hence, this method can be challenging when there are multiple closely spaced microseismic sources originating at different times. We in this work propose a sparsity promotion based method that can locate closely spaced microseismic events, with spatial separation as low as within half a wavelength, activating at small time intervals. We simultaneously estimate the origin time of microseismic events by estimating their source time functions. Our method exploits the fact that microseismic events are localized in space and have finite energy. We use accelerated Linearized Bregman algorithm with a preconditioning operator to arrive at a computationally feasible scheme.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/sharan2017SINBADFtts/sharan2017SINBADFtts.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/sharan2017SINBADFtts/sharan2017SINBADFtts.mov},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@conference {wang2016SEGfde,
	title = {Frequency down-extrapolation with TV norm minimization},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {1380-1384},
	abstract = {In this work, we present a total-variation (TV)-norm minimization method to perform model-free low-frequency data extrapolation for the purpose of assisting full-waveform inversion. Low-frequency extrapolation is the process of extending reliable frequency bands of the raw data towards the lower end of the spectrum. To this end, we propose to solve a TV-norm based convex optimization problem that has a global minimum and is equipped with a fast solver. The approach takes into account both the sparsity of the reflectivity series associated with a single trace, as well as the inter-trace correlations. A favorable byproduct of this approach is that it allows one to work with coarsely sampled trace along time (as coarse as 0.02s per sample), hence substantially reduces the size of the proposed optimization problem. Last, we show the effectiveness of the method for frequency-domain FWI on the Marmousi model.},
	keywords = {extrapolation, FWI, SEG, sparse},
	doi = {10.1190/segam2016-13879674.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/wang2016SEGfde/wang2016SEGfde.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/wang2016SEGfde/wang2016SEGfde_pres.pdf},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@presentation {sharan2016SINBADFhrm,
	title = {High resolution microseismic source collocation},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Microseismic waves generated by hydraulic fracturing are used to locate the fractures created in the reservoir. Most of the existing source collocation methods have their limitations in terms of resolution of very closely spaced microseismic events. In this work we propose a method to collocate microseismic events with high resolution along with simultaneous estimation of source time function. In the context of waves, the finite-difference modelling operator acts as an analysis operator for wavefields. This implies that wavefields under the action of a finite-difference modelling kernel can be focused to sources generating these wavefields. We exploit this property of wavefields to locate microseismic events and estimate the source time function jointly. To arrive at a computational feasible scheme, we use the accelerated version of Linearized Bregman (LBR) algorithm to solve the above mentioned problem. This is joint work with Rongrong Wang.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/sharan2016SINBADFhrm/sharan2016SINBADFhrm.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/sharan2016SINBADFhrm/sharan2016SINBADFhrm.mov},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@conference {chai2016SEGlbm,
	title = {A linearized Bregman method for compressive waveform inversion},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {1449-1454},
	abstract = {We present an implementation of a recently developed and relatively simple linearized Bregman method to solve the large-scale $\ell_1$-norm sparsity-promoting Gauss-Newton (GN) full-waveform inversion (FWI) problem. Numerical experiments demonstrate that: the much simpler linearized Bregman method does a comparable and even superior job compared to a state-of-the-art $\ell_1$-norm solver SPGL1, which is previously used in the modified Gauss-Newton FWI; the linearized Bregman method is also more efficient (faster) than SPGL1; the FWI result with the linearized Bregman method solving $\ell_1$-norm sparsity-promoting problems to get the model updates is better than that obtained by solving $\ell_2$-norm least-squares problems to get the model updates.},
	keywords = {Bregman method, SEG, waveform inversion},
	doi = {10.1190/segam2016-13848105.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/chai2016SEGlbm/chai2016SEGlbm.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/chai2016SEGlbm/chai2016SEGlbm_pres.pdf},
	author = {Xintao Chai and Mengmeng Yang and Philipp A. Witte and Rongrong Wang and Zhilong Fang and Felix J. Herrmann}
}
@conference {sharan2016SEGspj,
	title = {Sparsity-promoting joint microseismic source collocation and source-time function estimation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2016},
	note = {(SEG, Dallas)},
	month = {10},
	pages = {2574-2579},
	abstract = {In this work, we propose a new method to simultaneously locate microseismic events (e.g., induced by hydraulic fracturing) and estimate the source signature of these events. We use the method of linearized Bregman. This algorithm focuses unknown sources at their true locations by promoting sparsity along space and at the same time keeping the energy along time in check. We are particularly interested in situations where the microseismic data is noisy, sources have different signatures and we only have access to the smooth background-velocity model. We perform numerical experiments to demonstrate the usability of the proposed method. We also compare our results with full-waveform inversion based microseismic event collocation methods. Our method gives flexibility to simultaneously get a more accurate source image along with an estimate of the source-time function, which carries important information on the rupturing process and source mechanism.},
	keywords = {microseismic, noise, SEG, sources, sparse, wave equation},
	doi = {10.1190/segam2016-13871022.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/sharan2016SEGspj/sharan2016SEGspj.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2016/sharan2016SEGspj/sharan2016SEGspj_pres.pdf},
	author = {Shashin Sharan and Rongrong Wang and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {wang2016SINBADFtmf,
	title = {Two methods for frequency down-extrapolation},
	year = {2016},
	publisher = {SINBAD},
	abstract = {It is well known that low frequency data are essential to the success of the Full-Waveform Inversion. Field data from a usual marine acquisition, however, typically lack reliable low frequencies due to the existing physical limitations. One research direction focuses on studying low-frequency extrapolation, the process of extending reliable frequency bands of the raw data towards the lower end of the spectrum. We propose two optimization problems, by solving which low frequencies can be made out of the high SNR high frequency component. We present both theoretical explanation and numerical experiments to demonstrate the efficacy as well as limitation of these methods.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/wang2016SINBADFtmf/wang2016SINBADFtmf.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/wang2016SINBADFtmf/wang2016SINBADFtmf.mov},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@presentation {wang2015SINBADFifw,
	title = {Improving full waveform inversion with spectral extrapolation},
	year = {2015},
	publisher = {SINBAD},
	abstract = {We develop a sparsity based frequency extrapolation method to obtain low frequency data from the high frequency ones. As other frequency extrapolation methods, the extrapolation introduces noise due to dispersion and possibly densely distributed reflectors. To mitigate such noise, we incorporate Curvlet coefficient regularization in the extrapolation algorithm. Numerical results demonstrate the effectiveness of the extrapolation in frequency domain FWI in models of moderate complexity.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/wang2015SINBADFifw/wang2015SINBADFifw.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/wang2015SINBADFifw/wang2015SINBADFifw.mov},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@conference {esser2015EAGElcs,
	title = {A lifted $\ell_1/\ell_2$ constraint for sparse blind deconvolution},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2015},
	note = {(EAGE, Madrid)},
	month = {06},
	abstract = {We propose a modification to a sparsity constraint based on the ratio of $\ell_1$ and $\ell_2$ norms for solving blind seismic deconvolution problems in which the data consist of linear convolutions of different sparse reflectivities with the same source wavelet. We also extend the approach to the Estimation of Primaries by Sparse Inversion (EPSI) model, which includes surface related multiples. Minimizing the ratio of $\ell_1$ and $\ell_2$ norms has been previously shown to promote sparsity in a variety of applications including blind deconvolution. Most existing implementations are heuristic or require smoothing the $\ell_1/\ell_2$ penalty. Lifted versions of $\ell_1/\ell_2$ constraints have also been proposed but are challenging to implement. Inspired by the lifting approach, we propose to split the sparse signals into positive and negative components and apply an $\ell_1/\ell_2$ constraint to the difference, thereby obtaining a constraint that is easy to implement without smoothing the $\ell_1$ or $\ell_2$ norms. We show that a method of multipliers implementation of the resulting model can recover source wavelets that are not necessarily minimum phase and approximately reconstruct the sparse reflectivities. Numerical experiments demonstrate robustness to the initialization as well as to noise in the data.},
	keywords = {blind deconvolution, EAGE, EPSI},
	doi = {10.3997/2214-4609.201413420},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/esser2015EAGElcs/esser2015EAGElcs.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/esser2015EAGElcs/esser2015EAGElcs_poster.pdf},
	author = {Ernie Esser and Tim T.Y. Lin and Rongrong Wang and Felix J. Herrmann}
}
@conference {esser2015CAMSAPrsa,
	title = {Resolving scaling ambiguities with the $\ell_1/\ell_2$ norm in a blind deconvolution problem with feedback},
	booktitle = {IEEE International Workshop on Computational Advances in Multi-Sensor Adaptive Processing},
	year = {2015},
	note = {(IEEE CAMSAP Workshop, Canc{\'u}n, Mexico)},
	month = {12},
	pages = {365-368},
	abstract = {Compared to more mundane blind deconvolution problems, blind deconvolution in seismic applications involves a feedback mechanism related to the free surface. The presence of this feedback mechanism gives us an unique opportunity to remove ambiguities that have plagued blind deconvolution for a long time. While beneficial, this feedback by itself is insufficient to remove the ambiguities even with $\ell_1$ constraints. However, when paired with an $\ell_1/\ell_2$ constraint the feedback allows us to resolve the scaling ambiguity under relatively mild assumptions. Inspired by lifting approaches, we propose to split the sparse signal into positive and negative components and apply an $\ell_1/\ell_2$ constraint to the difference, thereby obtaining a constraint that is easy to implement. Numerical experiments demonstrate robustness to the initialization as well as to noise in the data.},
	keywords = {$\ell_1/\ell_2$, blind deconvolution, CAMSAP, lifting, method of multipliers},
	doi = {10.1109/CAMSAP.2015.7383812},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/CAMSAP/2015/esser2015CAMSAPrsa/esser2015CAMSAPrsa.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/CAMSAP/2015/esser2015CAMSAPrsa/esser2015CAMSAPrsa_poster.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/CAMSAP/2015/esser2015CAMSAPrsa/esser2015CAMSAPrsa_ext.pdf},
	author = {Ernie Esser and Rongrong Wang and Tim T.Y. Lin and Felix J. Herrmann}
}
@presentation {esser2015SINBADFrsa,
	title = {Resolving scaling ambiguities with the L1/L2 norm in a blind deconvolution problem with feedback},
	year = {2015},
	publisher = {SINBAD},
	abstract = {Compared to more mundane blind deconvolution problems, blind deconvolution in seismic applications involves a feedback mechanism related to the free surface. The presence of this feedback mechanism gives us an unique opportunity to remove ambiguities that have plagued blind deconvolution for a long time. While beneficial, this feedback by itself is insufficient to remove the ambiguities even with L1 constraints. However, when paired with an L1/L2 constraint the feedback allows us to resolve the scaling ambiguity under relatively mild assumptions. Inspired by lifting approaches, we propose to split the sparse signal into positive and negative components and apply an \ell_1/\ell_2 constraint to the difference, thereby obtaining a constraint that is easy to implement. Numerical experiments demonstrate robustness to the initialization as well as to noise in the data.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/esser2015SINBADFrsa/esser2015SINBADFrsa.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/esser2015SINBADFrsa/esser2015SINBADFrsa.mov},
	author = {Ernie Esser and Rongrong Wang and Tim T.Y. Lin and Felix J. Herrmann}
}
@presentation {sharan2015SINBADFscu,
	title = {Source collocation using the method of Linearized Bregman},
	year = {2015},
	publisher = {SINBAD},
	abstract = {In this work, we present a method to collocate sources from seismic reflections recorded at receiver locations. We use the method of linearized Bregman. This algorithm focuses unknown sources by promoting sparsity in the {\textquotedblleft}Helmholtz domain{\textquotedblright}{\textemdash}i.e., the wavefield under the action of the the Helmholtz system{\textemdash}under the constraint of a data misfit within a particular tolerance. We extend this method to noisy measurements. We are in particular interested in situations where the noise is coherent because it is in the range of the Helmholtz equation.},
	keywords = {Presentation, SINBAD, SINBADFALL2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/sharan2015SINBADFscu/sharan2015SINBADFscu.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Fall/sharan2015SINBADFscu/sharan2015SINBADFscu.mov},
	author = {Shashin Sharan and Rongrong Wang and Felix J. Herrmann}
}
@conference {wang2015SIAMwds,
	title = {Wavefield-denoising and source encoding},
	booktitle = {SIAM Conference on Mathematical and Computational Issues in the Geosciences},
	year = {2015},
	note = {(SIAM Conference on Mathematical and Computational Issues in the Geosciences, Stanford University, California)},
	month = {06-07},
	keywords = {SIAM, source encoding, WRI},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SIAM/2015/wang2015SIAMwds/wang2015SIAMwds_pres.pdf},
	author = {Rongrong Wang and Ozgur Yilmaz and Felix J. Herrmann}
}
@presentation {wang2014SINBADdwi,
	title = {Denoising the wavefield inversion problem through source blending and penalty method},
	year = {2014},
	publisher = {SINBAD},
	abstract = {Previous FWI denoising techniques, such as source blending and stacking, are designed applied to the seismic traces. Hence no physical information is used to improve the fidelity of the reconstruction, even when a good initial model is available. Moreover, traditional source blending requires the existence of common source gathers, which might be missing from the data. In this talk, we incorporate the source blending idea into the penalty method, and extract the synthetic sources corresponding to the principle directions of the gradient update. By blending the sources, we can not only achieve an acceleration and storage efficiency of the penalty method, but also make the reconstruction more robust to the random noise. This method is applicable to any source receiver position but is more efficient if common source gathers are available.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/wang2014SINBADdwi/wang2014SINBADdwi.pdf},
	author = {Rongrong Wang and Felix J. Herrmann}
}
@techreport {wang2014SEGfwi,
	title = {Full waveform inversion with interferometric measurements},
	number = {TR-EOAS-2014-5},
	year = {2014},
	month = {04},
	institution = {UBC},
	abstract = {In this note, we design new misfit functions for full-waveform inversion by using interferometric measurements to reduce sensitivity to phase errors. Though established within a completely different setting from the linear case, we obtain a similar observation: the interferometry can improve robustness under certain modeling errors. Moreover, in order to deal with errors on both source and receiver sides, we propose a higher order interferometry, which, as a generalization of the usual definition, involves the cross correlation of four traces. A proof of principle simulations is included on a stylized example.},
	keywords = {FWI},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2014/wang2014SEGfwi/wang2014SEGfwi.html},
	author = {Rongrong Wang and Ozgur Yilmaz and Felix J. Herrmann}
}
@presentation {esser2014SINBADlcs,
	title = {A lifted $\ell_1$/$\ell_2$ constraint for sparse blind deconvolution},
	year = {2014},
	publisher = {SINBAD},
	abstract = {We propose a modification to a sparsity constraint based on the ratio of l1 and l2 norms for solving blind seismic deconvolution problems in which the data consist of linear convolutions of different sparse reflectivities with the same source wavelet. No assumptions are made about the location of the support of either the wavelet or the sparse signals. Minimizing the ratio of l1 and l2 norms has been previously shown to promote sparsity in a variety of applications including blind deconvolution. Most existing implementations are heuristic or require smoothing the l1/l2 penalty. Lifted versions of l1/l2 constraints have also been proposed but are still challenging to implement. Inspired by the lifting approach, we propose to split the sparse signals into positive and negative components and apply an l1/l2 constraint to the difference, thereby obtaining a constraint that is easy to implement without smoothing the l1 or l2 norms. We show that a method of multipliers implementation of the resulting model can recover source wavelets that are not necessarily minimum phase and approximately reconstruct the sparse reflectivities. It appears to be robust to the initialization and to small amounts of noise in the data. We also discuss extensions to the Estimation of Primaries by Sparse Inversion (EPSI) convolution model.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/esser2014SINBADlcs/esser2014SINBADlcs.pdf},
	author = {Ernie Esser and Tim T.Y. Lin and Rongrong Wang and Felix J. Herrmann}
}
@presentation {wang2013SINBADnru,
	title = {Noise reduction by using interferometric measurements},
	year = {2013},
	publisher = {SINBAD},
	abstract = {The interferometric formulation of linear wave-based inversion problem was proposed by Demanet and Jugnon recently. Instead of directly fitting the data, they proposed to fit a subset of the data{\textquoteright}s cross-correlation. It can be verified that if the full cross-correlation is used, then the problem is equivalent to the usual least square problem. The subsampling, which is usually considered to cause instability to the solution, is surprisingly useful in this setting. Numerical experiments for the inverse source problem and the inverse scatting problem have both suggested that a {\textquoteright}good{\textquoteright} sampling strategy can actually increase the stability under modeling error caused by the uncertainty of a kinetic nature. We will study the mathematical mechanism behind this phenomenon, and try to see whether or not there exists a universally {\textquoteright}good{\textquoteright} sampling strategy independent of the types of forward operators we use.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/wang2013SINBADnru/wang2013SINBADnru.pdf},
	author = {Rongrong Wang and Bas Peters and Felix J. Herrmann}
}
