CN110287615B - 基于遥感解译和降雨实验的雨水径流污染负荷测算方法 - Google Patents

基于遥感解译和降雨实验的雨水径流污染负荷测算方法 Download PDF

Info

Publication number
CN110287615B
CN110287615B CN201910574359.2A CN201910574359A CN110287615B CN 110287615 B CN110287615 B CN 110287615B CN 201910574359 A CN201910574359 A CN 201910574359A CN 110287615 B CN110287615 B CN 110287615B
Authority
CN
China
Prior art keywords
rainfall
runoff
remote sensing
experiment
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910574359.2A
Other languages
English (en)
Other versions
CN110287615A (zh
Inventor
袁增伟
季佳莹
张玲
盛虎
刘欣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University
Original Assignee
Nanjing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University filed Critical Nanjing University
Priority to CN201910574359.2A priority Critical patent/CN110287615B/zh
Publication of CN110287615A publication Critical patent/CN110287615A/zh
Application granted granted Critical
Publication of CN110287615B publication Critical patent/CN110287615B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，属于污染负荷测算领域。本发明的方法包括设计降雨径流污染物总负荷核算公式
Figure DDA0002111683020000011
Figure DDA0002111683020000012
收集处理基础数据；开展遥感解译；通过降雨观测实验进行本地化降雨产流模拟；采集分析径流污染物浓度；进行污染负荷测算并开展分析；构建研究区域高空间分辨率污染物排放清单，精准把握径流污染源和污染负荷信息，为研究区精准治污工程及海绵城市建设提供依据。

Description

基于遥感解译和降雨实验的雨水径流污染负荷测算方法
技术领域
本发明属于污染负荷测算领域，更具体地说，涉及一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法。
背景技术
我国地表水环境质量呈逐步改善趋势，但城市黑臭水体现象却仍然突出，不仅损害了城市人居环境，也严重影响城市形象。城市黑臭水体整治工作可以从点源污染和面源污染两个角度开展。相对来说，点源污染的排查和核算比较简单且容易控制，通过几年的治理已经卓有成效，而面源来源广、成分复杂，其中地表径流污染在很多地区已经成为水体水质恶化、河流生态系统退化的重要原因，随着城市化进程加速以及暴雨现象加剧，越来越多地区正在面临同样的问题，为切实有效地改善表水环境污染状况，有必要对径流污染重点控制。
目前，国内的非点源测算大多直接使用国外模型，使用观测值进行参数率定以及模型校正。经检索发现，申请公布号CN107480812A，公布日为2017年12月15日的中国发明专利申请公开了一种城市小流域初期雨水污染负荷预测方法，该方法包括划分小流域河网与汇水区域；沿主河道与上下游子汇水区域的节点处，设置监测断面；监测旱季断面水质、水量；在监测断面处，监测降雨初期t时段内平均水质、水量，计算得各汇水区域的初期雨水污染负荷及污染物浓度；划分土地利用类型，对汇水区域遥感图进行土地解译，得到各汇水区域不同地域类别面积占比；建立初雨污染物浓度与地域类别和降雨强度的评估模型，并通过拟合得到污染系数和雨水冲刷系数；预测城市小流域区域污染趋势变化。该发明虽然获得的实测结果与模拟结果相对一致，但该发明构建的是经验模型，影响模拟结果的未知因素较多，无法保证更多事件下的适应性。
也有研究者根据区域的特性开发本地化的模型，但在模型的应用过程中存在参数相关性、不稳定性和不唯一性的问题，且空间分辨率不高，通常只能识别子流域和用地板块，无法根据模型结果很好的追溯污染物的来源，进而制定有效的治理措施。
此外，在产流模型中，SCS径流曲线模型是基于美国的自然地理状况和水文气象条件，通过大量的实验研究和统计分析得到，涉及参数较少、资料易获取、计算过程不复杂，且有效考虑到了不同下垫面的区别，但由于我国与美国的土壤地形条件及降水特征的明显差异，该SCS径流曲线模型不能直接运用到我国的流域产流模拟。一些研究者也通过以下途径对模型进行修正：包括重新率定径流曲线数CN值(许彦等，基于Arc View的SCS模型在流域径流计算中的应用，《水土保持研究》，2006，4：176-179；罗利芳等，径流曲线数法在黄土高原地表径流量计算中的应用，《水土保持通报》，2002，3：58-61)、寻求AMC更合理的划分 (王少平等，基于GIS的苏州河非点源污染的总量控制，《中国环境科学》，2002，22(6)： 520-524)、将SCS模型改进为栅格式水文分布模型应用于流域研究(罗鹏等，基于栅格式SCS 模型的分布式水文模型研究，《武汉大学学报(工学版)》，2011，44(2)：156-160)等，但这些研究最终没能有效地应用于径流污染源解析工作中。
发明内容
1.要解决的问题
针对现有径流污染的研究中模型构建不准确、无法有效追溯污染源的问题，本发明提供了一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，该方法包括设计降雨径流污染物总负荷核算公式
Figure BDA0002111681000000021
通过高精度遥感解译技术进行土地利用分析；开展降雨观测实验进行本地化降雨产流模拟；监测分析径流污染物浓度；测算并分析雨水径流污染负荷；构建研究区域高空间分辨率污染物排放清单，精准把握径流污染源和污染负荷信息，为研究区域精准治污工程设计及海绵城市建设提供科学依据。
2.技术方案
为了解决上述问题，本发明所采用的技术方案如下：
本发明的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，包括以下步骤：
(1)设计降雨径流污染物总负荷核算公式，
Figure BDA0002111681000000022
(2)收集处理研究区域降雨数据；
(3)使用遥感解译技术进行土地利用分析；
(4)开展降雨观测实验进行本地化降雨产流模拟；
(5)采集分析径流污染物浓度；
(6)雨水径流污染物负荷测算与分析；
(7)构建研究区域高空间分辨率污染物排放清单；
其中：Wj为第j类污染物径流总负荷；Si为第i种土地利用类型的面积；Qi为第i种土地利用类型的降雨径流深；EMCij为第i种土地利用类型上第j类污染物的降雨事件平均浓度； n为流域内土地利用类型数。
优选地，所述步骤(2)的具体步骤为：收集降雨数据，其中降雨数据包括降雨强度、前期干燥天数及研究区域内降雨站经纬度信息，而后将数据导入ARCGIS软件，使用泰森多边形插值分析得到研究区域内雨量分布图层。
优选地，所述步骤(3)的具体步骤为：使用卫星遥感图开展土地利用类型解译，获取各种土地利用类型空间分布，并在ARCGIS软件中分析各土地利用类型面积。
优选地，所述步骤(4)的具体步骤为：开展降雨观测实验，获取各土地利用类型的降雨初损值参数，用降雨初损值参数校正SCS径流曲线模型，而后用校正后的SCS径流曲线模型模拟各降雨事件下的径流产生情况。
优选地，所述步骤(5)的具体步骤为：采集径流样品并监测污染物浓度，通过降雨事件产流时采样的监测数据，结合同步得到的雨强或雨量数据，计算事件平均浓度EMC，
Figure BDA0002111681000000031
其中，EMC为降雨事件中某一污染物的平均浓度(mg/L)；M是整个降雨过程中该类污染物冲刷总量(mg)；V是径流产生总量(L)；t为产流时间(min)；Ct为污染物浓度随时间变化情况(mg/L)；Qt为产流量随时间变化情况(L)；n为整个降雨事件中径流取样次数；Qi为第i次取样时的产流量(L)；Ci为第i次取样时的污染物含量(mg/L)，并且针对各土地利用类型，径流浓度分为夏半年的径流浓度和冬半年的径流浓度。
优选地，所述步骤(6)的具体步骤为：基于编程技术，对污染负荷开展测算，根据测算结果分析污染负荷的分布情况。
优选地，收集降雨数据后，还对收集到的降雨数据进行核查，核查包括验证数据是否完整，查找异常数据，随机选取站点数据较差比较。
优选地，在进行解译之前对影像进行预处理，预处理包括控制点选择、影像配准、影像镶嵌、影像裁剪，而后确定土地利用的划分体系，确定解译的标志，比较解译中的分类方法，在ENVI软件中进行分类，最后在ARCGIS软件中对解译结果目视解译调整，并且分析遥感解译结果。
优选地，开展降雨观测实验，其中在雨季时，搭建径流场来开展天然降雨观测实验，观察天然降雨的产流情况；在非雨季时，使用便携式人工模拟降雨系统来开展模拟降雨观测实验。
3.有益效果
相比于现有技术，本发明的有益效果为：
(1)本发明的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，包括设计了降雨径流污染物总负荷核算公式，
Figure BDA0002111681000000032
该公式使用SCS径流曲线模型开展模拟，进而模拟径流量，此外本发明还包括开展降雨观测实验进行本地化降雨产流模拟，在获得径流污染物浓度EMC的同时，观测各用地类型上的降雨初损值，用以校正SCS 径流曲线模型，由此得到的模拟结果可靠性更强；
(2)本发明的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，使用编程技术实现污染负荷的自动测算，并且通过构建降雨径流高分辨率污染排放清单能够更直观的模拟研究区域径流产生情况，为日后研究区域的精准治污工程及海绵城市建设提供基础数据和指导；
(3)本发明的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，针对用地类型多样的综合性流域设计，不仅可以应用于多种用地类型的流域，也可以应用于单一用地类型的流域，本发明方法的应用范围较广。
附图说明
图1为本发明的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法的流程示意图；
图2为本发明的实施例的降雨量空间分布图；
图3为本发明的实施例的遥感解译后的土地利用类型分布图；
图4为本发明的实施例的流域内径流量日变化图。
具体实施方式
下面结合具体实施例对本发明进一步进行描述。
本发明的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，包括以下步骤：
(1)设计降雨径流污染物总负荷核算公式，
Figure BDA0002111681000000041
其中：Wj为第j类污染物径流总负荷；Si为第i种土地利用类型的面积；Qi为第i种土地利用类型的降雨径流深；EMCij为第i种土地利用类型上第j类污染物的降雨事件平均浓度；n为流域内土地利用类型数。
(2)收集处理研究区域降雨数据，获取由气象局提供的各降雨站点每日降雨信息，其中降雨数据包括降雨深度、前期干燥天数以及研究区域内降雨站经纬度信息；收集降雨数据后，还对收集到的降雨数据进行核查，核查包括验证数据是否完整，查找异常数据，随机选取站点数据较差比较，以确保数据的质量以及后续计算的准确性；而后将数据导入ARCGIS软件，使用泰森多边形插值分析得到研究区域内雨量分布图层。
需要说明的是，验证数据是否完整是指对获取的数据进行排查，查看数据是否有缺失，是否规律呈现；查找异常数据指是指以往年当地当场降雨雨量范围为依据，查找超出正常降雨量范围的降雨量数据；随机选取站点数据较差比较是指为了进一步核实数据质量，挑选降雨量较大的数据进行核实，一方面，参照气象网提供的降雨量预测数据进行判断，另一方面，所有站点监测的日降雨数据相互之间虽有差异，但必然有一致性，所以用站点的日降雨量数据进行相互验证。
(3)使用卫星遥感图开展土地利用类型的解译，获取土地利用类型，并在ARCGIS软件中分析汇总解译结果。
需要说明的是，国土规划部门的土地利用现状图以及开源的Landsat TM影像解译结果图分辨率一般较低，因此需要高分辨率商业遥感影像图开展土地利用类型的解译，保证最终解译后图层的空间分辨率不低于2m×2m。
还需要进一步说明的是，在进行解译之前需要对影像进行预处理，预处理包括控制点选择、影像配准、影像镶嵌、影像裁剪，而后确定土地利用的划分体系，确定解译的标志，比较解译中的分类方法，在ENVI软件中进行分类，最后在ARCGIS软件中对解译结果目视解译调整，并且分析遥感解译结果。
(4)开展降雨观测实验，获取各土地利用类型的降雨初损值Ia参数，用降雨初损值Ia 参数校正SCS径流曲线模型，而后用校正后的SCS径流曲线模型模拟各降雨事件下的径流产生情况。
需要指出的是，在开展降雨观测实验的过程中，观测开始降雨至开始产流这一段时间中的降雨量作为地面截留量，以代表特定前期土壤湿度条件下降雨产流SCS曲线模型中各用地类型的降雨初损值Ia，其中，径流产生的标志为径流场内地表有明显的雨水成股从堰口流出。
此外，开展降雨观测实验，其中在雨季时，搭建径流场来开展天然降雨观测实验，观察天然降雨的产流情况；在非雨季时，为了缩短实验周期，使用便携式人工模拟降雨系统来开展模拟降雨观测实验。
值得说明的是，为了更好地观测径流的形成，需要在非硬质土地利用类型的实验场地设置径流场，径流场的围建采用不锈钢铁皮材料，将四块铁皮首尾相接，只留坡度最低处的一口用于径流的流出，用铁锤将其砸入地表，用于阻断地表径流向周围地面的外流，保证有明显的地表径流从径流场出口流出，由于硬质地面由于降雨产流较为明显，易于肉眼判断，并且径流场搭建会破坏地表状态，不易恢复，因此不设径流场。
此外，值得进一步说明的是，在天然降雨观测实验中，使用玻璃材质的雨量杯来收集降雨，选用的雨量杯精度不低于0.1mm，使用雨量杯的具体操作为：在降雨前挑选周围开阔，无建筑物遮挡的地方安置雨量杯于实验场地，待判断地表产流后，立即停止收集降雨，对雨量杯进行读数，该读数即为该种用地类型在此前期土壤湿度条件下的降雨初损值。
在模拟降雨观测实验中，由于模拟降雨器自带雨量感测系统，系统可以导出连续的雨强数据，只需要在模拟降雨实验中记录地表开始产流的时刻，再与导出的降雨时刻数据相对应，即可通过这部分数据处理获得该类用地上的降雨初损值。
(5)采集径流样品并监测污染物浓度，通过各场降雨事件产流时连续采样的监测数据，结合同步得到的雨强或雨量数据，计算事件平均浓度EMC，
Figure BDA0002111681000000061
其中，EMC为降雨事件中某一污染物的平均浓度(mg/L)；M是整个降雨过程中该类污染物冲刷总量(mg)；V是径流产生总量(L)；t为产流时间(min)；Ct为污染物浓度随时间变化情况(mg/L)；Qt为产流量随时间变化情况(L)；n为整个降雨事件中径流取样次数；Qi为第i次取样时的产流量(L)；Ci为第i次取样时的污染物含量(mg/L)，并且针对各土地利用类型，径流浓度分为夏半年的径流浓度和冬半年的径流浓度。
需要说明的是，径流污染物的浓度随降雨时间变化波动较大，存在一次或二次冲刷现象，对单场降雨事件引起的径流污染负荷评价中需要对径流浓度做出总体评价，因此以每个样品对应的径流量为权重使用降雨事件平均浓度EMC来表示在单场降雨径流全过程中排放的某种污染物的平均浓度。
还需要说明的是，由于污染物干湿沉降以及人类活动的差异会导致各用地类型在不同季度的污染物累积不同，因此需要划分径流浓度代表的季度，然后设计监测时间和采样点位，选取降雨事件开展降雨监测，并根据研究区域的污染现状选择采样监测指标。
径流采样点位需要按照以下原则设置：①该点所代表的土地利用类型明确且单一，该土地类型上的生长物不与其他土地利用类型上的物质相混合；②该点所代表的土地利用类型能确保在遥感解译中正确解译出来，在采样点设置后中需要在解译获得的土地利用图层中进行验证；③采样点的设置过程中需要提前前往现场进行踩点，确保点位上没有异常人为活动，如管道检修，现场施工活动等；④确保采样的硬质土地利用类型上有能打开的雨水井，非硬质的土地利用类型上可以依照天然坡度挖渠引流。
(6)绘制污染负荷核算分析的逻辑框架，运用编程技术实现测算过程自动化并开展分析，获取径流污染负荷核算结果以及时空分布图。
(7)构建研究区域高空间分辨率污染物排放清单，在时间精度上实现日尺度的分析，空间精度上实现高分辨率网格的展现，为研究区精准治污工程及海绵城市建设提供依据。
实施例1
本实施例以2016年十五里河流域径流污染源时空解析为研究对象。
十五里河是合肥市区重要的城市内部景观河流，位于合肥市西南郊，自西北流向东南汇入巢湖，加上塘西河转输雨水汇水面积，流域总面积135.8km2，河道全长约22.9km，为合肥市西南部主要行洪通道之一。该流域内无矿山开采，无水田排水，农村养殖也几乎完全取缔，流域的径流污染主要来自于降雨径流。
(1)设计降雨径流污染物总负荷核算公式
根据十五里河前期预调研的情况了解，流域内各土地用地类型产流情况差别较大，径流污染负荷应分别计算，流域内径流污染总负荷可由各用地类型的降雨“径流深”、“污染物浓度”与“用地面积”这三部分相乘获得，即：
Figure BDA0002111681000000071
其中：Wj为第j类污染物径流总负荷；Si为第i种用地类型的面积；Qi为第i种用地类型的降雨径流深；EMCij为第i种用地类型上第j类污染物的事件平均浓度；n为流域内用地类型数。
(2)收集处理基础资料数据
如图1所示，收集2016年流域内日降雨数据，该数据由合肥市气象局提供，包括合肥市内所有降雨观测站的降雨信息以及点位信息，流域内共计6个降雨观测站，分别为十五里河站、G小西冲站、包河区政府站、合肥站、宿松路站、义城站。通过将这些点位数据导入ARCGIS 软件，与流域范围图层做相交处理，获取了流域内的观测点位。
在获取降雨观测站的日降雨数据后，对这些数据进行核查，以确保数据的质量以及后续计算的准确性，核查包括：①验证数据是否完整；②查找异常数据；③随机选取站点数据较差比较，核实选择的数据。本实施例所获取的十五里河流域2016年6个站点日降雨数据质量高，可用于后续的分析和计算。
由于6个降雨观测站在流域内呈散点分布，而降雨实际上呈面状分布，因此需要对流域以这6个降雨观测站所代表的范围进行划分，在ARCGIS工具中插值分析模块可以根据本流域内6个降雨站点的经纬度，直接生成每个站点代表范围的泰森多边形，得到十五里河流域内2016年降雨量分布，如图2所示。
(3)开展遥感解译
使用商业化的高分辨率卫星遥感图参数(见表1)，开展解译，获取土地利用类型，在 ARCGIS软件中分析汇总解译结果。
表1遥感影像图参数
Figure BDA0002111681000000072
Figure BDA0002111681000000081
获得遥感影像图之后，需要对影像进行预处理，然后参考国家最新土地利用现状分类标准，结合流域的实际情况和研究需要，将土地类型分为水域、道路、工业区、商业区、小区及行政区、农田、林地及果园、未建设裸地、灌木及草地这9个用地类型；制定解译标志，并参照Google地图进行验证；比较各种分类方法，探索优缺点后，在ENVI5.2软件中采用监督分类的方法划分土地利用类型；进行目视解译调整，最终获得十五里河流域土地利用类型图，如图3所示。
(4)开展降雨观测实验，校正SCS径流曲线模型并开展模拟
开展降雨观测实验，获取各用地类型的降雨初损值参数，用校正后的SCS径流曲线模型模拟各场降雨下径流产生情况。
结合本实施例中降雨径流观测实验开展的进度安排，以及各类用地类型的地表径流特征，不同用地类型的降雨径流观测方式如表2所示：
表2不同用地类型径流观测方式
Figure BDA0002111681000000082
本实施例中购买NLJY-10型便携式人工模拟降雨系统，各性能参数如表3所示，其性能在降雨均匀度系数以及降雨调节精度等参数方面，优于市场上其他模拟降雨系统，能保证降雨观测数据的可靠性。
表3模拟降雨器性能参数
Figure BDA0002111681000000091
实验结果发现，流域内降雨初损最大的用地类型为农田，在前期湿度为I级、II级、III 级的条件下，降雨初损分别达到了62.36mm、26.19mm和11.39mm，在6个降雨观测站点，2016年分别只有12、17、7、15、13、10日的降雨在农田用地上产生径流。而降雨初损相对最小的工业区和商业区，在三个前期湿度级别下，降雨初损分别为5mm、2.1mm和2.91mm。
本实施例中基于上述数据，流域内径流量日变化情况如图4所示。
(5)采集径流样品并监测污染物浓度
通过各场降雨事件产流时连续采样的监测数据计算事件平均浓度(EMC)，各用地类型区分夏半年与冬半年径流浓度。各用地类型的代表性降雨径流事件如表4所示，各降雨事件采样点位及所代表的土地利用类型如表5所示，同步监测流量与污染物浓度，获得夏半年与冬半年的各用地类型的EMC，根据研究区域的污染现状，本实施例中选择了COD、NH4 +-N、TP、TN四类污染物进行监测。
表4各用地类型代表性径流事件时间
Figure BDA0002111681000000092
Figure BDA0002111681000000101
表5径流采样点位经纬度
Figure BDA0002111681000000102
Figure BDA0002111681000000111
(6)污染负荷测算及分析
绘制核算分析的逻辑框架，使用R语言在RStudio平台中编程，开展径流污染负荷测算，输入1、径流浓度EMC.xlsx，2、降雨空间站空间插值.shp，3、分区日降雨量.xlsx，4、各用地类型初损值.xlsx，5、土地利用图层.tif这五个文件，运行程序，获取十五里河流域2016年径流污染负荷核算结果以及时空分布图。
(7)构建高空间分辨率降雨径流污染物排放清单
构建研究区域高空间分辨率降雨径流污染物排放清单，在时间精度上实现日尺度的分析，空间细度上实现高分辨率网格的展现，为十五里河流域径流污染后续治理提供支撑，也为海绵城市建设提供基础数据和指导。
本发明的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法结合SCS径流曲线模型的本地化降雨观测实验、高精度遥感解译及径流污染采样监测，开展污染负荷测算和分析，构建高空间分辨率径流污染源排放清单，为研究区精准治污工程及海绵城市建设提供依据。

Claims (9)

1.一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，包括以下步骤：
(1)设计降雨径流污染物总负荷核算公式，
Figure FDA0002574627110000011
(2)收集处理研究区域降雨数据；
(3)通过遥感解译进行土地利用分析；
(4)开展降雨观测实验进行本地化降雨产流模拟；
(5)监测分析雨水径流污染物浓度；
(6)雨水径流污染物负荷测算与分析；
(7)构建研究区域高空间分辨率污染物排放清单；
其中：Wj为第j类污染物径流总负荷；Si为第i种土地利用类型的面积；Qi为第i种土地利用类型的降雨径流深；EMCij为第i种土地利用类型上第j类污染物的降雨事件平均浓度；n为流域内土地利用类型数。
2.根据权利要求1所述的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，其特征在于：所述步骤(2)的具体步骤为：收集降雨数据，其中降雨数据包括降雨强度、前期干燥天数及研究区域内降雨站经纬度信息，而后将数据导入ARCGIS软件，使用泰森多边形插值分析得到研究区域内雨量分布图层。
3.根据权利要求1所述的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，其特征在于：所述步骤(3)的具体步骤为：使用卫星遥感图开展土地利用类型解译，获取研究区各种土地利用类型分布，并在ARCGIS软件中分析各类土地利用面积。
4.根据权利要求1所述的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，其特征在于：所述步骤(4)的具体步骤为：开展降雨观测实验，获取各土地利用类型的降雨初损值参数，用降雨初损值参数校正SCS径流曲线模型，而后用校正后的SCS径流曲线模拟各降雨事件下的径流产生情况。
5.根据权利要求1所述的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，其特征在于：所述步骤(5)的具体步骤为：采集径流样品并监测污染物浓度，通过降雨事件产流时采样的监测数据，结合同步得到的雨强或雨量数据，计算事件平均浓度EMC，
Figure FDA0002574627110000012
其中，EMC为降雨事件中某一污染物的平均浓度，其单位为mg/L；M是整个降雨过程中该类污染物冲刷总量，其单位为mg；V是径流产生总量，其单位为L；t为产流时间，其单位为min；Ct为污染物浓度随时间变化情况，其单位为mg/L；Qt为产流量随时间变化情况，其单位为L；n为整个降雨事件中径流取样次数；Qi为第i次取样时的产流量，其单位为L；Ci为第i次取样时的污染物含量，其单位为mg/L。
6.根据权利要求1所述的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，其特征在于：所述步骤(6)的具体步骤为：基于编程技术，对污染负荷开展测算，根据测算结果分析污染负荷的分布情况。
7.根据权利要求2所述的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，其特征在于：收集降雨数据后，还对收集到的降雨数据进行核查，核查包括验证数据是否完整，查找异常数据，随机选取站点数据较差比较。
8.根据权利要求3所述的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，其特征在于：在进行解译之前对影像进行预处理，预处理包括控制点选择、影像配准、影像镶嵌、影像裁剪，而后确定土地利用的划分体系，确定解译的标志，在ENVI软件中进行分类，最后在ARCGIS软件中对解译结果目视解译调整，并且分析遥感解译结果。
9.根据权利要求4所述的一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法，其特征在于：开展降雨观测实验，
其中在雨季时，搭建径流场来开展天然降雨观测实验，观察天然降雨的产流情况；
在非雨季时，使用便携式人工模拟降雨系统来开展模拟降雨观测实验。
CN201910574359.2A 2019-06-28 2019-06-28 基于遥感解译和降雨实验的雨水径流污染负荷测算方法 Active CN110287615B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910574359.2A CN110287615B (zh) 2019-06-28 2019-06-28 基于遥感解译和降雨实验的雨水径流污染负荷测算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910574359.2A CN110287615B (zh) 2019-06-28 2019-06-28 基于遥感解译和降雨实验的雨水径流污染负荷测算方法

Publications (2)

Publication Number Publication Date
CN110287615A CN110287615A (zh) 2019-09-27
CN110287615B true CN110287615B (zh) 2020-08-28

Family

ID=68019489

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910574359.2A Active CN110287615B (zh) 2019-06-28 2019-06-28 基于遥感解译和降雨实验的雨水径流污染负荷测算方法

Country Status (1)

Country Link
CN (1) CN110287615B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110675091B (zh) * 2019-10-16 2022-08-09 中国水利水电科学研究院 一种城市地表径流污染物负荷计算方法
CN111310124A (zh) * 2020-02-13 2020-06-19 北京市环境保护科学研究院 一种城市径流数据处理方法
CN112163347B (zh) * 2020-10-21 2021-09-07 中国科学院生态环境研究中心 一种城市地块尺度面源污染动态输出特征评估的方法
CN114970388B (zh) * 2022-05-12 2023-07-25 河海大学 一种自适应时间步长的场地地表污染扩散方法
CN115168446B (zh) * 2022-06-30 2023-04-07 中科三清科技有限公司 污染物溯源方法、装置及电子设备
CN115267950B (zh) * 2022-08-26 2023-07-14 中国长江三峡集团有限公司 一种初期雨水量化方法、装置、存储介质及电子设备

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7062406B2 (en) * 2003-11-13 2006-06-13 Ch2M Hill, Inc. Method and system for water quality trading
CN102867120B (zh) * 2012-09-05 2016-02-24 环境保护部卫星环境应用中心 一种基于遥感像元的非点源污染计算方法
CN103020424B (zh) * 2012-11-22 2016-04-27 北京师范大学 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法
CN106528788B (zh) * 2016-11-09 2022-04-01 上海市政工程设计研究总院（集团）有限公司 基于gis技术分析地面雨径流污染空间分布特征的方法
CN106682272A (zh) * 2016-12-05 2017-05-17 北京工业大学 一种确定swmm中含不同下垫面地表径流污染物参数w的方法
CN108733915B (zh) * 2018-05-17 2019-08-20 河海大学 基于栅格化处理的平原河网区河道面源污染负荷确定方法
CN108717453B (zh) * 2018-05-21 2020-05-05 河海大学 一种基于gis平台的平原河网区污染负荷计算方法

Also Published As

Publication number Publication date
CN110287615A (zh) 2019-09-27

Similar Documents

Publication Publication Date Title
CN110287615B (zh) 基于遥感解译和降雨实验的雨水径流污染负荷测算方法
Demirel et al. Combining satellite data and appropriate objective functions for improved spatial pattern performance of a distributed hydrologic model
Singh et al. Hydrological modeling of the Iroquois river watershed using HSPF and SWAT 1
Hasenauer et al. Validating daily climate interpolations over complex terrain in Austria
Lin et al. Assessing temporal trends of soil erosion and sediment redistribution in the Hexi Corridor region using the integrated RUSLE-TLSD model
Lasserre et al. A GIS-linked model for the assessment of nitrate contamination in groundwater
CN108280789B (zh) 一种精细时空尺度二氧化碳排放要素间影响的空间分析方法
Zhang et al. Runoff variation and its response to climate change in the Three Rivers Source Region
CN111428972A (zh) 基于swmm与efdc耦合模型的调蓄工程环境效应评估方法及装置
CN112785024B (zh) 一种基于流域水文模型的径流计算和预测方法
Caiqiong et al. Application of HYDRUS-1D model to provide antecedent soil water contents for analysis of runoff and soil erosion from a slope on the Loess Plateau
CN112101693B (zh) 基于正交分析的现状城区河道水质达标分析方法
Wu et al. The simulation research of dissolved nitrogen and phosphorus non-point source pollution in Xiao-Jiang watershed of Three Gorges Reservoir area
Hammouri et al. Hydrological modeling of ungauged wadis in arid environments using GIS: a case study of Wadi Madoneh in Jordan
Jiang et al. A framework to develop a watershed pollution load model for semiarid and semihumid areas
CN115859823A (zh) 一种基于滑动窗口子时期划分的时变参数水文模拟方法
CN116205136A (zh) 基于径流滞后信息的大尺度流域深度学习洪水预报方法
Liu et al. Establishing the relationship between the integrated multidimensional landscape pattern and stream water quality in subtropical agricultural catchments
Wu et al. Simulation of nitrogen and phosphorus loads in the Dongjiang River basin in South China using SWAT
Jia et al. Coupling simulation of water and energy budgets and analysis of urban development impact
Hassan et al. Statistical downscaling for climate change scenarios of rainfall and temperature
Coville et al. Loss of street trees predicted to cause 6000 L/tree increase in leaf-on stormwater runoff for Great Lakes urban sewershed
Farzana et al. Application of swat model for assessing water availability in Surma River basin
Mostafazadeh et al. Plot-size runoff and sediment yield affected by the type of geological formations and implications for water erosion in a semi-arid region
CN113361114B (zh) 一种基于径流路径的多尺度面源污染物入河系数测算方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant