随着科技进步和人口增长，农业生态系统所产生的环境影响日益显著，已成为国内外关注热点。黑龙江省依靠其较高的农业机械化程度、独特的自然环境和优越的地理条件优势，已称为全国水稻种植面积、生产量和商品量大省之一，在全国稻米市场上具有举足轻重的地位[1]。然而，水稻作为全球重要的农作物之一，其生产全过程涉及大量能源消耗和污染物排放。因此，研究稻米生产过程中能源消耗与环境评价，对优化种植方式、降低能耗、改善环境、促进可持续农业具有重要意义。

生命周期评价（Life Cycle Assessment，简称LCA）是国际上通行的评估一个产品或是整体活动的、贯穿它的整个生命的环境后果的一种工具，为可持续生产提供途径和方法[2]。1997年Audsley等[3]将LCA应用涉及到重要的农产品的生命周期。之后，Williams等[4]提供了关于10种农业和园艺重要产品的详细报告，其中水稻是重点部分。大量研究围绕有关稻田温室气体排放，及相应的缓解建议与方案[5-8]，也有学者对美国、日本、意大利、瑞典、泰国等地的绿色稻米生产的LCA及常规管理优化等方面进行研究[9-15]。我国学者基于LCA对黑龙江、吉林、上海、江苏、湖南、台湾等地的稻米生产能耗及其温室气体排放进行了研究[16-24]，但目前国内稻米生产系统研究主要集中在碳、氮排放方面的影响，对于寒地稻米主要开展品种筛选、产量与环境互作方面的研究[25]，有关寒地稻米生产的环境影响方面的研究未见报道。

研究拟通过调查黑龙江寒地稻米生产全过程的环境排放，利用生命周期评价方法对稻米生产的6种环境影响进行评价，综合分析寒地稻米生产环境影响结构，为寒地稻米生产体系全程低碳管理和黑龙江生态农业发展提供科学依据。

1　材料和方法

1.1　研究地区概况

研究调查区域位于黑龙江省友谊县友谊农场（北纬 46°28′15″至 46°58′39″，东经 137°27′50″至132°15′38″），地处佳木斯市东南部，双鸭山市东北部，三江平原大片沼泽地边缘，土地面积18.88万hm2，是中国第一大农场。友谊农场位于黑龙江第二积温带，气候属亚寒带大陆性季风气候，春季多风少雨、干旱，夏季炎热湿润短暂，秋季凉爽多雨，昼夜温差较大，冬季严寒漫长，全年结冰期180 d左右，冻层厚度在1.5～2 m之间，无霜期120～135 d。常年平均气温3℃左右，最低气温平均0.8℃，最高气温平均4.3℃，极端最高气温37.3℃，极端最低气温-37.2℃。境内西南部积温高于东北，常年大于10℃有效积温2 170～2 700℃之间。友谊农场耕地面积为10.67万hm2，水稻种植面积6.67万hm2，占黑龙江水稻种植面积的1.66%，综合农业机械化率达98.2%，2014年友谊水稻产量57.33万吨，占黑龙江水稻产量的2.58%。

1.2　基本原理与研究方法

生命周期评价方法（LCA）已应用于选择性稻米耕种与食品工艺方法的环境效益方面的评价。根据ISO14040[26]，LCA主要包括四个阶段：目标与范围定义，生命周期清单，生命周期影响评价和结果解释。

今后也应加强对典型村域进行长期定位观测研究，实时掌握寒地不同区域不同农业生产模式、不同农业产业结构、不同能源消费模式对环境排放的影响，进一步拓展绿色生态农业系统研究的深度和广度。

GWP的主要来源是生产过程中田间释放、化肥生产以及碾米包装所排放的温室气体。前人研究表明，灌溉、土壤水分、施肥对温室气体排放影响较大[18]。化肥生产需要大量能耗，而我国的电能和热能生产主要来源煤炭的燃烧[30]。而且，过量的化肥使用所产生的N、P会侵入水体引起地表水富营养化和酸化[31-32]，而偏酸性、质地轻的水稻土反硝化也会引起温室气体排放[33]。因此，减少用水量与肥料施用量，是减少温室气体排放和富营养化的最有效措施。

 

表1　CML2001方法的分类、特征化与权重因子Table 1　Classification，characterization and weighed factor in the CML 2001 method

  

环境影响类型ADP GWP ODP AP EP POCP单位kg Sbeqkg CO2 eqkg R11eqkg SO2 eqkg PO4-3eqkg C2H2 eq权重1.5 10.0 3.0 2.0 7.0 3.0

如表1所示，研究考虑的6类环境影响类型为：非生物资源耗竭（Abiotic Depletion，ADP），是指所用的不可再生资源；全球变暖潜能（Global Warming Potential，GWP），用来指示温室效应影响；臭氧层消耗潜能（Ozone Depletion Potential，ODP），用来指示平流层臭氧损耗；酸化效应（Acidification Potential，AP），用来指示酸雨现象；富营养化（Eutrophication Potential，EP），用来指示地表水富营养化；光化学臭氧潜力（Photochemical Ozone Creation Potential，POCP），用来指示光烟雾形成。

因此，研究系统包括4个子系统，农资生产、水稻种植、稻米加工、运输阶段。系统边界和输入释放种类如图1所示。

1.3　系统边界与功能单位

由图5可知，田间释放环节的GWP、AP、EP和POCP的影响值为寒地稻米生产各环节的最高值，分别占68.50%、39.18%、72.27%和89.21%.化肥生产阶段具有最高的ADP和ODP影响值，肥料生产中所需要大量煤炭等化石燃料的消耗，与有机肥相比，合成化肥需要非生物资源耗竭量最大，会排放了较多的碳氢化合物和氟烃类化合物，易造成了臭氧层消耗[34]。

就其发展轨迹来看，我国高等职业教育具有非常明显的“高等性”与“职业性”，与普通高校存在着较大差别，高职院校的校园文化建设，不能固守中等职业教育的文化塑造模式，也不能一味模仿和全搬全套普通高校的文化建设。作为辅导员，一方面，要充分认识到高职校园文化的“高等性”，重视优秀中华传统文化、社会主义核心价值体系，公民日常行为规范等的重要作用，追求符合规范，教会学生做人，突出人文导向；另一方面，要充分认识到高职校园文化的“职业性”，重视效益至上理念、行业企业文化、优秀企业标准、精神等的重要作用，追求符合专业，教会学生做事，突出就业导向。

  

图1　寒地稻米生命周期系统边界Fig.1　System boundaries for rice production in cold region

1.4　数据来源与清单分析

研究稻米生产的环境影响的初级活动数据包括田间各项作业活动，加工贮存各项作业活动，深加工与产品包装以及田间温室气体散逸等活动水平数据。其中，稻米生产中的物料投入来自调查数据，为了数据可靠性和完整性以及随后产品生命周期的比较分析需要，重复计算了3次，以便基础数据进行初步分析、筛选、计算和修正，见表2。农资生产的物耗、能耗和运输系统的数据采用GaBi 6.0数据库和E－coinvent数据中中国来源数据。

 

表2　1吨寒地稻米生产的生命周期主要物料投入清单Table 2　The main materials inputs during 1 ton rice production in cold region

  

项目电能柴油化肥农药单位kWh L kg kg投入量146.36 13.47 72.37 1.55

1.5　数据简化及缺失数据处理

研究稻米主要考虑系统边界内消耗量较大的原材料，对于用量较小的原材料且对环境影响较低的物质流（如生产设备所用润滑油等）、工厂设备制造以及人员活动等方面可不予考虑[28]。稻米生产中的所需农资原料的开采、收割后稻秆剩余物还田或燃烧过程气体及固体排放均为缺失数据，因此研究暂不考虑。

2　结果与分析

2.1　特征化

假定每一种环境影响类型总值为100%，通过每一个生命周期阶段在各环境影响类型中所占比例可知各阶段对每一种环境影响类型的贡献率相对值，寒地稻米生产的生命周期的农资生产、水稻种植、稻米加工和运输阶段4个子系统相对贡献率如图2所示。

提出使用BIF特征以及加入正则化和核化相关成分分析的方法实现跨年龄人脸识别。对BIF特征采用PCA降维后，采用KR-RCA实现跨年龄人脸识别。

 

表3　1吨寒地稻米生产的生命周期特征化结果Table 3　Characteristic results of LCA for 1 ton rice production in cold region

  

环境影响类型ADP GWP ODP AP EP POCP单位kg Sbeqkg CO2 eqkg R11eqkg SO2 eqkg PO4-3eqkg C2H2 eq化肥生产8.33 183.39 3.55×10-53.44×10-11.07 5.67×10-3除草剂生产0.51 8.74 1.40×10-68.39×10-31.47×10-29.72×10-5耕种2.52 52.42 1.07×10-51.56×10-12.28×10-11.89×10-3灌溉0.39 26.21 3.90×10-85.03×10-42.94×10-31.22×10-4智能育苗0.74 78.62 3.10×10-66.21×10-23.09×10-11.13×10-2田间释放0 1 496 7.30×10-86.58×10-15.32 2.41×10-1农资运输1.76 58.97 4.50×10-65.37×10-28.83×10-28.10×10-4收粮5.61 135.41 2.11×10-52.20×10-11.69×10-24.59×10-3干燥贮存2.11 17.47 8.32×10-64.87×10-26.62×10-21.35×10-3碾米包装6.37 126.67 1.09×10-51.28×10-18.83×10-23.24×10-3合计28.34 2 184 9.56×10-51.68 7.36 0.27

2.2　归一化

根据ISO14040[26]，世界环境负荷总当量数是指全世界在1年时间内消耗的资源、能源总量和污染物排放量。通过GaBi软件中2015年全世界环境负荷总当量数归一化标准值，对特征化结果进行标准化处理后，可得到寒地稻米生产过程中各个环境负荷归一化结果（表4）。

 

表4　1吨寒地稻米生产的生命周期归一化结果Table 4　Normalization results of LCA for 1 ton rice production in cold region

  

各环节化肥生产除草剂生产耕种灌溉智能育苗田间释放农资运输收粮干燥贮存碾米包装ADP 2.19×10-141.34×10-156.63×10-151.03×10-151.95×10-150 4.63×10-151.48×10-145.55×10-151.68×10-14GWP 4.38×10-122.09×10-131.25×10-126.26×10-131.88×10-123.58×10-111.41×10-123.24×10-124.18×10-133.03×10-12ODP 1.57×10-136.18×10-154.72×10-141.72×10-161.37×10-143.22×10-161.99×10-149.31×10-143.67×10-144.81×10-14AP 1.44×10-123.51×10-146.53×10-132.11×10-152.60×10-132.75×10-122.25×10-139.20×10-132.04×10-135.34×10-13EP 6.79×10-129.30×10-141.44×10-121.86×10-141.95×10-123.36×10-115.58×10-131.07×10-124.18×10-135.58×10-13POCP 1.54×10-132.64×10-155.13×10-143.30×10-153.08×10-136.54×10-122.20×10-141.25×10-133.67×10-148.80×10-14

2.3　加权

为了比较寒地稻米生产的各环节中不同环境影响类型的相对大小，需对归一化结果进行加权处理（表1），得到总体环境的最终量化结果，如表5所示。

 

表5　1吨寒地稻米生产的生命周期加权结果Table 5　Weighted results of LCA for 1 ton rice production in cold region

  

各环节化肥生产除草剂生产耕种灌溉智能育苗田间释放农资运输收粮干燥贮存碾米包装ADP 3.29×10-142.01×10-159.95×10-151.54×10-152.92×10-150 6.95×10-152.21×10-148.33×10-152.51×10-14GWP 4.38×10-112.09×10-121.25×10-116.26×10-121.88×10-113.58×10-101.41×10-113.24×10-114.18×10-123.03×10-11ODP 4.70×10-131.85×10-141.42×10-135.16×10-164.10×10-149.67×10-165.96×10-141.40×10-131.10×10-131.44×10-13AP 2.88×10-127.02×10-141.31×10-124.21×10-155.20×10-135.50×10-124.49×10-131.84×10-124.07×10-131.07×10-12EP 4.75×10-116.51×10-131.01×10-111.30×10-131.37×10-112.35×10-103.90×10-127.48×10-122.93×10-123.90×10-12POCP 4.62×10-137.92×10-151.54×10-139.90×10-159.24×10-131.96×10-116.60×10-143.74×10-131.10×10-132.64×10-13

3　讨论

3.1　子系统对各环境类型的影响

在生命周期评价中，特征化是利用不同影响类型的参数结果来共同体现寒地稻米生产系统的生命周期评价特征，它是反映对环境影响绝对量大小的参数，即绝对总量。利用GaBi 6.0软件计算得出寒地稻米生产生命周期的特征化结果，见表3。

由图可知，在寒地稻米生产中，对于GWP、AP、EP、POCP，水稻种植子系统的贡献值最大，其次是农资生产子系统，稻米加工子系统对环境影响最小，这与前人研究相一致[14,18,29].对于ADP和ODP，农资生产子系统贡献值最大，其次是稻米加工子系统，水稻种植子系统对环境影响最小。

其中，对于ADP，农资生产、水稻种植、稻米加工和运输阶段各占31.2%，12.9%、29.9%和26.0%；对于GWP，农资生产、水稻种植、稻米加工和运输阶段各占8.8%，75.7%、6.6%和8.9%；对于ODP，农资生产、水稻种植、稻米加工和运输阶段各占38.4%，14.5%、20.1%和27.0%；对于AP，农资生产、水稻种植、稻米加工和运输阶段各占21.0%，52.2%、10.6%和16.3%；对于EP，农资生产、水稻种植、稻米加工和运输阶段各占 14.8%，79.6%、2.1%和 3.5%；对于 POCP，农资生产、水稻种植、稻米加工和运输阶段各占2.1%，94.1%、1.7%和2.0%。

  

图2　寒地稻米生产过程中子系统对环境影响贡献值Fig.2　Relative contributions to each impact category for four subsystems from rice production in cold region

3.2　生产环节对各环境类型的影响

如图3所示，稻米生产中各环节的环境影响加权值。从图中可以看出，寒地稻米生产生命周期环境影响潜力大小依次是 GWP、EP、POCP、AP、ODP 和ADP。GWP是寒地稻米生产环境影响的主要类型，占总环境影响值的58.99%；EP次之，占总环境值的36.79%；AP、ODP和ADP总的环境影响值占总值的4.22%。

研究利用GaBi 6.0软件建立具体生产模型，采用常用的CML2001方法[27]，对稻米生产各生命周期的环境负荷清单数据进行特征化、标准化、归一化和加权后得到所对应的环境影响值。

另外，AEO可以渗入晶圆，托起污染颗粒，与晶圆隔离，从而将之去除[10]。因此，AEO对物理吸附造成的颗粒残留也有明显的去除作用。但过多的 AEO却会明显地增加缺陷。这一方面是因为 FAOA具有稳定泡沫的特性，无论是AEO过量还是FAOA过量，都会使抛光液中的泡沫高度增大，并且难以消除；另一方面是因为过多的AEO覆盖在铜表面会降低FAOI与双氧水对铜的化学作用，令铜的去除速率下降[11]。

  

图3　寒地稻米生产的环境影响加权值Fig.3　Weighted value of environment impact of rice production in cold region

对于寒地稻米生产10个环节（化肥生产、除草剂生产、耕种、灌溉、智能育苗、田间释放、农资运输、收粮、干燥贮存、碾米包装），建立环境影响柱状图，结果如图4所示。

由图4可知，寒地稻米生命周期各个环节的环境影响大小依次为田间释放、化肥生产、碾米包装、智能育苗、耕种、收粮、干燥贮存、除草剂生产和农资运输、田间释放和化肥生产的环境影响最为严重，分别占寒地稻米生产的生命周期总环境影响的69.83%和10.76%，其主要影响因素为GWP和EP；收粮、碾米包装和智能育苗环节次之，分别占总值的4.79%、4.03%和3.84%；而其他各个环节对环境总的影响比例为6.75%，其中除草剂生产影响值最小。

  

图4　寒地稻米生产各环节的环境影响总值Fig.4　The total value of environmental impact of each process of rice production in cold region

  

图5　寒地稻米生产各环节的环境影响贡献值Fig.5　Contributions to each process of rice production in cold region

研究以寒地稻米为研究对象，平均年产9.5 t·hm-2，功能单位为1 t稻米生命周期内所有物资投入，产出过程的环境影响。水稻种植原料主要有种子、化肥、农药等，种植过程中气体排放，农机使用，稻米加工贮藏，废弃物回收利用，产品销售，消费者使用，再生利用。基于不确定因素，包装好的精选稻米后续过程、消费者垃圾和包装垃圾管理阶段不作考虑。

3.3　敏感度分析

研究为了判断对LCA结果质量影响程度的关键清单数据，需对清单数据进行敏感度分析。根据敏感度分析方法，对寒地稻米生产的各项原材料及能源输入进行敏感度分析，结果见表6。

根据表6所示，对于GWP和ADP结果的敏感度分析，电能投入和柴油投入是对寒地稻米LCA结果影响最大的关键清单数据，其次是化肥投入，因此要减少寒地稻米生产中对GWP和ADP的影响，着重需要减少能源、电能和化肥的输入量。电能的投入是输入清单中影响GWP、AP、ADP、EP和POCP指标的最关键清单数据，同时柴油投入是输入清单中影响ADP和AP的关键清单数据，这说明电能和柴油的输入会对多个环境指标产生明显的影响，要减少寒地稻米生产的环境影响，需要优化寒地稻米生产的农艺参数，优化水稻品种[35]，优化栽培方式，如免耕[36]、旱播[37]、旱平垄作双侧双深高效栽培[38]等，着重减少用水量、化肥施用量和电能输入量。

针对国外媒体认为中国共产党开展群众路线活动更多的是基于政治经济层面的考量这一观点，剖析其研究方法我们就会发现，这一观点走入了一个大部分研究中国共产党的国外学者都存在的误区，那就是用世界政党发展的一般规律来看待中国共产党的发展。诚然，世界上大部分执政党进行自身改革的最大动因必然存在于政治、经济方面。但是，由这一观点出发的国外学者往往忽略了一个问题，那就是中国共产党并不是一个一般性政党，它当然有着一般政党具有的特点，但是其特殊性也是显而易见的。因此，如果忽略了这一点，仅仅认为群众路线教育实践活动是中国共产党在面临政治经济困局的情况下的应时之作，这样得出的结论就难免会流于主观，缺乏说服力。

综上所述，高校档案馆只有拥有丰富的馆藏资源，把档案利用的触角延伸到高校的每个院系，让档案服务覆盖到全校各部门，让师生认识到高校档案馆是学校蕴藏着宝藏的地方，才能更好地为学校的发展做出更大贡献。

 

表6　各环境影响指标对应的原料、能源、电力的敏感度分析Table 6　Sensitivity analysis of environmental impacts on raw materials，energy，and electricity

  

项目投入化肥投入除草剂投入电能投入柴油投入农资运输收粮干燥贮存碾米包装ADP 0.007 194 0.001 565 0.049 323 0.088 778 6.95×10-152.21×10-148.33×10-152.51×10-14GWP 0.078 379 0.027 360 0.729 293 0.092 744 1.41×10-113.24×10-114.18×10-123.03×10-11敏感度ODP 0.000 000 0.000 000 0.000 001 0.000 000 5.96×10-141.40×10-131.10×10-131.44×10-13AP 0.000 928 0.000 594 0.003 688 0.002 067 4.49×10-131.84×10-124.07×10-131.07×10-12EP 0.000 928 0.000 594 0.003 688 0.002 067 3.90×10-127.48×10-122.93×10-123.90×10-12POCP 0.000 005 0.000 004 0.000 086 0.000 002 6.60×10-143.74×10-131.10×10-132.64×10-13

4　结论

（1）在1 t寒地稻米的生产中，6种环境影响类型 ADP、GWP、ODP、AP、EP 和 POCP 潜值分别为1.12 ×10-13、5.22 ×10-10、1.27 ×10-12、1.40 ×10-11、3.26 ×10-10、2.20×10-11。

家长是孩子的启蒙老师，是孩子的第一任教育者。温馨的家庭氛围和良好的家庭文化将对留守儿童产生巨大的影响。但是父母长期在外，导致家长角色弱化，父母与孩子产生隔阂，儿童应有的心灵归属缺失。

（2）GWP和EP是寒地稻米生产中环境影响的主要类型。在4个子系统中，对于GWP、AP、EP、POCP，水稻种植子系统的贡献值最大；对于ADP和ODP，农资生产子系统贡献值最大。

（3）生产1 t寒地稻米，化肥生产、除草剂生产、耕种、灌溉、智能育苗、田间释放、农资运输、收粮、干燥贮存、碾米包装的影响总值分别为9.52×10-11、2.84 ×10-12、2.42 ×10-11、6.41 ×10-12、3.39 ×10-11、6.18 ×10-10、1.86 ×10-11、4.24 ×10-11、7.74 ×10-12、3.57 ×10-11、8.85×10-10。

护理安全是在实施护理服务全过程中，患者不发生法律和法定的规章制度允许范围以外的心理人体结构或功能上的损害、障碍或死亡。无论是在任何年代，护理安全问题都是学术界和社会共同关心的问题，护理工作必须以保障患者的基本安全问前提进行，这不仅是道德上的要求，更是护理工作的职业道德要求。在这样的背景下，有关部门和护理专业本身，加强对护理工作的管理，就显得十分必要。

（4）田间释放和化肥生产是寒地稻米生产中环境影响的主要环节。由于煤炭燃烧产生的污染物排放，田间释放环节的GWP、AP、EP和POCP的影响值为寒地稻米生产各环节的最高值。同时，化肥生产阶段具有最高的ADP和ODP影响值。

沿着network/cohesion/density计算网络密度。从表1可以看出，2010-2016年河南省旅游经济网络密度从0.26上升为0.54，意味着河南省旅游经济网络结构发育程度不断提高，城市之间的旅游经济的互动性增强。但是从理论上看，有18个城市节点组成的河南省旅游经济网络最大可能的关系数应为324条，然而2010年实际存在的联系数仅为40条，虽然逐年增加，但直到2016年实际存在的关系数也只有82条。可见河南省城市旅游经济网络结构的整体发育程度偏低，网络结构尚不健全，城市间旅游经济总体互动性不强，旅游经济联系不紧密。

参考文献：

［1］ 黑龙江省统计局.黑龙江省食品主导产业发展状况［EB/OL］.（2015-03-18）［2015.12.18］.http://www.hlj.stats.gov.cn/tjfx/sjtjfx/201503/t20150318_31088.htm.

［2］ 邓南圣，王小兵.生命周期评价［M］.北京：化学工业出版社，2003.

［3］　Audsley A，Alber S，Clift R，et al.Harmonisation of environmental life cycle assessment for agriculture.Final report，Concerted Action AIR3-CT94-2028 ［R］.Brussels：European Commission DG VI，1997.

［4］ Williams A G，Audsley E，Jones R，et al.Determining the environmental burdens and resource use in the production of agricultural and horticultural commodities.Main Report，Defra Research Project IS0205 ［R］.Bedford：Cranfield University and Defra，2006.

［5］ Majumdar D.Methane and nitrous oxide emission from irrigated rice fields：Proposed mitigation strategies［J］.Current Science，2003，84（10）：1317-1326.

［6］　Sahrawat K L.Nitrification inhibitors for controlling methane emission from submerged rice soils［J］.Current Science，2004，87（8）：1084-1087.

［7］　Leip A，Bocchi S.Measurements of greenhouse gas emissions from rice cultivation in Italy［C］//Idiap.Novara：Procedings of the Fourth Temperate Rice Conference，2007:25-28.

［8］ Neue H U.Fluxes of methane from rice fields and potential for mitigation ［J］.Soil Use and Management，1997，13（54）：258-267.

［9］ Blengini G A，Busto M.The life cycle of rice：LCA of alternative agri-food chain management systems in Vercelli（Italy）［J］.JournalofEnvironmentalManagement，2009，90（3）：1512-1522.

［10］ Kasmaprapruet S，Paengjuntuek W，Saikhwan P，et al.Life cycle assessment of milled rice production：Case study in Thailand ［J］.European journal of scientific research，2009，30（2）：195-203.

［11］ Hokazono S，Hayashi K.Variability in environmental impacts during conversion from conventional to organic farming：a comparison among three rice production systems in Japan ［J］.Journal of Cleaner Production，2012，28：101-112.

［12］　Mohammadi A，Rafiee S，Jafari A，et al.Joint life cycle assessmentand　data　envelopmentanalysis　forthe benchmarking of environmental impacts in rice paddy production［J］.Journal of Cleaner Production，2015，106：521-532.

［13］ Breiling M，Hashimoto S，Sato Y，et al.Rice-related greenhouse gases in Japan，variations in scale and time and significance for the Kyoto Protocol［J］.Paddy and Water Environment，2005，3（1）：39-46.

［14］　Harada H，Kobayashi H，Shindo H.Reduction in greenhouse gas emissions by no-tilling rice cultivation in Hachirogata polder，northern Japan：Life-cycle inventory analysis［J］.Soil Science and Plant Nutrition，2007，53（5）：668-677.

［15］ Roy P，Shimizu N，Okadome H，et al.Life cycle of rice：challenges and choices for Bangladesh ［J］.Journal of Food Engineering，2007，79（4）：1250-1255.

［16］ 董珑丽，李小平，姜雪，等.上海崇明岛水稻生产能耗与碳足迹生命周期评价［J］.农业环境科学学报，2014，33（6）：1254-1260.

［17］ 曹黎明，李茂柏，王新其，等.基于生命周期评价的上海市水稻生产的碳足迹［J］.生态学报，2014，34（2）：491-499.

［18］ Xue J F，Pu C，Liu S L，et al.Carbon and nitrogen footprint of double rice production in Southern China［J］.E－cological Indicators，2016，64：249-257.

［19］　Liu Q，Liu B，Ambus P，et al.Carbon footprint of rice production under biochar amendment-a case study in a Chinese rice cropping system［J］.Global Change Biology Bioenergy，2016，8（1）：148-159.

［20］ 徐小明.吉林西部水田土壤碳库时空模拟及水稻生产的碳足迹研究［D］.长春：吉林大学，2011.

［21］ 温和.黑龙江省村域农业生态系统碳平衡及低碳农业对策研究［D］.哈尔滨：东北农业大学，2011.

［22］ 王明新，夏训峰，刘建国，等.太湖地区高产水稻生命周期评价［J］.农业环境科学学报，2009，28（2）：420-424.

［23］ 梁龙，陈源泉，高旺盛.两种水稻生产方式的生命周期环境影响评价［J］.农业环境科学学报，2009，28（9）：222-226.

［24］　Yang S S，Lai C M，Chang H L，et al.Estimation of methane and nitrous oxide emissions from paddy fields in Taiwan［J］.Renewable Energy，2009，34（8）：1916-1922.

［25］ 郝宁，李红宇，王淼，等.利用AMMI模型分析寒地水稻区试品种产量的基因型与环境互作［J］.黑龙江科技信息，2015，26（27）:5-8.

［26］ ISO 14040.Environmental management-life cycle assessment-principles and framework ［S］.International Standard，2006.

［27］ Guinee J.Handbook on life cycle assessment operational guide to the ISO standards［J］.The International Journal of Life Cycle Assessment，2002，7（5）：311-313.

［28］ Silva D，Lahr F，Garcia R，et al.Life cycle assessment of medium density particleboard（MDP）produced in Brazil［J］.The International Journal of Life Cycle Assessment，2013，18（7）：1404-1411.

［29］　程竞卉，柳伟先，李燕花，等.早熟水稻新品种“延稻 1号”选育［J］.延边大学农学学报，2016（1）：40-43.

［30］ Tong L，Liu X，Liu X，et al.Life cycle assessment of water reuse systems in an industrial park［J］.Journal of Environmental Management，2013，129（18）：471-478.

［31］ Wang M，Xia X，Zhang Q，et al.Life cycle assessment of a rice production system in Taihu region，China［J］.International Journal of Sustainable Development and World Ecology，2010，17（2）：157-161.

［32］ 黄树辉，蒋文伟，吕军，等.氮肥和磷肥对稻田N2O排放的影响［J］.中国环境科学，2005，25（5）：540-543.

［33］ 曹娜，王睿，廖婷婷，等.厌氧条件下砂壤水稻土N2、N2O、NO、CO2和 CH4排放特征［J］.环境科学，2015（9）：3373-3382.

［34］ Billen P，Costa J，Van der Aa L，et al.Electricity from poultry manure：a cleaner alternative to direct land application［J］.Journal of Cleaner Production，2015，96：467-475.

［35］ 姜树坤，张凤鸣，白良明，等.黑龙江省第一积温带水稻育种趋势分析［J］.作物杂志，2015（4）：61-64.

［36］ 钱明媚，肖永良，彭文涛，等.免耕水稻土固定CO2自养微生物多样性［J］.中国环境科学，2015，35（12）：3754-3761.

［37］ 邹建文，秦艳梅，刘树伟.不同水分管理方式下水稻生长季N2O排放量估算：模型建立［J］.环境科学，2009，30（2）：313-321.

［38］ 郑桂萍，韩东来，吕艳东，等.水稻旱平垄作双侧双深高效栽培新模式与技术规范［J］.黑龙江八一农垦大学学报，2016，28（1）：1-3.