Operation simulation method, device and equipment of integrated energy base

Abstract

The invention discloses an operation simulation method, device and equipment of an integrated energy base, and relates to the technical field of planning design and operation scheduling of the integrated energy base. In the method, on one hand, a hydroelectric output function of a hydroelectric power station is divided into a plurality of linear equations, so that a quadratic function problem can be converted into a linear function problem, and the linear equation is divided into a plurality of linear equations; and furthermore, in the operation process of simulating the integrated energy base based on the hydroelectric output function, the solving difficulty can be greatly reduced, and the simulation difficulty of the integrated energy base can be further reduced. On the other hand, the output weight of the hydroelectric power plant is determined based on the output fluctuation coefficient and the power abandoning coefficient of the wind and light power plant, and the generating capacity of the hydroelectric power plant is corrected based on the output weight. Therefore, the finally obtained generating capacity of the hydroelectric power station can be relatively accurate, and the problem that the hydroelectric power station has poor compensation and regulation performance on the wind-solar power station can be solved.

Description

Technical Field

[0001] This application relates to the field of integrated energy base planning, design and operation scheduling technology, and in particular to integrated energy base operation simulation methods, devices and equipment. Background Technology

[0002] Integrated energy bases typically include power plants such as hydroelectric, wind, and solar power stations. During the design or operation phase of an integrated energy base, its operation is usually simulated (i.e., the operation of the power plants themselves). This is done to determine the rationality of the integrated energy base design or to determine its dispatching strategy.

[0003] Currently, the hydropower output function of a hydroelectric power station is used in the operation simulation of an integrated energy base. Since the hydropower output function is a quadratic function of head and power generation flow, solving the simulation model of an integrated energy base using current solvers presents practical problems such as low solution efficiency or even the inability to find a feasible solution, resulting in high simulation difficulty. Summary of the Invention

[0004] This application provides a method, device, electronic equipment, computer-readable storage medium, and computer program product for simulating the operation of an integrated energy base, in order to at least solve the problem of high simulation difficulty in related technologies.

[0005] This application provides a method for simulating the operation of an integrated energy base, which includes a hydroelectric power station and a wind and solar power station. The method includes: Acquire hydropower and wind and solar power information. The hydropower information includes the operation constraints of the hydropower station, maximum reservoir capacity, minimum reservoir capacity, maximum power generation flow, and minimum power generation flow. The wind and solar power information includes the wind and solar power generation sequence obtained with a first duration as the statistical dimension and the maximum power transmission capacity supported by the transmission channel. In the interval between the maximum reservoir capacity and the minimum reservoir capacity, multiple discrete reservoir capacities are divided, and in the interval between the maximum power generation flow and the minimum power generation flow, multiple discrete power generation flows are divided, and based on the discrete reservoir capacities and the discrete power generation flows, multiple discrete power generation capacities of the hydropower station are determined. Based on the discrete reservoir capacity, the discrete power generation flow rate, and the discrete power generation output, the hydropower output function of the hydropower station is divided into multiple linear equations; Using the second duration as a statistical dimension, the wind and solar power generation sequence is statistically analyzed to obtain the power output fluctuation coefficient of the wind and solar power station. The second duration is greater than or equal to the first duration. Using the second duration as the statistical dimension and the maximum transmission capacity supported by the transmission channel as the constraint, the wind and solar power generation sequence is statistically analyzed to obtain the curtailment coefficient of the wind and solar power station. Based on the power output fluctuation coefficient and the curtailment coefficient of the wind and solar power station, the power output weight of the hydropower station is determined. Based on the linear equation and the operating constraint information of the hydropower station, the operation of the hydropower station is simulated, and the power generation of the hydropower station is corrected based on the output weight. Based on the wind and solar power generation information, the operation of the wind and solar power station is simulated.

[0006] This application also provides an operation simulation device for an integrated energy base, the integrated energy base including a hydroelectric power station and a wind and solar power station, the device comprising: The information acquisition module is used to acquire hydropower generation information and wind and solar power generation information. The hydropower generation information includes the operation constraint information of the hydropower station, the maximum reservoir capacity, the minimum reservoir capacity, the maximum power generation flow, and the minimum power generation flow. The wind and solar power generation information includes the wind and solar power generation sequence obtained with the first duration as the statistical dimension and the maximum power transmission capacity supported by the transmission channel. The interval division module is used to divide the interval between the maximum reservoir capacity and the minimum reservoir capacity into multiple discrete reservoir capacities, and to divide the interval between the maximum power generation flow and the minimum power generation flow into multiple discrete power generation flow rates, and to determine multiple discrete power generation capacities of the hydropower station based on the discrete reservoir capacities and the discrete power generation flow rates. An equation construction module is used to divide the hydropower output function of the hydropower station into multiple linear equations based on the discrete reservoir capacity, the discrete power generation flow rate, and the discrete power generation output. The fluctuation coefficient determination module is used to statistically analyze the wind and solar power generation sequence using the second duration as a statistical dimension to obtain the power output fluctuation coefficient of the wind and solar power station, wherein the second duration is greater than or equal to the first duration. The curtailment coefficient determination module is used to statistically analyze the wind and solar power generation sequence using the second duration as a statistical dimension and the maximum power transmission capacity supported by the transmission channel as a constraint condition, so as to obtain the curtailment coefficient of the wind and solar power station. The weight determination module is used to determine the output weight of the hydropower station based on the output fluctuation coefficient and the curtailment coefficient of the wind and solar power station. The simulation module is used to simulate the operation of the hydropower station based on the linear equation and the operating constraint information of the hydropower station, and to correct the power generation of the hydropower station based on the output weight, and to simulate the operation of the wind and solar power station based on the wind and solar power generation information.

[0007] This application also provides an electronic device, including: a memory for storing a computer program; and a processor for executing the computer program to implement the steps of the above-described integrated energy base operation simulation method.

[0008] This application also provides a computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor, implements the steps of the above-described integrated energy base operation simulation method.

[0009] In some embodiments of this application, on the one hand, by decomposing the hydropower output function of a hydropower station into multiple linear equations, a quadratic function problem can be transformed into a linear function problem. This significantly reduces the difficulty of solving the problem during the simulation of the operation of a comprehensive energy base based on the hydropower output function, thereby reducing the overall simulation difficulty. On the other hand, based on the output fluctuation coefficient and curtailment coefficient of wind and solar power stations, the output weight of the hydropower station is determined, and the power generation of the hydropower station is corrected based on this output weight. Thus, the final power generation of the hydropower station can be more accurate, thereby solving the problem of poor compensation and regulation performance of hydropower stations for wind and solar power stations. Attached Figure Description

[0010] To more clearly illustrate the embodiments of this application, the accompanying drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0011] Figure 1 A flowchart illustrating the operational simulation method provided for some embodiments of this application; Figure 2 A flowchart illustrating the construction of linear equations provided for some embodiments of this application; Figure 3 A schematic diagram of a three-dimensional coordinate system provided for some embodiments of this application; Figure 4 Schematic diagram of a running simulation device provided for some embodiments of this application; Figure 5 A schematic diagram of the modules of an electronic device provided for some embodiments of this application. Detailed Implementation

[0012] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the protection scope of this application.

[0013] It should be noted that, in the description of this application, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. The terms "first," "second," etc., in this application are used to distinguish similar objects and are not used to describe a specific order or sequence.

[0014] To enable those skilled in the art to better understand the present application, the present application will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0015] Integrated energy bases typically include various types of power plants, such as hydroelectric power stations, wind and solar power stations (i.e., wind power stations and photovoltaic power stations), and thermal power stations. During the design or operation phases of an integrated energy base, simulating the operation of each type of power plant can determine the rationality of the design or the appropriate dispatch strategy. This allows for the maximization of natural resource utilization and minimization of environmental pollution while meeting power supply demands. For example, since hydroelectric resources can be stored, when wind and solar resources are abundant, the power generation flow of hydroelectric power stations can be reduced while that of wind and solar power stations is increased. This maximizes the utilization of wind and solar resources while storing hydroelectric resources. When wind and solar resources are insufficient, the power generation flow of hydroelectric power stations can be increased based on stored hydroelectric resources, while the power generation flow of wind and solar power stations can be reduced. This ensures that normal power supply demands are met.

[0016] Currently, the hydropower output function of a hydroelectric power station is used in the simulation of the operation of an integrated energy base. Since the hydropower output function is a quadratic function of the water head and the power generation flow rate, solving the simulation model of the integrated energy base using current solvers presents practical problems such as low solution efficiency or even the inability to find a feasible solution, resulting in high simulation difficulty.

[0017] Therefore, this application provides an operation simulation method for integrated energy bases, which can solve the problem of high difficulty in simulating the operation of integrated energy bases. The operation simulation method can be applied to electronic devices. Electronic devices can include, but are not limited to, tablet computers, laptops, desktop computers, servers, etc. (See also...) Figure 1 This is a flowchart illustrating the operation simulation method provided in some embodiments of this application. Figure 1 The simulation method includes steps S101 to S104.

[0018] Step S101: Obtain hydropower generation information and wind and solar power generation information. The hydropower generation information includes the operation constraint information of the hydropower station, the maximum reservoir capacity, the minimum reservoir capacity, the maximum power generation flow, and the minimum power generation flow. The wind and solar power generation information includes the wind and solar power generation sequence obtained with the first duration as the statistical dimension and the maximum power transmission capacity supported by the transmission channel.

[0019] Specifically, hydropower information refers to the design parameters of a hydropower station. The operational constraints of a hydropower station may include, but are not limited to, the following: 1) Characteristic parameters such as installed capacity, characteristic water levels, minimum discharge flow, and overall power output coefficient; 2) Characteristic curves such as water level-reservoir capacity relationship curve, discharge flow-tailwater level curve, head-expected output curve, and water level-maximum discharge flow curve; 3) The water balance of the hydropower station and the water level, flow rate, and output at different times. Among them, the output is the power generation capacity of the hydropower station.

[0020] Wind and solar power generation information may include, but is not limited to, wind power generation sequences and photovoltaic power generation sequences obtained by statistically analyzing typical years using the first time period as the statistical dimension. A typical year refers to a specific historical year that can represent long-term climate or resource fluctuation characteristics. The first time period is in hours; that is, the wind power generation sequence includes 8760 wind power generation data points obtained by statistical analysis in hourly units, and the photovoltaic power generation sequence includes 8760 photovoltaic power generation data points obtained by statistical analysis in hourly units.

[0021] Step S102: Divide the reservoir into multiple discrete reservoir capacities within the range of maximum and minimum reservoir capacity, and divide the hydropower station into multiple discrete power generation flows within the range of maximum and minimum power generation flows. Based on the discrete reservoir capacities and discrete power generation flows, determine the multiple discrete power generation capacities of the hydropower station.

[0022] Specifically, within the interval between the maximum and minimum storage capacity, discretizing the interval according to the first step length yields a storage capacity vector comprising multiple discrete storage capacities. Similarly, within the interval between the maximum and minimum power generation flow rates, discretizing this interval according to the second step size yields a power generation flow vector comprising multiple discrete power generation flow rates. .

[0023] Based on storage capacity vector and power generation flow vector The discrete power output P of a hydroelectric power station under different combinations of reservoir capacity and power generation flow can be calculated by following these steps: 1) Based on the reservoir capacity of the hydropower station at the beginning of a certain period Given the power generation flow rate q, calculate the reservoir capacity of the hydropower station at the end of the specified period. ; 2) Based on storage capacity The initial water level in front of the dam of the hydropower station was calculated using the water level-reservoir capacity curve. and based on storage capacity Based on the water level-reservoir capacity curve of the hydropower station, the final water level in front of the dam is calculated. ; 3) Based on the power generation flow rate q and the discharge flow-tailwater level curve of the hydropower station, the tailwater level of the hydropower station is calculated. ; 4) Based on expression The head h of the hydroelectric power station is calculated. 5) Based on the expression P=kqh, the discrete power generation P of the hydropower station can be calculated, where k is the comprehensive output efficiency coefficient of the hydropower station and q is the power generation flow of the hydropower station.

[0024] Step S103: Based on discrete reservoir capacity, discrete power generation flow rate, and discrete power generation output, the hydropower output function of the hydropower station is divided into multiple linear equations.

[0025] See also Figure 2 The above is a flowchart illustrating the construction of linear equations for some embodiments of this application. Figure 2 This includes steps S201 to S203.

[0026] Step S201: The discrete reservoir capacity, discrete power generation flow, and discrete power generation are used as discrete points on the three coordinate axes in a three-dimensional coordinate system, and multiple intersection points of the discrete reservoir capacity, discrete power generation flow, and discrete power generation are determined in the three-dimensional coordinate system.

[0027] Step S202: Based on multiple intersection points, construct multiple convex hull planes and determine the plane equations of each convex hull plane.

[0028] Step S203: The plane equations are used as multiple linear equations for the hydropower output function.

[0029] For ease of understanding, please refer to the following: Figure 3 This is a schematic diagram of a three-dimensional coordinate system provided for some embodiments of this application. Figure 3In this model, discrete reservoir capacity, discrete power generation flow, and discrete power generation are used as the three axes of a three-dimensional coordinate system. After labeling each discrete reservoir capacity, discrete power generation flow, and discrete power generation on its corresponding coordinate axis, connecting these discrete points yields the surface shown by the dashed line. The surface equation of this surface is the hydropower output function of the hydropower station, and the hydropower output function is a quadratic function.

[0030] To reduce the difficulty of solving the simulation of the integrated energy base operation, multiple convex hull planes can be used to approximate the surface. Thus, the surface equation can be approximately represented by the plane equations of these multiple convex hull planes. Since the plane equations are linear functions, the difficulty of solving the problem can be greatly reduced.

[0031] Specifically, the above involves constructing multiple convex hull planes and determining the plane equations of each convex hull plane, including: Construct the plane equations for each convex hull plane, as shown in expression (1).

[0032] (1) in, V represents discrete power generation flow, and V represents discrete reservoir capacity. Let represent the convex hull plane formed by the intersection points of the m-th and m+1-th discrete power generation flows and the n-th and n+1-th discrete reservoir capacities. Describes the convex hull plane The plane equations are given, m=1,…,M-1, n=1,…,N-1, where M is the number of discrete power generation flows and N is the number of discrete reservoir capacities. , , Let be the plane parameters to be solved.

[0033] Using the following expression (2) as the objective function and the following expressions (3) and (4) as constraints, the planar parameters of each convex hull plane are obtained by solving.

[0034] (2) (3)

[0035] in, This represents the fitted value of the theoretical power generation calculated based on the m-th discrete power generation flow and the n-th discrete reservoir capacity. Represents the plane with respect to the convex hull. Adjacent convex hull planes, Represents the convex hull plane The plane equation fits the discrete power generation value at the center point. Represents the convex hull plane The plane equation fits the discrete power generation value at the center point. This represents the intersection point formed by the m-th discrete power generation flow and the n-th discrete reservoir capacity. Indicates the intersection point Fitted values ​​of discrete power generation at the location, , , , Indicates the intersection point Four connected convex hull planes, , , , The equations of the planes representing the four convex hull planes at their intersection points The fitted value of discrete power generation at the location.

[0036] Specifically, the meaning of the above expression (2) is: in the process of solving the plane parameters, the difference between the power generation obtained based on the plane equation and the theoretical power generation needs to be minimized, so as to ensure the approximation degree between the convex hull plane represented by the plane equation and the original surface.

[0037] The meaning of the above expression (3) is: in the plane At the center point, the plane The value of needs to be smaller than the values ​​of other planes. This ensures the overall consistency of the convex hull plane.

[0038] The meaning of the above expression (4) is that at the intersection of multiple convex hull planes, the discrete power generation needs to be the minimum value among the multiple convex hull planes. In this way, the problem of overestimation of power output caused by the unsmooth transition between planes is prevented.

[0039] Step S104: Using the second duration as the statistical dimension, perform statistics on the wind and solar power generation sequence to obtain the power output fluctuation coefficient of the wind and solar power station. The second duration is greater than or equal to the first duration.

[0040] Specifically, the power output fluctuation coefficient represents the magnitude of fluctuations in wind and solar power generation. For any second time period, the standard deviation and mean of wind and solar power generation within that second time period can be determined, and based on the standard deviation and mean, the power output fluctuation coefficient of the wind and solar power station within that second time period can be determined.

[0041] In this embodiment, the second duration is 10 days. Based on the wind power generation sequence and photovoltaic power generation sequence obtained in step S101, the wind power generation and photovoltaic power generation at corresponding time points can be added together to obtain the wind and solar power generation at the corresponding time points. For example, the wind power generation and photovoltaic power generation in the first hour are added together to obtain the wind and solar power generation in the first hour; the wind power generation and photovoltaic power generation in the second hour are added together to obtain the wind and solar power generation in the second hour.

[0042] Based on the hourly wind and solar power generation, the standard deviation and mean of wind and solar power generation within the t-th ten-day period can be calculated. Based on the standard deviation and mean of wind and solar power generation within the t-th ten-day period, the output fluctuation coefficient of the wind and solar power station within the t-th ten-day period can be determined. Here, t is an integer between 1 and 36 (inclusive). Specifically, the output fluctuation coefficient of the wind and solar power station within the t-th ten-day period can be calculated based on expression (5).

[0043] (5) in, This represents the output fluctuation coefficient of a wind and solar power station within ten days t. This represents the standard deviation of wind and solar power generation within the t-th period. This represents the average wind and solar power generation within the t-th ten-day period. This represents the nth wind and solar power generation within the t-th ten-day period. This represents the amount of wind and solar power generated within the t-th ten-day period, where n is an integer between 1 and T3 (inclusive).

[0044] After normalizing the output fluctuation coefficients for each ten-day period, the output fluctuation coefficient sequence for a typical year can be obtained. .

[0045] Understandably, in practical applications, the second duration can be determined according to actual needs. For example, the second duration could also be one month. This application does not impose any restrictions on the specific value of the second duration.

[0046] Step S105: Using the second duration as the statistical dimension and the maximum transmission capacity supported by the transmission channel as the constraint, the wind and solar power generation sequence is statistically analyzed to obtain the curtailment coefficient of the wind and solar power station.

[0047] Specifically, the curtailment factor includes a characterization of the amount of electricity curtailed by wind and solar power plants. For any second time period, the sum of wind and solar power generation within that second time period can be determined, and based on the sum of wind and solar power generation and the maximum transmission capacity supported by the transmission channel, the curtailment factor of the wind and solar power plant within that second time period can be determined.

[0048] In this embodiment, the curtailment factor of the wind and solar power station in the t-th ten-day period can be calculated based on expression (6).

[0049] (6) in, This represents the curtailment coefficient of wind and solar power stations within the t-th ten-day period. This represents the amount of electricity wasted by wind and solar power plants within the t-th ten-day period. This represents the amount of wind and solar power generated by the wind and solar power station within the t-th ten-day period.

[0050] After normalizing the curtailment coefficients for each ten-day period, the curtailment coefficient sequence for wind and solar power plants in a typical year can be obtained. .

[0051] Step S106: Determine the output weight of the hydropower station based on the output fluctuation coefficient and curtailment coefficient of the wind and solar power station.

[0052] Specifically, the output weight is used to correct the power generation of a hydroelectric power station. The corrected power generation is less than the uncorrected power generation. The larger the output weight, the greater the correction to the power generation. In other words, the larger the output weight, the smaller the corrected power generation.

[0053] The output fluctuation coefficient and curtailment coefficient of wind and solar power plants can be directly proportional to the output weight of hydropower plants. The reasons for this proportionality are explained below.

[0054] 1) A large output fluctuation coefficient indicates that the power generation of wind and solar power plants is relatively high. In this case, the power generation of hydropower plants can be further reduced (i.e., the output weight needs to be greater). This avoids excessive power waste and natural resource waste when the power generation of wind and solar power plants is high, as this would also lead to high power generation of hydropower plants. Of course, it is understandable that when the output fluctuation coefficient is large, the power generation of wind and solar power plants may also be relatively low. When the power generation of wind and solar power plants is low, in order to ensure normal power supply, the power generation of other power plants in the power system (such as thermal power plants) can be controlled. In this way, on the one hand, the normal power supply of the grid can be guaranteed, and on the other hand, natural energy sources such as wind, solar and hydropower can be fully utilized, avoiding the waste of natural energy.

[0055] 2) When the curtailment coefficient is relatively large, it indicates that the power generation of wind and solar power stations is already excessive. In this case, the power generation of hydropower stations can be further reduced (i.e., the output weight needs to be greater). In this way, we can avoid the situation where the power generation of wind and solar power stations is excessive, and the power generation of hydropower stations is also high, thus causing excessive curtailment and waste of natural resources.

[0056] In this embodiment, for any second time period, the output weight of the wind and solar power station within that second time period can be determined based on the output fluctuation coefficient and curtailment coefficient of the wind and solar power station within that second time period. That is, based on the output fluctuation coefficient and curtailment coefficient of the wind and solar power station within ten-day periods t, the output weight of the wind and solar power station within ten-day periods t can be determined. Thus, the output weight sequence of the hydropower station in a typical year can be obtained. .

[0057] Based on the output weight within ten days t, the power generation of the hydroelectric power station within ten days t can be adjusted.

[0058] Furthermore, the output fluctuation coefficient and curtailment coefficient of wind and solar power plants can each have their corresponding weights. These weights characterize the influence of the output fluctuation coefficient and curtailment coefficient of wind and solar power plants on the output weight of hydropower plants. When determining the output weight of a hydropower plant based on the output fluctuation coefficient and curtailment coefficient of wind and solar power plants, a first weight for the output fluctuation coefficient and a second weight for the curtailment coefficient can be obtained. Based on the first and second weights, a weighted fusion calculation can be performed on the output fluctuation coefficient and curtailment coefficient of wind and solar power plants to obtain the output weight.

[0059] Specifically, the output weight of the hydroelectric power station in the t-th ten-day period can be determined based on expression (7). .

[0060] (7) in, Indicates the first weight. Indicates the second weight. .

[0061] Step S107: Based on the linear equations and the operating constraints of the hydropower station, simulate the operation of the hydropower station, correct the power generation of the hydropower station based on the output weight, and simulate the operation of the wind and solar power station based on the wind and solar power generation information.

[0062] Specifically, in simulating the operation of a wind and solar power station based on wind and solar power generation information, the actual process involves having the wind and solar power generation model simulate power generation according to the wind power generation sequence and the photovoltaic power generation sequence. This process is a conventional technique and will not be elaborated upon here.

[0063] Furthermore, based on linear equations, the operation of a hydroelectric power station is simulated, including... Obtain the power generation flow and reservoir capacity used to simulate the operation of a hydroelectric power station; Substituting the power generation flow and reservoir capacity into the plane equations of each convex hull plane, the multiple power generation capacities are obtained by solving the equations. Among multiple power generation capacities, the smallest power generation capacity is selected as the power generation capacity of the hydroelectric power station.

[0064] Specifically, this process can be represented by expression (8): (8) in, Let represent the power generation of the hydroelectric power station at time t, and let a, b, and c represent the planar parameters of each convex hull plane. This represents the power generation flow rate of the hydroelectric power station at time t. This represents the reservoir capacity of the hydroelectric power station at time t. Let represent the reservoir capacity of the i-th sub-hydropower station at time t+1.

[0065] The meaning of expression (8) is as follows: Substitute the power generation flow and reservoir capacity into the plane equations of each convex hull plane to obtain multiple power generation capacities; among the multiple power generation capacities, select the smallest power generation capacities as the power generation capacities of the hydroelectric power station. In this way, in the case of errors in the linear equations leading to overestimation of the calculation results, it can effectively prevent the actual power generation capacities of the hydroelectric power station from exceeding its actual capacity.

[0066] In summary, in the technical solutions of some embodiments of this application, on the one hand, by decomposing the hydropower output function of a hydropower station into multiple linear equations, a quadratic function problem can be transformed into a linear function problem. This significantly reduces the difficulty of solving the problem during the simulation of the operation of a comprehensive energy base based on the hydropower output function, thereby reducing the overall simulation difficulty. On the other hand, based on the output fluctuation coefficient and curtailment coefficient of wind and solar power stations, the output weight of the hydropower station is determined, and the power generation of the hydropower station is corrected based on this output weight. Thus, the final power generation of the hydropower station can be obtained more accurately, thereby solving the problem of poor compensation and regulation performance of hydropower stations for wind and solar power stations.

[0067] Furthermore, the multiple convex hull planes in this application can only approximate the original surface composed of discrete reservoir capacity, discrete power generation flow, and discrete power generation, but cannot completely coincide with the original surface. Therefore, the power generation calculated based on the plane equation of the convex hull plane should have an error compared with the theoretical power generation of the hydropower station. To improve the simulation accuracy, in some embodiments, when the grid-connected power of the hydropower station and the wind and solar power station reaches the maximum transmission power allowed by the transmission channel and there is power curtailment in the integrated energy base, if the power generation of the hydropower station in any target time period obtained based on the plane equation of the convex hull plane is greater than the theoretical power generation of the hydropower station in the target time period, it indicates that the power generation capacity of the hydropower station is overestimated. In this case, the power generation of the hydropower station and the power curtailment of the wind and solar power station can be reduced. Conversely, if the grid-connected power of hydropower stations and wind and solar power stations reaches the maximum transmission capacity allowed by the transmission channel and there is curtailment of power in the integrated energy base, and the power generation of the hydropower station in any target time period obtained based on the plane equation of the convex hull plane is less than the theoretical power generation of the hydropower station in the target time period, it indicates that the power generation capacity of the hydropower station is underestimated. In this case, the power generation of the hydropower station and the curtailment power of the wind and solar power stations can be increased. In this way, the simulation accuracy can be guaranteed when simulating the operation of hydropower stations and wind and solar power stations.

[0068] In some embodiments, the hydroelectric power station includes multiple sub-hydroelectric power stations, each sub-hydroelectric power station having its own corresponding power generation flow rate, reservoir capacity, and convex hull plane; the method of this application further includes: Using expression (9) as a constraint, the power generation flow and reservoir capacity of the i-th sub-hydropower station are substituted into the convex hull plane corresponding to the i-th sub-hydropower station to obtain the power generation of the i-th sub-hydropower station at multiple time points.

[0069] (9) in, Let represent the power generation of the i-th sub-hydropower station at time t, and let a, b, and c represent the planar parameters of the convex hull plane corresponding to the i-th sub-hydropower station. This represents the power generation flow rate of the i-th sub-hydropower station at time t. This represents the reservoir capacity of the i-th hydroelectric power station at time t. Let represent the reservoir capacity of the i-th sub-hydropower station at time t+1.

[0070] Specifically, when a hydropower station includes multiple sub-hydropower stations, these multiple sub-hydropower stations can be cascade hydropower stations. The following uses a cascade hydropower station as an example to illustrate the complete process of the method in this application.

[0071] 1) Based on the above steps S101 to S103, the hydropower output function of each sub-hydropower station is decomposed into multiple linear equations.

[0072] 2) Based on the linear equations and operational constraints of each sub-hydropower station, the operation of the sub-hydropower station is simulated, and based on the wind and solar power generation information, the operation of the wind and solar power station is simulated, so as to obtain the power generation and curtailment of the sub-hydropower station and the wind and solar power station at multiple time points.

[0073] 3) Based on the power generation and theoretical power generation of each sub-hydropower station obtained from the solution, the power curtailment of at least some of the sub-hydropower stations obtained from the solution is corrected.

[0074] 4) The total grid-connected power of the hydro-wind-solar power station in the target time period, calculated as follows. Less than the maximum allowable power of the transmission channel In the case where the power generation of the hydropower station obtained based on the plane equation of the convex hull plane in the target time period is less than the theoretical power generation of the hydropower station in the target time period, the total power generation of the hydropower station, wind power station, and solar power station in the target time period can be recalculated. If the calculated total grid-connected power of the hydropower station, wind power station, and solar power station is greater than the maximum power allowed by the transmission channel, the total curtailment power of the integrated energy base in the target time period is corrected to... - If the power generation of the hydropower station in the target time period obtained based on the plane equation of the convex hull plane is greater than the theoretical power generation of the hydropower station in the target time period, then the power generation of the hydro-wind-solar power station in the target time period can be corrected.

[0075] Furthermore, in some embodiments, when simulating the operation of hydroelectric power stations and wind and solar power stations, the method of this application may further include: using a third duration as a first scheduling cycle, using a second duration as a first scheduling period, using the maximum grid-connected power of the integrated energy base as the target, using at least the operational constraint information of the hydroelectric power station as a constraint, correcting the power generation of the hydroelectric power station based on the output weight to obtain the target water level of the hydroelectric power station at the beginning and end of each first scheduling period, the third duration being greater than or equal to the second duration, and using the second duration as a second scheduling cycle, using the first duration as a second scheduling period, using the target water level of the hydroelectric power station at the beginning and end of each first scheduling period as a constraint to determine the power generation of the hydroelectric power station in each second scheduling period.

[0076] Specifically, the third time period can be a year, the second time period can be a ten-day period, and the first time period can be an hour. This simulation method can also be called multi-timescale layered nested simulation. In this way, the number of parameters to be solved can be reduced in each simulation calculation, thus reducing the difficulty of the solution. For example, if the scheduling cycle is directly set to a year and the scheduling period is set to an hour, it is necessary to solve for the power generation and curtailment of wind, solar and hydropower stations over 8760 hours. The number of parameters to be solved is too large, which may lead to problems that cannot be solved. However, the multi-timescale layered nested simulation method can effectively reduce the number of variables to be solved. Simply put, if a continuous simulation of 8760 hours is performed, at least 8760 variables need to be solved (i.e., the power generation of the hydropower station over 8760 hours), which is extremely difficult to solve. However, by using the third time period as the first scheduling cycle and the second time period as the first scheduling period, and solving for the target water level of the hydropower station at the beginning and end of each ten-day period, the number of variables to be solved is significantly reduced. Similarly, when using the second duration as the second scheduling cycle and the third duration as the second scheduling period, to solve for the hourly power generation of the hydroelectric power station within each ten-day period, since the number of hours in one ten-day period is much lower than the number of hours in one year, the number of variables to be solved is also significantly reduced. Therefore, the number of variables to be solved can also be greatly reduced.

[0077] In some embodiments, when the third duration is used as the first scheduling cycle and the second duration is used as the first scheduling period, expression (10) can be used as the first objective function, and at least expressions (11) to (17) can be used as constraints to solve for the target water level of the hydropower station at the beginning and end of each first scheduling period. It should be noted that, in addition to expressions (11) to (17), other constraints may also be included, such as the sum of the grid-connected power of hydropower, wind power and solar power stations not exceeding the maximum transmission power allowed by the transmission channel. This application does not impose any restrictions on the constraints.

[0078] (10) in, This represents the grid-connected electricity volume of the integrated energy base. T1 is the number of ten-day periods within a year (default is 36), and N1 is the number of sub-hydropower stations. Let n be the power generation capacity of the nth sub-hydropower station during time period t. This represents the total power generation of n sub-hydropower stations during the t-th first scheduling period. The power output weight of the hydroelectric power station Let be the grid-connected power of the wind power station in the t-th first dispatch period. Let be the grid-connected power of the photovoltaic power station in the t-th first dispatch period; Let be the power curtailment of the wind and solar power station during the t-th first dispatch period. This represents the output weight of n sub-hydropower stations in the t-th first scheduling period. By introducing output weights, the scheduling accuracy of hydropower stations and wind and solar power stations can be improved.

[0079] (11) (12) (13) (14) (15) (16) (17) In the above expression (11), This represents the inflow rate into the reservoir of the nth sub-hydropower station during the tth first scheduling period. Indicates upstream The inflow between the downstream nth sub-hydropower station and the tth first scheduling period. express The water flow time delay to the nth sub-hydropower station , They represent the upstream The power generation flow and water abandonment flow during the t-th first scheduling period.

[0080] In the above expression (12), , Let these represent the reservoir capacity of the nth sub-hydropower station during the t-th and (t-1)th first scheduling periods, respectively. It is the discharge flow of the nth sub-hydropower station during the tth first scheduling period.

[0081] In the above expression (13), , These represent the reservoirs of the nth sub-hydropower station in the nth... Time period and the Water level in front of the dam during the specified time period , These represent the reservoirs of the nth sub-hydropower station. Upper and lower limits of water level for a given period of time This represents the water level fluctuation limit of the reservoir of the nth sub-hydropower station in adjacent time periods.

[0082] In the above expression (14), , These represent the time periods of the nth sub-hydropower station. The minimum and maximum discharge flow rates.

[0083] In the above expression (15), This indicates that the nth sub-hydropower station is in the time period. Power generation capacity, , These represent the time periods of the nth sub-hydropower station. The minimum and maximum permissible power generation capacity.

[0084] In the above expression (16), This represents the total power generation of the cascade hydropower stations during time period t.

[0085] In the above expression (17), , Let represent the water level-storage capacity function and the tailwater level-discharge flow function corresponding to the reservoir of the nth sub-hydropower station, respectively.

[0086] In some embodiments, when the second duration is used as the second scheduling period and the first duration is used as the second scheduling time period, expression (18) can be used as the second objective function and expressions (19) to (20) can be used as constraints.

[0087] (18) in, This indicates the amount of electricity generated by the integrated energy base within a ten-day period. This indicates the number of hours within a ten-day period. The total grid-connected power of the hydro-wind-solar power station during time period t; The total power generation of the cascade hydropower stations during time period t; This represents the grid-connected power output of the wind power station during time period t. The grid-connected power of the photovoltaic power station during time period t; This represents the amount of power curtailed by the wind and solar power station during time period t.

[0088] (19) (20) In the above expression (20), The minimum allowable power for the transmission channels of the integrated energy base. The maximum allowable transmission power for the power transmission channels of the integrated energy base.

[0089] Through the above description of the embodiments, those skilled in the art can clearly understand that the methods according to the above embodiments can be implemented by means of software plus necessary general-purpose hardware platforms. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation method.

[0090] Corresponding to the method, this application also provides an operation simulation device for an integrated energy base. (See also...) Figure 4 This is a schematic diagram of a running simulation device provided in some embodiments of this application. Figure 4 The simulation device includes: The information acquisition module 401 is used to acquire hydropower generation information and wind and solar power generation information. The hydropower generation information includes the operation constraint information of the hydropower station, the maximum reservoir capacity, the minimum reservoir capacity, the maximum power generation flow, and the minimum power generation flow. The wind and solar power generation information includes the wind and solar power generation sequence obtained with the first duration as the statistical dimension and the maximum power transmission capacity supported by the transmission channel. The interval division module 402 is used to divide the interval between the maximum and minimum reservoir capacity into multiple discrete reservoir capacities, and to divide the interval between the maximum and minimum power generation flow into multiple discrete power generation flow, and to determine multiple discrete power generation capacities of the hydropower station based on the discrete reservoir capacities and discrete power generation flow. Equation building module 403 is used to divide the hydropower output function of a hydropower station into multiple linear equations based on discrete reservoir capacity, discrete power generation flow, and discrete power generation. The fluctuation coefficient determination module 404 is used to statistically analyze the wind and solar power generation sequence using the second duration as the statistical dimension to obtain the power output fluctuation coefficient of the wind and solar power station. The second duration is greater than or equal to the first duration. The curtailment coefficient determination module 405 is used to statistically analyze the wind and solar power generation sequence using the second duration as the statistical dimension and the maximum transmission capacity supported by the transmission channel as the constraint condition, so as to obtain the curtailment coefficient of the wind and solar power station. The weight determination module 406 is used to determine the output weight of hydropower stations based on the output fluctuation coefficient and curtailment coefficient of wind and solar power stations. The simulation module 407 is used to simulate the operation of a hydropower station based on linear equations and the operating constraint information of the hydropower station, and to correct the power generation of the hydropower station based on the output weight, and to simulate the operation of a wind and solar power station based on wind and solar power generation information.

[0091] See also Figure 5 The embodiments of this application also provide an electronic device, including a memory 10 and a processor 20, wherein the memory 10 stores a computer program and the processor 20 is configured to run the computer program to perform the steps in any of the above embodiments of the integrated energy base operation simulation method.

[0092] Embodiments of this application also provide a computer-readable storage medium storing a computer program, wherein the computer program is configured to execute the steps in any of the above embodiments of the integrated energy base operation simulation method.

[0093] In one exemplary embodiment, the aforementioned computer-readable storage medium may include, but is not limited to, various media capable of storing computer programs, such as a USB flash drive, read-only memory (ROM), random access memory (RAM), portable hard disk, magnetic disk, or optical disk.

[0094] The embodiments of this application also provide a computer program product, which includes a computer program that, when executed by a processor, implements the steps in any of the above embodiments of the integrated energy base operation simulation method.

[0095] Embodiments of this application also provide another computer program product, including a non-volatile computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps in any of the above embodiments of the integrated energy base operation simulation method.

[0096] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0097] The above provides a detailed description of the operation simulation method, apparatus, and equipment for an integrated energy base provided in this application. Specific examples have been used to illustrate the principles and implementation methods of this application. The descriptions of the embodiments above are merely for the purpose of helping to understand the method and its core ideas. It should be noted that those skilled in the art can make various improvements and modifications to this application without departing from its principles, and these improvements and modifications also fall within the protection scope of the claims of this application.

Claims

1. A method for simulating the operation of an integrated energy base, characterized in that, The integrated energy base includes hydropower stations and wind and solar power stations, and the method includes: Acquire hydropower and wind and solar power information. The hydropower information includes the operation constraints of the hydropower station, maximum reservoir capacity, minimum reservoir capacity, maximum power generation flow, and minimum power generation flow. The wind and solar power information includes the wind and solar power generation sequence obtained with a first duration as the statistical dimension and the maximum power transmission capacity supported by the transmission channel. In the interval between the maximum reservoir capacity and the minimum reservoir capacity, multiple discrete reservoir capacities are divided, and in the interval between the maximum power generation flow and the minimum power generation flow, multiple discrete power generation flows are divided, and based on the discrete reservoir capacities and the discrete power generation flows, multiple discrete power generation capacities of the hydropower station are determined. Based on the discrete reservoir capacity, the discrete power generation flow rate, and the discrete power generation output, the hydropower output function of the hydropower station is divided into multiple linear equations; Using the second duration as a statistical dimension, the wind and solar power generation sequence is statistically analyzed to obtain the power output fluctuation coefficient of the wind and solar power station. The second duration is greater than or equal to the first duration. Using the second duration as the statistical dimension and the maximum transmission capacity supported by the transmission channel as the constraint, the wind and solar power generation sequence is statistically analyzed to obtain the curtailment coefficient of the wind and solar power station. Based on the power output fluctuation coefficient and the curtailment coefficient of the wind and solar power station, the power output weight of the hydropower station is determined. Based on the linear equation and the operating constraint information of the hydropower station, the operation of the hydropower station is simulated, and the power generation of the hydropower station is corrected based on the output weight. Based on the wind and solar power generation information, the operation of the wind and solar power station is simulated.

2. The method according to claim 1, characterized in that, Based on the discrete reservoir capacity, the discrete power generation flow rate, and the discrete power generation output, the output function of the hydropower station is divided into multiple linear equations, including: The discrete reservoir capacity, the discrete power generation flow, and the discrete power generation are used as discrete points on three coordinate axes in a three-dimensional coordinate system, and multiple intersection points of the discrete reservoir capacity, the discrete power generation flow, and the discrete power generation are determined in the three-dimensional coordinate system. Based on the multiple intersection points, multiple convex hull planes are constructed, and the plane equations of each of the convex hull planes are determined; The plane equation is used as multiple linear equations for the hydropower output function.

3. The method according to claim 2, characterized in that, The construction of multiple convex hull planes and the determination of the plane equations of each of the convex hull planes include: Construct the plane equations of each of the convex hull planes ,in, V represents discrete power generation flow, and V represents discrete reservoir capacity. Let represent the convex hull plane formed by the intersection points of the m-th and m+1-th discrete power generation flows and the n-th and n+1-th discrete reservoir capacities. Describes the convex hull plane The plane equations are given, m=1,…,M-1, n=1,…,N-1, where M is the number of discrete power generation flows and N is the number of discrete reservoir capacities. , , Let be the plane parameters to be solved; In expression Using the following expression as the objective function and constraints, the planar parameters of each convex hull plane are obtained by solving: in, This represents the fitted value of the theoretical power generation calculated based on the m-th discrete power generation flow and the n-th discrete reservoir capacity. Represents the plane with respect to the convex hull. Adjacent convex hull planes, Describes the convex hull plane The plane equation fits the discrete power generation value at the center point. Describes the convex hull plane The plane equation fits the discrete power generation value at the center point. This represents the intersection point formed by the m-th discrete power generation flow and the n-th discrete reservoir capacity. Indicates the intersection point The fitted value of the discrete power generation at the location, , , , Indicates the intersection point Four connected convex hull planes, , , , The plane equations representing the four convex hull planes at their intersection points The fitted value of discrete power generation at the location.

4. The method according to claim 3, characterized in that, Based on the aforementioned linear equations, the operation of the hydroelectric power station is simulated, including: Obtain the power generation flow and reservoir capacity used to simulate the operation of a hydroelectric power station; Substituting the power generation flow rate and the reservoir capacity into the plane equations of each of the convex hull planes, the multiple power generation capacities are obtained by solving the equations. Among the multiple power generation capacities, the smallest power generation capacity is selected as the power generation capacity of the hydroelectric power station.

5. The method according to claim 4, characterized in that, The integrated energy base also includes power transmission channels; the method further includes: If the grid-connected power of the hydropower station and the wind and solar power station reaches the maximum transmission power allowed by the transmission channel and there is power curtailment in the integrated energy base, and if the power generation of the hydropower station in any target time period, obtained based on the plane equation of the convex hull plane, is greater than the theoretical power generation of the hydropower station in the target time period, then the power generation of the hydropower station and the power curtailment of the wind and solar power station shall be reduced.

6. The method according to claim 4, characterized in that, The method further includes: If the grid-connected power of the hydropower station and the wind and solar power station reaches the maximum transmission power allowed by the transmission channel and there is curtailment of power in the integrated energy base, and the power generation of the hydropower station in any target time period obtained based on the plane equation of the convex hull plane is less than the theoretical power generation of the hydropower station in the target time period, then the power generation of the hydropower station and the curtailment of power of the wind and solar power station shall be increased.

7. The method according to claim 1, characterized in that, The method further includes, during the simulation of the operation of the hydroelectric power station and the wind and solar power station: The third duration is used as the first scheduling cycle, the second duration is used as the first scheduling period, the maximum on-grid power generation of the integrated energy base is used as the target, and the operation constraint information of the hydropower station is used as the constraint. The power generation of the hydropower station is corrected based on the output weight to obtain the target water level of the hydropower station at the beginning and end of each first scheduling period. The third duration is greater than or equal to the second duration. Using the second duration as the second scheduling cycle, the first duration as the second scheduling period, and the target water level at the beginning and end of each first scheduling period as a constraint, the power generation of the hydropower station in each second scheduling period is determined.

8. An operation simulation device for an integrated energy base, characterized in that, The integrated energy base includes a hydroelectric power station and a wind and solar power station, and the equipment includes: The information acquisition module is used to acquire hydropower generation information and wind and solar power generation information. The hydropower generation information includes the operation constraint information of the hydropower station, the maximum reservoir capacity, the minimum reservoir capacity, the maximum power generation flow, and the minimum power generation flow. The wind and solar power generation information includes the wind and solar power generation sequence obtained with the first duration as the statistical dimension and the maximum power transmission capacity supported by the transmission channel. The interval division module is used to divide the interval between the maximum reservoir capacity and the minimum reservoir capacity into multiple discrete reservoir capacities, and to divide the interval between the maximum power generation flow and the minimum power generation flow into multiple discrete power generation flow rates, and to determine multiple discrete power generation capacities of the hydropower station based on the discrete reservoir capacities and the discrete power generation flow rates. An equation construction module is used to divide the hydropower output function of the hydropower station into multiple linear equations based on the discrete reservoir capacity, the discrete power generation flow rate, and the discrete power generation output. The fluctuation coefficient determination module is used to statistically analyze the wind and solar power generation sequence using the second duration as a statistical dimension to obtain the power output fluctuation coefficient of the wind and solar power station, wherein the second duration is greater than or equal to the first duration. The curtailment coefficient determination module is used to statistically analyze the wind and solar power generation sequence using the second duration as a statistical dimension and the maximum power transmission capacity supported by the transmission channel as a constraint condition, so as to obtain the curtailment coefficient of the wind and solar power station. The weight determination module is used to determine the output weight of the hydropower station based on the output fluctuation coefficient and the curtailment coefficient of the wind and solar power station. The simulation module is used to simulate the operation of the hydropower station based on the linear equation and the operating constraint information of the hydropower station, and to correct the power generation of the hydropower station based on the output weight, and to simulate the operation of the wind and solar power station based on the wind and solar power generation information.

9. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the operation simulation method of the integrated energy base as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, wherein the computer program, when executed by a processor, implements the operation simulation method of the integrated energy base as described in any one of claims 1 to 7.

Images