WO2012071103A1 - System and method for estimating fluid distribution in a subterranean reservoir - Google Patents

System and method for estimating fluid distribution in a subterranean reservoir Download PDF

Info

Publication number
WO2012071103A1
WO2012071103A1 PCT/US2011/051908 US2011051908W WO2012071103A1 WO 2012071103 A1 WO2012071103 A1 WO 2012071103A1 US 2011051908 W US2011051908 W US 2011051908W WO 2012071103 A1 WO2012071103 A1 WO 2012071103A1
Authority
WO
WIPO (PCT)
Prior art keywords
microporosity
capillary pressure
pressure data
macroporosity
water saturation
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.)
Ceased
Application number
PCT/US2011/051908
Other languages
French (fr)
Inventor
Scott Hanson
Jeroen Brantjes
Katherine Trigg
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chevron USA Inc
Original Assignee
Chevron USA Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chevron USA Inc filed Critical Chevron USA Inc
Priority to CA2815115A priority Critical patent/CA2815115A1/en
Priority to EA201390754A priority patent/EA201390754A1/en
Priority to CN2011800546496A priority patent/CN103210180A/en
Priority to EP11843660.9A priority patent/EP2643713A4/en
Priority to AU2011332287A priority patent/AU2011332287B2/en
Publication of WO2012071103A1 publication Critical patent/WO2012071103A1/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Definitions

  • the present invention relates generally to methods and systems for estimating fluid distribution in a subterranean reservoir, and in particular methods and systems for calculating water saturation within the macroporosity and microporosity of the rock formations to estimate fluid distribution in a subterranean reservoir.
  • microporosity will hold most of the water in the formation and most of the water will not flow out of the microporosity, however some gas could be produced and included in the Gas In Place values so as to not underestimate reserves.
  • alternate recovery techniques such as horizontal drilling and hydraulic fracturing, could be designed to better recover the oil.
  • a computer- implemented method for estimating fluid distribution in subterranean reservoirs may include receiving capillary pressure data from at least one representative rock sample, then setting a threshold dividing the capillary pressure data representative of microporosity and the capillary pressure data representative of macroporosity.
  • the water saturation in the macroporosity may be determined from the capillary pressure data related to the macroporosity using a saturation height function.
  • the capillary pressure data related to the microporosity may be corrected to have an entry pore value equivalent to the pore size defining the microporosity and the corrected data may be used to determine the water saturation in the microporosity using a saturation height function.
  • the macroporosity water saturation and the microporosity water saturation may then be used to estimate fluid distribution within the subterranean reservoir.
  • the fluid distribution in the reservoir may also be used to estimate hydrocarbon reserves.
  • the present invention may also be practiced as a system including a data source containing capillary pressure data that is input to at least one computer processor configured to execute computer program modules.
  • the computer program modules may include an input module to receive the capillary pressure data, a thresholding module to set a threshold between the capillary pressure data representative of microporosity and the capillary pressure data representative of macroporosity, a correction module to correct or normalize the entry pore value of the capillary pressure data representative of microporosity to a pore throat size defining microporosity, a water saturation module to calculate the water saturation in the microporosity and macroporosity, and a fluid distribution module to estimate the fluid distribution in a subterranean reservoir.
  • the computer program modules may also include a a hydrocarbon reserves module to calculate hydrocarbon reserves and an output module to store or display fluid distribution, hydrocarbon reserves, water saturations, or corrected capillary data.
  • the system may also include a user interface to allow interaction with the computer program modules and/or observe results of the computer program modules.
  • the present invention encompasses an article of manufacture including a computer readable medium having computer readable code on it, which will allow a computer to implement a method for estimating fluid distribution in a subterranean reservoir.
  • the method may include determining the macroporosity water saturation from capillary pressure data representative of the macroporosity, correcting the entry pore value of the capillary pressure data representative of microporosity to a pore throat size defining microporosity, using the corrected data to determine the water saturation in the microporosity, and using the macroporosity water saturation and microporosity water saturation to calculate the fluid distribution in the subterranean reservoir.
  • the method may also include setting a threshold to separate the capillary pressure data representative of macroporosity from the capillary pressure data representative of microporosity and estimating hydrocarbon reserves.
  • Figure 1 is a flowchart illustrating a method in accordance with an embodiment of the invention
  • Figure 2 is a graph of capillary pressure data
  • Figure 3 is a graph of a saturation height function with microporosity data displayed before and after correction.
  • FIG. 4 schematically illustrates a system for performing a method in accordance with an embodiment of the invention.
  • the present invention may be described and implemented in the general context of a system and computer methods to be executed by a computer.
  • Such computer- executable instructions may include programs, routines, objects, components, data structures, and computer software technologies that can be used to perform particular tasks and process abstract data types.
  • Software implementations of the present invention may be coded in different languages for application in a variety of computing platforms and environments. It will be appreciated that the scope and underlying principles of the present invention are not limited to any particular computer software technology.
  • the present invention may be practiced using any one or combination of hardware and software configurations, including but not limited to a system having single and/or multiple processor computers, hand-held devices, programmable consumer electronics, mini-computers, mainframe computers, and the like.
  • the invention may also be practiced in distributed computing environments where tasks are performed by servers or other processing devices that are linked through a one or more data communications network.
  • program modules may be located in both local and remote computer storage media including memory storage devices.
  • the present invention may also be practiced as part of a down-hole sensor or measuring device or as part of a laboratory measuring device.
  • CD pre-recorded disk or other equivalent devices
  • CD may include a computer program storage medium and program means recorded thereon for directing the computer processor to facilitate the implementation and practice of the present invention.
  • Such devices and articles of manufacture also fall within the spirit and scope of the present invention.
  • the present invention relates to determining fluid distribution in subterranean reservoirs by calculating and combining water saturation in macroporosity and microporosity within the reservoir.
  • the inventor has found that by determining water saturation for macroporosity and microporosity separately, fluid distribution within the subterranean reservoir may be accurately modeled.
  • the macroporosity and microporosity in the reservoir can be determined by analyzing capillary pressure data from representative rock samples and the water saturation for each can be calculated using a saturation height function.
  • the capillary pressure data related to the microporosity may be corrected prior to calculating the water saturation.
  • capillary pressure data is received.
  • the capillary pressure data may be acquired by any of a variety of laboratory methods including Mercury Capillary Pressure Injection (MICP) and the air-brine centrifuge method.
  • MICP Mercury Capillary Pressure Injection
  • the capillary pressure data (P c ) can be related to the pore throat radius (r), with the following equation:
  • ⁇ is the interfacial tension in units of dynes/cm
  • ⁇ is the contact angle of the fluids in the rock
  • C is a constant determined for the rock type and is approximately l .
  • a pressure or pore throat size threshold between the macroporosity and microporosity can be set in step 12. In one embodiment, it is set based on a known threshold of pore throat size or capillary pressure. It may also be based on a graph of the capillary pressure data. The distribution of macroporosity and microporosity is observed based on graphing the pore throat size on the x axis and the cumulative porosity from core plugs on the y axis. Microporosity is determined by the lower bimodal distribution and the macroporosity is determined by the higher bimodal distribution.
  • Figure 2 shows a graph of capillary pressure data 200 that is being displayed in terms of the pore throat diameter versus the mercury saturation with the threshold 201 separating the microporosity 205 and the macroporosity 210.
  • the capillary pressure data that is above the pressure threshold or below the pore throat size threshold is related to the microporosity.
  • the present invention processes the capillary pressure data related to the microporosity differently from the capillary pressure data related to the macroporosity.
  • the capillary pressure data related to the microporosity is corrected, or normalized, to reflect that the entry pore value is equivalent to the pore throat size defining the microporosity, which means that the corrected data will reflect that the water saturation in the microporosity is 100% at the pressure or pore throat size threshold. This is done by subtracting the water saturation in the macroporosity and resetting the water saturation in the microporosity back to 100% as indicated in the following equation:
  • Sw (l-0— Swmic co ) + Swmic Eqn. 2 where S Wmic is the saturation above the microporosity cut off and S WmiCco is the saturation at the microporosity cut off.
  • the water saturation of the microporosity may be calculated at step 16.
  • this calculation may be done using a saturation height function such as a Leverett J-function (J) as shown here: where P c is the capillary pressure measurement in units of pounds per square inch (psi), ⁇ is the interfacial tension in units of dynes/cm, K is the permeability and ⁇ is the porosity of the rock Porosity may be obtained from log data.
  • J Leverett J-function
  • one method that may be used is to obtain the permeability core data from MICP data, helium or air injection data of the core plugs, then plot the log porosity data against the core permeability data. This provides a line from which an equation may be derived. That equation is then used to create permeability data for the entire depth where porosity data is known.
  • the J-function is a dimensionless value for a rock- fluid system that accounts for the effects of fluid and pore geometry.
  • the capillary pressure data below the capillary pressure threshold or above the pore throat size threshold determined in step 12 is related to the macroporosity and can be used directly to calculate the water saturation in the macroporosity in step 18. In one embodiment, this calculation may also be done using a saturation height function such as a Leverett J-function as shown in Equation 3 and Equation 4.
  • Figure 3 shows a water saturation vs. J function graph 300 with regions of microporosity 305 and macroporosity 310 indicated.
  • the data displayed includes the data related to the macroporosity and the data related to the microporosity without the correction 315 in step 14 and the data related to the microporosity after the correction 320 done in step 14. It will be appreciated that step 14 and step 16 are usually done in the order they are shown and step 18 may be done before, after, or concurrently with steps 14 and 16.
  • the water saturation in the macroporosity and the water saturation in the microporosity are combined to calculate the fluid distribution in the reservoir.
  • a map of the macroporosity and microporosity of the reservoir is needed.
  • the macroporosity and microporosity can be determined, for example, from MICP measurements and then mapped throughout the reservoir during formation evaluation and subsurface modeling.
  • One technique of populating a reservoir model with macroporosity and microporosity would be to use the relationship between acoustic impedance and each porosity component. Acoustic impedance may be available throughout the reservoir based on seismic data.
  • the water saturation of the macroporosity and microporosity are then added to the subsurface model and are used to make the fluid distribution model.
  • hydrocarbon reserves can be estimated. For example, a one acre foot contains 7,758 barrels of oil. If the porosity is 30%, then the acre foot contains 30% of that. If the porosity is 30% and the water saturation is 50%, it means that half of the porosity is filled up with water, so there is 15% of available porosity for oil. Using macro- and microporosity and water saturation for each of those porosity components, it would be possible to further determine hydrocarbon distribution within a reservoir.
  • the system includes a data source 42 which may include, among others, a data storage device or memory.
  • the stored capillary pressure data may be made available to a processor 44, such as a programmable general purpose computer.
  • the processor 44 is configured to execute a correction module to correct the capillary pressure data corresponding to the microporosity; a water saturation module to calculate a water saturation for the macroporosity and microporosity; a fluid distribution module to estimate the fluid distribution in the subterranean reservoir; and an output module to store or display at least one of the following: corrected capillary pressure data, the water saturation of the microporosity, the water saturation of the macroporosity, and/or the fluid distribution.
  • the processor 44 may include interface components such as user interface 46, and is used to implement the above-described transforms in accordance with embodiments of the invention.
  • the user interface 46 may be used both to display data and processed data products and to allow the user to select among options for implementing aspects of the method.
  • the water saturations and fluid distribution computed on the processor 44 may be displayed on the user interface 46, stored on the data storage device or memory 42, or both displayed and stored.
  • the fluid distribution estimated here may be further used to calculate hydrocarbon reserves.

Landscapes

  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Measuring Fluid Pressure (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

A system and method for determining fluid distribution in subterranean reservoirs including determining a water saturation in macroporosity from the capillary pressure data representative of the macroporosity using a saturation height function, correcting capillary pressure data representative of microporosity to have an entry pore value equivalent to a pore size defining the microporosity, determining a water saturation in the microporosity from the corrected capillary pressure data representative of the microporosity, and using the macroporosity water saturation and the microporosity water saturation to estimate fluid distribution within the subterranean reservoir. The system and method may also include the estimation of hydrocarbon reserves.

Description

SYSTEM AND METHOD FOR ESTIMATING FLUID DISTRIBUTION IN A
SUBTERRANEAN RESERVOIR
FIELD OF THE INVENTION
[0001] The present invention relates generally to methods and systems for estimating fluid distribution in a subterranean reservoir, and in particular methods and systems for calculating water saturation within the macroporosity and microporosity of the rock formations to estimate fluid distribution in a subterranean reservoir.
BACKGROUND OF THE INVENTION
[0002] Calculating fluid distribution in subterranean reservoirs is an important step in determining potential hydrocarbon reserves. As hydrocarbon exploration and production moves to unconventional reservoirs such as complex carbonate formations and shale gas formations, the calculation of the fluid distribution becomes more difficult because of the varying porosity of the rocks in the formation. In particular, microporosity within the formation may cause fluid distribution calculations to be inaccurate.
[0003] In gas reservoirs, microporosity will hold most of the water in the formation and most of the water will not flow out of the microporosity, however some gas could be produced and included in the Gas In Place values so as to not underestimate reserves. In oil reservoirs, if a significant amount of oil is held in microporosity, then alternate recovery techniques, such as horizontal drilling and hydraulic fracturing, could be designed to better recover the oil.
[0004] Existing methods for calculating fluid distribution in subterranean reservoirs do not take into account the differences between the fluids in macropores and the fluids in micropores. These existing methods may use an average porosity that combines the microporosity and macroporosity or simply ignore the microporosity. When the reservoir under consideration has significant microporosity, such as oil shale reservoirs and shale gas reservoirs, the existing methods may not accurately calculate the fluid distribution. SUMMARY OF THE INVENTION
[0005] Described herein are implementations of various approaches for determining fluid distribution in subterranean reservoirs and, more particularly, for using water saturation in macroporosity and water saturation in microporosity to determine fluid distribution in subterranean reservoirs. According to one aspect of the present invention, a computer- implemented method for estimating fluid distribution in subterranean reservoirs may include receiving capillary pressure data from at least one representative rock sample, then setting a threshold dividing the capillary pressure data representative of microporosity and the capillary pressure data representative of macroporosity. The water saturation in the macroporosity may be determined from the capillary pressure data related to the macroporosity using a saturation height function. The capillary pressure data related to the microporosity may be corrected to have an entry pore value equivalent to the pore size defining the microporosity and the corrected data may be used to determine the water saturation in the microporosity using a saturation height function. The macroporosity water saturation and the microporosity water saturation may then be used to estimate fluid distribution within the subterranean reservoir. The fluid distribution in the reservoir may also be used to estimate hydrocarbon reserves.
[0006] The present invention may also be practiced as a system including a data source containing capillary pressure data that is input to at least one computer processor configured to execute computer program modules. The computer program modules may include an input module to receive the capillary pressure data, a thresholding module to set a threshold between the capillary pressure data representative of microporosity and the capillary pressure data representative of macroporosity, a correction module to correct or normalize the entry pore value of the capillary pressure data representative of microporosity to a pore throat size defining microporosity, a water saturation module to calculate the water saturation in the microporosity and macroporosity, and a fluid distribution module to estimate the fluid distribution in a subterranean reservoir. The computer program modules may also include a a hydrocarbon reserves module to calculate hydrocarbon reserves and an output module to store or display fluid distribution, hydrocarbon reserves, water saturations, or corrected capillary data. The system may also include a user interface to allow interaction with the computer program modules and/or observe results of the computer program modules. [0007] In addition, the present invention encompasses an article of manufacture including a computer readable medium having computer readable code on it, which will allow a computer to implement a method for estimating fluid distribution in a subterranean reservoir. The method may include determining the macroporosity water saturation from capillary pressure data representative of the macroporosity, correcting the entry pore value of the capillary pressure data representative of microporosity to a pore throat size defining microporosity, using the corrected data to determine the water saturation in the microporosity, and using the macroporosity water saturation and microporosity water saturation to calculate the fluid distribution in the subterranean reservoir. The method may also include setting a threshold to separate the capillary pressure data representative of macroporosity from the capillary pressure data representative of microporosity and estimating hydrocarbon reserves.
[0008] The above summary section is provided to introduce a selection of concepts in a simplified form that are further described below in the detailed description section. The summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used to limit the scope of the claimed subject matter. Furthermore, the claimed subject matter is not limited to implementations that solve any or all disadvantages noted in any part of this disclosure.
BRIEF DESCRIPTION OF THE DRAWINGS
[0009] These and other features of the present invention will become better understood with regard to the following description, claims and accompanying drawings where:
[0010] Figure 1 is a flowchart illustrating a method in accordance with an embodiment of the invention;
[0011] Figure 2 is a graph of capillary pressure data;
[0012] Figure 3 is a graph of a saturation height function with microporosity data displayed before and after correction; and
[0013] Figure 4 schematically illustrates a system for performing a method in accordance with an embodiment of the invention. DETAILED DESCRIPTION OF THE INVENTION
[0014] The present invention may be described and implemented in the general context of a system and computer methods to be executed by a computer. Such computer- executable instructions may include programs, routines, objects, components, data structures, and computer software technologies that can be used to perform particular tasks and process abstract data types. Software implementations of the present invention may be coded in different languages for application in a variety of computing platforms and environments. It will be appreciated that the scope and underlying principles of the present invention are not limited to any particular computer software technology.
[0015] Moreover, those skilled in the art will appreciate that the present invention may be practiced using any one or combination of hardware and software configurations, including but not limited to a system having single and/or multiple processor computers, hand-held devices, programmable consumer electronics, mini-computers, mainframe computers, and the like. The invention may also be practiced in distributed computing environments where tasks are performed by servers or other processing devices that are linked through a one or more data communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices. The present invention may also be practiced as part of a down-hole sensor or measuring device or as part of a laboratory measuring device.
[0016] Also, an article of manufacture for use with a computer processor, such as a
CD, pre-recorded disk or other equivalent devices, may include a computer program storage medium and program means recorded thereon for directing the computer processor to facilitate the implementation and practice of the present invention. Such devices and articles of manufacture also fall within the spirit and scope of the present invention.
[0017] Referring now to the drawings, embodiments of the present invention will be described. The invention can be implemented in numerous ways, including for example as a system (including a computer processing system), a method (including a computer implemented method), an apparatus, a computer readable medium, a computer program product, a graphical user interface, a web portal, or a data structure tangibly fixed in a computer readable memory. Several embodiments of the present invention are discussed below. The appended drawings illustrate only typical embodiments of the present invention and therefore are not to be considered limiting of its scope and breadth.
[0018] The present invention relates to determining fluid distribution in subterranean reservoirs by calculating and combining water saturation in macroporosity and microporosity within the reservoir. The inventor has found that by determining water saturation for macroporosity and microporosity separately, fluid distribution within the subterranean reservoir may be accurately modeled. The macroporosity and microporosity in the reservoir can be determined by analyzing capillary pressure data from representative rock samples and the water saturation for each can be calculated using a saturation height function. The capillary pressure data related to the microporosity may be corrected prior to calculating the water saturation.
[0019] In this regard, an example of a method 100 in accordance with the present invention is illustrated in the flowchart of Figure 1. At step 10, capillary pressure data is received. As will be appreciated, the capillary pressure data may be acquired by any of a variety of laboratory methods including Mercury Capillary Pressure Injection (MICP) and the air-brine centrifuge method.
[0020] The capillary pressure data (Pc) can be related to the pore throat radius (r), with the following equation:
2*a*cos 0 _ 1
r = Eqn. I where σ is the interfacial tension in units of dynes/cm, Θ is the contact angle of the fluids in the rock, and C is a constant determined for the rock type and is approximately l . From this equation, it is clear that the large pore throats that may be considered macroporosity in the rock are related to low pressure measurements and, conversely, high pressure measurements will be related to small pore throats or microporosity. Microporosity may be defined as having pore throat radii of less than l μιη, typically less than 0.75 μιη, and often less than 0.5 μιη. Macroporosity often has pore throat radii greater than 0.5 μιη but may have pore throat radii as small as O. l μιη.
[0021] From the capillary pressure data, a pressure or pore throat size threshold between the macroporosity and microporosity can be set in step 12. In one embodiment, it is set based on a known threshold of pore throat size or capillary pressure. It may also be based on a graph of the capillary pressure data. The distribution of macroporosity and microporosity is observed based on graphing the pore throat size on the x axis and the cumulative porosity from core plugs on the y axis. Microporosity is determined by the lower bimodal distribution and the macroporosity is determined by the higher bimodal distribution.
[0022] Figure 2 shows a graph of capillary pressure data 200 that is being displayed in terms of the pore throat diameter versus the mercury saturation with the threshold 201 separating the microporosity 205 and the macroporosity 210. There are eight curves, 220- 227, representing capillary pressure data from eight representative rock samples. The capillary pressure data that is above the pressure threshold or below the pore throat size threshold is related to the microporosity. The present invention processes the capillary pressure data related to the microporosity differently from the capillary pressure data related to the macroporosity. Referring again to Figure 1, at step 14 the capillary pressure data related to the microporosity is corrected, or normalized, to reflect that the entry pore value is equivalent to the pore throat size defining the microporosity, which means that the corrected data will reflect that the water saturation in the microporosity is 100% at the pressure or pore throat size threshold. This is done by subtracting the water saturation in the macroporosity and resetting the water saturation in the microporosity back to 100% as indicated in the following equation:
Sw = (l-0— Swmicco) + Swmic Eqn. 2 where SWmic is the saturation above the microporosity cut off and SWmiCco is the saturation at the microporosity cut off.
[0023] After the capillary pressure data related to the microporosity is corrected, the water saturation of the microporosity may be calculated at step 16. In one embodiment, this calculation may be done using a saturation height function such as a Leverett J-function (J) as shown here:
Figure imgf000008_0001
where Pc is the capillary pressure measurement in units of pounds per square inch (psi), σ is the interfacial tension in units of dynes/cm, K is the permeability and φ is the porosity of the rock Porosity may be obtained from log data. To determine permeability, one method that may be used is to obtain the permeability core data from MICP data, helium or air injection data of the core plugs, then plot the log porosity data against the core permeability data. This provides a line from which an equation may be derived. That equation is then used to create permeability data for the entire depth where porosity data is known.
The J-function is a dimensionless value for a rock- fluid system that accounts for the effects of fluid and pore geometry. Once the J-function is calculated, a relationship between J and Sw must be determined. This may be done, for example, using linear regression analysis of logio(Sw) and logio(J) from air-brine capillary pressure data or from MICP data that has been corrected to simulate data from an air-brine system. After a regression line is established, a slope (b) and intercept (a) is noted from cross-plot data as input into a SW(J) function as: logio Sw (J = b * log10 (/) + a Eqn. 4
[0024] The capillary pressure data below the capillary pressure threshold or above the pore throat size threshold determined in step 12 is related to the macroporosity and can be used directly to calculate the water saturation in the macroporosity in step 18. In one embodiment, this calculation may also be done using a saturation height function such as a Leverett J-function as shown in Equation 3 and Equation 4. Figure 3 shows a water saturation vs. J function graph 300 with regions of microporosity 305 and macroporosity 310 indicated. The data displayed includes the data related to the macroporosity and the data related to the microporosity without the correction 315 in step 14 and the data related to the microporosity after the correction 320 done in step 14. It will be appreciated that step 14 and step 16 are usually done in the order they are shown and step 18 may be done before, after, or concurrently with steps 14 and 16.
[0025] Referring again to Figure 1, at step 20 the water saturation in the macroporosity and the water saturation in the microporosity are combined to calculate the fluid distribution in the reservoir. For this step, a map of the macroporosity and microporosity of the reservoir is needed. The macroporosity and microporosity can be determined, for example, from MICP measurements and then mapped throughout the reservoir during formation evaluation and subsurface modeling. One technique of populating a reservoir model with macroporosity and microporosity would be to use the relationship between acoustic impedance and each porosity component. Acoustic impedance may be available throughout the reservoir based on seismic data. The water saturation of the macroporosity and microporosity are then added to the subsurface model and are used to make the fluid distribution model. Once a fluid distribution model is established, hydrocarbon reserves can be estimated. For example, a one acre foot contains 7,758 barrels of oil. If the porosity is 30%, then the acre foot contains 30% of that. If the porosity is 30% and the water saturation is 50%, it means that half of the porosity is filled up with water, so there is 15% of available porosity for oil. Using macro- and microporosity and water saturation for each of those porosity components, it would be possible to further determine hydrocarbon distribution within a reservoir.
[0026] A system 400 for performing the method is schematically illustrated in Figure
4. The system includes a data source 42 which may include, among others, a data storage device or memory. The stored capillary pressure data may be made available to a processor 44, such as a programmable general purpose computer. The processor 44 is configured to execute a correction module to correct the capillary pressure data corresponding to the microporosity; a water saturation module to calculate a water saturation for the macroporosity and microporosity; a fluid distribution module to estimate the fluid distribution in the subterranean reservoir; and an output module to store or display at least one of the following: corrected capillary pressure data, the water saturation of the microporosity, the water saturation of the macroporosity, and/or the fluid distribution. The processor 44 may include interface components such as user interface 46, and is used to implement the above-described transforms in accordance with embodiments of the invention. The user interface 46 may be used both to display data and processed data products and to allow the user to select among options for implementing aspects of the method. The water saturations and fluid distribution computed on the processor 44 may be displayed on the user interface 46, stored on the data storage device or memory 42, or both displayed and stored. The fluid distribution estimated here may be further used to calculate hydrocarbon reserves.
[0027] While in the foregoing specification this invention has been described in relation to certain preferred embodiments thereof, and many details have been set forth for purpose of illustration, it will be apparent to those skilled in the art that the invention is susceptible to alteration and that certain other details described herein can vary considerably without departing from the basic principles of the invention. In addition, it should be appreciated that structural features or method steps shown or described in any one embodiment herein can be used in other embodiments as well.

Claims

WHAT IS CLAIMED IS:
A computer-implemented method for estimating fluid distribution in a subterranean reservoir comprising:
a. determining a macroporosity water saturation from capillary pressure data representative of a macroporosity;
b. correcting an entry pore value of capillary pressure data representative of a microporosity;
c. determining a microporosity water saturation from the corrected capillary pressure data representative of the microporosity; and
d. using the macroporosity water saturation and the microporosity water
saturation to estimate a fluid distribution within the subterranean reservoir.
The method of claim 1 , wherein the entry pore value of the capillary pressure data representative of the microporosity is set to a pore throat size defining the microporosity.
The method of claim 1 , further comprising estimating hydrocarbon reserves from the fluid distribution.
The method of claim 1 , further comprising analysis of capillary pressure data to determine a threshold dividing the capillary pressure data representative of the microporosity and the capillary pressure data representative of the macroporosity.
The method of claim 1 , wherein the capillary pressure data is from an air-brine centrifuge method.
The method of claim 1 , wherein the capillary pressure data is from a Mercury Capillary Pressure Injection.
The method of claim 6, wherein the capillary pressure data from the Mercury Capillary Pressure Injection is transformed to simulate capillary pressure data from the air-brine centrifuge method. 8) The method of claim 1 , wherein the macroporosity water saturation is calculated using a saturation height function.
9) The method of claim 8, wherein the saturation height function is a Leverett J- function.
10) The method of claim 1, wherein the microporosity water saturation is calculated using a saturation height function.
11) The method of claim 10, wherein the saturation height function is a Leverett J- function.
12) A system for estimating fluid distribution in a subterranean reservoir comprising: a. a data source containing capillary pressure data;
b. at least one computer processor being configured to communicate with the data source and to execute computer program modules, the computer modules comprising:
an input module to receive the capillary pressure data from the data source;
a correction module to correct the capillary pressure data
representative of a microporosity;
a water saturation module to calculate a water saturation for macroporosity and a water saturation for the microporosity; and a fluid distribution module to estimate a fluid distribution in the subterranean reservoir.
13) The system of claim 12, wherein the correction module sets the entry pore value of the capillary pressure data representative of the microporosity to a pore throat size defining the microporosity. 14) The system of claim 12, further comprising an output module to store or display at least one of the corrected capillary data, the water saturation for the macroporosity, the water saturation for the microporosity, and/or the fluid distribution.
15) The system of claim 12, further comprising a user interface device to allow
interaction with the computer program modules and/or observe results of the computer program modules.
16) The system of claim 12, further comprising a hydrocarbon reserves module for estimating hydrocarbon reserves from the fluid distribution.
17) The system of claim 12, further comprising a thresholding module to graph the capillary pressure data and determine a threshold between the capillary pressure data representative of the microporosity and the capillary pressure data representative of the macroporosity.
18) An article of manufacture comprising a computer readable medium having a
computer readable code embodied therein, the computer readable program code adapted to be executed to implement a method for estimating fluid distribution in a subterranean reservoir, the method comprising:
a. determining a macroporosity water saturation from capillary pressure data representative of a macroporosity;
b. correcting an entry pore value of capillary pressure data representative of a microporosity;
c. determining a microporosity water saturation from the corrected capillary pressure data representative of the microporosity; and
d. using the macroporosity water saturation and the microporosity water
saturation to estimate a fluid distribution within the subterranean reservoir.
19) The method of claim 18, wherein the correcting the entry value of the capillary pressure data representative of the microporosity sets the entry pore value to the pore throat size defining microporosity. 20) The method of claim 18, further comprising estimating hydrocarbon reserves from the fluid distribution.
21) The method of claim 18, further comprising setting a threshold dividing the capillary pressure data representative of the microporosity and the capillary pressure data representative of the macroporosity.
PCT/US2011/051908 2010-11-24 2011-09-16 System and method for estimating fluid distribution in a subterranean reservoir Ceased WO2012071103A1 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CA2815115A CA2815115A1 (en) 2010-11-24 2011-09-16 System and method for estimating fluid distribution in a subterranean reservoir
EA201390754A EA201390754A1 (en) 2010-11-24 2011-09-16 SYSTEM AND METHOD FOR ASSESSING FLUID DISTRIBUTION IN THE UNDERGROUND COLLECTOR
CN2011800546496A CN103210180A (en) 2010-11-24 2011-09-16 System and method for estimating fluid distribution in a subterranean reservoir
EP11843660.9A EP2643713A4 (en) 2010-11-24 2011-09-16 System and method for estimating fluid distribution in a subterranean reservoir
AU2011332287A AU2011332287B2 (en) 2010-11-24 2011-09-16 System and method for estimating fluid distribution in a subterranean reservoir

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US12/954,374 US8645070B2 (en) 2010-11-24 2010-11-24 System and method for estimating fluid distribution in a subterranean reservoir
US12/954,374 2010-11-24

Publications (1)

Publication Number Publication Date
WO2012071103A1 true WO2012071103A1 (en) 2012-05-31

Family

ID=46065120

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2011/051908 Ceased WO2012071103A1 (en) 2010-11-24 2011-09-16 System and method for estimating fluid distribution in a subterranean reservoir

Country Status (7)

Country Link
US (1) US8645070B2 (en)
EP (1) EP2643713A4 (en)
CN (1) CN103210180A (en)
AU (1) AU2011332287B2 (en)
CA (1) CA2815115A1 (en)
EA (1) EA201390754A1 (en)
WO (1) WO2012071103A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110029989A (en) * 2018-01-11 2019-07-19 中国石油化工股份有限公司 A kind of unconventional oil and gas recovery percent of reserves calculation method and system
US11360004B2 (en) 2017-09-15 2022-06-14 Korea Institute Of Geoscience And Mineral Resourues Shale gas extracting device and extracting method therefor

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2862951C (en) 2011-10-21 2018-10-30 Saudi Arabian Oil Company Methods, computer readable medium, and apparatus for determining well characteristics and pore architecture utilizing conventional well logs
US20140088878A1 (en) * 2012-09-27 2014-03-27 Jinhong Chen Isotherm and gas-in-place estimation considering capillary condensation in shale gas reservoir
WO2015021088A1 (en) * 2013-08-06 2015-02-12 Schlumberger Canada Limited Methods for determining a saturation-height function in oil and gas reservoirs
US10495774B2 (en) * 2014-01-13 2019-12-03 Schlumberger Technology Corporation Method for estimating irreducible water saturation from mercury injection capillary pressure
FR3036820B1 (en) * 2015-06-01 2021-12-31 Services Petroliers Schlumberger MODELING OF OIL FIELD RESERVOIR SATURATION AND PERMEABILITY
US10280722B2 (en) 2015-06-02 2019-05-07 Baker Hughes, A Ge Company, Llc System and method for real-time monitoring and estimation of intelligent well system production performance
FR3038408B1 (en) * 2015-06-30 2017-08-25 Services Petroliers Schlumrberger MODELING THE SATURATION AND PERMEABILITY OF PETROLEUM FIELD RESERVOIR
US10552553B2 (en) 2015-08-17 2020-02-04 Saudi Arabian Oil Company Capillary pressure analysis for petrophysical statistical modeling
CN106405679B (en) * 2016-11-09 2018-05-01 中国地质大学（北京） A kind of Quantitative Monitoring subsurface flow to well device and pull out the well lantern ring
AU2018212812A1 (en) 2017-01-26 2019-08-15 Dassault Systemes Simulia Corp. Multi-phase flow visualizations based on fluid occupation time
US11714040B2 (en) * 2018-01-10 2023-08-01 Dassault Systemes Simulia Corp. Determining fluid flow characteristics of porous mediums
US11530598B2 (en) 2018-08-21 2022-12-20 Dassault Systemes Simulia Corp. Determination of oil removed by gas via miscible displacement in reservoir rock
CN109709301B (en) * 2018-11-30 2021-09-28 中国石油天然气股份有限公司 Method, device and system for determining distribution of fracture pore type tight sandstone reservoir
CN111561312B (en) * 2020-04-21 2023-09-26 中国石油天然气股份有限公司 Saturation calculation method and device based on conventional well logging data
US11847391B2 (en) 2020-06-29 2023-12-19 Dassault Systemes Simulia Corp. Computer system for simulating physical processes using surface algorithm
US11907625B2 (en) 2020-12-29 2024-02-20 Dassault Systemes Americas Corp. Computer simulation of multi-phase and multi-component fluid flows including physics of under-resolved porous structures
CN112733388A (en) * 2021-01-26 2021-04-30 新奥数能科技有限公司 Method, apparatus, electronic device and medium for determining saturation temperature of water vapor
CN112966380B (en) * 2021-03-10 2021-12-21 东北石油大学 Method for determining steep slope sand body type gathering ridge oil and gas reserve
CN113740515B (en) * 2021-11-05 2022-03-01 中国科学院地质与地球物理研究所 Method for comprehensively representing occurrence lower limit critical condition of deep sea phase carbonate reservoir
US12578500B2 (en) 2022-03-10 2026-03-17 Saudi Arabian Oil Company Hydrocarbon reservoir saturation logging
US20250044473A1 (en) * 2023-07-31 2025-02-06 Saudi Arabian Oil Company Method of identifying a productive interval along a well
CN120847898B (en) * 2025-04-23 2026-03-20 辽宁工程技术大学 An Adaptive Weighted Fusion Method for Inverting Terrestrial Water Storage Anomalies

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4926128A (en) * 1989-02-13 1990-05-15 Mobil Oil Corporation Method for utilizing measured resistivities of porous rock under differing fluid saturations to identify fluid distribution equilibrium
US20090248309A1 (en) * 2008-03-28 2009-10-01 Schlumberger Technology Corporation Evaluating a reservoir formation
WO2009138934A1 (en) * 2008-05-15 2009-11-19 Koninklijke Philips Electronics N.V. Method and system for detecting a fluid distribution in an object of interest
US20100277167A1 (en) * 2008-01-30 2010-11-04 Baker Hughes Incorporated Method and Apparatus for Determining Multiscale Similarity Between NMR Measurements and a Reference Well Log

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5621169A (en) * 1994-01-18 1997-04-15 Restech, Inc. Method for determining hydrocarbon/water contact level for oil and gas wells
US6484102B1 (en) * 2000-08-24 2002-11-19 Digital Formation Inc. System for evaluating fluid distributions of subsurface reservoirs
WO2003025625A1 (en) 2001-09-19 2003-03-27 Halliburton Energy Services, Inc. Method and system for using conventional core data to calibrate bound water volumes derived from nmr logs
US7277795B2 (en) 2004-04-07 2007-10-02 New England Research, Inc. Method for estimating pore structure of porous materials and its application to determining physical properties of the materials
US7567079B2 (en) * 2007-06-08 2009-07-28 University Of New Brunswick Methods suitable for measuring capillary pressure and relative permeability curves of porous rocks
US8676556B2 (en) 2007-11-27 2014-03-18 Exxonmobil Upstream Research Company Method for determining the properties of hydrocarbon reservoirs from geophysical data
CN101413388A (en) * 2008-12-02 2009-04-22 大庆油田有限责任公司 Method for obtaining oil-water common-layer original oil-containing saturation degree and method for estimating non-test oil-water common-layer original oil-containing saturation degree
CN101892837B (en) * 2010-04-29 2013-03-20 中国石油天然气股份有限公司 Formation factor determining method and oil saturation determining method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4926128A (en) * 1989-02-13 1990-05-15 Mobil Oil Corporation Method for utilizing measured resistivities of porous rock under differing fluid saturations to identify fluid distribution equilibrium
US20100277167A1 (en) * 2008-01-30 2010-11-04 Baker Hughes Incorporated Method and Apparatus for Determining Multiscale Similarity Between NMR Measurements and a Reference Well Log
US20090248309A1 (en) * 2008-03-28 2009-10-01 Schlumberger Technology Corporation Evaluating a reservoir formation
WO2009138934A1 (en) * 2008-05-15 2009-11-19 Koninklijke Philips Electronics N.V. Method and system for detecting a fluid distribution in an object of interest

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP2643713A4 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11360004B2 (en) 2017-09-15 2022-06-14 Korea Institute Of Geoscience And Mineral Resourues Shale gas extracting device and extracting method therefor
CN110029989A (en) * 2018-01-11 2019-07-19 中国石油化工股份有限公司 A kind of unconventional oil and gas recovery percent of reserves calculation method and system
CN110029989B (en) * 2018-01-11 2021-11-02 中国石油化工股份有限公司 Unconventional oil and gas extraction degree calculation method and system

Also Published As

Publication number Publication date
EA201390754A1 (en) 2013-09-30
AU2011332287A1 (en) 2013-03-21
EP2643713A4 (en) 2017-11-15
US20120130639A1 (en) 2012-05-24
EP2643713A1 (en) 2013-10-02
CN103210180A (en) 2013-07-17
AU2011332287B2 (en) 2013-10-17
CA2815115A1 (en) 2012-05-31
US8645070B2 (en) 2014-02-04

Similar Documents

Publication Publication Date Title
US8645070B2 (en) System and method for estimating fluid distribution in a subterranean reservoir
CN104879103B (en) A method for analyzing the effect of layered water injection
US8332155B2 (en) System and method for hydrocarbon gas pay zone characterization in a subterranean reservoir
CN112287508B (en) Method and device for determining oil-gas-water interface based on equivalent density of formation pressure
US20180032356A1 (en) Enhanced Oil Recovery (EOR) Chemical Coreflood Simulation Study Workflow
CA2809969C (en) System and method for sweet zone identification in shale gas reservoirs
FR3043481A1 (en)
EP2103964A1 (en) Method for calculating the density of the fractures in a rocky medium
Thomas et al. The fate of submarine fresh groundwater reservoirs at the New Jersey shelf, USA
US10787902B2 (en) Method and system for correcting a capillary pressure curve
CN110632658A (en) Analysis method and device for lateral sealing of faults
Wang et al. Numerical analysis of the effect of heterogeneity on CO2 dissolution enhanced by gravity-driven convection
CN115144929A (en) Oil and gas reservoir reserves determination method, device, storage medium and electronic equipment
CN112394416A (en) Heterogeneous fracture control reservoir prediction method and device
CN118349801A (en) Determination method, device, equipment and medium for karst action influence depth
CN115075800A (en) Hydrocarbon reservoir identification method, apparatus, medium and product
US20230259670A1 (en) Permeability modeling in a reservoir simulation model using dynamic pressure transient analysis
CN109425904B (en) Method for acquiring pore pressure of carbonate rock stratum
Potocki Resolving differences between core and welltest permeability in basal Colorado sandstones, Canada: The role of rock heterogeneity, relative permeability and formation damage
WO2025183686A1 (en) Stochastic workflow for geological fault seal analysis in carbonate-hosted reservoirs
CN120906541A (en) Method and related equipment for measuring saturation of natural gas hydrate and associated gas co-storage layer
Han et al. Investigating the geochemical alterations in an aquifer due to long-term sequestration of CO2 using time-lapse seismic information
Deshenenkov et al. The Estimation of Granular Reservoirs Relative Permeability with Log Data
CN110017126A (en) Method and device for determining wellbore orientation of horizontal wells in fractured reservoirs
CN110159265A (en) Method and apparatus for determining oil saturation of a heterogeneous reservoir

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 11843660

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2011332287

Country of ref document: AU

Date of ref document: 20110916

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2011843660

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2815115

Country of ref document: CA

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 201390754

Country of ref document: EA

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112013006925

Country of ref document: BR

REG Reference to national code

Ref country code: BR

Ref legal event code: B01E

Ref document number: 112013006925

Country of ref document: BR

ENPW Started to enter national phase and was withdrawn or failed for other reasons

Ref document number: 112013006925

Country of ref document: BR