Deep convolutional neural networks for estimating porous material parameters with ultrasound tomography

We study the feasibility of data based machine learning applied to ultrasound tomography to estimate water-saturated porous material parameters. In this work, the data to train the neural networks is simulated by solving wave propagation in coupled poroviscoelastic-viscoelastic-acoustic media. As the forward model, we consider a high-order discontinuous Galerkin method while deep convolutional neural networks are used to solve the parameter estimation problem. In the numerical experiment, we estimate the material porosity and tortuosity while the remaining parameters which are of less interest are successfully marginalized in the neural networks-based inversion. Computational examples confirms the feasibility and accuracy of this approach

Wave propagation in a fractional viscoelastic Andrade medium: diffusive approximation and numerical modeling

This study focuses on the numerical modeling of wave propagation in fractionally-dissipative media. These viscoelastic models are such that the attenuation is frequency dependent and follows a power law with non-integer exponent. As a prototypical example, the Andrade model is chosen for its simplicity and its satisfactory fits of experimental flow laws in rocks and metals. The corresponding constitutive equation features a fractional derivative in time, a non-local term that can be expressed as a convolution product which direct implementation bears substantial memory cost. To circumvent this limitation, a diffusive representation approach is deployed, replacing the convolution product by an integral of a function satisfying a local time-domain ordinary differential equation. An associated quadrature formula yields a local-in-time system of partial differential equations, which is then proven to be well-posed. The properties of the resulting model are also compared to those of the original Andrade model. The quadrature scheme associated with the diffusive approximation, and constructed either from a classical polynomial approach or from a constrained optimization method, is investigated to finally highlight the benefits of using the latter approach. Wave propagation simulations in homogeneous domains are performed within a split formulation framework that yields an optimal stability condition and which features a joint fourth-order time-marching scheme coupled with an exact integration step. A set of numerical experiments is presented to assess the efficiency of the diffusive approximation method for such wave propagation problems.Comment: submitted to Wave Motio

Biot-JKD model: simulation of 1D transient poroelastic waves with fractional derivatives

A time-domain numerical modeling of Biot poroelastic waves is presented. The viscous dissipation occurring in the pores is described using the dynamic permeability model developed by Johnson-Koplik-Dashen (JKD). Some of the coefficients in the Biot-JKD model are proportional to the square root of the frequency: in the time-domain, these coefficients introduce order 1/2 shifted fractional derivatives involving a convolution product. Based on a diffusive representation, the convolution kernel is replaced by a finite number of memory variables that satisfy local-in-time ordinary differential equations. Thanks to the dispersion relation, the coefficients in the diffusive representation are obtained by performing an optimization procedure in the frequency range of interest. A splitting strategy is then applied numerically: the propagative part of Biot-JKD equations is discretized using a fourth-order ADER scheme on a Cartesian grid, whereas the diffusive part is solved exactly. Comparisons with analytical solutions show the efficiency and the accuracy of this approach.Comment: arXiv admin note: substantial text overlap with arXiv:1210.036

Rigorous bounds for seismic dispersion and attenuation due to wave-induced fluid flow in porous rocks

The Hashin-Shtrikman (HS) bounds define the range of bulk and shear moduli of an elastic composite, given the moduli of the constituents and their volume fractions. Recently, the HS bounds have been ex tended to the quasi-static moduli of composite viscoelastic media. Because viscoelastic moduli are complex, the viscoelastic bounds form a closed curve on the complex plane. We analyze these general viscoelastic bounds for a particular case of a porous solid saturated with a Newtonian fluid. In our analysis, for poroelastic media, the viscoelastic bounds for the bulk modulus are represented by a semicircle and a segment of the real axis, connecting formal HS bounds that are computed for an inviscid fluid. Importantly, viscoelastic bounds for poroelastic media turn out to be independent of frequency. However, because the bounds are quasi-static, the frequency must be much lower than Biot’s characteristic frequency. Furthermore, we find that the bounds for the bulk modulus are attainable (realizable). We also find that these viscoelastic bounds account for viscous shear relaxation and squirt-flow dispersion, but do not account for Biot’s global flow dispersion, because the latter strongly depends on inertial forces