水质模型软件：MIKE SHE二次开发_（1）.MIKE SHE软件基础理论.docx

PAGE1

PAGE1

MIKESHE软件基础理论

1.MIKESHE简介

MIKESHE是由丹麦水力研究所（DHI）开发的一款综合性的水文和水动力模型软件。该软件主要用于模拟地表水和地下水的流动、水质传输和转化过程，适用于水资源管理、洪水预测、水质评估等多种应用场景。MIKESHE通过将地表水、土壤水和地下水的流动过程有机地结合在一起，提供了一个全面的水文循环模型。

1.1MIKESHE的基本结构

MIKESHE的基本结构可以分为以下几个主要部分：

地表水模块：模拟地表径流、蒸发、降水等过程。

土壤水模块：模拟土壤中的水分传输、入渗、土壤蒸发等过程。

地下水模块：模拟地下水的流动、回补、开采等过程。

水质模块：模拟水质传输和转化过程，包括污染物的迁移、降解等。

界面模块：实现地表水、土壤水和地下水之间的相互作用。

1.2MIKESHE的应用领域

MIKESHE可以广泛应用于以下领域：

水资源管理：评估水资源的可用性和分配，优化水资源利用。

洪水预测：模拟流域内的洪水过程，提供预警和管理建议。

水质评估：评估水体的水质状况，模拟污染物的传输和转化。

生态水文：研究水文过程对生态系统的影响。

城市排水：模拟城市排水系统，优化设计和管理。

2.MIKESHE的数学模型

2.1水文循环的基本方程

MIKESHE的核心是基于水文循环的基本方程，主要包括以下几个部分：

地表水方程：描述地表水的流动和传输过程。

土壤水方程：描述土壤中的水分传输和转化过程。

地下水方程：描述地下水流的运动和传输过程。

水质传输方程：描述水质在水体中的传输和转化过程。

2.1.1地表水方程

地表水方程主要基于圣维南方程组（Saint-Venantequations），该方程组描述了地表水的连续性和动量守恒。具体方程如下：

?

?

其中：

A是过水断面面积。

Q是流量。

R是净流入量（包括降水、蒸发、入渗等）。

g是重力加速度。

z是水位。

Sf

S0

2.1.2土壤水方程

土壤水方程主要基于理查兹方程（Richardsequation），该方程描述了土壤中的水分传输过程。具体方程如下：

?

其中：

θ是土壤体积含水量。

Kθ

ψ是土壤基质势。

2.1.3地下水方程

地下水方程主要基于达西定律（Darcy’slaw）和连续性方程，描述了地下水的流动过程。具体方程如下：

?

其中：

h是地下水位。

K是渗透系数。

R是净流入量（包括回补、开采等）。

2.1.4水质传输方程

水质传输方程主要基于对流-弥散方程（Advection-Dispersionequation），描述了水质在水体中的传输和转化过程。具体方程如下：

?

其中：

C是水质浓度。

u是水流速度。

D是弥散系数。

R是水质源汇项（包括降解、沉积等）。

2.2MIKESHE的模型求解方法

MIKESHE采用数值方法求解上述方程，主要包括有限差分法（FiniteDifferenceMethod,FDM）和有限元法（FiniteElementMethod,FEM）。

2.2.1有限差分法

有限差分法通过将连续的偏微分方程离散化为差分方程，从而在计算网格上进行数值求解。例如，对于一维地表水方程，可以采用以下离散化形式：

A

Q

2.2.2有限元法

有限元法通过将连续的偏微分方程转化为积分方程，然后在计算网格上进行数值求解。例如，对于二维地下水方程，可以采用以下有限元形式：

Ω

其中：

?i

Ω是计算区域。

3.MIKESHE的数据输入与输出

3.1数据输入

MIKESHE需要多种类型的数据输入，主要包括：

地形数据：包括数字高程模型（DigitalElevationModel,DEM）和流域边界。

气象数据：包括降水量、蒸发量、气温等。

土壤数据：包括土壤类型、土壤渗透系数等。

地下水数据：包括地下水位、渗透系数等。

水质数据：包括水质浓度、水质源汇项等。

边界条件：包括地表水边界条件、地下水边界条件等。

3.1.1地形数据

地形数据是MIKESHE建模的基础，通常以数字高程模型（DEM）的形式提供。DEM可以用于生成流域边界、河流网络等。例如，可以使用GIS软件生成DEM文件，然后导入MIKESHE中。以下是生成DEM文件的Python代码示例：

importnumpyasnp

importrasterio

#生成一个示例的DEM数据

defgenerate_dem_data(shape,filename):

生成一个示例的DEM数据并保存为TIF文件

参数:

shape(tup