1. Field of Invention
This invention is related to mapping and highlighting on a seismic reflection image, reservoir formations that have higher-viscosity pore fluids like oil, in comparison with reservoir formations that have lower viscosity pore fluids like water and gas. Oil, the higher-viscosity pore fluid, requires more time to go back to its original equilibrium state after the pore fluid has been disturbed by a seismic compressional impulse. The time required to recover to the state of equilibrium after the reservoir rock has been deformed as a result of being exposed to an outside seismic impulse can be identified as a time-constant. The recovery time response, known as time-constant, required to reach the state of equilibrium is greater in the case of oil as a pore fluid, versus water or gas. The time-constant, which represents the decay of the stored potential energy in the pore fluids, generates lower frequencies. More specifically, this invention describes a method of identifying oil-saturated rocks by mapping the lower frequencies in the seismic spectrum. These lower frequencies are generated in the reservoir rocks due to the relative movement between the oil and the rock matrix. These lower frequencies will not be present in the case of water or gas, which have a much lower time-constant. Mapping the presence or absence of these lower frequencies becomes a powerful seismic attribute to locate the oil-bearing formations.
2. Description of the Related Art
At present, seismic reflection methods are universally being used to map the structure of the subsurface formations. The mapped structural information delineates the possible hydrocarbon traps that can be drilled. However, due to the complexity of the sedimentary rocks, it becomes difficult to rely on current seismic data to identify and map the reservoir properties prior to drilling the location of interest. Ideally, one would like to identify the presence and the type of pore fluids, and to determine their location and their extent. Using current seismic methods, it is difficult to extract that type of reservoir information since the only known parameters are velocity, attenuation and modulus. To extract information like reservoir pore fluids, porosity and permeability, one finds that there are more unknown parameters than known parameters in the equations used to resolve them, which causes a great deal of ambiguity in the seismic results. The common product of any set of seismic reflection data is acoustic impedance, which in turn depends on several rock properties including velocity, density, porosity, differential pressure, rock matrix, pore fluids, and fluid compressibility. With that many unknowns and not enough equations to solve them, it becomes obvious that there are problems that one faces when one tries to extract reservoir properties from the seismic reflection data. Subjective decisions are made, which may not represent the reality of the complex subsurface.
At times, S-wave (shear wave) recordings are carried out in addition to P-wave (compressional wave) recordings. This provides additional information to reduce the ambiguity of the results. However, S-wave seismic sources are not available in the market, and industry has to rely on mode-converted shear waves generated at the boundaries of the formations, which have different elastic properties. In areas that have complex geology and structural geometry, it becomes difficult to identify and process mode converted S-waves, since mode conversions from P to S-wave and vice-versa take place in a very unpredictable manner. Even in the case where data is well behaved and P and S-wave velocities can be derived, it is difficult to map certain rock properties like porosity and clay content due to a certain amount of self-similarity between the behaviors of the two waves.
To reduce the level of ambiguity and the risk associated with hydrocarbon exploration, seismic technology known as amplitude-variance-with-offset (AVO) has been introduced. The AVO anomaly is most commonly expressed as an increase in the reflection amplitude with the increase in offset distance. AVO has a reasonable chance of success in mapping gaseous hydrocarbon traps. However, like every technology, it has its own shortcomings. AVO has not been successful in Gulf of Mexico, where the sands are calcium carbonate-cemented, or in the Denver-Julesburg basin and Sacramento Valley. At present the success rate of AVO is claimed to be better than fifty percent.
Recently electromagnetic (EM) methods have been re-introduced to reduce the uncertainty of the seismic results prior to drilling commitment. Lately, controlled source electromagnetic technology (CSEM) has been used in deep water to reduce the risk of drilling expensive wells. The value of CSEM has been to complement the seismic data with direct detection of hydrocarbons in the subsurface. The basic concept of the technology is to measure the resistivity of the subsurface structure. Hydrocarbon-filled reservoirs show higher resistivity compared to water-filled reservoirs. This technology is still in its development stage and may prove to be very useful in the future. However, it is not a stand-alone technology and has to be integrated with reflection seismic to provide meaningful results. The resolution of EM technology is limited; the individual formations, which are oil-saturated, cannot be mapped nor their areal extent accurately defined.
At present, using currently available subsurface imaging methods, seismic or otherwise, more than 60% of the producible oil has been left behind in the reservoirs. The extraction technologies exist; however, the weakness in detecting oil is the lack of reliable subsurface imaging methods to identify the pockets of residual oil left unproduced. The oil industry needs a paradigm shift and needs to recognize that there has to be a greater focus on mapping and distinguishing the actual reservoir properties instead of the rock properties. Rock properties constitute mapping changes in velocity and amplitude anomalies due to changes in reservoir characteristics. Reservoir properties are more closely related to pore structure, pore connections, pore fluids and their flow characteristics. These properties can be better mapped by inventing and introducing new seismic methods that relate to the nonlinear characteristics of the reservoir rocks, which are related physically to the interaction between the rock matrix and the reservoir fluids as a result of the rock being excited by the controlled seismic signals.
The current seismic practice has made an incorrect assumption of elastic linearity, when dealing with the seismic wave propagation in the reservoir rocks. At present, the seismic industry generally ignores the effects of elastic nonlinearity in the reservoir rocks, which is the most important and differentiating characteristic in identifying the porous permeable and fluid-saturated subsurface formations. This simplified assumption of linearity is quite often based on convenience and to avoid complex and cumbersome mathematics necessary to deal with the nonlinear behavior. Implicit in the assumption of elastic linearity is the fact that the seismic wave or pulse recorded after being reflected or refracted will contain only those frequencies that were part of the original seismic input signal. In the assumption of an elastically linear system, no new frequencies can be generated while the seismic signal is traveling through different earth formations. This approach has shut the door for new developments that could resolve many of the current seismic imaging problems.
Most seismic recordings for hydrocarbon exploration are made with a useable bandwidth of 6 Hz to 8 Hz on the lower end of the seismic frequency spectrum and 70 Hz to 80 Hz on the upper end of the frequency spectrum. The current seismic recording practices are designed for a linear earth model that has been successful for mapping subsurface structure but has not been very effective in mapping reservoir properties like porosity, permeability and pore fluids. Elastic nonlinearity effects in a porous and permeable reservoir rock generate new frequencies that may be as low as 1 Hz or 2 Hz, and could be used as an indicator of higher-viscosity pore fluid like oil. To achieve that, seismic data has to be recorded with a bandwidth that has lower frequencies all the way down to 0 Hz. The current seismic is more suitable for mapping structural details of the subsurface formations, but has serious limitations when trying to extract reservoir properties from the seismic data.
For years, the oil industry has been searching for a technology that will enable it to map subsurface reservoirs that hold oil rather than water or gas. So far, there has not been any significant breakthrough. One of the physical characteristics of the oil is its higher viscosity, which differentiates it from water and gas. This invention uses this higher viscosity, which is oil's differential physical characteristic, to seismically map the oil-bearing subsurface reservoir rocks and highlight them on seismic subsurface reflection images in comparison with lower viscosity water-saturated or gas-saturated rocks.
The presence of the pores and pore fluids, which occupy a significant volume of the reservoir rocks, have a pronounced effect on their dynamic elastic properties. When a seismic impulse is applied to the reservoir rock, there is a local change in its volume, density, and modulus. During the compression cycle of the impulse, the volume of the rock is reduced. In response to the sudden pressure and volume change, the internal flow of the pore fluids takes place in the case of the porous and permeable reservoir rock. This internal flow of the pore fluids causes a part of the energy to be dissipated due to friction, and the remainder is stored in the form of potential energy due to the displacement of the pore fluids. This component of the potential energy, which is stored in the displaced pore fluid, tends to restore the original state of equilibrium, which has been disturbed. The process of restoration, which involves the pore fluids to go back to their original state of equilibrium, requires a certain period of time interval. The recovery time response is related to the viscosity of the pore fluid and the physical geometry of the pores and their interconnections. In the case of water and gas, which have lower viscosity, the restoring time-constant is lower in comparison with oil, which has higher viscosity. Higher-viscosity fluid like oil, shows relatively higher resistance to flow and consequently requires more time to find its state of equilibrium after the pore fluids have been disturbed.
Once an object displaying an elastic characteristic is distorted, its shape tends to be restored to its original equilibrium state. In the same manner, when the pore fluids in the reservoir rocks are forced to move in response to a seismic impulse, the fluids tend to go back to their original state of equilibrium. This recovery to their original state of equilibrium is time-dependent on the flow characteristics of the pore fluids, mainly its viscosity and the shape of the pores and their interconnections. The elastic behavior of the rock matrix and the pore fluids is physically coupled. The energy dissipation and the potential energy stored in the pore fluids due to their movement caused by a seismic impulse directly affects the elastic property of the reservoir rock. The elastic behavior of the reservoir rock, saturated with high-viscosity pore fluid like oil, will be different than the reservoir rock saturated with water or gas. After a seismic impulse has been applied, the restoring time-constant of the oil-saturated reservoir rock will be larger compared to the restoring time-constant of the reservoir rock saturated with water or gas.
The larger time-constant stretches the waveform of the seismic impulse and generates lower frequencies that are not part of the seismic input signal. The higher-viscosity pore fluids like oil cause more pulse stretching since they exhibit a larger time-constant. Consequently, more pronounced lower frequencies are generated when a seismic pulse propagates through an oil-saturated reservoir compared to when the same formation is saturated with water or gas.
This invention provides a method for seismic reflection mapping the location and the extent of the oil-bearing reservoir formations. The measurement of the relative amplitudes of the lower frequencies in the seismic spectrum of a particular subsurface horizon, which are in the order of 2 Hz to 6 Hz, is a strong indicator of the presence of high-viscosity pore fluid like oil. The dominance of the relative amplitudes of these lower frequencies in the seismic-reflected signals of a particular subsurface formation is indicative that these frequencies are being generated in that particular reservoir formation. By mapping and displaying the subsurface seismic reflection image using various frequency filtering methods available in the industry today, one can highlight the oil-bearing rocks and identify their location and extent. This method of mapping oil-bearing subsurface formations is equally applicable for land or marine operations. According to the current and accepted art, surface or borehole seismic can be used for the method described in this invention. To provide more distinct identification of the frequencies being generated in the oil-saturated reservoir rock against the seismic source-generated frequencies, the seismic source has to generate a short duration impulse that is rich in frequencies above 10 Hz, and weak in generating frequencies below 6 Hz. This is not a major problem since most land and marine seismic sources currently being used are designed to generate higher frequencies to provide a higher resolution subsurface image and lack the lower-frequency component in their output.
This invention will help simplify the whole process of oil exploration by improving the reliability of the seismic results. This invention provides a method of direct oil detection in the subsurface formations, and of mapping the location and extent of the oil-bearing rocks.
In the drawings,
At present, seismic surveys are recorded to map the subsurface structural anomalies and decisions are made based on this information. To improve the reliability of the results, further analysis is carried out based on the variations in amplitude and velocities. The major shortcoming of the current seismic imaging methods that are currently being used is that none of them can identify the presence or absence and the type of reservoir fluids with any certainty. Seismic data are acquired using surface sources and surface receivers, surface sources and downhole receivers, downhole sources and downhole receivers or any combination of the data acquisition techniques.
This invention analyzes the modifications in the frequency spectrum of the seismic input impulse in the form of additional frequencies being generated and added to the spectrum. The generation of these new frequencies takes place during the propagation of the seismic impulse through a reservoir rock that is oil-saturated or water-saturated. In response to the compressional impulse, which causes a volume change of the reservoir rock, pore fluids are squirted out of the pore interconnections and the equilibrium state of the reservoir rock gets disturbed. During the time that follows the compressional cycle of the seismic impulse, the pore fluids and the rock matrix tend to go back to their original state of equilibrium. Due to the different viscosities of the pore fluids, the time taken to recover back to the original state varies. This time of recovery can be identified as a ‘time-constant’. Higher-viscosity pore fluids like oil require a larger recovery time-constant compared to water or gas that have lower viscosity, hence a lower recovery time-constant. A larger time-constant causes the seismic pulse to stretch in time and generates lower frequencies that are not present in the original seismic input signal.
In
In
This invention introduces a new method to map and highlight the oil-bearing subsurface formations in comparison with formations that are dry or water-saturated. It uses a seismic short duration impulse as a source for 2D, 3D, or borehole seismic, to perform seismic reflection recording according to the current art. As shown in
For further analysis of the relative amplitudes of the lower frequencies in a particular seismic horizon of interest, there are many available software algorithms that are readily available and are being used as a routine in the industry today. A time or depth window, which covers the zone of interest, can be analyzed using Fourier Transform or Spectral Decomposition, which will highlight the differences in the relative amplitudes of the lower frequencies in a particular time or depth window to show the exact location and the extent of the oil-bearing subsurface rocks. Any other suitable frequency filtering methods that have been available for many years can be used to show the presence of oil-bearing formations on a seismic reflection cross-section image. This method of mapping oil-bearing subsurface rock formations is valid for all type of seismic (2D, 3D, Borehole, offshore or onshore).
Johnson, P. A.; Dynamic Elastic Nonlinearity in Earth Materials. Research Thesis for Doctorate, Paris University, December 1997.
O'Connell, R. J.; Viscoelastic Properties of Fluid-Saturated Cracked Solids. American Geophysical Union, 1977.
Winkler, K., Nur, A.; Pore Fluids and Seismic Attenuation in Rocks. American Geophysical Union, 1979.
Mavko, G. M., Nur, A.; Wave Attenuation in Partially Saturated Rocks. SEG Meeting, 1976.
Dvorkin, J.; Self-Similarity in Rock Physics.
Journal: The Leading Edge, August 2007.
Khan, T. A., Khan McGuire, S.; Elastic Nonlinearity of Reservoir Rocks—A Paradigm Shift;
Journal: CSEG Recorder, May 2005.
Khan, T. A., Khan McGuire, S.; Nonlinear seismic data can improve hydrocarbon recovery from reservoirs;
Journal: Oil & Gas Journal, June 2006.