共查询到20条相似文献，搜索用时 15 毫秒
1.
Numerical modeling of hydraulic fracture propagation,closure and reopening using XFEM with application to in‐situ stress estimation 
下载免费PDF全文
下载免费PDF全文 In this paper, a fully coupled model is developed for numerical modeling of hydraulic fracturing in partially saturated weak porous formations using the extended finite element method, which provides an effective means to simulate the coupled hydro‐mechanical processes occurring during hydraulic fracturing. The developed model is for short fractures where plane strain assumptions are valid. The propagation of the hydraulic fracture is governed by the cohesive crack model, which accounts for crack closure and reopening. The developed model allows for fluid flow within the open part of the crack and crack face contact resulting from fracture closure. To prevent the unphysical crack face interpenetration during the closing mode, the crack face contact or self‐contact condition is enforced using the penalty method. Along the open part of the crack, the leakage flux through the crack faces is obtained directly as a part of the solution without introducing any simplifying assumption. If the crack undergoes the closing mode, zero leakage flux condition is imposed along the contact zone. An application of the developed model is shown in numerical modeling of pump‐in/shut‐in test. It is illustrated that the developed model is able to capture the salient features bottomhole pressure/time records exhibit and can extract the confining stress perpendicular to the direction of the hydraulic fracture propagation from the fracture closure pressure. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
2.
Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock 下载免费PDF全文
下载免费PDF全文 The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution. 相似文献
3.
The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
4.
Coupled hydromechanical‐fracture simulations of nonplanar three‐dimensional hydraulic fracture propagation 下载免费PDF全文
下载免费PDF全文 This paper presents an algorithm and a fully coupled hydromechanical‐fracture formulation for the simulation of three‐dimensional nonplanar hydraulic fracture propagation. The propagation algorithm automatically estimates the magnitude of time steps such that a regularized form of Irwin's criterion is satisfied along the predicted 3‐D fracture front at every fracture propagation step. A generalized finite element method is used for the discretization of elasticity equations governing the deformation of the rock, and a finite element method is adopted for the solution of the fluid flow equation on the basis of Poiseuille's cubic law. Adaptive mesh refinement is used for discretization error control, leading to significantly fewer degrees of freedom than available nonadaptive methods. An efficient computational scheme to handle nonlinear time‐dependent problems with adaptive mesh refinement is presented. Explicit fracture surface representations are used to avoid mapping of 3‐D solutions between generalized finite element method meshes. Examples demonstrating the accuracy, robustness, and computational efficiency of the proposed formulation, regularized Irwin's criterion, and propagation algorithm are presented. 相似文献
5.
In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
6.
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. 相似文献
7.
Coupled formulation and algorithms for the simulation of non‐planar three‐dimensional hydraulic fractures using the generalized finite element method 下载免费PDF全文
下载免费PDF全文 This paper presents a coupled hydro‐mechanical formulation for the simulation of non‐planar three‐dimensional hydraulic fractures. Deformation in the rock is modeled using linear elasticity, and the lubrication theory is adopted for the fluid flow in the fracture. The governing equations of the fluid flow and elasticity and the subsequent discretization are fully coupled. A Generalized/eXtended Finite Element Method (G/XFEM) is adopted for the discretization of the coupled system of equations. A Newton–Raphson method is used to solve the resulting system of nonlinear equations. A discretization strategy for the fluid flow problem on non‐planar three‐dimensional surfaces and a computationally efficient strategy for handling time integration combined with mesh adaptivity are also presented. Several three‐dimensional numerical verification examples are solved. The examples illustrate the generality and accuracy of the proposed coupled formulation and discretization strategies. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
8.
水力压裂技术广泛使用于页岩气开采工程中。为了分析压裂过程中多裂缝扩展形成复杂裂缝网的机制,尝试将态型近场动力学理论引入页岩水平井水力压裂过程的力学建模与数值仿真,在物质点间相互作用力模型中加入等效水压力项以实现在新生裂缝面上跟踪施加水压力,建立了水力压裂过程的近场动力学分析模型。通过模拟页岩储层的水力压裂过程,可得到复杂的裂缝扩展路径、裂缝网络的形成过程以及裂缝扩展受射孔间距及页岩天然裂缝和层理的影响。研究结果表明:射孔间距过小会造成起裂干扰,使中心射孔的裂缝扩展受到抑制;在压裂压力一定的情况下适当增大射孔间距,可以显著增强页岩压裂形成裂缝网的能力;压裂过程中水平层理面可能张开形成水平裂缝,且天然裂缝会诱导形成更复杂的垂直裂缝。模型和方法可为页岩水力压裂过程和机制研究及工程实践提供参考。 相似文献
9.
水力压裂作为煤层强化增透技术的一种,其应力演化特征及裂隙形态与扩展范围的判断尤为重要。采用离散元数值方法,以导向压裂为背景,建立水力压裂流固耦合模型;通过应力路径、裂纹热点图等手段,探究水力压裂过程中压裂排量、泊松比、天然裂隙密度对应力演化和裂隙演化的影响及其细观规律。结果表明:不同压裂排量下的应力演化方向及最终应力路径曲线形状有着明显的不同,低排量下裂隙附近的应力比值逐渐增大,而在高排量下先增大后减小;煤层泊松比越大,平均压裂半径越低,但对起裂时间及裂隙的扩展形态影响不明显;天然裂隙的发育情况对水力裂隙的扩展起着关键性作用,高裂隙发育煤层水力裂隙扩展的方向性无法预测,应力演化方向会出现反转现象;压裂过程中不同区域的应力演化特征能够反映出裂隙的扩展状态,现场可通过监测压裂区域附近应力变化,判断水力压裂缝网的扩展范围。 相似文献
10.
Propagation of fractures, especially those emanating from wellbores and closed natural fractures, often involves Mode I and Mode II, and at times Mode III, posing significant challenges to its numerical simulation. When an embedded inclined fracture is subjected to compression, the fracture edge is constrained by the surrounding materials so that its true propagation pattern cannot be simulated by 2D models. In this article, a virtual multidimensional internal bond (VMIB) model is presented to simulate three‐dimensional (3D) fracture propagation. The VMIB model bridges the processes of macro fracture and micro bond rupture. The macro 3D constitutive relation in VMIB is derived from the 1D bond in the micro scale and is implemented in a 3D finite element method. To represent the contact and friction between fracture surfaces, a 3D element partition method is employed. The model is applied to simulate fracture propagation and coalescence in typical laboratory experiments and is used to analyze the propagation of an embedded fracture. Simulation results for single and multiple fractures illustrate 3D features of the tensile and compressive fracture propagation, especially the propagation of a Mode III fracture. The results match well with the experimental observation, suggesting that the presented method can capture the main features of 3D fracture propagation and coalescence. Moreover, by developing an algorithm for applying pressure on the fracture surfaces, propagation of a natural fracture is also simulated. The result illustrates an interesting and important phenomenon of Mode III fracture propagation, namely the fracture front segmentation. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
11.
A fully coupled element‐free Galerkin model for hydro‐mechanical analysis of advancement of fluid‐driven fractures in porous media 下载免费PDF全文
下载免费PDF全文 Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi‐analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element‐free Galerkin (EFG) mesh‐less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity equations, considers the hydro‐mechanical coupling between the crack and its surrounding porous medium. Therefore, the developed EFG model is capable of simulating fluid leak‐off and fluid lag phenomena. To create the discrete equation system, the Galerkin technique is applied, and the essential boundary conditions are imposed via penalty method. Then, the resultant constrained integral equations are discretized in space using EFG shape functions. For temporal discretization, a fully implicit scheme is employed. The final set of algebraic equations that forms a non‐linear equation system is solved using the direct iterative procedure. Modeling of cracks is performed on the basis of linear elastic fracture mechanics, and for this purpose, the so‐called diffraction method is employed. For verification of the model, a number of problems are solved. According to the obtained results, the developed EFG computer program can successfully be applied for simulating the complex process of hydraulic fracture propagation in porous media. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
12.
储层岩体中的天然结构面对水力压裂缝网改造具有重要的影响。本文采用真实破裂过程分析软件RFPA2D-Flow,在考虑岩体非均质性和岩体渗流-应力-损伤破裂特性的基础上,对不同尺度天然结构面影响的水力压裂裂缝扩展与演化行为进行了模拟分析和讨论,研究结果表明:（1）当水力裂缝遇天然非闭合裂隙时,在水力裂缝靠近非闭合裂隙区间形成拉张应力区,水力裂缝与区间非闭合裂隙间微元体累进性张拉破坏是导致水力裂缝与非闭合裂隙贯通的主要机制;（2）层理等优势结构对水力压裂裂缝扩展及缝网形态影响十分显著,当最大主应力方向与层理面走向小角度相交时,层理结构面对水力裂隙的扩展起主要作用,当最大主应力方向与层理面走向大角度相交时,最大主压应力与层理面共同对缝网扩展起主导作用,随着优势结构面的增多和差应力的增大,水力压裂形成的缝网范围和复杂性程度随之增大;（3）储层水力压裂是一种局部范围内的短暂动力扰动过程,尽管断层的存在可以极大地影响水力裂缝的扩展模式,增大水力裂隙扩展高度,但相比于储层埋深,水力压裂对断层封闭性的破坏范围和断层活动性的扰动程度十分有限。 相似文献
13.
水力压裂是青海共和盆地干热岩地热资源开发的难点技术问题之一。本文基于升级改造的大尺寸真三轴水力压裂物理模拟实验系统模拟干热岩储层高温高压环境,利用青海共和盆地露头岩心进行水力压裂物理模拟实验,揭示干热岩储层水力裂缝的起裂和扩展规律。通过物理模拟实验发现：干热岩储层裂缝起裂可以通过文中提出的起裂模型判断起裂方式和预测起裂压力;水力裂缝在岩石基质中的扩展形态简单,仅沿最大主应力方向延伸;但是水力裂缝会受到岩石中弱面的影响,发生转向沿弱面延伸,形成较复杂的裂缝形态。因此,建议在干热岩储层实际施工中,在天然裂缝发育较丰富的层段开展水力压裂,以实现复杂裂缝网络提取地热能。 相似文献
14.
建立了求解自然裂纹和水力裂纹扩展的扩展有限元法,对裂纹附近区域的节点采用广义形函数,并采用线增函数消除混合单元,以提高裂纹附近的精度。引入水力劈裂的非耦合模型,即假设裂纹中的水压力为均布力;用砂浆法（线段-线段接触法）结合增广型拉格朗日乘子法处理受压裂纹段的接触条件。并通过算例分析了以下内容：计算了受压裂纹和裂纹面分布均布水压力的水力裂纹的应力强度因子,并与解析解进行了比较,结果表明,提出的方法具有很高的精度;模拟了水力裂纹对自然裂纹面的影响,并分析了自然裂纹面上的接触力和接触状态。 相似文献
15.
Yao Yao 《Rock Mechanics and Rock Engineering》2012,45(3):375-387
Hydraulic fracturing technology is being widely used within the oil and gas industry for both waste injection and unconventional
gas production wells. It is essential to predict the behavior of hydraulic fractures accurately based on understanding the
fundamental mechanism(s). The prevailing approach for hydraulic fracture modeling continues to rely on computational methods
based on Linear Elastic Fracture Mechanics (LEFM). Generally, these methods give reasonable predictions for hard rock hydraulic
fracture processes, but still have inherent limitations, especially when fluid injection is performed in soft rock/sand or
other non-conventional formations. These methods typically give very conservative predictions on fracture geometry and inaccurate
estimation of required fracture pressure. One of the reasons the LEFM-based methods fail to give accurate predictions for
these materials is that the fracture process zone ahead of the crack tip and softening effect should not be neglected in ductile
rock fracture analysis. A 3D pore pressure cohesive zone model has been developed and applied to predict hydraulic fracturing
under fluid injection. The cohesive zone method is a numerical tool developed to model crack initiation and growth in quasi-brittle
materials considering the material softening effect. The pore pressure cohesive zone model has been applied to investigate
the hydraulic fracture with different rock properties. The hydraulic fracture predictions of a three-layer water injection
case have been compared using the pore pressure cohesive zone model with revised parameters, LEFM-based pseudo 3D model, a
Perkins-Kern–Nordgren (PKN) model, and an analytical solution. Based on the size of the fracture process zone and its effect
on crack extension in ductile rock, the fundamental mechanical difference of LEFM and cohesive fracture mechanics-based methods
is discussed. An effective fracture toughness method has been proposed to consider the fracture process zone effect on the
ductile rock fracture. 相似文献
16.
Fully coupled simulation of a hydraulic fracture interacting with natural fractures with a hybrid discrete‐continuum method 下载免费PDF全文
下载免费PDF全文 A hybrid discrete‐continuum numerical scheme is developed to study the behavior of a hydraulic fracture crossing natural fractures. The fully coupled hybrid scheme utilizes a discrete element model for an inner domain, within which the hydraulic fracture propagates and interacts with natural fractures. The inner domain is embedded in an outer continuum domain that is implemented to extend the length of the hydraulic fracture and to better approximate the boundary effects. The fracture is identified to propagate initially in the viscosity‐dominated regime, and the numerical scheme is calibrated by using the theoretical plane strain hydraulic fracture solution. The simulation results for orthogonal crossing indicate three fundamental crossing scenarios, which occur for various stress ratios and friction coefficients of the natural fracture: (i) no crossing, that is, the hydraulic fracture is arrested by the natural fracture and makes a T‐shape intersection; (ii) offset crossing, that is, the hydraulic fracture crosses the natural fracture with an offset; and (iii) direct crossing, that is, the hydraulic fracture directly crosses the natural fracture without diversion. Each crossing scenario is associated with a distinct net pressure history. Additionally, the effects of strength contrast and stiffness contrast of rock materials and intersection angle between the hydraulic fracture and the natural fracture are also investigated. The simulations also illustrate that the level of fracturing complexity increases as the number and extent of the natural fractures increase. As a result, we can conclude that complex hydraulic fracture propagation patterns occur because of complicated crossing behavior during the stimulation of naturally fractured reservoirs. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
17.
随着扩展有限元理论的深入研究,利用扩展有限元方法模拟水力压裂具有了一定的可操作性。相比于常规有限元方法,XFEM方法具有计算结果精度高和计算量小的优点。但是,如何模拟射孔孔眼、如何模拟流体与岩石相互作用以及分析水力裂缝的扩展规律仍然是难题。以研究水力压裂裂缝扩展规律为目的,建立了岩石多孔介质应力平衡方程、流体渗流连续性方程和边界条件。通过有限元离散化方法对耦合方程矩阵进行处理。通过富集函数定义初始裂缝（射孔孔眼）,选择最大主应力及损伤变量D分别作为裂缝起裂和扩展判定准则,利用水平集方法模拟水力裂缝扩展过程。数值模拟结果显示：增加射孔方位角、压裂液排量和减小水平地应力差,起裂压力上升;黏度对起裂压力无明显影响。增加射孔方位角、压裂液排量、黏度和减小水平地应力差值有助于裂缝宽度的增加。增加水平地应力差值、压裂液排量和减小射孔方位角以及压裂液黏度有助于裂缝长度增加,反之亦然。基于ABAQUS的水力裂缝扩展有限元法可对不同井型和诸多储层物性参数及压裂施工参数进行分析,且裂缝形态逼真,裂缝面凹凸程度清晰,结果准确。此研究可作为一种简便有效研究水力压裂裂缝扩展规律的方法为油田水力压裂设计与施工提供参考与依据。 相似文献
18.
Analysis of refracturing in horizontal wells: Insights from the poroelastic displacement discontinuity method 下载免费PDF全文
下载免费PDF全文 Ali Rezaei Giorgio Bornia Mehdi Rafiee Mohamed Soliman Stephen Morse 《国际地质力学数值与分析法杂志》2018,42(11):1306-1327
In this paper, a fully coupled 2‐dimensional poroelastic displacement discontinuity method is used to investigate the refracturing process in horizontal wells. One of the objectives of refracturing is to access new reserves by adding new hydraulic fractures in zones that were bypassed in the initial fracturing attempt. Pore pressure depletion in the vicinity of old fractures directly affects the state of stress and eventually the propagation of newly created hydraulic fractures. Thus, a poroelastic analysis is required to identify guidelines for the refracturing process, in particular to understand the extension of the pore pressure depletion, and eventually, the orientation of new as well as old fractures. We propose a fully coupled approach to model the whole process of child fracture propagation in a depleted area between 2 parent fractures in the same wellbore. This approach omits the need of using multistep workflow that is regularly used to model the process. The maximum tensile stress criterion (σ criterion) is used for hydraulic fracture propagation. The proposed method is verified using available analytical solutions for total stress and pore pressure loading modes on a line fracture in drained and undrained conditions. Then, test cases of multifractured horizontal wells are studied to calculate the time evolution of the stress and pore pressure fields around old fractures and to understand the effect of these fields on the propagation path of newly created fractures. Finally, the effect of the pore pressure depletion on the propagation path of the newly created fractures in the bypassed area of the same wellbore is studied. The results show that the depleted areas around old fractures are highly affected by the extent and severity of the stress redistribution and pore pressure depletion. It is observed that a successful creation of new fractures may only happen in certain time frames. The results of this study provide new insights on the behavior of newly created fractures in depleted zones. They also clarify the relationship between stress change and pore pressure depletion in horizontal wells. 相似文献
19.
Two-dimensional hydraulic fracturing simulations using the cohesive zone model (CZM) can be readily found in the literature; however, to our knowledge, verified 3D cohesive zone modeling is not available. We present the development of a 3D fully coupled hydro-mechanical finite element method (FEM) model (with parallel computation framework) and its application to hydraulic fracturing. A special zero-thickness interface element based on the CZM is developed for modeling fracture propagation and fluid flow. A local traction-separation law with strain softening is used to capture tensile cracking. The model is verified by considering penny-shaped hydraulic fracture and plain strain Kristianovich‑Geertsma‑de Klerk hydraulic fracture (in 3D) in the viscosity- and toughness-dominated regimes. Good agreement between numerical results and analytical solutions has been achieved. The model is used to investigate the influence of rock and fluid properties on hydraulic fracturing. Lower stiffness tip cohesive elements tend to yield a larger elastic deformation around the fracture tips before the tensile strength is reached, generating a larger fracture length and lower fracture pressure compared with higher stiffness elements. It is found that the energy release rate has almost no influence on hydraulic fracturing in the viscosity-dominated regime because the energy spent in creating new fractures is too small when compared with the total input energy. For the toughness-dominated regime, the released energy during fracturing should be accurately captured; relatively large tensile strength should be used in order to match numerical results to the asymptotic analytical solutions. It requires smaller elements when compared with those used in the viscosity-dominated regime. 相似文献
20.
水力压裂技术是成功实现干热岩资源开发利用的重要手段之一,数值模拟技术能够精准预测水力裂缝扩展。针对典型花岗岩,借助黏性单元法,分别模拟了致密花岗岩和天然裂缝存在情况下的水力裂缝扩展特征,得出以下结论：致密花岗岩的水力裂缝形态单一,天然裂缝的存在增加了压裂后裂缝的复杂性;致密花岗岩水力裂缝拓展主要分为憋压和拓展两个交替往复的阶段,当存在天然裂隙时,水力压裂过程会变得复杂;天然裂缝存在时,水力裂缝的缝长和缝宽分别为致密花岗岩的5. 7倍和1. 7倍;缝网的形成需要借助复杂的压裂工艺实现。研究结果可以为增强型地热系统(EGS)储层水力刺激工作提供理论支持。 相似文献