CA2869825C - Three-dimensional multi-modal core and geological modeling for optimal field development - Google Patents
Three-dimensional multi-modal core and geological modeling for optimal field development Download PDFInfo
- Publication number
- CA2869825C CA2869825C CA2869825A CA2869825A CA2869825C CA 2869825 C CA2869825 C CA 2869825C CA 2869825 A CA2869825 A CA 2869825A CA 2869825 A CA2869825 A CA 2869825A CA 2869825 C CA2869825 C CA 2869825C
- Authority
- CA
- Canada
- Prior art keywords
- reservoir
- model
- subsurface
- well
- lithofacies
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011161 development Methods 0.000 title claims abstract description 16
- 208000035126 Facies Diseases 0.000 claims abstract description 42
- 239000011435 rock Substances 0.000 claims abstract description 35
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000004458 analytical method Methods 0.000 claims abstract description 21
- 238000012545 processing Methods 0.000 claims description 79
- 230000015572 biosynthetic process Effects 0.000 claims description 33
- 238000005755 formation reaction Methods 0.000 claims description 33
- 230000003068 static effect Effects 0.000 claims description 23
- 238000005259 measurement Methods 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 238000013500 data storage Methods 0.000 claims description 8
- 230000006870 function Effects 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 3
- 239000011162 core material Substances 0.000 abstract description 59
- 239000011148 porous material Substances 0.000 abstract description 30
- 238000011084 recovery Methods 0.000 abstract description 18
- 229930195733 hydrocarbon Natural products 0.000 abstract description 17
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 17
- 239000004215 Carbon black (E152) Substances 0.000 abstract description 11
- 239000012530 fluid Substances 0.000 abstract description 4
- 210000000746 body region Anatomy 0.000 abstract description 2
- 230000035699 permeability Effects 0.000 description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 235000019738 Limestone Nutrition 0.000 description 4
- 239000006028 limestone Substances 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 4
- 239000000203 mixture Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- RXZBMPWDPOLZGW-HEWSMUCTSA-N (Z)-roxithromycin Chemical group O([C@@H]1[C@@H](C)C(=O)O[C@@H]([C@@]([C@H](O)[C@@H](C)C(=N\OCOCCOC)/[C@H](C)C[C@@](C)(O)[C@H](O[C@H]2[C@@H]([C@H](C[C@@H](C)O2)N(C)C)O)[C@H]1C)(C)O)CC)[C@H]1C[C@@](C)(OC)[C@@H](O)[C@H](C)O1 RXZBMPWDPOLZGW-HEWSMUCTSA-N 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000002902 bimodal effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000011545 laboratory measurement Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 239000012267 brine Substances 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 1
- 238000012066 statistical methodology Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
- Earth Drilling (AREA)
Abstract
A roadmap for a field development strategy for optimal recovery is provided in a high quality 3D geological model. This geological model combines geological attributes, pore and rock properties for an optimum 3D representation of the reservoir thousands of feet beneath the surface. The model is based on the pertinent geological facies, derived from well core description and detailed studies of rock, as well as fluid and pore properties (Full Pore System) obtained from laboratory analyses of core material and well log data. These data differentiate various important pore throat and pore body regions and relationships, i.e., macroporosity and microporosity. Understanding hydrocarbon volumes in the various pore type groups and then establishing proper recovery techniques through focused laboratory studies yields a field development strategy that can significantly increase hydrocarbon recovery from a reservoir.
Description
THREE-DIMENSIONAL MULTI-MODAL CORE AND GEOLOGICAL MODELING
FOR OPTIMAL FIELD DEVELOPMENT
Cross Reference To Related Applications:
[0001] This application claims priority from U.S. Provisional Application No. 61/640,790, filed May 1, 2012.
BACKGROUND OF THE INVENTION
1. Field of the Invention
FOR OPTIMAL FIELD DEVELOPMENT
Cross Reference To Related Applications:
[0001] This application claims priority from U.S. Provisional Application No. 61/640,790, filed May 1, 2012.
BACKGROUND OF THE INVENTION
1. Field of the Invention
[0002] The present invention relates to computerized simulation of petrophysical structure of hydrocarbon reservoirs in the earth, and in particular to determination of petrophysical models of the reservoir based on core samples from wells in the reservoir and other field data regarding the reservoir.
2. Description of the Related Art
2. Description of the Related Art
[0003] In the oil and gas industries, the development of underground hydrocarbon reservoirs typically includes development and analysis of computer simulation models. These underground hydrocarbon reservoirs are typically complex rock formations which contain both a petroleum fluid mixture and water. One form of such reservoir model has been what is known as reservoir attributes models. For these models, the reservoir has been organized into a 3-dimensional grid of cells.
[0004] Oil and gas companies have also come to depend on geological models as an important tool to enhance the ability to exploit a petroleum reserve.
Geological models of
Geological models of
5 reservoirs and oil/gas fields have become increasingly large and complex. In geological models, the reservoir is also organized into a 3-dimensional grid of a number of individual cells.
100051 Seismic data with increasing accuracy has permitted the cells in reservoir models of both types to be on the order of 25 meters areal (x and y axis) intervals.
For what are known as giant reservoirs, the number of cells are the least hundreds of millions, and reservoirs of what is known as giga-cell size (a billion cells or more) are encountered.
100061 Days or weeks of computer processing have are usually been spent on 3D
model attributes computations. However, the guiding geological facies were not well digitally integrated into the model. Billions of cell attributes calculations could be wasted due to inappropriate algorithms from the initial model guide.
100071 When characterizing and developing a reservoir field, a 3D
geological model of the reservoir covering the entire 3D reservoir has been required to be built to provide an accurate model for reservoir planning. Accurate indications of lithofacies are an essential input in a 3D geological modeling system. Lithofacies is based on data from wells and establishes as a basis to predict reservoir properties in cells with no well data. One of the sets of data available has been what is known are formation rock descriptions or characterizations which are made by analysts based on observations made from well core samples.
10008] Traditionally, well core data has been assembled and analyzed based on measurements and observations taken from well core samples obtained from formation rock adjacent well boreholes at a number of depths of interest in a well. So far as is known, the core description has been hand drawn and was thereafter available for use in the form of a paper copy or at best a scanned graphical image of the hand drawing. The important geological information for a reservoir model came in a wide variety of forma and formats.
This information included numerical information, text strings, graphical images, relationships and drawings. These have posed a challenge to incorporation into the reservoir model in a form that could be consistently manipulated with other data.
[0009] In addition, manual preparation of these types of well description data presented a tedious challenge in drawing standardized patterns of mineral composition, texture, sedimentary structures based on well core samples repeatedly throughout the length of the well bore. Correlating core description data among different wells which had been described by different people with different scale and size hand drawings was also often problematic.
An accurate representation of the exact shape and size of formation rock features of interest in the drawing for tens of occurrences and over several wells in a reservoir was difficult to obtain in a reliable and representative manner.
[0010] The images or sketches of well core data descriptions did not lend themselves to digital interpretation and modeling applications. These images were frequently scanned as graphics image files for subsequent use, but an analyst was still provided with a visual image indicating well core data descriptions based on original analysis and formatting of such information. The images could only be displayed in the form of a visual reference picture or image.
[0011] It has been recognized that formation rock in hydrocarbon reservoirs exhibits two levels of porosity, which have been identified as macroporosity and microporosity. Oil flow characteristics in the rock are markedly different based on the relative presence of macroporosity and microporosity. The identification and relative presence of each of these levels of porosity in the reservoir rock has been based on analysis and laboratory measurement of core samples from wells in the reservoir, especially Type 1 microporosity which is charged with oil. This is described by Clerke et al. GeoArabia 2008, Vol. 13, No. 4, Application of Thomeer Hyperbolas to Decode the Pore Systems, Facies, and Reservoir Properties of the Upper Jurassic Arab D Limestone, Ghawar Field, Saudi Arabia:
A "Rosetta Stone" Approach; and Clerke, SPE Journal 2009, Permeability, Relative Permeability, Microscopic Displacement Efficiency, and Pore Geometry of M_1 Bimodal Pore Systems in Arab D Limestone. The identification of the porosity level presence was linked to an understanding of the rock space architecture.
100121 Micron level core analysis, while accurate, generally lacked the integration to the rest of the digital interpretation and modeling applications and systems.
Hydrocarbons recovered from macroporosity have proven to be usually much larger than the microporosity-recovered volumes in the early years of field recovery. This has led to an imprecise and inaccurate forecast of the field ultimate recovery, i. e., the performance of the microporosity hydrocarbons has not been properly included.
100131 Existing geological modeling processes and applications have not satisfactorily taken into account the formation rock characteristic data and the porosity level presence indicated by core samples. Specifically, proper pore system and recovery process data has, so far as is known, rarely been acquired in sufficient statistical quantities and in appropriate coordination with the geological facies.
SUMMARY OF THE INVENTION
[0014] Briefly, the present provides a computer implemented method of forming with a computer system a model of petrophysical structure of a subsurface reservoir. The model is formed based on well core description data obtained from analysis of well core samples about subsurface features of rock formations adjacent well bores in the subsurface reservoir. An initial geostatistical model of spatial distribution of the petrophysical structure of the subsurface reservoir is formed based on the well core description data. A geological model of lithofacies for wells in the reservoir in which core have not been obtained is also formed. A digital core description of lithofacies of the subsurface reservoir at the well bores is obtained based on the geological model. Identifications of the presence of macroporosity and microporosity in the subsurface lithofacies are received.
A model is then formed of the petrophysical structure of the reservoir based on the well core description data, the digital core description and the identified presence of macroporosity and microporosity in the subsurface lithofacies.
[0015] The present invention also provides a data processing system for computerized simulation of a model of lithofacies of a subsurface reservoir. The data processing system includes a processor which forms an initial geostatistical model of spatial distribution of the petrophysical structure of the subsurface reservoir based on the well core description data. The processor also forms a geological model of lithofacies for wells in the reservoir in which core have not been obtained, and obtains a digital core description of lithofacies of the subsurface reservoir at the well bores based on the geological model. The processor receives identifications of the presence of macroporosity and microporosity in the subsurface lithofacies, and forms a model of the petrophysical structure of the reservoir based on the well core description data, the digital core description and the identified presence of macroporosity and microporosity in the subsurface lithofacies.
[0016] The present invention further provides a data storage device having stored in a non-transitory computer readable medium computer operable instructions for causing a data processing system to form a model of petrophysical structure of a subsurface reservoir based on well core description data obtained from analysis of well core samples about subsurface features of rock formations adjacent well bores in the subsurface reservoir. The instructions stored in the computer readable medium causing the data processing system to form an initial geostatistical model of spatial distribution of the petrophysical structure of the subsurface reservoir based on the well core description data.
The instructions in the computer readable medium also cause the data processing system to form a geological model of lithofacies for wells in the reservoir in which core have not been obtained, and to obtain a digital core description of lithofacies of the subsurface reservoir at the well bores based on the geological model. The instructions in the computer readable medium also cause the data processing system to receive identifications of the presence of macroporosity and microporosity in the subsurface lithofacies, and to form a model of the petrophysical structure of the reservoir based on the well core description data, the digital core description and the identified presence of macroporosity and microporosity in the subsurface lithofacies.
BRIEF DESCRIPTION OF THE DRAWINGS
[0017] Figure 1 is a functional block diagram of a set of data processing steps performed in a data processing system for reservoir geological modeling of subsurface earth formations according to the present invention.
[0018] Figure 2 is a further functional block diagram of a portion of the data processing steps of the functional block diagram of Figure 1 for reservoir geological modeling of subsurface earth formations according to the present invention.
[0019] Figure 2A is a display of data resulting from one of the processing steps of Figure 2.
[0020] Figure 2B is a display of data resulting from another of the processing steps of Figure 2.
[0021] Figure 2C is a geostatistical variogram of results obtained during the processing steps of Figure 2.
100051 Seismic data with increasing accuracy has permitted the cells in reservoir models of both types to be on the order of 25 meters areal (x and y axis) intervals.
For what are known as giant reservoirs, the number of cells are the least hundreds of millions, and reservoirs of what is known as giga-cell size (a billion cells or more) are encountered.
100061 Days or weeks of computer processing have are usually been spent on 3D
model attributes computations. However, the guiding geological facies were not well digitally integrated into the model. Billions of cell attributes calculations could be wasted due to inappropriate algorithms from the initial model guide.
100071 When characterizing and developing a reservoir field, a 3D
geological model of the reservoir covering the entire 3D reservoir has been required to be built to provide an accurate model for reservoir planning. Accurate indications of lithofacies are an essential input in a 3D geological modeling system. Lithofacies is based on data from wells and establishes as a basis to predict reservoir properties in cells with no well data. One of the sets of data available has been what is known are formation rock descriptions or characterizations which are made by analysts based on observations made from well core samples.
10008] Traditionally, well core data has been assembled and analyzed based on measurements and observations taken from well core samples obtained from formation rock adjacent well boreholes at a number of depths of interest in a well. So far as is known, the core description has been hand drawn and was thereafter available for use in the form of a paper copy or at best a scanned graphical image of the hand drawing. The important geological information for a reservoir model came in a wide variety of forma and formats.
This information included numerical information, text strings, graphical images, relationships and drawings. These have posed a challenge to incorporation into the reservoir model in a form that could be consistently manipulated with other data.
[0009] In addition, manual preparation of these types of well description data presented a tedious challenge in drawing standardized patterns of mineral composition, texture, sedimentary structures based on well core samples repeatedly throughout the length of the well bore. Correlating core description data among different wells which had been described by different people with different scale and size hand drawings was also often problematic.
An accurate representation of the exact shape and size of formation rock features of interest in the drawing for tens of occurrences and over several wells in a reservoir was difficult to obtain in a reliable and representative manner.
[0010] The images or sketches of well core data descriptions did not lend themselves to digital interpretation and modeling applications. These images were frequently scanned as graphics image files for subsequent use, but an analyst was still provided with a visual image indicating well core data descriptions based on original analysis and formatting of such information. The images could only be displayed in the form of a visual reference picture or image.
[0011] It has been recognized that formation rock in hydrocarbon reservoirs exhibits two levels of porosity, which have been identified as macroporosity and microporosity. Oil flow characteristics in the rock are markedly different based on the relative presence of macroporosity and microporosity. The identification and relative presence of each of these levels of porosity in the reservoir rock has been based on analysis and laboratory measurement of core samples from wells in the reservoir, especially Type 1 microporosity which is charged with oil. This is described by Clerke et al. GeoArabia 2008, Vol. 13, No. 4, Application of Thomeer Hyperbolas to Decode the Pore Systems, Facies, and Reservoir Properties of the Upper Jurassic Arab D Limestone, Ghawar Field, Saudi Arabia:
A "Rosetta Stone" Approach; and Clerke, SPE Journal 2009, Permeability, Relative Permeability, Microscopic Displacement Efficiency, and Pore Geometry of M_1 Bimodal Pore Systems in Arab D Limestone. The identification of the porosity level presence was linked to an understanding of the rock space architecture.
100121 Micron level core analysis, while accurate, generally lacked the integration to the rest of the digital interpretation and modeling applications and systems.
Hydrocarbons recovered from macroporosity have proven to be usually much larger than the microporosity-recovered volumes in the early years of field recovery. This has led to an imprecise and inaccurate forecast of the field ultimate recovery, i. e., the performance of the microporosity hydrocarbons has not been properly included.
100131 Existing geological modeling processes and applications have not satisfactorily taken into account the formation rock characteristic data and the porosity level presence indicated by core samples. Specifically, proper pore system and recovery process data has, so far as is known, rarely been acquired in sufficient statistical quantities and in appropriate coordination with the geological facies.
SUMMARY OF THE INVENTION
[0014] Briefly, the present provides a computer implemented method of forming with a computer system a model of petrophysical structure of a subsurface reservoir. The model is formed based on well core description data obtained from analysis of well core samples about subsurface features of rock formations adjacent well bores in the subsurface reservoir. An initial geostatistical model of spatial distribution of the petrophysical structure of the subsurface reservoir is formed based on the well core description data. A geological model of lithofacies for wells in the reservoir in which core have not been obtained is also formed. A digital core description of lithofacies of the subsurface reservoir at the well bores is obtained based on the geological model. Identifications of the presence of macroporosity and microporosity in the subsurface lithofacies are received.
A model is then formed of the petrophysical structure of the reservoir based on the well core description data, the digital core description and the identified presence of macroporosity and microporosity in the subsurface lithofacies.
[0015] The present invention also provides a data processing system for computerized simulation of a model of lithofacies of a subsurface reservoir. The data processing system includes a processor which forms an initial geostatistical model of spatial distribution of the petrophysical structure of the subsurface reservoir based on the well core description data. The processor also forms a geological model of lithofacies for wells in the reservoir in which core have not been obtained, and obtains a digital core description of lithofacies of the subsurface reservoir at the well bores based on the geological model. The processor receives identifications of the presence of macroporosity and microporosity in the subsurface lithofacies, and forms a model of the petrophysical structure of the reservoir based on the well core description data, the digital core description and the identified presence of macroporosity and microporosity in the subsurface lithofacies.
[0016] The present invention further provides a data storage device having stored in a non-transitory computer readable medium computer operable instructions for causing a data processing system to form a model of petrophysical structure of a subsurface reservoir based on well core description data obtained from analysis of well core samples about subsurface features of rock formations adjacent well bores in the subsurface reservoir. The instructions stored in the computer readable medium causing the data processing system to form an initial geostatistical model of spatial distribution of the petrophysical structure of the subsurface reservoir based on the well core description data.
The instructions in the computer readable medium also cause the data processing system to form a geological model of lithofacies for wells in the reservoir in which core have not been obtained, and to obtain a digital core description of lithofacies of the subsurface reservoir at the well bores based on the geological model. The instructions in the computer readable medium also cause the data processing system to receive identifications of the presence of macroporosity and microporosity in the subsurface lithofacies, and to form a model of the petrophysical structure of the reservoir based on the well core description data, the digital core description and the identified presence of macroporosity and microporosity in the subsurface lithofacies.
BRIEF DESCRIPTION OF THE DRAWINGS
[0017] Figure 1 is a functional block diagram of a set of data processing steps performed in a data processing system for reservoir geological modeling of subsurface earth formations according to the present invention.
[0018] Figure 2 is a further functional block diagram of a portion of the data processing steps of the functional block diagram of Figure 1 for reservoir geological modeling of subsurface earth formations according to the present invention.
[0019] Figure 2A is a display of data resulting from one of the processing steps of Figure 2.
[0020] Figure 2B is a display of data resulting from another of the processing steps of Figure 2.
[0021] Figure 2C is a geostatistical variogram of results obtained during the processing steps of Figure 2.
-6-[0022] Figure 2D example display of a 3-D geological model with facies distribution obtained during the processing steps of Figure 2.
[0023] Figure 3 is a display of data resulting from the digital core description processing step of Figure 1.
[0024] Figures 3A, 3B, 3C and 3D are enlarged displays of portions of the data display of Figure 3.
[0025] Figure 4 is a schematic diagram of inputs used in upscaling water saturation in processing according to the present invention.
[0026] Figure 4A is a display of plots in well log format of example modeling data from processing according to the present invention.
[0027] Figure 4B is a display of 3-dimensional data plots of example modeling data from processing according to the present invention.
[0028] Figure 5 is a display like that of Figure 4A of a plot in well log format of example modeling data from processing according to the present invention.
[0029] Figures 5A is an enlarged display of portions of the data display of Figure 5.
[0030] Figures 6A, 6B, 6C and 6D are further displays of example modeling data from processing according to the present invention.
[0031] Figure 7 is an enlarged thin section image of a portion of a core sample showing actual macropore and micropore structure.
[0032] Figure 8 is a capillary pressure curve display of measurements from laboratory of well core plugs to investigate the relative presence of macroporosity and microporosity.
[0023] Figure 3 is a display of data resulting from the digital core description processing step of Figure 1.
[0024] Figures 3A, 3B, 3C and 3D are enlarged displays of portions of the data display of Figure 3.
[0025] Figure 4 is a schematic diagram of inputs used in upscaling water saturation in processing according to the present invention.
[0026] Figure 4A is a display of plots in well log format of example modeling data from processing according to the present invention.
[0027] Figure 4B is a display of 3-dimensional data plots of example modeling data from processing according to the present invention.
[0028] Figure 5 is a display like that of Figure 4A of a plot in well log format of example modeling data from processing according to the present invention.
[0029] Figures 5A is an enlarged display of portions of the data display of Figure 5.
[0030] Figures 6A, 6B, 6C and 6D are further displays of example modeling data from processing according to the present invention.
[0031] Figure 7 is an enlarged thin section image of a portion of a core sample showing actual macropore and micropore structure.
[0032] Figure 8 is a capillary pressure curve display of measurements from laboratory of well core plugs to investigate the relative presence of macroporosity and microporosity.
-7-[0033] Figure 9 is a schematic block diagram of a data processing system for three-dimensional multi-modal core and geophysical modeling according to the present invention.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
[0034] With the present invention, petrophysical properties and attributes of subsurface reservoirs are determined and modeled. The petrophysical properties provide a cumulative representation of the nature and structure of the formations in the subsurface reservoir. The model so formed is based on the pertinent geological facies, derived from well core description and detailed studies of rock, as well as fluid and pore properties (Full Pore System) obtained from laboratory analyses of core material and well log data.
These data differentiate various important pore throat and pore body regions and relationships, i.e., Macroporosity, and Type 1, 2 and 3 micropores.
[0035] The present invention thus provides an understanding of hydrocarbon volumes in the various pore type groups, which allows the establishment of proper recovery techniques through focused laboratory studies. This yields a field development strategy that can significantly increase hydrocarbon recovery from a reservoir.
[0036] As will be set forth, the present invention begins with 3-dimensional geo-cellular models of giant reservoirs with facies and full pore system attributes.
Determinations are then made using pore system type guidance from digitally described cores and rock core plugs to the sub- micron level, and identification and categorization of the macroporosity and multiple microporosity types. Multi-modal 3-dimensional modeling programs can then be used to form output models of the petrophysical properties and attributes of subsurface reservoirs. The processing results of the present invention contribute to optimal macro-micro reservoir recovery strategy and can lead potentially to billions of barrels of production in improved reservoir recovery.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
[0034] With the present invention, petrophysical properties and attributes of subsurface reservoirs are determined and modeled. The petrophysical properties provide a cumulative representation of the nature and structure of the formations in the subsurface reservoir. The model so formed is based on the pertinent geological facies, derived from well core description and detailed studies of rock, as well as fluid and pore properties (Full Pore System) obtained from laboratory analyses of core material and well log data.
These data differentiate various important pore throat and pore body regions and relationships, i.e., Macroporosity, and Type 1, 2 and 3 micropores.
[0035] The present invention thus provides an understanding of hydrocarbon volumes in the various pore type groups, which allows the establishment of proper recovery techniques through focused laboratory studies. This yields a field development strategy that can significantly increase hydrocarbon recovery from a reservoir.
[0036] As will be set forth, the present invention begins with 3-dimensional geo-cellular models of giant reservoirs with facies and full pore system attributes.
Determinations are then made using pore system type guidance from digitally described cores and rock core plugs to the sub- micron level, and identification and categorization of the macroporosity and multiple microporosity types. Multi-modal 3-dimensional modeling programs can then be used to form output models of the petrophysical properties and attributes of subsurface reservoirs. The processing results of the present invention contribute to optimal macro-micro reservoir recovery strategy and can lead potentially to billions of barrels of production in improved reservoir recovery.
-8-[0037] The attributes may be for example water saturation, porosity and permeability, and are determined directly at each cell of a three dimensional grid of a subsurface hydrocarbon reservoir. Input parameters based on original data from petrophysical measurements and other attributes of the reservoir and its subsurface features are generated for the individual cells. The processing of inputs for petrophysical properties and attributes is implemented by a suitable data processing system. The data processing system can be a mainframe computer of any conventional type of suitable processing capacity, or a cluster computer of a suitable number of processor nodes. An example of such a data processing system is a Linux Cluster arrangement which is commercially available. Other digital processors, however, may also be used, such as a laptop computer, or any other suitable processing apparatus.
It should thus be understood that a number of commercially available data processing systems may be used for this purpose.
[0038] A flowchart F (Figure 1) indicates the basic computer processing sequence of the present invention and the cluster computation taking place for a reservoir property or attribute according to the present invention. The processing sequence of the flow chart F is performed separately for each reservoir petrophysical property or attribute of the formations of interest in the reservoir being modeled.
[0039] 3D Geological Model (Step 100): Step 100 is performed by assembling an existing geological model for processing according to the present invention. The general parameters read in during step 100 initially include the following: rock type, porosity for rock type, well pore entrance size, individual cell dimensions and locations in the x, y and z directions;
petrophysical measurements and known data values of parameters and attributes from core sample data; and data available from well logs in areas of the reservoir where actual data have been obtained. The reservoir data is organized into a 3-dimensional (x, y, z) grid of cells according to the dimensions and volume of the reservoir. Depending on the size of the
It should thus be understood that a number of commercially available data processing systems may be used for this purpose.
[0038] A flowchart F (Figure 1) indicates the basic computer processing sequence of the present invention and the cluster computation taking place for a reservoir property or attribute according to the present invention. The processing sequence of the flow chart F is performed separately for each reservoir petrophysical property or attribute of the formations of interest in the reservoir being modeled.
[0039] 3D Geological Model (Step 100): Step 100 is performed by assembling an existing geological model for processing according to the present invention. The general parameters read in during step 100 initially include the following: rock type, porosity for rock type, well pore entrance size, individual cell dimensions and locations in the x, y and z directions;
petrophysical measurements and known data values of parameters and attributes from core sample data; and data available from well logs in areas of the reservoir where actual data have been obtained. The reservoir data is organized into a 3-dimensional (x, y, z) grid of cells according to the dimensions and volume of the reservoir. Depending on the size of the
-9-reservoir, the number of cells in the model can be a billion or more. Once data modeling begins, and additional data is obtained from newer wells in the reservoir and other data sources, the existing 3-dimensional petrophysical model data is supplemented by the additional data which is then included and used as part of the existing geological model.
[0040] Geologically-Controlled Geostatistics (Step 102): During step 102, geostatistical modeling is performed. The focus of the processing is to provide a geostatistical model of the postulated spatial distribution of presence of lithofacies over the entire reservoir based on the existing lithofacies data from cored wells in the geological model from step 100. The processing during step 102 uses the existing lithofacies data and forms a geostatistical model.
Selected values for the geological features of interest chosen by a geologist or analyst from the existing lithofacies data are used in the geostatistical processing. For different geological facies in the reservoir, the geologist may apply different geostatistical modeling criteria or a different processing schema In this way the geologist may control the modeling based on geological considerations rather than having the processing governed solely on statistical or analytical methodology. There are a number of presently commercially available methods for geostatistical modeling for this purpose. An example is Roxar Object-based Modeling available from Roxar ASA of Stavanger, Norway, through Emerson Process Management of St. Louis, Missouri.
[0041] Lithofacies Modeling (Step 104): During step 104, information regarding the lithofacies based on actual measurements and observations, which are termed ground truth, taken from core sample is developed. Details of the processing during step 102 are illustrated in Figure 2. In step 104A shown in Figure 2, processing begins with assembling petrophysical log data such as that shown in the display of Figure 2A from the core description data for an existing well in the reservoir. During step 104B, the cored well
[0040] Geologically-Controlled Geostatistics (Step 102): During step 102, geostatistical modeling is performed. The focus of the processing is to provide a geostatistical model of the postulated spatial distribution of presence of lithofacies over the entire reservoir based on the existing lithofacies data from cored wells in the geological model from step 100. The processing during step 102 uses the existing lithofacies data and forms a geostatistical model.
Selected values for the geological features of interest chosen by a geologist or analyst from the existing lithofacies data are used in the geostatistical processing. For different geological facies in the reservoir, the geologist may apply different geostatistical modeling criteria or a different processing schema In this way the geologist may control the modeling based on geological considerations rather than having the processing governed solely on statistical or analytical methodology. There are a number of presently commercially available methods for geostatistical modeling for this purpose. An example is Roxar Object-based Modeling available from Roxar ASA of Stavanger, Norway, through Emerson Process Management of St. Louis, Missouri.
[0041] Lithofacies Modeling (Step 104): During step 104, information regarding the lithofacies based on actual measurements and observations, which are termed ground truth, taken from core sample is developed. Details of the processing during step 102 are illustrated in Figure 2. In step 104A shown in Figure 2, processing begins with assembling petrophysical log data such as that shown in the display of Figure 2A from the core description data for an existing well in the reservoir. During step 104B, the cored well
-10-petrophysical log description data from step 104A is subjected to neural network or comparable processing of a suitable type to postulate lithofacies data for wells from which no core samples have been taken. Figure 2B is an example geostatistical variogram developed in the model building processing of step 104B. In step 104C facies logs are formed of the neural network processing results of step 104B. Figure 2C is a display resulting from step 104C of facies logs from uncored wells obtained during processing by step 104B. During step 104D, the facies logs developed during steps 104B are upscaled from well log scale to 3-dimensional model scale for the reservoir. During step 104E a facies model of the results from upscaling during step 104D is formed. Figure 2D is a visualization of a 3-dimensional geological model with facies distribution formed during step 104E.
[0042] Step 106 Digital Core Description: During step 106 (Figure 1), the processing is performed according to the procedure described in U. S. Patent Application Publication No.
2013/0080133 entitled "Core Plugs to Giga-Cells Lithological Modeling", published March 28, 2013, naming one of applicants as co-inventor. The subject matter of this publication may be referred to for further details. During processing in step 106 the reservoir characterization and interpretation developed during steps 102 and 104 are integrated into the digital data description.
[0043] Figure 3 is an example display formed as result of digital core description processing during step 106. The left seven columns in Figure 3 (Stylolites &
Fractures, Pore Type, Composition & Porosity, Sedimentary Structures, Texture, Cycle/Sequence Hierarchy, Fossils) are digital core descriptions of the formation lithology. Figure 3A
is an enlarged portion of the display of Figure 3 indicating numerical values for designated areas. Figure 3B
is an enlarged portion of a portion of Figure 3 indicating schematically sedimentary structures of the subsurface lithology at the indicated area. Figure 3C is an enlarged portion of a portion
[0042] Step 106 Digital Core Description: During step 106 (Figure 1), the processing is performed according to the procedure described in U. S. Patent Application Publication No.
2013/0080133 entitled "Core Plugs to Giga-Cells Lithological Modeling", published March 28, 2013, naming one of applicants as co-inventor. The subject matter of this publication may be referred to for further details. During processing in step 106 the reservoir characterization and interpretation developed during steps 102 and 104 are integrated into the digital data description.
[0043] Figure 3 is an example display formed as result of digital core description processing during step 106. The left seven columns in Figure 3 (Stylolites &
Fractures, Pore Type, Composition & Porosity, Sedimentary Structures, Texture, Cycle/Sequence Hierarchy, Fossils) are digital core descriptions of the formation lithology. Figure 3A
is an enlarged portion of the display of Figure 3 indicating numerical values for designated areas. Figure 3B
is an enlarged portion of a portion of Figure 3 indicating schematically sedimentary structures of the subsurface lithology at the indicated area. Figure 3C is an enlarged portion of a portion
-11-of Figure 3 which indicates graphically the carbonate texture and grain size of the subsurface lithology in the indicated region. Figure 3D is an enlarged portion of a portion of Figure 3 indicating an example thin section image obtained from a core sample taken at an indicated from a well in the reservoir being modeled.
10044] Figure 5 is a portion of a log plot of petrophysical data from a facies log formed as a result of the processing during step 104B (Figure 2) and Figure 5A is an enlarged plot of a portion indicated at 5A in the plot of Figure 5. Table 1 below is a replica of the legend shown at the top of Figures 5 and 5A identifying the facies indicated in these log plots. It should be understood that other indicia of facies may be used or included, as well.
TABLE I
Environ. Facies Fluv Ch (Fluvial Channel) Fluv Bar (Fluvial Bar) Shaly SS Bar (Shaly Sand Stone Bar) Flood Plain Coal Swamp Lake Mud Clean Coastal Tidal Sands Coastal/tidal Sand Stone < 15% Mud Heterolithic Marg. Marine/Deltaic Sand Heterolithic Marine Mud Pro-deltaic Mud Bay/Offshore Mud-Shale Ironstone Glauconitic Sands Carbonate
10044] Figure 5 is a portion of a log plot of petrophysical data from a facies log formed as a result of the processing during step 104B (Figure 2) and Figure 5A is an enlarged plot of a portion indicated at 5A in the plot of Figure 5. Table 1 below is a replica of the legend shown at the top of Figures 5 and 5A identifying the facies indicated in these log plots. It should be understood that other indicia of facies may be used or included, as well.
TABLE I
Environ. Facies Fluv Ch (Fluvial Channel) Fluv Bar (Fluvial Bar) Shaly SS Bar (Shaly Sand Stone Bar) Flood Plain Coal Swamp Lake Mud Clean Coastal Tidal Sands Coastal/tidal Sand Stone < 15% Mud Heterolithic Marg. Marine/Deltaic Sand Heterolithic Marine Mud Pro-deltaic Mud Bay/Offshore Mud-Shale Ironstone Glauconitic Sands Carbonate
-12-[0045] Figures 6A, 6B, 6C and 6D are well log plots of permeability as an attribute from horizontal wells formed as a result of the processing during step 104B. These plots are wireline and fluid flow meter logs measured at wells as verification with the present invention. Plots of such attributes according to the present invention are useful and helpful information to analysts and geologists in planning for improved recovery of hydrocarbons from the reservoir.
[0046] Macroporosity and Microporosity Identification (Step 108): During step 108, the pore system architecture of the 3-dimensional geological model is analyzed and identified based on data obtained from testing of core samples. The identification and analysis may be performed, for example, according to the techniques described in a publication of the Society of Petroleum Engineers: "Permeability, Relative Permeability, Microscopic Displacement Efficiency and Pore Geometry of M_1 Bimodal Pore Systems in Arab D Limestone"
Society of Petroleum Engineers Journal, SPE J. 105259 (2009) (Clerke); "Application of Thomeer Hyperbolas to Decode the Pore Systems, Facies, and Reservoir Properties of the Upper Jurassic Arab D Limestone, Ghawar Field, Saudi Arabia: A 'Rosetta Stone' Approach", GeoArabia 2008, Vol. 13, No. 4, (Clerke, et al). Investigation and analysis according to these techniques providing information about pore geometrical parameters for the geological model is conducted down to the sub-micron level. The pore system architecture analysis performed yields data regarding macroporosity (Type M) and microporosity (1, 2, and 3) properties of facies in the geological model.
[0047] In Figure 7 macroporosity in the form of a Type M macropore in a thin section image of a lithological sample is indicated by an arrow 70, while microporosity in the form of Type 1 micropores in the sample is indicated by arrows 72 and 74. An enlarged arrow 76
[0046] Macroporosity and Microporosity Identification (Step 108): During step 108, the pore system architecture of the 3-dimensional geological model is analyzed and identified based on data obtained from testing of core samples. The identification and analysis may be performed, for example, according to the techniques described in a publication of the Society of Petroleum Engineers: "Permeability, Relative Permeability, Microscopic Displacement Efficiency and Pore Geometry of M_1 Bimodal Pore Systems in Arab D Limestone"
Society of Petroleum Engineers Journal, SPE J. 105259 (2009) (Clerke); "Application of Thomeer Hyperbolas to Decode the Pore Systems, Facies, and Reservoir Properties of the Upper Jurassic Arab D Limestone, Ghawar Field, Saudi Arabia: A 'Rosetta Stone' Approach", GeoArabia 2008, Vol. 13, No. 4, (Clerke, et al). Investigation and analysis according to these techniques providing information about pore geometrical parameters for the geological model is conducted down to the sub-micron level. The pore system architecture analysis performed yields data regarding macroporosity (Type M) and microporosity (1, 2, and 3) properties of facies in the geological model.
[0047] In Figure 7 macroporosity in the form of a Type M macropore in a thin section image of a lithological sample is indicated by an arrow 70, while microporosity in the form of Type 1 micropores in the sample is indicated by arrows 72 and 74. An enlarged arrow 76
-13-indicates schematically entry of water comes in to the pore space of the sample via the macropore and paths of the water into the micropores. Arrows 78 indicate hydrocarbons movement from micropore to macropore. Figure 8 shows the standard capillary pressure curve display of the type described in Clerke et al. GeoArabia 2008, Vol. 13, No. 4, identified above. Multiple lines indicated generally at 80A on the left portion of the display are laboratory measurements from multiple well core plugs in a reservoir. The right portion of the display in Figure 8 shows the occurrence frequency statistical display identifying at 80 macroporosity M and at 82, 84, and 86 the respective microporosity 1, 2, and 3 properties.
[0048] Multi-Modal 3-D Petrophysical Modeling (Step 110): Petrophysical modeling of the reservoir based on the input data, the digital core description data and the identified presence of macroporosity and microporosity in the cells of the reservoir. The modeling methodology utilized is preferably a multi-modal petrophysical algorithm. The choice of modeling utilized is based in part on the attributes of interest. The modeling may be performed, for example, according to the 3-dimensional modeling methodology of the type described in co-pending U. S. Patent Application No. 13/913,086, "Cluster Petrophysical Uncertainty Modeling", filed July 28, 2011, naming one of applicants as co-inventor.
[0049] The petrophysical models formed during step 110 may be of a number of forms.
The models may be models of petrophysical structure, such as facies logs of the type shown in Figures 5 and 5A; and facies models of the type shown in Figure 2D. The models of petrophysical properties may also be of reservoir attributes, either in log plot form (Figures 6A through 6D, inclusive) or in 3-dimensional form as shown in Figure 4B
described below.
[0050] Figure 4 illustrates schematically an example procedure for determining values for upscale water saturation S,UP as a reservoir attribute according to the present invention.
Input parameters including: crude, brine, and rock properties in-situ as indicated at 110A;
[0048] Multi-Modal 3-D Petrophysical Modeling (Step 110): Petrophysical modeling of the reservoir based on the input data, the digital core description data and the identified presence of macroporosity and microporosity in the cells of the reservoir. The modeling methodology utilized is preferably a multi-modal petrophysical algorithm. The choice of modeling utilized is based in part on the attributes of interest. The modeling may be performed, for example, according to the 3-dimensional modeling methodology of the type described in co-pending U. S. Patent Application No. 13/913,086, "Cluster Petrophysical Uncertainty Modeling", filed July 28, 2011, naming one of applicants as co-inventor.
[0049] The petrophysical models formed during step 110 may be of a number of forms.
The models may be models of petrophysical structure, such as facies logs of the type shown in Figures 5 and 5A; and facies models of the type shown in Figure 2D. The models of petrophysical properties may also be of reservoir attributes, either in log plot form (Figures 6A through 6D, inclusive) or in 3-dimensional form as shown in Figure 4B
described below.
[0050] Figure 4 illustrates schematically an example procedure for determining values for upscale water saturation S,UP as a reservoir attribute according to the present invention.
Input parameters including: crude, brine, and rock properties in-situ as indicated at 110A;
-14-porosity architecture 110B; historical free water level and residual oil data 110C; wettability alterations 110D; as well as parameters indicating the nature of upscaling 110E; and the nature and presence of multiple formation pore systems 110F are provided to the data processor, along with the petrophysical processing algorithm for determining upscale water saturation. A display of modeling results of a reservoir attribute, such as upscale water saturation SwUP, is then formed.
[0051] Figure 4A is a display of modeling the results of processing according to Figure 4 at an individual well. Figure 4A is an example log plot display as a function of depth for well data obtained as a result of such processing. As indicated, porosity, modeled permeability, petrophysical rock types (PRT's) and cap pressure (Pc) attributes are displayed as parameters.
Figure 4B is an example display of modeling the results of processing according to Figure 4 at a 3-dimensional model level. As indicated, porosity 4), entry cap pressure Pd and pore-size variability curvature G are displayed as attributes. SwUP determination in Figure 4 indicates the various input data types which go into the calculations with the present invention. Figure 4A shows the input well logs and output calculated log visualization at a single well. Figure 4B shows the same input and output but visualized in the 3D geological model format.
[0052] Modify and Adjust Field Development Strategy (Step 112): This step involves geologists and analysts evaluating the processing results for the purpose of modifying and adjusting field development strategy based upon the best indicative model generated from steps 100 through 110. By doing such optimization, the ultimate recovery of potentially billions of additional barrels of hydrocarbon reserves is possible.
[0053] Processing according the present invention as illustrated in Figure 1 can be considered to be a continuous cycle or loop. The "loop" means it is to be repeated and is a continuous process. Starting with an initial geological model, the various processing steps
[0051] Figure 4A is a display of modeling the results of processing according to Figure 4 at an individual well. Figure 4A is an example log plot display as a function of depth for well data obtained as a result of such processing. As indicated, porosity, modeled permeability, petrophysical rock types (PRT's) and cap pressure (Pc) attributes are displayed as parameters.
Figure 4B is an example display of modeling the results of processing according to Figure 4 at a 3-dimensional model level. As indicated, porosity 4), entry cap pressure Pd and pore-size variability curvature G are displayed as attributes. SwUP determination in Figure 4 indicates the various input data types which go into the calculations with the present invention. Figure 4A shows the input well logs and output calculated log visualization at a single well. Figure 4B shows the same input and output but visualized in the 3D geological model format.
[0052] Modify and Adjust Field Development Strategy (Step 112): This step involves geologists and analysts evaluating the processing results for the purpose of modifying and adjusting field development strategy based upon the best indicative model generated from steps 100 through 110. By doing such optimization, the ultimate recovery of potentially billions of additional barrels of hydrocarbon reserves is possible.
[0053] Processing according the present invention as illustrated in Figure 1 can be considered to be a continuous cycle or loop. The "loop" means it is to be repeated and is a continuous process. Starting with an initial geological model, the various processing steps
-15-are performed to produce a new geological model. Field development in the form of a field development strategy is based on the model to drill additional wells for production. As new wells are being drilled, new well logs and cores are being acquired. With the new information in hand, the current best 3D geological model is being continuously updated.
The continuous field development furnishes a roadmap for further development of the reservoir in the models being formed. The new data and interpretation is incorporated into the existing model and processing loops back to Step 100 as schematically indicated in Figure 1 to update the model and repeat the process.
[0054] As illustrated in Fig. 9, a data processing system D according to the present invention includes a computer 150 having a processor 152 and memory 154 coupled to processor 152 to store operating instructions, control information and database records therein. The computer 150 may, if desired, be a portable digital processor, such as a personal computer in the form of a laptop computer, notebook computer or other suitable programmed or programmable digital data processing apparatus, such as a desktop computer.
It should also be understood that the computer 150 may be a multicore processor with nodes such as those from Intel Corporation or Advanced Micro Devices (A1VID), an HPC Linux cluster computer or a mainframe computer of any conventional type of suitable processing capacity such as those available from International Business Machines (IBM) of Armonk, N.Y. or other source.
[0055] The computer 150 has a user interface 156 and an output data or graphical user display 158 for displaying output data or records of lithological facies and reservoir attributes according to the present invention. The output display 158 includes components such as a printer and an output display screen capable of providing printed output information or
The continuous field development furnishes a roadmap for further development of the reservoir in the models being formed. The new data and interpretation is incorporated into the existing model and processing loops back to Step 100 as schematically indicated in Figure 1 to update the model and repeat the process.
[0054] As illustrated in Fig. 9, a data processing system D according to the present invention includes a computer 150 having a processor 152 and memory 154 coupled to processor 152 to store operating instructions, control information and database records therein. The computer 150 may, if desired, be a portable digital processor, such as a personal computer in the form of a laptop computer, notebook computer or other suitable programmed or programmable digital data processing apparatus, such as a desktop computer.
It should also be understood that the computer 150 may be a multicore processor with nodes such as those from Intel Corporation or Advanced Micro Devices (A1VID), an HPC Linux cluster computer or a mainframe computer of any conventional type of suitable processing capacity such as those available from International Business Machines (IBM) of Armonk, N.Y. or other source.
[0055] The computer 150 has a user interface 156 and an output data or graphical user display 158 for displaying output data or records of lithological facies and reservoir attributes according to the present invention. The output display 158 includes components such as a printer and an output display screen capable of providing printed output information or
-16-visible displays in the form of graphs, data sheets, graphical images, data plots and the like as output records or images, [0056] The user interface 156 of computer 150 also includes a suitable user input device or input/output control unit 160 to provide a user access to control or access information and database records and operate the computer 150. Data processing system D
further includes a database 162 stored in computer memory, which may be internal memory 154, or an external, networked, or non-networked memory as indicated at 166 in an associated database server 168.
[0057] The data processing system D includes program code 170 stored in memory 154 of the computer 150. The program code 170, according to the present invention is in the form of non-transitory computer operable instructions causing the data processor 152 to perform the computer implemented method of the present invention in the manner described above.
[0058] It should be noted that program code 170 may be in the form of microcode, programs, routines, or symbolic computer operable languages that provide a specific set of ordered operations that control the functioning of the data processing system D and direct its operation. The instructions of program code 170 may be stored in non-transitory form in memory 154 of the computer 150, or on computer diskette, magnetic tape, conventional hard disk drive, electronic read-only memory, optical storage device, or other appropriate data storage device having a non-transitory computer usable medium stored thereon.
Program code 170 may also be contained in non-transitory form on a data storage device such as server 168 as a computer readable medium.
[0059] The method of the present invention performed in the computer 150 can be implemented utilizing the computer program steps of Figures 1, 2 and 4 stored in memory
further includes a database 162 stored in computer memory, which may be internal memory 154, or an external, networked, or non-networked memory as indicated at 166 in an associated database server 168.
[0057] The data processing system D includes program code 170 stored in memory 154 of the computer 150. The program code 170, according to the present invention is in the form of non-transitory computer operable instructions causing the data processor 152 to perform the computer implemented method of the present invention in the manner described above.
[0058] It should be noted that program code 170 may be in the form of microcode, programs, routines, or symbolic computer operable languages that provide a specific set of ordered operations that control the functioning of the data processing system D and direct its operation. The instructions of program code 170 may be stored in non-transitory form in memory 154 of the computer 150, or on computer diskette, magnetic tape, conventional hard disk drive, electronic read-only memory, optical storage device, or other appropriate data storage device having a non-transitory computer usable medium stored thereon.
Program code 170 may also be contained in non-transitory form on a data storage device such as server 168 as a computer readable medium.
[0059] The method of the present invention performed in the computer 150 can be implemented utilizing the computer program steps of Figures 1, 2 and 4 stored in memory
-17-154 and executable by system processor 152 of computer 150. The input data to processing system D are the input parameters and reservoir data of the types described above.
[0060] It can thus be seen that the present invention provides a loop processing methodology for determining geological model attributes of billions of reservoir cells based on geological facies guidance from well core, and then digitally describing associated pore systems to the sub-micron level, identifying the pertinent pore types e.g., macroporosity and multiple microporosity types, and then incorporating multi-modal petrophysical algorithms into the 3D modeling application. These 3-dimensional static models can then be meshed with the facies-based recovery processes to establish proper recovery strategy, which can lead to potentially billions of barrels of production increases through optimal recovery practices.
[0061] The present invention thus integrates traditionally independent processes and creates new applications and workflows to link and digitally calibrate different pieces to generate a sound scientific and business solution. Geological model attributes like the permeability and the relative permeability are not assigned but calculated for each geocell using transform equations acting upon pore system parameters. The reservoir attributes are far more strongly related to their rock types over their spatial location distribution. The processing determinations also incorporate the various geological facies interpreted at each model cell location.
[0062] Therefore, the geological facies play an important role in determining model computation algorithms. The geological facies are guided from the well log data and core descriptions which contain the ground truth. The process requires that all data to be fully digital. Core description is generally done in a manual drawing fashion. While accurate, it lacks the digital format which this new process requires. Macroporosity and microporosity data from core plug analyses or specific well log analyses are captured digitally for input.
[0060] It can thus be seen that the present invention provides a loop processing methodology for determining geological model attributes of billions of reservoir cells based on geological facies guidance from well core, and then digitally describing associated pore systems to the sub-micron level, identifying the pertinent pore types e.g., macroporosity and multiple microporosity types, and then incorporating multi-modal petrophysical algorithms into the 3D modeling application. These 3-dimensional static models can then be meshed with the facies-based recovery processes to establish proper recovery strategy, which can lead to potentially billions of barrels of production increases through optimal recovery practices.
[0061] The present invention thus integrates traditionally independent processes and creates new applications and workflows to link and digitally calibrate different pieces to generate a sound scientific and business solution. Geological model attributes like the permeability and the relative permeability are not assigned but calculated for each geocell using transform equations acting upon pore system parameters. The reservoir attributes are far more strongly related to their rock types over their spatial location distribution. The processing determinations also incorporate the various geological facies interpreted at each model cell location.
[0062] Therefore, the geological facies play an important role in determining model computation algorithms. The geological facies are guided from the well log data and core descriptions which contain the ground truth. The process requires that all data to be fully digital. Core description is generally done in a manual drawing fashion. While accurate, it lacks the digital format which this new process requires. Macroporosity and microporosity data from core plug analyses or specific well log analyses are captured digitally for input.
-18-Hydrocarbons contained in macroporosity and microporosity has separate flow behaviors in terms of permeability, relative permeability and flow rates. It is generally found that a sizable hydrocarbon volume is contained in some microporosity (Type 1) alongside the hydrocarbons in the macroporosity porosity. Most recovery practice focus only on the macro-porosity due to its large volume and rapid flow especially in the early stages of reservoir production.
Lack of understanding and proper analyses have limited the ability to derive specific production characteristic algorithms from the microporosity.
[0063] The present invention thus incorporates fully multi-modal petrophysical algorithms into a 3-dimensional geological modeling application. The value of this process is that proper modeling of the static and dynamic aspects of the hydrocarbons in the microporosity leads to improved optimal reservoir recovery strategies that could result in billions of barrels in additional recovery.
[0064] The invention has been sufficiently described so that a person with average knowledge in the matter may reproduce and obtain the results mentioned in the invention herein Nonetheless, any skilled person in the field of technique, subject of the invention herein, may carry out modifications not described in the request herein, to apply these modifications to a determined methodology, or in the performance of the same, requires the claimed matter in the following claims; such techniques and procedures shall be covered within the scope of the invention.
[0065] It should be noted and understood that there can be improvements and modifications made of the present invention described in detail above without departing from the spirit or scope of the invention as set forth in the accompanying claims.
Lack of understanding and proper analyses have limited the ability to derive specific production characteristic algorithms from the microporosity.
[0063] The present invention thus incorporates fully multi-modal petrophysical algorithms into a 3-dimensional geological modeling application. The value of this process is that proper modeling of the static and dynamic aspects of the hydrocarbons in the microporosity leads to improved optimal reservoir recovery strategies that could result in billions of barrels in additional recovery.
[0064] The invention has been sufficiently described so that a person with average knowledge in the matter may reproduce and obtain the results mentioned in the invention herein Nonetheless, any skilled person in the field of technique, subject of the invention herein, may carry out modifications not described in the request herein, to apply these modifications to a determined methodology, or in the performance of the same, requires the claimed matter in the following claims; such techniques and procedures shall be covered within the scope of the invention.
[0065] It should be noted and understood that there can be improvements and modifications made of the present invention described in detail above without departing from the spirit or scope of the invention as set forth in the accompanying claims.
-19-
Claims (13)
1. A computer implemented method of forming with a computer system a static geological model of petrophysical properties representing geological formations of a subsurface reservoir, including a presence of microporosity and macroporosity, the static geological model being organized into a number of cells arranged in an organized three-dimensional grid of cells according to dimensions and volume of the reservoir, the static geological model being formed based on well logs obtained from wells in the subsurface reservoir, rock types of the subsurface rock formations, and well core description data comprising petrophysical measurements obtained from analysis of well core samples of rock formations adjacent cored wells of the wells in the subsurface reservoir in which core samples have been obtained, the method comprising computer processing steps of:
(a) forming a geostatistical model of the presence of lithofacies over the subsurface reservoir based on the well logs obtained from wells in the subsurface reservoir and the well core samples from the cored wells in the reservoir in which well core samples have been obtained;
(b) forming a lithofacies model of postulated facies logs for the wells in the reservoir in which core samples have not been obtained, based on the petrophysical measurements from the well core samples for the cored wells in the reservoir;
(c) forming a digital core description of a lithology of cells of the subsurface reservoir as functions of depth in well bores based on the formed initial geostatistical model of spatial distribution of the presence of lithofacies and the formed lithofacies model;
(d) receiving identifications of the presence of macroporosity and microporosity in the subsurface lithofacies in the three-dimensional grid of cells of the static geological model based on the testing of the well core samples obtained from the cored wells in the subsurface reservoir;
(e) forming the static geological model of the petrophysical properties of the reservoir representing the geological formations of the subsurface reservoir, including the presence of microporosity and macroporosity, by petrophysical uncertainty modeling based on: the well logs obtained from wells in the subsurface reservoir; the rock types of the subsurface rock formations; the well core description data obtained from analysis of the well core samples; and the identified presence of macroporosity and microporosity in the subsurface lithofacies; and (f) forming an output display of the static geological model of petrophysical properties representing in the cells, and including the presence of microporosity and macroporosity, the geological formations of the static geological model for development of the reservoir.
(a) forming a geostatistical model of the presence of lithofacies over the subsurface reservoir based on the well logs obtained from wells in the subsurface reservoir and the well core samples from the cored wells in the reservoir in which well core samples have been obtained;
(b) forming a lithofacies model of postulated facies logs for the wells in the reservoir in which core samples have not been obtained, based on the petrophysical measurements from the well core samples for the cored wells in the reservoir;
(c) forming a digital core description of a lithology of cells of the subsurface reservoir as functions of depth in well bores based on the formed initial geostatistical model of spatial distribution of the presence of lithofacies and the formed lithofacies model;
(d) receiving identifications of the presence of macroporosity and microporosity in the subsurface lithofacies in the three-dimensional grid of cells of the static geological model based on the testing of the well core samples obtained from the cored wells in the subsurface reservoir;
(e) forming the static geological model of the petrophysical properties of the reservoir representing the geological formations of the subsurface reservoir, including the presence of microporosity and macroporosity, by petrophysical uncertainty modeling based on: the well logs obtained from wells in the subsurface reservoir; the rock types of the subsurface rock formations; the well core description data obtained from analysis of the well core samples; and the identified presence of macroporosity and microporosity in the subsurface lithofacies; and (f) forming an output display of the static geological model of petrophysical properties representing in the cells, and including the presence of microporosity and macroporosity, the geological formations of the static geological model for development of the reservoir.
2. The computer implemented method of claim 1, wherein the step of forming a geostatistical model comprises the step of processing the well logs and the well core description data to determine lithofacies data for the cells of the reservoir model adjacent the well bores in the reservoir.
3. The computer implemented method of claim 2, wherein the step of forming a lithofacies model further includes the step of upscaling the postulated facies logs for wells in the reservoir in which core samples have not been obtained to a three-dimensional model of facies of the reservoir.
4. The computer implemented method of claim 3, wherein the step of forming a lithofacies model comprises the step of forming a facies model of results from the step of upscaling the postulated facies logs.
5. A data processing system for computerized simulation of a static geological model of petrophysical properties representing geological formations of a subsurface reservoir, including a presence of microporosity and macroporosity, the static geological model being organized into a number of cells arranged in an organized three-dimensional grid of cells according to dimensions and volume of the reservoir, the static geological model being formed based on well logs obtained from wells in the subsurface reservoir, rock types of the subsurface rock formations, and well core description data comprising petrophysical measurements obtained from analysis of well core samples of rock formations adjacent cored wells of the wells in the subsurface reservoir in which core samples have been obtained, the data processing system comprising:
(a) a processor performing the steps of:
(1) forming a geostatistical model of the presence of lithofacies over the subsurface reservoir based on the well logs obtained from wells in the subsurface reservoir and the well core samples from the cored wells in the reservoir in which well core samples have been obtained;
(2) forming a lithofacies model of postulated facies logs based on the petrophysical measurements from the well core samples for the cored wells in the reservoir;
(3) forming a digital core description of a lithology of cells of the subsurface reservoir as functions of depth in well bores based on the formed initial geostatistical model of spatial distribution of the presence of lithofacies and the formed lithofacies model;
(4) receiving identifications of the presence of macroporosity and microporosity in the subsurface lithofacies in the three-dimensional grid of cells of the static geological model based on the testing of the well core samples obtained from the cored wells in the subsurface reservoir;
(5) forming the static geological model of the petrophysical properties of the reservoir representing the geological formations of the subsurface reservoir, including the presence of microporosity and macroporosity, by petrophysical uncertainty modeling based on: the well logs obtained from wells in the subsurface reservoir; the rock types of the subsurface rock formations; the well core description data obtained from analysis of well core samples; and the identified presence of macroporosity and microporosity in the subsurface lithofacies; and (b) a display forming an output display of the static geological model of petrophysical properties representing in the cells, and including the presence of microporosity and macroporosity, the geological formations of the static geological model for development of the reservoir.
(a) a processor performing the steps of:
(1) forming a geostatistical model of the presence of lithofacies over the subsurface reservoir based on the well logs obtained from wells in the subsurface reservoir and the well core samples from the cored wells in the reservoir in which well core samples have been obtained;
(2) forming a lithofacies model of postulated facies logs based on the petrophysical measurements from the well core samples for the cored wells in the reservoir;
(3) forming a digital core description of a lithology of cells of the subsurface reservoir as functions of depth in well bores based on the formed initial geostatistical model of spatial distribution of the presence of lithofacies and the formed lithofacies model;
(4) receiving identifications of the presence of macroporosity and microporosity in the subsurface lithofacies in the three-dimensional grid of cells of the static geological model based on the testing of the well core samples obtained from the cored wells in the subsurface reservoir;
(5) forming the static geological model of the petrophysical properties of the reservoir representing the geological formations of the subsurface reservoir, including the presence of microporosity and macroporosity, by petrophysical uncertainty modeling based on: the well logs obtained from wells in the subsurface reservoir; the rock types of the subsurface rock formations; the well core description data obtained from analysis of well core samples; and the identified presence of macroporosity and microporosity in the subsurface lithofacies; and (b) a display forming an output display of the static geological model of petrophysical properties representing in the cells, and including the presence of microporosity and macroporosity, the geological formations of the static geological model for development of the reservoir.
6. The data processing system of claim 5, wherein the processor in forming a geostatistical model processes the well logs and the well core description data to determine lithofacies data for the cells of the reservoir model adjacent the well bores in the reservoir.
7. The data processing system of claim 6, wherein the processor in forming a lithofacies model further upscales the facies logs for wells in the reservoir in which core samples have not been obtained to a three-dimensional model of facies of the reservoir.
8. The data processing system of claim 5, wherein the processor in forming a lithofacies model forms a facies model of the upscaled facies logs.
9. A data storage device having stored in a non-transitory computer readable medium non-transitory computer operable instructions for causing a data processing system to form a static geological model of petrophysical properties representing geological formations of a subsurface reservoir, including a presence of microporosity and macroporosity, the static geological model being organized into a number of cells arranged in an organized three-dimensional grid of cells according to dimensions and volume of the reservoir, the static geological model being formed based on well logs obtained from wells in the subsurface reservoir, rock types of the subsurface rock formations, and well core description data comprising petrophysical measurements obtained from analysis of well core samples of rock formations adjacent cored wells in the subsurface reservoir, the instructions stored in the computer readable medium causing the data processing system to perform the steps of:
(a) forming a geostatistical model of the presence of lithofacies over the subsurface reservoir based on logs obtained from wells in the subsurface reservoir and the well core samples from the cored wells in the reservoir in which well core samples have been obtained;
(b) forming a lithofacies model of postulated facies logs for the wells in the reservoir in which core samples have not been obtained, based on the petrophysical measurements from the well core samples for the cored wells in the reservoir;
(c) forming a digital core description of a lithology of cells of the subsurface reservoir as functions of depth in well bores based on the formed initial geostatistical model of spatial distribution of the presence of lithofacies and the formed lithofacies model;
(d) receiving identifications of the presence of macroporosity and microporosity in the subsurface lithofacies in the three-dimensional grid of cells of the static geological model based on the testing of the well core samples obtained from the cored wells in the subsurface reservoir;
(e) forming the static geological model of the petrophysical properties of the reservoir representing the geological formations of the subsurface reservoir, including the presence of microporosity and macroporosity, by petrophysical uncertainty modeling based on: the well logs obtained from wells in the subsurface reservoir; the rock types of the subsurface rock formations; the well core description data obtained from analysis of the well core samples; and the identified presence of macroporosity and microporosity in the subsurface lithofacies; and (f) forming an output display of the static geological model of petrophysical properties representing in the cells, and including the presence of microporosity and macroporosity, the geological formations of the static geological model for development of the reservoir.
(a) forming a geostatistical model of the presence of lithofacies over the subsurface reservoir based on logs obtained from wells in the subsurface reservoir and the well core samples from the cored wells in the reservoir in which well core samples have been obtained;
(b) forming a lithofacies model of postulated facies logs for the wells in the reservoir in which core samples have not been obtained, based on the petrophysical measurements from the well core samples for the cored wells in the reservoir;
(c) forming a digital core description of a lithology of cells of the subsurface reservoir as functions of depth in well bores based on the formed initial geostatistical model of spatial distribution of the presence of lithofacies and the formed lithofacies model;
(d) receiving identifications of the presence of macroporosity and microporosity in the subsurface lithofacies in the three-dimensional grid of cells of the static geological model based on the testing of the well core samples obtained from the cored wells in the subsurface reservoir;
(e) forming the static geological model of the petrophysical properties of the reservoir representing the geological formations of the subsurface reservoir, including the presence of microporosity and macroporosity, by petrophysical uncertainty modeling based on: the well logs obtained from wells in the subsurface reservoir; the rock types of the subsurface rock formations; the well core description data obtained from analysis of the well core samples; and the identified presence of macroporosity and microporosity in the subsurface lithofacies; and (f) forming an output display of the static geological model of petrophysical properties representing in the cells, and including the presence of microporosity and macroporosity, the geological formations of the static geological model for development of the reservoir.
10. The data storage device of claim 9, wherein the instructions for forming an initial geostatistical model further cause the data processing system to process the well logs and the well core description data to determine lithofacies data for the cells of the reservoir model adjacent the well bores in the reservoir.
11. The data storage device of claim 10, wherein the instructions for forming a lithofacies model further cause the data processing system to upscale the facies logs for wells in the reservoir in which core samples have not been obtained to a three-dimensional model of facies of the reservoir.
12. The data storage device of claim 10, wherein the instructions for forming a lithofacies model further cause the data processing system to form a facies model of the upscaled facies logs.
13. The data storage device of claim 9, wherein the data processing system includes a data display and the instructions include instructions causing the data processing system to form an output display of the nature and structure of formations of the reservoir.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/IB2013/000978 WO2013164685A2 (en) | 2012-05-01 | 2013-05-13 | Three-dimensional multi-modal core and geological modeling for optimal field development |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2869825A1 CA2869825A1 (en) | 2013-11-07 |
| CA2869825C true CA2869825C (en) | 2018-07-24 |
Family
ID=51932135
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA2869825A Active CA2869825C (en) | 2013-05-13 | 2013-05-13 | Three-dimensional multi-modal core and geological modeling for optimal field development |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN104380144B (en) |
| CA (1) | CA2869825C (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9939548B2 (en) | 2014-02-24 | 2018-04-10 | Saudi Arabian Oil Company | Systems, methods, and computer medium to produce efficient, consistent, and high-confidence image-based electrofacies analysis in stratigraphic interpretations across multiple wells |
| EP3329091B1 (en) * | 2015-07-28 | 2022-12-28 | Services Pétroliers Schlumberger | Method and system for generating a virtual core |
| US20210165937A1 (en) * | 2017-12-14 | 2021-06-03 | Schlumberger Technology Corporation | System and Method for Simulating Reservoir Models |
| CN109215029B (en) * | 2018-08-29 | 2021-10-15 | 电子科技大学 | Segmentation and extraction method of three-dimensional geological abnormal body based on convolutional neural network |
| CN111260791B (en) * | 2018-11-30 | 2023-06-09 | 中国石油化工股份有限公司 | Method for updating geosteering model |
| CN111259508B (en) * | 2018-11-30 | 2023-04-11 | 中国石油化工股份有限公司 | Drilling process simulation method and system |
| CN112147680A (en) * | 2019-06-27 | 2020-12-29 | 中国石油化工股份有限公司 | Rock internal micro-crack extraction and characterization method and system |
| CN112287532B (en) * | 2020-10-20 | 2024-03-05 | 中海石油(中国)有限公司 | Edge control earthquake driving modeling method |
| CN113047830B (en) * | 2021-04-22 | 2022-05-03 | 中国石油大学(北京) | Reservoir flow unit identification method and device, computer and storage medium |
| CN117784737B (en) * | 2024-02-23 | 2024-05-10 | 西安众望能源科技有限公司 | Intelligent RTU control system for oilfield well site |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB9925957D0 (en) * | 1999-11-02 | 1999-12-29 | Geco As | Method and apparatus for generating a cross plot in attribute space from a plurality of attribute maps and generating a class from the cross plot. |
| WO2009029133A1 (en) * | 2007-08-24 | 2009-03-05 | Exxonmobil Upstream Research Company | Method for multi-scale geomechanical model analysis by computer simulation |
| US9134457B2 (en) * | 2009-04-08 | 2015-09-15 | Schlumberger Technology Corporation | Multiscale digital rock modeling for reservoir simulation |
| WO2013040281A2 (en) * | 2011-09-15 | 2013-03-21 | Saudi Arabian Oil Company | Core-plug to giga-cells lithological modeling |
-
2013
- 2013-05-13 CA CA2869825A patent/CA2869825C/en active Active
- 2013-05-13 CN CN201380023176.2A patent/CN104380144B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CA2869825A1 (en) | 2013-11-07 |
| CN104380144B (en) | 2018-01-16 |
| CN104380144A (en) | 2015-02-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9684084B2 (en) | Three-dimensional multi-modal core and geological modeling for optimal field development | |
| CA2869825C (en) | Three-dimensional multi-modal core and geological modeling for optimal field development | |
| Skalinski et al. | Carbonate petrophysical rock typing: integrating geological attributes and petrophysical properties while linking with dynamic behaviour | |
| US8217337B2 (en) | Evaluating a reservoir formation | |
| CA2847693C (en) | Core-plug to giga-cells lithological modeling | |
| US8359184B2 (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
| US9134457B2 (en) | Multiscale digital rock modeling for reservoir simulation | |
| US20200095858A1 (en) | Modeling reservoir permeability through estimating natural fracture distribution and properties | |
| US7054749B1 (en) | Method for determining reservoir fluid volumes, fluid contacts, compartmentalization, and permeability in geological subsurface models | |
| Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
| EP4042211B1 (en) | Modeling reservoir permeability through estimating natural fracture distribution and properties | |
| Harris | The role of geology in reservoir simulation studies | |
| Oeth et al. | Characterization of small scale heterogeneity to predict acid fracture performance | |
| Kamali et al. | 3D geostatistical modeling and uncertainty analysis in a carbonate reservoir, SW Iran | |
| US8718992B2 (en) | Method for history matching of a geological model comprising a sub-seismic fault network | |
| Benetatos et al. | Coping with uncertainties through an automated workflow for 3D reservoir modelling of carbonate reservoirs | |
| Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
| Orellana et al. | Influence of variograms in 3D reservoir-modeling outcomes: An example | |
| Holden et al. | Integration of production logs helps to understand heterogeneity of Mishrif reservoir in Rumaila | |
| Liu et al. | Geostatistical modeling for fine reservoir description of Wei2 block of Weicheng oilfield, Dongpu depression, China | |
| Regli et al. | GEOSSAV: a simulation tool for subsurface applications | |
| Penna et al. | Geostatistical seismic inversion and 3D modelling of metric flow units, porosity and permeability in Brazilian presalt reservoir | |
| Li et al. | Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China | |
| Arizabalo et al. | Lacunarity of geophysical well logs in the Cantarell oil field, Gulf of Mexico | |
| US20240302564A1 (en) | A method and system for modeling carbonate cementation in forward depositional models |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request |
Effective date: 20170802 |