地下爆炸各向同性源和补偿线性矢量偶极源共同作用对Rg波低谷点的影响

本文利用地震面波激发理论计算了两种震源不同相对强度、相对埋深及源时间函数情况下的Rg波频谱,研究地下爆炸中各向同性(ISO)源和补偿线性矢量偶极(CLVD)源共同激发的Rg波信号低谷点特性.结果表明,由于两种震源相对强度、相对埋深及源频谱的影响,混合信号低谷点频率与单一CLVD源激发的信号有很大差异,因此实际地下爆炸中...     

地下爆炸Rg波低谷点激发机理

目前普遍认为Rg_S散射是地下核爆炸短周期Lg波的主要能量来源,其主要依据为由线性偶极补偿源（CLVD）激发的Rg波谱的低谷点随深度变化的特征与区域震相Lg波谱的低谷点随深度变化的特征具有一致性. 本文利用简正振型理论,分析了Rg波谱中低谷点的形成机理,给出了各种速度模型下CLVD源深度与低谷点频率之间的关系,并分析了利用公式hCLVD=V/(16fNull)估计震源深度的局限性. 本文的研究结果进一步支持了CLVD源是激发Lg波的主要因素的观点,并得到重要结论：CLVD源形式下的本征位移函数及其导数的叠加所形成的极小值,即源的空间分布特征,是形成Rg波低谷点的原因,而不能仅用水平向基频本征位移函数过零点来解释.     

层裂对区域震相Lg波的影响

地下核爆炸Lg波的激发机制, 目前尚未十分清楚. 普遍接受的观点是: 伴随层裂源的补偿线性矢量偶极源激发的Rg波的散射形成的S波, 是Lg波的最主要贡献因素. 本文利用频谱比对方法, 结合合成理论地震图方法, 基于东哈萨克斯坦地区速度模型, 分析了东哈萨克斯坦地壳速度模型下的层裂对区域震相Lg波的影响, 并进一步分析了该地区的实际观测资料. 结果表明, 层裂时间函数对Lg波具有扇贝形调制作用; Lg波振幅谱的低谷点源自于补偿线性矢量偶极（CLVD）源激发的Rg波. 这一结果也表明, 东哈萨克斯坦速度模型下的Lg波受层裂的调制作用, 与内华达试验区的Lg波是一样的, 进一步支持了Patton和Taylor的观点, 对于进一步理解及更好地利用Lg波具有重要理论指导意义.      

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.     

用于弹性波模拟的交错网格预处理和多次校正的Chebyshev谱元法

本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。     

基于混合阶矢量基函数的海洋可控源电磁三维谱元法数值模拟

在前人工作的基础上,本文推导了电导率任意各向异性介质的海洋可控源电磁三维谱元法正演方程.采用一次场/二次场分离算法结合混合阶矢量基函数,可以有效避免源点的奇异性的影响,从而提高数值解的精度.采用任意六面体单元离散研究区域,有利于模拟复杂地形和地电结构.利用不完全LU分解的Induced Dimension Reduction（IDR（s））迭代算法求解线性方程组,有效地提高了求解的效率.设计典型的地电模型进行正演计算,并将计算结果与有限元解进行对比,对比结果表明本文提出的基于混合阶矢量基函数的海洋可控源电磁三维谱元数值模拟算法是正确的、有效的.本文算法具有良好的通用性,可推广用于电导率呈任意各向异性的陆地电磁、井中电磁等数值模拟研究.

     

垂直和倾斜断层围陷波的谱元法数值模拟及其波场特征对比

低速断层破碎带内传播的断层围陷波是断层内部信息的重要载体,研究其波场特征是反,演断层内部结构的关键。前人对围陷波的模拟主要针对垂直断层,但实际断层大部分是倾斜的。本文利用谱元法数值模拟地震波在垂直和倾斜断层区域的传播,并分别单独改变断层破碎带或断层两盘速度差异等参数,对比分析地震动在时域和频域发生的变化,研究这些参数对波场特征的影响。对比研究表明,当存在一定宽度的低速破碎带时,倾斜断层对近断层区域波场的影响和垂直断层相比在幅值分布上有比较显著的差异,尤其是当断层两盘同时还存在速度差异时;在倾斜和垂直断层破碎带上方,均能产生断层围陷波,其幅值特征差异不大,但在跨断层平面上的质点极化方式可能会存在差异:垂直断层围陷波质点极化方向表现为单一,而倾斜断层可能会出现多个方向。这些模拟结果有助于反演真实断层破碎带结构,也能够加深对近断层的强地面运动特征的认识。     

可控源海洋电磁勘探中空气波影响的理论分析和数值模拟

可控源海洋电磁勘探(MCSEM)中空气波对海底电磁响应的影响已为业界所熟知,且已提出了多种压制方法.然而,由于空气波与海水层中其它信号相互作用的复杂性,至今其作用机理及对有效信号的影响方式尚未完全清楚,这也阻碍了浅水域MCSEM油气勘探的应用.本文在前人研究的基础上, 基于一维层状介质的电磁位和电磁场分析,采用电磁场的模式分解理论导出了空气、海水层和海底地层之间的相互作用关系.分析了MCSEM电磁响应的影响因素和空气波的作用机理.以无限水深假想模型、三层介质模型和四层介质模型为基础,导出了水平电偶极子(HED)的横电(TE)和横磁(TM)模式的电磁场关系,分析显示了前人提出的上下行波分离进行空气波压制所存在的缺陷.通过理论和数值模拟知道, 利用横电(TE)和横磁(TM)模式受空气相互作用影响的差异和电磁场水平分量与垂直分量受空气层相互作用影响程度的不同可减弱空气波对有效信号的影响,这将有助于实现浅水域海洋电磁勘探数据的有效利用.     

可控源海洋电磁勘探中空气波影响的理论分析和数值模拟

可控源海洋电磁勘探(MCSEM)中空气波对海底电磁响应的影响已为业界所熟知,且已提出了多种压制方法.然而,由于空气波与海水层中其它信号相互作用的复杂性,至今其作用机理及对有效信号的影响方式尚未完全清楚,这也阻碍了浅水域MCSEM油气勘探的应用.本文在前人研究的基础上, 基于一维层状介质的电磁位和电磁场分析,采用电磁场的模式分解理论导出了空气、海水层和海底地层之间的相互作用关系.分析了MCSEM电磁响应的影响因素和空气波的作用机理.以无限水深假想模型、三层介质模型和四层介质模型为基础,导出了水平电偶极子(HED)的横电(TE)和横磁(TM)模式的电磁场关系,分析显示了前人提出的上下行波分离进行空气波压制所存在的缺陷.通过理论和数值模拟知道, 利用横电(TE)和横磁(TM)模式受空气相互作用影响的差异和电磁场水平分量与垂直分量受空气层相互作用影响程度的不同可减弱空气波对有效信号的影响,这将有助于实现浅水域海洋电磁勘探数据的有效利用.     

热源激发重力波特征以及波流作用的数值模拟研究

本文在二维等温可压大气中引入了一个随时间和空间变化的热源扰动,分别以静止风和中纬1月份月平均向东的纬向风急流为背景,对不同背景下热源激发的重力波的传播详细过程及其特性进行了数值模拟研究.热源激发出来的重力波在初始阶段有很宽的频谱范围,随后由于重力波的传播效应,水平波长和垂直波长分布范围随时间都有所减小.顺风传播的重力波的小尺度和低频部分会容易被急流吸收,从而加强了对流层急流,而逆风传播的重力波更容易上传,会导致中间层区域向西的背景风增强.这体现了低层大气急流对中间层大气风场结构的影响.热源的尺度直接决定激发波的尺度;激发波的垂直尺度和时间尺度对热源的变化比其水平尺度更敏感.     

利用谱元法研究震源参数对于近断层地震动上盘效应的影响

国内外已有诸多学者基于强震记录或数值模拟方法开展了近断层地震动的上盘效应分析,并获得了初步的认识,但上盘效应随震源参数的变化规律仍有待深入探讨。本文使用运动学有限断层震源模型,运用高精度谱元法模拟逆断层破裂模式下近断层区域地震动,对比上、下盘地震动差异进而分析近断层地震动上盘效应。通过改变单一震源参数,研究断层上界埋深、震级大小以及断层倾角对地震动上盘效应的影响规律。研究结果表明：（1）使用运动学有限断层震源的谱元法数值模拟可以较好地反映近断层上盘效应;（2）上盘效应随测点到断层上界在地表投影的水平距离Rx的增加先增大后减小,最终趋于稳定。（3）上盘效应随断层上界埋深的增加先增大后减小,且随着断层上界埋深的增加,上盘效应峰值区会向远离断层破裂迹线的方向扩展。（4）上盘效应受矩震级影响较小,但总体随矩震级的增大而增大;上盘效应峰值区到断层破裂迹线的距离不随矩震级的变化而变化。（5）随着断层倾角的增加,上盘效应先增大后减小,当断层倾角为45°时上盘效应最为明显;随着断层倾角增加,上盘效应峰值区会向远离断层破裂迹线的方向扩展;当断层倾角为90°时不存在上盘效应。     

Rectangular footing on soil with depth-degrading stiffness: Vertical and rocking impedances under conditional existence of surface waves

The response of rectangular rigid footings resting on an elastic soil of shear modulus decreasing monotonically with depth is studied. Such profiles are typically encountered after ground improvement. The propagation characteristics of SV/P surface waves are investigated, showing the appearance of cut-off frequencies above which surface waves do not exist. The semi-analytical method of the subdivision of the footing/soil contact area is then used for solving the boundary value problem, whereby the influence functions for the sub-regions are determined by integration of the corresponding surface-to-surface Green׳s functions. Impedance functions are presented over a wide range of frequencies for typical values of the non-homogeneity parameters, the Poisson׳s ratio and the foundation geometry. The salient features that are associated with the non-homogeneity and the appearance of cut-off frequencies are elucidated.     

基于南水北调西线工程岩性特征的CSAMT法有限元三维数值模拟研究

有限元法是地球物理数值模拟中常用的方法,本文采用三维可控源音频大地电磁法(CSAMT)有限元数值模拟的程序,根据南水北调西线工程岩性的地球物理特征,设置了不同的三维模型,并对其进行了有限元数值计算分析,从三维空间中模拟场的规律,探索了不同地质异常体的特征,为提高可控源音频大地电磁法在南水北调西线工程深埋藏隧洞探测地质异常体的精度奠定了基础.