Power grid adjustment method based on loads of variable frequency air conditioner

Abstract

Provided is a power grid adjustment method based on a load of a variable frequency air conditioner. The method includes: establishing a mathematical simulation model of a virtual synchronous motor in a variable frequency air conditioner controller; establishing a virtual inertia control segment and a droop control segment of power grid adjustment according to the mathematical simulation model; obtaining a reference value Δω of rotation speed variation of a compressor by the virtual inertia control segment; obtaining a reference value ωref of rotational angular frequency of the compressor by the droop control segment; and inputting a sum of the reference value ωref of rotational angular frequency of the compressor and the reference value Δω of rotation speed variation of the compressor into a field-oriented controller (FOC) to control rotation of a motor.

Description

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims priority to Chinese patent applications No. 201810262888.4 filed on Mar. 28, 2018 and No. 201910107263.5 field on Feb. 2, 2019, disclosure of which is incorporated herein by reference in its entirety.

TECHNICAL FIELD

Embodiments of the present disclosure relate to the technical field of a demand side response (DSR) of a power grid, and particularly relate to a power grid adjustment method based on a load of a variable frequency air conditioner.

BACKGROUND

With access of more and more wind and photovoltaic power generation and electric vehicles, a randomness of power supplies and loads in the power grid is greatly increased so that traditional means for frequency modulation and peak regulation cannot meet demands. Meanwhile, the development of an Internet of Things (IOT) technology and a smart home technology gives a two-way communication capacity to the power grid and the loads. Therefore, a DSR is regarded as one of the most promising means for frequency modulation in the future.

In many controllable loads, a load of an air conditioner, especially a load of a central air conditioner, has a greatest potential to participate in the peak regulation and the frequency modulation. Firstly, the electricity consumption of air-conditioner loads in a city is huge and may account for 30%-40% of a peak load at the peak of electricity consumption. Secondly, a closed building space has relatively strong heat and cold storage capacities so that a comfort level of users will not be affected when adjusting the power of the air conditioner in a short time. Thirdly, a response speed of the load of the air conditioner mainly depends on transmission of a control signal, and test results of an existing engineering indicate that the response speed meets requirements of the frequency modulation and peak regulation.

SUMMARY

The present disclosure provides a power grid adjustment method based on a load of a variable frequency air conditioner, so as to solve a problem of insufficient moment of inertia of the power grid caused after a new energy resource accesses to the power grid. A control strategy for the load of the variable frequency air conditioner may be optimized by using the DSR of the power grid, the moment of inertia of a microgrid may be improved, and a fluctuation range of an alternating current (AC) frequency may be reduced when the load fluctuates.

In a first aspect, embodiments of the present disclosure provide a power grid adjustment method based on a load of a variable frequency air conditioner, including:

establishing a mathematical simulation model of a virtual synchronous motor in the variable frequency air conditioner controller;

establishing a virtual inertia control segment and a droop control segment of power grid adjustment according to the mathematical simulation model;

obtaining a reference value Δω of a rotation speed variation of a compressor by the virtual inertia control segment;

obtaining a reference value ω_ref of rotational angular frequency of the compressor by the droop control segment; and

inputting a sum of the reference value ω_ref of rotational angular frequency of the compressor and the reference value Δω of rotation speed variation of the compressor into a field-oriented controller (FOC) to control rotation of a motor.

In an embodiment, the power grid adjustment method further includes:

establishing a signal model of the virtual synchronous motor, and obtaining a cutoff frequency ω of each variable frequency air conditioner;

acquiring a phase margin (PM) of the virtual synchronous motor, and calculating a moment of inertia J of the virtual synchronous motor according to the PM and the cutoff frequency ω of the variable frequency air conditioner; and

obtaining a maximum value of the moment of inertia J of the virtual synchronous motor according to a value range of the PM.

In an embodiment, the mathematical simulation model is

$\{\begin{array}{c}J\frac{d\Delta {\omega }_{\mathrm{sg}}}{\mathrm{dt}}=\frac{-K_D\Delta {\omega }_{\mathrm{sg}}-\Delta P_{\mathrm{sg}}}{{\omega }_N}\\ \Delta \delta =\int \left(\Delta {\omega }_{\mathrm{sg}}-\Delta {\omega }_g\right)\mathrm{dt}\end{array}.$

Where Δω_sg is an angular frequency disturbance quantity of the virtual synchronous motor, Δω_g is an angular frequency disturbance quantity of a voltage of the power grid, K_D is a droop coefficient, ω_N is a rated rotation speed, Δδ is a power angle disturbance quantity, and ΔP_sg is a power disturbance quantity.

In an embodiment, a step of obtaining the reference value Δω of rotation speed variation of the compressor by the virtual inertia control segment includes:

acquiring the voltage U_g of the power grid, generating a virtual voltage orthogonal to the voltage of the power grid by a differential method, forming an αβ static orthogonal coordinate system, and performing Park transformation with a phase angle θ_v of a virtual electromotive force vector of the virtual synchronous motor to obtain components of the voltage of the power grid relative to V_d and V_q, where a power angle δ_v of the virtual synchronous motor is:δ_v=arctan(V_q/V_d).

A virtual power P_v of the virtual synchronous motor is:

$P_v=\frac{U_{g}E}{X_a}\sin {\delta }_v\approx \frac{U_{g}E}{X_a}{\delta }_v.$

The reference value Δω of rotation speed variation of the compressor is:

$\Delta \omega =\frac{1}{n}\Delta P_v.$

Wherein U_g is the voltage of the power grid; E is the virtual electromotive force, X_a is a virtual reactance, n is a constant, and ΔP_v is the amount of power consumption variation.

In an embodiment, a step of obtaining the reference value ω_ref of rotational angular frequency of the compressor by the droop control segment includes:

acquiring a constant operating power P_N of the variable frequency air conditioner according to a set temperature of the variable frequency air conditioner and an indoor temperature, and calculating by a formula P_D=P_N+K_D(ω_N−ω_g) to obtain a steady-state output power P_D of the variable frequency air conditioner, where ω_g is an angular frequency of the voltage of the power grid, and K_D is a droop coefficient of the power grid; and

inputting a difference between an actually measured power of the variable frequency air conditioner and the steady-state output power P_D of the variable frequency air conditioner into a PI controller to obtain the reference value ω_ref of rotational angular frequency of the compressor.

In an embodiment, a step of obtaining the reference value ω_ref of rotational angular frequency of the compressor by the droop control segment further includes:

setting a range of the power of the variable frequency air conditioner as {0, P_max}, then an adjustable range of the power is P_N±min{(P_max−P_N), PN} during a frequency modulation of the power grid, and P_max is a maximum power of the variable frequency air conditioner, and the power of the variable frequency air conditioner reaches an upper limit when a frequency error is 0.3 Hz; and

calculating the droop coefficient K_D of the power grid by a formula:

$K_D=\frac{0.8\min \left\{\left(P_{\max }-P_N\right),P_N\right\}}{0.3\times 2\pi }.$

In an embodiment, the power grid adjustment method further includes:

${\omega }_c=\sqrt{\sqrt{\frac{K_p^2}{J^2{\omega }_N^2}+\frac{K_D^4}{4J^4}}-\frac{K_D^4}{2J^2}},{\omega }_c\le \frac{K_p}{{\omega }_N{K}_D},$

obtaining the cutoff frequency ω_c of each variable frequency air conditioner by the above formula, where K_p is a constant;

$\mathrm{PM}=180°+{\mathrm{\angle T}}_p\left({\omega }_c\right)=90°-\mathrm{arc}\tan \left(\frac{{\omega }_{c}J}{K_D}\right)$

calculating to obtain a maximum value of the moment of inertia J according to the formula of the PM when the PM is in the range of 30° to 50°, wherein ∠T_p(ω_c) is a loop gain amplitude of a virtual inertia stage.

In an embodiment, a step of inputting the sum of the reference value ω_ref of the rotational angular frequency of the compressor and the reference value Δω of rotation speed variation of the compressor into the FOC to control the rotation of the motor includes:

collecting two-phase currents i_a and i_b;

performing clarke transformation on the two-phase currents i_a and i_b to obtain magnitudes of two-axis orthogonal currents i_α and i_β;

performing rotation transformation on the orthogonal current magnitudes i_α and i_β to obtain orthogonal current magnitudes i_d and i_q, where i_q is related to a torque, and i_d is related to a magnetic flux;

inputting the obtained i_d and i_q respectively into the PI controller to obtain corresponding outputs v_d and v_q;

obtaining a rotation angle θ of the motor by a sensor, and performing an inverse park transformation with v_d and v_q to obtain magnitudes of two-axis currents v_α and v_β; and

performing an inverse clarke transformation on v_α and v_β to obtain a required three-phase voltage, and inputting the three-phase voltage to an inverter bridge to drive the motor to rotate.

In a second aspect, embodiments of the present disclosure further provide a power grid adjustment method based on a load of a variable frequency air conditioner, including:

establishing a mathematical simulation model of a virtual synchronous motor of a variable frequency air conditioner;

establishing an adjustment structure of the virtual synchronous motor according to the mathematical simulation model of the virtual synchronous motor;

acquiring parameters of a power grid in real time;

inputting the parameters of the power grid into the adjustment structure for calculation to obtain a reference value Δω of a rotation speed variation of the variable frequency air conditioner and a reference value ω_ref of a rotational angular frequency which is needed to be changed; and

adjusting a working point of the variable frequency air conditioner according to the reference value Δω of the rotation speed variation of the variable frequency air conditioner and the reference value ω_ref of the rotational angular frequency.

In an embodiment, before the step of establishing the mathematical simulation model of the virtual synchronous motor of the variable frequency air conditioner, the method further includes:

establishing a signal model of the virtual synchronous motor, and obtaining a cutoff frequency ω_c of each variable frequency air conditioner;

acquiring a phase margin (PM) of the virtual synchronous motor, and calculating a moment of inertia J of the virtual synchronous motor according to the PM and the cutoff frequency ω_c of the variable frequency air conditioner; and

obtaining a maximum value of the moment of inertia J of the virtual synchronous motor according to a value range of the PM.

In an embodiment, the mathematical simulation model is:

$\{\begin{array}{c}J\frac{d\Delta {\omega }_{\mathrm{sg}}}{\mathrm{dt}}=\frac{-K_D\Delta {\omega }_{\mathrm{sg}}-\Delta P_{\mathrm{sg}}}{{\omega }_N}\\ \Delta \delta =\int \left(\Delta {\omega }_{\mathrm{sg}}-\Delta {\omega }_g\right)\mathrm{dt}\end{array}.$

Where Δω_sg is an angular frequency disturbance quantity of the virtual synchronous motor, Δω_g is the angular frequency disturbance quantity of a voltage of the power grid, K_D is a droop coefficient, ω_N is a rated rotation speed, Δδ is a power angle disturbance quantity, ΔP_sg is a power disturbance quantity, and J is the moment of inertia of the virtual synchronous motor.

In an embodiment, the adjustment structure includes a virtual inertia control segment and a droop control segment, and the step of establishing the adjustment structure of the virtual synchronous motor according to the mathematical simulation model of the virtual synchronous motor includes:

determining an inertia adjustment proportion parameter of the virtual inertia control segment and a droop adjustment proportion parameter of the droop control segment according to the mathematical simulation model of the virtual synchronous motor.

In an embodiment, a step of inputting the parameters of the power grid into the adjustment structure for calculation to obtain the reference value Δω of the rotation speed variation of the variable frequency air conditioner and the reference value ω_ref of the rotational angular frequency which is needed to be changed includes:

inputting a voltage of the power grid into the virtual inertia control segment for calculation to obtain the reference value Δω of rotation speed variation of the variable frequency air conditioner; and

inputting a frequency of the power grid into the droop control segment for calculation to obtain the reference value ω_ref of rotational angular frequency which is needed to be changed.

In an embodiment, the step of inputting the voltage of the power grid into the virtual inertia control segment for calculation to obtain the reference value Δω of the rotation speed variation of the variable frequency air conditioner includes:

acquiring the voltage U_g of the power grid and generating a virtual voltage orthogonal to the voltage of the power grid by a differential method to form an αβ static orthogonal coordinate system, and performing park transformation with a phase angle θ_v of a virtual electromotive force vector of the virtual synchronous motor to obtain components of the voltage of the power grid relative to V_d and V_q, where a power angle δ_v of the virtual synchronous motor isδ_v=arctan(V_q/V_d).

A virtual power P_v of the virtual synchronous motor is:

$P_v=\frac{U_{g}E}{X_a}\sin {\delta }_v\approx \frac{U_{g}E}{X_a}{\delta }_v.$

The reference value Δω of rotation speed variation of the compressor is:

$\Delta \omega =\frac{1}{n}\Delta P_v.$

Where U_g is the voltage of the power grid; E is the virtual electromotive force, X_a is a virtual reactance, n is a constant, and ΔP_v is the amount of power consumption variation.

In an embodiment, the step of inputting the frequency of the power grid into the droop control segment for calculation to obtain the reference value ω_ref of the rotational angular frequency which is needed to be changed includes:

acquiring a constant operating power P_N of the variable frequency air conditioner according to a set temperature of the variable frequency air conditioner and an indoor temperature, and calculating by a formula P_D=P_N±K_D(ω_N−ω_g) to obtain a steady-state output power P_D of the variable frequency air conditioner, where ω_g is an angular frequency of the voltage of the power grid, and K_D is a droop coefficient of the power grid; and

inputting a difference between an actually measured power of the variable frequency air conditioner and the steady-state output power P_D of the variable frequency air conditioner into a PI controller to obtain the reference value ω_ref of the rotational angular frequency of a compressor.

In an embodiment, before inputting the frequency of the power grid into the droop control segment for calculation to obtain the reference value ω_ref of the rotational angular frequency which is needed to be changed, the method further includes: acquiring the droop coefficient of the droop control segment; the droop coefficient K_D is:

$K_D=\frac{0.8\min \left\{\left(P_{\max }-P_N\right),P_N\right\}}{0.3\times 2\pi }.$

Where P_max is a maximum power of the variable frequency air conditioner, P_N is a rated power of the variable frequency air conditioner, and an adjustable range of the power of the variable frequency air conditioner is P_N±min {(P_max−P_N), P_N}.

In an embodiment, the cutoff frequency ω_c of the variable frequency air conditioner is:

${\omega }_c=\sqrt{\sqrt{\frac{K_p^2}{J^2{\omega }_N^2}+\frac{K_D^4}{4J^4}}-\frac{K_D^2}{2J^2}},{\omega }_2\le \frac{K_p}{{\omega }_N{K}_D},$

where K_p is a constant.

In an embodiment, a relationship between a value range of the PM and the maximum value of the moment of inertia J of the virtual synchronous motor is

$\mathrm{PM}=180°+\angle T_p\left({\omega }_c\right)=90°-\mathrm{arc}\tan \left(\frac{{\omega }_{c}J}{K_D}\right).$

wherein PM is the phase margin, and ∠T_p(ω_c) is a loop gain amplitude of a virtual inertia stage.

In an embodiment, the step of adjusting the working point of the variable frequency air conditioner according to the reference value Δω of the rotation speed variation of the variable frequency air conditioner and the reference value ω_ref of the rotational angular frequency includes:

collecting two-phase currents i_a and i_b of the variable frequency air conditioner;

performing clarke transformation on the two-phase currents i_a and i_b to obtain magnitudes of two-axis orthogonal currents i_α and i_β;

performing rotation transformation on the magnitudes of the orthogonal currents i_a and i_β to obtain orthogonal current magnitudes i_d and i_q, where, i_q is related to a torque, and i_d is related to a magnetic flux;

inputting the obtained i_d and i_q respectively into the PI controller to obtain corresponding outputs v_d and v_q;

obtaining a rotation angle θ of the motor of the variable frequency air conditioner by a sensor, and performing an inverse park transformation with v_d and v_q to obtain magnitudes of two-axis currents v_α and v_β; and

performing an inverse clarke transformation on v_α and v_β to obtain a required three-phase voltage, and inputting the three-phase voltage to an inverter bridge to drive the variable frequency air conditioner to work.

According to a technical solution of the present disclosure, the output power of the variable frequency air conditioner can simulate variation processes of the output power and the rotation speed of a no-load synchronous motor with the same capacity by establishing the mathematical simulation model of the virtual synchronous motor when the frequency of the power grid fluctuates. The technical solution can be adopted to solve the problem of insufficient moment of inertia of the power grid caused after accessing a new energy source to the power grid, improve the moment of inertia of the microgrid, and reduce the fluctuation range of AC frequency when the load fluctuates.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 is a schematic diagram illustrating a power grid adjustment method based on a load of a variable frequency air conditioner provided by an embodiment of the present disclosure;

FIG. 2 is a schematic diagram illustrating a control structure of a variable frequency air conditioner based on a virtual motor provided by an embodiment of the present disclosure;

FIG. 3 is a schematic diagram illustrating simulation results of a total power of a single group of variable frequency air conditioners under different control methods provided by an embodiment of the present disclosure;

FIG. 4 is a schematic diagram illustrating simulation results of frequencies of a power grid under different control methods provided by an embodiment of the present disclosure;

FIG. 5 is a schematic diagram illustrating a structure of a field-oriented controller (FOC) provided by an embodiment of the present disclosure; and

FIG. 6 is a schematic diagram illustrating another power grid adjustment method based on the load of the variable frequency air conditioner provided by an embodiment of the present disclosure.

DETAILED DESCRIPTION

The present disclosure will be further described below in detail in combination with drawings and embodiments. It can be understood that specific embodiments described herein are only used for explaining the present disclosure, not used for limiting the present disclosure. In addition, it shall be illustrated that for ease of description, drawings only show some structures related to the present disclosure rather than all structures.

In view of a problem of insufficient moment of inertia of a power grid caused after a new energy resource is accessed to the power grid, the present disclosure provides a power grid adjustment method based on a load of a variable frequency air conditioner by optimizing a control strategy for the load of the variable frequency air conditioner with a DSR of the power grid. By establishing a mathematical simulation model of a virtual synchronous motor, an output power of the variable frequency air conditioner can simulate an output power and a change process of a rotation speed of a no-load synchronous motor with the same capacity when frequencies of the power grid fluctuates. The power grid adjustment method can improve a moment of inertia of a microgrid and reduce a fluctuation range of AC frequency when a load fluctuates.

As shown in FIG. 1, which is a schematic diagram illustrating a power grid adjustment method based on a load of a variable frequency air conditioner provided by the present disclosure, the method includes the following steps.

In S1, a mathematical simulation model of a virtual synchronous motor is established in a variable frequency air conditioner controller.

In S2, a virtual inertia control segment and a droop control segment of the variable frequency air conditioner is established according to the mathematical simulation model.

In S3, a reference value Δω of rotation speed variation of a compressor is obtained by the virtual inertia control segment.

In S4, a reference value ω_ref of rotational angular frequency of the compressor is obtained by the droop control segment.

In S5, a sum of the reference value ω_ref of rotational angular frequency of the compressor and the reference value Δω of rotation speed variation of the compressor is inputted into a field-oriented controller (FOC) to control rotation of the motor.

The method further includes the following steps. A signal model of the virtual synchronous motor is established, and a cutoff frequency ω_c of each variable frequency air conditioner is obtained. A phase margin (PM) of the virtual synchronous motor is acquired, and a moment of inertia J of the virtual synchronous motor is calculated according to the PM and the cutoff frequency ω_c of the variable frequency air conditioner. A maximum value of the moment of inertia J of the virtual synchronous motor is obtained according to a value range of the PM.

As shown in FIG. 2, it is a schematic diagram illustrating a control structure of the variable frequency air conditioner based on the virtual motor provided by the present disclosure. a voltage U_g of the power grid is collected when a frequency of the power grid changes. The voltage of the power grid is subjected to the virtual inertia control segment to obtain the reference value Δω of rotation speed variation of the compressor, and meanwhile, the voltage of the power grid is subjected to the droop control segment to obtain the reference value ω_ref of rotational angular frequency of the compressor. Then, the sum of the two is inputted into the FOC for control. A fluctuation range of AC frequencies is reduced when the load fluctuates, and the moment of inertia of the power grid is improved.

The mathematical simulation model is:

$\hspace{1em}\{\begin{array}{c}J\frac{d\Delta {\omega }_{\mathrm{sg}}}{\mathrm{dt}}=\frac{-K_D\Delta {\omega }_{\mathrm{sg}}-\Delta P_{\mathrm{sg}}}{{\omega }_N}\\ \Delta \delta =\int \left(\Delta {\omega }_{\mathrm{sg}}-\Delta {\omega }_g\right)\mathrm{dt}\end{array};$ where, J is the moment of inertia of the synchronous motor, Δω_sg is an angular frequency disturbance quantity of the synchronous motor, Δω_g is an angular frequency disturbance quantity of the voltage of the power grid, K_D is a droop coefficient, ω_N is a rated rotation speed, Δδ is a power angle disturbance quantity, and ΔP_sg is a power disturbance quantity.

The step of obtaining the reference value Δω of rotation speed variation of the compressor by the virtual inertia control segment includes the following steps.

The voltage U_g of the power grid is acquired and a virtual voltage orthogonal to the voltage of the power grid is generated by a differential method to form an αβ static orthogonal coordinate system. A Park transformation is performed with a phase angle θ_v of a virtual electromotive force vector of the virtual synchronous motor to obtain components of the voltage of the power grid relative to V_d and V_q. Then, a power angle δ_v of the virtual synchronous motor is:δ_v=arctan(V_q/V_d);

A virtual power P_v of the virtual synchronous motor is:

$P_v=\frac{U_{g}E}{X_a}\sin {\delta }_v\approx \frac{U_{g}E}{X_a}{\delta }_v.$

The reference value Δω of rotation speed variation of the compressor is:

$\Delta \omega =\frac{1}{n}\Delta P_v.$

Where, U_g is the voltage of the power grid, E is the virtual electromotive force, X_a is a virtual reactance, n is a constant, and ΔP_v is a power consumption variation.

In practical application, the voltage of a power distribution network generally fluctuates between 0.95 and 1.05 times of a rated voltage. Meanwhile, E is a fixed value since a power factor of the variable frequency air conditioner is fixedly and approximately equal to 1.0 under PFC control and cannot participate in adjustment of the voltage of the power grid. Therefore, the constant K_p=U_gE/X_a.

A single-phase AC power is an average value within one power frequency period, which has a certain lag for measuring the power. It may be considered that the output power of the variable frequency air conditioner approximately has a linear relationship with a rotation speed of the compressor. In order to make the power of the variable frequency air conditioner rapidly track a reference power of the virtual inertia stage, Δω may be set as:

$\Delta \omega =\frac{1}{n}\Delta P_v.$

The step of obtaining the reference value ω_ref of rotational angular frequency of the compressor by the droop control segment includes the following steps. A constant operating power P_N of the variable frequency air conditioner is acquired according to a set temperature of the variable frequency air conditioner and an indoor temperature. A steady-state output power P_D of the variable frequency air conditioner is calculated by a formula P_D=P_N+K_D(ω_N−ω_g), where, ω_g is an angular frequency of the voltage of the power grid, and K_D is a droop coefficient of the power grid. A difference between an actually measured power of the variable frequency air conditioner and the steady-state output power P_D of the variable frequency air conditioner is inputted into a PI controller to obtain the reference value ω_ref of rotational angular frequency of the compressor.

In summary, the Δω is added with the ω_ref obtained by the droop control segment, and then the result is inputted into the FOC for a frequency modulation for the voltage.

The step of obtaining the reference value ω_ref of rotational angular frequency of the compressor by the droop control segment further includes the following step(s).

A range of the power of the variable frequency air conditioner is assumed to be {0, P_max}, then an adjustable range of the power is P_N±min{(P_max−P_N), P_N} during frequency modulation of the power grid, where P_max is a maximum power of the variable frequency air conditioner. When a frequency error is 0.3 Hz, the power of the variable frequency air conditioner reaches an upper limit, and the droop coefficient K_D of the power grid is calculated by the following formula:

$K_D=\frac{0.8\min \left\{\left(P_{\max }-P_N\right),P_N\right\}}{0.3\times 2\pi };$

A small signal model of the virtual synchronous motor is established. An expression of the cutoff frequency ω_c of each variable frequency air conditioner is obtained by the following formula:

${\omega }_c=\sqrt{\sqrt{\frac{K_p^2}{J^2{\omega }_N^2}+\frac{K_D^4}{4J^4}}-\frac{K_D^2}{2J^2}},{\omega }_c\le \frac{K_p}{{\omega }_N{K}_D};$

When a certain phase angle margin (30° to 50°) is retained, the maximum value of the moment of inertia J is calculated according to the following expression of the PM:

$\mathrm{PM}=180°+\angle T_p\left({\omega }_c\right)=90°-\mathrm{arc}\tan \left(\frac{{\omega }_{c}J}{K_D}\right);$ where, ω_g is the angular frequency of the voltage of the power grid, P_max is the maximum power of the variable frequency air conditioner, K_p is a constant, and ∠T_p(ω_c) is a loop gain amplitude of the virtual inertia stage.

As shown in FIG. 3, the three-phase AC accesses a system for verification, in which the system includes three groups of variable frequency air conditioners and each group includes seven variable frequency air conditioners. Operating points, droop coefficients and other parameters for each group of variable frequency air conditioners are the same.

It is assumed that a set temperature of each variable frequency air conditioner and a room temperature of an ambient environment are known, and then a rated operating point of each variable frequency air conditioner is obtained, as shown in a second column of Table 1; the droop coefficient K_D is obtained according to the above formula for K_D, as shown in a third column of Table 1; and the moment of inertia of each variable frequency air conditioner is obtained by the PM, as shown in a fourth column of Table 1. Table 1 is shown as follows:

TABLE 1
Parameters of the variable frequency air conditioners in one group
	Rated operating	Droop coefficient	Moment of
No.	point (W)	(W/(rad/s))	inertia (kg · m2)
1	500	212	0.122
2	750	318	0.108
3	1000	425	0.092
4	1250	530	0.076
5	1500	425	0.092
6	1750	318	0.108
7	2000	212	0.122

In practical application, the variable frequency air conditioners can be controlled by a normal droop control or the method of the present application when the load fluctuates. FIG. 3 and FIG. 4 show simulation results of dynamic processes of the system of the variable frequency air conditioners under the normal droop control and the proposed control strategy respectively when the load fluctuates. It should be illustrated that the proposed algorithms illustrated in FIG. 3 and FIG. 4 are the load control of the variable frequency air conditioner based on the virtual synchronous motor. It can be seen from FIG. 3 that at the moment of impact loads accessing a speed of the power reducing of the variable frequency air conditioner is faster than that under the control strategy proposed in droop control, and a power imbalance of the microgrid is reduced more rapidly. It can be seen from FIG. 4 that a frequency drop rate of the microgrid is greatly slowed down after applying the proposed control strategy.

As shown in FIG. 5, which is a schematic diagram illustrating a structure of the FOC provided by the present disclosure A working process of the FOC is as follows.

In step 1, two-phase currents i_a and i_b are collected at first.

In step 2, a clarke transformation is performed to obtain two-axis orthogonal current magnitudes i_α and i_β.

In step 3, a rotation transformation is performed to obtain orthogonal current magnitudes i_d and i_q, where i_q is related to a torque and i_d is related to a magnetic flux. In an actual control, i_d is often set as zero. the two obtained magnitudes are not time-varying, so the two magnitudes can be controlled separately, which is similar to a direct current (DC) quantity control. However, It is unnecessary to know the specific voltages applied to the three phases of the motor.

In step 4, the magnitudes i_d and i_q obtained in the step 3 are respectively sent into a PI adjuster to obtain the corresponding outputs v_d and v_q.

In step 5, a rotation angle of the motor is obtained by a sensor.

In step 6, an inverse park transformation is performed to obtain two-axis current magnitudes v_α and v_β.

In step 7, an inverse clarke transformation is performed on v_α and v_β in the step 6 to obtain a three-phase voltage required actually for being inputted to an inverter bridge to drive the motor to rotate.

According to a technical solution of the present embodiment, when the frequency of the power grid fluctuates, the output power of the variable frequency air conditioner can simulate a change process of the output power and the rotation speed of a by establishing the mathematical simulation model of the virtual synchronous motor with the same capacity as the no-load synchronous motor. The solution can be adopted to solve the problem of insufficient moment of inertia of the power grid caused after a new energy resource accesses to the power grid. Thus, the moment of inertia of the microgrid can be improved, and the fluctuation range of AC frequency can be reduced when the load fluctuates.

Embodiments of the present disclosure further provide a power grid adjustment method based on a load of a variable frequency air conditioner. FIG. 6 is a schematic diagram illustrating another power grid adjustment method based on the load of the variable frequency air conditioner provided by an embodiment of the present disclosure. As shown in FIG. 6, the method includes the following steps.

In S10, a mathematical simulation model of a virtual synchronous motor of the variable frequency air conditioner is established.

In S20, an adjustment structure of the virtual synchronous motor is established according to the mathematical simulation model of the virtual synchronous motor.

With reference to FIG. 2, the adjustment structure may be a virtual inertia control segment and a droop control segment. An inertia adjustment proportion parameter of the virtual inertia control segment and a droop adjustment proportion parameter of the droop control segment can be determined according to the mathematical simulation model of the virtual

$\begin{array}{cc}J\frac{d\Delta {\omega }_{\mathrm{sg}}}{\mathrm{dt}}=\frac{-K_D\Delta {\omega }_{\mathrm{sg}}-\Delta P_{\mathrm{sg}}}{{\omega }_N}& \left(1\right)\end{array}$

synchronous motor. For example, in of the mathematical simulation model of the virtual synchronous motor, ΔP_sg is a power disturbance quantity of the power grid and measured by the power grid in real time. When the frequency of the power grid changes, the power disturbance quantity ΔP_sg of the power grid is measured. It can be known from formula (1) that the power disturbance quantity ΔP_sg of the power grid includes two parts: one part is a product of K_DΔω_sg which is expressed as the droop control segment of the adjustment structure, and the other part is

$J\frac{d\Delta {\omega }_{\mathrm{sg}}}{\mathrm{dt}}$ which is expressed as the virtual inertia control segment. Therefore, the adjustment structure of the virtual synchronous motor can be established according to the mathematical simulation model of the virtual synchronous motor. In general, the formula (1) is represented in the adjustment structure by determining the parameters in the adjustment structure. For example, a first-order system parameter in the virtual inertia control segment is determined by J in the mathematical simulation model of the virtual synchronous motor, and a droop coefficient in the droop control segment is determined by K_D of the mathematical simulation model of the virtual synchronous motor.

In S30, the parameters of the power grid are acquired in real time.

With reference to FIG. 2, the adjustment structure acquires a real-time frequency, a voltage, a current and the like of the power grid from the power grid before calculation, thereby determining frequency change of the power grid and the power acquired by the variable frequency air conditioner from the power grid in real time.

In S40, the parameters of the power grid are inputted into the adjustment structure for calculation to obtain a reference value Δω of rotation speed variation of the variable frequency air conditioner and a reference value ω_ref of rotational angular frequency which is needed to be changed.

With reference to FIG. 2, the parameters of the power grid may be inputted into the adjustment structure for calculation to obtain the reference value Δω of rotation speed variation of the variable frequency air conditioner and the reference value ω_ref of rotational angular frequency which is needed to be changed after determining the parameters in the virtual inertia control segment and the droop control segment of the adjustment structure and acquiring the parameters of the power grid in real time. In an exemplary embodiment, the voltage of the power grid is inputted into the virtual inertia control segment for calculation to obtain the reference value Δω of rotation speed variation of the variable frequency air conditioner, and the frequency of the power grid is inputted into the droop control segment for calculation to obtain the reference value ω_ref of rotational angular frequency which is needed to be changed.

In S50, a working point of the variable frequency air conditioner is adjusted according to the reference value Δω of rotation speed variation of the variable frequency air conditioner and the reference value ω_ref of rotational angular frequency.

With reference to FIG. 2, an angular frequency value needed to be adjusted is obtained by a difference value between the reference value Δω of rotation speed variation of the variable frequency air conditioner and the reference value ω_ref of rotational angular frequency, and the angular frequency value is inputted into the FOC for controlling a PWM circuit, so as to adjust the working point of the variable frequency air conditioner. Thus, the working point of the variable frequency air conditioner is adaptively adjusted according to changes of the parameters of the power grid. When the frequency of the power grid is reduced, the variable frequency air conditioner adaptively reduces the working frequency to reduce the power acquired from the power grid, further to reduce the moment of inertia acquired by the variable frequency air conditioner from the power grid and reduce the fluctuation range of the frequency of the power grid when the power grid fluctuates.

Before the mathematical simulation model of the virtual synchronous motor of the variable frequency air conditioner is established, the method further includes the following steps. A signal model of the virtual synchronous motor is established, and a cutoff frequency ω_c of each variable frequency air conditioner is obtained; a PM of the virtual synchronous motor is acquired, and a moment of inertia J of the virtual synchronous motor is calculated according to the PM and the cutoff frequency ω_c of the variable frequency air conditioner; and a maximum value of the moment of inertia J of the virtual synchronous motor is obtained according to a value range of the PM.

The mathematical simulation model of the virtual synchronous motor can be conveniently calculated by determining the value range of the moment of inertia J of the virtual synchronous motor before establishing the mathematical simulation model of the virtual synchronous motor of the variable frequency air conditioner.

The mathematical simulation model is:

$\hspace{1em}\{\begin{array}{c}J\frac{d\Delta {\omega }_{\mathrm{sg}}}{\mathrm{dt}}=\frac{-K_D\Delta {\omega }_{\mathrm{sg}}-\Delta P_{\mathrm{sg}}}{{\omega }_N}\\ \Delta \delta =\int \left(\Delta {\omega }_{\mathrm{sg}}-\Delta {\omega }_g\right)\mathrm{dt}\end{array},$ Δω_sg is an angular frequency disturbance quantity of the virtual synchronous motor, Δω_g is an angular frequency disturbance quantity of a voltage of the power grid, K_D is a droop coefficient, ω_N is a rated rotation speed, Δδ is a power angle disturbance quantity, ΔP_sg is a power disturbance quantity, and J is the moment of inertia of the virtual synchronous motor.

The step of inputting the voltage of the power grid into the virtual inertia control segment for calculation to obtain the reference value Δω of rotation speed variation of the variable frequency air conditioner includes the following steps.

The voltage U_g of the power grid is acquired and a virtual voltage orthogonal to the voltage of the power grid is generated by a differential method to form an αβ static orthogonal coordinate system. A Park transformation is performed with a phase angle θ_v of a virtual electromotive force vector of the virtual synchronous motor to obtain components of the voltage of the power grid relative to V_d and V_q. Then, a power angle δ_v of the virtual synchronous motor is:δ_v=arctan(V_q/V_d).

A virtual power P_v of the virtual synchronous motor is:

$P_v=\frac{U_{g}E}{X_a}\sin {\delta }_v\approx \frac{U_{g}E}{X_a}{\delta }_v$

The reference value Δω of rotation speed variation of the compressor is:

$\Delta \omega =\frac{1}{n}\Delta P_v.$

Where, U_g is the voltage of the power grid, E is the virtual electromotive force, X_a is a virtual reactance, n is a constant, and ΔP_v is a power consumption variation.

In practical application, the voltage of a power distribution network generally fluctuates between 0.95 and 1.05 times of a rated voltage. Meanwhile, E is a fixed value since a power factor of the variable frequency air conditioner is fixedly and approximately equal to 1.0 under PFC control and cannot participate in adjustment of the voltage of the power grid. Therefore, the constant K_p=U_gE/X_a.

A single-phase AC power is an average value within one power frequency period, which has a certain lag for measuring the power. It may be considered that the output power of the variable frequency air conditioner approximately has a linear relationship with a rotation speed of the compressor. In order to make the power of the variable frequency air conditioner rapidly track a reference power of the virtual inertia stage, Δω may be set as:

$\Delta \omega =\frac{1}{n}\Delta P_v.$

The step of inputting the frequency of the power grid into the droop control segment for calculation to obtain the reference value ω_ref of rotational angular frequency to be changed includes the following steps. A constant operating power P_N of the variable frequency air conditioner is acquired according to a set temperature of the variable frequency air conditioner and an indoor temperature. A steady-state output power P_D of the variable frequency air conditioner is calculated by a formula P_D=P_N±K_D(ω_N−ω_g), where, ω_g is an angular frequency of the voltage of the power grid, and K_D is a droop coefficient of the power grid. A difference value between an actually measured power of the variable frequency air conditioner and the steady-state output power P_D of the variable frequency air conditioner is inputted into a PI controller to obtain the reference value ω_ref of rotational angular frequency of the compressor.

Specifically, before the frequency of the power grid is inputted into the droop control stage for calculation to obtain the reference value ω_ref of rotational angular frequency which is needed to be changed, the method further includes: obtaining the droop coefficient of the droop control stage. The droop coefficient K_D is:

$K_D=\frac{0.8\min \left\{\left(P_{\max }-P_N\right),P_N\right\}}{0.3\times 2\pi };$ where, P_max is a maximum power of the variable frequency air conditioner, P_N is a rated power of the variable frequency air conditioner, and an adjustable range of the power of the variable frequency air conditioner is P_N±min{(P_max−P_N), P_N}.

The cutoff frequency ω_c of the variable frequency air conditioner is:

${\omega }_c=\sqrt{\sqrt{\frac{K_p^2}{J^2{\omega }_N^2}+\frac{K_D^4}{4J^4}}-\frac{K_D^2}{2J^2}},{\omega }_c\le \frac{K_p}{{\omega }_N{K}_D};$ where K_p is a constant.

A relationship between a value range of the PM and the maximum value of the moment of inertia J of the virtual synchronous motor is:

$\mathrm{PM}=180°+\angle T_p\left({\omega }_c\right)=90°-\arctan \left(\frac{{\omega }_{c}J}{K_D}\right);$

where, PM is the phase margin, and ∠T_p(ω_c) is a loop gain amplitude of the virtual inertia stage.

As shown in FIG. 3, the three-phase AC accesses a system for verification, in which the system includes three groups of variable frequency air conditioners and each group includes seven variable frequency air conditioners. Operating points, droop coefficients and other parameters for each group of variable frequency air conditioners are the same.

It is assumed that a set temperature of each variable frequency air conditioner and a room temperature of an ambient environment are known, and then a rated operating point of each variable frequency air conditioner is obtained, with continued reference to Table 1, as shown in a second column of Table 1; the droop coefficient K_D is obtained according to the above formula for K_D, as shown in a third column of Table 1; and the moment of inertia of each variable frequency air conditioner is calculated and obtained by the PM, as shown in a fourth column of Table 1.

In practical application, the variable frequency air conditioners can be controlled by a normal droop control or the method of the present application when the load fluctuates. FIG. 3 and FIG. 4 show simulation results of dynamic processes of the system of the variable frequency air conditioners under the normal droop control and the proposed control strategy respectively when the load fluctuates. It should be illustrated that, FIG. 3 is a schematic diagram illustrating the simulation results of a total power of a single group of variable frequency air conditioners under different control methods provided by embodiments of the present invention, and FIG. 4 is a schematic diagram illustrating the simulation results of the frequency of the power grid under different control methods provided by embodiments of the present invention. The proposed algorithm illustrated in FIG. 3 and FIG. 4 are an algorithm for the load control of the variable frequency air conditioner based on the virtual synchronous motor. It can be seen from FIG. 3 that at the moment of impact loads accessing a speed of the power reducing of the variable frequency air conditioner is faster than that under the control strategy proposed in droop control, and a power imbalance of the power grid is reduced more rapidly. It can be seen from FIG. 4 that a frequency drop rate of the power grid is greatly slowed down after applying the proposed control strategy.

As shown in FIG. 5, which is a schematic diagram illustrating a structure of the FOC provided by an embodiment of the present invention, a working process of the FOC is as follows. The working point of the variable frequency air conditioner is adjusted according to the reference value Δω of rotation speed variation of the variable frequency air conditioner and the reference value ω_ref of rotational angular frequency, and this step includes the following steps.

Two-phase currents i_a and i_b of the variable frequency air conditioner are collected.

The clarke transformation is performed on the two-phase currents i_a and i_b to obtain two-axis orthogonal current magnitudes i_α and i_β.

The rotation transformation is performed on the orthogonal current magnitudes i_α and i_β to obtain orthogonal current magnitudes i_d and i_q, where i_q is related to a torque and i_d is related to a magnetic flux.

The obtained i_d and i_q are respectively inputted into the PI controller to obtain corresponding outputs v_d and v_q.

A rotation angle θ of the motor of the variable frequency air conditioner is obtained by a sensor, and the inverse park transformation is performed with v_d and v_q to obtain two-axis current magnitudes v_α and v_β.

The inverse clarke transformation is performed on v_α and v_β to obtain the required three-phase voltage for being inputted to an inverter bridge to drive the variable frequency air conditioner to work.

According to a technical solution of the present embodiment, when the frequency of the power grid fluctuates, the output power of the variable frequency air conditioner can simulate a change process of the output power and the rotation speed of a no-load synchronous motor by establishing the mathematical simulation model of the virtual synchronous motor with the same capacity as the no-load synchronous motor. The solution can be adopted to solve the problem of insufficient moment of inertia of the power grid caused after a new energy resource accesses to the power grid. Thus, the moment of inertia of the microgrid can be improved, and the fluctuation range of AC frequency can be reduced when the load fluctuates.

It should be noted that, the above only describes preferred embodiments of the present disclosure and applied technical principles. Those skilled in the art shall understand that the present disclosure is not limited to specific embodiments described herein. For those skilled in the art, various apparent variations, readjustments and replacements can be made without departing from a protection scope of the present disclosure. Therefore, although the present disclosure is described in detail through the above embodiments, the present disclosure is not limited to the above embodiments and may further include more other equivalent embodiments without departing from the concept of the present disclosure, while the scope of the present disclosure is decided by a scope of attached claims.

Claims

1. A power grid adjustment method based on a load of a variable frequency air conditioner, comprising: establishing a mathematical simulation model of a virtual synchronous motor in the variable frequency air conditioner controller;establishing a virtual inertia control segment and a droop control segment of power grid adjustment according to the mathematical simulation model;obtaining a reference value Δω of a rotation speed variation of a compressor by the virtual inertia control segment;obtaining a reference value ωref of a rotational angular frequency of the compressor by the droop control segment; andinputting a sum of the reference value ωref of the rotational angular frequency of the compressor and the reference value Δω of the rotation speed variation of the compressor into a field-oriented controller to control a rotation of the motor.

2. The power grid adjustment method according to claim 1, further comprising: establishing a signal model of the virtual synchronous motor, and obtaining a cutoff frequency ωc of each variable frequency air conditioner;acquiring a phase margin (PM) of the virtual synchronous motor, and calculating a moment of inertia J of the virtual synchronous motor according to the PM and the cutoff frequency ωc of the variable frequency air conditioner; andobtaining a maximum value of the moment of inertia J of the virtual synchronous motor according to a value range of the PM.

3. The power grid adjustment method according to claim 2, wherein the mathematical simulation model is

4. The power grid adjustment method according to claim 3, wherein the step of obtaining a reference value Δω of a rotation speed variation of a compressor by the virtual inertia control segment comprises: acquiring the voltage Ug of the power grid and generating a virtual voltage orthogonal to the voltage of the power grid by a differential method to form an αβ static orthogonal coordinate system, and performing a Park transformation with a phase angle θ of a virtual electromotive force vector of the virtual synchronous motor to obtain components of the voltage of the power grid relative to Vd and Vq; and then,a power angle δv of the virtual synchronous motor is:δv=arctan(Vq/Vd);a virtual power Pv of the virtual synchronous motor is:

5. The power grid adjustment method according to claim 4, wherein the step of obtaining a reference value ωref of a rotational angular frequency of the compressor by the droop control segment comprises: acquiring a constant operating power PN of the variable frequency air conditioner according to a set temperature of the variable frequency air conditioner and an indoor temperature, and calculating by a formula PD=PN+KD(ωN−ωg) to obtain a steady-state output power PD of the variable frequency air conditioner, wherein ωg is an angular frequency of the voltage of the power grid, and KD is a droop coefficient of the power grid; andinputting a difference between an actually measured power of the variable frequency air conditioner and the steady-state output power PD of the variable frequency air conditioner into a PI controller to obtain the reference value ωref of the rotational angular frequency of the compressor.

6. The power grid adjustment method according to claim 5, wherein the step of obtaining a reference value ωref of a rotational angular frequency of the compressor by the droop control segment further comprises: setting a range of the power of the variable frequency air conditioner as {0, Pmax}, then an adjustable range of the power is PN±min{(Pmax-PN), PN} during a frequency modulation of the power grid, where, Pmax is a maximum power of the variable frequency air conditioner, and the power of the variable frequency air conditioner reaches an upper limit when a frequency error is 0.3 Hz; andcalculating the droop coefficient KD of the power grid by a formula:

7. The power grid adjustment method according to claim 6, further comprising:

8. The power grid adjustment method according to claim 7, wherein the step of inputting a sum of the reference value ωref of the rotational angular frequency of the compressor and the reference value Δω of the rotation speed variation of the compressor into a field-oriented controller to control a rotation of the motor comprises: collecting two-phase currents ia and ib;performing clarke transformation on the two-phase currents ia and ib to obtain magnitudes of two-axis orthogonal currents iα and iβ;performing rotation transformation on the magnitudes of the orthogonal currents iα and iβ to obtain orthogonal current magnitudes id and iq, wherein iq is related to a torque, and id is related to a magnetic flux;inputting the obtained id and iq respectively into the PI controller to obtain corresponding outputs vd and vq;obtaining a rotation angle θ of the motor by a sensor, and performing an inverse park transformation with vd and vq to obtain magnitudes of two-axis currents vα and vβ; andperforming an inverse clarke transformation on vα and vβ to obtain a required three-phase voltage, and inputting the three-phase voltage to an inverter bridge to drive the motor to rotate.

9. A power grid adjustment method based on a load of a variable frequency air conditioner, comprising: establishing a mathematical simulation model of a virtual synchronous motor of the variable frequency air conditioner;establishing an adjustment structure of the virtual synchronous motor according to the mathematical simulation model of the virtual synchronous motor;acquiring parameters of a power grid in real time;inputting the parameters of the power grid into the adjustment structure for calculation to obtain a reference value Δω of a rotation speed variation of the variable frequency air conditioner and a reference value ωref of a rotational angular frequency of the variable frequency air conditioner; andadjusting a working point of the variable frequency air conditioner according to the reference value Δω of the rotation speed variation of the variable frequency air conditioner and the reference value ωref of the rotational angular frequency.

10. The method according to claim 9, wherein before the step of establishing a mathematical simulation model of a virtual synchronous motor of a variable frequency air conditioner, the method further comprises: establishing a signal model of the virtual synchronous motor, and obtaining a cutoff frequency ωc of each variable frequency air conditioner;acquiring a phase margin (PM) of the virtual synchronous motor, and calculating a moment of inertia J of the virtual synchronous motor according to the PM and the cutoff frequency ωc of the variable frequency air conditioner; andobtaining a maximum value of the moment of inertia J of the virtual synchronous motor according to a value range of the PM.

11. The method according to claim 10, wherein the mathematical simulation model is:

12. The method according to claim 11, wherein the adjustment structure comprises a virtual inertia control segment and a droop control segment; and the step of establishing an adjustment structure of the virtual synchronous motor according to the mathematical simulation model of the virtual synchronous motor comprises: determining an inertia adjustment proportion parameter of the virtual inertia control segment and a droop adjustment proportion parameter of the droop control segment according to the mathematical simulation model of the virtual synchronous motor.

13. The method according to claim 12, wherein the step of inputting the parameters of the power grid into the adjustment structure for calculation to obtain a reference value Δω of a rotation speed variation of the variable frequency air conditioner and a reference value ωref of a rotational angular frequency of the variable frequency air conditioner comprises: inputting a voltage of the power grid into the virtual inertia control segment for calculation to obtain the reference value Δω of the rotation speed variation of the variable frequency air conditioner; andinputting a frequency of the power grid into the droop control segment for calculation to obtain the reference value ωref of the rotational angular frequency of the variable frequency air conditioner.

14. The method according to claim 13, wherein the step of inputting the voltage of the power grid into the virtual inertia control segment for calculation to obtain the reference value Δω of the rotation speed variation of the variable frequency air conditioner comprises: acquiring the voltage Ug of the power grid and generating a virtual voltage orthogonal to the voltage of the power grid by a differential method to form an αβ static orthogonal coordinate system, and performing a Park transformation with a phase angle θv of a virtual electromotive force vector of the virtual synchronous motor to obtain components of the voltage of the power grid relative to Vd and Vq; and then,a power angle δv of the virtual synchronous motor is: δv=arctan(Vq/Vd);a virtual power Pv of the virtual synchronous motor is:

15. The method according to claim 13, wherein the step of inputting the frequency of the power grid into the droop control segment for calculation to obtain the reference value ωref of the rotational angular frequency of the variable frequency air conditioner which is needed to be changed comprises: acquiring a constant operating power PN of the variable frequency air conditioner according to a set temperature of the variable frequency air conditioner and an indoor temperature, and calculating by a formula PD=PN+KD(ωN−ωg) to obtain a steady-state output power PD of the variable frequency air conditioner, wherein ωg is an angular frequency of the voltage of the power grid, and KD is a droop coefficient of the power grid; andinputting a difference between an actually measured power of the variable frequency air conditioner and the steady-state output power PD of the variable frequency air conditioner into a PI controller to obtain the reference value ωref of the rotational angular frequency of a compressor.

16. The method according to claim 15, wherein before inputting the frequency of the power grid into the droop control segment for calculation to obtain the reference value ωref of the rotational angular frequency of the variable frequency air conditioner, the method further comprises: acquiring the droop coefficient of the droop control segment; the droop coefficient KD is:

17. The method according to claim 13, wherein the cutoff frequency ωc of the variable frequency air conditioner is:

18. The method according to claim 17, wherein a relationship between a value range of the PM and the maximum value of the moment of inertia J of the virtual synchronous motor is:

19. The method according to claim 9, wherein the step of adjusting a working point of the variable frequency air conditioner according to the reference value Δω of the rotation speed variation of the variable frequency air conditioner and the reference value ωref of the rotational angular frequency comprises: collecting two-phase currents ia and ib of the variable frequency air conditioner;performing clarke transformation on the two-phase currents ia and ib to obtain magnitudes of two-axis orthogonal currents iα and iβ;performing rotation transformation on the magnitudes of the orthogonal currents iα and iβ to obtain orthogonal current magnitudes id and iq, where, iq is related to a torque, and id is related to a magnetic flux;inputting the obtained id and iq respectively into the PI controller to obtain corresponding outputs vd and vq;obtaining a rotation angle θ of the motor of the variable frequency air conditioner by a sensor, and performing an inverse park transformation with vd and vq to obtain magnitudes of two-axis currents vα and vβ; andperforming an inverse clarke transformation on vα and vβ to obtain a required three-phase voltage, and inputting the three-phase voltage to an inverter bridge to drive the variable frequency air conditioner to work.