储罐呼吸排放量计算SH3002(2018)

1、储罐呼吸排放量的计算储罐呼吸排放量的计算 -设计设计 VOCs 治理设备治理设备-估算排放量的计算估算排放量的计算 一、一、 SH/T30022000石油库节能设计导则石油库节能设计导则中关于中关于拱顶罐、拱顶罐、 内浮顶罐“大呼吸” “小呼吸”排放量计算内浮顶罐“大呼吸” “小呼吸”排放量计算公式（公式（第第 2.2.1） 1，拱顶罐的拱顶罐的“大呼吸” “小呼吸”“大呼吸” “小呼吸”呼吸排放量呼吸排放量 “大呼吸”“大呼吸”排放量见公式排放量见公式（4） （） （5） （） （6） （） （7） （） （8） 式中：式中： L LDW DW: : 拱顶罐“大呼吸”排放量拱顶罐“大呼吸”排

2、放量 m /a；K KT T为周转系数，无量纲；为周转系数，无量纲； K K1 1油品系数，汽油油品系数，汽油=1=1、原油、原油=0.75=0.75，无量纲；，无量纲； P PY Y为油品平均温度下为油品平均温度下蒸气压蒸气压 KPaKPa；y y油蒸气相对分子量，无量纲；油蒸气相对分子量，无量纲； K K为单位换算常数为单位换算常数 5151.6.6，无量纲；，无量纲；V V1 1为泵送液体输入量为泵送液体输入量 m m； N N油罐周转次数，无量纲；油罐周转次数，无量纲；Q Q为油罐周转量，为油罐周转量，m m /a/a；V V为油罐体积，为油罐体积，m m； P PY1 Y1为油罐内液

3、面最低温度对应的蒸气压为油罐内液面最低温度对应的蒸气压 KPaKPa； P PY2 Y2为油罐内液面最高温度对应的蒸气压为油罐内液面最高温度对应的蒸气压 KPaKPa “小呼吸”排放量计算公式（“小呼吸”排放量计算公式（9） 式中：式中： L LDS DS: : 拱顶罐“拱顶罐“小小呼吸”排放量呼吸”排放量 m /a；K K2 2为单位换算为单位换算系数系数 3.053.05，无量纲；，无量纲； K K3 3油品系数，汽油油品系数，汽油=1=1、原油、原油= =0.580.58，无量纲；，无量纲； P P为为油罐内油品本体油罐内油品本体温度下蒸气压温度下蒸气压 KPaKPa；P Pa a为当地

4、大气压为当地大气压，KPaKPa； D D为为罐直径罐直径 m m；H H 为油罐内气体空间高度为油罐内气体空间高度 m m；T T大气温度的平均日温差大气温度的平均日温差，；，； F FP P为为涂料吸收，无量纲涂料吸收，无量纲；C C1 1为为小直径修正系数，无量纲。小直径修正系数，无量纲。 2 内浮顶内浮顶油罐的油罐的“大呼吸” “小呼吸”“大呼吸” “小呼吸”呼吸排放量呼吸排放量 “大呼吸”“大呼吸”排放量见公式（排放量见公式（10） 式中：式中： L LWF WF: : 为内浮顶为内浮顶罐“罐“大大呼吸”排放量呼吸”排放量 kg/a；Q Q1 1为油罐周转量为油罐周转量，1010 3

5、 3m m /a /a； C C 为油罐壁的黏附系数，为油罐壁的黏附系数，1010 - -3 3m m / /； ；P Py y为油品密度，为油品密度，kg/mkg/m； N NC C为为支柱个数，无量纲；支柱个数，无量纲；F FC C为支柱有效直径。为支柱有效直径。 “小呼吸”排放量见公式（“小呼吸”排放量见公式（11） （） （12） （） （13） 式中：式中： L LS SF F: : 内浮内浮顶罐“顶罐“小小呼吸”排放量呼吸”排放量 kg/a；K K8 8为单位换算为单位换算系数系数 0.450.45，无量纲；，无量纲； K Ke e为边缘密为边缘密封排放系数封排放系数，无量纲；，无

6、量纲；F FM M为浮盘附件总排放系数，无量纲；为浮盘附件总排放系数，无量纲； F Fd d为为顶板接缝长度顶板接缝长度系数，无量纲；系数，无量纲； K Kd d为为顶板接缝排放系数，焊接顶板顶板接缝排放系数，焊接顶板=0=0、非焊接顶板、非焊接顶板=3.66=3.66，无量纲；，无量纲； P P * *为蒸气压函数，无量纲； 为蒸气压函数，无量纲；m mV V为油气相对分子量，无量纲；为油气相对分子量，无量纲； K Kc c为为油品系数，汽油油品系数，汽油=1=1、原油、原油=0.4=0.4，无量纲；，无量纲； N Nmj mj为某种附件个数；为某种附件个数；K Kmjmj为某种附件排放系数

7、，无量纲。为某种附件排放系数，无量纲。 二二、其他其他的计算方法的计算方法和设计治理设备处理规模的思考和设计治理设备处理规模的思考 关于储罐关于储罐“大呼吸大呼吸” “” “小呼吸小呼吸”排放的计算，”排放的计算，国内外都开展了许国内外都开展了许 多研究工作， 其中较为得到公认并适合我国情况的有多研究工作， 其中较为得到公认并适合我国情况的有： 美国石油学会美国石油学会 （API）和西方石油）和西方石油学会（学会（WPI）的经验公式、）的经验公式、中石化系统（中石化系统（SHJ）也）也 编制了相应的经验型计算公式、编制了相应的经验型计算公式、 我国我国环保部有关部门编制的环保部有关部门编制的

8、石油化石油化 工行业工行业 VOCsVOCs 排放量计算办法 、排放量计算办法 、 美国美国 EPAEPA 等多个等多个计算方法计算方法。 这些计算。 这些计算 方法方法的特点，一是专业性很强，二是公式的特点，一是专业性很强，二是公式比较复杂比较复杂，三是需要引入三是需要引入不不 少因子数据很少因子数据很难查找难查找。美国环境保护署美国环境保护署(EPA)办公室的大气质量规划办公室的大气质量规划 和标准部门和标准部门(OAQPS)开发开发了了 TANKS 程序旨程序旨，用于用于估算油罐中有机液体估算油罐中有机液体 的大气排放。的大气排放。 对于我们环保一线工作者对于我们环保一线工作者， 可操作

9、性较低，可操作性较低， 也难以运用也难以运用。 环保部有关部门下发的环保部有关部门下发的石油化工行业石油化工行业 VOCs 排放量计算办法排放量计算办法 “二、二、 有机液体储存与调有机液体储存与调和挥发损失和挥发损失” ， 有有对固定顶对固定顶罐和罐和浮顶罐的 “呼浮顶罐的 “呼 吸” 排放的分别计算方法吸” 排放的分别计算方法。 固定顶罐固定顶罐 “呼吸” 损失包括 （“呼吸” 损失包括 （E固 固=ES+EW） ,ES 即即静置过程静置过程呼吸损失（ “小呼吸” ）呼吸损失（ “小呼吸” ）和和 EW即工作过程呼吸损失（ “大呼即工作过程呼吸损失（ “大呼 吸” ） 。浮顶罐“呼吸”损失

10、包括（吸” ） 。浮顶罐“呼吸”损失包括（E浮 浮=ER+EWD+EP+ED）, ER即边缘密封即边缘密封 损失、损失、EWD即挂壁损失、即挂壁损失、EP即浮盘附件损失、即浮盘附件损失、ED即缝隙损失。即缝隙损失。 综合各种计算方法， 考量 “大呼吸” 排放量综合各种计算方法， 考量 “大呼吸” 排放量， 储罐， 储罐的的储存油品储存油品 （或或 化工品化工品）的周转量（吞吐量）的周转量（吞吐量） 、周转次数等参数在计算中是、周转次数等参数在计算中是起起主要的主要的 影响参数；考量“小呼吸”损失，储罐尺寸、液体所占储罐体积和剩影响参数；考量“小呼吸”损失，储罐尺寸、液体所占储罐体积和剩 余空间

11、的比例余空间的比例大小、温度变化等参数在计算中是主要的影响参数。大小、温度变化等参数在计算中是主要的影响参数。 三、三、考量考量储罐“大呼吸”呼出量储罐“大呼吸”呼出量的相关标准的相关标准和简化公式和简化公式 确定确定 VOCs 治理治理治理设备的处理能力治理设备的处理能力，主要主要根据储罐在单位时间根据储罐在单位时间 内“大呼吸” “小呼吸”呼出内“大呼吸” “小呼吸”呼出 VOCs 的体的体最大最大积量积量（m /h） 。） 。至于“年至于“年 损耗”多少、或是多么大的储罐区等等参数，实际损耗”多少、或是多么大的储罐区等等参数，实际用不上。用不上。 SH/T30022000石油库节能设计导

12、则中关于拱顶罐、浮顶罐石油库节能设计导则中关于拱顶罐、浮顶罐 “大呼吸”呼出排放量的计算，以公式（“大呼吸”呼出排放量的计算，以公式（4）为例，）为例， V V1 1 前面前面K KT T、K K1 1、P PY Y、y y、K K、P PY Y等为等为无量纲无量纲类类系数系数，汽油，汽油=1=1、原油、原油=0.75=0.75， ；， ； 为油品平均温度下蒸气压为油品平均温度下蒸气压 K KPaPa；油蒸气相对分子量，无量纲；油蒸气相对分子量，无量纲； 为单位换算常数为单位换算常数 51.651.6，无量纲；，无量纲；V V1 1为泵送液体输入量为泵送液体输入量 m m； N N油罐周转次数

13、，无量纲；油罐周转次数，无量纲；Q Q为油罐周转量，为油罐周转量，m m /a/a；V V为油罐体积，为油罐体积，m m； 为油罐内液面最低温度对应的蒸气压为油罐内液面最低温度对应的蒸气压 KPaKPa； P PY2 Y2为油罐内液面最高温度对应的蒸气压为油罐内液面最高温度对应的蒸气压 KPaKPa 国家标准国家标准 油品装载系统油气回收设施设计规范油品装载系统油气回收设施设计规范 GB50759-2012 规定规定“油气回收装置的设计规模为最大装车体积流量的油气回收装置的设计规模为最大装车体积流量的1.01.1倍倍”。 石化石化标准石油化工储运系统罐区设计规范标准石油化工储运系统罐区设计规范

14、SH/T30072007， “液体进罐时的最大进液量所造成的罐体内液体蒸气呼出量液体进罐时的最大进液量所造成的罐体内液体蒸气呼出量，当液体，当液体 闪点（闭口）高于闪点（闭口）高于45时，应按最大进液量的时，应按最大进液量的1.07倍考虑；当液体闪倍考虑；当液体闪 点（闭口）低于或等于点（闭口）低于或等于45时，应按最大进液量的时，应按最大进液量的2.14倍倍。 估算估算储罐“大呼吸”储罐“大呼吸”油气排放量：按照统计得到的油罐车灌装汽油气排放量：按照统计得到的油罐车灌装汽 油（或其它轻质油品、易挥发化学品）的体积流量，参照相关标准确油（或其它轻质油品、易挥发化学品）的体积流量，参照相关标准确

15、 定估算蒸发气体体积的量（会大于进入定估算蒸发气体体积的量（会大于进入油罐内液体体积的量） ，确认油罐内液体体积的量） ，确认 一个系数。一个系数。实际运用中，实际运用中，不同不同环保科研环保科研、设计院、设计院、设备设备供应商对这个供应商对这个 系数的认知有所系数的认知有所不同，不同，业内业内有有 1.07、1.1、1.14、1.17、2.14 多个系数多个系数 的数据。的数据。我们确定我们确定 VOCs 治理设备处理能力，也须留有余量（设定一治理设备处理能力，也须留有余量（设定一 定操作弹性） 。定操作弹性） 。对于做环保装备而言，实际估算储罐呼吸排放量时，对于做环保装备而言，实际估算储罐

16、呼吸排放量时， 允许有一定允许有一定误差。误差。我们取我们取 1.2（参阅本文附三）（参阅本文附三） 可以简化排放量和确定可以简化排放量和确定油气油气回收处理装置的处理能力的估算回收处理装置的处理能力的估算， 并以此， 并以此 作为编制方案时选择油气回收处理规模的依据。作为编制方案时选择油气回收处理规模的依据。因此，我们认为，估因此，我们认为，估 算储罐算储罐“大呼吸”的呼出“大呼吸”的呼出排放量，确定治理设备处理能力的参数，可排放量，确定治理设备处理能力的参数，可 以以采用较为采用较为简化简化的的方法方法，主要确认储罐“呼出” （体积）量，主要确认储罐“呼出” （体积）量进行估算。进行估算。

17、 （参阅（参阅SH3002公式 （公式 （4） 中进料泵流量中进料泵流量V V1 1为主要为主要影响影响计算结果的计算结果的数据）数据） ， 因此，我们简化公式为因此，我们简化公式为： L LW W= =V V1 11.21.2 L LW W “大呼吸大呼吸”呼出量（气态）呼出量（气态）m /h， V V1 1见见 SH3002SH3002 公式公式（4 4） ，V V1 1为泵送液体输入量为泵送液体输入量 m m；我们定义我们定义为储罐为储罐单位时间单位时间进进 料量体积料量体积（液态）（液态）m /h 参考：下面为美国参考：下面为美国 EPA 对固定顶罐大呼吸损耗的计算。对固定顶罐大呼吸损

18、耗的计算。 （其中，周转量其中，周转量 Q Q 为动态数据，且是为动态数据，且是 影响呼吸量的主要因素影响呼吸量的主要因素） 设计设计 VOCs 治理设备的处理能力，无论储罐体积有多大、无论储罐治理设备的处理能力，无论储罐体积有多大、无论储罐 区有多少储罐、无论周转次数（及年吞吐量）有多少，我们只考量区有多少储罐、无论周转次数（及年吞吐量）有多少，我们只考量在在 单位时间内其进料作业造成“大呼吸”呼出排放量（单位时间内其进料作业造成“大呼吸”呼出排放量（m /h）和气温）和气温 温升造成“小呼吸”呼出排放量的最大数据，主要还是储罐进料泵的温升造成“小呼吸”呼出排放量的最大数据，主要还是储罐进料

19、泵的 流量。流量。 四四、储罐储罐“小呼吸”呼出量“小呼吸”呼出量 储罐储罐“小“小呼吸呼吸”排放是由于温度和大气压力的变化引起蒸气的排放是由于温度和大气压力的变化引起蒸气的 膨胀和收缩而产生的蒸气排出，它出现在罐内液面无任何变化的情膨胀和收缩而产生的蒸气排出，它出现在罐内液面无任何变化的情 况，是非人为干扰的自然排放方式。况，是非人为干扰的自然排放方式。参照下列公式参照下列公式： LB=0.191M （P/ （100910-P） ） ） 0.68D1.73H0.51 T0.45FPCKC 式中：式中： LB固定顶罐的呼吸排放量（固定顶罐的呼吸排放量（Kg/a） ；） ； M储罐内蒸气的分子量

20、；储罐内蒸气的分子量； P在大量液体状态下，真实的蒸气压力（在大量液体状态下，真实的蒸气压力（Pa） ；） ； D罐的直径（罐的直径（m） ；） ； H平均蒸气空间高度（平均蒸气空间高度（m） ；） ； T一天之内的平均温度差（一天之内的平均温度差（） ；） ； FP涂层因子（无量纲） ，涂层因子（无量纲） ，根据油漆状况取值在根据油漆状况取值在 11.5 之间；之间； C用于小直径罐的调节因子（无量纲）用于小直径罐的调节因子（无量纲）; KC产品因子（石油原油产品因子（石油原油 KC 取取 0.65，其他的有机液体取，其他的有机液体取 1.0） 。 以下附一、 二、 三， 分别补充对 “小呼

21、吸” 计算示例、 对以下附一、 二、 三， 分别补充对 “小呼吸” 计算示例、 对 SH/T3007 “热呼出”数据被误用为“小呼吸”呼出的分析、 “大呼吸”呼“热呼出”数据被误用为“小呼吸”呼出的分析、 “大呼吸”呼 出计算对取系数出计算对取系数 1.2 的思考。供参考。抛砖引玉，期望得到同行的思考。供参考。抛砖引玉，期望得到同行 经验指导。经验指导。请请赐赐给给 。 附附一一：案例：案例（计算储罐“小呼吸”排放量的（计算储罐“小呼吸”排放量的示例示例） 以四个规格甲醇储罐（10000、5000、3000、1000m）为例,估算“小呼吸”量： 标称容积 m D 直径 m H 高度 m 实际容

22、积 m 备注 10000 30 16 11300 5000 24 12.5 5600 3000 19 11.5 3200 1000 11.5 11 1100 LB=0.191M（P/（100910-P） ）0.68D1.73H0.51T0.45FPCKC M储罐内蒸气的分子量；甲醇取 32 P在大量液体状态下，真实的蒸气压力（Pa） ；用 30状态下饱和蒸汽压 22kpa D罐的直径（m） ；按照上表 H平均蒸气空间高度（m） ；取罐高度 1/2 T一天之内的平均温度差（） ；假设平均 15 FP涂层因子（无量纲） ，根据油漆状况取值在 11.5 之间； C用于小直径罐的调节因子（无量纲）;直

23、径在 09m 之间的罐体， C=1-0.0123(D-9)2 ; 罐径大于 9m 的 C=1； KC产品因子（石油原油 KC 取 0.65，其他的有机液体取 1.0） 计算取值： LB M P D H T FP C KC 10000 18434.14 32 22kpa 30 16 15 1.5 1 1 5000 11521.34 32 22kpa 24 12.5 15 1.5 1 1 3000 7370.771 32 22kpa 19 11.5 15 1.5 1 1 1000 3022.949 32 22kpa 11.5 11 15 1.5 1 1 计算结果： 【10000m储罐】 LB=0.

24、19132（22000/（100910-22000） ）0.68301.73160.51150.451.511 =6.1120.278798630.68241.7312.50.51150.451.511=18434.14 =18434kg（23042.68L/a） （注，这是（注，这是年排放量，若按照年排放量，若按照 360 天、每天天、每天 10 小时小时考虑，则有：考虑，则有： 3.2kg/h 的排放量，以密度的排放量，以密度 0.8 考虑，考虑，排放量为排放量为 6.4L/h） 【5000m储罐】 LB=0.19132 （22000/ （100910-22000） ） 0.68241.7

25、312.50.51150.451.511 =6.1120.278798630.68241.7312.50.51150.451.511=11521.34 =11521kg(14401L/a) （注，这是（注，这是年排放量，若按照年排放量，若按照 360 天、每天天、每天 10 小时小时考虑，则有：考虑，则有： 3.2kg/h 的排放量，以密度的排放量，以密度 0.8 考虑，考虑，排放量为排放量为 4L/h） 【3000m储罐】 LB=0.19132 （22000/ （100910-22000） ） 0.68191.7311.50.51150.451.511 =6.1120.278798630.6

26、8191.7311.50.51150.451.511=7370.771 =7371kg（9213.75L/a） （注，这（注，这是是年排放量，若按照年排放量，若按照 360 天、每天天、每天 10 小时小时考虑，则有：考虑，则有： 2.04kg/h 的排放量，以密度的排放量，以密度 0.8 考虑，考虑，排放量为排放量为 2.56L/h） 【1000m储罐】 LB=0.19132 （22000/ （100910-22000） ） 0.6811.51.7311.0.51150.451.511 =6.1120.278798630.6811.51.73110.51150.451.511=3022.94

27、9 =3023kg（3778.75L/a） （注，这是每年排放量，若按照（注，这是每年排放量，若按照 360 天、每天天、每天 10 小时小时考虑，则考虑，则 有：有：0.839kg/h 的排放量，以密度的排放量，以密度 0.8 考虑，考虑，排放量为排放量为 1.05L/h） 甲醇密度 0.7913，取 0.8 罐容积 LB（t/a） 排放量 m/a 每小时排放每小时排放 L/h 10000 1183414 14.79 6.4 5000 11.52134 14.56 4 3000 7.370771 9.315 2.56 1000 3.022949 3.82 1.05 储罐发生温升时的呼吸排放量

28、： 若按照若按照年年 360 天、每天天、每天 10 小时小时考虑， “小呼吸”排考虑， “小呼吸”排放量放量 10000m罐为 6.4L/h、5000m罐为 4L/h、3000m罐为 2.56L/h、1000m罐为 1.05L/h， 与“大呼吸”排放量相比较， “小呼吸”排放量很小。与“大呼吸”排放量相比较， “小呼吸”排放量很小。考量回收处理装置处考量回收处理装置处 理能力设计时，理能力设计时，适度适度考虑“小呼吸”排放量考虑“小呼吸”排放量，比较切合实际。，比较切合实际。 附二：关于附二：关于 SH30072014石油化工储运系统罐区设计规范石油化工储运系统罐区设计规范 5.1.6 条款

29、条款相关内容的理解和运用相关内容的理解和运用，（别别误用误用“热呼吸”“热呼吸” ） 关于如何考量储罐“大呼吸” “小呼吸”排放量，一段时间，我们在与设计关于如何考量储罐“大呼吸” “小呼吸”排放量，一段时间，我们在与设计 院或用户交流中，遇到一些将院或用户交流中，遇到一些将 SH30072014石油化工储运系统罐区设计规范石油化工储运系统罐区设计规范 5.1.6 条款作为估算排放量条款作为估算排放量， 因而给出“小呼吸” （， 因而给出“小呼吸” （热呼吸热呼吸）排放数据超大）排放数据超大的情况。的情况。 该条款规定：该条款规定： 5.1.6 通气管或呼吸阀的通气管或呼吸阀的通气量通气量，不

30、得小于下列各项的呼出量之和或吸入量之和：，不得小于下列各项的呼出量之和或吸入量之和： A）液体出罐时的最大出液量所造成的空气最大出液量所造成的空气吸入量吸入量，应按液体最大出液量考虑； B） “液体进入固定顶储罐时所造成的罐内液体气体罐内液体气体呼出量呼出量，当液体闪点（闭口）高于 45 时，应按最大进液量的 1.07 倍考虑；当液体闪点（闭口）低于或等于 45时，应按最大进 液量的 2.14 倍考虑。液体进入采用氮气或其他惰性气体密封保护系统的内浮顶储罐时所造 成的罐体内气体呼出量，应按最大进液量考虑； C）因大气最大温降导致罐内气体收缩所造成储罐吸入温降导致罐内气体收缩所造成储罐吸入的空气

31、量和因大气最大温升导致罐内温升导致罐内 气体膨胀而呼出气体膨胀而呼出的气体，宜按表 5.1.6 确定。 A) 吸入量，按液体最大出液量（体积）考虑； B) 呼出量，根据闪点温度，按最大进液量（体积）考虑； C) 因温降收缩造成吸入量和温升膨胀造成呼出量，即“热呼吸通气量” ，按上表。 （因有温 度变化的“小呼吸”特征，被当成小呼吸排 从表从表 5.1.6 看，热呼吸看，热呼吸通气量较大，通气量较大，如如 10000m储罐，吸入量和闪点储罐，吸入量和闪点37.8的呼出的呼出 量高达量高达 1210m /h，有的设计院和用户将此数据当成“小呼吸”排放量的依据，这是不确切，有的设计院和用户将此数据当成“小呼吸”排放量的依据，这是不确切 的。的。该表实际是该表实际是 呼吸阀通气管规格选取的规范呼吸阀通气管规格选取的规范（依据极限通气量，这是一个安全的要求）（依据极限通气量，这是一个安全的要求） ， 不能作为估算“小呼吸”排放量的依据。不能作为估算“小呼吸”排放量的依据。 附三：公式（附三：公式（4）至至（8）试算示例）试算示例 L LDW DW: : 拱顶罐“大呼吸”排放量拱顶罐“大呼吸”排放量 m /a；K KT T为周