水质模型软件：EFDC二次开发_（5）.EFDC中的水质参数设置.docx

PAGE1

PAGE1

EFDC中的水质参数设置

在上一节中，我们介绍了EFDC的基本概念和模型结构。本节将详细探讨EFDC中的水质参数设置，这是水质模拟中非常重要的一部分，直接关系到模拟结果的准确性和可靠性。水质参数的设置涉及到模型中各种物质的物理、化学和生物过程，包括物质的传输、转化、沉降和积累等。通过合理设置水质参数，可以有效地反映水质变化的复杂机制，从而为环境管理和决策提供科学依据。

1.水质参数的基本概念

水质参数是指描述水体中各种物质（如溶解氧、营养盐、有机物、重金属等）的物理、化学和生物特性的参数。在EFDC模型中，这些参数主要包括：

溶解氧（DissolvedOxygen,DO）：表示水体中溶解的氧的浓度，单位通常为mg/L。

营养盐（Nutrients）：包括氮（N）、磷（P）等，主要影响水体中的藻类生长。

有机物（OrganicMatter）：包括BOD（生化需氧量）、COD（化学需氧量）等，主要影响水体的耗氧过程。

重金属（HeavyMetals）：如铅、汞、镉等，主要影响水体的生态健康。

悬浮颗粒物（SuspendedSolids,SS）：主要影响水体的浑浊度和沉积过程。

2.溶解氧参数设置

溶解氧是水体中最重要的水质参数之一，直接影响水生生物的生存和水体的自净能力。在EFDC中，溶解氧的模拟主要涉及以下几个过程：

大气复氧：水体与大气之间的氧气交换过程。

光合作用：藻类通过光合作用产生氧气。

呼吸作用：水生生物和藻类的呼吸作用消耗氧气。

生化耗氧：有机物的分解过程消耗氧气。

沉降和沉积：溶解氧在水体中的垂直分布受到沉降和沉积过程的影响。

2.1大气复氧

大气复氧是指水体与大气之间的氧气交换过程。在EFDC中，大气复氧的速率通常用复氧系数（K2）来表示。复氧系数的计算公式如下：

K

其中：

K2

T：水体温度，单位为°C。

W：风速，单位为m/s。

U：水体表面流速，单位为m/s。

Wref和Uref

2.1.1代码示例

以下是一个Python代码示例，用于计算大气复氧系数：

#计算大气复氧系数

defcalculate_reaeration_coefficient(K2_ref,T,W,U,W_ref=1,U_ref=1):

计算大气复氧系数

参数:

K2_ref(float):参考温度下的复氧系数（1/d）

T(float):水体温度（°C）

W(float):风速（m/s）

U(float):水体表面流速（m/s）

W_ref(float):参考风速（m/s），默认为1m/s

U_ref(float):参考流速（m/s），默认为1m/s

返回:

float:大气复氧系数（1/d）

K2=(K2_ref*(T+273.15)**1.45*W**0.5*U**0.6)/((273.15+20)**1.45*W_ref**0.5*U_ref**0.6)

returnK2

#示例数据

K2_ref=0.5#参考温度下的复氧系数

T=25#水体温度

W=3#风速

U=0.5#水体表面流速

#计算复氧系数

K2=calculate_reaeration_coefficient(K2_ref,T,W,U)

print(f大气复氧系数K2:{K2:.4f}1/d)

2.2光合作用和呼吸作用

光合作用和呼吸作用是影响溶解氧浓度的重要生物过程。EFDC中通常用以下公式来描述这些过程：

Δ

其中：

P：光合作用产生的氧气量，单位为mg/L/d。

R：呼吸作用消耗的氧气量，单位为mg/L/d。

2.2.1代码示例

以下是一个Python代码示例，用于计算光合作用和呼吸作用对溶解氧的影响：

#计算光合作用和呼吸作用对溶解氧的影响

defcalculate_do_change(P,R):

计算光合作用和呼吸作用对溶解氧的影响

参数:

P(float):光合作用产生的氧气量（mg/L/d）

R(float):呼吸作用消耗的氧气量（mg/L/d）

返回:

float:溶解氧的变化量（mg/L/d）

delta_do=P-R

returndelta_do

#示例数据

P=5.0#光合