Patent Application
20090259404
The present invention relates generally to methods and systems for investigating subsurface formations using acoustic measurements made in a borehole. More particularly, the present disclosure is directed to methods and systems that provide improved processing of acoustic waveform data for purposes of characterizing subsurface formations.
Acoustic measurements using currently available tools and methods provide key information in oilfield borehole logging. The currently available acoustic tools are useful in providing a wide range of information regarding the surrounding formation and the borehole parameters. Surface seismic and vertical seismic profiling (VSP) methods are utilized to provide imaging of the overall geological structure of a hydrocarbon reservoir. Sonic and other well logging methods provide good resolution imaging in the immediate vicinity of the borehole. Sonic imaging is yet another technique that bridges the gap in spatial resolution between these seismic and well logging methods.
Some tools include a single source of sonic waves and two or more receivers, however, most of the tools now include two or more acoustic sources and many receivers arranged in an array. A primary use of acoustic borehole measurements is the estimation of compressional (P) wave and/or shear (S) wave formation slowness. The estimation of compressional and/or shear wave formation slowness is often expressed as an S-T (slowness vs. time) plane, and can be visualized at the wellsite with current technology.
While sonic imaging generally has been successful, waveform data acquired for sonic imaging purposes typically contain many types of arrivals in addition to the desired reflected arrivals, such as tool-borne noise and borehole-borne noise. For example, direct compressional and shear headwaves and tube waves, which are a part of the raw data that are acquired during downhole acoustic measurements, tend to obscure the reflected arrival waves relating to acoustic reflectors in subterranean formations.
The assignee of the present application, Schlumberger, has developed a wireline sonic imaging tool (referred to as the Borehole Acoustic Reflection Survey (BARS) tool) that allows reservoir features such as reflectors and fractures to be imaged.
Borehole acoustic reflection surveys are performed, for example, for the characterization of local geology, validation of geo-steering, and detection of fractures. The acquired sonic waveforms are processed and migrated to obtain images up to about 10 to 20 meters away from the borehole. In this, it is desirable to have improved techniques for modeling the dominant wavefields, such as the direct compressional, shear and Stoneley waves, for example, so that they may be removed, making subsurface imaging and characterization efficient and reliable.
Typically, in sonic data processing the reflected (p-to-p, p-to-s and s-to-p) and refracted (p-to-s and s-to-p) waves are usually weak compared to the direct compressional, shear and Stoneley waves. One important aspect of sonic data processing involves the efficient separation of wavefields to reliably extract the reflected and refracted signals. Known techniques for wavefield separation include the adaptive interference canceller (AIC) filter and median velocity filters in common offset gathers, i.e., gathers of traces recorded by common receiver stations of the sonic array tool. However, in typical sonic data processing to extract the reflected and refracted signals event signals whose apparent velocity is slower than the direct P-waves may be removed and the amplitudes, especially for events from layers which are parallel to the well, tend to become weaker. Therefore, it is also desirable to have techniques for modeling the dominant wavefields with greater precision so that extraction of the reflected and refracted wavefields conserves amplitudes of event signals.
The limitations of conventional data processing noted in the preceding are not intended to be exhaustive but rather are among many which may reduce the effectiveness of previously known techniques. The above should be sufficient, however, to demonstrate that acoustic data processing techniques existing in the past will admit to worthwhile improvement.
The present invention meets the above-described needs and others. Specifically, the present disclosure provides methods and systems for taking acoustic measurements relating to subterranean formations, and, in particular, for logging subterranean formations.
The present applicants recognized that sonic data have characteristic properties such as amplitude variations due to geometrical spreading, dispersion, and attenuation. Such amplitude variations within an array may be significant and could be utilized for purposes of improved, efficient modeling of dominant waveforms in acoustic data measurements.
The present applicants further recognized that it is possible to adapt and apply wavefield separation techniques using parametric inversion for purposes of sonic data processing. Applicants further recognized that in sonic data processing the dominant waves to be modeled and removed are three, viz., P-, S- and Stoneley. All three waves have the same sign of moveout and the input data are not multi-component. Additionally, amplitude variations within the array are significant and the waves may be dispersive.
Applicants have developed novel and improved data processing techniques that are applicable to acoustic logging data, for example, sonic logging data, and overcome, or at least reduce the effects of, one or more of the problems outlined above.
Accordingly, an object of the present disclosure is to provide improved systems and methods for processing acoustic waveform data acquired in a borehole in which the dominant waves are accurately and precisely modeled. A further object of certain embodiments herein is to provide improved systems and methods so that reflected and refracted arrival waves are easily identified and selected from tool-borne and borehole-borne arrivals. A further object of certain embodiments herein is to provide methods and systems for characterizing subsurface formations, validating geo-steering, identifying fractures, among other applications, using sonic waveform data acquired in a borehole to derive novel and improved results relating to the subsurface formations.
In one aspect of the present disclosure, a method of processing acoustic waveform data relating to a subterranean formation is provided. The method comprises taking acoustic measurements at one or more depths of a borehole traversing a subterranean formation; deriving slowness of the formation based on wavefields in the acoustic measurements; and modeling dominant wavefields in the acoustic measurements based on the formation slowness, wherein deriving the formation slowness comprises parametric inversion for complex and frequency dependent slowness. The derived slowness is complex, having real and imaginary parts. The dominant wavefields in the acoustic measurements comprise compressional, shear, Stoneley and casing mode waveforms.
In certain embodiments herein, the method comprises decomposing the dominant wavefields in the acoustic measurements. In other embodiments, the method comprises deriving geometrical spreading of direct compressional and shear waveforms in the acoustic measurements. In further embodiments of the present disclosure, the geometrical spreading of the direct compressional and shear waveforms is derived before deriving slowness of the formation.
Aspects of the present disclosure provide iteratively minimizing misfit between input and modeled waveforms. In certain embodiments, the misfit between input and modeled waveforms is minimized by determining a set of geometrical spreading and slowness parameters that minimize the misfit. Dominant wavefields are generated based on the set of geometrical spreading and slowness parameters. Further aspects herein provide subtracting the dominant wavefields from the input waveforms and extracting event signals from the acoustic measurements comprising reflected and/or refracted waveforms.
In some embodiments, a method of processing acoustic waveform data comprises deriving attenuation parameter Q for the subterranean formation. In other embodiments, a method of processing acoustic waveform data comprises analyzing dispersion of Stoneley waveforms for the subterranean formation.
In aspects of the present disclosure, the measurements taken may be sonic measurements. In other aspects, the measurements taken may be logging measurements. In yet other aspects, the measurements taken may be while-drilling measurements. In further aspects herein, the measurements taken may be borehole acoustic reflection survey measurements.
Aspects of the present disclosure provide iteratively minimizing misfit between input and modeled waveforms; subtracting the modeled waveforms from the input waveforms; extracting event signals from the acoustic measurements comprising reflected and/or refracted waveforms; and generating, in real-time with the acoustic measurements, an indication or imaging of acoustic reflectors in the formation.
A method of logging a subterranean formation is provided comprising generating acoustic waves with a source; receiving acoustic waveforms with a plurality of receivers; deriving slowness of the formation based on the acoustic waveforms; modeling dominant waveforms in the acoustic waveforms based on the formation slowness, wherein deriving the formation slowness comprises parametric inversion for complex and frequency dependent slowness and the derived complex slowness has real and imaginary parts; iteratively minimizing misfit between input and modeled waveforms; subtracting the modeled waveforms from the input waveforms; and extracting event signals from the acoustic waveforms comprising reflected and/or refracted waveforms. In aspects herein, the method of logging a subterranean formation further comprises providing at a wellsite, in real-time with receiving acoustic waveforms, an indication or imaging of acoustic reflectors in the formation; and utilizing the indication or imaging of acoustic reflectors for adjusting parameters of the subterranean logging. The indication or imaging of acoustic reflectors in the formation may be utilized for validation of geo-steering parameters.
A system for taking sonic measurements relating to a subterranean formation is provided. The system comprises an acoustic tool comprising at least one source and a plurality of receivers mounted thereon; a computer in communication with the acoustic tool; and a set of instructions executable by the computer that, when executed, receive acoustic waveforms with the plurality of receivers; derive slowness of the formation based on the acoustic waveforms; model dominant waveforms in the acoustic waveforms based on the formation slowness, wherein deriving the formation slowness comprises parametric inversion for complex and frequency dependent slowness and the derived complex slowness has real and imaginary parts; iteratively minimize misfit between input and modeled waveforms; subtract the modeled waveforms from the input waveforms; and extract event signals from the acoustic waveforms comprising reflected and/or refracted waveforms.
Additional advantages and novel features will be set forth in the description which follows or may be learned by those skilled in the art through reading these materials or practicing the invention. The advantages may be achieved through the means recited in the attached claims.
The accompanying drawings illustrate certain embodiments and are a part of the specification. Together with the following description, the drawings demonstrate and explain some of the principles of the present invention.
Throughout the drawings, identical reference numbers and descriptions indicate similar, but not necessarily identical, elements. While the principles described herein are susceptible to various modifications and alternative forms, specific embodiments have been shown by way of example in the drawings and will be described in detail herein. However, it should be understood that the invention is not intended to be limited to the particular forms disclosed. Rather, the invention includes all modifications, equivalents and alternatives falling within the scope of the appended claims.
Illustrative embodiments and aspects of the invention are described below. It will of course be appreciated that in the development of any such actual embodiment, numerous implementation-specific decisions must be made to achieve the developers' specific goals, such as compliance with system-related and business-related constraints, that will vary from one implementation to another. Moreover, it will be appreciated that such development effort might be complex and time-consuming, but would nevertheless be a routine undertaking for those of ordinary skill in the art having the benefit of this disclosure.
Reference throughout the specification to “one embodiment,” “an embodiment,” “some embodiments,” “one aspect,” “an aspect,” or “some aspects” means that a particular feature, structure, method, or characteristic described in connection with the embodiment or aspect is included in at least one embodiment of the present invention. Thus, the appearance of the phrases “in one embodiment” or “in an embodiment” or “in some embodiments” in various places throughout the specification are not necessarily all referring to the same embodiment. Furthermore, the particular features, structures, methods, or characteristics may be combined in any suitable manner in one or more embodiments. The words “including” and “having” shall have the same meaning as the word “comprising.”
Moreover, inventive aspects lie in less than all features of a single disclosed embodiment. Thus, the claims following the Detailed Description are hereby expressly incorporated into this Detailed Description, with each claim standing on its own as a separate embodiment of this invention.
Turning now to the drawings, wherein like numerals indicate like parts, the disclosure herein is directed to the concept of various techniques that may be utilized to facilitate and improve acoustic data processing, and the characterization of subsurface formations based on results that are derived from the improved data processing. The present disclosure contemplates applicability of the disclosed techniques to sonic data that are acquired in, for example, sonic logging operations for purposes of characterization and monitoring of underground reservoirs. The acoustic data may be acquired in exploration and/or production wells. The present disclosure further contemplates application of the disclosed techniques to single well seismic measurements.
The methods and systems may be applied to acoustic wellbore logging with a wireline logging tool. However, the methods and systems presented herein are not so limited. For example, the methods and systems may be applied to other applications such as logging while drilling (LWD), measurement while drilling (MWD), production logging. The present disclosure contemplates the use of the methods and systems herein in combination with various acoustic measurements, for example, with other seismic logging operations.
As described in greater detail below, the present disclosure provides various techniques which may be used to facilitate and improve wavefield estimation and separation. In one aspect of the present disclosure, wavefields recorded by an array of receivers are assumed to be expressed by the phase shift due to the difference of arrival times and a function of geometrical spreading. The wavefields are decomposed as estimated amplitudes of angular frequency components. In further aspects herein, acoustic data are reconstructed using the estimated amplitudes, which contain the decomposed wavefields. If the decomposed wavefields are accurate extractions, reflected and refracted event signals are separated out by subtraction of the reconstructed data from the input data.
Applicants discovered that novel and improved results are obtained when parameters of geometrical spreading and complex-valued slowness of the direct compressional, shear and Stoneley waves are inverted by minimizing the errors between the input and reconstructed waveforms. Aspects of the present disclosure teach that the inverted results may be used for purposes of parametric wavefield decomposition and reconstruction of waveforms using, for example, the three wavefields.
The present disclosure also contemplates reconstruction of casing mode waveforms. In this, casing modes contain the P-, S- and Stoneley waves that are characterized by the P- and S-slowness of the casing.
In further aspects, reconstructed waveforms may be subtracted from the input waveforms, and signals of reflected or refracted events may be separated out. Residuals in the separated waveforms may be removed further by median filters using the estimated apparent velocities.
The present disclosure contemplates sonic data processing to derive reconstructed dominant wavefields. Further, precise estimation of the dominant wavefields provides the capability for reliable and efficient wavefield separation of reflected and/or refracted waveforms. In addition, inverting for slowness according to the techniques of the present disclosure provides enhanced answer products, such as quality factor Q and dispersion analysis of Stoneley waves.
As mentioned above, in the past acoustic reflectors extracted from data acquired during a logging operation could be compromised by subtracting estimates of the dominant waveforms that under certain circumstances could be inaccurate or imprecise.
The present applicants recognized that sonic data have characteristic properties such as amplitude variations due to geometrical spreading, dispersion, and attenuation. Such amplitude variations within an array may be significant and could be utilized for purposes of improved, efficient modeling of dominant waveforms in acoustic data measurements for characterization of subsurface formations, validation of geo-steering, identification of fractures, among other applications that are known to those skilled in the art.
The applicants have proposed a novel technique in which direct P-, S- and Stoneley waves are decomposed. In one embodiment of the present disclosure, a geometrical spreading function may be utilized to solve parameters for direct P- and S-waves. The present disclosure contemplates parametric inversions for both real and frequency independent parameters, i.e., non-dispersive, and complex and frequency dependent slowness, i.e., dispersive.
In order to gain a better understanding of the various techniques and features described in this application, a brief description of the disclosed techniques will now be provided. The present disclosure provides methods and systems for wavefield separation using parametric inversion and decomposition for data acquired, for example, in a borehole acoustic survey. Parameters describing the propagation along the wellbore of compressional, shear, and Stoneley waves are estimated, including geometrical spreading and complex-valued velocities. The maximum portion of the measured wavefields, compatible with the estimated parameters, are identified and subtracted. A further attenuation of the direct waves is obtained using median filters along the estimated propagation velocities. After the direct waves are removed, a relatively larger presence remains of components of the wavefield that are possibly reflected in the formation. Applicants conducted tests, as described in further details hereinafter, to confirm the effectiveness of the techniques by application to real field data.
In some embodiments of the present disclosure, assuming each wavefield can be expressed by a phase shift in the frequency domain due to the differences in arrival times and a function of geometrical spreading, the parametric decomposition is expressed as:
Geometrical spreading is derived for P- and S-waves (it is constant for Stoneley waves) in one embodiment of the present disclosure by inversion:
g
k(z)=(1+ck(z−zmin))−1 Equation 2
The present disclosure also contemplates use of a geometrical spreading function such as:
g
p(z)=(
g
s(z)=
Amplitudes may be estimated by:
Gm=d
m=(G+G+εI)−1G+d Equation 4
The present disclosure contemplates parametric inversion of dispersive (complex) equations:
{circumflex over (d)}−d=0
H δp=({circumflex over (d)}−d) Equation 5
δp=−(H+H+εI)−1H+({circumflex over (d)}−d)
The present disclosure also provides parametric inversion of non-dispersive (real) equations:
{circumflex over (d)}=Gm.
E=({circumflex over (d)}−d)+({circumflex over (d)}−d) Equation 6
After shifting the waveforms by times appropriate for each of the three dominant direct waves—compressional, shear and Stoneley—a spatial median filter is used to remove residuals effectively left by the processing discussed above.
The present disclosure provides inverting for slowness according to the techniques herein to derive quality factor Q as follows:
Referring to
Referring again to
In the processing workflow of
In the processing workflow of
As described in further detail below, the techniques illustrated in the flowcharts of
Referring to
An array tool of a sonic logging system may be used for the acquisition of acoustic data. In this, recording time may be set to be sufficiently long to record the reflected waves.
The principles described herein may be applied to various types of data acquisition systems, such as wireline measurements, drilling measurements, production measurements, and so on.
The techniques described above may be used with logging tools and systems of any type. For example, referring again to
As one possibility, Schlumberger's wireline sonic imaging tool referred to as the Borehole Acoustic Reflection Survey (BARS) tool may be utilized for the purposes described herein.
The receivers R and source S are coupled to a computer processor (generally designated as 236 in
Measurements collected from the sonic tool 230, which may include waveforms generated by the receivers R over time, may be sent to the surface via a cable (not shown) from which the sonic tool is suspended or by any other convenient communication technique (e.g., mud telemetry or wired drillpipe). The computer processor 236 is commercially available from a wide variety of sources. The sonic data taken by the sonic tool 230 and received by the computer processor 236 (in-situ or at the surface) may be processed according to instructions accessible to the computer processor 236 to perform the processing described above. The imaging of acoustic reflectors may be performed at the wellsite, although in other embodiments described herein computations may be done at an off-site processing location using the sonic data collected by the sonic tool.
In
In one example of application of the acoustic data processing techniques of the present disclosure, Generalized Radon Transform migration (Miller, D., Oristaglio, M., and Beylkin, G., 1987, A new slant of seismic imaging: Migration and integral geometry, Geophysics, 52, 943-964) was used for PrP (P-reflected-P) and SrS (S-reflected-S) waves. Waveforms that were separated and processed using semblance weighted deconvolution (SWD) were used as input waveforms. 1-D velocity model (no dip) was created from inverted slowness by parametric inversion. The migration options were 0 degree dip and 60 degree aperture. The results are depicted in
As mentioned above, according to some aspects of the present disclosure, the methods and systems described herein provide robust, accurate, usable estimation or modeling of dominant acoustic wavefields with a relatively simple process that requires very little computer resources. The techniques described above enable an operator to see real-time indications of possible reflectors in the surrounding subsurface formations, and to act appropriately during the course of a logging operation.
Generally, the techniques disclosed herein may be implemented on software and/or hardware. For example, they can be implemented in an operating system kernel, in a separate user process, in a library package bound into network applications, on a specially constructed machine, or on a network interface card. In one embodiment, the techniques disclosed herein may be implemented in software such as an operating system or in an application running on an operating system.
A software or software/hardware hybrid implementation of the present techniques may be implemented on a general-purpose programmable machine selectively activated or reconfigured by a computer program stored in memory. Such a programmable machine may be implemented on a general-purpose network host machine such as a personal computer or workstation. Further, the techniques disclosed herein may be at least partially implemented on a card (e.g., an interface card) for a network device or a general-purpose computing device.
The embodiments and aspects were chosen and described in order to best explain the principles of the invention and its practical applications. The preceding description is intended to enable others skilled in the art to best utilize the principles described herein in various embodiments and with various modifications as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the following claims.
This application claims priority of U.S. Provisional Patent Application Ser. No. 60/043,401, filed 9 Apr. 2008, the entire contents of which are incorporated herein by reference.
| Number | Date | Country | |
|---|---|---|---|
| 61043401 | Apr 2008 | US |
