随着九洲江流域生猪等养殖业的迅猛发展, 畜禽养殖污染日趋严重, 对当地流域水环境造成巨大压力, 畜禽养殖业成为流域污染的主要来源之一。流域水资源生态补偿是平衡上下游发展与保护关系、解决跨区域污染的有效途径, 亦是环境经济学、环境管理学领域研究的热点问题。自20世纪90年代末生态补偿理论被引入到流域管理领域, 国内外专家、学者对流域生态补偿进行了积极的研究探索^[1-3]。目前, 国内外学者针对流域生态补偿的研究主要集中在生态补偿理论内涵^[4-6]、补偿模式^[7-9]、补偿标准^[10-13]、支付意愿^[14-19]等方面。流域生态补偿机制理论的研究主要从跨省流域生态补偿、省内跨市流域生态补偿、市域内跨县流域生态补偿等3个层次展开, 剖析涵义、原则、标准、驱动力与障碍、政策服务等^[20]。在补偿模式与途径设计方面, 主要提出政府主导型和市场交易型模式, 且提出了各自应用的条件与准则, 剖析了现有模式的局限问题, 建立了生态补偿制度选择和操作范式^[21]。现阶段, 流域生态补偿标准的测算主要从投入成本、环境效益、补偿意愿等角度切入^[22], 主要方法包括条件价值法、机会成本法、价值和费用分析法等^[23], 这些核算方法为流域生态补偿标准的确定提供依据。而有关养殖污染主控流域管理方面研究主要集中在养殖污染治理成本效益、养殖户偏好、管理模式、养殖模式、受偿意愿等^[24-26]。有些学者探讨了畜禽养殖治理和管理模式优化相结合的定量研究^[27]；孙若梅^[28]探讨了畜禽养殖业生态补偿内涵的界定, 并提出具体的生态补偿措施；张郁等^[29]分析了生态补偿情景下家庭资源禀赋对养猪户的环境行为影响；宾幕容等^[24]以湘江流域367位生猪养殖户的实地调查数据为基础, 对农户生猪养殖污染治理意愿的影响因素及层次结构进行了分析；昌苗苗和陈洋^[30]对九洲江流域生态补偿实践进展进行了介绍；王西琴等^[31]以九洲江为研究区域, 计算了限制生猪养殖的机会成本与环境效益, 建立了基于畜禽养殖模式转变的生态补偿空间优化模型；唐潥等^[32]以九洲江流域生猪养殖户为研究对象, 采用意愿调查法分别给出禁养与升级改造受偿意愿。这些理论与实践成果为本研究的开展奠定了很好的方法学基础, 但在养殖污染主控流域水资源生态补偿量的分析框架与技术方法上仍有较大的发展空间。本文以国家生态补偿示范区九洲江流域作为典型案例, 基于Ⅲ类水质目标要求计算污染物削减目标量, 设置禁养区关闭拆迁养猪场、高架床养殖模式改造、养殖粪尿废水集中处理三种典型控制方式, 以污染物削减目标为约束条件, 基于费效比构建线性规划模型, 求算用于水污染控制的最小政府补偿总额, 并以取水量、用水量、GDP为主要约束条件, 建立上下游地方政府及中央政府生态补偿金额分配模型。以期为流域生态补偿政策的制定提供参考和借鉴, 对其他流域生态补偿量的确定提供方法学参考。

1 研究区域概况

九洲江是一条跨粤桂两省(区)的独流入海河流, 发源于广西陆川县沙坡镇秦镜村的文龙径, 流经广西玉林市陆川、博白两县共10个镇进入广东湛江市的鹤地水库, 最后注入北部湾, 全长168 km, 总集雨面积3396 km²。鹤地水库位于湛江廉江市河唇镇, 是九洲江流域下游湛江市重要的饮用水源地, 水面面积共122 km²。近年来, 以生猪养殖主控的畜禽养殖业是流域上游广西玉林市的主要发展产业, 养殖源成为九洲江流域水体污染的主要来源, 养殖污染排放中化学需氧量、氨氮、总氮、总磷占总量比例分别为60.69%、57.83%、66.00%、78.51%^[25]。粤桂两省(区)党委、政府和社会各界高度重视九洲江流域水环境安全, 2015年3月, 粤桂两省联合出台《九洲江流域水环境补偿实施方案》；同年9月, 中共中央、国务院明确将九洲江流域环境综合治理列为国家《生态文明体制改革总体方案》跨地区生态补偿三个试点之一。2016年3月, 两省签订《九洲江流域上下游横向生态补偿协议》, 在两省及中央的共同努力下, 流域水环境已经得到初步改善。基于流域各乡镇生猪养殖现状, 欲通过3个方面进行流域畜禽污染治理工作(图 1)：一是根据《粤桂两省九洲江水污染防治规划》, 划定禁养区(干流沿岸500 m及支流沿岸200 m范围内)、限养区(干流500—2000 m及支流沿岸200—2000 m范围内)及适养区(其他区域), 实施生猪养殖精细化分区管理；二是推广生态养殖新模式, 实现高架床养殖；三是进行养殖废弃物处理设施建设, 实现废弃物集中处理。因此, 如何用最小的生态补偿稳定达到水质要求, 取得环境效益目标的同时兼顾九洲江上游广西玉林片区为保护水质所作的经济牺牲, 是九洲江流域畜禽养殖污染治理和生态补偿政策制定面临的重要问题。

2 流域生态补偿量计算方法

流域生态补偿额度的研究, 从最初定性确定补偿理论机制, 到定量研究具体补偿量, 方法逐渐明确科学, 也更加具有合理性。当前, 投入、效益、补偿意愿等3方面是绝大多数流域生态补偿额度测算的主题方向。在投入方面, 通过计算为保护或改善水资源环境状况, 所进行的各项投入, 包括保护环境而丧失的机会成本。主要方法包括机会成本法^[33]、生态重建成本分摊法^[34]、损失价值核算法^[35]等；在效益方面, 通过估算保护环境投入在经济、社会、生态等方面产生的外部效应来确定补偿值, 主要有生态系统服务价值法、水质水量法、污染权等^[36]方法；在补偿意愿方面, 通过对消费者进行直接调查和询问, 了解消费者的对于改善或保护环境的支付意愿, 主要有条件价值法、博弈法等^[37-38]。

补偿资金额度的测算是实现生态补偿的前提, 也是生态补偿的关键环节。基于现有研究成果, 本研究基于水质水量的跨界区域生态补偿量计算方法, 结合九洲江流域的主要特点, 建立上下游地区及中央生态补偿标准, 推动流域生态补偿具体实施与贯彻执行。

2.1 水环境目标约束下的流域生态补偿量 2.1.1 流域水环境污染物削减目标总量计算

水环境容量是在特定水环境目标约束下, 综合考虑水体自净能力下能够容纳的最大排污量^[34-35]。模型选取COD、总磷、总氮为指标, 在Ⅲ类水质要求下, 基于水环境容量为约束条件计算畜禽污染物削减目标总量。以下为削减目标总量模型函数：

(1)

(2)

(3)

式中, Q_削为污染物削减目标量；Q_入河为畜禽养殖污染物现状排放入河量；W为水环境容量；C_i为畜禽污水污染物浓度(见表 1), mg/L；T₁、T₂为禁养区、限养区生猪存栏量, 头；q₁为传统养殖模式日排水量, 0.03 t头^-1 d^-1(实地养殖场调研求均值所得, 方差δ²为0.001)；r₁、r₂为禁养区、限养区污染物入河系数(根据文献查阅^[39-41], 流域内禁养区范围污染物入河系数在0.60—0.70, 限养区范围内污染物入河系数在0.30—0.60。基于九洲江流域实地调查(图 2), 禁养区范围内养殖密度大, 取最大入河系数0.70, 限养区范围内取均值0.50)；Q₀为流域多年平均径流量, 9.28亿m³/a；C_s为污染物的环境标准值(Ⅲ类)；C₀为污染物的环境背景值, 这里忽略本底值的影响, 假设C₀ = 0；V为流域水资源体积, 11.05亿m³；K为降解系数, 根据文献现有研究^[42-43], 确定COD降解系数为0.20 d^-1；TP、TN降解系数为0.10 d^-1。

2.1.2 污染物削减分配优化模型构建

在污染物削减目标确定的基础上, 总削减量的分配采取整型规划及优先次序递推法, 禁养区内禁养清拆为第一步；限养区处理方式为第二步, 并设高架床养殖模式改造和污染物集中处理两种情景。按照先后顺序对禁养区、限养区进行可削减量估算, 比较总目标削减量和各部分可削减量, 若禁养区可满足总削减量要求, 则完成削减目标；若禁养区不能满足总削减量要求, 则递推至限养区, 直到累积可削减量达到总目标削减量值为止^[44](图 3)。

2.1.3 生态补偿量计算

限养区内有两种畜禽养殖污染治理方案, 不同治理方案的治理成本存在差异, 因此本文基于限养区不同污染治理方案设置不同成本计算情景, 从而选取最小生态补偿量。根据玉林市环保局出台的《粤桂两省九洲江水污染防治规划》, 其限养区污染治理分高架床养殖模式改造、养殖粪污集中式处理点建设两种方案, 因此生态补偿成本计算可设置两种情景, 即“禁养+高架床生态养殖改造”为情景一, “禁养+集中式处理设施建设”为情景二, 以最小生态补偿金额为目标函数, 以污染物削减目标量为主要约束条件, 得到目标函数如下：

(4)

(5)

(6)

(7)

(8)

式中, Y为最小目标生态补偿量, 元；Y₁、Y₂为不同情景下的生态补偿成本, 元；Q_禁为禁养区内可削减污染物量, t；Q_高架床为高架床养殖模式改造可削减污染物量, t头^-1 a^-1; M为高架床养殖模式猪舍容积比, 1.10 m²/头(两次实地调研取均值所得, 方差δ²为0.11)；q₂为高架床养殖模式下生猪养殖排水量, 0.0018 t头^-1 d^-1(两次调研取均值所得, 方差δ²为0.0002)；F₁、F₂分别为禁养区内清拆禁养、限养区内高架床养殖改造治理策略下的补偿金额, 分别为500.00、440.00元/ m², 数据参考于当地政府治理规划文本报告^[45], 且实地调查进行验证；F₃为集中式养殖粪污处理模式下补偿金额, 100.00元/t, 数据来自博世科陆川县固废新能源项目实地考察, 此项目是流域上游地区唯一的养殖粪污资源化处理的项目点, 数据参考于当地政府对此项目的补偿标准。

2.2 九洲江流域生态补偿量分担模型构建

围绕九洲江水环境污染的现状和特点, 粤桂双方多次协商。2014年8月, 广东省人民政府和广西壮族自治区人民政府联合签署《九洲江流域跨界水环境保护合作协议》；2016年3月21日, 在环保部、财政部的指导下, 广东、广西两省区在福建省签订了《九洲江流域上下游横向生态补偿协议》。两省虽有较好的合作基础, 然而九洲江流域上下游及中央补偿比例还是存在纷争。本研究结合流域上下游用水量以及经济承受能力, 引入中央政府财政纵向转移支付, 建立地区及中央生态补偿成本分摊模型, 以期为九洲江养殖水污染生态补偿的实施提供科学建议。具体模型如下：

(9)

(10)

(11)

式中, Y为流域生态补偿量, 元；Y_中央为中央政府纵向财政转移支付, 元；Y_广西为广西省生态补偿成本分摊量, 元；Y_广东为广东省生态补偿成本分摊量, 元；Q_广西为上游广西取水量, m³；Q_供水量为九洲江流域供水量, m³；G_上游为广西玉林市GDP总值(2016年), 亿元；G_上下游为上下游(广西玉林、广东湛江)GDP总值, 亿元。

2.3 参数灵敏度计算

考虑到模型参数的取值容易受到外界环境的影响, 且具有不确定性。本文利用Crystal Ball软件进行模型参数的灵敏度分析。首先采用蒙特卡罗模拟方法对模型参数的概率分布函数进行大量的重复抽样, 代入式(2)、式(3)、式(5)、式(6)建立生态补偿量的分布函数, 从而获得模型中重要的数学特征, 如数学期望、方差、区间估计等。在此步骤基础上, 利用方差分析法对参数进行方差分析, 根据参数输入对生态补偿量输出的秩相关系数的平方归一化后所占的比例来确定参数的灵敏度^[46]。具体公式如下所示：

(12)

(13)

式中, γ_k为k参数变量的方差贡献率；α_k0表示k参数变量对生态补偿量的秩相关系数；n表示样本数；R_ik表示参数变量k第i次抽样在样本中的排序；R_i0表示第i次抽样对应的生态补偿量的结果排序。

3 结果分析 3.1 污染物削减目标总量

根据公式(3)计算可知, COD、TP、TN水环境容量分别为20682.00、186.93、934.65 t/a。2016年九洲江流域生猪养殖头数为311.00万头, 其中禁养区范围内有66.90万头, 限养区生猪存栏量为244.10万头。在Ⅲ类水质要求下, 流域COD、TN、TP污染物削减目标总量分别为378752.98、13951.70 t和2161.60 t。

3.2 生态补偿量及分摊比例

根据公式(5)、(6)计算可知, 其情景一(“禁养+高架床生态养殖改造”模式)测算出生态补偿金额为15.03亿元；情景二(“禁养+集中式污水处理设施建设”模式)测算出生态补偿金额为15.55亿元。按照成本最小化补偿原则确定流域生态补偿量为15.03亿元。九洲江流域多年平均总供水量为13.97亿m³, 其向广西、广东境内供水量分别为3.00亿m³和10.97亿m³。据2016年统计公报显示, 广西玉林市GDP为1553.91亿元, 广东湛江市GDP为2584.78亿元。经计算, 中央、广西、广东分担率分别为0.25、0.22和0.53, 中央、广西、广东生态补偿分担金额分别为3.76、3.31、7.96亿元。

3.3 参数灵敏度分析

本文采用Crystal Ball软件进行蒙特卡罗模拟和灵敏度分析, 模型的输入设置了高架床养殖容积比(k1)、限养区污水入河系数(k2)、高架床养殖模式日排水量(k3)及传统养殖模式日排水量(k4)等4个参数变量。预测变量为生态补偿量。采用蒙特卡罗法进行5000次重复抽样, 得到生态补偿量预测值得直方图。95%置信区间水平下, 情景一生态补偿量置信区间为12.51—17.42亿元, 平均值为14.95亿元(图 4)；情景二生态补偿量置信区间为9.34—21.93亿元, 平均值为15.53亿元(图 4)。生态补偿量预测值的置信区间包含了上述计算结果, 代表了4个模型参数不确定条件下的所有可能结果。情景一下的生态补偿量大于15.00亿元的确定性内概率为45.95%；情景二下的生态补偿量大于15.00亿元的确定性内概率为52.66%(图 5)。说明情景一下的流域水污染治理在大部分情况下较情景二的方案更为经济, 与本文上述生态补偿量具体化计算结果一致。

情景一下引起生态补偿量变化的灵敏参数为高架床养殖容积比(k1)、限养区污水入河系数(k2)、高架床养殖模式日排水量(k3), 方差贡献率分别为79.50%、18.90%和1.60%；情景二下引起生态补偿量变化的灵敏参数为限养区污水入河系数(k2)、传统养殖模式日排水量(k4), 方差贡献率依次为88.80%、11.20%(图 6)。总体来看, 情景一下的高架床养殖容积比(k1)参数对流域生态补偿量最敏感, 其余参数影响较少。情景二模型中涉及两个参数, 其限养区污水入河系数(k2)为主要影响生态补偿量的参数因子。

综上所述, 不同流域水污染治理方案的选择下, 其需要的生态补偿成本存在差异, 影响其生态补偿成本的参数贡献率也各有不同。本文通过模拟计算不确定条件下的生态补偿量的所有可能值以及进行模型参数灵敏度分析, 为决策者提供了更全面的信息。

4 结论与展望 4.1 结论

流域生态补偿是一个跨行政区域的综合生态环境问题。生态补偿量的计算和测定是流域生态补偿的前提, 也是决定能否顺利实施补偿的关键环节。针对九洲江流域水生态补偿政策的问题与不足, 本文尝试构建以水质达标为前提的、以最小投入为目标的、兼顾责任与公平的跨区域流域上下游政府水环境资源生态补偿量计算模型, 主要结论如下：

在Ⅲ类水质目标要求下, 计算得出水环境治理生态补偿资金最小投入为15.03亿元。其禁养区内生猪养殖拆迁补偿为3.35亿元, 限养区内高架床生态养殖模式改造补偿为11.68亿元。中央及广西广东地方政府生态补偿分摊金额分别为3.76、3.30亿元及7.96亿元。此外, 模型参数的蒙特卡罗模拟和灵敏度分析结果显示, 95%置信水平下, 情景一生态补偿量置信区间为12.51—17.42亿元, 平均值为14.95亿元；情景二生态补偿量置信区间为9.34—21.93亿元, 平均值为15.53亿元。情景一模型中高架床养殖容积比(k1)、限养区污水入河系数(k2)、高架床养殖模式日排水量(k3)等参数对生态补偿量具有一定的影响, 方差贡献比分别为79.50%、18.90%和1.60%；情景二模型中限养区污水入河系数(k2)参数对生态补偿量最敏感, 方差贡献率达到88.80%, 其传统养殖模式日排水量(k4)参数对生态补偿量影响较小, 方差贡献率为11.20%。可见灵敏度分析结果代表了模型中参数不确定条件下生态补偿量的所有可能值, 为决策者提供了更全面的信息。

4.2 不足与展望

通过上文研究分析可以看出, 该文以保持水质投入的成本为补偿金额, 构建了基于水质水量的生态补偿量测算模型, 使得补偿量得到了合理量化, 但如何形成一整套完备的评价理论、指标体系、实施原则与计算方法等许多科学问题依然存在模糊和不确定性, 生态补偿成本构成方面仍然有待更深入的探讨和分析, 以避免统计的不完全。此外, 养殖规模, 养殖空间差异等也是生态补偿研究中重要内容, 生态补偿标准在不同地区, 不同养殖规模与方式存在差异性, 与地区养殖特点有一定的相关性, 在今后补偿标准制定中还应考虑到此异质性。

对于未来还可开展的研究工作, 一方面, 本文未考虑生态补偿方式的可行性, 本文生态补偿模式是资金的纵向转移和横向支付, 由于存在诸多的制度限制与障碍, 现实中单一的资金补偿模式无法满足实际需要, 也不容易操作, 因此, 还可考虑到补偿方式的不同, 采取多种方式优势互补, 将产业扶持、生态旅游、绿色农业等新式产业纳入到生态补偿机制中。另一方面, 畜禽养殖过程中会使用到一定比例的抗生素药剂, 畜禽粪便污水中会携带一部分抗生素残留, 未来的生态补偿研究还应适当的考虑到抗生素等难降解污染物的控制。因此该研究的生态补偿模型根据实际情况还需要进一步的修正。建立适合我国流域的生态补偿还需各利益相关方的共同参与, 各种制度的完善, 还有待进一步的研究。