Hydraulic fracturing then fluid circulation in enhanced geothermal system(EGS)reservoirs have been shown to induce seismicity remote from the stimulation-potentially generated by the distal projection of thermoporoela...Hydraulic fracturing then fluid circulation in enhanced geothermal system(EGS)reservoirs have been shown to induce seismicity remote from the stimulation-potentially generated by the distal projection of thermoporoelastic stresses.We explore this phenomenon by evaluating stress perturbations resulting from stimulation of a single stage of hydraulic fracturing that is followed by thermal depletion of a prismatic zone adjacent to the hydraulic fracture.We use Coulomb failure stress to assess the effect of resulting stress perturbations on instability on adjacent critically-stressed faults.Results show that hydraulic fracturing in a single stage is capable of creating stress perturbations at distances to 1000 m that reach 10^(-5)-10^(-4)MPa.At a closer distance,the magnitude of stress perturbations increases even further.The stress perturbation induced by temperature depletion could also reach 10^(-3)-10^(-2)MPa within 1000 m-much higher than that by hydraulic fracturing.Considering that a critical change in Coulomb failure stress for fault instability is 10^(-2)MPa,a single stage of hydraulic fracturing and thermal drawdown are capable of reactivating critically-stressed faults at distances within 200 m and 1000 m,respectively.These results have important implications for understanding the distribution and magnitudes of stress perturbations driven by thermoporoelastic effects and the associated seismicity during the simulation and early production of EGS reservoirs.展开更多
Mining-related seismicity poses significant challenges in underground coal mining due to its complex rupture mechanisms and associated hazards.To bridge gaps in understanding these intricate processes,this study emplo...Mining-related seismicity poses significant challenges in underground coal mining due to its complex rupture mechanisms and associated hazards.To bridge gaps in understanding these intricate processes,this study employed a multi-local seismic monitoring network,integrating both in-mine and local instruments at overlapping length scales.We specifically focused on a damaging local magnitude(ML)2.6 event and its aftershocks that occurred on 10 September 2022 in the vicinity of the 3308 working face of the Yangcheng coal mine in Shandong Province,China.Moment tensor(MT)inversion revealed a complex cascading rupture mechanism:an initial moment magnitude(M_(w))2.2 normal fault slip along the DF60 fault in an ESEeWNW direction,transitioning to a M_(w)3.0 event as the FD24 and DF60 faults unclamped.The scale-independent self-similarity and stress heterogeneity of mining-related seismicity were investigated through source parameter calculations,providing valuable insights into the driving mechanism of these seismic sequences.The in-mine network,constrained by its low dynamic changes,captured only the nucleation phase of the DF60 fault.Furthermore,standard decomposition of the MT solution from the seismic network proved inadequate for accurately identifying the complex nature of the rupture.To enhance safety and risk management in mining environments,we examined the implications of source reactivation within the cluster area post-stress-adjustment.This comprehensive multiscale analysis offers crucial insights into the complex rupture mechanisms and hazards associated with mining-related seismicity.The results underscore the importance of continuous multi-local network monitoring and advanced analytical techniques for improved disaster assessment and risk mitigation in mining operations.展开更多
Triggered seismicity is a key hazard where fluids are injected or withdrawn from the subsurface and may impact permeability. Understanding the mechanisms that control fluid injection-triggered seismicity allows its mi...Triggered seismicity is a key hazard where fluids are injected or withdrawn from the subsurface and may impact permeability. Understanding the mechanisms that control fluid injection-triggered seismicity allows its mitigation. Key controls on seismicity are defined in terms of fault and fracture strength, second-order frictional response and stability, and competing fluid-driven mechanisms for arrest. We desire to constrain maximum event magnitudes in triggered earthquakes by relating pre-existing critical stresses to fluid injection volume to explain why some recorded events are significantly larger than anticipated seismic moment thresholds. This formalism is consistent with several uncharacteristically large fluid injection-triggered earthquakes. Such methods of reactivating fractures and faults by hydraulic stimulation in shear or tensile fracturing are routinely used to create permeability in the subsurface. Microearthquakes (MEQs) generated by such stimulations can be used to diagnose permeability evolution. Although high-fidelity data sets are scarce, the EGS-Collab and Utah FORGE hydraulic stimulation field demonstration projects provide high-fidelity data sets that concurrently track permeability evolution and triggered seismicity. Machine learning deciphers the principal features of MEQs and the resulting permeability evolution that best track permeability changes – with transfer learning methods allowing robust predictions across multiple eological settings. Changes in permeability at reactivated fractures in both shear and extensional modes suggest that permeability change (Δk) scales with the seismic moment (M) of individual MEQs as Δk∝M. This scaling relation is exact at early times but degrades with successive MEQs, but provides a method for characterizing crustal permeability evolution using MEQs, alone. Importantly, we quantify for the first time the role of prestress in defining the elevated magnitude and seismic moment of fluid injection-triggered events, and demonstrate that such MEQs can also be used as diagnostic in quantifying permeability evolution in the crust.展开更多
Understanding the relationship between normal stiffness and permeability in rock fractures under high and true-triaxial in situ stress conditions is critical to assess hydro-mechanical coupling in the Earth's crus...Understanding the relationship between normal stiffness and permeability in rock fractures under high and true-triaxial in situ stress conditions is critical to assess hydro-mechanical coupling in the Earth's crust.Previous data on stiffness–permeability relations are measured under uniaxial stress states as well as under normal stress.However,many projects involve faulted formations with complex three-dimensional(3D)stress states or significant changes to the original stress state.We rectified this by following the permeability evolution using a true-triaxial stress-permeability apparatus as well as independently applying a spectrum of triaxial stresses from low to high.The relationship between permeability and fracture normal stiffness was quantified using constraints based on the principle of virtual work.The impacts of fracture-lateral and fracture-normal stresses on permeability and normal stiffness evolution were measured.It was found that permeability decreases with increasing fracture-lateral and fracture-normal stresses as a result of Poisson confinement,independent of the orientation of the fracture relative to the stresses.The lateral stresses dominated the evolution of normal stiffness at lower normal stresses(σ_(3)=10 MPa)and played a supplementary role at higher normal stresses(σ_(3)>10 MPa).Moreover,correlations between the evolution of permeability and normal stiffness were extended beyond the low-stiffness,high-permeability region to the high-stiffness,low-permeability region under high fracture-lateral stresses(10–80 MPa)with fracture-normal stress(10–50 MPa)conditions.Again,high lateral stresses further confined the fracture and therefore reduced permeability and increased normal stiffness,which exceeded the previous reported stiffness under no lateral stress conditions.This process enabled us to identify a fundamental change in the flow regime from multi-channel to isolated channelized flow.These results provide important characterizations of fracture permeability in the deep crust,including recovery from deep shale-gas reservoirs.展开更多
Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristic...Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristics complicate the gas-water two-phase flow process in porous media following hydrate decomposition, posing challenges for efficient development. This study examines the transport response of clayey-silt reservoir samples from the Shenhu area using gas-water two-phase flow experiments and CT scanning to explore changes in pore structure, gas-water distribution, and relative permeability under varying flow conditions. The results indicate that pore heterogeneity significantly influences flow characteristics. Gas preferentially displaces water in larger pores, forming fracture-like pores, which serve as preferential flow channels for gas migration. The preferential flow channels enhance gas-phase permeability up to 19 times that of the water phase when fluid pressures exceed total stresses. However,small pores retain liquid, leading to a high residual water saturation of 0.561. CT imaging reveals that these hydro-fractures improve gas permeability but also confine gas flow to specific channels. Pore network analysis shows that gas injection expands the pore-throat network, enhancing connectivity and forming fracture-like pores. Residual water remains trapped in smaller pores and throats, while structural changes, including new fractures, improve gas flow pathways and overall connectivity. Relative permeability curves demonstrate a narrow gas-water cocurrent-flow zone, a right-shifted iso-permeability point and high reservoir capillary pressure, indicating a strong "water-blocking" effect. The findings suggest that optimizing reservoir stimulation techniques to enhance fracture formation, reduce residual water saturation, and improve gas flow capacity is critical for efficient hydrate reservoir development.展开更多
In this study,we employed Bayesian inversion coupled with the summation-by-parts and simultaneousapproximation-term(SBP-SAT)forward simulation method to elucidate the mechanisms behind mininginduced seismic events cau...In this study,we employed Bayesian inversion coupled with the summation-by-parts and simultaneousapproximation-term(SBP-SAT)forward simulation method to elucidate the mechanisms behind mininginduced seismic events caused by fault slip and their potential effects on rockbursts.Through Bayesian inversion,it is determined that the sources near fault FQ14 have a significant shear component.Additionally,we analyzed the stress and displacement fields of high-energy events,along with the hypocenter distribution of aftershocks,which aided in identifying the slip direction of the critically stressed fault FQ14.We also performed forward modeling to capture the complex dynamics of fault slip under varying friction laws and shear fracture modes.The selection of specific friction laws for fault slip models was based on their ability to accurately replicate observed slip behavior under various external loading conditions,thereby enhancing the applicability of our findings.Our results suggest that the slip behavior of fault FQ14 can be effectively understood by comparing different scenarios.展开更多
In-situ upgrading by heating is feasible for low-maturity shale oil,where the pore space dynamically evolves.We characterize this response for a heated substrate concurrently imaged by SEM.We systematically follow the...In-situ upgrading by heating is feasible for low-maturity shale oil,where the pore space dynamically evolves.We characterize this response for a heated substrate concurrently imaged by SEM.We systematically follow the evolution of pore quantity,size(length,width and cross-sectional area),orientation,shape(aspect ratio,roundness and solidity)and their anisotropy—interpreted by machine learning.Results indicate that heating generates new pores in both organic matter and inorganic minerals.However,the newly formed pores are smaller than the original pores and thus reduce average lengths and widths of the bedding-parallel pore system.Conversely,the average pore lengths and widths are increased in the bedding-perpendicular direction.Besides,heating increases the cross-sectional area of pores in low-maturity oil shales,where this growth tendency fluctuates at<300℃ but becomes steady at>300℃.In addition,the orientation and shape of the newly-formed heating-induced pores follow the habit of the original pores and follow the initial probability distributions of pore orientation and shape.Herein,limited anisotropy is detected in pore direction and shape,indicating similar modes of evolution both bedding-parallel and bedding-normal.We propose a straightforward but robust model to describe evolution of pore system in low-maturity oil shales during heating.展开更多
The injection of large volumes of natural gas into geological formations,as is required for underground gas storage,leads to alterations in the effective stress exerted on adjacent faults.This increases the potential ...The injection of large volumes of natural gas into geological formations,as is required for underground gas storage,leads to alterations in the effective stress exerted on adjacent faults.This increases the potential for their reactivation and subsequent earthquake triggering.Most measurements of the frictional properties of rock fractures have been conducted under normal and shear stresses.However,faults in gas storage facilities exist within a true three-dimensional(3D)stress state.A double-direct shear experiment on rock fractures under both lateral and normal stresses was conducted using a true triaxial loading system.It was observed that the friction coefficient increases with increasing lateral stress,but decreases with increasing normal stress.The impact of lateral and normal stresses on the response is primarily mediated through their influence on the initial friction coefficient.This allows for an empirical modification of the rate-state friction model that considers the influence of lateral and normal stresses.The impact of lateral and normal stresses on observed friction coefficients is related to the propensity for the production of wear products on the fracture surfaces.Lateral stresses enhance the shear strength of rock(e.g.Mogi criterion).This reduces asperity breakage and the generation of wear products,and consequently augments the friction coefficient of the surface.Conversely,increased normal stresses inhibit dilatancy on the fracture surface,increasing the breakage of asperities and the concomitant production of wear products that promote rolling deformation.This ultimately reduces the friction coefficient.展开更多
Basalt is a major component of the earth and moon crust.Mineral composition and temperature influence frictional instability and thus the potential for seismicity on basaltic faults.We performed velocitystepping shear...Basalt is a major component of the earth and moon crust.Mineral composition and temperature influence frictional instability and thus the potential for seismicity on basaltic faults.We performed velocitystepping shear experiments on basalt gouges at a confining pressure of 100 MPa,temperatures in the range of 100-400℃ and with varied obsidian mass fractions of 0-100%under wet/dry conditions to investigate the frictional strength and stability of basaltic faults.We observe a transition from velocity-neutral to velocity-weakening behaviors with increasing obsidian content.The frictional stability response of the mixed obsidian/basalt gouges is characterized by a transition from velocitystrengthening to velocity-weakening at 200℃ and another transition to velocity-strengthening at temperatures>300℃.Conversely,frictional strengths of the obsidian-bearing gouges are insensitive to temperature and wet/dry conditions.These results suggest that obsidian content dominates the potential seismic response of basaltic faults with the effect of temperature controlling the range of seismogenic depths.Thus,shallow moonquakes tend to occur in the lower lunar crust due to the corresponding anticipated higher glass content and a projected temperature range conducive to velocity-weakening behavior.These observations contribute to a better understanding of the nucleation mechanism of shallow seismicity in basaltic faults.展开更多
In-situ thermal upgrading is used to tune the pore system in low-maturity oil shales. We introduce fractal dimension(D), form factor(ff) and stochastic entropy(H) to quantify the heating-induced evolution of pore morp...In-situ thermal upgrading is used to tune the pore system in low-maturity oil shales. We introduce fractal dimension(D), form factor(ff) and stochastic entropy(H) to quantify the heating-induced evolution of pore morphological complexity and azimuthal disorder and develop a model to estimate the impact on seepage capacity via permeability. Experiments are conducted under recreated in-situ temperatures and consider anisotropic properties—both parallel and perpendicular to bedding. Results indicate that azimuthal distribution of pores in the bedding-parallel direction are dispersed, while those in the bedding-perpendicular direction are concentrated. D values indicate that higher temperatures reduce the uniformity of the pore size distribution(PSD) in the bedding-parallel direction but narrow the PSD in the bedding-perpendicular direction. The greater ff(> 0.7) values in the bedding-parallel direction account for a large proportion, while the dominated in the bedding-perpendicular direction locates within 0.2-0.7, for all temperatures. The H value of the bedding-parallel sample remains stable at ~0.925 during heating, but gradually increases from 0.808 at 25℃ to 0.879 at 500℃ for the beddingperpendicular sample. Congruent with a mechanistic model, the permeability at 500℃ is elevated~1.83 times(bedding-parallel) and ~6.08 times(bedding-perpendicular) relative to that at 25℃—confirming the effectiveness of thermal treatment in potentially enhancing production from low-maturity oil shales.展开更多
In current dual porosity/permeability models,there exists a fundamental assumption that the adsorption-induced swelling is distributed uniformly within the representative elementary volume (REV),irrespective of its in...In current dual porosity/permeability models,there exists a fundamental assumption that the adsorption-induced swelling is distributed uniformly within the representative elementary volume (REV),irrespective of its internal structures and transient processes.However,both internal structures and transient processes can lead to the non-uniform swelling.In this study,we hypothesize that the non-uniform swelling is responsible for why coal permeability in experimental measurements is not only controlled by the effective stress but also is affected by the adsorption-induced swelling.We propose a concept of the swelling triangle composed of swelling paths to characterize the evolution of the non-uniform swelling and serve as a core link in coupled multiphysics.A swelling path is determined by a dimensionless volumetric ratio and a dimensionless swelling ratio.Different swelling paths have the same start and end point,and each swelling path represents a unique swelling case.The swelling path as the diagonal of the triangle represents the case of the uniform swelling while that as the two perpendicular boundaries represents the case of the localized swelling.The paths of all intermediate cases populate inside the triangle.The corresponding relations between the swelling path and the response of coal multiphysics are established by a non-uniform swelling coefficient.We define this method as the triangle approach and corresponding models as swelling path-based ones.The proposed concept and models are verified against a long-term experimental measurement of permeability and strains under constant effective stress.Our results demonstrate that during gas injection,coal multiphysics responses have a close dependence on the swelling path,and that in both future experiments and field predictions,this dependence must be considered.展开更多
Coal measure source rocks, located in the Xihu Sag of the East China Sea Shelf Basin, were analyzed to define the hydrocarbon generation potential, organic geochemistry/petrology characteristics, and coal preservation...Coal measure source rocks, located in the Xihu Sag of the East China Sea Shelf Basin, were analyzed to define the hydrocarbon generation potential, organic geochemistry/petrology characteristics, and coal preservation conditions. The Pinghu source rocks in the Xihu Sag are mainly gas-prone accompany with condensate oil generation. The coals and shales of the Pinghu Formation are classified from "fair" to "excellent" source rocks with total organic carbon(TOC) contents ranging from 25.2% to 77.2% and 1.29% to 20.9%, respectively. The coals are richer in TOC and S1+S2 than the shales, indicating that the coals have more generation potential per unit mass. Moreover, the kerogen type of the organic matter consists of types Ⅱ-Ⅲ and Ⅲ, which the maturity Ro ranges from 0.59% to 0.83%. Petrographically, the coals and shales are dominated by vitrinite macerals(69.1%–96.8%) with minor proportions of liptinite(2.5%–17.55%) and inertinite(0.2%–6.2%). The correlation between maceral composition and S1+S2 indicates that the main contributor to the generation potential is vitrinite. Therefore, the coals and shales of the Pinghu Formation has good hydrocarbon generation potential, which provided a good foundation for coal measure gas accumulation. Furthermore, coal facies models indicates that the Pinghu coal was deposited in limno-telmatic environment under high water levels, with low tree density(mainly herbaceous) and with low-moderate nutrient supply. Fluctuating water levels and intermittent flooding during the deposition of peat resulted in the inter-layering of coal, shale and sandstone, which potentially providing favorable preservation conditions for coal measure gas.展开更多
Coal permeability measurements are normally conducted under the assumption that gas pressure in the matrix is equalized with that in fracture and that gas sorption-induced swelling/shrinking strain is uniformly distri...Coal permeability measurements are normally conducted under the assumption that gas pressure in the matrix is equalized with that in fracture and that gas sorption-induced swelling/shrinking strain is uniformly distributed within the coal.However,the validity of this assumption has long been questioned and differential strain between the fracture strain and the bulk strain has long been considered as the primary reason for the inconsistency between experimental data and poroelasticity solutions.Although efforts have been made to incorporate the impact into coal permeability models,the fundamental nature of those efforts to split the matrix strain between fracture and coal bulk remains questionable.In this study,a new concept of differential swelling index(DSI)was derived to theoretically define the relation among sorption-induced strains of the coal bulk,fracture,and coal matrix at the equilibrium state.DSI was a function of the equilibrium pressure and its magnitudes were regulated by the Langmuir constants of both the matrix and the coal bulk.Furthermore,a spectrum of DSI-based coal permeability models was developed to explicitly consider the effect of differential strains.These models were verified with the experimental data under the conditions of uniaxial strain,constant confining pressure,and constant effective stress.展开更多
Heat production from geothermal reservoirs is a typical heat transfer process involving a cold working fluid contacting a hot rock formation.Compared to the thermal-physical characteristics of water,supercritical CO_(...Heat production from geothermal reservoirs is a typical heat transfer process involving a cold working fluid contacting a hot rock formation.Compared to the thermal-physical characteristics of water,supercritical CO_(2)(scCO_(2))has a higher heat storage capacity over a wide temperature-pressure range and may be favored as a heat transfer fluid.Singularly characteristic of scCO_(2)-based heat extraction is that the hydraulic-thermal properties of the scCO_(2) vary dramatically and dynamically with the spatial pressure gradient during unsteady-state flow along fracture.This highly nonlinear behavior presents a challenge in the accurate estimation of heat extraction efficiency in scCO_(2)-based EGS.In this paper,a thermal-h ydraulic-mechanical(THM)coupled model is developed by considering deformation of the fractured reservoir,non-Darcy flow and the varying thermal-physical properties of scCO_(2).The proposed model is validated by matching the modeling temperature distribution with published data.The results show that during continuous injection of scCO_(2),the fracture first widens and then narrows,ultimately reopening over the long term.The sequential fracture deformation behaviors are in response to the combined impacts of mechanical compression and thermally-induced deformation.By controlling the injection parameters of the scCO_(2),it is found that the heat extraction rate is positively correlated to its pore pressure or mass flow rate.The heat extraction rate can be significantly enhanced,when the inlet temperature of scCO_(2) is below its critical temperature.As a result,the heat increment recovered per unit mass of scCO_(2) decreases as the hot rock is gradually cooled.Meanwhile,the heat increment recovered per unit mass of scCO_(2) decreases by increasing the inlet temperature of scCO_(2) or its mass flow rate,but increases as the outlet pressure rises.Furthermore,multi-linear regression indicates that controlling the inlet temperature of the scCO_(2) can significantly improve the thermodynamic efficiency of heat extraction.展开更多
Subsurface fluid injections can disturb the effective stress regime by elevating pore pressure and potentially reactivate faults and fractures.Laboratory studies indicate that fracture rheology and permeability in suc...Subsurface fluid injections can disturb the effective stress regime by elevating pore pressure and potentially reactivate faults and fractures.Laboratory studies indicate that fracture rheology and permeability in such reactivation events are linked to the roughness of the fracture surfaces.In this study,we construct numerical models using discrete element method(DEM)to explore the influence of fracture surface roughness on the shear strength,slip stability,and permeability evolution during such slip events.For each simulation,a pair of analog rock coupons(three-dimensional bonded quartz particle analogs)representing a mated fracture is sheared under a velocity-stepping scheme.The roughness of the fracture is defined in terms of asperity height and asperity wavelength.Results show that(1)Samples with larger asperity heights(rougher),when sheared,exhibit a higher peak strength which quickly devolves to a residual strength after reaching a threshold shear displacement;(2)These rougher samples also exhibit greater slip stability due to a high degree of asperity wear and resultant production of wear products;(3)Long-term suppression of permeability is observed with rougher fractures,possibly due to the removal of asperities and redistribution of wear products,which locally reduces porosity in the dilating fracture;and(4)Increasing shear-parallel asperity wavelength reduces magnitudes of stress drops after peak strength and enhances fracture permeability,while increasing shear-perpendicular asperity wavelength results in sequential stress drops and a delay in permeability enhancement.This study provides insights into understanding of the mechanisms of frictional and rheological evolution of rough fractures anticipated during reactivation events.展开更多
Optimization of fracturing perforation is of great importance to the commingling gas production in coal measure strata.In this paper,a 3 D lattice algorithm hydraulic fracturing simulator was employed to study the eff...Optimization of fracturing perforation is of great importance to the commingling gas production in coal measure strata.In this paper,a 3 D lattice algorithm hydraulic fracturing simulator was employed to study the effects of perforation position and length on hydraulic fracture propagation in coal measures of the Lin-Xing block,China.Based on field data,three lithologic combinations are simulated:1)a thick section of coal seam sandwiched by sandstones;2)a thin coal seam layer overlay by gas-bearing tight sandstone;3)two coal seams separated by a thin layer of sandstone.Our simulation shows that perforation position and length in multi-layer reservoirs play a major role in hydraulic fracture propagation.Achieving maximum stimulated volume requires consideration of lithologic sequence,coal seam thickness,stress states,and rock properties.To improve the combined gas production in coal measure strata,it is possible to simultaneously stimulate multiple coal seams or adjacent gas-bearing sandstones.In these cases,perforation location and length also significantly impact fracture propagation,and therefore should be carefully designed.Our simulation results using 3 D lattice algorithm are qualitatively consistent with laboratory physical simulation.3 D lattice models can be used to effectively simulate the fracture propagation through layers in coal measure strata.The numerical results provide guidance for perforation optimization in the hydraulic fracturing of coal measure strata.展开更多
Multi-seam mining often leads to the retention of a significant number of coal pillars for purposes such as protection,safety,or water isolation.However,stress concentration beneath these residual coal pillars can sig...Multi-seam mining often leads to the retention of a significant number of coal pillars for purposes such as protection,safety,or water isolation.However,stress concentration beneath these residual coal pillars can significantly impact their strength and stability when mining below them,potentially leading to hydraulic support failure,surface subsidence,and rock bursting.To address this issue,the linkage between the failure and instability of residual coal pillars and rock strata during multi-seam mining is examined in this study.Key controls include residual pillar spalling,safety factor(f.),local mine stiffness(LMS),and the post-peak stiffness(k)of the residual coal pillar.Limits separating the two forms of failure,progressive versus dynamic,are defined.Progressive failure results at lower stresses when the coal pillar transitions from indefinitely stable(f,>1.5)to failing(f,<1.5)when the coal pillar can no longer remain stable for an extended duration,whereas sud-den(unstable)failure results when the strength of the pillar is further degraded and fails.The transition in mode of failure is defined by the LMS/k ratio.Failure transitions from quiescent to dynamic as LMS/k.<1,which can cause chain pillar instability propagating throughout the mine.This study provides theoretical guidance to define this limit to instability of residual coal pillars for multi-seam mining in similar mines.展开更多
Stimulated shale reservoirs consist of kerogen,inorganic matter,secondary and hydraulic fractures.The dispersed distribution of kerogen within matrices and complex gas flow mechanisms make production evaluation challe...Stimulated shale reservoirs consist of kerogen,inorganic matter,secondary and hydraulic fractures.The dispersed distribution of kerogen within matrices and complex gas flow mechanisms make production evaluation challenging.Here we establish an analytical method that addresses kerogen-inorganic matter gas transfer,dispersed kerogen distribution,and complex gas flow mechanisms to facilitate evaluating gas production.The matrix element is defined as a kerogen core with an exterior inorganic sphere.Unlike most previous models,we merely use boundary conditions to describe kerogen-inorganic matter gas transfer without the instantaneous kerogen gas source term.It is closer to real inter-porosity flow conditions between kerogen and inorganic matter.Knudsen diffusion,surface diffusion,adsorption/desorption,and slip corrected flow are involved in matrix gas flow.Matrix-fracture coupling is realized by using a seven-region linear flow model.The model is verified against a published model and field data.Results reveal that inorganic matrices serve as a major gas source especially at early times.Kerogen provides limited contributions to production even under a pseudo-steady state.Kerogen properties’influence starts from the late matrix-fracture inter-porosity flow regime,while inorganic matter properties control almost all flow regimes except the early-mid time fracture linear flow regime.The contribution of different linear flow regions is also documented.展开更多
This paper experimentally explores the frictional sliding behavior of two simulated gouges:one,a series of quartz–smectite mixtures,and the other,powdered natural rocks,aiming to evaluate and codify the effect of min...This paper experimentally explores the frictional sliding behavior of two simulated gouges:one,a series of quartz–smectite mixtures,and the other,powdered natural rocks,aiming to evaluate and codify the effect of mineralogy on gouge dilation and frictional strength,stability,and healing.Specifically,velocity-stepping and slide-hold-slide experiments were performed in a double direct shear configuration to analyze frictional constitutive parameters at room temperature,under normal stresses of 10,20,and 40 MPa.Gouge dilation was measured based on the applied step-wise changes in shear velocity.The frictional response of the quartz–smectite mixtures and powdered natural rocks are affected by their phyllosilicate content.Frictional strength and healing rates decrease with increasing phyllosilicate content,and at 20 wt.%a transition from velocity-weakening to velocity-strengthening behavior was noted.For both suites of gouges,dilation is positively correlated with frictional strength and healing rates,and negatively correlated with frictional stability.Changes in the permeability of gouge-filled faults were estimated from changes in mean porosity,indexed through measured magnitudes of gouge dilation.This combined analysis implies that the reactivation of caprock faults filled with phyllosilicaterich gouges may have a strong influence on permeability evolution in caprock faults.展开更多
基金funded by the National Natural Science Foundation of China(Grant Nos.42107163 and 42320104003)support from the G.Albert Shoemaker endowment.
文摘Hydraulic fracturing then fluid circulation in enhanced geothermal system(EGS)reservoirs have been shown to induce seismicity remote from the stimulation-potentially generated by the distal projection of thermoporoelastic stresses.We explore this phenomenon by evaluating stress perturbations resulting from stimulation of a single stage of hydraulic fracturing that is followed by thermal depletion of a prismatic zone adjacent to the hydraulic fracture.We use Coulomb failure stress to assess the effect of resulting stress perturbations on instability on adjacent critically-stressed faults.Results show that hydraulic fracturing in a single stage is capable of creating stress perturbations at distances to 1000 m that reach 10^(-5)-10^(-4)MPa.At a closer distance,the magnitude of stress perturbations increases even further.The stress perturbation induced by temperature depletion could also reach 10^(-3)-10^(-2)MPa within 1000 m-much higher than that by hydraulic fracturing.Considering that a critical change in Coulomb failure stress for fault instability is 10^(-2)MPa,a single stage of hydraulic fracturing and thermal drawdown are capable of reactivating critically-stressed faults at distances within 200 m and 1000 m,respectively.These results have important implications for understanding the distribution and magnitudes of stress perturbations driven by thermoporoelastic effects and the associated seismicity during the simulation and early production of EGS reservoirs.
基金funded by the National Natural Science Foundation of China(Grant No.51574225)Shandong Energy Group(Grant No.SNKJ2022BJ03-R28)for Caiping Lu+1 种基金the Research Team on MonitoringActivity Mechanisms of Unnatural Earthquakes of Shandong Earthquake Agency(Grant No.TD202301)for Chengyu Liu.
文摘Mining-related seismicity poses significant challenges in underground coal mining due to its complex rupture mechanisms and associated hazards.To bridge gaps in understanding these intricate processes,this study employed a multi-local seismic monitoring network,integrating both in-mine and local instruments at overlapping length scales.We specifically focused on a damaging local magnitude(ML)2.6 event and its aftershocks that occurred on 10 September 2022 in the vicinity of the 3308 working face of the Yangcheng coal mine in Shandong Province,China.Moment tensor(MT)inversion revealed a complex cascading rupture mechanism:an initial moment magnitude(M_(w))2.2 normal fault slip along the DF60 fault in an ESEeWNW direction,transitioning to a M_(w)3.0 event as the FD24 and DF60 faults unclamped.The scale-independent self-similarity and stress heterogeneity of mining-related seismicity were investigated through source parameter calculations,providing valuable insights into the driving mechanism of these seismic sequences.The in-mine network,constrained by its low dynamic changes,captured only the nucleation phase of the DF60 fault.Furthermore,standard decomposition of the MT solution from the seismic network proved inadequate for accurately identifying the complex nature of the rupture.To enhance safety and risk management in mining environments,we examined the implications of source reactivation within the cluster area post-stress-adjustment.This comprehensive multiscale analysis offers crucial insights into the complex rupture mechanisms and hazards associated with mining-related seismicity.The results underscore the importance of continuous multi-local network monitoring and advanced analytical techniques for improved disaster assessment and risk mitigation in mining operations.
基金Derek Elsworth acknowledges the support from a Gledden Visiting Fellowship from the Institute of Advanced Studies at the University of Western Australia,Australia,and the G.Albert Shoemaker Endowment at Pennsylvania State University,USA.
文摘Triggered seismicity is a key hazard where fluids are injected or withdrawn from the subsurface and may impact permeability. Understanding the mechanisms that control fluid injection-triggered seismicity allows its mitigation. Key controls on seismicity are defined in terms of fault and fracture strength, second-order frictional response and stability, and competing fluid-driven mechanisms for arrest. We desire to constrain maximum event magnitudes in triggered earthquakes by relating pre-existing critical stresses to fluid injection volume to explain why some recorded events are significantly larger than anticipated seismic moment thresholds. This formalism is consistent with several uncharacteristically large fluid injection-triggered earthquakes. Such methods of reactivating fractures and faults by hydraulic stimulation in shear or tensile fracturing are routinely used to create permeability in the subsurface. Microearthquakes (MEQs) generated by such stimulations can be used to diagnose permeability evolution. Although high-fidelity data sets are scarce, the EGS-Collab and Utah FORGE hydraulic stimulation field demonstration projects provide high-fidelity data sets that concurrently track permeability evolution and triggered seismicity. Machine learning deciphers the principal features of MEQs and the resulting permeability evolution that best track permeability changes – with transfer learning methods allowing robust predictions across multiple eological settings. Changes in permeability at reactivated fractures in both shear and extensional modes suggest that permeability change (Δk) scales with the seismic moment (M) of individual MEQs as Δk∝M. This scaling relation is exact at early times but degrades with successive MEQs, but provides a method for characterizing crustal permeability evolution using MEQs, alone. Importantly, we quantify for the first time the role of prestress in defining the elevated magnitude and seismic moment of fluid injection-triggered events, and demonstrate that such MEQs can also be used as diagnostic in quantifying permeability evolution in the crust.
基金funded by the joint fund of the National Key Research and Development Program of China(Grant No.2021YFC2902101)National Natural Science Foundation of China(Grant No.52374084)+1 种基金the 111 Project(Grant No.B17009)DE acknowledges support from the G.Albert Shoemaker endowment.
文摘Understanding the relationship between normal stiffness and permeability in rock fractures under high and true-triaxial in situ stress conditions is critical to assess hydro-mechanical coupling in the Earth's crust.Previous data on stiffness–permeability relations are measured under uniaxial stress states as well as under normal stress.However,many projects involve faulted formations with complex three-dimensional(3D)stress states or significant changes to the original stress state.We rectified this by following the permeability evolution using a true-triaxial stress-permeability apparatus as well as independently applying a spectrum of triaxial stresses from low to high.The relationship between permeability and fracture normal stiffness was quantified using constraints based on the principle of virtual work.The impacts of fracture-lateral and fracture-normal stresses on permeability and normal stiffness evolution were measured.It was found that permeability decreases with increasing fracture-lateral and fracture-normal stresses as a result of Poisson confinement,independent of the orientation of the fracture relative to the stresses.The lateral stresses dominated the evolution of normal stiffness at lower normal stresses(σ_(3)=10 MPa)and played a supplementary role at higher normal stresses(σ_(3)>10 MPa).Moreover,correlations between the evolution of permeability and normal stiffness were extended beyond the low-stiffness,high-permeability region to the high-stiffness,low-permeability region under high fracture-lateral stresses(10–80 MPa)with fracture-normal stress(10–50 MPa)conditions.Again,high lateral stresses further confined the fracture and therefore reduced permeability and increased normal stiffness,which exceeded the previous reported stiffness under no lateral stress conditions.This process enabled us to identify a fundamental change in the flow regime from multi-channel to isolated channelized flow.These results provide important characterizations of fracture permeability in the deep crust,including recovery from deep shale-gas reservoirs.
基金the National Natural Science Foundation of China (Nos. 42302143, 42172159)China Geological Survey Project (No. DD20211350)support from the G. Albert Shoemaker endowment
文摘Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristics complicate the gas-water two-phase flow process in porous media following hydrate decomposition, posing challenges for efficient development. This study examines the transport response of clayey-silt reservoir samples from the Shenhu area using gas-water two-phase flow experiments and CT scanning to explore changes in pore structure, gas-water distribution, and relative permeability under varying flow conditions. The results indicate that pore heterogeneity significantly influences flow characteristics. Gas preferentially displaces water in larger pores, forming fracture-like pores, which serve as preferential flow channels for gas migration. The preferential flow channels enhance gas-phase permeability up to 19 times that of the water phase when fluid pressures exceed total stresses. However,small pores retain liquid, leading to a high residual water saturation of 0.561. CT imaging reveals that these hydro-fractures improve gas permeability but also confine gas flow to specific channels. Pore network analysis shows that gas injection expands the pore-throat network, enhancing connectivity and forming fracture-like pores. Residual water remains trapped in smaller pores and throats, while structural changes, including new fractures, improve gas flow pathways and overall connectivity. Relative permeability curves demonstrate a narrow gas-water cocurrent-flow zone, a right-shifted iso-permeability point and high reservoir capillary pressure, indicating a strong "water-blocking" effect. The findings suggest that optimizing reservoir stimulation techniques to enhance fracture formation, reduce residual water saturation, and improve gas flow capacity is critical for efficient hydrate reservoir development.
基金the Graduate Innovation Program of China University of Mining and Technology,the Fundamental Research Funds for the Central Universities(Grant No.2023WLKXJ017)the Postgraduate Research&Practice Innovation Program of Jiangsu Province(Grant No.KYCX23_2776)the Shandong Energy Group(Grant No.SNKJ2022BJ03-R28)。
文摘In this study,we employed Bayesian inversion coupled with the summation-by-parts and simultaneousapproximation-term(SBP-SAT)forward simulation method to elucidate the mechanisms behind mininginduced seismic events caused by fault slip and their potential effects on rockbursts.Through Bayesian inversion,it is determined that the sources near fault FQ14 have a significant shear component.Additionally,we analyzed the stress and displacement fields of high-energy events,along with the hypocenter distribution of aftershocks,which aided in identifying the slip direction of the critically stressed fault FQ14.We also performed forward modeling to capture the complex dynamics of fault slip under varying friction laws and shear fracture modes.The selection of specific friction laws for fault slip models was based on their ability to accurately replicate observed slip behavior under various external loading conditions,thereby enhancing the applicability of our findings.Our results suggest that the slip behavior of fault FQ14 can be effectively understood by comparing different scenarios.
基金financially supported by the National Key Research and Development Program of China(Grant No.2022YFE0129800)the National Natural Science Foundation of China(Grant No.42202204)。
文摘In-situ upgrading by heating is feasible for low-maturity shale oil,where the pore space dynamically evolves.We characterize this response for a heated substrate concurrently imaged by SEM.We systematically follow the evolution of pore quantity,size(length,width and cross-sectional area),orientation,shape(aspect ratio,roundness and solidity)and their anisotropy—interpreted by machine learning.Results indicate that heating generates new pores in both organic matter and inorganic minerals.However,the newly formed pores are smaller than the original pores and thus reduce average lengths and widths of the bedding-parallel pore system.Conversely,the average pore lengths and widths are increased in the bedding-perpendicular direction.Besides,heating increases the cross-sectional area of pores in low-maturity oil shales,where this growth tendency fluctuates at<300℃ but becomes steady at>300℃.In addition,the orientation and shape of the newly-formed heating-induced pores follow the habit of the original pores and follow the initial probability distributions of pore orientation and shape.Herein,limited anisotropy is detected in pore direction and shape,indicating similar modes of evolution both bedding-parallel and bedding-normal.We propose a straightforward but robust model to describe evolution of pore system in low-maturity oil shales during heating.
基金supported by National Nature Science Foundation of China (Grant No.42177157)the Science and Technology Program of Liaoning Province (Grant No.2023JH1/10400003)the Applied Basic Research Programof Liaoning Province (Grant No.2023JH2/101300153).
文摘The injection of large volumes of natural gas into geological formations,as is required for underground gas storage,leads to alterations in the effective stress exerted on adjacent faults.This increases the potential for their reactivation and subsequent earthquake triggering.Most measurements of the frictional properties of rock fractures have been conducted under normal and shear stresses.However,faults in gas storage facilities exist within a true three-dimensional(3D)stress state.A double-direct shear experiment on rock fractures under both lateral and normal stresses was conducted using a true triaxial loading system.It was observed that the friction coefficient increases with increasing lateral stress,but decreases with increasing normal stress.The impact of lateral and normal stresses on the response is primarily mediated through their influence on the initial friction coefficient.This allows for an empirical modification of the rate-state friction model that considers the influence of lateral and normal stresses.The impact of lateral and normal stresses on observed friction coefficients is related to the propensity for the production of wear products on the fracture surfaces.Lateral stresses enhance the shear strength of rock(e.g.Mogi criterion).This reduces asperity breakage and the generation of wear products,and consequently augments the friction coefficient of the surface.Conversely,increased normal stresses inhibit dilatancy on the fracture surface,increasing the breakage of asperities and the concomitant production of wear products that promote rolling deformation.This ultimately reduces the friction coefficient.
基金funded by the National Natural Science Foundation of China(Nos.42320104003 and 42107163)the Funda mental Research Funds for the Central Universities.Derek Elsworth acknowledges support from the G.Albert Shoemaker endowment.
文摘Basalt is a major component of the earth and moon crust.Mineral composition and temperature influence frictional instability and thus the potential for seismicity on basaltic faults.We performed velocitystepping shear experiments on basalt gouges at a confining pressure of 100 MPa,temperatures in the range of 100-400℃ and with varied obsidian mass fractions of 0-100%under wet/dry conditions to investigate the frictional strength and stability of basaltic faults.We observe a transition from velocity-neutral to velocity-weakening behaviors with increasing obsidian content.The frictional stability response of the mixed obsidian/basalt gouges is characterized by a transition from velocitystrengthening to velocity-weakening at 200℃ and another transition to velocity-strengthening at temperatures>300℃.Conversely,frictional strengths of the obsidian-bearing gouges are insensitive to temperature and wet/dry conditions.These results suggest that obsidian content dominates the potential seismic response of basaltic faults with the effect of temperature controlling the range of seismogenic depths.Thus,shallow moonquakes tend to occur in the lower lunar crust due to the corresponding anticipated higher glass content and a projected temperature range conducive to velocity-weakening behavior.These observations contribute to a better understanding of the nucleation mechanism of shallow seismicity in basaltic faults.
基金financially supported by the National Key Research and Development Program of China (Grant No. 2022YFE0129800)the National Natural Science Foundation of China (Grant No. 42202204)support from the G. Albert Shoemaker endowment。
文摘In-situ thermal upgrading is used to tune the pore system in low-maturity oil shales. We introduce fractal dimension(D), form factor(ff) and stochastic entropy(H) to quantify the heating-induced evolution of pore morphological complexity and azimuthal disorder and develop a model to estimate the impact on seepage capacity via permeability. Experiments are conducted under recreated in-situ temperatures and consider anisotropic properties—both parallel and perpendicular to bedding. Results indicate that azimuthal distribution of pores in the bedding-parallel direction are dispersed, while those in the bedding-perpendicular direction are concentrated. D values indicate that higher temperatures reduce the uniformity of the pore size distribution(PSD) in the bedding-parallel direction but narrow the PSD in the bedding-perpendicular direction. The greater ff(> 0.7) values in the bedding-parallel direction account for a large proportion, while the dominated in the bedding-perpendicular direction locates within 0.2-0.7, for all temperatures. The H value of the bedding-parallel sample remains stable at ~0.925 during heating, but gradually increases from 0.808 at 25℃ to 0.879 at 500℃ for the beddingperpendicular sample. Congruent with a mechanistic model, the permeability at 500℃ is elevated~1.83 times(bedding-parallel) and ~6.08 times(bedding-perpendicular) relative to that at 25℃—confirming the effectiveness of thermal treatment in potentially enhancing production from low-maturity oil shales.
基金supported by the Australian Research Council(Grant No.DP200101293)supported by the UWA-China Joint Scholarships(201906430030).
文摘In current dual porosity/permeability models,there exists a fundamental assumption that the adsorption-induced swelling is distributed uniformly within the representative elementary volume (REV),irrespective of its internal structures and transient processes.However,both internal structures and transient processes can lead to the non-uniform swelling.In this study,we hypothesize that the non-uniform swelling is responsible for why coal permeability in experimental measurements is not only controlled by the effective stress but also is affected by the adsorption-induced swelling.We propose a concept of the swelling triangle composed of swelling paths to characterize the evolution of the non-uniform swelling and serve as a core link in coupled multiphysics.A swelling path is determined by a dimensionless volumetric ratio and a dimensionless swelling ratio.Different swelling paths have the same start and end point,and each swelling path represents a unique swelling case.The swelling path as the diagonal of the triangle represents the case of the uniform swelling while that as the two perpendicular boundaries represents the case of the localized swelling.The paths of all intermediate cases populate inside the triangle.The corresponding relations between the swelling path and the response of coal multiphysics are established by a non-uniform swelling coefficient.We define this method as the triangle approach and corresponding models as swelling path-based ones.The proposed concept and models are verified against a long-term experimental measurement of permeability and strains under constant effective stress.Our results demonstrate that during gas injection,coal multiphysics responses have a close dependence on the swelling path,and that in both future experiments and field predictions,this dependence must be considered.
基金National Science and Technology Major Project(2017ZX05035004-002)the National Natural Science Foundation of China(41802183)+1 种基金the National Postdoctoral Program for Innovative Talents(BX201700282)the China Postdoctoral Science Foundation(2017M621870),the project of investigation and evaluation of special coal resources(DD20I60187-1).
文摘Coal measure source rocks, located in the Xihu Sag of the East China Sea Shelf Basin, were analyzed to define the hydrocarbon generation potential, organic geochemistry/petrology characteristics, and coal preservation conditions. The Pinghu source rocks in the Xihu Sag are mainly gas-prone accompany with condensate oil generation. The coals and shales of the Pinghu Formation are classified from "fair" to "excellent" source rocks with total organic carbon(TOC) contents ranging from 25.2% to 77.2% and 1.29% to 20.9%, respectively. The coals are richer in TOC and S1+S2 than the shales, indicating that the coals have more generation potential per unit mass. Moreover, the kerogen type of the organic matter consists of types Ⅱ-Ⅲ and Ⅲ, which the maturity Ro ranges from 0.59% to 0.83%. Petrographically, the coals and shales are dominated by vitrinite macerals(69.1%–96.8%) with minor proportions of liptinite(2.5%–17.55%) and inertinite(0.2%–6.2%). The correlation between maceral composition and S1+S2 indicates that the main contributor to the generation potential is vitrinite. Therefore, the coals and shales of the Pinghu Formation has good hydrocarbon generation potential, which provided a good foundation for coal measure gas accumulation. Furthermore, coal facies models indicates that the Pinghu coal was deposited in limno-telmatic environment under high water levels, with low tree density(mainly herbaceous) and with low-moderate nutrient supply. Fluctuating water levels and intermittent flooding during the deposition of peat resulted in the inter-layering of coal, shale and sandstone, which potentially providing favorable preservation conditions for coal measure gas.
基金supported by National Key R&D Program of China(Grant No.2018YFC0407006)the 111 Project(Grant No.B17009)the Australian Research Council(Grant No.DP200101293)。
文摘Coal permeability measurements are normally conducted under the assumption that gas pressure in the matrix is equalized with that in fracture and that gas sorption-induced swelling/shrinking strain is uniformly distributed within the coal.However,the validity of this assumption has long been questioned and differential strain between the fracture strain and the bulk strain has long been considered as the primary reason for the inconsistency between experimental data and poroelasticity solutions.Although efforts have been made to incorporate the impact into coal permeability models,the fundamental nature of those efforts to split the matrix strain between fracture and coal bulk remains questionable.In this study,a new concept of differential swelling index(DSI)was derived to theoretically define the relation among sorption-induced strains of the coal bulk,fracture,and coal matrix at the equilibrium state.DSI was a function of the equilibrium pressure and its magnitudes were regulated by the Langmuir constants of both the matrix and the coal bulk.Furthermore,a spectrum of DSI-based coal permeability models was developed to explicitly consider the effect of differential strains.These models were verified with the experimental data under the conditions of uniaxial strain,constant confining pressure,and constant effective stress.
基金The financial support from the National Natural Science Foundation of China(Nos.41772154 and 42102338)Natural Science Foundation of Shandong Province(Nos.ZR2019MA009 and ZR2020QE115)SDUST Research Fund of China(No.2018TDJH102)。
文摘Heat production from geothermal reservoirs is a typical heat transfer process involving a cold working fluid contacting a hot rock formation.Compared to the thermal-physical characteristics of water,supercritical CO_(2)(scCO_(2))has a higher heat storage capacity over a wide temperature-pressure range and may be favored as a heat transfer fluid.Singularly characteristic of scCO_(2)-based heat extraction is that the hydraulic-thermal properties of the scCO_(2) vary dramatically and dynamically with the spatial pressure gradient during unsteady-state flow along fracture.This highly nonlinear behavior presents a challenge in the accurate estimation of heat extraction efficiency in scCO_(2)-based EGS.In this paper,a thermal-h ydraulic-mechanical(THM)coupled model is developed by considering deformation of the fractured reservoir,non-Darcy flow and the varying thermal-physical properties of scCO_(2).The proposed model is validated by matching the modeling temperature distribution with published data.The results show that during continuous injection of scCO_(2),the fracture first widens and then narrows,ultimately reopening over the long term.The sequential fracture deformation behaviors are in response to the combined impacts of mechanical compression and thermally-induced deformation.By controlling the injection parameters of the scCO_(2),it is found that the heat extraction rate is positively correlated to its pore pressure or mass flow rate.The heat extraction rate can be significantly enhanced,when the inlet temperature of scCO_(2) is below its critical temperature.As a result,the heat increment recovered per unit mass of scCO_(2) decreases as the hot rock is gradually cooled.Meanwhile,the heat increment recovered per unit mass of scCO_(2) decreases by increasing the inlet temperature of scCO_(2) or its mass flow rate,but increases as the outlet pressure rises.Furthermore,multi-linear regression indicates that controlling the inlet temperature of the scCO_(2) can significantly improve the thermodynamic efficiency of heat extraction.
基金support provided by United States Department of Energy Grant DE-FE0023354。
文摘Subsurface fluid injections can disturb the effective stress regime by elevating pore pressure and potentially reactivate faults and fractures.Laboratory studies indicate that fracture rheology and permeability in such reactivation events are linked to the roughness of the fracture surfaces.In this study,we construct numerical models using discrete element method(DEM)to explore the influence of fracture surface roughness on the shear strength,slip stability,and permeability evolution during such slip events.For each simulation,a pair of analog rock coupons(three-dimensional bonded quartz particle analogs)representing a mated fracture is sheared under a velocity-stepping scheme.The roughness of the fracture is defined in terms of asperity height and asperity wavelength.Results show that(1)Samples with larger asperity heights(rougher),when sheared,exhibit a higher peak strength which quickly devolves to a residual strength after reaching a threshold shear displacement;(2)These rougher samples also exhibit greater slip stability due to a high degree of asperity wear and resultant production of wear products;(3)Long-term suppression of permeability is observed with rougher fractures,possibly due to the removal of asperities and redistribution of wear products,which locally reduces porosity in the dilating fracture;and(4)Increasing shear-parallel asperity wavelength reduces magnitudes of stress drops after peak strength and enhances fracture permeability,while increasing shear-perpendicular asperity wavelength results in sequential stress drops and a delay in permeability enhancement.This study provides insights into understanding of the mechanisms of frictional and rheological evolution of rough fractures anticipated during reactivation events.
基金the financial support by the National Key Research and Development Program of China(Grant No.2020YFC1808102)the Natural Science Foundation of China(No.51874328 and No.52074311)。
文摘Optimization of fracturing perforation is of great importance to the commingling gas production in coal measure strata.In this paper,a 3 D lattice algorithm hydraulic fracturing simulator was employed to study the effects of perforation position and length on hydraulic fracture propagation in coal measures of the Lin-Xing block,China.Based on field data,three lithologic combinations are simulated:1)a thick section of coal seam sandwiched by sandstones;2)a thin coal seam layer overlay by gas-bearing tight sandstone;3)two coal seams separated by a thin layer of sandstone.Our simulation shows that perforation position and length in multi-layer reservoirs play a major role in hydraulic fracture propagation.Achieving maximum stimulated volume requires consideration of lithologic sequence,coal seam thickness,stress states,and rock properties.To improve the combined gas production in coal measure strata,it is possible to simultaneously stimulate multiple coal seams or adjacent gas-bearing sandstones.In these cases,perforation location and length also significantly impact fracture propagation,and therefore should be carefully designed.Our simulation results using 3 D lattice algorithm are qualitatively consistent with laboratory physical simulation.3 D lattice models can be used to effectively simulate the fracture propagation through layers in coal measure strata.The numerical results provide guidance for perforation optimization in the hydraulic fracturing of coal measure strata.
基金supported by the Climbling Project of Taishan Scholar in Shandong Province (No.tspd20210313)National Natural Science Foundation of China (Grant No.51874190,52079068,41941019,52090081 and 52074168)+3 种基金Taishan Scholar in Shandong Province (No.tsqn202211150)Outstanding Youth Fund Project in Shandong Province (No.ZQ2022YQ49)the State Key Laboratory of Hydroscience and Engineering,China (No.2021-KY-04)support from the G.Albert Shoemaker endowment.
文摘Multi-seam mining often leads to the retention of a significant number of coal pillars for purposes such as protection,safety,or water isolation.However,stress concentration beneath these residual coal pillars can significantly impact their strength and stability when mining below them,potentially leading to hydraulic support failure,surface subsidence,and rock bursting.To address this issue,the linkage between the failure and instability of residual coal pillars and rock strata during multi-seam mining is examined in this study.Key controls include residual pillar spalling,safety factor(f.),local mine stiffness(LMS),and the post-peak stiffness(k)of the residual coal pillar.Limits separating the two forms of failure,progressive versus dynamic,are defined.Progressive failure results at lower stresses when the coal pillar transitions from indefinitely stable(f,>1.5)to failing(f,<1.5)when the coal pillar can no longer remain stable for an extended duration,whereas sud-den(unstable)failure results when the strength of the pillar is further degraded and fails.The transition in mode of failure is defined by the LMS/k ratio.Failure transitions from quiescent to dynamic as LMS/k.<1,which can cause chain pillar instability propagating throughout the mine.This study provides theoretical guidance to define this limit to instability of residual coal pillars for multi-seam mining in similar mines.
基金supported by the Australian Research Council under Grant DP200101293UWA China Scholarshipsthe China Scholarship Council(CSC No.201707970011)。
文摘Stimulated shale reservoirs consist of kerogen,inorganic matter,secondary and hydraulic fractures.The dispersed distribution of kerogen within matrices and complex gas flow mechanisms make production evaluation challenging.Here we establish an analytical method that addresses kerogen-inorganic matter gas transfer,dispersed kerogen distribution,and complex gas flow mechanisms to facilitate evaluating gas production.The matrix element is defined as a kerogen core with an exterior inorganic sphere.Unlike most previous models,we merely use boundary conditions to describe kerogen-inorganic matter gas transfer without the instantaneous kerogen gas source term.It is closer to real inter-porosity flow conditions between kerogen and inorganic matter.Knudsen diffusion,surface diffusion,adsorption/desorption,and slip corrected flow are involved in matrix gas flow.Matrix-fracture coupling is realized by using a seven-region linear flow model.The model is verified against a published model and field data.Results reveal that inorganic matrices serve as a major gas source especially at early times.Kerogen provides limited contributions to production even under a pseudo-steady state.Kerogen properties’influence starts from the late matrix-fracture inter-porosity flow regime,while inorganic matter properties control almost all flow regimes except the early-mid time fracture linear flow regime.The contribution of different linear flow regions is also documented.
基金supported by the National Natural Science Foundation of China (Nos. 41672268 and 41772286)the Petro China Innovation Foundation (2018D-5007-0202)the US Department of Energy (DEFE0023354)
文摘This paper experimentally explores the frictional sliding behavior of two simulated gouges:one,a series of quartz–smectite mixtures,and the other,powdered natural rocks,aiming to evaluate and codify the effect of mineralogy on gouge dilation and frictional strength,stability,and healing.Specifically,velocity-stepping and slide-hold-slide experiments were performed in a double direct shear configuration to analyze frictional constitutive parameters at room temperature,under normal stresses of 10,20,and 40 MPa.Gouge dilation was measured based on the applied step-wise changes in shear velocity.The frictional response of the quartz–smectite mixtures and powdered natural rocks are affected by their phyllosilicate content.Frictional strength and healing rates decrease with increasing phyllosilicate content,and at 20 wt.%a transition from velocity-weakening to velocity-strengthening behavior was noted.For both suites of gouges,dilation is positively correlated with frictional strength and healing rates,and negatively correlated with frictional stability.Changes in the permeability of gouge-filled faults were estimated from changes in mean porosity,indexed through measured magnitudes of gouge dilation.This combined analysis implies that the reactivation of caprock faults filled with phyllosilicaterich gouges may have a strong influence on permeability evolution in caprock faults.
