一种农区种植业碳氮足迹评价方法及系统

本发明涉及环境核算，具体涉及一种农区种植业碳氮足迹评价方法及系统。

背景技术：

1、我国农业生产过程产生的温室气体占全国温室气体排放总量的近1/5，其中以二氧化碳、甲烷、氧化亚氮为主。农业源排放的甲烷占全国甲烷排放总量的近80％，农业源排放的氧化亚氮占全国氧化亚氮排放的近90％。

2、我国农业因高强度翻耕、作物连作、有机投入不足、秸秆焚烧、大量农用品投入导致氮磷面源污染和碳排放加剧，为积极稳步推进农业可持续发展，明晰农业碳氮足迹时空特征，有针对性行的提出减污降碳措施是实现农业可持续发展亟待解决的关键问题。

3、目前，国内对于种植业碳氮足迹仅限于收支核算阶段，并且没有形成一套系统、全面、科学的核算方法，导致核算结果无法体现区域农业生态环境现状。

技术实现思路

1、针对现有技术中的上述不足，本发明的目的在于提供一种农区种植业碳氮足迹评价方法及系统，旨在建立标准化农区种植业碳氮数据库，对农区种植业碳氮的输入-循环-输出进行精准核算，利用核算结果定量化评价农区种植生产过程中的碳氮排放及原位消纳能力。

2、为了达到上述发明目的，本发明采用的技术方案为：

3、第一方面，提供一种农区种植业碳氮足迹评价方法，其包括以下步骤：

4、s1，农区种植业生产过程碳氮源汇识别与数据获取；

5、s2，数据结构化入库与优化存储；

6、s3，数据调取与碳氮足迹核算；

7、s4，农区种植业碳氮足迹评价。

8、进一步地，步骤s1具体包括以下子步骤：

9、s11，确定农区种植业生产过程中是否有农用品投入，包括化肥、农药、农膜、有机肥、种子；获取农用品投入量，包括氮肥施用量、复合肥施用量、尿素施用量、磷肥施用量、钾肥施用量、农药使用量、农膜使用量、有机肥使用量、种子用量；

10、s12，确定农区种植业生产过程中是否有能源投入，包括农机柴油、灌溉电力；获取能源投入量，包括农机柴油消耗量、灌溉电力消耗量；

11、s13，确定农区是否种植豆科作物；获取豆科作物种植面积；

12、s14，确定农区是否种植水稻作物；获取水稻种植面积；

13、s15，确定农区是否存在盐碱地；获取盐碱地面积；

14、s16，确定农区种植业生产过程中是否有役畜投入；获取役畜投入数量；

15、s17，确定农区种植业秸秆处置方式；获取作物产量、秸秆还田比例、秸秆露天焚烧比例、秸秆能源化利用比例；

16、s18，获取农区种植业基本信息数据，包括种植面积、土壤有机质含量。

17、进一步地，步骤s2具体包括以下子步骤：

18、s21，建立农区种植业碳氮足迹结构化数据库，包括基础数据库、参数数据库、足迹评价数据库；

19、s22，将s1步骤中获取的数据导入至基础数据库，导入的数据信息包括数据名称、区域位置、年份、数据值、数据量纲、备注信息；

20、s23，建立参数数据库，包括农区种植业碳氮足迹核算过程中所需的排放因子、含量比例、等温潜势转换系数数据；

21、s24，建立足迹评价数据库，包括核算指标数据库、衍生指标数据库、评价指标数据库；

22、s25，对s21步骤中结构化数据库中的数据进行存储编码，优化存储结构。

23、进一步地，步骤s3具体包括以下子步骤：

24、s31，用投入产出和物质循环的概念将碳氮足迹核算划分为能源消耗碳排放量、役畜投入碳氮排放量、物质投入碳氮源汇量、元素循环碳氮排放量、作物产出碳氮消纳量、秸秆处置碳氮源汇量，计算结果存储于s24所述的足迹评价数据库中的核算指标数据库；

25、s32，根据表1调取数据，计算能源消耗碳排放量，包括农机柴油消耗碳排放量和灌溉电力碳排放量，计算方法为：

26、能源消耗碳排放量＝农机柴油消耗碳排放量+灌溉电力碳排放量

27、农机柴油消耗碳排放量＝农机柴油消耗量×柴油碳排放因子

28、灌溉电力碳排放量＝灌溉电力消耗量×电力碳排放因子

29、表1能源消耗碳排放量计算数据调取

30、 调取指标 对应存储编码 调取数据库位置 数据库编码 农机柴油消耗量 bd_diesel 基础数据库 basic data 灌溉电力消耗量 bd_ele 基础数据库 basic data 柴油碳排放因子 pd_diesel 参数数据库 parameter data 电力碳排放因子 pd_ele 参数数据库 parameter data

31、s33，根据表2调取数据，计算役畜投入碳氮排放量，包括役畜投入碳排放量和役畜投入氮排放量；其中，役畜投入碳排放量，包括肠道发酵甲烷排放量和粪便管理甲烷排放量；役畜投入氮排放量为粪便管理氧化亚氮排放量；计算方法为：

32、役畜投入碳排放量＝肠道发酵甲烷排放量+粪便管理甲烷排放量

33、肠道发酵甲烷排放量＝役畜投入数量×肠道发酵甲烷排放因子

34、粪便管理甲烷排放量＝役畜投入数量×粪便管理甲烷排放因子

35、役畜投入氮排放量＝粪便管理氧化亚氮排放量＝役畜投入数量×粪便管理氧化亚氮排放因子

36、表2役畜投入碳氮排放量计算数据调取

37、s34，根据表3调取数据，计算物质投入碳氮源汇量，包括物质投入碳源汇量和物质投入氮汇量；其中物质投入碳源汇量，包括有机肥碳汇量、种子碳汇量、盐碱地碳汇量、氮肥生产碳排放量、磷肥生产碳排放量、钾肥生产碳排放量、复合肥生产碳排放量、农药生产碳排放量、农膜生产碳排放量；物质投入氮汇量，包括有机肥氮汇量、氮肥氮汇量、复合肥氮汇量、种子氮汇量、大气沉降氮汇量、生物固氮量；计算方法为：

38、物质投入碳源汇量＝有机肥碳汇量+种子碳汇量+盐碱地碳汇量-农用品生产碳排放量

39、农用品生产碳排放量＝氮肥生产碳排放量+磷肥生产碳排放量+钾肥生产碳排放量+复合肥生产碳排放量+农药生产碳排放量+农膜生产碳排放量

40、物质投入氮汇量＝有机肥氮汇量+氮肥氮汇量+复合肥氮汇量+种子氮汇量+大气沉降氮汇量+生物固氮量

41、氮肥生产碳排放量＝氮肥施用量×氮肥生产碳排放因子

42、氮肥氮汇量＝氮肥施用量

43、磷肥生产碳排放量＝磷肥施用量×磷肥生产碳排放因子

44、钾肥生产碳排放量＝钾肥施用量×钾肥生产碳排放因子

45、复合肥生产碳排放量＝复合肥施用量×复合肥生产碳排放因子

46、复合肥氮汇量＝复合肥施用量×复合肥含氮量

47、农药生产碳排放量＝农药使用量×农药生产碳排放因子

48、农膜生产碳排放量＝农膜使用量×农膜生产碳排放因子

49、有机肥碳汇量＝有机肥使用量×有机肥含碳量

50、有机肥氮汇量＝有机肥使用量×有机肥含氮量

51、种子碳汇量＝种子用量×种子含碳量

52、种子氮汇量＝种子用量×种子含氮量

53、盐碱地碳汇量＝盐碱地面积×盐碱地碳汇因子

54、大气沉降氮汇量＝种植面积×氮沉降通量

55、生物固氮量＝豆科作物种植面积×生物固氮系数

56、表3物质投入碳氮源汇量计算数据调取

57、 调取指标 对应存储编码 调取数据库位置 数据库编码 氮肥施用量 bd_fern 基础数据库 basic data 磷肥施用量 bd_ferp 基础数据库 basic data 钾肥施用量 bd_ferk 基础数据库 basic data 复合肥施用量 bd_ferm 基础数据库 basic data 农药使用量 bd_pest 基础数据库 basic data 农膜使用量 bd_film 基础数据库 basic data 有机肥使用量 bd_om 基础数据库 basic data 种子用量 bd_seed 基础数据库 basic data 盐碱地面积 bd_salt 基础数据库 basic data 种植面积 bd_area 基础数据库 basic data 豆科作物种植面积 bd_bean 基础数据库 basic data 氮肥生产碳排放因子 pd_fernc 参数数据库 parameter data 磷肥生产碳排放因子 pd_ferpc 参数数据库 parameter data 钾肥生产碳排放因子 pd_ferkc 参数数据库 parameter data 复合肥生产碳排放因子 pd_fermc 参数数据库 parameter data 复合肥含氮量 pd_fermn 参数数据库 parameter data 农药生产碳排放因子 pd_pest 参数数据库 parameter data 农膜生产碳排放因子 pd_film 参数数据库 parameter data 有机肥含碳量 pd_omc 参数数据库 parameter data 有机肥含氮量 pd_omn 参数数据库 parameter data 种子含碳量 pd_seedc 参数数据库 parameter data 种子含氮量 pd_seedn 参数数据库 parameter data 盐碱地碳汇因子 pd_salt 参数数据库 parameter data 氮沉降通量 pd_dep 参数数据库 parameter data 生物固氮系数 pd_bio 参数数据库 parameter data

58、s35，根据表4调取数据，计算元素循环碳氮排放量，包括元素循环碳排放量和元素循环氮排放量；其中元素循环碳排放量，包括有机质矿化碳排放量、尿素施用碳排放量、稻田甲烷排放量；元素循环氮排放量，包括有机质矿化氮排放量、施肥导致氧化亚氮排放量、施肥导致氨气排放量、施肥导致硝酸盐排放量；计算方法为：

59、元素循环碳排放量＝有机质矿化碳排放量+尿素施用碳排放量+稻田甲烷排放量

60、元素循环氮排放量＝有机质矿化氮排放量+施肥导致氧化亚氮排放量+施肥导致氨气排放量+施肥导致硝酸盐排放量

61、有机质矿化碳排放量＝土壤有机质含量×土壤矿化率

62、有机质矿化氮排放量＝有机质矿化氮排放量＝[0.914×log(土壤有机质含量×0.58)+1.058]×种植面积

63、尿素施用碳排放量＝尿素施用量×尿素碳排放因子

64、稻田甲烷排放量＝水稻种植面积×稻田甲烷排放因子

65、施肥导致氧化亚氮排放量＝氮肥施用量×施肥氧化亚氮排放因子

66、施肥导致氨气排放量＝氮肥施用量×施肥氨气排放因子

67、施肥导致硝酸盐排放量＝氮肥施用量×施肥硝酸盐排放因子

68、表4元素循环碳氮排放量计算数据调取

69、 调取指标 对应存储编码 调取数据库位置 数据库编码 土壤有机质含量 bd_som 基础数据库 basic data 尿素施用量 bd_urea 基础数据库 basic data 水稻种植面积 bd_rice 基础数据库 basic data 氮肥施用量 bd_fern 基础数据库 basic data 土壤矿化率 pd_som 参数数据库 parameter data 尿素碳排放因子 pd_urea 参数数据库 parameter data 稻田甲烷排放因子 pd_rice 参数数据库 parameter data 施肥氧化亚氮排放因子 pd_n2o 参数数据库 parameter data 施肥氨气排放因子 pd_nh3 参数数据库 parameter data 施肥硝酸盐排放因子 pd_no3 参数数据库 parameter data

70、s36，根据表5调取数据，计算作物产出碳氮消纳量，包括作物产出碳消纳量和作物产出氮消纳量；其中，作物产出碳消纳量，包括籽粒碳消纳量、秸秆碳消纳量；作物产出氮消纳量，包括籽粒氮消纳量、秸秆氮消纳量；计算方法为：

71、作物产出碳消纳量＝籽粒碳消纳量+秸秆碳消纳量

72、作物产出氮消纳量＝籽粒氮消纳量+秸秆氮消纳量

73、籽粒碳消纳量＝作物产量×种子含碳量

74、籽粒氮消纳量＝作物产量×种子含氮量

75、秸秆碳消纳量＝作物产量×草谷比×秸秆含碳量

76、秸秆氮消纳量＝作物产量×草谷比×秸秆含氮量

77、表5作物产出碳氮消纳量计算数据调取

78、 调取指标 对应存储编码 调取数据库位置 数据库编码 作物产量 bd_yield 基础数据库 basic data 种子含碳量 pd_seedc 参数数据库 parameter data 种子含氮量 pd_seedn 参数数据库 parameter data 秸秆含碳量 pd_strawc 参数数据库 parameter data 秸秆含氮量 pd_strawn 参数数据库 parameter data 草谷比 pd_ratio 参数数据库 parameter data

79、s37，根据表6调取数据，计算秸秆处置碳氮源汇量，包括秸秆处置碳源汇量和秸秆处置氮源汇量；其中，秸秆处置碳源汇量，包括秸秆还田碳汇量、秸秆露天焚烧甲烷排放量、秸秆能源化碳减排量；秸秆处置氮源汇量，包括秸秆还田氮汇量、秸秆露天焚烧氧化亚氮排放量；计算方法为：

80、秸秆处置碳源汇量＝秸秆露天焚烧甲烷排放量×甲烷等温潜势-秸秆还田碳汇量-秸秆能源化碳减排量

81、秸秆处置氮源汇量＝秸秆露天焚烧氧化亚氮排放量-秸秆还田氮汇量

82、秸秆露天焚烧甲烷排放量＝作物产量×草谷比×秸秆露天焚烧比例×秸秆露天焚烧甲烷排放修正因子

83、秸秆露天焚烧氧化亚氮排放量＝作物产量×草谷比×秸秆露天焚烧比例×秸秆露天焚烧氧化亚氮排放因子

84、秸秆还田碳汇量＝作物产量×草谷比×秸秆还田比例×秸秆含碳量

85、秸秆还田氮汇量＝作物产量×草谷比×秸秆还田比例×秸秆含氮量

86、秸秆能源化碳减排量＝作物产量×草谷比×秸秆能源化比例×秸秆能源化碳减排系数

87、表6秸秆处置碳氮源汇量计算数据调取

88、

89、进一步地，步骤s4具体包括以下子步骤：

90、s41，根据s3步骤计算结果，计算农区种植业碳氮足迹衍生指标，包括单位面积碳消纳强度、单位面积氮消纳强度、碳排放量、氮排放量、温室气体净排放量，计算结果存储于s24所述的足迹评价数据库中的衍生指标数据库；计算方法为：

91、单位面积碳消纳强度＝作物产出碳消纳量÷种植面积

92、单位面积氮消纳强度＝作物产出氮消纳量÷种植面积

93、碳排放量＝能源消耗碳排放量+役畜投入碳排放量+元素循环碳排放量

94、氮排放量＝役畜投入氮排放量+元素循环氮排放量

95、温室气体净排放量＝能源消耗碳排放量+役畜投入碳排放量+役畜投入氮排放量×氮等温潜势+元素循环碳排放量+元素循环氮排放量×氮等温潜势+秸秆处置碳源汇量+农用品生产碳排放量；

96、s42，根据s3和s41步骤计算结果，计算农区种植业碳氮足迹评价指标，包括碳足迹、氮足迹、碳承载力、氮承载力、碳足迹指数、氮足迹指数，计算结果存储于s24所述的足迹评价数据库中的评价指标数据库；计算方法为：

97、碳足迹＝碳排放量÷单位面积碳消纳强度

98、氮足迹＝氮排放量÷单位面积氮消纳强度

99、碳承载力＝(作物产出碳消纳量-物质投入碳源汇量)÷单位面积碳消纳强度

100、氮承载力＝(作物产出氮消纳量-物质投入氮汇量)÷单位面积氮消纳强度

101、碳足迹指数＝碳足迹-碳承载力

102、氮足迹指数＝氮足迹-氮承载力；

103、s43，根据s42步骤计算结果，对农区种植业碳氮足迹进行评价，评价方法为：

104、

105、

106、第二方面，提供一种农区种植业碳氮足迹评价系统，其包括：

107、碳氮足迹源汇识别与数据获取模块，用于农区种植业生产过程碳氮源汇识别与数据获取；

108、碳氮足迹数据结构化入库与优化存储模块，用于数据结构化入库与优化存储；

109、碳氮足迹数据调取与碳氮足迹核算模块，用于数据调取与碳氮足迹核算；

110、农区种植业碳氮足迹评价模块，用于农区种植业碳氮足迹评价。

111、进一步地，碳氮足迹源汇识别与数据获取模块包括：

112、农用品投入确定与投入量获取子模块，用于确定农区种植业生产过程中是否有农用品投入，包括化肥、农药、农膜、有机肥、种子；获取农用品投入量，包括氮肥施用量、复合肥施用量、尿素施用量、磷肥施用量、钾肥施用量、农药使用量、农膜使用量、有机肥使用量、种子用量；

113、能源投入确定与投入量获取子模块，用于确定农区种植业生产过程中是否有能源投入，包括农机柴油、灌溉电力；获取能源投入量，包括农机柴油消耗量、灌溉电力消耗量；

114、豆科作物确定与获取子模块，用于确定农区是否种植豆科作物；获取豆科作物种植面积；

115、水稻作物确定与获取子模块，用于确定农区是否种植水稻作物；获取水稻种植面积；

116、盐碱地确定与获取子模块，用于确定农区是否存在盐碱地；获取盐碱地面积；

117、役畜投入确定与获取子模块，用于确定农区种植业生产过程中是否有役畜投入；获取役畜投入数量；

118、秸秆处置方式确定与获取子模块，用于确定农区种植业秸秆处置方式；获取作物产量、秸秆还田比例、秸秆露天焚烧比例、秸秆能源化利用比例；

119、种植业基本信息数据获取子模块，用于获取农区种植业基本信息数据，包括种植面积、土壤有机质含量。

120、进一步地，碳氮足迹数据结构化入库与优化存储模块包括：

121、种植业碳氮足迹结构化数据库建立子模块，用于建立农区种植业碳氮足迹结构化数据库，包括基础数据库、参数数据库、足迹评价数据库；

122、将碳氮足迹源汇识别与数据获取模块中获取的数据导入至基础数据库，导入的数据信息包括数据名称、区域位置、年份、数据值、数据量纲、备注信息；

123、参数数据库，包括农区种植业碳氮足迹核算过程中所需的排放因子、含量比例、等温潜势转换系数数据；

124、足迹评价数据库，包括核算指标数据库、衍生指标数据库、评价指标数据库；

125、数据存储编码与存储结构优化子模块，用于对结构化数据库中的数据进行存储编码，优化存储结构。

126、进一步地，碳氮足迹数据结构化入库与优化存储模块包括：

127、碳氮足迹核算划分子模块，用于将碳氮足迹核算划分为能源消耗碳排放量、役畜投入碳氮排放量、物质投入碳氮源汇量、元素循环碳氮排放量、作物产出碳氮消纳量、秸秆处置碳氮源汇量，计算结果存储于足迹评价数据库中的核算指标数据库；

128、能源消耗碳排放量计算子模块，用于根据表1调取数据，计算能源消耗碳排放量，包括农机柴油消耗碳排放量和灌溉电力碳排放量，计算方法为：

129、能源消耗碳排放量＝农机柴油消耗碳排放量+灌溉电力碳排放量

130、农机柴油消耗碳排放量＝农机柴油消耗量×柴油碳排放因子

131、灌溉电力碳排放量＝灌溉电力消耗量×电力碳排放因子

132、表1能源消耗碳排放量计算数据调取

133、 调取指标 对应存储编码 调取数据库位置 数据库编码 农机柴油消耗量 bd_diesel 基础数据库 basic data 灌溉电力消耗量 bd_ele 基础数据库 basic data 柴油碳排放因子 pd_diesel 参数数据库 parameter data 电力碳排放因子 pd_ele 参数数据库 parameter data

134、役畜投入碳氮排放量计算子模块，用于根据表2调取数据，计算役畜投入碳氮排放量，包括役畜投入碳排放量和役畜投入氮排放量；其中，役畜投入碳排放量，包括肠道发酵甲烷排放量和粪便管理甲烷排放量；役畜投入氮排放量为粪便管理氧化亚氮排放量；计算方法为：

135、役畜投入碳排放量＝肠道发酵甲烷排放量+粪便管理甲烷排放量

136、肠道发酵甲烷排放量＝役畜投入数量×肠道发酵甲烷排放因子

137、粪便管理甲烷排放量＝役畜投入数量×粪便管理甲烷排放因子

138、役畜投入氮排放量＝粪便管理氧化亚氮排放量＝役畜投入数量×粪便管理氧化亚氮排放因子

139、表2役畜投入碳氮排放量计算数据调取

140、

141、物质投入碳氮源汇量计算子模块，用于根据表3调取数据，计算物质投入碳氮源汇量，包括物质投入碳源汇量和物质投入氮汇量；其中物质投入碳源汇量，包括有机肥碳汇量、种子碳汇量、盐碱地碳汇量、氮肥生产碳排放量、磷肥生产碳排放量、钾肥生产碳排放量、复合肥生产碳排放量、农药生产碳排放量、农膜生产碳排放量；物质投入氮汇量，包括有机肥氮汇量、氮肥氮汇量、复合肥氮汇量、种子氮汇量、大气沉降氮汇量、生物固氮量；计算方法为：

142、物质投入碳源汇量＝有机肥碳汇量+种子碳汇量+盐碱地碳汇量-农用品生产碳排放量

143、农用品生产碳排放量＝氮肥生产碳排放量+磷肥生产碳排放量+钾肥生产碳排放量+复合肥生产碳排放量+农药生产碳排放量+农膜生产碳排放量

144、物质投入氮汇量＝有机肥氮汇量+氮肥氮汇量+复合肥氮汇量+种子氮汇量+大气沉降氮汇量+生物固氮量

145、氮肥生产碳排放量＝氮肥施用量×氮肥生产碳排放因子

146、氮肥氮汇量＝氮肥施用量

147、磷肥生产碳排放量＝磷肥施用量×磷肥生产碳排放因子

148、钾肥生产碳排放量＝钾肥施用量×钾肥生产碳排放因子

149、复合肥生产碳排放量＝复合肥施用量×复合肥生产碳排放因子

150、复合肥氮汇量＝复合肥施用量×复合肥含氮量

151、农药生产碳排放量＝农药使用量×农药生产碳排放因子

152、农膜生产碳排放量＝农膜使用量×农膜生产碳排放因子

153、有机肥碳汇量＝有机肥使用量×有机肥含碳量

154、有机肥氮汇量＝有机肥使用量×有机肥含氮量

155、种子碳汇量＝种子用量×种子含碳量

156、种子氮汇量＝种子用量×种子含氮量

157、盐碱地碳汇量＝盐碱地面积×盐碱地碳汇因子

158、大气沉降氮汇量＝种植面积×氮沉降通量

159、生物固氮量＝豆科作物种植面积×生物固氮系数

160、表3物质投入碳氮源汇量计算数据调取

161、 调取指标 对应存储编码 调取数据库位置 数据库编码 氮肥施用量 bd_fern 基础数据库 basic data 磷肥施用量 bd_ferp 基础数据库 basic data 钾肥施用量 bd_ferk 基础数据库 basic data 复合肥施用量 bd_ferm 基础数据库 basic data 农药使用量 bd_pest 基础数据库 basic data 农膜使用量 bd_film 基础数据库 basic data 有机肥使用量 bd_om 基础数据库 basic data 种子用量 bd_seed 基础数据库 basic data 盐碱地面积 bd_salt 基础数据库 basic data 种植面积 bd_area 基础数据库 basic data 豆科作物种植面积 bd_bean 基础数据库 basic data 氮肥生产碳排放因子 pd_fernc 参数数据库 parameter data 磷肥生产碳排放因子 pd_ferpc 参数数据库 parameter data 钾肥生产碳排放因子 pd_ferkc 参数数据库 parameter data 复合肥生产碳排放因子 pd_fermc 参数数据库 parameter data 复合肥含氮量 pd_fermn 参数数据库 parameter data 农药生产碳排放因子 pd_pest 参数数据库 parameter data 农膜生产碳排放因子 pd_film 参数数据库 parameter data 有机肥含碳量 pd_omc 参数数据库 parameter data 有机肥含氮量 pd_omn 参数数据库 parameter data 种子含碳量 pd_seedc 参数数据库 parameter data 种子含氮量 pd_seedn 参数数据库 parameter data 盐碱地碳汇因子 pd_salt 参数数据库 parameter data 氮沉降通量 pd_dep 参数数据库 parameter data 生物固氮系数 pd_bio 参数数据库 parameter data

162、元素循环碳氮排放量计算子模块，用于根据表4调取数据，计算元素循环碳氮排放量，包括元素循环碳排放量和元素循环氮排放量；其中元素循环碳排放量，包括有机质矿化碳排放量、尿素施用碳排放量、稻田甲烷排放量；元素循环氮排放量，包括有机质矿化氮排放量、施肥导致氧化亚氮排放量、施肥导致氨气排放量、施肥导致硝酸盐排放量；计算方法为：

163、元素循环碳排放量＝有机质矿化碳排放量+尿素施用碳排放量+稻田甲烷排放量

164、元素循环氮排放量＝有机质矿化氮排放量+施肥导致氧化亚氮排放量+施肥导致氨气排放量+施肥导致硝酸盐排放量

165、有机质矿化碳排放量＝土壤有机质含量×土壤矿化率

166、有机质矿化氮排放量＝有机质矿化氮排放量＝[0.914×log(土壤有机质含量×0.58)+1.058]×种植面积

167、尿素施用碳排放量＝尿素施用量×尿素碳排放因子

168、稻田甲烷排放量＝水稻种植面积×稻田甲烷排放因子

169、施肥导致氧化亚氮排放量＝氮肥施用量×施肥氧化亚氮排放因子

170、施肥导致氨气排放量＝氮肥施用量×施肥氨气排放因子

171、施肥导致硝酸盐排放量＝氮肥施用量×施肥硝酸盐排放因子

172、表4元素循环碳氮排放量计算数据调取

173、 调取指标 对应存储编码 调取数据库位置 数据库编码 土壤有机质含量 bd_som 基础数据库 basic data 尿素施用量 bd_urea 基础数据库 basic data 水稻种植面积 bd_rice 基础数据库 basic data 氮肥施用量 bd_fern 基础数据库 basic data 土壤矿化率 pd_som 参数数据库 parameter data 尿素碳排放因子 pd_urea 参数数据库 parameter data 稻田甲烷排放因子 pd_rice 参数数据库 parameter data 施肥氧化亚氮排放因子 pd_n2o 参数数据库 parameter data 施肥氨气排放因子 pd_nh3 参数数据库 parameter data 施肥硝酸盐排放因子 pd_no3 参数数据库 parameter data

174、作物产出碳氮消纳量计算子模块，用于根据表5调取数据，计算作物产出碳氮消纳量，包括作物产出碳消纳量和作物产出氮消纳量；其中，作物产出碳消纳量，包括籽粒碳消纳量、秸秆碳消纳量；作物产出氮消纳量，包括籽粒氮消纳量、秸秆氮消纳量；计算方法为：

175、作物产出碳消纳量＝籽粒碳消纳量+秸秆碳消纳量

176、作物产出氮消纳量＝籽粒氮消纳量+秸秆氮消纳量

177、籽粒碳消纳量＝作物产量×种子含碳量

178、籽粒氮消纳量＝作物产量×种子含氮量

179、秸秆碳消纳量＝作物产量×草谷比×秸秆含碳量

180、秸秆氮消纳量＝作物产量×草谷比×秸秆含氮量

181、表5作物产出碳氮消纳量计算数据调取

182、 调取指标 对应存储编码 调取数据库位置 数据库编码 作物产量 bd_yield 基础数据库 basic data 种子含碳量 pd_seedc 参数数据库 parameter data 种子含氮量 pd_seedn 参数数据库 parameter data 秸秆含碳量 pd_strawc 参数数据库 parameter data 秸秆含氮量 pd_strawn 参数数据库 parameter data 草谷比 pd_ratio 参数数据库 parameter data

183、秸秆处置碳氮源汇量计算子模块，用根据表6调取数据，计算秸秆处置碳氮源汇量，包括秸秆处置碳源汇量和秸秆处置氮源汇量；其中，秸秆处置碳源汇量，包括秸秆还田碳汇量、秸秆露天焚烧甲烷排放量、秸秆能源化碳减排量；秸秆处置氮源汇量，包括秸秆还田氮汇量、秸秆露天焚烧氧化亚氮排放量；所述计算方法为：

184、秸秆处置碳源汇量＝秸秆露天焚烧甲烷排放量×甲烷等温潜势-秸秆还田碳汇量-秸秆能源化碳减排量

185、秸秆处置氮源汇量＝秸秆露天焚烧氧化亚氮排放量-秸秆还田氮汇量

186、秸秆露天焚烧甲烷排放量＝作物产量×草谷比×秸秆露天焚烧比例×秸秆露天焚烧甲烷排放修正因子

187、秸秆露天焚烧氧化亚氮排放量＝作物产量×草谷比×秸秆露天焚烧比例×秸秆露天焚烧氧化亚氮排放因子

188、秸秆还田碳汇量＝作物产量×草谷比×秸秆还田比例×秸秆含碳量

189、秸秆还田氮汇量＝作物产量×草谷比×秸秆还田比例×秸秆含氮量

190、秸秆能源化碳减排量＝作物产量×草谷比×秸秆能源化比例×秸秆能源化碳减排系数

191、表6秸秆处置碳氮源汇量计算数据调取

192、

193、进一步地，农区种植业碳氮足迹评价模块包括：

194、农区种植业碳氮足迹衍生指标计算子模块，用于根据碳氮足迹数据调取与碳氮足迹核算模块计算结果，计算农区种植业碳氮足迹衍生指标，包括单位面积碳消纳强度、单位面积氮消纳强度、碳排放量、氮排放量、温室气体净排放量，计算结果存储于足迹评价数据库中的衍生指标数据库；计算方法为：

195、单位面积碳消纳强度＝作物产出碳消纳量÷种植面积

196、单位面积氮消纳强度＝作物产出氮消纳量÷种植面积

197、碳排放量＝能源消耗碳排放量+役畜投入碳排放量+元素循环碳排放量

198、氮排放量＝役畜投入氮排放量+元素循环氮排放量

199、温室气体净排放量＝能源消耗碳排放量+役畜投入碳排放量+役畜投入氮排放量×氮等温潜势+元素循环碳排放量+元素循环氮排放量×氮等温潜势+秸秆处置碳源汇量+农用品生产碳排放量

200、农区种植业碳氮足迹评价指标计算子模块，用于根据碳氮足迹数据调取与碳氮足迹核算模块和农区种植业碳氮足迹衍生指标计算子模块计算结果，计算农区种植业碳氮足迹评价指标，包括碳足迹、氮足迹、碳承载力、氮承载力、碳足迹指数、氮足迹指数，计算结果存储于足迹评价数据库中的评价指标数据库；计算方法为：

201、碳足迹＝碳排放量÷单位面积碳消纳强度

202、氮足迹＝氮排放量÷单位面积氮消纳强度

203、碳承载力＝(作物产出碳消纳量-物质投入碳源汇量)÷单位面积碳消纳强度

204、氮承载力＝(作物产出氮消纳量-物质投入氮汇量)÷单位面积氮消纳强度

205、碳足迹指数＝碳足迹-碳承载力

206、氮足迹指数＝氮足迹-氮承载力

207、农区种植业碳氮足迹评价子模块，用于根据农区种植业碳氮足迹评价指标计算子模块计算结果，对农区种植业碳氮足迹进行评价，评价方法为：

208、

209、

210、本发明的有益效果为：

211、1、本发明结合生命周期评价法、投入产出法、物质循环法等方法，对农区种植业碳氮源汇量进行计算，并对农区碳氮足迹与原位消纳能力进行评价。

212、2、本发明可实现农区种植业碳氮足迹评价方法中的单一指标或多个指标的核算，并实现核算指标的时空特征分析。

213、3、利用本发明的种植业碳氮足迹评价结果，可为种植业减污降碳治理提供数据支撑，并提高治理措施的针对性和准确性。