@article {tu2015GJIsem,
	title = {Source estimation with surface-related multiples{\textendash}-fast ambiguity-resolved seismic imaging},
	journal = {Geophysical Journal International},
	volume = {205},
	year = {2016},
	note = {(published online in Geophysical Journal International)},
	month = {03},
	pages = {1492-1511},
	abstract = {We address the problem of obtaining a reliable seismic image without prior knowledge of the source wavelet, especially from data that contain strong surface-related multiples. Conventional reverse-time migration requires prior knowledge of the source wavelet, which is either technically or computationally challenging to accurately determine; inaccurate estimates of the source wavelet can result in seriously degraded reverse-time migrated images, and therefore wrong geological interpretations. To solve this problem, we present a "wavelet-free" imaging procedure that simultaneously inverts for the source wavelet and the seismic image, by tightly integrating source estimation into a fast least-squares imaging framework, namely compressive imaging, given a reasonably accurate background velocity model. However, this joint inversion problem is difficult to solve as it is plagued with local minima and the ambiguity with respect to amplitude scalings because of the multiplicative, and therefore nonlinear, appearance of the source wavelet in the otherwise linear formalism. We have found a way to solve this nonlinear joint-inversion problem using a technique called variable projection, and a way to overcome the scaling ambiguity by including surface-related multiples in our imaging procedure following recent developments in surface-related multiple prediction by sparse inversion. As a result, we obtain without prior knowledge of the source wavelet high-resolution seismic images, comparable in quality to images obtained assuming the true source wavelet is known. By leveraging the computationally efficient compressive-imaging methodology, these results are obtained at affordable computational costs compared with conventional processing work flows that include surface-related multiple removal and reverse-time migration.},
	keywords = {computational seismology, free surface, Inverse theory, multiples, seismic imaging, time series analysis, wave propagation},
	url = {https://slim.gatech.edu/Publications/Public/Journals/GeophysicalJournalInternational/2016/tu2015GJIsem/tu2015GJIsem.pdf},
	url2 = {http://gji.oxfordjournals.org/content/205/3/1492.full.pdf?keytype=ref\&ijkey=zFydRGQpiJzeCS9},
	author = {Ning Tu and Aleksandr Y. Aravkin and Tristan van Leeuwen and Tim T.Y. Lin and Felix J. Herrmann}
}
@conference {tu2015PIMSsls,
	title = {Sparse least-squares seismic imaging with source estimation utilizing multiples},
	booktitle = {PIMS Workshop on Advances in Seismic Imaging and Inversion},
	year = {2015},
	note = {(PIMS Workshop on Advances in Seismic Imaging and Inversion, University of Alberta, Edmonton)},
	month = {05},
	abstract = {We present an least-squares seismic imaging method that (1) is efficient by subsampling the monochromatic source experiments; (2) does not require the prior knowledge of the source wavelet, by inverting for the seismic image and the source wavelet simultaneously using variable projection; and (3) makes active use of surface-related multiples to increase the illumination and to resolve the scaling ambiguity in source estimation, by tightly integrating the well-established surface-related-multiple-prediction into the wave-equation based modelling. As a result, we are able to obtain ambiguity-resolved, high-fidelity least-squares inverted seismic images with computational costs that are comparable to conventional reverse-time migrations with all the data. We demonstrate the efficacy of the proposed method using synthetic examples.},
	keywords = {inversion, least-squares migration, multiples, PIMS, seismic imaging, source estimation, workshop},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/PIMS/2015/tu2015PIMSsls/tu2015PIMSsls_pres.pdf},
	author = {Ning Tu and Felix J. Herrmann}
}
@conference {witte2015IIPFWIspl,
	title = {Sparsity-promoting least-square migration with linearized Bregman and compressive sensing},
	booktitle = {Inaugural Full-Waveform Inversion Workshop},
	year = {2015},
	note = {(Natal, Brazil)},
	month = {08-09},
	abstract = {We present a novel adaptation of a recently developed relatively simple iterative algorithm to solve large-scale sparsity-promoting optimization problems. Our algorithm is particularly suitable to large-scale geophysical inversion problems, such as sparse least-squares reverse-time migration or Kirchoff migration since it allows for a tradeoff between parallel computations, memory allocation, and turnaround times, by working on subsets of the data with different sizes. Comparison of the proposed method for sparse least-squares imaging shows a performance that rivals and even exceeds the performance of state-of-the art one-norm solvers that are able to carry out least-squares migration at the cost of a single migration with all data.},
	keywords = {Bregman, compressed sensing, FWI, migration, sparsity promotion, time domain},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/IIPFWI/witte2015IIPFWIspl/witte2015IIPFWIspl_pres.pdf},
	author = {Philipp A. Witte and Ning Tu and Ernie Esser and Mengmeng Yang and Mathias Louboutin and Felix J. Herrmann}
}
@presentation {tu2014SINBADsis,
	title = {Sparse inversion simplified},
	year = {2014},
	publisher = {SINBAD},
	abstract = {Implementation of the SPGl1 solver in low-level computer languages such as C or Fortran may sound daunting to many. This talk is tailored for you if you would like to reap benefits from using sparse constraints to regularize your inversion problem, and at the same time want some easy-to-implement algorithm to solve it. We discuss in this talk the linearized Bregman projection method that embraces both ingredients, and show promising applications in fast least-squares imaging and seismic data interpolation problems.},
	keywords = {Presentation, SINBAD, SINBADFALL2014, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2014/Fall/tu2014SINBADsis/tu2014SINBADsis.pdf},
	author = {Ning Tu and Navid Ghadermarzy and Rajiv Kumar and Felix J. Herrmann}
}
@presentation {tu2013SINBADmachar,
	title = {SLIM{\textquoteright}s findings on the Machar dataset},
	year = {2013},
	publisher = {SINBAD},
	abstract = {By courtesy of BP, we are able to get some hands-on experience imaging seismic data from the Machar field in North Sea. We performed reverse-time migration and sparsity-promoting migration with source-estimation to this dataset (or a cropped section of the dataset in many cases due to computational constraints), and had some interesting findings. In this presentation, we will show some conclusive results we have got, and explain the key techniques in imaging this dataset.},
	keywords = {Presentation, SINBAD, SINBADFALL2013, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2013/Fall/tu2013SINBADmachar/tu2013SINBADmachar.pdf},
	author = {Ning Tu and Felix J. Herrmann}
}
@conference {aravkin2013ICASSPssi,
	title = {Sparse seismic imaging using variable projection},
	booktitle = {ICASSP},
	year = {2013},
	month = {05},
	abstract = {We consider an important class of signal processing problems where the signal of interest is known to be sparse, and can be recovered from data given auxiliary information about how this data was generated. For example, a sparse green{\textquoteright}s function may be recovered from seismic experimental data using sparsity optimization when the source signature is known. Unfortunately, in practice this information is often missing, and must be recovered from data along with the signal using deconvolution techniques. In this paper, we present a novel methodology to simulta- neously solve for the sparse signal and auxiliary parameters using a recently proposed variable projection technique. Our main contribution is to combine variable projection with spar- sity promoting optimization, obtaining an efficient algorithm for large-scale sparse deconvolution problems. We demon- strate the algorithm on a seismic imaging example.},
	keywords = {Imaging, Optimization, sparsity, variable projection},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/ICASSP/2013/aravkin2013ICASSPssi/aravkin2013ICASSPssi.pdf},
	author = {Aleksandr Y. Aravkin and Tristan van Leeuwen and Ning Tu}
}
@conference {tu2011EAGEspmsrm,
	title = {Sparsity-promoting migration with surface-related multiples},
	booktitle = {EAGE Annual Conference Proceedings},
	year = {2011},
	month = {05},
	abstract = {Multiples, especially the surface-related multiples, form a significant part of the total up-going wave- field. If not properly dealt with, they can lead to false reflectors in the final image. So conventionally practitioners remove them prior to migration. Recently research has revealed that multiples can actually provide extra illumination so different methods are proposed to address the issue that how to use multiples in seismic imaging, but with various kinds of limitations. In this abstract, we combine primary estimation and sparsity-promoting migration into one convex-optimization process to include information from multiples. Synthetic examples show that multiples do make active contributions to seismic migration. Also by this combination, we can benefit from better recoveries of the Greens function by using sparsity-promoting algorithms since reflectivity is sparser than the Greens function.},
	keywords = {EAGE, Imaging, Presentation, Processing},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/tu11EAGEspmsrm/tu11EAGEspmsrm.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/EAGE/2011/tu11EAGEspmsrm/tu11EAGEspmsrm_pres.pdf},
	url2 = {http://earthdoc.eage.org/publication/publicationdetails/?publication=50369},
	author = {Ning Tu and Tim T.Y. Lin and Felix J. Herrmann}
}
@presentation {tu2010SINBADspm,
	title = {Sparsity promoting migration with surface-related multiples},
	year = {2010},
	publisher = {SINBAD},
	keywords = {Presentation, SINBAD, SINBADFALL2010, SLIM},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SINBAD/2010/Fall/tu2010SINBADspm/tu2010SINBADspm_pres.pdf},
	author = {Ning Tu}
}
@conference {lin2010SEGspm,
	title = {Sparsity-promoting migration from surface-related multiples},
	booktitle = {SEG Technical Program Expanded Abstracts},
	volume = {29},
	year = {2010},
	month = {10},
	pages = {3333-3337},
	publisher = {SEG},
	organization = {SEG},
	abstract = {Seismic imaging typically begins with the removal of multiple energy in the data, out of fear that it may introduce erroneous structure. However, seismic multiples have effectively seen more of the earth{\textquoteright}s structure, and if treated correctly can potential supply more information to a seismic image compared to primaries. Past approaches to accomplish this leave ample room for improvement; they either require extensive modification to standard migration techniques, rely too much on prior information, require extensive pre-processing, or resort to full-waveform inversion. We take some valuable lessons from these efforts and present a new approach balanced in terms of ease of implementation, robustness, efficiency and well-posedness, involving a sparsity-promoting inversion procedure using standard Born migration and a data-driven multiple modeling approach based on the focal transform.},
	keywords = {Presentation, Processing, SEG},
	doi = {10.1190/1.3513540},
	url = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2010/lin10SEGspm/lin10SEGspm.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Conferences/SEG/2010/lin10SEGspm/lin10SEGspm_pres.pdf},
	author = {Tim T.Y. Lin and Ning Tu and Felix J. Herrmann}
}
