Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : Production rates in horizontal shale gas wells display declines that are controlled by the low permeability and the intrinsic heterogeneity of the shale matrix. We present an original multi-continuum approach that yields a physical model able to reproduce the complexity of the decreasing gas rates. The model describes the dynamics of gas rate as function of the physical reservoir parameters and geometry, while the shale matrix heterogeneity is accounted for by a stochastic description of transmissivity field. From the 3D (Dimensional) problem setting, including the heterogeneous shale matrix, the fractures generated by the hydrofracking operations, as well as the production well characteristics, we establish an effective upscaled 1D model for the gas pressures in fracture and matrix as well as the volumetric flux. We analyse the decline curves behaviour, and we identify the time scales that characterize the dynamics of the gas rate decline using explicit analytical Laplace space solutions of the upscaled process model. Asymptotically, the flux curves decrease exponentially, while in an intermediate regime we find a power-law behaviour, in which the flux scales with a power law in time as (t^{-beta }), where (beta ) reflects the medium heterogeneity. We use this solution to fit a set of real data displaying distinctly different decline trends and study the sensitivity of the model to the reservoir parameters in order to identify their respective controls at the different stages of the decline curve dynamics. Results indicate that the initial value of the gas rate is determined by the transmissivity of the fractures and the initial pressure of the gas in the shale matrix. The latter causes mainly a shift in the entire decline curve. The early time of decline curve shape is primarily controlled by the fracture properties (compressibility and transmissivity). During the main part of the economically valuable production times, i.e. before the production rate drops exponentially, the decline curve is strongly controlled by the properties of the shale rocks including their heterogeneity, which is modelled by two parameters describing the non-Fickian pressure diffusion effects in a stochastic framework.
ISSN: 0169-3913

hal-01217247
https://hal.archives-ouvertes.fr/hal-01217247
DOI : 10.1007/s11242-015-0504-y

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In petroleum industry, the difference between pore pressure (Pp) and minimum horizontal stress Sh (termed the seal or retention capacity) is of major consideration because it is often assumed to represent how close a system is to hydraulic failure and thus the maximum hydrocarbon column height that can be maintained. While Sh and Pp are often considered to be independent parameters, several studies in the last decade have demonstrated that Sh and Pp are in fact coupled. However, the nature of this coupling relationship remains poorly understood. In this paper, we explore the influences of the spatial pore pressure distribution on Sh/Pp coupling and then on failure pressure predictions and trap integrity evaluation. With analytical models, we predict the fluid pressure sustainable within a reservoir before failure of its overpressured shale cover. We verify our analytical predictions with experiments involving analogue materials and fluids. We show that hydraulic fracturing and seal breach occur for fluid pressure greater than it would be expected from conventional retention capacity. This can be explained by the impact of the fluid overpressure field in the overburden and the pressure diffusion around the reservoir on the principal stresses. We calculate that supralithostatic pressure could locally be reached in overpressured covers. We also define the retention capacity of a cover (RC) surrounding a fluid source or reservoir as the difference between the failure pressure and the fluid overpressure prevailing in shale at the same depth. In response to a localized fluid pressure rise, we show that the retention capacity does not only depend on the pore fluid overpressure of the overburden but also on the tensile strength of the cover, its Poisson's ratio, and the depth and width of the fluid source.
ISSN: 0264-8172

hal-00590497
https://hal.archives-ouvertes.fr/hal-00590497
DOI : 10.1016/j.marpetgeo.2010.08.007

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Sedimentary rocks and detrital minerals sample large areas of the continental crust, and they are increasingly seen as a reliable archive for its global evolution. This study presents two approaches to model the growth of the continental crust through the sedimentary archive. The first builds on the variations in U-Pb, Hf and O isotopes in global databases of detrital zircons. We show that uncertainty in the Hf isotope composition of the mantle reservoir from which new crust separated, in the 176Lu/177Hf ratio of that new crust, and in the contribution in the databases of zircons that experienced ancient Pb loss(es), adds some uncertainty to the individual Hf model ages, but not to the overall shape of the calculated continental growth curves. The second approach is based on the variation of Nd isotopes in 645 worldwide fine-grained continental sedimentary rocks with different deposition ages, which requires a correction of the bias induced by preferential erosion of younger rocks through an erosion parameter referred to as K. This dimensionless parameter relates the proportions of younger to older source rocks in the sediment, to the proportions of younger to older source rocks present in the crust from which the sediment was derived. We suggest that a Hadean/Archaean value of K = 1 (i.e., no preferential erosion), and that post-Archaean values of K = 4–6, may be reasonable for the global Earth system. Models built on the detrital zircon and the fine-grained sediment records independently suggest that at least 65% of the present volume of continental crust was established by 3 Ga. The continental crust has been generated continuously, but with a marked decrease in the growth rate at ~ 3 Ga. The period from > 4 Ga to ~ 3 Ga is characterised by relatively high net rates of continental growth (2.9–3.4 km3 yr− 1 on average), which are similar to the rates at which new crust is generated (and destroyed) at the present time. Net growth rates are much lower since 3 Ga (0.6–0.9 km3 yr− 1 on average), which can be attributed to higher rates of destruction of continental crust. The change in slope in the continental growth curve at ~ 3 Ga is taken to indicate a global change in the way bulk crust was generated and preserved, and this change has been linked to the onset of subduction-driven plate tectonics. At least 100% of the present volume of the continental crust has been destroyed and recycled back into the mantle since ~ 3 Ga, and this time marks a transition in the average composition of new continental crust. Continental crust generated before 3 Ga was on average mafic, dense, relatively thin (< 20 km) and therefore different from the calc-alkaline andesitic crust that dominates the continental record today. Continental crust that formed after 3 Ga gradually became more intermediate in composition, buoyant and thicker. The increase in crustal thickness is accompanied by increasing rates of crustal reworking and increasing input of sediment to the ocean. These changes may have been accommodated by a change in lithospheric strength at around 3 Ga, as it became strong enough to support high-relief crust. This time period therefore indicates when significant volumes of continental crust started to become emergent and were available for erosion and weathering, thus impacting on the composition of the atmosphere and the oceans.
ISSN: 0037-0738

hal-01622696
https://hal.archives-ouvertes.fr/hal-01622696
DOI : 10.1016/j.sedgeo.2017.06.001

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Phyllosilicate content and related permeability of fault zones form primary controls on hydraulic and mechanical behavior of the brittle crust. Hence, understanding and predicting the localization of these ubiquitous minerals is a major issue for fundamental and applied geosciences. We describe normal fault zones cutting a foreland arkosic turbiditic formation suffering high-T diagenesis and formed under conditions (~200°C) typical of deeply buried reservoirs and common within the seismogenic interval. Microstructural analyses show a large proportion of phyllosilicates (up to 34%) in the fault rock, derived from near-complete feldspar alteration and disaggregation during deformation. This study shows that even faults with offsets (~20 cm) much lower than bed thickness can have such large feldspar-to-phyllosilicate transformation ratios, implying that the origin of the phyllosilicates is purely transformation related. These results imply that the potential sealing capacity and strength of faults could be predicted from the host rock feldspar content. Where sealing capacity and fault strength can be related to phyllosilicate content, these properties can then also be inferred from the predicted phyllosilicate content: this opens up new horizons concerning the hydraulic and the mechanical behavior of the upper crust.
ISSN: 0094-8276

hal-00903574
https://hal.archives-ouvertes.fr/hal-00903574
DOI : 10.1002/grl.50829