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

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

     

巴颜喀拉块体地壳结构多样性探测

青藏高原内部地壳岩性的改造、岩性随深度变化及形变构造是探索研究地壳增厚、物质运动问题的关键.巴颜喀拉块体位于青藏高原中北部,地域广袤,通过对块体内中、东部不同区域的深地震广角反射/折射震相的综合分析,利用反射率理论地震图方法对不同性质震相走时及振幅特征进行细致的模拟计算,进一步研究巴颜喀拉块体内部不同区域地壳精细结构.结果显示：巴颜喀拉块体地壳厚度50~60 km、整体向西逐渐增厚,结晶地壳平均速度6.07~6.18 km·s-1、岩层速度大幅降低,壳内多强反射界面结构、但不同区域差异明显;东部若尔盖盆地地壳介质速度整体低速、壳幔边界较为清晰;中部玉树—玛多段下部地壳发现约6.8 km·s-1的"高速度"介质结构,壳幔边界不清、被改造为2~4 km厚的高速度梯度层,显示了巴颜喀拉块体内部地壳增厚、介质岩性结构被改造的差异性.地壳内部多组强反射、低视速度走时震相揭示了介质岩性的低速破碎、弱化蠕变以及可能的壳内解耦构造.局部地区下地壳的高视速度震相特征显示了青藏高原地壳改造增厚大背景下可能存在稳定的"原始地壳"结构残留或是与上地幔物质的浸入交流.巴颜喀拉块体内不同区域地壳增厚、岩性结构、结晶基底及壳幔边界性质被改造的多样性为深入认识青藏高原地壳形变及动力学过程带来新的启迪.     

四川盆地深部地壳结构——深地震反射剖面探测

四川盆地位于扬子地块的西北部,被褶皱构造带所围绕,受周缘构造带的侧向挤压作用,盆地卷入了多期次和多边界的构造变形,为开展盆山耦合作用及多边界、多期次构造叠加与复合关系的研究提供了不可多得的理想野外实验室.为揭示四川盆地地壳结构,本文通过对3条不同时间采集的深地震反射剖面数据进行拼接联线处理,获得跨越四川盆地的330 km深地震反射偏移成果剖面,揭示了四川盆地地壳上地幔细结构：沉积层从西北向东南逐渐变薄,在龙门山前沉积层厚度超过15 km,在华蓥山下沉积层减薄到~8 km,且褶皱变形形成华蓥山薄皮褶皱冲断带;莫霍面出现在13~15 s（双程走时）,埋深约40~45 km;并发现从下地壳延伸至地幔的东南向的倾斜反射,从13 s向下延伸至18 s,结合四川盆地及其周边地区其他地球物理和地球化学花岗岩同位素年龄等资料,我们认为这些倾斜反射层是扬子克拉通地台西北缘发生的新元古代俯冲的遗迹.

     

2017年四川九寨沟MS7.0地震前区域重力场变化

利用南北地震带2014-2017年期间的流动重力观测资料,系统分析了区域重力场变化及其与2017年8月8日四川九寨沟7.0级地震发生的关系.结果表明：① 区域重力场异常变化与北西西向塔藏断裂和南北向岷江断裂带在空间上关系密切,反映了沿控震断裂在2013-2017年期间发生了引起地表重力变化效应的地壳变形和构造活动.② 九寨沟地震前,测区内出现了大空间范围的区域性重力异常,而震源区附近产生了局部重力异常,沿塔藏断裂带形成了重力变化高梯度带,其中,甘肃玛曲、迭部、青海河南蒙古族自治县、四川若尔盖、九寨沟一带重力差异变化达100×10^-8m·s^-2以上;这些可能反映九寨沟地震前,区域及震源区附近均产生与该地震孕育、发生有关的构造运动或应力增强作用.③ 九寨沟地震震中位于重力差异运动剧烈的鞍部等值线附近,与断裂走向基本一致的重力变化高梯度带零值线上,这一观测事实进一步佐证了重力场动态变化图像对强震地点预测具有重要的指示意义.

     

一种数字三分量地震微测井仪器设计

针对当前地震微测井仪器在施工和数据采集方面的难题,从提高数据采集质量及仪器功能性方面提出了几种改进方法.主要利用MEMS加速度传感器结合多路并行处理采集方法设计了数字检波器,并集成姿态检测技术直接实现水平分量数据的校正功能,保证系统采集的质量;提出了Ethernet-RS485分布式总线自动控制技术,实现各级联探管的随机组合,同时采用电动方式控制探管操作,提高系统的集成化程度.文中详细分析了利用姿态检测技术实现水平分量数据校正的方法,解决了级联探管随机组合控制和震源激发方位的判定问题.仪器单节探管支持三分量并行采集,有效分辨率达18位,参数可远程配置,同时具备快速显示探管姿态的功能,且支持远程控制,非常适于复杂环境下的工程勘探,整个系统实现全数字化操作且仅需四芯电缆互连,从而缩小了设备的体积和重量,降低施工布线的复杂性.

     

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明：新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.

     

汶川地震断裂带科学深钻WFSD-3附近上地壳S波分裂特征

为配合汶川地震断裂带科学钻探工程（WFSD）,中国地震局地球物理研究所在四川省绵竹市天池乡和灌县—安县断裂附近分别架设了15套微震仪器和17个短周期地震台.基于WFSD-3附近的微震、短周期和区域台网的固定台站记录的近震数据,通过横波窗内S波分裂计算,得到其上地壳各向异性参数,即快波偏振方向和慢波的时间延迟,并分析了研究区的上地壳各向异性特征.结果显示,研究区大部分区域的快波偏振方向为NE向,与龙门山断裂带走向一致,但在研究区微震台阵布设小区域内,快波偏振方向表现出东西分区特征,东部为NE向,西部为NW向.上地壳各向异性主要是受到岩层中随应力分布排列的微裂隙和岩石或矿物结构的影响,研究区内快波偏振方向主要表现为NE方向,与断裂走向一致,反映了研究区上地壳各向异性主要受控于结构控制的各向异性,局部区域的快波偏振方向为NW向,与区域最大主压应力方向一致,说明区域应力场对研究区上地壳各向异性也有影响.通过分析微震台阵的归一化时间延迟随时间的变化情况可以反映区域应力场的变化情况.微震台阵的慢S波时间延迟在2012、2013年较为离散,在2014年有收敛的趋势,反映了强震后区域应力场逐渐稳定的趋势.

     

随采地震监测数据采集控制软件开发

随采地震能够对工作面前方地质异常体进行连续探测和实时预报,成为近几年的研究热点,但是目前还没有能够在煤矿井下开展随采地震长期连续监测的装备及配套软件。为了解决这个问题,基于Microsoft Foundation Classes（MFC）开发框架,开发了一套随采地震监测数据采集软件,在室内、野外进行了为期3个月的联调测试,并且在贵州岩脚煤矿与井下随采地震监测设备开展了为期3个月的全面试运行。测试表明,软件实现了随采地震信号的高效采集、完全存储和处理软件的实时通信功能,具有运行稳定、操作便捷、处理高效、便于维护、无人值守等优点。     

南海海盆区莫霍面分布规律及其对深部钻探的意义

钻遇莫霍面是人类一直以来的梦想。深海海底是地球上离莫霍面最近的地方,目前有研究推测南海是世界上莫霍面深度最浅的海域之一,但缺乏足够的直接证据。深反射地震探测可以直接揭示岩石圈的构造形态,是莫霍面探测的重要手段。本文基于长达15000 km的深反射多道地震剖面的解释、处理、制图和分析,结合前人的研究,形成了南海海盆区莫霍面反射特征和空间分布的初步认识。① 南海东部次海盆南部早期经历了较快速扩张,岩浆供应充足,受扩张停止后岩浆活动影响较小,基底平坦,地质构造相对简单,同时洋壳地震速度结构不存在异常,且有较强的广角莫霍面反射波和可识别的地幔顶部折射波,具备莫霍面钻探的基本条件。② 南海海盆不同区域的莫霍面反射强度存在较大差异。其中东部次海盆莫霍面反射最为强烈且清晰,西北次海盆次之,西南次海盆仅有零星出现的清晰莫霍面反射且可信度不高。③ 识别南海海盆区莫霍面地震反射长度超过3500 km,首次形成了海盆区深度域莫霍面地震反射空间分布图。与重力反演的莫霍面深度相比,利用深反射多道地震计算的莫霍面深度细节更为丰富,并且可以在垂向上清晰刻画莫霍面的结构。整体上,南海海盆区莫霍面地震反射强烈和可信度高的区域中,深度较浅的区域之一是东部次海盆南部,最浅处仅约9. 5 km,其中水深4. 01 km,洋壳厚度仅5. 54 km。综合判断,东部次海盆南部是南海重要的莫霍面钻探备选区,这对南海莫霍面钻探选址具有重要意义。     

带有横向约束的全局优化波阻抗反演方法及应用

地震波阻抗反演是一种基于地震数据得到波阻抗参数的非线性优化问题。不依赖于目标函数梯度信息的全局优化算法是求解地震波阻抗反演问题的有效方法。然而这类方法采用逐道反演的策略,忽略了相邻地震道的空间相关性,导致反演结果的横向连续性差。针对该问题,提出了一种融入旁道最优解的模型空间初始化方法约束波阻抗反演的搜索空间范围,以改善反演结果的横向连续性,进而提出了一种带有横向约束的多组变异差分进化地震波阻抗反演方法。合成地震记录算例表明,该方法相比于传统方法不仅具有更快的收敛速度,而且反演结果具有更好的横向连续性。此外,该方法被应用于胜利油田某区块的储层波阻抗参数反演,反演结果与测井资料有很好的一致性,且有效地刻画了储层砂岩厚度。