CN112636380A - Alternating current and direct current power distribution system characteristic value analysis method based on Gehr circle theory - Google Patents

Alternating current and direct current power distribution system characteristic value analysis method based on Gehr circle theory Download PDF

Info

Publication number
CN112636380A
CN112636380A CN202011534064.1A CN202011534064A CN112636380A CN 112636380 A CN112636380 A CN 112636380A CN 202011534064 A CN202011534064 A CN 202011534064A CN 112636380 A CN112636380 A CN 112636380A
Authority
CN
China
Prior art keywords
circle
gehr
gerr
ith
characteristic value
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
CN202011534064.1A
Other languages
Chinese (zh)
Other versions
CN112636380B (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.)
Institute of Electrical Engineering of CAS
Original Assignee
Institute of Electrical Engineering of CAS
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 Institute of Electrical Engineering of CAS filed Critical Institute of Electrical Engineering of CAS
Priority to CN202011534064.1A priority Critical patent/CN112636380B/en
Publication of CN112636380A publication Critical patent/CN112636380A/en
Application granted granted Critical
Publication of CN112636380B publication Critical patent/CN112636380B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • H02J3/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
    • 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
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/22—The renewable source being solar energy
    • H02J2300/24—The renewable source being solar energy of photovoltaic origin
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • 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
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/50—Photovoltaic [PV] energy
    • Y02E10/56—Power conversion systems, e.g. maximum power point trackers
    • 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
    • Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Landscapes

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

Abstract

The invention relates to a characteristic value analysis method of an alternating current and direct current power distribution system based on a Gehr circle theory, which comprises the following steps: step 1, acquiring a system matrix A of a multi-terminal AC/DC power distribution system; judging system characteristic root of any Gerr circle; when the system characteristic root of the Geoer circle exists in the left half plane or the right half plane, the whole process is exited; otherwise, the virtual shaft is covered by the Gerr circular disc and the step 2 is carried out; step 2, performing similarity transformation on the Geoer circle, and judging the stability of the Geoer circle; if the disc of the Geer circle still covers the virtual axis, performing the overlapping operation of the Geer circle; step 3, finding a Gerr circle covering the virtual axis and updating the Gerr circle; and 4, after updating the Gerr circle, obtaining a reduced system characteristic value range for analysis. The invention can rapidly analyze the range of the characteristic value when the complexity of the AC/DC distribution network system is higher, thereby judging the stability of the system and ensuring the safe and stable operation of the system.

Description

Alternating current and direct current power distribution system characteristic value analysis method based on Gehr circle theory
Technical Field
The invention relates to the field of electric power, in particular to a method for quickly analyzing characteristic values of an alternating current and direct current power distribution and utilization system based on a Gerr circle theory.
Background
In recent years, with the rapid development and wide application of new energy, new materials and new power electronic devices, the requirements of users on power supply quality, reliability and the like are increasingly increased. Meanwhile, with the gradual improvement of the permeability of distributed energy sources such as fans, photovoltaics and fuel cells, and the increase of direct current loads such as direct current lighting and electric automobile charging piles, the traditional alternating current power distribution network gradually has difficulty in meeting new requirements of direct current equipment. According to the multi-terminal alternating current and direct current hybrid power distribution technology, a converter station link and a DC/DC converter are added in a traditional alternating current power distribution network, voltage adaptation is achieved, and local consumption of distributed energy can be achieved through control strategies such as master-slave control and droop control. By means of the characteristics of flexible networking mode, various control modes and the like, the multi-terminal alternating current and direct current hybrid power distribution technology can solve the increasing problems of distributed energy access and direct current load power supply, and is a trend of future power distribution network development. By integrating a multi-terminal direct current device (comprising a plurality of converter devices and a direct current network) in an alternating current system, a multi-terminal mutual-aid interconnected power distribution and utilization system is formed. A typical structure corresponding to a multi-terminal ac/dc distribution power system is shown in fig. 1, and integrates three voltage-source converters (VSCs), wherein ac sides of the VSC1, the VSC2, and the VSC3 are respectively connected to an ac system 1, an ac system 2, and an ac system 3, and a dc side is connected to a dc bus via a line of a certain length. The direct current network can be integrated to be connected with distributed energy sources such as wind power and photovoltaic, an energy storage system, a direct current load and the like, and when the voltage level of the equipment is not matched with the voltage level of the direct current bus, the DC/DC converter can be configured to carry out conversion.
In a multi-terminal alternating current and direct current power distribution and utilization system, a characteristic value analysis method is generally adopted to analyze the small disturbance stability of the system, and then a dominant characteristic value of the system is obtained and stability analysis is carried out. Meanwhile, along with the increase of distributed power sources, direct current loads and converter stations in the alternating current and direct current hybrid system, the dimension of a system matrix is increased, the complexity of the system is correspondingly and rapidly increased, and the difficulty is brought to the traditional eigenvalue analysis method for solving the system matrix.
Disclosure of Invention
In order to solve the technical problems, the invention provides a matrix eigenvalue analysis method of an alternating current and direct current distribution system based on the Gehr circle theorem, which can be used for rapidly analyzing the range of the eigenvalue when the complexity of the alternating current and direct current distribution network system is high, so that the stability of the system is judged, and the safe and stable operation of the system is ensured.
The technical scheme of the invention is as follows: a characteristic value analysis method for an alternating current and direct current power distribution system based on a Gehr circle theory comprises the following steps:
step 1, acquiring a system matrix A of a multi-terminal AC/DC power distribution system; judging system characteristic root of any Gerr circle; when the system characteristic root of the Geoer circle exists in the left half plane or the right half plane, the whole process is exited; otherwise, the virtual shaft is covered by the Gerr circular disc and the step 2 is carried out;
step 2, performing similarity transformation on the Geoer circle, and judging the stability of the Geoer circle; if the disc of the Geer circle still covers the virtual axis, performing the overlapping operation of the Geer circle;
step 3, finding a Gerr circle covering the virtual axis and updating the Gerr circle;
and 4, after updating the Gerr circle, obtaining a reduced system characteristic value range for analysis.
Further, the step 1 specifically includes:
in all G1,G2,…,Gn,…,GNIn the cover circle, wherein GnDenotes the nth Geer circle, AnnIs the center of the n-th Gaur circle and is expressed by the position of the real axis, rnnIs the radius of the nth bell circle,
Figure BDA0002849897940000021
j is an intermediate variable, and N represents the rank of the matrix a, i.e. the total number of all the geuer circles.
Further, in the step 1, system characteristic root judgment is carried out on any Gerr circle; when the system characteristic root of the Geoer circle exists in the left half plane or the right half plane, the whole process is exited; otherwise, the virtual shaft is covered by the Gerr circular disc and the step 2 is carried out; the method specifically comprises the following steps:
if any of the Gerr circles satisfies:
Annis less than or equal to 0, and Ann+rnn<0, the system characteristic roots are all in the left half plane, and the system is stable; the whole process is exited, and the calculation is finished;
if at least one Gehr circle occurs:
Ann>0, and Ann-rnn>0, the system characteristic root is in the right half plane, and the system is unstable; the whole process is exited, and the calculation is finished;
if at least one of the Gehr circles satisfies:
Annis less than or equal to 0, and Ann+rnn>And 0, indicating that the virtual axis is covered by the Geiger circle disc, and switching to the step 2 if the judgment is not intuitive.
Further, the step 2 performs similarity transformation on the Gerr circle and judges the stability of the Gerr circle; the method specifically comprises the following steps:
step 2.1, selecting an initial Gerr circle which is covered by the circular disc and meets the preset conditions:
searching all Gerr circles, selecting N as 1,2, …, and satisfying A in NnnIs less than or equal to 0 and Ann+rnnThe Gaur circle corresponding to the largest numerical value is assumed to be the ith Gaur circle, and GiDenotes the ith Gehr circle, AiiIs the center of the ith Gaur circle and is expressed by the position of the real axis, riiIs the radius of the ith Gaur circle
Figure BDA0002849897940000031
j is an intermediate variable).
Step 2.2 changing B to diag (B)11,B22,...,Bnn,...,BNN) The ith column is updated accordingly and other data is normalized, i.e. expressed as
Figure BDA0002849897940000032
And performing a diagonal transformation, i.e.:
Figure BDA0002849897940000033
at this time, the ith circle is reduced to be within the virtual axis, other discs are correspondingly increased, at this time, the judgment is continued, if some discs of the circle still cover the virtual axis, the overlapping operation of the circles of the circle is needed, and the step 3 is switched to.
Further, the step 2.2 continues to judge, including:
if for N-1, 2, …, any of the Gerl circles in N satisfies:
Annis less than or equal to 0, and Ann+rnn<And 0, the characteristic root of the system is in the left half plane, and the system is stable. And exiting the whole flow and finishing the calculation.
If at least one Gehr circle occurs:
Ann>0, and Ann-rnn>0, the system characteristic root is in the right half-plane, and the system is unstable. And exiting the whole flow and finishing the calculation.
If at least one of the Gehr circles satisfies:
Annis less than or equal to 0, and Ann+rnn>0, indicating that the disc with the Geer circle still covers the virtual axis, and performing Geer circle overlapping operation to reduce the cross area among different discs to the maximum and reduce further characteristic root range search; and (5) turning to the step 3.
Further, the step 3 finds a bell circle covering the virtual axis and updates the bell circle, which specifically includes:
step 3.1 searches all the gehr circles, selects N as 1,2, …, N satisfies annIs less than or equal to 0 and Ann+rnnThe Gaur circle with the largest value, assumed to be the k-th Gaur circle, is designated by GkDenotes the k-th Gehr circle, AkkThe center of the k-th Gaier circle is located at the position of the real axisLine represents, rkkIs the radius of the k-th bell circle,
Figure BDA0002849897940000034
j is an intermediate variable;
step 3.2 changing B to diag (B)11,B22,...,Bkk,...,BNN) The kth column is updated accordingly and the other data is normalized, i.e. denoted as B ═ diag (B)11＝1,B22＝1,...,Bkk＝x,...,BNN1) and performs a diagonal transformation, i.e.:
Figure BDA0002849897940000041
Figure BDA0002849897940000042
wherein the content of the first and second substances,
C＝Aii 2Aik 2-4Aii 2Aikrii+4Aii 2rii 2+2AiiAik 2rii+2AiiAikAkkrii-6AiiAikrii 2+4rkkAiiAikrii-4AiiAkkrii 2+4Aiirii 3+Aik 2rii 2+2AikAkkrii 2-2Aikrii 3+4rkkAikrii 2+Akk 2rii 2-2Akkrii 3+rii 4
Aiirepresenting the ith row and ith column element in the system matrix, AkkRepresenting the kth row and kth column elements, A, in the system matrixikRow ith and column kth elements, riiRepresents the ith coverRadius of the circle, rkkRepresenting the k-th radius of the bell circle.
And 4, after updating the Gerr circle, obtaining a reduced system characteristic value range for analysis.
And ending and exiting the process.
Has the advantages that:
the alternating current and direct current hybrid power distribution technology provides an effective technical means for the large-scale access and the optimized operation control of distributed energy sources, and becomes one of the important forms of the future power distribution network. When the system safety and stability is analyzed, a characteristic value analysis method is generally adopted to analyze the small disturbance stability of the system, so as to obtain a leading characteristic value and a transmission power boundary of the system, and further determine the system safety margin. Therefore, the invention provides the quick analysis method for the characteristic value of the alternating current and direct current power distribution system based on the Gehr circle theorem, solves the problem that the accurate value of the characteristic value is difficult to obtain when the complexity of the alternating current and direct current power distribution system is high, fills up the blank of the related technology, and has wide application prospect.
Drawings
FIG. 1 is a typical structure diagram of a multi-terminal AC/DC distribution system;
fig. 2 is a flow chart of a method for analyzing a characteristic value of an ac/dc distribution system based on the gehr circle theory.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, rather than all embodiments, and all other embodiments obtained by a person skilled in the art based on the embodiments of the present invention belong to the protection scope of the present invention without creative efforts.
According to the embodiment of the invention, the alternating current and direct current power distribution system characteristic value analysis method based on the Gehr circle theory specifically comprises the following steps:
step 1, acquiring a system matrix A of a multi-terminal AC/DC power distribution system; judging system characteristic root of any Gerr circle; when the system characteristic root of the Geoer circle exists in the left half plane or the right half plane, the whole process is exited; otherwise, the virtual shaft is covered by the Gerr circular disc and the step 2 is carried out;
step 2, performing similarity transformation on the Geoer circle, and judging the stability of the Geoer circle; if the disc of the Geer circle still covers the virtual axis, performing the overlapping operation of the Geer circle;
step 3, finding a Gerr circle covering the virtual axis and updating the Gerr circle;
and 4, after updating the Gerr circle, obtaining a reduced system characteristic value range for analysis.
The principle of the invention is as follows: after acquiring the system matrix a of the multi-terminal dc system, assume that:
Figure BDA0002849897940000051
Aijan element representing the ith row and the jth column in the system matrix A;
considering that the order of the system matrix is higher and higher, the calculation amount is larger and larger when each detailed characteristic root is solved, so that the approximate position of the characteristic root is quickly judged by means of the Gerr circle theory to accelerate the analysis speed.
With GnRepresents the nth Gehr circle of the matrix A, AnnIs the center of the n-th Gehr circle, rnnIs the radius of the nth Gehr circle (
Figure BDA0002849897940000052
j is an intermediate variable), N represents the rank of the matrix A, and according to the Gehr circle theory, the characteristic roots of the matrix A are all in the union of N Gehr circles.
Then each Gerr circle on the complex plane is represented as:
G1：{|z-A11|≤r11}
G2：{|z-A22|≤r22}
...
Gnn：{|z-Ann|≤rnn}
...
GNN：{|z-ANN|≤rNN}
and further finishing to obtain:
Figure BDA0002849897940000061
constructing a diagonal matrix B, and setting:
Figure BDA0002849897940000062
at this time, diagonal transformation is performed:
Figure BDA0002849897940000063
after the similarity transformation of the diagonal matrix is carried out, the characteristic root of the matrix can be not changed, the radius of each circle can be reduced or enlarged according to the principle, several originally coincident circles are effectively separated, or the circle near the virtual axis is further reduced, so that the characteristic value is rapidly judged on the leftmost semi-plane or the right semi-plane, and the stability analysis is facilitated.
According to the specific embodiment of the invention, as shown in fig. 2, a method for rapidly analyzing a characteristic value of an ac/dc power distribution system based on a bell circle theory includes the following specific implementation steps:
step 1, acquiring a system matrix A of a multi-terminal AC/DC power distribution system; judging system characteristic root of any Gerr circle; when the system characteristic root of the Geoer circle exists in the left half plane or the right half plane, the whole process is exited; otherwise, the virtual shaft is covered by the Gerr circular disc and the step 2 is carried out;
all G1,G2,…,Gn,…,GNIn the cover circle, wherein GnDenotes the nth Geer circle, AnnIs the center of the n-th Gaur circle and is expressed by the position of the real axis, rnnIs the radius of the nth Gehr circle (
Figure BDA0002849897940000064
j is an intermediate variable), N represents the rank of matrix a, i.e. the total number of all geuer circles. If any of the Gerr circles satisfies:
Annis less than or equal to 0, and Ann+rnn<And 0, the characteristic root of the system is in the left half plane, and the system is stable. And exiting the whole flow and finishing the calculation.
If at least one Gehr circle occurs:
Ann>0, and Ann-rnn>0, the system characteristic root is in the right half-plane, and the system is unstable. And exiting the whole flow and finishing the calculation.
If at least one of the Gehr circles satisfies:
Annis less than or equal to 0, and Ann+rnn>And 0, indicating that the virtual axis is covered by the Geiger circle disc, and turning to the step 2 if the virtual axis is not visually judged.
Step 2: performing similarity transformation on the Geer circle, and judging the stability of the Geer circle;
step 2.1, selecting an initial Gerr circle which is covered by the circular disc and meets the preset conditions:
searching all Gerr circles, selecting N as 1,2, …, and satisfying A in NnnIs less than or equal to 0 and Ann+rnnThe Gaur circle corresponding to the largest numerical value is assumed to be the ith Gaur circle, and GiDenotes the ith Gehr circle, AiiIs the center of the ith Gaur circle and is expressed by the position of the real axis, riiIs the radius of the ith Gaur circle
Figure BDA0002849897940000071
j is an intermediate variable).
Step 2.2 changing B to diag (B)11,B22,...,Bnn,...,BNN) The ith column is updated accordingly and other data is normalized, i.e. expressed as
Figure BDA0002849897940000072
And performing a diagonal transformation, i.e.:
Figure BDA0002849897940000073
at this moment, the ith Gehr circle is reduced to be within the virtual axis, other discs are correspondingly increased, and the judgment is continued:
if for N-1, 2, …, any of the Gerl circles in N satisfies:
Annis less than or equal to 0, and Ann+rnn<And 0, the characteristic root of the system is in the left half plane, and the system is stable. And exiting the whole flow and finishing the calculation.
If at least one Gehr circle occurs:
Ann>0, and Ann-rnn>0, the system characteristic root is in the right half-plane, and the system is unstable. And exiting the whole flow and finishing the calculation.
If at least one of the Gehr circles satisfies:
Annis less than or equal to 0, and Ann+rnn>0, indicating that the disc with the Geer circle still covers the virtual axis, and performing Geer circle overlapping operation to reduce the cross area among different discs to the maximum and reduce further characteristic root range search; turning to the step 3;
step 3, finding the Gerr circle covering the virtual axis and updating the Gerr circle;
step 3.1 searches all the gehr circles, selects N as 1,2, …, N satisfies annIs less than or equal to 0 and Ann+rnnThe Gaur circle with the largest value, assumed to be the k-th Gaur circle, is designated by GkDenotes the k-th Gehr circle, AkkIs the center of the kth Gaier circle, and is expressed by the position of the real axis, rkkIs the radius of the kth Gehr circle (
Figure BDA0002849897940000081
j is an intermediate variable).
Step 3.2 changing B to diag (B)11,B22,...,Bkk,...,BNN) The kth column is updated accordingly and the other data is normalized, i.e. denoted as B ═ diag (B)11＝1,B22＝1,...,Bkk＝x,...,BNN1) and performs a diagonal transformation, i.e.:
Figure BDA0002849897940000082
Figure BDA0002849897940000083
wherein the content of the first and second substances,
C＝Aii 2Aik 2-4Aii 2Aikrii+4Aii 2rii 2+2AiiAik 2rii+2AiiAikAkkrii-6AiiAikrii 2+4rkkAiiAikrii-4AiiAkkrii 2+4Aiirii 3+Aik 2rii 2+2AikAkkrii 2-2Aikrii 3+4rkkAikrii 2+Akk 2rii 2-2Akkrii 3+rii 4
Aiirepresenting the ith row and ith column element in the system matrix, AkkRepresenting the kth row and kth column elements, A, in the system matrixikRow ith and column kth elements, riiRepresents the radius of the ith Gehr circle, rkkRepresenting the k-th radius of the bell circle.
And 4, after updating the Gerr circle, obtaining a reduced system characteristic value range for analysis.
And ending and exiting the process.
Although illustrative embodiments of the present invention have been described above to facilitate the understanding of the present invention by those skilled in the art, it should be understood that the present invention is not limited to the scope of the embodiments, but various changes may be apparent to those skilled in the art, and it is intended that all inventive concepts utilizing the inventive concepts set forth herein be protected without departing from the spirit and scope of the present invention as defined and limited by the appended claims.

Claims (6)

1. A method for analyzing a characteristic value of an alternating current and direct current power distribution system based on a Gehr circle theory is characterized by comprising the following steps:
step 1, acquiring a system matrix A of a multi-terminal AC/DC power distribution system; judging system characteristic root of any Gerr circle; when the system characteristic root of the Geoer circle exists in the left half plane or the right half plane, the whole process is exited; otherwise, the virtual shaft is covered by the Gerr circular disc, and the step 2 is carried out;
step 2, performing similarity transformation on the Geoer circle, and judging the stability of the Geoer circle; if the disc of the Geer circle still covers the virtual axis, performing the overlapping operation of the Geer circle;
step 3, finding a Gerr circle covering the virtual axis and updating the Gerr circle;
and 4, after updating the Gerr circle, obtaining a reduced system characteristic value range for analysis.
2. The method for analyzing the characteristic value of the alternating current and direct current distribution power system based on the Gehr circle theory as claimed in claim 1, wherein the step 1 specifically comprises:
in all G1,G2,…,Gn,…,GNIn the cover circle, wherein GnDenotes the nth Geer circle, AnnIs the center of the n-th Gaur circle and is expressed by the position of the real axis, rnnIs the radius of the nth bell circle,
Figure FDA0002849897930000011
j is an intermediate variable, N representsThe rank of the matrix a is the total number of all the geuer circles.
3. The method for analyzing the characteristic value of the alternating current and direct current power distribution system based on the Gehr circle theory as claimed in claim 1, wherein the step 1 is to perform system characteristic root judgment on any Gehr circle; when the system characteristic root of the Geoer circle exists in the left half plane or the right half plane, the whole process is exited; otherwise, the virtual shaft is covered by the Gerr circular disc, and the step 2 is carried out; the method specifically comprises the following steps:
if any of the Gerr circles satisfies:
Annis less than or equal to 0, and Ann+rnn<0, the system characteristic roots are all in the left half plane, and the system is stable; the whole process is exited, and the calculation is finished;
if at least one Gehr circle occurs:
Ann>0, and Ann-rnn>0, the system characteristic root is in the right half plane, and the system is unstable; the whole process is exited, and the calculation is finished;
if at least one of the Gehr circles satisfies:
Annis less than or equal to 0, and Ann+rnn>And 0, indicating that the virtual axis is covered by the Geiger circle disc, and turning to the step 2 if the virtual axis is not visually judged.
4. The method for analyzing the characteristic value of the alternating current and direct current distribution system based on the Gehr circle theory as claimed in claim 1, wherein the step 2 is to perform similarity transformation on the Gehr circle and judge the stability of the Gehr circle; the method specifically comprises the following steps:
step 2.1, selecting an initial Gerr circle which is covered by the circular disc and meets the preset conditions:
searching all Gerr circles, selecting N as 1,2, …, and satisfying A in NnnIs less than or equal to 0 and Ann+rnnThe Gaur circle corresponding to the largest numerical value is assumed to be the ith Gaur circle, and GiDenotes the ith Gehr circle, AiiIs the center of the ith Gaur circle and is expressed by the position of the real axis, riiIs the radius of the ith bell circle,
Figure FDA0002849897930000021
j is an intermediate variable;
step 2.2 changing B to diag (B)11,B22,...,Bnn,...,BNN) The ith column is updated accordingly and other data is normalized, i.e. expressed as
Figure FDA0002849897930000022
And performing a diagonal transformation, i.e.:
Figure FDA0002849897930000023
at this time, the ith circle is reduced to be within the virtual axis, other discs are correspondingly increased, at this time, the judgment is continued, if some discs of the circle still cover the virtual axis, the overlapping operation of the circles of the circle is needed, and the step 3 is switched to.
5. The method for analyzing the characteristic value of the ac/dc distribution system based on the bell circle theory as claimed in claim 2, wherein the step 2.2 of continuing the judgment comprises:
if for N ═ 1,2, …, in N, the Gerr circle satisfies:
Annis less than or equal to 0, and Ann+rnn<0, the system characteristic roots are all in the left half plane, and the system is stable; the whole process is exited, and the calculation is finished;
if at least one Gehr circle occurs:
Ann>0, and Ann-rnn>0, if the system characteristic root is in the right half plane, the system is unstable, the whole process is exited, and the calculation is finished;
if at least one of the Gehr circles satisfies:
Annis less than or equal to 0, and Ann+rnn>0, it indicates that the disc with the cover circle still covers the virtual axis, and the cover circle overlapping operation is needed to minimize the defectThe cross area between the same disc reduces further characteristic root range search; and (5) turning to the step 3.
6. The method for analyzing the characteristic value of the ac/dc distribution system based on the bell circle theory as claimed in claim 1, wherein the step 3 of finding the bell circle covering the virtual axis and updating the bell circle comprises:
step 3.1 searches all the gehr circles, selects N as 1,2, …, N satisfies annIs less than or equal to 0 and Ann+rnnThe Gaur circle with the largest value, assumed to be the k-th Gaur circle, is designated by GkDenotes the k-th Gehr circle, AkkIs the center of the kth Gaier circle, and is expressed by the position of the real axis, rkkIs the radius of the k-th bell circle,
Figure FDA0002849897930000031
j is an intermediate variable;
step 3.2 changing B to diag (B)11,B22,...,Bkk,...,BNN) The kth column is updated accordingly and the other data is normalized, i.e. denoted as B ═ diag (B)11＝1,B22＝1,...,Bkk＝x,...,BNN1) and performs a diagonal transformation, i.e.:
Figure FDA0002849897930000032
Figure FDA0002849897930000033
wherein the content of the first and second substances,
Figure FDA0002849897930000034
Aiirepresenting the ith row and ith column element in the system matrix, AkkRepresenting the kth row and kth column elements, A, in the system matrixikRow ith and column kth elements, riiRepresents the radius of the ith Gehr circle, rkkRepresenting the k-th radius of the bell circle.
CN202011534064.1A 2020-12-21 2020-12-21 Alternating current and direct current power distribution system characteristic value analysis method based on Gehr circle theory Active CN112636380B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011534064.1A CN112636380B (en) 2020-12-21 2020-12-21 Alternating current and direct current power distribution system characteristic value analysis method based on Gehr circle theory

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011534064.1A CN112636380B (en) 2020-12-21 2020-12-21 Alternating current and direct current power distribution system characteristic value analysis method based on Gehr circle theory

Publications (2)

Publication Number Publication Date
CN112636380A true CN112636380A (en) 2021-04-09
CN112636380B CN112636380B (en) 2022-11-18

Family

ID=75321214

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011534064.1A Active CN112636380B (en) 2020-12-21 2020-12-21 Alternating current and direct current power distribution system characteristic value analysis method based on Gehr circle theory

Country Status (1)

Country Link
CN (1) CN112636380B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117405989A (en) * 2023-09-25 2024-01-16 国家电网有限公司华东分部 Detection method and device for broadband oscillation phenomenon, storage medium and computer equipment

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106125715A (en) * 2016-06-21 2016-11-16 西南交通大学 A kind of MIMO cascade system method for analyzing stability based on impedance return ratio matrix
CN107565545A (en) * 2017-07-28 2018-01-09 西南交通大学 A kind of tractor net coupled system stability criteria computational methods based on impedance return ratio matrix

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106125715A (en) * 2016-06-21 2016-11-16 西南交通大学 A kind of MIMO cascade system method for analyzing stability based on impedance return ratio matrix
CN107565545A (en) * 2017-07-28 2018-01-09 西南交通大学 A kind of tractor net coupled system stability criteria computational methods based on impedance return ratio matrix

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
雍龙泉: "基于盖尔圆定理的矩阵特征值估计", 《陕西理工大学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117405989A (en) * 2023-09-25 2024-01-16 国家电网有限公司华东分部 Detection method and device for broadband oscillation phenomenon, storage medium and computer equipment
CN117405989B (en) * 2023-09-25 2024-03-22 国家电网有限公司华东分部 Detection method and device for broadband oscillation phenomenon, storage medium and computer equipment

Also Published As

Publication number Publication date
CN112636380B (en) 2022-11-18

Similar Documents

Publication Publication Date Title
CN111539566B (en) Power distribution network multi-fault first-aid repair recovery method and system considering pre-scheduling before disaster
Cai et al. Distributed control scheme for package-level state-of-charge balancing of grid-connected battery energy storage system
CN112993979A (en) Power distribution network reactive power optimization method and device, electronic equipment and storage medium
CN114094600B (en) Collaborative operation control method and system for multi-optical storage VSG system
CN112636380B (en) Alternating current and direct current power distribution system characteristic value analysis method based on Gehr circle theory
CN112803434A (en) Reactive power optimization method, device, equipment and storage medium for active power distribution network
CN110808600A (en) Method and system for calculating access of battery energy storage station in transformer substation
CN114919752A (en) ECMS-MPC-based energy management method for hydrogen fuel hybrid unmanned aerial vehicle
CN114094599A (en) Multi-station fusion participating peak-shaving frequency-modulation voltage-regulation potential evaluation method and device
CN113794193A (en) Decision-making method for hydrogen production of renewable energy direct-current micro-grid
CN113162053A (en) Power distribution network voltage control method and system based on load transfer and energy storage regulation
CN112467717A (en) Hybrid energy system real-time load distribution method based on fuzzy control
CN111478335A (en) Power distribution network load flow calculation method and system considering distributed photovoltaic
CN113824161B (en) Control method for AC/DC micro-grid cluster system under DoS attack
CN114336756B (en) Camera adjustment configuration method and system of new energy island Direct Current (DC) outgoing system
CN114180023B (en) Multi-energy ship control management method and device based on load prediction algorithm
CN107104429B (en) power system load scheduling method considering distributed energy storage system
CN110729759B (en) Method and device for determining distributed power supply configuration scheme in micro-grid
CN111626591B (en) Risk assessment method for flexible platform area
CN110504725B (en) Multi-cell-stack rapid equalization control method and device for energy storage power station
Nguyen et al. Power loss minimization by optimal placement of distributed generation considering the distribution network configuration based on artificial ecosystem optimization
CN112736948A (en) Power adjusting method and device for energy storage system in charging station
Li et al. Development of control strategy to increase the lifetime of grid-connected Li-ion battery
CN111311032A (en) Micro-grid system capacity optimal configuration method based on sector radar map model
Jaya et al. Mode Transition in DC Microgrids with Non-Dispatchable Sources

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