Reference power grid planning method, system and medium based on load moment and acceptance domain

Abstract

The present application discloses a reference power grid planning method, system and medium based on load moment and acceptance domain, which mainly relates to the technical field of reference power grid planning, and is used to solve the problem that the prediction interval calculated by the existing traditional model is sometimes too conservative or radical and does not take into account the line length and flow data. It includes: obtaining a preset upper acceptance domain maximum value calculation function, transforming the nonlinear terms in the preset upper acceptance domain maximum value calculation function to obtain a linear expression of the calculation function; obtaining a preset constraint formula corresponding to the linear expression; determining the acceptance domain maximum value corresponding to the linear expression based on the priority target planning method and the preset constraint formula; and then obtaining the acceptance domain value range constraint; obtaining the calculation function of the sum of the preset lower minimum load moments and the acceptance domain value range constraint corresponding to the reference power grid, and then determining the specific value of the sum of the minimum load moments.

Description

Technical Field

[0001] The present application relates to the technical field of reference power grid planning, and in particular to a reference power grid planning method, system and medium based on load moment and acceptance domain. Background Art

[0002] The large-scale grid connection of renewable energy has become an important trend in the development of the power system. Due to the uncertainty factors such as randomness, intermittency, volatility and uncontrollability of renewable energy power generation, the difficulty of transmission network planning has increased. While taking into account the uncertainty factors brought by the integration of renewable energy into the transmission network, it is particularly urgent and important to establish a safe, economical and reliable transmission network planning model.

[0003] Reference grid planning is one of the basic contents in the power system. An important difference from traditional transmission network planning is that reference grid planning does not optimize the capacity of one or several new lines, but is applied to the entire transmission network system, including new and existing lines. The main goal of reference grid planning is to reasonably determine the capacity of transmission lines based on the accurate prediction of load growth during the planning period and the existing transmission network topology, quantify the optimal investment cost and line congestion cost, and compare with the real network to find new investment needs and stranded investments, so as to ensure that the power system develops in a safer, more economical and more reliable direction.

[0004] In order to quantify the power system's ability to accept renewable energy, traditional models use a deterministic prediction interval to describe the penetration rate of renewable energy in the power system. However, the deterministic prediction interval is sometimes too conservative or radical. In addition, when planning lines in actual projects, the goal is always to "power supply nearby" and the objective function is to minimize the sum of line investment cost and operating cost. The line length and flow data are not taken into account, which makes the design engineering value of the reference power grid planning line still imperfect. Summary of the invention

[0005] In view of the above-mentioned deficiencies in the prior art, the present application provides a reference power grid planning method, system and medium based on load moment and acceptance domain to solve the problem that the prediction interval calculated by traditional models is sometimes too conservative or radical and does not take into account the line length and flow data.

[0006] In the first aspect, the present application provides a reference power grid planning method based on load moment and acceptance domain, the method comprising: obtaining a preset upper acceptance domain maximum value calculation function, transforming the nonlinear terms in the preset upper acceptance domain maximum value calculation function to obtain a linear expression of the calculation function; wherein the parameters involved in the preset upper acceptance domain maximum value calculation function include the acceptable domain and active net load of each node; obtaining a preset constraint formula corresponding to the linear expression; wherein the preset constraint formula includes a branch flow constraint formula; based on the priority target planning method and the preset constraint formula, determining the acceptance domain maximum value corresponding to the linear expression; and then obtaining the acceptance domain value range constraint; obtaining the calculation function and acceptance domain value range constraint of the sum of the preset lower minimum load moments corresponding to the reference power grid, and then determining the specific value of the sum of the minimum load moments; wherein the reference power grid planning includes branch lines and planned sections, and the parameters involved in the calculation function include branch lines, planned sections, and branch flow upper limits.

[0007] Further, obtaining a preset upper layer acceptance domain maximum value calculation function specifically includes: obtaining a preset upper layer acceptance domain maximum value calculation function:

[0008] ]; where Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0009] Further, the nonlinear term in the preset upper acceptance domain maximum calculation function is transformed to obtain a linear expression of the calculation function, specifically including: on the basis of the preset upper acceptance domain maximum calculation function, two non-negative continuous variables and variable constraints are introduced, and the preset upper acceptance domain maximum calculation function is corrected to obtain a linear expression of the calculation function; wherein, the introduction of the two non-negative continuous variables is: and ;and , Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the total number of nodes; the variable constraints are: and .

[0010] Furthermore, the preset constraint formula includes at least: operating base point power balance constraint, generator set maximum adjustment constraint, generator set preset constraint; obtaining the preset constraint formula corresponding to the linear expression specifically includes: obtaining the operating base point power balance constraint:

[0011] ; Where, f is the f-th generator set, NG is the number of all generator sets; i represents the i-th node; NI represents the number of all nodes; represents the expected value of active power output of the fth generator set, and The value range of is [the minimum value of the net load downward acceptable domain, the maximum value of the net load upward acceptable domain]; represents the expected net load of the ith node, and The value range of is [lower limit of active net load disturbance, upper limit of active net load disturbance]; obtain the maximum adjustment constraint of the generator set:

[0012] and ;in, represents the participation factor, and ; It represents the maximum upward adjustment of the load of the jth generator set in the tth time interval; It represents the maximum downward adjustment of the load of the jth generator set in the tth time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; Represents the lower limit of the active net load disturbance at the i-th node in the t-th time interval; obtain the preset constraints of the generator set:

[0013] ,

[0014] ,

[0015] ;in, Indicates the lower limit of active power generated by the f-th generator set; It represents the upper limit of active power generated by the f-th generator set; Indicates the output active power; represents the upward climbing rate of the fth generator set; represents the downward climbing rate of the fth generator set; represents the duration of the tth time interval; represents the upward reserve rate of the fth generator in the tth time interval; It represents the downward reserve rate of the f-th generator in the t-th time interval.

[0016] Furthermore, the preset constraint formula corresponding to the linear expression is obtained, which specifically includes:

[0017] ,h=1,2,......,NL;

[0018] ,h=1,2,......,NL; where, represents the net load disturbance of the ith node, and i=1,2,...,NI; represents the expected value of the net load of the ith node; h represents the hth branch line; NL represents the total number of branch lines; represents the power transfer distribution factor of the net load of the ith node to the hth branch line; represents the power transfer distribution factor of the net load of the f-th generator set to the h-th branch line; represents the upper limit of the branch power flow of the h-th branch line.

[0019] Further, obtaining the acceptance domain value range constraint specifically includes: obtaining the acceptance domain value range constraint based on the acceptance domain maximum value corresponding to the linear expression:

[0020] ,in, represents the maximum value of the acceptance domain, t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0021] Further, obtaining a calculation function of the sum of the preset lower layer minimum load moments corresponding to the reference power grid specifically includes: obtaining a calculation function of the sum of the preset lower layer minimum load moments corresponding to the reference power grid: , where s represents the sth planned road section, NS represents the number of all planned road sections; h represents the hth branch road; NL represents the total number of branch roads; represents the length of the sth planned road section, represents the upper limit of the branch power flow of the h-th branch line.

[0022] Furthermore, the method also includes: calculating the specific numerical values ​​of the maximum value of the acceptance domain and the sum of the minimum load moments corresponding to the linear expression through a CPLEX solver.

[0023] In the second aspect, the present application provides a reference power grid planning system based on load moment and acceptance domain, and the system includes: a linear expression acquisition module, which is used to obtain a preset upper acceptance domain maximum value calculation function, and transform the nonlinear terms in the preset upper acceptance domain maximum value calculation function to obtain a linear expression of the calculation function; wherein, the parameters involved in the preset upper acceptance domain maximum value calculation function include the acceptable domain and active net load of each node; an acceptance domain value range acquisition module, which is used to obtain a preset constraint formula corresponding to the linear expression; wherein the preset constraint formula includes a branch flow constraint formula; based on the priority target planning method and the preset constraint formula, determine the acceptance domain maximum value corresponding to the linear expression; and then obtain the acceptance domain value range constraint; a minimum load moment sum calculation module, which is used to obtain the calculation function and acceptance domain value range constraint of the preset lower minimum load moment sum corresponding to the reference power grid, and then determine the specific value of the minimum load moment sum; wherein, the reference power grid planning includes branch lines and planned sections, and the parameters involved in the calculation function include branch lines, planned sections, and branch flow upper limits.

[0024] Furthermore, the linear expression acquisition module includes a function acquisition unit, which is used to obtain a preset upper layer acceptance domain maximum value calculation function: ]; where Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0025] Furthermore, the linear expression acquisition module includes an expression acquisition unit, which is used to introduce two non-negative continuous variables and variable constraints on the basis of a preset upper layer acceptance domain maximum calculation function, and to modify the preset upper layer acceptance domain maximum calculation function to obtain a linear expression of the calculation function; wherein the two non-negative continuous variables are introduced as follows: and ;and , Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the total number of nodes; the variable constraints are:

[0026] and .

[0027] Furthermore, the preset constraint formula includes at least: operating base power balance constraint, generator set maximum adjustment constraint, generator set preset constraint; the acceptance domain value range acquisition module includes a preset acquisition unit for acquiring the operating base power balance constraint: ; Where, f is the f-th generator set, NG is the number of all generator sets; i represents the i-th node; NI represents the number of all nodes; represents the expected value of active power output of the fth generator set, and The value range of is [the minimum value of the net load downward acceptable domain, the maximum value of the net load upward acceptable domain]; represents the expected net load of the ith node, and The value range of is [lower limit of active net load disturbance, upper limit of active net load disturbance]; obtain the maximum adjustment constraint of the generator set: and ;in, represents the participation factor, and ; It represents the maximum upward adjustment of the load of the jth generator set in the tth time interval; It represents the maximum downward adjustment of the load of the jth generator set in the tth time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; Represents the lower limit of the active net load disturbance at the i-th node in the t-th time interval; obtain the preset constraints of the generator set:

[0028] , , ;in, Indicates the lower limit of active power generated by the f-th generator set; It represents the upper limit of active power generated by the f-th generator set; Indicates the output active power; represents the upward climbing rate of the fth generator set; represents the downward climbing rate of the fth generator set; represents the duration of the tth time interval; represents the upward reserve rate of the fth generator in the tth time interval; It represents the downward reserve rate of the f-th generator in the t-th time interval.

[0029] Furthermore, the acceptance domain value range acquisition module includes a branch power flow constraint acquisition unit, which is used to obtain the branch power flow constraint formula:

[0030] ,h=1,2,......,NL;

[0031] ,h=1,2,......,NL; where, represents the net load disturbance of the ith node, and i=1,2,...,NI; represents the expected value of the net load of the ith node; h represents the hth branch line; NL represents the total number of branch lines; represents the power transfer distribution factor of the net load of the ith node to the hth branch line; represents the power transfer distribution factor of the net load of the f-th generator set to the h-th branch line; represents the upper limit of the branch power flow of the h-th branch line.

[0032] Furthermore, the acceptance domain value range obtaining module includes a range determination unit, which is used to obtain the acceptance domain value range constraint based on the acceptance domain maximum value corresponding to the linear expression:

[0033] ,in, represents the maximum value of the acceptance domain, t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0034] Furthermore, the minimum load moment sum calculation module includes a calculation function obtaining unit, which is used to obtain a calculation function of the preset lower layer minimum load moment sum corresponding to the reference power grid: , where s represents the sth planned road section, NS represents the number of all planned road sections; h represents the hth branch road; NL represents the total number of branch roads; represents the length of the sth planned road section, represents the upper limit of the branch power flow of the h-th branch line.

[0035] In a third aspect, the present application provides a non-volatile computer storage medium having computer instructions stored thereon, which when executed implement a reference power grid planning method based on load moment and acceptance domain as described above.

[0036] Those skilled in the art can understand that the present application has at least the following beneficial effects:

[0037] This application divides the reference grid planning of the transmission network into two levels based on the method of priority target planning and effective safety domain analysis: first, the output of the automatic generation control unit is adjusted through the decision (preset constraint formula) of the upper model (preset upper acceptance domain maximum value calculation function), that is, the maximum acceptable domain of new energy is obtained through the automatic generation control unit, which makes up for the problem of insufficient absorption capacity of the external new energy after the grid connection in the lower model (preset calculation function of the sum of the minimum load moment of the lower layer), and at the same time, the output of the generator at each node is determined to reduce the uncertainty in the lower model; secondly, the lower model (preset calculation function of the sum of the minimum load moment of the lower layer) takes the minimum value of the sum of the line capacity and the line length product planned by the reference grid as the target, that is, the load moment is minimized as the objective function. The priority target planning method is used to solve the upper and lower models, so that the upper and lower models have a very clear order of precedence. The maximum acceptance domain calculated in the upper model will become a constraint condition in the lower model after relaxation. After solving, the reference grid planning result of the maximum absorption domain of new energy will be obtained. The upper and lower models reflect the interaction and mutual influence of the upper and lower decisions by adding constraints.

[0038] Based on the reference grid planning of the maximum acceptance domain of new energy, this application starts from the perspective of optimizing the economic efficiency of power system operation, while taking into account a certain new energy acceptance capacity, and integrates the robust optimization method with the actual planning method to obtain a reference grid planning method for optimal consumption of new energy. This application adopts the concept of load moment, while considering power supply stability and line loss, reducing the number of variables, being more in line with engineering practice and reducing computational complexity, and being closer to real life. That is, it solves the problem that the prediction interval calculated by the traditional model is sometimes too conservative or radical and does not take into account the line length and flow data. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] Some embodiments of the present disclosure are described below with reference to the accompanying drawings, in which:

[0040] Figure 1 It is a flow chart of a reference power grid planning method based on load moment and acceptance domain provided in an embodiment of the present application.

[0041] Figure 2It is a schematic diagram of the internal structure of a reference power grid planning system based on load moment and acceptance domain provided in an embodiment of the present application. DETAILED DESCRIPTION

[0042] It should be understood by those skilled in the art that the embodiments described below are only preferred embodiments of the present disclosure, and do not mean that the present disclosure can only be implemented through the preferred embodiments. The preferred embodiments are only used to explain the technical principles of the present disclosure, and are not used to limit the protection scope of the present disclosure. Based on the preferred embodiments provided by the present disclosure, all other embodiments obtained by ordinary technicians in this field without creative work should still fall within the protection scope of the present disclosure.

[0043] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0044] The technical solution proposed in the embodiments of the present application is described in detail below with reference to the accompanying drawings.

[0045] The present application embodiment provides a reference power grid planning method based on load moment and acceptance domain, such as Figure 1 As shown, the method provided in the embodiment of the present application mainly includes the following steps:

[0046] Step 110: Obtain a preset upper layer acceptance domain maximum value calculation function, and transform the nonlinear terms in the preset upper layer acceptance domain maximum value calculation function to obtain a linear expression of the calculation function.

[0047] It should be noted that the parameters involved in the preset upper acceptance domain maximum value calculation function include the acceptable domain and active net load of each node. The acceptable domain is the acceptance capacity of the reference power grid system for uncertain disturbances, which is an optimizable interval quantity and changes with the change of the uncertainty fluctuation size in the reference power grid system. The unit of the acceptable domain can be a power capacity unit.

[0048] The step of obtaining the preset upper layer acceptance domain maximum value calculation function specifically includes:

[0049] Get the preset upper acceptance domain maximum calculation function:

[0050] ]; where Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0051] The nonlinear term in the preset upper acceptance domain maximum value calculation function is converted to obtain a linear expression of the calculation function, which can be specifically:

[0052] On the basis of the preset upper layer acceptance domain maximum calculation function, two non-negative continuous variables and variable constraints are introduced, and the preset upper layer acceptance domain maximum calculation function is modified to obtain the linear expression of the calculation function; wherein, the two non-negative continuous variables are introduced as follows:

[0053] and ;and , Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the total number of nodes; the variable constraints are:

[0054] and .

[0055] Step 120, obtaining a preset constraint formula corresponding to the linear expression; determining the maximum value of the acceptance domain corresponding to the linear expression based on the priority goal programming method and the preset constraint formula; and then obtaining the acceptance domain value range constraint.

[0056] It should be noted that the preset constraint formula includes branch power flow constraint formula, operating base point power balance constraint, generator set maximum adjustment constraint, and generator set preset constraint.

[0057] As an example, obtaining a preset constraint formula corresponding to a linear expression may specifically include:

[0058] Get the operating base point power balance constraint:

[0059] ; Where, f is the f-th generator set, NG is the number of all generator sets; i represents the i-th node; NI represents the number of all nodes; represents the expected value of active power output of the fth generator set, and The value range of is [the minimum value of the net load downward acceptable domain, the maximum value of the net load upward acceptable domain]; represents the expected value of the net load of the ith node, and The value range of is [lower limit of active net load disturbance, upper limit of active net load disturbance].

[0060] Get the maximum adjustment constraint of the generator set:

[0061] and ;in, represents the participation factor, and ; It represents the maximum upward adjustment of the load of the jth generator set in the tth time interval; It represents the maximum downward adjustment of the load of the jth generator set in the tth time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0062] Get the preset constraints of the generator set:

[0063] , , .

[0064] in, Indicates the lower limit of active power generated by the f-th generator set; It represents the upper limit of active power generated by the f-th generator set; Indicates the output active power; represents the upward climbing rate of the fth generator set; represents the downward climbing rate of the fth generator set; represents the duration of the tth time interval; represents the upward reserve rate of the fth generator in the tth time interval; It represents the downward reserve rate of the f-th generator in the t-th time interval.

[0065] As a second example, obtaining a preset constraint formula corresponding to a linear expression may further include:

[0066] Get the branch power flow constraint formula:

[0067] ,h=1,2,......,NL,

[0068] ,h=1,2,......,NL.

[0069] in, represents the net load disturbance of the ith node, and i=1,2,...,NI; represents the expected value of the net load of the ith node; h represents the hth branch line; NL represents the total number of branch lines; represents the power transfer distribution factor of the net load of the ith node to the hth branch line; represents the power transfer distribution factor of the net load of the f-th generator set to the h-th branch line; represents the upper limit of the branch power flow of the h-th branch line.

[0070] The step of obtaining the acceptance domain value range constraint may specifically be:

[0071] Based on the maximum value of the acceptance domain corresponding to the linear expression, the acceptance domain value range constraint is obtained:

[0072] .

[0073] in, represents the maximum value of the acceptance domain, t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0074] Step 130: Obtain a calculation function of the sum of the preset lower-layer minimum load moments and an acceptance domain value range constraint corresponding to the reference power grid, and then determine a specific value of the sum of the minimum load moments.

[0075] It should be noted that the reference power grid planning includes branch lines and planned sections, and the parameters involved in the calculation function include branch lines, planned sections, and branch power flow upper limits.

[0076] The calculation function for obtaining the sum of the preset lower-layer minimum load moments corresponding to the reference power grid may specifically be:

[0077] Get the calculation function of the sum of the preset lower-layer minimum load moments corresponding to the reference power grid: , where s represents the sth planned road section, NS represents the number of all planned road sections; h represents the hth branch road; NL represents the total number of branch roads; represents the length of the sth planned road section, represents the upper limit of the branch power flow of the h-th branch line.

[0078] In addition, the present application can also solve and calculate the specific values ​​of the maximum value of the acceptance domain and the sum of the minimum load moments corresponding to the linear expression through the CPLEX solver.

[0079] Based on the above description, it can be understood by those skilled in the art that, based on the reference grid planning of the maximum acceptance domain of new energy (acceptance domain maximum value), the present application starts from the perspective of optimizing the economic efficiency of power system operation, while taking into account a certain new energy acceptance capacity, and integrates the robust optimization method with the actual planning method to obtain a reference grid planning method for optimal consumption of new energy. This model retains the advantages of high robust optimization calculation efficiency and ensuring that the operation constraints (preset constraint formulas) within the predetermined disturbance set are satisfied, while also maintaining the statistical superiority of the actual planning results. The upper layer is a decision layer for maximizing the acceptable domain of new energy (preset upper acceptance domain maximum value calculation function), and the decision result is passed to the lower layer (preset lower minimum load moment sum calculation function). The lower layer constrains the system operation cost and line investment cost, and obtains the reference grid based on the maximum acceptable domain constructed by the upper model, so as to optimize the network topology and find stranded investments, thereby improving the power supply capacity of the transmission network to cope with severe fluctuations in net load under high wind power penetration.

[0080] besides, Figure 2 A reference power grid planning system based on load moment and acceptance domain is provided in an embodiment of the present application. Figure 2 As shown, the system provided in the embodiment of the present application mainly includes:

[0081] The linear expression obtaining module 210 is used to obtain a preset upper layer acceptance domain maximum calculation function, and transform the nonlinear terms in the preset upper layer acceptance domain maximum calculation function to obtain a linear expression of the calculation function.

[0082] It should be noted that the linear expression acquisition module 210 can be any feasible device / component / device that can obtain the linear expression of the calculation function. The parameters involved in the preset upper layer acceptance domain maximum value calculation function include the acceptance domain and active net load of each node;

[0083] The step of obtaining the preset upper layer acceptance domain maximum value calculation function specifically includes:

[0084] The function acquisition unit in the linear expression acquisition module 210 obtains the preset upper layer acceptance domain maximum value calculation function:

[0085] ]; where Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0086] The nonlinear term in the preset upper acceptance domain maximum value calculation function is converted to obtain a linear expression of the calculation function, which can be specifically:

[0087] The expression obtaining unit in the linear expression obtaining module 210 introduces two non-negative continuous variables and variable constraints on the basis of the preset upper layer acceptance domain maximum calculation function, and modifies the preset upper layer acceptance domain maximum calculation function to obtain the linear expression of the calculation function; wherein, the two non-negative continuous variables introduced are:

[0088] and ;and , Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the total number of nodes; the variable constraints are:

[0089] and .

[0090] The acceptance domain value range obtaining module 220 is used to obtain a preset constraint formula corresponding to the linear expression; based on the priority goal programming method and the preset constraint formula, determine the maximum value of the acceptance domain corresponding to the linear expression; and then obtain the acceptance domain value range constraint;

[0091] It should be noted that the acceptance domain value range acquisition module 220 can be any feasible device / component / device that can determine the maximum value of the acceptance domain corresponding to the linear expression; and then obtain the acceptance domain value range constraint. Among them, the preset constraint formula includes the branch flow constraint formula. The preset constraint formula at least includes: operating base point power balance constraint, generator set maximum adjustment constraint, and generator set preset constraint.

[0092] As an example, obtaining a preset constraint formula corresponding to a linear expression may specifically include:

[0093] The power balance constraint of the operating base point is obtained by obtaining the preset acquisition unit in the module 220 through the acceptance domain value range: ; Where, f is the f-th generator set, NG is the number of all generator sets; i represents the i-th node; NI represents the number of all nodes; represents the expected value of active power output of the fth generator set, and The value range of is [the minimum value of the net load downward acceptable domain, the maximum value of the net load upward acceptable domain]; represents the expected value of the net load of the ith node, and The value range of is [lower limit of active net load disturbance, upper limit of active net load disturbance].

[0094] Get the maximum adjustment constraint of the generator set:

[0095] and .

[0096] in, represents the participation factor, and ; It represents the maximum upward adjustment of the load of the jth generator set in the tth time interval; It represents the maximum downward adjustment of the load of the jth generator set in the tth time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0097] Get the preset constraints of the generator set:

[0098] , , .

[0099] in, Indicates the lower limit of active power generated by the f-th generator set; It represents the upper limit of active power generated by the f-th generator set; Indicates the output active power; represents the upward climbing rate of the fth generator set; represents the downward climbing rate of the fth generator set; represents the duration of the tth time interval; represents the upward reserve rate of the fth generator in the tth time interval; It represents the downward reserve rate of the f-th generator in the t-th time interval.

[0100] As a second example, obtaining a preset constraint formula corresponding to a linear expression may further include:

[0101] The branch power flow constraint obtaining unit in the module 220 is obtained by accepting the value range of the domain to obtain the branch power flow constraint formula:

[0102] ,h=1,2,......,NL;

[0103] ,h=1,2,......,NL.

[0104] in, represents the net load disturbance of the ith node, and i=1,2,...,NI; represents the expected value of the net load of the ith node; h represents the hth branch line; NL represents the total number of branch lines; represents the power transfer distribution factor of the net load of the ith node to the hth branch line; represents the power transfer distribution factor of the net load of the f-th generator set to the h-th branch line; represents the upper limit of the branch power flow of the h-th branch line.

[0105] The acceptance domain value range constraint is obtained, which can be specifically:

[0106] The range determination unit in the acceptance domain value range acquisition module 220 obtains the acceptance domain value range constraint based on the acceptance domain maximum value corresponding to the linear expression:

[0107] .

[0108] in, represents the maximum value of the acceptance domain, t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

[0109] The minimum load moment sum calculation module 230 is used to obtain the calculation function of the preset lower layer minimum load moment sum corresponding to the reference power grid and the acceptance domain value range constraint, and then determine the specific value of the minimum load moment sum.

[0110] It should be noted that the minimum load moment sum calculation module 230 can be any feasible device / component / device that can determine the specific value of the minimum load moment sum. Among them, the reference power grid planning includes branch lines and planned sections, and the parameters involved in the calculation function include branch lines, planned sections, and branch flow upper limits.

[0111] The calculation function for obtaining the sum of the preset lower-layer minimum load moments corresponding to the reference power grid may specifically be:

[0112] The calculation function obtaining unit in the minimum load moment sum calculation module 230 obtains the calculation function of the preset lower layer minimum load moment sum corresponding to the reference power grid:

[0113] , where s represents the sth planned road section, NS represents the number of all planned road sections; h represents the hth branch road; NL represents the total number of branch roads; represents the length of the sth planned road section, represents the upper limit of the branch power flow of the h-th branch line.

[0114] In addition, an embodiment of the present application further provides a non-volatile computer storage medium having executable instructions stored thereon, and when the executable instructions are executed, a reference power grid planning method based on load moment and acceptance domain as described above is implemented.

[0115] So far, the technical solutions of the present disclosure have been described in combination with the above multiple embodiments, but it is easy for those skilled in the art to understand that the protection scope of the present disclosure is not limited to these specific embodiments. Without departing from the technical principles of the present disclosure, those skilled in the art can split and combine the technical solutions in the above-mentioned various embodiments, and can also make equivalent changes or replacements to the relevant technical features. Any changes, equivalent replacements, improvements, etc. made within the technical concept and / or technical principle of the present disclosure will fall within the protection scope of the present disclosure.

Claims

1. A reference power grid planning method based on load moment and acceptance domain, characterized in that: The method comprises: Obtaining a preset upper-layer acceptance domain maximum value calculation function, transforming the nonlinear terms in the preset upper-layer acceptance domain maximum value calculation function to obtain a linear expression of the calculation function; wherein the parameters involved in the preset upper-layer acceptance domain maximum value calculation function include the acceptance domain and active net load of each node; Obtaining a preset constraint formula corresponding to the linear expression; wherein the preset constraint formula includes a branch power flow constraint formula, and the preset constraint formula includes at least: an operating base point power balance constraint, a generator set maximum adjustment constraint, and a generator set preset constraint; Get the preset constraint formula corresponding to the linear expression, including: Get the operating base point power balance constraint: ; Where, f is the f-th generator set, NG is the number of all generator sets; i is the i-th node; NI is the number of all nodes; represents the expected value of active power output of the fth generator set, and The value range of is from the minimum value of the net load downward acceptable domain to the maximum value of the net load upward acceptable domain; represents the expected value of the net load of the ith node, and The value range of is from the lower limit of active net load disturbance to the upper limit of active net load disturbance; Get the maximum adjustment constraint of the generator set: and ; in, represents the participation factor, and ; It represents the maximum upward adjustment of the load of the jth generator set in the tth time interval; It represents the maximum downward adjustment of the load of the jth generator set in the tth time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; represents the lower limit of the active net load disturbance on the ith node in the tth time interval; Get the preset constraints of the generator set: 、 、 ； in, Indicates the lower limit of active power generated by the f-th generator set; It represents the upper limit of active power generated by the f-th generator set; Indicates the output active power; represents the upward climbing rate of the fth generator set; represents the downward climbing rate of the fth generator set; represents the duration of the tth time interval; represents the upward reserve rate of the fth generator in the tth time interval; represents the downward reserve rate of the fth generator in the tth time interval; Get the preset constraint formula corresponding to the linear expression, including: Get the branch power flow constraint formula: ,h=1,2,......,NL; ,h=1,2,......,NL; in, represents the net load disturbance of the ith node, and i=1,2,...,NI; represents the expected value of the net load of the ith node; h represents the hth branch line; NL represents the total number of branch lines; represents the power transfer distribution factor of the net load of the ith node to the hth branch line; represents the power transfer distribution factor of the net load of the f-th generator set to the h-th branch line; represents the upper limit of the branch power flow of the h-th branch line; Based on the priority goal programming method and the preset constraint formula, the maximum value of the acceptance domain corresponding to the linear expression is determined; then the value range constraint of the acceptance domain is obtained; Among them, obtaining the acceptance domain value range constraints specifically includes: Based on the maximum value of the acceptance domain corresponding to the linear expression, the acceptance domain value range constraint is obtained: ， in, represents the maximum value of the acceptance domain, t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; represents the lower limit of the active net load disturbance on the ith node in the tth time interval; Obtain the calculation function of the sum of the preset lower-layer minimum load moments corresponding to the reference power grid and the acceptance domain value range constraint, and then determine the specific value of the sum of the minimum load moments; wherein the reference power grid planning includes branch lines and planned sections, and the parameters involved in the calculation function include branch lines, planned sections, and branch power flow upper limits; The calculation function for obtaining the sum of the preset lower-layer minimum load moments corresponding to the reference power grid specifically includes: Get the calculation function of the sum of the preset lower-layer minimum load moments corresponding to the reference power grid: , where s represents the sth planned road section, NS represents the number of all planned road sections; h represents the hth branch road; NL represents the total number of branch roads; represents the length of the sth planned road section, represents the upper limit of the branch power flow of the h-th branch line.

2. The reference power grid planning method based on load moment and acceptance domain according to claim 1 is characterized in that: Get the preset upper acceptance domain maximum calculation function, including: Get the preset upper acceptance domain maximum calculation function: ]； Where Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

3. The reference power grid planning method based on load moment and acceptance domain according to claim 1 or 2, characterized in that: The nonlinear terms in the preset upper acceptance domain maximum calculation function are converted to obtain the linear expression of the calculation function, including: On the basis of the preset upper layer acceptance domain maximum value calculation function, two non-negative continuous variables and variable constraints are introduced, and the preset upper layer acceptance domain maximum value calculation function is modified to obtain a linear expression of the calculation function; Among them, two non-negative continuous variables are introduced as follows: and ; and , Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; The constraints on the variables are: and .

4. The reference power grid planning method based on load moment and acceptance domain according to claim 1 is characterized in that: The method further comprises: The CPLEX solver is used to calculate the specific values ​​of the maximum value of the acceptance domain and the sum of the minimum load moments corresponding to the linear expression.

5. A reference power grid planning system based on load moment and acceptance domain, characterized in that: The system comprises: A linear expression acquisition module is used to obtain a preset upper-layer acceptance domain maximum value calculation function, and transform the nonlinear terms in the preset upper-layer acceptance domain maximum value calculation function to obtain a linear expression of the calculation function; wherein the parameters involved in the preset upper-layer acceptance domain maximum value calculation function include the acceptance domain and active net load of each node; The module for obtaining the value range of the acceptance domain is used to obtain the preset constraint formula corresponding to the linear expression; wherein the preset constraint formula includes the branch power flow constraint formula, and the preset constraint formula at least includes: the power balance constraint of the operating base point, the maximum adjustment constraint of the generator set, and the preset constraint of the generator set; based on the priority target programming method and the preset constraint formula, the maximum value of the acceptance domain corresponding to the linear expression is determined; and then the value range constraint of the acceptance domain is obtained; The acceptance domain value range acquisition module includes a preset acquisition unit, Used to obtain the operating base point power balance constraints: ； Where, f is the f-th generator set, NG is the number of all generator sets; i is the i-th node; NI is the number of all nodes; represents the expected value of active power output of the fth generator set, and The value range of is from the minimum value of the net load downward acceptable domain to the maximum value of the net load upward acceptable domain; represents the expected value of the net load of the ith node, and The value range of is from the lower limit of active net load disturbance to the upper limit of active net load disturbance; Get the maximum adjustment constraint of the generator set: and ; in, represents the participation factor, and ; It represents the maximum upward adjustment of the load of the jth generator set in the tth time interval; It represents the maximum downward adjustment of the load of the jth generator set in the tth time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; represents the lower limit of the active net load disturbance on the ith node in the tth time interval; Get the preset constraints of the generator set: 、 、 ； in, Indicates the lower limit of active power generated by the f-th generator set; It represents the upper limit of active power generated by the f-th generator set; Indicates the output active power; represents the upward climbing rate of the fth generator set; represents the downward climbing rate of the fth generator set; represents the duration of the tth time interval; represents the upward reserve rate of the fth generator in the tth time interval; represents the downward reserve rate of the fth generator in the tth time interval; The module for obtaining the value range of the acceptance domain includes a branch flow constraint obtaining unit. Used to obtain branch power flow constraint formula: ,h=1,2,......,NL; ,h=1,2,......,NL; in, represents the net load disturbance of the ith node, and i=1,2,...,NI; represents the expected value of the net load of the ith node; h represents the hth branch line; NL represents the total number of branch lines; represents the power transfer distribution factor of the net load of the ith node to the hth branch line; represents the power transfer distribution factor of the net load of the f-th generator set to the h-th branch line; represents the upper limit of the branch power flow of the h-th branch line; The module for obtaining the value range of the acceptance domain includes a range determination unit, Used to obtain the acceptance domain value range constraint based on the maximum value of the acceptance domain corresponding to the linear expression: ， in, represents the maximum value of the acceptance domain, t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; represents the lower limit of the active net load disturbance on the ith node in the tth time interval; The minimum load moment sum calculation module is used to obtain the calculation function of the preset lower-level minimum load moment sum corresponding to the reference power grid and the acceptance domain value range constraint, and then determine the specific value of the minimum load moment sum; wherein the reference power grid planning includes branch lines and planned sections, and the parameters involved in the calculation function include branch lines, planned sections, and branch power flow upper limit; The minimum load moment sum calculation module includes a calculation function acquisition unit, The calculation function used to obtain the sum of the preset lower-level minimum load moments corresponding to the reference grid: , where s represents the sth planned road section, NS represents the number of all planned road sections; h represents the hth branch road; NL represents the total number of branch roads; represents the length of the sth planned road section, represents the upper limit of the branch power flow of the h-th branch line.

6. The reference power grid planning system based on load moment and acceptance domain according to claim 5, characterized in that: The linear expression acquisition module includes a function acquisition unit, Used to obtain the preset upper acceptance domain maximum calculation function: ]； Where Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; It represents the length of the upward acceptable domain of the net load of the i-th node in the t-th time interval; It represents the length of the downward acceptable domain of the net load of the i-th node in the t-th time interval; represents the upper limit of the active net load disturbance on the ith node in the tth time interval; It represents the lower limit of the active net load disturbance on the i-th node in the t-th time interval.

7. The reference power grid planning system based on load moment and acceptance domain according to claim 5 or 6, characterized in that: The linear expression acquisition module includes an expression acquisition unit, It is used to introduce two non-negative continuous variables and variable constraints on the basis of the preset upper acceptance domain maximum value calculation function, and to modify the preset upper acceptance domain maximum value calculation function to obtain a linear expression of the calculation function; Among them, two non-negative continuous variables are introduced as follows: and ; and , Z represents the maximum value of the upper acceptance domain; t represents the tth time interval, NT represents the number of all time intervals within the preset optimization time scale; i represents the i-th node; NI represents the number of all nodes; The constraints on the variables are: and .

8. A non-volatile computer storage medium, characterized in that: Computer instructions are stored thereon, and when the computer instructions are executed, the reference power grid planning method based on load moment and acceptance domain as described in any one of claims 1 to 4 is implemented.

Images