一种插值相似方法作温度预测的应用研究

将插值方法运用到短期气象预测中，其基本原理是按照相似原理将产生相似结果来作预报。对云南昆明等四站温度进行了九年定性和定量预报检验，结果表明：这种方法具有良好的预报能力，可以作要素预报应用。     

基于二维三次卷积插值算法的辛几何射线追踪

射线追踪是地震波走时层析成像的基础,射线空间位置的准确性及射线走时的精度决定了层析成像的可靠性.本文根据哈密尔顿系统可以有效提高程函方程解稳定性的特性,采用辛几何算法（SAM-Symplectic Algorithm Method）及二维三次卷积插值技术进行地震波射线追踪.由于采用了SAM算法,保证了地震波波前精度,提高了射线空间位置的准确性.数值模拟结果表明SAM既能保证哈密尔顿系统的稳定性又具有运算速度快的特点,提高了射线追踪的计算精度.     

地层格架正则化约束下的二维立体层析反演

通过把地层格架信息作用于立体层析Fréchet导数矩阵,使得更新后的速度模型呈现出符合地质规律的块状特征.地层格架信息基于立体层析反演中得到的反射点位置进行非规则B样条插值拟合得到,因此在反演中它将会随着反射点位置的更新自然得到更新.与前人提出的保边缘层析算法或多层立体层析算法相比,本文提出的地层格架正则化无需引入混合正则化项或定义某种复杂的混合速度格式,更为直接也更容易实现.理论和实际数据算例证实了该正则化技巧的稳健性和可靠性,能够得到与实际地质构造特征更为一致的地质一致性反演结果.

     

赋予地质信息的采空区三维联合反演

浅部采空区探测一般采用高密度电法,但该方法的探测结果一般都是剖面形式.在圈定的大范围采空区内,高密度法也只能确定测点周边电阻或某一条剖面的电阻分布情况,并不能很好的利用其他周边测线的结果和采集到的地质标本以及其他地质信息.地质信息也只是作为反演结果的判定条件,不参与计算过程,同时以往采集的大量的二维高密度数据也需要得到进一步开发计算.因此,对地质信息进行赋值,通过标准偏差确定赋值的合理性,最后对联合地质资料进行反演迭代计算,得到三维铁矿采空区结果,并与二维反演结果、实地际踏勘情况、未添加赋值信息的三维反演结果相对比,赋予地质信息的三维反演结果对采空区、山体裂隙的识别更精细直观.     

基于多参量的模拟退火全局优化傅里叶有限差分算子

傅里叶有限差分法(FFD)能够处理复杂地质构造中的波传播问题,但对陡倾角成像仍有明显的误差.优化参数的方法能够在保持计算效率的前提下进一步提高陡倾角的成像精度.本文在有理近似的基础上,将FFD算子展开式中的常系数由两个拓展为四个,然后采用模拟退火算法对这四个参数进行全局优化.本方法除了考虑速度对比度以外,还考虑了频率和延拓步长等参量的影响.理论误差分析和脉冲响应测试均表明该方法能极大地提高FFD算子的精确传播角度.二维SEG/EAGE盐丘模型实验表明本文方法对陡倾角以及盐下构造的成像精度明显高于未优化的FFD法.将本文的方法与交替方向加插值的方法结合应用于三维脉冲响应测试更进一步证实了本文方法的有效性.     

基于回声探测监测沉水植物盖度的插值方法研究:以杭州西湖沉水植物盖度监测为例

沉水植物作为水生态系统的重要组成成分,在水生态系统物质循环和能量流动中发挥着重要作用,其覆盖度和生物量是评价湖泊等浅水水体系统稳定性的关键参数随着高效和无损伤监测的回声探测仪在沉水植物盖度监测中的应用,其精确度算法也受到了越来越多的关注本研究以成功恢复沉水植物的浅水湖泊杭州西湖为研究对象,利用BioSonics便携型回声探测仪——MX采集沉水植物回声样本同时结合人工样方设置,采集与回声探测对应位点的沉水植物样本,验证回声探测结果的精确性通过建立回归模型分析回声探测得到的沉水植物体积百分比(PVI)与人工样方获得的对应平均鲜重关系,结果表明二者具有较好的相关性分别采用普通克里金法、反距离权重法、径向基函数法3种插值方法对同一季节的不同湖泊和同一湖泊的不同季节未采集区域沉水植物的盖度数据进行插值分析,并对插值结果进行交叉验证,以确定方法的精确度交叉验证结果表明,插值精确度反距离权重法径向基函数法普通克里金法研究结果为回声探测与插值分析方法结合在大尺度浅水水体中沉水植物监测应用提供了技术支撑.     

气候模式模拟中海量数据的并行抽取与可视化分析诊断

在地球动力学和气候模拟等领域, 数值模拟产生的数据规模达到Tb至Pb量级。 实现这些海量数据的实时可视化和实时诊断分析面临很大的困难, 时空多尺度数据抽取可以解决这一瓶颈。 高精度数据的可视化结果展示需要高分辨率的显示设备, 并行的大屏幕显示技术是解决这一问题的有效手段。 地球科学数值模拟、 并行数据抽取和高分辨率显示都需要搭建高性能计算机集群。 本文在基于Lagrange插值的多维度、 多尺度、 多分辨率并行数据抽取算法的同时, 利用并行计算节点及LCD显示器, 基于Rocks cluster系统搭建起一个176核, 4×10×1024×1280分辨率的高性能计算模拟、 数据抽取和并行显示输出的集成平台, 并将该平台成功应用到气候模式模拟产生的海量数据的并行抽取和并行显示。     

汶川8.0级地震前后全国倾斜场变化初步分析

由于岩石介质的不均匀型,在易破区域会出现应力积累,伴随着倾斜会出现较大的变化。由岩石破裂实验可以看出破裂区在破裂前出现了体积变大（鼓包）或变小（凹坑）现象,如果在岩石表面放置倾斜仪器,我们不难得出各仪器的倾斜方向指向或背向破裂点。本文采用全国前兆台网地倾斜（水管）数据,使用5次曲面函数区域插值的方法,拟合出地倾斜场（图1）。     

汶川8．O级地震前后全国倾斜场变化初步分析

由于岩石介质的不均匀型,在易破区域会出现应力积累,伴随着倾斜会出现较大的变化。由岩石破裂实验可以看出破裂区在破裂前出现了体积变大（鼓包）或变小（凹坑）现象,如果在岩石表面放置倾斜仪器,我们不难得出各仪器的倾斜方向指向或背向破裂点。本文采用全国前兆台网地倾斜（水管）数据,使用5次曲面函数区域插值的方法,拟合出地倾斜场（图1）。     

重力时变的数学逼近

研究了根据时域内离散的重复重力测量数据获取连续重力变化的方法。分析比较了多种不同的插值、拟事和滤波方法的优缺点，认为Akima样条插值法，正交多项式拟合法和卡尔曼滤波法对重力数据处理效果较好，以华北63号相对重力点和香山绝对点为例给出了处理结果，并对三种方法的意义和适用范围进行了讨论。