[go: up one dir, main page]

CN109245097B - A method and system for calculating AC-DC system voltage coupling action factor based on extended Jacobian matrix - Google Patents

A method and system for calculating AC-DC system voltage coupling action factor based on extended Jacobian matrix Download PDF

Info

Publication number
CN109245097B
CN109245097B CN201811237027.7A CN201811237027A CN109245097B CN 109245097 B CN109245097 B CN 109245097B CN 201811237027 A CN201811237027 A CN 201811237027A CN 109245097 B CN109245097 B CN 109245097B
Authority
CN
China
Prior art keywords
voltage
equation
bus
power
load
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
CN201811237027.7A
Other languages
Chinese (zh)
Other versions
CN109245097A (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.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Anhui Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Anhui 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 State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Anhui Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201811237027.7A priority Critical patent/CN109245097B/en
Publication of CN109245097A publication Critical patent/CN109245097A/en
Application granted granted Critical
Publication of CN109245097B publication Critical patent/CN109245097B/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
    • 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

本发明提供一种基于扩展雅克比矩阵计算交直流系统电压耦合作用因子的方法和系统，所述方法和系统基于交直流系统的每条母线的功率平衡方程，建立基于扩展雅克比矩阵的确定电压耦合作用因子的方程式X，再根据采集的所述多馈入交直流系统每条母线的网络参数值和动态元件参数值求解扩展雅克比矩阵中每个元素的值，并基于扩展雅克比矩阵中每个元素的值，采用稀疏技术对方程式X求解，确定每条直流输电系统换流母线对所述故障的受端交流母线的电压耦合作用因子，所述方法和系统使多直流落点受端电网支撑能力的快速评估结果更为准确，将本发明应用于交直流系统规划、设计及运行等领域，对于指导大电网的规划和维护大电网的安全稳定运行具有重要意义。

Figure 201811237027

The present invention provides a method and system for calculating the voltage coupling action factor of an AC-DC system based on an extended Jacobian matrix. Equation X of the coupling action factor, and then solve the value of each element in the extended Jacobian matrix according to the collected network parameter values and dynamic element parameter values of each bus of the multi-feed AC-DC system, and based on the extended Jacobian matrix For the value of each element, the sparse technique is used to solve the equation X to determine the voltage coupling action factor of each DC transmission system commutation bus to the AC bus at the receiving end of the fault. The quick evaluation result of the support capacity of the power grid is more accurate, and the application of the present invention in the fields of AC and DC system planning, design and operation is of great significance for guiding the planning of the large power grid and maintaining the safe and stable operation of the large power grid.

Figure 201811237027

Description

Method and system for calculating voltage coupling action factor of alternating current-direct current system based on extended Jacobian matrix
Technical Field
The present invention relates to the field of power system planning and operation, and more particularly, to a method and system for calculating a voltage coupling effect factor of an ac/dc system based on an extended jacobian matrix.
Background
In theoretical research and engineering application of an alternating current-direct current system, concepts of a short-circuit ratio and an effective short-circuit ratio are generally adopted to evaluate a relative strength relation between an alternating current system and a direct current system, a direct current maximum transmission power, an overvoltage level and a possible resonant frequency. In a multi-direct-current-drop-point system, direct-current systems are mutually coupled and influence an alternating-current power grid, and the problem of mutual influence among multi-loop direct-current systems cannot be considered in the traditional method. Aiming at the defects, indexes such as a multi-feed short-circuit ratio, a multi-feed interaction factor and the like are provided in the field of electric power research, and the method can be used for preliminary and rapid evaluation of the supporting capability of a multi-direct-current-drop-point receiving-end power grid. In the prior art, a calculation formula for deriving a critical multi-feed interaction factor based on a minimum extinction angle judgment standard provides a new idea for quickly judging whether a multi-feed direct-current system can simultaneously generate phase change failure, but due to the limitation of the definition of the multi-feed interaction factor, the method is only suitable for the condition that a direct-current inversion side current conversion bus has failure. Therefore, under the initiation of the multi-feeding interaction factor definition, the prior art also proposes the concept of the ac/dc system voltage coupling interaction factor, which is defined as:
if a three-phase symmetric inductive ground fault occurs in the ac bus k of the receiving side system and the voltage on the bus decreases, the voltage variation of the inverter-side converter bus j is Δ UjDefining AC/DC voltage coupling factor ADVCFjkThe calculation formula of (2) is as follows:
ADVCFjk＝△Uj/△Uk
the method for calculating the AC/DC voltage coupling effect factor replaces the ratio of the voltage amplitude variation by the ratio of elements in the equivalent impedance matrix of the receiving end AC system, and is equivalent under the condition that the power network is pure and the load is also pure. However, the method makes certain assumptions and does not consider the influence of the characteristics of dynamic elements of the power system.
Disclosure of Invention
Aiming at the technical problem that the method for calculating the AC/DC voltage coupling effect factor in the prior art does not consider the influence of the characteristics of dynamic elements of a power system, the invention provides a method for calculating the AC/DC voltage coupling effect factor based on an expanded Jacobian matrix, which comprises the following steps:
determining a power balance equation of each bus in a multi-feed-in alternating current and direct current system, and when one receiving end alternating current bus of the multi-feed-in alternating current and direct current system has fault disturbance, deriving the voltage of the faulted bus on two sides of the power balance equation, and establishing an equation X for determining a voltage coupling action factor of each direct current transmission system converter bus in the multi-feed-in alternating current and direct current system to the faulted receiving end alternating current bus based on an extended Jacobian matrix;
acquiring a network parameter value and a dynamic element parameter value of each bus of the multi-feed AC/DC system, and solving a value of each element in an extended Jacobian matrix according to the parameter values, wherein the network parameter comprises an effective voltage value of each bus in the multi-feed AC/DC system, and conductance, susceptance and voltage angle difference between the buses, the dynamic element parameter is a parameter for determining power of a dynamic element in the multi-feed AC/DC system, and the dynamic element comprises a generator, a load, a DC converter and a dynamic reactive power compensation device;
and solving the equation X by adopting a sparse technology based on the value of each element in the extended Jacobian matrix determined by calculation, and determining the voltage coupling action factor of each direct current transmission system converter bus on the receiving end alternating current bus with the fault.
Further, determining a power balance equation of each bus in the multi-feed-in alternating current and direct current system, and when one receiving end alternating current bus of the multi-feed-in alternating current and direct current system has fault disturbance, deriving the voltage of the faulted bus at two sides of the power balance equation, and establishing an equation X for determining a voltage coupling action factor of each direct current transmission system converter bus in the multi-feed-in alternating current and direct current system to the faulted receiving end alternating current bus based on an extended Jacobi matrix comprises:
when the multi-infeed alternating current-direct current system comprises m loops of direct current and n buses in total, the power balance equation of each bus is expressed as:
Figure GDA0003232112830000031
in the formula, delta Pi、△QiRespectively representing the active power variation and the reactive power variation injected by the node i, wherein the equations in the formula (1) are an active power equation and a reactive power equation of the node i, PGi、QGiRespectively representing the active and reactive power output, P, of the generator injection node iLi、QLiRespectively representing the active and reactive loads, P, of node iDiRepresenting the DC power, Q, of node iDiRepresenting reactive power, U, injected into node i of the DC converteri、UjRespectively representing the voltages of nodes i, j, QSiRepresenting reactive output, G, of the injection node i of the dynamic reactive power compensatorij、BijRespectively representing the conductance and susceptance between nodes i, j, thetaijRepresenting the voltage angle difference between the nodes i and j, and taking a negative sign at a rectification side and a positive sign at an inversion side of the direct-current active power in the formula (1) for the direct-current power transmission system in the multi-feed alternating-current and direct-current system;
when the k-th alternating current bus of the receiving end alternating current system breaks down, for the formula (1), the unbalance quantity delta Q appears on the left side of the reactive power equation corresponding to the k-th alternating current buskAnd Δ P of the buskAnd the power equations of other buses still keep the left term zero, and the voltage of the x-th bus of the direct-current transmission system in the multi-feed alternating-current and direct-current system is set to be effectiveValue of UxThe variation of the effective voltage value is DeltaUxThe effective value of the voltage of the kth bus of the receiving end alternating current system is UkThe variation of the effective value of voltage is DeltaUkCoefficient of delta Ux/△UkVoltage coupling action factor ADVCF of x-th conversion bus of direct current transmission system relative to k-th bus of receiving end alternating current systemxkWhen x is more than or equal to 1 and less than or equal to m and k is more than or equal to m +1 and less than or equal to n, based on the power balance equation of the formula (1), considering the influence of a dynamic element generator set, a direct current system, load characteristics and a dynamic reactive power compensation device, and enabling two sides of the formula (1) to act on UkDerivation, generating equation (2), whose expression is:
Figure GDA0003232112830000041
let the elements in the right-hand vector of equation (2)
Figure GDA0003232112830000042
Generating an expression formula (3), wherein the expression formula is as follows:
Figure GDA0003232112830000043
in formula (3), the extended jacobian matrix is:
Figure GDA0003232112830000051
when i ≠ j, the element H in Jacobian matrixij、Nij、Mij、LijThe calculation formula of (a) is as follows:
Figure GDA0003232112830000052
when i ═ j, element H in Jacobian matrixii、Nii、Mii、LiiThe calculation formula of (a) is as follows:
Figure GDA0003232112830000061
as shown in the formula (3), the Jacobian matrix is a 2 n-dimensional square matrix, the number of variables to be solved is 2n-1, wherein the 2 k-th row corresponding to the reactive power equation of the bus k is a redundant row,
Figure GDA0003232112830000062
for unknowns, column 2k corresponds
Figure GDA0003232112830000063
Deleting the 2k row and moving the 2k column to the left of equation (3) to obtain equation X, which is expressed as:
Figure GDA0003232112830000064
further, the acquiring a network parameter value and a dynamic element parameter value of each bus of the multi-infeed alternating current-direct current system, and solving a value of each element in an extended jacobian matrix according to the parameter values includes:
when i is not equal to j, determining an element H in the Jacobian matrix according to the collected network parameter value of each bus of the multi-feed AC/DC system and the formula (4)ij、Nij、Mij、LijA value of (d);
when i is j, determining an element H in the Jacobian matrix according to the acquired network parameter value, dynamic element parameter value and formula (5) of each bus of the multi-feed AC/DC systemii、Nii、Mii、LiiA value of (a), wherein:
derivative terms of the associated generator in equation (5) for non-generator nodes
Figure GDA0003232112830000071
Zero, at the generator node, approximately considers the generator sub-transient reactance X at the disturbance instant "dThe back electromotive force E' is kept constant, the expression of the output power of the generator is an expression (6),wherein theta isδFor the generator internal potential E' and terminal voltage UiThe derivative term of the generator power to voltage is the equation (7):
Figure GDA0003232112830000072
Figure GDA0003232112830000073
for the non-load node, the derivative term of the load power to the voltage in equation (5)
Figure GDA0003232112830000074
Zero, and in addition the load power is dependent only on the effective value of the feed point voltage and not on its angle, so the derivative of the load power with respect to the voltage angle
Figure GDA0003232112830000075
Is zero;
for a load node, when the load is a constant power load, the constant current load power expression and its derivative to voltage are equation (8):
Figure GDA0003232112830000081
the direct current converter voltage and current equation expressed by the named value is an equation (9), the converter power equation derived from the equation (9) is an equation (10), wherein the voltage variable of the converter bus in the equation (10) is only UiEquation (9) and equation (10) are as follows, excluding the voltage angle:
Figure GDA0003232112830000082
Figure GDA0003232112830000083
in the formula of UdRepresenting a direct voltage, ntRepresents the number of the six-pulse current converters connected in series, kTRepresenting the transformer ratio, theta, of the converterdRepresenting the DC commutation angle of the rectifier or the extinction angle of the inverter, XcDenotes the equivalent commutation reactance, IdRepresenting a direct current, kγWhich represents the equivalent transformation ratio of the converter transformer,
Figure GDA0003232112830000084
representing the equivalent power factor angle, IiRepresents the current injected into the ac system by dc;
when the rectification side of the direct current transmission system adopts constant current control and the inversion side adopts constant extinction angle control, the derivative of the power, which is extracted from the alternating current system by the inverter side converter of the direct current system, to the voltage is an equation (11):
Figure GDA0003232112830000091
when the dynamic reactive power compensation device is a static reactive power compensator, the static reactive power compensator adopts the voltage deviation of a controlled bus as an input signal, controls the equivalent susceptance of the compensation device through a proportional amplification link, neglects a delay link, and outputs the relation between the equivalent susceptance and the voltage deviation as an expression (12), and the derivative of the voltage as an expression (13):
Bi＝-K△Ui＝-K(Ui-Ui0) (12)
Figure GDA0003232112830000092
further, when i is equal to j, determining an element H in a Jacobian matrix according to the collected network parameter value, dynamic element parameter value and formula (5) of each bus of the multi-feed AC/DC systemii、Nii、Mii、LiiThe values of (a) further include:
for a load node, when the load is a constant impedance load, the expression of the load power and its derivative to voltage is equation (14):
Figure GDA0003232112830000101
when the rectification side of the direct current transmission system adopts constant power control and the inversion side adopts constant extinction angle control, applying small voltage fluctuation on a converter bus, applying a control strategy, calculating power change of converter stations at two sides, and applying difference to replace partial differential;
when the dynamic reactive power compensation device is a static synchronous compensator, the voltage deviation of the controlled bus is used as an input signal, and the steady state equation can be expressed as follows:
△U＝UREF-U＝KDIS (15)
IS＝△U/KD＝BS△U (16)
the equivalent susceptance of the static synchronous compensator is denoted BS＝1/KDWhen K isDWhen zero is taken, the control node of the static synchronous compensator is in no-difference control, but is limited by the output current of the static synchronous compensator, and the derivative of the voltage is expressed by the formula (17):
Figure GDA0003232112830000102
according to another aspect of the present invention, there is provided a system for calculating a voltage coupling effect factor of an ac/dc system based on an extended jacobian matrix, the system comprising:
the system comprises an equation determining unit, a fault detection unit and a fault detection unit, wherein the equation determining unit is used for determining a power balance equation of each bus in the multi-feed-in alternating current and direct current system, and when one receiving end alternating current bus of the multi-feed-in alternating current and direct current system has fault disturbance, deriving the voltage of the faulted bus on two sides of the power balance equation, and establishing an equation X for determining a voltage coupling action factor of each direct current transmission system converter bus in the multi-feed-in alternating current and direct current system to the faulted receiving end alternating current bus based on an extended Jacobi matrix;
the element calculation unit is used for acquiring a network parameter value and a dynamic element parameter value of each bus of the multi-feed AC/DC system, and solving a value of each element in an extended Jacobian matrix according to the parameter values, wherein the network parameter comprises an effective voltage value of each bus in the multi-feed AC/DC system, and conductance, susceptance and voltage angle difference between the buses, the dynamic element parameter is a parameter for determining power of a dynamic element in the multi-feed AC/DC system, and the dynamic element comprises a generator, a load, a DC converter and a dynamic reactive power compensation device;
and the factor determination unit is used for solving the equation X by adopting a sparse technology based on the value of each element in the extended Jacobian matrix determined by calculation, and determining the voltage coupling action factor of each direct current transmission system conversion bus on the receiving end alternating current bus of the fault.
Further, the equation determining unit determines a power balance equation of each bus in the multi-feed ac/dc system, and when a receiving-end ac bus of the multi-feed ac/dc system has a fault disturbance, derives the voltage of the faulty bus on both sides of the power balance equation, and establishes an equation X for determining a voltage coupling effect factor of each dc power transmission system converter bus in the multi-feed ac/dc system on the faulty receiving-end ac bus based on an extended jacobian matrix, including:
when the multi-infeed alternating current-direct current system comprises m loops of direct current and n buses in total, the power balance equation of each bus is expressed as:
Figure GDA0003232112830000111
in the formula, delta Pi、△QiRespectively representing the active power variation and the reactive power variation injected by the node i, wherein the equations in the formula (1) are an active power equation and a reactive power equation of the node i, PGi、QGiRespectively representing the active output of the generator injection node iForce and reactive force, PLi、QLiRespectively representing the active and reactive loads, P, of node iDiRepresenting the DC power, Q, of node iDiRepresenting reactive power, U, injected into node i of the DC converteri、UjRespectively representing the voltages of nodes i, j, QSiRepresenting reactive output, G, of the injection node i of the dynamic reactive power compensatorij、BijRespectively representing the conductance and susceptance between nodes i, j, thetaijRepresenting the voltage angle difference between the nodes i and j, and taking a negative sign at a rectification side and a positive sign at an inversion side of the direct-current active power in the formula (1) for the direct-current power transmission system in the multi-feed alternating-current and direct-current system;
when the k-th alternating current bus of the receiving end alternating current system breaks down, for the formula (1), the unbalance quantity delta Q appears on the left side of the reactive power equation corresponding to the k-th alternating current buskAnd Δ P of the buskAnd the power equations of other buses still keep the left term zero, and the voltage effective value of the x-th bus of the direct-current power transmission system in the multi-feed alternating-current and direct-current system is set to be UxWhen the voltage is applied, the change quantity of the effective voltage value is delta UxThe effective value of the voltage of the kth bus of the receiving end alternating current system is UkThe variation of the effective value of voltage is DeltaUkCoefficient of delta Ux/△UkVoltage coupling action factor ADVCF of x-th conversion bus of direct current transmission system relative to k-th bus of receiving end alternating current systemxkWhen x is more than or equal to 1 and less than or equal to m and k is more than or equal to m +1 and less than or equal to n, based on the power balance equation of the formula (1), considering the influence of a dynamic element generator set, a direct current system, load characteristics and a dynamic reactive power compensation device, and enabling two sides of the formula (1) to act on UkDerivation, generating equation (2), whose expression is:
Figure GDA0003232112830000121
let the elements in the right-hand vector of equation (2)
Figure GDA0003232112830000122
The resultant formula (3), table thereofThe expression is as follows:
Figure GDA0003232112830000131
in formula (3), the extended jacobian matrix is:
Figure GDA0003232112830000132
when i ≠ j, the element H in Jacobian matrixij、Nij、Mij、LijThe calculation formula of (a) is as follows:
Figure GDA0003232112830000141
when i ═ j, element H in Jacobian matrixii、Nii、Mii、LiiThe calculation formula of (a) is as follows:
Figure GDA0003232112830000142
as shown in the formula (3), the Jacobian matrix is a 2 n-dimensional square matrix, the number of variables to be solved is 2n-1, wherein the 2 k-th row corresponding to the reactive power equation of the bus k is a redundant row,
Figure GDA0003232112830000143
for unknowns, column 2k corresponds
Figure GDA0003232112830000144
Deleting the 2k row and moving the 2k column to the left of equation (3) to obtain equation X, which is expressed as:
Figure GDA0003232112830000151
further, the acquiring, by the element calculation unit, a network parameter value and a dynamic element parameter value of each bus of the multi-feed ac/dc system, and solving a value of each element in the extended jacobian matrix according to the parameter values includes:
when i is not equal to j, determining an element H in the Jacobian matrix according to the collected network parameter value of each bus of the multi-feed AC/DC system and the formula (4)ij、Nij、Mij、LijA value of (d);
when i is j, determining an element H in the Jacobian matrix according to the acquired network parameter value, dynamic element parameter value and formula (5) of each bus of the multi-feed AC/DC systemii、Nii、Mii、LiiA value of (a), wherein:
derivative terms of the associated generator in equation (5) for non-generator nodes
Figure GDA0003232112830000152
Zero, at the generator node, approximately considers the generator sub-transient reactance X at the disturbance instant "dThe back electromotive force E' is kept constant, and the output power of the generator is expressed by the formula (6), wherein thetaδFor the generator internal potential E' and terminal voltage UiThe derivative term of the generator power to voltage is the equation (7):
Figure GDA0003232112830000153
Figure GDA0003232112830000161
for the non-load node, the derivative term of the load power to the voltage in equation (5)
Figure GDA0003232112830000162
Zero, and in addition the load power is dependent only on the effective value of the feed point voltage and not on its angle, so the derivative of the load power with respect to the voltage angle
Figure GDA0003232112830000163
Is zero;
for a load node, when the load is a constant power load, the constant current load power expression and its derivative to voltage are equation (8):
Figure GDA0003232112830000164
the direct current converter voltage and current equation expressed by the named value is an equation (9), the converter power equation derived from the equation (9) is an equation (10), wherein the voltage variable of the converter bus in the equation (10) is only UiEquation (9) and equation (10) are as follows, excluding the voltage angle:
Figure GDA0003232112830000171
Figure GDA0003232112830000172
in the formula of UdRepresenting a direct voltage, ntRepresents the number of the six-pulse current converters connected in series, kTRepresenting the transformer ratio, theta, of the converterdRepresenting the DC commutation angle of the rectifier or the extinction angle of the inverter, XcDenotes the equivalent commutation reactance, IdRepresenting a direct current, kγWhich represents the equivalent transformation ratio of the converter transformer,
Figure GDA0003232112830000174
representing the equivalent power factor angle, IiRepresents the current injected into the ac system by dc;
when the rectification side of the direct current transmission system adopts constant current control and the inversion side adopts constant extinction angle control, the derivative of the power, which is extracted from the alternating current system by the inverter side converter of the direct current system, to the voltage is an equation (11):
Figure GDA0003232112830000173
when the dynamic reactive power compensation device is a static reactive power compensator, the static reactive power compensator adopts the voltage deviation of a controlled bus as an input signal, controls the equivalent susceptance of the compensation device through a proportional amplification link, neglects a delay link, and outputs the relation between the equivalent susceptance and the voltage deviation as an expression (12), and the derivative of the voltage as an expression (13):
Bi＝-K△Ui＝-K(Ui-Ui0) (12)
Figure GDA0003232112830000181
further, when i is j, the element calculation unit determines an element H in the jacobian matrix according to the collected network parameter value and dynamic element parameter value of each bus of the multi-feed ac/dc system and the formula (5)ii、Nii、Mii、LiiThe values of (a) further include:
for a load node, when the load is a constant impedance load, the expression of the load power and its derivative to voltage is equation (14):
Figure GDA0003232112830000182
when the rectification side of the direct current transmission system adopts constant power control and the inversion side adopts constant extinction angle control, applying small voltage fluctuation on a converter bus, applying a control strategy, calculating power change of converter stations at two sides, and applying difference to replace partial differential;
when the dynamic reactive power compensation device is a static synchronous compensator, the voltage deviation of the controlled bus is used as an input signal, and the steady state equation can be expressed as follows:
△U＝UREF-U＝KDIS (15)
IS＝△U/KD＝BS△U (16)
to static synchronise the compensatorThe equivalent susceptance is represented as BS＝1/KDWhen K isDWhen zero is taken, the control node of the static synchronous compensator is in no-difference control, but is limited by the output current of the static synchronous compensator, and the derivative of the voltage is expressed by the formula (17):
Figure GDA0003232112830000191
the method and the system for calculating the voltage coupling effect factor of the alternating current and direct current system based on the extended Jacobian matrix comprehensively consider the characteristics of dynamic element models in the power system, incorporate dynamic element models such as a generator, a load, a direct current and a dynamic reactive power compensation device into the calculation of the voltage coupling effect factor of the alternating current and direct current system, and determine the voltage coupling effect factor of each direct current transmission system conversion current bus on the faulty receiving end alternating current bus by calculating the value of each element in the extended Jacobian matrix, so that the quick evaluation result of the supporting capability of the multi-direct current drop receiving end power grid is more accurate.
Drawings
A more complete understanding of exemplary embodiments of the present invention may be had by reference to the following drawings in which:
FIG. 1 is a flow chart of a method for calculating a voltage coupling effect factor of an AC/DC system based on an extended Jacobian matrix according to a preferred embodiment of the present invention;
fig. 2 is a schematic structural diagram of a system for calculating a voltage coupling effect factor of an ac/dc system based on an extended jacobian matrix according to a preferred embodiment of the present invention.
Detailed Description
The exemplary embodiments of the present invention will now be described with reference to the accompanying drawings, however, the present invention may be embodied in many different forms and is not limited to the embodiments described herein, which are provided for complete and complete disclosure of the present invention and to fully convey the scope of the present invention to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the accompanying drawings is not intended to be limiting of the invention. In the drawings, the same units/elements are denoted by the same reference numerals.
Unless otherwise defined, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Further, it will be understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
Fig. 1 is a flowchart of a method for calculating a voltage coupling effect factor of an ac/dc system based on an extended jacobian matrix according to a preferred embodiment of the present invention. As shown in fig. 1, the method 100 for calculating the ac/dc system voltage coupling factor based on the extended jacobian matrix according to the preferred embodiment starts with step 101.
In step 101, determining a power balance equation of each bus in a multi-feed-in alternating current and direct current system, and when one receiving end alternating current bus of the multi-feed-in alternating current and direct current system has fault disturbance, deriving the voltage of the faulted bus at two sides of the power balance equation, and establishing an equation X for determining a voltage coupling action factor of each direct current transmission system converter bus in the multi-feed-in alternating current and direct current system to the faulted receiving end alternating current bus based on an extended Jacobi matrix;
acquiring a network parameter value and a dynamic element parameter value of each bus of the multi-feed alternating current and direct current system, and solving a value of each element in an extended Jacobian matrix according to the parameter values, wherein the network parameter comprises an effective voltage value of each bus in the multi-feed alternating current and direct current system, and conductance, susceptance and voltage phase difference between the buses, the dynamic element parameter is a parameter for determining power of a dynamic element in the multi-feed alternating current and direct current system, and the dynamic element comprises a generator, a load, a direct current converter and a dynamic reactive power compensation device;
in step 103, based on the calculated value of each element in the extended jacobian matrix, solving equation X by using a sparse technique, and determining a voltage coupling action factor of each dc transmission system converter bus on the faulty receiving-end ac bus.
Preferably, determining a power balance equation of each bus in the multi-feed-in ac/dc system, and when a fault disturbance occurs to one receiving-end ac bus of the multi-feed-in ac/dc system, deriving the voltage of the faulty bus at two sides of the power balance equation, and establishing an equation X for determining a voltage coupling effect factor of each dc power transmission system converter bus in the multi-feed-in ac/dc system to the faulty receiving-end ac bus based on an extended jacobian matrix includes:
when the multi-infeed alternating current-direct current system comprises m loops of direct current and n buses in total, the power balance equation of each bus is expressed as:
Figure GDA0003232112830000211
in the formula, delta Pi、△QiRespectively representing the active power variation and the reactive power variation injected by the node i, wherein the equations in the formula (1) are an active power equation and a reactive power equation of the node i, PGi、QGiRespectively representing the active and reactive power output, P, of the generator injection node iLi、QLiRespectively representing the active and reactive loads, P, of node iDiRepresenting the DC power, Q, of node iDiRepresenting reactive power, U, injected into node i of the DC converteri、UjRespectively representing the voltages of nodes i, j, QSiRepresenting reactive output, G, of the injection node i of the dynamic reactive power compensatorij、BijRespectively representing the conductance and susceptance between nodes i, j, thetaijRepresenting the voltage angle difference between the nodes i and j, and taking a negative sign at a rectification side and a positive sign at an inversion side of the direct-current active power in the formula (1) for the direct-current power transmission system in the multi-feed alternating-current and direct-current system;
when the k-th AC bus of the receiving-end AC system is in fault, regarding the formula (1), the k-th AC bus is in faultUnbalance quantity delta Q appears on the left side of reactive power equation corresponding to current buskAnd Δ P of the buskAnd the power equations of other buses still keep the left term zero, and the voltage effective value of the x-th bus of the direct-current power transmission system in the multi-feed alternating-current and direct-current system is set to be UxWhen the voltage is applied, the change quantity of the effective voltage value is delta UxThe effective value of the voltage of the kth bus of the receiving end alternating current system is UkThe variation of the effective value of voltage is DeltaUkCoefficient of delta Ux/△UkVoltage coupling action factor ADVCF of x-th conversion bus of direct current transmission system relative to k-th bus of receiving end alternating current systemxkWhen x is more than or equal to 1 and less than or equal to m and k is more than or equal to m +1 and less than or equal to n, based on the power balance equation of the formula (1), considering the influence of a dynamic element generator set, a direct current system, load characteristics and a dynamic reactive power compensation device, and enabling two sides of the formula (1) to act on UkDerivation, generating equation (2), whose expression is:
Figure GDA0003232112830000221
let the elements in the right-hand vector of equation (2)
Figure GDA0003232112830000222
Generating an expression formula (3), wherein the expression formula is as follows:
Figure GDA0003232112830000223
in formula (3), the extended jacobian matrix is:
Figure GDA0003232112830000231
when i ≠ j, the element H in Jacobian matrixij、Nij、Mij、LijThe calculation formula of (a) is as follows:
Figure GDA0003232112830000232
when i ═ j, element H in Jacobian matrixii、Nii、Mii、LiiThe calculation formula of (a) is as follows:
Figure GDA0003232112830000241
as shown in the formula (3), the Jacobian matrix is a 2 n-dimensional square matrix, the number of variables to be solved is 2n-1, wherein the 2 k-th row corresponding to the reactive power equation of the bus k is a redundant row,
Figure GDA0003232112830000242
for unknowns, column 2k corresponds
Figure GDA0003232112830000243
Deleting the 2k row and moving the 2k column to the left of equation (3) to obtain equation X, which is expressed as:
Figure GDA0003232112830000244
preferably, the acquiring a network parameter value and a dynamic element parameter value of each bus of the multi-infeed alternating current-direct current system, and solving a value of each element in an extended jacobian matrix according to the parameter values includes:
when i is not equal to j, determining an element H in the Jacobian matrix according to the collected network parameter value of each bus of the multi-feed AC/DC system and the formula (4)ij、Nij、Mij、LijA value of (d);
when i is j, determining an element H in the Jacobian matrix according to the acquired network parameter value, dynamic element parameter value and formula (5) of each bus of the multi-feed AC/DC systemii、Nii、Mii、LiiA value of (a), wherein:
derivative terms of the associated generator in equation (5) for non-generator nodes
Figure GDA0003232112830000251
Zero, at the generator node, approximately considers the generator sub-transient reactance X at the disturbance instant "dThe back electromotive force E' is kept constant, and the output power of the generator is expressed by the formula (6), wherein thetaδFor the generator internal potential E' and terminal voltage UiThe derivative term of the generator power to voltage is the equation (7):
Figure GDA0003232112830000252
Figure GDA0003232112830000253
for the non-load node, the derivative term of the load power to the voltage in equation (5)
Figure GDA0003232112830000254
Zero, and in addition the load power is dependent only on the effective value of the feed point voltage and not on its angle, so the derivative of the load power with respect to the voltage angle
Figure GDA0003232112830000255
Is zero;
for a load node, when the load is a constant power load, the constant current load power expression and its derivative to voltage are equation (8):
Figure GDA0003232112830000261
the direct current converter voltage and current equation expressed by the named value is an equation (9), the converter power equation derived from the equation (9) is an equation (10), wherein the voltage variable of the converter bus in the equation (10) is only UiEquation (9) and equation (10) are as follows, excluding the voltage angle:
Figure GDA0003232112830000262
Figure GDA0003232112830000263
in the formula of UdRepresenting a direct voltage, ntRepresents the number of the six-pulse current converters connected in series, kTRepresenting the transformer ratio, theta, of the converterdRepresenting the DC commutation angle of the rectifier or the extinction angle of the inverter, XcDenotes the equivalent commutation reactance, IdRepresenting a direct current, kγWhich represents the equivalent transformation ratio of the converter transformer,
Figure GDA0003232112830000264
representing the equivalent power factor angle, IiRepresents the current injected into the ac system by dc;
when the rectification side of the direct current transmission system adopts constant current control and the inversion side adopts constant extinction angle control, the derivative of the power, which is extracted from the alternating current system by the inverter side converter of the direct current system, to the voltage is an equation (11):
Figure GDA0003232112830000271
when the dynamic reactive power compensation device is a static reactive power compensator, the static reactive power compensator adopts the voltage deviation of a controlled bus as an input signal, controls the equivalent susceptance of the compensation device through a proportional amplification link, neglects a delay link, and outputs the relation between the equivalent susceptance and the voltage deviation as an expression (12), and the derivative of the voltage as an expression (13):
Bi＝-K△Ui＝-K(Ui-Ui0) (12)
Figure GDA0003232112830000272
preferably, when i is j, each multi-feed AC/DC system is collectedThe values of the network parameters, the dynamic element parameters of the bus, and equation (5) determine the element H in the Jacobian matrixii、Nii、Mii、LiiThe values of (a) further include:
for a load node, when the load is a constant impedance load, the expression of the load power and its derivative to voltage is equation (14):
Figure GDA0003232112830000281
when the rectification side of the direct current transmission system adopts constant power control and the inversion side adopts constant extinction angle control, applying small voltage fluctuation on a converter bus, applying a control strategy, calculating power change of converter stations at two sides, and applying difference to replace partial differential;
when the dynamic reactive power compensation device is a static synchronous compensator, the voltage deviation of the controlled bus is used as an input signal, and the steady state equation can be expressed as follows:
△U＝UREF-U＝KDIS (15)
IS＝△U/KD＝BS△U (16)
the equivalent susceptance of the static synchronous compensator is denoted BS＝1/KDWhen K isDWhen zero is taken, the control node of the static synchronous compensator is in no-difference control, but is limited by the output current of the static synchronous compensator, and the derivative of the voltage is expressed by the formula (17):
Figure GDA0003232112830000282
fig. 2 is a schematic structural diagram of a system for calculating a voltage coupling effect factor of an ac/dc system based on an extended jacobian matrix according to a preferred embodiment of the present invention. As shown in fig. 2, the system 200 for calculating the voltage coupling factor of the ac/dc system based on the extended jacobian matrix according to the preferred embodiment includes:
an equation determining unit 201, configured to determine a power balance equation of each bus in the multi-feed ac/dc system, and when a receiving-end ac bus of the multi-feed ac/dc system has a fault disturbance, derive voltages of the faulty bus on two sides of the power balance equation, and establish an equation X based on an extended jacobian matrix to determine a voltage coupling effect factor of each dc power transmission system converter bus in the multi-feed ac/dc system on the faulty receiving-end ac bus;
an element calculation unit 202, configured to collect a network parameter value and a dynamic element parameter value of each bus of the multi-feed ac/dc system, and solve a value of each element in an extended jacobian matrix according to the parameter values, where the network parameter includes an effective voltage value of each bus in the multi-feed ac/dc system, and a conductance, susceptance, and voltage angle difference between the bus and the bus, and the dynamic element parameter is a parameter that determines power of a dynamic element in the multi-feed ac/dc system, where the dynamic element includes a generator, a load, a dc converter, and a dynamic reactive power compensation device;
and the factor determining unit 203 is configured to solve equation X by using a sparse technique based on the calculated value of each element in the extended jacobian matrix, and determine a voltage coupling action factor of each dc transmission system converter bus to the faulty receiving-end ac bus.
Preferably, the equation determining unit 201 determines a power balance equation of each bus in the multi-feed ac/dc system, and when a receiving-end ac bus of the multi-feed ac/dc system has a fault disturbance, derives the voltage of the faulty bus on both sides of the power balance equation, and establishes an equation X based on an extended jacobian matrix to determine a voltage coupling effect factor of each dc power transmission system converter bus in the multi-feed ac/dc system to the faulty receiving-end ac bus, including:
when the multi-infeed alternating current-direct current system comprises m loops of direct current and n buses in total, the power balance equation of each bus is expressed as:
Figure GDA0003232112830000291
in the formula, delta Pi、△QiRespectively representing the active power variation and the reactive power variation injected by the node i, wherein the equations in the formula (1) are an active power equation and a reactive power equation of the node i, PGi、QGiRespectively representing the active and reactive power output, P, of the generator injection node iLi、QLiRespectively representing the active and reactive loads, P, of node iDiRepresenting the DC power, Q, of node iDiRepresenting reactive power, U, injected into node i of the DC converteri、UjRespectively representing the voltages of nodes i, j, QSiRepresenting reactive output, G, of the injection node i of the dynamic reactive power compensatorij、BijRespectively representing the conductance and susceptance between nodes i, j, thetaijRepresenting the voltage angle difference between the nodes i and j, and taking a negative sign at a rectification side and a positive sign at an inversion side of the direct-current active power in the formula (1) for the direct-current power transmission system in the multi-feed alternating-current and direct-current system;
when the k-th alternating current bus of the receiving end alternating current system breaks down, for the formula (1), the unbalance quantity delta Q appears on the left side of the reactive power equation corresponding to the k-th alternating current buskAnd Δ P of the buskAnd the power equations of other buses still keep the left term zero, and the voltage effective value of the x-th bus of the direct-current power transmission system in the multi-feed alternating-current and direct-current system is set to be UxWhen the voltage is applied, the change quantity of the effective voltage value is delta UxThe effective value of the voltage of the kth bus of the receiving end alternating current system is UkThe variation of the effective value of voltage is DeltaUkCoefficient of delta Ux/△UkVoltage coupling action factor ADVCF of x-th conversion bus of direct current transmission system relative to k-th bus of receiving end alternating current systemxkWhen x is more than or equal to 1 and less than or equal to m and k is more than or equal to m +1 and less than or equal to n, based on the power balance equation of the formula (1), considering the influence of a dynamic element generator set, a direct current system, load characteristics and a dynamic reactive power compensation device, and enabling two sides of the formula (1) to act on UkDerivation, generating equation (2), whose expression is:
Figure GDA0003232112830000301
let the elements in the right-hand vector of equation (2)
Figure GDA0003232112830000302
Generating an expression formula (3), wherein the expression formula is as follows:
Figure GDA0003232112830000311
in formula (3), the extended jacobian matrix is:
Figure GDA0003232112830000312
when i ≠ j, the element H in Jacobian matrixij、Nij、Mij、LijThe calculation formula of (a) is as follows:
Figure GDA0003232112830000321
when i ═ j, element H in Jacobian matrixii、Nii、Mii、LiiThe calculation formula of (a) is as follows:
Figure GDA0003232112830000322
as shown in the formula (3), the Jacobian matrix is a 2 n-dimensional square matrix, the number of variables to be solved is 2n-1, wherein the 2 k-th row corresponding to the reactive power equation of the bus k is a redundant row,
Figure GDA0003232112830000323
for unknowns, column 2k corresponds
Figure GDA0003232112830000324
Deleting the 2k row and moving the 2k column to the left of equation (3) to obtain equation X, which is expressed as:
Figure GDA0003232112830000331
preferably, the acquiring, by the element calculating unit 202, a network parameter value and a dynamic element parameter value of each bus of the multi-infeed ac/dc system, and solving a value of each element in the extended jacobian matrix according to the parameter values includes:
when i is not equal to j, determining an element H in the Jacobian matrix according to the collected network parameter value of each bus of the multi-feed AC/DC system and the formula (4)ij、Nij、Mij、LijA value of (d);
when i is j, determining an element H in the Jacobian matrix according to the acquired network parameter value, dynamic element parameter value and formula (5) of each bus of the multi-feed AC/DC systemii、Nii、Mii、LiiA value of (a), wherein:
derivative terms of the associated generator in equation (5) for non-generator nodes
Figure GDA0003232112830000332
Zero, at the generator node, approximately considers the generator sub-transient reactance X at the disturbance instant "dThe back electromotive force E' is kept constant, and the output power of the generator is expressed by the formula (6), wherein thetaδFor the generator internal potential E' and terminal voltage UiThe derivative term of the generator power to voltage is the equation (7):
Figure GDA0003232112830000333
Figure GDA0003232112830000341
for non-loaded nodes, load in equation (5)Derivative term of power versus voltage
Figure GDA0003232112830000342
Zero, and in addition the load power is dependent only on the effective value of the feed point voltage and not on its angle, so the derivative of the load power with respect to the voltage angle
Figure GDA0003232112830000343
Is zero;
for a load node, when the load is a constant power load, the constant current load power expression and its derivative to voltage are equation (8):
Figure GDA0003232112830000344
the direct current converter voltage and current equation expressed by the named value is an equation (9), the converter power equation derived from the equation (9) is an equation (10), wherein the voltage variable of the converter bus in the equation (10) is only UiEquation (9) and equation (10) are as follows, excluding the voltage angle:
Figure GDA0003232112830000351
Figure GDA0003232112830000352
in the formula of UdRepresenting a direct voltage, ntRepresents the number of the six-pulse current converters connected in series, kTRepresenting the transformer ratio, theta, of the converterdRepresenting the DC commutation angle of the rectifier or the extinction angle of the inverter, XcDenotes the equivalent commutation reactance, IdRepresenting a direct current, kγWhich represents the equivalent transformation ratio of the converter transformer,
Figure GDA0003232112830000354
representing the equivalent power factor angle, IiRepresenting dc injected into ac systemA stream;
when the rectification side of the direct current transmission system adopts constant current control and the inversion side adopts constant extinction angle control, the derivative of the power, which is extracted from the alternating current system by the inverter side converter of the direct current system, to the voltage is an equation (11):
Figure GDA0003232112830000353
when the dynamic reactive power compensation device is a static reactive power compensator, the static reactive power compensator adopts the voltage deviation of a controlled bus as an input signal, controls the equivalent susceptance of the compensation device through a proportional amplification link, neglects a delay link, and outputs the relation between the equivalent susceptance and the voltage deviation as an expression (12), and the derivative of the voltage as an expression (13):
Bi＝-K△Ui＝-K(Ui-Ui0) (12)
Figure GDA0003232112830000361
preferably, the element calculating unit 202 determines the element H in the jacobian matrix according to the collected network parameter value, dynamic element parameter value and equation (5) of each bus of the multi-feed ac/dc system when i is equal to jii、Nii、Mii、LiiThe values of (a) further include:
for a load node, when the load is a constant impedance load, the expression of the load power and its derivative to voltage is equation (14):
Figure GDA0003232112830000362
when the rectification side of the direct current transmission system adopts constant power control and the inversion side adopts constant extinction angle control, applying small voltage fluctuation on a converter bus, applying a control strategy, calculating power change of converter stations at two sides, and applying difference to replace partial differential;
when the dynamic reactive power compensation device is a static synchronous compensator, the voltage deviation of the controlled bus is used as an input signal, and the steady state equation can be expressed as follows:
△U＝UREF-U＝KDIS (15)
IS＝△U/KD＝BS△U (16)
the equivalent susceptance of the static synchronous compensator is denoted BS＝1/KDWhen K isDWhen zero is taken, the control node of the static synchronous compensator is in no-difference control, but is limited by the output current of the static synchronous compensator, and the derivative of the voltage is expressed by the formula (17):
Figure GDA0003232112830000371
the invention has been described with reference to a few embodiments. However, other embodiments of the invention than the one disclosed above are equally possible within the scope of the invention, as would be apparent to a person skilled in the art from the appended patent claims.
Generally, all terms used in the claims are to be interpreted according to their ordinary meaning in the technical field, unless explicitly defined otherwise herein. All references to "a/an/the [ device, component, etc ]" are to be interpreted openly as referring to at least one instance of said device, component, etc., unless explicitly stated otherwise. The steps of any method disclosed herein do not have to be performed in the exact order disclosed, unless explicitly stated.

Claims (6)

1.一种基于扩展雅克比矩阵计算交直流系统电压耦合作用因子的方法，其特征在于，所述方法包括：1. a method for calculating AC-DC system voltage coupling action factor based on extended Jacobian matrix, is characterized in that, described method comprises: 确定多馈入交直流系统中每条母线的功率平衡方程，并在所述多馈入交直流系统的一条受端交流母线出现故障扰动时，将所述功率平衡方程两侧对发生故障的母线的电压求导，建立基于扩展雅克比矩阵的确定多馈入交直流系统中送端直流输电系统的每条换流母线对故障的受端交流母线的电压耦合作用因子的方程式X，其中：Determine the power balance equation of each bus in the multi-feed AC/DC system, and when a receiving end AC bus of the multi-feed AC/DC system is disturbed by a fault, compare the two sides of the power balance equation to the faulty bus The voltage derivation based on the extended Jacobian matrix is established to determine the equation X of the voltage coupling action factor of each commutation busbar of the sending-end DC transmission system in the multi-infeed AC-DC transmission system to the fault-receiving AC busbar, where: 当所述多馈入交直流系统包括m回直流，一共n条母线时，每条母线的功率平衡方程表示为：When the multi-feed AC/DC system includes m circuits of DC and n buses in total, the power balance equation of each bus is expressed as:
Figure FDA0003232112820000011
Figure FDA0003232112820000011
 式中，△Pi、△Qi分别表示节点i注入的有功功率变化量和无功功率变化量，式(1)中的方程分别为节点i的有功功率方程和无功功率方程，PGi、QGi分别表示发电机注入节点i的有功出力和无功出力，PLi、QLi分别表示节点i的有功负荷和无功负荷，PDi表示节点i的直流功率，QDi表示直流换流器注入节点i的无功功率，Ui、Uj分别表示节点i、j的电压有效值，QSi表示动态无功补偿装置注入节点i的无功出力，Gij、Bij分别表示节点i、j之间的电导和电纳，θij表示节点i、j之间的电压相角差，对于所述多馈入交直流系统中的送端直流输电系统，式(1)中的直流有功功率在整流侧取负号，在逆变侧取正号；In the formula, ΔP i and ΔQ i represent the variation of active power and reactive power injected by node i, respectively, the equations in formula (1) are the active power equation and reactive power equation of node i, respectively, P Gi , Q Gi represent the active and reactive output of the generator injected into node i, respectively, P Li and Q Li represent the active load and reactive load of node i, respectively, P Di represents the DC power of node i, and Q Di represents the DC commutation is the reactive power injected by the device into node i, U i and U j represent the RMS voltages of nodes i and j respectively, Q Si represents the reactive power injected into node i by the dynamic reactive power compensation device, G ij and B ij represent node i respectively , the conductance and susceptance between j, θ ij represents the voltage phase angle difference between nodes i and j, and for the sending-end DC transmission system in the multi-feed AC-DC system, the DC active power in Eq. (1) The power takes a negative sign on the rectifier side and a positive sign on the inverter side; 当受端交流系统第k条交流母线发生故障的瞬间，对于式(1)，所述第k条交流母线对应的无功功率方程左侧出现不平衡量△Qk，而该母线的△Pk以及其他母线的功率方程仍然保持左侧项为零，设所述多馈入交直流系统中直流输电系统的第x条母线的电压有效值为Ux时，所述电压有效值的变化量为△Ux，受端交流系统的第k条母线的电压有效值为Uk，电压有效值的变化量为△Uk，系数△Ux/△Uk为直流输电系统第x条换流母线相对受端交流系统的第k条母线的电压耦合作用因子ADVCFxk，当1≤x≤m，m+1≤k≤n时，基于式(1)功率平衡方程，并考虑动态元件发电机组、直流系统、负荷特性以及动态无功补偿装置的影响，将式(1)两侧对Uk求导，生成式(2)，其表达式为：When the k-th AC bus of the receiving side AC system fails, for formula (1), the left side of the reactive power equation corresponding to the k-th AC bus appears unbalanced ΔQ k , while the ΔP k of the bus and the power equations of other busbars still keep the left-hand term as zero. When the RMS voltage of the xth busbar of the DC transmission system in the multi-feed AC/DC system is U x , the variation of the RMS voltage is △U x , the RMS voltage of the k-th bus of the receiving AC system is U k , the variation of the voltage RMS is △U k , and the coefficient △U x /△U k is the x-th commutation bus of the DC transmission system The voltage coupling factor ADVCF xk of the kth busbar of the AC system at the receiving end, when 1≤x≤m, m+1≤k≤n, is based on the power balance equation of formula (1), and considers the dynamic element generator set, For the influence of DC system, load characteristics and dynamic reactive power compensation device, the two sides of equation (1) are derived from U k to generate equation (2), and its expression is:
Figure FDA0003232112820000021
Figure FDA0003232112820000021
 令式(2)右端向量中的元素
Figure FDA0003232112820000022
生成式(3)，其表达式为：Elements in the vector on the right-hand side of (2)
Figure FDA0003232112820000022
Generating formula (3), its expression is:
Figure FDA0003232112820000031
Figure FDA0003232112820000031
 式(3)中，所述扩展雅克比矩阵为：In formula (3), the extended Jacobian matrix is:
Figure FDA0003232112820000032
Figure FDA0003232112820000032
 当i≠j时，雅克比矩阵中的元素Hij、Nij、Mij、Lij计算公式如下：When i≠j, the calculation formulas of elements Hi ij , N ij , M ij , and L ij in the Jacobian matrix are as follows:
Figure FDA0003232112820000041
Figure FDA0003232112820000041
 当i＝j时，雅克比矩阵中的元素Hii、Nii、Mii、Lii的计算公式如下：When i=j, the calculation formulas of elements H ii , N ii , M ii , and Li ii in the Jacobian matrix are as follows:
Figure FDA0003232112820000042
Figure FDA0003232112820000042
 由式(3)可知，雅克比矩阵为2n维方阵，待求解变量共2n-1个，其中母线k的无功功率方程对应的第2k行为冗余行，
Figure FDA0003232112820000043
为未知量，第2k列对应
Figure FDA0003232112820000044
删去第2k行，将第2k列移到式(3)左侧，得到方程式X，其表达式为：It can be seen from equation (3) that the Jacobian matrix is a 2n-dimensional square matrix, and there are 2n-1 variables to be solved. Among them, the 2kth row corresponding to the reactive power equation of the bus k is a redundant row,
Figure FDA0003232112820000043
is an unknown quantity, the 2kth column corresponds to
Figure FDA0003232112820000044
Delete the 2kth row and move the 2kth column to the left side of equation (3) to obtain equation X, whose expression is:
Figure FDA0003232112820000051
Figure FDA0003232112820000051
 采集所述多馈入交直流系统每条母线的网络参数值和动态元件参数值，并根据所述网络参数值和动态元件参数值求解扩展雅克比矩阵中每个元素的值，其中，所述网络参数值包括所述多馈入交直流系统中每条母线的有效电压值，以及母线与母线之间的电导、电纳和电压相角差，所述动态元件参数值是确定所述多馈入交直流系统中的动态元件的功率的参数值，所述动态元件包括发电机、负荷、直流换流器和动态无功补偿装置；Collect network parameter values and dynamic element parameter values of each bus of the multi-feed AC/DC system, and solve the value of each element in the extended Jacobian matrix according to the network parameter values and dynamic element parameter values, wherein the The network parameter value includes the effective voltage value of each bus in the multi-feed AC/DC system, and the conductance, susceptance and voltage phase angle difference between the bus and the bus, and the dynamic element parameter value is used to determine the multi-feed The parameter value of the power of the dynamic elements entering the AC-DC system, the dynamic elements include generators, loads, DC converters and dynamic reactive power compensation devices; 基于所述扩展雅克比矩阵中每个元素的值，采用稀疏技术对方程式X求解，确定每条直流输电系统换流母线对所述故障的受端交流母线的电压耦合作用因子。Based on the value of each element in the extended Jacobian matrix, equation X is solved by sparse technology to determine the voltage coupling action factor of each commutation bus in the DC transmission system to the AC bus at the receiving end of the fault. 2.根据权利要求1所述的方法，其特征在于，所述采集所述多馈入交直流系统每条母线的网络参数值和动态元件参数值，并根据所述网络参数值和动态元件参数值求解扩展雅克比矩阵中每个元素的值包括：2 . The method according to claim 1 , wherein the network parameter value and dynamic element parameter value of each bus of the multi-feed AC-DC system are collected, and the network parameter value and dynamic element parameter value are collected according to the network parameter value and the dynamic element parameter. 3 . Value Solving The value of each element in the extended Jacobian matrix consists of: 当i≠j时，根据采集的多馈入交直流系统每条母线的网络参数值和式(4)确定雅克比矩阵中的元素Hij、Nij、Mij、Lij的值；When i≠j, the values of elements Hij , Nij, Mij , Lij in the Jacobian matrix are determined according to the collected network parameter values of each busbar of the multi-feed AC-DC system and formula (4); 当i＝j时，根据采集的多馈入交直流系统每条母线的网络参数值、动态元件参数值和式(5)确定雅克比矩阵中的元素Hii、Nii、Mii、Lii的值，其中：When i=j, the elements H ii , N ii , M ii , Li ii in the Jacobian matrix are determined according to the collected network parameter values and dynamic element parameter values of each bus of the multi-feed AC-DC system and equation (5). value, where: 对于非发电机节点，式(5)中相关发电机的导数项
Figure FDA0003232112820000061
为零，在发电机节点，近似认为发电机在扰动瞬间次暂态电抗X"di后的电动势E"i保持恒定，发电机输出功率表达式为式(6)，其中θδi为发电机内电势E"i与机端电压UGi的角度之差，发电机功率对电压的导数项为式(7)：For non-generator nodes, the derivative term of the relevant generator in Eq. (5)
Figure FDA0003232112820000061
At the generator node, it is approximately considered that the electromotive force E" i of the generator after the sub-transient reactance X" di remains constant at the moment of disturbance, and the generator output power is expressed as formula (6), where θ δi is the internal value of the generator. The difference between the angle of the potential E" i and the terminal voltage U Gi , the derivative term of the generator power to the voltage is formula (7):
Figure FDA0003232112820000062
Figure FDA0003232112820000062
Figure FDA0003232112820000063
Figure FDA0003232112820000063
 对于非负荷节点，式(5)中负荷功率对电压的导数项
Figure FDA0003232112820000064
为零，另外，负荷功率仅与馈电点电压有效值有关，而与其角度无关，因此负荷功率对电压角度的导数
Figure FDA0003232112820000065
为零；For non-load nodes, the derivative term of load power to voltage in Eq. (5)
Figure FDA0003232112820000064
is zero, in addition, the load power is only related to the rms value of the feeding point voltage, not its angle, so the derivative of the load power to the voltage angle
Figure FDA0003232112820000065
zero; 对于负荷节点，当所述负荷为恒定功率负荷时，恒定电流负荷功率表达式及其对电压的导数为式(8)：For a load node, when the load is a constant power load, the power expression of the constant current load and its derivative to the voltage are formula (8):
Figure FDA0003232112820000071
Figure FDA0003232112820000071
 用有名值表示的直流换流器电压、电流方程为式(9)，由式(9)推导出的换流器功率方程式为式(10)，其中，方程(10)中关于换流母线的电压变量仅有电压有效值Ui，不含电压角度，式(9)和式(10)如下：The voltage and current equations of the DC converter represented by the famous values are equation (9), and the power equation of the converter derived from equation (9) is equation (10), where the equation (10) about the commutation bus is The voltage variable has only the RMS voltage U i , without the voltage angle. Equations (9) and (10) are as follows:
Figure FDA0003232112820000072
Figure FDA0003232112820000072
Figure FDA0003232112820000073
Figure FDA0003232112820000073
 式中，Ud表示直流线路的传输电压，nt表示六脉动换流器串联个数，kT表示换流器变压器变比，θd表示整流器的直流换相角或逆变器的熄弧角，Xc表示等效换相电抗，Id表示直流线路的传输电流，kγ表示换流变压器等效变比，
Figure FDA0003232112820000074
表示等效功率因数角，Ii表示直流注入交流系统的电流；In the formula, U d represents the transmission voltage of the DC line, n t represents the number of six-pulse converters in series, k T represents the converter transformer transformation ratio, and θ d represents the DC commutation angle of the rectifier or the arc extinguishing of the inverter. angle, X c represents the equivalent commutation reactance, I d represents the transmission current of the DC line, k γ represents the equivalent transformation ratio of the converter transformer,
Figure FDA0003232112820000074
represents the equivalent power factor angle, and I i represents the current injected by the DC into the AC system; 当直流输电系统整流侧采用定电流控制，逆变侧采用定熄弧角控制时，直流系统逆变侧换流器从交流系统吸取的功率对电压的导数为方程式(11)：When the rectifier side of the DC transmission system adopts constant current control and the inverter side adopts constant arc extinguishing angle control, the derivative of the power drawn by the converter on the inverter side of the DC system from the AC system to the voltage is equation (11):
Figure FDA0003232112820000081
Figure FDA0003232112820000081
 当动态无功补偿装置为静止无功补偿器时，所述静止无功补偿器采用被控母线的电压偏差作为输入信号，经比例放大环节控制补偿装置的等效电纳，忽略延迟环节，其输出等效电纳与电压偏差的关系为式(12)，对电压的导数为式(13)：When the dynamic reactive power compensation device is a static reactive power compensator, the static reactive power compensator uses the voltage deviation of the controlled bus as the input signal, and controls the equivalent susceptance of the compensation device through the proportional amplification link, ignoring the delay link, its The relationship between the output equivalent susceptance and the voltage deviation is equation (12), and the derivative to the voltage is equation (13): Bi＝-K△Ui＝-K(Ui-Ui0) (12)B i =-K△U i =-K(U i -U i0 ) (12)
Figure FDA0003232112820000082
Figure FDA0003232112820000082
 3.根据权利要求2所述的方法，其特征在于，所述当i＝j时，根据采集的多馈入交直流系统每条母线的网络参数值、动态元件参数值和式(5)确定雅克比矩阵中的元素Hii、Nii、Mii、Lii的值还包括：3 . The method according to claim 2 , wherein when i=j, it is determined according to the collected network parameter value, dynamic element parameter value and formula (5) of each busbar of the multi-feed AC-DC system. 4 . The values of elements H ii , N ii , M ii , and L ii in the Jacobian matrix also include: 对于负荷节点，当所述负荷为恒定阻抗负荷时，负荷功率的表达式及其对电压的导数为式(14)：For a load node, when the load is a constant impedance load, the expression of the load power and its derivative to the voltage are Equation (14):
Figure FDA0003232112820000091
Figure FDA0003232112820000091
 当直流输电系统整流侧采用定功率控制，逆变侧采用定熄弧角控制时，在换流母线上施加微小电压波动，应用控制策略，计算两侧换流站的功率变化，并应用差分替代偏微分；When the rectifier side of the DC transmission system adopts constant power control and the inverter side adopts constant arc extinguishing angle control, a small voltage fluctuation is applied to the converter bus, the control strategy is applied, the power changes of the converter stations on both sides are calculated, and the differential substitution is applied. partial differential; 当动态无功补偿装置为静止同步补偿器时，将被控母线的电压偏差作为输入信号，其稳态方程可表示为：When the dynamic reactive power compensation device is a static synchronous compensator, the voltage deviation of the controlled bus is used as the input signal, and its steady-state equation can be expressed as: △Ui＝UREF-Ui＝KDIS (15)△U i =U REF -U i =K D I S (15) IS＝△Ui/KD＝BS△Ui (16)I S =△U i /K D =B S △U i (16) 将静止同步补偿器的等效电纳表示为BS＝1/KD，当KD取零时，静止同步补偿器的控制节点为无差控制，但受到静止同步补偿器输出电流的限制，对电压的导数为式(17)：The equivalent susceptance of the static synchronous compensator is expressed as B S = 1/K D , when K D is zero, the control node of the static synchronous compensator is in differential control, but is limited by the output current of the static synchronous compensator, The derivative with respect to voltage is equation (17):
Figure FDA0003232112820000092
Figure FDA0003232112820000092
 4.一种基于扩展雅克比矩阵计算交直流系统电压耦合作用因子的系统，其特征在于，所述系统包括：4. A system for calculating a voltage coupling action factor of an AC-DC system based on an extended Jacobian matrix, wherein the system comprises: 方程确定单元，其用于确定多馈入交直流系统中每条母线的功率平衡方程，并在所述多馈入交直流系统的一条受端交流母线出现故障扰动时，将所述功率平衡方程两侧对发生故障的母线的电压求导，建立基于扩展雅克比矩阵的确定多馈入交直流系统中送端直流输电系统的每条换流母线对故障的受端交流母线的电压耦合作用因子的方程式X，其中：an equation determination unit, which is used to determine the power balance equation of each bus in the multi-feed AC-DC system, and when a receiving end AC bus of the multi-feed AC-DC system has a fault disturbance, the power balance equation The voltage of the faulted bus is derived on both sides, and the voltage coupling action factor of each commutation bus of the sending-end DC transmission system in the multi-infeed AC-DC transmission system based on the extended Jacobian matrix is established to the faulty receiving-end AC bus. Equation X of , where: 当所述多馈入交直流系统包括m回直流，一共n条母线时，每条母线的功率平衡方程表示为：When the multi-feed AC/DC system includes m circuits of DC and n buses in total, the power balance equation of each bus is expressed as:
Figure FDA0003232112820000101
Figure FDA0003232112820000101
 式中，△Pi、△Qi分别表示节点i注入的有功功率变化量和无功功率变化量，式(1)中的方程分别为节点i的有功功率方程和无功功率方程，PGi、QGi分别表示发电机注入节点i的有功出力和无功出力，PLi、QLi分别表示节点i的有功负荷和无功负荷，PDi表示节点i的直流功率，QDi表示直流换流器注入节点i的无功功率，Ui、Uj分别表示节点i、j的电压有效值，QSi表示动态无功补偿装置注入节点i的无功出力，Gij、Bij分别表示节点i、j之间的电导和电纳，θij表示节点i、j之间的电压相角差，对于所述多馈入交直流系统中的送端直流输电系统，式(1)中的直流有功功率在整流侧取负号，在逆变侧取正号；In the formula, ΔP i and ΔQ i represent the variation of active power and reactive power injected by node i, respectively, the equations in formula (1) are the active power equation and reactive power equation of node i, respectively, P Gi , Q Gi represent the active and reactive output of the generator injected into node i, respectively, P Li and Q Li represent the active load and reactive load of node i, respectively, P Di represents the DC power of node i, and Q Di represents the DC commutation is the reactive power injected by the device into node i, U i and U j represent the RMS voltages of nodes i and j respectively, Q Si represents the reactive power injected into node i by the dynamic reactive power compensation device, G ij and B ij represent node i respectively , the conductance and susceptance between j, θ ij represents the voltage phase angle difference between nodes i and j, and for the sending-end DC transmission system in the multi-feed AC-DC system, the DC active power in Eq. (1) The power takes a negative sign on the rectifier side and a positive sign on the inverter side; 当受端交流系统第k条交流母线发生故障的瞬间，对于式(1)，所述第k条交流母线对应的无功功率方程左侧出现不平衡量△Qk，而该母线的△Pk以及其他母线的功率方程仍然保持左侧项为零，设所述多馈入交直流系统中直流输电系统的第x条母线的电压有效值为Ux时，所述电压有效值的变化量为△Ux，受端交流系统的第k条母线的电压有效值为Uk，电压有效值的变化量为△Uk，系数△Ux/△Uk为直流输电系统第x条换流母线相对受端交流系统的第k条母线的电压耦合作用因子ADVCFxk，当1≤x≤m，m+1≤k≤n时，基于式(1)功率平衡方程，并考虑动态元件发电机组、直流系统、负荷特性以及动态无功补偿装置的影响，将式(1)两侧对Uk求导，生成式(2)，其表达式为：When the k-th AC bus of the receiving side AC system fails, for formula (1), the left side of the reactive power equation corresponding to the k-th AC bus appears unbalanced ΔQ k , while the ΔP k of the bus and the power equations of other busbars still keep the left-hand term as zero. When the RMS voltage of the xth busbar of the DC transmission system in the multi-feed AC/DC system is U x , the variation of the RMS voltage is △U x , the RMS voltage of the k-th bus of the receiving AC system is U k , the variation of the voltage RMS is △U k , and the coefficient △U x /△U k is the x-th commutation bus of the DC transmission system The voltage coupling factor ADVCF xk of the kth busbar of the AC system at the receiving end, when 1≤x≤m, m+1≤k≤n, is based on the power balance equation of formula (1), and considers the dynamic element generator set, For the influence of DC system, load characteristics and dynamic reactive power compensation device, the two sides of equation (1) are derived from U k to generate equation (2), and its expression is:
Figure FDA0003232112820000111
Figure FDA0003232112820000111
 令式(2)右端向量中的元素
Figure FDA0003232112820000112
生成式(3)，其表达式为：Elements in the vector on the right-hand side of (2)
Figure FDA0003232112820000112
Generating formula (3), its expression is:
Figure FDA0003232112820000113
Figure FDA0003232112820000113
 式(3)中，所述扩展雅克比矩阵为：In formula (3), the extended Jacobian matrix is:
Figure FDA0003232112820000121
Figure FDA0003232112820000121
 当i≠j时，雅克比矩阵中的元素Hij、Nij、Mij、Lij计算公式如下：When i≠j, the calculation formulas of elements Hi ij , N ij , M ij , and L ij in the Jacobian matrix are as follows:
Figure FDA0003232112820000122
Figure FDA0003232112820000122
 当i＝j时，雅克比矩阵中的元素Hii、Nii、Mii、Lii的计算公式如下：When i=j, the calculation formulas of elements H ii , N ii , M ii , and Li ii in the Jacobian matrix are as follows:
Figure FDA0003232112820000131
Figure FDA0003232112820000131
 由式(3)可知，雅克比矩阵为2n维方阵，待求解变量共2n-1个，其中母线k的无功功率方程对应的第2k行为冗余行，
Figure FDA0003232112820000132
为未知量，第2k列对应
Figure FDA0003232112820000133
删去第2k行，将第2k列移到式(3)左侧，得到方程式X，其表达式为：It can be seen from equation (3) that the Jacobian matrix is a 2n-dimensional square matrix, and there are 2n-1 variables to be solved. Among them, the 2kth row corresponding to the reactive power equation of the bus k is a redundant row,
Figure FDA0003232112820000132
is an unknown quantity, the 2kth column corresponds to
Figure FDA0003232112820000133
Delete the 2kth row and move the 2kth column to the left side of equation (3) to obtain equation X, whose expression is:
Figure FDA0003232112820000134
Figure FDA0003232112820000134
 元素计算单元，其用于采集所述多馈入交直流系统每条母线的网络参数值和动态元件参数值，并根据所述网络参数值和动态元件参数值求解扩展雅克比矩阵中每个元素的值，其中，所述网络参数值包括所述多馈入交直流系统中每条母线的有效电压值，以及母线与母线之间的电导、电纳和电压相角差，所述动态元件参数值是确定所述多馈入交直流系统中的动态元件的功率的参数值，所述动态元件包括发电机、负荷、直流换流器和动态无功补偿装置；an element calculation unit, which is used to collect the network parameter value and dynamic element parameter value of each bus of the multi-feed AC/DC system, and solve each element in the extended Jacobian matrix according to the network parameter value and the dynamic element parameter value , wherein the network parameter value includes the effective voltage value of each busbar in the multi-feed AC/DC system, and the conductance, susceptance and voltage phase angle difference between the busbar and the busbar, and the dynamic element parameter The value is a parameter value that determines the power of dynamic elements in the multi-feed AC-DC system, the dynamic elements including generators, loads, DC converters, and dynamic reactive power compensation devices; 因子确定单元，其用于基于所述扩展雅克比矩阵中每个元素的值，采用稀疏技术对方程式X求解，确定每条直流输电系统换流母线对所述故障的受端交流母线的电压耦合作用因子。A factor determination unit, which is used to solve equation X by using a sparse technique based on the value of each element in the extended Jacobian matrix, and determine the voltage coupling of each DC transmission system commutation bus to the faulty receiving-end AC bus Action factor. 5.根据权利要求4所述的系统，其特征在于，所述元素计算单元采集所述多馈入交直流系统每条母线的网络参数值和动态元件参数值，并根据所述网络参数值和动态元件参数值求解扩展雅克比矩阵中每个元素的值包括：5 . The system according to claim 4 , wherein the element calculation unit collects the network parameter value and dynamic element parameter value of each bus of the multi-feed AC/DC system, and calculates the network parameter value according to the network parameter value and the dynamic element parameter value. 6 . Dynamic Element Parameter Value Solving The value of each element in the extended Jacobian matrix consists of: 当i≠j时，根据采集的多馈入交直流系统每条母线的网络参数值和式(4)确定雅克比矩阵中的元素Hij、Nij、Mij、Lij的值；When i≠j, the values of elements Hij , Nij, Mij , Lij in the Jacobian matrix are determined according to the collected network parameter values of each busbar of the multi-feed AC-DC system and formula (4); 当i＝j时，根据采集的多馈入交直流系统每条母线的网络参数值、动态元件参数值和式(5)确定雅克比矩阵中的元素Hii、Nii、Mii、Lii的值，其中：When i=j, the elements H ii , N ii , M ii , Li ii in the Jacobian matrix are determined according to the collected network parameter values and dynamic element parameter values of each bus of the multi-feed AC-DC system and equation (5). value, where: 对于非发电机节点，式(5)中相关发电机的导数项
Figure FDA0003232112820000141
为零，在发电机节点，近似认为发电机在扰动瞬间次暂态电抗X"di后的电动势E"i保持恒定，发电机输出功率表达式为式(6)，其中θδi为发电机内电势E"i与机端电压UGi的角度之差，发电机功率对电压的导数项为式(7)：For non-generator nodes, the derivative term of the relevant generator in Eq. (5)
Figure FDA0003232112820000141
At the generator node, it is approximately considered that the electromotive force E" i of the generator after the sub-transient reactance X" di remains constant at the moment of disturbance, and the generator output power is expressed as formula (6), where θ δi is the internal value of the generator. The difference between the angle of the potential E" i and the terminal voltage U Gi , the derivative term of the generator power to the voltage is formula (7):
Figure FDA0003232112820000142
Figure FDA0003232112820000142
Figure FDA0003232112820000151
Figure FDA0003232112820000151
 对于非负荷节点，式(5)中负荷功率对电压的导数项
Figure FDA0003232112820000152
为零，另外，负荷功率仅与馈电点电压有效值有关，而与其角度无关，因此负荷功率对电压角度的导数
Figure FDA0003232112820000153
为零；For non-load nodes, the derivative term of load power to voltage in Eq. (5)
Figure FDA0003232112820000152
is zero, in addition, the load power is only related to the rms value of the feeding point voltage, not its angle, so the derivative of the load power to the voltage angle
Figure FDA0003232112820000153
zero; 对于负荷节点，当所述负荷为恒定功率负荷时，恒定电流负荷功率表达式及其对电压的导数为式(8)：For a load node, when the load is a constant power load, the power expression of the constant current load and its derivative to the voltage are formula (8):
Figure FDA0003232112820000154
Figure FDA0003232112820000154
 用有名值表示的直流换流器电压、电流方程为式(9)，由式(9)推导出的换流器功率方程式为式(10)，其中，方程(10)中关于换流母线的电压变量仅有电压有效值Ui，不含电压角度，式(9)和式(10)如下：The voltage and current equations of the DC converter represented by the famous values are equation (9), and the power equation of the converter derived from equation (9) is equation (10), where the equation (10) about the commutation bus is The voltage variable has only the RMS voltage U i , without the voltage angle. Equations (9) and (10) are as follows:
Figure FDA0003232112820000161
Figure FDA0003232112820000161
Figure FDA0003232112820000162
Figure FDA0003232112820000162
 式中，Ud表示直流线路的传输电压，nt表示六脉动换流器串联个数，kT表示换流器变压器变比，θd表示整流器的直流换相角或逆变器的熄弧角，Xc表示等效换相电抗，Id表示直流线路的传输电流，kγ表示换流变压器等效变比，
Figure FDA0003232112820000163
表示等效功率因数角，Ii表示直流注入交流系统的电流；In the formula, U d represents the transmission voltage of the DC line, n t represents the number of six-pulse converters in series, k T represents the converter transformer transformation ratio, and θ d represents the DC commutation angle of the rectifier or the arc extinguishing of the inverter. angle, X c represents the equivalent commutation reactance, I d represents the transmission current of the DC line, k γ represents the equivalent transformation ratio of the converter transformer,
Figure FDA0003232112820000163
represents the equivalent power factor angle, and I i represents the current injected by the DC into the AC system; 当直流输电系统整流侧采用定电流控制，逆变侧采用定熄弧角控制时，直流系统逆变侧换流器从交流系统吸取的功率对电压的导数为方程式(11)：When the rectifier side of the DC transmission system adopts constant current control and the inverter side adopts constant arc extinguishing angle control, the derivative of the power drawn by the converter on the inverter side of the DC system from the AC system to the voltage is equation (11):
Figure FDA0003232112820000164
Figure FDA0003232112820000164
 当动态无功补偿装置为静止无功补偿器时，所述静止无功补偿器采用被控母线的电压偏差作为输入信号，经比例放大环节控制补偿装置的等效电纳，忽略延迟环节，其输出等效电纳与电压偏差的关系为式(12)，对电压的导数为式(13)：When the dynamic reactive power compensation device is a static reactive power compensator, the static reactive power compensator uses the voltage deviation of the controlled bus as the input signal, and controls the equivalent susceptance of the compensation device through the proportional amplification link, ignoring the delay link, its The relationship between the output equivalent susceptance and the voltage deviation is equation (12), and the derivative to the voltage is equation (13): Bi＝-K△Ui＝-K(Ui-Ui0) (12)B i =-K△U i =-K(U i -U i0 ) (12)
Figure FDA0003232112820000171
Figure FDA0003232112820000171
 6.根据权利要求5所述的系统，其特征在于，所述元素计算单元在i＝j时，根据采集的多馈入交直流系统每条母线的网络参数值、动态元件参数值和式(5)确定雅克比矩阵中的元素Hii、Nii、Mii、Lii的值还包括：6. The system according to claim 5, characterized in that, when i=j, the element calculation unit collects network parameter values, dynamic element parameter values and formula ( 5) Determining the values of elements H ii , N ii , M ii , and L ii in the Jacobian matrix also includes: 对于负荷节点，当所述负荷为恒定阻抗负荷时，负荷功率的表达式及其对电压的导数为式(14)：For a load node, when the load is a constant impedance load, the expression of the load power and its derivative to the voltage are Equation (14):
Figure FDA0003232112820000172
Figure FDA0003232112820000172
 当直流输电系统整流侧采用定功率控制，逆变侧采用定熄弧角控制时，在换流母线上施加微小电压波动，应用控制策略，计算两侧换流站的功率变化，并应用差分替代偏微分；When the rectifier side of the DC transmission system adopts constant power control and the inverter side adopts constant arc extinguishing angle control, a small voltage fluctuation is applied to the converter bus, the control strategy is applied, the power changes of the converter stations on both sides are calculated, and the differential substitution is applied. partial differential; 当动态无功补偿装置为静止同步补偿器时，将被控母线的电压偏差作为输入信号，其稳态方程可表示为：When the dynamic reactive power compensation device is a static synchronous compensator, the voltage deviation of the controlled bus is used as the input signal, and its steady-state equation can be expressed as: △Ui＝UREF-Ui＝KDIS (15)△U i =U REF -U i =K D I S (15) IS＝△Ui/KD＝BS△Ui (16)I S =△U i /K D =B S △U i (16) 将静止同步补偿器的等效电纳表示为BS＝1/KD，当KD取零时，静止同步补偿器的控制节点为无差控制，但受到静止同步补偿器输出电流的限制，对电压的导数为式(17)：The equivalent susceptance of the static synchronous compensator is expressed as B S = 1/K D , when K D is zero, the control node of the static synchronous compensator is in differential control, but is limited by the output current of the static synchronous compensator, The derivative with respect to voltage is equation (17):
Figure FDA0003232112820000181
Figure FDA0003232112820000181
CN201811237027.7A 2018-10-23 2018-10-23 A method and system for calculating AC-DC system voltage coupling action factor based on extended Jacobian matrix Active CN109245097B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811237027.7A CN109245097B (en) 2018-10-23 2018-10-23 A method and system for calculating AC-DC system voltage coupling action factor based on extended Jacobian matrix

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811237027.7A CN109245097B (en) 2018-10-23 2018-10-23 A method and system for calculating AC-DC system voltage coupling action factor based on extended Jacobian matrix

Publications (2)

Publication Number Publication Date
CN109245097A CN109245097A (en) 2019-01-18
CN109245097B true CN109245097B (en) 2021-11-02

Family

ID=65081402

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811237027.7A Active CN109245097B (en) 2018-10-23 2018-10-23 A method and system for calculating AC-DC system voltage coupling action factor based on extended Jacobian matrix

Country Status (1)

Country Link
CN (1) CN109245097B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308200B (en) * 2020-04-09 2021-05-18 国网江苏省电力有限公司电力科学研究院 A method and system for analyzing harmonic conduction in a power grid
CN112350595B (en) * 2020-11-19 2021-09-03 中南大学 Analog impedance control method for inhibiting input unbalance influence of AC/DC matrix converter
CN115313396A (en) * 2021-12-21 2022-11-08 中国电力科学研究院有限公司 Method and device for acquiring voltage supporting capability of power system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103887810A (en) * 2014-03-12 2014-06-25 国家电网公司 Method for judging whether continuous phase changing of direct-current system fails on basis of short-circuit ratio dynamic tracking
CN106803153A (en) * 2017-01-23 2017-06-06 中国电力科学研究院 A kind of appraisal procedure and system of many feed-in ac and dc systemses commutation failure risks

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103887810A (en) * 2014-03-12 2014-06-25 国家电网公司 Method for judging whether continuous phase changing of direct-current system fails on basis of short-circuit ratio dynamic tracking
CN106803153A (en) * 2017-01-23 2017-06-06 中国电力科学研究院 A kind of appraisal procedure and system of many feed-in ac and dc systemses commutation failure risks

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Comprehensive Optimization for Power System with Multiple HVDC Infeed;Dong Yang;《2016 IEEE PES Asia-Pacific Power and Energy Conference》;20161230;全文 *
一种快速评估多馈入直流系统换相失败风险的方法;邵瑶;《中国电机工程学报》;20170620;第2-3页 *
周勤勇.多直流馈入受端电网直流受电规模研究.《中国博士学位论文全文数据库•工程科技Ⅱ辑》.2016, *
多直流馈入受端电网直流受电规模研究;周勤勇;《中国博士学位论文全文数据库•工程科技Ⅱ辑》;20160115;第71-76页 *
直流控制方式对多馈入交直流系统电压相互作用的影响;陈修宇;《电力系统自动化》;20120125;全文 *
西北750 kV 电网大容量新型FACTS设备应用研究;左玉玺;《电网技术》;20130830;附图7 *

Also Published As

Publication number Publication date
CN109245097A (en) 2019-01-18

Similar Documents

Publication Publication Date Title
CN102651551B (en) Method for controlling power factor of three-phase converter, method for controlling reactive power of three-phase converter, and controller of three-phase converter
CN109245097B (en) A method and system for calculating AC-DC system voltage coupling action factor based on extended Jacobian matrix
CN102088249A (en) Unbalanced voltage compensation method, unbalanced voltage compensator, three-phase converter control method, and controller of three-phase converter
CN109412191B (en) Phase locking method, device and equipment for high-voltage direct-current power transmission system
CN111740454B (en) A unified AC and DC voltage control method for hybrid microgrid based on bus interface converter
Hoseinnia et al. A control scheme for voltage unbalance compensation in an islanded microgrid
CN106451410B (en) A kind of DC dynamic voltage restorer and restoration method
Huber et al. Analysis, design, and evaluation of three-phase three-wire isolated ac-dc converter implemented with three single-phase converter modules
CN105576646A (en) Additional active signal-based balance control system and method for VSC-MTDC system
CN109256789B (en) A three-phase unbalance adjusting device and method for limiting current thereof
Wang et al. A novel hybrid directional comparison pilot protection scheme for the LCC-VSC hybrid HVDC transmission lines
Bassey et al. Modeling single-phase PQ inverter for unbalanced power dispatch in islanded microgrid
TW201834347A (en) Bidirectional power conversion device enables normal one within the first and second micro-grids to generate the inversion zero-sequence current for compensating the zero-sequence current through the control module generating the set of pulse width modulation signal
Mohammed et al. Performance evaluation of R-UPQC and L-UPQC based on a novel voltage detection algorithm
US20040052015A1 (en) Device and a method for voltage control in an electric transmission network
CN106533238A (en) Voltage compensation-based grid-connected inverter control method for marine electrical system
Wu et al. Frequency characteristic and impedance analysis on three-phase grid-connected inverters based on DDSRF-PLL
CN113039696A (en) Power conversion device
CN113517713B (en) Static voltage safety domain analysis method and device suitable for alternating current-direct current hybrid system
CN114268120A (en) MMC alternating-current side near-end asymmetric fault short-circuit current calculation method
CN115173441A (en) Voltage control method and device based on direct-current power distribution mode and storage medium
CN111969643A (en) Differential flat control method for MMC-HVDC (modular multilevel converter-high voltage direct current) supplying power to passive network under asymmetric fault
Yiauw et al. A novel three-phase active power filter
Hu et al. Virtual impedance based parallel operation of multi-converters in low voltage microgrids
CN114172172B (en) A non-isolated dynamic voltage restorer three-phase load imbalance control method

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