CN109472406A - 一种超低排放多污染物协同脱除系统的全局优化方法 - Google Patents

一种超低排放多污染物协同脱除系统的全局优化方法 Download PDF

Info

Publication number
CN109472406A
CN109472406A CN201811300698.3A CN201811300698A CN109472406A CN 109472406 A CN109472406 A CN 109472406A CN 201811300698 A CN201811300698 A CN 201811300698A CN 109472406 A CN109472406 A CN 109472406A
Authority
CN
China
Prior art keywords
cost
pollutant
concentration
flue gas
inlet
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201811300698.3A
Other languages
English (en)
Other versions
CN109472406B (zh
Inventor
高翔
郑成航
黄悦琪
郭杉
郭一杉
张涌新
翁卫国
吴卫红
曲瑞阳
刘少俊
赵海涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201811300698.3A priority Critical patent/CN109472406B/zh
Publication of CN109472406A publication Critical patent/CN109472406A/zh
Priority to PCT/CN2019/114203 priority patent/WO2020088485A1/zh
Priority to US17/263,869 priority patent/US20210319373A1/en
Application granted granted Critical
Publication of CN109472406B publication Critical patent/CN109472406B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0633—Workflow analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
    • B01D53/34—Chemical or biological purification of waste gases
    • B01D53/46—Removing components of defined structure
    • B01D53/48—Sulfur compounds
    • B01D53/50—Sulfur oxides
    • B01D53/501—Sulfur oxides by treating the gases with a solution or a suspension of an alkali or earth-alkali or ammonium compound
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
    • B01D53/34—Chemical or biological purification of waste gases
    • B01D53/46—Removing components of defined structure
    • B01D53/60—Simultaneously removing sulfur oxides and nitrogen oxides
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
    • B01D53/34—Chemical or biological purification of waste gases
    • B01D53/46—Removing components of defined structure
    • B01D53/64—Heavy metals or compounds thereof, e.g. mercury
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
    • B01D53/34—Chemical or biological purification of waste gases
    • B01D53/74—General processes for purification of waste gases; Apparatus or devices specially adapted therefor
    • B01D53/75—Multi-step processes
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
    • B01D53/34—Chemical or biological purification of waste gases
    • B01D53/74—General processes for purification of waste gases; Apparatus or devices specially adapted therefor
    • B01D53/86—Catalytic processes
    • B01D53/8621—Removing nitrogen compounds
    • B01D53/8625—Nitrogen oxides
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
    • B01D53/34—Chemical or biological purification of waste gases
    • B01D53/74—General processes for purification of waste gases; Apparatus or devices specially adapted therefor
    • B01D53/86—Catalytic processes
    • B01D53/8665—Removing heavy metals or compounds thereof, e.g. mercury
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D2251/00—Reactants
    • B01D2251/40—Alkaline earth metal or magnesium compounds
    • B01D2251/404—Alkaline earth metal or magnesium compounds of calcium
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D2251/00—Reactants
    • B01D2251/60—Inorganic bases or salts
    • B01D2251/606—Carbonates
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D2257/00—Components to be removed
    • B01D2257/30—Sulfur compounds
    • B01D2257/302—Sulfur oxides
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D2257/00—Components to be removed
    • B01D2257/40—Nitrogen compounds
    • B01D2257/404—Nitrogen oxides other than dinitrogen oxide
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D2257/00—Components to be removed
    • B01D2257/60—Heavy metals or heavy metal compounds
    • B01D2257/602—Mercury or mercury compounds
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01D—SEPARATION
    • B01D2258/00—Sources of waste gases
    • B01D2258/02—Other waste gases
    • B01D2258/0283—Flue gases

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Health & Medical Sciences (AREA)
  • Oil, Petroleum & Natural Gas (AREA)
  • General Chemical & Material Sciences (AREA)
  • Biomedical Technology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Analytical Chemistry (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Educational Administration (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Treating Waste Gases (AREA)

Abstract

本发明涉及一种超低排放多污染物协同脱除系统的全局优化方法，通过超低排放系统多装置多污染物协同脱除过程的精确建模，获得多污染物在多装置生成、迁移、转化和脱除过程的精确描述；通过超低排放系统中全局运行的成本评价方法精确评估不同负荷、煤质、污染物浓度和运行参数下的多污染物减排成本；通过多污染物多目标多工况的全局运行优化方法实现不同排放目标下全局污染减排装置减排量的分钟级规划和优化；通过低排放多污染物排放可靠达标的先进控制方法保证污染物的可靠减排和卡边控制。本发明从模型、架构和算法层面对超低排放系统的建模、优化和控制方法进行了创新，提高了系统的可控和可调能力，指导超低排放系统的高效、可靠和经济运行。

Description

一种超低排放多污染物协同脱除系统的全局优化方法
技术领域
本发明属于能源环境工程污染减排装置控制技术领域，具体涉及一种超低排放多污染物协同脱除系统的全局优化方法。
背景技术
我国煤炭资源丰富，能够稳定满足发展对电力的需求，决定了我国发电能源以煤为主的格局。这种能源结构和传统的煤利用方式产生了大量的大气污染物，随着我国经济的不断发展，对环保的要求也在逐步提高，燃煤机组实施清洁化生产大势所趋，燃煤电厂达到超低排放是一个发展趋势。
为满足环保要求，我国燃煤电厂都安装了除尘装置，SCR和WFGD也成为了脱硝和脱硫的主流技术。然而这些技术大多是单独开发的，各自考虑各自的边界条件，形成相对独立的工艺技术和技术装备，各污染物单独脱除，若设计运行不当，会导致一系列问题，如系统复杂，总占地面积偏大，总成本偏高，总阻力偏大，运行和维护难度大等。因此，建立一种超低排放多污染物协同脱除系统的全局优化方法意义重大，能够充分考虑燃煤电厂现有污染物脱除设备性能和协同关系，最大限度地减小环保改造对锅炉热效率，厂用电率和机组效率的影响，实现环保方案可行，技术路线优化，整体投资减少，节能降耗显著的目标，实现超低排放系统的高效绿色运行。
发明内容
本发明针对现有超低排放系统污染物脱除装备相互独立、运行波动大、污染物脱除成本高的问题，提出了一种多变燃料、多变负荷的超低排放多污染物协同脱除系统的全局优化方法。本发明从模型、架构和算法层面对超低排放系统的建模、优化和控制方法进行了创新，突破多工况多污染物多装置多参数的复杂系统建模和成本评估，多输入多输出多装置的多目标优化及大时延、非线性和大范围变负荷装置的控制难题，从而提高系统的可控和可调能力，指导超低排放系统的高效、可靠和经济运行。
为实现上述功能，本发明的技术方案是：
一种超低排放多污染物协同脱除系统的全局优化方法，通过超低排放系统多装置多污染物协同脱除过程的精确建模，获得多污染物(包括氮氧化物(NOx)、硫氧化物(SO2和SO3)、颗粒物(PM)和重金属(Hg)等)在多装置生成、迁移、转化和脱除过程的精确描述；通过超低排放系统中全局运行的成本评价方法精确评估不同负荷、煤质、污染物浓度和运行参数下的多污染物减排成本；通过多污染物多目标多工况的全局运行优化方法实现不同排放目标下全局污染减排装置减排量的分钟级规划和优化；通过低排放多污染物排放可靠达标的先进控制方法保证污染物的可靠减排和卡边控制。
作为优选，所述超低排放系统多装置多污染物协同脱除过程的精确建模包括多污染物生成过程精确预测，脱硫系统多污染物协同脱除过程的精确建模，脱硝系统多污染物协同脱除与转化过程的精确建模，除尘系统多污染物脱除过程精确建模；
多污染物生成过程精确预测：通过采集锅炉侧和污染物连续在线监测系统(CEMS)长时间运行数据，结合煤质检测数据和NOx、SO2、SO3、PM、Hg多种污染物测试报告，建立多种工况、多变煤质下的锅炉运行数据库和多污染物出口浓度数据库；通过数据建模的方法，基于锅炉运行数据库和出口污染物数据库，建立描述多种煤质、负荷等工况下锅炉及煤质参数与锅炉出口多种污染物浓度对应关系的模型；
脱硫系统多污染物协同脱除过程的精确建模：基于多种污染物生成、迁移转化及在脱硫系统的脱除机理，结合煤质、锅炉运行参数、脱硫系统入口烟气参数，与脱硫系统运行参数历史数据结合，建立：入口SO2浓度、入口烟气温度、液气比、浆液密度和pH值与脱硫系统出口SO2浓度的对应关系；入口SO3浓度、入口SO2浓度、入口含尘量、入口烟气温度、液气比和烟气流速与脱硫系统出口SO3浓度的对应关系；入口Hg浓度、入口含尘量、负荷、烟气流速、入口烟气温度和出口Hg浓度的对应关系；进而获得脱硫系统内多污染物的协同脱除精确描述模型；
脱硝系统多污染物协同脱除与转化过程的精确建模：基于选择性催化还原脱硝系统的运行机理，结合入口烟气参数，CEMS在线检测结果、脱硝系统运行参数的历史数据，建立：入口NOx浓度、烟气参数、反应条件和还原剂供给量与脱硝系统出口NOx浓度的对应关系；入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝系统内部SO2到SO3转化率的对应关系；入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝系统内SO2到SO3转化率的对应关系；入口Hg浓度、多种烟气参数、反应条件和还原剂供给量与脱硝系统内颗粒态汞(HgP)、氧化态汞(Hg2+、Hg+)以及单质汞(Hg0)之间转化比例的对应关系；基于这些对应关系，获得多种工况、多种运行参数下脱硝系统多污染物协同脱除与转化过程精确描述模型；
除尘系统多污染物协同脱除过程精确建模：基于电除尘系统的颗粒物脱除机理，结合煤质、锅炉运行参数和电除尘系统入口烟气参数，与运行历史数据结合，建立：干式电除尘装置煤质、负荷、烟气温度、烟气流量、运行电压与出口颗粒物浓度的对应关系；干式电除尘装置煤质、负荷、烟气温度、烟气流量、脱硝系统理论SO3浓度、运行电压与出口SO3浓度的对应关系；干式电除尘装置煤质、负荷、烟气温度、烟气流量、脱硝系统理论Hg浓度与出口Hg浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、循环水量、运行电压与出口颗粒物浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口NOx浓度、循环水量、循环水pH、运行电压与出口NOx浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口SO2浓度、循环水量、循环水pH、运行电压与出口SO2浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口SO3浓度、循环水量、循环水pH、运行电压与出口SO2浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口Hg浓度、循环水量、循环水pH、运行电压与出口Hg浓度的对应关系；进而获得除尘系统多污染物协同脱出精确描述模型。
作为优选，所述多污染物生成过程精确预测的方法是根据实时锅炉运行数据，结合定期更新的煤质报告，实现锅炉侧出口NOx、SO2、SO3、PM、Hg多种污染物浓度的精确预测；
按照以下步骤建立锅炉侧模型：
S101：采集锅炉测参数、煤质参数、锅炉出口检测报告；
S102：采集长期运行历史数据，包括锅炉测给煤量、烟风、水，及CEMS检测的锅炉出口污染物浓度；
S103：根据检测报告和在线运行数据，划分典型锅炉负荷区间，获得典型锅炉负荷区间内锅炉运行参数、煤质与出口污染物浓度的对应数据集；
S104：根据所划分的数据集，对于有CEMS在线监测数据的污染物，利用数据建模方法，获得锅炉测运行数据及煤质参数与出口污染物浓度的对应关系；
S105：根据所划分的数据集，对于没有CEMS在线监测数据的污染物，根据煤质数据和测试报告，利用数据建模方法，获得锅炉侧运行数据及煤质参数与出口污染物浓度的对应关系；
S106：将步骤S104和S105所得到的模型结合，获得锅炉侧的多污染物浓度预测模型。
作为优选，所述脱硫系统多污染物协同脱除精确描述模型的建立具体包括下述步骤：
S201：采集脱硫系统各装置设计参数，测试报告；
S202：采集脱硫系统在线运行的历史数据，包括锅炉负荷、入口烟气参数、烟气流量、pH值、液气比、浆液密度、入口和出口浓度；
S203：根据数据分布，划分典型负荷区间和污染物浓度区间，获得典型锅炉负荷区间和污染物区间内锅炉负荷、煤质与出口污染物浓度的对应数据集；
S204：根据所划分的数据集，基于脱硫塔装置脱除SO2机理，在典型锅炉负荷区间和污染物区间建立SO2脱除机理模型，并利用历史数据修正模型，获得典型锅炉负荷区间和污染物区间入口SO2浓度、入口烟气参数、液气比、浆液密度和pH值与出口SO2浓度的对应关系；
S205：根据所划分的数据集，结合测试报告结果，利用数据建模方法，获得典型锅炉负荷区间和污染物区间入口SO3浓度、入口SO2浓度、入口含尘量、入口烟气温度、液气比和烟气流速与脱硫系统出口SO3浓度的对应关系；入口Hg浓度、入口含尘量、负荷、烟气流速、入口烟气温度和出口Hg浓度的对应关系；
S206：基于这些对应关系，获得多种工况、多种运行参数下脱硫系统多污染物协同脱除与转化过程精确描述模型。
作为优选，所述脱硝系统多污染物协同脱除与转化过程精确描述模型的建立具体包括下述步骤：
S301：采集脱硝装置设计参数，催化剂参数及其测试报告，还原剂测试报告，测试报告；
S302：采集影响脱硝装置在线运行的历史数据，包括锅炉负荷、运行温度、烟气流量、还原剂供给量、CEMS检测的各污染物入口和出口浓度；
S303：根据数据分布，划分典型负荷区间和污染物浓度区间，获得典型锅炉负荷区间和污染物区间内锅炉负荷、煤质与出口污染物浓度的对应数据集；
S304：根据所划分的数据集，基于脱硝系统脱除NOx机理，在典型锅炉负荷区间和污染物区间建立脱硝装置机理模型，并利用历史数据修正模型，获得典型锅炉负荷区间和污染物区间入口NOx浓度、烟气参数、反应条件和还原剂供给量与脱硝装置出口NOx浓度的对应关系；
S305：根据所划分的数据集，结合测试报告结果，利用数据建模方法，获得典型锅炉负荷区间和污染物区间入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝装置内部SO2到SO3转化率的对应关系；入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝装置内SO2到SO3转化率的对应关系；入口Hg浓度、多种烟气参数、反应条件和还原剂供给量与脱硝装置内颗粒态汞(HgP)、氧化态汞(Hg2+、Hg+)以及单质汞(Hg0)之间转化比例的对应关系；
S306：基于这些对应关系，获得多种工况、多种运行参数下脱硝系统多污染物协同脱除与转化过程精确描述模型；
所述除尘系统多污染物协同脱出精确描述模型的建立具体包括下述步骤：
S401：采集电除尘装置设计参数，测试报告；
S402：采集影响电除尘装置在线运行的历史数据，包括锅炉负荷、运行温度、烟气流量、二次电压、二次电流、CEMS检测的污染物入口和出口浓度；
S403：根据数据分布，划分典型负荷区间和污染物浓度区间，获得典型锅炉负荷区间和污染物区间内锅炉负荷、煤质与出口污染物浓度的对应数据集；
S404：根据所划分的数据集，基于电除尘装置脱除PM机理，在典型锅炉负荷区间和污染物区间建立电除尘装置PM脱除机理模型，并利用历史数据修正模型，获得典型锅炉负荷区间和污染物区间入口PM浓度、烟气参数、运行电压与电除尘装置出口PM浓度的对应关系；
S405：根据所划分的数据集，结合测试报告结果，利用数据建模方法(如采用神经网络算法等)，获得典型锅炉负荷区间和污染物区间入口SO3浓度、入口PM浓度、锅炉负荷、运行温度、烟气流量、运行电压与出口SO3浓度的对应关系；入口Hg浓度、入口PM浓度、锅炉负荷、运行温度、烟气流量、运行电压与出口Hg浓度的对应关系；
S406：基于这些对应关系，获得多种工况、多种运行参数下除尘系统多污染物协同脱除与转化过程精确描述模型。
所述脱硝系统多污染物脱除与转化过程精确建模的方法根据脱硝系统实时运行状态，实现NOx、SO2、SO3和Hg动态转化过程的实时描述和浓度的精确预测；所述脱硫系统多污染物协同脱除过程的精确建模中，入口烟气参数包括入口烟温、入口SO2浓度、入口SO3浓度、入口颗粒物浓度、入口Hg浓度；出口烟气参数包括出口烟温、出口SO2浓度、出口SO3浓度、出口颗粒物浓度、出口Hg浓度；对建立的机理模型进行修正时采取耦合叠加的方式。
所述超低排放系统中全局运行的成本评价方法，具体包括各污染物减排系统的运行成本和固定成本；
针对石灰石-石膏湿法脱硫系统，总体运行成本和固定成本表示为：
COSTWFGD_fix＝COSTWFGD_d+COSTWFGD_r+COSTWFGD_m
其中包括增压风机电耗COSTbf、氧化风机电耗COSTsa、浆液循环泵电耗COSTscp、浆液搅拌器电耗COSToab、石灰石浆液成本工艺水成本COSTWFGD_w、副产物石膏收益折旧成本COSTWFGD_d、维修成本COSTWFGD_r、人工成本COSTWFGD_m。
脱硫系统的运行成本包括能耗与物耗，能耗主要由脱硫系统电机设备产生，所述能耗包括增压风机电耗、氧化风机电耗、浆液循环泵电耗、浆液搅拌器电耗，成本计算方法分别如下：
其中，q为锅炉实时负荷，nbf,nsa,nscp,noab分别表示增压风机、氧化风机、浆液循环泵、浆液搅拌器的运行台数，Ui,Ii分别为第i台设备的电压与电流，为功率因数，一般取0.8，PE为电价；
αWFGD表示脱硫塔阻力占后半段总阻力的比例，用下式计算得到：
其中，pdt是脱硫塔压降，pWESP是湿式电除尘器的阻力压降，pgd2是烟道部分阻力压降；
脱硫系统的物耗包括石灰石消耗和工艺水消耗，石灰石-石膏湿法脱硫系统的脱硫吸收剂为石灰石浆液，根据物料平衡，其单位发电量成本消耗，即石灰石消耗的计算方法为：
其中，与分别表示脱硫系统入口烟气与出口烟气所含SO2浓度，与分别为CaCO3与SO2的分子量，λ为石灰石纯度，为石灰石价格，δ1为钙硫比，一般设定为1.02～1.05，V表示烟气流量，与锅炉负荷成正相关，用下式计算得到：
V＝m×q×Vtc
其中，Vtc为单位燃煤产生的烟气量，m为燃煤质量；
所述工艺水消耗成的计算方法为：
其中，为水分子量，Pw为水价格；
因为在脱除烟气中SO2的同时产生副产物石膏，石膏作为脱硫系统运行过程中的收益部分被纳入成本计算石灰石-石膏湿法脱硫系统，收益计算方法为：
其中，为石膏分子量，为石膏价格；
脱硫系统的固定成本包括折旧成本COSTWFGD_d、维修成本COSTWFGD_r、人工成本COSTWFGD_m，其计算方法分别为：
其中，Q为机组容量，H为机组年运行小时数，PWFGD_init为脱硫系统初始投资成本，ηWFGD为脱硫系统固定资产形成率，YWFGD为脱硫系统折旧年限，ηWFGD_r为脱硫系统维修费用占投资成本比重，nWFGD为工作人员总人数，为脱硫系统第i个工作人员年薪。
针对选择性催化还原脱硝系统，总体运行成本和固定成本表示为：
COSTSCR＝COSTSCR_idf+COSTsb+COSTadf+COSTNH3+COSTC
COSTSCR_fix＝COSTSCR_d+COSTSCR_r+COSTSCR_m
所述运行成本包括能耗与物耗，所述能耗包括引风机电耗、吹灰风机电耗与稀释风机电耗，其计算方法分别为：
其中，nidf,nsb,nadf分别为引风机、吹灰电机与稀释风机的运行数目；αSCR表示脱硝反应器阻力占前半段总阻力的比例，计算方法为：
吹灰系统按吹灰方式不同，运行成本的计算方式也不同，其中，Psteam是经验蒸汽能耗，CVs是经验参考催化剂用量，CV是催化剂实际用量。
脱硝系统的物耗包括液氨与催化剂成本，根据物料平衡，液氨成本的计算方式为：
其中，分别表示脱硝系统入口烟气与出口烟气所含NOx浓度，与MNO分别为NH3与NO的分子量，表示液氨价格。
催化剂损耗成本的计算方法为：
其中，Pc为催化剂价格，Q为机组容量，h为机组年运行小时数；
脱硝系统的固定成本包括折旧成本COSTSCR_d、维修成本COSTSCR_r、人工成本COSTSCR_m，其计算方法为：
其中，PSCR_init为脱硝系统初始投资成本，ηSCR为脱硝系统固定资产形成率，YSCR为脱硝系统折旧年限，ηSCR_r为脱硝系统维修费用占投资成本比重，nSCR为脱硝系统工作人员总人数，为第i个工作人员年薪。
针对电除尘系统，干式静电除尘系统和湿式电除尘系统的总体运行成本和固定成本，其分别表示为：
COSTESP＝COSTESP_idf+COSTESP_e
COSTESP_fix＝COSTESP_d+COSTESP_r+COSTESP_m
COSTWESP＝COSTWESP_idf+COSTWESP_e+COSTWESP_w+COSTNa+COSTwc
COSTWESP_fix＝COSTWESP_d+COSTWESP_r+COSTWESP_m
所述运行成本包括电耗，干式静电除尘的电耗包括引风机电耗与电源电耗，计算方法分别为：
式中，ne表示电场数量；αESP为静电除尘器阻力占前半段总阻力的比例，计算方法为：
湿式电除尘器的电耗包括引风机电耗、电源电耗、电耗成本以及物料成本，湿式电除尘器的引风机电耗计算方法为：
湿式电除尘器的电源电耗为：
相较于干式静电除尘器，湿式静电除尘器增加了电耗成本以及物料成本，增加的电耗成本以水循环系统电耗为主，其计算方法为：
湿式电除尘器的物料成本包括工艺水成本以及碱消耗成本，计算方法为：
干式静电除尘系统与湿式静电除尘系统的固定成本均包括折旧成本COSTESP_d与COSTWESP_d、维修成本COSTESP_r与COSTWESP_r、人工成本COSTESP_m与COSTWESP_m，其计算方法分别为：
其中，PESP_init与PWESP_init分别为干式静电除尘系统与湿式静电除尘系统的初始投资成本，ηESP与ηWESP分别为干式静电除尘系统与湿式静电除尘系统的固定资产形成率，YESP与YWESP分别为干式静电除尘系统与湿式静电除尘系统的折旧年限，ηESP_r与ηWESP_r分别为干式静电除尘系统与湿式静电除尘系统的维修费用占投资成本比重，nESP与nWESP分别为干式静电除尘系统与湿式静电除尘系统的工作人员总人数，与分别为干式静电除尘系统与湿式静电除尘系统第i个工作人员年薪。
多污染物多目标多工况的全局运行优化方法，是基于多装置多污染物协同脱除过程的精确建模及全局运行成本评价，针对由选择性催化还原脱硝系统(SCR)、干式静电除尘系统(ESP)、湿法烟气脱硫系统(WFGD)与湿式静电除尘系统(WESP)四个主要部分组成的典型超低排放环保岛系统，利用各装置对多种污染物具有协同脱除的作用，以及各污染组分间存在的相互耦合与竞争关系，通过多污染物子学科分解和全局及局部的群智能算法，实现不同排放目标下的全局污染减排装置减排量的分钟级规划和优化。
作为优选，本发明采用基于协同优化算法的分散决策方法作为燃煤电厂超低排放环保岛多模型复杂系统运行优化问题的求解策略；该方法具有高度的系统自治性及良好的适应性，多个模型各自求解，为保证决策结果的一致性，各模型间需反复协调差异，有一定的沟通成本，但求解过程从各模型自身的最优解出发，决策结果较优；
协同优化算法首先由系统级传递耦合变量的期望值给各学科级，学科级在满足自身约束的条件下各自优化，使优化结果尽可能的靠近系统级给予的期望值，然后将各自的优化结果传递给系统级；系统级接收到各学科级的优化结果后，协调耦合变量并产生新的耦合变量期望值，然后再次传递给各学科级；经过系统级和各学科级之间的反复迭代，最终得到耦合变量保持一致且满足优化目标的可行解。
所述超低排放多污染物排放可靠达标的先进控制方法，是利用实时数据建立污染脱除装置操纵量及各扰动参数对污染物脱除动态性能响应的模型，在全局优化所给出的装置减排设定值的条件下，通过模型预测控制方法对污染物脱除装置的控制量进行实时优化调控；同时通过模型的在线更新，使得动态性能响应的模型更加符合系统所自带的大延迟、非线性、变负荷的特点，在系统参数的变化过程中，也能够实现污染物排放的卡边控制，进一步的在系统变化的情况下减小污染物的脱除成本；
需要建立起在高、中、低负荷情况下的，升负荷、降负荷、平稳负荷的九种污染物脱除动态性能响应模型，在负荷变化不同的阶段，运用不同的动态响应模型来计算污染物脱除装置控制量的大小，更符合系统实际变化的规律；
模型的在线更新是每经过某一时间段，对不满足动态响应性能的模型进行参数的替换，同时时间段的长度会根据脱硫脱硝除尘三大系统性能的不同而不同，脱硝系统喷氨迅速，变化较快，每1小时更新一次模型，脱硫系统中pH调节缓慢，每1天更新一次模型，除尘系统属于物理脱除过程，反应最快，每10分钟更新一次模型。
本发明的有益效果为：
相较于单一污染物脱除装置的控制研究，超低排放环保岛系统工艺流程长；从控制角度来看，超低排放环保岛系统是一个具有多输入多输出、强耦合、强非线性、变参数、多工况、大范围变负荷等特征的复杂流程工业过程。本发明在现有超低排放系统基础上，解决了环保岛中关键装置对多种烟气污染物的协同脱除问题，提升包括氮氧化物、硫氧化物、颗粒物和汞等多种污染物的协同脱除，有效提升烟气污染减排装置的多种污染物高效脱除效果，降低运行成本，提升多种工况下的运行适应性，具有较高的技术与应用价值，对超低排放系统的高效稳定可靠低成本绿色运行具有重大意义。
附图说明
图1为本发明实施案例提供的燃煤电厂环保岛污染物协同脱除示意图；
图2为本发明实施案例提供的协同优化算法优化过程示意图；
图3为本发明实施案例提供的燃煤电厂环保岛协同优化算法框架示意图；
图4为本发明实施案例提供的不同工况条件下环保岛运行成本构成图一；
图5为本发明实施案例提供的不同工况条件下环保岛运行成本构成图二。
具体实施方式
下面结合附图和具体实施例，对本发明作进一步说明，但本发明所要保护的范围并不限于此。
参照图1～3，一种超低排放多污染物协同脱除系统的全局优化方法，通过超低排放系统多装置多污染物协同脱除过程的精确建模，获得多污染物(包括氮氧化物(NOx)、硫氧化物(SO2和SO3)、颗粒物(PM)和重金属(Hg)等)在多装置生成、迁移、转化和脱除过程的精确描述；通过超低排放系统中全局运行的成本评价方法精确评估不同负荷、煤质、污染物浓度和运行参数下的多污染物减排成本；通过多污染物多目标多工况的全局运行优化方法实现不同排放目标下全局污染减排装置减排量的分钟级规划和优化；通过低排放多污染物排放可靠达标的先进控制方法保证污染物的可靠减排和卡边控制。本发明通过深度挖掘超低排放系统多装置多污染物协同脱除机理，结合运行数据与全局优化算法，在保证污染物高效脱除的前提下，充分发挥超低排放系统节能降耗潜力，突破多工况多污染物多装置多参数的复杂系统建模和成本评估，多输入多输出多装置的多目标优化及大时延、非线性和大范围变负荷装置的控制难题，从而提高系统的可控和可调能力，实现超低排放系统的高效、可靠和经济运行。
所述超低排放系统多装置多污染物协同脱除过程的精确建模包括多污染物生成过程精确预测，脱硫系统多污染物协同脱除过程的精确建模，脱硝系统多污染物协同脱除与转化过程的精确建模，除尘系统多污染物脱除过程精确建模。
多污染物生成过程精确预测：通过采集锅炉侧和污染物连续在线监测系统(CEMS)长时间运行数据，结合煤质检测数据和NOx、SO2、SO3、PM、Hg多种污染物测试报告，建立多种工况、多变煤质下的锅炉运行数据库和多污染物出口浓度数据库；通过数据建模的方法，基于锅炉运行数据库和出口污染物数据库，建立描述多种煤质、负荷等工况下锅炉及煤质参数与锅炉出口多种污染物浓度对应关系的模型。
脱硫系统多污染物协同脱除过程的精确建模：基于多种污染物生成、迁移转化及在脱硫系统的脱除机理，结合煤质、锅炉运行参数、脱硫系统入口烟气参数，与脱硫系统运行参数历史数据结合，建立：入口SO2浓度、入口烟气温度、液气比、浆液密度和pH值与脱硫系统出口SO2浓度的对应关系；入口SO3浓度、入口SO2浓度、入口含尘量、入口烟气温度、液气比和烟气流速与脱硫系统出口SO3浓度的对应关系；入口Hg浓度、入口含尘量、负荷、烟气流速、入口烟气温度和出口Hg浓度的对应关系；进而获得脱硫系统内多污染物的协同脱除精确描述模型。
脱硝系统多污染物协同脱除与转化过程的精确建模：基于选择性催化还原脱硝系统的运行机理，结合入口烟气参数，CEMS在线检测结果、脱硝系统运行参数的历史数据，建立：入口NOx浓度、烟气参数、反应条件和还原剂供给量与脱硝系统出口NOx浓度的对应关系；入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝系统内部SO2到SO3转化率的对应关系；入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝系统内SO2到SO3转化率的对应关系；入口Hg浓度、多种烟气参数、反应条件和还原剂供给量与脱硝系统内颗粒态汞(HgP)、氧化态汞(Hg2+、Hg+)以及单质汞(Hg0)之间转化比例的对应关系；基于这些对应关系，获得多种工况、多种运行参数下脱硝系统多污染物协同脱除与转化过程精确描述模型。
除尘系统多污染物协同脱除过程精确建模：基于电除尘系统的颗粒物脱除机理，结合煤质、锅炉运行参数和电除尘系统入口烟气参数，与运行历史数据结合，建立：干式电除尘装置煤质、负荷、烟气温度、烟气流量、运行电压与出口颗粒物浓度的对应关系；干式电除尘装置煤质、负荷、烟气温度、烟气流量、脱硝系统理论SO3浓度、运行电压与出口SO3浓度的对应关系；干式电除尘装置煤质、负荷、烟气温度、烟气流量、脱硝系统理论Hg浓度与出口Hg浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、循环水量、运行电压与出口颗粒物浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口NOx浓度、循环水量、循环水pH、运行电压与出口NOx浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口SO2浓度、循环水量、循环水pH、运行电压与出口SO2浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口SO3浓度、循环水量、循环水pH、运行电压与出口SO2浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口Hg浓度、循环水量、循环水pH、运行电压与出口Hg浓度的对应关系；进而获得除尘系统多污染物协同脱出精确描述模型。
作为优选，所述多污染物生成过程精确预测的方法是根据实时锅炉运行数据，结合定期更新的煤质报告，实现锅炉侧出口NOx、SO2、SO3、PM、Hg多种污染物浓度的精确预测；
按照以下步骤建立锅炉侧模型：
S101：采集锅炉测参数、煤质参数、锅炉出口检测报告；
S102：采集长期运行历史数据，包括锅炉测给煤量、烟风、水，及CEMS检测的锅炉出口污染物浓度；
S103：根据检测报告和在线运行数据，划分典型锅炉负荷区间(100％，75％，50％负荷等)，获得典型锅炉负荷区间内锅炉运行参数、煤质与出口污染物浓度的对应数据集；
S104：根据所划分的数据集，对于有CEMS在线监测数据的污染物，利用数据建模方法(本实施例采用神经网络算法)，获得锅炉测运行数据及煤质参数与出口污染物浓度的对应关系；
S105：根据所划分的数据集，对于没有CEMS在线监测数据的污染物，根据煤质数据和测试报告，利用数据建模方法(本实施例采用神经网络算法)，获得锅炉侧运行数据及煤质参数与出口污染物浓度的对应关系；
S106：将步骤S104和S105所得到的模型结合，获得锅炉侧的多污染物浓度预测模型。
所述脱硫系统多污染物协同脱除精确描述模型的建立具体包括下述步骤：
S201：采集脱硫系统各装置设计参数，测试报告；
S202：采集脱硫系统在线运行的历史数据，包括锅炉负荷、入口烟气参数、烟气流量、pH值、液气比、浆液密度、入口和出口浓度；
S203：根据数据分布，划分典型负荷区间和污染物浓度区间，获得典型锅炉负荷区间和污染物区间内锅炉负荷、煤质与出口污染物浓度的对应数据集；
S204：根据所划分的数据集，基于脱硫塔装置脱除SO2机理，在典型锅炉负荷区间和污染物区间建立SO2脱除机理模型，并利用历史数据修正模型，获得典型锅炉负荷区间和污染物区间入口SO2浓度、入口烟气参数、液气比、浆液密度和pH值与出口SO2浓度的对应关系；
S205：根据所划分的数据集，结合测试报告结果，利用数据建模方法，获得典型锅炉负荷区间和污染物区间入口SO3浓度、入口SO2浓度、入口含尘量、入口烟气温度、液气比和烟气流速与脱硫系统出口SO3浓度的对应关系；入口Hg浓度、入口含尘量、负荷、烟气流速、入口烟气温度和出口Hg浓度的对应关系；
S206：基于这些对应关系，获得多种工况、多种运行参数下脱硫系统多污染物协同脱除与转化过程精确描述模型。
所述脱硝系统多污染物协同脱除与转化过程精确描述模型的建立具体包括下述步骤：
S301：采集脱硝装置设计参数，催化剂参数及其测试报告，还原剂测试报告，测试报告；
S302：采集影响脱硝装置在线运行的历史数据，包括锅炉负荷、运行温度、烟气流量、还原剂供给量、CEMS检测的各污染物入口和出口浓度；
S303：根据数据分布，划分典型负荷区间和污染物浓度区间，获得典型锅炉负荷区间和污染物区间内锅炉负荷、煤质与出口污染物浓度的对应数据集；
S304：根据所划分的数据集，基于脱硝系统脱除NOx机理，在典型锅炉负荷区间和污染物区间建立脱硝装置机理模型，并利用历史数据修正模型，获得典型锅炉负荷区间和污染物区间入口NOx浓度、烟气参数、反应条件和还原剂供给量与脱硝装置出口NOx浓度的对应关系；
S305：根据所划分的数据集，结合测试报告结果，利用数据建模方法(本实施例采用神经网络算法)，获得典型锅炉负荷区间和污染物区间入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝装置内部SO2到SO3转化率的对应关系；入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝装置内SO2到SO3转化率的对应关系；入口Hg浓度、多种烟气参数、反应条件和还原剂供给量与脱硝装置内颗粒态汞、氧化态汞以及单质汞之间转化比例的对应关系；
S306：基于这些对应关系，获得多种工况、多种运行参数下脱硝系统多污染物协同脱除与转化过程精确描述模型；
所述除尘系统多污染物协同脱出精确描述模型的建立具体包括下述步骤：
S401：采集电除尘装置设计参数，测试报告；
S402：采集影响电除尘装置在线运行的历史数据，包括锅炉负荷、运行温度、烟气流量、二次电压、二次电流、CEMS检测的污染物入口和出口浓度；
S403：根据数据分布，划分典型负荷区间和污染物浓度区间，获得典型锅炉负荷区间和污染物区间内锅炉负荷、煤质与出口污染物浓度的对应数据集；
S404：根据所划分的数据集，基于电除尘装置脱除PM机理，在典型锅炉负荷区间和污染物区间建立电除尘装置PM脱除机理模型，并利用历史数据修正模型，获得典型锅炉负荷区间和污染物区间入口PM浓度、烟气参数、运行电压与电除尘装置出口PM浓度的对应关系；
S405：根据所划分的数据集，结合测试报告结果，利用数据建模方法，获得典型锅炉负荷区间和污染物区间入口SO3浓度、入口PM浓度、锅炉负荷、运行温度、烟气流量、运行电压与出口SO3浓度的对应关系；入口Hg浓度、入口PM浓度、锅炉负荷、运行温度、烟气流量、运行电压与出口Hg浓度的对应关系；
S406：基于这些对应关系，获得多种工况、多种运行参数下除尘系统多污染物协同脱除与转化过程精确描述模型。
所述SO2脱除机理建模包括下述部分：
SO2脱除速率计算：
其中，为二氧化硫摩尔浓度的微分，dt为时间的微分，a为单位体积内气液传质面积。
在喷淋塔中，从烟气中除去的SO2量等于浆液吸收的SO2量，对喷淋塔取微元高度dz，得到传质方程：
其中，KOG为气液两相的总传质系数，kmol/(m2·s)；a为单位体积内传质面积，m2/m3；P为总压力，Pa；y为SO2在气相的摩尔分数，y*为SO2在气膜侧的平衡摩尔分数，G为烟气流量，kmol/(m2·s)；dy气相的摩尔分数的微分；；
对脱硫系统进行简化：假设沿垂直方向流动恒定、总传质系数恒定不变，对上式dz进行积分得到吸收高度Z：
其中，(1-y)LM为(1-y*)和(1-y)的对数平均值，对于脱硫塔内SO2吸收过程可近似为1；Δym为脱硫塔气膜侧的平均推动力；y1为入口气相的摩尔分数；y2为出口气相的摩尔分数。
机理模型的基本形式为：
cout＝cin*exp(-NTU)
其中，cout、cin分别为SO2进口与出口浓度，NTU为传质单元数，其物理意义是表示吸收塔吸收SO2的困难程度，数学表达式为：
影响NTU的参数包括：浆液pH值、Ca/S比、液气比；
喷淋区总传质洗系数为：
根据非稳态渗透理论，kG为气膜吸收系数，kmol/(m2·s·kPa)；为化学反应增强因子；kL为液膜吸收系数，m/s；Sh为无量纲的施伍德数；为SO2的气相扩散系数，m2/s；p为塔内的总压，Pa；Sc为施密特数；Mair为空气的摩尔质量，g/mol；为SO2的摩尔质量，g/mol；Vair为空气的摩尔体积，cm3/mol；为SO2的摩尔体积，cm3/mol；DL为SO2在浆液中的扩散系数，m2/s；tp为渗透时间，s。
所述颗粒物脱除机理建模具体包括以下部分：
分级数目除尘效率ηn(dp,t)定义为t时刻、尺度为dp的烟尘颗粒被收集的数目总和与初始数目之间的比值：
分别代表尺度dp的颗粒区间的下限和上限，Np(dp,t)表示t时刻、尺度为dp的烟尘颗粒的总数目浓度；ddp代表尺度dp的微分，np(dp,t)表示t时刻、尺度为dp对应分级的烟尘颗粒数。
针对变量dp进行积分可得：
Np(dp,t)＝Np(dp,0)exp(-Λ(dp)×t)
Λ(dp)为沉积核(清除系数)，表示尺度为dp的颗粒被浆液液滴所脱除的速率，单位为s-1，表示为：
沉积核与分级数目质量除尘效率ηm(dp,t)的关系为：
ηm(dp,t)≈ηn(dp,t)＝1-exp(-Λ(dp)×t)
整体质量除尘效率ηtotal(t)定义为时刻t所有尺度颗粒被收集的质量总和与初始质量之间的比值：
式中ρp表示密度；np(dp,0)表示初始时刻、尺度为dp对应分级的烟尘颗粒数。
作为优选，所述脱硝系统多污染物脱除与转化过程精确建模的方法根据脱硝系统实时运行状态，实现NOx、SO2、SO3和Hg动态转化过程的实时描述和浓度的精确预测。
作为优选，所述脱硫系统多污染物协同脱除过程的精确建模中，入口烟气参数包括入口烟温、入口SO2浓度、入口SO3浓度、入口颗粒物(PM)浓度、入口Hg浓度；出口烟气参数包括出口烟温、出口SO2浓度、出口SO3浓度、出口颗粒物(PM)浓度、出口Hg浓度；对建立的机理模型进行修正时采取耦合叠加的方式。
采取耦合叠加的方式，形式如下：
机理模型主要考虑主要因素的影响，数据修正考虑次要因素的影响。
作为优选，所述超低排放系统中全局运行的成本评价方法，具体包括各污染物减排系统的运行成本和固定成本；
针对石灰石-石膏湿法脱硫系统，总体运行成本和固定成本表示为：
COSTWFGD_fix＝COSTWFGD_d+COSTWFGD_r+COSTWFGD_m
其中包括增压风机电耗COSTbf、氧化风机电耗COSTsa、浆液循环泵电耗COSTscp、浆液搅拌器电耗COSToab、石灰石浆液成本工艺水成本COSTWFGD_w、副产物石膏收益折旧成本COSTWFGD_d、维修成本COSTWFGD_r、人工成本COSTWFGD_m。
脱硫系统的运行成本包括能耗与物耗，能耗主要由脱硫系统电机设备产生，所述能耗包括增压风机电耗、氧化风机电耗、浆液循环泵电耗、浆液搅拌器电耗，成本计算方法分别如下：
其中，q为锅炉实时负荷，nbf,nsa,nscp,noab分别表示增压风机、氧化风机、浆液循环泵、浆液搅拌器的运行台数，Ui,Ii分别为第i台设备的电压与电流，为功率因数，一般取0.8，PE为电价；
αWFGD表示脱硫塔阻力占后半段总阻力的比例，用下式计算得到：
其中，pdt是脱硫塔压降，pWESP是湿式电除尘器的阻力压降，pgd2是烟道部分阻力压降；
脱硫系统的物耗包括石灰石消耗和工艺水消耗，石灰石-石膏湿法脱硫系统的脱硫吸收剂为石灰石浆液，根据物料平衡，其单位发电量成本消耗，即石灰石消耗的计算方法为：
其中，与分别表示脱硫系统入口烟气与出口烟气所含SO2浓度，与分别为CaCO3与SO2的分子量，λ为石灰石纯度，为石灰石价格，δ1为钙硫比，一般设定为1.02～1.05，V表示烟气流量，与锅炉负荷成正相关，用下式计算得到：
V＝m×q×Vtc
其中，Vtc为单位燃煤产生的烟气量，m为燃煤质量。
所述工艺水消耗成的计算方法为：
其中，为水分子量，Pw为水价格；
因为在脱除烟气中SO2的同时产生副产物石膏，石膏作为脱硫系统运行过程中的收益部分被纳入成本计算石灰石-石膏湿法脱硫系统，收益计算方法为：
其中，为石膏分子量，为石膏价格；
脱硫系统的固定成本包括折旧成本COSTWFGD_d、维修成本COSTWFGD_r、人工成本COSTWFGD_m，其计算方法分别为：
其中，Q为机组容量，H为机组年运行小时数，PWFGD_init为脱硫系统初始投资成本，ηWFGD为脱硫系统固定资产形成率，YWFGD为脱硫系统折旧年限，ηWFGD_r为脱硫系统维修费用占投资成本比重，nWFGD为工作人员总人数，SWFGDi为脱硫系统第i个工作人员年薪。
针对选择性催化还原脱硝系统，总体运行成本和固定成本表示为：
COSTSCR_fix＝COSTSCR_d+COSTSCR_r+COSTSCR_m
所述运行成本包括能耗与物耗，所述能耗包括引风机电耗、吹灰风机电耗与稀释风机电耗，其计算方法分别为：
其中，nidf,nsb,nadf分别为引风机、吹灰电机与稀释风机的运行数目；Psteam是经验蒸汽能耗，CVs是经验参考催化剂用量，CV是催化剂实际用量
αSCR表示脱硝反应器阻力占前半段总阻力的比例，计算方法为：
脱硝系统的物耗包括液氨与催化剂成本，根据物料平衡，液氨成本的计算方式为：
其中，其中，分别表示脱硝系统入口烟气与出口烟气所含NOx浓度，与MNO分别为NH3与NO的分子量，表示液氨价格。
催化剂损耗成本的计算方法为：
其中，Pc为催化剂价格，Q为机组容量，h为机组年运行小时数；
脱硝系统的固定成本包括折旧成本COSTSCR_d、维修成本COSTSCR_r、人工成本COSTSCR_m，其计算方法为：
其中，PSCR_init为脱硝系统初始投资成本，ηSCR为脱硝系统固定资产形成率，YSCR为脱硝系统折旧年限，ηSCR_r为脱硝系统维修费用占投资成本比重，nSCR为脱硝系统工作人员总人数，SSCRi为第i个工作人员年薪；
针对电除尘系统，干式静电除尘系统和湿式电除尘系统的总体运行成本和固定成本，其分别表示为：
COSTESP＝COSTESP_idf+COSTESP_e
COSTESP_fix＝COSTESP_d+COSTESP_r+COSTESP_m
COSTWESP＝COSTWESP_idf+COSTWESP_e+COSTWESP_w+COSTNa+COSTwc
COSTWESP_fix＝COSTWESP_d+COSTWESP_r+COSTWESP_m
所述运行成本包括电耗，
干式静电除尘的电耗包括引风机电耗与电源电耗，计算方法分别为：
式中，ne表示电场数量；αESP为静电除尘器阻力占前半段总阻力的比例，计算方法为：
湿式电除尘器的电耗包括引风机电耗、电源电耗、电耗成本以及物料成本，湿式电除尘器的引风机电耗计算方法为：
湿式电除尘器的电源电耗为：
相较于干式静电除尘器，湿式静电除尘器增加了电耗成本以及物料成本，增加的电耗成本以水循环系统电耗为主，其计算方法为：
湿式电除尘器的物料成本包括工艺水成本以及碱消耗成本，计算方法为：
干式静电除尘系统与湿式静电除尘系统的固定成本均包括折旧成本COSTESP_d与COSTWESP_d、维修成本COSTESP_r与COSTWESP_r、人工成本COSTESP_m与COSTWESP_m，其计算方法分别为：
其中，PESP_init与PWESP_init分别为干式静电除尘系统与湿式静电除尘系统的初始投资成本，ηESP与ηWESP分别为干式静电除尘系统与湿式静电除尘系统的固定资产形成率，YESP与YWESP分别为干式静电除尘系统与湿式静电除尘系统的折旧年限，ηESP_r与ηWESP_r分别为干式静电除尘系统与湿式静电除尘系统的维修费用占投资成本比重，nESP与nWESP分别为干式静电除尘系统与湿式静电除尘系统的工作人员总人数，与分别为干式静电除尘系统与湿式静电除尘系统第i个工作人员年薪。
图4、5展示了不同工况条件下环保岛系统的运行成本构成。在不同负荷与不同污染物浓度条件下，各个系统的运行成本占环保岛系统的整体运行成本比例基本不变。其中，湿法脱硫系统运行成本占比最大，约占环保岛整体运行成本的50％，其次是脱硝系统，约占30％的运行成本。干式静电除尘系统与湿式静电除尘系统的运行成本相当，均占据环保岛整体运行成本的10％左右。
本发明在高负荷高污染物浓度条件下，各污染物在环保岛系统内的脱除，NOx、PM、SO2的其实浓度分别为：250mg/m3、15356.6mg/m3、1200mg/m3。NOx的浓度在经过SCR系统后降低至40mg/m3左右，在WFGD系统的协同脱除作用下被脱除至35mg/m3以下。绝大部分PM在经过ESP系统时被脱除，ESP的PM脱除效率达到99％以上，其出口处的PM浓度仅为12.07mg/m3，最终在WFGD与WESP的脱除作用下，烟气中PM浓度被控制在1.25mg/m3。SO2主要在WFGD中被脱除，当烟气通过WFGD系统后，SO2浓度已低于20mg/m3，后续在WESP的协同脱除作用下，SO2浓度被控制在12.69mg/m3。
多污染物多目标多工况的全局运行优化方法，是基于多装置多污染物协同脱除过程的精确建模及全局运行成本评价，针对由选择性催化还原脱硝系统(SCR)、干式静电除尘系统(ESP)、湿法烟气脱硫系统(WFGD)与湿式静电除尘系统(WESP)四个主要部分组成的典型超低排放环保岛系统，利用各装置对多种污染物具有协同脱除的作用，以及各污染组分间存在的相互耦合与竞争关系，通过多污染物子学科分解和全局及局部的群智能算法，实现不同排放目标下的全局污染减排装置减排量的分钟级规划和优化。
NOx除了在SCR脱硝装置内发生选择性催化还原反应以外，还可被脱硫系统所吸收，主要反应为：
NO+NO2+H2O→HNO2
SO2在石灰石-石膏湿法脱硫装置内未被脱除的残留部分，能够在烟气经过湿式电除尘装置时随冲刷液的流动被脱除。PM除了干式电除尘装置与湿式电除尘装置这两种主要脱除装置以外，其余部分还会在脱硫塔内通过浆液洗涤作用被协同脱除。
作为优选，本发明采用基于协同优化算法的分散决策方法作为燃煤电厂超低排放环保岛多模型复杂系统运行优化问题的求解策略；该方法具有高度的系统自治性及良好的适应性，多个模型各自求解，为保证决策结果的一致性，各模型间需反复协调差异，有一定的沟通成本，但求解过程从各模型自身的最优解出发，决策结果较优。
协同优化算法首先由系统级传递耦合变量的期望值给各学科级，学科级在满足自身约束的条件下各自优化，使优化结果尽可能的靠近系统级给予的期望值，然后将各自的优化结果传递给系统级；系统级接收到各学科级的优化结果后，协调耦合变量并产生新的耦合变量期望值，然后再次传递给各学科级；经过系统级和各学科级之间的反复迭代，最终得到耦合变量保持一致且满足优化目标的可行解。
作为优选，本发明按不同污染物将问题分解为三个子学科分别进行优化，系统级优化的目标为最小化污染物减排的总成本，其约束为各子学科的耦合变量与系统级期望值之间的不一致性小于其对应的松弛因子；学科级的优化目标均为最小化自身耦合变量与系统期望值之间的不一致性，其约束均为污染物排放浓度低于相应污染物的排放限值；
系统级与学科级的优化问题均选用收敛速度快，适用范围广的粒子群算法求解；
学科级优化问题的设计变量有耦合变量y和非耦合变量x，其优化问题表示为：
s.t.gi(xi,yi)≤0
其中，yij表示第i个子学科的第j个耦合变量，表示系统级分配给子学科的第j个变量的期望值，si表示第i个子学科耦合变量的个数，n表示子学科的个数，xi表示第i个子学科的非耦合变量，gi表示第i个子学科的约束条件；
系统级优化问题的设计变量有耦合变量z和非耦合变量w，其优化问题表示为：
minF(z,w)
其中，F表示系统级目标函数，zj表示系统级第j个耦合变量，表示第i个学科的第j个耦合变量的优化结果，ε为松弛因子。
采用动态松弛因子能够较好的满足优化条件，以较快的收敛速度找到最优解。优化初期将松弛因子置一较大的初始值，构造较大的可行域，以确保可行解的存在。随着迭代的进行，缩小松弛因子，减小可行域，逐步提升耦合变量的一致性。
设在第k次迭代，第i个子学科的最优解为定义子学科间最大不一致性D为：设在第k-1次迭代，系统级的最优解为定义协同不一致性Pi为：则第k次迭代过程第i个子学科的松弛因子的构造为：其中，d取较小正数，以保证松弛作用，定义如下：
以各污染物脱除过程的RBF神经网络模型为基础，将协同优化算法的系统级与学科级之间的最大迭代次数设置为50次，将粒子群算法的最大迭代次数设置为50，粒子数目设置为100，对典型负荷及入口浓度条件下的工况进行优化。
所述超低排放多污染物排放可靠达标的先进控制方法，是利用实时数据建立污染脱除装置操纵量及各扰动参数对污染物脱除动态性能响应的模型，在全局优化所给出的装置减排设定值的条件下，通过模型预测控制方法对污染物脱除装置的控制量进行实时优化调控；同时通过模型的在线更新，使得动态性能响应的模型更加符合系统所自带的大延迟、非线性、变负荷的特点，在系统参数的变化过程中，也能够实现污染物排放的卡边控制，进一步的在系统变化的情况下减小污染物的脱除成本；
全局优化会向控制系统给出装置减排的设定值，控制系统会将此设定值作为一个输入变量，使得污染物脱除的控制效果尽可能的接近该设定值。利用实时数据建立污染脱除装置操纵量及各扰动参数对污染物脱除动态性能响应的模型，该模型具有在高、中、低负荷情况下，包含升负荷、降负荷、平稳负荷共九种子模型，在负荷变动的不同阶段，采用模型预测控制方法，并应用不同的子模型对污染物脱除装置的控制量进行实时优化调控。同时每间隔一段时间将会对动态性能响应的模型进行更新，使得在实际脱除模型发生变化的时候也能够很好的控制操作参数，考虑到脱硝、脱硫及除尘系统的特性，设置脱硝系统每1小时更新一次模型，脱硫系统每1天更新一次模型，除尘系统每10分钟更新一次模型。
通过以上实施方案，本发明可以实现在现有超低排放系统基础上，解决环保岛中关键装置对多种烟气污染物的协同脱除问题，提升包括氮氧化物、硫氧化物、颗粒物和汞等多种污染物的协同脱除，有效提升烟气污染减排装置的多种污染物高效脱除效果，大幅降低运行成本，提升多种工况下的运行适应性。

Claims (10)

1.一种超低排放多污染物协同脱除系统的全局优化方法，其特征在于：通过超低排放系统多装置多污染物协同脱除过程的精确建模，获得多污染物在多装置生成、迁移、转化和脱除过程的精确描述；通过超低排放系统中全局运行的成本评价方法精确评估不同负荷、煤质、污染物浓度和运行参数下的多污染物减排成本；通过多污染物多目标多工况的全局运行优化方法实现不同排放目标下全局污染减排装置减排量的分钟级规划和优化；通过低排放多污染物排放可靠达标的先进控制方法保证污染物的可靠减排和卡边控制。
2.根据权利要求1所述超低排放多污染物协同脱除系统的全局优化方法，其特征在于：所述超低排放系统多装置多污染物协同脱除过程的精确建模包括多污染物生成过程精确预测，脱硫系统多污染物协同脱除过程的精确建模，脱硝系统多污染物协同脱除与转化过程的精确建模，除尘系统多污染物脱除过程精确建模；
多污染物生成过程精确预测：通过采集锅炉侧和污染物连续在线监测系统长时间运行数据，结合煤质检测数据和NOx、SO2、SO3、PM、Hg多种污染物测试报告，建立多种工况、多变煤质下的锅炉运行数据库和多污染物出口浓度数据库；通过数据建模的方法，基于锅炉运行数据库和出口污染物数据库，建立描述多种煤质、负荷工况下锅炉及煤质参数与锅炉出口多种污染物浓度对应关系的模型；
脱硫系统多污染物协同脱除过程的精确建模：基于多种污染物生成、迁移转化及在脱硫系统的脱除机理，结合煤质、锅炉运行参数、脱硫系统入口烟气参数，与脱硫系统运行参数历史数据结合，建立：入口SO2浓度、入口烟气温度、液气比、浆液密度和pH值与脱硫系统出口SO2浓度的对应关系；入口SO3浓度、入口SO2浓度、入口含尘量、入口烟气温度、液气比和烟气流速与脱硫系统出口SO3浓度的对应关系；入口Hg浓度、入口含尘量、负荷、烟气流速、入口烟气温度与出口Hg浓度的对应关系；进而获得脱硫系统内多污染物的协同脱除精确描述模型；脱硝系统多污染物协同脱除与转化过程的精确建模：基于选择性催化还原脱硝系统的运行机理，结合入口烟气参数，CEMS在线检测结果、脱硝系统运行参数的历史数据，建立：入口NOx浓度、烟气参数、反应条件和还原剂供给量与脱硝系统出口NOx浓度的对应关系；入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝系统内部SO2到SO3转化率的对应关系；入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝系统内SO2到SO3转化率的对应关系；入口Hg浓度、多种烟气参数、反应条件和还原剂供给量与脱硝系统内颗粒态汞、氧化态汞以及单质汞之间转化比例的对应关系；基于这些对应关系，获得多种工况、多种运行参数下脱硝系统多污染物协同脱除与转化过程精确描述模型；
除尘系统多污染物协同脱除过程精确建模：基于电除尘系统的颗粒物脱除机理，结合煤质、锅炉运行参数和电除尘系统入口烟气参数，与运行历史数据结合，建立：干式电除尘装置煤质、负荷、烟气温度、烟气流量、运行电压与出口颗粒物浓度的对应关系；干式电除尘装置煤质、负荷、烟气温度、烟气流量、脱硝系统理论SO3浓度、运行电压与出口SO3浓度的对应关系；干式电除尘装置煤质、负荷、烟气温度、烟气流量、脱硝系统理论Hg浓度与出口Hg浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、循环水量、运行电压与出口颗粒物浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口NOx浓度、循环水量、循环水pH、运行电压与出口NOx浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口SO2浓度、循环水量、循环水pH、运行电压与出口SO2浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口SO3浓度、循环水量、循环水pH、运行电压与出口SO2浓度的对应关系；湿式静电除尘装置煤质、负荷、烟气温度、烟气流量、进口Hg浓度、循环水量、循环水pH、运行电压与出口Hg浓度的对应关系；进而获得除尘系统多污染物协同脱出精确描述模型。
3.根据权利要求2所述超低排放多污染物协同脱除系统的全局优化方法，其特征在于：所述多污染物生成过程精确预测的方法是根据实时锅炉运行数据，结合定期更新的煤质报告，实现锅炉侧出口NOx、SO2、SO3、PM、Hg多种污染物浓度的精确预测；
按照以下步骤建立锅炉侧模型：
S101：采集锅炉测参数、煤质参数、锅炉出口检测报告；
S102：采集长期运行历史数据，包括锅炉测给煤量、烟风、水，及CEMS检测的锅炉出口污染物浓度；
S103：根据检测报告和在线运行数据，划分典型锅炉负荷区间，获得典型锅炉负荷区间内锅炉运行参数、煤质与出口污染物浓度的对应数据集；
S104：根据所划分的数据集，对于有CEMS在线监测数据的污染物，利用数据建模方法，获得锅炉测运行数据及煤质参数与出口污染物浓度的对应关系；
S105：根据所划分的数据集，对于没有CEMS在线监测数据的污染物，根据煤质数据和测试报告，利用数据建模方法，获得锅炉侧运行数据及煤质参数与出口污染物浓度的对应关系；
S106：将步骤S104和S105所得到的模型结合，获得锅炉侧的多污染物浓度预测模型。
4.根据权利要求2超低排放多污染物协同脱除系统的全局优化方法，其特征在于所述脱硫系统多污染物协同脱除精确描述模型的建立具体包括下述步骤：
S201：采集脱硫系统各装置设计参数，测试报告；
S202：采集脱硫系统在线运行的历史数据，包括锅炉负荷、入口烟气参数、烟气流量、pH值、液气比、浆液密度、入口和出口浓度；
S203：根据数据分布，划分典型负荷区间和污染物浓度区间，获得典型锅炉负荷区间和污染物区间内锅炉负荷、煤质与出口污染物浓度的对应数据集；
S204：根据所划分的数据集，基于脱硫塔装置脱除SO2机理，在典型锅炉负荷区间和污染物区间建立SO2脱除机理模型，并利用历史数据修正模型，获得典型锅炉负荷区间和污染物区间入口SO2浓度、入口烟气参数、液气比、浆液密度和pH值与出口SO2浓度的对应关系；
S205：根据所划分的数据集，结合测试报告结果，利用数据建模方法，获得典型锅炉负荷区间和污染物区间入口SO3浓度、入口SO2浓度、入口含尘量、入口烟气温度、液气比和烟气流速与脱硫系统出口SO3浓度的对应关系；入口Hg浓度、入口含尘量、负荷、烟气流速、入口烟气温度和出口Hg浓度的对应关系；
S206：基于这些对应关系，获得多种工况、多种运行参数下脱硫系统多污染物协同脱除与转化过程精确描述模型；
所述脱硝系统多污染物协同脱除与转化过程精确描述模型的建立具体包括下述步骤：
S301：采集脱硝装置设计参数，催化剂参数及其测试报告，还原剂测试报告，测试报告；
S302：采集影响脱硝装置在线运行的历史数据，包括锅炉负荷、运行温度、烟气流量、还原剂供给量、CEMS检测的各污染物入口和出口浓度；
S303：根据数据分布，划分典型负荷区间和污染物浓度区间，获得典型锅炉负荷区间和污染物区间内锅炉负荷、煤质与出口污染物浓度的对应数据集；
S304：根据所划分的数据集，基于脱硝系统脱除NOx机理，在典型锅炉负荷区间和污染物区间建立脱硝装置机理模型，并利用历史数据修正模型，获得典型锅炉负荷区间和污染物区间入口NOx浓度、烟气参数、反应条件和还原剂供给量与脱硝装置出口NOx浓度的对应关系；
S305：根据所划分的数据集，结合测试报告结果，利用数据建模方法，获得典型锅炉负荷区间和污染物区间入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝装置内部SO2到SO3转化率的对应关系；入口SO2浓度、多种烟气参数、反应条件和还原剂供给量与脱硝装置内SO2到SO3转化率的对应关系；入口Hg浓度、多种烟气参数、反应条件和还原剂供给量与脱硝装置内颗粒态汞、氧化态汞以及单质汞之间转化比例的对应关系；
S306：基于这些对应关系，获得多种工况、多种运行参数下脱硝系统多污染物协同脱除与转化过程精确描述模型；
所述除尘系统多污染物协同脱出精确描述模型的建立具体包括下述步骤：
S401：采集电除尘装置设计参数，测试报告；
S402：采集影响电除尘装置在线运行的历史数据，包括锅炉负荷、运行温度、烟气流量、二次电压、二次电流、CEMS检测的污染物入口和出口浓度；
S403：根据数据分布，划分典型负荷区间和污染物浓度区间，获得典型锅炉负荷区间和污染物区间内锅炉负荷、煤质与出口污染物浓度的对应数据集；
S404：根据所划分的数据集，基于电除尘装置脱除颗粒物机理，在典型锅炉负荷区间和污染物区间建立颗粒物脱除机理模型，并利用历史数据修正模型，获得典型锅炉负荷区间和污染物区间入口PM浓度、烟气参数、运行电压与电除尘装置出口PM浓度的对应关系；
S405：根据所划分的数据集，结合测试报告结果，利用数据建模方法，获得典型锅炉负荷区间和污染物区间入口SO3浓度、入口PM浓度、锅炉负荷、运行温度、烟气流量、运行电压与出口SO3浓度的对应关系；入口Hg浓度、入口PM浓度、锅炉负荷、运行温度、烟气流量、运行电压与出口Hg浓度的对应关系；
S406：基于这些对应关系，获得多种工况、多种运行参数下除尘系统多污染物协同脱除与转化过程精确描述模型。
5.根据权利要求4所述超低排放多污染物协同脱除系统的全局优化方法，其特征在于：所述脱硝系统多污染物脱除与转化过程精确建模的方法根据脱硝系统实时运行状态，实现NOx、SO2、SO3和Hg动态转化过程的实时描述和浓度的精确预测；所述脱硫系统多污染物协同脱除过程的精确建模中，入口烟气参数包括入口烟温、入口SO2浓度、入口SO3浓度、入口颗粒物浓度、入口Hg浓度；出口烟气参数包括出口烟温、出口SO2浓度、出口SO3浓度、出口颗粒物浓度、出口Hg浓度；对建立的机理模型进行修正时采取耦合叠加的方式。
6.根据权利要求1所述超低排放多污染物协同脱除系统的全局优化方法，其特征在于：所述超低排放系统中全局运行的成本评价方法，具体包括各污染物减排系统的运行成本和固定成本；
针对石灰石-石膏湿法脱硫系统，总体运行成本和固定成本表示为：
COSTWFGD_fix＝COSTWFGD_d+COSTWFGD_r+COSTWFGD_m
其中包括增压风机电耗COSTbf、氧化风机电耗COSTsa、浆液循环泵电耗COSTscp、浆液搅拌器电耗COSToab、石灰石浆液成本工艺水成本COSTWFGD_w、副产物石膏收益折旧成本COSTWFGD_d、维修成本COSTWFGD_r、人工成本COSTWFGD_m；
脱硫系统的运行成本包括能耗与物耗，所述能耗包括增压风机电耗、氧化风机电耗、浆液循环泵电耗、浆液搅拌器电耗，成本计算方法分别如下：
其中，q为锅炉实时负荷，nbf,nsa,nscp,noab分别表示增压风机、氧化风机、浆液循环泵、浆液搅拌器的运行台数，Ui,Ii分别为第i台设备的电压与电流，为功率因数，PE为电价；
αWFGD表示脱硫塔阻力占后半段总阻力的比例，用下式计算得到：
其中，pdt是脱硫塔压降，pWESP是湿式电除尘器的阻力压降，pgd2是烟道部分阻力压降；
脱硫系统的物耗包括石灰石消耗和工艺水消耗，所述石灰石消耗的计算方法为：
其中，与分别表示脱硫系统入口烟气与出口烟气所含SO2浓度，与分别为CaCO3与SO2的分子量，λ为石灰石纯度，为石灰石价格，δ1为钙硫比，V表示烟气流量，与锅炉负荷成正相关，用下式计算得到：
V＝m×q×Vtc
其中，Vtc为单位燃煤产生的烟气量，m为燃煤质量；
所述工艺水消耗成的计算方法为：
其中，为水分子量，Pw为水价格；
石膏作为脱硫系统运行过程中的收益部分被纳入成本计算石灰石-石膏湿法脱硫系统，收益计算方法为：
其中，为石膏分子量，为石膏价格；
脱硫系统的固定成本包括折旧成本COSTWFGD_d、维修成本COSTWFGD_r、人工成本COSTWFGD_m，其计算方法分别为：
其中，Q为机组容量，H为机组年运行小时数，PWFGD_init为脱硫系统初始投资成本，ηWFGD为脱硫系统固定资产形成率，YWFGD为脱硫系统折旧年限，ηWFGD_r为脱硫系统维修费用占投资成本比重，nWFGD为工作人员总人数，为脱硫系统第i个工作人员年薪；
针对选择性催化还原脱硝系统，总体运行成本和固定成本表示为：
COSTSCR_fix＝COSTSCR_d+COSTSCR_r+COSTSCR_m
所述运行成本包括能耗与物耗，所述能耗包括引风机电耗、吹灰风机电耗与稀释风机电耗，其计算方法分别为：
其中，nidf,nsb,nadf分别为引风机、吹灰电机与稀释风机的运行数目；Psteam是经验蒸汽能耗，CVs是经验参考催化剂用量，CV是催化剂实际用量，
αSCR表示脱硝反应器阻力占前半段总阻力的比例，计算方法为：
脱硝系统的物耗包括液氨与催化剂成本，液氨成本的计算方式为：
其中，分别表示脱硝系统入口烟气与出口烟气所含NOx浓度，与MNO分别为NH3与NO的分子量，表示液氨价格；
催化剂损耗成本的计算方法为：
其中，Pc为催化剂价格，Q为机组容量，h为机组年运行小时数；
脱硝系统的固定成本包括折旧成本COSTSCR_d、维修成本COSTSCR_r、人工成本COSTSCR_m，其计算方法为：
其中，PSCR_init为脱硝系统初始投资成本，ηSCR为脱硝系统固定资产形成率，YSCR为脱硝系统折旧年限，ηSCR_r为脱硝系统维修费用占投资成本比重，nSCR为脱硝系统工作人员总人数，为第i个工作人员年薪；
针对电除尘系统，干式静电除尘系统和湿式电除尘系统的总体运行成本和固定成本，其分别表示为：
COSTESP＝COSTESP_idf+COSTESP_e
COSTESP_fix＝COSTESP_d+COSTESP_r+COSTESP_m
COSTWESP＝COSTWESP_idf+COSTWESP_e+COSTWESP_w+COSTNa+COSTwc
COSTWESP_fix＝COSTWESP_d+COSTWESP_r+COSTWESP_m
所述运行成本包括电耗，
干式静电除尘的电耗包括引风机电耗与电源电耗，计算方法分别为：
式中，ne表示电场数量；αESP为静电除尘器阻力占前半段总阻力的比例，计算方法为：
湿式电除尘器的电耗包括引风机电耗、电源电耗、电耗成本以及物料成本，湿式电除尘器的引风机电耗计算方法为：
湿式电除尘器的电源电耗为：
电耗成本计算方法为：
物料成本包括工艺水成本以及碱消耗成本，计算方法为：
干式静电除尘系统与湿式静电除尘系统的固定成本均包括折旧成本COSTESP_d与COSTWESP_d、维修成本COSTESP_r与COSTWESP_r、人工成本COSTESP_m与COSTWESP_m，其计算方法分别为：
其中，PESP_init与PWESP_init分别为干式静电除尘系统与湿式静电除尘系统的初始投资成本，ηESP与ηWESP分别为干式静电除尘系统与湿式静电除尘系统的固定资产形成率，YESP与YWESP分别为干式静电除尘系统与湿式静电除尘系统的折旧年限，ηESP_r与ηWESP_r分别为干式静电除尘系统与湿式静电除尘系统的维修费用占投资成本比重，nESP与nWESP分别为干式静电除尘系统与湿式静电除尘系统的工作人员总人数，与分别为干式静电除尘系统与湿式静电除尘系统第i个工作人员年薪。
7.根据权利要求1所述超低排放多污染物协同脱除系统的全局优化方法，其特征在于：所述多污染物多目标多工况的全局运行优化方法，是基于多装置多污染物协同脱除过程的精确建模及全局运行成本评价，针对由选择性催化还原脱硝系统、干式静电除尘系统、湿法烟气脱硫系统与湿式静电除尘系统四个部分组成的超低排放环保岛系统，利用各装置对多种污染物具有协同脱除的作用，以及各污染组分间存在的相互耦合与竞争关系，通过多污染物子学科分解和全局及局部的群智能算法，实现不同排放目标下的全局污染减排装置减排量的分钟级规划和优化。
8.根据权利要求7所述超低排放多污染物协同脱除系统的全局优化方法，其特征在于：所述超低排放环保岛系统的运行优化问题的求解策略采用基于协同优化算法的分散决策方法；
协同优化算法首先由系统级传递耦合变量的期望值给各学科级，学科级在满足自身约束的条件下各自优化，使优化结果尽可能的靠近系统级给予的期望值，然后将各自的优化结果传递给系统级；系统级接收到各学科级的优化结果后，协调耦合变量并产生新的耦合变量期望值，然后再次传递给各学科级；经过系统级和各学科级之间的反复迭代，最终得到耦合变量保持一致且满足优化目标的可行解。
9.根据权利要求1所述超低排放多污染物协同脱除系统的全局优化方法，其特征在于：所述超低排放多污染物排放可靠达标的先进控制方法，是利用实时数据建立污染脱除装置操纵量及各扰动参数对污染物脱除动态性能响应的模型，在全局优化所给出的装置减排设定值的条件下，通过模型预测控制方法对污染物脱除装置的控制量进行实时优化调控；同时通过模型的在线更新，使得动态性能响应的模型更加符合系统所自带的大延迟、非线性、变负荷的特点，在系统参数的变化过程中，也能够实现污染物排放的卡边控制，进一步的在系统变化的情况下减小污染物的脱除成本。
10.根据权利要求9所述超低排放多污染物协同脱除系统的全局优化方法，其特征在于：建立起在高、中、低负荷情况下的，升负荷、降负荷、平稳负荷的九种污染物脱除动态性能响应模型，在负荷变化不同的阶段，运用不同的动态响应模型来计算污染物脱除装置控制量的大小；
模型的在线更新是每经过某一时间段，对不满足动态响应性能的模型进行参数的替换，同时时间段的长度会根据脱硫脱硝除尘三大系统性能的不同而不同，脱硝系统每1小时更新一次模型，脱硫系统每1天更新一次模型，除尘系统每10分钟更新一次模型。
CN201811300698.3A 2018-11-02 2018-11-02 一种超低排放多污染物协同脱除系统的全局优化方法 Active CN109472406B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201811300698.3A CN109472406B (zh) 2018-11-02 2018-11-02 一种超低排放多污染物协同脱除系统的全局优化方法
PCT/CN2019/114203 WO2020088485A1 (zh) 2018-11-02 2019-10-30 一种智能化多种污染物超低排放系统及全局优化方法
US17/263,869 US20210319373A1 (en) 2018-11-02 2019-10-30 Intelligent multi-pollutant ultra-low emission system and global optimization method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811300698.3A CN109472406B (zh) 2018-11-02 2018-11-02 一种超低排放多污染物协同脱除系统的全局优化方法

Publications (2)

Publication Number Publication Date
CN109472406A true CN109472406A (zh) 2019-03-15
CN109472406B CN109472406B (zh) 2020-08-11

Family

ID=65672062

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811300698.3A Active CN109472406B (zh) 2018-11-02 2018-11-02 一种超低排放多污染物协同脱除系统的全局优化方法

Country Status (3)

Country Link
US (1) US20210319373A1 (zh)
CN (1) CN109472406B (zh)
WO (1) WO2020088485A1 (zh)

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110935312A (zh) * 2019-12-16 2020-03-31 广州珠江电力有限公司 一种scr烟气脱硝系统的动态监控装置及动态监测方法
CN110966611A (zh) * 2019-12-09 2020-04-07 苏州水云境科技有限公司 一种炉排炉垃圾焚烧排放全流程监管系统
CN110991787A (zh) * 2019-10-16 2020-04-10 华电电力科学研究院有限公司 一种火电厂超低排放后脱硫装置性能的评价方法
CN111018025A (zh) * 2019-12-31 2020-04-17 浙江大学 烟气余热利用的智能化废水零排放耦合脱氯系统
CN111027890A (zh) * 2019-12-26 2020-04-17 浙江天地环保科技有限公司 一种燃煤机组超低排放脱硫系统评估方法
WO2020088485A1 (zh) * 2018-11-02 2020-05-07 浙江大学 一种智能化多种污染物超低排放系统及全局优化方法
CN111178745A (zh) * 2019-12-26 2020-05-19 浙江天地环保科技有限公司 一种燃煤机组超低排放脱硝系统评估方法
CN111461388A (zh) * 2020-01-13 2020-07-28 杭州电子科技大学 一种基于燃煤锅炉超低排放系统的多目标优化方法
CN111467934A (zh) * 2020-04-15 2020-07-31 华电电力科学研究院有限公司 一种燃煤电厂湿式静电除尘器so3协同脱除效率预测方法
CN113011684A (zh) * 2021-04-26 2021-06-22 新奥数能科技有限公司 火电机组协同优化方法、装置、电子设备及可读介质
CN113031552A (zh) * 2021-03-09 2021-06-25 浙江菲达环保科技股份有限公司 炉后环保设备协同控制方法及系统
CN113058393A (zh) * 2021-04-16 2021-07-02 新奥数能科技有限公司 一种脱硫系统中浆液pH值的控制方法以及装置
CN113082954A (zh) * 2021-04-07 2021-07-09 浙江大学 一种湿法脱硫装置的全过程智能运行调控系统
CN113245060A (zh) * 2021-05-08 2021-08-13 西安热工研究院有限公司 除尘、脱硫系统协同节能运行优化方法、系统、设备及存储介质
CN113408905A (zh) * 2021-06-21 2021-09-17 西安热工研究院有限公司 湿法脱硫系统吸收塔能效评估方法、系统、设备及存储介质
CN113868876A (zh) * 2021-09-30 2021-12-31 湖南工商大学 一种烟气污染物控制设备的集成式设计方法
CN114159968A (zh) * 2021-12-10 2022-03-11 山西大学 一种电厂烟气重金属多污染物协同控制的预测方法
CN114307578A (zh) * 2021-12-31 2022-04-12 安徽省中易环保新材料有限公司 基于智能控制的工业烟气净化方法及系统
CN114397814A (zh) * 2021-12-06 2022-04-26 中国电建集团贵州电力设计研究院有限公司 一种基于bp神经网络的火电机组最佳运行参数寻找方法
CN114492961A (zh) * 2022-01-12 2022-05-13 河北大学 基于智能化KNN机制的燃煤锅炉NOx排放预测方法
CN114924889A (zh) * 2022-06-21 2022-08-19 浙江大学 一种基于云边端协同的超低排放智能调控系统及方法
CN115309117A (zh) * 2022-08-04 2022-11-08 浙江大学 一种基于数据驱动的wfgd出口so2浓度预测及智能优化方法
CN116339153A (zh) * 2023-05-22 2023-06-27 科大智能物联技术股份有限公司 基于粒子群优化的石灰炉窑控制方法
CN117435949A (zh) * 2023-08-29 2024-01-23 苏州方舟环保科技有限公司 基于人工智能的工业废水处理效果检测评价方法
CN118286838A (zh) * 2024-04-01 2024-07-05 浙江浩普智能科技有限公司 燃烧-污染物治理全流程智能调控减污降碳方法及系统

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11406933B1 (en) * 2019-06-03 2022-08-09 Henry Krigmont Flue gas conditioning system controller
CN111871164B (zh) * 2020-06-15 2023-03-21 浙江工业大学 一种低温等离子体催化降解CVOCs的智能调控方法
CN112099343B (zh) * 2020-07-29 2022-06-17 福建龙净环保股份有限公司 基于神经网络的电除尘系统智能节能优化方法及其介质
CN113082943B (zh) * 2021-03-22 2023-03-28 苏州新盛超净化系统有限公司 一种voc去除系统的节能方法
CN113239565B (zh) * 2021-05-27 2022-07-12 中南大学 一种沸腾焙烧炉产品质量软测量方法及装置
CN113457396B (zh) * 2021-06-29 2023-12-22 华能伊敏煤电有限责任公司 基于火力发电厂scr脱硝控制系统喷氨自动控制方法
CN113468813B (zh) * 2021-07-07 2022-04-15 大唐环境产业集团股份有限公司 一种脱硫系统入口so2浓度预测方法、装置、电子设备
CN113609684A (zh) * 2021-08-09 2021-11-05 工数科技(广州)有限公司 基于工业数据和工艺机理的锅炉吨煤产汽优化方法
CN113976322A (zh) * 2021-10-29 2022-01-28 西安热工研究院有限公司 适用于燃煤电厂电除尘系统的节能潜力评估方法及系统
CN113941251A (zh) * 2021-11-08 2022-01-18 西安热工研究院有限公司 一种烟气脱硝流场模拟校核方法和系统
CN114114921A (zh) * 2021-11-26 2022-03-01 华能平凉发电有限责任公司 一种除尘电源的控制方法及装置
CN114037181B (zh) * 2021-11-26 2024-08-06 江苏科技大学 一种基于pso优化bp神经网络的耙吸挖泥船艏吹浓度预测方法
CN114191953B (zh) * 2021-12-07 2024-02-20 国网河北能源技术服务有限公司 一种基于卷积神经网络和XGBoost的烟气脱硫脱硝控制方法
CN114191982A (zh) * 2021-12-10 2022-03-18 广州卓博机电科技有限公司 一种分散氨空混合器及氨氮一体分区轮测智能喷氨系统
CN114472491B (zh) * 2021-12-31 2022-12-20 中国科学院沈阳应用生态研究所 一种分段热脱附工艺处置油污土壤的物料衡算与工艺优化方法
CN114740713B (zh) * 2021-12-31 2024-08-09 北京工业大学 一种湿法烟气脱硫过程的多目标优化控制方法
CN114444240B (zh) * 2022-01-28 2022-09-09 暨南大学 一种面向信息物理融合系统的延迟和寿命优化方法
CN114791102B (zh) * 2022-04-21 2023-09-22 中国矿业大学 一种基于动态运行数据分析的燃烧优化控制方法
CN114935627A (zh) * 2022-05-07 2022-08-23 金川集团信息与自动化工程有限公司 一种用于预测合成炉出口氯化氢纯度的方法
CN115111601B (zh) * 2022-07-04 2024-05-14 浙江大学 多变负荷下内嵌算法融合的多目标锅炉燃烧优化控制方法
CN115495879B (zh) * 2022-08-18 2024-04-30 山东科技大学 一种工作面近正断层开采保护煤柱宽度设计方法
CN115173490B (zh) * 2022-09-05 2022-11-25 北京邮电大学 综合能源站供能方法和装置、电子设备
TWI839196B (zh) * 2023-04-27 2024-04-11 台灣電力股份有限公司 長期電力負載預測系統
CN116727107A (zh) * 2023-05-12 2023-09-12 江苏南通发电有限公司 一种除灰排放控制的方法
CN116362164B (zh) * 2023-06-02 2023-08-11 华北电力科学研究院有限责任公司 空气预热器堵塞状况确定方法及装置
CN117313428B (zh) * 2023-11-28 2024-02-23 西安中核核仪器股份有限公司 一种火电厂污染物浓度在线监测修正方法
CN117948601A (zh) * 2024-01-29 2024-04-30 光大环保（中国）有限公司 一种垃圾焚烧污染物排放总量的afc控制方法和系统
CN118113968B (zh) * 2024-04-30 2024-07-12 潍坊鑫山环保重工科技有限公司 一种污水处理流量在线监测系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104504498A (zh) * 2014-12-04 2015-04-08 国家电网公司 一种燃煤发电机组超低排放环保电价监控方法
CN109472406B (zh) * 2018-11-02 2020-08-11 浙江大学 一种超低排放多污染物协同脱除系统的全局优化方法

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020088485A1 (zh) * 2018-11-02 2020-05-07 浙江大学 一种智能化多种污染物超低排放系统及全局优化方法
CN110991787A (zh) * 2019-10-16 2020-04-10 华电电力科学研究院有限公司 一种火电厂超低排放后脱硫装置性能的评价方法
CN110966611A (zh) * 2019-12-09 2020-04-07 苏州水云境科技有限公司 一种炉排炉垃圾焚烧排放全流程监管系统
CN110935312A (zh) * 2019-12-16 2020-03-31 广州珠江电力有限公司 一种scr烟气脱硝系统的动态监控装置及动态监测方法
CN111027890A (zh) * 2019-12-26 2020-04-17 浙江天地环保科技有限公司 一种燃煤机组超低排放脱硫系统评估方法
CN111178745A (zh) * 2019-12-26 2020-05-19 浙江天地环保科技有限公司 一种燃煤机组超低排放脱硝系统评估方法
CN111018025B (zh) * 2019-12-31 2024-01-02 浙江大学 烟气余热利用的智能化废水零排放耦合脱氯系统
CN111018025A (zh) * 2019-12-31 2020-04-17 浙江大学 烟气余热利用的智能化废水零排放耦合脱氯系统
CN111461388A (zh) * 2020-01-13 2020-07-28 杭州电子科技大学 一种基于燃煤锅炉超低排放系统的多目标优化方法
CN111467934A (zh) * 2020-04-15 2020-07-31 华电电力科学研究院有限公司 一种燃煤电厂湿式静电除尘器so3协同脱除效率预测方法
CN113031552A (zh) * 2021-03-09 2021-06-25 浙江菲达环保科技股份有限公司 炉后环保设备协同控制方法及系统
CN113031552B (zh) * 2021-03-09 2022-10-25 浙江菲达环保科技股份有限公司 炉后环保设备协同控制方法及系统
CN113082954A (zh) * 2021-04-07 2021-07-09 浙江大学 一种湿法脱硫装置的全过程智能运行调控系统
CN113058393A (zh) * 2021-04-16 2021-07-02 新奥数能科技有限公司 一种脱硫系统中浆液pH值的控制方法以及装置
CN113011684A (zh) * 2021-04-26 2021-06-22 新奥数能科技有限公司 火电机组协同优化方法、装置、电子设备及可读介质
CN113245060A (zh) * 2021-05-08 2021-08-13 西安热工研究院有限公司 除尘、脱硫系统协同节能运行优化方法、系统、设备及存储介质
WO2022237011A1 (zh) * 2021-05-08 2022-11-17 西安热工研究院有限公司 除尘、脱硫系统协同节能运行优化方法、系统、设备及存储介质
CN113408905B (zh) * 2021-06-21 2023-03-14 西安热工研究院有限公司 湿法脱硫系统吸收塔能效评估方法、系统、设备及存储介质
CN113408905A (zh) * 2021-06-21 2021-09-17 西安热工研究院有限公司 湿法脱硫系统吸收塔能效评估方法、系统、设备及存储介质
CN113868876A (zh) * 2021-09-30 2021-12-31 湖南工商大学 一种烟气污染物控制设备的集成式设计方法
CN114397814A (zh) * 2021-12-06 2022-04-26 中国电建集团贵州电力设计研究院有限公司 一种基于bp神经网络的火电机组最佳运行参数寻找方法
CN114159968A (zh) * 2021-12-10 2022-03-11 山西大学 一种电厂烟气重金属多污染物协同控制的预测方法
CN114159968B (zh) * 2021-12-10 2023-09-26 山西大学 一种电厂烟气重金属多污染物协同控制的预测方法
CN114307578A (zh) * 2021-12-31 2022-04-12 安徽省中易环保新材料有限公司 基于智能控制的工业烟气净化方法及系统
CN114492961A (zh) * 2022-01-12 2022-05-13 河北大学 基于智能化KNN机制的燃煤锅炉NOx排放预测方法
CN114492961B (zh) * 2022-01-12 2024-05-14 河北大学 基于智能化KNN机制的燃煤锅炉NOx排放预测方法
CN114924889A (zh) * 2022-06-21 2022-08-19 浙江大学 一种基于云边端协同的超低排放智能调控系统及方法
CN114924889B (zh) * 2022-06-21 2024-05-14 浙江大学 一种基于云边端协同的超低排放智能调控系统及方法
CN115309117A (zh) * 2022-08-04 2022-11-08 浙江大学 一种基于数据驱动的wfgd出口so2浓度预测及智能优化方法
CN115309117B (zh) * 2022-08-04 2024-05-14 浙江大学 一种基于数据驱动的wfgd出口so2浓度预测及智能优化方法
CN116339153B (zh) * 2023-05-22 2023-09-01 科大智能物联技术股份有限公司 基于粒子群优化的石灰炉窑控制方法
CN116339153A (zh) * 2023-05-22 2023-06-27 科大智能物联技术股份有限公司 基于粒子群优化的石灰炉窑控制方法
CN117435949A (zh) * 2023-08-29 2024-01-23 苏州方舟环保科技有限公司 基于人工智能的工业废水处理效果检测评价方法
CN118286838A (zh) * 2024-04-01 2024-07-05 浙江浩普智能科技有限公司 燃烧-污染物治理全流程智能调控减污降碳方法及系统

Also Published As

Publication number Publication date
US20210319373A1 (en) 2021-10-14
CN109472406B (zh) 2020-08-11
WO2020088485A1 (zh) 2020-05-07

Similar Documents

Publication Publication Date Title
CN109472406B (zh) 一种超低排放多污染物协同脱除系统的全局优化方法
CN113082954B (zh) 一种湿法脱硫装置的全过程智能运行调控系统
CN107243257B (zh) 适合全负荷的智能喷氨控制系统
CN112580250A (zh) 基于深度学习的火电机组脱硝系统及优化控制方法
CN102728217B (zh) 移动床活性焦联合脱硫脱硝脱汞的方法及系统
CN105404145B (zh) 基于指数预测和时滞预估补偿的脱硝新型串级控制方法
CN107748955A (zh) 一种燃煤电厂超低排放环保岛能效评估方法
CN204543996U (zh) 燃煤电厂烟气多污染物深度净化控制耦合系统
CN102921287B (zh) 双塔双循环石灰石湿法脱硫装置及方法
CN104226110A (zh) 一种燃煤锅炉scr脱硝控制方法与控制系统
CN113050559B (zh) 燃煤电厂脱硫系统与电除尘系统协同控制方法及系统
CN110263452A (zh) 一种烟道内烟气时间分布特性分析方法、系统及脱硝系统
US20060121616A1 (en) Method and system for increasing efficiency of FGD operation in fossil fuel boilers
CN112364562A (zh) 一种烟气环保岛协同控制方法及系统
CN109876658B (zh) 分区涡流卷吸喷氨系统及喷氨调节方法
CN106557027B (zh) 优化焦化烟气脱硫脱硝过程中臭氧运行的系统及方法
CN105233657A (zh) 一种烟气脱硫除尘深度处理工艺及设备
CN203777877U (zh) 烟气复合除尘除雾装置
CN205216550U (zh) 一种基于动态反应区的高效半干法脱硫塔及脱硫系统
CN110865623A (zh) 一种SCR脱硝控制中NOx测量信号替代系统及其控制方法
CN209254469U (zh) 一种烟气净化装置
CN202683069U (zh) 移动床活性焦联合脱硫脱硝脱汞的系统
CN115509189A (zh) 一种应用于湿法脱硫系统的大数据调控方法及系统
CN108126517A (zh) 一种电厂用的水解供氨脱硝系统及其控制方法
CN109933884A (zh) 一种面向燃煤机组scr脱硝系统神经网络逆控制的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant