Geostatistics applied to best-fit interpolation of orientation data

Automatic methods used in geosciences to interpolate between orientation data have often limited applicability and strength, in particular where large ranges of orientations occur. In this paper, we show that geostatistical methods yield rather strong and powerful results when applied to directional data. The procedure involves the calculation of variograms, followed by a kriging interpolation of the data. In order to free from the circular property of directional data, the treatment of initial angular data sets is performed using scalar values provided by the direction cosines of double-angle values. The strength and application of the method are demonstrated by the analysis of theoretical and natural data sets. Natural examples are focused on the calculation and the analysis of cleavage trajectory maps.     

一种大区域地貌矢量数据的DEM快速生成与粗差检测算法

本文论述了一种大区域范围从矢量化地貌数据中快速形成ＤＥＭ以及可视化粗差检测与修正的方法。实践表明，该算法具有实用、直观、快速等特点，可获得稳定而可靠的计算结果。     

区域高程异常内插方法研究

介绍了部分高程异常内插的方法,结合一个具体的实例,在2km×2km的小范围内采用多项式拟合内插的方法进行高程异常的逼近,运用GPS水准获得的数据进行检查,通过分析,应该采用5参数内插的方法,取得了2.8cm的精度;并据此进行高程异常内插分析,运用MATLAB编写程序对几种插值方法进行讨论。     

基于修正πHW变分的板非协调动态有限元法

本文本文讨论了基于修正πＨＷ变分的杂交应变板弯曲非协调动态有限元的分析列式,建立了四节点单元,进行了特征分析。算例表明,该单元是简单、有效的。     

数字高程模型的随机分形算法研究

本文研究了基于随机分形的数字高程模型算法。实践证明，用此方法进行形状逼近内插，所需数据量小，可节省大量内存及计算时间，较容易获得所需要的分辨率。     

复杂山区地形条件下ERA5再分析地表气温降尺度方法

高时空分辨率的气温栅格数据是多种地学模型和气候模型的重要输入。山区地形复杂,气温空间异质性强,如何获取高时空分辨率的山区地表气温数据一直是研究热点与难点。本文选择地形复杂的河北省张家口市作为试验区,基于局部薄盘样条函数对ERA5再分析日均近地表气温（2 m高度）进行空间插值,并利用随机森林算法,结合少量气象站观测气温数据、地形地表参数数据构建日均气温订正模型和气温逐时化模型,实现空间分辨率由0.1 °（约11 km）到30 m的逐时气温降尺度,最后将该模型拓展应用于其他时间与区域,检验本文发展的降尺度方法在没有站点观测数据条件下的时空移植性。结果显示,本文降尺度方法得到的高时空分辨率山区气温数据精度较高,1月均方根误差（RMSE）平均值为2.4 ℃,明显优于气象站点插值结果,且气温相对高低的空间分布更为合理、纹理更加丰富;将该方法应用到其他时间与区域的RMSE平均值分别为2.9 ℃与2.5 ℃,均小于再分析资料直接插值所产生的误差。研究结果总体表明,在气象站点较少甚至没有时,可利用本文方法通过ERA5再分析气温准确获取复杂地形条件下的山区高时空分辨率气温数据。     

混采数据分离中插值与去噪的同步处理

近年来,由于新兴的混合采集观测系统在很大程度上提高了采集效率,因此得到了很多学者和石油公司的青睐.但在实际应用中,这种特殊观测系统的采集质量却受到很多因素的影响.一方面,该观测系统采集到的炮记录会受到相邻炮记录的干扰;另一方面,复杂的采集环境使得地震记录中包含部分空道;另外,采集过程中的场地环境干扰会不可避免地带入随机噪音,它们都会影响采集质量.虽然已有一些学者对这些因素做了相关研究,但都是单独分析,未能综合考虑各种干扰因素.本文基于稀疏约束反演的基本原理,将混合炮数据的分离、缺失道集的插值以及对随机噪音的压制问题整合在一起,通过一步处理同时减小如上三方面因素的不利影响,在改善信噪比的同时极大地提高了地震资料的处理效率.文章利用模拟数据和实际数据对这种新方法进行了验证,均得到了比较满意的效果.     

基于插值原理的检查点技术波场重构与叠前逆时偏移

叠前逆时偏移等基于波场互相关原理的地球物理方法存在极大的计算与存储需求,因此采用合适的波场重构方法显得尤为重要.常规的随机边界法容易产生成像噪声,而有效边界法在三维情况仍难以实现,检查点技术具有内存要求小的特点,但存在较高的重算率,因此本文提出了插值原理的检查点技术波场重构方法.在满足Nyquist采样定理的前提下对相邻检查点间的波场进行规则抽样,将抽样波场作为插值节点,运用多项式插值算法重构任意时刻的波场,从而避免优化检查点技术反复递推造成的计算效率问题.数值实验表明：插值检查点重构算法能有效的恢复波场,其中三次样条插值重构精度最高,而牛顿法插值法计算代价较小适合于快速重构.经Sigsbee模型的叠前逆时偏移证明了插值算法的可行性,并且极大的提高了波场重构的计算效率.三维模型分析得出在增加少量存储的情况下插值重构法的重算率大幅度降低,存储量减少为有效边界法的7.1%,对于三维尺度的叠前逆时偏移有实际意义.     

海域流动点外部扰动引力无奇异计算模型

针对海域重力场变化特征和远程飞行器机动发射保障应用需求,本文分析研究了地球外部空间扰动引力三类传统计算模型的技术特点及其适用性,指出了采用表层法作为海域流动点扰动引力计算模型的合理性及需要解决的关键问题,分析论证了空中扰动引力计算对地面观测数据的分辨率和精度要求,提出通过引入局部积分域恒等式变换、局域泰勒级数展开和非网格点内插方法,消除表层法计算模型积分奇异性固有缺陷的研究思路,进而推出了适合于海域流动点应用的扰动引力无奇异计算模型,较好地满足了全海域和全高度段对局部扰动重力场快速赋值的实际需求.以超高阶全球位模型EGM2008作为标准场,通过数值计算验证了无奇异计算模型的可行性和有效性,在重力场变化比较剧烈的海沟区,该模型的计算精度优于2×10^-5m·s^-2.