WO2014207197A1

WO2014207197A1 - Method to generate scenarios of hydrocarbon reservoirs based on limited amount of information on a target hydrocarbon reservoir

Info

Publication number: WO2014207197A1
Application number: PCT/EP2014/063689
Authority: WO
Inventors: Sonia Embid Droz; Rubén RODRÍGUEZ TORRADO; Giorgio De Paola; Mohamed Hegazy; David Echeverría Ciaurri; Bruno FLACH; Ulisses Mello
Original assignee: Repsol, S.A.; International Business Machines Corporation
Priority date: 2013-06-28
Filing date: 2014-06-27
Publication date: 2014-12-31
Also published as: CN105453125A; CA2914465A1; BR112015032537A2; EP3014560A1

Abstract

The present invention is related to generating scenarios of hydrocarbon reservoirs based on limited amount of information on a target hydrocarbon reservoir, and more particularly to automatically supplying missing parameters and an uncertainty associated with each supplied parameter allowing to valuating the target hydrocarbon reservoir.

Description

METHOD TO GENERATE SCENARIOS OF HYDROCARBON RESERVOIRS

BASED ON LIMITED AMOUNT OF INFORMATION ON A TARGET HYDROCARBON RESERVOIR

DESCRIPTION

OBJECT OF THE INVENTION The present invention is related to generating scenarios of hydrocarbon reservoirs based on limited amount of information on a target hydrocarbon reservoir, and more particularly to automatically supplying missing parameters and an uncertainty associated with each supplied parameter allowing to valuating the target hydrocarbon reservoir.

PRIOR ART

Determining whether investing in a new hydrocarbon reservoir candidate is a good business decision depends on the inherent value of the reservoir. Factors determining the inherent value of the reservoir include, for example, the total amount of material that is ultimately recoverable from each new hydrocarbon reservoir (production potential), market prices (oil and/or natural gas prices) and the cost of recovering that material, or capture difficulty. Until the material is actually recovered, however, that inherent value can only be estimated from different, primarily known, reservoir prop- erties.

However, investment opportunities frequently include a number of candidate reservoirs with spotty available information on each. Even with a dearth of available information, however, it is imperative to assess the value of each candidate accurately. This assessment may be even more complex when types of available information vary from candidate to candidate. Arriving at an optimal project portfolio value requires a uniform assessment that consistently evaluates all alternatives uniformly.

Some investment candidate estimates have been based on identifying existing reser- voirs with properties that match or closely match (i.e., are similar to) the new, candidate reservoir. The closest existing reservoirs are known as "analogous reservoirs." Typically, one or more experts (e.g., geologists and reservoir engineers) identified and selected analogous reservoirs, based solely on experience and known or available investment candidate properties. Relevant information of such analogous reservoirs are stored in data bases having records with information on volumetric, facies and properties distribution, wherein all records not necessarily have all properties.

The specialized literature shows continuously increased interest in analogues and its predictions. Sophisticated data mining and machine learning algorithms allow estima- tion of properties and such estimations also define the degree of uncertainty which can be reliably used within geostatistical workflows. Example can be found on the usage of reservoir analogous in the context of reservoir drive mechanism, recovery factor or reserve estimation and classification. Insufficient information, however, can make selecting analogous reservoirs guesswork at best, and make estimating the value error prone and uncertain. A mis-valuation could lead to wasted resources, e.g., from passing on an undervalued reservoir to exploit an overvalued reservoir. Missing parameters increase uncertainty and the likelihood of a mis-valuation.

Thus, there is a need for generating scenarios compatible with the limited amount of information available on a target hydrocarbon reservoir, while including an associated uncertainty for the assessment; and more particularly, for estimating geological and petrophysical properties.

SUMARY OF THE INVENTION

The present invention solves the problems providing a method to generate scenarios of hydrocarbon reservoirs when only limited amount of information on a target is available. Ranking reservoirs with different level of associated uncertainty within the same framework becomes key to make robust and unbiased decisions. The method according to the present invention allows assessing not only the geological parameters not known a priory with their associated uncertainty but also the geological and structural uncertainty of the reservoir with an accurate propagation of the inherent uncer- tainty. Normally, geological information comes from well logs, cores analysis, seismic data, etc. and production data can also be used to mitigate the uncertainty however in absence of these the geological scenario becomes hard to predict. In most of cases, geological information comprises scalar values lacking of spatial distributions as a three-dimensional structure. The method according to present invention generates scenarios of hydrocarbon reservoirs having a three-dimensional structure even if departing from scalar values. According to a first aspect of the invention, the method to generate scenarios of hydrocarbon reservoirs compiles information on the target hydrocarbon reservoir. Such information provides a concept model for the whole domain of the scenario. This information comprises scalar properties and may also comprise uncertainty data such as minimum and maximum values of some scalar data.

The method requires the proposal of an integer number n of facies. For each facie, the method provides automatically the partial rate f_it and a normal distribution of the porosity having a mean value φι and the cutoff values φ™¹⁷¹ and φ ^αχ for a predetermined threshold such that the area of the tail of the normal distribution is deemed to be neglectible beyond such threshold. Such variables are provided by solving an optimization problem that distributes the normal distribution of each facies in such a way the concept model is reproduced. The concept model is characterized by the number of facies and φ™ψ, φ™^^{χ '}, and the φ_τε^ values of porosity of the target hydrocarbon reservoir.

The method mainly uses the computed partial rate ^■ to populate a plurality of three- dimensional scenarios having a three-dimensional domain and the shape of the n facies by means of a multipoint statistical method, preferably and not limited to a Sequential Indicator Simulation.

Once the plurality of three-dimensional scenarios with the n facies are generated, each facie with its partial rate ^■ and the normal distribution of the porosity defined by the mean value φι and the cutoff values φ™¹⁷¹ and φ™¹¹*, further embodiments of the invention allow to populate the spatial distribution of the permeability, the saturation and the net-to-gross variables. DESCRIPTION OF THE DRAWINGS

These and other features and advantages of the invention will be seen more clearly from the following detailed description of a preferred embodiment provided only by way of illustrative and non-limiting example in reference to the attached drawings.

Figure 1 This figure shows an example of a domain having three facies located within the domain. Figure 2 This figure shows an example of a schematic representation of an analogous data base wherein selected records have been selected when comparing the information of the target according to a predefined criterion. Figure 3 This figure shows an example of the distribution of the porosity for the concept model and how the normal distribution for each facies is distributed reproducing the concept model by solving an optimization problem. Figure 4 This figure shows an example of the set of scenarios generated in each property calculation providing a first generation and a second generation of scenarios.

DETAILLED DESCRIPTION OF THE INVENTION

As will be appreciated by one skilled in the art, aspects of the present invention may be embodied as a system, method or computer program product. Accordingly, aspects of the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment (including firmware, resident software, micro-code, etc.) or an embodiment combining software and hardware aspects that may all generally be referred to herein as a "circuit," "module" or "system." Furthermore, aspects of the present invention may take the form of a computer program product embodied in one or more computer readable medium(s) having computer readable program code embodied thereon.

Any combination of one or more computer readable medium(s) may be utilized. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. A computer readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the fore- going. More specific examples (a non-exhaustive list) of the computer readable storage medium would include the following: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. A computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A computer readable signal medium may be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.

Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.

Computer program code for carrying out operations for aspects of the present invention may be written in any combination of one or more programming languages, in- eluding an object oriented programming language such as Java, Smalltalk, C++ or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider). Aspects of the present invention are described below with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer readable medium that can direct a computer, other programmable data processing apparatus, or other devices to function in a particular manner, such that the instructions stored in the computer readable medium produce an article of manufacture including instructions which implement the function/act specified in the flowchart and/or block diagram block or blocks.

The computer program instructions may also be loaded onto a computer, other pro- grammable data processing apparatus, or other devices to cause a series of operational steps to be performed on the computer, other programmable apparatus or other devices to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide processes for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.

Turning now to the drawings and more particularly, Figure 1 shows an example of an a prospective field for evaluation according to a preferred embodiment of the present invention, which begins by compiling information on the target hydrocarbon. This in- formation may be interpreted by a skilled person who may provide the number n of facies from this collected information. The number of facies is preferably between 1 and 8. Figure 1 shows a field having a domain (D) under certain portion of its surface, the area identified with letter A. This figure shows a domain (D) having three facies F_lt F₂ F₃ distributed within the domain. In absence of any information from seismic or maps the reservoir domain (D) is considered as a hexahedral shape. Structural uncertainty is considered on both shape and volume. In the present case since the conceptual model is a regular shape the uncertainty is represented by the reservoir area and thickness and the width-over-length ratio.

An interpolative function has been developed to estimate the joint distribution of reservoir area and thickness from the analogues data base whereas the width-over-length ratio parameter is assumed to be uniformly varying between its boundaries. Numerical analysis has shown that the reservoir area and thickness could be approximated with a lognormal bivariate distribution with a high confidence interval.

It is important to underline that the joint probability of area and thickness its key in order to preserve the interval of confidence given from the analogues in terms of volumetric and reserves. If there is information about the total amount of the oil of the reservoir or such amount may be estimated then the domain must be chosen big enough to contain that amount of oil.

According to an embodiment of the invention the number of facies and other properties may be obtained comparing the target information with the analog data base.

Each facies F_it i = l. . n, is modeled by its partial rate f_it and a normal distribution of the porosity having a mean value φι and the cutoff values φ™¹⁷¹ and φ ^αχ for a predetermined threshold such that the area of the tail of the normal distribution is deemed to be neglectible beyond such threshold. Once the threshold is predetermined, the normal distribution is taken to be limited between φ™¹⁷¹ and φ ^αχ and the mean value satisfies φι = (φ ¹⁰- — 0-^nm)/2. The partial rate ^■ satisfies

= 1 being n the number of facies already fixed.

As is shown in figure 3, the whole domain is modeled by a φ™ψ, φ™^^{χ '}, and φ_τε^ val- ues which may be provided by a skilled person or the estimation may be obtained from the analogues data base, or more briefly "the analogues". Next step calculates f_it φι, φ™¹¹*, φγ¹¹⁷¹, i = l.. n as the solution of a minimization problem wherein the cost function is defined by the expression

wherein φ™/^ι1 and φ™αΗί ^iS defined by the condition g£ = Mm( φ?^ίη); i = l, ... , n

Φ _Μ ^χ = Μαχί φ™^αχ); ί = 1, ... , η

and φι are taken such that φι < φ_ι+1 ί = 1. . n— 1. In an embodiment, the norm ||- || has been implemented as the Euclidean norm.

The objective function to be minimized represents the error between the estimated property as a function of facies type and the value provided for instance by the analogues. Figure 3 shows a four facies case (n = 4) wherein the four normal distributions of porosity are distributed according to the minimum of the cost function; and, the partial rate f_it i = 1, ... ,4 is simultaneously determined.

This minimization problem may be solved using Mixed Integer Linear Programming (MILP). The solution jointly computes the relative proportion of each facies and the conditional distribution of the properties per facies. This estimates will be consistent with the value of min, max and average used as the reference: φ™ψ, φ™^^{χ '}, and <$>_ref .

In an embodiment of the invention the reference values φ™ψ, φ™^^{χ '}, and φ_τε^ are obtained from the analogues data base (2). When compiling information on the target (1) hydrocarbon reservoir, said information at least have a material property among the material properties comprised in the records of the analog data base. The reference values are determined as follows: 1. When providing an analogues data base (2) comprising records of known reservoirs, at least one record further comprises the porosity of the material, that is, the minimum value φ™™ of the porosity, the mean value φ_αη of the porosity, and the maximum value φ£ ^αχ of the porosity.

These three values allow assuming a property distribution shape such as a Normal distribution. The porosity will be used to compute the reference value as the mean value of a subset of records of the analogs data having such property.

When compiling information on the target hydrocarbon reservoir in step a), such information has at least a material property among the material properties comprised in the records of the analog data base.

Common properties in the target hydrocarbon reservoir and in the analogs data base allows comparing the target (1) information and closely related analogues records (2.1). Common properties may be different to the porosity.

Provide a similarity criterion for the comparison of reservoirs sharing at least a common material property.

The similarity criterion is the criteria to define which records of the analogues data base are close to the target (1). This criterion may even combine properties or correlation between two or more properties.

Select the set of m records (2.1) in the analogues data base (2) meeting the similarity criterion with the target reservoir, wherein at least one record of the selected record comprises the porosity of the material,

Determining φ™ψ, φ™^^{χ '} , and the φ_τβ values wherein said values are obtained for the set of records of the analog data base selected in step 4) as:

- ¾ⁿ = Mean( ¾) ; = l ... m

^" 0re/ ⁼ ^^{e a} (0anj) 7 = 1— ^m>^{' a n}d,

- = Μβαη(φ™ ) 7 = 1 ... m, for all records having the porosity information.

In an embodiment of the present invention, the plurality of three-dimensional scenarios having a three-dimensional domain is obtained generating the dimensions of said domain randomly. If the analogues have volumetric data such as the reservoir area and thickness, NTG (net-to-gross) and saturation information, then the maximum capacity to store oil within the domain is calculated. If the randomly generated domain (D) is not big enough to store the calculated volume of oil then said domain (D) is discarded.

For example, the OOIP (Original oil in place) may be estimated as:

V φ NTG S

B wherein V is the volume of the domain, φ is the mean value of the porosity, NTG is the mean value of the net-to-gross value, S is the saturation and, B is the expansion factor of the oil when it is extracted from the reservoir and stored at ambient pressure.

Once the facies has been generated and f_it φι, φ™¹¹*, φγ¹¹⁷¹, ί = l. . n variables are determined, the porosity is spatially distributed based on its statistical distribution φι,

_{mea ns 0}f geo-statistical algorithm providing a predetermined longitudinal correlation equation and a predetermined orientation. In this case, the spatial distribution of the porosity is defined in each scenario of the plurality of scenarios.

After populating the porosity according to its statistical distribution, the longitudinal correlation equation and the orientation, the permeability property (K) is determined for each scenario.

To carry out this further step, the records of the analogues data base (2) have to comprise the permeability property of the materials and the porosity. From the selected records the correlation between the permeability versus the porosity is calculated wherein most of the cases the law between both variables is expressed as the loga- rithmic value of the permeability versus the porosity according to the Carman-Kozeny law. The selected records are those selected when the reference values φ™ψ, φ™^^{χ '}, and <$>_ref are obtained from the analogues data base (2).

According to the Carman-Kozeny law the equation

φ³ is used to fit the scatter from the relation between permeability (K) and porosity (ø). The coefficient C can be estimated using a least square regression. To allow low scatter in the distribution further filters on the analogues results can be applied such as exploit possible correlation with additional variables as could be depth. The permeability variance has also been parametrized using a linear regression model as a function of porosity in order to correctly propagate the uncertainty in the three/dimensional realizations. To have a reliable simulation it is important to separate the effects of the mean (the trend) which is deterministic with porosity and the variance which is simulated by geo-statistical techniques. The correlation law provides the scalar distribution of the permeability as a function of - li the porosity. By means of a geo statistical algorithm, for each domain of the generated scenarios the permeability as a scalar field is defined over said three-dimensional domain. For each scenario generated when the porosity has been populated, the geo- statistical algorithm generates a plurality of scenarios taking into account the variabil- ity of the permeability property. Figure 4 shows the initial scalar data, the first group of scenarios generated when populating the porosity property taking into account uncertainty; and, according to the last step, the plurality of scenarios departing from each former scenario wherein this second generation of scenarios provides the variability of the permeability taking into account how the uncertainty for said permeability proper- ty is propagated.

Further, once the porosity and the permeability have been distributed in the second generation of scenarios, according to a further embodiment of the invention the water saturation is calculated. To compute initial water saturation and later the net-to-gross, the results for analogous have been further exploit.

In this case, the selected records, those selected when the reference values φ™ψ, φ™^^{χ '}, and 0_ref are obtained from the analogues data base (2), also comprises the water saturation, or data on water saturation allowing to calculate the S for instance by means of correlations taking into account several properties different from the water saturation.

For the m selected records, being m > n, using the mean porosity conditioned on fades computed in the optimization algorithm together with the mean per each ana- logues j we can express the following over determined system of m equations:

n

is solved preferably by least square

t φι = φ_αη_,ί ^iS defined for all j = l, ... , m selected records of the analogues data base (2) and the n unknowns, the facies proportion per analogue, _t i = 1, ... , n.

Analogously the system of equatio

and n

NTGanj =∑f_l NTG_i

i=l

are specify for water saturation and net-to-gross calculation. In this case, the facies proportion per analogue, f_t i = 1, ... , n is known from the former system and the mean of the water saturation for each facies and the mean of the net-to-gross also for each facies is obtained by solving the over determined system, for instance, by least square method.

Claims

1.- Method to generate scenarios of hydrocarbon reservoirs based on limited amount of information on a target hydrocarbon reservoir, comprising the steps:

a. compile information on the target hydrocarbon reservoir,

b. providing an integer number n of facies, preferably between 1 and 8, comprised in the domain of the scenario modeling the target hydrocarbon reservoir, wherein each facies F_it i = l. . n, is modeled by its partial rate f_it and a normal distribution of the porosity having a mean value φι and the cutoff values φ™¹⁷¹ and φ ^αχ for a predetermined threshold such that the area of the tail of the normal distribution is deemed to be neglectible beyond such threshold, subject to the following restrictions:

c. determining φ™ψ, φ™^^{χ '}, and the φ_τε^ values,

d. determining f_{i t} φι, φ™¹¹*, ψγ¹¹⁷¹, ί = 1. . η as the solution of a minimization problem wherein the cost function is defined by the expression

wherein φ™/^ι 1 and is defined by the condition g£ = Min( ΦΥ¹ ) i = 1, ... , n

Φ _Μ ^χ = Μαχί φ™^αχ); ί = 1, ... , η generating a plurality of three-dimensional scenarios having a three-dimensional domain and the shape of the n facies by means of a multipoint statistical method, preferably a Sequential Indicator Simulation, wherein multipoint statistical method is restricted to the number of facies n, the partial rates of facies f_{i t} i = 1. . n.

2.- Method according to claim 1 wherein it further comprises providing an analog data base comprising records of known reservoirs wherein each record comprises scalar data of material properties for determining values such as facies properties.

3.- Method according to claim 1 and 2 wherein the φ™ψ, φ™^^{χ '}, and φ_τε^ values are obtained as follows: when providing an analog data base comprising records of known reservoirs, at least one record further comprises the porosity of the material, that is, the minimum value φ™™ of the porosity, the mean value φ_αη of the porosity, and the maximum value φ^^χ of the porosity,

when compiling information on the target hydrocarbon reservoir in step a), such information has at least a material property among the material properties comprised in the records of the analog data base,

provide a similarity criterion for the comparison of reservoirs sharing at least a common material property,

select the set of m records in the analog data base meeting the similarity criterion with the target reservoir, wherein at least one record of the selected record comprises the porosity of the material,

when determining φ™ψ, φ™^^{χ '}, and the φ_τε^ values in step c), such values are obtained for the set of records of the analog data base selected in step 4) as:

- ¾ⁿ = Mean( ¾) ; = l ... m

^" 0re/ ⁼ ^^{e a} (0anj) 7 = 1— ^m>^{' a n}d,

- φ™_β = Μβαη(φ™ 7 = 1 ... m for all records having the porosity in formation.

4.- Method according to claim 1 and 2, and any of previous claims wherein:

- the records of the analog data base further comprises volumetric data such as the reservoir area and thickness, NTG (net-to-gross) and saturation information, - when generating a plurality of three-dimensional scenarios on step e) having a three-dimensional domain, the domain is generated randomly from a multidimensional correlation between area and thickness and has a volume verifying that is big enough to store the oil calculated according to the volumetric data.

5.- Method according to any of previous claims wherein the multipoint statistical method for generating the plurality of three-dimensional scenarios is restricted in each scenario to the number of facies n, the partial rates of facies f_it i = l. . n, and the property φ is spatially distributed based on its statistical distribution φι, φ™¹⁷¹, ^ ^by means of geo-statistical algorithm and also to a predetermined longitudinal correlation equation and a predetermined orientation.

6. - Method according to any previous claim and in particular claim 2, wherein a plurality of records of the analog data base further comprises the permeability (K) property of the materials and the porosity,

- determining a correlation between the permeability versus the porosity, preferably the Carman-Kozeny law expressed as the logarithmic value of the permeability versus the porosity, obtained from the set of analog data base records selected in step 4) in claim 3,

- the dispersion of the former law is obtained from the same set of analog data base records selected in step 4) in claim 3,

- the permeability as a scalar field defined over the three-dimensional domain for each three-dimensional scenario is generated based on the correlation between permeability and porosity and its variability by means of geo-statistical algorithm .

7. - Method according claim 1, 2 and 3; and to any previous claim wherein the m se- lected records on step 4 of claim 3 of the analog data base further comprises the water saturation values of the materials S, or, data on water saturation allowing to calculate the S, form such data by postprocessing and the average water saturation for each facies is calculated as follows:

- for the m selected records, being m > n, solving the over determined system of m equations, preferably by mean square method:

n

for all j = l, ... , m records of the analogues data base and the n unknowns

being φι the mean porosity for each facie F_it i = 1. . n calculated on step e); and, Φαη_,ί ί ⁼ 1_>—, τη the mean values of the m selected records obtained from the analogues data base,

- for each facies F_it i = 1. . n, the average value of the saturation is calculated by solving the over determined system of equations, preferably the by least square method:

8. - Method according claim 1, 2 and 3; and to any previous claim wherein the m se- lected records on step 4 of claim 3 of the analog data base further comprises the net- to-gross NTG values of the materials NTG,; or, data on net-to-gross NTG allowing to calculate the NTG, form such data by postprocessing and the average NTG for each facies is calculated as follows:

n

for all j = l, ... , m records of the analogues data base and the n unknowns

being φι the mean porosity for each facie F_it i = 1. . n calculated on step e); and, Φαη_,ί i ⁼ 1)—, τη the mean values of the m selected records obtained from the analogues data base,

- for each facies F_it i = l. . n, the average value of the net-to-gross NTG_t is calculated by solving the over determined system of equations, preferably the by least square method: n

^f_l NTG_i = NTG_anj

i=l

9. - A data processing system comprising means adapted to carry out a method according to any of claims 1 to 8.

10. - A computer program product adapted to perform a method according to any of claims 1 to 8.