CN113836748B - Geomagnetic disturbance-oriented disaster prevention and management method, device, storage medium and equipment - Google Patents
Geomagnetic disturbance-oriented disaster prevention and management method, device, storage medium and equipment Download PDFInfo
- Publication number
- CN113836748B CN113836748B CN202111210979.1A CN202111210979A CN113836748B CN 113836748 B CN113836748 B CN 113836748B CN 202111210979 A CN202111210979 A CN 202111210979A CN 113836748 B CN113836748 B CN 113836748B
- Authority
- CN
- China
- Prior art keywords
- geomagnetic
- power system
- bias
- governance
- current
- 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.)
- Active
Links
- 230000002265 prevention Effects 0.000 title claims abstract description 29
- 238000007726 management method Methods 0.000 title claims abstract description 28
- 230000006698 induction Effects 0.000 claims abstract description 94
- 238000000034 method Methods 0.000 claims abstract description 56
- 230000008569 process Effects 0.000 claims abstract description 33
- 238000013178 mathematical model Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 238000005315 distribution function Methods 0.000 claims abstract description 29
- 230000006870 function Effects 0.000 claims description 55
- 238000005457 optimization Methods 0.000 claims description 37
- 238000005070 sampling Methods 0.000 claims description 19
- 230000008859 change Effects 0.000 claims description 18
- 230000005684 electric field Effects 0.000 claims description 18
- 238000004088 simulation Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000007935 neutral effect Effects 0.000 description 6
- 239000002699 waste material Substances 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000003990 capacitor Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005358 geomagnetic field Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Classifications
-
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
Abstract
The application discloses a geomagnetic disturbance-oriented disaster prevention and management method, a geomagnetic disturbance-oriented disaster prevention and management device, a storage medium and geomagnetic induction current samples are calculated according to line induction voltages by simulating line induction voltages of all nodes based on geoelectric conductivity data. Based on each geomagnetic induction current sample, a probability distribution function of the geomagnetic induction current is determined. And analyzing the bias magnetic data to obtain bias magnetic current of the node. Substituting each bias current into the probability distribution function to obtain each actual magnetic induction current. Substituting each actual geomagnetic induction current into the magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model. The states of all the governance devices in the power system are adjusted so that the total number of governance devices in an open state is equal to the optimal solution shown by the calculation result, and therefore, under the condition that the number of governance devices participating in the governance process of the power system is ensured to meet governance requirements, the excessive governance devices are prevented from participating in the governance process.
Description
Technical Field
The present application relates to the field of electric power, and in particular, to a method, an apparatus, a storage medium, and a device for disaster prevention and management for geomagnetic disturbance.
Background
The increase in reactive power loss of the transformer caused by the geomagnetic induction current (Geomagnetically Induced Current, GIC) will destroy the reactive power balance of the system, causing damage to the transformer and even a blackout accident. The global synchronism of the geomagnetic storm gives rise to the transformer GIC reactive losses to full network mass-transmission, and in addition, the uncertainty of the geomagnetic field strength at the time of the magnetic storm gives rise to the geomagnetic storm to randomness. In real life, the running state of the power system gradually approaches the stability limit of the power system, and reactive power fluctuation caused by the magnetotelluric storm with strong randomness can bring great impact to the power system, so that large-scale accidents of the power system are induced.
Currently, regarding the possibility that the GIC may threaten the safe and stable operation of the system when the geomagnetic storm occurs, researchers propose three main treatment methods of installing a resistor, a capacitor and a resistor and a capacitor at the neutral point of the transformer. However, the configuration of these abatement devices often depends on the maximum value of the GIC in a certain geomagnetic storm event, and the volatility of the GIC is not considered, which causes the problem that the abatement device has too strong pertinence, and the problems of "oversupply" and "undersupply" are easy to occur, so that the stable requirement of the power system cannot be satisfied.
Disclosure of Invention
The application provides a disaster prevention management method, device, storage medium and equipment for geomagnetic disturbance, and aims to avoid excessive management devices to participate in a management process under the condition that the number of management devices participating in the management process of a power system is ensured to meet management requirements.
In order to achieve the above object, the present application provides the following technical solutions:
a disaster prevention and treatment method for geomagnetic disturbance comprises the following steps:
based on the pre-collected earth conductivity data, simulating to obtain line induction voltage of each node in the power system when the earth storm is influenced;
according to the line induced voltage of each node, calculating to obtain a geomagnetic induced current sample of each node;
determining a probability distribution function of the geomagnetic induction current based on each geomagnetic induction current sample;
analyzing the bias magnetic data of the power system to obtain bias magnetic current of each node;
substituting each bias current into the probability distribution function to obtain each actual magnetic induction current;
substituting each actual geomagnetic induction current into a pre-constructed magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model; the magnetic bias governance mathematical model comprises a multi-objective optimization function taking the number of governance devices participating in the governance process of the power system as an optimization objective; the calculation result comprises an optimal solution of the multi-objective optimization function;
and adjusting the states of all the governance devices in the power system so that the total number of governance devices in an on state is equal to the optimal solution.
Optionally, based on the pre-collected earth conductivity data, the simulating obtains the line induction voltage of each node in the power system when the line induction voltage is affected by the geomagnetic storm, including:
constructing a ground conductivity model based on the pre-collected ground conductivity data;
simulating and calculating an induced ground electric field, a line path and a line length of geomagnetic storms generated under the ground conductivity model by using a plane wave method;
and calculating to obtain the line induction voltage of each node in the power system when the power system is influenced by the geomagnetic storm according to the induction ground electric field, the line path and the line length.
Optionally, the determining a probability distribution function of the geomagnetic induction current based on each geomagnetic induction current sample includes:
calculating the difference value between each geomagnetic induction current sample and a preset threshold value to obtain each magnetic bias current sample;
fitting the occurrence frequency of each bias current sample according to various probability density functions to obtain fitting indexes of various probability density functions;
and marking the probability density function with the minimum fitting index as a probability distribution function of the geomagnetic induction current.
Optionally, the analyzing the bias magnetic data of the power system to obtain bias magnetic currents of the nodes includes:
collecting bias magnetic data of the power system according to a sampling time scale; the sampling time scale is determined according to the reference change rate and the relative information loss rate of the geomagnetic induction current sample under different time scales;
and analyzing the bias magnetic data to obtain bias magnetic current of each node.
A disaster prevention and management device for geomagnetic disturbance, comprising:
the voltage simulation unit is used for simulating and obtaining line induction voltage of each node in the power system when the power system is influenced by geomagnetic storms based on the pre-collected geodetic conductivity data;
the sample calculation unit is used for calculating geomagnetic induction current samples of the nodes according to the line induction voltage of each node;
a function determining unit, configured to determine a probability distribution function of the geomagnetic induction current based on each geomagnetic induction current sample;
the current analysis unit is used for analyzing the magnetic bias data of the power system to obtain the magnetic bias current of each node;
the current calculation unit is used for substituting each bias current into the probability distribution function to obtain each actual magnetically induced current;
the target optimization unit is used for substituting each actual geomagnetic induction current into a pre-constructed magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model; the magnetic bias governance mathematical model comprises a multi-objective optimization function taking the number of governance devices participating in the governance process of the power system as an optimization objective; the calculation result comprises an optimal solution of the multi-objective optimization function;
and the quantity adjusting unit is used for adjusting the states of all the governance devices in the power system so that the total number of governance devices in an on state is equal to the optimal solution.
Optionally, the voltage simulation unit is specifically configured to:
constructing a ground conductivity model based on the pre-collected ground conductivity data;
simulating and calculating an induced ground electric field, a line path and a line length of geomagnetic storms generated under the ground conductivity model by using a plane wave method;
and calculating to obtain the line induction voltage of each node in the power system when the power system is influenced by the geomagnetic storm according to the induction ground electric field, the line path and the line length.
Optionally, the function determining unit is specifically configured to:
calculating the difference value between each geomagnetic induction current sample and a preset threshold value to obtain each magnetic bias current sample;
fitting the occurrence frequency of each bias current sample according to various probability density functions to obtain fitting indexes of various probability density functions;
and marking the probability density function with the minimum fitting index as a probability distribution function of the geomagnetic induction current.
Optionally, the current analysis unit is specifically configured to:
collecting bias magnetic data of the power system according to a sampling time scale; the sampling time scale is determined according to the reference change rate and the relative information loss rate of the geomagnetic induction current sample under different time scales;
and analyzing the bias magnetic data to obtain bias magnetic current of each node.
A computer-readable storage medium comprising a stored program, wherein the program performs the method of disaster prevention management for geomagnetic disturbances.
A geomagnetic disturbance-oriented disaster prevention and management device, comprising: a processor, a memory, and a bus; the processor is connected with the memory through the bus;
the memory is used for storing a program, and the processor is used for running the program, wherein the disaster prevention and management method facing geomagnetic disturbance is executed when the program runs.
According to the technical scheme, based on the pre-collected earth conductivity data, line induction voltage of each node in the power system is obtained through simulation when the earth storm is influenced. And calculating geomagnetic induction current samples of the nodes according to the line induction voltage of the nodes. Based on each geomagnetic induction current sample, a probability distribution function of the geomagnetic induction current is determined. And analyzing the bias magnetic data of the power system to obtain bias magnetic current of each node. Substituting each bias current into the probability distribution function to obtain each actual magnetic induction current. Substituting each actual geomagnetic induction current into a pre-constructed magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model. The bias magnetic governance mathematical model comprises a multi-objective optimization function taking the number of governance devices participating in the governance process of the power system as an optimization target, and the calculation result comprises an optimal solution of the multi-objective optimization function. The states of all the governance devices in the power system are adjusted so that the total number of governance devices in the on state is equal to the optimal solution. Based on the scheme shown in the application, the bias current is substituted into the probability distribution function to obtain the actual magnetic induction current, and the actual magnetic induction current is substituted into the multi-objective optimization function to obtain the optimal solution of the number of the treatment devices participating in the treatment process of the power system, and the states of the treatment devices in the power system are adjusted to ensure that the total number of the treatment devices in the on state is equal to the optimal solution, so that under the condition that the number of the treatment devices participating in the treatment process of the power system is ensured to meet the treatment requirement, the excessive treatment devices are prevented from participating in the treatment process, and the resource waste is prevented.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic flow chart of a disaster prevention and treatment method for geomagnetic disturbance, which is provided by an embodiment of the present application;
FIG. 2 is a frequency distribution histogram according to an embodiment of the present application;
FIG. 3 is a schematic flow chart of another disaster prevention and control method for geomagnetic disturbance according to an embodiment of the present application;
fig. 4 is a schematic architecture diagram of a disaster prevention and management device for geomagnetic disturbance according to an embodiment of the present application.
Detailed Description
The following description of the embodiments of the present application will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all, of the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
As shown in fig. 1, a flow chart of a disaster prevention and treatment method for geomagnetic disturbance provided by an embodiment of the present application includes the following steps:
s101: and constructing a ground conductivity model based on the pre-collected ground conductivity data.
The specific implementation process of constructing the earth conductivity model based on the earth conductivity data is common knowledge familiar to those skilled in the art, and will not be described herein.
S102: and (3) using a plane wave method to simulate and calculate the induction ground electric field, the line path and the line length of the geomagnetic storm which occurs under the ground conductivity model.
The geomagnetic storm induced ground electric field generated under the ground conductivity model is calculated by using a plane wave method, which is known as a person skilled in the art, and is not described herein.
S103: and calculating to obtain the line induction voltage of each node in the power system when the power system is influenced by geomagnetic storm according to the induction ground electric field of the geomagnetic storm, the line path and the line length.
The specific calculation process of the line induction voltage of each node in the power system when the power system is affected by the geomagnetic storm is calculated according to the induction ground electric field of the geomagnetic storm, the line path and the line length, as shown in a formula (1).
In formula (1), U represents a line induced voltage, R represents a line path, E represents an induced ground field, and L represents a line length.
As can be seen from the above formula (1), the induced voltage of the line is substantially an integral value of the induced field of the local magnetic storm along the preset line, and since the induced field of the local magnetic storm is constant, the integral value is considered to be related to the geographic positions of the two ends of the preset line (i.e. the line start point and the line end point), i.e. the induced voltage of the line is the sum of the first voltage (i.e. the induced voltage of the line in the north-south direction) and the second voltage (i.e. the induced voltage of the line in the east-west direction), so that the above formula (1) can be converted into the formula (2) to facilitate calculation of the induced voltage of the line.
U＝E N L N +E E L E (2)
In formula (2), E N Represents the value of the north-oriented earth electric field, E E Represents the eastern earth electric field value, L N Represents the north length of the line L E Representing the east length of the line.
S104: and calculating a GIC sample of each node according to the line induced voltage of each node.
The specific calculation process of the GIC sample of each node is calculated according to the line induced voltage of each node, as shown in formula (3).
I node ＝(I+YZ) -1 J (3)
In formula (3), I node GIC samples representing any node in the power system, node being the index of the node, I being an N identity matrix, N representing the total number of nodes, Y representing the admittance matrix, Z representing the ground impedance matrix, J representing the equivalent current source in the line, and J= [ J ] 1 ,J 2 ,...,J N ]，[J 1 ,J 2 ,...,J N ]Any element in the above is shown in formula (4).
In formula (4), R ij Represents the equivalent resistance between node i and node j, and i=1, 2,..n, and j=1, 2,..n, V ij Representing the equivalent voltage between node i and node j, and the equivalent voltage of the node is equal to the line induced voltage of the node.
S105: and calculating the difference value between each GIC sample and a preset threshold value to obtain each bias current sample.
Wherein each bias current sample should be greater than 0 under the influence of a geomagnetic storm.
S106: fitting the occurrence frequency of each bias current sample according to various probability density functions to obtain fitting indexes of various probability density functions.
The probability density function includes a gamma distribution, a lognormal distribution, a Rayleigh distribution, an exponential distribution, and a Weber distribution.
The expression of the fitting index is shown in formula (5).
In formula (5), I Fitting Representing a fitting index, M represents the number of groupings in the frequency distribution histogram of the bias current samples (i.e., the total number of bias current samples of different values), generally, the frequency distribution histogram is generated based on the frequency of occurrence of the bias current samples of different values, i=1, 2, M, for probability density used in fitting processDegree function, y i Represents->Fitting the value obtained by the probability density function, < >>Representing the height of the ith square column (i.e. the frequency of occurrence of arbitrarily valued bias current samples) in the frequency histogram>Representing the center position of the frequency distribution histogram. A specific frequency distribution histogram can be seen in fig. 2.
S107: and marking the probability density function with the minimum fitting index as the probability distribution function of the GIC.
The smaller the fitting index is, the higher the fitting accuracy of the representative probability density function is, and the probability distribution of the real GIC is fitted.
S108: and collecting bias magnetic data of the power system according to the sampling time scale, and analyzing the bias magnetic data to obtain bias magnetic current of each node.
The sampling time scale is determined according to the reference change rate and the relative information loss rate of the GIC sample under different time scales.
It should be noted that, according to different time scales, part of bias magnetic data is lost in the process of collecting bias magnetic data of the power system, in order to prevent excessive amount of lost data, the lost bias magnetic data under different time scales needs to be measured, so as to obtain data loss under different time scales, and the time scale with the least data loss is selected as the sampling time scale, so that accuracy of bias magnetic data sampling is improved.
In the embodiment of the present application, the data loss amount may be determined by two indexes, namely, a reference change rate and a relative information loss rate, and the smaller the reference change rate and the smaller the relative information loss rate, the smaller the representative data loss amount. The reference change rate and the relative information loss rate of the GIC sample under each time scale are calculated in advance, and the time scale with the minimum reference change rate and the minimum relative information loss rate is selected from each time scale to be used as the sampling time scale.
Generally, the increase of the time scale can cause the information loss of the stability and the continuity of the GIC variation, so in order to measure the influence of the increase of the time scale on the stability of the GIC fluctuation, an index of a reference change rate is provided, and the smaller the reference change rate is, the smaller the stability of the time scale on the GIC fluctuation is, and specifically, the expression of the reference change rate is shown as a formula (6).
In the formula (6), S k Represents the reference rate of change, k is the index of the time scale, S B Representing a preset reference load, ΔGIC is the first order difference of the GIC sample, U B Representing a preset reference voltage.
In the bias magnetic data obtained under different time scale sampling, the average first-order difference component and the probability that the difference component is 0 can reflect the fluctuation condition of the bias magnetic data. Therefore, comparing the bias magnetic data under different time scales with the standard data, and calculating to obtain the relative information loss rate under different time scales, wherein the smaller the relative information loss rate is, the smaller the fluctuation of the bias magnetic data under the time scales is, and specifically, taking the time scales of 5s and 1min as an example, the expression of the relative information loss rate is shown as a formula (7).
In formula (7), λ represents the relative information loss rate, Δ 5s Represents the average first order difference, delta, at a time scale of 5S 1min Represents the average first order difference, delta, at a time scale of 1min 0 Representing a preset reference first-order differential quantity,representing the probability of a time scale of 5s at a time difference component of 0,/for>GIC represents the probability that the time difference component is 0 at a time scale of 1min 0 Representing a preset reference difference amount.
Specifically, taking three time scales of 5s, 1min and 5min as examples, the reference change rate and the relative information loss rate of the GIC sample at the three time scales are shown in table 1. As can be seen from table 1, the probability that the reference change rate is less than 1% is maximum when the time scale is 5S, that is, the reference change rate is minimum when the time scale is 5S, and the relative information loss rate is minimum when the time scale is 5S. For this, 5s is determined as the sampling time scale.
TABLE 1
The foregoing table 1 is presented for illustration only.
S109: substituting each bias current into the probability distribution function to obtain each actual GIC.
S110: substituting each actual GIC into a pre-constructed magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model.
The bias magnetic governance mathematical model comprises a multi-objective optimization function taking the number of governance devices participating in the governance process of the power system as an optimization objective. In addition, the calculation result includes an optimal solution of the multi-objective optimization function.
Specifically, the expression of the multi-objective optimization function is shown in formula (8).
In formula (8), I ALL Representing the neutral point of the transformerThe sum of bias currents, E represents a preset expected value, N represents the number of neutral point grounding transformers, I i Represents the actual GIC flowing through the neutral point of the ith transformer, R s Represents the sum of resistance values of adjustable resistors in the treatment device, M represents the number of treatment devices participating in the treatment process of the power system, R i Representing the resistance of the ith abatement device.
In terms of the effect, the purpose of disaster prevention and management is to hope that the smaller and the better the neutral point magnetic bias current, the more management devices are used, and in terms of the operation and maintenance, the fewer and the better the management devices are, the lower the magnetic bias current cannot be reduced. Therefore, the two objective functions shown in the formula (8) are contradictory and mutually restricted.
In addition, comprehensively considering the magnetic bias influence of the GIC and the grounding electrode, the proposed transformer magnetic bias treatment scheme cannot use the unified treatment standard because the action range and the influence degree of the GIC and the grounding electrode are different, and the magnetic bias current of the power grid caused by the two conditions needs to be treated separately, so that the constraint conditions of the two objective functions in the formula (8) also need to be set respectively. Specifically, I i The constraint of (2) is as shown in formula (9), R i The constraint of (2) is shown in formula (10).
I i ≤18A (9)
0＜R i ≤5Ω (10)
Optionally, in order to improve the accuracy of the calculation result, a non-dominant ranking genetic algorithm (Nondominated Sorting Genetic Algorithm, NSGA-ii) may be used to optimize the bias treatment mathematical model. Specifically, the mathematical model of bias magnetic treatment after NSGA-II optimization is shown as a formula (11).
In formula (11), I max Represents the maximum value of the sum of bias currents, |A i I represents I ALL Is a value of (2).
Specifically, the process for optimizing the bias magnetic treatment mathematical model by using NSGA-II comprises the following steps:
1. inputting grid bias data (i.e., the actual GIC flowing through the neutral point of each transformer);
2. constructing a target optimization model (namely, formula (11));
3. setting relevant parameters of NSGA-II, and initializing a population (namely, the number of the treatment devices represented by individuals in the population);
4. non-dominant ordering is carried out on each individual in the population to obtain a first generation subgroup;
5. iteratively executing a preset step to obtain a non-dominant solution; wherein, the preset steps are as follows: combining parent individuals and offspring individuals to obtain new individuals; performing non-dominant sorting and crowding degree calculation on each new individual; generating a new parent population; selecting, crossing and mutating individuals in the population;
6. parameters of the overall target optimization model are demodulated according to non-dominance.
It should be noted that the foregoing specific implementation is merely illustrative.
S111: the states of all the governance devices in the power system are adjusted so that the total number of governance devices in an on state is equal to the optimal solution.
The state of each governance device in the power system is adjusted, so that the total number of governance devices in an open state is equal to an optimal solution, the number of governance devices participating in the governance process of the power system can be ensured to meet governance requirements (namely, the sum of bias currents of the power system reaches a minimum value), and excessive governance devices are prevented from participating in the governance process, so that the problems of over governance and under governance in the governance process are effectively solved.
In summary, by using the scheme shown in this embodiment, under the condition that the number of treatment devices participating in the treatment process of the power system is ensured to meet the treatment requirement, the excessive treatment devices are prevented from participating in the treatment process, and resource waste is prevented.
It should be noted that S101 mentioned in the foregoing embodiment is an optional implementation manner of the disaster prevention and management method for geomagnetic disturbance shown in the application. In addition, S102 mentioned in the foregoing embodiment is also an optional implementation manner of the disaster prevention and management method for geomagnetic disturbance shown in the application. For this reason, the flow mentioned in the above embodiment can be summarized as the method shown in fig. 3.
As shown in fig. 3, a flow chart of another disaster prevention and treatment method for geomagnetic disturbance is provided in the embodiment of the present application, including the following steps:
s301: based on the pre-collected earth conductivity data, line induction voltage of each node in the power system when the earth storm is influenced is simulated.
S302: and calculating geomagnetic induction current samples of the nodes according to the line induction voltage of the nodes.
S303: based on each geomagnetic induction current sample, a probability distribution function of the geomagnetic induction current is determined.
S304: and analyzing the bias magnetic data of the power system to obtain bias magnetic current of each node.
S305: substituting each bias current into the probability distribution function to obtain each actual magnetic induction current.
S306: substituting each actual geomagnetic induction current into a pre-constructed magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model.
The bias magnetic governance mathematical model comprises a multi-objective optimization function taking the number of governance devices participating in the governance process of the power system as an optimization objective. The calculation result includes an optimal solution of the multi-objective optimization function.
S307: the states of all the governance devices in the power system are adjusted so that the total number of governance devices in the on state is equal to the optimal solution.
In summary, the bias current is substituted into the probability distribution function to obtain the actual magnetically induced current, and the actual magnetically induced current is substituted into the multi-objective optimization function to obtain the optimal solution of the number of the abatement devices participating in the abatement process of the electric power system, and the states of the abatement devices in the electric power system are adjusted so that the total number of the abatement devices in the on state is equal to the optimal solution, thereby avoiding the participation of excessive abatement devices in the abatement process and preventing the resource waste under the condition that the number of abatement devices participating in the abatement process of the electric power system meets the abatement requirement.
Corresponding to the method for disaster prevention and management for geomagnetic disturbance provided in the embodiment of the present application, the embodiment of the present application further provides a device for disaster prevention and management for geomagnetic disturbance.
As shown in fig. 4, an architecture diagram of a disaster prevention and management device for geomagnetic disturbance according to an embodiment of the present application includes:
the voltage simulation unit 100 is configured to simulate and obtain a line induced voltage of each node in the power system when the power system is affected by geomagnetic storms based on the collected data of the earth conductivity in advance.
The voltage simulation unit 100 is specifically configured to: constructing a ground conductivity model based on the pre-collected ground conductivity data; simulating and calculating an induced ground electric field, a line path and a line length of geomagnetic storms generated under a ground conductivity model by using a plane wave method; and calculating to obtain the line induction voltage of each node in the power system when the power system is influenced by geomagnetic storms according to the induction ground electric field, the line path and the line length.
The sample calculation unit 200 is configured to calculate a geomagnetic induction current sample of each node according to the line induction voltage of each node.
The function determining unit 300 is configured to determine a probability distribution function of the geomagnetic induction current based on the respective geomagnetic induction current samples.
The function determining unit 300 specifically is configured to: calculating the difference value between each geomagnetic induction current sample and a preset threshold value to obtain each magnetic bias current sample; fitting the occurrence frequency of each bias current sample according to various probability density functions to obtain fitting indexes of various probability density functions; and marking the probability density function with the minimum fitting index as a probability distribution function of the geomagnetic induction current.
The current analysis unit 400 is configured to analyze the bias magnetic data of the power system to obtain bias magnetic currents of each node.
The current analysis unit 400 specifically is configured to: collecting bias magnetic data of a power system according to a sampling time scale; the sampling time scale is determined according to the reference change rate and the relative information loss rate of the geomagnetic induction current sample under different time scales; and analyzing the bias magnetic data to obtain bias magnetic current of each node.
The current calculation unit 500 is configured to substitute each bias current into the probability distribution function to obtain each actual magnetically induced current.
The target optimization unit 600 is configured to substitute each actual geomagnetic induction current into a pre-constructed magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model; the magnetic bias governance mathematical model comprises a multi-objective optimization function taking the number of governance devices participating in the governance process of the power system as an optimization objective; the calculation result includes an optimal solution of the multi-objective optimization function.
And a quantity adjusting unit 700 for adjusting the states of the respective abatement devices in the power system so that the total number of abatement devices in the on state is equal to the optimal solution.
In summary, the bias current is substituted into the probability distribution function to obtain the actual magnetically induced current, and the actual magnetically induced current is substituted into the multi-objective optimization function to obtain the optimal solution of the number of the abatement devices participating in the abatement process of the electric power system, and the states of the abatement devices in the electric power system are adjusted so that the total number of the abatement devices in the on state is equal to the optimal solution, thereby avoiding the participation of excessive abatement devices in the abatement process and preventing the resource waste under the condition that the number of abatement devices participating in the abatement process of the electric power system meets the abatement requirement.
The application also provides a computer readable storage medium, wherein the computer readable storage medium comprises a stored program, and the program executes the geomagnetic disturbance-oriented disaster prevention and management method.
The application also provides a disaster prevention and treatment device for geomagnetic disturbance, which comprises: a processor, a memory, and a bus. The processor is connected with the memory through a bus, the memory is used for storing a program, and the processor is used for running the program, wherein the method for preventing and managing geomagnetic disturbance, provided by the application, is executed when the program runs, and comprises the following steps:
based on the pre-collected earth conductivity data, simulating to obtain line induction voltage of each node in the power system when the earth storm is influenced;
according to the line induced voltage of each node, calculating to obtain a geomagnetic induced current sample of each node;
determining a probability distribution function of the geomagnetic induction current based on each geomagnetic induction current sample;
analyzing the bias magnetic data of the power system to obtain bias magnetic current of each node;
substituting each bias current into the probability distribution function to obtain each actual magnetic induction current;
substituting each actual geomagnetic induction current into a pre-constructed magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model; the magnetic bias governance mathematical model comprises a multi-objective optimization function taking the number of governance devices participating in the governance process of the power system as an optimization objective; the calculation result comprises an optimal solution of the multi-objective optimization function;
and adjusting the states of all the governance devices in the power system so that the total number of governance devices in an on state is equal to the optimal solution.
Specifically, based on the above embodiment, the simulating, based on the pre-collected earth conductivity data, the line induction voltage of each node in the power system when being affected by the geomagnetic storm includes:
constructing a ground conductivity model based on the pre-collected ground conductivity data;
simulating and calculating an induced ground electric field, a line path and a line length of geomagnetic storms generated under the ground conductivity model by using a plane wave method;
and calculating to obtain the line induction voltage of each node in the power system when the power system is influenced by the geomagnetic storm according to the induction ground electric field, the line path and the line length.
Specifically, on the basis of the foregoing embodiment, the determining, based on each geomagnetic induction current sample, a probability distribution function of the geomagnetic induction current includes:
calculating the difference value between each geomagnetic induction current sample and a preset threshold value to obtain each magnetic bias current sample;
fitting the occurrence frequency of each bias current sample according to various probability density functions to obtain fitting indexes of various probability density functions;
and marking the probability density function with the minimum fitting index as a probability distribution function of the geomagnetic induction current.
Specifically, on the basis of the foregoing embodiment, the analyzing the bias magnetic data of the power system to obtain bias magnetic currents of the nodes includes:
collecting bias magnetic data of the power system according to a sampling time scale; the sampling time scale is determined according to the reference change rate and the relative information loss rate of the geomagnetic induction current sample under different time scales;
and analyzing the bias magnetic data to obtain bias magnetic current of each node.
The functions described in the methods of the present application, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computing device readable storage medium. Based on such understanding, a portion of the embodiments of the present application that contributes to the prior art or a portion of the technical solution may be embodied in the form of a software product stored in a storage medium, comprising several instructions for causing a computing device (which may be a personal computer, a server, a mobile computing device or a network device, etc.) to perform all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
In this specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, so that the same or similar parts between the embodiments are referred to each other.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the application. Thus, the present application is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (6)
1. The disaster prevention and treatment method for geomagnetic disturbance is characterized by comprising the following steps:
based on the pre-collected earth conductivity data, simulating to obtain line induction voltage of each node in the power system when the earth storm is influenced;
according to the line induced voltage of each node, calculating to obtain a geomagnetic induced current sample of each node;
calculating the difference value between each geomagnetic induction current sample and a preset threshold value to obtain each magnetic bias current sample;
fitting the occurrence frequency of each bias current sample according to various probability density functions to obtain fitting indexes of various probability density functions;
marking the probability density function with the minimum fitting index as a probability distribution function of geomagnetic induction current;
collecting bias magnetic data of the power system according to a sampling time scale; the sampling time scale is determined according to the reference change rate and the relative information loss rate of the geomagnetic induction current sample under different time scales;
analyzing the bias magnetic data to obtain bias magnetic current of each node;
substituting each bias current into the probability distribution function to obtain each actual magnetic induction current;
substituting each actual geomagnetic induction current into a pre-constructed magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model; the magnetic bias governance mathematical model comprises a multi-objective optimization function taking the number of governance devices participating in the governance process of the power system as an optimization objective; the calculation result comprises an optimal solution of the multi-objective optimization function;
and adjusting the states of all the governance devices in the power system so that the total number of governance devices in an on state is equal to the optimal solution.
2. The method of claim 1, wherein simulating the line induced voltages at each node in the power system when affected by the geomagnetic storm based on the pre-collected earth conductivity data, comprises:
constructing a ground conductivity model based on the pre-collected ground conductivity data;
simulating and calculating an induced ground electric field, a line path and a line length of geomagnetic storms generated under the ground conductivity model by using a plane wave method;
and calculating to obtain the line induction voltage of each node in the power system when the power system is influenced by the geomagnetic storm according to the induction ground electric field, the line path and the line length.
3. The utility model provides a disaster prevention and treatment device towards geomagnetic disturbance which characterized in that includes:
the voltage simulation unit is used for simulating and obtaining line induction voltage of each node in the power system when the power system is influenced by geomagnetic storms based on the pre-collected geodetic conductivity data;
the sample calculation unit is used for calculating geomagnetic induction current samples of the nodes according to the line induction voltage of each node;
the function determining unit is used for calculating the difference value between each geomagnetic induction current sample and a preset threshold value to obtain each bias current sample; fitting the occurrence frequency of each bias current sample according to various probability density functions to obtain fitting indexes of various probability density functions; marking the probability density function with the minimum fitting index as a probability distribution function of geomagnetic induction current;
the current analysis unit is used for collecting magnetic bias data of the power system according to a sampling time scale; the sampling time scale is determined according to the reference change rate and the relative information loss rate of the geomagnetic induction current sample under different time scales; analyzing the bias magnetic data to obtain bias magnetic current of each node;
the current calculation unit is used for substituting each bias current into the probability distribution function to obtain each actual magnetically induced current;
the target optimization unit is used for substituting each actual geomagnetic induction current into a pre-constructed magnetic bias treatment mathematical model to obtain a calculation result output by the magnetic bias treatment mathematical model; the magnetic bias governance mathematical model comprises a multi-objective optimization function taking the number of governance devices participating in the governance process of the power system as an optimization objective; the calculation result comprises an optimal solution of the multi-objective optimization function;
and the quantity adjusting unit is used for adjusting the states of all the governance devices in the power system so that the total number of governance devices in an on state is equal to the optimal solution.
4. A device according to claim 3, characterized in that the voltage simulation unit is specifically adapted to:
constructing a ground conductivity model based on the pre-collected ground conductivity data;
simulating and calculating an induced ground electric field, a line path and a line length of geomagnetic storms generated under the ground conductivity model by using a plane wave method;
and calculating to obtain the line induction voltage of each node in the power system when the power system is influenced by the geomagnetic storm according to the induction ground electric field, the line path and the line length.
5. A computer-readable storage medium, characterized in that the computer-readable storage medium comprises a stored program, wherein the program executes the geomagnetic disturbance oriented disaster prevention management method of any of claims 1 to 2.
6. A disaster prevention and management device for geomagnetic disturbance, comprising: a processor, a memory, and a bus; the processor is connected with the memory through the bus;
the memory is used for storing a program, and the processor is used for running the program, wherein the program executes the disaster prevention and management method for geomagnetic disturbance according to any one of claims 1-2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111210979.1A CN113836748B (en) | 2021-10-18 | 2021-10-18 | Geomagnetic disturbance-oriented disaster prevention and management method, device, storage medium and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111210979.1A CN113836748B (en) | 2021-10-18 | 2021-10-18 | Geomagnetic disturbance-oriented disaster prevention and management method, device, storage medium and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113836748A CN113836748A (en) | 2021-12-24 |
| CN113836748B true CN113836748B (en) | 2024-01-30 |
Family
ID=78965216
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111210979.1A Active CN113836748B (en) | 2021-10-18 | 2021-10-18 | Geomagnetic disturbance-oriented disaster prevention and management method, device, storage medium and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113836748B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117350102A (en) * | 2023-09-20 | 2024-01-05 | 国网上海市电力公司 | Subway electric power system management method, device, equipment and readable storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20150075907A (en) * | 2013-12-26 | 2015-07-06 | (주)에스이티시스템 | Method and apparatus for predicting geomagnetically induced current |
| CN104821560A (en) * | 2015-04-14 | 2015-08-05 | 华北电力大学 | Novel interphase distance III segment protecting method in consideration of influence of geomagnetically induced current |
| CN108010693A (en) * | 2017-12-26 | 2018-05-08 | 华北电力大学 | The suppressing method and device of transformer DC magnetic bias |
| CN112653135A (en) * | 2020-12-14 | 2021-04-13 | 国网内蒙古东部电力有限公司检修分公司 | Optimization method for governing geomagnetic storm power grid disasters by adopting small resistors |
| CN113034002A (en) * | 2021-03-26 | 2021-06-25 | 国网江苏省电力有限公司电力科学研究院 | Method for analyzing small-disturbance voltage stability of power system by geomagnetic storm |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11422179B2 (en) * | 2019-06-20 | 2022-08-23 | Afshin REZAEI ZARE | System and method for geomagnetic disturbance determination for power systems |
-
2021
- 2021-10-18 CN CN202111210979.1A patent/CN113836748B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20150075907A (en) * | 2013-12-26 | 2015-07-06 | (주)에스이티시스템 | Method and apparatus for predicting geomagnetically induced current |
| CN104821560A (en) * | 2015-04-14 | 2015-08-05 | 华北电力大学 | Novel interphase distance III segment protecting method in consideration of influence of geomagnetically induced current |
| CN108010693A (en) * | 2017-12-26 | 2018-05-08 | 华北电力大学 | The suppressing method and device of transformer DC magnetic bias |
| CN112653135A (en) * | 2020-12-14 | 2021-04-13 | 国网内蒙古东部电力有限公司检修分公司 | Optimization method for governing geomagnetic storm power grid disasters by adopting small resistors |
| CN113034002A (en) * | 2021-03-26 | 2021-06-25 | 国网江苏省电力有限公司电力科学研究院 | Method for analyzing small-disturbance voltage stability of power system by geomagnetic storm |
Non-Patent Citations (2)
| Title |
|---|
| 抑制变压器地磁感应电流的电容隔直装置安装位置优化;刘春明 等;电力系统自动化(16);全文 * |
| 抑制电网地磁感应电流的电容电阻优化布置方法;刘春明 等;《电网技术》;第42卷(第5期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113836748A (en) | 2021-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Visconti et al. | Measurement-based load modeling using transfer functions for dynamic simulations | |
| CN107453897B (en) | Node importance evaluation method and device, electronic equipment and storage medium | |
| CN110717619A (en) | Multi-scale space-time load prediction method and system for bottom-up power distribution network | |
| CN113836748B (en) | Geomagnetic disturbance-oriented disaster prevention and management method, device, storage medium and equipment | |
| Mikic | Variance-based energy loss computation in low voltage distribution networks | |
| Raposo et al. | Optimal meter placement algorithm for state estimation in power distribution networks | |
| Deb et al. | Modified spider monkey optimization-based optimal placement of distributed generators in radial distribution system for voltage security improvement | |
| Xing et al. | Exact solutions for average trapping time of random walks on weighted scale-free networks | |
| Chihota | Extending the Herman-Beta transform for probabilistic load flow analysis of radial feeders | |
| Noghabi et al. | Parameter uncertainty in the optimal coordination of overcurrent relays | |
| CN109829652B (en) | Long-time scale dynamic harmonic responsibility division method | |
| Paruta et al. | Greedy placement of measurement devices on distribution grids based on enhanced distflow state estimation | |
| CN112633665A (en) | Lightning protection decision method for power distribution network based on analytic hierarchy process | |
| Lin et al. | Voltage sag severity analysis based on improved FP-Growth algorithm and AHP algorithm | |
| Wu et al. | Extreme bifurcation interval with deep PV penetration considering correlated net bus load forecast errors | |
| Sushma et al. | Observability constrained GA approach for optimal PMU placement considering zero injection modeling | |
| Arya et al. | Determination of probabilistic risk of voltage collapse using radial basis function (RBF) network | |
| Biscaro et al. | Optimal phasor measurement units placement for fault location on overhead electric power distribution feeders | |
| Guerrero et al. | Adaptation of the clonal selection algorithm to the real-time coordinated Volt/VAr control through a software-in-the-loop strategy | |
| CN108334990B (en) | Reactive power compensation site selection and capacity optimization method and system for large power grid | |
| Gajjala et al. | Congestion management using grey wolf optimization in a deregulated power market | |
| Wu et al. | Voltage stability in a long-distance power transmission system impacted by the geoelectric field due to a geomagnetic disturbance | |
| Bendale et al. | Evaluation of different grounding grid designs for Microgrid | |
| Blagajac et al. | A computational tool to support spatial electric load forecasting | |
| CN111697575B (en) | Flexibility-improved quantitative analysis method, device and system |
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 |