水质模型软件：SWMM二次开发_（8）.SWMM扩展模块开发.docx

PAGE1

PAGE1

SWMM扩展模块开发

在上一节中，我们介绍了SWMM的基本功能和使用方法，包括其在城市排水系统和水质模拟中的应用。本节将重点探讨如何通过二次开发来扩展SWMM的功能，使其能够更好地满足特定项目的需求。SWMM的扩展模块开发主要涉及以下几个方面：自定义水质参数、开发新的水文过程、集成外部数据源以及优化模型性能。我们将会通过具体的例子来详细说明这些内容。

1.自定义水质参数

SWMM默认提供了一些水质参数，如COD、BOD、SS等，但在实际应用中，往往需要模拟更多特定的水质参数。自定义水质参数允许用户根据项目的具体需求添加新的水质变量，并在模型中进行计算和处理。

1.1原理

SWMM的水质模块是基于质量守恒原理的，即流入系统的水质质量等于流出系统的水质质量和系统内部的质量变化之和。自定义水质参数的原理与默认参数相同，但需要在源代码中添加新的变量和计算逻辑。

1.2内容

定义新的水质参数：首先需要在SWMM的输入文件中定义新的水质参数。这通常通过修改pollutants部分来实现。

修改源代码：在SWMM的源代码中，需要添加新的水质变量的定义，并在相应的计算函数中实现新的逻辑。

编译和测试：修改源代码后，需要重新编译SWMM，并进行测试以确保新的水质参数能够正确地参与模拟计算。

1.3例子

假设我们需要添加一个新的水质参数：硝酸盐（NO3）。

1.3.1定义新的水质参数

在SWMM的输入文件中，添加新的水质参数NO3：

[OPTIONS]

...

StartDate01/01/2023

StartTime00:00:00

EndDate01/02/2023

EndTime23:59:59

...

[POLLUTANTS]

NameUnitsCoeffExpKcKfLimit

NO3mg/L0.50.30.010.0210.0

1.3.2修改源代码

在SWMM的源代码中，找到处理水质参数的部分，例如Swmm5.c文件中的InitPollutants函数，添加新的NO3参数：

//Swmm5.c

voidInitPollutants(void)

{

intj,k;

for(j=0;jNobjects[POLLUT];j++)

{

pollut[j].index=j;

pollut[j].units=GetUnits(j);

pollut[j].coeff=getDouble(Coeff,j);

pollut[j].exp=getDouble(Exp,j);

pollut[j].kc=getDouble(Kc,j);

pollut[j].kf=getDouble(Kf,j);

pollut[j].limit=getDouble(Limit,j);

//新增NO3参数

if(strcmp(pollut[j].name,NO3)==0)

{

pollut[j].type=NO3;

}

}

}

在处理水质计算的函数中，例如Routing.c文件中的RoutePollutants函数，添加NO3的计算逻辑：

//Routing.c

voidRoutePollutants(doubledt)

{

inti,j,k;

for(i=0;iNobjects[SUBCATCH];i++)

{

for(j=0;jNobjects[POLLUT];j++)

{

doubleconc=subcatch[i].pollut[j].conc;

doubleload=subcatch[i].pollut[j].load;

//计算默认参数

if(pollut[j].type!=NO3)

{

//默认