Click to view article in PDF format.
GC3-D Seismic Data in Imaging Fracture Properties for Reservoir Development*
By
Bob Parney1 and Paul LaPointe2
Search and Discovery Article #40098 (2003)
*Adapted for online presentation from the Geophysical Corner column in AAPG Explorer October and November, 2002, entitled “Fractures Can Come Into Focus”and “Simple Seismic, Complex Fractures,” respectively, and prepared by the authors. Appreciation is expressed to the authors, to R. Randy Ray, Chairman of the AAPG Geophysical Integration Committee, and to Larry Nation, AAPG Communications Director, for their support of this online version.
1Axis Geophysics, Denver, Colorado ([email protected])
2Golder Associates, Redmond, Washington ([email protected])
It has long been recognized that the presence of naturally occurring fracture networks can lead to unpredictable heterogeneity within many reservoirs. Conversely, fractures provide high permeability pathways that can be exploited to extract reserves stored in otherwise low permeability matrix rock. One of the primary difficulties in managing fracture heterogeneity and the consequent uncertainty is that production rates and volumes are controlled by fracture network connectivity between the producing wells, while the primary sources of data on fracture properties are measured only in the vicinity of wells. In some ways this is like trying to predict the size of a schoolyard by close examination of a single link in the surrounding fence.
Recent advances in the processing of 3-D seismic data, however, are providing valuable new tools for the imaging of fracture properties between wells. Those tools are the analysis of seismic velocities as affected by raypath direction and offset distance. Specifically, adjusting velocities as a function of azimuth (velocity anisotropy) to improve reflection imaging has produced by-product data volumes of seismic velocity anisotropy (ANMO) and improved data volumes of azimuthal changes in amplitude as a function of offset (AVAZ).
These seismic advances raise the following questions:
-
How do fractures influence these data?
-
Geologically, what should this newly imageable level of fracture heterogeneity look like?
-
How do we interpret this new data for fracture properties?
-
How do we then make the link between fracture properties and reservoir performance?
uTheory of seismic response to fractures uFracture orientation in Rocky Mountains uInterpreting multiple fracture set properties
uTheory of seismic response to fractures uFracture orientation in Rocky Mountains uInterpreting multiple fracture set properties
uTheory of seismic response to fractures uFracture orientation in Rocky Mountains uInterpreting multiple fracture set properties
uTheory of seismic response to fractures uFracture orientation in Rocky Mountains uInterpreting multiple fracture set properties
uTheory of seismic response to fractures uFracture orientation in Rocky Mountains uInterpreting multiple fracture set properties
uTheory of seismic response to fractures uFracture orientation in Rocky Mountains uInterpreting multiple fracture set properties
|
Theory of Seismic Response To Fractures
The underlying theory behind the ANMO and AVAZ
processing is quite simple: Most geophysical processing algorithms
assume that all fractures are approximately vertical, and are locally
oriented in a single dominant direction (Figure 1). The maximum
detectable seismic effect is when the seismic raypath travels
perpendicular to the open fractures, crossing the slow velocity,
possibly fluid-filled, open fracture. A maximum and minimum direction of
fracture influence on P- The difference between the maximum and minimum effect gives some measure of the fracture intensity. This same process can be applied in a number of data volumes where the change in Vp or Vs as a function of azimuth is measured by the change in stacking velocities (azimuthal NMO) or the change in reflection coefficients (azimuthal AVO).
Fracture Orientation in Rocky MountainsA critical feature of recently processed AVAZ and ANMO data volumes has been that the dominant fracture orientation can change dramatically over short distances. Recent work on a project sponsored by the U.S. Department of Energy (www.fracturedreservoirs.com) shows that these changes are not only possible, but also highly likely in a Rocky Mountain compressional setting where the stress field is complex. The Circle Ridge Field, in Wyoming’s Wind River Reservation, Wind River Basin (Figure 2), was characterized through a combination of 2-D cross-sections and 3-D structural reconstructions based on well and surface data, and fracture data from surface outcrops and subsurface image logs. The fracture and structural data were supplemented with data from several transient well tests, a bromide tracer test and a nitrogen injection test. The structure is primarily determined by NE-SW compression, which caused the formation of a series of imbricate fault blocks along the Red Gully Fault, including several imbricates to the north (Figure 3). The entire structure has been characterized as a fault-breached, fault-propagation fold. Development of the structure is likely to have produced the fracturing within the reservoir units. Fracture development was predicted using strain calculated through a 3-D palinspastic reconstruction of the field. Figure 4 shows differences in extensional strain magnitude and orientation throughout a block of the Tensleep Formation in the hanging wall of the field’s Red Gully Fault. The contours and line lengths represent the magnitude of the maximum extensional strain due to the initial folding of the reservoir formations. The figure’s red lines represent the strike orientation of extensional fractures that would develop perpendicular to the local direction of maximum extensional strain. The red lines also show the dominant set; it is likely that a secondary joint set perpendicular to the set shown might also develop. Ninety-degree changes in dominant fracture orientation across fracture fairways seen in Figure 4 are consistent with orientation patterns predicted by AVAZ data in nearby reservoirs. These orientation variations arise due to inhomogeneities in the stress field and the resulting fracture networks are consistent with well image log and tracer data. Similar changes in fracture orientation occur in nearby outcrop at a much smaller scale (Figure 5). The black fractures occur only on the left portion of the outcrop, nowhere else. Red fractures dominate over blue fractures in the left portion, while blue fracture intensity increases markedly on the right hand side. Since seismic anisotropy can be influenced by the presence of natural fractures – and that a high degree of variability in fracture orientation and intensity is to be expected in a Rocky Mountain compressional setting – interpretation of seismic data requires a sound link with knowledge of the fracture geology in a region.
Interpreting Multiple Fracture Set PropertiesThe determination of fracture azimuth and intensity is usually based on the assumption that there is a single dominant fracture orientation, typically vertical. Frequently, fractures occur in several sets with cross-cutting orientations (Figure 6), and generally multiple sets are necessary in order to get well-connected plumbing for long-term productivity in the absence of high matrix permeability.
A number of attributes can be extracted from the seismic data. They can be grouped into two major categories:
Orientation attributes such as the fast P or S
Anisotropy: Fractures or No-Fractures
In the early development of anisotropic
seismic analysis it was thought that high levels of anisotropy, as
measured by the difference between the fast and slow P and S
Thus, the magnitude of the anisotropy does not
in itself differentiate between regions of high fracture intensity and
low fracture intensity. Other attributes such as interval velocity must
be
The Next Step: Calculating Permeability
Once the attributes of the natural fracture
system have been mapped, the next step is to take these attributes and
use them as a predictive tool. This process, however, is not as simple
as identifying fracture properties at a In order to assess the connectivity of a reservoir, the next step after obtaining the fracture attributes from the seismic data is to use DFN models to understand the consequences of fracture orientation and intensity on permeability. The DFN approach models fractures as two-dimensional polygonal planar objects, like playing cards, located in three-dimensional space (Figure 8a). Each fracture is characterized by its surface area and shape and has flow properties such as permeability, compressibility and aperture. Network models can be formed based on an interpretation of seismic attribute data, engineering data, and image log data as available. Once fractures are generated, a finite element mesh can be constructed according to the fracture geometry (Figure 8b), and a flow solution can be obtained that takes into account the connectedness of the fracture system. Figure 9 shows an example of a pressure pulse spreading through a fractured reservoir in response to injection.
Seismic + Fractures = Permeability Prediction
Discrete Fracture Network (DFN) models have
provided an important tool to make the connection between seismic
properties and reservoir. The DFN approach can be combined with seismic
attribute mapping by first developing an interpretation of the link
between attributes and fracturing. For example, the difference between
the fast and slow P or S velocities can be
Once the DFN model has been generated, a grid
can be placed over the model and a finite element mesh
SummaryRecent advances in the processing of 3-D seismic data are providing valuable new tools for quantifying fracture properties between wells. In order to make use of this new information, it is necessary to:
Although uncertainties abound, these attributes provide new insight into notoriously difficult reservoirs, and promise to enhance recovery through focused engineering efforts.
ReferenceKeefer, W.R., 1969, Geology of petroleum in Wind River Basin, Central Wyoming: AAPG Bulletin, v. 53, p. 1839-1865. |