黑龍江省土地利用碳平衡分區(qū)研究

關(guān)鍵詞土地利用；碳排放；生態(tài)承載力；碳平衡分區(qū)

中圖分類號K903 文獻(xiàn)標(biāo)識碼A 文章編號 1007-7731（2025）17-0069-07

DOI號10.16377/j.cnki.issn1007-7731.2025.17.018

Research on the zoning of land use carbon balance in Heilongjiang Province

LIU Jianing1 SUN Lina1ZHAO Xingcai2 HAN Liu3 （1School of Government， Heilongjiang University， Harbin 15oo80， China; （204號 ² Heilongjiang Provincial Institute of Natural Researches Rightsand Interests Surveyand Monitoring，Harbin 150000， China; 3School of Heilongjiang Provincial Civil Affairs Vocational and Technical，Harbin 15oooo，China）

AbstractThe land use data and socio-economic data were utilized of Heilongjiang Province，and based on the ArcGIS spatial analysis method，the carbon emission Gini coeffcient，carbonemission risk，carbon emission intensity， economiccontributioncoeffcientandecological carryingcoefficient were calculatedof itsinvestigatedregions.These were used as zoning indicators toconduct carbon balance zoning of the study area and propose targeted carbon emission reduction strategies foreach type ofregion.Theresults showed thatthecarbon compensation rateofthestudyarea from 1980 to 2023 was generallylowand decreased year byyear.From 2O00 to 2O23，thecarbonemission Gini coefficint of theinvestigatedregionsallexceeded0.6.From 2O00to2023，thecarbon emissonrisk valuesof the investigatedregions generally showedan increasing trend year byyear，and the spatial diferences became increasinglyobvious.From 1980 to 2023，thecarbon emission intensityof the studyarea showedadecreasing trend yearbyyear;from 2OoO to 2023，the carbon emisson intensityof each city decreased to varying degrees.From 20O to 2O23，the economic contribution coeffcient wasinarelativelyfluctuating state，and theecological carrying coefficient was generallylow.Basedon the relevantanalysis results，the carbon balanceof the studyarea was divided into low-carbon optimization areas，carbon total control areasandcarbon sink function areas.The low-carbonoptimization area includesarea A，where the carbon emissionrisk isatamedium tohighlevel，thecarbonemissonintensityislow，theeconomiccontributioncoeficientis atarelativelyhighlevel，andtheecologicalcarryingcoeficientisatarelativelylowlevel.Thecarbontotalcontrolareas includeareasB，F(xiàn)，M，D，I，E，JandC.Intheseareas，thecarbonemissionriskisbasicallyatamedium tohighlevel， thecarbonemissionintensityishigh，theeconomiccontributioncoeficientisatarelativelylowlevel，andtheecological carying coeffcient is ata medium to low level.The carbon sink function areas include areas N，L，Gand K.In these areas，thecarbon emission risk is atamedium tolow level;thecarbonemision intensityis low，theeconomic contributioncoeffcientisatarelativelyhigh level;；andtheecological carrying coefficient ishigh.Implementing differentiated carbon emission reduction strategies based on the actual situation ofeach city provides a reference for the low-carbon development of land use in the study area.

Keywords land use; carbon emissions; ecological capacity; carbon balance zoning

碳排放量增加引發(fā)的全球氣候變化已受到各界的廣泛關(guān)注，碳達(dá)峰、碳中和等目標(biāo)的實施是促進(jìn)綠色低碳可持續(xù)發(fā)展的重要舉措。土地是生態(tài)系統(tǒng)碳循環(huán)過程中碳排放的重要來源，對維持碳平衡起關(guān)鍵作用，利用科學(xué)的土地規(guī)劃管理策略以有效控制碳排放是當(dāng)前研究的重點，了解土地利用格局，進(jìn)行碳平衡分區(qū)研究，因地制宜提出不同區(qū)域碳減排策略，對碳減排的實施具有重要推動作用。相關(guān)學(xué)者從不同角度開展了土地利用碳排放研究，肖磊等1Heil等2基于碳排放系數(shù)法，測算不同用地類型碳排放及碳吸收量，并結(jié)合ArcGIS軟件進(jìn)行時空分析。相關(guān)學(xué)者3-4基于對數(shù)平均迪氏指數(shù)分解法（LMDI、地理加權(quán)回歸（GWR）等模型，分析了能源消費、地區(qū)生產(chǎn)總值、人口規(guī)模、土地開發(fā)度等指標(biāo)對土地利用碳排放量的影響。運用碳排放系數(shù)法、馬爾可夫模型及灰色預(yù)測模型等，借助IDRISI軟件可對碳排放特征進(jìn)行分析及預(yù)測[5-6。利用基尼系數(shù)、經(jīng)濟(jì)貢獻(xiàn)系數(shù)及生態(tài)承載系數(shù)可對研究區(qū)域進(jìn)行分區(qū)[7-9]

土地利用碳平衡分區(qū)研究有助于研究區(qū)域制訂有針對性的減排策略，實現(xiàn)土地利用綠色低碳發(fā)展。本文以黑龍江省為研究區(qū)，依據(jù)碳排放基尼系數(shù)、碳排放風(fēng)險、碳排放強(qiáng)度、經(jīng)濟(jì)貢獻(xiàn)系數(shù)和生態(tài)承載系數(shù)等指標(biāo)，對該地區(qū)進(jìn)行分區(qū)，包括低碳優(yōu)化區(qū)、碳總量控制區(qū)及碳匯功能區(qū)，并提出不同碳平衡區(qū)域的減排優(yōu)化策略，探尋實現(xiàn)低碳土地利用的有效途徑，為實現(xiàn)“雙碳”減排目標(biāo)提供參考。

1材料與方法

1.1 研究區(qū)基本情況

研究區(qū)位于 43^°26^′-53^°33^′N，121^°11^′-135^°05^′ E。轄13個市（A～N）;是重點林區(qū)之一，森林面積、森林總蓄積量和木材生產(chǎn)量均較高，森林面積2150.6萬 hm² ，森林覆蓋率 47.3% ，森林蓄積量22.38億 m³ 。

1.2數(shù)據(jù)來源

本研究土地利用數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http：//www.resdc.cn/），該數(shù)據(jù)以Landsat-TM圖像為信息來源，通過人工目視解譯生成，空間分辨率為 1km 。研究區(qū)土地利用碳排放量及各地市土地利用碳排放量來源于相關(guān)研究成果[0]，具體數(shù)據(jù)見表1～2。社會經(jīng)濟(jì)數(shù)據(jù)來源于《中國統(tǒng)計年鑒》《黑龍江省統(tǒng)計年鑒》及各市統(tǒng)計年鑒中的地區(qū)生產(chǎn)總值，主要用于計算研究區(qū)及各地市土地利用碳排放強(qiáng)度、經(jīng)濟(jì)貢獻(xiàn)系數(shù)及生態(tài)承載系數(shù)。

表11980—2023年研究區(qū)土地利用碳排放量

表22000—2023年研究區(qū)各地市土地利用碳排放量

單位：×104t

1.3試驗方法

1.3.1碳排放基尼系數(shù) 以基尼系數(shù)為指標(biāo)，評估各地區(qū)間碳補(bǔ)償比率的變化情況，較高的碳補(bǔ)償率表明該地區(qū)具有較高的固碳潛力。計算如式（1）。

式中， Z_i，Z_j 為 i，j 區(qū)域的碳補(bǔ)償率， n 為研究區(qū)監(jiān)測點樣本數(shù)， 為研究區(qū)監(jiān)測點碳補(bǔ)償率平均值。基尼系數(shù) lt;0.4 為相對合理狀態(tài)，基尼系數(shù) gt;0.6 為不協(xié)調(diào)狀態(tài)。

1.3.2碳排放風(fēng)險 參考張杰等[12]、湯峰等[13]的研究方法，計算研究區(qū)各地區(qū)碳排放風(fēng)險指數(shù)（C_Ri），C_Ri 越高表示其研究區(qū)域碳排放風(fēng)險越高，反之則越低。計算如式（2）。

式中， C_i 為i區(qū)土地利用碳排放量；S為研究區(qū)總面積。

1.3.3碳排放強(qiáng)度 利用碳排放強(qiáng)度表示單位地區(qū)生產(chǎn)總值帶來的碳排放量。計算如式（3）。

式中， C 為研究區(qū)土地利用碳排放量; G 為折價后的地區(qū)生產(chǎn)總值，碳排放強(qiáng)度越小，表示該區(qū)域單位地區(qū)生產(chǎn)總值所產(chǎn)生的二氧化碳量越小。根據(jù)其強(qiáng)度大小分為4個等級， 0～0.50L，0.51～1.00L，1.01～ 2.00、2.00以上。

1.3.4經(jīng)濟(jì)貢獻(xiàn)系數(shù) 利用經(jīng)濟(jì)貢獻(xiàn)系數(shù) （ECC） （20表示單位地區(qū)碳生產(chǎn)力大小。ECC表示某一地區(qū)生產(chǎn)總值占全省生產(chǎn)總值的比例與該區(qū)域碳排放量占全省比值。計算如式（4）。

式中， G_i 和 G 分別表示 i 區(qū)及研究區(qū)的生產(chǎn)總值， C_i 和 C 分別表示 i 區(qū)及研究區(qū)的碳排放量。ECCgt;2 時，表明該地區(qū)經(jīng)濟(jì)貢獻(xiàn)比率大于碳排放的比率，說明碳生產(chǎn)力較高； ECClt;2 ，說明碳生產(chǎn)力較弱。

1.3.5生態(tài)承載系數(shù) 碳生態(tài)承載系數(shù)（ESC）表示某一地區(qū)碳吸收量占研究區(qū)比例與該區(qū)域碳排放量占研究區(qū)比例的比值，反映了該區(qū)域碳匯能力的大小[14]。計算如式（5）。

式中， C_Ai 和 C_A 分別表示 i 區(qū)及研究區(qū)的碳吸收量。ESC值越高，表明該地區(qū)碳匯能力越強(qiáng)；反之，說明碳匯能力較弱。

1.3.6碳平衡分區(qū) 基于研究區(qū)各地區(qū)碳排放基尼系數(shù)、碳排放強(qiáng)度、經(jīng)濟(jì)貢獻(xiàn)及生態(tài)貢獻(xiàn)等指標(biāo)的計算，參考趙榮欽等14的方法，綜合考量研究區(qū)實際情況、區(qū)域碳平衡及綠色發(fā)展，對其進(jìn)行碳平衡劃分。

2 結(jié)果與分析

2.1 碳排放基尼系數(shù)

由圖1可知，研究區(qū)碳補(bǔ)償率整體處于低位，且隨著時間的推移逐年下降，從1980年的0.30下降到2023年的0.07，說明研究區(qū)的碳匯總量無法滿足碳匯需求?；嵯禂?shù)計算表明，2000、2010、2022、2023年的基尼系數(shù)分別為0.76、0.83、0.85、0.88，均超過0.60，處于不協(xié)調(diào)狀態(tài)，說明研究區(qū)碳補(bǔ)償率存在較大的空間分異，并呈現(xiàn)出不斷擴(kuò)大的趨勢?？傊?，研究區(qū)在減少碳排放方面應(yīng)進(jìn)一步加強(qiáng)管理。

圖11980—2023年研究區(qū)碳補(bǔ)償率變化

2.2 碳排放風(fēng)險

參考湯峰等[13的自然斷點法，結(jié)合研究區(qū)實際情況將碳排放風(fēng)險劃分為4個等級： -6.51～0.11 、0.12～12.26、12.27～27.05、27.06以上，分別代表低度風(fēng)險、中度風(fēng)險、高度風(fēng)險和重度風(fēng)險。

由表3可知，2010—2023年，除B區(qū)、C區(qū)、K區(qū)和L區(qū)外，各地市土地利用碳排放風(fēng)險值均呈不同程度的增加。其中，F(xiàn)區(qū)、G區(qū)和M區(qū)碳排放風(fēng)險值增速較快，2000—2023年F區(qū)碳排放風(fēng)險值由32.02增加至85.06，始終處于重度風(fēng)險狀態(tài)；G區(qū)由0.11增加至11.04，由低度風(fēng)險轉(zhuǎn)變?yōu)橹卸蕊L(fēng)險，其碳排放風(fēng)險轉(zhuǎn)變與耕地和建設(shè)面積增加以及林地面積下降有關(guān)；M區(qū)由4.41增加至15.37，由中度風(fēng)險變?yōu)楦叨蕊L(fēng)險。G區(qū)和M區(qū)在2000—2010年增速較慢，在2010—2023年增速較快，預(yù)測未來碳排放風(fēng)險仍有增加趨勢。N區(qū)和L區(qū)的碳排放風(fēng)險值始終處于較低水平，但風(fēng)險值逐年升高。A區(qū)、D區(qū)、E區(qū)、I區(qū)和J區(qū)的碳排放風(fēng)險值逐年升高，在2000—2010年增速較快，但在2010—2023年增速變緩，預(yù)測未來會有降低的趨勢。

2000年，各地市碳排放風(fēng)險等級空間差異較不明顯，N區(qū)、L區(qū)及G區(qū)為低度風(fēng)險地區(qū)，中度風(fēng)險地區(qū)7個，占比 53.85% ，B區(qū)和A區(qū)為高度風(fēng)險地區(qū)，F(xiàn)區(qū)為重度風(fēng)險地區(qū)。2010年，空間差異較為明顯，F(xiàn)區(qū)仍為重度風(fēng)險地區(qū)，中度風(fēng)險地區(qū)由7個降為5個，其中G區(qū)由低度風(fēng)險升為中度風(fēng)險;高度風(fēng)險地區(qū)由2個升至5個，其中J區(qū)、C區(qū)和E區(qū)由中度風(fēng)險升為高度風(fēng)險；低度風(fēng)險區(qū)2個，分別為N區(qū)和L區(qū)。2020年，各區(qū)域空間差異顯著，重度風(fēng)險區(qū)2個，分別為F區(qū)和A區(qū)，其中F區(qū)始終為重度風(fēng)險區(qū)；高度風(fēng)險區(qū)5個，高度風(fēng)險和重度風(fēng)險區(qū)占比 53.85% ;N區(qū)始終為低度風(fēng)險區(qū)，主要是由于該區(qū)生態(tài)環(huán)境較好，林地為主要的用地類型，碳吸收能力強(qiáng)，因此碳排放風(fēng)險等級低。2023年，各區(qū)域空間差異明顯，重度風(fēng)險區(qū)2個，分別為F區(qū)和A區(qū)，其中F區(qū)始終為重度風(fēng)險區(qū)；高度風(fēng)險區(qū)7個，高度風(fēng)險和重度風(fēng)險地區(qū)占比較高；N區(qū)始終為低度風(fēng)險區(qū)?？傮w來看，2000—2023年各地市碳排放風(fēng)險等級逐年升高，空間差異性逐漸明顯。

表32000一2023年研究區(qū)各地市土地利用碳排放風(fēng)險值

2.3 碳排放強(qiáng)度

由圖2可知，研究區(qū)碳排放強(qiáng)度呈逐年下降趨勢，其中1980—2000年碳排放強(qiáng)度降速較快，2000—2005年間碳排放強(qiáng)度有所回升，2005年開始平穩(wěn)下降，之后維持在較低碳排放強(qiáng)度水平。

圖21980—2023年研究區(qū)碳排放強(qiáng)度變化

由表4可知，2000年，研究區(qū)各區(qū)域大多處于較高的碳排放強(qiáng)度水平，N區(qū)， .M 區(qū)和A區(qū)碳排放強(qiáng)度在 0.50～1.00 ，其余城市碳排放強(qiáng)度均高于1.00，占比 76.92% ，其中L區(qū)、F區(qū)、I區(qū)和K區(qū)土地利用碳排放強(qiáng)度在 1.01～2.00 ，處于較高水平，而B區(qū)、G區(qū)、D區(qū)、E區(qū)、J區(qū)及C區(qū)土地利用碳排放強(qiáng)度均大于2.00 。2010年，大部分地區(qū)碳排放強(qiáng)度呈下降趨勢，且處于較低水平；N區(qū)和A區(qū)碳排放強(qiáng)度較低，在0.5以下，而M區(qū)、D區(qū)和J區(qū)碳排放強(qiáng)度較為穩(wěn)定，仍大于2.00，L區(qū)、F區(qū)、I區(qū)及K區(qū)的碳排放強(qiáng)度由1.00～2.00 降低到 0.50～1.00 ，B區(qū)、G區(qū)、E區(qū)及C區(qū)的碳排放強(qiáng)度由2.00以上降低到 1.00～2.00 。2020年，各區(qū)域碳排放強(qiáng)度進(jìn)一步降低，其中N區(qū)、L區(qū)、A區(qū)和K區(qū)碳排放強(qiáng)度小于0.50，D區(qū)碳排放強(qiáng)度在 1.00～2.00 ，而G區(qū)碳排放強(qiáng)度有所升高。2023年僅F區(qū)發(fā)生了變化，碳排放強(qiáng)度降低，其他城市的碳排放級別均未發(fā)生變化。

總體來看，2000—2023年各區(qū)域碳排放強(qiáng)度呈不同程度的降低，其中，L區(qū)、B區(qū)、K區(qū)從碳排放強(qiáng)度較高水平轉(zhuǎn)變?yōu)檩^低水平；N區(qū)和A區(qū)一直處于較低碳排放強(qiáng)度水平，其中N區(qū)生產(chǎn)總值較低，但單位生產(chǎn)總值所帶來的碳排放量極低，因此碳排放強(qiáng)度低，而A區(qū)碳排放量較高，但其生產(chǎn)總值較高，因此，A區(qū)碳排放強(qiáng)度相對較低；G區(qū)、D區(qū)、J區(qū)、E區(qū)在研究期內(nèi)一直處于較高碳排放強(qiáng)度水平。

表42000—2023年研究區(qū)各地市土地利用排放強(qiáng)度

2.4經(jīng)濟(jì)貢獻(xiàn)系數(shù)

由表5可知，2000—2023年各區(qū)域經(jīng)濟(jì)貢獻(xiàn)系數(shù)處于相對波動的狀態(tài)。2000年，除N區(qū)、M區(qū)和A區(qū)經(jīng)濟(jì)貢獻(xiàn)系數(shù)在2以上外，其余地區(qū)經(jīng)濟(jì)貢獻(xiàn)系數(shù)均小于2，B區(qū)、D區(qū)、G區(qū)、E區(qū)、J區(qū)及C區(qū)經(jīng)濟(jì)貢獻(xiàn)系數(shù)在 0～1.00 ，水平較低；L區(qū)、I區(qū)、F區(qū)和K區(qū)的經(jīng)濟(jì)貢獻(xiàn)系數(shù)在 1.00～2.00 ：2010年，B區(qū)、G區(qū)、A區(qū)、K區(qū)和M區(qū)經(jīng)濟(jì)貢獻(xiàn)水平有所提高，其中A區(qū)經(jīng)濟(jì)貢獻(xiàn)率提升最快，達(dá)到4.00以上；2020年，A區(qū)、D區(qū)、E區(qū)、F區(qū)、G區(qū)、J區(qū)、M區(qū)經(jīng)濟(jì)貢獻(xiàn)水平下降，B區(qū)、C區(qū)、I區(qū)、K區(qū)和N區(qū)經(jīng)濟(jì)水平提高。2023年研究區(qū)未發(fā)生明顯變化。整體來看，N區(qū)以及A區(qū)在研究期處于較高經(jīng)濟(jì)貢獻(xiàn)水平，N區(qū)生產(chǎn)總值較低，主要依靠較低碳排放量以及較高碳匯量維持較高經(jīng)濟(jì)貢獻(xiàn)率，A區(qū)碳排放量高，但生產(chǎn)總值遠(yuǎn)高于其他地區(qū)，因此能維持較高的經(jīng)濟(jì)貢獻(xiàn)水平。F區(qū)是重要的石油生產(chǎn)和石化工業(yè)基地，其雖地區(qū)生產(chǎn)總值較高，但由于能源消費帶來的碳排放量遠(yuǎn)高于其他地區(qū)，F(xiàn)區(qū)經(jīng)濟(jì)貢獻(xiàn)水平較低。

表52000一2023年研究區(qū)各地市土地利用經(jīng)濟(jì)貢獻(xiàn)系數(shù)

2.5生態(tài)承載系數(shù)

由表6可知，研究區(qū)生態(tài)承載系數(shù)呈現(xiàn)北部地區(qū)偏高，南部地區(qū)偏低的態(tài)勢。2000年，N區(qū)、L區(qū)以及G區(qū)生態(tài)承載系數(shù)較高，其余地區(qū)均低于5；A區(qū)、E區(qū)和K區(qū)生態(tài)承載系數(shù)在 1.00～5.00;B 區(qū)、F區(qū)、M區(qū)、D區(qū)、I區(qū)、J區(qū)和C區(qū)生態(tài)承載系數(shù)在 0～ 1.00。2010年，L區(qū)、I區(qū)生態(tài)承載系數(shù)明顯上升，E區(qū)生態(tài)承載系數(shù)有所下降，其余地區(qū)較為穩(wěn)定。2020年，L區(qū)、G區(qū)、I區(qū)生態(tài)承載系數(shù)均下降；K區(qū)生態(tài)承載系數(shù)上升，其余地區(qū)較為穩(wěn)定。2023年，各區(qū)域生態(tài)承載系數(shù)均未發(fā)生明顯改變。整體來看，2000—2023年，N區(qū)處于較高生態(tài)承載水平，與其林地資源豐富，碳吸收能力強(qiáng)有關(guān)；而B區(qū)、F區(qū)、M區(qū)、D區(qū)、I區(qū)、E區(qū)、J區(qū)、C區(qū)常年處于低生態(tài)承載水平，因此需加強(qiáng)其碳匯能力，控制碳排放量，以提高其生態(tài)承載能力。

表62000一2023年研究區(qū)各地市土地利用生態(tài)承載系數(shù)

2.6碳平衡分區(qū)

結(jié)合各區(qū)域碳排放基尼系數(shù)、碳排放強(qiáng)度、經(jīng)濟(jì)貢獻(xiàn)及生態(tài)貢獻(xiàn)等指標(biāo)計算結(jié)果，考量研究區(qū)實際情況，研究區(qū)碳平衡分為3類：低碳優(yōu)化區(qū)、碳總量控制區(qū)、碳匯功能區(qū)。由表7可知，低碳優(yōu)化區(qū)包括A區(qū)，該地區(qū)碳排放風(fēng)險在研究期內(nèi)始終呈中高風(fēng)險，且碳排放風(fēng)險有繼續(xù)增高趨勢，但該地區(qū)由于經(jīng)濟(jì)發(fā)展水平較高，因此碳排放強(qiáng)度較低；該地區(qū)經(jīng)濟(jì)貢獻(xiàn)系數(shù)始終處于較高水平，位于各區(qū)域前列，而生態(tài)承載系數(shù)始終處于較低水平。碳總量控制區(qū)包括B區(qū)、F區(qū)、M區(qū)、D區(qū)、I區(qū)、E區(qū)、J區(qū)、C區(qū)，該區(qū)域能源型城市較多，碳排放風(fēng)險基本處于中、高風(fēng)險水平，碳排放強(qiáng)度水平基本居研究區(qū)前列，經(jīng)濟(jì)貢獻(xiàn)系數(shù)基本處于較低水平，表明該區(qū)域碳排放的經(jīng)濟(jì)效益較低；生態(tài)承載系數(shù)處于中等及較低水平。碳匯功能區(qū)包括N區(qū)、L區(qū)、G區(qū)、K區(qū)，該區(qū)域碳匯功能用地較多，如森林資源豐富，因此該區(qū)域碳排放風(fēng)險在研究期基本呈中、低風(fēng)險水平；碳排放強(qiáng)度因碳匯功能較高而處于較低水平，經(jīng)濟(jì)貢獻(xiàn)系數(shù)也因碳吸收能力較好而處于中等偏上水平；該類區(qū)域生態(tài)承載系數(shù)較高。

表7研究區(qū)土地利用碳平衡分區(qū)

3討論

本研究將研究區(qū)碳平衡分為低碳優(yōu)化區(qū)、碳總量控制區(qū)和碳匯功能區(qū)。針對不同區(qū)域特點提出針對性碳減排策略。

3.1土地利用低碳優(yōu)化區(qū)減排策略

低碳優(yōu)化區(qū)包括A區(qū)，盡管該區(qū)的土地利用碳排放量、碳排放風(fēng)險較高，但其良好的經(jīng)濟(jì)基礎(chǔ)和發(fā)達(dá)的重工業(yè)為土地低碳利用綠色轉(zhuǎn)型提供了契機(jī)。因此，可依托其經(jīng)濟(jì)和技術(shù)實力，引導(dǎo)產(chǎn)業(yè)向高科技、新材料和清潔能源領(lǐng)域轉(zhuǎn)型，強(qiáng)化服務(wù)業(yè)競爭力。在確保自身經(jīng)濟(jì)增長的同時，憑借其影響力和經(jīng)濟(jì)輻射力，激發(fā)周邊區(qū)域的經(jīng)濟(jì)與產(chǎn)業(yè)結(jié)構(gòu)協(xié)同優(yōu)化，從而增強(qiáng)碳排放效益的持久性。此外，可逐步拓展城市土地建設(shè)以促進(jìn)工業(yè)化和城市化發(fā)展。在低碳優(yōu)化區(qū)內(nèi)，主要任務(wù)是強(qiáng)化工業(yè)品和服務(wù)產(chǎn)品的供應(yīng)，鼓勵生產(chǎn)要素聚集，適當(dāng)增加重點發(fā)展區(qū)域?qū)Ω叨酥圃鞓I(yè)和現(xiàn)代服務(wù)業(yè)的空間支持。發(fā)展過程中水資源、生活能源等資源消耗上升，對環(huán)境造成較大壓力，因此，需科學(xué)規(guī)劃居民居住區(qū)，適當(dāng)增加綠地面積，開展城市綠地、公園建設(shè)，以提高城市的生態(tài)承載力。

3.2土地利用碳總量控制區(qū)減排策略

碳總量控制區(qū)包括B區(qū)、F區(qū)、M區(qū)、D區(qū)、I區(qū)、E區(qū)、J區(qū)、C區(qū)。此類區(qū)域固碳能力相對較弱，經(jīng)濟(jì)發(fā)展優(yōu)勢不足。這些地區(qū)大多為能源密集型城市，面對高能耗帶來的高碳排放問題，應(yīng)加大科技創(chuàng)新投入，提高能源利用效率，推動產(chǎn)業(yè)向低碳、高效方向轉(zhuǎn)型。碳總量控制區(qū)重點在于優(yōu)化碳排放強(qiáng)度，提升經(jīng)濟(jì)效能。為此，應(yīng)科學(xué)配置工業(yè)和居民用地，確保其與產(chǎn)業(yè)規(guī)劃相協(xié)調(diào)，推動化石能源的高效再利用，避免生產(chǎn)過程中的資源損耗。同時大幅提升清潔能源的使用比例，推廣低碳生產(chǎn)工藝，降低資源消耗、污染及碳排放強(qiáng)度。此外，發(fā)展循環(huán)經(jīng)濟(jì)理念，普及低碳技術(shù)，以降低二氧化碳的排放，且生活能源的循環(huán)利用也應(yīng)得到重視。在追求經(jīng)濟(jì)增長的同時，城市中的碳匯資源也不容忽視，需堅決保護(hù)農(nóng)田以及濕地、林地等生態(tài)環(huán)境用地。

3.3土地利用碳匯功能區(qū)減排策略

碳匯功能區(qū)包括N區(qū)、L區(qū)、G區(qū)、K區(qū)。此類區(qū)域表現(xiàn)出較低的碳排放水平和風(fēng)險，其碳匯潛力較大，生態(tài)環(huán)境承受的壓力相對較小。此外，該區(qū)域森林資源豐富，是主要的林業(yè)基地；擁有優(yōu)質(zhì)的生態(tài)環(huán)境，有利于發(fā)展旅游業(yè)和高新技術(shù)產(chǎn)業(yè)，應(yīng)重視維護(hù)碳匯資源，保護(hù)這些關(guān)鍵的碳匯區(qū)域，持續(xù)促進(jìn)綠色產(chǎn)業(yè)可持續(xù)發(fā)展。碳匯功能區(qū)在水源保護(hù)、土壤穩(wěn)定、防風(fēng)固沙及生物多樣性保護(hù)等生態(tài)服務(wù)中發(fā)揮重要作用。在確保經(jīng)濟(jì)持續(xù)發(fā)展的同時，應(yīng)注重環(huán)境保護(hù)和修復(fù)工作。作為維持地區(qū)生態(tài)系統(tǒng)平衡的關(guān)鍵地帶，嚴(yán)格管理建設(shè)用地面積的增長，提供生態(tài)服務(wù)。同時，應(yīng)強(qiáng)化生態(tài)修復(fù)措施，保護(hù)并復(fù)原自然環(huán)境。充分利用當(dāng)?shù)鬲毺氐淖匀毁Y源和地理優(yōu)勢，推動小規(guī)模、高附加值的特色產(chǎn)業(yè)發(fā)展，以實現(xiàn)經(jīng)濟(jì)、社會與生態(tài)的可持續(xù)共生。

4結(jié)論

本研究綜合分析了研究區(qū)的碳排放基尼系數(shù)、碳排放風(fēng)險、碳排放強(qiáng)度、經(jīng)濟(jì)貢獻(xiàn)系數(shù)及生態(tài)承載系數(shù)，并以此為分區(qū)指標(biāo)對其進(jìn)行碳平衡分區(qū)，提出不同區(qū)域差異化減排策略，主要結(jié)論如下。

（1）研究區(qū)的碳補(bǔ)償率整體上處于低位，并且隨時間的推移逐年下降，說明研究區(qū)碳匯總量無法完全滿足其碳匯需求。2000—2023年基尼系數(shù)均超過0.60，處于不協(xié)調(diào)狀態(tài)，研究區(qū)的碳補(bǔ)償率存在較大的空間分異，并呈現(xiàn)出不斷擴(kuò)大的趨勢。

（2）2000—2023年，各地區(qū)土地利用碳排放風(fēng)險呈不同程度的增加，其中F區(qū)、G區(qū)和M區(qū)碳排放風(fēng)險值增速較快，B區(qū)、C區(qū)和K區(qū)碳排放風(fēng)險在2010—2023年有所下降。N區(qū)雖始終處于低度風(fēng)險，但風(fēng)險逐年升高。A區(qū)、D區(qū)、E區(qū)、I區(qū)和J區(qū)雖碳排放風(fēng)險逐年升高，但在2010一2023年增速變緩，預(yù)測未來有降低的趨勢。G區(qū)和M區(qū)碳排放風(fēng)險逐年遞增，未來仍有上升趨勢。在2000年，各地區(qū)碳排放風(fēng)險等級差異不明顯，在2010年，空間差異較明顯，在2020年和2023年，各地市空間差異十分明顯，總體來看，2000—2023年各地市碳排放風(fēng)險等級逐年升高，空間差異性逐漸明顯。

（3）1980—2023年，研究區(qū)碳排放強(qiáng)度呈逐年下降趨勢，其中1980—2000年間降速最快，2000—2005年有所回升，而后2005年開始平穩(wěn)下降，維持在較低碳排放強(qiáng)度水平?？傮w來看，2000—2023年研究區(qū)各地區(qū)碳排放強(qiáng)度呈現(xiàn)不同程度的降低；N區(qū)和A區(qū)一直處于低碳排放強(qiáng)度水平，G區(qū)、D區(qū)、J區(qū)一直處于較高碳排放強(qiáng)度水平。2000—2023年，研究區(qū)各地區(qū)經(jīng)濟(jì)貢獻(xiàn)值處于相對波動的狀態(tài)。研究期間N區(qū)以及A區(qū)處于較高經(jīng)濟(jì)貢獻(xiàn)水平，其余地區(qū)經(jīng)濟(jì)貢獻(xiàn)水平均在較低水平。研究區(qū)生態(tài)承載系數(shù)呈現(xiàn)北部地區(qū)偏高，南部地區(qū)偏低的態(tài)勢，整體呈偏低態(tài)勢。

（4）研究區(qū)碳平衡分為低碳優(yōu)化區(qū)、碳總量控制區(qū)和碳匯功能區(qū)。實施分類施策：低碳優(yōu)化區(qū)依托經(jīng)濟(jì)與技術(shù)優(yōu)勢，推動產(chǎn)業(yè)綠色轉(zhuǎn)型和服務(wù)業(yè)升級，通過經(jīng)濟(jì)輻射帶動區(qū)域協(xié)同發(fā)展，同時優(yōu)化城市空間布局提升生態(tài)承載力；碳總量控制區(qū)重點加強(qiáng)科技創(chuàng)新和能效提升，推進(jìn)產(chǎn)業(yè)低碳化改造和清潔能源替代，優(yōu)化土地利用結(jié)構(gòu)，發(fā)展循環(huán)經(jīng)濟(jì)并嚴(yán)格保護(hù)生態(tài)空間；碳匯功能區(qū)則強(qiáng)化森林、濕地等核心碳匯區(qū)保護(hù)，嚴(yán)格控制開發(fā)強(qiáng)度，重點發(fā)展生態(tài)旅游、有機(jī)農(nóng)業(yè)等環(huán)境友好型產(chǎn)業(yè)，同步推進(jìn)生態(tài)修復(fù)工程。通過精準(zhǔn)施策，構(gòu)建經(jīng)濟(jì)發(fā)展與碳減排協(xié)同推進(jìn)的區(qū)域可持續(xù)發(fā)展格局。

參考文獻(xiàn)

[1]肖磊，孫貴艷.基于土地利用變化的重慶碳排放時空格局研究[J].國土與自然資源研究，2023（3）：28-32.

[2] HEIL M T，WODON Q T. Inequality in CO₂ emissionsbetweenpoor and rich countries[J]. The journal ofenvironmentamp;development，1997，6（4）：426-452.

[3]馮杰，張勝，王濤.中國省際土地利用碳排放及其影響因素分析[J].統(tǒng)計與決策，2019，35（5）：141-145.

[4]鄭欣，程久苗，鄭碩.基于土地利用結(jié)構(gòu)變化的蕪湖市碳排放及其影響因素研究[J].水土保持研究，2012，19（3）：259-262，268.

[5]魏媛，簡小玉.基于碳達(dá)峰碳中和的貴州省土地利用碳排放演變及預(yù)測研究[J].生態(tài)經(jīng)濟(jì)，2022（4）：108-114.

[6]陳景.石家莊西部太行山區(qū)土地利用變化碳排放演變及預(yù)測研究[D].北京：中國地質(zhì)大學(xué)（北京），2015.

[7]魏燕茹，陳松林.福建省土地利用碳排放空間關(guān)聯(lián)性與碳平衡分區(qū)[J].生態(tài)學(xué)報，2021，41（14）：5814-5824.

[8]楊靜媛，張明，多玲花，等.江西省土地利用碳排放空間格局及碳平衡分區(qū)[J].環(huán)境科學(xué)研究，2022，35（10）：2312-2321.

[9]張正峰，張棟.基于社會網(wǎng)絡(luò)分析的京津冀地區(qū)碳排放空間關(guān)聯(lián)與碳平衡分區(qū)[J].中國環(huán)境科學(xué)，2023，43（4）：2057-2068.

[10]劉佳寧，孫麗娜.1980—2020年黑龍江省土地利用碳排放時空分異研究[J].中國農(nóng)學(xué)通報，2023，39（26）：86-97.

[11]張樂勤.基于基尼系數(shù)方法的長三角區(qū)域綠色發(fā)展?fàn)顩r均衡性探析[J].環(huán)境與可持續(xù)發(fā)展，2021，46（1）：102-108.

[12]張杰，陳海，劉迪，等.基于縣域尺度土地利用碳排放的時空分異及影響因素研究[J].西北大學(xué)學(xué)報（自然科學(xué)版），2022，52（1）：21-31.

[13]湯峰，殷彥玲，張蓬濤，等.基于土地利用變化的高陽縣生態(tài)系統(tǒng)服務(wù)價值和碳排放風(fēng)險研究[J].水土保持研究，2019，26（4）：265-271.

[14]趙榮欽，張帥，黃賢金，等.中原經(jīng)濟(jì)區(qū)縣域碳收支空間分異及碳平衡分區(qū)[J].地理學(xué)報，2014，69（10）：1425-1437.

（責(zé)任編輯：胡立萍）