
An Analytic Solution to a Coupled System of Equations for Modeling Photoacoustic Trace Gas Sensors and a Full Waveform Inversion Approach to Microseismic Source Estimation 
Author(s): 
Kaderli, Jordan 

Minkoff, Susan E. 

201805 

Dissertation 
Keywords:  Show Keywords 

We discuss two different problems in applied mathematics. The first is related to the modeling of a certain class of trace gas sensors. The second is the geophysical inverse problem of estimating the source of a microseismic event. Quartzenhanced photoacoustic spectroscopy and resonant optothermoacoustic detection are two promising techniques used in trace gas sensing. Both methods use a quartz tuning fork and modulated laser source to detect trace gases. We discuss a coupled system of equations for the pressure, temperature, and velocity of a fluid that accounts for both thermal effects and viscous damping. We derive an analytical solution to a pressuretemperature subsystem of the MorseIngard equations in the special case of cylindrical symmetry. We solve for the pressure and temperature in an infinitely long cylindrical fluid domain with a source function given by a constantwidth Gaussian beam that is aligned with the axis of the cylinder. In addition, we surround this cylinder with an infinitely long annular solid domain, and we couple the pressure and temperature in the fluid domain to the temperature in the solid. We show that the temperature in the solid near the fluidsolid interface can be an order of magnitude larger than that computed using a simpler model in which the temperature is governed by the heat equation. We also verify that the temperature solution of the coupled system exhibits a thermal boundary layer. These results suggest that for computational modeling of resonant optothermoacoustic detection sensors, the temperature in the fluid should be computed by solving the MorseIngard equations rather than the heat equation. In the second problem, we use full waveform inversion to estimate the full spatial and temporal description of a microseismic source which includes not only the location and origin time of the source but also the waveform itself. Beginning with the simplifying assumption of twodimensional acoustic wave propagation, we compute the gradient via the adjointstate method for both the spatial radiation pattern and the temporal waveform of the source. This approach identifies multiple sources, handles extremely low signaltonoise ratio data, and produces accurate results in the absence of a good initial estimate. Encouraged by the promising results of the two dimensional acoustic case, we apply the approach to the case where wave propagation is modeled via the velocitystress formulation of the elastic wave equation in three dimensions. A change of variables applied to the velocitystress formulation ensures that the system is selfadjoint. Thus full waveform inversion can be effectively tailored to use this transformed velocitystress system to estimate microseismic events with limited modifications to the forward wave solver. The inversion produces either a spatial source description in the form of a volume in which each point in the domain indicates the magnitude of stress perturbation and therefore deformation of the medium or a temporal source description which gives the time evolution of the deformation. The inversion does not require any a priori assumptions about the form of the source and does not require a good starting guess for accurate source recovery. 

PHD 

Doctoral 

http://hdl.handle.net/10735.1/5850 

Copyright ©2018 is held by the author. Digital access to this material is made possible by the Eugene McDermott Library. Further transmission, reproduction or presentation (such as public display or performance) of protected items is prohibited except with permission of the author. 

text 

Mathematics 
Files  Size  Format  View  

ETD5608022KADERLI7864.76.pdf  3.844Mb 
View/ 
Dissertation 