CN117638867A - Method and system for evaluating interaction strength between new energy station units - Google Patents
Method and system for evaluating interaction strength between new energy station units Download PDFInfo
- Publication number
- CN117638867A CN117638867A CN202311372695.1A CN202311372695A CN117638867A CN 117638867 A CN117638867 A CN 117638867A CN 202311372695 A CN202311372695 A CN 202311372695A CN 117638867 A CN117638867 A CN 117638867A
- Authority
- CN
- China
- Prior art keywords
- new energy
- equipment
- energy station
- matrix
- model
- 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
- 230000003993 interaction Effects 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000006870 function Effects 0.000 claims abstract description 58
- 239000011159 matrix material Substances 0.000 claims abstract description 49
- 238000012546 transfer Methods 0.000 claims abstract description 32
- 238000010276 construction Methods 0.000 claims abstract description 16
- 238000004590 computer program Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 5
- 239000000758 substrate Substances 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract description 6
- 230000010355 oscillation Effects 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 description 11
- 238000005259 measurement Methods 0.000 description 4
- 230000006855 networking Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Abstract
The invention discloses a method and a system for evaluating interaction strength among new energy station units, which relate to the technical field of new energy station unit strength evaluation and comprise the steps of collecting input and output of equipment grid connection points, and establishing a power current model of electromechanical or electromagnetic scale of follow-up equipment and a power voltage model of electromechanical or electromagnetic scale of network construction equipment in a frequency domain; forming a vector U by the input variables of all following type equipment and network construction type equipment, forming a vector Y by all output variables, constructing transfer functions of a connecting line and a transformer in a new energy station according to the U and the Y, and finding out a dominant characteristic function of the system under a specific frequency; and fitting an interaction intensity matrix P according to each model and the function, and calculating the P, wherein the calculated result is the interaction intensity of the new energy station units under the specific frequency. The result of the invention can guide the planning department of the electric power system to plan the new energy station better so as to avoid the oscillation problem caused by strong interaction between units.
Description
Technical Field
The invention relates to the technical field of strength evaluation among new energy station units, in particular to a method and a system for evaluating interaction strength among new energy station units.
Background
As more and more power electronics are connected to the power system, conventional synchronous generator-based power systems are gradually becoming power electronics-based power systems. In recent years, accidents accompanied by broadband oscillation frequently occur, the safe operation of the system is threatened, and great challenges are brought to system operators. The rapid dynamics and high-order nonlinearity of the power electronics, as well as the complex interactions between the various power electronics, make system analysis and cognition more difficult than traditional synchronous machine-dominated power systems. Therefore, the mature experience and method based on synchronous machines in conventional power systems may no longer be applicable, and development of new methods to evaluate interactions between diverse power electronics in new energy sites is urgently needed.
Power electronics are broadly classified into two types, follower-type devices and mesh-type devices, according to different synchronization implementations. For small signal stability analysis of an electric power system, there are two main linearization methods, including time domain state space eigenvalue analysis and frequency domain analysis. The state space analysis, as a detailed system-wide description, can accurately predict stability, but it requires detailed control information. In frequency domain analysis, impedance (or admittance) models employing voltage or current as input or output are commonly used in power electronics-based power systems.
The novel power system stability analysis is not only satisfactory for small signal stability calculated by the method, but the power grid operators more hopefully know how different units interact. The invention provides a method for evaluating interaction strength between new energy station units under specific frequency. The invention aims to guide a power system planning department to better plan a new energy station so as to avoid the oscillation problem caused by strong interaction between units. Reference comments may also be provided for accident analysis.
Disclosure of Invention
The present invention has been made in view of the above-described problems.
Accordingly, the present invention aims to solve the problems: the prior art has been to evaluate interactions between units by data driven analysis. This method requires the measured data of the station. Some new energy stations do not have the conditions of experimental measurement, and even if the new energy stations have the conditions, the problem of high experimental measurement cost exists. And the results of the methods depend on the data of each experiment, and when the site control parameters and the connection method are changed, accurate results can be obtained only by re-experiment measurement, so that the cost is greatly increased.
In order to solve the technical problems, the invention provides the following technical scheme: the method and system for evaluating interaction strength between new energy station units comprise the steps of collecting input and output at a device grid-connected point, and establishing a power current model of electromechanical or electromagnetic scale of a following device and a power voltage model of electromechanical or electromagnetic scale of a network-structured device in a frequency domain; forming a vector U by the input variables of all following type equipment and network construction type equipment, forming a vector Y by all output variables, constructing transfer functions of a connecting line and a transformer in a new energy station according to the U and the Y, and finding out a dominant characteristic function of the system under a specific frequency; and fitting an interaction intensity matrix P according to each model and the function, and calculating the P, wherein the calculated result is the interaction intensity of the new energy station units under the specific frequency.
As a preferable scheme of the method for evaluating the interaction strength between the new energy station units, the method comprises the following steps: the power current model of the electromechanical or electromagnetic scale of the following type equipment comprises the steps of numbering all following type equipment in a new energy station from 1 to m in sequence, numbering all network-structured equipment in the new energy station from m+1 to n in sequence, n is the total number of the following type equipment and the network-structured equipment, taking the active power and the reactive power at the point of equipment connection as input variables, taking the amplitude and the phase of the output current at the point of equipment connection as output variables, establishing the power current model of the electromechanical or electromagnetic scale of the following type equipment in a frequency domain, and the power current model of the ith following type equipment is expressed as:
where s is denoted as laplace operator, i=1, …, m, Δ is denoted as the variable after linearization, Δp i And DeltaQ i Respectively expressed as the output active power and reactive power of the device, delta theta I,i And DeltaI i Respectively expressed as the phase and amplitude of the device output current; the arrangement into a matrix form is expressed as:
wherein G is GML,i (s) is represented as G 1-4,i (s) transfer functions respectively corresponding thereto.
As a preferable scheme of the method for evaluating the interaction strength between the new energy station units, the method comprises the following steps: the electromechanical or electromagnetic power voltage model of the network-structured equipment comprises linearizing the network-structured equipment at a working point according to the electromechanical or electromagnetic model of the network-structured equipment, taking the active power and the reactive power at the point of equipment connection as input variables, taking the amplitude and the phase of the output voltage at the point of equipment connection as output variables, and establishing the electromechanical or electromagnetic power voltage model of the network-structured equipment under the frequency domain, wherein the power voltage model of the k network-structured equipment is expressed as:
where s is denoted as laplace operator, k=m+1, …, n, Δ is denoted as the variable after linearization, Δp k And DeltaQ k Respectively indicated as the input of the deviceActive power and reactive power, delta theta t,k And DeltaU k Respectively expressed as the phase and amplitude of the device output voltage; the arrangement into a matrix form is expressed as:
wherein G is GMF,k (s) is represented as G 5-8,k (s) transfer functions respectively corresponding thereto.
As a preferable scheme of the method for evaluating the interaction strength between the new energy station units, the method comprises the following steps: the construction of the transfer function of the connecting line and the transformer in the new energy station comprises the steps of forming a vector U by the input variables of all following type equipment and network construction type equipment, forming a vector Y by all output variables, and specifically representing:
by G GML,i (s) and G GMF,k (s) is a diagonal element, and the obtained block diagonal transfer function matrix G(s) is expressed as:
G(s)＝diag{…,G GML,i (s),…,G GMF,k (s),…}
wherein G(s) is a 2 n-dimensional square matrix representing the current or voltage variation of all devices driven by active and reactive power, diag being represented as a diagonal matrix;
the vector Y is used as an input variable, the vector U is used as an output variable, and the transfer function N(s) of the connecting line of the new energy station and the transformer is constructed as follows:
U＝N(s)Y。
as a preferable scheme of the method for evaluating the interaction strength between the new energy station units, the method comprises the following steps: the dominant characteristic function comprises the steps of establishing an open-loop transfer function matrix H(s) according to N(s) and G(s) and expressing the open-loop transfer function matrix H(s) as follows according to a multivariate frequency domain theory:
H(s)＝-N(s)G(s)
for frequency omega 1 Will bejω 1 Substituting H(s), obtaining the j-th eigenvalue of H(s) to be expressed as lambda j (jω 1 ) J=1, …,2n, and calculating a corresponding right eigenvector matrix V (jω) according to matrix theory 1 ) And a left eigenvector matrix W (jω) 1 ) Satisfy H (j omega) 1 ) Expressed as:
H(jω 1 )＝V(jω 1 )diag{...,λ j (jω 1 ),...}W T (jω 1 )
for each j=1, …,2n, |λ is calculated j (jω 1 ) +1| and sorting, the feature function corresponding to the minimum value is taken as the dominant feature function, lambda is used c (jω 1 ) And (3) representing.
As a preferable scheme of the method for evaluating the interaction strength between the new energy station units, the method comprises the following steps: the interaction intensity matrix P comprises a 2 n-dimensional square matrix defined as the interaction intensity matrix P, and the a-th row and the b-th column of the P are P ab ，a，b＝1，…，2n；
When a is not equal to b, p ab Expressed as:
p ab ＝|w 2a-1,c (jω 1 )v 2b-1c (jω 1 )+w 2a,c (jω 1 )v 2b-1c (jω 1 )
+w 2a-1,c (jω 1 )v 2b,c (jω 1 )+w 2a,c (jω 1 )v 2b,c (jω 1 )|
when a=b, p ab Expressed as:
p ab ＝0
wherein w is 2a-1,c Represented as W (jω) 1 ) Lines 2a-1 and column c, v 2b-1,c Represented as V (jω) 1 ) Row 2b-1 and column c, w 2a,c Represented as W (jω) 1 ) Row 2a and column c, v 2b,c Represented as V (jω) 1 ) 2b and c.
As a preferable scheme of the method for evaluating the interaction strength between the new energy station units, the method comprises the following steps: the interaction strength includes the a-th row and b-th list of the P matrixShowing the interaction strength of the a-th device to the b-th device, p ab Is proportional to the effect of the a-th device on the b-th device.
Another object of the present invention is to provide a system for evaluating the interaction strength between units of a new energy station, which can construct an open loop transfer function of the system, find a dominant characteristic function under a frequency of interest, and finally evaluate the interaction strength between units by calculating an interaction strength matrix P.
In order to solve the technical problems, the invention provides the following technical scheme: the system of the new energy station unit interaction strength evaluation method comprises the following steps: the model building module and the interaction strength module; the model construction module is used for constructing a model and a transfer function, collecting input and output at the equipment grid-connected point, establishing a power current model of a following equipment electromechanical or electromagnetic scale and a power voltage model of a network-structured equipment electromechanical or electromagnetic scale under a frequency domain, constructing the transfer function of a connecting line and a transformer in a new energy station according to the output and the output, and finding out a dominant characteristic function of a system under a specific frequency; the interaction strength module is used for calculating interaction strength among the units, fitting an interaction strength matrix P according to each model and the function, and calculating the P, wherein the calculated result is the interaction strength among the units of the new energy station under specific frequency.
A computer device comprising a memory and a processor, said memory storing a computer program, characterized in that said processor, when executing said computer program, carries out the steps of the method for evaluating the strength of interaction between sets of new energy stations as described above.
A computer-readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method for evaluating the interaction strength between sets of new energy stations as described above.
The invention has the beneficial effects that: the invention provides an evaluation method for interaction strength between new energy station units formed by following type and network-structured power electronic equipment under specific frequency. According to the invention, the dominant characteristic function under the concerned frequency is found by constructing the open loop transfer function of the system, and finally the mutual strength among the units is estimated by calculating the interaction strength matrix P. The method can obtain the interaction strength among the units under different frequency bands by implementing the method for multiple times, and the result of the method can guide the power system planning department to better plan the new energy station so as to avoid the oscillation problem caused by strong interaction among the units. Reference comments may also be provided for accident analysis.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art. Wherein:
fig. 1 is a flowchart of a method for evaluating interaction strength between new energy station sets in embodiment 1.
FIG. 2 is a graph of a transfer function matrix model of a method for evaluating interaction strength between sets of new energy field stations in example 1
Fig. 3 is a block diagram of an evaluation system for interaction strength between new energy station units in example 3.
Detailed Description
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 specific embodiments thereof which are illustrated in the appended drawings.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, but the present invention may be practiced in other ways other than those described herein, and persons skilled in the art will readily appreciate that the present invention is not limited to the specific embodiments disclosed below.
Example 1
Referring to fig. 1 and 2, a first embodiment of the present invention provides a method for evaluating interaction strength between sets of new energy stations, which includes the steps of, as shown in fig. 1:
step 1: and collecting input and output of the equipment at the grid-connected point, and establishing a power current model of electromechanical or electromagnetic scale of the following equipment and a power voltage model of electromechanical or electromagnetic scale of the net-structured equipment in a frequency domain.
The photovoltaic equipment, the wind power equipment, the energy storage equipment and the Statcom device based on phase-locked loop control in the new energy station are uniformly considered as follow-up equipment, all the follow-up equipment in the new energy station are numbered from 1 to m in sequence, all the net-structured equipment in the new energy station are numbered from m+1 to n in sequence, and n is the total number of the follow-up equipment and the net-structured equipment. The following type equipment is generally equivalent to a controlled current source, is linearized at a working point according to an electromechanical model or an electromagnetic model of the networking type equipment, takes active power and reactive power at the equipment networking point as input variables, takes the amplitude and the phase of output current at the equipment networking point as output variables, establishes a power current model of electromechanical or electromagnetic dimensions of the following type equipment in a frequency domain, and the power current model of the i-th following type equipment is expressed as:
where s is denoted as laplace operator, i=1, …, m, Δ is denoted as the variable after linearization, Δp i And DeltaQ i Respectively expressed as the output active power and reactive power of the device, delta theta I,i And DeltaI i Respectively expressed as the phase and amplitude of the device output current; the arrangement into a matrix form is expressed as:
wherein G is GML,i (s) is represented as G 1-4,i (s) transfer functions respectively corresponding to
A mesh-type device is generally equivalent to a controlled voltage source. Linearizing the network-structured equipment at a working point according to an electromechanical model or an electromagnetic model of the network-structured equipment, taking active power and reactive power at a network-connected point of the equipment as input variables, taking amplitude and phase of output voltage at the network-connected point of the equipment as output variables, and establishing a power voltage model of the network-structured equipment in electromechanical or electromagnetic scale under a frequency domain, wherein the power voltage model of the k-th network-structured equipment is expressed as:
where s is denoted as laplace operator, k=m+1, …, n, Δ is denoted as the variable after linearization, Δp k And DeltaQ k Respectively expressed as the output active power and reactive power of the device, delta theta t,k And DeltaU k Respectively expressed as the phase and amplitude of the device output voltage; the arrangement into a matrix form is expressed as:
wherein G is GMF,k (s) is represented as G 5-8,k (s) transfer functions respectively corresponding thereto.
Step 2: and forming a vector U by the input variables of all the following type equipment and the network construction type equipment, forming a vector Y by the output variables, constructing transfer functions of a connecting line and a transformer in the new energy station according to the U and the Y, and finding out a dominant characteristic function of the system under a specific frequency.
Besides the following type equipment and the network construction type equipment, the new energy station also comprises a connecting line and a transformer, wherein the input variables of all the following type equipment and the network construction type equipment form a vector U, and the output variables form a vector Y which is specifically expressed as:
by G GML,i (s) and G GMF,k (s) is a diagonal member, resulting in a block diagonal transfer functionThe matrix G(s) is expressed as:
G(s)＝diag{…,G GML,i (s),…,G GMF,k (s),…}
where G(s) is a 2 n-dimensional square matrix representing the change in current or voltage of all devices driven by active and reactive power, diag being represented as a diagonal matrix.
The vector Y is used as an input variable, the vector U is used as an output variable, and the transfer function N(s) of the connecting line of the new energy station and the transformer is constructed as follows:
U＝N(s)Y。
as shown in fig. 2, according to the multivariate frequency domain theory, establishing an open loop transfer function matrix H(s) from N(s) and G(s) is expressed as:
H(s)＝-N(s)G(s)
for frequency omega 1 Will j omega 1 Substituting H(s), obtaining the j-th eigenvalue of H(s) to be expressed as lambda j (jω 1 ) J=1, …,2n, and calculating a corresponding right eigenvector matrix V (jω) according to matrix theory 1 ) And a left eigenvector matrix W (jω) 1 ) Satisfy H (j omega) 1 ) Expressed as:
H(jω 1 )＝V(jω 1 )diag{...,λ j (jω 1 ),...}W T (jω 1 )
for each j=1, …,2n, |λ is calculated j (jω 1 ) +1| and sorting, the feature function corresponding to the minimum value is taken as the dominant feature function, lambda is used c (jω 1 ) And (3) representing.
Step 3: and fitting an interaction intensity matrix P according to each model and the function, and calculating the P, wherein the calculated result is the interaction intensity of the new energy station units under the specific frequency.
Defining the interaction strength matrix P as a 2 n-dimensional square matrix, and the a-th row and the b-th column of the P as P ab A, b=1, …,2n; when a is not equal to b, p ab Expressed as:
p ab ＝|w 2a-1,c (jω 1 )v 2b-1c (jω 1 )+w 2a,c (jω 1 )v 2b-1c (jω 1 )
+w 2a-1,c (jω 1 )v 2b,c (jω 1 )+w 2a,c (jω 1 )v 2b,c (jω 1 )|
when a=b, p ab Expressed as:
p ab ＝0
wherein w is 2a-1,c Represented as W (jω) 1 ) Lines 2a-1 and column c, v 2b-1,c Represented as V (jω) 1 ) Row 2b-1 and column c, w 2a,c Represented as W (jω) 1 ) Row 2a and column c, v 2b,c Represented as V (jω) 1 ) 2b and c.
Row a and column b of the P matrix show the strength of interaction of the a-th device with the b-th device (including the following type device and the net-structured device), P ab The larger the device, the more obvious the effect of the device on the b.
Example 2
A second embodiment of the present invention, which is different from the first embodiment, is: the method for evaluating the interaction strength between the new energy station units further comprises the step of comparing test results by means of scientific demonstration by adopting the traditional technical scheme with the release of the method in order to verify and explain the technical effects adopted in the method, so as to verify the true effects of the method.
Table 1: the method of the invention is compared with the traditional method in small signal stability analysis under different scenes
As can be seen from the table above, the inventive method exhibits a more negative real part in all scenarios, which suggests that the system is more prone to return to steady state with better stability.
The method provides quantitative measurement for the interaction between each unit in the system, so that the interaction of power electronic equipment is processed more effectively, and the impedance value is lower, so that the system has lower overall resistance.
The method of the invention has faster response time, which means that the system can reach a new steady state more quickly when disturbance occurs.
Example 3
Referring to fig. 3, a third embodiment of the present invention is shown, which is different from the first two embodiments: a system of a new energy station unit interaction strength evaluation method comprises a model construction module and an interaction strength module; the model construction module is used for constructing a model and a transfer function, collecting input and output at the grid-connected point of the equipment, establishing a power current model of electromechanical or electromagnetic scale of the following equipment and a power voltage model of electromechanical or electromagnetic scale of the grid-structured equipment in a frequency domain, constructing the transfer function of a connecting line and a transformer in the new energy station according to the output and the output, and finding out a dominant characteristic function of the system under specific frequency; the interaction strength module is used for calculating interaction strength among the units, fitting an interaction strength matrix P according to each model and the function, and calculating the P, wherein the calculated result is the interaction strength among the units of the new energy station under the specific frequency.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-only memory (ROM), a random access memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
Logic and/or steps represented in the flowcharts or otherwise described herein, e.g., a ordered listing of executable instructions for implementing logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device.
More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or more wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). In addition, the computer readable medium may even be paper or other suitable medium on which the program is printed, as the program may be electronically captured, via, for instance, optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner, if necessary, and then stored in a computer memory.
It is to be understood that portions of the present invention may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as is well known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
It should be noted that the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that the technical solution of the present invention may be modified or substituted without departing from the spirit and scope of the technical solution of the present invention, which is intended to be covered in the scope of the claims of the present invention.
Claims (10)
1. The method for evaluating the interaction strength between the new energy station units is characterized by comprising the following steps of: comprising the steps of (a) a step of,
collecting input and output of equipment grid connection points, and establishing a power current model of electromechanical or electromagnetic scale of the following equipment and a power voltage model of electromechanical or electromagnetic scale of the net-structured equipment in a frequency domain;
forming a vector U by the input variables of all following type equipment and network construction type equipment, forming a vector Y by all output variables, constructing transfer functions of a connecting line and a transformer in a new energy station according to the U and the Y, and finding out a dominant characteristic function of the system under a specific frequency;
and fitting an interaction intensity matrix P according to each model and the function, and calculating the P, wherein the calculated result is the interaction intensity of the new energy station units under the specific frequency.
2. The method for evaluating interaction strength between new energy station units according to claim 1, wherein: the power current model of the electromechanical or electromagnetic scale of the following type equipment comprises the steps of numbering all following type equipment in a new energy station from 1 to m in sequence, numbering all network-structured equipment in the new energy station from m+1 to n in sequence, n is the total number of the following type equipment and the network-structured equipment, taking the active power and the reactive power at the point of equipment connection as input variables, taking the amplitude and the phase of the output current at the point of equipment connection as output variables, establishing the power current model of the electromechanical or electromagnetic scale of the following type equipment in a frequency domain, and the power current model of the ith following type equipment is expressed as:
where s is denoted as laplace operator, i=1, …, m, Δ is denoted as the variable after linearization, Δp i And DeltaQ i Respectively expressed as the output active power and reactive power of the device, delta theta I,i And DeltaI i Respectively expressed as the phase and amplitude of the device output current;
the arrangement into a matrix form is expressed as:
wherein G is GML,i (s) is represented as G 1-4,i (s) transfer functions respectively corresponding thereto.
3. The method for evaluating interaction strength between new energy station units according to claim 2, wherein: the electromechanical or electromagnetic power voltage model of the network-structured equipment comprises linearizing the network-structured equipment at a working point according to the electromechanical or electromagnetic model of the network-structured equipment, taking the active power and the reactive power at the point of equipment connection as input variables, taking the amplitude and the phase of the output voltage at the point of equipment connection as output variables, and establishing the electromechanical or electromagnetic power voltage model of the network-structured equipment under the frequency domain, wherein the power voltage model of the k network-structured equipment is expressed as:
where s is denoted as laplace operator, k=m+1, …, n, Δ is denoted as the variable after linearization, Δp k And DeltaQ k Respectively expressed as the output active power and reactive power of the device, delta theta t,k And DeltaU k Respectively expressed as the phase and amplitude of the device output voltage;
the arrangement into a matrix form is expressed as:
wherein G is GMF,k (s) is represented as G 5-8,k (s) transfer functions respectively corresponding thereto.
4. The method for evaluating interaction strength between new energy station units according to claim 3, wherein: the construction of the transfer function of the connecting line and the transformer in the new energy station comprises the steps of forming a vector U by the input variables of all following type equipment and network construction type equipment, forming a vector Y by all output variables, and specifically representing:
by G GML,i (s) and G GMF,k (s) is a diagonal element, and the obtained block diagonal transfer function matrix G(s) is expressed as:
G(s)＝diag{…,G GML,i (s),…,G GMF,k (s),…}
wherein G(s) is a 2 n-dimensional square matrix representing the current or voltage variation of all devices driven by active and reactive power, diag being represented as a diagonal matrix;
the vector Y is used as an input variable, the vector U is used as an output variable, and the transfer function N(s) of the connecting line of the new energy station and the transformer is constructed as follows:
U＝N(s)Y。
5. the method for evaluating interaction strength between new energy station units according to claim 4, wherein: the dominant characteristic function comprises the steps of establishing an open-loop transfer function matrix H(s) according to N(s) and G(s) and expressing the open-loop transfer function matrix H(s) as follows according to a multivariate frequency domain theory:
H(s)＝-N(s)G(s)
for frequency omega 1 Will j omega 1 Substituting H(s), obtaining the j-th eigenvalue of H(s) to be expressed as lambda j (jω 1 ) J=1, …,2n, and calculating a corresponding right eigenvector matrix V (jω) according to matrix theory 1 ) And a left eigenvector matrix W (jω) 1 ) Satisfy H (j omega) 1 ) Expressed as:
H(jω 1 )＝V(jω 1 )diag{...,λ j (jω 1 ),...}W T (jω 1 )
for each j=1, …,2n, |λ is calculated j (jω 1 ) +1| and sorting, the feature function corresponding to the minimum value is taken as the dominant feature function, lambda is used c (jω 1 ) And (3) representing.
6. The method for evaluating interaction strength between new energy station units according to claim 5, wherein: the interaction intensity matrix P comprises a 2 n-dimensional square matrix defined as the interaction intensity matrix P, and the a-th row and the b-th column of the P are P ab ，a，b＝1，…，2n；
When a is not equal to b, p ab Expressed as:
p ab ＝|w 2a-1,c (jω 1 )v 2b-1c (jω 1 )+w 2a,c (jω 1 )v 2b-1c (jω 1 )+w 2a-1,c (jω 1 )v 2b,c (jω 1 )+w 2a,c (jω 1 )v 2b,c (jω 1 )|
when a=b, p ab Expressed as:
p ab ＝0
wherein w is 2a-1,c Represented as W (jω) 1 ) Lines 2a-1 and column c, v 2b-1,c Represented as V (jω) 1 ) Row 2b-1 and column c, w 2a,c Represented as W (jω) 1 ) Row 2a and column c, v 2b,c Represented as V (jω) 1 ) 2b and c.
7. The new energy station set of claim 6A method for evaluating the strength of interaction between two substrates, characterized by: the interaction strength comprises that the a-th row and the b-th column of the P matrix represent the interaction strength of the a-th device to the b-th device, and P ab Is proportional to the effect of the a-th device on the b-th device.
8. A system employing the method for evaluating interaction strength between new energy station sets according to any one of claims 1 to 7, characterized in that: the system comprises a model construction module and an interaction strength module;
the model construction module is used for constructing a model and a transfer function, collecting input and output at the equipment grid-connected point, establishing a power current model of a following equipment electromechanical or electromagnetic scale and a power voltage model of a network-structured equipment electromechanical or electromagnetic scale under a frequency domain, constructing the transfer function of a connecting line and a transformer in a new energy station according to the output and the output, and finding out a dominant characteristic function of a system under a specific frequency;
the interaction strength module is used for calculating interaction strength among the units, fitting an interaction strength matrix P according to each model and the function, and calculating the P, wherein the calculated result is the interaction strength among the units of the new energy station under specific frequency.
9. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that: the processor, when executing the computer program, implements the steps of the method for evaluating the interaction strength between the new energy station sets of any one of claims 1 to 7.
10. A computer-readable storage medium having stored thereon a computer program, characterized by: the computer program, when executed by a processor, implements the steps of the method for evaluating the interaction strength between sets of new energy field stations of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311372695.1A CN117638867A (en) | 2023-10-23 | 2023-10-23 | Method and system for evaluating interaction strength between new energy station units |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311372695.1A CN117638867A (en) | 2023-10-23 | 2023-10-23 | Method and system for evaluating interaction strength between new energy station units |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117638867A true CN117638867A (en) | 2024-03-01 |
Family
ID=90020654
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311372695.1A Pending CN117638867A (en) | 2023-10-23 | 2023-10-23 | Method and system for evaluating interaction strength between new energy station units |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117638867A (en) |
-
2023
- 2023-10-23 CN CN202311372695.1A patent/CN117638867A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ma et al. | Partially‐coupled gradient‐based iterative algorithms for multivariable output‐error‐like systems with autoregressive moving average noises | |
| Abdallah et al. | Stator winding inter‐turn short‐circuit detection in induction motors by parameter identification | |
| Yu et al. | Frequency domain modelling and control of fractional‐order system for permanent magnet synchronous motor velocity servo system | |
| Zhang et al. | Artificial neural network based identification of multi-operating-point impedance model | |
| Rojas-Dueñas et al. | Black-box modelling of a dc-dc buck converter based on a recurrent neural network | |
| CN112330488B (en) | Power grid frequency situation prediction method based on transfer learning | |
| CN115238980A (en) | Micro-grid transient power angle trajectory prediction method based on attention mechanism and space-time diagram convolutional neural network | |
| Picardi et al. | Parameter identification of induction motor based on particle swarm optimization | |
| Norrlöf | Disturbance rejection using an ILC algorithm with iteration varying filters | |
| Pasqualotto et al. | Increasing feasibility of neural network-based early fault detection in induction motor drives | |
| CN114583767A (en) | Data-driven wind power plant frequency modulation response characteristic modeling method and system | |
| CN117638867A (en) | Method and system for evaluating interaction strength between new energy station units | |
| Ashar et al. | ARX model identification for the real-time temperature process with Matlab-arduino implementation | |
| Xiao et al. | A novel impedance measurement technique for power electronic systems | |
| CN109613347B (en) | Excitation characteristic detection method and device of transformer | |
| Peter et al. | Voltage stability assessment in power systems using Artificial Neural Networks | |
| CN107730002A (en) | A kind of communication network shutdown remote control parameter intelligent fuzzy comparison method | |
| Fetzer | Development of a MATLAB/Simulink Framework for Phasor-Based Power System Simulation and Component Modeling Based on State Machines | |
| Romero Ugalde et al. | Robust estimation of balanced simplicity-accuracy neural networks-based models | |
| CN113221064A (en) | Model processing method, model processing device, computer equipment and storage medium | |
| Meher et al. | Relay‐based identification of Wiener model | |
| Barocio et al. | Analysis of factors affecting power system normal form results | |
| Marchiori et al. | Development and validation of an electromagnetic model of the active control system of MHD modes in RFX-mod | |
| Birbir et al. | Prediction of current harmonics in induction motors with Artificial Neural Network | |
| Zaman et al. | Critical equivalent series resistance estimation for voltage regulator stability using hybrid system identification and neural network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |