围压下泥岩断裂韧性测试与解释方法

根据围压条件下的断裂力学理论 ,采用双翼裂纹厚壁圆筒 ,对泥岩进行了不同围压和不同泥质含量条件下的人工岩样断裂韧性测试。建立了有限元方法解释断裂韧性的数值模型。在进行大量实验测试的基础上 ,通过回归分析表明 ,断裂韧性与围压、泥质含量有较好的线性统计关系。解决了水力压裂设计中断裂韧性参数的准确获取问题 ,有利于提高压裂设计、预测的精度     

考虑断裂韧性影响的页岩气储层可压性评价方法

可压性是衡量页岩气储层对压裂改造响应程度的重要依据,目前针对页岩可压性常用脆性指数对其进行评价。但在实际运用中存在脆性指数相近,实际压裂效果却差别甚远的问题。为了弥补利用脆性指数进行页岩气储层可压性评价的不足,在考虑断裂韧性影响的基础上引入临界机械能释放率,综合了Ⅰ型、Ⅱ型断裂韧性对岩石裂缝发育的影响,并定义平均临界机械能释放率。最后结合平均临界机械能释放率和脆性指数,初步建立了以裂缝发育指数为量化指标的页岩气储层可压性评价模型。基于焦石坝龙马溪组焦页48-2HF井测井资料,计算出其弹性模量、泊松比、Ⅰ型及Ⅱ型断裂韧性,并根据构建的裂缝发育指数模型,得出页岩储层可压性随埋深的变化关系。结合现场微震资料,验证了该模型能有效预测页岩气储层可压性随埋深的变化趋势,为页岩气实际开发过程中优选页岩储层压裂段提供可靠方法。      

软煤储层顶板水平井穿层工况下压裂缝扩展模型

为了探索煤层顶板中水平井向目标层穿层压裂的裂缝扩展规律,以华北石炭-二叠纪煤田为例,运用断裂力学、损伤力学以及流体力学等经典理论并结合现场实测资料,开展了压裂缝延伸距离与压裂时间时空演化规律的建模与验证。首先,基于原生裂缝特性、渗透特性以及压裂射孔段附加应力等因素,提出了顶板水平井垂向造缝的起裂压力计算公式;其次,在考虑裂缝性煤岩体损伤效应的基础上,引入Dougill损伤因子,将该计算模型拓展为延伸压力计算模型;最后,基于改进的经典PKN裂缝模型和压裂液滤失理论,建立了连续穿层工况下压裂缝延伸距离与压裂施工时间的函数关系。实践验证表明,根据理论模型合理调配时间参数,可以控制穿层裂缝的延伸距离。      

页岩I型断裂韧性的分形计算方法与应用

水平井分段体积压裂是页岩气商业规模开发的重要工艺措施,如何评价页岩储层可压性是该工艺成功的关键。页岩断裂韧性是可压性评价的重要支撑参数,从I型断裂裂纹（裂缝）微观形态入手,结合断裂力学理论和分形理论,建立了页岩I型断裂韧性分形计算方法;借助晶体劈裂功法计算页岩表面能,采用密度、声波时差2种参数获取计算结果,对比分析新的分形方法计算数据、传统方法预测数据与试验测试数据,新的分形方法计算平均误差为3.63%,传统预测方法平均误差为 %,验证了方法的准确性;参考实例水平井测井解释数据,计算了水平段I型断裂韧性指数的全井筒连续性剖面,结合可压性级别与较大储层改造体积概率关系,优选可压性级别为III级及II级改造。建立的页岩I型断裂韧性分形计算方法对定量评价页岩储层可压性具有重要意义。     

深层致密砂岩储层可压裂性评价新方法

岩石可压性评价是储层压裂改造层位优选、压后产能评估的重要基础工作。准中4区块致密砂岩储层埋藏深、物性差,亟需通过压裂改造提高工业产能。因此,以董2井北三维区侏罗系致密砂岩为例,基于岩石三轴实验建立了致密砂岩断裂能密度—弹性模量的拟合公式,采用矿物成分法和弹模-泊松比法确定了研究区不同深度岩石脆性指数,采用岩石破裂准则确定了研究区不同深度的裂缝发育指数。以断裂能密度表征致密砂岩断裂韧性,以裂缝发育指数表征储层天然裂缝发育程度,综合考虑岩石脆性、断裂韧性、地应力环境和天然裂缝发育程度的影响,采用层次分析法计算了各因素权重,建立了适合深层致密砂岩的可压性评价方法。研究结果表明,可压裂性指数大于0.55时,可压性好;可压裂性指数介于0.50~0.55之间时,可压性一般;可压裂性指数小于0.50时,可压性差;研究区D7井的最佳压裂层位为4145~4160 m、4470~4480 m、5290~5330 m,D8井的最佳压裂层位为5120~5330 m、5350~5365 m,D701井的最佳压裂层位为3900~3910 m、4430~4440 m、4455~4465 m、5125~5135 m。      

Geomechanical Characteristics of Gas Shales: A Case Study in the North Perth Basin

Gas shales are one type of unconventional reservoirs which have attracted significant attention for gas production in recent years. Gas production from very tight shales requires employment of hydraulic fracturing as a stimulation technique. To design hydraulic fracture operation the mechanical properties of the targeted and surrounding formations should be estimated. Also, the magnitude and orientation of in situ stresses in the field need to be known to estimate the fracture initiation and propagation pressures. This study focuses on gas shale characteristics in the North Perth Basin and uses data corresponding to well Arrowsmith-2 (AS-2) which is the first dedicated shale gas well drilled in Western Australia. A log-based analysis was used to build the rock mechanical model (RMM). The RMM results were used to set up a hydraulic fracturing laboratory experiment. The test was done in the presence of three principal stresses to mimic the real field stress conditions. The test results include the pressure–time curve which was used to estimate the initiation and propagation pressure at that depth. The results were used to draw some practical conclusions related to hydraulic fracturing operation in the field.     

地应力分层技术在压裂设计优化中的应用

纵向连续地应力分布对井中储层的压裂施工设计具有重要的指导意义。以卫星油田葡萄花油藏储层为研究对象,结合实际储层压裂数据,优化建立由连续测井数据计算地层最小主应力的数学模型,给出合理的沿井轴连续的地应力剖面。以自然电位曲线的分层为基准,按照一定规则将应力剖面曲线、自然伽马曲线分层取值,并计算待压裂层与上(下)隔层间的"隔层厚度"、"储隔层应力差"、"隔层品质"等参数。利用研究区实际压裂设计数据,以避免压裂过程中待压裂岩层发生"窜槽"为目的,研究给出了不同参数岩层,在压裂排量、压砂强度、单位厚度用液量及压裂方式等施工参数方面的约束条件。实际井例的应用说明,这些条件可有效保障勘探开发压裂工艺的成功率。     

Atomistic model of diopside-K-jadeite (CaMgSi<Subscript>2</Subscript>O<Subscript>6</Subscript>-KAlSi<Subscript>2</Subscript>O<Subscript>6</Subscript>) solid solution

Atomistic model was proposed to describe the thermodynamics of mixing in the diopside-K-jadeite solid solution (CaMgSi₂O6-KAlSi₂O₆). The simulations were based on minimization of the lattice energies of 800 structures within a 2 × 2 × 4 supercell of C2/c diopside with the compositions between CaMgSi₂O₆ and KAlSi₂O₆ and with variable degrees of order/disorder in the arrangement of Ca/K cations in M2 site and Mg/Al in Ml site. The energy minimization was performed with the help of a force-field model. The results of the calculations were used to define a generalized Ising model, which included 37 pair interaction parameters. Isotherms of the enthalpy of mixing within the range of 273–2023 K were calculated with a Monte Carlo algorithm, while the Gibbs free energies of mixing were obtained by thermodynamic integration of the enthalpies of mixing. The calculated T-X diagram for the system CaMgSi₂O₆-KAlSi₂O₆ at temperatures below 1000 K shows several miscibility gaps, which are separated by intervals of stability of intermediate ordered compounds. At temperatures above 1000 K a homogeneous solid solution is formed. The standard thermodynamic properties of K-jadeite (KAlSi₂O₆) evaluated from quantum mechanical calculations were used to determine location of several mineral reactions with the participation of the diopside-K-jadeite solid solution. The results of the simulations suggest that the low content of KalSi₂O₆ in natural clinopyroxenes is not related to crystal chemical factors preventing isomorphism, but is determined by relatively high standard enthalpy of this end member.     

On the risk of hydraulic fracturing in CO2 geological storage

We present a contribution on the risk of hydraulic fracturing in CO₂ geological storage using an analytical model of hydraulic fracturing in weak formations. The work is based on a Mohr–Coulomb dislocation model that is extended to account for material with fracture toughness. The complete slip process that is distributed around the crack tip is replaced by superdislocations that are placed in the effective centers. The analytical model enables the identification of a dominant parameter, which defines the regimes of brittle to ductile propagation and the limit at which a mode‐1 fracture cannot advance. We examine also how the corrosive effect of CO₂ on rock strength may affect hydraulic fracture propagation. We found that a hydraulically induced vertical fracture from CO₂ injection is more likely to propagate horizontally than vertically, remaining contained in the storage zone. The horizontal fracture propagation will have a positive effect on the injectivity and storage capacity of the formation. The containment in the vertical direction will mitigate the risk of fracturing and migration of CO₂ to upper layers and back to the atmosphere. Although the corrosive effect of CO₂ is expected to decrease the rock toughness and the resistance to fracturing, the overall decrease of rock strength promotes ductile behavior with the energy dissipated in plastic deformation and hence mitigates the mode‐1 fracture propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

原状冻土非线性断裂韧度测试与修正

在以往原状冻土现场试验方法的基础上,针对原状冻土的非线性断裂力学特性进行了Ⅰ型和Ⅱ型断裂力学试验。在不同深度的冻胀量、冻深与时间变化规律下严格控制试样制作及试验时的温度,利用动态数据采集系统测出试样加力点处位移与力的关系曲线,并计算出相应的非线性参数。提出了利用非线性修正因子计算原状冻土断裂韧度的方法,并算出Ⅰ型、Ⅱ型直裂纹加卸载试样的非线性断裂韧度,将其与岩石断裂力学中的塑性修正因子方法的计算结果进行了对比,二者基本一致,说明该方法可以用来计算、修正原状冻土非线性断裂韧度,从而为冻土非线性理论的研究拓宽了思路。     

Parameters controlling pressure and fracture behaviors in field injectivity tests: A numerical investigation using coupled flow and geomechanics model

Field injectivity tests are widely used in the oil and gas industry to obtain key formation characteristics. The prevailing approaches for injectivity test interpretation rely on traditional analytical models. A number of parameters may affect the test results and lead to interpretation difficulties. Understanding their impacts on pressure response and fracture geometry of the test is essential for accurate test interpretation. In this work, a coupled flow and geomechanics model is developed for numerical simulation of field injectivity tests. The coupled model combines a cohesive zone model for simulating fluid-driven fracture and a poro-elastic/plastic model for simulating formation behavior. The model can capture fracture propagation, fluid flow within the fracture and formation, deformation of the formation, and evolution of pore pressure and stress around the wellbore and fracture during the tests. Numerical simulations are carried out to investigate the impacts of a multitude of parameters on test behaviors. The parameters include rock permeability, the leak-off coefficient of the fracture, rock stiffness, rock toughness, rock strength, plasticity deformation, and injection rate. The sensitivity of pressure response and fracture geometry on each parameter is reported and discussed. The coupled flow and geomechanics model provides additional advantages in the understanding of the fundamental mechanisms of field injectivity tests.     

Pore pressure during metamorphism of carbonate rock: effect of volumetric properties of H<Subscript>2</Subscript>O–CO<Subscript>2</Subscript> mixtures

The molar volume of mixtures of CO₂ and H₂O is a strong function of the fluid composition. Both CO₂ and H₂O participate in the metamorphism of carbonate rocks, resulting in a change in the fluid composition during reaction. One of the effects of the change in composition is the increase in pore-fluid pressure with only small increases in extent of reaction, ;. Pressure calculated from the volumetric properties of CO₂-H₂O mixtures at 400 °C increases greatly with small increases of ; but drops at greater values because of the increase in pore volume as a result of (V_solid. The pore pressure increase at small values of ;, though, readily exceeds the reported tensile strength of carbonate rocks, and the rock cannot sustain significant reaction without fracturing. The result of a small amount of reaction is a fractured rock with increased permeability, which promotes fluid transport.     

运用水压破裂法测试岩石的抗拉强度

劈裂法受试验条件的影响很难得出理想的结果,而水压破裂试验有试件受力载荷不受人为因素干扰的优点,测得的抗拉强度具有较高的可信度,并且在工程实际中有广泛的用途。本文在研究水压破裂法测试岩石抗拉强度的方法及原理的基础上,用统计断裂力学整理试验结果,并与巴西法作了比较。试验表明,用威布尔二参数模型的统计断裂力学方法研究岩石抗拉强度是可行的,并具有良好的线性关系。     

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

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

水平井中多条裂缝同步扩展时裂缝竞争机制

水平井中多条水力裂缝间的应力干扰行为,造成了压裂液排量的非均匀分配,影响了水力裂缝的几何形态。采用边界元法研究岩体在压裂液作用下的变形程度,以幂律流体泊肃叶平板流动方程来研究水力裂缝内部的压裂液流场,考虑了多条裂缝间应力干扰和压裂液流量分配,建立了流-固耦合的水平井多条水力裂缝同步扩展模型。模型可模拟水平井多条水力裂缝几何形态、应力干扰情况和压裂液排量的分配情况,可解释水力裂缝之间的竞争机制。多条裂缝同步扩展时,压裂液排量并非均等地分配到各个裂缝之中,进入到内部裂缝的压裂液流量最小,内部裂缝宽度最小;内部的水力裂缝增长一定长度后,停止增长,并且在应力干扰下逐渐闭合。     

Temperature-driven meltwater production and hydrochemical variations at a glaciated alpine karst aquifer: implication for the atmospheric CO<Subscript>2</Subscript> sink under global warming

About two hydrological years of continuous data of discharge, temperature, electrical conductivity and pH have been recorded at the Glarey spring in the Tsanfleuron glaciated karst area in the Swiss Alps, to understand how glaciated karst aquifer systems respond hydrochemically to diurnal and seasonal recharge variations, and how calcite dissolution by glacial meltwater contributes to the atmospheric CO₂ sink. A thermodynamic model was used to link the continuous data to monthly water quality data allowing the calculation of CO₂ partial pressures and calcite saturation indexes. The results show diurnal and seasonal hydrochemical variations controlled chiefly by air temperature, the latter influencing karst aquifer recharge by ice and snowmelt. Karst process-related atmospheric CO₂ sinks were more than four times higher in the melting season than those in the freezing season. This finding has implication for understanding the atmospheric CO₂ sink in glaciated carbonate rock terrains: the carbon sink will increase with increasing runoff caused by global warming, i.e., carbonate weathering provides a negative feedback for anthropogenic CO₂ release. However, this is a transient regulation effect that is most efficient when glacial meltwater production is highest, which in turn depends on the future climatic evolution.     

An explicitly coupled hydro‐geomechanical model for simulating hydraulic fracturing in arbitrary discrete fracture networks

Modeling hydraulic fracturing in the presence of a natural fracture network is a challenging task, owing to the complex interactions between fluid, rock matrix, and rock interfaces, as well as the interactions between propagating fractures and existing natural interfaces. Understanding these complex interactions through numerical modeling is critical to the design of optimum stimulation strategies. In this paper, we present an explicitly integrated, fully coupled discrete‐finite element approach for the simulation of hydraulic fracturing in arbitrary fracture networks. The individual physical processes involved in hydraulic fracturing are identified and addressed as separate modules: a finite element approach for geomechanics in the rock matrix, a finite volume approach for resolving hydrodynamics, a geomechanical joint model for interfacial resolution, and an adaptive remeshing module. The model is verified against the Khristianovich–Geertsma–DeKlerk closed‐form solution for the propagation of a single hydraulic fracture and validated against laboratory testing results on the interaction between a propagating hydraulic fracture and an existing fracture. Preliminary results of simulating hydraulic fracturing in a natural fracture system consisting of multiple fractures are also presented. Copyright © 2012 John Wiley & Sons, Ltd.     

测试压裂分析技术在渤海B油田的应用

测试压裂分析技术是水力压裂施工前获取储层物性参数和分析压裂液性质的手段,是压裂施工顺利进行的技术保证。结合渤海某油田沙河街油藏压裂工程,介绍了测试压裂应用过程和理论基础,并简单介绍几种求取地层闭合压力的数学方法。闭合压力是压裂施工关键参数之一,不同数学计算方法推导出的闭合压力不尽相同,现场压裂工程师要根据实际测试结果选取可靠的数值。压降曲线拟合是在模型中反推出地层关键参数的方法,拟合完毕即得到针对施工地层的压裂模型,在此基础上进行主压裂设计和裂缝预测是较为准确的,以此进行施工也是安全的。     

Transient solution for a plane‐strain fracture driven by a shear‐thinning,power‐law fluid

This paper analyses the problem of a fluid‐driven fracture propagating in an impermeable, linear elastic rock with finite toughness. The fracture is driven by injection of an incompressible viscous fluid with power‐law rheology. The relation between the fracture opening and the internal fluid pressure and the fracture propagation in mobile equilibrium are described by equations of linear elastic fracture mechanics (LEFM), and the flow of fluid inside the fracture is governed by the lubrication theory. It is shown that for shear‐thinning fracturing fluids, the fracture propagation regime evolves in time from the toughness‐ to the viscosity‐dominated regime. In the former, dissipation in the viscous fluid flow is negligible compared to the dissipation in extending the fracture in the rock, and in the later, the opposite holds. Corresponding self‐similar asymptotic solutions are given by the zero‐viscosity and zero‐toughness (J. Numer. Anal. Meth. Geomech. 2002; 26 :579–604) solutions, respectively. A transient solution in terms of the crack length, the fracture opening, and the net fluid pressure, which describes the fracture evolution from the early‐time (toughness‐dominated) to the large‐time (viscosity‐dominated) asymptote is presented and some of the implications for the practical range of parameters are discussed. Copyright © 2006 John Wiley & Sons, Ltd.     

射孔水平井分段压裂破裂点优化方法

为确定水平井裂缝破裂点位置,指导射孔位置的优选,建立了综合考虑多条人工裂缝诱导应力、井筒内压、原地应力、热应力、压裂液渗流效应和射孔联合作用下的水平井分段压裂应力场模型,模型重点考虑了先压开裂缝产生的诱导应力场影响,更符合水平井分段压裂的实际特点,在此基础上结合岩石破裂准则从降低破裂压力角度形成了破裂点优化方法。研究结果表明先压开裂缝产生的诱导应力场影响使得水平井筒周围应力场更为复杂,研究裂缝起裂问题时必须考虑这一因素。对川西水平井分段压裂的破裂点进行了优化,压裂后输气求产,实际破裂压力梯度低于邻井同层的破裂压力梯度,大大降低了施工风险。该破裂点优化方法对水平井分段压裂优化设计有一定的指导意义。