基于水压致裂法的岩石抗拉强度研究

为了揭示水压致裂法和巴西劈裂法测量岩石抗拉强度的关系,开展了理论和现场试验研究。基于经典的水压致裂法理论,推导了不同围压下钻孔破裂压力和抗拉强度。利用断裂力学理论建立了水压致裂法和巴西劈裂法测得抗拉强度的关系。利用预制切槽方法模拟天然裂纹,对水力裂缝的起裂压力进行了研究。结果表明：围压为最大主应力等于3倍最小主应力测得的抗拉强度大于围压为0测得的抗拉强度;水压致裂法和巴西劈裂法测量抗拉强度关系与应力场、裂纹长度、断裂韧度3个变量有关;通过在晋城矿区王台铺矿的预制切槽试验,运用断裂力学建立的抗拉强度计算式更为符合现场实际。研究结果可为坚硬难垮落顶板预制切槽的水力压裂设计提供参考。     

单孔岩样水压致裂的数值分析

水压致裂是改变岩体结构的一种天然行为和人为手段。采用F－RFPA2D软件,对水压致裂过程、裂缝扩展形态及注水孔形状和大小、应力条件和岩样强度等影响因素进行了研究。将开始出现声发射的水压称为微裂压力,将声发射急剧增多、裂缝非稳定扩展直至岩样破坏的水压称为破裂压力。岩样尺寸一定时,微裂压力和破裂压力随内孔面积增加而降低,方形孔岩样的微裂压力和破裂压力均小于同面积的圆形孔。微裂压力和破裂压力随围压或岩样强度增加而增加,且其差值随岩样强度增加而增加,理论破裂压力与模拟值趋势基本一致。方形孔的宏观裂纹起裂位置多在角点附近,而圆形孔比较随机。无围压时,宏观裂纹的延伸方向随机;有围压时,宏观裂纹扩展方向大致与主应力方向一致,且沿较大主应力方向的宏观裂纹扩展至岩样破坏,较小主应力方向宏观裂纹不完全发育。研究结果对水压致裂试验和工程实践有一定参考意义。     

煤层水压致裂后煤岩应力解析

水压致裂后煤岩应力分布规律对水压致裂防冲效果起关键性作用。采用理论研究方法得出高压注水压致裂后及卸水后水区和气区的孔隙、瓦斯压力和煤体应力解析解。结果表明,致裂后水区孔隙压力沿径向变化不大,与注水压力接近;气区瓦斯压力沿径向呈递减趋势;在水区外围一定范围内形成瓦斯压力升高区;水区煤体环向应力将会减小,直到变为拉应力;气区煤体径向应力沿径向递减。卸水后水区孔隙压力、煤体径向应力沿径向呈递增趋势;气区煤体径向应力沿径向呈递增趋势,趋近于原始煤体应力;气区煤体环向应力沿径向呈递减趋势;气区孔隙压力沿径向呈递减趋势。这为煤层水压致裂预防冲击地压提供理论基础。      

水压致裂试验过程中诱发的声发射

为了研究水压致裂过程中裂缝的扩展机制,在北京房山花岗岩体中开展了大型水压致裂试验。试验是在一个深 301 m的岩石新鲜完整的FR钻孔中进行的,其周围半径约40 m的范围内布置了4个120 m深的声发射观测孔。在FR孔中深度60 m至140 m的范围内选取了7段没有天然节理的部位进行了水压致裂。在所有的试验中,发现水压致裂形成的裂缝通过天然节理与AE孔连通,只在深度118.5 m进行的试验,观测到由水压致裂产生的AE事件。由声发射的震源机制解得到的P轴和T轴的方向与水压致裂应力测量的方向一致。     

极地深部冰层取心钻探孔壁水压致裂研究

为了维护极地深部冰层取心钻探工作中钻孔的稳定,避免孔内事故的发生,冰层孔壁的水压致裂问题是亟待解决的重要科学问题之一。本文在深入分析国内外冰层钻探资料的基础上,结合冰盖动力学相关理论,计算得出了钻孔所在区域冰层的密度、温度及内部应力随深度的变化规律。在此基础上,结合油气资源勘探水力压裂技术与冰层钻进钻井液等相关理论,建立了合理的孔壁压差计算方法,深入探讨了适用于深部冰层钻探孔壁水压致裂机理。研制了冰层钻孔水压致裂模拟实验装置,可分析研究不同围压条件下冰样的脆性变形机理。结合我国Dome A深冰心钻探工程实际,提出钻孔可能发生水压致裂的深度区域及孔壁所需的起始裂纹长度判定,以期为后续的安全高效冰层钻进提供重要的理论依据。     

一种用于水压致裂试验的水压供给装置研制及应用

为研究水压致裂过程中水力裂缝的起裂、扩展规律及机制,研制了一种新型水压供给控制装置。该装置由水泵、增压泵、空气压缩机和导水钢垫板等组成。该装置的主要特点是易于与目前广泛存在的压力伺服机相配合,共同完成水压致裂试验,具有操作简便、效率高、成本低、组装方便等优点,使水压致裂试验更易大面积推广。利用该试验装置对类岩石试块进行了水压致裂试验,研究了不同围压下含预置裂隙试件的水力裂缝扩展规律,并与无预置裂隙试件水力裂缝扩展规律进行了对比。试验结果表明,对于含预置裂隙试块,当水平应力差较小时,预置裂隙裂尖对水力裂缝扩展影响较大,水力裂缝起裂后向预置裂隙裂尖方向扩展;随着水平应力差的增大,水力裂缝向最大水平应力方向偏转。试验证明该装置可靠性高,试验过程稳定。该研究对水压致裂试验具有促进意义。     

热应力影响下干热岩水压致裂数值模拟

干热岩水压致裂过程中低温诱导热应力与注入水压共同影响裂缝的萌生与扩展。首先通过THM耦合分析了低温压裂液注入过程中注入水压与热应力的相互作用及其对裂缝萌生的影响,随后建立描述岩石细观结构的THMD耦合模型对热应力影响下高温岩石水压致裂过程进行初探。结果表明：低温压裂液注入高温岩石产生的热应力包括岩石自身温度梯度形成的热应力与岩石颗粒非均匀膨胀导致的热应力,并在井筒周围呈现为拉应力。高注入压力将抑制热应力导致的多裂缝萌生,井筒附近热应力的存在对注入压力也具有削弱作用。基岩温度升高,裂缝萌生阶段更多裂缝在井筒附近起裂,缝网沿最大地应力方向的扩展速度减慢,但改造规模增加,同时多裂缝的存在也使得裂缝延伸压力增加。     

倾斜煤层水力压裂起裂压力计算模型及判断准则

针对目前倾斜煤层起裂机制不明确,导致水力压裂时盲目升高压力或增大注水量来增加煤层透气性的问题,根据最大拉应力理论,分析真实环境下倾斜煤层压裂钻孔周围应力状态,建立压裂钻孔周围煤岩体起裂压力计算模型及判断准则,并在重庆松藻煤电公司同华煤矿进行了验证,结果表明,根据实测压裂区域地应力状态,起裂压力随煤层倾角增大而增大,钻孔起裂位置随煤层倾角增大逐渐向走向方向偏转;现场试验起裂压力与理论计算相符,随煤层倾角增大而增大,从而验证了计算模型的正确性。     

水压致裂法测量裂隙岩体的地应力研究

通过对水压致裂法形成的钻孔印模分析,运用Taylor公式和最小二乘法迭代算法,计算纵向和横向裂隙岩体地应力大小及其误差。首先确定大地、裂隙和钻孔坐标系,根据印模曲线确定裂隙方位;然后用Taylor公式和最小二乘法迭代算法,计算应力分量;最后计算出主应力大小及其误差。将此方法应用于雪峰山隧址区ZK6钻孔的地应力场计算,得出其最大水平应力为13 MPa,方位113°左右,这与地质构造法实测的应力场相吻合,说明这一方法合理可行,具有较好的实用价值。      

基于曲线滑裂面的挡墙主动土压力分析

土体滑裂面形状对挡土墙主动土压力有重要影响。以无黏性填土挡墙为研究对象,假设在考虑土拱效应时,极限状态下墙后土体的滑裂面为曲线,基于水平微分单元法推导出平动模式下挡土墙主动土压力的分布。首先将计算与模型试验结果及已有理论研究成果进行对比分析,验证了方法的可靠性;其次,研究土体内摩擦角和墙-土摩擦角对主动土压力分布、合力大小和作用点高度的影响。结果表明:基于曲线滑裂面假设得到的滑动楔体范围略大于采用直线滑裂面的假设;对于不同高度的挡土墙,建议的计算结果与模型试验结果更为符合;对于不同的土体内摩擦角和墙-土摩擦角,土压力的分布形式和合力作用点与Paik解较为接近,但合力略大于Paik解。     

基于离散元法的水力压裂数值模拟

不同条件下水压裂隙的发展特性对有效开采页岩气具有重要的指导作用。针对岩体在微观上为颗粒和孔隙的结构系统,提出离散元水力压裂数值模拟方法,离散元能量转化和能量守恒计算方法,建立了相应的三维离散元模型。采用自主研发的三维离散元模拟软件Mat DEM3D,通过控制模型的竖向应变与颗粒直径,来模拟地层中的应力与压裂速率的变化。模拟结果表明:(1)水力压裂产生裂隙的数量和方向受岩石的各向异性,压力状态和变化速率所影响。(2)裂隙在压缩波传播时发展,当水压力高速增加时,诱发的裂隙数量增多,并且有效能量(断裂热)百分比也随之增加,压裂作用也变得更明显。(3)当竖向应变为零时,50%的裂隙呈垂直状态,当竖向应变为-1×10-4时,裂隙趋于沿着最大压力方向发展,竖向裂隙的百分率增大。数值模拟和能量分析为定量地研究岩石水力压裂过程提供了一个新的方法。     

Microcrack-based coupled damage and flow modeling of fracturing evolution in permeable brittle rocks

Microcracks in brittle rocks affect not only the local mechanical properties, but also the poroelastic behavior and permeability. A continuum coupled hydro-mechanical modeling approach is presented using a two-scale conceptual model representing realistic rock material containing micro-fractures. This approach combines a microcrack-based continuous damage model within generalized Biot poroelasticity, in which the tensors of macroscopic elastic stiffness, Biot effective stress coefficient and of overall permeability are directly related to microcrack growth. Heterogeneity in both mechanical and hydraulic properties evolves from an initially random distribution of damage to produce localized failure and fluid transmission. A significant advantage of the approach is the ability to accurately predict the evolution of realistic fracturing and associated fluid flow in permeable rocks where pre-existing fractures exert significant control. The model is validated for biaxial failure of rock in compression and replicates typical pre- and post-peak strength metrics of stress drop, AE event counts, permeability evolution and failure modes. The model is applied to the simulation of hydraulic fracturing in permeable rocks to examine the effects of heterogeneities, permeability and borehole pressurization rate on the initiation of fracturing. The results indicate that more homogenous rocks require higher hydraulic pressure to initiate fracturing and breakdown. Moreover, both the fracturing initiation pressure and breakdown pressure decrease with permeability but increase with borehole pressurization rate, and the upper and lower limit of the initiation pressure are seen to be given by the impermeable (Hubbert–Willis) and permeable (Haimson–Fairhurst) borehole wall solutions, respectively. The numerical results are shown to be in good agreement with the experimental observations and theoretical results. This coupled damage and flow modeling approach provides an alternative way to solve a variety of complicated hydro-mechanical problems in practical rock engineering with the process coupling strictly enforced.     

基于孔壁应变发展规律的压裂孔三阶段起裂特征试验研究

利用自行研制的煤岩体水力压裂试验系统,开展了配比型煤与原煤水力压裂试验,测试并分析了水力压裂过程中压裂孔孔壁应变-水压曲线,并基于孔壁应变的发展规律,分析了压裂孔的三阶段起裂特征。结果表明,在压裂孔起裂过程中,钻孔孔壁呈现拉伸与压缩应变两种类型,并呈现拉伸破裂区与压缩变形区,其中压缩型应变具有较好的可恢复性,其应变恢复比远大于拉伸型应变;钻孔起裂过程分为3个阶段,即水气作用诱导微损伤形成阶段,孔壁内形成气流通道并产生初始损伤;局部损伤带形成阶段,孔壁形成拉伸破裂区和压缩变形区;试件失稳破坏阶段,裂缝不断延伸直至试件破裂,拉伸破裂区依然保持拉伸变形并较好地保持残余变形,而压缩变形区则由于作用力转向而得到一定程度恢复。研究成果对于揭示钻孔起裂行为及能量的演化规律具有重要理论意义。     

Hydraulic fracturing modeling using the discontinuous deformation analysis (DDA) method

Recent attempts to couple discontinuous deformation analysis (DDA) blocks with finite element meshes have been regarded as fruitful. The hybrid model has been proven to be reliable in modeling continuum–discontinuum problems. In this paper, a hydraulic crack initiation–propagation extension complete with a coupled hydro-mechanical analysis algorithm is introduced to the hybrid model for the simulation of fracturing problems initiated by hydraulic pressure. Several simulation cases, of which the results appear to be consistent with existing theories of hydraulic fracturing mechanisms, are presented to demonstrate the capability of the proposed algorithm in modeling hydraulic fracturing processes.     

Numerical Modeling of Natural Fracture Distributions in Shale

The production efficiency of shale gas is affected by the interaction between hydraulic and natural fractures. This study presents a simulation of natural fractures in shale reservoirs, based on a discrete fracture network (DFN) method for hydraulic fracturing engineering. Fracture properties of the model are calculated from core fracture data, according to statistical mathematical analysis. The calculation results make full use of the quantitative information of core fracture orientation, density, opening and length, which constitute the direct and extensive data of mining engineering. The reliability and applicability of the model are analyzed with regard to model size and density, a calculation method for dominant size and density being proposed. Then, finite element analysis is applied to a hydraulic fracturing numerical simulation of a shale fractured reservoir in southeastern Chongqing. The hydraulic pressure distribution, fracture propagation, acoustic emission information and in situ stress changes during fracturing are analyzed. The results show the application of fracture statistics in fracture modeling and the influence of fracture distribution on hydraulic fracturing engineering. The present analysis may provide a reference for shale gas exploitation.     

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr–Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils.     

Investigation of directional hydraulic fracturing based on true tri-axial experiment and finite element modeling

The initiation and propagation of directional hydraulic fracturing (DHF) was investigated based on true tri-axial experiment and finite element modeling. The influences of notch angle, notch length and injection rate on the DHF were investigated. The initiation and propagation of DHF was modeled by a 3D nonlinear finite element method. A comparison between experimental investigation and numerical modeling results indicates that there is a good correlation between unbalanced force (UF) and fracturing. UF can be used to predict the hydraulic fracture initiation and propagation.     

页岩水平井多段分簇压裂起裂压力数值模拟

页岩储层孔隙度和渗透率极低,天然裂缝和水平层理发育,常规压裂增产措施无法满足页岩气的开发要求,水平井多段分簇压裂是页岩气开发的关键技术之一,该技术能够大幅度提升压裂改造的体积、产气量和最终采收率。为确定页岩储层水平井多段分簇射孔压裂的起裂点和起裂压力,采用有限元方法建立了水平井套管完井（考虑水泥环和套管的存在）多段分簇射孔的全三维起裂模型。数值模型的起裂压力与室内试验结果吻合较好,证明了数值模型的准确性和可靠性。利用数值模型研究了页岩水平井多段分簇射孔压裂的起裂点和起裂压力的影响因素,研究发现：射孔孔眼附近无天然裂缝或水平层理影响,起裂点发生在射孔簇孔眼的根部;射孔簇间距越小,中间射孔簇的干扰越大,可能造成中间的射孔簇无法起裂;射孔密度和孔眼长度增大,起裂压力降低;天然裂缝的存在,在某些情况能够降低起裂压力且改变起裂位置,主要与天然裂缝的分布方位及水平主应力差有关;水平层理可能会降低起裂压力,但与垂向主应力与水平最小主应力的差值有关。获得的起裂压力变化规律,可作为进一步研究水平井多段分簇射孔条件下的裂缝扩展规律的基础,可以为压裂设计和施工的射孔参数确定及优化给出具体建议。