@article {yang2020lrpo,
	title = {Low-rank representation of omnidirectional subsurface extended image volumes},
	journal = {Geophysics},
	volume = {86},
	number = {3},
	year = {2021},
	note = {(Geophysics)},
	month = {1},
	pages = {1MJ-WA152},
	abstract = {Subsurface-offset gathers play an increasingly important role in seismic imaging. These gathers are used during velocity model building and inversion of rock properties from amplitude variations. While powerful, these gathers come with high computational and storage demands to form and manipulate these high dimensional objects. This explains why only limited numbers of image gathers are computed over a limited offset range. We avoid these high costs by working with highly compressed low-rank factorizations. We arrive at these factorizations via a combination of probings with the double two-way wave equation and randomized singular value decompositions. In turn, the resulting factorizations give us access to all subsurface offsets without having to form the full extended image volumes. The latter is computationally prohibitive because extended image volumes are quadratic in image size. As a result, we can easily handle situations where conventional horizontal offset gathers are no longer focused. More importantly, the factorization also provides a mechanism to use the invariance relation of extended image volumes for velocity continuation. With this technique, extended image volumes for one background velocity model can directly be mapped to those of another background velocity model. Our low-rank factorization inherits this invariance property so we incur factorization costs only once when examining different imaging scenarios. Because all imaging experiments only involve the factors, they are computationally cheap with costs that scale with the rank of the factorization. We validate our methodology on 2D synthetics including a challenging imaging example with salt. Our experiments show that our low-rank factorization parameterizes extended image volumes naturally. Instead of brute force explicit cross-correlations between shifted source and receiver wavefields, our approach relies on the underlying linear-algebra structure that enables us to work with these objects without incurring unfeasible demands on computation and storage.},
	keywords = {extended image volumes, invariance relationship, low rank, power schemes, randomized linear algebra},
	doi = {10.1190/geo2020-0152.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2021/yang2020lrpo/Paper_final.html},
	author = {Mengmeng Yang and Marie Graff and Rajiv Kumar and Felix J. Herrmann}
}
@conference {yang2019SEGlrr,
	title = {Low-rank representation of subsurface extended image volumes with power iterations},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2019},
	note = {(SEG, San Antonio)},
	month = {09},
	pages = {4470-4474},
	abstract = {Full subsurface extended image volumes (EIVs) contain abundant information such as subsurface image gathers used in imaging, interpretation of rock properties or velocity analysis, but are expensive in term of computational cost and storage requirement due to their large size. Yet, due to their redundant information, monochromatic EIVs exhibit a low-rank structure that allows us to get a good approximation at a computational cost proportional to the rank k , while conventional techniques require at least the number of sources. Such low-rank factorized EIVs are computed thanks to the randomized singular values decomposition (SVD) algorithm. However, recent developments on low-rank approximations of EIVs raise two major questions. First, monochromatic EIVs rely on time-harmonic wave equation solvers, which do not scale well to realistic 3D models. Second, the rank of the monochromatic EIVs increases with the frequency, which yields increasing computational costs and storage. Here, we propose a construction approach based on a time-domain finite-difference wave-equation solver with time stepping, which is combined to power iteration schemes in the randomized SVD algorithm to accelerate the decay of the singular values. We compare the performances of simultaneous iterations, block-Krylov iterations and the original randomized SVD on the computation of the EIV for a small section of the Marmousi model. Then, we show further results on the full Marmousi model to validate our approach.},
	keywords = {full extended image volumes, power iterations, randomized SVD, SEG, time domain},
	doi = {10.1190/segam2019-3215961.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/yang2019SEGlrr/yang2019SEGlrr.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2019/yang2019SEGlrr/yang2019SEGlrr_pres.pdf},
	author = {Mengmeng Yang and Marie Graff and Rajiv Kumar and Felix J. Herrmann}
}
@conference {kumar2018SEGlreiv,
	title = {Low-rank representation of extended image volumes{\textendash}{\textendash}applications to imaging and velocity continuation},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2018},
	note = {(SEG, Anaheim)},
	month = {10},
	pages = {4418-4422},
	abstract = {In this paper, we present a cost-effective low-rank factorized form of computing the full subsurface extended image volumes without explicitly calculating the receiver wavefields. The propose approach is computationally feasible, which exploits the low-rank structure of full subsurface extended image volumes organized as a matrix, thus avoiding the customary loop over sources. Using carefully selected stylized examples, we show how conventional migration images can be extracted from the low-rank factorized form. We also present a velocity continuation procedure that uses the same low-rank form and allows us to map the extended image for one velocity model to another velocity model without recomputing all the source wavefields.},
	keywords = {image volume, low-rank, SEG, velocity continuation},
	doi = {10.1190/segam2018-2998404.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2018/kumar2018SEGlreiv/kumar2018SEGlreiv.html},
	author = {Rajiv Kumar and Marie Graff-Kray and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {graff2017SINBADFlrp,
	title = {Low-rank representation of omnidirectional subsurface extended image volumes},
	year = {2017},
	publisher = {SINBAD},
	abstract = {Extended image volumes are an important migration tool in seismic exploration. However the computation and the storage of omnidirectional subsurface extended image volumes are usually prohibitive. That is why some solutions have been already proposed for instance by focusing on horizontal offsets only. In our work, we will consider a linear algebra approach to deal with the low-rank representation of extended image volumes with full offsets. We will never build entirely the resulting matrix but get only actions of it on well-chosen probing vectors, based on Low-Rank decomposition or randomized SVD. This representation allows us to have access to all the energy of the extended image volume matrix and still limits the storage of the information and the computational cost.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/graff2017SINBADFlrp/graff2017SINBADFlrp.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/graff2017SINBADFlrp/graff2017SINBADFlrp.mov},
	author = {Marie Graff-Kray and Rajiv Kumar and Felix J. Herrmann}
}
@presentation {wason2016SINBADFlcr,
	title = {Low-cost, randomized 3D towed-marine time-lapse seismic acquisition},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Most often conventional 3D towed-marine seismic data acquisition suffers from limited illumination of subsurface due to narrow azimuth. Although, acquisition techniques such as Ocean Bottom Cable (OBC), Ocean Bottom Node (OBN), Rich-Azimuth (RAZ) acquisition, Multi-Azimuth (MAZ) acquisition have been developed to illuminate the subsurface from all possible angles, these acquisition techniques are expensive. Motivated by the design principles of compressed sensing, we acquire randomly subsampled and simultaneous (e.g., simultaneous long offset) 3D towed-marine time-lapse seismic data for few (random) azimuths, and recover densely sampled multiple azimuth baseline and monitor data via sparsity promotion and rank minimization. Our acquisition is low cost since we have subsampled measurements.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/wason2016SINBADFlcr/wason2016SINBADFlcr.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/wason2016SINBADFlcr/marine_acq_mask.pdf},
	author = {Rajiv Kumar and Felix Oghenekohwo and Shashin Sharan and Haneet Wason and Felix J. Herrmann}
}
@presentation {lopez2016SINBADFlrm,
	title = {Low-rank matrix recovery for parallel architectures},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Low-rank matrix recovery (LRMR) techniques offer potential tools for frugal seismic data acquisition, where dense acquisition is replaced by large-scale optimization. This shift of focus means that efficient numerical methods are critical for successful LRMR implementation in seismic applications. In this talk, we extend the rank-penalization methodology to a parallelizable framework. We adopt a factorization-based alternating minimization scheme and decouple it into an independent system of simpler sub-problems that can be handled in parallel. The methodology is flexible, where the approach can be adapted to the number of workers available. Numerical experiments are conducted to demonstrate that the quality of reconstruction is comparable to existing methods at a fraction of the computational time.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/lopez2016SINBADFlrm/lopez2016SINBADFlrm.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/lopez2016SINBADFlrm/lopez2016SINBADFlrm.mov},
	author = {Oscar Lopez and Rajiv Kumar and Felix J. Herrmann}
}
@presentation {kumar2016SINBADFlrm,
	title = {Low-rank methods for on-the-fly slicing \& dicing of seismic data \& image volumes},
	year = {2016},
	publisher = {SINBAD},
	abstract = {Conventional oil and gas fields are increasingly difficult to explore and produce, which calls for more complex wave-equation based inversion (WEI) algorithms that require dense long-offset samplings. These requirements result in an exponential growth in data volumes and prohibitive demands on computational resources. In this work, we propose a fast, resilient, and scalable wave-equation based inversion methodology for both the data and image-domain, which can handle complicated wave physics. First we show that both the data and image domains exhibit low-rank structure in a transform-domain, which can be exploited to compress the dense data or image volumes. Then, by accessing information from the compressed volumes on-the fly, we devise a scalable computational inversion framework driven by gradient calculations, which works with small subsets of source experiments. In the full-waveform inversion context, we demonstrate the efficacy of low-rank interpolation to improve downstream inversion results compared to merely inverting the velocity model using the subsampled data volume directly. We demonstrate the effectiveness of the proposed framework on two different waveform-inversion formulations, specifically performing full-waveform inversion on 5D data set generated using the overthrust model, and wave-equation based migration velocity analysis on the Marmousi model. This is joint work with Curt Da Silva and Yiming Zhang.},
	keywords = {Presentation, SINBAD, SINBADFALL2016, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/kumar2016SINBADFlrm/kumar2016SINBADFlrm.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2016/Fall/kumar2016SINBADFlrm/kumar2016SINBADFlrm.mov},
	author = {Rajiv Kumar and Curt Da Silva and Yiming Zhang and Felix J. Herrmann}
}
@presentation {wason2015SINBADldr,
	title = {Latest developments in randomized marine acquisition and source separation},
	year = {2015},
	publisher = {SINBAD},
	abstract = {We present the latest developments in randomized (simultaneous) marine acquisition for time-lapse surveys. We also show a comparison between two source separation algorithms {\textemdash} sparsity-promotion and rank-minimization {\textemdash} for simultaneous data.},
	keywords = {Presentation, SINBAD, SINBADSPRING2015, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2015/Spring/wason2015SINBADldr/wason2015SINBADldr.pdf},
	author = {Haneet Wason and Rajiv Kumar and Felix Oghenekohwo and Felix J. Herrmann}
}
@conference {kumar2015CSEGlse,
	title = {Least-squares extended imaging with surface-related multiples},
	booktitle = {CSEG Annual Conference Proceedings},
	year = {2015},
	note = {(CSEG, Calgary)},
	month = {05},
	abstract = {Common image gathers are used in building velocity models, inverting anisotropy parameters, and analyzing reservoir attributes. Often primary reflections are used to form image gathers and multiples are typically attenuated in processing to remove strong coherent artifacts generated by multiples that interfere with the imaged reflectors. However, researchers have shown that, if correctly used, multiples can actually provide extra illumination of the subsurface in seismic imaging, especially for delineating the near-surface features. In this work, we borrow ideas from literatures on imaging with surface-related multiples, and apply these ideas to extended imaging. This way we save the massive computation cost in separating multiples from the data before using them during the formation of image gathers. Also, we mitigate the strong coherent artifacts generated by multiples which can send the migration velocity analysis type algorithms in wrong direction. Synthetic examples on a three-layer model show the efficacy of the proposed formulation.},
	keywords = {CSEG, image-gather, surface-related multiples},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2015/kumar2015CSEGlse/kumar2015CSEGlse.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/CSEG/2015/kumar2015CSEGlse/kumar2015CSEGlse_pres.pdf},
	author = {Rajiv Kumar and Ning Tu and Tristan van Leeuwen and Felix J. Herrmann}
}
@techreport {kumar2015EAGElse,
	title = {Least-squares extended imaging with surface-related multiples},
	number = {TR-EOAS-2015-1},
	year = {2015},
	month = {01},
	institution = {UBC},
	abstract = {Common image gathers are used in building velocity models, inverting for anisotropy parameters, and analyzing reservoir attributes. Typically, only primary reflections are used to form image gathers as multiples can cause artifacts that interfere with the events of interest. However, it has been shown that multiples can actually provide extra illumination of the subsurface, especially for delineating the near-surface features. In this paper, we aim to form common image gathers directly from the data with surface related multiples by applying concepts that have been used to successfully deal with surface-related multiples in imaging. We achieve this by effectively inverting an extended migration operator. This results in extended images with better near-surface illumination that are free of artifacts that can hamper velocity analysis. In addition, being able to generate extended images directly from the total data avoids the need for (time-consuming) pre-processing. Synthetic examples on a layered model show that the proposed formulation is promising.},
	keywords = {image gathers, surface-related multiples},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2015/kumar2015EAGElse/kumar2015EAGElse.html},
	author = {Rajiv Kumar and Ning Tu and Tristan van Leeuwen and Felix J. Herrmann}
}
@presentation {li2013SINBADgom,
	title = {Lessons learned from Chevron Gulf of Mexico data set},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Chevron Gulf of Mexico data set is very challenging for FWI because of elastic phases, limited offset, lack of low frequencies and salt structure. To overcome these issue, we first use ray-based tomography on the hand-picked first breaks to generate initial model for FWI, and then we apply curvelet-denosing techniques to improve the poor signal-to-noise ratio of the observed data at low frequencies. Finally, Curvelet domain sparsity promoting Gauss-Newton FWI helps to suppress model space artifacts caused by elastic phases. This is joint work with Andrew J. Calvert, Ian Hanlon, Mostafa Javanmehri, Rajiv Kumar, Tristan van Leeuwen, Brendan R. Smithyman, Haneet Wason and Felix J. Herrmann},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/li2013SINBADgom/li2013SINBADgom.pdf},
	author = {Xiang Li and Andrew J. Calvert and Ian Hanlon and Mostafa Javanmehri and Rajiv Kumar and Tristan van Leeuwen and Brendan R. Smithyman and Haneet Wason and Felix J. Herrmann}
}
