Method for Optimizing Operation of Combined Cycle Gas Turbine System

Abstract

The present disclosure provides a method for optimizing operation of a combined cycle gas turbine system, which includes the following steps: firstly, building a process flow model of a gas-fired power generation system as well as a process flow model of a steam power generation system; then, determining energy efficiency indexes, an environmental evaluation index, and thermoeconomic evaluation indexes of the system; next, building an overall evaluation model by analyzing, through an entropy weight method, weight indexes such as a primary energy ratio, exergy efficiency, a per-unit emission amount of CO.sub.2, and a per-unit thermoeconomic cost of the system; and finally, building an optimization model by means of particle swarm optimization.

Claims

1. A method for optimizing operation of a combined cycle gas turbine system, comprising the following steps: S1, building a process flow model of a gas-fired power generation system as well as a process flow model of a steam power generation system; S2, determining energy efficiency indexes and an environmental evaluation index of a combined cycle gas turbine system, wherein a primary energy ratio and exergy efficiency of the system are served as the energy efficiency indexes of the system, and mass of CO.sub.2 emitted by the system to generate per-unit electricity is served as the environmental evaluation index; S3, determining thermoeconomic evaluation indexes of the combined cycle gas turbine system; S4, building an overall evaluation model by analyzing, through an entropy weight method, weight indexes such as the primary energy ratio, the exergy efficiency, a per-unit emission amount of the CO.sub.2, and a per-unit thermoeconomic cost of the system; particularly: S41, normalization of the indexes firstly, totally numbering m operating conditions, participating in evaluation, of the system as M, wherein M=(m.sub.1, m.sub.2, m.sub.3 m.sub.m); totally numbering n evaluation indexes of the system as D, wherein D=(d.sub.1, d.sub.2, d.sub.3 d.sub.n); and recording a value of the i.sup.th evaluation index of the evaluated operating condition m.sub.i as x.sub.ij to form an evaluation index matrix X=[x.sub.ij].sub.m*n composed of m*n indexes; $\begin{array}{cc}X=\left[\begin{array}{cccc}{x}_{11}& x_{12}& L& x_{1n}\\ x_{21}& x_{22}& L& x_{2n}\\ M& M& O& M\\ x_{m1}& x_{m2}& L& x_{mn}\end{array}\right]& \left(4-1\right)\end{array}$ then, normalizing the indexes based on types of the indexes, wherein the indexes expressing performance improved with an increase in values of evaluation results are normalized according to formula (4-2), and indexes expressing the performance improved with a decrease in the values of the evaluation results are normalized according to formula (4-3); $\begin{array}{cc}{V}_{\mathrm{ij}}=\frac{x_{\mathrm{ij}}-\min \left(x_j\right)}{\max \left(x_j\right)-\min \left(x_j\right)}& \left(4-2\right)\end{array}$ $\begin{array}{cc}{V}_{\mathrm{ij}}=\frac{\max \left(x_j\right)-x_{\mathrm{ij}}}{\max \left(x_j\right)-\min \left(x_j\right)}& \left(4-3\right)\end{array}$ in formula (4-2) and formula (4-3), min(x.sub.j) represents the minimum value of the j.sup.th evaluation index under the operating conditions; and max(x.sub.j) represents the maximum value of the j.sup.th evaluation index under the operating conditions; and finally, calculating a proportion of features of the i.sup.th load condition in the presence of the j.sup.th evaluation index to form a normalized matrix P expressed by formula (4-4); $\begin{array}{cc}{P}_{ij}=\frac{V_{ij}}{\stackrel{m}{\underset{i=1}{.Math.}}{V}_{ij}}& \left(4-4\right)\end{array}$ in formula (4-4), V.sub.ij represents a value of a normalized and dimensionless index x.sub.ij; and P.sub.ij represents the proportion of the features; S42, information entropy calculation on the indexes working out a value of information entropy corresponding to the j.sup.th evaluation index according to formula ((4-5); $\begin{array}{cc}{e}_j=-1/\ln \left(m\right)\underset{i=1}{\overset{m}{.Math.}}{p}_{\mathrm{ij}}.Math.\ln p_{\mathrm{ij}}& \left(4-5\right)\end{array}$ in formula (4-5), e.sub.j represents the value of the information entropy of the i.sup.th evaluation index; and P.sub.ij represents the proportion of the features; S43, weight calculation on the indexes working out a difference coefficient of the evaluation index X.sub.j according to formula (4-6), and working out an entropy weight w.sub.j of the j.sup.th evaluation index according to formula (4-7): $\begin{array}{cc}{d}_j=1-e_j& \left(4-6\right)\end{array}$ $\begin{array}{cc}{w}_j=\frac{d_j}{\underset{j=1}{\overset{n}{.Math.}}{d}_j}& \left(4-7\right)\end{array}$ in formula (4-6) and formula (4-7), d.sub.j represents a difference of the j.sup.th evaluation index; and w.sub.j represents a weight ratio of the j.sup.th evaluation index; S44, calculation on overall evaluation indexes wherein, an overall effectiveness evaluation index K.sub.i under the i.sup.th operating condition is as follows: $\begin{array}{cc}{K}_i=\underset{j=1}{\overset{n}{.Math.}}{w}_j{V}_{\mathrm{ij}}& \left(4-8\right)\end{array}$ in formula (4-8), V.sub.ij represents a value of a normalized and dimensionless index x.sub.ij; and S5, building an optimization model by means of particle swarm optimization.

2. The method for optimizing operation of a combined cycle gas turbine system according to claim 1, wherein step S5 particularly comprises: setting the overall evaluation model as an optimization objective; establishing constraint conditions of the system; establishing an adaptive function group; and after the optimization objective, the constraint conditions, and the adaptive function group are determined, building the operation optimization model according to a calculation process of the particle swarm optimization.

3. The method for optimizing operation of a combined cycle gas turbine system according to claim 2, wherein in the adaptive function group, independent variables include inlet guide vane (IGV) opening to be optimized, natural gas flow, and a natural gas price having an influence on the per-unit thermoeconomic cost of the system; and dependent variables include the primary energy ratio, the exergy efficiency, the per-unit emission amount of the CO.sub.2, and the per-unit thermoeconomic cost which are related to the optimization objective, as well as an operating load of the system and an outlet flue gas temperature of a gas turbine, which are related to the constraint conditions.

4. The method for optimizing operation of a combined cycle gas turbine system according to claim 1, wherein in step S1, the process flow model of the gas-fired power generation system as well as the process flow model of the steam power generation system are built based on an actual production process of the combined cycle gas turbine system by means of process simulation software, namely Aspen Plus, and thermodynamic models of devices of the combined cycle gas turbine system.

5. The method for optimizing operation of a combined cycle gas turbine system according to claim 1, wherein in step S2, a primary energy ratio index is established by analyzing, based on energy analysis, an energy balance among a gas turbine system, a waste heat boiler system, and a steam turbine system; an exergy efficiency index is established by analyzing, based on energy analysis, an exergy balance among main devices of the system; and components of a flue gas from the system is analyzed, and the mass of the CO.sub.2 emitted by the system to generate the per-unit electricity is served as the environmental evaluation index.

6. The method for optimizing operation of a combined cycle gas turbine system according to claim 1, wherein in step S3, thermoeconomic models are built through the following steps: S31, drawing a productive structure diagram of the system according to a productive consumption relationship between fuels and the devices of the system and between products and the devices of the system; S32, building fuel-product calculation models of the devices of the system, to determine the fuels and the products; and S33, building the thermoeconomic models of the devices of the system to analyze a thermoeconomic cost of the system.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0062] FIG. 1 shows a flow chart of a method for optimizing operation of a combined cycle gas turbine system of the present disclosure;

[0063] FIG. 2 shows a process flow model of a gas-fired power generation system in an embodiment of the present disclosure;

[0064] Reverence numerals: COMPRESS-air compressor; COMBUST-combustion chamber; TURBINE-turbine;

[0065] FIG. 3 shows a process flow model of a steam power generation system in the embodiment of the present disclosure;

[0066] Reference numerals: RHEAT2—intermediate-pressure reheater 2; HSUP2—high-pressure superheater 2; RHEAT1—intermediate-pressure reheater 1; HSUP1—high-pressure superheater 1; HVAPOR—high-pressure evaporator; HECONOMI—high-pressure economizer; MSUP—intermediate-pressure superheater; MVAPOR—intermediate-pressure evaporator; MECONO MI—intermediate-pressure economizer; LSUP—low-pressure superheater; LVAPOR—low-pressure evaporator; HEAT—feedwater heater; HDRUM—high-pressure steam drum; IDRUM—intermediate-pressure steam drum; LDRUM—low-pressure steam drum; HPC—high-pressure cylinder of a steam turbine; IPC—intermediate-pressure cylinder of the steam turbine; LPC—low-pressure cylinder of the steam turbine; COND—condenser; CPUMP—condensate pump; IPUMP—intermediate-pressure water-delivery pump; HPUMP—high-pressure water-delivery pump;

[0067] FIG. 4 shows a productive structure diagram a combined cycle gas turbine system in a city in the embodiment of the present disclosure;

[0068] FIG. 5 shows per-unit thermoeconomic costs of main productive devices;

[0069] FIG. 6 shows composition of the per-unit thermoeconomic costs of the main productive devices;

[0070] FIG. 7 shows a change of air flow along with that of IGV opening of the combined cycle gas turbine system;

[0071] FIG. 8 shows a calculation process of particle swarm optimization;

[0072] FIG. 9 shows overall evaluation results of the system before and after optimization in spring;

[0073] FIG. 10 shows overall evaluation results of the system before and after the optimization in summer;

[0074] FIG. 11 shows overall evaluation results of the system before and after the optimization in autumn; and

[0075] FIG. 12 shows overall evaluation results of the system before and after the optimization in winter.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0076] The preferred embodiments of the present disclosure are described below with reference to the drawings. It should be understood that the preferred embodiments described herein are only used to illustrate the present disclosure, rather than to limit the present disclosure.

[0077] A method for optimizing operation of a combined cycle gas turbine system of the present disclosure is explained in detail with a combined cycle gas turbine system in a city as an example. As shown in FIG. 1 to FIG. 12, the method includes steps S1-S5.

[0078] S1, Build process flow models;

[0079] Particularly, build the process flow models of the combined cycle gas turbine system in the city by means of Aspen Plus; Where, the process flow model of a gas-fired power generation system is shown in FIG. 2; and the process flow model of a steam power generation system is shown in FIG. 3.

[0080] S2, Determine energy efficiency indexes and an environmental evaluation index of the combined cycle gas turbine system;

[0081] Particularly, establish a primary energy ratio index by analyzing, based on energy analysis, an energy balance among a gas turbine system, a waste heat boiler system, and a steam turbine system;

[0082] Where, the primary energy ratio of the combined cycle gas turbine system is expressed as follows:

[00010]$\begin{array}{cc}{\eta }_Q=1-\frac{.Math.Q_{\mathrm{si}}}{Q_{\mathrm{fuel}}}=\frac{W_1+W_2}{Q_{\mathrm{fuel}}};& \left(2-1\right)\end{array}$

[0083] In formula (2-1), Q.sub.si represents an energy loss of each part, which is measured in kJ/s;

[0084] Q.sub.fuel represents a lower heating value of a fuel entering the gas turbine system, which is measured in kJ/s;

[0085] W.sub.1 represents electric energy generated by the gas turbine system, which is measured in kJ/s; and

[0086] W.sub.2 represents electric energy generated by the steam turbine system, which is measured in kJ/s;

[0087] establishing an exergy efficiency index by analyzing, based on exergy analysis, an exergy balance among main devices of the combined cycle gas turbine system;

[00011]$\begin{array}{cc}{\eta }_{\mathrm{Ex}}=\frac{\underset{x}{\overset{n}{.Math.}}{E}_{\mathrm{in},x}-.Math.I}{\underset{x}{\overset{n}{.Math.}}{E}_{\mathrm{in},x}};& \left(2-2\right)\end{array}$

[0088] In formula (2-2), E.sub.in,x represents a value of an exergy flow entering the system, which is measured in kJ/s; and

[0089] I represents an exergy loss of the system, which is measured in kJ/s;

[0090] Where, the primary energy ratio and exergy efficiency of the system are served as the energy efficiency indexes of the system;

[0091] Analyze components of a flue gas from the system, where the mass of CO.sub.2 emitted by the system to generate per-unit electricity is served as the environmental evaluation index;

[00012]$\begin{array}{cc}{\lambda }_{{\mathrm{CO}}_2}=\frac{m_{{\mathrm{CO}}_2}}{W_e};& \left(2-3\right)\end{array}$$\begin{array}{cc}{m}_{{\mathrm{CO}}_2}=m_{gas}.Math.\frac{M_{{\mathrm{CO}}_2}}{M_{gas}};& \left(2-4\right)\end{array}$

[0092] In formula (2-3) and formula (2-4), λ.sub.CO2 represents an amount of the CO.sub.2 emitted by the system to generate the per-unit electricity, which is measured in g/(kW.Math.h);

[0093] M.sub.CO2 represents an amount of CO.sub.2 in the flue gas, which is measured in g/kg; and

[0094] M.sub.CO2 represents molar mass of the CO.sub.2 in the flue gas, which is measured in kg/mol; and M.sub.gas represents molar mass the flue gas, which is measured in kg/mol.

[0095] The primary energy ratio of the combined cycle gas turbine system in the city is 55.56%. The exergy efficiency of the system is 52.84%, and the amount of the CO.sub.2 emitted by the system to generate the per-unit electricity is calculated as 1287.31 g/(kW.Math.h).

[0096] S3, Determine thermoeconomic evaluation indexes of the combined cycle gas turbine system;

[0097] (1) draw a productive structure diagram, as shown in FIG. 4, of the system according to a productive consumption relationship between fuels and the devices of the system and between products and the devices of the system;

[0098] (2) Build fuel-product calculation models of the devices of the system, as shown in table 1;

TABLE-US-00001 TABLE 1 Fuel-product calculation models Device in the system Fuel Product Gas turbine system Combustion chamber FB = E.sub.1 PB = E.sub.3 − E.sub.4 FS = T.sub.0(S.sub.3 − S.sub.4) Air compressor FB = E.sub.23 PB = E.sub.3 − E.sub.2 FS = T.sub.0(S.sub.3 − S.sub.2) Turbine FB = E.sub.4 − E.sub.5 PB = E.sub.22 + E.sub.23 FS = T.sub.0(S.sub.4 − S.sub.5) Waste heat boiler system FB = E.sub.5 − E.sub.6 PB = E.sub.13 + E.sub.15 + T.sub.17 + E.sub.7 − FS = T.sub.0(S.sub.13 + S.sub.15 + S.sub.17 + S.sub.7 − E.sub.8 − E.sub.11 − E.sub.12 − E.sub.16 − E.sub.21 S.sub.8 − S.sub.11 − S.sub.12 − S.sub.16 − S.sub.21) PS = T.sub.0 (S.sub.5 − S.sub.6) Steam turbine system High-pressure cylinder FB = E.sub.15 − E.sub.16 PB = E.sub.24 FS = T.sub.0(S.sub.15 − S.sub.16) Intermediate-pressure cylinder FB = E.sub.13 − E.sub.14 PB = E.sub.25 FS = T.sub.0(S.sub.13 − S.sub.14) Low-pressure cylinder FB = E.sub.18 − S.sub.14 PB = E.sub.26 FS = T.sub.0(S.sub.18 − S.sub.19) Condenser FB = E.sub.19 − E.sub.20 FS = T.sub.0(S.sub.19 − S.sub.20) Pump system Low-pressure pump FB = E.sub.27 PB = E.sub.21 − E.sub.20 FS = T.sub.0(S.sub.21 − S.sub.20) Intermediate-pressure pump FB = E.sub.29 PB = E.sub.11 − E.sub.9 FS = T.sub.0(S.sub.11 − S.sub.9) High-pressure pump FB = E.sub.28 PB = E.sub.12 − E.sub.10 FS = T.sub.0(S.sub.12 − S.sub.10) Electric generator FB = E.sub.22 + E.sub.24 + E.sub.25 + E.sub.26 PB = E.sub.30 Chimney FB = E.sub.6 FS = T.sub.0S.sub.6

[0099] (3) Build thermoeconomic models (as shown in table 2) of the devices of the combined cycle gas turbine system to analyze a thermoeconomic cost (as shown in FIG. 5) of the system; and

TABLE-US-00002 TABLE 2 Thermoeconomic models of the devices of the combined cycle gas turbine system Device in the system Thermoeconomic model Gas turbine system Combustion chamber PB.sub.1 .Math. C.sub.PB, 1 = FB.sub.1 .Math. C.sub.FB, 1 + FS.sub.1 .Math. C.sub.FS, 1 + Z.sub.1 Air compressor PB.sub.2 .Math. C.sub.PB, 2 = FB.sub.2 .Math. C.sub.FB, 2 + FS.sub.2 .Math. C.sub.FS, 2 + Z.sub.2 Turbine PB.sub.3 .Math. C.sub.PB, 3 = FB.sub.3 .Math. C.sub.FB, 3 + FS.sub.3 .Math. C.sub.FS, 3 + Z.sub.3 Waste heat boiler system PB.sub.4 .Math. C.sub.PB, 4 + PS.sub.4 .Math. C.sub.PS, 4 = FB.sub.4 .Math. C.sub.FB, 4 + FS.sub.4 .Math. C.sub.FS, 4 + Z.sub.4 Steam turbine system High-pressure cylinder PB.sub.5 .Math. C.sub.PB, 5 = FB.sub.5 .Math. C.sub.FB, 5 + FS.sub.5 .Math. C.sub.FS, 5 + Z.sub.5 Intermediate-pressure cylinder PB.sub.6 .Math. C.sub.PB, 6 = FB.sub.6 .Math. C.sub.FB, 6 + FS.sub.6 .Math. C.sub.FS, 6 + Z.sub.6 Low-pressure cylinder PB.sub.7 .Math. C.sub.PB, 7 = FB.sub.7 .Math. C.sub.FB, 7 + FS.sub.7 .Math. C.sub.FS, 7 + Z.sub.7 Condenser PS.sub.12 .Math. C.sub.PS, 12 = FB.sub.12 .Math. C.sub.FB, 12 + Z.sub.12 Pump system Low-pressure pump PB.sub.9 .Math. C.sub.PB, 9 = FB.sub.9 .Math. C.sub.FB, 9 + F.sub.9 .Math. C.sub.FS, 9 + Z.sub.9 Intermediate-pressure pump PB.sub.10 .Math. C.sub.PB, 10 = FB.sub.10 .Math. C.sub.FB, 10 + FS.sub.10 .Math. C.sub.FS, 10 + Z.sub.10 High-pressure pump PB.sub.11 .Math. C.sub.PB, 11 = FB.sub.11 .Math. C.sub.FB, 11 + FS.sub.11 .Math. C.sub.FS, 11 + Z.sub.11 Electric generator PB.sub.8 .Math. C.sub.PB, 8 = FB.sub.8 .Math. C.sub.FB, 8 + Z.sub.8 Flue gas PS.sub.13 .Math. C.sub.PS, 13 = FB.sub.13 .Math. C.sub.FB, 13 + Z.sub.13 J1 C.sub.PB, 14 = Σr.sub.i .Math. C.sub.PB, i(i = 1, 2) J2 C.sub.PB, 15 = Σr.sub.i .Math. C.sub.PB, i(i = 4, 9, 10, 11) J3 C.sub.FB, 8 = Σr.sub.i .Math. C.sub.PB, i(i = 5, 6, 7, 16) J4 C.sub.PS, 17 = Σr.sub.i .Math. C.sub.PS, i(i = 4, 12, 13) B1 C.sub.FB, j = C.sub.PB, 14(j = 3, 4, 13) B2 C.sub.PB, j = C.sub.FB, 2(j = 3, 16) B3 C.sub.FB, j = C.sub.PB, 15(j = 5, 6, 7, 12) B4 C.sub.FB, j = C.sub.PB, 8(j = 9, 10, 11) B5 C.sub.FS, j = C.sub.PS, 17(j = 1, 2, 3, 4, 5, 6, 7, 9, 10, 11)

[0100] (4) Through analysis on composition (as shown in FIG. 6) of the thermoeconomic cost of the system, analyze, based on operating parameters of the system under a basic operating condition, the thermoeconomic cost of the system under the basic operating condition by means of the thermoeconomic models to evaluate the economy of the system.

[0101] With respect to the combined cycle gas turbine system in the city, the low-pressure cylinder has the highest per-unit thermoeconomic cost of 0.5567 yuan/(kW.Math.h); and the combustion chamber has the lowest per-unit thermoeconomic cost of 0.2714 yuan/(kW.Math.h). A product of the electric generator is equivalent to the electric energy generated by the system. Therefore, the per-unit thermoeconomic cost of 0.4848 yuan/(kW.Math.h) is equivalent to the per-unit power generation cost of the system.

[0102] S4, Build an overall evaluation model of the combined cycle gas turbine system;

[0103] Where, a method for building the overall evaluation model is put forward to overall evaluate the energy efficiency, environmental friendliness, and economy of the system; Particularly, build the overall evaluation model by analyzing, through an entropy weight method, weight indexes such as the primary energy ratio, the exergy efficiency, the per-unit emission amount of the CO.sub.2, and the per-unit thermoeconomic cost of the system; where, detailed steps are as follows:

[0104] (1) Normalization of the indexes, as shown in table 3;

TABLE-US-00003 TABLE 3 Normalization results of a feature proportion matrix Per-unit Thermal Exergy emission Per-unit Operating efficiency efficiency amount P.sub.i3 thermoeconomic condition P.sub.i1 P.sub.i2 of the CO2 cost P.sub.i4 m.sub.1 0.0545 0.0553 0.0546 0.0821 m.sub.2 0.0342 0.0379 0.0347 0.0653 m.sub.3 0.0327 0.0361 0.0333 0.0571 m.sub.4 0.0125 0.0188 0.0306 0.0401 m.sub.5 0.0122 0.0178 0.0126 0.0231 m.sub.6 0.0034 0.0027 0.0032 0.0007 m.sub.7 0.0519 0.0540 0.0520 0.0812 m.sub.8 0.0316 0.0341 0.0323 0.0641 m.sub.9 0.0275 0.0300 0.0308 0.0553 m.sub.10 0.0093 0.0145 0.0305 0.0389 m.sub.11 0.0029 0.0099 0.0097 0.0207 m.sub.12 0.0000 0.0000 0.0000 0.0000 m.sub.13 0.1395 0.1181 0.1189 0.0869 m.sub.14 0.0918 0.0885 0.0893 0.0641 m.sub.15 0.0703 0.0709 0.0688 0.0481 m.sub.16 0.0560 0.0580 0.0561 0.0339 m.sub.17 0.1442 0.1237 0.1229 0.0884 m.sub.18 0.0941 0.0921 0.0915 0.0653 m.sub.19 0.0726 0.0750 0.0694 0.0496 m.sub.20 0.0589 0.0623 0.0587 0.0353

[0105] (2) Information entropy calculation on the indexes based on formula (4-5), where calculation results are shown in table 4;

[0106] (3) Weight calculation on the indexes based on formula (4-6) and formula (4-7), where calculation results are shown in table 4; and

TABLE-US-00004 TABLE 4 Calculation results of the entropy weight method Per-unit Primary Exergy emission Per-unit energy efficiency amount V.sub.3 thermoeconomic ratio V.sub.1 V.sub.2 of CO.sub.2 cost V.sub.4 Information entropy e.sub.j 0.3793 0.3680 0.3656 0.3514 Difference d.sub.j 0.6207 0.6320 0.63442 0.6486 Weight w.sub.j 0.2448 0.2492 0.2502 0.2558

[0107] (4) Weight calculation on the indexes;

[0108] Particularly, substitute the weight of each evaluation index into formula (4-8) to build the following overall effectiveness evaluation model of the combined cycle gas turbine system in the city;

[00013]$\begin{array}{cc}\begin{array}{c}{K}_i=0.2448V_{i1}+0.2492V_{i2}+0.2502V_{i3}+0.2558V_{i4}\\ =0.2448\times \frac{{\eta }_Q-50.18}{5.56}+0.2492\times \frac{{\eta }_{Ex}-47.61}{5.48}+0.2502\times \frac{1405.16-m_{{\mathrm{CO}}_2}}{127.48}\\ +0.2558\times \frac{0.6417-C_8}{0.1076}\end{array}.& \left(4-9\right)\end{array}$

[0109] S5, Build an optimization model by means of particle swarm optimization;

[0110] Particularly, in order to obtain the optimal operating parameters in a case of variable loads of the system, build, by means of the particle swarm optimization, the optimization model of the system with IGV opening of an air compressor and natural gas flow as variables to obtain the highest overall evaluation of the system. A change of air flow along with that of the IGV opening is shown in FIG. 7. A process of the particle swarm optimization is shown in FIG. 8.

The overall evaluation model is particularly built through the following steps:

[0111] (1) Optimization Objective

[0112] In order to improve the primary energy ratio and exergy efficiency of the system and reduce the per-unit emission amount of the CO.sub.2 and per-unit thermoeconomic cost of the system under different operating conditions, set the built overall evaluation model (4-9) of the combined cycle gas turbine system as the optimization objective;

[0113] (2) Constraint Conditions

[0114] In order to guarantee safe operation of the system and satisfy the demand of users for electric loads, establish the constraint conditions of the system, which are expressed by formula 5-1; where, with a combined cycle gas turbine system in Dazhou as an example, an outlet flue gas temperature of a gas turbine should not be higher than 600° C., and the IGV opening ranges from 12% to 98%; only in this case, the system can operate safely; and in order to make the electricity generated by the system be adequate for the electric loads needed by the users, a power generation load of the system is equalized to a needed power generation load;

[00014]$\begin{array}{cc}s.t.\{\begin{array}{c}\max \left[K_i\right]\\ T_6⩽T_{6,\max }\\ {\alpha }_{\min }⩽\alpha ⩽{\alpha }_{\max }\\ {\mathrm{Load}}_e={\mathrm{Load}}_{\mathrm{need}}\end{array}& \left(5-1\right)\end{array}$

[0115] In formula (5-1), T.sub.6,max represents the maximum allowable outlet flue gas temperature, namely 600° C., of the gas turbine;

[0116] α.sub.min represents the minimum value, namely 12%, of the IGV opening, and α.sub.max represents the maximum value, namely 98%, of the IGV opening; and

[0117] Laod.sub.e represents the power generation load of the system, and Load.sub.need represents the needed power generation load;

[0118] (3) Establishment of an Adaptive Function Group

[0119] In the particle swarm optimization, independent variables are required to be substituted into the adaptive function group to determine current “positions” of particles. The independent variables of the adaptive function group include the IGV opening to be optimized, the natural gas flow, and a natural gas price having an influence on the per-unit thermoeconomic cost of the system. Dependent variables include the primary energy ratio, the exergy efficiency, the per-unit emission amount of the CO.sub.2, and the per-unit thermoeconomic cost which are related to the optimization objective, as well as the operating load of the system and the outlet flue gas temperature of the gas turbine, which are related to the constraint conditions.

[0120] The adaptive function group is particularly established through the following steps:

[0121] Firstly, simulate, by means of the process flow model, operating conditions of the combined cycle gas turbine system in Dazhou in a case where the IGV opening ranges from 12% to 98% and the natural gas flow ranges from 8.16 kg/s to 12.95 kg/s, and calculate the primary energy ratio, the exergy efficiency, the per-unit emission amount of the CO.sub.2, the per-unit thermoeconomic cost, the operating load, and the outlet flue gas temperature of the gas turbine of the system under the operating conditions;

[0122] Then, according to simulation results, perform fitting, by means of a fitting analysis tool of a matrix laboratory (matlab), on the operating load (f.sub.Load), the primary energy ratio (f.sub.Q), the exergy efficiency (f.sub.Exergy), the per-unit emission amount (f.sub.CO2) of the CO.sub.2, the per-unit thermoeconomic cost (f.sub.Cost), and the outlet flue gas temperature (f.sub.T) of the gas turbine to obtain the adaptive function group related to the IGV opening (x), the natural gas flow (y), and the natural gas price (m); and

[0123] After the optimization objective, the constraint conditions, and the adaptive function group are determined, build the operation optimization model according to a calculation process of the particle swarm optimization by writing calculation codes through the matlab.

[0124] Where, the adaptive function group f.sub.i of the combined cycle gas turbine system in spring, summer, autumn, and winter is respectively denoted by f.sub.1, f.sub.2, f.sub.3, and f.sub.4.

[0125] The adaptive function group f.sub.1 of the combined cycle gas turbine system in the city in spring is expressed by formula (5-2) to formula (5-7).

f(x,y).sub.Load,1=−39.82−1.299x+16.71y−3.76×10.sup.−3x.sup.2+1.352×10.sup.−1xy−5.432×10.sup.−1y.sup.2   (5-2)

f(x,y).sub.Q,1=−11.14−0.1561x+4.156y−4.172x.sup.2+3.559×10.sup.−2xy−4.933y.sup.2+4.894×10.sup.−5x.sup.2y−2.055×10.sup.−3xy.sup.2+1.953×10.sup.−2y.sup.3   (5-3)

f(x,y).sub.Exergy,1=−11.18−0.1519x+4.176y−3.759×10.sup.−4x.sup.2+0.03502xy−0.4983y.sup.2−2.055×10.sup.−3xy.sup.2+1.988×10.sup.−2y.sup.3   (5-4)

f(x,y).sub.cO.sub.2.sub.,1=7.322×10.sup.4+869.3x−2.554×10.sup.4y+2.163x.sup.2−201.2xy+3023y.sup.2−0.266x.sup.2y+11.77xy.sup.2−119.3y.sup.3   (5-5)

f(x,y,m).sub.Cost,1=7.781×10.sup.−4x−4.6×10.sup.−2y+2.94×10.sup.−3m+0.98307   (5-6)

f(x,y).sub.T,1=−747.1+11.37x+283.2y+0.5083x.sup.2−6.306xy−13.81y.sup.2+1.353×10.sup.−3x.sup.3−5.5923×10.sup.−2x.sup.2y+0.4919xy.sup.2   (5-7)

[0126] The adaptive function group f.sub.2 of the combined cycle gas turbine system in the city in summer is expressed by formula (5-8) to formula (5-13).

f(x,y).sub.Load,2=−4.629−0.791x+8.157y−2.646×10.sup.−3x.sup.2+8.492×10.sup.−2xy−3.372×10.sup.−2y.sup.2   (5-8)

f(x,y).sub.Q,2=0.4564−4.809×10.sup.3x+7.441×10.sup.−3y−1.372×10.sup.−5×.sup.2+4.736×10.sup.−4xy   (5-9)

f(x,y).sub.Exergy,2=3.344−0.256x+0.7119y−3.362×10.sup.−3x.sup.2+0.06578xy−0.04461y.sup.2−4.452×10.sup.−4x.sup.2y−4.222×10.sup.−3xy.sup.2   (5-10)

f(x,y).sub.CO.sub.2.sub.,2=1.612×10.sup.4+98.33x−4862y−19.3xy+532.4y.sup.2−19.45y.sup.3+0.9695xy.sup.2   (5-11)

f(x,y,m).sub.Cost,2=8.8905×10.sup.−4x−4.789×10.sup.−2y+3.12×10.sup.−3m+0.99539   (5-12)

f(x,y).sub.T,2=−4300+67.76x+1675y+0.009236x.sup.2+14.13xy−195.4y.sup.2+0.0119x.sup.2y−0.8473xy.sup.2+7.923y.sup.3   (5-13)

[0127] The adaptive function group f.sub.3 of the combined cycle gas turbine system in the city in autumn is expressed by formula (5-14) to formula (5-19).

f(X,y).sub.Load,3=−66.1−1.642x+23.3y−4.866×10.sup.−3x.sup.2+0.1824xy−0.9448y.sup.2   (5-14)

f(x,y).sub.Q,3=−1.081−0.01798x+0.5048y+2.988×10.sup.−3xy−0.05288y.sup.2+1.294×10.sup.−4xy.sup.2+1.873×10.sup.−3y.sup.2   (5-15)

f(x,y).sub.Energy,3=−1.79−0.02354x+0.7238y+4.063×10.sup.−3xy0.07589y.sup.2−1.809×10.sup.−4xy.sup.2−2.669×10.sup.−3y.sup.3   (5-16)

f(x,y).sub.cO.sub.2.sub.,3=1330+19.58x+2.957y−0.05214x.sup.2−2.204xy−6.154×10.sup.−4x.sup.3+0.01502x.sup.2y   (5-17)

f(x,y,m).sub.cost,3=8.6461×10.sup.−4x−4.697×10.sup.−2y+3.12×10.sup.−3m+0.9854   (5-18)

f(x,y).sub.T,3=149.2−9.069x+5957y+0.1118x.sup.2+0.2542xy+0.007161x.sup.2y   (5-19)

[0128] The adaptive function group f.sub.4 of the combined cycle gas turbine system in the city in winter is expressed by formula (5-20) to formula (5-25).

f(x,y).sub.Load,4234.4−2.574x+82.17y+0.4274xy−7.924y.sup.2−0.01855xy.sup.2+0.2834y.sup.3   (5-20)

f(x,y).sub.Q,4=−0.9701—0.04801x+0.3477y+8.818×10.sup.−4x.sup.2+0.01463xy−0.0198y.sup.2+1.017×10.sup.−3Xy.sup.2+1.075×10.sup.−4x.sup.2y   (5-21)

f(x,y).sub.Exergy,4=−1.032−0.0612x+0.3541y+1.07x−0.01786xy−0.01995y.sup.2−1.283×10.sup.−4x.sup.2y+1.206×10.sup.−3xy.sup.2   (5-22)

f(x,y).sub.CO.sub.2.sub.,4=13490+105.8x−3862y−18.77xy+404.1y.sup.2+0.8465xy.sup.2−14.02y.sup.3   (5-23)

f(x,y,m).sub.Cost,4=7.37246×10.sup.−4x−0.04496y+0.0031m+0.96814   (5-24)

f(x,y).sub.T,4=−823.3+25.74x+298.6y+0.7508x.sup.2−10.1xy−14.69y.sup.2+2.184×10.sup.−3x.sup.3−0.08908x.sup.2y+0.7366xy.sup.2   (5-25)

[0129] The goodness of fit R2 of adaptive functions of the system in spring is 0.999, 0.983, 0.971, 0.949, 0.991, and 0.998 respectively; and if all values of the R2 are approximate to 1, the adaptive function group can commendably reflect a functional relationship between optimized parameters and the optimization objective and between the optimized parameters and the constraint conditions.

[0130] The IGV opening and natural gas flow of the combined cycle gas turbine system in the city in different seasons are optimized by means of the optimization model (as shown in FIG. 9 to FIG. 12). The primary energy ratio, the exergy efficiency, the per-unit emission amount of the CO.sub.2, and the per-unit thermoeconomic cost of an optimized system are analyzed by comparing overall evaluation results of the optimized system with overall evaluation results of a non-optimized system.

[0131] The overall evaluation results of the optimized system under different load conditions are higher than those of the non-optimized system. When the load of the system is 80%, the system is optimized to the greatest extent and has the overall evaluation result increased by 0.1576.

[0132] The above embodiments are only preferred ones of the present disclosure, and are not intended to limit the present disclosure in any form. Although the present disclosure has been disclosed by the foregoing embodiments, these embodiments are not intended to limit the present disclosure. Any person skilled in the art may make some changes or modifications to implement equivalent embodiments with equivalent changes by using the technical contents disclosed above without departing from the scope of the technical solution of the present disclosure. Any simple modification, equivalent change and modification made to the foregoing embodiments according to the technical essence of the present disclosure without departing from the content of the technical solution of the present disclosure shall fall within the scope of the technical solution of the present disclosure.