CN116992776B - Voltage source converter stability domain construction method and system based on piecewise affine - Google Patents
Voltage source converter stability domain construction method and system based on piecewise affine Download PDFInfo
- Publication number
- CN116992776B CN116992776B CN202311091069.5A CN202311091069A CN116992776B CN 116992776 B CN116992776 B CN 116992776B CN 202311091069 A CN202311091069 A CN 202311091069A CN 116992776 B CN116992776 B CN 116992776B
- Authority
- CN
- China
- Prior art keywords
- model
- voltage source
- source converter
- vsc
- impedance
- 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
Links
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 title claims abstract description 116
- 238000010276 construction Methods 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 claims abstract description 94
- 238000005457 optimization Methods 0.000 claims abstract description 71
- 230000008569 process Effects 0.000 claims abstract description 10
- 238000005192 partition Methods 0.000 claims description 27
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 230000014509 gene expression Effects 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 4
- 239000002131 composite material Substances 0.000 claims 2
- 150000001875 compounds Chemical group 0.000 claims 2
- 238000004590 computer program Methods 0.000 claims 2
- 238000004364 calculation method Methods 0.000 description 27
- 238000010586 diagram Methods 0.000 description 9
- 230000010355 oscillation Effects 0.000 description 9
- 238000011160 research Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 7
- 238000012937 correction Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明属于变流器技术领域，具体涉及一种基于分段仿射的电压源变流器稳定域构造方法及系统，包括：获取电压源变流器的阻抗模型；基于分段仿射和所获取的阻抗模型，构建电压源变流器的分段仿射阻抗模型；求解所构建的分段仿射阻抗模型，完成电压源变流器稳定域的构造；其中，在求解所构建的分段仿射阻抗模型的过程中，采用双层优化求解模型，双层优化求解模型中的内层模型结合所构建的分段仿射阻抗模型和奈奎斯特判据判断电压源变流器小干扰稳定域中工况点的稳定性，将内层模型作为双层优化求解模型中的外层模型的非线性约束，完成模型的求解，实现电压源变流器稳定域的构建。
The invention belongs to the technical field of converters, and specifically relates to a method and system for constructing a stable domain of a voltage source converter based on segmented affine, which includes: obtaining an impedance model of the voltage source converter; based on segmented affine and the obtained Obtain the impedance model to construct the segmented affine impedance model of the voltage source converter; solve the constructed segmented affine impedance model to complete the construction of the stable domain of the voltage source converter; among them, after solving the constructed segmented affine impedance model In the process of affine impedance model, a double-layer optimization solution model is used. The inner layer model in the double-layer optimization solution model is combined with the constructed piecewise affine impedance model and the Nyquist criterion to determine the small interference of the voltage source converter. For the stability of the operating point in the stable domain, the inner model is used as the nonlinear constraint of the outer model in the double-layer optimization solution model to complete the solution of the model and realize the construction of the stable domain of the voltage source converter.
Description
Technical Field
The invention belongs to the technical field of converters, and particularly relates to a voltage source converter stability domain construction method and system based on piecewise affine.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
With the rapid development of new energy, the access of high-proportion power electronic equipment easily causes the problem of small disturbance instability, restricts the transportation and the absorption of the new energy, and seriously threatens the safe operation of a power system.
New energy oscillation accidents often result from small disturbance negative damping instability. The existing stability analysis method aiming at the problems mainly adopts a state space method and an impedance method; however, the designated working conditions need to be analyzed, and it is difficult to comprehensively evaluate the overall operation stability of the system, so that a method capable of reflecting the stability of all the working conditions of the system needs to be established. Due to the time-varying nature of the grid, the stability domain needs to be updated in time according to the changes in the grid parameters. There is therefore a need for a stable domain construction method that can meet the requirements of online applications. At present, the research on the new energy full-working-condition stable domain is less, and two main challenges exist: the new energy station model has high order and complex impedance calculation; secondly, the construction of the stability domain boundaries is inefficient.
Aiming at the challenge of high order of the new energy station model, most researches adopt a reduced order or equivalent model when analyzing the new energy multi-machine model. Related researches adopt a balance realization theory to discard a non-dominant oscillation mode of a multi-voltage source converter (Voltage Source Converter, VSC for short) grid-connected system, but the non-dominant oscillation mode may cause partial coupling characteristic loss in a field; and when the operation working condition changes, the dominant oscillation mode reserved by the reduced order model is difficult to ensure to have consistent small disturbance characteristics. The related research adopts an equivalent method to equivalent the new energy station model to a single machine or multiple machine model; however, the equivalent method only focuses on the interaction mechanism of the whole new energy station and the power grid, and ignores the coupling characteristic between single machines.
The key to constructing a small perturbation stability domain is to determine its stability domain boundaries. The parameter stability domain is analyzed by related research, but the constructed stability domain changes along with the change of working conditions, so that guidance is difficult to provide in practical application. The related research adopts a point-by-point method to calculate the system stability under different serial complement values, and establishes a stable domain boundary of single parameters related to the serial complement values; however, the point-by-point method has the problem of large calculation amount, and is difficult to apply online. Related researches aim at micro-grid application scenes, active power stability domain boundaries of two VSCs are fitted based on kernel ridge regression, but the fitting method does not consider time-varying characteristics of the power grid. The related research is based on the offline construction of the stable domain boundary of the polymerization impedance and the prediction correction method, and compared with a traversal method, the calculation efficiency is improved; however, the calculation time of the adopted method is still long, and the method is difficult to be suitable for online application of stable domain solving.
Disclosure of Invention
In order to solve the problems, the invention provides a voltage source converter stable domain construction method and a system based on piecewise affine.
According to some embodiments, the first scheme of the invention provides a voltage source converter stability domain construction method based on piecewise affine, which adopts the following technical scheme:
a voltage source converter stable domain construction method based on piecewise affine comprises the following steps:
obtaining an impedance model of the voltage source converter;
constructing a segmented affine impedance model of the voltage source converter based on the segmented affine and the acquired impedance model;
solving the constructed piecewise affine impedance model to finish the construction of the stable domain of the voltage source converter;
in the process of solving the constructed segmented affine impedance model, a double-layer optimization solving model is adopted, an inner layer model in the double-layer optimization solving model is combined with the constructed segmented affine impedance model and a Nyquist criterion to judge the stability of the working point in the small interference stable domain of the voltage source converter, the inner layer model is used as the nonlinear constraint of an outer layer model in the double-layer optimization solving model, so that the solving of the model is completed, and the construction of the stable domain of the voltage source converter is realized.
As a further technical definition, in the process of constructing the piecewise affine impedance model of the voltage source converter, the acquired impedance model of the voltage source converter is converted into a piecewise low-order linear sub-model based on the piecewise affine theory.
As a further technical definition, the piecewise affine comprises constructing a voltage source converter impedance dataset, constructing dataset partitions according to a clustering algorithm, and implementing parameter estimation of an affine sub-model based on a least squares method.
As a further technical definition, the constructed piecewise affine impedance model of the voltage source converter is a dq coordinate system based on the voltage orientation of the grid-connected point of the voltage source converter itself; when a plurality of voltage source converters are connected in a grid, each constructed voltage source converter impedance model is transformed to the same dq coordinate system, and a multi-machine impedance model is obtained.
As a further technical limitation, in the process of judging the stability of the operating point in the small-interference stable domain of the voltage source converter, analyzing whether the return rate matrix meets the generalized nyquist criterion or not, namely, when the root tracks of two characteristic roots of the return rate matrix are not surrounded by (0, 1) on a complex plane, the operating point in the small-interference stable domain of the voltage source converter is stable; otherwise unstable.
As a further technical limitation, in the process of carrying out model solving, an inner layer model in the double-layer optimization solving model is solved by adopting an improved Nyquist stability criterion, and an outer layer model in the double-layer optimization solving model is solved by adopting a complex shape method.
As a further technical definition, the nonlinear constraint of the outer layer model in the double-layer optimization solution model is thatWherein P is VSC Is the active output of the voltage source converter, Q VSC Reactive power output of the voltage source converter; f (P) VSC ,Q VSC ) =0 satisfies the tide equation for the running operating point; p (P) min And P max The upper limit and the lower limit of the active output of the voltage source converter are respectively set; q (Q) min And Q max The reactive power output upper limit and the reactive power output lower limit of the voltage source converter are respectively; v rate Output rated voltage for voltage source converter port, v s The actual output voltage of the voltage source converter port; alpha is the power growth direction.
According to some embodiments, the second scheme of the invention provides a voltage source converter stability domain construction system based on piecewise affine, which adopts the following technical scheme:
a piecewise affine-based voltage source converter stability domain construction system comprising:
an acquisition module configured to acquire an impedance model of the voltage source converter;
a building block configured to build a segmented affine impedance model of the voltage source converter based on the segmented affine and the acquired impedance model;
the solving module is configured to solve the constructed piecewise affine impedance model to complete the construction of the stable domain of the voltage source converter;
In the process of solving the constructed segmented affine impedance model, a double-layer optimization solving model is adopted, an inner layer model in the double-layer optimization solving model is combined with the constructed segmented affine impedance model and a Nyquist criterion to judge the stability of the working point in the small interference stable domain of the voltage source converter, the inner layer model is used as the nonlinear constraint of an outer layer model in the double-layer optimization solving model, so that the solving of the model is completed, and the construction of the stable domain of the voltage source converter is realized.
According to some embodiments, a third aspect of the present invention provides a computer-readable storage medium, which adopts the following technical solutions:
a computer readable storage medium having stored thereon a program which when executed by a processor implements the steps in a method for piecewise affine-based voltage source converter stability domain construction according to the first aspect of the present invention.
According to some embodiments, a fourth aspect of the present invention provides an electronic device, which adopts the following technical solutions:
an electronic device comprising a memory, a processor and a program stored on the memory and executable on the processor, the processor implementing the steps in the piecewise affine-based voltage source converter stability domain construction method according to the first aspect of the present invention when the program is executed.
Compared with the prior art, the invention has the beneficial effects that:
the method is based on a piecewise affine theory and establishes a VSC impedance piecewise affine model; and a method for establishing a data set and a partition is provided for the VSC small disturbance model, so that the searching efficiency of affine coefficients and the accuracy of the affine model are ensured. And establishing a multi-machine VSC segmented affine model, retaining the frequency domain impedance information of the original system, and greatly reducing the model order.
The invention provides a construction method of a boundary of a VSC stable domain, which comprises the steps of constructing a double-layer optimization model for solving boundary points of the VSC stable domain, and solving the optimization model based on an improved generalized Nyquist stability criterion and a complex shape method; and then searching boundary points by changing the solving direction of the double-layer optimization model and performing segmentation fitting, so that the accuracy of the boundary of the stable domain is ensured and the construction efficiency is improved.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments and are incorporated in and constitute a part of this specification, illustrate and explain the embodiments and together with the description serve to explain the embodiments.
Fig. 1 is a diagram of a single VSC grid-connected structure in a first embodiment of the present invention;
FIG. 2 is a schematic illustration of a piecewise affine in accordance with one embodiment of the present invention;
FIG. 3 is a schematic diagram of a segmented affine partition according to a first embodiment of the invention;
fig. 4 is a diagram of a multi-VSC grid-connected structure in a first embodiment of the present invention;
FIG. 5 is a schematic diagram of a VSC small disturbance stability domain boundary according to a first embodiment of the present invention;
FIG. 6 is a schematic diagram of an improved generalized Nyquist stability criterion in accordance with a first embodiment of the present invention;
FIG. 7 is a flow chart of a solution for boundary conditions of a VSC stability domain based on a complex shape method according to a first embodiment of the invention;
FIG. 8 is a schematic diagram of iterative solution of a complex method in accordance with a first embodiment of the present invention;
FIG. 9 is a schematic diagram of a VSC small disturbance stability domain boundary solution in accordance with a first embodiment of the present invention;
FIG. 10 is a flow chart of an online quick construction of a VSC stability domain boundary in accordance with a first embodiment of the present invention;
FIG. 11 is a graph showing the stability domain accuracy contrast of the method and the spread point traversal method according to the first embodiment of the invention at different SCRs;
fig. 12 is a schematic diagram of single VSC stability domain boundary time domain verification at scr=1.7 in the first embodiment of the present invention;
FIG. 13 is a comparison of multiple VSC stability domains under different SCRs in accordance with an embodiment of the present invention
Fig. 14 is a schematic diagram of multi-machine VSC stability domain boundary time domain verification at scr=1.7 in the first embodiment of the present invention.
Detailed Description
The invention will be further described with reference to the drawings and examples.
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the present application. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present invention. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
Embodiments of the invention and features of the embodiments may be combined with each other without conflict.
Example 1
The embodiment of the invention introduces a voltage source converter stability domain construction method based on piecewise affine.
The embodiment adopts a piecewise affine theory to convert VSC impedance into a linear form, and simultaneously provides a stable domain construction method capable of meeting online application; for complex high-order VSC impedance small disturbance models, the segmentation method of segmented affines and affine sub-model parameter identification are challenging. And, how to further improve the stability domain boundary rapid construction efficiency based on the piecewise affine impedance model is also challenging.
According to the VSC model shown in fig. 1, which mainly comprises a current loop control and a phase-locked loop control of the current transformer, an impedance model is established, which analyzes small disturbance stability. In the dq coordinate system, the two-dimensional matrix of VSC impedance is as follows:
wherein P is VSC And Q VSC Active and reactive power of VSC output, v s Representing the grid tie point voltage. In the frequency domain s=j2pi f is a complex state variable, where f is frequency.
The VSC impedance characteristics are affected by the operating point, and are related to active, reactive and port voltages. In order to construct a VSC segmented affine model more suitable for practical applications, the running operating point state variables are defined as:
x＝[P VSC Q VSC v s ] (2)
as shown in fig. 2, the complex high-order VSC impedance model with respect to the operating points x and s is converted into a piecewise low-order linear submodel based on piecewise affine theory. The piecewise affine specifically comprises three steps: constructing a VSC impedance data set, constructing a data set partition according to a clustering algorithm, and realizing parameter estimation of an affine submodel based on a least square method.
1) Impedance data set construction method
Since the input variables have different units and ranges, normalization processing is performed on all the input variables x. And calculating the corresponding VSC output impedance under different x and s within the operable range, and constructing a sample data set.
The present embodiment proposes to construct the VSC impedance dataset by linearity analysis to avoid the influence of sample set density on the accuracy and computational efficiency of the segmented affine model.
Formula (1) relates to P VSC ，Q VSC ，v s And s have different linearities. With impedance element Z dd For the example of linearity of s, a linear correlation coefficient r is defined as:
where Cov is the covariance calculation symbol and Var is the standard deviation calculation symbol. The molecules of the linear correlation coefficient are Z dd (s) real and imaginary covariance. The denominator of the linear correlation coefficient is Z dd (s) the product of the real and imaginary standard deviations. The linear correlation coefficient is between 0 and 1, the closer it is to 1, the higher the linearity of the two variables.
According to the linearity of the VSC working points x and s with respect to the impedance, scattering is performed when constructing a sample set, and the sample space is preliminarily divided into different data sets.
2) Data clustering and partition space solving
The data set obtained in the step 1) contains a large number of sample points, and a data clustering method is adopted to divide a sample space. A data set space may be divided into disjoint subspaces, each of which may be represented by a low-order linear analytical formula.
In the embodiment, a hierarchical clustering method is adopted to cluster the data of the data set. The main idea is that initially each data is considered as a class, the Euclidean distance is used to calculate the sample class S 1 ，S 2 Is a distance of (3). In Z dd (x, s) for example, the Euclidean distance e between its samples is:
the maximum sample distance between the two classes is taken as the true distance E (S 1 ,S 2 )：
When E (S) 1 ,S 2 ) Beyond the set threshold value psi, aggregation between the two classes is not performed, and the number and the accuracy of clusters are adjusted by adjusting the size of the threshold value psi.
After determining the sample set partition based on the clustering algorithm, solving the interface coefficient of the partition by adopting support vector classification (support vector classification, SVC) [21] . With χ u The partitions of the running points x and s of the piecewise affine model are represented, the subscript u represents the number of the partition, and the expression of the interface can be expressed as:
χ u ＝{F u [x s] T +g u ，，0} (6)
wherein F is u And g u Is an interface coefficient matrix.
As shown in FIG. 3, in Z dd Visual partition indications are given for the examples. When the parameters of the power grid are determined, the port voltage can be determined according to the active and reactive power flow, s is taken to be 5-15Hz, P VSC At 0-25kW, Q VSC Between 0 and 10kVar as a data set. The segmented affine partitioning result is shown in fig. 3, and the hierarchical clustering algorithm divides the sample space into 4 pieces, and the interface coefficient matrix of each subspace is obtained according to SVC.
3) Affine sub-model parameter estimation
Based on the data in the partition of step 2), the sub-model parameter beta is identified by utilizing a parameter estimation method u . After decomposing the dataset according to step 1) and clustering the partitions according to step 2), the sub-partitions are then partitionedThe linearity of the internal data is high, and a first-order linear expression fitting can be adopted. In Z dd For example, its piecewise affine model is expressed as:
wherein beta is u ＝[β u,1 β u,2 β u,3 β u,4 β u,5 ]Coefficient parameter of the u-th affine sub-model, χ u And the U state variable partition space is represented, and U is the total number of partitions in the operation point and the full frequency domain space.
Realizing parameter beta in (7) based on least square method u Is a single-chip microcomputer. For the (u) th χ u The v-th sample point x within the partition u,v Defining an error distance as follows:
where v is the sample point number. By making J u The sum of squares of the distances from each sample point to the fitting equation is minimum, and the optimal affine sub-model parameter beta is obtained u ：
Each partition executes least square fitting, and a segmented affine impedance model of the VSC under the full working condition and the full frequency band can be obtained, which is expressed as follows:
wherein the VSCThe other three terms of the impedance matrix are equally available from (7).
And establishing a piecewise affine impedance model about the operation working points x and s for the VSC, converting the high-order nonlinear model into a partitioned first-order piecewise affine model, and greatly reducing the impedance order of the VSC. In practical application, the port characteristic of the VSC can be obtained according to the load flow operation, and then the partition submodel is searched to obtain a low-order impedance expression.
When the VSCs are connected, each VSC impedance model is built based on a dq coordinate system oriented by the voltage of the self-connected point. Because of the influence of the current collecting circuit, the voltages of the grid connection points are different, and the grid connection points cannot be directly connected to form an impedance network. The multi-machine impedance model is built by transforming each VSC impedance model into a unified dq coordinate system. The multi-machine VSC grid-connected structure is shown in fig. 4, wherein grid nodes are used as reference nodes, and the rest VSC nodes are used as PQ nodes. The operating point x of each VSC is obtained through load flow calculation, and meanwhile, the angle difference theta between the VSC node and the reference node can also be obtained VSC,k 。
According to (10), the impedance model of the kth VSC transformed into the unified dq coordinate system is:
wherein x is VSC,k Is an operating point state variable, T, of the kth VSC defined according to equation (2) dq,k Is defined as follows:
the line impedance and the output impedance of the VSC are adopted for series-parallel connection to obtain the aggregate impedance Z of the multi-VSC field farm Finally, a multi-VSC segmented affine impedance model can be obtained,
wherein x is farm ＝[x VSC,1 ,…,x VSC,k ,…,x VSC,K ] (14)
The impedance characteristics are accurately reflected according to x and s, so that the effectiveness of actual operation is greatly improved.
In a multi-VSC application, the impedance characteristics of each VSC may directly determine the low-order impedance expression based on the port characteristics. The problem that when the working condition changes, the high-order impedance of each VSC needs to be recalculated is avoided. And as the number of VSCs becomes larger, the advantage of the piecewise affine model is more obvious, so that the computational complexity can be greatly reduced.
The key to constructing a stability domain is to quickly determine a stability boundary, which consists of a series of operating points at critical steady state. Therefore, the embodiment provides a double-layer optimization model for quickly solving the boundary points of the stable domain, so that the construction efficiency and accuracy of the boundary of the small-interference stable domain of the VSC system can be improved.
As shown in fig. 5, the VSC small disturbance stability domain is surrounded by a series of discrete operating condition critical points (black dots) on the boundary. And judging the stability of each working condition point according to the generalized Nyquist criterion. VSC grid-connection stability can be determined by analyzing whether the return rate matrix L(s) satisfies the generalized nyquist criterion: on the complex plane, when two characteristic roots lambda of L(s) 1 ，λ 2 When none of the trajectory curves of (1, 0) is enclosed, the system is stable; otherwise, the system is unstable. Wherein the recovery matrix is expressed as:
wherein the equivalent impedance at the grid side is denoted as Z grid 。
In fig. 5 d represents the power limit for safe and stable operation in the power increase direction α. Solving the stability domain boundary problem may translate into solving a critical stable operating point for each power increase direction.
Aiming at the problems, the embodiment provides a double-layer optimization solving model which is divided into an inner layer model and an outer layer model: solving the stability of each working point by the inner layer model based on a generalized Nyquist criterion; and solving the power limit which can stably run under the power increasing direction by the outer layer model. The inner layer model needs to determine the operation condition for solving, and the power increase limit of the outer layer model is limited by the stability margin of the inner layer model.
1) Inner layer optimization model
g RE For when the imaginary part of the characteristic root trace is equal to 0, i.e. Im (lambda) 1,2 (s,P VSC ,Q VSC ) For example, =0), the abscissa of the smallest intersection point on the real axis. When g RE When the total weight of the total weight is greater than-1, the system is stable, g RE And the system is critically stable when being equal to-1. Otherwise, the system is unstable. The decision variable of the inner layer model is P VSC ，Q VSC And s, the optimization target is according to the characteristic root lambda 1 ，λ 2 Solving g by trajectory curve of (2) RE . The inner layer optimization problem thus constructed is:
wherein f (P VSC ,Q VSC ) =0 satisfies the tide equation for the running operating point; p (P) min And P max The upper and lower limits of active force of the VSC; q (Q) min And Q max The reactive output upper and lower limits of the VSC are respectively set; v rate The rated voltage is output to the VSC port, and the VSC is required to meet the condition that the port voltage is not out of limit during operation.
According to the optimization target g RE Defining stability margin M of a system G ：
M G (s,P VSC ,Q VSC )＝-20log(|g RE (s,P VSC ,Q VSC )|) (17)
2) Outer layer optimization model
The optimization goal of the stability domain boundary outer layer optimization model is the maximum power limit d along a certain power increase direction, and the decision variable is P VSC And Q VSC . The outer layer optimization problem thus constructed is:
wherein M is G Is coming fromStability margin constraints solved from the inner layer model. According to equation (18), after the power growth direction is determined, the problem of solving the boundary point of the critical stable domain can be converted into the maximum active and reactive output under the condition that the solution meets the constraint condition.
In solving the double-layer optimization problem, the inner layer and the outer layer are mutually limited, which greatly increases the calculation time. The section proposes to construct an improved generalized Nyquist criterion suitable for double-layer optimization model calculation, which mainly solves the problem of long solving time of an inner-layer model.
As shown in fig. 6, an improved generalized nyquist stabilization criterion is proposed, delineating forbidden viable domains. And a rectangular area with a real axis less than xi and a height epsilon taking the real axis as a center on the complex plane is shown as the following formula:
compared with the Middlebrook criterion and the GMPM criterion, the criterion method provided by the embodiment has lower conservation in stability judgment and is more suitable for practical application. When ε=0, the forbidden feasible region becomes the forbidden region of generalized Nyquist. The improved nyquist criterion is shown in formula (20):
when the nyquist curve passes through the red region, it is noted as-1, indicating instability. Otherwise, it is denoted as 1, indicating stability.
Latin hypercube sampling is carried out on s, and g is avoided by judging whether sampling points are located in forbidden domains or not RE Is a solution to the problem. And the symbol operation is converted into the numerical value calculation by utilizing a sampling point mode, so that the calculation efficiency can be remarkably improved.
Introducing an improved generalized Nyquist criterion according to the double-layer optimization model and the inner-layer optimization model, and introducing M label (P VSC ,Q VSC ) Substituting the constraints into the outer layer optimization model. Thereby changing the inner layer optimizing model into the outer layer optimizing modelNonlinear constraint of the form, formula (18) becomes the following formula:
the complex method is a nonlinear global optimization algorithm and does not need an explicit analytical expression of an optimization target or constraint condition. In solving the problem of equation (21), the complex method may better incorporate equation (20) into the constraint. The algorithm flow chart shown in fig. 7 mainly takes n points, and finds the worst solution and the best solution of the target result during iteration. The worst solution is then mapped to the center of the remaining n-1 points. If the mapped target result does not meet the limit of the optimization problem or the target result is worse, the mapping distance is reduced. When the errors of the n points are all iteratively converged to be smaller than the set threshold value, the iteration is stopped. n is typically 10 times the number of variables.
P with a constant power growth direction in the application of the present embodiment VSC And Q VSC Is fixed, thus taking the variable P in the complex method VSC And n=10. The principle of iterative solution of the complex method shown in FIG. 8, P w Representing the worst solution in this iteration, P b Represents the optimal solution, P c Representing the centroid, calculated according to the following formula:
P map is P w According to P c The mapped points are calculated according to the following equation:
P map ＝P c +γ(P c -P w ) (23)
where γ represents a mapping coefficient, and is set to 1.3 in this embodiment.
For P in 10 points in each iteration w And mapping, wherein when the 10 points are subjected to iterative mapping and all meet the set threshold range, iteration is stopped, and an optimal solution is obtained.
As shown in FIG. 9, in the construction of VWhen SC stabilizes a domain boundary, a series of discrete points on the boundary are required, each of which can be found by the complex method. The search is performed by a step size of delta alpha over a range of angles alpha (-90 deg., 90 deg.). Where Δα is searched in initially large angular increments. As indicated by blue arrows, each two adjacent discrete points are calculated at alpha i Boundary point (P) in +Δα/2 power increase direction check ,Q check ) As an error detection point. The linearization analytical formula for obtaining the boundary according to the ith and the (i+1) th points on the boundary is as follows:
calculation (P) check ,Q check ) Distance d to line segment shown in (23) error The method comprises the following steps:
when 0 is<d error <d limit When the boundary line segment is considered to meet the set boundary precision; when d error >d limit At this time, the angle increment Δα=Δα/2 is adjusted, and then d is calculated error Until it meets the boundary error accuracy.
Error of the j-th line segment as shown in fig. 9. When d error <At 0, the boundary of the stable region formed at this time has a small unstable region inside. At this time, an improved generalized Nyquist stability criterion is utilized for the inner layer optimization model. By increasing the rectangular forbidden domain shown in FIG. 6, i.e., increasing the values of ε and ζ, the stable domain is made more conservative, up to 0 <d error <d limit The occurrence of unstable small areas is avoided.
As shown in fig. 10, specific steps of an online rapid construction method for a stability domain boundary of a VSC piecewise affine impedance model are given: first, data is initialized, and a starting angle and an increment thereof are set. And then solving boundary points based on a complex shape method, and performing error detection between every two points. And correcting the points which are not satisfied by the errors. And finally, obtaining an analytic expression by segment fitting to obtain a complete boundary. By reasonably changing the growth direction of the search power, the stable domain can be ensured to have higher calculation efficiency and more accurate boundary, and the online construction is completed.
Calculation case analysis
The single-machine VSC simulation model shown in the figure 1 is built, and compared with the existing method from the two angles of accuracy and calculation efficiency of single-machine VSC stable domain construction respectively; firstly, comparing a point-by-point traversal method with a double-layer optimization method proposed by the embodiment, and comparing the boundary differences of stable domains under different short circuit ratios (Short circuit ratio, SCR); and when scr=1.7, stability domain boundary validity verification is performed. Secondly, comparing the calculation time of different methods, verifying the stability domain construction method based on the VSC segmented affine model has high calculation efficiency.
The method and the point-by-point traversal method provided by the embodiment are adopted to establish the stable domains of the active power and the reactive power of the VSC. As shown in fig. 11, the boundary errors are compared in three cases of scr=1.4, scr=1.7, and scr=2.0, respectively. As can be seen, the point-by-point traversal method and the method proposed in this example yields consistent boundaries with minimal errors. Therefore, the method provided by the embodiment can accurately search the stability domain boundary. By comparing the stability domain boundaries under different SCRs, the smaller the grid SCR, the smaller the stability domain the system can operate. For a given example, the change rule of active power and reactive power is observed, which shows that the VSC under the weak current network can appear under the unstable condition along with the increase of active power, and the stability of the system can be improved by properly increasing reactive power.
Three operating points, operating 1, are shown in FIG. 11: p (P) VSC ＝22kW，Q VSC =0 kVar; working condition 2: p (P) VSC ＝23kW，Q VSC =0 kVar; working condition 3: p (P) VSC ＝23kW，Q VSC =2kvar. As shown in fig. 12, the stability under these three conditions was compared in the simulation. Before t is less than 5s, the system is operated under the working condition 1 and is in stable operation. The system increases the active power output at 5 s. At 6.5s the system is at regime 2 and the system begins to oscillate at 9Hz, with a subsequent gradual increase in oscillation amplitude. Reactive output increase at 7.5sUp to 2kVar, the system is at regime 3. The oscillation amplitude then begins to decrease and the system gradually resumes steady state operation. The accuracy of the boundary of FIG. 11 is verified by the time domain simulation result, and in practical application, the operation condition can be adjusted according to the stability domain, so that the stable operation of the system is ensured.
Aiming at four methods of a point-by-point traversal method, a prediction correction method, a double-layer optimization method of a VSC full-order model and a double-layer optimization method of a VSC segmented affine model, the search efficiency of a comparison boundary is ensured under the same search precision. The calculated time comparisons for the four methods are shown in Table 1:
TABLE 1 grid-connected stability Domain boundary calculation time comparison
The boundary search method has an advantage in the calculation speed of each operating point over the point-by-point traversal method and the prediction correction method. And in the number of working condition points to be calculated, the traversal method is the most, the prediction correction is less, and the embodiment has the advantage in calculation speed by the least calculation of the optimization method. The method 4 greatly improves the speed of solving the impedance during calculation by adopting the piecewise affine model, improves the calculation efficiency by more than 2 times compared with the third method, and proves the effectiveness of adopting the piecewise affine model. The calculation efficiency (inverse proportion of calculation time) of the comparison traversal method and the prediction correction method in the method proposed by the embodiment is 54 and 6.45, and the calculation time meets the requirement of on-line operation.
The method of the embodiment supports the online construction of the stable domains of a plurality of VSCs. And building a four-machine VSC simulation model according to the VSC parallel structure shown in FIG. 4. Stability between VSCs is analyzed under a weak current network, and for convenience and intuitionistic presentation, the stability domain of active power between the branch 1 and the branch 2 is focused in verification of the embodiment. Setting the reactive power output of 4 VSCs to 0, fixing the active power output of the third VSC to 23kW, and the active power output of the fourth VSC to 20kW. As shown in fig. 13, P at scr=1.6, scr=1.7, and scr=1.8 was compared VSC1 And P VSC2 Is a stable domain of operation of (a). The smaller the SCR, the smaller the operational stability domain. Three working conditions are set, namely, working condition 1: p (P) VSC1 ＝14kW，P VSC2 =20 kW; working condition 2: p (P) VSC1 ＝16kW，P VSC2 =20 kW; working condition 3: p (P) VSC1 ＝16kW，P VSC2 ＝17kW。
As shown in FIG. 14, the accuracy of the stability domain boundary was verified by simulation. The system is in a steady state at condition 1, the active power of the VSC connected to branch 1 begins to increase at 5s, condition 2 is reached at 7.5s, and the system begins to generate an oscillation frequency of 8 Hz. And at 8s, the VSC connected through the branch circuit 2 reduces the power output to the working condition 3, and the system gradually recovers to be stable. The time domain simulation result verifies the effectiveness and accuracy of the method in the multi-machine VSC stability domain structure. In actual operation, a stable domain among the multi-machine interconnected VSCs can be constructed, and the output of each VSC is distributed through field control, so that the system is in a stable state. If oscillation occurs, the output of each VSC can be reasonably regulated and controlled to enable the system to restrain the oscillation.
As shown in table 2, the computational efficiency of the four-machine VSC stability domain boundary configuration was compared with method 3 and method 4. The efficiency ratio is calculated to be 5.98 by a double-layer optimization method based on the segmented affine impedance model and the impedance full-order model. The necessity of a segmented affine model to be applied to the stable domain boundary construction is demonstrated. The application of the piecewise affine impedance model can increase the computational efficiency by a factor of 2 compared to a single VSC stability domain configuration, indicating that the piecewise affine impedance model has more significant advantages in a larger scale system.
Table 2 4 VSC grid-connected stability area boundary calculation time comparison
The embodiment establishes a VSC impedance piecewise affine model based on piecewise affine theory; providing a method for establishing a data set and a partition for the VSC small disturbance model, and ensuring the searching efficiency of affine coefficients and the accuracy of the affine model; and establishing a multi-machine VSC segmented affine model, retaining the frequency domain impedance information of the original system, and greatly reducing the model order.
The embodiment provides a construction method of a boundary of a VSC stability domain, a double-layer optimization model for solving boundary points of the VSC stability domain is constructed, and the optimization model is solved based on an improved generalized Nyquist stability criterion and a complex method; and then searching boundary points by changing the solving direction of the double-layer optimization model and performing segmentation fitting, so that the accuracy of the boundary of the stable domain is ensured and the construction efficiency is improved.
The embodiment verifies the accuracy of the piecewise affine model by comparing the errors with the impedance full-order model. And the accuracy and the effectiveness of the method in the aspect of constructing the stable domain are verified through single VSC and multi-VSC simulation. Compared with the existing method, the method in the embodiment can remarkably improve the calculation efficiency.
Example two
The second embodiment of the invention introduces a voltage source converter stability domain construction system based on piecewise affine.
A piecewise affine-based voltage source converter stability domain construction system comprising:
an acquisition module configured to acquire an impedance model of the voltage source converter;
a building block configured to build a segmented affine impedance model of the voltage source converter based on the segmented affine and the acquired impedance model;
the solving module is configured to solve the constructed piecewise affine impedance model to complete the construction of the stable domain of the voltage source converter;
in the process of solving the constructed segmented affine impedance model, a double-layer optimization solving model is adopted, an inner layer model in the double-layer optimization solving model is combined with the constructed segmented affine impedance model and a Nyquist criterion to judge the stability of the working point in the small interference stable domain of the voltage source converter, the inner layer model is used as the nonlinear constraint of an outer layer model in the double-layer optimization solving model, so that the solving of the model is completed, and the construction of the stable domain of the voltage source converter is realized.
The detailed steps are the same as those of the voltage source converter stability domain construction method based on piecewise affine provided in the first embodiment, and are not described herein again.
Example III
The third embodiment of the invention provides a computer readable storage medium.
A computer readable storage medium having stored thereon a program which when executed by a processor performs the steps in a method for piecewise affine-based voltage source converter stability domain construction according to an embodiment of the present invention.
The detailed steps are the same as those of the voltage source converter stability domain construction method based on piecewise affine provided in the first embodiment, and are not described herein again.
Example IV
The fourth embodiment of the invention provides electronic equipment.
An electronic device comprising a memory, a processor and a program stored on the memory and executable on the processor, the processor implementing the steps in a method for constructing a voltage source converter stability domain based on piecewise affine according to an embodiment of the invention when the program is executed.
The detailed steps are the same as those of the voltage source converter stability domain construction method based on piecewise affine provided in the first embodiment, and are not described herein again.
The above description is only a preferred embodiment of the present embodiment, and is not intended to limit the present embodiment, and various modifications and variations can be made to the present embodiment by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present embodiment should be included in the protection scope of the present embodiment.
Claims (5)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311091069.5A CN116992776B (en) | 2023-08-28 | 2023-08-28 | Voltage source converter stability domain construction method and system based on piecewise affine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311091069.5A CN116992776B (en) | 2023-08-28 | 2023-08-28 | Voltage source converter stability domain construction method and system based on piecewise affine |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116992776A CN116992776A (en) | 2023-11-03 |
| CN116992776B true CN116992776B (en) | 2024-03-26 |
Family
ID=88524774
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311091069.5A Active CN116992776B (en) | 2023-08-28 | 2023-08-28 | Voltage source converter stability domain construction method and system based on piecewise affine |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116992776B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103050970A (en) * | 2013-01-15 | 2013-04-17 | 华北电力大学 | Stability analyzing and optimizing method suitable for layering and zoning of ultra-high voltage electric network |
| CN108681640A (en) * | 2018-05-17 | 2018-10-19 | 武汉大学 | A kind of hybrid simulation system, topology and emulation mode based on virtual impedance compensation |
| CN110601244A (en) * | 2019-07-19 | 2019-12-20 | 华北电力大学（保定） | Compensation control method for reducing influence of reactive power fluctuation on system stability |
| WO2021012298A1 (en) * | 2019-07-25 | 2021-01-28 | 东北大学 | Self-mutual-group multi-level stability identification and stability recovery method for multi-port energy router |
| CN112736926A (en) * | 2020-12-25 | 2021-04-30 | 华北电力大学 | Interval affine power flow dynamic optimization method for distributed new energy access power distribution network |
| CN112909924A (en) * | 2021-01-22 | 2021-06-04 | 武汉大学 | A new energy power system small disturbance stability domain acquisition method in the combination space of operation and control parameters |
| CN113011083A (en) * | 2021-02-25 | 2021-06-22 | 中国科学院电工研究所 | Simulation evaluation method for island operation time length of comprehensive energy system |
| CN113435134A (en) * | 2021-06-23 | 2021-09-24 | 山东大学 | Wind power grid connection stability determination method and system based on full dynamic impedance model |
| CN114389282A (en) * | 2022-01-21 | 2022-04-22 | 福州大学 | Affine-based uncertainty modal analysis method |
| CN114997505A (en) * | 2022-06-16 | 2022-09-02 | 华北电力大学 | A method and system for stability analysis of grid-connected power supply via converter |
| CN115036931A (en) * | 2022-05-20 | 2022-09-09 | 武汉工程大学 | Active power grid reactive voltage affine adjustable robust optimization method and device |
-
2023
- 2023-08-28 CN CN202311091069.5A patent/CN116992776B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103050970A (en) * | 2013-01-15 | 2013-04-17 | 华北电力大学 | Stability analyzing and optimizing method suitable for layering and zoning of ultra-high voltage electric network |
| CN108681640A (en) * | 2018-05-17 | 2018-10-19 | 武汉大学 | A kind of hybrid simulation system, topology and emulation mode based on virtual impedance compensation |
| CN110601244A (en) * | 2019-07-19 | 2019-12-20 | 华北电力大学（保定） | Compensation control method for reducing influence of reactive power fluctuation on system stability |
| WO2021012298A1 (en) * | 2019-07-25 | 2021-01-28 | 东北大学 | Self-mutual-group multi-level stability identification and stability recovery method for multi-port energy router |
| CN112736926A (en) * | 2020-12-25 | 2021-04-30 | 华北电力大学 | Interval affine power flow dynamic optimization method for distributed new energy access power distribution network |
| CN112909924A (en) * | 2021-01-22 | 2021-06-04 | 武汉大学 | A new energy power system small disturbance stability domain acquisition method in the combination space of operation and control parameters |
| CN113011083A (en) * | 2021-02-25 | 2021-06-22 | 中国科学院电工研究所 | Simulation evaluation method for island operation time length of comprehensive energy system |
| CN113435134A (en) * | 2021-06-23 | 2021-09-24 | 山东大学 | Wind power grid connection stability determination method and system based on full dynamic impedance model |
| CN114389282A (en) * | 2022-01-21 | 2022-04-22 | 福州大学 | Affine-based uncertainty modal analysis method |
| CN115036931A (en) * | 2022-05-20 | 2022-09-09 | 武汉工程大学 | Active power grid reactive voltage affine adjustable robust optimization method and device |
| CN114997505A (en) * | 2022-06-16 | 2022-09-02 | 华北电力大学 | A method and system for stability analysis of grid-connected power supply via converter |
Non-Patent Citations (3)
| Title |
|---|
| Mohammad Momeni ; Homayoun Meshgin Kelk ; Heidarali Talebi.Rotating Switching Surface Control of SeriesResonant Converter Based on a Piecewise Affine Model.IEEE Transactions on Power Electronics.2014,第1762-1772页. * |
| 基于P-Q控制模式的VSC-HVDC小干扰稳定性研究;张志军;李帆;孟庆波;;电气传动(08);第55-60+66页 * |
| 多变流器并网发电系统控制策略与稳定性研究;陈浩;中国优秀硕士学位论文全文数据库;第C042-1010页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116992776A (en) | 2023-11-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | A review of applications of artificial intelligent algorithms in wind farms | |
| Gionfra et al. | Wind farm distributed PSO-based control for constrained power generation maximization | |
| Li et al. | A sequential approach for small signal stability enhancement with optimizing generation cost | |
| CN109981103B (en) | Parameter optimization method and system for bi-second order generalized integrator frequency locking loop | |
| CN106408135A (en) | Power system optimal power flow method based on feedback learning cuckoo algorithm | |
| Sibtain et al. | Frequency stabilization for interconnected renewable based power system using cascaded model predictive controller with fractional order PID controller | |
| CN114759590A (en) | Converter station steady-state reactive power coordination strategy optimization method and device based on data driving | |
| CN116646976A (en) | A distributed photovoltaic cluster division method and system based on agglomerative hierarchical clustering algorithm | |
| Wang et al. | Wind power curve modeling with large-scale generalized kernel-based regression model | |
| Bahmani et al. | Wind turbines power regulation using a low-complexity linear parameter varying-model predictive control approach | |
| Yang et al. | Intelligent computing of positive switching impulse breakdown voltage of rod-plane air gap based on extremely randomized trees algorithm | |
| Guo et al. | Parallel multi-fidelity expected improvement method for efficient global optimization | |
| CN116992776B (en) | Voltage source converter stability domain construction method and system based on piecewise affine | |
| Chen et al. | Generalized predictive control application scheme for nonlinear hydro-turbine regulation system: Based on a precise novel control structure | |
| CN109611274B (en) | LQG (Linear quadratic glass) optimization control method for high wind speed area of wind generating set | |
| CN110556878A (en) | Distributed voltage control optimization method and system applied to wind power plant | |
| Ma et al. | Sensorless intelligent second-order integral sliding mode maximum power point tracking control for wind turbine system based on wind speed estimation | |
| CN110277783A (en) | An online active de-loading control method and system for a power system | |
| CN104408531A (en) | Uniform dynamic planning method for multi-dimensional multi-stage complex decision problem | |
| Shalbafian et al. | A novel robust nonlinear optimal second-order sliding mode control scheme for power optimization of wind energy conversion systems | |
| Xu et al. | Advanced control techniques for PV maximum power point tracking | |
| CN116861321A (en) | Wind farm stability judging method and system based on stable domain quantification | |
| CN117150898A (en) | Power transmission section limit transmission capacity assessment method based on parameter optimization ensemble learning | |
| CN114678888A (en) | New energy external characteristic drive source network coordination stability judgment method based on passivity index | |
| Qiao et al. | Transient stability assessment for ac-dc hybrid systems based on bayesian optimization xgboost |
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 |