Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
An analytical solution in the Laplace domain is derived for modeling anomalous pressure diffusion during pumping tests in aquifer displaying hierarchical fractal fracture networks. The proposed solution generalizes all of the analytical models for fractal flow published previously by combining multifractal diffusion and nested multiporosity with transient exchanges, interface skin effects, and well storage effect. Solutions are derived for fracture-delimited blocks with planar, cylindrical, and spherical shapes, as well as with any fractional dimensional shapes. Any combinations of these shapes can be defined in order to model a large range of situations. Within each permeability level, the fractal properties of the fracture network can be specified and the fractal dimension can be distinct from the shape dimension of the block.
ISSN: 0043-1397

hal-00411549
https://hal.archives-ouvertes.fr/hal-00411549
DOI : 10.1029/2008WR007125

Éditeur(s) : HAL CCSD
Résumé : International audience
The Corinth Rift is the most seismically active zone in Europe. The area is characterized by very localized NSextension at a rate of1.5cm/year, the occurrence of frequent and intensive microseismic crises and occasionalmoderate to large earthquakes like in 1995 (Mw=6.1).Since the year 2000, the Corinth Rift Laboratory (CRL, http://crlab.eu) consisting in a multidisciplinary naturalobservatory, aims at understanding the mechanics of faulting and earthquake nucleation in the Rift. Recent studieshave improved our view about fault geometry and mechanics within CRL, but there is still a critical need for a betterknowledge of the structure at depth both for the accuracy of earthquake locations and for mechanical interpretationof the seismicity.In this project, we aim to analyze the complete seismological database (13 years of recordings) of CRL by usingrecently developed methodologies of structural imaging, in order to determine at the same time and with highresolution, the local 3D structure and the earthquake locations.We perform an iterative joint determination of 3D velocity model and earthquake coordinates. In a first step, Pand S velocity models are determined using first arrival time tomography method proposed by Taillandier et al.(2009). It consists in the minimization of the cost function between observed and theoretical arrival times which isachieved by the steepest descent method (e.g. Tarantola 1987). This latter requires computing the gradient of thecost function by using the adjoint state method (Chavent 1974). In a second step, earthquakes are located in thenew velocity model with a non-linear inversion method based on a Bayesian formulation (Gesret et al. 2015). Step1 and 2 are repeated until the cost function no longer decreases.We present preliminary results consisting in: (1) the adjustement of a 1D velocity model that is used as initialmodel of the 3D tomography and (2) a first attempt of the joint determination of 3D velocity model and earthquakelocation.In addition to the tomographic imaging, we perform a preliminary receiver function analysis of teleseismic datarecorded by the broadband stations of the CRL network. The RF analysis should provide the interface depthsbeneath seismometers and increase the imaging resolution of the upper crustal structures provided by the 3Dtomography. In this first attempt, we adjust the 1D velocity model that produces a synthetic RF as similar aspossible to the observed RF for a subset of data. We compare the identified interfaces with structures imaged bythe tomography.
EGU general assembly 2016
Vienne, Austria

hal-01306224
https://hal-mines-paristech.archives-ouvertes.fr/hal-01306224

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The main goal of the monitoring stage of Carbon Capture and Storage (CCS) is to obtain an accurate estimation of the subsurface CO2 accumulation and to detect any possible leakage. Laboratory experiments are necessary to investigate the small scale processes governing the CO2–brine–rock interaction. They also provide a means to calibrate the results coming from field scale geophysical methods. In this work we set up an experimental system which is able to perform Electrical Resistivity Tomography (ERT) measurements on centimeter-scale rock samples at various P–T conditions. We present the results of two new experiments related to CO2 monitoring, performed on a cylindrical (4 × 8 cm) Fontainebleau rock sample. In the first one, we have quantified the CO2 saturation at different volume fractions, representing zones from a deep saline aquifer with varying degrees of saturation. In the second one, we have monitored and quantified the effect of CO2 dissolution in the brine at a pressure of 40 bar during eight days, emulating the invasion of CO2 into a shallow aquifer. Results highlight the importance of accounting for the contribution of surface conductivity in highly CO2-saturated regions, even in clay-free rocks, and also for brine conductivity variation due to CO2 dissolution. Ignoring any of these effects will end up in a CO2 saturation underestimation. We present a modified CO2 saturation equation to account for these two influences.
ISSN: 0926-9851

hal-01356323
https://hal.archives-ouvertes.fr/hal-01356323
DOI : 10.1016/j.jappgeo.2016.03.039

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We report measurements of long-period P wave polarization (P-pol) in Australia and Tahiti made by combining modeling of the polarization deviation and harmonic analysis. The analysis of the deviation of the horizontal polarization of the P wave as a function of event back azimuth may be used to obtain information about (1) sensor misorientation, (2) dipping discontinuities, (3) seismic anisotropy, and (4) velocity heterogeneities beneath a seismic station. The results from harmonic analysis and a grid search using Snell's law suggest the presence of a dipping seismic discontinuity beneath stations CTAO and CAN in Australia. These results are consistent with published receiver function studies for these stations. The P-pol fast axis orientation is close to the N-S absolute plate motion direction at station TAU (Tasmania), which may be due to plate-motion-driven alignment of olivine crystals in the asthenosphere. Interestingly, measurements of SKS splitting at Tahiti (French Polynesia) show an apparent isotropy, whereas an inversion of P-pol observations at PPTL seismic station located in Tahiti suggests the presence of two anisotropic layers. The fast axis azimuth is oriented E-Win the upper layer, and it is close to the NW-SE orientation in the lower layer. Since P-pol orientations are used for real-time earthquake locations, especially in poorly instrumented areas such as the South Pacific, we show that the bias from anisotropy and sensor misorientation determined here can be corrected to improve the location accuracy, which yields fundamental data for rapid location necessary for effective tsunami warning.
ISSN: 2169-9313

hal-00413141
https://hal.archives-ouvertes.fr/hal-00413141
https://hal.archives-ouvertes.fr/hal-00413141/document
https://hal.archives-ouvertes.fr/hal-00413141/file/FontaineetalJGR2008.pdf
DOI : 10.1029/2008JB005709