Red-bed mudstone, prevalent in southwest China, poses a formidable challenge due to its hydrophilic clay minerals, resulting in expansion, deformation, and cracking upon exposure to moisture. This study addresses upli...Red-bed mudstone, prevalent in southwest China, poses a formidable challenge due to its hydrophilic clay minerals, resulting in expansion, deformation, and cracking upon exposure to moisture. This study addresses uplift deformation disasters in high-speed railways by employing a moisture diffusion-deformation-fracture coupling model based on the finite-discrete element method(FDEM). The model integrates the influence of cracks on moisture diffusion. The investigation into various excavation depths reveals a direct correlation between depth and the formation of tensile cracks at the bottom of the railway cutting. These cracks expedite moisture migration, significantly impacting the temporal and spatial evolution of the moisture field. Additionally, crack expansion dominates hygroscopic deformation, with the lateral coordinate of the crack zone determining peak vertical displacement. Furthermore, key factors influencing deformation in railway cuttings, including the swelling factor and initial moisture content at the bottom of the cutting, are explored. The number of tensile and shear cracks increases with greater excavation depth, particularly concerning shear cracks. Higher swelling factors and initial moisture contents result in an increased total number of cracks, predominantly shear cracks. Numerical calculations provide valuable insights, offering a scientific foundation and directional guidance for the precise prevention, control, prediction, and comprehensive treatment of mudstone-related issues in high-speed railways.展开更多
In this study,the axial swelling strain of red-bed mudstone under different vertical stresses are measured by swell-under-load method,and the microstructure of mudstone after hygroscopic swelling is studied by mercury...In this study,the axial swelling strain of red-bed mudstone under different vertical stresses are measured by swell-under-load method,and the microstructure of mudstone after hygroscopic swelling is studied by mercury intrusion porosimetry(MIP).The weakening coefficient and Weibull distribution function are introduced into the coupling model of mudstone moisture diffusion-swelling deformation-fracture based on finite-discrete element method(FDEM).The weakening effect of moisture on mudstone's mechanical parameters,as well as the heterogeneity of swelling deformation and stress distribution,is considered.The microcrack behavior and energy evolution of mudstone during hygroscopic swelling deformation under different vertical stresses are studied.The results show that the axial swelling strain of mudstone decreases with increase of the vertical stress.At low vertical stresses,moisture absorption in mudstone leads to formation of cracks caused by hydration-induced expansion.Under high vertical stresses,a muddy sealing zone forms on the mudstone surface,preventing further water infiltration.The simulation results of mudstone swelling deformation also demonstrate that it involves both swelling of the mudstone matrix and swelling caused by crack expansion.Notably,crack expansion plays a dominant role in mudstone swelling.With increasing vertical stress,the cracks in mudstone change from tensile cracks to shear cracks,resulting in a significant reduction in the total number of cracks.While the evolution of mudstone kinetic energy shows similarities under different vertical stresses,the evolution of strain energy varies significantly due to the presence of different types of cracks in the mudstone.The findings provide a theoretical basis for understanding the hygroscopic swelling deformation mechanism of red-bed mudstone at various depths.展开更多
The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(R...The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.展开更多
Uniaxial compression tests and cyclic loading acoustic emission tests were conducted on 20%,40%,60%,80%,dry and saturated muddy sandstone by using a creep impact loading system to investigate the mechanical properties...Uniaxial compression tests and cyclic loading acoustic emission tests were conducted on 20%,40%,60%,80%,dry and saturated muddy sandstone by using a creep impact loading system to investigate the mechanical properties and acoustic emission characteristics of soft rocks with different water contents under dynamic disturbance.The mechanical properties and acoustic emission characteristics of muddy sandstones at different water contents were analysed.Results of experimental studies show that water is a key factor in the mechanical properties of rocks,softening them,increasing their porosity,reducing their brittleness and increasing their plasticity.Under uniaxial compression,the macroscopic damage characteristics of the muddy sandstone change from mono-bevel shear damage and‘X’type conjugate bevel shear damage to a roadway bottom-drum type damage as the water content increases.Dynamic perturbation has a strengthening effect on the mechanical properties of samples with 60%and less water content,and a weakening effect on samples with 80%and more water content,but the weakening effect is not obvious.Macroscopic damage characteristics of dry samples remain unchanged,water samples from shear damage and tensile–shear composite damage gradually transformed into cleavage damage,until saturation transformation monoclinic shear damage.The evolution of acoustic emission energy and event number is mainly divided into four stages:loading stage(Ⅰ),dynamic loading stage(Ⅱ),yield failure stage(Ⅲ),and post-peak stage(Ⅳ),the acoustic emission characteristics of the stages were different for different water contents.The characteristic value of acoustic emission key point frequency gradually decreases,and the damage degree of the specimen increases,corresponding to low water content—high main frequency—low damage and high water content—low main frequency—high damage.展开更多
Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines...Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.展开更多
Kerogen types exert a decisive effect on the onset and capacity of hydrocarbon generation of source rocks.Lacustrine source rocks in the Liaohe Western Depression are characterized by thick deposition,high total organ...Kerogen types exert a decisive effect on the onset and capacity of hydrocarbon generation of source rocks.Lacustrine source rocks in the Liaohe Western Depression are characterized by thick deposition,high total organic carbon(TOC)content,various kerogen types,and a wide range of thermal maturity.Consequently,their hydrocarbon generation potential and resource estimation can be misinterpreted.In this study,geochemical tests,numerical analysis,hydrocarbon generation kinetics,and basin modeling were integrated to investigate the differential effects of kerogen types on the hydrocarbon generation potential of lacustrine source rocks.Optimized hydrocarbon generation and expulsion(HGE)models of different kerogen types were established quantitatively upon abundant Rock-Eval/TOC/vitrinite reflectance(R_(o))datasets.Three sets of good-excellent source rocks deposited in the fourth(Es4),third(Es3),and first(Es1)members of Paleogene Shahejie Formation,are predominantly types I-II_(1),II_(1)-II_(2),and II-III,respectively.The activation energy of types I-II_(2)kerogen is concentrated(180-230 kcal/mol),whereas that of type III kerogen is widely distributed(150-280 kcal/mol).The original hydrocarbon generation potentials of types I,II_(1),II_(2),and III kerogens are 790,510,270,and 85 mg/g TOC,respectively.The Ro values of the hydrocarbon generation threshold for type I-III source rocks gradually increase from 0.42%to 0.74%,and Ro values of the hydrocarbon expulsion threshold increase from 0.49%to 0.87%.Types I and II_(1)source rocks are characterized by earlier hydrocarbon generation,more rapid hydrocarbon expulsion,and narrower hydrocarbon generation windows than types II_(2)and III source rocks.The kerogen types also affect the HGE history and resource potential.Three types(conventional,tight,and shale oil/gas)and three levels(realistic,expected,and prospective)of hydrocarbon resources of different members in the Liaohe Western Depression are evaluated.Findings suggest that the Es3 member has considerable conventional and unconventional hydrocarbon resources.This study can quantitatively characterize the hydrocarbon generation potential of source rocks with different kerogen types,and facilitate a quick and accurate assessment of hydrocarbon resources,providing strategies for future oil and gas exploration.展开更多
In the context of repositories for nuclear waste,understanding the behavior of gas migration through clayey rocks with inherent anisotropy is crucial for assessing the safety of geological disposal facilities.The prim...In the context of repositories for nuclear waste,understanding the behavior of gas migration through clayey rocks with inherent anisotropy is crucial for assessing the safety of geological disposal facilities.The primary mechanism for gas breakthrough is the opening of micro-fractures due to high gas pressure.This occurs at gas pressures lower than the combined strength of the rock and its minimum principal stress under external loading conditions.To investigate the mechanism of microscale mode-I ruptures,it is essential to incorporate a multiscale approach that includes subcritical microcracks in the modeling framework.In this contribution,we derive the model from microstructures that contain periodically distributed microcracks within a porous material.The damage evolution law is coupled with the macroscopic poroelastic system by employing the asymptotic homogenization method and considering the inherent hydro-mechanical(HM)anisotropy at the microscale.The resulting permeability change induced by fracture opening is implicitly integrated into the gas flow equation.Verification examples are presented to validate the developed model step by step.An analysis of local macroscopic response is undertaken to underscore the influence of factors such as strain rate,initial damage,and applied stress,on the gas migration process.Numerical examples of direct tension tests are used to demonstrate the model’s efficacy in describing localized failure characteristics.Finally,the simulation results for preferential gas flow reveal the robustness of the two-scale model in explicitly depicting gas-induced fracturing in anisotropic clayey rocks.The model successfully captures the common behaviors observed in laboratory experiments,such as a sudden drop in gas injection pressure,rapid build-up of downstream gas pressure,and steady-state gas flow following gas breakthrough.展开更多
基金funded by the National Natural Science Foundation of China (No. 42172308, No.51779018)the Youth Innovation Promotion Association CAS (No. 2022331)the Science and Technology Research and Development Program of China State Railway Group Co., Ltd. (No. J2022G002)。
文摘Red-bed mudstone, prevalent in southwest China, poses a formidable challenge due to its hydrophilic clay minerals, resulting in expansion, deformation, and cracking upon exposure to moisture. This study addresses uplift deformation disasters in high-speed railways by employing a moisture diffusion-deformation-fracture coupling model based on the finite-discrete element method(FDEM). The model integrates the influence of cracks on moisture diffusion. The investigation into various excavation depths reveals a direct correlation between depth and the formation of tensile cracks at the bottom of the railway cutting. These cracks expedite moisture migration, significantly impacting the temporal and spatial evolution of the moisture field. Additionally, crack expansion dominates hygroscopic deformation, with the lateral coordinate of the crack zone determining peak vertical displacement. Furthermore, key factors influencing deformation in railway cuttings, including the swelling factor and initial moisture content at the bottom of the cutting, are explored. The number of tensile and shear cracks increases with greater excavation depth, particularly concerning shear cracks. Higher swelling factors and initial moisture contents result in an increased total number of cracks, predominantly shear cracks. Numerical calculations provide valuable insights, offering a scientific foundation and directional guidance for the precise prevention, control, prediction, and comprehensive treatment of mudstone-related issues in high-speed railways.
基金funded by the National Natural Science Foundation of China(No.42172308)the Youth Innovation Promotion Association CAS(No.2022331).
文摘In this study,the axial swelling strain of red-bed mudstone under different vertical stresses are measured by swell-under-load method,and the microstructure of mudstone after hygroscopic swelling is studied by mercury intrusion porosimetry(MIP).The weakening coefficient and Weibull distribution function are introduced into the coupling model of mudstone moisture diffusion-swelling deformation-fracture based on finite-discrete element method(FDEM).The weakening effect of moisture on mudstone's mechanical parameters,as well as the heterogeneity of swelling deformation and stress distribution,is considered.The microcrack behavior and energy evolution of mudstone during hygroscopic swelling deformation under different vertical stresses are studied.The results show that the axial swelling strain of mudstone decreases with increase of the vertical stress.At low vertical stresses,moisture absorption in mudstone leads to formation of cracks caused by hydration-induced expansion.Under high vertical stresses,a muddy sealing zone forms on the mudstone surface,preventing further water infiltration.The simulation results of mudstone swelling deformation also demonstrate that it involves both swelling of the mudstone matrix and swelling caused by crack expansion.Notably,crack expansion plays a dominant role in mudstone swelling.With increasing vertical stress,the cracks in mudstone change from tensile cracks to shear cracks,resulting in a significant reduction in the total number of cracks.While the evolution of mudstone kinetic energy shows similarities under different vertical stresses,the evolution of strain energy varies significantly due to the presence of different types of cracks in the mudstone.The findings provide a theoretical basis for understanding the hygroscopic swelling deformation mechanism of red-bed mudstone at various depths.
基金supported by the National Key Research and Development Program of China(No.2023YFC2907600)the National Natural Science Foundation of China(Nos.42077267,42277174 and 52074164)+2 种基金the Natural Science Foundation of Shandong Province,China(No.ZR2020JQ23)the Opening Project of State Key Laboratory of Explosion Science and Technology,Beijing Institute of Technology(No.KFJJ21-02Z)the Fundamental Research Funds for the Central Universities,China(No.2022JCCXSB03).
文摘The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.
基金National Natural Science Foundation of China (No. 52204101)Natural Science Foundation of Shandong Province (No. ZR2022QE137)Open Project of State Key Laboratory for Geomechanics and Deep Underground Engineering in CUMTB (No. SKLGDUEK2023).
文摘Uniaxial compression tests and cyclic loading acoustic emission tests were conducted on 20%,40%,60%,80%,dry and saturated muddy sandstone by using a creep impact loading system to investigate the mechanical properties and acoustic emission characteristics of soft rocks with different water contents under dynamic disturbance.The mechanical properties and acoustic emission characteristics of muddy sandstones at different water contents were analysed.Results of experimental studies show that water is a key factor in the mechanical properties of rocks,softening them,increasing their porosity,reducing their brittleness and increasing their plasticity.Under uniaxial compression,the macroscopic damage characteristics of the muddy sandstone change from mono-bevel shear damage and‘X’type conjugate bevel shear damage to a roadway bottom-drum type damage as the water content increases.Dynamic perturbation has a strengthening effect on the mechanical properties of samples with 60%and less water content,and a weakening effect on samples with 80%and more water content,but the weakening effect is not obvious.Macroscopic damage characteristics of dry samples remain unchanged,water samples from shear damage and tensile–shear composite damage gradually transformed into cleavage damage,until saturation transformation monoclinic shear damage.The evolution of acoustic emission energy and event number is mainly divided into four stages:loading stage(Ⅰ),dynamic loading stage(Ⅱ),yield failure stage(Ⅲ),and post-peak stage(Ⅳ),the acoustic emission characteristics of the stages were different for different water contents.The characteristic value of acoustic emission key point frequency gradually decreases,and the damage degree of the specimen increases,corresponding to low water content—high main frequency—low damage and high water content—low main frequency—high damage.
基金supported by Science and Technology Project of Yunnan Provincial Transportation Department(Grant No.25 of 2018)the National Natural Science Foundation of China(Grant No.52279107)The authors are grateful for the support by the China Scholarship Council(CSC No.202206260203 and No.201906690049).
文摘Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.
基金This research is supported by the Joint Fund of the National Natural Science Foundation of China(grant number U19B6003-02)the Cooperation Program of PetroChina Liaohe Oilfield Company(grant Number HX20180604)the AAPG Foundation Grants-in-Aid Program(grant number 22269437).This study has benefited considerably from PetroChina Liaohe Oilfield Company for data support.We also thank the editor and the anonymous reviewers for their professional suggestions and comments.
文摘Kerogen types exert a decisive effect on the onset and capacity of hydrocarbon generation of source rocks.Lacustrine source rocks in the Liaohe Western Depression are characterized by thick deposition,high total organic carbon(TOC)content,various kerogen types,and a wide range of thermal maturity.Consequently,their hydrocarbon generation potential and resource estimation can be misinterpreted.In this study,geochemical tests,numerical analysis,hydrocarbon generation kinetics,and basin modeling were integrated to investigate the differential effects of kerogen types on the hydrocarbon generation potential of lacustrine source rocks.Optimized hydrocarbon generation and expulsion(HGE)models of different kerogen types were established quantitatively upon abundant Rock-Eval/TOC/vitrinite reflectance(R_(o))datasets.Three sets of good-excellent source rocks deposited in the fourth(Es4),third(Es3),and first(Es1)members of Paleogene Shahejie Formation,are predominantly types I-II_(1),II_(1)-II_(2),and II-III,respectively.The activation energy of types I-II_(2)kerogen is concentrated(180-230 kcal/mol),whereas that of type III kerogen is widely distributed(150-280 kcal/mol).The original hydrocarbon generation potentials of types I,II_(1),II_(2),and III kerogens are 790,510,270,and 85 mg/g TOC,respectively.The Ro values of the hydrocarbon generation threshold for type I-III source rocks gradually increase from 0.42%to 0.74%,and Ro values of the hydrocarbon expulsion threshold increase from 0.49%to 0.87%.Types I and II_(1)source rocks are characterized by earlier hydrocarbon generation,more rapid hydrocarbon expulsion,and narrower hydrocarbon generation windows than types II_(2)and III source rocks.The kerogen types also affect the HGE history and resource potential.Three types(conventional,tight,and shale oil/gas)and three levels(realistic,expected,and prospective)of hydrocarbon resources of different members in the Liaohe Western Depression are evaluated.Findings suggest that the Es3 member has considerable conventional and unconventional hydrocarbon resources.This study can quantitatively characterize the hydrocarbon generation potential of source rocks with different kerogen types,and facilitate a quick and accurate assessment of hydrocarbon resources,providing strategies for future oil and gas exploration.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.12302503 and U20A20266)Scientific and Technological Research Projects in Sichuan Province,China(Grant No.2023ZYD0154).
文摘In the context of repositories for nuclear waste,understanding the behavior of gas migration through clayey rocks with inherent anisotropy is crucial for assessing the safety of geological disposal facilities.The primary mechanism for gas breakthrough is the opening of micro-fractures due to high gas pressure.This occurs at gas pressures lower than the combined strength of the rock and its minimum principal stress under external loading conditions.To investigate the mechanism of microscale mode-I ruptures,it is essential to incorporate a multiscale approach that includes subcritical microcracks in the modeling framework.In this contribution,we derive the model from microstructures that contain periodically distributed microcracks within a porous material.The damage evolution law is coupled with the macroscopic poroelastic system by employing the asymptotic homogenization method and considering the inherent hydro-mechanical(HM)anisotropy at the microscale.The resulting permeability change induced by fracture opening is implicitly integrated into the gas flow equation.Verification examples are presented to validate the developed model step by step.An analysis of local macroscopic response is undertaken to underscore the influence of factors such as strain rate,initial damage,and applied stress,on the gas migration process.Numerical examples of direct tension tests are used to demonstrate the model’s efficacy in describing localized failure characteristics.Finally,the simulation results for preferential gas flow reveal the robustness of the two-scale model in explicitly depicting gas-induced fracturing in anisotropic clayey rocks.The model successfully captures the common behaviors observed in laboratory experiments,such as a sudden drop in gas injection pressure,rapid build-up of downstream gas pressure,and steady-state gas flow following gas breakthrough.
