由于沉积环境、构造运动和上覆压力等地质因素，地下介质的空间结构通常表现出各向异性特征，传统的基于全局空间范围的克里金插值方法难以准确描述地下各空间点的各向异性空间结构特征。本文由Marmousi2波阻抗模型生成合成地震数据，并从中估计随机介质的统计特征参数（图1-2），表征地下介质中各空间点的各向异性空间结构特征，然后计算地下介质在不同方向上的空间相关长度，得到各向异性空间协方差函数，进而得到各向异性克里金权重系数，最后得到最终的各向异性空间结构克里金插值结果（图3）。本文方法得到的克里金插值结果更可靠，更符合地下介质的实际各向异性空间结构特征。
图1. 模型数据 (a)Marmousi2 波阻抗部分模型数据；(b)模型合成地震记录
图2. 随机介质的自相关参数 (a)横向自相关长度；(b)纵向自相关长度；(c)自相关角度
图3. 三口伪井的克里金插值 (a)-(c) 基于使用全局空间范围的各向同性克里格插值方法，分别显示克里金插值结果、方差和克里金插值结果与模型数据之间的误差；(d)-(f) 基于使用全局空间范围的各向异性克里金插值方法；(g)-(i) 基于随机介质统计特征参数的各向异性克里金插值方法
随机介质的统计特征参数可以与地下介质各空间点的各向异性空间结构特征一一对应，本文方法可以得到更符合地下介质实际空间结构特征、误差更小、插值不确定性更低的插值结果。
研究成果发表在地球物理国际权威期刊IEEE Transactions on Geoscience and Remote Sensing上，论文第一作者为太阳成集团tyc151cc刘龙东，通讯作者为太阳成集团tyc151cc王保丽副教授。该研究得到了国家自然科学基金42174144和U23B6010的资助。
论文信息：Liu L, Wang B, Lin Y, Zhang G, Yin X, Zhu Z. Anisotropic Spatial Structure Kriging Interpolation Method Based on Statistical Characteristic Parameters of Random Media. IEEE Transactions on Geoscience and Remote Sensing, 2025, 63, 3578443.https://doi.org/10.1109/TGRS.2025.3578443.
