Feed-forward control method and system for adjusting rotating speed of turbine

Abstract

The invention provides a feedforward control method and system for turbine rotating speed regulation, and the method comprises the steps: collecting a real-time interference parameter of a turbine, and obtaining a historical interference parameter of the same type as the real-time interference parameter and a corresponding historical rotating speed variable quantity through matching; constructing a rotating speed influence coefficient sequence based on the historical interference parameters and the corresponding historical rotating speed variation; performing rotating speed influence analysis based on each real-time interference parameter and the rotating speed influence coefficient sequence, determining a real-time rotating speed influence value, and determining a total interference estimated value based on the real-time rotating speed influence value corresponding to each real-time interference parameter; based on an initial rotating speed difference value of a target rotating speed and a rated rotating speed of the turbine, determining a standard adjusting quantity and a compensation adjusting quantity, and based on the standard adjusting quantity and the compensation adjusting quantity, performing compensation adjustment on the total interference estimated value to obtain a feedforward adjusting quantity; and converting the feedforward adjusting quantity into a driving signal to drive an execution mechanism to act according to the driving signal. According to the invention, the operation stability and anti-interference performance of the turbine are improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of turbine regulation control, and in particular to a feedforward control method and system for turbine speed regulation. BACKGROUND

[0002] In the scenarios involving turbine operation such as ship navigation and power generation, stable regulation of turbine speed is a core link to ensure efficient and safe operation of equipment. The existing turbine speed regulation methods mostly adopt feedback control strategies, that is, by detecting the deviation between the current speed and the target speed in real time, the adjustment amount is calculated based on the deviation value to drive the actuator (such as throttle, guide vane, etc.) to make corrections.

[0003] However, the existing method has obvious regulation lag. When the turbine is running and is disturbed by external interference (such as sudden change of load, fluctuation of fuel quality, sudden change of water flow resistance, etc.), the speed will deviate first, and the feedback system can only start the adjustment action after detecting the deviation, resulting in large speed fluctuation and long time to recover to a stable state. For example, when a ship suddenly encounters strong water flow impact, the load of the propeller will instantaneously increase, and the turbine speed will rapidly decrease. The feedback control needs to go through the process of "detecting deviation-calculating adjustment amount-executing adjustment" before the speed starts to rise. The large fluctuation of the speed in this process not only reduces the power output efficiency, but also may cause additional wear to the turbine mechanical structure, and even cause equipment failure, thereby resulting in poor stability and anti-interference ability of the turbine operation. SUMMARY

[0004] The present application provides a feedforward control method and system for turbine speed regulation to improve the stability and anti-interference ability of the turbine operation.

[0005] In a first aspect, the present application provides a feedforward control method for turbine speed regulation, comprising: collecting real-time interference parameters in the turbine operation process, and matching historical interference parameters of the same type as the real-time interference parameters and corresponding historical speed change amounts in a historical operation database; based on the historical interference parameters of each interference parameter and the corresponding historical speed change amount, constructing a speed influence coefficient sequence; each element in the speed influence coefficient sequence represents a speed influence coefficient caused by each interference parameter under a unit change amount; based on each real-time interference parameter and the speed influence coefficient sequence, performing speed influence analysis to determine the real-time speed influence value corresponding to each real-time interference parameter, and based on the real-time speed influence value corresponding to each real-time interference parameter, determining a total interference estimated value; determining a standard adjustment amount and a compensation adjustment amount based on an initial speed difference between the target speed and the rated speed of the turbine, and compensating the total disturbance estimated value based on the standard adjustment amount and the compensation adjustment amount to obtain a feedforward adjustment amount; converting the feedforward adjustment amount into a driving signal of an actuator to drive the actuator to act according to the driving signal; the actuator includes a throttle regulator and a guide vane angle regulator.

[0006] According to the feedforward control method for turbine speed regulation provided by the application, the determination of the standard adjustment amount and the compensation adjustment amount based on the initial speed difference between the target speed and the rated speed of the turbine includes: determining a dynamic correction coefficient based on the initial speed difference and a speed change rate during turbine operation; the dynamic correction coefficient represents the influence of the change trend of the initial speed difference over time on the adjustment; performing product calculation based on the initial speed difference and the dynamic correction coefficient to obtain a target speed difference after dynamic trend correction; determining a compensation adjustment coefficient based on the target speed difference and the cumulative operation time of the turbine, and determining a basic adjustment amount based on the target speed difference and the speed regulation sensitivity during turbine operation; adjusting the basic adjustment amount based on the compensation adjustment coefficient to obtain the standard adjustment amount and the compensation adjustment amount.

[0007] According to the feedforward control method for turbine speed regulation provided by the application, the adjustment of the basic adjustment amount based on the compensation adjustment coefficient to obtain the standard adjustment amount and the compensation adjustment amount includes: determining an environmental influence coefficient based on the environmental temperature and the environmental air pressure during turbine operation; adjusting the basic adjustment amount based on the environmental influence coefficient to obtain the standard adjustment amount; compensating the standard adjustment amount based on the compensation adjustment coefficient to obtain the compensation adjustment amount.

[0008] According to the feedforward control method for turbine speed regulation provided by the application, the determination of the real-time speed influence value corresponding to each real-time disturbance parameter based on the speed influence coefficient sequence and each real-time disturbance parameter includes: dividing the parameter value of each disturbance parameter into a plurality of continuous and non-overlapping dynamic intervals based on the physical characteristics of each disturbance parameter, and determining the target dynamic interval corresponding to each real-time disturbance parameter; Based on the target dynamic interval and the parameter value of the historical interference parameter corresponding to each speed influence coefficient in the speed influence coefficient sequence, a speed influence coefficient sub-sequence corresponding to the target dynamic interval is matched in the speed influence coefficient sequence; Based on the speed influence coefficient sub-sequence, trend analysis is performed to determine an interval coefficient trend degree; the interval coefficient trend degree represents the overall increasing or decreasing trend of the coefficient in the dynamic interval with the historical interference parameter; Based on the interval coefficient trend degree, speed influence analysis is performed to determine the real-time speed influence value of each real-time interference parameter.

[0009] According to the propeller speed regulation feedforward control method provided by the application, the speed influence analysis based on the interval coefficient trend degree to determine the real-time speed influence value of each real-time interference parameter includes: Based on the speed influence coefficient sub-sequence, fluctuation analysis is performed to determine a fluctuation characteristic value; the fluctuation characteristic value represents the discrete degree of the coefficient in the dynamic interval; Based on the relative position proportion of the parameter value corresponding to each real-time interference parameter in the target dynamic interval, a parameter deviation degree is determined; Based on the interval coefficient trend degree, the fluctuation characteristic value and the parameter deviation degree, a trend correction coefficient of each real-time interference parameter is determined; Based on the trend correction coefficient, speed influence analysis is performed on each real-time interference parameter to determine the real-time speed influence value of each real-time interference parameter.

[0010] According to the propeller speed regulation feedforward control method provided by the application, the historical interference parameter includes a plurality of groups of historical observation values; and the speed influence coefficient sequence is constructed based on the historical interference parameter and the corresponding historical speed change amount of each interference parameter, which includes: For any first interference parameter, based on the difference between adjacent two groups of historical observation values in the historical interference parameter when the first interference parameter changes alone, an observation index sequence is determined; Based on the observation index sequence, a historical observation value sequence of the first interference parameter when the first interference parameter changes alone and a speed change amount sub-sequence corresponding to the historical observation value sequence are extracted; Based on the difference between adjacent two historical observation values in the historical observation value sequence, a parameter change difference value sequence is determined, and based on the difference between adjacent two speed change amounts in the speed change amount sub-sequence, a speed change difference value sequence is determined; Based on the parameter change difference value sequence and the speed change difference value sequence, a speed influence difference value sequence of the first interference parameter under a unit change amount when the first interference parameter changes alone is determined; Based on the speed influence difference value sequence of the first interference parameter, a speed influence coefficient sequence is constructed.

[0011] According to the present application, the method for regulating the rotating speed of the wheel machine by the feedforward control comprises the following steps: The rotating speed influence difference sequence based on the first interference parameter is aggregated in multiple intervals to obtain a first aggregated rotating speed influence value; The first aggregated rotating speed influence value and a second aggregated rotating speed influence value based on a second interference parameter are interactively influenced and fused to obtain an interactive influence stripping coefficient between the first interference parameter and the second interference parameter; the second interference parameter is different from the first interference parameter in type; The rotating speed influence coefficient of the first interference parameter is determined based on the first aggregated rotating speed influence value and the interactive influence stripping coefficient; The rotating speed influence coefficients of the first interference parameter are fused according to the types of the interference parameters to obtain the rotating speed influence coefficient sequence.

[0012] In the second aspect, the present application further provides a feedforward control system for regulating the rotating speed of a wheel machine, which is applied to the method for regulating the rotating speed of the wheel machine by the feedforward control as described in any one of the first aspect, and comprises the following components: The acquisition module is configured to acquire real-time interference parameters in the operation process of the wheel machine, and match historical interference parameters of the same type as the real-time interference parameters and corresponding historical rotating speed change amounts in a historical operation database; The sequence construction module is configured to construct a rotating speed influence coefficient sequence based on the historical interference parameters of each type of interference parameter and the corresponding historical rotating speed change amounts; each element in the rotating speed influence coefficient sequence represents a rotating speed influence coefficient caused by a unit change amount of each type of interference parameter; The interference estimation module is configured to perform rotating speed influence analysis based on each real-time interference parameter and the rotating speed influence coefficient sequence, to determine a real-time rotating speed influence value corresponding to each real-time interference parameter, and to determine a total interference estimated value based on the real-time rotating speed influence value corresponding to each real-time interference parameter; The feedforward regulation amount estimation module is configured to determine a standard regulation amount and a compensation regulation amount based on an initial rotating speed difference value between the target rotating speed and the rated rotating speed of the wheel machine, and to compensate the total interference estimated value based on the standard regulation amount and the compensation regulation amount to obtain a feedforward regulation amount; The feedforward regulation execution module is configured to convert the feedforward regulation amount into a driving signal of an execution mechanism, to drive the execution mechanism to act according to the driving signal; the execution mechanism comprises a throttle regulator and a guide vane angle regulator.

[0013] The application further provides an electronic device, comprising a memory for storing a computer software program; and a processor for reading and executing the computer software program, thereby realizing the feedforward control method for turbine speed regulation.

[0014] The application further provides a non-transitory computer readable storage medium, wherein the storage medium stores a computer software program, and the computer software program is executed by a processor to realize the feedforward control method for turbine speed regulation.

[0015] The application further provides a computer program product, comprising a computer program, and the computer program is executed by a processor to realize the feedforward control method for turbine speed regulation.

[0016] The feedforward control method for turbine speed regulation provided by the application embodiment can obtain factors that may affect the speed in advance before the speed has changed due to the interference, and avoid the problem that the feedback control is passively waiting for the deviation to appear, by collecting real-time interference parameters in the turbine operation process in real time. BRIEF DESCRIPTION OF DRAWINGS

[0017] Figure 1 is a flowchart of the feedforward control method for turbine speed regulation provided by the application; Figure 2 is a structural diagram of the feedforward control system for turbine speed regulation provided by the application; Figure 3 is an embodiment diagram of the electronic device provided by the application embodiment. DETAILED DESCRIPTION

[0018] The technical solutions in the application embodiments will be described clearly and completely below with reference to the drawings in the application embodiments. Obviously, the described embodiments are only part of the embodiments of the application, rather than all the embodiments of the application. Based on the embodiments in the application, all other embodiments obtained by those skilled in the art without creative work fall within the protection scope of the application.

[0019] In the description of the present application, the terms "first", "second" are only for descriptive purposes and cannot be understood to indicate or imply relative importance or imply the number of the indicated technical features. Therefore, the features defined with "first", "second" can explicitly or implicitly include one or more of the features. In the description of the present application, the meaning of "a plurality of" is two or more, unless otherwise specifically limited.

[0020] In the description of the present application, the term "for example" is used to indicate "as an example, illustration or description". Any embodiment described as "for example" in the present application is not necessarily interpreted as more preferred or more advantageous than other embodiments. The following description is given in order to enable any person skilled in the art to implement and use the present application. In the following description, details are listed for the purpose of explanation. It should be understood that those skilled in the art can recognize that the present application can be implemented without using these specific details. In other examples, well-known structures and processes will not be described in detail to avoid unnecessary details making the description of the present application obscure. Therefore, the present application is not intended to be limited to the embodiments shown, but is consistent with the broadest scope consistent with the principles and characteristics disclosed.

[0021] Optionally, referring to Figure 1 , Figure 1 is a flowchart of the feedforward control method for turbine speed regulation provided by the present application. The execution subject of the feedforward control method for turbine speed regulation in the embodiments of the present application is a feedforward control system. Therefore, the feedforward control method for turbine speed regulation comprises: Step 10, collecting real-time interference parameters in the running process of the turbine, and matching the historical interference parameters of the same type as the real-time interference parameters and the corresponding historical speed change amount in the historical running database.

[0022] Optionally, the feedforward control system collects real-time interference parameters in the running process of the turbine. The real-time interference parameters include the current load value, the real-time fuel pressure, the intake flow, the environmental temperature and the ship draft. During the collection process, the feedforward control system acquires data through sensors installed at key parts of the turbine, ensuring the real-time and accuracy of the data. The type of the sensor needs to match the type of the parameter, for example, the current load value can be collected by a load sensor, the real-time fuel pressure is collected by a pressure sensor, etc.

[0023] Further, after the collection is completed, the feedforward control system matches the real-time interference parameters with the historical interference parameters of the same type stored in the historical running database and the corresponding historical speed change amount. The historical running database stores a large amount of running data of the turbine under different working conditions. When matching, the principle of consistent parameter type needs to be followed, that is, the current load value corresponds to the historical load value, the real-time fuel pressure corresponds to the historical fuel pressure, etc.

[0024] In an embodiment, the real-time interference parameters collected by the feedforward control system at a certain time are: current load value 800 kW, real-time fuel pressure 5 MPa, intake flow rate 200 m³ / h, ambient temperature 30 °C, and ship draft 5 m. The feedforward control system starts the matching program and searches the historical running database for historical interference parameters of the same type as these parameters and corresponding historical speed variation amounts. After searching, a set of matching data is found: the historical speed variation amount corresponding to the historical load value 750 kW is -10 r / min; the historical speed variation amount corresponding to the historical fuel pressure 4.8 MPa is -5 r / min; the historical speed variation amount corresponding to the historical intake flow rate 190 m³ / h is +8 r / min; the historical speed variation amount corresponding to the historical ambient temperature 28 °C is +3 r / min; and the historical speed variation amount corresponding to the historical ship draft 4.8 m is -2 r / min. This set of data is the historical data that matches the real-time interference parameters successfully.

[0025] Step 20, based on the historical interference parameters of each interference parameter and the corresponding historical speed variation amount, a speed influence coefficient sequence is constructed. Each element in the speed influence coefficient sequence represents the speed influence coefficient caused by a unit variation of each interference parameter.

[0026] Further, the feedforward control system constructs a speed influence coefficient sequence based on the historical interference parameters of each interference parameter and the corresponding historical speed variation amount obtained in step 10. For each interference parameter, the speed influence coefficient caused by a unit variation of the interference parameter needs to be calculated, and the speed influence coefficient is an element in the speed influence coefficient sequence, specifically as in the processes of steps 201 to 205.

[0027] Step 30, based on each real-time interference parameter and the speed influence coefficient sequence, a speed influence analysis is performed to determine the real-time speed influence value corresponding to each real-time interference parameter, and based on the real-time speed influence value corresponding to each real-time interference parameter, a total interference estimated value is determined.

[0028] Further, the feedforward control system performs a speed influence analysis based on each real-time interference parameter and the corresponding speed influence coefficient sequence to determine the real-time speed influence value corresponding to each real-time interference parameter, specifically as in the processes of steps 301 to 304.

[0029] Further, the feedforward control system adds the real-time speed influence values corresponding to each real-time interference parameter, and the sum is the total interference estimated value. The total interference estimated value reflects the predicted influence amount of the current all real-time interference parameters on the speed of the turbine.

[0030] Step 40, based on the initial speed difference between the target speed and the rated speed of the turbine, determine the standard adjustment amount and the compensation adjustment amount, and compensate the total disturbance estimate based on the standard adjustment amount and the compensation adjustment amount to obtain the feedforward adjustment amount.

[0031] Further, the feedforward control system calculates the initial speed difference between the target speed and the rated speed of the turbine, i.e. initial speed difference = target speed - rated speed, and determines the standard adjustment amount and the compensation adjustment amount based on the initial speed difference. The standard adjustment amount is the basic adjustment amount calculated only according to the initial speed difference, and its calculation formula can be set as standard adjustment amount = initial speed difference * standard adjustment coefficient, and the standard adjustment coefficient is determined according to the characteristics of the turbine. The compensation adjustment amount is an additional adjustment amount set considering the adjustment deviation caused by various interference factors, which can be set according to the size and direction of the initial speed difference, for example, when the initial speed difference is larger, the compensation adjustment amount takes a larger value, and the specific process is shown in steps 401 to 404.

[0032] Further, the feedforward control system adds the standard adjustment amount and the compensation adjustment amount to obtain the total adjustment compensation amount, and then compensates the total disturbance estimate with the total adjustment compensation amount, i.e. feedforward adjustment amount = total adjustment compensation amount - total disturbance estimate, to obtain the final feedforward adjustment amount.

[0033] Step 50, convert the feedforward adjustment amount into a driving signal of the actuator based on the feedforward adjustment amount to drive the actuator to act according to the driving signal. The actuator includes a throttle regulator and a guide vane angle regulator.

[0034] Further, the feedforward control system needs to convert the feedforward adjustment amount into a driving signal of the actuator. During the conversion process, the conversion relationship needs to be determined according to the type and characteristics of the actuator, and different actuators correspond to different conversion parameters. For the throttle regulator, its driving signal is usually a current signal or a voltage signal, and the conversion relationship between the feedforward adjustment amount and the driving signal can be calibrated through experiments, for example, 1 r / min of feedforward adjustment amount corresponds to 0.02 A of current change. For the guide vane angle regulator, the driving signal may be an angle control signal, and the conversion relationship is the corresponding relationship between the feedforward adjustment amount and the angle change, for example, 1 r / min of feedforward adjustment amount corresponds to 0.5° of angle change. After the conversion is completed, the feedforward control system sends the driving signal to the corresponding actuator to drive the actuator to act according to the driving signal, thereby realizing the adjustment of the turbine speed.

[0035] In an embodiment, the known feedforward adjustment amount is 60.5 r / min. For the throttle regulator, setting a feedforward adjustment amount of 1 r / min corresponds to a current signal change of 0.02 A, so its drive signal = 60.5 r / min * 0.02 A / (r / min) = 1.21 A. For the guide vane angle regulator, setting a feedforward adjustment amount of 1 r / min corresponds to an angle signal change of 0.5°, so its drive signal = 60.5 r / min * 0.5° / (r / min) = 30.25°. The feedforward control system sends a current signal of 1.21 A to the throttle regulator and an angle signal of 30.25° to the guide vane angle regulator to drive the two to take corresponding actions to adjust the turbine speed to the target value.

[0036] The embodiment of the present application can obtain factors that may affect the speed in advance before the speed has changed due to the interference by collecting real-time interference parameters in the turbine operation process in real time, and avoids the problem of passive waiting for deviation of feedback control. The influence of the interference parameters on the speed is calculated based on the historical data, the total estimated value of the real-time interference caused by the real-time interference parameters is estimated in advance, and the influence of the interference is predicted. The total estimated value of the interference is compensated by generating a compensation adjustment amount containing a target speed adjustment amount and a compensation adjustment amount required to offset the interference based on the speed of the turbine, and a feedforward adjustment amount is obtained, so that the adjustment action starts before the speed deviates. Finally, the actuator is driven to act in advance, so that the adjustment amount starts to work while the interference factors act on the turbine, directly suppresses the fluctuation of the speed, solves the problem of hysteresis, and improves the stability and anti-interference ability of the turbine operation.

[0037] In an embodiment, the process of steps 201 to 205 includes: Step 201, for any first interference parameter, determining an observation index sequence based on the difference between adjacent two groups of historical observation values in the historical interference parameters when the first interference parameter changes alone.

[0038] Optionally, the first interference parameter in the embodiment of the present application is any one of the current load value, real-time fuel pressure, intake flow, ambient temperature and ship draft, and the feedforward control system selects the historical interference parameters when the first interference parameter changes alone from the historical operation database for the first interference parameter, the historical interference parameters include multiple groups of historical observation values, and then calculates the difference between adjacent two groups of historical observation values, wherein the historical observation value refers to the specific value of the first interference parameter recorded in the historical operation database. Further, the feedforward control system determines an observation index sequence capable of representing the change order of the first interference parameter according to the difference between adjacent two groups of historical observation values, and the elements in the index sequence are serial number identifiers of adjacent two groups of historical observation values.

[0039] In an embodiment, the first interference parameter is selected as the current load value, and the historical interference parameter data is screened from the historical operation database, in which the current load value varies alone and the real-time fuel pressure, the intake flow, the ambient temperature and the ship draft depth remain stable. Assuming that the screened historical observation values are 600 kW (serial number 1), 650 kW (serial number 2), 700 kW (serial number 3), 750 kW (serial number 4) and 800 kW (serial number 5) in turn. The difference between adjacent two groups of historical observation values is calculated, 650 kW-600 kW = 50 kW, 700 kW-650 kW = 50 kW, 750 kW-700 kW = 50 kW and 800 kW-750 kW = 50 kW, and these differences are all effective changes. The observation index sequence determined therefrom is [1, 2, 3, 4, 5], which reflects the historical observation order when the current load value varies alone.

[0040] In step 202, the historical observation value sequence of the first interference parameter when varying alone is extracted based on the observation index sequence, and the rotational speed change quantum sequence corresponding to the historical observation value sequence is extracted.

[0041] Further, the feedforward control system extracts the historical observation value sequence of the first interference parameter when varying alone and the corresponding rotational speed change quantum sequence based on the observation index sequence. In the extraction process, the historical observation value of the first interference parameter corresponding to the serial number in the observation index sequence is called from the historical operation database to form the historical observation value sequence. At the same time, the historical rotational speed change corresponding to these historical observation values is called to form the rotational speed change quantum sequence.

[0042] Optionally, the embodiment of the present application needs to ensure that the elements in the historical observation value sequence and the rotational speed change quantum sequence correspond to each other one by one, that is, the rotational speed change corresponding to the nth element in the historical observation value sequence is the nth element in the rotational speed change quantum sequence, so as to ensure the accuracy of subsequent analysis.

[0043] Continuing the above embodiment, the observation index sequence is [1, 2, 3, 4, 5], the feedforward control system extracts the historical observation value of the current load value from the historical operation database to obtain the historical observation value sequence [600 kW, 650 kW, 700 kW, 750 kW, 800 kW]. At the same time, the historical rotational speed change corresponding to the serial number is extracted, which is assumed to be 10 r / min, 5 r / min, 0 r / min, -5 r / min and -10 r / min respectively, so that the rotational speed change quantum sequence is [10 r / min, 5 r / min, 0 r / min, -5 r / min, -10 r / min]. The elements in these two sequences correspond to each other one by one, such as the historical observation value 600 kW corresponding to the rotational speed change 10 r / min, 650 kW corresponding to 5 r / min, and so on.

[0044] Step 203, based on the difference between the adjacent two historical observation values in the historical observation value sequence, determine the parameter change difference value sequence, and based on the difference between the adjacent two speed change quantities in the speed change quantity sequence, determine the speed change difference value sequence.

[0045] Further, the feedforward control system calculates the difference between adjacent elements for the obtained historical observation value sequence and speed change quantity sequence, respectively, to obtain the parameter change difference value sequence and the speed change difference value sequence. For the historical observation value sequence, the nth element in the parameter change difference value sequence is equal to the (n+1) th element minus the nth element in the historical observation value sequence (n starts from 1). For the speed change quantity sequence, the nth element in the speed change difference value sequence is equal to the (n+1) th element minus the nth element in the speed change quantity sequence (n starts from 1). By calculating the difference between adjacent elements, the change amplitude of the first disturbance parameter and the corresponding speed change amplitude can be clearly reflected.

[0046] Continuing the above example, the historical observation value sequence of the current load value is [600kW, 650kW, 700kW, 750kW, 800kW], then the parameter change difference value sequence is calculated: 650kW-600kW=50kW, 700kW-650kW=50kW, 750kW-700kW=50kW, 800kW-750kW=50kW, the parameter change difference value sequence is [50kW, 50kW, 50kW, 50kW]. The speed change quantity sequence is [10r / min, 5r / min, 0r / min, -5r / min, -10r / min], the speed change difference value sequence is calculated: 5r / min-10r / min=-5r / min, 0r / min-5r / min=-5r / min, -5r / min-0r / min=-5r / min, -10r / min-(-5r / min)=-5r / min, the speed change difference value sequence is [-5r / min, -5r / min, -5r / min, -5r / min].

[0047] Step 204, based on the parameter change difference value sequence and the speed change difference value sequence, determine the speed influence difference value sequence of the first disturbance parameter under the unit change quantity when the first disturbance parameter changes alone.

[0048] Further, the feedforward control system calculates a rotation speed influence difference sequence of the first interference parameter under a unit change based on the parameter change difference sequence and the rotation speed change difference sequence. The calculation method is that the nth element in the rotation speed influence difference sequence is equal to the nth element in the rotation speed change difference sequence divided by the nth element in the parameter change difference sequence. The element represents the rotation speed change difference caused by the unit change of the first interference parameter in the change process of the adjacent two groups of historical observation values, and can reflect the change of the rotation speed influence of the first interference parameter in different change stages.

[0049] Continuing the above embodiment, the parameter change difference sequence of the current load value is [50kW, 50kW, 50kW, 50kW], and the rotation speed change difference sequence is [-5r / min, -5r / min, -5r / min, -5r / min]. Then, the elements in the rotation speed influence difference sequence are calculated as follows: -5r / min÷50kW=-0.1r / (min·kW), -5r / min÷50kW=-0.1r / (min·kW), -5r / min÷50kW=-0.1r / (min·kW), -5r / min÷50kW=-0.1r / (min·kW), that is, the rotation speed influence difference sequence is [-0.1r / (min·kW), -0.1r / (min·kW), -0.1r / (min·kW), -0.1r / (min·kW)].

[0050] In step 205, a rotation speed influence coefficient sequence is constructed based on the rotation speed influence difference sequence of the first interference parameter.

[0051] Further, the feedforward control system constructs the rotation speed influence coefficient sequence according to the rotation speed influence difference sequence of the first interference parameter, specifically as the processes of steps 2051 to 2054.

[0052] The embodiment of the present application constructs the rotation speed influence coefficient sequence of each interference parameter, and the rotation speed influence coefficient sequence accurately represents the rotation speed influence coefficient caused by a unit change of each interference parameter, thereby providing reliable rotation speed influence coefficients of each interference parameter for determination of the total interference estimated value and final adjustment control, and effectively improving the accuracy and foresight of turbine rotation speed adjustment.

[0053] In an embodiment, the processes of steps 2051 to 2054 include: In step 2051, a first aggregated rotation speed influence value is obtained by multi-interval aggregation based on the rotation speed influence difference sequence and the parameter change difference sequence of the first interference parameter.

[0054] Optionally, the feedforward control system aggregates the speed influence difference sequence and the parameter change difference sequence for the first interference parameter in multiple intervals to obtain a first aggregated speed influence value, specifically: The feedforward control system divides the parameter change difference sequence into multiple intervals, and the division basis can be the size range of the parameter change difference, for example, elements with similar parameter change differences are classified into the same interval. For each interval, elements in the corresponding speed influence difference sequence are extracted, and then the average value of the speed influence differences in the interval is calculated (the weight can be set according to the size of the parameter change difference) to obtain the aggregated speed influence value of the interval. Finally, the aggregated speed influence values of all intervals are aggregated again, and the arithmetic average method is adopted in the embodiment of the application to obtain the first aggregated speed influence value.

[0055] Continue with the above embodiment, taking the first interference parameter as an example, the current load value is, the speed influence difference sequence is [-0.1r / (min·kW), -0.1r / (min·kW), -0.1r / (min·kW), -0.1r / (min·kW)], and the parameter change difference sequence is [50kW, 50kW, 50kW, 50kW]. The parameter change difference sequence is divided into one interval (because all elements are 50kW), and the speed influence difference sequence elements corresponding to the interval are all -0.1r / (min·kW). The aggregated speed influence value of the interval is (-0.1-0.1-0.1-0.1)÷4=-0.1r / (min·kW). Since there is only one interval, the first aggregated speed influence value is -0.1r / (min·kW).

[0056] In another embodiment, the parameter change difference sequence of the current load value is [40kW, 60kW, 40kW, 60kW], and the corresponding speed influence difference sequence is [-0.12r / (min·kW), -0.08r / (min·kW), -0.11r / (min·kW), -0.09r / (min·kW)]. It is divided into two intervals: interval 1 is the parameter change difference 40kW, and the corresponding speed influence difference is [-0.12, -0.11]; interval 2 is the parameter change difference 60kW, and the corresponding speed influence difference is [-0.08, -0.09]. The aggregated value of interval 1 is (-0.12-0.11)÷2=-0.115r / (min·kW), and the aggregated value of interval 2 is (-0.08-0.09)÷2=-0.085r / (min·kW). The first aggregated speed influence value is (-0.115-0.085)÷2=-0.1r / (min·kW).

[0057] Step 2052, based on the first aggregate speed influence value and the second aggregate speed influence value of the second interference parameter, the interaction influence fusion is carried out to obtain the interaction influence stripping coefficient between the first interference parameter and the second interference parameter. The second interference parameter is a different type of interference parameter from the first interference parameter.

[0058] Further, the feedforward control system carries out interaction influence fusion based on the first aggregate speed influence value and the second aggregate speed influence value of the second interference parameter to obtain the interaction influence stripping coefficient between the first interference parameter and the second interference parameter. The second interference parameter is a different type of interference parameter from the first interference parameter, for example, when the first interference parameter is the current load value, the second interference parameter can be the real-time fuel pressure.

[0059] Optionally, the actual interaction influence fusion of the application is realized by analyzing the historical data when the two interference parameters change together, specifically: the feedforward control system extracts the historical observation value and the corresponding speed change amount when the two interference parameters change together from the historical operation database, and calculates the actual total speed influence value under the joint action of the two; then according to the first aggregate speed influence value and the second aggregate speed influence value, the theoretical total speed influence value when the two act alone is calculated; finally, the interaction influence amount is determined by the difference between the actual total speed influence value and the theoretical total speed influence value, and the interaction influence stripping coefficient is obtained, and the calculation formula is: interaction influence stripping coefficient = interaction influence amount ÷ (first aggregate speed influence value * second aggregate speed influence value), and the interaction influence stripping coefficient is used to represent the degree of interaction between the two interference parameters on the speed influence.

[0060] Continue the above embodiment, the first interference parameter is the current load value, and the first aggregate speed influence value is -0.1 r / (min·kW); the second interference parameter is the real-time fuel pressure, and the second aggregate speed influence value is -0.05 r / (min·MPa). A group of data is extracted from the historical operation database when the two change together: the current load value changes by 50 kW, the real-time fuel pressure changes by 10 MPa, and the corresponding actual total speed influence value is -5.5 r / min. The theoretical total speed influence value = (50 kW * -0.1 r / (min·kW)) + (10 MPa * -0.05 r / (min·MPa)) = -5 r / min - 0.5 r / min = -5.5 r / min. The difference between the actual total speed influence value and the theoretical total speed influence value (interaction influence amount) is -5.5-(-5.5)=0 r / min. Then the interaction influence stripping coefficient = 0 ÷ (-0.1 * -0.05) = 0. In another embodiment, the actual total speed influence value is -6 r / min, the theoretical total speed influence value is -5.5 r / min, and the interaction influence amount is -0.5 r / min. The interaction influence stripping coefficient = (-0.5) ÷ (-0.1 * -0.05) = (-0.5) ÷ 0.005 = -100.

[0061] In step 2053, the first interference parameter is determined based on the first aggregated speed influence value and the interaction influence stripping coefficient.

[0062] Further, the feedforward control system determines the speed influence coefficient caused by the unit change of the first interference parameter based on the first aggregated speed influence value and the interaction influence stripping coefficient. The calculation method of the embodiment of the present application is: speed influence coefficient = first aggregated speed influence value - (interaction influence stripping coefficient * first aggregated speed influence value * second aggregated speed influence value ÷ unit change reference value of the second interference parameter). Wherein, the unit change reference value of the second interference parameter is the unit change commonly seen in historical data (such as the unit change reference value of the real-time fuel pressure being 1 MPa).

[0063] The embodiment of the present application strips the interaction between the two interference parameters to obtain the speed influence coefficient under the unit change when the first interference parameter acts alone, ensuring the independence and accuracy of the coefficient.

[0064] Continuing the above embodiment, the first aggregated speed influence value is -0.1 r / (min·kW), the interaction influence stripping coefficient is 0, the second aggregated speed influence value is -0.05 r / (min·MPa), and the unit change reference value of the second interference parameter (real-time fuel pressure) is 1 MPa. Then the speed influence coefficient of the current load value is -0.1 - (0 * -0.1 * -0.05 ÷ 1) = -0.1 r / (min·kW). If the interaction influence stripping coefficient is -100 and other parameters remain unchanged: the speed influence coefficient is -0.1 - (-100 * -0.1 * -0.05 ÷ 1) = -0.05 r / (min·kW).

[0065] In step 2054, the speed influence coefficients of the first interference parameter are fused according to the types of the interference parameters to obtain a speed influence coefficient sequence.

[0066] Further, the feedforward control system fuses the speed influence coefficients of the first interference parameter according to the types of the interference parameters to obtain a speed influence coefficient sequence. The fusion process of the embodiment of the present application collects multiple speed influence coefficients calculated through different historical data sets for the same type of first interference parameter, and then processes these coefficients using statistical methods (such as arithmetic mean, median, etc.) to obtain the final speed influence coefficient sequence. The length of the speed influence coefficient sequence is determined according to the number of historical data sets, and each element represents the unit change speed influence coefficient of the interference parameter in different historical scenarios, ensuring that the sequence can fully reflect the influence characteristics of the interference parameter.

[0067] Continuing the above embodiment, the first interference parameter is the current load value, and the rotation speed influence coefficients calculated by multiple sets of historical data are -0.1 r / (min·kW), -0.11 r / (min·kW), -0.09 r / (min·kW), and -0.1 r / (min·kW) respectively. When fusion is performed by type, the arithmetic mean is used to obtain the rotation speed influence coefficient sequence as [-0.1, -0.11, -0.09, -0.1] r / (min·kW) (the original data is retained as the sequence elements, or the average value is taken as a single element sequence, and here the original data is taken as an example).

[0068] The rotation speed influence coefficient sequence constructed in the embodiment of the application not only comprehensively aggregates the influences of multiple intervals, but also strips the interaction influences with other interference parameters, ensures that each element can truly reflect the influence of the first interference parameter on the rotation speed when the first interference parameter acts alone under a unit change, provides a reliable data basis for subsequent calculation of a real-time rotation speed influence value based on a real-time interference parameter and determination of a total interference estimated value, improves the precision of turbine rotation speed adjustment, and can more effectively cope with various interference factors in advance to ensure the stability of turbine operation.

[0069] In an embodiment, the process of steps 301 to 304 includes: Step 301, based on the physical characteristics of each interference parameter, the parameter values are divided into multiple continuous and non-overlapping dynamic intervals, and the target dynamic interval corresponding to each real-time interference parameter is determined.

[0070] Optionally, the feedforward control system divides the parameter values of each interference parameter into dynamic intervals according to the physical characteristics of the interference parameter. The physical characteristics of different interference parameters are different, and the division basis is also different. For example, the physical characteristics of the current load value are related to the carrying capacity of the turbine, and the intervals can be divided in proportion according to the maximum load and the minimum load of the turbine; the physical characteristics of the real-time fuel pressure are related to the working pressure range of the fuel delivery system, and the intervals can be divided according to the measurement range of the pressure sensor.

[0071] Optionally, the multiple dynamic intervals divided in the embodiment of the application need to meet the requirements of continuity and non-overlapping, to ensure that each parameter value can correspond to a unique dynamic interval.

[0072] Further, after the division is completed, the feedforward control system compares each real-time interference parameter value collected in real time with each dynamic interval to determine the target dynamic interval corresponding to the real-time interference parameter.

[0073] Continuing with the current load value as an example, its physical characteristics determine that the parameter value range is 0kW to 2000kW. The feedforward control system divides it into 5 continuous and non-overlapping dynamic intervals, namely [0kW, 400kW), [400kW, 800kW), [800kW, 1200kW), [1200kW, 1600kW), [1600kW, 2000kW]. If the current load value collected in real time is 700kW, it can be known through comparison that the value falls within the [400kW, 800kW) interval, so the target dynamic interval corresponding to the current load value is [400kW, 800kW). Taking the ambient temperature as an example, its physical characteristics determine that the parameter value range is -10°C to 50°C, which is divided into [-10°C, 10°C), [10°C, 30°C), [30°C, 50°C] three dynamic intervals, and the real-time ambient temperature is 25°C. The corresponding target dynamic interval is [10°C, 30°C).

[0074] Step 302, based on the target dynamic interval, the parameter value of the historical interference parameter corresponding to each speed influence coefficient in the speed influence coefficient sequence is matched to the speed influence coefficient sub-sequence corresponding to the target dynamic interval in the speed influence coefficient sequence.

[0075] Further, the feedforward control system extracts the parameter value of the historical interference parameter corresponding to each speed influence coefficient in the speed influence coefficient sequence. The parameter value of the historical interference parameter is the historical data relied on when constructing the speed influence coefficient sequence.

[0076] Further, the feedforward control system compares the parameter value of the historical interference parameter with the target dynamic interval, and selects the speed influence coefficient corresponding to the parameter value of the historical interference parameter falling within the target dynamic interval. The sequence of coefficients is the speed influence coefficient sub-sequence corresponding to the target dynamic interval.

[0077] Continuing with the above example, the target dynamic interval of the current load value is [400kW, 800kW), and its speed influence coefficient sequence is [-0.1r / (min·kW), -0.11r / (min·kW), -0.09r / (min·kW), -0.1r / (min·kW)], and the historical load value corresponding to each coefficient is 300kW, 500kW, 600kW, 900kW respectively. Comparing the historical load value with the target dynamic interval [400kW, 800kW), 500kW and 600kW fall within the interval, and the corresponding speed influence coefficients are -0.11r / (min·kW) and -0.09r / (min·kW), so the matched speed influence coefficient sub-sequence is [-0.11r / (min·kW), -0.09r / (min·kW)].

[0078] For the ambient temperature, the target dynamic interval is [10℃, 30℃), the rotational speed influence coefficient sequence is [0.02r / (min·℃), 0.03r / (min·℃), 0.025r / (min·℃), 0.035r / (min·℃)], and the corresponding historical ambient temperatures are 5℃, 15℃, 25℃ and 35℃ respectively. The matched rotational speed influence coefficient sub-sequence is [0.03r / (min·℃), 0.025r / (min·℃)].

[0079] In step 303, trend analysis is performed based on the rotational speed influence coefficient sub-sequence to determine the interval coefficient trend degree. The interval coefficient trend degree represents the overall increasing or decreasing trend of the coefficient in the dynamic interval with the change of the historical interference parameter.

[0080] Further, the feedforward control system performs trend analysis based on the rotational speed influence coefficient sub-sequence to determine the interval coefficient trend degree. The trend analysis can be achieved by calculating the difference between adjacent rotational speed influence coefficients in the sub-sequence. If most of the differences are positive, it indicates that the coefficient as a whole is in an upward trend. If most of the differences are negative, it indicates that the overall trend is downward. If the differences alternate between positive and negative and have similar absolute values, it indicates that the overall trend is stable. The interval coefficient trend degree is represented by a specific numerical value or descriptive words, such as "upward", "downward" and "stable", or by a trend slope. A positive slope indicates an upward trend, a negative slope indicates a downward trend, and a slope close to zero indicates a stable trend.

[0081] Continuing with the above example, the rotational speed influence coefficient sub-sequence of the current load value is [-0.11r / (min·kW), -0.09r / (min·kW)], and the difference between adjacent coefficients is (-0.09)-(-0.11)=0.02r / (min·kW). The difference is positive, indicating that the interval coefficient trend degree is upward. The rotational speed influence coefficient sub-sequence of the ambient temperature is [0.03r / (min·℃), 0.025r / (min·℃)], and the difference between adjacent coefficients is 0.025-0.03=-0.005r / (min·℃). The difference is negative, indicating that the interval coefficient trend degree is downward.

[0082] In step 304, rotational speed influence analysis is performed based on the interval coefficient trend degree to determine the real-time rotational speed influence value of each real-time interference parameter.

[0083] Further, the feedforward control system performs rotational speed influence analysis based on the interval coefficient trend degree to determine the real-time rotational speed influence value of each real-time interference parameter, as described in steps 3041 to 3044.

[0084] The embodiment of the present application divides the dynamic interval according to the physical characteristics, ensures the pertinence of parameter analysis, matches the coefficient subsequence of the target interval, improves the relevance of the coefficient and the real-time parameter, considers the change rule of the coefficient in the interval in the trend analysis, makes the calculated real-time speed influence value more in line with the actual situation, provides reliable basic data for subsequent determination of the total estimated value of the interference, and thus can more accurately predict the comprehensive influence of various interference factors on the speed, and improves the accuracy and stability of the turbine speed control.

[0085] In an embodiment, the process of steps 3041 to 3044 includes: Step 3041, performing fluctuation analysis based on the speed influence coefficient subsequence to determine a fluctuation characteristic value. The fluctuation characteristic value represents the discrete degree of the coefficient in the dynamic interval.

[0086] Optionally, the feedforward control system performs fluctuation analysis on the speed influence coefficient subsequence to determine a fluctuation characteristic value. The fluctuation characteristic value is used to quantify the discrete degree of the coefficient in the dynamic interval, and reflects the stability of the coefficient in the interval. Optionally, in the analysis process of the embodiment of the present application, firstly, the difference sequence of adjacent coefficients in the subsequence is calculated, and then statistical analysis is performed on the difference sequence. The combination index of the mean square error (MSE) and the range (R) of the difference sequence is used to represent the fluctuation characteristic value. The mean square error reflects the average deviation degree of the coefficient change, and the range reflects the maximum fluctuation range.

[0087] Continuing the above embodiment, the speed influence coefficient subsequence of the current load value is [-0.11, -0.09 r / (min·kW)], and the difference of adjacent coefficients is 0.02 r / (min·kW). The mean square error of the difference sequence is 0.02²=0.0004, and the range is 0.02. The fluctuation characteristic value calculation formula is: fluctuation characteristic value , is the mean value of the subsequence, that is, (-0.11-0.09) / 2=-0.1 r / (min·kW). Substituting the calculation gives: fluctuation characteristic=0.024.

[0088] Step 3042, determining the parameter deviation degree based on the relative position proportion of the parameter value corresponding to each real-time interference parameter in the target dynamic interval.

[0089] Further, the feedforward control system calculates the relative position proportion of the parameter value corresponding to each real-time interference parameter in the target dynamic interval to determine the parameter deviation degree.

[0090] Optionally, the parameter deviation degree in the embodiment of the present application measures the degree of deviation of the real-time parameter from the center value of the interval, which is obtained by calculating the proportion of the difference between the real-time parameter value and the lower limit value of the interval to the total width of the interval. The proportion value is subjected to nonlinear transformation processing through a preset transformation calculation formula to enhance the sensitivity to the edge area of the interval. The parameter deviation degree calculation formula in the embodiment of the present application is: parameter deviation degree=(real-time parameter value-lower limit value of interval) / (upper limit value of interval-lower limit value of interval), and the preset transformation calculation formula is: nonlinear parameter deviation degree= .

[0091] In an embodiment, the real-time parameter value of the current load value is 700kW, the target dynamic interval is [400kW, 800kW), and the interval width is 400kW. Therefore, the parameter deviation degree=0.75, and the nonlinear parameter deviation degree= ≈0.4621.

[0092] Step 3043, based on the interval coefficient trend degree, the fluctuation characteristic value and the parameter deviation degree, the trend correction coefficient of each real-time interference parameter is determined.

[0093] Further, the feedforward control system calculates the trend correction coefficient based on the interval coefficient trend degree, the fluctuation characteristic value and the parameter deviation degree. The trend correction coefficient is used to adjust the basic speed influence coefficient to more accurately reflect the speed influence under the real-time working condition. In the embodiment of the present application, the interval coefficient trend degree is first converted into a numerical representation (rising as 1, stable as 0, and falling as-1), and then the correction coefficient is calculated through a composite function which considers the interaction among the three. The specific formula of the composite function is: trend correction coefficient= .

[0094] Continuing the above embodiment, the interval coefficient trend degree of the current load value is rising (numerical representation is 1), the fluctuation characteristic value is 0.024, and the nonlinear parameter deviation degree is 0.4621. Therefore, the trend correction coefficient= ≈1.570.

[0095] Step 3044, based on the trend correction coefficient, the speed influence analysis is performed on each real-time interference parameter to determine the real-time speed influence value of each real-time interference parameter.

[0096] Further, the feedforward control system corrects the basic speed influence coefficient based on the trend correction coefficient to determine the real-time speed influence value. In the embodiment of the present application, the basic coefficient is obtained by geometric weighted average of the subsequence coefficients, and the weight is determined by the parameter deviation degree. Then the basic coefficient is multiplied by the trend correction coefficient to obtain the final real-time speed influence coefficient, and then multiplied by the real-time parameter change amount to obtain the real-time speed influence value.

[0097] Optionally, the calculation formula of the base coefficient is: base coefficient = (current load value - base load value) / base load value. .

[0098] Continuing the above example, the current load value's speed influence coefficient subsequence is [-0.11 r / (min·kW), -0.09 r / (min·kW)], and the parameter shift degree is 0.75. The base coefficient = -0.1061 r / (min·kW). ≈-0.1061 r / (min·kW). The real-time speed influence coefficient = base coefficient * trend correction coefficient ≈-0.1666 r / (min·kW), the real-time parameter change amount is 700 kW-500 kW = 200 kW, and the real-time speed influence value = real-time speed influence coefficient * parameter change amount ≈-33.32 r / min.

[0099] The embodiment of the present application captures the discrete characteristics of the coefficient through fluctuation analysis, quantifies the relative position of the real-time parameter in the interval through the parameter shift degree, and comprehensively considers the nonlinear interaction of the three through the trend correction coefficient, so that the influence of the interference parameter on the speed under different working conditions can be accurately reflected, especially when the parameter fluctuation is large or close to the interval boundary, the prediction accuracy is improved, and therefore the real-time speed influence value obtained is more in line with the actual operation, providing a reliable basis for subsequent interference total estimation value calculation and feedforward regulation, effectively improving the precision and stability of the turbine speed control.

[0100] In an embodiment, the process of steps 401 to 404 includes: Step 401, based on the initial speed difference value, the speed change rate during turbine operation is determined to determine the dynamic correction coefficient. The dynamic correction coefficient represents the influence of the change trend of the initial speed difference value over time on the adjustment.

[0101] Optionally, the feedforward control system determines the dynamic correction coefficient based on the initial speed difference value and the speed change rate during turbine operation. The initial speed difference value is the difference between the target speed of the turbine and the rated speed, and the speed change rate refers to the change amount of the speed per unit time, which is calculated by continuously collecting speed data through a speed sensor. The dynamic correction coefficient in the embodiment of the present application needs to consider the positive and negative of the initial speed difference value and the trend of the speed change rate: when the speed change rate and the initial speed difference value are in the same direction (such as the initial speed difference value is positive and the speed is still rising), the dynamic correction coefficient is greater than 1, to enhance the adjustment strength; when they are in opposite directions (such as the initial speed difference value is positive and the speed is falling), the dynamic correction coefficient is less than 1, to weaken the adjustment strength; when the speed change rate is zero, the dynamic correction coefficient is 1.

[0102] In an embodiment, the target speed of the turbine is 1500 r / min, the rated speed is 1450 r / min, and the initial speed difference value = 1500 r / min-1450 r / min = 50 r / min (positive value).

[0103] The rotation speed collected by the rotation speed sensor rises from 1440 r / min to 1450 r / min in the past 10 seconds, the rotation speed change rate = (1450 r / min-1440 r / min) ÷ 10 s = 1 r / (min・s) (positive value, same direction with the initial rotation speed difference). The calculation formula of the dynamic correction coefficient is: dynamic correction coefficient = 1 + (rotation speed change rate ÷ 10), and the data is substituted to get: dynamic correction coefficient = 1 + (1 ÷ 10) = 1.1. If the rotation speed change rate is -0.5 r / (min・s) (negative value, opposite direction with the initial rotation speed difference), then the dynamic correction coefficient = 1 + (-0.5 ÷ 10) = 0.95.

[0104] Step 402, based on the initial rotation speed difference and the dynamic correction coefficient, the product is calculated to get the target rotation speed difference after dynamic trend correction.

[0105] Further, the feedforward control system based on the initial rotation speed difference and the dynamic correction coefficient, the product is calculated to get the target rotation speed difference after dynamic trend correction, target rotation speed difference = initial rotation speed difference * dynamic correction coefficient, so that the corrected target rotation speed difference is more in line with the current dynamic operation state of the turbine.

[0106] Continue the above embodiment, the initial rotation speed difference is 50 r / min, and the dynamic correction coefficient is 1.1, so that the target rotation speed difference after dynamic trend correction = 50 r / min * 1.1 = 55 r / min. If the dynamic correction coefficient is 0.95, then the target rotation speed difference = 50 r / min * 0.95 = 47.5 r / min.

[0107] Step 403, based on the target rotation speed difference and the cumulative running time of the turbine, the compensation adjustment coefficient is determined, and based on the target rotation speed difference and the rotation speed adjustment sensitivity of the turbine running, the basic adjustment amount is determined.

[0108] Further, the feedforward control system based on the target rotation speed difference and the cumulative running time of the turbine, the compensation adjustment coefficient is determined. The cumulative running time refers to the total running time of the turbine from start to the current time, which is recorded by the timer. The compensation adjustment coefficient needs to consider the influence of the cumulative running time on the performance of the turbine: the longer the cumulative running time, the more serious the wear of the turbine components, and the larger the compensation adjustment coefficient to offset the adjustment deviation caused by performance degradation. At the same time, based on the target rotation speed difference and the rotation speed adjustment sensitivity of the turbine running, the basic adjustment amount is determined. The rotation speed adjustment sensitivity refers to the rotation speed change caused by unit adjustment amount, which is fitted by historical adjustment data, and the basic adjustment amount = target rotation speed difference ÷ rotation speed adjustment sensitivity.

[0109] Continuing the above example, the target speed difference value is 55 r / min, and the cumulative running time of the turbine is 1000 hours. The calculation formula for setting the compensation adjustment coefficient is: compensation adjustment coefficient = 1 + (cumulative running time ÷ 10000), and the data is substituted to obtain: compensation adjustment coefficient = 1 + (1000 ÷ 10000) = 1.1. If the cumulative running time is 2000 hours, then the compensation adjustment coefficient = 1 + (2000 ÷ 10000) = 1.2.

[0110] The speed regulation sensitivity of the turbine during operation is known to be 0.5 r / (min·unit adjustment amount) (i.e., for every 1 unit adjustment amount, the speed increases by 0.5 r / min), so the basic adjustment amount = 55 r / min ÷ 0.5 r / (min·unit adjustment amount) = 110 unit adjustment amounts.

[0111] Step 404, adjusting the basic adjustment amount based on the compensation adjustment coefficient to obtain a standard adjustment amount and a compensation adjustment amount.

[0112] Further, the feedforward control system adjusts the basic adjustment amount based on the compensation adjustment coefficient to obtain a standard adjustment amount and a compensation adjustment amount, specifically as the process of steps 4041 to 4043.

[0113] The embodiments of the present application comprehensively consider the dynamic change trend of the initial speed difference value (through the dynamic correction coefficient), the influence of the cumulative running time of the turbine (through the compensation adjustment coefficient), and the speed regulation sensitivity, so that the obtained adjustment amount not only meets the current speed regulation basic demand, but also adapts to the dynamic operation state of the turbine and the performance change after long-term use, providing a reliable data basis for subsequent feedforward adjustment amount and improving the adaptability and reliability of the turbine speed regulation.

[0114] In an embodiment, the process of steps 4041 to 4043 includes: Step 4041, determining an environmental influence coefficient based on the environmental temperature and the environmental air pressure during operation of the turbine.

[0115] Optionally, the feedforward control system determines the environmental influence coefficient based on the environmental temperature and the environmental air pressure during operation of the turbine. The environmental temperature and the air pressure are collected in real time by sensors, and both have a significant influence on the operating efficiency of the turbine. The system first calculates a temperature influence factor and an air pressure influence factor, respectively, and then fuses the two through a nonlinear function. The temperature influence factor is based on the principle of thermodynamics and considers the change of air density with temperature; the air pressure influence factor reflects the influence of atmospheric pressure on the intake amount. The calculation formula of the environmental influence coefficient is: Environmental influence coefficient = .

[0116] Wherein, is the real-time environmental temperature, T0 is the standard ambient temperature (usually taken as 298K), T1 is the temperature influence coefficient (about -0.0035 / K); P0 is the real-time ambient pressure, P1 is the standard atmospheric pressure (101.3kPa), P2 is the pressure influence coefficient (about 0.002 / kPa).

[0117] In an embodiment, the real-time ambient temperature T = 308K, and the ambient pressure P = 100kPa, thus the ambient influence coefficient is calculated as .

[0118] Step 4042, adjust the base adjustment amount based on the ambient influence coefficient to obtain a standard adjustment amount.

[0119] Further, the feedforward control system adjusts the base adjustment amount based on the ambient influence coefficient to obtain a standard adjustment amount. The actual operating efficiency of the turbine is changed by environmental factors, so the base adjustment amount needs to be corrected. The standard adjustment amount calculation formula is: standard adjustment amount = base adjustment amount / ambient influence coefficient, which ensures that the adjustment amount is increased in harsh environmental conditions (such as high temperature and low pressure), and the adjustment amount is reduced in favorable environmental conditions, maintaining the stability of the turbine output.

[0120] Continuing the above embodiment, the base adjustment amount is 110 units of adjustment amount, and the ambient influence coefficient is 0.9987, so the standard adjustment amount = 110 / 0.9987 ≈ 110.15 units of adjustment amount.

[0121] Step 4043, compensate the standard adjustment amount based on the compensation adjustment coefficient to obtain a compensation adjustment amount.

[0122] Further, the feedforward control system compensates the standard adjustment amount based on the compensation adjustment coefficient to obtain a compensation adjustment amount. The compensation adjustment coefficient reflects the influence of the cumulative running time of the turbine on the performance. As the running time increases, the wear and tear of the turbine components leads to a decrease in efficiency, and additional adjustment amount is needed to compensate. The compensation adjustment amount calculation formula is: compensation adjustment amount = standard adjustment amount*(compensation adjustment coefficient-1), which ensures that the compensation amount is proportional to the standard adjustment amount and increases with the increase of the running time.

[0123] Continuing the above embodiment, the compensation adjustment coefficient is 1.1, and the standard adjustment amount is 110.15 units of adjustment amount, so the compensation adjustment amount = 110.15*(1.1-1) ≈ 11.02 units of adjustment amount.

[0124] The embodiment of the present application can adapt the adjustment amount to different environmental conditions through the environmental influence coefficient, ensure the stable operation of the turbine under various working conditions, compensate the adjustment amount for the dynamic correction of the long-term performance attenuation of the turbine, prolong the effective service life of the turbine, and finally improve the precision and reliability of the turbine speed control through the combined action of the standard adjustment amount and the compensation adjustment amount, especially when the environmental conditions fluctuate greatly or the turbine is operated for a long time, the target speed can be maintained more effectively, and the fluctuation and error are reduced.

[0125] Further, the turbine speed regulation feedforward control system provided by the present application is described below, and the turbine speed regulation feedforward control system described below can be referred to in correspondence with the turbine speed regulation feedforward control method described above.

[0126] Optionally, referring to Figure 2 , Figure 2 is a structural schematic diagram of the turbine speed regulation feedforward control system provided by the present application, and the turbine speed regulation feedforward control system comprises: The acquisition module 210 acquires real-time interference parameters in the turbine operation process, and matches historical interference parameters of the same type as the real-time interference parameters and corresponding historical speed change amounts in the historical operation database; The sequence construction module 220 is configured to construct a speed influence coefficient sequence based on the historical interference parameters of each interference parameter and the corresponding historical speed change amount; each element in the speed influence coefficient sequence represents a speed influence coefficient caused by a unit change amount of each interference parameter; The interference estimation module 230 is configured to perform speed influence analysis based on each real-time interference parameter and the speed influence coefficient sequence, determine a real-time speed influence value corresponding to each real-time interference parameter, and determine a total estimated value of the interference based on the real-time speed influence value corresponding to each real-time interference parameter; The feedforward adjustment amount estimation module 240 is configured to determine a standard adjustment amount and a compensation adjustment amount based on an initial speed difference between the target speed of the turbine and the rated speed, and compensate the total estimated value of the interference based on the standard adjustment amount and the compensation adjustment amount, to obtain a feedforward adjustment amount; The feedforward adjustment execution module 250 is configured to convert the feedforward adjustment amount into a driving signal of an execution mechanism, to drive the execution mechanism to act according to the driving signal; the execution mechanism comprises a throttle regulator and a guide vane angle regulator.

[0127] This invention, through real-time acquisition of disturbance parameters during turbine operation, can identify potential factors affecting speed before the speed changes due to disturbances, avoiding the problem of passive waiting for deviations in feedback control. By calculating the impact of disturbance parameters on speed using historical data, the total estimated value of real-time disturbances caused by these parameters is predicted, enabling anticipation of disturbance impacts. Based on the turbine's speed, a standard adjustment amount is generated to reach the target speed, along with a compensation adjustment amount required to offset the disturbance. This compensation adjustment is applied to the total estimated disturbance value, resulting in a feedforward adjustment amount. This ensures that the adjustment action is initiated before speed deviation occurs. Finally, the actuator is driven to act in advance, so that the adjustment amount begins to take effect simultaneously with the disturbance, directly suppressing speed fluctuations, solving the lag problem, and improving the stability and anti-interference capability of turbine operation.

[0128] Please see Figure 3 , Figure 3 This is a schematic diagram illustrating an embodiment of the electronic device provided in this invention. For example... Figure 3 As shown, this embodiment of the invention provides an electronic device 300, including a memory 310, a processor 320, and a computer program 311 stored in the memory 310 and executable on the processor 320. When the processor 320 executes the computer program 311, it performs the following steps: Collect real-time disturbance parameters during turbine operation and match them with historical disturbance parameters of the same type as the real-time disturbance parameters and their corresponding historical speed changes in the historical operation database; Based on the historical interference parameters and corresponding historical speed changes for each interference parameter, a speed influence coefficient sequence is constructed; each element in the speed influence coefficient sequence represents the speed influence coefficient caused by each interference parameter under a unit change. Based on the sequence of each real-time interference parameter and the speed influence coefficient, the speed influence is analyzed to determine the real-time speed influence value corresponding to each real-time interference parameter, and based on the real-time speed influence value corresponding to each real-time interference parameter, the total interference estimate is determined. Based on the initial speed difference between the target speed and the rated speed of the turbine, the standard adjustment and the compensation adjustment are determined, and the total estimated disturbance value is compensated and adjusted based on the standard adjustment and the compensation adjustment to obtain the feedforward adjustment. The feedforward adjustment is converted into a drive signal for the actuator to drive the actuator to perform actions according to the drive signal; the actuator includes a throttle regulator and a guide vane angle regulator.

[0129] It should be noted that the descriptions of each embodiment in the above embodiments have different focuses. For parts that are not described in detail in a certain embodiment, please refer to the relevant descriptions in other embodiments.

[0130] Those skilled in the art will appreciate that embodiments of the present application can be devised for a variety of applications. It is intended that the present application be limited only by the scope of the appended claims, and it is intended that various modifications and alterations made by those skilled in the art be considered as within the scope of the present application. The embodiments of the present application will be described with reference to the attached drawings identified below.

[0131] The present application is described in reference to the drawings using a flowchart and / or a block diagram of the method, apparatus (system) and computer program product according to embodiments of the application. It will be understood that each block of the flowchart and / or block diagram, and combinations of blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded computer, 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 and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks.

[0132] These computer program instructions can 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 instructions which implement the function specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks.

[0133] The computer program instructions can also be loaded into 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 and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks. Figure 1 one or more functions specified in the flowchart and / or block diagram block or blocks.

[0134] While the preferred embodiments of the application have been described, additional variations and modifications can be made to the embodiments by those skilled in the art once they learn of the basic inventive concepts. Therefore, the appended claims are intended to cover all such modifications and variations as fall within the scope of the present application.

[0135] Obviously, many modifications and variations of the present application are possible in light of the above teachings. It is, therefore, to be understood that within the scope of the appended claims and their equivalents, the application can be practiced otherwise than as specifically described.

Claims

1. A feedforward control method for turbine speed regulation, characterized in that, include; Collect real-time disturbance parameters during turbine operation and match them with historical disturbance parameters of the same type as the real-time disturbance parameters and their corresponding historical speed changes in the historical operation database; Based on the historical interference parameters and corresponding historical speed changes for each interference parameter, a speed influence coefficient sequence is constructed; each element in the speed influence coefficient sequence represents the speed influence coefficient caused by each interference parameter under a unit change. Based on each real-time interference parameter and the speed influence coefficient sequence, a speed influence analysis is performed to determine the real-time speed influence value corresponding to each real-time interference parameter, and based on the real-time speed influence value corresponding to each real-time interference parameter, the total interference estimate is determined. Based on the initial speed difference between the target speed and the rated speed of the turbine, a standard adjustment amount and a compensation adjustment amount are determined, and the total estimated value of the disturbance is compensated and adjusted based on the standard adjustment amount and the compensation adjustment amount to obtain the feedforward adjustment amount; The feedforward adjustment amount is converted into a drive signal for the actuator to drive the actuator to perform actions according to the drive signal; the actuator includes a throttle regulator and a guide vane angle regulator.

2. The feedforward control method for turbine speed regulation according to claim 1, characterized in that, The determination of the standard adjustment and compensation adjustment based on the initial speed difference between the target speed and the rated speed of the turbine includes: Based on the initial speed difference and the rate of change of turbine speed during operation, a dynamic correction coefficient is determined; the dynamic correction coefficient characterizes the influence of the change trend of the initial speed difference over time on the regulation. The target speed difference after dynamic trend correction is obtained by multiplying the initial speed difference and the dynamic correction coefficient. Based on the target speed difference and the turbine's cumulative operating time, a compensation adjustment coefficient is determined, and based on the target speed difference and the turbine's speed adjustment sensitivity during operation, a basic adjustment amount is determined. The basic adjustment amount is adjusted based on the compensation adjustment coefficient to obtain the standard adjustment amount and the compensation adjustment amount.

3. The feedforward control method for turbine speed regulation according to claim 2, characterized in that, The step of adjusting the basic adjustment amount based on the compensation adjustment coefficient to obtain the standard adjustment amount and the compensation adjustment amount includes: The environmental impact coefficient is determined based on the ambient temperature and air pressure during turbine operation. The standard adjustment amount is obtained by adjusting the basic adjustment amount based on the environmental impact coefficient. The standard adjustment amount is compensated based on the compensation adjustment coefficient to obtain the compensation adjustment amount.

4. The feedforward control method for turbine speed regulation according to claim 1, characterized in that, The step of performing speed influence analysis based on each real-time interference parameter and the speed influence coefficient sequence to determine the real-time speed influence value corresponding to each real-time interference parameter includes: Based on the physical characteristics of each interference parameter, its parameter value is divided into multiple continuous and non-overlapping dynamic intervals, and the target dynamic interval corresponding to each real-time interference parameter is determined. Based on the target dynamic range and the parameter values ​​of the historical interference parameters corresponding to each speed influence coefficient in the speed influence coefficient sequence, a subsequence of speed influence coefficients corresponding to the target dynamic range is matched in the speed influence coefficient sequence. Trend analysis is performed on the subsequence of the rotational speed influence coefficient to determine the interval coefficient trend degree; the interval coefficient trend degree characterizes the overall increase or decrease trend of the coefficient within the dynamic interval as the historical disturbance parameter changes. Based on the trend of the interval coefficient, the rotational speed influence analysis is performed to determine the real-time rotational speed influence value of each real-time disturbance parameter.

5. The feedforward control method for turbine speed regulation according to claim 4, characterized in that, The step of performing speed influence analysis based on the trend degree of the interval coefficient to determine the real-time speed influence value of each real-time disturbance parameter includes: Fluctuation analysis is performed based on the subsequence of the rotational speed influence coefficient to determine the fluctuation characteristic value; the fluctuation characteristic value characterizes the degree of dispersion of the system within the dynamic range; The degree of parameter offset is determined based on the relative position ratio of the parameter value corresponding to each real-time interference parameter within the target dynamic range. Based on the trend degree of the interval coefficient, the fluctuation characteristic value, and the parameter offset degree, a trend correction coefficient for each real-time disturbance parameter is determined; Based on the trend correction coefficient, the rotational speed influence of each real-time disturbance parameter is analyzed to determine the real-time rotational speed influence value of each real-time disturbance parameter.

6. The feedforward control method for turbine speed regulation according to any one of claims 1 to 5, characterized in that, The historical interference parameters include multiple sets of historical observations; the sequence of speed influence coefficients, constructed based on the historical interference parameters and corresponding historical speed changes for each interference parameter, includes: For any first disturbance parameter, the observation index sequence is determined based on the difference between two adjacent sets of historical observations in the historical disturbance parameter when it changes individually; Based on the observation index sequence, extract the historical observation value sequence when the first disturbance parameter changes individually, and the corresponding rotation speed change quantum sequence; Based on the difference between two adjacent historical observations in the historical observation sequence, a parameter change difference sequence is determined, and based on the difference between two adjacent speed changes in the speed change quantum sequence, a speed change difference sequence is determined. Based on the parameter change difference sequence and the rotational speed change difference sequence, determine the rotational speed influence difference sequence when the first disturbance parameter changes alone under a unit change amount; A sequence of speed influence coefficients is constructed based on the speed influence difference sequence of the first interference parameter.

7. The feedforward control method for turbine speed regulation according to claim 6, characterized in that, The step of constructing a speed influence coefficient sequence based on the speed influence difference sequence of the first interference parameter includes: Based on the speed influence difference sequence and parameter change difference sequence of the first interference parameter, multi-interval aggregation is performed to obtain the first aggregated speed influence value; Based on the first polymerization speed influence value and the second interference parameter's second polymerization speed influence value, an interactive influence fusion is performed to obtain the interactive influence stripping coefficient between the first interference parameter and the second interference parameter; the second interference parameter is an interference parameter of a different type from the first interference parameter. Based on the first aggregation speed influence value and the interaction influence stripping coefficient, determine the speed influence coefficient caused by the first interference parameter under a unit change. The rotational speed influence coefficients of the first interference parameter are fused according to the type of interference parameter to obtain the rotational speed influence coefficient sequence.

8. A feedforward control system for turbine speed regulation, characterized in that, The feedforward control method for turbine speed regulation as described in any one of claims 1 to 7, wherein the feedforward control system for turbine speed regulation comprises: The acquisition module is used to collect real-time disturbance parameters during the operation of the turbine and match them with historical disturbance parameters of the same type as the real-time disturbance parameters and the corresponding historical speed changes in the historical operation database. The sequence construction module is used to construct a speed influence coefficient sequence based on the historical interference parameters and corresponding historical speed changes of each interference parameter; each element in the speed influence coefficient sequence represents the speed influence coefficient caused by each interference parameter under a unit change. The interference prediction module is used to perform speed influence analysis based on each real-time interference parameter and the speed influence coefficient sequence, determine the real-time speed influence value corresponding to each real-time interference parameter, and determine the total interference prediction value based on the real-time speed influence value corresponding to each real-time interference parameter. The feedforward regulation estimation module is used to determine the standard regulation and the compensation regulation based on the initial speed difference between the target speed and the rated speed of the turbine, and to compensate and adjust the total disturbance estimate based on the standard regulation and the compensation regulation to obtain the feedforward regulation. The feedforward adjustment execution module is used to convert the feedforward adjustment amount into a drive signal for the actuator, so as to drive the actuator to perform actions according to the drive signal; the actuator includes a throttle adjuster and a guide vane angle adjuster.

9. An electronic device, characterized in that, The electronic device includes a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor, when executing the program, implements the feedforward control method for turbine speed regulation as described in any one of claims 1 to 7.

10. A non-transitory computer-readable storage medium, wherein a computer software program is stored therein, characterized in that, When the processor executes the program, it implements the feedforward control method for turbine speed regulation as described in any one of claims 1 to 7.

Images