CN114661849B - 污染溯源方法及装置 - Google Patents
污染溯源方法及装置 Download PDFInfo
- Publication number
- CN114661849B CN114661849B CN202210288649.2A CN202210288649A CN114661849B CN 114661849 B CN114661849 B CN 114661849B CN 202210288649 A CN202210288649 A CN 202210288649A CN 114661849 B CN114661849 B CN 114661849B
- Authority
- CN
- China
- Prior art keywords
- pollution
- model
- track
- diffusion
- tracing
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000009792 diffusion process Methods 0.000 claims abstract description 73
- 238000012216 screening Methods 0.000 claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 239000003344 environmental pollutant Substances 0.000 claims description 16
- 231100000719 pollutant Toxicity 0.000 claims description 16
- 238000011109 contamination Methods 0.000 claims description 5
- 239000002245 particle Substances 0.000 claims description 5
- 238000012163 sequencing technique Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 3
- 230000002265 prevention Effects 0.000 abstract description 3
- 238000011835 investigation Methods 0.000 description 5
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000013618 particulate matter Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 2
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 2
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 229910002091 carbon monoxide Inorganic materials 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000779 smoke Substances 0.000 description 2
- 239000012855 volatile organic compound Substances 0.000 description 2
- 239000002028 Biomass Substances 0.000 description 1
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 1
- 239000000809 air pollutant Substances 0.000 description 1
- 231100001243 air pollutant Toxicity 0.000 description 1
- 229910021529 ammonia Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
- Y02P90/845—Inventory and reporting systems for greenhouse gases [GHG]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Databases & Information Systems (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Economics (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Fuzzy Systems (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种污染溯源方法及装置。该方法包括:将空气质量预报数据输入气团轨迹追溯模型,计算得到预设时刻所选点位经过预设时长在不同高度的后向扩散轨迹;提取后向扩散轨迹上各高度的经纬度数据组,根据各高度以及对应的经纬度数据组绘制轨迹线;根据来自污染源排放清单和/或污染源普查得到的各污染源的排放高度和位置信息,筛选出与所绘制的轨迹线间的距离小于预设长度的污染源,得到轨迹周边的污染源排放信息。该装置包括扩散轨迹计算单元、轨迹线绘制单元、以及污染源筛选单元。这种方法及装置结合后向轨迹法和污染源信息法各自的优势,将结果联系起来,进一步精细化科学性地聚焦了污染溯源结果,为后续污染防控工作提供了指导和方向。
Description
技术领域
本发明涉及大气污染溯源和管控技术领域,更为具体来说,本发明涉及一种污染溯源方法及装置。
背景技术
近年来,我国空气质量明显改善,随着污染物浓度的下降,污染管控工作逐步向精细化、差异化发展,对于污染来源的追溯方式和准确率有了更高要求。
当前,常见的污染来源追溯方法主要有:污染源排放清单法、数值模型模拟法、受体模型解析法等,其中,数值模型模拟法以二代大气扩散模型和三代空气质量数值模型为常用方式,另有忽略污染传输中物理化学转化过程而简化应用的气团扩散模型法。
大气扩散模型是大气污染治理的重要工具之一,其目的是研究排放进入大气的空气污染物的扩散、转化、迁移和清除的规律,目前描述大气扩散规律的模式主要分为3类,包括:经修正的高斯烟流扩散模型、欧拉扩散模型和拉格朗日扩散模型。
尽管现有的污染来源追溯方法多样,但各种方法各占一席,且解析结果对于污染来源的地区、行业等指示不尽相同。目前,缺少多类型方法结合使用来追溯污染来源的实践,因此不能很好地结合各种方法的优势,精细化溯源的发展也未能快速前进。
发明内容
本发明创新地提供了一种污染溯源方法及装置,将后向轨迹结果与污染源信息联用,实现更精细化更准确的污染溯源。
为实现上述的技术目的,一方面,本发明公开了一种污染溯源方法,所述污染溯源方法包括:将空气质量预报数据输入气团轨迹追溯模型,计算得到预设时刻的所选点位经过预设时长在不同高度的后向扩散轨迹;提取所述后向扩散轨迹上各个高度的经度和纬度数据组,根据所述各个高度以及对应的所述经度和纬度数据组绘制轨迹线;根据来自污染源排放清单和/或污染源普查得到的各个污染源的排放高度和位置信息,筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源,从而得到所述轨迹周边的污染源排放信息。
进一步地,所述的污染溯源方法,在筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源之后,还包括:将筛选出的所述污染源按照排放量从大到小进行排序。
进一步地,所述的污染溯源方法,在筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源之后,还包括:将筛选出不同地区、不同行业的所述污染源按照排放量从大到小进行排序以及根据所述空气质量预报数据的初始条件和边界条件对所述污染源的污染物生成的贡献从大到小进行排序。
进一步地,对于所述污染溯源方法,所述空气质量预报数据包括天气预报模式的输出结果和/或全球预报系统的预报数据。
进一步地,对于所述的污染溯源方法,所述气团轨迹追溯模型可以包括高斯烟流扩散模型、欧拉扩散模型、拉格朗日扩散模型、或者以高斯烟流扩散模型、欧拉扩散模型和/或拉格朗日扩散模型为原理的衍生模型。
进一步地,对于所述污染溯源方法,所述拉格朗日扩散模型的衍生模型包括混合粒子拉格朗日积分轨迹模型。
进一步地,对于所述污染溯源方法,所述气团轨迹追溯模型还包括烟羽动力模型或者以烟羽动力模型为原理的衍生模型。
实现上述的技术目的,另一方面,本发明公开了一种污染溯源装置,所述污染溯源装置包括:扩散轨迹计算单元,用于将空气质量预报数据输入气团轨迹追溯模型,计算得到预设时刻的所选点位经过预设时长在不同高度的后向扩散轨迹;轨迹线绘制单元,用于提取所述后向扩散轨迹上各个高度的经度和纬度数据组,根据所述各个高度以及对应的所述经度和纬度数据组绘制轨迹线;污染源筛选单元,用于根据来自污染源排放清单和/或污染源普查得到的各个污染源的排放高度和位置信息,筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源,从而得到所述轨迹周边的污染源排放信息。
进一步地,所述污染溯源装置还包括:排序单元,用于将筛选出的所述污染源按照排放量从大到小进行排序。
进一步地,对于所述污染溯源装置,所述气团轨迹追溯模型可以包括高斯烟流扩散模型、欧拉扩散模型、拉格朗日扩散模型、或者以高斯烟流扩散模型、欧拉扩散模型和/或拉格朗日扩散模型为原理的衍生模型。
为实现上述的技术目的,又一方面,本发明公开了一种计算设备。所述计算设备包括:一个或多个处理器,以及与所述一个或多个处理器耦合的存储器,所述存储器存储指令,当所述指令被所述一个或多个处理器执行时,使得所述一个或多个处理器执行上述方法。
为实现上述的技术目的,再一方面,本发明公开了一种机器可读存储介质。所述机器可读存储介质存储有可执行指令,所述指令当被执行时使得所述机器执行上述方法。
本发明的有益效果为:
现有常用的污染溯源方法中,Hysplit后向轨迹法仅能给出污染气团的来源路径,管控措施只能通过这一结果确定大致管控方位,但并未做到精准管控的目的。而通过污染源排放清单和/或污染源普查数据得到的信息,仅能整理得到整个区域的重点大排量污染源结果,但并非这些污染源均对某次污染产生了影响,因此仅通过污染源排放清单和/或污染源普查数据也未能做到科学治污的目的。本发明实施例提供的污染溯源方法及装置结合后向轨迹法和污染源信息法各自的优势,将结果联系起来,进一步精细化、科学性地聚焦了污染溯源结果,进而为后续的污染防控工作提供了指导和方向。
附图说明
图中,
图1为本发明一个实施例提供的污染溯源方法的流程图;
图2为本发明一个示例提供的污染溯源方法的流程图;
图3为本发明另一个实施例提供的污染溯源装置的结构示意图;
图4为根据本发明实施例的用于污染溯源处理的计算设备的结构框图。
具体实施方式
下面结合说明书附图对本发明提供的污染溯源方法及装置进行详细的解释和说明。
图1为本发明一个实施例提供的污染溯源方法的流程图。图2为本发明一个示例提供的污染溯源方法的流程图。
如图1和图2所示,在步骤S110,将空气质量预报数据输入气团轨迹追溯模型,计算得到预设时刻的所选点位经过预设时长在不同高度的后向扩散轨迹。其中,所选点位的气团所在的高度和预设时长可以在模型输入参数中任意配置。
空气质量预报数据可以包括天气预报模式(WRF,The Weather Research andForecasting Model)的输出结果和/或全球预报系统(GFS,Global Forecast System)的预报数据。
气团轨迹追溯模型可以包括高斯烟流扩散模型、欧拉扩散模型、拉格朗日扩散模型、或者以高斯烟流扩散模型、欧拉扩散模型和/或拉格朗日扩散模型为原理的各类衍生模型。输入数据可以随着模型的选用而有所差异。其中,拉格朗日扩散模型的衍生模型可以包括混合粒子拉格朗日积分轨迹模型(Hysplit,Hybrid Single-Particle LagrangianIntegrated Trajectory Model)。HYSPLIT模型又称为后向轨迹模型,不仅可以计算污染物浓度扩散的沉积,也能计算单个污染物粒子的运动轨迹,是当前大气科学界使用最广泛的大气传输和扩散模型之一。可以利用HYSPLIT模型计算得到具体某点不同高度气团的来源与去向轨迹结果,对污染溯源具有相当的指导作用。
作为一种可选实施方式,上述气团轨迹追溯模型还包括烟羽动力模型或者以烟羽动力模型为原理的衍生模型。
在步骤S120,提取后向扩散轨迹上各个高度的经度和纬度数据组,根据各个高度以及对应的经度和纬度数据组绘制轨迹线。
在步骤S130,根据来自污染源排放清单和/或污染源普查得到的各个污染源的排放高度和位置信息,筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源,从而得到轨迹周边的污染源排放信息。其中,所述污染源与所绘制的轨迹线之间的距离,为该污染源与所绘制的轨迹线之间的最小距离。可以通过整理污染源排放清单和/或污染源普查得到的各类污染源数据,将信息按照一定格式存入数据库中。其中,预设长度可以自定义设置。
污染源排放清单是对某一地区一种或几种污染物排放源的排放量进行估算,一套完整的大气污染物排放清单应当覆盖化石燃料固定燃烧、工艺过程、移动源、溶剂使用、开放扬尘、生物质燃烧和农业等排放源,包含二氧化硫(SO2)、氮氧化物(NOx)、一氧化碳(CO)、挥发性有机物(VOCs)、氨(NH3)、一次颗粒物(PM2.5和PM10)和臭氧(O3)等大气污染物,并具备动态更新机制。编制污染源排放清单需要进行大量的前期调研工作,包括所编制区域内各类污染源的基础信息、生产活动方式/环境以及生产活动产量等数据,从而计算得到各污染源对不同污染物质的排放量,这些调研可以在污染源普查结果的基础上进一步展开。结合编制排放清单时对点源、线源和面源的基础调查资料,可以建立污染源排放基础数据库。
作为一种可选实施方式,在筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源之后,该实施例的污染溯源方法还可以包括步骤:将筛选出的污染源按照排放量从大到小进行排序。将这些污染源按排放量进行排序之后,可以根据排放量筛选得到应重点管控的排放源,比如某企业、某工地等信息,从而来达到精准管控、科学治污的目的。
作为另一种可选实施方式,在筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源之后,该实施例的污染溯源方法还可以包括如下步骤:将筛选出不同地区、不同行业的污染源按照排放量从大到小进行排序以及根据空气质量预报数据的初始条件和边界条件对污染源的污染物生成的贡献从大到小进行排序。由于输入气团轨迹追溯模型的空气质量预报数据存在一定的区域范围,而该区域范围之外的区域的相关空气质量数据(如污染物浓度、风向等)对该区域范围的空气质量也存在一定的贡献,因此,可以将该区域范围对应的空气质量预报数据作为初始条件,该区域范围的边界对应的相关数据作为边界条件。进而,在上述将筛选出的污染源按照排放量从大到小进行排序之外,还可以计算初始条件和边界条件对上述筛选出的污染源的的污染物生成的贡献,并从大到小进行排序。从而,通过排放量和贡献大小的排序,筛选得到应重点管控的排放源,进一步提高筛选的精细化程度,提高得到的排放源的精确性。
图3为本发明另一个实施例提供的污染溯源装置的结构示意图。如图3所示,该实施例提供的污染溯源装置300包括扩散轨迹计算单元310、轨迹线绘制单元320、以及污染源筛选单元330。
扩散轨迹计算单元310用于将空气质量预报数据输入气团轨迹追溯模型,计算得到预设时刻的所选点位经过预设时长在不同高度的后向扩散轨迹。其中,空气质量预报数据可以包括WRF的输出结果和/或GFS的预报数据。气团轨迹追溯模型可以包括高斯烟流扩散模型、欧拉扩散模型、拉格朗日扩散模型、或者以高斯烟流扩散模型、欧拉扩散模型和/或拉格朗日扩散模型为原理的各类衍生模型。输入数据可以随着模型的选用而有所差异。其中,拉格朗日扩散模型的衍生模型可以包括混合粒子拉格朗日积分轨迹模型(Hysplit,Hybrid Single-Particle Lagrangian Integrated Trajectory Model)。作为一种可选实施方式,气团轨迹追溯模型还包括烟羽动力模型或者以烟羽动力模型为原理的衍生模型。扩散轨迹计算单元310的操作可以参照上面图1描述的步骤S110的操作。
轨迹线绘制单元320用于提取后向扩散轨迹上各个高度的经度和纬度数据组,根据各个高度以及对应的经度和纬度数据组绘制轨迹线。扩散轨迹线绘制单元320的操作可以参照上面图1描述的步骤S120的操作。
污染源筛选单元330用于根据来自污染源排放清单和/或污染源普查得到的各个污染源的排放高度和位置信息,筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源,从而得到轨迹周边的污染源排放信息。污染源筛选单元330的操作可以参照上面图1描述的步骤S130的操作。
作为一种可选实施方式,该实施例的污染溯源装置300还可以包括排序单元,用于将筛选出的污染源按照排放量从大到小进行排序。将这些污染源按排放量进行排序之后,可以根据排放量筛选得到应重点管控的排放源,比如某企业、某工地等信息,从而来达到精准管控、科学治污的目的。
作为另一种可选实施方式,排序单元用于将筛选出不同地区、不同行业的所述污染源按照排放量从大到小进行排序以及根据所述空气质量预报数据的初始条件和边界条件对所述污染源的污染物生成的贡献从大到小进行排序。
现有常用的污染溯源方法中,Hysplit后向轨迹法仅能给出污染气团的来源路径,管控措施只能通过这一结果确定大致管控方位,但并未做到精准管控的目的。而通过污染源排放清单和/或污染源普查数据得到的信息,仅能整理得到整个区域的重点大排量污染源结果,但并非这些污染源均对某次污染产生了影响,因此仅通过污染源排放清单和/或污染源普查数据也未能做到科学治污的目的。本发明实施例提供的污染溯源方法及装置结合后向轨迹法和污染源信息法各自的优势,将结果联系起来,进一步精细化、科学性地聚焦了污染溯源结果,进而为后续的污染防控工作提供了指导和方向。
图4为根据本发明实施例的用于污染溯源处理的计算设备的结构框图。
如图4所示,计算设备400可以包括至少一个处理器410、存储器420、内存430、通信接口440以及内部总线450,并且至少一个处理器410、存储器420、内存430和通信接口440经由总线450连接在一起。该至少一个处理器410执行在计算机可读存储介质(即,存储器420)中存储或编码的至少一个计算机可读指令(即,上述以软件形式实现的元素)。
在一个实施例中,在存储器420中存储有计算机可执行指令,当其执行时使得至少一个处理器410执行:将空气质量预报数据输入气团轨迹追溯模型,计算得到预设时刻的所选点位经过预设时长在不同高度的后向扩散轨迹;提取所述后向扩散轨迹上各个高度的经度和纬度数据组,根据所述各个高度以及对应的所述经度和纬度数据组绘制轨迹线;根据来自污染源排放清单和/或污染源普查得到的各个污染源的排放高度和位置信息,筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源,从而得到轨迹周边的污染源排放信息。
应该理解的是,在存储器420中存储的计算机可执行指令当执行时使得至少一个处理器410进行本公开的各个实施例中以上结合图1-3描述的各种操作和功能。
在本公开中,计算设备400可以包括但不限于:个人计算机、服务器计算机、工作站、桌面型计算机、膝上型计算机、笔记本计算机、移动计算设备、智能电话、平板计算机、蜂窝电话、个人数字助理(PDA)、手持装置、消息收发设备、可佩戴计算设备、消费电子设备等等。
根据一个实施例,提供了一种例如非暂时性机器可读介质的程序产品。非暂时性机器可读介质可以具有指令(即,上述以软件形式实现的元素),该指令当被机器执行时,使得机器执行本公开的各个实施例中以上结合图1-3描述的各种操作和功能。
具体地,可以提供配有可读存储介质的系统或者装置,在该可读存储介质上存储着实现上述实施例中任一实施例的功能的软件程序代码,且使该系统或者装置的计算机或处理器读出并执行存储在该可读存储介质中的指令。
在这种情况下,从可读介质读取的程序代码本身可实现上述实施例中任何一项实施例的功能,因此机器可读代码和存储机器可读代码的可读存储介质构成了本发明的一部分。
可读存储介质的实施例包括软盘、硬盘、磁光盘、光盘(如CD-ROM、CD-R、CD-RW、DVD-ROM、DVD-RAM、DVD-RW、DVD-RW)、磁带、非易失性存储卡和ROM。可选择地,可以由通信网络从服务器计算机上或云上下载程序代码。
以上所述仅为本发明的实施例,并非因此限制本发明的权利要求保护范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明权利要求的保护范围内。
Claims (10)
1.一种污染溯源方法,其特征在于,包括:
将区域范围内的天气预报模式的输出结果和/或全球预报系统的预报数据输入气团轨迹追溯模型,计算得到预设时刻的所选点位经过预设时长在不同高度的后向扩散轨迹;
提取所述后向扩散轨迹上各个高度的经度和纬度数据组,根据所述各个高度以及对应的所述经度和纬度数据组绘制轨迹线;
根据来自污染源排放清单和/或污染源普查得到的各个污染源的排放高度和位置信息,筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源;
通过排放量和贡献大小的排序,筛选得到应重点管控的排放源;
其中,排放量的排序包括:将筛选出不同地区、不同行业的所述污染源按照排放量从大到小进行排序;
其中,贡献大小的排序,包括:将区域范围内的天气预报模式的输出结果和/或全球预报系统的预报数据作为初始条件,将该区域范围的边界对应的相关数据作为边界条件,其中,所述相关数据是指对该区域范围的空气质量存在一定的贡献的数据,包括污染物浓度和风向;根据所述初始条件和边界条件对所述污染源的污染物生成的贡献从大到小进行排序。
2.根据权利要求1所述的污染溯源方法,其特征在于,在筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源之后,还包括:
将筛选出的所述污染源按照排放量从大到小进行排序。
3.根据权利要求1-2中任一项所述的污染溯源方法,其特征在于,所述气团轨迹追溯模型包括高斯烟流扩散模型、欧拉扩散模型、拉格朗日扩散模型、或者以高斯烟流扩散模型、欧拉扩散模型和/或拉格朗日扩散模型为原理的衍生模型。
4.根据权利要求3所述的污染溯源方法,其特征在于,所述拉格朗日扩散模型的衍生模型包括混合粒子拉格朗日积分轨迹模型。
5.根据权利要求3所述的污染溯源方法,其特征在于,所述气团轨迹追溯模型还包括烟羽动力模型或者以烟羽动力模型为原理的衍生模型。
6.一种污染溯源装置,其特征在于,包括:
扩散轨迹计算单元,用于将区域范围内的天气预报模式的输出结果和/或全球预报系统的预报数据输入气团轨迹追溯模型,计算得到预设时刻的所选点位经过预设时长在不同高度的后向扩散轨迹;
轨迹线绘制单元,用于提取所述后向扩散轨迹上各个高度的经度和纬度数据组,根据所述各个高度以及对应的所述经度和纬度数据组绘制轨迹线;
污染源筛选单元,用于:
根据来自污染源排放清单和/或污染源普查得到的各个污染源的排放高度和位置信息,筛选出与所绘制的轨迹线之间的距离小于预设长度的污染源;
通过排放量和贡献大小的排序,筛选得到应重点管控的排放源;
其中,排放量的排序包括:将筛选出不同地区、不同行业的所述污染源按照排放量从大到小进行排序;
其中,贡献大小的排序,包括:将区域范围内的天气预报模式的输出结果和/或全球预报系统的预报数据作为初始条件,将该区域范围的边界对应的相关数据作为边界条件,其中,所述相关数据是指对该区域范围的空气质量存在一定的贡献的数据,包括污染物浓度和风向;根据所述初始条件和边界条件对所述污染源的污染物生成的贡献从大到小进行排序。
7.根据权利要求6所述的污染溯源装置,其特征在于,还包括:排序单元,用于将筛选出的所述污染源按照排放量从大到小进行排序。
8.根据权利要求6或7所述的污染溯源装置,其特征在于,所述气团轨迹追溯模型可以包括高斯烟流扩散模型、欧拉扩散模型、拉格朗日扩散模型、或者以高斯烟流扩散模型、欧拉扩散模型和/或拉格朗日扩散模型为原理的衍生模型。
9.一种计算设备,其特征在于,包括:
一个或多个处理器,以及
与所述一个或多个处理器耦合的存储器,所述存储器存储指令,当所述指令被所述一个或多个处理器执行时,使得所述一个或多个处理器执行如权利要求1到5中任一项所述的方法。
10.一种机器可读存储介质,其特征在于,所述机器可读存储介质存储有可执行指令,所述指令当被执行时使得所述机器执行如权利要求1到5中任一项所述的方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110362245.9A CN112988940A (zh) | 2021-04-02 | 2021-04-02 | 污染溯源方法及装置 |
CN2021103622459 | 2021-04-02 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114661849A CN114661849A (zh) | 2022-06-24 |
CN114661849B true CN114661849B (zh) | 2023-07-14 |
Family
ID=76339056
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110362245.9A Pending CN112988940A (zh) | 2021-04-02 | 2021-04-02 | 污染溯源方法及装置 |
CN202210288649.2A Active CN114661849B (zh) | 2021-04-02 | 2022-03-23 | 污染溯源方法及装置 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110362245.9A Pending CN112988940A (zh) | 2021-04-02 | 2021-04-02 | 污染溯源方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN112988940A (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113763509B (zh) * | 2021-09-08 | 2022-09-02 | 中科三清科技有限公司 | 轨迹图绘制方法和装置 |
CN114155129B (zh) * | 2021-11-29 | 2023-08-29 | 中科三清科技有限公司 | 一种基于工业园区的大气环境溯评方法和系统 |
CN114359002B (zh) * | 2022-03-21 | 2022-05-20 | 四川国蓝中天环境科技集团有限公司 | 基于移动监测趋势挖掘的大气污染小尺度溯源方法及系统 |
CN116109323B (zh) * | 2022-12-26 | 2023-08-01 | 北京中科三清环境技术有限公司 | 臭氧污染溯源方法、装置、设备及存储介质 |
CN115792137B (zh) * | 2023-01-17 | 2023-06-06 | 河北先河环保科技股份有限公司 | 大气污染溯源方法及装置、终端 |
CN117970527A (zh) * | 2024-04-02 | 2024-05-03 | 南昌云宜然科技有限公司 | 一种大气污染物的组网溯源监测方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006105805A (ja) * | 2004-10-06 | 2006-04-20 | Matsushita Electric Ind Co Ltd | 軌跡データ管理装置およびプログラム |
CN106295905A (zh) * | 2016-08-22 | 2017-01-04 | 南京大学 | 一种基于拉格朗日输送模型的空气质量快速溯源预报方法 |
KR101948546B1 (ko) * | 2017-11-30 | 2019-05-20 | 대한민국(환경부 국립환경과학원장) | 대기오염물질 배출영향분석 시스템 및 그 방법 |
CN110807725A (zh) * | 2019-11-06 | 2020-02-18 | 软通动力信息技术有限公司 | 一种大气污染溯源方法、装置、计算设备和介质 |
CN110824110A (zh) * | 2019-10-30 | 2020-02-21 | 山东大学 | 一种基于拉格朗日轨迹模式和化学盒子模式的区域臭氧污染溯源系统 |
CN111157688A (zh) * | 2020-03-06 | 2020-05-15 | 北京市环境保护监测中心 | 评价污染源对空气质量监测站影响的方法和装置 |
CN111506574A (zh) * | 2020-03-19 | 2020-08-07 | 平安国际智慧城市科技股份有限公司 | 基于r树的污染物溯源方法、装置及其相关设备 |
CN112131739A (zh) * | 2020-09-22 | 2020-12-25 | 南京创蓝科技有限公司 | 一种乡镇级别大气污染溯源的预报方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11226323B2 (en) * | 2018-04-27 | 2022-01-18 | International Business Machines Corporation | Air-pollution emission source monitoring |
WO2021026792A1 (zh) * | 2019-08-14 | 2021-02-18 | 柯灵爱尔(北京)环境技术中心 | 一种选择污染物治理措施的方法 |
CN112182064B (zh) * | 2020-09-25 | 2021-07-20 | 中科三清科技有限公司 | 一种污染物来源分析方法、装置、电子设备及存储介质 |
-
2021
- 2021-04-02 CN CN202110362245.9A patent/CN112988940A/zh active Pending
-
2022
- 2022-03-23 CN CN202210288649.2A patent/CN114661849B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006105805A (ja) * | 2004-10-06 | 2006-04-20 | Matsushita Electric Ind Co Ltd | 軌跡データ管理装置およびプログラム |
CN106295905A (zh) * | 2016-08-22 | 2017-01-04 | 南京大学 | 一种基于拉格朗日输送模型的空气质量快速溯源预报方法 |
KR101948546B1 (ko) * | 2017-11-30 | 2019-05-20 | 대한민국(환경부 국립환경과학원장) | 대기오염물질 배출영향분석 시스템 및 그 방법 |
CN110824110A (zh) * | 2019-10-30 | 2020-02-21 | 山东大学 | 一种基于拉格朗日轨迹模式和化学盒子模式的区域臭氧污染溯源系统 |
CN110807725A (zh) * | 2019-11-06 | 2020-02-18 | 软通动力信息技术有限公司 | 一种大气污染溯源方法、装置、计算设备和介质 |
CN111157688A (zh) * | 2020-03-06 | 2020-05-15 | 北京市环境保护监测中心 | 评价污染源对空气质量监测站影响的方法和装置 |
CN111506574A (zh) * | 2020-03-19 | 2020-08-07 | 平安国际智慧城市科技股份有限公司 | 基于r树的污染物溯源方法、装置及其相关设备 |
CN112131739A (zh) * | 2020-09-22 | 2020-12-25 | 南京创蓝科技有限公司 | 一种乡镇级别大气污染溯源的预报方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114661849A (zh) | 2022-06-24 |
CN112988940A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114661849B (zh) | 污染溯源方法及装置 | |
CN110824110B (zh) | 一种基于拉格朗日轨迹模式和化学盒子模式的区域臭氧污染溯源系统 | |
Carslaw et al. | Openair—an R package for air quality data analysis | |
KR102324392B1 (ko) | 지역별 미세먼지 농도 예측방법 및 시스템 | |
CN114707831A (zh) | 臭氧污染防治的动态决策方法及装置 | |
CN108489875B (zh) | 一种基于时段统计分析的污染物溯源系统及方法 | |
CN111753426B (zh) | 颗粒物污染的源解析方法和装置 | |
CN110348746B (zh) | 基于单个污染源的空气质量影响评估方法和装置 | |
CN115271547B (zh) | 臭氧污染的源解析方法、装置和电子设备 | |
CN109376443A (zh) | 区域大气污染环境仿真动态模拟系统 | |
CN109115949B (zh) | 基于大数据的污染溯源方法及计算机可读介质 | |
CN107767081B (zh) | 一种燃煤电厂大气环境影响评价方法及装置 | |
CN114971031A (zh) | 大气污染排放的溯源方法、装置、电子设备及存储介质 | |
CN117092297A (zh) | 一种工业园区大气污染物的溯源方法、系统、设备及介质 | |
Havrylenko et al. | Criteria for optimising air quality monitoring in Ukrainian cities (by example of Kyiv) | |
CN116485191A (zh) | 基于拉格朗日扩散模型的环境影响预警方法及系统 | |
KR20170079552A (ko) | 환경감시 업무지원 시스템 및 그를 이용한 환경감시 방법 | |
Xu et al. | How to reach haze control targets by air pollutants emission reduction in the Beijing-Tianjin-Hebei region of China? | |
CN107766429A (zh) | 一种环保大数据自适应精准分类处理方法及系统 | |
CN114896952B (zh) | 一种大气污染源排放清单编制方法及装置 | |
Henderson et al. | Comparison of Lagrangian Process Analysis tools for Eulerian air quality models | |
JP6672943B2 (ja) | 物質発生源推定方法及び物質発生源推定装置 | |
CN113672873A (zh) | 污染物源解析方法、装置、电子设备及存储介质 | |
CN115705510A (zh) | 厂区气体污染分布预测方法、系统、电子设备及存储介质 | |
Tang et al. | A new circulating accumulation emission model for assessing dust emission from open pit mine |
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 |