共查询到20条相似文献，搜索用时 500 毫秒
1.
2.
柴达木盆地西部古近系储层裂缝有效性的影响因素 总被引：2，自引：0，他引：2
古近系盐湖相裂缝性泥灰岩储层是柴达木盆地西部地区的重要储层类型,储层中发育剪切裂缝、纵张裂缝、顺层滑动裂缝和成岩裂缝4种类型。裂缝普遍被方解石、石膏或钙芒硝等充填,裂缝的有效性对储层评价十分重要。有效裂缝不仅提供了储层的重要储集空间和主要渗流通道,还决定了油气的富集规律和储层的潜在生产能力。裂缝的有效性主要受裂缝形成时间、溶蚀作用、第四纪以来的构造抬升剥蚀作用、异常流体高压作用、膏盐层分布和断层的活动性等因素控制。晚期形成的裂缝有效性更好;溶蚀作用可明显地增加裂缝的有效性;第四纪以来的构造抬升剥蚀作用和异常流体高压作用可使裂缝重新裂开,使裂缝的有效性变好;离膏盐层越远,裂缝的有效性越好。裂缝的有效性与现今地应力方向关系不密切,主要受断层活动的影响,与主断层平行的裂缝有效性最好。 相似文献
3.
Hydraulic fracturing is the method of choice to enhance reservoir permeability and well efficiency for extraction of shale gas. Multi‐stranded non‐planar hydraulic fractures are often observed in stimulation sites. Non‐planar fractures propagating from wellbores inclined from the direction of maximum horizontal stress have also been reported. The pressure required to propagate non‐planar fractures is in general higher than in the case of planar fractures. Current computational methods for the simulation of hydraulic fractures generally assume single, symmetric, and planar crack geometries. In order to better understand hydraulic fracturing in complex‐layered naturally fractured reservoirs, fully 3D models need to be developed. In this paper, we present simulations of 3D non‐planar fracture propagation using an adaptive generalized FEM. This method greatly facilitates the discretization of complex 3D fractures, as finite element faces are not required to fit the crack surfaces. A solution strategy for fully automatic propagation of arbitrary 3D cracks is presented. The fracture surface on which pressure is applied is also automatically updated at each step. An efficient technique to numerically integrate boundary conditions on crack surfaces is also proposed and implemented. Strongly graded localized refinement and analytical asymptotic expansions are used as enrichment functions in the neighborhood of fracture fronts to increase the computational accuracy and efficiency of the method. Stress intensity factors with pressure on crack faces are extracted using the contour integral method. Various non‐planar crack geometries are investigated to demonstrate the robustness and flexibility of the proposed simulation methodology. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
4.
N. B. Nagel M. A. Sanchez-Nagel F. Zhang X. Garcia B. Lee 《Rock Mechanics and Rock Engineering》2013,46(3):581-609
Due to the low permeability of many shale reservoirs, multi-stage hydraulic fracturing in horizontal wells is used to increase the productive, stimulated reservoir volume. However, each created hydraulic fracture alters the stress field around it, and subsequent fractures are affected by the stress field from previous fractures. The results of a numerical evaluation of the effect of stress field changes (stress shadowing), as a function of natural fracture and geomechanical properties, are presented, including a detailed evaluation of natural fracture shear failure (and, by analogy, the generated microseismicity) due to a created hydraulic fracture. The numerical simulations were performed using continuum and discrete element modeling approaches in both mechanical-only and fully coupled, hydro-mechanical modes. The results show the critical impacts that the stress field changes from a created hydraulic fracture have on the shear of the natural fracture system, which in-turn, significantly affects the success of the hydraulic fracture stimulation. Furthermore, the results provide important insight into: the role of completion design (stage spacing) and operational parameters (rate, viscosity, etc.) on the possibility of enhancing the stimulation of the natural fracture network (‘complexity’); the mechanisms that generate the microseismicity that occurs during a hydraulic fracture stimulation; and the interpretation of the generated microseismicity in relation to the volume of stimulated reservoir formation. 相似文献
5.
Liao Jianxing Gou Yang Feng Wentao Mehmood Faisal Xie Yachen Hou Zhengmeng 《Acta Geotechnica》2020,15(2):279-295
Although hydraulic fracturing has been massively studied and applied as a key technique to enhance the gas production from tight formations, some problems and uncertainties exist to accurately predict and analyze the fracture behavior in complex reservoirs, especially in the naturally fractured reservoirs like shale reservoirs. This paper presents a full 3D numerical model (FLAC3D) to study hydraulic fracturing behavior under the impact of preexisting orthogonal natural fractures. In this numerical model, the hydraulic fracture propagation direction is assumed perpendicular to the minimum principal stress and activated only by tensile failure, whereas the preexisting natural fractures can be activated by tensile or shear failure or a combination of them, and only tensile failure can open the natural fracture as well. The newly developed model was used to study the impact of preexisting orthogonal natural fractures on hydraulic fracturing behavior, based on a multistage hydraulic fracturing operation in a naturally fractured reservoir from the Barnett Shale formation, northwest of Texas in USA. In this multistage operation, two more representative stages, i.e., stage 1 with a relatively large horizontal stress anisotropy of 3.3 MPa and stage 4 with a comparatively small one of 1.3 MPa, were selected to conduct the simulation. Based on the numerical results, one can observe that the interaction between hydraulic and natural fracture is driven mainly by induced stress around fracture tip. Besides, the horizontal stress anisotropy plays a key role in opening the natural fracture. Thus, no significant opened fracture is activated on natural fracture in stage 1, while in stage 4 an opened fracture invades to about 90 m into the first natural fracture. Conversely, the hydraulic fracture length in stage 1 is much longer than in stage 4, as some fluid volume is stored in the opened natural fracture in stage 4. In this work, the shear failure on natural fractures is treated as the main factor for inducing the seismic events. And the simulated seismic events, i.e., shear failure on natural fractures, are very comparable with the measured seismic events.
相似文献6.
A virtual multidimensional internal bond (VMIB) model developed to simulate the propagation of hydraulic fractures using the finite‐element method is formulated within the framework of the virtual internal bond theory (VIB) that considers a solid as randomized material particles in the micro scale, and derives the macro constitutive relation from the cohesive law between the material particles with an implicit fracture criterion. Hydraulic pressure is applied using a new scheme that enables simulation of hydraulically driven cracks. When the model is applied to study hydraulic fracture propagation in the presence of a natural fracture, the results show the method to be very effective. It shows that although the in situ stress ratio is the dominant factor governing the propagation direction, a natural fault can also strongly influence the hydraulic fracture behavior. This influence is conditioned by the shear stiffness of the fault and the distance to the original hydraulic fracture. The model results show that when the fault is strong in shear, its impact on hydraulic fracture trajectory is weak and the hydraulic fracture will likely penetrate the fault. For a weak fault, however, the fracture tends to be arrested at the natural fault. The distance between the fault and the hydraulic fracture is also important; the fault influence increases with decreasing distance. The VMIB does not require selection of a fracture criterion and remeshing when the fracture propagates. Therefore, it is advantageous for modeling fracture initiation and propagation in naturally fractured rock. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
7.
人工压裂是获取干热岩型地热资源的关键环节,压裂后的人工裂隙带结构对开采条件下水热传递过程具有重要控制作用。结合我国共和盆地干热岩储层地质条件,采用数值模拟方法着重分析干热岩不同产状人工裂隙带的渗透率与宽度对热储中水热传递过程的影响机理,明确不同人工裂隙结构条件下水热产出能力,进而优化井间距。结果表明:当人工裂隙带渗透率较小时（小于5 D）,裂隙带规模越大,开采井温度越高;当渗透率较大时（大于10 D）,在水平裂隙带中,随着裂隙带规模的增加,由于注入冷水的快速扩散导致整体低温区域增加,开采井温度反而降低。在水平裂隙带中注入冷水主要为水平向流动,随着渗透率的增加,开采井温度更易受注入冷水的影响而降低;但在垂直裂隙带及倾斜裂隙带中,随着渗透率的增加,垂向自由对流增强,注入冷水更易于向储层底部高温区域流动,经加热后到达开采井,使得开采温度提升。综合比较,同一井间距条件下,低渗水平裂隙带以及高渗垂直裂隙带的产热能力较其他裂隙带更强。 相似文献
8.
Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method 
下载免费PDF全文
下载免费PDF全文 An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
9.
Thermal recovery from a hot dry rock (HDR) reservoir viewed as a deformable fractured medium is investigated with a focus on the assumption of local thermal non‐equilibrium (LTNE). Hydraulic diffusion, thermal diffusion, forced convection and deformation are considered in a two‐phase framework, the solid phase being made by impermeable solid blocks separated by saturated fractures. The finite element approximation of the constitutive and field equations is formulated and applied to obtain the response of a generic HDR reservoir to circulation tests. A change of time profile of the outlet fluid temperature is observed as the fracture spacing increases, switching from a single‐step pattern to a double‐step pattern, a feature which is viewed as characteristic of established LTNE. A dimensionless number is proposed to delineate between local thermal equilibrium (LTE) and non‐equilibrium. This number embodies local physical properties of the mixture, elements of the geometry of the reservoir and the production flow rate. All the above properties being fixed, the resulting fracture spacing threshold between LTNE and LTE is found to decrease with increasing porosity or fluid velocity. The thermally induced effective stress is tensile near the injection well, illustrating the thermal contraction of the rock, while the pressure contribution of the fracture fluid is negligible during the late period. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
10.
Yongquan Hu Zhiqiang Li Jinzhou Zhao Zhengwu Tao Pan Gao 《Environmental Earth Sciences》2017,76(15):546
The ultra-low-permeability shale gas reservoir has a lot of well-developed natural fractures. It has been proven that hydraulic fracture growth pattern is usually a complex network fracture rather than conventional single planar fractures by micro-seismic monitoring, which can be explained as the shear and tensile failure of natural fractures or creation of new cracks due to the increase in reservoir pore pressure caused by fluid injection during the process of hydraulic fracturing. In order to simulate the network fracture growth, a mathematical model was established based on full tensor permeability, continuum method and fluid mass conservation equation. Firstly, the governing equation of fluid diffusivity based on permeability tensor was solved to obtain the reservoir pressure distribution. Then Mohr–Coulomb shear failure criterion and tensile failure criterion were used to decide whether the rock failed or not in any block on the basis of the calculated reservoir pressure. The grid-block permeability was modified according to the change of fracture aperture once any type of rock failure criterion was met within a grid block. Finally, the stimulated reservoir volume (SRV) zone was represented by an enhancement permeability zone. After calibrating the numerical solution of the model with the field micro-seismic information, a sensitivity study was performed to analyze the effects of some factors including initial reservoir pressure, injection fluid volume, natural fracture azimuth angle and horizontal stress difference on the SRV (shape, size, bandwidth and length). The results show that the SRV size increases with the increasing initial pore reservoir and injection fluid volume, but decreases with the increase in the horizontal principal stress difference and natural fracture azimuth angle. The SRV shape is always similar for different initial pore reservoir and injection fluid volume. The SRV is observed to become shorter in length and wider in bandwidth with the decrease in natural fracture azimuth angle and horizontal principal stress difference. 相似文献
11.
Stress-Dependent Fluid Flow and Permeability in Fractured Media: from Lab Experiments to Engineering Applications 总被引：5，自引：3，他引：5
J. Zhang W. B. Standifird J.-C. Roegiers Y. Zhang 《Rock Mechanics and Rock Engineering》2007,40(1):3-21
Summary. Permeability is a physical property in rocks of extreme importance in energy engineering, civil and environmental engineering,
and various areas of geology. Early on, fractures in fluid flow models were assumed to be rigid. However, experimental research
and field data confirmed that stress-deformation behavior in fractures is a key factor governing their permeability tensor.
Although extensive research was conducted in the past, the three-dimensional stress-permeability relationships, particularly
in the inelastic deformation stage, still remain unclear. In this paper, laboratory experiments conducted on large concrete
blocks with randomly distributed fractures and rock core samples are reported to investigate fluid flow and permeability variations
under uniaxial, biaxial and triaxial complete stress-strain process. Experimental relationships among flowrate, permeability
and fracture aperture in the fractured media are investigated. Results show that the flowrate and stress/aperture exhibit
“cubic law” relationship for the randomly distributed fractures.
A permeability-aperture relationship is proposed according to the experimental results. Based on this relationship, stress-dependent
permeability in a set of fractures is derived in a three-dimensional domain by using a coupled stress and matrix-fracture
interactive model. A double porosity finite element model is extended by incorporating such stress-dependent permeability
effects. The proposed model is applied to examine permeability variations induced by stress redistributions for an inclined
borehole excavated in a naturally fractured formation. The results indicate that permeability around underground openings
depends strongly on stress changes and orientations of the natural fractures. 相似文献
12.
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. 相似文献
13.
Researchers have recently realized that the non-tectonic natural fractures are developed in shale formations and significant for the exploitation of shale gas. Studies have shown that the tectonic fractures in naturally fractured reservoirs have influences on the maximization of stimulated reservoir volume (SRV) during hydraulic fracturing. However, the effect of the non-tectonic randomly natural fractures on the fracturing network propagation is not well understood. Laboratory experiments are proposed to study the evolution of fracturing network in naturally fractured formations with specimens that contain non-tectonic random fractures. The influences of the dominating factors were studied and analyzed, with an emphasis on natural fracture density, stress ratio, and injection rate. The response surface methodology was employed to perform the multiple-factor analysis and optimization in the maximization of the SRV. A sensitivity study reveals a number of interesting observations resulting from these parameters on the fracturing network evaluation. It is suggested from the geometry morphology of fracturing network that high natural fracture density and injection rate tend to maximize the fracturing network. The influence of stress contrast on fracturing network is nonlinear; an optimal value exists resulting in the best hydraulic fracturing effectiveness. 相似文献
14.
干热岩(HDR)是指不含或仅含少量流体,温度高于180 ℃,其热能在当前技术经济条件下可以利用的岩体。作为一种重要的非常规地热资源,干热岩的开发利用可以借鉴页岩油气的成功经验,采用相似的技术发展路径,找到“地热甜点”,开发出低成本且高效的钻完井技术,逐步形成和完善技术体系,建立与对象相适应的生产运行模式,以期实现对这种巨大资源的有效开发利用。增强型地热系统(EGS)被认为是干热岩资源开采的一种重要方式。EGS最初被称为工程型地热系统,后来才统称为增强型地热系统,是指通过实施特殊的工程工艺,改善地层储集性能或(和)向地层中注入流体,以实现对地热资源的有效开发。其基本方法原理为在干热岩体内钻两口或多口井,将低温流体通过注入井注入干热岩体的天然裂缝系统,或注入通过压裂技术在钻井之间建立的具有水力联系的人工裂缝中加热,通过吸收干热岩内所蕴含的热能,将流体温度提高到一定程度后从生产井采出至地表或近地表进行利用,形成人工热交换系统,用于发电或取暖等。采用EGS技术开发干热岩地热资源,选区选址恰当与否是能否取得成功的最关键环节之一。中深层地热资源可分为水热型和干热岩型两大类、五亚类。其中,干热岩根据其热储孔渗条件差异又可分为无水优储、无水差储和无水无储三亚类,适合作为EGS开发对象的干热岩资源为其中的无水优储和无水差储两种类型。五类地热资源规模呈金字塔形,开发技术难度逐渐增加。基于由热储埋深、热储温度、热储岩性、热储物性、盖层厚度、盖层断裂发育条件等组成的地质资源条件,由钻探成井技术、储层改造技术、管理运营技术组成的工程技术条件,以及由地热需求和资源经济性组成的经济市场条件三个因素,本文建立了三因素分析与多层次指标分解法相结合的干热岩EGS选区评价方法和关键指标,在国内干热岩资源条件较好的17个候选区中,优选出西藏羊八井高温地热区和渤海湾盆地济阳坳陷潜山分布带作为EGS试验有利区。 相似文献
15.
水力压裂是青海共和盆地干热岩地热资源开发的难点技术问题之一。本文基于升级改造的大尺寸真三轴水力压裂物理模拟实验系统模拟干热岩储层高温高压环境,利用青海共和盆地露头岩心进行水力压裂物理模拟实验,揭示干热岩储层水力裂缝的起裂和扩展规律。通过物理模拟实验发现：干热岩储层裂缝起裂可以通过文中提出的起裂模型判断起裂方式和预测起裂压力;水力裂缝在岩石基质中的扩展形态简单,仅沿最大主应力方向延伸;但是水力裂缝会受到岩石中弱面的影响,发生转向沿弱面延伸,形成较复杂的裂缝形态。因此,建议在干热岩储层实际施工中,在天然裂缝发育较丰富的层段开展水力压裂,以实现复杂裂缝网络提取地热能。 相似文献
16.
Most multiple-fractured horizontal wells experience long-term linear flow due to the ultralow permeability of shale gas reservoirs. Considering the existence of natural fractures caused by compression and shear stresses during the process of tectonic movement or the expansion of high-pressure gas, a shale gas reservoir can be more appropriately described by dual-porosity medium. Based on the assumption of slab dual-porosity, this paper uses the trilinear flow model to simulate the transient production behavior of multiple-fractured horizontal wells in shale gas reservoirs, which takes the desorption of adsorbed gas, Knudsen diffusion and gas slippage flow in the shale matrix into consideration. Production decline curves are plotted with the Stehfest numerical inversion algorithm, and sensitivity analysis is done to identify the most influential reservoir and hydraulic fracture parameters. It was found that the density and permeability of the natural fracture network are the most important parameters affecting the production dynamics of multiple-fractured horizontal wells in shale gas reservoirs. The higher the density and permeability of the natural fractures are, the shorter the time is required to exploit the same amount of reserve, which means a faster investment payoff period. The analytical model presented in this paper can provide some insight into the reserve evaluation and production prediction for shale gas reservoirs. 相似文献
17.
为有效模拟裂缝性页岩储层中水力裂缝随机扩展过程,基于单元节点的拓扑数据结构,利用网格节点分裂方式,建立了一种基于有限元网格嵌入零厚度内聚力单元的水力裂缝随机扩展新方法。利用KGD模型解析解和2种室内试验,验证了新方法的准确性和有效性。同时,通过数值算例研究了水平地应力差和储层非均质性对水力裂缝随机扩展过程的影响。研究表明:(1)该方法弥补了ABAQUS平台内置的内聚力单元无法有效模拟水力裂缝随机扩展的不足;(2)在较高水平地应力差下页岩储层非均质性越强,与水力裂缝相交的高角度天然裂缝越容易开启。所建方法能准确地描述复杂水力裂缝的随机扩展行为,可为裂缝性页岩储层的数值模拟提供新手段。 相似文献
18.
在低渗透油藏开发中普遍采用菱形反九点井网,但采油井等缝长压裂将导致非均匀驱替,当油藏物性各向异性时,非均匀驱替现象更加明显。以往针对菱形反九点井网的研究主要基于等缝长压裂,未考虑利用不等缝长设计改善开发效果。根据某各向异性低渗透油藏基本参数建立数值模型,进行不等缝长设计,分析油藏渗透性各向异性及油水井压裂情况对菱形反九点井网水驱效果的影响。随着储层Kx∶Ky的增加,角井裂缝优化穿透比减小,而边井裂缝优化穿透比增大。在相等Kx∶Ky条件下,注水井是否压裂对角井裂缝优化穿透比影响较小,对边井裂缝优化穿透比影响较大。数值模拟结果表明,在各向异性低渗透油藏中进行不等缝长优化设计,能有效改善菱形反九点井网的开发效果,优化结果对于各向异性低渗透油藏的水力压裂设计以及井网布置具有一定的指导意义。 相似文献
19.
裂隙岩体流固耦合问题是目前国内外研究热点之一,采用离散元软件UDEC对裂隙岩体发生节理剪胀的渗透性变化规律进行了模拟分析。基于现场调查的裂隙信息统计生成裂隙网络岩体模型。 通过固定垂直应力、不断增加应力比RS(RS=水平应力/垂直应力)使岩体出现剪胀,采用库伦滑移节理模式对岩体在剪胀过程中的渗透性变化情况进行模拟。结果发现:当应力比较小(RS3.1)时,节理水力隙宽、流速、渗透系数等参数都随着应力比的增加表现出明显的降低; 而当岩体出现剪胀现象之后(应力比大于3.1),发生剪切滑移和剪胀现象的节理控制着裂隙岩体的总体渗流行为,与不考虑节理剪胀的计算结果相比,岩体渗透能力出现了显著增长。这一结果表明,剪胀对裂隙岩体渗透性的影响是显著而不可忽视的。 相似文献
20.
Poro‐mechanical and thermo‐mechanical processes change the fracture aperture and thus affect the water flow pattern in the fracture during the cold water injection into enhanced geothermal systems (EGS). In addition, the stresses generated by these processes contribute to the phenomenon of reservoir seismicity. In this paper, we present a three‐dimensional (3D) partially coupled poro‐thermoelastic model to investigate the poroelastic and thermoelastic effects of cold water injection in EGS. In the model, the lubrication fluid flow and the convective heat transfer in the fracture are modeled by the finite element method, while the pore fluid diffusion and heat conductive transfer in the reservoir matrix are assumed to be 3D and modeled by the boundary integral equation method without the need to discretize the reservoir. The stresses at the fracture surface and in the reservoir matrix are obtained from the numerical model and can be used to assess the variation of in situ stress and induced seismicty with injection/extraction. Application of the model shows that rock cooling induces large tensile stresses and increases fracture conductivity, whereas the rock dilation caused by fluid leakoff decreases fracture aperture and increases compressive total stresses around the injection zone. However, increases in pore pressure reduce the effective stresses and can contribute to rock failure, fracture slip, and microseismic activity. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献