CN112968454A - Parameter setting method of energy storage system and energy storage system - Google Patents
Parameter setting method of energy storage system and energy storage system Download PDFInfo
- Publication number
- CN112968454A CN112968454A CN202110347994.4A CN202110347994A CN112968454A CN 112968454 A CN112968454 A CN 112968454A CN 202110347994 A CN202110347994 A CN 202110347994A CN 112968454 A CN112968454 A CN 112968454A
- Authority
- CN
- China
- Prior art keywords
- loop
- virtual
- reactive
- active
- transfer function
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000004146 energy storage Methods 0.000 title claims abstract description 50
- 238000012546 transfer Methods 0.000 claims abstract description 73
- 230000001360 synchronised effect Effects 0.000 claims abstract description 38
- 230000008878 coupling Effects 0.000 claims abstract description 35
- 238000010168 coupling process Methods 0.000 claims abstract description 35
- 238000005859 coupling reaction Methods 0.000 claims abstract description 35
- 230000001939 inductive effect Effects 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 8
- 239000000470 constituent Substances 0.000 claims 1
- 230000010355 oscillation Effects 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 61
- 238000010586 diagram Methods 0.000 description 16
- 238000003860 storage Methods 0.000 description 15
- 238000004590 computer program Methods 0.000 description 11
- 238000012545 processing Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000013016 damping Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a parameter setting method of an energy storage system and the energy storage system, wherein the energy storage system is merged into a power grid through a virtual synchronous machine control system, and the parameter setting method comprises the following steps: step 1, introducing a virtual impedance loop in the control of a virtual synchronous machine, and determining the composition parameters of the virtual impedance loop according to the power coupling degree; step 2, establishing a reactive loop open-loop transfer function and an active loop open-loop transfer function under the active loop and reactive loop decoupling condition; and 3, determining the inertia coefficient and the rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, the reactive loop open-loop transfer function and the active loop open-loop transfer function. The virtual impedance loop is introduced into the control of the virtual synchronous machine, so that the power oscillation of the energy storage system is effectively avoided, and the parameters of the energy storage system are directly and quantitatively calculated by establishing the open-loop transfer function of the virtual synchronous machine under the decoupling condition of the active loop and the reactive loop.
Description
Technical Field
The invention belongs to the technical field of energy storage control, and particularly relates to a parameter setting method of an energy storage system and the energy storage system.
Background
The conventional energy storage grid-connected inverter is high in response speed, almost has no rotational inertia, is difficult to participate in power grid regulation, cannot provide necessary voltage and frequency support for an active power distribution network containing a distributed power supply, cannot provide necessary damping action for a micro-grid with relatively poor stability, and lacks a mechanism for effectively synchronizing with the distribution network and the micro-grid. Therefore, a student provides a VSG (virtual synchronous generator) control method, so that a grid-connected converter simulates the characteristics of a synchronous generator, and further has primary frequency modulation characteristics and primary voltage regulation characteristics, and the stability of the frequency of a microgrid is enhanced by introducing virtual inertia and damping links, so that an energy storage system can be seamlessly switched no matter in a planned or unplanned island situation, and the normal work of an important load is guaranteed to the maximum extent.
However, since the mathematical model of the VSG control system is a second-order system, there is a pair of poles close to the virtual axis, when the setting of the relevant core parameters of the control system is not reasonable, the effect of the control system may be poor, and the system may oscillate and diverge when the system is severe. The existing scheme for designing the VSG parameters mainly comprises the following steps: 1. the closed-loop parameters of the VSG are designed by a root track method, but the optimized parameters can be obtained only by simulation and repeated trial and error. And the quality of the output voltage and current waveform of the VSG control system designed by the method is difficult to ensure. 2. Considering the VSG control system as a typical second order system, the goal is to ensure that the damping ratio ζ of the system is 0.707. However, the active loop and the reactive loop of the VSG control system each have two pending nodesControl parameters (active loop: damping coefficient D, virtual moment of inertia J; reactive loop: reactive regulation coefficient K, voltage regulation coefficient D)q) Only the constraint of ζ being 0.707 cannot directly derive all the control parameters, and therefore, further adjustment and trial and error are still needed to obtain an optimized result. 3. The open-loop transfer functions of an active loop and a reactive loop of the VSG control system are established, and the active loop and the reactive loop are proved to be approximately decoupled when the line impedance is inductive, and then the performance index is controlled to carry out quantitative setting on the control parameters. However, the method has the precondition that the line impedance is ensured to be inductive, and the line inductance needs to be increased in certain application occasions, so that the cost is high.
Disclosure of Invention
In view of the above problems, the present invention discloses a parameter setting method for an energy storage system and an energy storage system, so as to overcome the above problems or at least partially solve the above problems.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention discloses a parameter setting method of an energy storage system, wherein the energy storage system is merged into a power grid through a virtual synchronous machine control system, and the parameter setting method comprises the following steps:
and 3, determining the inertia coefficient and the rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, the reactive loop open-loop transfer function and the active loop open-loop transfer function.
Further, the determining the composition parameters of the virtual impedance loop according to the power coupling degree specifically includes:
defining the ratio of the partial derivative of active power P to potential E to the partial derivative of active power to power angle delta as the coupling coefficient of reactive power control to active power controlKpDefining the ratio of the partial derivative of the reactive power Q to the power angle delta to the partial derivative of the reactive power Q to the electric potential E as the coupling coefficient K of the active control to the reactive controlqThe calculation formula is as follows:
where s is the differential in the time domain, LvIs a virtual inductor, RvIs a virtual resistance, δ is the active power vs. power angle, XvIs a virtual inductive reactance, E is a potential;
i.e. the power coupling coefficient K1The calculation formula is as follows:
when the power coupling coefficient K1When the minimum value is taken, the virtual resistance R is obtainedvValue of (2) and the virtual inductive reactance XvThe value of (c).
Further, the determining the composition parameters of the virtual impedance loop according to the power coupling degree further includes:
subjecting the virtual inductive reactance XvSubstituting the following formula to calculate the virtual inductance:
wherein, ω isnThe rated angular frequency of the power grid.
Further, the establishing of the reactive loop open-loop transfer function and the active loop open-loop transfer function specifically includes:
establishing the open loop transfer function of the reactive loop:
wherein E isnIs the bridge arm midpoint voltage, UgTo output a voltage amplitude, ZsAs total impedance of the line, DqIs a reactive-voltage droop coefficient, K2Is the coefficient of inertia;
establishing the active loop open loop transfer function:
wherein, ω isnFor rating the angular frequency of the grid, DpThe active-voltage droop coefficient is shown, J is the virtual moment of inertia, and s is the differential of the time domain.
Further, the determining an inertia coefficient and a rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, the reactive loop open-loop transfer function and the active loop open-loop transfer function specifically includes:
and calculating to obtain the total line impedance according to the virtual resistor and the virtual inductive reactance, and substituting the total line impedance into the open loop transfer function of the reactive loop and the open loop transfer function of the active loop to determine the inertia coefficient and the rotational inertia of the virtual synchronous machine.
Further, the total line impedance is calculated according to the following formula:
Zs＝(R+Rv)+j(X+Xv)
wherein R is the actual line equivalent resistance, RvIs a virtual resistor, X is an actual line equivalent inductive reactance, XvIs a virtual inductive reactance.
Further, determining the inertia coefficient specifically includes:
in the reactive open-loop transfer function, let at the cut-off frequency fcqAnd the loop gain amplitude of the reactive loop is 1, then:
and solving to obtain an inertia coefficient:
further, determining the moment of inertia specifically includes:
in the active open-circuit transfer function, let at the cut-off frequency fcpAnd the loop gain amplitude of the reactive loop is 1, then:
solving to obtain:
and determining the value of the moment of inertia J according to the phase angle margin.
Further, the reactive-voltage droop coefficient DqAnd the active-voltage droop coefficient DpCalculated according to the following formula:
P＝Pset+Dp(ωn-ω)
Q＝Qset+Dq(En-Ug)
wherein P is the feedback value of the active power of the system, PsetIs the active power set value of the system, omeganFor the nominal angular frequency of the grid, omega being the angular frequency of the system output, EnIs the midpoint voltage of the bridge arm, Q is the reactive power feedback value of the system, QsetFor a system reactive power setpoint, UgIs the output voltage amplitude.
The invention also discloses an energy storage system, which adopts the parameter setting method to carry out parameter setting.
The invention has the advantages and beneficial effects that:
in the parameter setting method, the virtual impedance loop is introduced into the control of the virtual synchronous machine, so that the power oscillation of the energy storage system is effectively avoided, and the method directly and quantitatively calculates the parameters of the energy storage system by establishing the open-loop transfer function of the virtual synchronous machine under the decoupling condition of the active loop and the reactive loop, and has the advantages of simplicity and strong applicability.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
FIG. 1 is a diagram illustrating steps performed in a method for tuning parameters of an energy storage system according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of VSG control in an energy storage system in accordance with an embodiment of the present invention;
FIG. 3 shows the power coupling coefficient K according to an embodiment of the present invention1A functional relationship diagram of (1);
FIG. 4 is a diagram of a loop small signal structure in accordance with an embodiment of the present invention;
FIG. 5 is a plot of moment of inertia J as a function of time in one embodiment of the present invention;
FIG. 6 is a simulation waveform diagram of a simulation test in an embodiment of the present invention;
FIG. 7 is a schematic diagram of an energy storage system according to an embodiment of the present disclosure;
fig. 8 is a schematic structural diagram of an energy storage device according to an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be described in detail and fully with reference to the accompanying drawings. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention takes a cascade H-bridge energy storage system as a control object, adopts a virtual synchronous generator control strategy (VSG) to control, and aims at the problems existing in the prior scheme: firstly, a virtual impedance loop is introduced in the control of a virtual synchronous machine, then an open-loop transfer function is established under the decoupling condition of an active loop and a reactive loop, and the damping coefficient D, the virtual moment of inertia J and a virtual resistor R in an energy storage system are realizedvVirtual inductor LvQuantitative design was performed in equal amounts.
The technical solutions provided by the embodiments of the present invention are described in detail below with reference to the accompanying drawings.
In an embodiment of the present invention, a method for setting a parameter of an energy storage system is disclosed, the energy storage system is incorporated into a power grid through a virtual synchronous machine control system, as shown in fig. 1, the method for setting the parameter includes the following steps:
And 2, establishing a reactive ring open-loop transfer function and an active ring open-loop transfer function under the active ring and reactive ring decoupling condition.
And 3, determining the inertia coefficient and the rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, the reactive loop open-loop transfer function and the active loop open-loop transfer function.
In summary, in the parameter setting method of this embodiment, the virtual impedance loop is introduced into the control of the virtual synchronous machine, so that power oscillation of the energy storage system is effectively avoided, and the method directly and quantitatively calculates the parameters of the energy storage system by establishing the open-loop transfer function of the virtual synchronous machine under the decoupling condition of the active loop and the reactive loop, and is simple and highly applicable.
In one embodiment, determining the composition parameters of the virtual impedance loop according to the power coupling degree specifically includes:
defining the ratio of the partial derivative of active power P to potential E to the partial derivative of active power to power angle delta as the coupling coefficient K of reactive power control to active power controlpDefining the ratio of the partial derivative of the reactive power Q to the power angle delta to the partial derivative of the reactive power Q to the electric potential E as the coupling coefficient K of the active control to the reactive controlqThe calculation formula is as follows:
where s is the differential in the time domain, LvIs a virtual inductor, RvIs a virtual resistance, δ is the active power vs. power angle, XvIs a virtual inductive reactance, and E is a potential.
I.e. the power coupling coefficient K1The calculation formula is as follows:
wherein R isvIs a virtual resistance, δ is the active power vs. power angle, XvIs a virtual inductive reactance.
FIG. 3 shows the power coupling coefficient K obtained according to the above formula1As can be seen from fig. 3, when the virtual inductive reactance X of the line is reachedvThe larger the virtual resistance RvThe smaller the power coupling coefficient K1The smaller the value of (c), i.e. the smaller the degree of power coupling. When power coupling coefficient K1When the minimum value is taken, obtaining the virtual resistance RvValue of (2) and virtual inductive reactance XvThe value of (c).
Further, determining the composition parameters of the virtual impedance loop according to the power coupling degree, further comprising:
will sense the reactance X virtuallyvSubstituting the following formula to calculate the virtual inductance:
wherein, ω isnIs the grid voltage angular frequency.
In one embodiment, as shown in fig. 4, the establishing of the reactive loop open-loop transfer function and the active loop open-loop transfer function specifically includes:
establishing a reactive loop open-loop transfer function:
wherein E isnIs the bridge arm midpoint voltage, UgTo output a voltage amplitude, ZsAs total impedance of the line, DqIs a reactive-voltage droop coefficient, K2Is the coefficient of inertia and s is the differential in the time domain.
Establishing an active loop open loop transfer function:
wherein, UgTo output a voltage amplitude, ZsAs total impedance of the line, ωnAt a nominal angular frequency, DpThe active-voltage droop coefficient is shown, J is the virtual moment of inertia, and s is the differential of the time domain.
In one embodiment, determining an inertia coefficient and a rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, and the reactive loop open-loop transfer function and the active loop open-loop transfer function specifically includes:
and calculating according to the virtual resistance and the virtual inductive reactance to obtain the total impedance of the line, and introducing the total impedance of the line into a reactive loop open-loop transfer function and an active loop open-loop transfer function to determine the inertia coefficient and the rotational inertia of the virtual synchronous machine.
In one embodiment, the total line impedance is calculated according to the following equation:
Zs＝(R+Rv)+j(X+Xv)
wherein R is the actual line resistance value, RvIs a virtual resistor, X is an actual line inductance, XvIs a virtual inductive reactance.
In one embodiment, determining the coefficient of inertia is embodied as:
in the reactive open-loop transfer function, let at the cut-off frequency fcqAnd the loop gain amplitude of the reactive loop is 1, then:
and solving to obtain an inertia coefficient:
in one embodiment, the determining the moment of inertia is specifically:
in the active open-circuit transfer function, let at the cut-off frequency fcpAnd the loop gain amplitude of the reactive loop is 1, then:
solving to obtain:
the value of the moment of inertia J is determined according to the phase angle margin, fig. 5 is a functional relation graph of the moment of inertia J obtained according to the formula, and it can be seen from fig. 5 that after a virtual impedance loop is introduced in the control of the virtual synchronous machine, the selectable range of the moment of inertia J is far away from 50Hz, and power oscillation near the potential fundamental frequency of the energy storage system is avoided.
The simulation test is performed on the parameter setting method of the energy storage system, as can be seen from fig. 6, the energy storage system is pre-synchronized, then the active power (per unit value pu) is given to be 0.8pu from 2s, after 3.5s, the active power (per unit value pu) is given to be 0.5pu, the corresponding output power and output current are also 0.8pu and 0.5pu, and no oscillation occurs in the transition process, so that the correctness and effectiveness of the method are verified.
In one embodiment, the reactive-voltage droop coefficient DqAnd active-voltage droop coefficient DpCalculated according to the following formula:
P＝Pset+Dp(ωn-ω)
Q＝Qset+Dq(En-Ug)
wherein P is the feedback value of the active power of the system, PsetIs the active power set value of the system, omeganFor the nominal angular frequency of the grid, omega being the angular frequency of the system output, EnIs the midpoint voltage of the bridge arm, Q is the reactive power feedback value of the system, QsetFor a system reactive power setpoint, UgIs the output voltage amplitude.
In an embodiment of the present invention, an energy storage system 700 is disclosed, and the energy storage system 700 performs parameter setting by using the parameter setting method. The energy storage system 700 includes:
the virtual impedance introducing unit 710 is configured to introduce a virtual impedance loop in the control of the virtual synchronous machine, and determine a composition parameter of the virtual impedance loop according to the power coupling degree.
And an open-loop function establishing unit 720, configured to establish a reactive loop open-loop transfer function and an active loop open-loop transfer function under the active loop and reactive loop decoupling condition.
And an inertia coefficient and rotational inertia determining unit 730, configured to determine an inertia coefficient and a rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, and the reactive loop open-loop transfer function and the active loop open-loop transfer function.
In an embodiment, the virtual impedance introducing unit 710 is configured to determine a composition parameter of a virtual impedance loop according to the power coupling degree, and specifically includes:
defining the ratio of the partial derivative of active power P to potential E to the partial derivative of active power to power angle delta as the coupling coefficient K of reactive power control to active power controlpDefining the ratio of the partial derivative of the reactive power Q to the power angle delta to the partial derivative of the reactive power Q to the electric potential E as the coupling coefficient K of the active control to the reactive controlqThe calculation formula is as follows:
where s is the differential in the time domain, LvIs a virtual inductor, RvIs a virtual resistance, δ is the active power vs. power angle, XvIs a virtual inductive reactance, and E is a potential.
I.e. the power coupling coefficient K1The calculation formula is as follows:
when power coupling coefficient K1When the minimum value is taken, obtaining the virtual resistance RvValue of (2) and virtual inductive reactance XvThe value of (c).
In one embodiment, the virtual impedance introducing unit 710 is configured to determine a composition parameter of a virtual impedance loop according to the power coupling degree, and further includes:
will sense the reactance X virtuallyvSubstituting the following formula to calculate the virtual inductance:
wherein, ω isnThe rated angular frequency of the power grid.
In an embodiment, the open-loop function establishing unit 720 is configured to establish a reactive loop open-loop transfer function and an active loop open-loop transfer function specifically as follows:
establishing a reactive loop open-loop transfer function:
wherein E isnIs the bridge arm midpoint voltage, UgTo output a voltage amplitude, ZsAs total impedance of the line, DqIs a reactive-voltage droop coefficient, K2Is the coefficient of inertia.
Establishing an active loop open loop transfer function:
wherein, ω isnAt a nominal angular frequency, DpThe active-voltage droop coefficient is shown, J is the virtual moment of inertia, and s is the differential of the time domain.
In an embodiment, the inertia coefficient and rotational inertia determining unit 730 is configured to determine the inertia coefficient and the rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, and the reactive loop open-loop transfer function and the active loop open-loop transfer function, and specifically includes:
and calculating according to the virtual resistance and the virtual inductive reactance to obtain the total impedance of the line, and introducing the total impedance of the line into a reactive loop open-loop transfer function and an active loop open-loop transfer function to determine the inertia coefficient and the rotational inertia of the virtual synchronous machine.
Further, the total line impedance is calculated according to the following formula:
Zs＝(R+Rv)+j(X+Xv)
wherein R is the actual line resistance, RvIs a virtual resistor, X is an actual line inductance, XvIs a virtual inductive reactance.
In one embodiment, determining the coefficient of inertia is embodied as:
in the reactive open-loop transfer function, let be at cut-offFrequency fcqAnd the loop gain amplitude of the reactive loop is 1, then:
and solving to obtain an inertia coefficient:
in one embodiment, the determining the moment of inertia is specifically:
in the active open-circuit transfer function, let at the cut-off frequency fcpAnd the loop gain amplitude of the reactive loop is 1, then:
solving to obtain:
and determining the value of the moment of inertia J according to the phase angle margin.
In one embodiment, the reactive-voltage droop coefficient DqAnd active-voltage droop coefficient DpCalculated according to the following formula:
P＝Pset+Dp(ωn-ω)
Q＝Qset+Dq(En-Ug)
wherein P is the feedback value of the active power of the system, PsetIs the active power set value of the system, omeganFor the nominal angular frequency of the grid, omega being the angular frequency of the system output, EnIs the midpoint voltage of the bridge arm, Q is the reactive power feedback value of the system, QsetFor a system reactive power setpoint, UgIs the output voltage amplitude.
One embodiment of the present invention discloses an energy storage device, including: a processor; and a memory arranged to store computer executable instructions that, when executed, cause the processor to perform any of the energy storage system parameter tuning methods described above.
Fig. 8 is a schematic structural diagram of an energy storage device according to an embodiment of the present application. Referring to fig. 8, in the hardware level, the energy storage device includes a processor, and optionally an internal bus, a network interface, and a memory. The Memory may include a Memory, such as a Random-Access Memory (RAM), and may further include a non-volatile Memory, such as at least 1 disk Memory. Of course, the energy storage device may also include hardware required for other services.
The processor, the network interface, and the memory may be connected to each other via an internal bus, which may be an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus, an EISA (Extended Industry Standard Architecture) bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one double-headed arrow is shown in FIG. 8, but that does not indicate only one bus or one type of bus.
And the memory is used for storing programs. In particular, the program may include program code including computer operating instructions. The memory may include both memory and non-volatile storage and provides instructions and data to the processor.
The processor reads the corresponding computer program from the nonvolatile memory into the memory and then runs the computer program to form the target detection device on a logic level. The processor is used for executing the program stored in the memory and is specifically used for executing the following operations:
introducing a virtual impedance loop in the control of the virtual synchronous machine, and determining the composition parameters of the virtual impedance loop according to the power coupling degree; establishing a reactive loop open-loop transfer function and an active loop open-loop transfer function under the active loop and reactive loop decoupling condition; and determining the inertia coefficient and the rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, the open-loop transfer function of the reactive loop and the open-loop transfer function of the active loop.
The parameter tuning method disclosed in the embodiment of fig. 1 of the present application may be applied to a processor, or implemented by the processor. The processor may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware in a processor or instructions in the form of software. The Processor may be a general-purpose Processor, including a Central Processing Unit (CPU), a Network Processor (NP), and the like; but also Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field Programmable Gate Arrays (FPGAs) or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components. The various methods, steps, and logic blocks disclosed in the embodiments of the present application may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present application may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The storage medium is located in a memory, and a processor reads information in the memory and completes the steps of the method in combination with hardware of the processor.
An embodiment of the present application further provides a computer-readable storage medium storing one or more programs, where the one or more programs include instructions, which when executed by an energy storage device including a plurality of application programs, enable the energy storage device to perform a parameter setting method of an energy storage system in the embodiment shown in fig. 1, and are specifically configured to perform:
introducing a virtual impedance loop in the control of the virtual synchronous machine, and determining the composition parameters of the virtual impedance loop according to the power coupling degree; establishing a reactive loop open-loop transfer function and an active loop open-loop transfer function under the active loop and reactive loop decoupling condition; and determining the inertia coefficient and the rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, the open-loop transfer function of the reactive loop and the open-loop transfer function of the active loop.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include transitory computer readable media (transmyedia) such as modulated data signals and carrier waves.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in the process, method, article, or apparatus that comprises the element.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
While the foregoing is directed to embodiments of the present invention, other modifications and variations of the present invention may be devised by those skilled in the art in light of the above teachings. It should be understood by those skilled in the art that the foregoing detailed description is for the purpose of better explaining the present invention, and the scope of the present invention should be determined by the scope of the appended claims.
Claims (10)
1. A parameter setting method of an energy storage system, wherein the energy storage system is merged into a power grid through a virtual synchronous machine control system, is characterized by comprising the following steps:
step 1, introducing a virtual impedance loop in the control of a virtual synchronous machine, and determining the composition parameters of the virtual impedance loop according to the power coupling degree;
step 2, establishing a reactive loop open-loop transfer function and an active loop open-loop transfer function under the active loop and reactive loop decoupling condition;
and 3, determining the inertia coefficient and the rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, the reactive loop open-loop transfer function and the active loop open-loop transfer function.
2. The parameter tuning method according to claim 1, wherein the determining the composition parameters of the virtual impedance loop according to the power coupling degree specifically includes:
defining the ratio of the partial derivative of active power P to potential E to the partial derivative of active power to power angle delta as the coupling coefficient K of reactive power control to active power controlpDefining the ratio of the partial derivative of the reactive power Q to the power angle delta to the partial derivative of the reactive power Q to the electric potential E as the coupling coefficient K of the active control to the reactive controlqThe calculation formula is as follows:
where s is the differential in the time domain, LvIs a virtual inductor, RvIs a virtual resistance, δ is the active power vs. power angle, XvIs a virtual inductive reactance, E is a potential;
i.e. the power coupling coefficient K1The calculation formula is as follows:
when the power coupling coefficient K1When the minimum value is taken, the virtual resistance R is obtainedvValue of (2) and the virtual inductive reactance XvThe value of (c).
3. The parameter tuning method of claim 2, wherein the determining the constituent parameters of the virtual impedance loop according to the degree of power coupling further comprises:
subjecting the virtual inductive reactance XvSubstituting the following formula to calculate the virtual inductance:
wherein, ω isnThe rated angular frequency of the power grid.
4. The parameter setting method according to claim 2, wherein the establishing of the reactive loop open-loop transfer function and the active loop open-loop transfer function is specifically:
establishing the open loop transfer function of the reactive loop:
wherein E isnIs the bridge arm midpoint voltage, UgTo output a voltage amplitude, ZsAs total impedance of the line, DqIs a reactive-voltage droop coefficient, K2Is the coefficient of inertia;
establishing the active loop open loop transfer function:
wherein, ω isnFor rating the angular frequency of the grid, DpThe active-voltage droop coefficient is shown, J is the virtual moment of inertia, and s is the differential of the time domain.
5. The parameter setting method according to claim 4, wherein the determining an inertia coefficient and a rotational inertia of the virtual synchronous machine according to the composition parameters of the virtual impedance loop, the reactive loop open-loop transfer function and the active loop open-loop transfer function specifically comprises:
and calculating to obtain the total line impedance according to the virtual resistor and the virtual inductive reactance, and substituting the total line impedance into the open loop transfer function of the reactive loop and the open loop transfer function of the active loop to determine the inertia coefficient and the rotational inertia of the virtual synchronous machine.
6. The parameter tuning method of claim 5, wherein the total line impedance is calculated according to the following formula:
Zs＝(R+Rv)+j(X+Xv)
wherein R is the actual line equivalent resistance, RvIs a virtual resistor, X is an actual line inductance, XvIs a virtual inductive reactance.
7. The parameter tuning method according to claim 6, wherein determining the inertia coefficient specifically comprises:
in the reactive open-loop transfer function, let at the cut-off frequency fcqAnd the loop gain amplitude of the reactive loop is 1, then:
and solving to obtain an inertia coefficient:
8. the parameter tuning method according to claim 6, wherein determining the moment of inertia specifically comprises:
in the active open-circuit transfer function, let at the cut-off frequency fcpAnd the loop gain amplitude of the reactive loop is 1, then:
solving to obtain:
and determining the value of the moment of inertia J according to the phase angle margin.
9. The parameter tuning method of claim 8, wherein the reactive-voltage droop coefficient DqAnd the active-voltage droop coefficient DpCalculated according to the following formula:
P＝Pset+Dp(ωn-ω)
Q＝Qset+Dq(En-Ug)
wherein P is the feedback value of the active power of the system, PsetIs the active power set value of the system, omeganFor the nominal angular frequency of the grid, omega being the angular frequency of the system output, EnIs the midpoint voltage of the bridge arm, Q is the reactive power feedback value of the system, QsetFor a system reactive power setpoint, UgIs the output voltage amplitude.
10. An energy storage system, characterized in that the energy storage system is parameter-tuned using the parameter-tuning method of any of claims 1-9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110347994.4A CN112968454A (en) | 2021-03-31 | 2021-03-31 | Parameter setting method of energy storage system and energy storage system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110347994.4A CN112968454A (en) | 2021-03-31 | 2021-03-31 | Parameter setting method of energy storage system and energy storage system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112968454A true CN112968454A (en) | 2021-06-15 |
Family
ID=76280571
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110347994.4A Pending CN112968454A (en) | 2021-03-31 | 2021-03-31 | Parameter setting method of energy storage system and energy storage system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112968454A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113469454A (en) * | 2021-07-20 | 2021-10-01 | 华能陇东能源有限责任公司 | Integrated energy base direct current modulation parameter configuration method and system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3376627A1 (en) * | 2017-03-14 | 2018-09-19 | ABB Schweiz AG | Method and control system for controlling a power converter |
| CN112290603A (en) * | 2019-07-22 | 2021-01-29 | 华北电力大学（保定） | VSG power decoupling control method of virtual steady-state synchronous negative impedance |
-
2021
- 2021-03-31 CN CN202110347994.4A patent/CN112968454A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3376627A1 (en) * | 2017-03-14 | 2018-09-19 | ABB Schweiz AG | Method and control system for controlling a power converter |
| CN112290603A (en) * | 2019-07-22 | 2021-01-29 | 华北电力大学（保定） | VSG power decoupling control method of virtual steady-state synchronous negative impedance |
Non-Patent Citations (3)
| Title |
|---|
| 吴恒等: ""虚拟同步发电机功率 环的建模与参数设计"", 《中国电机工程学报》, vol. 35, no. 24, pages 6508 - 6518 * |
| 吴恒等: "虚拟同步发电机功率环的建模与参数设计", 《中国电机工程学报》 * |
| 吴恒等: "虚拟同步发电机功率环的建模与参数设计", 《中国电机工程学报》, vol. 35, no. 24, 20 December 2015 (2015-12-20), pages 6508 - 6518 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113469454A (en) * | 2021-07-20 | 2021-10-01 | 华能陇东能源有限责任公司 | Integrated energy base direct current modulation parameter configuration method and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9379546B2 (en) | Vector control of grid-connected power electronic converter using artificial neural networks | |
| Avelar et al. | A state equation model of a single-phase grid-connected inverter using a droop control scheme with extra phase shift control action | |
| CN108462206B (en) | Method and device for determining selectable range of virtual inertia and damping coefficient of VSG | |
| US20200081044A1 (en) | Power flow calculation method and device for ac-dc interconnected power system, storage medium and terminal | |
| Gonzalez-Cajigas et al. | Design and analysis of parallel-connected grid-forming virtual synchronous machines for island and grid-connected applications | |
| Neumann et al. | Modelling and control of photovoltaic inverter systems with respect to German grid code requirements | |
| CN109638881B (en) | Virtual synchronization method and system for energy storage inverter with power grid strength self-adaptive optimization | |
| CN113595123B (en) | High-frequency impedance calculation method and device for modularized multi-level converter | |
| Lacerda et al. | Phasor modeling approaches and simulation guidelines of voltage-source converters in grid-integration studies | |
| CN112968454A (en) | Parameter setting method of energy storage system and energy storage system | |
| Ebrahim et al. | Optimal design of proportional‐resonant controller and its harmonic compensators for grid‐integrated renewable energy sources based three‐phase voltage source inverters | |
| Oraa et al. | Modeling of a droop-controlled grid-connected DFIG wind turbine | |
| Malkhandi et al. | Return ratio shaping approach to stabilize inverter-weak grid system | |
| CN112736896B (en) | Three-stage initialization method and system suitable for full-electromagnetic transient simulation | |
| Sigrist et al. | An induction motor model for system frequency response models | |
| Chen et al. | High penetration of inverter-based power sources with VSG control impact on electromechanical oscillation of power system | |
| Feng et al. | Damping forced oscillations in power system via interline power flow controller with additional repetitive control | |
| CN111190074B (en) | Power grid synchronous detection method based on single-phase-locked loop | |
| CN112928764B (en) | Parallel energy storage system control method, parallel energy storage system and device | |
| CN116073437A (en) | Three-phase grid-connected inverter droop control strategy based on Lyapunov equation | |
| Liao et al. | Enhanced compensation strategy of voltage harmonics for doubly-fed induction generator | |
| Jiang et al. | Secondary frequency regulation scheme based on improved virtual synchronous generator in an islanded microgrid | |
| Serra et al. | Application of the modified IDA‐PBC for shunt active power filters control | |
| Urtasun et al. | Dynamic Enhancement of the Droop Control for Grid-Forming Inverters | |
| CN117269654B (en) | Method and device for testing and evaluating phase jump active response characteristics of network-structured converter |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210615 |
|
| RJ01 | Rejection of invention patent application after publication |