Computer implemented method for characterizing a target fluid of a hydrocarbon reservoir under uncertainty

Abstract

The present invention is related to computer implemented method for characterizing a target fluid of a hydrocarbon reservoir under uncertainty providing a set of fluid properties from a limited amount of information.

Claims

1. A computer implemented method for assessing the viability of fluid production from a hydrocarbon reservoir using characteristics of a target fluid present within the hydrocarbon reservoir, wherein characterization of the target fluid, at least in part, is determined based on properties of the hydrocarbon reservoir in which the target fluid is located, or on properties of the hydrocarbon reservoir rock that stores the target fluid, and the viability of fluid production is simulated through facilities of the reservoir having a predetermined field architecture including at least one pipe connecting a first point of fluid extraction located within the reservoir and a final end point of fluid extraction located at surface facilities of the reservoir, the method comprising the following steps: a) establishing m predetermined fluid properties (p.sub.1, p.sub.2, . . . , p.sub.m) including at least a first fluid property (p.sub.1) directed to one of the phases of the target fluid present in the hydrocarbon reservoir, either liquid or gas, and a second fluid property (p.sub.2) directed to a gas/liquid ratio of the target fluid; b) assigning a probability density function (ƒ.sub.1) for the first fluid property (p.sub.1) and a probability density function (ƒ.sub.2) for the second fluid property (p.sub.2); c) providing m−2 correlations c.sub.3, . . . , c.sub.m dependent at least on p.sub.1 and/or p.sub.2 properties such that p.sub.i=c.sub.i(p.sub.1, p.sub.2); i=3, . . . , m; d) for each property p.sub.i, i=3, . . . , m generating, automatically by a computer system, r sample values c.sub.i({tilde over (p)}.sub.1, {tilde over (p)}.sub.2) wherein {tilde over (p)}.sub.1 is a random sample according to its probability density function ƒ.sub.1 and {tilde over (p)}.sub.2 is a random sample according to its probability density function ƒ.sub.2; generating, automatically by a computer system, a probability density function ƒ.sub.i for property p.sub.i from the r sample values c.sub.i({tilde over (p)}.sub.1, {tilde over (p)}.sub.2); e) based on the probability density functions (ƒ.sub.1, ƒ.sub.2, . . . , ƒ.sub.m), determining the characterizing properties p.sub.i, i=1, . . . , m of the target fluid present in the hydrocarbon reservoir; f) using the characterizing properties of the target fluid determined in step e), assessing the fluid production viability of the reservoir by simulating the evolution of the target fluid along the piping of the predetermined field architecture.

2. The method according to claim 1 wherein the first fluid property (p.sub.1) directed to one of the phases of the target fluid is the American Petroleum Institute (API) gravity and the second fluid property (p.sub.2) directed to the gas/liquid ratio of the target fluid is the gas-oil ratio (GOR) ratio.

3. The method according to claim 1, wherein correlation c.sub.i(p.sub.1, p.sub.2) also depends on at least one additional fluid property p.sub.j, being j<i.

4. The method according to claim 3, wherein the fluid properties p.sub.i, i=3, . . . , m are sorted such that dependency of correlations is c.sub.3=c.sub.3 (p.sub.1, p.sub.2), c.sub.4=c.sub.4(p.sub.1, p.sub.2, p.sub.3), . . . , c.sub.i=c.sub.i(p.sub.1, p.sub.2, . . . , p.sub.i-1) wherein not necessarily all properties p.sub.3, . . . , p.sub.i-1 appear as dependent parameters.

5. The method according to claim 1, wherein it further comprises providing an analog data base, said analog data base executable in a computer system, comprising records of fluids with fluid properties from known reservoirs, wherein at least one record of the analog data base further comprises the first fluid property (p.sub.1), the second fluid property (p.sub.2) or both; providing a similarity module for the comparison of fluid properties of the target fluid and fluid properties of records of the analog data base, the module providing a similarity value as a result of said comparison, wherein before determining a property p.sub.i by generating r sample values in step c), i being an integer value in the interval [3,m], the following steps are carried out: i. determining a threshold parameter (sim); ii. selecting the set of records of the analog data base, by means of the similarity module, having a similarity value higher than the threshold parameter (sim) when compared to a target fluid, iii. determining a modified correlation c′.sub.i as c′.sub.i=c.sub.i*C+Δc.sub.i, being C and Δc.sub.i constant values, said constant values being determined by imposing that the sum of the errors between the value provided by the new correlation c′.sub.i and the property value provided by each of the selected records of the analog data base is minimum; iv. in step c) using the modified correlation c′.sub.i when the r sample values c.sub.i({tilde over (p)}.sub.1, {tilde over (p)}.sub.2) are computed.

6. The method according to claim 5, wherein the properties compared between the target fluid and the records of the analog data base comprises geochemical properties, thermodynamic properties or both, said properties being quantitative, multiple binary or binary variables or any combination thereof.

7. The method according to claim 5, wherein the target fluid and at least one record of the analogous data base comprise at least one binary property wherein each type of the binary properties is set to 0 or 1 and comprises a weighing value and the similarity module comprises a binary property contribution value, said contribution value being weighed with the weighing value, said weighing value determined according to the following rules for each binary property: if the target fluid does not have the binary property and/or the record does not have the binary property, then the contribution value to the similarity value is set to 0 and the weighing value of the binary property is not taken into account; if the target fluid and the record have the binary property but the target fluid value is different to the record value the contribution value to the similarity value is set to 0 so it does not account the weighing value of the property; if the target fluid and the record have the binary property and the target fluid value is equal to the record value the weighing value is set to 1 so the full weighing value of the property is the contribution value to the similarity value.

8. The method according to claim 5, wherein the target fluid and at least one record of the analogous data base comprise at least one multiple binary property wherein each type of the multiple binary properties is set to 0 or 1 and comprises a weighing value and the similarity module comprises a multiple binary property contribution value, said contribution value being weighed with the weighing value, said weighing value determined according to the following rules for each multiple binary property: if the target fluid does not have the multiple binary property and/or the record does not have the multiple binary property, then the contribution value to the similarity value is set to 0 and the weighing value of the multiple binary property is not taken into account; if the target fluid and the record have not in common any of the multiple integer combination of values, the contribution value to the similarity value is set to 0 and the weighing value of the property is not taken into account; if the target fluid and the record have in common at least one of the multiple integer combination of values, the weighing value is set to 1 so the full weighing value of the property is the contribution value to the similarity value.

9. The method according to claim 5, wherein the target fluid and at least one record of the analogous data base comprise at least one quantitative property wherein each type of the quantitative properties has a continuous value and comprises a weighing value and, the similarity module comprises a quantitative property contribution value, said contribution value being weighed with the weighing value, said weighing value determined according to the following rules for each quantitative property: if the target fluid does not have the continuous value property and/or the record does not have the continuous value property, then the contribution value to the similarity value is set to 0 and the weighing value of the quantitative property is not taken into account; if the target fluid and the record have the same quantitative property the contribution value to the similarity value is weighed by the weighing value: $\left(1-\frac{.Math.{\mathrm{Vc}}_{\mathrm{ih}}-{\mathrm{Vc}}_{\mathrm{jh}}.Math.}{G_h}\right)$ wherein Vc.sub.ih is the h.sup.th fluid quantitative property Vc.sub.i of the target fluid among the Vc quantitative properties of the target fluid, being Vc the number of quantitative variables of the target fluid; Vc.sub.jh is the h.sup.th fluid quantitative property Vc.sub.j of the analog record among the Vc quantitative properties of the target fluid; and G.sub.h is the common range of the h.sup.th quantitative variables Vc.sub.i and Vc.sub.j.

10. The method according to claim 5, wherein the similarity value s.sub.ij between a target fluid and a record from the analog data base is computed as $s_{\mathrm{ij}}=\frac{\begin{array}{c}\underset{h=1}{\overset{\mathrm{Vc}}{.Math.}}\left(\left(1-\frac{.Math.{\mathrm{Vc}}_{\mathrm{ih}}-{\mathrm{Vc}}_{\mathrm{jh}}.Math.}{G_h}\right){\mathrm{Pc}}_h\right)+\\ \underset{k=1}{\overset{\mathrm{Vb}}{.Math.}}\left({\mathrm{Vb}}_k\star {\mathrm{Pb}}_k\right)+\underset{l=1}{\overset{\mathrm{Vmb}}{.Math.}}\left({\mathrm{Vmb}}_l\star {\mathrm{Pmb}}_l\right)\end{array}}{{\mathrm{Pc}}_t+{\mathrm{Pb}}_t+{\mathrm{Pmb}}_t}$ wherein Vc is the number of quantitative variables of the target fluid; Vc.sub.ih is the h.sup.th fluid quantitative property Vc.sub.i of the target fluid among the Vc quantitative properties of the target fluid; Vc.sub.jh is the h.sup.th fluid quantitative property Vc.sub.j of the analog record among the Vc quantitative properties of the target fluid; Pc.sub.h is a weighing value of the h.sup.th quantitative variable; G.sub.h is the common range of the h.sup.th quantitative variables Vc.sub.i and Vc.sub.j; Vb is the number of binary variables of the target fluid; Vb.sub.k is a coincidence value, 1 if there is coincidence or 0 otherwise, for the k.sup.th binary variable of the target fluid; Pb.sub.k is a weighing value for the k.sup.th binary variable Vmb is the number of multiple binary variables of the target fluid; Vmb.sub.l is a coincidence value, 1 if there is coincidence or 0 otherwise, for the l.sup.h multiple binary variable of the target fluid; Pmb.sub.l is a weighing value of the l.sup.th multiple binary variable; Pc.sub.t is the total weight of the quantitative variables; Pb.sub.t is the total weight of the binary variables; Pmb.sub.t is the total weight of the multiple binary variables.

11. The method according to claim 1, wherein the probability density function characterizing a determined property is expressed as: a piecewise probability density function, two or more percentile measures or, both.

12. The method according to claim 1 wherein ƒ.sub.1 and ƒ.sub.2 are uniform distributions.

13. The method according to claim 1 wherein at least one correlation c.sub.i comprises thermodynamic properties.

14. The method according to claim 1, wherein the p.sub.j property is defined by three percentiles, a p.sub.low-percentile, a p.sub.mean-percentile and a p.sub.high-percentile.

15. The method according to claim 1, further comprising: determining a temperature and pressure (T,P) domain for representing the behavior of the fluid along the piping of the facilities; establishing regions as risk regions in the (T,P) domain where the temperature and pressure variables (T,P) are such as the fluid shows a phase transition; determining a fluid defined by a set of properties p.sub.1, p.sub.2, . . . , p.sub.m according to claim 1, wherein each property j∈[1, m] can take values in a range [a.sub.j, b.sub.j] and it is defined from at least a discrete value of its probability density function, and wherein at least one of properties is the first fluid property p.sub.1; generating two samples of fluid, a first sample of low density ƒ.sup.l and a second sample of high density ƒ.sup.h; providing a module for the simulation in a computer system of the evolution of a fluid along the piping of the facilities for a predetermined flow rate, said piping comprising one or more pipes defined at least by the diameter, length, degree of thermal insulation and ambient temperature; simulating, by means of the module for the simulation, the evolution of the first sample of fluid ƒ.sup.l and the evolution of the second sample of fluid ƒ.sup.h along the piping of the facilities providing a first path and a second path respectively of the evolution in the (T,P) domain; if the first path and the second path do not fall within any of the risk regions then the field architecture is determined as viable for any fluid with densities between the low density and the high density chosen on step c), otherwise the field architecture is determined as not being viable.

16. A computer program product stored on a non-transitory computer-readable medium and comprising computer-implementable instructions which, when executed by a computer, cause the computer to carry out the method according to claim 1.

Description

DESCRIPTION OF THE DRAWINGS

(1) These and other features and advantages of the invention will be more clearly seen 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.

(2) FIG. 1 This figure shows an example of a preferred appraisal system for automatically selecting a target fluid with existing analogous members stored in a analog data base according to a preferred embodiment of the present invention.

(3) FIG. 2 This figure shows an example of table with a record storing the properties of a fluid, a target fluid or an analogous fluid. The property is disclosed by means of the name of the “parameter”, column 1, and the property value shown in column 2.

(4) FIG. 3 This figure shows a flowchart with the main steps of an embodiment of the present invention.

(5) FIG. 4 This figure shows three binary properties being defined as a Boolean value, (yes or no, true or false, for example) and, in column 2, the assigned value in this specific embodiment.

(6) FIG. 5 This figure shows two multiple binary properties where column 2 shows the assigned value in this specific embodiment.

(7) FIGS. 6A-6D These figures shows examples of quantitative values, Fluid ID, source Rock Period, Maturation Level or the Biodegradation (Moldowan Index). Column 2 shows the quantitative value for a set of the status of the continuous parameter.

(8) FIG. 7 This figure shows a graphic displaying the Molecular Weight obtained in three identified analogue records, the Molecular Weight obtained by using the correlation and the Molecular Weight obtained by the correlation once it has been modified by shifting and scaling.

DETAILED DESCRIPTION OF THE INVENTION

(9) 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.

(10) 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 foregoing. 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.

(11) 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.

(12) 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.

(13) 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, including 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).

(14) 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.

(15) 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.

(16) The computer program instructions may also be loaded onto a computer, other programmable 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.

(17) Turning now to the drawings and more particularly, FIG. 1 shows an example of an appraisal system (100) determining the fluid properties of a target fluid (e.g., in a newly discovered hydrocarbon reservoir), based on automatic pairing with cataloged and characterized existing population analogue fluid members according to a preferred embodiment of the present invention. The appraisal system (100) includes computers (102, 104, 106) (3 in this example) coupled, wired or wirelessly to, and communicate with each other over, network 108, e.g., a local area network (LAN), the Internet, an intranet or a combination thereof. Typically, the computers (102, 104, 106) include one or more processors, e.g., central processing unit (CPU) (110), memory (112) and local storage (114) with a catalog listing known or existing analogue fluid members and fluid or reservoir properties e.g., a reservoirs database.

(18) In this example, the fluid property generation will be based on the proposed correlations and the knowledge provided by the analogue data base stored in the local storage (114) that will be used for correcting the proposed correlations c.sub.i in order to meet as much as possible with the fluid properties stored within the stored analogue data base.

(19) Target fluid (116) properties and properties of analogue data base (114) may comprise binary properties, multiple binary properties or quantitative properties.

(20) FIG. 2 displays a table showing a single record, shortened with a few properties for this example, with two columns, a first column with the parameters of a fluid, the target fluid (116) or an analogue fluid (114), and a second column with the value of the parameter.

(21) Among the parameters of the record represented in FIG. 2 are the Fluid ID, the reservoir ID, using integer values for identifying the fluid and the reservoir respectively, the API Gravity value being a quantitative value, the GOR value being a quantitative value, the Kerogen Type being a multiple binary value, the Source Rock being a multiple binary value, the Source Rock Period being a multiple binary variable, the Maturation Level being a quantitative value and, a binary variable indicating if there is a fluid mixing Evidence and a binary variable indicating if washing evidence exists. These last two binary variables can be represented by a Boolean value and also by using an integer value, typically 1 and 0 according to a pre-established convention.

(22) Maturation Level shows a question mark as the shown record does not contain the Maturation Level value. Records of the analogous data base may contain one or more parameters identifying a specific property with no value. Because of this, the rules for comparing a target fluid (116) and an analogous fluid first require to identify if a determined common parameter exists in both fluids.

(23) Even if along the entire description the set of fluid properties are identified as such, fluid properties, said properties must be interpreted in a general form wherein some of these variables are not necessarily physical properties of a fluid but properties of the reservoir where the fluid is located or the rock storing the fluid as all of these parameters allow the characterization of the fluid.

(24) According to this definition, the fluid characterization by generating the fluid properties takes into account the geochemistry.

(25) In this example and turning now to FIG. 3, the preferred system (100) first determines a first fluid property (p.sub.1), the API gravity of the liquid phase of the target fluid (116), and a second fluid property (p.sub.2) directed to the gas/liquid ratio of the target fluid (116), the GOR property. In an embodiment, the API gravity of the liquid phase of the target fluid (116) and GOR property is obtained from the information provided by the Basin Modeling or any other predicting model.

(26) As it is shown in the first step (1000), system (100) also comprises the following steps: a) As it is shown in FIG. 3, in the first step (1000) the target fluid (1000) is identified with the available input data (1001), the Fluid ID, the reservoir ID, and other parameters with properties that may be obtained for instance by laboratory assays or from samples obtained when drilling a wellbore in the reservoir. b) In a second step (1002) and in this embodiment the system (100) determines a probability density function (ƒ.sub.1) for the API gravity (p.sub.1) and a probability density function (ƒ.sub.2) for the GOR property (p.sub.2) characterizing these two first properties of the target fluid. In this particular case, a uniform probability density function is used for ƒ.sub.1 and ƒ.sub.2. c) In a third step (1003), m−2 correlations c.sub.3, . . . , c.sub.m are provided, said correlations dependent at least on p.sub.1 and/or p.sub.2 properties, wherein additional dependencies may exist, that is, p.sub.i=c.sub.i(p.sub.1, p.sub.2 . . . , p.sub.i-1); i=3, . . . , m. FIG. 3 shows an iterator i value ranging from 3 to m wherein for each iterator value three sub-steps (1004, 1005) are executed: a first sub-step (1004) wherein the identified c.sub.i({tilde over (p)}.sub.1, {tilde over (p)}.sub.2, . . . , {tilde over (p)}.sub.i-i) correlations are provided, a second sub-step (1005) generating, automatically by the computer system (100), r sample values c.sub.i({tilde over (p)}.sub.1, {tilde over (p)}.sub.2, . . . , {tilde over (p)}.sub.i-1) wherein {tilde over (p)}.sub.1, {tilde over (p)}.sub.2, . . . , {tilde over (p)}.sub.i-1 are random samples according to their corresponding probability density functions ƒ.sub.1, ƒ.sub.2, . . . , ƒ.sub.i-1. That is, each property p.sub.j, among those properties {tilde over (p)}.sub.1, {tilde over (p)}.sub.2, . . . , {tilde over (p)}.sub.i-i appearing explicitly in the dependency of c.sub.i, is sampled using the corresponding probability density function ƒ.sub.j. Those samples are introduced in the correlation c.sub.i for calculating each sample of c.sub.i({tilde over (p)}.sub.1, {tilde over (p)}.sub.2, . . . , {tilde over (p)}.sub.i-1) and, a third sub-step (1006) generating, automatically by the computer system (100), a probability density function ƒ.sub.i for property p.sub.i from the r sample values c.sub.i({tilde over (p)}.sub.1, {tilde over (p)}.sub.2). In this particular case, the r sample values allow to compute frequencies of sub-ranges according to a discretization of the domain for the property values providing a numerical approximation of the probability density function. d) In a forth step (1007) the method makes available the characterizing properties p.sub.i, i=1, . . . , m of the fluid as the probability density functions (ƒ.sub.1, ƒ.sub.2, . . . , ƒ.sub.m) being determined in step b) (1002) and the third sub-step (1006).

(27) In this example correlations c.sub.i({tilde over (p)}.sub.1, {tilde over (p)}.sub.2, . . . , {tilde over (p)}.sub.i-1) are dependent on {tilde over (p)}.sub.1 and/or {tilde over (p)}.sub.2 and other dependencies {tilde over (p)}.sub.j, j=3, . . . , i−1 may exist.

(28) According to this preferred embodiment, correlations c.sub.i provided in the first sub-step (1004) are modified by a shift value and a scaling value in order to force the correlation being a good representation of the analog fluids.

(29) In a first step, given the target fluid (116), a set of analogue fluids are identified and selected by implementing a similarity module which selects those records of the stored analogue data base (114) complying a similarity criteria, that is, given a threshold value (sim).

(30) Once a sub-set of records of the analogue data base are selected, the modified correlation c′.sub.i being expressed as c′.sub.i=c.sub.i*C+Δc.sub.i, being C and Δc.sub.i constant values, said constant values are determined by imposing that the sum of the errors between the value provided by the new correlation c′.sub.i and the property value provided by each of the selected records of the analog data base is minimum.

(31) In this adjustment process the selecting criterion is very relevant as it allows determining only those records that provide valuable information of the target fluid (116). A certain correlation c.sub.i involves properties p.sub.j j=1 . . . i but selected records of the analogues data base, when computing the modified correlation c′.sub.i, comprises additional properties such as geological, chemical and physical properties. These additional properties provide a contribution on the adjusting process improving the final expression c′.sub.i=c.sub.i*C+Δc.sub.i because of the influence of these additional parameters when selecting the analogous records.

(32) In this selection process, geological, chemical and physical properties of the fluid and the reservoir are assessed by identifying binary properties, multiple binary properties and quantitative properties. Each of these three kind of properties provide a contribution value and, all contribution values are weighed for computing the similarity value being used for assessing the distance between the target fluid and each of the records of the fluid of the analogue data base.

(33) Once all distances identified as a similarity index has been computed the analogue fluids are ranked by said similarity index. The selected analogue fluids are those satisfying that the absolute value of the similarity index is higher than the threshold value (sim).

(34) In this particular example the similarity value s.sub.ij between a target fluid and a record from the analog data base, taking into account the contribution of the three kinds of properties is computed as

(35) $s_{\mathrm{ij}}=\frac{\begin{array}{c}\underset{h=1}{\overset{\mathrm{Vc}}{.Math.}}\left(\left(1-\frac{.Math.{\mathrm{Vc}}_{\mathrm{ih}}-{\mathrm{Vc}}_{\mathrm{jh}}.Math.}{G_h}\right){\mathrm{Pc}}_h\right)+\\ \underset{k=1}{\overset{\mathrm{Vb}}{.Math.}}\left({\mathrm{Vb}}_k\star {\mathrm{Pb}}_k\right)+\underset{l=1}{\overset{\mathrm{Vmb}}{.Math.}}\left({\mathrm{Vmb}}_l\star {\mathrm{Pmb}}_l\right)\end{array}}{{\mathrm{Pc}}_t+{\mathrm{Pb}}_t+{\mathrm{Pmb}}_t}$
wherein
Vc is the number of quantitative variables of the target fluid;
Vc.sub.ih is the h.sup.th fluid quantitative property Vc.sub.i of the target fluid among the
Vc quantitative properties of the target fluid;
Vc.sub.jh is the h.sup.th fluid quantitative property Vc.sub.j of the analog record among the Vc quantitative properties of the target fluid;
Pc.sub.h is a weighing value of the h.sup.th quantitative variable;
G.sub.h is the common range of the h.sup.th quantitative variables Vc.sub.i and Vc.sub.j;
Vb is the number of binary variables of the target fluid;
Vb.sub.k is a coincidence value, 1 if there is coincidence or 0 otherwise, for the k.sup.th binary variable of the target fluid;
Pb.sub.k is a weighing value for the k.sup.th binary variable
Vmb is the number of multiple binary variables of the target fluid;
Vmb.sub.l is a coincidence value, 1 if there is coincidence or 0 otherwise, for the l.sup.th multiple binary variable of the target fluid;
Pmb.sub.l is a weighing value of the l.sup.th multiple binary variable;
Pc.sub.t is the total weight of the quantitative variables;
Pb.sub.t is the total weight of the binary variables;
Pmb.sub.t is the total weight of the multiple binary variables.

(36) FIG. 4 shows a table with three binary variables, a property indicating if there is TSR evidence, a property indicating if there is washing evidence and a property indicating if there is fluid mixing evidence.

(37) Each type of binary variable has assigned an integer value, 0 or 1, identifying each Boolean value allowing computing the property when compared to other properties such as the quantitative properties.

(38) These three binary properties are compared with records having a parameter identifying said properties.

(39) The numerator of the expression used for the computation of the similarity value s.sub.ij shows the term Σ.sub.k=1.sup.Vb(Vb.sub.k*Pb.sub.k) as the binary contribution value. This contribution value is being weighed by Pb.sub.k, the weighing value for the k.sup.th binary variable.

(40) FIG. 5 shows a table with two multiple binary variables, a property indicating the Kerogen Type of the fluid and property indicating the Source Rock.

(41) Each type of multiple binary variables has assigned an integer value, 0 or 1, identifying each type allowing computing the property when compared to other properties such as the quantitative properties.

(42) The same numerator of the expression used for the computation of the similarity value s.sub.ij shows the term Σ.sub.l=1.sup.Vmb(Vmb.sub.l*Pmb.sub.l) as the multiple binary contribution value. This contribution value is being weighed by Pmb.sub.l, the weighing value of the l.sup.th multiple binary variable.

(43) And finally, the same numerator of the expression used for the computation of the similarity value s.sub.ij shows the term

(44) $\underset{h=1}{\stackrel{\mathrm{Vc}}{.Math.}}\left(\left(1-\frac{.Math.{\mathrm{Vc}}_{\mathrm{ih}}-{\mathrm{Vc}}_{\mathrm{jh}}.Math.}{G_h}\right){\mathrm{Pc}}_h\right)$
as the quantitative contribution value. This contribution value is being weighed by Pc.sub.h, the weighing value of the h.sup.th quantitative variable.

(45) FIGS. 6A-6D show four tables with four quantitative properties respectively, only as an example, the Fluid ID, the Source Rock Period and the Maturation Level. The right column shows the value of each parameter.

(46) The computed similarity value s.sub.ij is being divided by three total weights,
Pc.sub.t+Pb.sub.t+Pmb.sub.t
normalizing the similarity value s.sub.ij.

(47) In each type of variables the weighing value is being penalized if there is no datum in the analogue record or for the target fluid.

(48) Each property has an associated value (being 0 or 1 for binary properties and multiple binary properties) that then is multiplied by a weighing value.

(49) For quantitative values, the value of each quantitative variable is determined by the following factor:

(50) $\left(1-\frac{.Math.{\mathrm{Vc}}_{\mathrm{ih}}-{\mathrm{Vc}}_{\mathrm{jh}}.Math.}{G_h}\right)$
and being weighed by the weighing value Pc.sub.h resulting the term

(51) $\underset{h=1}{\stackrel{\mathrm{Vc}}{.Math.}}\left(\left(1-\frac{.Math.{\mathrm{Vc}}_{\mathrm{ih}}-{\mathrm{Vc}}_{\mathrm{jh}}.Math.}{G_h}\right){\mathrm{Pc}}_h\right)$
already indicated above.

(52) Factor (1−|Vc.sub.ih−Vc.sub.jh|/G.sub.h) is the difference of the analogue value and the target value when compared to the total range of the variable. This difference ranging [0,1] multiply the weighing value of each variable providing part of the value or the whole value.

(53) For instance, the variable biodegradation is quantitative and has a weight of 13%. The range G.sub.h is 9 because it ranges from 1 to 10. If the target fluid has a biodegradation value of 1 and the analogue fluid has a value of 3, the difference is the absolute value of Vc.sub.ih−Vc.sub.jh that will result as 2. This result divided by G.sub.h (9) results is 0.222. The difference 1−0.222 is 0.788, the factor multiplying the weight of the biodegradation variable (13%). Therefore, the final result of the contribution of the quantitative biodegradation variable for this specific analogue record is 10.11%.

(54) The denominator of the similarity value s.sub.ij has the weight of all variables being the maximum value 100%.

(55) FIG. 7 shows an embodiment wherein three different analogous records have been identified with the target fluid. These three records comprise the API gravity property (API) and the Molecular Weight property (MW), the i.sup.th property being generated according to an embodiment of the invention. The proposed correlation c.sub.i for the Molecular Weight property is MW=630−(10*API).

(56) The modified correlation may be expressed in the form
MW=C[630−(10*API)]+Δc.sub.i
being C and Δc.sub.i constants to be determined by minimizing the sum of relative errors between the result of the modified correlation and the values of the properties of the analogue fluid.

(57) In this embodiment, C results 0.83 and Δc.sub.i results 0 (no need of shifting); that is, the modified correlation is MW=0.83*[630−(10*API)]

(58) FIG. 7 shows three different results, the Molecular Weight calculated (small diamonds) by using the original correlation, the Molecular Weight of the selected analogue fluids (circles) and the Molecular Weight calculated (squares) by using the modified correlation.

(59) The modified correlation is strongly dependent on the selection of analogue records, being this selection being carried out by using the disclosed method depending on the type of properties and the definition of the similarity value s.sub.ij.

(60) The modified correlation for MW is then used in step c) of the invention when generating, automatically by a computer system, a probability density function ƒ.sub.i for property p.sub.i from the r sample values calculated from MW (API).

(61) A practical application of the characterizing method of a target fluid is the generation of fluid properties in a hydrocarbon reservoir for assessing the viability of fluid production from said reservoir throughout predetermined field architecture under uncertainty of the fluid properties.

(62) When the facilities to be deployed in a field are designed, parameters as the dimensions of the pipes to be used, the path of said pipes from the reservoir to the surface, the thickness of the thermal insulation must be defined. Such parameters are determined according to the fluid properties as when the fluid flows from the reservoir to the surface, temperature and pressure decreases bringing on the appearance of solid phases that may prevent the normal flow of the fluid.

(63) In an embodiment of the invention, a method for assessing the viability of fluid production in a hydrocarbon reservoir and for certain facilities comprising at least one pipe connecting a first point of extraction located within the reservoir and a final end point located at surface facilities is proposed. The method comprises the following steps: determining a temperature and pressure (T,P) domain for representing the behavior of the fluid along the piping of the facilities. The path of the piping to be assessed is already determined. Real measurements and numerical production models are examples of means providing the temperature and pressure along a predetermined domain, in particular along the piping production system. establishing regions as risk regions in the (T,P) domain where the temperature and pressure variables (T,P) are such as the fluid shows a phase transition. Examples of regions having a risk in the (T,P) domain are those regions where the fluid shows a change of phase for instance solid appearance. Additional risk regions may be established, for instance regions where the viscosity is too high to ensure the flow. determining a fluid defined by a set of properties p.sub.1, p.sub.2 . . . , p.sub.m under uncertainty according to any of the embodiments already disclosed, wherein each property j∈[1, m] can take values in a range [a.sub.j, b.sub.j] and it is defined from at least a discrete value of its probability density function, and wherein at least one of property is the first fluid property p.sub.1; generating two samples of fluid, a first sample of low density ƒ.sup.l and a second sample of high density ƒ.sup.h; In this step the fluid is being characterized by defining probability density functions of each fluid property departing from two (p.sub.1, p.sub.2) properties, in particular these properties may be the API gravity and GOR which allows to define a low density ƒ.sup.l and a high density flow ƒ.sup.h. providing a module for the simulation in a computer system of the evolution of a fluid along the piping of the facilities for a predetermined flow rate, said piping comprising one or more pipes defined at least by the diameter, length, degree of thermal insulation and ambient temperature; simulating, by means of the module for the simulation, the evolution of the first sample of fluid ƒ.sup.l and the evolution of the second sample of fluid ƒ.sup.h along the piping of the facilities providing a first path and a second path respectively of the evolution in the (T,P) domain; The evolution of the low density ƒ.sup.l along the piping of the facilities and the evolution of the high density flow ƒ.sup.l along the same piping of the facilities provides two separated trajectories of the flow in the temperature and pressure (T,P) domain and a region between both trajectories. The evolution of the two fluids (ƒ.sup.l, ƒ.sup.h) is computed by the module for the simulation as the piping and other parameters of the facilities are already determined, the facilities being assessed. if the first path and the second path do not fall within any of the risk regions then the field architecture is determined as viable for any fluid with densities between the low density and the high density chosen on step c), otherwise the field architecture is determined as not being viable. The proposed method avoids the assessing of the entire region defined between the two trajectories of the two fluids (ƒ.sup.l, ƒ.sup.h) and only the first trajectory or path in the (T,P) domain taking into account the evolution of the first sample of fluid ƒ.sup.l and the second trajectory or path in the (T,P) domain taking into account the evolution of the second sample of fluid ƒ.sup.h are checked, verifying if any of these trajectories or paths fall within any of the already determined risk regions. If this is the case, the field architecture is determined as not being viable.

(64) When a field architecture defined by the facilities designed to be deployed in the field is determined as not being viable then said facilities must be re-designed for instance at the locations corresponding to the (T,P) domain falling within a risk region or, changing the entire flow conditions like the pump system.

(65) Once the field architecture is redefined the proposed method may be used again for assessing the viability of fluid production with the new facilities.