原文传递
基于WRF模型的计算机程序化计算大气环境容量的方法
| 专利名称： | 基于WRF模型的计算机程序化计算大气环境容量的方法 |
| 摘要： | 本发明涉及一种基于WRF模型的计算机程序化计算大气环境容量的方法，本发明有效解决了在计算大气环境容量时采用人工估算，公式中相关参数采用较大地理范围内年平均值和经验数值所带来的计算结果时空精度差、经验依赖性大的问题；解决的技术方案包括：选定某一区域对其进行网格化处理通过查询相关文件计并结合WRF模型计算出该区域的大气稳定度等级、大气混合层厚度，经计算公式得出污染物的大气环境容量；WRF模型可利用网上公开的基础数据经计算输出精确到小时的各种气象参数，且计算地理范围可精确到数平方公里，通过编制计算机程序从WRF计算结果文件中提取计算大气环境容量所需的参数，可实现任意范围内时间精确到小时的大气环境容量。 |
| 专利类型： | 发明专利 |
| 国家地区组织代码： | 河南;41 |
| 申请人： | 河南省环境保护科学研究院 |
| 发明人： | 袁彩凤;时翔明;肖军仓;李祥华;张志;张清敏;王凯丽;张晓果;王婧 |
| 专利状态： | 有效 |
| 申请号： | CN201811274304.1 |
| 公开号： | CN109342658A |
| 代理机构： | 北京卓恒知识产权代理事务所(特殊普通合伙) 11394 |
| 代理人： | 轩文君 |
| 分类号： | G01N33/00(2006.01)I;G;G01;G01N;G01N33 |
| 申请人地址： | 450000 河南省郑州市金水区顺河路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级时： |
| 所属类别： | 发明专利 |
检索历史
应用推荐
