# Patent

## Title
Systems and computer programs for simultaneous process and utility systems synthesis in partially and fully decentralized environments

## Abstract

Systems and program code for providing a combined synthesis of cost-effective heat and power (CHP) utility systems and the key process systems are provided. An exemplary system includes a CHP utility and process subsystems analysis and design program computer and program code for identifying the best key design and operating conditions for both subsystems stored in the memory of the computer. The program code, when executed, can cause the computer to perform operations including determining an optimal allocation of steam and power between both process and utility systems for one or more industrial process facilities in both partially and totally decentralized environments.

## Background
RELATED APPLICATIONS
This non-provisional patent application claims priority to and the benefit of U.S. Provisional Patent Application No. 61/612,470, filed on Mar. 19, 2012, titled “System, Method, and Computer Program For Simultaneous Processing Utility Systems Synthesis in Partially and Fully Decentralized Environments,” and is related to U.S. Non-Provisional patent application titled “Methods for Simultaneous Process and Utility Systems Synthesis in Partially and Fully Decentralized Environments,” filed on Feb. 1, 2013, as U.S. application Ser. No. 13/757,467, each incorporated by reference in its entirety. See Appendix 1 for a list of additional related applications.


BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates in general to the field of energy management systems, and in particular, to systems, computer program code, and methods related to simultaneously synthesizing cost-effective combined heat and power utility systems and the process plant's heat exchangers network and distillation network sequence.
2. Description of the Related Art
Many different types of processes consume multiple steam levels and electricity to obtain an output result, or to produce a required product or compound. For large-scale processes that consume significant amounts of steam, it is preferable to optimize the consumption of energy through careful operation, design or reconfiguration of the plant and the equipment used. Further, in some industrial manufacturing processes, specific streams of material flows need to be supplied to different types of equipment and machinery at specific temperatures. These material flows may need to be heated or cooled from an original starting or supply temperature to a target temperature. This, in turn, will require the consumption of steam to heat specific streams and consumption of water, for example, to cool down specific streams. Energy to heat or cool down the specific streams can come from utility plant systems and/or through a network of heat exchange units which can recover potentially lost energy, reducing the amount of energy needed from the utility plant systems.
In nowadays mega complexes, the cost of fuel, power and its infrastructure as well as the site process heat exchangers network (HEN) and distillation network can be huge. A total site utility system in a chemical processing site, interlinked to the chemical plants production processes, for example, provides heat and power to chemical production processes. The economic performance of the total site utility systems directly influences the operating cost of the chemical processing site, hence it is recognized by the inventors that improvements in the synthesis of process subsystems and the total site utility systems can lead to significant savings.
Two basic approaches have been developed for the synthesis of utility systems (CHP). The first one is based on thermodynamic performance of the plant and aims in maximizing the plant's overall energy efficiency. Some optimization techniques may be included as part of the solution procedures. The second approach only uses optimization methods such as, for example, MILP and MINLP to satisfy the predefined targets/demands. In the methods that have been employed in both commercial software (e.g. Aspen Utilities Planner of Aspentech, Star software of university of Manchester/PI department, and ProSteam of KBC) and academic literatures, the task has been to design a utility plant that satisfies given/specific steam and power demands. This approach was later extended to account for uncertainty and disturbance in such one-time given/specific power and steam demand to multi-periods of again specific/given power and steam demand.
The state-of-the-art methodology presently adopted by the industrial community (e.g., process engineering departments) and academia, is also employed without systematic consideration of both distillation columns subsystem best sequence, design and operating conditions, and HEN synthesis, as well as best design conditions using high fidelity utility CHP details. The inventors have recognized that such approach is performed without considering the process degrees of freedom in the design of HEN, large distillation network sequence, or small but important groups of distillation network sequence or using the best process conditions to enable the best synthesis of the utility system. That is, recognized by the inventors is that the state-of-the-art fails to recognize or employ simultaneous syntheses of such subsystems, CHP, distillation columns, HEN, and/or the rest of process conditions. The numerous potential combinations of the degrees of freedom in both utility system and the key process systems such as, for example, HEN and process distillation network sequence and process design and operating conditions, however, can allow significant opportunity for both systems' syntheses optimization.
The combination of the two approaches as been demonstrated. Neither one of the combined approaches addressed the simultaneous synthesis of key process systems, such as, for example, the distillation sequence system, the HEN system, much less taking into consideration all possible combinations of process changes in the rest of the process such as temperature, pressure, reflux ratio, inter-coolers, and inter-heaters, in distillation columns, and so on, with the utility system synthesis.
In order to try to “optimize” energy recovery in both the process and the utility plants, the state-of-the-art has generally attempted to solve the synthesis problem by utilizing two entirely decomposed optimization problems: one for the chemical/process plant and one for the utility plant. Those two problems are currently solved in most of the companies by two different teams, either located in two separately decentralized departments within the same company, or in two different companies with the first being in the process plant and the second being in the utility plant. According to such conventional methodology, the process synthesis team/group dictates the process steam demand and process steam generation as well as equipment driver's type (motor or steam) to the utility subsystem (CHP) synthesis team/group, ahead of the CHP synthesis.
As such, the utility system synthesis problem is always addressed in both industry and academia as a “follower” objective problem. The process is always synthesized first and its process conditions are set. Thereafter the desired steam quantity and quality are defined to the utility system. The decision for process liquid and gas driving using motors driven or steam turbine driven pumps and compressors are decided ahead of the utility system synthesis. Such arguably archaic approach applied by managers in the decentralized company office emanates from an adherence to decisions made in the past, where the utility system was merely a bunch of boilers and the major capital in the facility was the process systems (e.g., distillation and HEN subsystems).
Recognized by the inventors, however, is that as a result of the new era of cogeneration, tri-generation and even quadra-generation where power and water emerged as new products, it may be more economically prudent to “change the chairs” such that the process systems, such as the HEN and distillation network, become the “follower” and the “utility system” becomes the “leader”. The inventors further recognize that the “utility system synthesis guy” can provide better results where his counterpart “process synthesis guy” cooperates, especially if the capital cost of the utility system is much higher than that of HEN and/or distillation network. Naturally, effective management of such utility systems syntheses can render big impact on projects' profitability and life cycle success.
Nevertheless, also recognized is that even in case of the closeness in capital cost between the utility system and the process system, the problem may become one of conflicting design objectives types. The problem can also become a little more complicated when the designer and/or the owner of the central utility system providing utilities to the chemical complex is a “third party” (other than the chemical company owner) who may have an investment plan of his own and limited project budget, and/or may represent the country national power grid. In such case, conflicting objectives may arise that may need to be addressed and solved between the chemical complex owner and the utility system owner in order to produce a result which satisfies both parties. Hence, it is recognized that there is a need for a different solution approach than the currently available in both industry and academia.
Correspondingly, recognized is the need for systems, computer program, and methods synthesizing cost-effective combined heat and power (CHP) utility systems and the process plant's heat exchangers network and distillation network sequence, which can identify its/the best key design and operating conditions in both partially and totally decentralized environments in order to address the problem of obtaining an optimal tradeoff among the CHP utility system synthesis and the key process units (HEN and distillation sequences) syntheses in mega facilities, which require many complex and huge interactions to be explored for best economic decisions.
SUMMARY OF THE INVENTION
In view of the foregoing, various embodiments of the present invention advantageously provide systems, computer programs, computer readable media, and methods of providing a combined synthesis of cost-effective heat and power (CHP) utility systems and the key process systems (e.g., process plant's heat exchangers network and distillation network sequence), which can include identifying its best key design and operating conditions in both partially and totally decentralized environments. Advantageously, various embodiments of the present invention address the problem of obtaining an optimal tradeoff among the CHP utility system synthesis and the key process units (heat exchanger network (HEN) and distillation sequences) syntheses in mega facilities, which require many complex and huge interactions to be explored for best economic decisions.
One particular aspect where current approaches have an important limitation is in accounting explicitly for the interactions that take place when synthesizing total processing systems that include several major components such as hot utility system, HEN and chemical plant processing units. Various embodiments of the present invention can advantageously overcome such limitation by allowing virtual simultaneous selection of structural, design and operating parameters of the hot utility system, HEN and the “key” process units (distillation columns).
According to one or more embodiments of the present invention, a unique formulation is presented by a CHP utility and process system/subsystem mathematical model, which can advantageously explicitly account for the interactions that take place in the syntheses of CHP, HEN, and distillation networks/key process units and systems. The solved problem's degrees of freedom can advantageously include: the utility system super structure, the HEN, the process operating parameter/conditions, the process design parameters/conditions, and the structural parameters for some key process units such as distillation, absorption, adsorption, stripping and its operating pressure. The degrees of freedom can also include the reflux ratio, number of stages, and the columns intra-heaters and intra-coolers, as well as the different potential sequences. The degrees of freedom can further include different sources of fuel, boilers, gas turbines, HRSGs, and steam turbines different configurations (each having minimum and maximum levels).
For example, various embodiments of the present invention advantageously provide systematic systems, computer programs, computer readable media, and methods for the simultaneous syntheses of combined heat and power (CHP) utility and process subsystems for a facility or cluster of process facilities, e.g., usable by chemical process industries. One or more exemplary embodiments include a combined heat and power utility and process subsystems analysis and design program computer having a processor, and memory coupled to the processor to store software and database records therein, and combined heat and power utility and process subsystems analysis and design program code defining a computer program stored in the memory of the combined heat and power utility and process subsystems analysis and design program computer to perform simultaneous synthesis of combined heat and power utility and process subsystems for a process facility or cluster of process facilities.
The computer program can include including instructions that when executed by the combined heat and power utility and process subsystems analysis and design program computer, cause the computer to perform the following the operations: determining an optimal allocation of steam and power between both process and utility systems for one or more industrial process facilities, and determining a plurality of at least substantially optimal process and utility subsystem conditions. The operations can include: defining a plurality of rigorously bounded process variables for one or more process plants/facilities as decision variables to thereby form model input data for a CHP utility and process subsystem mathematical program model, processing the model input data by the CHP utility and process subsystem mathematical program model, and determining the plurality of at least substantially optimal process and utility subsystem conditions. The decision variables can include a process steam demand interval at each of a plurality of different steam pressure levels and/or a process power generation interval, with the interval describing the respective upper and lower bounds of the respective variables. The process conditions can include a process steam demand load value at each of the plurality of different steam pressure levels and/or a process power generation value for the one or more process facility defining model output. The utility subsystem conditions include combined heat and power utility subsystem steam and/or power generation and allocation values. The CHP utility and process subsystem mathematical program model can include a cost function configured to minimize total cost of the utility and process subsystems of the one or more process facilities, to thereby determine the at least substantially optimal process and utility subsystem conditions.
According to one or more embodiments, the operation of determining a plurality of at least substantially optimal process and utility subsystem conditions includes one or more of the following operations: determining CHP utility subsystem steam and power generation and allocation values, determining a number of steam headers for the modeled utility subsystem and operating conditions of each steam header, determining a number of boilers for the modeled utility subsystem and designated capacities of each boiler, determining a number of steam turbine generators for the modeled utility subsystem and capacity of each steam turbine generator, determining a number of cogeneration units for the modeled utility subsystem and designated capacities of each cogeneration unit, determining a number of electrical motors and steam turbines for process equipment driving and size of each, determining a size and location of solar power system for the modeled utility subsystem, and determining an optimal allocation of steam and power between both process and utility subsystems, whereby the optimal utility subsystem solution is solved first and the process subsystem conditions solution is solved second utilizing portions of the optimal utility subsystem solution as a basis thereof, or alternatively, the process subsystem conditions solution is solved first.
The operation of determining a plurality of at least substantially optimal process and utility subsystem conditions can also or alternatively include one or more of the following operations: identifying process conditions which render a minimum number of HEN units for the one or more process facilities for all reasonably possible distillation columns sequences, and identifying process conditions which render a minimum total heat exchanger surface area of a HEN network design for the one or more facilities and define the surface area of each heat exchanger unit in the HEN. The heat exchanger surface area identification operation can be performed, for example, under/for one or more of the following: one or more reasonably possible combinations of steam generation and demand loads' levels and/or one or more given or desired combinations of steam generation and demand loads' levels, one or more reasonably possible distillation columns sequences, one or more reasonably possible distillation columns sequences in conjunction with one or more possible combinations of steam generation and demand loads' levels, one or more reasonably possible distillation columns sequences and one or more given or desired combinations of steam generation and demand loads' levels, one or more reasonably possible HEN configurations and one or more reasonably possible combinations of steam generation and demand loads' levels, and one or more reasonably possible HEN configurations and one or more given or desired combinations of steam generation and demand loads' levels.
One or more embodiments of the systems, computer program, and/or methods, can advantageously include/perform the following operations: providing rigorously bounded decision variables to a combined heat and power mixed integer linear program model for processing by a computer, and determining at least substantially simultaneously by the computer, an optimal allocation of steam and power between both process and combined heat and power utility systems for one or more industrial process facilities responsive to the rigorously bounded decision variables and the combined heat and power mixed integer linear program model. The operation of determining can include providing, by the combined heat and power mixed integer linear program model, combined heat and power utility subsystem steam and power generation and allocation values and process system steam generation, steam demand loads and levels, power generation, and power export values as output variables calculated responsive to the rigorously bounded decision variables.
One or more other exemplary embodiments of the systems, computer program, and/or methods can also or alternatively include/perform synthesizing preparation and execution operations. The operations can include selecting an optimal fuel type/energy source mix (e.g., coal, heavy fuel oil, natural gas, biomass, waste materials, solar etc.), and selecting equipment which generate steam (e.g., for heating purposes, pumps and compressors driving, heat carrying, cleaning, cooling) and power (e.g., for lighting and other applications). The operations can also or alternatively include defining the number of steam headers and its conditions (pressure or saturation temperature), determining the shaft work network configuration, and/or determining the optimal allocation of steam and power to both process and utility plants usage in both partially and totally decentralized environments, under/for various conditions or constraints.
The various conditions/constraints can include one or more of the following: all reasonably possible combinations of process design parametric changes (e.g., pre-identified potential parameter values) using all generally available sources of fuel/energy available to the utility system provider, all reasonably possible combinations of process design and operation parametric changes using all generally available sources of fuel or energy available to the utility system provider, all reasonably possible distillation sequences using all generally available sources of fuel or energy available to the utility system provider, all reasonably possible distillation sequences and under all reasonably possible combinations of process design and operation parametric changes using all generally available sources of fuel or energy available to the utility system provider, all reasonably possible HEN designs using all generally available sources of fuel or energy available to the utility system provider, for all reasonably possible HEN designs under all reasonably possible combinations of process design and operation parametric changes using all generally available sources of fuel or energy available to the utility system provider, all reasonably possible distillation columns sequences and HEN designs using all generally available sources of fuel or energy available to the utility system provider, and all reasonably possible distillation columns sequences and HEN designs under all reasonably possible combinations of process design and operation parametric changes using all generally available sources of fuel or energy available to the utility system provider.
According to one or more embodiments of the systems, computer program, and/or methods, the systematic system-based methodology and computer program can also advantageously further or alternatively include provisions for determining the process plant's HEN network design and distillation network sequence with its main design and operating conditions under all reasonably possible combinations of process operation parametric changes using, for example, all generally available sources of fuel/energy for fully decentralized and/or partially decentralized process and utility systems design environments responsive to receiving and/or determining a plurality of sets of values each defining a potential range of attribute values for defining a separate one of the corresponding plurality of process operation parametric changes.
One or more embodiments of the systems, computer program, and/or methods can also or alternatively provide systematic systems, computer programs, computer readable media, and methods for the identification of process conditions which render the minimum number of plant's HEN units for all reasonably possible distillation columns sequences.
One or more embodiments of the systems, computer program, and/or methods can also alternatively provide systematic systems, computer programs, computer readable media, and methods for the identification of process conditions which render minimum total surface area of plant's HEN network design and define the surface area of each unit in the network: for possible combinations of steam generation and demand loads' levels, for given/desired combinations of steam generation and demand loads' levels, for possible distillation columns sequences, for possible distillation columns sequences and possible combinations of steam generation and demand loads' levels, for possible distillation columns sequences and for given/desired combinations of steam generation and demand loads' levels, for possible HEN configurations and possible combinations of steam generation and demand loads' levels, and/or for possible HEN configurations and given/desired combinations of steam generation and demand loads' levels. According to an exemplary process flow, optimal allocation of steam and power for the utility system solution is solved first and the allocation of steam and power for the process system conditions solution is solved second utilizing portions of the optimal utility allocation of steam and power as a basis thereof, although the reverse can be alternatively employed.



BRIEF DESCRIPTION OF THE DRAWINGS
So that the manner in which the features and advantages of the invention, as well as others which will become apparent, may be understood in more detail, a more particular description of the invention briefly summarized above may be had by reference to the embodiments thereof which are illustrated in the appended drawings, which form a part of this specification. It is to be noted, however, that the drawings illustrate only various embodiments of the invention and are therefore not to be considered limiting of the invention's scope as it may include other effective embodiments as well.
FIG. 1 is a high-level block system diagram illustrating an embodiment of a system for providing combined synthesis of cost-effective combined heat and power (CHP) utility systems and key process systems according to an embodiment of the present invention;
FIG. 2 is a block flow diagram illustrating a process for simultaneous processing utility systems synthesis in partially and totally decentralized environments according to an embodiment of the present invention;
FIG. 3 is a graphical illustration in the form of a Pareto diagram illustrating a method/process of determining a trade-off between process and utility capital costs according to an embodiment of the present invention;
FIGS. 4A-4B are graphical illustrations showing a schematic method/process for columns integration under possible combinations of reflux ratio and operating pressure modifications according to an embodiment of the present invention;
FIGS. 5A-5B are graphical illustrations showing a sequence list for complex and simple distillation networks examples according to an embodiment of the present invention;
FIG. 6 is a graphical illustration of a graphical user interface illustrating input data entry for a CHP Mathematical Program according to an embodiment of the present invention;
FIG. 7 is a schematic block diagram of a group sequences list for a simple distillation network according to an embodiment of the present invention;
FIGS. 8A-8B are schematic diagrams of an interval cascade diagram and interval grand composite curves, respectively, according to an embodiment of the present invention;
FIG. 9 is a graphical illustration of a combined set of grand composite curves according to an embodiment of the present invention;
FIGS. 10A-10B are graphical illustrations of grand composite curves used for selecting an optimal number of steam headers according to an embodiment of the present invention;
FIG. 11 is a block flow diagram illustrating a CHP number of headers selection algorithm according to an embodiment of the present invention;
FIG. 12 is a graphical interconnection diagram showing the number of cases to be evaluated based upon headers selection according to an embodiment of the present invention;
FIG. 13 is a table illustrating tabular input fields that provide for entry of interval data according to an embodiment of the present invention;
FIG. 14 is a table illustrating tabular input fields that provide for entry of equipment pressure and steam pressure type data according to an embodiment of the present invention;
FIG. 15 is a table illustrating tabular input fields that provide for entry of steam turbines and motors data according to an embodiment of the present invention;
FIG. 16 is a graphical illustration of a CHP program output module according to an embodiment of the present invention;
FIG. 17 is a table providing abbreviations used in the CHP program output module of FIG. 16 according to an embodiment of the present invention;
FIGS. 18-19 are temperature interval diagram illustrating selection of optimal process conditions and HEN minimum number of units identification according to an embodiment of the present invention;
FIGS. 20-21 are temperature interval diagrams illustrating selection of optimal process conditions and HEN minimum number of units identification according to an embodiment of the present invention;
FIG. 22 is a graphical user interface illustrating an optimal process conditions and HEN area identification projects listing main page according to an embodiment of the present invention;
FIG. 23 is a graphical depiction illustrating an optimal process conditions and HEN area identification screen providing project details according to an embodiment of the present invention;
FIGS. 24A-B is a graphical illustration providing an example of project details editing according to an embodiment of the present invention;
FIG. 25 is a graphical illustration of a network topology according to an embodiment of the present invention;
FIG. 26 is a graphical illustration of an optimal process conditions and HEN area identification enter data set according to an embodiment of the present invention;
FIG. 27 is a graphical illustration of an optimal process conditions and HEN area identification model formulation according to an embodiment of the present invention;
FIG. 28 is a graphical illustration of solution results of the model formulation according to an embodiment of the present invention;
FIG. 29 is a graphical illustration of an exemplary configuration of a HEN solution example according to an embodiment of the present invention;
FIG. 30 is a graphical illustration of the HEN solution example of FIG. 29 having process conditions and HEN surface area solved employing the process as the “Leader” and the Utility as the “Follower” according to an embodiment of the present invention;
FIG. 31 is a graphical illustration of the HEN solution example of FIG. 29 having process conditions and HEN surface area solved employing utility as the “Leader” and the process is the “Follower” according to an embodiment of the present invention;
FIG. 32 is a graphical illustration of an exemplary configuration of a HEN solution example according to an embodiment of the present invention;
FIG. 33 is a graphical illustration of the HEN solution example of FIG. 32 having process conditions and HEN surface area solved employing the process as the “Leader” and the Utility as the “Follower” according to an embodiment of the present invention;
FIG. 34 is a graphical illustration of the HEN solution example of FIG. 32 having process conditions and HEN surface area solved employing utility as the “Leader” and the process is the “Follower” according to an embodiment of the present invention;
FIG. 35 is a graphical illustration of the exemplary HEN solution example of FIG. 29 having a larger ΔT_min_i, larger resulting in larger utility requirements according to an embodiment of the present invention;
FIG. 36 is a graphical illustration of the HEN solution example of FIG. 35 having process conditions and HEN surface area solved employing the process as the “Leader” and the Utility as the “Follower” according to an embodiment of the present invention;
FIG. 37 is a graphical illustration of the HEN solution example of FIG. 35 having process conditions and HEN surface area solved employing utility as the “Leader” and the process is the “Follower” according to an embodiment of the present invention;
FIG. 38 is a graphical illustration of the exemplary HEN solution example of FIG. 29 having an even larger ΔT_min_i, larger resulting in even larger utility requirements according to an embodiment of the present invention;
FIG. 39 is a graphical illustration of the HEN solution example of FIG. 38 having process conditions and HEN surface area solved employing the process as the “Leader” and the Utility as the “Follower” according to an embodiment of the present invention; and
FIG. 40 is a graphical illustration of the HEN solution example of FIG. 38 having process conditions and HEN surface area solved employing utility as the “Leader” and the process is the “Follower” according to an embodiment of the present invention.



DETAILED DESCRIPTION
The present invention will now be described more fully hereinafter with reference to the accompanying drawings, which illustrate embodiments of the invention. This invention may, however, be embodied in many different forms and should not be construed as limited to the illustrated embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art. Like numbers refer to like elements throughout. Prime notation, if used, indicates similar elements in alternative embodiments.
In the synthesis of total site utility systems, it is often the case that considerable changes exist between total site utility systems and chemical processes. Therefore, the integration of chemical processes and total site utility systems needs to be considered in the synthesis problem. Conventional wisdom in that field is to sequentially synthesize the process units/system first and the utility system second, in decentralized departments. Such approach can result in severe economic loss in both the process units and utility system design. One major technical reason behind such myopic approach is that, the synthesis of the combined heat and power utility system, alone, is a complex combinatorial problem that has multiple objectives. Adding the complexity of heat exchanger network (HEN) synthesis to the complexity of this synthesis problem, which is, itself, mathematically also considered NP-hard along with the synthesis problem of selecting optimal plant's distillation sequence, would result in the problem becoming intractable.
The synthesis of the combined heat and power (CHP) utility system involves four major tasks: (1) the selection of fuel type/energy source mix (e.g., coal, heavy fuel oil, natural gas, biomass, waste materials, renewable/solar etc.); (2) the hardware/devices/equipment which generate steam (e.g., mainly for heating, driving, heat carrying, cleaning, and even cooling) and power; (3) the number of steam headers and its conditions (e.g., pressure or saturation temperature); and (4) the determination of the shaft work generation structural configuration; as well as the optimal allocation of such two commodities (e.g., steam and power) to both process and utility plants usage (e.g., heating loads allocation, drivers selection, cooling tasks, and so on). The inventors have recognized that such decision variables interact with the in-process steam and power demand and generation. The number of steam levels, and its placement are also recognized to have significant impact, not only on the steam headers design and cogeneration potential of the synthesized CHP system including the steam turbines efficiencies variation and its size selection, but also on in-process key units such as the distillation and stripping network structural and parametric designs, as well as the heat exchangers network synthesis.
Conventionally, there are no systems, computer programs, computer readable media, or methods that simultaneously synthesize the CHP system with the interlinked key process systems such as heat exchangers network and distillation sequence. Although commercial software does exist to synthesize CHP systems, such as, for example, Star, Pro Steam, and others, they do not provide such capability. In academia, the big HEN problems are considered NP-hard. Making the problem even more combinatorial via the addition of distillation sequences and CHP synthesis, would make the problem an intractable problem to solve, especially for a mega industrial complex, which can have as many as 25 plants on the single site. Beneficially, various embodiments of the present invention provide a systematic methodology for simultaneous synthesis of: cost-effective combined heat and power utility system; process plant heat exchanger network; and distillation network sequence, with identification of substantially optimal key design and operating conditions. Various embodiments of the present invention also can account for the interactions between the CHP, HEN, and key process units including the distillation columns, can synthesize the utility system in consideration of the degrees of freedom and design of the process system (e.g., HEN and distillation network sequence), and can determine or facilitate the determination of an optimal trade-off between costs associated with utility system and costs associated with the process system.
System Information
Various system options according to exemplary embodiments of the present invention are described below.
FIG. 1, for example, illustrates a system 30 for providing a cost-effective combined heat and power (CHP) utility system, process plant heat exchanger network (HEN or HEN network), and distillation network sequence, with identification of substantially optimal key design and operating conditions under all possible combinations of process and utility parameters variations for an industrial process facility (individual process facility or cluster of industrial process facilities). The system 30 can include a CHP utility and process subsystems analysis and design computer 31 having a processor 33, memory 35 coupled to the processor 33 to store software and database records therein, and a user interface 37 which can include a graphical display 39 for displaying graphical images, and a user input device 41 as known to those skilled in the art, to provide a user access to manipulate the software and database records. Note, the computer 31 can be in the form of a personal computer or in the form of a server or server farm serving multiple user interfaces 37 or other configuration known to those skilled in the art. Accordingly, the user interface 37 can be either directly connected to the computer 31 or through a network 38 as known to those skilled in the art.
The system 30 can also include a database 43 stored in the memory 35 (internal or external) of the CHP utility and process subsystems analysis and design computer 31 and having a plurality of sets of values each separately defining a potential range of values for each operational attribute for each of a plurality of the utility and process subsections. Such attributes can include, for example, fuel type/energy source (coal, heavy fuel oil, natural gas, biomass, waste materials, solar etc.), equipment which generate steam (for heating purposes, pumps and compressors driving, heat carrying, cleaning, cooling) and power (for lighting and other applications), steam headers and its range of conditions (pressure or saturation temperature), shaft work network configuration, and the range of values with respect to allocation of steam and power to both process and utility plants usage, and discrete values identified as providing optimal and potentially optimal results.
Such attributes can also include a lower and an upper boundary value for the process supply temperature (Ts) of each of the hot process streams and each of the cold process streams, a lower and an upper boundary value for a target temperature (Tt) of each of the hot process streams and each of the cold process streams, and/or a lower and an upper boundary value for a heat capacity flow rate (FCp) of each of the hot process streams and each of the cold process streams, capital costs of various utility system components and various HEN equipment for the industrial facility according to the received streams conditions, along with one or more sets of stream-specific minimum temperature approach values between streams (ΔT_min_i), streams initial types, streams matching constraints, global utility consumption values [Qh], [Qc] where the “[ ]” denotes interval values. Such attributes can also include the interval and/or discrete locations of the pinch regions often referred to as a “pinch point” which describe a “region of minimum choice lower and upper temperature boundaries” at least for each pinch point controlling process stream/stream temperature, identification of the streams that control the pinch locations, data linking the pinch points define a map or maps of the pinch locations according to a progressive change in ΔT_min_i or process conditions, and the minimum number of HEN units required for a network condition at each pinch location, among others.
The system 30 can also include a CHP utility and process subsystems analysis and design program 51 stored in memory 35 of the CHP utility and process subsystems analysis and design computer 31 and adapted to provide a systematic process which includes various unique phases of analysis and design, which provide a revolutionary solution approach to providing an optimal solution for both utility and process systems/subsystems in both partially and fully decentralized environments.
The program 51 can incorporate one or more combinations of the following algorithms to include: a heating steam demand and generation upper and lower loads calculation algorithm, an in-process steam generation algorithm, an in-process steam demand algorithm, a distillation sequences list and its capital cost estimation algorithm(s), distillation columns reflux ratio, number of trays and pressure selection and costing algorithms, a HEN minimum area calculation and process conditions identification algorithm, a HEN minimum number of units and process conditions identification algorithm, a HEN minimum capital cost algorithm, a utility and process system cost comparison algorithm, an independent CHP utility system design and costing algorithms, a total steam generation and demand (level) lower and upper bounds calculations and/or targeting algorithm(s), a number of steam headers selection algorithm, a steam headers' conditions identifications algorithm, a process and distillation columns integration algorithm, and a CHP optimization algorithm. The set of algorithms can also include an energy targeting and process conditions identification algorithm described, for example, in U.S. Pat. No. 7,698,022, a {ΔT_min_i} marching module is described, for example, in U.S. Pat. No. 8,032,262 titled “System, Method, and Program Product For Synthesizing Non-Constrained Heat Exchanger Networks, a pinch points interval identification module is described, for example, in U.S. patent application Ser. No. 13/041,057, titled “Systems and Program Product For Heat Exchanger Network Energy Efficiency Assessment and Lifetime Retrofit, an interval-based site-wide steam demand and generation bounding module is described, for example, in U.S. patent application Ser. No. 12/480,415, titled “Method and Software For Global Targeting of Process Utilities Under Varying Conditions, and a HEN synthesis module is described, for example, in U.S. Pat. No. 7,729,809, titled “System, Method, and Program Product for Targeting and Identification of Optimal Process Variables in Constrained Energy Recovery Systems. FIG. 2 illustrates the interconnection between the major algorithms according to an embodiment of the present invention.
Note, the CHP utility and process subsystems analysis and design program 51 can be in the form of microcode, programs, routines, and symbolic languages that provide a specific set for sets of ordered operations that control the functioning of the hardware and direct its operation, as known and understood by those skilled in the art. Note also, the CHP utility and process subsystems analysis and design program 51, according to an embodiment of the present invention, need not reside in its entirety in volatile memory, but can be selectively loaded, as necessary, according to various methodologies as known and understood by those skilled in the art.
Serial Synthesis of Combined Heat and Power (CHP) Utility and Process Subsystems Having Total System Minimum Annualized Cost
Methodology for Partially Decentralized Application:
In such category we typically have a project with a single company objective of synthesizing both the process and the utility subsystems to have a total system minimum annualized cost, but with the syntheses tasks of process subsystems and utility subsystem (CHP) conducted in a decentralized approach. For example, the process design department and utility design department are working separately. Both departments are trying to achieve the best possible synthesis/design (minimum total annualized cost) for their respective subsystems, subject to their problem's constraints. According to the current state of art, the utility subsystem (CHP) designer gets his/her constraints from its own problem superstructure in addition to the decision variables of the process subsystems' optimization problem solution (i.e., for the distillation network and heat exchangers network (HEN) subsystems). The process problem is solved first by the process designer to provide the nominal design exact demand and generation values for steam levels and for power demand, as well as equipment to drive by power, and the utility subsystem problem (CHP) is solved later for a given desired process designer situation.
According to various embodiments of the present invention, three different analysis categories are recognized, which are dealt with differently. For example, according to an embodiment of the present invention, in order to determine the appropriate category for proper analysis, all viable process subsystems syntheses alternatives are first evaluated to determine the process subsystems' total annualized cost (PSTAC) values. The parameters of the promising alternatives are then used to calculate the steam demand and generation bounds for the utility subsystem synthesis (CHP) under all possible combinations of desired process subsystems synthesis design and operating conditions (distillation networks and HEN) using, for example, {ΔT_min−i} intervals. The CHP subsystem then starts synthesizing its subsystem using such steam bounds to get best total annualized cost (USTAC) of the utility subsystem.
Case I: Process Subsystems are the Leader and Utility Subsystem is the Follower
If the CHP USTAC is much less in dollar value compared to the PSTAC, the process subsystems become the problem objective leader and the CHP becomes the problem objective follower. According to the solution approach to this problem, the process designer solves the problem to his/her best interest first and dictates to the follower to find his/her best objective under his/her superstructure constraints and the decision variables, obtained from the process designer problem solution, given to him/her by the process subsystems designer. These constraints and decision variables include steam demand and generation levels and operating conditions and an ad hoc selection of process equipment drivers, such as turbines and/or motors, for pumps and compressors. This class of problem is currently implemented, at least to some extent, in CHP utility subsystem synthesis in both academia and industry.
Case II: Process Subsystems are the Follower and Utility Subsystem is the Leader
If the PSTAC is much less in dollar value compared to the CHP USTAC, the CHP becomes the problem objective leader and the process subsystems become the problem of objective follower. According to this solution, the process designer initially provides a broad set of constraints which the utility system designer utilizes to formulate its optimal solution. The output parameters regarding steam and power are then provided back to the process system designer to be utilized in designing its optimal solution given the utility system constraints.
Case III: Tie Cooperating Approach between Process and Utility Subsystems
If the PSTAC and the CHP USTAC are similar in costs, a Pareto-type diagram can be established to analyze the cheapest distillation sequence and heat exchange network versus CHP costs using a simple iteration between departments. A Pareto-type diagram for company-wide best decision-making can be employed using {ΔT_min_i} as a marching variable. An example of such a marching technique is described, for example, in U.S. Pat. No. 8,032,262.
Method for Totally Decentralized Application
In such category; the syntheses tasks of both the process subsystems and utility subsystem (CHP) can be conducted in a fully decentralized fashion. For example, in a typical scenario, we have two companies working on a project in which one company's objective is to synthesize the process plant (e.g., chemical complex) and the other company's objective is to synthesize a utility plant (CHP) to serve it. The utility plant synthesizer is normally a party with a certain investment capability. He or she wants to invest but with a limited budget. The process designer can reduce its capital investment at the expense of the utility system capital. However, extending this approach to its maximum capability might hinder the opportunity for finding a utility investor. During the early stages of the project, both parties discuss and negotiate the point until achieving a win-win situation.
A methodology employed by an embodiment of the present invention in order to achieve the win-win situation, systematically handles this category of application using Pareto-like curve that enables both process and utility companies to optimally find their best acceptable decision and smartly help the chemical/other industrial processes plant in its search for best financing opportunities for finding a builder to the chemical/other industrial complex utility plant. FIG. 3 illustrates utilization of a Pareto-like curve useful to aid in the determination of the optimal settling point between the conflicting goals of the two companies, ideally located at a point where movement in either direction along the curve would result in a larger increase in costs to one company than the savings costs afforded the other company.
Process and Distillation Columns Integration under Possible Combinations of Reflux Ratio and Operating Pressure Modifications Systematic Methodology.
The above described cases and/or scenarios rely on the ability to determine an optimal process subsystems' total annualized cost. In order to do so, various embodiments of the present invention provide for process and distillation columns integration under a wide range of parameters. According to an exemplary methodology, the process can include the following steps/operations:













Step 1:
Start a base case design with low column(s) operating pressure for above ambient



processes and base case design study with high columns operating pressure for sub



ambient processes.


Step 2:
Conceptually design the distillation column using low reflux ratio to high reflux



ratio.


Step3:
Register Column(s) cost in each case.


Step 4:
Simulate distillation column(s), in the pressure cases (base and higher) or (base and



lower) or low to high reflux ratio and extract for the heat integration application



streams data in form of intervals for FCp of cold stream to be evaporated (evap)



heated and hot stream to be condensed (cond) cooled, Ts and Tt are numbers



specific to the pressure used for each column.


Step 5:
Output and store [FCp_evap_P_base], ts, tt; [FCp_evap_P_new], ts, tt and



[FCp_cond_P_base], Ts, Tt; [FCp_cond_P_new], Ts, Tt.


Step 6:
Conduct energy targeting under all possible combinations of process changes and



reflux ratios for distillation column(s) at base pressure.


Step 7:
Conduct energy targeting under all possible combinations of process changes and



reflux ratios for distillation column(s) at desired high pressure (for sub-ambient



processes).


Step 8:
Find thermal pinch points interval, define different pinch points location and



identify pinch point controlling stream (this step is useful to help select which



column′ pressure need to be increased (if all of them are viable options, continue).


Step 9:
Register the difference in energy operating cost between the low pressure base



case design and high pressure column′ case and decide either to increase the



column (s) pressure or not to increase it.


Step 10:
Use energy targeting step result for distillation column(s) operating at desired high



pressure and start using for column-by-column, while other columns at high



pressure, the column′ stream data lying at low operating pressure base case design



to explore the possibility to reduce some of the columns′ pressures without having



a big impact on energy operating cost saving.


Step 11:
Define the level of integration in the plant between the distillation columns and the



rest of the process or among the distillation columns only under all possible



combinations of process design modifications for base case pressure and high



pressure using the simple algorithm illustrated in the two graphs in FIGS. 4A and 



4B which render the way to find best level of integration via conducting energy



targeting for all elements of the “SET” shown in graph.


Step 12:
Collapse the whole problem streams data from intervals to just numbers for the



desired heating and/or cooling utilities case requirement to identify utilities load



and level, best distillation column design parameters including reflux ratio,



operating pressure besides rest of process design and operating conditions



modifications.


Step 13:
Repeat the above steps for each distillation sequence in one group, for many



possible sequences, or combinations of several groups of distillation columns but



with small number of sequences, for instance as per the two graphs in FIGS. 5A and



5B.









According to Step 1, the process begins by starting a base case design with low column(s) operating pressure for above ambient processes and base case design study with high columns operating pressure for sub ambient processes. Step 2 involves conceptually designing the distillation column using a low reflux ratio to a high reflux ratio. In Step 3, the column cost is registered in each case.
Note, for above ambient temperature processes, increase the column(s) pressure from the base case design to higher pressure to decrease the column(s) cost. The highest column pressure is the one that is situated beyond the pressure where the column (cost) increases is due to its metal weight and thickness increase. For sub-ambient processes, decrease the pressure of the base case design until the increase in column (s) vapor loading/diameter increases the column(s) (cost).
In Step 4, simulate distillation column(s) according to the pressure cases (base and higher) or (base and lower), or low to high reflux ratio, and extract data for the heat integration application streams data in form of intervals for FCp of cold streams to be either evaporated (evap) heated and hot streams to be condensed (cond) cooled, where the hot and cold stream supply and target temperatures (Ts, ts, and Tt, tt, respectively) are numbers specific to the pressure used for each column. In Step 5, the following output data is recorded: [FCp_evap_P_base], ts, tt; [FCp_evap_P_new], ts, tt and [FCp_cond_P_base], Ts, Tt; [FCp_cond_P_new], Ts, Tt.
In Step 6, energy targeting is conducted under all possible combinations of process changes and reflux ratios for distillation column(s) at base pressure, and in Step 7, energy targeting is conducted, for example, under all possible combinations of process changes and reflux ratios for distillation column(s) at the desired high pressure (for above ambient processes).
In Step 8, the process determines the thermal pinch points interval, determines the different pinch points location, and identifies the pinch point controlling stream to help select which column's pressure need to be increased. If all of them are viable options, the process continues.
In Step 9, the difference in energy operating cost between the low pressure base case design and high pressure column's case registered and a decision is made to either to increase the column (s) pressure or not to increase it.
In Step 10, the results from the energy targeting step for the distillation column(s) operating at desired high pressure is utilized to perform a column-by-column examination of pressure. In this step, the column's stream data lying at low operating pressure base case design is applied to one of the columns with the other columns remaining at high pressure to explore the possibility to reduce some of the columns' pressures without having a big impact on energy operating cost savings.
In Step 11, the process defines the level of integration in the plant between the distillation columns and the rest of the process or among the distillation columns only under all possible combinations of process design modifications for base case pressure and high pressure. An example technique employs the simple algorithm illustrated in FIG. 4A, which render the roadmap to find the best level of integration via conducting energy targeting for all elements of the “SET” shown in graph. Note, the N=5 distillation column problem shown in FIG. 4A provides an example of a systematic method for “total process” columns integration under all possible combinations of reflux ratio and operating pressure modifications. Here, the four distillation columns, represented by items 2-5 in the table, are integrated with the “rest of the process” subsystem cold streams to be heated and hot streams to be cooled, represented by the item 1. The term “options” represents the total number of possible energy integration combinations of the four distillation columns with the rest of the process and each other, and “sets” represents a list of energy integration options among the four distillation columns with the rest of the process.
The N=4 distillation column problem shown in FIG. 4B provides an example of a systematic method for “distillation columns-only” columns “among-themselves-integration” under all possible combinations of reflux ratio and operating pressure modifications. Here, the four distillation columns, represented by items 1-4 in the table, are integrated with each other only. The term “options” represents the total number of possible energy integration combinations of the four distillation columns among themselves, and the “sets” represents the list of energy integration options among the four distillation columns. For example, the first element in the “set” represents the option of integrating the four columns all together. The third element in the “set” represents the option whereby columns 1 and 2 and 4 are integrated among themselves while column three is left on utility-usage only scheme.
In Step 12, the whole problem streams data is collapsed from intervals to just numbers for the desired heating and/or cooling utilities case requirement to identify utilities load and level, best distillation column design parameters including reflux ratio and operating pressure in addition to the rest of the process design and operating conditions modifications.
In Step 13, the above steps are repeated for each distillation sequence in one group for the many possible sequences, or for the combinations of several groups of distillation columns, but with small number of sequences, as shown, for example, in FIGS. 5A and 5B.
Total Steam levels Targeting Under Possible Combinations of Process Modifications
Various embodiments of the present invention provide both process steam generation interval targeting from a non-cooling task and process steam demand interval targeting from a non-heating task.
Process Steam Generation Interval Targeting from a Non-Cooling Task
According to an exemplary configuration, the in-process steam generation model for process unit (i) is calculated as follows:

[STM_GEN]ik=F([PDP],[POP])ik, 

wherein PDP=[x1,x2], and POP=[y1,y2].

PDP and POP are process design parameters and process operating parameters related/representing different process units including reflux ratio and pressure inter-heaters and inter-coolers, supply temperatures and target temperatures, catalyst type that affects streams' flowrate and/or supply and target temperatures, solvent type in absorption that affects stream's target and supply temperature, among others as understood by one of ordinary skill in the art.
Structural decisions such as complex distillation columns sequences and simple distillation sequences in multi-groups are handled via enumerating algorithms, described above. Another important application is the process steam demand for “Equipment Driving” which is handled in CHP algorithm, described later.
Process Steam Demand Interval Targeting from a Non-Heating Task
According to an exemplary configuration, the in-process steam demand model for process unit (j) is calculated as follows:

[Process_STM_Demand]jk=F([PDP],[POP])jk, 

wherein PDP=[D1, D2], and POP=[O1, O2].

Heating Steam Levels Targeting Under Possible Combinations of Process Modifications
According to an exemplary configuration, heating steam levels targeting utilizes the algorithms described, for example, in U.S. patent application Ser. No. 12/480,415, titled “System, Program Product, and Related Methods for Global Targeting of Process Utilities Under Varying Conditions.” According to the exemplary configuration, the process includes the various steps shown in the table illustrated below and described thereafter.













Step 1:
List hot and cold stream data in interval (range) form for the



problem to be solved in a data entry table.


Step 2:
Construct and interval-cascade diagram using the entered



problem data.


Step 3:
Balance the interval cascade diagram via adding a hot stream.


Step 4:
Draw the residuals from each step versus the average



temperature on a T-Q diagram to produce the



Interval-Grande-Composite-Curves.


Step 5:
Combine the low-pressure, medium pressure, and high pressure



steam intervals rendered by the Interval-Grand-Composite



Curves to those produced from other process applications



due to non-heating process applications to identify the total



steam levels under all possible combinations of process



modifications.









As shown in FIG. 6, the steps include listing the stream data in interval form for the problem to be solved in a data entry table 61. An interval-cascade diagram 71 as shown, for example, in FIG. 7 is then constructed using the problem data entered in the entry table 61. As illustrated in FIGS. 8A and 8B, the steps can include balancing the interval cascade diagram 71 (FIG. 8A) via adding a hot utility at 81 and drawing the residuals from each step versus the shown temperature (which is the average temperature of such step in the temperature interval diagram) on a T-Q diagram 83 (FIG. 8B) to produce the Interval-Grande-Composite-Curves shown therein.
Beneficially, the curves include all possible grand composite curves/process steam heating demands and steam generation due to process cooling tasks. Correspondingly, the diagram 83 renders the low-pressure steam interval [LPS], medium pressure steam interval [MPS], and high pressure steam interval [HPS] under all possible combinations of process design modifications.
As shown in FIG. 9, the steps further include combining the low-pressure steam interval [LPS], medium pressure steam interval [MPS], and high pressure steam interval [HPS] rendered by the Interval-Grand-Composite Curves to those produced from other process applications due to non-heating process applications to identify the total steam level interval values required by the problem.
Number of Steam Headers Loads Identification Algorithm
According to an exemplary configuration, steam headers loads identification utilizes the CHP algorithm to calculate the utility system total annualized cost (USTAC) for 2, 3, 4, 5, 6 and 7 headers under all possible combinations of process design changes. For example, we find the USTAC for 2 headers (FIG. 10A) and then 3 headers (FIG. 10B) to identify the most economical number of steam headers.
As shown in FIG. 11, according to a preferred configuration, the process begins with two headers (block 101), employs the total stream generation and demand (level) 2+i lower and upper bounds calculation algorithm (block 103), determines the USTAC through use of the CHP program/algorithm (block 105), records the annualized cost data for the respective number of headers (block 107) and continues a loop, increasing the number of headers one-by-one (block 109) according to CHP calculation results until the increase in the number of steam headers no longer decreases the total cost (block 110).
Number of Steam Headers and Operating Conditions Enumeration Algorithm
FIG. 12 provides an illustrative example for the number of steam headers and operating conditions enumeration algorithm. The algorithm identifies the number of cases based on the number of steam headers selected by the user in conjunction with the number of operating conditions (i.e., hot, medium, and low steam pressures).

Number of Cases=Number of operating conditions^(number of headers−1).

For example, for a high steam temperature of header H1=700° F. and number of headers N=3, the Number of Cases=3^(3−1)=9.
Simultaneous Synthesis of Combined Heat and Power (CHP) Utility and Process Subsystems Through Employment of a Combined Heat and Power Mathematical Program
According to various embodiments of the present invention, the simultaneous CHP and key process units syntheses alternatives can be broadly classified into two categories: process related options and utility related options. The process related options can include altering the amount of steam and power demand by the process and steam and power generated by the process due to the process' separation and HEN systems design. Note, according to the preferred configuration, the process related options do not include the process chemical reaction path, catalysts and reactors network as options for change for the sake of optimal utility; HEN; and separation systems' design. The utility/CHP related options normally concentrate on the optimal generation of heat and power (CHP) from the synthesized system to satisfy the process needs.
According to an exemplary configuration, the CHP mathematical program is a mixed integer linear program (MILP) employed to provide simultaneous synthesis of CHP utility and process subsystem. Although there exists commercial software available in the market that employs MILP, and although many of the academic models in public domain also use either MILP or mixed integer non-linear programs (MINLPs) due to nonlinearity in its model objective function and/or model constraints, various embodiments of the present invention beneficially employ a MILP in which the process steam demand and process steam generation is defined at different pressure levels, and in which process steam demand, process steam generation, power generation, and power export are defined as decision variables with rigorous range-type boundaries.
Further, beneficially, the optimal values for both process and utility subsystems syntheses are left to be decided by the CHP MILP model and not as input variables—which is the case in all other CHP MILP and MINLP in academia and industry, where such variables are specified/given ahead of the MILP solution. In contrast to the various embodiments of the present invention, such input data are conventionally normally fed to the utility system designer by the chemical process designer.
FIGS. 13-15 illustrates select examples of input data graphical user interface fields illustrating how the input data is loaded according to an exemplary configuration. The program input data includes major program sets, utility and process parameters, and utility and process variables.
Appendix 3 lists the various examples of the program sets, input parameters, and input variables according to the exemplary configuration. The major program sets can include: “original sets” and “derived sets.” The parameters can include: “steam headers” parameters, “boilers” parameters, “cogeneration” parameters, “steam turbine generator” parameters, condensing steam turbine generator parameters, breakover station parameters, condenser parameters, process steam turbine parameters, switchable motors to process steam turbines parameters, process steam generator parameters, solar thermal unit parameters, condensate system parameters, deaerator parameters, power parameters, fuel parameters, and economics parameters. The variables can include steam headers variables, boilers variables, cogeneration variables, solar steam generator variables, steam turbine generator variables, condensing steam turbine generator variables, switchable motors to process steam turbines variables, break over station variables, condenser variables, condensate system variables, deaerator variables, process steam generator variables, boiler feed water system variables, and process cooling demand variables.
Select calculations directed to steam headers, cogeneration, boilers, steam turbine generators, condensing steam turbine generators, solar thermal steam generation/heating, steam demands, breakover stations, condenser, skim tank, deaerators, boiler feed water, power, constraints, fuel, steam reserves, and economics, according to the exemplary configuration, are provided below:
Steam Headers:
Header Inlet Flow.



  
    
      HDRIn
      ⁡
      
        [
        h
        ]
      
    
    =
    
      
        HDRRet
        ⁡
        
          [
          h
          ]
        
      
      +
      
        HDRBFW
        ⁡
        
          [
          h
          ]
        
      
      +
      
        
          ∑
          
            
              (
              
                b
                ,
                h
              
              )
            
            ∈
            BLRHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            BLRSTM
            ⁡
            
              [
              b
              ]
            
          
          ·
          
            BLRNo
            ⁡
            
              [
              b
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                so
                ,
                h
              
              )
            
            ∈
            SOLARHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            SOLARSTM
            ⁡
            
              [
              so
              ]
            
          
          ·
          
            SOLARNo
            ⁡
            
              [
              so
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                c
                ,
                h
              
              )
            
            ∈
            COGENHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            COGENStm
            ⁡
            
              [
              c
              ]
            
          
          ·
          
            COGENNo
            ⁡
            
              [
              c
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                s
                ,
                h
              
              )
            
            ∈
            STGOutHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            STGStm
            ⁡
            
              [
              s
              ]
            
          
          ·
          
            (
            
              1
              -
              
                STGLoss
                ⁡
                
                  [
                  s
                  ]
                
              
            
            )
          
          ·
          
            STGNo
            ⁡
            
              [
              s
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                b
                ,
                h
              
              )
            
            ∈
            BOSOut
          
        
        ⁢
        
          
        
        ⁢
        
          
            BOSStm
            ⁡
            
              [
              b
              ]
            
          
          ·
          
            (
            
              1
              -
              
                BOSLoss
                ⁡
                
                  [
                  b
                  ]
                
              
            
            )
          
          ·
          
            BOSNo
            ⁡
            
              [
              b
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                p
                ,
                h
              
              )
            
            ∈
            PSTOutHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            PSTStm
            ⁡
            
              [
              p
              ]
            
          
          ·
          
            (
            
              1
              -
              
                PSTLoss
                ⁡
                
                  [
                  p
                  ]
                
              
            
            )
          
        
      
      +
      
        
          ∑
          
            
              (
              
                p
                ,
                h
              
              )
            
            ∈
            PSGenHDR
          
        
        ⁢
        
          
        
        ⁢
        
          PSGenProd
          ⁡
          
            [
            p
            ]
          
        
      
    
  




  
    
      
    
    ⁢
    
      ∀
      
        h
        ∈
        HDR
      
    
  



Header Outlet Flow.



  
    
      HDROut
      ⁡
      
        [
        h
        ]
      
    
    =
    
      
        
          ∑
          
            
              (
              
                s
                ,
                h
              
              )
            
            ∈
            SDemHDR
          
        
        ⁢
        
          
        
        ⁢
        
          HDRDem
          ⁡
          
            [
            h
            ]
          
        
      
      +
      
        
          ∑
          
            
              (
              
                s
                ,
                h
              
              )
            
            ∈
            STGInHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            STGStm
            ⁡
            
              [
              s
              ]
            
          
          ·
          
            STGNo
            ⁡
            
              [
              s
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                s
                ,
                h
              
              )
            
            ∈
            STGInHDR
          
        
        ⁢
        
          
        
        ⁢
        
          PSTSwitchStm
          ·
          
            PSTNo
            ⁡
            
              [
              s
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                c
                ,
                h
              
              )
            
            ∈
            CSTGHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            CSTGStm
            ⁡
            
              [
              c
              ]
            
          
          ·
          
            CSTGNo
            ⁡
            
              [
              c
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                s
                ,
                h
              
              )
            
            ∈
            PSTInHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            PSTStm
            ⁡
            
              [
              s
              ]
            
          
          ·
          
            PSTNo
            ⁡
            
              [
              s
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                b
                ,
                h
              
              )
            
            ∈
            BOSIn
          
        
        ⁢
        
          
        
        ⁢
        
          
            BOSStm
            ⁡
            
              [
              b
              ]
            
          
          ·
          
            BOSNo
            ⁡
            
              [
              b
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                c
                ,
                h
              
              )
            
            ∈
            CONDHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            CONDStm
            ⁡
            
              [
              c
              ]
            
          
          ·
          
            CONDNo
            ⁡
            
              [
              c
              ]
            
          
        
      
      +
      
        
          ∑
          
            h
            ∈
            DEAStmHDR
          
        
        ⁢
        
          
        
        ⁢
        DEAStmFlow
      
      +
      
        
          ∑
          
            
              (
              
                p
                ,
                h
              
              )
            
            ∈
            PSTInHDR
          
        
        ⁢
        
          
        
        ⁢
        
          PSTDem
          ⁡
          
            [
            p
            ]
          
        
      
    
  




  
    
      
    
    ⁢
    
      ∀
      
        h
        ∈
        HDR
      
    
  



Material Balance.

HDRIn[h]·(1−HDRLoss[h])=HDROut[h]∀hεHDR

Energy Balance.



  
    
      
        HDRIn
        ⁡
        
          [
          h
          ]
        
      
      ·
      
        HDREnth
        ⁡
        
          [
          h
          ]
        
      
    
    =
    
      
        
          HDRRet
          ⁡
          
            [
            h
            ]
          
        
        ·
        
          HDREnth
          ⁡
          
            [
            h
            ]
          
        
      
      +
      
        
          HDRBFW
          ⁡
          
            [
            h
            ]
          
        
        ·
        
          BFWEnth
          ⁡
          
            [
            h
            ]
          
        
      
      +
      
        
          HDREnth
          ⁡
          
            [
            h
            ]
          
        
        ·
        
          
            ∑
            
              
                (
                
                  so
                  ,
                  h
                
                )
              
              ∈
              SOLARHDR
            
          
          ⁢
          
            
          
          ⁢
          
            
              SOLARSTM
              ⁡
              
                [
                so
                ]
              
            
            ·
            
              SOLARNo
              ⁡
              
                [
                so
                ]
              
            
          
        
      
      +
      
        
          HDREnth
          ⁡
          
            [
            h
            ]
          
        
        ·
        
          
            ∑
            
              
                (
                
                  b
                  ,
                  h
                
                )
              
              ∈
              BLRHDR
            
          
          ⁢
          
            
          
          ⁢
          
            
              BLRSTM
              ⁡
              
                [
                b
                ]
              
            
            ·
            
              BLRNo
              ⁡
              
                [
                b
                ]
              
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                s
                ,
                h
              
              )
            
            ∈
            STGInHDR
          
        
        ⁢
        
          
        
        ⁢
        
          PSTSwitchStm
          ·
          
            PSTNo
            ⁡
            
              [
              s
              ]
            
          
          ·
          PSTEnthOut
        
      
      +
      
        
          
            HDREnth
            ⁡
            
              [
              h
              ]
            
          
          ·
          
            
              
                
              
              ∑
            
            
              
                (
                
                  c
                  ,
                  h
                
                )
              
              ∈
              COGENHDR
            
          
        
        ⁢
        
          
        
        ⁢
        
          
            COGENStm
            ⁡
            
              [
              c
              ]
            
          
          ·
          
            COGENNo
            ⁡
            
              [
              c
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                s
                ,
                h
              
              )
            
            ∈
            STGOutHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            STGStm
            ⁡
            
              [
              s
              ]
            
          
          ·
          
            (
            
              1
              -
              
                STGLoss
                ⁡
                
                  [
                  s
                  ]
                
              
            
            )
          
          ·
          
            STGNo
            ⁡
            
              [
              s
              ]
            
          
          ·
          
            STGEnthOut
            ⁡
            
              [
              s
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                b
                ,
                h
              
              )
            
            ∈
            BOSOut
          
        
        ⁢
        
          
        
        ⁢
        
          
            BOSStm
            ⁡
            
              [
              b
              ]
            
          
          ·
          
            (
            
              1
              -
              
                BOSLoss
                ⁡
                
                  [
                  b
                  ]
                
              
            
            )
          
          ·
          
            BOSNo
            ⁡
            
              [
              b
              ]
            
          
          ·
          
            BOSEnth
            ⁡
            
              [
              b
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                p
                ,
                h
              
              )
            
            ∈
            PSTOutHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            PSTDem
            ⁡
            
              [
              p
              ]
            
          
          ·
          
            (
            
              1
              -
              
                PSTLoss
                ⁡
                
                  [
                  p
                  ]
                
              
            
            )
          
          ·
          
            PSSTEnth
            ⁡
            
              [
              p
              ]
            
          
        
      
      +
      
        
          ∑
          
            
              (
              
                p
                ,
                h
              
              )
            
            ∈
            PSGenHDR
          
        
        ⁢
        
          
        
        ⁢
        
          
            PSGenProd
            ⁡
            
              [
              p
              ]
            
          
          ·
          
            PSGenEnth
            ⁡
            
              [
              p
              ]
            
          
        
      
    
  




  
    
      
    
    ⁢
    
      ∀
      
        h
        ∈
        HDR
      
    
  



Cogeneration:
Steam Production.

COGENStm[c]=COGENPwr[c]·COGENStmRatio[c]∀cεCOGEN

Material Balance.

COGENBFW[c]·(1−COGENLoss)=COGENStm[c]+COGENBD[c]∀cεCOGEN

Blow Down Calculation.

COGENBD[c]·(COGENCC[c]−1)=COGENStm[c]∀cεCOGEN

Cogeneration Fuel Use.




  
    COGENFuel
    ⁡
    
      [
      c
      ]
    
  
  =
  
    {
    
      
        
          
            
              
                
                  
                    
                      COGENc
                      ⁢
                      
                        
                      
                      ⁢
                      
                        0
                        ⁡
                        
                          [
                          c
                          ]
                        
                      
                    
                    +
                    
                      COGENc
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                          1
                          ⁡
                          
                            [
                            c
                            ]
                          
                        
                        ·
                        
                          COGENPwr
                          ⁡
                          
                            [
                            c
                            ]
                          
                        
                      
                    
                    +
                  
                
              
              
                
                  
                    
                      COGENc
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                          2
                          ⁡
                          
                            [
                            c
                            ]
                          
                        
                        ·
                        
                          
                            COGENPwr
                            ⁡
                            
                              [
                              c
                              ]
                            
                          
                          2
                        
                      
                    
                    ,
                  
                
              
            
          
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                
                  COGENStm
                  ⁡
                  
                    [
                    c
                    ]
                  
                
              
              >
              0
            
          
        
        
          
            
              0
              ,
            
          
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                
                  COGENStm
                  ⁡
                  
                    [
                    c
                    ]
                  
                
              
              =
              0
            
          
        
      
      ⁢
      
        

      
      ⁢
      
        
      
      ⁢
      
        ∀
        
          c
          ∈
          COGEN
        
      
    
  




Capacity Constraints.

COGENPwr[c]≧COGENMin[c]·COGENOn[c]∀cεCOGEN

COGENPwr[c]≦COGENCap[c]·COGENOn[c]∀cεCOGEN

Boilers:
Material Balance.

BLRBFW[b]·(1−BLRLoss)=BLRStm[b]+BLRBD[b]∀bεBLR

Boiler Blowdown Calculation.

BLRBD[b]·(BLRCC[b]−1)=BLRStm[b]∀bεBLR

Boiler Fuel Use Calculation.

BLRFuel[b]=BLRc0[b]+BLRc1[b]·BLRStm[b]+BLRc2·BLRStm[b]2∀bεBLR

Boiler Capacity Constraints.

BLRStm[b]≧BLRMin[b]·BLROn[b]∀bεBLR

BLRStm[b]≦BLRCap[b]·BLROn[b]∀bεBLR

Boiler Steam Reserves.




  
    ∑
    
      b
      ∈
      BLR
    
  
  ⁢
  
    
  
  ⁢
  
    
      BLROn
      ⁡
      
        [
        b
        ]
      
    
    ·
    
      (
      
        
          BLRCap
          ⁡
          
            [
            b
            ]
          
        
        -
        
          BLRStm
          ⁡
          
            [
            b
            ]
          
        
      
      )
    
    ·
    
      BLRNo
      ⁡
      
        [
        b
        ]
      
    
  




Steam Turbine Generator:
Power Generation Calculation.

STGPwr[s]·3412=STGStm[s]·(STGEnthIn[s]−STGEnthOut[s])∀sεSTG

Capacity.

STGStm[s]≧0∀sεSTG

STGStm[s]≦STGCap[s]∀sεSTG

Condensing Steam Turbine Generator:
Power Generation Calculation.

CSTGPwr[c]·3412=CSTGStm[c]·(CSTGEnthIn[c]−CSTGEnthOut[c])∀cεCSTG

Capacity.

CSTGStm[c]≧0∀cεCSTG

CSTGStm[c]≦CSTGCap[c]∀cεCSTG

Solar Thermal Steam Generation/Heating:
Material Balance:

SOLARBFW[so]·(1−SOLARLoss)=SOLARStm[so]∀soεSOLAR

Steam Production:

SOLARStm[SO]=SOLARon[so]·SOLAREff[so]·SOLARRad∀soεSOLAR

Steam Demands:
Demand Ranges in a Header:

PSDem[h]≧HDRMinDem∀hεHDR 

PSDem[h]≦HDRMaxDem∀hεHDR 

Break Over Stations:
Capacity.

BOSStm[b]≧0∀bεBOS

BOStm[b]≦BOSCap[b]∀bεBOS

Condenser:
Duty Calculation.

CONDDuty[c]·1000=CONDStm[c]·CONDEnthDuty[c]∀cεCOND

Capacity.

CONDStm[c]≧0∀cεCOND

CONDStm[c]≦CONDCap[c]∀cεCOND

Skim Tank:
Material Balance.



  
    
      SkimTankInFlow
      ·
      
        (
        
          1
          -
          SkimTankLoss
        
        )
      
    
    =
    SkimTankOutFlow
  




  
    SkimTankInFlow
    =
    
      CondRetFlow
      +
      MUWtrFlow
      +
      
        
          ∑
          
            c
            ∈
            CSTG
          
        
        ⁢
        
          
        
        ⁢
        
          
            CSTGStm
            ⁡
            
              [
              c
              ]
            
          
          ·
          
            (
            
              1
              -
              
                CSTGLoss
                ⁡
                
                  [
                  c
                  ]
                
              
            
            )
          
          ·
          
            CSTGNo
            ⁡
            
              [
              c
              ]
            
          
        
      
      +
      
        
          ∑
          
            c
            ∈
            COND
          
        
        ⁢
        
          
        
        ⁢
        
          
            CONDStm
            ⁡
            
              [
              c
              ]
            
          
          ·
          
            (
            
              1
              -
              
                CONDLoss
                ⁡
                
                  [
                  c
                  ]
                
              
            
            )
          
          ·
          
            CONDNo
            ⁡
            
              [
              c
              ]
            
          
        
      
    
  



Energy Balance.




  
    SkimTankOutFlow
    ·
    SkimTankOutEnth
  
  =
  
    
      CondRetFlow
      ·
      CondRetTargEnth
    
    +
    
      MUWtrFlow
      ·
      MUWtrEnth
    
    +
    
      
        ∑
        
          c
          ∈
          CSTG
        
      
      ⁢
      
        
          CSTGStm
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          (
          
            1
            -
            
              CSTGLoss
              ⁡
              
                [
                c
                ]
              
            
          
          )
        
        ·
        
          CSTGWtrEnth
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          CSTGNo
          ⁡
          
            [
            c
            ]
          
        
      
    
    +
    
      
        ∑
        
          c
          ∈
          COND
        
      
      ⁢
      
        
          CONDStm
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          (
          
            1
            -
            
              CONDLoss
              ⁡
              
                [
                c
                ]
              
            
          
          )
        
        ·
        
          CONDEnth
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          CONDNo
          ⁡
          
            [
            c
            ]
          
        
      
    
  




Deaerator:
Material Balance.

(SkimTankOutFlow+DEAStmFlow)·(1−DEALoss)=DEAVentFlow+DEAOutWtrFlow

Energy Balance.

SkimTankOutFlow·SkimTankOutEnth+DEAStmFlow·DEAStmEnth=DEAVentFlow·DEAVentEnth+DEAOutWtrFlow·DEAOutWtrEnth

Boiler Feed Water:
Material Balance.



  
    BFWFow
    =
    
      
        
          ∑
          
            h
            ∈
            HDR
          
        
        ⁢
        
          
        
        ⁢
        
          HDRBFW
          ⁡
          
            [
            h
            ]
          
        
      
      +
      
        
          ∑
          
            b
            ∈
            BLR
          
        
        ⁢
        
          
        
        ⁢
        
          
            BLRBFW
            ⁡
            
              [
              b
              ]
            
          
          ·
          
            BLRNo
            ⁡
            
              [
              b
              ]
            
          
        
      
      +
      
        
          ∑
          
            c
            ∈
            COGEN
          
        
        ⁢
        
          
        
        ⁢
        
          
            COGEN
            ⁡
            
              [
              c
              ]
            
          
          ·
          
            COGENNo
            ⁡
            
              [
              c
              ]
            
          
        
      
    
  




  
    
      
    
    ⁢
    
      BFWFow
      =
      
        DEAOutWtrFlow
        ·
        
          (
          
            1
            -
            BFWLoss
          
          )
        
      
    
  



Power:
Total Power Generation Calculation.




  TotPwrGen
  =
  
    
      
        ∑
        
          c
          ∈
          COGEN
        
      
      ⁢
      
        
      
      ⁢
      
        
          COGENPwr
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          COGENNo
          ⁡
          
            [
            c
            ]
          
        
      
    
    +
    
      
        ∑
        
          s
          ∈
          STG
        
      
      ⁢
      
        
      
      ⁢
      
        
          STGPwr
          ⁡
          
            [
            s
            ]
          
        
        ·
        
          STGNo
          ⁡
          
            [
            s
            ]
          
        
      
    
    +
    
      
        ∑
        
          ms
          ∈
          STG
        
      
      ⁢
      
        
      
      ⁢
      MSwitchMWToPST
    
    +
    
      
        ∑
        
          c
          ∈
          CSTG
        
      
      ⁢
      
        
      
      ⁢
      
        
          CSTGPwr
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          CSTGNo
          ⁡
          
            [
            c
            ]
          
        
      
    
  




Satisfy Power Requirements.

TotPwrGen+TotPwrImp−TotPwrExp≧ProcPwrDem

Exported Power Calculation.




  TotPwrExp
  =
  
    {
    
      
        
          
            
              
                
                  TotPwrGen
                  +
                  TotPwrImp
                  -
                
              
            
            
              
                
                  ProcPwrDem
                  ,
                
              
            
          
        
        
          
            
              
                
                  
                    
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      TotPwrGen
                    
                    +
                    TotPwrImp
                  
                  >
                
              
            
            
              
                ProcPwrDem
              
            
          
        
      
      
        
          
            0
            ,
          
        
        
          Otherwise
        
      
    
  




Constraints:
Process Cooling Demand:




  PCoolingDem
  =
  
    
      
        ∑
        
          
            (
            
              i
              ,
              PC
            
            )
          
          ∈
          PCoolingDem
        
      
      ⁢
      
        
      
      ⁢
      
        AirCDuty
        ⁡
        
          [
          i
          ]
        
      
    
    +
    
      
        ∑
        
          
            (
            
              i
              ,
              PC
            
            )
          
          ∈
          PCoolingDem
        
      
      ⁢
      
        
      
      ⁢
      
        WtrCDuty
        ⁡
        
          [
          i
          ]
        
      
    
    +
    
      
        ∑
        
          
            (
            
              i
              ,
              PC
            
            )
          
          ∈
          PCoolingDem
        
      
      ⁢
      
        
      
      ⁢
      
        RefDutyR
        ⁡
        
          [
          i
          ]
        
      
    
  




Fuel:
Total Fuel Use.




  TotFuelUse
  =
  
    
      
        ∑
        
          b
          ∈
          BLR
        
      
      ⁢
      
        
      
      ⁢
      
        
          BLRFuel
          ⁡
          
            [
            b
            ]
          
        
        ·
        
          BLRNo
          ⁡
          
            [
            b
            ]
          
        
      
    
    +
    
      
        ∑
        
          c
          ∈
          COGEN
        
      
      ⁢
      
        
      
      ⁢
      
        
          COGENFuel
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          COGENNo
          ⁡
          
            [
            c
            ]
          
        
      
    
  




Steam Reserves:
Maximum Steam Generation from a Boiler Unit.




  MaxBlrStm
  =
  
    
      max
      
        b
        ∈
        BLR
      
    
    ⁢
    
      {
      
        
          BLRStm
          ⁡
          
            [
            b
            ]
          
        
        ·
        
          BLRNo
          ⁡
          
            [
            b
            ]
          
        
      
      }
    
  




Maximum Steam Generation from a Cogeneration Unit.




  MaxCogenStm
  =
  
    
      max
      
        c
        ∈
        COGEN
      
    
    ⁢
    
      {
      
        
          COGENStm
          ⁡
          
            [
            b
            ]
          
        
        ·
        
          COGENNo
          ⁡
          
            [
            b
            ]
          
        
      
      }
    
  




Required Steam Reserves Calculation.

RequiredStmRes=max{MaxBlrStm,MaxCogenStm}

Condensate Steam Reserves Calculation.




  CondStmRes
  =
  
    
      
        ∑
        
          c
          ∈
          CSTG
        
      
      ⁢
      
        
      
      ⁢
      
        
          CSTGStm
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          (
          
            1
            -
            
              CSTGLoss
              ⁡
              
                [
                c
                ]
              
            
          
          )
        
        ·
        
          CSTGNo
          ⁡
          
            [
            c
            ]
          
        
      
    
    +
    
      
        ∑
        
          c
          ∈
          COND
        
      
      ⁢
      
        
      
      ⁢
      
        
          COND
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          (
          
            1
            -
            
              CONDLoss
              ⁡
              
                [
                c
                ]
              
            
          
          )
        
        ·
        
          CONDNo
          ⁡
          
            [
            c
            ]
          
        
      
    
  




Steam Reserves Calculation.

StmRes=max{RequiredStmRes−CondStmRes,0}

Economics:
Total Cost Calculation.




  TotCost
  =
  
    
      
        (
        
          
            TotFuelUse
            ·
            PriceFuel
          
          +
          
            MUWtrFlow
            ·
            PriceMUWtr
          
          +
          
            TotPwrImp
            ·
            PriceImpPwr
          
          -
          
            TotPwrExp
            ·
            PriceExpPwr
          
        
        )
      
      /
      1000
    
    +
    
      
        ∑
        
          b
          ∈
          BLR
        
      
      ⁢
      
        
      
      ⁢
      
        
          BLRNo
          ⁡
          
            [
            b
            ]
          
        
        ·
        
          BLRPrice
          ⁡
          
            [
            b
            ]
          
        
      
    
    +
    
      
        ∑
        
          h
          ∈
          hdr
        
      
      ⁢
      
        
      
      ⁢
      
        
          Hdr
          ⁡
          
            [
            h
            ]
          
        
        ·
        
          HdrPrice
          ⁡
          
            [
            h
            ]
          
        
      
    
    +
    
      
        ∑
        
          c
          ∈
          COGEN
        
      
      ⁢
      
        
      
      ⁢
      
        
          COGENNo
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          COGENPrice
          ⁡
          
            [
            c
            ]
          
        
      
    
    +
    
      
        ∑
        
          sO
          ∈
          SOLAR
        
      
      ⁢
      
        
      
      ⁢
      
        
          SOLARNo
          ⁡
          
            [
            s
            ]
          
        
        ·
        
          SOLARPrice
          ⁡
          
            [
            s
            ]
          
        
      
    
    +
    
      
        ∑
        
          s
          ∈
          STG
        
      
      ⁢
      
        
      
      ⁢
      
        
          STGNo
          ⁡
          
            [
            s
            ]
          
        
        ·
        
          STGPrice
          ⁡
          
            [
            s
            ]
          
        
      
    
    +
    
      
        ∑
        
          c
          ∈
          CSTG
        
      
      ⁢
      
        
      
      ⁢
      
        
          CSTGNo
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          CSTGPrice
          ⁡
          
            [
            c
            ]
          
        
      
    
    +
    
      
        ∑
        
          b
          ∈
          BOS
        
      
      ⁢
      
        
      
      ⁢
      
        
          BOSNo
          ⁡
          
            [
            b
            ]
          
        
        ·
        
          BOSPrice
          ⁡
          
            [
            b
            ]
          
        
      
    
    +
    
      
        ∑
        
          c
          ∈
          COND
        
      
      ⁢
      
        
      
      ⁢
      
        
          CONDNo
          ⁡
          
            [
            c
            ]
          
        
        ·
        
          CONDPrice
          ⁡
          
            [
            c
            ]
          
        
      
    
    +
    
      
        ∑
        
          
            (
            
              i
              ,
              PC
            
            )
          
          ∈
          PCoolingDem
        
      
      ⁢
      
        
      
      ⁢
      
        
          AirCDuty
          ⁡
          
            [
            i
            ]
          
        
        ·
        
          AirCPrice
          ⁡
          
            [
            i
            ]
          
        
      
    
    +
    
      
        ∑
        
          
            (
            
              i
              ,
              PC
            
            )
          
          ∈
          PCoolingDem
        
      
      ⁢
      
        
      
      ⁢
      
        
          WtrCDuty
          ⁡
          
            [
            i
            ]
          
        
        ·
        
          WtrCPrice
          ⁡
          
            [
            i
            ]
          
        
      
    
    +
    
      
        ∑
        
          
            (
            
              i
              ,
              PC
            
            )
          
          ∈
          PCoolingDem
        
      
      ⁢
      
        
      
      ⁢
      
        
          RefDutyR
          ⁡
          
            [
            i
            ]
          
        
        ·
        
          RefCPrice
          ⁡
          
            [
            i
            ]
          
        
      
    
  




Results provided by the CHP mathematical program output module, according to the exemplary configuration, include: the number of steam headers and operating conditions, the number of cogeneration units and its/their designated capacities, the number of boilers and its/their designated capacities, the number of steam turbine generators and its/their capacities, the number of electrical motors and steam turbines for process equipment driving and its/their size, the size and location of solar system, the CHP utility system steam and power generation and allocation, and the process plant process steam generation and demand loads and levels. FIG. 16 illustrates an output model diagram including the CHP utility system resulting from the entered input parameters/variables identified in FIGS. 13-15. FIG. 17 provides the respective abbreviations utilized to identify the various components in the output model.
Optimal Process Conditions and HEN Minimum Number of Units Identification.
As part of the above described methodology and mathematical model program, according to the exemplary configuration, the mathematical program can include a module for identifying the best pinch point resulting in or associated with the minimum number of heat exchange unit (HEN) units. The module includes program code which incorporates program code provided for/under all possible combinations of process changes and optimal driving force distribution in energy recovery systems described, for example, in U.S. Pat. No. 7,698,022. The steps/operations performed by the module, according to the exemplary configuration, are as follows:













Step 1:
Perform energy utility targeting under all possible combinations



of different process conditions using stream-specific minimum



approach temperatures (ΔTmini).


Step 2:
Develop a temperature interval diagram for each possible pinch



point temperature to enable the calculations of the associated



number of HEN units.


Step 3:
For each identified pinch point, estimate the minimum number



of HEN units for both the region above and the region below the



respective pinch point using a pinch heuristic (Nu = S − 1).


Step 4:
Compare the results and define the pinch point temperature at



which minimum number of units can be obtained.


Step 5:
Use an interval collapsing algorithm to find the discrete process



conditions rendering the desired pinch point temperature that



results in the targeted minimum number of HEN units.









Step 1: This step includes performing an energy utility targeting operation under all possible combinations of different process conditions using, for example, stream-specific minimum approach temperatures (ΔTmini). If the targeting step renders a pinch interval, rather than a single pinch point, then there is potential to explore best pinch point for minimum number of HENs. The Pinch points set includes all the streams possible start temperatures included in this set.
Step 2: This step includes developing a temperature interval diagram for each possible pinch point temperature to enable the calculations of the associated number of HEN units.
Step 3: This step includes estimating the minimum number of HEN units for both the region above and the region below the respective pinch point of each identified pinch point using a pinch heuristic: Nu=S−1, where Nu=Number of HEN units, and S=Number of process streams plus Number of utility levels.
Step 4: Compare the results and define the pinch point temperature at which the minimum number of units can be obtained.
Step 5: This step includes utilizing the interval collapsing algorithm described, for example, in U.S. Pat. No. 7,698,022, to find the discrete process conditions rendering the desired pinch point temperature that results in the targeted minimum number of HEN units.
FIGS. 18-21 provide two examples of interval diagrams showing the location of the pinch point and corresponding process conditions leading to a reduced number of HEN units.
FIG. 18 illustrates a first example which includes four hot streams H1-H4 and two cold streams C1-C2. In this example, hot stream H3 and cold stream C1 cross the pinch line located at 131. FIG. 19 illustrates how utilization of a pinch point located at 132 renders process conditions which provide for utilization of two less heat exchanger units. FIG. 20 similarly illustrates a second example, which includes two hot streams H1-H2 and two cold streams C1-C2 at their respective supply and target temperature process conditions. In this example, hot stream H1 and cold streams C1-C2 cross the pinch line located at 141. FIG. 21 illustrates how utilization of a pinch point located at 142 renders process conditions which provide for utilization of two less heat exchanger units.
Optimal Process Conditions and HEN Area Method and Mathematical Model Details (“LINGO” Solver Environment)
As part of the above described method and mathematical model program, the exemplary configuration includes steps for determining optimal process conditions, and the optimal heat exchanger network (HEN) surface area. The steps/operations performed by the module, according to the exemplary configuration, are as follows:













Step 0:
Select new project or access previously saved project.


Step 1:
Define heating and cooling utilities lower and upper bounds



utilizing an energy utilities targeting methodology.


Step 2:
Develop possible configurations for heat exchangers network



under possible combinations of process conditions for a



designated minimum temperature approach.


Step 3:
Develop a non-linear program (NLP) model in which each



stream supply temperature, target temperature and heat capacity



flowrate (FCp), as well as the HEN heating and cooling



duties are rigorously bounded variables.


Step 4:
Solve the NLP to identify best process conditions associated



with minimum total heat exchangers network surface area



(i.e., the surface area of each heat exchanger in the network).


Step 5:
Repeat steps 3 and 4 for each possible HEN configuration.


Step 6:
Repeat procedures 3, 4 and 5 for each of one or more but more



typically a plurality of desired or possible {ΔT_min_i} sets.









Step 0: This step includes accessing a “listing projects” webpage (see, e.g., FIG. 22). As shown in FIG. 23, selection of existing projects includes accessing project data sets. Selecting either an existing project or a new project provides access to edit/originate the project details (see, e.g., FIGS. 13-15) and show the project map.
Step 1: This step includes defining the heating and cooling utilities lower and upper bounds for the process utilizing an energy utilities targeting methodology such as, for example, that described in U.S. Pat. No. 7,698,022.
Step 2: This step includes developing possible configurations for the heat exchangers network under possible/probable combinations of process conditions for designated minimum temperature approach (e.g. {ΔT_min_i}=10) utilizing, for example, the methodology described, for example, in U.S. patent application Ser. Nos. 12/898,461 and/or 12/898,484.
Step 3: This step includes developing a non-linear program (NLP) model in which each stream supply temperature, target temperature and heat capacity flowrate (FCp), as well as the HEN heating and cooling duties are rigorously bounded variables. An exemplary nonlinear program model is provided, for example, in Appendix 3 and in FIGS. 22-40. FIGS. 24A-B together illustrate a portion of the program code for entering the network topology. FIG. 25 illustrates a network schematic example. FIG. 26 illustrates an “enter data” snapshot. FIG. 27, for example, illustrates a model formulation view.
Step 4: This step includes solving the NLP to identify best process conditions associated with minimum total heat exchangers network surface area (i.e., the surface area of each heat exchanger in the network). FIG. 28 illustrates a solution results view.
Step 5: This step includes repeating steps 3 and 4 for each possible HEN configuration.
Step 6: This step includes repeating steps 3, 4 and 5 for each of one or more, but more typically a plurality of desired or possible {ΔT_min_i} sets.
Optimistic Solution
As described above, various embodiments of the present invention provide a bi-level programming solution philosophy which includes or otherwise provides for: a close to optimal CHP, distillation sequence, and HEN as well as the process conditions leading to it; a close to optimal CHP, distillation sequence and minimum number of heat exchangers units, as well as the process conditions leading to it; and a close to optimal CHP, distillation sequence and minimum heat exchangers units surface area, as well as the process conditions leading to it.
FIG. 29 illustrates a portion of a first exemplary HEN configuration illustrating a problem for solving for optimal process conditions and HEN surface area (A or SA) of each heat exchanger forming the HEN network. FIG. 30 illustrates the results of solving for optimal process conditions and corresponding HEN surface area configuration whereby the optimal process conditions solution was solved first and the at least semi-optimal utility system solution was solved second in view of the requirements of the optimal process configuration, i.e., the process units/system was the “leader” and the utility system was the “follower.” FIG. 31 illustrates the reverse situation whereby the optimal utility system solution was solved first and the process conditions solution was solved second, i.e., the utility system was that “leader” and the process units/system was the “follower.” In this configuration, the heating duty (Qh) of utility #3 (U3) provided with the other utilities, was given from the utility CHP solution to be 275 kW—a 25 kW reduction over that of the process condition leader configuration, but with a trade-off of a higher total heat exchanger surface area requirement.
FIG. 32 illustrates a portion of a second exemplary HEN configuration illustrating a problem that includes solving for optimal process conditions and HEN surface area of each heat exchanger forming the HEN network. FIG. 33 illustrates the results of solving for optimal process conditions and corresponding HEN surface area configuration whereby the optimal process conditions solution was solved first and the utility system solution was solved second in view of the requirements of the optimal process configuration, i.e., the process units/system was the “leader” and the utility system was the “follower.” FIG. 34 illustrates the reverse situation whereby the optimal utility system solution was solved first and the process conditions solution was solved second, i.e., the utility system was that “leader” and the process units/system was the “follower.” In this configuration, the heating duty (Qh) of utility #3 (U3) and cooling duty (Qc) were given from the utility CHP solution to be Qh=275 kW and Qc=75 kW—a 25 kW reduction for each of the two utilities over that of the process condition leader configuration along with associated reductions in utility capital costs, but with a trade-off of a higher total heat exchanger surface area requirement and corresponding associated capital costs.
FIG. 35 illustrates a portion of a third exemplary HEN configuration illustrating a problem for solving for optimal process conditions and HEN surface area of each heat exchanger forming the HEN network of FIG. 29, but with a higher ΔTmin−i. FIG. 36 illustrates the results of solving for optimal process conditions and corresponding HEN surface area configuration whereby the optimal process conditions solution was solved first and the utility system solution was solved second in view of the requirements of the optimal process configuration, i.e., the process units/system was the “leader” and the utility system was the “follower.” FIG. 37 illustrates the reverse situation whereby the optimal utility system solution was solved first and the process conditions solution was solved second, i.e., the utility system was that “leader” and the process units/system was the “follower.” In this configuration, there is no appreciable change in total utility requirements, but the resultant configuration has a larger total heat exchanger surface area requirement, and thus, larger capital costs which would need to be more than offset by the reduction in utility system capital costs over that which would be required if the surface area was not increased over that of the previous problem for utility system leader configuration to be the optimal configuration.
FIG. 38 illustrates a portion of a third exemplary HEN configuration illustrating a problem for solving for optimal process conditions and HEN surface area of each heat exchanger forming the HEN network of FIGS. 29 and 37, but with an even higher ΔTmin−i. FIG. 39 illustrates the results of solving for optimal process conditions and corresponding HEN surface area configuration whereby the optimal process conditions solution was solved first and the utility system solution was solved second in view of the requirements of the optimal process configuration, i.e., the process units/system was the “leader” and the utility system was the “follower.” FIG. 40 illustrates the reverse situation whereby the optimal utility system solution was solved first and the process conditions solution was solved second, i.e., the utility system was that “leader” and the process units/system was the “follower.” In this configuration, there is no appreciable change in total utility requirements, but the resultant configuration has a larger total heat exchanger surface area requirement, and thus, larger capital costs which would need to be more than offset by the reduction in utility system capital costs over that which would be required if the surface area was not increased over that of the previous problem for utility system leader configuration to be the optimal configuration.
Bi-level Programming Solution Philosophy
Not only do various embodiments of the present invention beneficially provide the above-described optimistic solution and the pessimistic solution (currently adapted, whereby the CHP decision-maker finds its optimum under process giving conditions from the process engineering department), such embodiments can provide a “tie-cooperation approach” which can provide: information to consider the least capital intensive distillation sequence and/or HEN versus best CHP utility cost during the decision-making of the CHP (e.g., through simple iterations between departments); a Pareto diagram for a partially decentralized company for company-Y best decision-making using {ΔTmin−i} as a marching variable; and a Pareto diagram for both, e.g., chemical company and third-party utility/CHP company having a limited budget/investment proposal, for mutual decision-making.
Computer Readable Medium Implementation
It is important to note that while embodiments of the present invention have been described in the context of a fully functional system/apparatus, those skilled in the art will appreciate that the mechanism of at least portions of the present invention and/or aspects thereof are capable of being distributed in the form of a non-transitory computer readable medium storing/containing or otherwise embodying instructions in a variety of forms for execution on one or more processors, or the like, and that embodiments of the present invention apply equally regardless of the particular type of media used to actually carry out the distribution. Non-transitory computer readable medium or media which is understood to mean includes all forms of computer readable storage media that do not fall under the category of being non-statutory subject matter, in general, or take the form of a propagating signal per se, in particular. Examples of the non-transitory computer readable media include but are not limited to: nonvolatile, hard-coded type media such as read only memories (ROMs), CD-ROMs, and DVD-ROMs, or erasable, electrically programmable read only memories (EEPROMs), recordable type media such as floppy disks, hard disk drives, CD-R/RWs, DVD-RAMs, DVD-R/RWs, DVD+R/RWs, HD-DVDs, memory sticks, mini disks, laser disks, Blu-ray disks, flash drives, and other newer types of memories, and in certain circumstances, transmission type media such as digital and analog communication links capable of storing/containing or otherwise embodying the instructions, to the exclusion of a propagating signal per se. For example, such media can store or otherwise contain both operating instructions and operations instructions related to the various program code 51 and the method steps, described above.
This non-provisional patent application claims priority to and the benefit of U.S. Provisional Patent Application No. 61/612,470, filed on Mar. 19, 2012, titled “System, Method, and Computer Program For Simultaneous Processing Utility Systems Synthesis in Partially and Fully Decentralized Environments,” and is related to U.S. Non-Provisional patent application titled “Methods for Simultaneous Process and Utility Systems Synthesis in Partially and Fully Decentralized Environments,” filed on Feb. 1, 2013, as U.S. application Ser. No. 13/757,467, each incorporated by reference in its entirety. See Appendix 1 for a list of additional related applications.
In the drawings and specification, there have been disclosed a typical preferred embodiment of the invention, and although specific terms are employed, the terms are used in a descriptive sense only and not for purposes of limitation. The invention has been described in considerable detail with specific reference to these illustrated embodiments. It will be apparent, however, that various modifications and changes can be made within the spirit and scope of the invention as described in the foregoing specification.
APPENDIX 1
Related Applications
The following applications and patents are each Incorporated by reference in their entirety:

    U.S. application Ser. No. 12/767,217, filed Apr. 26, 2010, now U.S. Pat. No. 8,032,262 titled System, Method, and Program Product for Synthesizing Non-Constrained and Constrained Heat Exchanger Networks, is a continuation-in-part of U.S. application Ser. No. 12/575,743 filed on Oct. 8, 2009, which claims priority to and the benefit of U.S. Provisional Application No. 61/256,754 filed on Oct. 30, 2009;
    U.S. patent application Ser. No. 12/575,743, filed Oct. 8, 2009, now U.S. Pat. No. 7,729,809, titled “System, Method, and Program Product for Targeting and Identification of Optimal Process Variables in Constrained Energy Recovery Systems”;
    U.S. patent application Ser. No. 11/768,084, filed Jun. 25, 2007, titled “System, Method, and Program Product for Targeting an Optimal Driving Force Distribution in Energy Recovery Systems”, now U.S. Pat. No. 7,698,022, which claims priority to and the benefit of U.S. Provisional Application No. 60/816,234, filed Jun. 23, 2006, titled “Method and Program Product for Targeting and Optimal Driving Force Distribution in Energy Recovery Systems,” each incorporated herein by reference in its entirety.
    U.S. patent application Ser. No. 12/715,255, filed Mar. 1, 2010, now U.S. Pat. No. 7,873,443 titled System, Method and Program Product for Targeting and Optimal Driving Force Distribution in Energy Recovery Systems, which claims priority to and the benefit of U.S. application Ser. No. 11/768,084 filed on Jun. 25, 2007;
    U.S. application Ser. No. 12/480,415 filed Jun. 8, 2009, titled Method and Software for Global Targeting of Process Utilities under Varying Conditions, which claims priority to and the benefit of U.S. Provisional Application 61/208,874 filed on Jun. 6, 2008, U.S. Provisional Application No. 61/171,030, filed on Apr. 20, 2009; and U.S. Provisional Application 61/172,698 filed on Apr. 24, 2009.
    U.S. application Ser. No. 12/898,461, filed Oct. 5, 2010, titled Systems, Program Product, and Methods for Synthesizing Heat Exchanger Networks that Exhibit Life-Cycle Switchability and Flexibility under all Possible Combinations of Process Variations; which claims priority to and the benefit of U.S. Provisional Application No. 61/356,900, filed on Jun. 21, 2010;
    U.S. application Ser. No. 12/767,275, filed Apr. 26, 2010, titled System, Method, and Program Product For Synthesizing Non-Thermodynamically Constrained Heat Exchanger Networks, which claims priority to and the benefit of U.S. Provisional Application No. 61/256,754, filed on Oct. 30, 2009;
    U.S. application Ser. No. 12/767,315, filed Apr. 26, 2010, tiled System, Method, and Program Product for Synthesizing Heat Exchanger Networks and Identifying Optimal Topography for Future Retrofit, which claims priority to and the benefit of U.S. Provisional Application No. 61/256,754:
    U.S. application Ser. No. 12/898,484 filed Oct. 5, 2010, titled Systems, Program Product, and Methods for Targeting Optimal Process Conditions that Render an Optimal Heat Exchanger Network Design Under Varying Conditions, which claims priority to and the benefit of U.S. Provisional Application No. 61/356,900, filed on Jun. 21, 2010;
    U.S. application Ser. No. 12/898,475, filed Oct. 5, 2010, titled Systems, Program Product, and Methods for Synthesizing Heat Exchanger Networks that Account for Future Higher Levels of Disturbances and Uncertainty, and Identifying Optimal Topology for Future Retrofit, which claims priority to and benefit of U.S. Provisional Patent Application 61/356,900 filed on Jun. 21, 2010.
    U.S. application Ser. No. 13/041,057, filed Mar. 4, 2011, titled “Systems and Program Product For Heat Exchanger Network Energy Efficiency Assessment and Lifetime Retrofit”, which is a continuation-in-part of and claims priority to and the benefit of U.S. application Ser. No. 12/767,315 filed on Apr. 26, 2010.
    U.S. application Ser. No. 13/041,087, filed Mar. 4, 2011, titled “Methods For Heat Exchanger Network Energy Efficiency Assessment and Lifetime Retrofit”, which is a continuation-in-part of and claims priority to and the benefit of U.S. application Ser. Nos. 12/767,315, 12/767,275, 12/767,217, filed on Apr. 26, 2010.


APPENDIX 2
Original Sets:


    HDR: Steam headers.
    BLR: Boilers.
    COGEN: Cogeneration Units.
    SOLAR: Solar thermal units
    STG: Steam turbine generator units.
    CSTG: Condensing steam turbine generator units.
    BOS: Break over stations.
    COND: Condensers
    PST: Process steam turbines.
    MSwitch: Motors switchable to process steam turbines
    PSGen: Process steam generators.
    PCDem: Process cooling demand
    FUEL: Fuel types

Derived Sets:

    BLRHDR[BLR×HDR]: Matching between each boiler and the header it is connected to.
    COGENHDR[COGEN×HDR]: Matching between each cogen and the header connected to.
    SOLARHDR[SOLAR×HDR]: Matching between each Solar thermal unit and the header connected to.
    STGInHDR[STG×HDR]: Input steam header to the STG.
    STGOutHDR[STG×HDR]: Output steam header from the STG.
    CSTGHDR[CSTG×HDR]: Input steam header to the CSTG.
    BOSIn[BOS×HDR]: Input steam header to the breakover station.
    BOSOut[BOS×HDR]: Output steam header of the breakover station.
    CONDHDR[COND×HDR]: Condenser Input steam header.
    DEAStmHDR[HDR]: Deaerator steam header.
    PSTInHDR[PST×HDR]: Inlet steam header to process steam turbine.
    PSTOutHDR[PST×HDR]: Outlet header from process steam turbine.
    PSGenHDR[PsGen×HDR]: Outlet header from process steam generator.

Parameters



Seam Headers:

    HDRPres[iεHDR]: Header Pressure.
    HDRTemp[iεHDR]: Header Temperature.
    HDRLoss[iεHDR]: % loss of steam from header.
    HDRRet[iεHDR]: Steam returning from process to header (flash recovery).
    HDRBfw[iεHDR]: BFW to a steam header
    HDRMaxDem[iεHDR]: Maximum process steam demand on a header
    HDRMinDem[iεHDR]: Minimum process steam demand on a header
    HDREnth[iεHDR]: Enthalpy of a steam header.
    HDRCost[iεHDR]: Cost of header


Boilers:

    BLRCap[iεBLR]: Steam generation capacity of boiler.
    BLRMin[iεBLR]: Minimum allowable steam generation rate for an operating boiler.
    BLRc0[iεBLR]: Boiler fuel consumption constant coefficient.
    BLRc1[iεBLR]: Boiler fuel consumption linear coefficient.
    BLRc2[iεBLR]: Boiler fuel consumption quadratic coefficient.
    BLRCC[iεBLR]: Boiler's cycle of concentration.
    BLRLoss[iεBLR]: Boiler's % loss of boiler feed water.
    BLRPrice[tεBLR]: Annualized capital cost of boiler of the given specs.


Cogeneration:

    COGENCap[iεCOGEN]: Power generation capacity of the cogen unit.
    COGENMin[iεCOGEN]: Minimum power generation of an operational cogen unit.
    COGENStmRatio[iεCOGEN]: Steam to power ratio of the cogen unit.
    COGENc0[iεCOGEN]: Fuel consumption constant coefficient.
    COGENc1[iεCOGEN]: Fuel consumption linear coefficient.
    COGENc2[iεCOGEN]: Fuel consumption quadratic coefficient.
    COGENCC[iεCOGEN]: Cogen's cycle of concentration parameter for the HRSG.
    COGENLoss[iεCOGEN]: Water loss from cogen unit as a % of BFW.
    COGENPrice[iεCOGEN]: Annualized capital cost of a cogen units of the given specs.


Steam Turbine Generator:

    STGEff[iεSTG]: Isentropic efficiency of the steam turbine generator.
    STGCap[iεSTG]: Steam capacity of STG.
    STGLoss[iεSTG]: % of steam loss from STG.
    STGEnthIn[iεSTG]: STG input steam enthalpy from input steam header.
    STGEnthOut[iεSTG]: STG output steam enthalpy.
    STGPrice[iεSTG]: Annualized capital cost of an STG unit of the given specs.


Condensing Steam Turbine Generator:

    CSTGEff[iεCSTG]: Isentropic efficiency of the condensing steam turbine generator.
    CSTG[iεCSTG]: Steam flow capacity through the CSTG.
    CSTGLoss[iεCSTG]: % of steam loss from CSTG.
    CSTGStmEnthIn[iεCSTG]: Enthalpy of input steam to the CSTG.
    CSTGStmEnthOut[iεCSTG]: Enthalpy specs of the steam out of the CSTG (out of the STG stage).
    CSTGWtrEnth[iεCSTG]: Enthalpy specs of condensate water of CSTG.
    CSTGDutyEnth[iεCSTG]: Enthalpy of duty of condensed water of CSTG.
    CSTGPrice[iεCSTG]: Annualized capital cost of CSTG.


Breakover Station:

    BOSCap[iεBOS]: Seam capacity of the breakover station.
    BOSLoss[iεBOS]: % of steam loss from the BOS.
    BOSStm[iεBOS]: Steam flow to BOS
    BOSEnth[iεBOS]: Enthalpy of BOS steam.
    BOSPrice[iεBOS]: Annualized capital cost of a BOSS of the given specs.


Condenser:

    CONDCap[iεCOND]: Capacity of steam flow through the condenser.
    CONDLoss[iεCOND]: % loss of steam flow through the condenser.
    CONDStmx[εCOND]: Steam quality at the condenser output.
    CONDEnth[iεCOND]: Enthalpy of condenser output.
    CONDEnthDuty[iεCOND]: Enthalpy of condenser duty.
    CONDPrice[iεCOND]: Annualized capital cost of a condenser unit with the given specs.


Process Steam Turbine:

    PSTLoss[iεPST]: % loss of steam flow through PST.
    PSTDem[iεPST]: Steam demand for the PST.
    PSTEnth[iεPST]: Enthalpy of the outlet steam from the PST.


Switchable Motors to Process Steam Turbines:

    MSwitchHp[iεMSwitch]: HP required by process driven by switchable motors
    MSwitchEff[iεMSwitch]: Efficiency of switchable motors
    MSwitchMWToPST[iεMSwitch]: Power demand in motors switchable to Process Steam
    Turbines
    MSwitchStmToPST[iεMSwitch]: Steam flow to process steam turbines switchable with motor


Process Steam Generator:

    PSGenEnth[iεPST]: Enthalpy of outlet steam from process steam generator.
    PSGenMaxProd[iεHDR]: Maximum Steam production rate of process steam generator.
    PSGenMinProd[iεHDR]: Minimum Steam production rate of process steam generator.


Solar Thermal Unit:

    SOLARBfwEnth[iεSOLAR]: BFW Enthalpy to solar thermal unit
    SOLAROutEnth[iεSOLAR]: Enthalpy of Hot water or steam generated from solar thermal unit
    SOLARStmx[iεSOLAR]: Quality of steam generated from solar thermal unit
    SOLAREff[iεSOLAR]: Efficiency of Solar thermal unit
    SOLARRad[iεSOLAR]: Avg. Daily heat radiation from sun (kwh)
    SOLARCost[iεSOLAR]: Annualized cost of Solar thermal unit


Condensate System:

    SkimTankLoss: % losses from skim tank.
    MWWtrEnth: Enthalpy of the makeup water.
    CondRetFlow: Flow rate of the return condensate.
    CondCoolerLoss: % loss of condensate from condensate coolers.
    CondRetTargEnth: The target enthalpy of the return condensate.


Deaerator:

    DEAstm: Steam flow to the deaerator.
    DEAStmEnth: Enthalpy of the steam used in the deaerator.
    DEAOutWtrEnth: Enthalpy of the water leaving the deaerator.
    DEALoss: % loss of water from deaerator.
    DEAVentFlow: Flow rate of steam vented from the deaerator.
    DEAVentEnth: Enthalpy of the steam vented from the deaerator.


Power:

    ProcPwrDem: Power demand from the process (fixed and variable loads).
    FixedPwrDm: Power demand for other than process (buildings, etc).


Fuel:

    FuelType[iεFuel]: Type of fuels
    FuelMaxAv[iεFuel]: Maximum quantity of a fuel type
    FuelHV[iεFuel]: Heating value of a fuel type
    FuelPrice[iεFuel]: Price of a fuel type


Economics:

    PriceImpPwr: Imported power price.
    PriceExpPwr: Exported power price.
    PriceMUWtr: Makeup water price.

Variables:



Steam Headers:

    HDRIn[iεHDR]: Steam flow into header.
    HDROut[iεHDR]: Steam flow out of header.
    HDRBFW[iεHDR]: Boiler feed water injected into header to maintain header enthalpy specs.
    HDRDem[iεHDR]: Process steam demand on a header; Determine optimal values of Steam for heating demand at different headers within a given intervals


Boilers:

    BLROn[iεBLR]: Binary variable to determine whether the boiler is on or off (Binary).
    BLRBFW[iεBLR]: Boiler feed water into boiler.
    BLRSTM[iεBLR]: Boiler steam generation rate.
    BLRBD[iεBLR]: Blowdown from boiler.
    BLRFuel[iεBLR]: Boiler's fuel consumption rate.
    BLRNo[iεBLR]: Number of installed units of the boiler (Nonnegative integer).


Cogeneration:

    COGENOn[iεCOGEN]: Binary variable to determine whether a cogen unit is on or off (Binary).
    COGENPwr[iεCOGEN]: Power generation rate from cogen unit.
    COGENBFW[iεCOGEN]: Boiler feed water to cogen unit.
    COGENStm[iεCOGEN]: Steam generation rate from cogen unit.
    COGENBD[iεCOGEN]: Water blowdown rate from cogen unit.
    COGENFuel[iεCOGEN]: Fuel consumption rate of cogen unit.
    COGENNo[iεCOGEN]: Number of installed cogen units (Nonnegative integer)


Solar Steam Generator:

    SOLARStm[so]: Steam flow from solar thermal unit
    SOLARNo[so]: Binary variable of solar thermal unit


Steam Turbine Generator:

    STGStm[iεSTG]: Steam flow through the STG.
    STGPwr[iεSTG]: Power generation from the STG.
    STGNo[iεSTG]: Number of installations of the STG.


Condensing Steam Turbine Generator:

    CSTGStm[iεCSTG]: Steam flow rate through CSTG.
    CSTGPwr[iεCSTG]: Power generated from CSTG.
    CSTGDuty[iεCSTG]: Duty of CSTG condenser.
    CSTGNo[iεCSTG]: Number of installed units of CSTG of the given specs.


Switchable Motors to Process Steam Turbines:

    MSwitchOnOff[iεMSwitch]: binary variable (0/1) of a switchable motor
    PSTSwitchStm[iεMSwitch]: steam flow to process steam turbines switchable with motors


Break over Station:

    BOSStm[iεBOS]: Steam flow through the BOS.
    BOSNo[iεBOS]: Number of installed BOS units.


Condenser:

    CONDStm[iεCOND]: Steam flowrate through the condenser.
    CONDDuty[iεCOND]: Duty of the condenser.
    CONDNo[iεCOND]: Number of installed condenser units.


Condensate System:

    SkimTankInFlow: Input flowrate to skim tank.
    SkimTankOutFlow: Output flowrate from skim tank.
    SkimTankOutEnt: Enthalpy of skim tank outlet.
    MUWtrFlow: Flow rate of makeup water to the system.


Deaerator.

    DEAStmFlow: Deaerator steam flow rate.
    DEAOutWtrFlow: Deaerator output water flow rate.


Process Steam Generator:

    PSGenProd[iεPST]: Steam production rate of process steam generator.


Boiler Feed water System:

    BFWFlow: Boiler feedwater flow rate.


Process Cooling Demand:

    PCoolingDem[iεPCDem]: Determine optimal values of: steam for cooling demand, water cooling, air-cooling and refrigeration cooling demands


APPENDIX 3
Model:


    H1_in_temp=510.0;
    H1_in_fcb=1.0;
    H2_in_temp>=610.0; H2_in_temp>=630.0;
    H2_in_fcb=2.0;
    H1_out_temp=210.0;
    H1_out_fcb=1.0;
    H2_out_temp=230.0;
    C1_in_temp=400.0;
    C1_in_fcb=4.0;
    C2_out_temp=460.0;
    C2_out_fcb=2.0;
    CU1_temp=220.0;
    CU1_fcb>=60.0; CU1_fcb<=100.0;
    HU2_temp=480.0;
    HU2_fcb=120.0;
    HU3_temp=630.0;
    HU3_fcb>=260.0; HU3_fcb<=300.0;
    CU1_exiting_temp=221.0;
    CU1_exiting_fcb>=60.0; CU1_exiting_fcb<=100.0;
    HU2_exiting_tem=479.0;
    HU2_exiting_fc=120.0;
    HU3_exiting_tem=629.0;
    HU3_exiting_fcb>=260.0; HU3_exiting_fcb<=300.0;
    E1_q=100.0;
    E1_theta_1>=10.0;
    E1_theta_2>=10.0;
    E1_u=1.2;
    E2_q=200.0;
    E2_theta_1>=10.0;
    E2_theta_2>=10.0;
    E2_u=1.5;
    E3_q>=400.0; E3_q<=440.0;
    E3_theta_1>=10.0;
    E3_theta_2>=10.0;
    E3_u=1.0;
    E4_q>=260.0; E4_q<=300.0;
    E4_theta_1>=10.0;
    E4_theta_2>=10.0;
    E4_u=2.0;
    U1_q>=60.0; U1_q<=100.0;
    U1_theta_1>=10.0;
    U1_theta_2>=10.0;
    U1_u=1.5;
    U2_q=120.0;
    U2_theta_1>=10.0;
    U2_theta_2>=10.0;
    U2_u=2.0;
    U3_q>=260.0; U3_q<=300.0;
    U3_theta_1>=10.0;
    U3_theta_2>=10.0;
    U3_u=2.0;
    !J1;
    C1_in_fcb=J1_to_E3_c_in_fcb+J1_to_J2_fcb+J1_to_E1_c_in_fcb;
    J1_to_E3_c_in_temp=C1_in_temp;J1_to_J2_temp=C1_in_temp;J1_to_E1_c_in_temp=C1_in_temp;
    !J2;
    E3_c_out_to_J2_fcb+J1_to_J2_fcb+E1_c_out_to_J2_fcb=J2_to_U3_c_in_fcb;
    E3_c_out_to_J2_fcb*E3_c_out_to_J2_temp+J1_to_J2_fcb*J1_to_J2_temp+E1_c_out_to_J2_fcb*E1_c_out_to_J2_temp=J2_to_U3_c_in_fcb*J2_to_U3_c_in_temp;
    !J3;
    C2_in_fcb=J3_to_E2_c_in_fcb+J3_to_E4_c_in_fcb;
    J3_to_E2_c_in_temp=C2_in_temp;J3_to_E4_c_in_temp=C2_in_temp;
    !J4;
    E2_c_out_to_J4_fcb+E4_c_out_to_J4_fcb=J4_to_U2_c_in_fcb;
    E2_c_out_to_J4_fcb*E2_c_out_to_J4_temp+E4_c_out_to_J4_fcb*E4_c_out_to_J4_temp=J4_to_U2_c_in_fcb*J4_to_U2_c_in_temp;
    !E1_h_in;
    H1_in_fcb=E1_h_in_to_h_fcb;
    E1_h_in_to_h_temp=H1_in_temp;
    !E1_h_out;
    E1_h_to_h_out_fcb=E1_h_out_to_E2_h_in_fcb;
    E1_h_out_to_E2_h_in_temp=E1_h_to_h_out_temp;
    !E1_c_in;
    J1_to_E1_c_in_fcb=E1_c_in_to_c_fcb;
    E1_c_in_to_c_temp=J1_to_E1_c_in_temp;
    !E1_c_out;
    E1_c_to_c_out_fcb=E1_c_out_to_J2_fcb;
    E1_c_out_to_J2_temp=E1_c_to_c_out_temp;
    !E2_h_in;
    E1_h_out_to_E2_h_in_fcb=E2_h_in_to_h_fcb;
    E2_h_in_to_h_temp=E1_h_out_to_E2_h_in_temp;
    !E2_h_out;
    E2_h_to_h_out_fcb=H1_out_fcb;
    H1_out_temp=E2_h_to_h_out_temp;
    !E2_c_in;
    J3_to_E2_c_in_fcb=E2_c_in_to_c_fcb;
    E2_c_in_to_c_temp=J3_to_E2_c_in_temp;
    !E2_c_out;
    E2_c_to_c_out_fcb=E2_c_out_to_J4_fcb;
    E2_c_out_to_J4_temp=E2_c_to_c_out_temp;
    !E3_h_in;
    H2_in_fcb=E3_h_in_to_h_fcb;
    E3_h_in_to_h_temp=H2_in_temp;
    !E3_h_out;
    E3_h_to_h_out_fcb=E3_h_out_to_E4_h_in_fcb;
    E3_h_out_to_E4_h_in_temp=E3_h_to_h_out_temp;
    !E3_c_in;
    J1_to_E3_c_in_fcb=E3_c_in_to_c_fcb;
    E3_c_in_to_c_temp=J1_to_E3_c_in_temp;
    !E3_c_out;
    E3_c_to_c_out_fcb=E3_c_out_to_J2_fcb;
    E3_c_out_to_J2_temp=E3_c_to_c_out_temp;
    !E4_h_in;
    E3_h_out_to_E4_h_in_fcb=E4_h_in_to_h_fcb;
    E4_h_in_to_h_temp=E3_h_out_to_E4_h_in_temp;
    !E4_h_out;
    E4_h_to_h_out_fcb=E4_h_out_to_U1_h_in_fcb;
    E4_h_out_to_U1_h_in_temp=E4_h_to_h_out_temp;
    !E4_c_in;
    J3_to_E4_c_in_fcb=E4_c_in_to_c_fcb;
    E4_c_in_to_c_temp=J3_to_E4_c_in_temp;
    !E4_c_out;
    E4_c_to_c_out_fcb=E4_c_out_to_J4_fcb;
    E4_c_out_to_J4_temp=E4_c_to_c_out_temp;
    !U1_h_in;
    E4_h_out_to_U1_h_in_fcb=U1_h_in_to_h_fcb;
    U1_h_in_to_h_temp=E4_h_out_to_U1_h_in_temp;
    !U1_h_out;
    U1_h_to_h_out_fcb=H2_out_fcb;
    H2_out_temp=U1_h_to_h_out_temp;
    !U1_c_in;
    CU1_fcb=U1_c_in_to_c_fcb;
    U1_c_in_to_c_temp=CU1_temp;
    !U1_c_out;
    U1_c_to_c_out_fcb=CU1_exiting_fcb;
    CU1_exiting_temp=U1_c_to_c_out_temp;
    !U2_h_in;
    HU2_fcb=U2_h_in_to_h_fcb;
    U2_h_in_to_h_temp=HU2_temp;
    !U2_h_out;
    U2_h_to_h_out_fcb=HU2_exiting_fcb;
    HU2_exiting_temp=U2_h_to_h_out_temp;
    !U2_c_in;
    J4_to_U2_c_in_fcb=U2_c_in_to_c_fcb;
    U2_c_in_to_c_temp=J4_to_U2_c_in_temp;
    !U2_c_out;
    U2_c_to_c_out_fcb=C2_out_fcb;
    C2_out_temp=U2_c_to_c_out_temp;
    !U3_h_in;
    HU3_fcb=U3_h_in_to_h_fcb;
    U3_h_in_to_h_temp=HU3_temp;
    !U3_h_out;
    U3_h_to_h_out_fcb=HU3_exiting_fcb;
    HU3_exiting_temp=U3_h_to_h_out_temp;
    !U3_c_in;
    J2_to_U3_c_in_fcb=U3_c_in_to_c_fcb;
    U3_c_in_to_c_temp=J2_to_U3_c_in_temp;
    !U3_c_out;
    U3_c_to_c_out_fcb=C1_out_fcb;
    C1_out_temp=U3_c_to_c_out_temp;
    !E1 inner Junctions;
    E1_h_in_to_h_fcb=E1_h_to_h_out_fcb;
    E1_c_in_to_c_fcb=E1_c_to_c_out_fcb;
    !E2 inner Junctions;
    E2_h_in_to_h_fcb=E2_h_to_h_out_fcb;
    E2_c_in_to_c_fcb=E2_c_to_c_out_fcb;
    !E3 inner Junctions;
    E3_h_in_to_h_fcb=E3_h_to_h_out_fcb;
    E3_c_in_to_c_fcb=E3_c_to_c_out_fcb;
    !E4 inner Junctions;
    E4_h_in_to_h_fcb=E4_h_to_h_out_fcb;
    E4_c_in_to_c_fcb=E4_c_to_c_out_fcb;
    !U1 inner Junctions;
    U1_h_in_to_h_fcb=U1_h_to_h_out_fcb;
    U1_c_in_to_c_fcb=U1_c_to_c_out_fcb;
    !U2 inner Junctions;
    U2_h_in_to_h_fcb=U2_h_to_h_out_fcb;
    U2_c_in_to_c_fcb=U2_c_to_c_out_fcb;
    !U3 inner Junctions;
    U3_h_in_to_h_fcb=U3_h_to_h_out_fcb;
    U3_c_in_to_c_fcb=U3_c_to_c_out_fcb;
    !Heat Exchangers;
    !E1;
    E1_q=(E1_h_in_to_h_temp−E1_h_to_h_out_temp)*E1_h_in_to_h_fcb;
    E1_q=(E1_c_to_c_out_temp−E1_c_in_to_c_temp)*E113 c_in_to_c_fcb;
    E1_theta_1=E1_h_to_h_out_temp−E1_c_in_to_c_temp;
    E1_theta_2=E1_h_in_to_h_temp−E1_c_to_c_out_temp;
    E1_lmtd=((E1_theta_1*E1_theta_2)*(E1_theta_1+E1_theta_2)/2)^0.333;
    E1_area=E1_q/(E1_u*E1_lmtd);
    !E2;
    E2_q=(E2_h_in_to_h_temp−E2_h_to_h_out_temp)*E2_h_in_to_h_fcb;
    E2_q=(E2_c_to_c_out_temp−E2_c_in_to_c_temp)*E2_c_in_to_c_fcb;
    E2_theta_1=E2_h_to_h_out_temp−E2_c_in_to_c_temp;
    E2_theta_2=E2_h_in_to_h_temp−E2_c_to_c_out_temp;
    E2_lmtd=((E2_theta_1*E2_theta_2)*(E2_theta_1+E2_theta_2)/2)^0.333;
    E2_area=E2_q/(E2_u*E2_lmtd);
    !E3;
    E3_q=(E3_h_in_to_h_temp−E3_h_to_h_out_temp)*E3_h_in_to_h_fcb;
    E3_q=(E3_c_to_c_out_temp−E3_c_in_to_c_temp)*E3_c_in_to_c_fcb;
    E3_theta_1=E3_h_to_h_out_temp−E3_c_in_to_c_temp;
    E3_theta_2=E3_h_in_to_h_temp−E3_c_to_c_out_temp;
    E3_lmtd=((E3_theta_1*E3_theta_2)*(E3_theta_1+E3_theta_2)/2)^0.333;
    E3_area=E3_q/(E3_u*E3_lmtd);
    !E4;
    E4_q=(E4_h_in_to_h_temp−E4_h_to_h_out_temp)*E4_h_in_to_h_fcb;
    E4_q=(E4_c_to_c_out_temp−E4_c_in_to_c_temp)*E4_c_in_to_c_fcb;
    E4_theta_1=E4_h_to_h_out_temp−E4_c_in_to_c_temp;
    E4_theta_2=E4_h_in_to_h_temp−E4_c_to_c_out_temp;
    E4_lmtd=((E4_theta_1*E4_theta_2)*(E4_theta_1+E4_theta_2)/2)^0.333;
    E4_area=E4_q/(E4_u*E4_lmtd);
    !U1;
    U1_q=(U1_h_in_to_h_temp−U1_h_to_h_out_temp)*U1_h_in_to_h_fcb;
    U1_q=(U1_c_to_c_out_temp−U1_c_in_to_c_temp)*U1_c_in_to_c_fcb;
    U1_theta_1=U1_h_to_h_out_temp−U1_c_in_to_c_temp;
    U1_theta_2=U1_h_in_to_h_temp−U1_c_to_c_out_temp;
    U1_lmtd=((U1_theta_1*U1_theta_2)*(U1_theta_1+U1_theta_2)/2)^0.333;
    U1_area=U1_q/(U1_u*U1_lmtd);
    !U2;
    U2_q=(U2_h_in_to_h_temp−U2_h_to_h_out_temp)*U2_h_in_to_h_fcb;
    U2_q=(U2_c_to_c_out_temp−U2_c_in_to_c_temp)*U2_c_in_to_c_fcb;
    U2_theta_1=U2_h_to_h_out_temp−U2_c_in_to_c_temp;
    U2_theta_2=U2_h_in_to_h_temp−U2_c_to_c_out_temp;
    U2_lmtd=((U2_theta_1*U2_theta_2)*(U2_theta_1+U2_theta_2)/2)^0.333;
    U2_area=U2_q/(U2_u*U2_lmtd);
    !U3;
    U3_q=(U3_h_in_to_h_temp−U3_h_to_h_out_temp)*U3_h_in_to_h_fcb;
    U3_q=(U3_c_to_c_out_temp−U3_c_in_to_c_temp)*U3_c_in_to_c_fcb;
    U3_theta_1=U3_h_to_h_out_temp−U3_c_in_to_c_temp;
    U3_theta_2=U3_h_in_to_h_temp−U3_c_to_c_out_temp;
    U3_lmtd=((U3_theta_1*U3_theta_2)*(U3_theta_1+U3_theta_2)/2)^0.333;
    U3_area=U3_q/(U3_u*U3_lmtd);
    !Objective;
    MIN=E1_area+E2_area+E3_area+E4_area+U1_area+U2_area +U3_area;
    calc:
    @SET('GLOBAL',1);
    endcalc
    end
    Go
    DIVERT solution.txt
    SOLUTION
    RVRT
    QUIT;

## Claims

1. A computer program to perform simultaneous synthesis of combined heat and power utility and process subsystems for a process facility, the computer program comprising processor readable code stored on a non-transitory computer readable medium, the processor readable code comprising a set of instructions, that when executed by one or more processors, cause the one or more processors to perform operations comprising:
defining model input data for a combined heat and power utility and process subsystem mathematical program model, the model comprising a cost function configured to minimize a total cost of combined heat and power utility and process subsystems of a process facility, the combined heat and power utility and process subsystems comprising a utility subsystem configured to generate steam and a process subsystem configured to consume steam generated by the utility subsystem, and the model input data comprising a plurality of rigorously bounded decision variables for one or more key process of the process facility, the plurality of rigorously bounded process variables comprising, for each of the one or more key processes:
process steam demand intervals for the key process, the process steam demand intervals comprising a process steam demand interval for each of a plurality of different steam pressure levels, wherein each process steam demand interval for a steam pressure level comprises an upper bound and a lower bound of process steam demand for the key process at the steam pressure level; and
a process power generation interval for the key process, the process power generation interval comprising an upper bound and a lower bound of power generation for the key process; and

processing the model input data through application of the combined heat and power utility and process subsystem mathematical program model to determine model output data comprising process subsystem conditions and utility subsystem conditions, the process subsystem conditions comprising, for each of the one or more key processes:
process steam demand load values corresponding to the process steam demand intervals for the key process, the process steam demand load values comprising a process steam demand load value for each of the plurality of different steam pressure levels; and
a process power generation value corresponding to the process power generation interval.

2. A computer program as defined in claim 1,
wherein the process subsystem conditions and the utility subsystem conditions are determined simultaneously through application of the combined heat and power utility and process subsystem mathematical program model; and
wherein the utility subsystem conditions include combined heat and power utility subsystem steam and power generation and allocation values.

3. A computer program as defined in claim 1,
wherein the utility subsystem conditions comprise a number of steam headers for the utility subsystem and operating conditions of each steam header.

4. A computer program as defined in claim 1,
wherein the utility subsystem conditions comprise a number of boilers for the utility subsystem and designated capacities of each boiler.

5. A computer program as defined in claim 1,
wherein the utility subsystem conditions comprise a number of steam turbine generators for the utility subsystem and capacity of each steam turbine generator.

6. A computer program as defined in claim 1,
wherein the utility subsystem conditions comprise a number of cogeneration units for the utility subsystem and designated capacities of each cogeneration unit.

7. A computer program as defined in claim 1,
wherein the process subsystem conditions comprise a number of electrical motors and steam turbines for process equipment driving and size of each.

8. A computer program as defined in claim 1,
wherein the utility subsystem conditions comprise a size and location of solar power system for the utility subsystem.

9. A computer program as defined in claim 1, the operations further comprising:
determining a first cost associated with the utility subsystem;
determining a second cost associated with the process subsystem;
determining a greater of the first cost associated with the utility subsystem and the second cost second cost associated with the process subsystem;
in response to determining that the first cost associated with the utility subsystem is greater than the second cost associated with the process subsystem:
determining the utility subsystem conditions; and
determining the process subsystem conditions based on the utility subsystem conditions determined; and

in response to determining that the first cost associated with the utility subsystem is less than the second cost associated with the process subsystem:
determining the process subsystem conditions; and
determining the utility subsystem conditions based on the process subsystem conditions determined.

10. A computer program as defined in claim 1, wherein the operations further comprise one or more of the following:
defining a number of steam headers and conditions thereof of the utility subsystem; and
determining a shaft work network configuration of the utility subsystem, wherein the shaft work network configuration is indicative of a configuration of physical equipment.

11. A computer program as defined in claim 1,
wherein determining the process subsystem conditions and the utility subsystem conditions comprises determining a main design and operating conditions of a heat exchanger network (HEN) design and distillation network sequence for the process facility for possible combinations of process operation parametric changes for fully decentralized, partially decentralized, or both fully and partially decentralized process and utility systems design environments.

12. A computer program as defined in claim 11,
wherein the main design and operating conditions of a HEN network design and distillation network sequence is determined using sources of fuel or energy available to the utility subsystem; and
wherein the main design and operating conditions of a HEN network design and distillation network sequence is determined responsive to receiving, determining, or receiving and determining a plurality of sets of values each defining a potential range of attribute values of a separate one of a corresponding plurality of process operation parametric changes.

13. A computer program as defined in claim 1,
wherein the process subsystem conditions render a minimum number of heat exchanger network (HEN) units for the process facility for possible distillation columns sequences.

14. A computer program as defined in claim 1,
wherein the process subsystem conditions render a minimum total heat exchanger surface area of a heat exchanger network (HEN) design for the process facility and define the surface area of each heat exchanger unit in the HEN.

15. A computer program as defined in claim 14, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for:
one or more possible combinations of steam generation and demand loads' levels; and
one or more combinations of steam generation and demand loads' levels.

16. A computer program as defined in claim 14, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for:
one or more possible distillation columns sequences;
one or more possible distillation columns sequences in conjunction with one or more possible combinations of steam generation and demand loads' levels; and
one or more possible distillation columns sequences and one or more combinations of steam generation and demand loads' levels.

17. A computer program as defined in claim 14, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for:
one or more possible HEN configurations and one or more possible combinations of steam generation and demand loads' levels; and
one or more possible HEN configurations and one or more combinations of steam generation and demand loads' levels.

18. A computer program as defined in claim 1, the operations further comprising designing the process facility based on the process subsystem conditions and utility subsystem conditions.

19. A computer program to perform simultaneous synthesis of combined heat and power utility and process subsystems for chemical plant, the computer program comprising processor readable code stored on a non-transitory computer readable medium, the processor readable code comprising a set of instructions, that when executed by one or more processors, cause the one or more processors to perform operations comprising:
defining model input data for a combined heat and power utility and process subsystem mathematical program model, the model comprising a cost function configured to minimize a total cost of combined heat and power utility and process subsystems of a chemical plant, the combined heat and power utility and process subsystems comprising a utility subsystem configured to generate steam and a process subsystem configured to consume steam generated by the utility subsystem, and the model input data comprising a plurality of rigorously bounded process variables for one or more key process of the chemical plant, the plurality of rigorously bounded process variables comprising, for each of the one or more key processes:
process steam demand intervals for the key process, the process steam demand intervals comprising a process steam demand interval for each of a plurality of different steam pressure levels, wherein each process steam demand interval for a steam pressure level comprises an upper bound and a lower bound of process steam demand for the key process at the steam pressure level; and
a process power generation interval for the key process, the process power generation interval comprising an upper bound and a lower bound of power generation for the key process; and

processing, by a computer, the model input data through application of the combined heat and power utility and process subsystem mathematical program model to determine model output data comprising process subsystem conditions and utility subsystem conditions, the process subsystem conditions comprising, for each of the one or more key processes:
process steam demand load values corresponding to the process steam demand intervals for the key process, the process steam demand load values comprising a process steam demand load value for each of the plurality of different steam pressure levels; and
a process power generation value corresponding to the process power generation interval.

20. A computer program as defined in claim 19, wherein the process subsystem conditions and the utility subsystem conditions are determined simultaneously through application of the combined heat and power utility and process subsystem mathematical program model.

21. A computer program as defined in claim 19, the operations further comprising:
determining a first cost associated with the utility subsystem;
determining a second cost associated with the process subsystem;
determining a greater of the first cost associated with the utility subsystem and the second cost second cost associated with the process subsystem;
in response to determining that the first cost associated with the utility subsystem is greater than the second cost associated with the process subsystem;
determining the utility subsystem conditions; and
determining the process subsystem conditions based on the utility subsystem conditions determined; and

in response to determining that the first cost associated with the utility subsystem is less than the second cost associated with the process subsystem:
determining the process subsystem conditions; and
determining the utility subsystem conditions based on the process subsystem conditions determined.

22. A computer program as defined in claim 19, wherein the operations further comprise:
defining a number of steam headers and conditions thereof of the utility system; and
determining a shaft work network configuration of the utility system, wherein the shaft work network configuration is indicative of a configuration of physical equipment.

23. A computer program as defined in claim 19,
wherein determining the process subsystem conditions and the utility subsystem conditions comprises determining a main design and operating conditions of a heat exchanger network (HEN) design and distillation network sequence for the chemical plant under possible combinations of process operation parametric changes for fully decentralized, partially decentralized, or both fully and partially decentralized process and utility systems design environments.

24. A computer program as defined in claim 23, wherein the main design and operating conditions of a HEN network design and distillation network sequence is determined using sources of fuel or energy available to the chemical plant.

25. A computer program as defined in claim 23, wherein the main design and operating conditions of a HEN network design and distillation network sequence is determined responsive to receiving, determining, or receiving and determining a plurality of sets of values each defining a potential range of attribute values of a separate one of a corresponding plurality of process operation parametric changes.

26. A computer program as defined in claim 19, wherein the processing comprises determining the process subsystem conditions and the utility subsystem conditions based on possible combinations of process design parametric changes using sources of fuel or energy available to the utility system for both fully and partially decentralized process and utility systems design environments.

27. A computer program as defined in claim 19, wherein the processing comprises determining the process subsystem conditions and the utility subsystem conditions based on possible combinations of process design and operation parametric changes using sources of fuel or energy available to the utility system for both fully and partially decentralized process and utility systems design environments.

28. A computer program as defined in claim 19, wherein the processing comprises determining the process subsystem conditions and the utility subsystem conditions based on possible distillation sequences using sources of fuel or energy available to the utility system for both fully and partially decentralized process and utility systems design environments.

29. A computer program as defined in claim 19, wherein the processing comprises determining the process subsystem conditions and the utility subsystem conditions based on possible distillation sequences and under all reasonably possible combinations of process design and operation parametric changes using sources of fuel or energy available to the utility system for both fully and partially decentralized process and utility systems design environments.

30. A computer program as defined in claim 19, wherein the processing comprises determining the process subsystem conditions and the utility subsystem conditions based on possible heat exchanger network (HEN) designs using sources of fuel or energy available to the utility system for both fully and partially decentralized process and utility systems design environments.

31. A computer program as defined in claim 19, wherein the processing comprises determining the process subsystem conditions and the utility subsystem conditions based on possible heat exchanger network (HEN) designs and under possible combinations of process design and operation parametric changes using sources of fuel or energy available to the utility system for both fully and partially decentralized process and utility systems design environments.

32. A computer program as defined in claim 19, wherein the processing comprises determining the process subsystem conditions and the utility subsystem conditions based on possible distillation columns sequences and heat exchanger network (HEN) designs using sources of fuel or energy available to the utility system for both fully and partially decentralized process and utility systems design environments.

33. A computer program as defined in claim 19, wherein the processing comprises determining the process subsystem conditions and the utility subsystem conditions based on possible distillation columns sequences and heat exchanger network (HEN) designs under possible combinations of process design and operation parametric changes using sources of fuel or energy available to the utility system for both fully and partially decentralized process and utility systems design environments.

34. A computer program as defined in claim 19,
wherein the process subsystem conditions render a minimum number of heat exchanger network (HEN) units for the chemical plant for possible distillation columns sequences.

35. A computer program as defined in claim 19,
wherein the process subsystem conditions render a minimum total heat exchanger surface area of a heat exchanger network (HEN) design for the chemical plant and define the surface area of each heat exchanger unit in the HEN.

36. A computer program as defined in claim 35, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for one or more possible combinations of steam generation and demand loads' levels.

37. A computer program as defined in claim 35, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for one or more combinations of steam generation and demand loads' levels.

38. A computer program as defined in claim 35, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for one or more possible distillation columns sequences.

39. A computer program as defined in claim 35, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for one or more possible distillation columns sequences in conjunction with one or more possible combinations of steam generation and demand loads' levels.

40. A computer program as defined in claim 35, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for one or more possible distillation columns sequences and one or more combinations of steam generation and demand loads' levels.

41. A computer program as defined in claim 35, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for one or more possible HEN configurations and one or more possible combinations of steam generation and demand loads' levels.

42. A computer program as defined in claim 35, wherein determining the process subsystem conditions which render the minimum total heat exchanger surface area comprises determining process subsystem conditions for one or more possible HEN configurations and one or more combinations of steam generation and demand loads' levels.

43. A computer program as defined in claim 19, the operations further comprising further comprising designing the chemical plant based on the process subsystem conditions and utility subsystem conditions.

44. A system to perform simultaneous synthesis of combined heat and power utility and process subsystems for a process facility, the system comprising:
a combined heat and power utility and process subsystems analysis and design program computer having a processor, and memory coupled to the processor; and
combined heat and power utility and process subsystems analysis and design program code defining a computer program stored in the memory of the combined heat and power utility and process subsystems analysis and design program computer to perform simultaneous synthesis of combined heat and power utility and process subsystems for a process facility, the computer program including instructions that when executed by the combined heat and power utility and process subsystems analysis and design program computer, cause the computer to perform the operations of:
defining model input data for a combined heat and power utility and process subsystem mathematical program model, the model comprising a cost function configured to minimize a total cost of combined heat and power utility and process subsystems of a process facility, the combined heat and power utility and process subsystems comprising a utility subsystem configured to generate steam and a process subsystem configured to consume steam generated by the utility subsystem, and the model input data comprising a plurality of rigorously bounded decision variables for one or more key process of the process facility, the plurality of rigorously bounded process variables comprising, for each of the one or more key processes:
process steam demand intervals for the key process, the process steam demand intervals comprising a process steam demand interval for each of a plurality of different steam pressure levels, wherein each process steam demand interval for a steam pressure level comprises an upper bound and a lower bound of process steam demand for the key process at the steam pressure level; and
a process power generation interval for the key process, the process power generation interval comprising an upper bound and a lower bound of power generation for the key process; and

processing the model input data through application of the combined heat and power utility and process subsystem mathematical program model to determine model output data comprising process subsystem conditions and utility subsystem conditions, the process conditions comprising, for each of the one or more key processes:
process steam demand load values corresponding to the process steam demand intervals for the key process, the process steam demand load values comprising a process steam demand load value for each of the plurality of different steam pressure levels; and
a process power generation value corresponding to the process power generation interval.

45. A system as defined in claim 44,
wherein the process subsystem conditions and the utility subsystem conditions are determined simultaneously through application of the combined heat and power utility and process subsystem mathematical program model.

46. A system as defined in claim 44, wherein the process facility comprises a chemical plant.