Summary of the invention
Embodiment of the present invention technical matters to be solved is, a kind of method and system predicting power equipments defect is provided, the contingent defect of power equipment can be gone out by accurately predicting, power equipment operation maintenance personnel is helped to carry out emphasis tour and maintenance targetedly, decrease the labor capacity of power equipment operation maintenance personnel, improve power equipments defect analysis and effect is instructed to actual production work.
In order to solve the problems of the technologies described above, embodiments provide a kind of method predicting power equipments defect, it realizes on N number of of a sort power equipment, and described method comprises:
A, extract the historical data of N number of power equipment described in the time period, be met K defect classification of predetermined condition according to described historical data statistics, and using the described K obtained a defect classification all as defect analysis variable; Wherein, K, N are positive integer;
B, in described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, be combined into vectorial P (0) by a described K probability further;
C, obtain described K the total degree that defect analysis variable occurs respectively in described historical data, and the described time period is divided into M timeslice, sequentially determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode, and according to the total degree of K defect analysis occurrences of described acquisition and the number of times of K*K kind defect analysis variable transitions mode, determine to be combined into matrix P by a described K*K transition probability further by the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively; Wherein, M is positive integer;
D, according to formula P (1)=P (0) * P, determine the maximal value in P (1), and using defect classification corresponding for the described maximal value determined as the prediction defect before new power equipment comes into operation.
Wherein, described method comprises further:
Obtain a defect having come into operation the current generation of power equipment, and find corresponding to the maximal value in K kind conversion regime in described matrix P according to the classification of described current defect, and using defect classification corresponding for the described maximal value found next prediction defect as the power equipment of described current generation defect.
Wherein, the concrete steps of described step a comprise:
Extract the historical data of N number of power equipment described in the time period, and determine the total degree that in described historical data, each defect occurs;
The total degree sequential from large to small that each defect is occurred, filter out the defect that a front K-1 total degree is large, and classification corresponding for a described K-1 defect and classification that defect do not occur are set to defect analysis variable further as K the defect classification satisfied condition.
Wherein, the concrete steps of described step b comprise:
In described historical data, obtain the lasting time of origin of described N number of power equipment each defect analysis variable corresponding, under filtering out same defect analysis variable described N number of power equipment duration and add up, obtain the duration that K defect analysis variable is corresponding respectively;
The duration that K the defect analysis variable obtained described in cumulative is corresponding respectively, obtain total duration;
The duration that the described K obtained a defect analysis variable is corresponding respectively is all divided by with the total duration of described acquisition, obtains the probability that K defect analysis variable is corresponding respectively, and a described K probability is combined into vectorial P (0).
Wherein, the concrete steps of described step c comprise:
Obtain described K the total degree that defect analysis variable occurs respectively in described historical data;
Described K defect analysis variable is mapped between two, obtains K*K kind defect analysis variable transitions mode;
Is divided into M timeslice the described time period, sorts from small to large according to the time, determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode;
Using the number of times of each defect analysis variable transitions mode all as molecule, and determine the defect analysis variable that in every a part, main mapping pair is answered, and the total degree filtering out the defect analysis occurrences that main mapping pair is answered in described every a part is as corresponding denominator, obtain the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively, further a described K*K transition probability is combined into matrix P.
Wherein, described power equipment is transformer, described defect analysis variable comprise do not occur defect, leakage of oil, cooling system failure, instrument fault, operating mechanism exception and exterior mechanical damage.
The embodiment of the present invention further provides a kind of method predicting power equipments defect, and it realizes on N number of of a sort power equipment, and described method comprises:
S1, extract the historical data of N number of power equipment described in the time period, be met K defect classification of predetermined condition according to described historical data statistics, and using the described K obtained a defect classification all as defect analysis variable; Wherein, K, N are positive integer;
S2, obtain described K the total degree that defect analysis variable occurs respectively in described historical data, and the described time period is divided into M timeslice, sequentially determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode, and according to the total degree of K defect analysis occurrences of described acquisition and the number of times of K*K kind defect analysis variable transitions mode, determine to be combined into matrix P by a described K*K transition probability further by the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively; Wherein, M is positive integer;
S3, obtain a defect having come into operation the current generation of power equipment, and find corresponding to the maximal value in K kind conversion regime in described matrix P according to the classification of described current defect, and using defect classification corresponding for the described maximal value found next prediction defect as the power equipment of described current generation defect.
Wherein, the concrete steps of described step S1 comprise:
Extract the historical data of N number of power equipment described in the time period, and determine the total degree that in described historical data, each defect occurs;
The total degree sequential from large to small that each defect is occurred, filter out the defect that a front K-1 total degree is large, and classification corresponding for a described K-1 defect and classification that defect do not occur are set to defect analysis variable further as K the defect classification satisfied condition.
Wherein, the concrete steps of described step S2 comprise:
Obtain described K the total degree that defect analysis variable occurs respectively in described historical data;
Described K defect analysis variable is mapped between two, obtains K*K kind defect analysis variable transitions mode;
Is divided into M timeslice the described time period, sorts from small to large according to the time, determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode;
Using the number of times of each defect analysis variable transitions mode all as molecule, and determine the defect analysis variable that in every a part, main mapping pair is answered, and the total degree filtering out the defect analysis occurrences that main mapping pair is answered in described every a part is as corresponding denominator, obtain the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively, further a described K*K transition probability is combined into matrix P.
Wherein, described method comprises further:
In described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, be combined into vectorial P (0) by a described K probability further;
According to formula P (1)=P (0) * P, determine the maximal value in P (1), and using defect classification corresponding for the described maximal value determined as the prediction defect before new power equipment comes into operation.
Wherein, described in described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, the concrete steps being combined into vectorial P (0) by a described K probability further comprise:
In described historical data, obtain the lasting time of origin of described N number of power equipment each defect analysis variable corresponding, under filtering out same defect analysis variable described N number of power equipment duration and add up, obtain the duration that K defect analysis variable is corresponding respectively;
The duration that K the defect analysis variable obtained described in cumulative is corresponding respectively, obtain total duration;
The duration that the described K obtained a defect analysis variable is corresponding respectively is all divided by with the total duration of described acquisition, obtains the probability that K defect analysis variable is corresponding respectively, and a described K probability is combined into vectorial P (0).
The embodiment of the present invention additionally provides a kind of system predicting power equipments defect, and it realizes on N number of of a sort power equipment, and described system comprises:
Determine defect analysis variable cell, for extracting the historical data of N number of power equipment described in the time period, K defect classification of predetermined condition is met according to described historical data statistics, and using the described K obtained a defect classification all as defect analysis variable; Wherein, K, N are positive integer;
Obtain probability unit, for in described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, be combined into vectorial P (0) by a described K probability further;
Obtain transition probability unit, for obtaining described K the total degree that defect analysis variable occurs respectively in described historical data, and the described time period is divided into M timeslice, sequentially determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode, and according to the total degree of K defect analysis occurrences of described acquisition and the number of times of K*K kind defect analysis variable transitions mode, determine the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively, matrix P is combined into further by a described K*K transition probability, wherein, M is positive integer,
New equipment prediction defective unit, for according to formula P (1)=P (0) * P, determine the maximal value in P (1), and using defect classification corresponding for the described maximal value determined as the prediction defect before new power equipment comes into operation.
Wherein, described system also comprises:
Operational outfit prediction defective unit, for obtaining a defect having come into operation the current generation of power equipment, and find corresponding to the maximal value in K kind conversion regime in described matrix P according to the classification of described current defect, and using defect classification corresponding for the described maximal value found next prediction defect as the power equipment of described current generation defect.
The embodiment of the present invention further provides a kind of system predicting power equipments defect, and it realizes on N number of of a sort power equipment, and described system comprises:
Determine defect analysis variable cell, for extracting the historical data of N number of power equipment described in the time period, K defect classification of predetermined condition is met according to described historical data statistics, and using the described K obtained a defect classification all as defect analysis variable; Wherein, K, N are positive integer;
Obtain transition probability unit, for obtaining described K the total degree that defect analysis variable occurs respectively in described historical data, and the described time period is divided into M timeslice, sequentially determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode, and according to the total degree of K defect analysis occurrences of described acquisition and the number of times of K*K kind defect analysis variable transitions mode, determine the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively, matrix P is combined into further by a described K*K transition probability, wherein, M is positive integer,
Operational outfit prediction defective unit, for obtaining a defect having come into operation the current generation of power equipment, and find corresponding to the maximal value in K kind conversion regime in described matrix P according to the classification of described current defect, and using defect classification corresponding for the described maximal value found next prediction defect as the power equipment of described current generation defect.
Wherein, described system also comprises:
Obtain probability unit, for in described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, be combined into vectorial P (0) by a described K probability further;
New equipment prediction defective unit, for according to formula P (1)=P (0) * P, determine the maximal value in P (1), and using defect classification corresponding for the described maximal value determined as the prediction defect before new power equipment comes into operation.
Implement the embodiment of the present invention, there is following beneficial effect:
In embodiments of the present invention, due to according to the historic defects information of similar power equipment and current operating situation, adopt Prediction of Markov algorithm, accurately predicting goes out power equipment (comprising does not come into operation and come into operation, and defect occurred) contingent defect, power equipment operation maintenance personnel can be helped to carry out emphasis tour and maintenance targetedly, reduce the labor capacity of power equipment operation maintenance personnel, improve power equipments defect analysis and effect is instructed to actual production work.
Embodiment
For making the object, technical solutions and advantages of the present invention clearly, below in conjunction with accompanying drawing, the present invention is described in further detail.
As shown in Figure 1, the embodiment of the present invention provides a kind of method predicting power equipments defect, and it realizes on N number of of a sort power equipment, and described method comprises:
Step S101, extract the historical data of N number of power equipment described in the time period, be met K defect classification of predetermined condition according to described historical data statistics, and using the described K obtained a defect classification all as defect analysis variable; Wherein, K, N are positive integer;
Detailed process is, extracts the historical data of N number of power equipment in the time period, and determines the total degree that in historical data, each defect occurs;
The total degree sequential from large to small that each defect is occurred, filter out the defect that a front K-1 total degree is large, and classification corresponding for K-1 defect and classification that defect do not occur are set to defect analysis variable further as K the defect classification satisfied condition.
As an example, power equipment is transformer, its main cause that defect occurs includes but not limited to leakage of oil, cooling system failure, instrument fault, operating mechanism exception, exterior mechanical damage etc., therefore, when extracting the service data in multiple transformer 1 year, this service data comprises the defective data that there is above-mentioned defect.Count the total degree that in all transformers in the historical data, all defect occurs respectively, then all defect is sorted according to total degree is descending, filtering out and coming several, foremost defect is that leakage of oil, cooling system failure, instrument fault, operating mechanism exception and exterior mechanical damage, and using the defect classification filtered out and there is not defect classification as corresponding defect analysis variable, the order of defect analysis variable-definition for there is not defect, go out leakage of oil, cooling system failure, instrument fault, operating mechanism exception and exterior mechanical damage.
Step S102, in described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, be combined into vectorial P (0) by a described K probability further;
Detailed process is, in the historical data, obtain the lasting time of origin of N number of power equipment each defect analysis variable corresponding, filter out the duration of N number of power equipment under same defect analysis variable and add up, obtain the duration of K defect analysis variable correspondence respectively;
The duration that cumulative K the defect analysis variable obtained is corresponding respectively, obtain total duration;
Duration corresponding respectively for the K obtained a defect analysis variable is all divided by with the total duration obtained, obtains the probability that K defect analysis variable is corresponding respectively, and K probability is combined into vectorial P (0).
In embodiments of the present invention, K probability can be expressed as p
1(0), p
2(0) ... p
k(0), therefore vectorial P (0)=(p
1(0), p
2(0) ... p
k(0)).
For aforementioned transformer, p
1(0) probability that defect does not occur is represented, p
2(0) ... p
6(0) probability of the defects such as leakage of oil, cooling system failure, instrument fault, operating mechanism exception and exterior mechanical damage is represented respectively.The step calculating the probability of defect analysis variable is following (with p
1(0) be example):
There is not the duration of the duration of defect, the duration of cumulative all Leakage in Transformer oil, duration of cumulative all transformer cooling system accident defects, the duration of cumulative all transformer instrument faults, cumulative all transformer operating mechanisms exception in cumulative all transformers, till the result of the duration that all transformer exterior mechanicals that obtains adding up out damage;
The duration that above-mentioned six defect analysis variablees add up out is totally added again, thus obtains total duration;
All transformers be there is not the duration of defect as molecule, total duration is calculated as denominator, obtains p
1(0).
In like manner, p can be obtained according to said method
2(0) ... p
6, thus further form vectorial P (0)=(p (0)
1(0), p
2(0) ... p
6(0)).
Step S103, obtain described K the total degree that defect analysis variable occurs respectively in described historical data, and the described time period is divided into M timeslice, sequentially determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode, and according to the total degree of K defect analysis occurrences of described acquisition and the number of times of K*K kind defect analysis variable transitions mode, determine to be combined into matrix P by a described K*K transition probability further by the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively; Wherein, M is positive integer;
Detailed process is, obtains the total degree that K defect analysis variable occurs respectively in the historical data;
K defect analysis variable is mapped between two, obtains K*K kind defect analysis variable transitions mode;
M timeslice will be divided into the time period, sort from small to large according to the time, determine that between adjacent time sheet, N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode;
Using the number of times of each defect analysis variable transitions mode all as molecule, and determine the defect analysis variable that in every a part, main mapping pair is answered, and the total degree filtering out the defect analysis occurrences that main mapping pair in every a part is answered is as corresponding denominator, obtain the transition probability that K*K kind defect analysis variable transitions mode is corresponding respectively, further K*K transition probability is combined into matrix P.
It should be noted that the demand in order to calculate, K defect analysis variable maps between two, and it comprises the mutual mapping of same defect analysis variable, and the leakage of oil variable mappings as transformer is identical leakage of oil variable, the matrix therefore formed.
Should be noted that, M the timeslice be divided into time period is little as far as possible, adjacent time sheet (a upper timeslice is to future time sheet) is occurred, and the number of times of defect class transitions meets for once, as above there is leakage of oil in a timeslice A transformer, but only there will be primary cooling system accident defect at future time sheet A transformer, instead of once above cooling system failure defect.
In embodiments of the present invention, the transition probability that K*K kind defect analysis variable transitions mode is corresponding respectively can be expressed as (p
11, p
12... p
1K; p
21, p
22... p
2K; ... p
k1, p
k2... p
kK), thus form matrix
Wherein, p
nnthe defect classification being expressed as the defect analysis variable of the same power equipment of adjacent time sheet corresponding is identical, n=1, and 2 ... K.
For aforementioned transformer, 1 year was divided according to 12 months, obtain 12 timeslices; p
11represent that the equipment deficiency state between upper and lower two months is all the transition probability that defect does not occur, p
12represent that the transition probability that defect and next timeslice there occurs leakage of oil defect does not occur a timeslice, therefore by that analogy, the transition probability that the defect analysis variable transitions mode of transformer is corresponding respectively can be expressed as (p
11, p
12... p
16; p
21, p
22... p
26; ... p
61, p
62... p
66).The step calculating the transition probability of defect analysis variable is following (there is not defect (p
11, p
12... p
16) be example):
Obtain the total degree that six defect analysis variablees occur respectively in the historical data, certainly now for (p
11, p
12... p
16), only need to determine that the total degree of defect appearance does not occur all transformers in the historical data;
There is not defect will map out six kinds of conversion regimes and be: defect does not occur and maps and defect does not occur, defect does not occur map leakage of oil, defect does not occur map cooling system failure, do not occur that defect maps instrument fault, that defect map operation mechanism does not occur is abnormal and defect does not occur map exterior mechanical and damage;
From sequence in one to Dec, determine that between adjacent time sheet, all transformers correspond to the number of times of above-mentioned six kinds of defect analysis variable transitions modes;
Using the number of times of above-mentioned six kinds of defect analysis variable transitions modes successively all as molecule, in historical data there is not the total degree of defect appearance as denominator in all transformers, obtains (p
11, p
12... p
16).
In like manner, (p can be obtained according to said method
21, p
22... p
26; ... p
61, p
62... p
66), thus form matrix further
Step S104, according to formula P (1)=P (0) * P, determine the maximal value in P (1), and using defect classification corresponding for the described maximal value determined as the prediction defect before new power equipment comes into operation.
By predict new power equipment come into operation before defect power equipment maintainer can be made to carry out emphasis tour and maintenance targetedly, reduce the labor capacity of power equipment operation maintenance personnel, improve power equipments defect analysis and effect is instructed to actual production work.
For aforementioned transformer,
The maximal value obtained in P (1) corresponds to cooling system failure defect, using the defect of cooling system failure defect as the transformer of prediction.
Predict in the embodiment of the present invention that transition matrix that the method for power equipments defect produces also has the failure prediction of the power equipment come into operation and effect is instructed to actual production work, and do not need the calculated amount increasing data to carry out failure prediction to current power equipment under certain condition, therefore described method comprises further:
Obtain a defect having come into operation the current generation of power equipment, and find in matrix P corresponding to the maximal value in K kind conversion regime according to the classification of this current defect, and using defect classification corresponding for the maximal value that finds next prediction defect as the power equipment of current generation defect.
As shown in Figure 2, embodiments provide the method for another kind of prediction power equipments defect, it realizes on N number of of a sort power equipment, and described method comprises:
Step S201, extract the historical data of N number of power equipment described in the time period, be met K defect classification of predetermined condition according to described historical data statistics, and using the described K obtained a defect classification all as defect analysis variable; Wherein, K, N are positive integer;
Detailed process is, extracts the historical data of N number of power equipment in the time period, and determines the total degree that in historical data, each defect occurs;
The total degree sequential from large to small that each defect is occurred, filter out the defect that a front K-1 total degree is large, and classification corresponding for K-1 defect and classification that defect do not occur are set to defect analysis variable further as K the defect classification satisfied condition.
Step S202, obtain described K the total degree that defect analysis variable occurs respectively in described historical data, and the described time period is divided into M timeslice, sequentially determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode, and according to the total degree of K defect analysis occurrences of described acquisition and the number of times of K*K kind defect analysis variable transitions mode, determine to be combined into matrix P by a described K*K transition probability further by the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively; Wherein, M is positive integer;
Detailed process is, obtains the total degree that K defect analysis variable occurs respectively in the historical data;
K defect analysis variable is mapped between two, obtains K*K kind defect analysis variable transitions mode;
M timeslice will be divided into the time period, sort from small to large according to the time, determine that between adjacent time sheet, N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode;
Using the number of times of each defect analysis variable transitions mode all as molecule, and determine the defect analysis variable that in every a part, main mapping pair is answered, and the total degree filtering out the defect analysis occurrences that main mapping pair in every a part is answered is as corresponding denominator, obtain the transition probability that K*K kind defect analysis variable transitions mode is corresponding respectively, further K*K transition probability is combined into matrix P.
It should be noted that the demand in order to calculate, K defect analysis variable maps between two, and it comprises the mutual mapping of same defect analysis variable, and the leakage of oil variable mappings as transformer is identical leakage of oil variable, the matrix therefore formed.
Should be noted that, M the timeslice be divided into time period is little as far as possible, adjacent time sheet (a upper timeslice is to future time sheet) is occurred, and the number of times of defect class transitions meets for once, as above there is leakage of oil in a timeslice A transformer, but only there will be primary cooling system accident defect at future time sheet A transformer, instead of once above cooling system failure defect.
In embodiments of the present invention, the transition probability that K*K kind defect analysis variable transitions mode is corresponding respectively can be expressed as (p
11, p
12... p
1K; p
21, p
22... p
2K; ... p
k1, p
k2... p
kK), thus form matrix
Wherein, p
nnthe defect classification being expressed as the defect analysis variable of the same power equipment of adjacent time sheet corresponding is identical, n=1, and 2 ... K.
Step S203, obtain a defect having come into operation the current generation of power equipment, and find corresponding to the maximal value in K kind conversion regime in described matrix P according to the classification of described current defect, and using defect classification corresponding for the described maximal value found next prediction defect as the power equipment of described current generation defect.
By prediction come into operation power equipment come into operation before defect power equipment maintainer can be made to carry out emphasis tour and maintenance targetedly, reduce the labor capacity of power equipment operation maintenance personnel, improve power equipments defect analysis and effect is instructed to actual production work.
Predict in the embodiment of the present invention that new power equipment defect before being taken into use also has and instruct effect to actual production work, therefore described method comprises further:
Step S21, in described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, be combined into vectorial P (0) by a described K probability further;
Detailed process is, in the historical data, obtain the lasting time of origin of N number of power equipment each defect analysis variable corresponding, filter out the duration of N number of power equipment under same defect analysis variable and add up, obtain the duration of K defect analysis variable correspondence respectively;
The duration that cumulative K the defect analysis variable obtained is corresponding respectively, obtain total duration;
Duration corresponding respectively for the K obtained a defect analysis variable is all divided by with the total duration obtained, obtains the probability that K defect analysis variable is corresponding respectively, and K probability is combined into vectorial P (0).
In embodiments of the present invention, K probability can be expressed as p
1(0), p
2(0) ... p
k(0), therefore vectorial P (0)=(p
1(0), p
2(0) ... p
k(0)).
Step S22, according to formula P (1)=P (0) * P, determine the maximal value in P (1), and using defect classification corresponding for the described maximal value determined as the prediction defect before new power equipment comes into operation.
Wherein, described in described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, the concrete steps being combined into vectorial P (0) by a described K probability further comprise:
As shown in Figure 3, the embodiment of the present invention additionally provides a kind of system predicting power equipments defect, and it realizes on N number of of a sort power equipment, and described system comprises:
Determine defect analysis variable cell 110, for extracting the historical data of N number of power equipment described in the time period, K defect classification of predetermined condition is met according to described historical data statistics, and using the described K obtained a defect classification all as defect analysis variable; Wherein, K, N are positive integer;
Obtain probability unit 120, for in described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, be combined into vectorial P (0) by a described K probability further;
Obtain transition probability unit 130, for obtaining described K the total degree that defect analysis variable occurs respectively in described historical data, and the described time period is divided into M timeslice, sequentially determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode, and according to the total degree of K defect analysis occurrences of described acquisition and the number of times of K*K kind defect analysis variable transitions mode, determine the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively, matrix P is combined into further by a described K*K transition probability, wherein, M is positive integer,
New equipment prediction defective unit 140, for according to formula P (1)=P (0) * P, determine the maximal value in P (1), and using defect classification corresponding for the described maximal value determined as the prediction defect before new power equipment comes into operation.
Wherein, system also comprises:
Operational outfit prediction defective unit 150, for the defect of the current generation of power equipment of having come into operation when acquisition one, and find corresponding to the maximal value in K kind conversion regime in described matrix P according to the classification of described current defect, and using defect classification corresponding for the described maximal value found next prediction defect as the power equipment of described current generation defect.
Wherein, obtain probability unit 120 to comprise:
Obtain duration module 1201, for in described historical data, obtain the lasting time of origin of described N number of power equipment each defect analysis variable corresponding, under filtering out same defect analysis variable described N number of power equipment duration and add up, obtain the duration that K defect analysis variable is corresponding respectively;
Obtain total duration module 1202, the duration that K the defect analysis variable for obtaining described in cumulative is corresponding respectively, acquisition total duration;
Determine probability module 1203, all be divided by with the total duration of described acquisition for the duration that the described K obtained a defect analysis variable is corresponding respectively, obtain the probability that K defect analysis variable is corresponding respectively, and a described K probability is combined into vectorial P (0).
Wherein, obtain transition probability unit 130 to comprise:
Obtain defect total degree module 1301, for obtaining described K the total degree that defect analysis variable occurs respectively in described historical data;
Defect mapping block 1302, for mapping between two described K defect analysis variable, obtains K*K kind defect analysis variable transitions mode;
Obtain defect conversion times module 1303, for the described time period is divided into M timeslice, sort from small to large according to the time, determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode;
Determine transition probability module 1304, for using the number of times of each defect analysis variable transitions mode all as molecule, and determine the defect analysis variable that in every a part, main mapping pair is answered, and the total degree filtering out the defect analysis occurrences that main mapping pair is answered in described every a part is as corresponding denominator, obtain the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively, further a described K*K transition probability is combined into matrix P.
As shown in Figure 4, the embodiment of the present invention additionally provides the system of another kind of prediction power equipments defect, and it realizes on N number of of a sort power equipment, and described system comprises:
Determine defect analysis variable cell 210, for extracting the historical data of N number of power equipment described in the time period, K defect classification of predetermined condition is met according to described historical data statistics, and using the described K obtained a defect classification all as defect analysis variable; Wherein, K, N are positive integer;
Obtain transition probability unit 220, for obtaining described K the total degree that defect analysis variable occurs respectively in described historical data, and the described time period is divided into M timeslice, sequentially determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode, and according to the total degree of K defect analysis occurrences of described acquisition and the number of times of K*K kind defect analysis variable transitions mode, determine the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively, matrix P is combined into further by a described K*K transition probability, wherein, M is positive integer,
Operational outfit prediction defective unit 230, for obtaining a defect having come into operation the current generation of power equipment, and find corresponding to the maximal value in K kind conversion regime in described matrix P according to the classification of described current defect, and using defect classification corresponding for the described maximal value found next prediction defect as the power equipment of described current generation defect.
Wherein, described acquisition transition probability unit 220 comprises:
Obtain defect total degree module 2201, for obtaining described K the total degree that defect analysis variable occurs respectively in described historical data;
Defect mapping block 2202, for mapping between two described K defect analysis variable, obtains K*K kind defect analysis variable transitions mode;
Obtain defect conversion times module 2203, for the described time period is divided into M timeslice, sort from small to large according to the time, determine that between adjacent time sheet, described N number of power equipment corresponds to the number of times of K*K kind defect analysis variable transitions mode;
Determine transition probability module 2204, for using the number of times of each defect analysis variable transitions mode all as molecule, and determine the defect analysis variable that in every a part, main mapping pair is answered, and the total degree filtering out the defect analysis occurrences that main mapping pair is answered in described every a part is as corresponding denominator, obtain the transition probability that described K*K kind defect analysis variable transitions mode is corresponding respectively, further a described K*K transition probability is combined into matrix P.
Wherein, described system also comprises:
Obtain probability unit 240, for in described historical data, all add up out the duration of described N number of power equipment under same defect analysis variable for each defect analysis variable, and the duration of cumulative K described N number of power equipment under same defect analysis variable and obtain total duration, and according to the described total duration that adds up out and duration corresponding to described each the defect analysis variable added up out, determine the probability of each defect state variable, be combined into vectorial P (0) by a described K probability further;
New equipment prediction defective unit 250, for according to formula P (1)=P (0) * P, determine the maximal value in P (1), and using defect classification corresponding for the described maximal value determined as the prediction defect before new power equipment comes into operation.
Wherein, described acquisition probability unit 240 comprises:
Obtain duration module 2401, for in described historical data, obtain the lasting time of origin of described N number of power equipment each defect analysis variable corresponding, under filtering out same defect analysis variable described N number of power equipment duration and add up, obtain the duration that K defect analysis variable is corresponding respectively;
Obtain total duration module 2402, the duration that K the defect analysis variable for obtaining described in cumulative is corresponding respectively, acquisition total duration;
Determine probability module 2403, all be divided by with the total duration of described acquisition for the duration that the described K obtained a defect analysis variable is corresponding respectively, obtain the probability that K defect analysis variable is corresponding respectively, and a described K probability is combined into vectorial P (0).
Implement the embodiment of the present invention, there is following beneficial effect:
In embodiments of the present invention, due to according to the historic defects information of similar power equipment and current operating situation, adopt Prediction of Markov algorithm, accurately predicting goes out power equipment (comprising does not come into operation and come into operation, and defect occurred) contingent defect, power equipment operation maintenance personnel can be helped to carry out emphasis tour and maintenance targetedly, reduce the labor capacity of power equipment operation maintenance personnel, improve power equipments defect analysis and effect is instructed to actual production work.
It should be noted that in said system embodiment, each included system unit is carry out dividing according to function logic, but is not limited to above-mentioned division, as long as can realize corresponding function; In addition, the concrete title of each functional unit, also just for the ease of mutual differentiation, is not limited to protection scope of the present invention.
One of ordinary skill in the art will appreciate that all or part of step realized in above-described embodiment method is that the hardware that can carry out instruction relevant by program has come, described program can be stored in a computer read/write memory medium, described storage medium, as ROM/RAM, disk, CD etc.
Above disclosedly be only present pre-ferred embodiments, certainly can not limit the interest field of the present invention with this, therefore according to the equivalent variations that the claims in the present invention are done, still belong to the scope that the present invention is contained.