CN112819648A - Power dispatching operation instruction safety verification method and system - Google Patents
Power dispatching operation instruction safety verification method and system Download PDFInfo
- Publication number
- CN112819648A CN112819648A CN202011590065.8A CN202011590065A CN112819648A CN 112819648 A CN112819648 A CN 112819648A CN 202011590065 A CN202011590065 A CN 202011590065A CN 112819648 A CN112819648 A CN 112819648A
- Authority
- CN
- China
- Prior art keywords
- line
- power
- transformer substation
- substation
- operation instruction
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000012795 verification Methods 0.000 title claims description 15
- 230000006854 communication Effects 0.000 claims abstract description 93
- 238000004891 communication Methods 0.000 claims abstract description 92
- 230000005540 biological transmission Effects 0.000 claims abstract description 48
- 238000004364 calculation method Methods 0.000 claims abstract description 41
- 230000008859 change Effects 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 19
- 108010074506 Transfer Factor Proteins 0.000 claims description 14
- 230000035945 sensitivity Effects 0.000 claims description 12
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 3
- 230000001502 supplementing effect Effects 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000011158 quantitative evaluation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 8
- 238000002347 injection Methods 0.000 description 5
- 239000007924 injection Substances 0.000 description 5
- 238000012423 maintenance Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000013524 data verification Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Abstract
The invention discloses a method and a system for safely checking a power dispatching operation instruction, wherein a power grid local power flow calculation model is established according to a local communication network between transformer substations, and the transmission power of a circuit in a communication coil after the dispatching operation instruction is executed is calculated; the change relation of the line power flow before and after the scheduling instruction is executed is researched, the electric power scheduling instruction credibility model is established, and quantitative evaluation of the scheduling instruction safety analysis is realized. The invention overcomes the defect that the existing method can not cope with network attack, and improves the safety level of power dispatching.
Description
Technical Field
The invention relates to the technical field of power systems, in particular to power system safety analysis, and specifically relates to a power dispatching operation instruction safety verification method and system.
Background
With the continuous improvement of the intelligent level of the power system, information and communication technologies are widely applied to the power system. In recent years, power system network attack events frequently occur, and malicious attackers perform data injection attacks by using security holes existing in power information systems. For example, a malicious attacker can destroy the integrity of the data of the power system by modifying the data of the measurement device, and mislead a dispatcher to issue a wrong dispatching operation instruction, thereby causing safety accidents such as disconnection of a line, load shedding and the like.
At present, the execution mode of the power dispatching operation instruction is simple from top to bottom, and an effective verification mechanism does not exist. The dispatcher judges and compiles an operation instruction ticket according to actual measurement data uploaded by the transformer substation, and sends the operation instruction ticket to the transformer substation in advance; and after receiving the operation instruction ticket, the transformer operation and maintenance personnel carry out field operation according to the operation instruction ticket.
The existing power dispatching operation instruction operation method has a plurality of defects: 1) the current power dispatching operation instruction mode does not consider the inaccuracy of judgment of a dispatcher on wrong or incomplete measurement data; 2) in addition, the current scheduling mode does not consider the possibility that the scheduling instruction is tampered in the power communication process; 3) the conventional power dispatching mode is simple from top to bottom, dispatching operation instructions are made on system macro data only by dispatching operators, and risk verification is not carried out on the dispatching operation instructions by substation operation instruction operators through field data. Therefore, the current power dispatching operation instruction operation method is not sufficient for coping with information attacks, has serious potential safety hazards, is easy to cause large-scale power accidents, and urgently needs to invent a new power dispatching operation instruction safety verification method and system.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a method and a system for safely checking an electric power dispatching operation instruction aiming at the defects of the prior art, so that the safety and stability of an electric power system are improved.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows: a safety verification method for power dispatching operation instructions comprises the following steps:
s1, supplementing operation order time and execution instruction reference transformer substation on the basis of related information of the existing power dispatching operation instructionActual load DM,iAnd the actual output G of the power stationM,i(ii) a Establishing a local communication network between the transformer substation and the adjacent transformer substations to realize the actual load D of the transformer substations in the communication circleiAnd generator output GiSharing data;
s2, calculating the change power delta F of the external line of the transformer substation at the boundary of the communication ring according to the current topological structure psi of the power system and the line load flow sensitivity characteristicm；
S3, changing power delta F according to the external line of the boundary transformer substationmActual load of transformer substation DiEstablishing an equivalent load model of the transformer substation and calculating the equivalent load D of the transformer substation i in the communication circlei′；
S4, according to the equivalent load D of the transformer substation i in the communication circlei' local topology psi of power systemloaclCalculating the transmission power F of the communication coil after the electric power operation command is executedl；
S5, according to the transmission power F of the communication coil inner circuit after executing the electric power operation instructionlEstablishing a power dispatching operation instruction credibility model, and calculating a dispatching operation instruction credibility interval IbAnd judging whether the scheduling operation instruction is executed or not.
The existing transformer substation lacks safety verification on scheduling operation instructions, and safety accidents are easily caused by blind or wrong execution of risk instructions. The safety verification of the scheduling operation instruction is realized through the steps S1-S5, the prior study and judgment and early warning of the scheduling instruction safety risk are realized through calculating the physical consequence after the instruction is executed, and the serious potential safety hazard caused by the execution of the high-risk instruction is avoided.
In step S1, a local communication network of the transformer substation and the adjacent transformer substations is established to realize the actual load D of the transformer substation in the communication circleiAnd generator output GiThe data sharing method comprises the following steps: determining a sub-node transformer substation j according to a connection line n-j of the transformer substation n by taking the transformer substation n as a root node; a local communication network is established between the transformer substation n and all the child nodes j, namely a communication ring of the transformer substation n is formed, and actual data D between the transformer substation n and the transformer substation j is realizedn、Dj、Gn、GjSharing of (2); and if the transformer substation j and the transformer substation outside the communication circle have a connecting line, the transformer substation j is a boundary transformer substation.
The conventional load flow calculation method depends on the working condition data of the whole network and is difficult to realize in the weak communication environment of the transformer substation. Aiming at the problem, step S1 realizes high-precision load flow calculation based on local network parameter information and local working condition data by exchanging working condition data of the local communication circle by using the local communication circle of the substation, and solves the problem of the weak communication environment between substations at present.
After step S1, before step S2, the following operations are also performed: according to the actual load D of the transformer substation in the communication circleiAnd message transformer station load data DM,iDifference of test data Δ Di：ΔDi＝DM,i-Di(ii) a When | Δ DiWhen | ≠ 0, will DiUploading the power dispatching center to wait for a second dispatching operation instruction; otherwise, or the second message | Δ DiIf ≠ 0, it proceeds to step S3. The step carries out pre-verification aiming at the working condition data received by the transformer substation, and only when the actual load of the transformer substation is different from the load data in the message, the subsequent safety verification step is started.
The specific implementation process of step S2 includes:
1) calculating reactance sensitivity matrix H of external circuit m of boundary transformer substationm：
Hm＝[Xk-Xh]1×f；
Wherein X is a node reactance matrix; k is a head end node of the line m, namely a boundary transformer substation k in the communication circle of the transformer substation n; h is the end node of line m; xkIs the k-th row of the matrix X, XhF is the number of nodes of the power system, namely the number of substations governed by the power dispatching center;
2) according to the reactance sensitivity matrix HmAnd the actual transmission power F of line aaAnd calculating the power flow transfer factor tau of the external line m after the line a is disconnectedm；
3) According to the power flow transfer factor taumAnd a reactance xmCalculating the transfer power DeltaF of the line mm：
The calculation process is simple, the calculation speed is high, and the reliability is high. In step 2), the power flow transfer factor tau is calculated by the following methodm：
A) According to the actual transmission power F of the line aaSystem node reactance matrix X, line reactance XaCalculating the injected power delta P at two ends of the line an，ΔPj；
B) According to injected power delta P at two ends of the line an、ΔPjAnd a reactance sensitivity matrix H of an external circuit m of the communication coil boundary transformer substationmCalculating the power flow transfer factor of the external line m after the line a is disconnected
Step S2 realizes the load flow calculation without changing the power grid topological structure through the calculation of the load flow transfer factor, avoids the multiple load flow calculation caused by the change of the power grid topological structure, and reduces the calculation complexity of the whole method.
The specific implementation process of step S3 includes:
a) according to boundary substation k communication circle in-connection circuit z transmission power FzK generator set output G of connection boundary transformer substationkBoundary substation actual load DkCalculating transmission power of external line of boundary transformer substation
Wherein omegal,kA set of connection lines for boundary substation k;
b) according to the actual load D of the boundary substation kkBoundary changeTransmission power F of external line of power stationmThe variable power delta F of the external line m after the line a is disconnectedmCalculating the equivalent load D of the boundary substation kk′：Dk′＝Dk+Fm+ΔFm；
c) Calculating equivalent load D of transformer substation j in n communication circle of transformer substationj′：Dj′＝Dj(ii) a Wherein D isjIs the actual load of substation j.
And step S3, calculating the equivalent load of the substation node, and the method has the advantages of realizing independent calculation of the load flow of the communication circle, avoiding the dependence of a load flow calculation model on the working condition data of the whole network, ensuring the accuracy of the load flow calculation and accelerating the speed of the load flow calculation.
In step S4, the in-coil transmission power F after the execution of the operation command is calculated by the following formulal：
Is the phase angle of the head end node of the line l,
is the terminal node phase angle of line l; x is the number oflIs the impedance of line i.
The specific implementation process of step S5 includes:
I) according to the transmission power F of the line in the communication coil after the line a is disconnectedlAnd line transmission capacity Fr,lCalculating a load factor R of each linel；
II) according to the load factor R of each linelDetermining a maximum load rate line b in a communication circle after the instruction is executed: rb＝max(Rl)；l∈Ωnl(ii) a Wherein omeganlThe method comprises the steps of collecting lines in a communication circle of a transformer substation n;
III) calculating a 95% safety confidence interval I of the line b according to the line transmission power data of the previous T days, namely the historical transmission power of the line b of T1b。
Step S5, the load rate of the circuit in the communication circle after the instruction execution is calculated, the historical data is fully mined, and the safety confidence interval of the circuit is determined, so that the safety risk assessment aiming at the operation instruction is more reliable, and the identification degree of the risk instruction is improved.
After step III), the following operations are performed: according to I calculated in step III)bWhen R isb∈IbWhen the operation instruction is executed; otherwise, the execution is refused; wherein R isbIs the load factor of line b. According to the step, pre-warning of the risk instruction is realized according to the pre-calculated safety confidence interval, and serious safety risk caused by wrong execution of the risk instruction is prevented.
The invention also provides a power dispatching operation instruction safety verification system, which comprises computer equipment; the computer device is configured or programmed for performing the steps of the above-described method.
Compared with the prior art, the invention has the beneficial effects that:
(1) the invention considers the attack characteristic of an attacker on the malicious tampering of the dispatching instruction, provides a power dispatching instruction safety multi-check communication mechanism of a power dispatching center-transformer substation, overcomes the defect that the existing dispatching instruction safety analysis method cannot cope with information attack, and improves the safety level of the power dispatching operation instruction.
(2) The invention researches a system local power flow calculation method under incomplete information, establishes a local power flow calculation model, ensures the accuracy of safety risk of local data verification scheduling operation instructions, overcomes the defect of overhigh large-scale point-to-point communication construction cost among traditional transformer substations, and improves the investment proficiency of power companies.
(3) The method analyzes the change relation of the line power flow before and after the scheduling operation instruction is executed, establishes the electric power scheduling operation instruction credibility model by utilizing historical data, realizes quantitative evaluation rationality of the scheduling operation instruction safety, and overcomes the serious defect that the existing method lacks physical consequence correlation analysis.
Drawings
Fig. 1 is a flowchart of a method for verifying safety of an operating instruction in power dispatching according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating a format of a scheduling operation order message according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a local communication network between substations in an embodiment of the present invention;
FIG. 4 is a block diagram of a system for verifying safety of power dispatching operation commands in an embodiment of the present invention;
FIG. 5 is a system block diagram of a substation load calculation module in an embodiment of the present invention;
fig. 6 is a system unit diagram of a local power flow calculation module in the embodiment of the invention.
Detailed Description
Fig. 1 is a flowchart of a method for verifying safety of an electric power scheduling operation instruction according to an embodiment of the present invention, which includes the following specific steps:
step S1: on the basis of the related information such as the command number, the serial number, the ticket filling date, the operation ending date, the operation task and the like of the existing power dispatching operation command message, the operation ticket time and the execution command are supplemented to refer to the actual load D of the transformer substationM,iAnd the actual output G of the power stationM,iFig. 2 shows a specific diagram of a new message format.
Step S2: establishing local communication network between the transformer substation and the adjacent transformer substations to realize the transformer substation load D in the communication circleiAnd generator output GiAnd (4) sharing data. The operation and maintenance personnel of the transformer substation change the load data D of the transformer substation according to the scheduling operation instructionM,iSharing data with communication circle DiAnd (6) carrying out data verification. And when the data is unreasonable, returning the real data to the dispatching center.
Step S3: calculating the change power delta F of the external line of the communication circle boundary transformer substation (hereinafter referred to as boundary transformer substation) according to the topological structure psi of the current system and the line tidal current sensitivity characteristicm；
Step S4: power delta F is changed according to external lines of boundary transformer substationmActual load of transformer substation DiEstablishing an equivalent load model of the transformer substation and calculating the equivalent load D of the transformer substation i in the communication circlei′；
Step S5:according to equivalent load D of transformer substation i in communication circlei' local topology of sum system psiloaclCalculating the transmission power F of the communication coil after the electric power operation command is executedl；
Step S6: according to the transmission power F of the communication coil internal circuit after the execution of the electric power operation instructionlEstablishing a power dispatching operation instruction credibility model, and calculating a dispatching operation instruction credibility interval IbAnd judging whether the scheduling operation instruction is executed or not.
The transformer substation is taken as an example for explanation, and the specific modeling process and the meaning of parameters in each model are shown in the following examples, which are not specifically explained here.
The following are more specific embodiments of the present invention:
further as a preferred embodiment, step S2 includes:
step S21: and according to the topological structure of the power network, establishing a local communication network between the transformer substations, namely a communication ring of the transformer substation. And selecting the transformer substation n as a root node, and determining a sub-node transformer substation j according to the connection circuit n-j of the transformer substation n. The substation n and all the sub-nodes j establish a local communication network through a special optical fiber network, namely a communication ring of the substation n is formed, and actual data D between the substation n and the substation j is realizedn、Dj、Gn、GjTo share. And if the transformer substation j and the transformer substation outside the communication circle have a connecting line, the transformer substation j is a boundary transformer substation.
In one embodiment, a schematic diagram of a substation 2 communication network is shown in fig. 3. And selecting the transformer substation 2 as a root node, and determining the transformer substations 1, 5, 3 and 4 as child nodes according to the lines 2-1, 2-5, 2-3 and 2-4. The transformer substation 2, the transformer substation 1, the transformer substation 5, the transformer substation 3 and the transformer substation 4 form a local communication network, namely a communication ring of the transformer substation 2, and real-time data sharing of the transformer substation 2, the transformer substation 1, the transformer substation 5, the transformer substation 3 and the transformer substation 4 is realized. Wherein, the transformer substation 5 and the transformer substation 4 have external circuits, and the transformer substation 5 and the transformer substation 4 are boundary transformer substations of the communication ring of the transformer substation 2.
Step S22: according to the fact that the transformer substation is in the communication circleInter load DiAnd message actual load data D of transformer substationM,iAnd checking data difference, wherein the specific formula is as follows:
ΔDi＝DM,i-Di
step S23: when | Δ DiWhen | ≠ 0, the substation operation and maintenance personnel will DiUploading the power dispatching center to wait for a second dispatching operation instruction; otherwise or second message | Δ DiIf | ≠ 0, proceed to step 3.
Further as a preferred embodiment, step S3 includes:
step S31: calculating a reactance sensitivity matrix H of the external circuit m of the boundary transformer substation according to the topological structure psi of the system and the node reactance matrix XmThe concrete formula is as follows:
Hm＝[Xk-Xh]1×f
wherein k is a head end node of the line m, namely a boundary substation k in the communication circle of the substation n; h is the end node of line m; xkIs the k-th row of the matrix X, XhAnd f is the number of system nodes, namely the number of substations governed by the power dispatching center.
Step S32: according to the reactance sensitivity matrix H obtained in S31mAnd the actual transmission power F of line aaAnd calculating the power flow transfer factor tau of the external line m after the line a is disconnectedm。
In one embodiment, as shown in fig. 2, after the lines 2-3 are disconnected, the load flow transfer factor tau of the external lines 4-7, 4-9 of the substation 4 at the boundary of the communication circle is calculated47、τ49And a power flow transfer factor tau of the external lines 5-6 of the boundary substation 556。
Further as a preferred embodiment, step S32 includes:
step S320: according to the actual transmission power F of the line aaSystem node reactance matrix X, line reactance XaCalculating the injected power delta P at two ends of the line an，ΔPjThe concrete formula is as follows:
wherein n is a head end node of the line a, namely a transformer substation n; j is the end node of line a, namely substation j; xnnIs the nth element in the nth row and column of XjjIs the jth row and jth column element in X, XnjIs the nth row and the jth column element in X.
Step S321: according to the injected power delta P between the two ends of the line a in S320n、ΔPjAnd a reactance sensitivity matrix H of the external circuit m of the communication coil boundary substation in S31mCalculating the power flow transfer factor tau of the external line m after the line a is disconnectedmThe concrete formula is as follows:
step S33: according to the power flow transfer factor taumAnd a reactance xmCalculating the transfer power DeltaF of the line mmThe concrete formula is as follows:
further as a preferred embodiment, step S4 includes:
step S41: according to boundary substation k communication circle in-connection circuit z transmission power FzK generator set output G of connection boundary transformer substationkBoundary substation actual load DkTransmission power F of external line of boundary transformer substationmThe calculation formula is as follows:
wherein omegal,kIs the set of connection lines of the border substation k.
Step S42: according to the actual load D of the boundary substation kkOutside of boundary transformer substationTransmission power F of connection linemThe variable power delta F of the external line m after the line a is disconnectedmCalculating the equivalent load D of the boundary substation kk', the specific formula is as follows:
Dk′＝Dk+Fm+ΔFm
step S43: calculating equivalent load D of transformer substation j in n communication circle of transformer substationj', the specific formula is as follows:
Dj′＝Dj
wherein D isjIs the actual load of substation j.
Further as a preferred embodiment, step S5 includes:
step S51: according to the equivalent load D of the transformer substation i obtained in S4i' actual output G of generator setiAnd calculating the injection power of a transformer substation i in the communication circle, wherein the specific formula is as follows:
Pi＝Gi-Di′
step S52: according to the injection power P of the substation i in S51iExecuting local system node reactance matrix X after operation instructionlocalAnd line impedance xlCalculating the transmission power F of the communication coil after the operation command is executedlThe concrete formula is as follows.
θ＝XlocalPlocal
Wherein the content of the first and second substances,
is the phase angle of the head end node of the line l,
is the terminal nodal phase angle of line l, θ is the vector of the composition of all nodal phase angles, xlIs the impedance of line l; plocalInjecting power P for local intra-network substation iiThe vectors of the components.
The steps S2-S5 study the local power flow calculation characteristics in the communication circle and establish a local power flow calculation model in the communication circle. The method avoids large-scale inter-substation communication while ensuring the calculation accuracy, and greatly reduces the communication cost.
Further as a preferred embodiment, step S6 includes:
step S61: according to the transmission power F of the line in the communication coil after the line a is disconnectedlAnd line transmission capacity Fr,lCalculating a load factor R of each linelThe concrete formula is as follows:
step S62: according to the load factor R of each linelDetermining a maximum load rate line b in a communication circle after the instruction is executed, wherein the specific formula is as follows:
Rb＝max(Rl)l∈Ωnl
wherein omeganlIs a line set in a communication circle of the transformer substation n.
In one embodiment, the load factor of the line 2-5 is the greatest after the operating command to disconnect the line 2-3 is performed, so that the line b is the line 2-5.
Step S63: from the line transmission power data of the previous 90 days, i.e. the historical transmission power of 8640 groups of lines b, a 95% safety confidence interval for line b is calculated.
Further preferably, step S63 includes:
step S630: transmitting load rate data R according to line b at 8640 time points in the first 90 daysb,tCalculating the mean value mu of the transmission load rate of the line bbThe concrete formula is as follows:
where t is the time, and t-1 represents the initial time 90 days ago.
In one embodiment, the dispatcher gives the dispatching instruction operation time: when 10:00 is given at 6/20/2020, time t1 is 10:00 given at 3/20/2020.
Step S631: transmitting load rate data R according to line b at 8640 time points in the first 90 daysb,tLine b transmission load factor mean μbCalculating the standard deviation sigma of the transmission load rate of the line bbThe concrete formula is as follows:
step S632: according to the mean value mu of the transmission load rate of the line b in S630bAnd the line b transmission load rate standard deviation σ in S631bCalculating the confidence coefficient of the transmission load rate of the line b as a 95% confidence interval IbThe concrete formula is as follows:
step S64: confidence interval I of 95% of load rate of transmission of line b according to S63bAnd load factor R of line b of S62bAnd judging the rationality of the scheduling instruction. When R isb∈IbWhen the operation instruction is executed; otherwise, execution is denied.
The step S6 researches the incidence relation between the scheduling instruction and the line power flow after the scheduling instruction is executed, establishes a power scheduling instruction confidence model, and overcomes the serious shortcoming that the existing evaluation method lacks instruction analysis by calculating the load rate confidence of the maximum load rate line after the scheduling operation instruction is executed to be a 95% confidence interval and judging the validity of the power scheduling instruction.
Fig. 4 is a schematic structural diagram of a power scheduling operation instruction execution system according to an embodiment of the present invention, where the system is adapted to execute a method according to any embodiment of the present invention, and includes: the system comprises a data reading module 100, a transformer substation load calculation module 200, a local load flow calculation module 300 and an operation instruction judgment module 400.
In the embodiment of the present invention, the data reading module 100, the transformer substation load calculation module 200, the local load flow calculation module 300, and the operation instruction determination module 400 are all modules in a computer device.
The computer equipment can be a microprocessor, and can also be equipment such as an upper computer and the like.
The data reading module 100 is configured to read data, including data such as a power topology structure, data of an existing power generator set of a power grid, an actual load of a substation, an actual load of an execution instruction reference substation, a system node reactance, and a local line reactance.
The transformer substation load calculation module 200 is used for calculating equivalent load of a transformer substation in the communication circle.
The local power flow calculation module 300 is configured to calculate a transmission power of a line in the communication coil after the transformer substation executes the operation instruction.
The operation instruction judging module 400 is configured to judge the reliability of the scheduling operation instruction.
The output end of the data reading module 100 is connected with the input end of the transformer substation load calculation module 200, and is used for inputting a power topological structure, a transformer substation actual load, a system node reactance and a local line reactance.
The output end of the data reading module 100 is connected to the input end of the local power flow calculation module 300, and is used for inputting the power topology and the generator data.
The output end of the transformer substation load calculation module 200 is connected with the input end of the local load flow calculation module 300 and is used for inputting equivalent load of a transformer substation.
The output end of the local power flow calculation module 300 is connected to the input end of the operation instruction judgment module 400, and is used for inputting and calculating the transmission power of the line in the communication coil after the transformer substation executes the operation instruction.
As shown in fig. 5, further, the substation load calculation module 200 includes: a data reading unit 201, a first calculating unit 202, a second calculating unit 203, and a third calculating unit 204.
The output end of the data reading unit 201 is connected to the input end of the first calculating unit 202, and is used for inputting a power line topology structure, a power line reactance, a line actual transmission power, and a system substation node reactance.
The output end of the data reading unit 201 is connected with the input end of the second calculating unit 203, and is used for inputting data of the existing generator set of the power grid, the actual load of the transformer substation, the topological structure of the power line, and the actual transmission power of the line.
The output end of the data reading unit 201 is connected with the input end of the third calculating unit 204, and is used for inputting the actual load of the electric substation and the topological structure of the electric power line.
The output end of the first calculating unit 202 is connected with the input end of the third calculating unit 204, and is used for inputting the boundary substation external line change power.
The output of the second calculating unit 203 is connected to the input of the third calculating unit 204, and is used for inputting the transmission power of the external line of the interface substation.
In one embodiment, the data reading unit 201 reads x, Fl、X、ψ、Di、Gi. After the data reading is finished, the unit transmits the information of each parameter to the first calculating unit 202, the second calculating unit 203 and the third calculating unit 204.
The first calculating unit 202 is configured to calculate a change power Δ F of an external line of the boundary substationm。
In one embodiment, according to x, FlX, calculating Δ Fm. After solving, the unit will compare Δ FmTo the third calculation unit 204.
The second calculating unit 203 is used for calculating the transmission power F of the external line of the boundary substationm。
In one embodiment, according to Gk、Dk、FzCalculating Fm. After solving, the unit will FmTo the third calculation unit 204.
The third calculating unit 204 is configured to calculate an equivalent load D of the substationi′。
In one embodiment, according to Di、Fm、ΔFmCalculating Di'. After solving, the unit will Di' as an output of the substation load calculation module 200.
As shown in fig. 6, further, the local power flow calculation module 300 includes: a data reading unit 301, a first calculating unit 302, a second calculating unit 303. Wherein the content of the first and second substances,
the output end of the data reading unit 301 is connected to the input end of the first calculating unit 302, and is used for equivalent load of a transformer substation, power topology, data of the existing generator set of the power grid, and reactance of a power line.
The output end of the first calculating unit 302 is connected with the input end of the second calculating unit 303, and is used for calculating the injection power of the transformer substation in the communication circle.
In one embodiment, data read unit 301, read Gi、Di′、XlocalAnd x. After the data reading is finished, the unit transmits the information of each parameter to the first calculating unit 302 and the second calculating unit 303.
The first calculating unit 302 is configured to calculate the injection power of the substation i.
In one embodiment, according to Gi、Di', calculating Pi. After solving, the unit will PiTo the second calculation unit 303.
The second calculating unit 303 is configured to calculate a transmission power F of a line in the communication coill。
In one embodiment, according to Xlocal、PiX, calculating Fl. After solving, the unit will FlAs an output of the local power flow calculation module 300.
The operation instruction judging module 400 judges according to FlCalculation of RbAnd judging whether the line load rate confidence coefficient is within a confidence interval of 95 percent or notbAnd the judgment result is used as the output end of the operation instruction judgment module 400.
The power dispatching operation instruction operation system forms a power dispatching instruction safety multi-check communication mechanism of a power dispatching center-transformer substation, realizes safety risk evaluation of the power dispatching operation instruction according to actual data, message contents and a large amount of historical data of the transformer substation in a power communication circle, and overcomes the limitation that the existing method is insufficient in dealing with information physical attack.
Claims (10)
1. A safety verification method for power dispatching operation instructions is characterized by comprising the following steps:
s1, supplementing operation order time and execution order reference substation actual load D on the basis of the related information of the existing power dispatching operation orderM，iAnd the actual output G of the power stationM，i(ii) a Establishing a local communication network between the transformer substation and the adjacent transformer substations to realize the actual load D of the transformer substations in the communication circleiAnd generator output GiSharing data;
s2, calculating the change power delta F of the external line of the transformer substation at the boundary of the communication ring according to the current topological structure psi of the power system and the line load flow sensitivity characteristicm；
S3, changing power delta F according to the external line of the boundary transformer substationmActual load of transformer substation DiEstablishing an equivalent load model of the transformer substation and calculating the equivalent load D of the transformer substation i in the communication circlei′；
S4, according to the equivalent load D of the transformer substation i in the communication circlei' local topology psi of power systemloaclCalculating the transmission power F of the communication coil after the electric power operation command is executedl；
S5, according to the transmission power F of the communication coil inner circuit after executing the electric power operation instructionlEstablishing a power dispatching operation instruction credibility model, and calculating a dispatching operation instruction credibility interval IbAnd judging whether the scheduling operation instruction is executed or not.
2. The method for safely verifying the power dispatching operation instruction as claimed in claim 1, wherein in step S1, the substation and the adjacent substations are establishedLocal communication network for realizing actual load D of transformer substation in communication circleiAnd generator output GiThe data sharing method comprises the following steps: determining a sub-node transformer substation j according to a connection line n-j of the transformer substation n by taking the transformer substation n as a root node; a local communication network is established between the transformer substation n and all the child nodes j, namely a communication ring of the transformer substation n is formed, and actual data D between the transformer substation n and the transformer substation j is realizedn、Dj、Gn、GjSharing of (2); and if the transformer substation j and the transformer substation outside the communication circle have a connecting line, the transformer substation j is a boundary transformer substation.
3. The power dispatching operation instruction safety check method according to claim 1, wherein after step S1 and before step S2, the following operations are further performed: according to the actual load D of the transformer substation in the communication circleiAnd message actual load data D of transformer substationM，iDifference of test data Δ Di：ΔDi＝DM，i-Di(ii) a When | Δ DiWhen | ≠ 0, will DiUploading the power dispatching center to wait for a second dispatching operation instruction; otherwise, or the second message | Δ DiIf ≠ 0, it proceeds to step S3.
4. The method for verifying the safety of the power dispatching operation instruction as claimed in claim 1, wherein the step S2 is implemented by:
1) calculating reactance sensitivity matrix H of external circuit m of boundary transformer substationm：
Hm＝[Xk-Xh]1×f；
Wherein X is a node reactance matrix; k is a head end node of the line m, namely a boundary transformer substation k in the communication circle of the transformer substation n; h is the end node of line m; xkIs the k-th row of the matrix X, XhF is the number of nodes of the power system, namely the number of substations governed by the power dispatching center;
2) according to the reactance sensitivity matrix HmAnd the actual transmission power F of line aaCalculating for out after disconnection of line aPower flow transfer factor tau of connection line mm；
3) According to the power flow transfer factor taumAnd a reactance xmCalculating the transfer power DeltaF of the line mm：
5. The power dispatching operation instruction safety checking method as claimed in claim 4, wherein in step 2), the power flow transfer factor τ is calculated by the following methodm：
A) According to the actual transmission power F of the line aaSystem node reactance matrix X, line reactance XaCalculating the injected power delta P at two ends of the line an，ΔPj；
B) According to injected power delta P at two ends of the line an、ΔPjAnd a reactance sensitivity matrix H of an external circuit m of the communication coil boundary transformer substationmCalculating the power flow transfer factor tau of the external line m after the line a is disconnectedm：
6. The method for verifying the safety of the power dispatching operation instruction as claimed in claim 1, wherein the step S3 is implemented by:
a) according to boundary substation k communication circle in-connection circuit z transmission power FzK generator set output G of connection boundary transformer substationkBoundary substation actual load DkCalculating transmission power F of external line of boundary transformer substationm：
Wherein omegal，kA set of connection lines for boundary substation k;
b) according to the actual load D of the boundary substation kkBoundary transformer substation external line transmission power FmThe variable power delta F of the external line m after the line a is disconnectedmCalculating the equivalent load D of the boundary substation kk′：Dk′＝Dk+Fm+ΔFm；
c) Calculating equivalent load D of transformer substation j in n communication circle of transformer substationj′：Dj′＝Dj(ii) a Wherein D isjIs the actual load of substation j.
7. The method for verifying safety of power dispatching operation command as claimed in claim 1, wherein in step S4, the transmission power F of the communication coil is calculated according to the following formulal：
Is the phase angle of the head end node of the line l,
is the end node phase angle, x, of line llIs the impedance of line i.
8. The safety verification method for the power dispatching operation instruction according to any one of claims 1 to 7, wherein the specific implementation process of the step S5 includes:
I) according to the transmission power F of the line in the communication coil after the line a is disconnectedlAnd line transmission capacity Fr，lCalculating a load factor R of each linel；
II) according to the load factor R of each linelDetermining a maximum load rate line b in a communication circle after the instruction is executed: rb＝max(Rl)；l∈Ωnl(ii) a Wherein omeganlThe method comprises the steps of collecting lines in a communication circle of a transformer substation n;
III) calculation based on the line transmission power data of the previous T days, i.e. the historical transmission power of the line b of T195% safety confidence interval I for line bb。
9. The power dispatching operation instruction safety checking method according to claim 8, wherein after the step III), the following operations are performed: according to I calculated in step III)bWhen R isb∈IbWhen the operation instruction is executed; otherwise, the execution is refused; wherein R isbIs the load factor of line b.
10. A power dispatching operation instruction safety verification system is characterized by comprising computer equipment; the computer device is configured or programmed for carrying out the steps of the method according to one of claims 1 to 9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011590065.8A CN112819648B (en) | 2020-12-29 | 2020-12-29 | Power dispatching operation instruction safety verification method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011590065.8A CN112819648B (en) | 2020-12-29 | 2020-12-29 | Power dispatching operation instruction safety verification method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112819648A true CN112819648A (en) | 2021-05-18 |
| CN112819648B CN112819648B (en) | 2024-02-13 |
Family
ID=75854461
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011590065.8A Active CN112819648B (en) | 2020-12-29 | 2020-12-29 | Power dispatching operation instruction safety verification method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112819648B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014173131A1 (en) * | 2013-04-23 | 2014-10-30 | 国家电网公司 | Large power grid overall situation on-line integrated quantitative evaluation method based on response |
| CN110543696A (en) * | 2019-08-15 | 2019-12-06 | 国电南瑞科技股份有限公司 | Method for enabling unmodeled small unit to participate in clearing and safety check of power market |
| CN110867970A (en) * | 2019-11-25 | 2020-03-06 | 广州供电局有限公司 | Power grid dispatching operation command operation method, device and system |
| CN111130116A (en) * | 2019-12-14 | 2020-05-08 | 贵州电网有限责任公司 | Scheduling operation power flow checking method based on key topology change item identification |
| CN111641596A (en) * | 2020-05-11 | 2020-09-08 | 湖南大学 | Power network information physical potential safety hazard assessment method and system and power system |
| US20200293627A1 (en) * | 2019-03-13 | 2020-09-17 | General Electric Company | Method and apparatus for composite load calibration for a power system |
| CN111950900A (en) * | 2020-08-12 | 2020-11-17 | 国网山东省电力公司烟台供电公司 | Power system source network load storage intrinsic safety risk assessment method |
-
2020
- 2020-12-29 CN CN202011590065.8A patent/CN112819648B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014173131A1 (en) * | 2013-04-23 | 2014-10-30 | 国家电网公司 | Large power grid overall situation on-line integrated quantitative evaluation method based on response |
| US20200293627A1 (en) * | 2019-03-13 | 2020-09-17 | General Electric Company | Method and apparatus for composite load calibration for a power system |
| CN110543696A (en) * | 2019-08-15 | 2019-12-06 | 国电南瑞科技股份有限公司 | Method for enabling unmodeled small unit to participate in clearing and safety check of power market |
| CN110867970A (en) * | 2019-11-25 | 2020-03-06 | 广州供电局有限公司 | Power grid dispatching operation command operation method, device and system |
| CN111130116A (en) * | 2019-12-14 | 2020-05-08 | 贵州电网有限责任公司 | Scheduling operation power flow checking method based on key topology change item identification |
| CN111641596A (en) * | 2020-05-11 | 2020-09-08 | 湖南大学 | Power network information physical potential safety hazard assessment method and system and power system |
| CN111950900A (en) * | 2020-08-12 | 2020-11-17 | 国网山东省电力公司烟台供电公司 | Power system source network load storage intrinsic safety risk assessment method |
Non-Patent Citations (4)
| Title |
|---|
| 严康等: "大规模空调负荷在恶意控制下的配电网风险分析", 智慧电力, vol. 48, no. 8 * |
| 李辉等: "基于J2EE的节能发电调度信息管理系统", 微计算机信息 * |
| 王艳等: "具有校验功能的地区调度操作票专家系统的研制", 电力系统保护与控制, vol. 37, no. 21 * |
| 陈俊全;叶航超;宋弦;肖倩宏;张耀;安成;: "基于关键拓扑变动项辨识的调度操作潮流校核方法", 机械与电子, no. 01 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112819648B (en) | 2024-02-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103631240B (en) | Heating industry heat death theory Information application service management platform based on Internet of Things | |
| CN101425681A (en) | Relay protection on-line check method for electric grid system | |
| CN104811433B (en) | The distributed Internet of things system and implementation method of a kind of C/S frameworks | |
| CN107346917A (en) | Secondary device based on virtual machine integrates joint debugging method | |
| CN107819323B (en) | The online O＆M method and system of secondary equipment of intelligent converting station based on visual presentation | |
| CN107706889B (en) | Based on the online O＆M method and system of secondary equipment of intelligent converting station intelligently compared | |
| CN106257367A (en) | A kind of nuclear power DCS platform test device based on analogue system and method for testing | |
| CN109919451B (en) | Digital power grid system and method based on real power grid topology analysis and safety knowledge base | |
| CN108470386A (en) | A kind of cable run intelligent inspection system based on GPS and cloud platform | |
| CN102591223A (en) | Intelligentized diagnostic system of substation-operation operation tickets | |
| CN107705547A (en) | Treat crossing site intelligent operation mode in a kind of passenger-cargo rolling dress harbour based on BIM technology | |
| CN107368972A (en) | It is a kind of that work platform and its application method are sent based on BIM technology | |
| CN109639756A (en) | A kind of terminal network incidence relation is shown and equipment accesses real-time monitoring system | |
| CN107563520A (en) | A kind of intelligent substation secondary safety measure method for visualizing | |
| CN105680441B (en) | Simulation based on stimulated replacement technology is to a device | |
| CN103426041A (en) | Radiation source real-time monitoring and managing system and method based on Internet of things | |
| CN110717736A (en) | Construction progress intelligent tracking method and system for building | |
| CN110611312A (en) | On-line simulation system and method for panoramic data of substation bay level | |
| CN112819648A (en) | Power dispatching operation instruction safety verification method and system | |
| CN103716207A (en) | Detection method and device for IEC61850 protocol | |
| CN106682742A (en) | Real-time data acquisition and analysis method based on Internet of Things | |
| CN108803508A (en) | Power Plant Thermal control algolithm application system and method | |
| CN206331274U (en) | A kind of equipment long-distance service system | |
| CN207249677U (en) | A kind of construction site safety monitoring system integrated based on RFID and BIM | |
| CN108322334A (en) | A kind of power communication system reliability estimation method |
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 |