简介：目前在研究地球物理参数的空间变异特性时常常是基于几何各向异性假设,然而实际中几何各向异性假设常常不成立。本文提出了一种求解地球物理参数各向异性的新方法,首先对各方向上的变异值作归一化处理,然后用拟合的方法得到各方向单位步长上的变异性强度,进而实现对全局各向异性的求解。以上海地区合成孔径雷达干涉测量（InSAR）中差分大气延迟样本数据为例,验证了新方法求解各向异性的优越性：新方法的偏差只有6.4%,而基于几何各向异性假设方法的偏差达到了21.2%,而且新方法的运算效率有非常显著的提高。进一步,通过克里金插值交叉验证实验,证明了基于新算法得到的结构函数的插值效果最好,从而从另外一个角度验证了新算法求解的各向异性更准确,更好地表征了区域化变量的空间结构特征。因此,与几何各向异性假设的方法相比,新方法不仅能更准确的描述地球物理参数的各向异性,而且有更高的运算效率,为更准确地估计出所需的地球物理参数奠定了基础。

简介：有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程，高阶显式有限差分方法通常用来提高求解精度，已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数，发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明：忽略这些小的差分系数能够近似维持有限差分的精度，但是显著减小了计算量。然后，引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性，采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后，给出了针对均匀和非均匀介质模型的弹性波模拟例子，表明了本文方法的优点。