WO2021002464A1 - 組合せ解決定システム - Google Patents

組合せ解決定システム Download PDF

Info

Publication number
WO2021002464A1
WO2021002464A1 PCT/JP2020/026247 JP2020026247W WO2021002464A1 WO 2021002464 A1 WO2021002464 A1 WO 2021002464A1 JP 2020026247 W JP2020026247 W JP 2020026247W WO 2021002464 A1 WO2021002464 A1 WO 2021002464A1
Authority
WO
WIPO (PCT)
Prior art keywords
series data
solution
combination solution
unit
essential
Prior art date
Application number
PCT/JP2020/026247
Other languages
English (en)
French (fr)
Inventor
裕太 泉
西村　忠史
森田　浩
晴貴 井上
Original Assignee
ダイキン工業株式会社
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 ダイキン工業株式会社 filed Critical ダイキン工業株式会社
Priority to CN202080048321.2A priority Critical patent/CN114041031B/zh
Priority to EP20834972.0A priority patent/EP3996011A4/en
Priority to US17/623,167 priority patent/US11568108B2/en
Publication of WO2021002464A1 publication Critical patent/WO2021002464A1/ja

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24—HEATING; RANGES; VENTILATING
    • F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00—Control or safety arrangements
    • F24F11/30—Control or safety arrangements for purposes related to the operation of the system, e.g. for safety or monitoring
    • F24F11/46—Improving electric energy efficiency or saving
    • F24F11/47—Responding to energy costs
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N5/00—Computing arrangements using knowledge-based models
    • G06N5/01—Dynamic search techniques; Heuristics; Dynamic trees; Branch-and-bound
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N99/00—Subject matter not provided for in other groups of this subclass
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Definitions

  • Patent Document 1 Japanese Unexamined Patent Publication No. 2006-48475
  • Patent Document 1 Japanese Unexamined Patent Publication No. 2006-48475
  • the mathematical programming problem is a black box optimization problem in which the objective function cannot be formulated, there is a problem that a lot of calculation time is required to derive the optimum solution.
  • the combination solution determination system of the first viewpoint includes a solution candidate generation unit, a simulation unit, an evaluation index calculation unit, a solution determination unit, and an essential series data extraction unit.
  • the solution candidate generation unit generates candidates for combination solutions.
  • the simulation unit calculates simulation data using the information about the combination solution candidates generated by the solution candidate generation unit and the series data for evaluating the combination solution.
  • the evaluation index calculation unit calculates the evaluation index based on the simulation data.
  • the solution determination unit determines a highly evaluated combination solution from a plurality of combination solution candidates based on the evaluation index calculated by the evaluation index calculation unit from each of the plurality of combination solution candidates.
  • the essential series data extraction unit contains information on m (m ⁇ 2) combination solution candidates among the n (n ⁇ 2) combination solution candidates generated by the solution candidate generation unit, and the first series data. Based on the 1-series data and the evaluation index calculated from the evaluation index calculation unit, the second-series data necessary for calculating the evaluation index with a predetermined accuracy is extracted from the first-series data.
  • the essential series data extraction unit extracts the third series data necessary for verifying the suitability of a predetermined constraint condition for the m combination solution candidates from the first series data.
  • the essential series data extraction unit acquires the essential series data by combining the second series data and the third series data. After the essential series data extraction unit extracts the essential series data, the simulation unit must provide at least information on combination solution candidates other than the m combination solution candidates among the n combination solution candidates, and essential information.
  • the simulation data is calculated using the series data including the series data.
  • the combination solution determination system of the first viewpoint can efficiently derive the solution of the black box optimization problem.
  • the combination solution determination system of the second viewpoint is the combination solution determination system of the first viewpoint, and the essential series data extraction unit converts the first series data, which is likely not to satisfy a predetermined constraint condition, into the third series data. Extract as.
  • the combination solution determination system of the second viewpoint can efficiently derive the solution of the black box optimization problem.
  • the combination solution determination system of the third viewpoint is a combination solution determination system of the first viewpoint or the second viewpoint, and the essential series data extraction unit extracts the second series data by sparse estimation and the third by the extreme value statistics. Extract series data.
  • the combination solution determination system of the third viewpoint can efficiently derive the solution of the black box optimization problem.
  • the combination solution determination system of the fourth viewpoint is a combination solution determination system of any one of the first to third viewpoints, and the essential series data extraction unit is unnecessary data for calculating the evaluation index with a predetermined accuracy.
  • the essential series data is acquired by excluding the data that is unnecessary when considering the predetermined constraint condition from the first series data.
  • the combination solution determination system of the fourth viewpoint can efficiently derive the solution of the black box optimization problem.
  • the combination solution determination system of the fifth viewpoint is a combination solution determination system of any one of the first to fourth viewpoints, and the combination solution is a combination solution of objects or a combination solution of methods.
  • the combination solution determination system of the fifth viewpoint can efficiently derive the solution of the black box optimization problem.
  • the combination solution determination system of the sixth viewpoint is a combination solution determination system of any one of the first to fifth viewpoints, and the combination solution improves the situation of the application target of the combination solution, and is the first series.
  • the data is historical data about the situation that affects the improvement of the situation to which the combination solution is applied.
  • the combination solution determination system of the sixth viewpoint can efficiently derive the solution of the black box optimization problem.
  • the combination solution determination system determines the solution of the mathematical programming problem that derives the air conditioning system that consumes the least amount of power.
  • FIG. 1 is a diagram showing an example of the configuration of the air conditioning system 100.
  • the air conditioning system 100 is composed of two outdoor units 10a and 10b, six indoor units 20a to 20f, and two ventilation devices 30a and 30b.
  • the outdoor units 10a and 10b are connected to one or a plurality of indoor units 20a to 20f via a refrigerant pipe.
  • the indoor units 20a to 20f are connected to any one of the two outdoor units 10a and 10b, and are installed in any one of the zones 40a to 40c of the building in which the air conditioner system 100 is installed. There is.
  • FIG. 1 is a diagram showing an example of the configuration of the air conditioning system 100.
  • the air conditioning system 100 is composed of two outdoor units 10a and 10b, six indoor units 20a to 20f, and two ventilation devices 30a and 30b.
  • the outdoor units 10a and 10b are connected to one or a plurality of indoor units 20a to 20f via a refrigerant pipe.
  • the outdoor units 10a and 10b are connected to three indoor units 20a to 20f, and two indoor units 20a to 20f are installed in each zone 40a to 40c.
  • Zones 40a to 40c are spaces subject to air conditioning by the air conditioning system 100.
  • the indoor units 20a to 20f maintain the zones 40a to 40c in a comfortable state by removing the sensible heat of the zones 40a to 40c.
  • the ventilators 30a, 30b ventilate one or more zones 40a-40c to keep the zones 40a-40c in a comfortable state.
  • the parameters (air conditioning selection parameters) determined by the selection of the air conditioning system 100 are, for example, the number of indoor units 20a to 20f in each zone 40a to 40c, the model and performance (capacity, etc.) of the indoor units 20a to 20f, and the outdoor unit 10a.
  • the refrigerant system is, for example, information on a set of outdoor units 10a and 10b and indoor units 20a to 20f connected to the outdoor units 10a and 10b.
  • the control parameters are, for example, the set temperature and the set humidity of each zone 40a to 40c.
  • the combination solution determination system solves the above problems, and efficiently selects the air conditioning system that consumes the least power from a large number of options. As a result, the combination solution determination system can select an air conditioning system that minimizes the total cost.
  • the total cost is the sum of the equipment cost and the electricity cost.
  • the equipment cost includes, for example, the cost of the equipment itself, the installation construction cost of the equipment, and the maintenance cost of the equipment.
  • the electricity cost includes the electricity charge required to operate the air conditioning system and is calculated from the power consumption of the equipment.
  • the mathematical programming problem whose solution is determined by the combination solution determination system according to the embodiment corresponds to a black box optimization problem in which a general-purpose solver cannot be used because the objective function cannot be formulated.
  • the combination solution determination system according to the embodiment can reduce the calculation time by using the iterative local search method, which is an application method of the local search method that can also be used for the black box optimization problem.
  • FIG. 2 is a schematic block diagram of the combination solution determination system 190.
  • the combination solution determination system 190 selects an air conditioning system that consumes the least amount of power.
  • the combination solution determination system 190 is composed of, for example, one or more computers. When the combination solution determination system 190 is composed of a plurality of computers, the plurality of computers may be connected to each other via a network.
  • the combination solution determination system 190 mainly includes a solution candidate generation unit 110, a simulation unit 120, an evaluation index calculation unit 130, a solution determination unit 140, and an essential series data extraction unit 150.
  • the solution candidate generation unit 110 to the essential sequence data extraction unit 150 are realized, for example, by the CPU of the computer executing a program stored in the storage device of the computer constituting the combination solution determination system 190.
  • the solution candidate generation unit 110 generates candidates for a combination solution.
  • the combination solution is a solution determined by the combination solution determination system 190.
  • the combination solution is information on the configuration of the air conditioning system 100, and specifically includes the above-mentioned air conditioning selection parameters.
  • the simulation unit 120 calculates the simulation data using the information about the candidate of the combination solution generated by the solution candidate generation unit 110 and the series data for evaluating the combination solution.
  • the simulation data includes the power consumption of the air conditioning system 100, which is calculated by simulation using the candidate of the combination solution as an input parameter.
  • the series data is data related to a predetermined time when simulation data (power consumption) is calculated.
  • the series data is, for example, the time of each hour for the one year. In this case, the series data is composed of 8760 time data (365 days ⁇ 24 hours / day).
  • the power consumption calculated by the simulation unit 120 is the power consumption of the outdoor units 10a and 10b, the indoor units 20a to 20f, and the ventilation devices 30a and 30b.
  • the power consumption of the ventilation devices 30a and 30b is calculated by a ventilation simulation.
  • the power consumption of the outdoor units 10a and 10b and the indoor units 20a to 20f is calculated by air conditioning simulation from the heat load (sensible heat load) of each zone 40a to 40c in consideration of the influence of the heat load by the ventilation devices 30a and 30b.
  • the evaluation index calculation unit 130 calculates the evaluation index based on the simulation data (power consumption) calculated by the simulation unit 120.
  • the evaluation index is a parameter calculated from the above-mentioned total cost, untreated heat load, untreated ventilation volume, and the like.
  • the untreated heat load is a sensible heat load that cannot be treated when the sensible heat load that can be processed by the indoor units 20a to 20f arranged in the zones 40a to 40c is less than the sensible heat load applied in the zones 40a to 40c.
  • the untreated ventilation volume corresponds to the ventilation load that cannot be treated when the ventilation load that can be processed by the ventilation devices 30a and 30b arranged in the zones 40a to 40c is less than the ventilation load applied in the zones 40a to 40c. ..
  • the smaller the total cost, untreated heat load, and untreated ventilation volume the smaller the evaluation index. From the viewpoint of maximizing the use of resources, the smaller the evaluation index is, the more preferable.
  • the solution determination unit 140 determines a highly evaluated combination solution from a plurality of combination solution candidates based on the evaluation index calculated by the evaluation index calculation unit 130 from each of the plurality of combination solution candidates.
  • the highly evaluated combination solution is the combination solution having the smallest evaluation index.
  • the combination solution determined by the solution determination unit 140 is the optimum air conditioning selection parameter from the viewpoint of maximizing the use of resources.
  • the required series data extraction unit 150 acquires the required series data.
  • the required series data is data extracted from the series data. More specifically, as described above, when the series data is composed of 8760 time data per hour, the essential series data is referred to as the 8760 time data (hereinafter referred to as "first series data"). .) Is a part of the time data.
  • the simulation unit 120 calculates the simulation data using the information on the candidates for the combination solution and the essential series data. Therefore, the simulation unit 120 can calculate the simulation data (power consumption) based on the time data that is smaller than the time data included in the first series data. As a result, the combination solution determination system 190 can reduce the calculation cost required for calculating the power consumption.
  • FIG. 3 is a flowchart of a process for selecting an air conditioning system by the combination solution determination system 190.
  • step S11 the system configuration is determined. Specifically, in step S11, the solution candidate generation unit 110 generates a combination solution candidate.
  • step S12 it is determined whether or not the required series data has been extracted. If the required series data has been extracted in step S16, the process proceeds to step S13, and if it has not been extracted, the process proceeds to step S18.
  • step S13 the simulation is executed. Specifically, the simulation unit 120 calculates the simulation data at each time included in the first series data by using the combination solution candidates and the first series data.
  • step S14 the simulation data (power consumption) at each time obtained in step S13 is saved as training data.
  • step S15 the evaluation index calculation unit 130 calculates and saves the evaluation index from the training data saved in step S14.
  • step S16 it is determined whether or not the first local solution has been derived based on the evaluation index calculated from the training data.
  • the first local solution is calculated by iterative local search.
  • the first local solution is, for example, a system configuration in which the evaluation index first shows the minimum value (local minimum value) from the start of simulation execution.
  • the essential series data is extracted from the first series data based on the training data obtained so far in step S17. Specifically, the essential series data extraction unit 150 extracts the time required for calculating the power consumption. Then, in step S11, a new system configuration is determined and a simulation is executed.
  • step S16 If the first local solution is not derived in step S16, the required series data is not extracted in step S17, a new system configuration is determined in step S11, and further training data is acquired.
  • step S18 the simulation is executed based on the time (essential series data) extracted in step S17. Specifically, the simulation unit 120 calculates the simulation data at each time included in the essential series data by using the combination solution candidates and the essential series data.
  • step S19 the evaluation index calculation unit 130 calculates and saves the evaluation index from the simulation data obtained in step S18.
  • step S11 the system configuration is newly determined in step S11, the simulation is executed in step S18, and the evaluation index is calculated and saved in step S19.
  • the above processing is performed on a predetermined number of combination solution candidates.
  • step S13 The simulation executed in step S13 (hereinafter referred to as “full calculation”) is performed using the first series data (8,760 time data) including all times.
  • step S18 the simulation executed in step S18 (hereinafter referred to as “extraction calculation”) is performed using the essential series data (less than 8760 time data) extracted from the first series data.
  • the calculation time required for the extraction calculation is shorter than the calculation time required for the full calculation. Therefore, until the first local solution is obtained, the training data is recorded by executing the full calculation for each system configuration (combination solution candidate), and after the first local solution is obtained, each system configuration is extracted.
  • the evaluation index can be obtained efficiently by executing the calculation. Therefore, by executing the simulation (extraction calculation) for a large number of system configurations, the calculation time required to obtain the system configuration (optimal solution) that minimizes the evaluation index can be suppressed.
  • the essential series data extraction unit 150 includes information on m (m ⁇ 2 and n> m) combination solution candidates among the n (n ⁇ 2) combination solution candidates generated by the solution candidate generation unit 110.
  • the second series data is extracted from the first series data based on the evaluation index using the first series data.
  • the second series data includes time data necessary for calculating the evaluation index with a predetermined accuracy. In other words, the second series data includes time data having a large contribution to the simulation data (power consumption).
  • the "m candidate for combinatorial solution" is information about the system configuration generated until the first local solution is obtained, in other words, information about the system configuration used for generating training data.
  • the essential series data extraction unit 150 selects the third series data necessary for verifying the suitability of the predetermined constraint conditions for the above-mentioned "m candidates for combination solutions" from the first series data. Extract.
  • the predetermined constraint condition is a condition relating to an untreated heat load.
  • the third series data includes, for example, time data in which the unprocessed heat load is likely to exceed a predetermined reference value. Even if the power consumption calculated by the simulation is small, if the unprocessed heat load is large, it is not preferable from the viewpoint of maximizing the use of resources. Therefore, an appropriate evaluation index can be calculated by adding a value corresponding to the untreated heat load as a penalty to the total cost.
  • a predetermined constraint condition includes a condition that the untreated heat load is equal to or less than a predetermined reference value.
  • the essential series data extraction unit 150 extracts the first series data, which is likely not to satisfy the predetermined constraint condition, as the third series data.
  • the essential series data extraction unit 150 acquires the essential series data by combining the second series data and the third series data.
  • the simulation unit 120 includes at least information on combination solution candidates other than the above-mentioned "m combination solution candidates" and series data including the essential series data. And, are used to calculate the simulation data.
  • the combinatorial solution candidates other than "m combinatorial solution candidates" are information about the system configuration generated after the first local solution is obtained.
  • the essential series data extraction unit 150 extracts the second series data by, for example, sparse estimation, and extracts the third series data by the extreme value statistics.
  • Sparse estimation is one of the methods for estimating the regression coefficient in regression analysis.
  • Sparse estimation is a method in which regression coefficient estimation and variable selection are performed at the same time in order to estimate the objective variable based on the value of the explanatory variable as few as possible when there are many explanatory variables.
  • Extreme value statistics are, in probability theory and statistics, samples X 1 , X 2 , of size n generated according to a certain cumulative distribution function. .. .. , X n , which represents how the number of x or more (or x or less) is distributed, is a method of estimating an extreme value based on a continuous probability distribution model.
  • various patterns can be considered.
  • the indoor unit maintains a comfortable temperature by removing the sensible heat applied in each zone, so an indoor unit that can handle more than the sensible heat load required in each zone is required. Will be done.
  • various combinations can be considered, such as arranging one indoor unit having a large sensible heat load that can be processed, or arranging a plurality of indoor units of a type having a small sensible heat load that can be processed.
  • the ventilation system must meet the ventilation volume required in one zone, and the introduction of the ventilation system changes the sensible heat load required in each zone.
  • various patterns can be considered depending on what kind of outdoor unit is to be introduced first, and then which zone the indoor unit is connected to form the air conditioning system.
  • the load that can be handled by the indoor units and ventilation equipment placed in each zone must exceed the load required in that zone, and the number of indoor units that can be placed in each zone. It is also added as a condition that a certain upper limit is not exceeded and that the number of indoor units that can be connected to the outdoor unit does not exceed a certain upper limit.
  • a certain upper limit is not exceeded and that the number of indoor units that can be connected to the outdoor unit does not exceed a certain upper limit.
  • an indoor unit with a load exceeding the required load in the zone is placed, there may be a time when an unprocessed heat load is generated depending on the connection state of the indoor unit, so it occurs here.
  • the condition that the untreated sensible heat load to be applied is less than or equal to the specified value is also added.
  • the sensible heat load that can be processed by the indoor units arranged in each zone is greater than or equal to the sensible heat load applied in that zone.
  • B Observe the upper limit of indoor units placed in each zone.
  • C Observe the upper limit of the number of indoor units that can be connected to the outdoor unit.
  • D Observe the upper limit of installation of outdoor units.
  • E The ventilation volume that can be processed by the ventilation system arranged in each zone is equal to or greater than the ventilation load applied in that zone.
  • the sensible heat load that could not be processed in each zone is less than or equal to the specified value.
  • the 0-1 integer programming problem is an optimization problem in which all variables take a value of 0 or 1.
  • constraints (A) to (F) to be considered in Section (4-2-1) (A), (E), and (F) are treated so that the conditions can be relaxed, and the amount that violates the constraints. Is minimized by the objective function.
  • y i, v , i ⁇ Z , v ⁇ V Variables that are 1 when the ventilator v is located in zone i, 0 otherwise x: x i, j, k vector notation y: y i, v vector notation ⁇ k , i ⁇ Z: Variable indicating whether an outdoor unit of type k is installed (variable represented by equation 1)
  • the objective functions f (x) and h (y) are the power consumption obtained as output by inputting the variables x i, j, k , y i, v to the simulation software, respectively, and C (f (x), h (y). h (y)) is a function that calculates the power charge from the power consumption.
  • G (y) in the constraint equation is the load obtained as an output by inputting the variables y i and v into the simulation software.
  • U t (x) in the constraint equation is the unprocessed sensible heat load at time t obtained as an output by inputting the variables x i, j, k to the simulation software.
  • a constrained black box optimization problem is transformed into an unconstrained optimization problem.
  • a method for changing to an unconstrained optimization problem a method of removing those that do not satisfy the constraint conditions during the search from the search space and a method of using the penalty function method are considered.
  • the penalty function method was applied to the constraints (A), (E), and (F), and those that violated the constraints (B) to (D) were excluded from the search space.
  • the penalty function method is a method of unconstrained optimization of a penalty function defined by adding a penalty term for not satisfying the constraint condition to the objective function.
  • the penalty function F (x) is defined as follows. ⁇ is a positive parameter.
  • the solution obtained by the equation of equation 6 is a function obtained by applying Lagrange's undetermined multiplier method to the equation of equation 6. Is equivalent to the solution obtained by minimizing with respect to the parameter ⁇ .
  • the parameters ( ⁇ , ⁇ , ⁇ ) of the generalized extreme value distribution GEV ( ⁇ , ⁇ , ⁇ ) are estimated by the maximum likelihood method by applying the generalized extreme value distribution to the block maximum data.
  • the generalized extreme value distribution does not satisfy the regular condition for the maximum likelihood estimator, but when ⁇ > 0.5, it has been shown that the maximum likelihood estimator is an asymptotic effective estimator, which has been applied from previous experiments. Since it is rare that the above ⁇ ⁇ 0.5, parameter estimation is generally performed by the maximum likelihood method.
  • the maximum value z that can cause the upper p p ( ⁇ , Y) probability point z p of the generalized extreme value distribution G that the period maximum data follows in order to suppress the probability that the annual maximum data exceeds z p in the period Y.
  • max the reference value M of untreated sensible heat
  • Table 1 shows the number of zones, the number of types of indoor units, and the number of types of outdoor units from the top.
  • Table 2 shows the ventilatory loads that can be processed for each ventilator type.
  • Table 3 shows the sensible heat load capacity that can be processed in each indoor unit type.
  • Table 4 shows the sensible heat load and the ventilation load applied to each zone.
  • the configuration of the air conditioner with the lowest power consumption is derived by using the iterative local search method, and at the same time, the calculation time is reduced by using the sparse estimation and the extreme value statistics.
  • the calculation environment is Intel (R) Core TM i7 CPU 2.5GHz, and the memory is 8GB.
  • the correlation coefficient between the estimated value and the realized value was 0.9999867. From this, it can be seen that the hierarchical relationship of the solution is hardly broken, and it can be confirmed that it is possible to perform a local search using the estimated value.
  • FIGS. 5 to 7 shows unprocessed sensible heat at all times for the air conditioner configuration searched for the 3642nd time.
  • the calculation time was reduced by 73% compared to the full calculation.
  • the air conditioner configuration problem is a black box optimization problem in which the objective function and constraints can be understood by simulation, and cannot be formulated. Therefore, in the present disclosure, optimization is performed using an iterative local search method.
  • iterative local search method By using the iterative local search method, it was possible to derive a configuration with a lower cost than before.
  • the calculation cost by the simulation software is high, and the simple iterative local search method takes an enormous amount of time.
  • the time to be calculated is changed from 8760 variables to 27 variables by estimating the objective function value from some dates and times using the data until the first local solution is reached as training data. It could be shortened, mae was 31.23, and the correlation coefficient was 0.9999867, which made a very accurate estimation.
  • the constraint conditions are not taken into consideration only by the time extraction by sparse estimation, there is a possibility that unprocessed sensible heat may occur at the time when the time is not extracted. Therefore, by extracting the time that may violate the constraint condition by the extremum statistics, the variables that may violate the constraint condition are extracted, and 172 variables that are combined with the time extracted by the sparse estimation to eliminate duplication. By performing the extraction calculation in, the calculation time could be reduced by 73%.
  • the combination solution determination system 190 can reduce the calculation cost required for the simulation. Therefore, the combination solution determination system 190 can reduce the calculation time and efficiently derive the solution of the black box optimization problem.
  • the essential series data extraction unit 150 may acquire the second series data by excluding the data unnecessary for calculating the evaluation index with a predetermined accuracy from the first series data.
  • the essential series data extraction unit 150 may acquire the third series data by excluding unnecessary data from the first series data when considering a predetermined constraint condition.
  • the combination solution determined by the combination solution determination system 190 is a combination solution of objects or a combination solution of methods.
  • the combination solution may be a combination solution of methods such as air conditioning control parameters and air conditioning control contents.
  • the combination solution determination system 190 may be used in addition to the air conditioning system.
  • the combination solution determination system 190 can be used as a plant design system.
  • the combination solution determination system 190 may be used to determine the optimum capacity of each device of the plant that realizes the minimization of the plant operating cost, or to determine the optimum operation method of each device. ..
  • the combinational solution determined by the combinational solution determination system 190 may improve the situation of the application target of the combinational solution.
  • the application target of the combination solution is the air conditioning target space.
  • the first series data is past data on the situation that affects the improvement of the situation to which the combination solution is applied.
  • the situation is, for example, the temperature, humidity and carbon dioxide concentration of the air-conditioned space.
  • Past data are, for example, heat load and ventilation load.
  • the essential series data extraction unit 150 extracts the second series data by sparse estimation.
  • the essential sequence data extraction unit 150 may extract the second sequence data by a method other than sparse estimation.
  • Other methods include multiple regression analysis, principal component analysis, discriminant analysis and random forest.
  • the essential series data extraction unit 150 extracts the third series data by the extreme value statistics. However, the essential series data extraction unit 150 may extract the third series data by a method other than the extremum statistics.
  • An example of another method is a method of extracting conditions exceeding a predetermined value at a certain frequency.
  • Another example of another method is to apply the generated values to a normal distribution, calculate the mean and variance, and extract conditions that exceed a predetermined value with a certain probability.
  • the combination solution determination system can efficiently derive the solution of the black box optimization problem.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Strategic Management (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Combustion & Propulsion (AREA)
  • Computational Linguistics (AREA)
  • Chemical & Material Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Mechanical Engineering (AREA)
  • Marketing (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Air Conditioning Control Device (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Abstract

ブラックボックス最適化問題の解を効率的に導出することができる組合せ解決定システムを提供する。組合せ解決定システム（１９０）は、解候補生成部（１１０）と、シミュレーション部（１２０）と、評価指標算出部（１３０）と、解決定部（１４０）と、必須系列データ抽出部（１５０）とを備える。シミュレーション部は、解候補生成部が生成した組合せ解の候補に関する情報と、系列データとを用いてシミュレーションデータを演算する。評価指標算出部は、シミュレーションデータに基づいて評価指標を算出する。解決定部は、評価指標に基づいて、複数の組合せ解の候補の中から評価の高い組合せ解を決定する。必須系列データ抽出部は、系列データ（第１系列データ）の中から、所定の精度で評価指標を算出するために必要な第２系列データを抽出する。必須系列データ抽出部は、第１系列データの中から、所定の制約条件の適否を検証するために必要な第３系列データを抽出する。必須系列データ抽出部は、第２系列データと第３系列データとを合わせて必須系列データを取得する。シミュレーション部は、必須系列データ抽出部が必須系列データを抽出した後に、組合せ解の候補に関する情報と、必須系列データとを用いてシミュレーションデータを演算する。

Description

組合せ解決定システム
　組合せ解決定システム
　従来、特許文献１（特開２００６－４８４７５号公報）に開示されるように、資源を最大限に活用するために、現実の様々な問題を数理計画問題として定式化する方法が用いられている。
　数理計画問題が、目的関数を定式化できないブラックボックス最適化問題である場合、最適解の導出に多くの計算時間が必要となる課題がある。
　第１観点の組合せ解決定システムは、解候補生成部と、シミュレーション部と、評価指標算出部と、解決定部と、必須系列データ抽出部とを備える。解候補生成部は、組合せ解の候補を生成する。シミュレーション部は、解候補生成部が生成した組合せ解の候補に関する情報と、組合せ解を評価するための系列データと、を用いてシミュレーションデータを演算する。評価指標算出部は、シミュレーションデータに基づいて評価指標を算出する。解決定部は、複数の組合せ解の候補のそれぞれから評価指標算出部が算出した評価指標に基づいて、複数の組合せ解の候補の中から評価の高い組合せ解を決定する。必須系列データ抽出部は、解候補生成部が生成したｎ個（ｎ≧２）の組合せ解の候補のうちのｍ個（ｍ≧２）の組合せ解の候補に関する情報と、系列データとしての第１系列データと、から評価指標算出部が算出した評価指標に基づいて、第１系列データの中から所定の精度で評価指標を算出するために必要な第２系列データを抽出する。必須系列データ抽出部は、第１系列データの中から、当該ｍ個の組合せ解の候補に対して所定の制約条件の適否を検証するために必要な第３系列データを抽出する。必須系列データ抽出部は、第２系列データと第３系列データとを合わせて必須系列データを取得する。シミュレーション部は、必須系列データ抽出部が必須系列データを抽出した後に、少なくとも、当該ｎ個の組合せ解の候補のうちの当該ｍ個の組合せ解の候補以外の組合せ解の候補に関する情報と、必須系列データを含む系列データと、を用いてシミュレーションデータを演算する。
　第１観点の組合せ解決定システムは、ブラックボックス最適化問題の解を効率的に導出することができる。
　第２観点の組合せ解決定システムは、第１観点の組合せ解決定システムであって、必須系列データ抽出部は、所定の制約条件を満たさなくなる可能性が高い第１系列データを、第３系列データとして抽出する。
　第２観点の組合せ解決定システムは、ブラックボックス最適化問題の解を効率的に導出することができる。
　第３観点の組合せ解決定システムは、第１観点又は第２観点の組合せ解決定システムであって、必須系列データ抽出部は、スパース推定によって第２系列データを抽出し、極値統計によって第３系列データを抽出する。
　第３観点の組合せ解決定システムは、ブラックボックス最適化問題の解を効率的に導出することができる。
　第４観点の組合せ解決定システムは、第１乃至第３観点のいずれか１つの組合せ解決定システムであって、必須系列データ抽出部は、所定の精度で評価指標を算出するために不要なデータと、所定の制約条件を考慮するときに不要なデータとを、第１系列データから排除することによって、必須系列データを取得する。
　第４観点の組合せ解決定システムは、ブラックボックス最適化問題の解を効率的に導出することができる。
　第５観点の組合せ解決定システムは、第１乃至第４観点のいずれか１つの組合せ解決定システムであって、組合せ解は、物の組合せ解、又は、方法の組合せ解である。
　第５観点の組合せ解決定システムは、ブラックボックス最適化問題の解を効率的に導出することができる。
　第６観点の組合せ解決定システムは、第１乃至第５観点のいずれか１つの組合せ解決定システムであって、組合せ解は、組合せ解の適用対象の状況を改善するものであり、第１系列データは、組合せ解の適用対象の当該状況の改善に影響を与える、当該状況に関する過去のデータである。
　第６観点の組合せ解決定システムは、ブラックボックス最適化問題の解を効率的に導出することができる。
空調システム１００の構成の一例を示す図である。 組合せ解決定システム１９０の概略的なブロック図である。 組合せ解決定システム１９０による空調システムを選定する処理のフローチャートである。 反復局所探索における時刻抽出を表すグラフである。 １１月１９日１０時における未処理顕熱を表すグラフである。 １１月１９日１１時における未処理顕熱を表すグラフである。 １１月１９日１２時における未処理顕熱を表すグラフである。 ３６４２回目に計算した構成における未処理顕熱を表すグラフである。
　実施形態に係る組合せ解決定システムについて、図面を参照しながら説明する。組合せ解決定システムは、消費電力が最小となる空調システムを導出する数理計画問題の解を決定する。
　（１）空調システムの構成
　空調システムは、建物に設置され、主として室外機、室内機及び換気装置から構成される。図１は、空調システム１００の構成の一例を示す図である。図１において、空調システム１００は、２台の室外機１０ａ，１０ｂと、６台の室内機２０ａ～２０ｆと、２台の換気装置３０ａ，３０ｂとから構成される。各室外機１０ａ，１０ｂは、１台または複数台の室内機２０ａ～２０ｆと冷媒配管を介して接続されている。各室内機２０ａ～２０ｆは、２台の室外機１０ａ，１０ｂのいずれか１つに接続され、かつ、空調システム１００が設置されている建物のゾーン４０ａ～４０ｃのいずれか１つに設置されている。図１では、各室外機１０ａ，１０ｂは３台の室内機２０ａ～２０ｆと接続され、各ゾーン４０ａ～４０ｃに２台の室内機２０ａ～２０ｆが設置されている。ゾーン４０ａ～４０ｃは、空調システム１００による空調対象空間である。室内機２０ａ～２０ｆは、ゾーン４０ａ～４０ｃの顕熱を取り除くことによって、ゾーン４０ａ～４０ｃを快適な状態に維持する。換気装置３０ａ，３０ｂは、１つ又は複数のゾーン４０ａ～４０ｃを換気して、ゾーン４０ａ～４０ｃを快適な状態に維持する。
　（２）空調システムの選定
　空調システム１００を建物に設置する前に、空調システム１００の選定を行う必要がある。空調システム１００の選定とは、ゾーン４０ａ～４０ｃの熱負荷、及び、消費電力等を考慮して、機器（室外機１０ａ，１０ｂ、室内機２０ａ～２０ｆ及び換気装置３０ａ，３０ｂ）を選択したり、機器の組み合わせを決定したりすることである。空調システム１００の選定で決定されるパラメータ（空調選定パラメータ）は、例えば、各ゾーン４０ａ～４０ｃの室内機２０ａ～２０ｆの台数、室内機２０ａ～２０ｆの機種及び性能（容量等）、室外機１０ａ，１０ｂの機種及び性能（容量等）、換気装置３０ａ，３０ｂの機種及び性能（換気量等）、冷媒系統、制御パラメータ等である。冷媒系統とは、例えば、室外機１０ａ，１０ｂと、当該室外機１０ａ，１０ｂに接続される室内機２０ａ～２０ｆとの組に関する情報である。制御パラメータとは、例えば、各ゾーン４０ａ～４０ｃの設定温度及び設定湿度である。
　従来、空調システムを選定する際には、ピーク時の熱負荷に対応できるように過大な性能を有する機器が選択されている。しかし、温度を固定値とした静的な熱負荷計算しか行われなかったり、現場の担当者の経験及び技量によって選定される機器にばらつきがあったりして、消費電力の観点からは適切でない機器が選定される課題がある。この課題を解決するためには、機器及び建物の特性、及び、空調対象空間の温度及び顕熱の変化を考慮して空調システムを選定する必要がある。しかし、その場合、熱負荷及び消費電力の計算コストが大きくなる。また、空調システムの選定時において、機器の数及び種類が多いほど、選定に要する計算時間が増大する。
　実施形態に係る組合せ解決定システムは、上記の課題を解決するものであり、多数の選択肢の中から消費電力が最小となる空調システムを効率的に選定する。これにより、組合せ解決定システムは、トータルコストを最小化する空調システムを選定することができる。トータルコストとは、機器コストと電気コストとの合計である。機器コストは、例えば、機器自体の費用、機器の設置工事費用、及び、機器の保守費用を含む。電気コストは、空調システムを稼動させるために必要な電気料金を含み、機器の消費電力等から算出される。
　実施形態に係る組合せ解決定システムによって解が決定される数理計画問題は、目的関数を定式化できないため汎用ソルバーを用いることができないブラックボックス最適化問題に該当する。ブラックボックス最適化問題の解を決定するためには、一般的に、目的関数の値をシミュレーションにより計算する必要があるため、計算コストがかかりやすい。例えば、消費電力が最小となる空調システムを選定する問題の場合、多数の時刻における消費電力を、シミュレーションソフトを用いて計算する必要があり、計算時間が長くなりやすい。実施形態に係る組合せ解決定システムは、ブラックボックス最適化問題にも用いることができる局所探索法の応用手法である反復局所探索法を用いることで、計算時間を削減することができる。
　（３）組合せ解決定システムの構成
　図２は、組合せ解決定システム１９０の概略的なブロック図である。組合せ解決定システム１９０は、消費電力が最小となる空調システムを選定する。組合せ解決定システム１９０は、例えば、１つ又は複数のコンピュータから構成される。組合せ解決定システム１９０が複数のコンピュータから構成される場合、当該複数のコンピュータは、ネットワークを介して互いに接続されてもよい。
　組合せ解決定システム１９０は、主として、解候補生成部１１０と、シミュレーション部１２０と、評価指標算出部１３０と、解決定部１４０と、必須系列データ抽出部１５０とを備える。解候補生成部１１０～必須系列データ抽出部１５０は、例えば、組合せ解決定システム１９０を構成するコンピュータの記憶装置に記憶されているプログラムを、当該コンピュータのＣＰＵが実行することにより実現される。
　解候補生成部１１０は、組合せ解の候補を生成する。組合せ解とは、組合せ解決定システム１９０によって決定される解である。組合せ解は、空調システム１００の構成に関する情報であり、具体的には、上述した空調選定パラメータを含む。
　シミュレーション部１２０は、解候補生成部１１０が生成した組合せ解の候補に関する情報と、組合せ解を評価するための系列データと、を用いてシミュレーションデータを演算する。シミュレーションデータは、組合せ解の候補を入力パラメータとしてシミュレーションにより算出される、空調システム１００の消費電力を含む。系列データとは、所定の時刻におけるシミュレーションデータ（消費電力）が算出される場合、当該時刻に関するデータである。シミュレーション部１２０が１年間の消費電力をシミュレーションにより算出する場合、系列データは、例えば、当該１年間の１時間ごとの時刻である。この場合、系列データは、８７６０個の時刻データから構成される（３６５日×２４時間／日）。シミュレーション部１２０が算出する消費電力は、室外機１０ａ，１０ｂ、室内機２０ａ～２０ｆ、及び、換気装置３０ａ，３０ｂの消費電力である。換気装置３０ａ，３０ｂの消費電力は、換気シミュレーションによって算出される。室外機１０ａ，１０ｂ及び室内機２０ａ～２０ｆの消費電力は、換気装置３０ａ，３０ｂによる熱負荷の影響を考慮して、各ゾーン４０ａ～４０ｃの熱負荷（顕熱負荷）から空調シミュレーションによって算出される。
　評価指標算出部１３０は、シミュレーション部１２０が算出したシミュレーションデータ（消費電力）に基づいて評価指標を算出する。評価指標とは、上述のトータルコスト、未処理熱負荷、及び、未処理換気量等から算出されるパラメータである。未処理熱負荷とは、ゾーン４０ａ～４０ｃに配置された室内機２０ａ～２０ｆが処理できる顕熱負荷がそのゾーン４０ａ～４０ｃでかかる顕熱負荷未満である場合に、処置できない分の顕熱負荷に相当する。未処理換気量とは、ゾーン４０ａ～４０ｃに配置された換気装置３０ａ，３０ｂが処理できる換気負荷がそのゾーン４０ａ～４０ｃでかかる換気負荷未満である場合に、処置できない分の換気負荷に相当する。トータルコスト、未処理熱負荷、及び、未処理換気量が小さいほど、評価指標は小さくなる。資源を最大限に活用する観点からは、評価指標は小さいほど好ましい。
　解決定部１４０は、複数の組合せ解の候補のそれぞれから評価指標算出部１３０が算出した評価指標に基づいて、複数の組合せ解の候補の中から評価の高い組合せ解を決定する。評価の高い組合せ解とは、評価指標が最も小さい組合せ解である。解決定部１４０によって決定された組合せ解は、資源を最大限に活用する観点から最適な空調選定パラメータである。
　必須系列データ抽出部１５０は、必須系列データを取得する。必須系列データとは、系列データから抽出されたデータである。具体的に説明すると、上述したように、系列データが１時間ごとの８７６０個の時刻データから構成される場合、必須系列データは、この８７６０個の時刻データ（以下「第１系列データ」と呼ぶ。）の一部の時刻データである。
　シミュレーション部１２０は、必須系列データ抽出部１５０が必須系列データを抽出した後に、組合せ解の候補に関する情報と、必須系列データと、を用いてシミュレーションデータを演算する。そのため、シミュレーション部１２０は、第１系列データに含まれる時刻データよりも少ない時刻データに基づいて、シミュレーションデータ（消費電力）を演算することができる。これにより、組合せ解決定システム１９０は、消費電力の演算に必要な計算コストを削減することができる。
　ここで、図３を参照しながら、組合せ解決定システム１９０が計算コストを削減することができる理由について説明する。図３は、組合せ解決定システム１９０による空調システムを選定する処理のフローチャートである。
　ステップＳ１１では、システム構成が決定される。具体的には、ステップＳ１１では、解候補生成部１１０によって組合せ解の候補が生成される。
　ステップＳ１２では、必須系列データが抽出済みであるか否かが判定される。ステップＳ１６において必須系列データが抽出済みである場合、ステップＳ１３に移行し、抽出済みでない場合、ステップＳ１８に移行する。
　ステップＳ１３では、シミュレーションが実行される。具体的には、シミュレーション部１２０は、組合せ解の候補と、第１系列データとを用いて、第１系列データに含まれる各時刻におけるシミュレーションデータを演算する。
　ステップＳ１４では、ステップＳ１３で得られた各時刻におけるシミュレーションデータ（消費電力）を、トレーニングデータとして保存する。
　ステップＳ１５では、評価指標算出部１３０によって、ステップＳ１４で保存されたトレーニングデータから評価指標が算出されて保存される。
　ステップＳ１６では、トレーニングデータから算出された評価指標に基づいて、最初の局所解が導出されたか否かが判定される。最初の局所解は、反復局所探索法によって算出される。最初の局所解は、例えば、シミュレーションの実行開始から、評価指標が最初に最小値（局所的な最小値）を示したシステム構成である。
　ステップＳ１６において最初の局所解が導出された場合、ステップＳ１７において、それまでに得られたトレーニングデータに基づいて、第１系列データから必須系列データが抽出される。具体的には、必須系列データ抽出部１５０は、消費電力の演算に必要な時刻を抽出する。その後、ステップＳ１１において新たなシステム構成が決定されて、シミュレーションが実行される。
　ステップＳ１６において最初の局所解が導出されなかった場合、ステップＳ１７において必須系列データが抽出されずに、ステップＳ１１において新たなシステム構成が決定されて、トレーニングデータがさらに取得される。
　ステップＳ１８では、ステップＳ１７で抽出された時刻（必須系列データ）に基づいて、シミュレーションが実行される。具体的には、シミュレーション部１２０は、組合せ解の候補と、必須系列データとを用いて、必須系列データに含まれる各時刻におけるシミュレーションデータを演算する。
　ステップＳ１９では、評価指標算出部１３０によって、ステップＳ１８で得られたシミュレーションデータから評価指標が算出されて保存される。
　以降、ステップＳ１１でシステム構成が新たに決定され、ステップＳ１８においてシミュレーションが実行され、ステップＳ１９において評価指標が算出されて保存される。以上の処理が、所定の数の組合せ解の候補について行われる。
　ステップＳ１３において実行されるシミュレーション（以下「フル計算」と呼ぶ。）は、全ての時刻を含む第１系列データ（８７６０個の時刻データ）を用いて行われる。一方、ステップＳ１８において実行されるシミュレーション（以下「抽出計算」と呼ぶ。）は、第１系列データから抽出された必須系列データ（８７６０個未満の時刻データ）を用いて行われる。抽出計算に要する計算時間は、フル計算に要する計算時間よりも短い。そのため、最初の局所解が得られるまでは、各システム構成（組合せ解の候補）についてフル計算を実行してトレーニングデータを記録し、最初の局所解が得られた後は、各システム構成について抽出計算を実行して効率的に評価指標を得ることができる。従って、多数のシステム構成についてシミュレーション（抽出計算）を実行することで、評価指標が最小となるシステム構成（最適解）を得るために必要な計算時間を抑えることができる。
　次に、図３のステップＳ１７において第１系列データの中から必須系列データを抽出する具体的な処理の概要について説明する。
　必須系列データ抽出部１５０は、解候補生成部１１０が生成したｎ個（ｎ≧２）の組合せ解の候補のうちのｍ個（ｍ≧２かつｎ＞ｍ）の組合せ解の候補に関する情報と、第１系列データとを用いて、評価指標に基づいて、第１系列データの中から第２系列データを抽出する。第２系列データとは、所定の精度で評価指標を算出するために必要な時刻データを含む。言い換えると、第２系列データとは、シミュレーションデータ（消費電力）に対する寄与が大きい時刻データを含む。「ｍ個の組合せ解の候補」とは、最初の局所解が得られるまでに生成されたシステム構成に関する情報であり、言い換えると、トレーニングデータの生成に利用されたシステム構成に関する情報である。
　さらに、必須系列データ抽出部１５０は、第１系列データの中から、上述の「ｍ個の組合せ解の候補」に対して所定の制約条件の適否を検証するために必要な第３系列データを抽出する。所定の制約条件とは、未処理熱負荷に関する条件である。第３系列データとは、例えば、未処理熱負荷が所定の基準値を超える可能性が高い時刻データを含む。シミュレーションにより算出された消費電力が小さくても、未処理熱負荷が大きい場合、資源を最大限に活用する観点からは好ましくない。そのため、未処理熱負荷に応じた値をペナルティーとしてトータルコストに加えることで、適切な評価指標を算出することができる。所定の制約条件としては、未処理熱負荷が所定の基準値以下であるという条件が挙げられる。この場合、必須系列データ抽出部１５０は、所定の制約条件を満たさなくなる可能性が高い第１系列データを、第３系列データとして抽出する。
　そして、必須系列データ抽出部１５０は、第２系列データと第３系列データとを合わせて必須系列データを取得する。シミュレーション部１２０は、必須系列データ抽出部１５０が必須系列データを抽出した後に、少なくとも、上述の「ｍ個の組合せ解の候補」以外の組合せ解の候補に関する情報と、必須系列データを含む系列データと、を用いてシミュレーションデータを演算する。「ｍ個の組合せ解の候補」以外の組合せ解の候補とは、最初の局所解が得られた後に生成されたシステム構成に関する情報である。
　必須系列データ抽出部１５０は、例えば、スパース推定によって第２系列データを抽出し、極値統計によって第３系列データを抽出する。
　スパース推定とは、回帰分析における回帰係数の推定方法の一つである。スパース推定は、説明変数が多い場合になるべく少ない説明変数の値をもとに目的変数を推定するため、回帰係数の推定と変数選択とを同時に行う手法である。
　極値統計とは、確率論および統計学において、ある累積分布関数にしたがって生じた大きさｎの標本Ｘ１，Ｘ２，．．．，Ｘｎのうち、ｘ以上（あるいはｘ以下）となるものの個数がどのように分布するかを表す、連続確率分布モデルに基づき極値を推定する手法である。
　（４）組合せ解決定システムの具体例
　（４－１）概要
　次に、具体例として、消費電力が最小となる空調システムを選定するための組合せ解決定システムについて詳細に説明する。この組合せ解決定システムは、スパース推定と極値統計により消費電力を計算すべき時刻を抽出し、抽出された時刻から一年間の消費電力の推定を行うことで、全体の計算コスト削減を目指す。
　（４－２）空調機構成問題
　ある建物に室内機を導入する際、必ずどの室内機も室外機と接続されている必要がある。一つの室外機に接続できる室内機の数には制限があるが、一つの室外機には複数台の室内機を接続することができる。このとき、建物に導入されている室外機と、それに接続されている室内機との組を空調システムという。
　一つの建物に空調機を設置する際、様々なパターンが考えられる。まず、室内機を設置する際、室内機は各ゾーンでかかる顕熱を取り除くことによって快適な温度を維持しているので、各ゾーンにおいて必要となる顕熱負荷以上を処理できる室内機が必要とされる。室内機の組み合わせ方としては、処理できる顕熱負荷が大きな室内機を一つ置くというパターンや、処理できる顕熱負荷が小さな種類の室内機を複数配置するなど様々な組み合わせが考えられる。また、換気装置は一つのゾーンで必要となる換気量を満たす必要があり、換気装置の導入により各ゾーンで必要となる顕熱負荷も変化する。次に、室外機を配置する際、まずどの種類の室外機を導入するのか、次にどのゾーンの室内機と接続し空調システムを構成するのかに対しても様々なパターンが考えられる。
　どの空調システムを導入するかによって消費電力が変化する。このとき消費電力が最小となる空調システムを導出する問題を空調機構成問題という。本開示では、消費電力を１５年分の電気料金に換算したランニングコストと、空調機の値段（イニシャルコスト）とを加えたものの最小化を目指す。
　（４－２－１）空調機構成において考慮すべき制約
　空調機を配置する上で必要な制約条件について説明する。
　まず、各ゾーンに配置されている室内機と換気装置が処理することのできる負荷が、そのゾーンで必要とされる負荷を上回ることが必要であり、また、各ゾーンに配置できる室内機の数がある上限を超えないことと、室外機に接続できる室内機の数もある上限を超えないことも条件として加える。また、ゾーンで必要となる負荷以上の室内機を配置しても、その室内機の接続状態などにより未処理の熱負荷が発生してしまう時刻が存在してしまうことがあるため、ここで発生する未処理顕熱負荷は規定値以下であるという条件も加える。上記の考慮すべき点を以下に整理する。
　（Ａ）各ゾーンに配置された室内機が処理することのできる顕熱負荷が、そのゾーンでかかる顕熱負荷以上である。
　（Ｂ）各ゾーンに配置される室内機の上限を守る。
　（Ｃ）室外機に接続できる室内機の台数の上限を守る。
　（Ｄ）室外機の設置上限を守る。
　（Ｅ）各ゾーンに配置された換気装置が処理することのできる換気量は、そのゾーンでかかる換気負荷以上である。
　（Ｆ）各ゾーンで発生した処理できなかった顕熱負荷は、規定値以下である。
　（４－２－２）０－１整数計画問題としての定式化
　空調機構成問題を、０－１整数計画問題として定式化する。
　０－１整数計画問題とは、最適化問題において、すべての変数が０または１の値を取るもののことである。（４－２－１）節で挙げた考慮すべき制約（Ａ）～（Ｆ）のうち、（Ａ）、（Ｅ）、（Ｆ）は条件を緩和できるように扱い、制約を違反した量を目的関数で最小化する。
　（４－２－２－１）記号説明
　定式化に利用する記号について説明する。
　定数
I：室内機の集合
O：室外機の集合
Z：ゾーンの集合
V：換気装置の集合
T={1, 2, …, 8760}：時刻の集合
aj, j∈I：室内機jの処理可能な顕熱負荷
bi, i∈Z：ゾーンiでかかる顕熱負荷
fi, i∈Z：ゾーンiでかかる換気負荷
gv, v∈V：換気装置vの処理可能な換気負荷
pj, j∈I：室内機jの値段
qv, v∈V：換気装置vの値段
c：ゾーンに配置される室内機の上限
d：室外機に接続できる室内機の台数の上限
e：室外機の設置上限
M：未処理顕熱の基準値
　変数
xi,j,k, i∈Z, j∈I, k∈O：室内機jがゾーンiに配置され種類kの室外機につながっているとき１、そうでないとき０となる変数
yi,v, i∈Z, v∈V：換気装置vがゾーンiに配置されているとき１、そうでないとき０となる　変数
x：xi,j,kのベクトル表記
y：yi,vのベクトル表記
δk, i∈Z：種類kの室外機が設置されているかを表す変数（数１で表される変数）
Figure JPOXMLDOC01-appb-M000001
　（４－２－２－２）定式化
　（４－２－２－１）節の記号を利用し、（４－２－１）節の（Ａ）～（Ｄ）の制約を元に定式化したものを以下に示す。
Figure JPOXMLDOC01-appb-M000002
　目的関数f(x), h(y)は、変数xi,j,k, yi,vをそれぞれシミュレーションソフトに入力することにより出力として得られる消費電力であり、C(f(x), h(y))は、消費電力から電力料金を算出する関数である。
　制約式中のg(y)は、変数yi,vをシミュレーションソフトに入力することにより出力として得られる負荷である。
　制約式中のut(x)は、変数xi,j,kをシミュレーションソフトに入力することにより出力として得られる、時刻tにおける未処理顕熱負荷である。
　（４－２－３）制約付きブラックボックス最適化問題の解法
　一般的に、制約条件付きブラックボックス最適化問題は無制約最適化問題に変形する。無制約最適化問題へ変更するための手法として、探索の際制約条件を満たさないものは探索空間から外す方法と、ペナルティー関数法を用いる方法が考えられている。
　本問題では、制約（Ａ）、（Ｅ）、（Ｆ）に関してペナルティー関数法を適用し、制約（Ｂ）～（Ｄ）に違反するものは探索空間から外した。ペナルティー関数法とは、制約条件を満たさないことに対するペナルティー項を目的関数に加えて定義されるペナルティー関数を無制約最適化する手法である。
　ペナルティー関数F(x)を次のように定義する。
Figure JPOXMLDOC01-appb-M000003
 　ρは、正のパラメータである。
　（４－２－４）ブラックボックス最適化問題の解法
　本節ではブラックボックス最適化問題に対する解法について紹介し、本開示で取り組んだ解法の方針を説明する。（４－２－３）節で用いた方法により問題を無制約ブラックボックス最適化問題に変形し、F(x)に対し基本的なブラックボックス最適化問題に対する解法として知られている反復局所探索法を用いることによって解の改善を目指した。
　（４－３）計算時刻の抽出
　計算時刻を抽出するにあたりどのようなデータ構造をしているかというトレーニングデータが必要になる。しかし、空調機構成問題では空調機を配置する建物の種類や立地条件、気候などによりシミュレーション結果も変動するため、一意なトレーニングデータを作成することができない。そこで本開示では、図４に示されるように、反復局所探索法において最初の局所解を導出するまでは時刻抽出を行わず全ての時刻で計算を行い、最初の局所解導出までに得られたデータをトレーニングデータとして時刻抽出を行い、それ以降の探索では計算時刻を抽出して計算を行うことにより計算時間の削減を実現した。
　（４－３－１）スパース推定
　スパース推定は、ここ１０数年、情報学、機械学習、統計学など、様々な分野から注目を集めているが、本開示では、L1正則化法の代表である、Tibshiraniの提案したLASSO（Least Absolute Shrinkage and Selection Operator）を用いた。LASSOは、回帰モデルの損失関数にパラメータのL1ノルムに基づく正則化項を加えた正則化損失関数を最小化することによってパラメータを推定する方法で、推定の安定化とともに変数選択も行うことができる手法である。以下ではスパース推定の理論と本開示での適用方法を説明する。
　（４－３－１－１）スパース推定の理論
　連続値をとる目的関数Yとｐ次元説明変数x=(x1, …, xp)Tに関して、ｎ個の観測によりデータ(xi, yi); i=1, …, nが得られたとする。xi=(xi1, …, xip)Tとする。
Figure JPOXMLDOC01-appb-M000004
 を満たしているように正規化する。
　X=(x(1), …, x(p)), xi=(xi1, …, xnj)T y=(y1, …, yn)Tとおくと回帰モデルはつぎのように書ける。
Figure JPOXMLDOC01-appb-M000005
　線形回帰モデルの回帰係数の推定法としてlassoと呼ばれる次の制約付き最小化関数を考える。
Figure JPOXMLDOC01-appb-M000006
　lassoに基づいてパラメータ推定を行うことにより、いくつかのパラメータの推定値はぴったり０に縮小する性質を持つ。数６の式で得られる解は、数６の式に対してラグランジュの未定乗数法を適用することにより得られる関数
Figure JPOXMLDOC01-appb-M000007
 をパラメータβに関して最小化することにより得られる解と同値になっている。
　（４－３－１－２）本開示におけるスパース推定
　ｉ回目のシミュレーション計算における時刻ｔにおける消費電力をxit(t=1, …, 8760)、全時刻の消費電力の合計値をyiとすると線形回帰モデルとして下式のように書くことができ、
Figure JPOXMLDOC01-appb-M000008
 このモデルに対してスパース推定を行うことにより、回帰係数(β1, …, β8760)のうち多くの時刻にかかる係数は０となり、係数が０でない時刻のみシミュレーションソフトで消費電力を導出し係数をかけることにより、全時刻における消費電力を推定した。
　（４－３－２）極値統計
　本節では、未処理負荷が発生する可能性のある時刻を抽出するために用いる極値統計について説明する。極値統計学とは、もともと自然災害の予測や評価の際に用いられてきた学問であり、限られた期間の観測データから将来どのような大きな値をもつ事象が発生するか予測することを目的に考えられた。観測されていない右裾の領域について推測するため、極値統計では大きな値をとるデータのみに分布を当てはめることを考える。極値統計にはいくつかの統計モデルが存在するが、本開示ではGEVモデルを活用する。GEVモデルでは、ある期間でのブロック最大データに対し一般極値分布をあてはめて解析を行った。以下ではGEVモデルとその解析方法について紹介する。
　（４－３－２－１）極値理論
　はじめに、独立で同一分布Fに従う確率変数X1, X2, …について考える。ｎ個の確率変数の最大値を
Figure JPOXMLDOC01-appb-M000009
 とおく。適当な尺度との変換により基準化するとほとんどの連続分布が満たす条件の下で、Znは退化していない分布に収束することが知られている。
（Fisher-Tippettの定理）
　ある定数an>0, bn∈Rと退化していない分布Gが存在して、
Figure JPOXMLDOC01-appb-M000010
 を満たすならば、分布Gは次の標準極値分布Ｇξで表すことができる。
Figure JPOXMLDOC01-appb-M000011
　（定義）
　次の分布を一般極値分布といいGEV(μ,σ,ξ)(－∞＜μ＜∞, σ＞０, －∞＜ξ＜∞)で表す。
Figure JPOXMLDOC01-appb-M000012
　一般極値分布GEV(μ,σ,ξ)のパラメータ(μ,σ,ξ)はブロック最大データに一般極値分布をあてはめ、最尤法により推定する。一般極値分布は最尤推定量に関する正則条件を満たしていないが、ξ＞０．５の場合、最尤推定量は漸近有効推定量になることが示されており、これまでの実験から応用上ξ≦０．５になることはまれであるため、一般にパラメータ推定は最尤法で行われる。
　一般極値分布GEV(μ,σ,ξ)の上側p確率点zp
Figure JPOXMLDOC01-appb-M000013
 において、この確率点zpは再現期間1/pの再現レベルと呼ぶ。ここでは再現レベルについて考える。分布Gを、ある期間の中（一年間など）の最大データが従う一般極値分布GEV(μ,σ,ξ)とする。期間最大データは、互いに独立で同一の分布Gに従う確率変数Z1, Z2, …の実現値とする。このときi番目の期間を考えるとZiは期間最大なので、事象Zi＞zpは最大データが値zpを超えることを、事象Zi≦zpは観測したすべてのデータが値zpを超えないことを表す。ここで２値の確率変数
Figure JPOXMLDOC01-appb-M000014
 を考える。このB1, B2, …はベルヌーイ試行になる。G(zp)=１－ｐよりBiは互いに独立で二項分布B(1,p)に従う。このとき
Figure JPOXMLDOC01-appb-M000015
 はT年間で各年の最大データzpを超えた年の数を表しており、JT～B(T,p)よりその平均E(JT)=Tpとなり、特にT=1/pとするとE(J1/p)=1より再現レベルzpを超える現象が観測されるのは1/p年で平均一回である。
　（４－３－２－２）本開示における極値統計
　本開示におけるブロックは局所探索法で探索した10点のうち未処理負荷が最大となるデータとし、各時刻において極値分布の作成を行った。作成した極値分布に基づいて各時刻で今後出る可能性のある未処理負荷の最大値zmaxを推定する。ここでは、より安全なレベル期間zpとして期間Yにおいてzpを超える確率を小さな確率αで抑える。すなわち十分小さいαに対して
Figure JPOXMLDOC01-appb-M000016
 を満たすｐを考える。ただしMY=max{Z1, Z2, … ZY}である。
ここで与えられたαに対して
Figure JPOXMLDOC01-appb-M000017
 となるｐを計算する。すなわち期間Yで年最大データがzpを超える確率をα以下に抑えるために期間最大データが従う一般極値分布Gの上側p=p(α,Y)確率点zpを起こりえる最大値zmaxと定義する。未処理顕熱の基準値Mとして、zmax>Mである時刻の抽出計算を行い、zmax<Mである時刻は間引きを行った。
　（４－４）問題例
　（４－４－１）基本情報
　空調機割り当ての上での基本的な情報を、以下の表１、表２、表３、表４に示す。
　表１は、上からゾーン数、室内機の種類数、室外機の種類数を表している。
　表２は、各換気装置タイプにおける処理可能な換気負荷を表している。
　表３は、各室内機タイプにおける処理可能な顕熱負荷容量を表している。
　表４は、各ゾーンにかかる顕熱負荷と換気負荷を表している。
Figure JPOXMLDOC01-appb-T000018
Figure JPOXMLDOC01-appb-T000019
Figure JPOXMLDOC01-appb-T000020
Figure JPOXMLDOC01-appb-T000021
　（４－４－２）各制約条件
　（４－２）節で述べた各制約に関して下に示す。
（Ｂ）ゾーンに配置できる室内機は２台まで。
（Ｃ）一つの室外機に接続できる室内機は６台まで。
（Ｄ）室外機の設置上限は４台まで。
（Ｆ）未処理顕熱負荷の規定値は5.5kWまでとする。
　本開示において、反復局所探索法を用いることによって消費電力の最も小さくなる空調機の構成を導出すると同時に、スパース推定と極値統計を用いることにより計算時間の削減を試みた。計算環境はIntel(R) CoreTM i7 CPU 2.5GHz、メモリは8GBである。
　スパース推定と極値統計により抽出計算した結果の評価と、従来人が選定していた構成と本開示で導出した最適解の比較を行った。
　（４－５）計算実験
　（４－５－１）終了条件
　反復局所探索法における終了時条件は、反復回数が10000点探索した時終了とした。また、一度に探索する近傍は10点とし、10点のうち最善の解に移動するものとした。単純反復局所探索法の過程で350点解を探索しても解が改善しない場合、その解を局所解とし反復局所探索のステップへ移行するものとした。
　（４－５－２）時刻抽出
　最初の局所解に至るまでに1100回の解計算が行われ、この1100回分の計算データをトレーニングデータとしてスパース推定、極値統計による時刻抽出を行った。
　（４－５－２－１）スパース推定による時刻抽出
　スパース推定を行うことにより24時間365日、8760変数のうち27変数の抽出が行われた。スパース推定による消費電力推定におけるmae（mean absolute error）は31.23であった。全体のmaeは小さく、また得られた全ての局所解における実現値と推測値は次の表５のようになり、局所解における評価値の推測も正しくできていると考えられる。
Figure JPOXMLDOC01-appb-T000022
　推測値と実現値の相関係数は、0.9999867であった。これよりほとんど解の上下関係が崩れていないことがわかり推測値による局所探索を行うことが可能であることを裏付けることができている。
　（４－５－２－２）極値統計による時刻抽出
　極値統計により時刻抽出における基準値Mは安全側に作用させるために本来の基準値（5.5kW）よりも低い4kWを用いた。本開示では10000回の計算を行っており、極値統計では10回の計算結果の最大値を用いていることからY=1000, α=0.05とし一般極値分布Gの上側p(0.05, 1000)=5.13×10-5確率点zpを計算することにより、計算過程において5パーセントの確率で起こりえる値zmaxを計算することができ、zmax≧4である時刻を抽出した。安全側に基準値を下げたにもかかわらず本来の基準値である5.5kWを超えた日時で極値統計による抽出では包括できなかった時刻が[11/19 10:00:00, 11/19 11:00:00, 11/19 12:00:00]の３点存在した。この３時刻における全探索過程における未処理顕熱の推移を図５～７に示す。図５～７より基準値5.5を超えたのは3642回目の計算によるものであることがわかる。3642回目に探索した空調機構成に対する全時刻の未処理顕熱を図８に示す。
　図８よりそれぞれの構成における未処理顕熱を確認すると、どの構成においても基準値を大きく超える未処理顕熱が発生している時刻があり、極値統計で選ばれなかった時刻で基準値を超える未処理顕熱が現れるのは、特殊な構成に対してのみであり局所解に近づいていく時に発生する未処理顕熱は極値統計による抽出でカバーできていることがわかる。
　（４－５－２－３）最適解と計算時間
　本開示において導出された最適構成と、従来人が選択していた構成に対する評価値は以下の表６のようになっている。
Figure JPOXMLDOC01-appb-T000023
　表６を見るとわかる通り最適設計がどの従来設計よりもトータルコストの良い結果が導出できており、ペナルティーの値も０にできている。
　計算時間に関する比較を以下の表７に示す。
Figure JPOXMLDOC01-appb-T000024
　フル計算するときと比較し計算時間を73％削減することができた。
　（４－６）結論
　空調機構成問題は、シミュレーションにより目的関数や制約条件がわかるブラックボックス最適化問題であり定式化することができない。そのため、本開示では反復局所探索法を用いて最適化を行った。反復局所探索法を用いることによって、従来よりもコストの低い構成を導出することができた。しかし、シミュレーションソフトによる計算コストが大きいという問題点があり、単純な反復局所探索法では莫大な時間がかかってしまった。
　そこで本開示では、最初に局所解に陥るまでのデータをトレーニングデータとし、スパース推定を用いて一部の日時から目的関数値の推定を行うことにより、計算すべき時刻を8760変数から27変数に短縮することができ、maeは31.23、相関係数が0.9999867と非常に精度の高い推定を行うことができた。しかし、スパース推定による時刻抽出のみでは制約条件を考慮していないため抽出されなかった時刻で未処理顕熱が発生してしまう可能性がある。そこで、極値統計による制約条件に違反する可能性のある時刻抽出を行うことにより、制約条件を違反する可能性のある変数を抽出し、スパース推定で抽出した時刻と合わせ重複をなくした172変数での抽出計算を行うことにより、計算時間を73％削減することができた。また、目的関数値よりも、従来設計よりもコストの小さいものを導出でき、フルで計算を行ったときと遜色ない結果を得ることができた。
　（５）効果
　組合せ解決定システム１９０は、シミュレーションに必要な計算コストを削減することができる。そのため、組合せ解決定システム１９０は、計算時間を削減して、ブラックボックス最適化問題の解を効率的に導出することができる。
　（６）変形例
　（６－１）変形例Ａ
　必須系列データ抽出部１５０は、所定の精度で評価指標を算出するために不要なデータを、第１系列データから排除することによって、第２系列データを取得してもよい。
　必須系列データ抽出部１５０は、所定の制約条件を考慮するときに不要なデータを、第１系列データから排除することによって、第３系列データを取得してもよい。
　（６－２）変形例Ｂ
　組合せ解決定システム１９０が決定する組合せ解は、物の組合せ解、又は、方法の組合せ解である。実施形態の場合、組合せ解は、空調制御パラメータ及び空調制御内容等の方法の組合せ解であってもよい。
　また、組合せ解決定システム１９０は、空調システム以外にも用いられてもよい。例えば、組合せ解決定システム１９０は、プラントの設計システムに用いることができる。この場合、組合せ解決定システム１９０は、プラント運用費用の最小化を実現するプラントの各機器の最適容量を決定したり、各機器の最適な運用方法を決定したりするために用いられてもよい。
　（６－３）変形例Ｃ
　組合せ解決定システム１９０が決定する組合せ解は、組合せ解の適用対象の状況を改善するものであってもよい。実施形態の場合、組合せ解の適用対象は、空調対象空間である。第１系列データは、組合せ解の適用対象の状況の改善に影響を与える、状況に関する過去のデータである。状況とは、例えば、空調対象空間の温度、湿度及び二酸化炭素濃度である。過去のデータとは、例えば、熱負荷及び換気負荷である。
　（６－４）変形例Ｄ
　必須系列データ抽出部１５０は、スパース推定によって第２系列データを抽出する。しかし、必須系列データ抽出部１５０は、スパース推定以外の他の方法によって第２系列データを抽出してもよい。他の方法とは、重回帰分析、主成分分析、判別分析及びランダムフォレストが挙げられる。
　（６－５）変形例Ｅ
　必須系列データ抽出部１５０は、極値統計によって第３系列データを抽出する。しかし、必須系列データ抽出部１５０は、極値統計以外の他の方法によって第３系列データを抽出してもよい。他の方法の一例は、ある一定の頻度で所定の値を超えた条件を抽出する方法である。他の方法の別の例は、発生した値を正規分布に当てはめて、平均と分散を算出し、ある確率で所定の値を超える条件を抽出する方法である。
　―むすび―
　以上、本開示の実施形態を説明したが、特許請求の範囲に記載された本開示の趣旨及び範囲から逸脱することなく、形態や詳細の多様な変更が可能なことが理解されるであろう。
　組合せ解決定システムは、ブラックボックス最適化問題の解を効率的に導出することができる。
１１０　　　解候補生成部
１２０　　　シミュレーション部
１３０　　　評価指標算出部
１４０　　　解決定部
１５０　　　必須系列データ抽出部
１９０　　　組合せ解決定システム
特開２００６－４８４７５号公報

Claims (6)

  1. 　組合せ解の候補を生成する解候補生成部（１１０）と、
    　前記解候補生成部が生成した組合せ解の候補に関する情報と、組合せ解を評価するための系列データと、を用いてシミュレーションデータを演算するシミュレーション部（１２０）と、
    　前記シミュレーションデータに基づいて評価指標を算出する評価指標算出部（１３０）と、
    　複数の組合せ解の候補のそれぞれから前記評価指標算出部が算出した前記評価指標に基づいて、前記複数の組合せ解の候補の中から評価の高い組合せ解を決定する解決定部（１４０）と、
    　必須系列データ抽出部（１５０）と、
    を備え、
    　前記必須系列データ抽出部は、
    　　前記解候補生成部が生成したｎ個（ｎ≧２）の組合せ解の候補のうちのｍ個（ｍ≧２）の組合せ解の候補に関する情報と、前記系列データとしての第１系列データと、から前記評価指標算出部が算出した前記評価指標に基づいて、前記第１系列データの中から所定の精度で前記評価指標を算出するために必要な第２系列データを抽出し、
    　　前記第１系列データの中から、前記ｍ個の組合せ解の候補に対して所定の制約条件の適否を検証するために必要な第３系列データを抽出し、
    　　前記第２系列データと前記第３系列データとを合わせて必須系列データを取得し、
    　前記シミュレーション部は、前記必須系列データ抽出部が前記必須系列データを抽出した後に、少なくとも、前記ｎ個の組合せ解の候補のうちの前記ｍ個の組合せ解の候補以外の組合せ解の候補に関する情報と、前記必須系列データを含む系列データと、を用いて前記シミュレーションデータを演算する、
    組合せ解決定システム（１９０）。
  2. 　前記必須系列データ抽出部は、前記所定の制約条件を満たさなくなる可能性が高い前記第１系列データを、前記第３系列データとして抽出する、
    請求項１に記載の組合せ解決定システム。
  3. 　前記必須系列データ抽出部は、スパース推定によって前記第２系列データを抽出し、極値統計によって前記第３系列データを抽出する、
    請求項１又は２に記載の組合せ解決定システム。
  4. 　前記必須系列データ抽出部は、所定の精度で前記評価指標を算出するために不要なデータと、前記所定の制約条件を考慮するときに不要なデータとを、前記第１系列データから排除することによって、前記必須系列データを取得する、
    請求項１から３のいずれか１項に記載の組合せ解決定システム。
  5. 　組合せ解は、物の組合せ解、又は、方法の組合せ解である、
    請求項１から４のいずれか１項に記載の組合せ解決定システム。
  6. 　組合せ解は、組合せ解の適用対象の状況を改善するものであり、
    　前記第１系列データは、組合せ解の適用対象の前記状況の改善に影響を与える、前記状況に関する過去のデータである、
    請求項１から５のいずれか１項に記載の組合せ解決定システム。
PCT/JP2020/026247 2019-07-04 2020-07-03 組合せ解決定システム WO2021002464A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202080048321.2A CN114041031B (zh) 2019-07-04 2020-07-03 组合解确定系统
EP20834972.0A EP3996011A4 (en) 2019-07-04 2020-07-03 COMBINATORY SOLUTION DETERMINATION SYSTEM
US17/623,167 US11568108B2 (en) 2019-07-04 2020-07-03 Combinatorial solution determination system

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2019-125449 2019-07-04
JP2019125449A JP6841305B2 (ja) 2019-07-04 2019-07-04 組合せ解決定システム

Publications (1)

Publication Number Publication Date
WO2021002464A1 true WO2021002464A1 (ja) 2021-01-07

Family

ID=74101366

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/026247 WO2021002464A1 (ja) 2019-07-04 2020-07-03 組合せ解決定システム

Country Status (5)

Country Link
US (1) US11568108B2 (ja)
EP (1) EP3996011A4 (ja)
JP (1) JP6841305B2 (ja)
CN (1) CN114041031B (ja)
WO (1) WO2021002464A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024047879A1 (ja) * 2022-09-02 2024-03-07 日本電信電話株式会社 特徴量選択装置、特徴量選択方法、及びプログラム

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7277803B2 (ja) * 2021-03-26 2023-05-19 ダイキン工業株式会社 決定システム
WO2023282348A1 (ja) * 2021-07-08 2023-01-12 ダイキン工業株式会社 出力システム、出力方法、出力プログラム、及び出力プログラムを記録したコンピュータ読み取り可能な記録媒体
JP2023010617A (ja) * 2021-07-08 2023-01-20 ダイキン工業株式会社 出力システム、出力方法、出力プログラム、及び出力プログラムを記録したコンピュータ読み取り可能な記録媒体

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006048475A (ja) 2004-08-06 2006-02-16 Fuji Electric Systems Co Ltd プラント最適設計方法
JP2009014300A (ja) * 2007-07-06 2009-01-22 Mitsubishi Electric Corp 空調熱源最適制御装置
JP2010033349A (ja) * 2008-07-29 2010-02-12 Toshiba Corp 接触点算出方法およびプログラム

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1345167A1 (en) * 2002-03-12 2003-09-17 BRITISH TELECOMMUNICATIONS public limited company Method of combinatorial multimodal optimisation
CA2554580C (en) * 2004-01-21 2014-03-25 Metalogic, Inc. Constraint condition solving method, constraint condition solving device, and constraint condition solving system
US8874477B2 (en) * 2005-10-04 2014-10-28 Steven Mark Hoffberg Multifactorial optimization system and method
CN101393570B (zh) * 2007-09-19 2013-05-29 深圳达实智能股份有限公司 中央空调运行仿真系统
JP4561890B2 (ja) * 2008-07-01 2010-10-13 ダイキン工業株式会社 群管理装置および群管理システム
JP6588800B2 (ja) * 2015-10-27 2019-10-09 アズビル株式会社 熱源運転支援装置および方法
JP2017102619A (ja) * 2015-11-30 2017-06-08 オムロン株式会社 制御パラメータ調整装置、制御パラメータ調整方法、制御パラメータ調整プログラム
CN110249344A (zh) * 2017-03-15 2019-09-17 富士胶片株式会社 最佳解判定方法、最佳解判定程序及最佳解判定装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006048475A (ja) 2004-08-06 2006-02-16 Fuji Electric Systems Co Ltd プラント最適設計方法
JP2009014300A (ja) * 2007-07-06 2009-01-22 Mitsubishi Electric Corp 空調熱源最適制御装置
JP2010033349A (ja) * 2008-07-29 2010-02-12 Toshiba Corp 接触点算出方法およびプログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3996011A4

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024047879A1 (ja) * 2022-09-02 2024-03-07 日本電信電話株式会社 特徴量選択装置、特徴量選択方法、及びプログラム

Also Published As

Publication number Publication date
US20220253572A1 (en) 2022-08-11
JP2021012479A (ja) 2021-02-04
EP3996011A1 (en) 2022-05-11
JP6841305B2 (ja) 2021-03-10
CN114041031B (zh) 2022-08-09
US11568108B2 (en) 2023-01-31
EP3996011A4 (en) 2023-08-09
CN114041031A (zh) 2022-02-11

Similar Documents

Publication Publication Date Title
WO2021002464A1 (ja) 組合せ解決定システム
Nguyen et al. A review on simulation-based optimization methods applied to building performance analysis
Arnaldo et al. Multiple regression genetic programming
He et al. Performance optimization of HVAC systems with computational intelligence algorithms
CN107392368B (zh) 一种基于气象预报的办公建筑动态热负荷组合预测方法
Tao et al. Application of SVR optimized by modified simulated annealing (MSA-SVR) air conditioning load prediction model
Tang et al. Modeling and short-term prediction of HVAC system with a clustering algorithm
Yang et al. Machine learning-based prognostics for central heating and cooling plant equipment health monitoring
Niu et al. A robust optimization model for designing the building cooling source under cooling load uncertainty
Sangireddy et al. Development of a surrogate model by extracting top characteristic feature vectors for building energy prediction
Wang et al. A data-driven uncertainty quantification method for stochastic economic dispatch
GB2584614A (en) Optimising building energy use
Araújo et al. AD-Based Surrogate Models for Simulation and Optimization of Large Urban Areas
Lux et al. Novel meshes for multivariate interpolation and approximation
Lu et al. Clustering compression-based computation-efficient calibration method for digital twin modeling of HVAC system
Haaland et al. A statistical approach to thermal management of data centers under steady state and system perturbations
Biemann et al. Data centre HVAC control harnessing flexibility potential via real-time pricing cost optimisation using reinforcement learning
Tena-García et al. Implementing data reduction strategies for the optimal design of renewable energy systems
Jiang et al. Monthly electricity consumption forecasting by the fruit fly optimization algorithm enhanced Holt-Winters smoothing method
Tian et al. Dynamic prediction of building HVAC energy consumption by ensemble learning approach
US20200278130A1 (en) Operation control method, storage medium, and operation control device
CN112990603A (zh) 考虑频域分解后数据特征的空调冷负荷预测方法及系统
Suzuki et al. Power Prediction for Sustainable HPC
Ninagawa AI time series control system modelling
Goliatt et al. Gradient boosting ensembles for predicting heating and cooling loads in building design

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: 20834972

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2020834972

Country of ref document: EP

Effective date: 20220204