CN107968402A - A kind of localization method of more current source type oscillation sources - Google Patents

A kind of localization method of more current source type oscillation sources Download PDF

Info

Publication number
CN107968402A
CN107968402A CN201711309964.4A CN201711309964A CN107968402A CN 107968402 A CN107968402 A CN 107968402A CN 201711309964 A CN201711309964 A CN 201711309964A CN 107968402 A CN107968402 A CN 107968402A
Authority
CN
China
Prior art keywords
node
msub
voltage
frequency
mrow
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
Application number
CN201711309964.4A
Other languages
Chinese (zh)
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.)
China Electric Power Research Institute Co Ltd CEPRI
North China Electric Power University
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
North China Electric Power University
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 China Electric Power Research Institute Co Ltd CEPRI, North China Electric Power University filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201711309964.4A priority Critical patent/CN107968402A/en
Publication of CN107968402A publication Critical patent/CN107968402A/en
Pending legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a kind of localization method of more current source type oscillation sources, the emulation data or actual operating data of each node voltage in electric system to be detected are obtained first, and obtain the frequency distribution of each node voltage using Fourier transformation according to the data obtained；Based on the frequency distribution, TLS ESPRIT analyses are carried out to each node voltage, obtain voltage magnitude, frequency and the phase of each node voltage at different frequencies；The line parameter circuit value of electric system to be detected is obtained, different frequency lower node admittance matrix is obtained according to the line parameter circuit value obtained；According to the node voltage and bus admittance matrix under the different frequency obtained, ask for the Injection Current of each node by nodal voltage equation, and more current source type oscillation sources are positioned with this.Above method calculation amount is small, can efficiently navigate to more current source type oscillation sources, and fault rate is low, is with a wide range of applications.

Description

A kind of localization method of more current source type oscillation sources
Technical field
The present invention relates to power grid security technical field, more particularly to a kind of localization method of more current source type oscillation sources.
Background technology
As electric power networks interconnection degree gradually steps up, system structure becomes more and more huger, each region connects each other, Each several part is closely related, interacts in the process of running so that the method for operation also becomes increasingly complex, and ensures that system is safe and reliable The difficulty of operation is gradually increased, and the power oscillation problem in power transmission process between system shows more and more prominent.It is this to shake Swing certain generator in often system and receive periodic disturbance starting of oscillation first, system is traveled to rapidly then along electric power networks In other nodes, wherein disturbing source can be equivalent to current source.
Prior art has focused largely on Energy Analysis for High to the localization method of this disturbing source, by judging certain biography Location perturbation source is come in the energy transmission direction of defeated line or certain node, although the disturbance source locating method based on energy spectrometer can Effectively positioning oscillation source, but computationally intensive, especially in the system of multiple disturbing source, such method is difficult to be accurately positioned to shaking Source is swung, or even the situation of positioning by mistake occurs.
The content of the invention
The object of the present invention is to provide a kind of localization method of more current source type oscillation sources, this method calculation amount is small, can be high Effect navigates to more current source type oscillation sources, and fault rate is low.
The purpose of the present invention is what is be achieved through the following technical solutions：
A kind of localization method of more current source type oscillation sources, the described method includes：
Step 1, the emulation data or actual operating data for obtaining each node voltage in electric system to be detected, and according to institute The data of acquisition obtain the frequency distribution of each node voltage using Fourier transformation；
Step 2, based on the frequency distribution, TLS-ESPRIT analyses are carried out to each node voltage, obtain each node voltage Voltage magnitude, frequency and phase at different frequencies；
Step 3, the line parameter circuit value for obtaining electric system to be detected, are obtained under different frequency according to the line parameter circuit value obtained Bus admittance matrix；
Step 4, according to the node voltage and bus admittance matrix under the different frequency obtained, pass through nodal voltage equation Ask for the Injection Current of each node, and more current source type oscillation sources are positioned with this.
As seen from the above technical solution provided by the invention, above method calculation amount is small, can efficiently navigate to more Current source type oscillation source, and fault rate is low, is with a wide range of applications.
Brief description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the present invention, for this For the those of ordinary skill in field, without creative efforts, other can also be obtained according to these attached drawings Attached drawing.
Fig. 1 is the localization method flow diagram of more current source type oscillation sources provided in an embodiment of the present invention；
The system construction drawing for the four machine system examples that Fig. 2 is enumerated by the embodiment of the present invention.
Embodiment
With reference to the attached drawing in the embodiment of the present invention, the technical solution in the embodiment of the present invention is carried out clear, complete Ground describes, it is clear that described embodiment is only part of the embodiment of the present invention, instead of all the embodiments.Based on this The embodiment of invention, the every other implementation that those of ordinary skill in the art are obtained without making creative work Example, belongs to protection scope of the present invention.
The embodiment of the present invention is described in further detail below in conjunction with attached drawing, is implemented as shown in Figure 1 for the present invention The localization method flow diagram for more current source type oscillation sources that example provides, the described method includes：
Step 1, the emulation data or actual operating data for obtaining each node voltage in electric system to be detected, and according to institute The data of acquisition obtain the frequency distribution of each node voltage using Fourier transformation；
In the step 1, in the frequency distribution of each node voltage obtained：
In addition to the fundamental component of 50Hz, other existing frequency components are expressed as：fk(k=1,2,3 ...).
For specific example, the system knot for the four machine system examples that the embodiment of the present invention is enumerated is illustrated in figure 2 Composition, with reference to figure 2, applies the current source type oscillation disturbances that frequency is 0.06Hz on the generator being connected with node 2, with section Apply the current source type oscillation disturbances that frequency is 40.25Hz on the connected generator of point 4.
The emulation data or actual operating data of each node voltage in the system are obtained first, then pass through Fourier transformation The frequency distribution for obtaining each node voltage is：50Hz, 40.25Hz, 0.06Hz, then it is existing in addition to the fundamental component of 50Hz Other frequency components are expressed as f1=40Hz, f2=0.06Hz.
Step 2, based on the frequency distribution, TLS-ESPRIT analyses are carried out to each node voltage, obtain each node voltage Voltage magnitude, frequency and phase at different frequencies；
In the step 2, detailed process is：
The frequency distribution is primarily based on, TLS-ESPRIT analyses are carried out to each node voltage, each node voltage is obtained and exists Voltage magnitude U under different frequencymk, frequency fk, phaseThen node i is in frequency fkComponent of voltage is expressed as：
Wherein, ωk=2 π fk；T is the moment；
Further being write above formula as vector form is：
Write each node voltage vector as column vector again, obtain the node voltage under different frequency
For the example of Fig. 2：
Node i is in frequency f1Component of voltage be：
Wherein, ω1=2 π f1；
Node i is in frequency f2Component of voltage be：
Wherein, ω2=2 π f2；
Then write as column vector, the node voltage obtained under different frequency is expressed as：
Step 3, the line parameter circuit value for obtaining electric system to be detected, are obtained under different frequency according to the line parameter circuit value obtained Bus admittance matrix；
In this step, the process that different frequency lower node admittance matrix is obtained according to the line parameter circuit value obtained is specific For：
The line parameter circuit value obtained includes resistance R, inductance L, direct-to-ground capacitance C, asks for first under different frequency between i, j node Line impedance Zij, then：
It is again π type equivalent circuits by line equivalent, obtains the bus admittance matrix Y under different frequencyk。
Or by taking the four machine systems of Fig. 2 as an example, the bus admittance matrix Y under final different frequencykIt is expressed as：
Step 4, according to the node voltage and bus admittance matrix under the different frequency obtained, pass through nodal voltage equation Ask for the Injection Current of each node, and more current source type oscillation sources are positioned with this.
In the step 4, the node voltage under obtained different frequency is expressed asNode admittance under different frequency Matrix is expressed as Yk, then the Injection Current of each node is asked for by nodal voltage equationIt is expressed as：
By taking the four machine systems of Fig. 2 as an example, the Injection Current of each node is expressed as：
Based on the expression formula, there is the node that electric current injects to be positioned as more current source type oscillation sources, therefore have frequency at node 2 Injected for the electric current of 40.25Hz, there is the electric current that frequency is 0.06Hz to inject at node 4, thus can position generator 2 and 4 to shake Swing source.
It is worth noting that, the content not being described in detail in the embodiment of the present invention belongs to professional and technical personnel in the field's public affairs The prior art known.
The foregoing is only a preferred embodiment of the present invention, but protection scope of the present invention be not limited thereto, Any one skilled in the art is in the technical scope of present disclosure, the change or replacement that can readily occur in, It should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with the protection model of claims Subject to enclosing.

Claims (5)

  1. A kind of 1. localization method of more current source type oscillation sources, it is characterised in that the described method includes：
    Step 1, the emulation data or actual operating data for obtaining each node voltage in electric system to be detected, and according to being obtained Data the frequency distribution of each node voltage is obtained using Fourier transformation；
    Step 2, based on the frequency distribution, TLS-ESPRIT analyses are carried out to each node voltage, obtain each node voltage not Voltage magnitude, frequency and phase under same frequency；
    Step 3, the line parameter circuit value for obtaining electric system to be detected, different frequency lower node is obtained according to the line parameter circuit value obtained Admittance matrix；
    Step 4, according to the node voltage and bus admittance matrix under the different frequency obtained, asked for by nodal voltage equation The Injection Current of each node, and more current source type oscillation sources are positioned with this.
  2. 2. the localization method of more current source type oscillation sources according to claim 1, it is characterised in that in the step 1, institute In the frequency distribution of each node voltage obtained：
    In addition to the fundamental component of 50Hz, other existing frequency components are expressed as：fk(k=1,2,3 ...).
  3. 3. the localization method of more current source type oscillation sources according to claim 1, it is characterised in that in the step 2, institute State and TLS-ESPRIT analyses carried out to each node voltage, obtain each node voltage voltage magnitude at different frequencies, frequency and The process of phase is specially：
    The frequency distribution is primarily based on, TLS-ESPRIT analyses are carried out to each node voltage, obtain each node voltage in difference Voltage magnitude U under frequencymk, frequency fk, phaseThen node i is in frequency fkComponent of voltage is expressed as：
    Wherein, ωk=2 π fk；T is the moment；
    Further being write above formula as vector form is：
    Write each node voltage vector as column vector again, obtain the node voltage under different frequency
  4. 4. the localization method of more current source type oscillation sources according to claim 1, it is characterised in that in the step 3, institute State and the process of different frequency lower node admittance matrix is obtained according to the line parameter circuit value obtained be specially：
    The line parameter circuit value obtained includes resistance R, inductance L, direct-to-ground capacitance C, asks for the line between i, j node under different frequency first Road impedance Zij, then：
    <mrow> <msub> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>j&amp;omega;</mi> <mi>k</mi> </msub> <msub> <mi>L</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mrow> <msub> <mi>j&amp;omega;</mi> <mi>k</mi> </msub> <msub> <mi>C</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> </mfrac> </mrow>
    It is again π type equivalent circuits by line equivalent, obtains the bus admittance matrix Y under different frequencyk。
  5. 5. the localization method of more current source type oscillation sources according to claim 1, it is characterised in that in the step 4, institute The Injection Current that each node is asked for by nodal voltage equation is stated, and it is specific to position the process of more current source type oscillation sources with this For：
    Node voltage under obtained different frequency is expressed asBus admittance matrix under different frequency is expressed as Yk, then The Injection Current of each node is asked for by nodal voltage equationIt is expressed as：
    <mrow> <msub> <mover> <mi>I</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>k</mi> </msub> <mo>=</mo> <msub> <mi>Y</mi> <mi>k</mi> </msub> <msub> <mover> <mi>U</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>k</mi> </msub> </mrow>
    Based on the expression formula, there is the node that electric current injects to be positioned as more current source type oscillation sources.
CN201711309964.4A 2017-12-11 2017-12-11 A kind of localization method of more current source type oscillation sources Pending CN107968402A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711309964.4A CN107968402A (en) 2017-12-11 2017-12-11 A kind of localization method of more current source type oscillation sources

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711309964.4A CN107968402A (en) 2017-12-11 2017-12-11 A kind of localization method of more current source type oscillation sources

Publications (1)

Publication Number Publication Date
CN107968402A true CN107968402A (en) 2018-04-27

Family

ID=61999602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711309964.4A Pending CN107968402A (en) 2017-12-11 2017-12-11 A kind of localization method of more current source type oscillation sources

Country Status (1)

Country Link
CN (1) CN107968402A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114301055A (en) * 2022-02-17 2022-04-08 河海大学 Method and system for obtaining inter-harmonic power flow of power system based on broadband measurement

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969552A (en) * 2014-05-22 2014-08-06 河海大学 Harmonic source positioning and analyzing method for distributed power generation system
CN103983851A (en) * 2014-05-22 2014-08-13 河海大学 Active power distribution network harmonic source positioning analysis method
CN104898019A (en) * 2015-05-11 2015-09-09 国家电网公司 Harmonic source positioning and harmonic responsibility division method applied to active distribution network
CN105021886A (en) * 2014-12-09 2015-11-04 中国矿业大学 Harmonic source positioning device based on latent root estimation
CN105259472A (en) * 2015-10-15 2016-01-20 广州供电局有限公司 Harmonic source positioning method of electric energy quality monitoring device, and electric energy quality monitoring device thereof

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969552A (en) * 2014-05-22 2014-08-06 河海大学 Harmonic source positioning and analyzing method for distributed power generation system
CN103983851A (en) * 2014-05-22 2014-08-13 河海大学 Active power distribution network harmonic source positioning analysis method
CN105021886A (en) * 2014-12-09 2015-11-04 中国矿业大学 Harmonic source positioning device based on latent root estimation
CN104898019A (en) * 2015-05-11 2015-09-09 国家电网公司 Harmonic source positioning and harmonic responsibility division method applied to active distribution network
CN105259472A (en) * 2015-10-15 2016-01-20 广州供电局有限公司 Harmonic source positioning method of electric energy quality monitoring device, and electric energy quality monitoring device thereof

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114301055A (en) * 2022-02-17 2022-04-08 河海大学 Method and system for obtaining inter-harmonic power flow of power system based on broadband measurement
CN114301055B (en) * 2022-02-17 2024-07-12 河海大学 Method and system for acquiring harmonic power flow between power systems based on broadband measurement

Similar Documents

Publication Publication Date Title
CN108667048B (en) Frequency domain stability judging method and device for oscillation stability of new energy grid-connected system
Wang et al. Small-signal stability analysis of inverter-fed power systems using component connection method
CN106532685B (en) For the generalized impedance criterion calculation method of gird-connected inverter stability analysis and application
Suul et al. Impedance‐compensated grid synchronisation for extending the stability range of weak grids with voltage source converters
Werho et al. Power system connectivity monitoring using a graph theory network flow algorithm
CN107834830A (en) Control method and control system for uninterrupted operation of mixed MMC (modular multilevel converter)
CN104135003B (en) APF control method based on active disturbance rejection and repetitive control
Zhao et al. Subsynchronous oscillation of PV plants integrated to weak AC networks
CN106936125A (en) A kind of broad sense Second Order Integral phaselocked loop small signal impedance modeling method
Beheshtaein et al. Fault location in microgrids: a communication‐based high‐frequency impedance approach
Deng et al. Digital differential protection technique of transmission line using instantaneous active current: theory, simulation and experiment
CN107425544B (en) A kind of the electro-magnetic transient stabilising arrangement and method of distributed generation system
Xu et al. MMC admittance model simplification based on signal-flow graph
Shehata Active and reactive power control of doubly fed induction generators for wind energy generation under unbalanced grid voltage conditions
Quan et al. A concise discrete adaptive filter for frequency estimation under distorted three-phase voltage
Yuan et al. Dynamic equivalent model of VSC based on singular perturbation
CN107968402A (en) A kind of localization method of more current source type oscillation sources
CN109031029B (en) Standard voltage obtaining method and device based on single-phase earth fault
Meng et al. Conversion and SISO equivalence of impedance model of single-phase converter in electric multiple units
Wang et al. Performance improvement of a three‐phase PLL under distorted grid conditions based on frequency adaptive hybrid pre‐filtering
CN106646136A (en) Method and system for locating single-phase earth fault of parallel damping turn-adjustable type arc-extinction coil
CN107453358B (en) A kind of electric network natural resonance structure analysis method based on complex frequency domain lower node admittance matrix
Kabiri et al. A discrete state-space model for SSR stabilizing controller design for TCSC compensated systems
Zhou et al. Small-Signal Stability Boundary of Heterogeneous Multi-Converter Power Systems Dominated by the Phase-Locked Loops’ Dynamics
Huang et al. Investigation on SSCI between PMSGs‐based wind farm and AC network

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20180427