The present invention generally relates to the field of petroleum geology and, more specifically, to a system and method for estimating mechanical properties of a subsurface rock formation backward through time.
This section is intended to introduce various aspects of the art, which may be associated with some embodiments of the present invention. This discussion is believed to assist in providing a framework to facilitate a better understanding of particular aspects of the present invention. Accordingly, it should be understood that this section should be read in this light, and not necessarily as admissions of prior art.
The prediction or estimation of the mechanical behaviors of a rock formation is of great importance to the petroleum industry. Such estimations allow, for example, for the strength of the rock formation to be modeled and the failures within the rock layers of a hydrocarbon system to be predicted. Rock failures can take a variety of forms, such as compaction (analogous to compressing a sponge), brittle fracture (i.e., creation of faults), or ductile behavior (i.e., flowage such as salt flows, or distributed deformation). A basic understanding of a typical hydrocarbon system aids in the appreciation of such failures.
Referring now to the simplified diagram shown in
As more layers of sediment are piled on, the temperature and pressure of the source rock 102 rises. If the temperature and pressure reach appropriate levels, often referred to as thermal windows, the organic material 104 within the source rock 102 is transformed into hydrocarbons 106, such as, but not limited to, natural petroleum and natural gas.
In the depicted example, adjacent to the source rock 102 is a porous, permeable layer of rock called reservoir rock 108. At least some of the hydrocarbons 106 generated from the source rock 102 have a lower density than the water in the pores of the source rock 102. As a result, the hydrocarbons 106 tend to migrate, or float, upward from the source rock 102 and into the porous reservoir rock 108.
As further depicted in
As noted above, the various layers of rock within the hydrocarbon system 100 are susceptible to failure. However, some rock failures can be thought of as “good fractures” while other are thought of as “bad fractures”. When drilling a hydrocarbon production well, it is desirable to drill the well bore in such a way as to take advantage of existing fractures or potential fractures to enhance the permeability of the reservoir and increase the flow of hydrocarbons through the reservoir to the well. Such fractures are considered “good fractures.”
Conversely, if a fracture were to occur in the seal rock 110 creating trap 114, then the hydrocarbons 106 could be released from the trap 114 and allowed to migrate to the surface 112. This would be an example of a “bad fracture”. In another example, if the source rock 102 was relatively impermeable and a fracture happened to establish a fluid connection between the reservoir rock 108 and an aquifer 116, then water would migrate from the aquifer 116 into the pores of the reservoir rock 108. The migrated water could plug the pores of the reservoir rock 108 and make it difficult, if not to impossible, for the hydrocarbons 106 to flow through its pore structure up to the well head.
The current state of the art in geology allows prediction of certain rock constitutive behaviors such as failure, and thus prediction of fracture, by using estimates of (i) the various stresses existing at different depths and locations in the formation and (ii) the strength of the rock, which generally increases with depth. Estimates of the first-mentioned values, namely stress patterns in a rock formation, are typically a function of three variables: the weight of the overburden above the formation, the degree of “push” or “pull” that the formation experiences due to tectonic forces, and the mechanical properties of the formation and adjacent rock units
One constraint to such a method of failure and fracture prediction is that it is not possible today to directly measure the values of the mechanical properties of the rock as they existed millions of years ago, as the rock was being transformed by compaction and other diagenetic processes. To address this constraint, present-day conditions of the rock are measured and by “working backwards” from that information, inferences can be drawn about the likely properties of the rock in the past.
A variety of generalized strategies for modeling ancient rock properties have been used to provide information for a fracture-prediction analysis. One approach makes use of a linear interpolation between (i) rock properties assumed to exist at the time of the rock's deposition or soon thereafter and (ii) properties measured or assumed at the present day. The interpolation is accomplished at the proportional geologic time interval that is being modeled.
Another approach is to assume that all rock properties, such as cohesion, internal angle of friction, consolidation pressure, etc., are a function of burial depth and compaction alone. In this approach, the estimation can be as simple as deriving a compaction history by numerically “stacking” layer upon layer and allowing gravity and the weight of the overlying layers to compress and densify the layer or layers of interest. This can also be calculated using a density profile with depth and computing the mean stress with which approximates the consolidation pressure.
Another approach is to measure present-day properties of the rocks of interest and to either use those values explicitly or modify those values for some time in the past. This technique could be described as an educated guesswork. Wholly ad hoc guesses have also been used.
None of these techniques provide a systematic and scientific way to estimate rock properties at some time in the past. Thus, there is a need for improvement in this field.
The present disclosure provides a system and method for mechanical property prediction of a subsurface rock formation.
One embodiment of the present disclosure is a method of estimating a mechanical property value of a rock formation comprising: generating a correlation model correlating at least one mechanical property with a porosity value and a composition value relative to a first rock type; generating an estimation of a past porosity value and a past composition value for the rock formation at a point of time in the past; and identifying in the correlation model a mechanical property value at the point of time in the past based on the estimated past porosity and past composition values.
The foregoing has broadly outlined the features of one embodiment of the present disclosure in order that the detailed description that follows may be better understood. Additional features and embodiments will also be described herein.
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the U.S. Patent and Trademark Office upon request and payment of the necessary fee.
The present invention and its advantages will be better understood by referring to the following detailed description and the attached drawings in which:
It should be noted that the figures are merely examples of several embodiments of the present invention and no limitations on the scope of the present invention are intended thereby. Further, the figures are generally not drawn to scale, but are drafted for purposes of convenience and clarity in illustrating various aspects of certain embodiments of the invention.
For the purpose of promoting an understanding of the principles of the invention, reference will now be made to the embodiments illustrated in the drawings and specific language will be used to describe the same. It will nevertheless be understood that no limitation of the scope of the invention is thereby intended. Any alterations and further modifications in the described embodiments, and any further applications of the principles of the invention as described herein are contemplated as would normally occur to one skilled in the art to which the invention relates. One embodiment of the invention is shown in great detail, although it will be apparent to those skilled in the relevant art that some features that are not relevant to the present invention may not be shown for the sake of clarity.
Persons skilled in the technical field will readily recognize that in practical applications of the disclosed methodology, it must be performed on a computer, typically a suitably programmed digital computer. Further, some portions of the detailed descriptions which follow are presented in terms of procedures, steps, logic blocks, processing and other symbolic representations of operations on data bits within a computer memory. 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 the present application, 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 or magnetic signals capable of being stored, transferred, combined, compared, and otherwise manipulated in a computer system.
It should be borne in mind, however, that all of these and similar terms are to be associated with the appropriate physical quantities and are merely convenient labels applied to these quantities. Unless specifically stated otherwise as apparent from the following discussions, it is appreciated that throughout the present application, discussions utilizing the terms such as “generating”, “identifying”, “processing”, “computing”, “calculating”, “determining”, or the like, refer to the action and processes of a computer system, or similar electronic computing device, that manipulates and transforms data represented as physical (electronic) quantities within the computer system's registers and memories into other data similarly represented as physical quantities within the computer system memories or registers or other such information storage, transmission or display devices.
The following statements are provided to aid the reader in understanding the present disclosure. While these statements are not a complete definition of term as understood to someone skilled in the relevant art(s), in general, the terms provided below have the following meanings in the context of the present disclosure.
Cement refers to fine-grained minerals (i.e., calcite, quartz, or clay) that bind together the more-coarse-grained matrix (i.e., sand or silt) of sedimentary rocks.
Chemical rock refers to rock formed by precipitation of minerals, such as from seawater. One common form of chemical rock is limestone.
Clastic rock refers to rock created from clasts. Clasts are fragments or grains that have been broken or eroded from pre-existing rock.
Compaction refers to the packing of later-deposited sediment grains into pores within an earlier-deposited layer under the weight of still-later-deposited sediment.
Composition: In the context of clastic rock, composition refers to the relative proportions of sand- to clay-sized particles. For chemical rocks, composition refers generally to the relative proportions of various precipitated minerals, such as, calcite, dolomite, ankerite, barite, etc. Composition can be determined through a variety of known techniques, such as, but not limited to, X-ray diffraction analysis of a whole-rock powder or petrographic point counting.
Diagenesis refers to physical, chemical, and/or biological changes to a given layer of sediment that occur over time as the layer is buried by the deposit of successive additional layers. Such changes may result from (i) compaction of the sediment, (ii) partial or complete dissolution of some minerals and/or organic materials in the given layer, resulting from compaction-generated increases in their solubility, and/or (iii) precipitation of new minerals into spaces left by such dissolution. The diagenesis process can be affected by heat generated by radioactive decay of elements that may be present in sediment, such as uranium, cesium, and the like.
Kerogen refers to certain organic materials found in sedimentary rock that are converted to hydrocarbons (e.g., oil and/or natural gas) when the rock's temperature reaches and remains at the required temperature for the required amount of time. See also Thermal maturity.
Lopatin-style modeling may involve, among things, (i) measuring the depth of different rock layers in a well bore and (ii) determining the ages of the various layers of rock.
The age determination may be accomplished by various techniques, such as, but not limited to, radiometric dating (i.e., carbon dating) and/or through comparisons of fossils found in the rock layers in what amounts to a forensic study of the rock layers penetrated by the well bore. Such modeling is commonly implemented via software, such as the TOUCHSTONE™ software commercially available from Geocosm, LLC in Austin, Tex., or the TEMIS software commercially available from the Institut Français du Petrole (IFP) in Rueil-Malmaison Cedex, France.
Mechanical properties of a rock that are of interest could include: Poisson's ratio, Young's modulus, angle of internal friction, and cohesion. However, other mechanical properties may be of interest for other applications.
Point counting generally involves mounting a slice of rock to a microscope slide and then polishing it down to a thickness of 30 microns. A microscope is used to count rock grains, pores, and/or other items of interest, in an area of the slice divided into a grid. The point count may be performed either by eye or by computer-controlled scanning.
Porosity generally refers to the extent to which rock is composed of spaces or voids. These spaces are commonly filled with air or water until other materials are packed into the spaces. Porosity is defined as the ratio of the volume of voids within a rock to the total volume of the rock. It can be determined through a variety of standard techniques, both laboratory and petrographically based.
Sedimentary rock may be formed by the accumulation of mineral and/or organic particles, typically carried by moving fluids such as wind or water flows. Sedimentary rock may also be formed by the precipitation of minerals from a solution. Sedimentary rock is often categorized as clastic rock or as chemical rock.
Thermal history, also known as temperature history, is typically obtained by analyzing two types of data. The first data type is the present-day temperature profile as measured for the various layers comprising a rock formation. Second is an analysis of the thermal maturity indicators for the various layers.
Thermal maturity refers to the extent to which kerogen in a rock has been converted into hydrocarbons by chemical or physical processes. Thermal maturity depends in part on the maximum temperature the rock has reached over time, and in part on the length of time elapsed at that temperature. Thermal maturity is sometimes assessed by measuring the light reflectance of a type of organic matter, known as vitrinite—akin to small lumps of coal—which increases with the maximum temperature reached and length of time at that temperature. Thermal maturity can also be assessed through a variety of other techniques, such as, but not limited to, measuring biomarker transformations, fluid inclusions, the color of spores and pollen, etc.
The disclosed methodology is a process to produce an estimate of one or more rock properties of interest at one or more times of interest.
Many known techniques combine an understanding of the age, thickness, and composition of rock layers deposited in a basin with estimates of the thermal stress experienced by the basin in order to estimate the timing and degree of compaction and cementation by minerals precipitated in the pore space, as well as the conversion of some mineral types to other mineral types. The evolution models may then be calibrated with present-day property of the rock(s) in order to make then more accurate. The present-day property data may include, but is not limited to, amount and types of cementation, thermal maturity indicators, and amounts and kinds of alteration products. The evolution model 202 provides a quantitative way to estimate how one or more of these rock properties have evolved over the lifetime of the sedimentary layers of interest.
The depicted process 200 also utilizes a correlation model 204. In one embodiment, the correlation model 204 is derived from a set of empirical calibrations that relate measured rock compositions, porosities, and mechanical properties for a particular rock type. In some embodiments, a plurality of correlation models is used. As will be appreciated by those of ordinary skill, rock porosity can be measured through any variety of known techniques, both laboratory and petrographically based.
Information from the evolution model 202 is then applied to the correlation model 204 to produce an estimate 206 of one or more rock mechanical properties of interest at one or more particular times of interest. In one embodiment, past rock composition and porosity data from the evolution model is applied to the correlation model. In some embodiments, the algorithms utilized by the correlation model 204 are directly applied to the evolution model 206 to generate, identify or estimate a historic mechanical property data set or volume. In other embodiments, porosity and composition values are determined for specific points of time in the past, thereby allowing mechanical properties to be calculated at each corresponding past time.
Another embodiment of the present disclosure is provided in
Process 300 begins by providing present day porosity data 302, present day composition data 304 and present day mechanical property data 306. Composition data 304 may comprise, but is not limited to, the amount of clay or clay-sized matrix in the rock. Data 302, 304, 306 may to pertain to sandstones, wackes, shales, or carbonates to name a few non-limiting examples.
As appreciated by those skilled in the art, data 302, 304, 306 may be obtained through a variety of known techniques and may be comprised of a single data set or separate data sets. In some embodiments, porosity data 302 may be obtained by applying standard thin-section point-counting techniques or helium porosimetry techniques to a given rock sample. In some embodiments, composition data 304 may be obtained by applying standard thin-section point-count techniques or quantitative measurement of the composition via x-ray diffraction techniques to a rock sample. It will be further appreciated by those of ordinary skill having the benefit of this disclosure that it may be efficient to simultaneously collect porosity data 302 and composition data 304 using a point-counting technique. In some embodiments, mechanical property data 306 may be measured directly using conventional techniques or ascertained from previously-stored data.
At step 308, the present day porosity data 302, the present day composition data 304, and the present day mechanical property data 306 are processed to produce a mechanical property correlation model. In one embodiment, the correlation model relates rock porosity and composition information with at least one mechanical property. In some embodiments, a correlation model can be produced for a variety of rock types, such as, but not limited to, sandstones, wackes, shales, or carbonates. A database of correlation models may also be provided. In such an embodiment, the correlation model database may comprise a wide spectrum of sedimentary lithologies most often encountered within a hydrocarbon reservoir environment.
In one embodiment, a correlation model consists of sets of related data points. In another embodiment, a correlation model consists of one or more equations describing one or more relationships between data points. While the depicted embodiment utilizes porosity, composition and mechanical property data, other embodiments may utilize fewer data types, additional forms of data, and/or different forms of data to produce the correlation model.
One approach to generating the mechanical property correlation model (step 308) is described in U.S. Patent Application Publication No. 2011/0015907 to Crawford et al. (hereafter the “'907 Publication”), which is incorporated by reference in its entirety. The techniques described in the '907 Publication are well suited for modeling sandstone rocks. For shale rocks, the correlation model may be derived using conventional techniques, such as those described by Nazmul H. Mondol, Knut Bjøorlykke, Jens Jahren, and Kaare Hoeg in “Experimental mechanical compaction of clay mineral aggregates—Changes in physical properties of mudstones during burial,” M
Carbonate rocks present different challenges for deriving a correlation model. Correlation of present-day porosity and composition to the mechanical properties of the rock can be done in accordance with known techniques. See, e.g., Elisabeth Bemer, O. Vinke and P. Longuemare, “Geomechanical Log Deduced from Porosity and Mineralogical Content,” in Oil & Gas Science and Technology—Rev. IFP Energies nouvelles, vol. 59, no. 4, pp. 405-426, July-August 2004.
As depicted, process 300 also includes providing thermal history data 310 and burial history data 312 for a given rock sample. Thermal history data 310 and burial history data 312 may be obtained using known historical modeling techniques, such as Lopatin style modeling and/or a time-temperature index (TTI) of thermal maturity. The thermal history data 310 and the burial history data 312 are then processed to generate an evolution model 314. In some embodiments, additional data may also be used to generate the evolution model. The evolution model provides a quantitative way to estimate how one or more of rock properties have evolved over the lifetime of a sedimentary layer of interest. As a result, estimates of past rock properties can be obtained (step 316) from the evolution model.
In one embodiment, a past porosity value of the rock of interest is obtained. In other embodiments, other rock property values are identified, such as, but not limited to, rock cementation or mineralogical composition values. Cementation refers to the amounts and types of cements present in a rock. Those of ordinary skill in the art will understand that for a given type of cement, an estimate of the amount and type of cement could be zero. In some embodiments, estimates of past rock properties are produced at a fixed time interval (i.e., every million years). In other embodiments, estimates of past rock properties are obtained only for a particular time, or times, of interest. For clastic rocks, the clay content of the rock (i.e., the amount of fine-grained matrix) is of particular interest in the composition modeling.
After the past rock property values are obtained, selected or identified, they are then applied to the correlation model(s) 318. In one embodiment of the disclosed methodology, past porosity and composition values identified via the evolution model are applied to the correlation model. These past values are then used to estimate a past mechanical property 320. In one embodiment, the past porosity and composition values are applied to a single correlation model for all past times. In another embodiment, past porosity and composition values corresponding to different times are applied to different correlation models. By shifting between different correlation models, one embodiment of the disclosed methodology can adjust for varying conditions (i.e., rock types) experienced during a rock's formation.
The estimated mechanical property information 320 may then be utilized in a variety of applications. Fracture sets can be formed from a variety of conditions, such as, but not limited to, weakening caused by increases in porosity or increases in clay caused by diagenesis. Therefore, in some embodiments, the disclosed methodology may be applied to a data set in order to reveal or predict non-tectonic ways in which regional fracture sets may have been formed (step 322). These past mechanical property estimates may also be used to model and predict rock formation strength (step 324).
The past mechanical property information may be further used to produce an estimate of the earth's stress conditions at a horizon of interest at the time of interest (step 326). The stress information may then be used, for example, to model the mechanical evolution of a geological structure. The historical rock mechanical property information is also helpful in predicting the occurrence, orientation and/or density of fractures that may have formed as a result of rock failure. As appreciated by those of skill in the art, steps 322, 324 and 326 are provided in
Having the benefit of the present disclosure, those of ordinary skill will readily appreciate the utility of the mechanical rock property estimates derived from the disclosed methodologies. In one example, the historical mechanical rock property estimates may be used to refine and improve evolution models of geological structures. Relevant evolution modeling may include, but is not limited to, assessing the stress history of a hydrocarbon reservoir, predicting the presence of faults in a rock formation and/or assisting in the interpretation of poorly-imaged regions within a hydrocarbon trap. In some embodiments, the improved estimates of mechanical rock property information allow workers to better determine where and how to drill a well into a particular rock formation.
The empirical correlations depicted in
From the empirical correlations depicted in
The relationships depicted on these graphs demonstrate how the strength of a rock increases systematically with decreasing porosity. In embodiments of the present disclosure, this relationship is combined with a system to predict porosity evolution through time to derive rock strength as a function of time. In general this is straightforward, except in cases where there is considerable creation of clay in the pore space of a rock. In these cases, rock strength can actually decrease through time owing to crossing the mineralogical thresholds between arenite and wacke, or wacke and shale.
normal stress relationship for two adjacent rock samples, from the time they were deposited to the present day.
An advantage of the disclosed techniques is that they provide an estimate of the earth's stress conditions, at the horizon of interest, at the time of interest. This is important, for example, when trying to model the mechanical evolution of a geological structure. It is also important in attempting to predict the occurrence, and thus the orientations and density, of fractures that may have formed as a result of rock failure.
Unless otherwise indicated, in “provide”-type operations described below, the data or other subject matter being provided may be acquired from a pre-existing store, or it may be created, measured, derived, etc., specifically for the purpose.
Those of ordinary skill having the benefit of this disclosure will appreciate that some of the operations described below need not necessarily be conducted in the precise sequence described. It is important to note that the steps depicted in
Embodiments of the present invention 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 stored in the computer. Such a computer program may be stored in a computer readable medium. A computer-readable medium includes any mechanism for storing or transmitting information in a form readable by a machine (e.g., a computer). For example, but not limited to, a computer-readable (e.g., machine-readable) medium includes a machine (e.g., 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.), and a machine (e.g., computer) readable transmission medium (electrical, optical, acoustical or other form of propagated signals (e.g., carrier waves, infrared signals, digital signals, etc.)).
Furthermore, as will be apparent to one of ordinary skill in the relevant art, the modules, features, attributes, methodologies, and other aspects of the invention can be implemented as software, hardware, firmware or any combination of the three. Of course, wherever a component of the present 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 present invention is in no way limited to implementation in any specific operating system or environment.
It should be understood that the preceding is merely a detailed description of specific embodiments of this invention and that numerous changes, modifications, and alternatives to the disclosed embodiments can be made in accordance with the disclosure here without departing from the scope of the invention. The preceding description, therefore, is not meant to limit the scope of the invention. Rather, the scope of the invention is to be determined only by the appended claims and their equivalents. It is also contemplated that structures and features embodied in the present examples can be altered, rearranged, substituted, deleted, duplicated, combined, or added to each other. The articles “the”, “a” and “an” are not necessarily limited to mean only one, but rather are inclusive and open ended so as to include, optionally, multiple such elements.
Nur et al., “Critical porosity: a key to relating physical properties to porosity in rocks,” The Leading Edge 17, 357-362 (1998).
This application claims the benefit of U.S. Provisional Patent Application 61/695,774, filed Aug. 31, 2012, entitled METHOD OF ESTIMATING ROCK MECHANICAL PROPERTIES, the entirety of which is incorporated by reference herein.
Number | Date | Country | |
---|---|---|---|
61695774 | Aug 2012 | US |