CN102156788A - Method for simulating transmission property of partial discharge signal in power cable - Google Patents

Method for simulating transmission property of partial discharge signal in power cable Download PDF

Info

Publication number
CN102156788A
CN102156788A CN2011101035748A CN201110103574A CN102156788A CN 102156788 A CN102156788 A CN 102156788A CN 2011101035748 A CN2011101035748 A CN 2011101035748A CN 201110103574 A CN201110103574 A CN 201110103574A CN 102156788 A CN102156788 A CN 102156788A
Authority
CN
China
Prior art keywords
power cable
finite
cable
discharge signal
moment
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
CN2011101035748A
Other languages
Chinese (zh)
Other versions
CN102156788B (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN 201110103574 priority Critical patent/CN102156788B/en
Publication of CN102156788A publication Critical patent/CN102156788A/en
Application granted granted Critical
Publication of CN102156788B publication Critical patent/CN102156788B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Relating To Insulation (AREA)

Abstract

The invention discloses a method for simulating transmission property of a partial discharge signal in a power cable, comprising the following steps of: step 1, creating an equivalent circuit model of a coaxial cable under a high frequency signal through an R-L-C network according to electrical parameters of the power cable, and calculating R, L and C distribution parameter matrixes of the equivalent circuit model through a finite element method; step 2, creating a frequency varying factor related finite-different time-domain method iterative model according to the R, L and C distribution parameter matrixes; and step 3, simulating the transmission property of the signal in the cable through the finite-different time-domain method iterative model. The method offers effective help for in-depth study on the transmission property of the partial discharge in the cable, particularly on high-frequency transmission mechanism, and implementing accurate detection and positioning of partial discharge.

Description

A kind of method of simulating local discharge signal propagation characteristic in power cable
 
Technical field
The present invention relates to the shelf depreciation field, relate in particular to a kind of method of simulating local discharge signal propagation characteristic in power cable.
Background technology
The inner shelf depreciation (PD) that takes place of cable can produce the pulse of frequency up to hundreds of MHz, can be subjected to serious decay and wave form distortion effect when in cable, propagating by the pulse of shelf depreciation generation, therefore for realizing the accurate detection and the location of shelf depreciation, must the high-frequency propagation mechanism of the further investigation shelf depreciation propagation characteristic in cable, especially shelf depreciation in cable.At present, all there is the high and high problem of detection cost of price in the equipment of detection local discharge signal propagation characteristic in cable.
Summary of the invention
At the deficiency that prior art exists, the present invention proposes a kind of method of simulation local discharge signal propagation characteristic in power cable with low cost, simple to operate.
For solving the problems of the technologies described above, the present invention adopts following technical scheme:
A kind of method of simulating local discharge signal propagation characteristic in power cable may further comprise the steps successively:
Step 1 according to the electrical parameter of power cable, adopts R- L- CNetwork is set up the equivalent-circuit model of concentric cable under the high-frequency signal, and the calculating equivalent-circuit model R, L, CThe distribution parameter matrix;
Step 2, according to R, L, CThe distribution parameter matrix is set up the Finite-Difference Time-Domain Method iterative model that relates to the factor that frequently becomes of power cable:
Figure 2011101035748100002DEST_PATH_IMAGE003
Wherein, ,
Figure 2011101035748100002DEST_PATH_IMAGE005
Expression respectively
Figure 726158DEST_PATH_IMAGE006
,
Figure 2011101035748100002DEST_PATH_IMAGE007
Moment power cable
Figure 436101DEST_PATH_IMAGE008
The voltage of position;
Figure 2011101035748100002DEST_PATH_IMAGE009
Expression
Figure 877577DEST_PATH_IMAGE006
Moment power cable
Figure 241563DEST_PATH_IMAGE010
The voltage of position;
Figure 2011101035748100002DEST_PATH_IMAGE011
,
Figure 415186DEST_PATH_IMAGE012
Expression respectively
Figure 2011101035748100002DEST_PATH_IMAGE013
Moment power cable
Figure 864622DEST_PATH_IMAGE008
,
Figure 47473DEST_PATH_IMAGE014
The electric current of position;
Figure 2011101035748100002DEST_PATH_IMAGE015
Expression
Figure 385DEST_PATH_IMAGE016
Moment power cable
Figure 613419DEST_PATH_IMAGE008
The electric current of position; LBe distributed inductance;
Figure 2011101035748100002DEST_PATH_IMAGE017
Be distributed capacitance; RBe resistance;
Figure 284572DEST_PATH_IMAGE018
Be the space interval of choosing;
Figure 2011101035748100002DEST_PATH_IMAGE019
Be the time interval of choosing;
The Finite-Difference Time-Domain Method iterative model of power cable head end is:
Figure 2011101035748100002DEST_PATH_IMAGE021
Wherein,
Figure 333430DEST_PATH_IMAGE022
,
Figure 2011101035748100002DEST_PATH_IMAGE023
Expression respectively
Figure 157161DEST_PATH_IMAGE006
,
Figure 187434DEST_PATH_IMAGE007
The voltage of moment power cable head end;
Figure 96615DEST_PATH_IMAGE024
Expression
Figure 808219DEST_PATH_IMAGE013
The electric current of moment power cable head end; ,
Figure 545844DEST_PATH_IMAGE026
For
Figure 684701DEST_PATH_IMAGE006
,
Figure 330446DEST_PATH_IMAGE007
The voltage of moment signal source;
Figure 2011101035748100002DEST_PATH_IMAGE027
Internal resistance for signal source;
Figure 393211DEST_PATH_IMAGE017
Be distributed capacitance;
Figure 440801DEST_PATH_IMAGE028
Be extra electric field intensity;
Figure 563609DEST_PATH_IMAGE018
Be the space interval of choosing;
Figure 634333DEST_PATH_IMAGE019
Be the time interval of choosing;
The Finite-Difference Time-Domain Method iterative model of power cable end is:
Wherein,
Figure 2011101035748100002DEST_PATH_IMAGE031
, Expression respectively
Figure 945863DEST_PATH_IMAGE006
,
Figure 238304DEST_PATH_IMAGE007
The voltage of moment power cable end;
Figure 2011101035748100002DEST_PATH_IMAGE033
Expression
Figure 961979DEST_PATH_IMAGE013
The electric current of moment power cable end; Internal impedance for load;
Figure 120876DEST_PATH_IMAGE017
Be distributed capacitance;
Figure 228509DEST_PATH_IMAGE028
Be extra electric field intensity;
Figure 436768DEST_PATH_IMAGE018
Be the space interval of choosing;
Figure 985561DEST_PATH_IMAGE019
Be the time interval of choosing;
Step 3 adopts the described Finite-Difference Time-Domain Method iterative model of step 2 to simulate the propagation characteristic of local discharge signal in power cable.
The above-mentioned equivalent-circuit model of concentric cable of setting up under the high-frequency signal is to adopt many conductor transmission line theory.
Above-mentioned R, L, CParameter matrix is to adopt finite element method to calculate.
Above-mentioned Finite-Difference Time-Domain Method electric current iteration mould
Figure 135919DEST_PATH_IMAGE003
Revise by introducing the time domain convolution.
Compared with prior art, the present invention has the following advantages and beneficial effect:
1, the inventive method has advantage with low cost, simple to operate;
2, the inventive method is applied to finding the solution of cablebreak movable model with Finite-Difference Time-Domain Method (FDTD) iterative model, has solved decoupling zero problem of difficult in the past;
3, the inventive method is further investigation the shelf depreciation propagation characteristic in power cable, particularly high-frequency propagation mechanism, and realizes that the accurate detection and the location of shelf depreciation provide effective help;
4, the inventive method adopts finite element method to calculate the distribution parameter of equivalent-circuit model, it is truer to the simulation of kelvin effect under the high frequency that finite element method calculates distribution parameter, the frequency that more can reflect distribution parameter becomes factor, has changed in the past to calculate the accurately deficiency of match frequency influence of distribution parameter by experimental formula.
Description of drawings
Fig. 1 is the process flow diagram of the embodiment of the invention;
Fig. 2 is the spatial spreading synoptic diagram of equivalent-circuit model of the present invention;
Fig. 3 is the Finite-Difference Time-Domain Method iterative process to the equivalent-circuit model of Fig. 2;
Fig. 4 is the demonstration test wiring diagram;
The analog result figure of Fig. 5 for adopting the inventive method to obtain;
The measurement result figure of Fig. 6 for adopting demonstration test shown in Figure 4 to obtain.
Embodiment
For local discharge signal, its frequency can reach more than hundreds of MHz, should not set up the lumped-parameter circuit model this moment under low frequency, and should set up the equivalent-circuit model of concentric cable under the high-frequency signal, according to the equivalent-circuit model of concentric cable under the high-frequency signal, set up proper model and come the propagation characteristic of simulating signal in cable then.
In conjunction with the accompanying drawings technical scheme of the present invention is described further below by embodiment.
Be illustrated in figure 1 as a kind of method of simulating local discharge signal propagation characteristic in power cable that the present invention proposes, specific as follows:
Step 1 according to the electrical parameter of power cable, adopts R- L- CNetwork is set up the equivalent-circuit model of concentric cable under the high-frequency signal, and the calculating equivalent-circuit model R, L, CThe distribution parameter matrix:
Be the propagation characteristic of research local discharge signal in power cable, can set up the equivalent-circuit model of concentric cable under the high-frequency signal.Concentric cable can be regarded as by many conductor transmission line and constitute under the high-frequency signal, according to many conductor transmission line theory, with heart yearn, inner semiconductor layer, outer semiconductor layer and the screen layer of cable inside all equivalence be transmission line, again since electricity lead GVery I is ignored, so adopt R- L- CNetwork is set up the equivalent-circuit model of concentric cable under the high-frequency signal.
The distribution parameter of equivalent transmission line can obtain by the cablebreak movable model under the following time domain:
Figure 402953DEST_PATH_IMAGE036
（1）
（2）
Wherein, CFor the distributed capacitance matrix, LFor the distributed inductance matrix, RFor resistor matrix, GBe conductance matrix, V F ( z, t) and I F ( z, t) be respectively voltage source and current source that the external electromagnetic field excitation produces, zRepresentation space, tExpress time.
Because electricity is led in the present embodiment GVery I is ignored, so formula (2) can be written as:
Figure 2011101035748100002DEST_PATH_IMAGE039
（3）
At present, adopt the distribution parameter matrix of the theoretical equivalent-circuit model of setting up of many conductor transmission line mainly to find the solution by experimental formula, but experimental formula can not accurately be simulated the influence of kelvin effect under the high frequency, greatly reduce the accuracy that distribution parameter calculates, make the difference of equivalent-circuit model and cable actual transmissions characteristic increase.The employing finite element method calculates equivalent-circuit model in the present embodiment R, L, CThe distribution parameter matrix, wherein, distributed capacitance is owing to be subjected to frequency influence less, so the distributed capacitance matrix CUnder electrostatic field, calculate; And distributed inductance and resistance are obvious owing to changed by frequency influence, so the distributed inductance matrix LAnd resistor matrix RThe time calculate under the humorous magnetic field.Finite element method more can the real simulated high frequency under the influence of kelvin effect, improved the accuracy that distribution parameter calculates.
The calculating of distributed capacitance is based on the energy principle shown in the formula (4) in the present embodiment:
Figure 2011101035748100002DEST_PATH_IMAGE041
（3）
In the formula (4),
W For calculating the energy in the field domain;
E Be electric field intensity;
D Be electric flux density;
u i Be iExcitation potential on the conductor;
q i For iThe quantity of electric charge on the conductor;
Wherein, q i Can be expressed as:
Figure 2011101035748100002DEST_PATH_IMAGE043
（5）
To obtain in formula (5) the substitution formula (3):
Figure 2011101035748100002DEST_PATH_IMAGE045
（6）
In the formula, KMatrix is the induction coefficient matrix, is the distributed capacitance matrix that need find the solution C
The calculating of distributed inductance and resistance then is based on complex impedance and applied power principle, adopts following formula:
Figure 2011101035748100002DEST_PATH_IMAGE047
（7）
Figure 2011101035748100002DEST_PATH_IMAGE049
（8）
Figure 2011101035748100002DEST_PATH_IMAGE051
（9）
Figure 2011101035748100002DEST_PATH_IMAGE053
（10）
In the formula,
ωBe angular frequency;
V IR , V II Be respectively coil iOn the real part and the imaginary part of current potential effective value;
V JR , V JI Be respectively coil jOn the real part and the imaginary part of current potential effective value;
I IR , I II Be respectively coil iIn the real part and the imaginary part of exciting current effective value;
L Ii , L Ij Be respectively conductor iSelf-induction and conductor i, jBetween mutual inductance;
R Ii , R Ij Be respectively conductor iResistance certainly and conductor i, jBetween mutual resistance.
In this enforcement heart yearn, inner semiconductor layer, outer semiconductor layer and the screen layer of cable inside are set up the transmission system model as 4 transmission lines, therefore calculate R, L, CParameter matrix is 4 rank matrixes.
Step 2, according to R, L, CThe distribution parameter matrix is set up the Finite-Difference Time-Domain Method iterative model that relates to the factor that frequently becomes of power cable:
Because the decoupling zero of the cablebreak movable model under the time domain is comparatively complicated, therefore according to step 2 gained R, L, CThe distribution parameter matrix is set up Finite-Difference Time-Domain Method (FDTD) iterative model of cable wave equation, is illustrated in figure 2 as the equivalent-circuit model according to step 1 gained, along transmission direction according to D z/ 2 space interval disperses to electric current and voltage, and whole like this transmission line can be divided into 2 * NDZ section.On whole transmission line, voltage VWith electric current IThe D of being separated by z/ 2, total NDZ+1 of voltage discrete point, total NDZ of electric current discrete point, simultaneously, the time is with D t/ 2 step-length disperses.According to the method described above, the voltage of many conductor transmission line, current wave process become one group of space, temporal discrete point, and its time-space relationship as shown in Figure 3.
According to the single order central-difference formula, the cablebreak movable model under the time domain shown in formula (1) and the formula (2) is dispersed, obtain following Finite-Difference Time-Domain Method (FDTD) iterative model:
Figure 2011101035748100002DEST_PATH_IMAGE055
（11）
Figure 2011101035748100002DEST_PATH_IMAGE057
（12）
K=1 in the formula (11), 2 ... NDZ; K=2 in the formula (12), 3 ... NDZ.
Need to prove, in order to guarantee the stable of model, △ tAnd △ zValue to satisfy condition: △ t≤ △ z/ v, wherein, vThe max model speed of in many conductor transmission line, propagating for electromagnetic wave, vMode conversion in the available frequency domain analytic approach obtains.Above-mentioned iterative model has been considered the coupling situation of external electromagnetic field, but the not influence of CONSIDERING BOUNDARY CONDITIONS is therefore invalid to frontier point.
According to Kirchhoff's law, suppose no external electromagnetic field excitation, and do not consider loss, promptly R= G=0, discrete and arrangement can obtain Finite-Difference Time-Domain Method (FDTD) iterative model at node 1 and node NDZ+1 place respectively on time domain, i.e. the Finite-Difference Time-Domain Method iterative model of frontier point is specific as follows:
Figure 503688DEST_PATH_IMAGE021
（13）
（14）
Formula (13) is Finite-Difference Time-Domain Method (FDTD) iterative model at node 1 place, wherein,
Figure 392326DEST_PATH_IMAGE022
,
Figure 129338DEST_PATH_IMAGE023
Expression respectively
Figure 468702DEST_PATH_IMAGE006
,
Figure 960863DEST_PATH_IMAGE007
The voltage of moment power cable head end;
Figure 936910DEST_PATH_IMAGE024
Expression The electric current of moment power cable head end;
Figure 996450DEST_PATH_IMAGE025
,
Figure 331616DEST_PATH_IMAGE026
For
Figure 873587DEST_PATH_IMAGE006
,
Figure 14718DEST_PATH_IMAGE007
The voltage of moment signal source;
Figure 512696DEST_PATH_IMAGE027
Internal resistance for signal source;
Figure 831813DEST_PATH_IMAGE017
Be distributed capacitance;
Figure 110347DEST_PATH_IMAGE028
Be extra electric field intensity;
Figure 992853DEST_PATH_IMAGE018
Be the space interval of choosing;
Figure 155456DEST_PATH_IMAGE019
Be the time interval of choosing.
Formula (14) is Finite-Difference Time-Domain Method (FDTD) iterative model at node NDZ+1 place, wherein,
Figure 160321DEST_PATH_IMAGE031
,
Figure 598256DEST_PATH_IMAGE032
Expression respectively
Figure 831922DEST_PATH_IMAGE006
,
Figure 179858DEST_PATH_IMAGE007
The voltage of moment power cable end;
Figure 27728DEST_PATH_IMAGE033
Expression
Figure 15276DEST_PATH_IMAGE013
The electric current of moment power cable end;
Figure 239584DEST_PATH_IMAGE034
Internal impedance for load; Be distributed capacitance;
Figure 726377DEST_PATH_IMAGE028
Be extra electric field intensity;
Figure 138904DEST_PATH_IMAGE018
Be the space interval of choosing; Be the time interval of choosing.
With Finite-Difference Time-Domain Method iterative model lead-in cable volatility model, if do not consider Finite-Difference Time-Domain Method iterative model that parameter becomes factor frequently into:
Figure 972921DEST_PATH_IMAGE021
（15）
Figure 428173DEST_PATH_IMAGE030
（16）
Figure 141045DEST_PATH_IMAGE002
（17）
Figure 2011101035748100002DEST_PATH_IMAGE059
（18）
Because marked change can take place impedance parameter under the high frequency, do not consider that therefore the Finite-Difference Time-Domain Method iterative model of the factor that frequently becomes can not satisfy the needs of cable transient analysis.After considering the skin effect influence, by introducing the above-mentioned Finite-Difference Time-Domain Method electric current of time domain convolution correction iterative model
Figure 769473DEST_PATH_IMAGE003
, obtain after the correction as drag:
?（19）
Wherein, BThe influence of skin effect during for high frequency;
Figure 2011101035748100002DEST_PATH_IMAGE061
,
Figure 307081DEST_PATH_IMAGE062
The recursion item
Figure 2011101035748100002DEST_PATH_IMAGE063
, coefficient a i , a i Value as shown in table 1.
Table 1 coefficient a i , a i Value
n a i a i
1 0.07909818 -0.001148443
2 0.11543423 -0.013818329
3 0.13435380 -0.05437596
4 0.21870422 -0.14216494
5 0.09822967 -0.30128437
6 0.51360484 -0.5614219
7 -0.209629 -0.9711713
8 1.1974447 -1.6338433
9 0.01122549 -2.8951329
10 0.74425260 -5.0410969
Step 3 adopts the described Finite-Difference Time-Domain Method iterative model of step 2 to simulate the propagation characteristic of local discharge signal in power cable.
Adopt the inventive method that long single core 110kV copper core crosslinked polyethylene (XLPE) cable transient state of a 15m is simulated, above-mentioned cable data specifically is shown in Table 2.At first, according to cable data in the table 2, adopt R- L- CNetwork is set up the equivalent-circuit model of concentric cable under the high-frequency signal; Then, the employing finite element method calculates the equivalent-circuit model of cable R, L, CThe distribution parameter matrix; Then, according to R, L, CThe distribution parameter matrix is set up the Finite-Difference Time-Domain Method iterative model that relates to the factor that frequently becomes; At last, inject the simulation local discharge signal in cable head, selected simulation local discharge signal is two exponential waves, wherein wavefront/the wave rear of two exponential waves is 10ns/100ns, simulate the response signal (being voltage responsive) of end shielding layer according to the Finite-Difference Time-Domain Method iterative model that obtains, analog result as shown in Figure 5.
The basic parameter of table 2 twisted polyethylene cable
The parameter amount Numerical value The parameter amount Numerical value
Conductor nominal section/mm 2 240 Outer semiconductor thickness/mm 1.0
Conductor diameter/mm 18.40 Dredge around copper wire shielding cross section/mm 2 95
Interior semiconductor thickness/mm 1.5 Oversheath thickness/mm 3.4
Insulation thickness/mm 19.0 The outside diameter of cable/mm 76.5
In order to verify the accuracy of the inventive method, adopt connection line chart shown in Figure 4 to verify the inventive method.As shown in Figure 4, the impulse source that pumping signal takes place in the contrast test is to adopt high pressure nanosecond signal generator, and it can produce the simulation local discharge signal; Used oscillograph bandwidth 100MHz, sampling rate has high-speed sampling and memory function greater than 1GS/s; When impedance ground is 100 Ω.During test, simulation shelf depreciation waveform is elected two exponential waves as, wherein wavefront/the wave rear of two exponential waves is 10ns/100ns, employing sensitivity is the response signal that the High Frequency Current Sensor of 7.65mV/mA is measured cable first and last end shield layer respectively, and the response signal (being voltage responsive) of the end shielding layer that test records as shown in Figure 6.
Comparison diagram 5 and Fig. 6, result show that the simulation result that obtains according to the present invention can both be corresponding with test findings preferably in voltage responsive amplitude and time delay variation, thereby verified correctness of the present invention.

Claims (5)

1. a method of simulating local discharge signal propagation characteristic in power cable is characterized in that, may further comprise the steps successively:
Step 1 according to the electrical parameter of power cable, adopts R- L- CNetwork is set up the equivalent-circuit model of concentric cable under the high-frequency signal, and the calculating equivalent-circuit model R, L, CThe distribution parameter matrix;
Step 2, according to R, L, CThe distribution parameter matrix is set up the Finite-Difference Time-Domain Method iterative model that relates to the factor that frequently becomes of power cable:
Figure 2011101035748100001DEST_PATH_IMAGE001
Wherein,
Figure 398235DEST_PATH_IMAGE002
,
Figure 253059DEST_PATH_IMAGE003
Expression respectively
Figure 973409DEST_PATH_IMAGE004
,
Figure 368618DEST_PATH_IMAGE005
Moment power cable
Figure 120674DEST_PATH_IMAGE006
The voltage of position;
Figure 880819DEST_PATH_IMAGE007
Expression
Figure 85536DEST_PATH_IMAGE004
Moment power cable
Figure 284436DEST_PATH_IMAGE008
The voltage of position;
Figure 890998DEST_PATH_IMAGE009
,
Figure 822044DEST_PATH_IMAGE010
Expression respectively Moment power cable
Figure 251069DEST_PATH_IMAGE006
,
Figure 712137DEST_PATH_IMAGE012
The electric current of position;
Figure 76735DEST_PATH_IMAGE013
Expression
Figure 318360DEST_PATH_IMAGE014
Moment power cable
Figure 531167DEST_PATH_IMAGE006
The electric current of position; LBe distributed inductance;
Figure 112321DEST_PATH_IMAGE015
Be distributed capacitance; RBe resistance;
Figure 447487DEST_PATH_IMAGE016
Be the space interval of choosing;
Figure 114092DEST_PATH_IMAGE017
Be the time interval of choosing;
The Finite-Difference Time-Domain Method iterative model of power cable head end is:
Figure 130590DEST_PATH_IMAGE018
Wherein,
Figure 2011101035748100001DEST_PATH_IMAGE019
,
Figure 566250DEST_PATH_IMAGE020
Expression respectively
Figure 275580DEST_PATH_IMAGE004
,
Figure 163902DEST_PATH_IMAGE005
The voltage of moment power cable head end; Expression
Figure 351223DEST_PATH_IMAGE011
The electric current of moment power cable head end;
Figure 231454DEST_PATH_IMAGE022
,
Figure 669389DEST_PATH_IMAGE023
For
Figure 27689DEST_PATH_IMAGE004
,
Figure 437941DEST_PATH_IMAGE005
The voltage of moment signal source;
Figure 489074DEST_PATH_IMAGE024
Internal resistance for signal source;
Figure 351988DEST_PATH_IMAGE015
Be distributed capacitance;
Figure 576296DEST_PATH_IMAGE025
Be extra electric field intensity;
Figure 575476DEST_PATH_IMAGE016
Be the space interval of choosing;
Figure 797510DEST_PATH_IMAGE017
Be the time interval of choosing;
The Finite-Difference Time-Domain Method iterative model of power cable end is:
Figure 144790DEST_PATH_IMAGE026
Wherein,
Figure 172789DEST_PATH_IMAGE027
,
Figure 292055DEST_PATH_IMAGE028
Expression respectively
Figure 419411DEST_PATH_IMAGE004
, The voltage of moment power cable end;
Figure 823027DEST_PATH_IMAGE029
Expression The electric current of moment power cable end;
Figure 360636DEST_PATH_IMAGE030
Internal impedance for load;
Figure 747755DEST_PATH_IMAGE015
Be distributed capacitance; Be extra electric field intensity; Be the space interval of choosing; Be the time interval of choosing;
Step 3 adopts the described Finite-Difference Time-Domain Method iterative model of step 2 to simulate the propagation characteristic of local discharge signal in power cable.
2. the method for simulation local discharge signal according to claim 1 propagation characteristic in power cable is characterized in that: the described equivalent-circuit model of concentric cable of setting up under the high-frequency signal is to adopt many conductor transmission line method.
3. the method for simulation local discharge signal according to claim 1 and 2 propagation characteristic in power cable is characterized in that:
Described R, L, CThe distribution parameter matrix is to adopt finite element method to obtain.
4. the method for simulation local discharge signal according to claim 1 and 2 propagation characteristic in power cable is characterized in that:
Described
Figure 230022DEST_PATH_IMAGE013
Revise by introducing the time domain convolution.
5. the method for simulation local discharge signal according to claim 3 propagation characteristic in power cable is characterized in that:
Described Revise by introducing the time domain convolution.
CN 201110103574 2011-04-25 2011-04-25 Method for simulating transmission property of partial discharge signal in power cable Expired - Fee Related CN102156788B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110103574 CN102156788B (en) 2011-04-25 2011-04-25 Method for simulating transmission property of partial discharge signal in power cable

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110103574 CN102156788B (en) 2011-04-25 2011-04-25 Method for simulating transmission property of partial discharge signal in power cable

Publications (2)

Publication Number Publication Date
CN102156788A true CN102156788A (en) 2011-08-17
CN102156788B CN102156788B (en) 2013-03-20

Family

ID=44438283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110103574 Expired - Fee Related CN102156788B (en) 2011-04-25 2011-04-25 Method for simulating transmission property of partial discharge signal in power cable

Country Status (1)

Country Link
CN (1) CN102156788B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102707211A (en) * 2012-07-03 2012-10-03 广州供电局有限公司 Method and system for simulating ultrahigh-frequency radiation characteristic of cable straight joint
CN102722629A (en) * 2012-07-03 2012-10-10 广州供电局有限公司 Cable insulation connector ultrahigh frequency radiation characteristic simulation method and cable insulation connector terminal ultrahigh frequency radiation characteristic simulation system
CN102721909A (en) * 2012-07-03 2012-10-10 广州供电局有限公司 Ultrahigh frequency radiation characteristic simulation method and system for cable body
CN110442978A (en) * 2019-08-08 2019-11-12 华北电力大学（保定） A kind of more conductor distribution capacity quick calculation methods based on FInite Element
CN110516360A (en) * 2019-08-28 2019-11-29 哈尔滨工程大学 A kind of long line rapid simulation method based on FDTD
CN110888030A (en) * 2019-11-26 2020-03-17 深圳供电局有限公司 Short air gap discharge identification circuit, method and system of transformer
CN112836407A (en) * 2021-02-01 2021-05-25 中国科学技术大学 Simulation method for obtaining frequency response characteristic of logging cable
CN113805135A (en) * 2021-08-23 2021-12-17 深圳供电局有限公司 Method, apparatus, storage medium, and electronic apparatus for correcting partial discharge output signal

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201373905Y (en) * 2009-03-04 2009-12-30 厦门红相电力设备股份有限公司 Safety detecting evaluation system of power cable
CN101819235A (en) * 2010-01-07 2010-09-01 南京大学 Short-wave radio set electromagnetic pulse test circuit based on finite-difference time-domain analytical method
JP2010276420A (en) * 2009-05-27 2010-12-09 Mitsubishi Electric Corp Device and method for detecting partial discharge

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201373905Y (en) * 2009-03-04 2009-12-30 厦门红相电力设备股份有限公司 Safety detecting evaluation system of power cable
JP2010276420A (en) * 2009-05-27 2010-12-09 Mitsubishi Electric Corp Device and method for detecting partial discharge
CN101819235A (en) * 2010-01-07 2010-09-01 南京大学 Short-wave radio set electromagnetic pulse test circuit based on finite-difference time-domain analytical method

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
《中国优秀博硕士学位论文全文数据库(博士)工程科技II辑》 20070415 杜学峰。 长电缆传输研究 摘要、第二、三、四章 1-5 , *
《中国优秀硕士学位论文全文数据库工程科技II辑》 20090215 衣斌。 空间电磁脉冲作用下屏蔽电缆的耦合规律研究 摘要、第二、三章，图3-4 1-5 , *
《高电压技术》 20090731 唐炬，等。 高压电缆附件局部放电超高频检测与分析 第1571～1577页 1-5 第35卷, 第7期 *
杜学峰。: "长电缆传输研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技II辑》 *
衣斌。: "空间电磁脉冲作用下屏蔽电缆的耦合规律研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *
阮羚，等。: "基于多导体传输线模型的变压器绕组分布参数计算", 《高压电器》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102707211A (en) * 2012-07-03 2012-10-03 广州供电局有限公司 Method and system for simulating ultrahigh-frequency radiation characteristic of cable straight joint
CN102722629A (en) * 2012-07-03 2012-10-10 广州供电局有限公司 Cable insulation connector ultrahigh frequency radiation characteristic simulation method and cable insulation connector terminal ultrahigh frequency radiation characteristic simulation system
CN102721909A (en) * 2012-07-03 2012-10-10 广州供电局有限公司 Ultrahigh frequency radiation characteristic simulation method and system for cable body
CN102722629B (en) * 2012-07-03 2014-08-13 广州供电局有限公司 Cable insulation connector ultrahigh frequency radiation characteristic simulation method and cable insulation connector terminal ultrahigh frequency radiation characteristic simulation system
CN102707211B (en) * 2012-07-03 2014-09-17 广州供电局有限公司 Method and system for simulating ultrahigh-frequency radiation characteristic of cable straight joint
CN102721909B (en) * 2012-07-03 2015-07-08 广州供电局有限公司 Ultrahigh frequency radiation characteristic simulation method and system for cable body
CN110442978A (en) * 2019-08-08 2019-11-12 华北电力大学（保定） A kind of more conductor distribution capacity quick calculation methods based on FInite Element
CN110442978B (en) * 2019-08-08 2023-04-28 华北电力大学（保定） Rapid calculation method for multi-conductor distributed capacitance based on finite element method
CN110516360A (en) * 2019-08-28 2019-11-29 哈尔滨工程大学 A kind of long line rapid simulation method based on FDTD
CN110888030A (en) * 2019-11-26 2020-03-17 深圳供电局有限公司 Short air gap discharge identification circuit, method and system of transformer
CN110888030B (en) * 2019-11-26 2022-09-02 深圳供电局有限公司 Short air gap discharge identification circuit, method and system of transformer
CN112836407A (en) * 2021-02-01 2021-05-25 中国科学技术大学 Simulation method for obtaining frequency response characteristic of logging cable
CN113805135A (en) * 2021-08-23 2021-12-17 深圳供电局有限公司 Method, apparatus, storage medium, and electronic apparatus for correcting partial discharge output signal

Also Published As

Publication number Publication date
CN102156788B (en) 2013-03-20

Similar Documents

Publication Publication Date Title
CN102156788B (en) Method for simulating transmission property of partial discharge signal in power cable
Dubickas et al. High-frequency model of the Rogowski coil with a small number of turns
Visacro et al. Experimental impulse response of grounding grids
Li et al. Generalized equivalent cable bundle method for modeling EMC issues of complex cable bundle terminated in arbitrary loads
CN110688816A (en) Method for calculating time domain macro of nanosecond pulse to electrically large-size coil coupling response
Musznicki et al. Adaptive estimation of the transformer stray capacitances for DC–DC converter modelling
CN202402198U (en) Ignition energy detection device for magneto
Hosseini et al. Transformer winding modeling based on multi-conductor transmission line model for partial discharge study
Xie et al. Analysis of EMP coupling to a device from a wire penetrating a cavity aperture using transient electromagnetic topology
Cui et al. Circuit modeling of the test setup for pulsed current injection
CN104569668A (en) Method for designing Gaussian even-pulse high-current and high-power broadband injection probe
CN102645574B (en) On-site transient overvoltage measuring method based on network transfer function calculation
CN103197126A (en) Lightning induction overvoltage simulation test platform
Modarres et al. Design and implementation of a resistive MV voltage divider
Hou et al. A novel high frequency circuit model of power cable joint
CN203012053U (en) Device for testing port impedance characteristic of high voltage electrical equipment
Wenlu et al. Simulation analysis of crosstalk and radiation effects of cables inside aerocraft
Zhou et al. Lightning inductive coupling characteristics of parallel cables: Computation and measurement validation
Skibin et al. High frequency cable models for system level simulations in power electronics applications
Ye et al. Design and optimization of high frequency current transformer
Sun et al. Analysis of FDTD Crosstalk of Multiconductor Overhead Transmission Lines with a Frequency-Varying Ground Impedance Loss
Tikhomirova Electromagnetic transient modelling of grounding structures
Dou et al. Simulation of electromagnetic interference coupling to a substation secondary cable
Xianghui et al. Rogowski sensor for power grid traveling wave based fault location
CN202998015U (en) Resistor array voltage signal attenuator

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130320

Termination date: 20150425

EXPY Termination of patent right or utility model