CN108875147A - A kind of VFTO mathematic(al) representation approximating method - Google Patents
A kind of VFTO mathematic(al) representation approximating method Download PDFInfo
- Publication number
- CN108875147A CN108875147A CN201810520125.5A CN201810520125A CN108875147A CN 108875147 A CN108875147 A CN 108875147A CN 201810520125 A CN201810520125 A CN 201810520125A CN 108875147 A CN108875147 A CN 108875147A
- Authority
- CN
- China
- Prior art keywords
- vfto
- waveform
- mathematic
- representation
- fitting
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
The present invention relates to technical field of electric power, more particularly to a kind of VFTO mathematic(al) representation approximating method, on the basis of calculating VFTO transient-wave and analyzing its shared rule, with Fourier's expression formula of Gaussian function, data fitting has been carried out to shared VFTO waveform segment, has obtained VFTO mathematic(al) representation；Include the following steps：VFTO waveform at step 1, selection VFTO single pulse initial time is as equivalent waveform；Step 2, according to the common ground at VFTO waveform at the beginning quarter, VFTO waveform is equivalent to attenuation oscillasion impulse；Step 3 is fitted attenuation oscillasion impulse described in step 2 using Gaussian function, obtains VFTO mathematic(al) representation.This method obtains typical VFTO waveform fitting expression formula, can very easily apply to VFTO as in the numerical value calculating of driving source.
Description
Technical field
The invention belongs to technical field of electric power more particularly to a kind of VFTO mathematic(al) representation approximating methods.
Background technique
Isolator operation in GIS substation is typical electromagnetic interference source, can generate electromagnetic interference, is insulated to GIS
It is damaged with associated electrical instrument.When GIS switch operation, contact gap can generate Hf gas discharge process, SF6Quickly hit
It wears, contact both end voltage moment is fallen (ns), and then forms step voltage traveling wave, and passes along disconnecting switch fracture to two sides
It broadcasts, while folding, reflex can be generated.GIS isolator operation speed is slow (5m/s), therefore during isolator operation
SF6Gas can puncture (primary complete operating process about multiple punch-through) repeatedly.Breakdown can all generate edge every time
The traveling wave that contact two sides are propagated.During GIS whole operation, these are because of SF6The traveling wave of breakdown and generation, passes through inside GIS
It rolls over repeatedly, ultimately form very fast transient overvoltage VFTO phenomenon after reflection and multiple stacking.
The typical overall process waveform of VFTO is as shown in Figure 1.Existing research shows that VFTO has the property that：1. rise
Between short, about 3ns~20ns, voltage build-up rate is high, up to 40MV/s；2. frequency is high, it is concentrated mainly on 0.5MHz~150MHz
In range, 300MHz is reached as high as；3. typical amplitude value is 1.5~2.0p.u., 2.5p.u. is reached as high as.
High frequency VFTO can not only destroy the insulation inside GIS, also have part VFTO wave and be coupled to outside GIS, jeopardize with
The primary equipment (such as transformer, overhead transmission line) of GIS connection causes GIS case body current potential to increase, and causes transient state ground potential
It rises (Transient ground potential rise Enclosure Voltage, TGPR), which is also known as transient enclosure
Current potential increases (Transient enclosure voltage, TEV).
Currently, researcher's mainly research of the expansion to VFTO transient characterisitics in terms of simulation calculation and field measurement two.
To VFTO mistake in terms of being mainly to establish gas with various discharging model, different circuit residual charges in VFTO simulation analysis
The influence of the amplitude, transient frequency, wave steepness of voltage；Under the conditions of mainly measuring different operation in VFTO field measurement, isolation
VFTO waveform caused by switching.At present for extra-high voltage GIS measure circuit, mainly using hand hole formula voltage sensor to VFTO into
Row field measurement, and existing research person VFTO has carried out on-the-spot test.In general, basic at present for the investigative technique of VFTO
Become mature.Since VFTO waveform is excessively complicated, the VFTO pulse of single also contains a large amount of data, causes not having also at present
There is the research being effectively fitted for VFTO waveform.
Summary of the invention
The object of the present invention is to provide the mathematical tables that a kind of GIS switching operation in substation causes fast transient overvoltage VFTO
Up to the method for formula.
To achieve the above object, the technical solution adopted by the present invention is that：A kind of VFTO mathematic(al) representation approximating method, is being counted
On the basis of calculating VFTO transient-wave and analyzing its shared rule, with Fourier's expression formula of Gaussian function, to shared
VFTO waveform segment has carried out data fitting, obtains VFTO mathematic(al) representation；Include the following steps：
VFTO waveform at step 1, selection VFTO single pulse initial time is as equivalent waveform；
Step 2, according to the common ground at VFTO waveform at the beginning quarter, VFTO waveform is equivalent to attenuation oscillasion impulse；
Step 3 is fitted attenuation oscillasion impulse described in step 2 using Gaussian function, obtains VFTO mathematical expression
Formula.
In above-mentioned VFTO mathematic(al) representation approximating method, at selection VFTO single pulse initial time described in step 1
VFTO waveform be the 0-1 period single pulse.
In above-mentioned VFTO mathematic(al) representation approximating method, attenuation oscillasion impulse described in step 2 is on power-frequency voltage component
Overlapped high-frequency Decay Rate wave, the waveform attenuating time of VFTO waveform fitting section are 35ns or so.
In above-mentioned VFTO mathematic(al) representation approximating method, using Gaussian function to described in step 2 described in step 3
Attenuation oscillasion impulse, which is fitted, to be included the following steps：
1) data around extreme value in matched curve are fitted using part Gaussian function first；
2) the connected local functions that fitting obtains are connected again, constitutes global fitting function.
The beneficial effects of the invention are as follows：Approximating method through the invention obtains typical VFTO waveform fitting expression formula, energy
Enough VFTO that very easily applies to are as in the numerical value calculating of driving source.
Detailed description of the invention
Fig. 1 is one embodiment of the invention VFTO mathematic(al) representation approximating method flow chart；
Fig. 2 (a) is VFTO transient state all-wave in one embodiment of the invention, and Fig. 2 (b) is selected by one embodiment of the invention
Partial enlargement transient-wave at different single pulse wave bands in VFTO all-wave；
Fig. 3 is one embodiment of the invention by the partial enlargement waveform of VFTO single pulse waveform different moments；Fig. 3 (a)
For 0.0195 moment waveform, Fig. 3 (b) is 0.0139 moment waveform, and Fig. 3 (c) is 0.0115 moment waveform, and Fig. 3 (d) is
0.0090015 moment waveform；
Fig. 4 is the selection schematic diagram that one embodiment of the invention VFTO is fitted waveform；
Fig. 4 (a) is the VFTO fitting waveform chosen, and Fig. 4 (b) is VFTO transient-wave.
Specific embodiment
Embodiments of the present invention are described in detail with reference to the accompanying drawing.
The present embodiment is difficult to accurately obtain the VFTO mathematic(al) representation with universal law for the current country, proposes one kind
GIS switching operation in substation causes the mathematic(al) representation approximating method of fast transient overvoltage.
The present embodiment is on the basis of accurately calculating VFTO transient-wave and analyzing its shared rule, with Gaussian high
Fourier's expression formula of this function has carried out data fitting to shared VFTO waveform segment, has obtained VFTO accurate mathematical expression formula
Method.A kind of VFTO mathematic(al) representation approximating method, includes the following steps：
1, the discovery of VFTO transient-wave rule
It is VFTO transient state all-wave in the present embodiment as shown in Fig. 2 (a), is the whole of VFTO caused by primary switch operates
Section waveform can be seen that VFTO waveform for decaying Sasser by Fig. 2 (a)；VFTO transient state all-wave is by cluster cluster single pulse waveform
It is formed by stacking.As shown in Fig. 2 (b), locally to be put at single pulse wave bands different in VFTO transient state all-wave selected by the present embodiment
Big transient-wave can be seen that the partial enlargement waveform of VFTO has certain regularity by Fig. 2 (b)；By to a large amount of VFTO
The processing discovery of single pulse waveform partial enlargement, transient-wave variation of the VFTO in " 0-1 " section have more certain regularity.Cause
This, chooses representative VFTO wave band and carries out numerical fitting, has certain representative.
Fig. 3 is partial enlargement waveform of the present embodiment in VFTO single pulse waveform different moments, is primary switch operation
Waveform of the generated VFTO at different moments therefrom chooses the waveform of some cycles, for being fitted to VFTO waveform；
Fig. 3 (a) is 0.0195 moment waveform, and Fig. 3 (b) is 0.0139 moment waveform, and Fig. 3 (c) is 0.0115 moment waveform, and Fig. 3 (d) is
0.0090015 moment waveform.By comparing the VFTO waveform of different moments, transient state wave when " 0-1 " section in single pulse is found
Shape, strong to the impact of GIS device, the VFTO transient characterisitics of selective analysis herein have stronger realistic meaning；And " 0-1 " section
The transient characterisitics for locating VFTO waveform are essentially identical, and waveform amplification section variation tendency is consistent, illustrates have in entire VFTO waveform
Certain representative, research gets up can more search out the expression formula for meeting changing rule.
It is equivalent to the progress of VFTO waveform, the common ground of single VFTO waveform is relied primarily on, in summary Consideration, chosen
VFTO waveform at single pulse " 0-1 " section is analyzed.
2, the equivalent waveform fitting process of VFTO
VFTO waveform carries out equivalent analysis to " 0-1 " section, relies primarily on the common ground of single VFTO waveform.According to
The waveform variation tendency at VFTO " 0-1 " moment, is segmented herein.It can be found that VFTO waveform base in each section of waveform
This has some cycles.It is big in view of data volume contained by selected VFTO waveform segment, it can not be effectively fitted, only choose tool herein
There is " 0-1 " of periodical wave band to carry out data fitting.VFTO waveform at " 0-1 " section is in gradually attenuation trend, waveform fitting
The die-away time of section is in 35ns or so.
3, the selection of VFTO fitting function
It is obtained by the above-mentioned concrete analysis to VFTO waveform, although VFTO waveform changing rule is not fairly obvious,
Variation tendency is closer to classical Decay Rate oscillation wave, therefore selects the characteristic of decaying wave to carry out waveform fitting to VFTO herein.Fig. 4 is
The present embodiment VFTO is fitted the selection schematic diagram of waveform, the waveform of selection and fitting waveform is compared, it can be seen that waveform is quasi-
Close effect；Fig. 3 (a) is the VFTO fitting waveform chosen, and Fig. 3 (b) is VFTO transient-wave.
It selects Gaussian function to carry out mathematic(al) representation fitting to VFTO waveform, is since the fitting segment of selection is presented
Irregular fluctuation reciprocal back and forth, and the extraction to Wave data are general using extraction curve discrete data.And Gaussian letter
Number is Fourier's variation characteristic function, and wherein fourier function feature is an advantage over the analysis to discrete data, and integrates very square
Just.
During Gaussian letter type function model of fit, first with local Gaussian function to extreme value week in matched curve
The data enclosed are fitted, then the connected local functions that fitting obtains are connected, and constitute global fitting function, so that quasi-
The expression formula of conjunction can preferably coincide with matched curve.
As shown in Figure 1, the VFTO mathematic(al) representation approximating method of the present embodiment includes the following steps：
(1) selection of VFTO matched curve
By the amplification analysis to VFTO single pulse waveform, when the VFTO waveform at discovery " 0-1 " moment changes compared with peak value
Carve more regular, and VFTO " 0-1 " moment waveform steepness is big, to the sensitive equipment and insulated part in GIS substation
Impact is strong, more worth research.Therefore the waveform amplification at VFTO single pulse " 0-1 " moment is chosen as equivalent waveform.
The present embodiment relies primarily on the common ground of single VFTO waveform, equivalent to the progress of VFTO waveform, carves at the beginning,
VFTO waveform is made of the waveform of same period substantially, according to variation tendency, can be classified as two sections, wherein each section all has
Certain period.But it is not obvious in view of VFTO wave forms attenuation trend and waveform changing rule, according to existing to operating
The decomposition of voltage waveform is screened, and waveform can be equivalent to attenuation oscillasion impulse, i.e., overlapped high-frequency Decay Rate wave on power-frequency voltage component,
The waveform attenuating time of VFTO waveform fitting section is 35ns or so.
(2) the equivalent waveform fitting process of VFTO
Gaussian function is selected in this VFTO mathematic(al) representation fit procedure, is based primarily upon the following：
First, the VFTO of selection is fitted segment, shows irregular fluctuation reciprocal back and forth, variation tendency rule is not
It is fairly obvious, and the extraction to Wave data, it is general using extraction curve discrete data.
Second, Gaussian function are Fourier's variation characteristic functions, be wherein an advantage over the characteristics of fourier function to from
The analysis of data is dissipated, and is integrated very convenient.
Third, in Gaussian letter type function model of fit, using first with local Gaussian function to fitting song
Data in line around extreme value are fitted, then the connected local functions that fitting obtains are connected, and constitute global fitting
Function enables the expression formula of fitting preferably to coincide with matched curve, and then reinforces the accuracy of fitting.
It should be understood that the part that this specification does not elaborate belongs to the prior art.
Although being described in conjunction with the accompanying a specific embodiment of the invention above, those of ordinary skill in the art should
Understand, these are merely examples, various deformation or modification can be made to these embodiments, without departing from original of the invention
Reason and essence.The scope of the present invention is only limited by the claims that follow.
Claims (4)
1. a kind of VFTO mathematic(al) representation approximating method, characterized in that calculating VFTO transient-wave and analyzing its shared rule
On the basis of, with Fourier's expression formula of Gaussian function, data fitting is carried out to shared VFTO waveform segment, has been obtained
VFTO mathematic(al) representation；Include the following steps：
VFTO waveform at step 1, selection VFTO single pulse initial time is as equivalent waveform；
Step 2, according to the common ground at VFTO waveform at the beginning quarter, VFTO waveform is equivalent to attenuation oscillasion impulse；
Step 3 is fitted attenuation oscillasion impulse described in step 2 using Gaussian function, obtains VFTO mathematic(al) representation.
2. VFTO mathematic(al) representation approximating method as described in claim 1, characterized in that choose VFTO single described in step 1
VFTO waveform at pulse initial time is the single pulse of 0-1 period.
3. VFTO mathematic(al) representation approximating method as described in claim 1, characterized in that attenuation oscillasion impulse described in step 2 is
Overlapped high-frequency Decay Rate wave on power-frequency voltage component, the waveform attenuating time of VFTO waveform fitting section are 35ns or so.
4. VFTO mathematic(al) representation approximating method as described in claim 1, characterized in that use Gaussian described in step 3
Function, which is fitted attenuation oscillasion impulse described in step 2, to be included the following steps：
1) data around extreme value in matched curve are fitted using part Gaussian function first；
2) the connected local functions that fitting obtains are connected again, constitutes global fitting function.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810520125.5A CN108875147B (en) | 2018-05-28 | 2018-05-28 | VFTO mathematical expression fitting method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810520125.5A CN108875147B (en) | 2018-05-28 | 2018-05-28 | VFTO mathematical expression fitting method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108875147A true CN108875147A (en) | 2018-11-23 |
| CN108875147B CN108875147B (en) | 2023-09-22 |
Family
ID=64335133
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810520125.5A Active CN108875147B (en) | 2018-05-28 | 2018-05-28 | VFTO mathematical expression fitting method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108875147B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112730951A (en) * | 2020-12-29 | 2021-04-30 | 国网陕西省电力公司电力科学研究院 | Method for measuring TEV and TGPR in gas insulated substation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011157047A1 (en) * | 2010-06-17 | 2011-12-22 | 西安交通大学 | Capacitive sensor calibration system for measuring fast transient over-voltages |
| WO2016082556A1 (en) * | 2014-11-26 | 2016-06-02 | 国家电网公司 | Gis isolation switch high frequency arc resistance model modeling method |
| CN106896306A (en) * | 2017-04-26 | 2017-06-27 | 国网上海市电力公司 | A kind of GIS oscillatory surges pressure test signal antinoise method |
| CN107153743A (en) * | 2017-05-27 | 2017-09-12 | 山东大学 | VFTO intelligent simulations analysis method and system based on intelligent search algorithm |
-
2018
- 2018-05-28 CN CN201810520125.5A patent/CN108875147B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011157047A1 (en) * | 2010-06-17 | 2011-12-22 | 西安交通大学 | Capacitive sensor calibration system for measuring fast transient over-voltages |
| WO2016082556A1 (en) * | 2014-11-26 | 2016-06-02 | 国家电网公司 | Gis isolation switch high frequency arc resistance model modeling method |
| CN106896306A (en) * | 2017-04-26 | 2017-06-27 | 国网上海市电力公司 | A kind of GIS oscillatory surges pressure test signal antinoise method |
| CN107153743A (en) * | 2017-05-27 | 2017-09-12 | 山东大学 | VFTO intelligent simulations analysis method and system based on intelligent search algorithm |
Non-Patent Citations (2)
| Title |
|---|
| 康 兵等: "几种时频分析方法在 VFTO 频谱分析中的应用", 《高电压技术》, 28 February 2015 (2015-02-28), pages 680 - 686 * |
| 汪剑鸣等: "基于特征参数的脉搏波高斯拟合", 《天津工业大学学报》, 31 October 2016 (2016-10-31), pages 59 - 65 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112730951A (en) * | 2020-12-29 | 2021-04-30 | 国网陕西省电力公司电力科学研究院 | Method for measuring TEV and TGPR in gas insulated substation |
| CN112730951B (en) * | 2020-12-29 | 2023-06-27 | 国网陕西省电力公司电力科学研究院 | Method for measuring TEV and TGPR in gas-insulated substation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108875147B (en) | 2023-09-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103698669B (en) | A kind of high-tension electricity apparatus local discharge pulse separation method | |
| WO2016082556A1 (en) | Gis isolation switch high frequency arc resistance model modeling method | |
| CN104535902A (en) | Partial discharge pulse detection system | |
| Judd et al. | Modelling partial discharge excitation of UHF signals in waveguide structures using Green's functions | |
| CN103424675A (en) | Ultrahigh frequency antenna array partial discharge detection system | |
| CN104796116A (en) | Shock waveform peak value measurement method based on parameter fitting | |
| CN104330599A (en) | Partial-discharge ultrahigh-frequency signal simulation source | |
| Isakova et al. | A differential high-voltage divider | |
| CN204188674U (en) | A kind of local discharge superhigh frequency signal imitation source | |
| CN108875147A (en) | A kind of VFTO mathematic(al) representation approximating method | |
| Ishida et al. | A measurement on electromagnetic noises from ESD generator just before and after ESD testing | |
| Itose et al. | Partial discharge detection by TEV sensors and signal propagation analysis in transformer model | |
| CN203224537U (en) | Linear numeral control pulse signal source for calibration of partial discharging UHF detection system | |
| Reid et al. | High bandwidth measurement of partial discharge pulses in SF6 | |
| Judd et al. | Transfer functions for UHF partial discharge signals in GIS | |
| CN204256113U (en) | Partial discharge pulse's detection system | |
| CN108647469A (en) | A kind of mathematic(al) representation approximating method of TEV mass datas | |
| Ji et al. | Influence of voltage waveforms on partial discharge characteristics of protrusion defect in GIS | |
| Oyama et al. | Development of detection and diagnostic techniques for partial discharges in GIS | |
| Wallace | 6 GHz time domain measurement of fast-transient events | |
| da Silva et al. | Mathematical constraints on the use of transmission line models to investigate the preliminary breakdown stage of lightning flashes | |
| Ghosh et al. | A new method for the estimation of time difference of arrival for localization of partial discharge sources using acoustic detection technique | |
| Patel et al. | Simulation and mathematical analysis of partial discharge measurement in transformer | |
| Wang et al. | An experimental investigation on the correlation between corona current pulse and the ground level ion flow field | |
| CN117031379B (en) | High-frequency characteristic verification circuit of current sensor and verification method thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |