Click to view article in PDF format
GCFracture Identification with Seismic Data, XinChang Field, China*
Peter Stewart1, John Tinnin1, James Hallin1, and Jim Granath1
Search and Discovery Article #40370 (2008)
Posted December 9, 2008
*Adapted from the Geophysical Corner column, prepared by the authors, as a two-part series in AAPG Explorer, October, 2008, entitled “China Study: Detecting Fractures”, and November, 2008, entitled “China Study: Fracture Data Integrated”. Editor of Geophysical Corner is Bob A. Hardage ([email protected]). Managing Editor of AAPG Explorer is Vern Stefanic; Larry Nation is Communications Director.
1 ION Geophysical/GX Technology, Houston, TX
This paper describes how Sinopec’s local operating company, Southwest Petroleum Branch (SWPB), utilized full-wave seismic data to improve production from a fractured tight-gas reservoir in XinChang Field, Sichuan Province, China. First we detail the data-acquisition technology and the data-processing workflow that produced high-resolution images and yielded fracture information that correlated with well production.
Historically, this region has been a prolific gas producer – shallow prospects were depleted early, and the reservoirs currently targeted are now at the base of a terrestrial sequence some 20,000 feet thick. These deeper Triassic reservoirs are low porosity (less than 4 percent) – but specific areas within the reservoir can be highly fractured. Production has been declining, and the region now needs an injection of new technology to sustain production.
Legacy seismic data correlate poorly with existing wells, and the quality of existing seismic data is insufficient to define reservoir targets. Attention was focused on implementing a seismic program that would allow the fracture network to be understood so future drilling locations could be determined. In this effort, a task force of ION and SWPB geoscientists found that the region produces high levels of coherent converted-shear (C-wave) energy. The team concluded that C-waves had the potential of providing stratigraphic, lithologic, and fracture detail that would be crucial for understanding the reservoir and for optimizing well placements and reducing drilling risk.
The design team recommended a data-acquisition program involving dense spatial sampling, full offset and azimuth distributions, and the adoption of 3C digital sensors. With the design approved, a new survey was acquired in 2004 using
an I/O System Four® recording system and VectorSeis® full-wave 3C sensors. It became apparent shortly after data-acquisition began that the new P-wave data were high quality and that bandwidth and signal-to-noise ratios were a step change improvement over legacy seismic data. In addition, high-quality, full-azimuth C-wave data were also recorded.
|
The data-processing workflow that was implemented resulted in high-resolution C-wave images. An important byproduct of this workflow was information related to fracture orientation and to fracture intensity. The first step in the data processing was to rotate horizontal components from their recorded field orientations to their source-receiver azimuths, or “radial,” directions. If seismic data are azimuthally isotropic, all C-wave reflection energy would be concentrated on the radial component, and data acquired by the transverse sensor could be discarded. However, after rotating the XinChang data to radial/transverse coordinates, a significant amount of C-wave energy remained on the transverse component. This data behavior confirmed the presence of shear-wave splitting, which occurs when a shear wave encounters an azimuthally anisotropic layer such as one of the XinChang fractured reservoir units. In S-wave splitting, the C-wave polarizes into two new waves, a phenomenon known as “birefringence.” One of the new split waves is polarized parallel to, and the second orthogonal to, the fracture orientation. The velocities of the two waves differ – the faster wave being polarized parallel to the fractures, and the slower wave polarized orthogonal to the fractures (Figure 1). In addition, each of the new waves splits again when it encounters a new anisotropic layer, resulting in a complicated mix of waves arriving at each sensor. Shear-wave splitting can yield valuable information regarding fractures; however, unless addressed correctly, wave-splitting reduces the bandwidth of stacked or migrated C-wave images. Thus, our data-processing workflow was modified to capitalize on the shear-splitting that was detected. After sensor rotation to radial/transverse coordinates, the main components of the workflow included surface wave attenuation, resolution of shear-wave statics, surface-consistent signal processing and Q compensation. These steps were performed independently on the radial and transverse components. Next, the data for each component were subdivided into 36 10-degree, azimuth sectors. Each azimuth sector was migrated separately via a prestack C-wave time migration. This migration step required a Following migration, each sector volume was subjected to residual move-out correction, muting and stacking. The azimuth volumes were then re-assembled into azimuth-sector The most important step in the data processing was our layer-stripping anisotropic correction. This procedure removed the effects of shear-wave spitting at each anisotropic boundary by simulating the effect of an isotropic medium. A single anisotropic layer was stripped to form an isotropic layer in the following manner: Step 1. Knowing that azimuths corresponding to polarity reversals observed on the transverse component define fracture orientation, these azimuth angles were used to rotate the data from radial/transverse to fast/slow directions. Step 2. A cross-correlation between fast and slow data determined the time lag between these two wave modes; a static correction was then done to time align slow and fast data. Step 3. An additional rotation back to radial and transverse coordinates concentrated all of the energy onto the radial component and produced azimuthally isotropic data. These adjusted data were stacked to form a high-resolution C-wave image (Figure 3b). An important byproduct of step one is fracture orientation. For any particular layer, maps of fracture orientation throughout the entire data volume were generated Post-Processing/Interpretation Upon completion of the 3D3C seismic data processing, we initiated post-processing and interpretation activities. This portion of the workflow involved integrating available well, outcrop, and core data with the processed 3D3C seismic data. The post-processing included:
These steps built better structural and stratigraphic models, mapped fracture patterns and intensity, and provided an improved understanding of the region’s geologic and tectonic history. The Sichuan Basin underwent dramatic subsidence rates during the early Mesozoic. Burial of the reservoir to depths of 20,000 feet or more occurred soon after deposition and before gas was generated in surrounding and underlying source rocks. The resulting compaction reduced reservoir porosities to less than 4 percent, causing reservoir rocks to be almost impermeable. As a result, production in XinChang Field is fracture-dependent. The fracture network:
Interbedded sand-shale sequences are the best exploration targets – these thinner-bedded, brittle layers fracture more easily and with higher density than do their thicker counterparts. The integration of geological history with production data resulted in a model showing that storage capacity in the reservoir depends on interconnection of fractures in fault-damage zones and on the connections that these damage zones make with naturally fractured sandstone reservoir beds. The densely sampled, full-azimuth P-wave data acquired in this study supplied higher frequencies than existed in legacy seismic data, resulting in an improved structural picture with excellent fault resolution. However, because sands and shales had similar acoustic impedances, lithology could not be determined from P-wave data alone. Fortunately, converted-shear (C-wave) data were valuable for discriminating lithology in these rocks and delivered vital insights into the stratigraphic architecture of targeted reservoirs. The C-wave data provided a means to define interbedded, sand-shale sequences that were areas of optimal fracture intensity and enabled delineation of source rock (Figure 4). Because fractures dominate storage and movement of gas in XinChang Field, developing a tool to identify and map the best fractured zones was a high priority. The C-wave dataset proved to be that needed tool, because the data provided azimuthal definitions of S-wave As illustrated in Figure 5, seismic wave propagation is minimally affected parallel to the dominant fracture trend in rocks that have a simple one-directional fracture system (the “fast direction,” Vfast), but a maximum As described above, fracture orientation was determined by deriving the azimuths of Vfast from shear-wave splitting analyses performed in a layered-Earth approach. In the XinChang survey area, regional geology, borehole breakout and FMI log results all indicate that the current maximum horizontal stress is oriented along azimuths of 80-110 degrees. This information validated the shear-wave splitting results for orientation in the uppermost Earth layer, which has an average orientation of 95 degrees (Figure 6). It also is important to note that fractures oriented close to this principal stress direction are more likely to be open, though the extreme overpressure in this area keeps other fractures open as well. Knowing fracture orientation in zones of higher anisotropy, where a single set of parallel fractures is more likely to exist, can help in designing directional or horizontal wells that will intersect more fractures, yielding higher production in these areas. Fracture density was determined from the Insights from similarity processing, curvature attributes and shear-wave splitting The authors thank the management of Sinopec and Southwest Petroleum Branch for granting permission to present this paper – especially Xu Xiangrong, president of Southwest Petroleum Branch Company, for his leadership and commitment to cutting-edge technologies. The authors also thank the other contributors to the interpretation project, including AAPG member Roger Palomino and Doug Allinson, Felix Diaz, Reinaldo Nossa, Santi Randazzo and Jim Simmons.
|