地下洞室围岩脆性破坏时的应力特征研究

在高应力作用下,岩爆、钻孔崩落、片帮都是地下空间硬脆围岩中常见的破坏现象,这三类现象本质上均可归于完整岩体的脆性破坏,它们分别反映了高应力作用下完整岩体不同的破坏程度。通过对前人关于岩爆判据、钻孔崩落判据和片帮应力强度比判据研究成果的类比分析可知,这些脆性破坏现象在破坏时具备相同的应力背景条件。脆性破坏的应力条件可以用地下空间周边切向最大应力与岩石单轴抗压强度之比（ / ）或者工程区最大主应力与岩石单轴抗压强度之比（ / ）来描述,两种指标本质上反映了相同的应力背景条件。对于 / , / = 0.4 ± 0.1是发生脆性破坏的应力临界条件;对于 / , / = 0.15 ± 0.05是发生脆性破坏的应力临界条件。大量的工程实例和基于Hoek-Brown强度准则的力学分析也证明了这一背景条件的正确性。这里两种指标都取了一个范围,主要是由于不同的岩体分级、岩性和工程地质条件会对指标的界定产生较为显著的影响。     

Inference of in situ stress from thermoporoelastic borehole breakouts based on artificial neural network

The determination of in situ stresses is very important in petroleum engineering. Hydraulic fracturing is a widely accepted technique for the determination of in situ stresses nowadays. Unfortunately, the hydraulic fracturing test is time-consuming and expensive. Taking advantage of the shape of borehole breakouts measured from widely available caliper and image logs to determine in situ stress in petroleum engineering is highly attractive. By finite element modeling of borehole breakouts considering thermoporoelasticity, the authors simulate the process of borehole breakouts in terms of initiation, development, and stabilization under Mogi-Coulomb criterion and end up with the shape of borehole breakouts. Artificial neural network provides such a tool to establish the relationship between in situ stress and shape of borehole breakouts, which can be used to determine in situ stress based on different shape of borehole breakouts by inverse analysis. In this paper, two steps are taken to determine in situ stress by inverse analysis. First, sets of finite element modeling provide sets of data on in situ stress and borehole breakout measures considering the influence of drilling fluid temperature and pore pressure, which will be used to train an artificial neural network that can eventually represent the relationship between the in situ stress and borehole breakout measures. Second, for a given measure of borehole breakouts in a certain drilling fluid temperature, the trained artificial neural network will be used to predict the corresponding in situ stress. Results of numerical experiments show that the inverse analysis based on finite element modeling of borehole breakouts and artificial neural network is a promising method to determine in situ stress.     

中国大陆科学钻主孔现今地应力状态

用钻孔崩落法确定了中国大陆科学钻探主钻孔5 047 m深度以上的现今地应力状态.由钻孔声波成像测井资料发现, 科学钻主钻孔在1 200 m深度以下出现了钻孔崩落现象.我们从1 216~5 047 m的深度范围内采集了143个钻孔成像测井图象资料, 对钻孔崩落椭圆长轴方位进行了统计, 结果表明崩落椭圆长轴平均方位为319.5°±3.5°, 最大水平主应力方位平均为49.5°±3.5°.利用崩落形状要素(崩落深度和崩落宽度) 以及岩石的内聚力和内摩擦角, 估算了1 269 m至5 047 m范围内52个深度上的最大和最小水平主应力的大小.结果表明, 在浅处1 216 m深度, 最大水平主应力为42 MPa, 最小水平主应力为30.3 MPa; 在深处5 000 mm深度, 最大水平主应力为160.5 MPa, 最小水平主应力为120 MPa; 地应力随深度近于线性增加.据岩石密度测井资料计算了各个深度上静负载应力.3个主应力的大小和方向反映出科学钻主孔位置的应力场处于走滑应力状态, 与临近地区地震震源机制解和其他方法得到的应力场一致.利用声发射法对岩心试件进行了声发射测量, 得到了最大水平主应力幅值, 并与崩落法测量结果进行了对比, 两者十分一致.      

Simulation of drilling‐induced compaction bands using discrete element method

Rock failure is observed around boreholes often with certain types of failure zones, which are called breakouts. Laboratory‐scale drilling tests in some high‐porosity quartz‐rich sandstone have shown breakouts in the form of narrow localized compacted zones in the minimum horizontal stress direction. They are called fracture‐like breakouts. Such compaction bands may affect hydrocarbon extraction by forming barriers that inhibit fluid flow and may also be a source of sand production. This paper presents the results of numerical simulations of borehole breakouts using 3D discrete element method to investigate the mechanism of the fracture‐like breakouts and to identify the role of far‐field stresses on the breakout dimensions. The numerical tool was first verified against analytical solutions. It was then utilized to investigate the failure mechanism and breakout geometry for drilled cubic rock samples of Castlegate sandstone subjected to different pre‐existing far‐field stresses. Results show that failure occurs in the zones of the highest concentration of tangential stress around the borehole. It is concluded that fracture‐like breakout develops as a result of a nondilatant failure mechanism consisting of localized grain debonding and repacking and grain crushing that lead to the formation of a compaction band in the minimum horizontal stress direction. In addition, it is found that the length of fracture‐like breakouts depends on both the mean stress and stress anisotropy. However, the width of the breakout is not significantly changed by the far‐field stresses. Copyright © 2013 John Wiley & Sons, Ltd.     

Stress Concentration at the Bottom of a Borehole and its Effect on Borehole Breakout Formation

Summary Stress concentration at the bottom of a borehole due to the corners with small radius of curvature in an axial section and its effect on the azimuth of breakout was studied. To this end, a 3-D finite element analysis was conducted and the stress around the borehole was examined for boreholes arbitrarily oriented to three principal axes of remote stress. Results show that, in the case of high strength rock, compressive failure resulting in spalling of a borehole may occur only at the bottom of the borehole. The spalling can occur continuously with drilling, and results in continuous spalling with depth, i.e., a breakout. This type of breakout tends to form on one side of the borehole and its orientation is approximately perpendicular to the orientation of standard breakouts, inferred from the stress concentration due to the cylindrical shape of the borehole.     

The influence of the intermediate principal stress on rock failure behaviour: A numerical study

The influence of the intermediate principal stress on rock strength has been studied comprehensively by previous researchers. However, the reason why rock strength firstly increases and subsequently decreases with the increase of intermediate principal stress is still unclear. In this paper, the mechanism of the intermediate principal stress effect on rock failure behaviour is revealed through a numerical method using the EPCA3D system (Elasto-Plastic Cellular Automaton). In this study, both homogeneous and heterogeneous rocks are considered. The heterogeneity of a rock specimen is modelled by introducing Weibull's statistical distribution. Two criteria, i.e. the Drucker–Prager and Mohr–Coulomb models, are used to determine whether a meso-scopic element in the rock specimen is in a failure state or not during the polyaxial stress loading process. The EPCA3D simulation reproduces the typical phenomenon of the intermediate principal stress effect that occurs in some rock experiments. By studying the EPCA3D simulated acoustic emission and complete stress–strain curves illustrating failure initiation, propagation and coalescence in the failure process of rocks, the essence of the intermediate principal stress effect is tracked. It is concluded that the heterogeneous stress distribution induced by the natural heterogeneity of rocks and the effect of the loading platen are two of the reasons producing the intermediate stress effect. Studies indicate that a moderate intermediate principal stress delays the onset of local failure, which in turn leads to an increase in the rock strength. However, once the intermediate principal stress reaches a certain value, local failure will be formed through the application of the intermediate principal stress. It is the number of failed elements in the pre-peak region that determines whether the rock strength decreases or not. The extent of rock strength reduction when the intermediate principal stress reaches a certain value is lessened with the increase in the minimum principal stress.     

复合岩盐层井眼蠕变缩径的数值模拟

盐岩在深部高温高压环境下具有很强的蠕变特性,在岩盐层钻井经常出现井眼缩径、卡钻等复杂事故。目前关于岩盐层井段安全钻进泥浆密度的计算方法主要适用于均匀地应力的情况,难以在非均匀地应力条件下应用。为此研究建立了三向地应力作用下盐岩和砂泥岩相互交错的复合岩盐层井眼蠕变缩径的三维有限元模型,并以塔里木油田羊塔克地区深井为例,研究了非均匀地应力条件下,复合岩盐层井眼随时间变化的蠕变缩径规律,得出了一定泥浆密度下岩盐层井眼在不同时刻的井径值。计算结果表明,安全泥浆密度有限元模拟值与羊塔克地区岩盐层井段实际采用的安全泥浆密度非常接近,验证了所建数值模型的合理性和计算结果的可靠性,研究成果为岩盐层井段的安全钻进提供了技术支撑。      

深部地质钻探钻遇弱面地层井孔围岩稳定分析

深部地层地质条件复杂,随钻探向纵深发展,钻遇弱面地层导致井孔围岩失稳将严重影响工程进程。本文总结分析具有弱面地层岩石物理力学特性,探讨强度各向异性孔壁围岩稳定模型的适用性;考虑流体向弱面渗透存在渗透各向异性和传热各向异性,发展了具有弱面地层井孔围岩稳定理论。研究发现,某些强度各向异性模型在特定情况下可能失去其适用性;在深部高温环境下忽略弱面地层渗透和热传的各向异性特征可能会导致井孔围岩压力、破坏区域以及坍塌压力产生较大误差。     

Numerical modeling of plastic deformation and failure around a wellbore in compaction and dilation modes

In many wellbore stability analyses, the ability to forecast both the occurrence and extent of plastic deformation and failure hinges upon a fundamental understanding of deformation mode and failure mechanism in the reservoir rock. This study focuses on analyzing plastic zones, localized deformations, and failures around a borehole drilled overbalanced or underbalanced through a highly porous rock formation. Based on several laboratory experiments, porous rocks are prone to deform under both shear-induced dilation and shear-enhanced compaction mechanisms depending on the stress state. The shapes of the deformation and failure patterns around the borehole are shown, depending on the initial stress state and the local stress paths. The inquiry of the local stress paths in the near-wellbore zone facilitates the understanding of the reasons for different types of failure mechanisms, including the mixed-mode and the plastic deformation structures. The modification of the 2D plane strain condition by imitating third stress in the numerical scheme helps us bring the stress paths closer to the real state of loading conditions. Our modeling reveals that the transition from isotropic to anisotropic stress state is accompanied by an increase in the deviatoric part of effective shear tensor that leads to the development of inelastic deformation, degradation, and subsequent rock failure. Particular interest is devoted to the modeling of strain localization especially in compaction mode around a wellbore and computing the amount of stress concentration at the tips of dog-eared breakouts. Stress concentration can result in a change in irreversible deformation mode from dilatancy to compaction, elucidating the formation of the shear-enhanced compaction phenomenon at the failure tips in the direction of the minimum horizontal stress.     

Borehole instability in high-porosity Berea sandstone and factors affecting dimensions and shape of fracture-like breakouts

Stress-induced breakouts in vertical boreholes are failure zones caused by excessive compressive stress concentration at the borehole wall along the springline of the least horizontal far-field stress. Wellbores are sometimes drilled into aquifers or oil reservoirs that are weak, poorly consolidated, and highly porous sandstone formations, which are often conducive to breakout formation. Breakouts are an expression of borehole instability and a potential source of sand production. On the other hand, the breakout phenomenon can be used advantageously in obtaining an estimate of the in situ stress condition. The average orientation of breakouts, as identified by borehole geophysical logging, is a reliable indicator of in situ stress directions. It has also been suggested that breakout dimensions could potentially be used as indicators of in situ stress magnitudes. The research reported here has concentrated on the unique type of breakouts observed for the first time in high-porosity Berea sandstone. Drilling experiments in rock blocks subjected to critical far-field true triaxial stress regimes, simulating in situ conditions, induced breakouts that were unlike the ‘dog-ear’ ones previously observed in granites, limestones, and low-porosity sandstones. The newly observed breakouts were thin, tabular, and very long, resembling fractures that counterintuitively extended perpendicular to the maximum principal stress. We found that a narrow zone ahead of a fracture-like breakout tip underwent apparent localized grain debonding and compaction. In the field, such zones have been termed ‘compaction bands’, and are a source of concern because in oil fields and aquifers they constitute curtains of low permeability that can impede the normal flow of oil or water. In order to determine whether a correlation exists between fracture-like breakouts and in situ stress, we conducted several series of tests in which the minimum horizontal and vertical stresses were held constant and the maximum horizontal stress (σ_H) was increased from test to test. These tests showed strong dependence of the breakout length on far-field stress, signaling that potentially the ability to assess fracture-like breakout length in the field could be used to estimate in situ stress magnitudes in conjunction with other indicators. Another series of tests revealed that breakout length increased substantially when borehole diameter was enlarged. This result suggested that in the field, where wellbore size is considerably larger, fracture-like breakout could extend to sizable distances, creating a sand production hazard. Two series of tests, one to evaluate the effect of drill-bit penetration rate, and the other to verify the drilling-fluid flow rate effect on breakout formation and dimensions yielded inconsistent results and showed no unique trends. Remarkably, fracture-like breakouts maintained a consistent narrow width of about 5–10 grain diameters, irrespective of the test conditions. This characteristic supports the suggestion that fracture-like breakouts are emptied compaction bands.     

Borehole failure analysis in a sandstone under anisotropic stresses

Borehole failure under anisotropic stresses in a sandstone is analyze numerically for various borehole sizes using a nonlinear elastic–plastic constitutive model for a Cosserat continuum. Borehole failure is identified as macroscopic failure of the borehole through the development of shear bands and breakouts. The results compare well both qualitatively and quantitatively with experimental results from polyaxial tests on Red Wildmoor sandstone. They show that the hole size effect of the borehole failure strength is independent of the far‐field stress anisotropy and follows a ? power law of the hole size. A similar scale effect equation with a ? power law is proposed for the scale effect of the maximum plastic shear strain at failure. This equation can be useful for better predicting hole‐size‐dependent failure with standard codes based on classical continua. The effect of stress anisotropy on the borehole failure stress is found to be independent of the hole size. The failure stress decreases linearly to 40% as the stress anisotropy increases. However, the maximum plastic shear strain at failure is stress anisotropy independent and therefore the critical plastic shear strain for failure is only hole‐size dependent. Copyright © 2009 John Wiley & Sons, Ltd.     

Compaction bands induced by borehole drilling

Drilling experiments in rock blocks subjected to pre-existing true triaxial far-field stresses simulating real in situ conditions often result in localized failure around the created borehole, which brings about the formation of borehole breakouts. In weakly bonded quartz-rich porous sandstones breakouts take the form of narrow tabular (slot-like) openings extending along a plane perpendicular to the maximum applied-stress direction. Scanning electron microscopes images of failed boreholes strongly suggest that these breakouts are compaction bands that have been emptied to different extents. The bands form as a result of the stress concentration accompanying the creation of the borehole. The evacuation of the compaction bands is brought about by the circulating drilling fluid flushing out debonded and often fragmented grains from within these bands (Haimson and co-workers, 2003–2007). The objective of this paper is to predict the conditions under which compaction bands are formed around boreholes. To this end, a new analytical model is formulated that enables prediction of the stress field around emptied and filled compaction bands, the various factors affecting the breakouts lengths, and their final length. Good agreement of the developed analytical model with experimental results obtained by Haimson and co-workers (Haimson and Klaetsch in Rock physics and geomechanics in the study of reservoirs and repositories, vol 284, pp 89–105, 2007; Haimson and Kovachich in Eng Geol 69:219–231, 2003; Klaetsch and Haimson in Mining and tunneling innovation and opportunity, University of Toronto press, pp 1365–1371, 2002; Sheets and Haimson in Proceedings, paper ARMA/NARMS 04-484, 2004) is demonstrated. The presented study is of practical relevance: boreholes are often drilled deep into weak porous sandstone formations for the purpose of extracting oil and gas, and the question of borehole stability is crucial. In addition, borehole breakouts are often used to estimate the state of stress in the Earth’s crust, and our new formulation will help improve these estimates.     

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.     

利用钻孔崩落数据再认识白鹤滩右岸地应力场特征

白鹤滩水电站是仅次于三峡水电站的第二大水电站,位于中国西南地区川滇菱形块体内的金沙江上。通常地壳应力状态是影响地下工程安全的重要地质因素,对地下硐室稳定性分析具有重要意义。在水电站右岸厂房建设过程中,为了水电站的长期安全运营,采用超声波井下电视录井测试系统对白鹤滩右岸厂房锚固洞内7处钻孔进行测试,基于钻孔崩落数据计算了现今白鹤滩右岸厂房区域上方工程岩体的主应力方向。研究结果表明:白鹤滩右岸厂房区域最大水平主应力（S_H）方向为北北东—南南西方向,主要受到构造应力、自重应力、河流剥蚀作用以及岸坡卸荷作用的共同影响,属于局部构造应力场。      

Effect of fluid pressure on rock compressive failure in a nearly impermeable crystalline rock: Implication on mechanism of borehole breakouts

We conducted laboratory true triaxial experiments in the nearly impermeable Pohang rhyolite to investigate failure mechanisms under ‘dry’ and ‘wet’ rock conditions. Under ‘dry’ conditions prismatic specimens were jacketed all around to prevent confining fluid penetration. Under ‘wet’ conditions one pair of the specimen faces was left unjacketed and in direct contact with the confining fluid (kerosene) applying the least principal stress in an attempt to simulate the case of an unlined borehole wall. In both testing setups the true triaxial compressive strength for a given least principal stress increases significantly as the intermediate principal stress rises. The unjacketed rhyolite strength is, however, only 60 to 85% of the strength under dry conditions, depending on the magnitude of the intermediate principal stress. In dry rhyolite the failure process begins upon dilatancy onset, followed by microcrack localization, and ending in a steeply dipping shear fracture. On the other hand, brittle fracture in wet specimens occurs almost immediately after the onset of dilatancy by the development of one or more through-going extensile fractures subparallel and adjacent to one of the unjacketed faces, resembling the extensile cracks leading to borehole breakouts in crystalline rocks. We infer that upon dilatancy the confining fluid intrudes and quickly propagates newly opened stress-induced microcracks subparallel to the unjacketed faces, leading to ‘early’ failure.     

Digital Image Based Approach for Three-Dimensional Mechanical Analysis of Heterogeneous Rocks

Summary This paper presents a digital image based approach for three-dimensional (3-D) numerical simulation and failure analysis of rocks by taking into account the actual 3-D heterogeneity. Digital image techniques are adopted to extract two-dimensional (2-D) material heterogeneity from material surface images. The 2-D image mesostructures are further extrapolated to 3-D cuboid mesostructures by assuming the material surface as a representation of the inner material heterogeneity within a very small depth. The iterative milling and scanning system is set up to generate the 3-D rock mesostructures. A Hong Kong granite specimen is used as an example to demonstrate the procedure of 3-D mesostructure establishment. The mechanical responses and failure process under the conventional Brazilian tensile test condition are examined through numerical analyses. The stress distribution, crack propagation process and failure model of heterogeneous material cases are simulated with a finite difference software. The numerical results indicate that material heterogeneity plays an important role in determining the failure behavior of rocks under external loading.     

三维梁-颗粒模型在岩石材料破坏模拟中的应用

用三维梁-颗粒模型BPM3D(beam-particlemodelinthreedimensions)对岩石类非均质脆性材料的力学性质和破坏过程进行了数值模拟。梁-颗粒模型是在离散单元法基础上,结合有限单元法中的网格模型提出的用于模拟岩石类材料损伤破坏过程的数值模型。在模型中,材料在细观层次上被离散为颗粒单元集合体,相邻颗粒单元由有限单元法中的弹脆性梁单元联结。梁单元的力学性质均按韦伯(Weibull)分布随机赋值,以模拟岩石类材料力学参数的空间变异性。材料内部裂纹通过断开梁单元来模拟。通过自动生成的非均质材料模型对岩石类材料的破坏机理进行研究。岩石类非均质脆性材料在单轴压缩状态下破坏过程细观数值模拟结果显示,岩石材料宏观破坏是由于其内部细观裂纹产生、扩展、贯通的结果。通过数值模拟结果之间的对比分析,揭示出岩石试样宏观破坏模式随细观层次上韦伯分布参数的变化而不同。与实际矿柱破坏形态的对比分析表明了模型的适用性。根据数值模拟结果对岩石类非均质材料的破坏机理进行了探讨。     

基于数字图像技术的岩土材料有限元-离散元分析

将数字图像技术与Munjiza提出的有限元法-离散元法（FDEM）耦合分析方法结合,研制了可表征岩石真实非均质性的FDEM分析系统,为从细观角度对岩体进行建模以及研究岩体的破裂机制提供了新途径。该系统借助数字图像技术从岩石断面的图像获取岩石材料的真实细观结构,借助成熟的网格剖分技术,将其映射到FDEM计算网格中,从而克服了原有FDEM在考虑材料非均质性所存在的不足。利用该系统,进行了巴西圆盘劈裂试验的数值模拟,再现了花岗岩在荷载作用下的真实破裂过程。数值模拟结果表明,考虑非均质性的岩样应力分布呈现出非对称性,岩石的细观结构对岩石中裂纹的扩展及应力分布有重要影响。     

A Combined Three-Dimensional Geological-Geostatistical-Numerical Model of Underground Excavations in Rock

Summary. This paper exploits geological and borehole geotechnical data obtained in the exploratory phase of a tunneling project to investigate in a first place if the kriging interpolation scheme may effectively reproduce the spatial variability of rock mass quality (Rock Mass Rating, RMR) in the vicinity of tunnels. For this purpose a quick solver in Fortran has been developed that performs variography analysis of 3D spatial data, fast kriging estimations of RMR between borehole sampling locations at the centroids of the elements of the numerical model, and model validation. For the purpose of an integrated underground excavation design, a step further is made by incorporating into the 3D mechanical numerical model of the rock mass, the three-dimensional (3D) solid geological model, thus coupling the geology with the ground (geotechnical) model (i.e. each element of the numerical model is assigned a geological material). The mechanical properties of each finite difference cell (or Representative Elementary Volume) of the ground model were then prescribed according to its geological type, the spatial heterogeneity of the rock mass expressed quantitatively with the kriging model, and the upscaling calculations of the mechanical properties of the intact rocks determined in the laboratory, based on the size-effect (strength dependence on size) and Damage Theory. Furthermore, a preliminary numerical simulation of the advance of unsupported tunnels in the model of the heterogeneous rock mass was performed for illustration purposes.     

井底应力场对气体钻井井斜的影响

建立了油气井钻井过程中的井底力学模型,该力学模型考虑了井眼内压力、岩石节理、井眼深度、地应力、岩石材料特性、井眼直径等影响因素,并讨论了各个因素对井底应力场的影响,阐明了气体钻井时井底应力场分布变化对井眼井斜的影响。气体钻井时井眼内的压力远小于井眼周围的围压,井底会出现很大的拉应力,岩石容易破坏,因此,气体钻井时钻井速度远远大于钻井液钻井,但拉应力越大,井眼与岩石节理面交接处产生的应力集中越大,应力集中处岩石更容易破坏,也更易产生井斜。气体钻井时井底的应力状态对地层倾角、材料特性、井眼直径、井眼深度及初始地应力的变化较钻井液钻井时更敏感。