Numerical modeling of blast-induced rock fragmentation in deep mining with 3D and 2D FEM method approaches

To optimize the excavation of rock using underground blasting techniques,a reliable and simplified approach for modeling rock fragmentation is desired.This paper presents a multistep experimentalnumerical methodology for simplifying the three-dimensional(3D)to two-dimensional(2D)quasiplane-strain problem and reducing computational costs by more than 100-fold.First,in situ tests were conducted involving single-hole and free-face blasting of a dolomite rock mass in a 1050-m-deep mine.The results were validated by laser scanning.The craters were then compared with four analytical models to calculate the radius of the crushing zone.Next,a full 3D model for single-hole blasting was prepared and validated by simulating the crack length and the radius of the crushing zone.Based on the stable crack propagation zones observed in the 3D model and experiments,a 2D model was prepared.The properties of the high explosive(HE)were slightly reduced to match the shape and number of radial cracks and crushing zone radius between the 3D and 2D models.The final methodology was used to reproduce various cut-hole blasting scenarios and observe the effects of residual cracks in the rock mass on further fragmentation.The presence of preexisting cracks was found to be crucial for fragmentation,particularly when the borehole was situated near a free rock face.Finally,an optimization study was performed to determine the possibility of losing rock continuity at different positions within the well in relation to the free rock face.

Parameter calibration of the tensile-shear interactive damage constitutive model for sandstone failure

The tensile-shear interactive damage(TSID)model is a novel and powerful constitutive model for rock-like materials.This study proposes a methodology to calibrate the TSID model parameters to simulate sandstone.The basic parameters of sandstone are determined through a series of static and dynamic tests,including uniaxial compression,Brazilian disc,triaxial compression under varying confining pressures,hydrostatic compression,and dynamic compression and tensile tests with a split Hopkinson pressure bar.Based on the sandstone test results from this study and previous research,a step-by-step procedure for parameter calibration is outlined,which accounts for the categories of the strength surface,equation of state(EOS),strain rate effect,and damage.The calibrated parameters are verified through numerical tests that correspond to the experimental loading conditions.Consistency between numerical results and experimental data indicates the precision and reliability of the calibrated parameters.The methodology presented in this study is scientifically sound,straightforward,and essential for improving the TSID model.Furthermore,it has the potential to contribute to other rock constitutive models,particularly new user-defined models.

黄芩苷通过ROS依赖性调节RhoA/ROCK通路保护氯化钴诱导心肌细胞损伤的实验研究

目的:本研究旨在探讨黄芩苷(Baicalin)对氯化钴(CoCl_(2))诱导H9c2心肌细胞缺氧损伤的作用机制。方法:采用CoCl_(2)建立心肌细胞缺氧损伤模型,并分别加入不同浓度的黄芩苷培养。将正常氧培养设置为对照组,CoCl_(2)培养H9c2心肌细胞设置为CoCl_(2)组;Baicalin、Y27632(Rho激酶抑制剂)预处理的H9c2缺氧心肌细胞设置为CoCl_(2)+Baicalin组、CoCl_(2)+Y27632组、CoCl_(2)+Baicalin+Y27632组。通过细胞毒性检测试剂盒(CCK8)检测细胞活性;荧光探针测定细胞中活性氧(ROS)的表达;WST-8检测心肌细胞超氧化歧化酶(SOD)活力;TBA法测定丙二醇(MDA)浓度;同时采用WesternBlot方法分析RhoA、ROCK1、ROCK2、TNF-α、IL-1β蛋白表达情况。结果:处理H9c2细胞24 h后,1 000μmol/L的CoCl_(2)和75μmol/L黄芩苷对治疗心肌细胞缺氧损伤具有较好的细胞活力。CoCl_(2)组细胞活性明显低于对照组,加入黄芩苷后可显著提高细胞活性(P<0.05);与对照组相比,CoCl_(2)组心肌细胞可上调ROS、MDA表达,下调SOD活力,升高RhoA、ROCK1、ROCK2、TNF-α、IL-1β蛋白表达水平,加入Baicalin后可逆转上诉蛋白的表达水平;与CoCl_(2)组相比,CoCl_(2)+Baicalin可抑制ROS、MDA的表达,上调SOD含量,下调RhoA、ROCK1、ROCK2、TNF-α、IL-1β蛋白表达水平,而加入Y27632(Rho激酶抑制剂)后则可显著增强Baicalin带来的保护作用。结论:Baicalin可减轻CoCl_(2)诱导H9c2心肌细胞缺氧损伤的炎症、氧化应激反应,其机制与抑制ROS依赖性RhoA/ROCK通路有关。

利多卡因调节RhoA/ROCK信号通路对缺血性脑卒中大鼠的神经保护作用

目的探讨利多卡因调节RhoA/ROCK信号通路对缺血性脑卒中(IS)大鼠神经保护作用。方法通过线栓法建立模型IS大鼠,随机分为模型组、利多卡因低(利多卡因-L,5 mg/kg)、中(利多卡因-M,10 mg/kg)、高剂量(利多卡因-H,20 mg/kg)组以及利多卡因-H+溶血磷脂酸(LPA)(20 mg/kg利多卡因+50μmoL/L LPA)组,以仅分离但不插尼龙线的大鼠为假手术组;干预结束后,对大鼠进行神经功能评分、脑梗死体积、脑含水量测定;分离脑组织,观察神经元变化、神经细胞凋亡以及RhoA、ROCK蛋白表达。结果与假手术组相比,模型组大鼠Zea-Longa评分(3.28±0.33、0.14±0.02)、脑梗死体积(38.09±3.81、0.00±0.00)、含水量(84.11±2.39、62.95±0.65)、神经细胞凋亡率(39.51±3.96、5.14±0.52)、RhoA(1.16±0.12、0.24±0.03)、ROCK蛋白(2.15±0.22、0.77±0.08)表达增加(P<0.05);与模型组相比,利多卡因组大鼠Zea-Longa评分(2.37±0.24、1.52±0.16、0.87±0.09)、脑梗死体积(27.42±2.76、16.64±1.68、8.22±0.83)、含水量(75.24±1.46、68.34±1.02、63.15±0.86)、神经细胞凋亡率(27.08±2.71、18.84±1.89、9.47±0.95)、RhoA(0.78±0.08、0.52±0.06、0.26±0.03)、ROCK蛋白(1.66±0.17、1.24±0.13、0.86±0.09)表达显著降低,组间差异有统计学意义(P<0.05);LPA逆转了利多卡因-H对IS大鼠的保护作用。结论利多卡因通过抑制RhoA/ROCK信号通路对IS大鼠发挥神经保护作用。

Development Model and Quantitative Prediction of Igneous Rock in S Oilfield, the Bohai Sea Area

The structure of S Oilfield of Bohai sea area in China is located in the low uplift slope zone of the Dongwa Laibei area in the Huanghekou depression. The target interval of the oilfield belongs to the Paleogene layered structural oil reservoir under the influence of igneous rocks. Due to the widespread development of igneous rocks, and the diversity of igneous rock types, the structure is complex, and the lateral changes are intense. The differences in volcanic eruption intensity and volcanic materials made the igneous rock masses overlap each other, with rapid changes in lithofacies, and complex seismic reflection structures. Therefore, it was difficult to carry out fine characterization of igneous rocks. Based on the lithofacies model, lithofacies combination types, and seismic response characteristics of volcanic mechanisms, this paper summarized three types of development models of volcanic mechanisms in the study area. At the same time, through technical means such as stratigraphic slicing, seismic attribute extraction, and attribute optimization analysis, the spatial distribution characteristics of various facies zones of igneous rocks at different stages were described in detail, achieving precise prediction of igneous rocks in Bozhong S Oilfield. During the development process of the oilfield, the technical research results guided the optimization of well locations and provided technical support for the efficient development of the oilfield.

Rock skeleton models and seismic porosity inversion

By substituting rock skeleton modulus expressions into Gassmann approximate fluid equation, we obtain a seismic porosity inversion equation. However, conventional rock skeleton models and their expressions are quite different from each other, resuling in different seismic porosity inversion equations, potentially leading to difficulties in correctly applying them and evaluating their results. In response to this, a uniform relation with two adjusting parameters suitable for all rock skeleton models is established from an analysis and comparison of various conventional rock skeleton models and their expressions including the Eshelby-Walsh, Pride, Geertsma, Nur, Keys-Xu, and Krief models. By giving the two adjusting parameters specific values, different rock skeleton models with specific physical characteristics can be generated. This allows us to select the most appropriate rock skeleton model based on geological and geophysical conditions, and to develop more wise seismic porosity inversion. As an example of using this method for hydrocarbon prediction and fluid identification, we apply this improved porosity inversion, associated with rock physical data and well log data, to the ZJ basin. Research shows that the existence of an abundant hydrocarbon reservoir is dependent on a moderate porosity range, which means we can use the results of seismic porosity inversion to identify oil reservoirs and dry or water-saturated reservoirs. The seismic inversion results are closely correspond to well log porosity curves in the ZJ area, indicating that the uniform relations and inversion methods proposed in this paper are reliable and effective.

Rock physics modeling of heterogeneous carbonatereservoirs： porosity estimation and hydrocarbon detection

In heterogeneous natural gas reservoirs, gas is generally present as small patchlike pockets embedded in the water-saturated host matrix. This type of heterogeneity, also called "patchy saturation", causes significant seismic velocity dispersion and attenuation. To establish the relation between seismic response and type of fluids, we designed a rock physics model for carbonates. First, we performed CT scanning and analysis of the fluid distribution in the partially saturated rocks. Then, we predicted the quantitative relation between the wave response at different frequency ranges and the basic lithological properties and pore fluids. A rock physics template was constructed based on thin section analysis of pore structures and seismic inversion. This approach was applied to the limestone gas reservoirs of the right bank block of the Amu Darya River. Based on poststack wave impedance and prestack elastic parameter inversions, the seismic data were used to estimate rock porosity and gas saturation. The model results were in good agreement with the production regime of the wells.

Failure characteristics of three-body model composed of rock and coal with different strength and stiffness

Four different types of three-body model composed of rock and coal with different strength and stiffness were established in order to study the failure characteristics of compound model such as roof-coal-floor. Through stress analysis of the element with variable strength and stiffness extracted from the strong-weak interface, the tri-axial compressive strength of the weak body and strong body near the interface as well as the areas away from the contact surface was found. Then, on the basis of three-dimensional fast Lagrangian method of continua and strain softening constitutive model composed of Coulomb-Mohr shear failure with tensile cut-off, stress and strain relationship of the four three-body combined models were analyzed under different confining pressures by numerical simulation. Finally, the different features of local shear zones and plastic failure areas of the four different models and their development trend with increasing confining pressure were discussed. The results show that additional stresses are derived due to the lateral deformation constraints near the strong-weak interface area, which results in the strength increasing in weak body and strength decreasing in strong body. The weakly consolidated soft rock and coal cementation exhibit significant strain softening behavior and bear compound tension-shear failure under uni-axial compression. With the increase of confining pressure, the tensile failure disappears from the model, and the failure type of composed model changes to local shear failure with different number of shearing bands and plastic failure zones. This work shows important guiding significance for the mechanism study of seismic, rock burst, and coal bump.

Brittleness index and seismic rock physics model for anisotropic tight-oil sandstone reservoirs

Brittleness analysis becomes important when looking for sweet spots in tightoil sandstone reservoirs. Hence, appropriate indices are required as accurate brittleness evaluation criteria. We construct a seismic rock physics model for tight-oil sandstone reservoirs with vertical fractures. Because of the complexities in lithology and pore structure and the anisotropic characteristics of tight-oil sandstone reservoirs, the proposed model is based on the solid components, pore connectivity, pore type, and fractures to better describe the sandstone reservoir microstructure. Using the model, we analyze the brittleness sensitivity of the elastic parameters in an anisotropic medium and establish a new brittleness index. We show the applicability of the proposed brittleness index for tight-oil sandstone reservoirs by considering the brittleness sensitivity, the rock physics response characteristics, and cross-plots. Compared with conventional brittleness indexes, the new brittleness index has high brittleness sensitivity and it is the highest in oil-bearing brittle zones with relatively high porosity. The results also suggest that the new brittleness index is much more sensitive to elastic properties variations, and thus can presumably better predict the brittleness characteristics of sweet spots in tight-oil sandstone reservoirs.

Thermal damage constitutive model for rock considering damage threshold and residual strength

With the gradual depletion of mineral resources in the shallow part of the earth,resource exploitation continues to move deeper into the earth,it becomes a hot topic to simulate the whole process of rock strain softening,deformation and failure in deep environment,especially under high temperature and high pressure.On the basis of Lemaitre’s strain-equivalent principle,combined with statistics and damage theory,a statistical constitutive model of rock thermal damage under triaxial compression condition is established.At the same time,taking into account the existing damage model is difficult to reflect residual strength after rock failure,the residual strength is considered in this paper by introducing correction factor of damage variable,the model rationality is also verified by experiments.Analysis of results indicates that the damage evolution curve reflects the whole process of rock micro-cracks enclosure,initiation,expansion,penetration,and the formation of macro-cracks under coupled effect of temperature and confining pressure.Rock thermal damage shows logistic growth function with the increase of temperature.Under the same strain condition,rock total damage decreases with the rise of confining pressure.By studying the electron microscope images(SEM)of rock fracture,it is inferred that 35.40 MPa is the critical confining pressure of brittle to plastic transition for this granite.The model parameter F reflects the average strength of rock,and M reflects the morphological characteristics of rock stress–strain curves.The physical meanings of model parameters are clear and the model is suitable for complex stress states,which provides valuable references for the study of rock deformation and stability in deep engineering.

Effects of confinement on rock mass modulus:A synthetic rock mass modelling(SRM) study

The main objective of this paper is to examine the influence of the applied confining stress on the rock mass modulus of moderately jointed rocks（well interlocked undisturbed rock mass with blocks formed by three or less intersecting joints）. A synthetic rock mass modelling（SRM） approach is employed to determine the mechanical properties of the rock mass. In this approach, the intact body of rock is represented by the discrete element method（DEM）-Voronoi grains with the ability of simulating the initiation and propagation of microcracks within the intact part of the model. The geometry of the preexisting joints is generated by employing discrete fracture network（DFN） modelling based on field joint data collected from the Brockville Tunnel using LiDAR scanning. The geometrical characteristics of the simulated joints at a representative sample size are first validated against the field data, and then used to measure the rock quality designation（RQD）, joint spacing, areal fracture intensity（P21）, and block volumes. These geometrical quantities are used to quantitatively determine a representative range of the geological strength index（GSI）. The results show that estimating the GSI using the RQD tends to make a closer estimate of the degree of blockiness that leads to GSI values corresponding to those obtained from direct visual observations of the rock mass conditions in the field. The use of joint spacing and block volume in order to quantify the GSI value range for the studied rock mass suggests a lower range compared to that evaluated in situ. Based on numerical modelling results and laboratory data of rock testing reported in the literature, a semi-empirical equation is proposed that relates the rock mass modulus to confinement as a function of the areal fracture intensity and joint stiffness.

Prediction of brittleness based on anisotropic rock physics model for kerogen-rich shale

The construction of a shale rock physics model and the selection of an appropriate brittleness index （B/） are two significant steps that can influence the accuracy of brittleness prediction. On one hand, the existing models of kerogen-rich shale are controversial, so a reasonable rock physics model needs to be built. On the other hand, several types of equations already exist for predicting the BI whose feasibility needs to be carefully considered. This study constructed a kerogen-rich rock physics model by performing the self- consistent approximation and the differential effective medium theory to model intercoupled clay and kerogen mixtures. The feasibility of our model was confirmed by comparison with classical models, showing better accuracy. Templates were constructed based on our model to link physical properties and the BL Different equations for the BI had different sensitivities, making them suitable for different types of formations. Equations based on Young＇s Modulus were sensitive to variations in lithology, while those using Lame＇s Coefficients were sensitive to porosity and pore fluids. Physical information must be considered to improve brittleness prediction.

The construction of shale rock physics model and brittleness prediction for high-porosity shale gas-bearing reservoir

Due to the huge differences between the unconventional shale and conventional sand reservoirs in many aspects such as the types and the characteristics of minerals,matrix pores and fluids,the construction of shale rock physics model is significant for the exploration and development of shale reservoirs.To make a better characterization of shale gas-bearing reservoirs,we first propose a new but more suitable rock physics model to characterize the reservoirs.We then use a well A to demonstrate the feasibility and reliability of the proposed rock physics model of shale gas-bearing reservoirs.Moreover,we propose a new brittleness indicator for the high-porosity and organic-rich shale gas-bearing reservoirs.Based on the parameter analysis using the constructed rock physics model,we finally compare the new brittleness indicator with the commonly used Young’s modulus in the content of quartz and organic matter,the matrix porosity,and the types of filled fluids.We also propose a new shale brittleness index by integrating the proposed new brittleness indicator and the Poisson’s ratio.Tests on real data sets demonstrate that the new brittleness indicator and index are more sensitive than the commonly used Young’s modulus and brittleness index for the high-porosity and high-brittleness shale gas-bearing reservoirs.

A review of discrete modeling techniques for fracturing processes in discontinuous rock masses

The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.

A New Example of Retrograde Solubility Model for Carbonate Rocks

Objective The dissolution and precipitation of carbonate during burial diagenetic process controls the reservoir property in deep buried strata. The geological process related with it has become a research focus during recent years. The most important dissolution fluids to carbonates are probably H2S and CO2 as byproducts of sulfate reduction in deep-buried setting with sulfate minerals, but carbonates are more soluble in relatively low temperature, which is the so-called retrograde solubility. Several geological processes can result in the decrease of temperature, including the upward migration of thermal fluids and tectonic uplift.

An analytical model for estimating rock strength parameters from small-scale drilling data

The small-scale drilling technique can be a fast and reliable method to estimate rock strength parameters. It needs to link the operational drilling parameters and strength properties of rock. The parameters such as bit geometry, bit movement, contact frictions and crushed zone affect the estimated parameters.An analytical model considering operational drilling data and effective parameters can be used for these purposes. In this research, an analytical model was developed based on limit equilibrium of forces in a Tshaped drag bit considering the effective parameters such as bit geometry, crushed zone and contact frictions in drilling process. Based on the model, a method was used to estimate rock strength parameters such as cohesion, internal friction angle and uniaxial compressive strength of different rock types from operational drilling data. Some drilling tests were conducted by a portable and powerful drilling machine which was developed for this work. The obtained results for strength properties of different rock types from the drilling experiments based on the proposed model are in good agreement with the results of standard tests. Experimental results show that the contact friction between the cutting face and rock is close to that between bit end wearing face and rock due to the same bit material. In this case,the strength parameters, especially internal friction angle and cohesion, are estimated only by using a blunt bit drilling data and the bit bluntness does not affect the estimated results.

Distinct element modelling of fracture plan control in continuum and jointed rock mass in presplitting method of surface mining

Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surface blast design. The purpose of presplitting is to form a fracture plane across which the radial cracks from the production blast cannot travel. The purpose of this study is to investigate of effect of presplitting on the generation of a smooth wall in continuum and jointed rock mass. The 2D distinct element code was used to simulate the presplitting in a rock slope. The blast load history as a function of time was applied to the inner wall of each blasthole. Important parameters that were considered in the analysis were stress tensor and fracturing pattern. The blast loading magnitude and blasthole spacing and jointing pattern were found to be very significant in the final results.

Modeling stress wave propagation in rocks by distinct lattice spring model

In this paper, the ability of the distinct lattice spring model （DLSM） for modeling stress wave propagation in rocks was fully investigated. The influence of particle size on simulation of different types of stress waves （e.g. one-dimensional （1D） P-wave, 1D S-wave and two-dimensional （2D） cylindrical wave） was studied through comparing results predicted by the DLSM with different mesh ratios （It） and those obtained from the corresponding analytical solutions. Suggested values of lr were obtained for modeling these stress waves accurately. Moreover, the weak material layer method and virtual joint plane method were used to model P-wave and S-wave propagating through a single discontinuity. The results were compared with the classical analytical solutions, indicating that the virtual joint plane method can give better results and is recommended. Finally, some remarks of the DLSM on modeling of stress wave propagation in rocks were provided.

System dynamics model of the support-surrounding rock system in fully mechanized mining with large mining height face and its application

Fully mechanized mining with large mining height(FMMLMH)is widely used in thick coal seam mining face for its higher recovery ratio,especially where the thickness is less than 7.0 m.However,because of the great mining height and intense rock pressure,the coal wall rib spalling,roof falling and the instability of support occur more likely in FMMLMH working face,and the above three types of disasters interact with each other with complicated relationships.In order to get the relationship between each two of coal wall,roof,floor and support,and reduce the occurrence probability of the three types of disasters,we established the system dynamics(SD)model of the support-surrounding rock system which is composed of"coal wall-roof-floor-support"(CW-R-F-S)in a FMMLMH working face based on the condition of No.15104 working face in Sijiazhuang coal mine.With the software of Vensim,we also simulated the interaction process between each two factors of roof,floor,coal wall and the support.The results show that the SD model of"CW-R-F-S"system can reveal the complicated and interactive relationship clearly between the support and surrounding rock in the FMMLMH working face.By increasing the advancing speed of working face,the support resistance or the length of support guard,or by decreasing the tipto-face distance,the stability of"CW-R-F-S"system will be higher and the happening probability of the disasters such as coal wall rib spalling,roof falling or the instability of support will be lower.These research findings have been testified in field application in No.15104 working face,which can provide a new approach for researching the interaction relationship of support and surrounding rock.

The three-dimension model for the rock-breaking mechanism of disc cutter andanalysis of rock-breaking forces

To study the rock deformation with three- dimensional model under rolling forces of disc cutter, by car- rying out the circular-grooving test with disc cutter rolling around on the rock, the rock mechanical behavior under rolling disc cutter is studied, the mechanical model of disc cutter rolling around the groove is established, and the the- ory of single-point and double-angle variables is proposed. Based on this theory, the physics equations and geometric equations of rock mechanical behavior under disc cutters of tunnel boring machine （TBM） are studied, and then the bal- ance equations of interactive forces between disc cutter and rock are established. Accordingly, formulas about normal force, rolling force and side force of a disc cutter are de- rived, and their validity is studied by tests. Therefore, a new method and theory is proposed to study rock- breaking mech- anism of disc cutters.