New energy station impedance network modeling method and device, electronic equipment and medium

Abstract

The invention discloses a new energy station impedance network modeling method and device, electronic equipment and a medium, which are used for solving the problem that the currently adopted frequency domain impedance modeling method mainly aims at single converter modeling and cannot be directly applied to a station-level system comprising multiple units and complex topology. Wherein the new energy station comprises a plurality of units; obtaining station line data of a new energy station and filtering parameter data of each unit; for each unit, constructing a single-unit admittance matrix according to the filtering parameter data in combination with dynamic link correction; constructing a coupling admittance matrix according to the station line data, and constructing a total admittance matrix according to the coupling admittance matrix and each single-machine admittance matrix; and based on the total admittance matrix, obtaining equivalent input impedance for evaluating the interaction stability of the new energy station and the power grid side impedance through matrix inversion.

Description

Technical Field

[0001] The present invention relates to the technical field of power system simulation, and in particular to a new energy station impedance network modeling method, device, electronic equipment and medium. Background Art

[0002] With the large-scale integration of renewable energy sources such as wind power and photovoltaics into the power grid, renewable energy sites (such as wind farms and photovoltaic power plants) have become an increasingly important component of the power system. However, these renewable energy devices are generally connected to the grid through power electronic converters. These converters exhibit complex control dynamics and diverse frequency domain characteristics, leading to new stability issues in the system, such as subsynchronous and supersynchronous oscillations and broadband resonance.

[0003] Traditional power system stability analysis methods are unable to effectively address high-frequency dynamic problems. Therefore, it is urgent to establish an accurate frequency-domain impedance network model to provide a theoretical tool for stability analysis and resonance suppression.

[0004] Frequency-domain impedance models, by characterizing the system's voltage-current response at specific frequencies, can intuitively reveal key information such as the location of the resonance point and stability boundary. Current frequency-domain impedance modeling methods are primarily designed for single-unit converter modeling and cannot be directly applied to station-level systems with multiple units and complex topologies. Summary of the Invention

[0005] The present invention provides a method, device, electronic device and medium for modeling the impedance network of a new energy station, which are used to solve or partially solve the problem that the current frequency domain impedance modeling method is mainly used for modeling a single converter and cannot be directly used for station-level systems containing multiple units and complex topologies.

[0006] The present invention provides a method for modeling an impedance network of a new energy station, wherein the new energy station includes multiple units; the method comprises:

[0007] Obtaining station line data of the new energy station and filter parameter data of each unit;

[0008] For each of the units, a single-unit admittance matrix is ​​constructed based on the filtering parameter data and combined with dynamic link correction;

[0009] Constructing a coupling admittance matrix according to the station line data, and constructing a total admittance matrix according to the coupling admittance matrix and each of the single-machine admittance matrices;

[0010] Based on the total admittance matrix, an equivalent input impedance is obtained by matrix inversion, and the equivalent input impedance is used to evaluate the interactive stability between the new energy station and the grid-side impedance.

[0011] Optionally, the filtering parameter data includes filtering inductance, filtering inductance equivalent resistance, and grid equivalent impedance; and constructing a single-machine admittance matrix based on the filtering parameter data and combined with dynamic link correction includes:

[0012] Based on the filter inductance and the grid equivalent impedance, a phase-locked loop correction is performed in combination with a frequency domain transfer function to obtain the dq-axis cross-coupling impedance and the qd-axis cross-coupling impedance of the converter in the unit;

[0013] Based on the filter inductor and the filter inductor equivalent resistance, the current control loop is corrected in combination with the current loop closed-loop transfer function to obtain the d-axis self-impedance and q-axis self-impedance of the converter in the unit;

[0014] constructing a frequency domain impedance matrix of the unit according to the dq-axis cross-coupling impedance, the qd-axis cross-coupling impedance, the d-axis self-impedance, and the q-axis self-impedance;

[0015] Based on the frequency domain impedance matrix, the single-machine admittance matrix of the unit is obtained by matrix inversion.

[0016] Optionally, the phase-locked loop correction process is performed by the following formula:

[0017]

[0018]

[0019] in, represents the dq axis cross-coupling impedance; represents the qd-axis cross-coupling impedance; represents the angular frequency; represents an imaginary unit; Indicates the DC bus voltage; represents the frequency domain transfer function; Represents the filter inductor; 、 PI controller parameters representing the frequency domain transfer function; Represents the equivalent impedance of the power grid.

[0020] Optionally, the current control loop correction process is performed by the following formula:

[0021]

[0022]

[0023]

[0024] in, represents the d-axis self-impedance; represents the q-axis self-impedance; represents the angular frequency; represents an imaginary unit; Represents the filter inductor; Indicates the equivalent resistance of the filter inductor; represents the current loop closed-loop transfer function; =1.5 Indicates the total delay time; represents the switching cycle; represents the frequency domain response of the PI controller in the current loop closed-loop transfer function; 、 PI controller parameters representing the current loop closed-loop transfer function.

[0025] Optionally, the station line data includes a station electrical wiring diagram, position coordinates of each unit, and line parameters of each line; and constructing a coupling admittance matrix based on the station line data includes:

[0026] Constructing a station topology connection matrix representing the line connection relationship in the new energy station based on the station electrical wiring diagram and the position coordinates of each unit;

[0027] For each of the lines, constructing a line impedance according to the line parameters;

[0028] A coupling admittance matrix characterizing multi-machine coupling in the new energy station is constructed based on the station topology connection matrix and each of the line impedances.

[0029] Optionally, constructing a total admittance matrix according to the coupled admittance matrix and each of the stand-alone admittance matrices includes:

[0030] Based on each of the single-machine admittance matrices, a block diagonal matrix is ​​constructed;

[0031] A total admittance matrix is ​​constructed by summing the block diagonal matrix and the coupled admittance matrix.

[0032] Optionally, the method further includes:

[0033] Obtaining a measured impedance, and dynamically adjusting model parameters of each of the stand-alone admittance matrices according to the measured impedance and the equivalent input impedance and in combination with an adaptive learning rate matrix;

[0034] The model parameters include filter inductance, filter inductance equivalent resistance, PI controller parameters of frequency domain transfer function, and PI controller parameters of current loop closed-loop transfer function.

[0035] The present invention also provides an impedance network modeling device for a new energy station, wherein the new energy station includes multiple units; the device comprises:

[0036] A data acquisition unit, configured to acquire station line data of the new energy station and filter parameter data of each unit;

[0037] A dynamic link correction unit is used to construct a single-machine admittance matrix for each of the units based on the filtering parameter data and combined with dynamic link correction;

[0038] An admittance matrix construction unit, configured to construct a coupling admittance matrix according to the station line data, and to construct a total admittance matrix according to the coupling admittance matrix and each of the stand-alone admittance matrices;

[0039] An equivalent input impedance solving unit is used to obtain an equivalent input impedance by matrix inversion based on the total admittance matrix, and the equivalent input impedance is used to evaluate the interactive stability between the new energy station and the grid-side impedance.

[0040] The present invention further provides an electronic device, comprising a processor and a memory:

[0041] The memory is used to store program code and transmit the program code to the processor;

[0042] The processor is used to execute the new energy station impedance network modeling method as described in any one of the above items according to the instructions in the program code.

[0043] The present invention also provides a computer-readable storage medium, which is used to store program code, and the program code is used to execute the new energy station impedance network modeling method as described in any one of the above items.

[0044] It can be seen from the above technical solutions that the present invention has the following advantages:

[0045] A method for modeling an impedance network for a new energy station is provided. The new energy station includes multiple generators. Station line data and filter parameter data for each generator are obtained. For each generator, a single-unit admittance matrix is ​​constructed based on the filter parameter data and dynamic link correction. A coupling admittance matrix is ​​constructed based on the station line data, and a total admittance matrix is ​​constructed based on the coupling admittance matrix and the individual generator admittance matrices. Based on the total admittance matrix, an equivalent input impedance for evaluating the interaction stability between the new energy station and the grid-side impedance is obtained through matrix inversion. For a new energy station containing wind power and photovoltaic converter clusters (for ease of explanation, the present invention uses generator units as the representation), this method uses dynamic link correction to progressively design single-unit impedance modeling, multi-unit coupling modeling, and station-level impedance aggregation. This method uses a multidimensional frequency-domain decoupling algorithm to establish a full-band impedance network model covering 0.1 Hz to 2 kHz. This model accurately characterizes the station's impedance characteristics across a wide frequency range (0.1 Hz to 2 kHz), meeting the requirements for collaborative analysis of subsynchronous / supersynchronous oscillations and high-frequency resonances. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0047] Figure 1 A flowchart of the steps of a new energy station impedance network modeling method;

[0048] Figure 2 This is a schematic diagram of the overall process of a new energy station impedance network modeling method;

[0049] Figure 3 This is a structural block diagram of a new energy station impedance network modeling device. DETAILED DESCRIPTION

[0050] Embodiments of the present invention provide a method, device, electronic device, and medium for modeling an impedance network of a new energy station, which are used to solve or partially solve the problem that the currently used frequency domain impedance modeling method is mainly used for modeling a single converter and cannot be directly used for station-level systems containing multiple units and complex topologies.

[0051] In order to make the purpose, features, and advantages of the present invention more obvious and easy to understand, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the embodiments described below are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0052] As an example, frequency-domain impedance models can intuitively reveal key information such as the location of resonance points and stability boundaries by characterizing the system's voltage-current response at specific frequencies. Current frequency-domain impedance modeling methods are primarily designed for single-unit converter modeling and cannot be directly applied to station-level systems with multiple units and complex topologies.

[0053] Therefore, one of the core inventive aspects of the present invention is to propose a new energy station impedance network modeling method based on dynamic coupling modeling to quantify the impact of inter-unit electrical connections and control interactions and achieve precise integration of station-level impedance networks. For new energy stations containing wind and photovoltaic converter clusters (for ease of illustration, this invention uses units as the representation), this method, combined with dynamic link correction, uses a layered approach to design single-unit impedance modeling, multi-unit coupling modeling, and station-level impedance aggregation. This method uses a multidimensional frequency-domain decoupling algorithm to establish a full-band impedance network model covering 0.1 Hz to 2 kHz. This accurately characterizes the station's impedance characteristics across a wide frequency range (0.1 Hz to 2 kHz), meeting the requirements for collaborative analysis of subsynchronous / supersynchronous oscillations and high-frequency resonances. Furthermore, based on this constructed full-band impedance network model, an online parameter correction algorithm is provided to dynamically update model parameters based on real-time operating data to maintain model accuracy under varying operating conditions, thereby supporting dynamic impedance characteristic analysis under multiple operating conditions.

[0054] Reference Figure 1 , shows a flowchart of the steps of a method for modeling an impedance network of a new energy station provided by an embodiment of the present invention. The new energy station includes multiple units; the method may specifically include the following steps:

[0055] Step 101: Acquire the station line data of the new energy station and the filter parameter data of each unit;

[0056] Among them, the station line data mainly includes the station electrical wiring diagram, the location coordinates of each unit, and the line parameters of each line (such as the resistance of the line). ,inductance ,capacitance The filter parameter data mainly includes the filter inductance corresponding to a single unit. , filter inductor equivalent resistance and grid equivalent impedance Thus, the station line data and the filter parameter data of each unit are collected for subsequent impedance network modeling.

[0057] Step 102: for each of the units, construct a single-unit admittance matrix based on the filtering parameter data and combined with dynamic link correction;

[0058] This step primarily implements impedance modeling for each unit. The modeling principle is to consider key dynamic links, such as the converter's phase-locked loop (PLL), current control loop, and filter, and then establish a frequency-domain impedance matrix in the dq-axis coordinate system based on these dynamic link corrections.

[0059] Among them, dynamic link correction can mainly include phase-locked loop correction and current control loop correction.

[0060] In some embodiments, constructing a single-machine admittance matrix based on the filtering parameter data and combined with dynamic link correction may mainly include the following sub-steps S01 to S04:

[0061] Step S01: Based on the filter inductance and the grid equivalent impedance, the phase-locked loop is corrected in combination with the frequency domain transfer function to obtain the dq-axis cross-coupling impedance and the qd-axis cross-coupling impedance of the converter in the unit;

[0062] Among them, the phase tracking error of the phase-locked loop is transferred through the frequency domain transfer function Corrects cross-impedance terms in the frequency-domain impedance matrix.

[0063] Furthermore, the phase-locked loop correction process can be performed by the following formula:

[0064]

[0065]

[0066] in, represents the dq axis cross-coupling impedance; represents the qd-axis cross-coupling impedance; and The unit is Ω, which represents the asymmetric characteristics introduced by the phase-locked loop (PLL) dynamic and decoupling control. represents the angular frequency; represents an imaginary unit; Indicates the DC bus voltage in V, reflecting the energy transmission capability of the DC side of the converter; represents the frequency domain transfer function; Indicates the filter inductance, in H, used to suppress high-frequency switching harmonics; 、 PI controller parameters representing the frequency domain transfer function; Represents the equivalent impedance of the power grid.

[0067] Step S02: Based on the filter inductance and the filter inductance equivalent resistance, the current control loop is corrected in combination with the current loop closed-loop transfer function to obtain the d-axis self-impedance and q-axis self-impedance of the converter in the unit;

[0068] Furthermore, the current control loop correction process can be performed by the following formula:

[0069]

[0070]

[0071]

[0072] in, represents the d-axis self-impedance; represents the q-axis self-impedance; and The unit is Ω, which represents the response of the converter to the coaxial disturbance in the synchronous rotating coordinate system; Indicates the equivalent resistance of the filter inductor, in Ω, which represents the line loss; represents the current loop closed-loop transfer function; =1.5 Indicates the total delay time; represents the switching cycle; and The unit of is s; represents the frequency domain response of the PI controller in the current loop closed-loop transfer function; 、 PI controller parameters representing the current loop closed-loop transfer function.

[0073] Step S03: constructing a frequency domain impedance matrix of the unit according to the dq-axis cross-coupling impedance, the qd-axis cross-coupling impedance, the d-axis self-impedance, and the q-axis self-impedance;

[0074] Specifically, the unit frequency domain impedance matrix constructed in step S03 is As shown in the following formula:

[0075]

[0076] To simplify the representation, It can also be written as .

[0077] Step S04: Based on the frequency domain impedance matrix, the single-machine admittance matrix of the unit is obtained by matrix inversion.

[0078] The following formula can be used to solve the Single unit admittance matrix of a unit :

[0079]

[0080] Among them, the single machine admittance matrix The unit is S.

[0081] Step 103: constructing a coupling admittance matrix according to the station line data, and constructing a total admittance matrix according to the coupling admittance matrix and each of the stand-alone admittance matrices;

[0082] This step mainly realizes the multi-machine coupling modeling of new energy stations.

[0083] In some embodiments, constructing a coupling admittance matrix based on station line data can mainly include: first, constructing a station topology connection matrix that characterizes the line connection relationship in the new energy station based on the station electrical wiring diagram and the position coordinates of each unit; then, for each line, constructing the line impedance based on the line parameters; finally, constructing a coupling admittance matrix that characterizes the coupling of multiple machines in the new energy station based on the station topology connection matrix and the impedances of each line.

[0084] Specifically, through the station electrical wiring diagram, unit location coordinates, line parameters (line resistance ,inductance ,capacitance Parameters such as ), the site topology connection matrix can be obtained Pass the exam Line impedance of the line .

[0085] Assume that the new energy station includes Units, combined with the site topology connection matrix and the impedance of each line, the coupling admittance matrix can be constructed as shown below :

[0086]

[0087] Where, is the station topology connection matrix (i.e. line connection relationship matrix), whose dimension is ; Chinese elements When it is 1, it means Lines connecting units With the crew , otherwise 0. is the total number of lines within the station.

[0088] Furthermore, a total admittance matrix is ​​constructed based on the coupled admittance matrix and each stand-alone admittance matrix. Specifically, a block diagonal matrix is ​​constructed based on each stand-alone admittance matrix; and a total admittance matrix is ​​constructed by summing the block diagonal matrix and the coupled admittance matrix.

[0089] To sum up, the total admittance matrix It is the sum of the single unit admittance matrix and the coupling admittance matrix. The specific formula is as follows:

[0090]

[0091] Step 104: Based on the total admittance matrix, obtain an equivalent input impedance by matrix inversion. The equivalent input impedance is used to evaluate the interaction stability between the new energy station and the grid-side impedance.

[0092] Equivalent input impedance of the station grid connection point It can be obtained by matrix inversion as shown below:

[0093]

[0094] in, Used to evaluate the impedance between new energy stations and the grid side To simplify the representation, It can also be written as .

[0095] In some embodiments, after the equivalent input impedance of the station grid connection point is calculated, the model parameters involved in each calculation model constructed in the above embodiments can be dynamically adjusted based on the equivalent input impedance and the impedance data obtained by actual measurement.

[0096] In a specific implementation, the measured impedance is obtained, and the model parameters of each single-machine admittance matrix are dynamically adjusted according to the measured impedance and the equivalent input impedance, while combining the adaptive learning rate matrix.

[0097] Among them, the model parameters include filter inductance , filter inductor equivalent resistance , PI controller parameters of frequency domain transfer function ( 、 ), PI controller parameters of the current loop closed-loop transfer function ( 、 ).

[0098] Furthermore, the following formula is used to dynamically adjust the model parameters based on real-time operation data:

[0099]

[0100] Where, represents the model parameter vector (e.g. 、 、 、 、 、 ), subscript Represents the adjusted updated model parameters, Indicates the model parameters before adjustment and update; is the adaptive learning rate matrix, optimized by the Kalman filter algorithm; Indicates the impedance data measured by the swept frequency method.

[0101] Therefore, through the above updating means, it is possible to achieve the following results based on the measured impedance data: Equivalent input impedance of the station grid connection point , get the adjusted and updated model parameters , so that in the subsequent calculation process, more accurate calculation accuracy can be maintained even under changing working conditions.

[0102] In an embodiment of the present invention, a method for modeling an impedance network for a new energy station based on dynamic coupling modeling is proposed. For new energy stations containing wind power and photovoltaic converter clusters (for ease of explanation, the present invention uses units as representation), on the one hand, combined with dynamic link correction, single-unit impedance modeling, multi-unit coupling modeling, and station-level impedance aggregation are designed in a layered and progressive manner. This achieves a full-band impedance network model covering 0.1Hz to 2kHz through a multi-dimensional frequency domain decoupling algorithm. This accurately characterizes the impedance characteristics of the station over a wide frequency band (0.1Hz to 2kHz), meeting the requirements for collaborative analysis of subsynchronous / supersynchronous oscillations and high-frequency resonances. On the other hand, based on the constructed full-band impedance network model, an online parameter correction algorithm is provided to dynamically update model parameters based on real-time operating data to maintain the model's accuracy under varying operating conditions, thereby supporting dynamic impedance characteristic analysis under multiple operating conditions.

[0103] For better explanation, refer to Figure 2 , showing a schematic diagram of the overall process of a method for modeling an impedance network of a new energy station provided by an embodiment of the present invention. It should be noted that this embodiment only briefly describes the general process of modeling an impedance network of a new energy station. The specific implementation process of each step can be understood by referring to the relevant content in the aforementioned embodiments. A detailed description is omitted here. It is understood that the present invention is not limited to this.

[0104] Step 201: Obtaining station line data of a new energy station and filter parameter data of each unit;

[0105] Step 202: For each unit, based on the filter parameter data, combined with the phase-locked loop correction and the current control loop correction, the frequency domain impedance matrix of the unit is constructed, and then the single unit admittance matrix is ​​obtained by matrix inversion;

[0106] Step 203: constructing a coupling admittance matrix representing multi-machine coupling in a new energy station based on the station line data;

[0107] Step 204: construct a block diagonal matrix based on each stand-alone admittance matrix, and construct a total admittance matrix by summing the block diagonal matrix and the coupled admittance matrix;

[0108] Step 205: Based on the total admittance matrix, obtain the equivalent input impedance by matrix inversion, and evaluate the interaction stability between the new energy station and the grid-side impedance by the equivalent input impedance;

[0109] Step 206: Obtain the measured impedance, and dynamically adjust the model parameters of each stand-alone admittance matrix according to the measured impedance and the equivalent input impedance, while combining the adaptive learning rate matrix.

[0110] Reference Figure 3 , shows a structural block diagram of an impedance network modeling device for a new energy station provided by an embodiment of the present invention, wherein the new energy station includes multiple units; the device may specifically include:

[0111] The data acquisition unit 301 is used to acquire the station line data of the new energy station and the filter parameter data of each unit;

[0112] A dynamic link correction unit 302 is configured to construct a single-machine admittance matrix for each of the units based on the filtering parameter data and combined with dynamic link correction;

[0113] An admittance matrix construction unit 303 is configured to construct a coupling admittance matrix according to the station line data, and to construct a total admittance matrix according to the coupling admittance matrix and each of the stand-alone admittance matrices;

[0114] The equivalent input impedance solving unit 304 is used to obtain the equivalent input impedance by matrix inversion based on the total admittance matrix, and the equivalent input impedance is used to evaluate the interactive stability between the new energy station and the grid-side impedance.

[0115] In an optional embodiment, the filtering parameter data includes filter inductance, filter inductance equivalent resistance and grid equivalent impedance; the dynamic link correction unit 302 includes:

[0116] a phase-locked loop correction unit, configured to perform phase-locked loop correction based on the filter inductance and the grid equivalent impedance in combination with a frequency domain transfer function to obtain a dq-axis cross-coupling impedance and a qd-axis cross-coupling impedance of the converter in the unit;

[0117] a current control loop correction unit, configured to perform current control loop correction based on the filter inductor and the filter inductor equivalent resistance in combination with a current loop closed-loop transfer function, to obtain a d-axis self-impedance and a q-axis self-impedance of the converter in the unit;

[0118] a frequency domain impedance matrix construction unit, configured to construct a frequency domain impedance matrix of the unit according to the dq-axis cross-coupling impedance, the qd-axis cross-coupling impedance, the d-axis self-impedance, and the q-axis self-impedance;

[0119] The single-machine admittance matrix solving unit is used to obtain the single-machine admittance matrix of the unit by matrix inversion based on the frequency domain impedance matrix.

[0120] In an optional embodiment, the phase-locked loop correction process is performed by the following formula:

[0121]

[0122]

[0123] in, represents the dq axis cross-coupling impedance; represents the qd-axis cross-coupling impedance; represents the angular frequency; represents an imaginary unit; Indicates the DC bus voltage; represents the frequency domain transfer function; Represents the filter inductor; 、 PI controller parameters representing the frequency domain transfer function; Represents the equivalent impedance of the power grid.

[0124] In an optional embodiment, the current control loop correction process is performed by the following formula:

[0125]

[0126]

[0127]

[0128] in, represents the d-axis self-impedance; represents the q-axis self-impedance; represents the angular frequency; represents an imaginary unit; Represents the filter inductor; Indicates the equivalent resistance of the filter inductor; represents the current loop closed-loop transfer function; =1.5 Indicates the total delay time; represents the switching cycle; represents the frequency domain response of the PI controller in the current loop closed-loop transfer function; 、 PI controller parameters representing the current loop closed-loop transfer function.

[0129] In an optional embodiment, the station line data includes a station electrical wiring diagram, position coordinates of each unit, and line parameters of each line; the admittance matrix construction unit 303 includes:

[0130] A station topology connection matrix construction unit, wherein the user constructs a station topology connection matrix representing the line connection relationship in the new energy station according to the station electrical wiring diagram and the position coordinates of each unit;

[0131] A line impedance construction unit, configured to construct a line impedance for each of the lines according to the line parameters;

[0132] A coupling admittance matrix construction unit is used to construct a coupling admittance matrix characterizing the coupling of multiple machines in the new energy station according to the station topology connection matrix and each of the line impedances.

[0133] In an optional embodiment, the admittance matrix construction unit 303 includes:

[0134] A block diagonal matrix construction unit, configured to construct a block diagonal matrix based on each of the single-machine admittance matrices;

[0135] A total admittance matrix construction unit is used to construct a total admittance matrix by summing the block diagonal matrix and the coupling admittance matrix.

[0136] In an optional embodiment, the device further includes:

[0137] A model parameter dynamic parameter adjustment unit, configured to obtain a measured impedance, and dynamically adjust the model parameters of each of the stand-alone admittance matrices according to the measured impedance and the equivalent input impedance, while combining an adaptive learning rate matrix;

[0138] The model parameters include filter inductance, filter inductance equivalent resistance, PI controller parameters of frequency domain transfer function, and PI controller parameters of current loop closed-loop transfer function.

[0139] As for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the aforementioned method embodiment.

[0140] An embodiment of the present invention further provides an electronic device, the device including a processor and a memory:

[0141] The memory is used to store program codes and transmit the program codes to the processor;

[0142] The processor is used to execute the new energy station impedance network modeling method of any embodiment of the present invention according to the instructions in the program code.

[0143] An embodiment of the present invention further provides a computer-readable storage medium, which is used to store program code, and the program code is used to execute the new energy station impedance network modeling method of any embodiment of the present invention.

[0144] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0145] In the several embodiments provided by the present invention, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the units is merely a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interface, device or unit, which can be electrical, mechanical or other forms.

[0146] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of these units may be selected to achieve the purpose of this embodiment according to actual needs.

[0147] In addition, the functional units in the various embodiments of the present invention may be integrated into a single processing unit, each unit may exist physically separately, or two or more units may be integrated into a single unit. The aforementioned integrated units may be implemented in the form of hardware or software functional units.

[0148] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the portion that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes various media that can store program code, such as a USB flash drive, a mobile hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0149] As described above, the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that the technical solutions described in the above embodiments can still be modified, or some of the technical features thereof can be replaced by equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A new energy station impedance network modeling method, characterized in that: The new energy station includes multiple units; the method includes: Obtaining station line data of the new energy station and filter parameter data of each unit; For each of the units, a single-unit admittance matrix is ​​constructed based on the filtering parameter data and combined with dynamic link correction; Constructing a coupling admittance matrix according to the station line data, and constructing a total admittance matrix according to the coupling admittance matrix and each of the single-machine admittance matrices; Based on the total admittance matrix, an equivalent input impedance is obtained by matrix inversion, and the equivalent input impedance is used to evaluate the interactive stability between the new energy station and the grid-side impedance.

2. The impedance network modeling method for a new energy station according to claim 1 is characterized in that: The filtering parameter data includes filtering inductance, filtering inductance equivalent resistance and grid equivalent impedance; the single-machine admittance matrix is ​​constructed based on the filtering parameter data and combined with dynamic link correction, including: Based on the filter inductance and the grid equivalent impedance, a phase-locked loop correction is performed in combination with a frequency domain transfer function to obtain the dq-axis cross-coupling impedance and the qd-axis cross-coupling impedance of the converter in the unit; Based on the filter inductor and the filter inductor equivalent resistance, the current control loop is corrected in combination with the current loop closed-loop transfer function to obtain the d-axis self-impedance and q-axis self-impedance of the converter in the unit; constructing a frequency domain impedance matrix of the unit according to the dq-axis cross-coupling impedance, the qd-axis cross-coupling impedance, the d-axis self-impedance, and the q-axis self-impedance; Based on the frequency domain impedance matrix, the single-machine admittance matrix of the unit is obtained by matrix inversion.

3. The impedance network modeling method for a new energy station according to claim 2 is characterized in that: The phase-locked loop correction process is performed as follows: in, represents the dq axis cross-coupling impedance; represents the qd-axis cross-coupling impedance; represents the angular frequency; represents an imaginary unit; Indicates the DC bus voltage; represents the frequency domain transfer function; Represents the filter inductor; 、 PI controller parameters representing the frequency domain transfer function; Represents the equivalent impedance of the power grid.

4. The impedance network modeling method for a new energy station according to claim 2 is characterized in that: The current control loop correction process is performed by the following formula: in, represents the d-axis self-impedance; represents the q-axis self-impedance; represents the angular frequency; represents an imaginary unit; Represents the filter inductor; Indicates the equivalent resistance of the filter inductor; represents the current loop closed-loop transfer function; =1.5 Indicates the total delay time; represents the switching cycle; represents the frequency domain response of the PI controller in the current loop closed-loop transfer function; 、 PI controller parameters representing the current loop closed-loop transfer function.

5. The new energy station impedance network modeling method according to claim 2 is characterized in that: The station line data includes the station electrical wiring diagram, the location coordinates of each unit, and the line parameters of each line; The constructing of a coupling admittance matrix according to the station line data includes: Constructing a station topology connection matrix representing the line connection relationship in the new energy station based on the station electrical wiring diagram and the position coordinates of each unit; For each of the lines, constructing a line impedance according to the line parameters; A coupling admittance matrix characterizing multi-machine coupling in the new energy station is constructed based on the station topology connection matrix and each of the line impedances.

6. The impedance network modeling method for a new energy station according to claim 2 is characterized in that: The constructing of a total admittance matrix according to the coupled admittance matrix and each of the stand-alone admittance matrices includes: Based on each of the single-machine admittance matrices, a block diagonal matrix is ​​constructed; A total admittance matrix is ​​constructed by summing the block diagonal matrix and the coupled admittance matrix.

7. The impedance network modeling method for a new energy station according to any one of claims 2 to 6, characterized in that: Also includes: Obtaining a measured impedance, and dynamically adjusting model parameters of each of the stand-alone admittance matrices according to the measured impedance and the equivalent input impedance and in combination with an adaptive learning rate matrix; The model parameters include filter inductance, filter inductance equivalent resistance, PI controller parameters of frequency domain transfer function, and PI controller parameters of current loop closed-loop transfer function.

8. A new energy station impedance network modeling device, characterized in that: The new energy station includes multiple units; the device includes: A data acquisition unit, configured to acquire station line data of the new energy station and filter parameter data of each unit; A dynamic link correction unit is used to construct a single-machine admittance matrix for each of the units based on the filtering parameter data and combined with dynamic link correction; An admittance matrix construction unit, configured to construct a coupling admittance matrix according to the station line data, and to construct a total admittance matrix according to the coupling admittance matrix and each of the stand-alone admittance matrices; An equivalent input impedance solving unit is used to obtain an equivalent input impedance by matrix inversion based on the total admittance matrix, and the equivalent input impedance is used to evaluate the interactive stability between the new energy station and the grid-side impedance.

9. An electronic device, characterized in that: The device includes a processor and a memory: The memory is used to store program code and transmit the program code to the processor; The processor is used to execute the new energy station impedance network modeling method according to any one of claims 1-7 according to the instructions in the program code.

10. A computer-readable storage medium, characterized in that The computer-readable storage medium is used to store program code, and the program code is used to execute the new energy station impedance network modeling method according to any one of claims 1 to 7.

Images