基于衛(wèi)星圖像的廠區(qū)尺度溫室氣體排放核算模型：以制漿造紙企業(yè)為例

摘要： 本研究以制漿造紙企業(yè)為例，構(gòu)建了一種融合了衛(wèi)星圖像和機(jī)器學(xué)習(xí)方法的廠區(qū)尺度溫室氣體排放核算模型。首先通過衛(wèi)星圖像將制漿造紙企業(yè)分成原生纖維制漿企業(yè)、高定量紙張生產(chǎn)企業(yè)、低定量紙張生產(chǎn)企業(yè)和特種紙企業(yè)4類；然后采用Lasso算法，將廠區(qū)面積與碳排放進(jìn)行擬合，4種不同類型企業(yè)的R2均在0. 74以上。采用碳排放核算模型對(duì)全國重點(diǎn)排放單位名單中的112家造紙企業(yè)，進(jìn)行了廠區(qū)尺度的溫室氣體排放核算。結(jié)果表明，112家造紙行業(yè)重點(diǎn)排放單位之間溫室氣體排放差異較大，最小排放量為4. 45萬t CO2eq，最大排放量達(dá)426. 01萬t CO2eq，該方法也為廠區(qū)尺度的碳排放核算提供了一種新思路。

關(guān)鍵詞：衛(wèi)星圖像；機(jī)器學(xué)習(xí)；廠區(qū)尺度；溫室氣體排放；制漿造紙企業(yè)

中圖分類號(hào)：TS7 文獻(xiàn)標(biāo)識(shí)碼：A DOI：10. 11980/j. issn. 0254-508X. 2025. 02. 008

在應(yīng)對(duì)全球氣候變化的過程中，溫室氣體排放核算是實(shí)現(xiàn)可持續(xù)發(fā)展的關(guān)鍵任務(wù)之一。在“雙碳”背景下，我國工業(yè)企業(yè)面臨節(jié)能減排壓力。明確排放責(zé)任并落實(shí)減排任務(wù)，首要步驟是開展廠區(qū)尺度溫室氣體排放核算。在工業(yè)生產(chǎn)過程中，由于生產(chǎn)流程的復(fù)雜性和多樣性，會(huì)直接和間接排放溫室氣體 [1]。因此，如何從廠區(qū)尺度核算工業(yè)溫室氣體排放，是亟待解決的問題。

目前，溫室氣體核算方法包括實(shí)測(cè)法和排放因子法[2]。其中，實(shí)測(cè)法依賴于先進(jìn)的儀器設(shè)備，直接測(cè)量溫室氣體通量。Faubert等[3]采取實(shí)測(cè)法，通過現(xiàn)場(chǎng)測(cè)量中試規(guī)模的CH4通量，量化廢物填埋廠的溫室氣體排放情況，結(jié)果表明，與直接土地利用相比，填埋的溫室氣體排放量增加了3倍。在測(cè)量CH4通量的基礎(chǔ)上，部分研究人員開發(fā)了同時(shí)測(cè)量CO2、CH4和N2O3種溫室氣體通量的方法，并更新碳排放因子，已在加拿大某大型工廠進(jìn)行實(shí)驗(yàn)[4]。盡管實(shí)測(cè)法具有實(shí)時(shí)性高、檢測(cè)準(zhǔn)確的優(yōu)勢(shì)，但是受限于高昂的設(shè)備成本以及生產(chǎn)現(xiàn)場(chǎng)復(fù)雜環(huán)境帶來的測(cè)量誤差等問題，實(shí)測(cè)法并未進(jìn)行大規(guī)模推廣。

排放因子法是我國工業(yè)企業(yè)碳排放核算過程中使用最普遍的一種方法，主要是通過活動(dòng)水平數(shù)據(jù)乘以碳排放因子，得到對(duì)應(yīng)的碳排放數(shù)據(jù)[5]。其中，活動(dòng)水平數(shù)據(jù)可以通過省市、行業(yè)的統(tǒng)計(jì)數(shù)據(jù)或者企業(yè)的實(shí)測(cè)數(shù)據(jù)獲得；碳排放因子可以通過企業(yè)實(shí)測(cè)或參考標(biāo)準(zhǔn)指南提供的默認(rèn)數(shù)據(jù)獲得。已有大量研究學(xué)者采用排放因子法，對(duì)煤炭企業(yè)[6]、制藥企業(yè)[7]、鋼鐵企業(yè)[8]、火力發(fā)電工廠[9]等不同類型工廠，進(jìn)行了碳排放核算。但在許多行業(yè)中，企業(yè)能耗數(shù)據(jù)收集困難，導(dǎo)致溫室氣體排放核算不準(zhǔn)確，且自報(bào)數(shù)據(jù)可能存在偏差，難以反映真實(shí)的溫室氣體排放情況。此外，傳統(tǒng)的排放因子法通?；谔囟ㄐ袠I(yè)或生產(chǎn)過程的統(tǒng)計(jì)數(shù)據(jù)，缺乏普適性，導(dǎo)致其在實(shí)際應(yīng)用中難以滿足不同行業(yè)或企業(yè)的需求[10]。

造紙工業(yè)作為能源密集型產(chǎn)業(yè)，在生產(chǎn)過程中會(huì)排放大量CO2，已有部分學(xué)者對(duì)造紙企業(yè)進(jìn)行了碳排放核算。Zhao等[11]依據(jù)造紙企業(yè)提供的生產(chǎn)過程的用能數(shù)據(jù)，核算了斯道拉恩索（廣西） 漿紙有限公司的備料、制漿、造紙、堿回收4個(gè)生產(chǎn)過程的溫室氣體排放量。Man等[12]收集了我國23家造紙企業(yè)42條生產(chǎn)線的實(shí)際能耗數(shù)據(jù)，核算了42條生產(chǎn)線的實(shí)際碳排放數(shù)據(jù)，結(jié)果表明，我國造紙企業(yè)溫室氣體排放水平遠(yuǎn)優(yōu)于國家標(biāo)準(zhǔn)。Kuparinen 等[13]從生命周期角度出發(fā)，核算了北歐硫酸鹽針葉木紙漿廠和南美硫酸鹽桉木紙漿廠的碳排放量，并根據(jù)全球紙漿生產(chǎn)數(shù)據(jù)估計(jì)了全球碳捕獲潛力。由于排放因子更新速度慢，研究人員將目光投向準(zhǔn)確性更高、可以動(dòng)態(tài)預(yù)測(cè)溫室氣體排放的機(jī)器學(xué)習(xí)法[14]。在造紙廢水處理過程中，降解污染物會(huì)產(chǎn)生CO2、CH4和N2O等溫室氣體并直接排放[15]，這引起了研究人員的廣泛關(guān)注。由于造紙廢水處理過程具備長期非線性動(dòng)態(tài)環(huán)境的特征，基于PSO-BP[16]和DNN-LSTM的深度神經(jīng)網(wǎng)絡(luò)模型[17]，已被開發(fā)用于預(yù)測(cè)造紙廢水處理過程中的溫室氣體排放情況，并識(shí)別影響溫室氣體排放量的關(guān)鍵變量。但深度神經(jīng)網(wǎng)絡(luò)建模過程需要大量的樣本點(diǎn)，并且可能比較耗時(shí)[18]，因此計(jì)算資源使用更少、計(jì)算速度更快的溫室氣體代理模型急需被開發(fā)出來，以用于預(yù)測(cè)造紙廢水處理的溫室氣體排放量[19]。綜上所述，造紙工業(yè)和其他工業(yè)均面臨著相同的碳排放核算問題，無論是排放因子法還是機(jī)器學(xué)習(xí)方法，由于數(shù)據(jù)收集困難，且均是針對(duì)單一企業(yè)，難以對(duì)企業(yè)進(jìn)行大規(guī)模的溫室氣體排放核算。

隨著地理信息技術(shù)的發(fā)展，通過衛(wèi)星圖像可以觀察城市或者工廠的空間布局、建筑分布等，完成對(duì)建筑形態(tài)方面的分析，并通過對(duì)建筑物的位置、數(shù)量及面積的分析，核算溫室氣體排放[20]。部分研究人員通過對(duì)比不同時(shí)間點(diǎn)的衛(wèi)星圖像，追蹤城市[21]和工業(yè)園區(qū)[22]內(nèi)建筑物的變化，并據(jù)此估算其碳排放量的動(dòng)態(tài)變化。以上研究表明，利用衛(wèi)星圖像進(jìn)行溫室氣體核算是可行的方法。造紙廠具有占地面積大、涉及多個(gè)生產(chǎn)和處理設(shè)施的特點(diǎn)。通過衛(wèi)星圖像，可以有效監(jiān)測(cè)和分析造紙廠的空間布局和設(shè)施分布[23]，從而估算其溫室氣體排放量。如原材料存儲(chǔ)區(qū)、生產(chǎn)車間和廢水處理設(shè)施等結(jié)構(gòu)或建筑的面積和位置數(shù)據(jù)，均可以通過衛(wèi)星圖像獲取和分析，為廠區(qū)尺度溫室氣體核算提供可靠數(shù)據(jù)支持。

本研究以制漿造紙企業(yè)為例，通過融合地理信息技術(shù)和機(jī)器學(xué)習(xí)算法，利用衛(wèi)星圖像探究造紙廠區(qū)不同建筑形態(tài)與溫室氣體排放之間的關(guān)聯(lián)關(guān)系，建立了造紙企業(yè)廠區(qū)尺度的溫室氣體核算模型。不僅利用建立的溫室氣體核算模型，核算多個(gè)造紙企業(yè)的廠區(qū)尺度溫室氣體排放，從廠區(qū)尺度估算2022年全國造紙行業(yè)重點(diǎn)排放單位的溫室氣體排放，也為其他工業(yè)進(jìn)行廠區(qū)尺度溫室氣體核算提供了新思路。

1 基于衛(wèi)星圖像的造紙廠區(qū)溫室氣體核算模型

1. 1 數(shù)據(jù)收集

本研究收集的數(shù)據(jù)包含4個(gè)部分，第1部分?jǐn)?shù)據(jù)為制漿造紙上市企業(yè)名錄數(shù)據(jù)，來源于《2022中國造紙年鑒》；第2部分?jǐn)?shù)據(jù)為企業(yè)POI數(shù)據(jù)，來源于百度地圖開放平臺(tái)；第3部分?jǐn)?shù)據(jù)為企業(yè)衛(wèi)星圖像數(shù)據(jù)，來源于國家地理信息公共服務(wù)平臺(tái)；第4部分為企業(yè)碳排放數(shù)據(jù)，來源于企業(yè)的環(huán)境、社會(huì)和公司治理報(bào)告、可持續(xù)性發(fā)展報(bào)告和年度報(bào)告。通過POI數(shù)據(jù)將碳排放數(shù)據(jù)、衛(wèi)星圖像數(shù)據(jù)進(jìn)行匹配。

1. 2 造紙企業(yè)建筑物特征分類

造紙企業(yè)的不同區(qū)域在衛(wèi)星圖像上具有明顯的特征差異，這些區(qū)域在建筑規(guī)模、形狀和顏色等方面各具特點(diǎn)，可將造紙廠區(qū)劃分為以下5個(gè)不同的區(qū)域，如圖1所示。

1） 木片堆垛區(qū)域。木片堆垛區(qū)域邊界較為清晰，在衛(wèi)星圖像中與周圍環(huán)境形成顯明對(duì)比，顏色為黃色或深棕色，形狀為圓形或方形。為滿足造紙企業(yè)長時(shí)間的原料需求，木片堆垛區(qū)域的規(guī)模通常較大，能達(dá)到廠區(qū)面積的30%左右。

2） 廢水處理區(qū)域。廢水處理區(qū)域有多個(gè)大型的沉淀池、氧化池和處理池，這些水體區(qū)域在衛(wèi)星圖像中通常呈現(xiàn)規(guī)則的圓形或方形；此外由于不同處理階段的廢水水體顏色不同，如未處理的廢水可能水體顏色較深，而經(jīng)過處理的水體顏色會(huì)變淺，水體顏色呈現(xiàn)階段性變化也是廢水處理區(qū)域明顯的特征之一。

3） 熱電站區(qū)域。熱電站區(qū)域通常具備高大的煙囪，這些煙囪在衛(wèi)星圖像中通常是細(xì)長的陰影或高亮的柱狀物；另一種常見的大型建筑冷卻塔，通常是圓形或橢圓形的結(jié)構(gòu)，頂部開口，這2種建筑冒出的蒸汽云，在衛(wèi)星圖像中通常呈現(xiàn)為白色或淺色的云狀物。

4） 其他堆垛區(qū)域。其他堆垛區(qū)域的物料通常以規(guī)則的方式、采用小包的類型堆放，在衛(wèi)星圖像中呈現(xiàn)出明顯的幾何形狀，一般為矩形，且不同類型的物料呈現(xiàn)出不同的顏色，物料之間存在明顯邊界，上有覆蓋物進(jìn)行雨水遮擋。

5） 建筑區(qū)域。建筑區(qū)域包含制漿車間、造紙車間、倉庫、辦公樓等建筑，在衛(wèi)星圖像上，表現(xiàn)為密集的建筑群，呈現(xiàn)為大片的連續(xù)色塊，與其他功能區(qū)域如其他堆垛區(qū)域和廢水處理區(qū)域等有明顯的分隔。

1. 3 造紙企業(yè)廠區(qū)形態(tài)特征分類

由于造紙企業(yè)生產(chǎn)原料、產(chǎn)品種類、生產(chǎn)規(guī)模和所處地區(qū)的不同，廠區(qū)包含的建筑存在差異，在衛(wèi)星圖像上呈現(xiàn)出不同的形態(tài)特征。根據(jù)造紙企業(yè)廠區(qū)的形態(tài)特征，可以將造紙企業(yè)劃分成4類。第1類造紙企業(yè)廠區(qū)為原生纖維制漿企業(yè)。原生纖維制漿企業(yè)的形態(tài)特征為存在較大的木片堆垛區(qū)域，如圖2（a）所示。第2類造紙企業(yè)廠區(qū)為高定量紙張生產(chǎn)企業(yè)。高定量紙張生產(chǎn)企業(yè)主要生產(chǎn)包裝用紙、文化用紙和工業(yè)用紙等產(chǎn)品，不包含原料制漿，而是以廢紙制漿或外購漿板為主，在衛(wèi)星圖像上展現(xiàn)較大的其他堆垛區(qū)域，且顏色較為多樣，擁有細(xì)長的生產(chǎn)線和較大的倉儲(chǔ)空間，如圖2（b）所示。第3類造紙企業(yè)廠區(qū)為低定量紙張生產(chǎn)企業(yè)。低定量紙張生產(chǎn)企業(yè)的主要產(chǎn)品為生活用紙，如圖2（c）所示，相比原生纖維制漿企業(yè)和高定量紙張生產(chǎn)企業(yè)，低定量紙張生產(chǎn)企業(yè)整體廠區(qū)規(guī)模較小，擁有一定的其他堆垛區(qū)域，在衛(wèi)星圖像中堆垛物體較為方正且多呈現(xiàn)為綠色和藍(lán)色，建筑物呈現(xiàn)網(wǎng)格狀。第4類企業(yè)為特種紙企業(yè)。特種紙企業(yè)主要產(chǎn)品為特種紙，如圖2（d）所示，其衛(wèi)星圖像上幾乎只包含建筑區(qū)域，建筑物同樣呈現(xiàn)網(wǎng)格狀，無其他堆垛區(qū)域。

1. 4 基于衛(wèi)星圖像的造紙廠區(qū)溫室氣體核算模型的建立

根據(jù)1.3 中的企業(yè)形態(tài)特征分類方式，將造紙企業(yè)分成4類，分別構(gòu)建造紙廠區(qū)碳排放核算模型?；?.2 中的建筑特征，采用地理信息系統(tǒng)（GIS） 工具進(jìn)行輔助分析，對(duì)于不同區(qū)域的占地面積進(jìn)行測(cè)算，為構(gòu)建溫室氣體排放模型提供空間數(shù)據(jù)。

本研究的碳排放模型是基于多元線性回歸法建立的，即采用機(jī)器學(xué)習(xí)算法中的最小絕對(duì)收縮和選擇算子（least absolute shrinkage and selection operator，LASSO） 算法，構(gòu)建碳排放量與廠區(qū)建筑面積之間的函數(shù)關(guān)系。LASSO算法的核心是在損失函數(shù)中引入1個(gè)正則化系數(shù)，來減少模型復(fù)雜度，同時(shí)保持預(yù)測(cè)的準(zhǔn)確性，損失函數(shù)如式（1）所示。

式中，J （θ ）為損失函數(shù)；θ 為回歸系數(shù)；m 為樣本數(shù)量；xi 和yi 分別為x 和y 的第i 個(gè)值；λ 為正則化系數(shù)。

2 全國造紙重點(diǎn)碳排放單位碳排放核算

通過收集各?。ㄖ陛犑?、自治區(qū)） 生態(tài)環(huán)境廳發(fā)布2022年度納入全國碳排放權(quán)交易配額管理的重點(diǎn)排放單位名單數(shù)據(jù)，除遼寧省、香港特別行政區(qū)、澳門特別行政區(qū)、臺(tái)灣地區(qū)、西藏自治區(qū)外，全國2022年度重點(diǎn)排放單位共2 206家，其中制漿造紙行業(yè)114家，約占5.17%，2家造紙企業(yè)未找到相關(guān)圖像數(shù)據(jù)。各?。ㄖ陛犑?、自治區(qū)） 造紙企業(yè)分布如圖3所示。由圖3可知，共計(jì)23個(gè)?。ㄖ陛犑?、自治區(qū)） 公布的重點(diǎn)排放單位名單中包含造紙企業(yè)，其中廣東省有24家造紙企業(yè)為重點(diǎn)排放單位，屬于全國最多的?。ㄖ陛犑?、自治區(qū)）。

收集上述112家造紙企業(yè)的衛(wèi)星圖像數(shù)據(jù)，根據(jù)1.2 中的方法對(duì)造紙企業(yè)廠區(qū)進(jìn)行建筑特征的劃分，隨后依據(jù)1.3 中的廠區(qū)分類方式，將112家造紙企業(yè)進(jìn)行分類并統(tǒng)計(jì)數(shù)量，結(jié)果如圖4 所示。由圖4 可知，原生纖維制漿企業(yè)和高定量紙張生產(chǎn)企業(yè)數(shù)量較多，分別有37家和62家，表明這2種類型造紙企業(yè)生產(chǎn)過程中的碳排放量較高。

通過天地圖平臺(tái)分別測(cè)算112家造紙企業(yè)的木片堆垛區(qū)域面積、其他堆垛區(qū)域面積、廢水處理區(qū)域面積、熱電站區(qū)域面積和建筑區(qū)域面積，并根據(jù)造紙企業(yè)分類輸入1.4 中所搭建好的碳排放核算模型。

3 結(jié)果與分析

3. 1 溫室氣體核算模型結(jié)果

基于Lasso回歸模型，得到部分造紙企業(yè)廠區(qū)尺度的碳排放擬合值，并建立散點(diǎn)圖，檢驗(yàn)碳排放擬合值與真實(shí)值的相關(guān)性，結(jié)果如圖5所示。由圖5可知，當(dāng)碳排放擬合值和真實(shí)值相同時(shí)，樣本點(diǎn)均分布在y=x的對(duì)角線上。基于Lasso回歸模型，建立的4種基于廠區(qū)面積的碳排放核算模型，決定系數(shù)（R2） 和百分比誤差（MAPE） 分別為0.99和9.10%（原生纖維制漿企業(yè)）；0.91和21.24%（高定量紙張生產(chǎn)企業(yè)）；0.74和25.17%（低定量紙張生產(chǎn)企業(yè)）；0.88和32.21%（特種紙企業(yè)）；表明碳排放核算模型擬合結(jié)果良好。

3. 2 全國造紙重點(diǎn)碳排放單位碳排放核算

3. 2. 1 ?。ㄖ陛犑?、自治區(qū)） 碳排放分析

根據(jù)本研究建立的碳排放核算模型，估算2022年全國造紙行業(yè)重點(diǎn)碳排放單位的碳排放總量為8 126.37萬t CO2eq，各?。ㄖ陛犑?、自治區(qū)） 的碳排放情況如表1所示。由表1可知，2022年碳排放總量前3的省份分別為廣東省、山東省和江蘇省，碳排放總量分別為1 417.71萬t CO2eq、1 027.70萬t CO2eq和894.78萬t CO2eq，分別占全國碳排放總量的17.45%、11.01%和8.95%。山東省作為全國紙漿、紙和紙板產(chǎn)量前2的造紙工業(yè)大省，其碳排放總量遠(yuǎn)低于廣東省，主要原因是山東省發(fā)布的全國重點(diǎn)碳排放單位名單中屬于造紙行業(yè)的企業(yè)僅有10 家，遠(yuǎn)低于廣東省的24家。但從省（直轄市、自治區(qū)） 單位企業(yè)碳排放量來看，山東省單位企業(yè)碳排放量為114.15萬t CO2eq，大于廣東省的59.07萬t CO2eq。

綜上所述，中國造紙行業(yè)的碳排放量具有顯著的區(qū)域差異，東部和南部沿海地區(qū)由于工業(yè)化程度高，企業(yè)密集，碳排放量也較高，這些地區(qū)需要更嚴(yán)格的環(huán)保政策和技術(shù)改進(jìn)來控制和減少碳排放量。

3. 2. 2 行業(yè)碳排放分析

從企業(yè)種類來看，原生纖維制漿企業(yè)碳排放量為5 095.93萬t CO2eq，占全國碳排放總量的59.01%；高定量紙張生產(chǎn)企業(yè)碳排放量為2 887.95萬t CO2eq，占全國碳排放總量的35.54%，這2種類型企業(yè)的碳排放總量約占全國重點(diǎn)排放單位名單中的制漿造紙企業(yè)碳排放總量的95%。

圖6為4種類型企業(yè)的碳排放量箱型圖，展示了不同分類企業(yè)的碳排放量的統(tǒng)計(jì)分布情況。由圖6可知，原生纖維制漿企業(yè)的碳排放量中位數(shù)最高，約為100 萬t CO2eq，且數(shù)據(jù)分布較廣，存在多個(gè)極端高值，最大值超過400萬t CO2eq。高定量紙張產(chǎn)品企業(yè)的碳排放量中位數(shù)約為50萬t CO2eq，分布較為集中，但也存在一些極端高值。低定量紙張生產(chǎn)企業(yè)的碳排放量中位數(shù)最低，其數(shù)據(jù)分布最為密集且無極端高值。特種紙企業(yè)的碳排放量中位數(shù)略高于低定量企業(yè)，約為20萬t CO2eq，數(shù)據(jù)分布相對(duì)集中。綜上所述，不同類型造紙企業(yè)在碳排放量上存在顯著差異，其中原生纖維制漿企業(yè)的碳排放量顯著高于其他類型企業(yè)。

3. 2. 3 企業(yè)碳排放分析

核算結(jié)果表明，核算中國東南部沿海地區(qū)高碳排放量的造紙企業(yè)較為集中，在這些高碳排放量企業(yè)中，原生纖維制漿企業(yè)和高定量企業(yè)的分布尤為顯著。如江蘇省、浙江省、山東省和廣東省等地區(qū)聚集了大量的原生纖維制漿企業(yè)和高定量紙張生產(chǎn)企業(yè)，存在大型造紙廠區(qū)，碳排放量大而聚集。同時(shí)，四川省和重慶市2個(gè)地區(qū)，也有大規(guī)模的造紙企業(yè)集群。此外，盡管特種紙企業(yè)的數(shù)量相對(duì)較少，但其碳排放量普遍在50萬t CO2eq以上，主要集中在江蘇省。低定量紙張生產(chǎn)企業(yè)的碳排放量相對(duì)較低，大部分在10 萬t CO2eq 以下，且分布較分散，無明顯的集中區(qū)域。

2022年全國重點(diǎn)排放單位造紙行業(yè)碳排放量前10企業(yè)如表2所示。由表2可知，經(jīng)估算，碳排放量前10的造紙企業(yè)，其碳排放量均在200萬t CO2eq以上，且企業(yè)之間碳排放量差異較大。其中，海南金海漿紙業(yè)有限公司的碳排放量最高，達(dá)426.01萬t CO2eq，而廣西金桂漿紙業(yè)有限公司的碳排放量最低，為206.80萬t CO2eq，2家企業(yè)碳排放量差值達(dá)219.21萬t CO2eq。

4 結(jié)論

本研究針對(duì)現(xiàn)有造紙工業(yè)碳排放核算方式存在的成本高昂、反饋滯后、操作復(fù)雜和數(shù)據(jù)收集困難等問題，提出一種基于衛(wèi)星圖像的溫室氣體核算模型，可以實(shí)現(xiàn)廠區(qū)尺度碳排放核算。

4. 1 通過衛(wèi)星圖像將廠區(qū)內(nèi)的建筑物進(jìn)行特征分類，根據(jù)不同廠區(qū)所包含建筑物的形態(tài)差異，可以將造紙企業(yè)分成原生纖維制漿企業(yè)、高定量紙張生產(chǎn)企業(yè)、低定量紙張生產(chǎn)企業(yè)、特種紙企業(yè)4類。

4. 2 針對(duì)以上4種不同企業(yè)，分別建立碳排放核算模型，決定系數(shù)（R2） 和百分比誤差（MAPE） 分別為0. 99和9. 10%（原生纖維制漿企業(yè)）；0. 91和21. 24%（高定量紙張生產(chǎn)企業(yè)）；0. 74和25. 17%（低定量紙張生產(chǎn)企業(yè)）；0. 88和32. 21%（特種紙企業(yè)）。

4. 3 利用碳排放核算模型，從廠區(qū)尺度對(duì)全國重點(diǎn)排放單位名單中的112家造紙企業(yè)進(jìn)行碳排放核算。經(jīng)估算，碳排放量前10的造紙企業(yè)，其碳排放量均在200萬t CO2eq以上，且企業(yè)之間碳排放量差異較大。

參考文獻(xiàn)

［1］ 黃 煒， 王逸倫， 肖相澤， 等. 工業(yè)生產(chǎn)過程碳排放核算方法研究［J］. 統(tǒng)計(jì)科學(xué)與實(shí)踐，2024（1）：50-53.

HUANG W， WANG Y L， XIAO X Z， et al. Research on CarbonEmission Accounting Methods in Industrial Production Processes［J］. Statistical Science and Practice， 2024（1）： 50-53.

［2］ 宋佳翼， 程正柏， 劉利琴， 等. 制漿造紙行業(yè)溫室氣體排放核算標(biāo)準(zhǔn)與方法淺析［J］. 中國造紙，2023，42（8）：11-18.

SONG J Y， CHENG Z B， LIU L Q， et al. Analysis on GreenhouseGas Emission Standards and Accounting Methods in Pulp and PaperIndustry［J］. China Pulp amp; Paper， 2023，42（8）：11-18.

［3］ FAUBERT P， BéLISLE C L， BERTRAND N， et al. Landapplication of pulp and paper mill sludge may reduce greenhouse gasemissions compared to landfilling［J］. Resources， Conservation andRecycling， DOI： 10. 1016/j. resconrec. 2019. 104415.

［4］ RIBEIRO PAULA R， CUSSON M， BERTRAND N， et al.Correction factors for large-scale greenhouse gas assessment frompulp and paper mill sludge landfill sites［J］. Waste Management，2024， 177： 178-181.

［5］ 葉佳蓉. 中國造紙和紙制品業(yè)碳排放時(shí)空分布及其驅(qū)動(dòng)因素研究［D］. 南京：南京林業(yè)大學(xué)，2017.

YE J R. A Study on the Spatiotemporal Distribution of Carbon Emis?sions and Its Driving Factors in the Chinese Papermaking and PaperProducts Industry［D］. Nanjing： Nanjing Forestry University，2017.

［6］ LI Y， JIN X Y， WANG G L， et al. Construction and Application ofa Carbon Emission Model for China’s Coal Production Enterprisesand Result Analysis［J］. Frontiers in Energy Research， DOI：10. 3389/fenrg. 2022. 889877.

［7］ SHANG G D， YANG L， WANG J， et al. Case study on carbonemission of a traditional Chinese medicine pharmaceutical enterprisetoward carbon neutrality［J］. Journal of Cleaner Production， DOI：10. 1016/j. jclepro. 2024. 141600.

［8］ XU R， TONG D， DAVIS S J， et al. Plant-by-plant decarbonizationstrategies for the global steel industry［J］. Nature Climate Change，2023， 13（10）：1067-1074.

［9］ 馮 欣，李 璇，曹安琪. 遼寧省火力發(fā)電企業(yè)碳排放核算及減排建議［J］. 環(huán)境保護(hù)與循環(huán)經(jīng)濟(jì)，2024，44（5）：1-6.

FENG X， LI X， CAO A Q. Carbon Emission Accounting and ReductionSuggestions for Thermal Power Enterprises in Liaoning Province［J］.Environmental Protection and Circular Economy， 2024， 44（5）： 1-6.

［10］ 何旭丹， 王煥松， 賈學(xué)樺， 等. 中國造紙和紙制品行業(yè)碳排放特征及減排路徑分析［J］. 中國造紙， 2023，42（11）：144-151.

HE X D， WANG H S， JIA X H， et al. Analysis of Carbon EmissionCharacteristics and Emission Reduction Path in China’s Paper andPaper Products Industry［J］. China Pulp amp; Paper， 2023，42（11）：144-151.

［11］ ZHAO Q， KUANG X. A Study on the Implementation of CorporateCarbon Emission Accounting in the Production Process［J］. PolishJournal of Environmental Studies， 2024， 33（3）：2941-2957.

［12］ MAN Y， YAN Y， WANG X， et al. Overestimated carbon emission ofthe pulp and paper industry in China［J］. Energy， DOI： 10. 1016/j.energy. 2023. 127279.

［13］ KUPARINEN K， VAKKILAINEN E， TYJALA T. Biomass-basedcarbon capture and utilization in kraft pulp mills［J］. Mitigation andAdaptation Strategies for Global Change， 2019， 24（7）：1213-1230.

［14］ ASSAEL J， HEURTEBIZE T， CARLIER L， et al. GreenhouseGases Emissions： Estimating Corporate Non-reported EmissionsUsing Interpretable Machine Learning［J］. Sustainability， DOI：10. 3390/su15043391.

［15］ 張 欣， 張 放， 蔡 慧， 等. 典型制漿造紙廠的CO2排放及碳強(qiáng)度的算法和分析［J］. 中國造紙，2019，34（1）：36-42.

ZHANG X， ZHANG F， CAI H， et al. Calculation and Analysis ofCO2 Emissions and Carbon Intensity of a Typical Integrated PaperMill in China［J］. China Pulp amp; Paper， 2019， 34（1）： 36-42.

［16］ 李文青. 造紙污水A/O處理過程N(yùn)2O排放動(dòng)力學(xué)模型構(gòu)建及減排優(yōu)化控制研究［D］. 廣州：華南理工大學(xué)，2023.

LI W Q. Construction of the N2O Emission Kinetic Model and Opti?mization of Emission Reduction Control in the A/O Treatment Pro?cess of Papermaking Wastewater［D］. Guangzhou： South ChinaUniversity of Technology， 2023.

［17］ 李世忠， 滿 奕， 何正磊. 基于DNN-LSTM 的造紙廢水處理過程溫室氣體排放分析模型［J］. 中國造紙，2024，43（4）：170-176.

LI S Z， MAN Y， HE Z L. Analysis Model of Greenhouse Gas Emis?sions from Papermaking Wastewater Treatment Process Based onDNN-LSTM［J］. China Pulp amp; Paper， 2024， 43（4）： 170-176.

［18］ SUN L， GAO H， PAN S， et al. Surrogate modeling for fluid flowsbased on physics-constrained deep learning without simulation data［J］. Computer Methods in Applied Mechanics and Engineering，DOI： 10. 1016/j. cma. 2019. 112732.

［19］ 鄭鴻澤， 何正磊， 洪蒙納， 等. 造紙廢水處理過程溫室氣體排放代理模型研究［J］. 中國造紙學(xué)報(bào)，2023，38（2）：59-68.

ZHENG H Z， HE Z L， HONG M N， et al. Study on SurrogateModel of Greenhouse Gas Emission in Papermaking WastewaterTreatment Process［J］. China Pulp amp; Paper， 2023， 38（2）： 59-68.

［20］ LEVIN N， KYBA C C M， ZHANG Q， et al. Remote sensing ofnight lights： A review and an outlook for the future［J］. RemoteSensing of Environment， DOI： 10. 1016/j. rse. 2019. 111443.

［21］ WANG M， WANG Y， WU Y， et al. Identifying the spatialheterogeneity in the effects of the construction land scale on carbonemissions： Case study of the Yangtze River Economic Belt， China［J］.Environmental Research， DOI： 10. 1016/j. envres. 2022. 113397.

［22］ HANNA J， BORTH D， MOMMERT M. Physics-guided MultitaskLearning for Estimating Power Generation and CO2 Emissions fromSatellite Imagery［J］. IEEE Transactions on Geoscience and RemoteSensing， DOI： 10. 1109/TGRS. 2023. 3286444.

［23］ BONDUR V G， CHIMITDORZHIEV T N， DMITRIEV A V， et al.Methods of Radar Interferometry and Optical Satellite ImageProcessing to Study Negative Effects on the Environment （a CaseStudy of the Baikalsk Pulp and Paper Mill）［J］. Izvestiya，Atmospheric and Ocean Physics， 2021， 57（12）：1527-1537.

（責(zé)任編輯：魏琳珊）

基金項(xiàng)目：制漿造紙工程國家重點(diǎn)實(shí)驗(yàn)室開放基金資助項(xiàng)目（202419）。