利用广义反射-透射系数方法求解含低速层水平层状介质模型中面波频散曲线问题

在计算面波频散问题的广义反射-透射系数方法基础上研究了含有低速层的均匀平层固体地球模型的面波的频散曲线以及简正振型的计算问题. 当不存在低速层时，仅利用陈晓非1993年提出的单一的久期函数对求解这两个问题已经足够. 但是当模型中存在低速层时，该单一久期函数在数值求解较低阶模式很困难. 本文提出用久期函数族代替该单一久期函数、对不同的简正振型模式采用对应的久期函数来求解的方法，从而更加容易有效地求解较低阶的模式，完善了利用广义反射-透射系数方法求解面波频散问题的理论和算法.     

由单台远震P波偏振资料反演北京台站邻域的速度结构

地震台网的一个基本目的是研究地球的结构, 鉴于核爆炸事件的监测以及地震减灾的需要, 台网中台站邻域的地下结构就成了远震位置精确测定中非常重要的课题. 利用CDSN的北京(BJI)台站的宽频带、高采样率的三分量地震记录计算远震的偏振, 试图用这些偏振资料反演北京台站邻域的速度结构, 特别是速度间断面的细节. 结果发现该台站西部地壳有一明显的低速带, 这同该地区已有的研究结果十分一致. 这使我们相信, 它不仅在实际上验证了, 在理论上利用偏振资料可以反演速度结构、特别是速度梯度变化的边界形状的可能性; 也从实践上表明了, 用单台高质量宽频带记录的偏振数据反演速度结构的可行性. 在这个意义上, 走时数据和偏振数据可以联合用来研究速度结构, 而偏振数据更适合于确定速度异常的边界. 而且, 如果联合使用偏振方法和接收函数方法, 充分利用它们的互补性, 就有可能既可得到台站邻域横向速度结构又可得到台站下方细致的垂向结构.     

各向异性ATI介质剪切位错源地震矩张量

考虑震源区为各向异性ATI介质情形下,给出了剪切位错源地震矩张量解析表达式并模拟了震源区各向异性对双力偶分量(DC)、补偿线性偶极子分量(CLVD)和各向同性分量(ISO)的影响,结果表明,即使剪切位错源仍能产生非双力偶各向同性分量,其导致沿着ATI介质对称轴方向体积变化;当断层面位于ATI介质对称平面或者震源区为各向...     

利用经验格林函数识别地下核爆炸与天然地震

利用目标区域内的小当量地下核爆炸事件作为经验格林函数(EGF), 来估计发生在该目标区域内的可疑地震事件的相对震源时间函数(RSTF). 天然地震和地下核爆炸的震源物理过程的本质区别会通过震源时间函数, 表现出不同的特征. 研究结果表明, 地下核爆炸的远场RSTF表现为简单的、与方位角无关的单脉冲或双脉冲形式, 宽度约为1.0 s; 而基于地下核爆炸事件作为经验格林函数反演得到的天然地震的RSTF通常表现出很复杂的形式; 不仅表现为周期不固定的多脉冲形式, 而且随方位角、震中距的变化而变化. 对目标区域范围内发生的可疑地震事件, 利用已有的经验格林函数来反演其RSTF, 并通过对RSTF的分析来研究是否为地下核爆炸, 对地下核爆炸检测及识别研究不失为一种有效的参考手段.     

利用τ-p波场反演法确定东沙群岛的地壳速度模型

本文利用τ-p波场反演法直接对OBS原始数据进行转换,得到一维速度深度函数.该方法包含倾斜迭加和向下延拓两个线性转换:倾斜迭加将原始数据由时间-偏移距域(t-x)转换为时间截距-慢度域(τ-p);向下延拓将时间截距-慢度域(τ-p)波场转换为深度-慢度域(z-p)波场.然后根据五个OBS台站得到的一维速度深度函数,合成东沙群岛地壳的二维速度结构模型.为了检验该速度结构模型,我们用拾取走时和利用模型计算的理论走时进行拟合,得到了满意的结果,从而为下一步进行层析成像提供了可靠的初始速度结构模型.     

三维层状孔隙介质中弹性波的一种积分表达式I：理论

孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen（LAC）广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段.     

地震学百科知识(八)——计算地震学

正引言研究理论地震学计算问题所涉及的数值算法以及计算机实现是一门分支学科。它的诞生得益于计算机的广泛应用和大量的高质量地震波形观测资料的获得,是介于理论地震学、计算科学和观测地震学之间的新兴交叉学科。基于基本的物理定律和严格的数学工具,理论地震学给出了定量描述地震震源自     

利用曲线网格有限差分方法研究三维倾斜断层的破裂动力学

利用曲线网格有限差分方法,研究了三维倾斜断层的破裂传播过程.基于断层面生成贴体曲线网格,并通过坐标变换将含曲线网格的物理空间转换到含均匀直角网格的计算空间,实现了有限差分方法对复杂界面的处理.通过模拟地震断层的自相似破裂和自发破裂,并与已有发表的结果对比,发现拟合程度较高,验证了本方法的有效性和精确性.重点研究了不同倾角的倾斜断层破裂,最后展望了今后用本方法对非均匀介质中和任意起伏地表下的任意非平面断层破裂动力学的进一步研究.     

一种有效的适合F-J方法从背景噪声资料中提取多阶频散曲线的台阵划分方法

近年发展起来的从地震背景噪声提取多阶频散曲线的F-J方法是一种台阵分析方法,对某一特定的研究区域内台阵的划分方式直接影响成像的效果.当研究区域地下结构较简单,且台站间距比较均匀时,可以使用简单的平移式移动阵列法;然而若研究区域地下结构复杂多变,尤其是台站间距疏密不均时,就需要根据具体的地域特征灵活划分台阵.本文介绍了一种适用于地下结构复杂及台站间距不均匀区域的台阵划分方法,具有一定的普适性.我们将其应用于美国阿拉斯加及加拿大西北部,验证了相比于简单的沿经纬度移动阵列,我们的分区方法确实大大提高了结果的分辨率,为后续反演研究区域三维速度结构提供了坚实的保障.     

城市微动高阶面波在浅层勘探中的应用：以苏州河地区为例

频率-贝塞尔变换方法(Frequency-Bessel Transform method,简称F-J方法)是一种分析微动信号的新方法,由于该方法采用频率矢量波数变换处理水平层状各向同性弹性模型中时空平稳随机分布的微动信号,所以从理论上可以提取出清晰的瑞利波基阶和高阶模态频散曲线,但是目前还没有相关的野外实验对此进行研究和应用.本文首先采用该方法对上海市苏州河地区采集的城市微动信号进行处理获得了频率-相速度谱,然后提取了多模态瑞利波频散曲线,最后通过粒子群算法对频散曲线进行联合反演,得到了浅地表0~70 m深度范围的S波速度结构,并且利用钻孔数据对反演的速度结构进行了验证.另外,本文还通过对比F-J方法和传统的SPAC(SPatial AutoCorrelation method)方法分别提取的频散曲线,展示了F-J方法在处理城市微动信号方面的优势.本文研究结果表明:(1)F-J方法可以从少量台站(21个台站)短时记录(1小时)的微动信号垂直分量中提取出清晰的基阶和高阶模态瑞利波频散曲线;(2)F-J方法提取的高阶模态频散曲线比传统SPAC方法提取的更加清晰,高频部分(>13 Hz)优势更为明显;(3)联合基阶和高阶模态频散曲线反演的浅地表速度结构更加精确,可以分辨出第四系沉积层中物性相差较小的速度界面和低速异常,在城市浅地表精细结构成像方面具有较好的应用前景.