Method, device, electronic device and medium for estimating parameters of synchronous condenser

Abstract

The present application discloses a method, device, electronic device and medium for estimating parameters of a synchronous condenser, which is applied to the field of synchronous condenser operation. The method provided by the present application includes: obtaining the operating data of the synchronous condenser before and after the load shedding respectively; preprocessing the operating data to obtain the incremental per-unit value corresponding to the operating data; determining the corresponding first-order integral value and second-order integral value according to the incremental per-unit value; determining the preliminary parameter value of the synchronous condenser according to the first-order integral value, the second-order integral value and the incremental per-unit value; and correcting the preliminary parameter value using the auxiliary variable method to obtain the target parameter value of the synchronous condenser. It can be seen that the step of determining the second-order integral value of the present application is applicable to the scenario where the excitation voltage changes during the load shedding period. At the same time, this method does not require a long measurement time, does not require providing the initial value of the parameter, does not require iteration, takes into account the change in the excitation voltage, has high parameter estimation accuracy, and is highly versatile and practical.

Description

Technical Field

[0001] The present application relates to the field of synchronous condenser operation, and in particular to a method, device, electronic equipment and medium for estimating synchronous condenser parameters. Background Art

[0002] Reasonable and reliable synchronous condenser parameters are essential for assessing their transient inoperative capacity, analyzing power system stability, and optimizing power system operation and control measures. Similar to the methods used to determine synchronous generator parameters, load-dump tests are often used to determine the direct-axis parameters of synchronous condensers. Currently, academic and engineering communities around the world have proposed various methods for estimating synchronous generator and / or synchronous condenser parameters based on load-dump tests. These methods can be broadly categorized into three types. One type is the traditional graphical method. This type of method requires the use of data when the generator and (or) phase regulator reaches steady state after load shedding, and the required measurement time is long. In addition, this method is not applicable to generators and (or) phase regulators whose excitation voltage cannot remain constant during load shedding. The second type is the numerical optimization method. This type of method treats the parameter estimation problem as a nonlinear optimization problem. When solving, the initial values ​​of the parameters need to be given, and then the estimation results are obtained through iterative calculation. When the initial values ​​of the parameters are unreasonable, it may lead to large estimation errors. The third type of method is the linear estimation method. This type of method uses the quadrature-axis voltage, excitation current, and direct-axis current to simultaneously identify two single-input and single-output transfer functions to obtain the direct-axis parameters of the generator and (or) phase regulator. However, this method is only applicable to scenarios where the excitation voltage remains constant during load shedding, and cannot handle scenarios where the excitation voltage changes during load shedding.

[0003] In view of the above technologies, seeking a method, device, electronic device and medium for estimating synchronous condenser parameters is an urgent problem to be solved by those skilled in the art. Summary of the Invention

[0004] The present application aims to provide a method, apparatus, electronic device, and medium for estimating synchronous condenser parameters. This approach addresses the problem in existing techniques where current synchronous condenser parameter estimation methods are only applicable to scenarios where the excitation voltage remains constant during load shedding, and are unable to address scenarios where the excitation voltage varies during load shedding.

[0005] To address the above technical issues, this application provides a method for estimating synchronous condenser parameters. This method is applicable to scenarios where the synchronous condenser's direct-axis equivalent circuit model includes one excitation winding and one damping winding. This model is widely accepted and used in power system stability analysis and simulation modeling. The method proposed in this application includes:

[0006] Acquiring operating data of the synchronous condenser before and after the load shedding respectively; wherein the operating data includes the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser before the load shedding, and the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser after the load shedding;

[0007] Preprocess the operating data to obtain the incremental per unit value corresponding to the operating data;

[0008] Determine the corresponding first-order integral value and second-order integral value according to the incremental per-unit value;

[0009] Determine the preliminary parameter values ​​of the synchronous condenser according to the first-order integral value, the second-order integral value and the incremental per unit value;

[0010] The auxiliary variable method is used to modify the preliminary parameter values ​​to obtain the target parameter values ​​of the synchronous condenser.

[0011] Preferably, the operating data is preprocessed to obtain incremental per-unit values ​​corresponding to the operating data;

[0012] Get the parameter reference values ​​corresponding to the synchronous condenser;

[0013] The corresponding incremental per unit value is determined based on the operating data before the load dump test, the operating data after the load dump test and the parameter reference value.

[0014] Preferably, obtaining the parameter reference value corresponding to the synchronous condenser includes:

[0015] According to the rated apparent power and rated voltage of the synchronous condenser, the armature current reference value, the excitation current reference value and the terminal voltage reference value are determined;

[0016] If the magnetic field voltage required to generate the rated stator terminal voltage of the synchronous condenser on the air gap line is known and reliable, the magnetic field voltage is used as the reference value of the excitation voltage;

[0017] If the magnetic field voltage required to generate the rated stator terminal voltage of the synchronous condenser on the air gap line is unknown or unreliable, the rated excitation voltage of the synchronous condenser is used as the excitation voltage reference value.

[0018] Preferably, when the magnetic field voltage required to generate the rated stator terminal voltage of the synchronous condenser on the air gap line is the excitation voltage reference value, determining the preliminary parameter value of the synchronous condenser according to the first-order integral value, the second-order integral value, and the incremental per-unit value includes:

[0019] Determine a first matrix according to the second-order integral value of the excitation voltage increment per unit value, the second-order integral value of the machine terminal voltage increment per unit value, and the second-order integral value of the excitation current increment per unit value;

[0020] Determine a second matrix according to the per-unit value of the machine-end voltage increment, the first-order integral value of the per-unit value of the machine-end voltage increment, the per-unit value of the excitation current increment, the first-order integral value of the per-unit value of the excitation current increment, the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the armature current increment, the second-order integral value of the per-unit value of the armature current increment, and the first-order integral value of the per-unit value of the excitation voltage increment;

[0021] A third matrix is ​​determined according to the first matrix and the second matrix, wherein the third matrix represents preliminary parameter values ​​of the synchronous condenser.

[0022] Preferably, when the rated excitation voltage is used as the excitation voltage reference value, determining the preliminary parameter value of the synchronous condenser according to the first-order integral value, the second-order integral value and the incremental per-unit value includes:

[0023] Determine a fourth matrix according to the second-order integral value of the per-unit value of the terminal voltage increment and the second-order integral value of the per-unit value of the excitation current increment;

[0024] Determine a fifth matrix according to the per-unit value of the machine-end voltage increment, the first-order integral value of the per-unit value of the machine-end voltage increment, the per-unit value of the excitation current increment, the first-order integral value of the per-unit value of the excitation current increment, the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the armature current increment, the second-order integral value of the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the excitation voltage increment, and the second-order integral value of the per-unit value of the excitation voltage increment;

[0025] A sixth matrix is ​​determined according to the fourth matrix and the fifth matrix, wherein the sixth matrix represents preliminary parameter values ​​of the synchronous condenser.

[0026] Preferably, the preliminary parameter values ​​are corrected using an auxiliary variable method to obtain target parameter values ​​of the synchronous condenser, including:

[0027] A first numerical simulation method is used to obtain a first time domain response of a per-unit value of a terminal voltage increment and a second time domain response of a per-unit value of an excitation current increment;

[0028] Obtain a first-order integral value of the first time domain response and a first-order integral value of the second time domain response;

[0029] determining a first auxiliary variable matrix according to the first time domain response, the first-order integral value of the first time domain response, the second time domain response, the first-order integral value of the second time domain response, the armature current increment per unit value, the first-order integral value of the armature current increment per unit value, the second-order integral value of the armature current increment per unit value, and the first-order integral value of the excitation voltage increment per unit value;

[0030] A first target matrix is ​​determined according to the first matrix, the first auxiliary variable matrix and the second matrix, where the first target matrix represents target parameter values ​​of the synchronous condenser.

[0031] Preferably, the preliminary parameter values ​​are corrected using an auxiliary variable method to obtain target parameter values ​​of the synchronous condenser, including:

[0032] A second numerical simulation method is used to obtain a third time domain response of the per-unit value of the terminal voltage increment and a fourth time domain response of the per-unit value of the excitation current increment;

[0033] determining a second auxiliary variable matrix according to the third time domain response, the first-order integral value of the third time domain response, the fourth time domain response, the first-order integral value of the fourth time domain response, the armature current increment per unit value, the first-order integral value of the armature current increment per unit value, the second-order integral value of the armature current increment per unit value, and the first-order integral value of the excitation voltage increment per unit value;

[0034] A second target matrix is ​​determined according to the fourth matrix, the second auxiliary variable matrix and the fifth matrix, where the second target matrix represents target parameter values ​​of the synchronous condenser.

[0035] To solve the above technical problems, the present application further provides a device for estimating parameters of a synchronous condenser, comprising:

[0036] an acquisition module, configured to acquire operating data of the synchronous condenser before and after the load dump test; wherein the operating data includes the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser before the load dump, and the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser after the load dump; and

[0037] A processing module, used for pre-processing the operating data to obtain the incremental per unit value corresponding to the operating data;

[0038] A first determining module is used to determine a corresponding first-order integral value and a second-order integral value according to the incremental per-unit value;

[0039] A second determining module is used to determine preliminary parameter values ​​of the synchronous condenser according to the first-order integral value, the second-order integral value and the incremental per-unit value;

[0040] The correction module is used to correct the preliminary parameter values ​​using the auxiliary variable method to obtain the target parameter values ​​of the synchronous condenser.

[0041] To solve the above technical problems, the present application further provides an electronic device, comprising a memory for storing a computer program;

[0042] A processor is configured to implement the steps of the method for estimating synchronous condenser parameters when executing a computer program.

[0043] In order to solve the above technical problems, the present application also provides an electronic device and a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the method for estimating synchronous phase condenser parameters as described above are implemented.

[0044] The present application provides a method for estimating parameters of a synchronous condenser, comprising: obtaining operating data of the synchronous condenser before and after load shedding, respectively; wherein the operating data include the armature current, the machine terminal voltage, the excitation current and the excitation voltage of the synchronous condenser before load shedding, and the armature current, the machine terminal voltage, the excitation current and the excitation voltage of the synchronous condenser after load shedding; preprocessing the operating data to obtain an incremental per-unit value corresponding to the operating data; determining a corresponding first-order integral value and a second-order integral value according to the incremental per-unit value; determining a preliminary parameter value of the synchronous condenser according to the first-order integral value, the second-order integral value and the incremental per-unit value; and correcting the preliminary parameter value using an auxiliary variable method to obtain a target parameter value of the synchronous condenser. It can be seen that the method provided in the present application is based on the load dump test, obtains the operating data before and after the load dump test, and processes the operating data to obtain the corresponding first-order integral values ​​and second-order integral values. On the basis of obtaining the preliminary parameter values ​​of the synchronous phase regulator, the preliminary parameter values ​​are corrected according to the auxiliary variable method to provide more prepared target parameter values. In the present application, the step of processing the operating data to obtain the corresponding second-order integral values ​​is applicable to the scenario where the excitation voltage changes during the load dump period. The present application expands the scope of use of the method. At the same time, the method does not require a long measurement time, does not require providing initial parameter values, does not require iteration, takes into account the excitation voltage changes, has high parameter estimation accuracy, and is highly universal and practical. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In order to more clearly illustrate the embodiments of the present application, the following is a brief introduction to the drawings required for use in the embodiments. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0046] Figure 1 A flowchart of a method for estimating synchronous condenser parameters provided in an embodiment of the present application;

[0047] Figure 2 A schematic diagram of a simulation model provided in an embodiment of the present application;

[0048] Figure 3A schematic diagram of the true values ​​and estimated values ​​of parameters provided in the embodiments of the present application;

[0049] Figure 4 A module diagram of an apparatus for estimating synchronous condenser parameters provided by another embodiment of the present application;

[0050] Figure 5 A structural diagram of an electronic device provided in another embodiment of the present application. DETAILED DESCRIPTION

[0051] The following will be combined with the accompanying drawings in the embodiments of this application to clearly and completely describe the technical solutions in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0052] The core of this application is to provide a method, device, electronic equipment and medium for estimating synchronous condenser parameters.

[0053] In order to enable those skilled in the art to better understand the present application, the present application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0054] Reasonable and reliable synchronous condenser parameters are essential for assessing their transient inoperative capacity, analyzing power system stability, and optimizing power system operation and control measures. Similar to the methods used to determine synchronous generator parameters, load-dump tests are often used to determine the direct-axis parameters of synchronous condensers. Currently, academic and engineering communities around the world have proposed various methods for estimating synchronous generator and / or synchronous condenser parameters based on load-dump tests. These methods can be broadly categorized into three types. One type is the traditional graphical method. This type of method requires the use of data when the generator and (or) phase regulator reaches steady state after load shedding, and the required measurement time is long. In addition, this method is not applicable to generators and (or) phase regulators whose excitation voltage cannot remain constant during load shedding. The second type is the numerical optimization method. This type of method treats the parameter estimation problem as a nonlinear optimization problem. When solving, the initial values ​​of the parameters need to be given, and then the estimation results are obtained through iterative calculation. When the initial values ​​of the parameters are unreasonable, it may lead to large estimation errors. The third type of method is the linear estimation method. This type of method uses the quadrature-axis voltage, excitation current, and direct-axis current to simultaneously identify two single-input and single-output transfer functions to obtain the direct-axis parameters of the generator and (or) phase regulator. However, this method is only applicable to scenarios where the excitation voltage remains constant during load shedding, and cannot handle scenarios where the excitation voltage changes during load shedding.

[0055] In order to solve the above technical problems, the present application provides a method for estimating synchronous condenser parameters, such as Figure 1As shown, the following steps are included:

[0056] S10: Acquire the operating data of the synchronous condenser before and after the load shedding respectively.

[0057] In a specific embodiment, the load dump test is to suddenly disconnect the load while the synchronous condenser is connected to the load. In this application, the operating data of the synchronous condenser before and after the load dump is obtained. Generally speaking, the operating data includes four types, namely armature current, machine terminal voltage, excitation current and excitation voltage. Among them, the armature current is: the current flowing through the armature winding of the synchronous condenser; the machine terminal voltage is: the output voltage of the synchronous condenser machine terminal, which can be a phase voltage or a line voltage; the excitation current is: the output current of the synchronous condenser exciter; the excitation voltage is: the output voltage of the synchronous condenser exciter.

[0058] Since the present application obtains the operating data before and after the load shedding, the operating data in the embodiment of the present application specifically include: the armature current before the load shedding, the machine terminal voltage before the load shedding, the excitation current before the load shedding and the excitation voltage before the load shedding, the armature current after the load shedding, the machine terminal voltage after the load shedding, the excitation current after the load shedding and the excitation voltage after the load shedding.

[0059] In a specific embodiment, the armature current, machine terminal voltage, excitation current and excitation voltage can be obtained by using corresponding sensors, etc., but this application is not limited to this. At the same time, this application does not limit the frequency of acquisition, which can be set according to user needs.

[0060] S11: Preprocess the operating data to obtain the incremental per-unit value corresponding to the operating data.

[0061] In a specific embodiment, the armature current before load shedding, the machine terminal voltage before load shedding, the excitation current before load shedding and the excitation voltage before load shedding, and the armature current, the machine terminal voltage after load shedding, the excitation current after load shedding and the excitation voltage after load shedding are preprocessed to obtain corresponding incremental per-unit values, that is, the incremental per-unit value of the armature current; the incremental per-unit value of the machine terminal voltage, the incremental per-unit value of the excitation current and the incremental per-unit value of the excitation voltage.

[0062] Wherein, as a preferred embodiment, the incremental per-unit value is calculated as follows: the ratio of the difference between the operating data after the load dump and the operating data before the load dump to the corresponding reference value is the incremental per-unit value. Taking the incremental per-unit value of the armature current as an example: the ratio of the difference between the armature current after the load dump and the armature current before the load dump to the armature current reference value is the incremental per-unit value of the armature current. Wherein, it should be noted that the reference value is determined by the technical parameters of the synchronous phase regulator. It should also be noted that the example of the present application is only one possible way of implementation, but is not limited to only this implementation method, and can be set according to the needs of the user.

[0063] S12: Determine the corresponding first-order integral value and second-order integral value according to the incremental per-unit value.

[0064] S13: Determine preliminary parameter values ​​of the synchronous condenser according to the first-order integral value, the second-order integral value, and the incremental per-unit value.

[0065] In a specific embodiment, the first-order integral value and the second-order integral value corresponding to the incremental per-unit value determined in the above steps are obtained according to the calculation formulas for the first-order integral and the second-order integral. Specific values ​​include: the first-order integral value of the incremental per-unit value of the armature current, the second-order integral value of the incremental per-unit value of the armature current; the first-order integral value of the incremental per-unit value of the generator-end voltage, the second-order integral value of the incremental per-unit value of the generator-end voltage; the first-order integral value of the incremental per-unit value of the excitation current, the second-order integral value of the incremental per-unit value of the excitation current; the first-order integral value of the incremental per-unit value of the excitation voltage, the second-order integral value of the incremental per-unit value of the excitation voltage.

[0066] Since the operating data and the parameters of the synchronous condenser have a certain corresponding relationship, the preliminary parameter values ​​of the synchronous condenser are determined according to the first-order integral value, the second-order integral value and the incremental per-unit value.

[0067] The specific correspondence between the operating data and the parameters of the synchronous condenser is not limited in this application and can be set according to the needs of the user.

[0068] S14: Using the auxiliary variable method to correct the preliminary parameter value to obtain the target parameter value of the synchronous condenser.

[0069] In a specific embodiment, the accuracy of the preliminary parameter value obtained in step S13 is low, so the auxiliary variable method is used to correct the preliminary parameter value to obtain the target parameter value of the synchronous condenser, that is, to obtain a more accurate target parameter value.

[0070] Among them, the auxiliary variable method is defined as: if a certain variable is highly correlated with the random explanatory variable in the model, but is not correlated with the random error term, then a consistent estimator can be obtained by using this variable and the corresponding regression coefficient in the model. This variable is called an auxiliary variable, and this estimation method is called the auxiliary variable method.

[0071] The present application provides a method for estimating parameters of a synchronous condenser, comprising: obtaining operating data of the synchronous condenser before and after load shedding, respectively; wherein the operating data include the armature current, the machine terminal voltage, the excitation current and the excitation voltage of the synchronous condenser before load shedding, and the armature current, the machine terminal voltage, the excitation current and the excitation voltage of the synchronous condenser after load shedding; preprocessing the operating data to obtain an incremental per-unit value corresponding to the operating data; determining a corresponding first-order integral value and a second-order integral value according to the incremental per-unit value; determining a preliminary parameter value of the synchronous condenser according to the first-order integral value, the second-order integral value and the incremental per-unit value; and correcting the preliminary parameter value using an auxiliary variable method to obtain a target parameter value of the synchronous condenser. It can be seen that the method provided in the present application is based on the load dump test, obtains the operating data before and after the load dump test, and processes the operating data to obtain the corresponding first-order integral values ​​and second-order integral values. On the basis of obtaining the preliminary parameter values ​​of the synchronous phase regulator, the preliminary parameter values ​​are corrected according to the auxiliary variable method to provide more prepared target parameter values. In the present application, the step of processing the operating data to obtain the corresponding second-order integral values ​​is applicable to the scenario where the excitation voltage changes during the load dump period. The present application expands the scope of use of the method. At the same time, the method does not require a long measurement time, does not require providing initial parameter values, does not require iteration, takes into account the excitation voltage changes, has high parameter estimation accuracy, and is highly universal and practical.

[0072] Based on the above embodiment, as a preferred embodiment, the operation data is pre-processed to obtain the incremental per unit value corresponding to the operation data;

[0073] Get the parameter reference values ​​corresponding to the synchronous condenser;

[0074] The corresponding incremental per unit value is determined based on the operating data before the load dump test, the operating data after the load dump test and the parameter reference value.

[0075] In a specific embodiment, in order to obtain the incremental per-unit value corresponding to the operating data, the parameter reference values ​​corresponding to the synchronous phase regulator are first obtained, that is, the armature current reference value, the machine terminal voltage reference value, the excitation current reference value and the excitation voltage reference value. Then, the corresponding incremental per-unit value is determined based on the operating data before the load dump test, the operating data after the load dump test and the parameter reference value. The specific method is: the ratio of the difference between the operating data after the load dump and the operating data before the load dump to the corresponding reference value is the incremental per-unit value. Taking the incremental per-unit value of the armature current as an example: the difference between the armature current after the load dump and the stable armature current before the load dump and the armature current reference value is the incremental per-unit value of the armature current. It should be noted that the reference value is determined by the technical parameters of the synchronous phase regulator. It should also be noted that the example of this application is only one way that can be implemented, but is not limited to only this implementation method, and can be set according to the needs of the user.

[0076] The armature current reference value, excitation current reference value, and generator terminal voltage reference value are determined based on the synchronous condenser's rated apparent power and rated voltage. However, there are two methods for determining the excitation voltage reference value. The first method is to use the field voltage as the excitation voltage reference value when the field voltage required to generate the synchronous condenser's rated stator terminal voltage on the air gap line is known and reliable. In other words, when the field voltage is reliable, the field voltage is used as the excitation voltage reference value (this reference value is also the reference value for the non-reversible per-unit system commonly used in excitation systems). The second method is to use the synchronous condenser's rated excitation voltage as the excitation voltage reference value when the field voltage required to generate the synchronous condenser's rated stator terminal voltage on the air gap line is unknown or unreliable.

[0077] According to the above content, there are two ways to select the excitation voltage reference value. Therefore, the subsequent steps of obtaining the preliminary parameter values ​​of the synchronous condenser and correcting the preliminary parameter values ​​will differ due to the different selection of the excitation voltage reference value.

[0078] When the magnetic field voltage is the excitation voltage reference value, the steps of obtaining the preliminary parameter value of the synchronous condenser include:

[0079] Determine a first matrix according to the second-order integral value of the excitation voltage increment per unit value, the second-order integral value of the machine terminal voltage increment per unit value, and the second-order integral value of the excitation current increment per unit value;

[0080] Determine a second matrix according to the per-unit value of the machine-end voltage increment, the first-order integral value of the per-unit value of the machine-end voltage increment, the per-unit value of the excitation current increment, the first-order integral value of the per-unit value of the excitation current increment, the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the armature current increment, the second-order integral value of the per-unit value of the armature current increment, and the first-order integral value of the per-unit value of the excitation voltage increment;

[0081] A third matrix is ​​determined according to the first matrix and the second matrix, wherein the third matrix represents preliminary parameter values ​​of the synchronous condenser.

[0082] In a specific embodiment, the armature current after load shedding, the terminal voltage after load shedding, the excitation voltage after load shedding, and the excitation current after load shedding are all obtained by taking N sampling points starting from the first sampling point after load shedding. This means that the data after load shedding are N groups, each group corresponding to one armature current after load shedding, terminal voltage after load shedding, excitation voltage after load shedding, and excitation current after load shedding. Wherein each sampling time is t i Indicates (i=1, 2, ..., N).

[0083] Among them, the expression of the first matrix is:

[0084] Γ=[y1(t1) … y1(t N ) y2(t1) … y2(t N )] T ;

[0085]

[0086]

[0087] in, is the second-order integral value of the per-unit value of the excitation voltage increment; It is the second-order integral value of the per-unit value of the terminal voltage increment; It is the second-order integral value of the per-unit value of the excitation current increment.

[0088] The expression of the second matrix is:

[0089] Φ=[AB];

[0090]

[0091]

[0092] Where Δv t (t i ) is the per-unit value of the terminal voltage increment; ΔI is the first-order integral value of the per-unit value of the terminal voltage increment; fd (t i ) is the per-unit value of the excitation current increment; is the first-order integral value of the excitation current increment per unit value; Δi d (t i ) is the per-unit value of the armature current increment; is the first-order integral value of the per-unit value of the armature current increment; is the first-order integral value of the per-unit value of the armature current increment; It is the first-order integral value of the per-unit value of the excitation voltage increment.

[0093] The third matrix is ​​determined based on the first matrix and the second matrix. The relationship between the first matrix, the second matrix and the third matrix is:

[0094] Γ = Φθ;

[0095] θ=[θ1 θ2 θ3 θ4 θ5 θ6 θ7 θ8 θ9];

[0096] Where θ is the third matrix. θ1=T′ d0 T″ d0 ,θ2=T′ d0 +T″ d0 ,θ3=L d T′ d T″ d ,θ4=L d (T′ d +T″ d ), θ5=L d , Among them, T′ d0 , T″ d0 , T′ d , T″ d They represent the direct-axis transient open-circuit time constant, the direct-axis subtransient open-circuit time constant, the direct-axis transient short-circuit time constant, and the direct-axis subtransient short-circuit time constant, respectively, all in seconds; L d , L md , L 1d , L f1d 、R 1d 、R fd They represent the direct-axis synchronous inductance, direct-axis synchronous mutual inductance, direct-axis damping winding leakage inductance, differential leakage inductance, direct-axis damping winding resistance, and excitation winding resistance, respectively, all of which are per-unit values; ω n Indicates the rated angular frequency in rad / s.

[0097] After the values ​​of the third matrix are determined by the first and second matrices, each element in the third matrix is ​​the preliminary parameter value obtained in this application.

[0098] Accordingly, the steps for correcting the preliminary parameter values ​​include:

[0099] A first numerical simulation method is used to obtain a first time domain response of a per-unit value of a terminal voltage increment and a second time domain response of a per-unit value of an excitation current increment;

[0100] Obtain a first-order integral value of the first time domain response and a first-order integral value of the second time domain response;

[0101] determining a first auxiliary variable matrix according to the first time domain response, the first-order integral value of the first time domain response, the second time domain response, the first-order integral value of the second time domain response, the armature current increment per unit value, the first-order integral value of the armature current increment per unit value, the second-order integral value of the armature current increment per unit value, and the first-order integral value of the excitation voltage increment per unit value;

[0102] A first target matrix is ​​determined according to the first matrix, the first auxiliary variable matrix and the second matrix, where the first target matrix represents target parameter values ​​of the synchronous condenser.

[0103] In a specific embodiment, the calculation formula of the first time domain response of the per-unit value of the terminal voltage increment is:

[0104]

[0105] The calculation formula for the second time domain response of the per-unit value of the excitation current is:

[0106]

[0107] in, is the first time domain response of the per-unit value of the terminal voltage increment, is the second time domain response of the per-unit value of the excitation current, p represents the differential operator,

[0108] The dimension of the first auxiliary variable matrix Z is the same as that of Φ, and the last 7 columns are the same, and the expressions of the first two columns of the first auxiliary variable matrix are:

[0109]

[0110] in, is the first-order integral value of the first time domain response; is the first-order integral value of the second time domain response,

[0111] Among them, the expression of the corresponding relationship between the first matrix, the first auxiliary variable matrix, the second matrix and the first target matrix is:

[0112]

[0113] in, is the first target matrix. According to the first target matrix and the above formula for parameters, the target parameter values ​​can be obtained, where the types of target parameter values ​​include: direct axis transient open circuit time constant T′ d0 ; Direct axis subtransient open circuit time constant T″ d0 ; Direct axis transient short circuit time constant T′ d ; Direct axis subtransient short circuit time constant T″ d; Direct-axis synchronous inductor L d ; Direct-axis transient inductance L′ d ; Direct-axis subtransient inductance L″ d .

[0114] When the rated excitation voltage is used as the excitation voltage reference value, the steps for obtaining the preliminary parameter values ​​include:

[0115] Determine a fourth matrix according to the second-order integral value of the per-unit value of the terminal voltage increment and the second-order integral value of the per-unit value of the excitation current increment;

[0116] Determine a fifth matrix according to the per-unit value of the machine-end voltage increment, the first-order integral value of the per-unit value of the machine-end voltage increment, the per-unit value of the excitation current increment, the first-order integral value of the per-unit value of the excitation current increment, the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the armature current increment, the second-order integral value of the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the excitation voltage increment, and the second-order integral value of the per-unit value of the excitation voltage increment;

[0117] A sixth matrix is ​​determined according to the fourth matrix and the fifth matrix, wherein the sixth matrix represents preliminary parameter values ​​of the synchronous condenser.

[0118] In a specific embodiment, the armature current after load shedding, the terminal voltage after load shedding, the excitation voltage after load shedding, and the excitation current after load shedding are obtained by taking N sampling points starting from the first sampling point after load shedding. In this case, the data after load shedding are N groups, each group corresponding to one armature current after load shedding, terminal voltage after load shedding, excitation voltage after load shedding, and excitation current after load shedding. The sampling time is t i Indicates (i=1, 2, ..., N).

[0119] Among them, the expression of the fourth matrix is:

[0120] Γ′=[y1′(t1) … y1′(t N ) y2′(t1) …y2′(t N )] T ;

[0121]

[0122]

[0123] in, It is the second-order integral value of the per-unit value of the terminal voltage increment; It is the second-order integral value of the per-unit value of the excitation current increment.

[0124] The expression of the fifth matrix is:

[0125] Φ′=[AB′];

[0126]

[0127]

[0128] Where Δv t (t i ) is the per-unit value of the terminal voltage increment; ΔI is the first-order integral value of the per-unit value of the terminal voltage increment; fd (t i ) is the per-unit value of the excitation current increment; is the first-order integral value of the excitation current increment per unit value; Δi d (t i ) is the per-unit value of the armature current increment; is the first-order integral value of the per-unit value of the armature current increment; is the first-order integral value of the per-unit value of the armature current increment; It is the first-order integral value of the per-unit value of the excitation voltage increment; It is the second-order integral value of the per-unit value of the excitation voltage increment.

[0129] The sixth matrix is ​​determined according to the fourth matrix and the fifth matrix. Then the relationship between the fourth matrix, the fifth matrix and the sixth matrix is:

[0130] Γ′=Φ′θ′;

[0131] θ′=[θ1 θ2 θ3 θ4 θ5 θ6′ θ7 θ8 θ9′ θ 10 ′ θ 11 ′]

[0132] Where θ′ is the sixth matrix. θ1=T′ d0 T″ d0 ,θ2=T′ d0 +T′ d0 ,θ3=L d T′ d T′ d ,θ4=L d (T′ d +T″ d ), θ5=L d , θ9′=K 2 , Among them, T′ d0 , T″ d0 , T′ d , T″ dThey represent the direct-axis transient open-circuit time constant, the direct-axis subtransient open-circuit time constant, the direct-axis transient short-circuit time constant, and the direct-axis subtransient short-circuit time constant, respectively, all in seconds; L d , L md , L 1d , L f1d 、R 1d 、R fd They represent the direct-axis synchronous inductance, direct-axis synchronous mutual inductance, direct-axis damping winding leakage inductance, differential leakage inductance, direct-axis damping winding resistance, and excitation winding resistance, respectively, all of which are per-unit values; ω n Indicates the rated angular frequency in rad / s.

[0133] After the values ​​of the sixth matrix are determined by the fourth and fifth matrices, the elements in the sixth matrix are the preliminary parameter values ​​obtained when the rated excitation voltage is used as the excitation voltage reference value in the present application.

[0134] Accordingly, the steps for correcting the preliminary parameter values ​​include:

[0135] A second numerical simulation method is used to obtain a third time domain response of the per-unit value of the terminal voltage increment and a fourth time domain response of the per-unit value of the excitation current increment;

[0136] determining a second auxiliary variable matrix according to the third time domain response, the first-order integral value of the third time domain response, the fourth time domain response, the first-order integral value of the fourth time domain response, the armature current increment per unit value, the first-order integral value of the armature current increment per unit value, the second-order integral value of the armature current increment per unit value, and the first-order integral value of the excitation voltage increment per unit value;

[0137] A second target matrix is ​​determined according to the fourth matrix, the second auxiliary variable matrix and the fifth matrix, where the second target matrix represents target parameter values ​​of the synchronous condenser.

[0138] In a specific embodiment, the expression of the third time domain response of the per-unit value of the terminal voltage increment is:

[0139]

[0140] The expression of the fourth time domain response of the per-unit value of the excitation current increment is:

[0141]

[0142] in, The third time domain response of the per-unit value of the terminal voltage increment; is the fourth time domain response of the per-unit value of the excitation current increment, p represents the differential operator,

[0143] The dimension of the second auxiliary variable matrix Z′ is the same as that of Φ′, and the last 9 columns are the same, and the expressions of the first two columns of the second auxiliary variable matrix are:

[0144]

[0145] in, is the first-order integral value of the third time domain response; is the first-order integral value corresponding to the fourth time domain,

[0146] The corresponding relationship between the fourth matrix, the second auxiliary variable matrix, the fifth matrix and the second target matrix is ​​expressed as follows:

[0147]

[0148] in, is the second target matrix. According to the second target matrix and the above formula for parameters, the target parameter values ​​can be obtained, where the types of target parameter values ​​include: direct axis transient open circuit time constant T′ d0 ; Direct axis subtransient open circuit time constant T″ d0 ; Direct axis transient short circuit time constant T′ d ; Direct axis subtransient short circuit time constant T″ d ; Direct-axis synchronous inductor L d ; Direct axis transient inductance L' d ; Direct-axis subtransient inductance L″ d .

[0149] The effectiveness of the method proposed in this invention is verified by simulation examples. The simulation model is as follows: Figure 2 In the simulation model, the rated capacity of the synchronous condenser is 300 Mvar. The steady-state operating conditions before load shedding are: the condenser absorbs 50 Mvar of reactive power. The per-unit values ​​of the excitation voltage and current in the non-reversible per-unit system are 0.7745 and 0.7745, respectively. In the reversible per-unit system, the per-unit values ​​of the direct-axis current and terminal voltage are -0.1168 and 0.9952, respectively. The sampling frequency is set to 2000 Hz.

[0150] First, assuming that the excitation voltage reference value under the irreversible per-unit system is known, the first target matrix is ​​obtained by solving the corresponding formula above: The value of L d , T′ d , T″ d , T′ d0 , T″ d0 , L′ d , L″ d The comparison between these parameter estimates and the true values ​​is shown in Table 1 below. It can be seen that the error between the parameter estimates and the true values ​​is small.d , T′ d , T″ d , T′ d0 , T″ d0 , L′ d , L″ d The relative error is less than 4%. Except for the relative error of θ7 which is slightly greater than 10%, the relative errors of the other parameters are all less than 7.5%.

[0151] Table 1

[0152]

[0153]

[0154] Secondly, assuming that the excitation voltage reference value of the irreversible per-unit system is unknown, θ 10 ′=2.0, solve according to the corresponding formula above, and get the first target matrix The value of L d , T′ d , T″ d , T′ d0 , T″ d0 , L′ d 、L d The comparison between these parameter estimates and the true values ​​is shown in Table 2 below. It can be seen that the error between the parameter estimates and the true values ​​is small, except for θ′ 10 Except for the relative error exceeding 10%, the relative errors of other parameters are all less than 3.5%. Based on the true value and estimated value of the parameters, the load dump test simulation is carried out, and the terminal voltage obtained is as follows Figure 3 As shown, it can be seen that the voltage response at the machine end tends to be consistent.

[0155] Table 2

[0156] parameter truth value Estimated value Relative error (%) <![CDATA[θ1]]> 0.1845 0.1882 2.01 <![CDATA[θ2]]> 5.7836 5.9118 2.22 <![CDATA[θ3]]> 0.0374 0.037 1.07 <![CDATA[θ4]]> 1.6793 1.7159 2.18 <![CDATA[θ5(L d )]]> 1.89 1.9302 2.13 <![CDATA[θ6′]]> 0.0062 0.006 3.23 <![CDATA[θ7′]]> 1.5 1.5326 2.17 <![CDATA[θ8′]]> 1.0944 1.113 1.70 <![CDATA[θ9′]]> 2.25 2.25 0.00 <![CDATA[θ 10 ′]]> 0.0522 0.0587 12.45 <![CDATA[θ 11 ′]]> 12.6478 12.9867 2.68 <![CDATA[T′ d ]]> 0.8656 0.8669 0.15 <![CDATA[T″ d ]]> 0.0229 0.0221 3.49 <![CDATA[T′ d0 ]]> 5.7515 5.8798 2.23 <![CDATA[T″ d0 ]]> 0.0321 0.032 0.31 <![CDATA[L′ d ]]> 0.2871 0.2874 0.10 <![CDATA[L″ d ]]> 0.2028 0.1967 3.01

[0157] The present application provides a method for estimating parameters of a synchronous condenser, comprising: obtaining operating data of the synchronous condenser before and after load shedding, respectively; wherein the operating data include the armature current, the machine terminal voltage, the excitation current and the excitation voltage of the synchronous condenser before load shedding, and the armature current, the machine terminal voltage, the excitation current and the excitation voltage of the synchronous condenser after load shedding; preprocessing the operating data to obtain an incremental per-unit value corresponding to the operating data; determining a corresponding first-order integral value and a second-order integral value according to the incremental per-unit value; determining a preliminary parameter value of the synchronous condenser according to the first-order integral value, the second-order integral value and the incremental per-unit value; and correcting the preliminary parameter value using an auxiliary variable method to obtain a target parameter value of the synchronous condenser. It can be seen that the method provided in the present application is based on the load dump test, obtains the operating data before and after the load dump test, and processes the operating data to obtain the corresponding first-order integral values ​​and second-order integral values. On the basis of obtaining the preliminary parameter values ​​of the synchronous phase regulator, the preliminary parameter values ​​are corrected according to the auxiliary variable method to provide more prepared target parameter values. In the present application, the step of processing the operating data to obtain the corresponding second-order integral values ​​is applicable to the scenario where the excitation voltage changes during the load dump period. The present application expands the scope of use of the method. At the same time, the method does not require a long measurement time, does not require providing initial parameter values, does not require iteration, takes into account the excitation voltage changes, has high parameter estimation accuracy, and is highly universal and practical.

[0158] In the above embodiments, a method for estimating synchronous condenser parameters is described in detail. This application also provides corresponding embodiments of an apparatus for estimating synchronous condenser parameters. It should be noted that this application describes embodiments of the apparatus from two perspectives: one based on functional modules and the other based on hardware.

[0159] Figure 4 A module diagram of an apparatus for estimating synchronous condenser parameters provided in another embodiment of the present application includes:

[0160] an acquisition module 11 for respectively acquiring operating data of the synchronous condenser before and after the load shedding; wherein the operating data includes the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser before the load shedding, and the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser after the load shedding; and

[0161] The processing module 12 is used to pre-process the operation data to obtain the incremental per unit value corresponding to the operation data;

[0162] A first determining module 13 is configured to determine a corresponding first-order integral value and a second-order integral value according to the incremental per-unit value;

[0163] A second determining module 14 is configured to determine preliminary parameter values ​​of the synchronous condenser according to the first-order integral value, the second-order integral value, and the incremental per-unit value;

[0164] The correction module 15 is used to correct the preliminary parameter values ​​using the auxiliary variable method to obtain the target parameter values ​​of the synchronous condenser.

[0165] Since the embodiments of the apparatus part correspond to the embodiments of the method part, please refer to the description of the embodiments of the method part for the embodiments of the apparatus part, and they will not be repeated here.

[0166] Figure 5 This is a structural diagram of an electronic device provided in another embodiment of the present application, such as Figure 5 As shown, the electronic device includes: a memory 20 for storing computer programs;

[0167] The processor 21 is configured to implement the steps of the method for estimating synchronous condenser parameters mentioned in the above embodiment when executing the computer program.

[0168] The electronic device provided in this embodiment may include but is not limited to a smart phone, a tablet computer, a laptop computer, or a desktop computer.

[0169] Among them, the processor 21 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 21 can be implemented in at least one hardware form of a digital signal processor (DSP), a field programmable gate array (FPGA), and a programmable logic array (PLA). The processor 21 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a central processing unit (CPU); the coprocessor is a low-power processor for processing data in the standby state. In some embodiments, the processor 21 may be integrated with a graphics processing unit (GPU), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 21 may also include an artificial intelligence (AI) processor, which is used to process computing operations related to machine learning.

[0170] The memory 20 may include one or more computer-readable storage media, which may be non-transitory. The memory 20 may also include high-speed random access memory and non-volatile memory, such as one or more disk storage devices or flash memory storage devices. In this embodiment, the memory 20 is used to store at least the following computer program 201. When loaded and executed by the processor 21, the computer program can implement the relevant steps of the method for estimating synchronous condenser parameters disclosed in any of the aforementioned embodiments. Furthermore, the resources stored in the memory 20 may include an operating system 202 and data 203, which may be stored in either a temporary or permanent manner. The operating system 202 may include Windows, Unix, Linux, etc.

[0171] In some embodiments, the electronic device may further include a display screen 22 , an input / output interface 23 , a communication interface 24 , a power supply 25 , and a communication bus 26 .

[0172] Those skilled in the art will understand that Figure 5 The structure shown in the figure does not constitute a limitation of the electronic device, and may include more or fewer components than shown in the figure.

[0173] The electronic device provided in an embodiment of the present application includes a memory and a processor. When the processor executes the program stored in the memory, it can implement the above-described method for estimating synchronous condenser parameters and has the same beneficial effects.

[0174] Finally, the present application also provides an embodiment corresponding to a computer-readable storage medium. The computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps described in the above method embodiment.

[0175] It is understandable that if the method in the above embodiment is implemented in the form of 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 application, or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, which is stored in a storage medium and executes all or part of the steps of the method described in each embodiment of the present application. The aforementioned storage medium includes various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0176] The above describes in detail the method, device, electronic device, and medium for estimating synchronous condenser parameters provided by the present application. The various embodiments in the specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. Reference can be made to the same or similar parts between the various embodiments. For the devices disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple, and for relevant parts, reference can be made to the method description. It should be noted that, for those skilled in the art, several improvements and modifications can be made to the present application without departing from the principles of the present application, and these improvements and modifications also fall within the scope of protection of the claims of the present application.

[0177] It should also be noted that, in this specification, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variants thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or apparatus. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article, or apparatus comprising the element.

Claims

1. A method for estimating parameters of a synchronous condenser, characterized in that: include: Acquiring operating data of the synchronous condenser before and after the load shedding respectively; wherein the operating data includes the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser before the load shedding, and the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser after the load shedding; Preprocessing the operating data to obtain an incremental per-unit value corresponding to the operating data; Determine the corresponding first-order integral value and second-order integral value according to the incremental per-unit value; Determining preliminary parameter values ​​of the synchronous condenser according to the first-order integral value, the second-order integral value, and the incremental per-unit value; Correcting the preliminary parameter value using an auxiliary variable method to obtain a target parameter value of the synchronous condenser; Different selections of the excitation voltage reference value result in different steps for obtaining the preliminary parameter values ​​of the synchronous condenser and for correcting the preliminary parameter values; Specifically, when the magnetic field voltage required by the rated stator terminal voltage of the synchronous condenser is the excitation voltage reference value, determining the preliminary parameter value of the synchronous condenser according to the first-order integral value, the second-order integral value, and the incremental per-unit value includes: Determine a first matrix according to the second-order integral value of the excitation voltage increment per unit value, the second-order integral value of the machine terminal voltage increment per unit value, and the second-order integral value of the excitation current increment per unit value; Determine a second matrix according to the per-unit value of the machine-end voltage increment, the first-order integral value of the per-unit value of the machine-end voltage increment, the per-unit value of the excitation current increment, the first-order integral value of the per-unit value of the excitation current increment, the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the armature current increment, the second-order integral value of the per-unit value of the armature current increment, and the first-order integral value of the per-unit value of the excitation voltage increment; Determining a third matrix based on the first matrix and the second matrix; wherein the third matrix represents the preliminary parameter values ​​of the synchronous condenser; The preliminary parameter value is corrected by using an auxiliary variable method to obtain a target parameter value of the synchronous condenser, including: Using a first numerical simulation method to obtain a first time domain response of the per-unit value of the terminal voltage increment and a second time domain response of the per-unit value of the excitation current increment; Obtaining a first-order integral value of the first time domain response and a first-order integral value of the second time domain response; determining a first auxiliary variable matrix according to the first time domain response, the first-order integral value of the first time domain response, the second time domain response, the first-order integral value of the second time domain response, the armature current increment per unit value, the first-order integral value of the armature current increment per unit value, the second-order integral value of the armature current increment per unit value, and the first-order integral value of the excitation voltage increment per unit value; determining a first target matrix according to the first matrix, the first auxiliary variable matrix, and the second matrix, wherein the first target matrix represents target parameter values ​​of the synchronous condenser; When the rated excitation voltage of the synchronous condenser is used as the excitation voltage reference value, determining the preliminary parameter value of the synchronous condenser according to the first-order integral value, the second-order integral value, and the incremental per-unit value includes: Determine a fourth matrix according to the second-order integral value of the per-unit value of the terminal voltage increment and the second-order integral value of the per-unit value of the excitation current increment; Determine a fifth matrix according to the per-unit value of the generator-end voltage increment, the first-order integral value of the per-unit value of the generator-end voltage increment, the per-unit value of the excitation current increment, the first-order integral value of the per-unit value of the excitation current increment, the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the armature current increment, the second-order integral value of the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the excitation voltage increment, and the second-order integral value of the per-unit value of the excitation voltage increment; Determining a sixth matrix based on the fourth matrix and the fifth matrix; wherein the sixth matrix represents the preliminary parameter values ​​of the synchronous condenser; The preliminary parameter value is corrected by using an auxiliary variable method to obtain a target parameter value of the synchronous condenser, including: Using a second numerical simulation method to obtain a third time domain response of the per-unit value of the machine-end voltage increment and a fourth time domain response of the per-unit value of the excitation current increment; determining a second auxiliary variable matrix according to the third time domain response, the first-order integral value of the third time domain response, the fourth time domain response, the first-order integral value of the fourth time domain response, the armature current increment per unit value, the first-order integral value of the armature current increment per unit value, the second-order integral value of the armature current increment per unit value, and the first-order integral value of the excitation voltage increment per unit value; A second target matrix is ​​determined according to the fourth matrix, the second auxiliary variable matrix, and the fifth matrix, where the second target matrix represents target parameter values ​​of the synchronous condenser.

2. The method for estimating synchronous condenser parameters according to claim 1, characterized in that: Preprocessing the operating data to obtain an incremental per-unit value corresponding to the operating data; Obtaining a parameter reference value corresponding to the synchronous condenser; The corresponding incremental per unit value is determined according to the operating data before the load dump test, the operating data after the load dump test, and the parameter reference value.

3. The method for estimating synchronous condenser parameters according to claim 2, characterized in that: The obtaining of the parameter reference value corresponding to the synchronous condenser includes: Determining an armature current reference value, an excitation current reference value, and a terminal voltage reference value according to the rated apparent power and rated voltage of the synchronous condenser; If the magnetic field voltage required to generate the rated stator terminal voltage of the synchronous condenser on the air gap line is known and reliable, the magnetic field voltage is used as the excitation voltage reference value; If the magnetic field voltage required to generate the rated stator terminal voltage of the synchronous condenser on the air gap line is unknown or unreliable, the rated excitation voltage of the synchronous condenser is used as the excitation voltage reference value.

4. A device for estimating parameters of a synchronous condenser, characterized in that: include: an acquisition module, configured to respectively acquire operating data of the synchronous condenser before and after the load shedding; wherein the operating data includes the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser before the load shedding, and the armature current, the terminal voltage, the excitation current, and the excitation voltage of the synchronous condenser after the load shedding; and a processing module, configured to pre-process the operating data to obtain an incremental per-unit value corresponding to the operating data; A first determining module is used to determine a corresponding first-order integral value and a second-order integral value according to the incremental per-unit value; A second determining module is configured to determine a preliminary parameter value of the synchronous condenser according to the first-order integral value, the second-order integral value, and the incremental per-unit value; a correction module, configured to correct the preliminary parameter value using an auxiliary variable method to obtain a target parameter value of the synchronous condenser; The selection of the excitation voltage reference value is different, and the steps executed by the second determination module and the steps executed by the correction module are different; Specifically, when the magnetic field voltage required by the rated stator terminal voltage of the synchronous condenser is the excitation voltage reference value, the second determining module is specifically configured to: Determine a first matrix according to the second-order integral value of the excitation voltage increment per unit value, the second-order integral value of the machine terminal voltage increment per unit value, and the second-order integral value of the excitation current increment per unit value; Determine a second matrix according to the per-unit value of the machine-end voltage increment, the first-order integral value of the per-unit value of the machine-end voltage increment, the per-unit value of the excitation current increment, the first-order integral value of the per-unit value of the excitation current increment, the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the armature current increment, the second-order integral value of the per-unit value of the armature current increment, and the first-order integral value of the per-unit value of the excitation voltage increment; Determining a third matrix based on the first matrix and the second matrix; wherein the third matrix represents the preliminary parameter values ​​of the synchronous condenser; Correction module, specifically used for: Using a first numerical simulation method to obtain a first time domain response of the per-unit value of the terminal voltage increment and a second time domain response of the per-unit value of the excitation current increment; Obtaining a first-order integral value of the first time domain response and a first-order integral value of the second time domain response; determining a first auxiliary variable matrix according to the first time domain response, the first-order integral value of the first time domain response, the second time domain response, the first-order integral value of the second time domain response, the armature current increment per unit value, the first-order integral value of the armature current increment per unit value, the second-order integral value of the armature current increment per unit value, and the first-order integral value of the excitation voltage increment per unit value; determining a first target matrix according to the first matrix, the first auxiliary variable matrix, and the second matrix, wherein the first target matrix represents target parameter values ​​of the synchronous condenser; When the rated excitation voltage of the synchronous condenser is used as the excitation voltage reference value, the second determining module is specifically configured to: Determine a fourth matrix according to the second-order integral value of the per-unit value of the terminal voltage increment and the second-order integral value of the per-unit value of the excitation current increment; Determine a fifth matrix according to the per-unit value of the generator-end voltage increment, the first-order integral value of the per-unit value of the generator-end voltage increment, the per-unit value of the excitation current increment, the first-order integral value of the per-unit value of the excitation current increment, the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the armature current increment, the second-order integral value of the per-unit value of the armature current increment, the first-order integral value of the per-unit value of the excitation voltage increment, and the second-order integral value of the per-unit value of the excitation voltage increment; Determining a sixth matrix based on the fourth matrix and the fifth matrix; wherein the sixth matrix represents the preliminary parameter values ​​of the synchronous condenser; Correction module, specifically used for: Using a second numerical simulation method to obtain a third time domain response of the per-unit value of the machine-end voltage increment and a fourth time domain response of the per-unit value of the excitation current increment; determining a second auxiliary variable matrix according to the third time domain response, the first-order integral value of the third time domain response, the fourth time domain response, the first-order integral value of the fourth time domain response, the armature current increment per unit value, the first-order integral value of the armature current increment per unit value, the second-order integral value of the armature current increment per unit value, and the first-order integral value of the excitation voltage increment per unit value; A second target matrix is ​​determined according to the fourth matrix, the second auxiliary variable matrix, and the fifth matrix, where the second target matrix represents target parameter values ​​of the synchronous condenser.

5. An electronic device, characterized in that: including a memory for storing a computer program; A processor is configured to implement the steps of the method for estimating synchronous condenser parameters according to any one of claims 1 to 4 when executing the computer program.

6. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method for estimating synchronous condenser parameters according to any one of claims 1 to 4 are implemented.

Images