@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}
}
@mastersthesis {tu2015THfis,
	title = {Fast imaging with surface-related multiples},
	year = {2015},
	note = {(PhD)},
	month = {08},
	school = {The University of British Columbia},
	type = {phd},
	address = {Vancouver},
	abstract = {Surface-related multiples, which are waves that bounce more than once between the water surface and the subsurface reflectors, constitute a significant part of the data acquired in marine seismic surveys. If left untreated, they can lead to misplaced phantom reflectors in the image, and result in erroneous interpretations of the subsurface structure. As a result, these multiples are removed before the imaging procedure in conventional seismic data processing. However, because they interact more with the subsurface medium, they may carry extra information that is not present in the primaries. Therefore instead of removing these multiples, a more desirable alternative is to make active use of them. We derive from the well-established "Surface-Related Multiple Elimination" relation, and arrive at a linearized expression of the wave-equation based modelling that incorporates the surface- related multiples. We then present a computationally efficient approach to iteratively invert this expression to obtain an image of the subsurface from data that contain multiples. We achieve the computational efficiency inside each iteration by (i) using the wave-equation solver to implicitly carry out the expensive multiple prediction; and (ii) reducing the number of wave-equation solves during data simulation by subsampling the monochromatic source experiments. We show that, compared with directly applying the cross-correlation/deconvolutional imaging conditions, the presented approach can suppress the coherent imaging artifacts from multiples more effectively. We also show that, by curvelet-domain sparsity promoting and occasionally drawing new data samples during the inversion, the proposed inversion method gains improved robustness to velocity errors in the background model, as well as modelling errors incurred during linearization of the wave-equation. To combine the information encoded in both the primaries and the multiples, we then propose a highly accurate source estimation method to jointly invert the total upgoing wavefield. We show with field data examples that we can reap benefits from both the relative noise-free primaries and the extra illumination coverage of the multiples. We also demonstrate that the inclusion of multiples help mitigate the amplitude ambiguity during source estimation. We conclude the thesis with an outlook for future research directions, as well as potential extensions of the proposed work.},
	keywords = {inversion, least-squares, multiples, seismic imaging},
	url = {https://slim.gatech.edu/Publications/Public/Thesis/2015/tu2015THfis/tu2015THfis.pdf},
	presentation = {https://slim.gatech.edu/Publications/Public/Thesis/2015/tu2015THfis/tu2015THfis_pres.pdf},
	author = {Ning Tu}
}
@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}
}
