Louis De Barros

Wave propagation in heterogeneous porous media formulated in the frequency-space domain using a discontinuous Galerkin method

Biphasic media with a dynamic interaction between fluid and solid phases must be taken into account to accurately describe seismic wave amplitudes in subsurface and reservoir geophysical applications. Consequently, the modeling of the wave propagation in heteregeneous porous media, which includes the frequency-dependent phenomena of the fluidsolid interaction, is considered for 2D geometries. From the Biot-Gassmann theory, we have deduced the discrete linear system in the frequency domain for a discontinuous finite-element method, known as the nodal discontinuous Galerkin method. Solving this system in the frequency domain allows accurate modeling of the Biot wave in the diffusive/propagative regimes, enhancing the importance of frequency effects. Because we had to consider finite numerical models, we implemented perfectly matched layer techniques. We found that waves are efficiently absorbed at the model boundaries, and that the discretization of the medium should follow the same rules as in the elastodynamic case, that is, 10 grids per minimum wavelength for a P0 interpolation order. The grid spreading of the sources, which could be stresses or forces applied on either the solid phase or the fluid phase, did not show any additional difficulties compared to the elastic problem. For a flat interface separating two media, we compared the numerical solution and a semianalytic solution obtained by a reflectivity method in the three regimes where the Biot wave is propagative, diffusive/propagative, and diffusive. In all cases, fluid-solid interactions were reconstructed accurately, proving that attenuation and dispersion of the waves were correctly accounted for. In addition to this validation in layered media, we have explored the capacities of modeling complex wave propagation in a laterally heterogeneous porous medium related to steam injection in a sand reservoir and the seismic response associated to a fluid substitution.

Mega-earthquakes rupture flat megathrusts

Mega-earthquakes go the flat way Megathrust faults in subduction zones cause large and damaging earthquakes. Bletery et al. argue that certain geometric features of the subduction zones relate to earthquake size. The key parameter is the curvature of the megathrust. Larger earthquakes occur where the subducting slab is flatter, providing a rough metric for estimating where mega-earthquakes may occur in the future. Science, this issue p. 1027 Large earthquakes in subduction zones are most likely to occur where the subducting slab is relatively flat. The 2004 Sumatra-Andaman and 2011 Tohoku-Oki earthquakes highlighted gaps in our understanding of mega-earthquake rupture processes and the factors controlling their global distribution: A fast convergence rate and young buoyant lithosphere are not required to produce mega-earthquakes. We calculated the curvature along the major subduction zones of the world, showing that mega-earthquakes preferentially rupture flat (low-curvature) interfaces. A simplified analytic model demonstrates that heterogeneity in shear strength increases with curvature. Shear strength on flat megathrusts is more homogeneous, and hence more likely to be exceeded simultaneously over large areas, than on highly curved faults.

Source geometry from exceptionally high resolution Long Period event observations at Mt Etna during the 2008 eruption.

During the second half of June, 2008, 50 broadband seismic stations were deployed on Mt Etna volcano in close proximity to the summit, allowing us to observe seismic activity with exceptionally high resolution. 129 long period events (LP) with dominant frequencies ranging between 0.3 and 1.2 Hz, were extracted from this dataset. These events form two families of similar waveforms with different temporal distributions. Event locations are performed by cross-correlating signals for all pairs of stations in a two-step scheme. In the first step, the absolute location of the centre of the clusters was found. In the second step, all events are located using this position. The hypocentres are found at shallow depths (20 to 700 m deep) below the summit craters. The very high location resolution allows us to detect the temporal migration of the events along a dike-like structure and 2 pipe shaped bodies, yielding an unprecedented view of some elements of the shallow plumbing system at Mount Etna. These events do not seem to be a direct indicator of the ongoing lava flow or magma upwelling.

Perturbations of the seismic reflectivity of a fluid-saturated depth-dependent poroelastic medium

Analytical formulas are derived to compute the first-order effects produced by plane inhomogeneities on the point source seismic response of a fluid-filled stratified porous medium. The derivation is achieved by a perturbation analysis of the poroelastic wave equations in the plane-wave domain using the Born approximation. This approach yields the Frechet derivatives of the P-SV- and SH-wave responses in terms of the Greens functions of the unperturbed medium. The accuracy and stability of the derived operators are checked by comparing, in the time-distance domain, differential seismograms computed from these analytical expressions with complete solutions obtained by introducing discrete perturbations into the model properties. For vertical and horizontal point forces, it is found that the Frechet derivative approach is remarkably accurate for small and localized perturbations of the medium properties which are consistent with the Born approximation requirements. Furthermore, the first-order formulation appears to be stable at all source-receiver offsets. The porosity, consolidation parameter, solid density, and mineral shear modulus emerge as the most sensitive parameters in forward and inverse modeling problems. Finally, the amplitude-versus-angle response of a thin layer shows strong coupling effects between several model parameters.

Crustal structure below Popocatépetl Volcano (Mexico) from analysis of Rayleigh waves

An array of ten broadband stations was installed on the Popocatepetl volcano (Mexico) for five months between October 2002 and February 2003. 26 regional and teleseismic earthquakes were selected and filtered in the frequency time domain to extract the fundamental mode of the Rayleigh wave. The average dispersion curve was obtained in two steps. Firstly, phase velocities were measured in the period range [2-50] s from the phase difference between pairs of stations, using Wiener filtering. Secondly, the average dispersion curve was calculated by combining observations from all events in order to reduce diffraction effects. The inversion of the mean phase velocity yielded a crustal model for the volcano which is consistent with previous models of the Mexican Volcanic Belt. The overall crustal structure beneath Popocatepetl is therefore not different from the surrounding area, and the velocities in the lower crust are confirmed to be relatively low. Lateral variations of the structure were also investigated by dividing the network into four parts and by applying the same procedure to each sub-array. No well-defined anomalies appeared for the two sub-arrays for which it was possible to measure a dispersion curve. However, dispersion curves associated with individual events reveal important diffraction for 6 s to 12 s periods which could correspond to strong lateral variations at 5 to 10 km depth.

A brittle failure model for long‐period seismic events recorded at Turrialba Volcano, Costa Rica

A temporary seismic network, consisting of 23 broadband and six short-period stations, was installed in a dense network at Turrialba Volcano, Costa Rica, between 8 March and 4 May 2011. During this time 513 long-period (LP) events were observed. Due to their pulse-like waveforms, the hypothesis that the events are generated by a slow-failure mechanism, based on a recent new model by Bean et al. (2014), is tested. A significant number (107) of the LPs are jointly inverted for their source locations and mechanisms, using full-waveform moment tensor inversion. The locations are mostly shallow, with depths < 800 m below the active Southwest Crater. The results of the decompositions of the obtained moment tensor solutions show complex source mechanisms, composed of high proportions of isotropic and low, but seemingly significant, proportions of compensated linear vector dipole and double-couple components. It is demonstrated that this can be explained as mode I tensile fracturing with a strong shear component. The source mechanism is further investigated by exploring scaling laws within the data. The LPs recorded follow relationships very similar to those of conventional earthquakes, exhibiting frequency-magnitude and corner frequency versus magnitude relationships that can be explained by brittle failure. All of these observations indicate that a slow-failure source model can successfully describe the generation of short-duration LP events at Turrialba Volcano.

Seismic velocity changes associated with aseismic deformations of a fault stimulated by fluid injection

Fluid pressure plays an important role in the stability of tectonic faults. However, the in situ mechanical response of faults to fluid-pressure variations is still poorly known. To address this question, we performed a fluid-injection experiment in a fault zone in shales while monitoring fault movements at the injection source and seismic velocity variations from a near-distance (<10 meters) monitoring network. We measured and located the P- and S-wave velocity perturbations in and around the fault using repetitive active sources. We observed that seismic velocity perturbations dramatically increase above 1.5 MPa of injection pressure. This is consistent with an increase of fluid flow associated with an aseismic dilatant shearing of the fault as shown by numerical modelling. We find that seismic velocity changes are sensitive both to fault opening by fluid invasion and effective stress variations, and can be an efficient measurement for monitoring fluid-driven aseismic deformations of faults.

Aseismic Motions Drive a Sparse Seismicity During Fluid Injections Into a Fractured Zone in a Carbonate Reservoir

Author(s): Duboeuf, L; De Barros, L; Cappa, F; Guglielmi, Y; Deschamps, A; Seguy, S | Abstract: ©2017. American Geophysical Union. All Rights Reserved. An increase in fluid pressure in faults can trigger seismicity and large aseismic motions. Understanding how fluid and faults interact is an essential goal for seismic hazard and reservoir monitoring, but this key relation remains unclear. We developed an in situ experiment of fluid injections at a 10 meter scale. Water was injected at high pressure in different geological structures inside a fault damaged zone, in limestone at 280 m depth in the Low Noise Underground Laboratory (France). Induced seismicity, as well as strains, pressure, and flow rate, was continuously monitored during the injections. Although nonreversible deformations related to fracture reactivations were observed for all injections, only a few tests generated seismicity. Events are characterized by a 0.5-to-4 kHz content and a small magnitude (approximately −3.5). They are located within 1.5 m accuracy between 1 and 12 m from the injections. Comparing strain measurements and seismicity shows that more than 96% of the deformation is aseismic. The seismic moment is also small compared to the one expected from the injected volume. Moreover, a dual seismic behavior is observed as (1) the spatiotemporal distribution of some cluster of events is clearly independent from the fluid diffusion (2) while a diffusion-type pattern can be observed for some others clusters. The seismicity might therefore appear as an indirect effect to the fluid pressure, driven by aseismic motion and related stress perturbation transferred through failure.

Using a Poroelastic Theory to Reconstruct Subsurface Properties: Numerical Investigation

The quantitative imaging of the Earth subsurface is a major challenge in geophysics. In oil and gas exploration and production, aquifer management and other applications such as the underground storage of CO2 , seismic imaging techniques are implemented to provide as much information as possible on fluid-filled reservoir rocks. Biot theory (Biot, 1956) and its extensions provide a convenient framework to connect the various parameters characterizing a porous medium to the wave properties, namely, their amplitudes, velocities and frequency contents. The poroelastic model involves more parameters than the elastodynamic theory, but on the other hand, the wave attenuation and dispersion characteristics at the macroscopic scale are determined by the intrinsic properties of the medium without having to resort to empirical relationships.

Relocation of Long Period (LP) seismic events reveals en echelon fractures in the upper edifice of Turrialba volcano, Costa Rica

Swarms of long-period (LP) events were recorded on Turrialba volcano, Costa Rica, during a seismic field experiment in 2009. Families of LP events were previously identified and located using a joint inversion for source location and mechanism; however the spatial resolution of the obtained locations was not sufficient for imaging the structures on which they occur. Using a waveform similarity-based location method, we take advantage of the joint location-mechanism inversion by relocating events around the obtained familial location. The location method is successfully tested on a synthetic data set, and is then applied to the Turrialba LP data set. The relocated events are jointly interpreted with their source mechanisms, and reveal an en echelon structure within the upper-edifice of the volcano. This can be interpreted as a response of a shearing band with high fluid pressure inducing tensile fractures at unconsolidated rock layer interfaces within the upper-edifice of the volcano.