CN108879754B - FCL optimal configuration method and system for reducing risk of secondary commutation failure - Google Patents

FCL optimal configuration method and system for reducing risk of secondary commutation failure Download PDF

Info

Publication number
CN108879754B
CN108879754B CN201810835675.6A CN201810835675A CN108879754B CN 108879754 B CN108879754 B CN 108879754B CN 201810835675 A CN201810835675 A CN 201810835675A CN 108879754 B CN108879754 B CN 108879754B
Authority
CN
China
Prior art keywords
fault current
node
determining
current limiter
branch
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.)
Active
Application number
CN201810835675.6A
Other languages
Chinese (zh)
Other versions
CN108879754A (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.)
Tsinghua University
State Grid Shanghai Electric Power Co Ltd
Original Assignee
Tsinghua University
State Grid Shanghai Electric Power Co Ltd
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 Tsinghua University, State Grid Shanghai Electric Power Co Ltd filed Critical Tsinghua University
Priority to CN201810835675.6A priority Critical patent/CN108879754B/en
Publication of CN108879754A publication Critical patent/CN108879754A/en
Application granted granted Critical
Publication of CN108879754B publication Critical patent/CN108879754B/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
    • 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]
    • 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 discloses a fault current limiter optimal configuration method and system for reducing the risk of secondary commutation failure. The method comprises the following steps: selecting voltage interaction factors among a plurality of adjacent direct current branches, and determining a target function by taking the minimum voltage interaction factor as a target; determining the installation position and the current-limiting reactance value of the fault current limiter by taking one fault current limiter installed at each time and the maximum number of fault current limiters installed at set times as constraint conditions; after the installation position and the current-limiting reactance value of the fault current limiter are determined, calculating voltage interaction factors between adjacent direct current branches of the multi-feed-in high-voltage direct-current power transmission system, jumping to the step of selecting the voltage interaction factors between the multiple adjacent direct current branches, determining a target function step by taking the minimum voltage interaction factor as a target, and stopping iteration until the fault current limiter reaches the set installation number. The method and the system provided by the invention can meet the FCL optimal configuration scheme of different commutation failure risk constraint conditions.

Description

FCL optimal configuration method and system for reducing risk of secondary commutation failure
Technical Field
The invention relates to the field of high-voltage power transmission, in particular to an FCL optimal configuration method and system for reducing the risk of secondary commutation failure.
Background
In a Multi-feed High-Voltage Direct Current (MI-HVDC) system, secondary commutation failure refers to a related record of the prevention of commutation failure phenomenon caused by commutation Voltage drop impact caused by that an inverter station in the MI-HVDC system suffers from alternating Current system failure and adjacent inverter stations suffer from commutation failure, and commutation Failure (FCL) occurs.
Disclosure of Invention
The invention aims to provide an FCL optimal configuration method and system for reducing the risk of secondary commutation failure, and an FCL optimal configuration scheme capable of meeting different commutation failure risk constraint conditions.
In order to achieve the purpose, the invention provides the following scheme:
a fault current limiter optimized configuration method to reduce the risk of a secondary commutation failure, the method comprising:
selecting voltage interaction factors among a plurality of adjacent direct current branches, and determining a target function by taking the minimum voltage interaction factor as a target;
determining the installation position and the current-limiting reactance value of the fault current limiter by taking one fault current limiter installed at each time and the maximum number of fault current limiters installed at set times as constraint conditions;
after the installation position and the current-limiting reactance value of the fault current limiter are determined, calculating voltage interaction factors between adjacent direct current branches of the multi-feed-in high-voltage direct-current power transmission system, jumping to the step of selecting the voltage interaction factors between the multiple adjacent direct current branches, determining a target function step by taking the minimum voltage interaction factor as a target, and stopping iteration until the fault current limiter reaches the set installation number.
Optionally, the determining the objective function with the minimum voltage interaction factor as a target specifically includes:
determining a multi-objective optimization function
Figure BDA0001744435440000021
Where i denotes the number of the optimization objective,
Figure BDA0001744435440000022
and
Figure BDA0001744435440000023
the numbers of two adjacent dc branches are indicated,
Figure BDA0001744435440000024
representing branches
Figure BDA0001744435440000025
And branch
Figure BDA0001744435440000026
Voltage interaction factor between;
determining an objective function from the multi-objective optimization function
Figure BDA0001744435440000027
Wherein, aiTo represent
Figure BDA0001744435440000028
The weight of (a) is determined,
Figure BDA0001744435440000029
indicates the branch after the t-1 iteration
Figure BDA00017444354400000210
And branch
Figure BDA00017444354400000211
Voltage interaction factor of.
Optionally, the determining the installation position and the current-limiting reactance value of the fault current limiter by using one fault current limiter installed at a time and the maximum number of fault current limiters installed as constraint conditions specifically includes:
determining the initial iteration number t as 1 sumInitial number n of barrier current limitersFCL＝1；
Determining a first constraint
Figure BDA00017444354400000212
Wherein N isLAs the number of branches, a branch information matrix is
Figure BDA00017444354400000213
Is a variable of 0, 1, lcA number indicating the preselected installation leg of the fault current limiter,
Figure BDA00017444354400000214
indicating fault current limiter installed in first branch setcOn the branch, only the installation position of one fault current limiter is determined in each iteration, only one fault current limiter is installed,
Figure BDA00017444354400000215
is the l-th of CLcLine, represents the lcThe information of the branch is then transmitted to the subscriber,
Figure BDA0001744435440000031
Figure BDA0001744435440000032
respectively represent the lcA head end node, a tail end node and a branch impedance of the branch,
Figure BDA0001744435440000033
Figure BDA0001744435440000034
node impedance matrixes Z respectively corresponding to t-th iterationtIs located at
Figure BDA0001744435440000035
Respectively represent nodes
Figure BDA0001744435440000036
Self-impedance, node
Figure BDA0001744435440000037
Self-impedance, node
Figure BDA0001744435440000038
And node
Figure BDA0001744435440000039
The mutual impedance between the two electrodes is high,
Figure BDA00017444354400000310
are each ZtIs located at
Figure BDA00017444354400000311
Figure BDA00017444354400000312
Respectively represent nodes
Figure BDA00017444354400000313
And node
Figure BDA00017444354400000314
Mutual impedance, node
Figure BDA00017444354400000315
And node
Figure BDA00017444354400000316
Mutual impedance, node
Figure BDA00017444354400000317
And node
Figure BDA00017444354400000318
Mutual impedance, node
Figure BDA00017444354400000319
And node
Figure BDA00017444354400000320
The mutual impedance of (a);
determining a second constraint
Figure BDA00017444354400000321
Wherein the content of the first and second substances,
Figure BDA00017444354400000322
determining a third constraint
Figure BDA00017444354400000323
Wherein the content of the first and second substances,
Figure BDA00017444354400000324
c＝z1+z2-2z3+zFCL，zFCL＝-(z4 2+z4jx)/(jx), x being the reactance value of the fault current limiter, j being the imaginary sign;
determining that the fourth constraint x is less than or equal to xmaxWhere x is the number of fault current limiters installed, xmaxA set number of fault current limiters is installed at the maximum.
The invention also provides a fault current limiter optimal configuration system for reducing the risk of secondary commutation failure, which comprises:
the target function determining module is used for selecting voltage interaction factors among a plurality of adjacent direct current branches, and determining a target function by taking the minimum voltage interaction factor as a target;
the fault current limiter determining module is used for determining the installation position and the current limiting reactance value of the fault current limiter by taking the condition that one fault current limiter is installed at a time and the set number of fault current limiters are installed at most as constraint conditions;
and the voltage interaction factor calculation module is used for calculating voltage interaction factors between adjacent direct current branches of the multi-feed-in high-voltage direct-current power transmission system after the installation position and the current-limiting reactance value of the fault current limiter are determined, jumping to the step of selecting the voltage interaction factors between the multiple adjacent direct current branches, and determining a target function step by taking the minimum voltage interaction factor as a target until the fault current limiter reaches the set installation number.
Optionally, the objective function determining module specifically includes:
a multi-objective optimization function determination unit for determining a multi-objective optimization function
Figure BDA0001744435440000041
Where i denotes the number of the optimization objective,
Figure BDA00017444354400000415
and
Figure BDA0001744435440000044
the numbers of two adjacent dc branches are indicated,
Figure BDA0001744435440000045
representing branches
Figure BDA0001744435440000046
And branch
Figure BDA0001744435440000047
Voltage interaction factor between;
an objective function determination unit for determining an objective function based on the multi-objective optimization function
Figure BDA0001744435440000048
Wherein, aiTo represent
Figure BDA0001744435440000049
The weight of (a) is determined,
Figure BDA00017444354400000410
indicates the branch after the t-1 iteration
Figure BDA00017444354400000411
And branch
Figure BDA00017444354400000412
Voltage interaction factor of.
Optionally, the voltage interaction factor calculation module specifically includes:
an initial condition determining unit for determining an initial iteration number t equal to 1 and an initial number n of fault current limitersFCL＝1；
A first constraint condition determination unit for determining a first constraint condition
Figure BDA00017444354400000413
Wherein N isLAs the number of branches, a branch information matrix is
Figure BDA00017444354400000414
Is a variable of 0, 1, lcA number indicating the preselected installation leg of the fault current limiter,
Figure BDA0001744435440000051
indicating fault current limiter installed in first branch setcOn the branch, only the installation position of one fault current limiter is determined in each iteration, only one fault current limiter is installed,
Figure BDA0001744435440000052
is the l-th of CLcLine, represents the lcThe information of the branch is then transmitted to the subscriber,
Figure BDA0001744435440000053
Figure BDA0001744435440000054
respectively represent the lcA head end node, a tail end node and a branch impedance of the branch,
Figure BDA0001744435440000055
Figure BDA0001744435440000056
node impedance matrixes Z respectively corresponding to t-th iterationtIs located at
Figure BDA0001744435440000057
Respectively represent nodes
Figure BDA0001744435440000058
Self-impedance, node
Figure BDA0001744435440000059
Self-impedance, node
Figure BDA00017444354400000510
And node
Figure BDA00017444354400000511
The mutual impedance between the two electrodes is high,
Figure BDA00017444354400000512
are each ZtIs located at
Figure BDA00017444354400000513
Figure BDA00017444354400000514
Respectively represent nodes
Figure BDA00017444354400000515
And node
Figure BDA00017444354400000516
Mutual impedance, node
Figure BDA00017444354400000517
And node
Figure BDA00017444354400000518
Mutual impedance, node
Figure BDA00017444354400000519
And node
Figure BDA00017444354400000520
Mutual impedance, node
Figure BDA00017444354400000521
And node
Figure BDA00017444354400000522
The mutual impedance of (a);
a second constraint condition determination unit for determining a second constraint condition
Figure BDA00017444354400000523
Wherein the content of the first and second substances,
Figure BDA00017444354400000524
a third constraint condition determination unit for determining a third constraint condition
Figure BDA00017444354400000525
Wherein the content of the first and second substances,
Figure BDA00017444354400000526
c＝z1+z2-2z3+zFCL，zFCL＝-(z4 2+z4jx)/(jx), x being the reactance value of the fault current limiter, j being the imaginary sign;
a fourth constraint condition determination unit for determining that the fourth constraint condition x is less than or equal to xmaxWhere x is the number of fault current limiters installed, xmaxA set number of fault current limiters is installed at the maximum.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects: according to the FCL optimal configuration method and the FCL optimal configuration system for reducing the risk of the secondary commutation failure, provided by the invention, the optimal configuration scheme of the FCL is solved by taking the minimum MIIF index as a target and the installation cost of the FCL as a constraint, and the risk of the secondary commutation failure is reduced by reducing the MIIF index between adjacent inversion stations after the failure.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
Fig. 1 is a flowchart of an FCL optimization configuration method for reducing the risk of a secondary commutation failure according to an embodiment of the present invention;
fig. 2 is a structural diagram of an FCL optimized configuration system for reducing the risk of a secondary commutation failure according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the 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, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide an FCL optimal configuration method and system for reducing the risk of secondary commutation failure, and an FCL optimal configuration scheme capable of meeting different commutation failure risk constraint conditions.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Fig. 1 is a flowchart of an FCL optimization configuration method for reducing the risk of the secondary commutation failure according to an embodiment of the present invention, and as shown in fig. 1, the risk of the secondary commutation failure occurring in the MI-HVDC system is positively correlated with the MIIF index, and the risk of the secondary commutation failure can be reduced by reducing the MIIF index between adjacent inverter stations after the failure. The method takes the minimum MIIF index as a target and the installation cost of the FCL as a constraint to solve the optimal configuration scheme of the FCL. The method comprises the following specific steps:
step 1: the number of initialization iterations t is 1, and the number n of FCLs is setFCL＝1；
Step 2: and solving an objective function of the t iteration. Calculating the proportion of the risk of the secondary commutation failure after the t-1 iteration is finished, and selecting according to the calculation condition
Figure BDA0001744435440000061
The MIIF index is the target of optimization. If it is
Figure BDA0001744435440000062
Optimization is not required; if it is
Figure BDA0001744435440000071
The optimization problem is a multi-objective optimization problem. The objective function is shown in equation (1).
Figure BDA0001744435440000072
In the formula, i represents the number of the optimization target,
Figure BDA0001744435440000073
and
Figure BDA0001744435440000074
representing multi-feed interaction factors
Figure BDA0001744435440000075
Number of associated two HVDC rounds.
And step 3: evaluating the t-th iteration according to the result of the t-1 th iteration
Figure BDA0001744435440000076
Branch information matrix CL of order. CL contains NLThe information of the branch is expressed as follows:
Figure BDA0001744435440000077
in the formula (I), the compound is shown in the specification,
Figure BDA0001744435440000078
is one
Figure BDA0001744435440000079
Row vector of orderIt is the l th in CLcLine, represents the lcInformation of the row branch. The information is whencWhen a row leg is selected as the FCL's installation location, it is necessary to compute multiple objective functions.
Figure BDA00017444354400000710
The expression of (a) is as follows:
Figure BDA00017444354400000711
in the formula (I), the compound is shown in the specification,
Figure BDA00017444354400000712
respectively represent the lcThe line branch's head-end node and branch impedance,
Figure BDA00017444354400000713
Figure BDA00017444354400000714
node impedance matrixes Z respectively corresponding to t-th iterationtIs located at
Figure BDA00017444354400000715
Figure BDA00017444354400000716
Respectively represent nodes
Figure BDA00017444354400000717
Self-impedance, node
Figure BDA00017444354400000718
Self-impedance, node
Figure BDA00017444354400000719
And node
Figure BDA00017444354400000720
The mutual impedance between the two electrodes is high,
Figure BDA0001744435440000081
are each ZtIs located at
Figure BDA0001744435440000082
Figure BDA0001744435440000083
Respectively represent nodes
Figure BDA0001744435440000084
And node
Figure BDA0001744435440000085
Mutual impedance, node
Figure BDA0001744435440000086
And node
Figure BDA0001744435440000087
Mutual impedance, node
Figure BDA0001744435440000088
And node
Figure BDA0001744435440000089
Mutual impedance, node
Figure BDA00017444354400000810
And node
Figure BDA00017444354400000811
The mutual impedance of (a).
And 4, step 4: and solving an optimization problem. Converting the multi-objective optimization problem into a single-objective optimization problem by a weighting coefficient method, wherein an objective function is as follows:
Figure BDA00017444354400000812
in the formula, aiTo represent
Figure BDA00017444354400000813
The weight of (a) is determined,
Figure BDA00017444354400000814
indicating after FCL startup
Figure BDA00017444354400000815
To pair
Figure BDA00017444354400000816
Voltage interaction factor of.
Constraint one:
Figure BDA00017444354400000817
in the formula (I), the compound is shown in the specification,
Figure BDA00017444354400000818
is a variable of 0, 1, lcA number representing a pre-selected branch,
Figure BDA00017444354400000819
indicating installation of FCL in preselected branch setcOn a bar, since only one FCL mounting location is determined per iteration, equation (5) should be satisfied.
Constraint two:
Figure BDA00017444354400000820
the left side of the formula (6) is one
Figure BDA00017444354400000821
The row vector of the order is obtained by the expression on the right side of the formula (6), wherein sigma is NLThe dimension-row vector is then calculated,
Figure BDA00017444354400000822
constraint condition three:
Figure BDA00017444354400000823
wherein the content of the first and second substances,
Figure BDA00017444354400000824
Figure BDA0001744435440000091
c＝z1+z2-2z3+zFCL(10)
zFCL＝-(z4 2+z4·jx)/(jx) (11)
where x is the reactance value of FCL and j is the imaginary symbol.
Constraint condition four:
x≤xmax(12)
in the formula, xmaxIs the upper limit value of x.
Solving the nonlinear mixed integer programming problem by using Bonmin in the OPTI-Toolbox to obtain the optimal installation position and the current-limiting reactance value of the FCL in the iterative process of the t step.
Solving the nonlinear mixed integer programming problem by using Bonmin in the OPTI-Toolbox to obtain the optimal installation position and the current-limiting reactance value of the FCL in the iterative process of the t step.
And 5: n isFCL＝nFCL+1, t ═ t + 1. If n isFCL≤nFCL,maxThen correct the nodal impedance matrix ZtAnd jumping to step 2.
When the power grid normally operates, the FCL presents zero impedance or minimum impedance, rated current passes through almost without loss, and the normal operation of the power grid is not influenced; when the power grid is in fault, the short-circuit current is larger than the critical current, and the FCL presents nonlinear high-impedance characteristics in millisecond-scale time.
The MIIF index is an index used to describe the voltage interaction between adjacent dc currents. The international large grid conference organization (CIGRE) defines MIIF indexes as: when the current conversion bus m is put into a symmetrical three-phase reactor, the voltage of the bus is reducedAt 1%, the voltage change rate of the converter bus n is MIIFn,m. By definition, MIIFn,mAs shown in formula (13)
Figure BDA0001744435440000092
In the formula,. DELTA.Vm，ΔVnRespectively showing the voltage drop amount of the converter buses m and n at the moment after the symmetrical three-phase reactor is put into operation. Δ Vm，ΔVnThe values of (a) are obtained by simulation testing. MIIFn,mHas a minimum value of 0 and a maximum value of 1.
According to the MIIF definition, the calculation formula of the MIIF index based on the node admittance matrix is as shown in formula (14).
Figure BDA0001744435440000101
In the formula, ZnmRepresenting the mutual impedance between node n and node m, ZmmRepresenting the self-impedance of node m, ZnmAnd ZmmAre elements of Z in the power saving impedance matrix, Z being the inverse of the node admittance matrix Y.
The FCL optimal configuration method for reducing the risk of the secondary commutation failure provided by the invention is used for solving the optimal configuration scheme of the FCL by taking the minimum MIIF index as a target and the installation cost of the FCL as a constraint, and reducing the risk of the secondary commutation failure by reducing the MIIF index between adjacent inversion stations after the failure.
The invention also provides a fault current limiter optimized configuration system for reducing the risk of secondary commutation failure, as shown in fig. 2, the system comprises:
the objective function determining module 201 is configured to select a voltage interaction factor between a plurality of adjacent dc branches, and determine an objective function with the minimum voltage interaction factor as a target;
a fault current limiter determining module 202, configured to determine an installation position and a current limiting reactance value of a fault current limiter under a constraint condition that one fault current limiter is installed at a time and a set number of fault current limiters are installed at most;
and the voltage interaction factor calculation module 203 is configured to calculate a voltage interaction factor between adjacent direct current branches of the multi-feed-in high-voltage direct-current power transmission system after the installation position and the current-limiting reactance value of the fault current limiter are determined, skip to selecting the voltage interaction factor between the multiple adjacent direct current branches, and determine a target function step with the minimum voltage interaction factor as a target until the number of the fault current limiters reaches a set installation number.
The objective function determining module 201 specifically includes:
a multi-objective optimization function determination unit for determining a multi-objective optimization function
Figure BDA0001744435440000102
Where i denotes the number of the optimization objective,
Figure BDA0001744435440000103
and
Figure BDA0001744435440000104
the numbers of two adjacent dc branches are indicated,
Figure BDA0001744435440000105
representing branches
Figure BDA0001744435440000106
And branch
Figure BDA0001744435440000107
Voltage interaction factor between;
an objective function determination unit for determining an objective function based on the multi-objective optimization function
Figure BDA0001744435440000111
Wherein, aiTo represent
Figure BDA0001744435440000112
The weight of (a) is determined,
Figure BDA0001744435440000113
indicates the branch after the t-1 iteration
Figure BDA0001744435440000114
And branch
Figure BDA0001744435440000115
Voltage interaction factor of.
The voltage interaction factor calculation module 202 specifically includes:
an initial condition determining unit for determining an initial iteration number t equal to 1 and an initial number n of fault current limitersFCL＝1；
A first constraint condition determination unit for determining a first constraint condition
Figure BDA0001744435440000116
Wherein N isLAs the number of branches, a branch information matrix is
Figure BDA0001744435440000117
Is a variable of 0, 1, lcA number indicating the preselected installation leg of the fault current limiter,
Figure BDA0001744435440000118
indicating fault current limiter installed in first branch setcOn the branch, only the installation position of one fault current limiter is determined in each iteration, only one fault current limiter is installed,
Figure BDA0001744435440000119
is the l-th of CLcLine, represents the lcThe information of the branch is then transmitted to the subscriber,
Figure BDA00017444354400001110
Figure BDA00017444354400001111
respectively represent the lcA head end node, a tail end node and a branch impedance of the branch,
Figure BDA00017444354400001112
Figure BDA00017444354400001113
node impedance matrixes Z respectively corresponding to t-th iterationtIs located at
Figure BDA00017444354400001114
Respectively represent nodes
Figure BDA00017444354400001115
Self-impedance, node
Figure BDA00017444354400001116
Self-impedance, node
Figure BDA00017444354400001117
And node
Figure BDA00017444354400001118
The mutual impedance between the two electrodes is high,
Figure BDA00017444354400001119
are each ZtIs located at
Figure BDA00017444354400001120
Figure BDA00017444354400001121
Respectively represent nodes
Figure BDA00017444354400001122
And node
Figure BDA00017444354400001123
Mutual impedance, node
Figure BDA00017444354400001124
And node
Figure BDA00017444354400001125
Mutual impedance, node
Figure BDA00017444354400001126
And node
Figure BDA00017444354400001127
Mutual impedance, node
Figure BDA00017444354400001128
And node
Figure BDA00017444354400001129
The mutual impedance of (a);
a second constraint condition determination unit for determining a second constraint condition
Figure BDA0001744435440000121
Wherein the content of the first and second substances,
Figure BDA0001744435440000122
a third constraint condition determination unit for determining a third constraint condition
Figure BDA0001744435440000123
Wherein the content of the first and second substances,
Figure BDA0001744435440000124
c＝z1+z2-2z3+zFCL，zFCL＝-(z4 2+z4jx)/(jx), x being the reactance value of the fault current limiter, j being the imaginary sign;
a fourth constraint condition determination unit for determining that the fourth constraint condition x is less than or equal to xmaxWhere x is the number of fault current limiters installed, xmaxA set number of fault current limiters is installed at the maximum.
The FCL optimal configuration system for reducing the risk of the secondary commutation failure provided by the invention solves the optimal configuration scheme of the FCL by taking the minimum MIIF index as a target and the installation cost of the FCL as a constraint, and reduces the risk of the secondary commutation failure by reducing the MIIF index between adjacent inversion stations after the failure.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.

Claims (6)

1. A method for optimally configuring a fault current limiter to reduce the risk of a secondary commutation failure, the method comprising:
selecting voltage interaction factors among a plurality of adjacent direct current branches, and determining a target function by taking the minimum voltage interaction factor as a target;
determining the installation position and the current-limiting reactance value of the fault current limiter by taking one fault current limiter installed at each time and the maximum number of fault current limiters installed at set times as constraint conditions;
after the installation position and the current-limiting reactance value of the fault current limiter are determined, calculating voltage interaction factors between adjacent direct current branches of the multi-feed-in high-voltage direct-current power transmission system, jumping to the step of selecting the voltage interaction factors between the multiple adjacent direct current branches, determining a target function step by taking the minimum voltage interaction factor as a target, and stopping iteration until the fault current limiter reaches the set installation number.
2. The method for optimally configuring a fault current limiter for reducing the risk of the secondary commutation failure according to claim 1, wherein the determining an objective function with the minimum voltage interaction factor as a target specifically comprises:
determining a multi-objective optimization function
Figure FDA0002405486150000011
Where i denotes the number of the optimization objective,
Figure FDA0002405486150000012
Figure FDA0002405486150000013
and
Figure FDA0002405486150000014
the numbers of two adjacent dc branches are indicated,
Figure FDA0002405486150000015
representing branches
Figure FDA0002405486150000016
And branch
Figure FDA0002405486150000017
Voltage interaction factor between;
determining an objective function from the multi-objective optimization function
Figure FDA0002405486150000018
Wherein, aiTo represent
Figure FDA0002405486150000019
The weight of (a) is determined,
Figure FDA00024054861500000110
indicates the branch after the t-1 iteration
Figure FDA00024054861500000111
And branch
Figure FDA00024054861500000112
Voltage interaction factor of.
3. The method for optimally configuring a fault current limiter to reduce the risk of a secondary commutation failure according to claim 1, wherein the determining the installation position and the current-limiting reactance value of the fault current limiter is performed under the constraint that one fault current limiter is installed at a time and a set number of fault current limiters are installed at most, and specifically comprises:
determining the initial iteration number t as 1 and the initial number n of fault current limitersFCL＝1；
Determining a first constraint
Figure FDA0002405486150000021
Wherein N isLAs the number of branches, a branch information matrix is
Figure FDA0002405486150000022
Figure FDA0002405486150000023
Is a variable of 0, 1, lcA number indicating the preselected installation leg of the fault current limiter,
Figure FDA0002405486150000024
indicating fault current limiter installed in first branch setcOn the branch, only the installation position of one fault current limiter is determined in each iteration, only one fault current limiter is installed,
Figure FDA0002405486150000025
is the l-th of CLcLine, represents the lcThe information of the branch is then transmitted to the subscriber,
Figure FDA0002405486150000026
Figure FDA0002405486150000027
respectively represent the lcA head end node, a tail end node and a branch impedance of the branch,
Figure FDA0002405486150000028
Figure FDA0002405486150000029
node impedance matrixes Z respectively corresponding to t-th iterationtIs located at
Figure FDA00024054861500000210
Figure FDA00024054861500000211
Respectively represent nodes
Figure FDA00024054861500000212
Self-impedance, node
Figure FDA00024054861500000213
Self-impedance, node
Figure FDA00024054861500000214
And node
Figure FDA00024054861500000215
The mutual impedance between the two electrodes is high,
Figure FDA00024054861500000216
are each ZtIs located at
Figure FDA00024054861500000217
Respectively represent nodes
Figure FDA00024054861500000218
And node
Figure FDA0002405486150000031
Mutual impedance, node
Figure FDA0002405486150000032
And node
Figure FDA0002405486150000033
Mutual impedance, node
Figure FDA0002405486150000034
And node
Figure FDA0002405486150000035
Mutual impedance, node
Figure FDA0002405486150000036
And node
Figure FDA0002405486150000037
The mutual impedance of (a);
determining a second constraint
Figure FDA0002405486150000038
Wherein, in the step (A),
Figure FDA0002405486150000039
determining a third constraint
Figure FDA00024054861500000310
Wherein the content of the first and second substances,
Figure FDA00024054861500000311
c＝z1+z2-2z3+zFCL，zFCL＝-(z4 2+z4jx)/(jx), x being the reactance value of the fault current limiter, j being the imaginary sign;
determining that the fourth constraint y is less than or equal to ymaxWherein y is fault current limitNumber of devices mounted, ymaxA set number of fault current limiters is installed at the maximum.
4. A fault current limiter optimized configuration system for reducing the risk of a secondary commutation failure, the system comprising:
the target function determining module is used for selecting voltage interaction factors among a plurality of adjacent direct current branches, and determining a target function by taking the minimum voltage interaction factor as a target;
the fault current limiter determining module is used for determining the installation position and the current limiting reactance value of the fault current limiter by taking the condition that one fault current limiter is installed at a time and the set number of fault current limiters are installed at most as constraint conditions;
and the voltage interaction factor calculation module is used for calculating voltage interaction factors between adjacent direct current branches of the multi-feed-in high-voltage direct-current power transmission system after the installation position and the current-limiting reactance value of the fault current limiter are determined, jumping to the step of selecting the voltage interaction factors between the multiple adjacent direct current branches, and determining a target function step by taking the minimum voltage interaction factor as a target until the fault current limiter reaches the set installation number.
5. The fault current limiter optimized configuration system for reducing the risk of the secondary commutation failure according to claim 4, wherein the objective function determination module specifically comprises:
a multi-objective optimization function determination unit for determining a multi-objective optimization function
Figure FDA0002405486150000041
Where i denotes the number of the optimization objective,
Figure FDA0002405486150000042
Figure FDA0002405486150000043
and
Figure FDA0002405486150000044
the numbers of two adjacent dc branches are indicated,
Figure FDA0002405486150000045
representing branches
Figure FDA0002405486150000046
And branch
Figure FDA0002405486150000047
Voltage interaction factor between;
an objective function determination unit for determining an objective function based on the multi-objective optimization function
Figure FDA0002405486150000048
Wherein, aiTo represent
Figure FDA0002405486150000049
The weight of (a) is determined,
Figure FDA00024054861500000410
indicates the branch after the t-1 iteration
Figure FDA00024054861500000411
And branch
Figure FDA00024054861500000412
Voltage interaction factor of.
6. The fault current limiter optimized configuration system for reducing the risk of the secondary commutation failure according to claim 4, wherein the voltage interaction factor calculation module specifically comprises:
an initial condition determining unit for determining an initial iteration number t equal to 1 and an initial number n of fault current limitersFCL＝1；
A first constraint condition determination unit for determiningFirst constraint condition
Figure FDA00024054861500000413
Wherein N isLAs the number of branches, a branch information matrix is
Figure FDA00024054861500000414
Figure FDA00024054861500000415
Is a variable of 0, 1, lcA number indicating the preselected installation leg of the fault current limiter,
Figure FDA00024054861500000416
indicating fault current limiter installed in first branch setcOn the branch, only the installation position of one fault current limiter is determined in each iteration, only one fault current limiter is installed,
Figure FDA00024054861500000417
is the l-th of CLcLine, represents the lcThe information of the branch is then transmitted to the subscriber,
Figure FDA0002405486150000051
Figure FDA0002405486150000052
respectively represent the lcA head end node, a tail end node and a branch impedance of the branch,
Figure FDA0002405486150000053
Figure FDA0002405486150000054
node impedance matrixes Z respectively corresponding to t-th iterationtIs located at
Figure FDA0002405486150000055
Figure FDA0002405486150000056
Respectively represent nodes
Figure FDA0002405486150000057
Self-impedance, node
Figure FDA0002405486150000058
Self-impedance, node
Figure FDA0002405486150000059
And node
Figure FDA00024054861500000510
The mutual impedance between the two electrodes is high,
Figure FDA00024054861500000511
are each ZtIs located at
Figure FDA00024054861500000512
Respectively represent nodes
Figure FDA00024054861500000513
And node
Figure FDA00024054861500000514
Mutual impedance, node
Figure FDA00024054861500000515
And node
Figure FDA00024054861500000516
Mutual impedance, node
Figure FDA00024054861500000517
And node
Figure FDA00024054861500000518
Mutual impedance, node
Figure FDA00024054861500000519
And node
Figure FDA00024054861500000520
The mutual impedance of (a);
a second constraint condition determination unit for determining a second constraint condition
Figure FDA00024054861500000521
Wherein, in the step (A),
Figure FDA00024054861500000522
a third constraint condition determination unit for determining a third constraint condition
Figure FDA00024054861500000523
Wherein the content of the first and second substances,
Figure FDA00024054861500000524
c＝z1+z2-2z3+zFCL，zFCL＝-(z4 2+z4jx)/(jx), x being the reactance value of the fault current limiter, j being the imaginary sign;
a fourth constraint condition determination unit for determining that the fourth constraint condition y is less than or equal to ymaxWhere y is the number of fault current limiters installed, ymaxA set number of fault current limiters is installed at the maximum.
CN201810835675.6A 2018-07-26 2018-07-26 FCL optimal configuration method and system for reducing risk of secondary commutation failure Active CN108879754B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810835675.6A CN108879754B (en) 2018-07-26 2018-07-26 FCL optimal configuration method and system for reducing risk of secondary commutation failure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810835675.6A CN108879754B (en) 2018-07-26 2018-07-26 FCL optimal configuration method and system for reducing risk of secondary commutation failure

Publications (2)

Publication Number Publication Date
CN108879754A CN108879754A (en) 2018-11-23
CN108879754B true CN108879754B (en) 2020-05-01

Family

ID=64305447

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810835675.6A Active CN108879754B (en) 2018-07-26 2018-07-26 FCL optimal configuration method and system for reducing risk of secondary commutation failure

Country Status (1)

Country Link
CN (1) CN108879754B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110176755B (en) * 2019-06-18 2021-03-02 广州供电局有限公司 Fault current limiter configuration method and device
CN111884206B (en) * 2020-07-15 2021-09-17 山东大学 Key line impedance adjusting method and system for inhibiting multi-direct-current commutation failure linkage

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102761123A (en) * 2012-07-11 2012-10-31 浙江大学 Multi-infeed HVDC (high-voltage direct current) system partitioning method based on mounting of fault current limiters
CN103812129A (en) * 2014-01-16 2014-05-21 东北电力大学 Method for determining commutation failure of multi-infeed direct current power transmission system
CN104934982A (en) * 2015-05-21 2015-09-23 广东电网有限责任公司电网规划研究中心 Direct current partitioning method for multi-direct current feed-in system
CN105207276A (en) * 2015-10-10 2015-12-30 国网技术学院 Multi-infeed power transmission system fault partition current limiter installation method based on graph theory
CN105429187A (en) * 2015-12-07 2016-03-23 中国电力科学研究院 Multi-infeed direct current evaluation method
CN105939020A (en) * 2016-05-14 2016-09-14 国家电网公司 Station-configuration method for dynamic reactive compensation apparatus capable of improving multi-feed direct current restoration capability
CN106845848A (en) * 2017-02-03 2017-06-13 许继电气股份有限公司 The appraisal procedure and device of reciprocal effect between Multi-infeed HVDC transmission system current conversion station
CN108173264A (en) * 2017-12-29 2018-06-15 国电南瑞科技股份有限公司 A kind of Optimal Configuration Method of fault current limiter

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102761123A (en) * 2012-07-11 2012-10-31 浙江大学 Multi-infeed HVDC (high-voltage direct current) system partitioning method based on mounting of fault current limiters
CN103812129A (en) * 2014-01-16 2014-05-21 东北电力大学 Method for determining commutation failure of multi-infeed direct current power transmission system
CN104934982A (en) * 2015-05-21 2015-09-23 广东电网有限责任公司电网规划研究中心 Direct current partitioning method for multi-direct current feed-in system
CN105207276A (en) * 2015-10-10 2015-12-30 国网技术学院 Multi-infeed power transmission system fault partition current limiter installation method based on graph theory
CN105429187A (en) * 2015-12-07 2016-03-23 中国电力科学研究院 Multi-infeed direct current evaluation method
CN105939020A (en) * 2016-05-14 2016-09-14 国家电网公司 Station-configuration method for dynamic reactive compensation apparatus capable of improving multi-feed direct current restoration capability
CN106845848A (en) * 2017-02-03 2017-06-13 许继电气股份有限公司 The appraisal procedure and device of reciprocal effect between Multi-infeed HVDC transmission system current conversion station
CN108173264A (en) * 2017-12-29 2018-06-15 国电南瑞科技股份有限公司 A kind of Optimal Configuration Method of fault current limiter

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Restoring recloser-fuse coordination by optimal fault current limiters planning in DG-integrated distribution systems;A. Elmitwally;《Electrical Power and Energy Systems》;20160531;第77卷;第9-18页 *
基于多馈入交互作用因子的直流分区原则和方法研究;徐蔚 等;《电力系统保护与控制》;20160501;第44卷(第9期);第108-116页 *
采用多馈入交互作用因子判断高压直流系统换相失败的方法;邵瑶，汤涌;《中国电机工程学报》;20120205;第32卷(第4期);第108-114页 *

Also Published As

Publication number Publication date
CN108879754A (en) 2018-11-23

Similar Documents

Publication Publication Date Title
CN103258299B (en) A kind of many direct currents concentrate the receiving end Net Frame of Electric Network optimization method of feed-in
Biswas et al. Optimum distributed generation placement with voltage sag effect minimization
Milanovic et al. Modeling of FACTS devices for voltage sag mitigation studies in large power systems
CN107394807B (en) Multi-feed-in direct current commutation failure rapid evaluation method
CN110212565B (en) Direct-current power control strategy for inhibiting continuous commutation failure
CN104598725B (en) The evaluation method of power distribution network adaptability after a kind of distributed power source is grid-connected
CN107658895B (en) Threshold value setting method for direct-current commutation failure prediction link
CN108879754B (en) FCL optimal configuration method and system for reducing risk of secondary commutation failure
Yoosefian et al. Protection of distribution network considering fault ride through requirements of wind parks
CN115566722A (en) Distributed photovoltaic access capacity calculation method considering distribution network protection adaptability
CN107546734A (en) A kind of voltage-type DC power flow controller no-load voltage ratio optimized calculation method
Nayak et al. AC/DC interactions in multi-infeed HVDC scheme: a case study
CN110176755B (en) Fault current limiter configuration method and device
CN104993466A (en) Cascading failure fast dynamic simulation method suitable for alternating current-direct current power grid
CN102542355A (en) Operation planning part load method
CN108988401B (en) Method for evaluating power transmission capacity and constraint conditions of multi-channel alternating current-direct current series-parallel connection section
CN112787351B (en) Coordination control method for low-voltage current limiting link of multi-loop direct current system
Kong et al. Optimization of current breaker and fault current limiter in DC micro-grid based on faulty transient analysis
CN108599226B (en) True bipolar MMC-HVDC system line overload emergency control method
Marvik et al. Analysis of fault detection and location in medium voltage radial networks with distributed generation
CN116054224A (en) Quick evaluation method for commutation failure of multi-feed direct-current transmission system
Haghifam et al. Enhancement circuit breaker reliability by using fault current limiter
Qi et al. Optimal-probabilistic coordination of directional overcurrent relays considering network topological uncertainties
CN114884035A (en) Discrete type grading optimization method for parameters of multi-terminal power grid current limiting equipment
CN107623328B (en) Method for reducing closed loop current of power distribution 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
GR01 Patent grant
GR01 Patent grant