@conference {zhang2017SEGmsd,
	title = {Massive 3D seismic data compression and inversion with hierarchical Tucker},
	booktitle = {SEG Technical Program Expanded Abstracts},
	year = {2017},
	note = {(SEG, Houston)},
	month = {09},
	pages = {1347-1352},
	abstract = {Modern-day oil and gas exploration, especially in areas of complex geology such as fault belts and sub-salt areas, is an increasingly expensive and risky endeavour. Typically long-offset and dense sampling seismic data are required for subsequent shot based processing procedures, e.g. wave-equation based inversion (WEI) and surface-related multiple elimination (SRME). However, these strict requirements result in an exponential growth in data volume size and prohibitive demands on computational resources, given the multidimensional nature of the data volumes. Moreover the physical constraints and cost limitations impose restrictions on acquiring fully sampled data. In this work, we propose to invert our large-scale data from a set of subsampled measurements, resulting in an estimate of the true volume in a compressed low-rank tensor format. Rather than expanding the data to its fully-sampled form for later downstream processes, we demonstrate how to use this compressed data directly via on-the-fly common shot or receiver gathers extraction. The combination of massive compression and fast on demand data reconstruction of 3D shot or receiver gathers leads to a substantial reduction in memory costs but with minimal effects on results in the subsequent processing procedures. We demonstrate the effective implementation of our proposed framework on full-waveform inversion on a 3D seismic synthetic data set generated from a Overthrust model.},
	keywords = {3D, FWI, SEG, tensor algebra},
	doi = {10.1190/segam2017-17742951.1},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/zhang2017SEGmsd/zhang2017SEGmsd.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2017/zhang2017SEGmsd/zhang2017SEGmsd_pres.pdf},
	author = {Yiming Zhang and Curt Da Silva and Rajiv Kumar and Felix J. Herrmann}
}
@presentation {zhang2017SINBADFmsd,
	title = {Massive seismic data compression \& recovery w/ on-the-fly data extraction},
	year = {2017},
	publisher = {SINBAD},
	abstract = {Industrial seismic exploration has moved towards complex geological areas, which requires typically long-offset and dense sampling data in order to avoid aliasing and inaccuracy in wave-equation based inversion algorithms. These strict requirements lead to massive data volume size and prohibitive demands on computational resources. In this work, we propose to compress our dense data in hierarchical Tucker tensor format by exploiting the low-rank structure of the data in a transformed domain. Then, we devise on-the-fly common shot or receiver gather extraction directly via the highly compressed factors. In subsampling scenarios, by interpolating this novel tensor format, we can also reconstruct the shot or receiver gather on a per-query basis rather than expanding the data to its fully-sampled form. We demonstrate the effective performance of our proposed technique on 3D stochastic full-waveform inversion, which allows the stochastic algorithm to extract shot gathers as it requires them throughout the inversion process. Moreover, we finally show how to computational effectively generate the CIGs from this compressed low-rank tensor representation of the data with the help of fast simultaneous shot or receiver gather generation.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/zhang2017SINBADFmsd/zhang2017SINBADFmsd.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/zhang2017SINBADFmsd/zhang2017SINBADFmsd.mov},
	author = {Yiming Zhang and Curt Da Silva and Rajiv Kumar and Felix J. Herrmann}
}
@presentation {lopez2017SINBADFmci,
	title = {Matrix Completion in Parallel Architectures: Julia Implementation},
	year = {2017},
	publisher = {SINBAD},
	abstract = {Matrix completion techniques offer potential tools for frugal seismic data acquisition, where dense acquisition is replaced by optimization. This shift of focus means that efficient numerical methods are critical for the implementation of these techniques in large-scale seismic applications. To this end, this talk modifies rank-penalization methodologies to suit parallel architectures. By adopting factorization-based alternating minimization schemes, each program can be decoupled into independent sub-problems handled in parallel. We showcase a distributed parallel execution in Julia and explore the scalability of the approach.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/lopez2017SINBADFmci/lopez2017SINBADFmci.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/lopez2017SINBADFmci/lopez2017SINBADFmci.mov},
	author = {Oscar Lopez and Keegan Lensink and Rajiv Kumar and Henryk Modzelewski}
}
@article {yarman2017GPmbd,
	title = {A model based data driven dictionary learning for seismic data representation},
	journal = {Geophysical Prospecting},
	year = {2017},
	note = {(published online in Geophysical Prospecting)},
	month = {04},
	abstract = {Planar waves events recorded in a seismic array can be represented as lines in the Fourier domain. However, in the real-world seismic events usually have curvature or amplitude variability that means their Fourier transforms are no longer strictly linear, but rather occupy conic regions of the Fourier domain that are narrow at low frequencies, but broaden at high frequencies where the effect of curvature becomes more pronounced. One can consider these regions as localised "signal cones". In this work, we consider a space-time variable signal cone to model the seismic data. The variability of the signal cone is obtained through scaling, slanting and translation of the kernel for cone-limited (C-limited) functions (functions whose Fourier transform lives within a cone) or C-Gaussian function (a multivariate function whose Fourier transform decays exponentially with respect to slowness and frequency) which constitutes our dictionary. We find a discrete number of scaling, slanting and translation parameters from a continuum by optimally matching the data. This is a nonlinear optimization problem which we address by a fixed point method which utilizes a variable projection method with e1 constraints on the linear parameters and bound constraints on the nonlinear parameters. We observe that slow decay and oscillatory behavior of the kernel for C-limited functions constitute bottlenecks for the optimization problem which we partially overcome by the C-Gaussian function. We demonstrate our method through an interpolation example. We present the interpolation result using the estimated parameters obtained from the proposed method and compare it with those obtained using sparsity promoting curvelet decomposition, matching pursuit Fourier interpolation and sparsity promoting plane wave decomposition methods.},
	keywords = {dictionary learning, kernel method, nonlinear optimization, reproducing kernel, variable projection, variational method},
	doi = {10.1111/1365-2478.12533},
	url = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalProspecting/2017/yarman2017GPmbd/yarman2017GPmbd.pdf},
	author = {Can Evren Yarman and Rajiv Kumar and James Rickett}
}
@presentation {kumar2017SINBADFmdt,
	title = {Multi-domain target-oriented imaging using extreme-scale matrix factorization},
	year = {2017},
	publisher = {SINBAD},
	abstract = {In this work, we present an alternative approach to redatum both source and receivers at depth, under the framework of reflectivity-based extended images with two-way wave propagation in the background medium. We propose a randomized svd based probing scheme that takes advantage of the algebraic structure of the extended imaging system to overcome the computational cost and memory usage associated with the number of wave-equation solutions and explicit storage employed by conventional migration methods. Experimental results on complex geological models demonstrate the efficacy of proposed methodology in performing multi-domain target imaging.},
	keywords = {Presentation, SINBAD, SINBADFALL2017, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/kumar2017SINBADFmdt/kumar2017SINBADFmdt.pdf},
	url2 = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2017/Fall/kumar2017SINBADFmdt/kumar2017SINBADFmdt.mov},
	author = {Rajiv Kumar and Marie Graff-Kray and Ivan Vasconcelos and Felix J. Herrmann}
}
@conference {kumar2015EAGEmcu,
	title = {Matrix completion on unstructured grids : 2-D seismic data regularization and interpolation},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2015},
	note = {(EAGE, Madrid)},
	month = {06},
	abstract = {Seismic data interpolation via rank-minimization techniques has been recently introduced in the seismic community. All the existing rank-minimization techniques assume the recording locations to be on a regular grid, e.g., sampled periodically, but seismic data are typically irregularly sampled along spatial axes. Other than the irregularity of the sampled grid, we often have missing data. In this paper, we study the effect of grid irregularity to conduct matrix completion on a regular grid for unstructured data. We propose an improvement of existing rank-minimization techniques to do regularization. We also demonstrate that we can perform seismic data regularization and interpolation simultaneously. We illustrate the advantages of the modification using a real seismic line from the Gulf of Suez to obtain high quality results for regularization and interpolation, a key application in exploration geophysics.},
	keywords = {EAGE, interpolation, matrix completion, NFFT, regularization},
	doi = {10.3997/2214-4609.201413448},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/kumar2015EAGEmcu/kumar2015EAGEmcu.html},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2015/kumar2015EAGEmcu/kumar2015EAGEmcu_poster.pdf},
	author = {Rajiv Kumar and Oscar Lopez and Ernie Esser and Felix J. Herrmann}
}
@techreport {kumar2014SEGmcu,
	title = {Matrix completion on unstructured grids : 2-D seismic data regularization and interpolation},
	number = {TR-EOAS-2014-3},
	year = {2014},
	month = {04},
	institution = {UBC},
	abstract = {Seismic data interpolation via rank-minimization techniques has been recently introduced in the seismic community. All the existing rank-minimization techniques assume the recording locations to be on a regular grid, e.g. sampled periodically, but seismic data are typically irregularly sampled along spatial axes. Other than the irregularity of the sampled grid, we often have missing data. In this paper, we study the effect of grid irregularity to conduct matrix completion on a regular grid for unstructured data. We propose an improvement of existing rank-minimization techniques to do regularization. We also demonstrate that we can perform seismic data regularization and interpolation simultaneously. We illustrate the advantages of the modification using a real seismic line from the Gulf of Suez to obtain high quality results for regularization and interpolation, a key application in exploration geophysics.},
	keywords = {interpolation, matrix completion, NFFT, regularization},
	url = {https://slim.gatech.edu/Publications/Public/TechReport/2014/kumar2014SEGmcu/kumar2014SEGmcu.html},
	author = {Rajiv Kumar and Oscar Lopez and Ernie Esser and Felix J. Herrmann}
}
@presentation {akalin2013SINBADmtc,
	title = {Matrix and tensor completion for large-scale seismic interpolation: a comparative study},
	year = {2013},
	publisher = {SINBAD},
	abstract = {Owing to their high dimensionality, interpolating 3D seismic data volumes remains a computationally daunting task. In this work, we outline a comprehensive framework for sampling and interpolating such volumes based on the well-understood theory of Matrix and Tensor completion. This interpolation theory consists of three components major signal structure, structure-destroying sampling, and structure-restoring optimization. By viewing interpolation in the context of this theory, we are able to specify exactly when these approaches are expected to perform well. We also introduce structure-revealing transformations that promote the inherent low-rank structure in seismic data as well as a factorization approach that scales to large problem sizes. Our methods are able to handle large-scale data volumes more accurately and more quickly compared to other more ad-hoc approaches, as we will demonstrate. This is joint work with Curt Da Silva, Rajiv Kumar, Ben Recht, and Felix J. Herrmann.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/akalin2013SINBADmtc/akalin2013SINBADmtc.pdf},
	author = {Okan Akalin and Curt Da Silva and Rajiv Kumar and Ben Recht and Felix J. Herrmann}
}
