基于地震道匹配的地层空间转换计算方法

地震数据是一个空间变化关系的数据载体,反映了地质体空间中的变化特征.在已知地震数据的前提下,如何刻画储层空间变化规律是精细储层研究的一个重要问题.笔者通过对地层空间转换函数进行研究,提出了一种基于地震道匹配求取地层空间转换函数的方法.该方法采用地震道多道匹配求取加权系数,使用多点的褶积加权的地层空间转换函数取代点对点的地层空间转换函数,从而准确反映储层物性变化,减小误差.基于Q油田的实际资料,通过模型测试及实际应用分析,结果显示:相对于传统方法,本方法建立的模型分辨率较高,包含信息较丰富,在有砂体分布的地方孔隙度值高,而且高值分布广,与测井孔隙度的分布基本一致.本方法作为一种地震信息驱动的确定性建模方法,可以更精确地进行储层表征,对地质建模有重要意义.     

基于高斯束理论的有限频核函数计算

相对于常规射线层析速度建模,基于波动理论的层析速度建模考虑了波的带限特性,反演分辨率更高.波动理论层析的核心在于波路径(有限频核函数)的计算.文中详细介绍了一种基于高斯束算子计算有限频核函数的方法,分析了初始束宽度和高斯束出射角度间隔对计算精度的影响;并针对高斯束近源处误差较大的缺陷,提出了改进的束参数以提高近源精度;详细分析了初始束宽度和角度间隔对改进高斯束方法的影响及改进高斯束的聚焦特性;数值算例验证了该方法在缓变介质中计算有限频核函数的可行性,计算效率较高且可处理回折波的核函数.     

基于非结构化有限元的三维井地电阻率法约束反演

电磁探测反演是典型的不适定问题,易造成反演结果的多解性,不适定性是反演自身固有的特征,没有求解的附加信息这一本质困难是很难克服的,解决该问题的有效方法是研究约束反演。本文采用目前较为主流的高斯牛顿—共轭梯度法(GN-CG),在反演目标函数中直接施加约束条件,将介质电阻率的取值范围作为先验信息和约束条件以外点罚函数法的方式引入到反演目标函数中,与常规三维电阻率反演目标函数相比,增加了不等式约束项的目标函数,理论上可以压制反演的多解性。通过多种理论模型的测试结果表明,本文基于不等式约束的三维井地电阻率反演算法有效地改善了反演结果的精度,以惩罚函数法施加不等式约束条件的方式是现实可行及有效的。     

基于SPEI的渭河流域干旱时空变化特征分析

以标准化降水蒸散指数(SPEI)作为评估指标,基于渭河流域28个气象站点1961—2017年实测降水量和气温数据,采用Mann-Kendall(M-K)趋势检验、经验正交函数以及小波变换等方法分析渭河流域干旱时空变化特征,并研究渭河流域干旱与6种大尺度气候因子之间的相关关系,进一步探讨主要气候因子对流域干旱时空分布特征的潜在影响。研究表明:渭河流域在1961—2017年间整体呈现出变旱的趋势。通过经验正交函数分解,渭河流域干旱分布场主要有3种典型模态类型,分别为全局型、西北—东南反向分布型以及东—西反向分布型;同时,大尺度气候因子南方涛动指数SOI与流域干旱分布场具有更好的相关关系,对该区域内干旱变化有较强的影响。     

顾及垂直递减率函数的中国区域大气加权平均温度模型

大气加权平均温度(T_m)是全球导航卫星系统(GNSS)水汽监测的关键参数。针对中国区域地形起伏较大的特点,本文构建了顾及精细季节变化的T_m垂直递减率函数模型,在此基础上,利用2007—2014年的Global Geodetic Observing System (GGOS) atmosphere格网数据建立了中国区域的T_m格网新模型(简称为CT_m模型)。以2015年GGOS格网数据和无线电探空资料为参考值,对CT_m模型进行精度检验,并与常用的Bevis公式和GPT2w模型进行比较分析。结果表明:①以GGOS格网数据为参考值,CT_m模型的年均偏差和均方根误差(RMS)分别为-0.52 K和3.28 K,相比于GPT2w-5和GPT2w-1模型,精度(RMS值)分别提高了27%和13%;②以探空数据为参考值,CT_m模型的年均偏差和RMS误差分别为0.26 K和3.75 K,相对于GPT2w-5和GPT2w-1模型,精度分别提高了21%和1...     

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

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

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

由于海底环境和海底地震仪(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波速度呈负梯度,我们认为扩张中心更深的地方存在热物质的供给.     

黄土高原表土磁化率与气候要素的定量关系研究

黄土高原地区已建立多个表土低频磁化率-气候转换函数,为黄土古气候定量重建提供了关键方法。但是表土磁化率变化的气候控制因素的系统研究尚未开展; 同时,低频磁化率部分受到沉积作用的影响,需要寻找气候意义更明确的指标。我们系统采集了黄土高原及周边地区的表土样品,运用相关、回归和因子分析等方法,研究了不同气候要素及其季节变化对表土磁化率和频率磁化率的控制作用。结果显示,湿度是控制黄土高原表土磁化率的主要因素,温度的影响相对较小; 降水的季节分配也有影响,月降水变率较小的地区磁化率较高。在此基础上,选择反映成壤磁性颗粒组分的频率磁化率,建立了频率磁化率-年均降水量的转换函数,为黄土古气候定量估算提供了新的途径。     

利用接收函数方法研究中国东北东南部地区不同构造体的地壳特征

基于我们布设的探测深俯冲的中国东北地震台阵NECsaids台阵和固定地震台长时段的观测记录及NECESSArray流动台阵共计152个台站数据提取得到的33752条P波接收函数,采用H-κ叠加分析和共转换点（CCP）叠加成像等方法进行统一分析处理,并汇集他人接收函数研究结果得到中国东北东南部地区迄今为止最高分辨率的地壳厚度和平均波速比分布图像.对中国东北东南部地区不同构造体的地壳特征综合分析研究表明：研究区不同陆块的地壳属性存在明显差别,张广才岭地块中南部的地壳厚度和波速比与华北克拉通东北缘相当,地壳厚度同地表地形之间显示有明显的正相关关系;松辽地块东南缘地壳最薄、波速比最高,地壳厚度同壳内波速比之间显示出明显的负相关关系;兴凯地块西部地区的地壳结构表现为稍厚的地壳厚度和研究区内最低的壳内波速比,其地壳厚度同壳内波速比之间亦显示出明显的负相关关系;佳木斯地块西南缘在具有"正常"的壳内波速比同时地壳最厚.研究区内的郯庐断裂带北延段在切穿其下Moho面的同时表现出南北分段的特征：北段（44.4°N—47°N）两分支之下的Moho面整体下凹,而南段（41.5°N—43.3°N）两分支之下的Moho面则整体上隆.长白山天池火山下方表现为Moho面下凹沉落及高壳内波速比特征,推测其壳内岩浆囊很可能存在于火山口东北隅至少10 km的范围内.     

内蒙古阿巴嘎地区地壳方位各向异性研究

本文利用架设在内蒙古阿巴嘎地区38个宽频带地震台站记录到的远震数据,通过拟合P波接收函数径向Pms转换波到时和叠加不同方位切向分量,确定了地壳各向异性参数.结果表明,大部分台站Pms延迟时间在0.35 s左右;而少数台站时差较大,推测可能受到索伦缝合带附近地壳残留倾斜界面影响.各向异性快波方向变化范围在N95°E到N180°E之间,平均为N130.6°E±19.1°,推测中下地壳矿物在ENE-WSW向区域主压应力作用下发生晶格定向排列可能是导致地壳各向异性的主要成因.研究区壳幔变形特征和机制不同,属于解耦变形.