CN113284244A - 一种基于格点化扩散模拟的大气污染物溯源方法 - Google Patents

一种基于格点化扩散模拟的大气污染物溯源方法 Download PDF

Info

Publication number
CN113284244A
CN113284244A CN202110426596.1A CN202110426596A CN113284244A CN 113284244 A CN113284244 A CN 113284244A CN 202110426596 A CN202110426596 A CN 202110426596A CN 113284244 A CN113284244 A CN 113284244A
Authority
CN
China
Prior art keywords
pollutant
grid
emission
target area
method based
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.)
Granted
Application number
CN202110426596.1A
Other languages
English (en)
Other versions
CN113284244B (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.)
Sichuan University
Original Assignee
Sichuan 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 Sichuan University filed Critical Sichuan University
Priority to CN202110426596.1A priority Critical patent/CN113284244B/zh
Publication of CN113284244A publication Critical patent/CN113284244A/zh
Application granted granted Critical
Publication of CN113284244B publication Critical patent/CN113284244B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205Re-meshing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29Geographical information databases
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Databases & Information Systems (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Software Systems (AREA)
  • Game Theory and Decision Science (AREA)
  • Geometry (AREA)
  • Remote Sensing (AREA)
  • General Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Computer Graphics (AREA)
  • Data Mining & Analysis (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

一种基于格点化扩散模拟的大气污染物溯源方法
技术领域
本发明涉及大气环境技术领域,具体涉及一种基于格点化扩散模拟的大气污染物溯源方法。
背景技术
扩散模拟是一种常用的大气污染物溯源方法,现有技术主要包括两个方面:受体点后向轨迹方法,可以解析出影响受体点的气团来源,从而实现污染物的溯源。该方法仅能识别气团路径,在空间上缺乏精细解析能力。基于后向轨迹衍生的PSCF(potential sourcecontribution function)方法,通过对一定时间内的多条后向轨迹进行统计分析,解析出这段时间内影响受体点的气团来源概率分布,从而实现污染物的溯源。该方法针对一段时间进行解析,在时间上缺乏精细解析能力。排放源前向扩散方法,在重点关注的排放源位置模拟排放,通过扩散模拟获取该排放源对受体点的贡献,从而实现污染物溯源。该方法主要针对一些重要排放源,不能实现全空间解析。大气污染物溯源在大气环境研究和空气质量管理中的发挥重要作用,通过污染物溯源,可以对受体点位的污染物来源进行解析,从而确定排放管控的重点区域,
发明内容
为解决现有技术中存在的问题,本发明提供了一种基于格点化扩散模拟的大气污染物溯源方法,针对所关心的区域,实现区域内任意受体点的大气污染物溯源,获取全空间格点化贡献比例解析,并可以实现实时解析,解决了上述背景技术中提到的问题。
为实现上述目的,本发明提供如下技术方案:一种基于格点化扩散模拟的大气污染物溯源方法,包括如下步骤:
一、利用中尺度气象模式计算高分辨率气象场数据;
二、针对目标区域划定网格,覆盖目标区域;
三、整理目标区域内的污染物排放数据,获得排放源强,利用G I S技术将污染物排放源强重采样到所划定的网格中;
四、在网格点上以面源排放的方式释放污染物,并开展扩散模拟,计算在目标区域内每个格点内的污染物浓度;
五、对每一个网格重复步骤四,保存计算结果;
六、对计算结果开展后处理,可以获取每个网格对目标区域内任意受体点的污染物浓度贡献值,同时可对每个网格的贡献值进行排名,识别重点管控网格。
优选的,所述步骤一中高分辨率气象场数据具体是指网格空间分辨率不低于3km×3km,时间分辨率是1小时。
优选的,所述步骤二中划分网格具体是通过地理信息系统将目标区域排放量进行1km×1km的空间进行划分网格。
优选的,所述的步骤三中整理目标区域内的污染物排放数据具体是指收集、获取辖区内的大气颗粒物排放相关的活动水平数据,建立区域颗粒物排放信息库。
优选的,所述的数据主要包括工业企业的化石燃料固定燃烧、工艺过程、交通道路扬尘、工地扬尘和非道路移动机械等。
优选的,所述步骤三中的获得排放源强具体是指获取区域内污染时段不同网格的气团滞留时间以及污染来源识别计算,获得排放源强。
优选的,所述步骤六中的后处理得到污染物浓度贡献值具体是指根据污染源并利用扩散模型模拟得到污染物浓度高值区域的污染气团的概率分布,扩散模型模拟的空间分辨率可以任意设置,模拟污染物浓度高值区前推几小时的气团轨迹,从而可得到污染气团的概率分布,将排放量与气团概率分布相乘得到污染源对不同网格的污染物浓度的贡献值。
本发明的有益效果是:本发明方法能够实现任意受体点的污染物来源全空间实时解析,识别重点源区,为空气质量网格化管理提供依据,可以确定每个控制区的最大排放量及允许排放指标,并确定环境承载力,根据环境承载力制定减排措施,确定容量分配并做好监督管理。
附图说明
图1为本发明实施例溯源方法流程示意图;
图2为本发明实施例网格划分示意图;
图3为本发明实施例各网格点PM2.5浓度贡献示意图;
图4为本发明实施例各网格点PM2.5浓度贡献排名示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
请参阅图1-4,本发明提供一种技术方案:一种基于格点化扩散模拟的大气污染物溯源方法,具体流程如图1所示,包括如下步骤:
1)利用中尺度气象模式计算高分辨率气象场数据;
高分辨率气象场数据具体是指网格空间分辨率不低于3km×3km,时间分辨率是1小时,气象领域将2km到2000km的天气系统称为中尺度系统。
2)针对目标区域划定网格,覆盖目标区域;
通过地理信息系统将目标区域排放量进行1km×1km的空间进行划分网格,并覆盖目标区域,如图2所示。
3)整理目标区域内的污染物排放数据,获得排放源强,利用GIS技术将污染物排放源强重采样到所划定的网格中;
收集、获取辖区内的大气颗粒物排放相关的活动水平数据,建立区域颗粒物排放信息库,主要包括工业企业的化石燃料固定燃烧、工艺过程、交通道路扬尘、工地扬尘和非道路移动机械等,并针对数据作面源处理,按照道路网络的空间分布确定排放分配的网格,不同道路排放的空间差异结合道路的车流量、车速等参数以体现排放的空间变异性,获取区域内污染时段不同网格的气团滞留时间以及污染来源识别计算,获得排放源强。
4)在网格点上以面源排放的方式释放污染物,建立扩散模型,并开展扩散模型模拟,由计算出的污染物排放源强,计算在目标区域内每个格点内的污染物浓度;
5)对每一个网格重复步骤4),保存计算结果;
6)对计算结果开展后处理,根据污染源并利用扩散模型模拟得到污染物浓度高值区域的污染气团的概率分布,扩散模型模拟的空间分辨率可以任意设置,模拟污染物浓度高值区前推几小时的气团轨迹,从而可得到污染气团的概率分布,将排放量与气团概率分布相乘得到污染源对不同网格的污染物浓度的贡献值,可以获取每个网格对目标区域内任意受体点的污染物浓度贡献值,同时可对每个网格的贡献值进行排名,识别重点管控网格。
以2020年10月24日四川省德阳市东山公园站点的PM2.5浓度解析为例(不限于东山公园站,区域内任意受体点均可解析),如图3和图4所示,可以获取每个网格实现对对目标区域内任意受体点的污染物浓度贡献值及排名情况,能识别重点管控网格,方法简单易行,解析全面快速,能有效地支持各地大气污染的靶向防治,具有极大的积极作用。
尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种基于格点化扩散模拟的大气污染物溯源方法,其特征在于,包括如下步骤:
一、利用中尺度气象模式计算高分辨率气象场数据;
二、针对目标区域划定网格,覆盖目标区域;
三、整理目标区域内的污染物排放数据,获得排放源强,利用GIS技术将污染物排放源强重采样到所划定的网格中;
四、在网格点上以面源排放的方式释放污染物,并开展扩散模拟,计算在目标区域内每个格点内的污染物浓度;
五、对每一个网格重复步骤四,保存计算结果;
六、对计算结果开展后处理,可以获取每个网格对目标区域内任意受体点的污染物浓度贡献值,同时可对每个网格的贡献值进行排名,识别重点管控网格。
2.根据权利要求1所述的基于格点化扩散模拟的大气污染物溯源方法,其特征在于:所述步骤一中高分辨率气象场数据具体指网格空间分辨率不低于3km×3km,时间分辨率是1小时。
3.根据权利要求1所述的基于格点化扩散模拟的大气污染物溯源方法,其特征在于:所述步骤二中划分网格具体是通过地理信息系统将目标区域排放量进行1km×1km的空间进行划分网格。
4.根据权利要求1所述的基于格点化扩散模拟的大气污染物溯源方法,其特征在于:所述的步骤三中整理目标区域内的污染物排放数据具体是指收集、获取辖区内的大气颗粒物排放相关的活动水平数据,建立区域颗粒物排放信息库。
5.根据权利要求4所述的基于格点化扩散模拟的大气污染物溯源方法,其特征在于:所述的数据主要包括工业企业的化石燃料固定燃烧、工艺过程、交通道路扬尘、工地扬尘和非道路移动机械等。
6.根据权利要求1所述的基于格点化扩散模拟的大气污染物溯源方法,其特征在于:所述步骤三中的获得排放源强具体是指获取区域内污染时段不同网格的气团滞留时间以及污染来源识别计算,获得排放源强。
7.根据权利要求1所述的基于格点化扩散模拟的大气污染物溯源方法,其特征在于:所述步骤六中的后处理得到污染物浓度贡献值具体是指根据污染源并利用扩散模型模拟得到污染物浓度高值区域的污染气团的概率分布,扩散模型模拟的空间分辨率可以任意设置,模拟污染物浓度高值区前推几小时的气团轨迹,从而可得到污染气团的概率分布,将排放量与气团概率分布相乘得到污染源对不同网格的污染物浓度的贡献值。
CN202110426596.1A 2021-04-20 2021-04-20 一种基于格点化扩散模拟的大气污染物溯源方法 Active CN113284244B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110426596.1A CN113284244B (zh) 2021-04-20 2021-04-20 一种基于格点化扩散模拟的大气污染物溯源方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110426596.1A CN113284244B (zh) 2021-04-20 2021-04-20 一种基于格点化扩散模拟的大气污染物溯源方法

Publications (2)

Publication Number Publication Date
CN113284244A true CN113284244A (zh) 2021-08-20
CN113284244B CN113284244B (zh) 2023-09-05

Family

ID=77276996

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110426596.1A Active CN113284244B (zh) 2021-04-20 2021-04-20 一种基于格点化扩散模拟的大气污染物溯源方法

Country Status (1)

Country Link
CN (1) CN113284244B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114818391A (zh) * 2022-06-28 2022-07-29 北京化工大学 一种基于多溯源模型的污染源浓度智能分析方法
CN115130831A (zh) * 2022-06-08 2022-09-30 北京工业大学 一种小尺度精细化大气污染溯源方法
CN116228502A (zh) * 2023-05-08 2023-06-06 中科三清科技有限公司 重点污染源排放量订正方法、装置、存储介质及电子设备

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150031577A (ko) * 2013-09-16 2015-03-25 주식회사 에니텍 대기오염물질 배출량 역산출 방법
CN105956664A (zh) * 2016-04-27 2016-09-21 浙江大学 一种河流点源突发污染事故溯源方法
CN109444232A (zh) * 2018-12-26 2019-03-08 苏州同阳科技发展有限公司 一种多通道智能化污染气体监测装置与扩散溯源方法
CN109522603A (zh) * 2018-10-18 2019-03-26 南京大学 基于云平台的车载拉格朗日实时大气污染溯源系统及方法
CN109878400A (zh) * 2018-10-19 2019-06-14 南京大学 模块化全功能拓展大气监测移动平台及监测使用方法
CN110095394A (zh) * 2019-05-27 2019-08-06 佛山市环境监测中心站 一种大气颗粒物污染快速溯源方法
CN110232198A (zh) * 2018-11-23 2019-09-13 深圳市城市交通规划设计研究中心有限公司 一种污染浓度的计算方法、计算装置及终端
CN112131739A (zh) * 2020-09-22 2020-12-25 南京创蓝科技有限公司 一种乡镇级别大气污染溯源的预报方法
CN112132450A (zh) * 2020-09-22 2020-12-25 南京创蓝科技有限公司 一种气态污染物的定位方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150031577A (ko) * 2013-09-16 2015-03-25 주식회사 에니텍 대기오염물질 배출량 역산출 방법
CN105956664A (zh) * 2016-04-27 2016-09-21 浙江大学 一种河流点源突发污染事故溯源方法
CN109522603A (zh) * 2018-10-18 2019-03-26 南京大学 基于云平台的车载拉格朗日实时大气污染溯源系统及方法
CN109878400A (zh) * 2018-10-19 2019-06-14 南京大学 模块化全功能拓展大气监测移动平台及监测使用方法
CN110232198A (zh) * 2018-11-23 2019-09-13 深圳市城市交通规划设计研究中心有限公司 一种污染浓度的计算方法、计算装置及终端
CN109444232A (zh) * 2018-12-26 2019-03-08 苏州同阳科技发展有限公司 一种多通道智能化污染气体监测装置与扩散溯源方法
CN110095394A (zh) * 2019-05-27 2019-08-06 佛山市环境监测中心站 一种大气颗粒物污染快速溯源方法
CN112131739A (zh) * 2020-09-22 2020-12-25 南京创蓝科技有限公司 一种乡镇级别大气污染溯源的预报方法
CN112132450A (zh) * 2020-09-22 2020-12-25 南京创蓝科技有限公司 一种气态污染物的定位方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
W. M. ANGEVINE 等: "Uncertainty in Lagrangian pollutant transport simulations due to meteorological uncertainty from a mesoscale WRF ensemble", GEOSCIENTIFIC MODEL DEVELOPMENT, vol. 07, no. 06, pages 2817 *
危浩 等: "特种污染物排放源强反演及可视化分析工具研发及Hg~(2+)分析示例", vol. 37, no. 02, pages 512 - 519 *
孙杰 等: "基于进化算法和水质模型的河口污染物溯源方法", 科学技术与工程, vol. 19, no. 14, pages 384 - 390 *
黄顺祥 等: "基于伴随方法的大气污染溯源", 科学通报, vol. 63, no. 16, pages 1594 - 1605 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115130831A (zh) * 2022-06-08 2022-09-30 北京工业大学 一种小尺度精细化大气污染溯源方法
CN115130831B (zh) * 2022-06-08 2024-02-27 北京工业大学 一种小尺度精细化大气污染溯源方法
CN114818391A (zh) * 2022-06-28 2022-07-29 北京化工大学 一种基于多溯源模型的污染源浓度智能分析方法
CN114818391B (zh) * 2022-06-28 2022-10-21 北京化工大学 一种基于多溯源模型的污染源浓度分析方法
CN116228502A (zh) * 2023-05-08 2023-06-06 中科三清科技有限公司 重点污染源排放量订正方法、装置、存储介质及电子设备
CN116228502B (zh) * 2023-05-08 2023-08-04 中科三清科技有限公司 重点污染源排放量订正方法、装置、存储介质及电子设备

Also Published As

Publication number Publication date
CN113284244B (zh) 2023-09-05

Similar Documents

Publication Publication Date Title
CN112418609B (zh) 一种基于面-网格-点的大气二次污染精准溯源方法
CN113284244A (zh) 一种基于格点化扩散模拟的大气污染物溯源方法
CN109522603B (zh) 基于云平台的车载拉格朗日实时大气污染溯源系统及方法
CN112749478B (zh) 一种基于高斯扩散模型大气污染溯源扩散分析系统及方法
CN110531030B (zh) 一种基于物联网大数据大气污染物分析的方法及装置
CN113624921B (zh) 多模式综合污染溯源方法
CN114356880A (zh) 基于数据驱动的小尺度区域大气污染物快速溯源方法
CN110909483A (zh) 一种基于网格化数据的点源大气污染物排放清单校验方法
CN115420854B (zh) 一种基于正后向模型结合的大气污染物溯源方法
CN104865353A (zh) 基于无人机的工业园区大气污染数据采集方法
CN110095394A (zh) 一种大气颗粒物污染快速溯源方法
CN111696369A (zh) 一种基于多源地理空间大数据的全市道路分时分车型交通流预测方法
CN112132450A (zh) 一种气态污染物的定位方法
KR101882874B1 (ko) 대기환경을 고려한 발전기 운영 계획 수립 방법 및 그 장치
CN115902114B (zh) 基于半定量法的小尺度大气污染溯源方法
CN112581107B (zh) 污染排放控制方法、装置及可存储介质
CN116626233A (zh) 基于多源数据融合的大气污染溯源方法、终端及存储介质
CN118013769B (zh) 一种基于WRF-Chem的大气污染物浓度预测方法
CN107767081B (zh) 一种燃煤电厂大气环境影响评价方法及装置
CN117291343A (zh) 一种基于空天地监测网络的大气污染防治管控系统
CN118169339B (zh) 基于多源数据的周边污染源溯源方法、装置、介质及设备
CN112541255A (zh) 一种大气污染排放来源自动识别系统及识别方法
CN115545565A (zh) 基于大气环境质量的园区排污总量管控方法和系统
CN116739222B (zh) 确定道路移动源的污染物浓度贡献值的方法与装置
CN117852764A (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