确定等效水力参数的方法

第一步

把王建军的模型复制为两份，见下表。

编号	实测值 (mm)	模型	模拟值 (mm)
B1	8.94	S1	7.61
B2	48.61	S2	45.45
B3	65.12	S5	63.73
B4	6.06	S6	5.75
B5	12.23	S3	12.42
B9	37.32	S11	36.21
B10	63.77	S8	58.44
B12	15.26	S13	16.68

B1模拟是365天

sum(vTop)[L]就是利用模型中这个值与实际蒸散发量进行对比的。

通过王建军的参数敏感性分析可知，

包气带岩性相对于大气降水、灌溉、水位埋深等因素，具有更高的空间变异性及随机性，土壤质地决定了包气带水力特征，是导致数值模拟结果不确定性的首要因素 (霍思远和靳孟贵, 2017)。

本次研究采用单因素扰动法（One factor at a time, OAT），选择±5%为变化幅度，对van Genuchten模型参数（θr、θs、Ks、α和n）的敏感性进行分析 (任理和薛静, 2017)。

通过对各参数在一定范围内变化（其他参数不变）的情况下，

以模拟计算的蒸发量作为指标分析出最为敏感的参数。

以S8试验点为例，土壤水力参数变化对蒸发量的影响曲线如图5.7所示。

从图5.7可以看出，土壤水力参数的敏感性由大到小排序为n、θs、α、Ks、θr。

因此在对模型进行参数识别时，应先调整参数n，再逐步优化模型。 

所以，将模型概化为一层以后，先调整n再逐步优化模型，

王建军（2020）将模型概化为三层，这些都是调整后的参数，我应该怎么去调整这个参数。

B.3. 模块 PROFILE

B.3.1. 土壤剖面离散化

模块 PROFILE（图 B.2）主要用于将一维土壤剖面离散化为离散节点，其中节点是通过将土壤剖面划分为小元素而生成的。

如果之前不存在节点，则程序会自动生成默认的等距点分布。

用户可以编辑节点的位置以优化不同元素的厚度。

有两种指定节点适当分布的方法，即通过

（1）编辑点的数量和

（2）指定固定点和节点密度。节点密度确定元素的相对长度，只能在固定点指定。

固定点可以在土壤剖面的任何位置插入或删除。

用户可以编辑固定点处的节点密度，以在该点周围局部细化节点分布。

仔细放置节点非常重要，因为节点分布在很大程度上决定了计算的最终质量和速度。

图 B.2. PROFILE 模块的主窗口；用于土壤剖面离散化。

B.3.2. 土壤剖面内土壤性质的规范

PROFILE 模块（图 B.3）帮助用户定义流域参数的空间分布（例如，土壤材料的空间分布、水力缩放因子、根吸水参数）和/或观测节点。

该模块中的所有参数都在图形环境中使用鼠标指定。

规定流域参数（初始条件、材料分布）相对简单。

用户首先必须选择要为其分配特定变量值的输运域的部分。

可以选择整个输运域、部分输运域或仅个别节点。

选择输运域的特定部分的步骤如下：用户首先单击“编辑条件”按钮，然后将鼠标移动到所选位置。

通过单击左按钮框定选择操作的起始和结束。

所选区域是由左按钮单击时鼠标位置定义的垂直区域。

选择完成后，将弹出“条件规范”窗口，用户必须指定特定变量的值。

然后将该值分配给所选区域。

在规定初始条件时，用户可以将常量值分配给所选域，也可以指定所选区域的顶部和底部的不同值，此时程序将在所选区域内线性插值变量。

变量始终分配给节点点，而不是元素。

英文术语	中文翻译	描述
Nodal Point	节点	节点是将土壤剖面离散化的节点， 并用绿色十字标记。这些节点从顶部（节点编号1） 到底部（节点NumNP）进行排序。
Elements	元素	元素是离散化土壤剖面的层，连接生成的节点。
Fixed Points	固定点	固定点是土壤剖面中用紫色星号标记的点。 这些点可用于调整节点点的局部离散化密度。 默认情况下，固定点位于土壤剖面的顶部和底部， 但也可以在土壤剖面的任何其他点插入或删除。
Nodal Density	节点密度	节点密度是一个实数，范围在<0.01, 100.>， 用于指定节点点的局部密度。密度只能在固定点指定。 程序区分顶部和底部密度。 在固定点处，顶部（底部）密度指定该点上方（或下方）元素的相对厚度。 如果顶部和底部密度相等，则节点密度在整个剖面上是连续的， 即两个元素具有相同的厚度。如果顶部和底部密度值不同， 则元素的厚度也将不同。 例如：如果DT = 3.，DB = 2.，则LT / LB = 1.5， 其中DT和DB分别是固定点处的顶部和底部密度， LT和LB分别是该固定点上方和下方元素的厚度。

图 B.3. PROFILE 模块的主窗口；用于规定土壤性质。

10、土壤剖面-图像编辑器

剖面-图像编辑器可对模拟土层的离散划分、土壤质地划分、土壤亚区划分、土壤初始的含水量（负压）等条件进行编辑，并用图像的形式显现出来。

在此步骤中选择菜单栏中的condition，其下有6 个选项。

Material Distribution	剖面离散
Root Distribution	根系分布
Scaling Factors	尺度因素
Initial Conditions	初始条件（初始的负压含水量分布等）
Subregion Distribution	子区划分
Observation Nodes	设置观察点

其中在2.模拟土壤的几何属性中，已经输入了土质的种类、土壤亚区的划分数目等，这里还要指定土壤整个剖面中，亚区的起止范围及对应的质地。

点击Edit condition，然后在所示土柱中拖动鼠标进行编辑。

剖面离散这一项根据自己的需要，将土壤划分有限元计算单元。

里面对Edit Node 进行编辑，可以编辑节点数目，节点划分密度等。

设置观察点时，直接用鼠标在土柱需要观察的点进行左击。

设置初始条件和根系分布时，需要设定范围边界和边界点的值，再使值线性分布。

此处设置往往比较麻烦，一般下一步还可以设置。