基于非结构化网格的最小二乘逆时偏移

相对于传统的逆时偏移,最小二乘逆时偏移具有更高的成像质量,这种改善是通过迭代反演来获得的,另外其精度与效率依赖于求解波动方程算法的精度与效率.本文给出了基于非结构化网格的最小二乘逆时偏移的方法,该方法能够充分结合最小二乘算法与非结构化网格精细刻画地下界面以及随速度自适应剖分的优点;并采用带补偿的拉普拉斯滤波算法,来消除梯度计算中的低频噪声,从而加速目标函数的收敛速度.通过简单倾斜模型以及复杂Marmousi模型测试,显示了该方法的有效性和潜力.     

基于动态图像变形的PP与PS波层位直接匹配

多分量地震资料的矢量偏移、多波地震资料的联合解释与反演均需要估算纵横波的速度比,实现纵波与转换横波在时间或深度域的匹配.基于DTW,本文实现了一种适用于PP与PS波直接匹配的动态图像变形算法.该算法分为三个部分:首先,使用二阶对称动态规划算法逐样点递归计算PP与PS波走时或深度的误差累积和;其次,在以误差累积和为目标函数的回溯阶段设定变形窗,并在纵横波速比约束的变形窗内递归回溯搜索匹配路径;最后,根据最大相关系数判定准则在匹配路径中确定最佳匹配路径,获得使PP与PS波匹配的拉伸或压缩时移量.利用所获得的拉伸压缩时移量计算纵横波速度比就可以实现PP与PS波之间的匹配.模型与实际陆上多分量地震资料测试结果表明:该方法具有较高的匹配精度,且对于信噪比、相似度较低的多分量地震资料,该方法也能产生较好的匹配效果.     

基于平面波加速的VTI介质最小二乘逆时偏移

地震各向异性集中表现为速度各向异性,势必影响地震波运动学特征.传统声波逆时偏移(RTM)和最小二乘逆时偏移(LSRTM)没有考虑介质各向异性特征,导致反射波不能正确归位、同相轴出现扭曲及寻优速度慢或不收敛等,VTI介质逆时偏移(VTI-RTM)矫正了声波成像的不足,但仍存在低频干扰严重、中深部成像不佳、振幅保持差等缺陷.为此,本文首先实现了VTI介质最小二乘逆时偏移(VTI-LSRTM)方法,为了节省I/O及内存需求并提高效率,进一步引入平面波编码技术,提出了一种基于平面波加速的VTI介质最小二乘逆时偏移(VTI-PLSRTM)策略.在此基础上开展了简单模型及复杂Marmousi模型成像试验,并与标准逆时偏移剖面对比表明:本方法能够校正各向异性造成的相位畸变,且在迭代中自动压制串扰及低频噪声、补偿中深部能量,是一种兼具质量与效率的保幅成像策略;对速度误差的敏感性测试说明该方法需要相对正确的偏移速度及Thomsen参数模型.     

基于卷积完全匹配层的非规则网格时域有限元探地雷达数值模拟

复频移完全匹配层(Complex Frequency-Shifted PML,CFS-PML)在长时间时域计算中对凋落波、倏失波具有好的吸收效果,并被广泛应用于时域有限差分模拟中.而本文采用卷积方法将CFS-PML应用于时域有限元求解GPR波动方程的数值模拟中.论文以TM波为例,推导了基于CPML(Convolutional PML)边界的时域有限元GPR波动方程求解公式,采用Newmark-β方法对时间导数进行离散,有效改善了时域有限元GPR数值计算程序的稳定性.并以狭长模型为例,开展了CPML边界中关键参数m、R和κ的选取实验,通过对比反射误差大小确定了综合最优参数组合.相同时刻UPML与CPML波场快照、3个检测点的反射误差比较,说明CPML较UPML具有更好的吸收效果.最后,采用非规则四边形网格对1个复杂GPR模型进行剖分,应用加载CPML边界条件的FETD程序对该模型进行了正演,得到了二维剖面法、宽角法正演GPR剖面图,说明非规则四边形对复杂模型的良好适应性,基于CPML边界条件的FETD可有效减少边界反射误差,能实现对任意复杂不规则模型的正演模拟.     

基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移

近年来,面向实际应用的TI介质准P波正演模拟与逆时偏移成像技术受到空前的关注.基于常规耦合型传播方程的正演模拟方法不仅存在伪横波及频散假象干扰,而且还遭受模型参数限制(η0)和不稳定影响;而纯qP波方程的推导繁琐,且由于方程中包含拟微分算子造成求解难度大且精度有限.为此,本文首先构建了一种适用于任意TI介质的纯qP波传播算子,然后借助Low-rank分解求取该算子中的空间-波数域矩阵,同时引入Cerjan衰减边界条件来压制边界反射干扰,最终实现了一种间接的纯qP波波场外推方案,并将其成功应用于复杂TI介质正演模拟与逆时偏移成像中.通过开展数值模拟,并与其他方法对比表明:①该方法既避免了纯qP波方程的繁琐推导,又克服了耦合型方程对模型参数的限制;②还彻底消除了残余伪横波噪音及数值频散;③且能适应较大时间或空间步长及高频震源,是一种相对准确且稳定的各向异性纵波正演与成像策略.     

深海6000 m拖曳式瞬变电磁系统及其应用

针对我国西南印度洋合同区热液硫化物矿快速、有效以及便捷地质勘探装备的要求,研发了由甲板控制系统、万米光电复合缆、仪器舱拖体和天线拖体组成的深海6000 m拖曳式瞬变电磁系统.为了便于拖体布放和快速发现异常,选择重叠回线收发装置类型;采用理论数值计算确定了在拖曳高度不大于50 m的前提下,观测1～100ms窗范围内的二次场响应,可以发现近海底深海热液硫化物矿堆;另外,通过不同拖延深度海上试验,研究了拖曳深度对瞬变电磁的影响规律及仪器性能.大洋第30航次第二航段在西南印度洋脊热液区,应用该系统发现了明显的瞬变电磁异常.印证了系统的有效性.     

甘肃地区一维速度模型计算研究

采用2009年1月1日—2014年6月甘肃省境内所有地震的震相数据,分类拟合计算该地区一维速度模型,计算结果显示,利用少量震级较大地震数据来拟合的速度结果近似于用所有地震拟合结果。拟合得出新模型初步结果后,再从所有地震资料中选定台站包围好、震相标识规范的142个地震,利用Hyposat定位方法来计算莫霍面深度,参考其他模型的莫霍面深度,综合分析后在一定的扰动范围内通过Hyposat批处理程序进行迭代定位,对8 640次计算模型做残差对比,最终选取残差最小模型作为新模型,此模型除部分参数有较小的差别外与甘青模型基本一致。     

高烈度区大跨刚构-连续梁桥地震响应分析

为研究大跨刚构-连续梁桥在高烈度地区的地震响应,选取羊记沟左线大桥为工程背景,利用Midas软件建立多个有限元模型,考虑纵、横桥向的水平地震波输入,比较反应谱与时程分析结果,获得结构动力响应特点,为大跨刚构-连续梁桥的抗震设计提供参考依据。结果表明:刚构-连续体系仅单墩刚构时不宜采用反应谱分析,高烈度地区桥梁进行地震响应时程分析时,选用与场地适应的地震波的同时,应考虑采用本地地震参数转化人工波作为地震激励源。桥墩不等高时,可通过改变墩截面形状及其与主梁连接方式调节内力分布,避免矮墩刚度大导致的内力集中。     

滇西南普洱强震群地区地壳形变的数值模拟

滇西南普洱一带地质构造复杂,中-强地震成群活动特征明显,地壳的形变特征与形变机制引起了地学界的广泛关注和深入研究。在参考、借鉴以往研究成果的基础上,本文构建了普洱一带包含活动断裂的二维板有限元模型,以GPS测量得到的块体年运动量作为模型区域的边界约束,模拟计算获得了普洱地区的构造形变空间分布,以NEE、NNW 2个方向、从数值模拟角度对普洱地区的水平形变特征和区域地壳垂直形变场进行了分析和探讨。结果表明:在青藏高原东南缘扩展、高原物质重力滑塌型南东向移动的动力学背景下,研究区内NEE向断裂的活动带动牵引NNW向断裂共轭运动,隆升性质的垂直形变场与60年的大区域水准测量结果相吻合,沿NNW向形变场的分布无规律,应是区域应力累积和震后调整双重作用的结果。沿NEE向形变场呈现了较明显的西部大于东部的分布图像,特别是西北部高值区对应了2014年景谷6.6级强震,最后对研究区的区域动力学问题进行了初步探讨。     

基于滑动窗口导纳技术反演西太平洋中部岩石圈的有效弹性厚度

西太平洋中部地区是西太平洋板块边缘沟-弧-盆体系构造演化的关键区域,其地质特征与构造演化一直是地学家关注的焦点问题之一。开展岩石圈有效弹性厚度的研究对于认识该区域的形成演化具有重要的科学意义。本文采用滑动窗口导纳技术,并在挠曲模型中考虑了表面荷载和内部荷载同时存在的情况,计算得到该区域的岩石圈有效弹性厚度(Te)。计算结果显示,研究区的Te值整体上为0~50 km,其变化基本上与构造单元相吻合,且与主要的构造边界密切相关。除海底火山具有相对较小的Te值(15~20 km)外,太平洋板块整体上具有较强的岩石圈强度(25~30 km)。马里亚纳海沟和菲律宾海沟的岩石圈强度从外隆起到海沟方向表现为明显的减弱,表明岩石圈由外隆起向海沟发生了弱化。帕里西维拉海盆西部相较于东部具有较弱的岩石圈强度,这可能与海盆的非对称扩张有关。卡罗琳板块的岩石圈整体上表现为相对均一的低Te值特征(<15 km)。欧里皮克海隆、卡罗琳海岭和索罗尔海槽的Te值为3 km,这可能是强烈的火山作用所导致的结果。