WO2007149582A2

WO2007149582A2 - System, method, and program product for optimizing heat transfer in energy recovery systems

Info

Publication number: WO2007149582A2
Application number: PCT/US2007/014687
Authority: WO
Inventors: Mahmoud Bahy Noureldin; Ahmed Saleh Aseeri
Original assignee: Saudi Arabian Oil Company; Aramco Services Company
Priority date: 2006-06-23
Filing date: 2007-06-25
Publication date: 2007-12-27
Also published as: EP2400359A1; CN101512448B; CA2655973C; WO2007149582A3; US7698022B2; CN101512448A; US7873443B2; CA2655973A1; US20100223198A1; EP2400359B1; NO20085348L; EP2044492B1; NO333862B1; EP2044492A2; US20080015839A1

Abstract

A system (30), methods, and user-friendly program product (51) to calculate global energy utility targets and define optimal driving force distribution for a process or cluster of processes under all possible process changes and streams specific minimum temperature approach values (ΔTmini), simultaneously, and without enumeration, are provided. The program product (51) can utilize stream-specific minimum temperature approach values (ΔTmini), where the superscript i represents the specific hot stream, as the optimization parameters instead of the single global ΔTmin currently used, in addition to identifying the optimal operating conditions. The program product (51) can define optimal process conditions and an optimal driving force distribution in heat recovery systems, and can produce an optimal Pareto-curve that shows the rigorous trade-off between energy cost and capital cost for any energy recovery system.

Description

SYSTEM, METHOD, AND PROGRAM PRODUCT FOR TARGETING AND OPTIMAL DRIVING FORCE DISTRIBUTION IN ENERGY

RECOVERY SYSTEMS

BACKGROUND OF THE INVENTION

1. Related Applications

This non-provisional application claims priority to and the benefit of US Patent Application No. 60/816,234 incorporated herein by reference in its entirety.

2. Field of the Invention

[0001] This present invention relates to the field of energy recovery systems and related methods.

3. Background Of The Invention

[0002] Many different types of processes consume energy to obtain an output result, or to produce a required product or compound. For example, chemical processes consume energy to provide a desired result. For large scale processes which consume significant amounts of energy, it is preferable to minimize the energy consumed where possible. In the case of, for example, electrical energy generation systems or relatively large manufacturing plant or factories, it is preferable to optimize and potentially minimize the consumption of energy through careful operation, design or reconfiguration of the plant and equipment used.

[0003] For example, 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 temperature to a target temperature. This in turn will require the consumption of energy to cool specific streams, and also to heat other specific streams. [0004] The total energy employed or consumed by the processes can be optimized to a global minimal level through, for example, careful placement and configuration of specific material streams with respect to one another. For example, there may be the potential for hot streams that require cooling to be placed in proximity with cold streams that require heating. Thermal energy already present in streams that needs to be removed, or streams that need to have heat added, can therefore be associated with one another to optimize the energy consumption of the process. In addition, the minimum temperature differences between hot streams and cold streams upon optimization can also result in huge savings in energy consumption.

[0005] These considerations can be taken into account during the energy targeting phase prior to design, or alternatively, during the reconfiguration or refitting of the plant or equipment. It is preferable to consider these optimization issues with a modeling system prior to the actual design, redesign, construction or modification of the actual plant and equipment.

[0006] The state-of-the-art software currently on the market includes AspenTech Inc. software known as Aspen Pinch, Hyprotech Inc. software known as HX-NET (acquired by AspenTech), Pinch Express of KBC and Sprint of UMIST. In the targeting phase, these software products allow specific stream conditions of a process to be tracked and individual operational attributes associated with these streams to be modeled and adjusted, if required. In general terms, such software products are normally employed to track the temperatures and heat capacity flows of specific material streams in a process. Although such software provides useful tools, they are not particularly flexible in application, and do not address some of the above problems systematically.

[0007] For example, in grassroots heat exchanger network design, the parameter known as global ΔT_min (minimum approach temperature) is typically used in the state-of-the-art commercial software to represent the desired level of heat recovery between hot and cold streams, such as, for example, the minimum temperature difference allowed to recover energy through a heat exchanger. A single value for a specific material stream attribute can only be adjusted at one time including the ΔT_m,_π. This forces a user of the system to employ a trial and error approach through 'tweaking' particular attributes of specific streams one at a time, to hopefully arrive at an optimized value for Q_c and/or Qh, which represent the total energy consumed for heating (Q_h) and the total energy consumed for cooling (Q_c) for the process model at optimal driving force distribution between the hot and the cold streams. This limitation becomes compounded and makes the existing software difficult to employ effectively in large- scale processes, which employ many material streams, where these material streams may have a number of operational attributes which can be modeled and adjusted. This is also the case when several sets of stream-specific minimum temperature approaches need to be analyzed for optimal driving force distribution determination in the studied energy system.

[0008] Further, process stream changes not only may result in reduction in energy utility, but also may bring reduction in ΔT_mj_n. Thus, in current state-of-the-art technology, the capital/energy trade-off in energy system synthesis must be readjusted after each process change. Recognized by the Applicants is that the change in the driving force distribution due to both stream-specific ΔT_mj_n and each process change also affects the utility level selection. The problem is interlinked and multi-dimensional in that process changes and stream-specific ΔT_mi_n selection are competing for optimal selection of utilities, optimal process conditions, and optimal energy recovery system synthesis.

[0009] Two main methods are currently in use to address such issues: mathematical programming and thermodynamic-heuristics based pinch technology. Both methods fail to solve the problem of finding optimal driving force distribution, systematically, due to both optimal process conditions and optimal stream-specific ΔT_mi_n in energy systems and without manual iteration. Optimal driving force distribution in energy systems comes from the combined effect of the system's process conditions and hot and cold stream minimum approach temperatures ΔT_mi_n. Optimal driving force distribution in energy systems can have a significant impact on energy consumption, utility selection, utility systems, and energy recovery systems capital investment. Therefore, recognized by the Applicants is that any proposed method for optimal energy recovery systems design/retrofit and optimization should address these issues systematically and without enumeration.

[0010] Currently there is no methods or program products that can handle the theoretical, practical, and economical energy targeting problems under variable driving force distribution to find optimal distribution without manual iterations (trial and error) and in a user-friendly manner.

[0011] NZ Patent No. 527,244 (July 2004) and WO Application No. 2005/010,783 (February 2005) have addressed the problem of energy targeting to find an optimal driving force distribution due to process conditions optimization, but only for a global ΔT_mj_n, and not for stream-specific ΔT_mm- Other prior teachings have suggested using heuristics to find the energy utility targets and an optimal driving force distribution at constant process conditions and stream- dependent ΔT_mi_n. Still other prior teachings have tried to find an optimal driving force distribution through process conditions optimization at a fixed global ΔT_mj_n using mathematical programming.

[0012] Recognized by the Applicants, however, is that in most industrial processes, it is at least inefficient, if not impractical, to require that all heat exchangers (and thus, all process streams and utilities) obey the same global minimum value for driving forces, since streams (and utilities) in general have very different heat transfer coefficients. Quite often, the difference in film heat transfer coefficients can be some orders of magnitude. Thus, some heat exchangers require large ΔT_min values in order to avoid requiring an excessive heat transfer area, while other units will manage well with much smaller ΔT_mj_n values. In addition, the hot stream-specific minimum approach temperature optimal set can lead to much better energy consumption targets.

[0013] When considering retrofitting, the same problems exist with still no practical solution to simultaneously finding energy targets under all possible combinations of different process conditions, while using stream-specific minimum approach temperatures (ΔT_mi_n')-

[0014] Some scientists has recognized the need for at least assigning individual contributions to the minimum driving forces for each stream and utility, based upon the heat transfer coefficient of matched streams. Typically, these ΔT_m;_n contributions not only reflect heat transfer conditions, but can be used to represent the need for expensive materials of construction, heat exchanger types, etc. These methods, however, are based upon heuristics and are iterative, lack systemization, and do not consider the possible changes in process conditions that can result in significant changes in the energy system driving force distribution. [0015] Recognized by the Applicants is that it would be beneficial to have a system, method and program product that utilizes both process conditions manipulation and stream-specific minimum approach temperatures ΔT_mi_n' to target for energy consumption, utility selection and design heat recovery systems at an optimal driving force distribution systematically, without iteration or enumeration, without customized modeling, and in a user friendly manner.

SUMMARY OF THE INVENTION

[0016] In view of the foregoing, embodiments of the present invention advantageously provide an improved method, system and program product for theoretical, practical and economical energy targeting for the design and retrofit of energy recovery systems to find global energy utility targets, establish a good starting point for energy system design/retrofit using desired levels of heat recovery, find optimal driving force distributions in the energy system which is a function of the combined effect of process conditions and stream-specific minimum approach temperatures ΔTπ_mΛ and establish high fidelity relationship between energy cost versus capital cost to design energy recovery systems systematically and without enumeration.

[0017] Specifically, embodiments of the present invention provide a system to optimize energy recovery for a process having a plurality of resource streams each having at least one operational attribute. According to an embodiment of the present invention, such a system can include an energy utility modeling computer having a processor, and memory coupled to the processor to store software and database records therein, and a database stored in the memory (volatile or nonvolatile, internal or external) of the energy modeling computer. The database can include a plurality of sets of values each separately defining a potential range of values for at least one operational attribute for each of a plurality of hot resource streams, and a plurality of sets of values each separately defining a potential range of values for at least one operational attribute for each of a plurality of cold resource streams. These attributes can include, for example, a lower and an upper boundary value for a supply temperature (Ts) of each of the plurality of hot resource streams and each of the plurality of cold resource streams, a lower and an upper boundary value for a target temperature (Tt) of each of the plurality of hot resource streams and each of the plurality of cold resource streams, a lower and an upper boundary value for a heat capacity flow rate (FCp) of each of the plurality of hot resource streams and each of the plurality of cold resource streams. The database can also include one or more sets of values including lower and upper stream specific minimum temperature approach boundary values ({ΔTmi_n ¹ [L:H]}) for each of the plurality of hot resource streams.

[0018] The system can also include energy utility modeling program product either on a separate deliverable computer readable medium, e.g., DVD, etc., or stored in the memory of the energy utility modeling computer and adapted to optimize energy recovery for the process. The energy utility modeling program product can include instructions that when executed, for example, by the energy utility modeling computer, can cause the computer to perform the operations of receiving a plurality of sets of values each defining a potential range of attribute values for a separate one of a corresponding plurality of hot resource streams; one or more sets of a plurality of values each defining a potential range of attribute values for a separate one of a corresponding plurality of cold resource streams; and/or a plurality of sets of stream specific minimum temperature approach values ({ΔT_mm'}) assignable to the plurality of hot resource streams, with at least two of the values within each set being different from each other. The operations can also include determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) and a global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]) responsive to the operations of receiving and assigning at least portions of the plurality of sets of stream specific minimum temperature approach values ({ΔT_mj_π ¹}) to the plurality of hot resource streams. Note, this calculation can be facilitated using a cascading algorithm.

[0019] The operations can also include determining from the plurality of sets of hot resource streams attribute values and from the plurality of sets of cold resource streams attribute values, a first set of discrete process conditions that provide a minimum global minimum heating energy utility (Qh(minimum)) defining a first driving force distribution, for example, with the aid of a collapsing intervals algorithm; and calculating a first energy cost requirement associated with implementing the first set of process conditions. Similarly, the operations can also include determining from the plurality of sets of hot resource streams attribute values and from the plurality of sets of cold resource streams attribute values, a second set of discrete process conditions that provide a minimum global minimum cooling energy utility (Qc(minimum)) defining a second driving force distribution, and calculating a second energy cost requirement associated with implementing the second set of process conditions. Correspondingly, the operations can, advantageously, include selecting the set of process conditions from the first and second sets of process conditions that provides a minimum energy cost requirement responsive to determining the first and the second energy cost requirements associated with the first and the second process conditions; designing a heat exchanger network responsive to the selected set of process conditions; and determining a capital cost associated with the heat exchanger network.

[0020] The operations can also include iteratively/incrementally determining the first and second sets of discrete process conditions, calculating the first and second energy cost requirements associated therewith, and selecting the set of process conditions providing the minimum energy cost requirement, for each of a plurality of combinations of stream specific minimum approach values (ΔT_mj_n') falling within the potential range of stream specific minimum approach values for each one of the plurality of hot resource streams, to thereby draw or otherwise form data for a capital versus energy cost analysis, e.g., using a Pareto optimal chart, to determine an appropriate heat exchanger network based on a trade-off between capital cost and energy cost, and to provide a decision-maker expected utility purchasing requirements based on that selection.

[0021] According to an embodiment of the present invention, provided is a method for modeling the energy consumption of a process, or cluster of processes, using a plurality of resource streams, with each resource stream having at least one but preferably at least three operational attributes. The method can include the steps of: receiving at least one set of a range of attribute values for at least one attribute of at least one resource stream used by the process, receiving at least one set of stream-specific minimum approach temperatures between the streams, and calculating at least one energy consumption value for the process using said received range or ranges of attribute values.

[0022] According to a further aspect of an embodiment of the present invention, provided is a method substantially as described above, wherein the method further includes the step of identifying all resource streams used within the process to be modeled. According to yet another aspect of an embodiment of the present invention, provided is a method substantially as described above, further including the step of identifying all operational attributes of resource streams used within the process which affect the energy consumption of the process. [0023] According to a further aspect of an embodiment of the present invention, provided is a method substantially as described above, further including the step of indicating the specific attribute value or values from the range or ranges supplied which result in the energy consumption value or values calculated. According to a further aspect of an embodiment of the present invention, provided is a method substantially as described above wherein the energy consumption value or values calculated can include global minimum energy consumed to heat resource streams, and/or global energy consumed to cool resource streams.

[0024] According to yet another aspect of an embodiment of the present invention there is provided an energy consumption modeling system that includes or incorporates a method as described above. According to yet another aspect of an embodiment of the present invention, there is provided a program product stored in a tangible computer medium to model the energy consumption of a process utilizing a plurality of resource streams, with each resource stream having at least one operational attribute; the program product being adapted to execute the steps of: receiving at least one set of a range of attribute values for at least one attribute of at least one resource stream used by the process; receiving at least one set of stream-specific minimum approach temperature between the streams; and calculating at least one energy consumption value for the process using said received at least one set of a range of attribute values.

[0025] Embodiments of an improved system, program product and method are provided which allow for flexibility in the ability of a user to input potential ranges of values or settings for the operational attributes of material streams within an industrial process. In an embodiment of the present invention, the improved system, program product and method also allows the calculation of the global minimum and maximum energy consumed for heating (Q_h) and the global minimum and maximum energy consumed for cooling (Q_c) for the process model at an optimal driving force distribution between the hot and the cold streams in the energy system systematically and without enumeration. An embodiment of the present invention also allows for defining the optimal driving force distribution that renders the global optimum of heating or cooling utility.

[0026] Embodiments of the present invention relate to methods and user-friendly program products to calculate global energy utility targets and define an optimal driving force distribution for a process or cluster of processes under all possible process changes and stream specific minimum temperature approaches, simultaneously, and without enumeration. These embodiments can utilize stream-specific minimum temperature approaches ΔT_mi_n', where the superscript i represents the specific hot stream, as the optimization parameters, instead of the single global ΔT_mi_n currently used by all the state-of-the-art software, in addition to identifying the optimal operating conditions. Embodiments of the program product define optimal process conditions, an optimal driving force distribution in heat recovery systems, and an optimal pareto- curve that shows the rigorous trade off between energy cost and capital cost for any energy recovery system, without subjective assumptions such as, for example, systems depreciation time and future value of money.

[0027] Embodiments of the present invention can have several commercial benefits. For example, embodiments of the present invention allow for application of a unique advanced methodology automated in a user-friendly program product to optimally design and efficiently operate process units and energy recovery system in industrial facilities. Industrial companies can improve energy efficiency consumption and pollution minimization points of view in designing and operating their facilities. An estimated 10% improvement in energy efficiency due to the application of embodiments of this invention, which represents optimization beyond what is currently available from the state-of-art tools and technology, can result in savings estimated in the tens of millions of dollars per year in energy consumption and project capital.

BRIEF DESCRIPTION OF THE DRAWINGS

[0028] 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.

[0029] FIG. 1 is a schematic block diagram of a system to determine global energy utility targets and to define an optimal driving force distribution for a process having a plurality of resource streams according to an embodiment of the present invention;

[0030] FIGS. 2A-B is a schematic flow diagram of an energy utility modeling program product illustrating high-level functions according to an embodiment of the present invention;

[0031] FIGS. 3-9 are graphs and tables illustrating current methods used in Pinch technology;

[0032] FIG. 10 is a graph illustrating assignment of a specific (distinct) minimum approach temperature (ΔT_mi_n') to each hot resource stream of a process in a calculation for finding a global heating energy utility target requirement (Q_h) and a global cooling energy utility target requirement (Q_c) according to an embodiment of the present invention;

[0033] FIG. 11 is a block flow diagram illustrating the temperature steps of a cascading algorithm used in a calculation for finding a global minimum heating energy utility requirement (Q_h) and a global minimum cooling energy utility requirement (Q₀), according to an embodiment of the present invention;

[0034] FIGS. 12 and 13 are comparative tables illustrating how Q_h and Q_c differ according to the selection of the minimum approach temperature ΔT_mi_n' for hot resource streams Hl and H2 according to an embodiment of the present invention; [0035] FIGS. 14-16 are comparative tables illustrating the effect of combining changes in heat capacity flowrate (FC_P) with different minimum approach temperatures ΔT_mi_n' for hot resource streams Hl and H2 according to an embodiment of the present invention;

[0036] FIG. 17 is a table illustrating comparative values of the Q_h and Q_c targets under all possible combinations of FC_P changes for minimum approach temperatures ΔT_mi_n' between [5: 15] ⁰K of for the Hl and H2 hot resource streams according to an embodiment of the present invention;

[0037] FIG. 18 is a table illustrating the Q_h, and Q_c targets under variable supply temperature (T_S), target temperature output (T_t), and FC_P conditions, but for same approach temperature of 10⁰K for both Hl and H2 hot resource streams according to an embodiment of the present invention;

[0038] FIG. 19 is a table illustrating the Q_h and Q₀ targets under variable process conditions and variable minimum approach temperatures ATm_Jn ¹ according to an embodiment of the present invention;

[0039] FIG. 20 is a schematic diagram of an exemplary graphical user interface according to an embodiment of the present invention;

[0040] FIG. 21 is a graph illustrating a pareto-curve for the objectives of minimum energy cost and minimum capital cost generated according to an embodiment of the present invention;

[0041] FIG. 22 is a graphical user interface illustrating the temperature steps of the first phase of a cascading algorithm implemented according to embodiments of the present invention; and

[0042] FIG. 23 is a graphical user interface illustrating the temperature steps of the second phase of a cascading algorithm implemented according to embodiments of the present invention. DETAILED DESCRIPTION

[0043] 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.

[0044] As shown in FIGS. 1-23, embodiments of the present invention provide a method, system and/or program product to be used to model the energy consumption of an industrial process to optimize heat/energy recovery within a process having a plurality of resource streams, through the determination of global energy utility targets and the determination of an optimal driving force distribution for the process. Those skilled in the art should appreciate that embodiments of the present invention may encompass specific hardware or apparatus used to implement the present invention in addition to a computer program product programmed into programmable logic or digital devices adapted to execute to a number of processing steps to achieve the aims of the invention.

[0045] As noted previously, and is perhaps best shown in FIGS. 3-6, the total energy employed or consumed by the processes can be optimized to a global minimal level through, for example, careful placement and configuration of specific material streams with respect to one another. Through analysis, such as through development and utilization of a set of composite curves (see, e.g., FIGS. 3-5) which provide temperature-enthalpy profiles of the availability in a process and the demands in the process, and through what is termed a "Pinch analysis" (see, e.g., FIG. 6), an optimal driving force distribution through process conditions optimization at a global ΔT_min. (see, e.g., FIG. 7).

[0046] As perhaps best shown in FIGS. 12 and 13, the global minimum heating energy utility Qh_mi_n and the global minimum cooling energy utility Qc_min, shown as 65 KW and 75 KW, respectively, in FIG. 7 for such system, however, may not be the optimal global minimum heating utility and/or global minimum cooling because not all heat exchangers (and thus, all process streams and utilities) need to obey the same global minimum value for driving forces, since streams (and utilities) in general have very different heat transfer coefficients.

[0047] Further, as perhaps best shown in FIGS. 14-19, process conditions need not, and usually are not, specific discrete values incapable of being adjusted or analyzed separately. Additional efficiencies, i.e., improved values for global minimum heating energy utility Qh_mi_n and the global minimum cooling energy utility Qc_min can be found by utilizing both process conditions manipulation (see, e.g., FIGS. 15, 17, 19) and stream-specific minimum approach temperatures ΔT_mj_n' (see, e.g., FIGS. 12, 13, 15, 17, 19) to target for energy consumption, utility selection, and design heat recovery (e.g., heat exchange) systems, at an optimal driving force distribution.

[0048] FIG. 1 illustrates a system 30 to determine global energy utility targets and to define an optimal driving force distribution for a process having a plurality of resource streams. The system 30 can include an energy utility modeling 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 serving multiple user interfaces 37. 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.

[0049] The system 30 can also include a database 43 stored in the memory 35 (internal or external) of energy modeling computer 31 and having a plurality of sets of values each separately defining a potential range of values for at least one operational attribute for each of a plurality of hot resource streams and a plurality of sets of values each separately defining a potential range of values for at least one operational attribute for each of a plurality of cold resource streams. Such attributes can include, for example, a lower and an upper boundary value for a supply temperature (Ts) of each of the hot resource streams and each of the cold resource streams, a lower and an upper boundary value for a target temperature (Tt) of each of the hot resource streams and each of the cold resource streams, and/or a lower and an upper boundary value for a heat capacity flow rate (FCp) of each of the hot resource streams and each of the cold resource streams (see, e.g., FIG. 20).

[0050] The system 30 can also include energy utility modeling program product 51 stored in memory 35 of the energy utility modeling computer 31 and adapted to provide theoretical, practical, and economical energy targeting for the design and retrofit of energy recovery systems to find global energy utility targets, establish a good starting point for energy system design/retrofit using desired levels of heat recovery, find optimal driving force distributions in the energy system which is a function of the combined effect of process conditions and stream- specific minimum approach temperatures ΔT_mi_n' to optimize heat recovery for the process, and establish high fidelity relationship between energy cost versus capital cost to design energy recovery systems systematically and without the need for manual (trial and error) enumeration, inherent in other prior systems.

[0051 J Note, the energy utility modeling program product 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 energy utility modeling program product 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.

[0052] FIGS. 2A-B illustrate a high-level flow diagram illustrating operation of the energy utility modeling program product 51 according to an embodiment of the present invention. As shown in block 101, the program product 51 can receive input data entered by a user or stored in the database 43. The data can include, for example, a lower and an upper boundary value for a supply temperature (Ts) of each of the hot resource streams and each of the cold resource streams, a lower and an upper boundary value for a target temperature (Tt) of each of the hot resource streams and each of the cold resource streams, and/or a lower and an upper boundary value for a heat capacity flow rate (FCp) of each of the hot resource streams and each of the cold resource streams (see, e.g., FIG. 20). The data can also include a set of lower and upper stream specific minimum temperature approach ΔT_mi_n' boundary values for each the plurality of hot resource streams. Note, in the flow diagram, the "L" refers to low or minimum value of the range or interval and the "U" refers to the high or maximum value of the range or interval.

[0053] As shown in block 103, each member of the set {ΔT_mi_n'} of stream specific minimum temperature approach boundary values can be assigned the lower boundary value {ΔT_mi_n ¹}={ΔT_min¹[L]}. As shown in block 105, using the input data and the ΔT_mj_n' value data, the program product 51 can then calculate the global minimum heating energy utility interval [Qh(minimum), Qh(maximum)] illustrated more simply as Qh [L:U] with the "L" meaning low or minimum value of the interval and the "U" meaning high or maximum value of the interval. The program product 51 can also determine the global minimum cooling energy utility interval [Qc(minimum), Qc(maximum)] illustrated more simply as Qc [L:U]. According to an embodiment of the program product 51, a cascade algorithm, as will be described in more detail later, is used to perform such operation.

[0054] As shown in block 107, the program product 51 then selects as its global minimum heating energy utility value Qh the minimum value (Qh(minimum) or simply Qh[L]) of those of the global minimum heating energy utility interval values, along with its corresponding global minimum cooling energy value Qc using, e.g., a collapsing intervals algorithm, described in more detail later, to find the exact process conditions within the received range boundaries that provide such minimums. The program product 51 then calculates an associated energy cost. Similarly, as shown in block 109, the program product 51 also selects as its global minimum cooling energy utility value Qc the minimum value (Qc(minimum) or simply Qc[L]) of those of the global minimum cooling energy utility interval values, along with its corresponding global minimum heating energy value Qh using, e.g., the collapsing intervals algorithm, to find the exact process conditions within the received range boundaries that provide such minimums; and calculates an associated energy cost. As shown in block 111, the program product 51 then selects the minimum energy cost of the two calculated energy costs.

[0055] As shown in block 113, using the specific minimum temperature approach AT_mj_n' values providing the process conditions which provide the minimum energy cost, a heat exchanger network (HEN) can then be designed. Further, a capital cost associated with the HEN can be calculated. As shown in block 115, the program product 51 can then use the calculated minimum energy cost along with its associated calculated capital cost to form a single point in a Pareto-optimal chart, such as, for example, that shown in FIG. 21.

[0056] As shown in block 117, if the members of the set {ΔT_mi_n'} of stream specific minimum temperature approach boundary values have not yet been assigned the upper boundary values such that {ΔT_mi_n ^l}⁼⁼{ΔT_min^l[U]}_> the analysis continues. If it has been assigned such value, the algorithm passes, or otherwise makes available, the completed Pareto curve to a user such as, for example, a decision maker, as shown in block 119, along with the corresponding definitions of the process conditions/variables, the finished sets of {ΔTmin¹}, and associated HEN design, as shown in block 121.

[0057] As shown in blocks 123 and 125, and as will be described in more detail later, incremental values are marched into ΔTmi_n ¹ in order to complete an analysis on the continuum of values. As shown in block 127, similar to block 105, but using the incremented ΔT_mi_n' value data, the program product 51 then calculates the global minimum heating energy utility interval Qh [L:U] and the global minimum cooling energy utility interval Qc [L:U], using, for example, the cascade algorithm. As shown in block 129, the global minimum heating utility interval, global minimum cooling utility interval, and set {ΔTmi_n ¹} of stream specific minimum temperature approach values, are then stored in, for example, database 43.

[0058] As shown in block 131, if the number "i" of hot streams (i) is not the maximum number, the program product 51 resets the set {ΔTmin¹} of stream specific minimum temperature approach values (block 133) increments the value of "i" and increments ΔTmin¹ (block 125) to continue the analysis and storage of datasets. If the number "i" of hot streams (i) is the maximum number, as shown in block 137, the program product 51 chooses the -set {ΔT_mi_n ¹} with the minimum increase in the more expensive global minimum energy utility values Qh[L] or Qc[L]. As shown in block 105, the program product 51 then calculates the global minimum heating energy utility interval Qh [L:U] and the global minimum cooling energy utility interval Qc [L:U] using this chosen set {ΔT_mi_n ¹} .

[0059] Additional details describing some of the more important functions of the program product 51, according to an embodiment of the present invention, follow: [0060] FIG. 20 is a screen shot of a graphical user interface for a program product, which illustrates an industrial process that incorporates four separate and distinct resource streams in an embodiment of the present invention. Resource streams Hl and H2 are hot streams, whereas streams Cl and C2 are cold streams. Resources streams Hl and H2, each have their own independent minimum temperature approach with reference to cold streams Cl and C2. Ih the embodiment shown in FIG. 20, three separate operational attributes for each resource stream are monitored and modelled. The supply temperature of each stream, its target temperature, and a heat capacity flow rate (FCp) for the stream are monitored and shown in FIG. 20. FIG. 20 also illustrates the situation provided when a range of attribute values have been received for each attribute of each resource stream. As this information is entered, the software continuously updates and calculates, under any possible set of minimum temperature difference between the hot and cold resources, values for both the global minimum energy required for cooling and the global minimum energy required for heating by the process employed. It also calculates both the global maximum energy required for cooling and the global maximum energy required for heating by the process employed. Aspects of the present invention have been described by way of example only and it should be appreciated that modifications and additions may be made thereto without departing from the scope thereof.

[0061] As this information is entered, the program product continuously updates and calculates, under any possible minimum temperature difference between the hot and cold resources, values for both the global minimum energy required for cooling and the global minimum energy required for heating by the process (Enthalpy Min). It also calculates both the global maximum energy required for cooling and the global maximum energy required for heating by the process (Enthalpy Max).

[0062] In a preferred embodiment, the program product provided is adapted to calculate global minimum energy consumption values required by hot and cold resource streams in any process or cluster of processes under a variable heat transfer driving force distribution, and determine for desired energy consumption targets the optimal process conditions and heat transfer driving force distribution among process streams using, for example, the mathematical formulation described below. [0063] As perhaps best shown in FIG. 10, the hot streams temperatures, received either via automated processes or by direct manual input, are shifted down one-by-one by the first set of desired minimum temperature differences,

between the hot and cold resource streams to form a set of possible discrete temperature values for a continuum of possible values for ΔT_mi_n ^ι. ΔT_mi_n ¹ represents the ΔT_mj_n of hot stream (i) which refers to the minimum temperature approach between a specific hot stream and all other cold streams. The first ΔT_mi_π' can be given a value equal to zero for all hot streams to establish the global theoretical energy targets of the energy system under a theoretical optimum driving force distribution from energy consumption point of view. The second value of AT_n,;,,¹ will be equal to, for example, one degree or some other small value as per the required temperature precision, to obtain the practical energy targets at an optimal driving force distribution from an energy point of view. This can continue until ΔT_mi_n' reaches an upper boundary. Note, the upper boundary for each hot stream ΔT_mi_n ¹ should be the stream-specific ΔTmin¹ that creates a situation of no possible integration between any specific hot stream and any cold stream, i.e., the difference between the temperature of the hot stream being analyzed and the coolest cold stream.

[0064] The shifted supply and target (output) temperatures of resource hot streams, and the actual supply and the target cold streams temperatures obtained above are then sorted in a descending order, with duplicates removed, with each successive temperature pair representing the boundaries of a "TEMPERATURE STEP" and defining a new temperature step "S". Each supply temperature and target temperature input, is in the form of intervals (e.g., a pair of range boundaries) and not single discrete numbers.

[0065] Each supply temperature and target temperature input is then divided, for example, on an equal basis according to the desired temperature precision. The completion of this step results in "N" number of process temperature steps. Note, FIGS. 9 and 11 illustrates a basic model showing the "temperature steps" for hot and cold streams having discrete attribute values, FIG. 9 illustrating such steps for an algorithm using a global ΔTmi_n of 1O⁰K (see FIG. 7) for a process having two hot resource streams Hl and H2 and two cold resource streams Cl and C2; and FIG. 11 illustrating such steps for the same process but instead using stream specific minimum temperature approaches ΔTmi_n', 5⁰K for Hl and 15⁰K for H2. [0066] The total number of temperature steps is "N+l," where S varies from 0,1, 2,..., N and the temperature step number "0" represents the external energy utility temperature step. In this step, known as the external energy step, energy output, as will be described below, is assigned two initial values: Q_s=0 ^Iow-^output= 0.0 in "energy units" and Q_s=o ^highjnitput = 0.0 in "energy units."

[0067] As perhaps best shown in FIG. 22, according to a cascading algorithm, each temperature step "S" 170 greater than 0, where S =1,2,..., N, has energy surplus Q_s ^surplus. Such energy surplus has two calculated values: Q_s ^tow-^surP^fus and Q_s ^his^h-^surP'"^s. _{It a}i_{so has energy output} q^^output g._{om one tem}p_{eramre ste}p ι₀ another. Such energy output Q_s ^output has also two calculated values: Q_s ^low-^outPut and Q₈^-⁰"**. These values for energy surplus Q_s ^suφIus and energy output Q_s ^output are calculated as follows for S=I ,2,- • -,N:

_Qs.^ow_^surplu^s _{= (} j£_FQ,J^»" -%FQ>?*)(lh, -Jb₁) ,

_Qs ^high_^s _Uφ ^lus -T_CS) ,

Q _Surpius

Q high output _ /-)Λ(gA_o«(paf i /iføgΛ Surplus

where n_s and m_s are the number of the resource hot and cold streams respectively represented in the s^th temperature step and Th₈ and Tc₈ are the higher and lower shifted temperatures for the hot streams, respectively, and actual temperature of the cold streams representing the temperature boundaries.

[0068] FCp^l0W _k : is the low value of the Heat Capacity Flowrate term resulted from the multiplication of the value of the flow F lower bound by the specific heat value Cp of the hot stream number k in flow-specific heat units.

[0069] FCp^hi9h _k : is the high value of the Heat Capacity Flowrate term resulted from the multiplication of the value of the flow F upper bound by the specific heat value Cp of the hot stream number k in flow-specific heat units. [0070] FCp'°^w _j : is the low value of the Heat Capacity Flowrate term resulted from the multiplication of the value of the flow F lower bound by the specific heat value Cp of the cold stream number j in flow-specific heat units.

[0071] FCp^hl9hj : is the high value of the Heat Capacity Flowrate term resulted from the multiplication of the value of the flow F upper bound by the specific heat value Cp of the cold stream number j in flow-specific heat units.

[0072] The global minimum heating energy utility (Qh) is then calculated as follows:

Qh(minimum) =ABS Min{ Min[ Q_s ^his^h-^outPut]₅ o.O}, where: S =0, 1, 2,..., N; and ABS is the absolute of the value between brackets.

Qh(maximum) =ABS Min{ Min[ Q^-⁰"^ O.O },

and the Qh interval becomes:

Qh= [Qh(minimum), Qh(maximum)], which represents utility, which must be obtained from outside the process.

[0073] The global minimum cooling energy utility (Qc) is then calculated using steps similar to that for calculating Qh. First, a set of intervals designated as zQh are calculated as follows:

Calculate ZQh interval

ZQh(minimum) = ABS Min{ Min[ Q_s ^his^h-^outPut] , 0.0}, where S =O₅ 1, 2,..., N;

ZQh(maximum) = ABS Min{ Min[ QJ™-™**], 0.0 }, and the ZQh interval becomes:

ZQh=[ZQh(minimum), ZQh(maximum)].

[0074] Next, as perhaps best shown in FIG. 23, according to a cascading algorithm, the calculations used to determine the Qh interval are implemented using, instead, the symmetric image of the above ZQh interval in the temperature step number "zero (0)" to represent the external hot energy utility in temperature step number "zero." This symmetric image of the ZQh interval renders the Qh interval elements, but in different order. That is, if ZQh = [100, 200], this symmetric image would be [200, 100]. Step number "zero," used in solving the Qh(minimum) and Qh(maximum) values for Qh, had two initial energy output values set as follows:

=0.0 in "energy units". For calculating Qc, we set these two initial energy output (Q_s ^output) values initially to be:

Q_s=o^low-^output=ZQh(maximum), and Q_s=o^high-^output=ZQh(minimum).

[0075] Then repeat the remaining calculations for surplus and output values of each temperature step of the cascading algorithm using these two new values instead of the "zero" values (used before), while using the left term of the interval as the Qs(low_output) and the right term of such interval as the Qs(high_output).

[0076] Upon completing all temperature intervals calculations, the output of the last interval will be the global minimum cooling utility (Qc(minimum)) and the global maximum cooling utility (Qc(maximum)).

[0077] The next steps taken after complaining execution of the cascade algorithm, i.e., completing calculations at the minimum and maximum heating and cooling utilities at minimum theoretical and then practical sets of ATm_1n ¹ based upon interval (range boundary) data, includes obtaining economic targets and calculating obtained energy cost and the corresponding heat recovery system capital cost used to draw a "Pareto" optimal graph (FIG. 21) that shows optimal energy operating cost and energy system capital cost for the one or more different driving force distributions, responsive to calculating the minimum and maximum heating and cooling utilities at the practical sets of ΔT_mi_n' values, as described above.

[0078] As perhaps best shown in FIG. 2A-B, economic targets are determined, for example, by selecting an optimal set of process conditions and ΔT_mj_n' by marching in ATm_1n ¹, for example, one-by-one degree, for each hot stream with its "stream-specific minimum approach temperature," and then using associated economic data to produce the Pareto-optimal for heat recovery system synthesis. The algorithm for analyzing economic data, according to an embodiment of the present invention, is as follows:

(1) Select desired utility as an energy objective;

(2) Collapse the process conditions interval for each specific hot stream being analyzed;

(3) Increase a ΔT_mi_n' of a specific hot stream by one degree; and (4)Determine the effect of change on the desired utility target.

[0079] Steps (3) and (4), above, are automatically repeated for all hot streams (one at a time) and define the change that results in minimum increase in desired energy consumption cost. The ΔT_mi_n' that results in such minimum increase is then assigned to be the new stream-specific ΔT_mi_n' for such hot stream, while other hot streams will continue to keep their old values.

[0080] These procedures can be repeated on sets of AT_1nJn ¹, starting from its lower bounds up to its upper bounds, to generate the Pareto-optimal graph.

Operational Conditions Interval Collapsing Algorithm:

[0081] To select the optimal set of process condition, the process conditions interval for each specific hot stream being analyzed is collapsed to render individual discrete boundary values while others are still in interval mode, to allow analysis of the specific process conditions of the subject hot stream at, for example, the range boundaries, or at some interim points therebetween, to thereby determine the desired minimum utility consumption cost and exact process conditions that lead to it. The criteria for interval collapsing, for example, is that we select the lower bound or the upper bound of the interval-based operational attribute that render global minimum of a desired energy utility target. Any operational attribute value that does not render the desired energy target will be ignored. The desired energy target can be the global minimum cooling or global minimum heating utilities or it can be the global minimum heating utilities subject to minimum cooling utilities and vice versa.

[0082] If both the lower bound value and the upper bound value of stream supply temperature or target temperature are rendering bad values for the specified energy target, i.e., the algorithm does not render the global minimum energy target, the interval-based supply temperature is tested one-by-one to select the optimal supply temperature that renders the desired energy target. The same procedures can be used for the interval-based stream target temperature.

[0083] Beneficially, insights obtained from each process application energy target calculation phase can be used to make the processing time of this algorithm consume less than a second of the computer time.

[0084] The new set of driving force distributions values and the set values extracted from the collapsed intervals of process conditions can then, for example, be fed into currently available commercial software, such as, for example, "SPRINT available through the University of Manchester, to design a suitable heat exchanger network and/or select a minimum capital cost heat exchanger network.

[0085] It is important to note that while embodiments of the present invention have been described in the context of a fully functional system, those skilled in the art will appreciate that the mechanism of the present invention and/or aspects thereof are capable of being distributed in the form of a computer readable medium of instructions in a variety of forms for execution on a processor, processors, or the like, and that the present invention applies equally regardless of the particular type of signal bearing media used to actually carry out the distribution. Examples of 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, flash drives, and other newer types of memories, and transmission type media such as digital and analog communication links.

[0086] Such media can include operating instructions related to both the method steps and the energy utility modeling program product 51, described above. For example, a computer readable medium that is readable by energy utility modeling computer 31 to optimize energy recovery for a process having a plurality of resource streams (Hl..Hn, Cl..Cn), can include instructions, that when executed, cause the computer 31 to perform the operation of receiving a plurality of sets of attribute values 161, 163, 165, (see, e.g., FIGS. 19 and 20) each defining a potential range of attribute values for a separate one of a corresponding plurality of hot resource streams (e.g., Hl, H2) and a plurality of sets of attribute values 161, 163, 165, each defining a potential range of attribute values for a separate one of a corresponding plurality of cold resource streams (e.g., Cl, C2). Particularly, the sets of attribute values 161, 163, 165, can include a lower and an upper boundary value 161 for a supply temperature (Ts) of each of the plurality of hot resource streams and each of the plurality of cold resource streams, a lower and an upper boundary value 163 for a target temperature (Tt) of each of the plurality of hot resource streams and each of the plurality of cold resource streams, and a lower and an upper boundary value 165 for a heat capacity flow rate (FCp) of each of the plurality of hot resource streams and each of the plurality of cold resource streams. The operation can also include receiving a set of lower and upper stream specific minimum temperature approach boundary values or intervals ({ΔT_mi_n'.L:!!]}), or specifically assigned stream specific minimum temperature approach values (AT_01Jn ¹) for each of the plurality of hot resource streams, e.g., 5°K for Hl, 15⁰K for H2, etc.

[0087] The operations can also include assigning one or more sets a plurality of stream specific minimum temperature approach values ({AT_m,-,,¹}) to the plurality of hot resource streams, with at least two of the values within each set of stream specific minimum temperature approach values ({ΔT_mj_n'}) being different from each other; and determining a global minimum heating energy utility interval (Oh=[Qh(minimum), Qh(maximum)]) and a global minimum cooling energy utility interval (Qc=[Qc(minimum), Qc(maximum)]) responsive to the operations of receiving and assigning.

[0088] The operations can also include determining from the plurality of sets of hot resource streams attribute values 161, 163, 165, and from the plurality of sets of cold resource streams attribute values 161, 163, 165, a first set of discrete process conditions that provide the minimum global minimum heating energy utility (Qh(minimum)) which define a first driving force distribution, and determining from the plurality of sets of hot resource streams attribute values 161, 163, 165, and from the plurality of sets of cold resource streams attribute values 161, 163, 165, a second set of discrete process conditions that provide the minimum global minimum cooling energy utility (Qc(minimum)) which define a second driving force distribution. [0089] The operations can further include calculating a first energy cost requirement associated with implementing the first set of process conditions, calculating a second energy cost requirement associated with implementing the second set of process conditions, selecting the set of process conditions from the first and second sets of process conditions that provides a minimum energy cost requirement responsive to determining the first and the second energy cost requirements associated with the first and the second process conditions, designing a heat exchanger network (HEN) responsive to the selected set of process conditions, and determining a capital cost associated with the heat exchanger network.

[0090] Also, for example, according to an embodiment of the present invention, provided is a computer readable medium that is readable by the energy utility modeling computer 31 to optimize energy recovery for a process having a plurality of resource streams, which can include instructions, that when executed, cause the computer 31 to perform the operations of receiving a plurality of sets of attribute values 161, 163, 165, (see, e.g., FIGS. 19 and 20) each defining a potential range of attribute values for a separate one of a corresponding plurality of hot resource streams (e.g., Hl, H2) and a plurality of sets of attribute values 161, 163, 165, each defining a potential range of attribute values for a separate one of a corresponding plurality of cold resource streams (e.g., Cl, C2), and receiving a plurality of sets of values ({ΔT_mj_n ¹}) ^^ach defining a potential range of stream specific minimum approach values for a separate one of the plurality of hot resource streams.

[0091] The operations can also include assigning at least one set of stream specific minimum temperature approach values ({ΔT_min¹}- {AT^'fL]} to {ΔT_mi_n'fU]}) to a plurality of hot resource streams, with at least two of the values within each set being different from each other; determining a global minimum heating energy utility (Qh(minirnum) or (Qh=Qh[L])) responsive to the operation of assigning the at least one set of stream specific minimum temperature approach values ({ΔT_mj_n'}); and determining a global minimum cooling energy utility (Qc(minimum) or (Qc=Qc[L])) responsive to the operation of assigning the at least one set of stream specific minimum temperature approach values ({ΔT_min ¹}).

[0092] The operation of determining a global minimum heating energy utility (Qh(minimum)) can also include the operation of determining a plurality of temperature step intervals 170 (see, e.g., FIG. 22) for one set of the plurality of stream specific minimum temperature approaches (AT_m1n ¹) assigned to the plurality of hot resource streams responsive to the operation of assigning at least one set of stream specific minimum temperature approach values ({ΔT_mj_n'}) to the plurality of hot resource streams; determining a minimum high output value (Q_s ^hi8h-^0UtPut) and a minimum low output value (Q,¹™-™**) _of the plurality of temperature step intervals 170; and determining a global minimum heating energy utility interval (Oh=[Qh(minimum), Qh(maximum)]) responsive to the operation of determining the minimum high output value

and the minimum low output value (Q_s'°^w-^0UtPut).

[0093] Note, as perhaps best show in FIG. 22, each temperature step interval 170 can have an input interval 171 (e.g., pair of inputs) indicating heat extracted collectively from the plurality of hot resource streams, an output interval 173

_t Q ^tά&^j»Wy) _{indicating heat} collectively applied to a plurality of cold resource streams, and an output interval 175 (Q₈SUrPiUS₁₃I-Q₈Io^W-SUrPiUS_{5 indicating surplus heat ava}ii_abie for a next of the plurality of

temperature step intervals. Note also, according to an embodiment of the present invention, the global minimum heating energy utility interval (Qh= [Qh(minimum), Qh(maximum)]) is determined according to the following calculation:

Qh= [Qh(minimum), Qh(maximum)]; where:

Qh(minimum) =ABS Min{ Min[ Q_s ^his^h-°^utPut], ₀.0};

Qh(maximum) =ABS Min{ Min[ Q_s ^low-°^utPut]₅ o.O }; where:

Q .

Q high output , s^high _Surplus .

where:

_Qsiow_^.ur_Pι_Us -Tc_s) ; and

_Qshi_Sh_suφius _{= (}j£_FC;.Λ/_gA -^FCp'°^w)(Th_s -Tc₁) .

*=1 J=I [0094] Also, according to an embodiment of the present invention, the plurality of hot resource streams can be assigned a plurality of sets of values ({ΔTmin¹ [L:U]}) each separately defining a potential range of stream specific minimum approach values (ΔTmi_n ¹) for each of the plurality of hot resource streams. The first value for the set of specific minimum temperature approaches ({ΔT_mi_n ^ι [L:U]}) utilized in the above calculations was the lowest range value (ΔT_mj_n' [L]) of each value range within the set of the plurality of stream specific minimum approach values ({ΔT_min' [L-U]}) assigned to each of the plurality of hot resource streams. Accordingly, the operations also include iteratively performing the operations of determining a plurality of temperature step intervals 170, determining a minimum high output value (Q_s ^hls^h-^outPut) Q_n^ _a minimum low output value (Q_s ^low-^Output), and determining a global minimum heating energy utility interval (Qh= [Qh(minimum), Qh(maximum)]), for each of the remaining received plurality of specific minimum temperature approach values (ΔT_min).

[0095] According to this embodiment of the computer readable medium, similar to that of determining the global minimum heating energy utility (Qh(minimum)), the operation of determining a global minimum cooling energy utility (Qc(minimum)) can include the operation of determining a global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]) which can further include the operations of first determining a symmetric image ([Qh(maximum), Qh(minimum)]) of a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]), and determining a plurality of temperature step intervals 170 for one set of the plurality of stream specific minimum temperature approaches {ΔT_mj_n'} assigned to the plurality of hot resource streams responsive to the operation of determining the symmetric image ([Qh(maxirnurn), Qh(minimum)]) which provides the initial input (see, e.g., FIG. 23) for determining ([Qc(minimum), Qc(maximum)]) using the same cascading algorithm, described previously.

[0096] The operations can also include determining a minimum high output value (Q_s ^high-^output) and a minimum low output value (Q_s ^low-°^utPut) _of each of the plurality of temperature step intervals 170, and determining the global cooling energy utility interval (Qc=[Qc(minimum), Qc(maximum)]) responsive to the determining of the minimum high output value (Q_s ^his^h-°^utPut) and the minimum low output value (Q_S ^IOW-<^>U<^PU'). The sub-operations of determining a plurality of temperature step intervals 170, determining a minimum high output value (Q_s ^Mβh-^0U|Pu«) and a minimum low output value (Q_s ^low-^output), and determining a global minimum cooling energy utility interval (Qc=[Qc(minimum), Qc(maximum)]), are then performed for each of the remaining plurality of specific minimum temperature approach values within the set {ΔTmm'}.

[0097] The operations can also further include determining from the plurality of sets of hot and cold resource streams attribute values 161, 163, 165, a first set of discrete process conditions that provide the minimum global minimum heating energy utility (Qh(minimum)) defining a first driving force distribution; and determining from the plurality of sets of hot and cold resource streams attribute values 161, 163, 165, a second set of discrete process conditions that provide the minimum global minimum cooling energy utility (Qc(minimum)) defining a second driving force distribution. The operations can further include calculating a first energy cost requirement associated with implementing the first set of process conditions, calculating a second energy cost requirement associated with implementing the second set of process conditions, selecting the set of process conditions from the first and second sets of process conditions that provides a minimum energy cost requirement responsive to determining the first and the second energy cost requirements associated with the first and the second process conditions, designing or otherwise providing instructions to design a heat exchanger network responsive to the selected set of process conditions, and determining a capital cost associated with the heat exchanger network.

[0098] The operations can also include iteratively/incrementally determining first and second sets of discrete process conditions, calculating first and second energy cost requirements associated therewith, and selecting the set of process conditions providing the minimum energy cost requirement, for each of a plurality of combinations of stream specific minimum approach values (ΔT_m,_n') falling within the potential range of stream specific minimum approach values for each one of the plurality of hot resource streams, to thereby draw a capital versus energy cost Pareto optimal chart (see, e.g., FIG. 21) used for determining an appropriate heat exchanger network based on a trade-off between capital cost and energy cost, and for providing a decision- maker expected utility purchasing requirements based on that selection.

[0099] The various embodiments of the present invention include numerous advantages and benefits. In a preferred embodiment, the program product 51 provided may also be adapted to indicate optimal cost settings or values of particular operational attributes which, when used, result in the optimal energy consumption total cost value calculated. These particular settings or values for operational attributes may then be used in the implementation of the actual process after modeling has occurred to potentially optimize the process.

[00100] Preferably the program product 51 provided can be adapted to model the energy consumption of a particular process. In general terms, the modeling operation or facility provided can be used in investigations or analyses of specific processes and their energy consumption characteristics or requirements. In a further preferred embodiment, the model or modeling facility provided may be used to optimize the energy consumption characteristics of a process, preferably to minimize the amount of energy consumed. Reference throughout this specification was also made to embodiments of the present invention being used as a tool to optimize the energy consumption of a process, but those skilled in the art should appreciate that embodiments of the present invention may also be used as a simple analytical or analysis tool, if required.

[00101] Furthermore, a process as referred to throughout this specification may be defined as any operation or set of operations which consumes energy to in turn produce a required result, or alternatively to produce a material or product. Processes can range, for example, from the domestic with the use of household hot water system to heat, cool or mechanically process or chop materials, to the industrial, to include electricity generation systems and plant as examples.

[00102] Processes that have their associated energy consumption modeled in conjunction with embodiments the present invention, may also include manufacturing plant, equipment or manufacturing production lines, or collections of production lines which may produce unrelated products or materials. Furthermore, a single process to be considered in conjunction with embodiments of the present invention may also consist of or incorporate a set of sub-processes which need not necessarily be related together, nor in some instances, need not require the physical apparatus or equipment employed to be located in close proximity or adjacent to other equipment used.

[00103] Reference throughout this specificationwas, however, made to a process modeled in conjunction with embodiments of the present invention being a single manufacturing production line for a specific type of product. However, those skilled in the art should appreciate that these references are made for the sake of simplicity only and the embodiments of the present invention may be applied to a wide range and number of different types of varied process, as discussed above.

[00104] Embodiments of the present invention may also be adapted to model energy consumption of a single process. The energy consumed may be in a variety of different forms of mechanical and thermal energy. However, for the sake of simplicity, reference throughout this specification will be made to the energy consumed being originally sourced from thermal energy. Again, those skilled in the art should appreciate that other types of energy may also be monitored or modeled in conjunction with the present invention and reference to the above only throughout this specification should in no way be seen as limiting.

[00105] In addition, the energy consumed in conjunction with embodiments of the present invention were, in general terms, referred to as being employed to heat or cool materials. Again, however, those skilled in the art should appreciate that other energy consumption applications may also be considered in conjunction with the embodiments of the present invention and reference to the above only should in no way be seen as limiting.

[00106] Preferably a process modeled in accordance with embodiments of the present invention may include a plurality of resource streams. In general terms, a resource stream may consist of a flow of material that is consumed, modified or employed in some way within the process. For example, resource streams may consist of flows of water or other types of fluid or gases, or alternatively two or three phase flows. Pressure as a state variable can be handled the same way energy and mass quantities are handled.

[00107] Reference throughout this specification was also made to the resource streams monitored using embodiments of the present invention and employed in the particular process modeled to include flows of one or more types of fluid. However, those skilled in the art should appreciate that other types of resources may also be employed and modeled in conjunction with embodiments of the present invention and reference to the above only throughout this specification should in no way be seen as limiting. Furthermore, those skilled in the art should also appreciate that a resource stream need not necessarily be made up of a substantially continuous flow of material. For example, in some instances a resource stream may be composed from a supply or flow or material delivered periodically in discreet blocks, if required.

[00108] In a preferred embodiment of the present invention, a resource stream may have at least one operational attribute associated with it. An operational and/or design attribute may be a specific characteristic or parameter associated with a particular stream which in turn will have an effect on the energy consumed by the process modeled. Those skilled in the art should appreciate that various different types of operational attributes may be monitored and modeled in conjunction with the embodiments of the present invention, depending on the particular type of resource stream which the operational attribute is associated with. For example, in some instances, examples of operational attributes could include a materials temperature, flow rate, composition, pressure, latent heat of vaporization, latent heat of condensation, latent heat of sublimation, and/or combinations of all these attributes.

[00109] For the sake of simplicity and in general terms, reference throughout this specification was made to the operational attributes monitored and modeled being the temperature of a fluid flow, the flow rate, and the specific heat of the resource stream of a fluid flow. Again however, those skilled in the art should appreciate that other configurations and implementations of the present invention are envisioned and reference to the above only throughout this specification should in no way be seen as limiting.

[00110] In a preferred embodiment of the present invention, the program product employed can record, hold or otherwise receive a list of identified resource streams within a process that affect the energy consumption of the process. Preferably all the resource streams involved may be monitored and modeled in conjunction with the embodiments of the present invention to provide the most accurate indication or calculation of the energy consumption value or values required.

[00111] hi alternative embodiments of the present invention, it is not necessary to identify and monitor all resource streams employed in a process. For example, in such an alternative embodiment, only a selected number of streams which are of the most importance or have the most significance with respect to energy consumption of the process may be monitored or identified. This can provide an approximate calculation of the energy consumption value or values required while using relatively low cost computer systems or to provide a fast approximate indication of an energy consumption value.

[00112] In a preferred embodiment of the present invention, all relevant operational and/or design attributes of the resource streams monitored or modeled may in turn be identified for or within the program product employed. Preferably all operational and/or design attributes of the streams monitored which have an effect on the energy consumption of a process may be tracked and modeled in conjunction with a preferred embodiment of the present invention. Modeling each and every operational and design attribute which affects energy consumption will in turn provide an accurate calculation of energy consumption for the process modeled.

[00113] In alternative embodiments of the present invention, each and every relevant operational and/or design attribute may not necessarily be monitored or modeled. For example, in one alternative embodiment, only a selected number of highly relevant or highly significant attributes with respect to energy consumption may be monitored or modeled in conjunction with the present invention. This program product allows the rigorous check of such significance because of its interactive data entry capability that enables the calculation of the energy consumption global minima and maxima upon the entry of the data of each resource stream. This approach can allow a fast and approximate energy consumption values to be calculated or alternatively allow the present invention to function using a relatively low capacity computer system or systems.

[00114] In a preferred embodiment of the present invention, the program product provided is adapted to receive at least one set of a range of attribute values from a user. The user of the program product may be a person, persons or organization wishing to model and potentially optimize the energy consumption of a particular process. The user involved can input at least one set of a range of attribute values for a selected operational or design attribute of a resource stream to in turn provide the base data or information required to set up a model of the process involved. [00115] In a preferred embodiment of the present invention, at least one set of a range of attribute values may be received for one or more operational attributes of a resource stream. This information can be indicative of an allowable range of values for a particular operational and/or design attribute that can vary the energy consumed by the process involved. Beneficially, by supplying a range of attribute values, this substantially simplifies the use of the program product provided for a user, who need not experiment with each and every potential case provided through all the discreet single values capable of being entered and received.

[00116] In a further preferred embodiment of the present invention, a single set of a range of attribute values may be received from a user and assigned to each and every operational and design attribute identified for each and every resource stream monitored and modeled in conjunction with the present invention. The range of attribute values provided can give the operational and design parameters within which each of the operational attributes can vary or fluctuate in the running of a process to be modeled as well as stream-specific minimum approach temperatures.

[00117] However, in alternative embodiments of the present invention, a single range of attribute values need not necessarily be received for any each and every operational and/or design attribute monitored and modeled in conjunction with the present invention. For example, in some alternative embodiments, a single range of operational or design attributes may be provided for a single operational or design attribute for a resource stream that has multiple attributes, as should be appreciated by those skilled in the art. Furthermore, it should also be appreciated that more than one set of ranges of attribute values may be received in relation to a single operational attribute. Several sets of continuous ranges of attribute values may be received, if required, in accordance with such an embodiment.

[00118] Reference throughout this specification was, however, made to the program product employed receiving a single set of a range of attribute values for each and every attribute of each and every resource stream monitored and modeled in conjunction with the present invention. However, those skilled in the art should appreciate that other configurations and implementations of the program product to be provided are envisioned and reference to the above only throughout this specification should in no way be seen as limiting. [00119] In a preferred embodiment of the present invention, the program product provided is adapted to calculate at least one energy consumption value using the received range or ranges of attribute values. Preferably the energy consumption value or values calculated may be an optimal value for the process involved, calculated through determining or selecting a specific collection of attribute values, which result in an optimized energy consumption result.

[00120] Ih a further preferred embodiment of the present invention, the program product employed may calculate the global minimum energy consumed by a process or a specifically optimized set of operational and/or design attribute values or settings. This calculation may then be used to fine tune the operation or design of equipment to be employed to implement the process or alternatively to adjust or assist in redesign and redevelopment or refitting work to be completed with respect to existing process system.

[00121] In a further preferred embodiment of the present invention, the program product employed may be used to calculate a value for Q_h, being the global energy utility employed to heat resource streams, and Q_c, being the global energy utility used to cool resource streams. It can also be used to calculate minimum heating utility subject to global minimum cooling utility and the minimum cooling utility subject to global minimum heating utility.

[00122] In a preferred embodiment of the present invention, the program product provided is adapted to calculate global minima energy consumption values required by hot and cold resource streams in any process or cluster of processes using unique mathematical formulation.

[00123] Those skilled in the art should appreciate that other implementations or aspects of the technology developed are also envisioned and reference to the above only throughout this specification should in no way be seen as limiting. Hence the invention is the methodology, the mathematical algorithm and the user-friendly program product that automatically execute the methodology. The calculation engine/mathematical algorithm in the program product is also another component of the invention since it can be utilized by others to develop their own program product.

[00124] This non-provisional application is related to US Patent Application No. 60/816,234 incorporated herein by reference in its entirety. [00125] 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. This invention is not to be construed as limited to the particular forms or embodiments disclosed, since these are regarded as illustrative rather than restrictive. 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. For example, in the exemplary embodiment, many points of the heat exchange analysis related to ΔT_mm¹ focused on hot resource streams with respect to cold resource streams. Alternatively, the heat exchange analysis could be directed to cold resource streams with respect to hot resource streams, and still be within the scope of the present invention.

Claims

THAT CLAIMED IS:

1. A system (30) to optimize energy recovery for a process having a plurality of resource streams each having at least one operational attribute (161, 163, 165), the system (30) comprising an energy utility modeling computer (31) having a processor (33), and memory (35) coupled to the processor (33) to store software and database records therein, the system (30) being characterized by: a database stored in the memory (35) of the energy modeling computer (31) and having a plurality of sets of values each separately defining a potential range of values for at least one operational attribute (161, 163, 165) for each of a plurality of hot resource streams (Hl..Hn) and a plurality of sets of values each separately defining a potential range of values for at least one operational attribute (161, 163, 165) for each of a plurality of cold resource streams (Cl..Cn); and energy utility modeling program product (51) stored in the memory (35) of the energy utility modeling computer (31) to optimize energy recovery for the process, the program product (51) including instructions that when executed by the energy utility modeling computer (31), cause the computer (31) to perform the operations of: receiving a plurality of sets of values each defining a potential range of attribute values (161, 163, 165) for a separate one of a corresponding plurality of hot resource streams (H L. Hn) and a plurality of sets of values each defining a potential range of attribute values (161, 163, 165) for a separate one of a corresponding plurality of cold resource streams (C L. Cn); assigning at least one set of a plurality of stream specific minimum temperature approach values ({ΔT_mi_π ¹}) to the plurality of hot resource streams (Hl..Hn), at least two of the values (ΔTmin¹) within each set of stream specific minimum temperature approach values ({ΔT_mi_n'}) being different from each other; and determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) and a global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]) responsive to the operations of receiving and assigning.

2. A system (30) as defined in Claim 1, wherein the energy utility modeling program product (51) further includes instructions to perform the operations of: determining from the plurality of sets of hot resource streams attribute values (161, 163, 165) and from the plurality of sets of cold resource streams attribute values (161, 163, 165), a first set of discrete process conditions that provide a minimum global minimum heating energy utility (Qh(minimum)) defining a first driving force distribution; and determining from the plurality of sets of hot resource streams attribute values (161, 163, 165) and from the plurality of sets of cold resource streams attribute values (161, 163, 165), a second set of discrete process conditions that provide a minimum global minimum cooling energy utility (Qc(minimum)) defining a second driving force distribution.

3. A system (30) as defined in Claim 2, wherein the energy utility modeling program product (51) further includes instructions to perform the operations of: calculating a first energy cost requirement associated with implementing the first set of process conditions; calculating a second energy cost requirement associated with implementing the second set of process conditions; and selecting the set of process conditions from the first and second sets of process conditions that provides a minimum energy cost requirement responsive to determining the first and the second energy cost requirements associated with the first and the second process conditions.

4. A system (30) as defined in Claim 3, wherein the energy utility modeling program product (51) further includes instructions to perform the operations of: designing a heat exchanger network responsive to the selected set of process conditions; and determining a capital cost associated with the heat exchanger network.

5. A system (30) as defined in Claim 4, wherein the energy utility modeling program product (51) further includes instructions to perform the operation of receiving a plurality of sets of values each separately defining a potential range of stream specific minimum approach values ({ΔT_mi_n [L:U]}) for each of the plurality of hot resource streams (Hl..Hn).

6. A system (30) as defined in Claim 5, wherein the operations of determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) and a global minimum cooling utility interval ([Qc(minimum), Qc(maximum)]), determining first and second sets of discrete process conditions, calculating first and second energy cost requirements associated therewith, and selecting the set of process conditions providing the minimum energy cost requirement, are performed incrementally for each of a plurality of combinations of stream specific minimum approach values falling within the potential range of stream specific minimum approach values ({AlWlI/U]}) for each one of the plurality of hot resource streams (Hl..Hn).

7. A system (30) as defined in any of Claims 1-6, wherein the operation of determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) and a global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]) is performed using a cascade algorithm.

8. A system (30) as defined any of Claims 1-7, wherein the operation of receiving includes the operations of receiving a lower and an upper boundary value (161) for a supply temperature (Ts) of each of the plurality of hot resource streams (Hl.. Hn) and each of the plurality of cold resource streams (Cl..Cn)₅ receiving a lower and an upper boundary value (163) for a target temperature (Tt) of each of the plurality of hot resource streams (Hl.. Hn) and each of the plurality of cold resource streams (C 1..Cn), receiving a lower and an upper boundary value (165) for a heat capacity flow rate (FCp) of each of the plurality of hot resource streams (Hl..Hn) and each of the plurality of cold resource streams (C 1..Cn); and receiving a set of lower and upper stream specific minimum temperature approach boundary values (ΔT_mm'P^U]) for each of the plurality of hot resource streams (H 1..Hn).

9. A system (30) as defined in any of Claims 1-8, wherein the operation of determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) includes the operations of: determining a plurality of temperature step intervals (170) for one set of the plurality of stream specific minimum temperature approach values ({ΔT_min'}) assigned to the plurality of hot resource streams (Hl..Hn) responsive to the operation of assigning the plurality of streams specific minimum temperature approach values (AT_m1n ¹), each temperature step interval (170) having an input interval (171) indicating heat extracted collectively from the plurality of hot resource streams (Hl..Hn), an output interval (173, Q_s ^low-^0UtPut , Q WJ^MW*) indicating heat collectively applied to the plurality of cold resource streams (C L. Cn), and an output interval (175, Q_s ^low-^surP'^us , Q_s ^high_^surpiu _{S) indicating surplus heat} available for a next of the plurality of temperature step intervals (170); determining a minimum high output value (min[Q_s ^hiεh-^0Utput]) and a minimum low output value (min[Q_s ^l0W-^0Utput]) of the plurality of temperature step intervals (170); and determining the global minimum heating energy utility interval ([Qh(minimum), Qh(rnaximum)]) responsive to the operation of determining the minimum high output value (min[Q_s ^high-^0Utput]) and the minimum low output value (min[Q_s ^low-^output]).

10. A system (30) as defined in Claim 9, wherein the one set of the plurality of specific minimum temperature approach values ({ΔTmin¹}) includes the lowest range value ({ΔTmin'[L]}) of each value range ({ΔTmnV[L:U]}) within the set of the plurality of stream specific minimum approach values ({ΔTmin¹}) assigned to each of the plurality of hot resource streams (H 1..Hn); wherein the energy utility modeling program product (51) further includes instructions to perform the operation of receiving a plurality of sets of values (ΔTmin'[L:U]) each separately defining a potential range of stream specific minimum approach values (ΔTmin¹) for each of the plurality of hot resource streams (H 1..Hn); and wherein the operations of determining a plurality of temperature step intervals (170), determining a minimum high output value (min[Q_s ^hl8h-^0Utput]) and the minimum low output value (min[Q_s ^l0W-^0Utput]), and determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) are performed for each of the remaining received plurality of sets of specific minimum temperature approach values ({ΔTmin¹}).

11. A system (30) as defined in Claim 10, wherein a first value for the specific minimum temperature approach ({ΔTmin¹}) used by the first temperature step (170) is initially set to zero for all hot resource streams (H 1..Hn) to establish a global theoretical energy target values (Qh, Qc) for the process.

12. A system (30) as defined in any of Claims 9-11, wherein the operation of determining the global minimum cooling energy utility interval (Qc(minimum), Qc(rnaximum)) includes the operations of: determining a symmetric image ([Qh(maximum)₅ Qh(mirdmum)]) of the global minimum heating energy utility interval ([Qh(minimum), Qh(rnaximurn)]); determining a plurality of temperature step intervals (170) for one set of the plurality of stream specific minimum temperature approach values ((AT_m1n ¹)) assigned to the plurality of hot resource streams (H 1..Hn) responsive to the operation of determining the symmetric image ([Qh(maximum), Qh(minimum)]), each temperature step interval (170) having an input (171) indicating heat extracted collectively from the plurality of hot resource streams (Hl..Hn), an output (173, Q_s ^low-^outP^ut , Q₅ ^WgIu^)UtP^Ut) _{indicating heat} collectively applied to the plurality of cold resource streams (CL.Cn), and an output (175,

_? Q_s ^high__S ^uπ>ⁱu_{S) indicating surplus heat} available for a next of the plurality of temperature step intervals (170); determining a minimum high output value (min[Q_s ^hlsh-^0Utput]) and a minimum low output value (min[Q_s ^l0W-^0Utput]) of the plurality of temperature step intervals (170); and determining the global cooling energy utility interval ([Qc(minimum), Qc(maximum)]) responsive to the determining of the minimum high output value (mm[Q_s ^hl8h-^0Utput]) and the minimum low output value (min[Q_s ^low-^output]).

13. A system (30) as defined in Claim 12, wherein the sub-operations of determining a plurality of temperature step intervals (170), determining a minimum high output value (min[Q_s ^high-^output]) and a minimum low output value (min[Q_s ^Iow-^output]), and determining a global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]) responsive thereto, of the operation of determining the global minimum cooling energy utility interval ([Qc(minirnum), Qc(maximum)]), are performed for each of a plurality of sets of specific minimum temperature approach values ({ΔTmin¹}).

14. A system (30) to optimize energy recovery for a process having a plurality of resource streams (Hl..Hn, Cl..Cn) each having at least one operational attribute (161, 163, 165), the system (30) comprising an energy utility modeling computer (31) having a processor (33), and memory (35) coupled to the processor (33) to store software and database records therein, the system (30) being characterized by: energy utility modeling program product (51) stored in the memory (35) of the energy utility modeling computer (31) to optimize energy recovery for the process, the program product (51) including instructions that when executed by the energy utility modeling computer (31), cause the computer (31) to perform the operations of: assigning at least one set of stream specific minimum temperature approach values ({ΔT_mi_n'}) to a plurality of hot resource streams (Hl.. Hn), at least two of the values (AT_01Jn ¹) within each at least one set ({ΔTmi_n ¹}) being different from each other; and determining a global minimum heating energy utility (Qh(minimum)) responsive to the operation of assigning the at least one set of stream specific minimum temperature approach values ({ΔT_mi_n'}).

15. A system (30) as defined in Claim 14, wherein the energy utility modeling program product (51) further includes instructions to perform the operations of: determining a global minimum cooling energy utility (Qc(minimum)) responsive to the operation of assigning the at least one set of stream specific minimum temperature approach values ((AT_nJ)); determining from a plurality of sets of hot resource streams attribute values (161, 163, 165) and from a plurality of sets of cold resource streams attribute values (161, 163, 165), a first set of discrete process conditions that provide the minimum global minimum heating energy utility (Qc(minimum)) defining a first driving force distribution; and determining from the plurality of sets of hot resource streams attribute values (161, 163, 165) and from the plurality of sets of cold resource streams attribute values (161, 163, 165), a second set of discrete process conditions that provide the minimum global minimum cooling energy utility (Qc(minimum)) defining a second driving force distribution.

16. A system (30) as defined in Claim 15, wherein the energy utility modeling program product (51) further includes instructions to perform the operations of: calculating a first energy cost requirement associated with implementing the first set of process conditions; ■ calculating a second energy cost requirement associated with implementing the second set of process conditions; and selecting the set of process conditions from the first and second sets of process conditions that provides a minimum energy cost requirement responsive to determining the first and the second energy cost requirements associated with the first and the second process conditions.

17. A system (30) as defined in Claim 16, wherein the energy utility modeling program product (51) further includes instructions to perform the operations of: designing a heat exchanger network responsive to the selected set of process conditions; and determining a capital cost associated with the heat exchanger network.

18. A system (30) as defined in either of Claims 15 or 16, wherein the energy utility modeling program product (51) further includes instructions to perform the operation of receiving a plurality of sets of values (ΔT_min'[L:U]) each defining a potential range of stream specific minimum approach values (ΔTmin¹) for a separate one of the plurality of hot resource streams (Hl.. Hn).

19. A system (30) as defined in Claim 18, wherein the operations of determining a global minimum heating energy utility (Qh(minimum)) and a global minimum cooling utility (Qc(minimum)), determining first and second sets of discrete process conditions, calculating first and second energy cost requirements associated therewith, and selecting the set of process conditions providing the minimum energy cost requirement, are performed incrementally for each of a plurality of combinations of stream specific minimum approach values (ΔT_mi_π') falling within the potential range of stream specific minimum approach values (ΔTm_Jn'tL:!!]) for each one of the plurality of hot resource streams (Hl..Hn).

20. A system (30) as defined in any of Claims 15-19, wherein the operation of determining a global minimum heating energy utility (Qh(minimum)) and a global minimum cooling energy utility (Qc(minimum)) is performed using a cascade algorithm.

21. A system (30) as defined in any of Claims 14-20, wherein the energy utility modeling program product (51) further includes instructions to perform the operation of receiving a plurality of sets of values (161, 163, 165) each defining a potential range of attribute values (161, 163, 165) for a separate one of a corresponding plurality of hot resource streams (Hl..Hn) and a plurality of sets of values (161, 163, 165) each defining a potential range of attribute values (161, 163, 165) for a separate one of a corresponding plurality of cold resource streams (Cl..Cn).

22. A system (30) as defined in Claim 21, wherein the operation of receiving includes the operations of receiving a lower and an upper boundary value (161) for a supply temperature (Ts) of each of the plurality of hot resource streams (H L. Hn) and each of the plurality of cold resource streams (C 1..Cn), receiving a lower and an upper boundary value (163) for a target temperature (Tt) of each of the plurality of hot resource streams (Hl..Hn) and each of the plurality of cold resource streams (Cl..Cn), receiving a lower and an upper boundary value (165) for a heat capacity flow rate (FCp) of each of the plurality of hot resource streams (Hl..Hn) and each of the plurality of cold resource streams (Cl..Cn); and receiving a set of lower and upper stream specific minimum temperature approach boundary values (ΔT_mj_n'[L:U]) for each of the plurality of hot resource streams (Hl.. Hn).

23. A system (30) as defined in any of Claims 14-22, wherein the operation of determining a global minimum heating energy utility (Qc(minimum)) includes the operations of: determining a plurality of temperature step intervals (170) for one set of the plurality of stream specific minimum temperature approach values (ΔT_mi_n ^ι) assigned to the plurality of hot resource streams (Hl..Hn) responsive to the operation of assigning at least one set of stream specific minimum temperature approach values ({ΔT_mi_n'}) to the plurality of hot resource streams (Hl..Hn), each temperature step interval (170) having an input interval (171) indicating heat extracted collectively from the plurality of hot resource streams (H 1..Hn), an output interval (173, Q_s ^low-^outP^ut _t Q_s ^high_^output) _indicati_ng _heat collectively applied to a plurality of cold resource streams (Cl..Cn), and an output interval (175, Q_s ^low-^sutP'^us ,

indicating surplus heat available for a next of the plurality of temperature step intervals (170); determining a minimum high output value (min[Q_s ^hl8h-^0Utput]) and a minimum low output value (min[Q_s ^l0W-^0Utput]) of the plurality of temperature step intervals (170); and determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) responsive to the operation of determining the minimum high output value (min[Q_s ^high-^0Utput]) and the minimum low output value (min[Q_s ^l0W-^output]).

24. A system (30) as defined in Claim 23, wherein the one set of the plurality of specific minimum temperature approach values ({ΔTmin¹}) includes the lowest range value ({ΔTmin'[L]}) of each value range within the set of the plurality of stream specific minimum approach values ({ΔTmin¹}) assigned to each of the plurality of hot resource streams (Hl..Hn); wherein the energy utility modeling program product (51) further includes instructions to perform the operation of receiving a plurality of sets of values (AT_m1n ¹ [L:U]) separately defining a potential range of stream specific minimum approach values (ΔTmi_n ¹) for each of the plurality of hot resource streams (Hl..Hn); and wherein the operations of determining a plurality of temperature step intervals (170), determining a minimum high output value (min[Q_s ^hlgh-^0Utput]) and a minimum low output value (min[Q_s ^low-^output]), and determining a global minimum heating energy utility interval ([Qh(minimurn), Qh(maximum)]) are performed for each of a plurality of sets of specific minimum temperature approach values ({ΔTmin¹}).

25. A system (30) as defined in Claim 24, wherein a first value for the specific minimum temperature approach (ΔTmin¹) used by the first temperature step is initially set to zero for all hot resource streams (Hl..Hn) to establish a global theoretical energy target values for the process.

26. A system (30) as defined in any of Claims 23-25, wherein the global minimum heating energy utility interval is determined according to the following calculation:

Qh= [Qh(minimum), Qh(maximum)]; where: Qh(minimum) =ABS Min{ Min[ Q_s ^his^h-^0UtPut]₅ o.O}; Qh(maximum) =ABS Min{ Min[ Q_s ^l0W-^0UtPut], o.O } ; where:

Q lew output _ /-dow output . s^low Surplus , s ~ ~y_s-ι ⁺y_s >

Q high output _ _/-thigh output , πhigh _Surplus .

where:

Q_s ^hi9h-^surplus -Tc₁).

27. A system (30) as defined in any of Claims 23-26, wherein the energy utility modeling program product (51) further includes instructions to perform the operation of determining the global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]), which further includes the operations of: determining a symmetric image ([Qh(maximum), Qh(minimurn)]) of a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]); determining a plurality of temperature step intervals (170) for one set of the plurality of stream specific minimum temperature approach values (ΔT_mj_n') assigned to the plurality of hot resource streams (H 1..Hn) responsive to the operation of determining the symmetric image ([Qh(maximum), Qh(minimum)]), each temperature step interval (170) having an input (171) indicating heat extracted collectively from the plurality of hot resource streams (Hl.. Hn), an output (173, Q_s ^low-^0UtPut , Q_sws^h_^output) _{indicating heat c0}Uectively applied to the plurality of cold resource streams (Cl..Cn), and an output (175, Q_s ^1ow-^suπ»u^s _> Q WU^MΦ*) indicating surplus heat available for a next of the plurality of temperature step intervals (170); determining a minimum high output value (min[Q_s ^hlgh-^0Utput]) and a minimum low output value (min[Q_s ^!ow-^output]) of the plurality of temperature step intervals (170); and determining the global cooling energy utility interval ([Qc(minimum), Qc(maximum)]) responsive to the determining of the minimum high output value (min[Q_s ^hl8h-^0Utput]) and the minimum low output value (min[Q_s ^Iow-^output]).

28. A system (30) as defined in Claim 27, wherein the operations of determining a plurality of temperature step intervals (170), determining a minimum high output value (mm[Q_s ^hl8h-^output]) and a minimum low output value (min[Q_s ^low-^output]), and determining a global minimum cooling energy utility interval of the operation of determining the global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]), are performed for each of a plurality of sets of specific minimum temperature approach values ({ΔTmin¹}).

29. A method for modeling the energy consumption of a process, said process using a plurality of resource streams (Hl..Hn, Cl..Cn), with each resource stream (Hl..Hn, Cl..Cn) having at least one operational attribute (161, 163, 165), the method being characterized by the steps of: receiving at least one set of a range of attribute values (161, 163, 165) for at least one attribute (161, 163, 165) of at least one of the resource streams (Hl..Hn, Cl..Cn) used by the process; receiving at least one set of stream-specific minimum approach temperature values between streams ({ΔT_mi_n ¹}), at least two of the values (ΔT_mi_n') within the at least one set of stream specific minimum temperature approach values ({ΔT_mi_n'}) being different from each other; and calculating at least one energy consumption value for the process using the received at least one set of a range of attribute values (161, 163, 165) and the at least one set of stream- specific minimum approach temperature values ({ΔTmin¹}).

30. A method as defined in Claim 29, being further characterized by the step of identifying one or more resource streams used within the process to be modeled.

31. A method as defined in either of Claims 29 or 30, being further characterized by the step of identifying one or more operational attributes (161, 163, 165) of resource streams (Hl..Hn, Cl..Cn) used within the process which affect the energy consumption of the process.

32. A method as defined in any of Claims 29-31 , being further characterized by the step of indicating a specific attribute value (161, 163, 165) from the received range supplied which result in the desired energy consumption value or values calculated.

33. A method as defined in any of Claims 29-32, wherein the energy consumption value or values calculated are selected from the group consisting of: global minimum energy Qh(minimum) consumed to heat resource streams (Hl..Hn), global minimum energy Qc(minimum) consumed to cool resource streams (Cl..Cn), global maximum energy Qh(maximurn) consumed to heat resource streams (H 1..Hn), and global maximum energy Qc (maximum) consumed to cool resource streams (Cl..Cn).

34. A method as defined in any of Claims 29-33, wherein each resource stream (Hl..Hn, Cl..Cn) has at least three user-defined operational attributes (161, 163, 165) each having at least one set of a range of attribute values (161, 163, 165).

35. A method as defined in Claim 34, wherein the at least three user-defined operational attributes (161, 163, 165) comprise the supply temperature (161, Ts) of the resource stream (Hl..Hn, CL.Cn), the target temperature (163, Tt) of the resource stream (HL.Hn, CL.Cn), and the heat capacity flow rate (165, FCp) of the resource stream (HL.Hn, CL.Cn).

36. A program product (51) stored in a tangible computer medium to model the energy consumption of a process, said process using a plurality of resource streams (HL.Hn, CL.Cn), with each resource stream having at least one operational attribute (161, 163, 165), said program product (51) including instructions that when executed by a computer cause the computer to perform operations being characterized by:

receiving at least one set of a range of attribute values (161, 163, 165) for at least one attribute (161, 163, 165) of at least one resource stream (HL.Hn, CL.Cn) used by the process;

receiving at least one set of stream-specific minimum approach temperature values ({ΔT_mi_n ¹}) between the streams, at least two of the values (ΔT_mj_n') within each at least one set of stream specific minimum temperature approach values ({ΔT_mm ¹}) being different from each other; and calculating at least one energy consumption value for the process using the received at least one set of a range of attribute values (161, 163, 165) and the at least one set of stream- specific minimum approach temperature values ({ΔT_mj_n'}).

37. A program product (51) as defined in Claim 36, the operations being further characterized by: identifying one or more resource streams (Hl..Hn, Cl..Cn) used within the process to be modeled.

38. A program product (51) as defined in either of Claims 36 or 37, the operations being further characterized by: identifying one or more operational attributes (161, 163, 165) of resource streams (Hl..Hn, Cl..Cn) used within the process which affect the energy consumption of the process.

39. A program product (51) as defined in any of Claims 36-38, the operations being further characterized by: indicating a specific attribute value (161, 163, 165) from the received range supplied which results in the energy consumption value or values calculated.

40. A program product (51) as defined in any of Claims 36-39, wherein the energy consumption value or values calculated are selected from the group consisting of: global minimum energy Qh(minimum) consumed to heat resource streams (Hl..Hn), global minimum energy Qc(minimum) consumed to cool resource streams (C L. Cn), global maximum energy Qh(maximum) consumed to heat resource streams (H 1..Hn), and global maximum energy Qc (maximum) consumed to cool resource streams (Cl..Cn).

41. A program product (51) as defined in any of Claims 36-40, wherein each resource stream (Hl..Hn, Cl..Cn) has at least three user-defined operational attributes (161, 163, 165) each having at least one set of a range of attribute values (161, 163, 165).

42. A program product (51) as defined in Claim 41, wherein the three user-defined operational attributes (161, 163, 165) comprise the supply temperature (161, Ts) of the resource stream (Hl..Hn, Cl..Cn), the target temperature (163, Tt) of the resource stream (Hl..Hn, Cl..Cn), and the heat capacity flow rate (165, FCp) of the resource stream (Hl..Hn, CL.Cn).

43. A computer readable medium that is readable by a computer (31) to optimize energy recovery for a process having a plurality of resource streams (H 1.. Hn, Cl-Cn), the computer readable medium comprising a set of instructions that, when executed by the computer (31), cause the computer (31) to perform operations being characterized by: receiving a plurality of sets of values each defining a potential range of attribute values (161, 163, 165) for a separate one of a corresponding plurality of hot resource streams (Hl..Hn) and aplurality of sets of values each defining a potential range of attribute values (161, 163, 165) for a separate one of a corresponding plurality of cold resource streams (Cl..Cn); assigning at least one set of a plurality of stream specific minimum temperature approach values ({ΔT_min'}) to the plurality of hot resource streams (Hl..Hn), at least two of the values (ΔT_mi_n ¹) within each set of stream specific minimum temperature approach values ({ΔT_mi_n'}) being different from each other; and determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) and a global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]) responsive to the operations of receiving and assigning.

44. A computer readable medium as defined in Claim 43, the operations being further characterized by: determining from the plurality of sets of hot resource streams attribute values (161, 163, 165) and from the plurality of sets of cold resource streams attribute values (161, 163, 165), a first set of discrete process conditions that provide the minimum global minimum heating energy utility (Qh(minimum)) defining a first driving force distribution; and determining from the plurality of sets of hot resource streams attribute values (161, 163, 165) and from the plurality of sets of cold resource streams attribute values (161, 163, 165), a second set of discrete process conditions that provide the minimum global minimum cooling energy utility (Qc(minimum)) defining a second driving force distribution.

45. A computer readable medium as defined in Claim 44, the operations being further characterized by: calculating a first energy cost requirement associated with implementing the first set of process conditions; calculating a second energy cost requirement associated with implementing the second set of process conditions; selecting the set of process conditions from the first and second sets of process conditions that provides a minimum energy cost requirement responsive to determining the first and the second energy cost requirements associated with the first and the second process conditions; designing a heat exchanger network responsive to the selected set of process conditions; and determining a capital cost associated with the heat exchanger network.

46. A computer readable medium that is readable by a computer (31) to optimize energy recovery for a process having a plurality of resource streams, the computer readable medium comprising a set of instructions that, when executed by the computer (31), cause the computer (31) to perform operations being characterized by: assigning at least one set of stream specific minimum temperature approach values ({ΔT_min¹}) to a plurality of hot resource streams (HL.Hn), at least two of the values within each set being different from each other; and determining a global minimum heating energy utility (Qh(minimum)) responsive to the operation of assigning the at least one set of stream specific minimum temperature approach values ((AT_nJ)).

47. A computer readable medium as defined in Claim 46, wherein each hot resource stream has a set of values defining a potential range of attribute values (161, 163, 165) for an attribute (161, 163, 165) of the respective hot resource stream (Hl..Hn), the computer readable medium further comprising a set of instructions that, when executed by the computer (31), cause the computer (31) to perform the operations of: determining a global minimum cooling energy utility (Qc(minimum)) responsive to the operation of assigning the at least one set of stream specific minimum temperature approach values ((AT₁J)); determining from the plurality of sets of hot resource streams attribute values (161, 163, 165) and from a plurality of sets of cold resource streams attribute values (161, 163, 165) a first set of discrete process conditions that provide the minimum global minimum heating energy utility (Qc(minimum)) defining a first driving force distribution; and determining from the plurality of sets of hot resource streams attribute values (161, 163, 165) and from the plurality of sets of cold resource streams attribute values (161, 163, 165), a second set of discrete process conditions that provide the minimum global minimum cooling energy utility (Qc(minimum)) defining a second driving force distribution.

48. A computer readable medium as defined in Claim 47, the operations being further characterized by: calculating a first energy cost requirement associated with implementing the first set of process conditions; calculating a second energy cost requirement associated with implementing the second set of process conditions; selecting the set of process conditions from the first and second sets of process conditions that provides a minimum energy cost requirement responsive to determining the first and the second energy cost requirements associated with the first and the second process conditions; designing a heat exchanger network responsive to the selected set of process conditions; and determining a capital cost associated with the heat exchanger network.

49. A computer readable medium as defined in either of Claims 47 or 48, the operations being further characterized by: receiving a plurality of sets of values

each defining a potential range of stream specific minimum approach values (ΔT_mi_n') for a separate one of the plurality of hot resource streams (Hl..Hn).

50. A computer readable medium as defined in Claim 49, wherein the operations of determining a global minimum heating energy utility Qh(minimum) and a global minimum cooling utility Qc(minimum), determining first and second sets of discrete process conditions, calculating first and second energy cost requirements associated therewith, and selecting the set of process conditions providing the minimum energy cost requirement, are performed incrementally for each of a plurality of combinations of stream specific minimum approach values (ΔT_m,n') falling within the potential range of stream specific minimum approach values (ΔT_mj_n'[L:U]) for each one of the plurality of hot resource streams (Hl..Hn).

51. A computer readable medium as defined in any of Claims 46-50, wherein the operation of determining a global minimum heating energy utility Qh(minimum) includes the operations of: determining a plurality of temperature step intervals (170) for one set of the plurality of stream specific minimum temperature approach values (ΔT_mi_n') assigned to the plurality of hot resource streams (Hl..Hn) responsive to the operation of assigning the at least one set of stream specific minimum temperature approach values ({ΔT_mi_n'}) to the plurality of hot resource streams (Hl..Hn), each temperature step interval (170) having an input interval (171) indicating heat extracted collectively from the plurality of hot resource streams (Hl.. Hn), an output interval (173, Q_s ^l0W-^0UtP^ut , Q_sw^gh_^outpu_{t) indicating heat} collectively applied to a plurality of cold resource streams (CL.Cn), and an output interval (175, Q_S ^1OW-^SUΓP^IUS ₅ (^-^¹^ _{indicating surplus heat} available for a next of the plurality of temperature step intervals (170); determining a minimum high output value (min[Q_s ^hlsh-^0Utput]) and a minimum low output value (min[Q_s ^l0W-^output]) of the plurality of temperature step intervals (170); and determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) responsive to the operation of determining the minimum high output value (min[Q_s ^hish-^0Utput]) and the minimum low output value (min[Q_s ^low-^output]).

52. A computer readable medium as defined in Claim 51, wherein the global minimum heating energy utility interval [Qh(minimum), Qh(maximum)] is determined according to the following calculation:

Qh= [Qh(minimum), Qh(maximum)]; where:

Qh(minimum) =ABS Min{ Min[ Q_s ^his^h-^0UtPut]₅ 0.0} ;

Qh(maximum) =ABS Min{ Min[ Q_s'^0W-^0UtPut], 0.0 }; where:

Q low output _ _/~_\low output , SΛIOW Surplus .

.

where:

_Qs ^!ow_^sur _P ^ius ₌ (£_FCp? -∑FQ>f* )(Th, -Tc₁) ; and

A=I y=i

_Qs ^high__Su_Φ ^lu^s _ -∑FCp'°^W){Th_s -Tc₁) .

53. A computer readable medium as defined in either of Claims 51 or 52, wherein the one set of the plurality of specific minimum temperature approach values ({ΔTmin¹}) includes the lowest range value ({ΔTmin'tL}}) of each value range within the set of the plurality of stream specific minimum approach values ({ΔTmin¹}) assigned to each of the plurality of hot resource streams (Hl.. Hn); wherein the operations are further characterized by: receiving a plurality of sets of values (ΔTmin¹ [L:U]) each separately defining a potential range of stream specific minimum approach values (ΔT_mi_n ¹) for each of the plurality of hot resource streams (Hl..Hn); and wherein the operations of determining a plurality of temperature step intervals (170), determining a minimum high output value (min[Q_s ^1S -^output]) and a minimum low output value (min[Q_s ^low-^output]), and determining a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]) are performed for each of a plurality of sets of specific minimum temperature approach values ({ΔTmin¹}).

54. A computer readable medium as defined in any of Claims 51-53, the operations being further characterized by: determining the global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]), which further includes the operations of: determining a symmetric image ([Qh(maximum), Qh(minimum)]) of a global minimum heating energy utility interval ([Qh(minimum), Qh(maximum)]); determining a plurality of temperature step intervals (170) for one set of the plurality of stream specific minimum temperature approach values ({ΔTm_Jn ¹}) assigned to the plurality of hot resource streams (H 1..Hn) responsive to the operation of determining the symmetric image ([Qh(maximum), Qh(minimum)]), each temperature step interval (170) having an input (171) indicating heat extracted collectively from the plurality of hot resource streams (Hl..Hn), an output (Q_s ^l0W-^0UtP^ut , Q^Ms^h.^outpui) _{indicating heat} collectively applied to the plurality of cold resource streams (Cl..Cn), and an output (Q_s'™-^surPlus , q mum^ indicating surplus heat available for a.next of the plurality of temperature step intervals (170); determining a minimum high output value (min[Q_s ^hl8h-^0Utput]) and a minimum low output value (min[Q_s ^l0W-^0Utput]) of the plurality of temperature step intervals (170); and determining the global cooling energy utility interval ([Qc(minimum), Qc(maximum)]) responsive to the determining of the minimum high output value (min[Q_s ^h'^8h-^0Utput]) and the minimum low output value (min[Q_s ^low-^output]).

55. A computer readable medium as defined in Claim 54, wherein the operations of determining a plurality of temperature step intervals (170), determining a minimum high output value (min[Q_s ^hi8h-^output]) and a minimum low output value (min[Q_s ^l0W-^0Utput]), and determining a global minimum cooling energy utility interval ([Qc(minimum), Qc(maximum)]) responsive thereto, of the operation of determining the global minimum cooling energy utility interval ([Qc(minimurn), Qc(maximum)]), are performed for each of a plurality of sets of specific minimum temperature approach values ((AT_n,;,,¹)).