Relative Impacts of Sea Ice Loss and Atmospheric Internal Variability on the Winter Arctic to East Asian Surface Air Temperature Based on Large-Ensemble Simulations with NorESM2

To quantify the relative contributions of Arctic sea ice and unforced atmospheric internal variability to the “warm Arctic, cold East Asia”(WACE) teleconnection, this study analyses three sets of large-ensemble simulations carried out by the Norwegian Earth System Model with a coupled atmosphere–land surface model, forced by seasonal sea ice conditions from preindustrial, present-day, and future periods. Each ensemble member within the same set uses the same forcing but with small perturbations to the atmospheric initial state. Hence, the difference between the present-day(or future) ensemble mean and the preindustrial ensemble mean provides the ice-loss-induced response, while the difference of the individual members within the present-day(or future) set is the effect of atmospheric internal variability. Results indicate that both present-day and future sea ice loss can force a negative phase of the Arctic Oscillation with a WACE pattern in winter. The magnitude of ice-induced Arctic warming is over four(ten) times larger than the ice-induced East Asian cooling in the present-day(future) experiment;the latter having a magnitude that is about 30% of the observed cooling. Sea ice loss contributes about 60%(80%) to the Arctic winter warming in the present-day(future) experiment. Atmospheric internal variability can also induce a WACE pattern with comparable magnitudes between the Arctic and East Asia. Ice-lossinduced East Asian cooling can easily be masked by atmospheric internal variability effects because random atmospheric internal variability may induce a larger magnitude warming. The observed WACE pattern occurs as a result of both Arctic sea ice loss and atmospheric internal variability, with the former dominating Arctic warming and the latter dominating East Asian cooling.

CAS-ESM2.0 Successfully Reproduces Historical Atmospheric CO_(2) in a Coupled Carbon−Climate Simulation

The atmospheric carbon dioxide(CO_(2))concentration has been increasing rapidly since the Industrial Revolution,which has led to unequivocal global warming and crucial environmental change.It is extremely important to investigate the interactions among atmospheric CO_(2),the physical climate system,and the carbon cycle of the underlying surface for a better understanding of the Earth system.Earth system models are widely used to investigate these interactions via coupled carbon-climate simulations.The Chinese Academy of Sciences Earth System Model version 2(CAS-ESM2.0)has successfully fixed a two-way coupling of atmospheric CO_(2)with the climate and carbon cycle on land and in the ocean.Using CAS-ESM2.0,we conducted a coupled carbon-climate simulation by following the CMIP6 proposal of a historical emissions-driven experiment.This paper examines the modeled CO_(2)by comparison with observed CO_(2)at the sites of Mauna Loa and Barrow,and the Greenhouse Gases Observing Satellite(GOSAT)CO_(2)product.The results showed that CAS-ESM2.0 agrees very well with observations in reproducing the increasing trend of annual CO_(2)during the period 1850-2014,and in capturing the seasonal cycle of CO_(2)at the two baseline sites,as well as over northern high latitudes.These agreements illustrate a good ability of CAS-ESM2.0 in simulating carbon-climate interactions,even though uncertainties remain in the processes involved.This paper reports an important stage of the development of CAS-ESM with the coupling of carbon and climate,which will provide significant scientific support for climate research and China’s goal of carbon neutrality.

Basalt Petrology, Water Chemistry, and Their Impact on the CO_(2) Mineralization Simulation at Leizhou Peninsula Sites, Southern China

Mineral carbonation, which precipitates dissolved carbon dioxide(CO_(2)) as carbonate minerals in basaltic groundwater environments, is a potential technique for negative emissions. The Leizhou Peninsula in southwest Guangdong province has extensive basalt, indicating a promising potential for CO_(2) storage through rapid mineralization. However, understanding of the basic geological setting, potential, and mechanisms of CO_(2) mineralization in the basalts of the Leizhou Peninsula is still limited. The mineralization processes associated with CO_(2)storage at two candidate sites in the area are investigated in this paper: Yongshi Farm and Tianyang Basin(of the dried maar lake). Petrography,rock geochemistry, basalt petrophysical properties, and groundwater hydrochemistry analyses are included in the study. Numerical simulation is used to examine the reaction process and its effects. The results show that basalts in the study areas mainly comprise plagioclase, pyroxene, and Fe–Ti oxides, revealing a total volume fraction exceeding 85%. Additionally, small amounts of quartz and fayalite are available, with volume fractions of 5.1% and 1.0%, respectively. The basalts are rich in divalent metal cations, which can form carbonate minerals, with an average of approximately 6.2 moles of metal cations per 1 kg of rock. The groundwater samples have a pH of 7.5–8.2 and are dominated by the Mg–Ca–HCO3 type. The basalts demonstrate a porosity range of 10.9% to 28.8%, with over 70% of interconnected pores. A 20-year geochemical simulation revealed that CO_(2) injection dissolves primary minerals, including anorthite, albite, and diopside, while CO_(2)mineralization dissolves precipitation secondary minerals, such as calcite, siderite, and dolomite. Furthermore, a substantial rise in pH from 7.6to 10.6 is observed in the vicinity of the injected well, accompanied by a slight reduction in porosity from 20% to 19.8%. Additionally, 36.8% of the injected CO_(2) underwent complete mineralization within five years, revealing an increasing percentage of 66.1% if the experimental period is extended to 20 years. The presence of abundant divalent metal cations in basalts and water-bearing permeable rocks in the Leizhou Peninsula supports the potential for mineral carbonation in basalts, as indicated by the geochemical simulation results. Additional research is necessary to identify the factors that influence the CO_(2) mineralization, storage, and sensitivity analysis of basalt in the Leizhou Peninsula.

Numerical Simulation of Asphaltene Precipitation and Deposition during Natural Gas and CO_(2) Injection

Asphaltene deposition is a significant problem during gas injection processes,as it can block the porous medium,the wellbore,and the involved facilities,significantly impacting reservoir productivity and ultimate oil recovery.Only a few studies have investigated the numerical modeling of this potential effect in porous media.This study focuses on asphaltene deposition due to natural gas and CO_(2) injection.Predictions of the effect of gas injection on asphaltene deposition behavior have been made using a 3D numerical simulation model.The results indicate that the injection of natural gas exacerbates asphaltene deposition,leading to a significant reduction in permeability near the injection well and throughout the reservoir.This reduction in permeability strongly affects the ability of gas toflow through the reservoir,resulting in an improvement of the displacement front.The displacement effi-ciency of the injection gas process increases by up to 1.40%when gas is injected at 5500 psi,compared to the scenario where the asphaltene model is not considered.CO_(2) injection leads to a miscible process with crude oil,extracting light and intermediate components,which intensifies asphaltene precipitation and increases the viscosity of the remaining crude oil,ultimately reducing the recovery rate.

Investigation of gravity influence on EOR and CO_(2) geological storage based on pore-scale simulation

Gravity assistance is a critical factor influencing CO_(2)-Oil mixing and miscible flow during EOR and CO_(2)geological storage.Based on the Navier-Stokes equation,component mass conservation equation,and fluid property-composition relationship,a mathematical model for pore-scale CO_(2) injection in oilsaturated porous media was developed in this study.The model can reflect the effects of gravity assistance,component diffusion,fluid density variation,and velocity change on EOR and CO_(2) storage.For nonhomogeneous porous media,the gravity influence and large density difference help to minimize the velocity difference between the main flow path and the surrounding area,thus improving the oil recovery and CO_(2) storage.Large CO_(2) injection angles and oil-CO_(2) density differences can increase the oil recovery by 22.6% and 4.2%,respectively,and increase CO_(2) storage by 37.9% and 4.7%,respectively.Component diffusion facilitates the transportation of the oil components from the low-velocity region to the main flow path,thereby reducing the oil/CO_(2) concentration difference within the porous media.Component diffusion can increase oil recovery and CO_(2) storage by 5.7% and 6.9%,respectively.In addition,combined with the component diffusion,a low CO_(2) injection rate creates a more uniform spatial distribution of the oil/CO_(2) component,resulting in increases of 9.5% oil recovery and 15.7% CO_(2) storage,respectively.This study provides theoretical support for improving the geological CO_(2) storage and EOR processes.

A review of interaction mechanisms and microscopic simulation methods for CO_(2)-water-rock system

This work systematically reviews the complex mechanisms of CO_(2)-water-rock interactions,microscopic simulations of reactive transport(dissolution,precipitation and precipitate migration)in porous media,and microscopic simulations of CO_(2)-water-rock system.The work points out the key issues in current research and provides suggestions for future research.After injection of CO_(2) into underground reservoirs,not only conventional pressure-driven flow and mass transfer processes occur,but also special physicochemical phenomena like dissolution,precipitation,and precipitate migration.The coupling of these processes causes complex changes in permeability and porosity parameters of the porous media.Pore-scale microscopic flow simulations can provide detailed information within the three-dimensional pore and throat space and explicitly observe changes in the fluid-solid interfaces of porous media during reactions.At present,the research has limitations in the decoupling of complex mechanisms,characterization of differential multi-mineral reactions,precipitation generation mechanisms and characterization(crystal nucleation and mineral detachment),simulation methods for precipitation-fluid interaction,and coupling mechanisms of multiple physicochemical processes.In future studies,it is essential to innovate experimental methods to decouple“dissolution-precipitation-precipitate migration”processes,improve the accuracy of experimental testing of minerals geochemical reaction-related parameters,build reliable characterization of various precipitation types,establish precipitation-fluid interaction simulation methods,coordinate the boundary conditions of different physicochemical processes,and,finally,achieve coupled flow simulation of“dissolution-precipitation-precipitate migration”within CO_(2)-water-rock systems.

Optical properties of La_(2)O_(3) and HfO_(2) for radiative cooling via multiscale simulations

Radiative cooling materials have gained prominence as a zero-energy solution for mitigating global warming.However,a comprehensive understanding of the atomic-scale optical properties and macroscopic optical performance of radiative cooling materials remains elusive,limiting insight into the underlying physics of their optical response and cooling efficacy.La_(2)O_(3)and HfO_(2),which represent rare earth and third/fourth subgroup inorganic oxides,respectively,show promise for radiative cooling applications.In this study,we used multiscale simulations to investigate the optical properties of La_(2)O_(3)and HfO_(2)across a broad spectrum.First-principles calculations revealed their dielectric functions and intrinsic refractive indices,and the results indicated that the slightly smaller bandgap of La_(2)O_(3)compared to HfO_(2)induces a higher refractive index in the solar band.Additionally,three-phonon scattering was found to provide more accurate infrared optical properties than two-phonon scattering,which enhanced the emissivity in the sky window.Monte Carlo simulations were also used to determine the macroscopic optical properties of La_(2)O_(3)and HfO_(2)coatings.Based on the simulated results,we identified that the particle size and particle volume fraction play a dominant role in the optical properties.Our findings underscore the potential of La_(2)O_(3)and HfO_(2)nanocomposites for environment-friendly cooling and offer a new approach for high-throughput screening of optical materials through multiscale simulations.

Formation mechanism and crystal simulation of Na_2O-doped calcium aluminate compounds

Calcium aluminate clinkers doped with Na2O were synthesized using analytically pure reagents CaCO3, Al2O3, SiO2 and Na2CO3. The effects of Na2O-doping on the formation mechanism of calcium aluminate compounds and the crystal property of 12CaO·7Al2O3 (C12A7) cell were studied. The results show that the minerals containing Na2O mainly include 2Na2O·3CaO·5Al2O3 and Na2O·Al2O3, when the Na2O content in clinkers is less than 4.26% (mass fraction). The rest of Na2O is mainly doped in 12CaO·7Al2O3, which results in the decrease of the crystallinity of 12CaO·7Al2O3. The crystallinity of 2Na2O·3CaO·5Al2O3 is also inversely proportional to the Na2O content in clinkers. The formation processes of 2Na2O·3CaO·5Al2O3 and 12CaO·7Al2O3 can be divided into two ways, which are the direct reactions of raw materials and the transformation of CaO·Al2O3, respectively. The simulation shows that the covalency of O-Na bond in Na2O-doped 12CaO·7Al2O3 cell is weaker than those of O-Ca and O-Al bonds. The free energy of the unit cell increases because of Na2O doping, which results in the improvement of chemical activity of 12CaO·7Al2O3. The leaching efficiency of Al2O3 in clinker is improved from 34.81% to 88.17% when the Na2O content in clinkers increases from 0 to 4.26%.

注CO_(2)驱提高致密砂岩气藏采收率实验及数值模拟

针对致密砂岩气藏地层压力下降快、稳产难度大、采收率较低等问题,利用长岩心驱替实验、可视化实验和数值模拟方法,从提高气藏采收率方面对注CO_(2)的驱替特征和扩散特征进行了研究。结果表明:地层中的超临界CO_(2)不是与天然气大范围地互溶,而是“部分互溶”;由于重力分异作用,CO_(2)在沉降过程中驱替天然气向上运移并被采出;CO_(2)驱替天然气可提高采收率17.8百分点。研究成果可为下一步致密砂岩气藏CO_(2)驱矿场试验提供理论依据。

盐水层CO_(2)封存稳定性评价指标建立及储层优选

盐水层中CO_(2)封存被认为是缓解温室效应的主要途径,而CO_(2)长期封存的稳定性受储层和注入参数的影响,目前评价盐水层CO_(2)长期封存稳定性仍缺乏统一的指标.为实现盐水层CO_(2)长期封存稳定性评价,优选能够实现CO_(2)稳定封存的储层和注气强度,提出将封存量和转化速率(或封存速率)作为CO_(2)长期封存稳定性的因素指标.考虑主控因素,基于数值模拟和非线性多元回归分析,结合变异系数法,建立CO_(2)不同封存方式长期封存稳定性的评价指标.利用专家赋权法确定构造、残余、溶解和矿化4种封存方式的权重,建立综合考虑4种封存方式的CO_(2)长期封存稳定性的评价指标,并对CO_(2)封存评价指标划分为较稳定、稳定和不稳定3个等级,实现矿场特定储层及注气强度下的CO_(2)封存稳定性评价.同时,建立储层参数及注入参数影响下的CO_(2)长期封存稳定性评价指标图版,对优选适合CO_(2)长期稳定封存的盐水层具有参考价值.

Mathematical modeling and simulations of stress mitigation by coating polycrystalline particles in lithium-ion batteries

A chemo-mechanical model is developed to investigate the effects on the stress development of the coating of polycrystalline Ni-rich LiNixMnyCo_(z)O_(2)(x≥0.8)(NMC)particles with poly(3,4-ethylenedioxythiophene)(PEDOT).The simulation results show that the coating of primary NMC particles significantly reduces the stress generation by efficiently accommodating the volume change associated with the lithium diffusion,and the coating layer plays roles both as a cushion against the volume change and a channel for the lithium transport,promoting the lithium distribution across the secondary particles more homogeneously.Besides,the lower stiffness,higher ionic conductivity,and larger thickness of the coating layer improve the stress mitigation.This paper provides a mathematical framework for calculating the chemo-mechanical responses of anisotropic electrode materials and fundamental insights into how the coating of NMC active particles mitigates stress levels.

Numerical simulation of enhancing coalbed methane recovery by injecting CO_2 with heat injection

The technology used to enhance coalbed methane(CBM) recovery by injecting CO_2(CO_2-ECBM) with heat, combining heat injection with CO_2 injection, is still in its infancy; therefore, theoretical studies of this CO_2-ECBM technology should be perused. First, the coupling equations of the di usion–adsorption–seepage–heat transfer fields of gas are established. The displacement processes under di erent pressures and temperatures are simulated by COMSOL. Finally, the displacement effects, a comparison of the CO_2 storage capacity with the CH_4 output and the e ective influencing radius of CO_2 injection are analyzed and discussed. The results show that(1) the displacement pressure and temperature are two key factors influencing the CH_4 output and the CO_2 storage capacity, and the increase in the CO_2 storage capacity is more sensitive to temperature and pressure than the CH_4 output.(2) The gas flow direction is from the injection hole to the discharge hole during the displacement process, and the regions with high velocity are concentrated at the injection hole and the discharge hole.(3) A reduction in the CH_4 concentration and an increase in the CO_2 concentration are obvious during the displacement process.(4) The e ective influencing radius of injecting CO_2 with heat increases with the increase in time and pressure. The relationship between the e ective influencing radius and the injection time of CO_2 has a power exponential function, and there is a linear relationship between the functional coe cient and the injection pressure of CO_2. This numerical simulation study on enhancing CBM recovery by injecting CO_2 with heat can further promote the implementation of CO_2-ECBM project in deep coal seams.

A Parallel Discrete Event Simulation Engine for the Low-Earth-Orbit Satellite Constellation Networks

Low-Earth-Orbit satellite constellation networks(LEO-SCN)can provide low-cost,largescale,flexible coverage wireless communication services.High dynamics and large topological sizes characterize LEO-SCN.Protocol development and application testing of LEO-SCN are challenging to carry out in a natural environment.Simulation platforms are a more effective means of technology demonstration.Currently available simulators have a single function and limited simulation scale.There needs to be a simulator for full-featured simulation.In this paper,we apply the parallel discrete-event simulation technique to the simulation of LEO-SCN to support large-scale complex system simulation at the packet level.To solve the problem that single-process programs cannot cope with complex simulations containing numerous entities,we propose a parallel mechanism and algorithms LP-NM and LP-YAWNS for synchronization.In the experiment,we use ns-3 to verify the acceleration ratio and efficiency of the above algorithms.The results show that our proposed mechanism can provide parallel simulation engine support for the LEO-SCN.

Implication of Geochemical Simulation for CO2 Storage Using Data of York Reservoir

In this paper, we build up a three-dimensional model for CO2 storage in the deep reservoir. And this paper gives the mathematical formalism of combined geochemical and multi-phase flow. The results give us the information about geochemical changing caused by CO2 injection into aqueous, the dissolution or precipitation of reservoir minerals caused by aqueous components change, the change of water density, also the differences between this model and the simulation model without considering geochemical. The basic data for simulation is from York Reservoir.

Simulation and optimization of an industrial gas condensate stabilization unit to modify LPG and NGL production with minimizing CO2 emission to the environment

In the present study, a great effort was made to improve the performance of an industrial liquefied petroleum gas(LPG) and natural gas liquid(NGL) production unit in one of the major gas refinery located at Pars special economic zone in Iran. To demonstrate and obtain the optimal condition, the unit was simulated by using a steady-state flowsheet simulator, i.e. Aspen Plus, under different operational conditions. According to the simulation results,the unit was not operational under its optimal conditions due to some defects in the cooling system at top stage of the debutanizer tower(DBT) during hot and humid seasons. Additionally, the vapor pressure of produced LPG and accordingly the amount of its flaring were decreased by reducing the temperature of debutanizer tower at top stages. In the optimization section, the DBT condenser and reboiler heat duty, temperature, and pressure were regulated as adjustable parameters. The simulation results demonstrated that by applying the optimum suggestion in the hot months, the reflux stream temperature was reached about 55 ℃ which caused an efficient increment in LPG production(about 4%) with adjusting the propane component in LPG, based on the standard range as the plant criteria. Moreover, after applying modifications, about 750 t of LPG product was saved from flaring during five hot months of the year, which resulted in 360000 USD extra annual income for the company.Finally, from environmental point of view, this optimization caused to reduce 81 t of CO_2 emission to the environment. Therefore, the current investigation must be introduced as a friendly environmentally process.

Two Improvements on Numerical Simulation of 2-DOF Vortex-Induced Vibration with Low Mass Ratio

Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vortex-induced vibration with low mass ratio, the accuracy is not satisfactory, especially for the maximum amplitudes. In Jauvtis and Williamson＇s work, the maximum amplitude of the cylinder with low mass ratio m＊=2.6 can reach as large as 1.5D to be called as the ＂super-upper branch＂, but from current literatures, few simulation results can achieve such value, even fail to capture the upper branch. Besides, it is found that the amplitude decays too fast in the lower branch with the RANS-based turbulence model. The reason is likely to be the defects of the turbulence model itself in the prediction of unsteady separated flows as well as the unreasonable setting of the numerical simulation parameters. Aiming at above issues, a modified turbulence model is proposed in this paper, and the effect of the acceleration of flow field on the response of vortex-induced vibration is studied based on OpenFOAM. By analyzing the responses of amplitude, phase and trajectory, frequency and vortex mode, it is proved that the vortex-induced vibration can be predicted accurately with the modified turbulence model under appropriate flow field acceleration.

EFFICIENT STRIP-MODE SAR RAW SIGNAL SIMULATION OF MIXED TARGETS BASED ON ACCURATE 2-D SPECTRUM

An accurate and efficient Synthetic Aperture Radar(SAR)raw data generator is of considerable value for testing system parameters and verifying imaging algorithms.Nevertheless,the existing simulator cannot exactly handle the case of the fast moving targets in high squint geometry.As for the issue,the analytical expression for the two Dimensional(2-D)signal spectrum of moving targets is derived and a fast raw echo simulation method is proposed in this study.The proposed simulator can accommodate the moving targets in the high squint geometry,whose processing steps of the simulation are given in detail and its computational complexity is analyzed.The simulation data for static and moving targets are processed and analyzed,and the results are given to validate the effectiveness of the proposed approach.

Using TOUGH2 numerical simulation to analyse the geothermal formation in Guide basin,China

The Guide sedimentary basin is located in the northeastern part of Qinghai-Xizang Plateau,which is rich in geothermal resources.However,exploitation of the geothermal resources has so far been limited,because of limited understanding of the resources quantity and storage gained from scientific researches.In this study,using a typical cross section across the basin and taking into account its geothermal and geological conditions,a new waterheat coupled model was built and associated modelling was done by the software TOUGH2.During modelling process,the accuracy and applicability of the model was confirmed through the calibration of relevant parameters for modelling the heat and water transport and the formation of geothermal reservoir across the basin,with particular focus on the Neogene geothermal field.Results show that the groundwater that flows from the basin margins to the center is heated by the Neogene and Paleogene sedimentary rocks with high geothermal gradients.Since the east-west extending fault F1 is conductive,it acts as preferential flow paths which on one hand provide additional and rapid flows to the thermal reservoir;and on the other hand,cool down the thermal water to a certain extent due to the infiltration of shallower water sources in the vicinity of the fault.Furthermore,the estimated geothermal resources quantity is close to that of previous studies.In comparison with the Paleogene rock formations,the Neogene geothermal reservoir shows a better nature in terms of water content,aquifer permeability and resources exploitability,although the resource quantity of the Paleogene reservoir is considerable.

Simulation of gas-solid flow characteristics of the circulating fluidized bed boiler under pure-oxygen combustion conditions

Under the pressure of carbon neutrality,many carbon capture,utilization and storage technologies have witnessed rapid development in the recent years,including oxy-fuel combustion(OFC)technology.However,the conventional OFC technology usually depends on the flue gas recirculation system,which faces significant investment,high energy consumption,and potential low-temperature corrosion problem.Considering these deficiencies,the direct utilization of pure oxygen to achieve particle fluidization and fuel combustion may reduce the overall energy consumption and CO_(2)-capture costs.In this paper,the fundamental structure of a self-designed 130 t·h^(-1) pure-oxygen combustion circulating fluidized bed(CFB)boiler was provided,and the computational particle fluid dynamics method was used to analyze the gas-solid flow characteristics of this new-concept boiler under different working conditions.The results indicate that through the careful selection of design or operational parameters,such as average bed-material size and fluidization velocity,the pure-oxygen combustion CFB system can maintain the ideal fluidization state,namely significant internal and external particle circulation.Besides,the contraction section of the boiler leads to the particle backflow in the lower furnace,resulting in the particle suspension concentration near the wall region being higher than that in the center region.Conversely,the upper furnace still retains the classic core-annulus flow structure.In addition to increasing solid circulation rate by reducing the average bed-material size,altering primary gas ratio and bed inventory can also exert varying degrees of influence on the gas-solid flow characteristics of the pure-oxygen combustion CFB boiler.

Geomechanics model test and numerical simulation of 2G-NPR bolt support effect in an active fault tunnel

Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surrounding the cross-fault tunnel with poor geological conditions has been chosen for the study.The support capacity of 2^(nd) Generation-Negative Poisson’s Ratio(2G-NPR)bolt in an active fault tunnel has been carried out on the basis of relevant results obtained from the geomechanical model test and numerical investigations of failure model for existing unsupported fault tunnel.The investigation shows that surrounding rock of the tunnel is prone to shear deformation and crack formation along the fault,as a result,the rock mass on the upper part of the fault slips as a whole.Furthermore,small-scale deformation and loss of blocks are observed around the tunnel;however,the 2G-NPR bolt support is found to be helpful in keeping the overall tunnel intact without any damage and instability.Due to the blocking effect of fault,the stress of the surrounding rock on the upper and lower parts of the fault is significantly different,and the stress at the left shoulder of the tunnel is greater than that at the right shoulder.The asymmetrical arrangement of 2G-NPR bolts can effectively control the asymmetric deformation and instability of the surrounding rock.The present numerical scheme is in good agreement with the model test results,and can reasonably reflect the stress and displacement characteristics of the surrounding rock of the tunnel.In comparison to unsupported and ordinary PR(Poisson’s Ratio)bolt support,2G-NPR bolt can effectively limit the fault slip and control the stability of the surrounding rock of the fault tunnel.The research findings may serve as a guideline for the use of 2G-NPR bolts in fault tunnel support engineering.