协同克里格方法在东海表面温度场数据插值中的应用 Application of Co-Kriging Method to Interpolation of Surface Temperature Data in East China Sea.pdfVIP

第3O卷第3期 海 岸 工 程 20l1年9月 文章编号 ：1002—3682(2011)03—0049—07 协同克里格方法在东海表面温度 场数据插值 中的应用 宋庆磊，杜德文 ，丁 明，闫仕娟 (国家海洋局 第一海洋研究所 ，山东 青岛 266061) 摘 要：首次将协 同克里格方法引入到海洋水文插值计算 中处理东海实验 区 2008—01 2。O8一 O3三个月抽稀后的月平均海表面温场数据 。将温场数据分为趋势分量和剩余分量 ，最终采 用两分量的估计值之和实现海表面温场数据的格网插值估计 。结果表 明：协 同克里格方法比 传 统方法 (趋势分析 、移动平均 、普通克里格插值)效果好 。提 出变差函数局部拟合的方法 ，解 决了在最大搜寻范 围内变差 函数无法收敛 ，各方向很难取得统一基 台值 的问题 。 关键词 ：普通克里格 ；协 同克里格 ；SST；变差函数局部拟合 中图分类号：P731．11 文献标识码 ：A 水温数据是水文研究的重要内容之一 ，具有空间位置和时间两大基础特征 。由于调 查条件限制造成的数据缺失，数据本身存在 的系统或 随机误差 ，采样密度过稀等原因，使 得数据在空间分布上往往具有不均匀性 ；而在调查时间上 ，也很难保证所有调查站位在统 一 时间内获取数据。将这些来源不 同、精度不一 的水温数据转化为某一时 间段平均 的格 网数据是水温原始调查数据转化为数据产品的基础性工作 。因此，有必要对原始测量数 据进行插补或插值处理 以及对数据进行时间段统计处理分析 ，以期在其客观性 、科学性和 完备性上作出提升 。 区域水文数据可看成 由物理机理控制的趋势场分量和受随机因素控制的剩余场分量 两部分组成 。物理趋势场分量可视为确定因素 ，可由理论模拟 ，随机场部分这里引入地质 统计方法进行处理 。将两个分量分别进行插值处理 ，然后求和得到最终插值结果 。本研 究主要讨论海表面温场数据的处理。 地质统计学是数学地质领域 中一 门发展迅速且有着广泛应用前景的新兴学科 ，它是 以区域化变量理论为基础 ，以变异函数为主要工具 ，研究在空间分布上既有随机性又有结 构性 ，或空间相关和依赖性的 自然现象 的科学_】]。其 中的克里格法是一种最优的无偏内 插估计方法 j，与传统插值方法相 比[4]，它不仅能对不规则格点进行插值 ，而且还能根据 多种要素场信息对某一要素场进行协同插值 ]。 目前地质统计学在海洋领域主要应用于 海洋地质 、海洋矿物勘探 ，另外在生物学、生态学、农业 、水文地质、气象学、环境 、林业 、医 学等许多方面都有成功应用的实例 ]。本研究首次将协同克里格方法引入到海洋水文插 收稿 日期 ：2011-06—26 作者简介：宋庆磊 (1981一)，男 ，硕士 ，主要从事海洋数据处理方面研究 ．E—mail：songql@fio．org．cn (陈 靖 编辑) 5O 海 岸 工 程 第 3期 值计算中，并在这一方向上初步尝试 ，分析比较 了趋势面插值 、移动平均法、普通克里格方 法 。和协同克里格方法 。在东海表面温场数据插值中的效果。 1 技术流程 技术流程分为以下几部分 ：实验数据准备 、平稳检验 、相关性检验 、实验变差函数计算 与拟合 、协N／普通克里格插值 、求得海表面温度均估计值、估计效果分析和分析讨论得出 结论 (图 1)。 相邻时间域的空间实验数据准备 是 ／-t一阶平稳、二阶平稳检验 }_— 否 剩余场 卜_— 去除物理机理 区域化变量 时间域相关性检验