Seismic characterization of reservoirs containing multiple fracture sets

Natural fractures in reservoirs play an important role in determining fluid flow during production and knowledge of the orientation and density of fractures is required to optimize production. Variations in reflection amplitude versus offset (AVO) are sensitive to the presence of fractures but current models used to invert the seismic response often make simplified assumptions that prevent fractured reservoirs from being characterized correctly. For example, many models assume a single set of perfectly aligned fractures, whereas most reservoirs contain several fracture sets with variable orientation within a given fracture set. In addition, many authors only consider the azimuthal variation in the small offset amplitude versus offset and azimuth response (the variation in AVO gradient with azimuth), while the effect of fractures on amplitude versus offset and azimuth increases with increasing offset. In this paper, the variation in the reflection coefficient of seismic P -waves as a function of azimuth and offset due to the presence of multiple sets of fractures with variable orientation within any fracture set is used to determine the components of a second-rank fracture compliance tensor  α _ij   . The variation in the trace of this tensor as a function of position in the reservoir can be used to estimate the variation in fracture density with position in the reservoir and may be used to choose the location of infill wells in the field. The principal axes of  α _ij   reveal the most compliant direction within the reservoir and may be used to optimize the trajectory of deviated wells. The determination of the principal axes of  α _ij   requires wide azimuth acquisition and the use of the small-offset amplitude versus offset and azimuth (the azimuthal variation of the AVO gradient) may give misleading results.     

Azimuthal variation in AVO response for fractured gas sands

Natural fractures in reservoirs play an important role in determining fluid flow during production, and hence the density and orientation of fractures is of great interest. In the presence of aligned vertical fractures, the reflection amplitude at finite offset varies with azimuth. The effect of natural fractures on the azimuthal AVO response from a gas-sandstone reservoir encased within shale is investigated. A simple expression for the difference in P-wave reflection coefficient from the top of the reservoir parallel and perpendicular to the strike of the fractures is obtained in terms of the normal and tangential compliances, Z_N and Z_T, of the fractures. This expression is valid for small anisotropy and material contrasts and is compared with the results of numerical modelling. For a given value of Z_T, the azimuthal variation in reflection coefficient at moderate offsets is found to increase with decreasing Z_N/Z_T. For gas-filled open fractures Z_N/Z_T ≈ 1, but a lower ratio of Z_N/Z_T may result from the presence of cement or clay within the fractures, or from the presence of a fluid with non-zero bulk modulus. For Z_N/Z_T = 1 and moderate offsets, the variation with offset of the reflection coefficient from the top of the fractured unit is dominated by the contrast in Poisson's ratio between the gas sand and the overlying shale, the effect of fractures only becoming noticeable as the critical angle for the unfractured sandstone is approached. However, for reflections from the base of the fractured unit, the variation in reflection amplitude with azimuth is much greater at conventional seismic offsets than for the reflection from the top. Azimuthal variations in the strength of the reflection from the top of the reservoir depend only on the variation in reflection coefficient, whereas the raypath is also a function of azimuth for reflections from the base of the fractured unit, leading to stronger, more visible, variations of AVO with azimuth. It follows that an azimuthal variation in AVO due to fractures in the overburden may be misinterpreted as due to the presence of aligned fractures in the reservoir.     

The effect of mineralization on the ratio of normal to tangential compliance of fractures

The effect of a fracture on the propagation of seismic waves can be represented in terms of the normal compliance B_N and tangential compliance B_T of the fracture. If   B_N / B_T = 1  for all fractures, the effective elastic stiffness tensor of an isotropic background containing an arbitrary orientation distribution of fractures is orthotropic (i.e., has three orthogonal planes of mirror symmetry) in the long-wave limit. However, deviations from orthotropy may occur if   B_N / B_T   differs significantly from unity and this can cause the principal axes of the P -wave NMO ellipse and of the variation in the PP -reflection amplitude as a function of azimuth, to deviate from the fast and slow polarization direction of a vertically propagating S -wave. Simple models of a fracture in terms of a planar distribution of cracks suggest that   B_N / B_T ≈ 1  for dry fractures. However, naturally occurring fractures often exhibit mineralization in the form of bridges between opposing faces of the fracture. The presence of such bridges leads to significant departures of   B_N / B_T   from unity.     

Chrome-rich clinopyroxene in orthopyroxenite from Maowu, Dabie Mountains, central China: A second record and its implications for petrogenesis

A second occurrence of chrome-rich clinopyroxene has been discovered as inclusions in orthopyroxene in orthopyroxenite from Maowu, the Dabie Mountains, Central China. The average formula for chrome-rich clinopyroxene can be expressed as (Na_0.39Ca_0.54)_0.93(Mg_0.57Fe²⁺_0.06Fe³⁺_0.01Cr_0.24Al_0.15)_1.03Si_2.02O₆, with a maximum amount of kosmochlor component of 28.52 mol%. The unit cell parameters obtained from a single-crystal are a  = 9.614 Å, b  = 8.800 Å, c =  5.240 Å, β = 106.59°, space group C2 / c . The indices of refraction are α = 1.697, β = 1.704, γ = 1.726. Chrome-rich clinopyroxene, which coexists with chromite, chromian rutile and chromian pyrope, crystallized at a temperature of 1025 °C and very high pressure, and therefore represents a mantle relic. Together with the appearance of low-pressure inclusion mineral assemblage and the estimation of physical–chemical conditions for matrix minerals, the Maowu eclogite–ultramafic complex is considered to be formed during ultrahigh-pressure metamorphism from the mantle-derived protolith.     

煤层气储层物性预测的AVO技术对地震纵波资料品质要求的探讨

煤层的含气性及渗透性可以利用纵波AVO技术及基于各向异性理论的方位AVO技术进行研究.煤层气属于吸附气,储层厚度较薄,地震异常响应较弱,对地震资料的品质及其携带的信息量提出了更高的要求.本文对AVO与方位AVO反演理论的近似条件及模型试算结果进行了分析,提出了对地震纵波数据信噪比、入射角、偏移距分布及方位角分布等参数的要求:1.观测系统设计及野外采集时应保证地震波的有效入射角在0.～30.或更大范围内,偏移距应大中小均匀分布,避免过于集中或缺失.2.地震数据应有足够宽分布的方位角,煤层气地震勘探中应选择宽方位观测系统.3、地震资料信噪比应足够高.     

Geochemistry and origin of Archean volcanic rocks from the Upper Keewatin assemblage (ca 2.7 Ga), Lake of the Woods Greenstone Belt, Western Wabigoon Subprovince, Superior Province, Canada

Abstract   Volcanic rocks from the Upper Keewatin assemblage ( ca 2720 Ma) were geochemically classified into five groups; komatiites, tholeiitic rocks having near-flat primitive mantle-normalized abundance patterns, Nb-enriched basalts and andesites (NEBA) plus normal calc-alkaline (NCA) rocks, adakites and shoshonites. The adakites having [La/Yb]_N >30 and <30 were probably derived from felsic magmas formed by partial melting of a subducted slab at relatively greater and smaller depths, respectively. Ascending adakite magmas, by interaction with the overlying mantle wedge, decreased in Al₂O₃ / Y ratio and selectively lost high-field strength elements, thereby forming mantle sources for both NEBA + NCA and shoshonite magmas. Under the influence of a mantle plume, the source of komatiites, the NEBA + NCA magmas were generated from that part of the mantle wedge metasomatized by adakite magmas having [La / Yb]_N <30, and tholeiitic magmas from unmetasomatized part of the same mantle wedge. Magmas of both adakites having [La / Yb]_N >30 and shoshonites were generated in a normal Archean Arc system setting.     

Hydraulic Conductivity Measurements in the Unsaturated Zone Using Improved Well Analyses

The flow of ponded water into and through the unsaturated zone depends on both the saturated and unsaturated components of the hydraulic conductivity. Recent studies indicate that the ratio of the saturated (K_fs) to the unsaturated (φ_m) components (K_fs/φ_m=α*) of flow lies within prescribed bounds for most field soils, i.e., 1m⁻¹≤α*≤ 100 m⁻¹. In addition, the fact that the calculation of K_fs and φ_m is not strongly dependent on the choice of α*, suggests that a site estimation of α* leads to reasonable "best estimates" of K_fs and φ_m when using the constant head well permeameter technique. As a consequence, measurement of the steady flow rate using only one ponded head may be all that is necessary for many practical applications. Multiple head measurements or independent measurements of α* or φ_m can be used, however, to give more accurate estimates of K_fs if required.     

HTI介质中的反射纵波方位属性

利用横波响应进行裂缝性各向异性介质的检测在实际应用中取得了很好的效果,但技术复杂、成本较高使该方法的广泛使用受到限制,而纵波资料采集和处理技术的精细有效保持了纵波的各种属性,这为直接利用纵波资料进行裂缝检测创造了条件.地下垂直定向裂缝通常用HTI介质模型来描述,为此,本文利用射线追踪和反射率法计算了层状各向同性介质背景下的HTI介质顶、底界面的反射纵波旅行时和反射系数,并分析了这些属性随观测方位的变化规律.研究表明,HTI介质底界面反射纵波旅行时和HTI介质顶界面反射系数表现出了明显的方位各向异性;旅行时、振幅和AVO梯度属性均在0°观测方位和90°观测方位上存在最大差异,可以用多种属性联合来精确判定裂缝的发育方向.     

Electromagnetic profiling interpretation across vertical faults and dikes

When a vertical fault or dike scatters a normally incident TE-mode plane wave, the horizontal components of the electric and the magnetic fields vary along the direction perpendicular to the strike. This scattering is also responsible for the formation of a vertical component of the magnetic field, which also varies along the direction perpendicular to the strike. Analysis of the lateral variations of the field components permits the identification of the type of geological structure, either fault or dike, as well as an estimate of the position of the geological contacts and the electrical conductivity of each medium. Previous exact analytical solutions have shown that two Fourier cosine integrals can express each field component. A series of functions specified for each model represent the kernel of each integral. From the field components obtained from the first non-zero term of each series, we calculated the following functions: the dip angle and the ellipticity of the vertical polarization ellipse in a plane perpendicular to the strike of the structure, and the azimuth and the ellipticity of the horizontal ellipse. We established master-curves of these functions for the interpretation of vertical faults and dikes for the polarization ellipsoid, for instance VLF or audio-frequency methods. This representation has as variable the induction number     , and as parameters σ ₂/ σ ₁ and     , where a denotes the half-thickness of the dike. The interpretation procedure using curve fitting is possible because the induction number is represented on a logarithmic scale. This procedure represents an important step before automatic interpretation is carried out. Two case histories using field data illustrate the effectiveness of the procedure and the type of quantitative interpretation obtained from it.     

Numerical analysis of seismic wave propagation in fluid-saturated porous multifractured media

Elastic wave propagation and attenuation in porous rock layers with oriented sets of fractures, especially in carbonate reservoirs, are anisotropic owing to fracture sealing, fracture size, fracture density, filling fluid, and fracture strike orientation. To address this problem, we adopt the Chapman effective medium model and carry out numerical experiments to assess the variation in P-wave velocity and attenuation, and the shear-wave splitting anisotropy with the frequency and azimuth of the incident wave. The results suggest that velocity, attenuation, and anisotropy vary as function of azimuth and frequency. The azimuths of the minimum attenuation and maximum P-wave velocity are nearly coincident with the average strike of the two sets of open fractures. P-wave velocity is greater in sealed fractures than open fractures, whereas the attenuation of energy and anisotropy is stronger in open fractures than sealed fractures. For fractures of different sizes, the maximum velocity together with the minimum attenuation correspond to the average orientation of the fracture sets. Small fractures affect the wave propagation less. Azimuth-dependent anisotropy is low and varies more than the other attributes. Fracture density strongly affects the P-wave velocity, attenuation, and shear-wave anisotropy. The attenuation is more sensitive to the variation of fracture size than that of velocity and anisotropy. In the seismic frequency band, the effect of oil and gas saturation on attenuation is very different from that for brine saturation and varies weakly over azimuth. It is demonstrated that for two sets of fractures with the same density, the fast shear-wave polarization angle is almost linearly related with the orientation of one of the fracture sets.     

基于二维三分量伪谱法模拟数据的EDA介质中横波分裂研究

为了检测定向裂隙介质中横波分裂的方位属性特征,分析地震属性随裂隙密度和方位变化,采用人工吸收边界和反周期扩展边界,用伪谱法获得不同裂隙密度和不同方位地质模型三分量地面记录;应用时频分析和剪切波偏振分析研究由于裂隙方位和密度引起的横波分裂.结果显示,裂隙密度和方位决定着横波分裂的时差和偏振.快慢横波的延迟时间随裂隙密度增大而增加,不同方位相同裂隙密度的横波分裂时差有微小的变化.在45°方位检测时间延迟时间最大.通过时频分析,可以看到不同方位的瞬时主频有显著的变化,在横波分裂处瞬时主频有明显变化.因此,瞬时主频和快横波的偏振以及延迟时间可以作为裂隙方位和密度的指示.     

Relationship between P- and S-wave velocities and geological properties of near-surface sediments of the continental slope of the Barents Sea

Seismic velocities ( V _p and V _s) of compressional (P-) and shear (S-) waves are important parameters for the characterization of marine sediments with respect to their sedimentological and geotechnical properties. P- and S-wave velocity data of near-surface marine sediments (upper 9 m) of the continental slope of the Barents Sea are analysed and correlated to sedimentological and geotechnical properties. The results show that the S-wave velocity is much more sensitive to changes in lithology and mechanical properties than the P-wave velocity, which is characterized by a narrow range of values. The correlation coefficients between S-wave velocity and silt and clay content, wet bulk density, porosity, water content and shear strength are higher than 0.5 while the correlation coefficients of P-wave velocity and the same parameters are always lower than 0.4. Although the relationship between V _s and clay content has been widely described, the data show that V _s is better correlated with silt content than with clay content for the sediments of the area investigated. However, they show different trends. While V _s increases with increasing clay content, it decreases with increasing silt content.     

Paleostress orientation from 3-D orientation distribution of microcracks in quartz from the Cretaceous granodiorite core samples drilled through the Nojima Fault, south-west Japan

Abstract The orientations of both healed extension microcracks and microcracks in quartz grains sealed mostly by carbonate minerals were measured from Cretaceous granodiorite core samples drilled along the Nojima Fault, southwest Japan. The preferred orientations of both healed and sealed microcracks consist of approximately three orthogonal sets, (components) A, B and C, in which A strikes NS–NW-SE and dips vertically, B strikes EW–NE-SW and dips vertically, and C is subhorizontal. Both the healed and sealed microcracks were possibly formed by hydraulic fracturing, and the successive release of tensile stress due to pore fluid overpressure in the principal stress directions could have caused this microcracking in mutually orthogonal directions. The quartz grains are also very moderately plastically deformed, which is indicated by the occurrence of kink bands and undulose extinction. The association of healed microcracks and kink bands in the quartz suggests that these microstructures formed under subgreenschist facies conditions (≈300°C) during hydrothermal activity that could have occurred immediately after the emplacement of granodiorite during the Late Cretaceous period. Based on both the preferred orientation of microcracks, and c-axis fabrics of kinked and unkinked grains (so called kink method), it is inferred that the σ₁-and σ₃-axis were oriented horizontally in NS–NW-SE and EW–NE-SW directions, respectively. The inferred paleostress field does not conform to the east–west-trending compression during the Quaternary period, but to the activation of EW–NE-SW-trending, left–lateral strike–slip faults during the Late Cretaceous period in southwest Japan.     

Origins of hydrocarbons in the Sagara oil field, central Japan

Abstract   We collected free-gas and in situ fluid samples up to a depth of 200.6 m from the Sagara oil field, central Japan (34°44'N, 138°15'E), during the Sagara Drilling Program (SDP) and measured the concentrations and stable carbon isotopic compositions of CH₄ and C₂H₆ in the samples. A combination of the CH₄/C₂H₆ ratios with the carbon isotope ratios of methane indicates that the hydrocarbon gases are predominantly of thermogenic origin at all depths. The isotope signature of hydrocarbon gases of δ¹³      < δ¹³     suggests that these gases in the Sagara oil field are not generated by polymerization, but by the decomposition of organic materials.     

Diamond-bearing and diamond-free metacarbonate rocks from Kumdy–Kol in the Kokchetav Massif, northern Kazakhstan

Abstract Dolomite marble from the Kumdy–Kol area of the Kokchetav Massif contains abundant microdiamond, mainly in garnet and a few in diopside. The mineral assemblage at peak metamorphic condition consists of dolomite + diopside + garnet (+ aragonite) ± diamond. Inclusions of very low MgCO₃ calcite and almost pure calcite occur in diopside and are interpreted as aragonite and/or aragonite + dolomite. Single-phase Mg–calcite in diopside with a very high MgCO₃ component (up to 21.7 mol%) was also found in diamond-free dolomitic marble, and is interpreted as a retrograde product from aragonite + dolomite to Mg–calcite. The dolomite stability constrains the maximum pressure (P) at < 7 GPa using previous experimental data, whereas the occurrence of diamond yields the minimum peak pressure–temperature (P–T) condition at 4.2 GPa and 980 °C at X co ₂ = 0.1. The highest MgCO₃ in Mg–calcite constrains the minimum P–T condition higher than 2.5 GPa and 800 °C for the exhumation stage. As these marbles were subjected to nearly identical P–T metamorphic conditions, the appearance of diamond in some carbonate rocks was explained by high X co ₂. A low X co ₂ condition refers to high oxidized conditions and diamond (and/or graphite) becomes unstable. Difference in X co ₂ for marble from the same area suggests local heterogeneity of fluid compositions during ultrahigh-pressure metamorphism.     

Model-based shear-wave velocity estimation versus empirical regressions

Modelling of AVO signatures for reservoir characterization requires V _S estimation from other available logs when shear-wave data are not available. We tested various models for predicting V _S from P-wave velocity, porosity and shale volume measured in well logs. Effective medium models which characterize the pore space in terms of ellipsoidal inclusions were compared with statistical V _P– V _S regressions. The inclusion models were calibrated by non-linear minimization of the difference between model-predicted velocities and actual measured velocities. The quality of the V _S prediction was quantified in terms of the rms error by comparison with shear-wave data in wells where both V _P and V _S were measured. The linear regressions were found to be more robust and the rms error in the prediction was comparable to effective medium model-based predictions.     

Moveout approximation for horizontal transversely isotropic and vertical transversely isotropic layered medium. Part I: 1D ray propagation‡

Anisotropy in subsurface geological models is primarily caused by two factors: sedimentation in shale/sand layers and fractures. The sedimentation factor is mainly modelled by vertical transverse isotropy (VTI), whereas the fractures are modelled by a horizontal transversely isotropic medium (HTI). In this paper we study hyperbolic and non‐hyperbolic normal reflection moveout for a package of HTI/VTI layers, considering arbitrary azimuthal orientation of the symmetry axis at each HTI layer. We consider a local 1D medium, whose properties change vertically, with flat interfaces between the layers. In this case, the horizontal slowness is preserved; thus, the azimuth of the phase velocity is the same for all layers of the package. In general, however, the azimuth of the ray velocity differs from the azimuth of the phase velocity. The ray azimuth depends on the layer properties and may be different for each layer. In this case, the use of the Dix equation requires projection of the moveout velocity of each layer on the phase plane. We derive an accurate equation for hyperbolic and high‐order terms of the normal moveout, relating the traveltime to the surface offset, or alternatively, to the subsurface reflection angle. We relate the azimuth of the surface offset to its magnitude (or to the reflection angle), considering short and long offsets. We compare the derived approximations with analytical ray tracing.     

Paleomagnetism and tectonics of Karaginsky Island, Bering Sea

Abstract Rocks from Karaginsky accretionary prism (Karaginsky Island, Bering Sea) yield both prefolding (close to original) and postfolding magnetic vectors. The prefolding vectors suggest that the Maastrichtian–Paleocene volcanic–terrigenous sequences of Karaginsky Island formed at approximately 40°N to 50°N ( n = 45, D _G = 325, I _G = 57, K _G = 6, α_95G = 8, F _G = 15.06, D _S = 332, I _S = 63, K _S = 20, α_95S = 4.5, F _S = 0.3297, F _cr = 2.64) and were not originally part of either Eurasia ( F = 19, Δ F = 6.5) or North America ( F = 17, Δ F = 4.4). The geologic blocks rotated insignificantly counterclockwise about the horizontal plane, suggesting that the structure of Karaginsky Island arose without major strike-slip motions. Analysis of secondary magnetizations (for example, n = 28, D _G = 311, I _G = − 50, K _G = 9, α_95G = 8.7, F _G = 2.44; D _S = 293, I _S = − 41, K _S = 5, α_95S = 11, F _S = 12.04, F _cr = 2.65) reveals that the development of this framework involved at least two stages of deformation. During the second stage the sequences must have been tilted to west-northwest and northwest directions at 45–65°. This agrees with the northwest vergence of the structure of Karaginsky Island.     

Geographical distribution of helium isotope ratios in northeastern Japan

In order to study the precise geographical distribution of helium isotope ratios in northeastern Japan and compare it with geophysical data, we collected 43 gas and water samples from hot and mineral springs in the region where the ratio had never been reported, and measured the ³He/⁴He and ⁴He/²⁰Ne ratios of these samples. It was found that the ³He/⁴He ratios show clear contrasts between the forearc and the back-arc regions in the Tohoku district in northeastern Japan. In the forearc region, the ratios are smaller than 1 R_A (1 R_A = 1.4 × 10⁻⁶; R_A means the ³He/⁴He ratio of the atmosphere). On the other hand, those along the volcanic front and in the back-arc region are apparently higher. Moreover, we found a variation in the ³He/⁴He ratios along the volcanic front. In Miyagi Prefecture (38–39°N), the ratios range from 2 to 5 R_A. On the other hand, the ratios are less than 1 R_A in and around the southern border between Iwate and Akita Prefectures (39–39.5°N). Comparing the distribution of helium isotope ratios to results of recent geophysical studies, we found that the features in geographical distribution of helium isotope ratios are similar to those of seismic low-velocity zone distributions and high Qp⁻¹ distributions in the uppermost mantle. These observations strongly suggest that the helium isotope ratios reflect the distribution of melts in the uppermost mantle and are a useful tool for investigating the origin, behavior, and distribution of deep fluids and melts.     

瞬时频域衰减属性及其在碳酸盐岩裂缝检测中的应用

基于AVO梯度、走时及速度的纵波方位角裂缝检测已有较多研究和应用.由于裂缝及所含流体的综合作用,碳酸岩裂缝性储层的地震响应具有高频衰减、频带低移及方位各向异性的特点,因此应用衰减属性进行裂缝预测具有重要价值.基于小波变换计算的瞬时频域振幅谱及以此为基础提取的衰减属性,由于在时域和频域都有较好的分辨率能力,能够准确地刻画出碳酸岩裂缝性储层的非均质性.结合三维地震纵波数据的分方位角处理和分析技术,可以由方位衰减属性估算裂缝相对密度.对Z区下古生界潜山碳酸岩储层进行的裂缝预测表明,由衰减梯度和85%能量比频率两种属性导出的裂缝发育密度具有很好的相关性,且与实际情况符合.