交错网格任意阶导数有限差分格式及差分系数推导

交错网格有限差分算法以其高效、精确、实用等优点在地震波数值模拟中得到广泛应用。目前交错网格有限差分的精度已达到时间4阶、空间2N阶;然而在求空间三次导数时,差分格式实际上并未达到所谓的2N阶精度,而是采用了低阶的差分格式及差分系数,这样有利于提高大尺度空间正演时的计算效率;但从计算精度的角度考虑,有必要推导出准确的满足2N阶精度的交错网格有限差分格式及差分系数,以得到更高精度的正演结果。笔者利用Taylor公式展开首次推导出了可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而完善了常规高精度交错网格有限差分算法。采用新推导的交错网格有限差分格式得到的正演波形与解析解进行了对比,证明了新推导的差分格式的正确性,并与常规差分格式的正演波形进行了比较,结果显示,新推导出的交错网格有限差分格式模拟结果稳定性好,精度更高。     

沉积盆地地区地壳结构估计——预测反褶积方法消除接收函数多次波混响

接收函数方法被广泛地应用于地壳上地幔结构的研究中,H-κ叠加方法是其中最常用的方法之一.对于布设在基岩区台站计算的接收函数,H-κ叠加方法可以准确地估计台站下方地壳厚度和平均波速比,但是对于沉积盆地地区计算的接收函数,由于低速沉积层内会产生多次波混响,干扰甚至覆盖接收函数中莫霍面的转换波和多次波震相,从而影响H-κ叠加结果的准确性.为准确估计沉积盆地地区地壳结构,本文提出使用预测反褶积方法去除接收函数中低速沉积层内多次波混响,其中预测步长由接收函数归一化自相关函数获得,物理意义为沉积层内S波双程走时.合成接收函数和实测接收函数试验表明,本文提出的预测反褶积方法可以有效地去除沉积层多次波混响,并结合改进的H-κ叠加方法可以准确地估计沉积层下覆地壳厚度和平均波速比.相比于其他去除接收函数多次波混响的方法,本文提出的预测反褶积方法具有参数设定简单、运算量小、震相幅值较大等特点,适用于大批量数据处理.     

金川矿区及周边地球物理特征和背景噪声成像的应用

为探明甘肃金川矿区地下深部矿产,利用矿区及周边共65个固定台站的连续地震背景噪声数据的Z分量,基于射线追踪的面波频散直接反演方法获得了金川矿区及周边地下0.9～5.4 km的三维S波速度结构.结合金川矿区及周边地球物理特征(重力及航磁异常)等,对反演成果图分析表明:金川矿区在地下浅部1.8 km以浅,表现为高速特征;而...     

黏弹各向异性介质中波的反射与透射问题分析

黏弹各向异性介质中传播不均匀波,其反射、透射模式不仅与介质分界面两侧速度对比有关,还与品质因子Q的对比有关. 用伪谱技术模拟黏弹各向异性介质分界面上波的反射、透射,并与弹性各向异性介质、黏弹各向同性介质和弹性各向同性介质的模拟结果做比较. 计算平面波的反射、透射系数,分析介质的黏弹性和各向异性对反射、透射系数的影响. 数值模拟了一个三层介质模型中的波场,分析两个分界面上产生的反射波的特征. 黏弹各向异性介质中,qS波比qP波衰减程度大.     

瑞雷波测试技术在岩土工程中的应用研究

本次研究利用瑞雷波测试技术对北京地区岩土工程目标体进行了研究.针对评估地基承载力方面的问题提出利用瑞雷波测试成果评价地基承载力及均匀性可以采用三种方法:直接应用瑞雷波波速评估地基的强度、在时域曲线判读群速度评价地基土的均匀性、应用瑞雷波测试综合信息评价地基强度.通过不同地区的工程实例分析,总结出一些特征及关系,对地基强度进行量化评估提供了一个合理、可靠、简便的途径.     

波动方程参数反演方法研究进展

根据波场信息重建波动方程系数项是地球物理反演中一个重要研究方向.波动方程参数反演方法按照求解方式的不同,可以分为直接反演方法、迭代反演方法、搜索类反演方法、混合优化反演方法和联合反演方法.目前发展的参数反演方法多数都处于理论计算阶段.研究能求解实际地震资料的参数反演方法具有重要意义.为取得可靠的反演结果,信息的转化开发...     

全极化探地雷达正演模拟

全极化探地雷达不但能获取当前探地雷达使用的共极化信息,而且能使用交叉极化信息,进而更好地分析目标属性.本文在时域有限差分法的基础上,构建了全极化探地雷达的正演模拟方法,通过水平正交双方向同时接收获取全极化信息.标准目标体金属板和角反射器的正演模拟与实验室物理模型测试一致的共极化和交叉极化响应信息,证实了该正演模拟方法的...     

大兴安岭西北部中新生代盆地群基底电性分带特征研究

位于大兴安岭西北部的中新生代盆地群(海拉尔、根河、漠河),其构造受到大兴安岭断裂、德尔布干断裂的控制,西北方向的蒙古-鄂霍茨克缝合带、南部西拉木伦河-延吉缝合带甚至更远的西太平洋板块运动、印度板块运动以及黑龙江中西部微板块间拼合等区域构造应力场叠加作用在该盆地群基底产生了复杂的深部构造特征.本文利用沿盆地群实施的4条(...     

基于Matcom 平台的VC + + 与Matlab 混合编程地震正演模拟

正演模拟是地震勘探的重要方法之一。为解决基于Matcom 4. 5 的VC + + 与Matlab 接口程序的地震正演模拟问题,通过用Matcom 4. 5 与VC + + 结合编程的新方法实现EDA 介质中的地震波正演模拟的实例,结果显示Matcom 实现VC + + 与Matlab 混合编程并提高了编程效率。     

Pearson 相关系数法快慢横波波场分离

横波分裂是各向异性介质的重要特征,当横波或转换波穿过各向异性介质到达地面时,地面三分量检波器的x 分量和y 分量接收到的地震记录中都会同时存在快横波和慢横波。将快横波和慢横波进行分离,进而计算介质的各向异性参数是多分量数据处理中重要的一步。将数学中的Pearson 相关系数引入到多分量地震勘探中,提出了Pearson 相关系数法进行旋转角度识别,进而分离快、慢横波波场。相比于传统的互相关法,Pearson 相关系数法从精度、抗噪性能和计算效率上都有提高。