简介:本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。
简介:Tilt—depth法可以用于快速反演磁源的上顶埋深,但其无法反演磁源底部深度,同时该方法的反演解过于单一。针对于此,本文在有限厚度台阶磁异常导数公式基础上,推导出了可以同时反演地质体上、下界面埋深的改进型Tilt—depth法计算公式,并采用在Tilt梯度图上选取多特征点进行联合反演的计算模式来提高反演解的可靠性。二维及三维组合模型试验均证实了改进型Tilt—depth法可以有效地反演出磁性体的上顶与下底深度,且上顶埋深的反演精度明显高于常规方法的。最后将改进型Tilt—depth法用于松辽盆地长岭断陷航磁数据反演,其上顶埋深的反演结果与钻孔钻遇火山岩的深度基本一致,证实改进型Tilt—depth法相对于常规方法的反演结果更加准确。
简介:目前在研究地球物理参数的空间变异特性时常常是基于几何各向异性假设,然而实际中几何各向异性假设常常不成立。本文提出了一种求解地球物理参数各向异性的新方法,首先对各方向上的变异值作归一化处理,然后用拟合的方法得到各方向单位步长上的变异性强度,进而实现对全局各向异性的求解。以上海地区合成孔径雷达干涉测量(InSAR)中差分大气延迟样本数据为例,验证了新方法求解各向异性的优越性:新方法的偏差只有6.4%,而基于几何各向异性假设方法的偏差达到了21.2%,而且新方法的运算效率有非常显著的提高。进一步,通过克里金插值交叉验证实验,证明了基于新算法得到的结构函数的插值效果最好,从而从另外一个角度验证了新算法求解的各向异性更准确,更好地表征了区域化变量的空间结构特征。因此,与几何各向异性假设的方法相比,新方法不仅能更准确的描述地球物理参数的各向异性,而且有更高的运算效率,为更准确地估计出所需的地球物理参数奠定了基础。