This invention relates generally to the field of interpretation of reservoir monitoring using time-lapse seismic measurements, and more specifically, to a method of predicting fluid and/or pressure change from 4D seismic data.
This section is intended to introduce various aspects of the art, which may be associated with embodiments of the disclosed techniques and methodologies. A list of references is provided at the end of this section and may be referred to hereinafter. This discussion, including the references, is believed to assist in providing a framework to facilitate a better understanding of particular aspects of the disclosure. Accordingly, this section should be read in this light and not necessarily as admissions of prior art.
Reservoir surveillance during hydrocarbon reservoir production is a key to meeting goals of reduced operating costs and maximized recovery. Differences between actual and predicted performance are typically used to update the reservoir's geological model and to revise its depletion strategy. The changes in reservoir fluid saturation, pressure and temperature that occur during production also induce changes in the reservoir acoustic properties of rocks that under favorable conditions may be detected by seismic methods. One way to observe changes in properties in a hydrocarbon reservoir is to compare seismic returns measured at different times. The resultant gathered data is known as 4D seismic data.
For porous siliciclastic sands, seismic response is generally strong enough to detect fluid movement (Johnston et al., many others). Using different forms of approximations of the reflectivity equation and a rock physics model, a set of coefficients can be estimated to make a combination of near and far difference amplitudes, which in turn may infer the saturation and pressure change (Tura, et al, 1998; Landro, 2001; Lumley, et al, 2003, Ribeiro & MacBeth, 2004, Angelov et al, 2004). This type of formulation works well when physical properties of the reservoir, such as porosity (phi) and volume of shale (vsh), are relatively constant.
Seismic velocities change as hydrocarbon saturation change caused by production. The amount of the velocity change for a given saturation change depends on a number of factors, i.e. nature of the rock frame; fluid property; pore pressure, temperature of the reservoir, etc. These dependencies are well understood and are implemented in currently available commercial software for fluid substitutions to model/explain the fluid effect on seismic data. The scales at which the two fluids (e.g. oil and water) are mixed have an impact on the velocity change from experiment data (Knight and Nolen-Hoeksema, 1990, Endres and Knight, 1989) and theoretical derivations (Biot 1956, Mavko and Nur 1975, Stoll, 1989, Dvorkin and Nur, 1993). The mixing pattern of the two fluids is in turn influenced by the permeability at a given frequency (e.g. seismic frequency) (Batzle et al. 2006). For low permeability rocks, fluids can not move easily as in the higher permeability rock. Fluid moving through low permeability rock therefore has a higher chance of showing unpredictable or patchy behavior at frequency ranges associated with seismic surveys.
Time shift is the difference in two-way seismic travel times that are observed when analyzing seismic surveys conducted at different times. Time shifts can be attributed to two sources: 1) pore-fluid property changes that alter the velocity at which seismic signals pass through a layer or interval, and 2) changes in seismic velocity and layer thickness that occur both inside and outside of the layer or overall reservoir because of reservoir compaction and stress-strain redistribution in the surrounding formations. For a layer with thickness, z, the change in relative seismic travel time is (Landro and Stammeijer, 2004)
Δt/t=Δz/z−Δv/v (1)
where t represents the two-way travel time across the layer and v is the velocity of the layer. Hatchell and Bourne (2005) has used time shifts to estimate the compaction of compacting reservoirs and layers within reservoirs by implicitly assuming the velocity change due to fluid is a minor component. However, compaction (as expressed by Δz) does not always contribute significantly to the time shifts, especially when reservoir rock is stiff or pressure is well maintained. In this case, the main contributing factor to an observed time shift will be the velocity change from either fluid change (such as saturation levels) or pressure change.
Separating the contributions from pressure change or saturation change is important to an effective 4D seismic model. He et al. (1998) has proposed to invert the 4D seismic data to impedance volumes and from there pressure and saturation change are inferred. Chu and Grant (2007) proposed a method to invert saturation and pressure directly from 4D seismic data by integrating well information, reserve simulation and 4D seismic data. However, when 4D signals are weak, such as in hard rock formations, a more model based methodology is desired.
An additional concern is that when a reservoir is thick, it may span a few cycles of seismic wave events, so that side-lobe energy may generate apparent difference events that can appear as real reservoir differences. This side-lobe energy complicates interpretation of multi-cycle hydrocarbon reservoirs where there is interference between reflectors. One proposed method to overcome this is to invert saturation and pressures changes at well locations using rock physics models and then populate the saturation and pressure changes in-between wells by statistical calibration methods. However, in hard rocks, such as carbonate rocks, seismic signals may not provide an adequate signal return to observe changes in hydrocarbon reservoir property values over time. What is needed is a method to optimally evaluate dynamic reservoir properties, such as saturation changes and pressure changes, under weak 4D signal scenarios generally encountered when dealing with hard rock formations.
The foregoing discussion of need in the art is intended to be representative rather than exhaustive. A technology addressing one or more such needs, or some other related shortcoming in the field, would benefit drilling and reservoir development planning, for example, providing decisions or plans for developing a reservoir more effectively and more profitably.
Reference material which may be relevant to the invention, and which may be referred to herein, include Angelov, P. et al, “Pore Pressure and Water Saturation Variations—Modification of Landro's AVO Approach,” Expanded Abstracts: 74th Annual Meeting of the SEG (2004); Hatchell, P. J. et al., Measuring Reservoir Compaction Using Time-Lapse Timeshifts,” Expanded Abstracts: 75th Annual Meeting of the SEG (2005); Landro, M., “Discrimination Between Pressure and Fluid Saturation Changes From Time-Lapse Seismic Data,” Geophysics, 66, 836-844 (2001); Lumley, D. E. et al., “Estimation of Reservoir Pressure and Saturations by Crossplot Inversion of 4D Seismic (2003); Ribeiro, C. et al., “A Petroelastic-Based Approach to Pressure and Saturation Estimation Using 4D Seismic,” Expanded Abstracts: 74th Annual Meeting of the SEG (2004); Tura, A. et al., “Subsurface Fluid Flow Properties From Time-Lapse Elastic Wave Reflection Data,” Proceedings of SPIE, Mathematical Methods in Geophysical Imaging V, V. 3453, 125-138 (1998); Knight, R. et al., “A Laboratory Study of the Dependence of Elastic Wave Velocities on Pore Scale Fluid Distribution,” Geopys. Res. Lett., 14, 1529-1532 (1990); Endres, A. L. et al., “A Theoretical Treatment of the Effect of Microscopic Fluid Distribution on the Dielectric Properties of Partially Saturated Rocks,” Geophys. Prospecting, 40, 307-324 (1992); Biot, M. A., “Theory of Propagation of Elastic Waves in Fluid Saturated Porous Solid, I. Low Frequency Range and II. Higher-Frequency Range,” J. Acoust. Soc. Am., 28, 168-191 (1956); Mavko, G. et al., “Melt Squirt in the Asthenosphere,” J. Geophys. Res., 80, 1444-1448 (1975); Stoll, R. D. “Sediment Acoustics,” Springer Verlag, Berlin, 154, (1989); Dvorkin, J. et al., “Dynamic Poroelasticity: A Unified Model With the Squirt and the Biot Mechanisms,” Geophys., 58, 524-533, (1993); Batzle, et al., “Fluid Mobility and Frequency-Dependent Seismic Velocity—Direct Measurements,” Geophysics, 71, N1-N9, (2006); Landro, M. et al., “Quantitative Estimation of Compaction of Velocity Changes Using 4D Impedance and Traveltime Changes,” Geophysics, 69, 949-957, (2004); U.S. Pat. No. 5,798,982; PCT Patent Application No. PCT/US08/03830, Inversion of 4D Seismic Data, with inventors Chu, D. et al.
A method is disclosed including analyzing seismic data related to a hydrocarbon reservoir containing hard rock. The seismic data includes a first set of seismic data taken at a first time and a second set of seismic data taken at a second, subsequent time. An interval in the hydrocarbon reservoir composed substantially of hard rock is identified. The interval has a permeability level associated therewith. 4D seismic attributes for the interval are calculated. A first set of rock physics relationships are applied to seismic data related to the interval when the permeability level associated therewith is in a first permeability range. A second set of rock physics relationships is applied to seismic data related to the interval when the permeability level associated therewith is in a second permeability range. The first and second ranges are non-overlapping. A fluid saturation change or a pressure change of the interval are inferred based on the rock physics relationships from 4D seismic attributes. The inferred fluid saturation change or pressure change of the interval is outputted.
According to aspects of the disclosed techniques and methodologies, the interval may be a first interval. An additional interval or intervals composed substantially of hard rock in the reservoir may be defined. Each additional interval or intervals has a permeability level associated therewith. Steps of calculating, applying, inferring, and outputting may be repeated for the additional interval or intervals. According to another aspect, an output of a reservoir simulation model may be verified using the fluid saturation change or the pressure change of the interval. The reservoir simulation model may be revised so that prediction from the simulation model is consistent with the fluid saturation change and the pressure change calculated from 4D seismic and selected rock physical relationships. The corrected output or the verified output of the reservoir simulation model may be outputted. Hydrocarbon extraction activity may be conducted based on the output of at least one of the fluid saturation change and the pressure change of the interval computed from the outputs of the one of the rock physics relationships and the 4D seismic attributes for the interval. The first set of rock physics relationships may include an equation to determine a bulk modulus of fluid in the interval using a harmonic (e.g., Voigt) average or an arithmetic (e.g., Reuss) average. The components of the fluid may include at least two of gas, oil, and water. Seismic data may include one or more of co-processing the seismic data to maintain maximum repeatability, preserving time shifts between the second and first sets of seismic data, making time alignments to the seismic data, cross-equalizing the seismic data, and retaining derivative volumes. At least one of the 4D seismic attributes may describe a change over time of a characteristic of the interval. One of the 4D seismic attributes may be an amplitude of a returning seismic signal, or a change over time of such an amplitude. One of the 4D seismic attributes may be a time shift. The second set of seismic data may relate to the hydrocarbon reservoir at least one additional time.
In another aspect, a method is provided including defining a plurality of non-overlapping subsurface intervals in a hydrocarbon reservoir, each subsurface interval being composed primarily of hard rock. 4D seismic data related to each of the subsurface intervals are analyzed separately. The 4D seismic data include base seismic data and monitor seismic data. Seismic attributes for each of the subsurface intervals are calculated. Rock physics relationships are applied to the 4D seismic data related to each of the subsurface intervals based upon a permeability level for each subsurface interval. A fluid saturation change of each of the subsurface intervals is inferred based on the rock physics relationships and seismic attributes of each of the subsurface intervals. Outputs of the hydrocarbon reservoir simulation model are verified using the inferred fluid saturation change of each of the subsurface intervals. The hydrocarbon reservoir simulation model is modified according to the inferred fluid saturation change of each of the subsurface intervals. The verified model or the modified model is outputted.
Aspects of the disclosed techniques and methodologies may also include conducting hydrocarbon extraction activity based on the output from the modified reservoir simulation model. The rock physics relationships may include an equation to determine a bulk modulus of fluid in the interval using a harmonic average of bulk moduli of components in the fluid when permeability in a subsurface interval is less than a first predetermined amount, and an arithmetic average of bulk moduli of the fluid components when the permeability in a subsurface interval is greater than a second predetermined amount. The seismic attributes may include an amplitude of a returning seismic signal, a change over time of an amplitude of a returning seismic signal, and a time shift.
Also disclosed is a method of extracting hydrocarbons from a subsurface reservoir. According to the method, Seismic data relating to a hydrocarbon reservoir containing hard rock is analyzed. The seismic data includes a first set of seismic data relating to the hydrocarbon reservoir at a first time and a second set of seismic data relating to the hydrocarbon reservoir at a second time. An interval composed substantially of hard rock in the hydrocarbon reservoir is defined. The interval has a permeability value associated therewith. 4D seismic attributes for the interval are calculated. A first set of rock physics relationships are applied to seismic data related to the interval when the permeability associated therewith is in a first range. A second set of rock physics relationships is applied to seismic data related to the interval when the permeability associated therewith is in a second range. The first and second ranges are non-overlapping. A fluid saturation change or a pressure change of the interval is inferred based on outputs of the first or second sets of rock physics relationships and the calculated 4D seismic attributes for the interval. An output of a reservoir simulation model is verified using the fluid saturation change or the pressure change of the interval. The reservoir simulation model is revised so that a prediction from the simulation model is consistent with the fluid saturation change and the pressure change calculated from 4D seismic and selected rock physical relationships. The revised output of the reservoir simulation model or the verified output of the reservoir simulation model is outputted. Hydrocarbons are extracted from the subsurface reservoir based on the outputted revised output or the verified output.
Further disclosed is a computer program product having computer executable logic recorded on a tangible computer readable medium. The computer program product includes: code for defining a subsurface interval as being composed primarily of hard rock; code for analyzing 4D seismic data related to the subsurface interval, the 4D seismic data including base seismic data and monitor seismic data; code for calculating seismic attributes for the subsurface interval; code for applying rock physics relationships to the 4D seismic data related to the subsurface interval based upon a permeability level of the subsurface interval; code for inferring a fluid saturation change of the subsurface interval based on at least one of the rock physics relationships and the seismic attributes of the subsurface interval; code for verifying an output of a hydrocarbon reservoir simulation model using the inferred fluid saturation change; code for correcting the hydrocarbon reservoir simulation model when said output is inconsistent with the inferred fluid saturation change; and code for outputting the verified output or the corrected output.
The foregoing and other advantages may become apparent upon reviewing the following detailed description and drawings of non-limiting examples of embodiments in which
To the extent that the following detailed description is specific to a particular embodiment or a particular use of the disclosed techniques, this is intended to be illustrative only, and is not to be construed as limiting the scope of the invention. On the contrary, it is intended to cover all alternatives, modifications and equivalents that may be included within the spirit and scope of the invention, as defined by the appended claims.
Some portions of the detailed description which follows are presented in terms of procedures, steps, logic blocks, processing and other symbolic representations of operations on data bits within a memory in a computing system or a computing device. These descriptions and representations are the means used by those skilled in the data processing arts to most effectively convey the substance of their work to others skilled in the art. In this detailed description, a procedure, step, logic block, process, or the like, is conceived to be a self-consistent sequence of steps or instructions leading to a desired result. The steps are those requiring physical manipulations of physical quantities. Usually, although not necessarily, these quantities take the form of electrical, magnetic, or optical signals capable of being stored, transferred, combined, compared, and otherwise manipulated. It has proven convenient at times, principally for reasons of common usage, to refer to these signals as bits, values, elements, symbols, characters, terms, numbers, or the like.
Unless specifically stated otherwise as apparent from the following discussions, terms such as “analyzing”, “identifying”, “calculating”, “applying”, “inferring”, “outputting”, “defining”, “completing”, “verifying”, “revising”, “co-processing”, “preserving”, “making”, “cross-equalizing”, “retaining”, “modifying”, “revising”, or the like, may refer to the action and processes of a computer system, or other electronic device, that transforms data represented as physical (electronic, magnetic, or optical) quantities within some electrical device's storage into other data similarly represented as physical quantities within the storage, or in transmission or display devices. These and similar terms are to be associated with the appropriate physical quantities and are merely convenient labels applied to these quantities.
Embodiments disclosed herein also relate to an apparatus for performing the operations herein. This apparatus may be specially constructed for the required purposes, or it may comprise a general-purpose computer selectively activated or reconfigured by a computer program or code stored in the computer. Such a computer program or code may be stored or encoded in a computer readable medium or implemented over some type of transmission medium. A computer-readable medium includes any medium or mechanism for storing or transmitting information in a form readable by a machine, such as a computer (‘machine’ and ‘computer’ are used synonymously herein). As a non-limiting example, a computer-readable medium may include a computer-readable storage medium (e.g., read only memory (“ROM”), random access memory (“RAM”), magnetic disk storage media, optical storage media, flash memory devices, etc.). A transmission medium may be twisted wire pairs, coaxial cable, optical fiber, or some other suitable transmission medium, for transmitting signals such as electrical, optical, acoustical or other form of propagated signals (e.g., carrier waves, infrared signals, digital signals, etc.)).
Furthermore, modules, features, attributes, methodologies, and other aspects can be implemented as software, hardware, firmware or any combination thereof. Wherever a component of the invention is implemented as software, the component can be implemented as a standalone program, as part of a larger program, as a plurality of separate programs, as a statically or dynamically linked library, as a kernel loadable module, as a device driver, and/or in every and any other way known now or in the future to those of skill in the art of computer programming. Additionally, the invention is not limited to implementation in any specific operating system or environment.
Various terms as used herein are defined below. To the extent a term used in a claim is not defined below, it should be given the broadest possible definition persons in the pertinent art have given that term as reflected in at least one printed publication or issued patent.
As used herein, “4D” refers to time-lapse seismic surveying that takes into account 4-dimensional seismic returns, wherein time between acquisitions represents a fourth data dimension. The three other dimensions relate to the spatial characteristics of the earth formation, two being horizontal length dimensions, and the third relating to depth in the earth formation, which can be represented by a length coordinate, or by a time coordinate such as the two-way travel time of a seismic wave from surface to a certain depth and back. A 4D seismic survey (or 4D seismic) models alterations in geological or geophysical characteristics of a given interval of the earth by comparing data corresponding to seismic surveys carried out at different times.
As used herein, “and/or” placed between a first entity and a second entity means one of (1) the first entity, (2) the second entity, and (3) the first entity and the second entity. Multiple elements listed with “and/or” should be construed in the same fashion, i.e., “one or more” of the elements so conjoined.
As used herein, “base seismic survey” and “monitor seismic survey” designate seismic surveys of the reservoir taken at different times. The assumption is that the base survey is carried out earlier in time than the monitor survey.
As used herein, “displaying” includes a direct act that causes displaying, as well as any indirect act that facilitates displaying. Indirect acts include providing software to an end user, maintaining a website through which a user is enabled to affect a display, hyperlinking to such a website, or cooperating or partnering with an entity who performs such direct or indirect acts. Thus, a first party may operate alone or in cooperation with a third party vendor to enable the reference signal to be generated on a display device. The display device may include any device suitable for displaying the reference image, such as without limitation a CRT monitor, a LCD monitor, a plasma device, a flat panel device, or printer. The display device may include a device which has been calibrated through the use of any conventional software intended to be used in evaluating, correcting, and/or improving display results (e.g., a color monitor that has been adjusted using monitor calibration software). Rather than (or in addition to) displaying the reference image on a display device, a method, consistent with the invention, may include providing a reference image to a subject. “Providing a reference image” may include creating or distributing the reference image to the subject by physical, telephonic, or electronic delivery, providing access over a network to the reference, or creating or distributing software to the subject configured to run on the subject's workstation or computer including the reference image. In one example, the providing of the reference image could involve enabling the subject to obtain the reference image in hard copy form via a printer. For example, information, software, and/or instructions could be transmitted (e.g., electronically or physically via a data storage device or hard copy) and/or otherwise made available (e.g., via a network) in order to facilitate the subject using a printer to print a hard copy form of reference image. In such an example, the printer may be a printer which has been calibrated through the use of any conventional software intended to be used in evaluating, correcting, and/or improving printing results (e.g., a color printer that has been adjusted using color correction software).
As used herein, “hard rock” means any stiff reservoir rock characterized by high velocity and high density. In other words, hard rock means rocks whose acoustic response is controlled primarily by the rock frame and is relatively insensitive to fluid saturation change and pressure change due to hydrocarbon production. Non-limiting examples of hard rock include various carbonate or carbonate-type rocks and deep siliciclastic rocks.
As used herein, “hydrocarbon reservoirs” include reservoirs containing any hydrocarbon substance, including for example one or more than one of any of the following: oil (often referred to as petroleum), natural gas, gas condensate, tar and bitumen. The following detailed description of various embodiments is presented with primary reference to subsurface oil reservoirs, but the principles discussed apply also to situations involving reservoirs containing other hydrocarbon materials, either in addition to oil or instead of oil.
As used herein, “impedance” is the product of seismic velocity and the density of a rock formation. Impedance, also called acoustic impedance, will typically vary among different rock layers. Impedance associated with compression waves, or P-waves, is called P-impedance. Impedance associated with shear waves, or S-waves, is called S-impedance.
As used herein, “interval” is a geologic region or layer.
As used herein, “machine-readable medium” refers to a medium that participates in directly or indirectly providing signals, instructions and/or data. A machine-readable medium may take forms, including, but not limited to, non-volatile media (e.g. ROM, disk) and volatile media (RAM). Common forms of a machine-readable medium include, but are not limited to, a floppy disk, a flexible disk, a hard disk, a magnetic tape, other magnetic medium, a CD-ROM, other optical medium, punch cards, paper tape, other physical medium with patterns of holes, a RAM, a ROM, an EPROM, a FLASH-EPROM, or other memory chip or card, a memory stick, and other media from which a computer, a processor or other electronic device can read.
As used herein, “permeability” is a measure of the ability of a rock formation to transmit fluids. It is related directly to the connections (pore throats) between the spaces (pore spaces) of rock grains of a formation. These connections allow the fluids to move through the formation.
As used herein, a “permeability level” may be a single permeability value or a range of permeability values.
As used herein, “rock” includes various geological materials that may be encountered during drilling operations, e.g., salt, clay, shale, sand and the like, in addition to those materials more formally classified as “rocks.”
As used herein, “rock physics model” relates petrophysical and production-related properties to the bulk elastic properties of a rock formation. Examples of petrophysical and production-related properties may include, but are not limited to, porosity, pore geometry, pore connectivity, volume of shale or clay, estimated overburden stress or related data, pore pressure, fluid type and content, fluid and/or gas saturation, clay content, mineralogy, temperature, and anisotropy.
As used herein, “seismic data” refers to information collected from a seismic survey. Seismic data is typically processed using procedures such as filtering, multiple removal, muting, stacking, and migration.
As used herein, “seismic survey” means a program for mapping geologic structure by observation of seismic waves, especially by creating seismic waves with artificial sources and observing the arrival times and amplitudes of the waves reflected from acoustic-impedance contrasts (representing rock formation interfaces, for example) or refracted using techniques known in the art.
As used herein, “siliciclastic rock” means rock composed of fragments or portions of silicates. For example, sandstone is a siliciclastic rock.
Different rock physics models may be employed according to the permeability of a hard rock formation to predict sensitivity of seismic attributes to saturation or pressure change before seismic attributes are used to predict saturation or pressure changes. Specifically, a feature of the disclosed methodologies is predicting where pressure changes or saturation changes dominate 4D signal returns, and using appropriate 4D signal attributes to infer prospective changes. Interpreting weak 4D signal returns is improved by using permeability-dependent rock physics models. The seismic attributes for detecting saturation or pressure change are optimized. Detection limits for a given 4D signal return can also be identified.
Since fluid response to seismic signals differs with different permeability and rock type, another aspect of the invention is dividing hard rock formations into different intervals according to relative permeability levels, and then extracting seismic attributes for each interval. The seismic attributes may include various descriptions of amplitudes, such as maximum/minimum/average amplitudes of seismic returns corresponding to the interval, etc. Attributes can be calculated separately using a base seismic survey and a monitor seismic survey, and the differences can be used to create a 4D attribute. For example, one attribute that may be used is the maximum amplitude (Amax) of the returns from a seismic survey. The difference between the returns (diffAmax) from the monitor seismic survey and the base seismic survey is a useful 4D attribute. Amax, as with other amplitude based attributes, is a measure of impedance, which in turn is a function of velocity and rock density, and can be used to demonstrate viability of aspects of disclosed techniques.
Another type of attribute that may be calculated from 4D seismic data is a time shift. A time shift variable ΔT caused by fluid saturation change and pressure change in a rock formation can be defined as follows:
ΔT=dt2−dt1 (2)
In Equation 2, dt1 and dt2 are time shift variables at different levels of the reservoir. The time shift variable dt2 represents the difference in time between the seismic returns representing a geologic feature at the base of a divided hydrocarbon reservoir interval as measured at the monitor seismic survey and the base seismic survey. The time shift variable dt1 represents the difference in time between the seismic returns representing a geologic feature at the top of the hydrocarbon reservoir interval as measured at the monitor seismic survey and the base seismic survey. ΔT is an approximate measure of velocity change of seismic signals within a hydrocarbon reservoir of interest.
Reservoir simulation may not always capture accurate fluid movement through a reservoir. However, reservoir models do have a certain level of accuracy after history matching to many wells in a mature hydrocarbon reservoir. If 4D seismic data and the seismic attributes derived from them accurately reflect the fluid saturation and/or pressure change in the reservoir, the data and attributes should correlate with reservoir simulation outputs to a certain degree.
It will be understood from the preceding description, including the Background section, that the present inventive method is not limited to rocks of at least a certain hardness. Instead, the significance of rock hardness to the present invention is that the value of the invention is greatest when applied to “hard” rocks because other methods have difficulties with the weak 4D seismic data characteristic of hard rocks.
Seismic modeling is used to guide the interpretation of attributes in response to pressure and saturation changes for analyzed intervals with different permeability levels.
To model the effect of pressure on seismic wave velocity within a hard rock formation, the effect of pressure on density is ignored and empirical relationships for P-wave velocity Vp and S-wave velocity Vs are used. For Vp in terms of pressure P,
Here, Vp,∞ is a high-pressure limiting velocity and is correlated with porosity and mineralogy; A is the fractional velocity change for pressure change from 0 to high pressure; A correlates with porosity; Pp,0 is a pressure scalar that is roughly constant for hard rock, i.e., independent of porosity. Similarly, the empirical relationship for Vs can be expressed as
To model the effect of saturation, the Gassmann fluid substitution equation is used to determine elastic parameters for a given saturation state with gas saturation Sg, oil saturation So and water saturation Sw. For a given porosity φ, an overall density of a rock formation can be calculated through
ρ=(1−φ)ρs+φρf and ρf=ρoSo+ρgSg+ρwSw (5)
where ρs is the density for solid rock, ρf is the density for mixed fluid, ρo is the density of oil, ρg is the density of gas and ρw is density of water. The bulk modulus K for a rock formation is calculated by
where Kdry is the bulk modulus of the dry rock frame, i.e., when the rock formation has no fluid in its pore space; Kdry is called the dry bulk modulus; Ks is the bulk modulus of the mineral material that makes up the rock; Kf is the bulk modulus of the fluid in the pore space of the rock formation. To obtain Kf for a given reservoir fluid, it is assumed that the saturation levels of the various fluid components of oil, gas and water (So, Sg and Sw, respectively) are known, and that the bulk moduli of oil, gas, and water (Kg, Ko and Kw, respectively) are also known. According to the invention, the bulk modulus for the reservoir fluid mixture is determined by the level of permeability of the rock formation. If the permeability level of the rock formation is high, a fluid mixing scheme known as Reuss averaging is used. With Reuss averaging, the bulk modulus of the mixed fluid Kf is given by a harmonic average of the component bulk moduli (Kg, Ko and Kw):
If the permeability level of the rock formation is low, a fluid mixing scheme known as the Voigt averaging is used. With Voigt averaging, the bulk modulus for the mixed fluid is given by an arithmetic average of the component bulk moduli:
K
f
=S
g
K
g
+S
o
K
o
+S
w
K
w (8)
The statistically meaningful correlations between 4D seismic attributes (e.g., ΔT and diffAmax) and pressure/saturation changes according to the invention indicate that the 4D seismic attributes may be used to infer pressure and/or saturation changes in a hydrocarbon reservoir. The inferred pressure and/or saturation changes may be used to verify or modify reservoir simulation models. Such verified and/or modified outputs may be advantageously used to predict changes in water/oil saturation in a rock formation. Accurate changes in saturation may then be used for various hydrocarbon extraction activities such as planning the location and timing of new wells, managing production from existing wells, predicting production lifetimes of wells or hydrocarbon reservoirs at various extraction rates, and other similar activities.
Example methods may be better appreciated with reference to flow diagrams. While for purposes of simplicity of explanation, the illustrated methodologies are shown and described as a series of blocks, it is to be appreciated that the methodologies are not limited by the order of the blocks, as some blocks can occur in different orders and/or concurrently with other blocks from that shown and described. Moreover, less than all the illustrated blocks may be required to implement an example methodology. Blocks may be combined or separated into multiple components. Furthermore, additional and/or alternative methodologies can employ additional blocks not shown herein. While the figures illustrate various actions occurring serially, it is to be appreciated that various actions could occur in series, substantially in parallel, and/or at substantially different points in time.
In one embodiment, the input data are stored in disk storage device 86. The system computer 82 may retrieve the appropriate data from the disk storage device 86 to perform the model development and reservoir performance prediction according to program instructions that correspond to the methods described herein. The program instructions may be written in a computer programming language, such as C++, Java and the like. The program instructions may be stored in a computer-readable memory, such as program disk storage device 88. System computer 82 presents output primarily onto a text/graphics display 90, or alternatively to a printer 92. The system computer 82 may store the results of the methods described above on disk storage 84, for later use and further analysis. A keyboard 94 and a pointing device (e.g., a mouse, trackball, or the like) 96 may be provided with the system computer 82 to enable interactive operation. The system computer 82 may be located at a data center remote from the reservoir. Additionally, while the description above is in the context of computer-executable instructions that may run on one or more computers, those skilled in the art will recognize that the subject matter as claimed also can be implemented in combination with other program modules and/or as a combination of hardware and software.
The disclosed embodiments and methodologies may be susceptible to various modifications and alternative forms and have been shown only by way of example. The disclosed embodiments and methodologies are not intended to be limited to the particular embodiments disclosed herein, but includes all alternatives, modifications, and equivalents falling within the spirit and scope of the appended claims.
This application claims the benefit of U.S. Provisional Patent Application Ser. No. 61/166,546 filed 3 Apr. 2009, the disclosures of which are hereby incorporated herein in their entirety by reference.
Number | Date | Country | |
---|---|---|---|
61166546 | Apr 2009 | US |