空间波数混合域磁异常场积分解三维数值模拟

本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点,沿水平方向进行二维傅里叶变换,把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域,优势之一在于便于浅层单元剖分可适当加密,随着深度增加,单元剖分适当稀疏,可以准确模拟任意复杂地形和磁性体的磁异常,兼顾了计算精度与计算效率;优势之二在于一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征磁化强度,可得出单元积分的解析表达式,计算精度高、效率高.该方法充分利用一维形函数积分的高效和高精度、快速傅里叶变换的高效性及算法高度并行性,实现了磁异常场高效、高精度的数值模拟.设计棱柱体模型,将模型解析解与空间波数混合域法的数值解对比,结果表明该方法计算精度高、效率高.设计了组合棱柱体复杂模型,对比分析了标准FFT扩边法与Gauss-FFT法的计算精度与计算效率,总结了标准FFT的扩边系数选取策略.针对任意复杂地形条件下的磁异常模拟问题,本文提出一种适用于起伏地形条件下的磁异常场快速计算方法,并对其有效性进行了验证.     

基于三维应变格林函数反演中小震震源机制

随着复杂速度结构反演的发展和高性能计算能力的提升,基于高精度3D介质模型计算格林函数反演震源机制更具可行性。中小地震因具有更好的覆盖和近似点源效应,在区域结构成像中有着广泛的应用。基于波形类的反演方法如波动方程层析成像\,全波形反演都需要震源机制解,而传统的震源机制反演方法不能很好地应用于中小地震。采用有限差分法构建应变格林张量(Strain Green Tensor,SGT)数据库,将合成波形和实际波形按震相截窗并滤波到不同的频带范围,先通过最小化互相关走时差来进行震源重定位,再通过最小化波形残差反演震源机制。通过合成数据测试验证方法的正确性,随后将该方法应用于青藏高原东部边缘龙门山地区,反演一系列M_W3.4～5.7的中小地震震源机制。由于应变格林张量数据库可预先构建,该方法可以应用于(近)实时震源机制解反演。     

元谋地电场观测数据小波阈值去噪

为了解决元谋地电场观测数据的噪声干扰问题,采用小波阈值去噪的方法,根据元谋地电场观测数据特征,选取适当的阈值及不同的小波函数、分解层数进行研究。综合各小波函数去噪效果,通过对相对误差、均方根误差、信噪比等指标进行评价,确定适合元谋地电场观测数据去噪的阈值为软阈值函数、分层阈值;最优小波函数分别为Sym4、Sym5、db6和db8;分解层数分别为5、7、9。最后应用研究结果处理实际地电场观测数据,取得了较好的去噪效果。     

上海区域量规函数研究

量规函数是用来校准补偿地震波随距离的衰减的,本文对上海区域台网的量规函数进行重新计算,同时和全国通用的量规函数加以对比,观察是否需要对上海区域量规函数进行修正。并使用统计平均震级和M_L震级定义两种方法来计算量规函数。      

碰撞非相干散射理论谱宏观法和微观法的比较研究

非相干散射理论谱是理解非相干散射雷达探测的理论基础,其求解过程分为微观法和宏观法.在考虑碰撞不考虑磁场的条件下,本文使用这两种方法对理论谱分别进行求解,给出了关键求解过程,并对等离子体的归一化导纳和极化率以及非相干散射的微分散射截面和谱密度函数进行分析,得到了它们各自之间的数学关系,概括性分析了两种方法在理论选择、最终结论及多种因素下求解的异同点.     

鄂尔多斯西缘北段的地壳结构和块体间变形关系

本文对喜马拉雅计划二期部分台站的远震波形数据进行接收函数提取,利用接收函数共转换点叠加方法获得阿拉善地块、鄂尔多斯地块以及银川—河套盆地下方0~80 km深度的速度间断面结构.结果表明:鄂尔多斯地块成层性好,地壳厚度为38~42 km,康拉德界面为18~22 km,阿拉善地区的Moho面深度为38~45 km.河套盆地地壳厚度约52 km,银川断陷盆地和贺兰山下方的Moho面最深为~55 km.鄂尔多斯西缘构造边界下方Moho面变化明显,且黄河断裂为深大断裂直接切割莫霍界面.根据本文的间断面成像结果我们进一步确定阿拉善地块与鄂尔多斯地块分属不同的大地构造单元.与此同时,我们推测贺兰山以西70~80 km范围内和鄂尔多斯地块西缘北段存在地壳增厚变形的可能.     

2013年芦山地震震源区地壳介质地震波速变化的特征分析

为探究芦山M7.0级地震后5年多来,震源区龙门山断裂带西南段介质波速的变化规律,本文基于2012年4月至2018年4月共6年的连续波形数据,运用移动窗互谱与频域偏振等分析方法,结合背景噪声源的特性,对不同深度范围内的相对波速变化以及震后的恢复过程与机制进行了研究.获得的主要认识包括：（1）年尺度而言,震源区周期为1~20 s的背景噪声场相对稳定,但成分复杂、2~10 s频带内至少存在2个能量相对稳定的噪声源;不同周期噪声的能量,在月变与季节性上的变化特征差异明显.（2）获得了长时间尺度、不同频带内介质相对波速的背景变化水平,1~2 s、2~4 s的波动幅度（约为±0.04%）与季节性变化规律强于4~10 s、10~20 s的,结合与降雨量相关的地下水位模型能很好地解释其变化规律.（3）震源区的同震波速降低现象清晰,降幅约为0.08%~0.1%;空间上,波速下降最为显著的区域主要集中在龙门山断裂带两侧约70 km范围内,其中四川盆地一侧平均约为0.1%,略高于青藏高原（0.08%）一侧;在断裂带内的降速不显著.对不同子频带进行测量的结果显示,震后除10~20 s外,其余3个子频带的相对波速在震后较短时间内（约20天左右）均出现较大幅度的波速降低现象,其中4~10 s的平均降速最大（约为0.08%）,分析认为主震及大量余震的松弛效应是引起介质波速下降的主要原因.（4）震后大约1年左右,波速变化基本恢复到震前水平,且至2018年4月前未观察到大幅的波速变化现象,总体上各频带内的结果均沿零线小幅波动.     

倾角函数展开及其在分析法轨道预报中的应用

田谐项摄动是分析法轨道预报中的重要部分,其中包含大量倾角函数及其偏导数的计算.由于具有精度更高、速度更快的优点,倾角函数一般通过递推方法计算.以文献中提出的改进Gooding方法为基础,将其给出的程序稍加改进,在计算2–50阶倾角函数时缩短了约24%的计算时间.考虑到分析法预报过程中轨道平倾角变化很小,以泰勒展开式计算倾角函数,可极大提高计算速度,较大程度地减小分析法预报耗时,且引力场阶次越高,减小幅度越大,取50阶时预报耗时缩短了48%.另一方面,以2阶展开式计算倾角函数时,与改进Gooding法相比,分析法预报星历偏差很小.对于500 km高度的低轨卫星,分别以改进Gooding法和2阶泰勒展开式计算倾角函数,预报3天,当地球引力场阶次不高于50时,二者预报星历偏差RMS (Root Mean Square)低于1 mm,且随着轨道高度的增加,预报星历偏差RMS逐渐减小.     

湖泊沉积物中元素相态的连续提取分析—以岱海为例

如何利用湖泊/海洋沉积物地球化学代用指标有效地指示环境很大程度取决于元素的来源及其在沉积物中赋存形态的确定.对位于半干旱区的内蒙古岱海14个沉积物样品的连续提取实验及18个元素在水可溶态、可交换态、碳酸盐结合态、铁锰氧化物结合态、有机物结合态及残留态等六种赋存相态中含量变化进行了测试和研究,结合沉积物的矿物组成,结果表明在岱海流域内只经历了较弱的化学风化过程,除了Ca和Sr外,其他元素主要赋存于残留态中,只发生了很小程度的活化溶解.受元素亲合性和湖水物理化学状态的制约,出溶的Sr与Ca的赋存相态基本一致,而出溶的Pb、Be、Co、Cu、K、Mn、Fe等元素则与自生碳酸盐、有机质表现出一定的亲合性.元素在不同赋存相态中分布为地球化学环境代用指标的提取和解释提供了机理性的实验依据.     

中尺度流域次降雨洪水过程模拟——以太湖上游西苕溪流域为例

古近系沙河街组是阳信洼陷最重要的生烃层段.其湖盆充填具阶段性和沉积演化的旋回性,形成了由沙四段、沙三段与沙一段组成的复合生烃系统和不同的地球化学特征.利用Rock-Eval生油岩评价仪、色谱-质谱仪等实验分析技术对不同层系样品进行了地球化学特征分析研究,其中生烃潜力指标包括有机碳含量(TOC)、残余生烃潜量(S1 S2)、氢指数IH、氯仿沥青"A";成熟度指标包括热解峰温Tmax、奇偶优势参数OEP、饱和烃轻重比(C21 C22)/(C28 C29)及∑C21-/∑C22 、镜质体反射率Ro、莫烷/藿烷及甾烷生物标志物参数C29ββ(ββ αα)、C29ααα20S(20S 20R);古环境指标包括异构烷烃参数Pr/Ph、Pr/n-C17、Ph/n-C18及伽马蜡烷等.结果表明,沙四段为弱还原-还原性的半深湖沉积,沉积了一套中等厚度、分布局限的烃源岩(TOC平均为1.5%),其较深位置的烃源岩基本进入成熟阶段,多形成成熟油;沙三段属弱还原-还原性的半深湖-深湖环境,其烃源岩中有机质丰度较高(TOC平均为3.5%),大部分烃源岩处于未成熟-低成熟状态,主要形成未熟油;沙一段为湖水咸化、还原性的半深湖相环境,其烃源岩中有机质丰度高(TOC平均为5%左右),但处于未成熟阶段,主要生成生物气.