Search NASA⌕ Search

Engineering topics

Masson, Yder

Publications and source records attributed to Masson, Yder.

Box tomography: localized imaging of remote targets buried in an unknown medium, a step forward for understanding key structures in the deep Earth

We discuss an application of tomographic methods to the imaging of localized geological structures buried inside the Earth, where neither seismic sources nor receivers are necessarily present. We call this ‘box tomography’. The essential difference between box-tomography and standard tomographic methods is that the numerical modelling (i.e. the ray tracing in traveltime tomography and the wave propagation in waveform tomography or full waveform inversion) is completely confined within the small box-region imaged. Thus, box tomography is a lot more efficient than global tomography (i.e. where we invert for the velocity in the larger volume that encompasses all the sources and receivers), for imaging localized objects. The principal difficulty in the box tomography problem comes from the medium outside the box-region imaged, that masks the target structure. In the data recorded at the surface, the information carried by the waves about the target structure is mixed with that about the surrounding medium. The main focus of this paper is to evaluate how well box tomography can extract useful information to produce accurate images of remote buried objects, especially in the difficult situation where the velocity distribution above the target structure is not known a priori. We implemented traveltime box tomography, where first arrival times are inverted using ray theory, and, waveform box tomography or full waveform box inversion, where seismograms are inverted using either scattering integrals or adjoint methods. Whatever the method employed, we found that, in the simple situation where the velocity structure is known around the target region, box-tomography is faster and more accurate than global tomography. Most importantly, through examples and simple analysis, we show that box-tomography can produce accurate images of localized structures even though the velocity distribution is not a priori known around the box-region imaged. In the last section of this paper, we evaluate the usage of box-tomography for imaging geological structures in the deep Earth. We performed 2-D synthetic experiments in a fluid planet and we argue that box tomography can be employed for refining our current global tomographic models within arbitrary regions inside the Earth. This has important implications for understanding the composition and dynamics of the Earth’s interior.

Masson, Yder↗

Fast computation of synthetic seismograms within a medium containing remote localized perturbations: a numerical solution to the scattering problem

We derive a fast discrete solution to the scattering problem. This solution allows us to compute accurate synthetic seismograms or waveforms for arbitrary locations of sources and receivers within a medium containing localized perturbations. The key to efficiency is that wave propagation modelling does not need to be carried out in the entire volume that encompasses the sources and the receivers but only within the subvolume containing the perturbations or scatterers. The proposed solution has important applications, for example, it permits the imaging of remote targets located in regions where no sources or receivers are present. Our solution relies on domain decomposition: within a small volume that contains the scatterers, wave propagation is modelled numerically, while in the surrounding volume, where the medium isn't perturbed, the response is obtained through wavefield extrapolation. The originality of this work is the derivation of discrete formulas for representation theorems and Kirchhof–Helmholtz integrals that naturally adapt to the numerical scheme employed for modelling wave propagation. Our solution applies, for example, to finite difference methods or finite/spectral elements methods. The synthetic seismograms obtained with our solution can be considered ‘exact’ as the total numerical error is comparable to that of the method employed for modelling wave propagation. We detail a basic implementation of our solution in the acoustic case using the finite difference method and present numerical examples that demonstrate the accuracy of the method. We show that ignoring some terms accounting for higher order scattering effects in our solution has a limited effect on the computed seismograms and significantly reduces the computational effort. Finally, we show that our solution can be used to compute localized sensitivity kernels and we discuss applications to target oriented imaging. Extension to the elastic case is straightforward and summarized in a dedicated section.

97 MATHEMATICS AND COMPUTING↗

Mapping the mechanical properties of rocks using automated microindentation tests

A microindentation scanner is constructed that measures the spatial fluctuation in the elastic properties of natural rocks. This novel instrument performs automated indentation tests on the surface of a rock slab and outputs 2-D maps of the indentation modulus at submillimeter resolution. Maps obtained for clean, well-consolidated, sandstone are presented and demonstrate the capabilities of the instrument. We postulate that the elastic structure of sandstones correlates well with their visual appearance. Moreover, we observe that the probability distribution of the indentation modulus fluctuations across the slab surfaces can be modeled using a lognormal probability density function. To illustrate possible use of the data obtained with the microindentation scanner, we use roughly 10 cm × 10 cm scans with millimeter resolution over four sandstone planar slabs to numerically compute the overall drained elastic moduli for each sandstone sample. We demonstrate that such numerically computed moduli are well modeled using the multicomponent form of the Hashin-Shtrikman lower bound that employs the observed lognormal probability distribution for the mesoscopic-scale moduli (the geometric mean works almost the same). We also compute the seismic attenuation versus frequency associated with wave-induced fluid flow between the heterogeneities in the scanned sandstones and report relatively small values for the inverse quality factor (Q -1 <10 -2 ) in the seismic frequency band 10 2 Hz 4 Hz. The numerically computed frequency dependence in the attenuation varies from one type of sandstone to another, and we determine vital anisotropy in the attenuation associated with waves propagating in different directions.

58 GEOSCIENCES↗