根据芬兰Valkea-kotinen湖2011年1月14日至2011年4月3日冰下溶解氧浓度和温度的现场观测数据，利用分布参数系统参数辨识方法，对Golosov等人建立的淡水湖泊冰下溶解氧浓度一维偏微分方程模型中溶解氧湍流垂直交换系数进行优化辨识，得到了适应该湖泊的溶解氧湍流垂直交换系数值，并利用该溶解氧湍流垂直交换系数值模拟了Valkea-kotinen湖冰下溶解氧浓度随时间变化的情况。数值模拟结果与实测数据吻合良好，说明辨识出的Valkea-kotinen湖冰下溶解氧湍流垂直交换系数是准确有效的，该结果能对寒区冰下溶解氧浓度的数值模拟研究提供参考依据。


水体中溶解氧是养殖环境中最重要的因素之一，在养殖水体中既是水体理化性质和生物学过程的综合反映，也是养殖池塘生产性能的重要参数[1]。水体中溶解氧不仅是各种水生生物呼吸代谢的基础，溶解氧水平的高低还直接反映了水体的质量。它与养殖生物的生存繁衍和水体的自净作用息息相关。

国内对湖泊、海洋近岸等水域溶解氧作了许多的研究[2-3]。对于冬季溶解氧问题，蒋国昌等[4]分析了东海冬季环流中溶解氧和营养盐分布特征;王海英等[5]介绍了海河三岔口段冬季溶解氧升高的原因;石强等[6]对渤海冬季溶解氧与表观耗氧量年际时空变化进行了分析;王臻[7]对罗源湾冬季表层水体中溶解氧分布特征进行讨论。然而国内对冰雪覆盖下淡水湖泊溶解氧浓度的研究还比较匮乏，仅见李明等[8]对芬兰淡水湖冰下溶解氧浓度变化规律进行了分析。国外对于冰雪覆盖下淡水湖泊溶解氧浓度问题已经进行了较深入的研究并取得一些成果。Terzhevik等[9]对浅层冰封湖水的温度与溶氧结构特点进行了研究;Robarts等[10]分析了汉堡湖冬季在冰雪覆盖下溶解氧并未大量减少的原因;Golosov等[11]对冰下溶解氧浓度问题给出了一种一维抛物型偏微分方程模型(即FlakeEco模型)，该模型中的溶解氧湍流垂直交换系数受不同湖泊环境因素影响而不尽相同，没有明确给出该交换系数的取值。


本文主要是根据坐落于(北纬62°10′-62°20′，东经33°10′-33°20′)的芬兰Valkea-kotinen湖2011年1月13日至5月17日不同站位不同深度每隔17 min溶解氧浓度和温度的现场观测数据，利用分布参数系统参数辨识方法对该湖溶解氧湍流垂直交换系数K进行优化辨识，并利用辨识结果数值模拟该湖冰下溶解氧浓度随时间变化的情况。数值模拟结果与实测数据相吻合，表明本文辨识出的Valkea-kotinen湖冰下溶解氧的湍流垂直交换系数是准确、有效的，该结果能对冰盖下溶解氧浓度的数值模拟研究提供参考依据。


参数辨识和数值模拟结果


本文依据芬兰Valkea-kotinen湖2011年1月13日至2011年5月17日3个站位采集的数据，采集的间隔为17min。测点位于湖泊中央及其两侧，深度级别的零点是冰水交界面。传感器部件被连接到浮漂上，在冰下随浮漂的起伏而震荡，因此传感器到冰底部的距离是固定的，但传感器到湖底的距离随着冰厚度变化而变化。湖泊冻结时间是2010年11月8日，湖面开始融化时间是2011年4月3日。利用NO1站位2011年1月14日至2011年4月3日的溶解氧浓度数据辨识计算出系数K，然后利用该K值对相同时间段的NO2站位溶解氧浓度随时间变化情况进行数值模拟。


辨识实例


以NO2站位该段时间冰下溶解氧浓度和水温数据为辨识依据。湖的深度为2.75m，测试溶解氧传感器的分布位置分别为0.20、0.45、0.70、0.95m。在所测时间段内4个测试点位置皆为有氧区域，符合用本文模型进行描述条件。


在实际计算时，时间节点间距Δt取为60s，空间网格节点间距Δz取为0.0lm。将辨识时间段起始时刻各深度实测溶解氧浓度值每隔0.01m经线性插值得到各节点的初始溶解氧浓度，以z1=0.20m及z2=0.95m处各时刻的实测溶氧浓度，经插值得到的各时间节点的溶解氧浓度数据作为上、下边界条件，进行参数K的辨识计算。γ[T(z，t)]是与温度有关的，根据该时间段冰下的实测温度，按深度每隔0.01m、时间每隔60s进行线性插值得到各节点处的温度值，γmin=10-8s-1，γmax=5×10-6s-1[11]。


经过大量的辨识计算，综合分析得到了Valkea-kotinen湖冰封期溶解氧湍流垂直交换系数K=0.013。


数值模拟结果


为了检验辨识得出的溶解氧湍流交换系数是否准确有效，利用该系数对Valkea-kotinen湖NO2站位冰下溶解氧浓度变化情况进行数值模拟。NO2站位湖的深度为5.5m，测试溶解氧及温度传感器的分布位置分别为0.20、0.45、0.70、0.95、1.95、2.95、3.95m。由于观测后期2.95、3.95m位置实测数据基本为0，即缺氧状态，不符合用该模型描述，因此以1.95m处观测数据为下边界。


根据NO2站位2011/01/14 0:11-2011/04/03 23:46实测数据，以z1=0.20m、z2=1.95m各时刻的实测数据作为上、下边界条件。与辨识K值相同方法求出各节点处初始和边界溶解氧浓度数据，温度也采用相同处理方法。然后对Valkea-kotinen湖NO2站位溶解氧浓度变化情况进行了数值模拟，图1绘出了NO2站位数值模拟和实测数据对比曲线，其中虚线表示数值模拟计算数据。


定义计算溶解氧浓度与实测浓度的误差为


当选定K=0.013时，计算与实测误差为3.32%。可以看出，计算结果与现场观测数据吻合较好，说明本文辨识出的溶解氧湍流垂直交换系数K=0.013用于计算Valkea-kotinen湖溶解氧浓度是可行的。


结论


本文根据芬兰Valkea-kotinen湖溶解氧浓度和温度的现场观测数据，利用分布参数系统参数辨识方法对该湖NO1站位溶解氧湍流垂直交换系数K进行优化辨识，然后用辨识出的系数K对该湖NO2站位溶解氧浓度变化情况进行数值模拟。数值模拟结果与实测数据吻合良好，表明本文辨识出的Valkea-kotinen湖冰下溶解氧湍流垂直交换系数是准确、有效的，该结果能对寒区冰下溶解氧浓度的数值模拟研究提供参考依据。