（二）初始条件

从1:20万地形图和2012年潜水和承压水水位等值线图上获取地表高层及地下水位数据，利用Surfer软件中的Kriging插值方法，确定模拟区各单元的地表高程值和地下水位数据。 根据所建立的数学模型，将研究区第四系地下水系统分为6个模拟层。数值模拟需要对6个模拟层进行初始水头赋值，而研究区地下水位数据主要集中在第1、4层，不能满足模型对于初始水头分布的要求，本文采用以下方法来确定各个模拟层的初始水头：

（1）根据钻孔岩性资料和抽水试验结果确定各分区的水文地质参数初始值，然后利用Kriging插值法确定各单元的初始水头，假定各模拟层初始水头分布相同，对各层同一平面位置的单元赋相同初始水头值；

（2）将水位观测孔作为一类边界点，依据降水、蒸发等气象数据，将2012年9月至2016年8月间的边界条件赋入模型，运行模型，并依据整体流场分布规律对模型参数进行调整，直至流场形态基本一致，并提取2016年8月各层的计算水头作为识别模型的初始条件；

（3）在（2）模型的基础上，以计算水位作为初始条件，将2016年8月至2017年7月的观测孔位作为已知水头边界，将期间的河流、渠系、水库、田间灌溉等入渗补给量、侧向径流补给、蒸散发、地下水开采量、侧向排泄、、等边界条件输入模型，计算1年后各单元的水头值；重复运算多次，得到一组水头分布；

（4）撤销观测孔所在的一类边界点，将（2）中获得的水头分布作为初始条件，赋值模型所需要的源汇项后运行模型，采取常规优选法初步确定水文地质参数；

（5）采用经过识别的水文地质参数代替初始估计值，多次循环迭代，当拟合达到一定程度后全面检验模型是否符合基本规律，否则重新调整直至满足要求。