CN113139757A - 污染源排污量化分析方法及装置 - Google Patents
污染源排污量化分析方法及装置 Download PDFInfo
- Publication number
- CN113139757A CN113139757A CN202110545811.XA CN202110545811A CN113139757A CN 113139757 A CN113139757 A CN 113139757A CN 202110545811 A CN202110545811 A CN 202110545811A CN 113139757 A CN113139757 A CN 113139757A
- Authority
- CN
- China
- Prior art keywords
- pollution source
- monitoring section
- pollution
- control area
- exit
- 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.)
- Pending
Links
- 238000004445 quantitative analysis Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012544 monitoring process Methods 0.000 claims abstract description 139
- 238000004364 calculation method Methods 0.000 claims abstract description 42
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 40
- 231100000719 pollutant Toxicity 0.000 claims abstract description 40
- 238000012163 sequencing technique Methods 0.000 claims abstract description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 63
- 238000003860 storage Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000011144 upstream manufacturing Methods 0.000 claims description 9
- 238000011002 quantification Methods 0.000 claims description 8
- 238000012821 model calculation Methods 0.000 claims description 4
- 238000009826 distribution Methods 0.000 abstract description 8
- 238000009792 diffusion process Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 230000007613 environmental effect Effects 0.000 description 9
- 239000010865 sewage Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 5
- 229910052698 phosphorus Inorganic materials 0.000 description 5
- 239000011574 phosphorus Substances 0.000 description 5
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 230000031018 biological processes and functions Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011109 contamination Methods 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 239000002893 slag Substances 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005538 encapsulation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000009432 framing Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000004856 soil analysis Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 239000002351 wastewater Substances 0.000 description 1
- 238000004065 wastewater treatment Methods 0.000 description 1
- 238000004457 water analysis Methods 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种污染源量化分析方法及装置。该方法包括基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分;获取目标区域中各污染源所在位置,并确定各污染源所属的控制区域;根据各监测断面相对于所在河流的上下游关系,确定各控制区域的入境监测断面和出境监测断面;计算各污染源与其所在控制区域的出境监测断面间的距离;对于各污染源,计算污染源的各种污染物到达出境监测断面处的污染物浓度;将污染源的各种污染物到达出境监测断面处的浓度进行加权平均计算并对计算结果排序,得到各污染源排污量化分析结果。在控制区域的划分上,结合流域自然地理特征和监测断面的位置,且将量化分析的特性更多地放在权重分配上。
Description
技术领域
本发明涉及水环境技术领域,更为具体来说,本发明涉及一种污染源排污量化分析方法及装置。
背景技术
在对污染源的评估和管理方面,目前主要采用一种点源突发性水污染源事故溯源贡献程度计算方法,通过计算环境容量和污染源排放贡献率的方式来辅佐决策。前者根据流域自然地理特征确定流域控制断面对应的控制区域,得出了目标区域内允许排放的污染负荷(质量单位)。后者通过直接利用污染源排放量监测数据或模拟某个污染源的排放扩散(或流量*浓度的结果)来评价该污染源对断面监测水质的贡献程度,即该污染源为目标区域内的水环境贡献了多少污染物。
现有的方案大多都是根据流域自然地理特征比如自然地势来确定控制区域(子流域或汇水范围)。这种划分在应用中存在的缺点是:环境容量范围的框定和污染源脱节。一般而言,环境容量是根据流域汇水区域、功能区划来划分的,常会导致一个环境容量范围内包含多个行政区划或者划分区域内缺少控制单元作为出入境监测断面。前者不便于政务管理工作的开展,后者会因缺少上游来水的水质数据导致结果失真。
以贡献率的方式能够计算出污染源排放在量上对环境污染的贡献程度,但却不能反映出其排放对环境的威胁性。比如某家大型污水厂的排水量大,但浓度符合规范,其整体上给该河流贡献了十分之一的污染物。某小型企业排水量小,但浓度很高,其排放对河流的贡献程度很低,但高浓度的污水却直接威胁到了水生态的安全。同时,从管控的角度来看,如果上述大型污水处理厂离出境断面较远,小型企业离得很近,则该小型企业的高浓度排放更容易导致出境断面处的监测水质频频超标;而大型企业由于距离远,排水中污染物衰减比例高,反而对下游断面的监测水质影响很小。不探明这些原因,会对该区域内经济生产活动的部署安排产生不便。
在水环境管理中,对于排污企业(污染源)的管理一直是重中之重。而如何摒弃一刀切的限排策略,在科学支撑下,对污染排放企业做到精准定位、合理管理,一直是该项工作中的难点。量化分析在传统意义上更多指的是将难以定量的指标通过数据或区间的方式为其赋值,再通过对量化值的统计来得出一个可量化比较的结论。
由此,目前亟需一种对污染源排污进行量化分析的方法,结合地理特征、监测断面位置部署以及行政区划的控制区域划分方式,保证每个控制区域在汇水上具有合理性的同时,具有上下游的监测断面和较少的行政单元,便于计算评定和应用管控,而且将量化分析的特性更多地放在了权重分配上,并且通过不同指标对应权重的分配来体现对各个指标关注程度的不同。
发明内容
本发明创新地提供了一种污染源排污量化分析方法及装置,在控制区域的划分上,结合了流域自然地理特征和监测断面的空间位置,而且将量化分析的特性更多地放在了权重分配上。
为实现上述的技术目的,一方面,本发明公开了一种污染源排污量化分析方法。所述污染源排污量化分析方法包括:基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分;获取目标区域中各个污染源所在的位置,并确定各个污染源所属的控制区域;根据各个监测断面相对于所在河流的上下游关系,确定各个控制区域的入境监测断面和出境监测断面;计算各个污染源与其所在的控制区域的出境监测断面之间的距离;对于各个污染源,根据所述污染源与其所在的控制区域的出境监测断面之间的距离,计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度;将所述污染源的各种污染物到达所述出境监测断面处的浓度进行加权平均计算并得到加权平均计算结果,对各个污染源的加权平均计算结果进行排序,从而得到各个污染源排污量化分析结果。
进一步地,对于所述污染源排污量化分析方法,基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分包括:基于目标区域的自然地理特征、监测断面所在位置和行政区域划分,对目标区域进行控制区域划分。
进一步地,对于所述污染源排污量化分析方法,根据所述污染源与其所在的控制区域的出境监测断面之间的距离,计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度,包括:计算到达所述污染源之前的混合水质;基于菲克定律的一维河流水质模型得到衰减模型;根据所述污染源与其所在的控制区域的出境监测断面之间的距离、以及到达所述污染源之前的混合水质,利用所述衰减模型计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度。
进一步地,对于所述污染源排污量化分析方法,计算到达所述污染源之前的混合水质包括:利用零维完全混合模型计算到达所述污染源之前的混合水质。
为实现上述的技术目的,另一方面,本发明公开了一种污染源排污量化分析装置。所述污染源排污量化分析装置包括:控制区域划分单元,用于基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分;污染源和控制区域关联单元,用于获取目标区域中各个污染源所在的位置,并确定各个污染源所属的控制区域;监测断面和控制区域关联单元,用于根据各个监测断面相对于所在河流的上下游关系,确定各个控制区域的入境监测断面和出境监测断面;距离计算单元,用于计算各个污染源与其所在的控制区域的出境监测断面之间的距离;污染源到出境处水质计算单元,用于对于各个污染源,根据所述污染源与其所在的控制区域的出境监测断面之间的距离,计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度;量化分析单元,用于将所述污染源的各种污染物到达所述出境监测断面处的浓度进行加权平均计算并得到加权平均计算结果,对各个污染源的加权平均计算结果进行排序,从而得到各个污染源排污量化分析结果。
进一步地,对于所述污染源排污量化分析装置,所述控制区域划分单元进一步用于基于目标区域的自然地理特征、监测断面所在位置和行政区域划分,对目标区域进行控制区域划分。
进一步地,对于所述污染源排污量化分析装置,所述污染源到出境处水质计算单元包括:到达污染源前水质计算模块,用于计算到达所述污染源之前的混合水质;衰减模型确定模块,用于基于菲克定律的一维河流水质模型得到衰减模型;衰减模型计算模块,用于根据所述污染源与其所在的控制区域的出境监测断面之间的距离、以及到达所述污染源之前的混合水质,利用所述衰减模型计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度。
进一步地,对于所述污染源排污量化分析装置,所述到达污染源前水质计算模块进一步用于利用零维完全混合模型计算到达所述污染源之前的混合水质。
为实现上述的技术目的,又一方面,本发明公开了一种计算设备。所述计算设备包括:一个或多个处理器,以及与所述一个或多个处理器耦合的存储器,所述存储器存储指令,当所述指令被所述一个或多个处理器执行时,使得所述一个或多个处理器执行上述方法。
为实现上述的技术目的,再一方面,本发明公开了一种机器可读存储介质。所述机器可读存储介质存储有可执行指令,所述指令当被执行时使得所述机器执行上述方法。
本发明的有益效果为:
本发明实施例的污染源排污量化分析方法及装置在控制区域的划分上,结合了流域自然地理特征和监测断面的空间位置,最大化实现每个控制区域具有出入境的监测断面,从而使每个控制单元都有有效的监测数据可用,便于计算评定和应用管控,结合政务实现精准管理;而且将量化分析的特性更多地放在了权重分配上,并且通过不同指标对应权重的分配来体现对各个指标关注程度的不同。另一方面,考虑了污染源排污浓度和与下游监测断面之间的距离,以污染物传递机理模型为内核,通过功能性的封装,建立了非稳态扩散模型,通过模拟控制区域内每个污染源排污到达该控制区域时衰减后的污染物浓度,来判断其排放对出境监测断面后的生态环境的威胁性,与环保业务紧密结合,适宜于环保工作。
附图说明
图中,
图1为本发明一个实施例提供的污染源排污量化分析方法的流程图;
图2为本发明另一个实施例提供的污染源排污量化分析装置的结构示意图;
图3为根据本发明实施例的用于污染源排污量化分析处理的计算设备的结构框图。
具体实施方式
下面结合说明书附图对本发明提供的污染源排污量化分析方法及装置进行详细的解释和说明。
图1为本发明一个实施例提供的污染源排污量化分析方法的流程图。
如图1所示,在步骤S110,基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分。
基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分。进一步来说,还可以结合行政区域划分,对控制区域划分进行修正。当然,也可以先考虑目标区域的自然地理特征和行政区划,对目标区域进行控制区域划分;然后再根据监测断面所在位置,对控制区域划分进行调整。
获取目标区域的数字高程模型(DEM,Digital Elevation Model)高程地图,确定监测断面的经纬度点位,利用SWAT(Soil and Water Assessment Tool)模型的水文响应单元划分功能或者ArcGIS的hydrology工具,制作控制区域划分图。在制作过程中,将监测断面所在的点位作为汇水节点,则划分出的控制区域的出入境会以监测断面所在位置区分。其中,ArcGIS是计算机制图应用,包含了全球范围内的底图、地图数据和应用程序,以及可配置的应用模板和开发人员使用的GIS工具和API,可用于创建Web地图、发布GIS服务、共享地图、数据和应用程序,以及管理组织的内容和多个用户。
控制区域划分完毕后,可以在ArcGIS中叠加行政区划图层,利用ArcGIS-ArcToolBox-Analysis Tools-Overlay中的Identity或Union功能,将行政区域划分矢量图和控制区域划分图进行特征合并,根据管理蓝图或河流特征将区域划分图做进一步的修正。比如,随着水环境监管程度不断完善,多地环保厅局都逐步在河流出入各个市县区处部署监测断面。由此,可以根据其未来宏观布局,利用行政区划图将控制区划图做进一步细分。
在步骤S120,获取目标区域中各个污染源所在的位置,并确定各个污染源所属的控制区域。
关联污染源和控制区域。获取包含经纬度信息的目标区域所有污染源清单,可以将整理后的清单数据导入ArcGIS,利用ArcGIS中的Display XY Data功能将其进行空间可视化布点。随后可以利用ArcGIS-ArcToolBox-Analysis Tools-Overlay中的SpatialJoin功能在污染源属性表中给每个污染源添加所属的控制区域字段。
在步骤S130,根据各个监测断面相对于所在河流的上下游关系,确定各个控制区域的入境监测断面和出境监测断面。
关联监测断面和控制区域。获取包含经纬度信息的目标区域中所有监测断面清单,首先可以利用Select Features工具全选所有监测断面,将监测断面向河流上游方向位移,可以利用SpatialJoin功能关联出控制区域的第一列监测断面字段;随后将监测断面向河流下游方向位移,可以利用SpatialJoin功能关联出控制区域的第二列字段。最后,利用控制区域中河流上下游的关系,区分各个监测断面是属于关联控制区域中的出流断面还是入流断面。经过增删改查,最终在控制区域属性表中创建各个控制区域的入境监测断面字段和出境监测断面字段。
在步骤S140,计算各个污染源与其所在的控制区域的出境监测断面之间的距离。
关联污染源和该污染源所在控制区域的上下游监测断面。以控制区域标识(ID)或名称为索引,为污染源信息表单(即上文中的污染源属性表)中添加所在控制区域的出境监测断面和入境监测断面相关字段。
计算并排序污染源与其下游的出境监测断面之间的概化空间距离。至此,在数据集中已经包括各个污染源的经纬度信息和与其相关的下游监测断面的经纬度信息。用这两组经纬度信息,计算出其概化的地表距离。由于污染源多由管道管网排污入河,与地表地形关联较少,所以此处可以忽略地球表面地形对计算带来的误差;同时由于地下管网数据难考,此处在宏观层面上暂时忽略了地下管网分布带来的距离误差。污染源与其下游的出境监测断面之间概化的地表距离的具体计算公式可以采用Google地图(Google Map)计算程式。
建立关联性数据集并接入时序数据。数据集可以包括:关联信息总表、污染源信息表单、污染源排放监测时序数据表、监测断面信息表单、监测断面水文与水质时序数据表、每个控制区域的环境参数方案表单、以及/或者权重方案表单。其中,关联信息总表中的核心字段可以包括:污染源名称、污染源标识(ID)、所属控制区域编号、上游监测断面标识(ID)、下游监测断面标识(ID)、该污染源距下游监测断面的距离、以及/或者排放类型等字段。污染源信息表单和监测断面信息表单包括其元数据字段,包括接入的业务信息和其经纬度。污染源与监测断面的时序数据表中分别接入污染源与监测断面处的实时时序水质监测值与流量值。
在步骤S150,对于各个污染源,根据污染源与其所在的控制区域的出境监测断面之间的距离,计算污染源的各种污染物到达出境监测断面处的污染物浓度。
步骤S150可以具体包括如下过程:计算到达污染源之前的混合水质;基于菲克(Fick)定律的一维河流水质模型得到衰减模型;根据污染源与其所在的控制区域的出境监测断面之间的距离、以及到达污染源之前的混合水质,利用衰减模型计算污染源的各种污染物到达出境监测断面处的污染物浓度。
其中,计算到达污染源之前的混合水质可以包括如下步骤:利用零维完全混合模型计算到达污染源之前的混合水质。
封装基于Fick定律的衰减模型。污染源排污到达目标断面的浓度的核心算法是根据Fick第一定律和第二定律得到的,并在此基础上可以考虑弥散扩散和转化降解的因素,包括沉淀作用、气液交换作用、吸附作用、化学转化和生物作用,分别以弥散系数(Dn)和各降解转换系数(k(Kp,Ke,Ks,Kc,Kb))来体现。Kp表示沉淀系数,Ke表示气液交换系数,Ks表示吸附作用系数,Kc表示化学转化系数,Kb表示生物作用系数。其中,Fick第一定律是在单位时间内通过垂直于扩散方向的单位截面积的扩散物质流量(Diffusion Flow)即扩散通量与该截面处的浓度梯度(Concentration gradient)成正比,也就是说浓度梯度越大,扩散通量越大。菲克(Fick)第二定律是在非稳态扩散过程中,在距离X处,浓度随时间的变化率等于该处的扩散通量随距离变化率的负值。
在步骤S160,将污染源的各种污染物到达出境监测断面处的浓度进行加权平均计算并得到加权平均计算结果,对各个污染源的加权平均计算结果进行排序,从而得到各个污染源排污量化分析结果。
假设污染源的污染物种类有N种,各种污染物达出境监测断面处的浓度为C1,C2,…,CN,假设M1,M2,…,MN分别是浓度C1,C2,…,CN对应的权重数值,各种污染物达出境监测断面处的浓度的加权平均值为L=(C1*M1+C2*M2…+CN*MN)/N。例如,如果对于各种污染物中的氮和磷更加关注,则将氮和磷分别对应的权重数值设置得大于其他污染物对应的权重数值。
在实际应用中,随着实测数据的更新推送,触发执行上述方法流程,来实现动态监管的目的。
图2为本发明另一个实施例提供的污染源排污量化分析装置的结构示意图。如图2所示,该实施例提供的污染源排污量化分析装置200包括控制区域划分单元210、污染源和控制区域关联单元220、监测断面和控制区域关联单元230、距离计算单元240、污染源到出境处水质计算单元250、以及量化分析单元260。
控制区域划分单元210用于基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分。控制区域划分单元210的操作可以参照上面图1描述的步骤S110的操作。作为一种可选实施方式,控制区域划分单元210进一步用于基于目标区域的自然地理特征、监测断面所在位置和行政区域划分,对目标区域进行控制区域划分。
污染源和控制区域关联单元220用于获取目标区域中各个污染源所在的位置,并确定各个污染源所属的控制区域。污染源和控制区域关联单元220的操作可以参照上面图1描述的步骤S120的操作。
监测断面和控制区域关联单元230用于根据各个监测断面相对于所在河流的上下游关系,确定各个控制区域的入境监测断面和出境监测断面。监测断面和控制区域关联单元230的操作可以参照上面图1描述的步骤S130的操作。
距离计算单元240用于计算各个污染源与其所在的控制区域的出境监测断面之间的距离。距离计算单元240的操作可以参照上面图1描述的步骤S140的操作。
污染源到出境处水质计算单元250用于对于各个污染源,根据污染源与其所在的控制区域的出境监测断面之间的距离,计算污染源的各种污染物到达出境监测断面处的污染物浓度。污染源到出境处水质计算单元250的操作可以参照上面图1描述的步骤S150的操作。污染源到出境处水质计算单元250可以包括到达污染源前水质计算模块、衰减模型确定模块和衰减模型计算模块。其中,到达污染源前水质计算模块用于计算到达污染源之前的混合水质。衰减模型确定模块用于基于菲克定律的一维河流水质模型得到衰减模型。衰减模型计算模块用于根据污染源与其所在的控制区域的出境监测断面之间的距离、以及到达污染源之前的混合水质,利用衰减模型计算污染源的各种污染物到达出境监测断面处的污染物浓度。其中,到达污染源前水质计算模块进一步用于利用零维完全混合模型计算到达污染源之前的混合水质。
量化分析单元260用于将污染源的各种污染物到达出境监测断面处的浓度进行加权平均计算并得到加权平均计算结果,对各个污染源的加权平均计算结果进行排序,从而得到各个污染源排污量化分析结果。量化分析单元260的操作可以参照上面图1描述的步骤S160的操作。
在现有技术的污染源排放贡献率技术中,现有的方案能够计算出某个管控区域内排放贡献最高的污染源。然而,通过该方案计算出的排放贡献最高的污染企业往往是当地规模较大的厂家,这些厂家由于废水处理技术成熟、管控严格,排放的废水在环境评定中已经达标,并不是造成断面水质超标的主要原因。相反,排放不达标的企业却未在贡献率评价的重点关注名单中。由此,这种现有的方案无法对污染管控工作提供支撑。
通过上述实施例计算出的量化分析结果,一些重点关注的超标排污企业都有了合理的得分和排名。比如,某省某磷矿厂由于存在历史遗留的磷矿渣储存不当的问题,导致排水中总磷浓度严重超标。但该渣场平时排水较少,排水多在降雨时发生,因此在现在的贡献率机制中便只占据了较小的贡献度。但上述实施例的量化分析方案中,每当降雨发生,该厂当天的污染浓度排名就会位居前三,表示其是导致监测断面水质超标的主要原因。
图3为根据本发明实施例的用于污染源排污量化分析处理的计算设备的结构框图。
如图3所示,计算设备300可以包括至少一个处理器310、存储器320、内存330、通信接口340以及内部总线350,并且至少一个处理器310、存储器320、内存330和通信接口340经由内部总线350连接在一起。该至少一个处理器310执行在计算机可读存储介质(即,存储器320)中存储或编码的至少一个计算机可读指令(即,上述以软件形式实现的元素)。
在一个实施例中,在存储器320中存储有计算机可执行指令,其当执行时使得至少一个处理器310:基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分;获取目标区域中各个污染源所在的位置,并确定各个污染源所属的控制区域;根据各个监测断面相对于所在河流的上下游关系,确定各个控制区域的入境监测断面和出境监测断面;计算各个污染源与其所在的控制区域的出境监测断面之间的距离;对于各个污染源,根据所述污染源与其所在的控制区域的出境监测断面之间的距离,计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度;将所述污染源的各种污染物到达所述出境监测断面处的浓度进行加权平均计算并得到加权平均计算结果,对各个污染源的加权平均计算结果进行排序,从而得到各个污染源排污量化分析结果。
应该理解的是,在存储器320中存储的计算机可执行指令当执行时使得至少一个处理器310进行本公开的各个实施例中以上结合图1-2描述的各种操作和功能。
在本公开中,计算设备300可以包括但不限于:个人计算机、服务器计算机、工作站、桌面型计算机、膝上型计算机、笔记本计算机、移动计算设备、智能电话、平板计算机、蜂窝电话、个人数字助理(PDA)、手持装置、消息收发设备、可佩戴计算设备、消费电子设备等等。
根据一个实施例,提供了一种例如非暂时性机器可读介质的程序产品。非暂时性机器可读介质可以具有指令(即,上述以软件形式实现的元素),该指令当被机器执行时,使得机器执行本公开的各个实施例中以上结合图1-2描述的各种操作和功能。
具体地,可以提供配有可读存储介质的***或者装置,在该可读存储介质上存储着实现上述实施例中任一实施例的功能的软件程序代码,且使该***或者装置的计算机或处理器读出并执行存储在该可读存储介质中的指令。
在这种情况下,从可读介质读取的程序代码本身可实现上述实施例中任何一项实施例的功能,因此机器可读代码和存储机器可读代码的可读存储介质构成了本发明的一部分。
可读存储介质的实施例包括软盘、硬盘、磁光盘、光盘(如CD-ROM、CD-R、CD-RW、DVD-ROM、DVD-RAM、DVD-RW、DVD-RW)、磁带、非易失性存储卡和ROM。可选择地,可以由通信网络从服务器计算机上或云上下载程序代码。
以上所述仅为本发明的实施例,并非因此限制本发明的权利要求保护范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明权利要求的保护范围内。
Claims (10)
1.一种污染源排污量化分析方法,其特征在于,包括:
基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分;
获取目标区域中各个污染源所在的位置,并确定各个污染源所属的控制区域;
根据各个监测断面相对于所在河流的上下游关系,确定各个控制区域的入境监测断面和出境监测断面;
计算各个污染源与其所在的控制区域的出境监测断面之间的距离;
对于各个污染源,根据所述污染源与其所在的控制区域的出境监测断面之间的距离,计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度;
将所述污染源的各种污染物到达所述出境监测断面处的浓度进行加权平均计算并得到加权平均计算结果,对各个污染源的加权平均计算结果进行排序,从而得到各个污染源排污量化分析结果。
2.根据权利要求1所述的污染源排污量化分析方法,其特征在于,基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分包括:
基于目标区域的自然地理特征、监测断面所在位置和行政区域划分,对目标区域进行控制区域划分。
3.根据权利要求1所述的污染源排污量化分析方法,其特征在于,根据所述污染源与其所在的控制区域的出境监测断面之间的距离,计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度,包括:
计算到达所述污染源之前的混合水质;
基于菲克定律的一维河流水质模型得到衰减模型;
根据所述污染源与其所在的控制区域的出境监测断面之间的距离、以及到达所述污染源之前的混合水质,利用所述衰减模型计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度。
4.根据权利要求3所述的污染源排污量化分析方法,其特征在于,计算到达所述污染源之前的混合水质包括:
利用零维完全混合模型计算到达所述污染源之前的混合水质。
5.一种污染源排污量化分析装置,其特征在于,包括:
控制区域划分单元,用于基于目标区域的自然地理特征和监测断面所在位置,对目标区域进行控制区域划分;
污染源和控制区域关联单元,用于获取目标区域中各个污染源所在的位置,并确定各个污染源所属的控制区域;
监测断面和控制区域关联单元,用于根据各个监测断面相对于所在河流的上下游关系,确定各个控制区域的入境监测断面和出境监测断面;
距离计算单元,用于计算各个污染源与其所在的控制区域的出境监测断面之间的距离;
污染源到出境处水质计算单元,用于对于各个污染源,根据所述污染源与其所在的控制区域的出境监测断面之间的距离,计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度;
量化分析单元,用于将所述污染源的各种污染物到达所述出境监测断面处的浓度进行加权平均计算并得到加权平均计算结果,对各个污染源的加权平均计算结果进行排序,从而得到各个污染源排污量化分析结果。
6.根据权利要求5所述的污染源排污量化分析装置,其特征在于,所述控制区域划分单元进一步用于基于目标区域的自然地理特征、监测断面所在位置和行政区域划分,对目标区域进行控制区域划分。
7.根据权利要求5所述的污染源排污量化分析装置,其特征在于,所述污染源到出境处水质计算单元包括:
到达污染源前水质计算模块,用于计算到达所述污染源之前的混合水质;
衰减模型确定模块,用于基于菲克定律的一维河流水质模型得到衰减模型;
衰减模型计算模块,用于根据所述污染源与其所在的控制区域的出境监测断面之间的距离、以及到达所述污染源之前的混合水质,利用所述衰减模型计算所述污染源的各种污染物到达所述出境监测断面处的污染物浓度。
8.根据权利要求7所述的污染源排污量化分析装置,其特征在于,所述到达污染源前水质计算模块进一步用于利用零维完全混合模型计算到达所述污染源之前的混合水质。
9.一种计算设备,其特征在于,包括:
一个或多个处理器,以及
与所述一个或多个处理器耦合的存储器,所述存储器存储指令,当所述指令被所述一个或多个处理器执行时,使得所述一个或多个处理器执行如权利要求1到4中任一项所述的方法。
10.一种机器可读存储介质,其特征在于,所述机器可读存储介质存储有可执行指令,所述指令当被执行时使得所述机器执行如权利要求1到4中任一项所述的方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011125454 | 2020-10-20 | ||
CN2020111254543 | 2020-10-20 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113139757A true CN113139757A (zh) | 2021-07-20 |
Family
ID=76817971
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110545811.XA Pending CN113139757A (zh) | 2020-10-20 | 2021-05-19 | 污染源排污量化分析方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113139757A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113777256A (zh) * | 2021-08-09 | 2021-12-10 | 力合科技(湖南)股份有限公司 | 环境监测点位的自动化布点方法及***、设备、存储介质 |
CN113806960A (zh) * | 2021-04-30 | 2021-12-17 | 中国环境科学研究院 | 一种工业污染源产排污核算模型及构建方法 |
CN113869798A (zh) * | 2021-11-03 | 2021-12-31 | 大连市勘察测绘研究院集团有限公司 | 一种自动筛选上游风险源的水污染快速溯源方法 |
CN114199314A (zh) * | 2021-12-25 | 2022-03-18 | 江西武大扬帆科技有限公司 | 基于5g和北斗技术的水文监测反馈*** |
CN114705249A (zh) * | 2022-04-11 | 2022-07-05 | 平安国际智慧城市科技股份有限公司 | 基于人工智能的污染物排放量监测方法及相关设备 |
CN115168446A (zh) * | 2022-06-30 | 2022-10-11 | 中科三清科技有限公司 | 污染物溯源方法、装置及电子设备 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103810532A (zh) * | 2012-11-06 | 2014-05-21 | 北京清华同衡规划设计研究院有限公司 | 优化城市排水***运行状况的方法 |
CN105675448A (zh) * | 2014-11-17 | 2016-06-15 | 中国科学院沈阳计算技术研究所有限公司 | 一种面向水环境质量监测的动态模拟污染物扩散的方法 |
CN106557029A (zh) * | 2016-11-11 | 2017-04-05 | 中国科学院生态环境研究中心 | 一种黑臭河流水污染控制与治理的方法 |
CN106600506A (zh) * | 2016-11-25 | 2017-04-26 | 中国科学院生态环境研究中心 | 用于消除湖库型黑臭水体的治理方法 |
CN110765213A (zh) * | 2019-09-07 | 2020-02-07 | 北京化工大学 | 一种地表水流域污染源排放清单(动态清单)的编制方法 |
CN110765419A (zh) * | 2019-10-14 | 2020-02-07 | 中国水利水电科学研究院 | 一种基于多源数据耦合的重要水功能区水质风险测评方法 |
CN111199347A (zh) * | 2019-12-30 | 2020-05-26 | 中国科学院南京地理与湖泊研究所 | 流域污染控制单元分区方法 |
CN111310327A (zh) * | 2020-02-10 | 2020-06-19 | 科盛环保科技股份有限公司 | 一种基于削减量模型的水环境达标定量化设计方法 |
CN111382940A (zh) * | 2020-03-09 | 2020-07-07 | 生态环境部环境规划院 | 面源污染风险的确定方法、装置、设备及存储介质 |
-
2021
- 2021-05-19 CN CN202110545811.XA patent/CN113139757A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103810532A (zh) * | 2012-11-06 | 2014-05-21 | 北京清华同衡规划设计研究院有限公司 | 优化城市排水***运行状况的方法 |
CN105675448A (zh) * | 2014-11-17 | 2016-06-15 | 中国科学院沈阳计算技术研究所有限公司 | 一种面向水环境质量监测的动态模拟污染物扩散的方法 |
CN106557029A (zh) * | 2016-11-11 | 2017-04-05 | 中国科学院生态环境研究中心 | 一种黑臭河流水污染控制与治理的方法 |
CN106600506A (zh) * | 2016-11-25 | 2017-04-26 | 中国科学院生态环境研究中心 | 用于消除湖库型黑臭水体的治理方法 |
CN110765213A (zh) * | 2019-09-07 | 2020-02-07 | 北京化工大学 | 一种地表水流域污染源排放清单(动态清单)的编制方法 |
CN110765419A (zh) * | 2019-10-14 | 2020-02-07 | 中国水利水电科学研究院 | 一种基于多源数据耦合的重要水功能区水质风险测评方法 |
CN111199347A (zh) * | 2019-12-30 | 2020-05-26 | 中国科学院南京地理与湖泊研究所 | 流域污染控制单元分区方法 |
CN111310327A (zh) * | 2020-02-10 | 2020-06-19 | 科盛环保科技股份有限公司 | 一种基于削减量模型的水环境达标定量化设计方法 |
CN111382940A (zh) * | 2020-03-09 | 2020-07-07 | 生态环境部环境规划院 | 面源污染风险的确定方法、装置、设备及存储介质 |
Non-Patent Citations (4)
Title |
---|
唐晓青,周旌主编: "《环境质量监测点位布设技术指南》", 30 June 2016, 河北科学技术出版社 * |
孙辰: "汉江襄阳段水环境容量及总量控制研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑》 * |
徐景航,傅国伟编: "《环境***工程》", 31 October 1990, 中国环境科学出版社 * |
生态环境部环境工程评估中心: "《建设项目环境影响评价全过程管理及高新技术研究与实践》", 31 October 2019, 中国环境出版集团 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113806960A (zh) * | 2021-04-30 | 2021-12-17 | 中国环境科学研究院 | 一种工业污染源产排污核算模型及构建方法 |
CN113806960B (zh) * | 2021-04-30 | 2023-02-21 | 中国环境科学研究院 | 一种工业污染源产排污核算模型的构建方法 |
CN113777256A (zh) * | 2021-08-09 | 2021-12-10 | 力合科技(湖南)股份有限公司 | 环境监测点位的自动化布点方法及***、设备、存储介质 |
CN113869798A (zh) * | 2021-11-03 | 2021-12-31 | 大连市勘察测绘研究院集团有限公司 | 一种自动筛选上游风险源的水污染快速溯源方法 |
CN114199314A (zh) * | 2021-12-25 | 2022-03-18 | 江西武大扬帆科技有限公司 | 基于5g和北斗技术的水文监测反馈*** |
CN114199314B (zh) * | 2021-12-25 | 2023-05-16 | 中铁水利信息科技有限公司 | 基于5g和北斗技术的水文监测反馈*** |
CN114705249A (zh) * | 2022-04-11 | 2022-07-05 | 平安国际智慧城市科技股份有限公司 | 基于人工智能的污染物排放量监测方法及相关设备 |
CN114705249B (zh) * | 2022-04-11 | 2024-04-30 | 平安国际智慧城市科技股份有限公司 | 基于人工智能的污染物排放量监测方法及相关设备 |
CN115168446A (zh) * | 2022-06-30 | 2022-10-11 | 中科三清科技有限公司 | 污染物溯源方法、装置及电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113139757A (zh) | 污染源排污量化分析方法及装置 | |
Kang et al. | Multicriteria Evaluation of Water Resources Sustainability in the Context of Watershed Management 1 | |
CN104392100A (zh) | 基于水质在线监测***的污染源扩散预警方法 | |
CN112001610A (zh) | 农业面源污染的处理方法及装置 | |
Varekar et al. | Design of sampling locations for river water quality monitoring considering seasonal variation of point and diffuse pollution loads | |
CN111444595B (zh) | 基于swmm与efdc耦合模型的截污工程环境效应评估方法及装置 | |
Zhang et al. | Determination of the distribution of flood forecasting error | |
Liang et al. | A pollutant load hierarchical allocation method integrated in an environmental capacity management system for Zhushan Bay, Taihu Lake | |
CN113158385B (zh) | 一种农业面源污染关键源区识别方法及*** | |
Yasarer et al. | Characterizing ponds in a watershed simulation and evaluating their influence on streamflow in a Mississippi watershed | |
Wang et al. | Accurately early warning to water quality pollutant risk by mobile model system with optimization technology | |
CN113139756A (zh) | 污染源排污敏感性评估方法及装置 | |
Tu et al. | Bivariate design of hydrological droughts and their alterations under a changing environment | |
Xiao et al. | Water quality assessment and pollution source apportionment using multivariate statistical techniques: a case study of the Laixi River Basin, China | |
Gomes et al. | Regionalization of precipitation with determination of homogeneous regions via fuzzy c-means | |
Habete et al. | Impact of forecasted land use change on design peak discharge at watershed and catchment scales: simple equation to predict changes | |
Zhou et al. | Conceptual planning approach of low impact development s for combined water quality‐quantity control at an urban scale: A case study in Southern China | |
Azadi et al. | Predicting sediment load using stochastic model and rating curves in a hydrological station | |
Asadi et al. | An integrated approach for prioritization of river water quality sampling points using modified Sanders, analytic network process, and hydrodynamic modeling | |
Osmi et al. | Development of Total Maximum Daily Load using water quality modelling as an approach for watershed management in Malaysia | |
CN112287613A (zh) | 针对流域水环境控制断面的污染物削减方法 | |
Chen et al. | Development and application of a multi-centre cloud platform architecture for water environment management | |
CN111353718A (zh) | 基于swmm与efdc的湿地、补水工程环境效应评估方法及装置 | |
Wadeea et al. | Effects of industrial wastewater on water quality Tigris river at Baghdad using (GIS) technique | |
Wang et al. | Application of WASP model and Gini coefficient in total mass control of water pollutants: a case study in Xicheng Canal, China |
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 |