附2
石油化工行业VOCs排放量计算办法
本办法所涉及监测和检测方法应符合相关标准规范要求。
石化行业的VOCs排放源分为:设备动静密封点泄漏;有机液体储存与调和挥发损失;有机液体装卸挥发损失;废水集输、储存、处理处臵过程逸散;燃烧烟气排放;工艺有组织排放;工艺无组织排放;采样过程排放;火炬排放;非正常工况(含开停工及维修)排放;冷却塔、循环水冷却系统释放;事故排放等12类源项。
企业某个核算周期(以年计)VOCs排放量为:
11E石化=?Em
m?1(公式1)
式中: Em 年。
各源项污染源的VOCs排放量应为该源项每一种污染物排放量的加和,见公式 2:
n石化行业各源项污染源VOCs排放量,千克/
Em??Ei
i?1(公式 2)
1
式中: Ei
式中:
某源项污染源排放的污染物i的排放量,千克/年
N?WFi?(公式3) Ei=??E排放源n,i??
WFVOCs?n?1?污染物i的排放量,千克/年;
Ei
E排放源n,i含污染物i的第n个排放源的VOCs排放量, N WFi
千克/年;
含污染物i的排放源总数;
流经或储存于排放源的物料中污染物i的平均 质量分数;
WFVOCs 流经或储存于排放源的物料中VOCs的平均质
量分数。
有机液体储存与调和、装卸过程中涉及附表1中单一物质的,应按本办法进行单一污染物排放量核算,并可在VOCs总量中予以扣除。鼓励有条件的企业进行全过程单一污染物排放量核算。
进入气相的VOCs,可按以下方法进行核算:VOCs排放量=废气处理设施未投用的排放量+废气处理设施投用但未收集的排放量+废气处理设施投用收集后未去除的排放量=VOCs产生量总量-废气处理设施投用收集且去除的量。
2
一、设备动静密封点泄漏
排放量核算结果的准确度从高到低排序为:实测法、相关方程法、筛选范围法、平均排放系数法。前三种方法是基于检测的核算方法,需获得检测仪器对物料的(合成)响应因子,见附录一。企业可根据自身LDAR开展情况,选择核算方法。
E设备式中: E设备 ti eTOCs,i
??WFVOCs,i???eTOC,i??ti?
??WFTOC,ii?1??n(公式4)
密封点的VOCs年排放量,千克/年; 密封点i的运行时间段,小时/年; 密封点i的TOCs排放速率,千克/小时;
WFVOCs,i 运行时间段内流经密封点i的物料中VOCs的
平均质量分数;
WFTOC,i 运行时间段内流经密封点i的物料中TOC的平
均质量分数
WFVOCs如未提供物料中VOCs的平均质量分数,则按1
WFTOC计。
(一)排放速率。 1.实测法。
采用包袋法和大体积采样法对密封点进行实测,所得排
3
放速率最接近真实排放情况,企业可选用该方法对密封点排放速率进行检测。
2.相关方程法。
eTOC?e0,in????ep,ii?1??ef,i?0?SV?1??SV?50000? ?1?SV?50000?(公式5)
式中: eTOC SV e0,i ep,i ef,i 时。
各类型密封点的排放速率按表1计算。
表1 石油炼制和石油化工设备组件的设备排放速率a
限定排放速率(千相关方程b(千克/小时/克/小时/排放源) 密封点类型 排放源) >50000 μmol/mol 石油炼制的排放速率(炼油、营销终端和油气生产) 7.8E-06 0.14 2.29E-06×SV0.746 阀门 2.4E-05 0.16 5.03E-05×SV0.610 泵 4.0E-06 0.11 1.36E-05×SV0.589 其它 7.5E-06 0.030 1.53E-06×SV0.735 连接件 3.1E-07 0.084 4.61E-06×SV0.703 法兰 2.0E-06 0.079 2.20E-06×SV0.704 开口阀或开口管线 石油化工的排放速率 6.6E-07 0.11 1.87E-06×SV0.873 气体阀门 4.9E-07 0.15 6.41E-06×SV0.797 液体阀门 7.5E-06 0.62 1.90E-05×SV0.824 轻液体泵c 6.1E-07 0.22 3.05E-06×SV0.885 连接件 默认零值排放速率(千克/小时/排放源) 注:表中涉及的千克/小时/排放源=每个排放源每小时的TOC排放量(千克)。 密封点的TOC排放速率,千克/小时; 修正后的净检测值,μmol/mol;
密封点i的默认零值排放速率,千克/小时; 密封点i的限定排放速率,千克/小时; 密封点i的相关方程核算排放速率,千克/小
4
a:美国环保署,1995b报告的数据。对于密闭式的采样点,如果采样瓶连在采样口,则使用“连接件”的排放系数;如采样瓶未与采样口连接,则使用“开口管线”的排放系数。
b:SV是检测设备测得的净检测值(SV,μmol/mol); c:轻液体泵系数也可用于压缩机、泄压设备和重液体泵。
3.筛选范围法。
石油炼制工业排放速率计算公式:
n??WFTOC,ieTOC???FA,i??WFTOC,i?Ni? (公式6)
??WFTOC,i?WF甲烷,ii?1??石油化学工业排放速率计算公式:
eTOC?
??Fi?1nA,i?WFTOC,i?Ni?
(公式7)
式中: eTOC FA,i WFTOC WF甲烷 Ni
密封点的TOC排放速率,千克/小时; 密封点i排放系数;
流经密封点i的物料中TOC的平均质量分数; 流经密封点i的物料中甲烷的平均质量分数, 最大取10%; 密封点的个数。
表2 筛选范围排放系数a 石油炼制系数 设备类型 ≥介质 10000μmol/mol排放系数(千克/小时/排放源) <10000μmol/mol排放系数(千克/小时/排放源) b石油化工系数 ≥10000μmol/mol排放系数(千克/小时/排放源) 0.113 <10000μmol/mol排放系数(千克/小时/排放源) 0.000081 c法兰、连0.0375 0.00006 所有 接件 注:a:EPA,1995b报告的数据。 b:这些系数是针对非甲烷有机化合物排放。 c:这些系数是针对总有机化合物排放。
筛选范围法用于核算某套装臵不可达法兰或连接件的VOCs排放速率,需至少检测50%该装臵的可达法兰或连接
5
件,并且至少包含1个净检测值大于等于10000μmol/mol的点,以10000μmol/mol为界,分析已检测法兰或连接件净检测值可能≥10000μmol/mol的数量比例,将该比例应用到同一装臵的不可达法兰或连接件,且按比例计算的大于等于10000μmol/mol的不可达点个数向上取整,采用表2系数并按公式6和公式7计算排放速率。
4.平均排放系数法。
未开展LDAR工作的企业,或不可达点(除符合筛选范围法适用范围的法兰和连接件外),应采用表3系数并按公式6和公式7计算排放速率。
表3 石油炼制和石油化工组件平均排放系数a
设备类型 介质 石油炼制排放系数 (千克/小时/排放源)b 0.0268 0.0109 0.00023 0.114 0.021 0.636 0.16 0.00025 0.0023 0.0150 石油化工排放系数 (千克/小时/排放源)c 0.00597 0.00403 0.00023 0.0199 0.00862 0.228 0.104 0.00183 0.0017 0.0150 气体 阀 轻液体 重液体 轻液体 泵d 重液体 压缩机 气体 泄压设备 气体 法兰、连接件 所有 开口阀或开口所有 管线 采样连接系统 所有 注:对于表中涉及的千克/小时/排放源=每个排放源每小时的TOC排放量(千克)。对于开放式的采样点,采用平均排放系数法计算排放量。如果采样过程中排出的臵换残液或气未经处理直接排入环境,按照“取样连接系统”和“开口管线”排放系数分别计算并加和;如果企业有收集处理设施收集管线冲洗的残液或气体,并且运行效果良好,可按“开口阀或开口管线”排放系数进行计算。
a:摘自EPA,1995b;
b:石油炼制排放系数用于非甲烷有机化合物排放速率; c:石油化工排放系数用于TOC(包括甲烷)排放速率;
d:轻液体泵密封的系数可以用于估算搅拌器密封的排放速率。
6
(二)排放时间。
采用中点法确定该密封点的排放时间,即第n次检测值代表时间段的起始点为第n-1次至第n次检测时间段的中点,终止点为第n次至第n+1次检测时间段的中点。发生泄漏修复的情况下,修复复测的时间点为泄漏时间段的终止点。
二、有机液体储存与调和挥发损失 (一)实测法。
设有VOCs有机气体控制设施的储罐或罐区,其排放量应采用实测法核算,监测频次不少于1次/月。
?6E储罐???E?10?(C入口,i?C出口,i)?Qi?ti??计算量,i?
i?1n (公式8)
式中: E储罐
含有机气体控制设施的储罐VOCs年排放量 千克/年;
E计算量,i 连接在有机气体控制设施i上的储罐的排放
C入口,i C出口,i
量,由公式法计算,千克/年;
有机气体控制设施i的入口VOCs浓度年平均值,毫克/标立方米;
有机气体控制设施i的出口VOCs浓度年平均值,毫克/标立方米;
7
Qi
ti
有机气体控制设施i的出口流量,标立方米/小时;
有机气体控制设施i的运行时间,小时/年。
(二)公式法。
该核算方法可应用于固定顶罐和浮顶罐。不适用于以下情况:所储物料组分不稳定或真实蒸汽压高于大气压、蒸气压未知或无法测量的;储罐浮盘设施失效的;其他不符合相关环保要求的。
公式法核算过程采用美制单位。完成核算后,可将排放量的美制单位(磅)转为国际单位制(千克)。
nE储罐???E固,i?E浮,i?
i?1(公式 9)
式中: E储罐 E固 E浮
储罐的VOCs年排放量,千克/年; 固定顶罐i的VOCs年排放量,千克/年; 浮顶罐i的VOCs年排放量,千克/年。
(公式10)
1.固定顶罐总损失。
E固?ES?EW 式中: E固 ES EW
固定顶罐总损失,磅/年; 静臵损失,磅/年,见公式11; 工作损失,磅/年,见公式28。
8
(1)静臵损失,ES。
ES?365VVWVKEKS 式中: ES VV WV KE KS
(公式11)
静臵损失(地下卧式罐的ES取0),磅/年; 气相空间容积,立方英尺,见公式12; 储藏气相密度,磅/立方英尺; 气相空间膨胀因子,无量纲量; 排放蒸气饱和因子,无量纲量。
立式罐气相空间容积VV,通过公式12计算:
??2?(公式12) VV??D?HVO
?4?式中: VV D HVO
气相空间容积,立方英尺; 罐径,英尺;
气相空间高度,英尺。
卧式罐气相空间容积VV,通过公式13核算: ?2VV?DEHVO (公式13)
4式中: VV HVO DE
固定顶罐蒸气空间体积,立方英尺; 蒸气实际空间高度(HVO=D),英尺; 卧式罐有效直径,英尺;
LD DE?0.7859
(公式14)
A.气相空间膨胀因子KE KE?0.0018?TV?0.0018??0.72?TAX?TAN??0.028?I??
(公式15)
式中: KE ΔTV TAX TAN α I 0.0018
0.72 0.028
气相空间膨胀因子,无量纲量; 日蒸气温度范围,兰氏度; 日最高环境温度,兰氏度; 日最低环境温度,兰氏度;
罐漆太阳能吸收率,无量纲量,见表4; 太阳辐射强度,英热/(平方英尺〃天);常数,(兰氏度)-1; 常数,无量纲量;
常数,兰氏度〃平方英尺〃天 / 英热。
表4 罐漆太阳能吸收率(α)
序号 1 2 3 罐漆颜色 白色 铝色 黑色 太阳能吸收因子 0.34 0.68 0.97 序号 4 5 6 罐漆颜色 浅灰色 中灰色 绿色 太阳能吸收因子 0.63 0.74 0.91
B.气相空间高度,HVO
HVO?HS?HL?HRO 式中: HVO HS HL
气相空间高度,英尺; 罐体高度,英尺; 液体高度,英尺;
10
(公式16)
HRO
罐顶计量高度,英尺,锥顶罐见注释a,拱顶 罐见注释b。
公式16注释:
a.对于锥顶罐,顶高度HRO核算方法如下:
HRO?1/3HR (公式17) 式中: HRO HR
罐顶计量高度,英尺; 罐顶高度,英尺; HR?SRRS
(公式18)
式中: SR 0.0625;
RS
罐壳半径,英尺。
罐锥顶斜率,英尺/英尺;如未知,则取
b.对于拱顶罐,罐顶计量高度HRO核算方法如下:
?11?H?2?HRO?HR???R?? (公式19)
26?RS?????式中: HRO RS HR
罐顶计量高度,英尺; 罐壳半径,英尺; 罐顶高度,英尺;
HR?RR??R?R2R20.5S
? (公式20)
RR RS
罐拱顶半径,英尺; 罐壳半径,英尺;
11
RR的值一般介于0.8D-1.2D之间,其中D=2RS。如果RR未知,则用罐体直径代替。
C.气相空间饱和因子,KS
1 Ks?1?0.053PVAHVO式中: KS PVA
气相空间饱和因子,无量纲量;
日平均液面温度下的饱和蒸气压,磅/平方英
(公式21)
寸(绝压),或参照公式26计算;
HVO 0.053
气相空间高度,英尺,见公式16; 常数,(磅/平方英寸(绝压)〃英尺)-1。
D.气相密度,WV
MVPVA WV?RTLA式中: WV MV R PVA
气相密度,磅/立方英尺; 气相分子质量,磅/磅-摩尔;
(公式22)
理想气体状态常数,10.741磅/(磅-摩尔〃英 尺〃兰氏度);
日平均液面温度下的饱和蒸气压,磅/平方英
寸(绝压),见公式26;
TLA
日平均液体表面温度,兰氏度,取年平均实
际储存温度,如无该数据,用公式23计算。
12
公式22注释:
a.日平均液体表面温度,TLA
TLA?0.44TAA?0.56TB?0.0079?I
?T?TAN?TAA??AX?
2??TB?TAA?6??1 式中: TLA TAA TAX TAN TB α I
日平均液体表面温度,兰氏度; 日平均环境温度,兰氏度;
(公式23) (公式24) (公式25)
计算月的日最高环境温度,兰氏度; 计算月的日最低环境温度,兰氏度。 储液主体温度,兰氏度;
罐漆太阳能吸收率,无量纲量,见表4; 太阳辐射强度,英热/(平方英尺〃天)。
表5 年平均储藏温度计算表 罐体颜色 白 铝 灰 黑 年平均储藏温度,TS(华氏度) TAA+0 TAA+2.5 TAA+3.5 TAA+5.0 当TLA值无法取得时,可用表5计算。
注:此表格中TAA为年平均环境温度(华氏度)。 E.真实蒸气压,PVA
对于石油液体储料的日平均液体表面蒸气压,可通过公式26计算:
??B??P?? VA?exp?A???TLA???13
(公式26)
式中: A B TLA PVA 压)。
对于油品:
蒸气压公式中的常数,无量纲量; 蒸气压公式中的常数,兰氏度; 日平均液体表面温度,兰氏度;
日平均液体表面蒸气压,磅/平方英寸(绝
A?15.64-1.854S0.5-?0.8742-0.3280S0.5?ln?RVP?B?8742-1042S0.5-?1049-179.4S0.5?ln?RVP?对于原油:
A?12.82-0.9672ln?RVP?
B?7261?1216ln?RVP?
式中: RVP S /vol%。
雷德蒸气压,磅/平方英寸;
10%蒸发量下ASTM蒸馏曲线斜率,°F
15%馏出温度-5%馏出温度S?
15?5单一物质(如苯、对二甲苯)的日平均液体表面蒸气压,
采用安托因方程计算。
?B?lgP?VA?A???TLA?C?
14
(公式27)
式中:
A、B、C 为安托因常数; TLA PVA
日平均液体表面温度,摄氏度;
日平均液面温度下的饱和蒸气压,毫米汞柱。
(2)工作损失,EW。
5.614EW?MVPVAQKNKPKB
RTLA式中: EW MV TLA PVA
工作损失,磅/年; 气相分子量,磅/磅-摩尔; 日平均液体表面温度,兰氏度;
(公式28)
真实蒸气压,磅/平方英寸(绝压),见公式
26;
Q KP KN KB
年周转量,桶/年;
工作损失产品因子,无量纲量;
对于原油KP=0.75; 对于其它有机液体KP=1;
工作排放周转(饱和)因子,无量纲量;
当周转数>36,KN=(180+N)/6N; 当周转数≤36,KN=1;
呼吸阀工作校正因子。
呼吸阀工作时的校正因子,KB可用公式29和公式30计
15
算:
当
?PBP?PA?KN???1.0
?PI?PA?(公式29)
时
?PI?PA??PVA??KN? KB???PBP?PA?PVA????? (公式30)
式中: KB PI PI PA KN PVA PBP
呼吸阀校正因子,无量纲量;
正常工况条件下气相空间压力,磅/平方英寸 (表压);
是一个实际压力(表压),如果处在大气压下 (不是真空或处在稳定压力下),PI为0; 大气压,磅/平方英寸(绝压);
工作排放周转(饱和)因子,无量纲量,见 公式28;
日平均液面温度下的蒸气压,磅/平方英寸(绝 压),见公式26;
呼吸阀压力设定,磅/平方英寸(表压)。
(公式31)
2.浮顶罐总损失。
E浮?ER?EWD?EF?ED 式中:
16
E浮 ER EWD EF ED
浮顶罐总损失,磅/年;
边缘密封损失,磅/年,见公式32; 挂壁损失,磅/年,见公式34; 浮盘附件损失,磅/年,见公式35;
浮盘缝隙损失(只限螺栓连接式的浮盘或浮顶),磅/年,见公式38。
(1)边缘密封损失,ER。
ER??KRa?KRbvn?DP*MVKC
(公式32)
式中: ER KRa
边缘密封损失,磅/年;
零风速边缘密封损失因子,磅-摩尔/英
尺〃年,见表6; KR
v n P*
n
有风时边缘密封损失因子, 磅-摩尔(迈/〃英
尺〃年),见表6; 罐点平均环境风速,迈;
密封相关风速指数,无量纲量,见表6; 蒸气压函数,无量纲量;
PVAPA* P?2??P?0.5??1??1?VA??PA??????(公式33)
式中: PVA
日平均液体表面蒸气压,磅/平方英寸(绝
17
压),见公式26;
PA D MV KC
1。
表6 浮顶罐边缘密封损失系数 罐体类型 密封 KRa KRb (磅-摩尔/英尺〃年) 磅-摩尔/(迈n〃英尺〃年) 机械式鞋形密封 5.8 1.6 0.6 液体镶嵌式(接触液面) 焊接 只有一级 挡雨板 边缘刮板 只有一级 挡雨板 边缘刮板 只有一级 边缘靴板 边缘刮板 1.6 0.7 0.3 气体镶嵌式(不接触液面) 6.7 3.3 2.2 机械式鞋形密封 铆接 10.8 9.2 1.1 0.4 0.2 0.3 2.0 1.9 1.5 0.2 0.1 0.003 3.0 3.0 4.3 0.3 0.3 0.6 1.5 1.2 0.3 0.3 0.3 0.4 n 2.1 1.6 1.0 大气压,磅/平方英寸(绝压); 罐体直径,英尺;
气相分子质量,磅/磅-摩尔;
产品因子,原油0.4,其它挥发性有机液体
只有一级 边缘靴板 边缘刮板 注:表中边缘密封损失因子kra,krb, n只适用于风速6.8米/秒以下
(2)挂壁损失,EWD。
0.943?QCSWL?EWD?D式中: EWD Q
?NCFC??1?D? ??(公式34)
挂壁损失,磅/年; 年周转量,桶/年;
18
相关推荐: