CN107478954B - A kind of super extra-high voltage system corona loss calculation method based on distributed parameter model - Google Patents
A kind of super extra-high voltage system corona loss calculation method based on distributed parameter model Download PDFInfo
- Publication number
- CN107478954B CN107478954B CN201710705591.6A CN201710705591A CN107478954B CN 107478954 B CN107478954 B CN 107478954B CN 201710705591 A CN201710705591 A CN 201710705591A CN 107478954 B CN107478954 B CN 107478954B
- Authority
- CN
- China
- Prior art keywords
- corona loss
- route
- voltage
- transmission line
- line
- 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.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 21
- 230000005540 biological transmission Effects 0.000 claims abstract description 36
- 230000005611 electricity Effects 0.000 claims description 13
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 239000004020 conductor Substances 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 2
- 230000005684 electric field Effects 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 4
- 238000013461 design Methods 0.000 abstract description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
-
- 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
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of methods for calculating ultra-high/extra-high voltage system corona loss based on Transmission Line Distributed Parameter model.The π shape Equivalent Circuit Parameter calculation formula of ultra-high/extra-high voltage long range transmission line of alternation current is first derived by the distributed parameter model of long range transmission line of alternation current, then by accurately considering that the voltage's distribiuting on long transmission line derives the annual corona loss formula in route unit length, and then the corona loss of system is acquired by corona loss power equation.The corona loss value that the method is found out is of great significance to the design, operation and economic analysis of ultra-high/extra-high voltage system closer to actual measured value.
Description
Technical field
The super extra-high voltage system corona loss calculation method based on distributed parameter model that the present invention relates to a kind of, belongs to electric power
System line loss technical field.
Background technique
In modern power systems, ultra-high/extra-high voltage transmission line of electricity has become main force's power grid in China.For general voltage etc.
The transmission line of electricity of grade, can use lumped parameter equivalent circuit.But ultra-high/extra-high voltage Transmission Lines distance is up to hundreds of kilometer
Even thousands of kilometers, when studying the transmission characteristic of this route, it is necessary to distributed parameter model is used, at this point, such as reusing concentration
Parameter model just has very big deviation.
Since the voltage class of ultra-high/extra-high voltage transmission line of electricity is high, circuit surface is also easy to produce electric discharge phenomena, thus super/extra-high
The corona loss of pressure transmission line is very important sufficiently complex phenomenon.According to a large amount of measurement research, 500KV and Above Transmission Lines
Corona loss is very sensitive to route actual motion voltage, and the line is busy the sizable specific gravity of road total losses, thus in computing system
Corona loss should not be ignored when route total losses.
Summary of the invention
Purpose: in order to overcome deficiency present in existing calculation method, the present invention provides a kind of based on distributed parameter model
Super extra-high voltage system corona loss calculation method.
Technical solution: in order to solve the above technical problems, the technical solution adopted by the present invention are as follows:
A kind of super extra-high voltage system corona loss calculation method based on distributed parameter model, specifically comprises the following steps:
S01: the distributed parameter model based on long range transmission line of alternation current derives the long range ac transmission of excess of export extra-high voltage
The π shape Equivalent Circuit Parameter calculation formula of route；
S02: by accurately considering that the voltage's distribiuting on super extra-high voltage system long transmission line is derived in route unit length
Annual corona loss formula；
S03: Newton-Raphson iteration is carried out based on corona loss power equation, and then finds out the corona on every route
The corona loss of loss and system.
The calculation formula of each parameter of long transmission line equivalent circuit is derived in step S01:
Characteristic impedance based on super UHV transmission line:
Wherein:R, x and b is respectively transmission line of electricity unit length
Resistance, reactance and susceptance；
Line propagation constant:
Wherein:α, β are respectively attenuation coefficient and phase coefficient, and α
=| γ | cos θ, β=| γ | sin θ；
The accurate characteristics of distributed parameters for considering transmission line of electricity:
Wherein: l is the length of transmission line of electricity；
U=α l, v=β l is enabled, hyperbolic functions are utilized
By derivation and abbreviation, enable
It can must finally surpass the parameter calculation formula of UHV transmission line π shape equivalent circuit:
The annual corona loss formula in route unit length is derived in step S02:
Under the long wire model of super UHV transmission line, it is known that when line end voltage, electric current, along route arbitrary point
Voltage, electric current calculation formula:
Wherein:Respectively indicating away from line end length is voltage, the electric current at x；Respectively line end
Voltage, electric current；
The then voltage expression of super extra-high voltage system middle or long line road every bit:
Wherein:For route head end voltage；Z is the impedance of route；
The first-class day of unit length, snowy day and the corona loss calculation formula of rainy day are as follows:
Wherein: Pf、Ps、PrRespectively indicate the first-class day of unit length, snowy day and the corona loss of rainy day (kw/km)；reqIt indicates
Split conductor equivalent redius (cm)；EMFor conductive line surfaces maximum field strength (kv/cm)；E0For the critical electric field calculated with peek's formula
Intensity (kv/cm)；
The surface field strength of middle phase are as follows:
The surface field strength of side phase are as follows:
Wherein: n is split conductor radical；R ' is wire radius (cm)；D is division spacing (cm)；Dm is geometric mean distance
(cm)；U is actual motion voltage (kv)；
K2=1.03K1/1.1；
Peek's formula indicates are as follows:
Wherein: m1For the conductive line surfaces coefficient of roughness；m2For meteorological coefficient；δ is the relative density of air；
If setting in 1 year good day as Tf, snowy day Ts, rainy day Tr, while t+1 point is uniformly generated in the line, from
Every voltage of head end to end is followed successively by U1、U2、……、Ut+1, then the corona loss of route is t small line segment corona loss
Summation, small line segment actual motion voltage is expressed as the average voltage of this section of head end and end, the small line of n, m point-to-point transmission in route
The annual mean of corona loss in section unit length are as follows:
Wherein, Un、UmThe voltage of n, m two o'clock respectively in route；
The annual mean of corona loss so in route unit length are as follows:
The corona loss of the corona loss and system on every route is found out in step S03:
For the sake of to simplify the calculation, the corona loss of route is divided into the both ends that equal two parts are suspended on transmission line of electricity,
Be considered " burden with power ", derive in this way for a N node system, the power equation of node i are as follows:
It willYij=Gij+jBijIt substitutes into, obtains:
Wherein: PiFor the injection active power of PQ and PV node；QiFor the injection reactive power of PQ node；Gij、BijRespectively
For branch conductance and susceptance；θijPhase difference of voltage between node i and j；lijFor the length of route i-j；Ui、UjRespectively save
The voltage of point i, j；
Every route unit length annual corona loss value P of system can be obtained after ox daraf(reciprocal of farad) iterationcij, and then acquire every
Corona loss P on routecijlij, the corona loss of system can be obtained by finally adding up.
The utility model has the advantages that a kind of super extra-high voltage system corona loss calculating side based on distributed parameter model provided by the invention
Method, this calculation method accurately considers the every bit voltage on route, thus result is more accurate, to ultra-high/extra-high voltage system
Design and analysis are of great significance.
Detailed description of the invention
Fig. 1 is eight node system figure of example of the invention.
Specific embodiment
The present invention will be further explained with reference to the accompanying drawing.
As shown in Figure 1, eight node systems wherein balance nodes 5, PV node 6, remaining is all PQ node；System shares 5
Route, wherein the wire length of L1, L2, L3, L4 are 400km, and the wire length of L5 is 800km；Tf=7471h, Ts=56h, Tr=1233h,
Segments t=20, air relative density δ=1, conducting wire coefficient of roughness m1=0.9, m when good day2=1.0, m when snowy day2=
0.997, the m when rainy day2=0.993,500kv route model 4*LGJ-300, wire radius r '=1.213cm divide spacing d
=45cm, three-phase is horizontally arranged, phase spacing 13m.
The value of each parameter of long transmission line equivalent circuit is found out in step S01:
R=9.88572 Ω, the X=109.126 Ω, B=0.00321S, the R=of route L5 of route L1, L2, L3, L4
16.279 Ω, X=199.285 Ω, B=0.00673S.
The annual mean of corona loss in route unit length is found out in step S02:
The corona loss of the corona loss and system on every route is found out in step S03:
For the sake of to simplify the calculation, the corona loss of route is divided into the both ends that equal two parts are suspended on transmission line of electricity,
Be considered " burden with power ", derive in this way for a N node system, the power equation of node i are as follows:
It willYij=Gij+jBijIt substitutes into, obtains:
Wherein: PiFor the injection active power of PQ and PV node；QiFor the injection reactive power of PQ node；lijFor route i-
The length of j.
Every route unit length annual corona loss value P of system can be obtained after ox daraf(reciprocal of farad) iterationcij, and then acquire every
Corona loss P on routecijlij: the corona loss on route L1, L2 is 3.55693MW, the corona damage on route L3, L4
Consumption is 3.90702MW, and the corona loss on route L5 is 7.94192MW；The corona loss that system can finally be obtained is
22.8698MW.The above is only a preferred embodiment of the present invention, it should be pointed out that: for the ordinary skill of the art
For personnel, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications
It should be regarded as protection scope of the present invention.
Claims (2)
1. a kind of super extra-high voltage system corona loss calculation method based on distributed parameter model, it is characterised in that: including as follows
Step:
S01: the distributed parameter model based on long range transmission line of alternation current derives excess of export extra-high voltage long range transmission line of alternation current
π shape Equivalent Circuit Parameter calculation formula；
S02: by accurately considering that the voltage's distribiuting on super extra-high voltage system long transmission line derives that the year in route unit length flat
Equal corona loss formula；
S03: Newton-Raphson iteration is carried out based on corona loss power equation, and then finds out the corona loss on every route
With the corona loss of system；
The step S01 includes:
Characteristic impedance based on super UHV transmission line:
Wherein:R, x and b be respectively transmission line of electricity resistance per unit length,
Reactance and susceptance；
Line propagation constant:
Wherein:α, β are respectively attenuation coefficient and phase coefficient, and α=| γ
| cos θ, β=| γ | sin θ；
The accurate characteristics of distributed parameters for considering transmission line of electricity:
Wherein: l is the length of transmission line of electricity；
U=α l, v=β l is enabled, hyperbolic functions are utilizedThrough derivation and
Abbreviation enablesIt can must finally surpass UHV transmission line π shape
The parameter calculation formula of equivalent circuit:
The step S02 includes:
Under the long wire model of super UHV transmission line, it is known that electricity when line end voltage, electric current, along route arbitrary point
Pressure, electric current calculation formula:
Wherein:Respectively indicating away from line end length is voltage, the electric current at x；Respectively line end voltage,
Electric current；
The then voltage expression of super extra-high voltage system middle or long line road every bit:
Wherein:For route head end voltage；Z is the impedance of route；B is the susceptance of route；
The first-class day of unit length, snowy day and the corona loss calculation formula of rainy day are as follows:
Wherein: Pf、Ps、PrRespectively indicate the corona loss in the first-class day of unit length, snowy day and rainy day, unit kw/km；reqIt indicates
Split conductor equivalent redius, unit cm；EMFor conductive line surfaces maximum field strength, unit kv/cm；E0Face for what is calculated with peek's formula
Boundary's electric field strength, unit kv/cm；
The surface field strength of middle phase are as follows:
The surface field strength of side phase are as follows:
Wherein: n is split conductor radical；R ' is wire radius, unit cm；D is division spacing, unit cm；Dm is geometric mean distance,
Unit cm；U is actual motion voltage, unit kv；
K2=1.03K1/1.1；
Peek's formula indicates are as follows:
Wherein: m1For the conductive line surfaces coefficient of roughness；m2For meteorological coefficient；δ is the relative density of air；
If setting in 1 year good day as Tf, snowy day Ts, rainy day Tr, while t+1 point is uniformly generated in the line, from head end
Every voltage to end is followed successively by U1、U2、……、Ut+1, then the corona loss of route is the total of t small line segment corona loss
With small line segment actual motion voltage is expressed as the average voltage of this section of head end and end, the small line segment list of n, m point-to-point transmission in route
The annual mean of corona loss on bit length are as follows:
Wherein, Un、UmThe voltage of n, m two o'clock respectively in route；
The annual mean of corona loss so in route unit length are as follows:
2. a kind of super extra-high voltage system corona loss calculation method based on distributed parameter model according to claim 1,
It is characterized by: the step S03 includes:
The corona loss of route is divided into the both ends that equal two parts are suspended on transmission line of electricity, " burden with power " is set as, pushes away
Export is for a N node system, the power equation of node i are as follows:
It willYij=Gij+jBijIt substitutes into, obtains:
Wherein: PiFor the injection active power of PQ and PV node；QiFor the injection reactive power of PQ node；Gij、BijRespectively prop up
Road conductance and susceptance；θijPhase difference of voltage between node i and j；lijFor the length of route i-j；Ui、UjRespectively node i,
The voltage of j；
Every route unit length annual corona loss value P of system is obtained after Newton-Raphson approach iterationcij, and then acquire every
Corona loss P on routecijlij, the corona loss of system can be obtained by finally adding up.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710705591.6A CN107478954B (en) | 2017-08-16 | 2017-08-16 | A kind of super extra-high voltage system corona loss calculation method based on distributed parameter model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710705591.6A CN107478954B (en) | 2017-08-16 | 2017-08-16 | A kind of super extra-high voltage system corona loss calculation method based on distributed parameter model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107478954A CN107478954A (en) | 2017-12-15 |
| CN107478954B true CN107478954B (en) | 2019-08-06 |
Family
ID=60598923
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710705591.6A Expired - Fee Related CN107478954B (en) | 2017-08-16 | 2017-08-16 | A kind of super extra-high voltage system corona loss calculation method based on distributed parameter model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107478954B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109254189B (en) * | 2018-11-02 | 2021-10-08 | 江苏方天电力技术有限公司 | Submarine cable voltage withstand test terminal voltage calculation method |
| CN109948252B (en) * | 2019-03-20 | 2022-08-26 | 河海大学 | Ultra/extra-high voltage alternating current corona loss numerical value calculation method based on streamer theory |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101714234A (en) * | 2009-10-23 | 2010-05-26 | 西北电网有限公司 | Grid line loss monitoring and analyzing system |
| CN101794977A (en) * | 2010-03-18 | 2010-08-04 | 华北电力大学 | Long-distance transmission line deicing method and device |
| CN102545252A (en) * | 2011-12-22 | 2012-07-04 | 河海大学 | Voltage source commutation-high voltage direct current (VSC-HVDC) power flow computing method based on three-stage convergence Newton method |
| CN103902815A (en) * | 2014-03-10 | 2014-07-02 | 马善娟 | Method for building model of transmission line of medium voltage distribution network |
| CN104155543A (en) * | 2014-07-17 | 2014-11-19 | 国家电网公司 | Transmission line parameter identification method based on power loss method |
| CN104810827A (en) * | 2015-05-12 | 2015-07-29 | 华北电力大学 | Line optimum transmission power calculation method considering parallel connection reactor and corona loss |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9851411B2 (en) * | 2012-06-28 | 2017-12-26 | Keith S. Champlin | Suppressing HF cable oscillations during dynamic measurements of cells and batteries |
| US9568516B2 (en) * | 2014-09-23 | 2017-02-14 | Schweitzer Engineering Laboratories, Inc. | Determining status of electric power transmission lines in an electric power transmission system |
-
2017
- 2017-08-16 CN CN201710705591.6A patent/CN107478954B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101714234A (en) * | 2009-10-23 | 2010-05-26 | 西北电网有限公司 | Grid line loss monitoring and analyzing system |
| CN101794977A (en) * | 2010-03-18 | 2010-08-04 | 华北电力大学 | Long-distance transmission line deicing method and device |
| CN102545252A (en) * | 2011-12-22 | 2012-07-04 | 河海大学 | Voltage source commutation-high voltage direct current (VSC-HVDC) power flow computing method based on three-stage convergence Newton method |
| CN103902815A (en) * | 2014-03-10 | 2014-07-02 | 马善娟 | Method for building model of transmission line of medium voltage distribution network |
| CN104155543A (en) * | 2014-07-17 | 2014-11-19 | 国家电网公司 | Transmission line parameter identification method based on power loss method |
| CN104810827A (en) * | 2015-05-12 | 2015-07-29 | 华北电力大学 | Line optimum transmission power calculation method considering parallel connection reactor and corona loss |
Non-Patent Citations (2)
| Title |
|---|
| ±800kV特高压直流线路采用不均匀***降低电晕损耗;肖颍涛 等;《高压电器》;20130430;第14-19页 |
| 计入电晕损耗的潮流模型及计算;李劲波 等;《电力***自动化》;19970131;第33-35页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107478954A (en) | 2017-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107478954B (en) | A kind of super extra-high voltage system corona loss calculation method based on distributed parameter model | |
| CN105548802B (en) | A kind of asynchronous fault distance-finding method in three end of T link based on distribution character along fault traveling wave | |
| CN103941151A (en) | Radiation net fault location method by means of zero mode and aerial mode time difference independent of double-end synchronization and with matching of magnitude of voltages and magnitude of currents | |
| CN106019080B (en) | A kind of common-tower double-return DC line Single Terminal Traveling Wave Fault Location method based on energy jump along the line | |
| CN110531216B (en) | Overhead line and cable hybrid transmission line fault section distinguishing method | |
| CN103675565B (en) | Ultra-high/extra-high voltage transmission line of alternation current power frequency parameters measuring method | |
| CN104052053B (en) | Based on the power distribution network double line terminal section impedance analysis method of Linear Regression Model in One Unknown | |
| CN108614190A (en) | Power distribution network T link method for parameter estimation based on μ PMU and AMI | |
| CN110011314A (en) | Tidal current computing method and system based on transmission line of electricity electro thermal coupling effect | |
| CN105071447A (en) | Alternating current/direct-current hybrid connection multi-infeed direct-current receiving end power grid strength obtaining method | |
| CN105676015B (en) | A kind of transmission line of electricity carrying current calculation method | |
| CN110333394A (en) | A kind of low-voltage network line impedance estimation method | |
| CN105021903B (en) | The measuring method of nearby HV Transmission Line electric-field intensity on a kind of hillside | |
| CN108362978A (en) | A kind of T-type power transmission line fault locating method based on range function phase characteristic | |
| CN111123027A (en) | Multi-end transmission line fault location method based on WAMS system | |
| CN108984966A (en) | It is a kind of meter and route of power transmission electro thermal coupling effect wind power grid Limit Calculation Method | |
| CN105911429A (en) | Particle swarm optimization based aerial conductor's double end out-synchronization fault location method | |
| CN105548783B (en) | Ultrahigh-voltage alternating-current transmission line with four-circuit on single tower Zero sequence parameter accurate measurement method | |
| CN105158555B (en) | A kind of distribution line booste operation saves measurement and the verification method of quantity of electricity | |
| CN110046322A (en) | A kind of long cable voltage test core voltage Analytic Calculation Method | |
| CN105548803B (en) | UHV transmission line Fault Locating Method based on weighted least square algorithm | |
| KR101205844B1 (en) | On-line calculation method of available capacity of transmission line | |
| Huang et al. | Optimization research on ampacity of underground high voltage cable based on interior point method | |
| Wang et al. | 3-D electric field computation with charge simulation method around buildings near HV transmission lines | |
| Chen et al. | Thermal grading around overhead transmission line under various environments and its influence to load capacity |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190806 |