CN116228501B - 排污超标区域行业确定方法、装置、存储介质及电子设备 - Google Patents
排污超标区域行业确定方法、装置、存储介质及电子设备 Download PDFInfo
- Publication number
- CN116228501B CN116228501B CN202310510709.5A CN202310510709A CN116228501B CN 116228501 B CN116228501 B CN 116228501B CN 202310510709 A CN202310510709 A CN 202310510709A CN 116228501 B CN116228501 B CN 116228501B
- Authority
- CN
- China
- Prior art keywords
- industry
- determining
- monitoring
- site
- area
- 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
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
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Abstract
本公开涉及一种排污超标区域行业确定方法、装置、存储介质及电子设备,该方法通过确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量;根据所述多个子区域内多个行业分别对每个所述监测站点的污染物浓度的模拟贡献量确定每个所述监测站点的所述污染物浓度的模拟值;根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型;根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,所述目标区域行业信息,用于表征排污超标的子区域内的排污超标行业,能够快速排查出排污超标的子区域内的排污超标行业。
Description
技术领域
本公开涉及环境监测技术领域,具体地,涉及一种排污超标区域行业确定方法、装置、存储介质及电子设备。
背景技术
在预计发生重污染时通常会采取应急减排措施以改善空气质量是世界通用做法。但是在应急减排时,总会存在不按规定实施应急减排的现象存在,这就容易造成空气质量改善幅度偏低、减排效果不明显的现象出现。因此在实施应急减排时,核查应急减排落实情况、找出污染物排放存在异常的对象,对改善空气质量尤为重要。
发明内容
本公开的目的是提供一种排污超标区域行业确定方法、装置、存储介质及电子设备。
本公开第一方面提供一种排污超标区域行业确定方法,所述方法包括:
确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量;
根据所述多个子区域内多个行业分别对每个所述监测站点的污染物浓度的模拟贡献量确定每个所述监测站点的所述污染物浓度的模拟值;
根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型,所述站点类型用于表征站点污染物浓度是否存在异常;
根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,所述目标区域行业信息,用于表征排污超标的子区域内的排污超标行业。
可选地,所述站点类型包括异常站点和非异常站点,所述根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型,包括:
获取所述监测站点对所述污染物浓度的监测值;
在所述监测值与所述模拟值的差值大于或者等于预设差值阈值的情况下,确定所述监测站点为所述异常站点;
在所述监测值与所述模拟值的差值小于预设差值阈值的情况下,确定所述监测站点为所述非异常站点。
可选地,所述根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,包括:
确定所述排查区域内所述异常站点的第一数量;
根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息。
可选地,所述根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息,包括:
在所述第一数量大于或者等于第一阈值的情况下,确定所述排查区域内多个监测站点中所述监测值与所述模拟值的差值最小的第一监测站点,以及所述监测值与所述模拟值的差值最大的第二监测站点;
确定对所述第二监测站点的所述污染物浓度的模拟贡献量大于对所述第一监测站点的所述污染物浓度的模拟贡献量的待定区域行业信息,所述待定区域行业信息包括一个或多个区域行业的标识,所述区域行业用于表征一个子区域下的一个行业;
根据所述待定区域行业信息确定所述目标区域行业信息。
可选地,所述根据所述待定区域行业信息确定所述目标区域行业信息,包括:
确定所述待定区域行业信息中每个区域行业对所述第二监测站点的所述污染物浓度的模拟贡献量与对所述第一监测站点的所述污染物浓度的模拟贡献量的贡献量差值;
根据所述排查区域中多个区域行业的所述贡献量差值对所述多个区域行业进行排序,以得到目标顺序的区域行业序列;
根据所述区域行业序列确定所述目标区域行业信息。
可选地,所述根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息,包括:
在所述第一数量小于第一阈值,且大于1的情况下,针对每个异常站点,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合;确定第一数量的异常站点对应的第一数量个所述区域行业集合的交集,以得到所述目标区域行业信息;其中,所述第一阈值大于1。
可选地,所述根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息,包括:
在所述第一数量为1的情况下,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合,以得到所述目标区域行业信息。
可选地,所述确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量,包括:
确定网格化的目标大气污染源排放清单和行业排放清单,所述行业排放清单包括不同行业的污染物排放信息;
获取所述排查区域中的子区域信息,所述子区域信息包括每个子区域的位置信息以及每个子区域的标识信息;
确定未来时间段内的气象预报数据;
将所述目标大气污染源排放清单,所述行业排放清单,所述气象预报数据和所述子区域信息输入预设空气质量模式,以获取所述预设空气质量模式输出的所述排查区域中每个子区域内的每个行业对每个监测站点的污染物浓度的所述模拟贡献量。
本公开的第二方面提供一种排污超标区域行业确定装置,所述装置包括:
第一确定模块,被配置为确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量;
第二确定模块,被配置为根据所述多个子区域内多个行业分别对每个所述监测站点的污染物浓度的模拟贡献量确定每个所述监测站点的所述污染物浓度的模拟值;
第三确定模块,被配置为根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型,所述站点类型用于表征站点污染物浓度是否存在异常;
第四确定模块,被配置为根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,所述目标区域行业信息,用于表征排污超标的子区域内的排污超标行业。
可选地,所述站点类型包括异常站点和非异常站点,所述第三确定模块,被配置为:
获取所述监测站点对所述污染物浓度的监测值;
在所述监测值与所述模拟值的差值大于或者等于预设差值阈值的情况下,确定所述监测站点为所述异常站点;
在所述监测值与所述模拟值的差值小于预设差值阈值的情况下,确定所述监测站点为所述非异常站点。
可选地,所述第四确定模块,被配置为:
确定所述排查区域内所述异常站点的第一数量;
根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息。
可选地,所述第四确定模块,被配置为:
在所述第一数量大于或者等于第一阈值的情况下,确定所述排查区域内多个监测站点中所述监测值与所述模拟值的差值最小的第一监测站点,以及所述监测值与所述模拟值的差值最大的第二监测站点;
确定对所述第二监测站点的所述污染物浓度的模拟贡献量大于对所述第一监测站点的所述污染物浓度的模拟贡献量的待定区域行业信息,所述待定区域行业信息包括一个或多个区域行业的标识,所述区域行业用于表征一个子区域下的一个行业;
根据所述待定区域行业信息确定所述目标区域行业信息。
可选地,所述第四确定模块,被配置为:
确定所述待定区域行业信息中每个区域行业对所述第二监测站点的所述污染物浓度的模拟贡献量与对所述第一监测站点的所述污染物浓度的模拟贡献量的贡献量差值;
根据所述排查区域中多个区域行业的所述贡献量差值对所述多个区域行业进行排序,以得到目标顺序的区域行业序列;
根据所述区域行业序列确定所述目标区域行业信息。
可选地,所述第四确定模块,被配置为:
在所述第一数量小于第一阈值,且大于1的情况下,针对每个异常站点,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合;确定第一数量的异常站点对应的第一数量个所述区域行业集合的交集,以得到所述目标区域行业信息;其中,所述第一阈值大于1。
可选地,所述第四确定模块,被配置为:
在所述第一数量为1的情况下,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合,以得到所述目标区域行业信息。
可选地,所述第一确定模块,被配置为:
确定网格化的目标大气污染源排放清单和行业排放清单,所述行业排放清单包括不同行业的污染物排放信息;
获取所述排查区域中的子区域信息,所述子区域信息包括每个子区域的位置信息以及每个子区域的标识信息;
确定未来时间段内的气象预报数据;
将所述目标大气污染源排放清单,所述行业排放清单,所述气象预报数据和所述子区域信息输入预设空气质量模式,以获取所述预设空气质量模式输出的所述排查区域中每个子区域内的每个行业对每个监测站点的污染物浓度的所述模拟贡献量。
本公开的第三方面提供一种非临时性计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现以上第一方面所述方法的步骤。
本公开的第四方面提供一种电子设备,包括:
存储器,其上存储有计算机程序;
处理器,用于执行所述存储器中的所述计算机程序,以实现以上第一方面所述方法的步骤。
上述技术方案,通过确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量;根据所述多个子区域内多个行业分别对每个所述监测站点的污染物浓度的模拟贡献量确定每个所述监测站点的所述污染物浓度的模拟值;根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型;根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,所述目标区域行业信息,用于表征排污超标的子区域内的排污超标行业,能够快速排查出排污超标的子区域内的排污超标行业,有利于提升排污超标行业的排查效率,降低污染排查所需的人力成本。
本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本公开的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本公开,但并不构成对本公开的限制。在附图中:
图1是本公开一示例性实施例示出的一种排污超标区域行业确定方法的流程图;
图2是根据图1所示实施例示出的一种排污超标区域行业确定方法的流程图;
图3是根据图2所示实施例示出的一种排污超标区域行业确定方法的流程图;
图4是根据图2所示实施例示出的另一种排污超标区域行业确定方法的流程图;
图5是根据图1所示实施例示出的另一种排污超标区域行业确定方法的流程图;
图6是本公开一示例性实施例示出的一种排污超标区域行业确定装置的框图;
图7是根据一示例性实施例示出的一种电子设备的框图;
图8是根据一示例性实施例示出的另一种电子设备的框图。
具体实施方式
以下结合附图对本公开的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本公开,并不用于限制本公开。
需要说明的是,本公开中所有获取信号、信息或数据的动作都是在遵照所在地国家相应的数据保护法规政策的前提下,并获得由相应装置所有者给予授权的情况下进行的。
在详细介绍本公开的具体实施方式之前,首先对本公开的应用场景进行以下说明,本公开可以应用于污染物异常排放的排查过程中,例如排查应急减排企业对减排措施的落实情况,或者对存在异常排放的区域或行业的排查,能够有效保证排查效率的同时保证排查结果的准确度。目前相关技术中的通常是对区域内每个行业下的多个企业进行逐个排查,而对于企业的异常排放,通常所采用的方法有企业用电量对比法、在线监测法和现场检查法。企业用电量对比法,主要将企业排查时期的用电量和正常生产时的用电量进行对比,通过用电量的差异判断企业在排查时期是否存在异常排放情况,该方法通常存在企业实时用电量信息获取困难的问题,因此应用有限。在线监测方法是在企业排口安装在线监测仪器,实时监控企业排口的排放,所述在线检测方法由于安装在线监测设备的企业有限,且近年来不断爆出企业在线监测造假的信息,给各级监管工作带来巨大困扰,因此应用也比较有限。现场检查法,主要是对企业进行现场检查,判断企业是否存在异常排放的问题,但需要耗费大量人力,排查效率较低。
为了解决以上技术问题,本公开提供了一种排污超标区域行业确定方法、装置、存储介质及电子设备,该方法通过确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量;根据所述多个子区域内多个行业分别对每个所述监测站点的污染物浓度的模拟贡献量确定每个所述监测站点的所述污染物浓度的模拟值;根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型;根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,所述目标区域行业信息,用于表征排污超标的子区域内的排污超标行业,能够快速排查出排污超标的子区域内的排污超标行业,有利于提升排污超标行业的排查效率,降低污染排查所需的人力成本。
下面结合具体实施例对本公开的实施方式进行详细阐述。
图1是本公开一示例性实施例示出的一种排污超标区域行业确定方法的流程图;如图1所示,该排污超标区域行业确定方法可以包括以下步骤:
步骤101,确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量。
其中,该排查区域可以省、市、县等区域,也可以是地图上指定的区域,该子区域可以是组成该排查区域的区域,例如在该排查区域为一个省对应的区域的情况下,该子区域可以是市对应的区域、县对应的区域或者乡镇级对应的区域,也可以是根据指定划分方式划分的多个子区域,该污染物可以是SO2、NOx、CO、PM10、PM2.5、VOCs中的任一种。该行业可以是指从事国民经济中同性质的生产、服务或其他经济社会的经营单位或者个体的组织结构体系,例如可以包括化工行业、服装行业、建筑行业、农业、旅游业等,该行业的划分方式可以预先设定,例如可以预先设定不同企业对应的行业属性。
步骤102,根据所述多个子区域内多个行业分别对每个所述监测站点的污染物浓度的模拟贡献量确定每个所述监测站点的所述污染物浓度的模拟值。
示例地,若一个排查区域内包括m个子区域,n个行业,用表示第/>个子区域内的第/>种行业对第/>个监测站点的模拟贡献量,则该第/>个监测站点的模拟值/>可以通过计算得到。
步骤103,根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型。
其中,所述站点类型用于表征站点污染物浓度是否存在异常。
一种实施方式中,可以获取所述监测站点对所述污染物浓度的监测值;在所述监测值与所述模拟值的差值大于或者等于预设差值阈值的情况下,确定所述监测站点为所述异常站点;在所述监测值与所述模拟值的差值小于预设差值阈值的情况下,确定所述监测站点为所述非异常站点。
另一种实施方式中,可以预先设置异常判定阈值,在该模拟值大于或者等于异常判定阈值的情况下,确定所述监测站点为所述异常站点;在所述模拟值小于异常判定阈值的情况下,确定所述监测站点为所述非异常站点。
步骤104,根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,所述目标区域行业信息,用于表征排污超标的子区域内的排污超标行业。
本步骤一种可能的实施方式中,可以确定所述排查区域内所述异常站点的第一数量;根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息。以上所述的根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息,可以包括:
在所述第一数量大于或者等于第一阈值的情况下,确定所述排查区域内多个监测站点中所述监测值与所述模拟值的差值最小的第一监测站点,以及所述监测值与所述模拟值的差值最大的第二监测站点;确定对所述第二监测站点的所述污染物浓度的模拟贡献量大于对所述第一监测站点的所述污染物浓度的模拟贡献量的待定区域行业信息,所述待定区域行业信息包括一个或多个区域行业的标识,所述区域行业用于表征一个子区域下的一个行业;根据所述待定区域行业信息确定所述目标区域行业信息。
在所述第一数量小于第一阈值,且大于1的情况下,针对每个异常站点,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合;确定第一数量的异常站点对应的第一数量个所述区域行业集合的交集,以得到所述目标区域行业信息;其中,所述第一阈值大于1。
在所述第一数量为1的情况下,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合,以得到所述目标区域行业信息。
本步骤另一种可能的实施方式中,可以确定每个子区域内每个行业对异常站点的模拟贡献量,将模拟贡献量最大的预设数量个区域行业确定为排污异常的区域行业,生成包含该排污异常的区域行业的目标区域行业信息。
其中,该目标区域行业信息可以包括一个或多个区域行业标识信息,该区域行业标识信息可以是X子区域中的Y行业, X为子区域的标识信息,Y为行业的标识信息。
以上技术方案,能够快速排查出排污超标的子区域内的排污超标行业,有利于提升排污超标行业的排查效率,降低污染排查所需的人力成本。
图2是根据图1所示实施例示出的一种排污超标区域行业确定方法的流程图;如图2所示,图1中步骤104所述的根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,可以包括:
步骤1041,确定所述排查区域内所述异常站点的第一数量。
步骤1042,根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息。
一种实施方式中,可以通过图3中S11至S13所示步骤实施,图3是根据图2所示实施例示出的一种排污超标区域行业确定方法的流程图;如图3所示,该步骤1042,可以包括:
S11,在所述第一数量大于或者等于第一阈值的情况下,确定所述排查区域内多个监测站点中所述监测值与所述模拟值的差值最小的第一监测站点,以及所述监测值与所述模拟值的差值最大的第二监测站点。
S12,确定对所述第二监测站点的所述污染物浓度的模拟贡献量大于对所述第一监测站点的所述污染物浓度的模拟贡献量的待定区域行业信息。
其中,所述待定区域行业信息包括一个或多个区域行业的标识,所述区域行业用于表征一个子区域下的一个行业。
示例地,若所述排查区域内多个监测站点中所述监测值与所述模拟值的差值最大的第二监测站点,例如为监测站点1,表明该站点排污超标情况最严重,或者受部分区域行业未严格落实减排的影响最大,需尽快排查排污超标的行业,以降低该监测站点的监测浓度。同时所述排查区域内多个监测站点中所述监测值与所述模拟值的差值最小的第一监测站点,例如为监测站点2,该监测站点2的所述监测值与所述模拟值的差值可以是超出预设差值阈值的,可以是未超出预设差值阈值的,即该监测站点2可以是异常站点,也可以是非异常站点。若该监测站点1的监测值为,该监测站点1的模拟值为/>,该监测站点2的监测值为/>,该监测站点2的模拟值为/>,则所述监测站点1的监测值与模拟值差值,与所述监测站点2的监测值与模拟值差值满足关系:/>,进行不等式变换即可得到:/>,即/>>/>,再进行不等式变换可以得到:/>>0。其中/>为第j个子区域第k个行业的减排系数,该减排系数为大于1的系数,该/>为第j个子区域第k个行业对监测站点1的污染物浓度的模拟贡献量,/>为第j个子区域第k个行业对监测站点2的污染物浓度的模拟贡献量。
>0可分解为三部分,/>、和/>部分,因/>,所以/>部分对上式仅是负贡献,即使该部分中的某些子区域或行业排污超标(未按规定减排),仅会使上式小于0部分增大;/>部分对上式无贡献,即使该部分中的某些区域或行业未按规定减排,对上式没有影响;而/>部分对上式是正贡献,只有该部分中某些区域或行业未按规定减排才会使上式大于0部分增大,进而使上式成立。这表明必有中的某些区域行业未按规定减排,并且/>值越大,其未按规定进行减排的可能性越大,因此,可以将该/>的区域行业作为该排污超标的区域行业,根据该排污超标的区域行业的标识生成该待定区域行业信息。
S13,根据所述待定区域行业信息确定所述目标区域行业信息。
本步骤中,可以确定所述待定区域行业信息中每个区域行业对所述第二监测站点的所述污染物浓度的模拟贡献量与对所述第一监测站点的所述污染物浓度的模拟贡献量的贡献量差值;根据所述排查区域中多个区域行业的所述贡献量差值对所述多个区域行业进行排序,以得到目标顺序的区域行业序列;根据所述区域行业序列确定所述目标区域行业信息。
其中,该目标顺序可以是由大至小的顺序,也可以是由小至大的顺序,在该目标顺序为由大至小的顺序的情况下,可以将排序的前R个区域行业作为该排污超标的区域行业,R为预先指定的数量。
这样,按照对监测站点1和监测站点2贡献量差值从高到低的顺序对R个区域行业进行重点巡查,可快速锁定排污异常的企业,进而能够有效提升排查效率,节约排查过程所需的人力物力。
另一种实施方式中,可以通过图4中S14至S15所示步骤实施,图4是根据图2所示实施例示出的另一种排污超标区域行业确定方法的流程图;如图4所示,该步骤1042,还可以包括:
S14,在所述第一数量小于第一阈值,且大于1的情况下,针对每个异常站点,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合。
示例地,若监测站点、监测站点b为异常站点,监测站点c、监测站点d和监测站点e为非异常站点。对于监测站点/>来说,排污超标(未按规定减排)的区域行业对监测站点c、监测站点d和监测站点e的贡献较小,但对监测站点/>贡献较大,即未按规定落实减排的区域行业且对监测站点/>影响较大的区域行业均存在于满足/>、/>、条件的区域行业中,其中,该/>为第j个子区域第k个行业对监测站点/>的污染物浓度的模拟贡献量,/>为第j个子区域第k个行业对监测站点c的污染物浓度的模拟贡献量,/>为第j个子区域第k个行业对监测站点d的污染物浓度的模拟贡献量,为第j个子区域第k个行业对监测站点e的污染物浓度的模拟贡献量。
S15,确定第一数量的异常站点对应的第一数量个所述区域行业集合的交集,以得到所述目标区域行业信息。
其中,所述第一阈值大于1。
仍以上述S14中的示例为例,可以对的区域行业集合,/>的区域行业集合以及/>的区域行业集合取交集,从而可以有效确定出排污超标且对监测站点/>影响较大的区域行业;同理可得到对监测站点b影响较大且排污超标的区域行业。
再一种实施方式中,图1中步骤1042,还可以包括:
在所述第一数量为1的情况下,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合,以得到所述目标区域行业信息。
示例地,若监测测站点1为异常站点,监测站点2、监测站点3、监测站点4和监测站点5为非异常站点。即排污超标的区域行业(未按规定减排的区域行业)对监测站点2、监测站点3、监测站点4和监测站点5的贡献均较小,但是对监测站点1献较大,即排污超标的区域行业为满足、/>、/>、/>条件的区域行业,对上述4种情景取交集,可以明显缩小排污超标且对监测站点1影响较大的区域或行业,进而得到该目标区域行业信息。
以上技术方案,能够快速排查出排污超标的子区域内的排污超标行业,有利于提升排污超标行业的排查效率,降低污染排查所需的人力成本。
图5是根据图1所示实施例示出的另一种排污超标区域行业确定方法的流程图;如图5所示,图1中步骤101所述的所述确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量,可以包括:
步骤1011,确定网格化的目标大气污染源排放清单和行业排放清单,所述行业排放清单包括不同行业的污染物排放信息。
其中,该行业的划分方式可以预先设定。
示例地,可以根据减排措施对大气污染源排放清单进行调整,生成按规定落实减排的大气污染源排放清单;使用大气污染源排放清单生成的落实减排后的大气污染源排放清单和各行业的排放清单,行业划分可根据实际需求设置;再使用SMOKE(大气排放源清单处理模型)或其他排放源处理软件生成网格化的目标大气污染源排放清单和行业排放清单。
步骤1012,获取所述排查区域中的子区域信息,所述子区域信息包括每个子区域的位置信息以及每个子区域的标识信息。
本步骤中,可以使用ARCgis(GIS平台中的一种软件名称)或其他软件根据该子区域信息对各区域进行标记生成区域标记文件,进行区域标记时可根据需求按行政区划进行标记,例如对各区县或者各乡镇进行标记。
步骤1013,确定未来时间段内的气象预报数据。
本步骤中,可以使用WRF(The Weather Research and Forecasting Model,天气预报模式)气象模式或GFS(Google Earth Engine,全球气象预报系统模式)气象预测数据生成该未来时间段内的气象预报数据。
步骤1014,将所述目标大气污染源排放清单,所述行业排放清单,所述气象预报数据和所述子区域信息输入预设空气质量模式,以获取所述预设空气质量模式输出的所述排查区域中每个子区域内的每个行业对每个监测站点的污染物浓度的所述模拟贡献量。
其中,该预设空气质量模型可以是CAMQ(描述大气污染物空气质量的复杂数值模型)、CAMx(综合空气质量模式)、WRF-CHEM(区域大气动力-化学耦合模式)、NAQPMS(嵌套网格空气质量预报系统)等三代空气质量模式,也可以是CALPUFF(非稳态拉格朗日烟团模型系统)、AERMOD(稳态烟羽模型)、ADMS(三维高斯型大气扩散模型)等二代空气质量模式。
这样,通过步骤1011至步骤1014,能够有效确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量,能够为排污超标的区域行业的确定提供可靠的数据依据。
图6是本公开一示例性实施例示出的一种排污超标区域行业确定装置的框图;如图6所示,所述装置可以包括:
第一确定模块601,被配置为确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量;
第二确定模块602,被配置为根据所述多个子区域内多个行业分别对每个所述监测站点的污染物浓度的模拟贡献量确定每个所述监测站点的所述污染物浓度的模拟值;
第三确定模块603,被配置为根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型,所述站点类型用于表征站点污染物浓度是否存在异常;
第四确定模块604,被配置为根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,所述目标区域行业信息,用于表征排污超标的子区域内的排污超标行业。
以上技术方案,通过第一确定模块601、第二确定模块602、第三确定模块603和第四确定模块604能够快速排查出排污超标的子区域内的排污超标行业,有利于提升排污超标行业的排查效率,降低污染排查所需的人力成本。
可选地,所述站点类型包括异常站点和非异常站点,所述第三确定模块603,被配置为:
获取所述监测站点对所述污染物浓度的监测值;
在所述监测值与所述模拟值的差值大于或者等于预设差值阈值的情况下,确定所述监测站点为所述异常站点;
在所述监测值与所述模拟值的差值小于预设差值阈值的情况下,确定所述监测站点为所述非异常站点。
可选地,所述第四确定模块604,被配置为:
确定所述排查区域内所述异常站点的第一数量;
根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息。
可选地,所述第四确定模块,被配置为:
在所述第一数量大于或者等于第一阈值的情况下,确定所述排查区域内多个监测站点中所述监测值与所述模拟值的差值最小的第一监测站点,以及所述监测值与所述模拟值的差值最大的第二监测站点;
确定对所述第二监测站点的所述污染物浓度的模拟贡献量大于对所述第一监测站点的所述污染物浓度的模拟贡献量的待定区域行业信息,所述待定区域行业信息包括一个或多个区域行业的标识,所述区域行业用于表征一个子区域下的一个行业;
根据所述待定区域行业信息确定所述目标区域行业信息。
可选地,所述第四确定模块604,被配置为:
确定所述待定区域行业信息中每个区域行业对所述第二监测站点的所述污染物浓度的模拟贡献量与对所述第一监测站点的所述污染物浓度的模拟贡献量的贡献量差值;
根据所述排查区域中多个区域行业的所述贡献量差值对所述多个区域行业进行排序,以得到目标顺序的区域行业序列;
根据所述区域行业序列确定所述目标区域行业信息。
可选地,所述第四确定模块604,被配置为:
在所述第一数量小于第一阈值,且大于1的情况下,针对每个异常站点,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合;确定第一数量的异常站点对应的第一数量个所述区域行业集合的交集,以得到所述目标区域行业信息;其中,所述第一阈值大于1。
可选地,所述第四确定模块604,被配置为:
在所述第一数量为1的情况下,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合,以得到所述目标区域行业信息。
可选地,所述第一确定模块601,被配置为:
确定网格化的目标大气污染源排放清单和行业排放清单,所述行业排放清单包括不同行业的污染物排放信息;
获取所述排查区域中的子区域信息,所述子区域信息包括每个子区域的位置信息以及每个子区域的标识信息;
确定未来时间段内的气象预报数据;
将所述目标大气污染源排放清单,所述行业排放清单,所述气象预报数据和所述子区域信息输入预设空气质量模式,以获取所述预设空气质量模式输出的所述排查区域中每个子区域内的每个行业对每个监测站点的污染物浓度的所述模拟贡献量。
以上技术方案,通过第一确定模块601能够有效确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量,能够为排污超标的区域行业的确定提供可靠的数据依据。
关于上述实施例中的装置,其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述,此处将不做详细阐述说明。
图7是根据一示例性实施例示出的一种电子设备的框图。如图7所示,该第一电子设备700可以包括:第一处理器701,第一存储器702。该第一电子设备700还可以包括多媒体组件703,第一输入/输出接口704,以及第一通信组件705中的一者或多者。
其中,第一处理器701用于控制该第一电子设备700的整体操作,以完成上述的排污超标区域行业确定方法中的全部或部分步骤。第一存储器702用于存储各种类型的数据以支持在该第一电子设备700的操作,这些数据例如可以包括用于在该第一电子设备700上操作的任何应用程序或方法的指令,以及应用程序相关的数据,例如联系人数据、收发的消息、图片、音频、视频等等。该第一存储器702可以由任何类型的易失性或非易失性存储设备或者它们的组合实现,例如静态随机存取存储器(Static Random Access Memory,简称SRAM),电可擦除可编程只读存储器(Electrically Erasable Programmable Read-OnlyMemory,简称EEPROM),可擦除可编程只读存储器(Erasable Programmable Read-OnlyMemory,简称EPROM),可编程只读存储器(Programmable Read-Only Memory,简称PROM),只读存储器(Read-Only Memory,简称ROM),磁存储器,快闪存储器,磁盘或光盘。多媒体组件703可以包括屏幕和音频组件。其中屏幕例如可以是触摸屏,音频组件用于输出和/或输入音频信号。例如,音频组件可以包括一个麦克风,麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在第一存储器702或通过第一通信组件705发送。音频组件还包括至少一个扬声器,用于输出音频信号。第一输入/输出接口704为第一处理器701和其他接口模块之间提供接口,上述其他接口模块可以是键盘,鼠标,按钮等。这些按钮可以是虚拟按钮或者实体按钮。第一通信组件705用于该第一电子设备700与其他设备之间进行有线或无线通信。无线通信,例如Wi-Fi,蓝牙,近场通信(Near Field Communication,简称NFC),2G、3G、4G、NB-IOT、eMTC、或其他5G等等,或它们中的一种或几种的组合,在此不做限定。因此相应的该第一通信组件705可以包括:Wi-Fi模块,蓝牙模块,NFC模块等等。
在一示例性实施例中,第一电子设备700可以被一个或多个应用专用集成电路(Application Specific Integrated Circuit,简称ASIC)、数字信号处理器(DigitalSignal Processor,简称DSP)、数字信号处理设备(Digital Signal Processing Device,简称DSPD)、可编程逻辑器件(Programmable Logic Device,简称PLD)、现场可编程门阵列(Field Programmable Gate Array,简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现,用于执行上述的排污超标区域行业确定方法。
在另一示例性实施例中,还提供了一种包括程序指令的计算机可读存储介质,该程序指令被处理器执行时实现上述的排污超标区域行业确定方法的步骤。例如,该计算机可读存储介质可以为上述包括程序指令的第一存储器702,上述程序指令可由第一电子设备700的第一处理器701执行以完成上述的排污超标区域行业确定方法。
图8是根据一示例性实施例示出的另一种电子设备的框图。例如,第二电子设备800可以被提供为一服务器。参照图8,第二电子设备800包括第二处理器822,其数量可以为一个或多个,以及第二存储器832,用于存储可由第二处理器822执行的计算机程序。第二存储器832中存储的计算机程序可以包括一个或一个以上的每一个对应于一组指令的模块。此外,第二处理器822可以被配置为执行该计算机程序,以执行上述的排污超标区域行业确定方法。
另外,第二电子设备800还可以包括电源组件826和第二通信组件850,该电源组件826可以被配置为执行第二电子设备800的电源管理,该第二通信组件850可以被配置为实现第二电子设备800的通信,例如,有线或无线通信。此外,该第二电子设备800还可以包括第二输入/输出接口858。第二电子设备800可以操作基于存储在第二存储器832的操作系统。
在另一示例性实施例中,还提供了一种包括程序指令的计算机可读存储介质,该程序指令被处理器执行时实现上述的排污超标区域行业确定方法的步骤。例如,该非临时性计算机可读存储介质可以为上述包括程序指令的第二存储器832,上述程序指令可由第二电子设备800的第二处理器822执行以完成上述的排污超标区域行业确定方法。
在另一示例性实施例中,还提供一种计算机程序产品,该计算机程序产品包含能够由可编程的装置执行的计算机程序,该计算机程序具有当由该可编程的装置执行时用于执行上述的排污超标区域行业确定方法的代码部分。
以上结合附图详细描述了本公开的优选实施方式,但是,本公开并不限于上述实施方式中的具体细节,在本公开的技术构思范围内,可以对本公开的技术方案进行多种简单变型,这些简单变型均属于本公开的保护范围。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合,为了避免不必要的重复,本公开对各种可能的组合方式不再另行说明。
此外,本公开的各种不同的实施方式之间也可以进行任意组合,只要其不违背本公开的思想,其同样应当视为本公开所公开的内容。
Claims (10)
1.一种排污超标区域行业确定方法,其特征在于,所述方法包括:
确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量;
根据所述多个子区域内多个行业分别对每个所述监测站点的污染物浓度的模拟贡献量确定每个所述监测站点的所述污染物浓度的模拟值;
根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型,所述站点类型用于表征站点污染物浓度是否存在异常;
根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,所述目标区域行业信息,用于表征排污超标的子区域内的排污超标行业;
所述站点类型包括异常站点和非异常站点,所述确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量,包括:
确定网格化的目标大气污染源排放清单和行业排放清单,所述行业排放清单包括不同行业的污染物排放信息;
获取所述排查区域中的子区域信息,所述子区域信息包括每个子区域的位置信息以及每个子区域的标识信息;
确定未来时间段内的气象预报数据;
将所述目标大气污染源排放清单,所述行业排放清单,所述气象预报数据和所述子区域信息输入预设空气质量模式,以获取所述预设空气质量模式输出的所述排查区域中每个子区域内的每个行业对每个监测站点的污染物浓度的所述模拟贡献量。
2.根据权利要求1所述的方法,其特征在于,所述根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型,包括:
获取所述监测站点对所述污染物浓度的监测值;
在所述监测值与所述模拟值的差值大于或者等于预设差值阈值的情况下,确定所述监测站点为所述异常站点;
在所述监测值与所述模拟值的差值小于预设差值阈值的情况下,确定所述监测站点为所述非异常站点。
3.根据权利要求2所述的方法,其特征在于,所述根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,包括:
确定所述排查区域内所述异常站点的第一数量;
根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息。
4.根据权利要求3所述的方法,其特征在于,所述根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息,包括:
在所述第一数量大于或者等于第一阈值的情况下,确定所述排查区域内多个监测站点中所述监测值与所述模拟值的差值最小的第一监测站点,以及所述监测值与所述模拟值的差值最大的第二监测站点;
确定对所述第二监测站点的所述污染物浓度的模拟贡献量大于对所述第一监测站点的所述污染物浓度的模拟贡献量的待定区域行业信息,所述待定区域行业信息包括一个或多个区域行业的标识,所述区域行业用于表征一个子区域下的一个行业;
根据所述待定区域行业信息确定所述目标区域行业信息。
5.根据权利要求4所述的方法,其特征在于,所述根据所述待定区域行业信息确定所述目标区域行业信息,包括:
确定所述待定区域行业信息中每个区域行业对所述第二监测站点的所述污染物浓度的模拟贡献量与对所述第一监测站点的所述污染物浓度的模拟贡献量的贡献量差值;
根据所述排查区域中多个区域行业的所述贡献量差值对所述多个区域行业进行排序,以得到目标顺序的区域行业序列;
根据所述区域行业序列确定所述目标区域行业信息。
6.根据权利要求3所述的方法,其特征在于,所述根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息,包括:
在所述第一数量小于第一阈值,且大于1的情况下,针对每个异常站点,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合;确定第一数量的异常站点对应的第一数量个所述区域行业集合的交集,以得到所述目标区域行业信息;其中,所述第一阈值大于1。
7.根据权利要求3所述的方法,其特征在于,所述根据所述第一数量确定所述排查区域中排污超标的目标区域行业信息,包括:
在所述第一数量为1的情况下,确定所述排查区域内的多个区域行业中对所述异常站点的污染物浓度的模拟贡献量大于对每个非异常站点的污染物浓度的模拟贡献量的区域行业集合,以得到所述目标区域行业信息。
8.一种排污超标区域行业确定装置,其特征在于,所述装置包括:
第一确定模块,被配置为确定排查区域中多个子区域内多个行业分别对所述排查区域内每个监测站点的污染物浓度的模拟贡献量;
第二确定模块,被配置为根据所述多个子区域内多个行业分别对每个所述监测站点的污染物浓度的模拟贡献量确定每个所述监测站点的所述污染物浓度的模拟值;
第三确定模块,被配置为根据每个所述监测站点的所述污染物浓度的模拟值确定所述监测站点的站点类型,所述站点类型用于表征站点污染物浓度是否存在异常;
第四确定模块,被配置为根据所述站点类型和每个所述监测站点的所述污染物浓度的模拟值确定所述排查区域中排污超标的目标区域行业信息,所述目标区域行业信息,用于表征排污超标的子区域内的排污超标行业;
所述站点类型包括异常站点和非异常站点,所述第一确定模块,被配置为:
确定网格化的目标大气污染源排放清单和行业排放清单,所述行业排放清单包括不同行业的污染物排放信息;
获取所述排查区域中的子区域信息,所述子区域信息包括每个子区域的位置信息以及每个子区域的标识信息;
确定未来时间段内的气象预报数据;
将所述目标大气污染源排放清单,所述行业排放清单,所述气象预报数据和所述子区域信息输入预设空气质量模式,以获取所述预设空气质量模式输出的所述排查区域中每个子区域内的每个行业对每个监测站点的污染物浓度的所述模拟贡献量。
9.一种非临时性计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现权利要求1-7中任一项所述方法的步骤。
10.一种电子设备,其特征在于,包括:
存储器,其上存储有计算机程序;
处理器,用于执行所述存储器中的所述计算机程序,以实现权利要求1-7中任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310510709.5A CN116228501B (zh) | 2023-05-08 | 2023-05-08 | 排污超标区域行业确定方法、装置、存储介质及电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310510709.5A CN116228501B (zh) | 2023-05-08 | 2023-05-08 | 排污超标区域行业确定方法、装置、存储介质及电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116228501A CN116228501A (zh) | 2023-06-06 |
CN116228501B true CN116228501B (zh) | 2023-07-21 |
Family
ID=86580937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310510709.5A Active CN116228501B (zh) | 2023-05-08 | 2023-05-08 | 排污超标区域行业确定方法、装置、存储介质及电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116228501B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117455708B (zh) * | 2023-12-22 | 2024-03-01 | 山西亚鑫新能科技有限公司 | 一种无组织排放智能治理系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10830922B2 (en) * | 2015-10-28 | 2020-11-10 | International Business Machines Corporation | Air quality forecast by adapting pollutant emission inventory |
CN111814111B (zh) * | 2020-05-25 | 2021-07-27 | 南京大学环境规划设计研究院集团股份公司 | 一种工业园区大气污染物溯源方法 |
CN112417788A (zh) * | 2020-11-30 | 2021-02-26 | 重庆市生态环境大数据应用中心 | 基于大数据的水环境污染分析系统及方法 |
CN113436045A (zh) * | 2021-06-15 | 2021-09-24 | 深圳中兴网信科技有限公司 | 大气污染多尺度融合溯源方法、装置和可读存储介质 |
CN115759345A (zh) * | 2022-10-21 | 2023-03-07 | 深圳市生态环境智能管控中心 | 智慧环保平台大气环境管理系统的质量预测模块及方法 |
-
2023
- 2023-05-08 CN CN202310510709.5A patent/CN116228501B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN116228501A (zh) | 2023-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111474307A (zh) | 污染物溯源方法、装置、计算机设备及存储介质 | |
Nguyen et al. | Development of the typical driving cycle for buses in Hanoi, Vietnam | |
CN116228501B (zh) | 排污超标区域行业确定方法、装置、存储介质及电子设备 | |
CN113777236B (zh) | 基于排放源的空气质量监控方法和装置 | |
CN110426493B (zh) | 空气质量监测数据校准方法、装置、设备和存储介质 | |
CN114329245A (zh) | 一种空气污染的溯源方法、装置、服务器及存储介质 | |
CN110807725A (zh) | 一种大气污染溯源方法、装置、计算设备和介质 | |
US20160370333A1 (en) | Generating fine resolution air pollution estimates | |
CN116227749A (zh) | 污染物排放量的确定方法、装置、存储介质及电子设备 | |
CN110348746B (zh) | 基于单个污染源的空气质量影响评估方法和装置 | |
CN112344990A (zh) | 环境异常监测方法、装置、设备及存储介质 | |
CN109115949B (zh) | 基于大数据的污染溯源方法及计算机可读介质 | |
CN112800165B (zh) | 一种基于聚类算法的产业集群定位方法、装置及电子设备 | |
CN116228502A (zh) | 重点污染源排放量订正方法、装置、存储介质及电子设备 | |
CN116757366B (zh) | 高污染排放污染源选址的方法、装置、介质及电子设备 | |
CN115238245B (zh) | 污染物监测方法、装置、存储介质及电子设备 | |
CN116148951B (zh) | 空气目标污染物组分预报模型训练方法、预报方法及装置 | |
US10288547B2 (en) | Facility state analysis device, analysis method for facility state, storage medium, and facility management system | |
CN110555586A (zh) | 一种基于热点网格的生态监测方法及装置 | |
CN116307258B (zh) | 污染源异常排放确定方法、装置、存储介质及电子设备 | |
CN117610895B (zh) | 重点污染源管控时间的确定方法、装置、电子设备及介质 | |
Zhu et al. | A data-driven approach for optimal design of integrated air quality monitoring network in a chemical cluster | |
CN113763221A (zh) | 一种大气污染的溯源方法、装置、电子设备及存储介质 | |
CN116739191B (zh) | 热点网格识别方法、装置、存储介质及电子设备 | |
CN116776645B (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 |