3. 污染源排放速率的时间维度处理
大气环境影响预测中,排放速率是模型最直接的输入之一,而它并不是一个固定不变的数。同一个烟囱,在正常生产、开停车、检修或事故状态下,排放的污染物浓度和流量可能相差数倍甚至数十倍。如果只用一个“平均排放速率”去代表全年,预测结果要么低估短时高浓度事件,要么高估长期平均影响。因此,排放速率必须按时间维度分层处理:先区分工况,再确定各工况下的排放速率,最后通过时间分配系数把年、日、小时量级衔接起来。
一、三类工况的界定与排放速率确定
正常工况指生产装置按设计参数稳定运行、污染治理设施同步正常投运的状态。此时排放速率通常由设计资料、物料衡算或实测数据确定。以某化工厂为例,其排气筒设计风量为 30000 m³/h,污染物 A 的排放浓度为 40 mg/m³,则正常工况排放速率为 30000×40×10⁻⁶ = 1.2 kg/h。这个数值是后续年排放量换算的基础。
非正常工况包括开停车、设备检修、治理设施故障但生产未停等情形。这类工况持续时间短但排放强度可能显著升高。确定方法一般是查阅装置操作规程中的开停车排放曲线,或按治理设施失效时的去除率下降来估算。例如上述装置在开停车阶段风量降至 20000 m³/h、浓度升至 120 mg/m³,则排放速率为 2.4 kg/h,是正常工况的 2 倍。非正常工况的排放速率应单独列出,不能与正常工况简单平均。
事故工况指泄漏、火灾、爆炸、安全阀起跳等突发情形。其排放速率往往难以精确预测,工程上常采用事故情景假设法:给定泄漏孔径、压力、温度,用源强估算模型(如泄漏速率公式或经验系数)计算短时释放速率。例如储罐接管断裂,按孔径 20 mm、压力 0.6 MPa 估算,初始泄漏速率可达数十 kg/h 量级。事故工况的排放速率用于风险预测和短时浓度评价,不参与年排放量统计。
二、年、日、小时排放量的换算关系
三种时间量级之间的换算,核心是“速率×时间×分配系数”。设正常工况小时排放速率为 Qh(kg/h),年运行小时数为 H(h/a),则年排放量 Qa = Qh×H。若装置年运行 8000 h,Qh=1.2 kg/h,则 Qa=9600 kg/a=9.6 t/a。
日排放量则由年排放量按生产日分配:Qd = Qa/D,D 为年生产天数。若年生产 330 天,则日均排放量约 29.1 kg/d。但要注意,日均排放量不等于日排放速率乘以 24 h,因为装置可能并非全天满负荷。更严谨的做法是:日排放量 = 当日各小时排放速率之和,即对小时序列积分。
小时排放量是模型输入的基本单位。对于稳态模型(如 AERMOD),通常输入的是年平均小时排放速率或逐时排放速率序列;对于短期浓度预测,则需要逐时排放速率。换算时须保持单位一致:浓度 mg/m³、风量 m³/h,相乘后除以 10⁶ 得 kg/h。
三、时间分配系数的实际用法
时间分配系数用于把年排放量拆分到月、日、小时,反映生产负荷和治理效率的季节性、昼夜性变化。常见做法是:先统计全年逐月产量或运行小时占比,得到月分配系数 fm(Σfm=1);再在月内按日分配系数 fd 拆分;最后按小时分配系数 fh 拆分到逐时。于是某小时排放速率 Qh,i = Qa×fm×fd×fh/Δt,其中 Δt 为该小时对应的时间长度(1 h)。
实际工作中,分配系数应来自企业生产报表、用电负荷曲线或在线监测数据,而不是凭空假设。若缺乏数据,可按“采暖季/非采暖季”“昼间/夜间”两段简化处理,但须在报告中说明假设依据。分配系数之和必须归一化,否则年排放量会失真。
四、一个可复用的计算示例
某水泥窑尾排气筒,正常工况风量 250000 m³/h,颗粒物浓度 20 mg/m³,年运行 7200 h,年生产 300 天。计算如下:
Qh = 250000 × 20 × 10⁻⁶ = 5.0 kg/h
Qa = 5.0 × 7200 = 36000 kg/a = 36 t/a
Qd = 36000 / 300 = 120 kg/d
若采暖季月系数 fm = 0.10(4个月合计0.40),
非采暖季 fm = 0.075(8个月合计0.60),
则采暖季某月排放量 = 36 × 0.10 = 3.6 t/月
把月排放量再除以该月生产天数得到日排放量,除以 24 得到平均小时排放速率,即可输入模型。若企业有逐时在线监测数据,应优先用实测序列替代平均分配,这样能保留夜间低负荷、昼间高负荷的真实波动。
五、处理原则小结
第一,工况先行:先判定正常、非正常、事故三类工况,分别确定源强,不混用。第二,量级对应:年排放量用于长期平均浓度和总量控制,日排放量用于日均浓度预测,小时排放量用于短期浓度和风险评价。第三,分配有据:时间分配系数应可追溯、可归一化,并在报告中说明来源。第四,单位统一:风量 m³/h、浓度 mg/m³、速率 kg/h 之间的 10⁻⁶ 换算不能漏。掌握这四点,排放速率的时间维度处理就能既满足模型输入要求,又经得起技术复核。