CN113484646A - Transformer substation secondary circuit complex fault positioning device and diagnosis method - Google Patents

Transformer substation secondary circuit complex fault positioning device and diagnosis method Download PDF

Info

Publication number
CN113484646A
CN113484646A CN202110823899.7A CN202110823899A CN113484646A CN 113484646 A CN113484646 A CN 113484646A CN 202110823899 A CN202110823899 A CN 202110823899A CN 113484646 A CN113484646 A CN 113484646A
Authority
CN
China
Prior art keywords
event
variable
data
sample
events
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.)
Pending
Application number
CN202110823899.7A
Other languages
Chinese (zh)
Inventor
张惠山
朱卫民
渠红涛
孟荣
王会增
刘晓飞
池威威
王昭雷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Maintenance Branch of State Grid Hebei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Maintenance Branch of State Grid Hebei Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Maintenance Branch of State Grid Hebei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202110823899.7A priority Critical patent/CN113484646A/en
Publication of CN113484646A publication Critical patent/CN113484646A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)

Abstract

The invention discloses a transformer substation secondary circuit complex fault positioning device and a diagnosis method, wherein the fault positioning device comprises a main control CPU, a measurement unit module and a data communication preposed unit; the main control CPU communicates data information with the data communication preposition unit; the function of the measurement unit module is the collection function of various information quantities of events, the function of the data communication preposition unit is to complete data communication with each measurement unit, and is responsible for temporary storage of data information, receiving the command of the device host control system and issuing the command to each measurement unit, completing time synchronization and clock correction of each measurement unit, and sending data to the host control system. And the fault positioning device acquires data information and performs logical diagnosis reasoning on different direct current determined data information acquisition points. The diagnosis method adopts a mathematical algorithm of pattern matching to carry out analysis, diagnosis and reasoning, and provides comprehensive data information for accurate reasoning and positioning of faults, thereby ensuring the accuracy and effectiveness of the faults.

Description

Transformer substation secondary circuit complex fault positioning device and diagnosis method
Technical Field
The invention belongs to the technical field of artificial intelligence, and particularly relates to a transformer substation secondary circuit complex fault positioning device and a diagnosis method.
Background
Some secondary circuits of the current transformer substation have complex formation mechanisms (typical complex faults include grounding of a direct current system of the transformer substation and series faults of an alternating current and direct current circuit), and the obvious characteristic of the complex formation mechanisms is uncertainty of individual characteristics presented by the faults, which is different from other faults. Besides some fault inherent characteristics which are clear, the method also has more individual unique characteristics which present uncertainty. The concrete aspects are as follows:
for a direct current ground fault, different fault points present more different fault characteristic information besides causing the voltage to ground of the substation direct current system to offset. For example, the dc ground faults of different branches have different ground resistances and different offset degrees of the dc bus voltage. And when the ground fault occurs to different branches, the individual characteristics of the branches show uncertainty. Namely, the grounding resistance and the direct current offset voltage of the grounding fault of different branches show the same fault characteristics. Meanwhile, the direct current grounding characteristics of the same branch are different under different events and different environmental conditions. For the direct-current grounding fault, the existing direct-current grounding branch line selection device carries out branch line confirmation according to the current balance principle of a direct-current branch line, the principle is simple and old, the detection sensitivity is low, the wrong selection, wrong selection and missing selection rates of the grounding fault branch line are high, the specific fault position of the grounding branch line cannot be determined, and therefore the auxiliary support effect of effective fault processing cannot be really achieved.
Secondly, for the series fault of the alternating current and direct current system, the individuality of the frequent fault characteristics is strong. Some faults present the characteristic phenomena of grounding and instantaneous recovery of a direct current system. Some of the monitoring systems present a large amount of irrelevant alarm information frequently sent by mistake by the instant monitoring system, fault characteristics of instant resetting and the like. At present, the fault can be judged only according to experience, and a reliable auxiliary supporting tool is not available.
Disclosure of Invention
The invention aims to solve the technical problem of providing a transformer substation secondary circuit complex fault positioning device and a diagnosis method, which adopt a mode matching mathematical algorithm to carry out analysis diagnosis and reasoning and provide comprehensive data information for accurate reasoning and positioning of faults so as to ensure the accuracy and effectiveness of the faults.
In order to solve the problems, the technical scheme adopted by the invention is as follows:
a transformer substation secondary circuit complex fault positioning device comprises a main control CPU, a measurement unit module and a data communication preposition unit;
the main control CPU communicates data information with the data communication preposition unit;
the function of the measurement unit module is the collection function of various information quantities of events, the function of the data communication preposition unit is to complete data communication with each measurement unit, and is responsible for temporary storage of data information, receiving the command of the device host control system and issuing the command to each measurement unit, completing time synchronization and clock correction of each measurement unit, and sending data to the host control system.
As a further improvement of the positioning device, the measuring unit module comprises a clamp type current measuring component, a contact type voltage and current measuring component and a temperature/humidity non-electrical quantity measuring component;
the fault occurrence of a secondary circuit is regarded as an event occurrence, and the change of the electrical quantity information and the change of the non-electrical quantity information which are accompanied by the fault occurrence of the secondary circuit of the measured detection point form the characteristic variable of the event;
the clamp type current measuring assembly, the contact type voltage and current measuring assembly and the temperature/humidity non-electrical quantity measuring assembly finish the real-time collection of the change of the characteristic variable, record data information and time information for generating the data, and send the data to the device main control system.
As a further improvement of the positioning device, the main control CPU performs information interaction with the transformer substation monitoring system, receives command detection information sent by the transformer substation monitoring system, simultaneously sends a command reasoning detection result to the transformer substation monitoring system for transformer substation operation and maintenance personnel to review, and receives a fault data updating command sent by the monitoring system to update the database of the monitoring system.
As a further improvement of the positioning device, the fault data refers to a certain fault event which is not identified by manual work, and after receiving the instruction of the monitoring system, the main control CPU records and accumulates the data measurement information corresponding to the fault event data to form a new sample of the type of event.
As a further improvement of the positioning device, the main control CPU receives a clock message issued by the monitoring system, and completes the clock correction of the host, the data communication front-end unit and the measurement unit module.
A diagnosis method using a transformer substation secondary circuit complex fault positioning device firstly determines the fact that: if the complex fault occurring at a certain determined position of a certain direct current branch is expressed as an event, the occurring event is called a sample; for a plurality of samples of the event, some data information presents a physical characteristic of normal distribution, namely information which is completely consistent with the data information cannot be achieved by means of a group of determined data, but because the data information presents a certain distribution characteristic on a macroscopic level, the data which conforms to the characteristic has the possibility of being a sample of the event to a certain extent, and the higher the matching degree is, the higher the probability of conforming to the event is; the complex fault of the secondary circuit is like the moving track of electrons in atoms, under the macroscopic condition, the moving track of the secondary circuit is represented as different electron clouds, and the probability distribution characteristics of the moving track of the secondary circuit are reflected by the shapes of the different electron clouds.
As a further improvement of the diagnosis method, based on the fact, the logical inference method of pattern matching of similarity recognition is as follows: for some event P, it contains the following feature vectors:
P=[p′1,p′2,......,p′n];
each vector contains the following variables:
p′i=(fi1,fi2,......,fiq);
then:
Figure BDA0003172927560000031
distribution characteristics of independent event random variables:
according to the statistical theory, for an event P, the events which occur are independent of each other, the variable values of the events are subject to the characteristic rule of normal distribution, and for n event sequence sets of the independent events P of the same type: { P1,P2,......,PnEach event PiVariable f in (1)i,jWherein { i belongs to (0-n), j belongs to (0-q) } obeys the characteristic of normal distribution; i.e. for the variable fi,jThe values obey the distribution characteristics of a normal probability density function:
Figure BDA0003172927560000032
wherein the variable x means fi,jμ is the mean and σ is the variance;
for a set of events P: { P1,P2,......,PnAny one of independent variables f ini,jWill become a sample space
Figure BDA0003172927560000033
The mean μ and variance σ of the variable can be calculated;
Figure BDA0003172927560000034
thus, for each variable f of the event PijAll will obtain an average value
Figure BDA0003172927560000035
A standard pattern of events P is thus formed:
Figure BDA0003172927560000036
Figure BDA0003172927560000041
as a further improvement of the diagnostic method, a weighted model of the event features: some of the variables of the event P are changed, and have a larger value; but the change of the variable values is small and stable; the same event, the larger change of the value and the higher distribution probability indicate that the value is higher in activity degree; similarly, a smaller change indicates a lower activity level for the value;
corresponding to the event P, the mean square error matrix of the sample is:
Figure BDA0003172927560000042
the mean square error of the sample reflects the degree of deviation of each variable of the event from the standard mean value of the variable;
defining a variable weight coefficient S:
s (x) σ · f (x), where σ is the sample mean square error of the variable, and f (x) is its probability distribution density function;
the probability distribution density function f (x) reflects the magnitude of probability value corresponding to the value of a variable with normal distribution characteristics under specific mean and variance, and the physical meaning of the product of the two indicates that: for any sample, the function s (x) reflects a quantitative index of the matching degree of the variables related to the specific event after the values of the variables refer to the probability distribution characteristics and the standard pattern under the condition that the standard pattern is taken as a reference.
As a further improvement of the diagnostic method, the pattern matching method is as follows:
for a specific event P' which occurs, judging whether the event P and the event P belong to the same type of event, and adopting the following modes:
s (X) ═ P' (X) · P (σ), that is:
Figure BDA0003172927560000043
using S (X) and standard PbzThe absolute difference of the corresponding variable is reflected by the difference degree between the standard value and the corresponding variable, and the matrix M is used forcz(X) is as follows:
Figure BDA0003172927560000051
in an ideal, perfect match situation, if a particular event P' is exactly one sample of event P, then MczThe value of each variable in (X) should be 0, however, in practice there must be a deviation, there is a certain amount of deviation, and it cannot be denied that event P' is not a sample of event P.
It is also necessary to consider the other case, namely the case where said event is for a specific event tkThe temporal events, but the real objective physical events, an event is generated and changed with the time, and different times correspond to different events P, i.e. a set of event sequences is formed:
Pt jh={Pt1,Pt2,...,Ptk....,Ptmthe pattern matching of events should not be limited to matching of variable values at specific times only, but also the development of events, i.e. the degree of matching at different times,
Pt jhwill generate different corresponding different time
Figure BDA0003172927560000052
Thereby forming a set of sequences of difference values:
Figure BDA0003172927560000053
set Mcz|T(X) each element in (X) is formed by the development of the same event over time; the event P development is continuous, so that the value of the event P development is continuously related; similarly, there is a correlation between the values of the events P 'at different times, and if the event P' is a sample event formed by the events P, the set M formed is in particularcz|TAnd (X) the corresponding variables contained in each element have linear correlation, if the correlation degree is higher, the P' is determined to be a sample of the event P and the event P is the same as the event P, otherwise, the event P cannot be considered to belong to the same event.
As a further improvement of the diagnostic method, the concept of correlation coefficient of statistical theory is utilized to detect the degree of correlation between the variables, and the analysis method is as follows:
for a corresponding set of difference values M with M time instantscz|T(X) two successive instants t of aggregation of the elementsiAnd ti+1Is an element of
Figure BDA0003172927560000054
And if the sample data size of the whole set is m, the correlation coefficient rho calculation formula is as follows:
Figure BDA0003172927560000061
whereby an independent variable f for an event PijWill generate an m correlation coefficients piThe data series of (2) is subjected to summation operation to obtain a data variable of a correlation coefficient:
Figure BDA0003172927560000062
it is obvious that
Figure BDA0003172927560000063
Each variable f corresponding to an event PijEventually, a matrix of correlation coefficients is formed:
Figure BDA0003172927560000064
the data correlation of the outstanding variable, its correlation coefficient is close to 1, and less than 0.5, consider its degree of correlation to be lower;
when close to 0, it can be considered that there is no correlation; therefore, for data with a lower degree of similarity, it is weakened; taking the above factors into consideration, for the set ρΣPerforming mathematical operation, each variable operation formula
Figure BDA0003172927560000065
Calculating to obtain a new correlation coefficient matrix rho'ΣIt is a q × n matrix;
for matrix ρ'ΣThe values of the elements in (a) are summed:
Figure BDA0003172927560000066
the condition that the correlation coefficient is greater than 0.6 can be considered to have correlation is changed into 0 quantity 6 after mathematical square operation, so when the correlation coefficient is greater than 0.6
Figure BDA0003172927560000067
When the condition is satisfied, it can be inferred that the event P' is a sample event of the event P.
Adopt the produced beneficial effect of above-mentioned technical scheme to lie in:
the invention provides a method and a device for diagnosing and reasoning the complex fault. And determining data information acquisition points for different direct currents. The collection point information contains electric quantity information such as current, voltage and the like, and also comprises non-electric quantity information such as temperature, humidity and the like. The mathematical algorithm matched with the mode is adopted for analysis, diagnosis and reasoning, and comprehensive data information is provided for accurate reasoning and positioning of the fault, so that the accuracy and the effectiveness of the fault are ensured.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed for the embodiments or the prior art descriptions will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
FIG. 1 is a schematic view of a measurement cell module;
FIG. 2 is a schematic diagram of a data communication head unit of the apparatus;
FIG. 3 is a schematic diagram of the device host hardware mechanism;
FIG. 4 is a host CPU control system software flow diagram;
fig. 5 is a variable diagram of a branch circuit.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the application, its application, or uses. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments according to the present application. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
The relative arrangement of the components and steps, the numerical expressions, and numerical values set forth in these embodiments do not limit the scope of the present application unless specifically stated otherwise. Meanwhile, it should be understood that the sizes of the respective portions shown in the drawings are not drawn in an actual proportional relationship for the convenience of description. Techniques, methods, and apparatus known to those of ordinary skill in the relevant art may not be discussed in detail but are intended to be part of the specification where appropriate. In all examples shown and discussed herein, any particular value should be construed as merely illustrative, and not limiting.
As shown in the figures 1-3 of the drawings,
a transformer substation secondary circuit complex fault positioning device comprises a main control CPU, a measurement unit module and a data communication preposition unit;
the main function of the measuring unit module is the acquisition function of various information quantities of events. A failure of a secondary loop may be considered an event. The occurrence of this event is accompanied by a change in the electrical quantity information (including current and voltage information) and a change in the non-electrical quantity information (including a change in the quantity of information such as temperature and humidity) at the monitoring point set manually. These quantities of information will be part of the characteristic variables of the event. The measuring unit module collects the change of the information quantity in real time, records data information, time information for generating the data and other data, and generates the data to the device main control system. The data communication preposition unit mainly completes data communication with each measurement unit, is responsible for temporarily storing data information, receives a command of a device host control system and issues the command to each measurement unit, completes time synchronization and clock correction of each measurement unit, sends data to the host control system and the like.
The device host system communicates data information with the data communication front-end unit. Data information and command information are transmitted. Meanwhile, the host CPU system and the transformer substation monitoring system perform information interaction, mainly receive command detection information sent by the transformer substation monitoring system, and simultaneously send a command reasoning detection result to the transformer substation monitoring system for transformer substation operation and maintenance personnel to review. And receiving a fault data updating command issued to the monitoring system to update the database of the monitoring system. The fault data is a certain fault event which is not manually identified, and after the host receives the instruction of the monitoring system, the data measurement information record corresponding to the fault event data is accumulated to become a new sample of the type of event. The device host provides a human-computer interaction interface to complete the functions of program debugging, information retrieval and the like. The host system receives the clock message issued by the monitoring system, and completes the clock correction of the host, the data preposition unit and the measurement unit.
The software flow of the host CPU control system is shown in FIG. 4, and the software system comprises a main program and an interrupt program. The main program completes clock correction, receiving of command information of the monitoring system and data interaction. The interrupt program mainly completes the diagnosis and reasoning function of the fault.
The invention provides a method and a device for diagnosing and reasoning the complex fault. Firstly, data information acquisition points are determined for different direct currents. The collection point information contains electric quantity information such as current, voltage and the like, and also comprises non-electric quantity information such as temperature, humidity and the like. The accurate reasoning and positioning of the fault provides comprehensive data information, thereby ensuring the accuracy and effectiveness of the fault. The main idea of reasoning is to use a mathematical algorithm of pattern matching to perform analysis, diagnosis and reasoning.
The pattern matching algorithm is based on the individuality that the data information generated by a fault event reflects the fault characteristics of the fault event belonging to the fault event. Individuality here refers to the characteristics of affiliation to the present fault event. The system refers to the fault characteristics of a certain specific position of a certain branch circuit.
The patent firstly determines the fact that the secondary loop complex fault is: if a complex fault occurring at a certain determined position of a certain direct current branch is expressed as an event, the occurring event is called a sample. And the value of the related data information caused by a plurality of such events is not completely determined. That is, an event is matched against a value having completely certain information, and the event of the type cannot be matched. However, for multiple samples of the event, some data information exhibits a normally distributed physical property. That is, information completely consistent with a certain set of data cannot be obtained by means of the certain set of data, but because the data information of the certain set of data shows a certain distribution characteristic on a macroscopic level, data conforming to the characteristic has a possibility of being a sample of the event to some extent, and the higher the matching degree is, the higher the probability of conforming to the event is. The secondary loop complex fault is like the trajectory of an electron in an atom. For a single electron, even if its motion parameters are known, it is not possible to determine where it will be at a particular time interval thereafter. But under macroscopic conditions, the trajectory of its motion appears as a distinct cloud of electrons. The different electron cloud shapes reflect the probability distribution characteristics of their motion trajectories.
Based on the facts, the main idea of pattern matching is to determine the related parameter information of the normal distribution characteristic of a certain event by using the determined sample data as accumulation after the certain event has occurred by means of a certain mathematical reasoning algorithm. For a currently occurring event, matching already existing events according to a mathematical reasoning algorithm, thereby deriving a reasoning for the event.
The similarity recognition pattern matching logic reasoning method comprises the following steps:
for some event P, it contains the following feature vectors:
P=[p′1,p′2,......,p′n];
each vector contains the following variables:
p′i=(fi1,fi2,......,fiq);
then:
Figure BDA0003172927560000091
distribution characteristics of independent event random variables:
according to the statistical theory, the probability of the occurrence of mutually independent random events in the natural physical world follows the regular characteristic of normal distribution. Based on this theory, for event P, the events that occur are independent of each other. The variable values contained in the normal distribution model obey the characteristic rule of normal distribution. Event sequence sets for n independent events P of the same type are noted: { P1,P2,......,PnEach event PiVariable f in (1)i,jWhere { i ∈ (0. about.n), j ∈ (0. about.q) } should obey the characteristics of a normal distribution. I.e. for the variable fi,jThe values obey the distribution characteristics of a normal probability density function:
Figure BDA0003172927560000101
wherein the variable x means fi,jμ is the mean and σ is the variance.
For a set of events P: { P1,P2,......,PnAny one of independent variables f ini,jWill become a sample space
Figure BDA0003172927560000102
The mean μ and variance σ of the variables can be calculated.
Figure BDA0003172927560000103
Thus, for each variable f of the event PijAll will obtain an average value
Figure BDA0003172927560000104
A standard pattern of events P is thus formed:
Figure BDA0003172927560000105
weight model of event features: of the variables of the event P, some variables are changed, and there is a change in a larger number. But the change in the value of some variables is small and stable. The same event, with a larger value change and a higher probability of distribution, indicates that the value is more active. Similarly, a smaller change indicates a lower activity level for the value.
Corresponding to the event P, the mean square error matrix of the sample is:
Figure BDA0003172927560000106
description of the drawings: the mean square error of the sample reflects the degree of deviation of each variable of the event from its standard mean.
Defining a variable weight coefficient S:
s (x) σ · f (x), where σ is the sample mean square error of the variable and f (x) is its probability distribution density function.
Description of the drawings: the probability distribution density function f (x) reflects the magnitude of probability value corresponding to the value of a variable with normal distribution characteristics under specific mean and variance, and the physical meaning of the product of the two indicates that: for any sample, the function s (x) reflects a quantitative index of the matching degree of the variables related to the specific event after the values of the variables refer to the probability distribution characteristics and the standard pattern under the condition that the standard pattern is taken as a reference.
Brief description of the pattern matching method:
for a specific event P' which occurs, judging whether the event P and the event P belong to the same type of event, and adopting the following modes:
s (X) ═ P' (X) · P (σ), that is:
Figure BDA0003172927560000111
using S (X) and standard PbzThe absolute difference of the corresponding variable is reflected by the difference degree between the standard value and the corresponding variable, and the matrix M is used forcz(X) is as follows:
Figure BDA0003172927560000112
in an ideal, perfect match situation, if a particular event P' is exactly one sample of event P, then MczThe value of each variable in (X) should be 0. However, the actual situation is certainly biased. There is a certain amount of deviation that cannot be negated that event P' is not a sample of event P.
It is also necessary to consider the other case, namely the case where said event is for a specific event tkThe event of time. However, in real objective physical events, the occurrence of an event changes with the time, and different times correspond to different events P, i.e. a set of event sequences is formed:
Pt jh={Pt1,Pt2,...,Ptk....,Ptm}. The pattern matching of events should not be limited to matching of variable values at a particular time, but should also take into account the development of the event, i.e. the degree of matching at different times.
Pt jhWill generate different corresponding different time
Figure BDA0003172927560000113
Thereby forming a set of sequences of difference values:
Figure BDA0003172927560000114
obviously, set Mcz|TEach element in (X) is formed by the development of the same event over time. The event P development is continuous, and therefore, the value is continuously correlated. Similarly, there is a correlation between the values of the event P' at different times. If the event P' is a sample event formed by the event P, the set M formed in particularcz|TThere must be a linear correlation between the corresponding variables contained in each element of (X). If the degree of correlation is high, then,it can be determined that P' is a sample of event P and that event P is homogeneous. Otherwise, it cannot be considered as belonging to the same event.
The degree of correlation between its variables can be detected using the concept of correlation coefficients of statistical theory. The analysis method is as follows: for a corresponding set of difference values M with M time instantscz|T(X) two successive instants t of aggregation of the elementsiAnd ti+1Is an element of
Figure BDA0003172927560000121
And
Figure BDA0003172927560000122
and if the sample data size of the whole set is m, the correlation coefficient rho calculation formula is as follows:
Figure BDA0003172927560000123
whereby an independent variable f for an event PijWill generate an m correlation coefficients piThe data series of (1). And carrying out summation operation on the data to obtain a data variable of a correlation coefficient:
Figure BDA0003172927560000124
it is obvious that
Figure BDA0003172927560000125
Each variable f corresponding to an event PijEventually, a matrix of correlation coefficients is formed:
Figure BDA0003172927560000126
the degree of correlation is considered to be low when the correlation coefficient of a variable whose data correlation is prominent is close to 1 and less than 0.5. When close to 0, it can be considered that there is no correlation. Thus, for numbers with a lower degree of similarityAccordingly, it can be weakened. Taking the above factors into consideration, for the set ρΣPerforming mathematical operation, each variable operation formula
Figure BDA0003172927560000127
Calculating to obtain a new correlation coefficient matrix rho'ΣIt is a q × n matrix.
For matrix ρ'ΣThe values of the elements in (a) are summed:
Figure BDA0003172927560000131
the condition that the correlation coefficient is greater than 0.6 can be considered to have correlation is changed into 0 quantity 6 after mathematical square operation, so when the correlation coefficient is greater than 0.6
Figure BDA0003172927560000132
When the condition is satisfied, it can be inferred that the event P' is a sample event of the event P.
The foregoing reasoning process for dc ground fault at one site of only one branch of the event P. In fact, there are various branches in the substation, including control, protection, remote signaling, wave recording, etc. For a ground fault occurring at a certain position of a first group of control loops of a certain 220kV circuit breaker, the formed sample can be transplanted to a ground fault occurring at a certain position of a first group of control loops of another 220kV circuit breaker, thereby forming a sample of a ground fault occurring at a certain position of another branch circuit. Thus, event P is given an attribute M, M ═ M that characterizes its legs1,m2,......,m3]. Other events with the same characteristic variables can be migrated to each other. Of course, the event P contains a variable that is the variable to which the rate belongs. After the event P is transplanted, the variables are also associated correspondingly and transformed into the variables of the transplanted event, which are called mapping of the event variables. Mapping description as shown in fig. 5, mapping transformation exists between the relevant variables to which the event P and the event Q belong:
Figure BDA0003172927560000133
the analysis is an intelligent reasoning method for a reasoning process of pattern matching of a secondary circuit complex secondary circuit fault event.

Claims (10)

1. The utility model provides a complicated fault positioning device of transformer substation's secondary circuit which characterized in that: the system comprises a main control CPU, a measuring unit module and a data communication preposition unit;
the main control CPU communicates data information with the data communication preposition unit;
the function of the measurement unit module is the collection function of various information quantities of events, the function of the data communication preposition unit is to complete data communication with each measurement unit, and is responsible for temporary storage of data information, receiving the command of the device host control system and issuing the command to each measurement unit, completing time synchronization and clock correction of each measurement unit, and sending data to the host control system.
2. The complex fault location device of the substation secondary circuit of claim 1, characterized in that: the measuring unit module comprises a clamp type current measuring assembly, a contact type voltage and current measuring assembly and a temperature/humidity non-electrical quantity measuring assembly; the fault occurrence of a secondary circuit is regarded as an event occurrence, and the change of the electrical quantity information and the change of the non-electrical quantity information which are accompanied by the fault occurrence of the secondary circuit of the measured detection point form the characteristic variable of the event;
the clamp type current measuring assembly, the contact type voltage and current measuring assembly and the temperature/humidity non-electrical quantity measuring assembly finish the real-time collection of the change of the characteristic variable, record data information and time information for generating the data, and send the data to the device main control system.
3. The complex fault location device of the substation secondary circuit according to claim 2, characterized in that: the main control CPU performs information interaction with the transformer substation monitoring system, receives command detection information sent by the transformer substation monitoring system, simultaneously sends a command reasoning detection result to the transformer substation monitoring system for transformer substation operation and maintenance personnel to review, and receives a fault data updating command sent by the monitoring system to update a database of the main control CPU.
4. The complex fault location device of the substation secondary circuit of claim 3, characterized in that: the fault data is a certain fault event which is not manually identified, and after the main control CPU receives the instruction of the monitoring system, the data measurement information corresponding to the fault event data is recorded and accumulated to form a new sample of the type of event.
5. The complex fault location device of the substation secondary circuit of claim 4, characterized in that: and the main control CPU receives a clock message issued by the monitoring system, and completes the clock correction of the host, the data communication preposition unit and the measurement unit module.
6. A diagnosis method using a transformer substation secondary circuit complex fault positioning device is characterized in that: such a fact is first determined for a complex fault of the secondary loop: if the complex fault occurring at a certain determined position of a certain direct current branch is expressed as an event, the occurring event is called a sample; for a plurality of samples of the event, some data information presents a physical characteristic of normal distribution, namely information which is completely consistent with the data information cannot be achieved by means of a group of determined data, but because the data information presents a certain distribution characteristic on a macroscopic level, the data which conforms to the characteristic has the possibility of being a sample of the event to a certain extent, and the higher the matching degree is, the higher the probability of conforming to the event is;
the complex fault of the secondary circuit is like the moving track of electrons in atoms, the moving characteristic of a single electron at a certain moment under the microscopic condition cannot be determined, the moving track of the single electron is represented as different electron clouds under the macroscopic condition, and the probability distribution characteristics of the moving track of the single electron cloud are reflected by the shapes of the different electron clouds.
7. The diagnosis method using the complex fault location device of the substation secondary circuit according to claim 6, characterized in that: based on the facts, the logical inference method of pattern matching of similarity recognition is as follows: for some event P, it contains the following feature vectors:
P=[p′1,p′2,......,p′n];
each vector contains the following variables:
p′i=(fi1,fi2,......,fiq);
then:
Figure FDA0003172927550000021
distribution characteristics of independent event random variables:
according to the statistical theory, for an event P, the events which occur are independent of each other, the variable values of the events are subject to the characteristic rule of normal distribution, and for n event sequence sets of the independent events P of the same type: { P1,P2,......,PnEach event PiVariable f in (1)i,jWherein { i belongs to (0-n), j belongs to (0-q) } obeys the characteristic of normal distribution; i.e. for the variable fi,jThe values obey the distribution characteristics of a normal probability density function:
Figure FDA0003172927550000022
wherein the variable x means fi,jμ is the mean and σ is the variance;
for a set of events P: { P1,P2,......,PnAny one of independent variables f ini,jWill become a sample space
Figure FDA0003172927550000023
The mean μ and variance σ of the variable can be calculated;
Figure FDA0003172927550000024
thus, for each variable f of the event PijAll will obtain an average value
Figure FDA0003172927550000025
A standard pattern of events P is thus formed:
Figure FDA0003172927550000026
8. the diagnosis method using the complex fault location device of the substation secondary circuit according to claim 7, characterized in that: weight model of event features: some of the variables of the event P are changed, and have a larger value; but the change of the variable values is small and stable; the same event, the larger change of the value and the higher distribution probability indicate that the value is higher in activity degree; similarly, a smaller change indicates a lower activity level for the value;
corresponding to the event P, the mean square error matrix of the sample is:
Figure FDA0003172927550000031
the mean square error of the sample reflects the degree of deviation of each variable of the event from the standard mean value of the variable;
defining a variable weight coefficient S:
s (x) σ · f (x), where σ is the sample mean square error of the variable, and f (x) is its probability distribution density function;
the probability distribution density function f (x) reflects the magnitude of probability value corresponding to the value of a variable with normal distribution characteristics under specific mean and variance, and the physical meaning of the product of the two indicates that: for any sample, the function s (x) reflects a quantitative index of the matching degree of the variables related to the specific event after the values of the variables refer to the probability distribution characteristics and the standard pattern under the condition that the standard pattern is taken as a reference.
9. The diagnosis method using the complex fault location device of the substation secondary circuit according to claim 8, characterized in that: the pattern matching method is as follows:
for a specific event P' which occurs, judging whether the event P and the event P belong to the same type of event, and adopting the following modes: s (X) ═ P' (X) · P (σ), that is:
Figure FDA0003172927550000032
using S (X) and standard PbzThe absolute difference of the corresponding variable is reflected by the difference degree between the standard value and the corresponding variable, and the matrix M is used forcz(X) is as follows:
Figure FDA0003172927550000033
in an ideal, perfect match situation, if a particular event P' is exactly one sample of event P, then Mcz(X) the value of each variable should be 0, however, in practice there must be a deviation, there is a certain amount of deviation, it cannot be denied that event P' is not a sample of event P;
it is also necessary to consider the other case, namely the case where said event is for a specific event tkThe temporal events, but the real objective physical events, an event is generated and changed with the time, and different times correspond to different events P, i.e. a set of event sequences is formed:
Pt jh={Pt1,Pt2,...,Ptk....,Ptmthe pattern matching of events should not be limited to matching of variable values at specific times only, but also the development of events, i.e. the degree of matching at different times,
Pt jhwill generate different corresponding different time
Figure FDA0003172927550000041
Thereby forming a set of sequences of difference values:
Figure FDA0003172927550000042
set Mcz|T(X) each element in (X) is formed by the development of the same event over time; the event P development is continuous, so that the value of the event P development is continuously related; similarly, there is a correlation between the values of the events P 'at different times, and if the event P' is a sample event formed by the events P, the set M formed is in particularcz|TAnd (X) the corresponding variables contained in each element have linear correlation, if the correlation degree is higher, the P' is determined to be a sample of the event P and the event P is the same as the event P, otherwise, the event P cannot be considered to belong to the same event.
10. The diagnosis method using the complex fault location device of the substation secondary circuit according to claim 9, characterized in that: the concept of a correlation coefficient of a statistical theory is utilized to detect the degree of correlation between variables, and the analysis method is as follows:
for a corresponding set of difference values M with M time instantscz|T(X) two successive instants t of aggregation of the elementsiAnd ti+1Is an element of
Figure FDA0003172927550000043
And
Figure FDA0003172927550000044
and if the sample data size of the whole set is m, the correlation coefficient rho calculation formula is as follows:
Figure FDA0003172927550000045
whereby an independent variable f for an event PijWill generate an m correlation coefficients piThe data series of (2) is subjected to summation operation to obtain a data variable of a correlation coefficient:
Figure FDA0003172927550000046
it is obvious that
Figure FDA0003172927550000047
Each variable f corresponding to an event PijEventually, a matrix of correlation coefficients is formed:
Figure FDA0003172927550000051
the data correlation of the outstanding variable, its correlation coefficient is close to 1, and less than 0.5, consider its degree of correlation to be lower; when close to 0, it can be considered that there is no correlation; therefore, for data with a lower degree of similarity, it is weakened; taking the above factors into consideration, for the set ρPerforming mathematical operation, each variable operation formula
Figure FDA0003172927550000052
Calculating to obtain a new correlation coefficient matrix rho'It is a q × n matrix;
for matrix ρ'The values of the elements in (a) are summed:
Figure FDA0003172927550000053
the condition that the correlation coefficient is greater than 0.6 can be considered to have correlation is changed into 0 quantity 6 after mathematical square operation, so when the correlation coefficient is greater than 0.6
Figure FDA0003172927550000054
When the condition is satisfied, it can be inferred that the event P' is a sample event of the event P.
CN202110823899.7A 2021-07-21 2021-07-21 Transformer substation secondary circuit complex fault positioning device and diagnosis method Pending CN113484646A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110823899.7A CN113484646A (en) 2021-07-21 2021-07-21 Transformer substation secondary circuit complex fault positioning device and diagnosis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110823899.7A CN113484646A (en) 2021-07-21 2021-07-21 Transformer substation secondary circuit complex fault positioning device and diagnosis method

Publications (1)

Publication Number Publication Date
CN113484646A true CN113484646A (en) 2021-10-08

Family

ID=77942761

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110823899.7A Pending CN113484646A (en) 2021-07-21 2021-07-21 Transformer substation secondary circuit complex fault positioning device and diagnosis method

Country Status (1)

Country Link
CN (1) CN113484646A (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102749549A (en) * 2012-06-21 2012-10-24 云南电力试验研究院(集团)有限公司电力研究院 Intelligent checking system for substation alternating current voltage secondary circuit
CN104699077A (en) * 2015-02-12 2015-06-10 浙江大学 Nested iterative fisher discriminant analysis-based fault diagnosis isolation method
CN105425768A (en) * 2015-11-06 2016-03-23 国网山东莒县供电公司 Electric power secondary equipment monitoring device and method
CN107958292A (en) * 2017-10-19 2018-04-24 山东科技大学 Transformer based on cost sensitive learning obscures careful reasoning method for diagnosing faults
CN108418304A (en) * 2018-03-09 2018-08-17 国家电网公司 Substation secondary circuit state monitoring method, apparatus and system
CN108563806A (en) * 2018-01-05 2018-09-21 哈尔滨工业大学(威海) Engine air passage parameter long-range forecast method based on similitude and system
CN208537620U (en) * 2018-05-25 2019-02-22 国网宁夏电力有限公司石嘴山供电公司 Secondary loop of mutual inductor earth current inspection device
CN111818636A (en) * 2020-06-03 2020-10-23 哈尔滨工业大学(威海) Vehicle-mounted Bluetooth positioning system and positioning method thereof
CN112910089A (en) * 2021-01-25 2021-06-04 国网山东省电力公司青岛供电公司 Transformer substation secondary equipment fault logic visualization method and system

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102749549A (en) * 2012-06-21 2012-10-24 云南电力试验研究院(集团)有限公司电力研究院 Intelligent checking system for substation alternating current voltage secondary circuit
CN104699077A (en) * 2015-02-12 2015-06-10 浙江大学 Nested iterative fisher discriminant analysis-based fault diagnosis isolation method
CN105425768A (en) * 2015-11-06 2016-03-23 国网山东莒县供电公司 Electric power secondary equipment monitoring device and method
CN107958292A (en) * 2017-10-19 2018-04-24 山东科技大学 Transformer based on cost sensitive learning obscures careful reasoning method for diagnosing faults
CN108563806A (en) * 2018-01-05 2018-09-21 哈尔滨工业大学(威海) Engine air passage parameter long-range forecast method based on similitude and system
CN108418304A (en) * 2018-03-09 2018-08-17 国家电网公司 Substation secondary circuit state monitoring method, apparatus and system
CN208537620U (en) * 2018-05-25 2019-02-22 国网宁夏电力有限公司石嘴山供电公司 Secondary loop of mutual inductor earth current inspection device
CN111818636A (en) * 2020-06-03 2020-10-23 哈尔滨工业大学(威海) Vehicle-mounted Bluetooth positioning system and positioning method thereof
CN112910089A (en) * 2021-01-25 2021-06-04 国网山东省电力公司青岛供电公司 Transformer substation secondary equipment fault logic visualization method and system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
林君豪等: "基于宏微观特征分层聚类的配电网拓扑相似性分析方法", 《电力系统自动化》 *
王红等: "故障事件严酷度类别的模糊模式识别法", 《故障事件严酷度类别的模糊模式识别法 *
赵渊等: "电网可靠性评估中计及加和特性的非参数解集负荷模型", 《中国电机工程学报》 *

Similar Documents

Publication Publication Date Title
Jarventausta et al. Using fuzzy sets to model the uncertainty in the fault location process of distribution networks
Shadi et al. A real-time hierarchical framework for fault detection, classification, and location in power systems using PMUs data and deep learning
US10422833B2 (en) Systems and methods for classifying power line events
CN111443259A (en) Active power distribution network fault diagnosis method and system based on local abnormal factor detection
CN109886475B (en) Information security situation perception system of measurement automation system based on AI
CN107121639B (en) A kind of multi-Dimensional parameters direct current system accumulator cell management method and device
CN104134999A (en) Power-distribution-network measurement effectiveness analysis practical calculation method based on multiple data sources
CN112098915A (en) Method for evaluating secondary errors of multiple voltage transformers under double-bus segmented wiring
US20190056436A1 (en) Similarity detection of abnormal waveforms using pmu measurement
CN112751418A (en) Intelligent power distribution network regional situation element sensing method and system
US20230072795A1 (en) Online State Estimation and Topology Identification Using Advanced Metering Infrastructure (AMI) Measurements
CN114781476B (en) System and method for analyzing faults of measuring equipment
CN115587673B (en) Voltage transformer error state prediction method and system
CN109752629A (en) A kind of power grid measurement problem intelligent diagnosing method and system
CN117220441A (en) Information monitoring intelligent early warning method based on panoramic monitoring
CN103311960A (en) Method for dividing coherent stability regions of strong-connection power grids
Chantler et al. The use of fault-recorder data for diagnosing timing and other related faults in electricity transmission networks
CN112348207A (en) Power grid disaster prevention early warning method and device
CN113484646A (en) Transformer substation secondary circuit complex fault positioning device and diagnosis method
CN115685045A (en) Online evaluation method for voltage transformer
CN115833387A (en) Automatic state inspection method for energy storage power station
CN115561550A (en) Data diagnosis method for direct current system of transformer substation
CN111313355B (en) Method for updating monitoring signal event rule under manual supervision
CN117559449B (en) Power grid power transmission stability control method, system and storage medium
Du et al. Fault diagnosis method of automation equipment in independent and controllable substation based on deep reinforcement learning

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20211008

RJ01 Rejection of invention patent application after publication