Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEffective elasticity of a medium with many parallel fractures
87
0
0.0
(
0
)
Authors
Filip P. Adamus
Ask ChatGPT about the research
No Arabic abstract
We consider an alternative way of obtaining the effective elastic properties of a cracked medium. Similarly, to the popular linear-slip model, we assume flat, parallel fractures, and long wavelengths. However, we do not treat fractures as weakness planes of displacement discontinuity. In contrast to the classical models, we represent fractures by a thin layer embedded in the background medium. In other words, we follow the Schoenberg-Douma matrix formalism for Backus averaging, but we relax their assumptions of infinite weakness and marginal thickness of a layer so that it does not correspond to the linear-slip plane. To represent the properties of a fracture, we need a fourth order elasticity tensor and a thickness parameter. The effective tensor becomes more complicated, but it may describe a higher concentration of parallel cracks more accurately. Apart from the derivations of the effective elasticity tensors, we perform numerical experiments in which we compare the performance of our approach with a linear-slip model in the context of highly fractured media. Our model becomes pertinent if filled-in cracks occupy more than one percent of the effective medium.
rate research
Read More
Phase-field modeling -- a continuous approach to discontinuities -- is gaining popularity for simulating rock fractures due to its ability to handle complex, discontinuous geometry without an explicit surface tracking algorithm. None of the existing phase-field models, however, incorporates the impact of surface roughness on the mechanical response of fractures -- such as elastic deformability and shear-induced dilation -- despite the importance of this behavior for subsurface systems. To fill this gap, here we introduce the first framework for phase-field modeling of rough rock fractures. The framework transforms a displacement-jump-based discrete constitutive model for discontinuities into a strain-based continuous model, and then casts it into a phase-field formulation for frictional interfaces. We illustrate the framework by constructing a particular phase-field form employing a rock joint model originally formulated for discrete modeling. The results obtained by the new formulation show excellent agreement with those of a well-established discrete method for a variety of problems ranging from shearing of a single discontinuity to compression of fractured rocks. Consequently, our phase-field framework provides an unprecedented bridge between a discrete constitutive model for rough discontinuities -- common in rock mechanics -- and the continuous finite element method -- standard in computational mechanics -- without any algorithm to explicitly represent discontinuity geometry.
In this paper, we consider a long-wave equivalent medium to a finely parallel-layered inhomogeneous medium, obtained using the Backus average. Following the work of Postma and Backus, we show explicitly the derivations of the conditions to obtain the equivalent isotropic medium. We demonstrate that there cannot exist a transversely isotropic (TI) equivalent medium with the coefficients $c^{overline{rm TI}}_{1212} eq c^{overline{rm TI}}_{2323}$, $c^{overline{rm TI}}_{1111} = c^{overline{rm TI}}_{3333}$ and $c^{overline{rm TI}}_{1122} = c^{overline{rm TI}}_{1133}$. Moreover, we consider a new parameter, $varphi$, indicating the anisotropy of the equivalent medium, and we show its range and properties. Subsequently, we compare $varphi$ to the Thomsen parameters, emphasizing its usefulness as a supportive parameter showing the anisotropy of the equivalent medium or as an alternative parameter to $delta$. We argue with certain Berryman et al. considerations regarding the properties of the anisotropy parameters $epsilon$ and $delta$. Additionally, we show an alternative way---to the one mentioned by Berryman et al.---of indicating changing fluid content in layered Earth.
Hydraulic tension fractures were produced in porous limestones using a specially designed hydraulic cell. The 3D geometry of the samples was imaged using X-ray computed microtomography before and after fracturation. Using these data, it was possible to estimate the permeability tensor of the core samples, extract the path of the rupture and compare it to the heterogeneities initially present in the rock.
While the hydraulic fracturing technology, aka fracking (or fraccing, frac), has become highly developed and astonishingly successful, a consistent formulation of the associated fracture mechanics that would not conflict with some observations is still unavailable. It is attempted here. Classical fracture mechanics, as well as the current commercial softwares, predict vertical cracks to propagate without branching from the perforations of the horizontal well casing, which are typically spaced at 10 m or more. However, to explain the gas production rate at the wellhead, the crack spacing would have to be only about 0.1 m, which would increase the overall gas permeability of shale mass about 10,000$times$. This permeability increase has generally been attributed to a preexisting system of orthogonal natural cracks, whose spacing is about 0.1 m. But their average age is about 100 million years, and a recent analysis indicated that these cracks must have been completely closed by secondary creep of shale in less than a million years. Here it is considered that the tectonic events that produced the natural cracks in shale must have also created weak layers with nano- or micro-cracking damage. It is numerically demonstrated that a greatly enhanced permeability along the weak layers, with a greatly increased transverse Biot coefficient, must cause the fracking to engender lateral branching and the opening of hydraulic cracks along the weak layers, even if these cracks are initially almost closed. A finite element crack band model, based on recently developed anisotropic spherocylindrical microplane constitutive law, demonstrates these findings.
We discuss several physical constraints imposed on elasticity parameters of a transversely isotropic (TI) tensor. There are three types of restrictions we investigate; a fundamental one of stability conditions, and two additional ones, commonly considered in seismology. The first commonly considered restriction comes from an assumption of a wave with a greater speed in the horizontal than vertical direction. The second constitute the assumption that quasi-P wave is faster than quasi-S waves. We show several numerical examples to examine how these restrictions affect a TI tensor with known values of certain elasticity constants that could be acquired from the vertical or horizontal measurements.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
