CN102680243A - Online judgment method for steam flow shock excitation fault of steam turbine generator unit - Google Patents

Abstract

The invention discloses an online judgment method for a steam flow shock excitation fault of a steam turbine generator unit in the technical field of state monitoring and fault diagnosis of rotary mechanical vibration. The online judgment method comprises the following steps of: acquiring relative vibration data of a shaft for supporting a bearing on one side of a high pressure rotor of the steam turbine generator unit, a rotating speed signal of a rotor, a key-phase signal of the rotor and power data of the unit; forming a unit power data sequence, a low-frequency vibration amplitude maximum sequence and a entropy sequence of the low-frequency vibration amplitude sequence; calculating an increasing trend parameter of the unit power data sequence, a maximum value of the low-frequency vibration amplitude maximum sequence, an entropy sequence skewness parameter of the low-frequency vibration amplitude sequence and Kendall relevant parameters of the entropy sequence of the low-frequency vibration amplitude sequence and the unit power data sequence; and finally, judging whether the steam flow shock excitation fault occurs or not by using the result. According to the online judgment method disclosed by the invention, the automatic real-time online monitoring and judgment of the steam flow shock excitation fault of the high pressure rotor are realized and the efficiency and accuracy of analysis and diagnosis are improved.

Description

Turbo-generator Set steam flow excitation On-line Fault method of discrimination

Technical field

The invention belongs to rotating machinery vibrating condition monitoring and fault diagnosis technical field, relate in particular to a kind of Turbo-generator Set steam flow excitation On-line Fault method of discrimination.

Background technology

Steam flow excitation be a kind of usually occur in the large-size steam turbine height (in) press epitrochanterian, the low-frequency vibration phenomenon of bringing out by the steam exciting force.Because high pressure high temperature turbosets are along with the raising of turbine steam condition; Can cause the increase of high pressure cylinder admission density, flow velocity to improve; Sensitivity to moderate improves the tangential force of vapor action on high pressure rotor to dynamic and static gaps, hermetically-sealed construction and rotor and cylinder; Increase the exciting force that acts on high pressure rotor, directly influenced the available rate of unit.

Unit generation steam flow excitation fault often shows as low-frequency vibration and becomes suddenly the high load capacity stage that big and steam flow unit fault often occurs in unit.Judge whether unit the steam flow excitation fault takes place; Usually accomplish by professional with certain field operation experiences and professional knowledge technical ability; Bring thus that analysis result is higher to personnel's subjectivity degree of dependence, the analytic process labor intensive resource time; Problems such as analytical work cost height, and can't accomplish real-time automatic on-line monitoring, analysis and the differentiation of steam flow excitation fault.Therefore, propose a kind of large turbo-type generator group steam flow excitation On-line Fault method of discrimination and just seem very important.

Large turbo-type generator group steam flow excitation On-line Fault method of discrimination provided by the invention; To unit operation rotor axle vibrate relatively, data such as the power of the assembling unit carry out real-time automatic on-line monitoring, analyze and differentiate; Judge whether high pressure rotor the steam flow excitation fault takes place, improve the efficient and the accuracy of high pressure rotor steam flow excitation fault analysis and diagnosis work.

Summary of the invention

The objective of the invention is to, propose a kind of Turbo-generator Set steam flow excitation On-line Fault method of discrimination, in order to realize real-time automatic on-line monitoring, analysis and the differentiation of large turbo-type generator group steam flow excitation fault.

For realizing above-mentioned purpose, technical scheme provided by the invention is that a kind of Turbo-generator Set steam flow excitation On-line Fault method of discrimination is characterized in that said method comprises:

Step 1: set initial moment T _S, very first time stepping length t ₁, second time stepping length t ₂, the first setting value D ₁With the second setting value D ₂

Step 2: gather the relative vibration data of axle of Turbo-generator Set high pressure rotor one side radial journal bearing, the tach signal of rotor, the key signal and the power of the assembling unit data of rotor in real time;

Step 3: from initial moment T _SBeginning, every at a distance from very first time stepping length t ₁The power of the assembling unit data P of storage current time ^UEvery at a distance from second time stepping length t ₂, the maximal value A of calculating current time low-frequency vibration amplitude sequence ^Fmax, the entropy E of current time low-frequency vibration amplitude sequence and the variation kurtosis parameter κ of current time low-frequency vibration amplitude sequence ^A, and the maximal value A of storage current time low-frequency vibration amplitude sequence ^FmaxEntropy E with current time low-frequency vibration amplitude sequence;

As the variation kurtosis parameter κ that satisfies current time low-frequency vibration amplitude sequence ^AGreater than the first setting value D1 and current time and initial moment T _SDifference more than or equal to the second setting value D ₂With the first stepping length t ₁Product the time, with current time as stopping T constantly _N, with the power of the assembling unit data P of the current time of storing ^UAs stopping T constantly _NPower of the assembling unit data

Execution in step 4;

Step 4: choose and stop T constantly _NIn two moment before, be designated as first T constantly formerly respectively ₁With second T constantly formerly ₂, and satisfy

Step 5: will be from first T constantly formerly ₁Rise to stopping T constantly _NEnd, every at a distance from very first time stepping length t ₁The power of the assembling unit data P of storage ^ULine up power of the assembling unit data sequence according to the sequencing of storage time

Will be from second T constantly formerly ₂Rise to stopping T constantly _NEnd, every at a distance from the second stepping length t ₂The maximal value A of the low-frequency vibration amplitude sequence of storage ^FmaxLine up low-frequency vibration amplitude maximal value sequence according to the sequencing of storage time

Every at a distance from the second stepping length t ₂The entropy E of the low-frequency vibration amplitude sequence of storage lines up the entropy sequence { E of low-frequency vibration amplitude sequence according to the sequencing of storage time _j, $j=\mathrm{1,2},...,\frac{T_N-T_2}{t_2};$

Step 6: computer set power parameter and low-frequency vibration parameter comprise:

1) computer set power data sequence increases progressively the trend parameter I ^P

2) maximal value

of calculating low-frequency vibration amplitude maximal value sequence

3) the variation degree of bias parameter S of the entropy sequence of calculating low-frequency vibration amplitude sequence ^E

4) calculate the entropy sequence of low-frequency vibration amplitude sequence and the Kendall related coefficient τ of power of the assembling unit data sequence;

Step 7: according to stopping T constantly _NPower of the assembling unit data

Power of the assembling unit data sequence increase progressively the trend parameter I ^P, low-frequency vibration amplitude maximal value sequence maximal value The variation degree of bias parameter S of the entropy sequence of low-frequency vibration amplitude sequence ^EKendall related coefficient τ with the entropy sequence and the power of the assembling unit data sequence of low-frequency vibration amplitude sequence judges whether high pressure rotor the steam flow excitation fault takes place.

The entropy E of said calculating current time low-frequency vibration amplitude sequence adopts formula

$E=\underset{k=1}{\overset{l}{\mathrm{\Σ}}}\left[{\left(A_k^{\mathrm{freq}}\right)}^2\ln \left({\left(A_k^{\mathrm{freq}}\right)}^2\right)\right],$

Wherein,

is k data of low-frequency vibration amplitude sequence; K=1; 2; ... l; L is the data number of predefined low-frequency vibration amplitude sequence; And when stipulating as

,

The variation kurtosis parameter κ of said current time low-frequency vibration amplitude sequence ^AAdopt formula

${\mathrm{\κ}}^A=1/l\underset{k=1}{\overset{l}{\mathrm{\Σ}}}{(A_k^{\mathrm{freq}}-{\mathrm{\μ}}^A)}^4/{\left({\mathrm{\σ}}^A\right)}^4,$

Wherein, Be k data of low-frequency vibration amplitude sequence, k=1,2 ... l, l are the data number of predefined low-frequency vibration amplitude sequence, μ ^ABe the average of low-frequency vibration amplitude sequence, promptly

σ ^ABe the standard deviation of low-frequency vibration amplitude sequence, promptly

Said computer set power data sequence increase progressively the trend parameter I ^PAdopt formula

I ^P=S ^P/[1/2×n×(n-1)]，

Wherein, n is the data number of power of the assembling unit data sequence, S ^PIt is the Ser.No. of power of the assembling unit data sequence; Said Ser.No. is meant the right sum of order in the data sequence; Said order is to being meant that the front and back position of a logarithm is identical with size order in a data sequence, and promptly the number of front is less than the number of back.

The variation degree of bias parameter S of the entropy sequence of said calculating low-frequency vibration amplitude sequence ^EAdopt formula

$S^E=1/n\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{(E_j-{\mathrm{\μ}}^E)}^3/{\left({\mathrm{\σ}}^E\right)}^3,$

Wherein, E _jBe j data of the entropy sequence of low-frequency vibration amplitude sequence, j=1,2 ..., n, μ ^EBe the average of the entropy sequence of low-frequency vibration amplitude sequence, promptly

σ ^EBe the standard deviation of the entropy sequence of low-frequency vibration amplitude sequence, promptly N is the data number of the entropy sequence of low-frequency vibration amplitude sequence.

The entropy sequence of said calculating low-frequency vibration amplitude sequence and the Kendall related coefficient τ of power of the assembling unit data sequence adopt formula

τ=(n _s-n _d)/[1/2×(n ²-1)]，

Wherein, n is the entropy sequence of low-frequency vibration amplitude sequence or the data number of power of the assembling unit data sequence, n _sBe in the entropy sequence and power of the assembling unit data sequence of low-frequency vibration amplitude sequence, the right sum of sequence that mediation is arranged, n _dBe in the entropy sequence and power of the assembling unit data sequence of low-frequency vibration amplitude sequence, the sum that the sequence of anharmonic arrangement is right;

The sequence that said mediation is arranged is to being meant, for 2 sequence { A _kAnd { B _kIn any 2 sequences to (A _m, B _m) and (A _n, B _n), q>=m, n>=1, q is sequence { A _kAnd sequence { B _kThe data number, if satisfy A simultaneously _m>A _nAnd B _m>B _nPerhaps satisfy A simultaneously _m<a _nAnd B _m<b _n, then 2 sequences are to (A _m, B _m) and (A _n, B _n) be that the sequence of be in harmonious proportion arranging is right;

The sequence of said anharmonic arrangement is to being meant, for 2 sequence { A _kAnd { B _kIn any 2 sequences to (A _m, B _m) and (A _n, B _n), q>=m, n>=1, q is sequence { A _kAnd sequence { B _kThe data number, if satisfy A simultaneously _m>A _nAnd B _m<b _nPerhaps satisfy A simultaneously _m<a _nAnd B _m>B _n, then 2 sequences are to (A _m, B _m) and (A _n, B _n) be that the sequence of anharmonic arrangement is right.

Said step 7 specifically is to satisfy following 5 conditions when simultaneously:

1) stops T constantly _NPower of the assembling unit data

Greater than the 3rd setting value;

2) power of the assembling unit data sequence increases progressively the trend parameter I ^PGreater than the 4th setting value;

The maximal value of 3) low-frequency vibration amplitude maximal value sequence

is greater than the 5th setting value;

The variation degree of bias parameter of the entropy sequence of 4) low-frequency vibration amplitude sequence is greater than S ^EGreater than the 6th setting value;

The Kendall related coefficient τ of the entropy sequence of 5) low-frequency vibration amplitude sequence and power of the assembling unit data sequence is greater than the 7th setting value;

Then judge high pressure rotor generation steam flow excitation fault; Otherwise, judge that the steam flow excitation fault does not take place high pressure rotor.

Method provided by the invention is utilized data such as unit operation rotor axle vibrates relatively, the power of the assembling unit, through the computational discrimination high pressure rotor whether the steam flow excitation fault takes place, and has realized automatic time on-line monitoring, analysis and the differentiation of steam flow excitation fault.

Description of drawings

Fig. 1 is a Turbo-generator Set steam flow excitation On-line Fault method of discrimination process flow diagram;

Fig. 2 is that Turbo-generator Set steam flow excitation On-line Fault is differentiated synoptic diagram.

Embodiment

Below in conjunction with accompanying drawing, preferred embodiment is elaborated.Should be emphasized that following explanation only is exemplary, rather than in order to limit scope of the present invention and application thereof.

Embodiment

Fig. 1 is a Turbo-generator Set steam flow excitation On-line Fault method of discrimination process flow diagram.Among Fig. 1, Turbo-generator Set steam flow excitation On-line Fault method of discrimination comprises:

Step 101: set initial moment T _S=0 second, very first time stepping length t ₁=3 seconds, second time stepping length t ₂=0.1 second, the first setting value D ₁=the 30 and second setting value D ₂=300.

Step 102: gather the relative vibration data of axle of Turbo-generator Set high pressure rotor one side radial journal bearing, the tach signal of rotor, the key signal and the power of the assembling unit data of rotor in real time.

The tach signal of the relative vibration data of armature spindle, rotor and key signal can obtain from the supervisory instrument (TSI) of configuration Turbo-generator Set, and power of the assembling unit data-signal can obtain from the dcs (DCS) of configuration Turbo-generator Set.In the present embodiment; The tach signal of the relative vibration data of armature spindle, rotor and key signal are supervisory instrument (TSI) acquisitions from the configuration Turbo-generator Set, and power of the assembling unit data-signal is dcs (DCS) acquisition from the configuration Turbo-generator Set.Fig. 2 is that Turbo-generator Set steam flow excitation On-line Fault is differentiated synoptic diagram, and is as shown in Figure 2, in the slot that data collecting card insertion industrial microcomputer (IPC) provides.Requirement according to data collecting card; The data acquisition conditioning device is handled the relative vibration signal of axle, the tach signal of rotor, the key signal from Turbo-generator Set supervisory instrument (TSI), the vibration at high speed data collecting card in the tach signal of the relative vibration signal of axle after treatment, rotor, the key signal input IPC.Each passage technology parameter of vibrating data collection card is 50ks/s, 24bit.Simultaneously, the data acquisition conditioning device is handled the power of the assembling unit data-signal from Turbo-generator Set dcs (DCS), the data collecting card in the bearing oil temperature data signal input IPC after treatment.Each passage technology parameter of data collecting card is 1ks/s, 16bit.

According to the concrete Turbo-generator Set steam flow excitation On-line Fault discriminating program of this method design, online discriminating program is installed in the industrial microcomputer (IPC).Once diagnosis cyclic process in the Turbo-generator Set steam flow excitation On-line Fault discriminating program comprises that the real time data acquisition that relates in the diagnostic method calculates that storage, real time discriminating, power of the assembling unit data parameters are calculated, the low-frequency vibration parameter calculates in real time, related coefficient is calculated in real time and series of computation such as fault verification is analyzed link.

Step 103: from initial moment T _CBeginning in=0 second, every at a distance from very first time stepping length t ₁The power of the assembling unit data P of=3 seconds storage current times ^UAs shown in Figure 2, power of the assembling unit data obtain from the dcs (DCS) of configuration Turbo-generator Set.

Every at a distance from second time stepping length t ₂=0.1 second, calculate current time low-frequency vibration amplitude sequence.As shown in Figure 2, the routine analyzer in the industrial microcomputer (IPC) is gathered the tach signal and the key signal of near the relative vibration data of the axle that records the Turbo-generator Set high pressure rotor A side radial journal bearing, rotor in real time through adopting the vibration at high speed data collecting card.Each passage technology parameter of vibrating data collection card is 50ks/s, 24bit.The relative vibration data of axle to unit high pressure rotor A side utilizes the FFT frequency spectrum analysis method, calculates the pairing vibration amplitude data sequence of the different vibration frequencies of current time from the low frequency to the high frequency (amplitude unit is μ m, micron).For the vibration at high speed data collecting card; Each constantly can both collect the pairing vibration amplitude data sequence of different vibration frequencies from the low frequency to the high frequency; Therefrom intercepting obtains the current time low-frequency vibration amplitude sequence of frequency less than unit working speed respective frequencies (50Hz); Be designated as and can set vibrating data collection frequency and image data amount, make low-frequency vibration amplitude sequence data number l=98.Utilize the relative vibration data of axle, the tach signal of rotor and the key signal of rotor of the Turbo-generator Set high pressure rotor one side radial journal bearing of current time collection; Calculating current time low-frequency vibration amplitude sequence

has been the technology that those skilled in the art use always, repeats no more in the present invention.

In the low-frequency vibration amplitude sequence that obtains current time

After, can calculate the maximal value A of current time low-frequency vibration amplitude sequence ^Fmax, its computing formula does

In the present embodiment, l=98.

Low-frequency vibration amplitude sequence according to current time

Can also calculate the entropy E of current time low-frequency vibration amplitude sequence and the variation kurtosis parameter κ of current time low-frequency vibration amplitude sequence ^A

The entropy E that calculates current time low-frequency vibration amplitude sequence according to the low-frequency vibration amplitude sequence

of current time adopts formula

$E=\underset{k=1}{\overset{l}{\mathrm{\&Sigma;}}}\left[{\left(A_k^{\mathrm{freq}}\right)}^2\ln \left({\left(A_k^{\mathrm{freq}}\right)}^2\right)\right]---\left(1\right)$

In the formula,

Be low-frequency vibration amplitude sequence

K data, k=1,2 ... l, l=98, and regulation is worked as ${\left(A_k^{\mathrm{Freq}}\right)}^2=0$ The time, $\mathrm{Ln}\left({\left(A_k^{\mathrm{Freq}}\right)}^2\right)=0.$

Low-frequency vibration amplitude sequence according to current time

Calculate the variation kurtosis parameter κ of current time low-frequency vibration amplitude sequence ^AAdopt formula

${\mathrm{\&kappa;}}^A=1/l\underset{k=1}{\overset{l}{\mathrm{\&Sigma;}}}{(A_k^{\mathrm{freq}}-{\mathrm{\&mu;}}^A)}^4/{\left({\mathrm{\&sigma;}}^A\right)}^4---\left(2\right)$

In the formula,

Be low-frequency vibration amplitude sequence K data, k=1,2 ... l, l=98, μ ^AIt is low-frequency vibration amplitude sequence

Average, promptly

σ ^AIt is low-frequency vibration amplitude sequence

Standard deviation, promptly ${\mathrm{\&sigma;}}^A=\sqrt{1/l\underset{k=1}{\overset{l}{\mathrm{\&Sigma;}}}{(A_k^{\mathrm{Freq}}-{\mathrm{\&mu;}}^A)}^2}.$

Every at a distance from second time stepping length t ₂=0.1 second, calculate the maximal value A of current time low-frequency vibration amplitude sequence ^FmaxBehind the entropy E of current time low-frequency vibration amplitude sequence, all to store, so that use in the subsequent calculations process.

Judge whether to satisfy condition simultaneously:

1) the variation kurtosis parameter κ of current time low-frequency vibration amplitude sequence ^AGreater than the first setting value D ₁, i.e. κ ^A>D ₁=30;

2) current time and initial moment T _SDifference more than or equal to the second setting value D ₂With the first stepping length t ₁Product, i.e. T _C-T _S>=D ₂* t ₁=10 * 30=300, T _CBe current time.

If satisfy condition 1 simultaneously) and 2), then with current time T _CAs stopping T constantly _N, with the current time T of storage _CPower of the assembling unit data P ^UAs stopping T constantly _NPower of the assembling unit data

In the present embodiment, suppose to work as T _CIn the time of=300 seconds, κ ^A>30, and T _C-T _S=300-0>=300 then stop T constantly _N=300 seconds, stop T constantly _N=300 seconds power of the assembling unit data do

Step 104: choose and stop T constantly _NIn two moment before, be designated as first T constantly formerly respectively ₁With second T constantly formerly ₂, and satisfy

In the present embodiment, choose first T constantly formerly ₁=0 second, second T constantly formerly ₂=290 seconds. $\frac{T_N-T_1}{t_1}=\frac{300-0}{3}=100,$ $\frac{T_N-T_2}{t_2}=\frac{300-290}{0.1}=100,$ Satisfy condition $\frac{T_N-T_1}{t_1}=\frac{T_N-T_2}{t_2}.$

Step 105: will be from first T constantly formerly ₁Rose in=0 second to stopping T constantly _NEnded in=300 seconds, every at a distance from very first time stepping length t ₁The power of the assembling unit data P of storage in=3 seconds ^ULine up power of the assembling unit data sequence according to the sequencing of storage time

Because whenever at a distance from 3 seconds storage one secondary data, therefore 100 power of the assembling unit data, power of the assembling unit data sequence were stored up in coexistence from 0 second to 300 seconds

The data number be 100, i.e. i=1,2 ... 100.

Will be from second T constantly formerly ₂Rose in=290 seconds to stopping T constantly _NEnded in=300 seconds, every at a distance from the second stepping length t ₂The maximal value A of the low-frequency vibration amplitude sequence of storage in=0.1 second ^FmaxLine up low-frequency vibration amplitude maximal value sequence according to the sequencing of storage time

Every at a distance from the second stepping length t ₂The entropy E of the low-frequency vibration amplitude sequence of storage in=0.1 second lines up the entropy sequence { E of low-frequency vibration amplitude sequence according to the sequencing of storage time _j,

Because whenever at a distance from 0.1 second storage one secondary data, therefore the entropy of 100 low-frequency vibration amplitude maximal values and 100 low-frequency vibration amplitude sequences was stored up in coexistence, so low-frequency vibration amplitude maximal value sequence from 290 seconds to 300 seconds

Entropy sequence { E with low-frequency vibration amplitude sequence _jThe data number all be 100, i.e. j=1,2 ... 100.

Step 106: computer set power parameter and low-frequency vibration parameter comprise:

1) computer set power data sequence increases progressively the trend parameter I ^P

Computer set power data sequence increase progressively the trend parameter I ^PAdopt formula

I ^P=S ^P/[1/2×n×(n-1)] （3）

In the formula, n is the data number of power of the assembling unit data sequence, in the present embodiment, and n=100.S ^PIt is the Ser.No. of power of the assembling unit data sequence.Ser.No. is meant the right sum of order in the data sequence.Order is to being meant that the front and back position of a logarithm is identical with size order in a data sequence, and promptly the number of front is less than the number of back.

2) maximal value

of calculating low-frequency vibration amplitude maximal value sequence

The maximal value of low-frequency vibration amplitude maximal value sequence is utilized in formula

present embodiment, n=100.

3) the variation degree of bias parameter S of the entropy sequence of calculating low-frequency vibration amplitude sequence ^E

Calculate the variation degree of bias parameter S of the entropy sequence of low-frequency vibration amplitude sequence ^EAdopt formula

$S^E=1/n\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{(E_j-{\mathrm{\&mu;}}^E)}^3/{\left({\mathrm{\&sigma;}}^E\right)}^3---\left(4\right)$

In the formula, E _jBe j data of the entropy sequence of low-frequency vibration amplitude sequence, j=1,2 ..., n, μ ^EBe the average of the entropy sequence of low-frequency vibration amplitude sequence, promptly

σ ^EBe the standard deviation of the entropy sequence of low-frequency vibration amplitude sequence, promptly

N is the data number of the entropy sequence of low-frequency vibration amplitude sequence.In the present embodiment, n=100.

4) calculate the entropy sequence of low-frequency vibration amplitude sequence and the Kendall related coefficient τ of power of the assembling unit data sequence.

Calculate the entropy sequence of low-frequency vibration amplitude sequence and the Kendall related coefficient τ of power of the assembling unit data sequence and adopt formula

τ=(n _s-n _d)/[1/2×(n ²-1)] （5）

In the formula, n is the entropy sequence of low-frequency vibration amplitude sequence or the data number of power of the assembling unit data sequence, in the present embodiment, and n=100.n _sBe in the entropy sequence and power of the assembling unit data sequence of low-frequency vibration amplitude sequence, the right sum of sequence that mediation is arranged, n _dBe in the entropy sequence and power of the assembling unit data sequence of low-frequency vibration amplitude sequence, the sum that the sequence of anharmonic arrangement is right.

The sequence that mediation is arranged is to being meant, for 2 sequence { A _kAnd { B _kIn any 2 sequences to (A _m, B _m) and (A _n, B _n), q>=m, n>=1, q is sequence { A _kAnd sequence { B _kThe data number, if satisfy A simultaneously _m>A _nAnd B _m>B _nPerhaps satisfy A simultaneously _m<a _nAnd B _m<b _n, then 2 sequences are to (A _m, B _m) and (A _n, B _n) be that the sequence of be in harmonious proportion arranging is right.The sequence of anharmonic arrangement is to being meant, for 2 sequence { A _kAnd { B _kIn any 2 sequences to (A _m, B _m) and (A _n, B _n), q>=m, n>=1, q is sequence { A _kAnd sequence { B _kThe data number, if satisfy A simultaneously _m>A _nAnd B _m<b _nPerhaps satisfy A simultaneously _m<a _nAnd B _m>B _n, then 2 sequences are to (A _m, B _m) and (A _n, B _n) be that the sequence of anharmonic arrangement is right.For 2 sequence { A _kAnd { B _kIn any 2 sequences to (A _m, B _m) and (A _n, B _n), Q>=m, n>=1, Q is sequence { A _kOr sequence { B _kThe data number, if satisfy A _m=A _nPerhaps B _m=B _n, then 2 sequences are to (A _m, B _m) and (A _n, B _n) neither to be in harmonious proportion collating sequence right, neither anharmonic collating sequence right.

In the present invention, for the entropy sequence { E of low-frequency vibration amplitude sequence<sub >i</sub>And power of the assembling unit data sequence<img file="BDA00001636897400121.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="95" />In any 2 sequences right<img file="BDA00001636897400122.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="172" />With<img file="BDA00001636897400123.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="193" />100>=m, n>=1 is if satisfy E simultaneously<sub >m</sub>>E<sub >n</sub>And<img file="BDA00001636897400124.GIF" he="57" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="169" />Perhaps satisfy E simultaneously<sub >m</sub><e<sub >n</sub>And<img file="BDA00001636897400125.GIF" he="57" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="198" />Then these 2 sequences are right<img file="BDA00001636897400126.GIF" he="57" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="173" />With<img file="BDA00001636897400127.GIF" he="57" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="166" />The sequence that is the mediation arrangement is right.Entropy sequence { E for low-frequency vibration amplitude sequence<sub >i</sub>And power of the assembling unit data sequence<img file="BDA00001636897400128.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="95" />In any 2 sequences right<img file="BDA00001636897400129.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="173" />With<img file="BDA000016368974001210.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="193" />100>=m, n>=1 is if satisfy E simultaneously<sub >m</sub>>E<sub >n</sub>And<img file="BDA000016368974001211.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="169" />Perhaps satisfy E simultaneously<sub >m</sub><e<sub >n</sub>And<img file="BDA000016368974001212.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="198" />Then these 2 sequences are right<img file="BDA000016368974001213.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="172" />With<img file="BDA000016368974001214.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="166" />The sequence that is anharmonic arrangement is right.Entropy sequence { E for low-frequency vibration amplitude sequence<sub >i</sub>And power of the assembling unit data sequence<img file="BDA000016368974001215.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="95" />In any 2 sequences right<img file="BDA000016368974001216.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="172" />With<img file="BDA000016368974001217.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="192" />100>=m, n>=1 is if satisfy E<sub >m</sub>=E<sub >n</sub>Perhaps<img file="BDA000016368974001218.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="198" />Then these 2 sequences are right<img file="BDA000016368974001219.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="173" />With<img file="BDA000016368974001220.GIF" he="56" img-content="drawing" img-format="GIF" inline="yes" orientation="portrait" wi="166" />Neither it is the collating sequence that is in harmonious proportion is right, neither anharmonic collating sequence right.

Step 107: according to stopping the power of the assembling unit data of TN constantly Power of the assembling unit data sequence increase progressively the trend parameter I ^P, low-frequency vibration amplitude maximal value sequence maximal value

The variation degree of bias parameter S of the entropy sequence of low-frequency vibration amplitude sequence ^EKendall related coefficient τ with the entropy sequence and the power of the assembling unit data sequence of low-frequency vibration amplitude sequence judges whether high pressure rotor the steam flow excitation fault takes place.

Set the 3rd setting value D respectively ₃=150MW (megawatt), the 4th setting value D ₄The=0.8, the 5th setting value D ₅=30 μ m (micron), the 6th setting value D ₆The=1.5 and the 7th setting value D ₇=0.7.Above-mentioned setting value is used for assisting to judge whether high pressure rotor the steam flow excitation fault takes place that each setting value is confirmed according to the high pressure rotor running requirements and the standard of Turbo-generator Set.

Satisfy following 5 conditions when simultaneously:

1) stops T constantly _NPower of the assembling unit data

Greater than the 3rd setting value, promptly

2) power of the assembling unit data sequence increases progressively the trend parameter I ^PGreater than the 4th setting value, i.e. I ^P>D ₄=0.8;

The maximal value of 3) low-frequency vibration amplitude maximal value sequence

is greater than the 5th setting value, promptly

$A_{\max }^{f\max }>D_5=30\mathrm{\&mu;m};$

The variation degree of bias parameter of the entropy sequence of 4) low-frequency vibration amplitude sequence is greater than S ^EGreater than the 6th setting value, i.e. S ^E>D ₆=1.5;

The Kendall related coefficient τ of the entropy sequence of 5) low-frequency vibration amplitude sequence and power of the assembling unit data sequence is greater than the 7th setting value, i.e. τ>D ₇=0.7;

Then judge high pressure rotor generation steam flow excitation fault; Otherwise, judge that the steam flow excitation fault does not take place high pressure rotor.

Turbo-generator Set steam flow excitation On-line Fault method of discrimination provided by the invention; To unit operation rotor axle vibrate relatively, data such as the power of the assembling unit carry out real-time automatic on-line monitoring, analyze and differentiate; Judge whether high pressure rotor the steam flow excitation fault takes place, improve the efficient and the accuracy of high pressure rotor steam flow excitation fault analysis and diagnosis work.

The above; Be merely the preferable embodiment of the present invention, but protection scope of the present invention is not limited thereto, any technician who is familiar with the present technique field is in the technical scope that the present invention discloses; The variation that can expect easily or replacement all should be encompassed within protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with the protection domain of claim.

Claims (7)

1. Turbo-generator Set steam flow excitation On-line Fault method of discrimination is characterized in that said method comprises:

Step 1: set initial moment T _S, very first time stepping length t ₁, second time stepping length t ₂, the first setting value D ₁With the second setting value D ₂

Step 2: gather the relative vibration data of axle of Turbo-generator Set high pressure rotor one side radial journal bearing, the tach signal of rotor, the key signal and the power of the assembling unit data of rotor in real time;

Step 3: from initial moment T _SBeginning, every at a distance from very first time stepping length t ₁The power of the assembling unit data P of storage current time ^UEvery at a distance from second time stepping length t ₂, the maximal value A of calculating current time low-frequency vibration amplitude sequence ^Fmax, the entropy E of current time low-frequency vibration amplitude sequence and the variation kurtosis parameter κ of current time low-frequency vibration amplitude sequence ^A, and the maximal value A of storage current time low-frequency vibration amplitude sequence ^FmaxEntropy E with current time low-frequency vibration amplitude sequence;

As the variation kurtosis parameter κ that satisfies current time low-frequency vibration amplitude sequence ^AGreater than the first setting value D ₁And current time and initial moment T _SDifference more than or equal to the second setting value D ₂With the first stepping length t ₁Product the time, with current time as stopping T constantly _N, with the power of the assembling unit data P of the current time of storing ^UAs stopping T constantly _NPower of the assembling unit data

Execution in step 4;

Step 4: choose and stop T constantly _NIn two moment before, be designated as first T constantly formerly respectively ₁With second T constantly formerly ₂, and satisfy

Step 5: will be from first T constantly formerly ₁Rise to stopping T constantly _NEnd, every at a distance from very first time stepping length t ₁The power of the assembling unit data P of storage ^ULine up power of the assembling unit data sequence according to the sequencing of storage time

Will be from second T constantly formerly ₂Rise to stopping T constantly _NEnd, every at a distance from the second stepping length t ₂The maximal value A of the low-frequency vibration amplitude sequence of storage ^FmaxLine up low-frequency vibration amplitude maximal value sequence according to the sequencing of storage time

Every at a distance from the second stepping length t ₂The entropy E of the low-frequency vibration amplitude sequence of storage lines up the entropy sequence { E of low-frequency vibration amplitude sequence according to the sequencing of storage time _j,

Step 6: computer set power parameter and low-frequency vibration parameter comprise:

1) computer set power data sequence increases progressively the trend parameter I ^P

2) maximal value

of calculating low-frequency vibration amplitude maximal value sequence

3) the variation degree of bias parameter S of the entropy sequence of calculating low-frequency vibration amplitude sequence ^E

4) calculate the entropy sequence of low-frequency vibration amplitude sequence and the Kendall related coefficient τ of power of the assembling unit data sequence;

Step 7: according to stopping T constantly _NPower of the assembling unit data

Power of the assembling unit data sequence increase progressively the trend parameter I ^P, low-frequency vibration amplitude maximal value sequence maximal value

The variation degree of bias parameter S of the entropy sequence of low-frequency vibration amplitude sequence ^EKendall related coefficient τ with the entropy sequence and the power of the assembling unit data sequence of low-frequency vibration amplitude sequence judges whether high pressure rotor the steam flow excitation fault takes place.

2. Turbo-generator Set steam flow excitation On-line Fault method of discrimination according to claim 1 is characterized in that the entropy E of said calculating current time low-frequency vibration amplitude sequence adopts formula

$E=\underset{k=1}{\overset{l}{\mathrm{\&Sigma;}}}\left[{\left(A_k^{\mathrm{freq}}\right)}^2\ln \left({\left(A_k^{\mathrm{freq}}\right)}^2\right)\right],$

Wherein,

is k data of low-frequency vibration amplitude sequence; K=1; 2; ... l; L is the data number of predefined low-frequency vibration amplitude sequence; And when stipulating as

,

3. Turbo-generator Set steam flow excitation On-line Fault method of discrimination according to claim 1 is characterized in that the variation kurtosis parameter κ of said current time low-frequency vibration amplitude sequence ^AAdopt formula

${\mathrm{\&kappa;}}^A=1/l\underset{k=1}{\overset{l}{\mathrm{\&Sigma;}}}{(A_k^{\mathrm{freq}}-{\mathrm{\&mu;}}^A)}^4/{\left({\mathrm{\&sigma;}}^A\right)}^4,$

Wherein,

Be k data of low-frequency vibration amplitude sequence, k=1,2 ... l, l are the data number of predefined low-frequency vibration amplitude sequence, μ ^ABe the average of low-frequency vibration amplitude sequence, promptly

σ ^ABe the standard deviation of low-frequency vibration amplitude sequence, promptly

4. Turbo-generator Set steam flow excitation On-line Fault method of discrimination according to claim 1, what it is characterized in that said computer set power data sequence increases progressively the trend parameter I ^PAdopt formula

I ^P=S ^P/[1/2×n×(n-1)]，

Wherein, n is the data number of power of the assembling unit data sequence, S ^PIt is the Ser.No. of power of the assembling unit data sequence; Said Ser.No. is meant the right sum of order in the data sequence; Said order is to being meant that the front and back position of a logarithm is identical with size order in a data sequence, and promptly the number of front is less than the number of back.

5. Turbo-generator Set steam flow excitation On-line Fault method of discrimination according to claim 1 is characterized in that the variation degree of bias parameter S of the entropy sequence of said calculating low-frequency vibration amplitude sequence ^EAdopt formula

$S^E=1/n\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{(E_j-{\mathrm{\&mu;}}^E)}^3/{\left({\mathrm{\&sigma;}}^E\right)}^3,$

Wherein, E _jBe j data of the entropy sequence of low-frequency vibration amplitude sequence, j=1,2 ..., n, μ ^EBe the average of the entropy sequence of low-frequency vibration amplitude sequence, promptly

σ ^EBe the standard deviation of the entropy sequence of low-frequency vibration amplitude sequence, promptly

N is the data number of the entropy sequence of low-frequency vibration amplitude sequence.

6. Turbo-generator Set steam flow excitation On-line Fault method of discrimination according to claim 1 is characterized in that the entropy sequence of said calculating low-frequency vibration amplitude sequence and the Kendall related coefficient τ of power of the assembling unit data sequence adopt formula

τ=(n _s-n _d)/[1/2×(n ²-1)]，

Wherein, n is the entropy sequence of low-frequency vibration amplitude sequence or the data number of power of the assembling unit data sequence, n _sBe in the entropy sequence and power of the assembling unit data sequence of low-frequency vibration amplitude sequence, the right sum of sequence that mediation is arranged, n _dBe in the entropy sequence and power of the assembling unit data sequence of low-frequency vibration amplitude sequence, the sum that the sequence of anharmonic arrangement is right;

The sequence that said mediation is arranged is to being meant, for 2 sequence { A _kAnd { B _kIn any 2 sequences to (A _m, B _m) and (A _n, B _n), q>=m, n>=1, q is sequence { A _kAnd sequence { B _kThe data number, if satisfy A simultaneously _m>A _nAnd B _m>B _nPerhaps satisfy A simultaneously _m<a _nAnd B _m<b _n, then 2 sequences are to (A _m, B _m) and (A _n, B _n) be that the sequence of be in harmonious proportion arranging is right;

The sequence of said anharmonic arrangement is to being meant, for 2 sequence { A _kAnd { B _kIn any 2 sequences to (A _m, B _m) and (A _n, B _n), q>=m, n>=1, q is sequence { A _kAnd sequence { B _kThe data number, if satisfy A simultaneously _m>A _nAnd B _m<b _nPerhaps satisfy A simultaneously _m<a _nAnd B _m>B _n, then 2 sequences are to (A _m, B _m) and (A _n, B _n) be that the sequence of anharmonic arrangement is right.

7. Turbo-generator Set steam flow excitation On-line Fault method of discrimination according to claim 1 is characterized in that said step 7 specifically is, satisfies following 5 conditions when simultaneously:

1) power of the assembling unit data

that stop moment TN are greater than the 3rd setting value;

2) power of the assembling unit data sequence increases progressively the trend parameter I ^PGreater than the 4th setting value;

The maximal value of 3) low-frequency vibration amplitude maximal value sequence

is greater than the 5th setting value;

The variation degree of bias parameter of the entropy sequence of 4) low-frequency vibration amplitude sequence is greater than S ^EGreater than the 6th setting value;

The Kendall related coefficient τ of the entropy sequence of 5) low-frequency vibration amplitude sequence and power of the assembling unit data sequence is greater than the 7th setting value;

Then judge high pressure rotor generation steam flow excitation fault; Otherwise, judge that the steam flow excitation fault does not take place high pressure rotor.