@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}
}
@mastersthesis {yang2020THsiw,
	title = {Seismic imaging with extended image volumes and source estimation},
	year = {2020},
	note = {(PhD)},
	month = {03},
	school = {Georgia Institute of Technology},
	type = {phd},
	address = {Atlanta},
	abstract = {Seismic imaging is an important tool for the exploration and production of oil \& gas, carbon sequestration, and the mitigation of geohazards. Through the process of seismic migration, images of subsurface geological structures are created from data collected at the surface. These images reflect changes in the physical rock properties such as wave speed and density. While significant progress has been made in the development of 3D imaging technology for complex geological areas, several challenges remain, some of which are addressed in this thesis. The first main contribution of this thesis is in the area of creating so-called subsurface-offset gathers, which play an increasingly important role in seismic imaging because they provide a multitude of information ranging from the reflection mechanism itself to information of the dips of specific reflectors and the accuracy of the background velocity model. Unfortunately, the formation and manipulation of these gathers come with exceedingly high computational and storage costs because extended image volumes are quadratic in the image size. These high costs are avoided by using techniques from modern randomized linear algebra that allow for compression of extended image volumes into low-rank factorized form{\textemdash}i.e., the image volume is approximately written as an outer product of a tall and a wide matrix. It is demonstrated that this factorization provides access to different types of sub-surface offset gathers, including common-image (point) gathers, without the need to explicitly form this outer product. As a result, challenging steep dip imaging situations, where conventional horizontal offset gathers no longer focus, can be handled. Moreover, extended image volumes for one background velocity model can directly be mapped to those of another background velocity model. As a result, factorization costs are incurred only once when examining imaging scenarios for different background velocity models. The second main contribution of this thesis is on the development of computationally efficient sparsity-promoting imaging techniques and on-the-fly source estimation. In this work, an adaptive technique is proposed where the unknown time signature of the sources is estimated during imaging. Without accurate knowledge of these source signatures, seismic images can be wrongly positioned and can have the wrong amplitudes hampering subsequent geophysical and geological interpretations. With the presented technique, this problem is mitigated. Finally, a contribution is made to address the detrimental effects of surface-related multiples. If not handled correctly, these multiples give rise to unwanted artifacts in the image. A new technique is introduced to address this issue in realistic settings where there is a strong density contrast at the ocean bottom. As a result, the surface-related multiples are mapped to the reflectors. Because bounce points at the surface can be considered as sources, this mapping of the multiples rather than removal increases the subsurface illumination.},
	keywords = {extended image volumes, low rank, multiples, PhD, randomized linear algebra, source estimation, sparsity-promoting inversion},
	url = {https://slim.gatech.edu/Publications/Public/Thesis/2020/yang2020THsiw/yang2020THsiw.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Thesis/2020/yang2020THsiw/yang2020THsiw_pres.pdf},
	author = {Mengmeng Yang}
}
@article {silva2018alr,
	title = {Applications of low-rank compressed seismic data to full waveform inversion and extended image volumes},
	journal = {Geophysics},
	year = {2019},
	note = {Accepted on January 22, 2019.},
	abstract = {Conventional oil and gas fields are increasingly difficult to explore and image, resulting in the call for more complex wave-equation based inversion algorithms that require dense long-offset samplings. Consequently, there is an exponential growth in the size of data volumes and prohibitive demands on computational resources. In this work, we propose a method to compress and process seismic data directly in a low-rank tensor format, which drastically reduces the amount of storage required to represent the data. We first outline how seismic data exhibits low-rank structure in a particular transform-domain, which can be exploited to compress the dense data in one extremely storage-efficient tensor format when the data is fully sampled. In the more realistic case of missing data, we can use efficient interpolation techniques to approximate the fully sampled volume in compressed form. In either case, once we have our data represented in its compressed tensor form, we design an algorithm to extract source or receiver gathers directly from the compressed parameters. This extraction process can be done on-the-fly directly on the compressed data, in a full waveform inversion context, and does not require scanning through the entire dataset in order to form shot gathers. To the best of our knowledge, this work is one of the first major contributions to working with seismic data applications directly in the compressed domain without reconstructing the entire data volume. We use a stochastic inversion approach, which works with small subsets of source experiments at each iteration, further reducing the computational and memory costs of full waveform inversion. We also demonstrate how this data compression and extraction technique can be applied to forming full subsurface image gathers through probing techniques.},
	keywords = {common shot/receiver gather, extended image volume, full waveform inversion, low rank, tensor},
	doi = {10.1190/geo2018-0116.1},
	url = {https://slim.gatech.edu/Publications/Public/Journals/Geophysics/2019/silva2018alr/silva2018alr.pdf},
	author = {Curt Da Silva and Yiming Zhang and Rajiv Kumar 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}
}
@techreport {esser2014sln,
	title = {Some lifting notes},
	number = {TR-EOAS-2014-1},
	year = {2014},
	note = {written on February 15, 2014},
	month = {02},
	institution = {UBC},
	keywords = {convex semidefinite programming, lifting, low rank, nonconvex quadratic problems},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2014/esser2014sln/esser2014sln.pdf},
	author = {Ernie Esser}
}
@techreport {kumar2013ICMLlr,
	title = {An SVD-free Pareto curve approach to rank minimization},
	number = {TR-EOAS-2013-2},
	year = {2013},
	month = {02},
	institution = {UBC},
	abstract = {Recent SVD-free matrix factorization formulations have enabled rank optimization for extremely large-scale systems (millions of rows and columns). In this paper, we consider rank-regularized formulations that only require a target data-fitting error level, and propose an algorithm for the corresponding problem. We illustrate the advantages of the new approach using the Netflix problem, and use it to obtain high quality results for seismic trace interpolation, a key application in exploration geophysics. We show that factor rank can be easily adjusted as the inversion proceeds, and propose a weighted extension that allows known subspace information to improve the results of matrix completion formulations. Using these methods, we obtain high-quality reconstructions for large scale seismic interpolation problems with real data.},
	keywords = {interpolation, low rank},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2013/kumar2013ICMLlr/kumar2013ICMLlr.pdf},
	author = {Aleksandr Y. Aravkin and Rajiv Kumar and Hassan Mansour and Ben Recht and Felix J. Herrmann}
}
