|Publication number||US3891836 A|
|Publication date||Jun 24, 1975|
|Filing date||Mar 25, 1974|
|Priority date||Apr 21, 1972|
|Publication number||US 3891836 A, US 3891836A, US-A-3891836, US3891836 A, US3891836A|
|Original Assignee||Mobil Oil Corp|
|Export Citation||BiBTeX, EndNote, RefMan|
|Patent Citations (5), Referenced by (51), Classifications (12)|
|External Links: USPTO, USPTO Assignment, Espacenet|
United States Patent Lee June 24, 1975 APPARATUS FOR OPIIMIZING I MULTIUNIT PROCESSING SYSTEMS m ry Examinerloseph F. Ruggiero  Inventor: Wooyoung Lee, Cherry Hill, NJ. Anomey' Agem or Flrm c' Huggeu  Assignee: all? Oil Corporation, New York.  ABSTRACT A processing system comprising a plurality of individl Filed: ual units is optimized by first estimating the yields of A L No; 454,620 the individual units at a standard set of operating con- I pp ditions and then establishing optimum flow rates using Related Application Data a linear programming model or similar mathematical I63] Continuation-in-part of Ser. No. 246,445, April 21, tehniques. Individual units are then controlled and lol972,abandoned. cally optimized consistent with a sensitivity analysis which is performed by treating a proposed change in  U-S- Cl. 3 5l-l2; 235/150 /1 operating conditions of an individual unit as a distur-  Int. Cl. G06f 15/46; 606g 7/58 bance in the unit yield column of the linear r0 ram- P 8  Field of Seard 235/1501. ming model. The overall system may then be opti- 235/l5l.l2; 444/1 mized for the changes in operating conditions by determining and establishing new flow rates. The steps of  References Cited changing operating conditions and establishing new UNlTED STATES PATENTS flow rates may be repeated until the sensitivity analy- 3 075 700 M963 Bishop H 235/150}l sis reveals that any further change in operating condi- 3:079:079 2/1963 p i et a] H 235/1501 tions will not further improve the profit of the overall 320L572 8/1965 Yetter 235/!5] y 3,594,559 7/1971 Pemberton 235M501 X 3,62l,2l7 11/1971 Carr 235/1501 x 18 C|a|m5- 19 D'awmg Figures OVERALL SYSTEM COMPUTER i817 2 T y-1a, r- 7" 7 -1 l l M I 1 I'm a f i LOCAL 2 l I LOAL LOCAL 12, come j I come COMP. UNIT l UNlT UNIT- 1 5 6 lS l 12;, "1| l I j l I 1212 I LOCAL I COMP. D3 l i umt '29 141 i I 3 l i I25 I213 l UNIT l UN" 2 10 I 7 i bg M ii 12 1 I UNIT 1 1 I I2 4 l2 I 4 t\ 5 I L f" J PATENTEI] JUN 2 4 I975 AM 5 mm W m 4 llllll... N M Wm U I, L I MWM W3 E 3 LC 44 Q a a 0 M f I 4 2 PATENTEDJUN24 I975 3.891, 836
SHEET 2 ESTIMATE UNIT COSTS AND YIELDS AT STANDARD OPERATING CONDITIONS Fig Z TRANSMIT COSTS AND YIELDS TO SYSTEM PLANNING /24 SOLVE LP FOR T TRANSMITTED YIELDS TRANSMIT LOCAL /26 OBJECTIVE FUNCTION AND LP SOLUTION TO LOCAL UNITS DETERMINE IF A CHANGE IN YIELD OF UNIT WILL CORRECT SATISFY LOCAL LP OBJECTIVE FUNCTION SOLUTION IF CHANGE SATISFIES LOCAL OBJECTIVE 30 FUNCTION, TRANSMIT TO SYSTEM PLANNING PATENTEDJUN24 I975 3.891.836
SHEET 3 ESTIMATE UNIT COSTS AND YIELDS AT STANDARD SYSTEM PLANNING i /24 SOLVE LP FOR J TRANSMITTED YIELDS TRANSMIT LOCAL 26 OBJECTIVE FUNCTION AND LP SOLUTION TO LOCAL UNITS T Ml FA HANGE I N YI IE LD g T T WILL TRANSMIT CORRECTED SATISFY LOCAL LP SOLUTION TO OBJECTIVE FUNCTION SYSTEM PLANNING IF CHANGE SATISFIES LOCAL OBJECTIVE FUNCTION, CORRECT my 30a LP SOLUTION SOLVE LP AT OVERALL SYSTEM COMPUTER USING ESTIMATED YIELDS CHANGE IN YIELDS WILL SATISFY LO OBJECTIVE FUNCTION DETERMINE IF IF CHANGE SATISFIES LOCAL OBJECTIVE FUNCTIONJRANSMIT CHANGE TO LOCAL UNIT DETERMINE IF LOCAL UNIT CAN MAKE CHANGE IF LOCAL UNIT CAN TO SYSTEM COMPUTER MAKE CHANGE, TRANSMIT CORRECT PATENTEDJuN24 I975 3. 891. 8 36 SHEET 5 OVERALL 54 SYSTEM COMPUTER I 'lg- 4 A 56 56 x! CA /48 i I GASOLINE CRuOE 1 --l FUEL CRUDE 2'---"""" PROCESSNG HEATING OIL CRUDE 3/i/v JET FUEL 3 L 52 CRuOE 4 L 50 LOCAL COMP.
LUBE I PROCESSING l I LUBE OIL Fig 5 A 24 CONSTRAINTMI) 2 x CONSTRAINTMZ) PATENTEI] JUN 2 4 1915 SHEEI SHEET PATENTEDJIJN 24 I975 NEW 300 SIMPLEX MQLTIPLI ERS NEW COSTS NEW BASIS MATRIX liq- I982 202 17 /2002 LT@Q1V SHEET PATENTED JUN 24 I975 lsp SHEET 2|4 CHECK IF c; 20
PATENTEDJUN 24 ms SW TCH NG SHEET PATENTED JUN 2 4 I975 W on w m m m 2 4 a w w w w a m a ll lllllllllllll ll c 8 PATENTEDJUN 24 ms sum 14 13 8 91; 8 36 o o o o c STORE l 11 1 soo STORE Ap 502 COMPUTE SENSITIVITY COEFFICIENTS 5' o 504 5= 1 6 A S1 Em COMPUTE NEW X's o sos J 4 1+ -m COMPUTE NEW INVERSE BASIS MATRIX 5, 50s I o I Q Q g m 0 I1 S -m COMPUTE NEW SIMPLEX MULTIPLIERS 51o 'm 2K S C COMPUTE NEW COST COEFFICIENTS 1 APPARATUS FOR OPTIMIZING MULTIUNIT PROCESSING SYSTEMS RELATED APPLICATION This is a continuation-in-part of copending application Ser. No. 246,445, filed Apr. 21, 1972, now abandoned which is incorporated herein by reference.
BACKGROUND OF THE INVENTION This invention relates to a method for optimizing large, complex processing systems comprising a plurality of individual units.
A number of techniques are available for optimization of such large processing systems. Linear programming is one of the most widely used techniques, and modern refinery complexes are optimized by linear programming almost without exception. However, a linear programming model is at best an approximation of a real physical system. For example, it is commonly assumed that the yield information incorporated into the LP (linear programming model) is relatively fixed and the coefficients of the constraint equations which represent these yields are also fixed. In reality, however, these yield coefficients are dependent upon the operating conditions of the units and any change in the operating conditions for a single unit will of necessity affect the operation of other units.
In many instances, the operating conditions may be closely controlled. Local unit managers have appropriate tools (i.e., mathematical models, optimizers, and process control computers) for the local optimization and control of their individual units. As a result, large amounts of detailed information about the processes are continuously generated which are valuable for accurate adjustment of operating conditions. Each process unit can thus be locally optimized and controlled continuously.
However, any change in operating conditions of a unit which locally optimizes the unit will not necessarily optimize the overall system. In fact, a change in operating conditions which optimizes the local unit may have an adverse effect on the overall system. As a result, changes in operating conditions of the individual units have been discouraged or avoided for fear of the adverse effect on the overall system.
Changes in the operating conditions of individual units have been avoided for another reason. Any attempt to change the operating conditions at an individual unit would necessarily render the LP and its solution obsolete since the yield coefficients in the LP for a particular unit would change. Such a change would therefore require the LP to be solved again and prior art computer techniques for solving LPs are extremely cumbersome and complex. In the case of a digital computer, a great deal of computer time is required to solve an LP, particularly where the LP described a complex processing system.
Because of the foregoing difficulties associated with making the changes in operating conditions in the prior art, individual local objective functions consistent with the overall objective function of the system have been assigned and every effort has been made to maintain those operating conditions at the individual units which will satisfy the local objective functions rather than change to better operating conditions. In other words, no effort is made to deliberately change the operating conditions at the units to optimize the overall system.
Sensitivity analyses such as that disclosed by C. S. Bightler and D. J. Wilde, Hydrocarbon Processing, 44, No. 2, 111 (1965) have been proposed to determine the effect of changes in LP constraints on the operation of a system. However, no specific method has been proposed to take advantage of the sensitivity analysis for changes in the yield coefficients of an LP which optimizes the overall system, i.e., making changes in the operating conditions which will optimize the overall system as indicated by the sensitivity analysis.
SUMMARY OF THE INVENTION It is an object of this invention to provide an improved method of and apparatus for controlling the operation of a processing system including a plurality of individual processing units supplied by a plurality of feed streams of materials being transformed by the units into a plurality of product streams flowing from the units where the marginal product yield and the marginal product cost for each unit are dependent upon the operating conditions.
It is a more specific object of this invention to provide an improved method of and apparatus for controlling the operation of the processing system in a manner so as to encourage and implement changes in operating conditions for the individual units even though those changes do require a change in the linear program model of the processing system.
In accordance with these objects, a preferred embodiment of the invention comprises feed stream computer means for generating initial feed stream signals representing the initial feed stream flow rates for a given yield and cost at each of the units under a given set of operating conditions so as to maximize the profitability of the overall system. Feed stream control means are coupled to the feed stream computer means for controlling the feed stream flow rates in response to the initial feed stream signals. Unit yield and cost computer means generate new yield and cost signals representing new yield and cost for one of the units corresponding to an increase in the profitability of the overall system. Unit control means are coupled to the unit yield and cost computer means for controlling the op erating conditions at the one unit in response to the new yield and cost signals. Another feed stream computer means is coupled to the unit yield and cost computer means for repeatedly and continuously generating new feed stream signals representing new feed stream flow rates in response to the new yield and cost signals. The other feed stream computer means is coupled to the feed stream control means for controlling the feed stream flow rates in response to the new feed stream signals.
It is also a specific object of this invention to provide an improved method of and apparatus for operating and utilizing a linear program model.
In accordance with this specific object, the feed stream computer means includes means for computing the feed stream flow rates x from the linear program model having an overall system objective function of maximizing the profitability.
subject to x +....p x =Q x 2 Ofor alli wherep .Q=
Q is the marginal profit coefficient of the jth unit, x is the flow rate of the ith stream, p j is the yield column of the jth unit and a function of operating conditions, Q is the demand constraint column. a is the yield coefficient of the ith feed stream producing the jth product stream unit, and b,- is the demand coefficient of the ith product stream.
It is another specific object of this invention to assure that all changes in operating conditions at the units improve the overall profitability of the system.
In accordance with this specific object, the unit yield and cost computer means includes means for computing a local objective function where Ac, is the change in the marginal profit coefficient for the mth unit, A is the change in the mth yield column for the mth unit and g" are the simplex multipliers [c. c,,] at the previous operating conditions where Q" is the inverse of the basis matrix in [p,,
- pill- Itis a further specific object of this invention to provide an improved method of and apparatus for correcting the linear program model after permitting a change in operating conditions at one of the operating units.
In accordacne with this specific object, the other feed stream computer means includes means for computing the feed stream flow rates where BRIEF DESCRIPTION OF THE DRAWINGS FIG. 4 is a simplified refinery system operated in accordance with the method of this invention;
FIG. 5 is a graphical solution of a local optimization problem;
FIG. 6 is a block diagram of a processing system operated in accordance with the method depicted in FIG.
FIGS. 7 (a-d) are schematic circuit diagrams of a local computer and process control shown in block form in FIG. 6 where FIG. 7a shows circuitry for computing unit yield and cost coefficients, FIG. 7b shows circuitry for computing the local constraints and the local objective function, and FIGS. 7(c and d) shows circuitry for correcting the linear program model",
FIGS. 8 (a-c) are schematic circuit diagrams of the system computer and control shown in block form in FIG. 6 where FIG. 8a shows circuitry for computing the solution to the linear program model, FIG. 8b shows circuitry for determining the bases in the linear program model and FIG. 8c shows circuitry for controlling the feed stream flow rates to the units;
FIG. 9 is a flow diagram for a programmed digital computer which performs the same correction of the linear program performed by the circuitry of FIGS. 7(0- FIG. 10 is a block diagram of a typical refinery system which is controlled in accordance with this invention;
FIG. 11 is a block diagram of a fluid catalytic cracking unit of the system shown in FIG. 6; and
FIG. 12 is a block diagram depicting the relation and interconnections between the circuit diagrams of FIGS. 7(a-d) and FIGS. 8(ac).
TABLE OF CONTENTS FOR THE DETAILED DESCRIPTION I. Complex Processing System Employing Invention II. Method of Operating Complex Processing System III. General Description of a Simple Processing System IV. Detailed Description of the Simple Processing System Including Analog Computer Control A. Computing Yield and Cost Coefficients at the Local Computer B. Computing an LP Solution at the System Computer C. Computing the Bases of the LP Solution at the System Computer D. Setting the Flow Rate Controls E. Computing and Storing the Inverse Basis Matrix F. Computing the Simplex Multipliers G. Checking Local Optimization H. Checking Local Constraints I. Checking the Local Objective Function j. Correcting the LP Solution 1. Computing Sensitivity Coefficients 2. Computing New X's 3. Computing New Basis Matrix 4. Computing New Simplex Multipliers 5. Computing New Cost Coefficients V. Control of the Simple Processing System with a Digital Computer A. Correction of the LP B. Computing the Local Objective Function VI. Numerical Examples VII. Method of Operating a Complex Refinery System. VIII. Modified Method of Operating a Complex System.
|Cited Patent||Filing date||Publication date||Applicant||Title|
|US3075700 *||Jun 6, 1960||Jan 29, 1963||Controller|
|US3079079 *||May 21, 1958||Feb 26, 1963||Thompson Ramo Wooldridge Inc||Process and control apparatus for improving operating efficiency|
|US3201572 *||May 8, 1961||Aug 17, 1965||Du Pont||Control method and apparatus|
|US3594559 *||Jun 26, 1968||Jul 20, 1971||Phillips Petroleum Co||Process control for polymerization control system having equation updating feedback networks|
|US3621217 *||Sep 9, 1969||Nov 16, 1971||Gulf Research Development Co||Special-purpose process control computer|
|Citing Patent||Filing date||Publication date||Applicant||Title|
|US4349869 *||Oct 1, 1979||Sep 14, 1982||Shell Oil Company||Dynamic matrix control method|
|US4486148 *||May 17, 1982||Dec 4, 1984||Michigan Consolidated Gas Company||Method of controlling a motive power and fluid driving system|
|US4578747 *||Oct 14, 1983||Mar 25, 1986||Ford Motor Company||Selective parametric self-calibrating control system|
|US4607325 *||Mar 30, 1984||Aug 19, 1986||Honeywell Inc.||Discontinuous optimization procedure modelling the run-idle status of plural process components|
|US4616308 *||Dec 2, 1985||Oct 7, 1986||Shell Oil Company||Dynamic process control|
|US4660166 *||Jan 22, 1985||Apr 21, 1987||Bell Telephone Laboratories, Incorporated||Electronic network for collective decision based on large number of connections between signals|
|US4704584 *||Jun 25, 1986||Nov 3, 1987||Fairfield Industries||Instantaneous floating point amplifier|
|US4779347 *||Aug 10, 1987||Oct 25, 1988||Shell Oil Company||Method of determining the shape of a flat fabric pattern for application to a general curved surface|
|US4888692 *||Nov 10, 1988||Dec 19, 1989||Texas Instruments Incorporated||Real-time scheduling system|
|US5038270 *||Apr 24, 1989||Aug 6, 1991||Mitsubishi Kasei Corporation||Method for controlling reactor system|
|US5047947 *||Jul 25, 1990||Sep 10, 1991||Grumman Aerospace Corporation||Method of modeling the assembly of products to increase production yield|
|US5133073 *||May 29, 1990||Jul 21, 1992||Wavetracer, Inc.||Processor array of N-dimensions which is physically reconfigurable into N-1|
|US5155677 *||Mar 30, 1992||Oct 13, 1992||International Business Machines Corporation||Manufacturing process optimizations|
|US5157785 *||May 29, 1990||Oct 20, 1992||Wavetracer, Inc.||Process cell for an n-dimensional processor array having a single input element with 2n data inputs, memory, and full function arithmetic logic unit|
|US5193202 *||May 29, 1990||Mar 9, 1993||Wavetracer, Inc.||Processor array with relocated operand physical address generator capable of data transfer to distant physical processor for each virtual processor while simulating dimensionally larger array processor|
|US5195041 *||Jul 24, 1989||Mar 16, 1993||Institute Of Business Technology||Method and apparatus for improving manufacturing processes|
|US5265034 *||May 13, 1991||Nov 23, 1993||The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration||Feedback controlled optics with wavefront compensation|
|US5351195 *||Mar 12, 1993||Sep 27, 1994||The George Group||Method for improving manufacturing processes|
|US5375448 *||Mar 1, 1993||Dec 27, 1994||Hitachi, Ltd.||Non-interference control method and device|
|US5465321 *||Apr 7, 1993||Nov 7, 1995||The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration||Hidden markov models for fault detection in dynamic systems|
|US5696696 *||Oct 11, 1994||Dec 9, 1997||Betzdearborn, Inc.||Apparatus and method for automatically achieving and maintaining congruent control in an industrial boiler|
|US5727128 *||May 8, 1996||Mar 10, 1998||Fisher-Rosemount Systems, Inc.||System and method for automatically determining a set of variables for use in creating a process model|
|US5732193 *||Jan 20, 1995||Mar 24, 1998||Aberson; Michael||Method and apparatus for behavioristic-format coding of quantitative resource data/distributed automation protocol|
|US5752007 *||Mar 11, 1996||May 12, 1998||Fisher-Rosemount Systems, Inc.||System and method using separators for developing training records for use in creating an empirical model of a process|
|US5923571 *||Oct 6, 1997||Jul 13, 1999||Betzdearborn, Inc.||Apparatus and method for automatic congruent control of multiple boilers sharing a common feedwater line and chemical feed point|
|US6064819 *||Apr 2, 1997||May 16, 2000||Imec||Control flow and memory management optimization|
|US6236894 *||Dec 19, 1997||May 22, 2001||Atlantic Richfield Company||Petroleum production optimization utilizing adaptive network and genetic algorithm techniques|
|US6936665||Sep 25, 2002||Aug 30, 2005||Bp Corporation North America Inc.||Integrated advanced chemical process control|
|US7092893||Jan 28, 2003||Aug 15, 2006||Praxair Technology, Inc.||Control of liquid production of air separation plant network|
|US7133784||Oct 28, 2003||Nov 7, 2006||Bp Corporation North America, Inc.||Process and method for chemical manufacturing using transformation of on-line instrumentation data|
|US8069122||Mar 20, 2008||Nov 29, 2011||Accenture Global Services Limited||Predictive cost reduction based on a thermodynamic model|
|US8099259 *||May 24, 2005||Jan 17, 2012||Intercat Equipment, Inc.||Method for monitoring catalyst requirements of a refinery|
|US8775138 *||Oct 15, 2009||Jul 8, 2014||Exxonmobil Chemical Patents Inc.||Methods for handling withdrawal of streams from a linear programming model developed from a thermodynamically-based reference tool|
|US20030073787 *||Sep 25, 2002||Apr 17, 2003||Stephens William D.||Integrated advanced chemical process control|
|US20040133363 *||Oct 28, 2003||Jul 8, 2004||Ramaswamy Vaidyanathan||Process and method for chemical manufacturing using transformation of on-line instrumentation data|
|US20040148264 *||Jan 28, 2003||Jul 29, 2004||Lawrence Megan||Control of liquid production of air separation plant network|
|US20050216209 *||May 24, 2005||Sep 29, 2005||Intercat Equipment, Inc.||Method for monitoring a FCC catalyst injection system|
|US20080235067 *||Mar 20, 2008||Sep 25, 2008||Accenture||Predictive Cost Reduction Based On A Thermodynamic Model|
|US20100101272 *||Mar 19, 2007||Apr 29, 2010||Coward Brian A||process of liquefying a gaseous methane-rich feed for obtaining liquid natural gas|
|US20100131250 *||Oct 15, 2009||May 27, 2010||Carpency Joseph F||Methods for Handling Withdrawal of Streams from a Linear Programming Model Developed from a Thermodynamically-Based Reference Tool|
|US20110213715 *||Mar 1, 2010||Sep 1, 2011||Accenture Global Services Gmbh||Lean Analytics|
|US20140095268 *||Sep 28, 2012||Apr 3, 2014||Avaya Inc.||System and method of improving contact center supervisor decision making|
|CN101125291B||Sep 25, 2002||Oct 13, 2010||伊内奥斯美国公司||Integrated advanced chemical process control|
|EP1921406A1 *||Nov 8, 2006||May 14, 2008||Honeywell Control Systems Ltd.||A process of liquefying a gaseous methane-rich feed for obtaining liquid natural gas|
|WO1985001807A1 *||Oct 14, 1983||Apr 25, 1985||Ford Motor Company||Selective parametric self-calibrating control system|
|WO1993022725A1 *||Apr 28, 1993||Nov 11, 1993||THE UNITED STATES GOVERNMENT as represented by THESECRETARY OF COMMERCE||Synthetic perturbation tuning of computer programs|
|WO2001008054A2 *||Jul 27, 2000||Feb 1, 2001||Raytheon Company||Method and system for process design|
|WO2001008054A3 *||Jul 27, 2000||Mar 14, 2002||Raytheon Co||Method and system for process design|
|WO2003026791A1 *||Sep 25, 2002||Apr 3, 2003||Bp Corporation North America Inc.||Integrated chemical process control|
|WO2008056269A1 *||Mar 19, 2007||May 15, 2008||Honeywell Control Systems Ltd.||A process of liquefying a gaseous methane-rich feed for obtaining liquid natural gas|
|WO2017001173A1 *||Jun 13, 2016||Jan 5, 2017||Siemens Aktiengesellschaft||Control device for a production module, production module having a control device, and method for operating the control device|
|U.S. Classification||700/3, 700/36, 700/38, 700/31, 700/29|
|International Classification||G05B13/02, B01J19/00, G06F15/18|
|Cooperative Classification||B01J19/0033, G05B13/021|
|European Classification||B01J19/00B6, G05B13/02A1|