JP7059697B2 - Optimal calculation unit, optimal calculation method, and program - Google Patents

Description

The present invention relates to an optimal calculation device, an optimum calculation method, and a program.

Optimal power flow calculation is widely used when formulating an operation plan for an electric power system. In the optimum power flow calculation, the parameters related to the power system are calculated so as to match the operation index of the power system within the range where the constraint conditions are satisfied. When the constraint condition includes a constraint condition over a plurality of times (time continuity constraint condition), the amount of calculation increases. Therefore, a technique for reducing the amount of calculation has been proposed (see, for example, Non-Patent Document 1). However, even with the technique shown in Non-Patent Document 1, it is difficult to reduce the amount of calculation when the situation contrary to the time continuity constraint condition occurs frequently.
[Prior Art Document]
[Non-patent literature]
[Non-Patent Document 1] Arisa Takehara, Ikuo Kurihara, "Development of Basic Logic for Multi-section Optimal Tide Flow Calculation Research Report: T03035", Report of Central Research Institute of Electric Power, Japan, Central Research Institute of Electric Power, March 31, 2004 , Research Report T03035

In the optimum calculation device that calculates the parameters within the range of the constraint condition, it is preferable to reduce the calculation amount.

In order to solve the above problems, in the first aspect of the present invention, an optimum calculation device is provided. The optimum calculation device may be an optimum calculation device that calculates a parameter value satisfying a preset constraint condition. The optimum calculation device includes a condition acquisition unit. The condition acquisition unit may acquire the first constraint condition and the second constraint condition. The first constraint may be a preset constraint independent for each position on the axis, and the second constraint may be a constraint set between positions on the axis. The optimum calculation device may include a first calculation unit. The first calculation unit may calculate the first provisional parameter value satisfying the first constraint condition for each position of the axis. The optimum calculation device may include an interposition determination unit. The position-to-position determination unit may determine whether or not the first provisional parameter value satisfies the second constraint condition. The optimum calculation device may include a second calculation unit. The second calculation unit may calculate the first correction parameter value having the smallest difference from the first provisional parameter value within the range where the second constraint condition is satisfied. The second calculation unit may calculate the first correction parameter value when the first provisional parameter value does not satisfy the second constraint condition.

The second calculation unit may use the first provisional parameter value as the definite parameter value when the first provisional parameter value satisfies the second constraint condition.

The optimum calculation device may further include a position-by-position determination unit. The position-by-position determination unit may determine whether or not the first correction parameter value satisfies the first constraint condition. The first calculation unit may calculate the second provisional parameter value based on the first correction parameter value when the first correction parameter value does not satisfy the first constraint condition. The first calculation unit may calculate the second provisional parameter value in the numerical range corresponding to the first correction parameter value.

One of the upper limit value and the lower limit value in the numerical range may be the first correction parameter value.

The first calculation unit and the second calculation unit may repeat the calculation of the provisional parameter value and the correction parameter value until the calculated parameter value satisfies both the first constraint condition and the second constraint condition.

The first calculation unit may gradually reduce the numerical range each time the calculation of the provisional parameter value is repeated.

The first calculation unit may calculate provisional parameter values at two or more positions by parallel processing.

The preset axis may be a time axis. The preset axis may be a spatial axis.

The parameter value may include the voltage value of at least one node included in the power system.

The program may cause the computer to function as an optimal computing unit.

In the second aspect of the present invention, an optimum calculation method is provided. The optimum calculation method may calculate a parameter value that satisfies a preset constraint condition. The optimum calculation method may include a condition acquisition stage. At the condition acquisition stage, the first constraint condition and the second constraint condition may be acquired. The first constraint may be a constraint that is independent for each position on the preset axis. The second constraint may be set between positions on the axis. The optimum calculation method may include a first calculation step. In the first calculation step, the first provisional parameter value satisfying the first constraint condition may be calculated for each position of the axis. The optimum calculation method may include a position-to-position determination step. In the position-to-position determination stage, it may be determined whether or not the first provisional parameter value satisfies the second constraint condition. The optimum calculation method may include a second calculation step. In the second calculation stage, when the first provisional parameter value does not satisfy the second constraint condition, the difference from the first provisional parameter value is the smallest in the range where the second constraint condition is satisfied. The correction parameter value may be calculated.

The outline of the above invention does not list all the necessary features of the present invention. A subcombination of these feature groups can also be an invention.

This is an example of using the optimum calculation device of the first embodiment of the present invention for the optimum power flow calculation. It is a figure which shows the schematic structure of the optimum calculation apparatus 100 of 1st Embodiment. It is an example of the flowchart which shows the processing content by the optimum calculation apparatus 100 of 1st Embodiment. It is an example of the flowchart which shows the processing content by the optimum calculation apparatus 100 of 2nd Embodiment. It is a figure which shows an example of the hardware composition of the computer which concerns on one Embodiment of this invention.

Hereinafter, the present invention will be described through embodiments of the invention, but the following embodiments do not limit the invention within the scope of the claims. Also, not all combinations of features described in the embodiments are essential to the means of solving the invention.

FIG. 1 is an example of using the optimum calculation device of the first embodiment of the present invention for optimum power flow calculation. The optimum calculation device 100 calculates a parameter value that satisfies a preset constraint condition. The optimum calculation device 100 may execute the optimum power flow calculation in the power system 2. The optimum calculation device 100 is a computer. In one example. The optimum calculation device 100 includes a CPU 102, a storage unit 104, and a communication interface 106. The CPU 102 may execute the optimum power flow calculation program stored in the storage unit 104. The optimum power flow calculation program is an example of a calculation program for calculating a parameter value satisfying a preset constraint condition.

The power system 2 subject to the optimum power flow calculation by the optimum calculation device 100 may include a plurality of power systems 200a, 200b, and 200c. In this example, the plurality of power systems 200a, 200b, and 200c correspond to power plants, respectively. Further, the electric power system 2 may include the electric power system of the consumer.

The plurality of power grids 200a, 200b, and 200c have nodes 212, 214, and 216, respectively. The plurality of power grids 200a, 200b, and 200c are connected to each other by nodes 212, 214, and 216, respectively. The power system 200a may include a generator 220 and a phase adjustment facility 230. The generator 220 may be connected to a transformer 222 whose tap ratio can be changed. The generator 220 powers the load via the node 212 and the line. The phase adjustment equipment 230 adjusts the reactive power. The phase adjusting equipment 230 may include an inductance component 232 and a capacitor component 234. Other power systems 200b and 200c have similar configurations.

The operation of the power system 2 has physical or operational constraints such that the voltage values at each node 212, 214, 216 and the output of each generator 220 must be kept within an appropriate range. The optimum calculation device 100 determines an operating state for optimizing the purpose such as minimizing the fuel cost while satisfying all the various constraints. Specifically, the optimum calculation device 100 calculates the variables x and τ as parameters while satisfying all the various constraints.

The optimization problem in the optimum power flow calculation is formulated as follows.

Here, f _k (x) in Eq. (1) is an evaluation function (objective function), x is a variable, τ is a variable, and Eqs. (2) and (3) are the first constraints. Yes, Eq. (4) is the second constraint. The evaluation function f _k (x) indicates an operational index of the power system 2. The evaluation function f _k (x) may be at least one of fuel cost minimization and transmission loss minimization. The variable x may be at least one of the voltage value and the power flow (active power P, reactive power Q) of each bus. The variable τ may be at least one of the transformer tap ratio and the phase adjustment equipment input.

Equations (2) and (3), which are the first constraints, have independent constraints (sometimes referred to as "individual constraints") for each time (time cross section) k = 1, ..., K. show. The individual constraint (g (x)) with x as a variable may include at least a part of the load bus power flow equation, the upper and lower limit constraints of the bus voltage, and the upper limit constraint of the power flow flowing through the transmission line. Further, the individual constraint condition (g (τ)) with τ as a variable is the upper / lower limit constraint of the active power of the generator 220, the upper / lower limit constraint of the invalid power of the generator 220, and the upper / lower limit constraint of the tap ratio of the transformer 222. , And at least one of the upper and lower limit constraints of the input amount of the phase adjustment equipment 230 may be included.

Equation (4), which is the second constraint condition, shows a constraint condition (time continuity constraint condition) set between a plurality of times (time cross sections). The second constraint condition (h) includes at least one of the upper limit constraint of the change in the active power output of the generator 220, the upper limit constraint of the tap ratio change of the transformer 222, and the upper limit constraint of the input amount change of the phase adjustment equipment 230. good. The second constraint condition (time continuity constraint condition) may be a constraint condition related to the amount of time change.

However, the first constraint condition (g) and the second constraint condition (h) are not limited to the above cases. The first constraint condition (g) may be an independent constraint condition for each position (time) on a predetermined time axis. The second constraint condition (h) may be any constraint condition set between a plurality of positions (time) on a predetermined time axis. The space between a plurality of positions (time) may be between two positions (time) and may be between three or more positions (time). Moreover, the evaluation function and the constraint condition itself may be the same as those used in the conventional optimum power flow calculation. Therefore, detailed description will be omitted.

In the optimum calculation device 100, the variables x and τ which are the calculated parameters may be instructed to the power systems 200a, 200b, and 200c via the communication interface 106 and the network. Alternatively, the calculated parameters may be transmitted to a terminal that executes the formulation of the operation plan of the electric power system 2 and used for the formulation of the operation plan.

In the example shown in FIG. 1, as the variable τ which is a parameter, the transformer tap ratio and the phase adjustment equipment input amount in the phase adjustment equipment 230 are transmitted to the power systems 200a, 200b, and 200c. Each power system 200a, 200b, and 200c may receive a command to control the tap ratio of the transformer 222 and the input amount of the phase adjustment equipment 230. As a result, it is possible to realize the operation of the power system 2 suitable for the operation indicators such as the minimization of fuel cost and the minimization of transmission loss.

FIG. 2 is a diagram showing a schematic configuration of the optimum calculation device 100 of the first embodiment. The optimum calculation device 100 includes a condition acquisition unit 10, a first calculation unit 12, an interposition determination unit 14, a second calculation unit 16, and a position-by-position determination unit 18. Each part may be realized as a function of the CPU 102. The condition acquisition unit 10 acquires the first constraint condition and the second constraint condition. The condition acquisition unit 10 may acquire the first constraint condition and the second constraint condition from the outside via the communication interface 106 or the like. In this example, the case where the optimum calculation device 100 calculates the variable τ as a parameter will be described.

As described above, the first constraint condition (g) is an independent constraint condition for each position (time) on a predetermined time axis. The second constraint condition (h) is a constraint condition set between a plurality of positions (time) on a predetermined time axis.

The first calculation unit 12 calculates a first provisional parameter value τ'that satisfies the first constraint condition (g) for each position on the time axis. That is, the first calculation unit 12 calculates the first provisional parameter value τ'for each time on the time axis. The first calculation unit 12 calculates the first provisional parameter value τ'as a relaxation problem excluding the second constraint condition (h) from the constraint conditions.

This mitigation problem can be divided into the optimum tidal current calculation problem for each time (single time cross section) as shown in the following equations (5), (6), and (7). Therefore, the first calculation unit 12 may calculate the provisional parameter value τ'at two or more times (positions on the time axis) by parallel processing.

The position-to-position determination unit 14 determines whether or not the first provisional parameter value τ'satisfies the second constraint condition (h). In the second calculation unit 16, when the first provisional parameter value τ'does not satisfy the second constraint condition (h), the difference from the first provisional parameter value τ'is within the range satisfying the second constraint condition (h). The smallest first correction parameter value τ "is calculated. The second calculation unit 16 calculates the first correction parameter value τ" using the following equations (8), (9), and (10). You can do it. That is, the second calculation unit 16 corrects the first provisional parameter value τ', which is the solution, and executes the solution correction calculation for calculating the first correction parameter value τ ".

The second calculation unit 16 may use the first provisional parameter value τ'as a definite parameter value when the first provisional parameter value τ'satisfies the second constraint condition (h). The position-by-position determination unit 18 determines whether or not the first correction parameter value τ "satisfies the first constraint condition (g). In the first calculation unit 12, the first correction parameter value τ" is the first constraint. When the condition (g) is not satisfied, the second provisional parameter value τ'is calculated based on the first correction parameter value τ ″. Specifically, the first calculation unit 12 calculates the first correction parameter value τ ″. The second provisional parameter value τ'is calculated so as to satisfy the first constraint condition (g) in the numerical range corresponding to.

In one example, when the first correction parameter value τ ″ is calculated so that τ _k ′> τ _k ″ at the time (time cross section) k, the first calculation unit 12 is the first at this time k. The numerical range of the constraint condition is changed from the following equation (11) to equation (12). Specifically, the first calculation unit 12 may set the upper limit value of the numerical range of the first constraint condition to the first correction parameter value τ ".

On the other hand, when the first correction parameter value τ ″ is calculated so that τ _k ″> τ _k ′ in the time cross section k, the first calculation unit 12 determines the numerical value of the first constraint condition at this time k. The range is changed from the above equation (11) to the following equation (13). Specifically, the first calculation unit 12 may set the lower limit value of the numerical range of the first constraint condition to the first correction parameter value τ ".

As described above, the first calculation unit 12 changes the numerical range in the constraint condition according to the first correction parameter value τ ″. Specifically, one of the upper limit value and the lower limit value in the numerical range is the first. Then, in the numerical range changed in this way, the first calculation unit 12 calculates the second provisional parameter value τ'so as to satisfy the first constraint condition (g). Therefore, the first calculation unit 12 solves the division relaxation problem in which the constraint range is narrowed by making the above modifications.

The first calculation unit 12 and the second calculation unit 16 indicate the provisional parameter value τ'and the correction parameter value until the calculated parameter value τ satisfies both the first constraint condition (g) and the second constraint condition (h). The calculation of _τk ″ may be repeated. The first calculation unit 12 may gradually reduce the numerical range each time the calculation of the provisional parameter value is repeated.

The optimum calculation device 100 of the present embodiment configured as described above performs the following processing. FIG. 3 is an example of a flowchart showing the processing contents by the optimum calculation device 100 of the first embodiment. FIG. 3 shows an optimum calculation method of parameters in this embodiment.

The condition acquisition unit 10 acquires the first constraint condition (g) and the second constraint condition (h) (step S10). The first constraint condition (g) is an independent constraint condition for each position (time) on a predetermined time axis. The second constraint condition (h) is a constraint condition set between a plurality of positions (time) on a predetermined time axis.

The first calculation unit 12 may calculate the first provisional parameter values τ'and x'that satisfy the first constraint condition for each time using parallel processing (step S11). In this example, the case where the variable τ is calculated as a parameter will be described as an example. In step S11, the second constraint condition (h) set at a plurality of times is not included as the constraint condition. Therefore, the constraints are relaxed. Under this relaxed constraint, the optimal tidal current calculation problem can be divided into optimal tidal current calculation problems for each time (single time section). Therefore, the first calculation unit 12 may calculate the provisional parameter value τ'at two or more times (positions on the time axis) by parallel processing.

The position-to-position determination unit 14 determines whether or not the first provisional parameter value τ'satisfies the second constraint condition (h) (step S12). When the first provisional parameter value τ'satisfies the second constraint condition (h) (step S12: YES), the second calculation unit 16 may use the first provisional parameter value τ'as a definite parameter value (step). S16). Then, the optimum calculation device 100 ends the process. On the other hand, when the first provisional parameter value τ'does not satisfy the second constraint condition (h) (step S12: NO), the second calculation unit 16 is in the range where the second constraint condition (h) is satisfied. 1 Calculate the first correction parameter value τ ″ with the smallest difference from the provisional parameter value τ ′ (step S13). The first correction parameter value τ ″ with the smallest difference is calculated so that the absolute value of the difference becomes smaller. It may be calculated by doing so, or it may be calculated so that the sum of the squares of the differences is the smallest.

The first calculation unit 12 changes the numerical range in the first constraint condition (g) according to the first correction parameter value τ "(step S14). The first calculation unit 12 has an upper limit value or a lower limit value in the numerical range. The numerical range may be changed so that one of them has the first correction parameter value τ ". Specifically, when the first correction parameter value τ ″ is calculated so that τ _k ′> τ _k ″ at the time (time cross section) k, for example, the first calculation unit 12 determines this time. The numerical range of the first constraint condition in k is changed from the above equation (11) to the equation (12). The first calculation unit 12 may set the upper limit value of the numerical range of the first constraint condition as the first correction parameter value τ ".

On the other hand, when the first correction parameter value τ ″ is calculated so that τ _k ″> τ _k ′ at time k, the numerical range of the first constraint condition at this time k is calculated from the above equation (11). Change to the following equation (13). The first calculation unit 12 may set the lower limit value of the numerical range of the first constraint condition as the first correction parameter value τ ″. However, the lower limit value or the upper limit value of the numerical range of the first constraint condition may be set as the first correction parameter. Instead of using the value τ "itself, a value obtained by adding or subtracting a correction value to the first correction parameter value τ" or a value obtained by multiplying the correction coefficient may be used.

The position-by-position determination unit 18 determines whether or not the first correction parameter value τ "satisfies the first constraint condition (g) (step S15). The first correction parameter value τ" is the first constraint condition (g). ) Is not satisfied (step S15: NO), the process returns to step S11. The first calculation unit 12 calculates the second provisional parameter value τ'so as to satisfy the first constraint condition (g) in the numerical range corresponding to the first correction parameter value τ "(step S11).

In the first calculation unit 12 and the second calculation unit 16, the calculated parameter values satisfy both the first constraint condition (g) and the second constraint condition (h) (step S12: YES, and step S15: YES). ), The process (step S11 to step S15) including the calculation of the provisional parameter value (step S11) and the calculation of the correction parameter value (step S13) may be repeated. A provisional parameter or a correction parameter that satisfies both the first constraint condition (g) and the second constraint condition (h) is terminated as a deterministic parameter (step S16).

As described above, according to the optimum calculation device 100 of the present embodiment, the optimum operation considering the dynamic change of the power system in the future is realized by performing the optimum power flow calculation considering a plurality of times (time cross sections). can do. In addition, regarding the optimum power flow calculation considering a plurality of times, the calculation for each time (time cross section) can be processed in parallel. Therefore, the solution can be calculated in almost the same calculation time as the optimum power flow calculation considering a single time (time cross section).

In the above first embodiment, the case where the preset axis is the time axis has been described. That is, the parameters were calculated using the first constraint condition independent for each position (time) on the time axis and the second constraint condition set between the times on the time axis. However, the present invention is not limited to this case. The preset axis may be the spatial axis. In this case, the first constraint condition is an independent constraint condition for each position on the spatial axis. The second constraint condition is a constraint condition set between a plurality of positions on the spatial axis. For example, a position on the spatial axis is a node. The first constraint condition may be an independent constraint condition for each node. The second constraint condition may be a constraint condition set between a plurality of nodes. For example, the second constraint condition may be an upper limit constraint or the like for a change in a physical quantity such as a voltage between a plurality of nodes.

FIG. 4 is an example of a flowchart showing the processing contents by the optimum calculation device 100 of the second embodiment. The configuration of the optimum calculation device 100 of the second embodiment is the same as the configuration shown in FIG. In the optimum calculation device 100 of the present embodiment, in FIG. 2, the condition acquisition unit 10 has a first constraint condition that is independent for each position on the spatial axis set in advance, and a second constraint that is set between the positions on the spatial axis. Get constraints. The "positional interval" may be between two or more positions and may be between three or more positions. The first constraint may be a spatial continuity constraint. The second constraint may be an individual constraint.

The first calculation unit 12 calculates the first provisional parameter value satisfying the first constraint condition for each position of the space axis. The position-to-position determination unit 14 determines whether or not the first provisional parameter value satisfies the second constraint condition. The second calculation unit 16 calculates the first correction parameter value having the smallest difference from the first provisional parameter value within the range where the second constraint condition is satisfied when the first provisional parameter value does not satisfy the second constraint condition. .. As described above, the optimum calculation device 100 of the second embodiment has the same structure as the optimum calculation device 100 of the first embodiment, except that the target axis is not the time axis but the space axis. good. Therefore, the repeated description will be omitted.

The condition acquisition unit 10 acquires the first constraint condition and the second constraint condition (step S20). The first constraint condition (g) is a constraint condition independent for each position on a predetermined spatial axis. The second constraint condition (h) is a constraint condition set between a plurality of positions on a predetermined spatial axis. The position on the spatial axis may be a node. The plurality of nodes may be between two nodes, or may be between three or more nodes.

The first calculation unit 12 calculates the first provisional parameter values τ'and x'that satisfy the first constraint condition for each node by using parallel processing (step S21). In this example, the case where the variable τ is calculated as a parameter will be described as an example. In step S21, the second constraint condition set in a plurality of spatial positions is not included as the constraint condition. Therefore, the constraints are relaxed. Under this relaxed constraint, the optimal power flow calculation problem can be divided into the optimal power flow calculation problem for each single node. Therefore, the first calculation unit 12 may calculate the provisional parameter value τ'at two or more nodes (positions on the spatial axis) by parallel processing.

The position-to-position determination unit 14 determines whether or not the first provisional parameter value τ'satisfies the second constraint condition (h) (step S22). When the first provisional parameter value τ'satisfies the second constraint condition (h) (step S22: YES), the second calculation unit 16 sets the first provisional parameter value τ'as a definite parameter value (step). S26). Then, the optimum calculation device 100 ends the process. On the other hand, when the first provisional parameter value τ'does not satisfy the second constraint condition (h) (step S22: NO), the second calculation unit 16 is in the range where the second constraint condition (h) is satisfied. 1 Calculate the first correction parameter value τ ″ with the smallest difference from the provisional parameter value τ ′ (step S23). The first correction parameter value τ ″ with the smallest difference is calculated so that the absolute value of the difference becomes smaller. It may be calculated by doing so, or it may be calculated so that the sum of the squares of the differences is the smallest.

The first calculation unit 12 changes the numerical range in the first constraint condition (g) according to the first correction parameter value τ ″ (step S24). The first calculation unit 12 changes the upper limit value or the lower limit value in the numerical range. The numerical range may be changed so that one of them has the first correction parameter value τ ". Specifically, the first calculation unit 12 is, for example, in the node k, when the first correction parameter value τ ″ is calculated so that τ _k ′> τ _k ″, the first calculation unit 12 in the node k. The numerical range of the constraint condition is changed from the above equation (11) to the equation (12). The first calculation unit 12 may set the upper limit value of the numerical range of the first constraint condition as the first correction parameter value τ ".

On the other hand, when the first correction parameter value τ ″ is calculated so that τ _k ″> τ _k ′ in the node k, the numerical range of the first constraint condition in this node k is set from the above equation (11). Change to the following equation (13). The first calculation unit 12 may set the lower limit value of the numerical range of the first constraint condition as the first correction parameter value τ ″. However, the lower limit value or the upper limit value of the numerical range of the first constraint condition may be set as the first correction parameter. Instead of using the value τ "itself, a value obtained by adding or subtracting a correction value to the first correction parameter value τ" or a value obtained by multiplying the correction coefficient may be used.

The position-by-position determination unit 18 determines whether or not the first correction parameter value τ "satisfies the first constraint condition (g) (step S25). The first correction parameter value τ" is the first constraint condition (g). ) Is not satisfied (step S25: NO), the process returns to step S21. The first calculation unit 12 calculates the second provisional parameter value τ'so as to satisfy the first constraint condition (g) in the numerical range corresponding to the first correction parameter value τ "(step S21). The optimum calculation device 100 may calculate the parameter x in the same manner as in the above τ.

As described above, according to the optimum calculation device 100 of the present embodiment, even when the optimum power flow calculation considering a plurality of nodes (spatial positions) is performed, the calculation can be parallelized for each node (spatial position). .. Therefore, the solution can be calculated in almost the same calculation time as the optimum power flow calculation considering a single node.

FIG. 5 shows an example of a computer 1200 in which a plurality of aspects of the present invention can be embodied in whole or in part. The program installed on the computer 1200 causes the computer 1200 to function as an operation associated with the device according to an embodiment of the present invention or as one or more "parts" of the device, or the operation or the one or more "parts". A unit can be run and / or a computer 1200 can be run a process according to an embodiment of the invention or a step in the process. Such a program may be executed by the CPU 1212 to cause the computer 1200 to perform a specific operation associated with some or all of the blocks of the flowcharts and block diagrams described herein.

The computer 1200 according to this embodiment includes a CPU 1212, a RAM 1214, a graphic controller 1216, and a display device 1218, which are interconnected by a host controller 1210. The computer 1200 also includes input / output units such as a communication interface 1222, a hard disk drive 1224, a DVD-ROM drive 1226, and an IC card drive, which are connected to the host controller 1210 via the input / output controller 1220. The computer also includes legacy input / output units such as the ROM 1230 and keyboard 1242, which are connected to the input / output controller 1220 via the input / output chip 1240.

The CPU 1212 operates according to the programs stored in the ROM 1230 and the RAM 1214, thereby controlling each unit. The graphic controller 1216 acquires image data generated by the CPU 1212 in a frame buffer or the like provided in the RAM 1214 or the graphic controller 1216 itself, and displays the image data on the display device 1218.

The communication interface 1222 communicates with other electronic devices via the network. The hard disk drive 1224 stores programs and data used by the CPU 1212 in the computer 1200. The DVD-ROM drive 1226 reads the program or data from the DVD-ROM 1201 and provides the program or data to the hard disk drive 1224 via the RAM 1214. The IC card drive reads the program and data from the IC card and / or writes the program and data to the IC card.

The ROM 1230 internally stores a boot program or the like executed by the computer 1200 at the time of activation, and / or a program depending on the hardware of the computer 1200. The input / output chip 1240 may also connect various input / output units to the input / output controller 1220 via a parallel port, a serial port, a keyboard port, a mouse port, and the like.

The program is provided by a computer-readable storage medium such as a DVD-ROM 1201 or an IC card. The program is read from a computer-readable storage medium, installed in a hard disk drive 1224, RAM1214, or ROM1230, which is also an example of a computer-readable storage medium, and executed by the CPU 1212. The information processing described in these programs is read by the computer 1200 and provides a link between the program and the various types of hardware resources described above. The device or method may be configured to implement the operation or processing of information in accordance with the use of the computer 1200.

For example, when communication is executed between the computer 1200 and an external device, the CPU 1212 executes a communication program loaded in the RAM 1214, and performs communication processing on the communication interface 1222 based on the processing described in the communication program. You may order. Under the control of the CPU 1212, the communication interface 1222 reads and reads transmission data stored in a transmission buffer area provided in a recording medium such as a RAM 1214, a hard disk drive 1224, a DVD-ROM 1201, or an IC card. The data is transmitted to the network, or the received data received from the network is written to the reception buffer area or the like provided on the recording medium.

Further, the CPU 1212 makes the RAM 1214 read all or necessary parts of the file or the database stored in the external recording medium such as the hard disk drive 1224, the DVD-ROM drive 1226 (DVD-ROM1201), and the IC card. Various types of processing may be performed on the data on the RAM 1214. The CPU 1212 may then write back the processed data to an external recording medium.

Various types of information, such as various types of programs, data, tables, and databases, may be stored on recording media for information processing. The CPU 1212 describes various types of operations, information processing, conditional judgment, conditional branching, unconditional branching, and information retrieval described in various parts of the present disclosure with respect to the data read from the RAM 1214. Various types of processing may be performed, including / replacement, etc., and the results are written back to the RAM 1214. Further, the CPU 1212 may search for information in a file, database, or the like in the recording medium. For example, when a plurality of entries each having an attribute value of the first attribute associated with the attribute value of the second attribute are stored in the recording medium, the CPU 1212 is the first of the plurality of entries. The attribute value of the attribute of is searched for the entry that matches the specified condition, the attribute value of the second attribute stored in the entry is read, and the attribute value of the second attribute is changed to the first attribute that satisfies the predetermined condition. You may get the attribute value of the associated second attribute.

The program or software module according to the above description may be stored on a computer 1200 or in a computer-readable storage medium near the computer 1200. Further, a recording medium such as a hard disk or RAM provided in a dedicated communication network or a server system connected to the Internet can be used as a computer-readable storage medium, whereby the program can be sent to the computer 1200 via the network. offer.

Although the present invention has been described above using the embodiments, the technical scope of the present invention is not limited to the scope described in the above embodiments. It will be apparent to those skilled in the art that various changes or improvements can be made to the above embodiments. It is clear from the description of the claims that the form with such changes or improvements may be included in the technical scope of the present invention.

The order of execution of each process such as operation, procedure, step, and step in the apparatus, system, program, and method shown in the claims, specification, and drawings is particularly "before" and "prior to". It should be noted that it can be realized in any order unless the output of the previous process is used in the subsequent process. Even if the scope of claims, the specification, and the operation flow in the drawings are explained using "first", "next", etc. for convenience, it means that it is essential to carry out in this order. It's not a thing.
[Item 1]
It is an optimal calculation device that calculates parameter values that satisfy preset constraint conditions.
A first constraint condition that is independent for each position on the preset axis, and a condition acquisition unit that acquires a second constraint condition that is set between the positions on the axis.
A first calculation unit that calculates a first provisional parameter value that satisfies the first constraint condition for each position of the axis, and a first calculation unit.
A position-to-position determination unit for determining whether or not the first provisional parameter value satisfies the second constraint condition,
When the first provisional parameter value does not satisfy the second constraint condition, the first correction parameter value having the smallest difference from the first provisional parameter value within the range satisfying the second constraint condition is used. The second calculation unit to calculate and
Optimal arithmetic unit equipped with.
[Item 2]
The optimum calculation device according to item 1, wherein the second calculation unit sets the first provisional parameter value as a definite parameter value when the first provisional parameter value satisfies the second constraint condition.
[Item 3]
Further, a position-by-position determination unit for determining whether or not the first correction parameter value satisfies the first constraint condition is provided.
The first calculation unit calculates a second provisional parameter value based on the first correction parameter value when the first correction parameter value does not satisfy the first constraint condition.
The optimal computing device according to item 1 or 2.
[Item 4]
The optimum calculation device according to item 3, wherein the first calculation unit calculates the second provisional parameter value in a numerical range corresponding to the first correction parameter value.
[Item 5]
The optimal calculation device according to item 4, wherein one of the upper limit value and the lower limit value in the numerical range is the first correction parameter value.
[Item 6]
The first calculation unit and the second calculation unit calculate the provisional parameter value and the correction parameter value until the calculated parameter value satisfies both the first constraint condition and the second constraint condition. The optimal calculation device according to item 4 or 5, which is repeated.
[Item 7]
The optimum calculation device according to item 6, wherein the first calculation unit gradually reduces the numerical range each time the calculation of the provisional parameter value is repeated.
[Item 8]
The optimum calculation device according to any one of items 1 to 7, wherein the first calculation unit calculates provisional parameter values at two or more positions by parallel processing.
[Item 9]
The optimal calculation device according to any one of items 1 to 8, wherein the preset axis is a time axis or a space axis.
[Item 10]
The optimal computing device according to any one of items 1 to 9, wherein the parameter value includes a voltage value of at least one node included in the power system.
[Item 11]
It is an optimal calculation method that calculates parameter values that satisfy preset constraint conditions.
A first constraint condition that is independent for each position on the preset axis, and a condition acquisition stage that acquires a second constraint condition that is set between the positions on the axis.
The first calculation step of calculating the first provisional parameter value satisfying the first constraint condition for each position of the axis, and
A position-to-position determination step for determining whether or not the first provisional parameter value satisfies the second constraint condition, and
When the first provisional parameter value does not satisfy the second constraint condition, the first correction parameter value having the smallest difference from the first provisional parameter value within the range satisfying the second constraint condition is used. The second calculation stage to calculate and
Optimal calculation method.
[Item 12]
A program for operating a computer as the optimum calculation device according to any one of items 1 to 10.

2 ... Power system, 10 ... Condition acquisition unit, 12 ... 1st calculation unit, 14 ... Position-to-position determination unit, 16 ... 2nd calculation unit, 18 ... Position-by-position determination unit, 100 ... Optimal calculation Device, 102 ... CPU, 104 ... Storage, 106 ... Communication interface, 200 ... Power system, 212 ... Node, 214 ... Node, 216 ... Node, 220 ... Generator, 222 ... Transformer Instrument, 230 ... Phase adjustment equipment, 232 ... Inductive component, 234 ... Condenser component, 1200 ... Computer, 1201 ... DVD-ROM, 1210 ... Host controller, 1212 ... CPU, 1214 ... RAM, 1216 Graphic controller, 1218, display device, 1220, input / output controller, 1222, communication interface, 1224, hard disk drive, 1226, DVD-ROM drive, 1230, ROM, 1240, input / output chip, 1242 ... Keyboard

Claims (11)

It is an optimal calculation device that calculates parameter values that satisfy preset constraint conditions in the power system .
A condition acquisition unit that acquires a first constraint condition that is independent for each position on the time axis or the space axis set in advance, and a second constraint condition that is set between the positions on the time axis or the space axis .
A first calculation unit that calculates a first provisional parameter value that satisfies the first constraint condition for each position on the time axis or the space axis .
A position-to-position determination unit for determining whether or not the first provisional parameter value satisfies the second constraint condition,
When the first provisional parameter value does not satisfy the second constraint condition, the first correction parameter value having the smallest difference from the first provisional parameter value within the range satisfying the second constraint condition is used. An optimal calculation device including a second calculation unit for calculation. The optimum calculation device according to claim 1, wherein the second calculation unit uses the first provisional parameter value as a definite parameter value when the first provisional parameter value satisfies the second constraint condition. Further, a position-by-position determination unit for determining whether or not the first correction parameter value satisfies the first constraint condition is provided.
The first calculation unit calculates a second provisional parameter value based on the first correction parameter value when the first correction parameter value does not satisfy the first constraint condition. Or the optimum calculation device according to 2. The optimum calculation device according to claim 3, wherein the first calculation unit calculates the second provisional parameter value in a numerical range corresponding to the first correction parameter value. The optimal calculation device according to claim 4, wherein one of the upper limit value and the lower limit value in the numerical range is the first correction parameter value. The first calculation unit and the second calculation unit calculate the provisional parameter value and the correction parameter value until the calculated parameter value satisfies both the first constraint condition and the second constraint condition. The optimal computing device according to claim 4 or 5. The optimum calculation device according to claim 6, wherein the first calculation unit gradually reduces the numerical range each time the calculation of the provisional parameter value is repeated. The optimum calculation device according to any one of claims 1 to 7, wherein the first calculation unit calculates provisional parameter values at two or more positions by parallel processing. The optimum calculation device according to any one of claims 1 to 8 , wherein the parameter value includes a voltage value of at least one node included in the power system. It is an optimal calculation method for calculating parameter values that satisfy preset constraint conditions in the power system .
A condition acquisition step for acquiring a first constraint condition that is independent for each position on the time axis or the space axis set in advance, and a second constraint condition that is set between the positions on the time axis or the space axis .
The first calculation step of calculating the first provisional parameter value satisfying the first constraint condition for each position of the time axis or the space axis , and
A position-to-position determination step for determining whether or not the first provisional parameter value satisfies the second constraint condition, and
When the first provisional parameter value does not satisfy the second constraint condition, the first correction parameter value having the smallest difference from the first provisional parameter value within the range satisfying the second constraint condition is used. Optimal calculation method including a second calculation stage for calculation. A program for operating a computer as the optimum calculation device according to any one of claims 1 to 9 .

Images