空气污染控制软件：AERMOD二次开发all.docx

PAGE1

PAGE1

AERMOD二次开发入门

1.AERMOD概述

AERMOD是一种广泛使用的空气污染扩散模型，由美国环保署（EPA）开发和维护。该模型主要用于评估点源、线源、面源和体源在不同气象条件下的污染物扩散情况。AERMOD基于Stable、Unstable和Neutral大气稳定度条件下的高斯扩散模型，能够模拟复杂地形和高层大气中的污染物传输。

2.AERMOD的基本输入文件

在进行AERMOD二次开发之前，了解其基本输入文件结构是非常重要的。AERMOD的输入文件主要包括以下几个部分：

源参数文件（.SRT）：定义排放源的位置、高度、排放速率等信息。

气象文件（.PRM）：包含气象参数，如风速、风向、大气稳定度等。

地形文件（.TDX）：描述地形信息，如地表高度、地表粗糙度等。

受体文件（.REC）：定义受体点的位置，用于计算污染物浓度。

2.1源参数文件（.SRT）

源参数文件定义了所有排放源的基本信息。文件格式如下：

SOURCEID,X,Y,Z,HR,MIN,SEC,Q,D,U,V,W

SOURCEID：源的标识符。

X,Y,Z：源的三维坐标（单位：米）。

HR,MIN,SEC：源的时间信息（小时、分钟、秒）。

Q：排放速率（单位：克/秒）。

D,U,V,W：源的直径和速度分量（单位：米，米/秒）。

示例代码：

#源参数文件示例

1,0,0,100,0,0,0,1000,2,0,0,0

2,1000,0,150,0,0,0,500,1,0,0,0

3.AERMOD的输出文件

AERMOD的输出文件通常包含污染物浓度的计算结果，可以是文本文件或图形文件。主要输出文件包括：

浓度文件（.OUT）：包含受体点的污染物浓度数据。

调试文件（.DBG）：包含模型运行的详细信息，用于调试和验证。

3.1浓度文件（.OUT）

浓度文件记录了每个受体点的污染物浓度。文件格式如下：

RECEPTOR,HOUR,MIN,SEC,CONC

RECEPTOR：受体点的标识符。

HOUR,MIN,SEC：时间信息（小时、分钟、秒）。

CONC：污染物浓度（单位：微克/立方米）。

示例代码：

#浓度文件示例

1,0,0,0,50.2

2,0,0,0,30.1

1,1,0,0,45.8

2,1,0,0,28.5

4.AERMOD二次开发的准备工作

进行AERMOD二次开发前，需要做好以下准备工作：

安装AERMOD：确保AERMOD软件已经正确安装在开发环境中。

熟悉AERMOD的API：了解AERMOD提供的API接口，以便进行编程调用。

准备输入数据：收集和整理所需的气象、地形、源参数和受体数据。

选择开发工具：根据个人偏好选择开发工具，如Python、MATLAB等。

4.1安装AERMOD

AERMOD的安装通常包括以下几个步骤：

下载安装包：从EPA官方网站下载最新版本的AERMOD安装包。

解压安装包：将下载的安装包解压到指定目录。

配置环境变量：将AERMOD的可执行文件路径添加到系统环境变量中。

示例代码：

#解压安装包

tar-xvfAERMOD_v19190.tar.gz

#配置环境变量（Linux示例）

exportPATH=$PATH:/path/to/AERMOD

5.使用Python进行AERMOD二次开发

Python是一种广泛使用的编程语言，具有丰富的库支持，非常适合进行AERMOD的二次开发。以下是一些常用的Python库和工具：

pandas：用于数据处理和分析。

numpy：用于数值计算。

matplotlib：用于绘图和可视化。

subprocess：用于调用外部程序，如AERMOD。

5.1读取和处理输入文件

使用Python读取和处理AERMOD的输入文件，可以简化数据准备过程。

示例代码：

importpandasaspd

#读取源参数文件

defread_source_file(file_path):

读取源参数文件并返回DataFrame

:paramfile_path:源参数文件路径

:return:pandasDataFrame

columns=[SOURCEID,X,Y,Z,HR,MIN,SEC,Q,D,U,V,W]

returnpd.read_csv(file_path,sep=,,header=None,names=columns)