WO2017060985A1 - エネルギー解析支援装置、エネルギー解析支援方法、エネルギー解析支援プログラム及びこれを記録した記録媒体 - Google Patents

エネルギー解析支援装置、エネルギー解析支援方法、エネルギー解析支援プログラム及びこれを記録した記録媒体 Download PDF

Info

Publication number
WO2017060985A1
WO2017060985A1 PCT/JP2015/078483 JP2015078483W WO2017060985A1 WO 2017060985 A1 WO2017060985 A1 WO 2017060985A1 JP 2015078483 W JP2015078483 W JP 2015078483W WO 2017060985 A1 WO2017060985 A1 WO 2017060985A1
Authority
WO
WIPO (PCT)
Prior art keywords
plant
representing
model
equation
equipment
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2015/078483
Other languages
English (en)
French (fr)
Inventor
智志 桐生
吉雄 丹下
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuji Electric Co Ltd
Original Assignee
Fuji Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fuji Electric Co Ltd filed Critical Fuji Electric Co Ltd
Priority to PCT/JP2015/078483 priority Critical patent/WO2017060985A1/ja
Priority to JP2017544112A priority patent/JP6477906B2/ja
Publication of WO2017060985A1 publication Critical patent/WO2017060985A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B17/00—Systems involving the use of models or simulators of said systems
    • G05B17/02—Systems involving the use of models or simulators of said systems electric

Definitions

  • the present invention relates to a technology for supporting analysis of a plant having resource supply facilities for supplying resources such as electric power and gas and resource demand facilities for consuming resources and consuming energy.
  • Patent Document 2 proposes a method for supporting decision-making at the time of configuration change by calculating an operation method for optimizing energy consumption etc. while changing the device configuration of the system.
  • Non-Patent Document 1 a universal symbol in which some variables are generically referred to as limiting symbols for logical expressions in which equations or inequalities of polynomials are combined by combining symbols such as logical product ( ⁇ ) or logical disjunction ( ⁇ ) Generate a first-order predicate logical expression given ⁇ ) and an existence symbol ( ⁇ ). Then, the system is optimized by eliminating the variable (bound variable) to which the quantifier is added in this first-order predicate logical expression and deriving a logical expression that the other variables (free variables) should satisfy.
  • ⁇ logical product
  • ⁇ logical disjunction
  • Patent Document 3 a control system analysis / design that converts a control system into a first-order predicate logical expression, erases a variable with a limiting symbol in the first-order predicate logical expression, and analyzes the control system from the obtained expression An apparatus is disclosed.
  • the prior art calculates only the optimum energy consumption when using a specific one of a plurality of facilities constituting the plant. I could not do it. That is, in the prior art, when the characteristics of a predetermined facility such as energy efficiency are changed, it has not been possible to directly determine how the optimum energy consumption changes accordingly.
  • the present invention makes it possible to present the user with information such as how to set the characteristics of equipment while making it unnecessary to repeat the calculation many times when optimizing the characteristics of the equipment constituting the plant. It aims to provide technology that can
  • An energy analysis support device is a plant energy analysis support device including a plurality of facilities, and has the following configuration.
  • An output unit that generates and outputs a graph in which one of the two axes represents the characteristic of the facility to be analyzed that is represented by the variable of.
  • An energy analysis support program is an energy analysis support program for causing an information processing apparatus to execute an energy analysis support process of a plant including a plurality of facilities, and has the following configuration.
  • a recording medium is a recording medium recording an energy analysis support program for causing an information processing apparatus to execute an energy analysis support process of a plant including a plurality of facilities, and the following configuration Have. -A formula representing the model of the plant, given the model of the plant, the equipment model representing the capacity of the equipment constituting the plant, and the external condition representing the energy consumption or load of the plant, the equipment model Generating a logical expression including at least one variable representing the characteristic of the equipment to be analyzed from the expression representing the external condition and the expression representing the external condition; and ⁇ said one variable represents based on the generated logical expression Generating and outputting a graph in which the characteristic of the equipment to be analyzed is one of two axes;
  • An energy analysis support method is an energy analysis support method for a plant including a plurality of facilities, and has the following configuration.
  • the present invention when optimizing the characteristics of the equipment constituting the plant, it is possible to present the user with information such as how to set the characteristics of the equipment while making it unnecessary to repeat the calculation many times. It becomes possible.
  • FIG. 1 It is a block diagram of the energy analysis supporting device which concerns on 1st Embodiment. It is a figure which shows the operation
  • FIG. 8 is a diagram exemplifying a group of mathematical expressions generated by the calculation unit from the system model of FIGS. 5 to 7, an equipment model, and an external condition. It is a figure explaining the method to produce
  • FIG. 18 is a diagram illustrating a graph displayed on the display unit in the first modification.
  • FIG. 18 is a diagram for explaining an external condition mathematical expression generated in the modification 2;
  • FIG. 10 is a diagram showing two first-order predicate logical expressions generated in modification 2;
  • FIG. 18 is a diagram showing a logical expression obtained by applying a elimination algorithm to the first-order logical expression of FIG.
  • FIG. 17 is a diagram illustrating a graph displayed on the display unit in the second modification; It is a figure explaining the equipment model formula generated in modification 3.
  • FIG. 18 is a diagram illustrating a graph displayed on the display unit in the first modification.
  • FIG. 18 is a diagram showing a first-order predicate logical expression generated in modification 3; It is a figure which shows the logical expression obtained by applying a limited-symbol elimination algorithm to the first-order logical expression of FIG.
  • FIG. 21 is a diagram illustrating a graph displayed on the display unit in the modification 3; It is a figure which shows the operation
  • FIG. 30 is a diagram illustrating a group of mathematical expressions generated by the calculation unit 11 from the system model of FIG. 5, the equipment model of FIG. 6, and the external conditions of FIG. It is a figure explaining the method to produce
  • FIG. 30 is a diagram showing a logical expression obtained by applying a elimination algorithm to the first-order logical expression of FIG. It is a figure which shows an example of the screen displayed on a display part in 2nd Embodiment. It is a figure which shows the other example of the screen displayed on a display part in 2nd Embodiment. It is a block diagram of an information processor. It is a figure which shows a recording medium.
  • FIG. 1 is a block diagram of an energy analysis support apparatus according to the present embodiment.
  • the energy analysis support device 1 of FIG. 1 has a calculation unit 11 and an output unit 12, and is connected to an input unit 2 such as a keyboard 2A and a pointing device 2B and a display unit 5 such as a monitor.
  • the energy analysis supporting device 1 receives an input of a system model of a plant, an equipment model, and an external condition via the input unit 2 (2A or 2B).
  • the calculation unit 11 of the energy analysis support device 1 generates a system model formula, an equipment model formula, and an external condition formula from the input information, and generates a first-order predicate formula using the generated formula.
  • the calculation unit 11 applies a definite mark elimination algorithm to the first-order predicate logical expression to obtain a logical expression representing the relationship between the desired variable and the total energy consumption.
  • the output unit 12 generates a graph based on the logical expression obtained by the calculation unit 11 by calculation, and outputs the graph to the display unit 5.
  • the display unit 5 displays a graph based on the information input from the output unit 5.
  • the energy analysis support device 1 may also start processing using the input information, triggered by the input of information such as a system model via the input unit 2. Alternatively, the process may be started according to the start instruction from the user using the information already stored in the device.
  • FIG. 2 is a diagram showing an operation flow in the calculation unit 11.
  • step S1 the calculation unit 11 acquires a system model equation, an equipment model equation, and an external condition equation.
  • the mathematical expression group may be generated by the calculation unit 11 using the system model, the equipment model and the external condition, or may be directly input through the input unit 2.
  • step S2 a first-order predicate logical expression relating to the total energy consumption is generated from the group of mathematical expressions acquired in step S1.
  • step S3 the first-order logical formula generated in step S2 is processed to obtain a total energy consumption formula.
  • step S4 information such as a total energy consumption formula is output to the output unit 12, and the process is ended.
  • FIG. 3 is a diagram for explaining a method of generating a first-order predicate logical expression in step S2 of FIG. 2 and a method of calculating a total energy consumption formula in step S3.
  • the calculation unit 11 combines the system model equation F1, the equipment model equation F2 and the external condition equation F3 with a logical product ( ⁇ ) to generate the equation (3-2) in FIG. Then, for the equation (3-2), the existence symbol ( ⁇ ) is further added to variables excluding the variable a representing the desired property of the facility and the variable P representing the total energy consumption of the plant (FIG. 3 (3-1)). Thus, a first-order predicate formula ⁇ is generated.
  • ⁇ : ⁇ P1, ..., ⁇ PN ⁇ L ⁇ L 1, ..., ⁇ LN (F1 ⁇ F2 ⁇ F3)
  • L is a variable that represents the total load of the plant.
  • N represents the number of facilities
  • Ln represents the load of the facility n.
  • the calculation unit 11 applies the restrictive symbol elimination algorithm to the generated first-order predicate logical expression ⁇ to obtain a logical expression ⁇ representing the relationship satisfied by the variable a representing the property of the facility and the variable P representing the total energy consumption. And outputs this to the output unit 12.
  • the output unit 12 generates a graph representing the relationship between the variable a representing the property of the facility and the variable P representing the total energy consumption of the plant from the logical expression ⁇ input from the calculation unit 11.
  • the output unit 12 outputs, to the display unit 5, information such as the generated graph, that is, information obtained as a result of performing the analysis support process in the energy analysis support apparatus 1.
  • the display unit 5 displays the calculation result in the energy analysis support device 1 based on the information input from the output unit 12.
  • FIG. 4 is a view schematically showing a graph generated by the output unit 12.
  • one axis (horizontal axis) takes a variable a that represents the characteristics of the facility, and the other axis (vertical axis) takes a variable P that represents the total energy consumption of the plant.
  • An area A satisfied by the logical expression ⁇ is visualized.
  • the energy analysis support device 1 performs calculation of various mathematical expressions using the input information such as various models, and obtains a graph of a calculation result using a specific example. explain. First, information such as various models input to the energy analysis supporting apparatus 1 will be described with reference to FIGS. 5 to 7.
  • FIG. 5 is a diagram illustrating a system model of a plant.
  • the air conditioning system of FIG. 5 has a plurality of refrigerators, and the plurality of refrigerators use the electric power supplied from the power supply to supply the thermal load required by the space to be air-conditioned.
  • the system model is composed of the equipment configuration of the air conditioning system shown in FIG. 5 and the heat load and power consumption of each equipment.
  • the refrigerator # 1 consumes the power P1 [kW] with respect to the thermal load L1 [kW].
  • the refrigerator # 2 consumes the power P2 [kW] for the heat load L2 [kW], and the refrigerator # 3 consumes the power P3 [kW] for the heat load L3 [kW].
  • the refrigerator # 4 consumes the power P4 [kW] with respect to the heat load L4 [kW].
  • the sum of the power consumption P1 to P4 consumed by each of the refrigerator # 1 to # 4 is the total power consumption P [kW] consumed in the entire system, and the total heat load L [kW] required by the air conditioning target space
  • the thermal loads L1 to L4 [kW] of the respective refrigerator # 1 to # 4 are distributed.
  • FIG. 6 is a diagram illustrating an installation model.
  • the facility model of FIG. 6 is configured from the power consumption characteristics with respect to the heat load of the existing machines, ie, the refrigerators # 1 to # 3.
  • Each refrigerator # 1 to # 3 has different power consumption characteristics and upper and lower limit constraints.
  • FIG. 7 is a diagram illustrating an external condition.
  • the total heat load L of the entire air conditioning system for responding to the request of the air conditioning target space is set to 500 [kW].
  • FIG. 8 is a diagram illustrating a group of mathematical expressions generated by the calculation unit 11 from the system model, the equipment model, and the external condition of FIGS. 5 to 7.
  • the system model equation F1 is an equation representing the air conditioning system model of FIG. 5, and is an equation obtained by combining two equations (8-1) and (8-2) with a logical product.
  • Formula (8-1) of FIG. 8 shows that the sum total of the power consumption P1 to P4 of each refrigerator # 1 to # 4 is equal to the total power consumption P which the power supply supplies to each refrigerator # 1 to # 4.
  • the equation (8-2) represents that the total heat load L of the space to be air-conditioned is distributed to the refrigerators # 1 to # 4.
  • the equipment model formula F2 is a formula that expresses the equipment model of FIG. 6, and is composed of a formula that combines the formulas representing the characteristics of each of the four refrigerators # 1 to # 4 with a logical product.
  • the 4th line is a numerical expression showing the characteristic of refrigerator # 4 newly introduced among the following expression F2.
  • the variable a included in the equation representing the characteristic of the refrigerator # 4 represents the power consumption characteristic with respect to the load of the refrigerator # 4.
  • the external condition formula F3 is a formula representing the external condition of FIG.
  • the calculation unit 11 holds the system model formula F1, the equipment model formula F2 and the external condition formula F3 in storage means such as a memory (not shown) in FIG. 1, and uses these formulas to hold the following formula Perform the process.
  • FIG. 9 is a diagram for explaining a method of generating a first-order predicate logical expression ⁇ from the mathematical expression groups F1 to F3 of FIG.
  • the calculation unit 11 of the energy analysis support device 1 combines the equation group of FIG. 8, that is, the system model equation F1, the equipment model equation F2, and the external condition equation F3 with a logical product to represent the entire system (F1). Generate ⁇ F2 ⁇ F3) (equation (9-2) in FIG. 9).
  • the calculation unit 11 presents the existence symbol ( ⁇ ) as a variable excluding the variable P representing the total energy consumption and the variable a representing the characteristic of the refrigerator # 4 with respect to the formula F1 ⁇ F2 ⁇ F3 representing the entire system. Is generated ((9-1) in FIG. 9) to generate a first-order predicate formula ⁇ regarding the total power consumption.
  • the calculation unit 11 applies the quantifier elimination algorithm (QE []) to the first-order predicate logical formula ⁇ in FIG.
  • the logical expression ⁇ is obtained by eliminating the variables.
  • the logical expression ⁇ is a relational expression that the variable a representing the characteristic of the refrigerator # 4 and the variable P representing the total power consumption of the air conditioning system satisfy.
  • the logical expression ⁇ represents the range of the total power consumption P that can be performed by changing the load distribution of each of the refrigerator # 1 to # 4 with respect to the power consumption characteristic a of the refrigerator # 4.
  • FIG. 11 is a diagram illustrating a graph displayed on the display unit 5.
  • the calculation unit 11 of the energy analysis support device 1 eliminates the limitation symbol from the first-order predicate logical expression ⁇ to obtain the logical expression ⁇ in FIG.
  • the obtained logical expression ⁇ is passed to the output unit 12 of the energy analysis support device 1.
  • the output unit 12 takes the variable P representing the total power consumption on the vertical axis and the variable a representing the power consumption characteristic of the refrigerator # 4 on the horizontal axis from the logical formula ⁇ , and represents the area A1 satisfied by the logical formula ⁇ Generate a graph. In FIG. 11, the area A1 is indicated by hatching.
  • the output unit 12 outputs the generated image information of the graph to the display unit 5.
  • the display unit 5 displays the graph of FIG. 11 based on the information input from the output unit 12.
  • the lower limit value in the vertical axis direction of the region A1 is the optimal power consumption with respect to the value of the power consumption characteristic a of the refrigerator # 4.
  • the lower limit LM_lw corresponds to the optimum total power consumption with respect to the value of the power consumption characteristic a of the refrigerator # 4.
  • the user refers to the lower limit LM_lw of the graph in FIG. 11 and grasps the total power consumption P of the whole air conditioning system in the case of the operation for realizing the most power saving, and the power consumption characteristic a of the refrigerator # 4 to be introduced Can be determined. Therefore, the user can appropriately determine the characteristic (power consumption characteristic in the example) of the newly introduced facility (refrigerator # 4) with reference to the graph, whereby the newly introduced facility can be appropriately determined. It is possible to make a decision.
  • the first-order predicate logical expression is generated from the input system model, the facility model, and the external condition. Then, a punctuation elimination algorithm is applied to the generated first-order predicate logical expression to obtain a logical expression representing the relationship between the characteristics of a certain facility of the plant and the total power consumption of the plant.
  • the first-order predicate formula is a formula that expresses the whole of the plant generated from the input system model etc., variables other than the variables that you want to finally obtain as the processing result, ie, the characteristics of a certain facility and total power consumption. It is generated by giving the presence symbol ( ⁇ ) to a variable.
  • the area where the characteristics of the equipment and the total energy consumption satisfy includes both the area when the equipment to be introduced is operated and the case where the equipment is not operated. .
  • this modification it differs in the point which draws the graph which distinguished the case of the case of operating a certain installation, and the case of not making it operate.
  • FIGS. 1 and 5 to 7 As described with reference to.
  • the energy analysis supporting device 1 When the information shown in FIG. 5 to FIG. 7 is input, the energy analysis supporting device 1 first generates a system model equation, an equipment model equation, and an external condition equation, as in the above embodiment.
  • the system model equation F1 and the external condition equation F3 are the same as in the previous embodiment, and are as described with reference to FIG.
  • FIG. 12 is a diagram illustrating an equipment model equation generated in the present modification.
  • the calculation unit 11 of the energy analysis support apparatus 1 uses, as an equipment model equation, equations F2 ′ and F2 when the refrigerator # 4 to be newly introduced to the air conditioning system is not operated and when it is operated.
  • the equipment model formula F2 ′ ′ is expressed by the formula of the equipment model formula F2 in FIG. 8 excluding the following formula (12-1) representing the condition that the refrigerator # 4 is not operated.
  • the calculation unit 11 generates two first-order predicate logical expressions using each of the generated two equipment model expressions F2 ′ and F2 ′ ′.
  • FIG. 13 is a diagram showing first-order predicate logical expressions' 2 'and' 2 'generated using the expressions F1 and F3 of FIG. 8 and the expressions F2' and F2 '' of FIG.
  • FIG. 14 is a diagram showing logical expressions ⁇ 2 ′ and ⁇ 2 ′ ′ obtained by applying the restrictive symbol elimination algorithm to the first-order predicate logical expressions ⁇ 2 ′ and ⁇ 2 ′ ′.
  • the logical formulas ⁇ 2 ′ and ⁇ 2 ′ ′ represent the relationship between the power consumption characteristic a of the refrigerator # 4 and the total power consumption P when the refrigerator # 4 newly introduced to the air conditioning system is not operated and operated, respectively.
  • the output unit 12 receives the equations (14-1) and (14-2) from the calculation unit 11, the output unit 12 generates a graph capable of comparing the regions satisfied by the two equations based thereon. Then, the output unit 12 outputs this to the display unit 5 as information obtained as a result of performing the analysis support process in the energy analysis support device 1.
  • the display unit 5 displays a graph based on the input information.
  • FIG. 15 is a diagram illustrating a graph displayed on the display unit 5. As in the graph of FIG. 11, the vertical axis represents a variable P representing the total power consumption, and the horizontal axis represents a variable a representing the power consumption characteristic of the refrigerator # 4.
  • the user performs analysis by comparing the respective regions when the newly introduced refrigerator # 4 is not in operation and in operation, with reference to the graph illustrated in FIG. 15.
  • the refrigerator # 4 when the value of the power consumption characteristic a of the refrigerator # 4 is small, the refrigerator # 4 is operated and the other inefficient refrigerators are stopped to stop the operation, that is, the refrigerator # 1. It can be seen that the total power consumption P can be significantly suppressed as compared with the case of operating # 3 (* 1 in FIG. 15).
  • the total power consumption P can be reduced by not operating the refrigerator # 4, that is, even if a new refrigerator # 4 is introduced, the total consumption It can be seen that the power P can not be suppressed (* 2 in FIG. 15).
  • Modification 2 In said embodiment and modification 1, the case where the air conditioning object space is relatively high heat load is set as an external condition. On the other hand, in this modification, the case where the air conditioning target space has a relatively low thermal load is set, and a logical expression ⁇ representing the relationship between the characteristics of the facility to be introduced at the low thermal load and the total energy consumption is determined. Draw on a graph.
  • the energy analysis supporting device 1 When the information shown in FIGS. 5 to 7 is input, the energy analysis supporting device 1 generates a system model equation, an equipment model equation, and an external condition equation, as in the above embodiment and the first modification.
  • facility model formulas F2 ′ and F2 ′ ′ are generated for each of the case where the newly introduced refrigerator # 4 is operated and the case where it is not operated, and the respective refrigerator # 4 is generated.
  • the case of obtaining the relational expression between the power consumption characteristic a and the total power consumption P will be described.
  • the system model formula F1 is the same as that of the above embodiment and is as described with reference to FIG.
  • the equipment model formulas F2 ′ and F2 ′ ′ are as described in the first modification with reference to FIG.
  • FIG. 16 is a diagram for explaining an external conditional expression generated in the present modification.
  • the thermal load L in the graph of the total thermal load of the air-conditioned space according to the time zone on the upper side of FIG. 16 is a time zone such as "6 o'clock to 7 o'clock" Set 400 [kW] which is the value of relatively low total heat load set in.
  • the calculation unit 11 generates two first-order predicate logical expressions using the system model expression F1, the two equipment model expressions F2 'and F2' ', and the external condition expression F3' at low heat load.
  • FIG. 17 is a diagram showing two generated first order predicate formulas '3' and '3'.
  • Two first-order predicate logic formulas ⁇ 3 'and ⁇ 3''in FIG. 17 are first-order predicate formulas relating to the total power consumption when the refrigerator # 4 is not operated and when it is operated.
  • the calculation unit 11 is an equation (F1 ⁇ F2′ ⁇ F3 ′) representing the entire system from the system model equation F1 of FIG. 8, the equipment model equation F2 ′ and F2 ′ ′ of FIG. 12, and the external condition equation F3 ′ of FIG. ), (F1 ⁇ F2′′ ⁇ F3 ′).
  • the variable to which the presence symbol ( ⁇ ) is assigned is the same as in the above-described embodiment and the first modification.
  • FIG. 18 is a diagram showing logical expressions ⁇ 3 ′ and ⁇ 3 ′ ′ obtained by applying the restrictive symbol elimination algorithm to the first-order predicate logical expressions ⁇ 3 ′ and ⁇ 3 ′ ′.
  • the logical expressions ⁇ 3 ′ and ⁇ 3 ′ ′ are relational expressions that the power consumption characteristic a of the refrigerator # 4 and the total power consumption P satisfy when the refrigerator # 4 is not operated and when the refrigerator # 4 is operated at low heat load, respectively.
  • the output unit 12 receives the logical expressions ⁇ 3 ′ and ⁇ 3 ′ ′ from the calculation unit 11, the output unit 12 generates a graph capable of comparing the regions satisfied by the two expressions based thereon, and outputs the graph to the display unit 5.
  • FIG. 19 is a diagram illustrating a graph displayed on the display unit 5. As in the graphs of FIGS. 11 and 15, the vertical axis represents the total power consumption P, and the horizontal axis represents the power consumption characteristic a of the refrigerator # 4.
  • a graph is created in which areas A3 'filled with the two logical expressions ⁇ 3' and ⁇ 3 '' are hatched to make it possible to distinguish A3 '' by oblique lines.
  • the user performs analysis by comparing the respective regions when the refrigerator # 4 to be newly introduced at low heat load is not operated and when it is operated.
  • the refrigerator # 4 is operated when the range of the power consumption characteristic a of the refrigerator # 4 is from 0.10 to 035 (the whole area of the graph shown in FIG. 15). It is understood that the total power consumption P can be suppressed if the user is allowed to do so.
  • FIG. 19 the relationship between the power consumption characteristic a of the refrigerator # 4 and the total power consumption P at the time of low heat load is illustrated as a graph which can be identified when the refrigerator # 4 is not operated and when it is operated. doing.
  • various modified examples can be considered such as a graph comparing high heat load time and low heat load time by using logical expressions ⁇ 2 ′ and ⁇ 2 ′ ′ of FIG. 14 obtained in the first modification, for example.
  • Modification 3 In the above embodiment and Modifications 1 and 2, the relationship between the variable a representing the power consumption characteristic which is one of the characteristics of the equipment (refrigerator) and the total power consumption P is obtained. On the other hand, in this modification, the relationship between the variable representing the characteristic other than the power consumption characteristic of the newly introduced refrigerator and the total power consumption P is determined.
  • the configuration of the energy analysis support device 1 and the information to be input are the same as those of the above-described embodiment and the like, and are as described with reference to FIG. 1 and FIG. 5 to FIG.
  • the energy analysis supporting device 1 When the information shown in FIGS. 5 to 7 is input, the energy analysis supporting device 1 generates a system model equation, an equipment model equation, and an external condition equation, as in the above embodiment and the like.
  • the system model equation F1 and the external condition equation F3 are as described above with reference to FIG.
  • FIG. 20 is a diagram for explaining an equipment model formula generated in the present modification.
  • a facility model formula F2 ′ ′ ′ in the case of operating the refrigerator # 4 is generated using a variable b that represents the upper limit of the heat load of the refrigerator # 4.
  • the power consumption characteristic of the refrigerator # 4 is fixed to 0.16.
  • a formula representing the characteristic of the refrigerator # 4 is as follows.
  • the calculation unit 11 generates a first-order predicate logical expression using the system model expression F1, the facility model expression F2 ′ ′ ′, and the external condition expression F3.
  • FIG. 21 is a diagram showing the generated first-order predicate logical expression '4'.
  • the first-order predicate logical expression ⁇ 4 ′ ′ in FIG. 21 is a first-order predicate logical expression in which the upper limit b of the heat load of the refrigerator # 4 is a variable.
  • the calculation unit 11 generates an equation (F1 ⁇ F2 ′ ′ ′ ⁇ F3) representing the entire system from the system model equation F1 of FIG. 8 and the facility model F2 ′ ′ ′ of FIG. 20 and the external condition equation of FIG.
  • the existence symbol ( ⁇ ) is added to variables other than the variables b and P to generate a logical expression ⁇ 4 ′ ′.
  • FIG. 22 is a diagram showing a logical expression ⁇ 4 ′ ′ obtained by applying the restrictive symbol elimination algorithm to the first-order predicate logical expression ⁇ 4 ′ ′.
  • the logical formula ⁇ 4 ′ ′ represents the relationship between the upper limit b of the heat load L4 of the refrigerator # 4 and the total power consumption P.
  • the output unit 12 receives the logical expression ⁇ 4 ′ ′ from the calculation unit 11, the output unit 12 generates a graph representing an area that the logical expression ⁇ 4 ′ ′ satisfies based on the logical expression ⁇ 4 ′ ′ and outputs the graph to the display unit 5.
  • FIG. 23 is a diagram illustrating a graph displayed on the display unit 5.
  • the vertical axis represents the total power consumption P, as in the graphs of FIGS.
  • the upper limit b of the heat load L4 of the refrigerator # 4 is taken.
  • the area A4 '' represents the area that the variables b and P of the logical expression [phi] 4 '' satisfy, by hatching.
  • the region A4 ′ ′ is the upper limit b of the thermal load L4 of the refrigerator # 4 and the total power consumption P when operating the refrigerator # 4. Is an area that the
  • an area A4 'in FIG. 23 is an area where the upper limit b of the heat load L4 of the refrigerator # 4 and the total power consumption P when the refrigerator # 4 is not operated.
  • Region A 4 ′ utilizes the logical expression ⁇ 2 ′ determined for the case where the refrigerator # 4 is not operated in the second modification.
  • the logical expression ⁇ 2 ′ is obtained by eliminating the limitation from the first-order predicate logical expression ⁇ 2 ′, and the first-order predicate logical expression ⁇ 2 ′ uses the system model expression F1, the facility model expression F2 ′ and the external condition expression F3. Is generated.
  • the equipment model formula F2 'does not depend on the upper limit b of the heat load L4 of the refrigerator # 4. Therefore, the logical expression ⁇ 2 ′ determined using the equipment model F2 ′ is also not based on the variable b. Therefore, the logical expression ⁇ 2 ′ is the upper limit b of the heat load L4 of the refrigerator # 4 of the present modification and the total b It can be said that the relationship with the power consumption P is expressed.
  • the user refers to the graph illustrated in FIG. 23 to compare and analyze the respective regions when the refrigerator # 4 is operated and not operated.
  • the energy analysis support device 1 eliminates the limitation sign from the generated first-order predicate logical expression, obtains a logical expression representing the relationship between the property of the facility and the total consumption energy, and obtains the graph from this. It is generated.
  • the logical expression obtained as a result of eliminating the limitation symbol differs in that it represents a relationship which is satisfied between the characteristics of the equipment.
  • the energy analysis supporting device 1 will be described focusing on differences from the first embodiment.
  • the configuration diagram of the energy analysis support device 1 according to the present embodiment is the same as that according to the first embodiment of FIG. 1.
  • the energy analysis support device 1 includes a calculation unit 11 and an output unit 12, and is connected to the input unit 2 such as the keyboard 2A and the pointing device 2B and the display unit 5 such as a monitor.
  • the energy analysis support device 1 receives inputs of a system model, an equipment model, and an external condition via the input unit 2 as in the first embodiment.
  • the information to be input is also the same as that of the above embodiment, and is as described with reference to FIGS. 5 to 7 and the like.
  • the energy analysis support device 1 generates a mathematical expression such as a system model mathematical expression from the input information to generate a first-order predicate logical expression, and erases a limiting symbol from the generated first-order predicate logical expression. And get the logical formula.
  • the logical expressions obtained by eliminating the limitation symbols in the present embodiment represent the relationships that the variables representing the characteristics of the plant equipment satisfy.
  • the energy analysis support device 1 generates a graph based on the obtained logical expression, outputs the graph to the display unit 5, and displays the graph.
  • FIG. 24 is a diagram showing an operation flow in the calculation unit 11.
  • Step S11 is the same as step S1 of FIG. 2 in comparison with the operation flow of the first embodiment of FIG.
  • step S12 the calculation unit 11 generates a first-order predicate logical expression relating to the operating state of the plant.
  • “the first-order predicate logical expression relating to the operating state of the plant” is different from the “first-order predicate logical expression relating to the total consumption energy” generated in step S2 of FIG. .
  • steps S13 and S14 processing similar to that of steps S3 and S4 in FIG. 2 is performed. That is, a logical expression in which the limiting symbol has been eliminated is calculated from the first-order predicate logical expression generated in step S12, the plant analysis mathematical expression obtained by the calculation is output to the output unit 12, and the processing is ended.
  • FIG. 25 is a diagram for explaining a first-order predicate logical expression generated by the calculation unit 11 in the present embodiment and a logical expression obtained by deleting a limitation symbol from the first-order predicate logical expression.
  • Calculation unit 11 like the first embodiment described above, the system model equation F1, the equipment model F2 and external conditions formula F3 i binds with logical ( ⁇ ) (equation 25 (25-2)) . Then, in the present embodiment, the calculation unit 11 assigns an existing symbol ( ⁇ ) to variables other than the two variables a and b representing the characteristic of the facility ((25-1) in FIG. 25). Generate a rank predicate formula ⁇ i .
  • the first-order predicate logical expression ⁇ ⁇ i generated in the present embodiment is referred to as “first-order predicate logic relating to the operating state of the plant” below. It is also called a formula.
  • the calculation unit 11 applies the limited symbol elimination algorithm to each of the first-order predicate logical expressions ⁇ i relating to the operating states of I plants, and represents a logical expression i i representing the relationship satisfied by the variables a and b representing the characteristics of the equipment. And output this to the output unit 12.
  • Generate The output unit 12 outputs this to the display unit 5 as information obtained as a result of performing the analysis support processing in the energy analysis support device 1.
  • the display unit 5 displays the calculation result in the energy analysis support device 1 based on the information input from the output unit 12.
  • FIG. 26 is a view schematically showing a graph generated by the output unit 12.
  • one axis (horizontal axis) takes a variable a representing the property of the facility
  • the other axis (vertical axis) takes another variable b representing the property of the facility.
  • the region satisfied by the logical expression ⁇ i obtained for each of the operating conditions of the plant, that is, each of the plurality of external condition expressions F3 i is visualized by color coding or the like.
  • the system model and the facility model of the plant are the same as in the first embodiment, as shown in FIGS. 5 and 6.
  • FIG. 27 is a diagram illustrating an external condition.
  • the first external condition is the total heat load L of the entire air conditioning system.
  • the total heat load L is set to 500 [kW].
  • the second external condition is the total power consumption P of the entire air conditioning system.
  • four conditions are set as the second external condition, and 150 kW or less, 140 kW or less, 130 kW or less, and 120 kW or less are set.
  • FIG. 28 is a diagram illustrating a group of mathematical expressions generated by the calculation unit 11 from the system model of FIG. 5, the equipment model of FIG. 6, and the external conditions of FIG.
  • the power consumption characteristic of the equipment model formula F2 is represented by “variable a” as in the first embodiment.
  • the upper limit of the heat load L4 of the refrigerator # 4 is represented by "variable b”.
  • FIG. 29 is a diagram for explaining a method of generating a first-order logic formulas [psi i from equation group F1, F2, F3 i in FIG 28.
  • calculation unit 11 the generated equation F1, F2, F3 i of FIG. 28 coupled with logical AND ( ⁇ ), the variable a, existential quantifiers on variables other than b ( ⁇ ) to generate a first-order predicate formula ⁇ i .
  • the external condition of i 1, specifically, the total heat load L of the air conditioning system is 500 kW, and the total power consumption P of the air conditioning system is 150 kW First-order predicate formula ⁇ 1 in the case where the following condition is an external condition is exemplified.
  • Calculation unit 11 outputs the generated first order predicate each quantifier elimination algorithm formulas ⁇ 1 ⁇ ⁇ 4 was (QE []) by applying the logical expression was erased granted variables existential quantifier phi 1 ⁇ phi Get four .
  • the logical expressions ⁇ 1 to ⁇ 4 are relational expressions that the power consumption characteristic a of the refrigerator # 4 and the upper limit b of the heat load L4 satisfy.
  • FIG. 31 is a view showing an example of a screen displayed on the display unit 5, and the regions A1 to A4 satisfied by the logical expressions ⁇ 1 to ⁇ 4 are shown by graphs.
  • the ordinate represents the upper limit b of the heat load L4 of the refrigerator # 4, and the abscissa represents the power consumption characteristic a of the refrigerator # 4.
  • Areas A1 to A4 in FIG. 31 represent the installed capacity of the refrigerator # 4 required to reduce the total power consumption P of the air conditioning system to 150, 140, 130, and 120 [kW] or less.
  • the ranges of the regions A1 to A4 satisfied by the power consumption characteristic a of the refrigerator # 4 and the upper limit b of the heat load L4 differ depending on the external conditions.
  • the user changes the total power consumption P according to the correlated characteristics such as the power consumption characteristic a of the refrigerator # 4 and the upper limit b of the heat load L4. Can be considered.
  • the regions A1 to A4 may be arranged side by side on the screen as illustrated in FIG. 31, or one graph may display visibly a portion of each region overlapping.
  • FIG. 32 is a view showing another example of the screen displayed on the display unit 5, and the areas A1 to A4 satisfied by the logical expressions ⁇ 1 to ⁇ 4 are displayed distinguishably by color coding or the like.
  • the areas A1 to A4 are displayed so as to be distinguishable by the user by hatching or hatching.
  • the user attempts to suppress the total power consumption P to a predetermined value or less, for example, when the power consumption characteristic a of the newly introduced refrigerator # 4 is 0.11.
  • the optimum value of the power consumption characteristic a of the refrigerator # 4 is determined when trying to suppress the total power consumption to a predetermined value or less. can do. It can be seen that the power consumption characteristic a should be 0.27 or less in order to reduce the total power consumption to 150 kW or less under the condition that the upper limit of the heat load is 200 kW or more. Similarly, it is understood that the power consumption characteristic a should be set to 0.18 or less in order to set the total power consumption to 140 [kW] or less under the same conditions.
  • a formula representing the relationship between the characteristics of a plurality of correlated facilities is obtained, and how the total power consumption of the plant changes. Can be visualized. For the user, using the visualized result, for example, it becomes easy to study how much the energy saving effect can be expected if equipment with which efficiency and constraints are introduced.
  • the constraint condition (upper limit of heat load in the above example) of the equipment to be newly introduced is set to one variable b, and a logical expression ⁇ i including this is obtained and provided to the user for analysis. doing.
  • energy saving may be realized when, for example, the constraint condition is alleviated (in the above example, the upper limit of the heat load is increased), such a state Can also be verified.
  • the energy analysis support device 1 of FIG. 1 can be configured, for example, using an information processing device (computer) as shown in FIG.
  • the information processing apparatus shown in FIG. 33 includes a CPU (central processing unit) 1001, a memory 1002, an input unit 1003, an output unit 1004, an external storage unit 1005, a medium drive unit 1006, and a network connection unit 1007. It is connected.
  • the memory 1002 includes, for example, a read only memory (ROM), a random access memory (RAM), and the like, and stores programs and data used for processing.
  • the CPU 1001 performs necessary processing by executing a program using the memory 1002.
  • the calculation unit 11 and the output unit 12 in FIG. 1 correspond to functions implemented by executing a program stored in the memory 1002.
  • the input device 1003 is, for example, a keyboard, a pointing device, a touch panel, etc., and is used to input an instruction or information from an operator such as a network administrator.
  • the output device 1004 is, for example, a display, a printer, a speaker, and the like, and is used to output an inquiry to a user, an alarm, a processing result, and the like.
  • the external storage device 1005 is, for example, a magnetic disk device, an optical disk device, a magneto-optical disk device, a tape device or the like.
  • the information processing apparatus stores the program and data in the external storage device 1005, loads them into the memory 1002 and uses them as needed.
  • the medium drive device 1006 drives a portable recording medium 1009 and accesses the recorded contents.
  • the portable recording medium 1009 is an arbitrary computer-readable recording medium such as a memory card, a flexible disk, a compact disk read only memory (CD-ROM), an optical disk, a magneto-optical disk and the like.
  • the operator stores the program and data in the portable storage medium 1009 and loads them into the memory 1002 for use as needed.
  • the network connection device 1007 is connected to an arbitrary communication network such as a local area network (LAN) or the Internet, and performs data conversion involved in communication.
  • the information processing apparatus receives the program and data from an external apparatus via the network connection apparatus 1007, loads them into the memory 1002, and uses them as necessary.
  • FIG. 34 shows a computer readable recording medium capable of supplying a program and data to the information processing apparatus of FIG.
  • Programs and data stored in the portable storage medium 1009 or the database 1103 of the server 1101 are loaded into the memory 1002 of the information processing apparatus 1102.
  • the server 1101 generates a carrier signal that carries the program and data, and transmits the carrier signal to the information processing apparatus 1102 via any transmission medium on the network.
  • the CPU 1001 executes the program using the data and performs necessary processing.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

プラントを構成する設備の特性を最適化する際に、何度も計算を繰り返すことを不要としつつ、どのように設備の特性を設定すべきか等の情報をユーザに提示することのできる技術を提供する。複数の設備からなるプラントのエネルギー解析支援装置１において、計算部１１は、プラントのモデルと、プラントを構成する設備の能力を表す設備モデルと、プラントの消費エネルギーまたは負荷を表す外部条件とが与えられると、プラントのモデルを表す数式、設備モデルを表す数式及び外部条件を表す数式より、解析対象の設備の特性を表した少なくとも一の変数を含む論理式を生成する。出力部１２は計算部１１により生成された論理式に基づいて、一の変数が表す解析対象の設備の特性を２軸のうちの一の軸とするグラフを生成して出力する。

Description

エネルギー解析支援装置、エネルギー解析支援方法、エネルギー解析支援プログラム及びこれを記録した記録媒体
　本発明は、電力、ガス等の資源を供給する資源供給設備と資源を需要してエネルギーを消費する資源需要設備とを有するプラントの解析を支援する技術に関する。
　複数の設備を備え、設備の運転時にはエネルギーを消費するプラントの運転を制御する際には、消費電力等の消費エネルギーを最小化するように作成された運転スケジュールに基づいて、設備ごとの負荷配分が決定されている。これに関して、例えば、特許文献１では、外気温度や日射量、在室人数予測等の外部条件から建物の熱負荷を予測して、プラント全体における総消費電力が最小となる運転スケジュールを立案するシステム制御装置が開示されている。
　また、特許文献２では、システムの装置構成を変更しながら消費エネルギー等を最適にする運用方法を計算することで、構成変更時の意思決定を支援する方法が提案されている。
　一方、システム制御等の問題を一階述語論理式で表現し、これを解くことで、システムの最適化を行う技術が知られている（例えば、非特許文献１を参照）。具体的には、多項式の等式や不等式を論理積（∧）や論理和（∨）等の結合記号で結合した論理式に対して、一部の変数に限定記号と総称される全称記号（∀）や存在記号（∃）を付与してなる一階述語論理式を生成する。そして、この一階述語論理式における限定記号の付与された変数（束縛変数）を消去して、それ以外の変数（自由変数）が満たすべき論理式を導くことで、システムの最適化を行う。
　例えば、特許文献３では、制御系を一階述語論理式に変換し、一階述語論理式における限定記号の付いた変数を消去し、得られた式から制御系を解析する制御系解析・設計装置について開示されている。
特開２０１１－２１４７９４号公報 特開２０１０－２５０４７３号公報 特開平１１－３２８２３９号公報
穴井宏和・横山和弘著、「ＱＥの計算アルゴリズムとその応用　数式処理による最適化」、東京大学出版会、２０１１年、ｐ．２１４－２２１
　上記の一階述語論理式を解くことによりシステムを最適化する技術を含め、従来の技術では、プラントを構成する複数の設備のうちの特定の設備を使用した場合における最適な消費エネルギーについてしか計算することができなかった。すなわち、従来技術では、エネルギー効率等の所定の設備の特性を変化させた場合に、これに応じて最適な消費エネルギーがどのように変化するかを直接求めることができなかった。
　このように、従来は、設備のエネルギー効率や制約条件等の特性を変更することによってどれだけ省エネが実現できるかを把握するためには、様々な設備の特性の値を具体的に設定してそれぞれについて最適化の計算を実施して、計算結果同士を比較していく必要があり、試行錯誤によるしか手段がなかった。
　本発明は、プラントを構成する設備の特性を最適化する際に、何度も計算を繰り返すことを不要としつつ、どのように設備の特性を設定すべきか等の情報をユーザに提示することのできる技術を提供することを目的とする。
　本発明の一態様に係るエネルギー解析支援装置は、複数の設備からなるプラントのエネルギー解析支援装置であって、下記の構成を有する。
・前記プラントのモデルと、前記プラントを構成する設備の能力を表す設備モデルと、前記プラントの消費エネルギーまたは負荷を表す外部条件とが与えられると、前記プラントのモデルを表す数式、前記設備モデルを表す数式及び外部条件を表す数式より、解析対象の設備の特性を表した少なくとも一の変数を含む論理式を生成する計算部；及び
・前記計算部により生成された論理式に基づいて、前記一の変数が表す前記解析対象の設備の特性を２軸のうちの一の軸とするグラフを生成して出力する出力部。
　また、本発明の一態様に係るエネルギー解析支援プログラムは、複数の設備からなるプラントのエネルギー解析支援処理を情報処理装置に実行させるためのエネルギー解析支援プログラムであって、下記の構成を有する。
・前記プラントのモデルと、前記プラントを構成する設備の能力を表す設備モデルと、前記プラントの消費エネルギーまたは負荷を表す外部条件とが与えられると、前記プラントのモデルを表す数式、前記設備モデルを表す数式及び外部条件を表す数式より、解析対象の設備の特性を表した少なくとも一の変数を含む論理式を生成するステップ；及び
・前記生成された論理式に基づいて、前記一の変数が表す前記解析対象の設備の特性を２軸のうちの一の軸とするグラフを生成して出力するステップ。
　また、本発明の一態様に係る記録媒体は、複数の設備からなるプラントのエネルギー解析支援処理を情報処理装置に実行させるためのエネルギー解析支援プログラムを記録した記録媒体であって、下記の構成を有する。
・前記プラントのモデルと、前記プラントを構成する設備の能力を表す設備モデルと、前記プラントの消費エネルギーまたは負荷を表す外部条件とが与えられると、前記プラントのモデルを表す数式、前記設備モデルを表す数式及び外部条件を表す数式より、解析対象の設備の特性を表した少なくとも一の変数を含む論理式を生成するステップ；及び
・前記生成された論理式に基づいて、前記一の変数が表す前記解析対象の設備の特性を２軸のうちの一の軸とするグラフを生成して出力するステップ。
　また、本発明の一態様に係るエネルギー解析支援方法は、複数の設備からなるプラントのエネルギー解析支援方法であって、下記の構成を有する。
・前記プラントのモデルと、前記プラントを構成する設備の能力を表す設備モデルと、前記プラントの消費エネルギーまたは負荷を表す外部条件とが与えられると、前記プラントのモデルを表す数式、前記設備モデルを表す数式及び外部条件を表す数式より、解析対象の設備の特性を表した少なくとも一の変数を含む論理式を生成するステップ；及び
・前記生成された論理式に基づいて、前記一の変数が表す前記解析対象の設備の特性を２軸のうちの一の軸とするグラフを生成して出力するステップ。
　本発明によれば、プラントを構成する設備の特性を最適化する際に、何度も計算を繰り返すことを不要としつつ、どのように設備の特性を設定すべきか等の情報をユーザに提示することが可能となる。
第１の実施形態に係るエネルギー解析支援装置の構成図である。 第１の実施形態に係る計算部における動作フローを示す図である。 第１の実施形態において一階述語論理式を生成する方法及び総消費エネルギー数式を計算する方法を説明する図である。 第１の実施形態に係る出力部が生成するグラフを模式的に示す図である。 プラントのシステムモデルを例示する図である。 設備モデルを例示する図である。 外部条件を例示する図である。 計算部が図５～図７のシステムモデル、設備モデル及び外部条件から生成する数式群を例示する図である。 第１の実施形態における一階述語論理式を生成する方法を説明する図である。 一階述語論理式から論理式を得る方法について説明する図である。 第１の実施形態において表示部に表示するグラフを例示する図である。 変形例１において生成する設備モデル数式を例示する図である。 変形例１において生成される一階述語論理式を示す図である。 図１３の一階述語論理式に限定記号消去アルゴリズムを適用して得られる論理式を示す図である。 変形例１において表示部に表示するグラフを例示する図である。 変形例２において生成する外部条件数式を説明する図である。 変形例２において生成される２つの一階述語論理式を示す図である。 図１７の一階述語論理式に限定記号消去アルゴリズムを適用して得られる論理式を示す図である。 変形例２において表示部に表示するグラフを例示する図である。 変形例３において生成する設備モデル数式を説明する図である。 変形例３において生成される一階述語論理式を示す図である。 図２１の一階述語論理式に限定記号消去アルゴリズムを適用して得られる論理式を示す図である。 変形例３において表示部に表示するグラフを例示する図である。 第２の実施形態に係る計算部における動作フローを示す図である。 第２の実施形態において計算部が生成する一階述語論理式及び一階述語論理式から限定記号を消去して得られる論理式について説明する図である。 第２の実施形態に係る出力部が生成するグラフを模式的に示す図である。 外部条件を例示する図である。 計算部１１が図５のシステムモデル、図６の設備モデル及び図２７の外部条件から生成する数式群を例示する図である。 第２の実施形態における一階述語論理式を生成する方法を説明する図である。 図２９の一階述語論理式に限定記号消去アルゴリズムを適用して得られる論理式を示す図である。 第２の実施形態において表示部に表示する画面の一例を示す図である。 第２の実施形態において表示部に表示する画面の他の例を示す図である。 情報処理装置の構成図である。 記録媒体を示す図である。
　以下、本発明の実施の形態について、図面を参照して詳細に説明する。
＜第１の実施形態＞
　図１は、本実施形態に係るエネルギー解析支援装置の構成図である。図１のエネルギー解析支援装置１は、計算部１１と、出力部１２とを有し、キーボード２Ａやポインティングデバイス２Ｂ等の入力部２やモニタ等の表示部５と接続される。
　エネルギー解析支援装置１は、入力部２（２Ａや２Ｂ）を介して、プラントのシステムモデル、設備モデル及び外部条件の入力を受け付ける。エネルギー解析支援装置１の計算部１１は、入力されたこれらの情報より、システムモデル数式、設備モデル数式及び外部条件数式を生成し、生成した数式を用いて一階述語論理式を生成する。計算部１１は、一階述語論理式に対し限定記号消去アルゴリズムを適用して、所望の変数と総消費エネルギーとの関係を表した論理式を得る。出力部１２は、計算部１１が計算により得た論理式に基づきグラフを生成してこれを表示部５に出力する。表示部５は、出力部５から入力された情報に基づき、グラフを表示する。
　こうして、エネルギー解析支援装置１は、解析対象の設備のエネルギー効率やプラント運転時の制約条件等の特性を変更した場合に、プラントの総消費エネルギーがどのように変化するかをユーザが解析する際にこれを支援するための情報を提供する。
　なお、エネルギー解析支援装置１は、入力部２を介してシステムモデル等の情報が入力されたことを契機として、入力された情報を用いて処理を開始することもある。あるいは、既に装置内に記憶されている情報を用いて、ユーザからの開始の指示にしたがって処理を開始することもある。
　図２は、計算部１１における動作フローを示す図である。
　計算部１１は、まず、ステップＳ１で、システムモデル数式、設備モデル数式及び外部条件数式を取得する。これらの数式群は、システムモデル、設備モデル及び外部条件を用いて計算部１１が生成してもよいし、入力部２を介して直接入力されてもよい。
　ステップＳ２では、ステップＳ１で取得した数式群より、総消費エネルギーに関する一階述語論理式を生成する。
　ステップＳ３では、ステップＳ２で生成した一階述語論理式を処理して総消費エネルギー数式を得る。最後に、ステップＳ４で、総消費エネルギー数式等の情報を出力部１２に出力し、処理を終了する。
　図３は、図２のステップＳ２の一階述語論理式を生成する方法及びステップＳ３の総消費エネルギー数式を計算する方法を説明する図である。
　計算部１１は、システムモデル数式Ｆ１、設備モデル数式Ｆ２及び外部条件数式Ｆ３を論理積（∧）で結合し、図３の数式（３－２）を生成する。そして、数式（３－２）に対し、更に、所望の設備の特性を表す変数ａ及びプラントの総消費エネルギーを表す変数Ｐを除く変数に対して存在記号（∃）を付与する（図３の（３－１））。こうして、一階述語論理式ψを生成する。
ψ:=∃P1,...,∃PN∃L∃L1,..., ∃LN(F1∧F2∧F3)
Ｌは、プラントの総負荷を表す変数である。また、Ｎは、設備の台数を表し、Ｐｎは、設備ｎ（ｎ＝１，２，…，Ｎ）の消費エネルギーを、Ｌｎは、設備ｎの負荷を表す変数である。
　計算部１１は、生成した一階述語論理式ψに限定記号消去アルゴリズムを適用して、設備の特性を表す変数ａと総消費エネルギーを表す変数Ｐとが満たす関係を表した論理式φを得て、これを出力部１２に出力する。
φ:=QE[ψ]
(φ:=(Pとaに関する数式))
　出力部１２は、計算部１１から入力された論理式φより、設備の特性を表す変数ａとプラントの総消費エネルギーを表す変数Ｐとの関係を表すグラフを生成する。そして、出力部１２は、生成したグラフ等の情報すなわちエネルギー解析支援装置１において解析支援処理を行った結果得られた情報を表示部５に出力する。表示部５は、出力部１２から入力された情報に基づき、エネルギー解析支援装置１における計算結果を表示する。
　図４は、出力部１２が生成するグラフを模式的に示す図である。図４に例示するグラフでは、一方の軸（横軸）には、設備の特性を表す変数ａをとり、他方の軸（縦軸）には、プラントの総消費エネルギーを表す変数Ｐをとる。論理式φが満たす領域Ａが可視化される。
　以下、本実施形態に係るエネルギー解析支援装置１が、入力された各種のモデル等の情報を用いてどのように各種数式の計算を行い、計算結果のグラフを得るかについて、具体例を用いて説明する。まず、図５～図７を参照して、エネルギー解析支援装置１に入力される各種のモデル等の情報について説明する。
　図５は、プラントのシステムモデルを例示する図である。
　図５の空調システムは、複数の冷凍機を有し、複数の冷凍機は、電源から供給される電力を用いて、空調対象空間が需要する熱負荷を供給する。システムモデルは、図５に示す空調システムの設備構成や各設備の熱負荷及び消費電力から構成される。
　具体的には、冷凍機＃１は、熱負荷Ｌ１［ｋＷ］に対して電力Ｐ１［ｋＷ］を消費する。冷凍機＃２は、熱負荷Ｌ２［ｋＷ］に対し電力Ｐ２［ｋＷ］を消費し、冷凍機＃３は、熱負荷Ｌ３［ｋＷ］に対し電力Ｐ３［ｋＷ］を消費する。
　設例では、更に１台の冷凍機＃４を新たに導入するとする。冷凍機＃４は、熱負荷Ｌ４［ｋＷ］に対して電力Ｐ４［ｋＷ］を消費する。
　また、各冷凍機＃１～＃４が消費する消費電力Ｐ１～Ｐ４の総和は、システム全体で消費する総消費電力Ｐ［ｋＷ］となり、空調対象空間が要求する総熱負荷Ｌ［ｋＷ］を各冷凍機＃１～＃４の熱負荷Ｌ１～Ｌ４［ｋＷ］として配分する。
　図６は、設備モデルを例示する図である。
　図６の設備モデルは、既存の設備である冷凍機＃１～＃３についての熱負荷に対する消費電力特性から構成される。各冷凍機＃１～＃３は、それぞれ異なる消費電力特性と上下限の制約をもつ。
　図７は、外部条件を例示する図である。
　本実施形態においては、外部条件として、空調対象空間の要求に対応するための空調システム全体の総熱負荷Ｌを、５００［ｋＷ］と設定する。
　次に、入力されたシステムモデル、設備モデル及び外部条件より一階述語論理式を生成し、これを処理して総消費エネルギー数式を得る方法について説明する。
　図８は、計算部１１が図５～図７のシステムモデル、設備モデル及び外部条件から生成する数式群を例示する図である。
　システムモデル数式Ｆ１は、図５の空調システムモデルを表現した式であり、２つの数式（８－１）、（８－２）を論理積で結合した式からなる。
F1:=P=P1+P2+P3+P4∧L=L1+L2+L3+L4
図８の式（８－１）は、各冷凍機＃１～＃４の消費電力Ｐ１～Ｐ４の総和が、電源が各冷凍機＃１～＃４に供給する総消費電力Ｐに等しいことを表す。式（８－２）は、空調対象空間の総熱負荷Ｌが各冷凍機＃１～＃４に配分されることを表す。
　設備モデル数式Ｆ２は、図６の設備モデルを表現した数式であり、４台の冷凍機＃１～＃４それぞれの特性を表す数式を、論理積で結合した式からなる。下記の数式Ｆ２のうち、４行目は、新たに導入する冷凍機＃４の特性を表す数式である。冷凍機＃４の特性を表す数式に含まれる変数ａは、冷凍機＃４の負荷に対する消費電力特性を表す。
F2:=(P1=0.12*L1+16∧69<=L1<=125∨P1=0∧L1=0)
∧(P2=0.27*L2+21∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212∨P3=0∧L3=0)
∧(P4=a*L4+15∧69<=L4<=212∧0.1<=a<=0.35∨P4=0∧L4=0)
　外部条件数式Ｆ３は、図７の外部条件を表現した数式であり、先に図７の説明においても述べたとおり、空調システム全体の総熱負荷Ｌを表す。
F3:=L=500
　計算部１１は、システムモデル数式Ｆ１、設備モデル数式Ｆ２及び外部条件数式Ｆ３を図１においては不図示のメモリ等の記憶手段に保持し、保持しているこれらの数式を用いて、以下の数式処理を実施する。
　図９は、図８の数式群Ｆ１～Ｆ３から一階述語論理式ψを生成する方法を説明する図である。
　エネルギー解析支援装置１の計算部１１は、まず、図８の数式群、すなわち、システムモデル数式Ｆ１、設備モデル数式Ｆ２及び外部条件数式Ｆ３を論理積で結合し、システム全体を表現した数式（Ｆ１∧Ｆ２∧Ｆ３）を生成する（図９の式（９－２））。
　そして、計算部１１は、システム全体を表現する数式Ｆ１∧Ｆ２∧Ｆ３に対して、総消費エネルギーを表す変数Ｐと冷凍機＃４の特性を表す変数ａとを除く変数に存在記号（∃）を付与し（図９の（９－１））、総消費電力に関する一階述語論理式ψを生成する。
ψ:=∃P1∃P2∃P3∃P4∃L∃L1∃L2∃L3∃L4（F1∧F2∧F3）
F1∧F2∧F3:=
(P=P1+P2+P3+P4∧L=L1+L2+L3+L4)
∧(P1=0.12*L1+16∧69<=L1<=125∨P1=0∧L1=0)
∧(P2=0.27*L2+21∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212∨P3=0∧L3=0)
∧(P4=a*L4+15∧69<=L4<=212∧0.1<=a<=0.35∨P4=0∧L4=0)
∧(L=500)
　図１０は、一階述語論理式ψから論理式φを得る方法について説明する図である。
　計算部１１は、先に図３を参照して説明したとおり、図９の一階述語論理式ψに限定記号消去アルゴリズム（ＱＥ［］）を適用して、論理式ψから存在記号の付与された変数を消去した論理式φを得る。論理式φは、冷凍機＃４の特性を表す変数ａと、空調システムの総消費電力を表す変数Ｐとが満たす関係式からなる。
φ:=QE[ψ]
(φ:=(-35+a<=0∧1/10-a<=0∧14263/100+212*a-P=0∧-P<0)∨(-35+a<=0∧1/10-a<=0∧14263/100+212*a-P=0∧12763/100-P<0)∨(-35+a<=0∧1/10-a<=0∧15083/100+212*a-P<=0∧-3888/25-212*a+P<=0∧-P<0)∨...(以下省略))
論理式φは、冷凍機＃４の消費電力特性ａに対して各冷凍機＃１～＃４の負荷配分を変えることにより実行可能な総消費電力Ｐの範囲を表している。
　図１１は、表示部５に表示するグラフを例示する図である。
　上記のとおり、エネルギー解析支援装置１の計算部１１は、一階述語論理式ψから限定記号を消去して、図１０の論理式φを得る。得られた論理式φは、エネルギー解析支援装置１の出力部１２に渡される。
　出力部１２は、論理式φより、縦軸に総消費電力を表す変数Ｐを、横軸に冷凍機＃４の消費電力特性を表す変数ａをとり、論理式φが満たす領域Ａ１を表したグラフを生成する。図１１においては、領域Ａ１を斜線で示す。出力部１２は、生成したグラフの画像情報を、表示部５に出力する。表示部５は、出力部１２から入力された情報に基づき図１１のグラフを表示する。
　図１１のグラフによれば、領域Ａ１の縦軸方向の下限値が、冷凍機＃４の消費電力特性ａの値に対する最適な消費電力である。下限ＬＭ＿ｌｗは、冷凍機＃４の消費電力特性ａの値に対する最適な総消費電力に相当する。ユーザは、図１１のグラフの下限ＬＭ＿ｌｗを参照して、最も省電力を実現する運用をした場合における空調システム全体の総消費電力Ｐを把握しつつ、導入する冷凍機＃４の消費電力特性ａを決定することができる。したがって、ユーザは、グラフを参照して新たに導入する設備（冷凍機＃４）の特性（設例では消費電力特性）を適切に決定することができ、これにより、新たに導入する設備を適切に決定することが可能となる。
　以上説明したように、本実施形態に係るエネルギー解析支援装置１によれば、入力されたシステムモデル、設備モデル及び外部条件より、一階述語論理式を生成する。そして、生成した一階述語論理式に限定記号消去アルゴリズムを適用して、プラントのある設備の特性とプラントの総消費電力との関係を表した論理式を得る。一階述語論理式は、入力されたシステムモデル等から生成したプラントの全体を表現した数式に、最終的に処理結果として得たい変数以外の変数、すなわち、ある設備の特性と総消費電力以外の変数に存在記号（∃）を付与することにより生成する。こうして生成した一階述語論理式から限定記号を消去して得られる論理式をグラフにしてモニタ等の表示部５に表示をすることで、ユーザにとっては、新たに導入する設備の特性と総消費エネルギーとが満たす関係を視覚的に把握することが可能となり、有効な解析に資する。
　ある設備の変更は、システムを構成する他の設備の運用にも影響を与える。このため、単に効率のよい設備を導入しても、十分な省エネによる投資回収効果が見込めるか否かが判然とせず、必要以上に高性能で効果な設備を導入しても、その能力が活かしきれずに無駄な投資となってしまう場合がある。しかし、本実施形態に係るエネルギー解析支援装置１を用いて解析を行うことで、設備の導入前の検証によりそのような事態の発生を回避することが可能となる。
（変形例１）
　上記の実施形態においては、ある設備を導入する場合に、その設備の特性と総消費エネルギーとの関係を表す論理式φを求めてこれをグラフに描画している。上記の実施形態により描画されるグラフによれば、設備の特性と総消費エネルギーとが満たす領域には、導入する設備を稼動させる場合に満たされる領域と稼動させない場合のそれとの両方を含んでいる。これに対し、本変形例では、ある設備を稼働させる場合及び稼働させない場合のそれぞれの場合を区別したグラフを描画する点で異なる。
　以下、本変形例について、上記の実施形態と異なる点を中心に説明する。なお、エネルギー解析支援装置１の構成、エネルギー解析支援装置１に入力される情報（システムモデル、設備モデル及び外部条件）については、上記の実施形態と同様であり、図１や図５～図７を参照して説明したとおりである。
　エネルギー解析支援装置１は、図５～図７の情報が入力されると、上記の実施形態と同様に、まず、システムモデル数式、設備モデル数式及び外部条件数式を生成する。システムモデル数式Ｆ１及び外部条件数式Ｆ３については、先の実施形態と同様であり、図８を参照して説明したとおりである。
　図１２は、本変形例において生成する設備モデル数式を例示する図である。
　本変形例においては、エネルギー解析支援装置１の計算部１１は、設備モデル数式として、新たに空調システムに導入する冷凍機＃４を稼働しない場合と、稼働させる場合とでそれぞれ数式Ｆ２’、Ｆ２’’を生成する。
F2’:=(P1=0.12*L1+16∧69<=L1<=125∨P1=0∧L1=0)
∧(P2=0.17*L2+21∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212∨P3=0∧L3=0)
∧(P4=0∧L4=0)
F2’’:=(P1=0.12*L1+16∧69<=L1<=125∨P1=0∧L1=0)
∧(P2=0.17*L2+21∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212∨P3=0∧L3=0)
∧(P4=a*L4+15∧69<=L4<=212∧0.1<=a<=0.35)
　設備モデル数式Ｆ２’は、図８の設備モデル数式Ｆ２から、冷凍機＃４を稼働する条件を表した以下の数式（１２－２）を除いた数式で表される。
(P4＝ａ*L4+15∧69<=L4<=212∧0.1<=a<=0.35)…（１２－２）
　設備モデル数式Ｆ２’’は、図８の設備モデル数式Ｆ２から、冷凍機＃４を稼働しない条件を表した以下の数式（１２－１）を除いた数式で表される。
(P4=0∧L4=0)…（１２－１）
　計算部１１は、生成した２つの設備モデル数式Ｆ２’、Ｆ２’’のそれぞれを用いて、２つの一階述語論理式を生成する。
　図１３は、図８の数式Ｆ１、Ｆ３と図１２の数式Ｆ２’、Ｆ２’’とを用いて生成される一階述語論理式ψ２’、ψ２’’を示す図である。
　冷凍機＃４が稼働しない場合の総消費電力に関する一階述語論理式ψ２’に関しては、システム全体を表現した数式として、（Ｆ１∧Ｆ２’∧Ｆ３）を生成する。冷凍機＃４が稼働する場合の一階述語論理式ψ２’’に関しては、（Ｆ１∧Ｆ２’’∧Ｆ３）を生成する。存在記号（∃）を付与する変数については、いずれの式ψ２’、ψ２’’についても、上記の実施形態と同様である。
ψ2’:= ∃P1∃P2∃P3∃P4 L∃L1∃L2∃L3∃L4(F1∧F2’∧F3)
(F1∧F2’∧F3:=
P=P1+P2+P3+P4∧L=L1+L2+L3+L4
∧(P1=0.12*L1+ 16∧69<=L1<=125 ∨P1=0∧L1=0)
∧(P2=0.17*L2+21 ∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212 ∨P3=0∧L3=0)
∧(P4=0∧L4=0)
∧L=500))
ψ2’’:=∃P1∃P2∃P3∃P4 ∃L∃L1∃L2∃L3∃L4(F1∧F2’’∧F3)
(F1∧F2’’∧F3:=
P=P1+P2+P3+P4∧L=L1+L2+L3+L4
∧(P1=0.12*L1+ 16∧69<=L1<=125 ∨P1=0∧L1=0)
∧(P2=0.17*L2+21 ∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212 ∨P3=0∧L3=0)
∧(P4=a*L4+15∧69<=L4<=212 ∧0.1<=a<=0.35)
∧L=500))
　計算部１１は、生成した２つの一階述語論理式ψ２’、ψ２’’のそれぞれについて限定記号消去アルゴリズムを適用し、冷凍機＃４の消費電力特性を表す変数ａと総消費電力を表す変数Ｐとの関係を表した論理式φ２’、φ２’’を得る。
　図１４は、一階述語論理式ψ２’、ψ２’’に限定記号消去アルゴリズムを適用して得られる論理式φ２’、φ２’’を示す図である。
φ2’:=QE[ψ2’]
(φ2’:= 623/4<=P<=15719/100)
φ2’’:=QE[ψ2’’]
(φ2’’:=(-35+a<=0∧1/10-a<=0∧3072/25+163*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧14263/100+113*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧12103/100+203*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧12903/100+193*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧387/4+200*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧2778/25+212*a-P=0∧-P<0)．．．(以下省略))
　論理式φ２’、φ２’’は、それぞれ空調システムに新たに導入する冷凍機＃４を稼働しない場合及び稼働する場合の冷凍機＃４の消費電力特性ａと総消費電力Ｐとの関係を表す。出力部１２は、計算部１１から式（１４－１）、（１４－２）を受け取ると、これに基づき２つの式が満たす領域を比較可能なグラフを生成する。そして、出力部１２は、これをエネルギー解析支援装置１において解析支援処理を実施した結果得られた情報として、表示部５に出力する。表示部５は、入力された情報に基づき、グラフを表示する。
　図１５は、表示部５に表示するグラフを例示する図である。図１１のグラフと同様に、縦軸には、総消費電力を表す変数Ｐを、横軸には、冷凍機＃４の消費電力特性を表す変数ａをとる。
　図１５に示すように、本変形例では、２つの論理式φ２’、φ２’’のそれぞれが満たす領域Ａ２’、Ａ２’’が色分け等により、互いに識別可能に表示される。図１５においては、領域Ａ２’、Ａ２’’を斜線、網掛けで識別可能とするグラフを作成する。
　ユーザは、図１５に例示するグラフを参照して、新たに導入する冷凍機＃４が稼働しない場合及び稼働する場合のそれぞれの領域を比較して、解析を行う。
　例えば、冷凍機＃４の消費電力特性ａの値が小さい場合には、冷凍機＃４を稼働させて他の効率の悪い冷凍機を停止させることで、稼働させない場合、すなわち冷凍機＃１～＃３を稼働させる場合よりも大幅に総消費電力Ｐを抑制できることがわかる（図１５の※１）。
　また、消費電力特性ａの値が約０．２９５以上の領域では、冷凍機＃４を稼働させない方が総消費電力Ｐ小さくできること、すなわち、新たに冷凍機＃４を導入しても、総消費電力Ｐを抑制できないことがわかる（図１５の※２）。
　このように、本変形例によれば、新たに導入する設備を稼働するか否かが総消費エネルギーに与える影響や、新たに導入する設備の特性に応じて、省エネの観点からは新たに導入する設備または既存の設備のいずれを稼働させるべきかを容易に判断可能となる。
（変形例２）
　上記の実施形態及び変形例１では、外部条件として、空調対象空間が相対的に高熱負荷である場合を設定している。これに対し、本変形例では、空調対象空間が相対的に低熱負荷である場合を設定して、低熱負荷時における導入する設備の特性と総消費エネルギーとの関係を表す論理式φを求めてグラフに描画する。
　以下、本変形例について、上記の実施形態及び変形例１と異なる点を中心に説明する。エネルギー解析支援装置１の構成、入力される情報（システムモデル、設備モデル及び外部条件）については、上記の実施形態及び変形例１と同様であり、図１や図５～図７を参照して説明したとおりである。
　エネルギー解析支援装置１は、図５～図７の情報が入力されると、上記の実施形態や変形例１と同様に、システムモデル数式、設備モデル数式及び外部条件数式を生成する。本変形例では、変形例１と同様に、新たに導入する冷凍機＃４を稼働する場合及び稼働しない場合のそれぞれについて設備モデル数式Ｆ２’、Ｆ２’’を生成して、それぞれ冷凍機＃４の消費電力特性ａと総消費電力Ｐとの関係式を求める場合について説明する。
　システムモデル数式Ｆ１については、上記の実施形態と同様であり、図８を参照して説明したとおりである。設備モデル数式Ｆ２’、Ｆ２’’については、変形例１で図１２を参照して説明したとおりである。
　図１６は、本変形例において生成する外部条件数式を説明する図である。
　本変形例では、エネルギー解析支援装置１の計算部１１は、低熱負荷時の外部条件数式として、
F3’:=L=400
を生成する。このように、本変形例では熱負荷Ｌには、図１６上側の時間帯別の空調対象空間の総熱負荷のグラフでは「６時～７時」「７時～８時」等の時間帯に設定されている、相対的に低い総熱負荷の値である４００［ｋＷ］を設定する。
　計算部１１は、システムモデル数式Ｆ１、２つの設備モデル数式Ｆ２’とＦ２’’、及び低熱負荷時の外部条件数式Ｆ３’を用いて、２つの一階述語論理式を生成する。
　図１７は、生成される２つの一階述語論理式ψ３’、ψ３’’を示す図である。
　図１７の２つの一階述語論理式ψ３’及びψ３’’は、それぞれ冷凍機＃４を稼働しない場合及び稼働する場合の総消費電力に関する一階述語論理式である。計算部１１は、図８のシステムモデル数式Ｆ１、図１２の設備モデル数式Ｆ２’とＦ２’’、図１６の外部条件数式Ｆ３’からそれぞれシステム全体を表現した式（Ｆ１∧Ｆ２’∧Ｆ３’）、（Ｆ１∧Ｆ２’’∧Ｆ３’）を生成する。存在記号（∃）を付与する変数は、上記の実施形態や変形例１と同様である。
ψ3’:= ∃P1∃P2∃P3∃P4 L∃L1∃L2∃L3∃L4((F1∧F2’∧F3’)
F1∧F2’∧F3’:=
P=P1+P2+P3+P4∧L=L1+L2+L3+L4
∧(P1=0.12*L1+ 16∧69<=L1<=125 ∨P1=0∧L1=0)
∧(P2=0.17*L2+21 ∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212 ∨P3=0∧L3=0)
∧(P4=0∧L4=0)
∧L=400)
ψ3’’:=∃P1∃P2∃P3∃P4 ∃L∃L1∃L2∃L3∃L4(F1∧F2’’∧F3’)
F1∧F2’’∧F3’:=
P=P1+P2+P3+P4∧L=L1+L2+L3+L4
∧(P1=0.12*L1+ 16∧69<=L1<=125 ∨P1=0∧L1=0)
∧(P2=0.17*L2+21 ∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212 ∨P3=0∧L3=0)
∧(P4=a*L4+15∧69<=L4<=212 ∧0.1<=a<=0.35)
∧L=400)
　先の実施形態等においても説明したように、計算部１１は、生成した一階述語論理式に限定記号消去アルゴリズムを適用し、冷凍機＃４の消費電力特性を表す変数ａと総消費電力を表す変数Ｐとの関係を表した論理式を得る。
　図１８は、一階述語論理式ψ３’、ψ３’’に限定記号消去アルゴリズムを適用して得られる論理式φ３’、φ３’’を示す図である。
φ3’:=QE[ψ3’]
（φ3’:=13329/100<=P<=14239/100）
φ3’’:=QE[ψ3’’]
(φ3:=(-35+a<=0∧1/10-a<=0∧2297/25+188*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧2532/25+153*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧2904/25+119*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧2364/25+209*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧12103/100+103*a-P=0∧-P<0)
∨(-35+a=0∧1/10-a<=0∧12903/100+93*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧10743/100+183*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧387/4+100*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧1663/20+180*a-P=0∧-P<0)
∨(-35+a<=0∧1/10-a<=0∧9003/100+156*a-P=0∧-P<0) ．．．(以下省略))
　論理式φ３’、φ３’’は、それぞれ低熱負荷時において冷凍機＃４を稼働しない場合と稼働する場合の冷凍機＃４の消費電力特性ａと総消費電力Ｐとが満たす関係式である。出力部１２は、計算部１１から論理式φ３’、φ３’’を受け取ると、これに基づき２つの式が満たす領域を比較可能なグラフを生成し、これを表示部５に出力する。
　図１９は、表示部５に表示するグラフを例示する図である。図１１や図１５のグラフと同様に、縦軸は、総消費電力Ｐを、横軸は、冷凍機＃４の消費電力特性ａをとる。
　図１９に示すように、本変形例では、２つの論理式φ３’、φ３’’のそれぞれが満たす領域Ａ３’を網掛けで、Ａ３’’を斜線で識別可能とするグラフを作成する。ユーザは、図１９に例示するグラフを参照して、低熱負荷時に新たに導入する冷凍機＃４を稼働しない場合及び稼働する場合のそれぞれの領域を比較して、解析を行う。
　例えば、図１９においては矢印で指し示す領域のように、冷凍機＃４の消費電力特性ａの範囲が０．１０から０３５の領域（図１５に示すグラフの全域）では、冷凍機＃４を稼働させた方が総消費電力Ｐを抑制できることがわかる。
　なお、図１９においては、低熱負荷時における冷凍機＃４の消費電力特性ａと総消費電力Ｐとの関係を、冷凍機＃４を稼働しない場合と稼働する場合とで識別可能なグラフを例示している。この他にも、例えば変形例１で求めた図１４の論理式φ２’、φ２’’を用いて、高熱負荷時と低熱負荷時とを比較するグラフ等、更に各種の変形例が考えられる。
　このように、本変形例によれば、外部条件を変更した場合であっても、新たに導入する設備の特性と総消費エネルギーとの関係がどのようになるかについて容易に判断可能となる。
（変形例３）
　上記の実施形態及び変形例１、２では、設備（冷凍機）の特性の一つである消費電力特性を表す変数ａと総消費電力Ｐとの関係を求めている。これに対し、本変形例では、新たに導入する冷凍機の消費電力特性以外の特性を表す変数と総消費電力Ｐとの関係を求める。
　以下、本変形例について、上記の実施形態等と異なる点を中心に説明する。エネルギー解析支援装置１の構成、入力される情報については、上記の実施形態等と同様であり、図１や図５～図７を参照して説明したとおりである。
　エネルギー解析支援装置１は、図５～図７の情報が入力されると、上記の実施形態等と同様に、システムモデル数式、設備モデル数式及び外部条件数式を生成する。システムモデル数式Ｆ１と外部条件数式Ｆ３については、いずれも図８を参照して説明したとおりである。
　図２０は、本変形例において生成する設備モデル数式を説明する図である。
　本変形例では、冷凍機＃４の熱負荷の上限を表す変数ｂを用いて、冷凍機＃４を稼働する場合における設備モデル数式Ｆ２’’’を生成する。冷凍機＃４の消費電力特性については、固定値０．１６とする。
F2’’’:=(P1=0.12*L1+ 16∧69<=L1<=125∨P1=0∧L1=0)
∧(P2=0.17*L2+21∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212∨P3=0∧L3=0)
∧(P4=0.16*L4+15∧69<=L4<=b∧69<=b<=212)
上の設備モデル数式Ｆ２’’’のうち、冷凍機＃４の特性を表す数式は、以下のとおりである。
(P4=0.16*L4+15∧69<=L4<=b∧69<=b<=212)
このうち、変数ｂを含む第２の式「69<=L4<=b」は、熱負荷Ｌ４が下限値６９以上且つ上限値ｂ以下であることを表す。第３の式「69<=b<=212」は、熱負荷Ｌ４の上限ｂが６９以上かつ２１２以下の値をとることを表す。
　計算部１１は、システムモデル数式Ｆ１、設備モデル数式Ｆ２’’’及び外部条件数式Ｆ３を用いて、一階述語論理式を生成する。
　図２１は、生成される一階述語論理式ψ４’’を示す図である。
　図２１の一階述語論理式ψ４’’は、冷凍機＃４の熱負荷の上限ｂを変数とする一階述語論理式である。計算部１１は、図８のシステムモデル数式Ｆ１、図２０の設備モデルＦ２’’’及び図８の外部条件数式からシステム全体を現した式（Ｆ１∧Ｆ２’’’∧Ｆ３）を生成し、変数ｂ及びＰ以外の変数に存在記号（∃）を付与して、論理式ψ４’’を生成する。
ψ4’’:=∃P1∃P2∃P3∃P4 ∃L∃L1∃L2∃L3∃L4(F1∧F2’’’∧F3)
(F1∧F2’’’∧F3:=
P=P1+P2+P3+P4∧L=L1+L2+L3+L4
∧(P1=0.12*L1+ 16∧69<=L1<=125∨P1=0∧L1=0)
∧(P2=0.17*L2+21 ∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212∨P3=0∧L3=0)
∧(P4=0.16*L4+15∧69<=L4<=b∧69<=b<=212)
∧L=500))
　計算部１１は、生成した一階述語論理式に限定記号消去アルゴリズムを適用し、冷凍機＃４の熱負荷Ｌ４の上限を表す変数ｂと総消費電力を表す変数Ｐとの関係を表した論理式を得る。
　図２２は、一階述語論理式ψ４’’に限定記号消去アルゴリズムを適用して得られる論理式φ４’’を示す図である。
φ4’’:=QE[ψ4’’]
（φ4’’:=(-212+b<=0∧163-b<=0∧-3724/25+P=0)
∨(-212+b<=0∧200-b<=0∧-515/4+P=0)
∨(-212+b<=0∧203-b<=0∧-15351/100+P=0)
∨(-212+b<=0∧113-b<=0∧-16071/100+P=0)
∨(-212+b<=0∧193-b<=0∧-15991/100+P=0)
∨(-212+b<=0∧69-b<=0∧-375-4*b+25*P<=0∧-95+P=0)
∨(-212+b<=0∧69-b<=0∧-3439/20+P<=0∧16523/100-P<=0)
∨(-212+b<=0∧69-b<=0∧-17083/100+P<=0∧16523/100-P<=0)
∨(-212+b<=0∧69-b<=0∧-1737/10+P<=0∧3419/20-P<=0)
∨(-212+b<=0∧69-b<=0∧-3419/20+P<=0∧17083/100-P<=0)
∨(-212+b<=0∧69-b<=0∧-1737/10+P<=0∧3439/20-P<=0)
∨(-212+b<=0∧163-b<=0∧3561+b-25*P=0∧-P<0)
∨(-212+b<=0∧163-b<=0∧-P<0∧-4050+2*b+25*P= 0)．．．（以下省略））
　論理式φ４’’は、冷凍機＃４の熱負荷Ｌ４の上限ｂと総消費電力Ｐとの関係を表す。出力部１２は、計算部１１から論理式φ４’’を受け取ると、これに基づき論理式φ４’’が満たす領域を表すグラフを生成し、これを表示部５に出力する。
　図２３は、表示部５に表示するグラフを例示する図である。図１１、１５及び図１９のグラフと同様に、縦軸は、総消費電力Ｐをとる。横軸では、冷凍機＃４の熱負荷Ｌ４の上限ｂをとる。
　図２３に示すグラフのうち、領域Ａ４’’は、論理式φ４’’の変数ｂとＰとが満たす領域を斜線で表している。図２０の設備モデル数式Ｆ２’’’を生成する方法において記載したとおり、領域Ａ４’’は、冷凍機＃４を稼働する場合における冷凍機＃４の熱負荷Ｌ４の上限ｂと総消費電力Ｐとが満たす領域である。
　これに対し、図２３の領域Ａ４’は、冷凍機＃４を稼働しない場合における冷凍機＃４の熱負荷Ｌ４の上限ｂと総消費電力Ｐとが満たす領域である。領域Ａ４’は、上記の変形例２にて冷凍機＃４を稼働しない場合について求めた論理式φ２’を利用している。論理式φ２’は、一階述語論理式ψ２’から限定記号を消去することにより得られ、一階述語論理式ψ２’は、システムモデル数式Ｆ１、設備モデル数式Ｆ２’及び外部条件数式Ｆ３を用いて生成している。
　ここで、図１２の式（１２－１）に示すように、設備モデル数式Ｆ２’は、冷凍機＃４の熱負荷Ｌ４の上限ｂには依存しない。このため、設備モデルＦ２’を用いて求めた論理式φ２’もまた、変数ｂによらず、したがって、論理式φ２’は、本変形例の冷凍機＃４の熱負荷Ｌ４の上限ｂと総消費電力Ｐとの関係を表していると言うことができる。
　ユーザは、図２３に例示するグラフを参照して、冷凍機＃４を稼働させる場合及び稼働しない場合のそれぞれの領域を比較して、解析を行う。
　例えば、図中丸印で囲った領域を参照すると、冷凍機＃４の熱負荷Ｌ４の上限ｂが約１６０以上になると、冷凍機＃４を導入し稼働した方が、稼働しない場合よりも総消費電力Ｐを抑制できることがわかる。
　このように、本変形例によれば、新たに導入する設備の消費電力特性以外の各種の特性についても、総消費エネルギーとの関係を表した論理式を求めることができ、これにより、更にユーザのプラントの解析に資する。
　従来技術によれば、設備の制約条件（上記の例では解析対象の冷凍機＃４の総熱負荷の上限）を変更した場合には、最適な消費エネルギーがどのように変化するかが判然とせず、省エネを実現するためにはどのような制約条件下で設備を運用する必要があるかがわからない、ということがあった。このような場合であっても、本実施形態に係るエネルギー解析支援装置１を用いて解析を行うことで、設備を運用する制約条件を緩和することで省エネ化ができる場合等を検証することが可能となる。
＜第２の実施形態＞
　上記の第１の実施形態に係るエネルギー解析支援装置１は、生成した一階述語論理式から限定記号を消去し、設備の特性と総消費エネルギーとの関係を表す論理式を得てこれからグラフを生成している。これに対し、本実施形態においては、限定記号を消去した結果得られる論理式が、設備の特性同士の間で満たされる関係を表している点で異なる。
　以下、本実施形態に係るエネルギー解析支援装置１について、上記の第１の実施形態と異なる点を中心に説明する。
　本実施形態に係るエネルギー解析支援装置１の構成図は、図１の第１の実施形態に係るそれと同様である。エネルギー解析支援装置１は、計算部１１と出力部１２とを有し、キーボード２Ａやポインティングデバイス２Ｂ等の入力部２やモニタ等の表示部５と接続される。
　エネルギー解析支援装置１は、第１の実施形態と同様に、入力部２を介してシステムモデル、設備モデル及び外部条件の入力を受け付ける。入力されるこれらの情報についても、上記の実施形態と同様であり、図５～図７等を参照して説明したとおりである。
　本実施形態に係るエネルギー解析支援装置１は、入力された情報からシステムモデル数式等の数式を生成して、一階述語論理式を生成し、生成した一階述語論理式から限定記号を消去して論理式を得る。上記のとおり、本実施形態において限定記号を消去して得られる論理式は、それぞれがプラントの設備の特性を表す変数同士が満たす関係を表している。エネルギー解析支援装置１は、得られた論理式に基づきグラフを生成して表示部５に出力し、表示をさせる。
　図２４は、計算部１１における動作フローを示す図である。
　図２の第１の実施形態の動作フローと比較すると、ステップＳ１１については、図２のステップＳ１と同様である。
　ステップＳ１２で、計算部１１は、プラントの運転状態に関する一階述語論理式を生成する。ここで、「プラントの運転状態に関する一階述語論理式」は、図２のステップＳ２で生成する「総消費エネルギーに関する一階述語論理式」と比較すると、存在記号を付与する変数が異なっている。
　ステップＳ１３及びステップＳ１４では、図２のステップＳ３及びステップＳ４と同様の処理を実施する。すなわち、ステップＳ１２で生成した一階述語論理式から限定記号を消去した論理式を計算し、計算により得られたプラント解析数式を出力部１２に出力し、処理を終了する。
　図２５は、本実施形態において計算部１１が生成する一階述語論理式及び一階述語論理式から限定記号を消去して得られる論理式について説明する図である。
　計算部１１は、上記の第１の実施形態と同様に、システムモデル数式Ｆ１、設備モデルＦ２及び外部条件数式Ｆ３ｉを論理積（∧）で結合する（図２５の式（２５－２））。そして、本実施形態においては、計算部１１は、設備の特性を表す２つの変数ａ、ｂ以外の変数に対して存在記号（∃）を付与して（図２５の（２５－１））一階述語論理式ψｉを生成する。
ψi:=∃P∃P1,...,∃PN∃L∃L1,...,∃LN(F1∧F2∧F3i)
　外部条件数式Ｆ３ｉの添え字「ｉ」は、プラント解析のため予め用意されている複数の外部条件数式を互いに識別するために付されている。実施例では、複数の外部条件数式を生成し、その数を「Ｉ」とする。計算部１１は、外部条件数式Ｆ３ｉ（ｉ＝１，２，…，Ｉ）の数（Ｉ）だけ一階述語論理式ψｉを生成する。
　このように、互いに異なる外部条件数式は、それぞれプラントの運転状態を表すことから、本実施形態において生成する一階述語論理式ψｉを、以下においては、「プラントの運転状態に関する一階述語論理式」ともいう。
　計算部１１は、Ｉ個のプラントの運転状態に関する一階述語論理式ψｉのそれぞれに限定記号消去アルゴリズムを適用し、設備の特性を表す変数ａ、ｂが満たす関係を表した論理式φｉを得て、これを出力部１２に出力する。
φi:=QE[ψi]
(φi:=(aとbに関する数式))
　出力部１２は、計算部１１が得たＩ個の論理式φｉ（ｉ＝１，２，…，Ｉ）より、２つの設備の特性を表す変数ａ、ｂが満たす領域を表したグラフを生成する。出力部１２は、これをエネルギー解析支援装置１において解析支援処理を実施した結果得られた情報として表示部５に出力する。表示部５は、出力部１２から入力された情報に基づき、エネルギー解析支援装置１における計算結果を表示する。
　図２６は、出力部１２が生成するグラフを模式的に示す図である。本実施形態においては、一方の軸（横軸）には、設備の特性を表す変数ａを、他方の軸（縦軸）には、設備の特性を表す他の変数ｂをとる。そして、図２６に模式的に示すように、プラントの運転状態ごと、すなわち、複数の外部条件数式Ｆ３ｉのそれぞれについて求めた論理式φｉが満たす領域が、色分け等により可視化される。
　以下、本実施形態に係るエネルギー解析支援装置１が入力された各種モデルを用いてどのように各種数式の計算を行い、計算結果のグラフを得るかについて、第１の実施形態と同様に、具体例を用いつつ、上記第１の実施形態と異なる点を中心に説明する。
　プラントのシステムモデル及び設備モデルについては、第１の実施形態と同様であり、図５及び図６に示すとおりである。
　図２７は、外部条件を例示する図である。
　第１の外部条件は、空調システム全体の総熱負荷Ｌである。設例では、総熱負荷Ｌを５００［ｋＷ］に設定する。
　第２の外部条件は、空調システム全体の総消費電力Ｐである。設例では、第２の外部条件として、４とおりを設定し、それぞれ１５０［ｋＷ］以下、１４０［ｋＷ］以下、１３０［ｋＷ］以下及び１２０［ｋＷ］以下とする。
　図２８は、計算部１１が図５のシステムモデル、図６の設備モデル及び図２７の外部条件から生成する数式群を例示する図である。
　システムモデル数式Ｆ１については、第１の実施形態と同様であり、図８の数式Ｆ１と同様である。
F1:=P=P1+P2+P3+P4∧L=L1+L2+L3+L4
　設備モデル数式Ｆ２については、新たに導入する冷凍機＃４の特性を表す変数のうち、消費電力特性については、第１の実施形態と同様に、「変数ａ」で表す。冷凍機＃４の熱負荷Ｌ４の上限については、「変数ｂ」で表す。
F2:=(P1=0.12*L1+16∧69<=L1<=125∨P1=0∧L1=0)
∧(P2=0.17*L2+21∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212∨P3=0∧L3=0)
∧(P4=a*L4+15∧69<=L4<=b∧0.1<=a<=0.35∧69<=b<=212∨P4=0∧L4=0)
　外部条件数式Ｆ３ｉ（ｉ＝１，２，３，４）については、第１の外部条件及び第２の外部条件がいずれも成立する場合の数式を生成する。すなわち、第１の外部条件の式と第２の外部条件の式とを論理積（∧）で結合して、外部条件数式Ｆ３ｉとする。
F31:=L=500∧P<=150
F32:=L=500∧P<=140
F33:=L=500∧P<=130
F34:=L=500∧P<=120
　図２９は、図２８の数式群Ｆ１、Ｆ２、Ｆ３ｉから一階述語論理式ψｉを生成する方法を説明する図である。
　先に図２５の説明においても述べたとおり、計算部１１は、生成した図２８の数式Ｆ１、Ｆ２、Ｆ３ｉを論理積（∧）で結合し、変数ａ、ｂ以外の変数に存在記号（∃）を付与して、一階述語論理式ψｉを生成する。図２９においては、４とおりの外部条件のうち、ｉ＝１の外部条件、具体的には、空調システムの総熱負荷Ｌが５００［ｋＷ］で、空調システムの総消費電力Ｐが１５０［ｋＷ］以下であることを外部条件とする場合の一階述語論理式ψ１を例示する。他の３とおりの外部条件の場合の一階述語論理式ψ２～ψ４の具体的な数式については、外部条件数式Ｆ３ｉ以外は、ψ１と共通であるため、図２９においては記載を省略している。すなわち、図２９のψ１の式の外部条件数式Ｆ３１を上記のＦ３２～Ｆ３４に置き換えた式が、一階述語論理式ψ２～ψ４である。
ψ1:=∃P∃P1∃P2∃P3∃P4∃L∃L1∃L2∃L3∃L4(F1∧F2∧F31)
ψ2:=∃P∃P1∃P2∃P3∃P4∃L∃L1∃L2∃L3∃L4(F1∧F2∧F32)
ψ3:=∃P∃P1∃P2∃P3∃P4∃L∃L1∃L2∃L3∃L4(F1∧F2∧F33)
ψ4:=∃P∃P1∃P2∃P3∃P4∃L∃L1∃L2∃L3∃L4(F1∧F2∧F34)
F1∧F2∧F31:=
P=P1+P2+P3+P4∧L=L1+L2+L3+L4
∧(P1=0.12*L1+16∧69<=L1<=125∨P1=0∧L1=0)
∧(P2=0.17*L2+21∧95<=L2<=175∨P2=0∧L2=0)
∧(P3=0.24*L3+26∧122<=L3<=212∨P3=0∧L3=0)
∧(P4=a*L4+15∧69<=L4<=b∧0.1<=a<=0.35∧69<=b<= 212∨P4=0∧L4=0)
∧L=500∧P<=150)
　図３０は、一階述語論理式ψ１～ψ４から限定記号を消去して得られる論理式φ１～φ４を示す図である。
　計算部１１は、生成した一階述語論理式ψ１～ψ４のそれぞれに限定記号消去アルゴリズム（ＱＥ［］）を適用して、存在記号の付与された変数を消去した論理式φ１～φ４を得る。論理式φ１～φ４は、冷凍機＃４の消費電力特性ａと、熱負荷Ｌ４の上限ｂとが満たす関係式からなる。
φ1:=QE[ψ1]
(φ1:=(-213/800+a=0∧-212+b<=0∧200-b<=0)∨(-678/4075+a=0∧-212+b<=0∧163-b<=0)∨(-2897/20300+a=0∧-212+b<=0∧203-b<=0)∨(-2097/19300+a=0∧-212+b<=0∧193-b<=0)∨(-3/25+a<0∧1/10-a<=0∧-212+b<=0∧200-b<=0) ...以下省略)
φ2:=QE[ψ2]
(φ2:=(-173/800+a=0∧-212+b<=0∧200-b<=0)∨(-428/4075+a=0∧-212+b<=0∧163-b<=0)∨(-3/25+a<0∧1/10-a<=0∧-212+b<=0∧200-b<=0)∨(-173/800+a<=0∧1/10-a<=0∧-212+b<=0∧200-b<=0)∨(-428/4075+a<=0∧1/10-a<=0∧-212+b<=0∧163-b<=0)...以下省略)
φ3:=QE[ψ3]
(φ3:=(-133/800+a=0∧-212+b<=0∧200-b<=0)∨(-3/25+a<0∧1/10-a<=0∧-212+b<=0∧200-b<=0)∨(-133/800+a<=0∧1/10-a<=0∧-212+b<=0∧200-b<=0)∨(-133/800+a<=0∧3/25-a<0∧-212+b<=0∧200-b<=0)...以下省略)
φ4:=QE[ψ4]
(φ4:=(-93/800+a=0∧-212+b<=0∧200-b<=0)∨(-93/800+a<=0∧1/10-a<=0∧-212+b<=0∧200-b<=0)∨(-7/20+a<=0∧1/10-a<=0∧-212+b<=0∧163-b<0∧1050-6*b+25*a*b=0)∨(-7/20+a<=0∧1/10-a<=0∧-212+b<=0∧163-b<0∧1050-6*b+25*a*b<=0)...以下省略)
　図３１は、表示部５に表示する画面の一例を示す図であり、各論理式φ１～φ４が満たす領域Ａ１～Ａ４をそれぞれグラフで示している。縦軸には、冷凍機＃４の熱負荷Ｌ４の上限ｂを、横軸には冷凍機＃４の消費電力特性ａをとる。図３１の領域Ａ１～Ａ４は、それぞれ空調システムの総消費電力Ｐを１５０、１４０、１３０、１２０［ｋＷ］以下とするために必要な冷凍機＃４の設備能力を表している。
　図３１に示すように、外部条件により、冷凍機＃４の消費電力特性ａとその熱負荷Ｌ４の上限ｂとが満たす領域Ａ１～Ａ４は、範囲が異なる。ユーザは、例えば図３１の４つのグラフを比較しながら、冷凍機＃４の消費電力特性ａと熱負荷Ｌ４の上限ｂという互いに相関のある特性に応じて総消費電力Ｐがどのように変わるかを検討することができる。
　領域Ａ１～Ａ４は、図３１に例示するようにそれぞれを画面上に並べて配置させてもよいし、１つのグラフに各領域のどのあたりが重複しているかを視認可能に表示させてもよい。
　図３２は、表示部５に表示する画面の他の例を示す図であり、各論理式φ１～φ４が満たす領域Ａ１～Ａ４を色分け等により識別可能に表示している。図３２の例では、斜線や網掛け等により各領域Ａ１～Ａ４をユーザが識別可能に表示している。
　ユーザは、図３２に例示するグラフを参照して、例えば、新たに導入する冷凍機＃４の消費電力特性ａが０．１１のとき、総消費電力Ｐを所定値以下に抑制しようとする場合の熱負荷Ｌ４の上限値を確認する。総消費電力を１５０［ｋＷ］以下にするには、冷凍機＃４の熱負荷の上限を１５２［ｋＷ］以上にする必要があることがわかる。同様に、総消費電力を１４０［ｋＷ］以下にするには、当該上限を１６８［ｋＷ］以上に、総消費電力を１３０［ｋＷ］以下にするには、当該上限を２００［ｋＷ］以上にする必要があることがわかる。
　また、冷凍機＃４の熱負荷の上限ｂが２００［ｋＷ］以上の範囲で総消費電力を所定値以下に抑制しようとする場合の冷凍機＃４の消費電力特性ａの最適な値を決定することができる。熱負荷の上限が２００［ｋＷ］以上の条件の下で、総消費電力を１５０［ｋＷ］以下にするには、消費電力特性ａを０．２７以下にすればよいことがわかる。同様に、同条件の下で総消費電力を１４０［ｋＷ］以下にするには、消費電力特性ａを０．１８以下にすればよいことがわかる。
　このように、本実施形態に係るエネルギー解析支援装置１によれば、相関のある複数の設備の特性間の関係を表した式が得られ、これにより、プラントの総消費電力がどのように変わるかを可視化できる。ユーザにとっては、可視化した結果を用いて、例えばどのような効率及び制約条件を備えた設備を導入すれば省エネ効果がどの程度見込めるのかの検討が容易になる。
　本実施形態では、新たに導入する設備の制約条件（上記の例では熱負荷の上限）を一方の変数ｂに設定して、これを含む論理式φｉを得てこれを解析用にユーザに提供している。従来は、ある設備の制約条件を変えた場合には、最適な消費エネルギーがどのように変化するかがわからなかったため、省エネを実現するには、設備をどのような制約条件下で運用する必要があるかもわからなかった。しかし、本実施形態に係るエネルギー解析支援装置１によれば、例えば制約条件を緩和した（上記の例では熱負荷の上限を上げた）場合には省エネを実現できることがあるが、そのような状態についても検証することが可能となる。
　ところで、図１のエネルギー解析支援装置１は、例えば、図３３に示すような情報処理装置（コンピュータ）を用いて構成することができる。図３３の情報処理装置は、ＣＰＵ（中央処理装置）１００１、メモリ１００２、入力装置１００３、出力装置１００４、外部記憶装置１００５、媒体駆動装置１００６、ネットワーク接続装置１００７を備え、それらはバス１００８により互いに接続されている。
　メモリ１００２は、例えば、ＲＯＭ（read only memory）、ＲＡＭ（random access memory）等を含み、処理に用いられるプログラムおよびデータを格納する。ＣＰＵ１００１は、メモリ１００２を利用してプログラムを実行することにより、必要な処理を行う。
　図１の計算部１１及び出力部１２は、メモリ１００２に格納されたプログラムを実行することにより実現される機能に対応する。
　入力装置１００３は、例えば、キーボード、ポインティングデバイス、タッチパネル等であり、ネットワーク管理者等のオペレータからの指示や情報の入力に用いられる。出力装置１００４は、例えば、ディスプレイ、プリンタ、スピーカ等であり、ユーザへの問い合わせ、アラーム、処理結果等の出力に用いられる。
　外部記憶装置１００５は、例えば、磁気ディスク装置、光ディスク装置、光磁気ディスク装置、テープ装置等である。情報処理装置は、この外部記憶装置１００５に、上記プログラムおよびデータを格納しておき、必要に応じて、それらをメモリ１００２にロードして使用する。
　媒体駆動装置１００６は、可搬記録媒体１００９を駆動し、その記録内容にアクセスする。可搬記録媒体１００９は、メモリカード、フレキシブルディスク、ＣＤ－ＲＯＭ（compact disk read only memory ）、光ディスク、光磁気ディスク等の任意のコンピュータ読み取り可能な記録媒体である。オペレータは、この可搬記録媒体１００９に上記プログラムおよびデータを格納しておき、必要に応じて、それらをメモリ１００２にロードして使用する。
　ネットワーク接続装置１００７は、ＬＡＮ（local area network）、インターネット等の任意の通信ネットワークに接続され、通信に伴うデータ変換を行う。情報処理装置は、必要に応じて、上記プログラムおよびデータを外部の装置からネットワーク接続装置１００７を介して受け取り、それらをメモリ１００２にロードして使用する。
　図３４は、図３３の情報処理装置にプログラムおよびデータを供給することのできるコンピュータ読み取り可能な記録媒体を示している。可搬記録媒体１００９やサーバ１１０１のデータベース１１０３に格納されたプログラムおよびデータは、情報処理装置１１０２のメモリ１００２にロードされる。サーバ１１０１は、そのプログラムおよびデータを搬送する搬送信号を生成し、ネットワーク上の任意の伝送媒体を介して情報処理装置１１０２に送信する。ＣＰＵ１００１は、そのデータを用いてそのプログラムを実行し、必要な処理を行う。
　本発明は上述した実施形態そのままに限定されるものではく、実施段階でのその要旨を逸脱しない範囲で構成要素を変形して具体化することができる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成することができる。例えば、実施形態に示される全構成要素を適宜組み合わせても良い。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。このような、発明の趣旨を逸脱しない範囲内において種々の変形や応用が可能であることは言うまでもない。
　１　　　エネルギー解析支援装置
　２　　　入力部
　５　　　表示部
　１１　　計算部
　１２　　出力部

Claims (17)

  1. 　複数の設備からなるプラントのエネルギー解析支援装置であって、
    　前記プラントのモデルと、前記プラントを構成する設備の能力を表す設備モデルと、前記プラントの消費エネルギーまたは負荷を表す外部条件とが与えられると、前記プラントのモデルを表す数式、前記設備モデルを表す数式及び外部条件を表す数式より、解析対象の設備の特性を表した少なくとも一の変数を含む論理式を生成する計算部と、
    　前記計算部により生成された論理式に基づいて、前記一の変数が表す前記解析対象の設備の特性を２軸のうちの一の軸とするグラフを生成して出力する出力部と、
    　を有することを特徴とするエネルギー解析支援装置。
  2. 　前記計算部は、前記プラントのモデル、前記設備モデル及び前記外部条件より、一階述語論理式を生成し、生成した一階述語論理式に限定記号消去アルゴリズムを用いて限定記号を消去することにより、前記一の変数を含んだ論理式を算出する
    　ことを特徴とする請求項１記載のエネルギー解析支援装置。
  3. 　前記計算部は、
    　　前記プラントのモデル、前記設備モデル及び前記外部条件より、前記プラント全体の総負荷と前記プラントの各設備に配分される負荷との関係、及び前記プラントにおける各設備の消費エネルギーと前記プラント全体の総消費エネルギーとの関係を表すシステムモデル数式と、各設備の負荷と消費エネルギーの関係を表す設備モデル数式と、前記プラント全体の総負荷または前記プラント全体の消費エネルギーを表した外部条件数式とを生成し、
    　　前記システムモデル数式、前記設備モデル数式及び前記外部条件数式を論理積で結合し、前記限定記号消去アルゴリズムにて消去すべき変数に対して存在記号を付与する
    　ことにより、前記一階述語論理式を生成することを特徴とする請求項２記載のエネルギー解析支援装置。
  4. 　前記計算部は、前記システムモデル数式、前記設備モデル数式及び前記外部条件数式を論理積で結合した式に、前記一の変数として前記生成する論理式に含まれる前記解析対象の設備の特性と、前記プラントの総消費エネルギーを表す変数以外の全ての変数に限定記号を付与することにより、前記一階述語論理式を生成することを特徴とする請求項３記載のエネルギー解析支援装置。
  5. 　前記計算部は、前記解析対象の設備を稼働させた場合及び停止させた場合のそれぞれについて前記設備モデル数式を生成して、それぞれの場合についての前記解析対象の設備の特性と前記プラントの総消費エネルギーとの実行可能な関係を表した論理式を得る
    　ことを特徴とする請求項４記載のエネルギー解析支援装置。
  6. 　前記計算部において限定記号を消去して得られる前記論理式に含まれる前記一の変数は、前記解析対象の設備の消費エネルギー効率を表す
    　ことを特徴とする請求項４記載のエネルギー解析支援装置。
  7. 　前記計算部において限定記号を消去して得られる前記論理式に含まれる前記一の変数は、前記解析対象の設備の負荷についての制約条件を表す
    　ことを特徴とする請求項４記載のエネルギー解析支援装置。
  8. 　前記計算部は、前記システムモデル数式、前記設備モデル数式及び前記外部条件数式を論理積で結合した式に、前記解析対象の設備の特性を表す前記一の変数及び前記解析対象の設備の他の特性を表す他の変数以外の全ての変数に限定記号を付与することにより、前記一階述語論理式を生成することを特徴とする請求項３記載のエネルギー解析支援装置。
  9. 　前記計算部において限定記号を消去して得られる前記論理式に含まれる前記一の変数は、前記解析対象の設備の消費エネルギー効率を表し、前記他の変数は、前記解析対象の設備の負荷についての制約条件を表す
    　ことを特徴とする請求項８記載のエネルギー解析支援装置。
  10. 　前記出力部は、前記限定記号を消去して得られた論理式に含まれる変数のうち、前記プラントの総消費エネルギーを表す変数または前記解析対象の設備の他の特性を表す他の変数を前記２軸のうちの他の軸として、前記グラフを生成して出力する
    　ことを特徴とする請求項４から９のいずれか１項に記載のエネルギー解析支援装置。
  11. 　前記計算部は、前記外部条件数式として、前記プラントの総消費エネルギーが所定の範囲となることを表す式を生成する
    　ことを特徴とする請求項１から１０のいずれか１項に記載のエネルギー解析支援装置。
  12. 　前記計算部は、前記外部条件数式として、前記プラントの総負荷が所定の範囲となることを表す式を生成する
    　ことを特徴とする請求項１から１０のいずれか１項に記載のエネルギー解析支援装置。
  13. 　前記計算部は、複数とおりの前記外部条件が与えられると、前記複数とおりの外部条件に対応する前記複数とおりの論理式を算出し、
    　前記出力部は、前記複数とおりの外部条件のそれぞれについて求めた論理式に基づいて、前記グラフを生成して出力する
    　ことを特徴とする請求項１から１２のいずれか１項に記載のエネルギー解析支援装置。
  14. 　前記出力部は、前記複数とおりの外部条件のそれぞれについて求めた論理式が満たす領域を同時に表示させたグラフを生成して出力する
    　ことを特徴とする請求項１３記載のエネルギー解析支援装置。
  15. 　複数の設備からなるプラントのエネルギー解析支援処理を情報処理装置に実行させるためのエネルギー解析支援プログラムであって、
    　前記プラントのモデルと、前記プラントを構成する設備の能力を表す設備モデルと、前記プラントの消費エネルギーまたは負荷を表す外部条件とが与えられると、前記プラントのモデルを表す数式、前記設備モデルを表す数式及び外部条件を表す数式より、解析対象の設備の特性を表した少なくとも一の変数を含む論理式を生成し、
    　前記生成された論理式に基づいて、前記一の変数が表す前記解析対象の設備の特性を２軸のうちの一の軸とするグラフを生成して出力する、
    　ことを特徴とするエネルギー解析支援プログラム。
  16. 　複数の設備からなるプラントのエネルギー解析支援処理を情報処理装置に実行させるためのエネルギー解析支援プログラムを記録した記録媒体であって、
    　前記プラントのモデルと、前記プラントを構成する設備の能力を表す設備モデルと、前記プラントの消費エネルギーまたは負荷を表す外部条件とが与えられると、前記プラントのモデルを表す数式、前記設備モデルを表す数式及び外部条件を表す数式より、解析対象の設備の特性を表した少なくとも一の変数を含む論理式を生成し、
    　前記生成された論理式に基づいて、前記一の変数が表す前記解析対象の設備の特性を２軸のうちの一の軸とするグラフを生成して出力する、
    　ことを特徴とする記録媒体。
  17. 　複数の設備からなるプラントのエネルギー解析支援方法であって、
    　前記プラントのモデルと、前記プラントを構成する設備の能力を表す設備モデルと、前記プラントの消費エネルギーまたは負荷を表す外部条件とが与えられると、前記プラントのモデルを表す数式、前記設備モデルを表す数式及び外部条件を表す数式より、解析対象の設備の特性を表した少なくとも一の変数を含む論理式を生成し、
    　前記生成された論理式に基づいて、前記一の変数が表す前記解析対象の設備の特性を２軸のうちの一の軸とするグラフを生成して出力する、
    　ことを特徴とするエネルギー解析支援方法。
PCT/JP2015/078483 2015-10-07 2015-10-07 エネルギー解析支援装置、エネルギー解析支援方法、エネルギー解析支援プログラム及びこれを記録した記録媒体 Ceased WO2017060985A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/JP2015/078483 WO2017060985A1 (ja) 2015-10-07 2015-10-07 エネルギー解析支援装置、エネルギー解析支援方法、エネルギー解析支援プログラム及びこれを記録した記録媒体
JP2017544112A JP6477906B2 (ja) 2015-10-07 2015-10-07 エネルギー解析支援装置、エネルギー解析支援方法、エネルギー解析支援プログラム及びこれを記録した記録媒体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2015/078483 WO2017060985A1 (ja) 2015-10-07 2015-10-07 エネルギー解析支援装置、エネルギー解析支援方法、エネルギー解析支援プログラム及びこれを記録した記録媒体

Publications (1)

Publication Number Publication Date
WO2017060985A1 true WO2017060985A1 (ja) 2017-04-13

Family

ID=58488266

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/078483 Ceased WO2017060985A1 (ja) 2015-10-07 2015-10-07 エネルギー解析支援装置、エネルギー解析支援方法、エネルギー解析支援プログラム及びこれを記録した記録媒体

Country Status (2)

Country Link
JP (1) JP6477906B2 (ja)
WO (1) WO2017060985A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017169354A (ja) * 2016-03-16 2017-09-21 富士電機株式会社 プラント運用支援装置、プラント運用支援方法、プログラム及びプラント運用支援システム

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012037141A (ja) * 2010-08-06 2012-02-23 Mitsubishi Heavy Ind Ltd 冷凍機制御装置
JP2015102263A (ja) * 2013-11-22 2015-06-04 富士電機株式会社 空調制御支援装置、空調制御支援プログラムおよび空調制御支援方法
WO2015115312A1 (ja) * 2014-01-30 2015-08-06 富士電機株式会社 エネルギー解析装置、エネルギー解析プログラム及び記録媒体

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012037141A (ja) * 2010-08-06 2012-02-23 Mitsubishi Heavy Ind Ltd 冷凍機制御装置
JP2015102263A (ja) * 2013-11-22 2015-06-04 富士電機株式会社 空調制御支援装置、空調制御支援プログラムおよび空調制御支援方法
WO2015115312A1 (ja) * 2014-01-30 2015-08-06 富士電機株式会社 エネルギー解析装置、エネルギー解析プログラム及び記録媒体

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YOSHIO TANGE: "Visualization of Energy Optimization Problems under Supply and Demand Balance for Energy Saving Pre-verification", THE TRANSACTIONS OF THE INSTITUTE OF ELECTRICAL ENGINEERS OF JAPAN C, vol. 134, no. 1, 1 January 2014 (2014-01-01), pages 78 - 84, XP008184317 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017169354A (ja) * 2016-03-16 2017-09-21 富士電機株式会社 プラント運用支援装置、プラント運用支援方法、プログラム及びプラント運用支援システム

Also Published As

Publication number Publication date
JP6477906B2 (ja) 2019-03-06
JPWO2017060985A1 (ja) 2018-01-25

Similar Documents

Publication Publication Date Title
US11953865B2 (en) HVAC controller with predictive cost optimization
US9915178B2 (en) Online heurisitc algorithm for combined cooling heating and power plant optimization
Deng et al. A competitive memetic algorithm for carbon-efficient scheduling of distributed flow-shop
Keshavarz et al. Makespan minimisation in flexible flowshop sequence-dependent group scheduling problem
Lei et al. Hybrid heuristic search approach for deadlock-free scheduling of flexible manufacturing systems using Petri nets
Inkulu et al. Human-robot collaborative task planning for assembly system productivity enhancement
Li et al. Scheduling FMS problems with heuristic search function and transition-timed Petri nets
Lin et al. Minimising makespan in job-shops with deterministic machine availability constraints
Batur et al. Quantile regression metamodeling: Toward improved responsiveness in the high-tech electronics manufacturing industry
Alves et al. Planning on discrete event systems using parallelism maximization
Li et al. A fast algorithm for buffer allocation problem
Xu et al. An ordinal transformation framework for multi-fidelity simulation optimization
Cheng et al. Two-agent single-machine scheduling to minimize the weighted sum of the agents’ objective functions
Wang et al. A complexity model for assembly supply chains and its application to configuration design
Li et al. Quantitative analysis of a transfer production line with Andon
JP6477906B2 (ja) エネルギー解析支援装置、エネルギー解析支援方法、エネルギー解析支援プログラム及びこれを記録した記録媒体
JP6119574B2 (ja) 空調制御支援装置、空調制御支援プログラムおよび空調制御支援方法
Yuan et al. Digital-twin-based dynamic flexible job shop scheduling problem via multi-agent proximal policy optimisation
Dai et al. A GA-based heuristic approach for offshore structure construction spatial scheduling under uncertainty
JP6458561B2 (ja) 負荷配分決定支援装置、そのプログラム、負荷配分決定支援方法
Lin et al. Relocation scheduling subject to fixed processing sequences
JP6142959B1 (ja) プラント運用支援装置、表示装置、プラント運用システム、及びプラント運用支援方法
JP6477293B2 (ja) 系統運用支援装置、系統運用支援方法及びそのプログラム
Shi et al. A novel two stage neighborhood search for flexible job shop scheduling problem considering reconfigurable machine tools
WO2018143170A1 (ja) 設計プロセス支援システム、制御方法及びプログラム

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15905805

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2017544112

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15905805

Country of ref document: EP

Kind code of ref document: A1