A control method for a steam supply device in a thermal power production system

Abstract

The application relates to a control method of steam supply equipment in a thermoelectric production system, which comprises the following steps: collecting a set of main steam header pressures and a set of steam supply header pressures corresponding to a first specified moment t of an i-th group of steam supply equipment groups according to time intervals; obtaining a main steam header pressure trend and a steam supply header pressure trend of the moment t by adopting a pre-set strategy for the set of main steam header pressures and the set of steam supply header pressures respectively; obtaining a main steam header pressure prediction sequence based on the main steam header pressure corresponding to the moment t, the main steam header pressure trend of the moment t and a first recursive formula, and obtaining a steam supply header pressure prediction sequence based on the steam supply header pressure corresponding to the moment t, the steam supply header pressure trend of the moment t and a second recursive formula; and obtaining optimal adjustment values of the main steam header and the steam supply header based on the main steam header pressure prediction sequence, the steam supply header pressure prediction sequence and a control solution space corresponding to the i-th group of steam supply equipment groups.

Description

Technical Field

[0001] This invention relates to the field of energy production technology, and in particular to a control method for steam supply equipment in a thermal power production system. Background Technology

[0002] Fluctuations in the steam pressure of the main pipe in a thermal power production system not only affect the safe operation of the unit, but also cause instability in the external heat supply parameters. Due to the strong heat storage capacity of the steam main pipe, the main pipe pressure has a large time delay characteristic. At the same time, the delay of boiler combustion and the strong correlation between different boilers bring greater challenges to the control of the unit's main pipe pressure.

[0003] Existing control methods typically employ a combination of fixed-pressure regulating boilers, pressure regulating furnaces, and non-pressure regulating furnaces. In this method, the non-pressure regulating furnace controls the boiler outlet steam flow rate to maintain a constant boiler load. However, the boiler outlet steam flow rate fluctuates with changes in the main steam header pressure. Furthermore, controlling the main steam header pressure solely through the boiler cannot solve the problem of large boiler response delays. When external load disturbances occur, the main steam header pressure parameters are prone to significant short-term fluctuations, leading to a decline in the unit's economic operating performance and affecting its safe operation. Summary of the Invention

[0004] (a) Technical problems to be solved

[0005] In view of the above-mentioned shortcomings and deficiencies of the prior art, the present invention provides a control method for steam supply equipment in a thermal power production system, which solves the technical problems of insufficient control accuracy and control lag in the pressure control of the main pipe of the existing steam supply equipment.

[0006] (II) Technical Solution

[0007] To achieve the above objectives, the main technical solutions adopted by the present invention include:

[0008] A control method for steam supply equipment in a thermal power production system includes:

[0009] A1. Collect the set S of main steam header pressures corresponding to the i-th group of steam supply equipment in the thermal power production system and the first specified time t according to a preset time interval. ai-t The set of pressures S of the steam supply header bi-t The set S of main steam header pressures corresponding to the first specified time t ai-t This includes the main steam header pressure collected sequentially within a first specified time period before the first specified time t; and the set S of the steam supply header pressure corresponding to the first specified time t. bi-t This includes the steam supply main pipe pressure collected sequentially within the first specified time period before the first specified time t;

[0010] The first specified time period is [t-(T+1)Δ, t]; Δ is a pre-set time interval;

[0011] A2, for the main steam header pressure corresponding to the first specified time t, the set S ai-t and the set S of the main steam header pressure corresponding to the first specified time t bi-t respectively using a pre-set strategy, respectively obtaining the main steam header pressure trend trend a·t and the main steam header pressure trend trend b·t ;

[0012] A3, based on the main steam header pressure corresponding to the time t, the main steam header pressure trend trend a·t and the first recursive formula, obtaining the main steam header pressure prediction sequence S ai , and based on the main steam header pressure corresponding to the time t, the main steam header pressure trend trend b·t and the second recursive formula, obtaining the main steam header pressure prediction sequence S bi ;

[0013] A4, based on the main steam header pressure prediction sequence S ai , the main steam header pressure prediction sequence S bi and the control solution space C corresponding to the i-th group of steam supply equipment group in the thermal power production system obtained in advance, obtaining the optimal adjustment value of the main steam header and the optimal adjustment value of the steam supply header.

[0014] Preferably,

[0015] The thermal power production system includes N groups of steam supply equipment groups;

[0016] Any steam supply equipment group includes a first steam supply equipment or includes a plurality of second steam supply equipment meeting the first communication condition;

[0017] Wherein, the first communication condition is that all the main steam headers connected in front of the second steam supply equipment are in communication with each other, and all the steam supply headers connected behind the second steam supply equipment are in communication with each other;

[0018] The first steam supply equipment is a steam supply equipment meeting the second communication condition; the second communication condition is that the main steam header connected in front of the steam supply equipment is not in communication with the main steam header of other steam supply equipment, and the steam supply header connected behind the steam supply equipment is not in communication with the steam supply header of other steam supply equipment;

[0019] Wherein, the set S of the main steam header pressure corresponding to the first specified time t ai-t is:

[0020] Sai-t=[pa t-T+1 ,... pat-c ,...pa t-2 pa t-1 pa t ];

[0021] T is any positive integer between [3, 10];

[0022] pa t-c The main steam header pressure is defined as the pressure of the main steam header at c time intervals prior to the first specified time t.

[0023] The set S of steam supply header pressures corresponding to the first specified time t bi-t for:

[0024] Sbi-t = [pb] t-T+1 ...pb t-c ...pb t-2 ,pb t-1 [pbt];

[0025] pb t-c The steam supply header pressure is the pressure at c time intervals prior to the first specified time t.

[0026] Preferably,

[0027] In A2, the set S of main steam header pressures corresponding to the first specified time t ai-t Using a pre-set strategy, the main steam header pressure trend at time t is obtained. a·t Specifically, it includes:

[0028] Using a fixed index sequence [1, 2, 3...T] as the independent variable, and the set S of the main steam header pressure corresponding to the first specified time t... ai-t As the dependent variable, a pre-defined fitting method was used to obtain the main steam header pressure trend at time t. a·t Specifically, this includes: using a T-degree polynomial to fit the independent and dependent variables to obtain the first polynomial function corresponding to the dependent variable;

[0029] Based on the first polynomial function corresponding to the dependent variable, its first derivative Pa′(t) is obtained, and the value of this first derivative Pa′(t) at point T is taken as the main steam header pressure trend at time t. a·t ;

[0030] In A2, the set S of steam supply header pressures corresponding to the first specified time t bi-t Using a pre-set strategy, the steam supply header pressure trend at time t is obtained. b·t Specifically, it includes:

[0031] The fixed index sequence [1, 2, 3...T] is taken as an independent variable, and a set of the steam supply main pipe pressure corresponding to the first specified moment t is taken as a dependent variable bi-t The steam supply main pipe pressure trend trend at the moment t is obtained by using a preset fitting mode as a dependent variable b·t , specifically comprising: fitting the independent variable and the dependent variable by using a T-order polynomial to obtain a second polynomial function corresponding to the dependent variable;

[0032] According to the second polynomial function corresponding to the dependent variable, the first derivative Pb'(t) is obtained, and the value of the first derivative Pb'(t) at T is taken as the steam supply main pipe pressure trend trend at the moment t b.t .

[0033] Preferably, the A3 specifically comprises:

[0034] A3-1, based on the main steam supply main pipe pressure corresponding to the moment t, the main steam supply main pipe pressure trend trend at the moment t a.t and the first recursive formula, the main steam supply main pipe pressure prediction sequence S is obtained ai , specifically comprising:

[0035] A3-1-1, the main steam supply main pipe pressure corresponding to the moment t, the main steam supply main pipe pressure trend trend at the moment t a.t and the pre-obtained flow change amount dsteam in the main steam supply main pipe at the moment t a.t are substituted into the first recursive formula to obtain the main steam supply main pipe pressure pa t+1 corresponding to the moment t+Δ;

[0036] pa t+1 = pa t +trend a·t +f*dsteam a·t ;

[0037] The first recursive formula is:

[0038] pa t+k = pa t+k-1 +trend a·t+k-1 +f*dsteam a·t+k-1 ;

[0039] A3-1-2, based on the set S of the main steam supply main pipe pressure at the moment t ai-t and the main steam supply main pipe pressure pa t+1 corresponding to the moment t+Δ, the first recursive formula is used to obtain the main steam supply main pipe pressure of the second time interval to the k time interval after the moment t, respectively, and the main steam supply main pipe pressure of the k time interval after the moment t is taken as the main steam supply main pipe pressure prediction sequence S ai ;

[0040] wherein, trend a·t+k-1 is a set of main steam header pressure corresponding to k-1 time intervals after the first specified time t with fixed index sequence [1, 2, 3…T] as independent variable, and the collection S ai-t-k-1 is dependent variable, which is obtained by using a pre-set fitting method;

[0041] k is greater than or equal to 2, pa t+k-1 is obtained by using the first recursive formula;

[0042] f is a model coefficient obtained in advance in the first recursive formula;

[0043] dsteam a·t+k-1 = steam a·t+k-1 - steam a·t+k-2 ;

[0044] dsteam a·t+k-1 is the flow change amount in the main steam header in the k-1 time interval after the time t;

[0045] steam a·t+k-1 is the flow in the main steam header in the k-1 time interval after the time t;

[0046] The main steam header pressure prediction sequence S ai is:

[0047] S ai = [pa t+1 , pa t+2 , pa t+3 , …, pa t+k ];

[0048] A3-2, based on the time t corresponding to the steam header pressure, the steam header pressure trend trend b-t and the second recursive formula, the steam header pressure prediction sequence S bi is obtained, specifically including:

[0049] A3-2-1, the time t corresponding to the steam header pressure, the steam header pressure trend trend b·t and the flow change amount dsteam b·t in the steam header at time t are substituted into the second recursive formula to obtain t +Δtime corresponding to the steam header pressure pb t+1 ;

[0050] wherein, pb t+1 = pb t +trend b·t+ g * dsteam b·t ;

[0051] The second recursive formula is:

[0052] pb t+k = pb t+k-1 + trend b·t+k-1 + g * dsteam b·t+k-1 .

[0053] A3-2-2, a set S of the main steam supply pipe pressure based on time t bi-t and the main steam supply pipe pressure pb corresponding to time t + Δ t+1 , a second recursive formula is used to obtain the main steam supply pipe pressure after the second time interval to the kth time interval after the current time t, respectively, and the main steam supply pipe pressure after the k time intervals after time t forms a main steam supply pipe pressure prediction sequence S bi ;

[0054] Wherein, trend b·t+k-1 is obtained by using a predetermined fitting method with a fixed index sequence [1, 2, 3... T] as the independent variable and a set S of the main steam supply pipe pressure corresponding to the k-1 time intervals after the first specified time t bi-t-k-1 as the dependent variable;

[0055] pb t+k-1 is obtained by using the second recursive formula;

[0056] g is a model coefficient obtained in advance in the second recursive formula;

[0057] dsteam b·t+k-1 = steam b·t+k-1 - steam b.t+k-2 ;

[0058] dsteam b·t+k-1 is the flow change amount in the main steam supply pipe for the k-1 time intervals after time t;

[0059] steam b·t+k-1 is the flow in the main steam supply pipe for the k-1 time intervals after time t;

[0060] The main steam supply pipe pressure prediction sequence S bi is:

[0061] S bi = [pb t+1 , pb t+2 , pb t+3 , …, pb t+k ].

[0062] Preferably,

[0063] The model coefficient f in the first recursive formula is obtained in the following manner:

[0064] From the beginning to the end of the second specified time period, a time window with a length of L x Δ is translated sequentially H times from the beginning to the end of the second specified time period, and each time the translation is Δ;

[0065] The second specified time period is [t0-(H+L)Δ, t0];

[0066] Wherein, t0 is the second specified time; H is a pre-set number; L is a random integer between [6, 12];

[0067] After each translation of the time window, a set of main steam header steam flow at the time corresponding to the beginning of the time window and a set of main steam header pressure at the time corresponding to the beginning of the time window are obtained, and a set of main steam header steam flow at the time corresponding to the end of the time window and a set of main steam header pressure at the time corresponding to the end of the time window are obtained;

[0068] For the main steam header steam flow at the time corresponding to the beginning of the time window and the main steam header steam flow at the time corresponding to the end of the time window obtained after the dth translation of the time window, the difference value of the main steam header steam flow corresponding to the dth translation of the time window is obtained; wherein, 1≤d≤H;

[0069] Taking the fixed index sequence [1, 2, 3...T] as the independent variable and the set of main steam header pressure at the time corresponding to the beginning of the time window obtained after the dth translation of the time window as the dependent variable, a pre-set fitting method is used to obtain the main steam header pressure trend at the time corresponding to the beginning of the time window obtained after the dth translation of the time window;

[0070] Taking the fixed index sequence [1, 2, 3...T] as the independent variable and the set of main steam header pressure at the time corresponding to the end of the time window obtained after the dth translation of the time window as the dependent variable, a pre-set fitting method is used to obtain the main steam header pressure trend at the time corresponding to the end of the time window obtained after the dth translation of the time window;

[0071] Based on the main steam header pressure trend at the time corresponding to the beginning of the time window obtained after the dth translation of the time window and the main steam header pressure trend at the time corresponding to the end of the time window obtained after the dth translation of the time window, the difference value of the main steam header pressure trend corresponding to the dth translation of the time window is obtained;

[0072] According to the difference of the main steam header steam flow corresponding to each time window translation and the difference of the main steam header pressure trend, a first two-dimensional array is generated;

[0073] The first two-dimensional array is:

[0074]

[0075] delta_steam ad is the difference of the main steam header steam flow corresponding to the dth time window translation; delta_trend ad is the difference of the main steam header pressure trend corresponding to the dth time window translation.

[0076] The difference of the main steam header steam flow corresponding to each translation in the first two-dimensional array and the difference of the main steam header pressure trend are respectively substituted into the first non-intercept linear regression model obtained in advance to train the specific value of the regression model coefficient f0 in the first non-intercept linear regression model, and the specific value of the regression model coefficient f0 in the first non-intercept linear regression model is used as the model coefficient f in the first recursive formula;

[0077] The first non-intercept linear regression model is:

[0078] The difference of the main steam header pressure trend=f0×the difference of the main steam header steam flow.

[0079] Preferably,

[0080] The model coefficient g in the second recursive formula is obtained in the following manner:

[0081] From the beginning to the end of the second specified time period, a time window with a length of L×Δ is used to translate H times from the beginning to the end of the second specified time period, and each translation is by Δ.

[0082] After each translation of the time window, the set of the main steam header steam flow at the time corresponding to the beginning of the time window and the set of the main steam header pressure at the time corresponding to the beginning of the time window, and the set of the main steam header steam flow at the time corresponding to the end of the time window and the set of the main steam header pressure at the time corresponding to the end of the time window are obtained.

[0083] For the dth time window translation, the difference of the main steam header steam flow corresponding to the dth time window translation is obtained from the main steam header steam flow at the time corresponding to the beginning of the time window and the main steam header steam flow at the time corresponding to the end of the time window.

[0084] Taking the fixed index sequence [1, 2, 3...T] as the independent variable, taking the set of the steam supply main pressure corresponding to the time point at the start of the time window obtained after the dth translation of the time window as the dependent variable, and adopting a preset fitting mode, the steam supply main pressure trend corresponding to the time point at the start of the time window obtained after the dth translation of the time window is obtained;

[0085] Taking the fixed index sequence [1, 2, 3...T] as the independent variable, taking the set of the steam supply main pressure corresponding to the time point at the end of the time window obtained after the dth translation of the time window as the dependent variable, and adopting a preset fitting mode, the steam supply main pressure trend corresponding to the time point at the end of the time window obtained after the dth translation of the time window is obtained;

[0086] Based on the steam supply main pressure trend corresponding to the time point at the start of the time window obtained after the dth translation of the time window and the steam supply main pressure trend corresponding to the time point at the end of the time window obtained after the dth translation of the time window, the difference value of the steam supply main pressure trend corresponding to the dth translation of the time window is obtained;

[0087] According to the difference value of the steam supply main steam flow corresponding to each translation of the time window and the difference value of the steam supply main pressure trend, a second two-dimensional array is generated;

[0088] The second two-dimensional array is:

[0089]

[0090] delta_steam bd is the difference value of the steam supply main steam flow corresponding to the dth translation of the time window;

[0091] delta_trend bd is the difference value of the steam supply main pressure trend corresponding to the dth translation of the time window;

[0092] The difference value of the steam supply main steam flow corresponding to each translation in the second two-dimensional array and the difference value of the steam supply main pressure trend are respectively substituted into the second non-intercept linear regression model obtained in advance, and the specific value of the regression model coefficient g0 in the second non-intercept linear regression model is trained to obtain the model coefficient g in the second recursive formula;

[0093] The second non-intercept linear regression model is:

[0094] The difference value of the steam supply main pressure trend = g0x the difference value of the steam supply main steam flow.

[0095] Preferably,

[0096] The control solution space C in the thermal power production system is obtained based on the pre-acquired set of adjustable values ​​corresponding to each group of steam supply equipment in the thermal power production system, specifically including:

[0097] The Cartesian product of the adjustable value set of each group of steam supply equipment in the thermal power production system is taken as the final control solution space C;

[0098] Where C = C1 × C2 ... × C i ...×C N ;

[0099] C i This is the set of adjustable values ​​corresponding to the i-th group of steam supply equipment in the thermal power production system;

[0100] C i This includes M pre-set steam flow adjustment values ​​in the main steam header corresponding to the i-th group of steam supply equipment;

[0101] In this case, each row of data in the final control solution space C is considered as a set of adjusted feasible solutions.

[0102] Preferably, A4 specifically includes:

[0103] A41. For any feasible solution for adjusting the steam flow rate in the main steam header of any group e in the final control solution space C: [c1, c2, ..., c i , ...c N [; Obtain the steam flow adjustment solution in the steam supply header corresponding to the feasible solution for adjusting the steam flow in the main steam header of this group.]

[0104] in,

[0105] Among them, c i This is the steam flow adjustment value in the main steam header corresponding to the i-th group of steam supply equipment in the thermal power production system;

[0106] c i * This is the steam flow adjustment value in the steam supply header corresponding to the i-th group of steam supply equipment in the thermal power production system.

[0107] w i It is the i-th c in the pre-obtained e-th feasible solution group. i The corresponding coefficient;

[0108] A42. The main steam header pressure at time t and the main steam header pressure trend at time t. a-t And the flow adjustment value c from the main steam header iThe flow variation amount in the main steam header is substituted into the first recursive formula to obtain the new main steam header pressure pa t+1 * ;

[0109] pa t+1 * = pa t +trend a·t +f*c i ;

[0110] A43, based on the main steam header pressure set S at time t ai-t and the new main steam header pressure pa t+1 * , the third recursive formula is used to obtain the new main steam header pressure after the second time interval to the kth time interval after the current time t, and the new main steam header pressure after the k time intervals after the time t forms the final main steam header pressure prediction sequence S ai *;

[0111] wherein the third recursive formula is:

[0112]

[0113] wherein, is a set of new main steam header pressures corresponding to the k-1 time intervals after the first specified time t, which is obtained by using a predetermined fitting method with a fixed index sequence [1, 2, 3...T] as the independent variable and the set S ai-t-k-1 * as the dependent variable;

[0114] when k = 1, is obtained by the first recursive formula; when k > 1, is obtained by the third recursive formula;

[0115] The final main steam header pressure prediction sequence S ai * is:

[0116] S ai * = [pa* t+1 , pa* t+2 , pa* t+3 , …, pa* t+k ];

[0117] A44, the steam supply header pressure corresponding to time t, the steam supply header pressure trend trend a·t at time t, and the flow adjustment value c i *The flow variation amount in the steam supply main pipe is substituted into the second recursive formula to obtain the new steam supply main pipe pressure corresponding to the time t+Δ

[0118] wherein,

[0119] A45, based on the steam supply main pipe pressure set S at the time t bi-t and the new steam supply main pipe pressure corresponding to the time t+Δ The fourth recursive formula is used to obtain the new steam supply main pipe pressure in the second time interval to the kth time interval after the current time t, and the new steam supply main pipe pressure in the k time intervals after the time t forms the final steam supply main pipe pressure prediction sequence S bi *;

[0120] wherein, the fourth recursive formula is:

[0121]

[0122] wherein, is a set of new steam supply main pipe pressures corresponding to k-1 time intervals after the first specified time t with a fixed index sequence [1, 2, 3...T] as the independent variable, and is obtained by using a pre-set fitting method; bi-t-k-1 *;

[0123] when k=1, is obtained by the second recursive formula; when k>1, is obtained by the fourth recursive formula;

[0124] The final steam supply main pipe pressure prediction sequence S bi *is:

[0125] S bi *=[pb* t+1 , pb* t+2 , pb* t+3 , …, pb* t+k ].

[0126] A46, for any each group of main steam supply pipe steam flow adjustment feasible solution corresponding to the final main steam supply pipe pressure prediction sequence S ai *in the final control solution space C and the final steam supply main pipe pressure prediction sequence S bi *is screened to determine the optimal recommended value.

[0127] Preferably, the A46 specifically comprises:

[0128] A46-1, filtering the final main steam header pressure prediction sequence and the final steam header pressure prediction sequence corresponding to the i-th group of steam supply equipment corresponding to any each group of main steam header steam flow adjustment feasible solution in the final control solution space C, to obtain the final main steam header pressure prediction sequence and the final steam header pressure prediction sequence corresponding to the i-th group of steam supply equipment that meet the preset target condition;

[0129] Wherein, the target condition is:

[0130] Any new main steam header pressure value in the final main steam header pressure prediction sequence is greater than or equal to the first preset value lower i1 corresponding to the i-th group of steam supply equipment, and less than or equal to the second preset value upper i1 corresponding to the i-th group of steam supply equipment;

[0131] Any new steam header pressure value in the final steam header pressure prediction sequence is greater than or equal to the third preset value lower i2 corresponding to the i-th group of steam supply equipment, and less than or equal to the fourth preset value upper i2 corresponding to the i-th group of steam supply equipment;

[0132] A46-2, for the final main steam header pressure prediction sequence corresponding to the i-th group of steam supply equipment that meets the preset target condition, obtaining the main steam header pressure convergence factor μ i1 corresponding to the i-th group of steam supply equipment according to formula (1);

[0133] The formula (1) is:

[0134]

[0135] For the final steam header pressure prediction sequence corresponding to the i-th group of steam supply equipment that meets the preset target condition, obtaining the steam header pressure convergence factor μ i2 corresponding to the i-th group of steam supply equipment according to formula (2);

[0136] The formula (2) is:

[0137]

[0138] A46-3, according to the main steam header pressure convergence factor and the steam header pressure convergence factor corresponding to each group of steam supply equipment that meets the preset target condition, obtaining the first pressure convergence value μ corresponding to the thermal power production system by formula (3);

[0139] Wherein, the formula (3) is:

[0140] A46-4, the steam flow adjustment value in the main steam header and the steam flow adjustment value in the steam supply header corresponding to the final control solution space C corresponding to the first pressure convergence value μmin are taken as the optimal adjustment value of the main steam header and the optimal adjustment value of the steam supply header, respectively.

[0141] (III) Beneficial Effects

[0142] The beneficial effects of the present application are as follows: the control method of the steam supply equipment in the thermal power production system of the present application, since the main steam header pressure prediction sequence S a·t and the steam supply header pressure prediction sequence S ai are obtained based on the main steam header pressure corresponding to time t, the main steam header pressure trend trend b·t at time t, and the first recursive formula, and the steam supply header pressure corresponding to time t, the steam supply header pressure trend trend bi at time t, and the second recursive formula, and the optimal adjustment value of the main steam header and the optimal adjustment value of the steam supply header are obtained based on the main steam header pressure prediction sequence S ai , the steam supply header pressure prediction sequence S bi , and the control solution space C corresponding to the i-th group of steam supply equipment groups in the thermal power production system obtained in advance; the optimal adjustment value is obtained according to the control solution space C and the predicted header pressure trend, which improves the control accuracy of the header pressure; at the same time, the control accuracy of the header pressure is improved by adjusting the main steam header and the steam supply header. BRIEF DESCRIPTION OF DRAWINGS

[0143] Figure 1 The flow chart of the control method of the steam supply equipment in the thermal power production system of the present application;

[0144] Figure 2 The schematic diagram of one embodiment of the steam supply equipment group in the present application;

[0145] Figure 3 The schematic diagram of another embodiment of the steam supply equipment group in the present application.

[0146]

Explanation of reference signs

[0147] 1: main steam header;

[0148] 2: steam supply header. DETAILED DESCRIPTION

[0149] In order to better explain the present application and facilitate understanding, the present application is described in detail in combination with the drawings through specific embodiments.

[0150] To better understand the above technical solutions, exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that the present invention can be understood more clearly and thoroughly, and that the scope of the present invention can be fully conveyed to those skilled in the art.

[0151] See Figure 1 This embodiment provides a control method for steam supply equipment in a thermal power production system, including:

[0152] A1. Collect the set S of main steam header pressures corresponding to the i-th group of steam supply equipment in the thermal power production system and the first specified time t according to a preset time interval. ai-t The set of pressures S of the steam supply header bi-t The set S of main steam header pressures corresponding to the first specified time t ai-t This includes the main steam header pressure collected sequentially within a first specified time period before the first specified time t; and the set S of the steam supply header pressure corresponding to the first specified time t. bi-t This includes the steam supply main pipe pressure collected sequentially within the first specified time period before the first specified time t;

[0153] The first specified time period is [t-(T+1)Δ, t]; Δ is a preset time interval;

[0154] Specifically, the thermal power production system includes N groups of steam supply equipment;

[0155] Each steam supply equipment group includes a first steam supply equipment, see [link to relevant documentation]. Figure 2 Or it may include multiple second steam supply devices that meet the first connection conditions, such as Figure 3 As shown.

[0156] Among them, see Figure 3 The first connection condition is that the main steam headers 1 connected in front of all the second steam supply equipment are interconnected, and the steam supply headers 2 connected behind all the second steam supply equipment are interconnected.

[0157] See Figure 2 The first steam supply equipment is a steam supply equipment that meets the second connection condition; the second connection condition is that the main steam header 1 connected in front of the steam supply equipment is not connected to the main steam header 1 of other steam supply equipment, and the steam supply header 2 connected behind it is not connected to the steam supply header 2 of other steam supply equipment.

[0158] Wherein, S is the set of main steam header pressures corresponding to the first specified time t. ai-t for:

[0159] S ai-t = [pa t-T+1 ,... pa t-c ,... pa t-2 , pa t-1 , pa t ].

[0160] T is any positive integer between [3, 10].

[0161] pa t-c is the main steam header pressure at the c time intervals before the first specified time t.

[0162] The set S bi-t of the main steam header pressure corresponding to the first specified time t is:

[0163] S bi-t = [pb t-T+1 ,... pb t-c ,... pb t-2 , pb t-1 , pb t ].

[0164] pb t-c is the main steam header pressure at the c time intervals before the first specified time t.

[0165] A2, for the set S ai-t of the main steam header pressure corresponding to the first specified time t, and the set S bi-t of the main steam header pressure, respectively adopts a pre-set strategy to obtain the main steam header pressure trend trend a·t and the main steam header pressure trend trend b·t at time t, respectively.

[0166] In the specific application of the embodiment, in the A2, for the set S ai-t of the main steam header pressure corresponding to the first specified time t, a pre-set strategy is adopted to obtain the main steam header pressure trend trend a·t at time t, specifically including:

[0167] Taking the fixed index sequence [1, 2, 3...T] as the independent variable, and the set S ai-t of the main steam header pressure corresponding to the first specified time t as the dependent variable, a pre-set fitting method is adopted to obtain the main steam header pressure trend trend a·t at time t, specifically including: using a T-order polynomial to fit the independent variable and the dependent variable to obtain a first polynomial function corresponding to the dependent variable.

[0168] According to the first polynomial function corresponding to the dependent variable, a first derivative Pa'(t) is obtained, and a value at T in the first derivative Pa'(t) is taken as the main steam header pressure trend trend at the time t a·t .

[0169] In the A2, a set S of the main steam header pressure corresponding to the first specified time t bi-t , a pre-set strategy is used to obtain the main steam header pressure trend trend at the time t b·t , specifically including:

[0170] The fixed index sequence [1, 2, 3...T] is taken as the independent variable, and the set S of the main steam header pressure corresponding to the first specified time t bi-t is taken as the dependent variable, and a pre-set fitting method is used to obtain the main steam header pressure trend trend at the time t b·t , specifically including: using a T-order polynomial to fit the independent variable and the dependent variable to obtain a second polynomial function corresponding to the dependent variable.

[0171] According to the second polynomial function corresponding to the dependent variable, a first derivative Pb'(t) is obtained, and a value at T in the first derivative Pb'(t) is taken as the main steam header pressure trend trend at the time t b·t .

[0172] A3, based on the main steam header pressure corresponding to the time t, the main steam header pressure trend trend at the time t a·t and the first recursive formula, the main steam header pressure prediction sequence S ai is obtained, and based on the main steam header pressure corresponding to the time t, the main steam header pressure trend trend at the time t b·t and the second recursive formula, the main steam header pressure prediction sequence S bi .

[0173] Among them, the A3 specifically includes:

[0174] A3-1, based on the main steam header pressure corresponding to the time t, the main steam header pressure trend trend at the time t a·t and the first recursive formula, the main steam header pressure prediction sequence S ai , specifically including:

[0175] A3-1-1, the main steam header pressure corresponding to the time t, the main steam header pressure trend trend at the time t a·t and the pre-obtained flow change amount dsteam in the main steam header at the time t a·t are substituted into the first recursive formula to obtain the main steam header pressure pa corresponding to the time t+Δ t+1 .

[0176] wherein pa t+1 = pa t + trend a·t + f * dsteam a·t .

[0177] The first recursive formula is:

[0178] pa t = pa t+k-1 + trend a·t+k-1 + f * dsteam a·t+k-1 .

[0179] A3-1-2, a set S of main steam header pressures based on time t ai-t and the main steam header pressure pa t+1 corresponding to time t+Δ, using the first recursive formula, the main steam header pressures of the second time interval to the kth time interval after time t are obtained respectively, and the main steam header pressure prediction sequence S ai is composed of the main steam header pressures of the k time intervals after time t.

[0180] wherein trend a·t+k-1 is a set S ai-t-k-1 of main steam header pressures corresponding to the k-1 time intervals after the first specified time t, using a pre-set fitting method, with a fixed index sequence [1, 2, 3...T] as the independent variable.

[0181] When k is greater than or equal to 2, pa t+k-1 is obtained by using the first recursive formula.

[0182] f is a model coefficient obtained in advance in the first recursive formula.

[0183] dsteam a·t+k-1 = steam a·t+k-1 - steam a·t+k-2 .

[0184] dsteam a·t+k-1 is the flow change amount in the main steam header in the k-1 time interval after time t.

[0185] steam a·t+k-1 is the flow in the main steam header in the k-1 time interval after time t.

[0186] The main steam header pressure prediction sequence S ai is:

[0187] S ai = [pa t+1 , pat+2 , pa t+3 , …, pa t+k ].

[0188] In the embodiment, the model coefficient f in the first recursive formula is obtained in the following manner:

[0189] From the beginning to the end of the second specified time period, a time window with a length of L*Delta is used to sequentially shift H times from the beginning to the end of the second specified time period, and each time the time window is shifted by Delta.

[0190] The second specified time period is [t0-(H+L)*Delta, t0].

[0191] Wherein, t0 is the second specified time; H is a pre-set number; L is a random integer between [6, 12].

[0192] After each shift of the time window, a set of main steam header steam flow at the time corresponding to the beginning of the time window and a set of main steam header pressure at the time corresponding to the beginning of the time window are obtained, and a set of main steam header steam flow at the time corresponding to the end of the time window and a set of main steam header pressure at the time corresponding to the end of the time window are obtained.

[0193] For the main steam header steam flow at the time corresponding to the beginning of the time window and the main steam header steam flow at the time corresponding to the end of the time window obtained after the dth shift of the time window, the difference between the main steam header steam flow corresponding to the dth shift of the time window is obtained; wherein, 1≤d≤H.

[0194] Taking the fixed index sequence [1, 2, 3...T] as the independent variable and the set of main steam header pressure at the time corresponding to the beginning of the time window obtained after the dth shift of the time window as the dependent variable, a pre-set fitting method is used to obtain the main steam header pressure trend at the time corresponding to the beginning of the time window obtained after the dth shift of the time window.

[0195] Taking the fixed index sequence [1, 2, 3...T] as the independent variable and the set of main steam header pressure at the time corresponding to the end of the time window obtained after the dth shift of the time window as the dependent variable, a pre-set fitting method is used to obtain the main steam header pressure trend at the time corresponding to the end of the time window obtained after the dth shift of the time window.

[0196] Based on the main steam header pressure trend at the time corresponding to the beginning of the time window obtained after the dth shift of the time window and the main steam header pressure trend at the time corresponding to the end of the time window obtained after the dth shift of the time window, the difference between the main steam header pressure trend corresponding to the dth shift of the time window is obtained.

[0197] The first two-dimensional array is generated based on the difference in steam flow rate of the main steam header and the difference in pressure trend of the main steam header after each shift of the time window.

[0198] The first two-dimensional array is:

[0199]

[0200] Among them, delta_steam ad delta_trend represents the difference in steam flow rate of the main steam header after the d-th shift of the time window. ad It represents the difference in the pressure trend of the main steam header corresponding to the d-th shift of the time window.

[0201] The difference in steam flow rate of the main steam header and the difference in pressure trend of the main steam header corresponding to each translation in the first two-dimensional array are substituted into the pre-acquired first non-intercept linear regression model to train and obtain the specific value of the regression model coefficient f0 in the first non-intercept linear regression model. The specific value of the regression model coefficient f0 in the first non-intercept linear regression model is used as the model coefficient f in the first recursive formula.

[0202] The first non-intercept linear regression model is as follows:

[0203] The difference in pressure trend of the main steam header = f0 × the difference in steam flow rate of the main steam header.

[0204] A3-2. Based on the steam supply header pressure at time t, and the trend of steam supply header pressure at time t. b-t The second recursive formula is used to obtain the steam supply header pressure prediction sequence S. bi Specifically, it includes:

[0205] A3-2-1. The steam supply header pressure at time t and the steam supply header pressure trend at time t. b·t And the change in flow rate dsteam in the steam supply header at time t b·t Substituting into the second recursive formula, we obtain the steam supply header pressure pb at time t+Δ. t+1 .

[0206] Among them, pb t+1 =pb t +trend b·t +g*dsteam b·t .

[0207] The second recursive formula is:

[0208] pbt+k=pb t+k-1 +trendb·t+k-1 + g * dsteam b·t+k-1 .

[0209] A3-2-2, a set S of the steam main pressures at time t bi-t and the steam main pressure pb corresponding to time t + Δ t+1 , a second recursive formula is used to obtain the steam main pressures after the second time interval to the kth time interval after the current time t, and the steam main pressure prediction sequence S is composed of the steam main pressures after the k time intervals after time t bi .

[0210] wherein trend b·t+k-1 is obtained by using a predetermined fitting method with a fixed index sequence [1, 2, 3... T] as the independent variable and a set S of the steam main pressures corresponding to the k-1 time intervals after the first specified time t bi-t-k-1 as the dependent variable.

[0211] pb t+k-1 is obtained by using the second recursive formula.

[0212] g is a model coefficient obtained in advance in the second recursive formula.

[0213] dsteam b·t+k-1 = steam b·t+k-1 - steam b·t+k-2 .

[0214] dsteam b·t+k-1 is the flow rate change amount in the steam main after the k-1 time intervals after time t.

[0215] steam b·t+k-1 is the flow rate in the steam main after the k-1 time intervals after time t.

[0216] steam main pressure prediction sequence S bi is:

[0217] S bi = [pb t+1 , pb t+2 , pb t+3 , …, pb t+k ].

[0218] In the embodiment, the model coefficient g in the second recursive formula is obtained in the following manner:

[0219] From the beginning to the end of the second specified time period, a time window with a length of L x Δ is used to translate from the beginning to the end of the second specified time period H times in turn, and each translation is Δ.

[0220] After each time window translation, a set of the steam supply main steam flow at the time corresponding to the start of the time window and a set of the steam supply main steam pressure at the time corresponding to the start of the time window, and a set of the steam supply main steam flow at the time corresponding to the end of the time window and a set of the steam supply main steam pressure at the time corresponding to the end of the time window are obtained.

[0221] For the steam supply main steam flow at the time corresponding to the start of the time window and the steam supply main steam flow at the time corresponding to the end of the time window obtained after the dth time window translation, a difference value of the steam supply main steam flow corresponding to the dth time window translation is obtained.

[0222] Taking the fixed index sequence [1, 2, 3...T] as the independent variable and the set of the steam supply main steam pressure at the time corresponding to the start of the time window obtained after the dth time window translation as the dependent variable, a steam supply main steam pressure trend at the time corresponding to the start of the time window obtained after the dth time window translation is obtained by using a pre-set fitting mode.

[0223] Taking the fixed index sequence [1, 2, 3...T] as the independent variable and the set of the steam supply main steam pressure at the time corresponding to the end of the time window obtained after the dth time window translation as the dependent variable, a steam supply main steam pressure trend at the time corresponding to the end of the time window obtained after the dth time window translation is obtained by using a pre-set fitting mode.

[0224] Based on the steam supply main steam pressure trend at the time corresponding to the start of the time window obtained after the dth time window translation and the steam supply main steam pressure trend at the time corresponding to the end of the time window obtained after the dth time window translation, a difference value of the steam supply main steam pressure trend corresponding to the dth time window translation is obtained.

[0225] According to the difference value of the steam supply main steam flow and the difference value of the steam supply main steam pressure trend corresponding to each time window translation, a second two-dimensional array is generated.

[0226] The second two-dimensional array is as follows:

[0227]

[0228] delta_steam bd is the difference value of the steam supply main steam flow corresponding to the dth time window translation.

[0229] delta_trend bd is the difference value of the steam supply main steam pressure trend corresponding to the dth time window translation.

[0230] The difference value of the steam flow of the steam supply main corresponding to each translation in the second two-dimensional array and the difference value of the steam supply main pressure trend are substituted into the second linear regression model without intercept obtained in advance to train the specific value of the regression model coefficient g0 in the second linear regression model without intercept, and the specific value of the regression model coefficient g0 in the second linear regression model without intercept is used as the model coefficient g in the second recursive formula.

[0231] The second linear regression model without intercept is:

[0232] The difference value of the steam supply main pressure trend = g0 x the difference value of the steam flow of the steam supply main.

[0233] A4, based on the main steam supply main pressure prediction sequence S ai , the main steam supply main pressure prediction sequence S bi , and the control solution space C corresponding to the i-th group of steam supply equipment groups in the thermal power production system obtained in advance, the optimal adjustment value of the main steam supply main and the optimal adjustment value of the steam supply main are obtained.

[0234] The corresponding control solution space C in the thermal power production system is obtained according to the adjustable value set corresponding to each group of steam supply equipment groups in the thermal power production system, and specifically includes:

[0235] The Cartesian product of the adjustable value set of each group of steam supply equipment groups in the thermal power production system is used as the final control solution space C. Wherein, C = C1 x C2... x C i ... x C N .

[0236] C i is the adjustable value set corresponding to the i-th group of steam supply equipment groups in the thermal power production system. C i includes M steam flow adjustment values preset in the main steam supply main corresponding to the i-th group of steam supply equipment groups. Wherein, each row of data in the final control solution space C is used as a group of adjustment feasible solutions.

[0237] For example, assuming that there are 2 groups of steam supply equipment groups in the thermal power production system, C1 represents the adjustable value set corresponding to the first group of steam supply equipment groups, C2 represents the adjustable value set corresponding to the second group of steam supply equipment groups, and if the adjustable value set corresponding to the first group of steam supply equipment groups has a total of 3 steam flow adjustment values (M = 3 at this time), they are respectively: 11 , a 12 , a 13 If the adjustable value set corresponding to the second group of steam supply equipment groups has a total of 3 steam flow adjustment values, they are respectively: 21 , b 22 , b 23Therefore, the control solution space C corresponding to the thermoelectric production system at this time is:

[0238] In this system, each row of data in the corresponding control solution space C of the thermal power production system is considered as a set of adjustable feasible solutions. For example, if the first row of data in the final control solution space C is a... 11 b 21 Then a 11 and b 21 , which is a set of adjustable feasible solutions, where a 11 This corresponds to the steam flow adjustment value in the main steam header corresponding to the first group of steam supply equipment, b. 21 This corresponds to the steam flow adjustment value in the main steam header corresponding to the second group of steam supply equipment.

[0239] In the practical application of this embodiment, A4 specifically includes:

[0240] A41. For any feasible solution for adjusting the steam flow rate in the main steam header of any group e in the final control solution space C: [c1, c2, ... c i , ...c N [; Obtain the steam flow adjustment solution in the steam supply header corresponding to the feasible solution for adjusting the steam flow in the main steam header of this group.]

[0241] in,

[0242] Among them, c i This is the steam flow adjustment value in the main steam header corresponding to the i-th group of steam supply equipment in the thermal power production system.

[0243] c i * This is the steam flow adjustment value in the steam supply header corresponding to the i-th group of steam supply equipment in the thermal power production system.

[0244] w i It is the i-th c in the pre-obtained e-th feasible solution group. i The corresponding coefficient.

[0245] A42. The main steam header pressure at time t and the main steam header pressure trend at time t. a-t And the flow adjustment value c from the main steam header i Substituting the change in flow rate within the main steam header into the first recursive formula, we obtain the new main steam header pressure pa at time t+Δ. t+1 * .

[0246] Among them, pa t+1 *= pa t + trend a·t + f * c i .

[0247] A43, based on the main steam header pressure at time t set S ai-t And the new main steam header pressure pa corresponding to time t+Δ t+1 * , using the third recursive formula, respectively, to obtain the new main steam header pressure after the first 2 time interval to the k time interval, and by the k time interval after time t of the new main steam header pressure consisting of the final main steam header pressure prediction sequence S ai *.

[0248] Wherein, the third recursive formula is:

[0249]

[0250] Wherein, Is a fixed index sequence [1, 2, 3...T] as the independent variable, the first specified time t after k-1 time interval corresponding to the new main steam header pressure set S ai-t-k-1 * for the dependent variable, using the pre-set fitting method to obtain.

[0251] When k = 1, Is obtained by the first recursive formula; when k > 1, Is obtained by the third recursive formula.

[0252] The final main steam header pressure prediction sequence S ai * is:

[0253] S ai * = [pa* t+1 , pa* t+2 , pa* t+3 , …, pa* t+k ].

[0254] A44, the steam supply header pressure corresponding to time t, steam supply header pressure trend a·t At time t and the flow adjustment value c i * As the flow change in the steam supply header, into the second recursive formula, to obtain the new steam supply header pressure corresponding to time t+Δ

[0255] Wherein,

[0256] A45, based on the steam supply header pressure at time t set S bi-tThe new supply main pipe pressure corresponding to the time t+Δ The fourth recursive formula is used to obtain the new supply main pipe pressure of the second time interval to the kth time interval after the current time t, and the final supply main pipe pressure prediction sequence S is composed of the new supply main pipe pressure of the k time intervals after the time t bi

[0257] The fourth recursive formula is:

[0258]

[0259] The fourth recursive formula is: is a set of new supply main pipe pressures corresponding to the k-1 time intervals after the first specified time t, with a fixed index sequence [1, 2, 3...T] as the independent variable and the dependent variable obtained by using a pre-set fitting method. bi-t-k-1

[0260] When k=1, is obtained by the second recursive formula; when k>1, is obtained by the fourth recursive formula.

[0261] The final supply main pipe pressure prediction sequence S bi

[0262] bi t+1 t+2 t+3 t+k

[0263] A46, for any each group of the final control solution space C in the steam flow adjustment feasible solution of the main steam pipe, the final main steam pipe pressure prediction sequence S ai corresponding to the i-th group of the steam supply equipment group and the final supply main pipe pressure prediction sequence S bi are screened to determine the optimal recommended value.

[0264] Specifically, the A46 specifically includes:

[0265] A46-1, for any each group of the final control solution space C in the steam flow adjustment feasible solution of the main steam pipe, the final main steam pipe pressure prediction sequence S

[0266] corresponding to the i-th group of the steam supply equipment group and the final supply main pipe pressure prediction sequence S are screened to obtain the final main steam pipe pressure prediction sequence S ​​​​​​​​​

[0267] Any new value of the final main steam header pressure prediction sequence is greater than or equal to the first preset value lower corresponding to the i-th group of steam supply equipment i1 , and less than or equal to the second preset value upper corresponding to the i-th group of steam supply equipment i1 .

[0268] Any new value of the final main steam header pressure prediction sequence is greater than or equal to the first preset value lower corresponding to the i-th group of steam supply equipment i2 , and less than or equal to the second preset value upper corresponding to the i-th group of steam supply equipment i2 .

[0269] A46-2, for the final main steam header pressure prediction sequence corresponding to the i-th group of steam supply equipment meeting the preset target condition, the main steam header pressure convergence factor μ corresponding to the i-th group of steam supply equipment is obtained according to formula (1) i1 .

[0270] The formula (1) is:

[0271]

[0272] For the final main steam header pressure prediction sequence corresponding to the i-th group of steam supply equipment meeting the preset target condition, the main steam header pressure convergence factor μ corresponding to the i-th group of steam supply equipment is obtained according to formula (1) i2 .

[0273] The formula (2) is:

[0274]

[0275] A46-3, according to the main steam header pressure convergence factor and the steam supply header pressure convergence factor corresponding to each group of steam supply equipment meeting the preset target condition, the first pressure convergence value μ corresponding to the thermal power production system is obtained by using formula (3).

[0276] The formula (3) is:

[0277] A46-4, the main steam header steam flow adjustment value in the final control solution space C corresponding to the minimum first pressure convergence value μ and the steam supply header steam flow adjustment value corresponding thereto are respectively taken as the optimal adjustment value of the main steam header and the optimal adjustment value of the steam supply header.

[0278] The control method of the steam supply equipment in the thermal power production system in the embodiment is characterized in that: the main steam header pressure corresponding to the time t, the main steam header pressure trend trend a·t at the time t, and the first recursive formula are used to obtain the main steam header pressure prediction sequence S ai The steam supply header pressure corresponding to the time t, the steam supply header pressure trend trend b·t at the time t, and the second recursive formula are used to obtain the steam supply header pressure prediction sequence S bi The optimal adjustment value of the main steam header and the optimal adjustment value of the steam supply header are obtained based on the main steam header pressure prediction sequence S ai , the steam supply header pressure prediction sequence S bi , and the control solution space C corresponding to the i-th group of steam supply equipment groups in the thermal power production system obtained in advance.

[0279] The control system of the steam supply equipment in the thermal power production system in the embodiment is characterized in that: the control system comprises at least one processor and at least one memory in communication connection with the processor, wherein the memory stores program instructions executable by the processor, and the processor calling the program instructions can execute the control method of the steam supply equipment in the thermal power production system as described above.

[0280] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system or a computer program product. Therefore, the present application can be in the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can be in the form of a computer program product implemented on one or more computer usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer usable program codes.

[0281] The present application is described with reference to the flowcharts and / or block diagrams of the method, device (system) and computer program product according to the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of the flows and / or blocks in the flowcharts and / or block diagrams can be realized by computer program instructions.

[0282] It should be noted that the description using "comprises" or "comprising" does not exclude the presence of elements or steps other than those listed in a claim. The words "a" or "an" preceding the disclosure of a plurality of elements or steps do not exclude the presence of a plurality of such elements or steps. The application can be implemented by means of both hardware and software, and any combination thereof. In a claim reciting a means, the term "comprising" is to be interpreted as including the elements or steps recited in the claim, not as excluding any additional elements or steps. The term "first", "second", "third" and the like in the description do not necessarily denote any ordinal. The terms "first", "second", "third" and the like can be understood as names.

[0283] Furthermore, it is noted that the described features, structures, or characteristics can be combined in any suitable manner in one or more embodiments or examples. It will also be readily appreciated that the terms "comprise", "comprising", "include", "including", "has", "having" or variants thereof are not to be construed as limiting the described embodiments or examples to the precise configurations recited. In describing and claiming the embodiments or examples, "including" and "comprising" are open-ended terms that are intended to permit the inclusion of additional features or steps without materially departing from the spirit or scope of the invention. It will also be readily appreciated that the terms "first", "second", "third", etc. are used merely as labels, and are not intended to impose numerical requirements on their objects.

[0284] While the preferred embodiments of the application have been described, additional variations and modifications can be made to the embodiments without departing from the spirit and scope of the application. Therefore, it should be understood that the application is not limited by the preferred embodiments and examples described above. Rather, the scope of the application is defined by the appended claims and their equivalents.

[0285] 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 control method of a steam supply device in a thermoelectric power production system, characterized by, Comprising: A1, collecting a set S of main steam header pressures corresponding to a first specified time t in the thermal power production system according to a pre-set time interval ai-t and the set S of main steam header pressures corresponding to the first specified time t bi-t ; the set S of main steam header pressures corresponding to the first specified time t ai-t , including the main steam header pressures collected in sequence within a first specified time period before the first specified time t; and the set S of main steam header pressures corresponding to the first specified time t bi-t , including the main steam header pressures collected in sequence within a first specified time period before the first specified time t The first specified time period is [t-(T+1)Δ, t]; Δ is a pre-set time interval; A2, a set S of main steam header pressures corresponding to the first specified time t ai-t and a set S of boiler feedwater header pressures bi-t respectively using a pre-set strategy, to obtain a main steam header pressure trend trend a·t and a boiler feedwater header pressure trend trend b·t ; A3, based on the main steam header pressure corresponding to time t, the main steam header pressure trend trend a·t and the first recursive formula, obtain the main steam header pressure prediction sequence S ai , and based on the main steam header pressure corresponding to time t, the main steam header pressure trend trend b·t and the second recursive formula, obtain the main steam header pressure prediction sequence S bi ; A4、based on the main steam header pressure prediction sequence S ai , the steam supply header pressure prediction sequence S bi and the pre-acquired control solution space C corresponding to the i-th group of steam supply equipment groups in the thermal power generation system, the optimal adjustment value of the main steam header and the optimal adjustment value of the steam supply header are obtained.

2. The control method of the steam supply equipment in the thermal power production system according to claim 1, characterized in that, The thermal power production system comprises N groups of steam supply equipment groups; Any steam supply equipment group comprises one first steam supply equipment or a plurality of second steam supply equipment meeting the first communication condition; The first communication condition is that the main steam headers connected in front of all the second steam supply equipment are in communication with each other, and the steam supply headers connected behind all the second steam supply equipment are in communication with each other; The first steam supply equipment is a steam supply equipment meeting the second communication condition; the second communication condition is that the main steam header connected in front of the steam supply equipment is not in communication with the main steam header of other steam supply equipment, and the steam supply header connected behind the steam supply equipment is not in communication with the steam supply header of other steam supply equipment; wherein the set S of main steam header pressures corresponding to the first specified time instant t ai-t is: S ai-t = [pa t-T+1 ,... pa t-c ,... pa t-2 , pa t-1 , pa t ]; T is any positive integer between 3 and 10; pa t-c Pc is the main steam header pressure at the c time intervals before the first specified time t; A set S of steam main pressures corresponding to a first specified time t bi-t is: S bi-t = [pb t-T+1 ,... pb t-c ,... pb t-2 , pb t-1 , pb t ]; pb t-c Pb is the supply main pressure at the c time intervals before the first specified time t.

3. The control method of the steam supply equipment in the thermal power production system according to claim 2, characterized in that, In the A2, the main steam header pressure corresponding to the first specified time t is obtained ai-t The main steam header pressure trend trend at time t is obtained by using a preset strategy a.t , specifically comprising: Taking the fixed index sequence [1, 2, 3...T] as an independent variable, taking the set S of the main steam header pressure corresponding to the first specified time t ai-t Taking the set S of the main steam header pressure corresponding to the first specified time t as a dependent variable, a fitting method is used to obtain the main steam header pressure trend trend at time t a.t , specifically including: using a T-order polynomial to fit the independent variable and the dependent variable to obtain a first polynomial function corresponding to the dependent variable; According to the first polynomial function corresponding to the dependent variable, a first derivative Pa'(t) of the first polynomial function is obtained, and a value of the first derivative Pa'(t) at T is taken as a main steam header pressure trend trend a·t ; In the A2, for the first designated time t, the set S of the steam main pressure corresponding to the set S bi-t , using a pre-set strategy, the steam main pressure trend trend at time t is obtained b.t , specifically including: Taking the fixed index sequence [1, 2, 3...T] as an independent variable, taking the set S of the steam supply main pressure corresponding to the first specified moment t bi-t Taking the set S of the steam supply main pressure corresponding to the first specified moment t as a dependent variable, a fitting mode is adopted to obtain the steam supply main pressure trend trend of the moment t b.t , specifically including: using a T-order polynomial to fit the independent variable and the dependent variable to obtain a second polynomial function corresponding to the dependent variable; According to the second polynomial function corresponding to the dependent variable, a first derivative Pb'(t) is obtained, and a value at T in the first derivative Pb'(t) is taken as the steam main pressure trend trend b.t .

4. The control method of the steam supply device in a thermoelectric power production system according to claim 3, characterized by, The A3 specifically comprises: A3-1, based on the main steam header pressure corresponding to the time t, the main steam header pressure trend trend at time t a·t and the first recursive formula, obtain the main steam header pressure prediction sequence S ai , specifically comprising: A3-1-1, the main steam pipe pressure corresponding to time t, the main steam pipe pressure trend at time t trend a.t and the flow change amount dsteam in the main steam pipe at time t obtained in advance a.t Substituting into the first recursive formula, the main steam pipe pressure pa corresponding to time t+Δ is obtained t+1 ; where pa t+1 = pa t + trend a·t + f * dsteam a·t ; The first recursive formula is: pa t+k = pa t+k-1 + trend a·t+k-1 + f * dsteam a·t+k-1 ; A3-1-2, a set S of main steam header pressures based on time t ai-t and the main steam header pressure pa corresponding to time t+Δ t+1 , a first recursive formula is used to obtain the main steam header pressures of the second time interval to the kth time interval after time t, and the main steam header pressure prediction sequence S is composed of the main steam header pressures of the k time intervals after time t ai ; wherein trend a·t+k-1 is obtained by using a preset fitting method, with a fixed index sequence [1, 2, 3...T] as an independent variable and a set S of the main steam header pressure corresponding to k-1 time intervals after the first specified time t as a dependent variable ai-t-k-1 is obtained by using a preset fitting method, with a fixed index sequence [1, 2, 3...T] as an independent variable and a set S of the main steam header pressure corresponding to k-1 time intervals after the first specified time t as a dependent variable k is greater than or equal to 2, pa t+k-1 is recursively calculated using the first recursive formula; f is a model coefficient pre-acquired in the first recursive formula; dsteam a·t+k-1 = steam a·t+k-1 - steam a·t+k-2 ; dsteam a·t+k-1 is the flow rate variation of the main steam header for the k-1 time interval after the t time. steam a·t+k-1 Qk-1is the flow rate in the main steam header for the k-1th time interval after time t; Main steam header pressure prediction sequence S ai is: S ai = [pa t+1 , pa t+2 , pa t+3 ,..., pa t+k ]; A3-2, based on the time t corresponding to the main steam supply pipe pressure, the time t of the main steam supply pipe pressure trend trend b·t and the second recursive formula, obtain the main steam supply pipe pressure prediction sequence S bi , specifically comprising: A3-2-1, the steam supply main pressure corresponding to the time t, the steam supply main pressure trend trend at the time t b·t and the flow change amount dsteam in the steam supply main at the time t b·t Substituting the second recursive formula, the steam supply main pressure pb corresponding to the time t+Δ is obtained t+1 ; where pb t+1 = pb t + trend b·t + g * dsteam b·t ; The second recursive formula is: pb t+k = pb t+k-1 + trend b·t+k-1 + g * dsteam b·t+k-1 ; A3-2-2, a set S of the supply main pipe pressures based on time t bi-t and the supply main pipe pressure pb corresponding to time t+Δ t+1 , the second recursive formula is used to obtain the supply main pipe pressures of the second time interval to the kth time interval after the current time t, and the supply main pipe pressure prediction sequence S is composed of the supply main pipe pressures of the k time intervals after the time t bi ; wherein trend b·t+k-1 is a set of the supply main pressure corresponding to k-1 time intervals after the first specified time t with a fixed index sequence [1, 2, 3...T] as the independent variable, and bi-t-k-1 S is obtained by using a preset fitting method as the dependent variable. pb t+k-1 are recursively derived using a second recursive formula; g is a model coefficient pre-acquired in the second recursive formula; dsteam b·t+k-1 = steam b·t+k-1 - steam b·t+k-2 ; dsteam b·t+k-1 is the flow rate variation of the steam supply main for k-1 time intervals after time t; steam b·t+k-1 Qs(t) is the flow rate in the steam supply main for k-1 time intervals after time t; The sequence S of prediction of the pressure of the steam main bi is: S bi = [pb t+1 , pb t+2 , pb t+3 ,..., pb t+k ].

5. The control method of the steam supply equipment in the thermal power production system according to claim 4, characterized in that, The model coefficient f in the first recursive formula is acquired in the following manner: From the beginning to the end of the second specified time period, a time window with a length of L×Δ is translated H times from the beginning to the end of the second specified time period in the direction of the beginning to the end, and each time the translation is Δ; The second specified time period is [t0-(H+L)Δ, t0]; Wherein, t0 is the second specified time; H is a pre-set number; L is a random integer between 6 and 12; After each translation of the time window, a set of main steam header steam flow at the time corresponding to the beginning of the time window and a set of main steam header pressure at the time corresponding to the end of the time window are acquired; For the main steam header steam flow at the time corresponding to the beginning of the time window and the main steam header steam flow at the time corresponding to the end of the time window acquired after the dth translation of the time window, the difference of the main steam header steam flow corresponding to the dth translation of the time window is acquired; wherein, 1≤d≤H; Taking the fixed index sequence [1, 2, 3...T] as the independent variable and the set of main steam header pressure at the time corresponding to the beginning of the time window acquired after the dth translation of the time window as the dependent variable, a pre-set fitting method is used to obtain the main steam header pressure trend at the time corresponding to the beginning of the time window acquired after the dth translation of the time window. The fixed index sequence [1, 2, 3...T] is taken as an independent variable, the set of main steam header pressure at the time corresponding to the end of the time window obtained after the dth translation of the time window is taken as a dependent variable, a preset fitting mode is adopted, and a main steam header pressure trend at the time corresponding to the end of the time window obtained after the dth translation of the time window is obtained; Based on the main steam header pressure trend at the time corresponding to the start of the time window obtained after the dth translation of the time window and the main steam header pressure trend at the time corresponding to the end of the time window obtained after the dth translation of the time window, a difference value of the main steam header pressure trend corresponding to the dth translation of the time window is obtained; According to the difference value of the main steam header steam flow corresponding to each translation of the time window and the difference value of the main steam header pressure trend, a first two-dimensional array is generated; The first two-dimensional array is: wherein delta_steam ad is the difference of the main steam header steam flow corresponding to the dth translation of the time window; delta_trend ad is the difference of the main steam header pressure trend corresponding to the dth translation of the time window; The difference value of the main steam header steam flow corresponding to each translation in the first two-dimensional array and the difference value of the main steam header pressure trend are respectively substituted into the first non-intercept linear regression model obtained in advance, the specific value of the regression model coefficient f0 in the first non-intercept linear regression model is trained, and the specific value of the regression model coefficient f0 in the first non-intercept linear regression model is taken as the model coefficient f in the first recursive formula. The first non-intercept linear regression model is: The difference value of the main steam header pressure trend = f0*the difference value of the main steam header steam flow.

6. The control method of the steam supply equipment in the thermal power generation system according to claim 5, wherein The model coefficient g in the second recursive formula is obtained in the following manner: From the start to the end of the second specified time period, a time window with a length of L*Delta is translated H times from the start to the end of the second specified time period in the direction, and each translation is Delta; After each translation of the time window, a set of steam header steam flow at the time corresponding to the start of the time window and a set of steam header pressure at the time corresponding to the start of the time window are obtained, and a set of steam header steam flow at the time corresponding to the end of the time window and a set of steam header pressure at the time corresponding to the end of the time window are obtained; For the steam header steam flow at the time corresponding to the start of the time window and the steam header steam flow at the time corresponding to the end of the time window obtained after the dth translation of the time window, a difference value of the steam header steam flow corresponding to the dth translation of the time window is obtained; The fixed index sequence [1, 2, 3...T] is taken as an independent variable, the set of main steam header pressure at the time corresponding to the end of the time window obtained after the dth translation of the time window is taken as a dependent variable, a preset fitting mode is adopted, and a main steam header pressure trend at the time corresponding to the end of the time window obtained after the dth translation of the time window is obtained; Taking the fixed index sequence [1, 2, 3…T] as the independent variable, taking the set of the steam supply main pipe pressure corresponding to the time point at the end of the time window obtained after the dth translation of the time window as the dependent variable, and adopting a pre-set fitting mode, a steam supply main pipe pressure trend corresponding to the time point at the end of the time window obtained after the dth translation of the time window is obtained; Based on the steam supply main pipe pressure trend corresponding to the time point at the start of the time window obtained after the dth translation of the time window and the steam supply main pipe pressure trend corresponding to the time point at the end of the time window obtained after the dth translation of the time window, a difference value of the steam supply main pipe pressure trend corresponding to the dth translation of the time window is obtained; According to the difference value of the steam supply main pipe steam flow corresponding to each translation of the time window and the difference value of the steam supply main pipe pressure trend, a second two-dimensional array is generated; Wherein the second two-dimensional array is: Wherein, delta-steambd is the difference value of the steam supply main pipe steam flow corresponding to the dth translation of the time window; delta_trendbd is the difference value of the steam supply main pipe pressure trend corresponding to the dth translation of the time window; The difference value of the steam supply main pipe steam flow corresponding to each translation in the second two-dimensional array and the difference value of the steam supply main pipe pressure trend are respectively substituted into the second non-intercept linear regression model obtained in advance, and the specific value of the regression model coefficient g0 in the second non-intercept linear regression model is obtained by training, and the specific value of the regression model coefficient g0 in the second non-intercept linear regression model is taken as the model coefficient g in the second recursive formula; Wherein, the second non-intercept linear regression model is: The difference value of the steam supply main pipe pressure trend = g0 x the difference value of the steam supply main pipe steam flow.

7. The control method of the steam supply equipment in the thermal power generation system according to claim 6, characterized in that, wherein The corresponding control solution space C in the thermal power generation system is obtained according to the set of adjustable values corresponding to each group of steam supply equipment groups in the thermal power generation system, and specifically includes: The Cartesian product of the set of adjustable values of each group of steam supply equipment groups in the thermal power generation system is taken as the final control solution space C; C = C1 x C2... x C i ... x C N ; C i a set of adjustable values corresponding to the i-th group of steam supply equipment in the thermal power production system; C i comprises M steam flow adjustment values preset in the main steam header corresponding to the i-th group of steam supply equipment groups; Wherein, each row of data in the final control solution space C is taken as a group of adjustment feasible solutions.

8. The control method of the steam supply device in a thermoelectric power production system according to claim 7, characterized by, The A4 specifically includes: A41、for any e-th group of main steam pipe inner steam flow adjustment feasible solution in the final control solution space C: [c1, c2,... c i ]; obtaining the main steam pipe inner steam flow adjustment solution corresponding to the group of main steam pipe inner steam flow adjustment feasible solution N ] and the main steam pipe inner steam flow adjustment solution corresponding to the group of main steam pipe inner steam flow adjustment feasible solution wherein wherein c i is the steam flow adjustment value for the main steam header corresponding to the i-th group of steam supply equipment in the thermal power production system; c i * is the steam flow adjustment value of the steam supply header corresponding to the i-th group of steam supply equipment in the thermal power production system; w i is the i-th c i corresponding coefficient; A42、the main steam header pressure corresponding to time t, the main steam header pressure trend at time t trend a-t and the flow adjustment value c in the main steam header i As the flow change amount in the main steam header, substitute into the first recursive formula to obtain the new main steam header pressure pa corresponding to time t+Δ t+1 *; wherein pa t+1 * = pa t + trend a·t + f*c i ; A43, the set S of main steam header pressures at time t ai-t and the new main steam header pressure pa corresponding to time t+Δ t+1 * using a third recursive formula to obtain the new main steam header pressures for the second to kth time intervals after time t, and to form the final main steam header pressure prediction sequence S from the new main steam header pressures for the k time intervals after time t ai *; Wherein, the third recursive formula is: wherein, is a set of new main steam header pressures corresponding to k-1 time intervals after the first specified time t, with a fixed index sequence [1, 2, 3...T] as the independent variable ai-t-k-1 is the dependent variable, obtained using a pre-set fitting method; At k = 1, is obtained from the first recurrence formula; when k > 1, is obtained from the third recurrence formula; Final main steam header pressure prediction sequence S ai * is: S ai *=[pa t+1 、pa t+2 、pa t+3 、…、pa t+k ]; A44, the supply main pressure corresponding to time t, the supply main pressure trend trend at time t a·t and the flow adjustment value c in the supply main pipe i * as the flow change amount in the supply main pipe, substituted into the second recursive formula to obtain the new supply main pressure corresponding to time t+Δ wherein, A45, a set of steam main pressures at time t bi-t a new steam main pressure corresponding to time t+Δ using a fourth recursive formula, a new steam main pressure for each of the second to kth time intervals after time t is obtained, and the final steam main pressure prediction sequence S is composed of the new steam main pressures for the k time intervals after time t bi *; Wherein, the fourth recursive formula is: wherein, is a set of new main steam line pressures corresponding to k-1 time intervals after the first specified time t with a fixed index sequence [1, 2, 3...T] as the independent variable bi-t-k-1 is the dependent variable, obtained using a pre-set fitting method; At k = 1, is obtained from the second recurrence formula; when k > 1, is obtained from the fourth recurrence formula; Final supply main pressure prediction sequence S bi * is: S bi * = [pb t+1 , pb t+2 , pb t+3 ,..., pb t+k ]; A46. For any each group of main steam header steam flow adjustment feasible solution in the final control solution space C, the final main steam header pressure prediction sequence S corresponding to the i-th group of steam supply equipment group ai *and the final main steam header pressure prediction sequence S bi *For screening, the optimal recommended value is determined.

9. The control method of the steam supply device in a thermoelectric power production system according to claim 8, characterized by, The A46 specifically includes: A46-1, for any each group of main steam main pipe steam flow adjustment feasible solutions in the final control solution space C, the final main steam main pipe pressure prediction sequence and the final steam supply main pipe pressure prediction sequence corresponding to the i-th group of steam supply equipment groups are screened, and the final main steam main pipe pressure prediction sequence and the final steam supply main pipe pressure prediction sequence corresponding to the i-th group of steam supply equipment groups that meet the pre-set target condition are obtained; Wherein, the target condition is: Any new value of the main steam header pressure in the final main steam header pressure prediction sequence is greater than or equal to a first preset value lower corresponding to the ith group of steam supply equipment i1 , and less than or equal to a second preset value upper corresponding to the ith group of steam supply equipment i1 ; Any new steam main pressure value in the final steam main pressure prediction sequence is greater than or equal to the third preset value lower corresponding to the i-th group of steam supply equipment groups i2 , and less than or equal to the fourth preset value upper corresponding to the i-th group of steam supply equipment groups i2 ; A46-2, for the final main steam header pressure prediction sequence corresponding to the i-th group of steam supply equipment group meeting the preset target condition, the main steam header pressure convergence factor μ corresponding to the i-th group of steam supply equipment group is obtained according to formula (1) i1 ; The formula (1) is: According to formula (2), the steam supply main pipe pressure convergence factor μ corresponding to the i-th group of steam supply equipment groups is obtained for the final steam supply main pipe pressure prediction sequence corresponding to the i-th group of steam supply equipment groups that meet the preset target conditions i2 ; The formula (2) is: A46-3, according to the main steam main pipe pressure convergence factor and the steam supply main pipe pressure convergence factor corresponding to each group of steam supply equipment groups that meet the pre-set target condition, the formula (3) is adopted to obtain the first pressure convergence value μ corresponding to the thermal power generation system; wherein the formula (3) is: A46-4, the steam flow adjustment value in the main steam header and the steam flow adjustment value in the steam supply header corresponding to the steam flow adjustment value in the main steam header are taken as the optimal adjustment value of the main steam header and the optimal adjustment value of the steam supply header, respectively.

Images