New energy power flow calculation method, system, equipment and medium

By establishing mathematical calculation models and calculating the power flow in each period, the problem that the existing power system cannot optimize the power flow is solved, and the efficient operation of the power system and the economic utilization of new energy are achieved.

CN119787371APending Publication Date: 2025-04-08NORTHWEST ELECTRIC POWER DESIGN INST OF CHINA POWER ENG CONSULTING GRP
View PDF 0 Cites 0 Cited by

Patent Information

Application Number
CN202411969101.X
Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2024-12-30
Publication Date
2025-04-08

AI Technical Summary

Technical Problem

The existing power system cannot accurately calculate the power flow for each period, resulting in the inability to optimize the system, increasing the cost of grid investment and low operating efficiency.

Method used

By obtaining the new energy installed capacity, new energy output characteristics and equivalent consumption capacity data of each subsystem, a mathematical calculation model is established, and the power flow for each period is calculated according to the constraints, and the optimal power flow under the economical optimal situation is determined.

Benefits of technology

It realizes accurate calculation of the power flow of the power system in each period, optimizes the system operation, reduces the cost of power grid investment, and improves the operating efficiency and new energy utilization efficiency.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119787371A_ABST
    Figure CN119787371A_ABST
Patent Text Reader

Abstract

The invention relates to the technical field of power system planning, in particular to a new energy power flow calculation method, system and device and a medium, and the method comprises the steps: obtaining the new energy installed capacity, new energy output characteristics and equivalent consumption capability data of each subsystem in a researched energy system; establishing a mathematical calculation model according to the new energy installed capacity, the new energy output data and the equivalent consumption capability data of each subsystem, and considering constraint conditions of the mathematical calculation model; and starting from the first time period, calculating the mathematical calculation model according to the constraint condition, and obtaining the power flow corresponding to each time period. According to the method, the power flow can be effectively calculated, and reference is provided for planning and design.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to the technical field of power system planning, and in particular to a new energy power flow calculation method, system, equipment and medium. Background Art

[0002] Under the background of "dual carbon", new energy sources such as wind power and photovoltaics have been vigorously developed. New energy sources in different regions have natural complementarity, and it is necessary to exchange power by adding interconnection lines between different systems to achieve mutual assistance. Traditional methods often use a large number of production simulation calculations to indirectly judge the system section demand based on the power abandonment rate indicators in the calculation results of different schemes. On the other hand, there is a peak distribution of new energy power abandonment, and there are few moments when the exchange power is large. At this time, it is not economical to strengthen the grid. Moreover, the traditional method calculates the power flow demand under the premise of ignoring the economy, which may lead to inefficient operation of the power system and increase unnecessary power grid investment costs. Therefore, a new method is needed to accurately calculate the power flow demand under the condition of new energy consumption, while taking into account economic factors to ensure the maximum cost-effectiveness of power grid operation. Summary of the invention

[0003] The technical problem to be solved by the present invention is to provide a new energy power flow calculation method, system, equipment and medium to address the deficiencies in the above-mentioned prior art, so as to solve the technical problem that the existing power system cannot calculate the power flow in each time period and thus cannot optimize the system.

[0004] The objective of the present invention is achieved by the following technical solutions: In a first aspect, the present invention provides a method for calculating a new energy power flow, comprising: Obtain data on the installed capacity of renewable energy, renewable energy output characteristics, and equivalent absorption capacity of each subsystem in the energy system under study; Establish a mathematical calculation model based on the installed capacity of new energy, new energy output data and equivalent absorption capacity data of each subsystem, and consider the constraints of the mathematical calculation model; Starting from the first period, the mathematical calculation model is calculated according to the constraints, and the power flow corresponding to each period is obtained. As a further improvement of the present invention, the mathematical calculation model is calculated based on the new energy output data and equivalent absorption capacity data of the two mutual aid subsystems in the power system, and the mathematical calculation model is:

[0005]

[0006]

[0007] In the formula, is the power system output sequence in the i-th period; is the output sequence corresponding to subsystem A; is the output sequence corresponding to subsystem B; is the equivalent absorption capacity of subsystem A; is the equivalent absorption capacity of subsystem B; is the total equivalent absorption capacity, is the system curtailment rate of the power system.

[0008] As a further improvement of the present invention, in the mathematical calculation model, the power abandonment rate of the power system for:

[0009] in, For the The power flow into system A at all times; is the equivalent absorption capacity of subsystem A; is the equivalent absorption capacity of subsystem B; is the output sequence corresponding to subsystem A; is the output sequence corresponding to subsystem B.

[0010] As a further improvement of the present invention, by calculating the power abandonment rate in the power system, the exchange power flowing into the A subsystem at each moment is calculated, and the exchange power calculation formula is:

[0011] In the formula, For the The switching power flowing into system A at all times; is the equivalent absorption capacity of subsystem A; is the equivalent absorption capacity of subsystem B; is the output sequence corresponding to subsystem A; is the output sequence corresponding to subsystem B.

[0012] As a further improvement of the present invention, after obtaining the exchange power at each moment according to the mathematical calculation model, the power flow demand at each moment is obtained according to the exchange power, and the optimized power flow under the economic optimum is determined.

[0013] As a further improvement of the present invention, the corresponding optimized power flow under the economic optimum is:

[0014]

[0015] In the formula, is the number of economically utilized hours, is the exchange power between systems A and B at the i-th moment; The corresponding power flow under the economic optimum condition; It represents the sum of the hours when the exchange power is greater than the economic optimum.

[0016] As a further improvement of the present invention, the economic utilization hours are calculated based on the grid investment, a preset annual cost coefficient and the new energy cost per kilowatt-hour of the power system.

[0017] In a second aspect, the present invention provides a new energy power flow calculation system, which is used to implement the above-mentioned new energy power flow calculation method, including: The data acquisition module is used to obtain the installed capacity of new energy, the output characteristics of new energy and the equivalent absorption capacity data of each subsystem in the energy system under study; The model building module is used to establish a mathematical calculation model based on the installed capacity of new energy, new energy output data and equivalent absorption capacity data of each subsystem, taking into account the constraints of the mathematical calculation model; The power flow calculation module calculates the mathematical calculation model according to the constraint conditions starting from the first time period to obtain the power flow corresponding to each time period.

[0018] In a third aspect, the present invention provides a computer-readable storage medium storing one or more programs, wherein the one or more programs include instructions, which, when executed by a computing device, enable the computing device to execute the above-mentioned new energy power flow calculation method.

[0019] In a fourth aspect, the present invention provides a computing device, comprising: One or more processors, a memory and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs include steps for executing the above-mentioned new energy power flow calculation method.

[0020] The beneficial effect of the present invention is that by obtaining the installed capacity, output data and equivalent absorption capacity data of each subsystem, the distribution, scale and absorption capacity of the new energy resources in the entire power system can be clearly and comprehensively understood. Based on these detailed data and the mathematical calculation model established based on them, scientific predictions can be made on the degree of expansion and absorption capacity improvement direction required by the future power system when accepting more new energy, so as to formulate a practical and forward-looking system upgrade and transformation plan to ensure that the power system can continuously and stably accept the incorporation of new energy. Based on the established mathematical calculation model, the power flow corresponding to each time period is calculated, so that the power dispatcher can reasonably arrange the power exchange between different subsystems in each time period according to the accurate calculation results, so that the electric energy can achieve efficient flow and optimal distribution in the system, minimize the problems of power waste and congestion caused by unreasonable scheduling, and improve the efficiency and stability of the operation of the entire power system.

[0021] By accurately grasping the equivalent absorption capacity of each subsystem and rationally allocating power flow in each time period, each subsystem can be encouraged to absorb the output of renewable energy as fully as possible, reduce the phenomenon of power abandonment caused by insufficient absorption capacity, improve the utilization efficiency of renewable energy in the entire power system, give full play to the role of renewable energy in energy supply, and promote the optimization of energy structure towards clean and low-carbon direction. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0023] Figure 1 is a schematic diagram of a method for calculating a new energy power flow in an embodiment of the present invention; Figure 2 is a schematic diagram of determining economic power flow according to an embodiment of the present invention; Figure 3 It is a schematic diagram of the structure of an electronic device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0024] In order to make the purpose and technical solution of the present invention clearer and easier to understand, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. The specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0025] The technical solution of the present invention will be clearly and completely described below in conjunction with the accompanying drawings and specific embodiments, wherein the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments.

[0026] Example 1 like Figure 1-Figure 2 As shown, this embodiment provides a new energy power flow calculation method. By obtaining the installed capacity of new energy, the output characteristics of new energy and the equivalent absorption capacity data of each subsystem in the energy system under study; establishing a mathematical calculation model based on the installed capacity of new energy, the output characteristics of new energy and the equivalent absorption capacity data of each subsystem; starting from the first period, the mathematical calculation model is calculated according to the constraints of the mathematical calculation model, and the corresponding power flow of each period is obtained. This method is convenient and fast to calculate. The following is a specific implementation method.

[0027] Collect data such as the installed capacity of new energy, output characteristics of new energy, and equivalent absorption capacity of each subsystem in the system under study. Among them, the method of this embodiment measures two provincial power grid subsystems. The output of new energy at any time is the product of the installed capacity of new energy and the output characteristics of new energy at that time. The output sequence of new energy is shown in Table 1. The equivalent absorption capacities of subsystems A and B are 37.55 million kilowatts and 28.83 million kilowatts respectively.

[0028] Table 1. A and B new energy output series Unit: 10,000 kilowatts

[0029] As shown in Table 1, the unit kilowatt investment of the power grid is 1,000 yuan / kilowatt, the annual cost coefficient is 10%, and the cost of new energy electricity is 0.2 yuan / kWh.

[0030] A mathematical calculation model is established based on the new energy installed capacity, new energy output data and equivalent absorption capacity data of each subsystem, and the constraints of the mathematical calculation model are taken into consideration.

[0031] The mathematical calculation model is:

[0032]

[0033] In the formula, is the power system output sequence in the i-th period; is the output sequence corresponding to subsystem A; is the output sequence corresponding to subsystem B; is the equivalent absorption capacity of subsystem A; is the equivalent absorption capacity of subsystem B; is the total equivalent absorption capacity.

[0034] The system curtailment rate of the power system for:

[0035] in, Indicates the smaller value of the two.

[0036] Then, according to the definition, we can get:

[0037] In the formula, For the The power flow into system A at all times.

[0038] At any time, by calculating the power abandonment rate in the power system, the exchange power flowing into subsystem A at each time is calculated. The exchange power calculation formula is:

[0039] In the formula, For the The switching power flowing into system A at all times; is the equivalent absorption capacity of subsystem A; is the equivalent absorption capacity of subsystem B; is the output sequence corresponding to subsystem A; is the output sequence corresponding to subsystem B.

[0040] Specifically, when and At this moment, each subsystem can absorb its own new energy. ; when and At this moment, each subsystem cannot absorb its own new energy. ; when and At this moment, system A still has the ability to absorb the traffic and can assist system B in absorbing the traffic. At this time, there are two situations: i) Then there is

[0041] ii) Then there is

[0042] when and At this moment, system B still has the ability to absorb and can assist system A in absorbing. Is a negative value.

[0043] i) Then there is

[0044] ii) Then there is

[0045] According to the above formula, solve the exchange power at each moment; Starting from the first period, the output of system A 1088 is less than 3755, and the output of system B 1378 is less than 2883. At this moment, the power flow between systems A and B is 0.

[0046] In the second period, the output of system A 1126 is less than 3755, and the output of system B 1235 is less than 2883. At this moment, the power flow between systems A and B is 0.

[0047] In the 3rd to 10th period, the power flow between the AB system is 0 at this moment.

[0048] In the 11th period, the output of system A is 3734, which is less than 3755, and the output of system B is 3145, which is greater than 2883. The output of system A+B is greater than the equivalent absorption capacity of the combined system. According to the above formula, the power flow between system A and system B at this moment is 3755-3734=210,000 kilowatts.

[0049] And so on until the last time period is reached. This embodiment mainly targets the power flow within the 8760 time period.

[0050] After obtaining the exchange power at each moment according to the mathematical calculation model, the power flow demand at each moment is obtained according to the exchange power, and the optimized power flow under the economic optimal condition is determined.

[0051] In this embodiment, the power flow considering economic considerations is obtained in the following manner: The investment in the power grid should be less than or equal to its benefits, so it should meet the following conditions: Total power grid investment*annual cost coefficient Abandoned electricity* cost of new energy If the per unit value is used, the following conditions must be met: Grid unit kilowatt investment * annual cost coefficient Utilization hours*renewable energy cost per kilowatt-hour When the equal sign is taken, it is the economic utilization hours of power flow considering economy. Economic utilization hours = grid unit kilowatt investment * annual cost coefficient / new energy kilowatt-hour cost The economic power flow can be obtained as follows:

[0052]

[0053] In the formula, is the number of economically utilized hours, is the exchange power between systems A and B at the i-th moment; The corresponding power flow under the economic optimum condition; It represents the sum of the hours when the exchange power is greater than the economic optimum.

[0054] Finally, the power flow between the two subsystems in this embodiment is obtained as: =11.9 million kilowatts In summary, the power flow corresponding to each time period is calculated based on the established mathematical calculation model, so that the power dispatcher can reasonably arrange the power exchange between different subsystems in each time period according to the accurate calculation results, so that the electric energy can achieve efficient flow and optimal distribution within the system, minimize the problems of power waste and congestion caused by unreasonable scheduling, and improve the efficiency and stability of the entire power system. By accurately grasping the equivalent absorption capacity of each subsystem and the reasonable allocation of power flow in each time period, each subsystem can be prompted to absorb the output of new energy as fully as possible, reduce the phenomenon of power abandonment caused by insufficient absorption capacity, etc., improve the utilization efficiency of new energy in the entire power system, give full play to the role of new energy in energy supply, and promote the optimization of energy structure towards clean and low-carbon.

[0055] Optimized power flow distribution and efficient new energy consumption can reduce economic losses caused by power abandonment, unreasonable scheduling, etc., while also avoiding additional costs such as excessive construction of redundant facilities to cope with new energy consumption issues, thereby reducing the overall operating costs of the power system and improving economic benefits.

[0056] Example 2 This embodiment further provides a new energy power flow calculation system, which is used to implement the new energy power flow calculation method in Embodiment 1. The system includes: The data acquisition module is used to collect data on the installed capacity of new energy, output characteristics of new energy and equivalent absorption capacity of each subsystem in the energy system under study; The model building module is used to establish a mathematical calculation model based on the installed capacity of new energy, new energy output data and equivalent absorption capacity data of each subsystem, and determine the constraints of the mathematical calculation model; The power flow calculation module calculates the mathematical calculation model according to the constraint conditions starting from the first time period to obtain the power flow corresponding to each time period.

[0057] The specific implementation methods of each module have been described in Example 1 and will not be repeated here.

[0058] Example 3 In another embodiment of the present invention, the present invention also provides a storage medium, specifically a computer-readable storage medium (Memory), which is a memory device in a terminal device for storing programs and data. It is understandable that the computer-readable storage medium here can include both the built-in storage medium in the terminal device and the extended storage medium supported by the terminal device, and can be any tangible medium containing or storing a program, which can be used by an instruction execution system, device or device or used in combination with it. The computer-readable storage medium provides a storage space, which stores the operating system of the terminal. In addition, one or more instructions suitable for being loaded and executed by the processor are also stored in the storage space, and these instructions can be one or more computer programs (including program codes). It should be noted that more specific examples (non-exhaustive list) of the computer-readable storage medium here include: an electrical connection with one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above.

[0059] Computer readable storage media also include data signals propagated in baseband or as part of a carrier wave, which carry readable program codes. Such propagated data signals can take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. The readable storage medium can also be any readable medium other than a readable storage medium, which can send, propagate, or transmit a program for use by or in conjunction with an instruction execution system, device, or device. The program code contained on the readable storage medium can be transmitted using any appropriate medium, including but not limited to wireless, wired, optical cable, RF, etc., or any suitable combination of the above.

[0060] Program code for performing the operations of the present invention may be written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Java, C++, etc., and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user computing device, partially on the user device, as a separate software package, partially on the user computing device and partially on a remote computing device, or entirely on a remote computing device or server. In the case of a remote computing device, the remote computing device may be connected to the user computing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computing device (e.g., through the Internet using an Internet service provider).

[0061] The processor may load and execute one or more instructions stored in a computer-readable storage medium to implement the corresponding steps of the method for calculating the new energy power flow in the above embodiment; the processor may load and execute the following steps: Obtain data on the installed capacity of renewable energy, renewable energy output characteristics, and equivalent absorption capacity of each subsystem in the energy system under study; Establish a mathematical calculation model based on the installed capacity of new energy, new energy output data and equivalent absorption capacity data of each subsystem, and consider the constraints of the mathematical calculation model; Starting from the first period, the mathematical calculation model is calculated according to the constraints, and the power flow corresponding to each period is obtained.

[0062] Example 4 See also Figure 3 , the terminal device is a computer device, and the computer device 60 of this embodiment includes: a processor 61, a memory 62, and a computer program 63 stored in the memory 62 and executable on the processor 61. When the computer program 63 is executed by the processor 61, the new energy power flow calculation method in the embodiment is implemented. To avoid repetition, it is not described one by one here. Alternatively, when the computer program 63 is executed by the processor 61, the functions of each model / unit in the calculation system in the new energy power flow calculation method in Embodiment 1 are implemented. To avoid repetition, it is not described one by one here.

[0063] The computer device 60 may be a computing device such as a desktop computer, a notebook, a PDA, or a cloud server. The computer device 60 may include, but is not limited to, a processor 61 and a memory 62. Those skilled in the art will appreciate that Figure 3This is only an example of the computer device 60 and does not constitute a limitation of the computer device 60. It may include more or fewer components than shown in the figure, or a combination of certain components, or different components. For example, the computer device may also include input and output devices, network access devices, buses, etc.

[0064] The processor 61 may be a central processing unit (CPU), or other general-purpose processors, central processing units, graphics processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, data processing logic based on quantum computing, discrete hardware components, etc. A general-purpose processor may be a microprocessor or the processor may be any conventional processor, etc.

[0065] The memory 62 may be an internal storage unit of the computer device 60, such as a hard disk or memory of the computer device 60. The memory 62 may also be an external storage device of the computer device 60, such as a plug-in hard disk, a smart memory card (SmartMediaCard, SMC), a secure digital (SecureDigital, SD) card, a flash card (FlashCard), etc. equipped on the computer device 60.

[0066] Furthermore, the memory 62 may include both an internal storage unit of the computer device 60 and an external storage device. The memory 62 is used to store computer programs and other programs and data required by the computer device. The memory 62 may also be used to temporarily store data that has been output or is to be output.

[0067] Any reference to a memory, database or other medium used in the embodiments provided in this application may include at least one of a non-volatile and a volatile memory. Non-volatile memory may include a read-only memory (ROM), a magnetic tape, a floppy disk, a flash memory, an optical memory, a high-density embedded non-volatile memory, a resistive random access memory (ReRAM), a magnetic random access memory (MRAM), a ferroelectric random access memory (FRAM), a phase change memory (PCM), a graphene memory, etc. Volatile memory may include a random access memory (RAM) or an external cache memory, etc. As an illustration and not limitation, RAM may be in various forms, such as a static random access memory (SRAM) or a dynamic random access memory (DRAM), etc.

[0068] The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., but is not limited thereto. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but is not limited thereto.

Claims

1. A method for calculating the flow of renewable energy power, characterized in that: include: Obtain data on the installed capacity of renewable energy, renewable energy output characteristics, and equivalent absorption capacity of each subsystem in the energy system under study; Establish a mathematical calculation model based on the installed capacity of new energy, new energy output data and equivalent absorption capacity data of each subsystem, and consider the constraints of the mathematical calculation model; Starting from the first period, the mathematical calculation model is calculated according to the constraints, and the power flow corresponding to each period is obtained.

2. The method for calculating the new energy power flow according to claim 1, characterized in that: The mathematical calculation model is calculated based on the new energy output data and equivalent absorption capacity data of the two mutual aid subsystems in the power system. The mathematical calculation model is: In the formula, is the power system output sequence in the i-th period; is the output sequence corresponding to subsystem A; is the output sequence corresponding to subsystem B; is the equivalent absorption capacity of subsystem A; is the equivalent absorption capacity of subsystem B; is the total equivalent absorption capacity, is the system curtailment rate of the power system.

3. The method for calculating the new energy power flow according to claim 2, characterized in that: In the mathematical calculation model, the power abandonment rate of the power system for: in, For the The power flow into system A at all times; is the equivalent absorption capacity of subsystem A; is the equivalent absorption capacity of subsystem B; is the output sequence corresponding to subsystem A; is the output sequence corresponding to subsystem B.

4. The method for calculating the new energy power flow according to claim 3, characterized in that: By calculating the power abandonment rate in the power system, the exchange power flowing into subsystem A at each moment can be calculated. The exchange power calculation formula is: In the formula, For the The switching power flowing into system A at all times; is the equivalent absorption capacity of subsystem A; is the equivalent absorption capacity of subsystem B; is the output sequence corresponding to subsystem A; is the output sequence corresponding to subsystem B.

5. The method for calculating the new energy power flow according to claim 4, characterized in that: After obtaining the exchange power at each moment according to the mathematical calculation model, the power flow demand at each moment is obtained according to the exchange power, and the optimized power flow under the economic optimal condition is determined.

6. The method for calculating the new energy power flow according to claim 5, characterized in that: Under the economically optimal condition, the corresponding optimal power flow is: In the formula, is the number of economically utilized hours, is the exchange power between systems A and B at the i-th moment; The corresponding power flow under the economic optimum condition; It represents the sum of the hours when the exchange power is greater than the economic optimum.

7. The method for calculating the new energy power flow according to claim 6, characterized in that: The economic utilization hours are calculated based on the grid investment, a preset annual cost coefficient, and the new energy cost per kilowatt-hour of the power system.

8. A new energy power flow calculation system, used to implement the new energy power flow calculation method according to any one of claims 1 to 7, characterized in that: include: The data acquisition module is used to obtain the installed capacity of new energy, the output characteristics of new energy and the equivalent absorption capacity data of each subsystem in the energy system under study; Mathematical model module, used to establish mathematical calculation model according to the installed capacity of new energy, new energy output data and equivalent absorption capacity data of each subsystem, and consider the constraints of the mathematical calculation model; The power flow calculation module calculates the mathematical calculation model according to the constraint conditions starting from the first time period to obtain the power flow corresponding to each time period.

9. A computer-readable storage medium storing one or more programs, characterized in that: The one or more programs include instructions, which, when executed by a computing device, enable the computing device to execute the new energy power flow calculation method according to any one of claims 1 to 7.

10. A computing device, characterized in that include: One or more processors, a memory and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs include steps for executing the new energy power flow calculation method described in any one of claims 1 to 7.