CN111727447A - 一种选择污染物治理措施的方法 - Google Patents
一种选择污染物治理措施的方法 Download PDFInfo
- Publication number
- CN111727447A CN111727447A CN201980003511.XA CN201980003511A CN111727447A CN 111727447 A CN111727447 A CN 111727447A CN 201980003511 A CN201980003511 A CN 201980003511A CN 111727447 A CN111727447 A CN 111727447A
- Authority
- CN
- China
- Prior art keywords
- measure
- pollutant
- measures
- source
- standard
- 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
Links
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 113
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 111
- 238000000034 method Methods 0.000 title claims abstract description 63
- 230000000694 effects Effects 0.000 claims description 22
- 238000012544 monitoring process Methods 0.000 claims description 17
- 239000000356 contaminant Substances 0.000 claims description 16
- 238000012216 screening Methods 0.000 claims description 9
- 239000000203 mixture Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 239000000428 dust Substances 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims description 5
- 238000012163 sequencing technique Methods 0.000 claims 1
- 238000003915 air pollution Methods 0.000 abstract description 22
- 238000010187 selection method Methods 0.000 abstract description 10
- 238000005516 engineering process Methods 0.000 abstract description 5
- 230000007547 defect Effects 0.000 abstract 1
- 239000003570 air Substances 0.000 description 24
- 230000009471 action Effects 0.000 description 23
- 239000000809 air pollutant Substances 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 231100001243 air pollutant Toxicity 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 9
- 238000012937 correction Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 230000002265 prevention Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 239000000443 aerosol Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000003245 coal Substances 0.000 description 3
- 239000005431 greenhouse gas Substances 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 238000003911 water pollution Methods 0.000 description 3
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 230000003137 locomotive effect Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000005406 washing Methods 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 229910002651 NO3 Inorganic materials 0.000 description 1
- NHNBFGGVMKEFGY-UHFFFAOYSA-N Nitrate Chemical compound [O-][N+]([O-])=O NHNBFGGVMKEFGY-UHFFFAOYSA-N 0.000 description 1
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 239000012080 ambient air Substances 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 239000002283 diesel fuel Substances 0.000 description 1
- 239000012717 electrostatic precipitator Substances 0.000 description 1
- 239000010419 fine particle Substances 0.000 description 1
- 239000003502 gasoline Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000005571 horizontal transmission Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
- 230000005570 vertical transmission Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
- Y02A50/20—Air quality improvement or preservation, e.g. vehicle emission control or emission reduction by using catalytic converters
-
- 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
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Processing Of Solid Wastes (AREA)
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Sampling And Sample Adjustment (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
Abstract
针对背景技术中对污染治理措施选取方式的不足,本发明提供了一种高效、快速、便捷、科学的空气污染治理措施选择方法。该方法是一种根据城市空气污染特征为城市优选空气污染防治措施的动态选择方法。该方法依据城市PM2.5、O3、SO2、NO2等污染物浓度超标率和其相对应的排放来源构成(根据排放清单得出)对空气污染治理措施进行选取,通过为措施库中各措施进行评分,根据分数大小为各城市优选适应其空气质量现状的措施。
Description
技术领域
本发明涉及一种选择环境污染治理措施的方法,属于环境治理领域。
背景技术
2016年新大气法设立了城市空气质量限期达标机制,城市成为空气质量改善的主体,然而目前我国仍有64%以上的城市没有实现空气质量达标,大气污染防治形势依然严峻。随着我国大气污染防治的不断深入,许多治理工作都进入深水区,政策制定和执行的难度都将会不断升级,对城市空气质量的精细化管理要求将越来越高,因此以城市为单元的定制化措施需求也越来越迫切。
目前,“底数不清”的状况是制约我国大气污染防治工作的重要瓶颈。很多城市的决策者对主要污染物排放总量、时空分布、行业贡献、减排潜力等信息掌握不足,在选择大气污染防治相关措施时,脱离当地实际情况,缺乏有效的科学支撑,措施不能与当地实际的空气污染状况想匹配,难以实现精准治霾。
目前已有一些控制质量预报和评估***,用于给相关决策人员提供建议,如CMAQ(Community Multi-scale Air Quality),CMAQ是美国环保署(EPA)在拉格朗日轨迹模型和欧拉网格模型后提出的第三代空气质量预报和评估***,该模型是在“一个大气”理论的指导下,以中尺度气象模式和SMOKE(Spare Matrix Operator Kerenl Emission)等源排放模型为依托,考虑了大气污染过程中水平传输、垂直传输、扩散过程、源排放、化学反应和去除过程等对污染物浓度的影响,将复杂空气污染状况进行综合处理。模型可模拟平流传输、湍流扩散、气相化学反应、气溶胶动力学、排放过程、沉降过程、云过程和液相过程,可用于评价大气中细颗粒物、对流层臭氧、气溶胶以及酸沉降污染水平。但是该***主要用于空气质量预报,对措施选择有一定参考,但没有行程***化的评估方式。
CMAx(Comprehensive Air Quality Model with Extensions)模型同样是欧拉型化学传输模型,在“一个大气”理念的指导下考虑气-液-固多相化学机制,利用气象场模拟结果,通过SMOKE源排放模型对源排放清单进行处理,最后CAMx模型对污染物浓度进行模拟,与models3-CMAQ模型不同的是CAMx模型具有双向嵌套的网格结构,可以多重网格同时进行计算,在时间范围和空间范围内模拟的更精细。CAMx还包括多种分析工具如臭氧来源解析技术(OSAT)、颗粒物来源追踪技术(PAST)、网格烟羽模块(PiG)等。CAMx工具同样主要是用于污染预测,并辅助一些污染来源解析工具,也不能有效的选择污染治理的措施。
WRF-Chem(Weather Research and Forecasting model coupled withChemistry)模型也是大气环境治理领域常用的技术,该模型的气象模式和化学传输模块使用相同的格点、时间步长、传输方案和物理方案,避免因差值等造成的误差,同时二者为同步计算,在时间和空间分辨率上完成耦合,实现真正的在线传输,从而完成对太阳辐射、大气动力和气溶胶化学等多过程的耦合和反馈。WRF-Chem可以对大气质量进行预测,也同样无法高效科学的选择污染治理措施。
以上这些用于治理措施的技术,主要用于预测空气质量,不能根据污染特征针对性的为决策者选择合适的措施或措施组合。使用这些技术需要运用复杂的数学、物理以及化学模型,需要专业的工作人员来进行实施,一次预测过程还需要使用专用的大型计算机,耗时长成本高。
用于污染措施选取的相关论文也比较少,最相关的一篇是【深圳交通拥堵综合治理的模式选择及策略措施】,出自2016年中国城市交通规划年会论文集,这篇论文主要介绍的是在交通治理方面的模式选择及策略措施的方法,该方法首先进行调研,收集其他的治理方法,经过人为的评估最后制定并决策出交通拥堵治理策略措施,并没有量化的选取指标,以及科学的论证方法。
在相关专利方面,没有检索到用于污染治理措施选择的方法专利。
因此根据我们的检索,对目前的大气污染治理措施的选取问题上,主要依靠人工的经验进行,没有一个科学的、高效的技术手段对相关的污染治理措施进行评估和选取。
发明内容
相关术语
排放清单编号i:i=1,2,3,4……,代表排放清单中的类别编号,如道路移动源、非道路移动源、扬尘源、工业源等。
排放清单编号(i)示例表
i值 | 1 | 2 | 3 | 4 |
污染源类别 | 道路移动源 | 非道路移动源 | 扬尘源 | 工业源 |
措施清单编号j:j=A,B,C,D……代表措施编号,如升级发动机机车:柴油电动混合动力机车、飞机地面支持设备:电力替代燃料、岸基电力、洗煤措施可以治理的污染物编号k:k=1,2,3,4……,k值代表为该措施可以治理的污染物,如NOX、PM2.5、SO2、O3。
可治理的污染物编号(k)示例表
k值 | 1 | 2 | 3 | 4 |
措施 | NO<sub>X</sub> | PM<sub>2.5</sub> | SO<sub>2</sub> | O<sub>3</sub> |
措施的特征编号α:α=1,2,3……代表措施的特征,如实施周期特征,资金需求特征。
措施类别权重Wi系数:Wi表示一类排放措施的权重系数,这个权重系数与污染物超标情况和来源构成有关,城市污染物超标率越高、措施库所对应排放源类别(措施库的分类与排放清单中排放源分类一致)在污染物来源构成中占比越大,对本权重系数贡献越大;如W1为机动车污染治理措施权重,W2为燃煤污染治理措施权重。
Fj:措施评分的分数。
超标率εk:εk为污染物年均浓度监测值与国家标准中的污染物年均浓度限值的比值。
污染物浓度限值Sk:Sk表示国家空气质量标准中规定的空气污染物年均浓度限值。污染物浓度监测值Tk:Sk表示空气污染物浓度的监测值。
根据εk值判断各污染物超标情况,不同的超标情况对应的不同的超标系数ek。
目前很多城市的决策者对主要污染物排放总量、时空分布、行业贡献、减排潜力等信息掌握不足,在选择大气污染防治相关措施时,脱离当地实际情况,缺乏有效的科学支撑,措施不能与当地实际的空气污染状况相匹配。
减少大气污染的措施目前已经基本完善,但是一个城市或者一片区域为了实现精准有效的治理雾霾,需要在有限的资源下对治理措施进行合理、高效、准确的选取,同时选取的措施需要满足一定的时间资源、成本限制。目前的措施选择方式主要依靠人工的经验进行,没有科学依据,容易造成决策失误,资源使用效率低下,不能有针对性的有效解决问题。虽然有一些措施制定的辅助手段,如大气污染预测模型等,但这些模型基本都是用于预测空气质量,不能直接帮助决策者进行空气污染治理措施的选取,起到的辅助作用较小,而且这些模型基本都需要超级计算机运行,需要高水平的专业人员操作,耗时长、成本高、操作难度大。综上所述,目前还没有一个科学的、***的、便捷、高效的技术手段,用于相关的污染治理措施进行评估和选择。
针对背景技术中对污染治理措施选取方式的不足,本发明提供了一种高效、快速、便捷、科学的空气污染治理措施选择方法。该方法是一种根据城市空气污染特征为城市优选空气污染防治措施的动态选择方法。该方法依据城市PM2.5、O3、SO2、NO2等污染物浓度超标率和其相对应的排放来源构成(根据排放清单得出)对空气污染治理措施进行选取,通过为措施库中各措施进行评分,根据分数大小为各城市优选适应其空气质量现状的措施。
本发明所提供的空气污染治理措施选择方法可帮助决策者根据当地空气污染特征,科学的评判各项治理措施的适用性、有效性,进而有针对性的解决当地空气污染问题。该发明提供的空气污染治理措施选择方法一方面避免了依靠经验进行选择所带来的盲目性、低效性,提高了科学性、针对性;一方面与各大气污染预测模型等辅助手段相比,简单、高效、成本低、无需高级专业人员操作,可使决策者相对较快速地选择与当地实际情况相适宜的空气污染治理措施,提高措施选择效率。
本发明是一种根据城市空气污染特征为城市优选空气污染防治措施的动态选择方法。措施选取的基本依据是城市PM2.5、SO2、NO2的浓度超标率和其相对应的排放来源构成(根据排放清单得出)。该方法可通过为措施库中各措施进行评分,根据分数大小为各城市优选适应其空气质量现状的措施。
具体的选取流程包括建立措施库、分析污染物监测数据超标情况、分析污染物来源构成、计算措施权重、计算措施评分、选取措施。
需要建立当地的空气污染治理措施库并按照当地污染物排放清单排放源的分类方法,根据措施本身的适用范围进行分类,形成不同类别的措施库。措施库中每项措施有针对各项污染物的去除适用性系数;措施的污染物治理效果系数用来描述某项措施是否适用于去除某项污染物,具体值根据措施对污染物的去除效率确定。
需要分析污染物监测数据超标情况,可以根据相应的环境空气质量标准来计算污染物浓度的超标率,根据超标情况确定不同的超标系数。
需要分析污染物来源构成,根据城市各污染物的排放清单分析各污染物(如PM2.5、SO2、NO2)的排放来源构成。污染物排放清单指各种排放源在一定时间跨度和空间区域内向大气中排放的大气污染物的量的集合。根据排放清单计算得出各排放源对各污染物的贡献率。
根据各污染物的超标情况及来源构成为措施库中各措施类别分配不同权重。
根据各措施类别权重、污染物超标情况、措施的污染物去除效果系数计算各措施的评分,并根据评分按照一定的选择办法选取最优的措施。最简便的选择办法可以是评分排名的方法。
本方法还可以用于水污染、温室气体污染治理,相应的所需治理污染物为水污染的污染物或者温室气体污染物;治理措施为对应的水污染、温室气体治理措施即可。
建立措施库
措施库是所有可治理环境污染的有效措施(如:机动车加装DPF、燃煤电厂超低排放改造、清洁取暖等)的清单库。措施库中的措施按照城市、地理位置区域污染物排放清单排放源的分类方法,根据措施本身的适用范围进行分类,措施库的类别与排放清单中污染源的类别一致,形成不同类别的措施库,如工业源、道路移动源、非道路移动源、扬尘源、VOC相关源、火电厂、天然源等类别的措施。措施库中每项措施有针对各项污染物的措施的污染物治理效果系数,措施的污染物治理效果系数用来描述某项措施是否适用于去除某项污染物。措施的污染物治理效果系数可以直接由措施清单而来(如美国EPA措施清单);或者措施的污染物治理效果系数可以由措施的处理效果对应转化而来,取值范围为(0,1)。
措施的处理效果对应表
措施对污染物的实际去除效率 | 措施的污染物治理效果系数 |
0%-25% | 0 |
26%-75% | 0.5 |
76%-100% | 1 |
措施的污染物治理效果系数表格示例表
如污染物1可以为NOX,污染物2可以为PM2.5,污染物3可以为SO2;污染物编号可以根据城市的情况进行编辑选择。
分析城市空气质量监测数据超标情况
分析空气污染物超标率:根据空气污染物监测情况与空气质量标准得到各污染物超标率,空气质量标准可以为国标GB3095-2012、地方空气质量标准、国际空气质量标准等。
污染物浓度限值Sk:Sk表示空气质量标准中规定的空气污染物年均浓度限值,PM2.5、SO2、NO2的分别为浓度限值。
污染物浓度监测值Tk:Tk表示空气污染物浓度的监测值,TPM2.5、TSO2、TNO2、TO3分别为PM2.5、SO2、NO2、O3的浓度监测值(PM2.5、SO2、NO2可以为年均浓度监测值;如果选择O3,O3可以为日最大8小时滑动平均值的第90百分位数)根据环境空气质量标准计算PM2.5、O3、SO2、NO2的浓度超标率,计算方法如下。
εk为污染物年均浓度监测值与国家标准中的污染物年均浓度限值的比值,各污染物浓度限值可以由我国家标准而得,还可以由别的空气质量标准而来,如下表。
中国空气污染物浓度国家标准限制表
污染物 | 统计指标 | 国家标准 |
PM<sub>2.5</sub> | 年均浓度 | 35μg/m<sup>3</sup> |
SO<sub>2</sub> | 年均浓度 | 60μg/m<sup>3</sup> |
NO<sub>2</sub> | 年均浓度 | 40μg/m<sup>3</sup> |
O<sub>3</sub> | 日最大8小时滑动平均值的第90百分位数 | 160μg/m<sup>3</sup> |
当计算得到εk后便可以判断各污染物超标情况,不同的超标情况对应的得到的超标系数ek如下表,超标系数的对应方式可以根据当地的不同情况调整。
当PM2.5、SO2、NO2的超标情况处于以下几种情况时,还可以对其监测值进行修正:
1)PM2.5轻微超标,SO2和NO2至少有其一不超标且接近排放限值;
2)PM2.5不超标且接近排放限值,SO2和NO2至少有其一轻微超标;
3)PM2.5、SO2和NO2超标情况相同(各自的εk值相似);
修正方法如下:
分析来源构成
根据城市各污染物的排放清单分析各污染物(PM2.5、O3、SO2、NO2)的排放来源构成。污染物排放清单指各种排放源在一定时间跨度和空间区域内向大气中排放的大气污染物的量的集合。根据排放清单计算得出各排放源对各污染物的贡献率,具体计算方法如下:
Ek:污染物的排放量参数Ek表示污染物k的排放量。
措施类别权重计算
根据各污染物的超标情况及来源构成(贡献率)为措施库中各措施分配不同权重,权重计算方法如下:
进行措施初步筛选,可以选取εk大于1的污染物(即超标污染物)所对应的防治措施;
根据超标污染物超标率和其排放清单中各污染源的比重分配超标污染物涉及到的各类措施库的权重,最终将各措施所涉及措施库的权重相加,即为各措施的权重,具体计算方法方法如下。
Wi:i类措施库的权重。
i:超标污染物涉及到的措施库类别,即污染源类别。
k:超标污染物名称。
有一些措施可同时应用在多种污染源防治中,因此同一种防治措施可同时出现在多个类别的措施库中,这类措施权重的计算需要将其在不同的措施类别中的权重都考虑进来,在评分中进行计算。
结合基础数据和措施库进行措施分数计算
措施的评价根据评分公式进行,措施的评分公式如下。
Fj:措施评分
例如:
上述公式代表1类措施库中措施A的分数情况。
措施选择方式
根据上述公式计算每项措施的评分Fj,之后根据Fj选择适用的措施。具体选择办法有以下两种,即分数直接排名法与阈值筛选法。
直接排名法
将各措施按照其评分Fj大小进行排列,Fj越大,措施排名越靠前,最终选取排名前10%的措施。
例如:总共有10项措施,这些措施的评分Fj及排名情况见下表:
措施 | A | B | C | D | E | F | G | H | I | J |
F | F<sub>A</sub>=10 | F<sub>B</sub>=9 | F<sub>B</sub>=8 | F<sub>D</sub>=7 | F<sub>E</sub>=6 | F<sub>F</sub>=5 | F<sub>G</sub>=4 | F<sub>H</sub>=3 | F<sub>I</sub>=2 | F<sub>J</sub>=1 |
排名 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
根据上表,排名前10%的措施为排名第1的措施(10*10%=1),即措施A,最终选取措施A。
阈值筛选法:选取最大的评分Fj,将0.5*Fj定为阈值,最终选取评分Fj大于阈值的措施。
例如:总共有10项措施,这些措施的评分Ej见下表:
措施 | A | B | C | D | E | F | G | H | I | J |
F | F<sub>A</sub>=10 | F<sub>B</sub>=9 | F<sub>B</sub>=8 | F<sub>D</sub>=7 | F<sub>E</sub>=6 | F<sub>F</sub>=5 | F<sub>G</sub>=4 | F<sub>H</sub>=3 | F<sub>I</sub>=2 | F<sub>J</sub>=1 |
根据上表,阈值=0.5*10=5,评分Fj大于5的措施有A、B、C、D、E,则最终选取措施A、B、C、D、E。
阈值筛选法
通过直接排名法对措施进行初步筛选后,且初步筛选得出的措施数目大于1,综合考虑实际的资金预算与时间要求对措施进行二次筛选。
当优先考虑时间要求时,即要求在一定时间范围内达到一定的处理效果时:根据每项措施的实施周期选择措施,选择实施周期小于或等于所要求时间范围的措施。当优先考虑资金预算时,即需用一定额度的资金费用达到一定的处理效果时:根据每项措施的资金成本选择措施,选择资金成本小于或等于资金预算值的措施。当综合考虑时间要求与资金预算时:根据措施的实施周期与资金成本对应用基本选择方法选择出来的措施的评分Fj进行修正,修正方法如下:
F′j:修正后的措施评分
修正之后,剔除E′j值最小的措施,剩下的措施即为最终优选出的措施。
附图说明
图1:措施选择基础方案流程示意图。
图2:引入适用性条件改进的措施选择方案流程示意图。
图3:引入PM2.5源解析修正改进的措施选择方案流程示意图。
图4:结合适用条件和PM2.5源解析修改进的措施选择方案流程示意图。
图5:甲城市PM2.5源解析结果来源构成。
图6:甲城市SO2源解析结果来源构成。
图7:甲城市NO2源解析结果来源构成。
具体实施方式
实施例一
应用上述的措施选择方法为甲城市筛选出最优的污染治理措施。甲城市空气质量监测结果如下表。
甲市各空气污染物浓度监测结果
甲城市PM2.5、SO2、NO2、O3的浓度超标率计算如下:
甲城市各空气污染物源解析结果如下表。
甲城市污染物排放来源构成表
甲城市可用污染源防治措施如下表。
甲城市可用污染源防治措施表
由于SO2超标率小于1,说明该市SO2已达标,不需采取相关措施着重治理,因此在措施库中选取PM2.5、NO2所对应的防治措施。
各措施类别权重计算如下:
通过公式:
可以计算出每个措施j的评分情况如下表:
非道路移动源 | 道路移动源 | 工业源 | 火电厂 | 其他+自然源 | 扬尘源 | 措施总分 | |
措施1 | 0.6 | 0.6 | |||||
措施2 | 0.1 | 0.1 | |||||
措施3 | 0.6 | 0.6 | |||||
措施4 | 0 | 0 | |||||
措施5 | 0.1333 | 0.1333 | |||||
措施6 | 0 | 0 | |||||
措施7 | 0 | 0 | |||||
措施8 | 0.46 | 0.46 | |||||
措施9 | 4.3 | 4.3 | |||||
措施10 | 0 | 0 | |||||
措施11 | 1.2665 | 1.2665 |
计算示例:
由于措施4、6、7的主要去除物为SO2,而该实例中SO2未超标,所以以上三个措施的最终评分为0;或者由于是治理SO2的措施,SO2在本案例中没有超标,相关SO2措施不纳入考虑范围。由于措施10的四项污染物去除率不足25%,所以措施10对四项污染物的去除适应性系数均为0,也因此措施10的最终评分为0。
实施例二
应用上述的措施选择方法为甲城市筛选出最优的污染治理措施。甲城市空气质量监测结果如下表。
甲市各空气污染物浓度监测结果
甲城市PM2.5、SO2、NO2、O3的浓度超标率计算如下:
甲城市各空气污染物源解析结果如下表。
甲城市污染物排放来源构成表
甲城市可用污染源防治措施在选择中,同一种措施对有对多个污染源有治理效果,如下表。
甲城市可用污染源防治措施表
措施编号 | 措施名称 | 非道路移动源 | 道路移动源 | 工业源 | 火电厂 |
1 | 使用重组汽油(RFG)标准 | √ | √ | ||
2 | 飞机地面支持设备:电力替代燃料 | √ | |||
3 | 柴油改造 | √ | √ | ||
4 | 洗煤 | √ | |||
5 | 低NOx燃烧器 | √ | √ | ||
6 | 水泥窑:喷雾干燥器吸收器 | √ | |||
7 | ESP静电除尘器 | √ | √ |
各措施类别所占权重计算如下:
通过公式
可以计算出每个措施j,在源解析i的领域的评分情况如下表:
措施评分结果表
非道路移动源 | 道路移动源 | 工业源 | 火电厂 | 措施总分 | |
措施1 | 0 | 0 | 0 | 0 | 0 |
措施2 | 0.1 | 0 | 0 | 0 | 0.1 |
措施3 | 0.05 | 0.43 | 0 | 0 | 0.48 |
措施4 | 0 | 0 | 0 | 0 | 0 |
措施5 | 0 | 0 | 0.1333 | 0.23 | 0.3633 |
措施6 | 0 | 0 | 0 | 0 | 0 |
措施7 | 0 | 0 | 1.333 | 2.3 | 3.633 |
由于措施4、6的主要去除物为SO2,而该实例中SO2未超标,所以以上三个措施的最终评分为0;或者由于是治理SO2的措施,SO2在本案例中没有超标,相关SO2措施不纳入考虑范围。由于措施1的四项污染物去除率不足25%,所以措施1对四项污染物的去除适应性系数均为0,也因此措施1的最终评分为0。最终措施5的总得分最高,使用分数排名的方式,那么最终优选的是措施5。
Claims (11)
1.一种选择污染治理措施的方法,包含如下步骤:
1)建立污染治理措施库:所述污染治理措施库包含治理措施对各项污染物的去除效果;
2)计算污染物超标系数:所述超标系数的由环境质量标准限值和环境监测数值计算得出;
3)分析污染物来源构成:根据城市污染物的排放源清单,分析计算得出各排放源对各污染物的贡献率;
4)计算治理措施类别权重:根据各污染物的超标率和各污染物的贡献率,计算措施库类别权重;
5)计算措施评分:根据各措施类别权重、污染物超标情况、措施的污染物去除效果系数计算各措施分数;
6)对措施分数进行筛选和排序,选取最优的措施。
2.权利要求1所述的方法,其特征在于,措施库中的措施分类方法依据地理位置区域的污染物排放清单的排放源分类方法,所述分类方法的类别包括工业源、道路移动源、非道路移动源、扬尘源、VOC相关源、火电厂、天然源。
3.权利要求1所述的方法,其特征在于,污染物超标系数的由污染物的超标比例对应判断得来,对应判断方法如下:
4.权利要求1所述的方法,其特征在于,污染物超标系数还可以直接等于污染物的超标比例。
9.权利要求1所述的方法,其特征在于,所述措施分数进行筛选和排序的方式为分数直接排名法,将各措施按照其评分Fj从大到小进行排列,选取排名前一定比例的措施,所述一定比例可以为5%,10%,20%,30%,40%。
10.权利要求1所述的方法,其特征在于,所述措施分数进行筛选和排序的方式为阈值筛选法,所述阈值筛选法为,当使用直接排名法得到了多于一个措施的结果后,根据措施的特征参数进行对得分进行修正后再次进行从大到小的排名,选取排名前一定比例的措施,所述一定比例可以为5%,10%,20%,30%,40%。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210098699.4A CN114418434B (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
CN202210103938.0A CN114418442B (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2019/100493 WO2021026792A1 (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
Related Child Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210098699.4A Division CN114418434B (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
CN202210103938.0A Division CN114418442B (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111727447A true CN111727447A (zh) | 2020-09-29 |
CN111727447B CN111727447B (zh) | 2022-03-01 |
Family
ID=72563727
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210103938.0A Active CN114418442B (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
CN201980003511.XA Active CN111727447B (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
CN202210098699.4A Active CN114418434B (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210103938.0A Active CN114418442B (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210098699.4A Active CN114418434B (zh) | 2019-08-14 | 2019-08-14 | 一种选择污染物治理措施的方法 |
Country Status (3)
Country | Link |
---|---|
CN (3) | CN114418442B (zh) |
GB (2) | GB202201846D0 (zh) |
WO (1) | WO2021026792A1 (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113487178A (zh) * | 2021-07-02 | 2021-10-08 | 天津大学 | 用于企业大气污染物减排分级评价的处理方法 |
CN113487098A (zh) * | 2021-07-14 | 2021-10-08 | 清华苏州环境创新研究院 | 一种大气污染预警信息表达及展示方法 |
CN114139856A (zh) * | 2021-10-22 | 2022-03-04 | 河北先河环保科技股份有限公司 | 船舶污染物排放治理方案确定方法及终端设备 |
CN114240263A (zh) * | 2022-02-24 | 2022-03-25 | 天津天融环境科技发展有限公司 | 一种分级管理的方法、装置、电子设备及存储介质 |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112925256A (zh) * | 2019-10-12 | 2021-06-08 | 吴郁君 | 城市环境大数据综合协同管理运转平台 |
CN112801422B (zh) * | 2021-03-23 | 2024-01-23 | 南京智汇环境气象产业研究院有限公司 | 数值模型和线性规划模型相结合的大气污染协同控制方法 |
CN112988940A (zh) * | 2021-04-02 | 2021-06-18 | 中科三清科技有限公司 | 污染溯源方法及装置 |
CN113554052A (zh) * | 2021-05-10 | 2021-10-26 | 北京数汇通信息技术有限公司 | 一种城市空气质量评估诊断*** |
CN114166708A (zh) * | 2021-11-25 | 2022-03-11 | 泛测(北京)环境科技有限公司 | 颗粒物传输贡献分析方法、装置、设备及存储介质 |
CN114455715B (zh) * | 2022-03-03 | 2023-04-25 | 四川省建筑设计研究院有限公司 | 一种基于“药方式”的水体生态治理方法及*** |
CN114792057B (zh) * | 2022-05-30 | 2023-05-12 | 交通运输部规划研究院 | 一种基于ais数据的空气质量模型船舶排放清单处理方法 |
CN115358631B (zh) * | 2022-09-21 | 2023-09-05 | 张家港市艾尔环保工程有限公司 | 一种基于有害物质检测的废气定向处理方法及*** |
CN115389718B (zh) * | 2022-10-31 | 2023-01-17 | 山东公用环保集团检测运营有限公司 | 一种基于人工智能的环境检测预警方法 |
CN116432948B (zh) * | 2023-03-16 | 2024-05-10 | 中节能国祯环保科技股份有限公司 | 一种城市雨源型河流***分步治理方法 |
CN116757366B (zh) * | 2023-08-14 | 2023-11-10 | 中科三清科技有限公司 | 高污染排放污染源选址的方法、装置、介质及电子设备 |
CN117473398B (zh) * | 2023-12-26 | 2024-03-19 | 四川国蓝中天环境科技集团有限公司 | 一种基于运渣车活动的城市扬尘污染源分类方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102521494A (zh) * | 2011-12-02 | 2012-06-27 | 北京师范大学 | 一种湖泊氮磷营养盐削减达标的容量分配方法 |
CN102628852A (zh) * | 2012-03-13 | 2012-08-08 | 北京工业大学 | 基于污染物来源识别技术的大气污染源分级方法 |
CN102855401A (zh) * | 2012-09-10 | 2013-01-02 | 济南市环境保护科学研究院 | 一种水污染源优控污染物综合指标筛选的方法 |
KR20150011531A (ko) * | 2013-07-23 | 2015-02-02 | (주) 동우엔지니어링 | 공기 오염 검출 모니터링 시스템 |
CN205245460U (zh) * | 2015-12-30 | 2016-05-18 | 河南工业职业技术学院 | 一种用于智能家居的室内空气检测及净化*** |
CN106296023A (zh) * | 2016-08-18 | 2017-01-04 | 南华大学 | 基于三标度层次分析法的铀尾矿库退役环境治理效果评估方法 |
CN106339974A (zh) * | 2016-10-21 | 2017-01-18 | 常州市环境监测中心 | 一种市区大气污染物排放清单的建立方法 |
CN106649960A (zh) * | 2016-10-12 | 2017-05-10 | 环境保护部环境规划院 | 大气多污染物环境容量三维迭代计算方法 |
CN107194139A (zh) * | 2016-03-14 | 2017-09-22 | 日电(中国)有限公司 | 大气污染源分级方法及计算设备 |
WO2018070561A1 (ko) * | 2016-10-12 | 2018-04-19 | 주식회사 엔즈 | 휴대용 에어 블로워 및 이를 이용한 미세먼지 모니터링 시스템 |
CN108226046A (zh) * | 2017-12-26 | 2018-06-29 | 中国科学技术大学 | 一种机动车尾气遥测结果修正方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101706904A (zh) * | 2009-11-25 | 2010-05-12 | 国家海洋环境监测中心 | 我国近岸海域优先控制污染物名录确定方法 |
CN102830211B (zh) * | 2012-09-10 | 2015-08-12 | 济南市环境保护科学研究院 | 一种水污染源应急污染物综合指标筛选的方法 |
CN105628091B (zh) * | 2014-10-31 | 2018-05-04 | 青岛海尔空调器有限总公司 | 环境数据的检测方法和环境数据的检测*** |
CN106021898A (zh) * | 2016-05-13 | 2016-10-12 | 中国环境科学研究院 | 一种地下水污染源强分级评价方法 |
CN106324204B (zh) * | 2016-11-02 | 2018-09-28 | 扬州大学 | 城市大气pm2.5主要污染源的确定方法 |
CN106651036A (zh) * | 2016-12-26 | 2017-05-10 | 东莞理工学院 | 空气质量预报*** |
CN108896713B (zh) * | 2018-06-26 | 2021-04-23 | 湖南凯迪工程科技有限公司 | 针对气体污染源的污染治理方法及治理*** |
CN109604316B (zh) * | 2018-10-26 | 2020-11-27 | 中国辐射防护研究院 | 一种放射性污染土壤修复技术筛选方法 |
-
2019
- 2019-08-14 GB GBGB2201846.9D patent/GB202201846D0/en active Pending
- 2019-08-14 CN CN202210103938.0A patent/CN114418442B/zh active Active
- 2019-08-14 CN CN201980003511.XA patent/CN111727447B/zh active Active
- 2019-08-14 CN CN202210098699.4A patent/CN114418434B/zh active Active
- 2019-08-14 GB GB2201846.9A patent/GB2600072A/en not_active Withdrawn
- 2019-08-14 WO PCT/CN2019/100493 patent/WO2021026792A1/zh active Application Filing
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102521494A (zh) * | 2011-12-02 | 2012-06-27 | 北京师范大学 | 一种湖泊氮磷营养盐削减达标的容量分配方法 |
CN102628852A (zh) * | 2012-03-13 | 2012-08-08 | 北京工业大学 | 基于污染物来源识别技术的大气污染源分级方法 |
CN102855401A (zh) * | 2012-09-10 | 2013-01-02 | 济南市环境保护科学研究院 | 一种水污染源优控污染物综合指标筛选的方法 |
KR20150011531A (ko) * | 2013-07-23 | 2015-02-02 | (주) 동우엔지니어링 | 공기 오염 검출 모니터링 시스템 |
CN205245460U (zh) * | 2015-12-30 | 2016-05-18 | 河南工业职业技术学院 | 一种用于智能家居的室内空气检测及净化*** |
CN107194139A (zh) * | 2016-03-14 | 2017-09-22 | 日电(中国)有限公司 | 大气污染源分级方法及计算设备 |
CN106296023A (zh) * | 2016-08-18 | 2017-01-04 | 南华大学 | 基于三标度层次分析法的铀尾矿库退役环境治理效果评估方法 |
CN106649960A (zh) * | 2016-10-12 | 2017-05-10 | 环境保护部环境规划院 | 大气多污染物环境容量三维迭代计算方法 |
WO2018070561A1 (ko) * | 2016-10-12 | 2018-04-19 | 주식회사 엔즈 | 휴대용 에어 블로워 및 이를 이용한 미세먼지 모니터링 시스템 |
CN106339974A (zh) * | 2016-10-21 | 2017-01-18 | 常州市环境监测中心 | 一种市区大气污染物排放清单的建立方法 |
CN108226046A (zh) * | 2017-12-26 | 2018-06-29 | 中国科学技术大学 | 一种机动车尾气遥测结果修正方法 |
Non-Patent Citations (2)
Title |
---|
刘成伦 等: ""重庆市机动车尾气对大气环境的影响分析及减缓措施"", 《环境污染与防治》 * |
武娟妮 等: ""中国工业大气污染物排放的影响因素研究"", 《环境污染与防治》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113487178A (zh) * | 2021-07-02 | 2021-10-08 | 天津大学 | 用于企业大气污染物减排分级评价的处理方法 |
CN113487098A (zh) * | 2021-07-14 | 2021-10-08 | 清华苏州环境创新研究院 | 一种大气污染预警信息表达及展示方法 |
CN114139856A (zh) * | 2021-10-22 | 2022-03-04 | 河北先河环保科技股份有限公司 | 船舶污染物排放治理方案确定方法及终端设备 |
CN114240263A (zh) * | 2022-02-24 | 2022-03-25 | 天津天融环境科技发展有限公司 | 一种分级管理的方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN114418442A (zh) | 2022-04-29 |
CN114418442B (zh) | 2022-09-27 |
WO2021026792A1 (zh) | 2021-02-18 |
GB2600072A (en) | 2022-04-20 |
CN111727447B (zh) | 2022-03-01 |
GB202201846D0 (en) | 2022-03-30 |
CN114418434B (zh) | 2022-11-18 |
CN114418434A (zh) | 2022-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111727447B (zh) | 一种选择污染物治理措施的方法 | |
CN107871210B (zh) | 一种大气环境容量核算方法 | |
Zhu et al. | Daily air quality index forecasting with hybrid models: A case in China | |
Ying et al. | Source contributions to the regional distribution of secondary particulate matter in California | |
Duque et al. | Evaluating strategies to reduce urban air pollution | |
Memmesheimer et al. | Long-term simulations of particulate matter in Europe on different scales using sequential nesting of a regional model | |
Iodice et al. | Appraisal of pollutant emissions and air quality state in a critical I talian region: Methods and results | |
Parra et al. | Development of the high spatial resolution EMICAT2000 emission model for air pollutants from the north-eastern Iberian Peninsula (Catalonia, Spain) | |
CN116128421B (zh) | 一种耦合污染源排放和解析的大气污染管控方案生成方法 | |
KR20210086326A (ko) | 지역별 미세먼지 농도 예측방법 및 시스템 | |
Volta et al. | GAMES, a comprehensive gas aerosol modelling evaluation system | |
Cheng et al. | A new statistical modeling and optimization framework for establishing high-resolution PM10 emission inventory–I. Stepwise regression model development and application | |
Wang et al. | Capturing uncertainty in emission estimates related to vehicle electrification and implications for metropolitan greenhouse gas emission inventories | |
Goyal et al. | Development of strategic air quality improvement framework for urban hotspots | |
Iodice et al. | Air pollution and air quality state in an Italian national interest priority site. part 2: the pollutant dispersion | |
Iodice et al. | Environmental Assessment of a Wide Area Under Surveillance with Different Air Pollution Sources. | |
Liu et al. | Comprehensive evaluation of flue gas desulfurization and denitrification technologies of six typical enterprises in Chengdu, China | |
Bui et al. | Spatio-temporal modelling and related factors analysis of PM2. 5 pollution: A case study of Ho Chi Minh city, Vietnam | |
Chen et al. | A hybrid approach to forecast air quality during high-PM concentration pollution period | |
Dzhidzhelava et al. | The environmental component of construction projects in the aspect of sustainable development | |
Samaali et al. | A new tool for processing atmospheric emission inventories: Technical aspects and application to the ESCOMPTE study area | |
Sowden et al. | Developing a spatially and temporally resolved emission inventory for photochemical modeling in the City of Cape Town and assessing its uncertainty | |
Thunis et al. | Air pollution and emission reductions over the Po-valley: Air Quality Modelling and Integrated Assessment | |
Wang | Long-term Greenhouse Gas Emissions Mitigation in California and the Associated Regional Air Quality and Public Health Impacts | |
Ansari | Mitigation and Control of Urban Air Pollution in Beijing |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Xie Hongxing Inventor after: Zhang Yunpeng Inventor after: Wang Kan Inventor after: He Xin Inventor after: Tokitaka Inventor before: Xie Hongxing Inventor before: He Xin Inventor before: Tokitaka |
|
GR01 | Patent grant | ||
GR01 | Patent grant |