谱元法模拟地震波传播地震点源处理方法

1

作者 李孟洋 刘少林 +1 位作者 杨顶辉 李小凡 《地震研究》 北大核心 2025年第3期401-412,共12页

在利用谱元法求解地震波运动方程时需要对方程进行离散,然而由于震源项常包含δ函数,该函数是奇异函数,通过数值近似具有一定难度,因此需要采用特殊方式处理震源项。引入δ函数的2类离散形式,以谱元法求解地震波运动方程为例,详细给出... 在利用谱元法求解地震波运动方程时需要对方程进行离散,然而由于震源项常包含δ函数,该函数是奇异函数,通过数值近似具有一定难度,因此需要采用特殊方式处理震源项。引入δ函数的2类离散形式,以谱元法求解地震波运动方程为例,详细给出了震源项离散形式。为了验证这种震源处理方法的有效性,开展了多组数值模拟实验。首先在规则网格下,模拟了δ函数离散形式处理单极震源和矩张量震源的情况;其次在非规则网格下,模拟了δ函数离散形式处理地震点源情况。数值实验表明,在谱元法模拟地震波传播时,利用δ函数离散形式近似地震点源项,可得到高精度数值模拟结果。 展开更多

关键词 谱元法 地震波运动方程 数值模拟 Δ函数

在线阅读 下载PDF 职称材料

基于预条件广义逐次超松弛迭代法的数值格林函数计算方法 被引量：1

2

作者 徐杨杨 商耀达 孙建国 《吉林大学学报（地球科学版）》 CAS CSCD 北大核心 2024年第5期1696-1710,共15页

为了改善Born散射级数解决地震强散射问题时的收敛性,将带有虚部分量的复波数格林函数引入到求解格林函数Lippmann–Schwinger(L-S)积分方程数值解的广义逐次超松弛迭代法中,弱化格林函数的奇异性。引入预条件算子降低系数矩阵的条件数... 为了改善Born散射级数解决地震强散射问题时的收敛性,将带有虚部分量的复波数格林函数引入到求解格林函数Lippmann–Schwinger(L-S)积分方程数值解的广义逐次超松弛迭代法中,弱化格林函数的奇异性。引入预条件算子降低系数矩阵的条件数,加速迭代级数的收敛速度,给出了复波数L-S方程的预条件广义逐次超松弛(preconditioned generalized successive over-relaxation,Pre-GSOR)迭代格式。通过数值分析和收敛性分析重新选取合适的衰减因子和预条件算子,得到了满足地震强散射条件的收敛Born级数,并将其用于地震强散射问题中数值格林函数的计算。数值结果表明:复波数L-S方程Pre-GSOR迭代法可以得到与实波数L-S方程直接法相匹配的数值模拟结果;复波数L-S方程Pre-GSOR迭代法系数矩阵条件数在高频时仅为原系数矩阵条件数的10%,相同迭代次数下归一化收敛残差可降低3个数量级以上,且对高频适应性强,可有效改善实波数L-S方程广义超松弛迭代法在强散射介质中的收敛停滞问题。 展开更多

关键词 地震散射波场 格林函数 广义超松弛迭代 预条件算子

在线阅读 下载PDF 职称材料

基于CEEMDAN-INHT的地下洞室爆破振动时频分析应用研究 被引量：1

3

作者 孙苗 吴立 杨钧凯 《爆破》 CSCD 北大核心 2024年第1期14-20,共7页

爆破地震波信号采集会因监测环境、测试系统等因素导致实测信号中混有噪声,噪声的存在将导致信号希尔伯特-黄变换(Hilbert-Huang Transform,HHT)时频分析结果失真。原因有二:其一是经验模态分解(Ensemble Empirical Mode,EMD)处理含噪... 爆破地震波信号采集会因监测环境、测试系统等因素导致实测信号中混有噪声,噪声的存在将导致信号希尔伯特-黄变换(Hilbert-Huang Transform,HHT)时频分析结果失真。原因有二:其一是经验模态分解(Ensemble Empirical Mode,EMD)处理含噪爆破地震波信号会得到具有模态混淆现象的固有模态函数(Intrinsic mode function,IMF)分量;其二是Hilbert变换受Bedrosian定理的约束在处理模态混淆分量时会产生负值瞬时频率,从而造成巨大的分析误差。为获得真实的爆破振动属性需对HHT进行改进,在EMD中添加自适应噪声信号得到自适应补充集合经验模态分解(Complete Ensemble Empirical Mode Decomposition with Adaptive Noise,CEEMDAN)算法;再对CEEMDAN得到的IMF进行归一化Hilbert变换,得到改进归一化Hilbert变换(Improved Normalized Hilbert Transform,INHT)。通过上述两步可建立CEEMDAN-INHT时频分析算法,为验证该算法可有效提高含噪爆破地震波振动信号时频分析精度,进行HHT和CEEMDAN-INHT含噪仿真振动信号时频分析对比研究。最后将CEEMDAN-INHT用于某地下洞室爆破地震波信号时频分析中,发现该算法能有效克服EMD固有的模态混淆现象,同时得到反映真实爆破振动属性的时-频-能特征参数,对从频率-能量的角度进行洞室爆破开挖共振分析,实现爆破地震波危害控制具有一定的现实意义。 展开更多

关键词 爆破地震波信号 经验模态分解 HILBERT变换 固有模态函数

在线阅读 下载PDF 职称材料

反平面线源荷载作用下浅埋圆形非完全粘结隧道动力响应研究

4

作者 李志文 李潜 +2 位作者 徐斌 李晓锋 李海波 《爆炸与冲击》 EI CAS CSCD 北大核心 2024年第8期60-71,共12页

为加深理解波源距离和非完全粘结对地震波散射的影响规律,结合位移不连续模型、波函数展开法、Graf公式和镜像方法推导了反平面线源荷载下浅埋圆形非完全粘结隧道动力响应的级数解,并通过衬砌内外边界条件残余量与级数解截断项数的关系... 为加深理解波源距离和非完全粘结对地震波散射的影响规律,结合位移不连续模型、波函数展开法、Graf公式和镜像方法推导了反平面线源荷载下浅埋圆形非完全粘结隧道动力响应的级数解,并通过衬砌内外边界条件残余量与级数解截断项数的关系校验了该解的精度。通过对该级数解进行参数分析,系统地探讨了衬砌与围岩的接触刚度、衬砌模量、衬砌厚度、隧道埋深和波源距离等因素对衬砌内表面位移和周向剪应力的影响。结果表明:衬砌与围岩的接触刚度对隧道的动力响应具有显著的影响,尤其在某些较小接触刚度情况下隧道动力响应幅值可能非常大;增大衬砌模量会减小位移,但同时会导致周向剪应力增加;增大衬砌厚度能同时减小位移和周向剪应力;增大隧道埋深会使最大位移和周向剪应力向隧道拱顶附近移动;增大线源与隧道的水平距离会使隧道背波侧相对幅值增大。 展开更多

关键词 浅埋隧道 地震波散射 反平面线源 波函数展开法

在线阅读 下载PDF 职称材料

爆破地震波作用下隧道围岩动应力集中系数分析 被引量：21

5

作者 李兴华 龙源 +2 位作者 纪冲 周翔 路亮 《岩土工程学报》 EI CAS CSCD 北大核心 2013年第3期578-582,共5页

将爆破地震波进行合理近似简化,采用波函数展开法,推导了无限岩石介质中爆破地震波作用下隧道围岩动应力集中系数的表达式,结合具体算例,从理论上定量分析了围岩不同位置上爆破地震波主频和围岩物理力学指标对动应力集中系数的影响。分... 将爆破地震波进行合理近似简化,采用波函数展开法,推导了无限岩石介质中爆破地震波作用下隧道围岩动应力集中系数的表达式,结合具体算例,从理论上定量分析了围岩不同位置上爆破地震波主频和围岩物理力学指标对动应力集中系数的影响。分析结果表明:隧道围岩迎爆侧的动应力集中系数整体上大于背爆侧,动应力集中系数峰值偏向于迎爆侧;在不同主频爆破地震波作用下,动应力集中系数的影响存在最大值;隧道围岩的弹性模量对动应力集中系数的影响较大,而泊松比对动应力集中系数的影响较小。 展开更多

关键词 爆破地震波 波函数展开法 隧道围岩 动应力集中系数

在线阅读 下载PDF 职称材料

爆破地震波作用下既有圆形隧道衬砌动应力集中系数分析 被引量：22

6

作者 李兴华 龙源 +3 位作者 纪冲 周翔 何洋扬 路亮 《岩土力学》 EI CAS CSCD 北大核心 2013年第8期2218-2224,共7页

将爆破地震波进行合理近似简化,采用波函数展开法,推导出无限岩石介质中爆破地震波作用下隧道围岩和衬砌动应力集中系数的表达式。结合具体算例,分析隧道不同位置上爆破地震波主频、隧道几何参数和隧道围岩和衬砌物理力学指标对动应力... 将爆破地震波进行合理近似简化,采用波函数展开法,推导出无限岩石介质中爆破地震波作用下隧道围岩和衬砌动应力集中系数的表达式。结合具体算例,分析隧道不同位置上爆破地震波主频、隧道几何参数和隧道围岩和衬砌物理力学指标对动应力集中系数的影响。分析结果表明,爆破地震波主振频率对隧道围岩和衬砌中的动应力集中系数(DSCF)影响较大,隧道不同位置和不同方向应力的动应力集中系数随主振频率的变化趋势不同,低主振频率引起的动应力集中系数较大,不利于隧道的抗爆破振动;隧道几何参数对隧道围岩和衬砌内侧切向应力动应力集中系数影响趋势不同;隧道衬砌弹性模量对隧道结构动应力集中系数影响较大,泊松比次之;泊松比较大的隧道围岩具有一定的抗爆破振动作用,但效果有限。 展开更多

关键词 爆破地震波 波函数展开法 隧道围岩 隧道衬砌 动应力集中系数

在线阅读 下载PDF 职称材料

多孔岩石中应力波传播特性及震源函数 被引量：5

7

作者 周钟 王肖钧 +1 位作者 肖卫国 劳俊 《中国科学技术大学学报》 CAS CSCD 北大核心 2008年第11期1310-1316,共7页

建立了两种着重描述多孔岩石的塑性压缩(孔隙坍塌)机理的弹塑性本构模型,分别为将孔隙压实过程计入P-V关系的变模量理想塑性包络面模型和将孔隙演化行为计入帽盖的双屈服面帽盖模型.在此基础上,对地下球形药包起爆引起的弹塑性波在多孔... 建立了两种着重描述多孔岩石的塑性压缩(孔隙坍塌)机理的弹塑性本构模型,分别为将孔隙压实过程计入P-V关系的变模量理想塑性包络面模型和将孔隙演化行为计入帽盖的双屈服面帽盖模型.在此基础上,对地下球形药包起爆引起的弹塑性波在多孔凝灰岩中的传播过程进行计算,揭示了多孔岩石中应力波传播的一些基本特性.并采用改进的帽盖模型,研究了多孔凝灰岩中地下爆炸震源函数的特性,着重分析了孔隙度对震源函数的稳态值、拐角频率以及过冲特性的影响. 展开更多

关键词 多孔岩石 本构模型 孔隙演化 应力波 震源函数

在线阅读 下载PDF 职称材料

P波入射下浅埋圆形复合式衬砌隧道减震力学特性研究 被引量：10

8

作者 王帅帅 高波 申玉生 《铁道学报》 EI CAS CSCD 北大核心 2016年第6期108-116,共9页

本文利用波函数展开法,给出平面P波入射下浅埋圆形双层衬砌隧道动应力集中系数级数解,并分析P波入射角和入射频率对不同衬砌刚度组合隧道动应力集中系数的影响。研究表明:P波入射角度和入射频率对衬砌洞室动应力集中系数有明显影响,随着... 本文利用波函数展开法,给出平面P波入射下浅埋圆形双层衬砌隧道动应力集中系数级数解,并分析P波入射角和入射频率对不同衬砌刚度组合隧道动应力集中系数的影响。研究表明:P波入射角度和入射频率对衬砌洞室动应力集中系数有明显影响,随着P波入射角和入射频率增大,衬砌的动应力集中系数分布变复杂;在低频入射情况下,当组合式衬砌洞室中衬砌1为刚度较低的柔性衬砌和减震层时,刚性衬砌2的动应力集中系数明显降低,衬砌1刚度越低,衬砌2动应力集中系数越小;与高频入射波相比,柔性衬砌和减震层对低频入射波下衬砌2的动应力有更好的减震效果。本文结合某穿越Ⅸ度地震区浅埋管道隧道工程,研究减震层厚度和弹性模量对"围岩-减震层-衬砌"洞室动应力集中系数的影响,建议减震层弹性模量小于围岩弹性模量的1/20,减震层最优厚度宜取洞室内净空直径的1/25。 展开更多

关键词 P波 圆形洞室 组合式衬砌 柔性衬砌 减震层 动应力集中系数

在线阅读 下载PDF 职称材料

基于地震波衰减的特征属性重构方法 被引量：7

9

作者 魏文 王兴谋 +2 位作者 李红梅 王树刚 王红 《石油物探》 EI CSCD 北大核心 2012年第3期219-224,209,共6页

以储层流体识别为目的,提出了基于地震波衰减的特征属性重构方法。首先从地震波衰减理论出发,通过改变地层品质因子表达式,推导了地震特征属性函数,利用振幅类、频率类及衰减属性重构新的地震属性;然后设计粘弹性介质模型进行试算,结果... 以储层流体识别为目的,提出了基于地震波衰减的特征属性重构方法。首先从地震波衰减理论出发,通过改变地层品质因子表达式,推导了地震特征属性函数,利用振幅类、频率类及衰减属性重构新的地震属性;然后设计粘弹性介质模型进行试算,结果表明,利用特征属性函数重构的地震属性能较好地突出储层含流体变化特征,从而验证了基于地震波衰减的特征属性重构方法识别储层含油气性的有效性。将该方法应用于老河口油田北部老168井区河道砂岩储层油气识别,基于叠后地震资料精细刻画了储层流体的展布特征,取得了较好的应用效果。 展开更多

关键词 地震波衰减 特征属性函数 老168井区 河道砂岩 油气识别

在线阅读 下载PDF 职称材料

流体饱和层状半空间中沉积谷地对地震波的散射—IBIEM求解 被引量：8

10

作者 刘中宪 梁建文 赵瑞斌 《岩土工程学报》 EI CAS CSCD 北大核心 2013年第3期512-522,共11页

针对流体饱和层状半空间中任意形状沉积谷地对地震波的散射问题,基于单层位势理论,建立一种新的间接边界积分方程法(IBIEM):首先推导饱和层状半空间中两类膨胀波源和剪切波源动力格林函数,进而由分布在沉积谷地和层状半空间交界面附近... 针对流体饱和层状半空间中任意形状沉积谷地对地震波的散射问题,基于单层位势理论,建立一种新的间接边界积分方程法(IBIEM):首先推导饱和层状半空间中两类膨胀波源和剪切波源动力格林函数,进而由分布在沉积谷地和层状半空间交界面附近两虚拟波源面上的PI、PII和SV波源分别构造沉积内外的散射波场。由交界面连续性条件建立方程并求解确定虚拟波源密度,总波场即可由自由波场和散射波场叠加而得。然后通过边界条件验算、退化解答与现有结果的比较,验证了方法的计算精度。通过一组典型算例,初步探讨了平面SV波在饱和层状半空间中沉积谷地周围散射的基本规律。研究表明:饱和层状沉积对地震波的散射和相应干土情况存在显著差异,边界渗透条件对沉积内部地震响应具有重要影响;饱和层状场地中波的散射和均质半空间情况有着本质不同,总体反应体现了饱和土层的共振放大效应和沉积谷地对波的散射效应。实际场地中地震波传播的精确模拟需考虑岩土介质的多相性和层理特征。 展开更多

关键词 流体饱和层状半空间 沉积谷地 地震波散射 格林函数 间接边界积分方程法(IBIEM)

在线阅读 下载PDF 职称材料

基于Love波相速度反演南北地震带地壳上地幔结构 被引量：12

11

作者 徐小明 丁志峰 +1 位作者 叶庆东 吕苗苗 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2015年第11期3928-3940,共13页

收集了南北地震带区域地震台网中292个地震台站2008年1月至2011年3月期间的地震波形数据,由频时分析方法提取了Love波相速度频散曲线,经过反演得到了研究区内的Love波相速度分布.根据Love波纯路径频散,采用线性反演方法对0.25°... 收集了南北地震带区域地震台网中292个地震台站2008年1月至2011年3月期间的地震波形数据,由频时分析方法提取了Love波相速度频散曲线,经过反演得到了研究区内的Love波相速度分布.根据Love波纯路径频散,采用线性反演方法对0.25°×0.25°的网格点进行了一维S波速度结构反演,利用线性插值获取了南北地震带地区的三维S波速度结构.结果显示了松潘—甘孜地体和川滇菱形块体地区的下地壳具有明显的S波低速层分布,该异常分布特征支持解释青藏高原隆升及其地壳物质运移的下地壳流模型.在100至120km深度上,川滇菱形块体西北部呈现较强的S波高速异常,这可能是印度岩石圈板块沿喜马拉雅东构造结下插至该区域所致,该区域下地壳的低速软弱物质与上地幔的高速强硬物质形成了鲜明对比,暗示了地壳和上地幔可能具有不同的构造运动和变形方式,这为该区域的壳幔动力学解耦提供了条件. 展开更多

关键词 南北地震带 层析成像 LOVE波 互相关函数 相速度

在线阅读 下载PDF 职称材料

用接收函数反演甘肃测震台网下方的S波速度结构 被引量：6

12

作者 安张辉 吴庆举 周民都 《西北地震学报》 CSCD 北大核心 2006年第3期263-267,共5页

利用接收函数方法对甘肃测震台网下方的一维S波速度进行了研究。结果显示甘肃测震台网各个台站下方地壳内部可分为两层：第一层深度为20～25km之间；第二层为Moho界面，其平均深度约为50km，界面上表面的速度约为3．8km／s，界面底部的... 利用接收函数方法对甘肃测震台网下方的一维S波速度进行了研究。结果显示甘肃测震台网各个台站下方地壳内部可分为两层：第一层深度为20～25km之间；第二层为Moho界面，其平均深度约为50km，界面上表面的速度约为3．8km／s，界面底部的速度约为4．5km／s。 展开更多

关键词 甘肃测震台网 接收函数 S波速度结构

在线阅读 下载PDF 职称材料

上海地表软土层、细砂层的地震波反应 被引量：3

13

作者 徐永林 熊里军 《中国地震》 CSCD 北大核心 2003年第1期84-88,共5页

关键词 上海 软土层 细砂层 地震波 传递函数 地震波 地震动

在线阅读 下载PDF 职称材料

地下结构地震响应的计算模型 被引量：8

14

作者 林皋 《力学学报》 EI CSCD 北大核心 2017年第3期528-542,共15页

地震时地下结构在围岩的约束作用下发生变形,其动态特性与地面结构有很大不同.自二十世纪七八十年代以来,地下结构抗震设计与研究取得了很大进展.总的看来,工程设计中普遍采用的计算方法与设计导则大都建立在比较简单假定的基础上,实际... 地震时地下结构在围岩的约束作用下发生变形,其动态特性与地面结构有很大不同.自二十世纪七八十年代以来,地下结构抗震设计与研究取得了很大进展.总的看来,工程设计中普遍采用的计算方法与设计导则大都建立在比较简单假定的基础上,实际的岩土介质条件都是十分复杂的.地下结构抗震研究的近期成果则表现在对地下结构动力分析中的波动散射问题提出了波函数展开法以及边界积分方程方法等多种计算方法.但计算相对复杂,在工程设计中的推广应用有一定困难.此文致力于地下结构计算模型的改进,使之具有良好的计算精度与效率,又便于工程应用.为此,提出了一种地下结构抗震响应分析的新的计算模型.模型具有较广泛的适应性,可以进行河谷、孔洞、地下铁道、隧洞等地下结构的散射与绕射分析.对于复杂层状的地质条件,提出了格林函数求解简便而有效的方法.数值算例论证了方法的精度和效率. 展开更多

关键词 地下结构 地震响应 波的散射 格林函数 层状半空间地基

在线阅读 下载PDF 职称材料

背景噪声频散曲线测定及其在华北地区的应用 被引量：21

15

作者 房立华 吴建平 《地震学报》 CSCD 北大核心 2009年第5期544-554,共11页

利用华北地震台阵垂直分量的观测资料,采用滑动绝对平均方法对资料进行预处理,通过互相关方法从背景噪声中提取瑞利面波的格林函数,开发了群速度频散曲线的自动提取工具,测量了位于华北盆地、燕山隆起和太行山隆起的3条频散曲线,反演得... 利用华北地震台阵垂直分量的观测资料,采用滑动绝对平均方法对资料进行预处理,通过互相关方法从背景噪声中提取瑞利面波的格林函数,开发了群速度频散曲线的自动提取工具,测量了位于华北盆地、燕山隆起和太行山隆起的3条频散曲线,反演得到了3个区域的S波速度结构.研究分析表明,滑动绝对平均方法可以有效降低地震和台站附近干扰源的影响.为了得到可靠的层析成像结果,应计算格林函数的信噪比,选择高信噪比的格林函数测量其频散曲线,进行层析成像反演.当信噪比大于7时,一般都能得到稳定可靠的频散曲线.群速度频散曲线的最大可信周期(Tmax)与台站间距有关,华北地区最大可信周期以不超过台站间距的1/12为宜,周期大于Tmax时不同月份测得的频散曲线变化较大. 展开更多

关键词 噪声 互相关 格林函数 瑞利波群速度 频散曲线 华北

在线阅读 下载PDF 职称材料

面波频散与接收函数联合反演南北地震带北段壳幔速度结构 被引量：24

16

作者 郑晨 丁志峰 宋晓东 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2018年第4期1211-1224,共14页

通过对南北地震带北段区域所布设的676个流动地震台站观测资料进行处理,联合反演面波频散与接收函数数据,获得了研究区内地壳厚度、沉积层厚度的分布情况以及地壳上地幔高分辨率S波速度结构成像结果.反演结果显示研究区地壳厚度从青藏... 通过对南北地震带北段区域所布设的676个流动地震台站观测资料进行处理,联合反演面波频散与接收函数数据,获得了研究区内地壳厚度、沉积层厚度的分布情况以及地壳上地幔高分辨率S波速度结构成像结果.反演结果显示研究区地壳厚度从青藏高原东北缘向外总体逐渐变薄,秦岭造山带地壳厚度较同属青藏高原东北缘的北祁连块体明显减薄;鄂尔多斯盆地及河套盆地分布有非常厚的沉积层,阿拉善块体部分区域也有一定沉积层分布,沉积层与研究区内盆地位置较为一致;松潘—甘孜块体、北祁连造山带等青藏高原东北缘总体表现为S波低速异常;在中下地壳,松潘—甘孜块体下方的低速体比北祁连造山带下方的低速体S波速度值更小、分布深度更浅,更有可能对应于部分熔融的地壳;鄂尔多斯盆地在中下地壳以及上地幔内有着较大范围的高速异常一直延伸到120km以下,而河套盆地地幔只在80km以上部分有着高速异常的分布,此深度可能代表了河套盆地的岩石圈厚度,来自深部地幔的热物质上涌造成了该区域的岩石圈减薄;阿拉善块体在地壳和上地幔都表现出高低速共存的分布特征,暗示阿拉善块体西部岩石圈可能受青藏高原东北缘的挤压作用发生改造. 展开更多

关键词 联合反演 面波频散 接收函数 南北地震带 青藏高原东北缘

在线阅读 下载PDF 职称材料

地震衍射CT中的震源函数与入射场

17

作者 刁顺 杨慧珠 许云 《石油地球物理勘探》 EI CSCD 北大核心 1994年第4期434-440,共7页

本大应用爆炸力学和理论力学对爆炸震源（球状和往状）及集中力震源（落重、可控、敲击）进行了研究。给出了地震衍射CT中所涉及的震源函数形式和相应的入射场解析解与数值结果，并采用弹塑性理论、冲量定理及Laplace变换导出了位移Gree... 本大应用爆炸力学和理论力学对爆炸震源（球状和往状）及集中力震源（落重、可控、敲击）进行了研究。给出了地震衍射CT中所涉及的震源函数形式和相应的入射场解析解与数值结果，并采用弹塑性理论、冲量定理及Laplace变换导出了位移Green函数，进而作褶积运算，最后得到位移入射场表达式及其图表显示。 展开更多

关键词 震源 入射场 地震波 函数 CT

在线阅读 下载PDF 职称材料

地震波动强度变化的数学建模分析与仿真 被引量：2

18

作者 柴瑞帅 《地震工程学报》 CSCD 北大核心 2018年第3期549-554,共6页

传统利用灰色关联分析方法对地震波动强度变化进行数学建模分析与仿真时,对地震波动强度变化的数列进行仿真分析时,忽略了地震波动强度的时间属性对结果的影响,导致分析结果准确性较低。本论述提出新的地震波动强度变化数学建模分析与... 传统利用灰色关联分析方法对地震波动强度变化进行数学建模分析与仿真时,对地震波动强度变化的数列进行仿真分析时,忽略了地震波动强度的时间属性对结果的影响,导致分析结果准确性较低。本论述提出新的地震波动强度变化数学建模分析与仿真方法,通过地震波动强度序列的经验分布确定门限自回归模型的门限值,依据该门限值、AIC最小准则以及最小残差平方等方法获取地震波动强度序列的门限自回归模型,分析自回归模型的极限环和振荡的属性特点,得到地震波动强度变化的初步数值模拟结果。本论述构建了基于均生函数的地震波动强度序列的数学模型,通过均生函数数学建模方法拟合地震波动强度时间序列,依据时间序列基于双评分准则选取拟合周期,实现地震波动强度的数值仿真。实验结果表明,所提方法对地震波动强度变化模型具有较高的准确性和稳定性。 展开更多

关键词 地震波动强度 数学建模 地震波动强度序列 均生函数 双评分准则 门限自回归模型

在线阅读 下载PDF 职称材料

根据接收函数反演得到的首都圈地壳上地幔三维S波速度结构 被引量：20

19

作者 王峻 刘启元 +3 位作者 陈九辉 李顺成 郭飚 李昱 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2009年第10期2472-2482,共11页

利用2002-2003年中国地震局地质研究所台阵实验室以唐山大震区为中心布设的40个流动宽频带地震台站和首都圈数字台网的33个宽频带台站的远震数据,采用接收函数非线性反演方法得到其中72个宽频带台站下方60 km深度范围内的S波速度结构.... 利用2002-2003年中国地震局地质研究所台阵实验室以唐山大震区为中心布设的40个流动宽频带地震台站和首都圈数字台网的33个宽频带台站的远震数据,采用接收函数非线性反演方法得到其中72个宽频带台站下方60 km深度范围内的S波速度结构.根据得到的各台站下方地壳上地幔的S波速度结构,并综合刘启元等(1997)用接收函数非线性反演方法得到的延怀盆地15个宽频带流动台站下方的地壳上地幔S波速度结构模型,给出了39°N-41°N ,114°E-7.19. 5°E区域内沿不同走向、不同深度S波速度分布.由于综合了利用首都圈数字地震台网的宽频带台站以及流动地震台阵的观测数据,本文给出了较前人同类研究空间分辨率更好的结果.结果表明:(1)研究区的速度结构,特别是怀来以东的速度结构十分复杂.在10-20 km深度范围内,研究区地壳具有高速和低速异常块体的交错结构.研究区中上地壳速度结构主要被与张渤地震带大体重合的NW向高速条带和穿越唐山大震区的NE向高速条带所控制,而其中下地壳的速度结构主要为延怀一三河一唐山地区上地幔隆起所控制.(2)研究区内存在若干壳内S波低速体,它们主要分布在唐山,三河及延怀盆地等地区.在这些地区,壳内低速体伴随着壳幔界面的隆起和上地幔顶部速度结构的横向变化.(3)地表断层分布与地壳速度结构分区有较好的相关性,表明断层对不同块体有明显的控制作用.其中,宝坻断裂,香河断裂和唐山断裂均为超壳断裂.(4)首都圈内大地震的分布与壳内低速体及上地幔顶部的速度结构有密切关系.对于唐山大地震的成因,仅考虑板块作用引起的水平应力场是不够的,有必要充分重视由于上地幔变形引起的地壳垂直变形和上地幔物质侵人造成的热效应. 展开更多

关键词 首都圈 接收函数反演 流动地震台阵 地壳上地慢S波速度结构

在线阅读 下载PDF 职称材料

华北地区基于噪声的瑞利面波群速度层析成像 被引量：122

20

作者 房立华 吴建平 吕作勇 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2009年第3期663-671,共9页

本文利用华北地震科学台阵2007年1～4月份、190个宽频带和10个甚宽带垂直分量的地震噪声数据,通过互相关方法提取了瑞利面波的经验格林函数,用多重滤波方法测量了瑞利面波的群速度频散曲线.我们将研究区域划分为0.5°×0.5°... 本文利用华北地震科学台阵2007年1～4月份、190个宽频带和10个甚宽带垂直分量的地震噪声数据,通过互相关方法提取了瑞利面波的经验格林函数,用多重滤波方法测量了瑞利面波的群速度频散曲线.我们将研究区域划分为0.5°×0.5°的网格,利用噪声层析成像方法得到了研究区域7 s、12 s、16 s、23 s 的瑞利面波群速度分布图像,所得结果较好地揭示了地壳内部、尤其是浅部地壳的横向速度变化.研究表明,短周期的群速度分布同地表地质结构、地形密切相关;华北地区的地壳结构具有明显的横向不均匀性,华北盆地及山间的沉积盆地显示出低速异常,而基岩广泛出露的太行山和燕山隆起区,呈现高速异常;多数强震(M≥6.0级)都发生在高群速度与低群速度的过渡地带. 展开更多

关键词 地震噪声 互相关 格林函数 群速度频散 瑞利面波层析成像 华北

在线阅读 下载PDF 职称材料