CN118074058B - Intelligent distribution room safety control method and system - Google Patents
Intelligent distribution room safety control method and system Download PDFInfo
- Publication number
- CN118074058B CN118074058B CN202410452296.4A CN202410452296A CN118074058B CN 118074058 B CN118074058 B CN 118074058B CN 202410452296 A CN202410452296 A CN 202410452296A CN 118074058 B CN118074058 B CN 118074058B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- component
- period
- harmonic component
- residual
- 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
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 33
- 230000002950 deficient Effects 0.000 claims abstract description 12
- 238000001514 detection method Methods 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims description 18
- 238000000354 decomposition reaction Methods 0.000 claims description 13
- 230000002159 abnormal effect Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 13
- 230000000694 effects Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 230000001276 controlling effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004804 winding Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000007667 floating Methods 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 239000012774 insulation material Substances 0.000 description 1
- 230000036244 malformation Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013021 overheating Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
Abstract
The invention relates to the technical field of power distribution circuit management and control, in particular to an intelligent power distribution room safety management and control method and system, wherein the method comprises the following steps: collecting an original current signal; acquiring a fundamental component and a harmonic component; calculating harmonic distortion coefficients; further acquiring the jump disturbance degree of the harmonic component; acquiring a residual fundamental component and a residual harmonic component; acquiring a harmonic intensity coefficient; further obtaining the distortion index of the intensity of the defective harmonic of each original current signal in each period; acquiring harmonic intensity distortion consistency coefficients of all transformers; and detecting the working state of the transformer according to the harmonic intensity distortion consistency coefficient and a detection algorithm. The invention aims to solve the problem that the prior distribution room safety control method is too dependent on a safety threshold value, so that system false alarm or missing alarm is easy to occur.
Description
Technical Field
The invention relates to the technical field of power distribution circuit management and control, in particular to an intelligent power distribution room safety management and control method and system.
Background
The distribution room is a key element in the power system, its main function being to receive power from substations or power plants and distribute it to end users. The distribution room is typically composed of a transformer for converting a high voltage to a safe voltage level for home or business use, a distribution cabinet, control equipment, and safety equipment. The transformer is a part of the distribution room, which is easy to have potential safety hazards, such as short circuit of transformer windings, poor grounding of iron cores and the like, which can lead to unstable output voltage and influence the safe operation of the circuit system of the distribution room.
The existing distribution room safety control method is characterized in that a plurality of data parameters are synthesized for analysis through detecting the distribution room environment parameters and the equipment states, so that a system is too complicated, and then whether potential safety hazards exist or not is judged according to a preset safety threshold. The running environment and conditions of the distribution room may change with time, excessive parameter setting may increase the burden of processing data by the system, and incorrect setting of the safety threshold value may easily cause misinformation or missing report of the system, which affects the reliability and effectiveness of the system.
Disclosure of Invention
In order to solve the technical problems, the invention aims to provide an intelligent power distribution room safety management and control method and system, and the adopted technical scheme is as follows:
in a first aspect, an embodiment of the present invention provides a method for controlling security of an intelligent power distribution room, including the following steps:
collecting original current signals of each transformer in each period;
Acquiring fundamental wave components and harmonic wave components of each original current signal by using a sliding window iterative DFT algorithm; acquiring harmonic distortion coefficients of each subsequence in each period according to the harmonic components and the amplitude of the original current signal; acquiring the jump disturbance degree of the harmonic component of each subsequence in each period according to the harmonic distortion coefficient and the amplitude of the harmonic component; using a decomposition technique to obtain a residual fundamental component and a residual harmonic component of the harmonic component in each sub-sequence in each period; acquiring harmonic intensity coefficients of each subsequence of the harmonic component in each period according to the amplitude of the residual harmonic component; acquiring the defective harmonic intensity distortion index of each original current signal in each period according to the relation among the harmonic component jump disorder degree, the harmonic intensity coefficient and each subsequence in the harmonic component; acquiring harmonic intensity distortion consistency coefficients of all transformers according to the defective harmonic intensity distortion index, the residual fundamental wave component and the fundamental wave component of the transformers;
and detecting the working state of the transformer according to the harmonic intensity distortion consistency coefficient and a detection algorithm.
Further, the obtaining the fundamental component and the harmonic component of each original current signal by using the sliding window iterative DFT algorithm includes:
processing the original current signals of each period by using a sliding window iterative DFT algorithm to obtain fundamental wave components of the original current signals;
And subtracting the amplitude of the original current signal at the same moment from the amplitude of the fundamental component to obtain the harmonic component of the original current signal.
Further, the obtaining harmonic distortion coefficients of each sub-sequence in each period according to the harmonic component and the amplitude of the original current signal includes:
Dividing the original current signal, the fundamental wave component and the harmonic component of each period into a first preset number of subsequences with equal length respectively;
For each sub-sequence of each period, calculating the square of the amplitude of the harmonic component at each sampling time in the sub-sequence as a first square value, calculating the square of the amplitude of the original current signal at each sampling time in the sub-sequence as a second square value, calculating the sum of all the first square values in the sub-sequence as a first sum, calculating the sum of all the second square values in the sub-sequence as a second sum, and taking the square root of the ratio of the first sum to the second sum as the harmonic distortion coefficient of each sub-sequence of each period.
Further, the obtaining the harmonic component jump disorder degree of each subsequence in each period according to the harmonic distortion coefficient and the amplitude of the harmonic component includes:
For each subsequence of each period, calculating the difference value between each amplitude value in the subsequence of the harmonic component and the amplitude value of the first sampling time, calculating the average value of all amplitude values in the subsequence of the harmonic component, calculating the absolute value of the ratio of the difference value to the average value, calculating the sum value of all the absolute values in the subsequence of the harmonic component, calculating the ratio of the harmonic distortion coefficient of the subsequence to the number of sampling times in the subsequence as a first ratio, and taking the product of the first ratio and the sum value as the jump turbulence degree of the harmonic component of each subsequence in each period.
Further, the obtaining the residual fundamental component and the residual harmonic component of the harmonic component in each sub-sequence in each period by using a decomposition technology includes:
for each period, taking the subsequence of the harmonic component as the input of a decomposition algorithm, and acquiring a second preset number of subband components;
and calculating correlation coefficients of each sub-band component and the sub-sequence of the harmonic components, taking the sub-band component with the largest correlation number as a residual fundamental component of the harmonic components, and taking sub-band components except the residual fundamental component as residual harmonic components.
Further, the obtaining the harmonic intensity coefficient of each sub-sequence of the harmonic component in each period according to the amplitude of the residual harmonic component includes:
for each residual harmonic component of the harmonic component, taking the maximum amplitude value and the minimum amplitude value of the residual harmonic component as an interval upper limit and an interval lower limit to construct an interval, dividing the interval into a third preset number of equal-length subintervals, and taking the sum of the probabilities of the amplitudes of all sampling moments in each subinterval in the residual harmonic component as the interval amplitude probability of each subinterval;
For each period, calculating a correlation coefficient between the residual harmonic component and the harmonic component of each sub-sequence of the harmonic component, calculating a sum value of the correlation coefficient and a preset adjusting parameter as a third sum value, calculating a ratio of an average value of interval amplitude probabilities of all sub-intervals in the residual harmonic component of each sub-sequence of the harmonic component to the third sum value, and calculating an average value of all the ratios in each sub-sequence of the harmonic component as a harmonic intensity coefficient of each sub-sequence of the harmonic component in each period.
Further, the obtaining the distortion index of the defective harmonic intensity of each original current signal in each period according to the relation among the harmonic component jump disorder degree, the harmonic intensity coefficient and each subsequence in the harmonic component includes:
For each period, calculating the sum value of the harmonic component jump disorder degree and the harmonic intensity coefficient of each subsequence as a fourth sum value, calculating the DTW distance between each subsequence in the harmonic components and the residual fundamental wave component of each subsequence, calculating the sum value of the DTW distance and a preset adjusting parameter as a fifth sum value, calculating the ratio of the fourth sum value to the fifth sum value, and taking the average value of all the ratios in each period as the distortion index of the residual harmonic intensity of each original current signal in each period.
Further, the obtaining the harmonic intensity distortion consistency coefficient of each transformer according to the harmonic intensity distortion index, the residual fundamental component and the fundamental component of the transformer comprises the following steps:
The distortion indexes of the incomplete harmonic intensities of each period of each transformer are arranged according to the time sequence to form a distortion index sequence of each transformer;
Adding the residual fundamental wave components of the transformers and the amplitudes of the fundamental wave components at the corresponding sampling moments to obtain approximate fundamental wave components of the transformers, and obtaining fitting curves of the approximate fundamental wave components by using a curve fitting algorithm;
and calculating the mean square error between the approximate fundamental component of each transformer and the fitting curve of the approximate fundamental component, calculating the DTW distance between each transformer and the distortion index sequences of other transformers, calculating the sum of all the DTW distances of each transformer and the mean square error, and taking the absolute value of the reciprocal of the sum as the harmonic intensity distortion consistency coefficient of each transformer.
Further, the detecting the working state of the transformer according to the harmonic intensity distortion consistency coefficient and the detection algorithm comprises the following steps:
Taking the harmonic intensity distortion consistency coefficient of each transformer as the input of the COF algorithm, outputting an abnormal factor of the harmonic intensity distortion consistency coefficient, and if the abnormal factor of the harmonic intensity distortion consistency coefficient of the transformer is larger than a preset abnormal threshold value, indicating that the transformer has a safety problem; on the contrary, the transformer has no safety problem.
In a second aspect, an embodiment of the present invention further provides an intelligent power distribution room security management and control system, including a memory, a processor, and a computer program stored in the memory and running on the processor, where the processor implements the steps of any one of the methods described above when executing the computer program.
The invention has at least the following beneficial effects:
The invention separates fundamental wave component and harmonic wave component from original current signal by sliding window iteration DFT algorithm, calculates harmonic wave component jump disorder degree based on fluctuation characteristic of harmonic wave component, and has effect of judging whether original signal data sequence distribution is regular. Further, a wavelet packet decomposition is adopted to extract a residual fundamental component and a residual harmonic component from the harmonic component, and a defective harmonic intensity distortion index is constructed based on the residual harmonic component, so that the influence degree of harmonic on a transformer circuit system is reflected; further, the harmonic intensity distortion consistency coefficient is constructed by approximating the sine wave fitting degree of the fundamental wave component and the incomplete harmonic intensity distortion index sequence, so that the working state of the transformer in the distribution room can be accurately judged; finally, according to the harmonic intensity distortion consistency coefficient, the possible early faults are identified by combining a COF algorithm, and the defect that the system false alarm or omission is easily caused by the fact that the existing distribution room safety control method is too dependent on a safety threshold value is overcome.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions and advantages of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are only some embodiments of the invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart illustrating steps of a method for controlling security of an intelligent distribution room according to an embodiment of the present invention;
Fig. 2 is a flowchart of harmonic intensity distortion consistency coefficient acquisition.
Detailed Description
In order to further describe the technical means and effects adopted by the invention to achieve the preset aim, the following detailed description refers to specific implementation, structure, characteristics and effects of an intelligent power distribution room safety control method and system according to the invention by combining the accompanying drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" means that the embodiments are not necessarily the same. Furthermore, the particular features, structures, or characteristics of one or more embodiments may be combined in any suitable manner.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
The invention provides a method and a system for controlling safety of an intelligent power distribution room, which are specifically described below with reference to the accompanying drawings.
Referring to fig. 1, a flowchart of steps of a method for controlling security of an intelligent power distribution room according to an embodiment of the invention is shown, the method includes the following steps:
and S001, collecting current data of the transformer in the distribution room, and preprocessing the data.
Transformers are critical components of the distribution room, the safety of which is of paramount importance. It is responsible for regulating the voltage, ensuring efficient transfer of power between different voltage levels. However, in the process of handling high voltage and current, if internal materials are damaged, serious electric shock accidents may be caused. In addition, overheating of the transformer not only threatens the safety of the equipment, but also increases the risk of fire, and thus, this embodiment protects against these potential hazards by detecting and analyzing the current in the vicinity of the transformer.
In order to accurately collect current data of the output end of the transformer, a current transformer is installed near each transformer to collect the current data, the time interval of data collection is set to be 2ms, the current data of each transformer is obtained, and the duration of each collection period is set to be 1 second. The acquisition process is inevitably interfered by the instrument and external environment factors, and in order to eliminate the influence of the interference on the current data, a neighbor ordering algorithm (Sorted-Neighborhood Method) is adopted to repeatedly clean the data. The input of the algorithm is collected current data, and the output is cleaned current data. The neighbor ranking algorithm is a well-known technique, and the specific process is not described in detail.
So far, the cleaned current data is obtained, and the cleaned current data of each transformer is used as an original current signal of each transformer in each period.
Step S002, constructing a harmonic component jump disorder degree based on the up-and-down fluctuation characteristics of the harmonic component in the original current signal, analyzing the residual harmonic component obtained by further decomposition, constructing a defective harmonic intensity distortion index, and constructing a harmonic intensity distortion consistency coefficient based on the fundamental component and the working states of different transformers.
Transformers in distribution rooms are directly affected by nonlinear loads, such as rectifiers, frequency converters, switching power supplies, etc., which, when operated, produce currents of non-sinusoidal waveforms. The rectifying circuits of these devices typically employ large capacitors for filtering, resulting in severe malformation of the ac side input current waveform, thereby producing a large amount of harmonic current; furthermore, the electromagnetic coupling characteristics of the transformer may lead to the generation of harmonic voltages. The core saturation nonlinear characteristic of the transformer can generate harmonic voltage under specific conditions, thereby affecting the current waveform.
For transformers, harmonic currents can directly lead to increased losses in the transformer core and windings. The core and windings may overheat due to additional losses caused by harmonics, which may accelerate aging of the insulation material, reduce the life of the transformer, and may cause insulation failure. Therefore, deep analysis of harmonic currents in a distribution room is critical to ensure safe management.
The original current signal of the I-th transformer is denoted as the original current signal I. In the embodiment, a sliding window iterative DFT algorithm is adopted to extract harmonic components in an original current signal, the input of the sliding window iterative DFT algorithm is the original current signal I in each period, a window function is set as a Hamming window, the window length is set as 20ms, and the output of the algorithm is a fundamental component J. Since the method of extracting the harmonic component of the current by sliding window iterative DFT is a known technique, the specific process is not repeated.
Further, the amplitude I of the original current signal at the same moment is subtracted from the amplitude of the fundamental component J to obtain a harmonic component X of the original current signal.
The fundamental component J approximates a sine wave with regular periodicity, while the amplitude of the harmonic component X is in a fast and irregular floating up and down state. And the greater the degree of up-down floating of the harmonic component X, the greater the degree of distortion of the original current signal I. The original current signal I, the fundamental wave component J and the harmonic component X have the same data number in the same period, and the original current signal I, the fundamental wave component J and the harmonic component X in each period are divided into 10 subsequences with the same length. According to the characteristics, constructing the harmonic component jump disorder degree:
In the method, in the process of the invention, The harmonic component jump disorder degree of the m-th sub-sequence in each period is represented,Representing the harmonic distortion coefficients of the mth sub-sequence in each period,The number of sampling times of the mth sub-sequence in each period is represented,Representing the mth subsequence in the harmonic componentThe amplitude of the individual sampling instants,Representing the mth subsequence in the harmonic componentThe amplitude of the individual sampling instants,Representing the mean of the magnitudes of the mth subsequence in the harmonic component.Representing the mth subsequence of the original current signal IThe amplitude of the sampling instants.
Harmonic distortion coefficientIs an index for measuring the degree of deviation of a current waveform from a pure sine waveform, and is the ratio of the sum of squares of harmonic components to the sum of squares of original current signalsThe larger the weight occupied by the harmonic component in the original current signal, the larger the influence of the harmonic on the original current signal,The larger the value, the more severe the distortion of the original current signal. When the higher the difference of values within the harmonic component sub-sequence,The larger the absolute value of (2)Relatively stable and small, and calculate the jump disorder degree of the current harmonic waveThe larger the data sequence distribution, the more irregular the original current signal I.
Furthermore, the sliding window iterative DFT algorithm has the unavoidable defect of extracting harmonic components, and because the method adopted in the current harmonic analysis is to calculate the fundamental current component and then subtract the fundamental component from the original current signal to obtain the harmonic component, the harmonic component X is composed of the sum of all subharmonic components. However, in the distribution room line, the power frequency is not always kept at 50HZ but surrounds 50HZ, and when the power frequency deviates from 50HZ, leakage and a fence effect are caused, so that the amplitude of the fundamental component also deviates. In the transformer line of the distribution room, the harmonic component obtained by subtracting the fundamental component from the original current signal contains a residual fundamental component, resulting in an error in detection.
According to the above features, the present embodiment further analyzes the harmonic component by using a wavelet packet decomposition technique to accurately separate the residual fundamental component and the residual harmonic component, the input of the wavelet packet decomposition algorithm is each sub-sequence of the harmonic component X in each period, the number of decomposition layers is set to 2, the output of the algorithm is 4 sub-band components of each sub-sequence of the harmonic component X in each period, and the 4 sub-band components are respectively denoted as X1, X2, X3, X4. The spearman correlation coefficients of the 4 subband components and the corresponding subsequences of the harmonic component X in each period are calculated separately, the subband component with the spearman correlation coefficient closest to 1 is denoted as X1, since X1 is closest to the trend of the harmonic component X, as residual fundamental component, and X2, X3 and X4 as residual harmonic components. The wavelet packet decomposition algorithm and the spearman correlation coefficient are known techniques, and the specific process is not described in detail.
And then carrying out probability statistics on the amplitude of the residual harmonic component, wherein the calculation method comprises the following steps: taking the residual harmonic component X2 as an example, if the minimum amplitude and the maximum amplitude in the residual harmonic component X2 are a and b respectively, the interval (a, b) is divided intoSub-intervals of equal length, in this embodimentThe sum of the probabilities of the occurrence of the amplitudes of the sampling moments in each subinterval in the residual harmonic component X2 is taken as the interval amplitude probability of the subinterval, wherein the value of the sum is 20. To further analyze the effect of harmonics on the transformer line, a defective harmonic intensity distortion index is constructed:
In the method, in the process of the invention, Representing the harmonic power distortion index of the original current signal I at the period w,The number of sub-sequences within the harmonic component X is represented,The harmonic component jump disorder of the mth subsequence in the period w is represented,The harmonic intensity coefficient representing the m-th sub-sequence of harmonic components within the period w,The number of sub-sequences within the period w is indicated,Representing the mth subsequence in the harmonic component X within the period w,Representing the m-th sub-sequence of harmonic components within the period w,Representing the residual fundamental component resulting from the decomposition of the mth subsequence of the harmonic component X in period w and the mth subsequence of the harmonic component in period wThe calculation of the DTW distance is a well-known technique, and the specific process is not described again.Representing the number of residual harmonic components in the mth subsequence of the harmonic component X in the period w,Representing the number of subintervals in the residual harmonic component,The interval amplitude probability of the q-th subinterval of the p-th residual harmonic component in the m-th subsequence of the harmonic component X in the period w is represented,Representing the p-th residual harmonic component,Representing the harmonic component of the wave,Representing the spearman correlation coefficient between the p-th residual harmonic component and the harmonic component X,For adjusting parameters, the checked value is 1, and the denominator is avoided to be 0.
When the proportion of the high amplitude in the residual harmonic component is greater,The larger the value of (2) is at the same timeThe smaller the value of (c) is, the lower the similarity between the residual harmonic component and the harmonic component is, so that the harmonic intensity coefficient isThe larger the residual harmonic component is, the more complex it is. The smaller the DTW distance between the separated residual fundamental wave component and the harmonic wave component is, the higher the similarity of the two sequences is, which means that the fundamental wave information contained in the harmonic wave component is still more, and the distortion index of the strength of the defective harmonic wave is obtainedThe larger the effect of the harmonics on the transformer circuitry is further explained.
In a distribution room, a plurality of transformers are usually arranged to run simultaneously, and the stability of a power supply system is ensured under the working state of high load. However, the current data difference between different transformers is too large, so that the problems of unbalanced load and overload of part of transformers are easy to occur. And the more the transformer is overloaded and the stronger the harmonic influence on the line, the higher the safety risk. Therefore, the harmonic intensity distortion consistency coefficient is constructed based on the analysis of the working states of different transformers.
According to the above steps, for each original current signal of the transformer, each period is calculatedValues.
For each transformer, the distortion indexes of the harmonic intensities of the incomplete phases of the transformers are arranged according to the time sequence to form a distortion index sequence of each transformer. When the distribution room works normally, the difference between the corresponding distortion index sequences of the transformers is small. And adding the amplitudes of the residual fundamental wave components and the fundamental wave components of each transformer at the corresponding sampling time to obtain each approximate fundamental wave component JS of each transformer. And a fitting curve of the approximate fundamental component JS is obtained by adopting a least square method technology and is marked as F, the least square method is a known technology, and the specific process is not repeated. The harmonic intensity distortion consistency coefficient is calculated as follows:
In the method, in the process of the invention, Representing the harmonic intensity distortion uniformity coefficient of the transformer u,The number of the power transformers in the distribution room is represented,Representing the approximate fundamental component of the transformer u,A fitted curve representing the approximate fundamental component JS,The function is a mean square error function and,AndThe sequences of skew indices for transformers u and v are shown,Representing a sequence warping index sequenceAnd a skew index sequenceBetween (a) and (b)Distance. The harmonic intensity distortion consistency coefficient obtaining flow chart is shown in fig. 2.
If it isThe smaller the value, the greater the fitting degree of the approximate fundamental component and the sine waveform, which shows that the fundamental information of the current data of the transformer u appears more obviously,The smaller the value, the smaller the harmonic intensity distortion difference between the transformer u and the transformer v, the less likely the transformer u is overloaded. The higher the obtained harmonic intensity distortion consistency coefficient is, the better the working state of the transformer u is.
And step S003, detecting the working state of the transformer by adopting a COF detection algorithm to identify the safety problem of the transformer.
So far, after the harmonic intensity distortion consistency coefficient of each transformer is calculated, the potential safety faults in the distribution room are judged by combining a COF algorithm. The input of the COF algorithm is the corresponding harmonic intensity distortion consistency coefficient of each transformer, the output is an abnormal factor of each harmonic intensity distortion consistency coefficient, and the COF algorithm is a known technology and the specific process is not repeated. If the abnormal factor value corresponding to the harmonic intensity distortion consistency coefficient of a certain transformer is larger than the preset abnormal threshold value, the preset abnormal threshold value of the embodiment is 3, which indicates that the transformation has a safety problem; on the contrary, the transformer has no safety problem.
Further, because more harmonic current appears in the current of the transformer in the distribution room, the transformer load in the distribution room can be unbalanced, the line temperature rises, and then the safety problems such as fire disaster are caused. According to the embodiment, the harmonic intensity distortion consistency coefficient corresponding to each transformer is constructed by analyzing the output current data of each transformer in the distribution room, and the safety of the distribution room is controlled according to the harmonic intensity distortion consistency coefficient, wherein the specific control mode is as follows: collecting current data of transformers in a distribution room in the operation process, decomposing fundamental wave components and harmonic components in the transformers according to the collected current data of each transformer, analyzing harmonic characteristics to obtain harmonic intensity distortion consistency coefficients corresponding to each transformer, conveying the obtained harmonic intensity distortion consistency coefficients corresponding to each transformer to a distribution room safety control system, obtaining a transformer with potential safety hazard in the distribution room according to the harmonic intensity distortion consistency coefficients corresponding to each transformer in the distribution room and a preset abnormal threshold value, closing the transformer with potential safety hazard in the distribution room in time, calling an idle transformer to replace the transformer, and checking the fault transformer in time.
Based on the same inventive concept as the above method, the embodiment of the invention also provides an intelligent distribution room safety management and control system, which comprises a memory, a processor and a computer program stored in the memory and running on the processor, wherein the processor realizes the steps of any one of the above intelligent distribution room safety management and control methods when executing the computer program.
It should be noted that: the sequence of the embodiments of the present invention is only for description, and does not represent the advantages and disadvantages of the embodiments. And the foregoing description has been directed to specific embodiments of this specification. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments.
The foregoing description of the preferred embodiments of the present invention is not intended to be limiting, but rather, any modifications, equivalents, improvements, etc. that fall within the principles of the present invention are intended to be included within the scope of the present invention.
Claims (7)
1. The intelligent distribution room safety control method is characterized by comprising the following steps of:
collecting original current signals of each transformer in each period;
Acquiring fundamental wave components and harmonic wave components of each original current signal by using a sliding window iterative DFT algorithm; acquiring harmonic distortion coefficients of each subsequence in each period according to the harmonic components and the amplitude of the original current signal; acquiring the jump disturbance degree of the harmonic component of each subsequence in each period according to the harmonic distortion coefficient and the amplitude of the harmonic component; using a decomposition technique to obtain a residual fundamental component and a residual harmonic component of the harmonic component in each sub-sequence in each period; acquiring harmonic intensity coefficients of each subsequence of the harmonic component in each period according to the amplitude of the residual harmonic component; acquiring the defective harmonic intensity distortion index of each original current signal in each period according to the relation among the harmonic component jump disorder degree, the harmonic intensity coefficient and each subsequence in the harmonic component; acquiring harmonic intensity distortion consistency coefficients of all transformers according to the defective harmonic intensity distortion index, the residual fundamental wave component and the fundamental wave component of the transformers;
Detecting the working state of the transformer according to the harmonic intensity distortion consistency coefficient and a detection algorithm;
the method for obtaining the harmonic component jump disorder degree of each subsequence in each period according to the harmonic distortion coefficient and the amplitude of the harmonic component comprises the following steps:
For each subsequence of each period, calculating the difference value between each amplitude value in the subsequence of the harmonic component and the amplitude value of the first sampling moment, calculating the average value of all amplitude values in the subsequence of the harmonic component, calculating the absolute value of the ratio of the difference value to the average value, calculating the sum value of all the absolute values in the subsequence of the harmonic component, calculating the ratio of the harmonic distortion coefficient of the subsequence to the number of sampling moments in the subsequence as a first ratio, and taking the product of the first ratio and the sum value as the jump turbulence degree of the harmonic component of each subsequence in each period;
The obtaining the distortion index of the defective harmonic intensity of each original current signal in each period according to the relation among the jump disorder degree of the harmonic component, the harmonic intensity coefficient and each subsequence in the harmonic component comprises the following steps:
For each period, calculating the sum value of the harmonic component jump disorder degree and the harmonic intensity coefficient of each subsequence as a fourth sum value, calculating the DTW distance between each subsequence in the harmonic components and the residual fundamental wave component of each subsequence, calculating the sum value of the DTW distance and a preset adjusting parameter as a fifth sum value, calculating the ratio of the fourth sum value to the fifth sum value, and taking the average value of all the ratios in each period as the distortion index of the residual harmonic intensity of each original current signal in each period;
the obtaining the harmonic intensity distortion consistency coefficient of each transformer according to the defective harmonic intensity distortion index, the residual fundamental component and the fundamental component of the transformer comprises the following steps:
The distortion indexes of the incomplete harmonic intensities of each period of each transformer are arranged according to the time sequence to form a distortion index sequence of each transformer;
Adding the residual fundamental wave components of the transformers and the amplitudes of the fundamental wave components at the corresponding sampling moments to obtain approximate fundamental wave components of the transformers, and obtaining fitting curves of the approximate fundamental wave components by using a curve fitting algorithm;
and calculating the mean square error between the approximate fundamental component of each transformer and the fitting curve of the approximate fundamental component, calculating the DTW distance between each transformer and the distortion index sequences of other transformers, calculating the sum of all the DTW distances of each transformer and the mean square error, and taking the absolute value of the reciprocal of the sum as the harmonic intensity distortion consistency coefficient of each transformer.
2. The intelligent electrical room safety control method of claim 1, wherein the obtaining fundamental and harmonic components of each raw current signal using a sliding window iterative DFT algorithm comprises:
processing the original current signals of each period by using a sliding window iterative DFT algorithm to obtain fundamental wave components of the original current signals;
And subtracting the amplitude of the original current signal at the same moment from the amplitude of the fundamental component to obtain the harmonic component of the original current signal.
3. The intelligent power distribution room safety control method as claimed in claim 1, wherein the obtaining harmonic distortion coefficients of each sub-sequence in each period according to the harmonic component and the amplitude of the original current signal comprises:
Dividing the original current signal, the fundamental wave component and the harmonic component of each period into a first preset number of subsequences with equal length respectively;
For each sub-sequence of each period, calculating the square of the amplitude of the harmonic component at each sampling time in the sub-sequence as a first square value, calculating the square of the amplitude of the original current signal at each sampling time in the sub-sequence as a second square value, calculating the sum of all the first square values in the sub-sequence as a first sum, calculating the sum of all the second square values in the sub-sequence as a second sum, and taking the square root of the ratio of the first sum to the second sum as the harmonic distortion coefficient of each sub-sequence of each period.
4. The intelligent distribution room safety control method according to claim 1, wherein the acquiring the residual fundamental wave component and the residual harmonic wave component of the harmonic wave component in each sub-sequence in each period by using a decomposition technology comprises:
for each period, taking the subsequence of the harmonic component as the input of a decomposition algorithm, and acquiring a second preset number of subband components;
and calculating correlation coefficients of each sub-band component and the sub-sequence of the harmonic components, taking the sub-band component with the largest correlation number as a residual fundamental component of the harmonic components, and taking sub-band components except the residual fundamental component as residual harmonic components.
5. The intelligent power distribution room safety control method according to claim 1, wherein the obtaining the harmonic intensity coefficients of each sub-sequence of the harmonic component in each period according to the amplitude of the residual harmonic component comprises:
for each residual harmonic component of the harmonic component, taking the maximum amplitude value and the minimum amplitude value of the residual harmonic component as an interval upper limit and an interval lower limit to construct an interval, dividing the interval into a third preset number of equal-length subintervals, and taking the sum of the probabilities of the amplitudes of all sampling moments in each subinterval in the residual harmonic component as the interval amplitude probability of each subinterval;
For each period, calculating a correlation coefficient between the residual harmonic component and the harmonic component of each sub-sequence of the harmonic component, calculating a sum value of the correlation coefficient and a preset adjusting parameter as a third sum value, calculating a ratio of an average value of interval amplitude probabilities of all sub-intervals in the residual harmonic component of each sub-sequence of the harmonic component to the third sum value, and calculating an average value of all the ratios in each sub-sequence of the harmonic component as a harmonic intensity coefficient of each sub-sequence of the harmonic component in each period.
6. The intelligent power distribution room safety control method according to claim 1, wherein the detecting the working state of the transformer according to the harmonic intensity distortion consistency coefficient in combination with the detection algorithm comprises the following steps:
Taking the harmonic intensity distortion consistency coefficient of each transformer as the input of the COF algorithm, outputting an abnormal factor of the harmonic intensity distortion consistency coefficient, and if the abnormal factor of the harmonic intensity distortion consistency coefficient of the transformer is larger than a preset abnormal threshold value, indicating that the transformer has a safety problem; on the contrary, the transformer has no safety problem.
7. An intelligent distribution room security management system comprising a memory, a processor and a computer program stored in the memory and running on the processor, wherein the processor, when executing the computer program, implements the steps of the method of any of claims 1-6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410452296.4A CN118074058B (en) | 2024-04-16 | 2024-04-16 | Intelligent distribution room safety control method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410452296.4A CN118074058B (en) | 2024-04-16 | 2024-04-16 | Intelligent distribution room safety control method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118074058A CN118074058A (en) | 2024-05-24 |
| CN118074058B true CN118074058B (en) | 2024-07-02 |
Family
ID=91097453
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410452296.4A Active CN118074058B (en) | 2024-04-16 | 2024-04-16 | Intelligent distribution room safety control method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118074058B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118316202B (en) * | 2024-06-07 | 2024-08-09 | 国网黑龙江省电力有限公司伊春供电公司 | Real-time inspection device for power system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111337745A (en) * | 2020-03-18 | 2020-06-26 | 石家庄德堃电子科技有限公司 | Multifunctional three-phase power meter containing multiple harmonic measurement |
| CN117630558A (en) * | 2024-01-26 | 2024-03-01 | 山西鼎泰源科技有限公司 | Intelligent fault diagnosis system for electromechanical equipment |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3783639B2 (en) * | 2002-03-20 | 2006-06-07 | 株式会社デンソー | AC generator motor for vehicle |
| JP5850709B2 (en) * | 2011-11-07 | 2016-02-03 | 株式会社ダイヘン | Single operation detection device for grid-connected inverter device |
| CN117368799B (en) * | 2023-12-07 | 2024-02-23 | 山西思极科技有限公司 | Diagnosis method for short-circuit fault of power transmission line of power system |
-
2024
- 2024-04-16 CN CN202410452296.4A patent/CN118074058B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111337745A (en) * | 2020-03-18 | 2020-06-26 | 石家庄德堃电子科技有限公司 | Multifunctional three-phase power meter containing multiple harmonic measurement |
| CN117630558A (en) * | 2024-01-26 | 2024-03-01 | 山西鼎泰源科技有限公司 | Intelligent fault diagnosis system for electromechanical equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118074058A (en) | 2024-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN118074058B (en) | Intelligent distribution room safety control method and system | |
| CN105067966B (en) | The low-voltage alternating-current fault arc detection method of feature based modal components energy spectrometer | |
| Yılmaz et al. | An improved automated PQD classification method for distributed generators with hybrid SVM-based approach using un-decimated wavelet transform | |
| CN110068759B (en) | Fault type obtaining method and device | |
| Guo et al. | Detection and line selection of series arc fault in multi-load circuit | |
| CN117368799B (en) | Diagnosis method for short-circuit fault of power transmission line of power system | |
| CN109581270A (en) | Method based on multi-scale wavelet transform Fast Identification CVT failure | |
| CN117390408B (en) | Power transformer operation fault detection method and system | |
| CN110320436A (en) | Flexible direct current power distribution network high resistance earthing fault detection method based on color relationship classifier | |
| CN112183400A (en) | Novel latent fault feature extraction method and system for distribution transformer | |
| Chen et al. | Reliable detection method of variable series arc fault in building integrated photovoltaic systems based on nonstationary time series analysis | |
| CN115600754A (en) | Method, system, equipment and storage medium for predicting potential safety hazard of user power utilization | |
| CN113917294B (en) | Intelligent self-adaptive arc detection method based on wavelet decomposition and application device thereof | |
| CN117540235B (en) | Real-time monitoring method for temporary electricity utilization abnormal data | |
| CN113484658A (en) | Method, system, medium, and electronic device for diagnosing arc fault | |
| CN112748318A (en) | Series fault arc detection method based on improved PSO-BP neural network | |
| CN117575247A (en) | Province-station scheduling management method and system | |
| CN117074783A (en) | Real-time monitoring and early warning method for overheat state of power equipment | |
| CN110794210B (en) | Method and device for judging voltage harmonic isolation effect, power supply system, computer equipment and storage medium | |
| CN110017894B (en) | Method and device for filtering random noise in vibration and sound detection of transformer in running state | |
| Florkowski | Exploitation stresses and challenges in diagnostics of electrical industrial equipment | |
| CN112329628A (en) | Method and device for extracting multi-feature fusion features of converter transformer vibration signal | |
| CN113128396B (en) | Power quality composite disturbance classification method | |
| CN117767228B (en) | Bus protection method, device, equipment and medium based on transient frequency spectrum components | |
| CN118091526B (en) | On-line monitoring method and device for measuring error of current transformer |
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 |