大别-苏鲁及其邻近地区基于背景 噪声的勒夫波群速度成像

大别-苏鲁造山带是中国大陆东部地区最重要的构造之一. 为了研究该地区的地壳上地幔速度结构,本文收集了国家数字地震台网和中国区域地震台网的山东、河南、安徽、江苏和湖北等省的144个宽频带地震台连续两年(2009年5月—2011年5月)的水平向地震记录(E分量和N分量)数据,首先对台站对之间E-E,E-N,N-N和N-E分量进行互相关,然后分别对这4个互相关分量采用相位权重叠加法进行叠加,最后旋转到横向分量(T-T)获得勒夫波经验格林函数(EGF);用频时分析(FTAN)方法获得4000余条勒夫波群速度频散曲线,并进一步反演得到了周期为6—40s的勒夫波群速度分布图.结果表明,周期为6—10s的勒夫波群速度分布与地表构造特征相吻合.大别造山带、苏鲁造山带、湖北西部隆起均表现为高速;华北盆地发育,表现为大面积的低速;江汉盆地、南襄盆地、合肥盆地等因其规模不同而显示不同程度的低速.在周期为6—30s的勒夫波群速度分布图上,大别和苏鲁地区均显示高速,已有的研究结果中地壳的低速并没有得到反映.其原因一方面可能与勒夫波群速度纵向分辨不高有关,另一方面高压变质岩深度分布可能比已有研究结果给出的要深. 6—30 s的分布图上郯庐断裂带及其邻近地区表现为不同程度的高速,可能与该地区白垩纪以来处于拉张构造体制,地幔物质受到扰动,造成物质上涌有关.     

塔克拉玛干沙漠腹地风蚀起沙特征观测研究

利用塔克拉玛干沙漠腹地塔中气象站2009—2018年的地面观测资料,分析塔中地区的风蚀起沙特征,结果表明:(1)塔中地区临界起沙风速呈明显的季节变化特征,夏季>春季>秋季>冬季,临界起沙风速介于4.47～4.92 m/s;(2)春、夏季是风蚀事件的多发季节,春、夏季沙尘水平通量的平均值分别为2638.9、3298.9 ...     

上海TWP3型边界层风廓线雷达探测性能评估

基于2018年7月—2019年6月上海宝山移动式风廓线雷达（mobile wind profiler radar,MWPR）资料,从获取率、水平风误差和上海6类典型天气条件下测风数据的可靠性三方面,对其探测性能进行了评估与分析。结果发现：（1）获取率受湍流强度影响最大,受大气温度、湿度影响次之。获取率具有较明显的日变化和季节变化特征,边界层内夏秋季获取率高于冬春季,且有午后获取率高、深夜和清晨低的特点;自由大气中冬春季获取率高于夏秋季。以获取率80%为其业务准入标准,边界层内MWPR的获取率达标。（2）相比L波段探空测风数据,MWPR整层的水平风速均偏小,低模态时有34.4%的数据位于偏小2.0~6.0 m·s^-1区间内;风向偏差超过±15°以上的水平风向数据占比为48.9%。（3）定性分析结果显示,MWPR探测到的垂直速度的指向性是合理的。雨强较大时,MWPR探测的垂直速度实际上是降水粒子的下降速度与大气湍流速度之和。MWPR的水平风数据在台风天、雪天等大气均一性强的时段内可靠性高。总体上,上海宝山MWPR的探测性能基本能满足日常监测服务之需,但其测风精度与业务准入标准相比尚有一定差距。

     

海底地震仪远震记录接收函数反演——以南海西南次海盆为例

由于海底环境和海底地震仪(OBS)结构的特殊性,用OBS远震记录进行接收函数岩石圈反演研究因为存在一定的困难,所以还很少见.在深入分析问题的基础上,以国产I-4C型宽频带OBS在南海西南次海盆记录的天然地震为实例,我们将傅里叶变换和小波变换相结合以压制海底地震仪记录中的非平稳干扰,获得了信噪比较高、震相清晰的地震记录,进而成功开展了远震记录的岩石圈结构接收函数反演.主要结论是:(1)OBS接收函数的求取是可行的,关键是压制非平稳干扰.(2)西南次海盆的Moho面埋深为海底下10~12 km(地壳厚6~8 km),沉积物厚度为1~2 km,浅部地壳存在低速区,与沉积物和海底扩张停止后的岩浆喷发产生的岩石碎屑和裂隙有关.(3)在扩张脊中央Moho面上方6~12 km存在S波低速区,推测扩张中心可能存在下地壳熔融或岩浆房,在17~30 km区间S波速度呈负梯度,我们认为扩张中心更深的地方存在热物质的供给.     

利用接收函数反演青藏高原西部地壳S波速度结构

相对于宽阔的腹地,青藏高原西部南北向宽度仅约600 km,却记录了印度和欧亚板块汇聚的深部过程及其响应.本文用22台宽频带流动地震台站在西缘构建了一条南北向探测剖面(~80°E,TW-80试验).利用接收函数反演剖面下方S波速度结构,综合西部已有的宽频带探测结果,分析认为:印度板块向北俯冲可能已到达班公湖-怒江缝合带附近,俯冲过程中下地壳发生榴辉岩化;喀拉昆仑断裂带、班公湖-怒江缝合带、阿尔金断裂带均为切穿地壳的深断裂,莫霍面发生错断;喀拉昆仑断裂带和龙木错断裂带之间的中上地壳没有发现连续的S波低速体,说明可能缺乏解耦层,支持青藏高原西部地壳为整体缩短增厚模式.     

基于微动方法研究五大连池火山区尾山火山锥浅层剪切波速度结构

为了研究五大连池火山区尾山火山锥浅层三维波速结构特征,在尾山火山锥附近区域布设了无线地震检波器密集台阵,记录连续地震背景噪声波形数据.基于微动方法（拓展空间自相关方法）提取了台站间2~5 Hz频率范围的Rayleigh面波相速度频散曲线.利用面波层析成像方法反演获得2~5 Hz Rayleigh面波二维相速度图像,基于每一个网格节点的频散曲线,进一步反演获得了尾山火山锥附近区域地表至700 m深度的三维剪切波速度结构.成像结果显示：在0~150 m较浅深度,靠近尾山火山锥区域显示为相对高速异常,远离火山锥区域则显示为相对低速异常.而至150~700 m较深深度,波速异常特征与浅部相反,靠近尾山火山锥的区域显示为相对低速异常,远离火山锥的区域显示为相对高速异常.在远离尾山火山锥区域,浅层的相对低速异常可能与松散沉积层有关,深部的高速异常则反映了结晶变质岩的影响.在靠近尾山火山锥区域,浅部的相对高速异常应该反映了出露地表的玄武岩,而深部的相对低速异常则可能反映了火山通道周围广泛发育的破碎裂隙结构及其火山喷发后孔隙流体填充的影响.     

基于改进的散射积分算法的初至波走时层析

初至波走时层析是获取近地表速度结构的一种常用方法.随着采集技术的不断发展,可使用的数据量迅速增多,传统的基于射线追踪和解方程组的地震走时层析成像方法面临着内存占用大、方程求解不稳定等问题.为了解决这些问题,本文基于前人在波形反演研究中提出的一种改进的散射积分算法,提出了一种预条件最速下降法初至波走时层析.该方法无需存储核函数矩阵与Hessian矩阵即可方便地实现目标函数梯度的计算与预条件,且该方法计算效率高、求解稳定、易于并行.数值实验结果表明,该方法可以获得与传统方法精度相当的反演结果,但所占用的内存大幅减小.     

哈拉湖地区天然气水合物地震探测技术试验

哈拉湖位于青藏高原,在该区开展的地质调查、音频大地电磁测深和地球化学调查结果表明,该区具有良好的天然气水合物找矿前景。为寻找天然气水合物,在该区开展了高精度反射地震试验。反射地震采用1200道接收,道间距2 m,偏移距1 m,炮间距8 m,覆盖次数150次,排列中间激发的观测系统,每道采用6个60 Hz检波器单点组合接收;采样间隔0.5 ms,记录长度2 s,宽频带采集;激发震源为大型车载可控震源。采用该工作方法得到的地震剖面信噪比和分辨率较高,构造形态特征明显。根据试验区地形地貌和地震剖面上反映的永冻土层厚度、断裂构造和高速层分布,结合其他资料,确定了验证孔位。     

羌塘盆地中部地区地壳S波速度结构及构造意义

羌塘盆地是我国最大的海相盆地.本文根据在羌塘盆地内布设的27个宽频带地震观测台站记录的远震波形数据,利用非线性复谱比反演算法得到各台站下方100 km深度范围内S波速度结构.结果表明,羌塘地区Moho深度较为平缓,平均深度为61 km;北羌塘地壳内低速层广泛分布;北羌塘具有两个较大的沉积盆地,龙尾错和白滩湖坳陷,沉积厚度分别有10 km和15 km.尽管北羌塘下地壳受到强烈的新生代火山岩作用改造,但是这种深部岩浆热作用会加速烃源岩中有机质的热演化历程,北羌塘两个盆地具有很好的油气前景.与北羌塘低速层分布相比,南羌塘下低速层更深,可能与班公怒江洋于中生代的俯冲消减及拉萨地体北向俯冲有关.

     

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

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