废有机溶剂再生利用企业输转过程VOCs排放模型的建立

牛茹轩, 吴川福, 杨玉飞, 黄泽春, 汪群慧. 废有机溶剂再生利用企业输转过程VOCs排放模型的建立[J]. 环境工程学报, 2020, 14(12): 3495-3504. doi: 10.12030/j.cjee.202002034
引用本文: 牛茹轩, 吴川福, 杨玉飞, 黄泽春, 汪群慧. 废有机溶剂再生利用企业输转过程VOCs排放模型的建立[J]. 环境工程学报, 2020, 14(12): 3495-3504. doi: 10.12030/j.cjee.202002034
NIU Ruxuan, WU Chuanfu, YANG Yufei, HUANG Zechun, WANG Qunhui. Establishment of VOCs emission model during the transfer process of waste organic solvent recycling enterprise[J]. Chinese Journal of Environmental Engineering, 2020, 14(12): 3495-3504. doi: 10.12030/j.cjee.202002034
Citation: NIU Ruxuan, WU Chuanfu, YANG Yufei, HUANG Zechun, WANG Qunhui. Establishment of VOCs emission model during the transfer process of waste organic solvent recycling enterprise[J]. Chinese Journal of Environmental Engineering, 2020, 14(12): 3495-3504. doi: 10.12030/j.cjee.202002034

废有机溶剂再生利用企业输转过程VOCs排放模型的建立

    作者简介: 牛茹轩(1987—),女,博士研究生。研究方向:VOCs排放与固体废物管理。E-mail:nrx1987@163.com
    通讯作者: 汪群慧(1959—),女,博士,教授。研究方向:固体废物资源化。E-mail:wangqh59@sina.com
  • 基金项目:
    国家重点研发计划(2018YFC1900103)
  • 中图分类号: X701

Establishment of VOCs emission model during the transfer process of waste organic solvent recycling enterprise

    Corresponding author: WANG Qunhui, wangqh59@sina.com
  • 摘要: 为提出废有机溶剂再生利用过程中可挥发性有机物(VOCs)减排的合理化建议,通过建立3种废有机溶剂再生利用企业的VOCs排放模型,估算了18种常见废有机溶剂的排放系数。结果表明:排放系数的顺序为小型企业>中型企业>大型企业;小型再生利用企业的VOCs排放系数较高,为9.04~12.59 g·kg−1, 中型再生利用企业的VOCs排放系数为8.89~10.69 g·kg−1,大型再生利用企业的VOCs排放系数为0.09~3.04 g·kg−1。中、小型企业在装料管残留的环节排放系数最高;大型企业在装料过程排放系数最高;温度对VOCs的排放有影响,对大型企业影响较大。综合上述结果,建议小型企业采用密封管道输送物料且在进料与装料口加装VOCs吸收装置,大型企业的操作环境以及储罐保持在较低温度处。上述研究结果可以为VOCs排放控制提供参考。
  • 电化学厌氧消化(electrochemical anaerobic digestion,EAD)是一种利用微生物电解池(microbial electrolysis cell,MEC)的耦合厌氧消化(anaerobic digestion,AD)的新技术,其原理是利用MEC的制氢功能,对系统内的二氧化碳(CO2)进行固定和转化,生成更多的CH4,从而减少CO2的排放[1],并极大地提高了CH4的转化率和沼气的品质[2],所以,EAD作为一种新型的沼气制备技术,因其具有甲烷转化率高、能量回收率高、能耗低、适用于多种有机废弃物和废水的降解等优势而受到了广泛的关注,是目前AD领域的最新发现[3-5].

    但不管是AD,还是EAD,绝大多数相关研究仍然停留在提高消化能力或者副产物的产率等方面[6-8],很少涉及代谢过程的研究. 代谢通量分析(metabolic flux analysis,MFA)是一项研究胞内代谢状态和分析系统代谢能力的强有力技术,通过MFA:(1)可以确定细胞或者系统内存在的不同代谢途径,明确物质的流向和流量[8];(2)可以计算胞内的中间产物的通量,对未知途径加以辨别[9];(3)可以识别代谢节点的刚柔性[10],和对环境扰动的影响做出评价,从而对代谢途径中具有特殊地位的途径实施有效调控[11](如优化目标产物代谢途径、计算最大理论转化率等[12]). 因此,本文将利用MFA中的通过通量平衡分析(flux balance analysis,FBA)的方法,构建EAD代谢网络,研究EAD代谢网络中的相关代谢途径的通量分布信息.

    在代谢通量分析中,通量是反映各代谢产物在整个代谢网络中参与程度的重要参数[7],可以依据通量的大小评价各代谢途径在整个代谢网络中所发挥的作用[13]. 然而由于EAD代谢网络中的相关代谢途径,往往是在相应的酶催化下进行的,在EAD代谢途径中,乙酸的形成和转化是EAD产甲烷过程中最为关键的节点. 乙酸的形成主要的途径是乙酰辅酶A通过磷酸转乙酰酶(PTA)转化为乙酰基-P,然后再通过乙酸激酶(AK)从乙酰基-P转化为乙酰基-P,AK催化反应产生ATP[14]. 厌氧消化EMP途径所产生的丙酮酸转化为乙醛和二氧化碳,然后乙醛通过ADH还原成乙醇[15]. 丁酸的产生则与磷酸转丁酰酶(PTB)与丁酸激酶相关. 此外,ADP/ATP转运酶也发挥着重要的作用. ADP/ATP转运酶又被称为腺嘌呤核苷酸易位子或ADP/ATP载体蛋白,ADP/ATP转运酶通过ADP 与ATP相互转化产生能量从而为细胞提供能量[16]. 而在乙酸转化中,CoF420参与产甲烷菌中的氧化还原反应,充当电子载体. CoF430是甲基辅酶M还原酶的辅助因子,在甲烷生成的最后一步释放甲烷. 而辅酶M是古菌产甲烷菌代谢中甲基转移反应所需的辅酶.

    综上所述,利用FBA不仅可以对EAD系统的代谢过程进行表征,还能对代谢通量进行分析和直观地了解代谢网络中不同代谢途径的产物生成情况以及代谢关键节点处的酶活性变化水平.

    EAD实验装置由 电解电流及电压记录单元、 EAD反应单元和排水集气单元组成,EAD反应单元主要包含400 mL的发酵罐,其中电极距离为3.0 cm,一端的电极片完全地浸于发酵液中,另一端通过导线与宽屏无纸记录仪(SIN-R7000A,中国)相连;在EAD系统中,阳极为石墨电极,阴极为铂电极,参比电极为饱和的Ag/AgCl电极. EAD实验装置如图1所示.

    图 1  EAD实验装置
    Figure 1.  EAD experimental setup
    a. 电解电流及电压记录单元,b. EAD反应单元,c. 排水集气单元
    a. electrolysis current and voltage recording unit, b. EAD reaction unit, c. drainage gas gathering unit.

    厌氧活性污泥来自云南师范大学昆明实验室通过猪粪厌氧发酵. 经检测,pH值为8.07,TS为17.11%,VS为10.08%.

    缓冲营养液配比:NaH2PO4∙2H2O(5.54 g·L−1)、Na2HPO4∙12H2O(23.09 g·L−1)、KCl(0.26 g·L−1)和NH4Cl(0.62 g·L−1).

    微量元素液配比:Na2WO4∙2H2O(0.03 g·L−1)、NaCl(1 g·L−1)、FeCl2∙7H2O(0.07 g·L−1)、CuSO4∙5H2O(0.01 g·L−1)、NiCl2∙6H2O(0.02 g·L−1)、C6H9NO6(1.50 g·L−1)、MgSO4(3.00 g·L−1)、AlK(SO42∙12H2O(0.01 g·L−1)、ZnCl2(0.13 g·L−1)、CaCl2∙2H2O(0.10 g·L−1)、H3BO3(0.01 g·L−1)、MnCl2∙4H2O(0.60 g·L−1)、CoCl2∙6H2O(0.10 g·L−1)、Na2MoO4(0.03 g·L−1)和复合维生素(1 粒·L−1).

    EAD实验设置3个实验组,每组3个平行. 首先在每个实验组的发酵罐中添加:葡萄糖1.0 g,接种物120.0 g,缓冲营养液100.0 mL,微量元素液10.0 mL,再补加蒸馏水使物料的总质量为360.0 g;其次将密封好的发酵罐连接气体收集装置,再利用高纯氮排出体系内的溶解氧和空气(通气时间为3 min),使EAD体系处于厌氧状态;最后施加恒定电压1.0 V,进行培养.

    培养阶段:将发酵罐置于恒温磁力搅拌器中,设置温度为30.0 ℃;每天监测pH值变化,待其恢复到7.0左右,立即添加等量的葡萄糖将体系培养至稳定.

    扰动阶段:培养阶段结束后,测定发酵液中底物和代谢产物的浓度,同时记下该时间点为扰动的初始点,随后添加1.0 g的葡萄糖,施加扰动后,通过监测体系中底物的消耗量和代谢产物的生成量,然后进行代谢通量分析.

    日产气量:采用排水法法,每天定时记录产气量. TS和VS:在(105±5)℃下干燥并在(550±10)℃下马弗炉燃烧1 h后,测量污泥的. pH值:采用pH计(HAOSHI H-1008A,中国).

    VFAs:样品经离心分离,取上清液,采用气相色谱仪(FULI, GC9790II,中国),以FID为检测器进行测量[17]. 色谱柱为KB- FFAP(30 m×0. 32 mm×0. 25 μm)毛细管柱.

    气体含量:采用气相色谱仪(FULI, GC9790II,中国),以TCD 检测器对H2、CH4和CO2含量进行测量. 色谱柱为TDX-01填充柱. 葡萄糖含量:采用3,5-二硝基水杨酸比色法测定[18].

    酶活测量采用改良的双抗体夹心酶联免疫吸附测定(ELISA)方法[19-21]在测量过程中,以ADP/ATP转运酶、磷酸转乙酰酶(PTA)、乙酸激酶(AK)、磷酸转丁酰酶(PTB)、丁酸激酶(BK)、辅酶M(CoM)、甲基丙二酰CoA变位酶(MCM)、辅酶F420(F420)、辅酶F430(F430)为标准酶. 检测时在酶标包被板上设定空白孔、标准孔、待测样品孔. 空白孔不加样品及标准酶,作为对照,标准孔加标准酶,待测样品孔中先加酶试剂盒中的样品稀释液,然后再加待测样品. 按照酶试剂盒说明书进行操作,最后采用Epoch酶标仪在450 nm波长下进行测定,以空白孔调零依序测量各孔的吸光度(OD值),并通过标准曲线计算样品中酶的活性.

    通量分析采用通量平衡分析(FBA)的方法,采用代谢通量分析软件Cell Net Analyzer进行代谢网络的构建及分析. FBA是MFA中基于质量守恒的代谢过程研究常用的一种数学建模分析法[22],在拟稳态条件下,模型中代谢产物处于平衡状态,此时,可将代谢物的通量描述为一个线性方程,只需对底物和相关产物的浓度变化进行测量,再采用线性规划或者计量矩阵分析方法则可确定每种产物的代谢通量[23].

    EAD反应过程中产气特性的变化是表征反应器效能的重要因素. 图2为培养过程中的产气变化情况,从图2可以发现,在反应器启动阶段,气体中以CO2为主,随着反应器的运行,产甲烷菌逐渐生长富集. 在稳定阶段,其甲烷含量最高可以达到55%左右,葡萄糖厌氧消化中理论产甲烷含量在50%左右,也就是说,MEC的引入明显提升了葡萄糖的理论产甲烷率,而且此时CO2的含量降低到11%左右,明显的降低了CO2的排放. 在反应器稳定后,日产气量也得到了明显的提升,从最初的50 mL·d−1左右增加到了200 mL·d−1左右. 对于氢气含量,除启动阶段有氢气产生外,反应中后期基本没有氢气产生,这可能是由于前期pH较低,适于产氢菌的生长,中后期pH值恢复较快,可以迅速恢复至已6.5以上,不利于产氢菌的生长,导致氢气含量极低. 此外,也可能是由于电化学辅助促进了H2还原CO2途径产甲烷,使得反应器中的H2迅速被噬氢产甲烷菌消耗,从而提高了CH4的含量,降低了H2的含量.

    图 2  EAD培养阶段中产气体积及气体组分含量
    Figure 2.  Gas production volume and gas component content in EAD operation stage

    VFAs是厌氧消化中重要的中间产物,其产生和分解过程受多种因素影响,同时也关系到EAD产甲烷过程的稳定性. EAD体系中 VFAs的浓度一旦超过 200 mmol·L−1 就有可能发生酸抑制问题[24],并对产甲烷过程产生严重的影响,然而本次实验研究并未出现VFAs累积的现象. 图3显示了EAD培养过程中各短链有机酸的变化情况. 从图3可以看出,各短链有机酸的变化与底物补充时间相关,培养前期VFAs波动较大且浓度较高,但15 d后,因阳极上形成了电活性微生物膜,致使VFAs能够在阳极表面降解,尽管进料时间间隔缩短为3 d,也没有产生VFAs累积,表明EAD具有较强的抗酸化的能力.

    图 3  EAD运行阶段中VFAs的含量变化
    Figure 3.  Changes in the content of VFAs during the EAD operation stage

    在厌氧消化中丁酸累积是抑制厌氧消化的重要因素之一,由于丁酸的产氢产乙酸过程必须依赖氢还原CO2产甲烷过程,否则本身不能自发进行,一旦产氢产乙酸和产甲烷过程失衡,将会产生严重的丁酸抑制作用。然而在EAD中丁酸的降解除了依赖氢还原CO2产甲烷过程外,还可以在电活性微生物的作用下,在电极表面降解. 在这条途径中,丁酸分解为CO2、H+和电子, H+向阴极迁移在阴极上获得电子而形成H2,这部分H2同样被氢营养型产甲烷菌利用还原CO2,加强H2还原CO2产甲烷途径.

    对于纯培养体系而言,由于体系内菌种单一,代谢产物和途径比较明确,构建与之相对应的代谢网络相对比较容易. 但对于复杂的混合培养而言,特别是像EAD这种由多种微生物构成的体系,变量较多,对于代谢网络的构建及分析难度较大[25],为了使代谢通量分析能适用于混合体系,提出了“统一微生物”这一个概念[26],同时国内外一些学者已经在部分混合体系中作了应用上的尝试,并且成功地建立了发酵制氢体系的代谢网络[27],还有在新型多级厌氧甲烷反应器中也进行了FBA的研究[28],但还需进一步地完善,才能使其适用于EAD的混合培养体系中. EAD是在AD中引入MEC,由于MEC的引入,使H+向阴极迁移,在阴极上生成H2,这在AD基础上增加了一条重要的产氢途径. 根据相关研究表明[9, 12],在严格的厌氧条件下,丙酮酸几乎不会流向三羧酸循环,甲酸也不是所有严格厌氧菌都能产生. 结合相关研究[29-33],构建了如图4所示的EAD代谢网络该网络包含了37个主要的代谢反应式,每一步的代谢的具体生化反应过程如表1所示.

    图 4  电压扰动下EAD代谢网络及通量分布
    Figure 4.  EAD metabolic flux under voltage perturbation
    通量值是将表2数据输入CNA软件运行后,以100为基换算得到,数值大小只表示通量大小.
    The flux value is calculated based on 100 after inputting the data in Table 2 into the CNA software and running. The value only represents the flux size.
    表 1  EAD体系葡萄糖代谢过程包含的主要生化反应
    Table 1.  The main biochemical reactions included in the glucose metabolism process of the EAD system
    NO.反应式 EquationNO.反应式 Equation
    R1GLC+PEP ∀ G6P + PyrR20HPr + 2 H2O ∀ HAc + CO2 +3 H2
    R2NADH + Pyr ∀ HLa + NADR21HBu + 2 H2O ∀ 2 HAc + 2 H2
    R3Pyr + NADH ∀ NAD + HFo +AcCoAR22HPr = HPr (ext)
    R42 Fd + Pyr + CoA ∀ CO2 + AcCoA + 2 FdHR23EtOH+ HPr ∀ Hva (ext) + H2O
    R52 Fd + NADH = 2 FdH + NADR24HFo = HFo (ext)
    R6NADPH + NAD ⇌ NADH + NADPR25CO2 = CO2 (ext)
    R72 NADH = H2 + 2 NADR26CO2 + 4 H2 = CH4 + 2 H2O
    R82 FdH = H2 + 2 FdR27H2 = H2 (ext)
    R9HLa = HLa (ext)R282 AcCoA + H2 ∀ PrOH(ext) +2 CoA + CO2
    R10HLa + NADH ∀ HPr + NADR29HAc ∀ CO2 + CH4
    R11HFo ∀ CO2 + H2R30HAc = HAc (ext)
    R12AcCoA + ADP + iP ∀ATP + HAc + CoAR31EtOH = EtOH (ext)
    R13AcCoA + 2 NADH ∀ EtOH + CoA + 2 NADR32HBu + EtOH ∀ HCa (ext) + H2O
    R142 AcCoA + NADH ⇌ CoA + H2O + CroCoA + NADR33HBu = HBu (ext)
    R15CroCoA + 2 Fd + NADH ∀ ButCoA + 2 FdH + NADR34HBu + 2 NADH ∀ BuOH(ext) + CoA + 2 NAD
    R16CroCoA + NADH ∀ ButCoA + NADR35Pyr ⇌ Pyr (ext)
    R17ButCoA + ADP + iP ∀ HBu + ATP + CoAR36CH4 = CH4 (ext)
    R182 CO2 + 4 H2 ∀ HAc +2 H2OR372 H+ + 2 e- = H2
    R19EtOH + 2H2O ∀ 4 H+ + HAc
      注:iP为磷酸根离子,CoA为辅酶A,H+为氢离子,e-为转移电子或者电解协助装置提供的电子.
      Note: iP is the phosphate ion, CoA is the coenzyme A, H+ is the hydrogen ion, and e- is the electron provided by the transfer electron or electrolysis assistance device.
     | Show Table
    DownLoad: CSV

    在培养阶段中,当EAD运行稳定后,从取样口取样,检测残余葡萄糖和相关代谢物的含量,然后添加1 g的葡萄糖,以此时的状态为FBA的初始状态,实施电压扰动. 在对照组电压为1.0 V的基础上,扰动组1电压降低至0.6 V,而扰动组2增加至1.4 V. 14 h后,再次检测残余葡萄糖和相关代谢物,产气量和气体含量. 经过计算得到EAD扰动前后的发酵液中代谢产物的净生成速率,如表2所示.

    表 2  EAD体系中关键代谢物的净生成速率(mmol·L−1·h−1
    Table 2.  Net formation rate of key metabolites in EAD system(mmol·L−1·h−1
    代谢物Metabolites对照组1.0 VControl group 1.0 V扰动组0.6 VDisturbance group 0.6 V扰动组1.4 VDisturbance group 1.4 V
    葡萄糖1.1274 ± 0.01541.0958 ± 0.00551.1195 ± 0.0175
    甲 酸0.0032± 0.00000.0086 ± 0.00210.0093± 0.0026
    乙 酸0.7511 ± 0.14400.5229 ± 0.07700.7064 ± 0.1202
    丙 酸0.1061 ± 0.00960.0860 ± 0.00490.0902 ± 0.0023
    丁 酸0.4180± 0.03350.4218± 0.03910.4026 ± 0.0201
    乳 酸0.0107 ± 0.00460.0081 ± 0.00350.0128 ± 0.0056
    丙酮酸0.0156 ± 0.00370.0335± 0.00740.0261 ± 0.0171
    戊 酸0.0230 ± 0.01070.0107 ± 0.00460.0277 ± 0.0093
    己 酸0.0000 ± 0.00000.0000 ± 0.00000.0000 ± 0.0000
    乙 醇0.0173 ± 0.00140.0092 ± 0.00190.0160 ± 0.0067
    丙 醇0.0000 ± 0.00000.0000 ± 0.00000.0000 ± 0.0000
    丁 醇0.0000 ± 0.00000.0000 ± 0.00000.0000 ± 0.0000
    CO20.2258 ± 0.02290.2255 ± 0.03760.2716 ± 0.0071
    CH40.2947 ± 0.01270.4909 ± 0.05150.3924 ± 0.0293
    H20.0281 ± 0.02310.0228 ± 0.01370.0039 ± 0.0024
      注:表2中的实验数据为电压扰动实验组的主要代谢物的净生成速率;葡萄糖的测量值代表底物的净摄取速率;丁酸为正丁酸和异丁酸的总和,戊酸为正戊酸和异戊酸的总和.
      Note: The experimental data in Table 2 are the net production rates of major metabolites in the voltage disturbance group; The measured value of glucose represents the net uptake rate of substrate; Butyrate is the sum of n-butyrate and isobutyrate, valerate is the sum of n-valerate and isovalerate.
     | Show Table
    DownLoad: CSV

    根据表2中的净生成速率,采用运行Cell Net Analyzer软件,获得EAD代谢通量分布信息,如图4所示,并以H2、CO2、HAc和CH4代谢物进行评价.

    在0.6 V扰动下,NADH通量明显优于1.4 V电压扰动,且此条件下最终产物甲烷通量达到43.52,明显高于1.0 V和1.4 V下的26.14和34.81. 由图4可知,H2的代谢途径主要形成途径为:R7、R20、R21和R37. 由通量值可知,EAD中H2主要来源于电极反应R37,和丁酸降解为乙酸的途径R21,其次是NADH平衡调节产氢R7. CH4为最终的目标代谢物,主要由H2还原CO2途径R26和乙酸转化途径R29两条途径形成. EAD体系中CH4的合成途径是以乙酸转化途径占主导地位,其次为H2还原CO2途径,1.0 V下EAD体系的这两条途径的占CH4产生总量的64.95%和35.05%. 0.6 V和1.4 V扰动下,通过H2还原CO2产生的产甲烷(R26)分别占13.44%和12.76%,相比对照组5.97%均有所提高,表明电压扰动使得EAD产甲烷代谢途径偏向了氢营养型产甲烷. 值得注意的EAD产甲烷系统明显产生了甲酸,但根据代谢通量,CH4和H2并无一来源于甲酸,说明EAD反应体系中可能没有噬甲酸型产甲烷菌存在. 最终产物H2大多由EAD中的阴极产生,但在最终产物中,氢气通量明显下降,说明H2还原为CH4(R26)的途径较为活跃. 本次从发酵液代谢物中未检测到己酸、丙醇和己酸盐. 乙酸的来源有R12、R18、R19、R20和R21途径. 通量值表明乙酰辅酶A转化是乙酸最主要的来源,其次是丁酸转化、H2还原CO2、乙醇转化和丙酸转化. 其通量大小受多个中间产物的影响,但在电压扰动(0.6 V和1.4 V)下,其通量明显增大,其中以0.6 V扰动增幅最大,说明反应器中的产酸微生物在此时相对较活跃. 乳酸和乙醇的产生不利于H2的生产,图中0.6 V扰动下的乳酸通量较高,这也可能是导致了H2较低的原因之一.

    表3为葡萄糖代谢过程中关键酶活性的变化. 其中ADP/ATP转运酶的活性,0.6 V和1.4 V扰动时的酶活性均高于对照组1.0 V 的239.28 IU·L−1. ADP/ATP是能量代谢的载体,为微生物生长和繁殖提供能量基础. 当电压扰动时,需要大量能量来使微生物群落发生转变,因此ADP/ATP转运酶酶活性增加. 在乙酸代谢途径中首要的是乙酸的活化,它需要由乙酸激酶AK和磷酸转乙酰酶PTA的先后催化来进行. 在ATP参与下由乙酸激酶AK催化,乙酸被磷酸化为乙酰磷酸,然后乙酰磷酸再在磷酸转乙酰酶PTA催化下生成乙酰辅酶A,进行乙酸代谢[34-35]. 但从图5中可以看出,所产生的乙酸大部分都R12途径逆反应生成乙酰辅酶A,其中当进行0.6 V扰动时,R12逆反应通量最大,这与在0.6 V时PTA活性达到115.69 U·L−1检测结果相吻合. 当进行电压扰动时,EAD体系内在乙酸代谢中PTA和AK的酶活调节现象相似,有研究表明[36]当谷氨酸棒杆菌生长在以乙酸(取代葡萄糖)为碳源的培养基上时,PTA和AK酶活性明显提高,该现象是微生物碳代谢中葡萄糖效应的具体体现,表明酶在基因表达环节上可能存在葡萄糖基质上的负调控效应或乙酸基质上的正调控效应. 降低电压,PTA与AK酶活性增加. PTB和BK酶与丁酸产生有关,由图5代谢通量可知,丁酸的产生是从乙酰基通过PTB与BK的共同催化所形成. 在丁酸代谢中,PTB催化反应产生大量的丁酰辅酶A,但产生的丁酸大部分却通过R21转化为乙酸. 在0.6 V扰动下,PTB的酶活性352.19 U·L−1,明显高于对照组310.19 U·L−1和1.4 V的338.04 U·L−1,但BK酶活性,对照组1.0 V的酶活性52.59 U·L−1,要略高于0.6 V的50.92 U·L−1. PTB酶活性高于BK,这是由于PTB需要先将丁酰辅酶A转化为在丁酰磷酸,最终再由BK催生与ADP反应形成丁酸和ATP. 而BK在1.0 V与0.6 V酶活性差异不大,这是可能是由于电压BK影响不大. 辅酶M在甲烷形成的最终反应步骤中充当甲基的载体. CoM的化学结构虽然很简单,但在甲基还原酶的反应体系中却有高度的专一性,因此在进行电压扰动时,辅酶M的酶活性变化差异不大. 值得一提是,MCM酶活性变化与CoM酶活性变化相似. 在丙酮酸转化为丙酸代谢过程中,关键酶为MCM,在电压扰动下,酶活性变化差异较小. 结合图6中的产丙酸代谢通量图可以看出,丙酸通量变化也差异较小,说明酶活性与所构建的代谢通量有较好的相关性. CoF420在产甲烷途径中电子载体,而CoF430产甲烷古菌形成甲烷最后步骤作为甲基载体的四吡咯化合物,因此产甲烷的总通量应当取决于CoF420与CoF430共同酶活性,在0.6 V电压扰动下,总酶活性为63.01 IU·L−1,对照组为56.94 IU·L−1,1.4 V扰动为59.57 IU·L−1,这与图5中产甲烷代谢通量中,0.6 V时产甲烷最多,1.4 V次之,对照组最少结果相吻合.

    表 3  代谢途径中关键酶活性水平的变化
    Table 3.  Changes in activity levels of key enzymes in metabolic pathways
    酶名称Enzyme对照组1.0 VControl group 1.0 V扰动组0.6 VDisturbance group 0.6 V扰动组1.4 VDisturbance group 1.4 V单位Unit
    ADP/ATP转运酶239.28±23.54337.04±7.88344.23±7.54IU·L−1
    磷酸转乙酰酶(PTA)99.78±6.39115.69±11.3297.69±7.50U·L−1
    乙酸激酶(AK)106.26±24.23151.38±28.81120.25±28.93U·L−1
    磷酸转丁酰酶(PTB)310.19±1.29352.19±12.66338.04±22.08U·L−1
    丁酸激酶(BK)52.59±6.6350.92±0.9036.08±4.89U·L−1
    辅酶M(CoM)173.84±17.97191.01±23.53176.85±12.82U·L−1
    甲基丙二酰CoA变位酶(MCM)58.17±4.0651.71±2.7154.20±3.09U·L−1
    辅酶F420(CoF42040.38±2.2243.72±4.3736.86±6.48IU·L−1
    辅酶F430(CoF43016.56±2.0519.29±1.5522.71±2.14IU·L−1
     | Show Table
    DownLoad: CSV

    本研究通过外加电压对EAD体系扰动过程中代谢通量的影响,对EAD代谢网络通量上的变化结合关键酶活性的变化,进一步分析电压扰动对EAD系统代谢产物生成上的提高或者降低提供理论依据,由通量分析结果可知,EAD体系中并没有很高的H2生成,相反,EAD体系中产甲烷能力得到提高. 原因可能是由于电解协助装置的引入,极大地促进H2的产生,随后,所产生的H2被用来还原CO2产乙酸和产甲烷的代谢过程,进而提升了EAD系统的产甲烷能力. 结合关键节点处酶活性变化,酶活性与代谢通量有着较好的相关性. 在进行电压扰动时,ADP/ATP转运酶、PTA与AK、PTB与BK,这些酶活性均随电压变化而产生较大的差异. 而CoM与MCM酶活性,具有高度专一性的酶,随电压扰动变化不大. 在0.6 V电压扰动下,CoF420与CoF430酶活性与产甲烷代谢通量变化相吻合.

  • 图 1  小型废有机溶剂再生利用企业工艺流程及VOCs排放节点

    Figure 1.  Process flow and VOCs emission nodes in small-scale waste organic solvent recycling and utilization enterprise

    图 2  中型废有机溶剂再生利用企业工艺流程及VOCs排放节点

    Figure 2.  Process flow and VOCs emission nodes in medium-scale waste organic solvent recycling and utilization enterprise

    图 3  大型废有机溶剂再生利用企业工艺流程及VOCs排放节点

    Figure 3.  Process flow and VOCs emission nodes in large-scale waste organic solvent recycling and utilization enterprise

    图 4  3种废有机溶剂再生利用企业VOCs排放系数

    Figure 4.  VOCs emission coefficients of waste organic solvents in three types of recycling and utilization enterprises

    图 5  3种废有机溶剂再生利用企业VOCs的排放系数比例

    Figure 5.  VOCs emission coefficients ratio of waste organic solvents in three types of recycling and utilization enterprises

    图 6  温度对VOCs排放的影响

    Figure 6.  Effect of temperature on VOCs emissions

    表 1  3种废有机溶剂再生利用企业VOCs排放节点对比

    Table 1.  Comparison of VOCs emission nodes in three types of waste organic solvent utilization enterprises

    VOCs排放节点混装卸料进料进料管残留装料装料管残留
    小型企业
    中型企业
    大型企业
      注:“√”表示存在此排放点。
    VOCs排放节点混装卸料进料进料管残留装料装料管残留
    小型企业
    中型企业
    大型企业
      注:“√”表示存在此排放点。
    下载: 导出CSV

    表 2  装载损耗计算中的饱和因子

    Table 2.  Saturation factor in load loss calculation

    装载方式装载容器饱和因子
    底部/液下新桶罐或清洗后的桶罐0.5
    正常工况(普通)的桶罐0.6
    上次卸料采用蒸汽平衡装置1.0
    喷溅式新桶罐或清洗后的桶罐1.45
    正常工况(普通)的桶罐1.45
    上次卸料采用蒸气平衡装置1.0
      注:数据来源于《VOCs污染源排查指南》。
    装载方式装载容器饱和因子
    底部/液下新桶罐或清洗后的桶罐0.5
    正常工况(普通)的桶罐0.6
    上次卸料采用蒸汽平衡装置1.0
    喷溅式新桶罐或清洗后的桶罐1.45
    正常工况(普通)的桶罐1.45
    上次卸料采用蒸气平衡装置1.0
      注:数据来源于《VOCs污染源排查指南》。
    下载: 导出CSV

    表 3  18种物质蒸汽压

    Table 3.  Vapor Pressures of 18 types of substance

    名称25 ℃蒸汽压/kPa名称25 ℃蒸汽压/kPa名称25 ℃蒸汽压/kPa
    苯乙烯0.80异丙醇4.34正己烷20.17
    邻二甲苯0.88三氯乙烯9.91乙酸乙酯21.23
    间二甲苯1.11均二氯乙烷10.65乙酸甲酯28.46
    对二甲苯1.1712.69丙酮30.67
    乙苯1.27四氯化碳15.37丙醛42.46
    甲苯3.90甲醇16.85二氯甲烷58.12
      注:数据来源于《化学化工物性数据手册》[25]
    名称25 ℃蒸汽压/kPa名称25 ℃蒸汽压/kPa名称25 ℃蒸汽压/kPa
    苯乙烯0.80异丙醇4.34正己烷20.17
    邻二甲苯0.88三氯乙烯9.91乙酸乙酯21.23
    间二甲苯1.11均二氯乙烷10.65乙酸甲酯28.46
    对二甲苯1.1712.69丙酮30.67
    乙苯1.27四氯化碳15.37丙醛42.46
    甲苯3.90甲醇16.85二氯甲烷58.12
      注:数据来源于《化学化工物性数据手册》[25]
    下载: 导出CSV

    表 4  挥发性有机物各类源排放系数的推荐值

    Table 4.  Recommended values of VOCs Emission coefficients from various sources

    1级2级名称排放系数/(g·kg−1)
    工艺过程源基础化学原料制造乙烯0.097
    工艺过程源基础化学原料制造丙烯0.111
    工艺过程源基础化学原料制造丙烯腈0.988
    工艺过程源基础化学原料制造乙苯0.1
    工艺过程源基础化学原料制造丁二烯139.74
    工艺过程源基础化学原料制造苯乙烯0.223
      注:数据来源于《排放清单指南》。
    1级2级名称排放系数/(g·kg−1)
    工艺过程源基础化学原料制造乙烯0.097
    工艺过程源基础化学原料制造丙烯0.111
    工艺过程源基础化学原料制造丙烯腈0.988
    工艺过程源基础化学原料制造乙苯0.1
    工艺过程源基础化学原料制造丁二烯139.74
    工艺过程源基础化学原料制造苯乙烯0.223
      注:数据来源于《排放清单指南》。
    下载: 导出CSV
  • [1] 王扶潘, 朱乔, 冯凝, 等. 深圳大气中VOCs的二次有机气溶胶生成潜势[J]. 中国环境科学, 2014, 34(10): 2449-2457.
    [2] NGUYEN K, DABDUB D. NOx and VOC control and its effects on the formation of aerosols[J]. Aerosol Science and Technology, 2002, 36(5): 560-572. doi: 10.1080/02786820252883801
    [3] SHAO M, ZHANG Y, ZENG L, et al. Ground-level ozone in the Pearl River Delta and the roles of VOC and NOx in its production[J]. Journal of Environmental Management, 2009, 90(1): 512-518. doi: 10.1016/j.jenvman.2007.12.008
    [4] 国务院办公厅. 关于推进大气污染联防联控工作改善区域空气质量的指导意见: 国办发[2010]33号[Z]. 北京, 2010.
    [5] 环境保护部. 大气挥发性有机物源排放清单编制技术指南(试行): 公告 2014年 第55号[Z]. 北京, 2014.
    [6] WEI W, WANG S, CHATANI S, et al. Emission and speciation of non-methane volatile organic compounds from anthropogenic sources in China[J]. Atmospheric Environment, 2008, 42: 4976-4988. doi: 10.1016/j.atmosenv.2008.02.044
    [7] 环境保护部, 发展改革委. 国家危险废物名录: 部令第1号[Z]. 北京, 2016.
    [8] 郑俊. 上海市典型工业区大气中VOCs及其来源分析[J]. 上市师范大学学报(自然科学版), 2017, 46(2): 298-303.
    [9] 何家禧, 周伟, 黄红英. 工业用混合有机溶剂挥发性化学组分分析[J]. 中国职业医学, 2013, 40(3): 256-257.
    [10] 莫梓伟, 陆思华, 李悦, 等. 北京市典型溶剂使用企业VOCs排放成分特征[J]. 中国环境科学, 2015, 35(2): 374-380.
    [11] 田亮, 魏巍, 程水源, 等. 典型有机溶剂使用行业VOCs成分谱及臭氧生成潜势[J]. 安全与环境学报, 2017, 7(1): 314-320.
    [12] 赵锐, 黄络萍, 程军, 等. 成都市工业源重点VOC排放行业排放清单及空间分布特征[J]. 环境科学学报, 2018, 38(4): 1358-1367.
    [13] 陈浩, 林艳发, 钟学情. 2011-2013年深圳市龙岗区工作场所有机溶剂成分监测分析[J]. 实用预防医学, 2013, 22(4): 474-476. doi: 10.3969/j.issn.1006-3110.2013.04.033
    [14] 王肖丽, 金宇帆, 黄志炯, 等. 基于实地调研的广东省工业VOC排放清单改进研究[J]. 环境科学学报, 2019, 39(4): 1113-1124.
    [15] 余翔翔, 杨畅, 赵崇业, 等. 温州市木质家具行业VOCs排放特征及排放系数研究[J]. 现代商贸工业, 2019, 29(24): 194-195.
    [16] 梁小明, 陈来国, 孙西勃, 等. 基于原料类型及末端治理的典型溶剂使用源VOCs排放系数[J]. 环境科学, 2019, 40(10): 4382-4394.
    [17] 李建伟, 刘新宇, 修光利. VOCs无组织排放估算方法和控制标准初探[J]. 化学世界, 2010(10): 632-634. doi: 10.3969/j.issn.0367-6358.2010.10.017
    [18] 刘自力. 废气中有机溶剂的处理及回收[J]. 化学工程与设备, 2017(4): 221-222.
    [19] 牛茹轩, 高明, 蔡洪英, 等. 重庆市废有机溶剂产生特性[J]. 中国环境科学, 2019, 39(3): 1101-1108. doi: 10.3969/j.issn.1000-6923.2019.03.025
    [20] 谢添, 杨文, 郭婷, 等. 化工行业挥发性有机物无组织排放特征研究: 以天津化工企业为例[J]. 南开大学学报(自然科学版), 2017, 50(3): 79-83.
    [21] 苏艳明. 石化工业储罐无组织排放VOCs的定量探讨[J]. 石油化工安全环保技术, 2018, 34(3): 56-60. doi: 10.3969/j.issn.1673-8659.2018.03.018
    [22] 张秀青. 石化企业废气无组织排放源及排放量估算简介[J]. 装备环境工程, 2008, 5(5): 74-77. doi: 10.3969/j.issn.1672-9242.2008.05.020
    [23] 生态环境部. 石化行业VOCs污染源排查工作指南: 环办[2015]104号[Z]. 北京, 2014.
    [24] 王卓. 石化企业VOCs检测及排放量核算研究[D]. 青岛: 中国石油大学, 2016.
    [25] 刘光启, 马连湘, 项曙光. 化学化工物性数据手册: 有机卷[M]. 北京: 化学工业出版社, 2013.
  • 加载中
图( 6) 表( 4)
计量
  • 文章访问数:  3755
  • HTML全文浏览数:  3755
  • PDF下载数:  46
  • 施引文献:  0
出版历程
  • 收稿日期:  2020-02-07
  • 录用日期:  2020-03-29
  • 刊出日期:  2020-12-10
牛茹轩, 吴川福, 杨玉飞, 黄泽春, 汪群慧. 废有机溶剂再生利用企业输转过程VOCs排放模型的建立[J]. 环境工程学报, 2020, 14(12): 3495-3504. doi: 10.12030/j.cjee.202002034
引用本文: 牛茹轩, 吴川福, 杨玉飞, 黄泽春, 汪群慧. 废有机溶剂再生利用企业输转过程VOCs排放模型的建立[J]. 环境工程学报, 2020, 14(12): 3495-3504. doi: 10.12030/j.cjee.202002034
NIU Ruxuan, WU Chuanfu, YANG Yufei, HUANG Zechun, WANG Qunhui. Establishment of VOCs emission model during the transfer process of waste organic solvent recycling enterprise[J]. Chinese Journal of Environmental Engineering, 2020, 14(12): 3495-3504. doi: 10.12030/j.cjee.202002034
Citation: NIU Ruxuan, WU Chuanfu, YANG Yufei, HUANG Zechun, WANG Qunhui. Establishment of VOCs emission model during the transfer process of waste organic solvent recycling enterprise[J]. Chinese Journal of Environmental Engineering, 2020, 14(12): 3495-3504. doi: 10.12030/j.cjee.202002034

废有机溶剂再生利用企业输转过程VOCs排放模型的建立

    通讯作者: 汪群慧(1959—),女,博士,教授。研究方向:固体废物资源化。E-mail:wangqh59@sina.com
    作者简介: 牛茹轩(1987—),女,博士研究生。研究方向:VOCs排放与固体废物管理。E-mail:nrx1987@163.com
  • 1. 北京科技大学能源与环境工程学院,北京 100083
  • 2. 工业典型污染物资源化处理北京市重点实验室,北京 100083
  • 3. 中国环境科学研究院,土壤与固体废物环境研究所,北京 100012
基金项目:
国家重点研发计划(2018YFC1900103)

摘要: 为提出废有机溶剂再生利用过程中可挥发性有机物(VOCs)减排的合理化建议,通过建立3种废有机溶剂再生利用企业的VOCs排放模型,估算了18种常见废有机溶剂的排放系数。结果表明:排放系数的顺序为小型企业>中型企业>大型企业;小型再生利用企业的VOCs排放系数较高,为9.04~12.59 g·kg−1, 中型再生利用企业的VOCs排放系数为8.89~10.69 g·kg−1,大型再生利用企业的VOCs排放系数为0.09~3.04 g·kg−1。中、小型企业在装料管残留的环节排放系数最高;大型企业在装料过程排放系数最高;温度对VOCs的排放有影响,对大型企业影响较大。综合上述结果,建议小型企业采用密封管道输送物料且在进料与装料口加装VOCs吸收装置,大型企业的操作环境以及储罐保持在较低温度处。上述研究结果可以为VOCs排放控制提供参考。

English Abstract

  • 随着全球工业化的快速发展,城市大气中的污染物含量逐渐增多。挥发性有机物(VOCs)作为臭氧和二次有机气溶胶的重要前体物,极易发生光化学反应,形成光化学烟雾[1-3]。VOCs的排放不仅影响我国城市环境空气质量,同时对排放企业职工及周边居民的身体健康造成严重的损害。在2010年国务院发布的《关于推进大气污染联防联控工作改善区域空气质量的指导意见》[4]中,将VOCs列为继SO2、NO2和颗粒物之后拟重点防控的大气污染物。但由于VOCs种类繁多、来源复杂,因此,VOCs排放量及其地区和行业的分布是科学研究和决策支持都亟须的基本信息。2014年8月,生态环境部发布了《大气挥发性有机物源排放清单编制技术指南(试行)》[5](以下简称《排放清单指南》)。《排放清单指南》重点提出了VOCs排放源的分类方法、源排放量的计算方法、排放因子的取值等主要内容。

    有机溶剂是一大类在生活和生产中广泛应用的有机化合物,它存在于涂料、黏合剂、油墨和清洁剂中,其主要的排放过程包括表面涂层、表面喷漆、印刷印染、清洗、农药喷洒等,主要成分包括脂肪烃、芳香烃、氢化烃、萜烯烃、卤代烃、醇、醛酸酯、乙二醇及其衍生物、酮、醛、醚缩醛、含氮有机物及含硫有机物等。工业有机溶剂由于具有挥发性,会产生大量VOCs,是我国城市VOCs的重要来源,大约占全国人为源VOCs排放总量的20%以上[6]。工业有机溶剂由于在使用中不参与化学反应,因此,废弃后保持了原有成分的物理性质。废有机溶剂自2008年被列入《国家危险废物名录》[7]后,按照危险废物的管理要求执行,越来越受到环境管理部门的重视。

    虽然有许多学者已经开展了关于有机溶剂释放VOCs的特征研究,但是目前已有研究重点为有机溶剂使用过程中排放VOCs的特征。郑俊[8]在对上海市典型工业区大气中VOCs及其来源分析中提出,有机溶剂的使用是城市中VOCs的主要来源之一。何家禧等[9]对深圳市工业用混合溶剂进行采样分析后发现,常用工业有机溶剂多含有可挥发成分,使用范围较广,且危害度较高。莫梓伟等[10]开展了北京市使用典型溶剂企业VOCs排放成分特征。田亮等[11]研究了典型有机溶剂使用行业VOCs成分谱及臭氧生成潜势。赵锐等[12]研究了成都市工业源重点VOCs排放行业排放清单及空间分布特征。陈浩等[13]对深圳市龙岗区工作场所有机溶剂成分进行了监测分析。王肖丽等[14]基于实地调研,提出了广东省工业VOC排放清单改进建议。同时,也有针对VOCs的排放系数的研究。余翔翔等[15]进行了温州市木质家具行业VOCs排放特征及排放系数研究。梁小明等[16]对基于原料类型及末端治理的典型溶剂使用源VOCs排放系数进行了研究。李建伟等[17]进行了VOCs无组织排放估算方法和控制标准初探研究。

    但是,有机溶剂在废弃成为废有机溶剂后,仍然保持了强烈的挥发性。废有机溶剂主要的处置方式包括贮存、焚烧、再生利用、安全填埋以及水泥窑协同处置等方式[18]。以重庆市为例,每年有35.6%的废有机溶剂的处置方式为再生利用[19]。《排放清单指南》中仅列出了溶剂使用源排放因子,以计算溶剂使用过程中的排放源强,但是废有机溶剂在再生利用企业VOCs的排放情况鲜有报道[20-22],且当前研究的排放系数是由大量样本数量而得的经验值。因此,为了完善有机溶剂整个生命周期中VOCs的排放研究,本研究通过对废有机溶剂再生利用企业的调查分析,结合《排放清单指南》以及文献中提出的溶剂使用企业常见VOCs的成分,选取18种典型有机溶剂成分作为研究对象,主要包括正己烷、苯、甲苯、对二甲苯、间二甲苯、邻二甲苯、乙苯、苯乙烯、丙醛、甲醇、异丙醇、丙酮、乙酸甲酯、乙酸乙酯、二氯甲烷、均二氯乙烷、三氯乙烯、四氯化碳。本研究建立VOCs排放估算模型,以计算各情景废有机溶剂再生利用企业VOCs排放系数的差异,最终提出在不同情景模式下减少VOCs排放的合理化建议。

  • 目前,废有机溶剂再生利用技术以精蒸馏再生为主,工艺过程与石化行业的精蒸馏过程相似。因此,废有机溶剂再生利用企业的VOCs排放源可参考《石化行业VOCs污染源排查工作指南》[23](以下简称《VOCs污染源排查指南》)。《VOCs污染源排查指南》指出,本指南源项共有12类,包括设备动静密封点泄漏,有机液体储存于调和挥发损失,有机液体装卸挥发损失,废水集输、储存、处理处置过程逸散,工艺有组织排放,冷却塔、循环水冷却系统释放,非正常工况(含开停工及维修)排放,工艺无组织排放,火炬排放,燃烧烟气排放,采样过程排放以及事故排放。根据工艺环节污染源解析的理论基础,按照资料收集、源项解析、合规性检查、统计核算(包括监测/检测)、格式上报的原则对VOCs的排放量进行核算。

    本研究通过对我国废有机溶剂再生利用企业的实地调研,按照不同的操作模式以及处置规模,我国废有机溶剂再生利用企业分为大、中、小3类。为了研究3类废有机溶剂再生利用企业的排放区别,分别对3种情景进行VOCs排放源项差异性分析。

    小型再生利用企业的操作模式:由于处理量有限,回收对象为废有机溶剂产生量较少的企业,因此,收集到的废有机溶剂常以20 L塑料桶包装。进料前,需要先将20 L桶装废有机溶剂混装入200 L铁桶;在进料过程中,进料管伸入桶底以减少液面搅动而引起VOCs的挥发,由隔膜泵抽入蒸馏釜;经过精蒸馏后,产品再次通过软管伸入桶底,装料进入200 L铁桶中,具体流程及VOCs排放环节见图1

    小型再生利用企业的VOCs排放节点有5个。1) 20 L塑料桶原料混装倒入200 L铁桶的过程中,桶内废气被装入废液置换时所产生的VOCs排放溢出桶外,同时混装的过程跑冒滴漏、喷溅等意外造成废液的蒸发而产生的VOCs。2)进料过程中,进料软管伸入桶底,铁桶中废液液面下降,桶内气体空间压力降低,吸入桶外空气使桶内气体空间的蒸汽浓度下降,加速废液的蒸发。进料结束后桶内VOCs通过桶口溢出桶外,造成VOCs的排放。3)进料结束后,进料管壁外侧残留废液蒸发造成VOCs的排放。4)装料过程中,装料软管伸入桶底,桶中液面上升,桶内废气被装入的产品置换溢出桶外时所产生的VOCs排放。5)装料结束后,装料管内、外壁残留废液蒸发造成的VOCs的排放。

    中型再生利用企业的操作模式:废有机溶剂的包装为200 L铁桶,因此,不再有混装过程,其他流程同小型再生利用企业相同。具体流程及VOCs排放节点见图2

    中型再生利用企业的VOCs排放节点有4个。1)进料过程中,进料软管伸入桶底,铁桶中废液液面下降,桶内气体空间压力降低,吸入桶外空气使桶内气体空间的蒸汽浓度下降,加速废液的蒸发。进料结束后桶内VOCs通过桶口溢出桶外,造成VOCs的排放。2)进料结束后,进料管壁外侧残留废液蒸发造成VOCs的排放。3)装料过程中,装料软管伸入桶底,桶中液面上升,桶内废气被装入的产品置换时所产生的VOCs排放溢出桶外。4)装料结束后,装料管内、外壁残留废液蒸发造成VOCs的排放。

    大型再生利用企业的操作模式:大型再生利用企业回收废有机溶剂量大,常以槽罐车的形式进行包装及运输,企业内常以大型(如50 m3)原料储罐进行原料的储存,槽罐车进厂后卸料装入原料储罐中。进料时,进料管道与蒸馏釜封闭连接,无排放。装料过程使用自动装料机采用喷溅式装料。具体流程及VOCs排放节点见图3

    大型再生利用企业的VOCs排放节点有3个。1)槽罐车卸料进入50 m3原料罐时卸料造成的储罐大呼吸的蒸发损耗。2)装料过程中,装料机的装料管伸入桶口,采用喷溅式装料的方式,桶内蒸汽被装载的产品置换时产生的VOCs的排放。3)装料结束后,装料管口壁残留废液的蒸发。

    对3种再生利用企业VOCs排放节点进行分析,3种类型企业排放点的对比结果见表1

  • 《VOCs污染源排查指南》中提出企业可根据自身的情况选择核算方法[24]。排放量的核算方法通常包括实测法、公式法以及排放系数法。1)实测法。对于储罐或进料装料过程设置有机气体控制设施并有监测资料或可以进行实际检测时,推荐使用实测法进行估算。2)公式法。对于未设置有机气体控制设施的卸料、进料装料过程,按照核算方法的优先顺数,推荐使用美国环保署(EPA)发布的“污染物排放因子文件”(AP-42)(第5版)提供的评价公式,以我国有机液体理化参数和设备构造特点为基准进行估算。3)排放系数法。在企业收集的相关资料不全,无法采用实测法或公式法进行核算时,可使用系数法进行估算。为了建立全国范围的废有机溶剂再生利用企业的VOCs排放估算模型,根据实际情况,并不能对全国废有机溶剂再生利用企业进行实测,因此,公式法以及排放系数法是进行废有机溶剂再生利用企业VOCs排放估算的最佳选择方法。现对本研究提出假设,液体为理想状态下的纯液体,空气压为标准大气压,液体表面温度为25 ℃。污染源强计算式(1)~式(5)来源于《VOCs污染源排查指南》。

    工作排放的VOCs排放量LW与装料或卸料时所储蒸汽的排放有关,固定顶罐的工作排放计算方法见式(1)。

    式中:LW为工作排放的VOCs排放量;R为理想气体状态常数;TLA为日平均液体表面温度;MV为气相分子质量;PVA为真实蒸汽压;Q为周转量;KP为工作损耗产品因子,有机液体的KP=1;KN为工作排放周转因子,本研究KN=1;KB为呼吸阀工作校正因子,本研究KB=1。

    有机化学品(如苯、对二甲苯)的平均液体表面温度下的蒸汽压采用安托因方程计算,计算方法见式(2)。

    式中:PVA为平均液体表面温度下的蒸汽压;A、B、C为安托因常数;TLA为日平均液体表面温度,℃。

    在装卸排放时,计算方法见式(3)。

    式中:E装卸为装载过程的VOCs排放量,t;LL为装载损耗排放因子,kg·m−3V为物料周转量,m3η为总控制效率,%;当装卸系统未设蒸汽平衡/处理系统时,则总控制效率η取0。

    装载损耗排放因子的计算方法见式(4)。

    式中:S为饱和因子,代表排出的挥发物料接近饱和的程度,饱和因子的选取见表2C0为装载罐气、液相处于平衡状态,将挥发物料看做理想气体下的物料密度,kg·m−3C0计算方法见式(5)。

    式中:M为蒸汽的分子质量,g·mol−1PT为温度T时装载物料的真实蒸汽压,Pa;1.2×10−4为单位转换系数。

    进料/装料管残液排放量计算方法见式(6)。

    式中:E为进料管残液的排放量,kg;0.005为最低液面高度,m;S为沾染残液管壁面积,m2ρ为液体密度,kg·m−3

    在进行跑冒滴漏排放计算时,混装过程跑冒滴漏或喷溅的量采用系数法进行估算,可参考《散装液态石油产品损耗》(GB 11085-1989)中提出的灌桶损耗标准0.01%进行计算。

    排放系数的估算包括小型、中型、大型3种废有机溶剂再生利用企业排放系数的估算。

    1)小型再生利用企业排放系数。小型再生利用企业VOCs的排放包括混装、进料、进料管残留、装料以及装料管残料。其中混装排放包括喷溅式卸料和跑冒滴漏排放;进料时为工作排放量;装料时采用液下装料;进料管只有外壁有残留,排放量为E管1;而装料管与进料管尺寸相同,但内外壁均有残留,则装料管残留量为进料管2倍。因此,小型废有机溶剂再生利用企业VOCs的排放模型见式(7)。

    2)中型再生利用企业排放系数。中型再生利用企业VOCs的排放包括进料、进料管残留、装料以及装料管残料。其中进料时为工作排放量;进料管只有外壁有残留,排放量为E管1;而装料管内外壁均有残留,则装料管残留量为进料管2倍,排放系数为2E管1。因此小型再生利用企业VOCs的排放模型见式(8)。

    3)大型再生利用企业排放系数。大型再生利用企业VOCs的排放包括卸料时的工作损耗,装料机喷溅式装料以及装料机管残料。因此大型企业VOCs的排放模型见式(9)。

  • 本研究对大型、中型、小型3种废有机溶剂再生利用企业的不同VOCs排放情景模式的排放系数进行估算,结果见图4

    图4显示,3种再生利用企业中VOCs的排放系数差距较大。其中,小型再生利用企业的VOCs排放系数较高,为9.04~12.59 g·kg−1,中型再生利用企业的VOCs排放系数为8.89~10.69 g·kg−1,大型再生利用企业的VOCs排放系数为0.09~3.04 g·kg−1。比较了3种再生利用企业排放系数的比值,小型再生利用企业的排放系数为大型再生利用企业的3.64~99.93倍,中型再生利用企业的排放系数为大型再生利用企业的3.02~98.24倍。由此可见,与大型再生利用企业比,小型再生利用企业和大型再生利用企业VOCs的排放系数要高很多。其中二氯甲烷、丙醛、正己烷等VOCs的排放系数在3种类型企业中均较高,这是因为这些成分的蒸汽压较高,易挥发造成的,如表3所示。

    《排放清单指南》中列出了挥发性有机物各类源排放系数的推荐值,来自于工艺过程中有机物的制造源推荐值见表4。基础化学原料制造行业通常采用精蒸馏塔设备生产,与废有机溶剂再生利用企业的生产过程相似,而《排放清单指南》中不包括废有机溶剂再生利用过程VOCs的排放源。由表4可知,有机物的排放系数为0.097~139.74 g·kg−1,根据物质种类不同差距较大。本研究估算废有机溶剂再生利用企业的排放系数为0.090~13.35 g·kg−1,与表4中排放系数差异不大。因此,VOCs的排放清单中补充废有机溶剂再生利用过程是十分有必要的。

  • 小型再生利用企业的排放系数最高可达到大型再生利用企业的100倍左右。大、中、小型再生利用企业的排放节点不同,使得各类型企业的排放系数差距较大,图5分析了不同排放节点VOCs的排放系数。

    图5(a)所示,小型再生利用企业由于废有机溶剂的种类不同,在不同排放节点的VOCs排放系数差距较大。各排放环节的排放系数排序装料管残留>进料管残留>混装>进料>装料。在装料管残留环节的排放比例为44.09%~65.13%,在进料管残留环节的排放比例为22.04%~32.56%,在混装环节的排放比例为1.7%~17.03%,在进料环节的排放比例为0.41%~11.22%,在装料环节的排放比例为0.2%~5.61%。在选取的研究对象中,虽然排放系数的排序相同,但由于各物质蒸汽压的区别,如表4所示,在各环节的排放比例差距较大。蒸汽压越高的物质(如二氯甲烷,丙酮,正己烷等)在进料管残留以及装料管残留环节的排放比例越小,在混装、进料以及装料的排放比例越大。小型废有机溶剂再生利用企业VOCs的排放特点是,在装料管残留以及进料管残留环节的排放比例较高,而蒸汽压越低的物质在这2个环节的排放比例越高,因此小型废有机溶剂再生企业应重点关注进料管残留以及装料管的液体残留量。

    图5(b)所示,中型再生利用企业在装料管残留与出料管残留的环节排放系数较高,在装料管残留环节的排放比例最高,为53.14%~66.25%,其次为在进料管残留环节的排放,排放比例为26.57%~33.13%。在进料以及装料环节排放料较少,在进料环节的排放比例为0.42%~13.53%,在装料环节的排放比例为0.21%~6.77%。与小型企业相似,蒸汽压越低的物质在进料管残留以及装料管残留的排放比例越高。

    图5(c)所示,大型再生利用企业的主要排放环节为装料以及卸料环节。与中、小型企业的排放特别差别较大,虽然各物质的蒸汽压各不相同,但装料环节的排放比例为59.17%~59.19%,卸料环节的排放比例为40.79%~40.81%,装料管残留比例较小。可见大型再生利用企业各物质的排放比例受到蒸汽压的影响较小。

  • 图6所示,废有机溶剂再生利用企业VOCs的排放系数受到蒸汽压的影响较大,而液面温度是蒸汽压变化的重要原因。我国领土辽阔,南北跨越的经纬度近50°。冬季,我国黑龙江省漠河市的气温在−30 ℃以下,夏季新疆吐鲁番盆地平均气温高达32 ℃,各地温差较大。为了探究温度对VOCs排放的影响,本研究计算了温度为15、25和35 ℃时VOCs的排放系数。

    图6显示了大、中、小型再生利用企业温度改变对VOCs排放的影响。在温度由15 ℃升高到25 ℃的过程中,小型再生利用企业的排放系数增加了0.56%~12.04%,中型再生利用企业的排放系数增加了0.29%~7.22%,大型再生利用企业的排放系数增加了50.48%~118.33%;而由25 ℃增加到35 ℃的过程中,小型再生利用企业的排放系数增加了0.96%~15.49%,中型再生利用企业的排放系数增加了0.50%~9.49%,大型再生利用企业的排放系数增加了48.66%~102.55%。可见,大型VOCs的排放系数受到温度的影响较大。这是由于进料管残留与装料管残留排放是中、小型废有机溶剂再生利用企业重要的排放环节,而进料管残留于装料管残留VOCs的排放量与密度相关,而液体密度由于温度变化的变化量较小,可忽略不计,因此中、小型废有机溶剂再生利用企业VOCs的排放受到温度的影响较小。而大型废有机溶剂再生利用企业主要的VOCs排放来源于卸料与装料环节,这两个环节都受到温度的影响。因此大型废有机溶剂再生利用企业VOCs的排放更容易受到温度的影响。

  • 1)排放模型的不确定性分析。本研究采用了《石化行业VOCs污染源排查工作指南》[23]中质量蒸发估算计算蒸发速率,计算过程采用的是标准状态下的理想纯液体,由于废有机溶剂常常为混合物,因而与实际的平均排放因子存在一定的误差,但可以从整体上反应不同规模废有机溶剂再生利用企业VOCs的排放情况。

    2)参数的不确定性分析。本研究的参数通过调研的情况归纳总结,各企业的生产情况与生产设备,工艺环节存在一定的差异,各企业各地区的气象状况差异较大,工人操作的差异也较大,但模型可以根据实际情况,对参数进行修正。

    综上所述,由于行业活动水平存在一定的差异,无法避免地导致排放系数的估算存在一定的不确定性。但是,将废有机溶剂再生利用企业的VOCs排放系数加入VOCs排放清单是一项重要的工作。总体而言,通过本文模型的建立以及计算,可以大体反应不同规模废有机溶剂再生利用企业VOCs排放的情况。

  • 中、小型废有机溶剂再生利用企业的主要排放环节为进料管残留与装料管残留液体造成VOCs的排放。因此,建议在进料口与装料口加装废气吸收装置,同时尽量采用密闭管道进行液体的输送,以减少废气的排放。

    大型废有机溶剂再生利用企业的主要排放环节为卸料与装料。且VOCs的排放受到温度的影响较大,因此建议在卸料与装料的过程中尽量在避光通风处进行,温度较高时,可采用凉水喷淋的方式降低储罐的温度以减少VOCs的蒸发。

  • 1)大、中、小型3种再生利用企业中,排放系数的顺序为小型企业 > 中型企业 > 大型企业。中、小型再生利用企业排放系数最高的节点为装料管残留液的排放过程,大型再生利用企业排放系数最高的节点为装料的过程。物质的饱和蒸汽压越大,VOCs的排放系数也越大。

    2)温度升高,废有机溶剂中VOCs的排放系数增大,大型废有机溶剂再生利用企业受到温度的影响较大,中、小型废有机溶剂再生利用企业受到温度的影响相差不大。

    3)从排放节点来看,中、小型废有机溶剂再生利用企业因其工艺特点,工作环境,导致VOCs排放量较大,建议逐步取缔,如果作为过渡措施,可以要求企业在进料管以及装料管位置加装VOCs收集,且尽量采用密闭管道输送物料。对于大型废有机溶剂再生利用企业,应要求作业环境保持避光通风或采用喷淋等降温处理,同时在装料管位置加装VOCs收集系统。

参考文献 (25)

返回顶部

目录

/

返回文章
返回