Characterization of dipping fractures in a transverselyisotropic background

Although it is believed that natural fracture sets predominantly have near‐vertical orientation, oblique stresses and some other mechanisms may tilt fractures away from the vertical. Here, we examine an effective medium produced by a single system of obliquely dipping rotationally invariant fractures embedded in a transversely isotropic with a vertical symmetry axis (VTI) background rock. This model is monoclinic with a vertical symmetry plane that coincides with the dip plane of the fractures. Multicomponent seismic data acquired over such a medium possess several distinct features that make it possible to estimate the fracture orientation. For example, the vertically propagating fast shear wave (and the fast converted PS‐wave) is typically polarized in the direction of the fracture strike. The normal‐moveout (NMO) ellipses of horizontal reflection events are co‐orientated with the dip and strike directions of the fractures, which provides an independent estimate of the fracture azimuth. However, the polarization vector of the slow shear wave at vertical incidence does not lie in the horizontal plane – an unusual phenomenon that can be used to evaluate fracture dip. Also, for oblique fractures the shear‐wave splitting coefficient at vertical incidence becomes dependent on fracture infill (saturation). A complete medium‐characterization procedure includes estimating the fracture compliances and orientation (dip and azimuth), as well as the Thomsen parameters of the VTI background. We demonstrate that both the fracture and background parameters can be obtained from multicomponent wide‐azimuth data using the vertical velocities and NMO ellipses of PP‐waves and two split SS‐waves (or the traveltimes of PS‐waves) reflected from horizontal interfaces. Numerical tests corroborate the accuracy and stability of the inversion algorithm based on the exact expressions for the vertical and NMO velocities.     

Fluid‐dependent shear‐wave splitting in fractured media

Azimuthal anisotropy in rocks can result from the presence of one or more sets of partially aligned fractures with orientations determined by the stress history of the rock. A shear wave propagating in an azimuthally anisotropic medium splits into two components with different polarizations if the source polarization is not aligned with the principal axes of the medium. For vertical propagation of shear waves in a horizontally layered medium containing vertical fractures, the shear‐wave splitting depends on the shear compliance of the fractures, but is independent of their normal compliance. If the fractures are not perfectly vertical, the shear‐wave splitting also depends on the normal compliance of the fractures. The normal compliance of a fluid‐filled fracture decreases with increasing fluid bulk modulus. For dipping fractures, this results in a decrease in shear‐wave splitting and an increase in shear‐wave velocity with increasing fluid bulk modulus. The sensitivity of the shear‐wave splitting to fluid bulk modulus depends on the interconnectivity of the fracture network, the permeability of the background medium and on whether the fracture is fully or partially saturated.     

Sensitivity of time-lapse seismic to reservoir stress path

The change in reservoir pore pressure due to the production of hydrocarbons leads to anisotropic changes in the stress field acting on the reservoir. Reservoir stress path is defined as the ratio of the change in effective horizontal stress to the change in effective vertical stress from the initial reservoir conditions, and strongly influences the depletion‐induced compaction behaviour of the reservoir. Seismic velocities in sandstones vary with stress due to the presence of stress‐sensitive regions within the rock, such as grain boundaries, microcracks, fractures, etc. Since the response of any microcracks and grain boundaries to a change in stress depends on their orientation relative to the principal stress axes, elastic‐wave velocities are sensitive to reservoir stress path. The vertical P‐ and S‐wave velocities, the small‐offset P‐ and SV‐wave normal‐moveout (NMO) velocities, and the P‐wave amplitude‐versus‐offset (AVO) are sensitive to different combinations of vertical and horizontal stress. The relationships between these quantities and the change in stress can be calibrated using a repeat seismic, sonic log, checkshot or vertical seismic profile (VSP) at the location of a well at which the change in reservoir pressure has been measured. Alternatively, the variation of velocity with azimuth and distance from the borehole, obtained by dipole radial profiling, can be used. Having calibrated these relationships, the theory allows the reservoir stress path to be monitored using time‐lapse seismic by combining changes in the vertical P‐wave impedance, changes in the P‐wave NMO and AVO behaviour, and changes in the S‐wave impedance.     

Inversion of normal moveout for monoclinic media1

Multiple vertical fracture sets, possibly combined with horizontal fine layering, produce an equivalent medium of monoclinic symmetry with a horizontal symmetry plane. Although monoclinic models may be rather common for fractured formations, they have hardly been used in seismic methods of fracture detection due to the large number of independent elements in the stiffness tensor. Here, we show that multicomponent wide-azimuth reflection data (combined with known vertical velocity or reflector depth) or multi-azimuth walkaway VSP surveys provide enough information to invert for all but one anisotropic parameters of monoclinic media. In order to facilitate the inversion procedure, we introduce a Thomsen-style parametrization for monoclinic media that includes the vertical velocities of the P-wave and one of the split S-waves and a set of dimensionless anisotropic coefficients. Our notation, defined for the coordinate frame associated with the polarization directions of the vertically propagating shear waves, captures the combinations of the stiffnesses responsible for the normal-moveout (NMO) ellipses of all three pure modes. The first group of the anisotropic parameters contains seven coefficients (ε^(1,2), δ^(1,2,3) and γ^(1,2)) analogous to those defined by Tsvankin for the higher-symmetry orthorhombic model. The parameters ε^(1,2), δ^(1,2) and γ^(1,2) are primarily responsible for the pure-mode NMO velocities along the coordinate axes x₁ and x₂ (i.e. in the shear-wave polarization directions). The remaining coefficient δ⁽³⁾ is not constrained by conventional-spread reflection traveltimes in a horizontal monoclinic layer. The second parameter group consists of the newly introduced coefficients ζ^(1,2,3) which control the rotation of the P-, S₁- and S₂-wave NMO ellipses with respect to the horizontal coordinate axes. Misalignment of the P-wave NMO ellipse and shear-wave polarization directions was recently observed on field data by Pérez et al. Our parameter-estimation algorithm, based on NMO equations valid for any strength of the anisotropy, is designed to obtain anisotropic parameters of monoclinic media by inverting the vertical velocities and NMO ellipses of the P-, S₁- and S₂-waves. A Dix-type representation of the NMO velocity of mode-converted waves makes it possible to replace the pure shear modes in reflection surveys with the PS₁- and PS₂-waves. Numerical tests show that our method yields stable estimates of all relevant parameters for both a single layer and a horizontally stratified monoclinic medium.     

Seismic dispersion and attenuation in layered porous rocks with fractures of varying orientations

Wave‐induced fluid flow plays an important role in affecting the seismic dispersion and attenuation of fractured porous rocks. While numerous theoretical models have been proposed for the seismic dispersion and attenuation in fractured porous rocks, most of them neglect the wave‐induced fluid flow resulting from the background anisotropy (e.g. the interlayer fluid flow between different layers) that can be normal in real reservoirs. Here, according to the theories of poroelasticity, we present an approach to study the frequency‐dependent seismic properties of more realistic and complicated rocks, i.e. horizontally and periodically layered porous rock with horizontal and randomly orienting fractures, respectively, distributed in one of the two periodical layers. The approach accounts for the dual effects of the wave‐induced fluid flow between the fractures and the background pores and between different layers (the interlayer fluid flow). Because C₃₃ (i.e., the modulus of the normally incident P‐wave) is directly related to the P‐wave velocity widely measured in the seismic exploration, and its comprehensive dispersion and attenuation are found to be most significant, we study mainly the effects of fracture properties and the stiffness contrast between the different layers on the seismic dispersion and attenuation of C₃₃. The results show that the increasing stiffness contrast enhances the interlayer fluid flow of the layered porous rocks with both horizontal and randomly orienting fractures and weakens the wave‐induced fluid flow between the fractures and the background pores, especially for the layered porous rock with horizontal fractures. The modelling results also demonstrate that for the considered rock construction, the increasing fracture density reduces the interlayer fluid flow while improves the dispersion and attenuation in the fracture‐relevant frequency band. Increasing fracture aspect ratio is found to reduce the dispersion and attenuation in the fracture‐relevant frequency band only, especially for the layered porous rock with horizontal fractures.     

Seismic attenuation and velocity dispersion in fractured rocks: The role played by fracture contact areas

The presence of fractures in fluid‐saturated porous rocks is usually associated with strong seismic P‐wave attenuation and velocity dispersion. This energy dissipation can be caused by oscillatory wave‐induced fluid pressure diffusion between the fractures and the host rock, an intrinsic attenuation mechanism generally referred to as wave‐induced fluid flow. Geological observations suggest that fracture surfaces are highly irregular at the millimetre and sub‐millimetre scale, which finds its expression in geometrical and mechanical complexities of the contact area between the fracture faces. It is well known that contact areas strongly affect the overall mechanical fracture properties. However, existing models for seismic attenuation and velocity dispersion in fractured rocks neglect this complexity. In this work, we explore the effects of fracture contact areas on seismic P‐wave attenuation and velocity dispersion using oscillatory relaxation simulations based on quasi‐static poroelastic equations. We verify that the geometrical and mechanical details of fracture contact areas have a strong impact on seismic signatures. In addition, our numerical approach allows us to quantify the vertical solid displacement jump across fractures, the key quantity in the linear slip theory. We find that the displacement jump is strongly affected by the geometrical details of the fracture contact area and, due to the oscillatory fluid pressure diffusion process, is complex‐valued and frequency‐dependent. By using laboratory measurements of stress‐induced changes in the fracture contact area, we relate seismic attenuation and dispersion to the effective stress. The corresponding results do indeed indicate that seismic attenuation and phase velocity may constitute useful attributes to constrain the effective stress. Alternatively, knowledge of the effective stress may help to identify the regions in which wave induced fluid flow is expected to be the dominant attenuation mechanism.     

不同尺度裂缝对弹性波速度和各向异性影响的实验研究

裂缝广泛分布于地球介质中并且具有多尺度的特点,裂缝尺度对于油气勘探和开发有着重要的意义.本文制作了一组含不同长度裂缝的人工岩样,其中三块含裂缝岩样中的裂缝直径分别为2 mm、3 mm和4 mm,裂缝的厚度都约为0.06 mm,裂缝密度大致相同(分别为4.8%、4.86%和4.86%).在岩样含水的条件下测试不同方向上的纵横波速度,实验结果表明,虽然三块裂缝岩样中的裂缝密度大致相同,但是含不同直径裂缝岩样的纵横波速度存在差异.在各个方向上,含数量众多的小尺度裂缝的岩样中纵横波速度都明显低于含少量的大尺度裂缝的岩样中纵横波速度.尤其是对纵波速度和SV波速度,在不同尺度裂缝岩样中的差异更明显.在含数量多的小尺度裂缝的岩样中纵波各向异性和横波各向异性最高,而含少量的大尺度的裂缝的岩样中的纵波各向异性和横波各向异性较低.实验测量结果与Hudson理论模型预测结果进行了对比分析,结果发现Hudson理论考虑到了裂缝尺度对纵波速度和纵波各向异性的影响,但是忽略了其对横波速度和横波各向异性的影响.     

The effect of orthorhombic anisotropy and its implication for oil recovery and reservoir exploitation

Results of an experimental study of shear‐ and compressional‐wave propagation in an orthorhombically anisotropic medium are presented. The experiments were performed on a physical model consisting of two sets of fractures. The first set consisted of orientated rubber inclusions simulating weak material‐filled cracks. The second set consisted of a system of closely spaced parallel fractures simulated by thin plates of epoxy resin, superimposed orthogonally on the first set. Three cases of fracture orientations within the model were identified and studied. Case 1 is analogous to a jointed fracture reservoir with one vertical set of fluid‐filled cracks or fractures and one non‐filled horizontal set. This case is referred to as JFV. Case 2 is analogous to a double fracture reservoir with one horizontal set of fluid‐filled fractures or cracks and one non‐filled vertical set. This is referred to as DFH. Case 3 is analogous to a double fracture reservoir with two vertical sets of fractures or cracks, with only one fluid‐filled. Case 3 is referred to as DFV. A pulse transmission method was performed on all three modelled cases along the three principal axes. A directional variation in the compressional‐ and shear‐wave velocities, as well as distinct shear‐wave splitting, was observed. The elastic constants for each case were determined and differences between them were noted and compared with the controlled results of both layered (transverse isotropy, TI) and vertically fractured (azimuthally anisotropic models, VF) media. The differences in elastic moduli and velocities indicate the potential of recognizing the different fracture orientations and suggest an approach to designing a method of drilling to further enhance oil recovery and reservoir exploitation.     

Pore fluid viscosity effects on P‐ and S‐wave anisotropy in synthetic silica‐cemented sandstone with aligned fractures

Ultrasonic (500 kHz) P‐ and S‐wave velocity and attenuation anisotropy were measured in the laboratory on synthetic, octagonal‐shaped, silica‐cemented sandstone samples with aligned penny‐shaped voids as a function of pore fluid viscosity. One control (blank) sample was manufactured without fractures, another sample with a known fracture density (measured from X‐ray CT images). Velocity and attenuation were measured in four directions relative to the bedding fabric (introduced during packing of successive layers of sand grains during sample construction) and the coincident penny‐shaped voids (fractures). Both samples were measured when saturated with air, water (viscosity 1 cP) and glycerin (100 cP) to reveal poro‐visco‐elastic effects on velocity and attenuation, and their anisotropy. The blank sample was used to estimate the background anisotropy of the host rock in the fractured sample; the bedding fabric was found to show transverse isotropy with shear wave splitting (SWS) of 1.45 ± 1.18% (i.e. for S‐wave propagation along the bedding planes). In the fractured rock, maximum velocity and minimum attenuation of P‐waves was seen at 90° to the fracture normal. After correction for the background anisotropy, the fractured sample velocity anisotropy was expressed in terms of Thomsen's weak anisotropy parameters ε, γ & δ. A theory of frequency‐dependent seismic anisotropy in porous, fractured, media was able to predict the observed effect of viscosity and bulk modulus on ε and δ in water‐ and glycerin‐saturated samples, and the higher ε and δ values in air‐saturated samples. Theoretical predictions of fluid independent γ are also in agreement with the laboratory observations. We also observed the predicted polarisation cross‐over in shear‐wave splitting for wave propagation at 45° to the fracture normal as fluid viscosity and bulk modulus increases.     

OnP F and other phases in the early part of seismogram

Summary It has been found that when seismic energy propagates along the surface of the homogeneous crust beside usual Rayleigh waves, it produces certain instability in layers through which it propagates. In the light of this instability, a type of motion corresponding to longitudinal wave will be prominent in horizontal component compared to the vertical component; while transverse wave will be prominent in the vertical component but weak in the horizontal component, a contradiction with the existing knowledge. This has been identified withP _F phase. On taking the medium of propagation as slightly heterogeneous which allows existence of low velocity layer, a few larger number of such instabilities have been found. Velocity equation for Rayleigh waves for such media reveals existence of different velocities corresponding to vertical and horizontal components. Table for these velocities has been furnished.     

裂缝诱导双相HTI介质地震波场错格伪谱法模拟与波场特征分析

裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.     

Velocity dispersion in fractured rocks in a wide frequency range

Experimental measurements of fracture-induced seismic waves velocity variations at frequencies ~ 1 kHz, ~ 40 kHz and ~ 1 MHz were performed directly in the field at the rocky outcrop and in the laboratory on specific rock samples collected from the outcrops. The peridotite–lherzolite outcrop appeared macroscopically uniform and contained three systems of visible parallel sub-vertical fractures. This rock has substantial bulk density and higher than average value of seismic wave velocity. The presence of fracture systems gives rise to its velocity anisotropy. The seismic waves passing through the rock fractures are subject to velocity dispersion and frequency dependent attenuation. Our data, obtained from field and laboratory measurements, were compared with theoretical model predictions. In this model we successfully used displacement discontinuity approach. For the velocity dispersion evaluation we used multi-frequency measurements. The a priori observation of orientations and densities of fracture sets allowed evaluation of their stiffness. Our approach revealed that the first arrivals of seismic waves can be used for evaluation of P-wave group velocities, the specific case, in which we expect anomalous velocity dispersion. Our observations contribute to the issue of up-scaling of well-log derived velocities in fractured rock to the scale of standard seismic exploration frequencies.     

Field measurements of normal and shear fracture compliance

Field estimates of the normal and shear compliance of water saturated fractures were obtained from a seismic experiment carried out on a wave‐cut platform of Upper Caithness Flagstone on the North coast of Scotland. The rocks are cut by two orthogonal sets of vertical fractures. Vertical geophones were glued to the rock surface and seismic waves were generated by striking the surface with a sledge hammer. First arrival traveltimes were obtained for source‐receiver distances up to 30 m. Differences of compressional velocity with direction across the platform were interpreted as resulting from the compliance of the fractures intersected by the transmitted waves. The average normal compliance for fractures at this site is 4.0 × 10^?12± 1.0 × 10^?12 m/Pa. There is much less certainty about the precise value of shear compliance. However, the normal/shear compliance ratio is determined to be ≤0.1.     

Seismic anisotropy in a hydrocarbon field estimated from microseismic data

The study of seismic anisotropy in exploration seismology is gaining interest as it provides valuable information about reservoir properties and stress directions. In this study we estimate anisotropy in a petroleum field in Oman using observations of shear‐wave splitting from microseismic data. The data set was recorded by arrays of borehole geophones deployed in five wells. We analyse nearly 3400 microearthquakes, yielding around 8500 shear‐wave splitting measurements. Stringent quality control reduces the number of reliable measurements to 325. Shear‐wave splitting modelling in a range of rock models is then used to guide the interpretation. The difference between the fast and slow shear‐wave velocities along the raypath in the field ranges between 0–10% and it is controlled both by lithology and proximity to the NE‐SW trending graben fault system that cuts the field formations. The anisotropy is interpreted in terms of aligned fractures or cracks superimposed on an intrinsic vertical transversely isotropic (VTI) rock fabric. The highest magnitudes of anisotropy are within the highly fractured uppermost unit of the Natih carbonate reservoir. Anisotropy decreases with depth, with the lowest magnitudes found in the deep part of the Natih carbonate formation. Moderate amounts of anisotropy are found in the shale cap rock. Anisotropy also varies laterally with the highest anisotropy occurring either side of the south‐eastern graben fault. The predominant fracture strikes, inferred from the fast shear‐wave polarizations, are consistent with the trends of the main faults (NE‐SW and NW‐SE). The majority of observations indicate subvertical fracture dip (>70° ). Cumulatively, these observations show how studies of shear‐wave splitting using microseismic data can be used to characterize fractures, important information for the exploitation of many reservoirs.     

利用含可控裂缝人工岩样研究裂缝密度对各向异性的影响

利用新方法制作出含可控裂缝的双孔隙人工砂岩物理模型,具有与天然岩石更为接近的矿物成分、孔隙结构和胶结方式,其中裂缝密度、裂缝尺寸和裂缝张开度等裂缝参数可以控制以得到实验所需要的裂缝参数,岩样具有真实的孔隙和裂缝空间并可以在不同饱和流体状态下研究流体性质对于裂缝介质性质的影响.本次实验制作出一组具有不同裂缝密度的含裂缝人工岩样,对岩样利用SEM扫描电镜分析可以看到真实的孔隙结构和符合我们要求的裂缝参数,岩样被加工成八面棱柱以测量不同方向上弹性波传播的速度,用0.5 MHz的换能器使用透射法测量在饱和空气和饱和水条件下各个样品不同方向上的纵横波速度,并得出纵横波速度、横波分裂系数和纵横波各向异性强度受裂缝密度和饱和流体的影响.研究发现流体对于纵波速度和纵波各向异性强度的影响较强,而横波速度、横波分裂系数和横波各向异性强度受饱和流体的影响不大,但是对裂缝密度的变化更敏感.     

Elastic wave velocities in a granitic geothermal reservoir

Geothermal resources have potential for providing cost-effective and sustainable energy. Monitoring of production-induced changes in geothermal reservoirs using seismic waves requires understanding of the elastic properties of the rock and how they change due to injection of fluids and opening and closing of natural and hydraulic fractures. P- and S-wave velocities measured in a granitic geothermal reservoir using sonic logging are systematically lower than those predicted using the composition of the rock. Cracks may occur in granitic rocks from tectonic stresses and from the thermal expansion mismatch between differently oriented anisotropic crystals. An isotropic orientation distribution of microcracks causes a significant reduction in both the P- and S-velocities, consistent with the observed sonic P- and S-velocities. Vertical fractures cause a difference in the velocity of vertically propagating shear waves polarized parallel and perpendicular to the fractures. An assumption that the lower measured velocities are caused by the presence of vertical fractures is inconsistent with the sonic data. This is because vertical fractures cause a decrease in slow S-wave velocity that greatly exceeds the decrease in P-wave velocity, in contrast to the observed data. The growth of vertical fractures in the geothermal reservoir may be monitored using the difference in velocity of the fast and slow shear waves, while the change in P-velocity in a crossplot of measured P- and slow S-velocities is useful for estimating the ratio of the normal-to-shear compliance of the fractures.     

Analysis of different parameterisations of waveform inversion of compressional body waves in an elastic transverse isotropic Earth with a vertical axis of symmetry

In a multi‐parameter waveform inversion, the choice of the parameterisation influences the results and their interpretations because leakages and the tradeoff between parameters can cause artefacts. We review the parameterisation selection when the inversion focuses on the recovery of the intermediate‐to‐long wavenumbers of the compressional velocities from the compressional body (P) waves. Assuming a transverse isotropic medium with a vertical axis of symmetry and weak anisotropy, analytical formulas for the radiation patterns are developed to quantify the tradeoff between the shear velocity and the anisotropic parameters and the effects of setting to zero the shear velocity in the acoustic approach. Because, in an anisotropic medium, the radiation patterns depend on the angle of the incident wave with respect to the vertical axis, two particular patterns are discussed: a transmission pattern when the ingoing and outgoing slowness vectors are parallel and a reflection pattern when the ingoing and outgoing slowness vectors satisfy Snell's law. When the inversion aims at recovering the long‐to‐intermediate wavenumbers of the compressional velocities from the P‐waves, we propose to base the parameterisation choice on the transmission patterns. Since the P‐wave events in surface seismic data do not constrain the background (smooth) vertical velocity due to the velocity/depth ambiguity, the preferred parameterisation contains a parameter that has a transmission pattern concentrated along the vertical axis. This parameter can be fixed during the inversion which reduces the size of the model space. The review of several parameterisations shows that the vertical velocity, the Thomsen parameter δ, or the Thomsen parameter ε have a transmission pattern along the vertical axis depending on the parameterisation choice. The review of the reflection patterns of those selected parameterisations should be done in the elastic context. Indeed, when reflection data are also inverted, there are potential leakages of the shear parameter at intermediate angles when we carry out acoustic inversion.     

Estimation of reservoir fracturing from marine VSP using local shear-wave conversion

A marine VSP is designed to estimate the orientation and density of fracturing within a gas-producing dolomite layer in the southern North Sea. The overburden anisotropy is firstly estimated by analysing shear waves converted at or just below the sea-bed, from airgun sources at four fixed offset azimuths. Full-wave modelling helps confirm that the background has no more than 3% vertical birefringence, originating from TIH anisotropy with a symmetry axis orientated perpendicular to the maximum horizontal compressive stress of NW–SE. This finding concurs with current hypotheses regarding the background rock matrix in the upper crust. More detailed anisotropy estimates reveal two thin zones with possible polarization reversals and a stronger anisotropy. The seismic anisotropy of the dolomite is then determined from the behaviour of locally converted shear waves, providing a direct link with the physical properties of its fractures. It is possible to utilize this phenomenon due to the large seismic velocity contrast between the dolomite and the surrounding evaporites. Two walkaway VSPs at different azimuths, recorded on three-component receivers placed inside the target zone, provide the appropriate acquisition design to monitor this behaviour. Anisotropy in the dolomite generates a transverse component energy which scales in proportion to the degree of anisotropy. The relative amplitudes, for this component, between the different walkaway azimuths relate principally to the orientation of the anisotropy. Full-wave modelling confirms that a 50% vertical birefringence from TIH anisotropy with a similar orientation to the overburden is required to simulate the field observations. This amount of anisotropy is not entirely unexpected for a fine-grained brittle dolomite with a potentially high fracture intensity, particularly if the fractures contain fluid which renders them compliant to the shear-wave motion.