附2 石化行业VOCs排放量计算办法模板.docx

附2 石化行业VOCs排放量计算办法模板.docx

《附2 石化行业VOCs排放量计算办法模板.docx》由会员分享，可在线阅读，更多相关《附2 石化行业VOCs排放量计算办法模板.docx（58页珍藏版）》请在冰豆网上搜索。

附2石化行业VOCs排放量计算办法模板

附2

石油化工行业VOCs排放量计算办法

本办法所涉及监测和检测方法应符合相关标准规范要求。

石化行业的VOCs排放源分为：

设备动静密封点泄漏；有机液体储存与调和挥发损失；有机液体装卸挥发损失；废水集输、储存、处理处置过程逸散；燃烧烟气排放；工艺有组织排放；工艺无组织排放；采样过程排放；火炬排放；非正常工况（含开停工及维修）排放；冷却塔、循环水冷却系统释放；事故排放等12类源项。

企业某个核算周期（以年计）VOCs排放量为：

（

公式1）

式中：

Em石化行业各源项污染源VOCs排放量，千克/年。

各源项污染源的VOCs排放量应为该源项每一种污染物排放量的加和，见公式2：

（

公式2）

式中：

Ei某源项污染源排放的污染物i的排放量，千克/年

（

公式3）

式中：

Ei污染物i的排放量，千克/年；

E排放源n,i含污染物i的第n个排放源的VOCs排放量，

千克/年；

N含污染物i的排放源总数；

WFi流经或储存于排放源的物料中污染物i的平均

质量分数；

WFVOCs流经或储存于排放源的物料中VOCs的平均质

量分数。

有机液体储存与调和、装卸过程中涉及附表1中单一物质的，应按本办法进行单一污染物排放量核算，并可在VOCs总量中予以扣除。

鼓励有条件的企业进行全过程单一污染物排放量核算。

进入气相的VOCs，可按以下方法进行核算：

VOCs排放量=废气处理设施未投用的排放量+废气处理设施投用但未收集的排放量+废气处理设施投用收集后未去除的排放量=VOCs产生量总量-废气处理设施投用收集且去除的量。

一、设备动静密封点泄漏

排放量核算结果的准确度从高到低排序为：

实测法、相关方程法、筛选范围法、平均排放系数法。

前三种方法是基于检测的核算方法，需获得检测仪器对物料的（合成）响应因子，见附录一。

企业可根据自身LDAR开展情况，选择核算方法。

（

公式4）

式中：

E设备密封点的VOCs年排放量，千克/年；

ti密封点i的运行时间段，小时/年；

eTOCs,i密封点i的TOCs排放速率，千克/小时；

WFVOCs,i运行时间段内流经密封点i的物料中VOCs的

平均质量分数；

WFTOC,i运行时间段内流经密封点i的物料中TOC的平

均质量分数

如未提供物料中VOCs的平均质量分数，则

按1计。

（一）排放速率。

1.实测法。

采用包袋法和大体积采样法对密封点进行实测，所得排放速率最接近真实排放情况，企业可选用该方法对密封点排放速率进行检测。

2.相关方程法。

（

公式5）

式中：

eTOC密封点的TOC排放速率，千克/小时；

SV修正后的净检测值，μmol/mol；

e0,i密封点i的默认零值排放速率，千克/小时；

ep,i密封点i的限定排放速率，千克/小时；

ef,i密封点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

轻液体泵c

7.5E-06

0.62

1.90E-05×SV0.824

连接件

6.1E-07

0.22

3.05E-06×SV0.885

注：

表中涉及的千克/小时/排放源＝每个排放源每小时的TOC排放量（千克）。

a：

美国环保署，1995b报告的数据。

对于密闭式的采样点，如果采样瓶连在采样口，则使用“连接件”的排放系数；如采样瓶未与采样口连接，则使用“开口管线”的排放系数。

b：

SV是检测设备测得的净检测值（SV，μmol/mol）；

c：

轻液体泵系数也可用于压缩机、泄压设备和重液体泵。

3.筛选范围法。

石油炼制工业排放速率计算公式：

（

公式6）

石油化学工业排放速率计算公式：

（

公式7）

式中：

eTOC密封点的TOC排放速率，千克/小时；

FA,i密封点i排放系数；

WFTOC流经密封点i的物料中TOC的平均质量分数；

WF甲烷流经密封点i的物料中甲烷的平均质量分数，

最大取10%；

Ni密封点的个数。

表2筛选范围排放系数a

设备类型

介质

石油炼制系数b

石油化工系数c

≥10000μmol/mol排放系数（千克/小时/排放源）

<10000μmol/mol排放系数（千克/小时/排放源）

≥10000μmol/mol排放系数（千克/小时/排放源）

<10000μmol/mol排放系数（千克/小时/排放源）

法兰、连接件

所有

0.0375

0.00006

0.113

0.000081

注：

a：

EPA，1995b报告的数据。

b：

这些系数是针对非甲烷有机化合物排放。

c：

这些系数是针对总有机化合物排放。

筛选范围法用于核算某套装置不可达法兰或连接件的VOCs排放速率，需至少检测50%该装置的可达法兰或连接件，并且至少包含1个净检测值大于等于10000µmol/mol的点，以10000µmol/mol为界，分析已检测法兰或连接件净检测值可能≥10000µmol/mol的数量比例，将该比例应用到同一装置的不可达法兰或连接件，且按比例计算的大于等于10000µmol/mol的不可达点个数向上取整，采用表2系数并按公式6和公式7计算排放速率。

4.平均排放系数法。

未开展LDAR工作的企业，或不可达点（除符合筛选范围法适用范围的法兰和连接件外），应采用表3系数并按公式6和公式7计算排放速率。

表3石油炼制和石油化工组件平均排放系数a

设备类型

介质

石油炼制排放系数

（千克/小时/排放源）b

石油化工排放系数

（千克/小时/排放源）c

阀

气体

0.0268

0.00597

轻液体

0.0109

0.00403

重液体

0.00023

0.00023

泵d

轻液体

0.114

0.0199

重液体

0.021

0.00862

压缩机

气体

0.636

0.228

泄压设备

气体

0.16

0.104

法兰、连接件

所有

0.00025

0.00183

开口阀或开口管线

所有

0.0023

0.0017

采样连接系统

所有

0.0150

0.0150

注：

对于表中涉及的千克/小时/排放源＝每个排放源每小时的TOC排放量（千克）。

对于开放式的采样点，采用平均排放系数法计算排放量。

如果采样过程中排出的置换残液或气未经处理直接排入环境，按照“取样连接系统”和“开口管线”排放系数分别计算并加和；如果企业有收集处理设施收集管线冲洗的残液或气体，并且运行效果良好，可按“开口阀或开口管线”排放系数进行计算。

a：

摘自EPA，1995b；

b：

石油炼制排放系数用于非甲烷有机化合物排放速率；

c：

石油化工排放系数用于TOC（包括甲烷）排放速率；

d：

轻液体泵密封的系数可以用于估算搅拌器密封的排放速率。

（二）排放时间。

采用中点法确定该密封点的排放时间，即第n次检测值代表时间段的起始点为第n-1次至第n次检测时间段的中点，终止点为第n次至第n+1次检测时间段的中点。

发生泄漏修复的情况下，修复复测的时间点为泄漏时间段的终止点。

二、有机液体储存与调和挥发损失

（一）实测法。

设有VOCs有机气体控制设施的储罐或罐区，其排放量应采用实测法核算，监测频次不少于1次/月。

（

公式8）

式中：

E储罐含有机气体控制设施的储罐VOCs年排放量

千克/年；

E计算量，i连接在有机气体控制设施i上的储罐的排放量，由公式法计算，千克/年；

C入口,i有机气体控制设施i的入口VOCs浓度年平均值，毫克/标立方米；

C出口,i有机气体控制设施i的出口VOCs浓度年平均值，毫克/标立方米；

Qi有机气体控制设施i的出口流量，标立方米/小时；

ti有机气体控制设施i的运行时间，小时/年。

（二）公式法。

该核算方法可应用于固定顶罐和浮顶罐。

不适用于以下情况：

所储物料组分不稳定或真实蒸汽压高于大气压、蒸气压未知或无法测量的；储罐浮盘设施失效的；其他不符合相关环保要求的。

公式法核算过程采用美制单位。

完成核算后，可将排放量的美制单位（磅）转为国际单位制（千克）。

（

公式9）

式中：

E储罐储罐的VOCs年排放量，千克/年；

E固固定顶罐i的VOCs年排放量，千克/年；

E浮浮顶罐i的VOCs年排放量，千克/年。

1.固定顶罐总损失。

（

公式10

）

式中：

E固固定顶罐总损失，磅/年；

ES静置损失，磅/年，见公式11；

EW工作损失，磅/年，见公式28。

（1）静置损失，ES。

（

公式11）

式中：

ES静置损失（地下卧式罐的ES取0），磅/年；

VV气相空间容积，立方英尺，见公式12；

WV储藏气相密度，磅/立方英尺；

KE气相空间膨胀因子，无量纲量；

KS排放蒸气饱和因子，无量纲量。

立式罐气相空间容积VV，通过公式12计算：

（

公式12）

式中：

VV气相空间容积，立方英尺；

D罐径，英尺；

HVO气相空间高度，英尺。

卧式罐气相空间容积VV，通过公式13核算：

（

公式13）

式中：

VV固定顶罐蒸气空间体积，立方英尺；

HVO蒸气实际空间高度（HVO=D），英尺；

DE卧式罐有效直径，英尺；

（

公式14）

A.气相空间膨胀因子KE

（

公式15）

式中：

KE气相空间膨胀因子，无量纲量；

ΔTV日蒸气温度范围，兰氏度；

TAX日最高环境温度，兰氏度；

TAN日最低环境温度，兰氏度；

α罐漆太阳能吸收率，无量纲量，见表4；

I太阳辐射强度，英热/（平方英尺·天）；0.0018常数，（兰氏度）-1；

0.72常数，无量纲量；

0.028常数，兰氏度·平方英尺·天/英热。

表4罐漆太阳能吸收率（α）

序号

罐漆颜色

太阳能吸收因子

序号

罐漆颜色

太阳能吸收因子

1

白色

0.34

4

浅灰色

0.63

2

铝色

0.68

5

中灰色

0.74

3

黑色

0.97

6

绿色

0.91

B.气相空间高度，HVO

（

公式16）

式中：

HVO气相空间高度，英尺；

HS罐体高度，英尺；

HL液体高度，英尺；

HRO罐顶计量高度，英尺，锥顶罐见注释a，拱顶

罐见注释b。

公式16注释：

a.对于锥顶罐，顶高度HRO核算方法如下：

（

公式17）

式中：

HRO罐顶计量高度，英尺；

HR罐顶高度，英尺；

（

公式18）

式中：

SR罐锥顶斜率，英尺/英尺；如未知，则取0.0625；

RS罐壳半径，英尺。

b.对于拱顶罐，罐顶计量高度HRO核算方法如下：

（

公式19）

式中：

HRO罐顶计量高度，英尺；

RS罐壳半径，英尺；

HR罐顶高度，英尺；

（

公式20）

RR罐拱顶半径，英尺；

RS罐壳半径，英尺；

RR的值一般介于0.8D-1.2D之间，其中D=2RS。

如果RR未知，则用罐体直径代替。

C.气相空间饱和因子，KS

（

公式21）

式中：

KS气相空间饱和因子，无量纲量；

PVA日平均液面温度下的饱和蒸气压，磅/平方英寸（绝压），或参照公式26计算；

HVO气相空间高度，英尺，见公式16；

0.053常数，（磅/平方英寸（绝压）·英尺）-1。

D.气相密度，WV

（

公式22）

式中：

WV气相密度，磅/立方英尺；

MV气相分子质量，磅/磅-摩尔；

R理想气体状态常数，10.741磅/（磅-摩尔·英

尺·兰氏度）；

PVA日平均液面温度下的饱和蒸气压，磅/平方英寸（绝压），见公式26；

TLA日平均液体表面温度，兰氏度，取年平均实际储存温度，如无该数据，用公式23计算。

公式22注释：

a.日平均液体表面温度，TLA

（

公式23）

（

公式24）

（

公式25）

式中：

TLA日平均液体表面温度，兰氏度；

TAA日平均环境温度，兰氏度；

TAX计算月的日最高环境温度，兰氏度；

TAN计算月的日最低环境温度，兰氏度。

TB储液主体温度，兰氏度；

α罐漆太阳能吸收率，无量纲量，见表4；

I太阳辐射强度，英热/（平方英尺·天）。

当TLA值无法取得时，可用表5计算。

表5年平均储藏温度计算表

罐体颜色

年平均储藏温度，TS（华氏度）

白

TAA+0

铝

TAA+2.5

灰

TAA+3.5

黑

TAA+5.0

注：

此表格中TAA为年平均环境温度（华氏度）。

E.真实蒸气压，PVA

对于石油液体储料的日平均液体表面蒸气压，可通过公式26计算：

（

公式26）

式中：

A蒸气压公式中的常数，无量纲量；

B蒸气压公式中的常数，兰氏度；

TLA日平均液体表面温度，兰氏度；

PVA日平均液体表面蒸气压，磅/平方英寸（绝压）。

对于油品：

对于原油：

式中：

RVP雷德蒸气压，磅/平方英寸；

S10%蒸发量下ASTM蒸馏曲线斜率，°F/vol%。

单一物质（如苯、对二甲苯）的日平均液体表面蒸气压，采用安托因方程计算。

（

公式27）

式中：

A、B、C为安托因常数；

TLA日平均液体表面温度，摄氏度；

PVA日平均液面温度下的饱和蒸气压，毫米汞柱。

（2）工作损失，EW。

（

公式28）

式中：

EW工作损失，磅/年；

MV气相分子量，磅/磅-摩尔；

TLA日平均液体表面温度，兰氏度；

PVA真实蒸气压，磅/平方英寸（绝压），见公式26；

Q年周转量，桶/年；

KP工作损失产品因子，无量纲量；

对于原油KP=0.75；

对于其它有机液体KP=1；

KN工作排放周转（饱和）因子，无量纲量；

当周转数＞36，KN=（180+N）/6N；

当周转数≤36，KN=1；

KB呼吸阀工作校正因子。

呼吸阀工作时的校正因子，KB可用公式29和公式30计算：

当

（

公式29）

时

（

公式30）

式中：

KB呼吸阀校正因子，无量纲量；

PI正常工况条件下气相空间压力，磅/平方英寸

（表压）；

PI是一个实际压力（表压），如果处在大气压下

（不是真空或处在稳定压力下），PI为0；

PA大气压，磅/平方英寸（绝压）；

KN工作排放周转（饱和）因子，无量纲量，见

公式28；

PVA日平均液面温度下的蒸气压，磅/平方英寸（绝

压），见公式26；

PBP呼吸阀压力设定，磅/平方英寸（表压）。

2.浮顶罐总损失。

（

公式31）

式中：

E浮浮顶罐总损失，磅/年；

ER边缘密封损失，磅/年，见公式32；

EWD挂壁损失，磅/年，见公式34；

EF浮盘附件损失，磅/年，见公式35；

ED浮盘缝隙损失（只限螺栓连接式的浮盘或浮顶），磅/年，见公式38。

（1）边缘密封损失，ER。

（

公式32）

式中：

ER边缘密封损失，磅/年；

KRa零风速边缘密封损失因子，磅-摩尔/英尺·年，见表6；

KR有风时边缘密封损失因子，磅-摩尔/（迈n·英尺·年），见表6；

v罐点平均环境风速，迈；

n密封相关风速指数，无量纲量，见表6；

P*蒸气压函数，无量纲量；

（

公式33）

式中：

PVA日平均液体表面蒸气压，磅/平方英寸（绝压），见公式26；

PA大气压，磅/平方英寸（绝压）；

D罐体直径，英尺；

MV气相分子质量，磅/磅-摩尔；

KC产品因子，原油0.4，其它挥发性有机液体1。

表6浮顶罐边缘密封损失系数

罐体类型

密封

KRa

（磅-摩尔/英尺·年）

KRb

磅-摩尔/（迈n·英尺·年）

n

焊接

机械式鞋形密封

只有一级

5.8

0.3

2.1

边缘靴板

1.6

0.3

1.6

边缘刮板

0.6

0.4

1.0

液体镶嵌式（接触液面）

只有一级

1.6

0.3

1.5

挡雨板

0.7

0.3

1.2

边缘刮板

0.3

0.6

0.3

气体镶嵌式（不接触液面）

只有一级

6.7

0.2

3.0

挡雨板

3.3

0.1

3.0

边缘刮板

2.2

0.003

4.3

铆接

机械式鞋形密封

只有一级

10.8

0.4

2.0

边缘靴板

9.2

0.2

1.9

边缘刮板

1.1

0.3

1.5

注：

表中边缘密封损失因子kra，krb，n只适用于风速6.8米/秒以下

（2）挂壁损失，EWD。

（

公式34）

式中：

EWD挂壁损失，磅/年；

Q年周转量，桶/年；

CS罐体油垢因子，见表7；

WL有机液体密度，磅/加仑，部分物料参数见附表2；

D罐体直径，英尺；

0.943常数，1000立方英尺·加仑/桶2；

NC固定顶支撑柱数量（对于自支撑固定浮顶或外浮顶罐：

NC=0。

），无量纲量；

FC有效柱直径，取值1。

表7储罐罐壁油垢因子

介质

罐壁状况（桶/1000平方英尺）

轻锈

中锈

重锈

汽油

0.0015

0.0075

0.15

原油

0.006

0.03

0.6

其它油品

0.0015

0.0075

0.15

备注：

储罐内壁平均3年以上（包括3年）除锈一次，为重锈；平均两年除锈一次，为中锈；平均每年除锈一次，为轻绣。

（3）浮盘附件损失，EF。

（

公式35）

式中：

EF浮盘附件损失，磅/年；

FF总浮盘附件损失因子，磅-摩尔/年；

（

公式36）

式中：

NFi某类浮盘附件数，无量纲量；

KFi某类附件损失因子，磅-摩尔/年，见公式37；

nf某类的附件总数，无量纲量；

P*，MV，KC的定义见公式26。

FF的值可由罐体实际参数中附件种类数（NF）乘以每一种附件的损失因子（KF）计算。

对于浮盘附件，KFi可由公式37计算：

（

公式37）

式中：

KFi浮盘附件损失因子，磅-摩尔/年；

KFai无风情况下浮盘附件损失因子，磅-摩尔/年，

见表8；

KFbi有风情况下浮盘附件损失因子，

磅-摩尔/（迈m·年），见表8；

mi某类浮盘损失因子，无量纲量，见表8；

Kv附件风速修正因子，无量纲量（外浮顶罐，Kv=0.7；内浮顶罐和穹顶外浮顶罐，Kv=0）；

v平均气压平均风速，迈。

表8浮顶罐浮盘附件损失系数表

附件

状态

kfa

（磅-摩尔/年）

kfb

（磅-摩尔/

（迈n·年））

m

人孔

螺栓固定盖子，有密封件

1.6

0

0

无螺栓固定盖子，无密封件

36

5.9

1.2

无螺栓固定盖子，有密封件

31

5.2

1.3

计量井

螺栓固定盖子，有密封件

2.8

0

0

无螺栓固定盖子，无密封件

14

5.4

1.1

无螺栓固定盖子，有密封件

4.3

17

0.38

支柱井

内嵌式柱形滑盖，有密封件

33

/

/

内嵌式柱形滑盖，无密封件

51

/

/

管柱式滑盖，有密封件

25

/

/

管柱式挠性纤维衬套密封

10

/

/

取样管/井

有槽管式滑盖/重加权，有密封件

0.47

0.02

0.97

有槽管式滑盖/重加权，无密封件

2.3

0

0

切膜纤维密封（开度10%）

12

有槽导杆和取样井

无密封件滑盖（不带浮球）

43

270

1.4

有密封件滑盖（不带浮球）

无密封件滑盖（带浮球）

31

36

2.0

有密封件滑盖（带浮球）

有密封件滑盖（带导杆凸轮）

41

48

1.4

有密封件滑盖（带导杆衬套）

11

46

1.4

有密封件滑盖（带导杆衬套及凸轮）

8.3

4.4

1.6

有密封件滑盖（带浮球和导杆凸轮）

21

7.9

1.8

有密封件滑盖（带浮球、衬套和凸轮）

11

9.9

0.89

无槽导杆和取样井

无衬垫滑盖

31

150

1.4

无衬垫滑盖带导杆

25

2.2

2.1

衬套衬垫带滑盖

25

13

2.2

有衬垫滑盖带凸轮

14

3.7

0.78

有衬垫滑盖带衬套

8.6

12

0.81

呼吸阀

附重加权，未加密封件

7.8

0.01

4.0

附重加权，加密封件

6.2

1.2

0.94

浮盘支柱

可调式（浮筒区域）有密封件

1.3

0.08

0.65

可调式（浮筒区域）无密封件

2.0

0.37

0.91

可调式（中心区域）有密封件

0.53

0.11

0.13

可调式（中心区域）无密封件

0.82

0.53

0.14

可调式，双层浮顶

0.82

0.53

0.14

可调式（浮筒区域），衬垫

1.2

0.14

0.65

可调式（中心区域），衬垫

0.49

0.16

0.14

固定式

0

0

0

边缘通气阀

配重机械驱动机构，有密封件

0.71

0.1

1.0

配重机械驱动机构，无密封件

0.68

1.8

1.0

楼梯井

滑盖，有密封件

98

滑盖，无密封件

56

浮盘排水

/

1.2

注：

表中浮盘附件密封损失因子kra，krb，n只适用于风速6.8米/秒以下

（4）浮盘缝隙损失，ED。

螺栓固定的浮盘存在盘缝损失，由公式38计算：

（

公式38）

式中：

KD盘缝损耗单位缝长因子，0.14磅-摩尔/（英尺·年）；

SD盘缝长度因子，英尺/平方英尺，为浮盘缝隙长度与浮盘面积的比值，见表9；

D，P*，M