Gouge stability controlled by temperature elevation and obsidian addition in basaltic faults and implications for moonquakes

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.

Thermo-hydro-mechanical (THM) coupled simulation of the land subsidence due to aquifer thermal energy storage (ATES) system in soft soils

Aquifer thermal energy storage(ATES)system has received attention for heating or cooling buildings.However,it is well known that land subsidence becomes a major environmental concern for ATES projects.Yet,the effect of temperature on land subsidence has received practically no attention in the past.This paper presents a thermo-hydro-mechanical(THM)coupled numerical study on an ATES system in Shanghai,China.Four water wells were installed for seasonal heating and cooling of an agriculture greenhouse.The target aquifer at a depth of 74e104.5 m consisted of alternating layers of sand and silty sand and was covered with clay.Groundwater level,temperature,and land subsidence data from 2015 to 2017 were collected using field monitoring instruments.Constrained by data,we constructed a field scale three-dimensional(3D)model using TOUGH(Transport of Unsaturated Groundwater and Heat)and FLAC3D(Fast Lagrangian Analysis of Continua)equipped with a thermo-elastoplastic constitutive model.The effectiveness of the numerical model was validated by field data.The model was used to reproduce groundwater flow,heat transfer,and mechanical responses in porous media over three years and capture the thermo-and pressure-induced land subsidence.The results show that the maximum thermoinduced land subsidence accounts for about 60%of the total subsidence.The thermo-induced subsidence is slightly greater in winter than that in summer,and more pronounced near the cold well area than the hot well area.This study provides some valuable guidelines for controlling land subsidence caused by ATES systems installed in soft soils.

A review of methods,applications and limitations for incorporating fluid flow in the discrete element method

The past decade has witnessed the substantial growth in research interests and progress on the subject of coupled hydro-mechanical processes in rocks and soils,driven mainly by the surge of research in unconventional hydrocarbon reservoirs and associated hazards.Many coupling techniques have been developed to include the effects of fluid flow in the discrete element method(DEM),and the techniques have been applied to a variety of geomechanical problems.Although these coupling methods have been successfully applied in various engineering fields,no single fluid/DEM coupling method is universal due to the complexity of engineering problems and the limitations of the numerical methods.For researchers and engineers,the key to solve a specific problem is to select the most appropriate fluid/DEM coupling method among these modeling technologies.The purpose of this paper is to give a comprehensive review of fluid flow/DEM coupling methods and relevant research.Given their importance,the availability or unavailability of best practice guidelines is outlined.The theoretical background and current status of DEM are introduced first,and the principles,applications,and advantages and disadvantages of different fluid flow/DEM coupling methods are discussed.Finally,a summary with speculation on future development trends is given.

Three-dimensional distinct element modeling of fault reactivation and induced seismicity due to hydraulic fracturing injection and backflow

This paper presents a three-dimensional fully hydro-mechanical coupled distinct element study on fault reactivation and induced seismicity due to hydraulic fracturing injection and subsequent backflow process,based on the geological data in Horn River Basin,Northeast British Columbia,Canada.The modeling results indicate that the maximum magnitude of seismic events appears at the fracturing stage.The increment of fluid volume in the fault determines the cumulative moment and maximum fault slippage,both of which are essentially proportional to the fluid volume.After backflow starts,the fluid near the joint intersection keeps flowing into the critically stressed fault,rather than backflows to the wellbore.Although fault slippage is affected by the changes of both pore pressure and ambient rock stress,their contributions are different at fracturing and backflow stages.At fracturing stage,pore pressure change shows a dominant effect on induced fault slippage.While at backflow stage,because the fault plane is under a critical stress state,any minor disturbance would trigger a fault slippage.The energy analysis indicates that aseismic deformation takes up a majority of the total deformation energy during hydraulic fracturing.A common regularity is found in both fracturing-and backflow-induced seismicity that the cumulative moment and maximum fault slippage are nearly proportional to the injected fluid volume.This study shows some novel insights into interpreting fracturing-and backflowinduced seismicity,and provides useful information for controlling and mitigating seismic hazards due to hydraulic fracturing.

Influence of fracture roughness on shear strength,slip stability and permeability:A mechanistic analysis by three-dimensional digital rock modeling

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.

Developing cold-resistant high-adhesive electronic substrate for WIMPs detectors at CDEX

Herein we report a prototypical electronic substrate specifically designed to serve the weakly interacting massive particles(WIMPs)detectors at the China Dark Matter Experiment(CDEX).Because the bulky high-purity germanium(HPGe)detectors operate under liquid-nitrogen temperatures and ultralow radiation backgrounds,the desired electronic substrates must maintain high adhesivity across different layers in such cold environment and be free from any radioactive nuclides.To conquer these challenges,for the first time,we employed polytetrafluoroethylene((C2F4)n)foil as the base substrate,in conjunction with ion implantation and deposition techniques using an independently developed device at Beijing Normal University for surface modification prior to electroplating.The remarkable peeling strengths of 0.88±0.06 N/mm for as-prepared sample and 0.75±0.05 N/mm for that after 2.5-days of soaking inside the liquid nitrogen were observed,while the regular standards commonly require 0.4 N/mm^0.6 N/mm for electronic substrates.

Effect of mineralogy on friction-dilation relationships for simulated faults:Implications for permeability evolution in caprock faults

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.

Experimental investigation on frictional properties of stressed basalt fractures

The frictional strength and sliding stability of faults are crucial in interpreting earthquake mechanisms and cycles.Herein,we report friction experiments on basalt fractures,using a self-designed triaxial apparatus that allows direct shear of samples under coupled hydro-mechanical conditions.Velocitystepping(VS)and slide-hold-slide(SHS)experiments are performed on both bare and gouge-bearing surfaces of Xiashan basalt subjected to cyclic shear velocities at 1e30 mm/s,effective normal stresses of 1e5 MPa,and pore pressures of 70e300 kPa.The measured basalt friction coefficients are in the range of 0.67e0.74,which is sensitive to gouge thickness,normal stress,and water.Specifically,a reduction in friction coefficient is observed with an increment in gouge thickness,normal stress,and pore pressure.Based on the microscopic observation of the pre-and post-shearing sliding surfaces,this weakening effect in friction coefficient can be attributed to powder lubrication.Furthermore,the VS test results reveal predominantly velocity-strengthening behavior at investigated slip velocities,and this velocity strengthening behavior does not appear to be influenced by variations in normal stress,gouge thickness,and water.However,changes in sliding velocity and normal stress can lead to a shift between stable and unstable sliding.Specifically,stable sliding is favored by high sliding velocities and low normal stress applied in this study.Finally,we analyze the experimental data by calculating the rate-and-state parameters using the rate-and state-dependent friction(RSF)theory.Importantly,the calculated friction rate parameter(a-b)supports the velocity-strengthening behavior.Both frictional relaxation(Dmc)during hold periods and frictional healing(Dm)upon re-shearing are linearly proportional to the logarithmic hold time,which may be attributed to the growth in true contact area with hold time.This study sheds light on the roles of sliding velocity,and gouge thickness in controlling frictional strength and stability of basalt fractures.

Frictional stability of Longmaxi shale gouges and its implication for deep seismic potential in the southeastern Sichuan Basin

Microearthquakes accompanying shale gas recovery highlight the importance of exploring the frictional and stability properties of shale gouges.Aiming to reveal the influencing factors on fault stability,this paper explores the impact of mineral compositions,effective stress and temperature on the frictional stability of Longmaxi shale gouges in deep reservoirs located in the Luzhou area,southeastern Sichuan Basin.Eleven shear experiments were conducted to define the frictional strength and stability of five shale gouges.The specific experimental conditions were as follows:temperatures:90–270°C;a confining stress:95 MPa;and pore fluid pressures:25–55 MPa.The results show that all five shale gouges generally display high frictional strength with friction coefficients ranging from 0.60 to 0.70 at the aforementioned experiment condition of pressures,and temperatures.Frictional stability is significantly affected by temperature and mineral compositions,but is insensitive to variation in pore fluid pressures.Fault instability is enhanced at higher temperatures(especially at>200°C)and with higher tectosilicate/carbonate contents.The results demonstrate that the combined effect of mineral composition and temperature is particularly important for induced seismicity during hydraulic fracturing in deep shale reservoirs.

Biological Effects of Sunflower Seeds Implanted by Carbon Ion

[ Objective ] This study aimed to investigate the effects of carbon ion implantation and implantation times on growth and genetic variation of sunflowers. [ Method] Carbon ions were implanted into Bakui 138, Bakui i36 and Bakui 118 seeds at dose of 5 - 10is C/cm2, before they were planted. Their Fl-generation seeds were irradiated again. Seeds of the both generations were planted and the growth d the seedlings was observed in field tests. Finally, their genetic variation was analyzed through RAPD. [ Result] The germination rate and several agronomic traits like plant height, stem diameter, leaf number and yields of Bakui 138 of once-irradiated group were significantly improved, while that of twice-irradiated group showed opposite trend. The variation of Bakui 136 and Bakui 118 was insig- nificant. At the molecular level, the genetic distance with the control group of once and twice-irradiated groups was 0. 111 1, 0. 108 7 in Bakui 138; 0. 068 O, O. 030 3 in Bakui 136 and 0.062 5,0.043 5 in Bakui 118. [Conclusion] Carbon ion implantation had a significant effect on the growth and development of Bakui 138, and the effect varied with irradiation times. Moreover, it caused genomic variation in the three sunflower cuhivars.

Probing dark matter particles from evaporating primordial black holes via electron scattering in the CDEX-10 experiment

Dark matter(DM)is a major constituent of the Universe.However,no definite evidence of DM particles(denoted as“χ”)has been found in DM direct detection(DD)experiments to date.There is a novel concept of detectingχfrom evaporating primordial black holes(PBHs).We search forχemitted from PBHs by investigating their interaction with target electrons.The examined PBH masses range from 1×10^(15)to 7×10^(16)g under the current limits of PBH abundance fPBH.Using 205.4 kg·day data obtained from the CDEX-10 experiment conducted in the China Jinping Underground Laboratory,we exclude theχ-electron(χ-e)elastic-scattering cross sectionσ_(χe)~5×10^(-29)cm^(2)forχwith a mass■keV from our results.With the higher radiation background but lower energy threshold(160 eV),CDEX-10 fills a part of the gap in the previous work.If(m_(χ),σ_(χe))can be determined in the future,DD experiments are expected to impose strong constraints on fPBHfor large MPBHs.

Finite domain solution of a KGD hydraulic fracture in the viscosity-dominated regime

This paper describes a numerical algorithm for solving the classic problem of a plane strain(KGD)fracture propagating in an impermeable elastic medium with zero toughness.The method,which takes advantage of the self-similar nature of the solution,combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation.This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects,noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method.To enhance the efficiency and accuracy of the numerical scheme,the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening,while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation.The proposed technique yields results that closely match the analytical solution,even with a coarse mesh.This approach offers potential for addressing more complex hydraulic fracturing problems in the future.

Two-phase flow thermo-hydro-mechanical modeling for a water flooding field case

Simulation of subsurface energy system involves multi-physical processes such as thermal,hydraulical,and mechanical(THM)processes,and requires a so-called THM coupled modeling approach.THM coupled modeling is commonly performed in geothermal energy production.However,for hydrocarbon extraction,we need to consider multiphase flow additionally.In this paper,we describe a three-dimensional numerical model of non-isothermal two-phase flow in the deformable porous medium by integrating governing equations of two-phase mixture in the porous media flow in the reservoir.To account for inter-woven impacts in subsurface condi-tions,we introduced a temperature-dependent fluid viscosity and a fluid density along with a strain-dependent reservoir permeability.Subsequently,we performed numerical experiments of a ten-year water flooding pro-cess employing the open-source parallelized code,OpenGeoSys.We considered different well patterns with colder water injection in realistic scenarios.Our results demonstrate that our model can simulate complex interactions of temperature,pore pressure,subsurface stress and water saturation simultaneously to evaluate the recovery per-formance.High temperature can promote fluid flow while cold water injection under non-isothermal conditions causes the normal stress reduction by significant thermal stress.Under different well patterns the displacement efficiency will be changed by the relative location between injection and production wells.This finding has provided the important reference for fluid flow and induced stress evolution during hydrocarbon exploitation under the environment of large reservoir depth and high temperature.

First experimental constraints on WIMP couplings in the effective field theory framework from CDEX

We present weakly interacting massive particles(WIMPs) search results performed using two approaches of effective field theory from the China Dark Matter Experiment(CDEX), based on the data from both CDEX-1B and CDEX-10 stages. In the nonrelativistic effective field theory approach, both time-integrated and annual modulation analyses were used to set new limits for the coupling of WIMP-nucleon effective operators at 90% confidence level(C.L.) and improve over the current bounds in the low mχregion. In the chiral effective field theory approach, data from CDEX-10 were used to set an upper limit on WIMP-pion coupling at 90% C.L. We for the first time extended the limit to the m_(χ)<6 GeV/c^(2) region.

Performances of a prototype point-contact germanium detector immersed in liquid nitrogen for light dark matter search

The CDEX-10 experiment searches for light weakly interacting massive particles, a form of dark matter, at the China Jinping Underground Laboratory, where approximately 10 kg of germanium detectors are arranged in an array and immersed in liquid nitrogen. Herein, we report on the experimental apparatus, detector characterization, and spectrum analysis of one prototype detector. Owing to the higher rise-time resolution of the CDEX-10 prototype detector as compared with CDEX-1 B, we identified the origin of an observed category of extremely fast events. For data analysis of the CDEX-10 prototype detector, we introduced and applied an improved bulk/surface event discrimination method. The results of the new method were compared to those of the CDEX-1 B spectrum. Both sets of results showed good consistency in the 0-12 ke Vee energy range, except for the 8.0 keV K-shell X-ray peak from the external copper.

Evaluation of fault stability and seismic potential for Hutubi underground gas storage due to seasonal injection and extraction

The Hutubi gas field was put into production in 1998 and then converted into an underground gas storage(UGS)facility in 2013,and since then a cluster of earthquakes associated with seasonal injection and extraction activities have been recorded nearby.To evaluate the fault stability and seismic potential,we established a pseudo-3D geomechanical model to simulate the process of seasonal injection and extraction.Reservoir pore pressures from 1998 to 2019 were obtained through multiphase reservoir simulation and validated by history matching the field injection and production data.We then imported pore pressures into the geomechanical model to simulate the poroelastic perturbation on faults for over 20 years.The fidelity of this model was validated by comparing the simulated surface deformation with global positioning system(GPS)measured data.We used Coulomb failure stress(CFS)as the indicator for the likelihood of fault slippage.The simulation results show that the location of the induced earthquake cluster was within the positive Coulomb stress perturbation(DCFS)area,in which fault slippage was promoted.In addition,DCFS at the earthquake location kept increasing after the injection began.These findings could explain the induced earthquakes with the Coulomb failure stress theory.Furthermore,we conducted a parameter sensitivity study on the dominant factors such as the maximum operating pressure(MOP),frictional coefficient,and dip angle of the pre-existing fault.The results indicate that the magnitude of DCFS caused by seasonal injection and extraction decreases with distance;MOPs are constrained to 32.9,36.2,and 39.5 MPa according to different DCFS thresholds;the critical dip angle ranges are 0-20°and 80°-100°;and strengthening the fault friction can either increase or decrease the seismic potential.This study can help determine the MOP for Hutubi underground gas storage(HTB UGS)and provide a framework for simulating the potential causes of induced seismicity for other sites.