CA2442596A1 - Method for enhancing production allocation in an integrated reservoir and surface flow system - Google Patents
Method for enhancing production allocation in an integrated reservoir and surface flow system Download PDFInfo
- Publication number
- CA2442596A1 CA2442596A1 CA002442596A CA2442596A CA2442596A1 CA 2442596 A1 CA2442596 A1 CA 2442596A1 CA 002442596 A CA002442596 A CA 002442596A CA 2442596 A CA2442596 A CA 2442596A CA 2442596 A1 CA2442596 A1 CA 2442596A1
- Authority
- CA
- Canada
- Prior art keywords
- objective function
- constraint
- wellbores
- fluid flow
- recalculating
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 100
- 230000002708 enhancing effect Effects 0.000 title claims abstract description 6
- 239000012530 fluid Substances 0.000 claims abstract description 103
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 38
- 230000035945 sensitivity Effects 0.000 claims description 13
- 238000009826 distribution Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000011084 recovery Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 abstract description 15
- 239000007789 gas Substances 0.000 description 22
- 239000003921 oil Substances 0.000 description 16
- 239000003208 petroleum Substances 0.000 description 13
- 238000004088 simulation Methods 0.000 description 12
- 239000007788 liquid Substances 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 239000000243 solution Substances 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000002347 injection Methods 0.000 description 4
- 239000007924 injection Substances 0.000 description 4
- 239000013049 sediment Substances 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000037230 mobility Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/14—Obtaining from a multiple-zone well
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
Abstract
A method for enhancing allocation of fluid flow rates among a plurality of wellbores coupled to surface facilities is disclosed. The method includes modeling fluid flow characteristics of the wellbores and reservoirs penetrated by the wellbores. The method includes modeling fluid flow characteristics of the surface facilities (Fig. 2, item 42). An optimizer (44) adapted to determine an enhanced value of an objective function corresponding to the modeled fluid flow characteristics of the wellbores and the surface facilities is then operated. The objective function relates to at least one production system performance parameter. Fluid flow rates are then allocated according to the optimization ( item 60).
Description
METHOD FOR ENHANCING PRODUCTION ALLOCATION IN
AN INTEGRATED RESERVOIR AND SURFACE FLOW SYSTEM
CROSS-REFERENCE TO RELATED APPLICATION
This application claims priority benefit from U.S. provisional patent application number 601286,134 filed April 24, 2001.
FIELD OF THE INVENTION
The invention relates generally to the field of petroleum production equipment and production control systems. More specifically, the invention relates to methods and systems for controlling production from a plurality of petroleum wells and reservoirs coupled to a limited number of surface facilities so as to enhance use of the facilities and production from the reservoirs.
BACKGROUND OF THE INVENTION
Petroleum is generally produced by drilling wellbores through permeable earth formations having petroleum reservoirs therein, and causing petroleum fluids in the reservoir to move to the earth's surface through the wellbores. Movement is accomplished by creating a pressure difference between the reservoir and the wellbore. Produced fluids from the wells may include various quantities of crude oil, natural gas and/or water, depending on the conditions in the particular reservoir being produced. Depending on conditions in the particular reservoir, the amounts and rates at which the various fluids will be extracted from a particular well depend on factors which include pressure difference between the reservoir and the wellbore. As is known in the art, wellbore pressure may be adjusted by operating various devices such as chokes (orifices) disposed in the fluid flow path along the wellbore, pumps, compressors, fluid injection devices (which pump fluid into a reservoir to increase its pressure). Generally speaking, changing the rate at which a total volume of fluid is extracted from any particular wellbore may also affect relative rates at which oil, water and gas are produced from each wellbore.
AN INTEGRATED RESERVOIR AND SURFACE FLOW SYSTEM
CROSS-REFERENCE TO RELATED APPLICATION
This application claims priority benefit from U.S. provisional patent application number 601286,134 filed April 24, 2001.
FIELD OF THE INVENTION
The invention relates generally to the field of petroleum production equipment and production control systems. More specifically, the invention relates to methods and systems for controlling production from a plurality of petroleum wells and reservoirs coupled to a limited number of surface facilities so as to enhance use of the facilities and production from the reservoirs.
BACKGROUND OF THE INVENTION
Petroleum is generally produced by drilling wellbores through permeable earth formations having petroleum reservoirs therein, and causing petroleum fluids in the reservoir to move to the earth's surface through the wellbores. Movement is accomplished by creating a pressure difference between the reservoir and the wellbore. Produced fluids from the wells may include various quantities of crude oil, natural gas and/or water, depending on the conditions in the particular reservoir being produced. Depending on conditions in the particular reservoir, the amounts and rates at which the various fluids will be extracted from a particular well depend on factors which include pressure difference between the reservoir and the wellbore. As is known in the art, wellbore pressure may be adjusted by operating various devices such as chokes (orifices) disposed in the fluid flow path along the wellbore, pumps, compressors, fluid injection devices (which pump fluid into a reservoir to increase its pressure). Generally speaking, changing the rate at which a total volume of fluid is extracted from any particular wellbore may also affect relative rates at which oil, water and gas are produced from each wellbore.
Production processing equipment, known by a general term "surface facilities", includes various devices to separate oil and water in liquid form from gas in the produced petroleum. Extracted liquids may be temporarily stored or may be moved to a pipeline for transportation away from the location of the wellbore.
Gas may be transported by pipeline to a point of sale, or may be transported by pipe for further processing away from the location of the wellbore. The surface facilities are typically designed to process selected volumes or quantities of produced petroleum.
The selected volumes depend on what is believed to be likely volumes of production from various wellbores, and how many wellbores are to be coupled to a particular set of surface facilities. Depending on the physical location of the reservoir, such as below the ocean floor or other remote location, it is often economically advantageous to couple a substantial number of wells, and typically from a plurality of different reservoirs, to a single set of surface facilities. As for less complicated installations, the surface facilities coupled to multiple wells and reservoirs are typically selected to most efficiently process expected quantities of the various fluids produced from the wells. An important aspect of the economic performance of surface facilities' is appropriate selection of sizes and capacities of various components of the surface facilities. Equipment which is too small for actual quantities of fluids produced may limit the rate at which the various wellbores may be produced. Such condition may result in poor economic performance of the entire reservoir and surface facility combination. Conversely, equipment which has excess capacity may increase capital costs beyond those necessary, reducing overall rate of return on investment.
Still another problem in the efficient use of surface facilities can arise when some wellbores change fluid production rates. As is known in the art, such changes in rate may result from natural depletion of the reservoir, and from unforeseen problems with one or more wellbores in a reservoir, among others. Sometimes, it is possible to change production rates in other wellbores coupled to the surface facilities to maintain throughput in the surface facilities. As is known in the art, however, such production rate changes may be accompanied by changes in relative quantities of water, oil and gas produced from the affected wellbores. Such relative rate changes may affect the ability of the surface facilities to operate efficiently.
One way to determine expected quantities of produced fluids from each wellbore in each reservoir is to mathematically simulate the performance of each well in each reservoir to be coupled to the surface facilities. Typically this mathematical simulation is performed using a computer program. Such reservoir simulation computer programs are well known in the art. Reservoir simulation programs, however, typically do not include any means to couple the simulation result to a simulation of the operation of surface facilities. Therefore, there is no direct linkage between selective operation of the various wellbores and whether the surface facilities are being operated in an optimal way.
One system that attempts to couple reservoir simulation with surface facility simulation is described in, G. G. Hepguler et al, Integration of a field surface and production network with a reservoir simulator, SPE Computer Appl. vol. 9, p.
~~, Society of Petroleum Engineers, Richardson, TX (1997). A limitation to the system described in the Hepguler et al reference is that it is unable to generate a corrective action with respect to the surface facilities which may arise out of infeasibility.
Infeasibility is defined as the production system operating outside a constraint or limit, for example, defining a maximum allowable water production which is lower than an expected water production from reservoir simulation. Another limitation in the Hepulger et al system is that there is poor convergence in an optimization routine in the system. Other prior art optimization systems are described, for example in M.
R. Palke et al, Nonlinear optimization of well production considering gas lift and phase behavior, Proceedings, SPE production operations symposium, p. 341, Society of Petroleum Engineers, Richardson, TX (1995). This reference deals primarily with optimizing gas lift systems and does not describe any means for optimizing surface facility use in conjunction with optimizing reservoir production.
A method for optimizing production allocation between wellbores in a reservoir is described in, Zakirov et al, Optmizing reservoir performance by automatic allocation of well fates, Conference Proceedings, 5th Math of Oil Recovery, Europe, p. 375 (1996). The method described in this reference does not deal with optimizing the use of surface facilities in conjunction with optimizing reservoir production.
It is desirable to have a simulation system that can enhance or optimize, both reservoir production and surface facility operation simultaneously, while also being able to assist in isolating and rectifying causes of the production system operating outside constraints.
SUMMARY OF THE INVENTION
The invention generally is a method for enhancing allocation of fluid flow rates among a plurality of wellbores coupled to surface facilities. The method includes modeling fluid flow characteristics of the wellbores and reservoirs penetrated by the wellbores. The method includes modeling fluid flow characteristics of the surface facilities. An optimizer adapted to determine an optimal value of an objective function corresponding to the modeled fluid flow characteristics of the wellbores and the surface facilities is then operated. The objective function relates to at least one production system performance parameter. Fluid flow rates are then allocated among the plurality of wellbores as determined by the operating the optimizer.
In some embodiments, a constraint on the system is adjusted. The optimizer is again operated using the adjusted constraint. This is repeated until an enhanced fluid flow rate allocation is determined.
In some embodiments, non-convergence of the optimizer is determined. At least one system constraint is adjusted and the optimizer is again operated.
This is repeated until the optimizer converges.
In some embodiments, the optimizer includes successive quadratic programming. A value of a Lagrange multiplier associated with at least one system constraint is determined as a result of the successive quadratic programming.
The value of the Lagrange multiplier can be used to determine a sensitivity of the production system to the at least one constraint.
Other aspects and advantages of the invention will be apparent from the following description and the appended claims.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 shows an example of a plurality of wellbores coupled to various surface facilities.
Figure 2 is a flow chart showing operation of one embodiment of the invention.
DETAILED DESCRIPTION
Figure 1 shows one example of a petroleum production system. The production system in Figure 1 includes a plurality of wellbores W, which rnay penetrate the same reservoir, or a plurality of different subsurface petroleum reservoirs (not shown). The wellbores W are coupled in any manner known in the art to various surface facilities. Each wellbore W may be coupled to the various surface facilities using a flow control device C, such as a controllable choke, or similar fixed or variable flow restriction, in the fluid coupling between each wellbore W
and the surface facilities. The flow control device C may be locally or remotely operable.
The surface facilities may include, for example, production gathering platforms 22, 24, 26, 28, 30, 32 and 33, where production from one or more of the wellbores W may be collected, stored, commingled andlor remotely controlled.
Control in this context means having a fluid flow rate from each wellbore W
selectively adjusted or stopped. Fluid produced from each of the wellbores W
is coupled directly, or commingled with produced fluids from selected other ones of the wellbores W, to petroleum fluid processing devices which may include separators S.
The separators S may be of any type known in the art, and are generally used to separate gas, oil and sediment and water from the fluid extracted from the wellbores W. Each separator S may have a gas output 13, and outputs for liquid oil 10 and for water and sediment 12. The liquid oil 10 and water and sediment 12 outputs may be coupled to storage units or tanks (not shown) disposed on one or more of the platforms 22, 24, 26, 28, 30, 32 and 33, or the liquid outputs 10, 12 may be coupled to a pipeline (not shown) for transportation to a location away from the wellbore W
locations or the platforms 22, 24, 26, 28, 30, 32 and 33. The gas outputs 13 may be coupled directly, or commingled at one of the platforms, for example platform 26, to serial-connected compressors 14, 16, then to a terminal 18 for transport to a sales line (not shown) or to a gas processing plant 20, which may itself be on a platform or at a remote physical location. Gas processing plants are known in the art for removing impurities and gas liquids from "separated" gas (gas that is extracted from a device such as one of the separators S). Any one or all of the platforms 22, 24, 26, 28, 30, 32 and 33 may also include control devices (not shown) for regulating the total amount of fluid, including gas, delivered from the respective platform to the separator S, to the pipeline (not shown) or to the compressors 14, 16. It should be clearly noted that the production system shown in Figure 1 is only an example of the types of production systems and elements thereof than can be used with the method of the invention. The method of the invention only requires that the fluid flow characteristics of each component in any production system be able to be modeled or characterized so as to be representable by an equation or set of equations.
"Component" in this context means both the wellbores W and one or more components of the surface facilities. Accordingly, the invention is not intended to be limited to use with a production system that includes or excludes any one or more of the components of the system shown in Figure 1.
In a production system, such as the one shown in Figure 1, as some of the wellbores W are operated to extract particular amounts (at selected rates) of fluid from the one or more subsurface reservoirs (not shown), various quantities of gas, oil and/or water will flow into these wellbores W at rates which may be estimated by solution to reservoir mass and momentum balance equations. Such mass and _7_ momentum balance equations are well known in the art for estimating wellbore production. The fluid flow rates depend on relative fluid mobilities in the subsurface reservoir and on the pressure difference between the particular one of the wellbores W
and the reservoir (not shown). As is known in the art, as any one or more of the wellbores W is selectively controlled, such as by operating its associated flow control device C, the rates at which the various fluids are produced from each such wellbore W will change, both instantaneously and over time. The change over time, as is known in the art, is related to the change in pressure and fluid content distribution in the reservoir as fluids are extracted at known rates. These changes in fluid flow rates may also be calculated using mass and momentum balance equations known in the art. Such changes in fluid flow rates will have an effect on operation of the various components of the surface facilities, including for example, the compressors 14, 16, and the separators S. As will be further explained, a method according to the invention seeks to optimize one or more selected production system performance parameters with respect to both fluid extracted from the one or more subsurface reservoirs (not shown) and with respect to operation of the surface facilities.
It should be noted that in the example production system of Figure 1, any one or more of the wellbores W may be an injector well, meaning that fluid is not extracted from that wellbore, but that the fluid is pumped into that wellbore.
Fluid pumping into a wellbore, as is known in the axt, is generally either for disposal of fluid or for providing pressure to the subsurface reservoir (not shown). As a practical matter, the only difference between an injector well (where injection is into one of the reservoirs) and a producing (fluid extracting) wellbore is that for reservoir simulation purposes, an injector well will act as a source of pressure into the reservoir, rather than a pressure sink from the reservoir.
One aspect of the invention is to determine an allocation of fluid flow rates from each of the wellbores W in the production system so that a particular production performance parameter is optimized. The production performance parameter may be, for example, maximization' of oil production, minimization of gas and/or water _g_ production, or maximizing an economic value of the entire production system, such as by net present value or similar measure of value, or maximizing an ultimate oil or gas recovery from the one or more subsurface reservoirs (not shown). It should be noted that the foregoing are only examples of production performance parameters and that the invention is not limited to the foregoing parameters as the performance parameter wluch is to be enhanced or optimized.
In a method according to this aspect of the invention, fluid flow allocation is modeled mathematically by a non-linear optimization procedure. The non-linear optimization includes an objective function and a set of inequality and equality constraints. The objective function can be expressed as:
F - ~~k~k(W x) The objective function is subject to the following equality constraints represented by the expressions:
H(w,x) = 0 which represents the subsurface reservoir mass and momentum balance equations and S(w,x) = 0 which represents the surface facilities flow and pressure balance equations.
The objective function is also subject to inequality constraints:
a<_C(w,x)<_b where w represents subsurface reservoir variables such as fluid component mole number, fluid pressure, temperature, etc. x represents "decision" variables such as pressure in any wellbore W at the depth of the subsurface reservoir (known as "bottom hole pressure" - BHP), pressure at any surface "node" (a connection between any two elements of the surface facilities), and a and b represent lower and upper boundaries for each of the constraints C . Constraints may include system operating parameters such as gas/oil ratio (GOR), flow rate, pressure, water cut (fractional amount of produced liquid consisting of water), or any similar parameter which is affected by changing the fluid flow rate out of any of the wellbores W, or by changing any operating parameter of any element of the surface facilities, such as separators S
or compressors 14, 16. ' Variable wk in the above objective function represents a set of weighting factors, which can be applied individually to individual contribution variables, irk, in the objective function. The individual contribution variables may include flow rates of the various fluids from each of the wellbores W, although the individual contribution variables are not limited to flow rates. As previously explained, the flow rates can be calculated using well known mass and momentum balance equations.
In a method according to this aspect of the invention, any one of the wellbores W
or any surface device, including but not limited to the separators S and/or compressors 14, 16 may be represented as one of the reservoir variables or one of the decision variables.
Similarly, the objective function can be arranged to include any configuration of wellbores and surface facilities.
The ones of the constraints C which represent selected ("target") values of fluid production rates for the system, such as total water flow rate, GOR, or oil flow rate, for example, are preferably inequality constraints with the target values set as an upper or lower boundary, as is consistent with the particular target. Doing this enables the optimizer to converge under conditions where the actual system production rate is different from the target, but does not fall outside the limit set by the target.
An optimization system according to the invention enables production allocation with respect to a production performance parameter that includes reservoir variables in the calculation. Prior art systems that attempt to couple reservoir simulation with surface facility simulation, for example the one described in, G. G.
Hepguler et al, I~cteg~ation of a field surface and production hetwo~k with a reservoir simulator, SPE Computer Appl. vol. 9, p. 88, Society of Petroleum Engineers, Richardson, TX (1997) [referred to in the Background section herein], do not seek to optimize production allocation and reservoir calculations in a single executable program. One advantage that may be offered by a system according to the invention is a substantial saving in computation time.
In one embodiment of a method according to the invention, the objective function can be optimized by using successive quadratic programming (SQP). In SQP, the objective function is approximated as a quadratic function, and constraints are linearized. The SQP algorithm used in embodiments of the invention can be described as follows. Consider a general nonlinear optimization problem of the form:
Minimize F(x) x E Rn (1) subject to constraints h;(x) = 0 i = 1, .... , neq (2) g~(x) <_ 0 j = 1, .... , n;eq (3) If g~(x) = 0 then the constraint is active while the constraint is inactive if g~(x) < 0. A Lagrange function L(x, u, v) is defined so that:
L(x, u, v) ---- F(x) + ~ u1 h~ (x) + ~ v~ g~ (x) (4) minimizing L(x, u, v) also minimizes F(x) subject to the above constraints.
Here u;
and v~ represent the Lagrange multiplier for equality constraint i and inequality constraint j, respectively. v~ > 0 for active constraints, while v~ = 0 when the constraint is inactive. It can be shown that the following conditions are satisfied at the optimum:
OL(x, u, v) = OF(x) + ~ u;~h; (x) + ~ v~ ~g~ (x) = o (s) uth~ (x) = 0 (6) v~ g~ (x) = 0 (7) v; >- 0 (8) These conditions are called Karesh-Kuhn-Tucker (KKT) optimality criteria. It can be shown that applying Newton's method to solve the optimality criteria for the problem described in equations (1) - (4) is equivalent to solving the following quadratic problem:
Minimize OF(xo )Ox + z ~xT H(xo )fix (9) g(xo)+Og(xo)~ c 0 (10) h(xo) + Vh(xo )~ = 0 (11) where x0 represents the current guess or estimate as to the actual minimum value of the objective function, and H(x~ represents the Hessian at xo.
Here, as previously explained, the obj ective function is approximated quadratically wlule the constraints are linearly approximated. The minimum found for this approximate problem would be exact if the Hessian, (H(xo)), is also exact.
However, an inexact Hessian can be used in the foregoing formulation to save computation cost. By applying the above quadratic approximation successively, the real minimum of the objective function is obtained at convergence.
The terms "optimize" and "optimizing" as used with respect to this invention are intended to mean to determine or determining, respectively, an apparent optimum value of the objective function. As will be appreciated by those skilled in the art, in certain circumstances a localized optimum value of the objective function may be determined during any calculation procedure which seeks to determine the true ("global") optimum value of the objective function. Accordingly, the terms "optimize" and "optimizing" are intended to include within their scope any calculation procedure which seeks to determine an enhanced or optimum value of the objective function. Any allocation of fluid flow rates and/or surface facility operating parameters which result from such calculation procedure, whether the global optimum or a localized optimum value of the objective function is actually determined, are therefore also within the scope of this invention. In some instances, as will be readily appreciated by those skilled in the art, it may be desirable for a production system operator to intentionally select a fluid flow rate allocation among the wellbores that is less than optimal as determined by the optimizer. Accordingly, the invention shall not be limited in scope only to determining an optimal fluid flow rate allocation as a result of operating an optimization program according to the various embodiments of the invention.
In a particular embodiment of the invention, the Lagrange multipliers defined in equation (4) can be used to determine a sensitivity of the optimizer to any or all of the optimizer constraints. The values of one or more of the Lagrange multipliers are a measure of the sensitivity of the objective function to the associated constraints. The measure of sensitivity can be used to determine which of the constraints may be relaxed or otherwise adjusted to provide a substantial increase in the value of the system performance parameter that is to be optimized. As an example, a selected maximum total system water production may be a "bottleneck" to total oil production.
During optimization, the Lagrange multiplier associated with the maximum total system water production may indicate that a slight relaxation or adjustment of the selected maximum water production rate may provide the production system with the capacity to substantially increase maximum oil production rate, and correspondingly, the economic value (for example, net present value) of the production system.
The foregoing is meant to serve only as one example of use of the Lagrange multipliers calculated by the optimizer to determine constraint sensitivity. Any other constraint used in the optimizer may also undergo similar sensitivity analysis to determine production system "bottlenecks".
In one embodiment of a method according to the invention, a so-called "infeasible path" strategy is used, where the initial estimate or guess (x0) is allowed to be infeasible. "Infeasible" means that some or all of the constraints and variables are out of their respective minimum or maximum bounds. For example, one or more of the wellbores W may produce water at a rate which exceeds a maximum water production rate target for the entire system, or the total gas production, as another example, may exceed the capacity of the compressors. The optimization algorithm simultaneously tries to reach to an optimum as well as a feasible solution.
Thus feasibility is determined only at convergence. The advantage of this strategy is reduced objective and constraint function evaluation cost. How the infeasible solution strategy of the method of the invention is used will be further explained.
The solution of the optimization problem provides an optimal fluid flow rate and pressure distribution within the entire surface facility network. A part of this solution is then used in the reservoir simulator as the boundary conditions, while then solving the mass and momentum balance equations that describe the fluid flow in the reservoir.
A flow chart of how an optimization method according to the invention can be used in operating a production system is shown in Figure 2. After surface facility equations and reservoir equations are set up, and initial conditions in the surface facility and reservoir are set, at 40 the system time is incremented. If any surface facility operating parameters or structures have been changed from the previous calculation, shown at 42, such changes are entered into the conditions and/or equations for the surface facilities and reservoir. At 44, the conditions and constraints are entered into an optimization routine as previously described. At 46, the optimizer it is determined as to whether the optimizer has reached convergence. As previously explained, when the optimizer reaches convergence, an optimal value of the objective function is determined. When the optimal value of the objective function is determined, the system performance parameter which is represented by the objective function is at an optimal value. As previously explained, the performance parameter can be, for example, economic value, maximum oil production, minimized gas and/or water production, minimum operating cost, or any other parameter related to a measure of production and/or economic performance of the production system such as shown in Figure 1. The result of the optimization is an allocation of fluid production rates from each of the wellbores (W in Figure 1) which results in the optimization of the selected system performance parameters.
Referring again to Figure 2, the output of the optimizer includes fluid production rate allocation among the wellbores in the production system. In actual production and/or injection at the rates allocated by the optimizer, each wellbore (W
in Figure 1) will cause a pressure sink or pressure increase (depending on whether the wellbore is a producing well or injection well) at the reservoir. Such pressure changes propagate through the reservoir, and these pressure changes can be calculated using the mass and momentum balance equations referred to earlier. Therefore, as fluids are produced or injected into each wellbore W, a distribution of conditions in the subsurface reservoir changes. Using the output of the optimizer, the set of fluid flow rates for each wellbore as a set of boundary conditions, as shown at 62, a new distribution of conditions (particularly including but not limited to pressure) for the subsurface reservoir is calculated, at 64.
In some instances, the changes in reservoir conditions will result in changes in fluid flow rates from one or more of the wellbores (W in Figure 1). As these changes take place, they become part of the initial conditions for operating the optimizer, as indicated in Figure 2 by a line leading back to box 40.
In other cases, the optimizer will not converge. Failure of convergence, as explained earlier with reference to the description of the SQP aspect of the optimizer, is typically because at least one of the constraints is violated. The constraints may include operating parameters such as maximum acceptable water production in the system, maximum GOR, minimum inlet pressure to the compressor (14 in Figure 1), and others. In the event no system fluid production allocation will enable meeting all the constraints, the optimizer will not converge. In another aspect of the invention, a cause of the optimizer failing to converge may lead to isolation of one or more elements of the production system which cause the constraints to be violated.
At box 48 in Figure 2, one or more of the constraints may be relaxed or removed. For example a maximum acceptable water production may be increased, or removed as a constraint, or, alternatively, a minimum oil production may be reduced or removed as a constraint. Then, at box 50, the optimizer is run again. If convergence is achieved, then the violated constraint has been identified, at 52. At 54, corrective action can be taken to repair or correct the violated constraint. For example, if a maximum horsepower rating of the compressor (14 in Figure 1) is exceeded by a selected system gas flow rate, the compressor may be substituted by a higher rating compressor, and the optimizer run again, at 56. Any other physical change to the production system which alters or adjusts a system constraint can be detected and corrected by the method elements outlined in boxes 48, 50, 52 and 54, and the examples referred to herein should not be interpreted as limiting the types of system constraints that can be affected by the method of this invention. At box 58, if the optimizer has converged, then the flow rates axe allocated among the wellbores (W in Figure 1) according to the solution determined by the optimizer. At 60, these fluid flow rates are used as boundary conditions to perform a recalculation of the reservoir conditions, as in the earlier case where the initial run of the optimizer converged (at box 46).
While the invention has been described with respect to a limited number of embodiments, those skilled in the art, having benefit of this disclosure, will appreciate that other embodiments can be devised which do not depart from the scope of the invention as disclosed herein. Accordingly, the scope of the invention should be limited only by the attached claims.
Gas may be transported by pipeline to a point of sale, or may be transported by pipe for further processing away from the location of the wellbore. The surface facilities are typically designed to process selected volumes or quantities of produced petroleum.
The selected volumes depend on what is believed to be likely volumes of production from various wellbores, and how many wellbores are to be coupled to a particular set of surface facilities. Depending on the physical location of the reservoir, such as below the ocean floor or other remote location, it is often economically advantageous to couple a substantial number of wells, and typically from a plurality of different reservoirs, to a single set of surface facilities. As for less complicated installations, the surface facilities coupled to multiple wells and reservoirs are typically selected to most efficiently process expected quantities of the various fluids produced from the wells. An important aspect of the economic performance of surface facilities' is appropriate selection of sizes and capacities of various components of the surface facilities. Equipment which is too small for actual quantities of fluids produced may limit the rate at which the various wellbores may be produced. Such condition may result in poor economic performance of the entire reservoir and surface facility combination. Conversely, equipment which has excess capacity may increase capital costs beyond those necessary, reducing overall rate of return on investment.
Still another problem in the efficient use of surface facilities can arise when some wellbores change fluid production rates. As is known in the art, such changes in rate may result from natural depletion of the reservoir, and from unforeseen problems with one or more wellbores in a reservoir, among others. Sometimes, it is possible to change production rates in other wellbores coupled to the surface facilities to maintain throughput in the surface facilities. As is known in the art, however, such production rate changes may be accompanied by changes in relative quantities of water, oil and gas produced from the affected wellbores. Such relative rate changes may affect the ability of the surface facilities to operate efficiently.
One way to determine expected quantities of produced fluids from each wellbore in each reservoir is to mathematically simulate the performance of each well in each reservoir to be coupled to the surface facilities. Typically this mathematical simulation is performed using a computer program. Such reservoir simulation computer programs are well known in the art. Reservoir simulation programs, however, typically do not include any means to couple the simulation result to a simulation of the operation of surface facilities. Therefore, there is no direct linkage between selective operation of the various wellbores and whether the surface facilities are being operated in an optimal way.
One system that attempts to couple reservoir simulation with surface facility simulation is described in, G. G. Hepguler et al, Integration of a field surface and production network with a reservoir simulator, SPE Computer Appl. vol. 9, p.
~~, Society of Petroleum Engineers, Richardson, TX (1997). A limitation to the system described in the Hepguler et al reference is that it is unable to generate a corrective action with respect to the surface facilities which may arise out of infeasibility.
Infeasibility is defined as the production system operating outside a constraint or limit, for example, defining a maximum allowable water production which is lower than an expected water production from reservoir simulation. Another limitation in the Hepulger et al system is that there is poor convergence in an optimization routine in the system. Other prior art optimization systems are described, for example in M.
R. Palke et al, Nonlinear optimization of well production considering gas lift and phase behavior, Proceedings, SPE production operations symposium, p. 341, Society of Petroleum Engineers, Richardson, TX (1995). This reference deals primarily with optimizing gas lift systems and does not describe any means for optimizing surface facility use in conjunction with optimizing reservoir production.
A method for optimizing production allocation between wellbores in a reservoir is described in, Zakirov et al, Optmizing reservoir performance by automatic allocation of well fates, Conference Proceedings, 5th Math of Oil Recovery, Europe, p. 375 (1996). The method described in this reference does not deal with optimizing the use of surface facilities in conjunction with optimizing reservoir production.
It is desirable to have a simulation system that can enhance or optimize, both reservoir production and surface facility operation simultaneously, while also being able to assist in isolating and rectifying causes of the production system operating outside constraints.
SUMMARY OF THE INVENTION
The invention generally is a method for enhancing allocation of fluid flow rates among a plurality of wellbores coupled to surface facilities. The method includes modeling fluid flow characteristics of the wellbores and reservoirs penetrated by the wellbores. The method includes modeling fluid flow characteristics of the surface facilities. An optimizer adapted to determine an optimal value of an objective function corresponding to the modeled fluid flow characteristics of the wellbores and the surface facilities is then operated. The objective function relates to at least one production system performance parameter. Fluid flow rates are then allocated among the plurality of wellbores as determined by the operating the optimizer.
In some embodiments, a constraint on the system is adjusted. The optimizer is again operated using the adjusted constraint. This is repeated until an enhanced fluid flow rate allocation is determined.
In some embodiments, non-convergence of the optimizer is determined. At least one system constraint is adjusted and the optimizer is again operated.
This is repeated until the optimizer converges.
In some embodiments, the optimizer includes successive quadratic programming. A value of a Lagrange multiplier associated with at least one system constraint is determined as a result of the successive quadratic programming.
The value of the Lagrange multiplier can be used to determine a sensitivity of the production system to the at least one constraint.
Other aspects and advantages of the invention will be apparent from the following description and the appended claims.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 shows an example of a plurality of wellbores coupled to various surface facilities.
Figure 2 is a flow chart showing operation of one embodiment of the invention.
DETAILED DESCRIPTION
Figure 1 shows one example of a petroleum production system. The production system in Figure 1 includes a plurality of wellbores W, which rnay penetrate the same reservoir, or a plurality of different subsurface petroleum reservoirs (not shown). The wellbores W are coupled in any manner known in the art to various surface facilities. Each wellbore W may be coupled to the various surface facilities using a flow control device C, such as a controllable choke, or similar fixed or variable flow restriction, in the fluid coupling between each wellbore W
and the surface facilities. The flow control device C may be locally or remotely operable.
The surface facilities may include, for example, production gathering platforms 22, 24, 26, 28, 30, 32 and 33, where production from one or more of the wellbores W may be collected, stored, commingled andlor remotely controlled.
Control in this context means having a fluid flow rate from each wellbore W
selectively adjusted or stopped. Fluid produced from each of the wellbores W
is coupled directly, or commingled with produced fluids from selected other ones of the wellbores W, to petroleum fluid processing devices which may include separators S.
The separators S may be of any type known in the art, and are generally used to separate gas, oil and sediment and water from the fluid extracted from the wellbores W. Each separator S may have a gas output 13, and outputs for liquid oil 10 and for water and sediment 12. The liquid oil 10 and water and sediment 12 outputs may be coupled to storage units or tanks (not shown) disposed on one or more of the platforms 22, 24, 26, 28, 30, 32 and 33, or the liquid outputs 10, 12 may be coupled to a pipeline (not shown) for transportation to a location away from the wellbore W
locations or the platforms 22, 24, 26, 28, 30, 32 and 33. The gas outputs 13 may be coupled directly, or commingled at one of the platforms, for example platform 26, to serial-connected compressors 14, 16, then to a terminal 18 for transport to a sales line (not shown) or to a gas processing plant 20, which may itself be on a platform or at a remote physical location. Gas processing plants are known in the art for removing impurities and gas liquids from "separated" gas (gas that is extracted from a device such as one of the separators S). Any one or all of the platforms 22, 24, 26, 28, 30, 32 and 33 may also include control devices (not shown) for regulating the total amount of fluid, including gas, delivered from the respective platform to the separator S, to the pipeline (not shown) or to the compressors 14, 16. It should be clearly noted that the production system shown in Figure 1 is only an example of the types of production systems and elements thereof than can be used with the method of the invention. The method of the invention only requires that the fluid flow characteristics of each component in any production system be able to be modeled or characterized so as to be representable by an equation or set of equations.
"Component" in this context means both the wellbores W and one or more components of the surface facilities. Accordingly, the invention is not intended to be limited to use with a production system that includes or excludes any one or more of the components of the system shown in Figure 1.
In a production system, such as the one shown in Figure 1, as some of the wellbores W are operated to extract particular amounts (at selected rates) of fluid from the one or more subsurface reservoirs (not shown), various quantities of gas, oil and/or water will flow into these wellbores W at rates which may be estimated by solution to reservoir mass and momentum balance equations. Such mass and _7_ momentum balance equations are well known in the art for estimating wellbore production. The fluid flow rates depend on relative fluid mobilities in the subsurface reservoir and on the pressure difference between the particular one of the wellbores W
and the reservoir (not shown). As is known in the art, as any one or more of the wellbores W is selectively controlled, such as by operating its associated flow control device C, the rates at which the various fluids are produced from each such wellbore W will change, both instantaneously and over time. The change over time, as is known in the art, is related to the change in pressure and fluid content distribution in the reservoir as fluids are extracted at known rates. These changes in fluid flow rates may also be calculated using mass and momentum balance equations known in the art. Such changes in fluid flow rates will have an effect on operation of the various components of the surface facilities, including for example, the compressors 14, 16, and the separators S. As will be further explained, a method according to the invention seeks to optimize one or more selected production system performance parameters with respect to both fluid extracted from the one or more subsurface reservoirs (not shown) and with respect to operation of the surface facilities.
It should be noted that in the example production system of Figure 1, any one or more of the wellbores W may be an injector well, meaning that fluid is not extracted from that wellbore, but that the fluid is pumped into that wellbore.
Fluid pumping into a wellbore, as is known in the axt, is generally either for disposal of fluid or for providing pressure to the subsurface reservoir (not shown). As a practical matter, the only difference between an injector well (where injection is into one of the reservoirs) and a producing (fluid extracting) wellbore is that for reservoir simulation purposes, an injector well will act as a source of pressure into the reservoir, rather than a pressure sink from the reservoir.
One aspect of the invention is to determine an allocation of fluid flow rates from each of the wellbores W in the production system so that a particular production performance parameter is optimized. The production performance parameter may be, for example, maximization' of oil production, minimization of gas and/or water _g_ production, or maximizing an economic value of the entire production system, such as by net present value or similar measure of value, or maximizing an ultimate oil or gas recovery from the one or more subsurface reservoirs (not shown). It should be noted that the foregoing are only examples of production performance parameters and that the invention is not limited to the foregoing parameters as the performance parameter wluch is to be enhanced or optimized.
In a method according to this aspect of the invention, fluid flow allocation is modeled mathematically by a non-linear optimization procedure. The non-linear optimization includes an objective function and a set of inequality and equality constraints. The objective function can be expressed as:
F - ~~k~k(W x) The objective function is subject to the following equality constraints represented by the expressions:
H(w,x) = 0 which represents the subsurface reservoir mass and momentum balance equations and S(w,x) = 0 which represents the surface facilities flow and pressure balance equations.
The objective function is also subject to inequality constraints:
a<_C(w,x)<_b where w represents subsurface reservoir variables such as fluid component mole number, fluid pressure, temperature, etc. x represents "decision" variables such as pressure in any wellbore W at the depth of the subsurface reservoir (known as "bottom hole pressure" - BHP), pressure at any surface "node" (a connection between any two elements of the surface facilities), and a and b represent lower and upper boundaries for each of the constraints C . Constraints may include system operating parameters such as gas/oil ratio (GOR), flow rate, pressure, water cut (fractional amount of produced liquid consisting of water), or any similar parameter which is affected by changing the fluid flow rate out of any of the wellbores W, or by changing any operating parameter of any element of the surface facilities, such as separators S
or compressors 14, 16. ' Variable wk in the above objective function represents a set of weighting factors, which can be applied individually to individual contribution variables, irk, in the objective function. The individual contribution variables may include flow rates of the various fluids from each of the wellbores W, although the individual contribution variables are not limited to flow rates. As previously explained, the flow rates can be calculated using well known mass and momentum balance equations.
In a method according to this aspect of the invention, any one of the wellbores W
or any surface device, including but not limited to the separators S and/or compressors 14, 16 may be represented as one of the reservoir variables or one of the decision variables.
Similarly, the objective function can be arranged to include any configuration of wellbores and surface facilities.
The ones of the constraints C which represent selected ("target") values of fluid production rates for the system, such as total water flow rate, GOR, or oil flow rate, for example, are preferably inequality constraints with the target values set as an upper or lower boundary, as is consistent with the particular target. Doing this enables the optimizer to converge under conditions where the actual system production rate is different from the target, but does not fall outside the limit set by the target.
An optimization system according to the invention enables production allocation with respect to a production performance parameter that includes reservoir variables in the calculation. Prior art systems that attempt to couple reservoir simulation with surface facility simulation, for example the one described in, G. G.
Hepguler et al, I~cteg~ation of a field surface and production hetwo~k with a reservoir simulator, SPE Computer Appl. vol. 9, p. 88, Society of Petroleum Engineers, Richardson, TX (1997) [referred to in the Background section herein], do not seek to optimize production allocation and reservoir calculations in a single executable program. One advantage that may be offered by a system according to the invention is a substantial saving in computation time.
In one embodiment of a method according to the invention, the objective function can be optimized by using successive quadratic programming (SQP). In SQP, the objective function is approximated as a quadratic function, and constraints are linearized. The SQP algorithm used in embodiments of the invention can be described as follows. Consider a general nonlinear optimization problem of the form:
Minimize F(x) x E Rn (1) subject to constraints h;(x) = 0 i = 1, .... , neq (2) g~(x) <_ 0 j = 1, .... , n;eq (3) If g~(x) = 0 then the constraint is active while the constraint is inactive if g~(x) < 0. A Lagrange function L(x, u, v) is defined so that:
L(x, u, v) ---- F(x) + ~ u1 h~ (x) + ~ v~ g~ (x) (4) minimizing L(x, u, v) also minimizes F(x) subject to the above constraints.
Here u;
and v~ represent the Lagrange multiplier for equality constraint i and inequality constraint j, respectively. v~ > 0 for active constraints, while v~ = 0 when the constraint is inactive. It can be shown that the following conditions are satisfied at the optimum:
OL(x, u, v) = OF(x) + ~ u;~h; (x) + ~ v~ ~g~ (x) = o (s) uth~ (x) = 0 (6) v~ g~ (x) = 0 (7) v; >- 0 (8) These conditions are called Karesh-Kuhn-Tucker (KKT) optimality criteria. It can be shown that applying Newton's method to solve the optimality criteria for the problem described in equations (1) - (4) is equivalent to solving the following quadratic problem:
Minimize OF(xo )Ox + z ~xT H(xo )fix (9) g(xo)+Og(xo)~ c 0 (10) h(xo) + Vh(xo )~ = 0 (11) where x0 represents the current guess or estimate as to the actual minimum value of the objective function, and H(x~ represents the Hessian at xo.
Here, as previously explained, the obj ective function is approximated quadratically wlule the constraints are linearly approximated. The minimum found for this approximate problem would be exact if the Hessian, (H(xo)), is also exact.
However, an inexact Hessian can be used in the foregoing formulation to save computation cost. By applying the above quadratic approximation successively, the real minimum of the objective function is obtained at convergence.
The terms "optimize" and "optimizing" as used with respect to this invention are intended to mean to determine or determining, respectively, an apparent optimum value of the objective function. As will be appreciated by those skilled in the art, in certain circumstances a localized optimum value of the objective function may be determined during any calculation procedure which seeks to determine the true ("global") optimum value of the objective function. Accordingly, the terms "optimize" and "optimizing" are intended to include within their scope any calculation procedure which seeks to determine an enhanced or optimum value of the objective function. Any allocation of fluid flow rates and/or surface facility operating parameters which result from such calculation procedure, whether the global optimum or a localized optimum value of the objective function is actually determined, are therefore also within the scope of this invention. In some instances, as will be readily appreciated by those skilled in the art, it may be desirable for a production system operator to intentionally select a fluid flow rate allocation among the wellbores that is less than optimal as determined by the optimizer. Accordingly, the invention shall not be limited in scope only to determining an optimal fluid flow rate allocation as a result of operating an optimization program according to the various embodiments of the invention.
In a particular embodiment of the invention, the Lagrange multipliers defined in equation (4) can be used to determine a sensitivity of the optimizer to any or all of the optimizer constraints. The values of one or more of the Lagrange multipliers are a measure of the sensitivity of the objective function to the associated constraints. The measure of sensitivity can be used to determine which of the constraints may be relaxed or otherwise adjusted to provide a substantial increase in the value of the system performance parameter that is to be optimized. As an example, a selected maximum total system water production may be a "bottleneck" to total oil production.
During optimization, the Lagrange multiplier associated with the maximum total system water production may indicate that a slight relaxation or adjustment of the selected maximum water production rate may provide the production system with the capacity to substantially increase maximum oil production rate, and correspondingly, the economic value (for example, net present value) of the production system.
The foregoing is meant to serve only as one example of use of the Lagrange multipliers calculated by the optimizer to determine constraint sensitivity. Any other constraint used in the optimizer may also undergo similar sensitivity analysis to determine production system "bottlenecks".
In one embodiment of a method according to the invention, a so-called "infeasible path" strategy is used, where the initial estimate or guess (x0) is allowed to be infeasible. "Infeasible" means that some or all of the constraints and variables are out of their respective minimum or maximum bounds. For example, one or more of the wellbores W may produce water at a rate which exceeds a maximum water production rate target for the entire system, or the total gas production, as another example, may exceed the capacity of the compressors. The optimization algorithm simultaneously tries to reach to an optimum as well as a feasible solution.
Thus feasibility is determined only at convergence. The advantage of this strategy is reduced objective and constraint function evaluation cost. How the infeasible solution strategy of the method of the invention is used will be further explained.
The solution of the optimization problem provides an optimal fluid flow rate and pressure distribution within the entire surface facility network. A part of this solution is then used in the reservoir simulator as the boundary conditions, while then solving the mass and momentum balance equations that describe the fluid flow in the reservoir.
A flow chart of how an optimization method according to the invention can be used in operating a production system is shown in Figure 2. After surface facility equations and reservoir equations are set up, and initial conditions in the surface facility and reservoir are set, at 40 the system time is incremented. If any surface facility operating parameters or structures have been changed from the previous calculation, shown at 42, such changes are entered into the conditions and/or equations for the surface facilities and reservoir. At 44, the conditions and constraints are entered into an optimization routine as previously described. At 46, the optimizer it is determined as to whether the optimizer has reached convergence. As previously explained, when the optimizer reaches convergence, an optimal value of the objective function is determined. When the optimal value of the objective function is determined, the system performance parameter which is represented by the objective function is at an optimal value. As previously explained, the performance parameter can be, for example, economic value, maximum oil production, minimized gas and/or water production, minimum operating cost, or any other parameter related to a measure of production and/or economic performance of the production system such as shown in Figure 1. The result of the optimization is an allocation of fluid production rates from each of the wellbores (W in Figure 1) which results in the optimization of the selected system performance parameters.
Referring again to Figure 2, the output of the optimizer includes fluid production rate allocation among the wellbores in the production system. In actual production and/or injection at the rates allocated by the optimizer, each wellbore (W
in Figure 1) will cause a pressure sink or pressure increase (depending on whether the wellbore is a producing well or injection well) at the reservoir. Such pressure changes propagate through the reservoir, and these pressure changes can be calculated using the mass and momentum balance equations referred to earlier. Therefore, as fluids are produced or injected into each wellbore W, a distribution of conditions in the subsurface reservoir changes. Using the output of the optimizer, the set of fluid flow rates for each wellbore as a set of boundary conditions, as shown at 62, a new distribution of conditions (particularly including but not limited to pressure) for the subsurface reservoir is calculated, at 64.
In some instances, the changes in reservoir conditions will result in changes in fluid flow rates from one or more of the wellbores (W in Figure 1). As these changes take place, they become part of the initial conditions for operating the optimizer, as indicated in Figure 2 by a line leading back to box 40.
In other cases, the optimizer will not converge. Failure of convergence, as explained earlier with reference to the description of the SQP aspect of the optimizer, is typically because at least one of the constraints is violated. The constraints may include operating parameters such as maximum acceptable water production in the system, maximum GOR, minimum inlet pressure to the compressor (14 in Figure 1), and others. In the event no system fluid production allocation will enable meeting all the constraints, the optimizer will not converge. In another aspect of the invention, a cause of the optimizer failing to converge may lead to isolation of one or more elements of the production system which cause the constraints to be violated.
At box 48 in Figure 2, one or more of the constraints may be relaxed or removed. For example a maximum acceptable water production may be increased, or removed as a constraint, or, alternatively, a minimum oil production may be reduced or removed as a constraint. Then, at box 50, the optimizer is run again. If convergence is achieved, then the violated constraint has been identified, at 52. At 54, corrective action can be taken to repair or correct the violated constraint. For example, if a maximum horsepower rating of the compressor (14 in Figure 1) is exceeded by a selected system gas flow rate, the compressor may be substituted by a higher rating compressor, and the optimizer run again, at 56. Any other physical change to the production system which alters or adjusts a system constraint can be detected and corrected by the method elements outlined in boxes 48, 50, 52 and 54, and the examples referred to herein should not be interpreted as limiting the types of system constraints that can be affected by the method of this invention. At box 58, if the optimizer has converged, then the flow rates axe allocated among the wellbores (W in Figure 1) according to the solution determined by the optimizer. At 60, these fluid flow rates are used as boundary conditions to perform a recalculation of the reservoir conditions, as in the earlier case where the initial run of the optimizer converged (at box 46).
While the invention has been described with respect to a limited number of embodiments, those skilled in the art, having benefit of this disclosure, will appreciate that other embodiments can be devised which do not depart from the scope of the invention as disclosed herein. Accordingly, the scope of the invention should be limited only by the attached claims.
Claims (71)
1. A method for enhancing allocation of fluid flow rates among a plurality of wellbores coupled to surface facilities, comprising:
modeling fluid flow characteristics of the wellbores and at least one reservoir penetrated thereby;
modeling fluid flow characteristics of the surface facilities;
operating an optimizer adapted to determine an enhanced value of an objective function, the objective function corresponding simultaneously to the modeled fluid flow characteristics of the wellbores and the surface facilities, the objective function relating to at least one production system performance parameter; and allocating fluid flow rates among the plurality of wellbores as determined by the operating the optimizer.
modeling fluid flow characteristics of the wellbores and at least one reservoir penetrated thereby;
modeling fluid flow characteristics of the surface facilities;
operating an optimizer adapted to determine an enhanced value of an objective function, the objective function corresponding simultaneously to the modeled fluid flow characteristics of the wellbores and the surface facilities, the objective function relating to at least one production system performance parameter; and allocating fluid flow rates among the plurality of wellbores as determined by the operating the optimizer.
2. The method as defined in claim 1 wherein the at least one production system performance parameter comprises economic value.
3. The method as defined in claim 1 wherein the at least one production system performance parameter comprises minimum water production rate.
4. The method as defined in claim 1 wherein the at least one production system performance parameter comprises minimum gas/oil ratio.
5. The method as defined in claim 1 wherein the at least one production system performance parameter comprises maximum oil production rate.
6. The method as defined in claim 1 wherein the at least one production system performance parameter comprises maximum ultimate recovery.
7. The method as defined in claim 1 wherein the objective function is optimized by successive quadratic programming.
8. The method as defined in claim 1 further comprising:
determining non-convergence of the objective function;
adjusting at least one constraint on the objective function;
recalculating the objective function; and repeating the adjusting at least one constraint and recalculating until the objective function converges.
determining non-convergence of the objective function;
adjusting at least one constraint on the objective function;
recalculating the objective function; and repeating the adjusting at least one constraint and recalculating until the objective function converges.
9. The method as defined in claim 8 further comprising:
repeating determining non-convergence of the objective function;
adjusting at least one element of the surface facilities;
recalculating the objective function;
repeating the adjusting at least one element and recalculating the objective function until the objective function converges.
repeating determining non-convergence of the objective function;
adjusting at least one element of the surface facilities;
recalculating the objective function;
repeating the adjusting at least one element and recalculating the objective function until the objective function converges.
10. The method as defined in claim 8 wherein the at least one constraint comprises maximum water production rate.
11. The method as defined in claim 8 wherein the at least one constraint comprises maximum gas/oil ratio.
12. The method as defined in claim 8 wherein the at least one constraint comprises maximum water cut.
13. The method as defined in claim 1 further comprising:
calculating a fluid pressure distribution in the at least one reservoir after a selected time interval;
recalculating fluid flow rates from the wellbores in response to the fluid pressure distribution calculation; and repeating the operating the optimizer and reallocating fluid flow rates among the wellbores in response to the repeated operating the optimizer.
calculating a fluid pressure distribution in the at least one reservoir after a selected time interval;
recalculating fluid flow rates from the wellbores in response to the fluid pressure distribution calculation; and repeating the operating the optimizer and reallocating fluid flow rates among the wellbores in response to the repeated operating the optimizer.
14. The method as defined in claim 1, further comprising:
determining a sensitivity of the objective function to at least one system constraint;
adjusting the at least one constraint and recalculating the objective function using the adjusted constraint; and reallocating fluid flow rates among the plurality of wellbores as determined by the recalculated objective function.
determining a sensitivity of the objective function to at least one system constraint;
adjusting the at least one constraint and recalculating the objective function using the adjusted constraint; and reallocating fluid flow rates among the plurality of wellbores as determined by the recalculated objective function.
15. The method as defined in claim 14 wherein determining the sensitivity comprises determining an optimal value of the objective function by sequential quadratic approximating, and determining a value of a Lagrange multiplier associated with the at least one constraint.
16. The method as defined in claim 1 wherein the optimizer comprises at least one constraint corresponding to a target value of at least one system parameter, the optimizer adapted to converge when a value of the at least one constraint is within a range bounded by the target value.
17. The method as defined in claim 16 wherein the at least one system parameter comprises a minimum oil production rate.
18. The method as defined in claim 16 wherein the at least one system parameter comprises a maximum water production rate.
19. A method for enhancing allocation of fluid flow rates among a plurality of wellbores coupled to surface facilities, comprising:
modeling fluid flow characteristics of the wellbores and at least one reservoir penetrated thereby;
modeling fluid flow characteristics of the surface facilities;
operating an optimizer adapted to determine an optimal value of an objective function, the objective function corresponding to the modeled fluid flow characteristics of the wellbores and the surface facilities, the objective function relating to at least one production system performance parameter, the optimizing comprising at least one constraint corresponding to a target value of at least one system operating parameter, the optimizer adapted to converge when a value of the at least one constraint is within a range bounded by the target value; and allocating fluid flow rates among the plurality of wellbores as determined by the operating the optimizer.
modeling fluid flow characteristics of the wellbores and at least one reservoir penetrated thereby;
modeling fluid flow characteristics of the surface facilities;
operating an optimizer adapted to determine an optimal value of an objective function, the objective function corresponding to the modeled fluid flow characteristics of the wellbores and the surface facilities, the objective function relating to at least one production system performance parameter, the optimizing comprising at least one constraint corresponding to a target value of at least one system operating parameter, the optimizer adapted to converge when a value of the at least one constraint is within a range bounded by the target value; and allocating fluid flow rates among the plurality of wellbores as determined by the operating the optimizer.
20. The method as defined in claim 19 wherein the at least one production system performance parameter comprises economic value.
21. The method as defined in claim 19 wherein the at least one production system performance parameter comprises water production rate.
22. The method as defined in claim 19 wherein the at least one production system performance parameter comprises minimum gas/oil ratio.
23. The method as defined in claim 19 wherein the at least one production system performance parameter comprises oil production rate.
24. The method as defined in claim 19 wherein the at least one production system performance parameter comprises ultimate recovery.
25. The method as defined in claim 19 wherein the optimizer comprises successive quadratic programming.
26. The method as defined in claim 19 further comprising:
determining non-convergence of the objective function;
adjusting the value of the at least one constraint;
recalculating the objective function; and repeating the adjusting the value of the at least one constraint and recalculating until the objective function converges.
determining non-convergence of the objective function;
adjusting the value of the at least one constraint;
recalculating the objective function; and repeating the adjusting the value of the at least one constraint and recalculating until the objective function converges.
27. The method as defined in claim 26 further comprising:
repeating determining non-convergence of the objective function;
adjusting at least one element of the surface facilities;
recalculating the objective function;
repeating the adjusting at least one element and recalculating until the objective function converges.
repeating determining non-convergence of the objective function;
adjusting at least one element of the surface facilities;
recalculating the objective function;
repeating the adjusting at least one element and recalculating until the objective function converges.
28. The method as defined in claim 26 wherein the at least one constraint comprises a maximum water production.
29. The method as defined in claim 26 wherein the at least one constraint comprises a maximum gas/oil ratio.
30. The method as defined in claim 26 wherein the at least one constraint comprises a maximum water cut.
31. The method as defined in claim 19 further comprising:
calculating a fluid pressure distribution in the at least one reservoir after a selected time interval;
recalculating fluid flow rates from the wellbores in response to the fluid pressure distribution calculation;
repeating the operating the optimizer; and reallocating fluid flow among the plurality of wellbores in response to the repeated operation of the optimizer.
calculating a fluid pressure distribution in the at least one reservoir after a selected time interval;
recalculating fluid flow rates from the wellbores in response to the fluid pressure distribution calculation;
repeating the operating the optimizer; and reallocating fluid flow among the plurality of wellbores in response to the repeated operation of the optimizer.
32. The method as defined in claim 19, further comprising:
determining a sensitivity of the objective function to at least one system operating constraint in a plurality of system operating constraints;
adjusting the at least one system operating constraint and recalculating the objective function using the adjusted system operating constraint; and reallocating fluid flow rates among the plurality of wellbores as determined by the recalculated objective function.
-
determining a sensitivity of the objective function to at least one system operating constraint in a plurality of system operating constraints;
adjusting the at least one system operating constraint and recalculating the objective function using the adjusted system operating constraint; and reallocating fluid flow rates among the plurality of wellbores as determined by the recalculated objective function.
-
33. The method as defined in claim 32 wherein determining the sensitivity comprises calculating the objective function by sequential quadratic approximating, and determining a value of a Lagrange multiplier associated with the at least one system operating constraint.
34. The method as defined in claim 32 wherein the at least one system operating constraint comprises a maximum water production.
35. The method as defined in claim 32 wherein the at least one system operating constraint comprises a maximum gas/oil ratio.
36. The method as defined in claim 32 wherein the at least one system operating constraint comprises a maximum water cut.
37. A method for optimizing allocation of fluid flow rates among a plurality of wellbores coupled to surface facilities, comprising:
modeling fluid flow characteristics of the wellbores and at least one reservoir penetrated thereby;
modeling fluid flow characteristics of the surface facilities;
optimizing an objective function, the objective function corresponding simultaneously to the modeled fluid flow characteristics of the wellbores and the surface facilities, the objective function relating to at least one production system performance parameter; and allocating fluid flow rates among the plurality of wellbores as determined by the optimizing.
modeling fluid flow characteristics of the wellbores and at least one reservoir penetrated thereby;
modeling fluid flow characteristics of the surface facilities;
optimizing an objective function, the objective function corresponding simultaneously to the modeled fluid flow characteristics of the wellbores and the surface facilities, the objective function relating to at least one production system performance parameter; and allocating fluid flow rates among the plurality of wellbores as determined by the optimizing.
38. The method as defined in claim 37 wherein the at least one production system performance parameter comprises economic value.
39. The method as defined in claim 37 wherein the at least one production system performance parameter comprises water production rate.
40. The method as defined in claim 37 wherein the at least one production system performance parameter comprises gas/oil ratio.
41. The method as defined in claim 37 wherein the at least one production system performance parameter comprises oil production rate.
42. The method as defined in claim 37 wherein the at least one production system performance parameter comprises ultimate recovery.
43. The method as defined in claim 37 wherein the objective function is optimized by successive quadratic programming.
44. The method as defined in claim 37 further comprising:
determining non-convergence of the objective function;
adjusting at least one constraint on the objective function;
recalculating the objective function; and repeating the adjusting at least one constraint and recalculating until the objective function converges.
determining non-convergence of the objective function;
adjusting at least one constraint on the objective function;
recalculating the objective function; and repeating the adjusting at least one constraint and recalculating until the objective function converges.
45. The method as defined in claim 44 further comprising:
repeating determining non-convergence of the objective function;
adjusting at least one element of the surface facilities;
recalculating the objective function;
repeating the adjusting at least one element and recalculating until the objective function converges.
repeating determining non-convergence of the objective function;
adjusting at least one element of the surface facilities;
recalculating the objective function;
repeating the adjusting at least one element and recalculating until the objective function converges.
46. The method as defined in claim 44 wherein the at least one constraint comprises water production rate.
47. The method as defined in claim 44 wherein the at least one constraint comprises gas/oil ratio.
48. The method as defined in claim 44 wherein the at least one constraint comprises water cut.
49. The method as defined in claim 37 further comprising:
calculating a fluid pressure distribution in the at least one reservoir after a selected time interval;
recalculating fluid flow rates from the wellbores in response to the fluid pressure distribution calculation;
repeating the optimizing the objective function; and reallocating fluid flow among the plurality of wellbores in response to the repeated optimizing.
calculating a fluid pressure distribution in the at least one reservoir after a selected time interval;
recalculating fluid flow rates from the wellbores in response to the fluid pressure distribution calculation;
repeating the optimizing the objective function; and reallocating fluid flow among the plurality of wellbores in response to the repeated optimizing.
50. The method as defined in claim 37, further comprising:
determining a sensitivity of the objective function to at least one system constraint;
adjusting the at least one constraint and recalculating the objective function using the adjusted constraint; and reallocating fluid flow rates among the plurality of wellbores as determined by the recalculated objective function.
determining a sensitivity of the objective function to at least one system constraint;
adjusting the at least one constraint and recalculating the objective function using the adjusted constraint; and reallocating fluid flow rates among the plurality of wellbores as determined by the recalculated objective function.
51. The method as defined in claim 50 wherein determining the sensitivity comprises optimizing the objective function by sequential quadratic approximating, and determining a value of a Lagrange multiplier associated with the at least one constraint.
52. The method as defined in claim 37 wherein the optimizing comprises at least one constraint corresponding to a target value of at least one system parameter, the optimizing adapted to converge when a value of the at least one constraint is within a range bounded by the target value.
53. The method as defined in claim 52 wherein the at least one system parameter comprises a minimum oil production rate.
54. The method as defined in claim 52 wherein the at least one system parameter comprises a maximum water production rate.
55. A method for optimizing allocation of fluid flow among a plurality of wellbores coupled to surface facilities, comprising:
modeling fluid flow characteristics of the wellbores and at least one reservoir penetrated thereby;
modeling fluid flow characteristics of the surface facilities;
optimizing an objective function, the objective function corresponding to the modeled fluid flow characteristics of the wellbores and the surface facilities, the objective function relating to at least one production system performance parameter;
determining a sensitivity of the objective function to at least one system constraint;
adjusting the at least one system constraint and recalculating the objective function using the adjusted system constraint; and reallocating fluid flow rates among the plurality of wellbores as determined by the recalculated objective function; and allocating fluid flow rates among the plurality of wellbores as determined by the optimizing.
modeling fluid flow characteristics of the wellbores and at least one reservoir penetrated thereby;
modeling fluid flow characteristics of the surface facilities;
optimizing an objective function, the objective function corresponding to the modeled fluid flow characteristics of the wellbores and the surface facilities, the objective function relating to at least one production system performance parameter;
determining a sensitivity of the objective function to at least one system constraint;
adjusting the at least one system constraint and recalculating the objective function using the adjusted system constraint; and reallocating fluid flow rates among the plurality of wellbores as determined by the recalculated objective function; and allocating fluid flow rates among the plurality of wellbores as determined by the optimizing.
56. The method as defined in claim 55 wherein the at least one production system performance parameter comprises economic value.
57. The method as defined in claim 55 wherein the at least one production system performance parameter comprises water production rate.
58. The method as defined in claim 55 wherein the at least one production system performance parameter comprises gas/oil ratio.
59. The method as defined in claim 55 wherein the at least one production system performance parameter comprises oil production rate.
60. The method as defined in claim 55 wherein the at least one production system performance parameter comprises ultimate recovery.
61. The method as defined in claim 55 wherein the objective function is optimized by successive quadratic programming.
62. The method as defined in claim 55 further comprising:
determining non-convergence of the objective function;
adjusting at least one constraint on the objective function;
recalculating the objective function; and repeating the adjusting at least one constraint and recalculating until the objective function converges.
determining non-convergence of the objective function;
adjusting at least one constraint on the objective function;
recalculating the objective function; and repeating the adjusting at least one constraint and recalculating until the objective function converges.
63. The method as defined in claim 62 further comprising:
repeating determining non-convergence of the objective function;
adjusting at least one element of the surface facilities;
recalculating the objective function;
repeating the adjusting at least one element and recalculating until the objective function converges.
repeating determining non-convergence of the objective function;
adjusting at least one element of the surface facilities;
recalculating the objective function;
repeating the adjusting at least one element and recalculating until the objective function converges.
64. The method as defined in claim 62 wherein the at least one constraint comprises water production rate.
65. The method as defined in claim 62 wherein the at least one constraint comprises gas/oil ratio.
66. The method as defined in claim 62 wherein the at least one constraint comprises water cut.
67. The method as defined in claim 55 further comprising:
calculating a fluid pressure distribution in the at least one reservoir after a selected time interval;
recalculating fluid flow rates from the wellbores in response to the fluid pressure distribution calculation;
repeating the optimizing the objective function; and reallocating fluid flow among the plurality of wellbores in response to the repeated optimizing.
calculating a fluid pressure distribution in the at least one reservoir after a selected time interval;
recalculating fluid flow rates from the wellbores in response to the fluid pressure distribution calculation;
repeating the optimizing the objective function; and reallocating fluid flow among the plurality of wellbores in response to the repeated optimizing.
68. The method as defined in claim 55 wherein determining the sensitivity comprises optimizing the objective function by sequential quadratic approximating, and determining a value of a Lagrange multiplier associated with the at least one constraint.
69. The method as defined in claim 55 wherein the optimizing comprises at least one constraint corresponding to a target value of at least one system parameter, the optimizing adapted to converge when a value of the at least one constraint corresponding to the target value is within a range bounded by the target value.
70. The method as defined in claim 69 wherein the at least one system parameter comprises an oil production rate.
71. The method as defined in claim 69 wherein the at least one system parameter comprises a water production rate.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US28613401P | 2001-04-24 | 2001-04-24 | |
| US60/286,134 | 2001-04-24 | ||
| PCT/US2002/012287 WO2002086277A2 (en) | 2001-04-24 | 2002-04-19 | Method for enhancing production allocation in an integrated reservoir and surface flow system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2442596A1 true CA2442596A1 (en) | 2002-10-31 |
Family
ID=23097232
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002442596A Abandoned CA2442596A1 (en) | 2001-04-24 | 2002-04-19 | Method for enhancing production allocation in an integrated reservoir and surface flow system |
Country Status (7)
| Country | Link |
|---|---|
| US (2) | US7379853B2 (en) |
| EP (1) | EP1389259B1 (en) |
| AT (1) | ATE310890T1 (en) |
| CA (1) | CA2442596A1 (en) |
| DE (1) | DE60207549D1 (en) |
| NO (1) | NO20034745L (en) |
| WO (1) | WO2002086277A2 (en) |
Families Citing this family (84)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EA005604B1 (en) | 2001-02-05 | 2005-04-28 | Шлумбергер Холдингс Лимитид | Optimization of reservoir, well and surface network systems |
| US8401832B2 (en) * | 2002-11-23 | 2013-03-19 | Schlumberger Technology Corporation | Method and system for integrated reservoir and surface facility networks simulations |
| US7096172B2 (en) * | 2003-01-31 | 2006-08-22 | Landmark Graphics Corporation, A Division Of Halliburton Energy Services, Inc. | System and method for automated reservoir targeting |
| US7200540B2 (en) * | 2003-01-31 | 2007-04-03 | Landmark Graphics Corporation | System and method for automated platform generation |
| WO2006065270A1 (en) * | 2004-05-11 | 2006-06-22 | Quantm Ltd. | Path analysis system |
| US7627461B2 (en) * | 2004-05-25 | 2009-12-01 | Chevron U.S.A. Inc. | Method for field scale production optimization by enhancing the allocation of well flow rates |
| CA2569102A1 (en) * | 2004-06-07 | 2005-12-22 | Exxonmobil Upstream Research Company | Method for solving implicit reservoir simulation matrix equation |
| CA2570058C (en) | 2004-06-25 | 2013-07-30 | Shell Canada Limited | Closed loop control system for controlling production of hydrocarbon fluid from an underground formation |
| US7809537B2 (en) * | 2004-10-15 | 2010-10-05 | Saudi Arabian Oil Company | Generalized well management in parallel reservoir simulation |
| US7373285B2 (en) * | 2004-12-01 | 2008-05-13 | Bp Corporation North America Inc. | Application of phase behavior models in production allocation systems |
| US7596480B2 (en) * | 2005-04-14 | 2009-09-29 | Saudi Arabian Oil Company | Solution method and apparatus for large-scale simulation of layered formations |
| AU2006344398B2 (en) * | 2005-10-06 | 2011-05-19 | Logined B.V. | Method, system and apparatus for numerical black oil delumping |
| EA014140B1 (en) * | 2005-11-21 | 2010-10-29 | Шеврон Ю.Эс.Эй. Инк. | Method for field scale production optimization |
| US8195401B2 (en) | 2006-01-20 | 2012-06-05 | Landmark Graphics Corporation | Dynamic production system management |
| US8571688B2 (en) * | 2006-05-25 | 2013-10-29 | Honeywell International Inc. | System and method for optimization of gas lift rates on multiple wells |
| US10260329B2 (en) * | 2006-05-25 | 2019-04-16 | Honeywell International Inc. | System and method for multivariable control in three-phase separation oil and gas production |
| GB2456925B (en) * | 2006-10-30 | 2011-08-10 | Logined Bv | System and method for performing oilfield simulation operations |
| BRPI0720188B1 (en) * | 2006-10-31 | 2018-10-16 | Exxonmobil Upstream Res Co | modeling method of a computer-readable storage system and reservoir system |
| US8078444B2 (en) * | 2006-12-07 | 2011-12-13 | Schlumberger Technology Corporation | Method for performing oilfield production operations |
| US8190458B2 (en) * | 2007-01-17 | 2012-05-29 | Schlumberger Technology Corporation | Method of performing integrated oilfield operations |
| US7946127B2 (en) * | 2007-02-21 | 2011-05-24 | Honeywell International Inc. | Apparatus and method for optimizing a liquefied natural gas facility |
| US8775141B2 (en) * | 2007-07-02 | 2014-07-08 | Schlumberger Technology Corporation | System and method for performing oilfield simulation operations |
| US20090076632A1 (en) * | 2007-09-18 | 2009-03-19 | Groundswell Technologies, Inc. | Integrated resource monitoring system with interactive logic control |
| US8892221B2 (en) * | 2007-09-18 | 2014-11-18 | Groundswell Technologies, Inc. | Integrated resource monitoring system with interactive logic control for well water extraction |
| US8140310B2 (en) | 2007-11-01 | 2012-03-20 | Schlumberger Technology Corporation | Reservoir fracture simulation |
| US7668707B2 (en) * | 2007-11-28 | 2010-02-23 | Landmark Graphics Corporation | Systems and methods for the determination of active constraints in a network using slack variables and plurality of slack variable multipliers |
| EP2247820A4 (en) * | 2007-12-13 | 2016-02-24 | Exxonmobil Upstream Res Co | Parallel adaptive data partitioning on a reservoir simulation using an unstructured grid |
| EP2232301A1 (en) * | 2007-12-21 | 2010-09-29 | ExxonMobil Upstream Research Company | Modeling in sedimentary basins |
| US20090198505A1 (en) * | 2008-02-05 | 2009-08-06 | Peter Gipps | Interactive path planning with dynamic costing |
| EP2151540A1 (en) * | 2008-06-16 | 2010-02-10 | Bp Exploration Operating Company Limited | Method and apparatus for configuring oil and/or gas producing system |
| EP2161406A1 (en) * | 2008-09-03 | 2010-03-10 | BP Exploration Operating Company Limited | Method and apparatus for configuring oil and/or gas producing system |
| US20110087476A1 (en) * | 2008-06-16 | 2011-04-14 | Patrick James Calvert | Method and apparatus for configuring oil and/or gas production system |
| US8670966B2 (en) * | 2008-08-04 | 2014-03-11 | Schlumberger Technology Corporation | Methods and systems for performing oilfield production operations |
| EP2350915A4 (en) * | 2008-09-30 | 2013-06-05 | Exxonmobil Upstream Res Co | Method for solving reservoir simulation matrix equation using parallel multi-level incomplete factorizations |
| BRPI0919456A2 (en) * | 2008-09-30 | 2015-12-22 | Exxonmobil Upstream Res Co | method for modeling fluid flow in a hydrocarbon reservoir |
| US9228415B2 (en) | 2008-10-06 | 2016-01-05 | Schlumberger Technology Corporation | Multidimensional data repository for modeling oilfield operations |
| US8301426B2 (en) * | 2008-11-17 | 2012-10-30 | Landmark Graphics Corporation | Systems and methods for dynamically developing wellbore plans with a reservoir simulator |
| AU2010208105B2 (en) * | 2009-01-30 | 2015-01-22 | Chevron U.S.A. Inc. | System and method for predicting fluid flow in subterranean reservoirs |
| US8490696B2 (en) | 2009-04-20 | 2013-07-23 | David Randolph Smith | Method and apparatus to enhance oil recovery in wells |
| AU2014201895B2 (en) * | 2009-04-20 | 2015-09-03 | David Randolph Smith | Method and apparatus to enhance oil recovery in wells |
| US20120130696A1 (en) * | 2009-08-12 | 2012-05-24 | Exxonmobil Upstream Research Company | Optimizing Well Management Policy |
| US9134454B2 (en) | 2010-04-30 | 2015-09-15 | Exxonmobil Upstream Research Company | Method and system for finite volume simulation of flow |
| US8463586B2 (en) | 2010-06-22 | 2013-06-11 | Saudi Arabian Oil Company | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids |
| AU2011283190A1 (en) | 2010-07-29 | 2013-02-07 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
| EP2599023B1 (en) | 2010-07-29 | 2019-10-23 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
| EP2599032A4 (en) | 2010-07-29 | 2018-01-17 | Exxonmobil Upstream Research Company | Method and system for reservoir modeling |
| AU2011283192B2 (en) | 2010-07-29 | 2014-07-17 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
| US8433551B2 (en) | 2010-11-29 | 2013-04-30 | Saudi Arabian Oil Company | Machine, computer program product and method to carry out parallel reservoir simulation |
| US8386227B2 (en) | 2010-09-07 | 2013-02-26 | Saudi Arabian Oil Company | Machine, computer program product and method to generate unstructured grids and carry out parallel reservoir simulation |
| US9058446B2 (en) | 2010-09-20 | 2015-06-16 | Exxonmobil Upstream Research Company | Flexible and adaptive formulations for complex reservoir simulations |
| WO2012082273A1 (en) | 2010-12-13 | 2012-06-21 | Chevron U.S.A. Inc. | Method and system for coupling reservoir and surface facility simulations |
| US8972232B2 (en) * | 2011-02-17 | 2015-03-03 | Chevron U.S.A. Inc. | System and method for modeling a subterranean reservoir |
| US20120215364A1 (en) * | 2011-02-18 | 2012-08-23 | David John Rossi | Field lift optimization using distributed intelligence and single-variable slope control |
| WO2013039606A1 (en) | 2011-09-15 | 2013-03-21 | Exxonmobil Upstream Research Company | Optimized matrix and vector operations in instruction limited algorithms that perform eos calculations |
| US20130110474A1 (en) * | 2011-10-26 | 2013-05-02 | Nansen G. Saleri | Determining and considering a premium related to petroleum reserves and production characteristics when valuing petroleum production capital projects |
| US9767421B2 (en) * | 2011-10-26 | 2017-09-19 | QRI Group, LLC | Determining and considering petroleum reservoir reserves and production characteristics when valuing petroleum production capital projects |
| US9946986B1 (en) | 2011-10-26 | 2018-04-17 | QRI Group, LLC | Petroleum reservoir operation using geotechnical analysis |
| US20140326449A1 (en) * | 2012-02-24 | 2014-11-06 | Landmark Graphics Corporation | Determining optimal parameters for a downhole operation |
| US9803457B2 (en) | 2012-03-08 | 2017-10-31 | Schlumberger Technology Corporation | System and method for delivering treatment fluid |
| US9863228B2 (en) * | 2012-03-08 | 2018-01-09 | Schlumberger Technology Corporation | System and method for delivering treatment fluid |
| AU2013274731B2 (en) * | 2012-06-15 | 2016-08-25 | Landmark Graphics Corporation | Systems and methods for optimizing facility limited production and injection in an integrated reservoir and gathering network |
| US9031822B2 (en) | 2012-06-15 | 2015-05-12 | Chevron U.S.A. Inc. | System and method for use in simulating a subterranean reservoir |
| CA2874994C (en) * | 2012-06-15 | 2017-02-07 | Landmark Graphics Corporation | Systems and methods for solving a multi-reservoir system with heterogeneous fluids coupled to a common gathering network |
| EP2817734B1 (en) * | 2012-06-15 | 2018-05-23 | Landmark Graphics Corporation | Methods and systems for gas lift rate management |
| US10036829B2 (en) | 2012-09-28 | 2018-07-31 | Exxonmobil Upstream Research Company | Fault removal in geological models |
| US10012071B2 (en) | 2013-07-11 | 2018-07-03 | Laurie Sibbald | Differential method for equitable allocation of hydrocarbon component yields using phase behavior process models |
| US20150337631A1 (en) * | 2014-05-23 | 2015-11-26 | QRI Group, LLC | Integrated production simulator based on capacitance-resistance model |
| CA2948667A1 (en) | 2014-07-30 | 2016-02-04 | Exxonmobil Upstream Research Company | Method for volumetric grid generation in a domain with heterogeneous material properties |
| US10443358B2 (en) | 2014-08-22 | 2019-10-15 | Schlumberger Technology Corporation | Oilfield-wide production optimization |
| US9951601B2 (en) | 2014-08-22 | 2018-04-24 | Schlumberger Technology Corporation | Distributed real-time processing for gas lift optimization |
| AU2015339883B2 (en) | 2014-10-31 | 2018-03-29 | Exxonmobil Upstream Research Company | Methods to handle discontinuity in constructing design space for faulted subsurface model using moving least squares |
| WO2016069171A1 (en) | 2014-10-31 | 2016-05-06 | Exxonmobil Upstream Research Company | Handling domain discontinuity in a subsurface grid model with the help of grid optimization techniques |
| FR3054705B1 (en) * | 2016-07-29 | 2018-07-27 | Veolia Environnement-VE | TOOL FOR MANAGING MULTIPLE WATER RESOURCES |
| US11263370B2 (en) | 2016-08-25 | 2022-03-01 | Enverus, Inc. | Systems and methods for allocating hydrocarbon production values |
| US10303819B2 (en) | 2016-08-25 | 2019-05-28 | Drilling Info, Inc. | Systems and methods for allocating hydrocarbon production values |
| CA3043231C (en) | 2016-12-23 | 2022-06-14 | Exxonmobil Upstream Research Company | Method and system for stable and efficient reservoir simulation using stability proxies |
| CN108729911A (en) * | 2017-04-24 | 2018-11-02 | 通用电气公司 | Optimization devices, systems, and methods for resource production system |
| US20210165124A1 (en) * | 2017-06-14 | 2021-06-03 | Landmark Graphics Corporation | Modeling Geological Strata Using Weighted Parameters |
| US11466554B2 (en) | 2018-03-20 | 2022-10-11 | QRI Group, LLC | Data-driven methods and systems for improving oil and gas drilling and completion processes |
| US11506052B1 (en) | 2018-06-26 | 2022-11-22 | QRI Group, LLC | Framework and interface for assessing reservoir management competency |
| US11078773B2 (en) | 2018-12-03 | 2021-08-03 | Saudi Arabian Oil Company | Performing continuous daily production allocation |
| US11714393B2 (en) | 2019-07-12 | 2023-08-01 | Johnson Controls Tyco IP Holdings LLP | Building control system with load curtailment optimization |
| US11226600B2 (en) * | 2019-08-29 | 2022-01-18 | Johnson Controls Tyco IP Holdings LLP | Building control system with load curtailment optimization |
| CA3148100A1 (en) * | 2019-08-07 | 2021-02-11 | Drilling Info, Inc. | Systems and methods for allocating hydrocarbon production values |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5305209A (en) * | 1991-01-31 | 1994-04-19 | Amoco Corporation | Method for characterizing subterranean reservoirs |
| FR2742794B1 (en) * | 1995-12-22 | 1998-01-30 | Inst Francais Du Petrole | METHOD FOR MODELING THE EFFECTS OF WELL INTERACTIONS ON THE AQUEOUS FRACTION PRODUCED BY AN UNDERGROUND HYDROCARBON DEPOSIT |
| US6112126A (en) * | 1997-02-21 | 2000-08-29 | Baker Hughes Incorporated | Adaptive object-oriented optimization software system |
| US5992519A (en) * | 1997-09-29 | 1999-11-30 | Schlumberger Technology Corporation | Real time monitoring and control of downhole reservoirs |
| US6236894B1 (en) * | 1997-12-19 | 2001-05-22 | Atlantic Richfield Company | Petroleum production optimization utilizing adaptive network and genetic algorithm techniques |
| US6101447A (en) * | 1998-02-12 | 2000-08-08 | Schlumberger Technology Corporation | Oil and gas reservoir production analysis apparatus and method |
| US6266619B1 (en) * | 1999-07-20 | 2001-07-24 | Halliburton Energy Services, Inc. | System and method for real time reservoir management |
| US6928399B1 (en) | 1999-12-03 | 2005-08-09 | Exxonmobil Upstream Research Company | Method and program for simulating a physical system using object-oriented programming |
| US20020177955A1 (en) * | 2000-09-28 | 2002-11-28 | Younes Jalali | Completions architecture |
| EA005604B1 (en) | 2001-02-05 | 2005-04-28 | Шлумбергер Холдингс Лимитид | Optimization of reservoir, well and surface network systems |
| US7627461B2 (en) * | 2004-05-25 | 2009-12-01 | Chevron U.S.A. Inc. | Method for field scale production optimization by enhancing the allocation of well flow rates |
| EA014140B1 (en) | 2005-11-21 | 2010-10-29 | Шеврон Ю.Эс.Эй. Инк. | Method for field scale production optimization |
-
2002
- 2002-04-19 US US10/126,215 patent/US7379853B2/en active Active
- 2002-04-19 AT AT02728833T patent/ATE310890T1/en not_active IP Right Cessation
- 2002-04-19 WO PCT/US2002/012287 patent/WO2002086277A2/en not_active Application Discontinuation
- 2002-04-19 CA CA002442596A patent/CA2442596A1/en not_active Abandoned
- 2002-04-19 EP EP02728833A patent/EP1389259B1/en not_active Revoked
- 2002-04-19 DE DE60207549T patent/DE60207549D1/en not_active Expired - Lifetime
-
2003
- 2003-10-23 NO NO20034745A patent/NO20034745L/en unknown
-
2007
- 2007-10-31 US US11/981,411 patent/US7752023B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| NO20034745L (en) | 2003-12-23 |
| US20020165671A1 (en) | 2002-11-07 |
| EP1389259A2 (en) | 2004-02-18 |
| EP1389259A4 (en) | 2004-06-09 |
| US7379853B2 (en) | 2008-05-27 |
| WO2002086277A3 (en) | 2003-05-22 |
| ATE310890T1 (en) | 2005-12-15 |
| EP1389259B1 (en) | 2005-11-23 |
| US20080065363A1 (en) | 2008-03-13 |
| US7752023B2 (en) | 2010-07-06 |
| DE60207549D1 (en) | 2005-12-29 |
| WO2002086277A2 (en) | 2002-10-31 |
| NO20034745D0 (en) | 2003-10-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1389259B1 (en) | Method for enhancing production allocation in an integrated reservoir and surface flow system | |
| Kosmidis et al. | A mixed integer optimization formulation for the well scheduling problem on petroleum fields | |
| US20120215364A1 (en) | Field lift optimization using distributed intelligence and single-variable slope control | |
| US20120130696A1 (en) | Optimizing Well Management Policy | |
| US20180023373A1 (en) | Method and apparatus for configuring oil and/or gas producing system | |
| US20140094974A1 (en) | Lift and choke control | |
| WO2008055186A2 (en) | System and method for performing oilfield simulation operations | |
| Camponogara et al. | An automation system for gas-lifted oil wells: Model identification, control, and optimization | |
| EP3339565B1 (en) | Systems and methods for assessing production and/or injection system startup | |
| Wilson et al. | Geomechanical modeling of flowback scenarios to establish best practices in the Midland Basin horizontal program | |
| Ribeiro et al. | Model Predictive Control with quality requirements on petroleum production platforms | |
| EP2847708B1 (en) | Methods and systems for non-physical attribute management in reservoir simulation | |
| Lasdon et al. | Optimal hydrocarbon reservoir production policies | |
| Hohendorff Filho et al. | Effect of reservoir and production system integration on field production strategy selection | |
| AU2002258860A1 (en) | Method for enhancing production allocation in an integrated reservoir and surface flow system | |
| Plucenio et al. | Gas-lift optimization and control with nonlinear mpc | |
| Binder | Production optimization in a cluster of gas-lift wells | |
| US20100100354A1 (en) | Dynamic calculation of allocation factors for a producer well | |
| Khairy et al. | Gas Lift Optimization for a Sudanese Field | |
| EP2161406A1 (en) | Method and apparatus for configuring oil and/or gas producing system | |
| Storvold | Optimization of investment decisions and production planning in aging offshore petroleum fields | |
| US20230408043A1 (en) | Multiphase flow instability control | |
| Soares et al. | Using the Equal Slope Methods to Optimized Artificial Lift in K Field With Two Artificial Lift Types | |
| Georgiadis et al. | Integrated optimization of oil and gas production | |
| Goridko et al. | SPE-201878-MS |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| FZDE | Discontinued |
Effective date: 20160418 |