深圳市生态系统生产总值核算技术规范(DB4403-T 141-2021).pdf

ICS 13.020.20 CCS Z 04 DB4403 深圳市地方标准 DB4403/T 1412021 深圳市生态系统生产总值核算 技术规范 Technical specification for accounting gross ecosystem product of Shenzhen 2021-02-23 发布 2021-03-01 实施 深圳市市场监督管理局 发 布 DB4403/T 1412021 I 目 次 前言.II 引言.III 1 范围.1 2 规范性引用文件.1 3 术语和定义.1 4 总则.2 5 评估指标体系.3 6 调查与评估使用数据要求.4 7 评估使用的主要方法.5 附录 A（规范性）生态系统生产总值功能量核算方法.7 附录 B（规范性）生态系统生产总值价值量核算方法.16 附录 C（规范性）生态系统生产总值核算方法.17 附录 D（资料性）生态系统生产总值核算参数补充.18 参考文献.35 DB4403/T 1412021 II 前 言 本文件按照 GB/T 1.12020标准化工作导则 第 1 部分：标准化文件的结构和起草规则的规定起草。本文件由深圳市生态环境局提出并归口。本文件起草单位：中国科学院生态环境研究中心、深圳市环境科学研究院、中国标准化研究院。本文件主要起草人：欧阳志云、韩宝龙、孙芳芳、李思诗、潘雪莲、陈礼、肖燚、宋昌素、束承继、从睿、邓瑞民、黄爱兵、付岚、陈龙、董笑语。DB4403/T 1412021 III 引 言 为全面践行习近平生态文明思想和“绿水青山就是金山银山”理念，深入落实中共中央国务院关于支持深圳建设中国特色社会主义先行示范区的意见中对“探索实施生态系统服务价值核算制度”的重要指示，促进深圳城市社会经济与生态环境可持续发展，支撑深圳生态系统管理，制定本文件。本文件可为将生态效益纳入经济社会发展评价体系、完善发展成果考核评价体系提供重要支撑，为建立生态产品价值实现机制、区域生态补偿、自然资源资产审计、自然资产负债表编制等制度的实施提供科学依据。DB4403/T 1412021 1 深圳市生态系统生产总值核算技术规范 1 范围 本文件规定了深圳市生态系统生产总值核算的总则、评估指标体系、调查与评估使用数据要求、评估使用的主要方法等内容。本文件适用于深圳市和各区（新区、深汕特别合作区）的生态系统服务价值核算、深圳市生态产品价值核算等与深圳市生态系统生产总值核算内涵相同的工作。2 规范性引用文件 下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB 30952012 环境空气质量标准 GB 38382002 地表水环境质量标准 3 术语和定义 下列术语和定义适用于本文件。3.1 生态系统 ecosystem 地球表面生物及其环境通过能流、物流、信息流形成的功能整体。注：通常包括森林生态系统、草地生态系统、湿地生态系统、荒漠生态系统、农田生态系统、城市生态系统等类型。3.2 生态系统服务 ecosystem services 人类从生态系统中得到的惠益。注：通常包括物质产品、调节服务、文化旅游服务以及支持服务。支持服务属于中间过程服务，不属于最终服务，为了避免重复计算，在生态系统生产总值中不考虑支持服务。3.3 生态系统生产总值 gross ecosystem product,GEP 生态系统为人类福祉和经济社会可持续发展提供的各种最终产品与服务价值的总和。注：通常又称生态系统服务价值、生态产品价值，主要包括生态系统提供的物质产品、调节服务和文化旅游服务的价值。3.4 物质产品 provisioning services DB4403/T 1412021 2 人类从生态系统获取的可在市场交换的食物、纤维、燃料、遗传资源、天然药物、装饰材料与其他物质材料等。3.5 调节服务 regulating services 生态系统调节作用提供的调节气候、涵养水源、保持土壤、调蓄洪水、降解污染物、固碳等改善人类生存与生活环境的惠益。3.6 文化旅游服务 cultural&tourism services 人类通过精神感受、知识获取、休闲娱乐和美学体验从生态系统获得的非物质惠益。3.7 功能量 product amount 生态系统产品与服务的物理量。注：功能量可以是粮食产量、洪水调蓄量、土壤保持量、固碳量与景点旅游人数等。3.8 价值量 product value 生态系统产品与服务的货币价值。4 总则 4.1 生态系统生产总值核算目的 根据各类生态模型法定量评估生态系统生产总值的功能量，借助价格将不同计量单位的生态系统产品和服务的功能量货币化得到生态系统生产总值。生态系统生产总值可以为生态效益纳入经济社会发展评价体系、完善发展成果考核评价体系提供重要支撑，为区域生态补偿、自然资源资产审计等制度的制定提供科学依据，为评估区域生态资产及其变化状况提供科学方法。4.2 生态系统生产总值指标选择原则 4.2.1 来自自然生态过程 生态系统产品（含服务）应由自然生态过程产生，不以自然生态过程为主的人工生态环境设施产生的生态环境改善产品（含服务）不作为生态系统生产总值指标。4.2.2 对人类有益 生态系统生产总值核算中仅考虑对人类产生惠宜的生态系统服务。4.2.3 非生态系统服务中间过程 生物多样性维持、土壤保肥等中间过程服务不计入生态系统生产总值指标体系。4.2.4 技术可行数据可获取 DB4403/T 1412021 3 本文件选取的指标为当前技术成熟、数据可获取的指标，部分生态系统服务由于缺乏可行的价值化模型和可获取的关键数据，不在本文件考虑之内。5 评估指标体系 深圳市生态系统生产总值评估主要考虑的内容如表1所示。表1 深圳生态系统生产总值评估指标体系 一级指标一级指标 二级指标二级指标 指标说明指标说明 物质产品物质产品 农林牧渔产品农林牧渔产品 能够在市场交易，并由本地生态系统生产的一次农、林、牧、渔业产品。生态能源生态能源*在市场交易，并由本地生态系统生产的生态能源。水资源水资源*能够在市场交易，并由本地供给的水资源。调节服务调节服务 减少泥沙淤积减少泥沙淤积 本地生态系统通过林冠层、枯落物、根系等各个层次保护土壤、消减降雨侵蚀力，增加土壤抗蚀性，减少泥沙阻塞河道。减少面源污染减少面源污染 本地生态系统通过保持土壤，减少氮、磷等土壤营养物质进入下游水体（包括河流、湖泊、水库和海湾等），可降低下游水体的面源污染。调节气候调节气候 本地陆域生态系统蒸腾、蒸发带走的热量。固定二氧化碳固定二氧化碳 本地生态系统吸收大气中的二氧化碳合成有机质，将碳固定在植物或土壤中。削减洪涝削减洪涝 本地生态系统吸纳的降水、蓄积的径流和过境水。涵养水源涵养水源 由本地生态系统拦截滞蓄的降水，通过增强土壤下渗、蓄积，和补充地下水，净增加的本地水资源量。削减交通噪声削减交通噪声 本地道路两侧绿化带对噪声的削减。防护海岸带防护海岸带*本地海岸带生态系统如红树林、滨海盐沼等减低海浪，避免或减小海堤或海岸侵蚀的功能。净化空气净化空气 本地生态系统吸收、过滤、阻隔和分解大气污染物，改善大气环境（二氧化硫、氮氧化物、颗粒污染物）。净化水体净化水体 湖泊、河流、沼泽等水域湿地生态系统吸附、降解和转化水体污染物（COD、氨氮、总磷）。文化旅游服务文化旅游服务 旅游休闲服务旅游休闲服务 本地生态系统提供的文娱休闲服务，有益于丰富知识、愉悦身心。自然景观溢价自然景观溢价*由于本地生态系统的美丽自然景观使房产交易（转让、出租）过程中发生了溢价。康养服务康养服务*本地生态空间提供疗养、康复（如减少呼吸道疾病就医量）等有益健康的功能。注：表中*表示可选项，根据核算区生态系统特征和生态产品供给情况进行核算。DB4403/T 1412021 4 6 调查与评估使用数据要求 生态系统生产总值评估使用的数据主要包括遥感数据、监测数据、统计数据、实地调查数据以及参考文献数据等。关键结果数据的单位、时间等要求如表2所示。表2 深圳生态系统生产总值评估数据要求 编号编号 生态系统产品（服务）类型生态系统产品（服务）类型 功能量单位功能量单位 价值量单位价值量单位 时间时间 0 生态系统生产总值 亿元 评估时段 1 物质产品 亿元 评估时段 11 农林牧渔 亿元 亿元 评估时段 12 生态能源*千瓦时 亿元 评估时段 13 水资源*万立方米 亿元 评估时段 2 调节服务 亿元 评估时段 21 涵养水源 万立方米 亿元 评估时段 22 削减洪涝 万立方米 亿元 评估时段 23 调节气候 亿千瓦时 亿元 评估时段 24 固定二氧化碳 万吨 亿元 评估时段 25 净化空气 吨 亿元 评估时段 251 净化空气二氧化硫 吨 亿元 评估时段 252 净化空气氮氧化物 吨 亿元 评估时段 253 净化空气颗粒污染物 吨 亿元 评估时段 26 净化水体 吨 亿元 评估时段 261 净化水体氨氮 吨 亿元 评估时段 262 净化水体总磷 吨 亿元 评估时段 263 净化水体 COD 吨 亿元 评估时段 27 减少泥沙淤积 万立方米 亿元 评估时段 28 减少面源污染 万吨 亿元 评估时段 281 减少面源污染总氮 万吨 亿元 评估时段 282 减少面源污染总磷 万吨 亿元 评估时段 29 削减交通噪声 平均分贝 亿元 评估时段 210 防护海岸带*千米 亿元 评估时段 3 文化旅游服务 亿元 评估时段 31 旅游休闲服务 亿元 亿元 评估时段 32 自然景观溢价*亿元 亿元 评估时段 33 康养服务*人次 亿元 评估时段 注：表中*表示可选项，根据核算区生态系统特征和生态产品供给情况进行核算。DB4403/T 1412021 5 7 评估使用的主要方法 7.1 生态系统生产总值功能量核算方法 生态系统生产总值功能量核算包括三大类，即物质产品功能量核算、调节服务功能量核算、文化旅游服务功能量核算。生态系统生产总值功能量的核算项目、功能量指标和核算方法如表 3 所示。具体方法应符合附录 A 的规定。表3 生态系统生产总值功能量核算方法 一级指标一级指标 二级指标二级指标 功能量功能量 核算方法概述核算方法概述 物质产品物质产品 农林牧渔产品农林牧渔产品 初级农林牧渔产品产值 查阅当地统计年鉴或农业部门数据 生态能源生态能源*各类生态能源生产数量 查阅当地统计年鉴或能源部门数据 水资源水资源*本地自然水体供水量 总供水量中，由当地自然水资源供给的量；查阅当地水资源公报或水务部门数据。调节服务调节服务 减少泥沙淤积减少泥沙淤积 减少泥沙淤积量 在产流降雨条件下，由通用土壤流失方程计算得出土壤保持量，再乘以泥沙形成系数。减少面源污染减少面源污染 减少各类面源污染物量 根据土壤保持量和土壤中氮、磷的含量，计算因自然生态系统的土壤保持功能减少面源污染的量。调节气候调节气候 消耗的热量 在高于适宜温度时期，本地各类生态系统单位面积蒸散发消耗热量乘以面积，并加总。固定二氧化碳固定二氧化碳 二氧化碳固定量 根据两年间的植被生物量变化计算二氧化碳固定量，或根据不同生态系统类型固碳速率和面积计算。削减洪涝削减洪涝 区域调蓄洪水量 在城市范围内利用 SCS 模型计算指标削减径流量；以及利用监测数据计算湖泊、水库的滞留水量。涵养水源涵养水源 水源涵养量 本地降雨量减去径流量，再减去蒸散发量。削减交通噪声削减交通噪声 路侧绿化带平均消减噪声分贝量 根据在不同路段的典型样地监测数据，评估道路绿化（两侧及内部）平均消减噪声量。防护海岸带防护海岸带*达到安全岸线强度的自然岸线长度 根据 InVEST 风暴潮风险评估模型计算达到安全岸线强度的自然岸线长度。净化空气净化空气 净化各类大气污染物量 根据本地大气污染物达标水平，选择污染物排放量或者净化量（每类生态系统单位面积净化量乘面积）为实物量。净化水体净化水体 净化各类水体污染物量 根据本地水质达标水平，选择水体污染物排放量或者净化量（每类生态系统单位面积净化量乘面积）为实物量。文化旅游文化旅游服务服务 旅游休闲服务旅游休闲服务 旅游休闲人次与时间 根据抽样调查统计获取的自然风景旅游与休闲人数和平均滞留时间。自然景观溢价自然景观溢价*景观溢价价值 根据抽样调查统计获取当年房屋交易中的景观溢价价值，以及酒店交易中的景观溢价价值。康养服务康养服务*减少呼吸道疾病就医量和死亡量 根据中国城市大气污染健康终端效应时间序列的 Meta 分析和 WHO 的“污染物浓度死亡风险”变化曲线，计算暴露人口的变化量。注：表中*表示可选项，根据核算区生态系统特征和生态产品供给情况进行核算。DB4403/T 1412021 6 7.2 生态系统生产总值价值量核算方法 根据各服务功能指标的内涵，价值量的核算方法如表 4 所示，具体方法应符合附录 B 的规定。在所有类型的生态系统产品与服务功能量、价值量核算的基础上，将各类生态系统产品与服务价值加总，得到生态系统生产总值（GEP），具体方法应符合附录 C 的规定。鉴于价值量核算定价方法的单价变动性，建议定价比选方法和推荐定价方案见附录 D。表4 生态系统生产总值价值量核算方法 一级指标一级指标 二级指标二级指标 核算内容核算内容 价值量核算方法概述价值量核算方法概述 物质产品物质产品 农林牧渔产品农林牧渔产品 初级农林牧渔产品价值 根据深圳统计年鉴数据核算 生态能源生态能源*生态能源价值 根据深圳统计年鉴数据核算 水资源水资源*水资源价值 原水价格 调节服务调节服务 减少泥沙淤积减少泥沙淤积 减少泥沙淤积价值 土方清运成本 减少面源污染减少面源污染 减少面源氮价值 治理成本 减少面源磷价值 治理成本 调节气候调节气候 植被蒸腾调节温度价值 普通居民合表用户用电成本 水面蒸发调节温度价值 固定二氧化碳固定二氧化碳 固碳价值 碳交易价格 削减洪涝削减洪涝 城区植被暴雨径流调节价值 海绵城市蓄水池建设成本 郊野植被暴雨径流调节价值 水库单位库容造价 湖泊洪水调蓄价值 水库洪水调蓄价值 涵养水源涵养水源 水源涵养价值 水库单位库容造价 削减交通噪声削减交通噪声 交通噪声消减价值 人工降噪幕墙建设成本 防护海岸带防护海岸带*海岸带防护价值 人工岸线建设成本 净化空气净化空气 净化二氧化硫价值 治理成本 净化氮氧化物价值 治理成本 净化颗粒污染物价值 治理成本 净化水体净化水体 净化 COD 价值 治理成本 净化总氮价值 治理成本 净化总磷价值 治理成本 文化旅游服务文化旅游服务 旅游休闲服务旅游休闲服务 旅游休闲服务价值 根据调查问卷 自然景观溢价自然景观溢价*房产产权销售溢价 销售数据的拟合公式模拟 酒店销售溢价 酒店销售数据调查 康养服务康养服务*减少呼吸道疾病价值 门诊、住院成本 人力资本替代成本 注：表中*表示可选项，根据核算区生态系统特征和生态产品供给情况进行核算。DB4403/T 1412021 7 附 录 A（规范性）生态系统生产总值功能量核算方法 A.1 物质产品功能量核算方法 各类物质产品的产量通过统计资料获取。物质产品总产量按式（A.1）进行求和。Epro=1（A.1）式中：Epro 物质产品总产量（t/a）；i 核算区物质产品种类，i=1，2，3，n；Ei 第 i 种物质产品的产量（t/a）。注 1：物质产品包括自然生态系统野生的和人工集约化种养殖的农业产品、林业产品、畜牧业产品、渔业产品以及水资源与生态能源。注 2：产量数据来自林业、农业、渔业及统计部门相关资料，以及实地调查。注 3：生态能源功能量单位为 kWh，水资源功能量单位为吨。A.2 调节服务功能量核算方法 A.2.1 减少泥沙淤积服务 减少泥沙淤积服务的功能量核算按式（A.2）和式（A.3）。Qsd=(Qsr/)（A.2）Qsr=R K L S 1 C P （A.3）式中：Qsd 减少泥沙淤积量（m/a）；泥沙淤积系数；Qsr 土壤保持量（t/a）；土壤容重（t/m）；R 降雨侵蚀力因子（MJmmhm-2h-1a-1），用多年平均年降雨侵蚀力指数表示，核算见附录 D.1；K 土壤可蚀性因子（thm2hhm-2MJ-1mm-1），用标准样方上单位降雨侵蚀力所引起的土壤流失量来表示，核算见附录 D；L 坡长因子（无量纲），核算见附录 D；S 坡度因子（无量纲），核算见附录 D；DB4403/T 1412021 8 C 植被覆盖因子（无量纲），核算见附录 D；P 水土保持措施因子（无量纲），核算见附录 D。A.2.2 减少面源污染服务 减少面源污染服务的功能量核算按式（A.4）。Qdpd=Qsr=1（A.4）式中：Qdpd 减少面源污染量（t/a）；i 土壤中氮、磷等营养物质数量，i=1，2，n；Qsr 土壤保持量（t/a）（核算方法见式 A.3）；ci 土壤中氮、磷等营养物质的纯含量（%）。A.2.3 调节气候服务 采用生态系统蒸腾蒸发总消耗的能量作为调节气候的功能量，核算按式（A.5）、式（A.6）和式（A.7）。Ett=Ept+Ewe（A.5）Ept=EPPi Si D3 106/(3600)（A.6）Ewe=Ew 103/(3600)（A.7）式中：Ett 生态系统蒸腾蒸发消耗的总能量（kWh/a）；Ept 生态系统植被蒸腾消耗的能量（kWh/a）；Ewe 生态系统水面蒸发消耗的能量（kWh/a）；i 生态系统类型（无量纲）；EPPi第 i 类生态系统单位面积蒸腾消耗热量（kJm-2d-1）；Si 第 i 类生态系统面积（km2）；D 空调开放天数（天）；r 空调能效比：3.0（无量纲）；Ew 水面蒸发量（m3）；q 挥发潜热，即蒸发 1 克水所需要的热量（J/g）。DB4403/T 1412021 9 A.2.4 固定二氧化碳服务 采用生态系统二氧化碳固定量作为核算指标，核算方法有以下三种：如果净生态系统生产力（NEP,Net ecosystem productivity）数据可得，生态系统二氧化碳固定量核算按式（A.8）；QtCO2=MCO2/MC NEP（A.8）式中：QtCO2 生态系统二氧化碳固定量（tCO2/a）；MCO2/MCC 转化为 CO2的系数，44/12；NEP 净生态系统生产力（tC/a），核算见附录 D。如果净生态系统生产力（NEP）数据不可得，采用生物量法测算生态系统二氧化碳固定量，核算按式（A.9）；QtCO2=MCO2MC A CC (AGBt2 AGBt1)（A.9）式中：QtCO2 生态系统二氧化碳固定量（tCO2/a）；MCO2/MCC 转化为 CO2的系数，44/12；A 生态系统面积（ha）；CC 生物量-碳转换系数；AGBt2 第 t2年的生物量（t/ha）；AGBt1 第 t1年的生物量（t/ha）。如果净生态系统生产力（NEP）和生物量数据均不可得，采用固碳速率法计算生态系统二氧化碳固定量，核算按式（A.10）。QtCO2=MCO2MC (FCS+GSCS+WCS+CSCS)（A.10）式中：QtCO2 生态系统二氧化碳固定量（tCO2/a）；MCO2/MCC 转化为 CO2的系数，44/12；FCS 森林（及灌丛）固碳量（tC/a），核算见附录 D；DB4403/T 1412021 10 GSCS 草地固碳量（tC/a），核算见附录 D；WCS 湿地固碳量（tC/a），核算见附录 D；CSCS 农田固碳量（tC/a），核算见附录 D。A.2.5 削减洪涝服务 削减洪涝量包括林灌草削减洪涝水量和水库削减洪涝量两部分，核算按式（A.11）、式（A.12）和式（A.13）。Cfm=Cfc+Crc（A.11）Cfc=Pi Rfi Ai=1 1000（A.12）Crc=Crfm（A.13）式中：Cfm 削减洪涝量（m3/a）；Cfc 林灌草削减洪涝量（m3/a）；Crc 水库削减洪涝量（m3/a）；n 生态系统类型数（无量纲）；i 第 i 类生态系统类型（无量纲）；Pi 暴雨降雨量（mm），采用 50mm 日降雨划定暴雨日；Rfi 暴雨径流量（mm），核算见附录 D；Ai 第 i 类生态系统的面积（km2）；Crfm 水库防洪库容（m3/a），核算见附录 D。A.2.6 涵养水源服务 采用水量平衡方程，计算生态系统通过拦截滞蓄降水，增强土壤下渗、蓄积，涵养土壤水分、调节地表径流和补充地下水所增加的水资源总量，核算按式（A.14）。Qwr=Ai Pi Ri ETi 10-3=1（A.14）式中：Qwr 涵养水源量（m3/a）；n 生态系统类型总数（无量纲）；i 第 i 类生态系统类型（无量纲）；DB4403/T 1412021 11 Ai 第 i 类生态系统的面积（m2）；产流降雨量（mm/a）；Ri 地表径流量（mm/a），核算见附录 D；ETi 蒸散发量（mm/a）。A.2.7 削减交通噪声服务 将全市具有交通绿化带且非高架道路的路段按照宽度、流量、道路等级进行分类，每类道路两侧设置监测样点采集相关数据计算削减交通噪声量，核算按式（A.15）。=Ri NAini=1 Rini=1（A.15）式中：ANA 为全市交通噪音平均削减量（dB）；Ri 第 i 类道路的长度（km）；NAi 第 i 类道路两侧的平均降噪分贝（dB）。注：降噪分贝数由绿化带近路侧和远路侧噪声差值确定。A.2.8 防护海岸带服务 通过统计植被重要性指数 HRI0 的岸段数量计算海岸带的防护总长度，核算按式（A.16）。=d NHRI 0 n1（A.16）式中：D 防护海岸带长度（km）；n 植被对海岸带起防护作用的岸段数量；d 单位栅格长度（km）；N 植被重要性指数 HRI0 的岸段数量；HRI 植被重要性指数（无量纲），核算见附录 D。A.2.9 净化空气服务 采用大气污染物净化量作为核算指标，按照 GB 30952012 中对环境空气质量应控制项目的规定，选取二氧化硫、氮氧化物、颗粒污染物等污染物指标，当颗粒污染物数据不可得时可采用工业粉尘替代，核算方法有以下两种：如果污染物浓度未超过环境空气功能区质量标准，大气污染物净化量为污染物排放量，核算按DB4403/T 1412021 12 式（A.17）；Qap=Qi=1（A.17）式中：Qap 大气污染物排放总量（kg/a）；n 大气污染物类别的数量（无量纲）；i 污染物类别（无量纲）；Qi 第 i 类大气污染物排放量（kg/a）。如果污染物浓度超过环境空气功能区质量标准，大气污染物净化量为生态系统自净能力，核算按式（A.18）。Qap=Qij=1=1 Ai（A.18）式中：Qap 生态系统空气净化能力（kg/a）；m 生态系统类型的数量（无量纲）；i 生态系统类型（无量纲）；n 大气污染物类别的数量（无量纲）；j 大气污染物类别（无量纲）；Qij 第 i 类生态系统第 j 种大气污染物的单位面积净化量（kg/km2a）；Ai 第 i 类生态系统类型面积（km2）。A.2.10 净化水体服务 采用水体污染物净化量作为核算指标，按照 GB 38382002 对水环境质量应控制项目的规定，选取COD、总氮、总磷等污染物指标，核算方法有以下两种：如果污染物浓度超过地表水水域环境功能标准限值，水体污染物净化量为生态系统自净能力，核算按式（A.19）；Qwp=Qi A=1（A.19）式中：Qwp 水体污染物净化量（t/a）；n 水体污染物类别的数量（无量纲）；DB4403/T 1412021 13 i 污染物类别（无量纲）；Qi 水域对第 i 类水体污染物的单位面积年净化量（t/km2a）；A 水域面积（km2）。如果污染物排放浓度未超过地表水水域环境功能标准限值，水体污染物净化量为排放量与随水输送出境的污染物量之差，核算按式（A.20）。Qwp=(Qei+Qai Qdi)=1（A.20）式中：Qwp 污染物净化总量（kg/a）；n 水体污染物类别的数量（无量纲）；i 污染物类别（无量纲）；Qei 某种（类）污染物入境总量（kg/a）；Qai 区域内某种（类）污染物排入河流总量（kg/a）；Qdi 某种（类）污染物出境量（kg/a）。A.3 文化旅游服务核算方法 A.3.1 旅游休闲服务 采用区域内自然景观的年旅游总人次作为旅游休闲服务的功能量评价指标，核算按式（A.21）。运用旅行费用法核算人们通过休闲旅游活动体验生态系统与自然景观美学价值，并获得知识和精神愉悦的非物质价值，核算按式（A.22）、式（A.23）、式（A.24）和式（A.25）。Nt=Ntini=1（A.21）Vr=（CC+CS）Nt（A.22）CC=Nj CCj=1（A.23）Nj=nj N（A.24）CCj=Ctc,j+Clf,j+Cad,j（A.25）式中：Nt 游客总人数（人次）；i 旅游景观个数，i=1，2，n；Nti 第 i 个旅游景观的游览人数（人次）；DB4403/T 1412021 14 Vr 评估区域的旅游休闲价值（元）；CC 评估区域的消费者支出（元/人次）；CS 评估区域的消费者剩余（消费者的意愿支付与实际支付之差）（元/人次）；j 景区类型（如公园、绿道等）；Nj 来自 j 类景区的游客总人次（人次）；CCj 来自 j 类景区的游客每人次旅行成本或支付意愿（元/人次）；nj 调查样本中 j 类景区（适用非收费景区，如社区公园、绿道等）的游客人均访问次数（次/人）；N 被评估地常住人口数量（人）；Ctc,j 调查样本中 j 类景区的游客平均每人次交通费用（元/人次）；Clf,j 平均每人次食宿及其他费用（元/人次）；Cad,j 平均每人次门票费用（元/人次）。A.3.2 自然景观溢价 A.3.2.1 房屋交易自然景观溢价服务 总的房屋销售自然景观溢价核算按式（A.26）。=(PPi Pi)TE（A.26）式中：TP 城市总体销售房屋自然景观溢价（元/a）；PPi 第 i 个样本房屋的单位面积房屋的自然景观溢价（元/a），核算见附录 D；Pi 第 i 个样本房屋的单位面积成交价（元/a）；TE 城市总体销售房屋金额（元/a）。A.3.2.2 酒店租赁自然景观溢价服务 总的酒店租赁自然景观溢价核算按式（A.27）。=(PVi PiPViM PViM PVi+N Pi THi)（A.27）式中：THP 酒店租赁总的自然景观溢价（元/a）；PVi 第 i 片区的景观房平均价格（元/m2）；DB4403/T 1412021 15 Pi 第 i 片区的非景观房平均价格（元/m2）；M,N 调查获得景观房与非景观房按照 M：N 比例销售；THi 第 i 片区的总体酒店营业收入（元/a）。A.3.3 康养服务 A.3.3.1 减少呼吸道疾病服务 以 GB 30952012 中二级标准作为清新空气的服务基准，以 PM10、PM2.5和 O3作为代表污染物，通过健康效应计算暴露人口的变化量，进而计算出因空气质量改善带来的发病损失和致死损失的变化量。核算按式（A.28）。N=Nop+Nh+Nd（A.28）式中：N 因大气环境改善而带来的呼吸系统疾病患者人数总变化（人/a）；Nop 因大气环境改善而带来的呼吸系统疾病患者门诊人数变化（人/a），核算见附录 D；Nh 因大气环境改善而带来的呼吸系统疾病患者住院人数变化（人/a），核算见附录 D；Nd 因 PM10、PM2.5、O3浓度变化而带来的死亡人数变化（人/a），核算见附录 D。DB4403/T 1412021 16 附 录 B（规范性）生态系统生产总值价值量核算方法 生态系统生产总值功能量与其对应单价的乘积即为价值量。其中，部分类型生态系统产品功能量即为产值数据（或是直接可从统计年鉴查询到产值数据）。生态系统生产总值价值量核算按式（B.1）、式（B.2）和式（B.3）。=1 （B.1）=1 （B.2）=1 （B.3）式中：EPV 生态系统物质产品价值（元/a）；EPi 第 i 类物质产品的功能量（见附录 A）；Pi 第 i 类物质产品的价格（见附录 D）；ESV 生态系统调节服务价值（元/a）；ESj 第 j 类生态系统调节服务功能量（见附录 A）；Pj 第 j 类生态系统调节服务功能的价格（见附录 D）；ECV 生态系统文化旅游服务价值（元/a）；ECk 第 k 类生态系统文化旅游服务功能量（见附录 A）；Pk 第 k 类生态系统文化旅游服务的价格（见附录 D）。注 1：农林牧渔产品和生态能源产品：农牧渔业及生态能源产值的计算均直接采用深圳统计年鉴或深圳市农业年报中对应的产品产值数据表示。注 2：旅游休闲服务：运用旅行消费法，通过有效的调查问卷来调查统计消费者对于风景观赏、休闲游乐的实际支付（或支付意愿）计算旅游服务价值。注 3：房地产景观溢价：利用链接地产的房屋交易数据计算，回归方程因变量为房屋价格（单位为人民币），自变量为一定缓冲区内的自然景观结构与社会经济结构；通过模拟有无自然景观获得价格差即为房地产景观溢价价值，因此无需额外的价值量转换。注 4：酒店景观溢价：利用携程酒店抽样调查数据计算，调查每个样本酒店的景观房与非景观房价格差（单位为人民币），然后利用当年统计年鉴酒店收入等比例计算得到，因此无需额外价值转换。注 5：削减洪涝服务：发生在人口集聚区（不适宜建设水库地方，划定方法参考附录 D）的植被削减洪涝服务，其替代工程应采用城市海绵工程；发生在其他地方的植被削减洪涝服务，其替代工程应采用水库工程。DB4403/T 1412021 17 附 录 C（规范性）生态系统生产总值核算方法 所有类型的生态系统产品和服务价值之和，即为生态系统生产总值，核算按式（C.1）。GEP=EPV+ESV+ECV（C.1）式中：GEP 生态系统生产总值（元/a）；EPV 生态系统物质产品价值（元/a）；ESV 生态系统调节服务价值（元/a）；ECV 生态系统文化旅游服务价值（元/a）。DB4403/T 1412021 18 附 录 D（资料性）生态系统生产总值核算参数补充 D.1 减少泥沙淤积服务核算参数 D.1.1 降雨侵蚀力因子 降雨侵蚀力因子R核算方法有两种方法，可根据数据可得性选择：基于月降雨量的计算模型，核算按式（D.1）；=1.735 10 1.5lgpi2p 0.818812i=1（D.1）式中：R 降雨侵蚀力值（MJmmhm-2h-1a-1）；i 月份；Pi 月降雨量（mm）；P 年降雨量（mm）。基于日降雨量资料的半月降雨侵蚀力模型来估算降雨侵蚀力，核算按式（D.2）、式（D.3）、式（D.4）和式（D.5）。=Ri24n=1（D.2）=1（D.3）=0.8363+18.14412+24.45512（D.4）=21.586 7.1891（D.5）式中：R 降雨侵蚀力值（MJmmhm-2h-1a-1）；Ri 某半月时段的降雨侵蚀力值（MJmmhm-2h-1a-1）；、模型待定参数；Dj 半月时段内第 j 天的侵蚀性日雨量（mm）（要求日雨量大于等于12mm，否则以0计算）；k 半月时段内的天数，半月时段的划分以每月第 15 日为界，每月前 15 天作为一个半 月时段，该月剩下部分作为另一个半月时段，将全年依次划分为 24 个时段；Pd12 日雨量 12mm 以上（包括等于 12mm）的日平均雨量（mm）；DB4403/T 1412021 19 Py12 日雨量 12mm 以上（包括 12mm）的年平均雨量（mm）。D.1.2 土壤可蚀性因子 土壤可蚀性因子K核算按式（D.6）。K=2.1 10-4 12 OM M1.14+3.25 S 2 +2.5 3/100 0.1317（D.6）式中：K 土壤可蚀性值（thm2hhm-2MJ-1mm-1）；OM 土壤有机质含量百分比（%）；M 土壤颗粒级配参数，即粒径分级制中（粉粒+极细砂）与（100-粘粒）百分比之积；S 土壤结构系数，取值见表 D.1；P 渗透等级，取值见表 D.1。注：粘粒为（0.002 mm）；粉粒为（0.002-0.05 mm）；极细砂为（0.05-0.1 mm）；砂粒为（0.1-2.0 mm）。表 D.1 结构系数与渗透等级的定义 结构性指数结构性指数 S S 含义含义 可渗透性指数可渗透性指数 P P 含义含义 1 非常坚固 1 快速 2 很坚固 2 中快速 3 较坚固 3 中速 4 坚固 4 中慢速 5 慢速 6 极慢 D.1.3 坡长和坡度因子 坡长因子L的核算按式（D.7）、式（D.8）和式（D.9），坡度因子S的核算按式（D.10）。L=（22.13）m（D.7）=/1+（D.8）=sin0.0896/3.0 sin0.8+0.56（D.9）S=10.8sin()+0.03,18%（D.10）式中：L 坡长因子；DB4403/T 1412021 20 坡长（m）；m 无量纲常数，取决于坡度百分比值；S 坡度因子；坡度（弧度）。D.1.4 植被覆盖因子 植被覆盖因子C采用植被覆盖度计算得到，核算按式（D.11）。植被覆盖度 f 基于植被指数 NDVI数据计算得到，核算按式（D.12）。C=1 =00.6508 0.3436lgf 0 78.3%（D.11）=(NDVI NDVIsoil)(NDVImax NDVIsoil)（D.12）式中：C 植被覆盖因子；f 植被覆盖度（%）；NDVIsoil 纯裸土象元的 NDVI 值；NDVImax 纯植被象元的 NDVI 值。D.1.5 水土保持措施因子 在P不能通过实验测定的情况下，可以通过土地利用类型的典型P值替代，不同土地利用类型取值见表 D.2。表 D.2 不同土地利用类型的水土保持措施因子P值 土地利用类型土地利用类型 耕地 园地 林地 草地 裸地 建设用地 水体 P P值值 0.5 0.6 1 1 1 0 0 注 1：水土保持措施因子反映水土保持措施对土壤侵蚀的抑制作用，是指采取水土保持措施的土壤侵蚀量与采用顺坡种植时土壤侵蚀量的比值。注 2：P值的取值范围为 0-1，值越小，表示水土保持能力越好，土壤越不易被侵蚀。D.2 调节气候服务核算参数 D.2.1 不同生态系统类型单位面积蒸腾参数 可由城市生态监测数据获取或文献获取。D.2.2 水面蒸发量参数 DB4403/T 1412021 21 可由城市气象监测数据获取或文献获取。D.3 固定二氧化碳服务核算参数 D.3.1 净生态系统生产力法核算参数 净生态系统生产力 NEP 核算方法有以下两种：如果异氧呼吸消耗数据可得，NEP 由 NPP 减去异氧呼吸消耗得到，核算按式（D.13）；NEP=NPP RS（D.13）式中：NEP 净生态系统生产力（tC/a）；NPP 净初级生产力（tC/a）；RS 异养呼吸释放碳量（tC/a）。如果异氧呼吸消耗数据不可得，NEP 根据本地NEP 和 NPP 的转换系数计算得到，核算按式（D.14）。NEP=NPP MC6/MC6H10O5（D.14）式中：NEP 净生态系统生产力（tC/a）；NEP 和 NPP 的转换系数；NPP 净初级生产力（t干物质/a）；MC6/MC6H10O5 干物质转化为 C 的系数，72/162。D.3.2 固碳速率法核算参数 固碳速率法中森林（及灌丛）固碳量FCS核算按式（D.15），草地固碳量GSCS核算按式（D.16），湿地固碳量WCS核算按式（D.17），农田固碳量CSCS核算按式（D.18）。FCS=FCSR SF （1+）（D.15）GSCS=GSR SG（D.16）=SCSR SWi 10-2ni=1（D.17）CSCS=(BSS+SCSRN+PR SCSRS)SC（D.18）式中：FCS 森林（及灌丛）固碳量（tC/a）；FCSR 森林及灌丛的固碳速率（tCha-1a-1）；DB4403/T 1412021 22 SF 森林及灌丛面积（ha）；森林及灌丛土壤固碳系数；GSCS 草地固碳量（tC/a）；GSR 草地土壤的固碳速率（tCha-1a-1）；SG 草地面积（ha）；WCS 湿地固碳量（tC/a）；n 水域湿地的种类；i 水域湿地类别，i=1，2，n；SCSRi 第 i