火力发电厂脱硫系统液固两相流仿真建模方法

1.一种火力发电厂脱硫系统液固两相流仿真建模方法，其特征在于，包括如下步骤：
S1，建立火力发电厂脱硫系统液固两相流的等效物理模型；
S2，选择对火力发电厂脱硫系统液固两相流进行有效控制的工艺操作参数，设为等效模型的输入变量；
S3，根据所述输入变量和预设脱硫效率，计算浆液槽内平衡氢离子浓度；
S4，根据所述浆液槽内平衡氢离子浓度，计算循环浆液中二氧化硫摩尔分数；
S5，根据所述循环浆液中二氧化硫摩尔分数，计算火力发电厂脱硫系统液固两相流的脱硫效率；
S6，判断所述脱硫效率是否和预设脱硫效率一致；若是，执行S7；若否，将所述脱硫效率置为预设脱硫效率，重复执行S3～S6的步骤，直到所述脱硫效率达到预设脱硫效率为止，执行S7；
S7，将所述等效物理模型设为火力发电厂脱硫系统液固两相流的仿真模型，结束仿真操作。
2.根据权利要求1所述的火力发电厂脱硫系统液固两相流仿真建模方法，其特征在于，所述根据所述输入变量和预设脱硫效率，计算浆液槽内平衡氢离子浓度的公式如下：
其中， 表示浆液中平衡氢离子浓度， 表示从烟气中脱除的二氧化硫量， 表示CaCO3的消耗量， 表示CaCO3的摩尔分数，VC浆液槽中浆液体积，V表示原烟气的总体积， 表示脱硫效率，yin表示进口原烟气进口二氧化硫的体积分数，Wlm表示含石灰石浆液量， 表示含石灰石浆液质量浓度，mC浆液槽中浆液质量，ρc浆液槽中浆液密度。
3.根据权利要求1所述的火力发电厂脱硫系统液固两相流仿真建模方法，其特征在于，所述根据所述浆液槽内平衡氢离子浓度，计算循环浆液中二氧化硫的摩尔分数的公式如下所示：
其中，xin表示循环浆液中二氧化硫的摩尔分数，θ表示浆液中剩余二氧化硫的摩尔分数与平衡氢离子摩尔分数的平衡系数，VC浆液槽中浆液体积， 表示浆液中平衡氢离子浓度，mC浆液槽中浆液质量，MC表示浆液的摩尔质量。
4.根据权利要求1所述的火力发电厂脱硫系统液固两相流仿真建模方法，其特征在于，所述根据所述循环浆液中二氧化硫摩尔分数，计算火力发电厂脱硫系统液固两相流的脱硫效率的步骤包括：
计算火力发电厂脱硫系统的传质单元高度；
根据所述传质单元高度，计算火力发电厂脱硫系统的传质单元数；
利用亨利定律，计算出口净烟气二氧化硫的体积分数；
根据进口原烟气进口二氧化硫的体积分数和出口净烟气二氧化硫的体积分数，计算火力发电厂脱硫系统液固两相流的脱硫效率。
5.根据权利要求4所述的火力发电厂脱硫系统液固两相流仿真建模方法，其特征在于，所述计算火力发电厂脱硫系统的传质单元高度的公式如下所示：
其中，HOG表示传质单元高度，V表示原烟气的总体积，KY表示总吸收系数，a表示单位体积塔层提供的有效传质体积系数，Ω表示吸收塔的截面积。
6.根据权利要求4所述的火力发电厂脱硫系统液固两相流仿真建模方法，其特征在于，所述根据所述传质单元高度，计算火力发电厂脱硫系统的传质单元数的公式如下所示：
其中，NOG表示传质单元数，Z表示喷淋层的高度，HOG表示传质单元高度，s表示脱硫系数，m表示相平衡系数，V表示原烟气的总体积，yin表示进口原烟气进口二氧化硫的体积分数，yout表示出口净烟气二氧化硫的体积分数， 表示溶质在气相的摩尔分数，L表示浆液的摩尔浓度，E表示亨利系数，p表示总压，β表示修正系数。
7.根据权利要求4所述的脱硫系统液固两相流仿真建模方法，其特征在于，所述利用亨利定律，计算出口净烟气二氧化硫的体积分数的公式如下所示：
其中， 表示溶质在气相的摩尔分数，m表示相平衡系数，xin表示循环浆液中的摩尔分数，L表示浆液的摩尔浓度，xout表示喷淋层后循环浆液中的摩尔分数，V表示原烟气的总体积， 表示脱硫效率，yin表示进口原烟气进口二氧化硫的体积分数。
8.根据权利要求4所述的火力发电厂脱硫系统液固两相流仿真建模方法，其特征在于，根据所述进口原烟气进口二氧化硫的体积分数和出口净烟气二氧化硫的体积分数，计算火力发电厂脱硫系统液固两相流的脱硫效率的公式如下所示：
式中， 表示脱硫效率，yin表示进口原烟气进口二氧化硫的体积分数，yout表示出口净烟气二氧化硫的体积分数。
9.根据权利要求1所述的火力发电厂脱硫系统液固两相流仿真建模方法，其特征在于，在所述根据所述输入变量和预设脱硫效率，计算浆液槽内平衡氢离子浓度的步骤之前，还包括：计算脱硫系统制浆区浆液箱内的石膏密度。
10.根据权利要求9所述的火力发电厂脱硫系统液固两相流仿真建模方法，其特征在于，所述计算脱硫系统制浆区浆液箱内的石膏密度的公式如下所示：
RO＝KA*KB/(KB-RATIO*(KB-KA))；
其中，RO表示浆液箱内的石膏密度，RATIO表示浆液箱内的石膏质量浓度，KA表示液相密度，KB表示固相密度。

火力发电厂脱硫系统液固两相流仿真建模方法\n技术领域\n[0001] 本发明涉及火力发电厂污染物控制技术领域，特别是涉及一种火力发电厂脱硫系统液固两相流仿真建模方法。\n背景技术\n[0002] 近些年来，随着电力事业的飞速发展，火力发电厂的火电装机总容量也随之增加，这大大增加了火电厂的耗煤量。耗煤量的大量增加，由煤燃烧产生的污染物，像飞灰、氮氧化物和二氧化硫对环境的污染也变得愈来愈严重，因此，严格控制燃煤产生的污染物的排放成为电力事业发展的重要组成部分。\n[0003] 火力发电厂锅炉尾部烟气中的二氧化硫排放量控制是治理大气污染物的重要一环，而石灰石-石膏湿法脱硫工艺是当前应用范围最广、工艺技术最成熟的标准脱硫工艺技术。石灰石-石膏湿法脱硫工艺是湿法脱硫的一种，是当前大机组火电厂烟气脱硫的基本工艺，它采用价廉易得的石灰石或石灰作脱硫吸收剂，石灰石经破碎磨细成粉状与水混合搅拌成吸收浆液，当采用石灰为吸收剂时，石灰粉经消化处理后加水制成吸收剂浆液；在吸收塔内，吸收浆液与烟气接触混合，烟气中的二氧化硫与浆液中的碳酸钙以及鼓入的氧化空气进行化学反应被脱除，最终反应产物为石膏；脱硫后的烟气经除雾器除去带出的细小液滴，经换热器加热升温后排入烟囱；脱硫石膏浆经过脱水装置脱水后进行回收。通过上述操作，可以大大降低煤燃烧排放的烟气对环境的污染程度；电站仿真机可以对火力发电厂脱硫系统内部的实际运行情况进行模拟，并获取相应的数据，便于工程技术人员掌握火力发电厂脱硫系统内部的运行情况。\n[0004] 现有的电站仿真机对火力发电厂脱硫系统采用的仿真建模方法，现有的电站仿真机采用的仿真建模方法只能获取火力发电厂脱硫系统的静态数据，来对火力发电厂脱硫系统的脱硫效率进行计算，无法实时反映火力发电厂脱硫系统内部的运行情况。\n发明内容\n[0005] 基于此，有必要针对现有的电站仿真机采用的仿真建模方法只能获取火力发电厂脱硫系统的静态数据，无法实时反映火力发电厂脱硫系统内部的运行情况的技术问题，提供一种火力发电厂脱硫系统液固两相流仿真建模方法。\n[0006] 一种火力发电厂脱硫系统液固两相流仿真建模方法，包括如下步骤：\n[0007] 步骤S1，建立火力发电厂脱硫系统液固两相流的等效物理模型；\n[0008] 在上述步骤的等效物理模型中，火力发电厂脱硫系统包括：吸收塔、浆液泵、浆液箱、搅拌器、控制开关的阀门和连接各个组成部分的管道等；通过建立上述火力发电厂脱硫系统液固两相流的等效物理模型，将火力发电厂脱硫系统液固两相流现场的运行情况实时地反映在所述等效物理模型上，可以实时计算脱硫系统各环节浆液情况和脱硫效果，模拟火力发电厂脱硫系统机组的实际运行情况，并为工程技术人员提供无损实验数据。\n[0009] 步骤S2，选择对火力发电厂脱硫系统液固两相流进行有效控制的工艺操作参数，设为等效模型的输入变量；\n[0010] 在步骤S2中，所述等效模型的输入变量可以是浆液箱内的石膏质量浓度、通过脱硫系统的原烟气总体积等。\n[0011] 步骤S3，根据所述输入变量和预设脱硫效率，计算浆液槽内平衡氢离子浓度；\n[0012] 在步骤S3中，根据步骤S2的输入变量和预设脱硫效率，计算浆液槽内平衡氢离子浓度。在本步骤中，可以将浆液区看作全混流反应器，浆液槽内的浆液具有均一的浓度和温度分布，二氧化硫在石灰石浆液中吸收反应的反应机理如下所示：\n[0013]\n[0014]\n[0015]\n[0016] 上式中的平衡常数分别为：Ks＝4.25E39，Ke＝7.69E9，Ksp＝12692.5。\n[0017] 根据上述二氧化硫在石灰石浆液中吸收反应的反应机理，可知，浆液中氢离子浓度的增量等于浆液槽中生成的氢离子浓度减去在反映过程中消耗的氢离子浓度。因此，浆液槽内平衡氢离子浓度可以通过下式进行计算：\n[0018]\n[0019] 其中， 表示浆液中平衡氢离子浓度， 表示从烟气中脱除的二氧化硫体积量， 表示CaCO3的消耗量， 表示CaCO3的摩尔分数，VC浆液槽中浆液体积，V表示原烟气的总体积， 表示脱硫效率，yin表示进口原烟气进口二氧化硫的体积分数，Wlm表示含石灰石浆液量，mlm表示含石灰石浆液质量浓度，mC表示浆液槽中浆液质量，ρc表示浆液槽中浆液密度。\n[0020] 步骤S4，根据所述浆液槽内平衡氢离子浓度，计算循环浆液中二氧化硫摩尔分数；\n[0021] 在步骤S4中，根据步骤S3计算得到的浆液槽内平衡氢离子浓度，计算循环浆液中二氧化硫的摩尔分数，可以采用如下公式：\n[0022]\n[0023] 其中，xin表示循环浆液中二氧化硫的摩尔分数，θ表示浆液中剩余二氧化硫的摩尔分数与平衡氢离子摩尔分数的平衡系数，MC表示浆液的摩尔质量。\n[0024] 步骤S5，根据所述循环浆液中二氧化硫摩尔分数，计算脱硫系统液固两相流的脱硫效率；\n[0025] 在步骤S5中，根据步骤S4计算得到的循环浆液中二氧化硫摩尔分数，计算脱硫系统液固两相流的脱硫效率。\n[0026] 步骤S6，判断所述脱硫效率是否和预设脱硫效率一致；若是，执行步骤S7；若否，将所述脱硫效率替换预设脱硫效率，重复执行步骤S3～S6，直到所述脱硫效率达到预设脱硫效率为止，执行步骤S7；\n[0027] 步骤S7，将所述等效物理模型设为火力发电厂脱硫系统液固两相流的仿真模型，结束仿真操作。\n[0028] 上述火力发电厂脱硫系统液固两相流仿真建模方法，根据进入吸收塔的烟气中的二氧化硫与加入吸收塔内的石灰石在石灰石浆液中吸收反应的反应机理，计算浆液槽内平衡氢离子浓度；然后，根据浆液槽内平衡氢离子浓度和预设脱硫效率，计算得到循环浆液中二氧化硫摩尔分数；最后，根据进口原烟气进口二氧化硫的体积分数和出口净烟气二氧化硫的体积分数，计算出脱硫系统液固两相流的脱硫效率；若所述脱硫效率和预设脱硫效率不一致，则将所述脱硫效率替换预设脱硫效率，重复计算脱硫系统的脱硫效率的步骤，直到所述脱硫效率与预设脱硫效率一致为止，然后将所述等效物理模型设为火力发电厂脱硫系统液固两相流的仿真模型，结束仿真操作；所述火力发电厂脱硫系统液固两相流仿真建模方法，通过对火力发电厂脱硫系统中浆液槽内平衡氢离子浓度、循环浆液中二氧化硫摩尔分数以及脱硫效率的数据进行实时计算，使得工程技术人员能够得到无损实验数据，便于掌握火力发电厂脱硫系统内部的实时运行情况，也为提高脱硫系统的脱硫效率提供参考依据。\n附图说明\n[0029] 图1为本发明的其中一个实施例的火力发电厂脱硫系统液固两相流仿真建模方法流程图；\n[0030] 图2为本发明的另一个实施例的计算火力发电厂脱硫系统脱硫效率的方法流程图；\n[0031] 图3为应用本发明的其中一个实施例的火力发电厂脱硫系统液固两相流仿真建模方法对火力发电厂脱硫系统进行实时控制的方法流程图。\n具体实施方式\n[0032] 为了更进一步阐述本发明所采取的技术手段及取得的效果，下面结合附图及较佳实施例，对本发明的技术方案，进行清楚和完整的描述。\n[0033] 如图1所示，图1为本发明的其中一个实施例的火力发电厂脱硫系统液固两相流仿真建模方法流程图；\n[0034] 一种火力发电厂脱硫系统液固两相流仿真建模方法，包括如下步骤：\n[0035] 步骤S1，建立火力发电厂脱硫系统液固两相流的等效物理模型；\n[0036] 步骤S2，选择对火力发电厂脱硫系统液固两相流进行有效控制的工艺操作参数，设为等效模型的输入变量；\n[0037] 步骤S3，根据所述输入变量和预设脱硫效率，计算浆液槽内平衡氢离子浓度；\n[0038] 步骤S4，根据所述浆液槽内平衡氢离子浓度，计算循环浆液中二氧化硫摩尔分数；\n[0039] 步骤S5，根据所述循环浆液中二氧化硫摩尔分数，计算脱硫系统液固两相流的脱硫效率；\n[0040] 步骤S6，判断所述脱硫效率是否和预设脱硫效率一致；若是，执行步骤S7；若否，将所述脱硫效率置为预设脱硫效率，重复步骤S3～S6，直到所述脱硫效率达到预设脱硫效率为止，执行步骤S7；\n[0041] 步骤S7，将所述等效物理模型设为火力发电厂脱硫系统液固两相流的仿真模型，结束仿真操作。\n[0042] 通过上述步骤，本发明的火力发电厂脱硫系统液固两相流仿真建模方法，通过为工程技术人员提供无损实验数据，准确地反映火力发电厂脱硫系统的实时运行情况，更为火力发电厂脱硫系统的完善提供参考依据。\n[0043] 在其中一个实施例中，本发明的脱硫系统液固两相流仿真建模方法，所述根据所述输入变量和预设脱硫效率，计算浆液槽内平衡氢离子浓度，可以采用如下的公式进行计算：\n[0044]\n[0045] 其中， 表示浆液中平衡氢离子浓度， 表示从烟气中脱除的二氧化硫体积量， 表示CaCO3的消耗量， 表示CaCO3的摩尔分数，VC浆液槽中浆液体积，V表示原烟气的总体积， 表示脱硫效率，yin表示进口原烟气进口二氧化硫的体积分数，Wlm表示含石灰石浆液量，mlm表示含石灰石浆液质量浓度，mC表示浆液槽中浆液质量，ρc表示浆液槽中浆液密度。\n[0046] 在其中一个实施例中，本发明的脱硫系统液固两相流仿真建模方法，所述根据所述浆液槽内平衡氢离子浓度，计算循环浆液中二氧化硫的摩尔分数，可以采用如下的公式进行计算：\n[0047]\n[0048] 其中，xin表示循环浆液中二氧化硫的摩尔分数，θ表示浆液中剩余二氧化硫的摩尔分数与平衡氢离子摩尔分数的平衡系数，MC表示浆液的摩尔质量。\n[0049] 图2为本发明的另一个实施例的计算火力发电厂脱硫系统脱硫效率的方法流程图，由于气体吸收过程是传质控制过程，因此，气体吸收过程的脱硫效率的计算可以从上层到下层进行。\n[0050] 在其中一个实施例中，本发明的脱硫系统液固两相流仿真建模方法，所述根据所述循环浆液中二氧化硫摩尔分数，计算脱硫系统液固两相流的脱硫效率的步骤包括：\n[0051] S501：计算火力发电厂脱硫系统的传质单元高度；\n[0052] S502：根据所述传质单元高度，计算火力发电厂脱硫系统的传质单元数；\n[0053] S503：利用亨利定律，计算出口净烟气二氧化硫的体积分数；\n[0054] S504：根据进口原烟气进口二氧化硫的体积分数和出口净烟气二氧化硫的体积分数，计算脱硫系统液固两相流的脱硫效率。\n[0055] 在其中一个实施例中，本发明的脱硫系统液固两相流仿真建模方法，所述计算火力发电厂脱硫系统的传质单元高度，可以采用如下的公式进行计算：\n[0056]\n[0057] 其中，HOG表示传质单元高度，KY表示总吸收系数，a表示单位体积塔层提供的有效传质体积系数，Ω表示吸收塔的截面积。\n[0058] 在其中一个实施例中，本发明的脱硫系统液固两相流仿真建模方法，所述根据所述传质单元高度，计算火力发电厂脱硫系统的传质单元数，可以采用如下的公式进行计算：\n[0059]\n[0060] 其中，NOG表示传质单元数，Z表示喷淋层的高度，HOG表示传质单元高度，s表示脱硫系数，m表示相平衡系数，V表示原烟气的总体积，yin表示进口原烟气进口二氧化硫的体积分数，yout表示出口净烟气二氧化硫的体积分数，E表示亨利系数，p表示总压，β表示修正系数。\n[0061] 在其中一个实施例中，本发明的脱硫系统液固两相流仿真建模方法，所述利用亨利定律，计算出口净烟气二氧化硫的体积分数，可以采用如下的公式进行计算：\n[0062]\n[0063] 其中， 表示溶质在气相的摩尔分数，m表示相平衡系数，xin表示循环浆液中的摩尔分数，L表示浆液的摩尔浓度，xout表示喷淋层后循环浆液中的摩尔分数，V表示原烟气的总体积， 表示脱硫效率，yin表示进口原烟气进口二氧化硫的体积分数。\n[0064] 在其中一个实施例中，本发明的脱硫系统液固两相流仿真建模方法，所述进口原烟气进口二氧化硫的体积分数和出口净烟气二氧化硫的体积分数，计算脱硫系统液固两相流的脱硫效率，可以采用如下的公式进行计算：\n[0065]\n[0066] 式中， 表示脱硫效率，yin表示进口原烟气进口二氧化硫的体积分数，yout表示出口净烟气二氧化硫的体积分数。\n[0067] 在其中一个实施例中，本发明的脱硫系统液固两相流仿真建模方法，在所述根据所述输入变量和预设脱硫效率，计算浆液槽内平衡氢离子浓度的步骤之前，还可以包括：计算脱硫系统制浆区浆液箱内的石膏密度。\n[0068] 在其中一个实施例中，本发明的脱硫系统液固两相流仿真建模方法，所述计算脱硫系统制浆区浆液箱内的石膏密度，可以采用如下公式进行计算：\n[0069] RO＝KA*KB/(KB-RATIO*(KB-KA))；\n[0070] 其中，RO表示浆液箱内的石膏密度，RATIO表示浆液箱内的石膏质量浓度，KA表示液相密度，KB表示固相密度。\n[0071] 如图3所示，图3为应用本发明的其中一个实施例的火力发电厂脱硫系统液固两相流仿真建模方法对火力发电厂脱硫系统进行实时控制的方法流程图。本实施例的火力发电厂脱硫系统液固两相流仿真建模方法中对脱硫系统的模拟包括脱硫工艺的三种状态：运行状态、暂停状态和终止状态，而控制这三种状态之间的转换可以通过如下几种按钮：启动按钮、停止按钮、暂停按钮和继续按钮。\n[0072] 因此，这三种状态的转换关系如下所示：\n[0073] 1)终止状态――(启动)――运行状态；\n[0074] 2)暂停状态――(继续)――运行状态；\n[0075] 3)运行状态――(暂停)――暂停状态；\n[0076] 4)运行状态――(停止)――终止状态；\n[0077] 5)暂停状态――(停止)――终止状态。\n[0078] 通过上述状态之间的转换，可以模拟火力发电厂脱硫系统液固两相流的实时运行状况，能够为工程技术人员提供无损实验数据，准确地反映火力发电厂脱硫系统的实时运行情况，更为火力发电厂脱硫系统的改进提供参考依据。只有实时掌控火力发电厂脱硫系统内部的运行情况，才能有效控制火力发电厂脱硫系统的脱硫效率，才能有效地降低煤燃烧排放的烟气对环境的污染。\n[0079] 以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。\n[0080] 以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。