Summary of the invention
The invention provides a kind of isolating switch of on-line monitoring that is applicable to and just divided firm chalaza determination methods and device thereof.The present invention is simple and practical, computational accuracy is higher relatively; Its installation cost is lower, reliable and stable, can move continuously, can satisfy the requirement that on-line monitoring detects about the firm chalaza of firm branch and other mechanical parameters.
The present invention adopts following technical scheme:
A kind of isolating switch of the present invention has just divided the online determination methods of firm chalaza; Comprise: online determination methods of the firm chalaza of isolating switch and the firm online determination methods of branch; It is characterized in that; Utilize the online determination methods of the firm chalaza of isolating switch respectively the firm chalaza of A, B, C three-phase enforcement isolating switch to be carried out online judgement, utilize the online determination methods of firm branch respectively the firm branch of A, B, C three-phase enforcement isolating switch to be carried out online judgement
The online determination methods of the firm chalaza of said isolating switch may further comprise the steps:
Step 1 is according to sampling time of sampled point order, utilizes displacement transducer to obtain the shift simulation amount disA (i) of the moving contact of breaker making process of constant duration T successively, i=1,2,3 ...; Time interval T=1/f; F is a sampling rate, and i is the sequence number of sampled point and is positive integer
Step 2 placement algorithm parameter: select the breaker operation mechanism bottom panel as the measuring basis face; And with the distance between measuring basis face to the gate-dividing state lower moving contact as the actual least displacement value of moving contact Amin; With the distance between measuring basis face to the "on" position lower moving contact as the actual maximum displacement value Amax of moving contact; Obtain actual maximum displacement value Amax of moving contact and the actual least displacement value of moving contact Amin through measurement again
The moving contact displacement signal pre-service that step 3 pair step 1 is gathered: adopt the wavelet noise algorithm to carry out the one dimension noise reduction, eliminate noise, data are disA1 (i) after the de-noising, i=1,2,3 ....Utilize the limit filtration means to remove burr, obtaining pretreated discrete displacement data is disA2 (i), i=1,2,3 ...,
The calculating of step 4 actual discrete displacement disAs (i), total kilometres strA, close a floodgate initial moment nhA1 and discrete closing speed yA (i):
Step 4.1 actual discrete displacement disAs (i) calculates:
At first calculate total kilometres strA, total kilometres strA=Amax-Amin;
Next calculates linear scale factor KK; Said linear scale factor KK represent the actual displacement amount and the signal voltage value gathered between linear scaling; Computing formula is: KK=strA/ (av2-av1); Wherein av2 representes the displacement transducer sample voltage value of "on" position, and av1 representes the displacement transducer sample voltage value of gate-dividing state;
Obtain actual discrete displacement disAs (i) at last, the actual discrete displacement is:
disAs(i)=(disA2(i)-av1)×KK
The initial moment nhA1 of step 4.2 search: search for disAs (i) backward one by one, and with first greater than making process total kilometres strA 10% the point in the correspondence moment be nhA1,
The pairing discrete closing speed yA of step 4.3 actual discrete displacement disAs (i) (i) obtains through 3 formula of numerical differentiation; DisAs (l), disAs (l+1), disAs (l+2) are carried out segmentation as one group to be calculated; L=1,4,7 ..., 3n+1; N is a natural number, and its algorithm is following:
L=1,4,7 ..., 3n+1, n is a natural number,
Obtain the discrete closing speed yA (i) of making process, i=1,2,3 ...,
The firm chalaza of step 5 making process judges that it is specifically following step by step:
Step 5.1 is starting point with nhA1, search speed nhA2 at zero point backward, and its basis for estimation is:
DisAs (k+1)-disAs (k)<0 and disAs (k)-disAs (k-1)>0
Wherein, k successively value be nhA1, nhA1+1, nhA1+2 ..., be the nhA2 point until the k that occurs satisfying condition and with first k point that satisfies condition,
Step 5.2 is starting point with nhA2, searches for acceleration nhA3 at zero point forward, and its basis for estimation is:
YA (j+1)-yA (j)<0 and yA (j)-yA (j-1)>0
Wherein, j successively value be nhA2, nhA2-1, nhA2-2 ..., be the nhA3 point until the j that occurs satisfying condition and with first j point that satisfies condition,
Step 6 as firm chalaza, is asked for isolating switch contacting travel ConA with the nhA3 point,
ConA=strA-disAs(nhA3);
The online determination methods of the firm branch of said isolating switch may further comprise the steps:
Step 7 is according to sampling time of sampled point order, utilizes displacement transducer to obtain the shift simulation amount dispA (i) of the moving contact of breaker separating brake process of constant duration T successively, i=1,2,3 ...; Time interval T=1/f; F is a sampling rate, and i is the sequence number of sampled point and is positive integer
The moving contact displacement signal pre-service that step 8 pair step 7 is gathered: adopt the wavelet noise algorithm to carry out the one dimension noise reduction; Eliminate noise, utilize the limit filtration means to remove burr, obtaining pretreated discrete displacement data is dispA2 (i); I=1,2,3 ...
Step 9 actual discrete displacement dispAs (i) calculates, dispAs (i)=(dispA2 (i)-av1) * KK,
Step 10 is searched for dispAs (i) in order backward according to the principle that divide-shut brake process contacting travel ConA equates, i=1,2,3 ..., and be nfA1 with first point greater than the ConA institute correspondence moment, i.e. the firm branch of separating brake process.
A kind of device that isolating switch has just divided the online determination methods of firm chalaza that is used to implement of the present invention; Comprise: the operating mechanism that is provided with isolating switch; Base plate with operating mechanism is the measuring basis face; And on operating mechanism, be provided with the three-phase displacement transducer, the three-phase displacement transducer is connected through " Z " shape web member with insulated tension pole in the isolating switch.
Compared with prior art, the present invention has following advantage:
The present invention utilizes acceleration to judge the firm branch of breaker closing zero point, have higher confidence level, and this method computation process is simple, and operand is less, is applicable to on-line monitoring.Utilize the three-phase displacement transducer to gather stroke signal, can calculate most mechanical parameters of isolating switch operation, do not need extra data channel, make that apparatus structure is simple, lower, the stable and reliable operation of cost, satisfy the requirement of on-line monitoring.Utilize wavelet noise means and threshold value limit filtration burr removing method, eliminated the interference that site environment noise and the mechanical vibration of isolating switch own are produced, the accuracy that has improved determination methods.
Embodiment
Below in conjunction with accompanying drawing the present invention is done explanation further.
In order to obtain the exact position that the isolating switch three-phase has just divided firm chalaza, the at first just model and the installation site of necessary choose reasonable sensor respectively.The present invention selects miniature displacement transducer KTS-25, is characterized in advantages of small volume, is easy to install; Simple in structure, power consumption is lower; Higher linearity and repeatable accuracy are arranged, can guarantee sampling precision; Adopt resistance sensor, response is very fast, has guaranteed the real-time of data, is specially adapted to 10-20KV vacuum circuit breaker on-line monitoring system.The present invention adopt three independently displacement transducer measure the displacement of three-phase divide-shut brake respectively, overcome and adopted photoelectric encoder can only measure the defective of single-phase displacement in the past.Sensor is installed in moving contact of breaker insulated tension pole below, and sensor activity parts pull bar is fixed on the insulating pull rod for breaker, and shell is fixed on the isolating switch pedestal.Installation direction makes sensor parallel with insulated tension pole, so the displacement of sensor can reflect the displacement of moving contact itself fully.Concrete installation site is as shown in Figure 1.
The sampling rate 25KHz that the present invention adopts, sampling length arrives 2s afterwards for 200ms before receiving deciliter signal.The resulting signal of three-phase displacement transducer through the signal condition filtering circuit, is translated into digital signal by the A/D change-over circuit.Transfer to primary processor ARM through auxiliary processor FPGA, ARM receives displacement data, its sampled data with other monitoring variables is packed according to self-defined message format, and be sent to the up-computer data processing unit through Ethernet interface.Host computer receives message, unpacks and set up data database storing.Form the data stream of real-time synchronous high-speed, can guarantee the validity and the degree of accuracy of data.Software algorithm realizes through MATLAB software, generates com component through the compiling of MATLAB Builder for COM compiler, calls com component during isolating switch on-line monitoring running software and accomplishes data computation, and the form of result of calculation with data and figure shown.
Be illustrated in figure 2 as a kind of isolating switch of on-line monitoring that is applicable to and just divided firm chalaza determination methods and device software algorithm flow chart.Adopt identical algorithm flow following to A, B, C three-phase:
(1) read host computer storage three-phase making process displacement data disA (i) i=1,2,3 ..., and predefined configuration data comprises SF, the actual maximum displacement value Amax of moving contact and the actual least displacement value of moving contact Amin.
(2) raw data pre-service.Data unavoidably can receive some extraneous interference in collection, transmission, storing process, the interference of some sampled point even quite serious has produced very big influence to the computational accuracy of algorithm.In order to get rid of these interference, need take some de-noisings and means of filtering.Small echo has obtained using widely in signal noise silencing as the important means in the signal Processing.Adopt one dimension wavelet de-noising mode wden, its parameter is chosen as: select the sym5 small echo; The threshold treatment mode adopts the soft-threshold mode, and parameter is Sorh=s; Threshold value selects to adopt maximum minimum principle to produce threshold value, and parameter is tptr=minimaxi; Noise level according to each layer wavelet decomposition estimates to adjust Scal=mln; 5 layers of wavelet decomposition and reconstruct.The displacement discrete data is disA1 (i) after the de-noising, and its displacement curve is level and smooth, more near actual.In order to eliminate the formed large-scale sudden change of some bigger interference; Adopt deburring limit filtration algorithm following: the numerical value of more adjacent two sampled point disA1 (i) and disA1 (i+1); If the difference between them has surpassed the threshold value 0.1 that sets; Then think mutation disturbance has taken place, need give up the numerical value of this point, come this disA1 of reconstruct (i+1)=(disA1 (i)+disA1 (i+2))/2 with the linear interpolation of disA1 (i+2) with disA1 (i).Discrete data is disA2 (i) after the pre-service.
(3) combined floodgate actual displacement disAs (i), the initial moment nhA1 that closes a floodgate, the discrete speed yA (i) that closes a floodgate calculate.Because the simulating signal of displacement transducer collection is a magnitude of voltage, so must utilize configuration data to be translated into actual shift value.Algorithm is following: at first calculate linear scale factor KK, the linear scaling between expression actual displacement amount and the signal voltage value gathered.KK=strA/ (av2-av1), wherein av2 representes the displacement transducer sample voltage value of A phase "on" position, av1 representes the displacement transducer sample voltage value of A phase gate-dividing state, total kilometres strA=Amax-Amin.Proportion of utilization COEFFICIENT K K is converted into actual discrete displacement disAs (i) with all sampled datas again, disAs (i)=(disA2 (i)-av1) * KK.Obtain actual discrete displacement disAs (i) thus.In the calculating of mechanical parameter, usually can be with 10% total kilometres as starting point.Therefore, need ask for 10% travel point nhA1.Can be through searching for the actual discrete displacement backward one by one, like disAs (i), and with first greater than making process total kilometres strA 10% correspondence point constantly be nhA1.This algorithm is based on breaker closing process acceleration criterion at zero point, therefore must obtain the discrete velocity amplitude of each sample point.The pairing discrete closing speed yA of actual discrete displacement disAs (i) (i) obtains through 3 formula of numerical differentiation; DisAs (l), disAs (l+1), disAs (l+2) are carried out segmentation as one group to be calculated; L=1,4,7 ..., 3n+1, n is a natural number, its algorithm is following:
Wherein l=1,4,7 ..., 3n+1, n is a natural number, obtains the discrete closing speed yA (i) of making process.
(4) just chalaza is judged.Tradition utilize three phase shifts to firm chalaza judge be through speed zero point as foundation, obtained certain effect.Yet, exist certain error as firm chalaza zero point with moving contact of breaker speed.This is because dynamic/static contact itself has certain elasticity, can produce certain excess of stroke in the making process.In firm chalaza constantly, because the reacting force of static contact makes moving contact speed begin to descend, this brief acceleration reduces to zero, but velocity amplitude is non-vanishing.Therefore, acceleration zero point as criterion relative velocity zero point of just chalaza, to have higher precision relatively, more near actual conditions.Its algorithm is following: with nhA1 as starting point; Utilize formula disAs (k+1)-disAs (k)<0 and disAs (k)-disAs (k-1)>0 to search for disAs (i) backward as basic criterion; K successively value be nhA1, nhA1+1, nhA1+2 ..., be the nhA2 point until the k that occurs satisfying condition and with first k point that satisfies condition.Because backward acceleration is bigger, speed will reduce to zero very soon, so speed zero point is very approaching and at acceleration rear at zero point with acceleration zero point.So with speed nhA2 at zero point as starting point, reverse search acceleration zero point in its front certain limit.Judgment criterion: yA (j+1)-yA (j)<0 and yA (j)-yA (j-1)>0; J successively value be nhA2, nhA2-1, nhA2-2 ... Be the nhA3 point until the j that occurs satisfying condition and with first j point that satisfies condition, with this as the firm chalaza of breaker closing process.Its specific algorithm flow process is as shown in Figure 3.
(5) the breaker closing contacting travel is defined as the firm chalaza of isolating switch to the moving contact of breaker stop motion displacement of moving contact constantly, can easily obtain contacting travel ConA=strA-disAs (nhA3) according to firm chalaza nhA3.
(6) the firm branch of separating brake process is calculated.At first read the shift simulation amount dispA (i) of moving contact separating brake process, carry out de-noising and deburring pre-service, and image data dispA (i) is converted into actual displacement amount dispAs (i) through linear scale factor KK with above-mentioned identical method.According to the principle that divide-shut brake process contacting travel ConA equates, search for dispAs (i) in order backward, and be nfA1, be i.e. the firm branch of separating brake process with first point greater than the ConA institute correspondence moment.
Embodiment 1
A kind of isolating switch has just divided the online determination methods of firm chalaza; Comprise: online determination methods of the firm chalaza of isolating switch and the firm online determination methods of branch; It is characterized in that; Utilize the online determination methods of the firm chalaza of isolating switch respectively the firm chalaza of A, B, C three-phase enforcement isolating switch to be carried out online judgement, utilize the online determination methods of firm branch respectively the firm branch of A, B, C three-phase enforcement isolating switch to be carried out online judgement
The online determination methods of the firm chalaza of said isolating switch may further comprise the steps:
Step 1 is according to sampling time of sampled point order, utilizes displacement transducer to obtain the shift simulation amount disA (i) of the moving contact of breaker making process of constant duration T successively, i=1,2,3 ...; Time interval T=1/f; F is a sampling rate, and i is the sequence number of sampled point and is positive integer
Step 2 placement algorithm parameter: select the breaker operation mechanism bottom panel as the measuring basis face; And with the distance between measuring basis face to the gate-dividing state lower moving contact as the actual least displacement value of moving contact Amin; With the distance between measuring basis face to the "on" position lower moving contact as the actual maximum displacement value Amax of moving contact; Obtain actual maximum displacement value Amax of moving contact and the actual least displacement value of moving contact Amin through measurement again
The moving contact displacement signal pre-service that step 3 pair step 1 is gathered: adopt the wavelet noise algorithm to carry out the one dimension noise reduction, eliminate noise, data are disA1 (i) after the de-noising, i=1,2,3 ....Utilize the limit filtration means to remove burr, obtaining pretreated discrete displacement data is disA2 (i), i=1,2,3 ...,
The calculating of step 4 actual discrete displacement disAs (i), total kilometres strA, close a floodgate initial moment nhA1 and discrete closing speed yA (i):
Step 4.1 actual discrete displacement disAs (i) calculates:
At first calculate total kilometres strA, total kilometres strA=Amax-Amin;
Next calculates linear scale factor KK; Said linear scale factor KK represent the actual displacement amount and the signal voltage value gathered between linear scaling; Computing formula is: KK=strA/ (av2-av1); Wherein av2 representes the displacement transducer sample voltage value of "on" position, and av1 representes the displacement transducer sample voltage value of gate-dividing state;
Obtain actual discrete displacement disAs (i) at last, the actual discrete displacement is:
disAs(i)=(disA2(i)-av1)×KK
The initial moment nhA1 of step 4.2 search: search for disAs (i) backward one by one, and with first greater than making process total kilometres strA 10% the point in the correspondence moment be nhA1,
The pairing discrete closing speed yA of step 4.3 actual discrete displacement disAs (i) (i) obtains through 3 formula of numerical differentiation; DisAs (l), disAs (l+1), disAs (l+2) are carried out segmentation as one group to be calculated; L=1,4,7 ..., 3n+1; N is a natural number, and its algorithm is following:
L=1,4,7 ..., 3n+1, n is a natural number,
Obtain the discrete closing speed yA (i) of making process, i=1,2,3 ...,
The firm chalaza of step 5 making process judges that it is specifically following step by step:
Step 5.1 is starting point with nhA1, search speed nhA2 at zero point backward, and its basis for estimation is:
DisAs (k+1)-disAs (k)<0 and disAs (k)-disAs (k-1)>0
Wherein, k successively value be nhA1, nhA1+1, nhA1+2 ..., be the nhA2 point until the k that occurs satisfying condition and with first k point that satisfies condition,
Step 5.2 is starting point with nhA2, searches for acceleration nhA3 at zero point forward, and its basis for estimation is:
YA (j+1)-yA (j)<0 and yA (j)-yA (j-1)>0
Wherein, j successively value be nhA2, nhA2-1, nhA2-2 ..., be the nhA3 point until the j that occurs satisfying condition and with first j point that satisfies condition,
Step 6 as firm chalaza, is asked for isolating switch contacting travel ConA with the nhA3 point,
ConA=strA-disAs(nhA3);
The online determination methods of the firm branch of said isolating switch may further comprise the steps:
Step 7 is according to sampling time of sampled point order, utilizes displacement transducer to obtain the shift simulation amount dispA (i) of the moving contact of breaker separating brake process of constant duration T successively, i=1,2,3 ...; Time interval T=1/f; F is a sampling rate, and i is the sequence number of sampled point and is positive integer
The moving contact displacement signal pre-service that step 8 pair step 7 is gathered: adopt the wavelet noise algorithm to carry out the one dimension noise reduction; Eliminate noise, utilize the limit filtration means to remove burr, obtaining pretreated discrete displacement data is dispA2 (i); I=1,2,3 ...
Step 9 actual discrete displacement dispAs (i) calculates, dispAs (i)=(dispA2 (i)-av1) * KK,
Step 10 is searched for dispAs (i) in order backward according to the principle that divide-shut brake process contacting travel ConA equates, i=1,2,3 ..., and be nfA1 with first point greater than the ConA institute correspondence moment, i.e. the firm branch of separating brake process.
Said wavelet noise algorithm is to utilize Matlab that the moving contact displacement signal is carried out the one dimension wavelet de-noising, and its mode is wden, and parameter is chosen as: select the sym5 small echo; The threshold treatment mode adopts the soft-threshold mode, and parameter is Sorh=s; Threshold value selects to adopt maximum minimum principle to produce threshold value, and parameter is tptr=minimaxi; Noise level according to each layer wavelet decomposition estimates to adjust Scal=mln; 5 layers of wavelet decomposition and reconstruct.
3, isolating switch according to claim 1 has just divided the online determination methods of firm chalaza, it is characterized in that said deburring algorithm is following: the numerical value of more adjacent two sampled point dis A1 (i) and disAl (i+1); I=1,2,3 ..., if the difference between them has surpassed the threshold value 0.1 that sets, then think mutation disturbance has taken place; The numerical value that need give up this point; Threshold value can be set through manual, comes this idsA1 of reconstruct (i+1) with disAl (i) with the linear interpolation of disAl (i+2)
disAl(i+1)=(disA1(i)+disA1(i+2))/2。
Embodiment 2
A kind ofly be used to implement the device that isolating switch has just divided the online determination methods of firm chalaza; Comprise: the operating mechanism 1 that is provided with isolating switch 11; Base plate with operating mechanism 1 is a measuring basis face 4; And on operating mechanism 1, be provided with three-phase displacement transducer 2, three-phase displacement transducer 2 is connected through " Z " shape web member 3 with insulated tension pole 111 in the isolating switch 11.
Embodiment 3
Use three models to gather vacuum circuit breaker three phase displacement signals of model respectively in the practical implementation as VMB5-12 as the micro linear displacement transducer of KTS-25; The letter that the isolating switch on-line monitoring system that utilizes independent development to develop will be gathered is changed through A/D, is sent to host computer through Ethernet interface and accomplishes data processing, graphic presentation and data storage function.Be illustrated in figure 4 as the original signal figure and the preprocessed signal figure of vacuum circuit breaker three phase shifts.As can be seen from the figure, contain a large amount of undesired signals in the raw bits shifting signal, brought difficulty for follow-up calculating.After passing threshold limit filtration and the wavelet noise, its displacement curve is level and smooth relatively, and characteristic quantity is comparatively obvious, the requirement that algorithmic is calculated.
Combined floodgate and the separating brake process mechanism characterisitic parameter that is calculated for the isolating switch on-line monitoring software of independent development as shown in table 1.Through comparing with VMB5-12 type vacuum circuit breaker delivery test testing result, the result is as shown in table 1.Can find out that from table main mechanical parameter is compared comparatively approaching with the off-line test data.If as exact value, then three-phase clearance between open contacts precision is respectively: 1.23%, 2.10%, 1.22% with the offline inspection data.Three-phase contacting travel precision is respectively: 3.03%, 5.04%, 3.11%.Therefore can draw, the firm chalaza determination methods of firm branch according to the invention has higher degree of accuracy relatively, and simple, meets the requirement of on-line monitoring system.
The table 1VMB5-12 type vacuum circuit breaker testing result of dispatching from the factory contrasts with algorithm computation result according to the invention