CN1924595A - Virtual instrument technique based gas insulation combined electric appliances online detecting method - Google Patents

Virtual instrument technique based gas insulation combined electric appliances online detecting method Download PDF

Info

Publication number
CN1924595A
CN1924595A CN 200610095138 CN200610095138A CN1924595A CN 1924595 A CN1924595 A CN 1924595A CN 200610095138 CN200610095138 CN 200610095138 CN 200610095138 A CN200610095138 A CN 200610095138A CN 1924595 A CN1924595 A CN 1924595A
Authority
CN
China
Prior art keywords
carry out
signal
data
shelf depreciation
monitoring
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 200610095138
Other languages
Chinese (zh)
Other versions
CN100427961C (en
Inventor
张晓星
唐炬
唐世宇
杜林�
莫文强
李剑
刘明军
许中荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CNB200610095138XA priority Critical patent/CN100427961C/en
Publication of CN1924595A publication Critical patent/CN1924595A/en
Application granted granted Critical
Publication of CN100427961C publication Critical patent/CN100427961C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Relating To Insulation (AREA)

Abstract

This invention relates to gas insulation combination apparatus online monitor method based on virtual device in the device insulation monitor technique field. This invention uses the computer and gas insulation combination apparatus local online test and positioning device and through virtual device technique and comprises parameter setting, online monitor, data analysis, history research four units module computer program and all-day GIS local discharging for online monitoring with multiple data analysis, warning and fault finding and requiring function.

Description

Gas insulation combined electric appliances online detecting method based on virtual instrument technique
Technical field
The invention belongs to the insulation of electrical installation monitoring technical field, relate to local discharge of gas-insulator switchgear ultrahigh frequency on-line monitoring method, particularly based on the local discharge of gas-insulator switchgear ultrahigh frequency on-line monitoring method of virtual instrument technique.
Background technology
Gas insulated combined electrical equipment in the electrical equipment (GIS) is in operation, because the existence of the various defectives of GIS inside causes GIS to cause various forms of shelf depreciations, as corona, impurity discharge etc.The inner shelf depreciation of GIS continues for some time and can develop into Fault of Insulating Breakdown, causes power outage, influences the safe operation of electric system, brings about great losses to national economy.The shelf depreciation of the early stage insulation fault of GIS has the characteristic of ultrahigh frequency, and its equivalent frequency is very high, can reach 1GHz.Because GIS is good waveguiding structure, the ultra-high frequency signal of shelf depreciation generation can be propagated therein thus, and leak in the ring flange junction by its insulator, so at the inner local discharge signal that detects local discharge signal or leak of GIS in its external detection, all can in time find the inherent vice of GIS, take appropriate measures, can avoid the generation of electrical equipment malfunction, guarantee the safe operation of electric system.
As Chinese patent application number is 200310118942.1 disclosed " gas insulated transformer substation high-frequency wideband partial discharge monitoring methods ", its local discharge detection device is the enclosure that a bottom surface has the slit, being packaged with frequency band range in the box is the broadband detection antenna of 10~3000MHz, and antenna ends is useful on the SMA high-frequency coaxial connector of feed; The one end is connected with 50 ohm matched impedance.But because the frequency of corona discharge is usually at 100-200MHz in the air, this sensor low-frequency range is low excessively, can not overcome to detect on-the-spot strong corona interference, and influence easily is interfered.Its diagnosis unit application wavelet packet and fractal dimension come the local feature of quantitative test local discharge signal, realization is to the quantitative analysis of local discharge signal waveform character, but owing to always have various noises at the on-line monitoring scene, the existence of noise has great influence to fractal dimension, causes calculating inaccurate.
And for example the application number of the applicant's application is the localization method of 200610054229.9 disclosed " online detection and location device of local discharge of gas-insulator switchgear and localization methods ", though can accurately locate to Partial Discharge Sources, but it can not make accurate judgement to the fault type of shelf depreciation, can not make assessment to the development trend of shelf depreciation, be difficult to more references to the maintainer.
Summary of the invention
The objective of the invention is to weak point at existing GIS Partial Discharge Detection, a kind of combination of gases electrical equipment on-line monitoring method based on virtual instrument technique is provided, have characteristics such as efficient, that use is flexible, powerful, realized the virtual instrument of " software is exactly instrument ", in conjunction with digital processing technologies such as modern signal analyzing and processing and databases, help the GIS staff of transformer station to carry out the judgement of GIS discharge fault type, repair and maintenance decision-making, be convenient to those skilled in the art and easily use.
The application number that the present invention utilizes the applicant to apply for is " online detection and location device of local discharge of gas-insulator switchgear and localization method " disclosed device of 200610054229.9, based on virtual instrument technique, by LabVIEW software platform and high-speed figure collector, adopt the mode bus that is fit to, utilize virtual control to control bottom hardware, the collection, analyzing and processing and the result that are finished signal by powerful software show, have realized functions such as on-line monitoring, data analysis, data in real time preservation, historical query, fault alarm, demonstration.And use special programming software to save a lot of man power and materials' input, shorten product development cycle, develop the product of superior performance.
The object of the present invention is achieved like this: a kind of combination of gases electrical equipment on-line monitoring method based on virtual instrument technique, the application number that utilizes computing machine and the applicant application is the sensor array that the micro-strip paster antenna of locating device in 200610054229.9 " online detection and location device of local discharge of gas-insulator switchgear and localization method " is formed, the multi-way intelligence selector switch, amplifilter, high-speed figure collector and control module, carry out partial discharge monitoring by the GIS on-line monitoring computer program based on virtual instrument technique, this routine package containing parameter is provided with, on-line monitoring, data analysis, four unit modules of historical query.Its concrete method step is as follows:
(1) parameter setting
At first enter the parameter setting: mainly carry out sample frequency, sampling length, regularly sampling parameter setting such as length, the activation threshold value setting, the alarm threshold value setting, monitoring mode (comprise manually and automatically) is provided with, and filtering parameter such as is provided with at the software initial setting up.
(2) systemic-function is selected
After (1) step, the parameter setting was finished, enter the function selecting unit, and automatic selection function 1, enter the on-line monitoring module; When data acquisition finished, selection function 2 entered data analysis module automatically; When operating personnel patrolled and examined, selection function 3 entered the historical query module; In the time of need stopping to monitor, operating personnel's selection function 4, system withdraws from, EOP (end of program).
1) when selecting 1, enters the on-line monitoring module, carry out on-line monitoring.The on-line monitoring unit at first according to the sampling parameter setting in the parameter set unit, carries out initialization to the high-speed figure collector.Carry out monitoring mode then and judge, when monitoring mode was automated manner, the local discharge signal to the sensor acquisition at the tested disc insulator in full station place carried out circulatory monitoring automatically.When monitoring mode is manual mode, select the sensor signal at some tested disc insulators place to concentrate monitoring, manual mode is used to needing when local discharge signal surpasses early warning value to occur the GIS interval of primary part observation.When monitoring mode is automated manner: according to sampling timing the interval is set and starts timer, when interval time then, computing machine is at first by total line traffic control multi-way intelligence switch, micro-strip paster antenna in the gating sensor array in turn, start the high-speed figure collector simultaneously, conversion of signals with each micro-strip paster antenna is a digital signal successively, be sent to by bus again and carry out data processing in the computing machine, and data are shown and stores processor, patrol and examine and return wait when finishing and regularly arrive next time.When monitoring mode was manual mode: computing machine was at first by certain micro-strip paster antenna in total line traffic control multi-way intelligence switching gate sensor array, start the high-speed figure collector simultaneously, conversion of signals with this micro-strip paster antenna is a digital signal successively, be sent to by bus again and carry out data processing in the computing machine, and data are shown and stores processor.Memory module has realized preserving automatically for the real-time raw data of the shelf depreciation of measuring, and database adopts general relevant database, and memory contents comprises local discharge signal waveform and device number, acquisition time and maximum pd quantity etc.After on-line monitoring is finished, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, selection function 2 enters data analysis module automatically.
2) when selecting 2, enter data analysis module, carry out data analysis.Data analysis unit is used for multiple small echo denoising and the spectrum analysis to local discharge data, and fault type is carried out pattern-recognition and Partial Discharge Sources location.
1. small echo denoising again: the application number that utilizes the inventor to invent is the multiple Wavelet noise-eliminating method in 200510057243.x " online detector for partial discharge of gas-insulated substation and denoising method ", adopt multiple wavelet technique that the sensor signal of gathering is carried out denoising, undesired signals such as the on-the-spot white noise that may exist, arrowband are detected in the place to go.According to signal amplitude after the multiple small echo denoising, when greater than alarm threshold value, then carry out sound and light alarm.
2. spectrum analysis: the change procedure that joint time frequency analysis can careful portrayal discharge signal be taken place on time-frequency plane, local discharge signal is carried out spectrum analysis, can reflect the time-frequency characteristic of local discharge signal, be provided with the back and further analyze.
3. pattern-recognition: pattern-recognition is carried out fault diagnosis to the GIS shelf depreciation, so that accurately grasp the defect property and the guide maintenance of GIS inside.Pattern-recognition at first to carry out feature extraction, then the feature of extracting is carried out pattern-recognition by neural network.The concrete steps of its feature extraction are as follows:
1〉PD Signal Pretreatment: the local discharge signal sample is carried out multiple small echo denoising, and with the signal normalization after the denoising;
2〉the multiple wavelet decomposition of processing back signal: pretreated signal is carried out multiple wavelet transformation, and multiple wavelet decomposition number of plies n is the 3-5 layer, obtains the multiple wavelet coefficient of each yardstick;
3〉each scale coefficient fuzzy clustering: to the Fuzzy C-Means Clustering of carrying out of each scale coefficient (real part R, imaginary part I and the composite information R|I| of multiple wavelet coefficient), cluster numbers c is 2-4.
C={c I, 1, c I, 2..., c I, nFor waiting the i yardstick wavelet coefficient of classifying; V={v 1, v 2..., v cBe c cluster centre, v p∈ R PBe subordinate to matrix U=[u P, k] N * c, and satisfy u P, k∈ [0,1], Σ p = 1 c u p , k = 1 , Given preliminary classification matrix U C * n (l), wherein l is an iterations.
The cluster centre vector of preliminary classification then:
V p ( l ) = Σ k = 1 n [ ( u p , k ( l ) ) m C k ] Σ k = 1 n ( u p , k ( l ) ) m - - - ( 1 )
Renewal is subordinate to matrix:
u p , k ( l + 1 ) = 1 Σ j = 1 c ( | | C k - V p ( l ) | | | | C k - V j ( l ) | | ) 2 m - 1 - - - ( 2 )
Objective function J bBe defined as:
J b ( U , V ) = Σ k = 1 n Σ p = 1 c u pk | | c l , k - v p | | 2 - - - ( 3 )
With formula (1-3) repeated calculation cluster centre, classification matrix (being the degree of membership matrix), up to objective function J bReach minimum.
4〉form each scale feature amount: as characteristic quantity, each shelf depreciation sample generates n * c characteristic quantity with the energy of the wavelet coefficient of each cluster under each yardstick.
f i , k = d i , l 2 + · · · + d i , m 2 , k = 1,2 , · · · , n - - - ( 4 )
F wherein I, kBe the characteristic quantity of k cluster under the i yardstick, d I, l..., d I, mBe all wavelet coefficients in k cluster under the i yardstick.
When carrying out pattern-recognition, at first dissimilar shelf depreciation waveform training samples is formed the characteristic quantity storehouse by above-mentioned characteristic extraction step, and neural network training; The waveform sample that the shelf depreciation of on-line monitoring collection in worksite need be classified forms characteristic quantity by above-mentioned characteristic extraction step then, and input neural network carries out the identification of Partial Discharge Sources defect type, recognition result confession technician reference.
4. shelf depreciation source location: the application number that utilizes the inventor to invent carries out the shelf depreciation source location for the Partial Discharge Sources localization method in 200510057243.x " online detection and location device of local discharge of gas-insulator switchgear and localization method ".
In (2)-2) after data analysis finishes, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, selection function 1 enters the on-line monitoring module automatically.
3) when selecting 3, enter the historical query module, carry out historical query.For making things convenient for the operation conditions of power equipment in staff of transformer station and the upper management personnel long-term observation station, monitoring system provides query function.At first inquire about, generate historical data report and historical discharge signal amplitude changing trend diagram at historical shelf depreciation data, and formulation equipment year, the moon and day trend curve inquiry.Because the signal amplitude that uhf sensor is measured is closely related with travel path, so demarcate very big difficulty with discharge capacity PC.The inventor adopts the signal amplitude of measuring UHF local discharge signal waveform by on-line continuous, the historical variations trend of observation signal amplitude, judge the order of severity of shelf depreciation, when signal amplitude is changed significantly, expression built-in electrical insulation degradation trend is remarkable, report to the police to the monitoring personnel with acousto-optic, the technician should in time handle.After historical query is finished, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, proceed historical query.
4) when selecting 4, system withdraws from, EOP (end of program).
The present invention adopts technique scheme, based on graphical programming language LabView and Database Systems, realized a kind of based on virtual instrument technique GIS on-line monitoring computer program, the variation of can be round-the-clock automatic monitoring GIS transformer station multi-point partial discharge signal, carry out the several data analysis, and have early warning, fault diagnosis and a positioning function, multiple query function.Show that through trial run this software has stable, friendly interface is easy to operate, the characteristics that function is complete.The present invention can be widely used in the GIS partial discharge monitoring of generating plant, transformer station.
Description of drawings
Fig. 1 is the theory diagram of the hardware unit of the present invention's utilization;
Fig. 2 is the program flow chart of the inventive method;
Fig. 3 is the program flow diagram of the on-line monitoring of present embodiment;
Fig. 4 is the program flow diagram that the shelf depreciation waveform character of present embodiment extracts;
Fig. 5 is the program flow diagram of the historical query of present embodiment.
Embodiment
Below in conjunction with embodiment, further specify the present invention.
Shown in Fig. 1~5, a kind of gas insulation combined electric appliances online detecting method based on virtual instrument technique, the application number that utilizes computing machine and the applicant application is the device of 200610054229.9 " online detection and location device of local discharge of gas-insulator switchgear and localization method ", carry out the on-line monitoring of GIS shelf depreciation by program, concrete method step is as follows:
(1) parameter setting
At first enter the parameter setting: mainly carry out sample frequency, sampling length, regularly sampling parameter setting such as length, the activation threshold value setting, the alarm threshold value setting, monitoring mode (comprise manually and automatically) is provided with, and filtering parameter such as is provided with at the software initial setting up.
(2) systemic-function is selected
After (1) step, the parameter setting was finished, enter the function selecting unit, and automatic selection function 1, enter the on-line monitoring module; When data acquisition finished, selection function 2 entered data analysis module automatically; When operating personnel patrolled and examined, selection function 3 entered and goes through official's enquiry module; In the time of need stopping to monitor, operating personnel's selection function 4, system withdraws from, EOP (end of program).
1) when selecting 1, enters the on-line monitoring module, carry out on-line monitoring according to program circuit shown in Figure 3.The on-line monitoring unit at first according to the sampling parameter setting in the parameter set unit, carries out initialization to the high-speed figure collector.Carrying out monitoring mode then judges, when monitoring mode is automated manner: according to sampling timing the interval is set and starts timer, when interval time then, computing machine is at first by total line traffic control multi-way intelligence switch, first micro-strip paster antenna in the gating sensor array, start the high-speed figure collector simultaneously, conversion of signals with each micro-strip paster antenna is a digital signal successively, be sent to by bus again and carry out data processing in the computing machine, and data are shown and stores processor, this back that disposes judges once to patrol and examine whether finish, when patrolling and examining when not finishing, next micro-strip paster antenna in the gating sensor array carries out signals collecting, waveform shows and data storage processing; When patrolling and examining end, return and wait for that fixed time interval arrives next time.When monitoring mode was manual mode: computing machine was at first by certain micro-strip paster antenna in total line traffic control multi-way intelligence switching gate sensor array, start the high-speed figure collector simultaneously, conversion of signals with this micro-strip paster antenna is a digital signal successively, be sent to by bus again and carry out data processing in the computing machine, and data are shown and stores processor.Memory module has realized preserving automatically for the real-time raw data of the shelf depreciation of measuring, and database adopts general relevant database, and memory contents comprises local discharge signal waveform and device number, acquisition time and maximum pd quantity etc.After on-line monitoring is finished, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, selection function 2 enters data analysis module automatically.
2) when selecting 2, enter data analysis module, carry out data analysis.Data analysis unit is used for multiple small echo denoising and the spectrum analysis to local discharge data, and fault type is carried out pattern-recognition and Partial Discharge Sources location.
1. small echo denoising again: the application number that utilizes the inventor to invent is the multiple Wavelet noise-eliminating method in 200510057243.x " online detector for partial discharge of gas-insulated substation and denoising method ", adopt multiple wavelet technique that the sensor signal of gathering is carried out denoising, undesired signals such as the on-the-spot white noise that may exist, arrowband are detected in the place to go.According to signal amplitude after the multiple small echo denoising, when greater than alarm threshold value, then carry out sound and light alarm.
2. spectrum analysis: the change procedure that joint time frequency analysis can careful portrayal discharge signal be taken place on time-frequency plane, local discharge signal is carried out spectrum analysis, can reflect the time-frequency characteristic of local discharge signal, be provided with the back and further analyze.
3. pattern-recognition: pattern-recognition at first to carry out feature extraction, then the feature of extracting is carried out pattern-recognition by neural network.Wherein the concrete steps of feature extraction are as shown in Figure 4:
1〉PD Signal Pretreatment: the local discharge signal sample is carried out multiple small echo denoising, and with the signal normalization after the denoising;
2〉the multiple wavelet decomposition of processing back signal: pretreated signal is carried out multiple wavelet transformation, and multiple wavelet decomposition number of plies n is 5 layers, obtains the multiple wavelet coefficient of each yardstick;
3〉each scale coefficient fuzzy clustering: to the Fuzzy C-Means Clustering of carrying out of each scale coefficient (real part R, imaginary part I and the composite information R|I| of multiple wavelet coefficient), cluster numbers c is 3.
C={c I, 1, c I, 2..., c I, nFor waiting the i yardstick wavelet coefficient of classifying; V={v 1, v 2..., v cBe c cluster centre, v p∈ R PBe subordinate to matrix U=[u P, k] N * c, and satisfy u P, k∈ [0,1], Σ p = 1 c u p , k = 1 , Given preliminary classification matrix U C * n (l), wherein l is an iterations.
The cluster centre vector of preliminary classification then:
V p ( l ) = Σ k = 1 n [ ( u p , k ( l ) ) m C k ] Σ k = 1 n ( u p , k ( l ) ) m - - - ( 1 )
Renewal is subordinate to matrix:
u p , k ( l + 1 ) = 1 Σ j = 1 c ( | | C k - V p ( l ) | | | | C k - V j ( l ) | | ) 2 m - 1 - - - ( 2 )
Objective function J bBe defined as:
J b ( U , V ) = Σ k = 1 n Σ p = 1 c u pk | | c l , k - v p | | 2 - - - ( 3 )
With formula (1-3) repeated calculation cluster centre, classification matrix (being the degree of membership matrix), up to objective function J bReach minimum.
4〉form each scale feature amount: as characteristic quantity, each shelf depreciation sample generates 15 characteristic quantities with the energy of the wavelet coefficient of each cluster under each yardstick.
f i , k = d i , l 2 + · · · + d i , m 2 , k = 1,2 , 3 - - - ( 4 )
F wherein I, kBe the characteristic quantity of k cluster under the i yardstick, d I, l, d I, mBe all wavelet coefficients in k cluster under the i yardstick.
When carrying out pattern-recognition, at first dissimilar shelf depreciation waveform training samples is formed the characteristic quantity storehouse by above-mentioned characteristic extraction step, and neural network training; The waveform sample that the shelf depreciation of on-line monitoring collection in worksite need be classified forms characteristic quantity by above-mentioned characteristic extraction step then, and input neural network carries out the identification of Partial Discharge Sources defect type, recognition result confession technician reference.
4. shelf depreciation source location: the application number that utilizes the inventor to invent carries out the shelf depreciation source location for the Partial Discharge Sources localization method in 200510057243.x " online detection and location device of local discharge of gas-insulator switchgear and localization method ".
In (2)-2) after data analysis finishes, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, selection function 1 enters the on-line monitoring module automatically.
3) when selecting 3, enter the historical query module, carry out historical query according to program circuit shown in Figure 5, its concrete steps are as follows:
At first select year, the moon and day date of historical query, and need to select the sensor corresponding equipment number in the sensor array of inquiry, the judgement of whether inquiring about then, when inquiry, the historical shelf depreciation data query that carries out this equipment is handled and is shown; When not inquiring about, return and re-enter the Query Dates selection.Displaying contents comprises historical shelf depreciation data, historical data report and historical discharge signal amplitude changing trend diagram.Operating personnel judge the order of severity of shelf depreciation by observing the historical variations trend of discharge amplitude, and when the discharge amplitude was changed significantly, expression built-in electrical insulation degradation trend was remarkable, reports to the police to the monitoring personnel with acousto-optic, and the technician should in time handle.After historical query is finished, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, proceed historical query.
4) when selecting 4, system withdraws from, EOP (end of program).

Claims (2)

1, a kind of combination of gases electrical equipment on-line monitoring method based on virtual instrument technique, the application number that utilizes computing machine and the applicant application is " online detection and location device of local discharge of gas-insulator switchgear and localization method " middle locating device of 200610054229.9, carry out partial discharge monitoring by the GIS on-line monitoring computer program based on virtual instrument technique, its concrete method step is as follows:
(1) parameter setting
At first enter the parameter setting: mainly carry out sample frequency, sampling length, regularly sampling parameter setting such as length, the activation threshold value setting, the alarm threshold value setting, monitoring mode (comprise manually and automatically) is provided with, and filtering parameter such as is provided with at the software initial setting up;
(2) systemic-function is selected
After (1) step, the parameter setting was finished, enter the function selecting unit, and automatic selection function 1, enter the on-line monitoring module; When data acquisition finished, selection function 2 entered data analysis module automatically; When operating personnel patrolled and examined, selection function 3 entered the historical query module; In the time of need stopping to monitor, operating personnel's selection function 4, system withdraws from, EOP (end of program);
1) when selecting 1, enter the on-line monitoring module, carry out on-line monitoring, the on-line monitoring unit is at first according to the sampling parameter setting in the parameter set unit, the high-speed figure collector is carried out initialization, carrying out monitoring mode then judges, when monitoring mode is automated manner, automatically the local discharge signal to the sensor acquisition at the tested disc insulator in full station place carries out circulatory monitoring, when monitoring mode is manual mode, select the sensor signal at some tested disc insulators place to concentrate monitoring, manual mode is used to needing when local discharge signal surpasses early warning value to occur the GIS interval of primary part observation, when monitoring mode is automated manner: according to sampling timing the interval is set and starts timer, when interval time then, computing machine is at first by total line traffic control multi-way intelligence switch, micro-strip paster antenna in the gating sensor array in turn, start the high-speed figure collector simultaneously, conversion of signals with each micro-strip paster antenna is a digital signal successively, be sent to by bus again and carry out data processing in the computing machine, and data are shown and stores processor, patrol and examine and return wait regularly arrival next time when finishing, when monitoring mode was manual mode: computing machine was at first by certain micro-strip paster antenna in total line traffic control multi-way intelligence switching gate sensor array, start the high-speed figure collector simultaneously, conversion of signals with this micro-strip paster antenna is a digital signal successively, be sent to by bus again and carry out data processing in the computing machine, and data are shown and stores processor, memory module has realized preserving automatically for the real-time raw data of the shelf depreciation of measuring, database adopts general relevant database, memory contents comprises local discharge signal waveform and device number, acquisition time and maximum pd quantity etc., after on-line monitoring is finished, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, selection function 2 enters data analysis module automatically;
2) when selecting 2, enter data analysis module, carry out data analysis, data analysis unit is used for multiple small echo denoising and the spectrum analysis to local discharge data, and fault type is carried out pattern-recognition and Partial Discharge Sources location;
1. small echo denoising again: the application number that utilizes the inventor to invent is the multiple Wavelet noise-eliminating method in 200510057243.x " online detector for partial discharge of gas-insulated substation and denoising method ", adopt multiple wavelet technique that the sensor signal of gathering is carried out denoising, undesired signals such as the on-the-spot white noise that may exist, arrowband are detected in the place to go, according to signal amplitude after the multiple small echo denoising, when greater than alarm threshold value, then carry out sound and light alarm;
2. spectrum analysis: the change procedure that joint time frequency analysis can careful portrayal discharge signal be taken place on time-frequency plane, local discharge signal is carried out spectrum analysis, can reflect the time-frequency characteristic of local discharge signal, be provided with the back and further analyze;
3. pattern-recognition: pattern-recognition is carried out fault diagnosis to the GIS shelf depreciation, so that accurately grasp the defect property and the guide maintenance of GIS inside, pattern-recognition at first to carry out feature extraction, then the feature of extracting is carried out pattern-recognition by neural network, the concrete steps of its feature extraction are as follows:
1〉PD Signal Pretreatment: the local discharge signal sample is carried out multiple small echo denoising, and with the signal normalization after the denoising;
2〉the multiple wavelet decomposition of processing back signal: pretreated signal is carried out multiple wavelet transformation, and multiple wavelet decomposition number of plies n is the 3-5 layer, obtains the multiple wavelet coefficient of each yardstick;
3〉each scale coefficient fuzzy clustering: to the Fuzzy C-Means Clustering of carrying out of each scale coefficient (real part R, imaginary part I and the composite information R|I| of multiple wavelet coefficient), cluster numbers c is 2-4;
C={c I, 1, c I, 2..., c I, nFor waiting the i yardstick wavelet coefficient of classifying; V={v 1, v 2..., v cBe c cluster centre, v p∈ R PBe subordinate to matrix U=[u P, k] N * c, and satisfy u P, k∈ [0,1], Σ p = 1 c u p , k = 1 , Given preliminary classification matrix U C * n (l), wherein l is an iterations;
The cluster centre vector of preliminary classification then:
V p ( l ) = Σ k = 1 n [ ( u p , k ( l ) ) m C k ] Σ k = 1 n ( u p , k ( l ) ) m - - - ( 1 )
Renewal is subordinate to matrix:
u p , k ( l + 1 ) = 1 Σ j = 1 c ( | | C k - V p ( l ) | | | | C k - V j ( l ) | | ) 2 m - 1 - - - ( 2 )
Objective function J bBe defined as:
J b ( U , V ) = Σ k = 1 n Σ p = 1 c u pk | | c i , k - v p | | 2 - - - ( 3 )
With formula (1-3) repeated calculation cluster centre, classification matrix (being the degree of membership matrix), up to objective function J bReach minimum;
4〉form each scale feature amount: as characteristic quantity, each shelf depreciation sample generates n * c characteristic quantity with the energy of the wavelet coefficient of each cluster under each yardstick;
f i , k = d i , l 2 + · · · + d i , m 2 , k = 1,2 , · · · , n - - - ( 4 )
F wherein I, kBe the characteristic quantity of k cluster under the i yardstick, d I, l..., d I, mBe all wavelet coefficients in k cluster under the i yardstick;
When carrying out pattern-recognition, at first dissimilar shelf depreciation waveform training samples is formed the characteristic quantity storehouse by above-mentioned characteristic extraction step, and neural network training; The waveform sample that the shelf depreciation of on-line monitoring collection in worksite need be classified forms characteristic quantity by above-mentioned characteristic extraction step then, and input neural network carries out the identification of Partial Discharge Sources defect type, recognition result confession technician reference;
4. shelf depreciation source location: the application number that utilizes the inventor to invent carries out the shelf depreciation source location for the Partial Discharge Sources localization method in 200510057243.x " online detection and location device of local discharge of gas-insulator switchgear and localization method ";
In (2)-2) after data analysis finishes, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, selection function 1 enters the on-line monitoring module automatically;
3) when selecting 3, enter the historical query module, carry out historical query, for making things convenient for the operation conditions of power equipment in staff of transformer station and the upper management personnel long-term observation station, monitoring system provides query function, at first inquire about at historical shelf depreciation data, generate historical data report and historical discharge signal amplitude changing trend diagram, and formulation equipment year, the moon and day trend curve inquiry, because the signal amplitude that uhf sensor is measured is closely related with travel path, so demarcate very big difficulty with discharge capacity PC, the inventor adopts the signal amplitude of measuring UHF local discharge signal waveform by on-line continuous, the historical variations trend of observation signal amplitude, judge the order of severity of shelf depreciation, when signal amplitude was changed significantly, expression built-in electrical insulation degradation trend was remarkable, report to the police to the monitoring personnel with acousto-optic, the technician should in time handle, and after historical query is finished, withdraws from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, proceed historical query;
4) when selecting 4, system withdraws from, EOP (end of program).
2, according to the described gas insulation combined electric appliances online detecting method of claim 1 based on virtual instrument technique, it is characterized in that, a kind of gas insulation combined electric appliances online detecting method based on virtual instrument technique, the application number that utilizes computing machine and the applicant application is the device of 200610054229.9 " online detection and location device of local discharge of gas-insulator switchgear and localization method ", carry out the on-line monitoring of GIS shelf depreciation by program, concrete method step is as follows:
(1) parameter setting
At first enter the parameter setting: mainly carry out sample frequency, sampling length, regularly sampling parameter setting such as length, the activation threshold value setting, the alarm threshold value setting, monitoring mode (comprise manually and automatically) is provided with, and filtering parameter such as is provided with at the software initial setting up;
(2) systemic-function is selected
After (1) step, the parameter setting was finished, enter the function selecting unit, and automatic selection function 1, enter the on-line monitoring module; When data acquisition finished, selection function 2 entered data analysis module automatically; When operating personnel patrolled and examined, selection function 3 entered the historical query module; In the time of need stopping to monitor, operating personnel's selection function 4, system withdraws from, EOP (end of program);
1) when selecting 1, enter the on-line monitoring module, carry out on-line monitoring according to program circuit shown in Figure 3, the on-line monitoring unit is at first according to the sampling parameter setting in the parameter set unit, the high-speed figure collector is carried out initialization, carrying out monitoring mode then judges, when monitoring mode is automated manner: according to sampling timing the interval is set and starts timer, when interval time then, computing machine is at first by total line traffic control multi-way intelligence switch, first micro-strip paster antenna in the gating sensor array, start the high-speed figure collector simultaneously, conversion of signals with each micro-strip paster antenna is a digital signal successively, is sent to by bus and carries out data processing in the computing machine, and data are shown and stores processor, this back that disposes judges once to patrol and examine whether finish, when patrolling and examining when not finishing, next micro-strip paster antenna in the gating sensor array carries out signals collecting, waveform shows and data storage processing; When patrolling and examining end, return and wait for that fixed time interval arrives next time, when monitoring mode was manual mode: computing machine was at first by certain micro-strip paster antenna in total line traffic control multi-way intelligence switching gate sensor array, start the high-speed figure collector simultaneously, conversion of signals with this micro-strip paster antenna is a digital signal successively, be sent to by bus again and carry out data processing in the computing machine, and data are shown and stores processor, memory module has realized preserving automatically for the real-time raw data of the shelf depreciation of measuring, database adopts general relevant database, memory contents comprises local discharge signal waveform and device number, acquisition time and maximum pd quantity etc. after on-line monitoring is finished, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, selection function 2 enters data analysis module automatically;
2) when selecting 2, enter data analysis module, carry out data analysis, data analysis unit is used for multiple small echo denoising and the spectrum analysis to local discharge data, and fault type is carried out pattern-recognition and Partial Discharge Sources location;
1. small echo denoising again: the application number that utilizes the inventor to invent is the multiple Wavelet noise-eliminating method in 200510057243.x " online detector for partial discharge of gas-insulated substation and denoising method ", adopt multiple wavelet technique that the sensor signal of gathering is carried out denoising, undesired signals such as the on-the-spot white noise that may exist, arrowband are detected in the place to go, according to signal amplitude after the multiple small echo denoising, when greater than alarm threshold value, then carry out sound and light alarm;
2. spectrum analysis: the change procedure that joint time frequency analysis can careful portrayal discharge signal be taken place on time-frequency plane, local discharge signal is carried out spectrum analysis, can reflect the time-frequency characteristic of local discharge signal, be provided with the back and further analyze;
3. pattern-recognition: pattern-recognition at first to carry out feature extraction, then the feature of extracting is carried out pattern-recognition by neural network, wherein the concrete steps of feature extraction as shown in Figure 4:
1〉PD Signal Pretreatment: the local discharge signal sample is carried out multiple small echo denoising, and with the signal normalization after the denoising;
2〉the multiple wavelet decomposition of processing back signal: pretreated signal is carried out multiple wavelet transformation, and the multiple wavelet decomposition number of plies is that n is 5 layers, obtains the multiple wavelet coefficient of each yardstick;
3〉each scale coefficient fuzzy clustering: to the Fuzzy C-Means Clustering of carrying out of each scale coefficient (real part R, imaginary part I and the composite information R|I| of multiple wavelet coefficient), cluster numbers c is 3;
C={c I, 1, c I, 2..., c I, nFor waiting the i yardstick wavelet coefficient of classifying; V={v 1, v 2..., v cBe c cluster centre, v p∈ R PBe subordinate to matrix U=[u P, k] N * c, and satisfy u P, k∈ [0,1], Σ p = 1 c u p , k = 1 , Given preliminary classification matrix U C * n (l), wherein l is an iterations;
The cluster centre vector of preliminary classification then:
V p ( l ) = Σ k = 1 n [ ( u p , k ( l ) ) m C k ] Σ k = 1 n ( u p , k ( l ) ) m - - - ( 1 )
Renewal is subordinate to matrix:
u p , k ( l + 1 ) = 1 Σ j = 1 c ( | | C k - V p ( l ) | | | | C k - V j ( l ) | | ) 2 m - 1 - - - ( 2 )
Objective function J bBe defined as:
J b ( U , V ) = Σ k = 1 n Σ p = 1 c u pk | | c i , k - v p | | 2 - - - ( 3 )
With formula (1-3) repeated calculation cluster centre, classification matrix (being the degree of membership matrix), up to objective function J bReach minimum,
4〉form each scale feature amount: as characteristic quantity, each shelf depreciation sample generates 15 characteristic quantities with the energy of the wavelet coefficient of each cluster under each yardstick;
f i , k = d i , l 2 + · · · + d i , m 2 , k = 1,2 , · · · , n - - - ( 4 )
F wherein I, kBe the characteristic quantity of k cluster under the i yardstick, d I, l..., d I, mBe all wavelet coefficients in k cluster under the i yardstick;
When carrying out pattern-recognition, at first dissimilar shelf depreciation waveform training samples is formed the characteristic quantity storehouse by above-mentioned characteristic extraction step, and neural network training; The waveform sample that the shelf depreciation of on-line monitoring collection in worksite need be classified forms characteristic quantity by above-mentioned characteristic extraction step then, and input neural network carries out the identification of Partial Discharge Sources defect type, recognition result confession technician reference;
4. shelf depreciation source location: the application number that utilizes the inventor to invent carries out the shelf depreciation source location for the Partial Discharge Sources localization method in 200510057243.x " online detection and location device of local discharge of gas-insulator switchgear and localization method ";
In (2)-2) after data analysis finishes, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, selection function 1 enters the on-line monitoring module automatically;
3) when selecting 3, enter the historical query module, carry out historical query according to program circuit shown in Figure 5, its concrete steps are as follows:
At first select year, the moon and day date of historical query, and need to select the sensor corresponding equipment number in the sensor array of inquiry, the judgement of whether inquiring about then, when inquiry, the historical shelf depreciation data query that carries out this equipment is handled and is shown; When not inquiring about, return and re-enter the Query Dates selection, displaying contents comprises historical shelf depreciation data, historical data report and historical discharge signal amplitude changing trend diagram, operating personnel are by observing the historical variations trend of discharge amplitude, judge the order of severity of shelf depreciation, when the discharge amplitude is changed significantly, expression built-in electrical insulation degradation trend is remarkable, reports to the police to the monitoring personnel with acousto-optic, and the technician should in time handle, after historical query is finished, withdraw from judgement, when withdrawing from, returned for (2) step and carry out function selecting; When not withdrawing from, proceed historical query;
4) when selecting 4, system withdraws from, EOP (end of program).
CNB200610095138XA 2006-09-20 2006-09-20 Virtual instrument technique based gas insulation combined electric appliances online detecting method Expired - Fee Related CN100427961C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200610095138XA CN100427961C (en) 2006-09-20 2006-09-20 Virtual instrument technique based gas insulation combined electric appliances online detecting method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200610095138XA CN100427961C (en) 2006-09-20 2006-09-20 Virtual instrument technique based gas insulation combined electric appliances online detecting method

Publications (2)

Publication Number Publication Date
CN1924595A true CN1924595A (en) 2007-03-07
CN100427961C CN100427961C (en) 2008-10-22

Family

ID=37817318

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200610095138XA Expired - Fee Related CN100427961C (en) 2006-09-20 2006-09-20 Virtual instrument technique based gas insulation combined electric appliances online detecting method

Country Status (1)

Country Link
CN (1) CN100427961C (en)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101307756B (en) * 2008-07-08 2010-06-16 河北省电力研究院 Thermoelectric generating set auxiliary machinery fault fast-reducing load feed water pump trip behavior feasibility determination method
CN101852836A (en) * 2010-05-05 2010-10-06 华北电力大学 System of detecting partial discharge in equipment
CN102096025A (en) * 2009-12-10 2011-06-15 上海莫克电子技术有限公司 Local discharge detector of gas insulating switch
CN102096026A (en) * 2009-12-11 2011-06-15 上海莫克电子技术有限公司 Portable partial discharge detector for insulating switch
CN102135579A (en) * 2010-02-03 2011-07-27 乌鲁木齐希望电子有限公司 Online monitoring and managing system for power consumption information
CN102279354A (en) * 2011-07-01 2011-12-14 西安交通大学 Cruciform ultrasound array sensor and method for positioning transformer partial discharge
CN102288882A (en) * 2011-08-16 2011-12-21 重庆大学 Partial discharge monitoring system and method for multichannel electrical equipment
CN102353877A (en) * 2011-06-10 2012-02-15 河南电力试验研究院 Partial discharge detection method for realizing phase resolved partial discharge (PRPD) pattern and time resolved partial discharge (TRPD) pattern
CN102539527A (en) * 2010-12-27 2012-07-04 湖北省电力试验研究院 GIS (gas insulated substation) partial discharge mode identification method based on ultrasonic testing
CN102621434A (en) * 2012-04-19 2012-08-01 山东大学 Nonlinear fuzzy detection method for operating safety of power cable tunnel
CN102645621A (en) * 2012-05-17 2012-08-22 广东电网公司电力科学研究院 Multisource partial discharge detection method and device of transformer substation based on space characteristic parameters
CN103278789A (en) * 2013-04-26 2013-09-04 国家电网公司 Simulative discharger for GIS (gas-insulated switchgear) partial discharge detection
CN105652168A (en) * 2016-03-22 2016-06-08 黄旭 Transformer substation local discharge signal detection system with communication device
CN106646165A (en) * 2016-12-23 2017-05-10 西安交通大学 Method and system for classifying and positioning internal insulation defects of GIS
CN108333443A (en) * 2017-12-22 2018-07-27 国家电网公司 The alarm method of power equipment intermittent defect
CN109061426A (en) * 2018-11-02 2018-12-21 国网河北省电力有限公司电力科学研究院 Partial discharge of transformer method for diagnosing faults and on-Line Monitor Device
CN113933660A (en) * 2021-08-31 2022-01-14 华能澜沧江水电股份有限公司 Generator insulation fault degradation monitoring method based on partial discharge characteristics
CN116340312A (en) * 2022-11-01 2023-06-27 湖北能源集团新能源发展有限公司 Data storage method and system of GIS equipment partial discharge on-line monitoring system

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4648468B2 (en) * 2009-06-05 2011-03-09 ファナック株式会社 Wire-cut electrical discharge machine with machining status judgment function
CN102735999B (en) * 2012-07-02 2014-12-10 山东电力集团公司电力科学研究院 Gas insulated substation (GIS) partial discharge online monitoring system and fault mode identifying method thereof

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10210647A (en) * 1997-01-24 1998-08-07 Tohoku Electric Power Co Inc Insulation abnormality diagnosing equipment
JP4157636B2 (en) * 1999-01-28 2008-10-01 株式会社日立製作所 Partial discharge diagnostic equipment for gas insulation equipment
JP3901045B2 (en) * 2002-07-18 2007-04-04 日新電機株式会社 Sensor device and insulation diagnostic device
CN1230684C (en) * 2003-12-18 2005-12-07 西安交通大学 High-frequency wide-band local discharging on-line monitoring method in gas insulative converting station
CN100363747C (en) * 2005-04-26 2008-01-23 重庆大学 Gas insulation combined electric device local discharge superhigh frequency detection apparatus and method
CN100458454C (en) * 2005-08-26 2009-02-04 重庆大学 Online detector for partial discharge of gas-insulated substation and noise reduction method

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101307756B (en) * 2008-07-08 2010-06-16 河北省电力研究院 Thermoelectric generating set auxiliary machinery fault fast-reducing load feed water pump trip behavior feasibility determination method
CN102096025A (en) * 2009-12-10 2011-06-15 上海莫克电子技术有限公司 Local discharge detector of gas insulating switch
CN102096026A (en) * 2009-12-11 2011-06-15 上海莫克电子技术有限公司 Portable partial discharge detector for insulating switch
CN102135579A (en) * 2010-02-03 2011-07-27 乌鲁木齐希望电子有限公司 Online monitoring and managing system for power consumption information
CN101852836A (en) * 2010-05-05 2010-10-06 华北电力大学 System of detecting partial discharge in equipment
CN101852836B (en) * 2010-05-05 2012-07-25 华北电力大学 System of detecting partial discharge in equipment
CN102539527B (en) * 2010-12-27 2013-06-19 湖北省电力公司电力科学研究院 GIS (gas insulated substation) partial discharge mode identification method based on ultrasonic testing
CN102539527A (en) * 2010-12-27 2012-07-04 湖北省电力试验研究院 GIS (gas insulated substation) partial discharge mode identification method based on ultrasonic testing
CN102353877A (en) * 2011-06-10 2012-02-15 河南电力试验研究院 Partial discharge detection method for realizing phase resolved partial discharge (PRPD) pattern and time resolved partial discharge (TRPD) pattern
CN102279354A (en) * 2011-07-01 2011-12-14 西安交通大学 Cruciform ultrasound array sensor and method for positioning transformer partial discharge
CN102288882A (en) * 2011-08-16 2011-12-21 重庆大学 Partial discharge monitoring system and method for multichannel electrical equipment
CN102621434A (en) * 2012-04-19 2012-08-01 山东大学 Nonlinear fuzzy detection method for operating safety of power cable tunnel
CN102645621A (en) * 2012-05-17 2012-08-22 广东电网公司电力科学研究院 Multisource partial discharge detection method and device of transformer substation based on space characteristic parameters
CN102645621B (en) * 2012-05-17 2014-11-19 广东电网公司电力科学研究院 Multisource partial discharge detection method and device of transformer substation based on space characteristic parameters
CN103278789A (en) * 2013-04-26 2013-09-04 国家电网公司 Simulative discharger for GIS (gas-insulated switchgear) partial discharge detection
CN103278789B (en) * 2013-04-26 2015-09-23 国家电网公司 GIS office puts and detects simulation discharger
CN105652168A (en) * 2016-03-22 2016-06-08 黄旭 Transformer substation local discharge signal detection system with communication device
CN106646165A (en) * 2016-12-23 2017-05-10 西安交通大学 Method and system for classifying and positioning internal insulation defects of GIS
CN106646165B (en) * 2016-12-23 2020-06-09 西安交通大学 GIS internal insulation defect classification and positioning method and system
CN108333443A (en) * 2017-12-22 2018-07-27 国家电网公司 The alarm method of power equipment intermittent defect
CN109061426A (en) * 2018-11-02 2018-12-21 国网河北省电力有限公司电力科学研究院 Partial discharge of transformer method for diagnosing faults and on-Line Monitor Device
CN113933660A (en) * 2021-08-31 2022-01-14 华能澜沧江水电股份有限公司 Generator insulation fault degradation monitoring method based on partial discharge characteristics
CN116340312A (en) * 2022-11-01 2023-06-27 湖北能源集团新能源发展有限公司 Data storage method and system of GIS equipment partial discharge on-line monitoring system
CN116340312B (en) * 2022-11-01 2024-06-04 湖北能源集团麻城风电有限公司 Data storage method and system of GIS equipment partial discharge on-line monitoring system

Also Published As

Publication number Publication date
CN100427961C (en) 2008-10-22

Similar Documents

Publication Publication Date Title
CN100427961C (en) Virtual instrument technique based gas insulation combined electric appliances online detecting method
CN103645425B (en) High-voltage cable insulation defect partial discharge on-line monitoring diagnosis method
CN105606975B (en) A kind of orientable superfrequency cable local discharge detection method and device
CN1234014C (en) Power cable local discharge on-line monitoring method and device
CN107831413B (en) The electromagnetic wave detection device and method of transmission line of electricity corona discharge based on unmanned plane
CN100483144C (en) Corona detecting method and system for insulator remote surface on optical amplifier
CN108344931B (en) Power equipment safety analysis system based on ultraviolet spectrum technology
CN109649432B (en) System and method for monitoring integrity of steel rail of cloud platform based on guided wave technology
CN111428561B (en) Broadband electromagnetic interference source identification method and system, storage medium and equipment
CN101702002A (en) Local discharge electroacoustic united detection positioning system and positioning method thereof
CN109490730A (en) Cable discharge detection method, device, storage medium and processor
CN103134765A (en) Chinese medicine sample authenticity preliminary screening method based on terahertz time-domain spectrum
CN1743839A (en) Structure defect ultrasonic on-line intelligent identifying system and identifying method
CN106404399B (en) Method for Bearing Fault Diagnosis based on self-adaptive redundant Lifting Wavelet packet decomposition tree
WO2023284127A1 (en) Vibration-signal-based online monitoring system for gil defects, and method
CN112986870A (en) Distributed power transformer winding state monitoring method and system based on vibration method
CN102616156A (en) Device and method for detecting and elevating contact network state based on electrical capacity analysis
CN111999614A (en) Local discharge multi-source fusion detection early warning method for high-voltage switch cabinet
CN2809640Y (en) Ultrasonic online intelligent recognition system for structural defect
CN117761469A (en) Inspection method and system for diagnosing partial discharge of overhead line of power distribution network
CN206132913U (en) Special high frequency cable partial discharge detector device that can fix a position
CN113156278A (en) CNN network-based GIS partial discharge positioning method
CN109917245B (en) Ultrasonic detection partial discharge signal mode identification method considering phase difference
Zhang et al. A convolutional neural network-based UHF partial discharge atlas classification system for GIS
Bhukya et al. Wavelet Aided Deep Neural Network for Classification of Partial Discharge Sources in Air Insulated Substations

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20081022

Termination date: 20130920