CN111106603A - Method and system for identifying limit interval of thermal stability power of power transmission section - Google Patents

Method and system for identifying limit interval of thermal stability power of power transmission section Download PDF

Info

Publication number
CN111106603A
CN111106603A CN201911366425.3A CN201911366425A CN111106603A CN 111106603 A CN111106603 A CN 111106603A CN 201911366425 A CN201911366425 A CN 201911366425A CN 111106603 A CN111106603 A CN 111106603A
Authority
CN
China
Prior art keywords
power
fault
limit
thermal stability
node
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911366425.3A
Other languages
Chinese (zh)
Other versions
CN111106603B (en
Inventor
鲍颜红
张金龙
徐泰山
薛峰
徐伟
任先成
周海锋
宋东阔
洪姗姗
戴玉臣
袁震
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nari Technology Co Ltd
Original Assignee
Nari Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nari Technology Co Ltd filed Critical Nari Technology Co Ltd
Priority to CN201911366425.3A priority Critical patent/CN111106603B/en
Publication of CN111106603A publication Critical patent/CN111106603A/en
Priority to PCT/CN2020/114927 priority patent/WO2021128959A1/en
Application granted granted Critical
Publication of CN111106603B publication Critical patent/CN111106603B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a method and a system for identifying a limit interval of thermal stability power of a power transmission section. Solving an optimization model with the minimum power of the transmission section under the fault as a target, and solving a lower limit of a thermal stability power limit interval based on all optimization results; constructing a sub-optimization model taking minimum active power difference injected by a node before/after a fault as a target, constructing a main optimization model taking maximum ground state transmission section power as the target and considering the Benders cut constraint, and iteratively solving the main and sub-optimization models to determine the upper limit of a thermal stability power limit interval; the traditional thermal stability power limit calculation problem is converted into a nonlinear optimization problem based on Benders decomposition, the identified limit interval is more accurate, and the method has guiding significance for dispatching personnel to fully master the safety and stability boundary of the power grid, guarantee the safety and stability operation of the power grid and fully utilize the power transmission capability of the power transmission section.

Description

Method and system for identifying limit interval of thermal stability power of power transmission section
Technical Field
The invention relates to a method and a system for identifying a limit interval of thermal stability power of a power transmission section, and belongs to the technical field of power system automation.
Background
The power transmission section is used as a power transmission channel and an electrical communication corridor, and when the transmission power is too large, potential safety and stability hazards exist. The power of a transmission section is a key index for dispatching operators to monitor the safe operation level of a power grid, and the thermal stability constraint is a main factor for limiting the transmission power of the transmission section. For the transmission capacity of the transmission section, when the section reaches the thermal stability power limit, the corresponding power flow adjustment modes are different, the corresponding maximum transmission capacities are inconsistent, and the thermal stability power limit of the section is an interval formed by the maximum value and the minimum value. The existing method for identifying the thermal stability power limit interval of the section (such as the patent 201910264792.6) has poor identification accuracy.
Disclosure of Invention
The invention provides a method and a system for identifying a thermal stability power limit interval of a power transmission section, which solve the problem of poor identification accuracy of the existing method.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows:
a method for identifying a limit interval of a thermal stability power of a power transmission section comprises the following steps,
acquiring a section-related thermal stability mode;
screening out a key thermal stability mode from the cross section related thermal stability modes;
obtaining a lower limit of a thermal stability power limit interval according to the key thermal stability mode and a pre-constructed first optimization model; the first optimization model takes the minimum power of a transmission section under a fault as a target;
responding to the total number of the faults of all the key thermal stability modes as 1, and obtaining the upper limit of the thermal stability power limit interval according to a second pre-constructed optimization model; the second optimization model takes the maximum power of the transmission section under the fault as a target;
responding to the fact that the total number of faults of all key thermal stability modes is larger than 1, and obtaining the upper limit of a thermal stability power limit interval according to a main optimization model and a sub-optimization model which are constructed in advance; the main optimization model aims at enabling ground state transmission section power to be maximum and includes Benders cut constraint, and the sub-optimization model aims at enabling active power difference injected into nodes before/after faults to be minimum.
The cross-section-dependent thermostabilization mode comprises: the main power flow transfer element is a thermal stability mode of a section component, the thermal stability mode comprises a fault of a power grid and a corresponding main power flow transfer element, and the main power flow transfer element with the fault is a thermal stability investigation element with a power flow transfer ratio or a load rate variation larger than a threshold value under the fault.
The process of screening out the key thermal stability mode from the section-related thermal stability modes comprises the following steps,
solving a pre-constructed third optimization model aiming at each section-related thermal stability mode, and responding to the fact that the optimization result is not smaller than the current limit value of a section component, wherein the section-related thermal stability mode is a key thermal stability mode; wherein the third optimization model targets the fault lower section component current maximum.
The third optimization model is that,
objective function f 1:
Figure BDA0002338541890000021
Figure BDA0002338541890000022
constraint conditions are as follows:
Figure BDA0002338541890000031
Figure BDA0002338541890000032
Figure BDA0002338541890000033
Figure BDA0002338541890000034
Figure BDA0002338541890000035
Figure BDA0002338541890000036
Figure BDA0002338541890000037
wherein the content of the first and second substances,
Figure BDA0002338541890000038
is the branch current amplitude between the nodes I and j under the fault x, Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure BDA0002338541890000039
respectively the active power of a generator, the active power of a load, the reactive power of the generator and the reactive power of the load of a node i under the fault x,
Figure BDA00023385418900000310
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (x)、
Figure BDA00023385418900000311
The voltage at node i and the voltage at node j at fault x,
Figure BDA00023385418900000312
respectively the lower and upper voltage limits of node i at fault x,
Figure BDA00023385418900000313
respectively is the lower limit of the load active power, the upper limit of the load active power, the lower limit of the power factor angle and the upper limit of the power factor angle of the node i,
Figure BDA00023385418900000314
is the power factor angle of node i at fault x, N is the number of nodes except node i,
Figure BDA00023385418900000315
amplitude and phase angle of admittance between nodes i and j under fault x, respectively
Figure BDA00023385418900000316
Figure BDA00023385418900000317
The phase angles of the nodes i and j under the fault x are respectively; z is all main power flow transfer elements corresponding to the fault x, and y is a fault lower section component.
The obtaining of the lower limit of the thermal stability power limit interval according to the key thermal stability mode and the pre-constructed first optimization model comprises:
and traversing all key thermal stability modes, solving a pre-constructed first optimization model, solving the power of the power transmission section before the fault based on an optimization result, and taking the minimum value of the power transmission section before the fault as the lower limit of the limit interval of the thermal stability power.
The first optimization model is a model of,
objective function f 2:
Figure BDA0002338541890000041
Figure BDA0002338541890000042
constraint conditions are as follows:
Figure BDA0002338541890000043
Figure BDA0002338541890000044
Figure BDA0002338541890000045
Figure BDA0002338541890000046
Figure BDA0002338541890000047
Figure BDA0002338541890000048
Figure BDA0002338541890000049
Figure BDA00023385418900000410
wherein the content of the first and second substances,
Figure BDA00023385418900000411
for the active power of the branch between nodes i and j under fault x, SLThe cross section is formed into a branch set,
Figure BDA0002338541890000051
is the branch current amplitude between the nodes I and j under the fault x, Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure BDA0002338541890000052
respectively the active power of a generator, the active power of a load, the reactive power of the generator and the reactive power of the load of a node i under the fault x,
Figure BDA0002338541890000053
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (x)、
Figure BDA0002338541890000054
The voltage at node i and the voltage at node j at fault x,
Figure BDA0002338541890000055
respectively the lower and upper voltage limits of node i at fault x,
Figure BDA0002338541890000056
respectively is the lower limit of the load active power, the upper limit of the load active power, the lower limit of the power factor angle and the upper limit of the power factor angle of the node i,
Figure BDA0002338541890000057
is the power factor angle of node i at fault x, N is the number of nodes except node i,
Figure BDA0002338541890000058
is the amplitude, parameter, of admittance between nodes i and j under fault x
Figure BDA0002338541890000059
Figure BDA00023385418900000510
The phase angles of the nodes i and j under the fault x are respectively, and epsilon is a threshold value; and Z is all main power flow transfer elements corresponding to the fault x.
The step of obtaining the thermal stability power limit interval upper limit according to a second pre-constructed optimization model in response to the total number of faults of all key thermal stability modes being 1 includes:
and responding to the total number of the faults of all the key thermal stability modes as 1, solving a pre-constructed second optimization model, solving the power of the power transmission section before the fault based on an optimization result, and taking the power of the power transmission section before the fault as the upper limit of the limit interval of the thermal stability power.
The second optimization model is a model of,
objective function f 3:
Figure BDA00023385418900000511
Figure BDA00023385418900000512
constraint conditions are as follows:
Figure BDA0002338541890000061
Figure BDA0002338541890000062
Figure BDA0002338541890000063
Figure BDA0002338541890000064
Figure BDA0002338541890000065
Figure BDA0002338541890000066
Figure BDA0002338541890000067
wherein the content of the first and second substances,
Figure BDA0002338541890000068
for the active power of the branch between nodes i and j under fault x, SLThe cross section is formed into a branch set,
Figure BDA0002338541890000069
is the branch current amplitude between the nodes I and j under the fault x, Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure BDA00023385418900000610
respectively the active power of a generator, the active power of a load, the reactive power of the generator and the reactive power of the load of a node i under the fault x,
Figure BDA00023385418900000611
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (x)、
Figure BDA00023385418900000612
The voltage at node i and the voltage at node j at fault x,
Figure BDA00023385418900000613
respectively the lower and upper voltage limits of node i at fault x,
Figure BDA00023385418900000614
respectively is the lower limit of the load active power, the upper limit of the load active power, the lower limit of the power factor angle and the upper limit of the power factor angle of the node i,
Figure BDA00023385418900000615
is the power factor angle of node i at fault x, N is the number of nodes except node i,
Figure BDA00023385418900000616
is the amplitude, parameter, of admittance between nodes i and j under fault x
Figure BDA00023385418900000617
Figure BDA00023385418900000618
The phase angles of the nodes i and j under the fault x are respectively; and Z is all main power flow transfer elements corresponding to the fault x.
The sub-optimization model is that,
objective function fx.sub：
Figure BDA0002338541890000071
Constraint conditions are as follows:
Figure BDA0002338541890000072
Figure BDA0002338541890000073
Figure BDA0002338541890000074
Figure BDA0002338541890000075
Figure BDA0002338541890000076
Figure BDA0002338541890000077
Figure BDA0002338541890000078
Figure BDA0002338541890000079
wherein N isgIn order to adjust the number of generators and load nodes,
Figure BDA00023385418900000710
in order to introduce the virtual variables,
Figure BDA00023385418900000711
in order to be a function of the involved variables,
Figure BDA00023385418900000712
is the branch current amplitude between the nodes I and j under the fault x, Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure BDA00023385418900000713
respectively the active power of a generator, the active power of a load, the reactive power of the generator and the reactive power of the load of a node i under the fault x,
Figure BDA00023385418900000714
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (x)、
Figure BDA00023385418900000715
The voltage at node i and the voltage at node j at fault x,
Figure BDA00023385418900000716
respectively the lower and upper voltage limits of node i at fault x,
Figure BDA0002338541890000081
respectively is the lower limit of the load active power, the upper limit of the load active power, the lower limit of the power factor angle and the upper limit of the power factor angle of the node i,
Figure BDA0002338541890000082
is the power factor angle of node i at fault x, N is the number of nodes except node i,
Figure BDA0002338541890000083
is the amplitude, parameter, of admittance between nodes i and j under fault x
Figure BDA0002338541890000084
Figure BDA0002338541890000085
The phase angles of the nodes i and j under the fault x are respectively, and Z is all the main power flow transfer elements corresponding to the fault x.
The main optimization model is that,
objective function f 4:
Figure BDA0002338541890000086
Figure BDA0002338541890000087
constraint conditions are as follows:
Figure BDA0002338541890000088
Figure BDA0002338541890000089
Figure BDA00023385418900000810
Figure BDA00023385418900000811
Figure BDA00023385418900000812
Figure BDA00023385418900000813
Figure BDA00023385418900000814
responsive to the primary optimization model for the first timeIn the calculation, the calculation is carried out,
Figure BDA00023385418900000815
x∈SX；
Figure BDA00023385418900000816
in response to the primary optimization model being non-first calculated,
Figure BDA0002338541890000091
wherein the content of the first and second substances,
Figure BDA0002338541890000092
is the active power of the branch between the nodes i and j under the ground state, SLThe cross section is formed into a branch set,
Figure BDA0002338541890000093
the amplitude of the current of the branch between the nodes I and j under the ground state is Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure BDA0002338541890000094
respectively the active power, the load active power, the reactive power and the load reactive power of the generator at the node i under the ground state,
Figure BDA0002338541890000095
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (0)、
Figure BDA0002338541890000096
The voltage at node i and the voltage at node j in the ground state,
Figure BDA0002338541890000097
respectively a lower voltage limit and an upper voltage limit of the node i in the ground state,
Figure BDA0002338541890000098
respectively is the lower limit of the load active power, the upper limit of the load active power of the node i,A lower power factor angle limit and an upper power factor angle limit,
Figure BDA0002338541890000099
is the power factor angle of node i in the ground state, N is the number of nodes excluding node i,
Figure BDA00023385418900000910
amplitude, parameter, of admittance between nodes i, j under ground state
Figure BDA00023385418900000911
Figure BDA00023385418900000912
Phase angles, S, of nodes i, j, respectively, in the ground stateXIs a set of critical faults, the critical faults are faults in a critical thermal stability mode,
Figure BDA00023385418900000913
is a set of critical elements under the fault x, the critical elements are main power flow transfer elements corresponding to the fault in the critical thermal stability mode,
Figure BDA00023385418900000914
for the power flow transfer ratio of critical component y at fault x,
Figure BDA00023385418900000915
the power is estimated for the critical safety of the critical component y under fault x,
Figure BDA00023385418900000916
respectively active power of key element y under the ground state, active power of fault element corresponding to fault x and NgTo adjust the number of generator and load nodes, mu is the acceleration factor,
Figure BDA00023385418900000917
multiplier, f, corresponding to the containment constraint of the subproblem under fault xx.subIn order to sub-optimize the model objective function,
Figure BDA00023385418900000918
is the solution of the main problem in the previous round.
A power transmission section thermal stability power limit interval identification system comprises,
an acquisition module: the method comprises the steps of obtaining a section-related thermal stability mode;
a screening module: the method is used for screening out a key thermal stability mode from the related thermal stability modes of the section;
a lower limit calculation module: the lower limit of the thermal stability power limit interval is obtained according to the key thermal stability mode and a pre-constructed first optimization model; the first optimization model takes the minimum power of a transmission section under a fault as a target;
the single fault upper limit solving module: the method comprises the steps that the total number of faults responding to all key thermal stability modes is 1, and the upper limit of a thermal stability power limit interval is obtained according to a second pre-constructed optimization model; the second optimization model takes the maximum power of the transmission section under the fault as a target;
a multi-fault upper limit solving module: the method comprises the steps that the total number of faults of all key thermal stability modes is larger than 1, and the upper limit of a thermal stability power limit interval is obtained according to a pre-constructed main optimization model and a pre-constructed sub-optimization model; the main optimization model aims at enabling ground state transmission section power to be maximum and includes Benders cut constraint, and the sub-optimization model aims at enabling active power difference injected into nodes before/after faults to be minimum.
A computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform a power transmission profile thermally stable power limit interval identification method.
A computing device comprising one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs comprising instructions for performing a power transmission profile thermally stable power limit interval identification method.
The invention achieves the following beneficial effects: the method solves an optimization model with the minimum power of the transmission section under the fault as a target, and solves the lower limit of the limit interval of the thermal stability power based on all optimization results; constructing a sub-optimization model taking minimum active power difference injected by a node before/after a fault as a target, constructing a main optimization model taking maximum ground state transmission section power as the target and considering the Benders cut constraint, and iteratively solving the main and sub-optimization models to determine the upper limit of a thermal stability power limit interval; the traditional thermal stability power limit calculation problem is converted into a nonlinear optimization problem based on Benders decomposition, the identified limit interval is more accurate, and the method has guiding significance for dispatching personnel to fully master the safety and stability boundary of the power grid, guarantee the safety and stability operation of the power grid and fully utilize the power transmission capability of the power transmission section.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
As shown in fig. 1, a method for identifying a thermally stable power limit interval of a power transmission section includes the following steps:
step 1, a section-related thermal stability mode is obtained, wherein the section-related thermal stability mode comprises faults of a power grid and corresponding main power flow transfer elements.
Identifying a main power flow transfer element based on a ground state and a power flow after a fault, specifically, taking a thermal stability investigation element with a power flow transfer ratio or a load rate variation larger than a threshold value under the fault as the main power flow transfer element of the fault; wherein, the power flow transfer ratio formula of the element for thermal stability inspection under fault is as follows,
Figure BDA0002338541890000111
wherein the content of the first and second substances,
Figure BDA0002338541890000112
other than the initial modeThe active power of the element y' and the sum of the active power of the element corresponding to the fault x are examined in the lower thermal stability,
Figure BDA0002338541890000113
the active power of element y' is examined for thermal stability after the initial mode failure x.
Each thermal stability mode only comprises one main power flow transfer element, and if the same fault corresponds to a plurality of main power flow transfer elements, a plurality of thermal stability modes are formed.
And 2, screening out a key thermal stability mode from the cross section related thermal stability modes.
Aiming at each section-related thermal stability mode, solving a pre-constructed third optimization model based on an interior point method, and responding to the fact that the optimization result is not less than the current limit value of a section component, wherein the section-related thermal stability mode is a key thermal stability mode; in response to the optimization result being less than the current limit of the section component element, the section-related thermostabilization mode is set to be invalid; the third optimization model is used for forming the maximum target of the element current by using the fault lower section; and if the number of the screened key thermal stability modes is 0, the power transmission section is not limited.
The third optimization model is:
objective function f 1:
Figure BDA0002338541890000121
Figure BDA0002338541890000122
constraint conditions are as follows:
Figure BDA0002338541890000123
Figure BDA0002338541890000124
Figure BDA0002338541890000125
Figure BDA0002338541890000126
Figure BDA0002338541890000127
Figure BDA0002338541890000128
Figure BDA0002338541890000131
wherein the content of the first and second substances,
Figure BDA0002338541890000132
is the branch current amplitude between the nodes I and j under the fault x, Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure BDA0002338541890000133
respectively the active power of a generator, the active power of a load, the reactive power of the generator and the reactive power of the load of a node i under the fault x,
Figure BDA0002338541890000134
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (x)、
Figure BDA0002338541890000135
The voltage at node i and the voltage at node j at fault x,
Figure BDA0002338541890000136
respectively the lower and upper voltage limits of node i at fault x,
Figure BDA0002338541890000137
respectively is the lower limit of the load active power, the upper limit of the load active power of the node i,A lower power factor angle limit and an upper power factor angle limit,
Figure BDA0002338541890000138
is the power factor angle of node i at fault x, N is the number of nodes except node i,
Figure BDA0002338541890000139
amplitude and phase angle of admittance between nodes i and j under fault x, respectively
Figure BDA00023385418900001310
The phase angles of the nodes i and j under the fault x are respectively; z is all main power flow transfer elements corresponding to the fault x, and y is a fault lower section component.
And 3, constructing a key fault set by taking a fault union set of each key thermal stability mode, and constructing a key element set of each fault by taking a main power flow transfer element union set corresponding to each key fault.
Step 4, obtaining a lower limit of a thermal stability power limit interval according to the key thermal stability mode and a pre-constructed first optimization model; wherein the first optimization model aims at the minimum power of the transmission section under the fault.
The method specifically comprises the following steps: traversing all key thermal stability modes, solving a pre-constructed first optimization model based on an interior point method, solving power of a power transmission section before fault based on an optimization result, and taking the minimum value of the power transmission section before fault as the lower limit of a thermal stability power limit interval.
The first optimization model is:
objective function f 2:
Figure BDA0002338541890000141
Figure BDA0002338541890000142
constraint conditions are as follows:
formulae (a) to (g);
Figure BDA0002338541890000143
wherein the content of the first and second substances,
Figure BDA0002338541890000144
for the active power of the branch between nodes i and j under fault x, SLA set of branches is formed for the section, and epsilon is a threshold value.
Step 5, responding to the fact that the total number of faults of all key thermal stability modes is 1, and obtaining the upper limit of a thermal stability power limit interval according to a second optimization model which is constructed in advance; wherein the second optimization model targets the transmission section power under fault as maximum.
The method specifically comprises the following steps: and responding to the total number of the faults of all the key thermal stability modes as 1, solving a pre-constructed second optimization model, solving the power of the power transmission section before the fault based on an optimization result, and taking the power of the power transmission section before the fault as the upper limit of the limit interval of the thermal stability power.
The second optimization model is:
objective function f 3:
Figure BDA0002338541890000145
Figure BDA0002338541890000146
constraint conditions are as follows: formulae (a) to (f);
Figure BDA0002338541890000147
step 6, responding to the fact that the total number of faults of all key thermal stability modes is larger than 1, solving and continuously updating the involvement variables (Bendserkin constraint, active power injection of nodes before/after faults) of the main sub-problem according to a pre-constructed main optimization model and a sub-optimization model based on an interior point method until the objective function value of the sub-optimization model is smaller than a threshold value and the main optimization model meets a convergence condition, and taking the objective function value of the main optimization model as the upper limit of a thermal stability power limit interval; the main optimization model aims at enabling ground state transmission section power to be maximum and includes Benders cut constraint, and the sub-optimization model aims at enabling active power difference injected into nodes before/after faults to be minimum.
The sub-optimization model is:
objective function fx.sub：
Figure BDA0002338541890000151
Constraint conditions are as follows:
formulae (a) to (f);
Figure BDA0002338541890000152
Figure BDA0002338541890000153
wherein N isgIn order to adjust the number of generators and load nodes,
Figure BDA0002338541890000154
the introduced virtual variable has the function of temporarily relieving the constraint of the main and sub problems when the constraint of the involvement of the main and sub problems is acted, ensuring the solution of the sub problems and simultaneously having the correlation function of the connection and the constraint of the main problems,
Figure BDA0002338541890000155
to tie up variables, from the main problem, it remains unchanged during the sub-problem processing.
The main optimization model is as follows:
objective function f 4:
Figure BDA0002338541890000156
Figure BDA0002338541890000161
constraint conditions are as follows:
Figure BDA0002338541890000162
Figure BDA0002338541890000163
Figure BDA0002338541890000164
Figure BDA0002338541890000165
Figure BDA0002338541890000166
Figure BDA0002338541890000167
Figure BDA0002338541890000168
responding to the main optimization model, calculating for the first time, and estimating static overload constraint based on the power flow sensitivity:
Figure BDA0002338541890000169
x∈SX；
Figure BDA00023385418900001610
responding to the non-first calculation of the main optimization model, solving the Bendbiscut constraint of the subproblem conversion with the result larger than the specified threshold value:
Figure BDA00023385418900001611
wherein the content of the first and second substances,
Figure BDA00023385418900001612
is the active power of the branch between the nodes i and j under the ground state,
Figure BDA00023385418900001613
the amplitude of the branch current between the nodes i and j under the ground state,
Figure BDA00023385418900001614
the active power, the load active power, the reactive power and the load reactive power of the generator of a node i under the ground state respectively, Vi (0)、
Figure BDA00023385418900001615
The voltage at node i and the voltage at node j in the ground state,
Figure BDA00023385418900001616
respectively a lower voltage limit and an upper voltage limit of the node i in the ground state,
Figure BDA00023385418900001617
is the power factor angle of node i in the ground state, N is the number of nodes excluding node i,
Figure BDA0002338541890000171
amplitude, parameter, of admittance between nodes i, j under ground state
Figure BDA0002338541890000172
Figure BDA0002338541890000173
Phase angles, S, of nodes i, j, respectively, in the ground stateXIs a set of critical faults, the critical faults are faults in a critical thermal stability mode,
Figure BDA0002338541890000174
is a set of critical elements under the fault x, the critical elements are main power flow transfer elements corresponding to the fault in the critical thermal stability mode,
Figure BDA0002338541890000175
for fault x lower gateThe power flow transfer ratio of the key element y,
Figure BDA0002338541890000176
the power is estimated for the critical safety of the critical component y under fault x,
Figure BDA0002338541890000177
respectively the active power of the key element y in the ground state, the active power of the fault element corresponding to the fault x and the sum of the active powers,
Figure BDA0002338541890000178
mu is an acceleration factor, N is a multiplier corresponding to the containment constraint of the subproblem under the fault xgIn order to adjust the number of generator and load nodes,
Figure BDA0002338541890000179
in order to solve the main problem in the previous round,
Figure BDA00023385418900001710
respectively the active power and the current amplitude, I, of the key element y after the failure x in the initial modemax.yIs the upper current amplitude limit for the critical component y.
fx.sub、
Figure BDA00023385418900001711
From the sub-problem under the fault x,
Figure BDA00023385418900001712
the solution from the previous round of main problem is kept unchanged in the current round of main problem processing;
Figure BDA00023385418900001713
the multiplier corresponding to the tie constraint of the subproblem under the fault x corresponds to the inequality constraint condition in the subproblem under the fault x
Figure BDA00023385418900001714
When the representation subproblem obtains the optimal solution, the objective function value corresponds to
Figure BDA00023385418900001715
(ii) a sensitivity of change; the so-called Benders cut constraint is that
Figure BDA00023385418900001716
And the operation mode of the power grid is finely adjusted on the basis, so that the subproblems are feasible.
The method solves the optimization model with the minimum power of the transmission section under the fault as the target, and solves the lower limit of the limit interval of the thermal stability power based on all the optimization results; constructing a sub-optimization model taking minimum active power difference injected by a node before/after a fault as a target, constructing a main optimization model taking maximum ground state transmission section power as the target and considering the Benders cut constraint, and iteratively solving the main and sub-optimization models to determine the upper limit of a thermal stability power limit interval; the traditional thermal stability power limit calculation problem is converted into a nonlinear optimization problem based on Benders decomposition, the identified limit interval is more accurate, and the method has guiding significance for dispatching personnel to fully master the safety and stability boundary of the power grid, guarantee the safety and stability operation of the power grid and fully utilize the power transmission capability of the power transmission section.
A power transmission section thermal stability power limit interval identification system comprises,
an acquisition module: the method comprises the steps of obtaining a section-related thermal stability mode;
a screening module: the method is used for screening out a key thermal stability mode from the related thermal stability modes of the section;
a lower limit calculation module: the lower limit of the thermal stability power limit interval is obtained according to the key thermal stability mode and a pre-constructed first optimization model; the first optimization model takes the minimum power of a transmission section under a fault as a target;
the single fault upper limit solving module: the method comprises the steps that the total number of faults responding to all key thermal stability modes is 1, and the upper limit of a thermal stability power limit interval is obtained according to a second pre-constructed optimization model; the second optimization model takes the maximum power of the transmission section under the fault as a target;
a multi-fault upper limit solving module: the method comprises the steps that the total number of faults of all key thermal stability modes is larger than 1, and the upper limit of a thermal stability power limit interval is obtained according to a pre-constructed main optimization model and a pre-constructed sub-optimization model; the main optimization model aims at enabling ground state transmission section power to be maximum and includes Benders cut constraint, and the sub-optimization model aims at enabling active power difference injected into nodes before/after faults to be minimum.
A computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform a power transmission profile thermally stable power limit interval identification method.
A computing device comprising one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs comprising instructions for performing a power transmission profile thermally stable power limit interval identification method.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The present invention is not limited to the above embodiments, and any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention are included in the scope of the claims of the present invention which are filed as the application.

Claims (13)

1. A method for identifying a limit interval of thermal stability power of a power transmission section is characterized by comprising the following steps: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
acquiring a section-related thermal stability mode;
screening out a key thermal stability mode from the cross section related thermal stability modes;
obtaining a lower limit of a thermal stability power limit interval according to the key thermal stability mode and a pre-constructed first optimization model; the first optimization model takes the minimum power of a transmission section under a fault as a target;
responding to the total number of the faults of all the key thermal stability modes as 1, and obtaining the upper limit of the thermal stability power limit interval according to a second pre-constructed optimization model; the second optimization model takes the maximum power of the transmission section under the fault as a target;
responding to the fact that the total number of faults of all key thermal stability modes is larger than 1, and obtaining the upper limit of a thermal stability power limit interval according to a main optimization model and a sub-optimization model which are constructed in advance; the main optimization model aims at enabling ground state transmission section power to be maximum and includes Benders cut constraint, and the sub-optimization model aims at enabling active power difference injected into nodes before/after faults to be minimum.
2. The method for identifying a thermally stable power limit interval of a power transmission section according to claim 1, wherein: the cross-section-dependent thermostabilization mode comprises: the main power flow transfer element is a thermal stability mode of a section component, the thermal stability mode comprises a fault of a power grid and a corresponding main power flow transfer element, and the main power flow transfer element with the fault is a thermal stability investigation element with a power flow transfer ratio or a load rate variation larger than a threshold value under the fault.
3. A method for identifying a thermally stable power limit interval of a power transmission section according to claim 1 or 2, characterized in that: the process of screening out the key thermal stability mode from the section-related thermal stability modes comprises the following steps,
solving a pre-constructed third optimization model aiming at each section-related thermal stability mode, and responding to the fact that the optimization result is not smaller than the current limit value of a section component, wherein the section-related thermal stability mode is a key thermal stability mode; wherein the third optimization model targets the fault lower section component current maximum.
4. A method for identifying a thermally stable power limit interval of a power transmission section according to claim 3, characterized in that: the third optimization model is that,
objective function f 1:
Figure FDA0002338541880000021
Figure FDA0002338541880000022
constraint conditions are as follows:
Figure FDA0002338541880000023
Figure FDA0002338541880000024
Figure FDA0002338541880000025
Figure FDA00023385418800000211
Figure FDA0002338541880000026
Figure FDA0002338541880000027
Figure FDA0002338541880000028
wherein the content of the first and second substances,
Figure FDA0002338541880000029
is the branch current amplitude between the nodes I and j under the fault x, Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure FDA00023385418800000210
respectively the active power of a generator, the active power of a load, the reactive power of the generator and the reactive power of the load of a node i under the fault x,
Figure FDA0002338541880000031
the lower limit of the active power, the upper limit of the active power and the lower limit of the reactive power of the generator which are respectively a node iUpper limit of reactive power, Vi (x)、
Figure FDA0002338541880000032
The voltage at node i and the voltage at node j at fault x,
Figure FDA0002338541880000033
respectively the lower and upper voltage limits of node i at fault x,
Figure FDA0002338541880000034
respectively is the lower limit of the load active power, the upper limit of the load active power, the lower limit of the power factor angle and the upper limit of the power factor angle of the node i,
Figure FDA0002338541880000035
is the power factor angle of node i at fault x, N is the number of nodes except node i,
Figure FDA0002338541880000036
amplitude and phase angle of admittance between nodes i and j under fault x, respectively
Figure FDA0002338541880000037
The phase angles of the nodes i and j under the fault x are respectively; z is all main power flow transfer elements corresponding to the fault x, and y is a fault lower section component.
5. The method for identifying a thermally stable power limit interval of a power transmission section according to claim 1, wherein: the obtaining of the lower limit of the thermal stability power limit interval according to the key thermal stability mode and the pre-constructed first optimization model comprises:
and traversing all key thermal stability modes, solving a pre-constructed first optimization model, solving the power of the power transmission section before the fault based on an optimization result, and taking the minimum value of the power transmission section before the fault as the lower limit of the limit interval of the thermal stability power.
6. A method for identifying a thermally stable power limit interval of a power transmission section according to claim 1 or 5, characterized in that: the first optimization model is a model of,
objective function f 2:
Figure FDA0002338541880000038
Figure FDA0002338541880000039
constraint conditions are as follows:
Figure FDA00023385418800000417
Figure FDA0002338541880000041
Figure FDA00023385418800000418
Figure FDA0002338541880000042
Figure FDA0002338541880000043
Figure FDA0002338541880000044
Figure FDA0002338541880000045
Figure FDA0002338541880000046
wherein the content of the first and second substances,
Figure FDA0002338541880000047
for the active power of the branch between nodes i and j under fault x, SLThe cross section is formed into a branch set,
Figure FDA0002338541880000048
is the branch current amplitude between the nodes I and j under the fault x, Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure FDA0002338541880000049
respectively the active power of a generator, the active power of a load, the reactive power of the generator and the reactive power of the load of a node i under the fault x,
Figure FDA00023385418800000410
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (x)、
Figure FDA00023385418800000411
The voltage at node i and the voltage at node j at fault x,
Figure FDA00023385418800000412
respectively the lower and upper voltage limits of node i at fault x,
Figure FDA00023385418800000413
respectively is the lower limit of the load active power, the upper limit of the load active power, the lower limit of the power factor angle and the upper limit of the power factor angle of the node i,
Figure FDA00023385418800000414
is the power factor angle of node i at fault x, N is the number of nodes except node i,
Figure FDA00023385418800000415
is the amplitude, parameter, of admittance between nodes i and j under fault x
Figure FDA00023385418800000419
Figure FDA00023385418800000416
The phase angles of the nodes i and j under the fault x are respectively, and epsilon is a threshold value; and Z is all main power flow transfer elements corresponding to the fault x.
7. The method for identifying a thermally stable power limit interval of a power transmission section according to claim 1, wherein: the step of obtaining the thermal stability power limit interval upper limit according to a second pre-constructed optimization model in response to the total number of faults of all key thermal stability modes being 1 includes:
and responding to the total number of the faults of all the key thermal stability modes as 1, solving a pre-constructed second optimization model, solving the power of the power transmission section before the fault based on an optimization result, and taking the power of the power transmission section before the fault as the upper limit of the limit interval of the thermal stability power.
8. A method for identifying a thermally stable power limit interval of a power transmission section according to claim 1 or 7, characterized in that: the second optimization model is a model of,
objective function f 3:
Figure FDA0002338541880000051
Figure FDA0002338541880000052
constraint conditions are as follows:
Figure FDA0002338541880000053
Figure FDA0002338541880000054
Figure FDA0002338541880000055
Figure FDA0002338541880000056
Figure FDA0002338541880000057
Figure FDA0002338541880000058
Figure FDA0002338541880000059
wherein the content of the first and second substances,
Figure FDA0002338541880000061
for the active power of the branch between nodes i and j under fault x, SLThe cross section is formed into a branch set,
Figure FDA0002338541880000062
is the branch current amplitude between the nodes I and j under the fault x, Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure FDA0002338541880000063
respectively the active power of a generator, the active power of a load, the reactive power of the generator and the reactive power of the load of a node i under the fault x,
Figure FDA0002338541880000064
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (x)、
Figure FDA0002338541880000065
The voltage at node i and the voltage at node j at fault x,
Figure FDA0002338541880000066
respectively the lower and upper voltage limits of node i at fault x,
Figure FDA0002338541880000067
respectively is the lower limit of the load active power, the upper limit of the load active power, the lower limit of the power factor angle and the upper limit of the power factor angle of the node i,
Figure FDA0002338541880000068
is the power factor angle of node i at fault x, N is the number of nodes except node i,
Figure FDA0002338541880000069
is the amplitude, parameter, of admittance between nodes i and j under fault x
Figure FDA00023385418800000610
Figure FDA00023385418800000611
The phase angles of the nodes i and j under the fault x are respectively; and Z is all main power flow transfer elements corresponding to the fault x.
9. The method for identifying a thermally stable power limit interval of a power transmission section according to claim 1, wherein: the sub-optimization model is that,
objective function fx.sub：
Figure FDA00023385418800000612
Constraint conditions are as follows:
Figure FDA00023385418800000613
Figure FDA00023385418800000614
Figure FDA00023385418800000615
Figure FDA0002338541880000071
Figure FDA0002338541880000072
Figure FDA0002338541880000073
Figure FDA0002338541880000074
Figure FDA0002338541880000075
wherein N isgIn order to adjust the number of generators and load nodes,
Figure FDA0002338541880000076
in order to introduce the virtual variables,
Figure FDA0002338541880000077
in order to be a function of the involved variables,
Figure FDA0002338541880000078
is the branch current amplitude between the nodes I and j under the fault x, Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure FDA0002338541880000079
respectively the active power of a generator, the active power of a load, the reactive power of the generator and the reactive power of the load of a node i under the fault x,
Figure FDA00023385418800000710
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (x)、
Figure FDA00023385418800000711
The voltage at node i and the voltage at node j at fault x,
Figure FDA00023385418800000712
respectively the lower and upper voltage limits of node i at fault x,
Figure FDA00023385418800000713
respectively is the lower limit of the load active power, the upper limit of the load active power, the lower limit of the power factor angle and the upper limit of the power factor angle of the node i,
Figure FDA00023385418800000714
is the power factor angle of node i at fault x, N is the number of nodes except node i,
Figure FDA00023385418800000715
is the amplitude, parameter, of admittance between nodes i and j under fault x
Figure FDA00023385418800000716
Figure FDA00023385418800000717
The phase angles of the nodes i and j under the fault x are respectively, and Z is all the main power flow transfer elements corresponding to the fault x.
10. The method for identifying a thermally stable power limit interval of a power transmission section according to claim 1, wherein: the main optimization model is that,
objective function f 4:
Figure FDA00023385418800000718
Figure FDA0002338541880000081
constraint conditions are as follows:
Figure FDA0002338541880000082
Figure FDA0002338541880000083
Figure FDA0002338541880000084
Figure FDA0002338541880000085
Figure FDA0002338541880000086
Figure FDA0002338541880000087
Figure FDA0002338541880000088
in response to the primary optimization model being calculated for the first time,
Figure FDA0002338541880000089
x∈SX；
Figure FDA00023385418800000810
in response to the primary optimization model being non-first calculated,
Figure FDA00023385418800000811
wherein the content of the first and second substances,
Figure FDA00023385418800000812
is the active power of the branch between the nodes i and j under the ground state, SLThe cross section is formed into a branch set,
Figure FDA00023385418800000813
the amplitude of the current of the branch between the nodes I and j under the ground state is Imax.ijIs the upper limit of the branch current amplitude between the nodes i and j,
Figure FDA00023385418800000814
respectively the active power, the load active power, the reactive power and the load reactive power of the generator at the node i under the ground state,
Figure FDA00023385418800000815
the lower active limit, the upper active limit, the lower reactive limit and the upper reactive limit of the generator of the node i, Vi (0)、
Figure FDA00023385418800000816
The voltage at node i and the voltage at node j in the ground state,
Figure FDA00023385418800000817
respectively a lower voltage limit and an upper voltage limit of the node i in the ground state,
Figure FDA0002338541880000091
respectively is the lower limit of the load active power, the upper limit of the load active power, the lower limit of the power factor angle and the upper limit of the power factor angle of the node i,
Figure FDA0002338541880000092
is the power factor angle of node i in the ground state, N is the number of nodes excluding node i,
Figure FDA0002338541880000093
amplitude, parameter, of admittance between nodes i, j under ground state
Figure FDA0002338541880000094
Figure FDA0002338541880000095
Phase angles, S, of nodes i, j, respectively, in the ground stateXIs a set of critical faults, the critical faults are faults in a critical thermal stability mode,
Figure FDA0002338541880000096
is a set of critical elements under the fault x, the critical elements are main power flow transfer elements corresponding to the fault in the critical thermal stability mode,
Figure FDA0002338541880000097
for the power flow transfer ratio of critical component y at fault x,
Figure FDA0002338541880000098
the power is estimated for the critical safety of the critical component y under fault x,
Figure FDA0002338541880000099
respectively active power of key element y under the ground state, active power of fault element corresponding to fault x and NgTo adjust the number of generator and load nodes, mu is the acceleration factor,
Figure FDA00023385418800000910
multiplier, f, corresponding to the containment constraint of the subproblem under fault xx.subIn order to sub-optimize the model objective function,
Figure FDA00023385418800000911
is the solution of the main problem in the previous round.
11. A power transmission section thermal stability power limit interval identification system is characterized in that: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
an acquisition module: the method comprises the steps of obtaining a section-related thermal stability mode;
a screening module: the method is used for screening out a key thermal stability mode from the related thermal stability modes of the section;
a lower limit calculation module: the lower limit of the thermal stability power limit interval is obtained according to the key thermal stability mode and a pre-constructed first optimization model; the first optimization model takes the minimum power of a transmission section under a fault as a target;
the single fault upper limit solving module: the method comprises the steps that the total number of faults responding to all key thermal stability modes is 1, and the upper limit of a thermal stability power limit interval is obtained according to a second pre-constructed optimization model; the second optimization model takes the maximum power of the transmission section under the fault as a target;
a multi-fault upper limit solving module: the method comprises the steps that the total number of faults of all key thermal stability modes is larger than 1, and the upper limit of a thermal stability power limit interval is obtained according to a pre-constructed main optimization model and a pre-constructed sub-optimization model; the main optimization model aims at enabling ground state transmission section power to be maximum and includes Benders cut constraint, and the sub-optimization model aims at enabling active power difference injected into nodes before/after faults to be minimum.
12. A computer readable storage medium storing one or more programs, characterized in that: the one or more programs include instructions that, when executed by a computing device, cause the computing device to perform any of the methods of claims 1-10.
13. A computing device, characterized by: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for performing any of the methods of claims 1-10.
CN201911366425.3A 2019-12-26 2019-12-26 Method and system for identifying limit interval of thermal stability power of power transmission section Active CN111106603B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201911366425.3A CN111106603B (en) 2019-12-26 2019-12-26 Method and system for identifying limit interval of thermal stability power of power transmission section
PCT/CN2020/114927 WO2021128959A1 (en) 2019-12-26 2020-09-12 Power transmission section thermal stability power limit interval identification method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911366425.3A CN111106603B (en) 2019-12-26 2019-12-26 Method and system for identifying limit interval of thermal stability power of power transmission section

Publications (2)

Publication Number Publication Date
CN111106603A true CN111106603A (en) 2020-05-05
CN111106603B CN111106603B (en) 2021-11-30

Family

ID=70424851

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911366425.3A Active CN111106603B (en) 2019-12-26 2019-12-26 Method and system for identifying limit interval of thermal stability power of power transmission section

Country Status (2)

Country Link
CN (1) CN111106603B (en)
WO (1) WO2021128959A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111900730A (en) * 2020-07-20 2020-11-06 国电南瑞科技股份有限公司 Online calculation method and device for thermal stability quota interval of power transmission channel
WO2021128959A1 (en) * 2019-12-26 2021-07-01 国电南瑞科技股份有限公司 Power transmission section thermal stability power limit interval identification method and system

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113783236B (en) * 2021-09-03 2023-10-20 南方电网科学研究院有限责任公司 Real-time optimal scheduling method, device and medium for high-proportion new energy power system
CN114024312A (en) * 2021-11-16 2022-02-08 国网江苏省电力有限公司电力科学研究院 Method and system for improving reactive support capability, storage medium and computing equipment
CN114243800B (en) * 2021-11-17 2024-04-02 国网浙江省电力有限公司 Accident standby energy storage power station access point and access capacity optimization coordination configuration method
CN114759559B (en) * 2022-06-16 2022-09-06 中国能源建设集团湖南省电力设计院有限公司 Power grid stability limit analysis method
CN117353312A (en) * 2023-09-15 2024-01-05 国家电网有限公司华东分部 Power grid thermal stability quota checking method, device, equipment and medium
CN117094537B (en) * 2023-10-19 2024-01-05 国网浙江省电力有限公司丽水供电公司 Power grid planning method, device, electronic equipment and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413043A (en) * 2013-08-09 2013-11-27 中国南方电网有限责任公司 Solving method for power system multi-dimensional space cross section thermal stability limit boundary
CN107067126A (en) * 2016-11-25 2017-08-18 中国南方电网有限责任公司 It is a kind of based on power flow transfer than thermally-stabilised key transmission channel recognition method
US20180034317A1 (en) * 2016-07-26 2018-02-01 Schweitzer Engineering Laboratories, Inc. Microgrid Power Flow Monitoring and Control
CN109936134A (en) * 2019-04-03 2019-06-25 国电南瑞科技股份有限公司 A kind of thermostabilization passway for transmitting electricity power limit section recognition methods and equipment

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101969202B (en) * 2010-11-23 2012-10-17 国网电力科学研究院 Method for recognizing transmission section power limit interval of power system transient stability
CN108964012A (en) * 2017-05-22 2018-12-07 武汉大学 The security constraint Unit Combination dual blank-holder decomposed based on Benders
CN111106603B (en) * 2019-12-26 2021-11-30 国电南瑞科技股份有限公司 Method and system for identifying limit interval of thermal stability power of power transmission section

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413043A (en) * 2013-08-09 2013-11-27 中国南方电网有限责任公司 Solving method for power system multi-dimensional space cross section thermal stability limit boundary
US20180034317A1 (en) * 2016-07-26 2018-02-01 Schweitzer Engineering Laboratories, Inc. Microgrid Power Flow Monitoring and Control
CN107067126A (en) * 2016-11-25 2017-08-18 中国南方电网有限责任公司 It is a kind of based on power flow transfer than thermally-stabilised key transmission channel recognition method
CN109936134A (en) * 2019-04-03 2019-06-25 国电南瑞科技股份有限公司 A kind of thermostabilization passway for transmitting electricity power limit section recognition methods and equipment

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
徐岩等: "基于N-1原则的输电断面热稳定功率极限评估", 《电力系统及其自动化学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021128959A1 (en) * 2019-12-26 2021-07-01 国电南瑞科技股份有限公司 Power transmission section thermal stability power limit interval identification method and system
CN111900730A (en) * 2020-07-20 2020-11-06 国电南瑞科技股份有限公司 Online calculation method and device for thermal stability quota interval of power transmission channel
CN111900730B (en) * 2020-07-20 2021-09-28 国电南瑞科技股份有限公司 Online calculation method and device for thermal stability quota interval of power transmission channel

Also Published As

Publication number Publication date
WO2021128959A1 (en) 2021-07-01
CN111106603B (en) 2021-11-30

Similar Documents

Publication Publication Date Title
CN111106603B (en) Method and system for identifying limit interval of thermal stability power of power transmission section
Thangaiah et al. Delay‐dependent stability analysis of power system considering communication delays
CN104578411A (en) Multi-source data safety verification and failure locking method after PMU (Phasor Measurements Units) data's accessing to AGC (Automatic Generation Control)
CN110137948B (en) Method and device for evaluating sub-synchronous oscillation risk of doubly-fed wind power plant
CN109829246B (en) Line parameter identification method based on parameter comprehensive suspicion
CN112330488A (en) Power grid frequency situation prediction method based on transfer learning
CN109873418B (en) Strong correlation power transmission section identification method and system in transient power angle stable mode
CN106056305B (en) Power generation system reliability rapid evaluation method based on state clustering
CN114936544A (en) Data processing method and device, electronic equipment and storage medium
CN113094887B (en) Optimization method and device for frequency-shifting electromagnetic transient simulation and electronic equipment
CN111756031B (en) Power grid operation trend estimation method and system
CN107169644A (en) A kind of power distribution network safe operation management-control method
CN114977222B (en) Analysis method for broadband oscillation influence factors and propagation paths of power system
CN108551164B (en) Voltage stability control method and device for direct-current micro-grid
CN114386510A (en) Method and system for identifying measurement errors of power system
CN115907526A (en) Power grid voltage sag evaluation method and system, storage medium and computing device
CN113517693A (en) Weak voltage transient stability node analysis method and device
Mahmoudabadi et al. An application of CHA to concurrent short-term transmission expansion & reactive power planning
CN113113908A (en) Time domain analysis method and system suitable for frequency response of modern large power grid
CN113076628A (en) Analysis method and system suitable for frequency safety index of modern large power grid
CN104037810A (en) Improved risk assessment electricity generation and load power adjustment method
CN111162519A (en) Power grid topological structure node voltage sensing method and device
CN111049158A (en) Method and system for determining broadband oscillation stability of power system based on spectrum radius
CN113221301B (en) Method and device for generating stability criterion of power electronic power station
CN111900730B (en) Online calculation method and device for thermal stability quota interval of power transmission channel

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant