煤炭开发过程碳排放特征及碳中和发展的技术途径

任世华 谢亚辰 焦小淼 谢和平

任世华, 谢亚辰, 焦小淼, 等. 煤炭开发过程碳排放特征及碳中和发展的技术途径 [J]. 工程科学与技术, 2022, 54(1): 60-68. doi: 10.15961/j.jsuese.202100924
引用本文: 任世华, 谢亚辰, 焦小淼, 等. 煤炭开发过程碳排放特征及碳中和发展的技术途径 [J]. 工程科学与技术, 2022, 54(1): 60-68. doi: 10.15961/j.jsuese.202100924
REN Shihua, XIE Yachen, JIAO Xiaomiao, et al. Characteristics of Carbon Emissions During Coal Development and Technical Approaches for Carbon Neutral Development [J]. Advanced Engineering Sciences, 2022, 54(1): 60-68. doi: 10.15961/j.jsuese.202100924
Citation: REN Shihua, XIE Yachen, JIAO Xiaomiao, et al. Characteristics of Carbon Emissions During Coal Development and Technical Approaches for Carbon Neutral Development [J]. Advanced Engineering Sciences, 2022, 54(1): 60-68. doi: 10.15961/j.jsuese.202100924

煤炭开发过程碳排放特征及碳中和发展的技术途径

基金项目: 中国工程院咨询研究项目(2022–28–12)
详细信息
    • 收稿日期:  2021-09-13
    • 网络出版时间:  2021-12-29 12:46:00
  • 作者简介:

    任世华(1980—),男,研究员.研究方向:煤炭相关技术经济评价、能源战略. E-mail:ren@cct.org.cn
    任世华,煤炭科学研究总院科技支持中心主任,研究员,主要从事洁净煤技术研究、技术经济评价等工作。曾主持或参与国家科技重大专项、国家863项目、国家科技支撑、自然科学基金、国家发改委重点课题、国家能源局重点课题等50余项科研课题的研究。曾获省部级科学技术奖特等奖1项、一等奖2项、二等奖9项、三等奖13项,2014年获全国煤炭青年科学技术奖。发表论文30余篇,参与出版著作11部。

  • 中图分类号: TD8

Characteristics of Carbon Emissions During Coal Development and Technical Approaches for Carbon Neutral Development

  • 摘要: 理清煤炭开发过程碳排放特征,是推动煤炭开发过程碳达峰、碳中和的前提和基础。基于煤炭开发全生命周期碳排放清单分析方法,重点从生产用能、瓦斯排放及矿后活动3个环节,建立煤炭开发过程碳排放计算模型,测算煤炭开发过程碳排放量,并分析不同环节碳排放特征,提出煤炭开发过程碳减排技术途径。研究表明:1)生产用能碳排放强度呈下降趋势,由2010年81.5 kg/t,先快速降低到2015年的66.5 kg/t,而后缓慢降低到2020年的65.4 kg/t;2)瓦斯排放(碳排放)强度呈快速下降趋势,由2010年的123.7 kg/t,降低到2020年的67.6 kg/t;3)矿后活动碳排放强度呈缓慢降低趋势,由2010年的21.5 kg/t,逐渐降至2020年的18.0 kg/t。基于碳排放量和排放特征分析,提出了煤炭开发过程实现碳中和的5大技术途径,即大力推广应用煤炭开发节能提效技术,减少煤炭开发过程能源消耗带来的碳排放;持续攻关煤矿瓦斯抽采利用技术,减少煤炭开发过程甲烷排放;加快探索煤炭开发新技术,从原理上改变碳排放特性;有序研发和示范煤矿区煤与新能源耦合利用技术,降低单位产品碳排放强度;尽早谋划煤矿区二氧化碳捕集、利用与固化、封存技术,形成煤炭行业独有的碳中和实现途径。

     

    Abstract: Clarifying the carbon emission characteristics is the prerequisite and basis for promoting the carbon peak and carbon neutrality in the coal development process. Based on the carbon emission inventory analysis method of the whole life cycle of coal development, focusing on the three links of production energy, gas emission and post-mining activities, the carbon emission calculation model of the coal development process was established, and the carbon emission during the coal development process was measured. The characteristics of carbon emissions in different links was analyzed, and the technical approaches for carbon emission reduction in the coal development process was proposed. The results showed that the carbon emission intensity of energy used in production showed a fluctuating downward trend, from 81.5 kg/t in 2010, first reduced to 66.5 kg/t in 2015, and then slowly reduced to 65.4 kg/t in 2020. The carbon emission intensity of coal mine gas emissions (carbon emissions) shown a rapid downward trend, from 123.7 kg/t in 2010 to 67.6 kg/t in 2020. The carbon emission intensity of post-mining activities was slowly decreasing, from 21.5 kg/t in 2010 to 18.0 kg/t in 2020. Based on the above researches, five technical approaches for carbon neutral achievement in the coal development process were proposed, namely, vigorously promoting the application of coal development energy-saving and efficiency-increasing technologies to reduce carbon emissions caused by coal development energy consumption, continuously tackling key coal mine gas extraction and utilization technologies to reduce methane emissions during coal development, accelerating the exploration of disruptive technologies for coal development to change carbon emission characteristics in principle, developing and demonstrating coal and new energy coupling utilization technologies in coal mines to reduce the carbon emission coefficient per unit product, planning the way to realize the technology of carbon dioxide capture, utilization, solidification, and storage in coal mine areas to form a unique way for achieving carbon neutrality of the coal industry.

     

  • 中国政府提出“力争2030年前实现二氧化碳达峰、2060年前实现碳中和”的“双碳”目标,将推动一场广泛而深刻的经济社会系统性变革[1]。煤炭开发利用过程中产生的碳排放是中国碳排放的主要来源,约占全国碳排放总量的60%~70%。从煤炭开发和利用过程看,煤炭利用过程的碳排放量占比近90%[2],煤炭开发过程的碳排放量占比约10%[3]。煤炭开发过程的碳排放总量虽不多,但甲烷(CH4)排放量占能源活动CH4总排放量的80%以上,约占中国CH4总排放量的1/3[4],因CH4较CO2更为显著的增温效益,已引起广泛的关注[5-7]

    一些学者对煤炭开发过程CO2排放及CH4排放进行研究,取得了一系列的研究成果[8-9]。王宁[10]和李学武[11]等采用原煤产量法,以中国煤矿区CH4排放为例建立“产量–涌出量模型”,并结合矿井实际调研数据,得出区域CH4涌出系数。马翠梅等[12]对中国2010—2016年煤炭开采和矿后活动的CH4排放量进行估算,提出加强煤矿区CH4的回收利用。Zhu等[13]更新分省煤矿区CH4排放因子,对中国2006—2010年煤矿区CH4排放进行估算,预测了2011—2020年中国煤矿CH4排放量。张媛[14]以晋城长河流域的11个煤矿为研究区,更新矿区煤矿开采CO2排放强度和CH4排放强度,得出矿区开采吨煤炭排放量为44.5 kg。于胜民等[3]以中国7 575处井工煤矿CO2相对涌出量测定数据为基础,统计得井工煤矿开采过程CO2逃逸排放因子为每t原煤逃逸CO26.57 m3,不同瓦斯等级矿井之间均值差异显著,低瓦斯矿井、高瓦斯矿井、突出矿井均值分别为每吨原煤产CO25.41、9.38和10.50 m3

    煤炭开发过程中排放的CH4,既是温室气体,也是低碳的能源,为充分利用煤矿区CH4资源,近年来中国加大了对地面煤层气和井下瓦斯的抽采利用力度。依据中国煤炭工业发展报告及相关文献[12,15],中国煤层气抽采总量由2005年的2.33×109 m3快速增加到2020年的2.05×1010 m3,年均增加17.0%;对应的利用量由2005年的7.4×108m3快速增加到2020年的1.288×1010 m3,年均增加22.7%。煤层气的总体利用率则由2006年的31.6%增加到2020年的62.6%,年均增长2.1%(图1)。

    图  1  煤层气(煤矿瓦斯)抽采量及利用率
    Fig.  1  Coalbed methane (coal mine gas) extraction volume and utilization rate
    下载: 全尺寸图片

    “双碳”目标下,煤炭开发过程的低碳清洁化是必然要求[16-17]。理清煤炭开发过程碳排放量和排放特征对于寻找可持续利用煤炭资源方法和途径至关重要[18],碳核算及碳排放清单的建立是挖掘减排潜力的前提和基础[19]。然而,近年来多数学者在核算煤炭开发过程碳排放时,对煤矿区CH4的回收利用考虑不充分,使得对煤矿区CH4排放的估算偏高[13,20-22]。鉴于此,本文在项目组前期研究构建的煤炭开发利用碳排放清单[23-26]的基础上,对中国煤炭开发过程碳排放量进行估算,分析煤炭开发过程不同环节的碳排放特征,进而提出中国煤炭开发过程碳中和发展的技术途径。

    煤炭开发过程是指由井工或露天煤矿开采出原煤,并经洗选成为煤炭产品的过程,煤炭开发过程碳排放计算模型的范围、边界及输入输出如图2所示。该计算模型的输入主要包括煤炭开发过程原煤、油、气及电力、热力等能源的消耗量,而模型输出主要包含井工开采或露天开采主要环节的生产能耗、瓦斯排放和矿后活动的碳排放量及排放强度。模型输出还可根据需要输出井工开采单产品能耗(煤、油、气、电力)CO2排放量、井工开采单产品CH4排放量、井工开采单产品矿后活动CH4排放量、露天开采生产用能单产品能耗(煤、油、气、电力)CO2排放量及露天开采单产品CH4排放量。

    图  2  煤炭开发过程碳排放计算模型范围、边界及输入输出[24]
    Fig.  2  Scope, boundary and input and output of the carbon emission calculation model during the coal development process[24]
    下载: 全尺寸图片

    依据《中国煤炭生产企业温室气体排放核算方法与报告指南(试行)》[27],将煤炭开发过程中温室气体(GHG)排放分为3个环节:生产用能碳排放、瓦斯排放(碳排放)及矿后活动碳排放。生产用能碳排放主要包含输入化石燃料燃烧CO2排放量、火炬燃烧CO2排放量及净购入电力和热力隐含CO2排放量;瓦斯排放(碳排放)主要包含井工和露天开采前、开采中CH4逃逸量(折算为CO2排放量);矿后活动碳排放主要包含露天开采、废弃矿井,以及原煤在运输、洗选、储存过程CH4(折算为CO2排放量)和CO2的逃逸排放量。煤炭开发过程碳排放总量EGHG为:

    $$ {E}_{\text{GHG}}\text={E}_{瓦斯}\text+{E}_{矿后}\text+{E}_{用能} $$ (1)

    式中: ${E_{{\text{GHG}}}}$ 为煤炭开发过程碳排放总量,t; $ {E}_{用能} $ $ {E}_{瓦斯} $ $ {E}_{矿后} $ 分别为煤炭开发过程生产用能碳排放量、瓦斯排放(碳排放)量及矿后活动碳排放量,t,E用能E瓦斯E矿后表达式如下:

    $$ {E}_{用能}\text=\sum\limits_{{i=}1}^{4}{E}_{\text{C}{\text{O}}_{\text{2}}}^{{i}}+{E}_{\text{C}{\text{O}}_{\text{2}}}^{\text{p}}+{E}_{\text{C}{\text{O}}_{\text{2}}}^{\text{h}} $$ (2)

    式中: $E_{{{\rm{C}}}{{{\rm{O}}}_{2}}}^{{\rm{p}}}$ 为煤炭开发过程输入电力隐含碳排放量,t; $E_{{{\rm{C}}}{{{\rm{O}}}_{2}}}^{{\rm{h}}}$ 为煤炭开发过程输入热力隐含碳排放量,t; $E_{{{\rm{C}}}{{{\rm{O}}}_{2}}}^{i}$ 为煤炭开发过程化石燃料燃烧碳排放量,t,其中, $ {i} $ =1、2、3、4,分别表示煤炭开发过程输入煤、油、气燃烧及煤矿瓦斯火炬燃烧碳排放量; $E_{{{\rm{C}}}{{{\rm{O}}}_{2}}}^{i} $ $E_{{{\rm{C}}}{{{\rm{O}}}_{2}}}^{{\rm{p}}}$ 的表达式如下:

    $$ E_{{{\rm{C}}}{{{\rm{O}}}_{2}}}^{i} = {M_{i}} \times {\rm{E}}{{\rm{F}}_{i}} = {M_{i}} \times {Q_{i}} \times {C_{i}} \times {\beta _{i}} $$ (3)

    式中: $ {M_{i}} $ 为消耗化石燃料的质量,t; ${\rm{E}}{{\rm{F}}_{i}}$ 为消耗化石燃料对应碳排放因子,t/t; $ {Q_{i}} $ 为化石燃料的低位发热量,TJ/t; $ {C_{i}} $ 为化石燃料燃烧产生单位热量对应的碳排放量,t/TJ; $ {\;\beta _{i}} $ 为化石燃料的碳转化率,即燃料中的碳在燃烧过程转化成CO2并排放到大气中的比例,%。

    $$ {E}_{\text{C}{\text{O}}_{\text{2}}}^{\text{p}}=\gamma \times {W}_{电} $$ (4)

    式中: $ {W}_{电} $ 为煤炭开发过程输入电力,kW·h; $ \gamma $ 为区域电网年平均供电CO2排放因子,t/MW·h;实际计算中认为 $E_{{{\rm{C}}}{{{\rm{O}}}_{2}}}^{{\rm{h}}}$ $E_{{{\rm{C}}}{{{\rm{O}}}_{2}}}^{{\rm{p}}}$ 的计算方法一样。

    $$ {{E}}_{瓦斯}\text={\rm{GW}}{{\rm{P}}}_{\text{C}{\text{H}}_{\text{4}}}\times {E}_{\text{C}{\text{H}}_{\text{4}}}^{开发} $$ (5)

    式中: ${\rm{GW}}{{\rm{P}}_{{{\rm{C}}}{{{\rm{H}}}_{4}}}}$ 为CH4的温增指数,根据IPCC第2次评估报告,100 a时间尺度内1 t CH4相当于21 t CO2的增温能力,因此GWPCH4取值21; $ {E}_{\text{C}{\text{H}}_{\text{4}}}^{开发} $ 为煤炭开发过程CH4的逃逸碳排放量,计算式如式(6)所示:

    $$ {E}_{\text{C}{\text{H}}_{\text{4}}}^{开发}\text={\rho }_{{\rm{CH}}_{4}}\times (\alpha_{\rm{p}} \times W_{原煤}/10^4 -{Q}_{\text{C}{\text{H}}_{\text{4}}}^{利用}) $$ (6)

    式中: ${\;\rho _{{{\rm{C}}}{{{\rm{H}}}_{4}}}}$ 为标况下CH4的密度,取值为7.17 t /104 N·m3 $ {\alpha _{{\rm{p}}}} $ 为煤炭开采CH4排放因子,m3/t; $ {W}_{原煤} $ 为原煤产量,t; $ {Q}_{\text{C}{\text{H}}_{\text{4}}}^{利用} $ 为瓦斯回收利用量,104 N·m3

    $$ {E}_{矿后}\text={\rm{GW{P}}}_{\text{C}{\text{H}}_{\text{4}}}\times {E}_{\text{C}{\text{H}}_{\text{4}}}^{矿后}+{E}_{\text{C}{\text{O}}_{\text{2}}}^{逃逸} $$ (7)

    其中:

    $$ {E}_{\text{C}{\text{H}}_{\text{4}}}^{矿后}\text=\rho_{{\rm{CH}}_4}\times (\alpha_0 \times W_{原煤})/10^4 $$ (8)

    式中: ${E}_{\text{C}{\text{H}}_{\text{4}}}^{矿后}$ 为煤炭开发矿后活动CH4释放量,t; $ {E}_{\text{C}{\text{O}}_{\text{2}}}^{逃逸} $ 为煤炭开发过程CO2的逃逸量,t; $ \alpha _0 $ 为矿后活动CH4排放因子,m3/t。在实际计算时,基于数据的可获得性暂没有计算火炬燃烧CO2排放量和煤炭开发过程CO2的逃逸量。煤炭开发过程不同环节的碳排放强度Fi为:

    $$ {F}_{{i}}\text={{E}_{{i}}/W}_{原煤} $$ (9)

    式中: $ {F_{i}} $ 为煤炭开发过程不同环节的碳排放强度,kg/t,其中, $ {i} $ =1、2、3,分别表示生产用能、瓦斯排放及矿后活动; $ {E_{i}} $ 为煤炭开发不同环节碳排放量,kg,其中, $ {i} $ =1、2、3,分别表示生产用能、瓦斯排放及矿后活动碳排放量。

    依据煤炭开发过程碳排放量和排放强度的估算方法[20],选取《中国统计年鉴(2011—2020)》[28]和项目组测算得出的煤炭开发过程能源消耗数据,并参考《IPCC2006国家温室气体排放清单计算指南》[29-30]和相关文献[31-33],更新了煤炭开发过程不同能源消耗的碳排放因子,对中国2010—2020年煤炭开发过程生产用能的碳排放量和排放强度进行估算,结果见图3

    图  3  生产用能碳排放量及排放强度
    Fig.  3  Carbon emissions and emission intensity of energy used in production
    下载: 全尺寸图片

    图3可知,中国煤炭开发过程生产用能碳排放由2010年2.64×108 t,先增加到2011年的2.79×108 t,随后逐渐降低到2016年的2.25×108 t,而后逐渐增加到2020年2.57×108 t。煤炭开发过程生产用能碳排放主要包括煤炭、电力及油气消耗碳排放,其中:煤炭消耗碳排放量整体上呈现先降低而后增加的趋势,由2010年的1.80×108 t,先快速降低到2016年的1.32×108 t,年均降低4.9%,而后缓慢增加到2020年的1.40×108 t,年均增加1.5%。煤炭消耗碳排放占生产用能总碳排放的比例由2010年68.2%降低到2020年的54.7%,年均降低1.4%。而电力消耗碳排放量呈现增加趋势,由2010年的0.63×108 t,增加到2020年的0.98×108 t,年均增速4.8%,使得电力消耗碳排放占生产用能总碳排放的比例由2010年的23.8%增加到2020年的38.1%,年均增加1.4%。油气消耗碳排放量占比呈现先增加,而后降低的趋势,由2010年的8.0%先增加到2015年的9.1%,随后波动下降到2020年的7.3%。生产用能碳排放受原煤产量、单位产品能源消耗强度、能源消耗碳排放强度的影响,其中原煤产量是最主要影响因素,生产用能碳排放量的变化趋势基本与原煤产量变化趋势一致。

    近年来,随着中国煤炭开发机械化水平持续提高,大型煤炭企业采煤机械化水平已高达97.1%,达到发达国家水平[34],使得煤矿生产能耗逐年降低[35],结合煤矿区“电代煤”及“气代煤”的改造升级[36],显著改变了生产用能结构,使得中国煤炭开发过程生产用能碳排放强度由2010年81.5 kg/t,先快速降低到2015年的66.5 kg/t,年均降速4.0%,而后缓慢降低到2020年的65.4 kg/t,年均降速0.3%。由图3可知:吨煤煤耗碳排放强度持续降低,由2010年的55.6 kg/t快速降低到2016年的39.3 kg/t,年均降速5.6%,而后缓慢降低到2020年的35.8 kg/t,年均降低2.3%;吨煤电耗碳排放强度整体上呈现波动增加趋势,由2010年的19.4 kg/t,波动增加到2020年的24.9 kg/t,年均增加2.6%;吨煤油气消耗碳排放强度呈现降低趋势,由2010年的6.5 kg/t,波动降低到2020年的4.8 kg/t,年均降低3.0%。

    依据煤监部门统计[37]及相关文献[9,38],结合项目组研究成果[23-26],估算了2010—2020年中国煤矿瓦斯溢出量并与其他研究结果[4,13,20-22,39]对比(图4)。

    图  4  不同机构煤矿瓦斯溢出量估算量对比
    Fig.  4  Comparison of the estimated amount of coal mine gas emissions from different institutions
    下载: 全尺寸图片

    图4可知,近年来,中国煤矿瓦斯溢出量整体上处于下降趋势。由2010年的0.19×108 t,先缓慢降低到2016年的0.13×108 t,而后处于稳定状态,2020年排放量为0.12×108 t。近年煤矿瓦斯的抽采利用率持续增加,使得瓦斯排放量的变化与原煤产量变化规律不一致。通过对比已有相关研究结果(图4)可知,不同学者对煤矿瓦斯排放的估算结果差异较大,主要原因在于估算依据和数据来源不同,本文估算的瓦斯排放量综合考虑了煤矿瓦斯的抽采利用率。

    依据式(5)和(9),估算瓦斯排放折合碳排放量及排放强度(图5),中国瓦斯排放(碳排放)量由2010年的4.01×108 t先缓慢波动增加到2012年的4.03×108 t,而后逐渐降低到2016年的2.81×108 t;随着瓦斯抽采利用率的提高,瓦斯排放(碳排放)量逐渐降低到2020年的2.65×108 t。吨煤瓦斯碳排放呈现逐渐降低的趋势,由2010年的123.7 kg/t,逐渐降低到2020年的67.6 kg/t,年均降速为5.8%,由此可知,瓦斯排放(碳排放)量的变化与煤矿瓦斯溢出量变化规律一致,主要受原煤产量和瓦斯抽采利用率的波动影响,与吨煤瓦斯排放(碳排放)强度的降低对其影响较小。

    图  5  煤炭开发碳排放量及排放强度
    Fig.  5  Carbon emissions and emission intensity of coal development
    下载: 全尺寸图片

    依据中国煤炭工业统计资料整理了中国井工煤矿和露天煤矿原煤产量数据[33],露天煤矿瓦斯含量相对较低,因此矿后活动排放因子采用《IPCC2006指南》缺省值。按照《矿井瓦斯涌出量预测方法》(AQ 1018—2006)行业标准,井工矿的矿后活动CH4排放量依据式(10)计算得出[12]

    $$ {W_{{\rm{c}}}}{\text{ = }}10.385 \times {{\rm{e}}^{{{{ - 7}{.207}} \mathord{\left/ {\vphantom {{{ - 7}{.207}} {{W_{0}}}}} \right. } {{W_{0}}}}}} $$ (10)

    式中: ${W_{{\rm{c}}}}$ 为原煤的矿后活动CH4排放量,m3/t; $ {W_{0}} $ 为煤层的原始CH4含量,m3/t,通常取值1~5 m3/t[40],本文取值3 m3/t。结合Zhu等[13]相关研究成果估算了煤炭开发矿后活动碳排放量和排放强度(图5)。

    图5可知:中国煤炭开发矿后活动碳排放量由2010年的7.0×107 t,先波动增加到2014年的7.7×107 t,《随后降低到2016年的6.4×107 t,此后增加到2020年的7.1×107 t。吨煤矿后活动碳排放强度呈现逐渐降低的趋势,由2010年的21.5 kg/t,逐渐降低到2020年的18.0 kg/t,年均降速为1.7%。矿后活动碳排放量主要受原煤产量的波动影响,而吨煤矿后活动碳排放强度的降低对其影响较小。

    图5可知:中国煤炭开发过程碳排放总量由2010年的7.35×108 t,先增加到2011年的7.59×108 t,随后逐渐降低到2016年的5.70×108 t,而后再缓慢增加到2020年的5.93×108 t;而煤炭开发过程吨煤碳排放强度呈现逐年降低的趋势,由2010年的226.7 kg/t,快速降低到2016年的169.6 kg/t,年均降速4.7%,此后缓慢再降低到2020年的151.1 kg/t,年均降速2.9%,研究结论与Zhou等[41]用全生命周期方法估算的煤炭开发过程中碳排放强度164.4 kg/t的结果相近。

    图6为煤炭开发过程碳排放结构及变化趋势。

    图  6  煤炭开发过程碳排放结构及变化趋势
    Fig.  6  Carbon emission structure and change trend during coal development
    下载: 全尺寸图片

    图6可知:煤炭开发过程煤炭消耗碳排放量占比24%左右,随着煤炭开发效率的提高而下降,由2010年的24.5%,缓慢降低到2020年的23.7%,年均降低0.1%;随着采煤机械化程度的提高,煤炭开发过程电耗碳排放量占比呈现增加趋势,由2010年的8.6%,快速增加到2020年的16.5%,年均增加0.8%;煤炭开发过程油气消耗碳排放量占比较少,约占3%,呈现先增加,而后减少的变化趋势,由2010年的2.9%,先增加到2016年的3.5%,随后逐渐减少到2020年的3.1%;煤炭开发过程瓦斯排放(碳排放)量占比整体上呈现降低趋势,占比由2010年的54.6%,降低到2020年的44.8%,年均降低1%;而矿后活动碳排放量占比呈现增加趋势,由2010年的9.5%,增加到2020年的11.9%,年均增加0.2%。本文估算的2016年煤炭开发过程煤矿瓦斯排放(碳排放)量和矿后活动碳排放量的占比为60.5%,与Wang等[20]估算瓦斯排放占总排放62%的研究结果一致。2020年煤炭开发过程煤矿瓦斯排放(碳排放)量占总排放量56.7%,随着中国采煤机械化水平的提高,煤炭开发吨煤能耗、吨煤瓦斯排空总体呈现下降趋势,使得碳排放总量及强度下降趋势明显,未来煤矿瓦斯抽采利用是煤炭开发过程碳减排的最核心内容。

    当前的煤炭开发方式不可避免地消耗能源和带来碳排放,依靠其自身优化,可在一定程度上减少碳排放,但难以实现碳中和要求。通过技术创新,推进煤炭开发过程节能提效,提高煤矿瓦斯抽采利用率,探索低碳型煤炭开发新工艺技术,推动煤与新能源耦合利用,布局煤矿区二氧化碳捕集、利用与固化及封存技术,是煤炭开发过程实现碳中和的必然要求。

    煤炭开发过程机械设备运行等消耗的煤、油、气及电力、热力等生产用能产生的CO2排放占总排放量的20%左右。因此,大力推广应用节能提效技术,提高煤炭开发过程中能源利用效率,减少能源用量是降低碳排放的优先途径[35]

    推广基于节能降碳的煤炭开采优化设计技术,在确保安全的条件下,严格按生产规模优化配置装备和能力,减少“大马拉小车”的能源浪费;研发应用智能变频永磁驱动等技术,提高掘进机、采煤机等大型矿用设备能源利用效率,减少能源用量;加快研发应用煤矿智能化和矿山物联网技术,攻克自适应割煤、煤岩识别、超前支护自动化、智能放煤、装备智能定位及路径规划等技术难题,减少不必要的功率损失和能源消耗;全面应用余热、余压、节水、节材等综合利用技术,以能源和材料的再利用间接减少能源消耗等。

    近年来,瓦斯抽采利用率不断提高,但抽采瓦斯直接排空的比例依然在50%左右,使得煤矿瓦斯造成的碳排放占2020年煤炭开发过程碳排放的56.7%。因此,持续加强煤矿瓦斯抽采利用技术的研发和应用,进一步提高煤矿瓦斯抽采率和利用率,是降低煤炭开发过程碳排放量的必由之路。

    在国家科技重大专项“大型油气田及煤层气”资助的系列技术攻关引领下,集聚优势科技资源,加强协同,持续攻关低渗煤层抽采关键工艺技术,提高低渗煤层的煤层气渗出效率,解决煤炭生产过程中抽掘采在时间和空间上的匹配问题,推进煤矿区煤层气应抽尽抽;突破低浓度瓦斯提纯和利用关键工艺技术,提高CH4利用率和利用量,推进煤矿区煤层气(煤矿瓦斯)应用尽用,实现CH4零排放;攻克废弃(关闭)矿井煤层气资源评价和抽采技术,推进关闭矿井CH4高效抽采利用,减少CH4通过煤矿巷道和地层裂缝向大气中逸散;加强大气级、场地级和设备级CH4排放监测、统计、校验、模拟等基础技术研究,为煤矿CH4排放监督和管理提供基础手段。

    当前的煤炭开采工艺和方法,从原理上不可避免地消耗能源和引起CH4排空,必须加大探索新的开采方法和技术,从原理上减少煤炭开采过程中的能源消耗和瓦斯排空,支撑煤炭开采节能降耗和低碳化。

    加快探索煤炭深部原位流态化开采理论和技术,攻克煤炭资源流态化迴行开采工艺、煤炭资源原位物理流态化工艺和技术、煤炭资源原位气化工艺和技术、煤炭资源原位液化工艺和技术等[42-43],推进煤炭资源以液体、气体及电能的方式从地下输出,实现煤炭资源的清洁低碳高效开发利用;突破煤与瓦斯物理流态化同采方法和技术,通过井巷工程共用、复用,降低单一煤炭或煤层气开采的能源消耗,推进低碳、低生态损害的煤与瓦斯协同开采。

    煤矿区除煤炭资源外,还有大量的土地、风、光等资源,具有发展可再生能源的先天优势[44-45],推进煤矿区煤与新能源耦合利用,是降低单位产品碳排放强度的重要途径。

    研发和示范煤矿地下水库电力调峰技术、煤矿区煤与太阳能光热耦合发电技术、煤矿区煤与风能耦合发电技术、煤矿区煤与地热能耦合发电/供热技术、煤矿区风能/太阳能制氢与煤清洁转化耦合技术等,支撑新能源发展,同时,减少单一以煤为原料生产产品的能源单耗和碳排放。进一步推进煤炭开发过程使用煤与新能源耦合发电的电力、氢能,减少燃煤、燃油,降低碳排放。

    煤炭利用过程产生的CO2量大,难以实现全部利用,基于煤矿区的地下空间优势,应提早布局探索煤矿区碳处置与封存技术,为无法利用的CO2提供最终的处置方案。

    探索煤矿深部原位CO2与CH4制氢新原理和技术,将CO2与CH4转化为无碳的氢能;攻克CO2矿化发电新理论与技术,在煤矿区实现CO2能源化再利用和固碳;突破高效CO2电化学捕集新原理新技术,实现煤矿区煤炭利用的低成本高效碳捕集;研究采空区CO2封存原理与控制技术、煤炭开采与采空区CO2充填协同方法,推进在适宜的煤矿区进行大规模CO2封存。

    基于煤炭开发全生命周期碳排放清单分析方法,从生产用能、瓦斯排放及矿后活动3个方面,估算了煤炭开发过程碳排放特征,分析了不同环节单位煤炭产品碳排放量及其变化趋势,并提出了煤炭开发过程碳中和发展的技术路径。主要结论如下:

    1)生产用能碳排放由2010年的2.64×108 t,先增加到2011年的2.79×108 t,随后逐渐降低到2016年的2.25×108 t,而后逐渐增加到2020年的2.57×108 t;对应生产用能碳排放强度由2010年的81.5 kg/t,先降低到2016年的66.5 kg/t,而后呈现增加趋势,逐渐增加到2020年的65.4 kg/t。

    2)瓦斯排放(碳排放)量由2010年的4.01×108 t,先波动降低到2016年的2.81×108 t,而后随着瓦斯抽采利用率的提高,逐渐降低到2020年的2.65×108 t;对应的瓦斯排放(碳排放)强度呈现逐渐降低的趋势,由2010年的123.7 kg/t,逐渐降低到2020年的67.6 kg/t。

    3)矿后活动碳排放量由2010年的7.0×107 t,先波动增加到2014年的7.7×107 t,随后降低到2016年的6.4×107 t,此后增加到2020年的7.1×107 t;对应的矿后活动碳排放强度呈现逐渐降低的趋势,由2010年的21.5 kg/t,逐渐降低到2020年的18.0 kg/t。

    4)煤炭开发过程碳排放总量由2010年的7.35×108 t,先增加到2011年的7.59×108 t,随后逐渐降低到2016年的5.70×108t,而后再缓慢增加到2020年的5.93×108 t;对应煤炭开发碳排放强度呈现逐年降低的趋势,由2010年的226.7 kg/t,快速降低到2016年的169.6 kg/t,此后缓慢降低到2020年的151.1 kg/t。

    5)基于碳排放量和排放特征,煤炭开发过程低碳发展可采取5大技术途径,即大力推广应用煤炭开发节能提效技术,减少煤炭开发能源消耗带来的碳排放;持续攻关煤矿瓦斯抽采利用技术,减少煤炭开发过程CH4排放;加快探索煤炭开发新技术,从原理上改变碳排放特性;有序研发和示范煤矿区煤与新能源耦合利用技术,降低单位产品碳排放系数;尽早谋划煤矿区CO2捕集、利用与固化、封存技术,形成煤炭行业独有的碳中和实现途径。

  • 图  1   煤层气(煤矿瓦斯)抽采量及利用率

    Fig.  1   Coalbed methane (coal mine gas) extraction volume and utilization rate

    下载: 全尺寸图片

    图  2   煤炭开发过程碳排放计算模型范围、边界及输入输出[24]

    Fig.  2   Scope, boundary and input and output of the carbon emission calculation model during the coal development process[24]

    下载: 全尺寸图片

    图  3   生产用能碳排放量及排放强度

    Fig.  3   Carbon emissions and emission intensity of energy used in production

    下载: 全尺寸图片

    图  4   不同机构煤矿瓦斯溢出量估算量对比

    Fig.  4   Comparison of the estimated amount of coal mine gas emissions from different institutions

    下载: 全尺寸图片

    图  5   煤炭开发碳排放量及排放强度

    Fig.  5   Carbon emissions and emission intensity of coal development

    下载: 全尺寸图片

    图  6   煤炭开发过程碳排放结构及变化趋势

    Fig.  6   Carbon emission structure and change trend during coal development

    下载: 全尺寸图片
  • [1] 中共中央 国务院. 关于完整准确全面贯彻新发展理念做好碳达峰碳中和工作的意见[EB/OL].(2021–09–22)[2021–10–22].http://www.gov.cn/zhengce/2021-10/24/content_5644613.htm.

    The Central Committee of The State Council of China.Working guidance for carbon dioxide peaking and carbon neutrality in full and faithful implementation of the new development philosophy[EB/OL].(2021–09–22)[2021–10–22].http://www.gov.cn/zhengce/2021-10/24/content_5644613.htm.
    [2] Liu Junling,Wang Ke,Zou Ji,et al.The implications of coal consumption in the power sector for China’s CO2 peaking target[J].Applied Energy,2019,253:113518. doi: 10.1016/j.apenergy.2019.113518
    [3] 于胜民,朱松丽,张俊龙.中国井工煤矿开采过程的二氧化碳逃逸排放因子研究[J].中国能源,2018,40(5):10–16. doi: 10.3969/j.issn.1003-2355.2018.05.002

    Yu Shengmin,Zhu Songli,Zhang Junlong.Research on the emission factors of carbon dioxide escape during the mining process of China’s coal mines[J].Energy of China,2018,40(5):10–16 doi: 10.3969/j.issn.1003-2355.2018.05.002
    [4] Miller S M,Michalak A M,Detmers R G,et al.China’s coal mine methane regulations have not curbed growing emissions[J].Nature Communications,2019,10(1):1–8. doi: 10.1038/s41467-018-07882-8
    [5] Gao Junlian,Guan Chengshe,Zhang Bo.China’s CH4 emissions from coal mining:A review of current bottom-up inventories[J].Science of The Total Environment,2020,725:138295. doi: 10.1016/j.scitotenv.2020.138295
    [6] Wang Ke,Zhang Jianjun,Cai Bofeng,et al.Emission factors of fugitive methane from underground coal mines in China:Estimation and uncertainty[J].Applied Energy,2019,250:273–282. doi: 10.1016/j.apenergy.2019.05.024
    [7] Lin Jiang,Khanna N,Liu Xu,et al.China’s non-CO2 greenhouse gas emissions:Future trajectories and mitigation options and potential[J].Scientific Reports,2019,9(1):1–10.
    [8] 高俊莲,徐向阳,郑凤琴,等.基于全生命周期的煤炭碳排放清单计算与不确定性分析[J].中国煤炭,2017,43(6):22–26. doi: 10.3969/j.issn.1006-530X.2017.06.005

    Gao Junlian,Xu Xiangyang,Zheng Fengqin,et al.Coal carbon emission inventory calculation and uncertainties analysis based on lifecycle analysis[J].China Coal,2017,43(6):22–26 doi: 10.3969/j.issn.1006-530X.2017.06.005
    [9] 柳君波,徐向阳,霍志佳等.中国煤炭格局变化对煤矿甲烷排放的影响及原因[J].生态经济,2021,37(7):176–182. https://www.cnki.com.cn/Article/CJFDTOTAL-STJJ202107027.htm

    Liu Junbo,Xu Xiangyang,Huo Zhijia,et al.Influence of coal pattern change on coal mine methane emission in China[J].Ecological Economy,2021,37(7):176–182 https://www.cnki.com.cn/Article/CJFDTOTAL-STJJ202107027.htm
    [10] 王宁.基于区域因子分析的我国煤矿甲烷排放与控排对策研究[D].北京:中国矿业大学(北京),2013.

    Wang Ning.Research on China coal mine methane emissions and control measures based on regional factor analysis[D].Beijing:China University of Mining and Technology (Beijing),2013.
    [11] 李学武,张伟,李跃,等.基于涌出因子分析煤炭矿区甲烷利用及减排[J].煤炭学报,2014,39(增刊2):390–396. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB2014S2016.htm

    Li Xuewu,Zhang Wei,Li Yue,et al.CMM utilization and emission reduction based on gush factor analysis[J].Journal of China Coal Society,2014,39(Supp2):390–396 https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB2014S2016.htm
    [12] 马翠梅,戴尔阜,刘乙辰,等.中国煤炭开采和矿后活动甲烷逃逸排放研究[J].资源科学,2020,42(2):311–322. doi: 10.18402/resci.2020.02.10

    Ma Cuimei,Dai Erfu,Liu Yichen,et al.Methane fugitive emissions from coal mining and post-mining activities in China[J].Resources Science,2020,42(2):311–322 doi: 10.18402/resci.2020.02.10
    [13] Zhu Tao,Bian Wenjing,Zhang Shuqing,et al.An improved approach to estimate methane emissions from coal mining in China[J].Environmental Science & Technology,2017,51(21):12072–12080.
    [14] 张媛.煤炭开采对矿区生态系统的碳排放效应研究[D].晋中:山西农业大学,2019.

    Zhang Yuan.Study on carbon emission effect of coal mining in mining area ecosystem[D].Jinzhong:Shanxi Agricultural University,2019.
    [15] 司荣军,李润之.低浓度含氧瓦斯爆炸动力特性及防控关键技术[J].煤炭科学技术,2020,48(10):17–36. https://www.cnki.com.cn/Article/CJFDTOTAL-MTKJ202010002.htm

    Si Rongjun,Li Runzhi.Dynamic characteristics of low-concentration oxygen-containing gas explosion and key prevention and control key technologies[J].Coal Science and Technology,2020,48(10):17–36 https://www.cnki.com.cn/Article/CJFDTOTAL-MTKJ202010002.htm
    [16] 陈浮,于昊辰,卞正富,等.碳中和愿景下煤炭行业发展的危机与应对[J].煤炭学报,2021,46(6):1808–1820. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202106013.htm

    Chen Fu,Yu Haochen,Bian Zhengfu,et al.How to handles the crisis of coal industry in China under the vision of carbon neutrality[J].Journal of China Coal Society,2021,46(6):1808–1820 https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202106013.htm
    [17] 谢和平,任世华,谢亚辰,等.碳中和目标下煤炭行业发展机遇[J].煤炭学报,2021,46(6):1808–1820. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202107015.htm

    Xie Heping,Ren Shihua,Xie Yachen,et al.Development opportunities of the coal industry towards the goal of carbon neutrality[J].Journal of China Coal Society,2021,46(6):1808–1820 https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202107015.htm
    [18] 魏英敏.消费伦理与保护生态环境[J].中国矿业大学学报(社会科学版),2000,2(3):3–7. doi: 10.3969/j.issn.1009-105X.2000.03.001

    Wei Yingmin.The ethics of consumption and the protection of ecological environment[J].Journal of China University of Mining & Technology(Social Sciences),2000,2(3):3–7 doi: 10.3969/j.issn.1009-105X.2000.03.001
    [19] Mi Zhifu,Zheng Jiali,Meng Jing,et al.China’s energy consumption in the new normal[J].Earth’s Future,2018,6(7):1007–1016. doi: 10.1029/2018EF000840
    [20] Wang Bing,Cui Chaoqun,Zhao Yixin,et al.Carbon emissions accounting for China’s coal mining sector:Invisible sources of climate change[J].Natural Hazards,2019,99(3):1345–1364. doi: 10.1007/s11069-018-3526-2
    [21] Liu Gang,Peng Shushi,Lin Xin,et al.Recent slowdown of anthropogenic methane emissions in China driven by stabilized coal production[J].Environmental Science & Technology Letters,2021,8(9):739–746.
    [22] Sheng Jianxiong,Song Shaojie,Zhang Yuzhong,et al.Bottom-up estimates of coal mine methane emissions in China:A gridded inventory,emission factors,and trends[J].Environmental Science & Technology Letters,2019,6(8):473–478.
    [23] 任世华.煤炭资源开发利用效率分析评价模型研究[J].中国能源,2015,37(2):33–36. doi: 10.3969/j.issn.1003-2355.2015.02.006

    Ren Shihua.Research on analysis and evaluation model of coal resources development and utilization efficiency[J].Energy of China,2015,37(2):33–36 doi: 10.3969/j.issn.1003-2355.2015.02.006
    [24] 煤炭工业洁净煤工程技术研究中心.煤炭开发利用碳排放清单分析及减排潜力研究[R].北京:煤炭科学研究总院,2012.
    [25] 任世华,罗腾,赵路正.煤炭开发利用碳减排潜力分析[J].中国能源,2013,35(11):24–27. doi: 10.3969/j.issn.1003-2355.2013.11.006

    Ren Shihua,Luo Teng,Zhao Luzheng.Analysis on the potential of carbon emission reduction in coal development and utilization[J].Energy of China,2013,35(11):24–27 doi: 10.3969/j.issn.1003-2355.2013.11.006
    [26] 樊金璐,孙键,赵娆.我国煤炭行业全生命周期碳排放与碳流通图[J].煤炭经济研究,2017,37(9):34–37. https://www.cnki.com.cn/Article/CJFDTOTAL-MTJN201709007.htm

    Fan Jinlu,Sun Jian,Zhao Rao.Carbon emission and carbon flow chart for whole life cycle of China coal industry[J].Coal Economic Research,2017,37(9):34–37 https://www.cnki.com.cn/Article/CJFDTOTAL-MTJN201709007.htm
    [27] 国家发展和改革委员会.中国煤炭生产企业温室气体排放核算方法与报告指南(试行)[R].北京:国家发展与改革委员会,2014.
    [28] 国家统计局.中国统计年鉴(2021)[M].北京:中国统计出版社,2021.

    National Bureau of Statistics.China statistical yearbook (2021)[M].Beijing:China Statistics Press,2021
    [29] IPCC.2006 IPCC guidelines for national greenhouse gas inventories[M].Kanagawa:The Institute for Global Environmental Strategies,2006.
    [30] IPCC.2013 supplement to the 2006 IPCC guidelines for national greenhouse gas inventories[M].Geneva:Intergovernmental Panel on Climate Change,2013.2014
    [31] 杨永均,张绍良,侯湖平.煤炭开采的温室气体逸散排放估算研究[J].中国煤炭,2014,40(1):114–117. doi: 10.3969/j.issn.1006-530X.2014.01.032

    Yang Yongjun,Zhang Shaoliang,Hou Huping.The research on the estimation of fugitive greenhouse gas emissions from coal mining[J].China Coal,2014,40(1):114–117 doi: 10.3969/j.issn.1006-530X.2014.01.032
    [32] 才庆祥,刘福明,陈树召.露天煤矿温室气体排放计算方法[J].煤炭学报,2012,37(1):103–106. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201201025.htm

    Cai Qingxiang,Liu Fuming,Chen Shuzhao.Calculation method of greenhouse gas emission in open pit coal mines[J].Journal of China Coal Society,2012,37(1):103–106 https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201201025.htm
    [33] 张振芳.露天煤矿碳排放量核算及碳减排途径研究[D].徐州:中国矿业大学,2013.

    Zhang Zhenfang.Study on carbon emissions accounting and carbon emission reduction approach of surface coal mine[D].Xuzhou:China University of Mining and Technology,2013.
    [34] 张凯,李全生.碳约束条件下我国煤炭清洁高效开发利用战略研究[J].煤炭经济研究,2019,39(11):10–14. https://www.cnki.com.cn/Article/CJFDTOTAL-MTJN201911003.htm

    Zhang Kai,Li Quansheng.Strategic research on clean and efficient development and utilization of coal in China under carbon constraint conditions[J].Coal Economic Research,2019,39(11):10–14 https://www.cnki.com.cn/Article/CJFDTOTAL-MTJN201911003.htm
    [35] 葛世荣,刘洪涛,刘金龙,等.我国煤矿生产能耗现状分析及节能思路[J].中国矿业大学学报,2018,47(1):9–14. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD201801002.htm

    Ge Shirong,Liu Hongtao,Liu Jinlong,et al.Energy consumption and energy-saving strategies for coal mine production in China[J].Journal of China University of Mining & Technology,2018,47(1):9–14 https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD201801002.htm
    [36] 夏伦娣.“双代煤”工程温室气体和大气污染物的协同减排效益评估[D].邯郸:河北工程大学,2020.

    Xia Lundi.Evaluation of the benefit of the “coal to gas and electricity” project on the coordinated emission reduction of greenhouse gases and air pollutants[D].Handan:Hebei University of Engineering,2020.
    [37] 国家煤矿安全监察局.全国煤矿矿井瓦斯等级鉴定资料汇编(2011)[R].北京:国家煤矿安全检查局,2012.
    [38] 吴金刚,毛俊睿.中国废弃煤矿瓦斯资源评价与抽采利用研究进展[J].煤矿安全,2021,52(7):162–169. https://www.cnki.com.cn/Article/CJFDTOTAL-MKAQ202107031.htm

    Wu Jingang,Mao Junrui.Research progress on gas resource evaluation,extraction and utilization in abandoned coal mines in China[J].Safety in Coal Mines,2021,52(7):162–169 https://www.cnki.com.cn/Article/CJFDTOTAL-MKAQ202107031.htm
    [39] Peng Shushi,Piao Shilong,Bousquet P,et al.Inventory of anthropogenic methane emissions in mainland China from 1980 to 2010[J].Atmospheric Chemistry and Physics,2016,16(22):14545–14562. doi: 10.5194/acp-16-14545-2016
    [40] 国家发展和改革委员会应对气候变化司.2005年中国温室气体清单研究[M].北京:中国环境出版社,2014.
    [41] Zhou Aitao,Hu Jiaying,Wang Kai.Carbon emission assessment and control measures for coal mining in China[J].Environmental Earth Sciences,2020,79(19):1–15.
    [42] 王金华,谢和平,刘见中,等.煤炭近零生态环境影响开发利用理论和技术构想[J].煤炭学报,2018,43(5):1198–1209. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201805002.htm

    Wang Jinhua,Xie Heping,Liu Jianzhong,et al.Coal development and utilization theory and technical system of near-zero ecological environment impact[J].Journal of China Coal Society,2018,43(5):1198–1209 https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201805002.htm
    [43] 谢和平,鞠杨,高明忠,等.煤炭深部原位流态化开采的理论与技术体系[J].煤炭学报,2018,43(5):1210–1219. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201805003.htm

    Xie Heiping,Ju Yang,Gao Mingzhong,et al.Theories and technologies for in-situ fluidized mining of deep underground coal resources[J].Journal of China Coal Society,2018,43(5):1210–1219 https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201805003.htm
    [44] 谢和平,高明忠,刘见中,等.煤矿地下空间容量估算及开发利用研究[J].煤炭学报,2018,43(6):1487–1503. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201806001.htm

    Xie Heping,Gao Mingzhong,Liu Jianzhong,et al.Research on exploitation and volume estimation of underground space in coal mines[J].Journal of China Coal Society,2018,43(6):1487–1503 https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201806001.htm
    [45] 谢和平,侯正猛,高峰,等.煤矿井下抽水蓄能发电新技术:原理、现状及展望[J].煤炭学报,2015,40(5):965–972. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201505001.htm

    Xie Heping,Hou Zhengmeng,Gao Feng,et al.A new technology of pumped-storage power in underground coal mine:Principles,present situation and future[J].Journal of China Coal Society,2015,40(5):965–972 https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201505001.htm
图(6)

本文结构

    /

    返回文章
    返回