CN116885742A - Primary frequency modulation parameter adjustment method, system and medium based on model simulation - Google Patents

Abstract

The invention discloses a primary frequency modulation parameter adjustment method, a system and a medium based on model simulation, wherein the method comprises the steps of determining a primary frequency modulation parameter to be adjusted in a pre-established control system model, determining an empirical value of the primary frequency modulation parameter as a current setting value of the primary frequency modulation parameter to be assigned to the control system model, carrying out primary frequency modulation simulation by combining a pre-established main steam pressure model and a steam turbine model to obtain total primary frequency modulation action times and qualified primary frequency modulation times, calculating the qualified primary frequency modulation rate according to the total primary frequency modulation action times and the qualified primary frequency modulation times, judging whether the requirement is met, and if the requirement is not met, correcting the current setting value of the primary frequency modulation parameter to be adjusted, and continuing iterative simulation; otherwise, ending and exiting. The invention aims to enable the primary frequency modulation parameters of the thermal power unit to better overcome the influence of spontaneous fluctuation of the load and the main steam pressure of the thermal power unit, and finally enables the thermal power unit to have better primary frequency modulation power response performance.

Description

Primary frequency modulation parameter adjustment method, system and medium based on model simulation

Technical Field

The invention belongs to the technical field of primary frequency modulation of thermal power generating units, and particularly relates to a primary frequency modulation parameter adjustment method, a primary frequency modulation parameter adjustment system and a primary frequency modulation parameter adjustment medium based on model simulation.

Background

The primary frequency modulation performance of the thermal power generating unit is influenced by a primary frequency modulation control strategy and primary frequency modulation parameters. After the primary frequency modulation control strategy is determined, how to determine the primary frequency modulation parameters is an important factor affecting primary frequency modulation performance. In the conventional method, an empirical method is a main means for adjusting primary frequency modulation parameters. For working conditions such as step disturbance test, power grid accident frequency simulation test and the like, the empirical method can generally meet the requirement of primary frequency modulation parameter adjustment. However, when the power grid normally operates, the fluctuation amplitude of the grid frequency is small, and spontaneous fluctuation of the load and the main steam pressure of the thermal power unit can obviously influence the power response characteristic of primary frequency modulation. The empirical method is difficult to fully meet the requirements of primary frequency modulation parameter adjustment and performance optimization under the working condition. Therefore, how to fully consider the influence of spontaneous fluctuation of the load and the main steam pressure of the thermal power unit in the primary frequency modulation parameter adjustment process, and guaranteeing that the thermal power unit meets the requirement of primary frequency modulation performance assessment in the normal operation mode of the power grid is still a main difficulty of the conventional primary frequency modulation parameter adjustment.

Disclosure of Invention

The invention aims to solve the technical problems: aiming at the problems in the prior art, the invention provides a primary frequency modulation parameter adjustment method, a primary frequency modulation parameter adjustment system and a primary frequency modulation parameter adjustment medium based on model simulation, which aim to enable a primary frequency modulation parameter of a thermal power unit to better overcome the influence of spontaneous fluctuation of load and main steam pressure of the thermal power unit and finally enable the thermal power unit to have better primary frequency modulation power response performance.

In order to solve the technical problems, the invention adopts the following technical scheme:

a primary frequency modulation parameter adjustment method based on model simulation comprises the following steps:

s1, determining a primary frequency modulation parameter to be adjusted in a pre-established control system model, and determining an empirical value of the primary frequency modulation parameter as a current setting value of the primary frequency modulation parameter;

s2, assigning the current set value of the primary frequency modulation parameter to be adjusted to a control system model, carrying out primary frequency modulation simulation by combining a pre-established main steam pressure model and a steam turbine model through the control system model to obtain the total number of primary frequency modulation actions and the qualified number of primary frequency modulation, and calculating the qualified rate of primary frequency modulation according to the total number of primary frequency modulation actions and the qualified number of primary frequency modulation;

s3, judging whether the primary frequency modulation qualification rate meets the requirement, if not, correcting the current set value of the primary frequency modulation parameter to be adjusted, and jumping to the step S2; otherwise, ending and exiting.

Optionally, the control system model of primary frequency modulation of the thermal power generating unit, which is pre-established in step S1, includes:

the rotating speed unequal rate function module is used for forming a frequency modulation instruction by the input power grid frequency through a rotating speed unequal rate function;

the controller set value generation module is used for adding the frequency modulation instruction and the load setting to form a controller set value;

the load deviation calculation module is used for feeding back the set value of the controller and the input load to calculate the load deviation;

the closed loop controller is used for generating a controller output according to the load deviation;

the gain link is used for taking the frequency modulation instruction as frequency modulation feedforward after gain;

and the summation module is used for adding the controller output and the frequency modulation feedforward or directly outputting the controller output and the frequency modulation feedforward as the control system model output.

Optionally, the closed loop controller is a PI controller.

Optionally, the primary frequency modulation parameters to be adjusted in the control system model determined in step S1 include: gain factor K, PI controller of gain link _p And integral coefficient K _I 。

Optionally, the functional expression of the pre-established main steam pressure model in step S2 is:

in the above-mentioned method, the step of,p _T is the main steam pressure;Ais the amplitude;fis frequency;ttime is; θ ₀ Is the phase angle;bis biased.

Optionally, the turbine model pre-established in step S2 includes:

a valve flow characteristic link for outputting a control system model of the control system model by using a preset valve flow characteristic functionf(x) Processing to obtain a first intermediate value;

the multiplication module is used for multiplying the main steam pressure output by the main steam pressure model by a first intermediate value;

a regulating stage pressure-load link for applying a preset regulating stage pressure-load function to the output of the multiplication moduleg(x) Processing to obtain a second intermediate value;

a transfer function module for adopting a preset transfer function for the second intermediate valueG(s) Generating a third intermediate value;

a summing module for adding a third intermediate value, the input steady state power response deltaAnd generating active power P output by the steam turbine model.

Optionally, the preset valve flow characteristic functionf(x) The abscissa is the output of the control system model and the ordinate is the pressure ratio for the first intermediate value.

Optionally, the preset regulation stage pressure-load functiong(x) The abscissa of (2) is the regulation stage pressure and the ordinate is the load.

Optionally, the preset regulation stage pressure-load functionG(s) The functional expression of (2) is:

，

In the above-mentioned method, the step of,F _HP is the power ratio of the high-pressure cylinder,λis the power overshoot factor of the high-pressure cylinder,T _CH is the volume time constant of the high-pressure steam chamber,T _RH is the reheater volume time constant.

Optionally, the steady state power response delta of the inputThe functional expression of (2) is:

in the above-mentioned method, the step of,mfor the number of phase angles,is a steady-state power amplitude characteristic parameter +>I element of->Is a steady-state power frequency characteristic parameter->Is selected from the group consisting of the (i) th element,ttime is; />Is the firstiPhase angle, & gt>，randIs a random one between 0 and 1A machine number; wherein the steady-state power amplitude characteristic parameter->And steady state power frequency characteristic parameterF _f The acquisition of (1) comprises:

s101, extracting steady-state load data of the unit within a certain period of time，/>，/>Is the firstiSteady state load data at a point in time,nis steady-state load data length, unit steady-state load data +.>Time interval of adjacent elements->Equal, corresponding sampling rate is +.>；

S102, calculating steady-state load variables of the unit according to the following formula：

In the above-mentioned method, the step of,for steady state load data of the unit->Is the average value of (2); />；

S103, obtaining an intermediate variable through discrete Fourier calculation，/>Wherein, and have:

，

in the above-mentioned method, the step of,for intermediate variable +.>I element of->Is->The first of (3)kThe number of elements to be added to the composition,jin units of imaginary numbers,；

s104, calculating steady-state power amplitude characteristic parameters according to the following ，/>Elements of (a) and (b):

，/>，/>，

in the above-mentioned method, the step of,is a steady-state power amplitude characteristic parameter +>Element 1 of->Is a steady-state power amplitude characteristic parameter +>To the first part of (1)l2) +1 elements, ++>Is a steady-state power amplitude characteristic parameter +>At the time of the mth element of (1)mThe value range of (2) is +.>，/>For intermediate variable +.>Element 1 of->For intermediate variable +.>To the first part of (1)l2) +1 elements; calculating steady-state power frequency characteristic parameters according to the following formula>，/>：

，

In the above-mentioned method, the step of,is a steady-state power frequency characteristic parameter->M element of->Is the sampling rate.

Optionally, in step S2, when performing primary frequency modulation simulation by combining the control system model with the pre-established main steam pressure model and the steam turbine model to obtain the total number of primary frequency modulation actions and the qualified number of primary frequency modulation actions, the method includes:

s201, acquiring a power grid frequency fluctuation data set，/>Wherein the firstiIndividual element->，Wherein->The first of (3)jIndividual element->Represents->Grid frequency at time,/->Is the grid frequency sampling interval; initializing the total action times of primary frequency modulation and the qualified times of primary frequency modulation;

s202, randomly selecting a power grid frequency fluctuation data setAny elementF _i As the input of the control system model, the active power P output by the steam turbine model is generated by carrying out single primary frequency modulation simulation by combining the control system model with a pre-established main steam pressure model and the steam turbine model;

S203, by elementF _i And calculating the actual integral contribution electric quantity of the simulation according to the active power P output by the steam turbine modelH _i And theoretical integral contribution powerH _e And will actually integrate the contributing powerH _i Divided by theoretical integral contributionH _e Obtaining the primary frequency modulation contribution power rate of the simulationkThe method comprises the steps of carrying out a first treatment on the surface of the Adding 1 to the total action frequency of the primary frequency modulation, if the primary frequency modulation contributes to the electric quantity ratekIf the frequency is larger than or equal to the set value, adding 1 to the qualified frequency of primary frequency modulation;

and S204, judging whether the total frequency modulation action times are greater than or equal to a preset value, if so, judging that the total frequency modulation action times and the qualified frequency modulation times are obtained by carrying out primary frequency modulation simulation through a control system model in combination with a pre-established main steam pressure model and a steam turbine model, otherwise, continuing to iterate in the step S202.

Optionally, the actual integrated contribution powerH _i The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ the system frequency exceeds the dead zone of primary frequency modulation action of the unit,t _t the moment when the system frequency enters the dead zone of the primary frequency modulation action of the unit,P _t is thattThe time unit actually generates active power,P ₀ is thatt ₀ And the time unit actually generates active power, and t is time.

Optionally, the theoretical integration contribution power H _e The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ for the moment when the system frequency exceeds the dead zone of the primary frequency modulation action of the unit,t _t the time when the system frequency enters the primary frequency modulation action dead zone of the unit,is thattAt the moment when the power grid frequency exceeds the frequency modulation artificial dead zone value, < ->Rated active output for the unit, < >>Rated frequency for the power grid, < >>The speed change rate of the unit is set, and t is time.

Optionally, in step S3, calculating the qualified rate of primary frequency modulation according to the total number of primary frequency modulation actions and the qualified rate of primary frequency modulation refers to dividing the qualified rate of primary frequency modulation by the total number of primary frequency modulation actions to obtain the qualified rate of primary frequency modulation.

In addition, the invention also provides a primary frequency modulation parameter adjustment system based on model simulation, which comprises the following steps:

the primary frequency modulation parameter selecting and value taking program unit is used for determining primary frequency modulation parameters to be adjusted in a pre-established control system model and determining the experience value of the primary frequency modulation parameters as the current setting value of the primary frequency modulation parameters;

the primary frequency modulation simulation program unit is used for assigning the current set value of the primary frequency modulation parameter to be adjusted to the control system model, carrying out primary frequency modulation simulation by combining the control system model with a pre-established main steam pressure model and a steam turbine model to obtain the total number of primary frequency modulation actions and the qualified number of primary frequency modulation, and calculating the qualified rate of primary frequency modulation according to the total number of primary frequency modulation actions and the qualified number of primary frequency modulation;

The simulation result judging program unit is used for judging whether the primary frequency modulation qualification rate meets the requirement, if not, correcting the current set value of the primary frequency modulation parameter to be adjusted, and calling the primary frequency modulation simulation program unit to continuously execute a round of primary frequency modulation simulation; otherwise, ending and exiting.

Optionally, the primary frequency modulation parameter selection and value-taking program unit pre-establishes a control system model of primary frequency modulation of the thermal power generating unit, which comprises the following steps:

the rotating speed unequal rate function module is used for forming a frequency modulation instruction by the input power grid frequency through a rotating speed unequal rate function;

the controller set value generation module is used for adding the frequency modulation instruction and the load setting to form a controller set value;

the load deviation calculation module is used for feeding back the set value of the controller and the input load to calculate the load deviation;

the closed loop controller is used for generating a controller output according to the load deviation;

the gain link is used for taking the frequency modulation instruction as frequency modulation feedforward after gain;

and the summation module is used for adding the controller output and the frequency modulation feedforward or directly outputting the controller output and the frequency modulation feedforward as the control system model output.

Optionally, the closed loop controller is a PI controller.

Optionally, the primary frequency modulation parameters to be adjusted in the control system model determined by the primary frequency modulation parameter selection and value-taking program unit include: gain factor K, PI controller of gain link _p And integral coefficient K _I 。

Optionally, the functional expression of the primary steam pressure model pre-established in the primary frequency modulation simulation program unit is:

in the above-mentioned method, the step of,p _T is the main steam pressure;Ais the amplitude;fis frequency;ttime is;θ ₀ is the phase angle;bis biased.

Optionally, the pre-established turbine model in the primary frequency modulation simulation program unit includes:

a valve flow characteristic link for outputting a control system model of the control system model by using a preset valve flow characteristic functionf(x) Processing to obtain a first intermediate value;

the multiplication module is used for multiplying the main steam pressure output by the main steam pressure model by a first intermediate value;

a regulating stage pressure-load link for applying a preset regulating stage pressure-load function to the output of the multiplication moduleg(x) Processing to obtain a second intermediate value;

a transfer function module for adopting a preset transfer function for the second intermediate valueG(s) Generating a third intermediate value;

a summing module for adding a third intermediate value, the input steady state power response deltaAnd generating active power P output by the steam turbine model.

Optionally, the preset valve flow characteristic functionf(x) The abscissa is the output of the control system model and the ordinate is the pressure ratio for the first intermediate value.

Optionally, the preset regulation stage pressure-load functiong(x) The abscissa of (2) is the regulation stage pressure and the ordinate is the load.

Optionally, the preset regulation stage pressure-load functionG(s) The functional expression of (2) is:

，

in the above-mentioned method, the step of,F _HP is the power ratio of the high-pressure cylinder,λis the power overshoot factor of the high-pressure cylinder,T _CH is the volume time constant of the high-pressure steam chamber,T _RH is the reheater volume time constant.

Optionally, the steady state power response delta of the inputThe functional expression of (2) is:

in the above-mentioned method, the step of,mfor the number of phase angles,is a steady-state power amplitude characteristic parameter +>I element of->Is a steady-state power frequency characteristic parameter->Is selected from the group consisting of the (i) th element,ttime is; />Is the firstiPhase angle, & gt>，randIs a random number between 0 and 1; wherein the steady-state power amplitude characteristic parameter->And steady state power frequency characteristic parameterF _f The acquisition of (1) comprises:

s101, extracting steady-state load data of the unit within a certain period of time，/>，/>Is the firstiSteady state load data at a point in time,nis steady-state load data length, unit steady-state load data +.>Time interval of adjacent elements->Equal, corresponding sampling rate is +.>；

S102, calculating steady-state load variables of the unit according to the following formula：

In the above-mentioned method, the step of,for steady state load data of the unit- >Is the average value of (2); />；

S103, obtaining an intermediate variable through discrete Fourier calculation，/>Wherein, and have:

，

in the above-mentioned method, the step of,for intermediate variable +.>I element of->Is->The first of (3)kThe number of elements to be added to the composition,jin units of imaginary numbers,；

s104, calculating steady-state power amplitude characteristic parameters according to the following，/>Elements of (a) and (b):

，/>，/>，

in the above-mentioned method, the step of,is a steady-state power amplitude characteristic parameter +>Element 1 of->Is a steady-state power amplitude characteristic parameter +>To the first part of (1)l2) +1 elements, ++>Is a steady-state power amplitude characteristic parameter +>At the time of the mth element of (1)mThe value range of (2) is +.>，/>For intermediate variable +.>Element 1 of->For intermediate variable +.>To the first part of (1)l2) +1 elements; calculating steady-state power frequency characteristic parameters according to the following formula>，/>：

，

In the above-mentioned method, the step of,is a steady-state power frequency characteristic parameter->M element of->Is the sampling rate.

Optionally, when the primary frequency modulation simulation program unit performs primary frequency modulation simulation by combining a control system model with a pre-established main steam pressure model and a steam turbine model to obtain a total primary frequency modulation action frequency and a primary frequency modulation qualified frequency, the primary frequency modulation simulation program unit includes:

s201, acquiring a power grid frequency fluctuation data set，/>Wherein the firstiIndividual element->，Wherein->The first of (3)jIndividual element->Represents- >Grid frequency at time,/->Is the grid frequency sampling interval; initializing the total action times of primary frequency modulation and the qualified times of primary frequency modulation;

s202, randomly selecting a power grid frequency fluctuation data setAny elementF _i As the input of the control system model, the active power P output by the steam turbine model is generated by carrying out single primary frequency modulation simulation by combining the control system model with a pre-established main steam pressure model and the steam turbine model;

s203, by elementF _i And calculating the actual integral contribution electric quantity of the simulation according to the active power P output by the steam turbine modelH _i And theoretical integral contribution powerH _e And will actually integrate the contributing powerH _i Divided by theoretical integral contributionH _e Obtaining the primary frequency modulation contribution power rate of the simulationkThe method comprises the steps of carrying out a first treatment on the surface of the Adding 1 to the total action frequency of the primary frequency modulation, if the primary frequency modulation contributes to the electric quantity ratekIs more than or equal to the set value,then 1 is added to the qualified frequency of primary frequency modulation;

and S204, judging whether the total frequency modulation action times are greater than or equal to a preset value, if so, judging that the total frequency modulation action times and the qualified frequency modulation times are obtained by carrying out primary frequency modulation simulation through a control system model in combination with a pre-established main steam pressure model and a steam turbine model, otherwise, continuing to iterate in the step S202.

Optionally, the actual integrated contribution powerH _i The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ the system frequency exceeds the dead zone of primary frequency modulation action of the unit,t _t the moment when the system frequency enters the dead zone of the primary frequency modulation action of the unit,P _t is thattThe time unit actually generates active power,P ₀ is thatt ₀ And the time unit actually generates active power, and t is time.

Optionally, the theoretical integration contribution powerH _e The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ for the moment when the system frequency exceeds the dead zone of the primary frequency modulation action of the unit,t _t the time when the system frequency enters the primary frequency modulation action dead zone of the unit,is thattAt the moment when the power grid frequency exceeds the frequency modulation artificial dead zone value, < ->Rated active output for the unit, < >>Rated frequency for the power grid, < >>The speed change rate of the unit is set, and t is time.

Optionally, the step of calculating the primary frequency modulation qualified rate according to the total number of primary frequency modulation actions and the primary frequency modulation qualified number in the simulation result judging program unit refers to dividing the primary frequency modulation qualified number by the total number of primary frequency modulation actions to obtain the primary frequency modulation qualified rate.

In addition, the invention also provides a model simulation-based primary frequency modulation parameter adjustment system, which comprises a microprocessor and a memory which are connected with each other, wherein the microprocessor is programmed or configured to execute the model simulation-based primary frequency modulation parameter adjustment method.

Furthermore, the invention provides a computer readable storage medium having a computer program stored therein, the computer program being for being programmed or configured by a microprocessor to perform the model-simulation-based primary frequency modulation parameter adjustment method.

Compared with the prior art, the invention has the following advantages: determining a primary frequency modulation parameter to be adjusted in a pre-established control system model, and determining an empirical value of the primary frequency modulation parameter as a current setting value of the primary frequency modulation parameter; the method comprises the steps of assigning a current set value of a primary frequency modulation parameter to be adjusted to a control system model, carrying out primary frequency modulation simulation by combining a pre-established main steam pressure model and a steam turbine model through the control system model to obtain total primary frequency modulation action times and qualified primary frequency modulation times, calculating the primary frequency modulation qualification rate according to the total primary frequency modulation action times and the qualified primary frequency modulation times, and iteratively adjusting the primary frequency modulation parameter based on the primary frequency modulation qualification rate obtained by simulation, so that the primary frequency modulation parameter of a unit can better overcome the influence of load and spontaneous fluctuation of main steam pressure of the thermal power unit, and finally, the thermal power unit has better primary frequency modulation power response performance.

Drawings

FIG. 1 is a schematic diagram of a basic flow of a method according to an embodiment of the present invention.

Fig. 2 is a schematic structural diagram of a control system model according to an embodiment of the present invention.

Fig. 3 shows the frequency spectrum analysis result in the embodiment of the invention.

FIG. 4 is a schematic view of a turbine model in accordance with an embodiment of the present invention.

FIG. 5 is a graph showing the flow characteristics of a valve according to an embodiment of the present inventionf(x) Is a schematic diagram of the curve of (a).

FIG. 6 is a graph of the pressure-load function of the conditioning stage in an embodiment of the inventiong(x) Is a schematic diagram of the curve of (a).

FIG. 7 is a diagram of a grid frequency fluctuation dataset element F in an embodiment of the invention ₁ Is a schematic diagram of the curve of (a).

FIG. 8 is a grid frequency fluctuation dataset element F in an embodiment of the invention ₂ Is a schematic diagram of the curve of (a).

FIG. 9 is a diagram of a grid frequency fluctuation dataset element F in an embodiment of the invention ₃ Is a schematic diagram of the curve of (a).

FIG. 10 is a diagram of a grid frequency fluctuation dataset element F in an embodiment of the invention ₄ Is a schematic diagram of the curve of (a).

FIG. 11 is a diagram of a grid frequency fluctuation dataset element F in an embodiment of the invention ₅ Is a schematic diagram of the curve of (a).

FIG. 12 is a grid frequency fluctuation dataset element F in an embodiment of the invention ₆ Is a schematic diagram of the curve of (a).

FIG. 13 is a diagram of a grid frequency fluctuation dataset element F in an embodiment of the invention ₇ Is a schematic diagram of the curve of (a).

FIG. 14 is a grid frequency fluctuation dataset element F in an embodiment of the invention ₈ Is a schematic diagram of the curve of (a).

FIG. 15 is a diagram of a grid frequency fluctuation dataset element F in an embodiment of the invention ₉ Is a schematic diagram of the curve of (a).

FIG. 16 is a diagram of a grid frequency fluctuation dataset element F in an embodiment of the invention ₁₀ Is a schematic diagram of the curve of (a).

Detailed Description

The method for adjusting the primary frequency modulation parameters based on model simulation according to the invention is further described in detail below by taking a certain 300MW subcritical thermal power generating unit as an implementation object.

As shown in fig. 1, the primary frequency modulation parameter adjustment method based on model simulation in this embodiment includes:

s1, determining a primary frequency modulation parameter to be adjusted in a pre-established control system model, and determining an empirical value of the primary frequency modulation parameter as a current setting value of the primary frequency modulation parameter;

s2, assigning the current set value of the primary frequency modulation parameter to be adjusted to a control system model, carrying out primary frequency modulation simulation by combining a pre-established main steam pressure model and a steam turbine model through the control system model to obtain the total number of primary frequency modulation actions and the qualified number of primary frequency modulation, and calculating the qualified rate of primary frequency modulation according to the total number of primary frequency modulation actions and the qualified number of primary frequency modulation;

s3, judging whether the primary frequency modulation qualification rate meets the requirement, if not, correcting the current set value of the primary frequency modulation parameter to be adjusted, and jumping to the step S2; otherwise, ending and exiting.

As shown in fig. 2, the control system model of primary frequency modulation of the thermal power generating unit, which is pre-established in step S1 of the present embodiment, includes:

the rotating speed unequal rate function module is used for forming a frequency modulation instruction by the input power grid frequency through a rotating speed unequal rate function;

the controller set value generation module is used for adding the frequency modulation instruction and the load setting to form a controller set value;

the load deviation calculation module is used for feeding back the set value of the controller and the input load to calculate the load deviation;

the closed loop controller is used for generating a controller output according to the load deviation;

the gain link is used for taking the frequency modulation instruction as frequency modulation feedforward after gain;

and the summation module is used for adding the controller output and the frequency modulation feedforward or directly outputting the controller output and the frequency modulation feedforward as the control system model output.

The working process of the control system model is as follows: forming a frequency modulation instruction by the power grid frequency through a rotating speed unequal rate link; adding the frequency modulation instruction and the load setting to form a controller setting value; the controller receives a controller set value and load feedback operation to form a controller output; forming frequency modulation feedforward after the frequency modulation instruction passes through a gain link; the controller output may be added to the frequency modulated feedforward or directly output to form a control system model output.

It should be noted that the closed-loop controller may employ a desired closed-loop control algorithm as needed. For example, as an alternative implementation manner, the closed-loop controller in this embodiment is a PI controller.

The primary frequency modulation parameters to be adjusted in the control system model determined in step S1 of this embodiment include: gain factor K, PI controller of gain link _p And integral coefficient K _I . In this embodiment, the gain factor of the gain element K, PI is the scaling factor K of the controller _p And integral coefficient K _I The initial value is set empirically, specifically, the gain factor is set initially empiricallyK=0.5, scaling factorK _P =0.4, integral coefficientK _I =0.01。

The main steam pressure model is used for generating simulated main steam pressure, and a function can be adopted to generate the required main steam pressure or an actual wave recording signal of the main steam pressure according to the requirement. The functional expression of the pre-established main steam pressure model in step S2 of this embodiment is:

in the above-mentioned method, the step of,p _T is the main steam pressure;Ais the amplitude;fis frequency;ttime is;θ ₀ is the phase angle;bis biased. Taking typical main steam pressure fluctuation data of a thermal power generating unit under a stable operation condition, carrying out Fourier frequency amplitude analysis on the data, and respectively representing the maximum amplitude value and the corresponding frequency value in the frequency amplitude analysis AAndfthe method comprises the steps of carrying out a first treatment on the surface of the Bias ofbManually setting a value within a main steam pressure range under a normal operation condition of the thermal power generating unit;θ ₀ randomly taking values within the range of 0-2 pi. The frequency amplitude analysis in this embodiment is shown in FIG. 3, where the maximum amplitude isA= 0.1806, corresponding frequencyValue off= 0.002857Hz; bias ofbManually set to 12;θ ₀ and randomly taking values within the range of 0-2 pi.

In this embodiment, step S2 assigns the current setting value of the primary frequency modulation parameter to be adjusted to the control system model, and performs primary frequency modulation simulation by combining the control system model with the pre-established main steam pressure model and the steam turbine model, specifically, by using commercial software Matlab/Simulink, and may also use other simulation software as required.

As shown in fig. 4, the turbine model previously established in step S2 of the present embodiment includes:

a valve flow characteristic link for outputting a control system model of the control system model by using a preset valve flow characteristic functionf(x) Processing to obtain a first intermediate value;

the multiplication module is used for multiplying the main steam pressure output by the main steam pressure model by a first intermediate value;

a regulating stage pressure-load link for applying a preset regulating stage pressure-load function to the output of the multiplication module g(x) Processing to obtain a second intermediate value;

a transfer function module for adopting a preset transfer function for the second intermediate valueG(s) Generating a third intermediate value;

a summing module for adding a third intermediate value, the input steady state power response deltaAnd generating active power P output by the steam turbine model.

The working process of the steam turbine model is as follows: the input of the steam turbine model is respectively the output of the control system model and the output of the main steam pressure model, and the output of the control system model forms a first intermediate value after passing through a valve flow characteristic link; the first intermediate value is multiplied by the main steam pressure output and then forms a second intermediate value after passing through a regulating stage pressure-load link; the second intermediate value passes through a transfer functionG(s) Form a third intermediate value and respond to the input steady-state power by increasingAnd generating active power P output by the steam turbine model.

In this embodiment, the preset valve flow characteristic functionf(x) The abscissa of (2) is the output of the control system model and the ordinate is the pressure ratio for the first intermediate value, as shown in fig. 5.

In this embodiment, the preset regulation stage pressure-load functiong(x) The abscissa of (2) is the regulation stage pressure and the ordinate is the load, as shown in fig. 6.

In this embodiment, the preset regulation stage pressure-load functionG(s) The functional expression of (2) is:

，

in the above-mentioned method, the step of,F _HP is the power ratio of the high-pressure cylinder,λis the power overshoot factor of the high-pressure cylinder,T _CH is the volume time constant of the high-pressure steam chamber,T _RH is the reheater volume time constant. Parameters in the transfer function are determined according to actual measurement data of parameters of a steam turbine and a regulating system of the unit, and in the embodimentF _HP =0.3，λ=0.8，T _CH =0.1，T _RH =10.5。

In this embodiment, the steady state power response delta of the inputThe functional expression of (2) is:

in the above-mentioned method, the step of,mfor the number of phase angles,is a steady-state power amplitude characteristic parameter +>I element of->Is a steady-state power frequency characteristic parameter->Is selected from the group consisting of the (i) th element,ttime is; />Is the firstiPhase angle, & gt>，randIs a random number between 0 and 1; the summation module adds the third intermediate value, the input steady-state power response increment>The active power P that generates the turbine model output can be expressed as: />。

In this embodiment, the steady-state power amplitude characteristic parameterAnd steady state power frequency characteristic parameterF _f The acquisition of (1) comprises:

s101, extracting steady-state load data of the unit within a certain period of time，/>，/>Is the firstiSteady state load data at a point in time,nis steady-state load data length, unit steady-state load data +.>Time interval of adjacent elements- >Equal, corresponding sampling rate is +.>；

S102, calculating steady-state load variables of the unit according to the following formula：

In the above-mentioned method, the step of,for steady state load data of the unit->Is the average value of (2); />；

S103, obtaining an intermediate variable through discrete Fourier calculation，/>Wherein, and have:

，

in the above-mentioned method, the step of,for intermediate variable +.>I element of->Is->The first of (3)kThe number of elements to be added to the composition,jin units of imaginary numbers,；

s104, calculating steady-state power amplitude characteristic parameters according to the following，/>Elements of (a) and (b):

，/>，/>，

in the above-mentioned method, the step of,is a steady-state power amplitude characteristic parameter +>Element 1 of->Is a steady-state power amplitude characteristic parameter +>To the first part of (1)l2) +1 elements, ++>Is a steady-state power amplitude characteristic parameter +>At the time of the mth element of (1)mThe value range of (2) is +.>，/>For intermediate variable +.>Element 1 of->For intermediate variable +.>To the first part of (1)l2) +1 elements; calculating steady-state power frequency characteristic parameters according to the following formula>，/>：

，

In the above-mentioned method, the step of,is a steady-state power frequency characteristic parameter->M element of->Is the sampling rate.

In this embodiment, when performing primary frequency modulation simulation by combining the control system model with the pre-established main steam pressure model and the steam turbine model to obtain the total number of primary frequency modulation actions and the qualified number of primary frequency modulation in step S2, the method includes:

s201, acquiring a power grid frequency fluctuation data set ，/>Wherein the firstiIndividual element->，Wherein->The first of (3)jIndividual element->Represents->Grid frequency at time,/->Is the grid frequency sampling interval; initializing the total action times of primary frequency modulation and the qualified times of primary frequency modulation; in this embodiment, n=10, the grid frequency fluctuation dataset element F ₁ -F ₁₀ As shown in fig. 7-16, respectively.

S202, randomly selecting a power grid frequency fluctuation data setAny elementF _i As the input of the control system model, the active power P output by the steam turbine model is generated by carrying out single primary frequency modulation simulation by combining the control system model with a pre-established main steam pressure model and the steam turbine model; />

S203, by elementF _i And calculating the actual integral contribution electric quantity of the simulation according to the active power P output by the steam turbine modelH _i And theoretical integral contribution powerH _e And will actually integrate the contributing powerH _i Divided by theoretical integral contributionH _e Obtaining the primary frequency modulation contribution power rate of the simulationkThe method comprises the steps of carrying out a first treatment on the surface of the Adding 1 to the total action frequency of the primary frequency modulation, if the primary frequency modulation contributes to the electric quantity ratekIf the frequency is larger than or equal to the set value, adding 1 to the qualified frequency of primary frequency modulation;

and S204, judging whether the total frequency modulation action times are greater than or equal to a preset value, if so, judging that the total frequency modulation action times and the qualified frequency modulation times are obtained by carrying out primary frequency modulation simulation through a control system model in combination with a pre-established main steam pressure model and a steam turbine model, otherwise, continuing to iterate in the step S202.

In the present embodiment, the actual integrated contribution powerH _i The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ the system frequency exceeds the dead zone of primary frequency modulation action of the unit,t _t the moment when the system frequency enters the dead zone of the primary frequency modulation action of the unit,P _t is thattThe time unit actually generates active power,P ₀ is thatt ₀ And the time unit actually generates active power, and t is time.

In the present embodiment, the theoretical integration contributes to the electric quantityH _e The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ for the moment when the system frequency exceeds the dead zone of the primary frequency modulation action of the unit,t _t the time when the system frequency enters the primary frequency modulation action dead zone of the unit,is thattAt the moment when the power grid frequency exceeds the frequency modulation artificial dead zone value, < ->Rated active output for the unit, < >>Rated frequency (50 Hz) for the power grid, < >>The speed change rate of the unit is set, and t is time.

In this embodiment, in step S3, calculating the qualified primary frequency according to the total number of primary frequency modulation actions and the qualified primary frequency modulation actions refers to dividing the qualified primary frequency modulation actions by the total number of primary frequency modulation actions to obtain the qualified primary frequency modulation Q, which may be expressed as:

，

in this embodiment, when 1 primary frequency modulation occurs, the total number of primary frequency modulation is increased by 1; when the primary frequency modulation contributing power rate k is more than or equal to 50 during primary frequency modulation action, the number of times of primary frequency modulation qualification is increased by 1. A total of 200 sets of primary frequency modulation simulation data are obtained, and according to the primary frequency modulation qualified frequency=117, the total primary frequency modulation operation frequency=200, so that the primary frequency modulation qualified rate is obtained Q=58.5%。

Finally, in this embodiment, 200 primary frequency simulations are performed to obtain 200 sets of primary frequency simulation data, and according to the primary frequency qualified frequency=189, the total primary frequency action frequency=200, so that the primary frequency qualified rate is obtainedQ=94.5%. At this timeQAnd the qualification rate is more than 90%, the qualification rate threshold requirement of the 'two rules' of the Huazhong power grid for small disturbance examination-free is met, and the primary frequency modulation qualification rate meets the expected requirement. Therefore, the value of the primary frequency modulation parameter is not adjusted any more and repeated iterative simulation is not needed.

To verify the method of the present embodiment, the gain element coefficients obtained by the empirical method are used in the present embodimentK（K=1.5), ratio coefficientK _P （K _P =0.4), integral coefficientK _I （K _I =0.1) is set as the actual main control corresponding parameter of the turbine, and a small disturbance primary frequency modulation qualification rate test method, medium and system are used (publication number: 115436736A) after the test parameters of the disclosed primary frequency modulation active test system are adjusted, the actual primary frequency modulation qualification rate of the unit is calculated to carry out the primary frequency modulation test 174 times, and the primary frequency modulation contribution power rate is calculatedkThe frequency of more than or equal to 50 is 161 times, namely primary frequency modulationThe qualification rate is 92.5%, and the test result shows that the simulation qualification rate obtained by the primary frequency modulation parameter adjustment method based on the model design can be well matched with the primary frequency modulation qualification rate of the actual unit, and has a good guiding effect on the primary frequency modulation parameter adjustment of the actual unit.

In summary, according to the primary frequency modulation parameter adjustment method based on the model simulation, the model for primary frequency modulation parameter adjustment is established, the primary frequency modulation optimization index is set, mass primary frequency modulation simulation under typical working conditions is developed based on the design model, and the primary frequency modulation parameters are iteratively adjusted based on the simulation result, so that the primary frequency modulation parameters of the unit can better overcome the influence of spontaneous fluctuation of the load and the main steam pressure of the thermal power unit, and finally the thermal power unit has better primary frequency modulation power response performance. According to the primary frequency modulation parameter adjustment method based on model simulation, a model capable of simulating spontaneous fluctuation of power of the thermal power generating unit is established, and accuracy of simulation results is improved; the primary frequency modulation optimization index based on the primary frequency modulation qualification rate is designed, and the pertinency of parameter optimization adjustment is improved; the power grid frequency fluctuation data set is established, and the influence of the power grid frequency fluctuation characteristic on the simulation result is more truly considered by randomly selecting any element of the power grid frequency fluctuation data set as the power grid frequency. The primary frequency modulation parameter adjustment method based on the model design realizes the optimal adjustment of the primary frequency modulation parameter more comprehensively through massive and real simulation calculation, and reduces the time required by the optimal adjustment of the conventional method while obtaining better adjustment effect.

In addition, the embodiment also provides a primary frequency modulation parameter adjustment system based on model simulation, which comprises:

the primary frequency modulation parameter selecting and value taking program unit is used for determining primary frequency modulation parameters to be adjusted in a pre-established control system model and determining the experience value of the primary frequency modulation parameters as the current setting value of the primary frequency modulation parameters;

the primary frequency modulation simulation program unit is used for assigning the current set value of the primary frequency modulation parameter to be adjusted to the control system model, carrying out primary frequency modulation simulation by combining the control system model with a pre-established main steam pressure model and a steam turbine model to obtain the total number of primary frequency modulation actions and the qualified number of primary frequency modulation, and calculating the qualified rate of primary frequency modulation according to the total number of primary frequency modulation actions and the qualified number of primary frequency modulation;

the simulation result judging program unit is used for judging whether the primary frequency modulation qualification rate meets the requirement, if not, correcting the current set value of the primary frequency modulation parameter to be adjusted, and calling the primary frequency modulation simulation program unit to continuously execute a round of primary frequency modulation simulation; otherwise, ending and exiting.

In addition, the embodiment also provides a primary frequency modulation parameter adjustment system based on model simulation, which comprises a microprocessor and a memory which are connected with each other, wherein the microprocessor is programmed or configured to execute the primary frequency modulation parameter adjustment method based on model simulation.

In addition, the present embodiment also provides a computer readable storage medium having a computer program stored therein, the computer program being configured or programmed by a microprocessor to perform the model-simulation-based primary frequency modulation parameter adjustment method.

It will be appreciated by those skilled in the art that 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-readable 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 flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations 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.

The above description is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above examples, and all technical solutions belonging to the concept of the present invention belong to the protection scope of the present invention. It should be noted that modifications and adaptations to the present invention may occur to one skilled in the art without departing from the principles of the present invention and are intended to be within the scope of the present invention.

Claims (30)

1. The primary frequency modulation parameter adjustment method based on the model simulation is characterized by comprising the following steps of:

s1, determining a primary frequency modulation parameter to be adjusted in a pre-established control system model, and determining an empirical value of the primary frequency modulation parameter as a current setting value of the primary frequency modulation parameter;

s2, assigning the current set value of the primary frequency modulation parameter to be adjusted to a control system model, carrying out primary frequency modulation simulation by combining a pre-established main steam pressure model and a steam turbine model through the control system model to obtain the total number of primary frequency modulation actions and the qualified number of primary frequency modulation, and calculating the qualified rate of primary frequency modulation according to the total number of primary frequency modulation actions and the qualified number of primary frequency modulation;

s3, judging whether the primary frequency modulation qualification rate meets the requirement, if not, correcting the current set value of the primary frequency modulation parameter to be adjusted, and jumping to the step S2; otherwise, ending and exiting.

2. The method for adjusting primary frequency modulation parameters based on model simulation according to claim 1, wherein the control system model of primary frequency modulation of the thermal power generating unit, which is pre-established in step S1, comprises:

the rotating speed unequal rate function module is used for forming a frequency modulation instruction by the input power grid frequency through a rotating speed unequal rate function;

the controller set value generation module is used for adding the frequency modulation instruction and the load setting to form a controller set value;

the load deviation calculation module is used for feeding back the set value of the controller and the input load to calculate the load deviation;

the closed loop controller is used for generating a controller output according to the load deviation;

the gain link is used for taking the frequency modulation instruction as frequency modulation feedforward after gain;

and the summation module is used for adding the controller output and the frequency modulation feedforward or directly outputting the controller output and the frequency modulation feedforward as the control system model output.

3. The model simulation-based primary frequency modulation parameter adjustment method according to claim 2, wherein the closed-loop controller is a PI controller.

4. The method for adjusting a primary frequency modulation parameter based on model simulation according to claim 3, wherein the primary frequency modulation parameter to be adjusted in the control system model determined in step S1 comprises: gain factor K, PI controller of gain link _p And integral coefficient K _I 。

5. The method for adjusting primary frequency modulation parameters based on model simulation according to claim 1, wherein the functional expression of the primary steam pressure model pre-established in step S2 is:

in the above-mentioned method, the step of,p _T is the main steam pressure;Ais the amplitude;fis frequency;ttime is;θ ₀ is the phase angle;bis biased.

6. The model simulation-based primary frequency modulation parameter adjustment method according to claim 1, wherein the pre-established turbine model in step S2 comprises:

a valve flow characteristic link for outputting a control system model of the control system model by using a preset valve flow characteristic functionf(x) Processing to obtain a first intermediate value;

the multiplication module is used for multiplying the main steam pressure output by the main steam pressure model by a first intermediate value;

a regulating stage pressure-load link for applying a preset regulating stage pressure-load function to the output of the multiplication moduleg(x) Processing to obtain a second intermediate value;

a transfer function module for adopting a preset transfer function for the second intermediate valueG(s) Generating a third intermediate value;

a summing module for adding a third intermediate value, the input steady state power response deltaAnd generating active power P output by the steam turbine model.

7. The model simulation-based primary frequency modulation parameter adjustment method according to claim 6, wherein the preset valve flow characteristic functionf(x) The abscissa is the output of the control system model, and the ordinate is the output of the control system modelThe coordinates are for the pressure ratio as the first intermediate value.

8. The model simulation based primary frequency modulation parameter adjustment method according to claim 6, wherein the preset regulation level pressure-load functiong(x) The abscissa of (2) is the regulation stage pressure and the ordinate is the load.

9. The model simulation based primary frequency modulation parameter adjustment method according to claim 6, wherein the preset regulation level pressure-load functionG(s) The functional expression of (2) is:

，

in the above-mentioned method, the step of,F _HP is the power ratio of the high-pressure cylinder,λis the power overshoot factor of the high-pressure cylinder,T _CH is the volume time constant of the high-pressure steam chamber,T _RH is the reheater volume time constant.

10. The model simulation based primary frequency modulation parameter adjustment method according to claim 6, wherein the input steady-state power response incrementThe functional expression of (2) is:

in the above-mentioned method, the step of,mfor the number of phase angles,is a steady-state power amplitude characteristic parameter +>I element of->Is a steady-state power frequency characteristic parameter- >Is selected from the group consisting of the (i) th element,ttime is; />Is the firstiPhase angle, & gt>，randIs a random number between 0 and 1; wherein the steady-state power amplitude characteristic parameter->And steady state power frequency characteristic parameterF _f The acquisition of (1) comprises:

s101, extracting steady-state load data of the unit within a certain period of time，/>，/>Is the firstiSteady state load data at a point in time,nis steady-state load data length, unit steady-state load data +.>Time interval of adjacent elements->Equal, corresponding sampling rate is +.>；

S102, calculating steady-state load variables of the unit according to the following formula：

In the above-mentioned method, the step of,for steady state load data of the unit->Is the average value of (2); />；

S103, obtaining an intermediate variable through discrete Fourier calculation，/>Wherein, and have:

，

in the above-mentioned method, the step of,for intermediate variable +.>I element of->Is->The first of (3)kThe number of elements to be added to the composition,jin units of imaginary numbers,；

s104, calculating steady-state power amplitude characteristic parameters according to the following，/>Elements of (a) and (b):

，/>，/>，

in the above-mentioned method, the step of,is a steady-state power amplitude characteristic parameter +>Element 1 of->Is a steady-state power amplitude characteristic parameter +>To the first part of (1)l2) +1 elements, ++>Is a steady-state power amplitude characteristic parameter +>At the time of the mth element of (1)mThe range of the values is as follows，/>For intermediate variable +.>Element 1 of->For intermediate variable +.>To the first part of (1)l2) +1 elements; calculating steady-state power frequency characteristic parameters according to the following formula >，/>：

，

In the above-mentioned method, the step of,is a steady-state power frequency characteristic parameter->M element of->Is the sampling rate.

11. The method for adjusting a primary frequency modulation parameter based on model simulation according to claim 1, wherein in step S2, when performing primary frequency modulation simulation by combining a control system model with a pre-established main steam pressure model and a steam turbine model to obtain a total number of primary frequency modulation actions and a qualified number of primary frequency modulation actions, the method comprises:

s201, acquiring a power grid frequency fluctuation data set，/>Wherein the firstiIndividual element->，/>Wherein->The first of (3)jIndividual element->Represents->Grid frequency at time,/->Is the grid frequency sampling interval; initializing the total action times of primary frequency modulation and the qualified times of primary frequency modulation;

s202, randomly selecting a power grid frequency fluctuation data setAny elementF _i As the input of the control system model, the active power P output by the steam turbine model is generated by carrying out single primary frequency modulation simulation by combining the control system model with a pre-established main steam pressure model and the steam turbine model;

s203, by elementF _i And calculating the actual integral contribution electric quantity of the simulation according to the active power P output by the steam turbine modelH _i And theory ofIntegral contribution powerH _e And will actually integrate the contributing power H _i Divided by theoretical integral contributionH _e Obtaining the primary frequency modulation contribution power rate of the simulationkThe method comprises the steps of carrying out a first treatment on the surface of the Adding 1 to the total action frequency of the primary frequency modulation, if the primary frequency modulation contributes to the electric quantity ratekIf the frequency is larger than or equal to the set value, adding 1 to the qualified frequency of primary frequency modulation;

and S204, judging whether the total frequency modulation action times are greater than or equal to a preset value, if so, judging that the total frequency modulation action times and the qualified frequency modulation times are obtained by carrying out primary frequency modulation simulation through a control system model in combination with a pre-established main steam pressure model and a steam turbine model, otherwise, continuing to iterate in the step S202.

12. The model simulation-based primary frequency modulation parameter adjustment method according to claim 11, wherein the actual integrated contribution power isH _i The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ the system frequency exceeds the dead zone of primary frequency modulation action of the unit,t _t the moment when the system frequency enters the dead zone of the primary frequency modulation action of the unit,P _t is thattThe time unit actually generates active power,P ₀ is thatt ₀ And the time unit actually generates active power, and t is time.

13. The model simulation-based primary frequency modulation parameter adjustment method according to claim 11, wherein the theoretical integration contribution electric quantity H _e The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ for the moment when the system frequency exceeds the dead zone of the primary frequency modulation action of the unit,t _t the time when the system frequency enters the primary frequency modulation action dead zone of the unit,is thattAt the moment when the power grid frequency exceeds the frequency modulation artificial dead zone value, < ->Rated active output for the unit, < >>Rated frequency for the power grid, < >>The speed change rate of the unit is set, and t is time.

14. The method for adjusting a primary frequency modulation parameter based on model simulation as claimed in claim 1, wherein the step S3 of calculating the primary frequency modulation qualification rate according to the total number of primary frequency modulation actions and the number of primary frequency modulation qualification times is to divide the number of primary frequency modulation qualification times by the total number of primary frequency modulation actions to obtain the primary frequency modulation qualification rate.

15. A model simulation-based primary frequency modulation parameter adjustment system, comprising:

the primary frequency modulation parameter selecting and value taking program unit is used for determining primary frequency modulation parameters to be adjusted in a pre-established control system model and determining the experience value of the primary frequency modulation parameters as the current setting value of the primary frequency modulation parameters;

the primary frequency modulation simulation program unit is used for assigning the current set value of the primary frequency modulation parameter to be adjusted to the control system model, carrying out primary frequency modulation simulation by combining the control system model with a pre-established main steam pressure model and a steam turbine model to obtain the total number of primary frequency modulation actions and the qualified number of primary frequency modulation, and calculating the qualified rate of primary frequency modulation according to the total number of primary frequency modulation actions and the qualified number of primary frequency modulation;

The simulation result judging program unit is used for judging whether the primary frequency modulation qualification rate meets the requirement, if not, correcting the current set value of the primary frequency modulation parameter to be adjusted, and calling the primary frequency modulation simulation program unit to continuously execute a round of primary frequency modulation simulation; otherwise, ending and exiting.

16. The model emulation-based primary frequency modulation parameter adjustment system according to claim 15, wherein the primary frequency modulation parameter selection and value-taking program unit pre-establishes a primary frequency modulation control system model of the thermal power generating unit, comprising:

the rotating speed unequal rate function module is used for forming a frequency modulation instruction by the input power grid frequency through a rotating speed unequal rate function;

the controller set value generation module is used for adding the frequency modulation instruction and the load setting to form a controller set value;

the load deviation calculation module is used for feeding back the set value of the controller and the input load to calculate the load deviation;

the closed loop controller is used for generating a controller output according to the load deviation;

the gain link is used for taking the frequency modulation instruction as frequency modulation feedforward after gain;

and the summation module is used for adding the controller output and the frequency modulation feedforward or directly outputting the controller output and the frequency modulation feedforward as the control system model output.

17. The model simulation based primary frequency modulation parameter adjustment system of claim 16, wherein the closed loop controller is a PI controller.

18. The model emulation-based chirp parameter tuning system of claim 17, wherein the chirp parameter to be tuned in the control system model determined by the chirp parameter selection and valuation program unit comprises: gain factor K, PI controller of gain link _p And integral coefficient K _I 。

19. The model simulation-based primary frequency modulation parameter adjustment system according to claim 15, wherein the primary frequency modulation simulation program unit has a function expression of a primary steam pressure model established in advance, which is:

in the above-mentioned method, the step of,p _T is the main steam pressure;Ais the amplitude;fis frequency;ttime is;θ ₀ is the phase angle;bis biased.

20. The model simulation based chirp parameter tuning system of claim 15, wherein the pre-established turbine model in the chirp simulation program unit comprises:

a valve flow characteristic link for outputting a control system model of the control system model by using a preset valve flow characteristic function f(x) Processing to obtain a first intermediate value;

the multiplication module is used for multiplying the main steam pressure output by the main steam pressure model by a first intermediate value;

a regulating stage pressure-load link for applying a preset regulating stage pressure-load function to the output of the multiplication moduleg(x) Processing to obtain a second intermediate value;

a transfer function module for adopting a preset transfer function for the second intermediate valueG(s) Generating a third intermediate value;

a summing module for adding a third intermediate value, the input steady state power response deltaAnd generating active power P output by the steam turbine model.

21. Model simulation based primary frequency modulation parameter tuning in accordance with claim 20The whole system is characterized in that the preset valve flow characteristic functionf(x) The abscissa is the output of the control system model and the ordinate is the pressure ratio for the first intermediate value.

22. The model simulation based chirping parameter tuning system of claim 20 wherein the preset tuning stage pressure-load functiong(x) The abscissa of (2) is the regulation stage pressure and the ordinate is the load.

23. The model simulation based chirping parameter tuning system of claim 20 wherein the preset tuning stage pressure-load function G(s) The functional expression of (2) is:

，

in the above-mentioned method, the step of,F _HP is the power ratio of the high-pressure cylinder,λis the power overshoot factor of the high-pressure cylinder,T _CH is the volume time constant of the high-pressure steam chamber,T _RH is the reheater volume time constant.

24. The model simulation based primary frequency modulation parameter tuning system of claim 20, wherein the input steady state power response deltaThe functional expression of (2) is:

in the above-mentioned method, the step of,mfor the number of phase angles,is a steady-state power amplitude characteristic parameter +>I element of->Is a steady-state power frequency characteristic parameter->Is selected from the group consisting of the (i) th element,ttime is; />Is the firstiPhase angle, & gt>，randIs a random number between 0 and 1; wherein the steady-state power amplitude characteristic parameter->And steady state power frequency characteristic parameterF _f The acquisition of (1) comprises:

s101, extracting steady-state load data of the unit within a certain period of time，/>，/>Is the firstiSteady state load data at a point in time,nis steady-state load data length, unit steady-state load data +.>Time interval of adjacent elements->Equal, corresponding sampling rate is +.>；

S102, calculating steady-state load variables of the unit according to the following formula：

In the above-mentioned method, the step of,for steady state load data of the unit->Is the average value of (2); />；

S103, obtaining an intermediate variable through discrete Fourier calculation，/>Wherein, and have:

，

In the above-mentioned method, the step of,for intermediate variable +.>I element of->Is->The first of (3)kThe number of elements to be added to the composition,jin units of imaginary numbers,；

s104, calculating steady-state power amplitude characteristic parameters according to the following，/>Elements of (a) and (b):

，/>，/>，

in the above-mentioned method, the step of,is a steady-state power amplitude characteristic parameter +>Element 1 of->Is a steady-state power amplitude characteristic parameter +>To the first part of (1)l2) +1 elements, ++>Is a steady-state power amplitude characteristic parameter +>At the time of the mth element of (1)mThe range of the values is as follows，/>For intermediate variable +.>Element 1 of->For intermediate variable +.>To the first part of (1)l2) +1 elements; calculating steady-state power frequency characteristic parameters according to the following formula>，/>：

，

In the above-mentioned method, the step of,is a steady-state power frequency characteristic parameter->M element of->Is the sampling rate.

25. The model simulation-based primary frequency modulation parameter adjustment system according to claim 15, wherein when the primary frequency modulation simulation program unit performs primary frequency modulation simulation by combining a pre-established main steam pressure model and a steam turbine model through a control system model to obtain a total number of primary frequency modulation actions and a number of primary frequency modulation qualified times, the primary frequency modulation parameter adjustment system comprises:

s201, acquiring a power grid frequency fluctuation data set，/>Wherein the firstiIndividual element->，/>Wherein->The first of (3)jIndividual element->Represents- >Grid frequency at time,/->Is the grid frequency sampling interval; initializing the total action times of primary frequency modulation and the qualified times of primary frequency modulation;

s202, randomly selecting a power grid frequency fluctuation data setAny elementF _i As the input of the control system model, the control system model is combined with a pre-established main steam pressure model and a steam turbine model to carry out single primary frequency modulation simulation to generate the steam turbineActive power P output by the model;

s203, by elementF _i And calculating the actual integral contribution electric quantity of the simulation according to the active power P output by the steam turbine modelH _i And theoretical integral contribution powerH _e And will actually integrate the contributing powerH _i Divided by theoretical integral contributionH _e Obtaining the primary frequency modulation contribution power rate of the simulationkThe method comprises the steps of carrying out a first treatment on the surface of the Adding 1 to the total action frequency of the primary frequency modulation, if the primary frequency modulation contributes to the electric quantity ratekIf the frequency is larger than or equal to the set value, adding 1 to the qualified frequency of primary frequency modulation;

and S204, judging whether the total frequency modulation action times are greater than or equal to a preset value, if so, judging that the total frequency modulation action times and the qualified frequency modulation times are obtained by carrying out primary frequency modulation simulation through a control system model in combination with a pre-established main steam pressure model and a steam turbine model, otherwise, continuing to iterate in the step S202.

26. The model simulation based primary frequency modulation parameter adjustment system of claim 25, wherein the actual integrated contribution powerH _i The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ the system frequency exceeds the dead zone of primary frequency modulation action of the unit,t _t the moment when the system frequency enters the dead zone of the primary frequency modulation action of the unit,P _t is thattThe time unit actually generates active power,P ₀ is thatt ₀ And the time unit actually generates active power, and t is time.

27. The model simulation based primary frequency modulation parameter adjustment system of claim 25, wherein the theoretical integration contribution power isH _e The expression of the calculation function of (c) is:

in the above-mentioned method, the step of,t ₀ for the moment when the system frequency exceeds the dead zone of the primary frequency modulation action of the unit,t _t the time when the system frequency enters the primary frequency modulation action dead zone of the unit,is thattAt the moment when the power grid frequency exceeds the frequency modulation artificial dead zone value, < ->Rated active output for the unit, < >>Rated frequency for the power grid, < >>The speed change rate of the unit is set, and t is time.

28. The system for adjusting a primary frequency modulation parameter based on model simulation according to claim 15, wherein the step of calculating the primary frequency modulation qualification rate according to the total number of primary frequency modulation actions and the number of primary frequency modulation qualification times in the simulation result judging program unit is to divide the number of primary frequency modulation qualification times by the total number of primary frequency modulation actions to obtain the primary frequency modulation qualification rate.

29. A model simulation based primary frequency modulation parameter adjustment system comprising a microprocessor and a memory connected to each other, wherein the microprocessor is programmed or configured to perform the model simulation based primary frequency modulation parameter adjustment method of any one of claims 1 to 14.

30. A computer readable storage medium having a computer program stored therein, wherein the computer program is for programming or configuring by a microprocessor to perform the model simulation based primary frequency modulation parameter adjustment method of any one of claims 1 to 14.