CN101651363B - Loss sampling data processing method based on third-order spline interpolation method - Google Patents
Loss sampling data processing method based on third-order spline interpolation method Download PDFInfo
- Publication number
- CN101651363B CN101651363B CN2009100330087A CN200910033008A CN101651363B CN 101651363 B CN101651363 B CN 101651363B CN 2009100330087 A CN2009100330087 A CN 2009100330087A CN 200910033008 A CN200910033008 A CN 200910033008A CN 101651363 B CN101651363 B CN 101651363B
- Authority
- CN
- China
- Prior art keywords
- loss
- data
- sampling
- sampling data
- sampled data
- 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.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/16—Electric power substations
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Emergency Protection Circuit Devices (AREA)
- Detection And Prevention Of Errors In Transmission (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses a loss sampling data processing method based on the third-order spline interpolation method, which comprises the following steps: the communication situation of a digital protective device is detected in real time, if the frame loss or frame error occurs, the status bit of the corresponding loss sampling data is marked invalid; after the communication state recovers, if the number of lost sampling points is not more than the number of the already set preset lost sampling points which is allowed to be evaluated, the waveform fitting is carried out by using a third-order spline function on the basis of at least three effective sampling data before a loss data segment and at least three effective sampling data after the loss data segment, then the approximate value of the loss sampling data is evaluated in the interpolation method, and the corresponding state of the loss sampling data is marked estimated; if the number of lost sampling points is more than the number of the already set lost sampling points, no measure is taken or only the approximate value of the loss sampling data is evaluated, but the status bit of the corresponding loss sampling data, which is marked invalid, is not changed. The invention solves the problem of the delayed actuation or the locking of the digital protective device, which is resulted from the incomplete data source due to the flame loss or the frame error, thereby ensuring the rapid actuation and the accuracy of the digital protective device.
Description
Technical field
The present invention relates to a kind of sampled data value approximate solution method, be applicable to that digital transformer substation because of the loss sampling data processing that mistake frame/frame losing causes, belongs to technical field of power systems.
Background technology
In the digital transformer substation, after process layer data hop is subjected to serious electromagnetic interference, the generation of error code will be inevitable, mistake frame/frame losing phenomenon occur thereby cause transmitting data.When mistake frame/frame losing took place, the general measure taked of protective device at present was latch-up protection, reopens after the pending data transmission state is stable again.Obvious this way exchanges its reliability for the sacrificing protection responsiveness, if system is during because of the various fault of the claimed quick excision of stable problem, protective value can't meet the demands.Simultaneously, disturbed protective device can't move rapidly, maybe will influence the cooperation of protections at different levels in the system, and then causes the extension of accident impact.
Summary of the invention
Technical problem to be solved by this invention is, overcomes in the digitizing station protective device and guarantees rapidity and accuracy that device moves in the case because of tripping or delay voltage that the loss of small sample data point causes.
For solving the problems of the technologies described above, the invention provides a kind of loss sampling data processing method based on the third-order spline interpolation method, it is characterized in that may further comprise the steps:
1) counts according to the sampled data of the sample frequency of MU in the digital transformer substation and the transmission of every frame and determine the sampling number that allowing loss under the certain calculation accuracy condition satisfying;
2) to protection each sampling number in the buffering area according to setting up mode bit, when the authentic and valid value of corresponding sampled data, put " effectively " and indicate, when invalid, put the engineering noise sign, if estimated value, then put " estimation ", during initialization, all sampled data states all are set to engineering noise;
3) real-time monitoring device signal intelligence, if there is frame losing/mistake frame condition to take place, mark corresponding sample data mode bit is an engineering noise, after treating that communications status recovers, if the sampling number of loss is not more than allowing of having set when asking for the loss sampling number, then respectively to be no less than the efficiently sampling data of 3 points before and after the lost data section, use three rank spline functions to carry out waveform fitting, ask the approximation of loss sampling data then with interpolation method, and the approximation corresponding states of loss sampling data is put " estimation " mark; Otherwise do not take any measure, or only ask the loss sampling data approximation, but do not change the engineering noise sign of its mode bit;
4) when the protection computing, the protection timing element is screened statistics to the sampled data state that participates in computing, when the engineering noise sampled data of finding the participation computing is counted greater than M, and latch-up protection, described M 〉=3; When " estimations " sampled data that find to participate in computing is counted when being not more than N, can be with its protection time-delay corresponding time of counting of increase; When " estimations " sampled data that find to participate in computing is counted greater than N, then latch-up protection, up to the participation computing be " effectively " sampled data till, the scope of described N is 5～10.
The beneficial effect that the present invention reached is as follows: the method for asking for the loss sampling data approximation based on the third-order spline interpolation method provided by the invention; fully adaptive system is rich in the harmonic wave situation; solution is guaranteed the quick-action and the accuracy of protective device because of caused digital protection device action time-delay of the imperfect problem of data source that frame losing/the mistake frame causes or locking.
Description of drawings
Fig. 1 is a match oscillogram after the cubic spline interpolation of 4 point sampling data of 4 successive losses;
Fig. 2 is a match oscillogram after the cubic spline interpolation of 8 point sampling data of 4 successive losses;
Fig. 3 is a match oscillogram after the cubic spline interpolation of 12 point sampling data of 4 successive losses;
Fig. 4 is a match oscillogram after the cubic spline interpolation of 16 point sampling data of 4 successive losses;
Fig. 5 is a match oscillogram after the cubic spline interpolation of 20 point sampling data of 4 successive losses.
Embodiment
A kind of loss sampling data processing method based on the third-order spline interpolation method of the present invention may further comprise the steps:
1) counts according to the sampled data of the sample frequency of MU in the digital transformer substation and the transmission of every frame and determine the sampling number that allowing loss under the certain calculation accuracy condition satisfying.With the MU sample frequency was 4800 point/seconds, and every frame comprises 4 sampled point data instances, under one of percentage error of calculation requires, can ask for the loss sampling data approximation of nearly 12 points, promptly allowed continuous 3 frame traffic data error code or loss to occur;
2) to protection each sampling number in the buffering area according to setting up mode bit, when the authentic and valid value of corresponding sampled data, put " effectively " and indicate, when invalid, put the engineering noise sign, if estimated value, then put " estimation ", during initialization, all sampled data states all are set to engineering noise;
3) real-time monitoring device signal intelligence, if there is frame losing/mistake frame condition to take place, mark corresponding sample data mode bit is an engineering noise, after treating that communications status recovers, if the sampling number of loss is not more than allowing of having set when asking for the loss sampling number, then, use three rank spline functions to carry out waveform fitting respectively to be no less than the efficiently sampling data of 3 points before and after the lost data section, ask the approximation of loss sampling data then with interpolation method, and this data corresponding states is put " estimation " mark; Otherwise do not take any measure, or only ask the loss sampling data approximation, but do not change the engineering noise sign of its mode bit;
When 4) protecting computing, the protection timing element is screened statistics to the sampled data state that participates in computing.When the engineering noise sampled data of finding the participation computing is counted greater than M (desirable 3), latch-up protection; When " estimations " sampled data that find to participate in computing is counted when being not more than N (can be different to its value of dissimilar protections, as to the desirable smaller value of quick-action segment protect, the long delay segment protect is then got higher value), can be with its protection time-delay corresponding time of counting of increase; When " estimations " sampled data that find to participate in computing is counted greater than N, then latch-up protection, up to the participation computing be " effectively " sampled data till.
Fig. 1 to Fig. 5 is the waveform fitting situation ask for the sampled data approximation with the third-order spline interpolation method under different loss sampling datas are counted situation after.
Select for use line end to be provided with 330kV, the 250km high-tension line of serial compensation capacitance and shunt reactor, the current waveform when behind the capacitor three-phase shortcircuit taking place is as simulation object, and its equation is described as:
Be of the influence of check cubic spline interpolation to result of calculation, the 20ms internal loss varying number sampled data situation that begins constantly from t=0 has been carried out emulation, now will lose the one piece of data (comprise sampling number and be respectively 4,8,12,16,20, down together) and the simulation result of two segment datas (two segment datas are 4 effective sampling points at interval) is expressed as follows.
Data sampling rate was 4800 point/seconds in the simulation process, and Fu Shi calculates the DFT algorithm that adopts 96 points/cycle.When carrying out the comparison of DFT result of calculation, the interpolation calculation starting point is fixed on the 11st point (t ≈ 2.2917ms) and locates; When carrying out the waveform fitting comparison, the interpolation calculation starting point is fixed on the 1st point (t=0ms) and locates.
Wherein show 1-table 3 and write down the influence of under different loss sampling datas are counted situation third-order spline interpolation method the protection error of calculation.To its further analysis, can obtain following conclusion: 1) interpolation result is different to protection algorithm (DFT algorithm) the result of calculation influence of different times harmonics, and harmonic number is high more, influences big more; 2) generally speaking, the interpolation result error increases and increases with the loss sampling number, and for first-harmonic, when the loss sampling number was not more than 12, the relative error precision was in some thousandths of, and when equaling 16, the relative error precision is in a few percent.
The DFT result of calculation of table 1 each harmonic wave after using cubic spline interpolation under the loss multi-point sampling data cases
Table 1
Table 2 is DFT errors of calculation of each harmonic wave after using cubic spline interpolation under the loss single hop sampled data situation
Table 2
Table 3 is DFT errors of calculation of each harmonic wave after using cubic spline interpolation under two sections sampled data situations of loss
Table 3
Below announced the present invention with preferred embodiment, so it is not in order to restriction the present invention, and all technical schemes that scheme obtained of taking to be equal to replacement or equivalent transformation all drop in protection scope of the present invention.
Claims (1)
1. loss sampling data processing method based on the third-order spline interpolation method is characterized in that may further comprise the steps:
1) counts according to the sampled data of the sample frequency of MU in the digital transformer substation and the transmission of every frame and determine the sampling number that allowing loss under the certain calculation accuracy condition satisfying;
2) to protection each sampling number in the buffering area according to setting up mode bit, when corresponding sampled data is authentic and valid value, put " effectively " and indicate, when invalid, put the engineering noise sign, if estimated value, then put " estimation ", during initialization, all sampled data states all are set to engineering noise;
3) real-time monitoring device signal intelligence, if there is frame losing/mistake frame condition to take place, mark corresponding sample data mode bit is an engineering noise, after treating that communications status recovers, if the sampling number of loss is not more than allowing of having set when asking for the loss sampling number, then respectively to be no less than the efficiently sampling data of 3 points before and after the lost data section, use three rank spline functions to carry out waveform fitting, ask the approximation of loss sampling data then with interpolation method, and the approximation corresponding states of loss sampling data is put " estimation " mark; Otherwise do not take any measure, or only ask the loss sampling data approximation, but do not change the engineering noise sign of its mode bit;
4) when the protection computing, the protection timing element is screened statistics to the sampled data state that participates in computing, when the engineering noise sampled data of finding the participation computing is counted greater than M, and latch-up protection, described M 〉=3; When " estimations " sampled data that find to participate in computing is counted when being not more than N, with its protection time-delay corresponding time of counting of increase; When " estimations " sampled data that find to participate in computing is counted greater than N, then latch-up protection, up to the participation computing be " effectively " sampled data till, the scope of described N is 5～10.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100330087A CN101651363B (en) | 2009-06-08 | 2009-06-08 | Loss sampling data processing method based on third-order spline interpolation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100330087A CN101651363B (en) | 2009-06-08 | 2009-06-08 | Loss sampling data processing method based on third-order spline interpolation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101651363A CN101651363A (en) | 2010-02-17 |
| CN101651363B true CN101651363B (en) | 2011-07-20 |
Family
ID=41673518
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009100330087A Active CN101651363B (en) | 2009-06-08 | 2009-06-08 | Loss sampling data processing method based on third-order spline interpolation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101651363B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5203440B2 (en) * | 2010-11-24 | 2013-06-05 | 横河電機株式会社 | Harmonic component measuring device |
| CN102243272A (en) * | 2011-04-01 | 2011-11-16 | 重庆大学 | High-precision harmonic analysis method of sampling data synchronization |
| CN106291205B (en) * | 2016-10-14 | 2019-04-16 | 广州视源电子科技股份有限公司 | A kind of blood oxygen probe method for diagnosing faults and device |
| CN108596475B (en) * | 2018-04-23 | 2021-11-30 | 华北电力大学 | PMU data recovery method based on dynamic change of interpolation interval |
| CN109710449A (en) * | 2019-01-10 | 2019-05-03 | 福州芝麻智能科技有限公司 | A kind of scripting language basic data type supports the method and system of invalid value operation |
| CN116933022B (en) * | 2023-09-14 | 2024-01-23 | 国网江苏省电力有限公司苏州供电分公司 | Intelligent synchronous alignment preprocessing method and system for data of multi-source load equipment |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5587917A (en) * | 1994-10-17 | 1996-12-24 | Eaton Corporation | Data collection and processing for digital AC power system monitor/analyzer |
| CN101409463A (en) * | 2008-12-02 | 2009-04-15 | 江西省电力科学研究院 | Protection and video system gang control method for electric power system digitalization transforming plant |
-
2009
- 2009-06-08 CN CN2009100330087A patent/CN101651363B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5587917A (en) * | 1994-10-17 | 1996-12-24 | Eaton Corporation | Data collection and processing for digital AC power system monitor/analyzer |
| CN101409463A (en) * | 2008-12-02 | 2009-04-15 | 江西省电力科学研究院 | Protection and video system gang control method for electric power system digitalization transforming plant |
Non-Patent Citations (1)
| Title |
|---|
| 孙同明等.应用三次样条函数快速计算插值FFT算法.《电力自动化设备》.2007,第27卷(第6期),60-62. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101651363A (en) | 2010-02-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101651363B (en) | Loss sampling data processing method based on third-order spline interpolation method | |
| Biswal et al. | Hilbert Huang transform based online differential relay algorithm for a shunt-compensated transmission line | |
| US9720027B2 (en) | Fault point locating method of hybrid lines based on analysis of comprehensive characteristics of single-end electric quantity and transient travelling waves | |
| Costa et al. | Two-terminal traveling-wave-based transmission-line protection | |
| CN101551433B (en) | Distribution network feed out circuit fault circuit selection method by using HHT detection technology | |
| CN103424669B (en) | A kind of selection method utilizing fault feeder zero-sequence current matrix principal component analysis (PCA) first principal component | |
| CN111398733B (en) | Ocean nuclear power platform power grid ground fault line selection protection method and system | |
| da Silva França et al. | Traveling wave-based transmission line earth fault distance protection | |
| Lin et al. | A novel adaptive single-phase reclosure scheme using dual-window transient energy ratio and mathematical morphology | |
| CN103675602A (en) | Method for discriminating lightning flashover and non-flashover of power transmission lines | |
| CN101635446B (en) | Method for processing sample data based on waveform fitting and secondary sampling | |
| CN103633629A (en) | High-voltage direct current power transmission line protection method based on wavelet transformation and energy spectrum analysis | |
| CN111463764A (en) | Direct-current transmission line protection method based on initial voltage traveling wave frequency domain attenuation rate | |
| CN101707396B (en) | Lost sample data processing method based on mathematic models of current and voltage | |
| CN106159875A (en) | A kind of short circuit current zero predicting method and short circuit current phase selection disjunction control method | |
| CN104393581A (en) | Triangle looped network transient state quantity unit protection method utilizing voltage and current abrupt change quantity wavelet coefficient correlation analysis | |
| CN108808636B (en) | A kind of multi-terminal HVDC transmission line protection method based on transfer entropy | |
| CN102435896A (en) | Intermittent grounding fault rapid identification method of ship medium-voltage power system | |
| CN109193578B (en) | Rapid protection method for direct-current transmission line based on wave-front generalized Logistic function fitting | |
| CN108287286B (en) | Polarity verification method based on single-phase earth fault recording data | |
| CN104393580B (en) | A kind of T of utilization voltage x current Sudden Changing Rate wavelet coefficient correlation analysis gets access to grid transient cell protection method | |
| EP3618212A1 (en) | Method and apparatus for suppressing impact of compensator on line distance protection | |
| CN108847653B (en) | A kind of transformer excitation flow recognition method and device based on waveform change rate | |
| CN104702245A (en) | Method for suppressing interference of surging lightning waves | |
| Wang et al. | Fault location and phase selection for UHV six-phase transmission lines |
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 | ||
| ASS | Succession or assignment of patent right |
Owner name: NARI NANJING CONTROL SYSTEM LTD. Effective date: 20141128 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20141128 Address after: High road high tech Development Zone Nanjing city Jiangsu province 210061 No. 20 Patentee after: NARI Technology Development Co., Ltd. Patentee after: SGCC NARI Nanjing Control System Co., Ltd. Address before: High road high tech Development Zone Nanjing city Jiangsu province 210061 No. 20 Patentee before: NARI Technology Development Co., Ltd. |