• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    紅外掩日遙感監(jiān)測(cè)煉油廠非甲烷烷烴排放通量

    2022-07-19 03:37:04李凌波程夢(mèng)婷劉新宇
    中國環(huán)境科學(xué) 2022年7期
    關(guān)鍵詞:煉油廠烷烴通量

    李凌波,宮 超,程夢(mèng)婷,李 龍,劉新宇

    紅外掩日遙感監(jiān)測(cè)煉油廠非甲烷烷烴排放通量

    李凌波1*,宮 超2,程夢(mèng)婷1,李 龍1,劉新宇1

    (1.中國石油化工股份有限公司,大連石油化工研究院,遼寧 大連 116045;2.中國石油化工股份有限公司安慶分公司,安徽 安慶 246002)

    采用紅外掩日通量遙感監(jiān)測(cè)(SOF)技術(shù)分別在2014年5月~2015年12月和2021年10月監(jiān)測(cè)了我國7座大型煉油廠(其中6座原油年加工量均超過1000萬t)非甲烷烷烴排放通量(kg/h)及分布.每座煉油廠監(jiān)測(cè)3~8d,測(cè)量18~73次,總計(jì)獲得328個(gè)排放通量測(cè)量數(shù)據(jù),根據(jù)國內(nèi)煉油廠VOCs排放煙羽中非甲烷烷烴質(zhì)量分?jǐn)?shù)估算了VOCs排放量,結(jié)合監(jiān)測(cè)期間實(shí)際原油加工量計(jì)算了非甲烷烷烴和VOCs排放系數(shù).結(jié)果顯示: 7座煉油廠2014~2015年的非甲烷烷烴排放系數(shù)測(cè)定值為0.016%~0.11%,平均為0.081%;VOCs排放系數(shù)估算為0.020%~0.14%,平均為0.10%.無組織排放約占煉油廠非甲烷烷烴排放總量的70%以上,其中輕油貯罐排放占比過半.國內(nèi)7座煉油廠2014~2015年非甲烷烷烴排放系數(shù)的最好水平與美國南加州6座煉油廠同期SOF監(jiān)測(cè)的最好水平相當(dāng),但非甲烷烷烴排放系數(shù)的平均水平約為其平均水平的3.9倍,國內(nèi)煉油廠的VOCs排放控制水平更加參差不齊.國內(nèi)1座千萬噸級(jí)煉油廠2021年監(jiān)測(cè)的非甲烷烷烴排放通量和排放系數(shù)分別較2015年削減72.4%和74.2%.SOF可為石化VOCs無組織排放監(jiān)測(cè)、量化和排放清單修訂提供最佳實(shí)用技術(shù),本研究結(jié)果提供了國內(nèi)石化行業(yè)VOCs綜合整治初期典型煉油廠非甲烷烷烴和VOCs的基線排放實(shí)測(cè)數(shù)據(jù),以及1座千萬t級(jí)煉油廠6a后治理攻堅(jiān)效果.

    煉油廠;非甲烷烷烴;VOCs;無組織排放;排放通量監(jiān)測(cè);紅外掩日通量遙感監(jiān)測(cè)

    石化工業(yè)是僅次于溶劑使用的第二大揮發(fā)性有機(jī)物(VOCs)人為排放源,在國內(nèi)人為VOCs排放的占比逐年增加[1].歷史原因和城市化進(jìn)程的加快已使國內(nèi)的很多石化企業(yè)瀕臨城市或被城市包圍,VOCs及其伴生異味的排放控制已成為其可持續(xù)發(fā)展的關(guān)鍵.石化企業(yè)VOCs排放源包括設(shè)備密封隨機(jī)泄漏、貯罐呼吸與泄漏、產(chǎn)品裝卸揮發(fā)等十二類[2],其中大部分為無組織排放源,隨著VOCs有組織排放逐步有效控制,無組織排放占據(jù)VOCs排放的主導(dǎo),正成為VOCs排放控制的重點(diǎn)[3-4].石化VOCs無組織排放源數(shù)量多、分散且排放無規(guī)則,部分無規(guī)則排放口,部分為非穩(wěn)態(tài)和隨機(jī)泄漏排放,污染物排放的時(shí)間和空間分布波動(dòng)較大,總體表現(xiàn)為大型面源和體積源特征,與氣象條件復(fù)合形成復(fù)雜的污染特征,排放監(jiān)控難度很大,排放通量監(jiān)測(cè)與核算尤為困難和復(fù)雜.國際上廣泛應(yīng)用美國EPA AP-42排放因子/模型核算石化企業(yè)VOCs排放[5],基于該方法估算的美國煉油廠VOCs排放系數(shù)約為0.01%~0.02%[6-7],約為正常生產(chǎn)工況下實(shí)際監(jiān)測(cè)值的1/5~1/3[3,8],可能低估排放.目前國內(nèi)主要參考美國EPA AP42排放因子/模型核算石化企業(yè)VOCs排放,然而核算的煉油廠VOCs排放系數(shù)為0.182%[4,9],遠(yuǎn)高于國際主流數(shù)據(jù),約為美國的10倍,可能高估了實(shí)際排放.該方法主要基于計(jì)算,缺少實(shí)測(cè)檢驗(yàn)及本土化驗(yàn)證與修訂,還存在排放因子質(zhì)量不佳、計(jì)算參數(shù)選擇人為因素影響大、方法不確定度高等問題,難以得到客觀、準(zhǔn)確、可靠的VOCs排放數(shù)據(jù),導(dǎo)致排放清單不確定性較高,困擾VOCs排放控制決策、監(jiān)管、管控及治理效果量化評(píng)估[3].雖然衛(wèi)星遙感監(jiān)測(cè)和大氣環(huán)境監(jiān)測(cè)可在一定程度校驗(yàn)和改進(jìn)VOCs排放清單的精度和分辨,但受限于VOCs排放種類繁多、反應(yīng)活性和生命期不一、時(shí)空分布不均勻、模型精度等因素.有必要發(fā)展VOCs無組織排放通量廠界或近源監(jiān)測(cè)技術(shù),修訂完善排放因子/排放模型,提高VOCs排放清單的精度和分辨,推動(dòng)石化VOCs排放控制與治理由粗放向精準(zhǔn)化發(fā)展.

    VOCs無組織排放通量監(jiān)測(cè)技術(shù)有間接測(cè)量和遙感監(jiān)測(cè)兩類[3].間接測(cè)量技術(shù)主要包括示蹤技術(shù)(TCT)和反向擴(kuò)散模型(RDM)等[10],遙感監(jiān)測(cè)技術(shù)主要包括徑向羽映射(RPM)、差分吸收激光雷達(dá)(DIAL)和紅外掩日通量(SOF)[10-12]等技術(shù).間接測(cè)量技術(shù)相對(duì)簡(jiǎn)單,可用于小型無組織排放源通量測(cè)算,但難以準(zhǔn)確可靠地監(jiān)測(cè)石化企業(yè)這類大型面源或體積源排放通量,RDM通過下風(fēng)向多點(diǎn)監(jiān)測(cè)VOCs濃度并反推無組織排放源強(qiáng)的方法偏差較大,TCT的準(zhǔn)確性優(yōu)于RDM,但也受限于示蹤排放和混合條件,監(jiān)測(cè)不確定性較高.SOF、DIAL和RPM等遙感監(jiān)測(cè)技術(shù)可直接測(cè)量VOCs總體排放及分布,其中RPM光路覆蓋范圍有限,僅適于小型面源監(jiān)測(cè),SOF和DIAL可定量監(jiān)測(cè)石化企業(yè)等大型面源和體積源VOCs排放.SOF和DIAL監(jiān)測(cè)石化企業(yè)VOCs排放應(yīng)用超過20a,技術(shù)較為成熟,歐洲標(biāo)準(zhǔn)委員會(huì)(CEN)已頒布基于DIAL和SOF等技術(shù)的VOCs無組織排放監(jiān)測(cè)標(biāo)準(zhǔn)方法[13-14].DIAL可全天候監(jiān)測(cè),監(jiān)測(cè)半徑1km左右,測(cè)量精度更高,并定位排放源,但設(shè)備龐大,操作復(fù)雜,監(jiān)測(cè)費(fèi)用較高,上風(fēng)向輸入難以同步扣除[3].SOF監(jiān)測(cè)范圍可覆蓋石化企業(yè)廠區(qū)總體及廠內(nèi)功能區(qū)或石化園區(qū)VOCs排放,結(jié)合生產(chǎn)工況和氣象條件優(yōu)化并實(shí)施大量監(jiān)測(cè),通過統(tǒng)計(jì)分析獲取代表性排放數(shù)據(jù),監(jiān)測(cè)時(shí)需要陽光和適宜的風(fēng)場(chǎng)等天氣條件[3].

    非甲烷烷烴是國內(nèi)大氣含量最豐富的非甲烷烴(NMHC)或VOCs物種[15],也是煉油企業(yè)排放最多的VOCs,貢獻(xiàn)煉油廠VOCs排放總量的大部分[16-19].雖然非甲烷烷烴的化學(xué)活性不及同為煉油特征VOCs的苯系物和烯烴,但由于排放量遠(yuǎn)高于苯系物和烯烴,對(duì)煉油廠周邊區(qū)域大氣的臭氧生成潛勢(shì)也有重要貢獻(xiàn)[16-17].作為煉油廠VOCs排放的主體組分,非甲烷烷烴也是煉油廠VOCs排放總量測(cè)算的基礎(chǔ),可為其他特征VOCs的排放測(cè)算提供重要參比.

    本文采用SOF技術(shù)在2014~2015年監(jiān)測(cè)了7家大型煉油企業(yè)非甲烷烷烴無組織排放通量及分布,并在2021年復(fù)測(cè)了其中1家煉油企業(yè),估算了VOCs無組織排放通量,以及非甲烷烷烴、VOCs排放速率和排放系數(shù),并與美國加州煉油廠同期同類監(jiān)測(cè)結(jié)果進(jìn)行了對(duì)比,可用于評(píng)估煉油企業(yè)VOCs排放管控治理水平與效果,校正/修訂煉油廠VOCs排放清單,或創(chuàng)建基于實(shí)測(cè)的排放清單,也可作為國內(nèi)石化工業(yè)VOCs無組織排放控制初期典型煉油企業(yè)VOCs基線排放.

    1 材料與方法

    1.1 測(cè)量?jī)x器

    SOF儀器由瑞典FluxSense公司提供,主要由傅里葉變換紅外光譜儀、太陽跟蹤器、光路傳輸、GPS定位、數(shù)據(jù)處理、風(fēng)速儀、車載供電系統(tǒng)及監(jiān)測(cè)車構(gòu)成.

    1.2 測(cè)量原理

    SOF測(cè)量示意見圖1.SOF儀器安裝在監(jiān)測(cè)車上,以太陽的紅外輻射為光源,利用太陽跟蹤器跟蹤太陽,并將陽光導(dǎo)入傅立葉變換紅外光譜儀,移動(dòng)測(cè)量排放煙羽的烴類(非甲烷烷烴和烯烴)紅外吸收,反演烴類柱濃度(mg/m2)分布,結(jié)合風(fēng)向、風(fēng)速測(cè)量和GPS定位測(cè)算烴類質(zhì)量排放速率(kg/h).在石化園區(qū)、廠區(qū)、裝置區(qū)、儲(chǔ)運(yùn)或廢水處理設(shè)施界外上下風(fēng)向或封閉監(jiān)測(cè)一圈,可獲得上述區(qū)域烴類排放速率(kg/h).在典型排放工況和典型氣象條件下測(cè)量多次,通過統(tǒng)計(jì)分析得到烴類代表性排放速率,根據(jù)排放污染帶中烴類占VOCs總量的比例(烴類總濃度/非甲烷有機(jī)物總濃度),可估算VOCs的排放速率,進(jìn)而推算烴類或VOCs年排放總量.

    SOF測(cè)量過程具體可分解為:烴類特征紅外吸收區(qū)太陽光紅外吸收光譜測(cè)量;實(shí)測(cè)紅外吸收光譜擬合反演烴類路徑積分濃度(柱濃度);移動(dòng)測(cè)量污染帶吸收截面,并結(jié)合風(fēng)速風(fēng)向測(cè)量計(jì)算烴類排放通量.

    太陽光穿過排放氣擴(kuò)散層會(huì)被多種待測(cè)物吸收,通過太陽跟蹤器跟蹤并將太陽光導(dǎo)入傅里葉變換紅外光譜儀(FTIR),測(cè)量太陽光中紅外吸收光譜.在烴類特征紅外吸收區(qū),用已知烴類的標(biāo)準(zhǔn)紅外光譜通過非線性擬合實(shí)測(cè)吸收光譜,得到烴類的柱濃度.烷烴C-H鍵伸縮振動(dòng)的特征紅外吸收譜帶在2700~3100cm-1,在低分辨(8cm-1)下測(cè)量,不同鏈長(zhǎng)的烷烴在該特征吸收譜帶的紅外吸收光譜很接近,長(zhǎng)鏈烷烴基本一致,僅乙烷和丙烷等短鏈烷烴有輕微差別[12].可以用少數(shù)幾種烷烴的紅外光譜擬合烷烴混合物的紅外吸收光譜,并估算烷烴混合物的質(zhì)量柱濃度[12].煉油廠排放中環(huán)烷烴的占比很低,SOF測(cè)量選擇乙烷、正丙烷、正丁烷、異戊烷和正辛烷用于擬合,測(cè)量區(qū)間2700~ 3005cm-1,采用標(biāo)準(zhǔn)吸收譜庫,非線性搜索算法擬合[12].水蒸汽和甲烷在此區(qū)域有吸收,參與擬合,但不定量,水蒸汽和甲烷采用HITRAN Database標(biāo)準(zhǔn)譜.甲烷在此區(qū)域的吸收主要為窄線,其大氣背景濃度較高,吸收線通過大氣層后幾乎耗盡,在分辨8cm-1下,甲烷響應(yīng)靈敏度很低,基本不影響其他烷烴的柱濃度定量.

    圖1 SOF測(cè)量示意

    SOF監(jiān)測(cè)排放通量的計(jì)算可通過微體積元積分實(shí)現(xiàn),原理如圖2所示.

    圖2 SOF監(jiān)測(cè)排放通量的計(jì)算

    C=m/(x×b×h) (1)

    b=u×t(2)

    SOF,i=C×h(3)

    m/t=SOF,i×x×u(4)

    =Sm/t=S(SOF,i×x×u) (5)

    式中:C為時(shí)刻微體積元內(nèi)待測(cè)物的濃度,mg/m3;m為時(shí)刻微體積元內(nèi)待測(cè)物的質(zhì)量;x為時(shí)刻微體積元SOF監(jiān)測(cè)車運(yùn)動(dòng)位移,通過行車過程GPS記錄測(cè)量;b為時(shí)刻微體積元風(fēng)位移,校正至與SOF監(jiān)測(cè)車運(yùn)動(dòng)方向垂直;h為微體積元太陽光測(cè)量光束穿過排放煙羽擴(kuò)散層路徑長(zhǎng)度,結(jié)合時(shí)刻陽光與垂直方向的角度校正到垂直地面方向;u為時(shí)刻平均風(fēng)速,通過風(fēng)速儀測(cè)量并記錄,結(jié)合時(shí)刻風(fēng)向與SOF監(jiān)測(cè)車移動(dòng)方向的夾角,校正至與SOF監(jiān)測(cè)車移動(dòng)方向垂直方向;t為時(shí)間;SOF,i為時(shí)刻待測(cè)物的柱濃度,即待測(cè)物的路徑(太陽光測(cè)量光束穿過排放氣擴(kuò)散層路徑校正到垂直方向)積分濃度(mg/m3×m=mg/m2),在監(jiān)測(cè)車移動(dòng)過程中用車載SOF監(jiān)測(cè)儀連續(xù)測(cè)量,并通過SOF測(cè)量和反演擬合得到;為待測(cè)物排放通量,可通過微體積元待測(cè)物柱濃度SOF,i、位移微元x和風(fēng)速u3者乘積累加得到.

    1.3 測(cè)量方法

    測(cè)量條件為:晴天有陽光,風(fēng)速2~12m/s,車速20~30km/h,SOF監(jiān)測(cè)車距離排放源10~5000m移動(dòng)測(cè)量,圍繞或在測(cè)量區(qū)域上下風(fēng)向測(cè)量.風(fēng)向和風(fēng)速采用機(jī)械式風(fēng)速儀(美國Young 公司)測(cè)量并實(shí)時(shí)記錄(美國Campbell Sci logger CR200X記錄儀),風(fēng)速儀安裝在測(cè)量區(qū)域內(nèi)或附近,安裝高度30~40m,周圍無障礙,避開地表紊流風(fēng).測(cè)量數(shù)據(jù)采用瑞典FluxSense公司的SOF-Measure和SOFsuite軟件采集和處理,SOF測(cè)量數(shù)據(jù)、監(jiān)測(cè)車GPS和風(fēng)速風(fēng)向通過時(shí)間同步測(cè)量和記錄,在衛(wèi)星遙感影像(帶有GPS定位信息)上生成測(cè)量區(qū)域非甲烷烷烴和烯烴排放柱濃度分布,結(jié)合風(fēng)速和風(fēng)向同步測(cè)量數(shù)據(jù)計(jì)算烴類排放速率.

    2 結(jié)果與討論

    2.1 煉油廠非甲烷烷烴排放測(cè)量結(jié)果

    這些煉油廠分布在東北1座(原油加工量約800萬t/a)、華北1座(原油加工量超過千萬t/a)、華東5座(原油加工量均超過千萬t/a).測(cè)量時(shí)間基本覆蓋四季,在春季、秋季和冬季各有2座煉油廠,橫跨夏季和秋季有1座煉油廠.每座煉油廠測(cè)量3~8d(SOF測(cè)量需要適宜的陽光和風(fēng)速,測(cè)量日期不一定連續(xù)),每天一般從上午09:00測(cè)量到下午15:00,每天的測(cè)量次數(shù)取決于廠區(qū)的大小和陽光情況.由表1可知,每座煉油廠累計(jì)測(cè)量次數(shù)為18~73次,具備一定的統(tǒng)計(jì)分析樣本;G煉油廠和T煉油廠(2021年監(jiān)測(cè))非甲烷烷烴排放速率相對(duì)標(biāo)準(zhǔn)偏差(RSD)較高,生產(chǎn)和排放工況波動(dòng)較大,其余6座煉油廠的非甲烷烷烴排放速率RSD為21.6%~33.1%,生產(chǎn)和排放工況相對(duì)穩(wěn)定.煉油廠VOCs排放污染主要通過風(fēng)場(chǎng)和擴(kuò)散遷移,由圖3可知非甲烷烷烴排放污染帶位置隨風(fēng)向變化,柱濃度隨風(fēng)速和源強(qiáng)變化.

    表1 7座大型煉油廠非甲烷烷烴排放SOF監(jiān)測(cè)結(jié)果

    注:*F煉油廠監(jiān)測(cè)范圍包括煉油區(qū)和乙烯生產(chǎn)區(qū),各占非甲烷烷烴排放總量的58.4%和41.6%.**監(jiān)測(cè)到1次事故性瞬間高排放(8640kg/h),用Grubbs檢驗(yàn)法剔除該非正常生產(chǎn)工況的極端排放值,采用其余20次監(jiān)測(cè)數(shù)據(jù)計(jì)算.***括號(hào)內(nèi)百分?jǐn)?shù)為相對(duì)標(biāo)準(zhǔn)偏差.

    圖3 煉油廠非甲烷烷烴排放SOF測(cè)量柱濃度分布

    2.2 煉油廠非甲烷烷烴排放波動(dòng)分析

    石化行業(yè)VOCs排放可歸類為十二個(gè)源項(xiàng)[2],其中設(shè)備密封點(diǎn)泄漏數(shù)量多且排放部位隨機(jī)不固定,有機(jī)液體儲(chǔ)存與調(diào)和揮發(fā)、有機(jī)液體裝卸揮發(fā)、工藝無組織排放及采樣過程排放數(shù)量較多、分散且間歇排放,廢水集輸、儲(chǔ)存、處理處置過程逸散和冷卻塔、循環(huán)水冷卻系統(tǒng)釋放通常無規(guī)則排放口,工藝有組織排放、燃燒煙氣排放基本為連續(xù)穩(wěn)態(tài)排放,火炬排放為點(diǎn)源非穩(wěn)態(tài)排放.上述排放源大部分為無組織排放,即使生產(chǎn)工況正常時(shí)總體排放也是非穩(wěn)態(tài)的,三座煉油廠測(cè)量日不同時(shí)段非甲烷烷烴排放測(cè)量結(jié)果反映排放波動(dòng)較大且無明顯規(guī)律性(圖4).6座千萬t級(jí)煉油企業(yè)2015年每個(gè)測(cè)量日各次SOF測(cè)量值的RSD分布見表2.由表2可知,各測(cè)量日的RSD幾乎呈現(xiàn)正態(tài)分布,RSD在20%~30%的測(cè)量日最多,70%測(cè)量日的RSD不超過30%.這6座千萬t級(jí)煉油企業(yè)2015年不同測(cè)量日非甲烷烷烴排放日均值匯總見表3,由表3可知,G煉油廠不同測(cè)量日的非甲烷烷烴排放日均值RSD高達(dá)68.2%,反映SOF監(jiān)測(cè)期間的生產(chǎn)和排放工況穩(wěn)定性較差;除G煉油廠外,其余6座煉油廠不同測(cè)量日的非甲烷烷烴排放日均值RSD在11.2%~24.9%,表明在生產(chǎn)工況平穩(wěn)的情況下,即使僅測(cè)量1d,測(cè)量日均值的偏差不超過25%.

    表2 每個(gè)測(cè)量日非甲烷烷烴排放相對(duì)標(biāo)準(zhǔn)偏差分布

    表3 不同測(cè)量日非甲烷烷烴排放日均值相對(duì)標(biāo)準(zhǔn)偏差分布

    注:“—”表示該日未測(cè)量.

    2.3 煉油廠非甲烷烷烴排放監(jiān)測(cè)代表性分析

    VOCs排放計(jì)算的可靠性由監(jiān)測(cè)的代表性決定,在典型氣象、生產(chǎn)和排放工況獲取一定SOF監(jiān)測(cè)數(shù)據(jù)樣本,通過統(tǒng)計(jì)分析可得到氣象條件(環(huán)境溫度、風(fēng)向、風(fēng)速、太陽輻射等)擾動(dòng)、生產(chǎn)和排放工況波動(dòng)影響下的代表性監(jiān)測(cè)數(shù)據(jù).平均值在某種程度上可代表煉油廠排放水平,但某些異常數(shù)據(jù)(如工況不穩(wěn)定)可能影響平均值的代表性.中位值是處在代表性范圍內(nèi)的一次監(jiān)測(cè)值,仍可能有一定的偏差.為從監(jiān)測(cè)數(shù)據(jù)組中獲取最具代表性的非甲烷烷烴排放數(shù)據(jù),分別繪制了7座煉油廠及其中6座千萬噸級(jí)煉油廠(2015年)合計(jì)的非甲烷烷烴排放SOF監(jiān)測(cè)值頻率分布柱形圖,并進(jìn)行對(duì)數(shù)擬合,結(jié)果見圖5,圖中橫軸為排放分布區(qū)間,縱軸為各區(qū)間內(nèi)監(jiān)測(cè)到次數(shù)的統(tǒng)計(jì)值,曲線為概率分布擬合線.

    圖5 煉油廠非甲烷烷烴排放速率頻率分布與擬合

    由圖5可知:(1)煉油廠非甲烷烷烴排放速率頻率符合正偏態(tài)分布,對(duì)數(shù)曲線適于擬合非甲烷烷烴SOF監(jiān)測(cè)頻率分布,6座千萬t級(jí)煉油廠合計(jì)(2015年)、F、Q、T、G、Z、J和Y煉油廠(2014~2015年)及T煉油廠(2021年)非甲烷烷烴排放速率SOF監(jiān)測(cè)值擬合相關(guān)系數(shù)()分別為0.9832、0.9317、0.9462、0.8040、0.9138、0.6102、0.8199、0.9318、0.7134,擬合優(yōu)度總體良好;(2)F煉油廠監(jiān)測(cè)到1個(gè)瞬時(shí)高排放(8640kg/h),可能為事故排放,但頻次低、時(shí)間短,對(duì)排放貢獻(xiàn)有限,如果用平均值計(jì)算可能高估排放,擬合可降低此類小概率事件的影響;(3)最大概率密度處的擬合值基本可反映監(jiān)測(cè)期間代表性非甲烷烷烴排放,該值更接近中位值,而不是平均值;(4) F(包括煉油和乙烯)、Q、T、G、Z、J、Y 煉油廠2014~2015年監(jiān)測(cè)非甲烷烷烴代表性排放速率分別為1690,880,1335,596,406,1148和1242kg/h,T煉油廠2021年監(jiān)測(cè)非甲烷烷烴代表性排放速率為368kg/h.

    2.4 煉油廠非甲烷烷烴及VOCs排放系數(shù)估算

    煉油廠非甲烷烷烴或VOCs排放系數(shù)定義為非甲烷烷烴或VOCs排放總量與其原油加工量的比值,可通過非甲烷烷烴或VOCs排放速率監(jiān)測(cè)值與監(jiān)測(cè)期間原油實(shí)際加工量的比值計(jì)算.基于7座煉油廠SOF監(jiān)測(cè)值統(tǒng)計(jì)擬合得到的非甲烷烷烴代表性排放速率計(jì)算了非甲烷烷烴排放系數(shù)(表4).煉油廠排放的VOCs主要為非甲烷烷烴、烯烴和芳烴,環(huán)烷烴及其他含硫、氮、氧、鹵素VOCs的含量較低.SOF可測(cè)量非甲烷烷烴的質(zhì)量排放速率,通過測(cè)量排放煙羽中非甲烷烷烴占VOCs總質(zhì)量濃度的比值可計(jì)算其質(zhì)量分?jǐn)?shù),可推算VOCs排放速率.煉油廠VOCs排放中的非甲烷烷烴質(zhì)量分?jǐn)?shù)與其加工工藝、原油種類、生產(chǎn)和儲(chǔ)運(yùn)設(shè)施管理等有關(guān),一般為70%~ 90%.國內(nèi)幾座煉油廠VOCs排放中非甲烷烷烴的質(zhì)量分?jǐn)?shù)分別為76.2%[20]、77.6%[21]、78.9%[22]、82%[18]和83.8%[23],平均約為79.7%.采用國內(nèi)煉油廠VOCs排放中非甲烷烷烴質(zhì)量分?jǐn)?shù)均值估算了7家煉油廠的VOCs排放速率和排放系數(shù)(表4).由表4可知,7座煉油廠2014~2015年非甲烷烷烴和VOCs排放系數(shù)均值分別為0.081%和0.10%,Z煉化公司的非甲烷烷烴和VOCs排放系數(shù)最低,分別為0.016%和0.020%,其他企業(yè)的非甲烷烷烴和VOCs排放系數(shù)約為Z煉油廠的3~7倍.7座煉油廠非甲烷烷烴和VOCs排放系數(shù)的RSD約為45%,反映各煉油廠非甲烷烷烴和VOCs排放控制水平差異較大.(3)7座煉油廠的區(qū)位、規(guī)模、工藝流程和原油類型比較典型,監(jiān)測(cè)時(shí)間基本覆蓋四季(有研究顯示煉油廠烷烴排放隨季節(jié)波動(dòng)有限[24]),監(jiān)測(cè)期間生產(chǎn)和排放工況基本正常,監(jiān)測(cè)的非甲烷烷烴或VOCs排放系數(shù)在一定程度上反映了當(dāng)時(shí)(2015年)煉油廠VOCs排放控制水平.(4)T煉油廠2021年監(jiān)測(cè)的非甲烷烷烴(或VOCs)排放速率和非甲烷烷烴(或VOCs)排放系數(shù)較2015年分別削減72.4%和74.2%.

    表4 基于SOF監(jiān)測(cè)估算的7座煉油廠非甲烷烷烴和VOCs排放系數(shù)

    注:* 2014年監(jiān)測(cè);** 2015年監(jiān)測(cè);*** 2021年監(jiān)測(cè).

    2.5 國內(nèi)外煉油廠非甲烷烷烴排放對(duì)比

    美國加州南岸空氣質(zhì)量管理局(SCAQMD)采用SOF技術(shù)于2015年監(jiān)測(cè)了轄區(qū)內(nèi)6座煉油廠非甲烷烷烴排放,監(jiān)測(cè)的非甲烷烷烴排放系數(shù)為0.016%~0.039%,平均為0.021%[25].國內(nèi)7座煉油廠與美國加州煉油廠同期SOF監(jiān)測(cè)結(jié)果相比,非甲烷烷烴排放系數(shù)的最好水平與其頂尖水平相同(均為0.016%),但平均水平(0.081%)約為其3.9倍,排放系數(shù)差異較大,范圍(0.016%~0.11%)約為其4.1倍.美國煉油廠VOCs排放控制起步較早,控制難度較大的VOCs無組織排放源均納入新源標(biāo)準(zhǔn)(NSPS)和有害空氣污染物排放標(biāo)準(zhǔn)(NESHAP),應(yīng)用了最佳實(shí)用控制技術(shù)(BAT),如1983年就已在煉油廠實(shí)施泄漏檢測(cè)與修復(fù)(LDAR)控制設(shè)備密封泄漏VOCs排放[26-27].美國SCAQMD的煉油廠排放控制法規(guī)和標(biāo)準(zhǔn)比聯(lián)邦更為嚴(yán)格,轄區(qū)內(nèi)煉油廠VOCs排放控制代表美國乃至世界領(lǐng)先水平.國內(nèi)VOCs排放控制起步較晚,2015年才開始在煉油廠實(shí)施LDAR等VOCs無組織排放控制措施,雖然2015年監(jiān)測(cè)的7座國內(nèi)煉油廠中有1座的VOCs排放控制達(dá)到國際先進(jìn)水平,但平均排放距離國際先進(jìn)水平尚有差距,各企業(yè)的VOCs排放控制水平參差不齊.T煉油廠VOCs排放在2015~2021年間降幅較大,2021年SOF監(jiān)測(cè)的非甲烷烷烴排放系數(shù)接近SCAQMD轄區(qū)煉油廠2015年平均水平.

    2.6 煉油廠非甲烷烷烴排放分布監(jiān)測(cè)及溯源

    圖6 SOF監(jiān)測(cè)煉油廠內(nèi)烷烴排放分布

    T煉油廠、Q煉油廠和J煉油廠廠內(nèi)非甲烷烷烴排放SOF監(jiān)測(cè)見圖6.監(jiān)測(cè)發(fā)現(xiàn),T煉油廠區(qū)西南角火炬存在未點(diǎn)火間歇排放,對(duì)全廠非甲烷烷烴排放的貢獻(xiàn)約為18.9%,汽油、航煤、輕污油等輕油罐區(qū)占全廠非甲烷烷烴排放的38.7%.Q煉油廠內(nèi)罐區(qū)合計(jì)約占全廠非甲烷烷烴排放的55.2%(其中汽油罐區(qū)約占80%);循環(huán)水冷卻塔約占全廠烷烴排放的10.1%,疑似換熱器泄漏通過循環(huán)水冷卻塔排放;生產(chǎn)裝置中,焦化和重整的非甲烷烷烴排放相對(duì)較高,分別占全廠非甲烷烷烴排放的7.5%(未監(jiān)測(cè)到焦炭塔開塔除焦時(shí)排放,焦化裝置的監(jiān)測(cè)值可能偏低)和6.4%;油品裝車、火炬和焦炭裝車分別約占全廠烷烴排放的6.2%、2.0%、1.9%;上述排放源中除火炬外均為無組織排放,約占該煉油廠非甲烷烷烴排放總量的87.3%,扣除循環(huán)水冷卻塔非正常工況排放(可能為換熱器泄漏事故性排放),該煉油廠正常工況下無組織排放約占非甲烷烷烴排放總量的77.2%.J煉油廠汽油調(diào)和罐和污油罐排放較高,裝置區(qū)中重整裝置排放相對(duì)較高,疑似存在換熱器泄漏通過循環(huán)水冷卻塔排放.

    2.7 煉油廠VOCs排放分析及控制措施建議

    無組織排放源是煉油廠非甲烷烷烴排放的主體,占比超過70%,其中輕油貯罐是最主要的非甲烷烷烴排放源,約占煉油廠非甲烷烷烴排放總量的40%以上.生產(chǎn)和排放工況穩(wěn)定是煉油廠VOCs無組織排放控制的重要保障,2014~2015年SOF監(jiān)測(cè)期間有1座煉油廠無異常排放,監(jiān)測(cè)的非甲烷烷烴排放系數(shù)最低,其余6座煉油廠中不同程度存在事故性瞬態(tài)高排放、火炬未點(diǎn)火或焚燒不完全間歇高排放、換熱器泄漏并通過循環(huán)水冷卻塔排放、浮頂罐泄漏等非正常排放,造成非甲烷烷烴排放系數(shù)較高且差異較大.歐盟研究發(fā)現(xiàn)定期用SOF等技術(shù)監(jiān)測(cè)煉油廠VOCs排放總量與分布,并針對(duì)性采取控制措施,可有效降低煉油廠VOCs非正常排放[28].

    國內(nèi)煉油廠排放清單核算所用VOCs排放系數(shù)為0.182%(1.82g/kg)[4,9],約為美國和歐洲煉油廠VOCs排放系數(shù)的10倍[7,29],約為本文SOF實(shí)測(cè)值的1.3~9.1倍(平均1.8倍),可能高估實(shí)際排放.煉油廠VOCs排放清單一般參考EPA排放系數(shù)/排放模型核算,美國EPA報(bào)告的煉油廠VOCs排放系數(shù)為0.01%~0.02%[7],歐洲多個(gè)煉油廠VOCs排放系數(shù)平均為0.0188%[29].排放通量監(jiān)測(cè)值與排放清單核算值通常不一致[30-33],包括SOF在內(nèi)的多種實(shí)際監(jiān)測(cè)顯示煉油廠VOCs排放通量測(cè)量值約為排放清單值的3~15倍[3,7].由于煉油廠VOCs無組織排放復(fù)雜且核算難度大、EPA排放系數(shù)/模型總體質(zhì)量不高、驗(yàn)證與修訂不充分、假定煉油廠運(yùn)轉(zhuǎn)在良好維護(hù)狀態(tài)、某些排放源缺乏成熟估算方法、計(jì)算參數(shù)選擇和方法解讀人為因素影響大等因素,排放清單核算結(jié)果不確定性較高,往往與體量而不是控制水平相關(guān),不能反映企業(yè)VOCs排放實(shí)際控制水平和企業(yè)間生產(chǎn)及VOCs治理管控差異,難以有效支撐控制決策和指導(dǎo)控制實(shí)踐.

    國內(nèi)7座煉油廠2014~2015年SOF監(jiān)測(cè)結(jié)果可作為石化VOCs無組織排放控制初期典型煉油企業(yè)非甲烷烷烴(或VOCs)基線排放.2015年以來,國家強(qiáng)化了VOCs排放治理與管控,有組織排放基本有效控制,煉油廠VOCs排放總體應(yīng)為下降趨勢(shì),T煉油廠2021年的非甲烷烷烴排放較2015年削減70%以上,治理措施主要包括泄漏檢測(cè)與修復(fù)(LDAR)、輕油或污油罐區(qū)排放控制、火炬非正常排放管控等.未來控制重點(diǎn)為無組織排放,應(yīng)重點(diǎn)解決同類企業(yè)VOCs排放嚴(yán)重參差不齊問題.SOF監(jiān)測(cè)為石化VOCs無組織排放監(jiān)測(cè)、量化核算、排放清單修訂控制決策、政策規(guī)范制定、措施優(yōu)化、減排效果評(píng)估、合規(guī)性監(jiān)督提供了BAT.煉油廠應(yīng)全面建立覆蓋排放源頭、過程、廠界及周邊社區(qū)的VOCs無組織排放的實(shí)時(shí)監(jiān)控、預(yù)警與溯源體系,保障企業(yè)運(yùn)行工況正常和平穩(wěn).每年采用SOF等VOCs排放通量監(jiān)測(cè)技術(shù)測(cè)量煉油廠非甲烷烷烴排放總量及分布,篩查重點(diǎn)排放源,評(píng)估控制效果,針對(duì)性地落實(shí)控制措施,驗(yàn)證和修訂排放系數(shù)/排放模型,建立高精度、高分辨VOCs排放清單管理體系.可采用VOCs排放實(shí)測(cè)值與清單值的比值衡量煉油廠VOCs排放實(shí)際控制效果,只要落實(shí)了VOCs排放控制BAT且運(yùn)維良好,該比值將在較低水平,該比值超過一定水平說明存在泄漏等異常排放或VOCs排放控制效果不佳.

    美國SCAQMD轄區(qū)煉油廠2015~2020年非甲烷烷烴SOF監(jiān)測(cè)值總體穩(wěn)定, VOCs排放管控治理已基本達(dá)到削減緩慢的平臺(tái)階段[34].參考其SOF監(jiān)測(cè)數(shù)據(jù)推測(cè),煉油廠VOCs排放控制最佳水平的非甲烷烷烴排放系數(shù)SOF監(jiān)測(cè)值為0.01%左右.

    2.8 SOF測(cè)量不確定性分析與改進(jìn)

    煉油廠非甲烷烷烴排放通量SOF測(cè)量不確定性主要源于非甲烷烷烴柱濃度、非甲烷烷烴背景柱濃度、基線、風(fēng)速、風(fēng)向和GPS(記錄SOF測(cè)量路線、速度、柱濃度的位置等)等測(cè)量的隨機(jī)或系統(tǒng)誤差.非甲烷烷烴柱濃度測(cè)量不確定性主要來自太陽紅外光譜測(cè)量和反演,非甲烷烷烴SOF測(cè)量光譜吸收截面和柱濃度反演的不確定度(相對(duì)誤差)分別約為3.5%和12%[12,35-36].所有煉油廠監(jiān)測(cè)期間均未發(fā)現(xiàn)測(cè)量區(qū)域上風(fēng)向非甲烷烷烴輸入背景,非甲烷烷烴背景柱濃度不確定度幾乎可以忽略不計(jì).基線不確定性源自測(cè)量基線波動(dòng)、漂移、噪聲和光譜干擾等,采用非甲烷烷烴柱濃度測(cè)量基線波動(dòng)或噪聲引起的柱濃度測(cè)量偏差,以及基線漂移引起的通量計(jì)算偏差估算了煉油廠排放SOF測(cè)量的基線不確定度,約為0.9%~10%,其中基線漂移貢獻(xiàn)較大.風(fēng)場(chǎng)主導(dǎo)的垂直風(fēng)梯度、排放煙羽高度及混合速度對(duì)排放通量計(jì)算不確定性貢獻(xiàn)較大,這些不確定性主要計(jì)入風(fēng)速.具體而言,風(fēng)速的不確定度來自3個(gè)方面:(1)風(fēng)速采用機(jī)械測(cè)風(fēng)儀在測(cè)量區(qū)域內(nèi)典型位置高度30~40m測(cè)量,與SOF測(cè)量時(shí)間同步但空間位置不同步,用于排放煙羽常規(guī)高度范圍內(nèi)SOF測(cè)量柱濃度垂向平均風(fēng)速存在一定不確定性;(2)排放煙羽高度可能超過常規(guī)高度范圍引起垂向平均風(fēng)速不確定性;(3)排放煙羽移動(dòng)速度與風(fēng)可能不完全同步引起通量測(cè)算不確定性(系統(tǒng)誤差);三者合計(jì)的不確定度約為11%~30%[12,35-36].排放煙羽運(yùn)動(dòng)方向與測(cè)量截面夾角的不確定性主要由風(fēng)向決定,風(fēng)向不確定度約為4%~12%[12,35-36].所有煉油廠SOF測(cè)量過程中GPS信號(hào)連續(xù)穩(wěn)定,GPS測(cè)量的不確定度約為0.1%[37],幾乎可以忽略不計(jì).上述所有測(cè)量不確定性的復(fù)合不確定度(各項(xiàng)相對(duì)誤差平方和開方)約為17%~36%.總體上風(fēng)場(chǎng)是非甲烷烷烴排放通量SOF測(cè)量不確定性的主要來源.作為標(biāo)準(zhǔn)化工作的一部分,CEN在法國一座停產(chǎn)的煉油廠開展了丙烷控制釋放盲測(cè),驗(yàn)證包括SOF在內(nèi)的多種監(jiān)測(cè)方法有效性,該測(cè)試場(chǎng)地及多點(diǎn)控制釋放較為接近VOCs實(shí)際排放場(chǎng)景,第三方操作控制釋放裝置,釋放量對(duì)監(jiān)測(cè)方保密[13],結(jié)果顯示大部分SOF測(cè)試相對(duì)偏差(測(cè)量值與實(shí)際釋放值)及全部測(cè)試(15次)平均相對(duì)偏差的絕對(duì)值均在上述估算的不確定度上限范圍內(nèi)[13].

    車載SOF可測(cè)量小到煉化廠區(qū),大到工業(yè)園區(qū)或區(qū)域排放(如大型石化工業(yè)帶或油氣勘探開發(fā)區(qū)域等)[12,38-40],船載[39]或機(jī)載[41-42]SOF測(cè)量平臺(tái)可實(shí)現(xiàn)港口或臨海工業(yè)園[39]、海上或山火[41]等更復(fù)雜區(qū)域排放測(cè)量.為更為準(zhǔn)確地測(cè)量區(qū)域排放,可從SOF測(cè)量技術(shù)和地面/機(jī)載/星載測(cè)量與SOF測(cè)量交互校驗(yàn)/校正兩方面改進(jìn)SOF測(cè)量的不確定性. SOF測(cè)量技術(shù)方面,需改進(jìn)對(duì)測(cè)量不確定性影響最大的風(fēng)場(chǎng)測(cè)量技術(shù),研發(fā)與SOF監(jiān)測(cè)同步的車載測(cè)風(fēng)雷達(dá);同時(shí)可采用無人機(jī)載技術(shù)測(cè)量排放煙羽廓線及非甲烷烷烴混合比,提高非甲烷烷烴混合比的準(zhǔn)確性和代表性;拓展SOF測(cè)量組分范圍,測(cè)量主體組分并增加SOF測(cè)量的組分?jǐn)?shù)量,以降低VOCs排放測(cè)算的不確定性. SOF一般測(cè)量主體組分排放通量,并通過排放煙羽中主體組分在VOCs總質(zhì)量濃度的占比測(cè)算VOCs排放通量.煉油廠VOCs排放主體組分為非甲烷烷烴,其他類型工業(yè)園區(qū)或區(qū)域排放的VOCs主體組分可能不同.除非甲烷烷烴外,SOF還能測(cè)量烯烴[38]、氯乙烯[43]、氨[44]和一氧化碳[41-42]等有特征紅外吸收組分的排放通量.交互校驗(yàn)/校正兩方面,用機(jī)載采樣分析或地面大氣自動(dòng)站監(jiān)測(cè)到的排放煙羽VOCs組成數(shù)據(jù)校驗(yàn)SOF質(zhì)量柱濃度測(cè)量反演的準(zhǔn)確性[12],區(qū)域排放混合帶機(jī)載測(cè)量可與SOF測(cè)量相互驗(yàn)證排放通量的準(zhǔn)確性,美國休斯頓附近大型石化區(qū)烯烴排放SOF與機(jī)載測(cè)量數(shù)據(jù)基本一致[38,45].機(jī)載SOF測(cè)量已用于校正或改進(jìn)星載遙感監(jiān)測(cè)的不確定性[41].

    3 結(jié)論

    3.1 采用SOF技術(shù)在2014.05~2015.12和2021.10系統(tǒng)監(jiān)測(cè)了國內(nèi)7座大型煉油廠非甲烷烷烴排放通量及分布,監(jiān)測(cè)結(jié)果反映了國內(nèi)石化VOCs無組織排放控制初期的基線排放水平.監(jiān)測(cè)發(fā)現(xiàn)煉油廠非甲烷烷烴排放主要來自無組織排放源,其中輕油貯罐約占煉油廠非甲烷烷烴排放總量的40%以上.在生產(chǎn)和排放工況相對(duì)平穩(wěn)的情況下,煉油廠非甲烷烷烴排放波動(dòng)有限,總體呈偏正態(tài)分布,日均值的相對(duì)偏差不超過25%,非甲烷烷烴代表性排放更接近中位值而不是平均值.

    3.2 7座煉油廠2014~2015年非甲烷烷烴排放系數(shù)測(cè)定值為0.016%~0.11%,平均值為0.081%;VOC排放系數(shù)估算值為0.020%~0.14%,平均值為0.10%.與美國加州6座煉油廠同期SOF監(jiān)測(cè)相比,國內(nèi)煉油廠非甲烷烷烴排放系數(shù)的最好水平與其頂尖水平相當(dāng),但平均水平約為其平均水平的3.9倍,且排放系數(shù)的差異度較高,總體與國際先進(jìn)水平尚有差距.

    3.3 經(jīng)過6a的管控治理后,本文7座煉油廠中的1座千萬t級(jí)煉油廠2021年非甲烷烷烴排放系數(shù)較2015年削減74.2%.

    3.4 SOF技術(shù)可用于石化企業(yè)或園區(qū)VOCs無組織排放監(jiān)測(cè)、監(jiān)管、管控與治理評(píng)估、排放清單實(shí)測(cè)、校驗(yàn)與修訂,未來SOF技術(shù)將向規(guī)范化、多組分及高精度測(cè)量發(fā)展.

    [1] Li M, Zhang Q, Zheng B, et al. Persistent growth of anthropogenic non-methane volatile organic compound (NMVOC) emissions in China during 1990~2017: drivers, speciation and ozone formation potential [J]. Atmos Chem Phys, 2019,19:8897-8913.

    [2] 環(huán)境保護(hù)部辦公廳.環(huán)辦[2015]104號(hào)關(guān)于印發(fā)《石化行業(yè)VOCs污染源排查工作指南》及《石化企業(yè)泄漏檢測(cè)與修復(fù)工作指南》的通知[EB/OL].https://www.mee.gov.cn/gkml/hbb/bgt/201511/ t20151124_317577.htm. 2015-11-17.

    Office of the Ministry of Environmental Protection. Huan Ban [2015] No. 104Notice on issuing《Guidelines for investigation of VOCs Pollution Sources in petrochemical industry 》and 《Guidelines for leak detection and repair in petrochemical enterprises》[EB/OL].https://www.mee.gov.cn/gkml/hbb/bgt/201511/t20151124_317577.htm. 2015-11-17.

    [3] 李凌波,程夢(mèng)婷,李 龍,等.煉油企業(yè)揮發(fā)性有機(jī)物無組織排放通量監(jiān)測(cè)的現(xiàn)狀與發(fā)展[J]. 中國環(huán)境監(jiān)測(cè), 2020,36(3):19-28.

    Li L, Cheng M, Li L, et al. Current and future developments in fugitive volatile organic compounds emission flux monitoring in petroleum refining industry [J]. Environmental Monitoring in China, 2020,36(3):19-28.

    [4] Simayi M, Hao Y, Li J, et al. Historical volatile organic compounds emission performance and reduction potentials in China’s petroleum refining industry [J]. Journal of Cleaner Production, 2021,292:125810.

    [5] U.S.Environmental Protection Agency.AP-42:Compilation of Air Emissions Factors [EB/OL]. https://www.epa.gov/air-emissions- factors-and-quantification/ap-42-compilation-air-emissions-factors.

    [6] Nelson T P. An examination of historical air pollutant emissions from US petroleum refineries [J]. Environmental Progress & Sustainable Energy, 2013,32(2):425-432.

    [7] Cuclis A. Why emission factors don’t work at refineries and what to do about it [EB/OL]. https://www3.epa.gov/ttnchie1/conference/ei20/ session7/acuclis.pdf. 2012.

    [8] Karl T, Striednig M, Graus M, et al. Urban flux measurements reveal a large pool of oxygenated volatile organic compound emissions [J]. Proc Natl Acad Sci USA, 2018,115(6):1186-1191.

    [9] 環(huán)境保護(hù)部.公告2014年第55號(hào)關(guān)于發(fā)布《大氣細(xì)顆粒物一次源排放清單編制技術(shù)指南(試行)》等4項(xiàng)技術(shù)指南的公告附件2大氣揮發(fā)性有機(jī)物源排放清單編制技術(shù)指南(試行) [EB/OL].https: //www.mee.gov.cn/gkml/hbb/bgg/201408/t20140828_288364.htm. 2014-08-28.

    Ministry of Environmental Protection. Announcement [2014] No. 55 on Issuing four technical guidelines including《Technical guidelines for the preparation of emission inventories of atmospheric fine particulate matter from primary sources (trial implementation)》,Annex 2 《Technical guidelines for the preparation of emission inventories of atmospheric volatile organic compounds (trial implementation)》[EB/OL].https://www.mee.gov.cn/gkml/hbb/bgg/201408/t20140828_ 288364.htm. 2014-08-28.

    [10] U.S. Environmental Protection Agency Office of Air Quality Planning and Standards Air Quality Assessment Division Measurement Technology Group. EPA Handbook: Optical and remote sensing for measurement and monitoring of emissions flux of gases and particulate matter [EB/OL]. https://www.epa.gov/sites/production/ files/2018-08/documents/gd-52v.2.pdf.

    [11] Robinson R, Gardiner T, Innocenti F, et al. Infrared differential absorption Lidar (DIAL) measurements of hydrocarbon emissions [J]. Journal of Environmental Monitoring, 2011,13:2213-2220.

    [12] Johansson J K E, Mellqvist J, Samuelsson J, et al. Emission measurements of alkenes, alkanes, SO2, and NO2from stationary sources in Southeast Texas over a 5year period using SOF and mobile DOAS [J]. J Geophys Res Atmos, 2014,119:1973-1991.

    [13] Robinson R, Innocenti F, Helmore J, et al. The development and validation of a new standardized method for monitoring diffuse and fugitive VOC emissions from oil and gas facilities [C]. AWMA 2019 Air Quality Measurement Methods and Technology, April 2-4, 2019.Durham,NC.

    [14] European Committee for Standardization (CEN). EN 17628:2022 Fugitive and diffuse emissions of common concern to industry sectors-Standard method to determine diffuse emissions of volatile organic compounds into the atmosphere [S]. April 2022.

    [15] Wang F, Du W, Lv S, et al. Spatial and temporal distributions and sources of anthropogenic NMVOCs in the atmosphere of China: a review [J]. Advances in Atmosphere Science, 2021,38(7):1085-1100.

    [16] Wei W, Cheng S, Li G, et al. Characteristics of ozone and ozone precursors (VOCs and NO) around a petroleum refinery in Beijing, China [J]. Journal of Environmental Sciences, 2014,26:332-342.

    [17] Lv D, Lu S, Tan X, et al. Source profiles, emission factors and associated contributions to secondary pollution of volatile organic compounds (VOCs) emitted from a local petroleum refinery in Shandong [J]. Environmental Pollution, 2021,274:116589.

    [18] Zhang Z, Wang H, Chen D, et al. Emission characteristics of volatile organic compounds and their secondary organic aerosol formation potentials from a petroleum refinery in Pearl River Delta, China [J]. Science of the Total Environment, 2017,584-585:1162-1174.

    [19] Lv D, Lu S, He S, et al. Research on accounting and detection of volatile organic compounds from a typical petroleum refinery in Hebei, North China [J]. Chemosphere, 2021,281:130653.

    [20] Zhang G, Wang N, Jiang X, et al. Characterization of ambient volatile organic compounds (VOCs) in the area adjacent to a petroleum refinery in Jinan, China [J]. Aerosol and Air Quality Research, 2017, 17:944-950.

    [21] 李勤勤,張志娟,李 楊,等.石油煉化無組織VOCs的排放特征及臭氧生成潛力分析[J]. 中國環(huán)境科學(xué), 2016,36(5):1323-1331.

    Li Q, Zhang Z, Li Y, et al. Characteristics and ozone formation potential of fugitive volatile organic compounds(VOCs) emitted from petrochemical industry in Pearl River Delta [J]. China Environmental Science, 2016,36(5):1323-1331.

    [22] 呂兆豐,魏巍,楊 干,等.某石油煉制企業(yè)VOCs排放源強(qiáng)反演研究[J]. 中國環(huán)境科學(xué), 2015,35(10):2958-2963.

    Lv Z, Wei W, Yang G, et al. Inversion research in VOCs source emission of a petroleum refinery [J]. China Environmental Science, 2015,35(10):2958-2963.

    [23] 孫筱強(qiáng),程水源,魏 巍,等.煉油廠VOCs排放及其對(duì)O3生成的敏感性研究[J]. 安全與環(huán)境學(xué)報(bào), 2011,11(6):101-105.

    Sun X, Cheng S, Wei W, et al. Research on the refinery O3variation characteristics and VOCs reactivity [J]. Journal of Safety and Environment, 2011,11(6):101-105.

    [24] Pikelnaya O, Tsai C, Polidori A, et al. Understanding seasonal variations in refinery voc emissions using quarterly optical remote sensing mobile surveys [C]//AWMA 2019 Air Quality Measurement Methods and Technology. April 2-4, 2019.Durham,NC.

    [25] Mellqvist J, Samuelsson J, Isoz O, et al. Emission measurements of VOCs, NO2and SO2from the refineries in the south coast air basin using solar occultation flux and other optical remote sensing methods [EB/OL]. http://www.aqmd.gov/fenceline-monitoring/project-1.

    [26] 李凌波,劉忠生,方向晨,等.煉油廠VOC排放控制策略-設(shè)備與管閥件泄漏[J]. 當(dāng)代石油石化, 2013,21(9):1-9.

    Li L, Liu Z, Fang X, et al. Strategies of refinery VOC emission control-equipment leaks [J]. Petroleum & Petrochemical Today, 2013,21(9):1-9.

    [27] 李凌波,劉忠生,方向晨.煉油廠VOC排放控制策略-儲(chǔ)運(yùn)、廢水處理、工藝尾氣、冷卻塔及火炬[J]. 當(dāng)代石油石化, 2013,(10):4-12.

    Li L, Liu Z, Fang X. The strategies for refinery VOC emission control-storage tanks and transfer operations, wastewater treatment, process vents, cooling towers and flares [J]. Petroleum & Petrochemical Today, 2013,21(10):4-12.

    [28] Barthe P, Chaugny M, Roudier S, et al. Best available techniques (BAT) reference document for the refining of mineral oil and gas. Industrial Emissions Directive 2010/75/EU (Integrated Pollution Prevention and control) [EB/OL]. 2015. http://eippcb.jrc.ec.europa.eu/reference/ BREF/REF_BREF_2015.pdf.

    [29] Roveda L, Polvara E, Invernizzi M, et al. Definition of an emission factor for VOC emitted from Italian and European refineries [J]. Atmosphere, 2020,11(6):564-576.

    [30] Acton W J F, Huang Z, Davison B, et al. Surface–atmosphere fluxes of volatile organic compounds in Beijing [J]. Atmos. Chem. Phys., 2020,20: 15101-15125.

    [31] Ryerson T B, Trainer M, Angevine W M, et al. Effect of petrochemical industrial emissions of reactive alkenes and NOon tropospheric ozone formation in Houston, Texas [J]. J. Geophys. Res., 2003, 108(D8):4249.

    [32] Henry R C, Spiegelman C H, Collins J F, et al. Reported emissions of organic gases are not consistent with observations [J]. Proc Natl Acad Sci USA, 1997,94:6596-6599.

    [33] Li S M, Leithead A, Moussa S G, et al. Differences between measured and reported volatile organic compound emissions from oil sands facilities in Alberta, Canada [J]. Proc. Natl. Acad. Sci. USA, 2017,114 (19):E3756-E3765.

    [34] Porter J, Tsai C, Pikelnaya O, et al. Routine direct measurements of refinery emissions of total alkanes over the last 5years in the Los Angeles air basin [C]//AGU Fall Meeting 2020.December 1-17, 2020.

    [35] Johansson J K E, Mellqvist J, Samuelsson J, et al. Quantitative measurements and modeling of industrial formaldehyde emissions in the Greater Houston area during campaigns in 2009 and 2011 [J]. J. Geophys. Res. Atmos., 2014,119:4303-4322.

    [36] European Committee for Standardization (CEN).CEN/TC 264/WG 38 "Determination of diffuse VOC emissions" Final Report [EB/OL].https://www.vdi.de/fileadmin/pages/vdi_de/redakteure/ueber_uns/fachgesellschaften/KRdL/dateien/WG_38_Final_Report_SACEN2014-07.pdf.

    [37] Kille N. Development of solar occultation flux spectroscopy for ground and airborne applications [D]. Colorado: University of Colorado, 2020.

    [38] Mellqvist J, Samuelsson J, Johansson J, et al. Measurements of industrial emissions of alkenes in Texas using the solar occultation flux method [J]. J. Geophys. Res., 2010,115:D00F17.

    [39] Mellqvist J, Samuelsson J, Andersson P, et al. Using solar occultation flux and other optical remote sensing methods to measure VOC emissions from a variety of stationary sources in the South Coast Air Basin [EB/OL]. http://www.aqmd.gov/docs/default-source/fenceline_ monitroing/project_2/fluxsense_project2_2015_final_report.pdf.

    [40] Mellqvist J, Samuelsson J, Offerle B, et al. VOC emissions monitoring in the Lost Hills Vicinity, San Joaquin Valley [EB/OL]. https://ww2. arb.ca.gov/sites/default/files/2021-03/FluxSenseInc_CARB_18ISD023_SJV_final_update_201027.pdf.

    [41] Bela M M, Kille N, McKeen S A, et al. Quantifying carbon monoxide emissions on the scale of large wildfires [J]. Geophysical Research Letters, 2022,49:e2021GL095831.

    [42] Kille N, Zarzana K J, Alvarez J R, et al. The CU airborne solar occultation flux instrument: performance evaluation during BB-FLUX [J]. ACS Earth Space Chem., 2022,6:582-596.

    [43] Johansson J, Mellqvist J, Li L, et al.Quantification of industrial VOC emissions in China using solar occultation flux (SOF) [EB/OL]. https://meetingorganizer.copernicus.org/EGU2017/EGU2017-19062.pdf.

    [44] KilleN, BaidarS, HandleyP, et al. The CU mobile solar occultation flux instrument: structure functions and emission rates of NH3, NO2and C2H6[J].Atmos. Meas. Tech., 2017,10:373-392.

    [45] Degouw J A, Hekkert S T, Mellqvist J, et al. Airborne measurements of ethene from industrial sources using laser photo-acoustic spectroscopy [J]. Environ. Sci. Technol., 2009,43:2437-2442.

    Measurement of emission fluxes of total non-methane alkanes from refineries using solar occultation flux remote sensing technique.

    LILing-bo1*, GONG Chao2, CHENG Meng-ting1, LI Long1, LIU Xin-yu1

    (1.SINOPEC Dalian Research Institute of Petroleum and Petrochemicals, Dalian 116045, China;2.SINOPEC Anqing Petrochemical Company, Anqing 246002, China)., 2022,42(7):3046~3057

    The emission fluxes (kg/h) and distributions of total non-methane alkanes from seven major refineries (in which six refineries were processing crude oil more than ten million tons per year) in China were quantified during May 2014 and December 2015 and in October 2021 using the Solar Occultation Flux (SOF) method. At each refinery site, 18~73measurement transects were performed during three to eight individual measurement days. In total, 328measurements of total non-methane alkane emission flux were performed for the seven refineries. Total VOC emissions were estimated from the measured total non-methane alkane emissions based on the average mass fraction of non-methane alkanes in total volatile organic compounds (VOC) measured in VOC emission plumes of several refineries in China. Total non-methane alkanes and VOC emission factors were calculated by scaling measured emissions of total non-methane alkanes and estimated emissions of VOCs with crude oil processing rates during the measurement period. The monitoring results indicate that the measured total non-methane alkanes emission factors ranged from 0.016% to 0.11% for the seven studied refineries and were 0.081% on average. The estimated VOC emission factors for the seven studied refineries ranged from 0.020% to 0.14%, with an average value of 0.10% during 2014 and 2015. Fugitive emissions accounted for more than 70% of total non-methane alkane emissions from refineries, and light oil tanks contributed more than 50% of fugitive non-methane alkane emissions. The results were compared to SOF measurements at six refineries in the South Coast Air Quality Management District (SCAQMD) in 2015, and the lowest level of total non-methane alkane emission factors in the seven refineries in China measured during 2014 and 2015 was the same as that for all six refineries in SCAQMD. However, the average level of total non-methane alkane emission factors of the seven refineries in China was 3.9times higher than that of the six refineries in SCAQMD, and there was a higher variability in total non-methane alkane emission factors among the seven refineries in China. The emission flux and emission factor of non-methane alkanes monitored in a major Chinese refinery in 2021 showed a reduction by 72.4% and 74.2% respectively compared with 2015. SOF could be developed as the best available technique for monitoring and quantifying fugitive and diffuse VOC emissions from petroleum and petrochemical industrial sites, and for revising VOC emission inventories. This study has provided a baseline of measured emissions of total non-methane alkanes and VOCs from typical refineries in the initial stage of VOC emissions control actions implemented in China in 2015, and VOCs control effectiveness of a major refinery after 6 years.

    refinery;non-methane alkanes;volatile organic compounds;diffuse and fugitive emission;emission flux measurement;solar occultation flux (SOF)

    X831;X87

    A

    1000-6923(2022)07-3046-12

    李凌波(1969-),男,內(nèi)蒙莫旗人,教授級(jí)高級(jí)工程師,工學(xué)學(xué)士,主要從事石油石化環(huán)境監(jiān)測(cè)及污染物組學(xué)研究.發(fā)表論文40余篇.

    2021-12-23

    中國石油化工股份有限公司科技開發(fā)項(xiàng)目(CLY15070,319006-6)

    * 責(zé)任作者, 教授級(jí)高級(jí)工程師, lilingbo.fshy@sinopec.com

    猜你喜歡
    煉油廠烷烴通量
    冬小麥田N2O通量研究
    氣相色譜六通閥在正構(gòu)烷烴及碳數(shù)分布測(cè)定中的應(yīng)用
    云南化工(2021年11期)2022-01-12 06:06:30
    高苯原料油烷烴異構(gòu)化的MAX-ISOM技術(shù)
    ExxonMobil公司宣布將在2025年底投資更多的基礎(chǔ)油項(xiàng)目
    自動(dòng)點(diǎn)火系統(tǒng)FCS在煉油廠火炬系統(tǒng)中的應(yīng)用
    煉油廠節(jié)能減排儲(chǔ)氣柜設(shè)計(jì)及安裝
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    烷烴油滴在超臨界二氧化碳中溶解的分子動(dòng)力學(xué)模擬
    中國石油工程建設(shè)公司承建的蘇丹煉油廠
    春、夏季長(zhǎng)江口及鄰近海域溶解甲烷的分布與釋放通量
    19禁男女啪啪无遮挡网站| 桃花免费在线播放| 在线永久观看黄色视频| 精品人妻一区二区三区麻豆| 国产精品一区二区精品视频观看| 精品国产国语对白av| 日本av手机在线免费观看| 欧美黄色淫秽网站| 欧美日韩福利视频一区二区| 国产精品偷伦视频观看了| 日韩视频一区二区在线观看| 十八禁网站网址无遮挡| 日韩欧美一区二区三区在线观看 | 人人妻人人爽人人添夜夜欢视频| 欧美在线黄色| 97在线人人人人妻| 国产日韩欧美亚洲二区| videos熟女内射| 亚洲专区中文字幕在线| 不卡一级毛片| 超碰97精品在线观看| 性色av一级| 又紧又爽又黄一区二区| 国产精品免费视频内射| 亚洲欧美日韩另类电影网站| 国产精品久久久久久精品电影小说| 一级毛片女人18水好多| 亚洲 欧美一区二区三区| 黄色a级毛片大全视频| 91大片在线观看| 丰满迷人的少妇在线观看| 美女扒开内裤让男人捅视频| 亚洲熟女精品中文字幕| 大陆偷拍与自拍| 国产成人精品在线电影| 男女床上黄色一级片免费看| 日本五十路高清| 秋霞在线观看毛片| 亚洲专区字幕在线| 国产精品久久久久成人av| 国产亚洲精品一区二区www | 香蕉丝袜av| 国产在线视频一区二区| 国产伦理片在线播放av一区| 亚洲精品久久久久久婷婷小说| av超薄肉色丝袜交足视频| 久久天躁狠狠躁夜夜2o2o| 久9热在线精品视频| 久久久久国产一级毛片高清牌| www.熟女人妻精品国产| 人妻久久中文字幕网| 国产精品免费大片| 亚洲欧洲日产国产| 国产日韩一区二区三区精品不卡| 欧美大码av| 淫妇啪啪啪对白视频 | 国产精品99久久99久久久不卡| 一边摸一边抽搐一进一出视频| 女性被躁到高潮视频| 一二三四在线观看免费中文在| 老司机在亚洲福利影院| 高清欧美精品videossex| 热99国产精品久久久久久7| 咕卡用的链子| 高清欧美精品videossex| 国产精品久久久久久精品电影小说| 欧美国产精品va在线观看不卡| 免费日韩欧美在线观看| 下体分泌物呈黄色| 少妇 在线观看| 国产亚洲av高清不卡| 亚洲专区国产一区二区| 99热全是精品| 久久人人爽人人片av| 激情视频va一区二区三区| 一二三四在线观看免费中文在| 久久人人爽人人片av| 婷婷成人精品国产| 999久久久精品免费观看国产| 亚洲成人免费电影在线观看| 亚洲精品成人av观看孕妇| 麻豆乱淫一区二区| 女人被躁到高潮嗷嗷叫费观| 女人被躁到高潮嗷嗷叫费观| 国产高清视频在线播放一区 | 又大又爽又粗| 国产色视频综合| 久久久精品国产亚洲av高清涩受| 久久久久久亚洲精品国产蜜桃av| 日韩中文字幕视频在线看片| www.精华液| www.熟女人妻精品国产| 久久影院123| 免费在线观看日本一区| 老司机深夜福利视频在线观看 | 日韩欧美一区二区三区在线观看 | 中亚洲国语对白在线视频| 交换朋友夫妻互换小说| 色精品久久人妻99蜜桃| 老鸭窝网址在线观看| 日韩大码丰满熟妇| 如日韩欧美国产精品一区二区三区| 成年人黄色毛片网站| 国产熟女午夜一区二区三区| 操出白浆在线播放| 十分钟在线观看高清视频www| 国产97色在线日韩免费| av网站在线播放免费| 丰满迷人的少妇在线观看| 色视频在线一区二区三区| 国产av国产精品国产| 9191精品国产免费久久| 亚洲成人手机| 黄色毛片三级朝国网站| 国产99久久九九免费精品| 窝窝影院91人妻| 一本大道久久a久久精品| 人妻人人澡人人爽人人| 亚洲国产精品成人久久小说| 下体分泌物呈黄色| 国产日韩欧美视频二区| 黄色 视频免费看| 国产精品一区二区精品视频观看| 精品久久久精品久久久| 丝袜脚勾引网站| 黑人巨大精品欧美一区二区mp4| 久久久精品区二区三区| 亚洲熟女毛片儿| 丝瓜视频免费看黄片| 在线观看人妻少妇| 一区福利在线观看| 欧美激情 高清一区二区三区| 日本黄色日本黄色录像| 久久久国产精品麻豆| 天堂中文最新版在线下载| 中文欧美无线码| 夜夜夜夜夜久久久久| 丰满人妻熟妇乱又伦精品不卡| 青春草亚洲视频在线观看| 久久 成人 亚洲| 国产成+人综合+亚洲专区| 秋霞在线观看毛片| 亚洲国产欧美日韩在线播放| 国产1区2区3区精品| netflix在线观看网站| 欧美精品啪啪一区二区三区 | 欧美激情极品国产一区二区三区| 欧美在线一区亚洲| 国产精品亚洲av一区麻豆| 久久99热这里只频精品6学生| 两性午夜刺激爽爽歪歪视频在线观看 | 国产欧美日韩精品亚洲av| 黄片大片在线免费观看| 久久久久久久久久久久大奶| 久9热在线精品视频| 少妇的丰满在线观看| 亚洲精品国产一区二区精华液| 久久ye,这里只有精品| 久久国产亚洲av麻豆专区| a级毛片黄视频| 蜜桃国产av成人99| 日韩欧美一区视频在线观看| 国产精品自产拍在线观看55亚洲 | 久久人人爽av亚洲精品天堂| 亚洲国产欧美网| 777久久人妻少妇嫩草av网站| 在线亚洲精品国产二区图片欧美| 搡老岳熟女国产| 国产在线免费精品| 老鸭窝网址在线观看| 国产精品秋霞免费鲁丝片| 高清av免费在线| 午夜91福利影院| 国产精品一二三区在线看| 97在线人人人人妻| 国产精品国产av在线观看| 日本精品一区二区三区蜜桃| 狠狠婷婷综合久久久久久88av| 久久九九热精品免费| 人妻久久中文字幕网| 在线观看舔阴道视频| 久久人人爽人人片av| videos熟女内射| 99国产精品免费福利视频| cao死你这个sao货| 亚洲人成电影免费在线| 国产男女内射视频| 国产成人av激情在线播放| 在线观看舔阴道视频| 在线十欧美十亚洲十日本专区| 精品亚洲乱码少妇综合久久| videos熟女内射| 在线天堂中文资源库| 久久精品国产亚洲av高清一级| 亚洲人成电影免费在线| 久久久精品94久久精品| 日日摸夜夜添夜夜添小说| 久热爱精品视频在线9| 成人国产一区最新在线观看| 精品人妻在线不人妻| 成在线人永久免费视频| 波多野结衣一区麻豆| 欧美精品av麻豆av| 国产成人a∨麻豆精品| 国产日韩欧美视频二区| 日本五十路高清| 亚洲成av片中文字幕在线观看| 国产精品久久久久久人妻精品电影 | 午夜福利乱码中文字幕| 性高湖久久久久久久久免费观看| 国产黄色免费在线视频| 中文字幕色久视频| 久久 成人 亚洲| 亚洲国产欧美日韩在线播放| 午夜福利乱码中文字幕| 免费观看av网站的网址| 午夜老司机福利片| 日本av手机在线免费观看| 成年人午夜在线观看视频| 久久精品国产a三级三级三级| 国产一区二区三区综合在线观看| 男女免费视频国产| 免费不卡黄色视频| av在线播放精品| 久久人人爽人人片av| 一区二区日韩欧美中文字幕| 亚洲欧美精品自产自拍| 国产精品亚洲av一区麻豆| 久久久久网色| 久久久精品免费免费高清| 免费看十八禁软件| 免费av中文字幕在线| 在线十欧美十亚洲十日本专区| 少妇精品久久久久久久| 美女脱内裤让男人舔精品视频| 波多野结衣一区麻豆| 精品一区二区三区av网在线观看 | av在线播放精品| 午夜两性在线视频| h视频一区二区三区| 欧美精品人与动牲交sv欧美| 日韩中文字幕视频在线看片| 国产精品久久久久成人av| 国产日韩欧美亚洲二区| 久久女婷五月综合色啪小说| 午夜精品久久久久久毛片777| 亚洲情色 制服丝袜| 麻豆乱淫一区二区| 叶爱在线成人免费视频播放| 亚洲五月婷婷丁香| 丰满迷人的少妇在线观看| 一二三四在线观看免费中文在| 男女无遮挡免费网站观看| 亚洲精品中文字幕一二三四区 | 国产精品秋霞免费鲁丝片| 午夜福利视频精品| 999精品在线视频| 久久精品国产亚洲av香蕉五月 | 一级黄色大片毛片| 欧美97在线视频| 亚洲黑人精品在线| av天堂久久9| 又紧又爽又黄一区二区| 波多野结衣av一区二区av| 亚洲av电影在线进入| 欧美日韩视频精品一区| 日日爽夜夜爽网站| 欧美日韩亚洲高清精品| 日韩视频一区二区在线观看| 亚洲欧美成人综合另类久久久| 啦啦啦视频在线资源免费观看| 久久久精品区二区三区| 亚洲精品美女久久av网站| 麻豆国产av国片精品| 精品国产乱子伦一区二区三区 | 国产成人影院久久av| 色视频在线一区二区三区| 国产精品久久久久久精品古装| 日本精品一区二区三区蜜桃| 香蕉国产在线看| 人人澡人人妻人| 免费在线观看影片大全网站| 高清视频免费观看一区二区| 美女主播在线视频| 脱女人内裤的视频| 好男人电影高清在线观看| 99精国产麻豆久久婷婷| 亚洲免费av在线视频| av电影中文网址| 国产亚洲精品第一综合不卡| 精品国产一区二区三区久久久樱花| 桃花免费在线播放| 欧美国产精品va在线观看不卡| 极品少妇高潮喷水抽搐| 欧美xxⅹ黑人| 婷婷丁香在线五月| 97人妻天天添夜夜摸| 欧美老熟妇乱子伦牲交| 色精品久久人妻99蜜桃| 两个人看的免费小视频| 飞空精品影院首页| 大陆偷拍与自拍| 可以免费在线观看a视频的电影网站| 免费日韩欧美在线观看| 久久久精品免费免费高清| 国产精品成人在线| 99国产精品99久久久久| 亚洲精品一区蜜桃| 国产亚洲精品久久久久5区| 成人三级做爰电影| 亚洲精品第二区| 五月开心婷婷网| 十分钟在线观看高清视频www| 丰满饥渴人妻一区二区三| 久久久久精品人妻al黑| 一区在线观看完整版| 真人做人爱边吃奶动态| 精品人妻1区二区| 狠狠婷婷综合久久久久久88av| av又黄又爽大尺度在线免费看| 亚洲精品一卡2卡三卡4卡5卡 | 高清黄色对白视频在线免费看| 午夜成年电影在线免费观看| 日韩中文字幕欧美一区二区| 中文字幕人妻丝袜一区二区| 国产麻豆69| 美女国产高潮福利片在线看| 亚洲精品一二三| 如日韩欧美国产精品一区二区三区| 精品一区二区三区av网在线观看 | 人人妻人人添人人爽欧美一区卜| 热99国产精品久久久久久7| 久久精品aⅴ一区二区三区四区| 亚洲激情五月婷婷啪啪| 国产一区二区 视频在线| 亚洲国产成人一精品久久久| 黄色视频不卡| 各种免费的搞黄视频| 美女国产高潮福利片在线看| 午夜免费观看性视频| 十八禁人妻一区二区| 亚洲欧美日韩高清在线视频 | 精品人妻一区二区三区麻豆| 狠狠婷婷综合久久久久久88av| 男女国产视频网站| 超碰97精品在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 欧美国产精品一级二级三级| 久久久久国产精品人妻一区二区| videosex国产| 午夜成年电影在线免费观看| 老司机靠b影院| 国产男女超爽视频在线观看| 老司机午夜福利在线观看视频 | 18禁黄网站禁片午夜丰满| 久久人人爽av亚洲精品天堂| 亚洲成人免费电影在线观看| h视频一区二区三区| 国产国语露脸激情在线看| 欧美性长视频在线观看| a级片在线免费高清观看视频| 亚洲天堂av无毛| 无限看片的www在线观看| 少妇被粗大的猛进出69影院| 久久精品久久久久久噜噜老黄| 亚洲av成人一区二区三| 精品视频人人做人人爽| 美女主播在线视频| 欧美精品啪啪一区二区三区 | 中文欧美无线码| 国产免费视频播放在线视频| 丁香六月欧美| 午夜影院在线不卡| 嫁个100分男人电影在线观看| 精品少妇久久久久久888优播| 91字幕亚洲| 精品免费久久久久久久清纯 | 亚洲成人国产一区在线观看| 一进一出抽搐动态| 99国产精品一区二区三区| 日韩欧美一区二区三区在线观看 | 在线观看舔阴道视频| 日韩欧美一区视频在线观看| 国产熟女午夜一区二区三区| 美国免费a级毛片| 国产不卡av网站在线观看| 美女视频免费永久观看网站| 亚洲人成77777在线视频| 精品欧美一区二区三区在线| 亚洲精品成人av观看孕妇| 日韩熟女老妇一区二区性免费视频| 一本色道久久久久久精品综合| 亚洲中文字幕日韩| 亚洲av电影在线进入| www.熟女人妻精品国产| 狠狠精品人妻久久久久久综合| 亚洲精品中文字幕一二三四区 | 一级毛片女人18水好多| 国产高清视频在线播放一区 | 手机成人av网站| 亚洲男人天堂网一区| 极品人妻少妇av视频| 久热爱精品视频在线9| 久久精品人人爽人人爽视色| 欧美一级毛片孕妇| 亚洲精品中文字幕在线视频| 国产亚洲精品久久久久5区| 黄色a级毛片大全视频| av欧美777| 99香蕉大伊视频| 亚洲一区中文字幕在线| 国产精品二区激情视频| 老司机福利观看| 亚洲国产欧美日韩在线播放| 老鸭窝网址在线观看| 黄色毛片三级朝国网站| 日韩精品免费视频一区二区三区| 精品熟女少妇八av免费久了| av不卡在线播放| a级片在线免费高清观看视频| 精品欧美一区二区三区在线| 老汉色∧v一级毛片| 啪啪无遮挡十八禁网站| 脱女人内裤的视频| 午夜免费观看性视频| 精品久久久久久电影网| 天堂俺去俺来也www色官网| 狠狠婷婷综合久久久久久88av| 亚洲国产精品一区三区| 亚洲 国产 在线| 国产精品 国内视频| 国产在视频线精品| 深夜精品福利| 亚洲avbb在线观看| 亚洲精品一区蜜桃| 国产在视频线精品| 啦啦啦视频在线资源免费观看| 欧美在线黄色| 美女主播在线视频| 国产精品一区二区精品视频观看| 亚洲人成77777在线视频| 国产一区二区 视频在线| 欧美av亚洲av综合av国产av| 免费观看a级毛片全部| 曰老女人黄片| 成人免费观看视频高清| 久久ye,这里只有精品| 欧美精品高潮呻吟av久久| 午夜精品久久久久久毛片777| 国产精品99久久99久久久不卡| 视频区欧美日本亚洲| 精品第一国产精品| 精品一区二区三卡| 国产精品秋霞免费鲁丝片| 久久久久久免费高清国产稀缺| 69av精品久久久久久 | 黄色毛片三级朝国网站| av视频免费观看在线观看| 最新在线观看一区二区三区| 侵犯人妻中文字幕一二三四区| 成年美女黄网站色视频大全免费| av在线播放精品| 蜜桃在线观看..| 国产男人的电影天堂91| 人人妻,人人澡人人爽秒播| 国产色视频综合| 欧美精品啪啪一区二区三区 | 青草久久国产| 久久久久网色| 久久中文字幕一级| 日日夜夜操网爽| 久久中文看片网| 欧美另类亚洲清纯唯美| 女人高潮潮喷娇喘18禁视频| 亚洲色图 男人天堂 中文字幕| 国产精品免费视频内射| 成年美女黄网站色视频大全免费| 少妇猛男粗大的猛烈进出视频| 人人妻人人爽人人添夜夜欢视频| 法律面前人人平等表现在哪些方面 | 在线精品无人区一区二区三| 国产精品 欧美亚洲| 91成人精品电影| 亚洲人成电影免费在线| 欧美日韩黄片免| 十八禁人妻一区二区| 欧美精品亚洲一区二区| 精品国产乱码久久久久久男人| 91精品国产国语对白视频| 欧美黑人欧美精品刺激| 夫妻午夜视频| 精品亚洲成a人片在线观看| 日韩有码中文字幕| 妹子高潮喷水视频| 视频在线观看一区二区三区| 亚洲欧美清纯卡通| 免费高清在线观看日韩| 一级毛片精品| 高清视频免费观看一区二区| 国产精品一区二区在线不卡| 两个人免费观看高清视频| 亚洲精品乱久久久久久| 亚洲伊人色综图| 欧美日韩一级在线毛片| 久久久久久免费高清国产稀缺| 1024视频免费在线观看| 亚洲专区国产一区二区| 国产亚洲精品久久久久5区| 99久久99久久久精品蜜桃| 亚洲精品美女久久久久99蜜臀| 男女国产视频网站| 亚洲精品久久成人aⅴ小说| 91麻豆av在线| 国产在线观看jvid| 操美女的视频在线观看| 女人被躁到高潮嗷嗷叫费观| 亚洲五月婷婷丁香| 又紧又爽又黄一区二区| 欧美精品人与动牲交sv欧美| 免费日韩欧美在线观看| 国产精品自产拍在线观看55亚洲 | 精品久久久久久久毛片微露脸 | 亚洲精品国产av蜜桃| 如日韩欧美国产精品一区二区三区| 手机成人av网站| 日本a在线网址| 国产在视频线精品| 久久中文字幕一级| 亚洲国产欧美网| 国产三级黄色录像| 国产1区2区3区精品| 国产精品二区激情视频| 中文字幕最新亚洲高清| 久久久精品免费免费高清| 午夜久久久在线观看| 久久狼人影院| 国产男女内射视频| 日本91视频免费播放| 免费在线观看黄色视频的| 12—13女人毛片做爰片一| 免费高清在线观看视频在线观看| 手机成人av网站| 亚洲黑人精品在线| 手机成人av网站| 99国产精品一区二区蜜桃av | 天天操日日干夜夜撸| 在线观看舔阴道视频| 欧美日本中文国产一区发布| 亚洲av电影在线观看一区二区三区| 国产不卡av网站在线观看| 搡老乐熟女国产| 亚洲国产精品999| 久久综合国产亚洲精品| 深夜精品福利| 脱女人内裤的视频| 99re6热这里在线精品视频| 50天的宝宝边吃奶边哭怎么回事| 国产伦理片在线播放av一区| 久久热在线av| 99精品欧美一区二区三区四区| 久久综合国产亚洲精品| 精品人妻熟女毛片av久久网站| 国产99久久九九免费精品| av有码第一页| 亚洲av成人一区二区三| 日韩欧美免费精品| 在线十欧美十亚洲十日本专区| 69av精品久久久久久 | 老司机深夜福利视频在线观看 | 王馨瑶露胸无遮挡在线观看| 久久久国产精品麻豆| 国产成人a∨麻豆精品| 19禁男女啪啪无遮挡网站| 久久 成人 亚洲| 久久久国产精品麻豆| 亚洲欧美清纯卡通| 永久免费av网站大全| 欧美精品亚洲一区二区| 欧美日韩成人在线一区二区| 婷婷色av中文字幕| 18禁黄网站禁片午夜丰满| 国产深夜福利视频在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美97在线视频| 久久久精品94久久精品| 精品乱码久久久久久99久播| 国产精品一区二区免费欧美 | 黑人猛操日本美女一级片| kizo精华| 看免费av毛片| 国产一区二区三区在线臀色熟女 | 国产视频一区二区在线看| 黄片大片在线免费观看| 超碰97精品在线观看| bbb黄色大片| 老熟女久久久| e午夜精品久久久久久久| 国产真人三级小视频在线观看| 欧美另类亚洲清纯唯美| 国产91精品成人一区二区三区 | 精品第一国产精品| 99精品欧美一区二区三区四区| 国精品久久久久久国模美| 国产av精品麻豆| 日韩免费高清中文字幕av| 99国产精品一区二区蜜桃av | 99香蕉大伊视频| 成年av动漫网址| 99久久精品国产亚洲精品| 欧美+亚洲+日韩+国产| 久久久久视频综合| av视频免费观看在线观看| 丰满迷人的少妇在线观看| 亚洲中文av在线| 女人爽到高潮嗷嗷叫在线视频|