密封点挥发性有机物（VOCs）排放速率的数值模拟及实测应用

密封点挥发性有机物（VOCs）排放速率的数值模拟及实测应用摘要：目前对设备与管线组件密封点挥发性有机物（VOCs）排放速率的讨论尚停留阅历估算的阶段，本文以连接件密封点为典型代表，对乙醇生产装置中所涉密封点VOCs排放速率进行数值模拟和实测应用讨论。利用直接测量包扎法对连接件进行包扎密封，将其所形成的包袋内流通空间类比为石油化工反应器，引入化工单元混合时间和空时的概念，通过计算流体力学数值模拟的方法，对比包袋内速度场、浓度场和压力场的分布状况，定量分析出包扎法系统中惰性载气吹扫流量和包袋体积之间的关系，优选合理的操作参数，现场实测并推算出连接件VOCs排放速率。讨论结果表明：适当减小包袋体积可增加混合强度，缩短混合时间；在保持袋内微正压条件下，合理选择吹扫流量以匹配不同的包袋体积，吹扫时间大于袋内混合空时的2倍可基本达到近稳态；建立新的连接件密封点VOCs泄漏浓度和泄漏排放速率之间的幂函数关系式，用以评估核算VOCs年排放量更精确     。挥发性有机物（VOCs）作为二次气溶胶的重要参加物，导致了大气污染物PM2.5、O3的生成。VOCs来自人为和自然 ，工业源是主要人为排放源之一，而石化行业排放又占据肯定重要地位[1]。据美国环保署（EPA）估算，2000年美国石化设备泄漏产生70367t/a，占石化企业无组织VOCs总排放量的20%以上[2]。据国际石油工业环保协会估算，石化企业生产过程中工艺设备管线物料泄漏导致的VOCs排放量占全厂VOCs无组织排放量总量的40%～50%[3]。泄漏检测与修复技术（LDAR）是利用便携式监测设备，检测化工企业各类设备与管线组件上阀门、法兰及连接件等易产生泄漏的密封点，并修复超过肯定浓度的泄漏点，从而达到掌握原料泄漏损失和环境爱护的目的。随着排污许可证管理和环境爱护税的快速推动，VOCs排放核算的精确     性渐渐成为企业和政府关怀的焦点[4]。对石化行业VOCs进行收税的前提是评估核算出VOCs年排放量，目前多采纳单个泄漏点排放浓度的阅历公式法[5]，主要套用EPA、欧盟等国外石化无组织排放因子对石化企业装置组件的无组织排放进行估算，考虑到各国装置操作参数、运行条件、使用原料、运行管理等方面差异，这种借用国外无组织排放因子的做法无助于精确     核算VOCs排放量[6]。石化工业排污许可中对设备与管线组件密封点泄漏挥发性有机物排放量的填报则是结合国内装置密封点的统计状况及中国石化集团公司、中国石油集团公司及中国海油集团公司完成的第一轮LDAR的检测结果，推算出有机气体和挥发性有机液体流经的设备与管线组件总数的0.3%，以泄漏浓度大于等于10000μmol/mol对于各种密封点的定值泄漏系数计算泄漏挥发性有机物的量，认为与企业实际排放量基本吻合[7]。从目前看，对设备与管线组件密封点VOCs排放速率的讨论尚停留在数学统计和阅历估算的阶段，对组件泄漏VOCs排放过程及排放量尚缺乏实质性了解和深化讨论。本讨论通过计算流体力学数值模拟的方法，选取石化行业乙醇生产装置的连接件密封，定量分析VOCs排放速率，直接测量包扎法系统中惰性载气吹扫流量和包袋体积之间的关系，摸索直接测量法实测系统操作规律，利用包扎法进行现场试验，结合LDAR相关方程推导出VOCs泄漏排放速率，建立一套设备与管线组件密封点VOCs排放速率的系统讨论方法，为VOCs真实排放量的核算供应方法和参考依据。后续可结合数值模拟法和现场包扎法对设备与管线组件其他密封点类型进行VOCs泄漏讨论，特殊是结构简单的密封类型，如阀门、泵，可充分发挥数值模拟的技术特点和优势，并借鉴化工单元合适的参数来优化设计现场试验方案，推导各类密封点VOCs排放速率的核算公式，猎取更精确     的年排放量，为建立切合我国石化行业的VOCs排放因子计算标准和评估体系奠定基础。1泄漏排放速率的评估方法1.1评估方法石化行业泄漏损失评估的方法主要是依据EPA泄漏损失评估协议，包括直接测量法、LDAR相关方程法、筛分法和平均排放系数法[8,9]。这四种方法中，直接测量法是通过直接测量泄漏气量，收集并分析泄漏气体组成，其精确     性和牢靠性最高。但由于石化企业的设备与管线组件过多，直接测量法的现场工作量巨大，实际操作性不强，应用较多的是LDAR相关方程法。该方法基于LDAR定量检测每一受控密封点泄漏浓度值，将其净值SV（μmol/mol）代入相应关系方程得出排放速率eVOCs（kg/h）来估算年排放量[10]。相关方程法是将泄漏浓度值利用数学公式转化为排放速率，中间避开了时间变量，精确     性必定受到各装置操作参数和运行介质、设备和管线组件的密封材质以及设备运行管理等各个方面的影响，应依据石化企业的实际状况利用直接测量法加以修正，重点考虑我国石化生产装置的本地化特征，以削减评估核算误差。1.2包扎法直接测量法一般采纳包扎法，将排放密封点或排放口用薄膜包扎起来形成一个密封包袋，让已知流量的惰性载气通入，待载气达到平衡后，收集气体样品后测量总有机物浓度，也可以针对气样中的单个化合物浓度进行分析后加和，通过样品浓度和载气流量推算VOCs的排放速率。包扎法操作系统流程示意图见图1。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359346347.jpg”alt=“1.jpg”width=“915”height=“520”/设备和管线组件的密封点可分为动静密封点，共9大类[2]，因连接件结构规整，将连接密封点作为目标对象。惰性载气选用N2，现场用高纯氮气经减压后进入包扎系统。为削减材质对VOCs成分的影响，全部连接管均采纳特氟龙管。包袋前后分别设置一个流量计，如前后流量基本全都，说明包扎严密好。N2经三通阀分两侧进入袋内，便于对泄漏VOCs气体的混合和携带。通过稳定性检测推断N2和VOCs混合程度，现场选用便携式有机物检测仪进行VOCs浓度测量，数据变化在±5%时，认定袋内气体混合匀称、稳定[6]。测定袋内气体温度后进行采样，将混合气体充入空白采样袋中，样品即可送色谱分析，减去背景值得到净浓度值。1.3排放速率的计算依据物料衡算，包扎法估算密封点VOCs泄漏的排放速率如式(1)。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359341891.jpg”alt=“2.jpg”width=“668”height=“47”/式中，气体常数转换系数为1.219×10-5(K·kmol)/m3；体积浓度分数转换系数为106μmol/mol。2计算流体力学分析2.1数值模拟建模计算流体力学数值模拟是采纳有限体积法对多组分物质在流场空间内的质量守恒方程和动量守恒方程进行离散求解一种数学方法。数值模拟可为现场试验供应有效的数据关心，优化试验方案。依据式(1)，包扎法系统中仅吹扫流量是自变量，其他则是化工系统本身参数，和包扎操作无关。另外，包袋体积的大小与包扎方法、设备和管线组件密封点类型和大小有关，是包扎系统的一个操作变量。连接件包扎后外形规章，类似化工中的套管反应器，连接件在中部，VOCs向上发生泄漏。在同侧，N2分两支进入包袋，呈V字形。N2与VOCs充分混合后从对侧尾部出口流出，数值模拟三维建模的物理结构及相关设置见图2。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359352546.jpg”alt=“2.jpg”width=“915”height=“285”/对不同包袋体积均采纳三维六面体结构化网格建模，如图2所示。管道外径d1=20mm，包袋直径d2分别为24mm/30mm/40mm。内部流场中设三个监测点位：P1在上游N2进口对侧、P2在中部VOCs泄漏对侧、P3在下游出口处。2.2不同包袋体积的混合时间混合时间是进料在全反应器范围内分布匀称实际上所需要的时间，一般采纳示踪剂脉冲注入反应器后达到最终稳定±5%的时间即为试验测定的混合时间tm[12]。将VOCs类比为示踪剂，数值模拟过程中采纳物质输运方程求解，通过VOCs的对流集中方程计算其在混合气体中的摩尔分数[13]。通过非稳态时间计算，监测出口P3点VOCs浓度随时间的变化。在吹扫流量Q=6.03×10-2m3/h的条件下，不同包袋体积VOCs浓度随时间变化趋势及其混合时间见图3，其中对混合浓度进行无因次处理表示为混合程度。参考抱负反应器的设计方程，当量直径de=d2-d1，空时τ=V/Q，不同包袋体积对应的相关流体力学参数见表1。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359355377.jpg”alt=““width=“915”height=“556”/<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359357080.jpg”alt=““width=“912”height=“179”/由图3和表1可见，随着包袋体积的增加，VOCs在袋内混合时间加长。若紧贴近管道进行包扎，形成的包袋体乐观小，混合时间tm仅为1.3s。同样吹扫流量条件下，包袋体积比为6.8（9.42/1.38）倍时，混合时间的倍数略高为7.4（9.6/1.3）倍。随体积的增加，混合空间扩大，返混程度加剧，混合匀称需要更长的时间。空时是抱负反应器宏观混合状况的反映，宏观流场混合匀称后，不同的物料组分需要在微观尺度进一步对流集中才可达到最终的混合匀称。表1数据表明，混合时间均大于空时，而且体积越大，倍数越高，呈非线性增加。EPA建议包袋体积尽可能小，密封空间减小，可保证混合时间最短，且由于停留时间和传热表面积削减可削减或防止内部重组分有机化合物的凝聚[11]。由于现场设备与管线组件的简单状况，特殊是对于动密封点、阀门类简单结构和大设备的包扎，很难保证包袋体乐观小。如包袋体积较大，可适当增加吹扫流量，以降低空时保证宏观混合强度，并延长吹扫时间直到满意混合时间的要求，吹扫时间在空时的两倍以上基本可达到近稳态。2.3不同物料的混合状况考虑到现场很难具备紧贴设备和管线组件进行抱负包扎的实施条件，选择中体积为后续讨论对象，考察包袋内不同点位的物料浓度动态变化状况，O2和N2的体积分数随时间的变化规律分别见图4(a)和(b)。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359352664.jpg”alt=““width=“738”height=“516”/<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359352859.jpg”alt=““width=“792”height=“580”/包扎完成后，充入N2前，包袋内初始状态为空气，O2和N2的体积分数分别占21%和79%，VOCs为0。假设包扎密封好，不存在空气进入袋内的状况，氧气则为完全灭失气体，如图4(a)所示。其浓度在各个点位逐次被氮气吹扫出，浓度缓慢降低，各个点位的浓度最终逐步衰减为0。P1～P3在包袋内从上游到下游逐次分布，O2灭失也是依此挨次逐步衰减。如图4(b)所示，N2为主流气体，浓度逐次、逐步增大，挨次和O2灭失挨次全都。由于N2从VOCs泄漏的上游吹入，VOCs很难进入上游和其对侧绕流区域，这些区域基本被N2渐渐布满，N2体积分数接近100%。在下游区域，N2和VOCs进行混合、携带后流出，由于VOCs占有肯定比例，所以P3点的N2体积分数明显低于P1和P2点。2.4速度场、浓度场和压力场对比分析在中体积包袋中，在吹扫流量Q=6.03×10-2m3/h的条件下，选VOCs泄漏下游中间截面，速度场、浓度场（VOCs）和压力场（总压）如图5所示。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359366716.jpg”alt=“3.jpg”width=“917”height=“300”/如图5(a)速度场所示，因出口位于下方，中下游上部区域速度开头锐减，由于出口缩径的影响，速度在接近出口区域逐步增加；图5(b)压力场显示整个包袋内总压力约在800Pa，分布较匀称，为微正压状态。由于动压头的影响，在出口处总压力略高；图5(c)VOCs的浓度场显示，在中下游上部区域VOCs浓度相对较高，逐次往出口移动。对比三场，可以发觉速度场和压力场相关性极高，而浓度场主要和VOCs泄漏源的位置以及吹扫气体的出入口位置有关。在包扎法实测操作中，吹扫气体的出入口位置设置至关重要。2.5不同吹扫流量的VOCs浓度变化对中体积包袋，将吹扫流量依次增大4倍，即Q2=4×Q1、Q3=16×Q1，在这3个流量条件下，对应的VOCs浓度在P3点随时间的变化趋势见图6，为便利对比，对数据进行4倍处理，3条曲线数值分别为GC1、GC2×4、GC3×16。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359369556.jpg”alt=“5.jpg”width=“918”height=“540”/泄漏VOCs可视为连续型示踪剂，初始状态时浓度为0，泄漏开头后浓度渐渐增加，如图6所示。随着吹扫流量增至Q2，包袋内的混合强度明显增加，混合时间快速减小。流量进一步增至Q3后，袋内消失剧烈湍流，造成出口P3点区域的浓度脉动。依据式(1)，吹扫流量增大，相当于稀释了样品VOCs浓度，GC值则对应降低4倍。由图6可见，4倍处理后最终浓度值并不相同，而是逐步递减，流量增加后反而简单造成排放速率eVOCs的核算偏低。这主要是由于流量增大，造成袋内湍流脉动增加，袋内压力过高，转变袋内气体的密度。EPA特殊强调袋内压力不应超过6894Pa[11]，压力过高，简单造成包扎材料的破损和对包扎密封的破坏。泄漏点的正常泄漏状态时表压为0，假如包袋内压力过高，转变了包扎前密封点原始泄漏状态，简单造成实测数据失真。3包扎法实测分析在石化行业乙醇生产装置现场进行包扎法直接测量分析，密封点的泄漏值SV范围为0～8×103μmol/mol，可涵盖泄漏值2×103μmol/mol的认定标准[14]。现场实测中重复性有：①现场密封点的包扎次数；②每次包扎密封点的采样次数；③一个采样袋的测试次数。各组重复次数均为3，基本可保证明测数据的重复概率。为和LDAR数据及其相关方程曲线对比，选用幂函数形式拟合，详细拟合曲线及相关数据见图7，因相关方程法的