Energy storage unit output rapid solving method and system for time sequence production simulation

Abstract

The invention discloses an energy storage unit output rapid solving method and system for time sequence production simulation, and belongs to the technical field of power system planning. The method comprises the following steps: acquiring power data required by time sequence production simulation of a power system, and generating an annual equivalent load curve based on the power data; based on a daily equivalent load curve in the annual equivalent load curve, determining charging and discharging working positions of the energy storage unit; and in a preset period, the workloads of the charging and discharging working positions and the arrangement coefficient of the workloads are solved, and the energy storage unit output of the charging and discharging working positions is determined based on the arrangement coefficient. The method can efficiently and quickly determine the output of the energy storage unit.

Description

Technical Field

[0001] The present invention relates to the technical field of power system planning, and more particularly to a method and system for quickly calculating the output of an energy storage unit for sequential production simulation. Background Art

[0002] Power system power balance calculation refers to the production simulation calculation of the power supply and demand balance, which is a basic calculation in power system planning and system design. For new power systems with a high proportion of renewable energy access, time-series production simulation based on time-series curve data such as load and renewable energy generation can take into account the time-series characteristics of power system production and is suitable for simulating the actual dispatching and production process of the power system. With the continuous increase in the proportion of renewable energy access, the net load volatility and uncertainty faced by the power system have increased significantly, posing huge challenges to the power system supply and renewable energy consumption. In this context, energy storage units have become an important flexible regulation resource in the new power system due to their bidirectional regulation capabilities of charging and discharging.

[0003] However, the charging and discharging working conditions of the energy storage unit bring integer variables to its optimization model. Its production simulation optimization solution belongs to a large-scale mixed integer optimization problem, which is difficult to solve and time-consuming to calculate. It is difficult to meet the computational timeliness requirements of the multi-scheme analysis and planning work of power and electricity balance with a high proportion of new energy and energy storage in the context of the new power system. Summary of the Invention

[0004] To address the above problems, the present invention proposes a method for quickly calculating the output of an energy storage unit for time-series production simulation, comprising:

[0005] Acquire power data required for power system time-series production simulation, and generate an annual equivalent load curve based on the power data;

[0006] Determining charging and discharging operating positions of the energy storage unit based on the daily equivalent load curve in the annual equivalent load curve;

[0007] In a preset period, the workload and workload arrangement coefficient of the charging and discharging working positions are solved, and the output of the energy storage unit at the charging and discharging working positions is determined based on the arrangement coefficient.

[0008] Optional power data required for power system time series production simulation, including:

[0009] The curve data of the power system's load, wind power output, and photovoltaic output throughout the year;

[0010] Capacity and output characteristic data of various units;

[0011] The rated power, charging time and energy storage efficiency of the energy storage unit.

[0012] Optional expression for the annual equivalent load curve is as follows:

[0013] P Eq (t) = P Load (t)+P EX (t)-P Wind (t)-P Solar (t)

[0014] Among them, t is the time, P Eq (t) is the equivalent load power at time t, P Load (t) is the load power at time t, P EX (t) is the AC and DC power output at time t, P Wind (t), P Solar (t) are the wind power and photovoltaic power output at time t.

[0015] Optionally, the calculation formula for charging and discharging working positions is as follows:

[0016]

[0017] Among them, t is the time, u d (t) is the discharge working position at time t, P Eq (t) is the equivalent load power at time t, P ess,mid is the average value of the equivalent load curve within period T, u c (t) is the charging working position at time t.

[0018] Optional, the average value P of the equivalent load curve within the period T ess,mid The calculation formula is as follows:

[0019]

[0020] Among them, t is the time, P Eq (t) is the equivalent load power at time t;

[0021] The average value P of the equivalent load curve within the period T ess,mid , satisfying the following calculation formula:

[0022]

[0023] Among them, t is the time, P Eq (t) is the equivalent load power at time t, f(P ess,mid ) is P ess,mid A monotonically decreasing function, η ess is the efficiency of the energy storage unit.

[0024] Optionally, use the bisection method to solve P ess,mid The monotonically decreasing function f(P ess,mid ).

[0025] Optionally, the calculation formula for the arrangement coefficient is as follows:

[0026]

[0027] Among them, k seg is the arrangement coefficient, P Eq (t) is the equivalent load power at time t, P ess (t) is the output of the energy storage unit at time t, T is the period, and t is the time.

[0028] In another aspect, the present invention further proposes a system for rapidly calculating the output of an energy storage unit for sequential production simulation, comprising:

[0029] An initialization unit, configured to obtain power data required for a power system time-series production simulation and generate an annual equivalent load curve based on the power data;

[0030] a calculation unit, configured to determine a charging and discharging operating position of the energy storage unit based on a daily equivalent load curve in the annual equivalent load curve;

[0031] The solving unit is used to solve the workload and workload arrangement coefficient of the charging and discharging working positions within a preset period, and determine the output of the energy storage unit at the charging and discharging working positions based on the arrangement coefficient.

[0032] Optional power data required for power system time series production simulation, including:

[0033] The curve data of the power system's load, wind power output, and photovoltaic output throughout the year;

[0034] Capacity and output characteristic data of various units;

[0035] The rated power, charging time and energy storage efficiency of the energy storage unit.

[0036] Optional expression for the annual equivalent load curve is as follows:

[0037] P Eq (t) = P Load (t)+P EX (t)-P Wind (t)-P Solar (t)

[0038] Among them, t is the time, P Eq (t) is the equivalent load power at time t, P Load (t) is the load power at time t, PEX (t) is the AC and DC power output at time t, P Wind (t), P Solar (t) are the wind power and photovoltaic power output at time t.

[0039] Optionally, the calculation formula for charging and discharging working positions is as follows:

[0040]

[0041] Among them, t is the time, u d (t) is the discharge working position at time t, P Eq (t) is the equivalent load power at time t, P ess,mid is the average value of the equivalent load curve within period T, u c (t) is the charging working position at time t.

[0042] Optional, the average value P of the equivalent load curve within the period T ess,mid The calculation formula is as follows:

[0043]

[0044] Among them, t is the time, P Eq (t) is the equivalent load power at time t;

[0045] The average value P of the equivalent load curve within the period T ess,mid , satisfying the following calculation formula:

[0046]

[0047] Among them, t is the time, P Eq (t) is the equivalent load power at time t, f(P ess,mid ) is P ess,mid A monotonically decreasing function, η ess is the efficiency of the energy storage unit.

[0048] Optionally, use the bisection method to solve P ess,mid The monotonically decreasing function f(P ess,mid ).

[0049] Optionally, the calculation formula for the arrangement coefficient is as follows:

[0050]

[0051] Among them, k seg is the arrangement coefficient, P Eq (t) is the equivalent load power at time t, P ess (t) is the output of the energy storage unit at time t, T is the period, and t is the time.

[0052] In yet another aspect, the present invention further provides a computing device comprising: one or more processors;

[0053] a processor for executing one or more programs;

[0054] When the one or more programs are executed by the one or more processors, the above-described method is implemented.

[0055] In another aspect, the present invention further provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed, the method described above is implemented.

[0056] Compared with the prior art, the present invention has the following beneficial effects:

[0057] The present invention provides a method for rapidly calculating the output of an energy storage unit for time-series production simulation, comprising: obtaining power data required for power system time-series production simulation and generating a full-year equivalent load curve based on the power data; determining the charging and discharging operating positions of the energy storage unit based on the daily equivalent load curves within the full-year equivalent load curve; calculating the workload and workload scheduling coefficients for the charging and discharging operating positions within a preset period, and determining the output of the energy storage unit at the charging and discharging operating positions based on the scheduling coefficients. The present invention can efficiently and rapidly determine the output of the energy storage unit. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Figure 1 is a flow chart of the method of the present invention;

[0059] Figure 2 is a flow chart of an embodiment of the method of the present invention;

[0060] Figure 3 This is an equivalent load curve diagram obtained by subtracting the output of new energy from the load curve of the system in the example calculation of the method embodiment of the present invention;

[0061] Figure 4 This is an equivalent load curve diagram after the output of the energy storage unit of the system in the example of the method embodiment of the present invention is quickly arranged;

[0062] Figure 5 This is an equivalent load curve corresponding to the MIP optimization solution result of the energy storage unit output in the example system of the method embodiment of the present invention;

[0063] Figure 6 It is a structural diagram of the system of the present invention. DETAILED DESCRIPTION

[0064] Exemplary embodiments of the present invention will now be described with reference to the accompanying drawings. However, the present invention may be embodied in many different forms and is not limited to the embodiments described herein. These embodiments are provided to provide a thorough and complete disclosure of the present invention and to fully convey the scope of the present invention to those skilled in the art. The terminology used in the exemplary embodiments shown in the accompanying drawings is not intended to limit the present invention. In the accompanying drawings, identical elements are denoted by the same reference numerals.

[0065] Unless otherwise specified, the terms used herein (including technical terms) have the meanings commonly understood by those skilled in the art. In addition, it is understood that terms defined in commonly used dictionaries should be understood to have the same meanings as those in the context of the relevant fields, and should not be understood as idealized or overly formal meanings.

[0066] Example 1:

[0067] The present invention proposes a method for quickly calculating the output of energy storage units for time-series production simulation. Figure 1 As shown, including:

[0068] Step 1: Obtain power data required for power system time series production simulation, and generate an annual equivalent load curve based on the power data;

[0069] Step 2: Determine the charging and discharging operating positions of the energy storage unit based on the daily equivalent load curve in the annual equivalent load curve;

[0070] Step 3: Within a preset period, the workload and workload arrangement coefficient of the charging and discharging working positions are solved, and based on the arrangement coefficient, the output of the energy storage unit at the charging and discharging working positions is determined.

[0071] The power data required for power system time series production simulation includes:

[0072] The curve data of the power system's load, wind power output, and photovoltaic output throughout the year;

[0073] Capacity and output characteristic data of various units;

[0074] The rated power, charging time and energy storage efficiency of the energy storage unit.

[0075] The expression of the annual equivalent load curve is as follows:

[0076] P Eq (t) = P Load (t)+P EX (t)-P Wind (t)-P Solar (t)

[0077] Among them, t is the time, P Eq (t) is the equivalent load power at time t, P Load (t) is the load power at time t, P EX (t) is the AC and DC power output at time t, P Wind (t), P Solar (t) are the wind power and photovoltaic power output at time t.

[0078] The calculation formulas for the charging and discharging working positions are as follows:

[0079]

[0080] Among them, t is the time, u d (t) is the discharge working position at time t, P Eq (t) is the equivalent load power at time t, P ess,mid is the average value of the equivalent load curve within period T, u c (t) is the charging working position at time t.

[0081] Among them, the average value P of the equivalent load curve in period T ess,mid The calculation formula is as follows:

[0082]

[0083] Among them, t is the time, P Eq (t) is the equivalent load power at time t;

[0084] The average value P of the equivalent load curve within the period T ess,mid , satisfying the following calculation formula:

[0085]

[0086] Among them, t is the time, P Eq (t) is the equivalent load power at time t, f(P ess,mid ) is P ess,mid A monotonically decreasing function, η ess is the efficiency of the energy storage unit (between 0 and 1, usually 0.85).

[0087] Among them, the bisection method is used to solve P ess,mid The monotonically decreasing function f(P ess,mid ).

[0088] The calculation formula of the arrangement coefficient is as follows:

[0089]

[0090] Among them, kseg is the arrangement coefficient, P Eq (t) is the equivalent load power at time t, P ess (t) is the output of the energy storage unit at time t, T is the period, and t is the time.

[0091] The present invention will be further described below with reference to specific embodiments:

[0092] Specific case steps, such as Figure 2 Shown, including, as follows:

[0093] Step 1: Input the data required for power system time series production simulation calculation;

[0094] Step 2: Calculate and form the equivalent load curve for the whole year;

[0095] Step 3: Arrange the charging and discharging working positions of the energy storage unit according to the daily equivalent load curve;

[0096] Step 4: Solve the problem over two consecutive days, and solve the optimal working power for each charging and discharging position in chronological order from the beginning to the end. Arrange the output of the energy storage unit at the current charging and discharging position based on the optimal working power.

[0097] Step 5: Save the output arrangement result of the energy storage unit at the current charge and discharge working position, advance to the next charge and discharge working position, and perform rolling optimization until the output arrangement for the whole year is completed.

[0098] In step 1:

[0099] The data required for the power system time series production simulation calculation includes the annual load, wind power, photovoltaic output curve data, the capacity and output characteristics of various units, etc. Among them, the important parameters of the energy storage unit include the unit rated power P ess,n , Charging time T ess and energy storage efficiency η ess (between 0-1, usually 0.85). ess,n Defined as the discharge capacity of the energy storage unit, then E ess,n =η ess ·P ess,n ·T ess Optionally, the sampling duration of the curve data is 1 hour, that is, the curve data for the whole year contains 8760 time points.

[0100] In step 2:

[0101] The expression of the annual equivalent load curve is:

[0102] P Eq (t) = P Load (t)+PEX (t)-P Wind (t)-P Solar (t)

[0103] Where, t is the time from 1 to 8760 hours; P Eq is the equivalent load power; P Load is the load power; P EX P is the AC and DC output power; Wind 、P Solar are the output power of wind power and photovoltaic power respectively.

[0104] In step 3:

[0105] The basic principle of arranging the charging and discharging working positions of the energy storage unit according to the daily equivalent load curve is to arrange the energy storage unit to discharge at the peak position of the equivalent load curve and to arrange the energy storage unit to charge at the valley position of the equivalent load curve, thereby giving full play to the peak-shaving and valley-filling effect of the energy storage unit. ess,mid , then the equivalent load curve is higher than P ess,mid The position is the possible discharge position, and the corresponding discharge position mark u d =1; equivalent load curve is lower than P ess,mid The position is a possible charging position, and the corresponding charging position mark u c =1, the expression is:

[0106]

[0107] For a continuous u d (t) = 1, the interval corresponding to t is called a discharge working position; for a continuous u c (t)=1, the interval corresponding to t is called a charging working position.

[0108] Here is a method to select the middle value P ess,mid A simple method is to take the average value of the equivalent load curve in the period T corresponding to the current day as P ess,mid , the expression is:

[0109]

[0110] Considering the impact of the energy storage unit's charging and discharging efficiency, the charge capacity obtained by dividing the middle value should be greater than the discharge capacity. Therefore, the middle value should be slightly higher than the average value of the equivalent load curve to ensure that the charge capacity is greater than the discharge capacity.

[0111] The position of the intermediate value should satisfy the requirement that during the arrangement period, the total discharge energy is equal to the total charge energy multiplied by the efficiency of the energy storage unit, as shown in the following expression:

[0112]

[0113] Note that the function f(P ess,mid ) is P ess,mid The monotonically decreasing function is used to arrange the maximum and minimum values of the equivalent load curve within the cycle as P ess,mid The initial upper and lower limits are used to iteratively solve the above monotonically decreasing function f(P ess,mid ), we can get the intermediate value P ess,mid The solution.

[0114] Specifically, the bisection method is used to solve the monotonically decreasing function f(x) = 0. First, determine the initial lower limit a and upper limit b of x, satisfying f(a) > 0 and f(b) < 0. Then, take the midpoint c = (a + b) / 2 between a and b and calculate f(c). If f(c) > 0, the solution is within (c, b), and update a = c. If f(c) < 0, the solution is within (a, c), and update b = c. Then, take the midpoint c = (a + b) / 2 again and calculate f(c). Repeat the iterative calculation until the length of the interval (a, b) is less than the preset precision value. At this time, the value of c is the solution to x.

[0115] In step 4:

[0116] The output of the energy storage unit at the current charge and discharge working position is arranged according to the working power. First, the output of the energy storage unit at the discharge working position [t B ,t E ] Arrange the output of the energy storage unit, considering the rated power limit and power capacity limit of the energy storage unit, and the maximum possible discharge power E in the current discharge working section seg,max The expression is:

[0117]

[0118] If the actual working power is taken as E seg =k seg ·E seg,max , then the marker line P ess,flag Arrange the output of energy storage unit P ess The expression is as follows:

[0119] P ess (t) = min{max(P Eq (t)-P ess,flag ,0),P ess,n}t∈[t B ,t E ]

[0120] The corresponding working power needs to be E seg , the expression is as follows:

[0121]

[0122] It can be seen that f(P ess,flag ) is P ess,flag The monotonically decreasing function can be solved using the bisection method, P ess,flag The initial lower limit value is P ess,mid , the initial upper limit value is [t B ,t E ]The maximum value of the equivalent load curve. The output of the energy storage unit in the current discharge working section P ess (t) After the arrangement, the amount of electricity stored in the energy storage unit is updated to:

[0123]

[0124] Similarly, the following is a method for arranging the charging working period of the energy storage unit. B ,t E ], the maximum possible charging capacity E seg,max The expression is:

[0125]

[0126] If the actual working power is taken as E seg =k seg ·E seg,max , then the marker line P ess,flag Arrange the output of energy storage unit P ess The expression is as follows:

[0127] P ess (t) = -min{max(P ess,flag -P Eq (t),0),P ess,n}t∈[t B ,t E ]

[0128] Output P of energy storage unit ess (t) is a negative value, indicating that it is in the charging state, and the corresponding working power needs to be E seg , the expression is as follows:

[0129]

[0130] It can be seen that f(P ess,flag ) is P ess,flag The monotonically decreasing function can be solved using the bisection method, P ess,flag The initial lower limit value is [t B ,t E ]The minimum value of the equivalent load curve, the initial upper limit is P ess,midThe output of the energy storage unit in the current discharge working section is P ess (t) After the arrangement, the amount of electricity stored in the energy storage unit is updated to:

[0131]

[0132] According to the above method, the output of the energy storage unit at each charging and discharging position is arranged from front to back in chronological order. The parameter to be determined is the working power arrangement coefficient k of each section. seg A simple strategy is the greedy strategy, where the working power arrangement coefficient k seg All are set to 1, meaning that each charging and discharging location arranges charging and discharging power to the best of its ability. However, the greedy strategy has its limitations. It can easily cause excessive discharge at one discharging location and no power to be discharged at the next discharging location, or excessive charging at one charging location and no remaining capacity for charging at the next charging location.

[0133] In order to avoid the limitations of the greedy strategy, the present invention iteratively determines the optimal power arrangement coefficient k for the current charge and discharge working position. seg For other charging and discharging working positions within the next two consecutive days, the output of the energy storage unit is arranged according to the greedy strategy. For the current charging and discharging working position, k seg The upper limit is 1 and the lower limit is 0. The objective function is to minimize the variance of the energy storage unit output after the equivalent load curve is arranged within the current charge and discharge working position and the period T corresponding to the next two consecutive days. The expression is:

[0134]

[0135] The golden section algorithm is used to search for the above objective function to find the best k seg The golden section algorithm only needs to calculate the value of the objective function once in each iteration, and the optimization speed is fast. seg , which is set in the range [a, b]. In order to reuse the objective function value of the previous iteration in the next iterative calculation, the split ratio α satisfies:

[0136] (1-α) / α=α / 1

[0137] Solving this equation yields α≈0.618, which is the golden ratio, hence the name golden ratio algorithm. The specific steps of the golden ratio algorithm are as follows:

[0138] 1) Let p=a+(1-α)(ba), q=a+α(ba), and calculate f p =f(p),f q =f(q);

[0139] 2) If |ba| is less than the allowable error, the search result x = (a + b) / 2 is obtained, otherwise go to step 3);

[0140] 3) If f p <f q , then let b=q,q=p,f q =f p , p=a+(1-α)(ba), calculate f p =f(p), return to step 2), otherwise go to step 4);

[0141] 4) Let a = p, p = q, f p =f q , q=a+α(ba), calculate f q =f(q), return to step 2).

[0142] According to the optimal power arrangement coefficient k of the current charging and discharging working position seg , obtaining the energy storage unit output schedule for the current charge / discharge operating position. After saving this, proceed to the next charge / discharge operating position and use the same method to solve for its optimal power schedule coefficient. This rolling optimization continues until the energy storage unit output schedule for the entire year is complete.

[0143] In step 5:

[0144] The energy storage unit output arrangement for the whole year requires arranging the energy storage unit output arrangement results of each charging and discharging working position saved in the rolling optimization process in chronological order to obtain the energy storage unit output arrangement results for the whole year.

[0145] The planned power grid data of a certain region in a certain year is selected to verify the rationality and effectiveness of the present invention. The load curve, new energy output curve and equivalent load curve obtained by subtracting new energy output from the load in the selected region are as follows: Figure 3 shown.

[0146] In the example system, set the rated power P of the energy storage unit ess,n 100MW, charging time T ess For 4 hours, the energy storage efficiency η ess The output of the energy storage unit is arranged on the above equivalent load curve using the method for quickly arranging the output of the energy storage unit proposed by the present invention. The arrangement result is as follows: Figure 4 As shown in .

[0147] For comparison, a rigorous mixed integer programming (MIP) model is established to solve the energy storage unit output in the above example. The model expression is:

[0148]

[0149] -P ess,n ≤P ess (t)≤P ess,n

[0150] 0≤E ess (t)≤η ess P ess,n T ess

[0151] E ess (t) = E ess (t-1)-[1-u c (t)+η ess u c (t)]·P ess (t)

[0152]

[0153] By calling the optimization solver to solve the above mixed integer programming model, the energy storage output optimization solution curve is obtained as follows: Figure 5 As shown in .

[0154] observe Figure 5 and Figure 4 As is known, the present invention can obtain very similar energy storage unit output results compared with a strict mixed integer programming model, and achieves a good peak-shaving and valley-filling effect on the equivalent load curve.

[0155] Example 2:

[0156] The present invention also proposes a system 200 for quickly calculating the output of an energy storage unit for time-series production simulation, such as Figure 6 As shown, including:

[0157] The initialization unit 201 is used to obtain power data required for power system time series production simulation and generate an annual equivalent load curve based on the power data;

[0158] A calculation unit 202 is configured to determine a charging and discharging operating position of the energy storage unit based on the daily equivalent load curve in the annual equivalent load curve;

[0159] The solving unit 203 is used to solve the workload and workload arrangement coefficient of the charging and discharging working positions within a preset period, and determine the output of the energy storage unit at the charging and discharging working positions based on the arrangement coefficient.

[0160] The power data required for power system time series production simulation includes:

[0161] The curve data of the power system's load, wind power output, and photovoltaic output throughout the year;

[0162] Capacity and output characteristic data of various units;

[0163] The rated power, charging time and energy storage efficiency of the energy storage unit.

[0164] The expression of the annual equivalent load curve is as follows:

[0165] P Eq (t) = P Load (t)+P EX (t)-P Wind (t)-P Solar (t)

[0166] Among them, t is the time, P Eq (t) is the equivalent load power at time t, P Load (t) is the load power at time t, P EX (t) is the AC and DC power output at time t, P Wind (t), P Solar (t) are the wind power and photovoltaic power output at time t.

[0167] The calculation formulas for the charging and discharging working positions are as follows:

[0168]

[0169] Among them, t is the time, u d (t) is the discharge working position at time t, P Eq (t) is the equivalent load power at time t, P ess,mid is the average value of the equivalent load curve within period T, u c (t) is the charging working position at time t.

[0170] Among them, the average value P of the equivalent load curve in period T ess,mid The calculation formula is as follows:

[0171]

[0172] Among them, t is the time, P Eq (t) is the equivalent load power at time t;

[0173] The average value P of the equivalent load curve within the period T ess,mid , satisfying the following calculation formula:

[0174]

[0175] Among them, t is the time, P Eq (t) is the equivalent load power at time t, f(Pess,mid ) is P ess,mid A monotonically decreasing function, η ess is the efficiency of the energy storage unit (between 0 and 1, usually 0.85).

[0176] Among them, the bisection method is used to solve P ess,mid The monotonically decreasing function f(P ess,mid ).

[0177] The calculation formula of the arrangement coefficient is as follows:

[0178]

[0179] Among them, k seg is the arrangement coefficient, P Eq (t) is the equivalent load power at time t, P ess (t) is the output of the energy storage unit at time t, T is the period, and t is the time.

[0180] The present invention can efficiently and quickly determine the output of the energy storage unit.

[0181] Example 3:

[0182] Based on the same inventive concept, the present invention also provides a computer device, which includes a processor and a memory, wherein the memory is used to store a computer program, the computer program includes program instructions, and the processor is used to execute the program instructions stored in the computer storage medium. The processor may be a central processing unit (CPU), or may be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. It is the computing core and control core of the terminal, which is suitable for implementing one or more instructions, specifically suitable for loading and executing one or more instructions in the computer storage medium to implement the corresponding method flow or corresponding function, so as to implement the steps of the method in the above embodiment.

[0183] Example 4:

[0184] Based on the same inventive concept, the present invention also provides a storage medium, specifically a computer-readable storage medium (Memory), which is a memory device in a computer device for storing programs and data. It can be understood that the computer-readable storage medium here can include both built-in storage media in the computer device and, of course, extended storage media supported by the computer device. The computer-readable storage medium provides a storage space that 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. These instructions can be one or more computer programs (including program codes). It should be noted that the computer-readable storage medium here can be a high-speed RAM memory or a non-volatile memory, such as at least one disk memory. The processor can load and execute one or more instructions stored in the computer-readable storage medium to implement the steps of the method in the above embodiment.

[0185] It will be understood by those skilled in the art that the embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code. The solutions in the embodiments of the present invention may be implemented in various computer languages, for example, the object-oriented programming language Java and the interpreted scripting language JavaScript.

[0186] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowcharts and / or block diagrams. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0187] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0188] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 The steps for the function specified in one or more boxes.

[0189] Although the preferred embodiments of the present invention have been described, those skilled in the art may make additional changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present invention.

[0190] Obviously, those skilled in the art may make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if such changes and modifications fall within the scope of the claims and their equivalents, the present invention is intended to include such changes and modifications.

Claims

1. A method for quickly calculating the output of an energy storage unit for time-series production simulation, characterized in that: include: Acquire power data required for power system time-series production simulation, and generate an annual equivalent load curve based on the power data; Determining charging and discharging operating positions of the energy storage unit based on the daily equivalent load curve in the annual equivalent load curve; In a preset period, the workload and workload arrangement coefficient of the charging and discharging working positions are solved, and the output of the energy storage unit at the charging and discharging working positions is determined based on the arrangement coefficient.

2. The method for quickly calculating the output of an energy storage unit according to claim 1, characterized in that: The power data required for the power system time series production simulation includes: The curve data of the power system's load, wind power output, and photovoltaic output throughout the year; Capacity and output characteristic data of various units; The rated power, charging time and energy storage efficiency of the energy storage unit.

3. The method for quickly calculating the output of an energy storage unit according to claim 1, characterized in that: The expression of the annual equivalent load curve is as follows: P Eq (t)=P Load (t)+P EX (t)-P Wind (t)-P Solar (t) Among them, t is the time, P Eq (t) is the equivalent load power at time t, P Load (t) is the load power at time t, P EX (t) is the AC and DC power output at time t, P Wind (t), P Solar (t) are the wind power and photovoltaic power output at time t.

4. The method for quickly calculating the output of an energy storage unit according to claim 1, characterized in that: The calculation formulas for the charging and discharging working positions are as follows: Among them, t is the time, u d (t) is the discharge working position at time t, P Eq (t) is the equivalent load power at time t, P ess,mid is the average value of the equivalent load curve within period T, u c (t) is the charging working position at time t.

5. The method for quickly calculating the output of an energy storage unit according to claim 4, characterized in that: The average value P of the equivalent load curve within the period T ess,mid The calculation formula is as follows: Among them, t is the time, P Eq (t) is the equivalent load power at time t; The average value P of the equivalent load curve within the period T ess,mid , satisfying the following calculation formula: Among them, t is the time, P Eq (t) is the equivalent load power at time t, f(P ess,mid ) is P ess,mid A monotonically decreasing function, η ess is the efficiency of the energy storage unit.

6. The method for quickly calculating the output of an energy storage unit according to claim 5, characterized in that: Solve for P using the bisection method ess,mid The monotonically decreasing function f(P ess,mid ).

7. The method for quickly calculating the output of an energy storage unit according to claim 1, characterized in that: The calculation formula of the arrangement coefficient is as follows: Among them, k seg is the arrangement coefficient, P Eq (t) is the equivalent load power at time t, P ess (t) is the output of the energy storage unit at time t, T is the period, and t is the time.

8. A system for quickly calculating the output of an energy storage unit for time-series production simulation, characterized in that: include: An initialization unit, configured to obtain power data required for a power system time-series production simulation and generate an annual equivalent load curve based on the power data; a calculation unit, configured to determine a charging and discharging operating position of the energy storage unit based on a daily equivalent load curve in the annual equivalent load curve; The solving unit is used to solve the workload and workload arrangement coefficient of the charging and discharging working positions within a preset period, and determine the output of the energy storage unit at the charging and discharging working positions based on the arrangement coefficient.

9. The energy storage unit output rapid solution system according to claim 8, characterized in that: The power data required for the power system time series production simulation includes: The curve data of the power system's load, wind power output, and photovoltaic output throughout the year; Capacity and output characteristic data of various units; The rated power, charging time and energy storage efficiency of the energy storage unit.

10. The energy storage unit output rapid solution system according to claim 8, characterized in that: The expression of the annual equivalent load curve is as follows: P Eq (t)=P Load (t)+P EX (t)-P Wind (t)-P Solar (t) Among them, t is the time, P Eq (t) is the equivalent load power at time t, P Load (t) is the load power at time t, P EX (t) is the AC and DC power output at time t, P Wind (t), P Solar (t) are the wind power and photovoltaic power output at time t.

11. The energy storage unit output rapid solution system according to claim 8, characterized in that: The calculation formulas for the charging and discharging working positions are as follows: Among them, t is the time, u d (t) is the discharge working position at time t, P Eq (t) is the equivalent load power at time t, P ess,mid is the average value of the equivalent load curve within period T, u c (t) is the charging working position at time t.

12. The energy storage unit output rapid solution system according to claim 11, characterized in that: The average value P of the equivalent load curve within the period T ess,mid The calculation formula is as follows: Among them, t is the time, P Eq (t) is the equivalent load power at time t; The average value P of the equivalent load curve within the period T ess,mid , satisfying the following calculation formula: Among them, t is the time, P Eq (t) is the equivalent load power at time t, f(P ess,mid ) is P ess,mid A monotonically decreasing function, η ess is the efficiency of the energy storage unit.

13. The energy storage unit output rapid solution system according to claim 12, characterized in that: Solve for P using the bisection method ess,mid The monotonically decreasing function f(P ess,mid ).

14. The energy storage unit output rapid solution system according to claim 8, characterized in that: The calculation formula of the arrangement coefficient is as follows: Among them, k seg is the arrangement coefficient, P Eq (t) is the equivalent load power at time t, P ess (t) is the output of the energy storage unit at time t, T is the period, and t is the time.

15. A computer device, characterized in that: include: one or more processors; a processor for executing one or more programs; When the one or more programs are executed by the one or more processors, the method according to any one of claims 1 to 7 is implemented.

16. A computer-readable storage medium, characterized in that A computer program is stored thereon, and when the computer program is executed, the method according to any one of claims 1 to 7 is implemented.

Images