4.4.3 初始条件

第1步：

第2步：

通过计算2012年9月和2013年9月的两组观测井的数据的均方根误差，NRMSE= 0.05，认为该断时间的为近似均衡状态，利用这一段时间做稳定流模拟，将这一段时间的源汇项输入到模型中，不断调整参数，使得模拟值与观测值，接近，满足识别条件以后。以……时间段做非稳定流模拟，不断调整参数，得到给水度。

第3步：

（3）将水位观测孔作为constant head，将2012年9月-2013年9月其他的边界条件数据赋值到模型，运行模型，并根据整体流场分布特征不断地调整参数和初始水位。直至流场形态基本一致，即满足识别条件。此时，提取2016年6月各层的计算水头作为识别模型的初始条件。

第4步：

第5步：

第6步：

第7步：

 

由于收集到的地下水潜水等水位线数据为2012年9月。

包气带剖面土壤含水量的测试时间为2016年4月~6月。根据收集到资料的特征，2012年9月-2017年5月做饱和带模拟，并将识别后的2016年6月的水位数据作为饱和-非饱和水流耦合模拟的初始水位数据。对于饱和带水流模拟利用modflow-NWT模块进行模拟求解；对于饱和-非饱和水流耦合模拟利用modflow-surfact模块进行求解。

对研究区地下水等水位线在MAPGIS6.7软件中进行抽稀，得到研究区的初始水位数据，在visual modflow flex 6.1软件中利用Kringing方式进行插值，作为研究区的地下水模型的初始水位数据。由于研究区地下水位数据只集中在第一层，利用参数——初始水位迭代法（P-H0）法，在求取参数过程中随时调整初始水头分布，初始水头的分布是与模型识别同时完成的。具体操作如下：

（1）根据钻孔资料，在第三章对模型的结构进行概化，并根据各钻孔资料的分层结果，对各钻孔所在分层的综合水平渗透系数进行计算，根据计算的综合水平渗透系数的数值，结合观测孔的位置，对参数进行分区。

（2）利用Kringing插值法确定初始水头，并对每层的初始水头都赋相同的值。