密封點揮發(fā)性有機物（VOCs）排放速率的數(shù)值模擬及實測應(yīng)用

密封點揮發(fā)性有機物（VOCs）排放速率的數(shù)值模擬及實測應(yīng)用摘要：目前對設(shè)備與管線組件密封點揮發(fā)性有機物（VOCs）排放速率的討論尚停留閱歷估算的階段，本文以連接件密封點為典型代表，對乙醇生產(chǎn)裝置中所涉密封點VOCs排放速率進行數(shù)值模擬和實測應(yīng)用討論。利用直接測量包扎法對連接件進行包扎密封，將其所形成的包袋內(nèi)流通空間類比為石油化工反應(yīng)器，引入化工單元混合時間和空時的概念，通過計算流體力學數(shù)值模擬的方法，對比包袋內(nèi)速度場、濃度場和壓力場的分布狀況，定量分析出包扎法系統(tǒng)中惰性載氣吹掃流量和包袋體積之間的關(guān)系，優(yōu)選合理的操作參數(shù)，現(xiàn)場實測并推算出連接件VOCs排放速率。討論結(jié)果表明：適當減小包袋體積可增加混合強度，縮短混合時間；在保持袋內(nèi)微正壓條件下，合理選擇吹掃流量以匹配不同的包袋體積，吹掃時間大于袋內(nèi)混合空時的2倍可基本達到近穩(wěn)態(tài)；建立新的連接件密封點VOCs泄漏濃度和泄漏排放速率之間的冪函數(shù)關(guān)系式，用以評估核算VOCs年排放量更精確?????。揮發(fā)性有機物（VOCs）作為二次氣溶膠的重要參加物，導致了大氣污染物PM2.5、O3的生成。VOCs來自人為和自然?，工業(yè)源是主要人為排放源之一，而石化行業(yè)排放又占據(jù)肯定重要地位[1]。據(jù)美國環(huán)保署（EPA）估算，2000年美國石化設(shè)備泄漏產(chǎn)生70367t/a，占石化企業(yè)無組織VOCs總排放量的20%以上[2]。據(jù)國際石油工業(yè)環(huán)保協(xié)會估算，石化企業(yè)生產(chǎn)過程中工藝設(shè)備管線物料泄漏導致的VOCs排放量占全廠VOCs無組織排放量總量的40%～50%[3]。泄漏檢測與修復技術(shù)（LDAR）是利用便攜式監(jiān)測設(shè)備，檢測化工企業(yè)各類設(shè)備與管線組件上閥門、法蘭及連接件等易產(chǎn)生泄漏的密封點，并修復超過肯定濃度的泄漏點，從而達到掌握原料泄漏損失和環(huán)境愛護的目的。隨著排污許可證管理和環(huán)境愛護稅的快速推動，VOCs排放核算的精確?????性漸漸成為企業(yè)和政府關(guān)懷的焦點[4]。對石化行業(yè)VOCs進行收稅的前提是評估核算出VOCs年排放量，目前多采納單個泄漏點排放濃度的閱歷公式法[5]，主要套用EPA、歐盟等國外石化無組織排放因子對石化企業(yè)裝置組件的無組織排放進行估算，考慮到各國裝置操作參數(shù)、運行條件、使用原料、運行管理等方面差異，這種借用國外無組織排放因子的做法無助于精確?????核算VOCs排放量[6]。石化工業(yè)排污許可中對設(shè)備與管線組件密封點泄漏揮發(fā)性有機物排放量的填報則是結(jié)合國內(nèi)裝置密封點的統(tǒng)計狀況及中國石化集團公司、中國石油集團公司及中國海油集團公司完成的第一輪LDAR的檢測結(jié)果，推算出有機氣體和揮發(fā)性有機液體流經(jīng)的設(shè)備與管線組件總數(shù)的0.3%，以泄漏濃度大于等于10000μmol/mol對于各種密封點的定值泄漏系數(shù)計算泄漏揮發(fā)性有機物的量，認為與企業(yè)實際排放量基本吻合[7]。從目前看，對設(shè)備與管線組件密封點VOCs排放速率的討論尚停留在數(shù)學統(tǒng)計和閱歷估算的階段，對組件泄漏VOCs排放過程及排放量尚缺乏實質(zhì)性了解和深化討論。本討論通過計算流體力學數(shù)值模擬的方法，選取石化行業(yè)乙醇生產(chǎn)裝置的連接件密封，定量分析VOCs排放速率，直接測量包扎法系統(tǒng)中惰性載氣吹掃流量和包袋體積之間的關(guān)系，摸索直接測量法實測系統(tǒng)操作規(guī)律，利用包扎法進行現(xiàn)場試驗，結(jié)合LDAR相關(guān)方程推導出VOCs泄漏排放速率，建立一套設(shè)備與管線組件密封點VOCs排放速率的系統(tǒng)討論方法，為VOCs真實排放量的核算供應(yīng)方法和參考依據(jù)。后續(xù)可結(jié)合數(shù)值模擬法和現(xiàn)場包扎法對設(shè)備與管線組件其他密封點類型進行VOCs泄漏討論，特殊是結(jié)構(gòu)簡單的密封類型，如閥門、泵，可充分發(fā)揮數(shù)值模擬的技術(shù)特點和優(yōu)勢，并借鑒化工單元合適的參數(shù)來優(yōu)化設(shè)計現(xiàn)場試驗方案，推導各類密封點VOCs排放速率的核算公式，獵取更精確?????的年排放量，為建立切合我國石化行業(yè)的VOCs排放因子計算標準和評估體系奠定基礎(chǔ)。1泄漏排放速率的評估方法1.1評估方法石化行業(yè)泄漏損失評估的方法主要是依據(jù)EPA泄漏損失評估協(xié)議，包括直接測量法、LDAR相關(guān)方程法、篩分法和平均排放系數(shù)法[8,9]。這四種方法中，直接測量法是通過直接測量泄漏氣量，收集并分析泄漏氣體組成，其精確?????性和牢靠性最高。但由于石化企業(yè)的設(shè)備與管線組件過多，直接測量法的現(xiàn)場工作量巨大，實際操作性不強，應(yīng)用較多的是LDAR相關(guān)方程法。該方法基于LDAR定量檢測每一受控密封點泄漏濃度值，將其凈值SV（μmol/mol）代入相應(yīng)關(guān)系方程得出排放速率eVOCs（kg/h）來估算年排放量[10]。相關(guān)方程法是將泄漏濃度值利用數(shù)學公式轉(zhuǎn)化為排放速率，中間避開了時間變量，精確?????性必定受到各裝置操作參數(shù)和運行介質(zhì)、設(shè)備和管線組件的密封材質(zhì)以及設(shè)備運行管理等各個方面的影響，應(yīng)依據(jù)石化企業(yè)的實際狀況利用直接測量法加以修正，重點考慮我國石化生產(chǎn)裝置的本地化特征，以削減評估核算誤差。1.2包扎法直接測量法一般采納包扎法，將排放密封點或排放口用薄膜包扎起來形成一個密封包袋，讓已知流量的惰性載氣通入，待載氣達到平衡后，收集氣體樣品后測量總有機物濃度，也可以針對氣樣中的單個化合物濃度進行分析后加和，通過樣品濃度和載氣流量推算VOCs的排放速率。包扎法操作系統(tǒng)流程示意圖見圖1。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359346347.jpg”alt=“1.jpg”width=“915”height=“520”/設(shè)備和管線組件的密封點可分為動靜密封點，共9大類[2]，因連接件結(jié)構(gòu)規(guī)整，將連接密封點作為目標對象。惰性載氣選用N2，現(xiàn)場用高純氮氣經(jīng)減壓后進入包扎系統(tǒng)。為削減材質(zhì)對VOCs成分的影響，全部連接管均采納特氟龍管。包袋前后分別設(shè)置一個流量計，如前后流量基本全都，說明包扎嚴密好。N2經(jīng)三通閥分兩側(cè)進入袋內(nèi)，便于對泄漏VOCs氣體的混合和攜帶。通過穩(wěn)定性檢測推斷N2和VOCs混合程度，現(xiàn)場選用便攜式有機物檢測儀進行VOCs濃度測量，數(shù)據(jù)變化在±5%時，認定袋內(nèi)氣體混合勻稱、穩(wěn)定[6]。測定袋內(nèi)氣體溫度后進行采樣，將混合氣體充入空白采樣袋中，樣品即可送色譜分析，減去背景值得到凈濃度值。1.3排放速率的計算依據(jù)物料衡算，包扎法估算密封點VOCs泄漏的排放速率如式(1)。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359341891.jpg”alt=“2.jpg”width=“668”height=“47”/式中，氣體常數(shù)轉(zhuǎn)換系數(shù)為1.219×10-5(K·kmol)/m3；體積濃度分數(shù)轉(zhuǎn)換系數(shù)為106μmol/mol。2計算流體力學分析2.1數(shù)值模擬建模計算流體力學數(shù)值模擬是采納有限體積法對多組分物質(zhì)在流場空間內(nèi)的質(zhì)量守恒方程和動量守恒方程進行離散求解一種數(shù)學方法。數(shù)值模擬可為現(xiàn)場試驗供應(yīng)有效的數(shù)據(jù)關(guān)心，優(yōu)化試驗方案。依據(jù)式(1)，包扎法系統(tǒng)中僅吹掃流量是自變量，其他則是化工系統(tǒng)本身參數(shù)，和包扎操作無關(guān)。另外，包袋體積的大小與包扎方法、設(shè)備和管線組件密封點類型和大小有關(guān)，是包扎系統(tǒng)的一個操作變量。連接件包扎后外形規(guī)章，類似化工中的套管反應(yīng)器，連接件在中部，VOCs向上發(fā)生泄漏。在同側(cè)，N2分兩支進入包袋，呈V字形。N2與VOCs充分混合后從對側(cè)尾部出口流出，數(shù)值模擬三維建模的物理結(jié)構(gòu)及相關(guān)設(shè)置見圖2。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359352546.jpg”alt=“2.jpg”width=“915”height=“285”/對不同包袋體積均采納三維六面體結(jié)構(gòu)化網(wǎng)格建模，如圖2所示。管道外徑d1=20mm，包袋直徑d2分別為24mm/30mm/40mm。內(nèi)部流場中設(shè)三個監(jiān)測點位：P1在上游N2進口對側(cè)、P2在中部VOCs泄漏對側(cè)、P3在下游出口處。2.2不同包袋體積的混合時間混合時間是進料在全反應(yīng)器范圍內(nèi)分布勻稱實際上所需要的時間，一般采納示蹤劑脈沖注入反應(yīng)器后達到最終穩(wěn)定±5%的時間即為試驗測定的混合時間tm[12]。將VOCs類比為示蹤劑，數(shù)值模擬過程中采納物質(zhì)輸運方程求解，通過VOCs的對流集中方程計算其在混合氣體中的摩爾分數(shù)[13]。通過非穩(wěn)態(tài)時間計算，監(jiān)測出口P3點VOCs濃度隨時間的變化。在吹掃流量Q=6.03×10-2m3/h的條件下，不同包袋體積VOCs濃度隨時間變化趨勢及其混合時間見圖3，其中對混合濃度進行無因次處理表示為混合程度。參考抱負反應(yīng)器的設(shè)計方程，當量直徑de=d2-d1，空時τ=V/Q，不同包袋體積對應(yīng)的相關(guān)流體力學參數(shù)見表1。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359355377.jpg”alt=““width=“915”height=“556”/<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359357080.jpg”alt=““width=“912”height=“179”/由圖3和表1可見，隨著包袋體積的增加，VOCs在袋內(nèi)混合時間加長。若緊貼近管道進行包扎，形成的包袋體樂觀小，混合時間tm僅為1.3s。同樣吹掃流量條件下，包袋體積比為6.8（9.42/1.38）倍時，混合時間的倍數(shù)略高為7.4（9.6/1.3）倍。隨體積的增加，混合空間擴大，返混程度加劇，混合勻稱需要更長的時間。空時是抱負反應(yīng)器宏觀混合狀況的反映，宏觀流場混合勻稱后，不同的物料組分需要在微觀尺度進一步對流集中才可達到最終的混合勻稱。表1數(shù)據(jù)表明，混合時間均大于空時，而且體積越大，倍數(shù)越高，呈非線性增加。EPA建議包袋體積盡可能小，密封空間減小，可保證混合時間最短，且由于停留時間和傳熱表面積削減可削減或防止內(nèi)部重組分有機化合物的凝聚[11]。由于現(xiàn)場設(shè)備與管線組件的簡單狀況，特殊是對于動密封點、閥門類簡單結(jié)構(gòu)和大設(shè)備的包扎，很難保證包袋體樂觀小。如包袋體積較大，可適當增加吹掃流量，以降低空時保證宏觀混合強度，并延長吹掃時間直到滿意混合時間的要求，吹掃時間在空時的兩倍以上基本可達到近穩(wěn)態(tài)。2.3不同物料的混合狀況考慮到現(xiàn)場很難具備緊貼設(shè)備和管線組件進行抱負包扎的實施條件，選擇中體積為后續(xù)討論對象，考察包袋內(nèi)不同點位的物料濃度動態(tài)變化狀況，O2和N2的體積分數(shù)隨時間的變化規(guī)律分別見圖4(a)和(b)。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359352664.jpg”alt=““width=“738”height=“516”/<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359352859.jpg”alt=““width=“792”height=“580”/包扎完成后，充入N2前，包袋內(nèi)初始狀態(tài)為空氣，O2和N2的體積分數(shù)分別占21%和79%，VOCs為0。假設(shè)包扎密封好，不存在空氣進入袋內(nèi)的狀況，氧氣則為完全滅失氣體，如圖4(a)所示。其濃度在各個點位逐次被氮氣吹掃出，濃度緩慢降低，各個點位的濃度最終逐步衰減為0。P1～P3在包袋內(nèi)從上游到下游逐次分布，O2滅失也是依此挨次逐步衰減。如圖4(b)所示，N2為主流氣體，濃度逐次、逐步增大，挨次和O2滅失挨次全都。由于N2從VOCs泄漏的上游吹入，VOCs很難進入上游和其對側(cè)繞流區(qū)域，這些區(qū)域基本被N2漸漸布滿，N2體積分數(shù)接近100%。在下游區(qū)域，N2和VOCs進行混合、攜帶后流出，由于VOCs占有肯定比例，所以P3點的N2體積分數(shù)明顯低于P1和P2點。2.4速度場、濃度場和壓力場對比分析在中體積包袋中，在吹掃流量Q=6.03×10-2m3/h的條件下，選VOCs泄漏下游中間截面，速度場、濃度場（VOCs）和壓力場（總壓）如圖5所示。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359366716.jpg”alt=“3.jpg”width=“917”height=“300”/如圖5(a)速度場所示，因出口位于下方，中下游上部區(qū)域速度開頭銳減，由于出口縮徑的影響，速度在接近出口區(qū)域逐步增加；圖5(b)壓力場顯示整個包袋內(nèi)總壓力約在800Pa，分布較勻稱，為微正壓狀態(tài)。由于動壓頭的影響，在出口處總壓力略高；圖5(c)VOCs的濃度場顯示，在中下游上部區(qū)域VOCs濃度相對較高，逐次往出口移動。對比三場，可以發(fā)覺速度場和壓力場相關(guān)性極高，而濃度場主要和VOCs泄漏源的位置以及吹掃氣體的出入口位置有關(guān)。在包扎法實測操作中，吹掃氣體的出入口位置設(shè)置至關(guān)重要。2.5不同吹掃流量的VOCs濃度變化對中體積包袋，將吹掃流量依次增大4倍，即Q2=4×Q1、Q3=16×Q1，在這3個流量條件下，對應(yīng)的VOCs濃度在P3點隨時間的變化趨勢見圖6，為便利對比，對數(shù)據(jù)進行4倍處理，3條曲線數(shù)值分別為GC1、GC2×4、GC3×16。<imgsrc=“/UploadFiles/2020001/20209/Env/202009052359369556.jpg”alt=“5.jpg”width=“918”height=“540”/泄漏VOCs可視為連續(xù)型示蹤劑，初始狀態(tài)時濃度為0，泄漏開頭后濃度漸漸增加，如圖6所示。隨著吹掃流量增至Q2，包袋內(nèi)的混合強度明顯增加，混合時間快速減小。流量進一步增至Q3后，袋內(nèi)消失劇烈湍流，造成出口P3點區(qū)域的濃度脈動。依據(jù)式(1)，吹掃流量增大，相當于稀釋了樣品VOCs濃度，GC值則對應(yīng)降低4倍。由圖6可見，4倍處理后最終濃度值并不相同，而是逐步遞減，流量增加后反而簡單造成排放速率eVOCs的核算偏低。這主要是由于流量增大，造成袋內(nèi)湍流脈動增加，袋內(nèi)壓力過高，轉(zhuǎn)變袋內(nèi)氣體的密度。EPA特殊強調(diào)袋內(nèi)壓力不應(yīng)超過6894Pa[11]，壓力過高，簡單造成包扎材料的破損和對包扎密封的破壞。泄漏點的正常泄漏狀態(tài)時表壓為0，假如包袋內(nèi)壓力過高，轉(zhuǎn)變了包扎前密封點原始泄漏狀態(tài)，簡單造成實測數(shù)據(jù)失真。3包扎法實測分析在石化行業(yè)乙醇生產(chǎn)裝置現(xiàn)場進行包扎法直接測量分析，密封點的泄漏值SV范圍為0～8×103μmol/mol，可涵蓋泄漏值2×103μmol/mol的認定標準[14]?，F(xiàn)場實測中重復性有：①現(xiàn)場密封點的包扎次數(shù)；②每次包扎密封點的采樣次數(shù)；③一個采樣袋的測試次數(shù)。各組重復次數(shù)均為3，基本可保證明測數(shù)據(jù)的重復概率。為和LDAR數(shù)據(jù)及其相關(guān)方程曲線對比，選用冪函數(shù)形式擬合，詳細擬合曲線及相關(guān)數(shù)據(jù)見圖7，因相關(guān)方程法的