裂缝性碳酸盐岩裂缝的双侧向测井响应特征及解释方法

裂缝评价是裂缝性碳酸盐岩储层评价的关键, 其常规评价方法受到裂缝发育的不均匀性及储层各向异性的影响而存在诸多困难.采用三维数值算法, 利用宏观各向异性地层模型, 研究不同的裂缝参数条件下双侧向测井响应特征, 由此导出一种用于裂缝孔隙度计算的快速算法.分析表明, 裂缝的双侧向响应同裂缝孔隙度与孔隙流体电导率之间存在明显的线性关系, 裂缝的倾角造成双侧向测井曲线幅度差异的变化; 不同倾斜情况下, 将双侧向测井响应近似表示为岩石基岩电阻率、裂缝孔隙度、裂缝流体电导率的函数, 用于裂缝孔隙度的快速计算.实际资料处理表明, 利用双侧向依据该方法确定的裂缝参数同成像测井资料具有良好的对应性.      

煤层气储层裂隙阵列侧向测井响应数值模拟与分析

裂隙是煤层气储层重要的流体渗流通道,是煤层气有效储集体形成的重要条件。研究裂隙的测井响应特征是确定煤层裂隙参数的关键。利用三维有限元方法,确定不同裂隙条件的阵列侧向测井响应。研究表明:阵列侧向测井修正的视电导率与裂隙孔隙度、裂隙流体电导率基本为线性关系;随裂隙倾角增大,阵列侧向视电阻率升高;随基块电阻率增大,阵列侧向视电阻率增加;高角度裂隙(组)的阵列侧向测井深浅视电阻率呈现正差异,而基块电阻率较大的低角度裂隙(组)的响应可呈现负差异特征;相同裂隙孔隙度下,裂隙密度增加可能使裂隙组的阵列侧向视电阻率增加,当裂隙密度足够大时,测井响应基本不随裂隙密度发生变化;交叉裂隙的测井响应与平行裂隙组差异明显,且不同倾角组合的交叉裂隙的测井响应之间差异明显。      

裂缝性储层裂缝的阵列侧向测井响应数值模拟

为了研究裂缝性储层裂缝的阵列侧向测井响应, 采用裂缝组平板裂缝模型, 利用三维有限元方法, 进行不同裂缝参数条件下的高分辨率阵列侧向测井响应数值模拟.表明阵列侧向视电导率和裂缝孔隙度与裂缝流体电导率成线性关系; 低角度裂缝的阵列侧向深浅视电阻率呈现负差异特征, 高角度缝呈现正差异; 相同裂缝孔隙度下, 裂缝密度增加可使得阵列侧向视电阻率增加, 当裂缝密度足够大时, 测井响应特征基本不随裂缝密度发生变化; 交叉缝的测井响应与具有相同角度的平行裂缝组差异明显, 且不同倾角组合的交叉缝的测井响应之间差异明显.阵列侧向测井对裂缝有较强的敏感性, 根据其测井视电阻率大小可以评价裂缝孔隙度, 根据深浅视电阻率的幅度差异可以研究裂缝的产状, 这些为将阵列侧向测井用于裂缝性储层裂缝评价有指导意义.      

应用随裂缝倾角变化的孔隙模型改进储层饱和度计算

复杂的孔隙结构及多变的孔隙类型使得裂缝性储层成为测井评价的难点,通过假定不同孔隙类型间的导电方式来建立合理的孔隙模型,进而结合阿尔奇方程实现饱和度计算是目前裂缝性储层定量评价的主要方法和途径。在经典的双孔隙模型中求取岩石系统的总孔隙指数m时,未考虑裂缝与电流线的夹角变化带来的影响。通过在裂缝性储层中同时考虑裂缝曲折度、裂缝与电流线夹角双重因素,构建了改进双重孔隙模型,并结合电成像资料提供的倾角信息进行了实例处理。结果表明,改进模型的计算值与取心分析有着更好的一致性;若不考虑裂缝与电流线夹角等因素,将带来321%的相对误差。     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

基于RJM断续节理岩体强度与破坏特征的数值模拟研究

裂隙结构的存在对于工程岩体的强度和稳定性具有重要影响,岩石宏观裂隙的产生源自于微破裂的积累。针对岩体裂隙的粗糙特性,通过Matlab建立考虑粗糙度的节理模型（Roughness Joint Model）,采用简化的正弦曲线来表示粗糙节理,并将其导入到颗粒流试验模型中进行单轴压缩试验。对比完整岩体、直线型裂隙岩体、RJM岩体三者破坏的应力-应变曲线,改变裂隙倾角（与水平方向夹角）α,岩桥倾角β,裂隙密度γ,建立不同裂隙分布的断续节理岩体数值试样,开展一系列数值模拟试验。研究结果发现:（1）裂隙的存在明显降低了岩体的抗压强度,RJM模型峰值强度和峰值应变均高于直线型裂隙岩体;（2）岩体抗压强度总体上随裂隙倾角增大而增加,随裂隙密度增加而减小,但随岩桥倾角的改变呈非线性变化,岩桥倾角45°时峰值强度最低,峰值应变最小;（3）裂隙分布会影响岩体的破裂模式,微裂隙的扩展反映了岩体力学性质的各向异性;（4）不同倾角下增加裂隙密度,岩体强度下降程度不同,倾角75°时密度对强度影响最小,30°和60°时影响最大。     

井斜角对裂缝特征参数统计的影响及其意义

井轨迹对裂缝钻遇率有重要的影响。以某裂缝性碳酸盐岩油藏为例,利用100余口井的FMI资料按照小斜度井（包括直井）和大斜度井（包括水平井）分组分别统计裂缝的平均倾角、裂缝密度参数,并与不考虑井斜角利用所有FMI资料统计的结果相互对比,以分析井斜角对裂缝倾角、裂缝密度参数统计的影响。结果表明,按小斜度井钻遇裂缝统计的裂缝平均倾角和裂缝密度均明显低于按大斜度井和所有井未分组钻遇裂缝的统计结果;因此,小斜度井钻遇高角度裂缝的几率明显低于大斜度井。这一结论具有三方面的意义：对于本文中某裂缝性碳酸盐岩油藏,目前主要为直井和小斜度井钻遇的G-H层,应部署水平井进一步评价其油气勘探潜力;为更准确地统计高角度裂缝性储层的裂缝倾角并客观地评价裂缝密度,应尽量利用井轨迹垂直于裂缝走向和主要裂缝面的水平井钻遇的裂缝来分析统计;对于同类高角度裂缝性油藏的勘探评价和开发,应尽量采用大斜度井和水平井并确保井轨迹垂直于主要裂缝面。     

库车坳陷构造裂缝发育特征及分布规律

通过对野外露头和井下岩心构造裂缝的观察和描述,及其性质、产状、密度、强度、开度、充填程度和充填物等测 量要素的观测和统计,对库车前陆冲断带山前构造裂缝发育特征及分布规律进行了综合分析。该区主要发育三组近南北走 向的剪裂缝和张剪裂缝,裂缝倾角总体以大于45°的高角度缝和垂直缝为主,野外和岩心裂缝开度分别以0~5 mm和0~1mm 为主,裂缝多未充填或半充填,充填物以钙质充填为主。裂缝各参数的分布呈现出东西分段,南北分带的特点,以克拉苏 -依奇克里克构造带的裂缝密度和强度最高。具体到层位上,库车东部以下侏罗统阿合组裂缝最为发育,库车西部则以下 白垩统巴什基齐克组砂岩的裂缝密度最高。总体上,库车坳陷构造裂缝十分发育,有益于致密储层孔渗条件的改善,有利 于裂缝型油气藏的形成。     

Petrophysical characterization and natural fracturing in an olivine-dolerite aquifer

As part of a study investigating the naturally-occurring fractures in mafic rocks, two holes were drilled 450 m apart through the Palisades dolerite sill in New York. Well-2 is 229 m deep and Well-3 was drilled to 305 m, both penetrating through the sill and into the underlying Triassic sediments of the Newark Basin. Both holes were logged with downhole geophysical tools, including the BHTV, which acoustically images fractures intersecting the well. Understanding the fracture pattern, density, and porosity in the sill is essential for identifying possible zones of active fluid flow and high permeability. Using the BHTV logs, 96 and 203 fractures were digitally mapped within the sill in Well-2 and Well-3, respectively. Most fractures appear to dip steeply (76-78°). There is a shift in fracture orientation, however, and these fractures may or may not be continuous over the short lateral distance between Well-2 and Well-3. The lithology of the sill as identified by drill chips is nevertheless continuous between the holes. Both intersect a 7 m thick olivine-rich layer about 15 m above the bottom of the sill. Several fractures identified in Well-2 have large apparent aperture (>6cm) which correspond to high porosity zones (6-14%) observed in the logs. Resistivity logs were used to compute porosity using Archie's law and match well with the neutron porosity log in Well-2. We use the relationship between porosity and fracture aperture within the sill at Well-2 to infer the porosity in Well-3. High-porosity, large-aperture zones, including the target olivine layer, are identified in both holes. Changes in the temperature gradient log indicate active fluid flow in the sill, although flow appears to be most active in the sediments. Direct field measurements of bulk permeability, hydrologic modeling of fluid flow and calibration of fracture and log porosity will be undertaken in the future.     

Application of Stochastic Fracture Network with Numerical Fluid Flow Simulations to Groundwater Flow Modeling in Fractured Rocks

ＩＮＴＲＯＤＵＣＴＩＯＮＧｒｏｕｎｄｗａｔｅｒｏｒｆｌｕｉｄｆｌｏｗｍｏｄｅｌｉｎｇｉｎｆｒａｃｔｕｒｅｄｒｏｃｋｓｉｓａｃｏｍｐｌｉｃａｔｅｄｔｈｅｏｒｅｔｉｃａｌａｎｄａｐｐｌｉｅｄｔｏｐｉｃ.Ｂｏｔｈｔｈｅｏｒｅｔｉｃａｌｌｙａｎｄｏｐｅｒａｔｉｏｎａｌｌｙ ,ｉｔｉｓｉｍｐｏｒｔａｎｔｉｎｍａｎｙｆｉｅｌｄｓｓｕｃｈａｓｇｅｏｌｏｇｉｃａｌａｎｄｈｙｄｒｏｇｅｏｌｏｇｉｃａｌｅｎｇｉｎｅｅｒｉｎｇ ,ｅｎｖｉｒｏｎｍｅｎｔａｌｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｅｔｒｏｌｅｕｍｅｎｇｉｎｅｅｒｉｎ…     

基于CT-CA的天然裂隙岩石受拉破坏特性研究

岩石内部存在的裂隙、孔洞等天然损伤对岩石的力学性能和破坏过程有重要影响,依据细胞自动机理论结合CT无损识别技术实现了含天然裂隙岩石在劈裂条件下裂纹扩展和贯通全过程及其力学性能变化规律的研究。从裂隙砂岩的真实细观结构出发,构建了天然裂隙岩石的数值计算模型,运用CASRock数值计算软件完成了含不同裂隙倾角的砂岩劈裂破坏的数值试验,分析了裂隙倾角对砂岩的力学特性、裂纹扩展过程及能量演化的影响规律。研究表明：(1)天然裂隙砂岩的抗拉强度与裂隙倾角密切相关,随着裂隙倾角的增加,其抗拉强度呈现先减小后增加的趋势;(2)裂隙起裂于天然裂隙尖端,当裂隙倾角0°≤θ<48°时,岩样的破坏是由错开型裂纹引起,裂纹沿着与天然裂隙近垂直方向扩展;当裂隙倾角48°≤θ<94°时,岩样的破坏是由张开型裂纹引起,裂纹沿着与天然裂隙近平行方向扩展;(3)劈裂过程中裂纹尖端应力场存在拉应力区和压应力区,拉应力造成翼裂纹由天然裂隙尖端沿加载端方向萌生扩展,而压应力则引发次生裂纹沿天然裂隙方向扩展;(4)含天然裂隙砂岩劈裂破坏过程能量演化可划分为4个阶段,随裂隙倾角的增大,峰值点处的总能量密度、弹性能密度先缓...     

构造电阻率差比值法识别火山岩裂缝地层天然气层

在基岩电阻率较高的硬地层中,除去泥质、孔隙和其他矿物等因素的影响后,地层电阻率与致密围岩电阻率的差异就被认为是裂缝以及地层孔隙中储存的流体性质引起的。选择相同岩性含气储层与致密围岩层的电阻率值,利用数值反演的方法确定出计算含气层段消除裂缝以及孔隙中流体性质影响以后的地层真电阻率公式,定义地层真电阻率和深侧向电阻率的差值与地层深侧向电阻率的比值为构造电阻率差比值,该参数主要反映裂缝以及孔隙中的流体性质对电阻率降低幅度的影响。因此,利用差比值法可以识别火山岩裂缝地层的流体性质,进而制作油层、气层和水层的判别图版,并结合其他测井曲线、油藏动态资料以及气测信息综合识别火山岩裂缝地层天然气层。该方法在准噶尔盆地研究区火山岩天然气层的识别中取得了很好的应用效果,解释结论与试油结果基本吻合。     

东准噶尔蒙西斑岩铜钼矿床脉体特征及其形成机制

蒙西斑岩铜钼矿床位于东准噶尔伊吾县琼河坝花岗岩北侧,以发育细脉、网脉状矿化为特征。脉体类型多样,包括石英脉、石英硫化物脉和硫化物脉等。根据脉体力学成因机制,其又可划分为水压破裂充填脉和构造破裂充填脉。构造破裂充填脉体按破裂形成的位错特征有正断与逆断两种。地表石英脉产状陡立,明显受断裂控制。脉体体积分数统计结果显示流体富集区呈向北缓倾的带状分布于深100～400m范围内,并与矿化富集带有较好的对应关系。脉体富集带内,脉体以共轭形式存在,一组为倾角较小的逆断破裂充填脉,另一组为倾角较大的正断破裂充填脉,他们可能为缓倾逆冲剪切带的次级破裂与充填脉,即富矿带内脉体是沿矿区低角度逆冲断裂次级破裂面充填的。矿区地表及深部(400m以下)脉体以陡立为主,矿化较弱。矿区流体的运移具先沿水压直立破裂往上运移,进入剪切带后沿剪切带次级破裂侧向和向上运移,并在剪切带中富集成矿。低温矿物组合脉体穿切高温矿物组合脉体的特征说明脉体形成过程矿区处于隆升构造环境,这对斑岩铜矿成矿有利。     

A fracture mapping and extended finite element scheme for coupled deformation and fluid flow in fractured porous media

This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd.     

断层相关褶皱概念模型中的裂缝域

断层相关褶皱对裂缝的控制作用主要体现在岩层的几何形态对裂缝方位的影响和变形对岩石破裂的控制程度。本文在断层相关褶皱概念模型的基础上提出了裂缝域,并以Suppe的平行膝折概念模型为对象,利用几何学方法对模型中区域性裂缝的裂缝域进行分区描述,并提出了断层相关褶皱的曲线概念模型中区域性裂缝的分布规律。以新疆卡普沙良河一处断层转折褶皱剖面上裂缝的实测数据为依据,对模型进行验证,发现测量结果同模型基本一致。该方法为前陆冲断带地区断层相关褶皱对裂缝控制作用的定量预测提供了一种全新的思路。     

Stress-Dependent Fluid Flow and Permeability in Fractured Media: from Lab Experiments to Engineering Applications

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.     

Mechanical and hydraulic properties of rocks related to induced seismicity

Witherspoon, P.A. and Gale, J.E., 1977. Mechanical and hydraulic properties of rocks related to induced seismicity. Eng. Geol., 11(1): 23–55.The mechanical and hydraulic properties of fractured rocks are considered with regard to the role they play in induced seismicity. In many cases, the mechanical properties of fractures determine the stability of a rock mass. The problems of sampling and testing these rock discontinuities and interpreting their non-linear behavior are reviewed. Stick slip has been proposed as the failure mechanism in earthquake events. Because of the complex interactions that are inherent in the mechanical behavior of fractured rocks, there seems to be no simple way to combine the deformation characteristics of several sets of fractures when there are significant perturbations of existing conditions. Thus, the more important fractures must be treated as individual components in the rock mass.In considering the hydraulic properties, it has been customary to treat a fracture as a parallel-plate conduit and a number of mathematical models of fracture systems have adopted this approach. Non-steady flow in fractured systems has usually been based on a two-porosity model, which assumes the primary (intergranular) porosity contributes only to storage and the secondary (fracture) porosity contributes only to the overall conductivity. Using such a model, it has been found that the time required to achieve quasi-steady state flow in a fractured reservoir is one or two orders of magnitude greater than it is in a homogeneous system. In essentially all of this work, the assumption has generally been made that the fractures are rigid.However, it is clear from a review of the mechanical and hydraulic properties that not only are fractures easily deformed but they constitute the main flow paths in many rock masses. This means that one must consider the interaction of mechanical and hydraulic effects. A considerable amount of laboratory and field data is now available that clearly demonstrates this stress-flow behavior. Two approaches have been used in attempting to numerically model such behavior: (1) continuum models, and (2) discrete models. The continuum approach only needs information as to average values of fracture spacing and material properties. But because of the inherent complexity of fractured rock masses and the corresponding decrease in symmetry, it is difficult to develop an equivalent continuum that will simulate the behavior of the entire system. The discrete approach, on the other hand, requires details of the fracture geometry and material properties of both fractures and rock matrix. The difficulty in obtaining such information has been considered a serious limitation of discrete models, but improved borehole techniques can enable one to obtain the necessary data, at least in shallow systems. The possibility of extending these methods to deeper fracture systems needs more investigation. Such data must be considered when deciding whether to use a continuum or discrete model to represent the interaction of rock and fluid forces in a fractured rock system, especially with regard to the problem of induced seismicity. When one is attempting to alter the pressure distribution in a fault zone by injection or withdrawal of fluids, the extent to which this can be achieved will be controlled in large measure by the behavior of the fractures that communicate with the borehole. Since this is essentially a point phenomenon, i.e., the changes will propagate from a relatively small region around the borehole, the use of a discrete model would appear to be preferable.     

Naturally propped fractures caused by quartz cementation preserve oil reservoirs in basement rocks

Much silica precipitation in oil reservoirs occurred in the presence of hydrocarbons, evidenced by the entrapment of oil fluid inclusions in quartz. Also, silica in sedimentary basins is commonly precipitated at oil‐window temperatures. This spatial and temporal relationship between oil and quartz precipitation aids the entry of oil into fractured reservoirs, including fractured basement. Where quartz is precipitated as fracture linings, the fractures are propped open by bridging quartz crystals, creating high fracture porosity and permeability. Evidence from fossil fractured reservoirs shows a large proportion of oil residue is in such propped open fractures.     

高黏粒含量海洋土电阻率特征分析及模型构建

在高黏粒含量海洋土中,黏粒吸附孔液中阳离子形成的双电层会直接影响土体导电性,与孔液、土体结构共同形成土体导电性的影响因素,使导电结构变得复杂。对饱和的海洋土而言,固结程度不同代表了土体结构不同。为量化分析双电层对不同固结程度高黏粒含量海洋土导电性的影响,界定电阻率测试技术对海洋土结构变化的反映能力,分别以黄河三角洲和南海北部陆坡海洋土为研究对象,进行电阻率和其他物理性质测试。基于测试结果分析不同黏粒含量海洋土电阻率的变化规律和独有特征,借鉴黏性土多元导电理论构建适合的电阻率模型。研究表明,黏粒含量对不同固结程度海洋土电阻率的影响不同,根据其变化特征可分为3个阶段：正常固结的海洋土电阻率变化可用简化的二元模型公式表示,当孔隙度小于50%时,电阻率随孔隙度呈明显的幂函数降低,当孔隙度介于50%~60%时,电阻率随孔隙度降低速率明显减小;欠固结海洋土孔隙度均高于60%,电阻率随孔隙度线性减小,可用线性模型公式描述。电阻率测量技术对高黏粒含量正常固结海洋土结构变化（孔隙度变化）反映能力良好,探测灵敏度可达25 Ω·m;对欠固结海洋土结构变化反映困难。     

裂隙岩体渗透系数确定方法综述

总结了近年来裂隙岩体渗透系数的确定方法,主要有现场水力试验法、裂隙测量法和离散裂隙网络渗流数值试验法。每种方法都有自己的适用性和测量尺度。裂隙岩体渗透系数存在尺度效应,针对不同尺度的研究对象,应尽量选取与渗流模型网格剖分尺度匹配的测量方法。裂隙的延伸具有方向性,测试点的布设应合理科学,不同测量方法需要结合起来才能得到裂隙岩体真实的渗透系数。