一种确定流域土地开发利用承载力的方法

本发明属于环境保护
技术领域：
，具体涉及一种确定流域土地开发利用承载力的方法。
背景技术：
：为了加强长江流域生态环境保护和修复，促进资源合理高效利用，保障生态安全，《长江经济带生态环境保护规划》明确要求“长江沿线一切经济活动都要以不破坏生态环境为前提，抓紧制定产业准入负面清单，明确空间准入和环境准入的清单式管理要求。除在建项目外，严禁在干流及主要支流岸线1公里范围内布局新建重化工园区，严控在中上游沿岸地区新建石油化工和煤化工项目。”可见，国家对长江经济带，尤其是沿江1公里范围内的产业发展做了总体要求，长江经济带多个省市也在此基础上制定了沿江1公里范围内的产业调整细则。但排入河流的污染物除了沿江范围还来自整个流域的其他点源、面源等排放源，若不将整个流域的土地开发方式、开发规模进行科学安排，将不利于河流水质保护，无法达到预期保护目标。在此背景下，深入研究流域的土地开发方式、开发规模是有必要的。但目前尚没有合理的办法，以河流环境容量为限定前提，判断某流域范围内的土地开发利用承载力是否超过上限，进而对整个流域土地开发方式、开发规模的优化调整提供决策支持。所以，需要一个新的技术方案来实现上述效果。技术实现要素：发明目的：为了克服现有技术中存在的不足，提供一种确定流域土地开发利用承载力的方法，能够定量确定某个流域的土地开发利用承载力，可为流域土地利用开发决策提供参考，达到保护河流水体的目的。技术方案：为实现上述目的，本发明提供一种确定流域土地开发利用承载力的方法，包括如下步骤：s1：调查并计算河流的水环境容量；s2：确定该河流所处的流域范围；s3：基于步骤s2确定的流域范围，调查流域范围内土地利用现状及土地利用规划，确定现状各土地利用类型的单位土地面积的水污染物排放量，即入江排污系数；s4：基于步骤s1获得的水环境容量以及步骤s3获得的土地利用规划和入江排污系数，确定流域土地开发利用承载力，判断当地土地利用规划是否超过其开发利用承载力。进一步地，所述步骤s1中水环境容量的获取方法为：给出如下公式：w＝qp(cs-c0)+kvcs式中：w为单一河段(河段为两节点之间的河道)计算单元的环境容量，通常用t/a表示；qp为计算单元的河流设计流量，m3/s；cs为计算单元污染物控制因子的环境质量目标值，mg/l。根据国家地表水污染物总量控制指标，此处控制因子可选取化学需氧量(cod)和氨氮(nh3-n)。根据所在河流各河段地表水功能区划及环境保护目标的要求，cod和nh3-n的环境质量目标值执行《地表水环境质量标准》(gb3838-2002)中相应的污染物浓度标准；c0为计算单元入流断面控制因子的浓度值，即河流中污染物的实际浓度值，可取多个河流断面监测到的污染物浓度的平均值，mg/l；k为控制因子污染物在地表水体中的自净系数，d-1；v为计算单元蓄水量，m3。进一步地，所述步骤s2中流域范围的确定方法为：基于流域或区域dem数据，通过arcgis软件的水文分析工具，提取所需流域的边界线，边界线以内的区域即为流域范围，通过arcgis软件统计工具获取流域面积。进一步地，所述步骤s3中入江排污系数的确定方法为：给出如下公式：ci＝pi/ai(i＝1,2,3…n)式中：ci为第i种土地利用类型(如工业用地、居住用地、农田用地、商业用地，等)的入江排污系数，t/a.m-2；pi为第i种土地利用类型的污染物排放量，t/a，可通过区域排污统计资料获取；ai为第i种土地利用类型的面积，m2；可通过土地利用现状统计资料获取，进一步地，所述步骤s4中判断当地土地利用规划是否超过其开发利用承载力的方法为：若有：则说明土地利用规划已超过流域开发利用承载力，当地土地开发方式和开发规模需要优化调整；若有：则说明土地利用规划未超过流域开发利用承载力。式中：gai为该流域土地利用规划中第i种土地利用类型的面积，m2；可通过土地利用规划资料获取；ci为第i种土地利用类型(如工业用地、居住用地、农田用地、商业用地，等)的入江排污系数，t/a.m-2。有益效果：本发明与现有技术相比，能够定量确定某个流域的土地开发利用承载力，可快速判断流域土地开发方式、开发规模是否超过河流水环境承载力，能够为流域土地利用优化调整提供决策支撑，免于经验决策的盲目性，从而达到保护河流水体的目的。附图说明图1为本发明方法的流程示意图。具体实施方式下面结合附图和具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。如图1所示，本发明提供一种确定流域土地开发利用承载力的方法，包括如下步骤：s1：调查并计算河流的水环境容量。水环境容量的获取方法为：给出如下公式：w＝qp(cs-c0)+kvcs式中：w为单一河段(河段为两节点之间的河道)计算单元的环境容量，通常用t/a表示；qp为计算单元的河流设计流量，m3/s；cs为计算单元污染物控制因子的环境质量目标值，mg/l。根据国家地表水污染物总量控制指标，此处控制因子可选取化学需氧量(cod)和氨氮(nh3-n)。根据所在河流各河段地表水功能区划及环境保护目标的要求，cod和nh3-n的环境质量目标值执行《地表水环境质量标准》(gb3838-2002)中相应的污染物浓度标准；c0为计算单元入流断面控制因子的浓度值，即河流中污染物的实际浓度值，可取多个河流断面监测到的污染物浓度的平均值，mg/l；k为控制因子污染物在地表水体中的自净系数，d-1；v为计算单元蓄水量，m3。s2：确定该河流所处的流域范围。具体为：基于流域或区域dem数据，通过arcgis软件的水文分析工具，提取所需流域的边界线，边界线以内的区域即为流域范围，通过arcgis软件统计工具获取流域面积。s3：基于步骤s2确定的流域范围，调查流域范围内土地利用现状及土地利用规划，确定现状各土地利用类型的单位土地面积的水污染物排放量，即入江排污系数。入江排污系数的确定方法为：给出如下公式：ci＝pi/ai(i＝1,2,3…n)式中：ci为第i种土地利用类型(如工业用地、居住用地、农田用地、商业用地，等)的入江排污系数，t/a.m-2；pi为第i种土地利用类型的污染物排放量，t/a，可通过区域排污统计资料获取；ai为第i种土地利用类型的面积，m2；可通过土地利用现状统计资料获取，s4：基于步骤s1获得的水环境容量以及步骤s3获得的土地利用规划和入江排污系数，确定流域土地开发利用承载力，判断当地土地利用规划是否超过其开发利用承载力。若有：则说明土地利用规划已超过流域开发利用承载力，当地土地开发方式和开发规模需要优化调整；若有：则说明土地利用规划未超过流域开发利用承载力。式中：gai为该流域土地利用规划中第i种土地利用类型的面积，m2；可通过土地利用规划资料获取；ci为第i种土地利用类型(如工业用地、居住用地、农田用地、商业用地，等)的入江排污系数，t/a.m-2。基于上述技术方案，本实施例中将上述技术方案应用于长江某支流流域土地开发利用承载力的判断和确定，具体的过程为：1、调查并计算所研究的长江支流河段水环境容量。选取江苏省境内某长江一级支流，该支流长度为2000m。根据当地地表水功能区划，该段水体水质标准执行《地表水环境质量标准》(gb3838—2002)中ⅳ类水标准，此处选取cod为计算演示指标。综合考虑该河流水文参数、水质目标、上游来水水质和污染物降解能力等因素，水环境容量计算模型中，各计算参数选择及环境容量w计算结果如表1。表1环境容量计算参数选取名称qp(m3/s)cs(mg/l)c0(mg/l)k(d-1)v(m3)w(t/a)cod602018.20.2651236003644.992、确定该河流的流域范围。基于流域dem数据，在arcgis10.2软件中提取了流域边界线，经统计，流域面积a为4km2。3、调查流域范围内土地利用现状及土地利用规划，确定现状各土地利用类型的单位土地面积的水污染物排放量，即入江排污系数ci。根据统计资料，现状土地利用中：工业用地面积a1＝0.5km2，c1＝1.5×10-3t/a.m-2；居住用地面积a2＝0.5km2，c2＝0.6×10-3t/a.m-2；农业用地面积a3＝3km2，c3＝0.4×10-3t/a.m-2。根据该流域土地利用规划，未来将增加工业用地和居住用地面积，减少农田用地面积，规划实施后，工业用地面积ga1＝2.5km2，居住用地面积ga2＝1km2，农业用地面积ga3＝0.5km2。4、基于获得的土地利用规划及水环境容量，判断流域土地利用规划是否超过其开发利用承载力。ga1×c1+ga2×c2+ga3×c3＝(2.5×1.5+1×0.6+0.5×0.4)×103＝4550(t/a)4550t/a>w(3644.99t/a)说明该流域土地利用规划已超过其开发利用承载力，流域土地开发方式和开发规模需要进行优化调整。建议该流域在开发过程中，严格控制农业用地向工业用地和居住用地转化，土地开发类型和开发规模以不超过开发利用承载力为最低要求，达到保护长江水体的目的。当前第1页12