CN107328995A - A kind of high-voltage dc transmission electrical impedance scanning method - Google Patents
A kind of high-voltage dc transmission electrical impedance scanning method Download PDFInfo
- Publication number
- CN107328995A CN107328995A CN201710585634.1A CN201710585634A CN107328995A CN 107328995 A CN107328995 A CN 107328995A CN 201710585634 A CN201710585634 A CN 201710585634A CN 107328995 A CN107328995 A CN 107328995A
- Authority
- CN
- China
- Prior art keywords
- node
- voltage
- admittance
- scanning method
- files
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
Comprise the following steps the invention discloses a kind of high-voltage dc transmission electrical impedance scanning method：1) it is excel files by .dat files and .swi file translations；2) node to be measured in setting impedance scan scope and high-voltage dc transmission electric network；3) fundamental frequency model of transmission line, transformer, generator, load is converted into corresponding high frequency model；4) admittance matrix Yn；5) using Sparse technology storage node admittance matrix Yn；6) harmonic impedance is solved；The admittance of node to be measured is placed on last column of admittance matrix by the high-voltage dc transmission electrical impedance scanning method, saves backward steps, simplifies calculating process, while .dat files and .swi files can be received directly；In addition, the software interface close friend developed by the high-voltage dc transmission electrical impedance scanning method, it is easy to operate, it is easy to learn.
Description
Technical field
The present invention relates to technical field of HVDC transmission, more particularly to a kind of high-voltage dc transmission electrical impedance scanning method.
Background technology
Domestic designing unit generally carries out AC system harmonic impedance equivalence calculating using NIMSCAN softwares, but
NIMSCAN softwares have as a drawback that；
1) impedance scan function is single, and NIMSCAN softwares can only read the file of .nim forms in itself, needs and third party connects
Mouth program is used cooperatively, and is .nim files by .dat files and .swi file translations, otherwise when forming new .nim files, is needed
By set form .nim files are filled in each original document；And the quantity of documents scanned every time is huge, to staff with
Carry out very big workload, influence operating efficiency, and excessive troublesome operation may make staff's manual operations error, lead
Cause corrupt data.The minimax impedance Data-Statistics of multiple .out files are needed by one .nim files one .out file of correspondence
User oneself is wanted to handle or completed by third party software；NIMSCAN softwares do not generate system impedance figure, amplitude versus frequency characte in itself yet
Figure and phase-frequency characteristic figure.
2) it is restricted to scanning nodal point number, 12000 nodes of most scannings, and the nodal point number of the BPA files of real system
Huge, user need to first carry out equivalence, the big complex operation of workload in BPA to the distant place node overruned；
3) interface is unfriendly, in-convenience in use, NIMSCAN is based on DOS environment exploitations, no operation interface, knot to be measured
The input such as period information is required to write in .nim input files, it is desirable to operating personnel to the data modes of .nim files very
It is familiar with.
The content of the invention
It is an object of the invention to provide a kind of high-voltage dc transmission electrical impedance scanning method.
Therefore, technical solution of the present invention is as follows：
A kind of high-voltage dc transmission electrical impedance scanning method comprises the following steps：
1) will be comprising the trend and the .dat files of impedance perunit value and .swi file translations under each force device fundamental frequency
Node information and force device parameter information are included in excel files, excel files；
2) node to be measured in setting impedance scan scope and high-voltage dc transmission electric network；
3) fundamental frequency model of transmission line, transformer, generator, load is converted into corresponding high frequency model, and according to height
Frequency model inference goes out the relational expression I of admittance and voltage x currentn=YnUn；
4) information in excel files is read, according to step 3) in high frequency model calculate admittance between each node, it is and raw
Into the matrix of node admittance Yn,And the admittance of node to be measured is arranged on admittance matrix YnMost end
OK；
5) using Sparse technology storage node admittance matrix Yn；
6) to formula In=YnUnIn current vector InAssignment, puts 1, its aftercurrent by most footline node row electric current to be measured
Set to 0, the U tried to achievenThe harmonic impedance of node as to be measured.
Preferably, step 5) in admittance matrix YnIn storage, by admittance matrix YnDiagonal element leave concentratedly one-dimensional
In array, the nonzero element in its upper triangular matrix is stored in another one-dimension array successively by row, then correspondingly with two
One-dimension array indicates row, column position of each element in original matrix.
Preferably, the step 6) middle solution UnShi Caiyong Gaussian elimination methods.
Compared with prior art, the admittance of node to be measured is placed on admittance matrix by the high-voltage dc transmission electrical impedance scanning method
Last column, save backward steps, simplify calculating process, while .dat files and .swi files can be received directly；Remove
Outside this, the software interface developed by the high-voltage dc transmission electrical impedance scanning method is friendly, easy to operate, easy to learn.
Brief description of the drawings
The flow chart for the high-voltage dc transmission electrical impedance scanning method that Fig. 1 provides for the present invention.
Fig. 2 is generator Equivalent Model.
Fig. 3 is transformer Equivalent Model.
Fig. 4 is the π type Equivalent Models of two-winding transformer.
Fig. 5 is transmission line of electricity Equivalent Model.
Fig. 6 is series load model.
Fig. 7 is the HVDC transmission line model in embodiment.
Fig. 8 is input parameter excel tables.
Fig. 9 is the amplitude versus frequency characte figure under PSCAD is verified.
Figure 10 is the phase-frequency characteristic figure under PSCAD is verified.
Figure 11 is system impedance administrative division map.
Embodiment
Below in conjunction with the accompanying drawings and specific embodiment the present invention is described further, but following embodiments are absolutely not to this hair
It is bright to have any limitation.
As shown in figure 1, the high-voltage dc transmission electrical impedance scanning method comprises the following steps：1) each force device base will be included
The .dat files and .swi file translations of the lower trend of frequency and impedance perunit value are comprising tying in excel files, excel files
Point information and force device parameter information；2) node to be measured in setting impedance scan scope and HVDC transmission system；3)
The fundamental frequency model of transmission line, transformer, generator, load is converted into corresponding high frequency model；4) read in excel files
Information, according to step 3) in high frequency model calculate admittance between each node, generation the matrix of node admittance Yn,And the admittance of node to be measured is arranged on admittance matrix YnMost footline；5) deposited using Sparse technology
Store up the matrix of node admittance Yn；6) to admittance calculation formula In=YnUnIn current vector InAssignment, will most footline node row to be measured
Electric current puts 1, and its aftercurrent is set to 0, and U is solved using Gaussian elimination methodn, the U tried to achievenAs node to be measured harmonic impedance.
Admittance matrix is left the diagonal element of admittance matrix in one-dimension array concentratedly, thereon three angular moments in storage
Nonzero element in battle array is stored in another one-dimension array successively by row, then correspondingly with two one-dimension arrays indicates each member
Row, column position of the element in original matrix.
Fig. 2 is the Equivalent Model of generator, shown in the value such as formula (1) of resistance reactance, wherein, R (1) takes generator
Armature resistance Ra；X (1) cut-off axle subtranient reactance Xd”；H is overtone order；
Fig. 3 and Fig. 4 are transformer Equivalent Model, and parameter chooses such as formula (2) and formula (3).
Wherein, R (1) takes transformer fundamental frequency resistance RT；X (1) takes transformer fundamental frequency resistance XT。
Need progress π type equivalences to be calculated again on this basis, see Fig. 4.
Reduction is to each equivalent branch road conductance such as following formula of i sides in Fig. 4, Z in formulaT(h)=R (h)+jX (h), k are i sides to j
The voltage ratio of side, i.e. k=Ui/Uj。
Three-winding transformer is both needed to be converted into two two-winding transformer inputs, the present invention and turned with BPA systems in BPA
The excel files of change are not required to consider three-winding transformer as input.
Especially, the electric power networks containing transformer are calculated, k=Ui/UjWhen, calculated by formula (3), each branch admittance is with k
It is worth on the basis of the voltage of side, if measuring point is not at k, the perunit value conversion between different voltage class need to be carried out.
Fig. 4 is transmission line of electricity Equivalent Model, wherein shown in parameter value such as formula (4),
R, X, B therein take the famous value of the circuit lumped parameter under the fundamental frequency read in from BPA.
Fig. 5 is series load model,
1) load is perception
Work as Q>When 0, load is perception, shown in equivalent impedance such as formula (5),
2) load is capacitive
Work as Q<When 0, load is capacitive, shown in equivalent impedance such as formula (6),
H subharmonic impedance values ZhMethod for solving such as formula (7) shown in,
When electric current is unit harmonic current value, UhThe value showed is equal to harmonic wave equivalent impedance Z hereinh, resistance
The inverse of anti-value is admittance；
The high-voltage dc transmission electrical impedance scanning method provided using PSCAD softwares the present invention verifies that Fig. 7 is to test
Demonstrate,prove the circuit diagram of the high-voltage dc transmission power network used；Fig. 8 be BPA data conversions into input parameter excel forms；Set
Range of scanned frequencies is 50~2500Hz, and node to be measured is node 5；Model is built in PSCAD softwares and carries out simulating, verifying, is adjusted
With PSCAD harmonic impedance scan module, acquiescence output .out formal files, the data file of output can be by text document
Open.
According to the demand of impedance scan, circuit pack uses Bei Jielong models in PSCAD emulation；In addition, load is using permanent
Impedance model, generator uses normal impedance power source model.
Amplitude versus frequency characte figure, phase-frequency characteristic figure and the system impedance figure of result of calculation respectively as shown in Fig. 9, Figure 10 and Figure 11,
Plus sige "+" is program calculated result in figure, and circle is PSCAD simulation results, and its result and this software result of calculation are basically identical.
This embodiment demonstrates the correctness of the method for the invention, and the impedance scan available for actual DC engineering is designed
In.
Claims (3)
1. a kind of high-voltage dc transmission electrical impedance scanning method, it is characterised in that comprise the following steps：
1) will be excel comprising the trend and the .dat files of impedance perunit value and .swi file translations under each force device fundamental frequency
Node information and force device parameter information are included in file, excel files；
2) node to be measured in setting impedance scan scope and high-voltage dc transmission electric network；
3) fundamental frequency model of transmission line, transformer, generator, load is converted into corresponding high frequency model；
4) information in excel files is read, according to step 3) in high frequency model calculate admittance between each node, and generate knot
Point admittance matrix Yn,And the admittance of node to be measured is arranged on admittance matrix YnMost footline；
5) using Sparse technology storage node admittance matrix Yn；
6) to formula In=YnUnIn current vector InAssignment, puts 1, its aftercurrent is set to 0, and is asked by most footline node row electric current to be measured
The U obtainednThe harmonic impedance of node as to be measured.
2. high-voltage dc transmission electrical impedance scanning method according to claim 1, it is characterised in that step 5) in admittance matrix
YnIn storage, by admittance matrix YnDiagonal element leave concentratedly in one-dimension array, the nonzero element in its upper triangular matrix
It is stored in successively in another one-dimension array by row, then correspondingly with two one-dimension arrays indicate each element in original matrix
Row, column position.
3. high-voltage dc transmission electrical impedance scanning method according to claim 1 or 2, it is characterised in that the step 6) in
Solve UnShi Caiyong Gaussian elimination methods.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710585634.1A CN107328995A (en) | 2017-07-18 | 2017-07-18 | A kind of high-voltage dc transmission electrical impedance scanning method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710585634.1A CN107328995A (en) | 2017-07-18 | 2017-07-18 | A kind of high-voltage dc transmission electrical impedance scanning method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107328995A true CN107328995A (en) | 2017-11-07 |
Family
ID=60226449
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710585634.1A Pending CN107328995A (en) | 2017-07-18 | 2017-07-18 | A kind of high-voltage dc transmission electrical impedance scanning method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107328995A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109038534A (en) * | 2018-07-17 | 2018-12-18 | 中国南方电网有限责任公司 | A kind of novel DC transmission system harmonic impedance scanning means |
| CN110717145A (en) * | 2019-09-19 | 2020-01-21 | 南昌大学 | Piecewise symmetric reverse Gaussian-approximately-equal-elimination method based on symmetric sparse matrix technology |
| CN112231926A (en) * | 2020-10-29 | 2021-01-15 | 高小翎 | Node-optimized harmonic impedance measurement method for large power system |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101282042A (en) * | 2007-12-24 | 2008-10-08 | 国网直流工程建设有限公司 | Impedance frequency eqmivalent value novel method for designing AC filter |
| CN102914697A (en) * | 2012-10-30 | 2013-02-06 | 西安交通大学 | Micro-grid harmonic impedance measuring method based on three-phase symmetrical square wave current injection |
| CN103018534A (en) * | 2012-12-20 | 2013-04-03 | 广东电网公司电力科学研究院 | Method and system for determining harmonic voltage |
| CN103544378A (en) * | 2013-09-28 | 2014-01-29 | 南方电网科学研究院有限责任公司 | Method for calculating harmonic wave impedance of alternating current system for direct current transmission |
| CN104063519A (en) * | 2014-07-16 | 2014-09-24 | 国家电网公司 | BPA power grid data analyzing and managing method and system based on EXCEL |
| CN105335428A (en) * | 2014-08-08 | 2016-02-17 | 国家电网公司 | Data conversion method and apparatus |
| CN105529711A (en) * | 2016-01-28 | 2016-04-27 | 国家电网公司 | BPA data based alternating current system harmonic impedance scanning method |
-
2017
- 2017-07-18 CN CN201710585634.1A patent/CN107328995A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101282042A (en) * | 2007-12-24 | 2008-10-08 | 国网直流工程建设有限公司 | Impedance frequency eqmivalent value novel method for designing AC filter |
| CN102914697A (en) * | 2012-10-30 | 2013-02-06 | 西安交通大学 | Micro-grid harmonic impedance measuring method based on three-phase symmetrical square wave current injection |
| CN103018534A (en) * | 2012-12-20 | 2013-04-03 | 广东电网公司电力科学研究院 | Method and system for determining harmonic voltage |
| CN103544378A (en) * | 2013-09-28 | 2014-01-29 | 南方电网科学研究院有限责任公司 | Method for calculating harmonic wave impedance of alternating current system for direct current transmission |
| CN104063519A (en) * | 2014-07-16 | 2014-09-24 | 国家电网公司 | BPA power grid data analyzing and managing method and system based on EXCEL |
| CN105335428A (en) * | 2014-08-08 | 2016-02-17 | 国家电网公司 | Data conversion method and apparatus |
| CN105529711A (en) * | 2016-01-28 | 2016-04-27 | 国家电网公司 | BPA data based alternating current system harmonic impedance scanning method |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109038534A (en) * | 2018-07-17 | 2018-12-18 | 中国南方电网有限责任公司 | A kind of novel DC transmission system harmonic impedance scanning means |
| CN110717145A (en) * | 2019-09-19 | 2020-01-21 | 南昌大学 | Piecewise symmetric reverse Gaussian-approximately-equal-elimination method based on symmetric sparse matrix technology |
| CN112231926A (en) * | 2020-10-29 | 2021-01-15 | 高小翎 | Node-optimized harmonic impedance measurement method for large power system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lei et al. | A general unified AC/DC power flow algorithm with MTDC | |
| Orillaza et al. | Harmonic state-space model of a controlled TCR | |
| Dugan et al. | An open source platform for collaborating on smart grid research | |
| CN107328995A (en) | A kind of high-voltage dc transmission electrical impedance scanning method | |
| CN106655195B (en) | Active power distribution network high-frequency harmonic power flow calculation method | |
| Semlyen | S-domain methodology for assessing the small signal stability of complex systems in nonsinusoidal steady state | |
| CN106845041B (en) | MMC-based real-time simulation system, simulation method and MMC valve simulator | |
| Abbood et al. | Data‐Driven Modeling of a Commercial Photovoltaic Microinverter | |
| CN106159955A (en) | Based on the power system distributed optimal power flow method punishing Duality Decomposition continuously | |
| CN111628504A (en) | Inverter modeling method containing amplitude limiter and related device | |
| Mazhari et al. | Interfacing transient stability and extended harmonic domain for dynamic harmonic analysis of power systems | |
| Ou et al. | Research and application of small time-step simulation for MMC VSC-HVDC in RTDS | |
| von Jouanne et al. | Power hardware-in-the-loop (PHIL): A review to advance smart inverter-based grid-edge solutions | |
| Zhang et al. | Self‐adaptable reactive power‐voltage controller for virtual synchronous generators | |
| Gao et al. | Real-time model of the simplified australian 14-generator system in hypersim | |
| Eminoglu | A new model for wind turbine systems | |
| CN106953331A (en) | A kind of Harmonic Power Flow of Power Systems algorithm for considering uncertain and three-phase imbalance | |
| Williams et al. | Increasing green energy penetration and efficient utilization through a finite control set model predictive-based virtual inductance droop control | |
| CN110378020A (en) | Line commutation inverter multiband Dynamic Phasors electromagnetical transient emulation method and system | |
| CN109309385A (en) | Hybrid active filter Optimal Configuration Method in a kind of active power distribution network | |
| Rehman et al. | AC/DC multi‐infeed power flow solution | |
| CN102540907B (en) | Parallel digital and analog integrated simulation interface and physical simulation subsystem interface | |
| CN107563062A (en) | A kind of AC and DC hybrid transmission system modeling method based on Dynamic Phasors | |
| Tang et al. | Application of hardware‐in‐loop in teaching power electronic course based on a low‐cost platform | |
| Fetzer | Development of a MATLAB/Simulink Framework for Phasor-Based Power System Simulation and Component Modeling Based on State Machines |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171107 |