高少華，鄒 兵，崔積山，丁德武，朱勝杰

（1. 中國石化 青島安全工程研究院，山東 青島 266071；2. 環(huán)境保護部環(huán)境工程評估中心，北京 100012；3. 化學(xué)品安全控制國家重點實驗室，山東 青島 266071）

設(shè)備動靜密封點（以下簡稱“設(shè)備”）泄漏是石化企業(yè)主要揮發(fā)性有機物（VOCs）污染源之一［1-2］。國外實踐表明，泄漏檢測與修復(fù)（LDAR）是控制該污染源排放的有效手段。近期，隨著國家和各地VOCs排污收費辦法的相繼出臺和排污許可證管理工作［3］的快速推進，VOCs排放核算的準確性逐漸成為企業(yè)和政府關(guān)心的焦點［4］。目前多數(shù)企業(yè)采用相關(guān)方程法核算設(shè)備泄漏的排放量［5-7］。該方法與檢測數(shù)據(jù)緊密相關(guān)［8］。檢測VOCs的原理有多種，國家標準規(guī)定氫火焰離子化檢測儀（FID）為石化企業(yè)LDAR項目的檢測儀器［9-10］。FID可以檢測多數(shù)VOCs，具有基線穩(wěn)定、線性范圍寬、基本不受水蒸氣影響等特點。然而FID與其他原理的檢測儀器存在相似問題，即對不同組分的VOCs響應(yīng)并不一致，個別情況甚至差別較大。環(huán)保部2015年發(fā)布的《石化企業(yè)泄漏檢測與修復(fù)工作指南》（簡稱《工作指南》）給出了響應(yīng)因子的基本應(yīng)用方法［11］。然而，響應(yīng)因子隨被測氣體組分等多參數(shù)的變化及其在排放核算中的應(yīng)用至今未見詳盡報道。

本工作通過分析儀器型號、被測氣體種類和被測氣體濃度等多個因素與響應(yīng)因子的關(guān)系，給出應(yīng)用響應(yīng)因子的建議，并通過實例應(yīng)用說明響應(yīng)因子修正排放量的方法，希望對設(shè)備泄漏核算有一點借鑒作用。

1 響應(yīng)因子的定義

FID用氫火焰將進入檢測器的被測氣體電離，檢測器的極化電場使離子遷移而形成與氣體濃度成正比的微電流。該類儀器可以檢測多數(shù)VOCs，然而對不同的物質(zhì)，儀器靈敏度并不一致。為了表征儀器的這一特性，定義某儀器的響應(yīng)因子為：被測氣體摩爾分數(shù)的實際值與該儀器檢測值的比值，見式（1）。

式中：F為響應(yīng)因子，無量綱；xA為被測氣體的實際摩爾分數(shù)，10-6（μmol·mol-1）；xS為被測氣體摩爾分數(shù)的檢測值，10-6（μmol·mol-1）。

據(jù)美國環(huán)保署（EPA）給出的數(shù)據(jù)，煉油設(shè)備接觸的物料至少有66種VOCs（包括苯系物等38種有害空氣污染物），化工設(shè)備接觸的物料至少有378種為VOCs［12］。FID制造商提供了66種化合物的響應(yīng)因子［13］。針對部分化合物的氣體或蒸氣標準物質(zhì)在現(xiàn)階段難以得到的問題，Jalali-Heravi等［14］研究了運用神經(jīng)網(wǎng)絡(luò)和多元線性回歸預(yù)測FID響應(yīng)因子的方法，然而至今未見到實用報道。

2 響應(yīng)因子的相關(guān)因素分析

2.1 FID型號對響應(yīng)因子的影響

目前國內(nèi)石化企業(yè)使用的FID多數(shù)為進口儀器。表1為廠家給出的5種型號FID對正己烷、二甲苯和甲醇、甲苯、苯乙烯純物質(zhì)的響應(yīng)因子。

表1 不同型號FID的響應(yīng)因子

采用變異系數(shù)表征各FID的離散程度，計算方法見式（2）。

式中：ζj為檢測被測氣體j時，n種型號FID產(chǎn)生的變異系數(shù)，無量綱；Fij為i＃FID檢測被測氣體j的響應(yīng)因子，無量綱；為n種型號FID檢測被測氣體j的響應(yīng)因子的平均值。

將表1的數(shù)據(jù)分別代入式（2），可以計算出對于給定被測氣體j，響應(yīng)因子隨FID型號的變異系數(shù)。計算結(jié)果見表2。由表2可見：響應(yīng)因子隨檢測儀器的變異系數(shù)為49.42%～71.70%，說明對于表1中任何一種被測氣體，響應(yīng)因子均隨FID型號有顯著變化。

表2 響應(yīng)因子隨FID型號的變異系數(shù)

需要說明的是：來自同一制造商的FID的響應(yīng)因子差異并不大，可以計算出同一制造商的FID的變異系數(shù)范圍為0～19.17%；而不同制造商生產(chǎn)的FID的結(jié)構(gòu)及工作參數(shù)不同，各自的響應(yīng)因子差異也較大，這是響應(yīng)因子均隨FID型號有顯著變化的主要原因。如2＃FID檢測甲醇的響應(yīng)因子為3.19，3＃FID則為23.46。

2.2 被測氣體種類對響應(yīng)因子的影響

與2.1節(jié)類似，某FID檢測m種被測氣體時，響應(yīng)因子的變異系數(shù)可按式（3）計算。

式中：ζi為i＃FID分別檢測m種被測氣體時產(chǎn)生的變異系數(shù)，無量綱；為i＃FID分別檢測m種被測氣體時響應(yīng)因子的平均值。

將表1的數(shù)據(jù)代入式（3），得到響應(yīng)因子隨被測氣體種類的變異系數(shù)，見表3。由表3可見：FID檢測正己烷、二甲苯、甲醇、甲苯和苯乙烯等5種被測氣體組分的響應(yīng)因子變異系數(shù)為148.22%～175.46%，說明響應(yīng)因子隨被測氣體種類變化極為顯著。例如，3＃FID檢測甲醇的響應(yīng)因子高達23.46，檢測甲苯的響應(yīng)因子為0.95，前者是后者的24.69倍。

表3 響應(yīng)因子隨被測氣體種類的變異系數(shù)

2.3 被測氣體濃度對響應(yīng)因子的影響

1＃FID檢測苯系物、烯烴和烷烴的響應(yīng)因子隨被測氣體濃度的變化分別見圖1～圖3（數(shù)據(jù)由廠家提供）。

圖11#FID檢測苯系物的響應(yīng)因子隨被測氣體濃度的變化

由圖1～圖3可見：被測氣體的摩爾分數(shù)在（10～10000）×10-6范圍內(nèi)時，F(xiàn)ID檢測苯系物（二甲苯、甲苯和苯）的響應(yīng)因子從0.32～0.35降至0.10～0.12；檢測烯烴（C2～C4）的響應(yīng)因子從0.66～1.34降至0.24～0.74；檢測烷烴（C4，C6，C8）的響應(yīng)因子從0.35～0.58降至0.09～0.27。與2.1，2.2節(jié)類似，可以計算出：苯系物、烯烴和烷烴的響應(yīng)因子隨被測氣體濃度的變異系數(shù)分別為38.32%～40.04%、24.85%～37.72%和30.85%～49.58%。說明響應(yīng)因子不僅與被測氣體種類有關(guān)，而且與其濃度相關(guān)。對于上述3類物質(zhì)，F(xiàn)ID的響應(yīng)因子隨被測氣體濃度的增加而降低。

圖21#FID檢測烯烴的響應(yīng)因子隨被測氣體濃度的變化

圖31#FID檢測烷烴的響應(yīng)因子隨被測氣體濃度的變化

3 響應(yīng)因子的修正方法

根據(jù)《工作指南》，煉油裝置可不考慮FID響應(yīng)因子對檢測數(shù)據(jù)的影響，化工裝置則應(yīng)根據(jù)物料中VOCs的組分計算主要物料的響應(yīng)因子。因此，在實施LDAR項目前，首先需要辨識生產(chǎn)裝置的類別。

化工裝置在建立密封點數(shù)據(jù)庫過程中需要收集各主要單元的物料組分。按照標準規(guī)定的方法，獲取各物料的響應(yīng)因子［15］，并按《工作指南》要求，基于以下原則對檢測數(shù)據(jù)進行修正：1）F≤3時，檢測數(shù)據(jù)無需修正；2）3＜F≤10時，需要修正檢測數(shù)據(jù)；3）F ＞10時，選擇物料中F ＞10的氣體或響應(yīng)特性相近的氣體作為校準氣體，按標準規(guī)定的方法得出響應(yīng)因子。

化工裝置設(shè)備管線接觸的物料多數(shù)是由多種組分混合而成，對于多組分物料，可按式（4）計算物料的響應(yīng)因子（也稱為“復(fù)合響應(yīng)因子”）［16］。

式中：F＇為多組分物料的復(fù)合響應(yīng)因子；yk為被測氣體組分k的摩爾分數(shù)與物料中總烴的摩爾分數(shù)之比；Fk為被測氣體組分k的響應(yīng)因子。

4 響應(yīng)因子應(yīng)用實例

4.1 芳烴重整裝置

某石化企業(yè)芳烴重整裝置以直餾、加氫裂化石腦油或乙烯裂解汽油為原料，生產(chǎn)苯、對二甲苯和鄰二甲苯等芳烴產(chǎn)品，按照國標分類方法應(yīng)歸為化工。以檢出的145個接觸甲苯、二甲苯兩組分物料泄漏點為例，從圖1可以看出甲苯、二甲苯的響應(yīng)因子不超過3。按照式（4）計算，這兩種組分混合物的響應(yīng)因子也小于3。為了探討響應(yīng)因子對核算結(jié)果的影響，本文采用原始檢測數(shù)據(jù)（不考慮響應(yīng)因子）和響應(yīng)因子修正數(shù)據(jù)，分別核算甲苯及二甲苯混合物的排放量，并進行對比。

4.1.1 原始檢測數(shù)據(jù)法的核算結(jié)果

將145個泄漏點按照密封類型和物料狀態(tài)分類，閥門23個（氣體3個，輕液20個），法蘭97個（氣體6個，輕液91個），開口閥或開口管線13個（氣體1個，輕液12個），連接件12個（氣體1個，輕液11個），分別代入用《石化行業(yè)VOCs污染源排放排查工作指南》（簡稱《排查指南》）給出的相關(guān)方程，可以計算出145個泄漏點甲苯及二甲苯混合物的排放量為24525 kg/a。

4.1.2 響應(yīng)因子修正法的核算結(jié)果

為了簡化計算，假設(shè)物料中各組分的摩爾分數(shù)與泄漏后形成蒸氣的摩爾分數(shù)相同，即忽略兩組分揮發(fā)性差別，按式（4）計算物料的復(fù)合響應(yīng)因子。

由制造商給出的資料，檢測摩爾分數(shù)為500×10-6的甲苯和二甲苯的響應(yīng)因子分別為0.328和0.308，檢測摩爾分數(shù)10000×10-6的甲苯和二甲苯的響應(yīng)因子分別為0.114和0.104。為了簡化起見，采用《工作指南》給出方法，即選取較大的響應(yīng)因子進行修正。

將F1（甲苯）=0.328，F(xiàn)2（二甲苯）=0.308，以及甲苯和二甲苯的摩爾分數(shù)（各約為50%）代入式（4），可以求得F＇=0.318。

將檢測值和復(fù)合響應(yīng)因子計算值分別代入式（1），可得各密封點泄漏氣體的實際摩爾分數(shù)xA。將xA分別代入《排查指南》給出的相關(guān)方程，得出甲苯和二甲苯的排放量為4581 kg/a。響應(yīng)因子修正后的排放量僅為修正前的18.68%。

4.2 甲基叔丁基醚（MTBE）生產(chǎn)裝置

某企業(yè)MTBE裝置有387個接觸甲醇密封點。檢出泄漏點13個（以當(dāng)?shù)匦孤╅撝?00計算），其中閥門3個，法蘭和連接件4個，泵2個，開口閥4個。該MTBE裝置以混合C4和甲醇為原料，通過醚化合成MTBE。需要注意的是，該企業(yè)雖然以煉油生產(chǎn)為主，但按照國標分類方法［10］，其MTBE裝置為化工裝置。

從表1可以看出5種FID檢測甲醇的響應(yīng)因子均大于3，因此，需要對接觸甲醇設(shè)備和管線的檢測數(shù)據(jù)進行修正。表4給出檢測數(shù)據(jù)修正前后該裝置的甲醇排放量。由表4可見，響應(yīng)因子修正后的排放量是修正前的3倍左右。

表4 檢測數(shù)據(jù)修正前后某MTBE裝置的甲醇排放量

4.3 實例結(jié)果分析

從4.1節(jié)芳烴重整裝置145個甲苯和二甲苯泄漏點的FID測定響應(yīng)因子修正結(jié)果可以發(fā)現(xiàn)，修正后排放量計算值會降低。這是由于這兩種物質(zhì)在摩爾分數(shù)為500×10-6和10000×10-6時的響應(yīng)因子均小于1。4.2節(jié)MTBE裝置387個密封點，用響應(yīng)因子修正后的排放量計算值明顯增加，這主要是FID檢測甲醇響應(yīng)因子超過3造成的。

對比圖1～圖3可以發(fā)現(xiàn)，對于苯系物和烷烴，響應(yīng)因子修正后的排放量計算值會明顯下降，對于烯烴，下降幅度比前兩類化合物小。對于煉油設(shè)備特別是接觸苯系物、烷烴的設(shè)備，其排放量計算值用響應(yīng)因子修正后通常會有明顯下降；化工裝置涉及的物料組分較為繁雜，修正結(jié)果與物料具體組分有關(guān)。

4.4 應(yīng)用建議

為了應(yīng)用響應(yīng)因子修正核算結(jié)果，在LDAR項目建立過程中，除按照《工作指南》要求收集化工設(shè)裝置物料組分信息外，建議企業(yè)：1）通過收集工藝參數(shù)、現(xiàn)場采樣［17］和實驗室分析［18］等途徑全面收集煉油裝置涉苯系物、烷烴物料的組分信息；2）對于同一套裝置宜采用統(tǒng)一型號的FID，建立主要物料的響應(yīng)因子數(shù)據(jù)庫；3）采用逐點修正的方法，逐步實現(xiàn)設(shè)備泄漏排放的精細化核算。需要指出的是：至今制造商提供的FID基礎(chǔ)響應(yīng)因子數(shù)據(jù)庫有限，氣體標準物質(zhì)種類又難以覆蓋煉油和化工生產(chǎn)涉及的VOCs組分。建議研究機構(gòu)、儀器制造商和企業(yè)盡快聯(lián)手開展相關(guān)研究，開發(fā)VOCs響應(yīng)因子快速測定系統(tǒng)。

5 結(jié)語

a）從FID型號、被測氣體種類和被測氣體濃度等方面分析了影響FID響應(yīng)因子的多重因素。以某石化企業(yè)芳烴重整和MTBE裝置的部分設(shè)備為實例，計算了基于響應(yīng)因子的排放量。

b）FID檢測VOCs的響應(yīng)因子與FID型號、被測氣體種類和被測氣體濃度有關(guān)，變異系數(shù)分別為49.42%～71.70%、148.22%～175.46%、24.85%～49.58%。影響強度排序為：被測氣體種類＞FID型號＞被測氣體濃度。應(yīng)采用收集物料組分信息、建立基于檢測儀器的物料響應(yīng)因子數(shù)據(jù)庫和逐點修正的核算策略。

c）芳烴重整和MTBE裝置的實例應(yīng)用表明，響應(yīng)因子修正前后的排放量計算值變化明顯。響應(yīng)因子修正是設(shè)備泄漏排放的精細化核算的有效途徑之一。

致謝本研究得到環(huán)保部工程評估中心VOCs工作組童莉、周學(xué)雙、郭森、莊思源、李成福、方偉、沙沙、孫慧、賈萍、梁睿、羅霂、祝曉燕、冉麗君、何少林、王赫婧、楊一鳴等的大力支持，謹表謝意。

［1］ RTI International. Emissions Estimation Protocol for Petroleum Refineries［EB/OL］. 2015［2017-2-20］.https：//www3.epa.gov/ttn/chief/efpac/protocol/Protocol%20Report%202015.pdf.

［2］ 何少林，陳輝，于景琦，等. 石化廢水揮發(fā)性有機物逸散量的估算方法［J］. 化工環(huán)保，2015，35（6）：620 - 624.

［3］ 國務(wù)院辦公廳. 控制污染物排放許可制實施方案［EB/OL］. 2016［2017-3-25］. http：//zfs.mep.gov.cn/gz/bmhb/gwygf/201611/t20161124_368124.shtml.

［4］ 何少林，崔積山，王赫婧. 石化廢水揮發(fā)性有機物逸散管控分析［J］. 化工環(huán)保，2016，38（3）：65 - 69.

［5］ 朱勝杰，石昕，高少華，等. 加氫裂化裝置泄漏排放量評估研究［J］. 石油化工安全環(huán)保技術(shù)，2016，32（2）：29 - 31.

［6］ 丁德武，賈潤中，高少華，等. LDAR技術(shù)在石化企業(yè)應(yīng)用中常見問題解析［J］. 安全、健康和環(huán)境，2016，16（11）：39 - 42.

［7］ 劉桂霞，寧淼，楊清霞. 基于氣分裝置案例的LDAR減排效果影響因素識別［J］. 環(huán)境與可持續(xù)發(fā)展，2017，16（11）：166 - 169.

［8］ 環(huán)境保護部. 石化行業(yè)VOCs污染源排查工作指南［EB/OL］. 2015［2017-2-21］. http：//www.zhb.gov.cn/gkml/hbb/bgt/201511/W020151124546327744099.pdf.

［9］ 環(huán)境保護部，國家質(zhì)量監(jiān)督檢驗檢疫總局. GB 31570—2015石油煉制工業(yè)污染物排放標準［S］. 北京：中國環(huán)境科學(xué)出版社，2015.

［10］ 環(huán)境保護部，國家質(zhì)量監(jiān)督檢驗檢疫總局. GB 31571—2015石油化學(xué)工業(yè)污染物排放標準［S］. 北京：中國環(huán)境科學(xué)出版社，2015.

［11］ 環(huán)境保護部. 石化企業(yè)泄漏檢測與修復(fù)工作指南［EB/OL］. 2015［2017-2-22］. http：//www.zhb.gov.cn/gkml/hbb/bgt/201511/W020151124546328688845.pdf.

［12］ The U.S. Government Publishing Office. Standards of performance for equipment leaks of VOC in the synthetic organic chemicals manufacturing industry for which construction，reconstruction，or modification commenced after november 7，2006［EB/OL］. 2007［2017-2-25］. http：//www.ecfr.gov/cgi-bin/text-idx?S ID=f67100426c6dfabbaea2b7b9564a22c3&mc=true&n ode=sp40.7.60.vv_0a&rgn=div6.

［13］ Thermo Environmental Instruments Inc. Response factors［EB/OL］. 2000［2017-02-25］. http：//www.petersonenvironmental.com/ThermoTVA1000ResponseFactors.pdf.

［14］ Jalali-Heravi M，Garkani-Nejad Z. Prediction of relative response factors for flame ionization and photoionization detection using self-training artificial neural networks［J］. J Chromatogr A，9502002，950：183 -194

［15］ 環(huán)境保護部. HJ733—2014泄漏和敞開液面排放的揮發(fā)性有機物檢測技術(shù)導(dǎo)則［S］.北京：中國環(huán)境科學(xué)出版社，2014.

［16］ United States Environmental Protection Agency. Protocol for equipment leak emission estimates［EB/OL］.1995［2017-2-27］. https：//www3.epa.gov/ttn/chief/efdocs/equiplks.pdf.

［17］ 邵霞，李國傲，任培芳，等. 典型揮發(fā)性有機化合物在氣體采樣罐中的存儲穩(wěn)定性［J］. 化工環(huán)保，2017，37（1）：116 - 120.

［18］ 王興華，田英，任敏紅，等. 改進氣相色譜法測定總烴與非甲烷總烴［J］. 化工環(huán)保，2015，35（6）：662 - 666.