Source estimation is essential for all wave-equation-based seismic inversions, including full-waveform inversion (FWI) and the recently proposed wavefield-reconstruction inversion (WRI). When the source estimation is inaccurate, errors will propagate into the predicted data and introduce additional data misfit. As a consequence, inversion results that minimize this data misfit may become erroneous. To mitigate the errors introduced by the incorrect and preestimated sources, an embedded procedure that updates sources along with medium parameters is necessary for the inversion. So far, such a procedure is still missing in the context of WRI, a method that is, in many situations, less prone to local minima related to so-called cycle skipping, compared with FWI through exact data fitting. Although WRI indeed helps to mitigate issues related to cycle skipping by extending the search space with wavefields as auxiliary variables, it relies on having access to the correct source functions. To remove the requirement of having the accurate source functions, we have developed a source-estimation technique specifically designed for WRI. To achieve this task, we consider the source functions as unknown variables and arrive at an objective function that depends on the medium parameters, wavefields, and source functions. During each iteration, we apply the so-called variable projection method to simultaneously project out the source functions and wavefields. After the projection, we obtain a reduced objective function that only depends on the medium parameters and invert for the unknown medium parameters by minimizing this reduced objective. Numerical experiments illustrate that this approach can produce accurate estimates of the unknown medium parameters without any prior information of the source functions.
Acoustoelastic full-waveform inversion for transcranial ultrasound computed tomography
