污染物超排对广州土地资源规划的影响

林宜涛，李永平，王 佳

(华北电力大学 环境科学与工程学院,北京 102206)

摘 要： 以广州市城市生态系统为例，构建了基于土地交易机制的区间两阶段随机规划模型。模型考虑了工业排放的四种污染物和生态用地产生的八种生态服务价值，并分析它们对城市土地资源规划管理的影响。结果表明，随着污染物排放水平的增加，其土地交易量也随之增加。同时，土地交易机制也可以有效地控制污染物的排放量，以保证生态系统的可持续发展。最后，根据研究结果建议决策者对汽车制造业与电子制造业进行必要的产业升级，使其向着高技术、高土地利用效率、低污染物排放的方向进行产业转型。

关键词： 不确定性；区域生态系统；土地交易机制；区间两阶段随机规划；工业污染

0 引言

生态系统对人类至关重要，人类生存所必须的生活资料、经济生产资料、安全和公共健康环境都是由自然生态系统提供^[1]。同时，自然生态系统为人类和其他生物提供了各种直接和间接的服务和无形的利益^[2]。这些生态系统服务及其经济价值在过去几十年中成为科学家、决策者和利益相关者关注的焦点^[3-6]。例如，于开芹^[7]定量分析了城市化进程中生态用地变化对生态服务价值变化的影响；张天海^[8]等人研究了滨海城市的土地利用格局演变对生态系统服务价值的影响；吴后建^[9]等人对湿地生态系统服务价值进行定量评价，为区域湿地保护管理提供重要的决策依据。然而，诸如城市化进程加快、工业发展、人口增长和人口活动扩张等一系列问题导致生态系统和服务价值的巨大损失和恶化。

具体而言，人类活动加剧了自然资源与有限土地储备之间的冲突，并改变了土地利用模式。对于城市生态系统而言，为了最大限度地提高经济效益，高经济效益的用地类型就会对低经济效益的用地类型进行侵占，其中最主要的就是工业用地对生态用地的侵占。而工业用地的扩张将直接导致污染物的超标排放，其中主要的超标排放污染物为固废、SO₂、PM₁₀和工业废水。超标的污染物会对生态系统产生影响，同时工业用地对生态用地的直接侵占也进一步限制了生态系统的自我调节能力，最直接的结果就是导致生态系统的破坏。针对这些问题，陈兴旺^[10]基于城市复合生态系统理念，建立了城市规划方案的评价体系，通过建立模糊综合评价模型，提出了城市规划的发展之路；陈天鹏^[11]基于生态城市理念，用理论和实证相结合的方法，对生态城市的生态规划技术方法进行了研究，探讨了生态城市的产业生态化问题；石永林^[12]基于可持续发展的生态城市的理论，提出了基于可持续发展的生态城市的概念，并建立了相应的评价指标体系。但在区域经济发展和生态环境保护的过程中伴随着诸多不确定性，主要来自于它与随机的自然资源的管理过程、参数的估计误差、数据信息的不足等。而这些系统中存在的不确定性也影响着城市生态系统规划的实施。

1）需求表达和汲取机制。面对日益增长的公共服务需求，艾医卫等［50］提出建立需求表达机制和需求汲取机制。针对农民的利益表达机制缺失，胡志平［12］提出重构农民的利益表达机制与公共参与机制。王蔚等［51］认为，建立农民对公共服务需求的表达机制，应注重农民作为农村公共服务需求主体的作用，建立起以尊重农民个体权利基础上的多样的组织化公共服务需求表达机制。

因此，本研究结合区间参数规划(IPP)和两阶段随机规划(TSP)方法，采用区间两阶段随机规划方法^[13](ITSP)处理生态系统与城市系统中存在的不确定性；同时，为了解决工业用地对生态用地的侵占问题以及其引发的污染物超标排放问题，引入土地交易机制，以寻求最优的土地利用规划方案。最后通过建立区域生态系统的区间两阶段随机规划模型对广州市进行土地利用规划。

1 基于区域生态系统的区间两阶段随机规划模型

本研究目的在于将有限的土地资源进行合理分配，使系统收益达到最大的同时保证生态系统中各种污染物的浓度控制在国家所规定的标准范围内。本研究采用区间两阶段随机规划模型，第一阶段在保证区域土地资源承载力的前提下，确定各工业用地与生态用地的土地资源量的目标值。但是，由于区域城市生态系统中存在诸多不确定性，并且会对第一阶段的决策带来影响。例如，在发生工业污染超排事件时，超排的这部分污染物必然会对区域生态系统产生破坏。因此，在模型中引入惩罚系数(如表1)对第一阶段超排的污染物进行处罚，以补偿受到破坏的生态系统。而超排事件一旦发生则需要进行第二阶段的规划，以此调整土地规划决策，保证区域生态系统稳定。因研究区域的主要污染物为工业废水、固废、SO₂和PM₁₀，所以本研究只考虑以上四种污染物对生态系统的影响。具体的规划模型如下：

目标函数：

认知灵活性指顺应改变的情境而转换到另一种思维或行为，以符合新情境的需要的能力反映思维和行为的适应能力。Ajilore等[24]研究发现，双相障碍Ⅰ型稳定期患者(n=22)认知灵活性与健康对照人群(n=20)比较差异无统计学意义。刘传朋[25]研究显示青少年双相抑郁患者无论是在急性期还是缓解期均存在认知灵活性等方面的缺陷。

Maxf =PRO_IND^±+PRO_ECO^±-

(1)

式中，f 表示系统总收益；PRO_IND表示工业用地净收益(元)；PRO_ECO表示生态用地净收益(元)；p 表示随机事件发生的概率；i 表示不同工业类型，i =1表示汽车制造业，i =2表示电子制造业，i =3表示石油化工业；t 表示不同规划期，t =1，2，3分别表示第一，第二，第三规划期；e 表示随机事件发生的不同概率水平，e =1，2，3分别表示低，中，高三种概率水平；y 表示第二阶段土地交易量(ha)；L 表示规划期年限(年)；PW_it 表示i 工业在t 时期单位产品生产所排放的污水量(m³/ton)；PC_it 表示i 工业在t 时期单位面积产品生产量(ton/ha)；PDW_it 表示i 工业在t 时期单位面积对污水的处理量(m³/ha)；PPM_it 表示i 工业在t 时期单位产品生产所排放的PM₁₀的量(kg/ton)；PDPM_it 表示i 工业在t 时期单位面积对PM₁₀的处理量(kg/ha)；PSO_it 表示i 工业在t 时期单位产品生产所排放的SO₂的量(kg/ton)；PDSO_it 表示i 工业在t 时期单位面积对SO₂的处理量(kg/ha)；PS_it 表示i 工业在t 时期单位产品生产所排放的固废的量(kg/ton)；PDS_it 表示i 工业在t 时期单位面积对固废的处理量(kg/ha)；η _it 表示固废回收率；DTWC表示废水排放超标时工业所需要承担罚款的惩罚系数；DTPMC表示PM₁₀排放超标时工业所需要承担罚款的惩罚系数；DTSOC表示SO₂排放超标时工业所需要承担罚款的惩罚系数；DTGC表示固废排放超标时工业所需要承担罚款的惩罚系数；上标+和-表示上限和下限。

总体来讲,多重随机森林加权大数投票对于动态行为分类效果较好,分析原因是使用多组最佳参数组合使波动误差缩小,同时加权投票机制使最终决策更加稳定。但是在实际的应用中,随着基分类器个数增多,模型训练及行为识别过程中时间消耗也会越多,应当根据实际场景去选择基分类器的个数。

表 1各超排污染物的惩罚系数表

注：根据何军^[14]等人对生态补偿机制的研究，可以知道在生态遭到破坏时所需要补偿的金额，而将补偿

金额与生态系统面积相比则可以获得破坏单位面积生态用地所需要赔偿的金额，即惩罚系数。

(7) 固废处理能力约束：

1.1工业净收益核算

在本研究中，工业用地(即汽车制造业，电子制造业和石油化工业)净收益主要考虑了8种收益与成本类型，主要包括：工业产品收益、工业产品成本、污水处理成本、固废回收收益、固废处理成本、SO₂处理成本、PM₁₀处理成本和氮氧化物处理成本。以下是工业用地收益的具体计算公式：

(8) 植被覆盖率约束:

(2)

式中，PCI_it 表示i 工业在t 时期单位产品的收益(元/ton)；PCR_it 表示i 工业在t 时期单位产品的成本(元/ton)；TC表示t 时期污水处理的费用(元/m³)；SC表示t 时期回收固废产生的收益(元/kg)；SF表示t 时期固废处理的费用(元/kg)；SOC_it 表示i 工业在t 时期对SO₂的处理费用(元/kg)；PMC_it 表示i 工业在t 时期对PM₁₀的处理费用(元/kg)；NOC_it 表示i 工业在t 时期对氮氧化物的处理费用(元/kg)。

2018年11月11日下午，发网第五届双11媒体见面会在发网吴江嘉合仓举行。邀请40多家媒体记者走进仓库现场，见证双11征战。与往年不同，2018年的媒体见面会设置为体验式参观，嘉合仓设置了专门的团购流水线以供媒体记者实际上线实操，40余位记者被安排在不同环节，分别完成打单、分拣、复核、出库等操作。记者所在体验的发网吴江嘉合仓总面积6万平方米，2018年双11期间预计单量将超600万，雪花秀、完美日记、奥利奥、高露洁、古今内衣、施华蔻等近60余个知名品牌产品从该仓发货。

1.2生态服务价值核算

在本研究中，生态用地(即草地、林地和湿地)产生的服务价值只要包括了8种生态服务价值类型：固碳价值、释氧价值、涵养水源价值、净化水质价值、固土价值、保肥价值、净化空气价值和旅游价值。而在生态服务价值核算方法方面，本研究主要是采用土地利用类型面积与经验系数核算法并未采用能值核算法。原因主要包括：(1)土地利用类型面积与经验系数核算法与能值核算法相比其在数据获取方面更加便捷；(2)土地利用类型面积与经验系数核算法相比于能值核算法可以核算更多的服务价值，例如旅游价值就无法利用能值核算法进行核算。以下是生态用地收益的具体计算公式：

(3)

式中，CP_kt 表示k 类型生态用地在t 时期对二氧化碳吸收的能力(ton/ha)；C_t 表示t 时期的碳税(元/ton)；O_t 表示t 时期制氧费用(元/ton)；RP_kt 表示k 类用地在t 时期单位面积涵养水源的能力(m³/ha)；RO_kt 表示k 类用地在t 时期的径流系数；RCI表示广州市水价(元/m³)；WP_kt 表示k 类生态用地在t 时期水质净化的能力(m³/ha)；ZP_kt 表示k 类生态用地在t 时期固定土壤的能力(ton/ha)；ZCI表示河道清淤的费用(元/ton)；γ_qt 表示t 时期q 类元素在土壤中的占比(%)；PI_kt 表示t 时期q 元素肥料的单价(元/ton)；I_q 表示q 元素在肥料中的占比(%)；YC表示旅游的人均花费(元/人)；YP_kt 表示k 类生态用地t 时期单位面积的旅游人数(人/ha)；WGW_krt 表示k 类生态用地t 时期对r 类大气污染物的净化能力(kg/ha)；GWC_nt 表示t 时期工业处理n 类污染物的费用(元/kg)。

1.3模型约束条件

在本研究中，模型涉及11种约束条件，主要包括：土地资源约束、废水排放约束、废水处理能力约束、PM₁₀处理能力约束、SO₂处理能力约束、固废排放约束、固废处理能力约束、植被覆盖率约束、土壤侵蚀约束、生态系统涵养水源能力约束和释氧能力约束。具体约束条件如下：

(1) 土地资源约束：

∀t ，

(4)

式中，CSW表示区域生态系统水资源总量(m³)；α_kt 表示t 时期k 类生态用地径流系数；PF_it 表示t 时期i 工业生产单位产品所消耗的水量(m³/ ton)；GW表示废水回用率(%)；RJW表示人均用水量(m³/人)；RK表示区域人口数(人)。

(2) 废水排放约束：

≤∀t ，

(5)

式中，SW_et 表示t 时期e 概率水平下废水的排放量(m³)。

式中，ZP _kt表示t时期k类生态用地的土地侵蚀系数；SEA 表示区域被侵蚀的生态用地面积(ha )。

(3) 废水处理能力约束：

≤∀t ，

(6)

式中，TPC表示污水处理厂处理的废水量(m³)。

(4) PM₁₀处理能力约束：

≤∀t ，

(7)

(12) 非负约束:

(5) SO₂处理能力约束：

≤∀t，

(8)

式中，WSO_kt 表示t 时期k 类型生态用地对SO₂的净化系数(kg/ha)；SSO_et 表示t 时期e 概率水平下SO₂的排放量(kg)。

(6) 固废排放约束：

≤∀i ,t ，

(9)

式中，SG_et 表示t 时期e 概率水平下固废的排放量(kg)。

2017年2月22日，新华网每日电讯刊载一篇名为《东北一些农业合作社宁可多花十倍钱，也要买进口犁》的报道，提出国产犁具技术落后、使用效果差等问题，导致出现了进口高价高端犁具畅销的局面，这一现象折射出我国农机产业大而不强的突出问题，虽然我国是世界第一农机制造和使用大国，但是我国农机技术性能低、产品可靠性差，这一系列问题更加表明了农机供给侧改革的必要性与紧迫性。鉴于此，笔者介绍了兵团农机化现状及发展趋势、犁具技术发展及兵团应用现状、兵团犁具生产技术现状和存在的问题以及相应的措施及建议。

分别计算N=6, 8涡流发生器流场, 其余参数采用上文的最优组合: 安装角度α=12°, 后缘高度H=1.2δ, 安装位置距试验段入口20δ. 图10和表5分别为附面层速度分布和相应的附面层厚度, 图11 为尾涡诱导速度沿流向的分布.

≤∀t ,

(10)

式中，SP_kt 表示t 时期k 类型生态用地对固废的净化系数(kg/ha)；TSP表示处理的固废量(kg)。

采用蒲公英300克、大青叶100克、制半夏100克、连翘200克、黄连200克、金银花200克、黄芩150克、桔梗100克、甘草50克，水煎原汁或共制粗粉，煎汁拌料，为5000只雏鸡1日用量，连用3～5天，粉末每只雏鸡每天0.5～0.6克，开水浸后拌料饲喂，用于防治鸡传染性支气管炎、呼吸道严重感染的治疗，有效率高达99.03%，预防量减半，连续使用3日。

随着我国矿山开采机械化水平的提高，大型矿井和千万吨级别高产高效工作面的增加，矿井建设和矿井开拓接续问题凸显[1-3]。“十二五”期间，国家“863 计划”重点项目“煤炭智能化掘采技术与装备”的“矿山竖井掘进机研制” 课题，已开发出具有安全防护的竖井综合凿井装备——MSJ5.8/1.6D型矿山竖井掘进机。MSJ5.8/1.6D型矿山竖井掘进机是一种涉及多学科的技术密集型施工设备，主要应用于竖井井筒建设，能够实现快速机械破岩、支护，是大型现代化矿井建设机械化的发展方向，具有广阔的应用前景[4-8]。

≥∀t ,

(11)

式中，PFR _it表示t时期i工业的植被覆盖率(%)；LHS 表示工业用地应满足的最小绿地面积(ha )。

(9) 土壤侵蚀约束:

≥∀t

(12)

在诺基亚产品销售渠道上，并没有采用销售多元化的方式，甚至没有很好的利用突飞猛进的互联网销售平台，这些都在不同程度上制约了诺基亚的发展。诺基亚不得已选择了与微软合作，而主要以软件发展为主，并且与众多公司都有合作，这就决定了其不能成为强有力的上游供应商。

(10) 生态系统涵养水源能力约束：

≥∀t

(13)

式中，x _it 表示t 时期i 类型工业用地所规划的面积(ha)；x _kt 表示t 时期k 类型生态用地所规划的面积(ha)；TA表示研究区域总规划面积(ha)。

(11) 释氧能力约束：

≥∀t ，

(14)

式中，ZO_kt 表示t 时期单位面积k 类生态用地释放的氧气量(ton/ha)；MINO区域生态系统维持稳定所需要的最小氧气量(ton)。

式中，WPM_kt 表示t 时期k 类型生态用地对PM₁₀的净化系数(kg/ha)；SPM_et 表示t 时期e 概率水平下PM₁₀的排放量(kg)。

≥≥≥0,∀i ，k ，e ，t 。

(15)

2 研究区域概况及数据来源

2.1研究区域概况

广州位于中国广东省中部和南部，由10个市区(天河、越秀、海珠、番禺、南沙、荔湾、白云、萝岗、黄埔和花都)和两个县级城市(从化和增城)组成。在过去的三十年里，广州经历了快速的经济增长，见证了前所未有的城市化进程，城市化率从1990年的57.4%上升到2015年的89.4%。同时，国内生产总值(GDP)从2005年的5 154亿元人民币增加到2015年的13 619亿元人民币，增长了近三倍。在所有行业中，汽车制造业、电子制造业和石油化工业是三大支柱产业，其产品占国内生产总值(GDP)的47%。但是，随着经济的快速发展和人口的增长，广州的生态系统受到严重破坏，遇到了一系列生态问题，如城市绿化面积减少、生物多样性减少、城市热岛效应增加等^[15-18]。例如，广州市20世纪90年代城区中总悬浮颗粒物平均浓度为0.256mg/m³，是推荐标准0.06～0.09 mg/m³的4.2倍，是香港的2.6倍，新加坡的4.8倍^[19]。此外，自1960年以来，广州市的林地、草地和湿地等城市绿地面积急剧下降，由此带来的生态系统的退化会阻碍城市经济的增长，同时阻碍区域生态系统的可持续发展。因此，本研究在保证区域土地资源承载力的条件下，对区域生态系统的土地利用规划提供合理建议。

2.2数据来源

本研究根据2000—2015年《广东省统计年鉴》《广州市统计年鉴》等有关资料，得到工业方面的数据，例如，产量、收益、固废回收率及废水回用率。而生态方面数据，根据徐俏^[20]等人的研究，得到广州市林地土壤肥力保持系数。根据庄大昌^[21]等人的研究，得到广州市湿地单位面积旅游收益系数为382元/hm²。表2也给出了不同时期生态影响因子的经验系数值。参考方小林、高岚^[22]的研究，可知涵养水源系数为[9149.5，9374.35] ton/ha，净化水质系数为[27.33，27.68]m³/ha。参考王璐^[23]等人的研究，获得了广州市的平均水资源总量7.98×10⁹m³，其中地表水7.88×10⁹m³；地下水1.48×10⁹m³。而根据年鉴中污染物排放量的数据，假设预测年份的3种污染物排放量水平出现的概率分别为20%、60%和20%。表3给出了不同时期不同概率水平下的污染物排放总量。

表 2模型的部分生态影响因子经验系数

注：根据方小林、高岚^[22]对广东碳汇林生态服务功能价值的研究，获得净化水质系数、涵养水源系数和固碳系数。同时，依照中华人民共和国行业标准《土壤侵蚀分类分级标准》，获得土壤侵蚀模数。

表 3不同概率水平下的污染物排放总量

注：根据《广东省统计年鉴》及《广州市统计年鉴》可以获得相应时期内废水、固废、PM₁₀和SO₂的排放总量。

3 结果分析

为探究污染物在不同消减水平下对于区域生态土地利用的影响，本研究基于污染物的消减水平引入五个情景(如表4)，同时基于两种交易模式(交易与不交易)，研究将产生10种情景，以下将对所产生的情景与结果进行分析。

3.1不同情景下区域生态系统总收益分析

表 4本研究设置的情景

图1给出了在交易模式下不同情景不同时期下的系统总收益的变化情况。从图1可以看出随着时期的增加系统的总收益也会随之增加。例如，在情景S1下，系统的总收益从第一时期的[9.84×10⁹，10.82×10⁹]元，上升到第三时期的[12.45×10⁹，14.06×10⁹]元。与此同时，随着污染物消减水平的增加系统总收益随之减少。例如，在第三时期，系统的总收益从情景S1的[12.45×10⁹，14.06×10⁹]元，下降到情景S5的[11.31×10⁹，12.85×10⁹]元。而增加污染物消减水平意味着工业需要降低污染物排放，这也就说明在更高的污染物消减水平下决策者更加倾向于生态的安全，从而对污染物采取更加严格的排放政策，而这样的做法必将影响工业生产的效率，从而导致工业收益下降。虽然，与此同时生态服务价值所带来的收益会增加，但是与工业收益相比，生态服务价值所提供的收益较低，所以整个系统的收益呈下降趋势。图2给出了在不交易模式下不同情景不同时期下的系统总收益的变化情况。与图1相比可以看出系统收益的整体趋势是一致的，但是可以明显看出在交易模式下的系统总收益大于不交易模式下的系统总收益。例如，在第三时期情景S1时，交易模式下系统总收益为[12.45×10⁹，14.06×10⁹]元，而不交易模式下系统总收益为[9.79×10⁹，11.19×10⁹]元。这可能是因为在不进行土地交易的情况下，工业可以获得最大的经济收益，但系统需要承担超标污染物的巨额罚款；而在进行土地交易的情况下，工业由于受到土地交易的影响，其无法获得最大的经济效益，但与此同时，系统无需再承担超标污染物所带来的罚款，因此使得交易模式的系统收益要高于不交易模式的系统收益。这也就说明了通过土地交易模式可以维持生态系统的稳定性和可持续性，并且在这个基础上为系统获得最大的经济收益。

对井队主要技术人员技术方案的掌握情况进行考核；对井队在施工中技术方案的执行度进行考核；对井队主要技术负责人的业务能力进行考核；对钻井工程师、钻井液工程师负责的钻井工具、仪器、设备是否按规范使用进行考核；对施工井的井口质量、井身质量、取心收获率、施工进度进行考核。

3.2不同情景下土地交易情况分析

图3和4给出了不同情景在三种概率水平下各个用地的土地交易面积的变化情况。从图中可以看出随着污染物排放水平的提高，其土地交易面积增加，这也就意味着在高污染物排放的情况下更需要进行土地交易。例如，在情景S1下，石油化工业的土地交易面积从低水平的不进行土地交易，增加到高水平的[0.68×10³，3.62×10³]ha。这同样也说明了土地交易模式在维持生态系统可持续性方面是有效的。而随着污染物消减水平的增加，土地交易面积也随之增加。例如，在中等污染物排放水平下，汽车制造业的土地交易面积从情景S1的[0.94×10³，1.22×10³]ha，增加到情景S5的[3.53×10³，3.87×10³]ha。这可能是因为增加污染物的消减水平会增加工业排污的压力，从而导致部分工业用地与生态用地进行交易以此维持生态系统的稳定。而从图3和4中也可以看出在污染物消减处于低水平时，系统没有进行土地交易，说明在低水平时，工业排放的污染物在生态系统可接受范围内，所以这时候不需要进行土地交易。最后，根据图3和4的结果可以发现广州市现在急需进行土地交易的工业类型为汽车制造业与电子制造业，其中在中污染物排放水平时汽车制造业为主要的土地卖方，而在高污染物排放水平时电子制造业为主要的土地卖方，这也就说明汽车制造业与电子制造业应该寻求新的发展方式。

3.3土地交易机制对污染物排放的影响

图5和图6给出了不同交易模式下固废排放量的变化情况。从图中可以看出在进行土地交易模式下固废的排放量要低于不进行土地交易模式下固废排放量。例如，在第二时期地污染物排放水平下并且不进行土地交易模式时，汽车制造业和电子制造业以及石油化工业的固废排放量分别为[0.54×10⁶，1.42×10⁶]kg，[2.76×10⁶，3.22×10⁶] kg和[0.73×10⁶，1.89×10⁶] kg；而在进行土地交易模式时，汽车制造业和电子制造业以及石油化工业的固废排放量分别为[0.42×10⁶，1.14×10⁶] kg，[0.72×10⁶，1.93×10⁶] kg和[0.58×10⁶，1.03×10⁶] kg。从结果可以发现土地交易模式可以降低工业污染物的排放量，这同时也说明了土地交易模式可以有效的缓解工业对生态环境造成的破坏。

4 讨论与结论

区域生态系统是一个具有随机性、动态性和不确定性等一系列复杂特征的巨系统。而本文结合概率密度函数与区间数构建区间两阶段随机规划模型，并引入土地交易机制，对区域生态系统中的土地资源配置进行研究。在模型的第一阶段，当工业排放的污染物没有超过生态承受阈值时，工业用地与生态用地拥有其预定的目标值；而在模型的第二阶段，当工业排放的污染物超过生态承受阈值时，通过引入土地交易机制，使得超排严重的工业用地与生态用地进行土地交易，以减少工业污染排放对生态系统的危害。计算结果表明，随着污染物排放水平的增加，其土地交易量也随之增加，而在土地交易的模式下污染物的排放量要小于不进行土地交易模式下的污染物排放量。如在中污染物排放水平下，土地交易量在第1、2和3时期分别为0 ha、[3.53×10³，3.87×10³] ha和[0.56×10³，2.63×10³] ha；而在高污染物排放水平下，土地交易量在第1、2和3时期分别为0 ha、[5.54×10³，6.53×10³] ha和[2.07×10³，2.86×10³] ha。而对于污染物排放量来说，进行土地交易可以降低[2.31×10⁶，2.43×10⁶] kg固废排放量。这也说明了土地交易机制可以有效的控制污染物的排放量使其维持在生态系统可调节范围内，以保证生态系统的可持续发展。因此，建议管理决策者对汽车制造业与电子制造业进行必要的产业升级，使汽车制造业与电子制造业向着高技术，高土地利用效率，低污染物排放和更加生态的方向进行产业转型。

农副食品加工业是农业产业纵深发展的产物，对农产品进行一系列的加工活动，能够有效提高农产品的利用效率、增加农产品的附加值，从而提高农民收入、促进农业繁荣、推动农村发展。诚然，农副食品加工业带来的诸多益处是显而易见，但是，也应当注意到当前其发展仍存在一些问题，尤其是资源利用效率不高，投入产出不配比的现象。对于农副食品加工业来说，效率的高低直接影响到农民收入和农业的发展，因此，进一步探索农副食品加工业的效率变动，对于提高农业资源配置效率，解决“三农”问题具有重要现实意义。

本研究将土地交易机制运用到区域生态系统土地资源配置模型中，但是还有许多其他生态参数在区域生态系统中发挥重要作用，如水文连通性、生物连通性、生态阈值和生态幅度。通过模型难以评估这些生态参数的微小变化。因此，仍需要在未来的工作中对生态区域土地资源配置进行研究。

小学全科师范生党员学习党的时代性思想加上教师本身自主的学习教师新时代新思想，两者相互促进、相互影响，更好的践行社会主义核心价值观。

参考文献：

[1] 景永才,陈利顶,孙然好.基于生态系统服务供需的城市群生态安全格局构建框架[J].生态学报,2018,38(12):4121-4131.

[2] 谢高地,张彩霞,张昌顺,等.中国生态系统服务的价值[J].资源科学,2015,37(9):1740-1746.

[3] 邓小文,孙贻超,韩士杰.城市生态用地分类及其规划的一般原则[J].应用生态报,2006,16(10):2003-2006.

[4] 李洪远.基于遥感和GIS的天津滨海新区30年间生态用地变化分析[J].南水北调与水利科技,2013,11(1):75-80.

[5] 王振健,李如雪.城市生态用地分类、功能及其保护利用研究:以山东聊城市为例[J].水土保持研究,2006,13(6):306-308.

[6] 李 锋,叶亚平,宋博文,等.城市生态用地的空间结构及其生态系统服务动态演变:以常州市为例[J].生态学报,2011,31(19):5623-5631.

[7] 于开芹.区域生态用地变化对生态服务价值的影响[J].安徽农业科学,2018,46(23):53-60.

[8] 张天海,郭 蔚,田 野,等.滨海城市土地利用格局演变及对生态系统服务价值的影响[J].生态学报,2018,38(21):2-10.

[9] 吴后建,但新球,刘世好,等.湖南省湿地生态系统服务价值初步评价[J].湿地科学,2016,14(6):782-787.

[10] 陈兴旺.城市规划对城市复合生态系统的优化功能研究[D].西安:西北大学,2007.

[11] 陈天鹏.生态城市建设与评价研究[D].哈尔滨:哈尔滨工业大学,2008.

[12] 石永林.基于可持续发展的生态城市建设研究[D].哈尔滨:哈尔滨工业大学,2006.

[13] LI Y P，HUANG G H.Two-stage planning for sustainable water-quality management under uncertainty[J].J Environ Manage,2009(90):2402-2413.

[14] 何 军,马 娅,张昌顺,等.基于生态服务价值的广州市生态补偿研究[J].生态经济,2017,12(33):184-218.

[15] 周阳品,陈 军.广州市生态环境与经济发展阶段的实证研究[J].资源开发与市场,2009,25(8):732-735.

[16] 黄镇国,余汉豪,邹春洋,等.广州市的经济发展与环境保护[J].云南地理环境研究,2005,17(1):14-18.

[17] 陈伯通,罗建中,冯爱坤.广州地区酸雨状况及其影响因素探讨[J].环境污染与防治,2006,28(2):112-115.

[18] 曹兴兴.广州市城市生物多样性保护规划研究[D].广州:仲恺农业工程学院,2013.

[19] 张 捷,景守武.改革开放以来广州市生态文明建设经验总结[J].城市观察,2018(3):61-72.

[20] 徐 俏,何孟常,杨志峰,等.广州市生态系统服务功能价值评估[J].北京师范大学学报(自然科学版),2003,2(39):268-272.

[21] 庄大昌,杨青生.广州市城市湿地生态系统服务功能价值评估[J].热带地理,2009,5(29):407-411.

[22] 方小林,高 岚.广东碳汇林生态服务功能价值评价[J].中南林业科技大学学报(社会科学版),2016(10):1673-9272.

[23] 王 璐,杨 洁,胡月明,等.广州市土地利用生态服务价值测算研究[J].水土保持通报,2009(29):230-234.

Impact of Excessive Emission of Pollutants on Land Resource Planning in Guangzhou

LIN Yitao，LI Yongping，WANG Jia

(College of Environmental Science and Engineering，North China Electric Power University，Beijing 102206，China)

Abstract :Taking the urban ecosystem of Guangzhou as an example，an interval two-stage stochastic planning model based on land transaction mechanism is constructed.The model takes into consideration the four pollutants emitted by industry and eight ecological service values generated by ecological land，and analyzes their impact on urban land resource planning management.The results show that the land transaction volume increases as the pollutant emission level increases.Meanwhile，the land transaction mechanism can also effectively control the emission of pollutants to ensure the sustainable development of the ecosystem.In the end，based on the research results，it is recommended that policy makers should carry out the necessary industrial upgrading of automobile manufacturing industry and electronics manufacturing industry for transformation toward high technology，high land use efficiency and low pollutant emissions.

Keywords ：uncertainty;regional ecosystem;land transaction mechanism;interval two-stage stochastic programming;industrial pollution

DOI： 10.16246/ j.issn.1673- 5072.2019.02.014

收稿日期： 2018-09-17

基金项目： 国家重点研发计划重点专项(2016YFC0502803)

作者简介： 林宜涛(1994—)，男，福建三明人，硕士研究生，主要从事生态保护与恢复、城市生态系统管理与规划。

E-mail：18850543945@163.com

通信作者： 李永平(1971—)，女，湖北荆州人，博士，教授，主要从事环境系统不确定性研究。

E-mail:yongping-li33@163.com

引文格式： 林宜涛，李永平，王 佳.污染物超排对广州土地资源规划的影响[J].西华师范大学学报(自然科学版)，2019,40(2)：187-196.[LIN Yitao,LI Yongping,WANG Jia.Impact of Excessive Emission of Pollutants on Land Resource Planning in Guangzhou[J].Journal of China West Normal University (Natural Sciences),2019,40(2):187-196.]

中图分类号： X321

文献标志码： A

文章编号： 1673-5072(2019)02-0187-010

标签：不确定性论文;  区域生态系统论文;  土地交易机制论文;  区间两阶段随机规划论文;  工业污染论文;  华北电力大学环境科学与工程学院论文;