Although transcranial photoacoustic imaging has been previously investigated by several groups, there are many unknowns about the distorting effects of the skull due to the impedance mismatch between the skull and underlying layers. The current computational methods based on finite-element modeling are slow, especially in the cases where fine grids are defined for a large 3-D volume. We develop a very fast modeling/simulation framework based on deterministic ray-tracing. The framework considers a multilayer model of the medium, taking into account the frequency-dependent attenuation and dispersion effects that occur in wave reflection, refraction, and mode conversion at the skull surface. The speed of the proposed framework is evaluated. We validate the accuracy of the framework using numerical phantoms and compare its results to k-Wave simulation results. Analytical validation is also performed based on the longitudinal and shear wave transmission coefficients. We then simulated, using our method, the major skull-distorting effects including amplitude attenuation, time-domain signal broadening, and time shift, and confirmed the findings by comparing them to several ex vivo experimental results. It is expected that the proposed method speeds up modeling and quantification of skull tissue and allows the development of transcranial photoacoustic brain imaging.
Keywords: acoustic attenuation; acoustic dispersion; mode conversion; numerical simulation; photoacoustic imaging; skull phantom; transcranial brain imaging.