利用近震高频接收函数研究四川理县西山村滑坡体结构

远震接收函数已广泛用于反演台站下方的结构,然而由于地球的非弹性衰减作用,远震数据较难获得高频接收函数,对浅地表结构约束不足.为了克服这一问题,我们使用近震数据的高频接收函数来研究浅表速度结构,并应用于四川理县西山村滑坡体上3个宽频带地震仪记录到的近震事件.本文发展了接收函数V_P-k（V_P/V_S）叠加方法,结合接收函数H-k叠加和波形反演方法获得了台站下方滑坡体的厚度、S波速度和平均V_P/V_S比,并与钻孔得到的滑坡体厚度进行对比.结果表明,滑坡体具有小尺度的横向不均匀性,台站下方滑坡体的平均V_P/V_S比在2.4~3.1之间变化并且在底层存在78~143 m·s^-1左右的S波低速层.本文观测到的高V_P/V_S比和底层低的S波速度结构,与电磁法获得的滑坡体底层低的电阻率和底部富水特征一致,表明滑坡体h1底界面的抗剪强度相对较弱,是潜在的滑坡危险区域.本文研究结果表明,利用近震接收函数能有效约束浅表的速度结构,进而能为滑坡灾害治理提供一定的地震学参考.

     

瞬变电磁采样函数优化法降噪技术

常用的双极性同步采样方法对瞬变电磁系统的观测噪声具有一定的抑制作用,但在实际工作中其噪声抑制能力仍显不足.本文在双极性同步采样方法基础上提出一种优化方法,以提升瞬变电磁系统的噪声抑制能力,即：通过对取样道加高斯窗函数的方式,赋予取样道内各采样点合适的累加权重,提升系统对取样道内高频噪声的抑制能力;使用高斯窗函数与梳状滤波器构造复合窗函数,通过对包含所有叠加周期的完整观测序列加复合窗函数的形式优化各周期的叠加权重,提升系统对实际工频噪声的抑制能力.对优化前后不同采样函数的谱特性及野外实测数据的处理效果进行比较,结果表明本文提出的方法较传统双极性同步采样方法具有更强的噪声抑制能力.

     

坡度随水平分辨率变化及其空间格局研究

以黄土丘陵沟壑区的县南沟流域为研究区,基于1∶1万地形图,利用ANUDEM软件生成5m到200m分辨率DEM,并利用Arc/Info中计算坡度的方法提取了各种分辨率的坡度。研究表明,随着DEM分辨率的降低,单个样点坡度值表现出不确定性,但同一坡度级所有栅格点的坡度均值呈现一定的规律性,低坡度段表现为先升高后降低,中坡度段呈现微弱变化,陡坡度段呈现对数降低趋势;沟沿线上坡度值呈比较剧烈的下降趋势、分水线和流水线上坡度缓慢下降。     

正则化滤波在反演地表光谱反照率中的应用

利用有限的卫星观测数据,通过反演地物二向性反射分布函数(BRDF)模型提取复杂地表的真实信息,是一个不适定的地球科学反演问题。为了克服离散不适定性带来的困难,在分析引起模型解不稳定的原因和数据误差传递机制的基础上,使用构建滤波函数的方法实现BRDF模型的正则化约束反演。实验结果表明,正则化滤波算法与MODIS AMBRALS算法精度相当,完全适用于具有充足观测和稀疏观测情况下的地物参数反演。     

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

地震数据是一个空间变化关系的数据载体,反映了地质体空间中的变化特征.在已知地震数据的前提下,如何刻画储层空间变化规律是精细储层研究的一个重要问题.笔者通过对地层空间转换函数进行研究,提出了一种基于地震道匹配求取地层空间转换函数的方法.该方法采用地震道多道匹配求取加权系数,使用多点的褶积加权的地层空间转换函数取代点对点的地层空间转换函数,从而准确反映储层物性变化,减小误差.基于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波低速体,说明可能缺乏解耦层,支持青藏高原西部地壳为整体缩短增厚模式.