CN109342658A - 基于wrf模型的计算机程序化计算大气环境容量的方法 - Google Patents

基于wrf模型的计算机程序化计算大气环境容量的方法 Download PDF

Info

Publication number
CN109342658A
CN109342658A CN201811274304.1A CN201811274304A CN109342658A CN 109342658 A CN109342658 A CN 109342658A CN 201811274304 A CN201811274304 A CN 201811274304A CN 109342658 A CN109342658 A CN 109342658A
Authority
CN
China
Prior art keywords
formula
follows
east longitude
north latitude
atmospheric environment
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
Application number
CN201811274304.1A
Other languages
English (en)
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.)
Henan Provincial Ecological Environment Technology Center
Original Assignee
Henan Research Academy Of Environment Protection Sciences
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 Henan Research Academy Of Environment Protection Sciences filed Critical Henan Research Academy Of Environment Protection Sciences
Priority to CN201811274304.1A priority Critical patent/CN109342658A/zh
Publication of CN109342658A publication Critical patent/CN109342658A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004—Gaseous mixtures, e.g. polluted air

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medicinal Chemistry (AREA)
  • Food Science & Technology (AREA)
  • Combustion & Propulsion (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于WRF模型的计算机程序化计算大气环境容量的方法，本发明有效解决了在计算大气环境容量时采用人工估算，公式中相关参数采用较大地理范围内年平均值和经验数值所带来的计算结果时空精度差、经验依赖性大的问题；解决的技术方案包括：选定某一区域对其进行网格化处理通过查询相关文件计并结合WRF模型计算出该区域的大气稳定度等级、大气混合层厚度，经计算公式得出污染物的大气环境容量；WRF模型可利用网上公开的基础数据经计算输出精确到小时的各种气象参数，且计算地理范围可精确到数平方公里，通过编制计算机程序从WRF计算结果文件中提取计算大气环境容量所需的参数，可实现任意范围内时间精确到小时的大气环境容量。

Description

基于WRF模型的计算机程序化计算大气环境容量的方法
技术领域
本发明涉及环境影响评价技术领域，尤其是涉及一种基于WRF模型的计算机程序化计算大气环境容量的方法。
背景技术
环境容量(environment capacity)又称环境负载容量、地球环境承载容量或负荷量，是在人类生存和自然生态系统不致受害的前提下，某一环境所能容纳的污染物的最大负荷量。大气环境容量主要与大气混合层厚度、风速、沉降速率、降水、污染物标准浓度等因素相关，以往的大气环境容量研究和实际应用中大气环境容量计算主要采用A值法，依据国家颁布的相关标准和及再次基础上进行改良的公式进行计算，但由于采用人工估算，公式中相关参数的取值多采用较大地理范围内的年平均值和经验数值，因此计算结果存在时空精度差、经验依赖性较大的问题。本发明以WRF模型为例，该模型可利用在互联网上公开的基础数据经计算输出精确到小时的各种气象参数，且计算地理范围可精确到数平方公里，通过编制计算机程序实现自动化计算，最终可实现任意地理范围内时间精确到小时的大气环境容量计算；因此，本发明提供一种基于WRF模型的计算机程序化计算大气环境容量的方法。
发明内容
针对上述情况，为克服现有技术之缺陷，本发明提供一种基于WRF模型的计算机程序化计算大气环境容量的方法。有效的解决了在计算大气环境容量时采用人工估算，公式中相关参数采用较大地理范围内年平均值和经验数值所带来的计算结果时空精度差、经验依赖性大的问题。
基于WRF模型的计算机程序化计算大气环境容量的方法，其特征在于，计算步骤如下：
步骤一、计算某一选定区域W污染物的大气环境容量，选定所要计算区域中心坐标为北纬A，东经B，边长为X*Y的范围，对于所要计算区域进行网格化处理，将其格网化为x*y个网格，选取计算范围以单个网格为例进行计算，计算时间长度为K年L月M日0时-23时；
根据公式
δ＝[0.006918-0.399912cosθ0+0.070257sinθ0-0.006758cos2θ0+0.000907sin2θ0-0.002697cos3θ0+0.001480sin3θ0]×180/π
式中：θ0＝360dn/365,deg；
δ—太阳倾角，deg；
dn—一年中日期序数，0，1，2……，364
得出所要计算区域K年L月M日的太阳倾角δ
步骤二、将步骤一中求得的太阳倾角代入公式
h0＝arcsin[sinφsinδ+cosφcosδcos(15t+λ-300)]
式中：h0—太阳高度角，deg；
φ—当地纬度，deg；
t—北京时间，h；
λ—当地纬度，deg；
计算网格点中心坐标为：北纬A，东经B，得出0时-23时各个时间段内对应的太阳高度角h0，根据太阳高度角可知：K年L月M日昼间为c-d时，夜间为g-h时，e-f时；
步骤三、查询地面气象sam文件得出网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的总云量/低云量；
步骤四、由太阳高度角h0、总云量/低云量、以及昼夜情况查询太阳辐射等级表得出所要计算网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的太阳辐射等级；
步骤五、将以上数据参数由WRF模型运算完成输出所要计算网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内的风速U、干沉降速率Ud、降水速率R；
步骤六、由太阳辐射等级与风速的对应关系，查询大气稳定度等级表得出网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内的大气稳定度等级；
步骤七、结合地区序号表确定网格点中心坐标为：北纬A，东经B，所处的地区序号，并且结合大气稳定度等级查询我国不同地区和在不同大气稳定度等级下所对应的as/bs数值，求得网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的as/bs数值；
步骤八、综合网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的风速以及0时-23时各个时间段内对应的as/bs数值，根据公式：
在大气稳定度为A、B、C和D级时：
在大气稳定度为E和F级时：
式中：H—混合层厚度，m；
U—10m高度上平均风速，m/s；大于6m/s时取为6m/s；
as，bs—混合层系数；
f—地转参数；
Ω—地转角速度，由《制定地方大气污染物排放标准的技术方法》GB/T13201-91里规定取为7.29×10-5rad/s；
—地理纬度；
求得网格点中心坐标为：北纬A，东经B，地转参数f以及0时-23时各个时间段内对应的大气混合层厚度；
步骤九、大气环境容量Q的表达式表示为：
其中：H:大气混合层厚度(m)
U：风速(m/s)
Ud:干沉降速度(m/s)
R：降水速率(mm/a)
ωr：清洗比，由《制定地方大气污染物排放标准的技术方法》GB/T13201-91里规定取值为1.9×10-5
Cs：污染物标准浓度(mg/m3)
S：区域面积(km2)
Q：大气污染物容量(104t/a)
区域内存在不同大气功能区时各功能区Ai计算公式为：
Ai＝ACi
式中Ci为大气环境质量标准规定的第i功能区类别对应的年日平均浓度限值(mg/m3)
Qi＝Ai×Csi×Si/S1/2
式中：
Csi：大气环境质量标准规定的第i功能区类别对应的年平均浓度限值(mg/m3)
Si：第i功能区面积
Ci、Csi由查询气象文件及相关国家标准获得；
计算得出网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的Qi值，由于公式中Qi值量纲为104t/a，乘以系数1.14最终得出小时值Qhi对Qhi进行求和得出网格点中心坐标为：北纬A，东经B，0时-23时内W污染物的大气环境容量。
(1)本发明以WRF模型为基础，由于WRF模型可以利用互联网上公开的基础数据经计算输出精确到小时的各种气象参数，且计算地理范围可精确到数平方公里，通过编制计算机程序从WRF计算结果文件中提取计算大气环境容量所需的参数实现自动化计算，最终可实现任意地理范围内时间精确到小时的大气环境容量计算；
(2)本发明可以快速准确的计算所关注区域的大气环境容量，为企业生产和大气环境相关政策的制定提供科学依据。
附图说明
图1为本发明所选计算区域示意图。
具体实施方式
有关本发明的前述及其他技术内容、特点与功效，在以下配合参考附图至图对实施例的详细说明中，将可清楚的呈现。以下实施例中所提到的结构内容，均是以说明书附图为参考。
下面将参照附图描述本发明的各示例性的实施例。
实施例1，基于WRF模型的计算机程序化计算大气环境容量的方法，其特征在于，计算步骤如下：
步骤一、计算某一选定区域SO2污染物的大气环境容量，选定所要计算区域中心坐标为北纬33.3°，东经113.5°，边长为50km×50km的范围，对于所要计算区域进行网格化处理，将其格网化为50×50个网格，选取计算范围以单个网格为例进行计算，计算时间长度为2012年1月1日0时-23时；
计算区域示意图为选取计算网格
根据公式
δ＝[0.006918-0.399912cosθ0+0.070257sinθ0-0.006758cos2θ0+0.000907sin2θ0-0.002697cos3θ0+0.001480sin3θ0]×180/π
式中：θ0＝360dn/365,deg；
δ—太阳倾角，deg；
dn—一年中日期序数，0，1，2……，364
dn＝0,θ0＝0，计算得出δ＝-23.06
步骤二、将步骤一中求得的太阳倾角代入公式
h0＝arcsin[sinφsinδ+cosφcosδcos(15t+λ-300)]
式中：h0—太阳高度角，deg；
φ—当地纬度，deg；
t—北京时间，h；
λ—当地纬度，deg；
计算网格点中心坐标为：北纬33.3°，东经113.5°，得出0时-23时各个时间段内对应的太阳高度角h0
t 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
h<sub>0</sub> -78 -77 -67 -55 -42 -30 -18 -6 5 15 24 30 33 33 30 23 14
t 17 18 19 20 21 22 23
h<sub>0</sub> 4 -7 -19 -31 -44 -56 -68
根据太阳高度角可知：2012年1月1日昼间为8-17时，夜间为0-7时，18-23时；
步骤三、查询地面气象sam文件得出网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内对应的总云量/低云量
t 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
总云 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
低云 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
t 17 18 19 20 21 22 23
总云 5 5 5 5 5 5 5
低云 5 5 5 5 5 5 5
注：云量(全天空十分制)观测规则见中央气象局编定的《地面气象观测规范》第3.3节。
步骤四、由太阳高度角h0、总云量/低云量、以及昼夜情况查询太阳辐射等级表得出所要计算网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内对应的太阳辐射等级；
太阳辐射等级
网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内对应的太阳辐射等级：
t 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
b<sub>1</sub> 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
t 17 18 19 20 21 22 23
b<sub>1</sub> 0 0 0 0 0 0 0
步骤五、将以上数据参数由WRF模型运算完成输出所要计算网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内的风速U、干沉降速率Ud、降水速率R
T 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
U 2.4 3.0 3.5 3.6 0.5 4.5 4.9 5.0 4.6 4.8 5.3 5.2 0.5 4.9 4.3 3.5 2.5
T 17 18 19 20 21 22 23
U 1.8 1.4 1.0 0.7 0.4 0.3 0.6
T 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
R 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
T 17 18 19 20 21 22 23
R 0 0 0 0 0 0 0
步骤六、由太阳辐射等级与风速的对应关系，查询大气稳定度等级表得出网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内的大气稳定度等级
大气稳定度等级表
注：1)地面风速(m/s)系指离地面10m高度处10分钟平均风速
网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内的大气稳定度等级：
t 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
b<sub>2</sub> D D D D D D D D D D D D D D D D D
t 17 18 19 20 21 22 23
b<sub>2</sub> D D D D D D D
步骤七、结合地区序号表确定网格点中心坐标为：北纬33.3°，东经113.5°，所处的地区序号，并且结合大气稳定度等级查询我国不同地区和在不同大气稳定度等级下所对应的as/bs数值，求得网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内对应的as/bs数值；
由表1确定地区序号为“3”
表1
结合地区序号和大气稳定度等级b2，由表E1查出as,bs值：
表E1我国各地区as和bs值
网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内对应的as/bs数值：
T 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
a<sub>s</sub> 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019 0.019
T 17 18 19 20 21 22 23
a<sub>s</sub> 0.019 0.019 0.019 0.019 0.019 0.019 0.019
步骤八、综合网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内对应的风速以及0时-23时各个时间段内对应的as/bs数值，根据公式：
在大气稳定度为A、B、C和D级时：
在大气稳定度为E和F级时：
式中：H—混合层厚度，m；
U—10m高度上平均风速，m/s；大于6m/s时取为6m/s；
as，bs—混合层系数，见表E1；
f—地转参数；
Ω—地转角速度，由《制定地方大气污染物排放标准的技术方法》GB/T13201-91里规定取值为7.29×10-5rad/s；
—地理纬度；
得出网格点中心坐标为：北纬33.3°，东经113.5°的地转参数为f＝7.96×10-5；
求得网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内对应的大气混合层厚度：
t 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
H 573 716 835 859 119 1074 1170 1193 1098 1146 1265 1241 119 1170 1026 835 597
t 17 18 19 20 21 22 23
H 430 334 239 167 95 72 143
综合以上结果，A值计算公式中所需参数已全部求得，计算所需的各参数分别为：
步骤九、大气环境容量Q的表达式表示为：
其中：H:大气混合层厚度(m)
U：风速(m/s)
Ud:干沉降速度(m/s)
R：降水速率(mm/a)
ωr：清洗比，由《制定地方大气污染物排放标准的技术方法》GB/T13201-91里规定取值为1.9×10-5；
Cs：污染物标准浓度(mg/m3)
S：区域面积(km2)
Q：大气污染物容量(104t/a)
由此计算出2012年1月1日0-23时逐时A值：
t 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
A 3.84 5.99 8.16 8.64 0.17 13.49 15.99 16.65 14.10 15.35 18.71 18.01 0.17 15.99 12.31 8.16 4.16
t 17 18 19 20 21 22 23
A 2.16 1.31 0.67 0.33 0.11 0.06 0.24
区域内存在不同大气功能区时各功能区Ai计算公式为：
Ai＝ACi
式中Ci为大气环境质量标准规定的第i功能区类别对应的年日平均浓度限值(mg/m3)
Qi＝Ai×Csi×Si/S1/2
式中：
Csi：大气环境质量标准规定的第i功能区类别对应的年平均浓度限值(mg/m3)
Si：第i功能区面积
Ci、Csi由查询气象文件及相关国家标准获得；
根据环境空气质量标准(GB3095-2012)中相关规定：
环境空气质量功能区的分类和标准分级
环境空气质量功能区分类：
一类区为自然保护区、风景名胜区和其它需要特殊保护的地区；
二类区为居住区、商业交通居民混合区、文化区、工业区和农村地区；
环境空气质量标准分级：
一类区执行一级浓度限值
二类区执行二级浓度限值
各项污染物的浓度限值：
设定计算网格点执行二级标准Csi＝0.06mg/m3Ci＝0.15mg/m3
由公式
Qi＝Ai×Csi×Si/S1/2
Ai＝ACi
由此计算出2012年1月1日0-23时逐时Ai值：
t 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
A<sub>i</sub> 0.58 0.90 1.22 1.30 0.03 2.02 2.40 2.50 2.11 2.30 2.81 2.70 0.03 2.40 1.85 1.22 0.62
t 17 18 19 20 21 22 23
A<sub>i</sub> 0.32 0.20 0.10 0.05 0.02 0.009 0.04
若出现U＝0公式中分母为0的情况，该时间点计算结果无效舍弃
计算得出网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时各个时间段内对应的Qi值，由于公式中Qi值量纲为104t/a(104吨/年)，乘以系数1.14最终得出小时值Qhi对Qhi进行求和得出网格点中心坐标为：北纬33.3°，东经113.5°，0时-23时内W污染物的大气环境容量；
1.14为万吨/年与吨/小时的换算系数，由10000(吨)/8760(小时)计算得到；
求解Qi值，由于公式中Qi值量纲为104t/a，乘以系数1.14最终得出小时值Qhi(t/h)
对于Qhi进行求和，得出2012年1月1日一天时间内选定计算网格内SO2的环境容量为0.038t。
本发明以WRF模型为基础，由于WRF模型可以利用互联网上公开的基础数据经计算输出精确到小时的各种气象参数，且计算地理范围可精确到数平方公里，通过编制计算机程序从WRF计算结果文件中提取计算大气环境容量所需的参数实现自动化计算，最终可实现任意地理范围内时间精确到小时的大气环境容量计算；本发明可以快速准确的计算所关注区域的大气环境容量，为企业生产和大气环境相关政策的制定提供科学依据。
以上所述的实施例并非对本发明的范围进行限定，在不脱离本发明设计构思的前提下，本领域所述技术人员对本发明的技术方案作出的各种变形和改进，均应纳入本发明的权利要求书确定的保护范围内。

Claims (1)

1.基于WRF模型的计算机程序化计算大气环境容量的方法，其特征在于，计算步骤如下：
步骤一、计算某一选定区域W污染物的大气环境容量，选定所要计算区域中心坐标为北纬A，东经B，边长为X*Y的范围，对于所要计算区域进行网格化处理，将其格网化为x*y个网格，选取计算范围以单个网格为例进行计算，计算时间长度为K年L月M日0时-23时；
根据公式
δ＝[0.006918-0.399912cosθ0+0.070257sinθ0-0.006758cos2θ0+0.000907sin2θ0-0.002697cos3θ0+0.001480sin3θ0]×180/π
式中：θ0＝360dn/365,deg；
δ—太阳倾角，deg；
dn—一年中日期序数，0，1，2……，364
得出所要计算区域K年L月M日的太阳倾角δ
步骤二、将步骤一中求得的太阳倾角代入公式
h0＝arcsin[sinφsinδ+cosφcosδcos(15t+λ-300)]
式中：h0—太阳高度角，deg；
φ—当地纬度，deg；
t—北京时间，h；
λ—当地纬度，deg；
计算网格点中心坐标为：北纬A，东经B，得出0时-23时各个时间段内对应的太阳高度角h0，根据太阳高度角可知：K年L月M日昼间为c-d时，夜间为g-h时，e-f时；
步骤三、查询地面气象sam文件得出网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的总云量/低云量；
步骤四、由太阳高度角h0、总云量/低云量、以及昼夜情况查询太阳辐射等级表得出所要计算网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的太阳辐射等级；
步骤五、将以上数据参数由WRF模型运算完成输出所要计算网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内的风速U、干沉降速率Ud、降水速率R；
步骤六、由太阳辐射等级与风速的对应关系，查询大气稳定度等级表得出网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内的大气稳定度等级；
步骤七、结合地区序号表确定网格点中心坐标为：北纬A，东经B，所处的地区序号，并且结合大气稳定度等级查询我国不同地区和在不同大气稳定度等级下所对应的as/bs数值，求得网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的as/bs数值；
步骤八、综合网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的风速以及0时-23时各个时间段内对应的as/bs数值，根据公式：
在大气稳定度为A、B、C和D级时：
在大气稳定度为E和F级时：
式中：H—混合层厚度，m；
U—10m高度上平均风速，m/s；大于6m/s时取为6m/s；
as，bs—混合层系数；
f—地转参数；
Ω—地转角速度，由《制定地方大气污染物排放标准的技术方法》GB/T13201-91里规定取为7.29×10-5rad/s；
—地理纬度；
求得网格点中心坐标为：北纬A，东经B，地转参数f以及0时-23时各个时间段内对应的大气混合层厚度；
步骤九、大气环境容量Q的表达式表示为：
其中：H:大气混合层厚度(m)
U：风速(m/s)
Ud:干沉降速度(m/s)
R：降水速率(mm/a)
ωr：清洗比，由《制定地方大气污染物排放标准的技术方法》GB/T13201-91里规定取为(无量纲取值1.9×10-5)
Cs：污染物标准浓度(mg/m3)
S：区域面积(km2)
Q：大气污染物容量(104t/a)
区域内存在不同大气功能区时各功能区Ai计算公式为：
Ai＝ACi
式中Ci为大气环境质量标准规定的第i功能区类别对应的年日平均浓度限值(mg/m3)
Qi＝Ai×Csi×Si/S1/2
式中：
Csi：大气环境质量标准规定的第i功能区类别对应的年平均浓度限值(mg/m3)
Si：第i功能区面积
Ci、Csi由查询气象文件及相关国家标准获得；
计算得出网格点中心坐标为：北纬A，东经B，0时-23时各个时间段内对应的Qi值，由于公式中Qi值量纲为104t/a，乘以系数1.14最终得出小时值Qhi对Qhi进行求和得出网格点中心坐标为：北纬A，东经B，0时-23时内W污染物的大气环境容量。
CN201811274304.1A 2018-10-30 2018-10-30 基于wrf模型的计算机程序化计算大气环境容量的方法 Pending CN109342658A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811274304.1A CN109342658A (zh) 2018-10-30 2018-10-30 基于wrf模型的计算机程序化计算大气环境容量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811274304.1A CN109342658A (zh) 2018-10-30 2018-10-30 基于wrf模型的计算机程序化计算大气环境容量的方法

Publications (1)

Publication Number Publication Date
CN109342658A true CN109342658A (zh) 2019-02-15

Family

ID=65311117

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811274304.1A Pending CN109342658A (zh) 2018-10-30 2018-10-30 基于wrf模型的计算机程序化计算大气环境容量的方法

Country Status (1)

Country Link
CN (1) CN109342658A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109783980A (zh) * 2019-02-16 2019-05-21 河南省环境保护科学研究院 基于大气环境容量、环境质量及污染排放量制定大气污染物减排方案的方法
CN114417579A (zh) * 2021-12-31 2022-04-29 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种对流层电波环境数值模式结果的不同单位网格距转换方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105115859A (zh) * 2015-07-29 2015-12-02 北京千安哲信息技术有限公司 用于监测泄漏气体扩散浓度的方法及装置
CN106529746A (zh) * 2016-12-29 2017-03-22 南京恩瑞特实业有限公司 基于动力和热力因子的空气质量动态融合统计预报方法
CN106649960A (zh) * 2016-10-12 2017-05-10 环境保护部环境规划院 大气多污染物环境容量三维迭代计算方法
CN107871210A (zh) * 2017-11-03 2018-04-03 南开大学 一种大气环境容量核算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105115859A (zh) * 2015-07-29 2015-12-02 北京千安哲信息技术有限公司 用于监测泄漏气体扩散浓度的方法及装置
CN106649960A (zh) * 2016-10-12 2017-05-10 环境保护部环境规划院 大气多污染物环境容量三维迭代计算方法
CN106529746A (zh) * 2016-12-29 2017-03-22 南京恩瑞特实业有限公司 基于动力和热力因子的空气质量动态融合统计预报方法
CN107871210A (zh) * 2017-11-03 2018-04-03 南开大学 一种大气环境容量核算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张志 等: "郑州市大气污染原因及防治策略研究", 《环境与可持续发展》 *
环境保护部环境工程评估中心: "《环境影响评价技术方法》", 2 February 2011, 中国环境科学出版社 *
胡艺文: "黑龙江省大气PM2.5环境容量测算研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109783980A (zh) * 2019-02-16 2019-05-21 河南省环境保护科学研究院 基于大气环境容量、环境质量及污染排放量制定大气污染物减排方案的方法
CN114417579A (zh) * 2021-12-31 2022-04-29 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种对流层电波环境数值模式结果的不同单位网格距转换方法
CN114417579B (zh) * 2021-12-31 2023-01-24 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种对流层电波环境数值模式结果的不同单位网格距转换方法

Similar Documents

Publication Publication Date Title
Yu et al. Origin of air pollution during a weekly heavy haze episode in Hangzhou, China
Prospero Saharan dust transport over the North Atlantic Ocean and Mediterranean: An overview
Pillai et al. A study of PM, PM10 and PM2. 5 concentration at a tropical coastal station
Yu Water spray geoengineering to clean air pollution for mitigating haze in China’s cities
Ahmed et al. Geomorphological characteristics of the Um-Rimam depression in northern Kuwait
Green et al. Seasonal and interplanetary magnetic field dependence of the field-aligned currents for both Northern and Southern Hemispheres
Munir et al. Spatiotemporal analysis of fine particulate matter (PM2. 5) in Saudi Arabia using remote sensing data
CN109342658A (zh) 基于wrf模型的计算机程序化计算大气环境容量的方法
Itahashi et al. Long-term trends of satellite-based fine-mode aerosol optical depth over the Seto Inland Sea, Japan, over two decades (2001–2020)
Lenes et al. A model of Saharan dust deposition to the eastern Gulf of Mexico
Claiborn et al. Regional measurements and modeling of windblown agricultural dust: The Columbia Plateau PM10 Program
Yu et al. A climatological study of air pollution potential in China
Muter et al. Analysis of seasonal and annual relative humidity using GIS for selected stations over Iraq during the period (1980-2017)
Jin et al. Observation and simulation of near-surface wind and its variation with topography in Urumqi, West China
Kim et al. Spatio-temporal characteristics of wind observations over South Korea: 1982–2011
Eldamaty et al. GIS-Based Multi Criteria Analysis for Solar Power Plant Site Selection Support in Mecca
Renné et al. Solar resource assessment for Sri Lanka and Maldives
Kim et al. Atmospheric deposition of nitrogen and sulfur in the Yellow Sea region: Significance of long-range transport in East Asia
CN106203891A (zh) 基于多源遥感影像数据的国土开发密度三维综合评价方法
Teng et al. Study on the Temporal and Spatial Variation of PM 2.5 in Eight Main Cities of Yunnan Province
Sarr et al. Impact of temporal and spatial variability of solar resource on technical sizing of isolated solar installations in Senegal using satellite data
Tadić et al. Erosion intensity assessment using erosion potential method and geographic information systems: A case study of Beočin municipality, Serbia
Fang et al. A climatic environmental performance assessment method for ecological city construction: Application to Beijing Yanqi Lake
CN107766799B (zh) 基于尺度效应的多源遥感数据源汇景观的分析方法及系统
Xi et al. Spatiotemporal Evolution Characteristics of Ecosystem Service Values Based on NDVI Changes in Island Cities

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
CB02 Change of applicant information

Address after: 450000 No.1, Shunhe Road, Jinshui District, Zhengzhou City, Henan Province

Applicant after: Henan Academy of ecological and Environmental Sciences

Address before: 450000 No.1, Shunhe Road, Jinshui District, Zhengzhou City, Henan Province

Applicant before: HENAN ENVIRONMENTAL PROTECTION Research Institute

CB02 Change of applicant information
TA01 Transfer of patent application right

Effective date of registration: 20231207

Address after: 450000, No.10 Xueli Road, Boxue Road Street, Zhengdong New District, Zhengzhou City, Henan Province

Applicant after: Henan Provincial Ecological Environment Technology Center

Address before: 450000 No.1, Shunhe Road, Jinshui District, Zhengzhou City, Henan Province

Applicant before: Henan Academy of ecological and Environmental Sciences

TA01 Transfer of patent application right