US20140350860A1 - Systems, methods, and computer-readable media for continuous capillary pressure estimation - Google Patents
Systems, methods, and computer-readable media for continuous capillary pressure estimation Download PDFInfo
- Publication number
- US20140350860A1 US20140350860A1 US13/902,333 US201313902333A US2014350860A1 US 20140350860 A1 US20140350860 A1 US 20140350860A1 US 201313902333 A US201313902333 A US 201313902333A US 2014350860 A1 US2014350860 A1 US 2014350860A1
- Authority
- US
- United States
- Prior art keywords
- log data
- modeled
- computer
- thomeer
- log
- 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.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 239000011148 porous material Substances 0.000 claims abstract description 83
- 230000035699 permeability Effects 0.000 claims abstract description 65
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 44
- 230000006870 function Effects 0.000 claims description 64
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 29
- 229910052753 mercury Inorganic materials 0.000 claims description 29
- 230000015654 memory Effects 0.000 claims description 19
- 238000003860 storage Methods 0.000 claims description 14
- 238000004088 simulation Methods 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 11
- 239000012530 fluid Substances 0.000 claims description 10
- 230000001419 dependent effect Effects 0.000 claims 2
- 230000008569 process Effects 0.000 description 22
- 238000009826 distribution Methods 0.000 description 18
- 238000010200 validation analysis Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- 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
- E21B49/008—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 by injection test; by analysing pressure variations in an injection or production test, e.g. for estimating the skin factor
-
- 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
- This invention relates generally to sedimentary basin exploration and reservoir hydrocarbon production and, more particularly, to estimating capillary pressure for basin exploration and reservoir development and production.
- Capillary pressure is a property used for basin exploration and reservoir development and production.
- capillary pressure may be used to compute the original hydrocarbon in place (OHIP) and to estimate a recovery factor.
- OHIP original hydrocarbon in place
- Reservoir modeling and simulation enables characterization of static reservoir properties and prediction of reservoir dynamic behavior. Exploration and production activities and strategies are strongly impacted by reservoir modeling workflows and techniques. Accordingly, accurate capillary pressure estimation will contribute in reducing uncertainties on reservoir model predictions and in optimizing the exploration activities and the reservoir development plan.
- Capillary pressure may be used to determine the saturation distribution and the total in situ volumes of fluids (e.g., oil, water, and gas).
- Capillary pressure data is obtained using special core analysis (SCAL).
- SCAL core analysis
- Capillary pressure measurement support size used for SCAL is significantly smaller than the modeling support size used in reservoir simulation that uses the capillary pressure data.
- Capillary pressure is typically measured by mercury injection capillary pressure (MICP).
- Other techniques for measuring capillary pressure may include core analysis using computer tomography (CT) and nuclear magnetic resonance (NMR).
- a method includes accessing well log data from a well log for a well, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data and determining Thomeer parameters from the permeability log data, the porosity log data, the water saturation log data, and the oil saturation log data.
- the Thomeer parameters include, for each pore system, a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure.
- Determining the Thomeer parameters includes determining a modeled permeability, determining a modeled porosity, and determining a modeled water saturation. Additionally, determining the Thomeer parameters includes evaluating an objective based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints.
- the objective function includes
- T is the Thomeer parameter
- Sw FAL is the value of the water saturation data
- So(T) is a modeled oil saturation for the current Thomeer parameter T.
- the one or more linear equality constraints include:
- Bv i is a fractional bulk volume occupied by mercury
- Pc is an applied capillary pressure
- ⁇ is the conversion factor from mercury-air to oil-water
- n is the number of pore systems in the reservoir
- ⁇ FAL is the porosity data.
- the one or more linear inequality constraints include:
- G i is the pore geometrical factor
- K(T) is the modeled permeability and K FAL is the permeability from log data.
- the method further includes determining the capillary pressure of the basin/reservoir using a Thomeer model having the Thomeer parameters.
- a non-transitory tangible computer-readable storage medium having executable computer code stored thereon for determining capillary pressure in a basin reservoir is provided.
- the computer code includes a set of instructions that causes one or more processors to perform the following operations: accessing well log data from a well log for a well, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data and determining Thomeer parameters from the permeability log data, the porosity log data, the water saturation log data, and the oil saturation log data.
- the Thomeer parameters include, for each pore system, a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure.
- Determining the Thomeer parameters includes determining a modeled permeability, determining a modeled porosity, and determining a modeled water saturation. Additionally, determining the Thomeer parameters includes evaluating an objective based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints.
- the objective function includes
- T is the Thomeer parameter
- Sw FAL is the value of the water saturation data
- So(T) is a modeled oil saturation for the current Thomeer parameter T.
- the one or more linear equality constraints include:
- Bv i is a fractional bulk volume occupied by mercury
- Pc is an applied capillary pressure
- ⁇ is the conversion factor from mercury-air to oil-water
- n is the number of pore systems in the reservoir
- ⁇ FAL is the porosity data.
- the one or more linear inequality constraints include:
- G i is the pore geometrical factor
- K(T) is the modeled permeability and K FAL is the permeability log data.
- the computer code includes a set of instructions that causes one or more processors to perform the following operations: determining the capillary pressure of the reservoir using a Thomeer model having the Thomeer parameters.
- a system includes well log data, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data, one or more processors, and a tangible non-transitory computer-readable memory having executable computer code stored thereon for determining capillary pressure in a basin and a reservoir.
- the computer code includes a set of instructions that causes one or more processors to perform the following operations: accessing well log data from a well log for a well, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data and determining Thomeer parameters from the permeability log data, the porosity log data, the water saturation log data, and the oil saturation log data.
- the Thomeer parameters include, for each pore system, a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure. Determining the Thomeer parameters includes determining a modeled permeability, determining a modeled porosity, and determining a modeled water saturation for the current Thomeer parameter T. Additionally, determining the Thomeer parameters includes evaluating an objective based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints.
- the objective function includes
- T is the Thomeer parameter
- Sw FAL is the value of the water saturation data
- So(T) is a modeled oil saturation.
- Bv i is a fractional bulk volume occupied by mercury
- Pc is an applied capillary pressure
- ⁇ is the conversion factor from mercury-air to oil-water
- n is the number of pore systems in the reservoir
- ⁇ FAL is the porosity data.
- the one or more linear inequality constraints include:
- G i is the pore geometrical factor
- K(T) is the modeled permeability and K FAL is the permeability from log data.
- the computer code further includes a set of instructions that causes one or more processors to perform the following operations: determining the capillary pressure of the basin/reservoir using a Thomeer model having the Thomeer parameters.
- a computer-implemented method for determining capillary pressure includes accessing well log data from a well log for a well, the well log data including permeability log data, porosity log data, water saturation log data, and oil saturation log data.
- the method further includes evaluating an objective function measuring the different between the permeability log data and a modeled permeability, the porosity log data and a modeled porosity, and the oil saturation log data and a modeled oil saturation, the modeled permeability, the modeled porosity, and the modeled oil saturation each a function of Thomeer parameters.
- the method includes determining the capillary pressures of the basin and a reservoir using a Thomeer model having the Thomeer parameters, the Thomeer parameters comprising a fractional bulk volume, for each pore system, a pore geometrical factor, and a minimum entry pressure.
- FIG. 1 is a graph of a Thomeer function representing an MICP B v -curve for a mono-modal pore system according to the prior art
- FIG. 2 is a block diagram of a process for estimating Thomeer parameters from standard well log data and determining capillary pressure in accordance with an embodiment of the present invention
- FIG. 3 is a block diagram of an optimization process for the objective function and constraints for estimating Thomeer parameters in accordance with an embodiment of the present invention
- FIG. 4 is a block diagram of a system for estimating Thomeer parameters from standard well log data and determining capillary pressure in accordance with an embodiment of the present invention
- FIG. 5 is a block diagram of a computer in accordance with an embodiment of the present invention.
- FIGS. 6-10 are graphs depicting the results of validation testing of the estimation of Thomeer parameters in accordance with an embodiment of the present invention.
- FIGS. 11-15 are plots of modeled data and actual data for the five synthetic wells used for the validation testing of the estimation of Thomeer parameters in accordance with an embodiment of the present invention.
- Well log data is obtained from a well log for a well and used to determine Thomeer parameters for each pore system, i.e., minimum entry pressure, pore geometrical factor, and fractional bulk volume occupied, used in a Thomeer model for determining capillary pressure.
- the well log data may include permeability log data, porosity log data, water saturation log data, and oil saturation log data.
- the Thomeer parameters are determined by evaluating an objective function that measures the mismatch between the well log data and modeled data having the Thomeer parameters as input.
- the objective function is iteratively evaluated using linear equality constraints, linear inequality constraints, and nonlinear equality constraints until convergence criteria are met.
- the evaluation may be performed using sequential quadratic programming.
- the determined capillary pressures for a basin and a reservoir may be provided for basin exploration, prospect evaluation, and reservoir modeling and reservoir simulation.
- mercury injection capillary pressure may be used with core samples to determine capillary pressures and pore size distributions.
- the mercury-air systems used with MICP may be converted to oil-water systems that typically exist in basin and a reservoir.
- the Thomeer model is based on an observed hyperbolic relationship between the amount of mercury entering a pore system in an MICP experiment and the applied mercury pressures.
- the Thomeer model provides an empirical formula for the occupied fractional bulk volume B v and the mercury pressures P c , as shown in Equation 1 as follows:
- Bv i (Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc
- Bv, ⁇ is the fractional bulk volume occupied by mercury at infinitely high capillary pressure
- G is the pore geometrical factor
- Pc is the applied capillary pressure
- Pd is the minimum entry pressure
- Equation 1 the shape of the function described above in Equation 1 is determined by the three Thomeer parameters described above: P d , G, and B v, ⁇ .
- the Thomeer parameterization for MICP experiments may be used to describe the internal architecture of a basin and a reservoir pore system.
- FIG. 1 depicts a Thomeer function representing an MICP Bv-curve for a mono-modal pore system.
- the pore geometrical factor G determines the curvature of the illustrated MICP Bv-curve, such that a large value of G results in a gradual onset and small value of G results in a sudden and sharp onset.
- a summation of up to three different Thomeer hyperbolas may be used to constitute a complete MICP B v -curve. Because each such hyperbola would upscale independently, multi-modality is left out of the upscale derivation.
- pore systems employed in the techniques described below may be assumed to be mono-modal. However, the techniques described below are not restricted to mono-modal pore systems and may be used with multi-modal pore systems.
- Equation 2 Equation 2
- Bv(Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc
- ⁇ is the porosity
- G is the pore geometrical factor
- Pc is the applied capillary pressure
- Pd is the minimum entry pressure
- Equation 2 may be expressed in terms of porosity and mercury saturations, as shown below in Equations 3 and 4 as follows:
- Bv(Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc
- ⁇ is the porosity
- S Hg (Pc) is the mercury saturation at capillary pressure Pc
- Pc is the applied capillary pressure
- S Hg (Pc) is the mercury saturation at capillary pressure Pc
- Pc is the applied capillary pressure
- G is the pore geometrical factor
- Pd is the minimum entry pressure
- the conversion from mercury saturation to actual reservoir fluid saturation may be performed using the interfacial tension values for oil-brine-rock and for Hg-air-rock (or Hg-vapor-rock).
- Typical estimations using Thomeer analysis involve the fitting of the Thomeer hyperbola to core plug MICP data by determining the Thomeer parameters from such data. This results in the full characterization of the entire pore space of each core plug.
- core plug MICP data is very sparse and Thomeer parameters estimation from MICP data will generate capillary pressure data with a high degree of uncertainty.
- a core plug may have a typical volume of about 10 cm 3 and reservoir element, as probed by wire-line logs such as resistivity, may represent a reservoir volume equivalent billions of core plugs, while a typical reservoir model grid cell is about 250 m ⁇ 250 m ⁇ 1 m. Consequently, reservoir modeling and simulation (e.g., for reserve estimation and production forecast) based on such capillary pressure data may also be highly uncertain.
- FIG. 2 depicts a process 200 for estimating the Thomeer parameters from standard well logs, such as those obtained from fluid analysis logs. Some or all steps of the process 200 may be implemented as executable computer code stored on a non-transitory tangible computer-readable storage medium and executed by one or more processors of a special-purpose machine, e.g., a computer programmed to execute the code.
- well log data is obtained (block 202 ).
- the obtained well log data 204 may include porosity data, water saturation data, oil saturation data, and permeability data.
- the Thomeer parameters for the Thomeer model are estimated from the well log data (block 206 ) by minimizing an objective function measuring a mismatch between observed data and predicted data from a theoretical model having the Thomeer parameters as an input.
- the estimated Thomeer parameters are inputted into the Thomeer formula (block 208 ). Based on these parameters, capillary pressure is estimated for a reservoir using the Thomeer equation (block 210 ). The estimated capillary pressures may then be used in further basin exploration, prospect evaluation, and reservoir modeling and simulations (block 212 ).
- embodiments of the present invention may estimate Thomeer parameters from “standard” well logs having, for example, porosity data, water saturation data, oil saturation data, and permeability data.
- the estimation of Thomeer parameters from well log data is based on an inverse problem theory. Accordingly, the Thomeer parameters are estimated by minimizing an objective function measuring a mismatch between observed data and predicted data from a theoretical model having the Thomeer parameters as an input.
- the Thomeer parameters may be abbreviated using Equation 5 as follows:
- T is the abbreviation for the Thomeer parameters
- n is the number of pore systems in the reservoir
- Bv i (Pc) 1 ⁇ i ⁇ n is the fractional bulk volume occupied by mercury at infinitely high capillary pressure
- Pc is the applied capillary pressure
- G i ) 1 ⁇ i ⁇ n is the pre geometrical factor
- Pd i ) 1 ⁇ i ⁇ n is the minimum entry pressure.
- the well log data may include a porosity log, water saturation log, oil saturation log, and permeability log.
- a porosity model using the multi-pore system porosity definition may be represented by Equation 6 as follows:
- ⁇ ⁇ ( T ) ⁇ ⁇ ⁇ i - 1 n ⁇ Bv i ( 6 )
- ⁇ is the porosity
- ⁇ is the conversion factor from mercury-air to oil-water
- n is the number of pore systems in the reservoir
- Bv i is the fractional bulk volume occupied by mercury.
- Equation 7 An oil saturation model using the Thomeer equations may be represented by Equations 7 and 8 as follows:
- So i (G i , Pd i ) is the modeled oil saturation using Thomeer equations for the pore system i described by the Thomeer parameters (Bv i (Pc), G i , Pd i ), Bv ⁇ is the fractional bulk volume occupied by mercury at infinitely high capillary pressure, Pc is the applied capillary pressure, G i is the pore geometrical factor and Pd i is the minimum entry pressure.
- So(T) is the modeled oil saturation using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T at capillary pressure Pc
- ⁇ is the porosity
- So i (G i , Pd i ) is the modeled oil saturation using Thomeer equations for the pore system i described by the Thomeer parameters (Bv i (Pc), G i , Pd i )
- Bv i is the fractional bulk volume occupied by mercury
- G i the pore geometrical factor
- Pd i is the minimum entry pressure.
- An absolute permeability model may be represented using the Buiting-Clerke equation, as shown below in Equation 9:
- K(T) is the modeled permeability using the using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T at capillary pressure Pc
- Bv i (Pc) is the fractional bulk volume occupied by mercury at infinitely high capillary pressure
- Pc is the applied capillary pressure
- G i is the pore geometrical factor
- Pd i is the minimum entry pressure.
- the Thomeer parameters are estimated from well log data based on an inverse problem theory.
- the inverse problem formulation includes an objection function minimized using linear and nonlinear constraints.
- the objection function is shown in Equation 10 as follows:
- F(T) is the objective function using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T
- w is the
- Sw FAL is the water saturation log value (e.g., from a fluid analysis log)
- So(T) is the modeled oil saturation using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T at capillary pressure Pc
- T is the abbreviation for the Thomeer parameters.
- Equation 11 The linear equality constraints are shown by Equation 11 as follows:
- Bv i (Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc
- Pc is the applied capillary pressure
- ⁇ is the conversion factor from mercury-air to oil-water
- n is the number of pore systems in the reservoir
- ⁇ FAL is the porosity log data (e.g., from a fluid analysis log).
- Equations 12-16 The linear inequality constraints are shown by Equations 12-16 as follows:
- Bv i is the fractional bulk volume occupied by mercury
- Bv i (Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc
- Pc is the applied capillary pressure
- n is the number of pore systems in the reservoir.
- G i is the pore geometrical factor and n is the number of pore systems in the basin/reservoir.
- Pd i is the minimum entry pressure and n is the number of pore systems in the reservoir.
- Bv i is the fractional bulk volume occupied by mercury
- Bv i (Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc
- Pc is the applied capillary pressure
- n is the number of pore systems in the reservoir.
- Pd i is the minimum entry pressure and n is the number of pore systems in the reservoir.
- Equation 17 The nonlinear equality constraints are shown by Equation 17 as follows:
- K(T) is the modeled permeability using the using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T and K FAL is the permeability log data (e.g., from a fluid analysis log).
- the objective function is evaluated using sequential quadratic programming (SQP).
- SQP sequential quadratic programming
- Equations 18-22 the objective function is evaluated using sequential quadratic programming.
- SQP is the standard general purpose technique to solve smooth nonlinear optimization problems under the following assumptions: the problem is not large, functions and gradients can be evaluated with sufficiently high precision, and the problem is smooth and well-scaled. Accordingly, in such embodiments a quadratic programming sub-problem may at solved at any iteration by linearizing the constraints and quadratically approximating the Lagrangian function shown in Equation 23 as follows:
- FIG. 3 depicts the optimization process 300 for the objective function and constraints described above. Some or all steps of the process 300 may be implemented as executable computer code stored on a non-transitory tangible computer-readable storage medium and executed by one or more processors of a special-purpose machine, e.g., a computer programmed to execute the code.
- the process 300 begins at initialization block 302 and is an iterative process that is stopped after convergence criteria are met. The initialization may include obtaining the well log data and receiving any other data or parameters used in the subsequent determinations. As shown in FIG. 3 , Equation 5, the abbreviation for the Thomeer parameters, is depicted in block 304 and represents the Thomeer parameters to be estimated via the process 300 .
- Equation 7 for estimating the oil saturation is depicted in block 306 .
- Equation 8 for estimating the oil saturation is estimated, as shown by Equation 8 depicted in block 308 .
- the modeled porosity as a function of the Thomeer parameters is estimated, as shown by Equation 6 depicted in block 310 .
- the modeled permeability as a function of the Thomeer parameters is estimated, as shown by Equation 9 depicted in block 312 .
- the evaluation of the equations described above is performed by minimizing the objective function described above in Equation 10 using the linear equality constraints described above in Equation 11, the linear inequality constraints described above in Equations 12-16, and the nonlinear inequality constraints described above in Equation 17 (block 314 ).
- the objection function may be optimized using SQP techniques.
- Next convergence criteria are evaluated to determine if the criteria are met (decision block 316 ). If the convergence criteria are not met (line 318 ), the next iteration of the process 300 is executed (block 320 ). In this manner, the equations described above are evaluated until the convergence criteria of the minimization of the objective function are met. If the convergence criteria are met (line 322 ), the process may stop (block 324 ) and the values for the Thomeer parameters may be used in subsequent calculations of capillary pressure.
- FIG. 4 depicts a system 400 for determining capillary pressure in accordance with an embodiment of the present invention.
- the system 400 may include well log data 402 and server 404 having a Thomeer parameter estimation process 406 and a capillary pressure estimation process 408 .
- FIG. 4 also illustrates a network 410 and a reservoir modeling and simulation system 412 .
- well log data 402 may be obtained from “standard” well logging techniques, such as a fluid analysis log, Facimage, and other techniques.
- no additional procedures, such as core plugs are required to obtain data for estimating the Thomeer parameters and the capillary pressure.
- the well log data includes porosity data 414 , saturation data 416 (e.g., water saturation data and oil saturation data), and permeability data 418 .
- the server 404 may execute a Thomeer parameter estimation process 406 to estimate Thomeer parameters 420 from the well log data 402 , according to the equations and processes described above.
- the Thomeer parameters 420 are estimated from the well log data 402 by minimizing an objective function measuring the mismatch between the observed well log data 402 and the predicted data from a theoretical model having the Thomeer parameters as the input.
- the server 404 may also execute the capillary pressure estimation process 408 to estimate capillary pressures 422 from the estimated Thomeer parameters 420 using the Thomeer equations described above.
- the estimated capillary pressures 422 may be provided to the reservoir modeling and simulation system 412 , such as via the network 410 .
- the reservoir modeling and simulation system may produce a reservoir model 424 and a reservoir simulation 426 using the capillary pressure as an input, in addition to other inputs typically provided for such models and simulations.
- the determination of capillary pressure according to the techniques described above is relatively quick as they are based on numerical modeling and not laboratory experiences.
- the cost of estimating capillary pressure may be reduced by using standard well log data for the estimation of capillary pressure as described above instead of laboratory experiments.
- the capillary pressure for the entire well is estimated and is not limited to only the relatively small number of score samples having SCAL. Accordingly, the estimated capillary pressure reflects the capillary pressures in the well and not just for the small plug volume used for SCAL. Consequently, the capillary pressure is more appropriate for reservoir simulation and modeling wherein the modeling cell size may be closer to the well neighborhood scale than to the core scale.
- FIG. 5 depicts a computer 500 , such as a server, in accordance with an embodiment of the present invention.
- the computer depicted in FIG. 5 may be used in conjunction with the present techniques.
- the computer 500 may communicate over a network 502 , described further below.
- the computer 500 may include various internal and external components that contribute to the function of the device and which may allow the computer 500 to function in accordance with the techniques discussed herein.
- various components of computer 500 may be provided as internal or integral components of the computer 500 or may be provided as external or connectable components.
- FIG. 5 depicts merely one example of a particular implementation and is intended to illustrate the types of components and functionalities that may be present in computer 500 .
- the computer 500 may be a server, a desktop computer, a laptop computer, a tablet computer, a smartphone, or other types of computers. As shown in FIG. 5 , the computer 500 may include one or more processors 504 and memory 506 . Additionally, the computer 500 may include, for example, an interface 508 , a display 510 , an input device 512 , input/output ports 514 and a network interface 516 .
- the display 510 may include a cathode ray tube (CRT) display, a liquid crystal display (LCD), an organic light emitting diode (OLED) display, or other types of displays.
- the display 510 may display a user interface (e.g., a graphical user interface) and may display various function and system indicators to provide feedback to a user, such as power status, call status, memory status, etc.
- the display 510 may include a touch-sensitive display (referred to as a “touch screen). In such embodiments, the touch screen may enable interaction with the computer via a user interface displayed on the display 510 .
- the display 510 may display a user interface for implementing the techniques described above, such as, for example, selecting well log data, initiating determination of Thomeer parameters, viewing the status of the processes described above, viewing the determined Thomeer parameters, viewing the determined capillary pressures, and so on.
- the one or more processors 504 may provide the processing capability required to execute the operating system, programs, user interface, and functions of the computer 500 .
- the one or more processors 500 may include microprocessors, such as “general-purpose” microprocessors, a combination of general and special purpose microprocessors, and Application-Specific Integrated Circuits (ASICs).
- the computer 500 may thus be a single processor system or a multiple processor system.
- the one or more processors 500 may include single-core processors and multicore processors and may include graphics processors, video processors, and/or related chip sets.
- the memory 506 may be accessible by the processor 502 and other components of the computer 500 .
- the memory 506 (which may include tangible non-transitory computer readable storage mediums) may include volatile memory and non-volatile memory accessible by the processor 502 and other components of the computer 500 .
- the memory 506 may store a variety of information and may be used for a variety of purposes. For example, the memory 506 may store the firmware for the computer 500 , an operating system for the computer 500 , and any other programs or executable code necessary for the computer 500 to function.
- the memory 506 may include volatile memory, such as random access memory (RAM) and may also include non-volatile memory, such as ROM, a solid state drive (SSD), a hard drive, any other suitable optical, magnetic, or solid-state storage medium, or a combination thereof.
- volatile memory such as random access memory (RAM)
- non-volatile memory such as ROM, a solid state drive (SSD), a hard drive, any other suitable optical, magnetic, or solid-state storage medium, or a combination thereof.
- the memory may store executable computer code that includes program instructions 518 executable by the one or more processors 502 to implement one or more embodiments of the present invention.
- program instructions 518 may include instructions 520 for Thomeer parameter estimation (e.g., the Thomeer parameter estimation process 406 ) and instructions 522 for capillary pressure estimation (e.g., the capillary pressure estimation process 408 ).
- the program instructions 518 may include a computer program (which in certain forms is known as a program, software, software application, script, or code).
- a computer program may be written in a programming language, including compiled or interpreted languages, or declarative or procedural languages.
- a computer program may include a unit suitable for use in a computing environment, including as a stand-alone program, a module, a component, a subroutine, etc., that may or may not correspond to a file in a file system.
- the program instructions 518 may be deployed to be executed on computers located locally at one site or distributed across multiple remote sites and interconnected by a communication network (e.g., network 502 ).
- the interface 508 may include multiple interfaces and may couple various components of the computer 500 to the processor 502 and memory 504 .
- the interface 508 , the processor 502 , memory 504 , and one or more other components of the computer 500 may be implemented on a single chip. In other embodiments, these components and/or their functionalities may be implemented on separate chips.
- the computer 500 also includes a user input device 512 that may be used to interact with and control the computer 500 .
- a user input device 512 may include any number of user input devices 512 , such as a keyboard, a mouse, a trackball, a digital stylus or pen, buttons, switches, or any other suitable input device.
- the input device 512 may be operable with a user interface displayed on the computer 500 to control functions of the computer 500 or of other devices connected to or used by the computer 500 .
- the input device 500 may allow a user to navigate a user interface, input data to the computer 500 , select data provided by the computer 500 , and direct the output of data from the computer 500 .
- the computer 500 may also include an input and output port 514 to enable connection of devices to the computer 500 .
- the input and output 514 may include an audio port, universal serial bus (USB) ports, AC and DC power connectors, serial data ports, and so on. Further, the computer 500 may use the input and output ports to connect to and send or receive data with other devices, such as other computers, printers, and so on.
- USB universal serial bus
- the computer 500 depicted in FIG. 5 also includes a network interface 516 , such as a network interface card (NIC), wireless (e.g., radio frequency) receivers, etc.
- the network interface 516 may include known circuitry for communicating with communication networks via electromagnetic signals transmitted over a wired or wireless connection. Such circuitry may include, for example, antennas, amplifiers, transceivers, receivers, processors, and so on.
- the network interface 516 may communicate with various communication networks (e.g., network 502 ), such as the Internet, an intranet, a cellular telephone network, a wireless local area network (LAN) a metropolitan area network (MAN), or other suitable communication networks.
- network 502 such as the Internet, an intranet, a cellular telephone network, a wireless local area network (LAN) a metropolitan area network (MAN), or other suitable communication networks.
- the network interface 516 may implement any suitable communications standard, protocol and/or technology, including wired Ethernet, wireless Ethernet (Wi-Fi) ((e.g., IEEE 802.11a, IEEE 802.11b, IEEE 802.11g and/or IEEE 802.11n), a 3G network (e.g., based upon the IMT-2000 standard), high-speed downlink packet access (HSDPA), wideband code division multiple access (W-CDMA), code division multiple access (CDMA), time division multiple access (TDMA), a 4G network (e.g., IMT Advanced, Long-Term Evolution Advanced (LTE Advanced), etc.), and any other suitable communications standard, protocol, or technology.
- FIGS. 6-10 depict the results of various validation tests to estimate reservoir capillary pressure in accordance with the techniques described herein.
- the validation testing is based on reference sets of Thomeer parameters that are used to populate five synthetic wells. For each synthetic well, porosity, permeability, and saturation logs are generated using the models described above in Equations 6-10. For validation, the actual Thomeer parameters are assumed to be unknown and the synthetic porosity, permeability, and saturation log data is used to estimate Thomeer parameters and, thus, the capillary pressure.
- FIG. 6 depicts graphs of the model well data vs. the estimated well data and illustrates the efficiency of the optimization of the objective function using the linear equality constraints, linear inequality constraints, and nonlinear equality constraints.
- FIG. 6 depicts a first graph 600 illustrating estimated saturation data vs. modeled saturation data, a second graph 602 illustrating estimated permeability data vs. modeled permeability data, and a third graph 604 illustrating estimated porosity data vs. modeled porosity data. The process took less than one minute to reach the convergence criteria and determine a solution that matches the saturation data and honors the linear constraints on the porosity and the non-linear constraints on the permeability
- FIG. 7 depicts graphs of the estimated Thomeer parameters (G, Bv, and Pd) vs. the actual Thomeer parameters for a macro system.
- FIG. 7 depicts a first graph 700 illustrating estimated pore geometry factor G vs. actual pore geometry factor G, a second graph 702 illustrating estimated fractional bulk volume By vs. actual fractional bulk volume By, and a third graph 704 illustrating estimated minimum entry pressure Pd vs. actual minimum entry pressure Pd.
- FIG. 7 also includes bar graphs that depict the distribution of the estimated Thomeer parameters as compared to the distribution of the actual Thomeer parameters obtained from a model database. As shown in FIG.
- the figure includes a first bar graph 706 that depicts the distribution of estimated and actual pore geometry factor G, a second bar graph 708 that depicts the distribution of estimated and actual fractional bulk volume Bv, and a third bar graph 710 that depicts the distribution of estimated and actual minimum entry pressure Pd.
- the distribution of the estimated Thomeer parameters is consistent with the distribution of the actual Thomeer parameters. Additionally, as shown from these figures, the range for the Thomeer parameters is also consistent.
- FIG. 8 depicts a first graph 800 illustrating estimated pore geometry factor G vs. actual pore geometry factor G, a second graph 802 illustrating estimated fractional bulk volume Bv vs. actual fractional bulk volume Bv, and a third graph 804 illustrating estimated minimum entry pressure Pd vs. actual minimum entry pressure Pd.
- G estimated Thomeer parameters
- Pd estimated Thomeer parameters
- the figure also includes a first bar graph 806 that depicts the distribution of estimated and actual pore geometry factor G, a second bar graph 808 that depicts the distribution of estimated and actual fractional bulk volume Bv, and a third bar graph 810 that depicts the distribution of estimated and actual minimum entry pressure Pd.
- the figures illustrate a higher difficulty of estimation due to the insensitivity of the permeability to the pore system 2 characteristics. However, these figures show that a consistent trend for fractional bulk volume Bv is obtained as the porosity is sensitive to the Bv values for the three pore systems.
- the minimum entry pressure Pd and the pore geometry factor G the low sensitivity of the saturation to these two parameters increases the difficulty of their estimation and the optimization program may generate a solution consistent with the prior term included in the objective function. As shown in the bar graphs 806 , 808 , and 810 , the prior term will control the consistency of these parameters' distributions.
- FIG. 9 depicts a first graph 900 illustrating estimated pore geometry factor G vs. actual pore geometry factor G, a second graph 902 illustrating estimated fractional bulk volume Bv vs. actual fractional bulk volume Bv, and a third graph 904 illustrating estimated minimum entry pressure Pd vs. actual minimum entry pressure Pd.
- the figures illustrate a higher difficulty of estimation due to the insensitivity of the permeability to the pore system 23 characteristics.
- FIG. 9 also depicts a first bar graph 806 that depicts the distribution of estimated and actual pore geometry factor G, a second bar graph 808 that depicts the distribution of estimated and actual fractional bulk volume Bv, and a third bar graph 810 that depicts the distribution of estimated and actual minimum entry pressure Pd that depict the control the prior term of the objective function over the parameters' distributions.
- FIG. 10 depicts the Thomeer parameters' distributions from for the model parameters and database parameters for the three pore systems discussed above and in accordance with an embodiment of the present invention.
- FIG. 10 illustrates the consistency between the modeled and actual Thomeer parameters and thus validates the utility of the capillary pressures derived from the Thomeer equation using Thomeer estimated according the techniques described herein.
- FIGS. 11-15 are plots of modeled data and actual data for the five synthetic wells used for the validation testing described above and in accordance with an embodiment of the present invention.
- Each plot includes: the match of the modeled porosity vs. the actual porosity data, the match of the modeled saturation vs. the actual saturation data; the match of the modeled permeability vs. the actual permeability data; the geometry factor G for the three pore systems; the bulk volume Bv for the three pore systems, and the minimum entry pressure Pd for the three pore systems.
- the plots illustrate that the techniques described estimate the Thomeer parameters for a multi-pore system that is consistent with the actual Thomeer parameters and well log data used for validation testing.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Debugging And Monitoring (AREA)
Abstract
Provided are methods, systems, and computer-readable media for determining capillary pressure in a basin/reservoir. Well log data is obtained that includes permeability log data, porosity log data, water saturation log data, and oil saturation log data. Thomeer parameters for a multi-pore system of a Thomeer model are determined by evaluating an objective function that measures the mismatch between the well log data and modeled data having the Thomeer parameters as input. The objective function is iteratively evaluated using linear equality constraints, linear inequality constraints, and nonlinear equality constraints until convergence criteria are met.
Description
- 1. Field of the Invention
- This invention relates generally to sedimentary basin exploration and reservoir hydrocarbon production and, more particularly, to estimating capillary pressure for basin exploration and reservoir development and production.
- 2. Description of the Related Art
- Various exploration and production systems are employed to find and extract oil, natural gas, and other resources from natural basins and reservoirs in the earth. Capillary pressure is a property used for basin exploration and reservoir development and production. For example, capillary pressure may be used to compute the original hydrocarbon in place (OHIP) and to estimate a recovery factor. Reservoir modeling and simulation enables characterization of static reservoir properties and prediction of reservoir dynamic behavior. Exploration and production activities and strategies are strongly impacted by reservoir modeling workflows and techniques. Accordingly, accurate capillary pressure estimation will contribute in reducing uncertainties on reservoir model predictions and in optimizing the exploration activities and the reservoir development plan. Capillary pressure may be used to determine the saturation distribution and the total in situ volumes of fluids (e.g., oil, water, and gas).
- Capillary pressure data is obtained using special core analysis (SCAL). However, due to the expense and time associated with SCAL procedures, only a relatively few measurements of capillary pressure are typically performed. Additionally, the capillary pressure measurement support size used for SCAL is significantly smaller than the modeling support size used in reservoir simulation that uses the capillary pressure data. Capillary pressure is typically measured by mercury injection capillary pressure (MICP). Other techniques for measuring capillary pressure may include core analysis using computer tomography (CT) and nuclear magnetic resonance (NMR).
- Various embodiments of methods, computer-readable media, and systems for determining capillary pressure in a basin and a reservoir are provided herein. In some embodiments, a method is provided that includes accessing well log data from a well log for a well, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data and determining Thomeer parameters from the permeability log data, the porosity log data, the water saturation log data, and the oil saturation log data. The Thomeer parameters include, for each pore system, a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure. Determining the Thomeer parameters includes determining a modeled permeability, determining a modeled porosity, and determining a modeled water saturation. Additionally, determining the Thomeer parameters includes evaluating an objective based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints. The objective function includes
-
- wherein T is the Thomeer parameter, SwFAL is the value of the water saturation data, and So(T) is a modeled oil saturation for the current Thomeer parameter T. The one or more linear equality constraints include:
-
- wherein Bvi is a fractional bulk volume occupied by mercury, Pc is an applied capillary pressure, α is the conversion factor from mercury-air to oil-water, n is the number of pore systems in the reservoir, and φFAL is the porosity data. The one or more linear inequality constraints include:
-
Bv i min ≦Bv i(Pc)≦Bv i max for 1≦i≦n -
G i min ≦G i ≦G i max for 1≦i≦n - wherein Gi is the pore geometrical factor,
-
Pd i min ≦Pd i ≦Pd i max for 1≦i≦n - wherein Pdi is a minimum entry pressure,
-
If Bv i(Pc)≠0 then Bv i+1(Pc)≦Bv i(Pc) for 1≦i≦n−1 -
Pd i ≦Pd i+1 for 1≦i≦n−1, - and the one or more nonlinear equality constraints include:
-
K(T)=K FAL (20) - wherein K(T) is the modeled permeability and KFAL is the permeability from log data. The method further includes determining the capillary pressure of the basin/reservoir using a Thomeer model having the Thomeer parameters.
- In other embodiments, a non-transitory tangible computer-readable storage medium having executable computer code stored thereon for determining capillary pressure in a basin reservoir is provided. The computer code includes a set of instructions that causes one or more processors to perform the following operations: accessing well log data from a well log for a well, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data and determining Thomeer parameters from the permeability log data, the porosity log data, the water saturation log data, and the oil saturation log data. The Thomeer parameters include, for each pore system, a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure. Determining the Thomeer parameters includes determining a modeled permeability, determining a modeled porosity, and determining a modeled water saturation. Additionally, determining the Thomeer parameters includes evaluating an objective based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints. The objective function includes
-
- wherein T is the Thomeer parameter, SwFAL is the value of the water saturation data, and So(T) is a modeled oil saturation for the current Thomeer parameter T. The one or more linear equality constraints include:
-
- wherein Bvi is a fractional bulk volume occupied by mercury, Pc is an applied capillary pressure, α is the conversion factor from mercury-air to oil-water, n is the number of pore systems in the reservoir, and φFAL is the porosity data. The one or more linear inequality constraints include:
-
Bv i ≦Bc i(Pc)≦Bv i max for 1≦i≦n -
G i min ≦G i ≦G i max for 1≦i≦n - wherein Gi is the pore geometrical factor,
-
Pd i min ≦Pd i ≦Pd i max for 1≦i≦n - wherein Pdi is a minimum entry pressure,
-
If Bv i(Pc)≠0 then Bv i+1(Pc)Bv i(Pc) for 1≦i≦n−1, -
Pd i ≦Pd i+1 for 1≦i≦n−1 - and the one or more nonlinear equality constraints include:
-
K(T)=K FAL (20) - wherein K(T) is the modeled permeability and KFAL is the permeability log data. The computer code includes a set of instructions that causes one or more processors to perform the following operations: determining the capillary pressure of the reservoir using a Thomeer model having the Thomeer parameters.
- Additionally, in some embodiments, a system is provided that includes well log data, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data, one or more processors, and a tangible non-transitory computer-readable memory having executable computer code stored thereon for determining capillary pressure in a basin and a reservoir. The computer code includes a set of instructions that causes one or more processors to perform the following operations: accessing well log data from a well log for a well, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data and determining Thomeer parameters from the permeability log data, the porosity log data, the water saturation log data, and the oil saturation log data. The Thomeer parameters include, for each pore system, a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure. Determining the Thomeer parameters includes determining a modeled permeability, determining a modeled porosity, and determining a modeled water saturation for the current Thomeer parameter T. Additionally, determining the Thomeer parameters includes evaluating an objective based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints. The objective function includes
-
- wherein T is the Thomeer parameter, SwFAL is the value of the water saturation data, and So(T) is a modeled oil saturation. The one or more linear equality constraints include:
-
- wherein Bvi is a fractional bulk volume occupied by mercury, Pc is an applied capillary pressure, α is the conversion factor from mercury-air to oil-water, n is the number of pore systems in the reservoir, and φFAL is the porosity data. The one or more linear inequality constraints include:
-
Bv i ≦Bc i(Pc)≦Bv i max for 1≦i≦n -
G i min ≦G i ≦G i max for 1≦i≦n - wherein Gi is the pore geometrical factor,
-
Pd i min ≦Pd i ≦Pd i max for 1≦≦n - wherein Pdi is a minimum entry pressure,
-
If Bv i(Pc)≠0 then Bv i+1(Pc)Bv i(Pc) for 1≦i≦n−1, -
Pd i ≦Pd i+1 for 1≦i≦n−1 - and the one or more nonlinear equality constraints include:
-
K(T)=K FAL - wherein K(T) is the modeled permeability and KFAL is the permeability from log data. The computer code further includes a set of instructions that causes one or more processors to perform the following operations: determining the capillary pressure of the basin/reservoir using a Thomeer model having the Thomeer parameters.
- Further, in some embodiments, a computer-implemented method for determining capillary pressure is provided. The method includes accessing well log data from a well log for a well, the well log data including permeability log data, porosity log data, water saturation log data, and oil saturation log data. The method further includes evaluating an objective function measuring the different between the permeability log data and a modeled permeability, the porosity log data and a modeled porosity, and the oil saturation log data and a modeled oil saturation, the modeled permeability, the modeled porosity, and the modeled oil saturation each a function of Thomeer parameters. Additionally, the method includes determining the capillary pressures of the basin and a reservoir using a Thomeer model having the Thomeer parameters, the Thomeer parameters comprising a fractional bulk volume, for each pore system, a pore geometrical factor, and a minimum entry pressure.
-
FIG. 1 is a graph of a Thomeer function representing an MICP Bv-curve for a mono-modal pore system according to the prior art; -
FIG. 2 is a block diagram of a process for estimating Thomeer parameters from standard well log data and determining capillary pressure in accordance with an embodiment of the present invention; -
FIG. 3 is a block diagram of an optimization process for the objective function and constraints for estimating Thomeer parameters in accordance with an embodiment of the present invention; -
FIG. 4 is a block diagram of a system for estimating Thomeer parameters from standard well log data and determining capillary pressure in accordance with an embodiment of the present invention; -
FIG. 5 is a block diagram of a computer in accordance with an embodiment of the present invention; -
FIGS. 6-10 are graphs depicting the results of validation testing of the estimation of Thomeer parameters in accordance with an embodiment of the present invention; -
FIGS. 11-15 are plots of modeled data and actual data for the five synthetic wells used for the validation testing of the estimation of Thomeer parameters in accordance with an embodiment of the present invention. - While the invention is susceptible to various modifications and alternative forms, specific embodiments thereof are shown by way of example in the drawings and will herein be described in detail. The drawings may not be to scale. It should be understood, however, that the drawings and detailed description thereto are not intended to limit the invention to the particular form disclosed, but to the contrary, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the present invention as defined by the appended claims.
- As discussed in more detail below, provided in some embodiments are systems, methods, and computer-readable media for determining capillary pressure in a basin and a reservoir. Well log data is obtained from a well log for a well and used to determine Thomeer parameters for each pore system, i.e., minimum entry pressure, pore geometrical factor, and fractional bulk volume occupied, used in a Thomeer model for determining capillary pressure. The well log data may include permeability log data, porosity log data, water saturation log data, and oil saturation log data. The Thomeer parameters are determined by evaluating an objective function that measures the mismatch between the well log data and modeled data having the Thomeer parameters as input. The objective function is iteratively evaluated using linear equality constraints, linear inequality constraints, and nonlinear equality constraints until convergence criteria are met. In some embodiments, the evaluation may be performed using sequential quadratic programming. The determined capillary pressures for a basin and a reservoir may be provided for basin exploration, prospect evaluation, and reservoir modeling and reservoir simulation.
- As will be appreciated, mercury injection capillary pressure (MICP) may be used with core samples to determine capillary pressures and pore size distributions. The mercury-air systems used with MICP may be converted to oil-water systems that typically exist in basin and a reservoir. The Thomeer model is based on an observed hyperbolic relationship between the amount of mercury entering a pore system in an MICP experiment and the applied mercury pressures. The Thomeer model provides an empirical formula for the occupied fractional bulk volume Bv and the mercury pressures Pc, as shown in
Equation 1 as follows: -
- Where Bvi(Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc, Bv,∞ is the fractional bulk volume occupied by mercury at infinitely high capillary pressure, G is the pore geometrical factor, Pc is the applied capillary pressure, and Pd is the minimum entry pressure.
- Thus, the shape of the function described above in
Equation 1 is determined by the three Thomeer parameters described above: Pd, G, and Bv,∞. The Thomeer parameterization for MICP experiments may be used to describe the internal architecture of a basin and a reservoir pore system.FIG. 1 depicts a Thomeer function representing an MICP Bv-curve for a mono-modal pore system. - As shown in the figure, the asymptotes for the Thomeer parameters Bv and Pd are indicated. The pore geometrical factor G determines the curvature of the illustrated MICP Bv-curve, such that a large value of G results in a gradual onset and small value of G results in a sudden and sharp onset. For multi-modal systems (e.g., carbonates), a summation of up to three different Thomeer hyperbolas may be used to constitute a complete MICP Bv-curve. Because each such hyperbola would upscale independently, multi-modality is left out of the upscale derivation. Thus, unless otherwise noted, pore systems employed in the techniques described below may be assumed to be mono-modal. However, the techniques described below are not restricted to mono-modal pore systems and may be used with multi-modal pore systems.
- Additionally, in some embodiments Bv,∞ may be assumed to equal φ. Thus, applying this assumption to
Equation 1, in such embodiments the Thomeer formula is shown inEquation 2 as follows: -
- Where Bv(Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc, φ is the porosity, G is the pore geometrical factor, Pc is the applied capillary pressure, and Pd is the minimum entry pressure.
- Accordingly,
Equation 2 may be expressed in terms of porosity and mercury saturations, as shown below in Equations 3 and 4 as follows: -
B v(P c)≈φS Hg(P c) (3) - Where Bv(Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc, φ is the porosity, SHg(Pc) is the mercury saturation at capillary pressure Pc, and Pc is the applied capillary pressure.
-
- Wherein SHg(Pc) is the mercury saturation at capillary pressure Pc, Pc is the applied capillary pressure, G is the pore geometrical factor, and Pd is the minimum entry pressure.
- As will be appreciated, the conversion from mercury saturation to actual reservoir fluid saturation (i.e., SHG to Soil) may be performed using the interfacial tension values for oil-brine-rock and for Hg-air-rock (or Hg-vapor-rock).
- Typically estimations using Thomeer analysis involve the fitting of the Thomeer hyperbola to core plug MICP data by determining the Thomeer parameters from such data. This results in the full characterization of the entire pore space of each core plug. However, such core plug MICP data is very sparse and Thomeer parameters estimation from MICP data will generate capillary pressure data with a high degree of uncertainty. For example, a core plug may have a typical volume of about 10 cm3 and reservoir element, as probed by wire-line logs such as resistivity, may represent a reservoir volume equivalent billions of core plugs, while a typical reservoir model grid cell is about 250 m×250 m×1 m. Consequently, reservoir modeling and simulation (e.g., for reserve estimation and production forecast) based on such capillary pressure data may also be highly uncertain.
-
FIG. 2 depicts aprocess 200 for estimating the Thomeer parameters from standard well logs, such as those obtained from fluid analysis logs. Some or all steps of theprocess 200 may be implemented as executable computer code stored on a non-transitory tangible computer-readable storage medium and executed by one or more processors of a special-purpose machine, e.g., a computer programmed to execute the code. Initially, well log data is obtained (block 202). The obtained well logdata 204 may include porosity data, water saturation data, oil saturation data, and permeability data. Next, as described in detail below, the Thomeer parameters for the Thomeer model are estimated from the well log data (block 206) by minimizing an objective function measuring a mismatch between observed data and predicted data from a theoretical model having the Thomeer parameters as an input. Next, the estimated Thomeer parameters are inputted into the Thomeer formula (block 208). Based on these parameters, capillary pressure is estimated for a reservoir using the Thomeer equation (block 210). The estimated capillary pressures may then be used in further basin exploration, prospect evaluation, and reservoir modeling and simulations (block 212). - As mentioned above, embodiments of the present invention may estimate Thomeer parameters from “standard” well logs having, for example, porosity data, water saturation data, oil saturation data, and permeability data. As described in detail below, the estimation of Thomeer parameters from well log data is based on an inverse problem theory. Accordingly, the Thomeer parameters are estimated by minimizing an objective function measuring a mismatch between observed data and predicted data from a theoretical model having the Thomeer parameters as an input. In the techniques described below, the Thomeer parameters may be abbreviated using Equation 5 as follows:
-
T=(Bv i(Pc),G i ,Pd i)1≦i≦n (5) - wherein T is the abbreviation for the Thomeer parameters, n is the number of pore systems in the reservoir, Bvi(Pc)1≦i≦n is the fractional bulk volume occupied by mercury at infinitely high capillary pressure, Pc is the applied capillary pressure, (Gi)1≦i≦n is the pre geometrical factor and (Pdi)1≦i≦n is the minimum entry pressure.
- As explained above, the well log data may include a porosity log, water saturation log, oil saturation log, and permeability log. A porosity model using the multi-pore system porosity definition may be represented by Equation 6 as follows:
-
- Where φ is the porosity, α is the conversion factor from mercury-air to oil-water, n is the number of pore systems in the reservoir, and Bvi is the fractional bulk volume occupied by mercury.
- An oil saturation model using the Thomeer equations may be represented by
Equations 7 and 8 as follows: -
- Where Soi (Gi, Pdi) is the modeled oil saturation using Thomeer equations for the pore system i described by the Thomeer parameters (Bvi(Pc), Gi, Pdi), Bv∞ is the fractional bulk volume occupied by mercury at infinitely high capillary pressure, Pc is the applied capillary pressure, Gi is the pore geometrical factor and Pdi is the minimum entry pressure.
-
- Where So(T) is the modeled oil saturation using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T at capillary pressure Pc, φ is the porosity, Soi (Gi, Pdi) is the modeled oil saturation using Thomeer equations for the pore system i described by the Thomeer parameters (Bvi(Pc), Gi, Pdi), Bvi is the fractional bulk volume occupied by mercury, Gi the pore geometrical factor and Pdi is the minimum entry pressure.
- An absolute permeability model may be represented using the Buiting-Clerke equation, as shown below in Equation 9:
-
- Wherein K(T) is the modeled permeability using the using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T at capillary pressure Pc, Bvi(Pc) is the fractional bulk volume occupied by mercury at infinitely high capillary pressure, Pc is the applied capillary pressure, Gi is the pore geometrical factor and Pdi is the minimum entry pressure.
- As mentioned above, the Thomeer parameters are estimated from well log data based on an inverse problem theory. The inverse problem formulation includes an objection function minimized using linear and nonlinear constraints. The objection function is shown in
Equation 10 as follows: -
- Wherein F(T) is the objective function using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T, w is the, SwFAL is the water saturation log value (e.g., from a fluid analysis log), So(T) is the modeled oil saturation using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T at capillary pressure Pc, and T is the abbreviation for the Thomeer parameters.
- The linear equality constraints are shown by Equation 11 as follows:
-
- Wherein Bvi(Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc, Pc is the applied capillary pressure, α is the conversion factor from mercury-air to oil-water, n is the number of pore systems in the reservoir, and φFAL is the porosity log data (e.g., from a fluid analysis log).
- The linear inequality constraints are shown by Equations 12-16 as follows:
-
Bv i min ≦Bv i(Pc)≦Bv i max for 1≦i≦n (12) - Wherein Bvi is the fractional bulk volume occupied by mercury, Bvi(Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc, Pc is the applied capillary pressure, and n is the number of pore systems in the reservoir.
-
G i min ≦G i ≦G i max for 1≦i≦n (13) - Wherein Gi is the pore geometrical factor and n is the number of pore systems in the basin/reservoir.
-
Pd i min ≦Pd i ≦Pd i max for 1≦i≦n (14) - Wherein Pdi is the minimum entry pressure and n is the number of pore systems in the reservoir.
-
if Bv i(Pc)≠0 then Bv i+1(Pc)≦Bv i(Pc) for 1≦i≦n−1 (15) - Wherein Bvi is the fractional bulk volume occupied by mercury, Bvi(Pc) is the fractional bulk volume occupied by mercury at capillary pressure Pc, Pc is the applied capillary pressure, and n is the number of pore systems in the reservoir.
-
Pd i ≦Pd i+1 for 1≦i≦n−1 (16) - Wherein Pdi is the minimum entry pressure and n is the number of pore systems in the reservoir.
- The nonlinear equality constraints are shown by Equation 17 as follows:
-
K(T)=K FAL (17) - Wherein K(T) is the modeled permeability using the using Thomeer equations for the multi-modal pore system described by the Thomeer parameters T and KFAL is the permeability log data (e.g., from a fluid analysis log).
- In some embodiments, the objective function is evaluated using sequential quadratic programming (SQP). By way of background, SQP may be used to solve differentiable nonlinear programming problems having forms shown by Equations 18-22 as follows:
-
min f(x) (18) -
x i ≦x≦x u (20) -
g j(x)=0,j=1, . . . ,m e (21) -
g j(x)≦,j=m e+1, . . . ,m (22) - Wherein x is an n-dimensional parameter vector and all problems functions f(x) and gj(x), j=1, . . . , m are assumed to be continuously differentiable. As will be appreciated, SQP is the standard general purpose technique to solve smooth nonlinear optimization problems under the following assumptions: the problem is not large, functions and gradients can be evaluated with sufficiently high precision, and the problem is smooth and well-scaled. Accordingly, in such embodiments a quadratic programming sub-problem may at solved at any iteration by linearizing the constraints and quadratically approximating the Lagrangian function shown in Equation 23 as follows:
-
- Wherein x is the primal variable and u=(u1, . . . , um)T is the Lagrange multiplier vector.
-
FIG. 3 depicts theoptimization process 300 for the objective function and constraints described above. Some or all steps of theprocess 300 may be implemented as executable computer code stored on a non-transitory tangible computer-readable storage medium and executed by one or more processors of a special-purpose machine, e.g., a computer programmed to execute the code. Theprocess 300 begins atinitialization block 302 and is an iterative process that is stopped after convergence criteria are met. The initialization may include obtaining the well log data and receiving any other data or parameters used in the subsequent determinations. As shown inFIG. 3 , Equation 5, the abbreviation for the Thomeer parameters, is depicted inblock 304 and represents the Thomeer parameters to be estimated via theprocess 300. Next, the various equations for the well log data as a function of the Thomeer parameters are determined. As shown inFIG. 3 , Equation 7 for estimating the oil saturation is depicted inblock 306. Subsequently, the modeled oil saturation as a function of the Thomeer parameters is estimated, as shown byEquation 8 depicted inblock 308. Additionally, as mentioned above, the modeled porosity as a function of the Thomeer parameters is estimated, as shown by Equation 6 depicted inblock 310. Similarly, the modeled permeability as a function of the Thomeer parameters is estimated, as shown byEquation 9 depicted inblock 312. - The evaluation of the equations described above is performed by minimizing the objective function described above in
Equation 10 using the linear equality constraints described above in Equation 11, the linear inequality constraints described above in Equations 12-16, and the nonlinear inequality constraints described above in Equation 17 (block 314). As mentioned above, in some embodiments the objection function may be optimized using SQP techniques. Next convergence criteria are evaluated to determine if the criteria are met (decision block 316). If the convergence criteria are not met (line 318), the next iteration of theprocess 300 is executed (block 320). In this manner, the equations described above are evaluated until the convergence criteria of the minimization of the objective function are met. If the convergence criteria are met (line 322), the process may stop (block 324) and the values for the Thomeer parameters may be used in subsequent calculations of capillary pressure. -
FIG. 4 depicts asystem 400 for determining capillary pressure in accordance with an embodiment of the present invention. Thesystem 400 may include well logdata 402 andserver 404 having a Thomeerparameter estimation process 406 and a capillarypressure estimation process 408. Additionally,FIG. 4 also illustrates anetwork 410 and a reservoir modeling andsimulation system 412. As will appreciated, well logdata 402 may be obtained from “standard” well logging techniques, such as a fluid analysis log, Facimage, and other techniques. Advantageously, no additional procedures, such as core plugs, are required to obtain data for estimating the Thomeer parameters and the capillary pressure. As noted above, in some embodiments the well log data includesporosity data 414, saturation data 416 (e.g., water saturation data and oil saturation data), andpermeability data 418. - As described above, the
server 404 may execute a Thomeerparameter estimation process 406 to estimateThomeer parameters 420 from the well logdata 402, according to the equations and processes described above. As described above, theThomeer parameters 420 are estimated from the well logdata 402 by minimizing an objective function measuring the mismatch between the observed well logdata 402 and the predicted data from a theoretical model having the Thomeer parameters as the input. Theserver 404 may also execute the capillarypressure estimation process 408 to estimatecapillary pressures 422 from the estimatedThomeer parameters 420 using the Thomeer equations described above. The estimatedcapillary pressures 422 may be provided to the reservoir modeling andsimulation system 412, such as via thenetwork 410. The reservoir modeling and simulation system may produce areservoir model 424 and areservoir simulation 426 using the capillary pressure as an input, in addition to other inputs typically provided for such models and simulations. - Advantageously, the determination of capillary pressure according to the techniques described above is relatively quick as they are based on numerical modeling and not laboratory experiences. Additionally, the cost of estimating capillary pressure may be reduced by using standard well log data for the estimation of capillary pressure as described above instead of laboratory experiments. Moreover, the capillary pressure for the entire well is estimated and is not limited to only the relatively small number of score samples having SCAL. Accordingly, the estimated capillary pressure reflects the capillary pressures in the well and not just for the small plug volume used for SCAL. Consequently, the capillary pressure is more appropriate for reservoir simulation and modeling wherein the modeling cell size may be closer to the well neighborhood scale than to the core scale.
-
FIG. 5 depicts acomputer 500, such as a server, in accordance with an embodiment of the present invention. The computer depicted inFIG. 5 , and other computers providing comparable capabilities, may be used in conjunction with the present techniques. Thecomputer 500 may communicate over anetwork 502, described further below. Thecomputer 500 may include various internal and external components that contribute to the function of the device and which may allow thecomputer 500 to function in accordance with the techniques discussed herein. As will be appreciated, various components ofcomputer 500 may be provided as internal or integral components of thecomputer 500 or may be provided as external or connectable components. It should further be noted thatFIG. 5 depicts merely one example of a particular implementation and is intended to illustrate the types of components and functionalities that may be present incomputer 500. - In various embodiments, the
computer 500 may be a server, a desktop computer, a laptop computer, a tablet computer, a smartphone, or other types of computers. As shown inFIG. 5 , thecomputer 500 may include one ormore processors 504 andmemory 506. Additionally, thecomputer 500 may include, for example, aninterface 508, adisplay 510, aninput device 512, input/output ports 514 and anetwork interface 516. - The
display 510 may include a cathode ray tube (CRT) display, a liquid crystal display (LCD), an organic light emitting diode (OLED) display, or other types of displays. Thedisplay 510 may display a user interface (e.g., a graphical user interface) and may display various function and system indicators to provide feedback to a user, such as power status, call status, memory status, etc. In some embodiments, thedisplay 510 may include a touch-sensitive display (referred to as a “touch screen). In such embodiments, the touch screen may enable interaction with the computer via a user interface displayed on thedisplay 510. In some embodiments, thedisplay 510 may display a user interface for implementing the techniques described above, such as, for example, selecting well log data, initiating determination of Thomeer parameters, viewing the status of the processes described above, viewing the determined Thomeer parameters, viewing the determined capillary pressures, and so on. - The one or
more processors 504 may provide the processing capability required to execute the operating system, programs, user interface, and functions of thecomputer 500. The one ormore processors 500 may include microprocessors, such as “general-purpose” microprocessors, a combination of general and special purpose microprocessors, and Application-Specific Integrated Circuits (ASICs). Thecomputer 500 may thus be a single processor system or a multiple processor system. The one ormore processors 500 may include single-core processors and multicore processors and may include graphics processors, video processors, and/or related chip sets. - The
memory 506 may be accessible by theprocessor 502 and other components of thecomputer 500. The memory 506 (which may include tangible non-transitory computer readable storage mediums) may include volatile memory and non-volatile memory accessible by theprocessor 502 and other components of thecomputer 500. Thememory 506 may store a variety of information and may be used for a variety of purposes. For example, thememory 506 may store the firmware for thecomputer 500, an operating system for thecomputer 500, and any other programs or executable code necessary for thecomputer 500 to function. Thememory 506 may include volatile memory, such as random access memory (RAM) and may also include non-volatile memory, such as ROM, a solid state drive (SSD), a hard drive, any other suitable optical, magnetic, or solid-state storage medium, or a combination thereof. - The memory may store executable computer code that includes
program instructions 518 executable by the one ormore processors 502 to implement one or more embodiments of the present invention. For example, theprocesses program instructions 518. Thus, in someembodiments program instructions 518 may include instructions 520 for Thomeer parameter estimation (e.g., the Thomeer parameter estimation process 406) andinstructions 522 for capillary pressure estimation (e.g., the capillary pressure estimation process 408). Theprogram instructions 518 may include a computer program (which in certain forms is known as a program, software, software application, script, or code). A computer program may be written in a programming language, including compiled or interpreted languages, or declarative or procedural languages. A computer program may include a unit suitable for use in a computing environment, including as a stand-alone program, a module, a component, a subroutine, etc., that may or may not correspond to a file in a file system. Theprogram instructions 518 may be deployed to be executed on computers located locally at one site or distributed across multiple remote sites and interconnected by a communication network (e.g., network 502). - The
interface 508 may include multiple interfaces and may couple various components of thecomputer 500 to theprocessor 502 andmemory 504. In some embodiments, theinterface 508, theprocessor 502,memory 504, and one or more other components of thecomputer 500 may be implemented on a single chip. In other embodiments, these components and/or their functionalities may be implemented on separate chips. - The
computer 500 also includes auser input device 512 that may be used to interact with and control thecomputer 500. In general, embodiments of thecomputer 500 may include any number ofuser input devices 512, such as a keyboard, a mouse, a trackball, a digital stylus or pen, buttons, switches, or any other suitable input device. Theinput device 512 may be operable with a user interface displayed on thecomputer 500 to control functions of thecomputer 500 or of other devices connected to or used by thecomputer 500. For example, theinput device 500 may allow a user to navigate a user interface, input data to thecomputer 500, select data provided by thecomputer 500, and direct the output of data from thecomputer 500. - The
computer 500 may also include an input and output port 514 to enable connection of devices to thecomputer 500. The input and output 514 may include an audio port, universal serial bus (USB) ports, AC and DC power connectors, serial data ports, and so on. Further, thecomputer 500 may use the input and output ports to connect to and send or receive data with other devices, such as other computers, printers, and so on. - The
computer 500 depicted inFIG. 5 also includes anetwork interface 516, such as a network interface card (NIC), wireless (e.g., radio frequency) receivers, etc. For example, thenetwork interface 516 may include known circuitry for communicating with communication networks via electromagnetic signals transmitted over a wired or wireless connection. Such circuitry may include, for example, antennas, amplifiers, transceivers, receivers, processors, and so on. Thenetwork interface 516 may communicate with various communication networks (e.g., network 502), such as the Internet, an intranet, a cellular telephone network, a wireless local area network (LAN) a metropolitan area network (MAN), or other suitable communication networks. Thenetwork interface 516 may implement any suitable communications standard, protocol and/or technology, including wired Ethernet, wireless Ethernet (Wi-Fi) ((e.g., IEEE 802.11a, IEEE 802.11b, IEEE 802.11g and/or IEEE 802.11n), a 3G network (e.g., based upon the IMT-2000 standard), high-speed downlink packet access (HSDPA), wideband code division multiple access (W-CDMA), code division multiple access (CDMA), time division multiple access (TDMA), a 4G network (e.g., IMT Advanced, Long-Term Evolution Advanced (LTE Advanced), etc.), and any other suitable communications standard, protocol, or technology. -
FIGS. 6-10 depict the results of various validation tests to estimate reservoir capillary pressure in accordance with the techniques described herein. The validation testing is based on reference sets of Thomeer parameters that are used to populate five synthetic wells. For each synthetic well, porosity, permeability, and saturation logs are generated using the models described above in Equations 6-10. For validation, the actual Thomeer parameters are assumed to be unknown and the synthetic porosity, permeability, and saturation log data is used to estimate Thomeer parameters and, thus, the capillary pressure. -
FIG. 6 depicts graphs of the model well data vs. the estimated well data and illustrates the efficiency of the optimization of the objective function using the linear equality constraints, linear inequality constraints, and nonlinear equality constraints.FIG. 6 depicts afirst graph 600 illustrating estimated saturation data vs. modeled saturation data, a second graph 602 illustrating estimated permeability data vs. modeled permeability data, and a third graph 604 illustrating estimated porosity data vs. modeled porosity data. The process took less than one minute to reach the convergence criteria and determine a solution that matches the saturation data and honors the linear constraints on the porosity and the non-linear constraints on the permeability -
FIG. 7 depicts graphs of the estimated Thomeer parameters (G, Bv, and Pd) vs. the actual Thomeer parameters for a macro system.FIG. 7 depicts afirst graph 700 illustrating estimated pore geometry factor G vs. actual pore geometry factor G, a second graph 702 illustrating estimated fractional bulk volume By vs. actual fractional bulk volume By, and a third graph 704 illustrating estimated minimum entry pressure Pd vs. actual minimum entry pressure Pd.FIG. 7 also includes bar graphs that depict the distribution of the estimated Thomeer parameters as compared to the distribution of the actual Thomeer parameters obtained from a model database. As shown inFIG. 7 , the figure includes a first bar graph 706 that depicts the distribution of estimated and actual pore geometry factor G, a second bar graph 708 that depicts the distribution of estimated and actual fractional bulk volume Bv, and a third bar graph 710 that depicts the distribution of estimated and actual minimum entry pressure Pd. As shown from these figures, the distribution of the estimated Thomeer parameters is consistent with the distribution of the actual Thomeer parameters. Additionally, as shown from these figures, the range for the Thomeer parameters is also consistent. -
FIG. 8 depicts graphs of the estimated Thomeer parameters (G, By, and Pd) vs. the actual Thomeer parameters for a second pore system having n=2.FIG. 8 depicts afirst graph 800 illustrating estimated pore geometry factor G vs. actual pore geometry factor G, a second graph 802 illustrating estimated fractional bulk volume Bv vs. actual fractional bulk volume Bv, and a third graph 804 illustrating estimated minimum entry pressure Pd vs. actual minimum entry pressure Pd. As shown inFIG. 8 , the figure also includes a first bar graph 806 that depicts the distribution of estimated and actual pore geometry factor G, a second bar graph 808 that depicts the distribution of estimated and actual fractional bulk volume Bv, and a third bar graph 810 that depicts the distribution of estimated and actual minimum entry pressure Pd. The figures illustrate a higher difficulty of estimation due to the insensitivity of the permeability to thepore system 2 characteristics. However, these figures show that a consistent trend for fractional bulk volume Bv is obtained as the porosity is sensitive to the Bv values for the three pore systems. For the minimum entry pressure Pd and the pore geometry factor G, the low sensitivity of the saturation to these two parameters increases the difficulty of their estimation and the optimization program may generate a solution consistent with the prior term included in the objective function. As shown in the bar graphs 806, 808, and 810, the prior term will control the consistency of these parameters' distributions. -
FIG. 9 depicts graphs of the estimated Thomeer parameters (G, Bv, and Pd) vs. the actual Thomeer parameters for a second pore system having n=3.FIG. 9 depicts afirst graph 900 illustrating estimated pore geometry factor G vs. actual pore geometry factor G, a second graph 902 illustrating estimated fractional bulk volume Bv vs. actual fractional bulk volume Bv, and a third graph 904 illustrating estimated minimum entry pressure Pd vs. actual minimum entry pressure Pd. The figures illustrate a higher difficulty of estimation due to the insensitivity of the permeability to the pore system 23 characteristics. Here again, however, the graphs show a consistent trend for fractional bulk volume Bv as the porosity is sensitive to the Bv values for the three pore systems. Similarly, as noted above, the low sensitivity of the saturation to the minimum entry pressure Pd and the pore geometry factor G increases the difficulty of their estimation, and the optimization program may generate a solution consistent with the prior term included in the objective function.FIG. 9 also depicts a first bar graph 806 that depicts the distribution of estimated and actual pore geometry factor G, a second bar graph 808 that depicts the distribution of estimated and actual fractional bulk volume Bv, and a third bar graph 810 that depicts the distribution of estimated and actual minimum entry pressure Pd that depict the control the prior term of the objective function over the parameters' distributions. -
FIG. 10 depicts the Thomeer parameters' distributions from for the model parameters and database parameters for the three pore systems discussed above and in accordance with an embodiment of the present invention.FIG. 10 illustrates the consistency between the modeled and actual Thomeer parameters and thus validates the utility of the capillary pressures derived from the Thomeer equation using Thomeer estimated according the techniques described herein. -
FIGS. 11-15 are plots of modeled data and actual data for the five synthetic wells used for the validation testing described above and in accordance with an embodiment of the present invention. Each plot includes: the match of the modeled porosity vs. the actual porosity data, the match of the modeled saturation vs. the actual saturation data; the match of the modeled permeability vs. the actual permeability data; the geometry factor G for the three pore systems; the bulk volume Bv for the three pore systems, and the minimum entry pressure Pd for the three pore systems. The plots illustrate that the techniques described estimate the Thomeer parameters for a multi-pore system that is consistent with the actual Thomeer parameters and well log data used for validation testing. - Further modifications and alternative embodiments of various aspects of the invention will be apparent to those skilled in the art in view of this description. Accordingly, this description is to be construed as illustrative only and is for the purpose of teaching those skilled in the art the general manner of carrying out the invention. It is to be understood that the forms of the invention shown and described herein are to be taken as examples of embodiments. Elements and materials may be substituted for those illustrated and described herein, parts and processes may be reversed or omitted, and certain features of the invention may be utilized independently, all as would be apparent to one skilled in the art after having the benefit of this description of the invention. Changes may be made in the elements described herein without departing from the spirit and scope of the invention as described in the following claims. Headings used herein are for organizational purposes only and are not meant to be used to limit the scope of the description.
- As used throughout this application, the word “may” is used in a permissive sense (i.e., meaning having the potential to), rather than the mandatory sense (i.e., meaning must). The words “include”, “including”, and “includes” mean including, but not limited to. As used throughout this application, the singular forms “a”, “an” and “the” include plural referents unless the content clearly indicates otherwise. Thus, for example, reference to “an element” includes a combination of two or more elements. Unless specifically stated otherwise, as apparent from the discussion, it is appreciated that throughout this specification discussions utilizing terms such as “processing”, “computing”, “calculating”, “determining” or the like refer to actions or processes of a specific apparatus, such as a special purpose computer or a similar special purpose electronic processing/computing device. In the context of this specification, a special purpose computer or a similar special purpose electronic processing/computing device is capable of manipulating or transforming signals, typically represented as physical electronic or magnetic quantities within memories, registers, or other information storage devices, transmission devices, or display devices of the special purpose computer or similar special purpose electronic processing/computing device.
Claims (33)
1. A computer-implemented method for determining capillary pressure in a reservoir, the method comprising:
accessing well log data from a well log for a well, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data;
determining Thomeer parameters from the permeability log data, the porosity log data, the water saturation log data, and the oil saturation log data, the Thomeer parameters comprising a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure, the determining comprising:
determining a modeled permeability;
determining a modeled porosity;
determining a modeled water saturation, and
evaluating an objective function based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints, the objective function comprising:
wherein T is the Thomeer parameters,
SwFAL is the value of the water saturation data,
So(T) is a modeled oil saturation;
the one or more linear equality constraints comprising:
wherein Bvi is a fractional bulk volume occupied by mercury,
Pc is an applied capillary pressure;
α is the conversion factor from mercury-air to oil-water,
n is the number of pore systems in the reservoir,
φFAL is the porosity data;
the one or more linear inequality constraints comprising:
Bv i min ≦Bv i(Pc)≦Bv i max for 1≦1≦n
G i min ≦G i ≦G i max for 1≦i≦n
Bv i min ≦Bv i(Pc)≦Bv i max for 1≦1≦n
G i min ≦G i ≦G i max for 1≦i≦n
wherein Gi is the pore geometrical factor,
Pd i min ≦Pd i ≦Pd i max for 1≦i≦n
Pd i min ≦Pd i ≦Pd i max for 1≦i≦n
wherein Pdi is a minimum entry pressure,
If Bv i(Pc)≠0 then Bv i+1(Pc)Bv i(Pc) for 1≦i≦n−1,
Pd i ≦Pd i+1 for 1≦i≦n−1, and
If Bv i(Pc)≠0 then Bv i+1(Pc)Bv i(Pc) for 1≦i≦n−1,
Pd i ≦Pd i+1 for 1≦i≦n−1, and
the one or more nonlinear equality constraints comprising:
K(T)=K FAL
K(T)=K FAL
wherein K(T) is the modeled permeability,
KFAL is the permeability log data; and
determining the capillary pressure of the reservoir using a Thomeer model having the determined Thomeer parameters.
2. The computer-implemented method of claim 1 , wherein the modeled permeability comprises:
3. The computer-implemented method of claim 1 , wherein the modeled porosity comprises:
4. The computer-implemented method of claim 1 , wherein the modeled oil saturation comprises:
5. The computer implemented method of claim 1 , wherein the Thomeer model comprises:
6. The computer-implemented method of claim 1 , wherein evaluating the objective function based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints comprises iteratively evaluating the objective function until convergence criteria are met.
7. The computer-implemented method of claim 1 , wherein evaluating the objective function based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints comprises iteratively evaluating the objective function evaluating the objective function using sequential quadratic programming (SQP).
8. The computer-implemented method of claim 1 , wherein the well log comprises a fluid analysis log.
9. The computer-implemented method of claim 1 , wherein the reservoir comprises an oil reservoir.
10. The computer-implemented method of claim 1 , comprising providing the capillary pressures to a reservoir modeling system, a reservoir simulation system, or a combination thereof.
11. A non-transitory tangible computer-readable storage medium having executable computer code stored thereon for determining capillary pressure in a reservoir, the computer code comprising a set of instructions that causes one or more processors to perform the following operations:
accessing well log data from a well log for a well, the well log data including permeability log data, porosity log data, water saturation log data, and oil saturation log data;
determining Thomeer parameters from the permeability log data, the porosity log data, the water saturation log data, and the oil saturation log data, the Thomeer parameters comprising a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure, the determining comprising:
determining a modeled permeability;
determining a modeled porosity;
determining a modeled water saturation, and
evaluating an objective function based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints, the objective function comprising:
wherein T is the Thomeer parameters,
SwFAL is the value of the water saturation data,
So(T) is a modeled oil saturation;
the one or more linear equality constraints comprising:
wherein Bvi is a fractional bulk volume occupied by mercury,
Pc is an applied capillary pressure;
α is the conversion factor from mercury-air to oil-water,
n is the number of pore systems in the reservoir,
φFAL is the porosity data;
the one or more linear inequality constraints comprising:
Bv i min ≦Bv i(Pc)≦Bv i max for 1≦i≦n
G i min ≦G i ≦G i max for 1≦i≦n
Bv i min ≦Bv i(Pc)≦Bv i max for 1≦i≦n
G i min ≦G i ≦G i max for 1≦i≦n
wherein Gi is the pore geometrical factor,
Pd i min ≦Pd i ≦Pd i max for 1≦i≦n
Pd i min ≦Pd i ≦Pd i max for 1≦i≦n
wherein Pdi is a minimum entry pressure,
If Bv i(Pc)≠0 then Bv i+1(Pc)≦Bv i(Pc) for 1≦i≦n−1,
Pd i ≦Pd i+1 for 1≦i≦n−1, and
If Bv i(Pc)≠0 then Bv i+1(Pc)≦Bv i(Pc) for 1≦i≦n−1,
Pd i ≦Pd i+1 for 1≦i≦n−1, and
the one or more nonlinear equality constraints comprising:
K(T)=K FAL
K(T)=K FAL
wherein K(T) is the modeled permeability,
KFAL is the permeability log data; and
determining the capillary pressures of the reservoir using a Thomeer model having the determined Thomeer parameters.
12. The non-transitory tangible computer-readable storage medium of claim 12 , wherein the modeled permeability comprises:
13. The non-transitory tangible computer-readable storage medium of claim 12 , wherein the modeled porosity comprises:
14. The non-transitory tangible computer-readable storage medium of claim 12 , wherein the modeled oil saturation comprises:
15. The computer implemented method of claim 1 , wherein the Thomeer model comprises:
16. The non-transitory tangible computer-readable storage medium of claim 12 , wherein evaluating the objective function based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints comprises iteratively evaluating the objective function until convergence criteria are met.
17. The non-transitory tangible computer-readable storage medium of claim 12 , wherein evaluating the objective function based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints comprises iteratively evaluating the objective function evaluating the objective function using sequential quadratic programming (SQP).
18. The non-transitory tangible computer-readable storage medium of claim 12 , wherein the well log comprises a fluid analysis log.
19. The non-transitory tangible computer-readable storage medium of claim 12 , wherein the reservoir comprises an oil reservoir.
20. A system for determining capillary pressure in a basin and reservoir, the system comprising:
well log data, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data;
one or more processors;
a tangible non-transitory computer-readable memory having executable computer code stored thereon for determining capillary pressure in a reservoir, the computer code comprising a set of instructions that causes the one or more processors to perform the following operations:
determining Thomeer parameters from the permeability log data, the porosity log data, the water saturation log data, and the oil saturation log data, the Thomeer parameters comprising a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure, the determining comprising:
evaluating an objective function based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints, the objective function comprising:
wherein T is the Thomeer parameters,
SwFAL is the value of the water saturation data,
So(T) is a modeled oil saturation;
the one or more linear equality constraints comprising:
wherein Bvi is a fractional bulk volume occupied by mercury,
Pc is an applied capillary pressure;
α is the conversion factor from mercury-air to oil-water,
n is the number of pore systems in the reservoir,
φFAL is the porosity data;
the one or more linear inequality constraints comprising:
Bv i min ≦Bv i(Pc)≦Bv i max for 1≦i≦n
G i min ≦G i ≦G i max for 1≦i≦n
Bv i min ≦Bv i(Pc)≦Bv i max for 1≦i≦n
G i min ≦G i ≦G i max for 1≦i≦n
wherein Gi is the pore geometrical factor,
Pd i min ≦Pd i ≦Pd i max for 1≦i≦n
Pd i min ≦Pd i ≦Pd i max for 1≦i≦n
wherein Pdi is a minimum entry pressure,
If Bv i(Pc)≠0 then Bv i+1(Pc)≦Bv i(Pc) for 1≦i≦n−1,
Pd i ≦Pd i+1 for 1≦i≦n−1, and
If Bv i(Pc)≠0 then Bv i+1(Pc)≦Bv i(Pc) for 1≦i≦n−1,
Pd i ≦Pd i+1 for 1≦i≦n−1, and
the one or more nonlinear equality constraints comprising:
K(T)=K FAL
K(T)=K FAL
wherein K(T) is the modeled permeability,
KFAL is the permeability log data; and
determining the capillary pressures of the reservoir using a Thomeer model having the determined Thomeer parameters.
21. The system of claim 20 , tangible non-transitory computer-readable memory storing a modeled permeability, a modeled porosity, and a modeled water saturation.
22. The system of claim 21 , wherein the modeled permeability comprises:
23. The system of claim 21 , wherein the modeled porosity comprises:
24. The system of claim 21 , wherein the modeled oil saturation comprises:
25. The system of claim 20 , the tangible non-transitory computer-readable memory storing a the Thomeer model, the Thomeer model comprising:
26. The system of claim 20 , wherein evaluating the objective function based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints comprises iteratively evaluating the objective function until convergence criteria are met.
27. The system of claim 20 , wherein evaluating the objective function based on one or more linear equality constraints, one or more linear inequality constraints, and one or more nonlinear equality constraints comprises iteratively evaluating the objective function evaluating the objective function using sequential quadratic programming (SQP).
28. The system of claim 20 , wherein the well log comprises a fluid analysis log.
29. The system of claim 20 , wherein the reservoir comprises an oil reservoir.
30. The system of claim 20 , comprising a network coupled to the one or more processor.
31. The system of claim 20 , comprising providing, over the network, the capillary pressures to a reservoir modeling system, a reservoir simulation system, or a combination thereof.
32. A computer-implemented method for determining capillary pressure in a reservoir, the method comprising:
accessing well log data from a well log for a well, the well log data comprising permeability log data, porosity log data, water saturation log data, and oil saturation log data;
evaluating an objective function measuring the different between the permeability log data and a modeled permeability, the porosity log data and a modeled porosity, and the oil saturation log data and a modeled oil saturation, the modeled permeability, the modeled porosity, and the modeled oil saturation each a function of Thomeer parameters; and
determining the capillary pressures of the reservoir using a Thomeer model having the Thomeer parameters, the Thomeer parameters comprising a fractional bulk volume, a pore geometrical factor, and a minimum entry pressure for each pore system.
33. The computer-implemented method of claim 32 , wherein evaluating the objective function measuring the different between the permeability log data and a modeled permeability, the porosity log data and a modeled porosity, and the oil saturation log data and a modeled oil saturation comprises evaluating the objective function based on a linear equality constraint dependent on the porosity log data, a linear inequality constraint, and a nonlinear equality constraint dependent on the permeability log data.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/902,333 US20140350860A1 (en) | 2013-05-24 | 2013-05-24 | Systems, methods, and computer-readable media for continuous capillary pressure estimation |
CA2910955A CA2910955A1 (en) | 2013-05-24 | 2014-05-23 | Systems, methods, and computer-readable media for continuous capillary pressure estimation |
PCT/US2014/039339 WO2014190261A2 (en) | 2013-05-24 | 2014-05-23 | Systems, methods, and computer-readable media for continuous capillary pressure estimation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/902,333 US20140350860A1 (en) | 2013-05-24 | 2013-05-24 | Systems, methods, and computer-readable media for continuous capillary pressure estimation |
Publications (1)
Publication Number | Publication Date |
---|---|
US20140350860A1 true US20140350860A1 (en) | 2014-11-27 |
Family
ID=51022442
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/902,333 Abandoned US20140350860A1 (en) | 2013-05-24 | 2013-05-24 | Systems, methods, and computer-readable media for continuous capillary pressure estimation |
Country Status (3)
Country | Link |
---|---|
US (1) | US20140350860A1 (en) |
CA (1) | CA2910955A1 (en) |
WO (1) | WO2014190261A2 (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016164507A1 (en) * | 2015-04-09 | 2016-10-13 | Schlumberger Technology Corporation | Oilfield reservior saturation and permeability modeling |
WO2016196421A1 (en) * | 2015-06-01 | 2016-12-08 | Schlumberger Technology Corporation | Oilfield reservoir saturation and permeability modeling |
WO2017100565A1 (en) | 2015-12-10 | 2017-06-15 | Saudi Arabian Oil Company | Predicting and modeling changes in capillary pressure and relative permeabilities in a porous medium due to mineral precipitation and dissolution |
US11066905B2 (en) * | 2015-06-30 | 2021-07-20 | Schlumberger Technology Corporation | Oilfield reservoir saturation and permeability modeling |
US11314909B2 (en) | 2018-10-05 | 2022-04-26 | Qatar Foundation For Education, Science And Community Development | Methods and systems for simulating multiphase flow through porous media |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2785321T3 (en) | 2014-12-17 | 2020-10-06 | Reliance Worldwide Corp | Fluid flow meter, method of identifying a pipe type and computer-readable medium |
CN107102377B (en) * | 2017-06-27 | 2018-09-14 | 中国石油大学(华东) | The method of quantitative forecast tight sand favorable oil/gas exploration area |
US11231474B2 (en) * | 2019-04-08 | 2022-01-25 | Baker Hughes Oilfield Operations Llc | Inversion of NMR echo trains using a supplementary nonlinear equality constraint |
CN110826016B (en) * | 2019-11-04 | 2023-03-21 | 中国石油天然气股份有限公司 | Saturation and free water surface calculation and prediction method based on microscopic heterogeneity representation |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9081117B2 (en) * | 2010-09-15 | 2015-07-14 | Baker Hughes Incorporated | Method and apparatus for predicting petrophysical properties from NMR data in carbonate rocks |
US9671525B2 (en) * | 2013-03-08 | 2017-06-06 | Saudi Arabian Oil Company | Determining continuous capillary pressure curves for subsurface earth formations using saturation and NMR log data |
-
2013
- 2013-05-24 US US13/902,333 patent/US20140350860A1/en not_active Abandoned
-
2014
- 2014-05-23 WO PCT/US2014/039339 patent/WO2014190261A2/en active Application Filing
- 2014-05-23 CA CA2910955A patent/CA2910955A1/en not_active Abandoned
Non-Patent Citations (1)
Title |
---|
Xu et al., SATURATION-HEIGHT AND INVASION CONSISTENT HYDRAULIC ROCK TYPING USING MULTI-WELL CONVENTIONAL LOGS, 16-20 June 2012, SPWLA 53rd Annual Logging Symposium, pp 1-16 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016164507A1 (en) * | 2015-04-09 | 2016-10-13 | Schlumberger Technology Corporation | Oilfield reservior saturation and permeability modeling |
WO2016196421A1 (en) * | 2015-06-01 | 2016-12-08 | Schlumberger Technology Corporation | Oilfield reservoir saturation and permeability modeling |
US10787902B2 (en) | 2015-06-01 | 2020-09-29 | Schlumberger Technology Corporation | Method and system for correcting a capillary pressure curve |
US11066905B2 (en) * | 2015-06-30 | 2021-07-20 | Schlumberger Technology Corporation | Oilfield reservoir saturation and permeability modeling |
WO2017100565A1 (en) | 2015-12-10 | 2017-06-15 | Saudi Arabian Oil Company | Predicting and modeling changes in capillary pressure and relative permeabilities in a porous medium due to mineral precipitation and dissolution |
US10049172B2 (en) | 2015-12-10 | 2018-08-14 | Saudi Arabian Oil Company | Predicting and modeling changes in capillary pressure and relative permeabilities in a porous medium due to mineral precipitation and dissolution |
US11314909B2 (en) | 2018-10-05 | 2022-04-26 | Qatar Foundation For Education, Science And Community Development | Methods and systems for simulating multiphase flow through porous media |
Also Published As
Publication number | Publication date |
---|---|
WO2014190261A2 (en) | 2014-11-27 |
WO2014190261A3 (en) | 2015-01-15 |
CA2910955A1 (en) | 2014-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20140350860A1 (en) | Systems, methods, and computer-readable media for continuous capillary pressure estimation | |
Zhong et al. | Predicting CO2 plume migration in heterogeneous formations using conditional deep convolutional generative adversarial network | |
Chai et al. | An efficient method for fractured shale reservoir history matching: The embedded discrete fracture multi-continuum approach | |
Mariethoz et al. | Bayesian inverse problem and optimization with iterative spatial resampling | |
Oladyshkin et al. | An integrative approach to robust design and probabilistic risk assessment for CO 2 storage in geological formations | |
Oladyshkin et al. | Bayesian updating via bootstrap filtering combined with data-driven polynomial chaos expansions: methodology and application to history matching for carbon dioxide storage in geological formations | |
March et al. | Accurate early‐time and late‐time modeling of countercurrent spontaneous imbibition | |
Yang et al. | A fully implicit constraint-preserving simulator for the black oil model of petroleum reservoirs | |
US20150006081A1 (en) | Adaptive time-lapse sub-surface electrical resistivity monitoring | |
Liao et al. | Probabilistic collocation method for strongly nonlinear problems: 1. Transform by location | |
Dostert et al. | Coarse-gradient Langevin algorithms for dynamic data integration and uncertainty quantification | |
Mendes et al. | A new computational strategy for solving two‐phase flow in strongly heterogeneous poroelastic media of evolving scales | |
Bao et al. | Data-driven model reduction based on sparsity-promoting methods for multiphase flow in porous media | |
Lu et al. | Assisted history matching for fractured reservoirs by use of hough-transform-based parameterization | |
Tanaka et al. | Streamline-based history matching of arrival times and bottomhole pressure data for multicomponent compositional systems | |
Man et al. | Sequential ensemble‐based optimal design for parameter estimation | |
WO2020154317A1 (en) | Integration of seismic driven rock property into a geo-cellular model | |
Liao et al. | Probabilistic collocation method for strongly nonlinear problems: 2. Transform by displacement | |
Dai et al. | Spatio-temporal modelling of dam deformation using independent component analysis | |
Ibrahima et al. | An efficient distribution method for nonlinear two-phase flow in highly heterogeneous multidimensional stochastic porous media | |
AU2013400127B2 (en) | Relative permeability inversion from historical production data using viscosity ratio invariant step-function relative permeability approximations | |
US10167703B2 (en) | Optimal well placement under constraints | |
Nagao et al. | Data-Driven Discovery of Physics for Reservoir Surveillance | |
Ginting et al. | Rapid quantification of uncertainty in permeability and porosity of oil reservoirs for enabling predictive simulation | |
US20210349226A1 (en) | Construction of a high-resolution advanced 3d transient model with multiple wells by integrating pressure transient data into static geological model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SAUDI ARABIAN OIL COMPANY, SAUDI ARABIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:MEZGHANI, MOKHLES M.;AL IBRAHIM, MUSTAFA;BUITING, JOHANNES JACOBUS MARIA;SIGNING DATES FROM 20130605 TO 20130606;REEL/FRAME:030590/0529 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- AFTER EXAMINER'S ANSWER OR BOARD OF APPEALS DECISION |