成像测井缝洞解释评价研究进展

综合露头刻度井、数值模拟和物理模拟的实验结果对微电阻率扫描成像测井缝洞解释评价的研究进展进行了总结和探讨,分别包括岩芯归位、缝洞刻度率、岩芯和成像裂缝参数的差异性,以及裂缝宽度和地层孔隙度的计算。岩芯归位确定了缝洞体在地层中的深度和方位。电成像测井对过井眼未充填缝和孔洞发育段的刻度率为100%,而单个孔洞在岩芯和电成像测井中通常无法一一对应。单条裂缝的倾角或井周长度在岩芯和成像测井中具有较高的相关性,而裂缝宽度的相关性一般。数模和物模的结果表明电成像测井裂缝宽度的计算受多种因素影响,当裂缝的视宽度大于0.1 mm时,电成像测井计算的裂缝宽度和裂缝视宽度基本一致,当小于0.1 mm时,电成像测井计算的宽度值误差较大。目前电成像测井通过孔隙度谱法和图像分割法计算地层的孔隙度,但都受限于除孔洞的其他低阻地质体的影响,因此需要加强岩芯刻度的功能。     

成像测井缝洞解释评价研究进展

综合露头刻度井、数值模拟和物理模拟的实验结果对微电阻率扫描成像测井缝洞解释评价的研究进展进行了总结和探讨,分别包括岩芯归位、缝洞刻度率、岩芯和成像裂缝参数的差异性,以及裂缝宽度和地层孔隙度的计算.岩芯归位确定了缝洞体在地层中的深度和方位.电成像测井对过井眼未充填缝和孔洞发育段的刻度率为100％,而单个孔洞在岩芯和电成像...     

融合地震属性和成像测井信息优化的储层缝洞带评价与钻井轨迹

针对新疆某工区强非均质性岩溶缝洞型碳酸盐岩储层,提出了一种基于地震多属性优化与电成像测井信息融合的碳酸盐岩储层缝洞带识别评价与钻井轨迹优化设计方法。首先,利用相干体、蚂蚁体等多种地震属性与张量约束的稀疏脉冲反演结果综合实现缝洞发育带空间分布预测与雕刻,确定可能油气圈闭的分布范围;然后,基于数学形态学滤波方法,由电成像测井资料自动拾取井壁的裂缝和溶蚀孔洞分布,统计井壁裂缝倾向和倾角;再次,用工区内钻孔多臂井径资料计算的井筒崩落信息与电成像资料拾取的诱导缝信息预测工区的最大和最小主应力方向;最后,融合地震属性与电成像测井信息评价的缝洞分布结果,预测工区内走滑断裂带的最大主应力方向以NE—SW为主,局部为NW—SE,且以中高角度为主的裂缝走向与主应力近乎平行。采用预测主应力分布方法实现了沿NW—SE方向的大斜度井与水平井优化设计,钻遇了大规模优质储集体,获得了稳定油流,并以此为基础在周围打井若干口,极大地降低了深层勘探的钻井风险。     

电成像测井资料在裂缝成因分析中的应用

以塔里木盆地为例,探讨高分辨率电成像测井资料在裂缝成因分析中的应用。根据电成像测井图像解释,利用裂缝图像的颜色(或灰度)和产状特征,结合岩芯标定和常规测井资料,可以有效地区分构造裂缝与非构造裂缝、张裂缝与剪裂缝,确定裂缝充填与否和充填物的成分,以及裂缝的溶蚀改造程度。利用从电成像测井图像上提取的共轭裂缝的产状信息,还可以恢复裂缝形成时古构造应力场的方向。     

Estimation of 2-D and 3-D Fracture Densities from 1-D Data - Experimental and Field Results

2-D and 3-D densities of fractures are commonly used in mining safety design, natural gas and oil production in fractured reservoirs, and the characterization of subsurface flow and transportation systems in fractured rocks. However, many field data sets are collected in 1-D frequency (f) (e.g., scanlines and borehole data). We have developed an ARC/ INFO-based technology to calculate fracture frequency and densities for a given fracture network. A series of numerical simulations are performed in order to determine the optimal orientation of a scanline, along which the maximum fracture frequency of a fracture network can be obtained. We calculated the frequency (f) and densities (both D1 and D2) of 36 natural fracture trace maps, and investigated the statistical relationship between fracture frequency and fracture density D1, i.e. D1=1.340f+ 0.034. We derived analytical solutions for converting dimensional density (D1) to non-dimensional densities (D2 and D3) assuming that fracture length distribution f     

Fractal and Multifractal Properties of the Spatial Distribution of Natural Fractures--Analyses and Applications

The Cantor's dust theory is applied to investigate the scaling properites of the spatial distribution of natural fractures obtained from detailed scanline surveys of 27 field sites in the Appalachian Plateau of western New York, USA. The results obtained in this study indicate: 1) fracture spacing is characterized by fractal and multifractal properties. On small scales analyses yield an average fractal dimension of 0.15, which suggests a very high degree of clustering. In contrast, on large scales, fractures tend to be more regular and evenly distributed with an average fracture dimension of 0.52; 2) fractal dimension varies with different sets in different orientations, which can be attributed to interactions between pre-existing fractures and younger ones, as well as variations of the intensity of the stresses under which the fractures were formed; 3) a time sequence of fracture set formation can be proposed based on fractal and multifractal analyses, which consists of (from old to young): N-S, NW, ENE, and NE-striking sets. This time sequence is confirmed by the study of the abutting relationships of different fracture sets observed in the field.     

Fractal and Multifractal Propertiesof the Spatial Distribution of Natural Fractures——Analyses and Applications

The Cantor's dust theory is applied to investigate the scaling properties of the spatial distribution of natural fractures obtained from detailed scanline surveys of 27 field sites in the Appalachian Plateau of western New York, USA. The results obtained in this study indicate: 1) fracture spacing is characterized by fractal and multifractal properties. On small scales analyses yield an average fractal dimension of 0.15, which suggests a very high degree of clustering. In contrast, on large scales, fractures tend to be more regular and evenly distributed with an average fracture dimension of 0.52; 2) fractal dimension varies with different sets in different orientations, which can be attributed to interactions between pre-existing fractures and younger ones, as well as variations of the intensity of the stresses under which the fractures were formed; 3) a time sequence of fracture set formation can be proposed based on fractal and multifractal analyses, which consists of (from old to young): N-S, NW, EN     

Application of X-ray computed tomography for analyzing cleat spacing and cleat aperture in coal samples

New techniques to determine distributions of cleat aperture, cleat orientation and cleat spacing from CT scans have been developed. For cleat orientation and spacing distributions, two different coal blocks were scanned. The CT scans have been analyzed for the three orthogonal directions. Histograms of the cleat orientations are bimodal, expressing the typical cleat texture of face and butt cleats and bedding perpendicular relaxation fractures. Deviations up to 20° from the peak values in the cleat orientation distributions were used as input for automated image analysis of cleat spacing. Distributions of the cleat spacing measurements are related to the face and butt cleat directions. The term “relevant cleat length” is introduced as a measure to extract the amount of cleat length involved with the cleat spacing measurements. The ratio ranges from 0.03 to 0.38 and expresses the difference in cleat texture in both samples. Cleat spacing versus relevant cleat length shows sample specific patterns for face cleat, butt cleat and bedding. To describe cleat aperture quantitatively, peak height and missing attenuation have been used. The image of a cleat was seen as a convolution of a rectangular fracture profile with a Gaussian point spread function.     

Fracture sampling on a cylinder: From scanlines to boreholes and tunnels

Summary In this paper we develop correction factors for orientation bias when sampling rock fractures with a borehole. For a borehole with zero radius we obtain Terzaghi's result for scanline sampling: otherwise the appropriate correction factor may differ significantly from that obtained for the scanline model, depending on the relative size of the borehole and the rock fracture. Two separate cases are discussed: the sampling of fractures with a borehole and sampling with the surface of a borehole. Analytical expressions for the correction factors are graphed and tabulated.     

Numerical evaluation of the suitability of the equivalent porous medium model for characterizing the two-dimensional flow field in a fractured geologic medium

The equivalent porous medium (EPM) method is an efficient approximation method for groundwater yield analysis considering the equivalent permeability in a fractured geologic medium (FGM). The EPM method is widely used in many practical hydrogeological problems from local to regional scales. However, when calculating water head and velocity distributions, the suitability of the EPM model remains insufficiently evaluated. The suitability refers to the head error caused by the application of the EPM model. The smaller the error, the better the suitability. In this study, the influence of fracture geometric attributes on the suitability was quantitatively studied in numerical simulation experiments, and the EPM model simulation results were compared to those obtained with the discrete fracture network (DFN) model. The results indicated that the suitability decreased with increasing fracture spacing. When the fracture spacing was smaller than 0.6 m, the influence of an increase in the fracture spacing on the suitability was obvious. For the same fracture spacing, the suitability generally increased with increasing trace length. When the spatial variability of fracture aperture is not considered, the change of the fracture aperture did not impact the suitability. The fracture orientation slightly impacted the suitability, which can be ignored. An evaluation standard based on fracture parameters is proposed to estimate the suitability of the EPM model, which provides a scientific basis to ascertain whether this model can be applied to a given site to solve the head-related hydrogeological issues with FGMs.

     

Calculation of fracture parameters and their effect on porosity and permeability using image logs and petrophysical data in carbonate Asmari reservoir,SW Iran

Evaluation of fractures and their parameters, such as aperture and density, is necessary in the optimization of oil production and field development. The purpose of this study is the calculation of fracture parameters in the Asmari reservoir using two electrical image logs (FMI, EMI), and the determination of fracture parameters’ effect on the porosity and permeability using thin sections and velocity deviation log (VDL). The results indicate that production in the Asmari reservoir is a combination of fractures and rock matrix. Fracture aperture (VAH) and fracture porosity (VPA) are only measurable with core and image logs directly. However, regarding core limitations, the image log has been recognized as the best method for fracture parameter determination due to their high resolution (2.5 mm). In this study, VDL log and thin sections have been used as auxiliary methods which may be available in all wells. The VDL log provides a tool to obtain downhole information about the predominant pore type in carbonates. Results indicate that between fracture parameters, VAH is considered as the most important parameter for determining permeability. For well No. 3, VAH ranges from minimum 51 × 10^?5 mm to maximum 0. 047 mm and VPA changes from min 10^?5% to maximum 0.02056%. For well No. 6, VAH varies from 5 × 10^?4 to 0.0695 mm and VPA varies from 10^?5 to 0.015%. Therefore, due to high fracture density and fracture aperture, it seems that most of effective porosity originates from fractures especially in well No. 3. However, VDL for well No. 6 indicates that intercrystalline and vuggy porosity are the dominant porosity. This result may be an indication for fracture set diversity in the two studied wells. While in well No. 3, they related to the folding and active faults, in well No. 6 they are only of folding type. Furthermore, results indicate the high capability for both of EMI and FMI image logs for calculation of fracture and vug parameters in the carbonate reservoirs.     

致密低渗透砂岩储层多尺度天然裂缝表征方法

本文利用地表露头、岩芯、测井、三维地震和微观分析资料,根据多尺度裂缝的发育特征及不同资料的精度,探讨了多尺度裂缝的表征方法。致密低渗透砂岩储层发育多尺度天然裂缝,根据天然裂缝规模以及限制天然裂缝扩展延伸的岩石力学层界面,将天然裂缝分为大尺度裂缝、中尺度裂缝、小尺度裂缝和微尺度裂缝4级,不同尺度裂缝的发育特征存在较大差异。其中大尺度裂缝主要应用三维地震资料的叠后属性进行检测和表征;中小尺度裂缝在单井上可以应用岩芯、成像测井和常规测井相结合的方法进行表征,在平面上采用基于地质和测井约束的三维地震叠前各向异性分析以及储层地质力学相结合的方法进行综合预测和评价;微尺度裂缝一般应用铸体薄片、普通薄片、扫描电镜和三维CT扫描等微观分析方法进行观察描述,通过主控因素约束的方法进行预测。利用本文提出的多尺度裂缝表征方法实现了对鄂尔多斯盆地西南部致密砂岩储层多尺度裂缝表征,可为其它致密低渗透砂岩储层多尺度裂缝表征提供借鉴。     

Use of Descartes Folium Equation for Deriving a Relation between Total Aperture of Fractures after Uniaxial Compression and Strain Parameters of Different Rocks Exhibiting Negative Total Volumetric Strains

The axial, crack and total volumetric strains, porosity, elastic constants, crack damage stresses, uniaxial compressive strengths, as well as fracture apertures and number of fracture traces in rock samples surface after compression were defined for different chalk, basalt, dolomite, granite, limestone and sandstone samples exhibiting negative total volumetric strain at failure. It is established that the total (summed) aperture of vertical fractures obtained on the lateral surface of rock sample is related to three characteristic strain parameters: axial strain at the onset of negative total volumetric strain, axial failure strain and negative total volumetric strain at failure. The relation is based on Descartes folium equation, where the length of the loop of folium is equal to axial strain coordinate at the onset of negative total volumetric strain. This relation shows that the total aperture increases according to power law with increasing difference between axial failure strain and axial strain at the onset of negative total volumetric strain. Simultaneously, an increase in this difference leads to an increase in the value of negative total volumetric strain at failure. It is found that a direct correlation between total aperture of fractures and negative total volumetric strain at failure is relatively weak. Nevertheless, total aperture of fractures tends to increase with increasing absolute value of negative total volumetric strain at failure. It is revealed that there is no connection between the number of fracture traces and negative total volumetric strain at failure.     

致密低渗透砂岩储层多尺度天然裂缝表征方法

本文利用地表露头、岩芯、测井、三维地震和微观分析资料,根据多尺度裂缝的发育特征及不同资料的精度,探讨了多尺度裂缝的表征方法。致密低渗透砂岩储层发育多尺度天然裂缝,根据天然裂缝规模以及限制天然裂缝扩展延伸的岩石力学层界面,将天然裂缝分为大尺度裂缝、中尺度裂缝、小尺度裂缝和微尺度裂缝4级,不同尺度裂缝的发育特征存在较大差异。其中大尺度裂缝主要应用三维地震资料的叠后属性进行检测和表征;中小尺度裂缝在单井上可以应用岩芯、成像测井和常规测井相结合的方法进行表征,在平面上采用基于地质和测井约束的三维地震叠前各向异性分析以及储层地质力学相结合的方法进行综合预测和评价;微尺度裂缝一般应用铸体薄片、普通薄片、扫描电镜和三维CT扫描等微观分析方法进行观察描述,通过主控因素约束的方法进行预测。利用本文提出的多尺度裂缝表征方法实现了对鄂尔多斯盆地西南部致密砂岩储层多尺度裂缝表征,可为其它致密低渗透砂岩储层多尺度裂缝表征提供借鉴。     

Cubic law with aperture-length correlation: implications for network scale fluid flow

Previous studies have computed and modeled fluid flow through fractured rock with the parallel plate approach where the volumetric flow per unit width normal to the direction of flow is proportional to the cubed aperture between the plates, referred to as the traditional cubic law. When combined with the square root relationship of displacement to length scaling of opening-mode fractures, total flow rates through natural opening-mode fractures are found to be proportional to apertures to the fifth power. This new relationship was explored by examining a suite of flow simulations through fracture networks using the discrete fracture network model (DFN). Flow was modeled through fracture networks with the same spatial distribution of fractures for both correlated and uncorrelated fracture length-to-aperture relationships. Results indicate that flow rates are significantly higher for correlated DFNs. Furthermore, the length-to-aperture relations lead to power-law distributions of network hydraulic conductivity which greatly influence equivalent permeability tensor values. These results confirm the importance of the correlated square root relationship of displacement to length scaling for total flow through natural opening-mode fractures and, hence, emphasize the role of these correlations for flow modeling.     

Permeability tensor and representative elementary volume of fractured rock masses

Based on a simulation of three-dimensional fracture networks and a superposition principle of liquid dissipation energy for fractured rock masses, a model of the fracture permeability tensor is proposed. An elastic constitutive model of rock fractures, considering fracture closure and dilation during shearing, is also proposed, based on the dilation angle of the fracture. Algorithms of flow-path searching and calculation of the effective flow coefficients for fracture networks are presented, together with a discussion on the influence of geometric parameters of the fractures (trace length, spacing, aperture, orientation and the number of fracture sets) on magnitude, anisotropy of hydraulic permeability and the size of a representative elementary volume (REV). The anisotropy of hydraulic permeability of fractured rock masses is mainly affected by orientation and the number of fracture sets, and the REV size is mainly influenced by trace length, spacing and the number of fracture sets. The results of studies on REV size and the influence of in-situ stress on hydraulic conductivity of the rock mass on the slope of Jinping-I hydropower station, China, are presented using the developed models and methods. The simulation results agreed well with the results obtained from field water-pressure measurements, with an error of less than 10 %.     

Characterizing Fracture Spatial Patterns by Using Semivariograms

Semivariogram is applied to fracture data obtained from detailed scanline surveys of nine field sites in western New York, USA in order to investigate the spatial patterns of natural fractures. The length of the scanline is up to 36 m. How both fracture spacing and fracture length vary with distance is determined through semivariogram calculations. In this study, the authors developed a FORTRAN program to resample the fracture data from the scanline survey. By calculating experimental semivariogram, the authors found five different types of spatial patterns that can be described by linear, spherical, reversed spherical, polynomial I (for a<0) and polynomial II (for a>0) models, of which the last three are newly proposed in this study. The well-structured semivariograms of fracture spacing and length indicate that both the location of the fractures and the length distribution within their structure domains are not random. The results of this study also suggest that semivariograms can provide useful infor     

Isotopic effects in fracture-dominated reactive fluid-rock systems

A mathematical model is presented that describes the effects of pore fluid aqueous diffusion and reaction rate on the isotopic exchange between fluids and rocks in reactive geo-hydrological systems where flow is primarily through fractures. The model describes a simple system with parallel equidistant fractures, and chemical transport in the matrix slabs between fractures by aqueous diffusion through a stagnant pore fluid. The solid matrix exchanges isotopes with pore fluid by solution-precipitation at a rate characterized by a time constant, R (yr⁻¹), which is an adjustable parameter. The effects of reaction on the isotopes of a particular element in the fracture fluid are shown to depend on the ratio of the diffusive reaction length for that element (L) to the fracture spacing (b). The reaction length depends on the solid-fluid exchange rate within the matrix, the partitioning of the element between the matrix pore fluid and the matrix solid phase, the porosity and density of the matrix, and the aqueous diffusivity. For L/b < 0.3, fluid-rock isotopic exchange is effectively reduced by a factor of 2L/b relative to a standard porous flow (single porosity) model. For L/b > 1, the parallel fracture model is no different from a porous flow model. If isotopic data are available for two or more elements with different L values, it may be possible to use the model with appropriate isotopic measurements to estimate the spacing of the primary fluid-carrying fractures in natural fluid-rock systems. Examples are given using Sr and O isotopic data from mid-ocean ridge (MOR) hydrothermal vent fluids and Sr isotopes in groundwater aquifers hosted by fractured basalt. The available data for MOR systems are consistent with average fracture spacing of 1-4 m. The groundwater data suggest larger effective fracture spacing, in the range 50-500 m. In general, for fractured rock systems, the effects of fracture-matrix diffusive exchange must be considered when comparing isotopic exchange effects for different elements, as well as for estimating water age using radioactive and cosmogenic isotopes.     

利用成像测井资料分析汶川地震断裂带科学钻探3号孔(WFSD-3)裂缝特征

WFSD-3孔是汶川地震断裂带科学钻探主要钻孔之一,全井段(终孔深度1 502.30 m)实施了连续取心(累计取心进尺1 548.44 m)和测井作业。采集到的成像测井资料包含丰富的原位地质特征信息,对裂缝、破碎带识别和构造应力场分析具有重要作用。利用该钻孔的电阻率成像测井和超声成像测井资料,结合岩心资料进行了裂缝特征分析。结果表明,WFSD-3孔岩层高角度斜交缝最为发育,其次是低角度斜交缝,垂直缝和水平缝极少;25～200 m和900～1 000 m深度范围内裂缝尤为发育;不同深度的裂缝倾向存在明显差异:410 m之上主要分布于260°～290°,410～730m集中于330°～360°,730～960 m主要分布于210°～240°,960～1 185 m与410 m之上基本一致,1 410～1 450 m与前述各深度段不同,集中于180°～200°;裂缝与破碎带、层理密切相关,宏观分布受构造控制。