CN118916632A - Heterogeneous power supply power system correction participation analysis method based on impedance model - Google Patents
Heterogeneous power supply power system correction participation analysis method based on impedance model Download PDFInfo
- Publication number
- CN118916632A CN118916632A CN202410957283.2A CN202410957283A CN118916632A CN 118916632 A CN118916632 A CN 118916632A CN 202410957283 A CN202410957283 A CN 202410957283A CN 118916632 A CN118916632 A CN 118916632A
- Authority
- CN
- China
- Prior art keywords
- admittance
- matrix
- power supply
- sensitivity factor
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for AC mains or AC distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
-
- H02J2103/30—
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种基于阻抗模型的异构电源电力系统修正参与分析方法，涉及电力系统稳定控制领域。该方法生成异构电源电力系统的全系统节点阻抗矩阵和全系统节点导纳矩阵，从而获得元件的导纳灵敏度因子；根据元件的导纳灵敏度因子，确定参数灵敏度因子；考虑变压器支路变比对全系统节点导纳矩阵的影响，修正元件的导纳灵敏度因子和参数灵敏度因子；根据修正后的元件的导纳灵敏度因子和修正后的参数灵敏度因子，利用三层链式法则总体评价元件的参与程度。本发明能够准确地对异构电源电力系统电源或网络元件进行参与分析。
The present invention discloses a modified participation analysis method for a heterogeneous power system based on an impedance model, and relates to the field of power system stability control. The method generates a full-system node impedance matrix and a full-system node admittance matrix of a heterogeneous power system, thereby obtaining the admittance sensitivity factor of the element; according to the admittance sensitivity factor of the element, the parameter sensitivity factor is determined; considering the influence of the transformer branch ratio on the full-system node admittance matrix, the admittance sensitivity factor and the parameter sensitivity factor of the element are corrected; according to the corrected admittance sensitivity factor and the corrected parameter sensitivity factor of the element, the three-layer chain rule is used to evaluate the participation degree of the element as a whole. The present invention can accurately perform participation analysis on power supplies or network elements of a heterogeneous power system.
Description
Technical Field
The invention relates to the field of stable control of power systems, in particular to a heterogeneous power supply power system correction participation analysis method based on an impedance model.
Background
As the proportion of renewable energy sources in global energy structures continues to increase, the permeability of converter-INTERFACED GENERATION (CIG) interface power sources in electrical power systems increases significantly, which presents new stability challenges to the electrical power systems. In the last 10 years, high permeability converters have brought about multiple synchronous oscillation events to the actual grid. The power angle, frequency and voltage dynamics of the current transformer occur over a wider and faster time scale compared to synchronous machines of conventional power systems, which may further bring about new electromagnetic modes. When the permeability of the converter reaches 100%, the dynamic characteristics of the power system will change radically.
Current converters that are connected to ac power systems are mainly classified into grid-form (GFM) converters and grid-form (GFL) converters, which have similarities with synchronous machines in control modes, but have some important differences, which make a new modeling framework necessary for analyzing the stability of heterogeneous power supply power systems. Impedance is an attribute capable of reflecting dynamic characteristics of a power supply, and related scholars propose a method for judging the stability of small signals of a multi-equipment power system based on a generalized Nyquist stability criterion. In addition, the impedance network transformation method of the integrated subsystem can also realize the determination of stability, and can simplify the complexity of generating a network matrix for stability analysis. There is also a relevant scholars proposing a stability criterion for positive mode damping, but this only applies to weak damped oscillation modes. The method for embedding the dynamics of the coordinate transformation into the impedance model realizes the integral impedance modeling of the alternating current power grid. For a multi-converter system, these methods can only determine whether the system is stable, and cannot locate the element with the greatest influence on the system instability.
The impedance model participation analysis method provides an admittance sensitivity factor and a parameter sensitivity factor, and the admittance sensitivity factor can fairly participate in analysis on the admittance of the power supply and the network element, so that the impedance model participation analysis of the whole system is realized. In addition, the related literature further applies the impedance model participation analysis method to an AC/DC hybrid system, and illustrates the equivalence of the method in the AC/DC hybrid system for direct and indirect analysis. However, the impedance model has a large error in some scene applications, and the influence factors causing the error need to be discussed and the theoretical calculation process of the method needs to be corrected. Therefore, the impedance model-based heterogeneous power supply power system correction participation analysis method based on the impedance model is realized.
The prior researches mainly guide how to position the element with the largest participation degree in the system mode from the participation analysis result of the admittance sensitivity of the power supply or the network element, and the prior researches evaluate the misjudgment phenomenon of the total participation degree of the element based on the maximum value obtained by the Cauchy inequality, so that the mispositioning of important elements in the system can be caused. In addition, the mode analysis of the network element by the existing impedance model mode analysis method is only aimed at the circuit element, and the transformer branch can generate larger numerical error on the mode analysis result of the impedance model due to the transformation ratio, so that the stability of a small signal of the system can be weakened due to the adopted error parameter adjustment measures.
Disclosure of Invention
The invention aims to provide a heterogeneous power supply power system correction participation analysis method based on an impedance model, which can accurately perform participation analysis on a heterogeneous power supply power system power supply or a network element.
In order to achieve the above object, the present invention provides the following solutions:
An impedance model-based heterogeneous power supply power system correction participation analysis method comprises the following steps: generating a full-system node impedance matrix and a full-system node admittance matrix of the heterogeneous power supply power system; obtaining an admittance sensitivity factor of the element according to the full-system node impedance matrix and the full-system node admittance matrix; determining parameter sensitivity factors according to admittance sensitivity factors of the elements, and correcting the admittance sensitivity factors and the parameter sensitivity factors of the elements by considering the influence of transformer branch transformation ratios on an admittance matrix of all system nodes; obtaining a modal prediction variable quantity caused by parameter change according to the corrected parameter sensitivity factor; determining errors of a modal prediction variable quantity and a modal actual variable quantity caused by parameter changes; and according to the admittance sensitivity factor and the corrected parameter sensitivity factor of the element, the participation degree of the element is generally evaluated by using a three-layer chain rule.
Optionally, generating a system-wide node impedance matrix and a system-wide node admittance matrix of the heterogeneous power supply power system specifically includes:
the power supply node admittance matrix of the heterogeneous power supply power system is established as follows: wherein Y G represents the power supply node admittance matrix, Y G1 represents the admittance of the power supply connected with the bus 1, Y Gi represents the admittance of the power supply connected with the bus i, and Y Gn represents the admittance of the power supply connected with the bus n;
The network node admittance matrix of the heterogeneous power supply power system is established as follows: Wherein Y N represents a network node admittance matrix, Y N11 represents a self-admittance element of bus bar 1, Y N1j represents a transadmittance element between bus bar 1 and bus bar j, Y N1n represents a transadmittance element between bus bar 1 and bus bar n, Y Ni1 represents a transadmittance element between bus bar i and bus bar 1, Y Nij represents a transadmittance element between bus bar i and bus bar j, Y Nin represents a transadmittance element between bus bar i and bus bar n, Y Nn1 represents a transadmittance element between bus bar n and bus bar 1, Y Nnj represents a transadmittance element between bus bar n and bus bar j, and Y Nnn represents a self-admittance element of bus bar i;
Taking a disturbance value of injection current at each bus in a network as input, taking a voltage change value at each bus in the network as output, and obtaining a full-system node impedance matrix of the heterogeneous power supply power system by utilizing a reciprocal relation of a node admittance matrix and a node impedance matrix, wherein the full-system node impedance matrix comprises the following components: z= (I+Z NYG)-1ZN＝(YG+YN)-1＝Y-1) the total system node admittance matrix is Y=Y G+YN, wherein Z N is the node impedance matrix of the network, I is the identity matrix, Z represents the total system node impedance matrix taking the power supply impedance characteristic into account, and Y represents the total system node admittance matrix taking the power supply impedance characteristic into account.
Optionally, taking a disturbance value of injection current at each bus in the network as input, taking a voltage variation value at each bus in the network as output, and obtaining a node impedance matrix of the whole system of the heterogeneous power supply power system by utilizing a reciprocal relation of the node admittance matrix and the node impedance matrix, wherein the node impedance matrix is as follows: z= (i+z NYG)-1ZN＝(YG+YN)-1＝Y-1, full system node admittance matrix is y=y G+YN, specifically including:
Taking the disturbance value of the injected current at each bus in the network as an input, the input is: Δu= [ Δi 1…ΔIi…ΔIn]T; wherein Deltau is input, deltaI 1 is a disturbance value of injection current at the bus 1, deltaI i is a disturbance value of injection current at the bus I, deltaI n is a disturbance value of injection current at the bus n;
taking the voltage change value at each bus of the network as output, the output is: Δy= [ Δu 1…ΔUi…ΔUn]T; wherein Δy is the output, Δu 1 is the voltage change value at bus 1, Δu i is the voltage change value at bus i, and Δu n is the voltage change value at bus n;
According to the input and the output, an input-output transfer function matrix is obtained by utilizing a closed-loop formula, wherein the input-output transfer function matrix is delta y= (I+Z NYG)-1ZN delta u;
according to the input and output transfer function matrix, determining a node impedance matrix of the whole system of the heterogeneous power supply power system as follows: z= (i+z NYG)-1ZN＝(YG+YN)-1＝Y-1, full system node admittance matrix is y=y G+YN.
Optionally, the admittance sensitivity factor of the element is obtained according to the system-wide node impedance matrix and the system-wide node admittance matrix, which specifically comprises:
according to the total system node impedance matrix and the total system node admittance matrix, a partial derivative formula for taking partial derivatives of admittances of power supplies or network elements in the system by the mode is obtained as follows: wherein, Representing the partial derivative of the modal lambda to the admittance Y of the power supply or network element, Y ij representing the transadmittance between busbar i and busbar j in the node admittance matrix Y taking into account the power supply impedance characteristics, tr () representing the trace of the matrix, res λ Z representing the taking of the remainder of the overall system node impedance matrix Z at modal lambda;
Expanding the partial conductance formula to obtain an expansion formula of the partial conductance formula, wherein the expansion formula comprises the following steps: Wherein Z ii represents the self-admittance element of bus bar i, Z jj represents the self-admittance element of bus bar j, Z ji is the impedance between bus bar j and bus bar i, and Z ij is the impedance between bus bar i and bus bar j; res λZii represents taking the remainder at modality λ for Z ii, res λZjj represents taking the remainder at modality λ for Z jj, res λZij represents taking the remainder at modality λ for Z ij, res λZji represents taking the remainder at modality λ for Z ji;
according to the expansion of the partial conductivity formula, the admittance sensitivity factors of the obtained element are as follows: where s λ,y is the admittance sensitivity factor of the element, which represents the conjugate transpose.
Optionally, the parameter sensitivity factor is calculated according to the following formula: Wherein s λ,ρ is a parameter sensitivity factor, s λ,y is an admittance sensitivity factor of the element, For the sensitivity information of element admittance to parameters, < > represents the Frobenius inner product.
Optionally, considering the influence of the transformer branch transformation ratio on the admittance matrix of the nodes of the whole system, the admittance sensitivity factor and the parameter sensitivity factor of the correction element specifically comprise:
Considering the influence of the transformer branch transformation ratio on the partial conductance formula, the expansion of the correction partial conductance formula is as follows: wherein k is the transformer transformation ratio;
Correcting the admittance sensitivity factor of the element according to the expansion of the corrected partial conductance formula;
according to the admittance sensitivity factor of the element after correction, the formula is utilized Correcting the parameter sensitivity factor; where s 'λ,ρ is the modified parameter sensitivity factor, s' λ,y is the modified admittance sensitivity factor of the element,Sensitivity information for the element admittance to the parameters,Δρ is the small perturbation of the parameter ρ and y Δρ is the result of the change in admittance y of the power supply or network element after the change in the parameter ρ.
Optionally, considering the influence of the transformer branch transformation ratio on the admittance matrix of the all-system node, correcting the admittance sensitivity factor and the parameter sensitivity factor of the element, and then further including: according to the corrected parameter sensitivity factor and the formula delta lambda pr＝s′λ,ρ rho 5%, obtaining the modal prediction variable quantity caused by 5% of parameter change; wherein Δλ pr is the modal prediction variation caused by 5% change in parameter ρ, and s' λ,ρ is the modified parameter sensitivity factor; and determining errors of the modal prediction variation and the modal actual variation caused by the parameter variation of 5%.
Optionally, the error determination formula is: Wherein Error is the Error of the predicted mode change amount and the actual mode change amount caused by parameter change, deltalambda pr is the predicted mode change amount caused by 5% of parameter ρ change, deltalambda is the actual mode change amount after parameter change.
Optionally, the participation degree of the element is generally evaluated by using a three-layer chain rule according to the admittance sensitivity factor and the parameter sensitivity factor of the element after correction, and specifically includes:
the first layer defining the three-layer chain law is: Wherein Δλ max is the maximum variation of the mode corresponding to the admittance variation of the element, and σ 2 and ω 2 are the real part and the imaginary part of the second layer of the three-layer chain rule, respectively;
According to the admittance sensitivity factor of the corrected element, defining a second layer of the three-layer chain rule as follows: σ 2+jω2＝<s′λ,y, y >; wherein s' λ,y is the admittance sensitivity factor of the modified element, y is the admittance of the power supply or network element, and < represents the Frobenius inner product;
according to the corrected parameter sensitivity factor, a third layer of the three-layer chain rule is defined as follows: wherein, For the modified parameter sensitivity factor, σ 3 and ω 3 are the real and imaginary parts, respectively, of the third layer of the three-layer chain law;
the degree of participation of the elements was generally evaluated using the first, second and third layers of the three-layer chain law.
Optionally, the participation degree of the element is generally evaluated by using a three-layer chain rule according to the admittance sensitivity factor and the parameter sensitivity factor after correction, and then the method further comprises: and according to the corrected parameter sensitivity factor, aiming at enhancing the disturbance stability of the heterogeneous power supply power system, and determining the parameter adjustment direction.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
The invention discloses a heterogeneous power supply power system correction participation analysis method based on an impedance model, which considers the influence of transformer branch transformation ratio on a node admittance matrix of a whole system, corrects admittance sensitivity factors and parameter sensitivity factors of elements, further evaluates the participation degree of the elements by utilizing a three-layer chain rule overall according to the admittance sensitivity factors and the corrected parameter sensitivity factors of the elements, effectively solves the problem of unstable element judgment errors, and can accurately perform participation analysis on power supplies or network elements of the heterogeneous power supply power system.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed 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 other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a flow chart of a heterogeneous power supply power system correction participation analysis method based on an impedance model according to an embodiment of the present invention.
Fig. 2 is a schematic structural diagram of a heterogeneous power supply power system according to an embodiment of the present invention.
Fig. 3 is a schematic diagram of an ideal transformer branch according to an embodiment of the present invention.
Fig. 4 is a schematic diagram of a retrofit IEEE14 bus system provided by an embodiment of the present invention.
Fig. 5 is a schematic diagram of modal distribution of a modified IEEE14 bus system according to an embodiment of the present invention.
Fig. 6 is a schematic diagram of calculation results of a first layer and a second layer of a three-layer chain rule of a mode 1 according to an embodiment of the present invention.
Fig. 7 is a schematic diagram of calculation results of a first layer and a second layer of a three-layer chain rule of a mode 2 according to an embodiment of the present invention.
Fig. 8 is a schematic diagram of calculation results of the first layer and the second layer of the three-layer chain rule of the mode 3 according to the embodiment of the present invention.
Fig. 9 is a schematic diagram of calculation results of the first layer and the second layer of the three-layer chain rule of the modality 4 according to the embodiment of the present invention.
Fig. 10 is a schematic diagram of a first layer correction result of a three-layer chain rule of mode 2, mode 3 and mode 4 according to an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention aims to solve the technical problems that the impedance model participation analysis method is misjudged on the numerical error of transformer element participation analysis and the participation degree of evaluation element, and specifically comprises the following steps: according to the error influence of transformer transformation ratio and the erroneous judgment result of the participation degree of the overall evaluation element of the impedance model participation analysis method, the impedance model participation analysis method needs to be further corrected. The invention aims to provide an impedance model-based heterogeneous power supply power system correction participation analysis method based on an impedance model, which accurately performs participation analysis on a heterogeneous power supply power system power supply or a network element.
In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.
As shown in fig. 1, a heterogeneous power supply power system correction participation analysis method based on an impedance model in the present embodiment includes the following steps 1 to 7.
Step 1: and generating a full-system node impedance matrix and a full-system node admittance matrix of the heterogeneous power supply power system.
And (one) establishing an impedance model of the power supply and the network element.
And establishing a small signal model of the heterogeneous power supply power system, and respectively carrying out impedance modeling on a power supply and a power network of the heterogeneous power supply power system.
The small signal model of the power supply can be obtained through deduction of a state equation and measurement. Notably, in multi-machine power system impedance modeling, individual power supply impedance models derived from state equations need to ensure coordinate system alignment of all power supplies. Y Gi denotes the admittance of the power supply to which busbar i is connected:
Where Z Gi represents the impedance model of the power supply to which bus i is connected, AndRespectively the corresponding elements of the dq rotational coordinate system. If there is no power at bus i, there is Y Gi =0.
The power admittance matrix of the system can be obtained by each power impedance model:
where Y G represents the admittance matrix of the system power supply. Y G1 represents the admittance of the power supply to which bus 1 is connected, Y Gi represents the admittance of the power supply to which bus i is connected, and Y Gn represents the admittance of the power supply to which bus n is connected.
Establishing a network node admittance matrix Y N for the power network:
Where Y Nij represents the transadmittance element between bus bar i and j, and Y Nii represents the self-admittance element of bus bar i. Y N11 represents the self-admittance element of bus bar 1, Y N1j represents the transadmittance element between bus bar 1 and bus bar j, Y N1n represents the transadmittance element between bus bar 1 and bus bar n, Y Ni1 represents the transadmittance element between bus bar i and bus bar 1, Y Nij represents the transadmittance element between bus bar i and bus bar j, Y Nin represents the transadmittance element between bus bar i and bus bar n, Y Nn1 represents the transadmittance element between bus bar n and bus bar 1, Y Nnj represents the transadmittance element between bus bar n and bus bar j, and Y Nnn represents the self-admittance element of bus bar i.
And (II) generating a full system node impedance matrix.
The structure of the heterogeneous power supply power system is shown in fig. 2. Consider taking as input the disturbance value of the injected current at each busbar in the network:
Δu＝[ΔI1…ΔIi…ΔIn]T (4)
where ΔI i is the injection current disturbance value of bus I.
Take as output the voltage change value deltay at each bus of the network:
Δy＝[ΔU1…ΔUi…ΔUn]T (5)
where ΔU i is the voltage change value at bus i.
The input/output transfer function matrix can be obtained by using a closed-loop formula as follows:
Δy＝(I+ZNYG)-1ZNΔu (6)
Wherein the method comprises the steps of Is the network node impedance matrix and I is the identity matrix.
The transfer function matrix of the whole system can be obtained by the method:
Z＝(I+ZNYG)-1ZN＝(YG+YN)-1＝Y-1 (7)
Where Z and Y are a system-wide node impedance matrix and a system-wide node admittance matrix that take into account the power supply impedance characteristics. The total system node admittance matrix is: y=y G+YN.
Step 2: and obtaining the admittance sensitivity factor of the element according to the full-system node impedance matrix and the full-system node admittance matrix.
The modal λ takes the partial derivative of the admittance y of a certain power supply or network element in the system to obtain:
Wherein the method comprises the steps of Representing the partial derivative of the modal lambda with respect to the admittance Y of the power supply or network element, Y ij represents the transadmittance between bus i and bus j in the full system node admittance matrix Y. tr () represents the trace of the matrix. Res λ Z represents taking the remainder of matrix Z at modality λ, which can be identified from the spectrum by a rational approximation.
The above expansion can be obtained:
Wherein Z ii represents the self-admittance element of bus bar i, Z jj represents the self-admittance element of bus bar j, Z ji is the impedance between bus bar j and bus bar i, and Z ij is the impedance between bus bar i and bus bar j. Res λZii represents taking the remainder at modality λ for Z ii, res λZjj represents taking the remainder at modality λ for Z jj, res λZij represents taking the remainder at modality λ for Z ij, res λZji represents taking the remainder at modality λ for Z ji.
The admittance sensitivity factor defining y is:
Where s λ,y is the admittance sensitivity factor of the element, which represents the conjugate transpose.
Then there are:
Δλ＝<sλ,y,Δy> (11)
Where < > represents the Frobenius inner product, which describes the amount of change in the mode corresponding to the change in admittance of a certain element in the system.
Step 3: and determining a parameter sensitivity factor according to the admittance sensitivity factor of the element, and correcting the admittance sensitivity factor and the parameter sensitivity factor of the element by considering the influence of the transformer branch transformation ratio on the admittance matrix of the nodes of the whole system.
If the variations caused by the parameter ρ in admittance y are considered, there is a parameter sensitivity factor:
Wherein s λ,ρ is a parameter sensitivity factor, s λ,y is an admittance sensitivity factor of the element, Sensitivity information for the element admittance to the parameter. The parameter ρ may be any parameter describing the admittance y dynamic equation. The modal change amount caused by the change of the parameter ρ is:
Δλ＝sλ,ρ·Δρ (13)
The transformer branch transformation ratio can bring larger error to the impedance model modal analysis method. The calculation process needs to be modified. Consider an ideal transformer branch between bus i and bus j, with transformation ratio k on bus i side. The ideal transformer branch is shown in fig. 3.
The contribution of the ideal transformer branch to the admittance matrix of the all-system node is as follows:
the transformer transformation ratio k will have an effect on the resulting equation (8) and therefore the result of the series branch in equation (9) needs to be corrected:
that is, the expansion of the correction partial conductance formula is:
Correcting the admittance sensitivity factor of the element according to the expansion of the corrected partial conductance formula;
according to the admittance sensitivity factor of the element after correction, the formula is utilized Correcting the parameter sensitivity factor; where s 'λ,ρ is the modified parameter sensitivity factor, s' λ,y is the modified admittance sensitivity factor of the element,Sensitivity information for the element admittance to the parameters,Δρ is the small perturbation of the parameter ρ and y Δρ is the result of the change in admittance y of the power supply or network element after the change in the parameter ρ.
And 3, error analysis and error correction can be performed, and specifically, the method comprises the following steps:
And obtaining the modal prediction variable quantity caused by parameter change according to the corrected parameter sensitivity factor.
The predicted change in mode due to the 5% change in the parameter ρ is:
Δλpr＝s′λ,ρ·ρ·5％ (16)
Wherein Δλ pr is the modal prediction variation caused by 5% change in parameter ρ, and s' λ,ρ is the modified parameter sensitivity factor.
And determining errors of the modal prediction variable quantity and the modal actual variable quantity caused by parameter changes.
After the parameter is changed, the actual change value of the system is delta lambda, and the Error of the predicted value and the actual value is:
Step 4: and according to the admittance sensitivity factor and the corrected parameter sensitivity factor of the element, the participation degree of the element is generally evaluated by using a three-layer chain rule.
The three-layer chain rule is used for carrying out modal analysis of different transparency on each element of the system. Equation (11) is scaled using the cauchy inequality:
Δλmax＝||sλ,y||·||Δy||＝ε||sλ,y||·||y|| (18)
wherein|| representation of the Frobenius norm, epsilon represents the change proportion of y, and the value is generally smaller.
Then, λ,y ·y is defined as the first layer of the impedance model modality analysis method to measure the possible overall participation of y in the model λ.
The second layer and the third layer are:
σ2+jω2＝<sλ,y,y> (19)
Where σ 2 and ω 2 are the real and imaginary parts of the second layer, respectively, and σ 3 and ω 3 are the real and imaginary parts of the second layer, respectively.
The second and third layers describe the degree of participation in the modality from the element and parameter perspective, respectively. The first layer of the three-layer analysis method obtains what elements of the system may have great participation, the second layer obtains the influence of the elements on the modal damping (real part) and natural oscillation frequency (imaginary part), and the third layer requires sensitivity information of the elements on parameters compared with the second layer, so that the influence of specific parameters on the mode can be analyzed.
From the point of view of information acquisition, the information used by the first layer and the second layer is the same. The first layer, thus obtained by scaling the cauchy-inequality, is not suitable for measuring the degree of participation of all elements in its entirety. The first layer redefined here to measure the participation of all elements is:
Where Δλ max is the maximum amount of change in the mode corresponding to the change in admittance of the element.
According to the admittance sensitivity factor of the corrected element, defining a second layer of the three-layer chain rule as follows: σ 2+jω2＝<s′λ,y, y >; where s' λ,y is the admittance sensitivity factor of the modified element.
According to the corrected parameter sensitivity factor, a third layer of the three-layer chain rule is defined as follows: wherein, Is the modified parameter sensitivity factor.
The degree of participation of the elements was generally evaluated using the first, second and third layers of the three-layer chain law.
And after the step 4, key element parameters can be adjusted to stabilize the heterogeneous power supply power system:
Sensitivity information of element admittance to its internal parameters
Where Δρ is the small perturbation of the element parameter ρ and y Δρ is the result of the change in admittance y after the change in the element parameter ρ. Substituting equation (22) into equation (20) allows the overall evaluation of the degree of participation of the element.
Substituting equation (22) into equation (12) propagates the parameter perturbation to the impedance perturbation by calculating the parameter sensitivity factor of equation (12). Thereby guiding the selection of the adjustment direction of the internal parameters of the element and realizing the enhancement of the small disturbance stability of the heterogeneous power supply power system. The specific adjustment mode is as follows: and according to the corrected parameter sensitivity factor, aiming at enhancing the disturbance stability of the heterogeneous power supply power system, and determining the parameter adjustment direction.
The invention provides an impedance model-based heterogeneous power supply power system correction participation analysis method based on an impedance model, which comprises the steps of establishing an impedance model comprising all power supplies and network elements of a system, generating a node impedance matrix of the whole system according to the impedance models of the power supplies and the network elements, generating an admittance sensitivity matrix of the elements according to the node impedance matrix of the whole system, and obtaining a correction admittance sensitivity factor to participate in analysis on the system elements.
The method solves the problem that the unstable element judgment is wrong due to calculation errors in the impedance model participation analysis method in the actual heterogeneous power supply power system. The error influence of the transformer transformation ratio on the impedance model modal analysis method is considered, and the misjudgment phenomenon of the participation degree of the evaluation element is specifically analyzed through a plurality of modal results, and a correction result is given. The device can rapidly and conveniently position the element with the highest unstable participation degree to the system in the heterogeneous power supply power system.
Thus, the heterogeneous power supply power system correction participation analysis method based on the impedance model is designed.
The method of the present invention is described in detail below in conjunction with specific examples:
Table 1 gives the calculated results and errors without and with correction. It is pointed out here that the latter calculated values for the modes are all considered frequency values divided by 2pi, the modes lambda 1 and lambda 2 being-4.0935 +30.0633i and-2.8589 +7.3169i, respectively.
Table 1 error results of transformer transformation ratio correction versus model analysis
| λ,ρ | λ1,L56 | λ2,L56 |
| Uncorrected Δλ pr | -0.0741+0.2415i | 0.0084-0.1039i |
| Correction of Deltalambda pr | -0.0612+0.2129i | 0.0084-0.0901i |
| Actual Deltalambda | -0.0621+0.2059i | 0.0084-0.0871i |
| Unmodified Error | 14.854％ | 16.146％ |
| Correction Error | 3.165％ | 3.251％ |
From the results of table 1, it can be seen that, for the modes λ 1 and λ 2, the error calculated by the impedance model modal analysis method is only about 3% after the transformer transformation ratio correction is considered, and the error result without the transformer transformation ratio correction is about 15%. The error result is almost 5 times.
It is pointed out here that the line branch is only a special case of a transformation ratio of 1. The calculation results in table 1 are calculated when the original value of the transformation ratio is 0.932, and the further the transformation ratio deviates from 1, the larger the calculated error becomes. If the star connection of the transformer is considered, the transformation ratio is a complex number with an angular offset of 30 degrees. The error of the calculation result will be further increased if the transformer transformation ratio is not considered. Therefore, when the impedance model modal analysis method is actually used, it is necessary to consider the correction of the transformer transformation ratio.
The modal analysis of the series line parameters from previous studies was limited to only one 3-node passive system. The first layer of the chain law is scaled by the cauchy inequality and is considered as an index that can describe the degree of participation of the different elements in general. This is a significant error in the actual system analysis, since the cauchy inequality yields only one possible value, which needs to be corrected.
The modal distribution of the low frequency portion of the system shown in fig. 4 is shown in fig. 5. Consider the 4 modes marked red as in fig. 5, for which the first and second layers of the three-layer chain law are calculated. The calculation results of modality 1, modality 2, modality 3, and modality 4 are shown in fig. 6, fig. 7, fig. 8, and fig. 9. Parts (a), (b) and (c) in fig. 6, 7, 8 and 9 are the first layer of the three-layer chain law, the second real part of the three-layer chain law (representing damping) and the second imaginary part of the three-layer chain law (representing natural frequency), respectively.
Taking the mode 1 of fig. 6 as an example, the part (a) of fig. 6 represents the first layer of the three-layer chain rule, to measure the participation degree of all the elements, and the darker the color of the color represents the greater the participation degree of the element in the mode 1. The diagonal elements represent the power sources connected with the buses, and the non-diagonal elements represent the participation degree of the branches connected between the two buses in series. The white portion of the figure represents that no actual line or power source is present. Part (a) of fig. 6 shows that the mesh-type converter A6 connected to the busbar 6 in the mode 1 has the highest participation, and the second is that the participation of the branches 5-6 is also larger. The greater the degree of participation the greater the impact of the disturbance of the element on modality 1.
Part (b) of fig. 6 shows the change in the real part of the mode 1 when the admittances of the elements increase. The green color represents a movement away from the imaginary axis, which increases the modal damping. The red color represents a movement towards the direction close to the imaginary axis, which makes the modal damping smaller. As can be seen from fig. 6, if one considers increasing the number of parallels of the grid-connected converters at the bus bar 6 or decreasing their reference power, this will increase their admittance and thus the damping of this mode. If the branches 5-6 increase their admittance by connecting the same-scale branches in parallel, this will result in a reduction of the damping of the mode 1. It is noted that it is emphasized here that the direction of admittance change coincides with the direction of the original admittance, which brings about the damping change effect shown in fig. 6. When a certain parameter in an element changes, the mode does not necessarily change according to the direction in the figure, and this relates to a third layer of mode analysis. Part (c) of fig. 6 represents the change in the modal frequency as the admittance of each element increases.
As can be seen from fig. 6 and 7, when the admittances of the branches 5-6 or of the grid-connected converter A6 change, this brings about a diametrically opposite change in the modes 1 and 2. This means that adjusting the parameters of the various elements of the system does not increase the damping of all modes, i.e. the double-sided nature of the damping. If the system parameters are adjusted, adjustments need to be made taking into account the direction in which the damping of mode 2 increases. Since the mode 2 itself dampens less, nearer to the imaginary axis, it is easier to excite oscillations.
The maximum possible values for the four modal amplitude variations shown in part (a) of each of fig. 6 to 9 are approximately 120, 600, 5 and 120, respectively. The maximum amplitude values calculated here taking the maximum values of the damping part (b) and the frequency part (c) in fig. 6 to 9 are approximately 85, 28, 1.2 and 0.5, respectively. The difference in value from the first layer results calculated by the previous method is very large, and the modal 4 variation amplitude calculated by the previous method is about 240 times the actual variation.
Part (a) of fig. 8 and 9 shows that the participation degree of the mesh type converter A6 in the mode 3 and the mode 4 is the greatest. However, both parts (b) and (c) of fig. 8 show that the synchronous machine A8 and the branches 7-8 are involved to a greater extent than the grid-connected converter A6. Similarly, parts (b) and (c) of fig. 9 show that the participation degree of the synchronous machine A2, the synchronous machine A3, and the synchronous machine A8 is also high.
Fig. 10 shows the calculated first layer after modality 2, modality 3 and modality 4 correction. Comparing fig. 7-9, it can be seen that the results of the modified calculations more truly reflect the overall participation of each element in the system modality.
The scheme of the invention has the beneficial effects that: the invention provides a correction participation analysis method of a heterogeneous power supply power system based on an impedance model. The method effectively solves the problems that the impedance model participates in larger calculation errors of the analysis method when the large-scale heterogeneous power supply power system is applied and unstable elements are caused to judge wrongly. The error correction considering the transformer transformation ratio and the misjudgment reason of the participation degree of the overall evaluation element are given, and the impedance model participation analysis method is further corrected. And the theoretical analysis is verified through the numerical calculation result. The research result has certain reference significance for the participation analysis of the current heterogeneous power supply power system with the high-ratio converter.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The principles and embodiments of the present invention have been described herein with reference to specific examples, the description of which is intended only to assist in understanding the methods of the present invention and the core ideas thereof; also, it is within the scope of the present invention to be modified by those of ordinary skill in the art in light of the present teachings. In view of the foregoing, this description should not be construed as limiting the invention.
Claims (10)
1. The heterogeneous power supply power system correction participation analysis method based on the impedance model is characterized by comprising the following steps of:
Generating a full-system node impedance matrix and a full-system node admittance matrix of the heterogeneous power supply power system;
obtaining an admittance sensitivity factor of the element according to the full-system node impedance matrix and the full-system node admittance matrix;
Determining parameter sensitivity factors according to admittance sensitivity factors of the elements, and correcting the admittance sensitivity factors and the parameter sensitivity factors of the elements by considering the influence of transformer branch transformation ratios on an admittance matrix of all system nodes;
And according to the admittance sensitivity factor and the corrected parameter sensitivity factor of the element, the participation degree of the element is generally evaluated by using a three-layer chain rule.
2. The method for analyzing the revision participation of the heterogeneous power supply power system based on the impedance model according to claim 1, wherein the generating of the all-system node impedance matrix and the all-system node admittance matrix of the heterogeneous power supply power system specifically comprises:
the power supply node admittance matrix of the heterogeneous power supply power system is established as follows: wherein Y G represents the power supply node admittance matrix, Y G1 represents the admittance of the power supply connected with the bus 1, Y Gi represents the admittance of the power supply connected with the bus i, and Y Gn represents the admittance of the power supply connected with the bus n;
The network node admittance matrix of the heterogeneous power supply power system is established as follows: Wherein Y N represents a network node admittance matrix, Y N11 represents a self-admittance element of bus bar 1, Y N1j represents a transadmittance element between bus bar 1 and bus bar j, Y N1n represents a transadmittance element between bus bar 1 and bus bar n, Y Ni1 represents a transadmittance element between bus bar i and bus bar 1, Y Nij represents a transadmittance element between bus bar i and bus bar j, Y Nin represents a transadmittance element between bus bar i and bus bar n, Y Nn1 represents a transadmittance element between bus bar n and bus bar 1, Y Nnj represents a transadmittance element between bus bar n and bus bar j, and Y Nnn represents a self-admittance element of bus bar i;
Taking a disturbance value of injection current at each bus in a network as input, taking a voltage change value at each bus in the network as output, and obtaining a full-system node impedance matrix of the heterogeneous power supply power system by utilizing a reciprocal relation of a node admittance matrix and a node impedance matrix, wherein the full-system node impedance matrix comprises the following components: z= (I+Z NYG)-1ZN＝(YG+YN)-1＝Y-1) the total system node admittance matrix is Y=Y G+YN, wherein Z N is the node impedance matrix of the network, I is the identity matrix, Z represents the total system node impedance matrix taking the power supply impedance characteristic into account, and Y represents the total system node admittance matrix taking the power supply impedance characteristic into account.
3. The method for analyzing the correction participation of the heterogeneous power supply power system based on the impedance model according to claim 2, wherein the disturbance value of the injection current at each bus in the network is taken as an input, the voltage variation value at each bus in the network is taken as an output, and the node admittance matrix and the node impedance matrix reciprocal relationship are utilized to obtain the node impedance matrix of the whole system of the heterogeneous power supply power system, which is: z= (i+z NYG)-1ZN＝(YG+YN)-1＝Y-1, full system node admittance matrix is y=y G+YN, specifically including:
Taking the disturbance value of the injected current at each bus in the network as an input, the input is: Δu= [ Δi 1…ΔIi…ΔIn]T; wherein Deltau is input, deltaI 1 is a disturbance value of injection current at the bus 1, deltaI i is a disturbance value of injection current at the bus I, deltaI n is a disturbance value of injection current at the bus n;
taking the voltage change value at each bus of the network as output, the output is: Δy= [ Δu 1…ΔUi…ΔUn]T; wherein Δy is the output, Δu 1 is the voltage change value at bus 1, Δu i is the voltage change value at bus i, and Δu n is the voltage change value at bus n;
According to the input and the output, an input-output transfer function matrix is obtained by utilizing a closed-loop formula, wherein the input-output transfer function matrix is delta y= (I+Z NYG)-1ZN delta u;
according to the input and output transfer function matrix, determining a node impedance matrix of the whole system of the heterogeneous power supply power system as follows: z= (i+z NYG)-1ZN＝(YG+YN)-1＝Y-1, full system node admittance matrix is y=y G+YN.
4. The method for analyzing the correction participation of the heterogeneous power supply power system based on the impedance model according to claim 1, wherein the method for obtaining the admittance sensitivity factor of the element according to the impedance matrix of the nodes of the whole system and the admittance matrix of the nodes of the whole system specifically comprises the following steps:
according to the total system node impedance matrix and the total system node admittance matrix, a partial derivative formula for taking partial derivatives of admittances of power supplies or network elements in the system by the mode is obtained as follows: wherein, Representing the partial derivative of the modal lambda to the admittance Y of the power supply or network element, Y ij representing the transadmittance between busbar i and busbar j in the node admittance matrix Y taking into account the power supply impedance characteristics, tr () representing the trace of the matrix, res λ Z representing the taking of the remainder of the overall system node impedance matrix Z at modal lambda;
Expanding the partial conductance formula to obtain an expansion formula of the partial conductance formula, wherein the expansion formula comprises the following steps: Wherein Z ii represents the self-admittance element of bus bar i, Z jj represents the self-admittance element of bus bar j, Z ji is the impedance between bus bar j and bus bar i, and Z ij is the impedance between bus bar i and bus bar j; res λZii represents taking the remainder at modality λ for Z ii, res λZjj represents taking the remainder at modality λ for Z jj, res λZij represents taking the remainder at modality λ for Z ij, res λZji represents taking the remainder at modality λ for Z ji;
according to the expansion of the partial conductivity formula, the admittance sensitivity factors of the obtained element are as follows: where s λ,y is the admittance sensitivity factor of the element, which represents the conjugate transpose.
5. The heterogeneous power supply power system correction participation analysis method based on the impedance model according to claim 1, wherein the calculation formula of the parameter sensitivity factor is:
Wherein s λ,ρ is a parameter sensitivity factor, s λ,y is an admittance sensitivity factor of the element, For the sensitivity information of element admittance to parameters, < > represents the Frobenius inner product.
6. The method for analyzing the correction participation of the heterogeneous power supply power system based on the impedance model according to claim 4, wherein the admittance sensitivity factor and the parameter sensitivity factor of the correction element are taken into consideration in consideration of the influence of the transformer branch transformation ratio on the admittance matrix of the nodes of the whole system, and specifically comprise:
Considering the influence of the transformer branch transformation ratio on the partial conductance formula, the expansion of the correction partial conductance formula is as follows: wherein k is the transformer transformation ratio;
Correcting the admittance sensitivity factor of the element according to the expansion of the corrected partial conductance formula;
according to the admittance sensitivity factor of the element after correction, the formula is utilized Correcting the parameter sensitivity factor; where s 'λ,ρ is the modified parameter sensitivity factor, s' λ,y is the modified admittance sensitivity factor of the element,Sensitivity information for the element admittance to the parameters,Δρ is the small perturbation of the parameter ρ and y Δρ is the result of the change in admittance y of the power supply or network element after the change in the parameter ρ.
7. The method for analyzing the correction participation of the heterogeneous power supply power system based on the impedance model according to claim 1, wherein the correction element comprises an admittance sensitivity factor and a parameter sensitivity factor in consideration of the influence of the transformer branch transformation ratio on the admittance matrix of the whole system node, and further comprises:
According to the corrected parameter sensitivity factor and the formula delta lambda pr＝s′λ,ρ rho 5%, obtaining the modal prediction variable quantity caused by 5% of parameter change; wherein Δλ pr is the modal prediction variation caused by 5% change in parameter ρ, and s' λ,ρ is the modified parameter sensitivity factor;
and determining errors of the modal prediction variation and the modal actual variation caused by the parameter variation of 5%.
8. The heterogeneous power supply power system correction participation analysis method based on an impedance model according to claim 7, wherein the error determination formula is:
Wherein Error is the Error of the predicted mode change amount and the actual mode change amount caused by parameter change, deltalambda pr is the predicted mode change amount caused by 5% of parameter ρ change, deltalambda is the actual mode change amount after parameter change.
9. The method for analyzing the participation of the modified heterogeneous power supply power system based on the impedance model according to claim 1, wherein the participation degree of the elements is generally evaluated by using a three-layer chain rule according to the admittance sensitivity factor of the modified elements and the modified parameter sensitivity factor, and specifically comprises the following steps:
the first layer defining the three-layer chain law is: Wherein Δλ max is the maximum variation of the mode corresponding to the admittance variation of the element, and σ 2 and ω 2 are the real part and the imaginary part of the second layer of the three-layer chain rule, respectively;
According to the admittance sensitivity factor of the corrected element, defining a second layer of the three-layer chain rule as follows: σ 2+jω2＝<s′λ,y, y >; wherein s' λ,y is the admittance sensitivity factor of the modified element, y is the admittance of the power supply or network element, and < represents the Frobenius inner product;
according to the corrected parameter sensitivity factor, a third layer of the three-layer chain rule is defined as follows: wherein, For the modified parameter sensitivity factor, σ 3 and ω 3 are the real and imaginary parts, respectively, of the third layer of the three-layer chain law;
the degree of participation of the elements was generally evaluated using the first, second and third layers of the three-layer chain law.
10. The method for analysis of the modified participation of a heterogeneous power supply power system based on an impedance model according to claim 9, wherein the participation degree of the elements is generally evaluated by using a three-layer chain rule according to the admittance sensitivity factor of the modified elements and the modified parameter sensitivity factor, and further comprising:
And according to the corrected parameter sensitivity factor, aiming at enhancing the disturbance stability of the heterogeneous power supply power system, and determining the parameter adjustment direction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410957283.2A CN118916632A (en) | 2024-07-17 | 2024-07-17 | Heterogeneous power supply power system correction participation analysis method based on impedance model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410957283.2A CN118916632A (en) | 2024-07-17 | 2024-07-17 | Heterogeneous power supply power system correction participation analysis method based on impedance model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN118916632A true CN118916632A (en) | 2024-11-08 |
Family
ID=93304259
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410957283.2A Pending CN118916632A (en) | 2024-07-17 | 2024-07-17 | Heterogeneous power supply power system correction participation analysis method based on impedance model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118916632A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN119628086A (en) * | 2024-11-25 | 2025-03-14 | 华北电力大学 | A method, device, equipment, medium and product for dynamic gray box modal analysis of grid-following inverter |
| CN119787484A (en) * | 2024-12-24 | 2025-04-08 | 华北电力大学 | Extended impedance modal analysis method and device for internal dynamics of grid-connected inverter |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130226482A1 (en) * | 2012-02-24 | 2013-08-29 | Hongbo Sun | Decoupled Three-Phase Power Flow Analysis Method for Unbalanced Power Distribution Systems |
| CN112531731A (en) * | 2020-12-14 | 2021-03-19 | 国网河南省电力公司经济技术研究院 | Reactive compensation method for alternating current-direct current hybrid power distribution network based on sensitivity analysis |
| US20210247839A1 (en) * | 2018-09-03 | 2021-08-12 | Drexel University | On-chip voltage assignment through particle swarm optimization |
| CN113987795A (en) * | 2021-10-27 | 2022-01-28 | 华北电力大学 | A method for calculating modal sensitivity of DC/DC converters based on node admittance matrix |
| US20220155353A1 (en) * | 2019-03-13 | 2022-05-19 | Depsys Sa | Method for determining sensitivity coefficients of an electric power network using metering data |
| CN116937698A (en) * | 2023-08-21 | 2023-10-24 | 上海电力大学 | Small signal modeling method for power system with heterogeneous power supply |
| CN118100139A (en) * | 2024-01-22 | 2024-05-28 | 国电南瑞科技股份有限公司 | A method and system for dynamic partitioning of power grid considering electrical characteristics |
| CN118313320A (en) * | 2024-04-24 | 2024-07-09 | 上海交通大学 | Control-circuit hybrid network and frequency domain modal calculation method for multi-converter system |
-
2024
- 2024-07-17 CN CN202410957283.2A patent/CN118916632A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130226482A1 (en) * | 2012-02-24 | 2013-08-29 | Hongbo Sun | Decoupled Three-Phase Power Flow Analysis Method for Unbalanced Power Distribution Systems |
| US20210247839A1 (en) * | 2018-09-03 | 2021-08-12 | Drexel University | On-chip voltage assignment through particle swarm optimization |
| US20220155353A1 (en) * | 2019-03-13 | 2022-05-19 | Depsys Sa | Method for determining sensitivity coefficients of an electric power network using metering data |
| CN112531731A (en) * | 2020-12-14 | 2021-03-19 | 国网河南省电力公司经济技术研究院 | Reactive compensation method for alternating current-direct current hybrid power distribution network based on sensitivity analysis |
| CN113987795A (en) * | 2021-10-27 | 2022-01-28 | 华北电力大学 | A method for calculating modal sensitivity of DC/DC converters based on node admittance matrix |
| CN116937698A (en) * | 2023-08-21 | 2023-10-24 | 上海电力大学 | Small signal modeling method for power system with heterogeneous power supply |
| CN118100139A (en) * | 2024-01-22 | 2024-05-28 | 国电南瑞科技股份有限公司 | A method and system for dynamic partitioning of power grid considering electrical characteristics |
| CN118313320A (en) * | 2024-04-24 | 2024-07-09 | 上海交通大学 | Control-circuit hybrid network and frequency domain modal calculation method for multi-converter system |
Non-Patent Citations (2)
| Title |
|---|
| LE ZHENG 等: "Modal Analysis of Power System with High CIG Penetration B ased on Impedance Models", ARXIV, 1 June 2024 (2024-06-01), pages 1 - 8 * |
| 刘书铭;李陈莹;李琼林;崔雪;刘会金;: "电力系统串联谐波谐振的特性分析与灵敏度计算", 电力系统保护与控制, no. 09, 1 May 2015 (2015-05-01) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN119628086A (en) * | 2024-11-25 | 2025-03-14 | 华北电力大学 | A method, device, equipment, medium and product for dynamic gray box modal analysis of grid-following inverter |
| CN119787484A (en) * | 2024-12-24 | 2025-04-08 | 华北电力大学 | Extended impedance modal analysis method and device for internal dynamics of grid-connected inverter |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN118916632A (en) | Heterogeneous power supply power system correction participation analysis method based on impedance model | |
| CN110707686B (en) | The method of identifying the station area and the method of constructing the line topology of the station area | |
| CN107843810B (en) | A State Estimation-Based Online Fault Section Location Method for Active Distribution Networks | |
| CN109494724B (en) | LU decomposition-based large power grid Thevenin equivalent parameter online identification method | |
| CN114611676A (en) | New energy power generation system impedance model identification method and system based on neural network | |
| CN105490282B (en) | Consider the micro-capacitance sensor real-time voltage control method of the idle torque equilibrium degree of micro battery | |
| CN109993361A (en) | A kind of power distribution network operation trend prediction technique based on PMU | |
| CN108075480B (en) | State estimation method and system of AC/DC system | |
| CN108183484B (en) | A Decoupled Semilinear Optimal Power Flow Model Based on Hot Start Environment | |
| CN116031887B (en) | A method, system, device and medium for generating data of power grid simulation analysis examples | |
| CN115149582A (en) | A voltage stability analysis method of photovoltaic grid-connected power system based on monotonic system theory | |
| CN105071387B (en) | Based on the power distribution network Fast decoupled state estimation method that complex field standardization is theoretical | |
| CN107579525A (en) | A cold-start linearized optimal power flow calculation method capable of calculating complete power flow information | |
| CN109149583B (en) | A simple forward-backward generation power flow calculation method for active distribution network | |
| CN118885722A (en) | A Modal Analysis Method for Impedance Model of Power System | |
| Ma et al. | A two-step approach for multi-topology identification and parameter estimation of power distribution networks | |
| CN104392285B (en) | A kind of Optimal Power Flow Problems acquisition methods containing Hybrid HVDC | |
| CN113341275A (en) | Method for positioning single-phase earth fault of power distribution network | |
| CN110601167B (en) | Non-communication bus voltage non-deviation control method | |
| CN116718832A (en) | A high-precision electric energy measurement method and system for nationally produced devices | |
| CN110544960A (en) | A Distributed Control Method for Improving Reactive Power Sharing Capability of Island Microgrid | |
| CN113030540A (en) | Fundamental wave and harmonic wave electric energy bidirectional metering method of distributed new energy grid connection | |
| Yan et al. | Real-time carbon flow algorithm of electrical power systems based on network power decomposition | |
| CN105119301B (en) | A kind of micro-capacitance sensor theoretical based on equivalent feeder line is idle to divide equally control method | |
| CN112115568B (en) | An improved coupled least squares transmission line parameter identification method based on PMU measurement |
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 |