Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

Structure-from-motion photogrammetry analysis of historical aerial photography: Determining beach volumetric change over decadal scales

Historical aerial photographs are an invaluable tool in shoreline mapping and change detection in coastal landscapes. We evaluate the extent to which structure-from-motion (SfM) photogrammetric methods can be applied to quantify volumetric changes along sandy beaches, using archival imagery. We demonstrate the application of SfM-derived digital surface models (DSMs) at East Beach and Lady Bay in southwest Victoria, Australia, using photographic datasets taken in 1969, 1977 and 1986, and compare them to LiDAR-derived DSMs acquired at both sites in 2007. The SfM approaches resulted in two entire and two partial suitable DSMs out of six datasets. Good-quality DSMs were spatially continuous with a good spread of ground control points (GCPs) near the beach at Lady Bay, whereas unsuitable DSMs were mostly restricted by poor distribution and number of GCPs in spatially segmented areas of East Beach, due to limited overlapping of images, possible poor quality of GCPs and also the propagation of errors in the derived point clouds. A volume of approximately 223 000 ± 72 000 m³ was deposited at Lady Bay between 1969 and 2007, despite minimal erosion observed near the breakwater. The partially suitable dataset of East Beach indicated that beach erosion of at least 39 m³ m⁻¹ occurred immediately to the east of the seawall after 1977. We also discuss the drawbacks and strengths of SfM approaches as a benchmark of historical erosion assessments along sandy beaches. © 2020 John Wiley & Sons, Ltd.     

Styles of salt tectonics in the Sab’atayn basin,onshore Yemen

A variety of distinct salt tectonic features are present in the Sab’atayn Basin of western Yemen. Based on the interpretation of regional 2D seismic reflection data and exploration wells in the central part of the basin, an Upper Jurassic evaporite formation produced numerous salt rollers, salt pillows, reactive, flip-flop, and falling diapirs. Due to regional extension, halokinetics began as soon as the early Cretaceous, within just a few million years after the deposition of the Tithonian Sab’atayn evaporite sequence, by formation of salt rollers. The salt locally formed salt pillows which evolved to reactive and active salt diapirs and diapiric salt walls as the result of renewed, but low-strain extension in the basin. Some of the diapiric walls further evolved into falling diapirs due to ongoing extension. As the result of a prominent extensional episode at the end of the Cretaceous, many of the diapiric walls in the basin are controlled by large normal faults on their updip flanks. As the post-Cretaceous sedimentary cover is largely missing in the study area, the assumed reactivation of salt structures during the Cenozoic remains poorly constrained. The interpreted changes in the style of salt tectonics in the Sab’atayn Basin offer a better understanding of the regional-scale tectonic development of the Arabian plate during the late Jurassic and Cretaceous.     

Influence of damage on pore size distribution and permeability of rocks

The model proposed in this article relates permeability to porosity measurements that can easily be performed in the laboratory. The pore size distribution (PSD) curve is updated with strains and damage. The updated volumetric fractions of natural pores and cracks are introduced in the expression of permeability. Contrary to classical permeability models based on PSD integrations, the model proposed in this article accounts for possible changes in the porosity modes: one mode for undamaged samples and two modes for cracked samples. The proposed approach also accounts for varying states of damage, as opposed to classical fracture network models, in which the cracks pattern is fixed. The only material parameters that are required to describe the microstructure are the lower and upper bounds of the pores size for both natural pores and cracks. All the other PSD parameters involved in the model are related to macroscopic parameters that can easily be determined in the laboratory, such as the initial void ratio. The framework proposed in this article can be used in any damage constitutive model to determine the permeability of a brittle porous medium. Drained triaxial compression tests have been simulated. Before cracks initiation, permeability decreases while the larger natural pores are getting squeezed. After the occurrence of damage, permeability grows due to the increase of cracks density. The model performs well to represent the influence of the confining pressure on damage evolution and permeability variations. Copyright © 2012 John Wiley & Sons, Ltd.     

Identifying oceanographic conditions conducive to coastal impacts on temperate open coastal beaches

Natural Hazards - Warnings issued by meteorological or oceanographic agencies are a common means of allowing people to prepare for likely impactful events. Quantifying the relationships between...     

Micromechanical modeling for rate‐dependent behavior of salt rock under cyclic loading

The dependence of rock behavior on the deformation rate is still not well understood. In salt rock, the fundamental mechanisms that drive the accumulation of irreversible deformation, the reduction of stiffness, and the development of hysteresis during cyclic loading are usually attributed to intracrystalline plasticity and diffusion. We hypothesize that at low pressure and low temperature, the rate‐dependent behavior of salt rock is governed by water‐assisted diffusion along grain boundaries. Accordingly, a chemo‐mechanical homogenization framework is proposed in which the representative elementary volume (REV) is viewed as a homogeneous polycrystalline matrix that contains sliding grain‐boundary cracks. The slip is related to the mass of salt ions that diffuse along the crack surface. The relationship between fluid inclusion‐scale and REV‐scale stresses and strains is established by using the Mori–Tanaka homogenization scheme. It is noted from the model that a lower strain rate and a larger number of sliding cracks enhance stiffness reduction and hysteresis. Thinner sliding cracks (i.e., thinner brine films) promote stiffness reduction and accelerate stress redistributions. The larger the volume fraction of the crack inclusions, the larger the REV deformation and the larger the hysteresis. Results presented in this study shed light on the mechanical behavior of salt rock that is pertinent to the design of geological storage facilities that undergo cyclic unloading, which could help optimize the energy production cycle with low carbon emissions.