Anisotropic velocity analysis1

Results from walkaway VSP and shale laboratory experiments show that shale anisotropy can be significantly anelliptic. Heterogeneity and anellipticity both lead to non-hyperbolic moveout curves and the resulting ambiguity in velocity analysis is investigated for the case of a factorizable anisotropic medium with a linear dependence of velocity on depth. More information can be obtained if there are several reflectors. The method of Dellinger et al. for anisotropic velocity analysis in layered transversely isotropic media is examined and is shown to be restricted to media having relatively small anellipticity. A new scheme, based on an expansion of the inverse-squared group velocity in spherical harmonics, is presented. This scheme can be used for larger anellipticity, and is applicable for horizontal layers having monoclinic symmetry with the symmetry plane parallel to the layers. The method is applied to invert the results of anisotropic ray tracing on a model Sand/shale sequence. For transversely isotropic media with small anisotropy, the scheme reduces to the method of Byun et al. and Byun and Corrigan. The expansion in spherical harmonics allows the P-phase slowness surface of each layer to be determined in analytic form from the layer parameters obtained by inversion without the need to assume that the anisotropy is weak.     

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.     

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.     

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.     

Constructing a discrete fracture network constrained by seismic inversion data

Rock fractures are of great practical importance to petroleum reservoir engineering because they provide pathways for fluid flow, especially in reservoirs with low matrix permeability, where they constitute the primary flow conduits. Understanding the spatial distribution of natural fracture networks is thus key to optimising production. The impact of fracture systems on fluid flow patterns can be predicted using discrete fracture network models, which allow not only the 6 independent components of the second‐rank permeability tensor to be estimated, but also the 21 independent components of the fully anisotropic fourth‐rank elastic stiffness tensor, from which the elastic and seismic properties of the fractured rock medium can be predicted. As they are stochastically generated, discrete fracture network realisations are inherently non‐unique. It is thus important to constrain their construction, so as to reduce their range of variability and, hence, the uncertainty of fractured rock properties derived from them. This paper presents the underlying theory and implementation of a method for constructing a geologically realistic discrete fracture network, constrained by seismic amplitude variation with offset and azimuth data. Several different formulations are described, depending on the type of seismic data and prior geologic information available, and the relative strengths and weaknesses of each approach are compared. Potential applications of the method are numerous, including the prediction of fluid flow, elastic and seismic properties of fractured reservoirs, model‐based inversion of seismic amplitude variation with offset and azimuth data, and the optimal placement and orientation of infill wells to maximise production.     

Teleprogramming for subsea teleoperation using acousticcommunication

This paper considers the performance of subsea intervention tasks from an unmanned untethered submersible while using acoustic communications. It is argued that the low bandwidth and high delay imposed by acoustic modems makes it unwise to adopt conventional teleoperation techniques and a system is presented which permits subsea teleoperative tasks to be carried out using such limited communication resources. The described implementation employs active techniques to assist the operator both in performing actions and in recovering from those problems which will inevitably occur during real-world interaction. It provides the operator with both simulated and real visual imagery and is designed to adapt dynamically to changing bandwidth and computational resources. Experiments are described in which an operator in Philadelphia, PA, controlled a robot manipulator mounted on the JASON underwater vehicle submerged off the Massachusetts coast. All communication over this 500-km distance was via a combination of Internet and a simulated acoustic link. Analysis of the bandwidth requirements showed them to be consistent with those from acoustic subsea networks     

Concerning the electron temperature discrepancy between in situ and remote probes in the E-region

Rocket borne Langmuir probe measurements of electron temperature in the E-region are examined in relation to recent laboratory investigations of surface drift effects which can lead to erroneously high and time-dependent electron temperature measurements. The rocket data is consistent with the laboratory expectations thus supporting the suggested importance of surface effects in rocket measurements and in relation to the E-region discrepancy with simultaneous incoherent radar scatter measurements.     

Fluid-dependent shear-wave splitting in a poroelastic medium with conjugate fracture sets

The dependence of shear‐wave splitting in fractured reservoirs on the properties of the filling fluid may provide a useful attribute for identifying reservoir fluids. If the direction of wave propagation is not perpendicular or parallel to the plane of fracturing, the wave polarized in the plane perpendicular to the fractures is a quasi‐shear mode, and therefore the shear‐wave splitting will be sensitive to the fluid bulk modulus. The magnitude of this sensitivity depends upon the extent to which fluid pressure can equilibrate between pores and fractures during the period of the deformation. In this paper, we use the anisotropic Gassmann equations and existing formulations for the excess compliance due to fracturing to estimate the splitting of vertically propagating shear waves as a function of the fluid modulus for a porous medium with a single set of dipping fractures and with two conjugate fracture sets, dipping with opposite dips to the vertical. This is achieved using two alternative approaches. In the first approach, it is assumed that the deformation taking place is quasi‐static: that is, the frequency of the elastic disturbance is low enough to allow enough time for fluid to flow between both the fractures and the pore space throughout the medium. In the second approach, we assume that the frequency is low enough to allow fluid flow between a fracture set and the surrounding pore space, but high enough so that there is not enough time during the period of the elastic disturbance for fluid flow between different fracture sets to occur. It is found that the second approach yields a much stronger dependence of shear‐wave splitting on the fluid modulus than the first approach. This is a consequence of the fact that at higher wave frequencies there is not enough time for fluid pressure to equilibrate and therefore the elastic properties of the fluid have a greater effect on the magnitude of the shear‐wave splitting.