Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass...Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass is very important to elucidate the permeability evolution mechanisms, which is difficult to realize in laboratory, but easy to be achieved in numerical simulations. Therefore, the particle flow code （PFC）, a discrete element method, is used to simulate permeability behaviors of rock materials in this study. Owe to the limitation of the existed solid-fluid coupling algorithm in PFC, an improved flow-coupling algorithm is presented to better reflect the preferential flow in rock fractures. The comparative analysis is conducted between original and improved algorithm when simulating rock permeability evolution during triaxial compression, showing that the improved algorithm can better describe the experimental phenomenon. Furthermore, the evolution of pore pressure and flow rate distribution during the flow process are analyzed by using the improved algorithm. It is concluded that during the steady flow process in the fractured specimen, the pore pressure and flow rate both prefer transmitting through the fractures rather than rock matrix. Based on the results, fractures are divided into the following three types： I） fractures link to both the inlet and outlet, II） fractures only link to the inlet, and III） fractures only link to the outlet. The type I fracture is always the preferential propagating path for both the pore pressure and flow rate. For type II fractures, the pore pressure increases and then becomes steady. However, the flow rate increases first and begins to decrease after the flow reaches the stop end of the fracture and finally vanishes. There is no obvious pore pressure or flow rate concentration within type III fractures.展开更多
It is important to calibrate micro-parameters for applying partied flow code(PFC)to study mechanical characteristics and failure mechanism of rock materials.Uniform design method is firstly adopted to determine the mi...It is important to calibrate micro-parameters for applying partied flow code(PFC)to study mechanical characteristics and failure mechanism of rock materials.Uniform design method is firstly adopted to determine the microscopic parameters of parallel-bonded particle model for three-dimensional discrete element particle flow code(PFC3D).Variation ranges of microscopic of the microscopic parameters are created by analyzing the effects of microscopic parameters on macroscopic parameters(elastic modulus E,Poisson ratio v,uniaxial compressive strengthσc,and ratio of crack initial stress to uniaxial compressive strengthσci/σc)in order to obtain the actual uniform design talbe.The calculation equations of the microscopic and macroscopic parameters of rock materials can be established by the actual uniform design table and the regression analysis and thus the PFC3D microscopic parameters can be quantitatively determined.The PFC3D simulated results of the intact and pre-cracked rock specimens under uniaxial and triaxial compressions(including the macroscopic mechanical parameters,stress−strain curves and failure process)are in good agreement with experimental results,which can prove the validity of the calculation equations of microscopic and macroscopic parameters.展开更多
Hydraulic fracturing,as a key technology of deep energy exploitation,accelerates the rapid development of the modern petroleum industry.To study the mechanisms of hydraulic fracture propagation and rock failure mode o...Hydraulic fracturing,as a key technology of deep energy exploitation,accelerates the rapid development of the modern petroleum industry.To study the mechanisms of hydraulic fracture propagation and rock failure mode of the vertical well hydraulic fracturing,the true triaxial hydraulic fracturing test and numerical simulation are carried out,and the influence of the principal stress difference,water injection displacement,perforation angle and natural fracture on fracture propagation is analyzed.The results show that the fracture propagation mode of limestone is mainly divided into two types:the single vertical fracture and the transverse-longitudinal crossed complex fracture.Under high displacement,the fracturing pressure is larger,and the secondary fracture is more likely to occur,while variable displacement loading is more likely to induce fracture network.Meanwhile,the amplitude of acoustic emission(AE)waveform of limestone during fracturing is between 0.01 and 0.02 mV,and the main frequency is maintained in the range of 230−300 kHz.When perforation angleθ=45°,it is easy to produce the T-type fracture that connects with the natural fracture,while X-type cracks are generated whenθ=30°.The results can be used as a reference for further study on the mechanism of limestone hydraulic fracturing.展开更多
The complexity of a rock masses structure can lead to high uncertainties and risk during underground engineering construction.Laboratory tests on fractured rock-like materials containing a tunnel were conducted,and tw...The complexity of a rock masses structure can lead to high uncertainties and risk during underground engineering construction.Laboratory tests on fractured rock-like materials containing a tunnel were conducted,and twodimensional particle flow models were established.The principal stress and principal strain distributions surrounding the four-arc-shaped and inverted U-shaped tunnels were investigated,respectively.Numerical results indicated that the dip angle combination of preexisting fractures directly affects the principal stress,principal strain distribution and the failure characteristics around the tunnel.The larger the absolute value of the preexisting fracture inclination angle,the higher the crushing degree of compression splitting near the hance and the larger the V-shaped failure zone.With a decrease in the absolute value of the preexisting fracture inclination angle,the compressive stress concentration of the sidewall with preexisting fractures gradually increases.The types of cracks initiated around the four-arc-shaped tunnel and the inverted U-shape tunnel are different.When the fractures are almost vertical,they have a significant influence on the stress of the sidewall force of the four-arc-shaped tunnel.When the fractures are almost horizontal,they have a significant influence on the stress of the sidewall of the inverted U-shaped tunnel.The findings provide a theoretical support for the local strengthening design of the tunnel supporting structure.展开更多
基于颗粒离散元计算架构建立单矿物石英岩的颗粒簇模型(grain based model,GBM),以石英岩微裂纹发展规律为约束优化矿物模型的弹性模量及其边界强度参数,开展热-固耦合下石英岩的高温损伤及强度劣化细观机理研究,获得以下结论:(1)优化...基于颗粒离散元计算架构建立单矿物石英岩的颗粒簇模型(grain based model,GBM),以石英岩微裂纹发展规律为约束优化矿物模型的弹性模量及其边界强度参数,开展热-固耦合下石英岩的高温损伤及强度劣化细观机理研究,获得以下结论:(1)优化矿物弹性模量及其边界强度参数后可很好地再现石英岩高温作用下微裂纹演化规律:首先从石英晶粒边界萌生,逐渐扩展至晶粒内部,最终形成裂纹网络。(2)单轴压缩模拟中石英岩弹性模量随温度的升高总体上呈现出先增大后减小的趋势,泊松比表现为先减小后增大,且泊松比较弹性模量对热损伤表现出更高的敏感性。(3)随着温度的升高,石英岩单轴压缩过程中微裂纹的离散性增加,破坏模式从多斜面剪切破坏逐渐转变为无明显宏观裂纹的分散式拉伸(劈裂)破坏。(4)岩石的力学强度劣化规律与其本身的致密性及矿物成分、含量密切相关,与花岗岩相比石英岩热损伤温度阈值高,约为200℃。(5)3种不同石英含量的岩石的热损伤结果表明:矿物边界和力学强度较低的矿物内部首先发生热破裂,矿物之间的镶嵌结构对石英矿物的热损伤有抑制作用。展开更多
基金Project（BK20150005） supported by the Natural Science Foundation of Jiangsu Province for Distinguished Young Scholars, China Project（2015XKZD05） supported by the Fundamental Research Funds for the Central Universities, China
文摘Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass is very important to elucidate the permeability evolution mechanisms, which is difficult to realize in laboratory, but easy to be achieved in numerical simulations. Therefore, the particle flow code （PFC）, a discrete element method, is used to simulate permeability behaviors of rock materials in this study. Owe to the limitation of the existed solid-fluid coupling algorithm in PFC, an improved flow-coupling algorithm is presented to better reflect the preferential flow in rock fractures. The comparative analysis is conducted between original and improved algorithm when simulating rock permeability evolution during triaxial compression, showing that the improved algorithm can better describe the experimental phenomenon. Furthermore, the evolution of pore pressure and flow rate distribution during the flow process are analyzed by using the improved algorithm. It is concluded that during the steady flow process in the fractured specimen, the pore pressure and flow rate both prefer transmitting through the fractures rather than rock matrix. Based on the results, fractures are divided into the following three types： I） fractures link to both the inlet and outlet, II） fractures only link to the inlet, and III） fractures only link to the outlet. The type I fracture is always the preferential propagating path for both the pore pressure and flow rate. For type II fractures, the pore pressure increases and then becomes steady. However, the flow rate increases first and begins to decrease after the flow reaches the stop end of the fracture and finally vanishes. There is no obvious pore pressure or flow rate concentration within type III fractures.
基金Projects(51474251,51874351)supported by the National Natural Science Foundation,China。
文摘It is important to calibrate micro-parameters for applying partied flow code(PFC)to study mechanical characteristics and failure mechanism of rock materials.Uniform design method is firstly adopted to determine the microscopic parameters of parallel-bonded particle model for three-dimensional discrete element particle flow code(PFC3D).Variation ranges of microscopic of the microscopic parameters are created by analyzing the effects of microscopic parameters on macroscopic parameters(elastic modulus E,Poisson ratio v,uniaxial compressive strengthσc,and ratio of crack initial stress to uniaxial compressive strengthσci/σc)in order to obtain the actual uniform design talbe.The calculation equations of the microscopic and macroscopic parameters of rock materials can be established by the actual uniform design table and the regression analysis and thus the PFC3D microscopic parameters can be quantitatively determined.The PFC3D simulated results of the intact and pre-cracked rock specimens under uniaxial and triaxial compressions(including the macroscopic mechanical parameters,stress−strain curves and failure process)are in good agreement with experimental results,which can prove the validity of the calculation equations of microscopic and macroscopic parameters.
基金Projects(51879148,51709159,51911530214)supported by the National Natural Science Foundation of ChinaProject(2019GSF111030)supported by Shandong Provincial Key R&D Program of ChinaProject(KT201804)supported by the Project of Special Fund for Science and Technology of Water Resources Department of Guizhou Province,China。
文摘Hydraulic fracturing,as a key technology of deep energy exploitation,accelerates the rapid development of the modern petroleum industry.To study the mechanisms of hydraulic fracture propagation and rock failure mode of the vertical well hydraulic fracturing,the true triaxial hydraulic fracturing test and numerical simulation are carried out,and the influence of the principal stress difference,water injection displacement,perforation angle and natural fracture on fracture propagation is analyzed.The results show that the fracture propagation mode of limestone is mainly divided into two types:the single vertical fracture and the transverse-longitudinal crossed complex fracture.Under high displacement,the fracturing pressure is larger,and the secondary fracture is more likely to occur,while variable displacement loading is more likely to induce fracture network.Meanwhile,the amplitude of acoustic emission(AE)waveform of limestone during fracturing is between 0.01 and 0.02 mV,and the main frequency is maintained in the range of 230−300 kHz.When perforation angleθ=45°,it is easy to produce the T-type fracture that connects with the natural fracture,while X-type cracks are generated whenθ=30°.The results can be used as a reference for further study on the mechanism of limestone hydraulic fracturing.
基金Project(41807241) supported by the National Natural Science Foundation of ChinaProject(2021M693544) supported by China Postdoctoral Science FoundationProject(2022JM-160) supported by the Natural Science Basic Research Program of Shaanxi,China。
文摘The complexity of a rock masses structure can lead to high uncertainties and risk during underground engineering construction.Laboratory tests on fractured rock-like materials containing a tunnel were conducted,and twodimensional particle flow models were established.The principal stress and principal strain distributions surrounding the four-arc-shaped and inverted U-shaped tunnels were investigated,respectively.Numerical results indicated that the dip angle combination of preexisting fractures directly affects the principal stress,principal strain distribution and the failure characteristics around the tunnel.The larger the absolute value of the preexisting fracture inclination angle,the higher the crushing degree of compression splitting near the hance and the larger the V-shaped failure zone.With a decrease in the absolute value of the preexisting fracture inclination angle,the compressive stress concentration of the sidewall with preexisting fractures gradually increases.The types of cracks initiated around the four-arc-shaped tunnel and the inverted U-shape tunnel are different.When the fractures are almost vertical,they have a significant influence on the stress of the sidewall force of the four-arc-shaped tunnel.When the fractures are almost horizontal,they have a significant influence on the stress of the sidewall of the inverted U-shaped tunnel.The findings provide a theoretical support for the local strengthening design of the tunnel supporting structure.
文摘基于颗粒离散元计算架构建立单矿物石英岩的颗粒簇模型(grain based model,GBM),以石英岩微裂纹发展规律为约束优化矿物模型的弹性模量及其边界强度参数,开展热-固耦合下石英岩的高温损伤及强度劣化细观机理研究,获得以下结论:(1)优化矿物弹性模量及其边界强度参数后可很好地再现石英岩高温作用下微裂纹演化规律:首先从石英晶粒边界萌生,逐渐扩展至晶粒内部,最终形成裂纹网络。(2)单轴压缩模拟中石英岩弹性模量随温度的升高总体上呈现出先增大后减小的趋势,泊松比表现为先减小后增大,且泊松比较弹性模量对热损伤表现出更高的敏感性。(3)随着温度的升高,石英岩单轴压缩过程中微裂纹的离散性增加,破坏模式从多斜面剪切破坏逐渐转变为无明显宏观裂纹的分散式拉伸(劈裂)破坏。(4)岩石的力学强度劣化规律与其本身的致密性及矿物成分、含量密切相关,与花岗岩相比石英岩热损伤温度阈值高,约为200℃。(5)3种不同石英含量的岩石的热损伤结果表明:矿物边界和力学强度较低的矿物内部首先发生热破裂,矿物之间的镶嵌结构对石英矿物的热损伤有抑制作用。
