Stress initialization methods for dynamic numerical simulation of rock mass with high in-situ stress

In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.

A comprehensive review of in-situ polymer gel simulation for conformance control

Gel treatment has been widely applied in mature oilfields to improve sweep efficiency and control water production.Correct numerical simulation is of major importance to the optimization design and prediction of a successful gel treatment.However,there exist many problems in current simulation studies in published literature.This paper first presents a comprehensive review on the major factors that have been considered at different gelation stages during gel treatment,the models used in the commercial/inhouse simulators,and current numerical simulation studies on both laboratory and field scales.Then we classify the current in-situ gel numerical simulation problems as 1,deficient model problem that has published numerical model but has not been applied in simulator and application studies;2,missing model problem that does not have published quantitative model;and 3,inaccurate application problem that does not consider the major factors of gel performance,based on the reasons from some questionable results of current simulation studies.Finally,we point out the major research efforts that should be made in the future to better simulate in-situ gel treatment process.The review indicates that numerous simulation studies using commercial software packages intend to predigest the gel treatment,many of which,however,ignore important mechanisms and mislead the operation of gel treatment.In fact,a full assessment of simulating in-situ gels cannot be achieved unless the quantitative models can be qualified in terms of transport and plugging mechanisms based on the experimental results.

Numerical simulation of oil shale in-situ mining using fluid-thermo-solid coupling

For the in-situ gas-injection mining technology of oil shale by process,a numerical simulation method with flow-thermo-solid coupling is proposed in this paper.This method adopts separate simulations and step-by-step coupling simulation ideas combined with the advantages of the finite element method and the finite volume method.The numerical simulation of flow-thermo-solid coupling is decomposed into two parts:flow-solid coupling and thermo-solid coupling.Considering the Fuyu oil shale in-situ production test area in Songliao Basin as an example,it is concluded that the oil shale has undergone four heating stages:a rapid temperature rise,a steady temperature rise,a slow temperature rise,and heat preservation.It takes about 10 years for the stress-strain state of the oil shale layer to reach a steady-state through the thermo-solid coupling.The main strain zones of the oil shale layer are distributed near the fracturing fractures connected to the gas injection well and at the edge of the fracturing fractures.The areas with the plastic deformation in the oil shale layer predominantly appear near the gas injection wells,production wells,and fracturing channels.The areas with the largest fracture strength are mostly distributed near the edge of the fracturing fractures with low flow velocity and low temperature.

Solute drag-controlled grain growth in magnesium investigated by quasi in-situ orientation mapping and level-set simulations

Critical properties of metallic materials,such as the yield stress,corrosion resistance and ductility depend on the microstructure and its grain size and size distribution.Solute atoms that favorably segregate to grain boundaries produce a pinning atmosphere that exerts a drag pressure on the boundary motion,which strongly affects the grain growth behavior during annealing.In the current work,the characteristics of grain growth in an annealed Mg-1 wt.%Mn-1 wt.%Nd magnesium alloy were investigated by advanced experimental and modeling techniques.Systematic quasi in-situ orientation mappings with a scanning electron microscope were performed to track the evolution of local and global microstructural characteristics as a function of annealing time.Solute segregation at targeted grain boundaries was measured using three-dimensional atom probe tomography.Level-set computer simulations were carried with different setups of driving forces to explore their contribution to the microstructure development with and without solute drag.The results showed that the favorable growth advantage for some grains leading to a transient stage of abnormal grain growth is controlled by several drivers with varying importance at different stages of annealing.For longer annealing times,residual dislocation density gradients between large and smaller grains are no longer important,which leads to microstructure stability due to predominant solute drag.Local fluctuations in residual dislocation energy and solute concentration near grain boundaries cause different boundary segments to migrate at different rates,which affects the average growth rate of large grains and their evolved shape.

Kinetic simulation of hydrocarbon generation and its application to in-situ conversion of shale oil

The kinetic parameters of hydrocarbon generation are determined through experimental simulation and mathematical calculation using four typical samples selected from the Cretaceous Nenjiang Formation in the northwest of Songliao Basin,Chang 7 Member of Triassic Yanchang Formation in the southwest of Ordos Basin,Paleogene in the southwest of Qaidam Basin,and Lucaogou Formation of Jimusar Sag in the east of Junggar Basin.The results show that activation energy of hydrocarbon generation of organic matter is closely related to maturity and mainly ranges between 197 kJ/mol and 227 kJ/mol.On this basis,the temperature required for organic matter in shale to convert into oil was calculated.The ideal heating temperature is between 270℃and 300℃,and the conversation rate can reach 90%after 50-300 days of heating at constant temperature.When the temperature rises at a constant rate,the temperature corresponding to the major hydrocarbon generation period ranges from 225 to 350℃at the temperature rise rate of 1-150℃/month.In order to obtain higher economic benefits,it is suggested to adopt higher temperature rise rate(60-150℃/month).The more reliable kinetic parameters obtained can provide a basis for designing more reasonable scheme of in-situ heating conversion.

Effects of in-situ stress on the stability of a roadway excavated through a coal seam

Roadways excavated through a coal seam can exert an adverse effect on roadway stability. To investigate the effects of in-situ stress on roadway stability, numerical models were built and high horizontal stresses at varying orientations were applied. The results indicate that stress concentrations, roadway deformation and failure increase in magnitude and extent as the excavation angle with respect to the maximum horizontal stress increases. In addition, the stress adjacent to the coal-rock interface sharply varies in space and evolves with time; coal is much more vulnerable to deformation and failure than rock.The results provide insights into the layout of roadways excavated through a coal seam. Roadways should be designed parallel or at a narrow angle to the maximum horizontal stress. The concentrated stress at the top corner of the face-end should be reduced in advance, and the coal seam should be reinforced immediately after excavation.

Assessment of recoverable oil and gas resources by in-situ conversion of shale——Case study of extracting the Chang 7_(3) shale in the Ordos Basin

The purpose of this study is to investigate the entire evolution process of shales with various total organic contents(TOC)in order to build models for quantitative evaluation of oil and gas yields and establish methods for assessing recoverable oil and gas resources from in-situ conversion of organic matters in shale.Thermal simulation experiments under in-situ conversion conditions were conducted on Chang 7_(3) shales from the Ordos Basin in a semi-open system with large capacity.The results showed that TOC and R_(o) were the key factors affecting the in-situ transformation potential of shale.The remaining oil and gas yields increased linearly with TOC but inconsistently with R_(o).R_(o) ranged 0.75%—1.25%and 1.05%—2.3%,respectively,corresponding to the main oil generation stage and gas generation stage of shale in-situ transformation.Thus a model to evaluate the remaining oil/gas yield with TOC and R_(o) was obtained.The TOC of shale suitable for in-situ conversion should be greater than 6%,whereas its R_(o) should be less than 1.0%.Shales with 0.75%(R_(o))could obtain the best economic benefit.The results provided a theoretical basis and evaluation methodology for predicting the hydrocarbon resources from in-situ conversion of shale and for the identification of the optimum“sweet spots”.The assessment of the Chang 7_(3) shale in the Ordos Basin indicated that the recoverable oil and gas resources from in-situ conversion of organic matters in shale are substantial,with oil and gas resources reaching approximately 450×10^(8) t and 30×10^(12)m^(3),respectively,from an area of 4.27×10^(4) km^(2).

Numerical analysis on mechanical difference of sandstone under in-situ stress,pore pressure preserved environment at depth

Deep in-situ rock mechanics considers the influence of the in-situ environment on mechanical properties,differentiating it from traditional rock mechanics.To investigate the effect of in-situ stress,pore pressure preserved environment on the mechanical difference of sandstone,four tests are numerically modeled by COMSOL:conventional triaxial test,conventional pore pressure test,in-situ stress restoration and reconstruction test,and in-situ pore pressure-preserved test(not yet realized in the laboratory).The in-situ stress restoration parameter is introduced to characterize the recovery effect of in-situ stress on elastic modulus and heterogeneous distribution of sandstone at different depths.A random function and nonuniform pore pressure coefficient are employed to describe the non-uniform distribution of pore pressure in the in-situ environment.Numerical results are compared with existing experimental data to validate the models and calibrate the numerical parameters.By extracting mechanical parameters from numerical cores,the stress-strain curves of the four tests under different depths,in-situ stress and pore pressure are compared.The influence of non-uniform pore pressure coefficient and depth on the peak strength of sandstone is analyzed.The results show a strong linear relationship between the in-situ stress restoration parameter and depth,effectively characterizing the enhanced effect of stress restoration and reconstruction methods on the elastic modulus of conventional cores at different depths.The in-situ pore pressurepreserved test exhibits lower peak stress and peak strain compared to the other three tests,and sandstone subjected to non-uniform pore pressure is more prone to plastic damage and failure.Moreover,the influence of non-uniform pore pressure on peak strength gradually diminished with increasing depth.

Design and development of the deep-rock in-situ condition-preserved coring calibration platform

To systematically validate and calibrate the theory and technology of the deep in-situ conditionpreserved coring, the in-situ conditions at different depths should be simulated, and the full-size coring tests should be carried out in this simulated environment. Therefore, a deep-rock in-situ conditionpreserved coring calibration platform was designed and developed. The self-tightening sealing structure and the quick-disassembly structure were designed on the basis of an innovative segmented nonuniformdiameter structure, which was a breakthrough from the traditional high-pressure vessel frame and was verified by finite element simulation and actual testing under extreme working conditions, respectively.To simulate the actual deep in-situ environment with a temperature of 150℃ and pressure of 140 MPa for a large Φ450 mm×H1400 mm core, temperature and pressure control systems were designed by coupling, and a pre-embedded high-pressure-resistant temperature sensor was designed. Finally, highprecision assembly automation, complex movement coordination of the coring device with the platform,and rotary dynamic sealing were achieved by utilizing the combination of adaptive cabin body servo control and an adaptive mechanical structure in a limited space, laying a solid foundation for the calibration of in-situ condition-preserved coring.

In-situ remediation of deep petroleum-contaminated soil injection

A computational fluid dynamics(CFD)numerical simulation and field experiment were used to investigate optimal operating parameters of high-pressure jet grouting equipment and clarify the boundary law of the injection area in the remediation process.The response surface optimization design results show that the optimal injection pressure is 30 MPa,rotation speed is 23 r/min,commission speed is 30 cm/min,and the optimal injection diameter is 147.3 cm.Based on the CFD numerical simulation,the ratio of the injection core,turbulent zone,and seepage zone is approximately 1∶4∶2.The distribution law of jet core,turbulence zone and seepage zone at different cross-sections under 30 MPa operating conditions is as follows:The jet core radius is approximately 100 mm,the turbulence zone is mainly distributed at 100 to 500 mm,the seepage zone is mainly distributed at 500 to 700 mm,the seepage zone could be completed within 2 h,and the proportion of the three boundary zones in the injection zone is similar to that of the numerical simulation.This study provides theoretical parameters and practical reference for the remediation of deep pollution via in-situ chemical oxidation in the Loess Plateau soil environment.

Stability analysis of gravity anchor foundation of layered argillaceous sandstone under dry-wet cycles

To investigate the stability of gravity anchors of suspension bridges,in-situ tests of the vertical bearing capacity of the bedrock,shear resistance of the anchor-rock interface,shear resistance of the bedrock were conducted in a suspension bridge project.Under dry-wet cycles,the deterioration law of the mechanical properties of argillaceous sandstone was identified in laboratory tests:the elastic modulus,cohesion and friction of the argillaceous sandstone deteriorated significantly at first few dry-wet cycles and then declined slowly after 10 cycles,ultimately these three mechanical parameters were reduced to about 1/3,1/3,2/3 of the initial value respectively.Moreover,numerical simulation was used to restore in-situ shear tests and a good agreement was obtained.Base on the results of in-situ and laboratory tests,the stability of the gravity anchor foundation under natural conditions and drywet cycles was calculated and its failure modes were analyzed.The results demonstrated that the dry-wet cycles caused uneven settlement of the anchor foundation,resulting in more serious stress concentration in the substrate.The dry-wet cycles remarkably reduced the stability coefficient of the anchor foundation,whose failure mode shifted from overturning failure under natural conditions to sliding failure.When there was weak interlayer in the rock layer,the anti-sliding stability of the anchor foundation was affected drastically.

Simulation of Excess Pore Water Pressure During Deep Soil Mixing Columns Installing

Most of current studies of deep soil mixing (DSM) methods are focused on the soil strength improvement and soil treatment effectiveness. But the DSM installation leads to excess pore water pressure and soil disturbance, which will bring great harm to adjacent structures, such as shell tunnels and historic buildings. The procedure of excess pore water pressure buildup while large number DSM columns are installed is complicated. In order to find methods to predict and simulate the excess pore water pressure during DSM column installation, the complicated dissipation and buildup of excess pore water pressure through in-situ test are studied in this paper. In-situ test was conducted in soft clay near the Huangpu River in Shanghai. The pore water pressure was investigated by an automatic monitoring system. Test results indicate that the excess pore water pressure induced by one DSM column installation is composed of the compaction pressure and the reversing pressure. The empirical equations of excess pore water pressure dissipation and buildup were built by mathematical fitting methods. A compound method is proposed to simulate the excess pore water pressure due to DSM installation. Using this method to predict the excess pore water pressure in the situ test, results show a well agreement between the prediction and the measurements.

Mixing processes and patterns of fluids in alkane-CO_(2)-water systems under high temperature and high pressure——Microscopic visual physical thermal simulations and molecular dynamics simulations

Various types of geofluids exist in deep and ultra-deep layers in petroliferous basins.The geofluids are much more active under high-temperature and high-pressure(HTHP)conditions,but their properties are unclear.We simulated the mixing of different fluids in CH_(4)/C_(3)H_(8)/C_(6)H_(14)/C_(8)H_(18)-water systems and C_(6)H_(14)/C_(8)H_(18)-CO_(2)-H_(2)O systems at temperatures of 25℃ to 425℃ and pressures of 5 MPa to 105 MPa,using an in-situ micron quartz capillary tube thermal simulation system and molecular dynamics numerical simulation software.The mixing processes,patterns,and mechanisms of various fluids were analyzed at microscale under increasing temperature and pressure conditions.The results show that the miscibility of fluids in the different alkane-H_(2)O and alkane-CO_(2)-H_(2)O systems is not instantaneous,but the miscibility degree between different fluid phases increases as the temperature and pressure rise during the experiments.The physical thermal experiments(PTEs)show that the mixing process can be divided into three stages:initial miscibility,segmented dynamic miscibility,and complete miscibility.The molecular dynamics numerical simulations(MDNSs)indicate that the mixing process of fluids in the alkane-H_(2)O and alkane CO_(2)-H_(2)O systems can be divided into seven and eight stages,respectively.The carbon number affects the miscibility of alkanes and water,and the temperature and pressure required to reach the same miscibility stage with water increase with the carbon number(C_(3)H_(8),C_(6)H_(14),CH_(4),C_(8)H_(18)).CO_(2) has a critical bridge role in the miscibility of alkanes and water,and its presence significantly reduces the temperatures required to reach the initial,dynamic,and complete miscibility of alkanes and water.The results are of great significance for analyzing and understanding the miscibility of geofluids in deep and ultra-deep HTHP systems.

Efficient grain refinement of Al alloys induced by in-situ nanoparticles

Grain refinement of Al and Al alloys during the solidification remains a long-term challenge.Here we report a general strategy for efficient grain refinement through instantaneously in-situ forming abundant nanoparticles with concentrated size distribution.The Al-Zn-Mg-Cu alloy and commercial purity Al are respectively refined by using this proposed strategy and the traditional refinement method(adding AlTi-B alloy into Al melt),and the experimental and numerical results indicate that the TiBnanoparticles instantaneously form in the Al melt during the in-situ refinement process.Compared to the traditional refinement treatment method,the new approach not only can promote the high number density of the nucleation sites and narrow their size distribution,but also tremendously attenuates the agglomeration and settlement of the heterogeneous nucleation sites.It exhibits much better grain refinement ability and inhibits the decline of the grain refinement efficiency.This approach may have broad application prospect not only for the grain refinement of Al and Al alloys,but also for the grain refinement of other alloy systems,e.g.steel,magnesium,copper,and so on.

Numerical study on the effect of in-situ stress on smoothwall blasting in deep tunnelling

In the present study,a numerical model is first calibrated against the crack networks and pressure attenuation data in laboratory blasting test.Then,based on the calibrated numerical model,two-hole plane models are developed and used to perform a series of sim-ulations of smoothwall blasting in deep tunnelling subjected to in-situ stress.The evolutions of rock fracture and excavation damage zone in the roof/floor and sidewalls under different far-field hydrostatic pressure and anisotropic in-situ stress conditions are numerically investigated.The findings in numerical modelling are also analytically interpreted with the stress distributions around the designed tunnel perimeter and perimeter borehole.The numerical and analytical results show that the variations of rock cracking and excavation dam-aged zone induced by smoothwall blasting with in-situ stress are mainly attributed to the high tangential compressive stress concentration around the remaining rock after inner primary blasts and the tensile stress acting on the wall of perimeter hole,which control the crack propagation and initiation respectively.At last,the implications of findings for practical smoothwall blasting in deep tunnelling are discussed.

Room temperature in-situ preparation of hydrazine-linked covalent organic frameworks coated capillaries for separation and determination of polycyclic aromatic hydrocarbons

Covalent organic frameworks(COFs)have been increasingly used in capillary electrochromatography due to their excellent characteristics.In this work,hydrazine-linked TFPB-DHzDS(TFPB:1,3,5-tris(4-formylphenyl)benzene;DHzDS:2,5-bis(3-(ethylthio)propoxy)terephthalohydrazide)was first synthesized by a simpler and easier method at room temperature and introduced into capillary electrochromatography as coating material.The TFPB-DHzDS coated capillaries were prepared by an in-situ growth process at room temperature.After optimizing the coating concentration and experimental conditions of capillary electrochromatography,baseline separation of two groups of polycyclic aromatic hydrocarbons was achieved based on the TFPB-DHzDS coated capillary.And the established method was used successfully to determine PAHs in natural water and soil samples.The spiked recoveries of polycyclic aromatic hydrocarbons in these samples ranged from 90.01%to 111.0%,indicating that the method is reliable and could detect polycyclic aromatic hydrocarbons in natural samples.Finally,molecular simulation was applied to study and visualize the interaction between the analytes and coating COF materials to investigate the molecular level separation mechanism further.

A three-dimensional numerical study on the stability of layered rock spillway tunnels in alpine canyon areas

Rock masses in alpine canyon areas exhibit strong heterogeneity,discontinuity,and are subject to strong tectonic effects and stress unloading,leading to extremely complex distribution of in-situ stress.In addition,the occurrence of layered rock masses makes it more complex,with obvious anisotropic mechanical properties.This study proposes a comprehensive method for evaluating the stability of layered rock spillway tunnels in a hydropower station in an alpine canyon.First,the failure criterion and mechanical model of layered rock masses considering the anisotropy induced by the bedding plane and the true triaxial stress regime were established;an inversion theory and calculation procedure for in-situ stress in alpine canyon areas were then introduced.Finally,by using a self-developed numerical tool,i.e.CASRock,the stability of the layered rock spillway tunnel in a hydropower station was numerically analyzed.The results show that,affected by geological structure and stratigraphic lithology,there is significant differentiation in the in-situ stress in alpine canyons,with horizontal tectonic stress as the main factor.The occurrence of layered rock masses in the region has a significant impact on the stability of surrounding rock,and the angle between the bedding strike and the tunnel axis as well as the bedding dip both exert a significant influence on the failure characteristics of the surrounding rock.

PROCESS CONTROL OF THE SAND WAVE MIGRATION IN BEIBU GULF OF THE SOUTH CHINA SEA

Based on the environment characteristics of the Beibu Gulf of South China Sea, a quasi-three-dimensional physical model is built. By coupling the bottom boundary layer with the two-dimensional tidal current field near the seabed surface, the quasi-three-dimensional hydrodynamic numerical simulation is carried out. The sand wave migration process is dealt with by coupling the hydrodynamic model with the sediment transport model. The computational results are shown to be in good agreement with the observed data, which indicates that the quasi-three-dimensional physical model can be used to simulate the migration process for small scale sand waves. Then, based on measured data, the evolution of the sand wave migration is investigated. An effective formula is developed to predict the migration rate, in which not only the effects of the environment but also the features of sand waves are considered.

Impact of hydraulic perforation on fracture initiation and propagation in shale rocks

To enhance the oil and gas recovery rate, hydraulic fracturing techniques have been widely adopted for stimulation of low-permeability reservoirs. Pioneering work indicates that hydraulic perforation and layout could significantly affect fracture initiation and propagation in low-permeability reservoir rocks subjected to complex in-situ stresses. This paper reports on a novel numerical method that incorporates fracture mechanics principles and the numerical tools FRANC3D and ANSYS to investigate the three-dimensional initiation and propagation behavior of hydro-fracturing cracks in shale rock. Considering the transverse isotropic property of shale rocks, the mechanical parameters of reservoir rocks attained from laboratory tests were adopted in the simulation. The influence of perforation layouts on the 3D initiation of hydro-fracturing fractures in reservoir rocks under geo-stresses was quantitatively illuminated. The propagation and growth of fractures in three dimensions in different perforating azimuth values were illustrated. The results indicate that： 1） the optimal perforation direction should be parallel to the maximum horizontal principal stress, 2） the crack plane gradually turns toward the direction of the maximum horizontal principal stress when they are not in parallel, 3） compared with the linear and symmetric pattern, the staggered perforation is the optimal one, 4） the proper perforation density is four to six holes per meter, 5） the optimal perforation diameter in this model is 30 ram, and 6） the influence of the perforation depth on the fracture initiation pressure is low.

Numerical evaluation of the effectiveness of NO_2 and N_2O_5 generation during the NO ozonation process

Wet scrubbing combined with ozone oxidation has become a promising technology for simultaneous removal of SO_2 and NO_x in exhaust gas.In this paper,a new 20-species,76-step detailed kinetic mechanism was proposed between O_3 and NO_x.The concentration of N_2O_5 was measured using an in-situ IR spectrometer.The numerical evaluation results kept good pace with both the public experiment results and our experiment results.Key reaction parameters for the generation of NO_2 and N_2O_5 during the NO ozonation process were investigated by a numerical simulation method.The effect of temperature on producing NO_2 was found to be negligible.To produce NO_2,the optimal residence time was1.25 sec and the molar ratio of O_3/NO about 1.For the generation of N_2O_5,the residence time should be about 8 sec while the temperature of the exhaust gas should be strictly controlled and the molar ratio of O_3/NO about 1.75.This study provided detailed investigations on the reaction parameters of ozonation of NO_x by a numerical simulation method,and the results obtained should be helpful for the design and optimization of ozone oxidation combined with the wet flue gas desulfurization methods（WFGD） method for the removal of NO_x.