Operational simulation methods, devices, and equipment for integrated energy bases

Abstract

This application discloses a method, apparatus, and equipment for simulating the operation of integrated energy bases, relating to the fields of integrated energy base planning, design, and operation scheduling technology. The method of this application, on the one hand, decomposes the hydropower output function of a hydropower station into multiple linear equations, thus transforming a quadratic function problem into a linear function problem. This significantly reduces the difficulty of solving the problem during the simulation of the integrated energy base's operation 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.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of integrated energy base planning and design and operation scheduling, and in particular to an operation simulation method, device and equipment of an integrated energy base. BACKGROUND

[0002] An integrated energy base usually includes power stations such as water, wind and light. During the design or operation stage of the integrated energy base, the operation process of the integrated energy base (i.e., the operation process of the water, wind and light power stations) is usually simulated. In this way, it is determined whether the design of the integrated energy base is reasonable or the scheduling strategy of the integrated energy base is determined.

[0003] At present, when the operation simulation of the integrated energy base is performed, the water power output function of the hydroelectric power station is used. Since the water power output function is a quadratic function with respect to the water head and the power generation flow, when the current some solvers are used to solve the simulation model of the integrated energy base, there is a real problem that the solving efficiency is low or even no feasible solution can be found, and the simulation difficulty is high. SUMMARY

[0004] The present application provides an operation simulation method, device, electronic equipment, computer readable storage medium and computer program product of an integrated energy base to at least solve the problem of high simulation difficulty in the related art.

[0005] The present application provides an operation simulation method of an integrated energy base, the integrated energy base including a hydroelectric power station and a wind and light power station, and the method comprising:

[0006] Obtaining hydroelectric power information and wind and light power information, the hydroelectric power information including operation constraint information, maximum reservoir capacity, minimum reservoir capacity, maximum power generation flow and minimum power generation flow of the hydroelectric power station, and the wind and light power information including a wind and light power generation sequence obtained by taking a first time length as a statistical dimension and a maximum power transmission capacity supported by a power transmission channel;

[0007] In the interval of the maximum reservoir capacity and the minimum reservoir capacity, a plurality of discrete reservoir capacities are divided, and in the interval of the maximum power generation flow and the minimum power generation flow, a plurality of discrete power generation flows are divided, and based on the discrete reservoir capacities and the discrete power generation flows, a plurality of discrete power generation powers of the hydroelectric power station are determined;

[0008] Based on the discrete reservoir capacities, the discrete power generation flows and the discrete power generation powers, the water power output function of the hydroelectric power station is divided into a plurality of linear equations;

[0009] Taking a second time length as a statistical dimension, the wind and light power generation sequence is counted to obtain an output fluctuation coefficient of the wind and light power station, and the second time length is greater than or equal to the first time length.

[0010] statistically, to obtain an abandoned electricity coefficient of the wind-solar power station, the second time length being greater than or equal to the first time length;

[0011] determine an output weight of the hydraulic power station based on the output fluctuation coefficient and the abandoned electricity coefficient of the wind-solar power station;

[0012] simulate the operation of the hydraulic power station based on the linear equations and the operation constraint information of the hydraulic power station, correct the power generation of the hydraulic power station based on the output weight, and simulate the operation of the wind-solar power station based on the wind-solar power generation information.

[0013] The application further provides a device for simulating operation of a comprehensive energy base, the comprehensive energy base comprising a hydraulic power station and a wind-solar power station, the device comprising:

[0014] an information acquisition module, configured to acquire hydraulic power information and wind-solar power information, the hydraulic power information comprising operation constraint information, maximum reservoir capacity, minimum reservoir capacity, maximum power generation flow rate and minimum power generation flow rate of the hydraulic power station, the wind-solar power information comprising a wind-solar power generation sequence obtained by taking a first time length as a statistical dimension and maximum power transmission capacity supported by a power transmission channel;

[0015] an interval division module, configured to divide a plurality of discrete reservoir capacities in an interval of the maximum reservoir capacity and the minimum reservoir capacity, divide a plurality of discrete power generation flow rates in an interval of the maximum power generation flow rate and the minimum power generation flow rate, and determine a plurality of discrete power generation powers of the hydraulic power station based on the discrete reservoir capacities and the discrete power generation flow rates;

[0016] an equation construction module, configured to divide a hydroelectric output function of the hydraulic power station into a plurality of linear equations based on the discrete reservoir capacities, the discrete power generation flow rates and the discrete power generation powers;

[0017] a fluctuation coefficient determination module, configured to statistically obtain an output fluctuation coefficient of the wind-solar power station by taking a second time length as a statistical dimension, the second time length being greater than or equal to the first time length;

[0018] an abandoned electricity coefficient determination module, configured to statistically obtain an abandoned electricity coefficient of the wind-solar power station by taking the second time length as a statistical dimension and taking the maximum power transmission capacity supported by the power transmission channel as a constraint condition;

[0019] a weight determination module, configured to determine an output weight of the hydraulic power station based on the output fluctuation coefficient and the abandoned electricity coefficient of the wind-solar power station.

[0020] simulate operation of the hydropower station based on the linear equations and the operation constraint information of the hydropower station, and correct power generation of the hydropower station based on the output weight, and simulate operation of the wind-solar power station based on the wind-solar power generation information.

[0021] The application further provides an electronic device, comprising a memory for storing a computer program, and a processor for executing the computer program to implement the steps of the operation simulation method of the integrated energy base.

[0022] The application further provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the steps of the operation simulation method of the integrated energy base.

[0023] In the technical scheme of some embodiments of the application, on the one hand, the water power output function of the hydropower station is split into multiple linear equations, so that the quadratic function problem can be converted into a linear function problem, and then the difficulty of solving can be greatly reduced in simulating the operation of the integrated energy base based on the water power output function, and thus the simulation difficulty of the integrated energy base can be reduced. On the other hand, the output weight of the hydropower station is determined based on the output fluctuation coefficient and the power abandonment coefficient of the wind-solar power station, and the power generation of the hydropower station is corrected based on the output weight. In this way, the power generation of the hydropower station obtained finally can be more accurate, and thus the problem of poor compensation adjustment performance of the hydropower station to the wind-solar power station can be solved. BRIEF DESCRIPTION OF DRAWINGS

[0024] In order to more clearly illustrate the embodiments of the application, the drawings needed in the embodiments will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the application, and other drawings can be obtained by those skilled in the art without creative labor.

[0025] Figure 1 a flowchart of the operation simulation method provided by some embodiments of the application;

[0026] Figure 2 a flowchart of the construction of the linear equation provided by some embodiments of the application;

[0027] Figure 3 a schematic diagram of the three-dimensional coordinate system provided by some embodiments of the application;

[0028] Figure 4 a module schematic diagram of the operation simulation device provided by some embodiments of the application;

[0029] Figure 5 A module schematic diagram of an electronic device provided for some embodiments of the present application. DETAILED DESCRIPTION

[0030] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, but not all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work fall within the protection scope of the present application.

[0031] It should be noted that, in the description of the present application, the terms “comprise”, “contain” or any other variants thereof are intended to cover non-exclusive inclusion, so that the process, method, article or device comprising a series of elements not only includes those elements, but also includes other elements not explicitly listed or inherent to such process, method, article or device. The terms “first”, “second” and the like in the present application are used to distinguish similar objects, and are not used to describe a specific order or sequence.

[0032] In order for those skilled in the art to better understand the technical solutions of the present application, the present application will be further described in detail below with reference to the drawings and specific embodiments.

[0033] The integrated energy base usually includes various types of power stations such as hydropower stations, wind-solar power stations (i.e. wind power stations and photovoltaic power stations), thermal power stations, etc. In the design phase or operation phase of the integrated energy base, by simulating the operation process of each type of power station, it can be determined whether the design of the integrated energy base is reasonable or the dispatching strategy of the integrated energy base is determined. In this way, the natural resources can be utilized to the maximum extent and the environmental pollution can be reduced while meeting the power supply demand. For example, since the water resources in the hydropower station can be stored, when the wind-solar resources are sufficient, the power generation flow of the hydropower station can be reduced and the power generation flow of the wind-solar power station can be increased. In this way, the wind-solar resources are utilized to the maximum extent and the water resources are stored. When the wind-solar resources are insufficient, the power generation flow of the hydropower station can be increased and the power generation flow of the wind-solar power station can be reduced based on the stored water resources. In this way, the normal power supply demand can be met.

[0034] At present, in the operation simulation process of the integrated energy base, the water power output function of the hydropower station is used. Since the water power output function is a quadratic function about the water head and the power generation flow, when the current some solvers are used to solve the simulation model of the integrated energy base, there is a real problem that the solving efficiency is low or even no feasible solution can be found, and the simulation difficulty is high.

[0035] In view of this, the application provides a method for simulating operation of an integrated energy base, which can solve the problem of high difficulty in simulating operation of an integrated energy base. The method for simulating operation can be applied to an electronic device. The electronic device can include, but is not limited to, a tablet computer, a notebook computer, a desktop computer, a server, and the like. For a better understanding of the application, reference can be made to the accompanying drawings Figure 1 A flowchart of a method for simulating operation provided for some embodiments of the application is shown. Figure 1 In some embodiments, the method for simulating operation includes steps S101-S104.

[0036] In step S101, hydropower information and wind-solar power information are obtained. The hydropower information includes operation constraint information of a hydropower station, maximum reservoir capacity, minimum reservoir capacity, maximum power generation flow rate, and minimum power generation flow rate. The wind-solar power information includes a wind-solar power generation sequence obtained by taking a first time length as a statistical dimension, and a maximum power transmission capacity supported by a power transmission channel.

[0037] Specifically, the hydropower information refers to design parameters of the hydropower station. The operation constraint information of the hydropower station can include, but is not limited to, the following information:

[0038] 1) characteristic parameters such as installed capacity, each characteristic water level, minimum discharge flow rate, and comprehensive output coefficient;

[0039] 2) characteristic curves such as water level-reservoir capacity relationship curve, discharge flow rate-tail water level curve, water head-expected output curve, and water level-maximum discharge flow rate curve;

[0040] 3) water balance of the hydropower station, and water level, flow rate, and output of each period. The output refers to the power generation power of the hydropower station.

[0041] The wind-solar power information can include, but is not limited to, a wind power generation sequence and a photovoltaic power generation sequence of a typical year obtained by taking a first time length as a statistical dimension. The typical year refers to a specific historical year that can represent long-term climate or resource fluctuation characteristics. The first time length is an hour, that is, the wind power generation sequence includes 8760 wind power generations obtained by statistical calculation in units of hours, and the photovoltaic power generation sequence includes 8760 photovoltaic power generations obtained by statistical calculation in units of hours.

[0042] In step S102, a plurality of discrete reservoir capacities are divided in the interval of the maximum reservoir capacity and the minimum reservoir capacity, and a plurality of discrete power generation flow rates are divided in the interval of the maximum power generation flow rate and the minimum power generation flow rate, and based on the discrete reservoir capacities and the discrete power generation flow rates, a plurality of discrete power generation powers of the hydropower station are determined.

[0043] Specifically, in the interval of the maximum reservoir capacity and the minimum reservoir capacity, the interval is discretized according to a first step length, and a reservoir capacity vector including a plurality of discrete reservoir capacities can be obtained. Similarly, in the interval of the maximum power generation flow and the minimum power generation flow, the interval is discretized according to the second step size, and a power generation flow vector including a plurality of discrete power generation flows can be obtained .

[0044] Based on the reservoir capacity vector and the power generation flow vector , the plurality of discrete power generation powers P of the hydropower station under different combinations of reservoir capacity and power generation flow can be calculated according to the following steps:

[0045] 1) Based on the initial reservoir capacity of the hydropower station at the beginning of a time period and the power generation flow q, the reservoir capacity of the hydropower station at the end of the time period is calculated;

[0046] 2) Based on the reservoir capacity and the water level-reservoir capacity curve of the hydropower station, the initial water level before the dam of the hydropower station is calculated, and based on the reservoir capacity and the water level-reservoir capacity curve of the hydropower station, the final water level before the dam of the hydropower station is calculated;

[0047] 3) Based on the power generation flow q and the discharge-tailwater level curve of the hydropower station, the tailwater level of the hydropower station is calculated;

[0048] 4) Based on the expression , the power generation head h of the hydropower station is calculated;

[0049] 5) Based on the expression P=kqh, the discrete power generation power 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 reference flow of the hydropower station.

[0050] Step S103, based on the discrete reservoir capacity, the discrete power generation flow and the discrete power generation power, the hydropower station power generation function is divided into a plurality of linear equations.

[0051] For reference Figure 2 , a linear equation construction flowchart is provided for some embodiments of the present application. Figure 2 Including steps S201-S203.

[0052] Step S201, the discrete reservoir capacity, the discrete power generation flow and the discrete power generation power are taken as discrete points on three coordinate axes in a three-dimensional coordinate system, and a plurality of intersection points of the discrete reservoir capacity, the discrete power generation flow and the discrete power generation power are determined in the three-dimensional coordinate system.

[0053] Step S202, based on the plurality of intersection points, constructing a plurality of convex hull planes, and determining a plane equation of each convex hull plane.

[0054] Step S203, taking the plane equation as a plurality of linear equations of the water and electricity output function.

[0055] For the convenience of understanding, combined with reference to Figure 3 The schematic diagram of the three-dimensional coordinate system provided for some embodiments of the present application. Figure 3 In the figure, the discrete reservoir capacity, the discrete power generation flow and the discrete power generation power are taken as three coordinate axes of the three-dimensional coordinate system, after marking each discrete reservoir capacity, discrete power generation flow and discrete power generation power on the corresponding coordinate axis, connecting these discrete points, the curved surface shown by the dashed line can be obtained. The curved surface equation of the curved surface is the water and electricity output function of the hydropower station, and the water and electricity output function is a quadratic function.

[0056] In order to reduce the difficulty of solving in the operation simulation process of the comprehensive energy base, a plurality of convex hull planes can be used to approximate the curved surface. In this way, the plane equation of the plurality of convex hull planes can be used to approximately represent the curved surface equation of the curved surface. The plane equation is a linear function, so the difficulty of solving can be greatly reduced.

[0057] Specifically, the above constructing a plurality of convex hull planes and determining a plane equation of each convex hull plane comprises:

[0058] The plane equation of each convex hull plane is constructed, as shown in expression (1).

[0059] (1)

[0060] wherein, V represents the discrete reservoir capacity, represents the convex hull plane composed of the intersection point of the mth and (m+1)th discrete power generation flow and the nth and (n+1)th discrete reservoir capacity, represents the plane equation of the convex hull plane , m=1,…,M-1, n=1,…,N-1, M is the number of discrete power generation flows, N is the number of discrete reservoir capacities, , , is a plane parameter to be solved.

[0061] Taking the following expression (2) as the objective function, and taking the following expressions (3) and (4) as the constraints, the plane parameters of each convex hull plane are solved.

[0062] (2)

[0063] (3)

[0064]

[0065] wherein, represents a theoretical power generation fitting value calculated based on the mth discrete power generation flow and the nth discrete reservoir capacity, represents a discrete power generation fitting value of a plane equation of a convex hull plane adjacent to the convex hull plane, represents a discrete power generation fitting value of a plane equation of a convex hull plane at a center point, represents a discrete power generation fitting value of a plane equation of a convex hull plane at a center point, represents an intersection point formed by the mth discrete power generation flow and the nth discrete reservoir capacity, represents a discrete power generation fitting value at the intersection point , , , , represents four convex hull planes connected to the intersection point , , , , represents a discrete power generation fitting value of a plane equation of the four convex hull planes at the intersection point .

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

[0067] The meaning of the above expression (3) is that, at the center point of the plane , the value of the plane needs to be less than the values of other planes. In this way, the overall consistency of the convex hull plane is ensured.

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

[0069] In step S104, the wind-solar power generation sequence is counted with the second time length as the statistical dimension, to obtain the output fluctuation coefficient of the wind-solar power station, and the second time length is greater than or equal to the first time length.

[0070] Specifically, the output fluctuation coefficient represents the fluctuation size of the wind-solar power generation. For any second time length, the standard deviation and the mean of the wind-solar power generation in the second time length can be determined, and based on the standard deviation and the mean, the output fluctuation coefficient of the wind-solar power station in the second time length can be determined.

[0071] In this embodiment, the second time length is 1 week. Based on the wind power generation sequence and the photovoltaic power generation sequence obtained in step S101, the wind power generation and the photovoltaic power generation at the corresponding time point can be added to obtain the wind-solar power generation at the corresponding time point. For example, the wind power generation and the photovoltaic power generation at the first hour are added to obtain the wind-solar power generation at the first hour; the wind power generation and the photovoltaic power generation at the second hour are added to obtain the wind-solar power generation at the second hour.

[0072] Based on the wind-solar power generation of each hour, the standard deviation and the mean of the wind-solar power generation in the tth week can be counted. Based on the standard deviation and the mean of the wind-solar power generation in the tth week, the output fluctuation coefficient of the wind-solar power station in the tth week can be determined. Wherein, t is an integer between 1 and 36 (including 1 and 36). Specifically, the output fluctuation coefficient of the wind-solar power station in the tth week can be calculated based on expression (5).

[0073] (5)

[0074] wherein, represents the output fluctuation coefficient of the wind-solar power station in the tth week, represents the standard deviation of the wind-solar power generation in the tth week, represents the mean of the wind-solar power generation in the tth week, represents the nth wind-solar power generation in the tth week, represents the number of wind-solar power generations in the tth week, and n is an integer between 1 and T3 (including 1 and T3).

[0075] After the output fluctuation coefficients of each week are normalized, the output fluctuation coefficient sequence in a typical year can be obtained .

[0076] It can be understood that in actual application, the second time length can be determined according to actual needs. For example, the second time length can also be 1 month. The specific value of the second time length is not limited in the present application.

[0077] Step S105, taking the second time length as the statistical dimension and taking the maximum power transmission capacity supported by the power transmission channel as the constraint condition, the wind-solar power generation sequence is counted to obtain the curtailment coefficient of the wind-solar power station.

[0078] Specifically, the curtailment coefficient includes a value representing the curtailment amount of the wind-solar power station. For any second time length, the sum of the wind-solar power generation amount in the second time length can be determined, and based on the sum of the wind-solar power generation amount and the maximum power transmission amount supported by the power transmission channel, the curtailment coefficient of the wind-solar power station in the second time length can be determined.

[0079] In this embodiment, the curtailment coefficient of the wind-solar power station in the t th decade can be calculated based on expression (6).

[0080] (6)

[0081] wherein, represents the curtailment coefficient of the wind-solar power station in the t th decade, represents the curtailment amount of the wind-solar power station in the t th decade, represents the wind-solar power generation amount of the wind-solar power station in the t th decade.

[0082] After the curtailment coefficients of the decades are normalized, the curtailment coefficient sequence of the wind-solar power station in a typical year can be obtained as .

[0083] In step S106, the output weight of the hydropower station is determined based on the output fluctuation coefficient and the curtailment coefficient of the wind-solar power station.

[0084] Specifically, the output weight is used to correct the power generation amount of the hydropower station. The corrected power generation amount is less than the uncorrected power generation amount. The greater the output weight, the greater the correction degree of the power generation amount. That is, the greater the output weight, the smaller the corrected power generation amount.

[0085] The output fluctuation coefficient and the curtailment coefficient of the wind-solar power station can be proportional to the output weight of the hydropower station. The reasons for the proportionality are described below.

[0086] 1) When the output fluctuation coefficient is relatively large, it indicates that the power generation amount of the wind-solar power station has a relatively large value. At this time, the power generation amount of the hydropower station can be further reduced (i.e., the output weight needs to be larger). In this way, when the power generation amount of the wind-solar power station is high, the power generation amount of the hydropower station is also high, thereby causing excessive curtailment amount and waste of natural resources. Of course, it can be understood that when the output fluctuation coefficient is relatively large, the power generation amount of the wind-solar power station also has a relatively small value. When the power generation amount of the wind-solar power station is small, in order to ensure normal power supply, other power stations (such as thermal power stations) in the power system can be controlled to generate power. In this way, on the one hand, the power grid can be ensured to supply power normally, and on the other hand, the natural energy such as wind-solar-hydropower can be fully utilized to avoid waste of natural energy.

[0087] 2) when the abandoned electricity coefficient is large, it indicates that the power generation of the wind-solar power station is already excessive, at this time, the power generation of the hydraulic power station can be further reduced (i.e. the output weight needs to be larger), so as to avoid the power generation of the hydraulic power station being high when the power generation of the wind-solar power station is excessive, thereby causing excessive abandoned electricity and waste of natural resources.

[0088] In the embodiment, for any second time length, the output weight of the wind-solar power station in the second time length can be determined based on the output fluctuation coefficient and the abandoned electricity coefficient of the wind-solar power station in the second time length. That is, based on the output fluctuation coefficient and the abandoned electricity coefficient of the wind-solar power station in the t th decade, the output weight of the wind-solar power station in the t th decade can be determined, so that the output weight sequence of the hydraulic power station in the typical year can be obtained .

[0089] Based on the output weight in the t th decade, the power generation of the hydraulic power station in the t th decade can be corrected.

[0090] Further, the output fluctuation coefficient and the abandoned electricity coefficient of the wind-solar power station can have respective corresponding weights. The weight is used to represent the influence of the output fluctuation coefficient and the abandoned electricity coefficient of the wind-solar power station on the output weight of the hydraulic power station. When the output weight of the hydraulic power station is determined based on the output fluctuation coefficient and the abandoned electricity coefficient of the wind-solar power station, the first weight of the output fluctuation coefficient and the second weight of the abandoned electricity coefficient can be obtained. Based on the first weight and the second weight, the output fluctuation coefficient and the abandoned electricity coefficient of the wind-solar power station can be weighted and fused to obtain the output weight.

[0091] Specifically, the output weight of the hydraulic power station in the t th decade can be determined based on expression (7) .

[0092] (7)

[0093] wherein, represents the first weight, represents the second weight, .

[0094] In step S107, the operation of the hydraulic power station is simulated based on the linear equation and the operation constraint information of the hydraulic power station, and the power generation of the hydraulic power station is corrected based on the output weight, and the operation of the wind-solar power station is simulated based on the wind-solar power information.

[0095] Specifically, when the operation of the wind-solar power station is simulated based on the wind-solar power information, the wind-solar power model actually simulates 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 described here.

[0096] Further, based on the linear equations, the operation of the hydropower station is simulated, including

[0097] obtaining the power generation flow and the reservoir capacity for simulating the operation of the hydropower station;

[0098] substituting the power generation flow and the reservoir capacity into the plane equation of each convex hull plane respectively to obtain a plurality of power generations;

[0099] selecting the minimum power generation as the power generation of the hydropower station from the plurality of power generations.

[0100] Specifically, the process can be represented by expression (8):

[0101] (8)

[0102] wherein, represents the power generation of the hydropower station at the t th time, a, b, and c represent the plane parameters of each convex hull plane, represents the power generation flow of the hydropower station at the t th time, represents the reservoir capacity of the hydropower station at the t th time, represents the reservoir capacity of the i th sub-hydropower station at the t+1 th time.

[0103] The meaning of expression (8) is that the power generation flow and the reservoir capacity are substituted into the plane equation of each convex hull plane respectively to obtain a plurality of power generations; the minimum power generation is selected as the power generation of the hydropower station from the plurality of power generations. In this way, in the case that the linear equation has errors and leads to a high calculation result, the actual power generation of the hydropower station can be effectively prevented from exceeding its actual capacity.

[0104] To sum up, in the technical scheme of some embodiments of the present application, on the one hand, by splitting the water and electricity output function of the hydropower station into a plurality of linear equations, the quadratic function problem can be converted into a linear function problem, and then in the process of simulating the operation of the comprehensive energy base based on the water and electricity output function, the solving difficulty can be greatly reduced, and then the simulation difficulty of the comprehensive energy base can be reduced. On the other hand, based on the output fluctuation coefficient and the power abandonment coefficient of the wind and light power station, the output weight of the hydropower station is determined, and the power generation of the hydropower station is corrected based on the output weight. In this way, the final power generation of the hydropower station can be more accurate, and then the problem of poor compensation adjustment performance of the hydropower station to the wind and light power station can be solved.

[0105] Further, the plurality of convex hull planes of the present application can only approximate the original curved surface constituted by the discrete reservoir capacity, the discrete power generation flow and the discrete power generation power, but cannot be completely coincident with the original curved surface, thus the power generation power calculated based on the plane equation of the convex hull plane and the theoretical power generation power of the hydropower station should be with errors. In order to improve the simulation accuracy, in some embodiments, in the case that the grid-connected power of the hydropower station and the wind-solar power station reaches the maximum transmission power allowed by the transmission channel and the integrated energy base exists power abandonment, if the power generation power of the hydropower station at any target time period obtained based on the plane equation of the convex hull plane is greater than the theoretical power generation power of the hydropower station at the target time period, it indicates that the power generation capacity of the hydropower station is overestimated, at this time, the power generation power of the hydropower station and the power abandonment power of the wind-solar power station can be reduced. Conversely, in the case that the grid-connected power of the hydropower station and the wind-solar power station reaches the maximum transmission power allowed by the transmission channel and the integrated energy base exists power abandonment, if the power generation power of the hydropower station at any target time period obtained based on the plane equation of the convex hull plane is less than the theoretical power generation power of the hydropower station at the target time period, it indicates that the power generation capacity of the hydropower station is underestimated, at this time, the power generation power of the hydropower station and the power abandonment power of the wind-solar power station can be increased, thus the simulation accuracy can be ensured when simulating the operation of the hydropower station and the wind-solar power station.

[0106] In some embodiments, the hydropower station comprises a plurality of sub-hydropower stations, each of which has a respective corresponding power generation flow, reservoir capacity and convex hull plane; the method of the present application further comprises:

[0107] With expression (9) as a constraint condition, the power generation flow and the 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 power of the i-th sub-hydropower station at a plurality of time points.

[0108] (9)

[0109] wherein, Pti represents the power generation power of the i-th sub-hydropower station at the t-th time point, a, b, c represent the plane parameters of the convex hull plane corresponding to the i-th sub-hydropower station, Qti represents the power generation flow of the i-th sub-hydropower station at the t-th time point, Vti represents the reservoir capacity of the i-th sub-hydropower station at the t-th time point, Vti+1 represents the reservoir capacity of the i-th sub-hydropower station at the t+1-th time point.

[0110] Specifically, when the hydropower station comprises a plurality of sub-hydropower stations, the plurality of sub-hydropower stations can be cascade hydropower stations. The complete process of the method of the present application is described below taking the cascade hydropower stations as an example.

[0111] 1) Based on the above steps S101~S103, the water power output function of each sub-hydropower station is respectively split into a plurality of linear equations.

[0112] 2) Based on the linear equations and operation constraint information of each sub-hydropower station, the operation of the sub-hydropower station is simulated, and based on the wind and light power generation information, the operation of the wind and light power station is simulated, to obtain the power generation power and abandoned power of the sub-hydropower station and the wind and light power station at a plurality of time points.

[0113] 3) Based on the power generation power and theoretical power generation power of each sub-hydropower station obtained by solving, the abandoned power of at least part of the sub-hydropower stations obtained by solving is corrected.

[0114] 4) In the case that the total grid-connected power of the water, wind and light power station calculated in the target time period is less than the maximum power allowed by the power transmission channel , If the power generation power 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 power of the hydropower station in the target time period, the total power generation power of the water, wind and light power station in the target time period can be recalculated; if the total grid-connected power of the water, wind and light power station calculated is greater than the maximum power allowed by the power transmission channel, the total abandoned power of the comprehensive energy base in the target time period is corrected to - ; if the power generation power of the hydropower station obtained based on the plane equation of the convex hull plane in the target time period is greater than the theoretical power generation power of the hydropower station in the target time period, the power generation power of the water, wind and light power station in the target time period can be corrected.

[0115] Further, in some embodiments, when simulating the operation of the hydropower station and the wind and light power station, the method of the present application can further comprise: taking a third length as a first scheduling period, taking a second length as a first scheduling period, taking the maximum on-grid power of the comprehensive energy base as the target, taking at least the operation constraint information of the hydropower station as the constraint, correcting the power generation of the hydropower station based on the output weight, obtaining the target water level of the hydropower station at the beginning and end of each first scheduling period, the third length being greater than or equal to the second length, and taking the second length as a second scheduling period, taking the first length as a second scheduling period, taking the target water level of the hydropower station at the beginning and end of each first scheduling period as the constraint, determining the power generation of the hydropower station in each second scheduling period.

[0116] Specifically, the third time length can be a year, the second time length can be a decade, and the first time length can be an hour. This simulation method can also be referred to as a multi-time scale hierarchical nested simulation. Thus, the number of parameters to be solved can be reduced and the difficulty of solving can be reduced in the calculation process of each simulation. For example, if a year is directly used as a scheduling period and an hour is directly used as a scheduling period, the power generation and power abandonment of the wind-solar-hydro power station in 8760 hours need to be solved, and the number of parameters to be solved is too large, which can cause a problem of being unable to solve. If the method based on the multi-time scale hierarchical nested simulation is used, the number of variables to be solved can be effectively reduced. In simple terms, if 8760 hours of continuous simulation is performed, at least 8760 variables (i.e., the power generation of the hydro power station in 8760 hours) need to be solved, and the solving is extremely difficult. However, when the third time length is used as the first scheduling period and the second time length is used as the first scheduling period, the target water level of the hydro power station at the beginning and end of each decade is solved, and the number of variables to be solved is significantly reduced. Similarly, when the second time length is used as the second scheduling period and the third time length is used as the second scheduling period, the power generation of the hydro power station in each hour of each decade is solved, and since the number of hours in a decade is much lower than the number of hours in a year, the number of variables to be solved is also significantly reduced. Thus, the number of variables to be solved can be greatly reduced.

[0117] In some embodiments, when the third time length is used as the first scheduling period and the second time length is used as the first scheduling period, expression (10) can be used as the first objective function, and at least expressions (11)-(17) can be used as constraint conditions to solve the target water level of the hydro power station at the beginning and end of each first scheduling period. It should be noted that, in addition to expressions (11)-(17), other constraint conditions can also be included, such as the sum of the grid-connected powers of the wind-solar-hydro power stations cannot exceed the maximum transmission power allowed by the transmission channel. The present application does not limit the constraint conditions.

[0118] (10)

[0119] wherein, is the grid-connected power of the integrated energy base, T1 is the number of decades in a year (by default, 36), and N1 is the number of sub-hydro power stations, is the power generation of the nth sub-hydro power station in the tth period, represents the total power generation of the n sub-hydro power stations in the tth first scheduling period, is the output weight of the hydro power station, is the grid-connected power of the wind power station in the tth first scheduling period, is the grid-connected power of the photovoltaic power station in the tth first scheduling period; is the power abandonment of the wind-solar power station in the tth first scheduling 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.

[0120] (11)

[0121] (12)

[0122] (13)

[0123] (14)

[0124] (15)

[0125] (16)

[0126] (17)

[0127] 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.

[0128] 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.

[0129] 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.

[0130] In the above expression (14), , respectively represent the minimum and maximum discharge flow of the nth sub-hydropower station in the time period

[0131] In the above expression (15), represents the power generation of the nth sub-hydropower station in the time period , respectively represent the minimum and maximum power generation allowed for the nth sub-hydropower station in the time period

[0132] In the above expression (16), represents the total power generation of the cascade hydropower station in the time period t.

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

[0134] In some embodiments, when the second scheduling period is the second time length and the second scheduling time period is the first time length, expression (18) can be taken as the second objective function, and expressions (19)-(20) can be taken as the constraint conditions.

[0135] (18)

[0136] wherein, represents the online power generation of the comprehensive energy base in one ten-day period, represents the number of hours in one ten-day period, is the total grid-connected power of the water, wind and solar power station in the time period t; is the total power generation of the cascade hydropower station in the time period t; is the grid-connected power of the wind power station in the time period t; is the grid-connected power of the solar power station in the time period t; is the curtailed power of the wind and solar power station in the time period t.

[0137] (19)

[0138] (20)

[0139] In the above expression (20), is the minimum power allowed by the power transmission channel of the comprehensive energy base, is the maximum power transmission allowed by the power transmission channel of the comprehensive energy base. ​​​

[0140] Through the description of the above embodiments, those skilled in the art can clearly understand that the method according to the above embodiments can be realized by means of software on a general hardware platform as necessary, and of course can also be realized by hardware, but in many cases the former is a better embodiment.

[0141] Corresponding to the method, the application also provides an operation simulation device of an integrated energy base. For reference Figure 4 The module schematic diagram of the operation simulation device is provided for some embodiments of the application. Figure 4 In the application, the operation simulation device comprises:

[0142] The information acquisition module 401 is configured to acquire hydropower generation information and wind-solar generation information, the hydropower generation information comprising operation constraint information, maximum reservoir capacity, minimum reservoir capacity, maximum power generation flow rate and minimum power generation flow rate of a hydropower station, and the wind-solar generation information comprising a wind-solar power generation sequence obtained by taking a first time length as a statistical dimension and a maximum power transmission capacity supported by a power transmission channel;

[0143] The interval division module 402 is configured to divide a plurality of discrete reservoir capacities in the interval of the maximum reservoir capacity and the minimum reservoir capacity, divide a plurality of discrete power generation flow rates in the interval of the maximum power generation flow rate and the minimum power generation flow rate, and determine a plurality of discrete power generation powers of the hydropower station based on the discrete reservoir capacities and the discrete power generation flow rates;

[0144] The equation construction module 403 is configured to divide a hydropower output function of the hydropower station into a plurality of linear equations based on the discrete reservoir capacities, the discrete power generation flow rates and the discrete power generation powers;

[0145] The fluctuation coefficient determination module 404 is configured to take a second time length as a statistical dimension to statistically process the wind-solar power generation sequence, and obtain an output fluctuation coefficient of the wind-solar power station, the second time length being greater than or equal to the first time length;

[0146] The curtailment coefficient determination module 405 is configured to take the second time length as a statistical dimension, take the maximum power transmission capacity supported by the power transmission channel as a constraint condition, and statistically process the wind-solar power generation sequence to obtain a curtailment coefficient of the wind-solar power station;

[0147] The weight determination module 406 is configured to determine an output weight of the hydropower station based on the output fluctuation coefficient and the curtailment coefficient of the wind-solar power station;

[0148] The simulation module 407 is configured to simulate the operation of the hydropower station based on the linear equations and the operation constraint information of the hydropower station, correct the power generation capacity of the hydropower station based on the output weight, and simulate the operation of the wind-solar power station based on the wind-solar generation information.

[0149] For reference Figure 5The embodiment of the present application further provides an electronic device, comprising a memory 10 and a processor 20, 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-mentioned embodiments of the operation simulation method of the integrated energy base.

[0150] The embodiment of the present application further provides a computer readable storage medium, which stores a computer program, wherein the computer program is configured to perform the steps in any of the above-mentioned embodiments of the operation simulation method of the integrated energy base when running.

[0151] In an example embodiment, the above-mentioned computer readable storage medium can include, but is not limited to, a U disk, a read-only memory (ROM), a random access memory (RAM), a mobile hard disk, a magnetic disk or an optical disk, and various media that can store computer programs.

[0152] The embodiment of the present application further provides a computer program product, which comprises a computer program, and the computer program is executed by a processor to implement the steps in any of the above-mentioned embodiments of the operation simulation method of the integrated energy base.

[0153] The embodiment of the present application further provides another computer program product, which comprises a non-volatile computer readable storage medium, and the non-volatile computer readable storage medium stores a computer program, and the computer program is executed by a processor to implement the steps in any of the above-mentioned embodiments of the operation simulation method of the integrated energy base.

[0154] The skilled person can further realize that the units and algorithm steps of the examples described in combination with the embodiments disclosed herein can be realized by electronic hardware, computer software or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the components and steps of the examples have been described in the above description in general terms. Whether the functions are realized in hardware or software depends on the specific application and design constraints of the technical solution. The skilled person can use different methods to realize the described functions for each specific application, but such implementation should not be considered beyond the scope of the present application.

[0155] The operation simulation method, device and equipment of the comprehensive energy base provided by the application are described in detail. The principles and implementation modes of the application are described by applying specific examples. The above description of the embodiments is only used to help understand the method of the application and its core idea. It should be pointed out that, for ordinary skilled persons in the technical field, some improvements and modifications can be made to the application without departing from the principles of the application, and these improvements and modifications also fall within the protection scope of the claims of the application.

Claims

1. A method of operating simulation of an integrated energy hub, characterized by, The integrated energy base comprises a hydropower station and a wind-solar power station, and the method comprises: obtaining hydropower information and wind-solar power information, wherein the hydropower information comprises operation constraint information, maximum reservoir capacity, minimum reservoir capacity, maximum power generation flow and minimum power generation flow of the hydropower station, and the wind-solar power information comprises a wind-solar power generation sequence obtained by taking a first time length as a statistical dimension and a maximum power transmission capacity supported by a power transmission channel; dividing a plurality of discrete reservoir capacities in the interval of the maximum reservoir capacity and the minimum reservoir capacity, and dividing a plurality of discrete power generation flows in the interval of the maximum power generation flow and the minimum power generation flow, and determining a plurality of discrete power generation powers of the hydropower station based on the discrete reservoir capacities and the discrete power generation flows; dividing a hydropower output function of the hydropower station into a plurality of linear equations based on the discrete reservoir capacities, the discrete power generation flows and the discrete power generation powers; obtaining an output fluctuation coefficient of the wind-solar power station by taking a second time length as a statistical dimension to statistically process the wind-solar power generation sequence, wherein the second time length is greater than or equal to the first time length; obtaining an abandoned power coefficient of the wind-solar power station by taking the second time length as a statistical dimension and taking the maximum power transmission capacity supported by the power transmission channel as a constraint condition to statistically process the wind-solar power generation sequence; determining an output weight of the hydropower station based on the output fluctuation coefficient and the abandoned power coefficient of the wind-solar power station; simulating the operation of the hydropower station based on the linear equations and the operation constraint information of the hydropower station, correcting the power generation capacity of the hydropower station based on the output weight, and simulating the operation of the wind-solar power station based on the wind-solar power information.

2. The method of claim 1, wherein, The method of dividing the hydropower output function of the hydropower station into a plurality of linear equations based on the discrete reservoir capacities, the discrete power generation flows and the discrete power generation powers comprises: taking the discrete reservoir capacities, the discrete power generation flows and the discrete power generation powers as discrete points on three coordinate axes in a three-dimensional coordinate system, and determining a plurality of intersection points of the discrete reservoir capacities, the discrete power generation flows and the discrete power generation powers in the three-dimensional coordinate system; constructing a plurality of convex hull planes based on the plurality of intersection points, and determining a plane equation of each convex hull plane; taking the plane equation as a plurality of linear equations of the hydropower output function.

3. The method of claim 2, wherein, The method of constructing a plurality of convex hull planes and determining a plane equation of each convex hull plane comprises: constructing a plane equation of each of the convex hull planes wherein, represents a discrete power generation flow, V represents a discrete reservoir capacity, represents a convex hull plane formed by the intersection of the mth, m+1th discrete power generation flow and the nth, n+1th discrete reservoir capacity, represents a plane equation of the convex hull plane , m = 1, …, M-1, n = 1, …, N-1, M is the number of discrete power generation flows, and N is the number of discrete reservoir capacities, , , are plane parameters to be solved. With the expression As the objective function, with the following expression as a constraint, the plane parameters of each convex hull plane are solved: wherein, represents the theoretical power generation fitting value calculated based on the mth discrete power generation flow and the nth discrete reservoir capacity, represents the plane equation of the convex hull plane adjacent to the convex hull plane, represents the discrete power generation fitting value of the plane equation of the convex hull plane at the center point, represents the discrete power generation fitting value of the plane equation of the convex hull plane at the center point, represents the intersection point formed by the mth discrete power generation flow and the nth discrete reservoir capacity, represents the discrete power generation fitting value at the intersection point , , , , represents the four convex hull planes connected to the intersection point , , , , represents the discrete power generation fitting value of the plane equation of the four convex hull planes at the intersection point .

4. The method of claim 3, wherein, simulating the operation of the hydropower station based on the linear equations comprises: obtaining a power generation flow and a reservoir capacity used for simulating the operation of the hydropower station; substituting the power generation flow and the reservoir capacity into the plane equation of each convex hull plane respectively to obtain a plurality of power generation powers; selecting the smallest power generation power as the power generation power of the hydropower station from the plurality of power generation powers.

5. The method of claim 4, wherein, The integrated energy base further comprises a power transmission channel, and the method further comprises: In a case where the grid-connected power of the hydropower station and the wind-solar power station reaches the maximum transmission power allowed by the transmission channel and there is curtailed power in the integrated energy base, if the power generated by the hydropower station in any target time period based on the plane equation of the convex hull plane is greater than the theoretical power generated by the hydropower station in the target time period, the power generated by the hydropower station and the curtailed power of the wind-solar power station are reduced.

6. The method of claim 4, wherein, The method further comprises: In a case where the grid-connected power of the hydropower station and the wind-solar power station reaches the maximum transmission power allowed by the transmission channel and there is curtailed power in the integrated energy base, if the power generated by the hydropower station in any target time period based on the plane equation of the convex hull plane is less than the theoretical power generated by the hydropower station in the target time period, the power generated by the hydropower station and the curtailed power of the wind-solar power station are increased.

7. The method of claim 1, wherein, In simulating the operation of the hydropower station and the wind-solar power station, the method further comprises: taking a third time length as a first scheduling period, taking the second time length as a first scheduling time period, taking the maximum on-grid power of the integrated energy base as a target, taking at least the operation constraint information of the hydropower station as a constraint, correcting the power generation of the hydropower station based on the output weight, to obtain the target water level of the hydropower station at the beginning and end of each first scheduling time period, the third time length being greater than or equal to the second time length; taking the second time length as a second scheduling period, taking the first time length as a second scheduling time period, taking the target water level of the hydropower station at the beginning and end of each first scheduling time period as a constraint, determining the power generation of the hydropower station in each second scheduling time period.

8. An integrated energy base operation simulation device, characterized by, The integrated energy base comprises a hydropower station and a wind-solar power station, and the device comprises: an information acquisition module configured to acquire hydropower information and wind-solar power information, the hydropower information comprising operation constraint information, a maximum reservoir capacity, a minimum reservoir capacity, a maximum power generation flow rate and a minimum power generation flow rate of the hydropower station, and the wind-solar power information comprising a wind-solar power generation sequence obtained by taking a first time length as a statistical dimension and a maximum transmission power supported by a transmission channel; an interval division module configured to divide a plurality of discrete reservoir capacities in an interval of the maximum reservoir capacity and the minimum reservoir capacity, and divide a plurality of discrete power generation flow rates in an interval of the maximum power generation flow rate and the minimum power generation flow rate, and determine a plurality of discrete power generations of the hydropower station based on the discrete reservoir capacities and the discrete power generation flow rates; an equation construction module configured to divide a hydropower output function of the hydropower station into a plurality of linear equations based on the discrete reservoir capacities, the discrete power generation flow rates and the discrete power generations; a fluctuation coefficient determination module configured to take a second time length as a statistical dimension to statistically process the wind-solar power generation sequence, to obtain an output fluctuation coefficient of the wind-solar power station, the second time length being greater than or equal to the first time length; and a power generation determination module configured to take the second time length as a second scheduling period, take the first time length as a second scheduling time period, take the target water level of the hydropower station at the beginning and end of each first scheduling time period as a constraint, and determine the power generation of the hydropower station in each second scheduling time period. The abandoned electricity coefficient determination module is configured to take the second time length as a statistical dimension, take the maximum power transmission amount supported by the power transmission channel as a constraint condition, and statistically process the wind-solar power generation amount sequence to obtain the abandoned electricity coefficient of the wind-solar power station. The weight determination module is configured to determine the output weight of the hydroelectric power station based on the output fluctuation coefficient and the abandoned electricity coefficient of the wind-solar power station. The simulation module is configured to simulate the operation of the hydroelectric power station based on the linear equation and the operation constraint information of the hydroelectric power station, correct the power generation amount of the hydroelectric power station based on the output weight, and simulate the operation of the wind-solar power station based on the wind-solar power generation information.

9. An electronic device, comprising: The computer readable storage medium stores a computer program, and the computer program is executed by the processor to implement the operation simulation method of the integrated energy base. The computer readable storage medium stores a computer program, and the computer program is executed by the processor to implement the operation simulation method of the integrated energy base. The computer readable storage medium stores a computer program, and the computer program is executed by the processor to implement the operation simulation method of the integrated energy base.

10. A computer-readable storage medium, characterized in that, ​

Images