根據(jù)芬蘭Valkea-kotinen湖2011年1月14日至2011年4月3日冰下溶解氧濃度和溫度的現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)，利用分布參數(shù)系統(tǒng)參數(shù)辨識(shí)方法，對(duì)Golosov等人建立的淡水湖泊冰下溶解氧濃度一維偏微分方程模型中溶解氧湍流垂直交換系數(shù)進(jìn)行優(yōu)化辨識(shí)，得到了適應(yīng)該湖泊的溶解氧湍流垂直交換系數(shù)值，并利用該溶解氧湍流垂直交換系數(shù)值模擬了Valkea-kotinen湖冰下溶解氧濃度隨時(shí)間變化的情況。數(shù)值模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)吻合良好，說明辨識(shí)出的Valkea-kotinen湖冰下溶解氧湍流垂直交換系數(shù)是準(zhǔn)確有效的，該結(jié)果能對(duì)寒區(qū)冰下溶解氧濃度的數(shù)值模擬研究提供參考依據(jù)。

水體中溶解氧是養(yǎng)殖環(huán)境中最重要的因素之一，在養(yǎng)殖水體中既是水體理化性質(zhì)和生物學(xué)過程的綜合反映，也是養(yǎng)殖池塘生產(chǎn)性能的重要參數(shù)[1]。水體中溶解氧不僅是各種水生生物呼吸代謝的基礎(chǔ)，溶解氧水平的高低還直接反映了水體的質(zhì)量。它與養(yǎng)殖生物的生存繁衍和水體的自凈作用息息相關(guān)。

國(guó)內(nèi)對(duì)湖泊、海洋近岸等水域溶解氧作了許多的研究[2-3]。對(duì)于冬季溶解氧問題，蔣國(guó)昌等[4]分析了東海冬季環(huán)流中溶解氧和營(yíng)養(yǎng)鹽分布特征;王海英等[5]介紹了海河三岔口段冬季溶解氧升高的原因;石強(qiáng)等[6]對(duì)渤海冬季溶解氧與表觀耗氧量年際時(shí)空變化進(jìn)行了分析;王臻[7]對(duì)羅源灣冬季表層水體中溶解氧分布特征進(jìn)行討論。然而國(guó)內(nèi)對(duì)冰雪覆蓋下淡水湖泊溶解氧濃度的研究還比較匱乏，僅見李明等[8]對(duì)芬蘭淡水湖冰下溶解氧濃度變化規(guī)律進(jìn)行了分析。國(guó)外對(duì)于冰雪覆蓋下淡水湖泊溶解氧濃度問題已經(jīng)進(jìn)行了較深入的研究并取得一些成果。Terzhevik等[9]對(duì)淺層冰封湖水的溫度與溶氧結(jié)構(gòu)特點(diǎn)進(jìn)行了研究;Robarts等[10]分析了漢堡湖冬季在冰雪覆蓋下溶解氧并未大量減少的原因;Golosov等[11]對(duì)冰下溶解氧濃度問題給出了一種一維拋物型偏微分方程模型(即FlakeEco模型)，該模型中的溶解氧湍流垂直交換系數(shù)受不同湖泊環(huán)境因素影響而不盡相同，沒有明確給出該交換系數(shù)的取值。

本文主要是根據(jù)坐落于(北緯62°10′-62°20′，東經(jīng)33°10′-33°20′)的芬蘭Valkea-kotinen湖2011年1月13日至5月17日不同站位不同深度每隔17 min溶解氧濃度和溫度的現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)，利用分布參數(shù)系統(tǒng)參數(shù)辨識(shí)方法對(duì)該湖溶解氧湍流垂直交換系數(shù)K進(jìn)行優(yōu)化辨識(shí)，并利用辨識(shí)結(jié)果數(shù)值模擬該湖冰下溶解氧濃度隨時(shí)間變化的情況。數(shù)值模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)相吻合，表明本文辨識(shí)出的Valkea-kotinen湖冰下溶解氧的湍流垂直交換系數(shù)是準(zhǔn)確、有效的，該結(jié)果能對(duì)冰蓋下溶解氧濃度的數(shù)值模擬研究提供參考依據(jù)。

參數(shù)辨識(shí)和數(shù)值模擬結(jié)果

本文依據(jù)芬蘭Valkea-kotinen湖2011年1月13日至2011年5月17日3個(gè)站位采集的數(shù)據(jù)，采集的間隔為17min。測(cè)點(diǎn)位于湖泊中央及其兩側(cè)，深度級(jí)別的零點(diǎn)是冰水交界面。傳感器部件被連接到浮漂上，在冰下隨浮漂的起伏而震蕩，因此傳感器到冰底部的距離是固定的，但傳感器到湖底的距離隨著冰厚度變化而變化。湖泊凍結(jié)時(shí)間是2010年11月8日，湖面開始融化時(shí)間是2011年4月3日。利用NO1站位2011年1月14日至2011年4月3日的溶解氧濃度數(shù)據(jù)辨識(shí)計(jì)算出系數(shù)K，然后利用該K值對(duì)相同時(shí)間段的NO2站位溶解氧濃度隨時(shí)間變化情況進(jìn)行數(shù)值模擬。

辨識(shí)實(shí)例

以NO2站位該段時(shí)間冰下溶解氧濃度和水溫?cái)?shù)據(jù)為辨識(shí)依據(jù)。湖的深度為2.75m，測(cè)試溶解氧傳感器的分布位置分別為0.20、0.45、0.70、0.95m。在所測(cè)時(shí)間段內(nèi)4個(gè)測(cè)試點(diǎn)位置皆為有氧區(qū)域，符合用本文模型進(jìn)行描述條件。

在實(shí)際計(jì)算時(shí)，時(shí)間節(jié)點(diǎn)間距Δt取為60s，空間網(wǎng)格節(jié)點(diǎn)間距Δz取為0.0lm。將辨識(shí)時(shí)間段起始時(shí)刻各深度實(shí)測(cè)溶解氧濃度值每隔0.01m經(jīng)線性插值得到各節(jié)點(diǎn)的初始溶解氧濃度，以z1=0.20m及z2=0.95m處各時(shí)刻的實(shí)測(cè)溶氧濃度，經(jīng)插值得到的各時(shí)間節(jié)點(diǎn)的溶解氧濃度數(shù)據(jù)作為上、下邊界條件，進(jìn)行參數(shù)K的辨識(shí)計(jì)算。γ[T(z，t)]是與溫度有關(guān)的，根據(jù)該時(shí)間段冰下的實(shí)測(cè)溫度，按深度每隔0.01m、時(shí)間每隔60s進(jìn)行線性插值得到各節(jié)點(diǎn)處的溫度值，γmin=10-8s-1，γmax=5×10-6s-1[11]。

經(jīng)過大量的辨識(shí)計(jì)算，綜合分析得到了Valkea-kotinen湖冰封期溶解氧湍流垂直交換系數(shù)K=0.013。

數(shù)值模擬結(jié)果

為了檢驗(yàn)辨識(shí)得出的溶解氧湍流交換系數(shù)是否準(zhǔn)確有效，利用該系數(shù)對(duì)Valkea-kotinen湖NO2站位冰下溶解氧濃度變化情況進(jìn)行數(shù)值模擬。NO2站位湖的深度為5.5m，測(cè)試溶解氧及溫度傳感器的分布位置分別為0.20、0.45、0.70、0.95、1.95、2.95、3.95m。由于觀測(cè)后期2.95、3.95m位置實(shí)測(cè)數(shù)據(jù)基本為0，即缺氧狀態(tài)，不符合用該模型描述，因此以1.95m處觀測(cè)數(shù)據(jù)為下邊界。

根據(jù)NO2站位2011/01/14 0:11-2011/04/03 23:46實(shí)測(cè)數(shù)據(jù)，以z1=0.20m、z2=1.95m各時(shí)刻的實(shí)測(cè)數(shù)據(jù)作為上、下邊界條件。與辨識(shí)K值相同方法求出各節(jié)點(diǎn)處初始和邊界溶解氧濃度數(shù)據(jù)，溫度也采用相同處理方法。然后對(duì)Valkea-kotinen湖NO2站位溶解氧濃度變化情況進(jìn)行了數(shù)值模擬，圖1繪出了NO2站位數(shù)值模擬和實(shí)測(cè)數(shù)據(jù)對(duì)比曲線，其中虛線表示數(shù)值模擬計(jì)算數(shù)據(jù)。

定義計(jì)算溶解氧濃度與實(shí)測(cè)濃度的誤差為

當(dāng)選定K=0.013時(shí)，計(jì)算與實(shí)測(cè)誤差為3.32%?？梢钥闯?，計(jì)算結(jié)果與現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)吻合較好，說明本文辨識(shí)出的溶解氧湍流垂直交換系數(shù)K=0.013用于計(jì)算Valkea-kotinen湖溶解氧濃度是可行的。

結(jié)論

本文根據(jù)芬蘭Valkea-kotinen湖溶解氧濃度和溫度的現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)，利用分布參數(shù)系統(tǒng)參數(shù)辨識(shí)方法對(duì)該湖NO1站位溶解氧湍流垂直交換系數(shù)K進(jìn)行優(yōu)化辨識(shí)，然后用辨識(shí)出的系數(shù)K對(duì)該湖NO2站位溶解氧濃度變化情況進(jìn)行數(shù)值模擬。數(shù)值模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)吻合良好，表明本文辨識(shí)出的Valkea-kotinen湖冰下溶解氧湍流垂直交換系數(shù)是準(zhǔn)確、有效的，該結(jié)果能對(duì)寒區(qū)冰下溶解氧濃度的數(shù)值模擬研究提供參考依據(jù)。