Poroelastic Model Research Articles

4D Finite element modeling of stress distribution in depleted reservoir of south Iraq oilfield

The harvest of hydrocarbon from the depleted reservoir is crucial during field development. Therefore, drilling operations in the depleted reservoir faced several problems like partial and total lost circulation. Continuing production without an active water drive or water injection to support reservoir pressure will decrease the pore and fracture pressure. Moreover, this depletion will affect the distribution of stress and change the mud weight window. This study focused on vertical stress, maximum and minimum horizontal stress redistributions in the depleted reservoirs due to decreases in pore pressure and, consequently, the effect on the mud weight window. 1D and 4D robust geomechanical models are built based on all available data in a mature oil field. The 1D model was used to estimate all mechanical rock properties, stress, and pore pressure. The minimum and maximum horizontal stress were determined using the poroelastic horizontal strain model. Furthermore, the mechanical properties were calibrated using drained triaxial and uniaxial compression tests. The pore pressure was tested using modular dynamic tester log MDT. The Mohr–Coulomb model was applied in the 4D model to calculate the stress distribution in the depleted reservoir. According to study wells, the target area has been classified into four main groups in Mishrif reservoir based on depletion: highly, moderately, slightly, and no depleted region. Also, the results showed that the units had been classified into three main categories based on depletion state (from above to low depleted): L1.1, L1.2, and M1. The mean average reduction in minimum horizontal stress magnitude was 322 psi for L1.1, 183.86 psi for L1.2, and 115.56 psi for M1. Thus, the lower limit of fracture pressure dropped to a high value in L1.1, which is considered a weak point. As a result of changing horizontal stress, the mud weight window became narrow.

Comparative biomechanical analysis of rigid vs. flexible fixation devices for the lumbar spine: A geometrically patient-specific poroelastic finite element study

Background and ObjectiveLumbar spinal stenosis (LSS), or the narrowing of the spinal canal, continues to be the leading preoperative diagnosis for adults older than 65 years who undergo spine surgery. Although the treatment of LSS depends on its severity, the optimal surgical technique towards decreasing the risk of adjacent segment disease (ASD) remains elusive. This study aimed to comparatively analyze spinal biomechanics with rigid and flexible fixation devices (i.e., rigid and dynamic posterolateral fusion (PLF) and interspinous process (ISP) devices) during daily activities. MethodsUsing a validated parametric poroelastic finite element modeling approach, 8 subject-specific pre-operative models were developed, and their validity was evaluated. Parametric FE models of the lumbar spines were then regenerated based on post-operation images for (A) rigid PLF (B) dynamic PLF (C) rigid ISP device (Coflex) and (D) flexible ISP device (DIAM) at L4-L5 level. Biomechanical responses for instrumented and adjacent intervertebral discs (IVDs) were analyzed and compared subject to static and cyclic loading. ResultsThe preoperative models were well comparable with previous works in literature. The postoperative results for the PLF and Coflex rigid systems, demonstrated greater ROM; higher values of stress and strain in the AF region; and increased disc height and fluid loss at the adjacent levels, as compared with the pre-op models and the post-op results of the flexible systems (i.e., dynamic PLF and DIAM). The calculated forces on the facet joint were of smaller magnitude for the ISP devices as compared to the PLF, particularly during extension. ConclusionsThis study demonstrates that the dynamic PLF construct and DIAM implants could be effective to maintain the natural poroelastic characteristics of adjacent IVDs, which could be beneficial for enhancing long-term clinical outcomes. FEM provides clinicians with an invaluable patient-specific quantitative tool for informed surgical planning and discerning follow-up management.

The Human Meniscus Behaves as a Functionally Graded Fractional Porous Medium under Confined Compression Conditions

In this study, we observe that the poromechanical parameters in human meniscus vary spatially throughout the tissue. The response is anisotropic and the porosity is functionally graded. To draw these conclusions, we measured the anisotropic permeability and the “aggregate modulus” of the tissue, i.e., the stiffness of the material at equilibrium, after the interstitial fluid has ceased flowing. We estimated those parameters within the central portion of the meniscus in three directions (i.e., vertical, radial and circumferential) by fitting an enhanced model on stress relation confined compression tests. We noticed that a classical biphasic model was not sufficient to reproduce the observed experimental behaviour. We propose a poroelastic model based on the assumption that the fluid flow inside the human meniscus is described by a fractional porous medium equation analogous to Darcy’s law, which involves fractional operators. The fluid flux is then time-dependent for a constant applied pressure gradient (in contrast with the classical Darcy’s law, which describes a time independent fluid flux relation). We show that a fractional poroelastic model is well-suited to describe the flow within the meniscus and to identify the associated parameters (i.e., the order of the time derivative and the permeability). The results indicate that mean values of λβ,β in the central body are λβ=5.5443×10−10m4Ns1−β, β=0.0434, while, in the posterior and anterior regions, are λβ=2.851×10−10m4Ns1−β, β=0.0326 and λβ=1.2636×10−10m4Ns1−β, β=0.0232, respectively. Furthermore, numerical simulations show that the fluid flux diffusion is facilitated in the central part of the meniscus and hindered in the posterior and anterior regions.

Fluid induced deformation in porous media – Sensitivity analysis of a poroelastic model

The poromechanical behaviour of granular materials are influenced by the rheological properties and stress state of the injected fluid in addition to the state of the porous media. Fluid injection through a granular continuum generally results in the elastic or plastic deformation of the material which reflects as the change in porosity due to particle rearrangement. This phenomenon results in fluid induced instabilities in the porous media which is commonly observed during CO2 sequestration, oil and hydrocarbon recovery, high pressure grouting, hydraulic fracturing. Solid-fluid interaction in porous media is a fundamental multi-physics problem encountered in civil engineering, petroleum, and mining industries. Injecting fluid into the porous materials results in the deformation of the existing solid skeleton, especially when the flow rate is greater than the ability of the media to permeate the fluid. Deformations can be considered as poroelastic when the storage of reversible elastic energy controls the process. Among the existing poroelastic models which predict the fluid induced deformation in porous media, a coupled non-linear continuum-based model proposed by MacMinn is used in this study. This model takes into account various fluid properties and predicts the material response under different boundary conditions. The efficiency of this model to capture the deformation characteristics of porous media under different flow rate, fluid viscosity, porosity of the media and other geometric parameters will be carried out in this study.

Lactation alters fluid flow and solute transport in maternal skeleton: A multiscale modeling study on the effects of microstructural changes and loading frequency.

The female skeleton undergoes significant material and ultrastructural changes to meet high calcium demands during reproduction and lactation. Through the peri-lacunar/canalicular remodeling (PLR), osteocytes actively resorb surrounding matrix and enlarge their lacunae and canaliculi during lactation, which are quickly reversed after weaning. How these changes alter the physicochemical environment of osteocytes, the most abundant and primary mechanosensing cells in bone, are not well understood. In this study, we developed a multiscale poroelastic modeling technique to investigate lactation-induced changes in stress, fluid pressurization, fluid flow, and solute transport across multiple length scales (whole bone, porous midshaft cortex, lacunar-canalicular pore system (LCS), and pericellular matrix (PCM) around osteocytes) in murine tibiae subjected to axial compression at 3N peak load (~320με) at 0.5, 2, or 4Hz. Based on previously reported skeletal anatomical measurements from lactating and nulliparous mice, our models demonstrated that loading frequency, LCS porosity, and PCM density were major determinants of fluid and solute flows responsible for osteocyte mechanosensing, cell-cell signaling, and metabolism. When loaded at 0.5Hz, lactation-induced LCS expansion and potential PCM reduction promoted solute transport and osteocyte mechanosensing via primary cilia, but suppressed mechanosensing via fluid shear and/or drag force on the cell membrane. Interestingly, loading at 2 or 4Hz was found to overcome the mechanosensing deficits observed at 0.5Hz and these counter effects became more pronounced at 4Hz and with sparser PCM in the lactating bone. Synergistically, higher loading frequency (2, 4Hz) and sparser PCM enhanced flow-mediated mechanosensing and diffusion/convection of nutrients and signaling molecules for osteocytes. In summary, lactation-induced structural changes alter the local environment of osteocytes in ways that favor metabolism, mechanosensing, and post-weaning recovery of maternal bone. Thus, osteocytes play a role in balancing the metabolic and mechanical functions of female skeleton during reproduction and lactation.

Numerical simulation and microtomography study for drying a deformable isodiametric-cellular food

Abstract The aim of this work is the simulation of volumetric strain of tuberous crop during drying. We propose a poroelastic model for predicting the drying kinetics and volume loss of potato cubes during convective drying. The Biot’s theory of poroelasticity was used, which considers the Lamé parameters, Young’s modulus and Poisson’s ratio. Drying kinetics and volumetric strain were modeled and compared versus experimental data. An X-ray microtomograph coupled with image analysis was used to visualize the shape and size of the samples during drying. Drying experiments were conducted at 50, 60 and 70 °C, 20% RH, with an air velocity of 1 and 2 m/s. The drying process was interrupted several times to perform tomographic acquisitions. We found a period of ideal shrinkage, nevertheless, the volumetric strain reveals a kinetic behavior over time. The model computes the volumetric strain, which describes correctly the experimental data obtained by microtomography.

Poro-Viscoelastic Effects During Biomechanical Testing of Human Brain Tissue

Brain tissue is one of the softest tissues in the human body and the quantification of its mechanical properties has challenged scientists over the past decades. Associated experimental results in the literature have been contradictory as characterizing the mechanical response of brain tissue not only requires well-designed experimental setups that can record the ultrasoft response, but also appropriate approaches to analyze the corresponding data. Due to the extreme complexity of brain tissue behavior, nonlinear continuum mechanics has proven an expedient tool to analyze testing data and predict the mechanical response using a combination of hyper-, visco-, or poro-elastic models. Such models can not only allow for personalized predictions through finite element simulations, but also help to comprehensively understand the physical mechanisms underlying the tissue response. Here, we use a nonlinear poro-viscoelastic computational model to evaluate the effect of different intrinsic material properties (permeability, shear moduli, nonlinearity, viscosity) on the tissue response during different quasi-static biomechanical measurements, i.e., large-strain compression and tension as well as indentation experiments. We show that not only the permeability but also the properties of the viscoelastic solid largely control the fluid flow within and out of the sample. This reveals the close coupling between viscous and porous effects in brain tissue behavior. Strikingly, our simulations can explain why indentation experiments yield that white matter tissue in the human brain is stiffer than gray matter, while large-strain compression experiments show the opposite trend. These observations can be attributed to different experimental loading and boundary conditions as well as assumptions made during data analysis. The present study provides an important step to better understand experimental data previously published in the literature and can help to improve experimental setups and data analysis for biomechanical testing of brain tissue in the future.

Impact of adjacent pre-existing disc degeneration status on its biomechanics after single-level anterior cervical interbody fusion

Background and objectiveMechanics and biology may be interconnected and amplify each other during disc degeneration. It remains unknown the influence of pre-existing disc degeneration and its impact on adjacent segment degeneration (ASD) after anterior cervical discectomy and fusion (ACDF). This study aimed to discuss the necessity of including degenerated adjacent segments in single-level ACDF surgery from a biomechanical view. MethodsA poroelastic C2-T1 finite element model was created and validated. A C5-C6 ACDF model was developed based on this normal model. Moderate C4-C5 disc degeneration was created by appropriately modifying the morphology and tissue material properties in this fusion model. Degenerative morphology was modeled based on Thompson's grading system and Walraevens's scoring system for cervical spine, including disc height, whole disc area, nucleus pulposus (NP) area, endplate sclerosis and curvature. Stresses in disc and endplate and loads in facet joint were computed under moment loads in the fusion models with normal and pre-existing degenerative disc condition. ResultsAs for the adjacent disc, the stress values in degenerative condition were 7.41%, 5% and 5.26% larger than that in normal situation during extension, axial rotation and lateral bending motion, respectively. The disc stress changes mainly stemmed from annulus fibrosus (AF) tissue, but not NP. In the endplate, stress values of degeneration status were 4.17, 4.35 and 6.06% larger than that of normal condition under axial rotation, lateral bending and extension. The facet load magnitudes of pre-existing degeneration were 11.28, 11.57, 11.78 and 11.42% greater than that of normal condition in flexion, extension, axial rotation and lateral bending motion. ConclusionPre-existing degenerated disc experience increased biomechanical changes in adjacent segment after single-level ACDF. It may pose a long-term cumulative problem related to biomechanics in cervical spine after fusion. Before surgery, surgeons should be careful about selecting the fusion level.

Prediction of the natural frequencies of different degrees of degenerated human lumbar segments L2-L3 using dynamic finite element analysis

Background and ObjectiveChronic exposure to resonant environment may cause more serious injuries to human lumbar spine than other vibrations. On the condition that the natural frequency of human lumbar spine is known, excitation frequency from an external vibration source can be optimized to keep away from the natural frequency and thus avoid lumbar resonance. Therefore, this study aimed to present an approach to predict the natural frequency of the human lumbar spine. MethodsFour poroelastic finite element models of human L2–L3 spinal motion segments with different degrees of degeneration were established. Dynamic finite element analyses of these models during 1 h of vibration were then conducted. The mechanical parameters of these models under vibrations at different excitation frequencies were predicted. The excitation frequencies that resulted in the greatest changes in the lumbar mechanical parameters were identified as the natural frequencies of the established L2–L3 spinal motion segments. ResultsSimulation results showed that the natural frequencies of the healthy and mildly degenerated L2–L3 spinal motion segments, moderately degenerated L2–L3 spinal motion segments, and seriously degenerated L2–L3 spinal motion segments were in the range of 5–7, 3–5, and 1–3 Hz, respectively. ConclusionsThe predicted results indicated that the natural frequencies of the human L2–L3 spinal motion segments gradually decreased with the severity of degeneration. These phenomena may be related to changes in the lumbar structures and materials because of degeneration. This study provided a feasible method to predict the lumbar natural frequencies for different populations, which may be helpful in optimizing external vibration sources to avoid lumbar resonance.

Elastic, Dynamic Viscoelastic and Model-Derived Fibril-Reinforced Poroelastic Mechanical Properties of Normal and Osteoarthritic Human Femoral Condyle Cartilage

Osteoarthritis (OA) degrades articular cartilage and weakens its function. Modern fibril-reinforced poroelastic (FRPE) computational models can distinguish the mechanical properties of main cartilage constituents, namely collagen, proteoglycans, and fluid, thus, they can precisely characterize the complex mechanical behavior of the tissue. However, these properties are not known for human femoral condyle cartilage. Therefore, we aimed to characterize them from human subjects undergoing knee replacement and from deceased donors without known OA. Multi-step stress-relaxation measurements coupled with sample-specific finite element analyses were conducted to obtain the FRPE material properties. Samples were graded using OARSI scoring to determine the severity of histopathological cartilage degradation. The results suggest that alterations in the FRPE properties are not evident in the moderate stages of cartilage degradation (OARSI 2-3) as compared with normal tissue (OARSI 0-1). Drastic deterioration of the FRPE properties was observed in severely degraded cartilage (OARSI 4). We also found that the FRPE properties of femoral condyle cartilage related to the collagen network (initial fibril-network modulus) and proteoglycan matrix (non-fibrillar matrix modulus) were greater compared to tibial and patellar cartilage in OA. These findings may inform cartilage tissue-engineering efforts and help to improve the accuracy of cartilage representations in computational knee joint models.

Predicting acoustic performance of high surface area particle stacks with a poro-elastic model

Because of the high sound absorption they offer at low frequencies, there is a growing interest in high surface area particles and how they might be applied in noise control. Therefore, a model that can accurately predict the acoustic behavior of this type of materials will be useful in relevant applications. A poro-elastic model based on a combination of Biot theory and an existing rigid model of granular activated carbon (GAC) is introduced in the current work. The input parameters for this model consist of a certain number of properties that are known by measurement, and a set of values obtained by matching the model prediction with acoustic measurements. Measured absorption coefficients and surface impedance of stacks of several types of different activated carbon particles are shown in this paper. A fitting procedure that determines the unknown parameters is also described. It is shown that the model is able to predict the acoustic behavior of the particle stacks, and especially to capture the frame resonances at low frequencies, thus, validating the proposed model. Beyond the activated carbon used in the present tests, it is reasonable to generalize this model to stacks of other high surface area particles.

Numerical study of long-time growth of hydraulic fractures in a line drive

This work concerns detailed numerical modelling of injection fracturing in chalk formations under 2–3 km burial. The conditions and properties are chosen to resemble field conditions as typically found in oil fields subject to water injection in the Central Graben in the North Sea off-shore Denmark, Norway and the UK, but the studied phenomena are general to injection fracturing irrespective of geographical location and technological application. The complexity of reality is a challenging factor when mathematical models of hydraulic fractures in the subsurface are formulated. To this end the finite element method and a multi physics approach is instrumental, and data from fields developed and operated under conditions prudent for modelling are equally important. In the present study, a poroelastic finite element model is applied and used under plane strain conditions for a linear fracture mechanics investigation. This approach is suitable for analysing fracturing in arrays of parallel and horizontal wells as can be found in many oil fields undergoing water-flooding. Based on realistic field data, e.g. reservoir properties and rate and pressure scenarios representative for typical fields in the North Sea region, fracture initiation and fracture growth are analysed in details in a realistic field setting. Using a formulation of the J-integral, that includes Functionally Graded Material and poroelastic effects, it is demonstrated that formation fracture toughness determines fracture initiation and fracture height. The fracture propagation speed is controlled by continuity, e.g. the injection rate must balance the leak-off from the fracture wall. By combining the fields of geotechnical engineering, petroleum engineering and mechanical engineering with realistic data measurements, this study provides a novel realistic study of how the producer and injector pressure influence the fracturing process for induced hydraulic fractures in a line drive.

The influence of alkaline treatment on the mechanical and structural properties of bacterial cellulose

The unique mechanical properties of hydrated bacterial cellulose make it suitable for biomedical applications. This study evaluates the effect of concentrated sodium hydroxide treatment on the structural and mechanical properties of bacterial cellulose hydrogels using rheological, tensile, and compression tests combined with mathematical modelling. Bacterial cellulose hydrogels show a concentration-dependent and irreversible reduction in shear moduli, compression, and tensile strength after alkaline treatment. Applying a poroelastic biphasic model to through-thickness compressive stress-relaxation tests showed the alkaline treatment to induce no significant change in axial compression, an effect was observed in the radial direction, potentially due to the escape of water from within the hydrogel. Scanning electron microscopy showed a more porous structure of bacterial cellulose. These results show how concentration-dependent alkaline treatment induces selective weakening of intramolecular interactions between cellulose fibres, allowing the opportunity to precisely tune the mechanical properties for specific biomedical application, e.g., faster-degradable materials.

Poroelastic model in a vertically sealed gas storage: a case study from cyclic injection/production in a carbonate aquifer

SUMMARY Natural gas can be temporarily stored in a variety of underground facilities, such as depleted gas and oil fields, natural aquifers and caverns in salt rocks. Being extensively monitored during operations, these systems provide a favourable opportunity to investigate how pressure varies in time and space and possibly induces/triggers earthquakes on nearby faults. Elaborate and detailed numerical modelling techniques are often applied to study gas reservoirs. Here we show the possibilities and discuss the limitations of a flexible and easily formulated tool that can be straightforwardly applied to simulate temporal pore-pressure variations and study the relation with recorded microseismic events. We use the software POEL (POroELastic diffusion and deformation) which computes the poroelastic response to fluid injection/extraction in a horizontally layered poroelastic structure. We further develop its application to address the presence of vertical impermeable faults bounding the reservoir and of multiple injection/extraction sources. Exploiting available information on the reservoir geometry and physical parameters, and records of injection/extraction rates for a gas reservoir in southern Europe, we perform an extensive parametric study considering different model configurations. Comparing modelled spatiotemporal pore-pressure variations with in situ measurements, we show that the inclusion of vertical impermeable faults provides an improvement in reproducing the observations and results in pore-pressure accumulation near the faults and in a variation of the temporal pore-pressure diffusion pattern. To study the relation between gas storage activity and recorded local microseismicity, we applied different seismicity models based on the estimated pore-pressure distribution. This analysis helps to understand the spatial distribution of seismicity and its temporal modulation. The results show that the observed microseismicity could be partly linked to the storage activity, but the contribution of tectonic background seismicity cannot be excluded.

A numerical study on the effect of osmotic pressure on stress and strain in intercellular structures of tumor tissue in the poro-elastic model

A numerical study on the effect of osmotic pressure on stress and strain in intercellular structures of tumor tissue in the poro-elastic model

Runge-Kutta discontinuous Galerkin method for solving wave equations in 2D isotropic and anisotropic poroelastic media at low frequencies

We have developed a Runge-Kutta discontinuous Galerkin (RKDG) method for solving wave equations in isotropic and anisotropic poroelastic media at low frequencies. First, the 2D Biot’s two-phase equations are transformed into a first-order system with dissipation. Then, the system is discretized by using the discontinuous Galerkin method with a third-order Runge-Kutta time discretization. The numerical stability conditions for solving porous equations are also investigated. We test several examples to validate our method in isotropic and anisotropic poroelastic media. Comparisons of seismic responses with the finite-difference method on fine grids show the correctness of this method. Moreover, the numerical results indicate that the RKDG method can provide clear fast P-, slow P-, and S-waves for anisotropic poroelastic media on coarse meshes. Also, a two-layer porous model, a poroelastic-elastic model with horizontal interface, and an isotropic-anisotropic poroelastic model with sinusoidal interface demonstrate that our method can deal with complex wave propagation. Therefore, the simulation results show that RKDG is an accurate and stable method for solving Biot’s equations.

A poro-viscoelastic model for the subcutaneous injection of monoclonal antibodies

Subcutaneous delivery of therapeutic monoclonal antibodies (mAbs) has recently attracted exceptional interests in the pharmaceutical industry. The mechanical response of the tissue during the injection and the drug distribution after the injection are yet to be understood. We are motivated to study the subcutaneous injection using a nonlinear poro-elastic model with a viscoelastic constitutive law of the solid. The model is first validated against experimental data of relaxation tests. Then the poro-viscoelastic model is applied to subcutaneous injection, and is compared against a poro-hyperelastic model. The pressure at the injection site, which serves as an indicator of possible fracture initiation, differs significantly by taking into account the viscosity of the solid.

Impacts of Reservoir Water Level Fluctuation on Measuring Seasonal Seismic Travel Time Changes in the Binchuan Basin, Yunnan, China

An airgun source in a water reservoir has been developed in the past decade as a green active source that had been proven effective to derive short-term subsurface structural changes. However, seasonal water level fluctuation in the reservoir affects the airgun signal, and thus whether the airgun signals can be used to derive robust seasonal variation in subsurface structure remains unclear. We use the airgun data observed in the Binchuan basin to estimate the seasonal variation of seismic travel time and compare the results with those derived from ambient noise data in the same frequency band. Our main observation is that seasonal change δt/t from airgun is negatively correlated to the variation of dominant frequency and water table fluctuation in the reservoir. One possible explanation is that water table fluctuation in the reservoir affects the dominant frequency of the airgun signal and causes significant phase shift. We also compute the travel time changes in P-wave from the empirical Green’s function after deconvolving the waveforms from a reference station that is 50 m from the airgun source. The dominant frequency after deconvolution still shows seasonal variation and correlates inversely to the travel time changes, suggesting that deconvolution cannot completely eliminate the source effect on travel time changes. We also use ambient noise cross-correlation to retrieve coda waves and then derive travel time changes in monthly stacked cross-correlations relative to a yearly average cross-correlation. We observe that seismic travel time increases to its local maximum in the end of August. The travel time changes lag behind the precipitation for about one month. We apply a poroelastic physical model to explain seismic travel time changes and find that a combined effect from precipitation and evaporation might induce the seasonal changes as shown in the ambient noise data. However, the pattern of travel time changes from the airgun differs from that from ambient noise, reflecting the strong effects of airgun source property changes. Therefore, we should be cautious to derive long-term subsurface structural variation from the airgun source and put more attention on stabilizing the dominant frequency of each excitation in the future experiments.

A dual porosity poroelastic model for simulation of gas flow in saturated claystone as a potential host rock for deep geological repositories

Claystone is considered as a potential host rock for radioactive waste disposal at great depth in many countries. Considerable quantity of gas can be generated in a deep geological repository (DGR) due to several processes, which may affect the integrity of the host rock and the safety of the DGR. Thus, understanding the migration of gas through the host rock is an essential requirement in developing a safety case for a deep geological disposal of radioactive waste. In this paper, a fully coupled dual porosity poroelastic model is developed to simulate the gas transport process in initially saturated claystone. The model considers the hydro-mechanical behavior for the fractured porous medium consisting of both the porous matrix (represented by the porous continuum) and fractures (represented by the fractured continuum), which are two separate and overlapping porous media. Each continuum has its own constitutive law that imposes a superimposed effect on the behavior of the fractured porous medium. The mechanical coupling between the two continua is through total stress equilibrium and strain superimposition technique, while the fluid coupling in the two continua is through water exchange term and porosity interaction. The capability of the developed dual porosity poroelastic model is evaluated by comparing the simulated results with that recorded in the laboratory gas injection tests on claystone, in which the main experimental behaviors, i.e., the major gas breakthrough, sample volume dilation, gas induced fracturing are well represented.

Lateral fluid propagation and strike slip fault reactivation related to hydraulic fracturing and induced seismicity in the Duvernay Formation, Fox Creek area, Alberta

SUMMARY Fault shear slip potential is analysed in the area where induced earthquakes (up to 3.9 Mw) occurred in May–June 2015 approximately 30 km south of Fox Creek, Western Canada Sedimentary Basin, Canada. The induced earthquakes were generated by the hydraulic fracturing of the Upper Devonian Duvernay Formation. Interpretation of a 3-D seismic survey and analysis of the ant tracking attribute identifies a linear discontinuity that likely represents a subvertical fault with strike length of 1.4 km, which is aligned with the zone of induced earthquake hypocentres. 1-D–3-D geomechanical modelling is conducted to characterize mechanical rock properties, initial reservoir pressure and stress field. Hydraulic fracture propagation and reservoir pressure buildup simulations are run to analyse lateral fluid pressure diffusion during well treatment. The interaction of natural fractures introduced as Discrete Fracture Network and hydraulic fractures is tested. 3-D poroelastic reservoir geomechanical modelling is completed to simulate slip reactivation of the identified fault zone. The obtained results support that additional pressure buildup of 20 MPa in treatment wells can propagate laterally along hydraulic fractures (and potentially natural fracture network) for about 550 m and reach the fault zone. The increase of fluid pressure by 20 MPa in the fault zone results in dextral slip along the fault, mostly in the interval of the Duvernay and overlying Ireton Formations, corroborating prior focal mechanism results and hypocentral depths. The simulations indicate that lateral transmission of additional fluid pressure from the fracturing stimulation area to the fault zone could happen in a few days after the treatment of lateral wells that is supported by the observed induced earthquakes. This study helps to quantify changes in fluid pressure and stresses that may result in fault shear slip during hydraulic fracturing and predict the potential of induced seismicity connected to hydrocarbon production from the Duvernay Play.