Fracture detection by Gaussian beam imaging of seismic data and image spectrum analysis

Localization of fractured areas is of primary interest in the study of oil and gas geology in carbonate environments. Hydrocarbon reservoirs in these environments are embedded within an impenetrable rock matrix but possess a rich system of various microheterogeneities, i.e., cavities, cracks, and fractures. Cavities accumulate oil, but its flow is governed by a system of fractures. A distinctive feature of wave propagation in such media is the excitation of the scattered/diffracted waves by the microheterogeneities. This scattering could be a reliable attribute for characterization of the fine structure of reservoirs, but it has extremely low energy and any standard data processing renders them practically invisible in comparison with images produced by specular reflections. Therefore, any attempts to use these waves for image congestion of microheterogeneities should first have a preliminary separation of the scattering and specular reflections. In this paper, the approach to performing this separation is based on the asymmetric summation. It is implemented by double focusing of Gaussian beams. To do this, the special weights are computed by propagating Gaussian beams from the target area towards the acquisition system separately for sources and receivers. The different mutual positioning of beams in each pair introduces a variety of selective images that are destined to represent some selected singular primitives of the target objects such as fractures, cavities, and edges. In this way, one can construct various wave images of a target reservoir, particularly in scattered/diffracted waves. Additional removal of remnants of specular reflections is done by means of spectral analysis of the scattered/diffracted waves' images to recognize and cancel extended lineaments. Numerical experiments with Sigsbee 2A synthetic seismic data and some typical structures of the Yurubcheno‐Tokhomskoye oil field in East Siberia are presented and discussed.     

Wide-angle linear forward modelling of synthetic seismograms

We address the issue of linearity and scale dependence in forward modelling of seismic data from well logs, for large ray parameters, wide angles or large offsets. We present a forward model, within the context of seismic‐to‐well matching, that is linearized in the elastic properties of the earth. This model preserves linearity at large ray parameters and can handle fine‐layering effects such as induced anisotropy. Starting from a low‐contrast small‐ray‐parameter model, we extend it to a large‐ray‐parameter model by fully linearizing the elastic‐property contrasts. Overall linearity of the forward model is extended by partitioning the compressional‐wave and shear‐wave velocity fields into two fundamental scales: a kinematic scale that governs wavefield propagation effects and a dynamic scale that governs wavefield scattering effects. This analysis reveals that the standard practice in forward modelling of strongly filtering the ratio of compressional‐wave velocity to shear‐wave velocity is well founded in the underlying physics. The partitioning of the velocity fields also leads naturally to forward modelling that accounts fully for stretch effects, to resolution of the angle‐of‐incidence versus ray‐parameter dichotomy in seismic‐amplitude analysis, and to full accounting for induced anisotropy and dispersion effects due to fine‐layering of isotropic media. With the onset of routine long‐offset acquisition and the compelling need to optimize asset management in order to maximize reserve recovery, this forward model recognizes the physics of seismic wave propagation and enables a more complete exploitation of amplitude information in pre‐critical seismic data.     

Modelling the wave phenomena in acoustic and elastic media with sharp variations of physical properties using the grid‐characteristic method

This paper introduces a novel method of modelling acoustic and elastic wave propagation in inhomogeneous media with sharp variations of physical properties based on the recently developed grid‐characteristic method which considers different types of waves generated in inhomogeneous linear‐elastic media (e.g., longitudinal, transverse, Stoneley, Rayleigh, scattered PP‐, SS‐waves, and converted PS‐ and SP‐waves). In the framework of this method, the problem of solving acoustic or elastic wave equations is reduced to the interpolation of the solutions, determined at earlier time, thus avoiding a direct solution of the large systems of linear equations required by the FD or FE methods. We apply the grid‐characteristic method to compare wave phenomena computed using the acoustic and elastic wave equations in geological medium containing a hydrocarbon reservoir or a fracture zone. The results of this study demonstrate that the developed algorithm can be used as an effective technique for modelling wave phenomena in the models containing hydrocarbon reservoir and/or the fracture zones, which are important targets of seismic exploration.     

Elastic full waveform inversion of multi‐component OBC seismic data

Elastic full waveform inversion of seismic reflection data represents a data‐driven form of analysis leading to quantification of sub‐surface parameters in depth. In previous studies attention has been given to P‐wave data recorded in the marine environment, using either acoustic or elastic inversion schemes. In this paper we exploit both P‐waves and mode‐converted S‐waves in the marine environment in the inversion for both P‐ and S‐wave velocities by using wide‐angle, multi‐component, ocean‐bottom cable seismic data. An elastic waveform inversion scheme operating in the time domain was used, allowing accurate modelling of the full wavefield, including the elastic amplitude variation with offset response of reflected arrivals and mode‐converted events. A series of one‐ and two‐dimensional synthetic examples are presented, demonstrating the ability to invert for and thereby to quantify both P‐ and S‐wave velocities for different velocity models. In particular, for more realistic low velocity models, including a typically soft seabed, an effective strategy for inversion is proposed to exploit both P‐ and mode‐converted PS‐waves. Whilst P‐wave events are exploited for inversion for P‐wave velocity, examples show the contribution of both P‐ and PS‐waves to the successful recovery of S‐wave velocity.     

Full‐wavefield migration: using surface and internal multiples in imaging

Multiple scattering is usually ignored in migration algorithms, although it is a genuine part of the physical reflection response. When properly included, multiples can add to the illumination of the subsurface, although their crosstalk effects are removed. Therefore, we introduce full‐wavefield migration. It includes all multiples and transmission effects in deriving an image via an inversion approach. Since it tries to minimize the misfit between modeled and observed data, it may be considered a full waveform inversion process. However, full‐wavefield migration involves a forward modelling process that uses the estimated seismic image (i.e., the reflectivities) to generate the modelled full wavefield response, whereas a smooth migration velocity model can be used to describe the propagation effects. This separation of modelling in terms of scattering and propagation is not easily achievable when finite‐difference or finite‐element modelling is used. By this separation, a more linear inversion problem is obtained. Moreover, during the forward modelling, the wavefields are computed separately in the incident and scattered directions, which allows the implementation of various imaging conditions, such as imaging reflectors from below, and avoids low‐frequency image artefacts, such as typically observed during reverse‐time migration. The full wavefield modelling process also has the flexibility to image directly the total data (i.e., primaries and multiples together) or the primaries and the multiples separately. Based on various numerical data examples for the 2D and 3D cases, the advantages of this methodology are demonstrated.     

Analytical expressions for the borehole seismic signatures excited by an explosive source buried in an elastic half‐space

Borehole guided waves that are excited by explosive sources outside of the borehole are important for interpreting borehole seismic surveys and for rock property inversion workflows. Borehole seismograms are typically modelled using numerical methods of wave propagation. In order to benchmark such numerical algorithms and partially to interpret the results of modelling, an analytical methodology is presented here to compute synthetic seismograms. The specific setup is a wavefield emanating from a monopole point source embedded within a homogeneous elastic medium that interacts with a fluid‐filled borehole and a free surface. The methodology assumes that the wavelength of the seismic signal is much larger than the borehole radius. In this paper, it is supposed that there is no poroelastic coupling between the formation and the borehole. The total wavefield solution consists of P, PP, and PS body waves; the surface Rayleigh wave; and the low‐frequency guided Stoneley wave (often referred as the tube wave) within the borehole. In its turn, the tube wave consists of the partial responses generated by the incident P‐wave and the reflected PP and PS body waves at the borehole mouth and by the Rayleigh wave, as well as the Stoneley wave eigenmode. The Mach tube wave, which is a conic tube wave, additionally appears in the Mach cone in a slow formation with the tube‐wave velocity greater than the shear one. The conditions of appearance of the Mach wave in a slow formation are formulated. It is shown that the amplitude of the Mach tube wave strongly depends on Poisson's ratio of the slow surrounding formation. The amplitude of the Mach tube wave exponentially decreases when the source depth grows for weakly compressible elastic media with Poisson's ratio close to 0.5 (i.e., saturated clays and saturated clay soils). Asymptotic expressions are also provided to compute the wavefield amplitudes for different combinations of source depth and source‐well offset. These expressions allow an approximate solution of the wavefield to be computed much faster (within several seconds) than directly computing the implicit integrals arising from the analytical formulation.     

Detecting underground cavities using microtremor data: physical modelling and field experiment

In this paper, we study the possibilities of the use of microtremor records in the detection and delineation of near‐surface underground cavities. Three‐dimensional physical modelling data showed that the averaging amplitude spectra of a large number of microtremor records makes it possible to evaluate the frequencies and amplitudes of compressional standing waves generated by microtremor in the space between the ground surface and underground cavities. We illustrate how these parameters can be used to estimate the shape of the underground cavity horizontal projection. If the compressional wave velocity in the enclosing rock is known, it is possible to evaluate the depth to the cavity top using the frequencies of the standing waves. The results of the field experiment confirmed the possibility of underground cavities detection using microtremor data.     

Anisotropic elastic full‐waveform inversion of walkaway vertical seismic profiling data from the Arabian Gulf

Borehole seismic addresses the need for high‐resolution images and elastic parameters of the subsurface. Full‐waveform inversion of vertical seismic profile data is a promising technology with the potential to recover quantitative information about elastic properties of the medium. Full‐waveform inversion has the capability to process the entire wavefield and to address the wave propagation effects contained in the borehole data—multi‐component measurements; anisotropic effects; compressional and shear waves; and transmitted, converted, and reflected waves and multiples. Full‐waveform inversion, therefore, has the potential to provide a more accurate result compared with conventional processing methods. We present a feasibility study with results of the application of high‐frequency (up to 60 Hz) anisotropic elastic full‐waveform inversion to a walkaway vertical seismic profile data from the Arabian Gulf. Full‐waveform inversion has reproduced the majority of the wave events and recovered a geologically plausible layered model with physically meaningful values of the medium.     

Frequency‐dependent anisotropy of porous rocks with aligned fractures

Naturally fractured reservoirs are becoming increasingly important for oil and gas exploration in many areas of the world. Because fractures may control the permeability of a reservoir, it is important to be able to find and characterize fractured zones. In fractured reservoirs, the wave‐induced fluid flow between pores and fractures can cause significant dispersion and attenuation of seismic waves. For waves propagating normal to the fractures, this effect has been quantified in earlier studies. Here we extend normal incidence results to oblique incidence using known expressions for the stiffness tensors in the low‐ and high‐frequency limits. This allows us to quantify frequency‐dependent anisotropy due to the wave‐induced flow between pores and fractures and gives a simple recipe for computing phase velocities and attenuation factors of quasi‐P and SV waves as functions of frequency and angle. These frequency and angle dependencies are concisely expressed through dimensionless velocity anisotropy and attenuation anisotropy parameters. It is found that, although at low frequencies, the medium is close to elliptical (which is to be expected as a dry medium containing a distribution of penny‐shaped cracks is known to be close to elliptical); at high frequencies, the coupling between P‐wave and SV‐wave results in anisotropy due to the non‐vanishing excess tangential compliance.     

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.     

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

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

地震多级散射波正演模拟方法

地震波在穿越地下散射体群时会产生多级散射波,分析其地震响应特征,可推断散射体的分布情况和性质。本文从二维标量波动方程出发,结合地震散射理论和波恩近似理论,推导了多级散射波方程。在此基础上,采用高阶有限差分法对双点散射体模型和复杂散射体模型进行数值模拟,分析了多级散射波的传播规律和波场特征,并通过抽取多级散射记录和各级散射记录的单道记录与参考单道记录的对比,验证了本文推导散射波方程的准确性。      

Calibration of cross-line components for sea-bed 4C acquisition systems

Faithful recording of the elastic wavefield at the sea‐bed is required for quantitative applications of 4C seismic. The accuracy of the recorded vectorial wavefield depends on factors that vary from deployment to deployment. This paper focuses on one such factor: the interaction of the acquisition system with the sea‐bed, which is referred to here as coupling. We show, using multi‐azimuth data recorded with a cable‐based sea‐bed acquisition system, whose sensor housing is cylindrically shaped and with the in‐line geophone fixed to the cable, that coupling depends on the propagation direction and wave type (P‐ or S‐waves) of the incident wavefield. We show that coupling is more critical for S‐waves than for P‐waves. Detection of inconsistent coupling using both P‐ and S‐waves is therefore mandatory. A data‐driven processing method to compensate for the frequency‐dependent coupling response of the cross‐line geophone is derived. Its application to field data verifies the effectiveness of the method.     

地震干涉测量法近地表散射波分离技术

针对山地地震勘探数据低信噪比问题,近地表散射波分离意义显得尤为突出,地震干涉测量法为此提供了一种技术手段.本文将地震干涉测量理论和散射理论结合起来,导出了近地表散射波地震干涉测量表达式,分为互相关型和褶积型表达式,它们由实际波场和背景波场干涉测量构成.根据近地表散射波分离理论,结合陆上地震勘探实际观测系统,采用褶积和反褶积混合型地震干涉测量配置,用实际地震资料展示了近地表散射波分离技术的应用效果.经过理论分析和砾石区实际资料试验,表明地震干涉测量不仅能分离测线上散射源产生的散射波,而且能分离部分侧面散射波.该技术的优点在于它适应于起伏地形和不均匀近地表结构,并且不需要起伏地形和近地表速度信息.为了从实际资料中消除近地表散射波,本文采用多道匹配滤波自适应减法,在砾石区见到较好效果.     

Low-frequency microtremor anomalies at an oil and gas field in Voitsdorf, Austria

Results of a passive microtremor survey at an oil and gas field in Voitsdorf, Austria, are presented. The survey consists in six parallel profiles approximately 9 km long over two hydrocarbon reservoirs. For each profile the seismic wavefield was recorded synchronously at 11 in-line stations. The measurements were conducted with broadband seismometers and lasted, for each profile, at least 12 hours overnight. Data interpretation is based on a comprehensive data set and on the analysis of four different spectral attributes. These attributes quantify the characteristic features of the wavefield's Fourier spectra in the low-frequency range (<10Hz). One attribute quantifies the spectral energy in the vertical wavefield component, another attribute quantifies the maxima in vertical-to-horizontal spectral ratios and two attributes describe the frequency shifts of peaks within the spectra of vertical and horizontal wavefield components. Due to temporal variations of the signals we combine the long-term measurements (several hours of continuous records) of multiple profiles. This procedure considerably enhances the consistency of each spectral attribute and makes them suitable to quantify lateral variations of the wavefield. The results show that using a combination of several attributes significantly increases the reliability of detecting anomalies in the microtremor wavefield that are presumably caused by hydrocarbon reservoirs. A numerical study of two-dimensional seismic wave propagation is applied to investigate the peak frequency shift attributes. The results of the study indicate that the attributes may contain information on the depth of hydrocarbon reservoirs, assuming that the reservoir acts as a (secondary) source of low-frequency seismic waves.