卢龙地震动场的预测

分别用经验格林函数方法和随机方法对卢龙ML5.7地震的地震动进行模拟对比研究,结果显示:用经验格林函数方法模拟得到的波形和反应谱都与原始记录接近,用该方法模拟的地震动,无需掌握场地的具体情况;对于无小震记录的地区,如果知道了其场地情况与附近有小震记录地区相似,可以用该小震记录作为经验格林函数模拟地震动。对于不能用经验格林函数方法模拟地震动的地区,可以用随机方法模拟地震动。用随机方法模拟地震动具有简单、方便、无需小震记录的优点。根据卢龙场地具体条件,尝试用两种方法结合预测了卢龙地震加速度场,结果比较好地反映了断层周围加速度的峰值分布规律。     

用经验格林函数法和随机方法模拟卢龙地震的比较

在介绍用经验格林函数方法和随机方法模拟地震动的基础上,分别用这两种方法模拟了卢龙5．7级地震。模拟结果表明,用经验格林函数方法模拟地震动的结果与地震记录符合得很好。用随机方法模拟地震动在高频部分与地震动记录符合得比较好,而在中、低频部分符合得不太好。两种方法比较,经验格林函数方法模拟的结果更符合实际记录。但随机方法具有计算简单、无需小震记录的优点,因此发展随机的方法是有意义的。随机方法中中、低频部分不太符合的原因是随机方法中的很多经验关系是在美国大陆得到的。因此建议加强对中国大陆经验关系的研究,并把中国大陆的经验关系加入到随机方法中,使随机的方法成为中国大陆普适的模拟地震动的方法。     

用经验格林函数方法模拟汶川M_W7.9地震的强地震动

汶川M_W7.9地震的发生震惊了中国乃至全世界,不仅仅因为该地震震级巨大,还因为其可怕的伤亡。评估汶川M_W7.9地震的复杂性及破坏力需对其地震特征进行研究。本文运用经验格林函数方法(EGF)模拟了地震的强地震动。首先,选取了4个震级为M_W4.9和M_W5.0的余震,对这些余震的运动利用格林函数进行了模拟。通过模拟,可以估算出包括断层规模和应力降在内的余震震源参数。然后,通过试错法正演模拟得出M_W7.9地震凹凸体(强震动生成区)的震源模型。为验证该震源模型的准确性,运用经验格林函数方法模拟合成了9个台站的强地震动。将模拟得出的强地震动与观测记录到的加速度波形、速度波形和拟加速度反应谱相比较。尽管有一些差异,但在大部分台站,合成的强地震动与观测记录结果相吻合。在模拟得出的震源模型中,断层面上3个凹凸体的应力降约为10MPa,而凹凸体面积与破裂面积比为0.24,接近经验关系式得出的结果0.22。     

经验格林函数方法的理论探讨

通过比较经验格林函数和理论格林函数模拟地震动方法,认为经验格林函数方法能比理论格林函数更好地模拟地震动。同时亦指出用经验格林函数方法模拟地震动时出现两个问题:(1)模拟的结果高频不稳定。(2)模拟的结果有一个"人造周期"。地震工作者提出选择适当的频响函数,比较有效地解决了模拟结果高频不稳定的问题;还提出了完备经验格林函数集的概念,以提高模拟结果与记录的符合程度。但完备经验格林函数集提出了每个大震子断层上都有小震记录的过于苛刻的要求,使得完备经验格林函数集的获得只能靠理论和经验关系结合的办法。本文通过对大震断层滑动速度函数的理论分析认为:可以用多个不同滑动速度和不同上升时间的小震合成经验格林函数,用这种经验格林函数合成的地震动不会产生"人造周期"。但需要加强对小震记录的反演工作并获得丰富的小震参数。     

用改进的经验格林函数方法模拟唐山地震动

建筑物的抗震设防需要尽可能地掌握未来大地震强震动记录信息,但大地震强震动记录的匮乏阻碍了抗震设防实践的发展。经验格林函数方法作为模拟地震动的主要方法,可以提供可靠的大地震强震动记录,但也存在着许多问题,如缺乏对大地震断层滑动分布不均匀的描述、用经验确定小震数目、模拟方法受到大小地震相似条件的限制等。文中对上述经验格林函数方法存在的问题进行了研究,改进的经验格林函数方法,有效地解决了上述问题。并用其对唐山大地震进行了模拟,并把模拟的地震动时程和反应谱与实际记录相比较,发现用改进方法模拟的地震动加速度反应谱比用未改进方法模拟结果更接近实际的地震动记录加速度反应谱。由此说明改进的经验格林函数可更准确的模拟地震动。     

九寨沟MS7.0地震强地震动模拟及漳扎镇地震动强度预测

破坏性地震强度预测可用于工程领域抗震设防以及地震危险性分析评估,是防震减灾中一项很重要的基础工作.为了再现九寨沟地震的地震动强度,评估缺失强震记录的九寨章扎台站的地震动强度,本文用经验格林函数法对九寨沟地震进行了数值模拟.选取了震源周边地震动峰值加速度超过10 Gal的10个强震台站进行模拟.因未得到九寨沟地震的余震,初次尝试将汶川地震和定西地震的余震作为格林函数模拟九寨沟地震.模拟结果整体上可以反映各台站地震动的强度特征,尤其是地震动高频成份拟合较好.模拟值的地震动峰值加速度、时程数据、反应谱等与观测值拟合较好.预测结果显示漳扎镇的地震动峰值加速度值约为180~200 Gal.预测结果也表明在缺少大震的余震记录时,经验格林函数法使用其他大震的余震同样可以再现目标地震的强度特征.本研究也为经验格林函数方法在缺乏小震记录地区的使用积累了经验.最后总结了格林函数的选取标准,为经验格林函数方法来预测未来强震动特征积累了经验.     

用格林函数方法模拟卢龙地震

建筑物的抗震设防需要地震动的输入,但往往在抗震设防区缺少地震记录,这就需要地震动的模拟.在总结现有经验格林函数方法的基础上,提出了改进的经验格林函数方法,并用该方法对卢龙地震进行了模拟.模拟的结果表明:模拟的结果无论是地震动加速度时程还是反应谱与实际地震记录都符合得比以往好.说明改进的格林函数方法是模拟地震动的有效方法,可以在工程实践中应用推广.     

2021年5月21日云南漾濞MS6.4地震随机有限断层三维地震动模拟

2021年5月21日云南省大理州漾濞县发生MS6.4 地震,震中附近遭受了强烈地震破坏.为预测此次地震的地震动影响场,利用震源运动学破裂随机模型,基于随机有限断层三维地震动模拟方法,给出了此次地震中 28个触发强震动台站的三分量加速度时程模拟记录,并结合强震动观测记录,估计了此次地震的地震应力降及震源破裂过程,进一步模拟给出了此次地震中 2823 个虚拟观测点的三分量加速度时程.结果表明,模拟记录的峰值地面加速度(PGA)、峰值地面速度(PGV)与观测值接近,并体现了地震动峰值的衰减规律、近场饱和效应和破裂方向性效应;模拟与观测记录的 5%阻尼比拟加速度反应谱(PSA)的幅值接近、谱形相似,在0.05～10 s周期段,模拟记录可以很好地预测地震动.基于三分量模拟记录给出了漾濞MS 6.4 地震的仪器测定地震烈度图,与云南省地震局发布的烈度图接近,极震区烈度最高可达Ⅷ度,震源破裂方向性导致震中 SE方向的烈度普遍高于 NW方向,受局部场地条件影响沿洱海西侧出现高烈度异常区.     

经验格林函数方法模拟地震动研究

地震动模拟研究的目的是为建筑物的抗震设计提供未来地震动的时程或反应谱。本研究在总结了经验格林函数方法和随机方法模拟地震动的基础上,对所存在的问题进行了系统研究,得到如下成果:1改进了经验格林函数方法,提高了经验格林函数方法的严密性用经验格林函数方法模拟地震动的结果在高频部分偏低的部分原因是由于假设大小地震的滑动分布都是相     

基于改进的经验格林函数方法的地震动模拟

改进的经验格林函数方法弃掉以往完全用经验方法确定用于合成大震的小震数目的方法,考虑了大地震断层各个子断层位错的不均匀性,使得经验格林函数方法的物理意义更加清楚。用改进的经验格林函数方法,以汶川MS6.4余震记录作为经验格林函数模拟了九寨白河台、松潘安宏台、皋兰台和成都中和4个台站的加速度时程,与以往的经验格林函数方法比较,模拟结果的反应谱与记录符合得更好,证明了改进的经验格林函数方法可以反映震源位错的复杂分布,可以在工程地震中用来估计场地地震动过程。     

中强地震随机有限断层模型应力降参数的确定方法

随机有限断层模型是模拟地震动加速度时程的一个重要工具.但将其应用于中强地震时,由于震源信息的准确性较差从而使模型参数具有较大的不确定性.尤其针对其中最为关键的应力降参数,目前相关研究较为缺乏且尚未形成系统的确定方法.本文基于美国Little Skull Mountain Mw5.6级地震2个近场台站记录的地震动模拟,详细研究了采用随机有限断层法拟合中强地震地震动伪加速度反应谱（PSA）来确定应力降参数值的方法,并在计算应力降时引入了其它震源参数的不确定性,随后对此方法的可行性进行了验证.研究表明：采用不同频段反应谱残差和计算得到的应力降值差别较大,确定中强地震应力降较为合适的反应谱频段是中高频,采用该频段确定的应力降参数值模拟的反应谱和峰值加速度与实际记录较为符合;脉冲子断层百分比、断层长宽、倾角和深度等震源参数按截断的正态分布或均匀分布随机抽样,拟合得到的应力降参数值与通过实际震源模型参数得到的值相近.以上研究结果对确定一个区域中强地震应力降或中强地震近场强震动模拟研究提供了更进一步的研究方法和研究方向.     

利用经验格林函数谱比法估计2013年芦山地震序列的应力降

本文以芦山地震余震为例,分析了基于强震动观测记录的经验格林函数（EGF）谱比法估计地震的拐角频率和应力降的可行性。首先给出能够可靠估计地震拐角频率和应力降的经验格林函数谱比曲线的质量标准;然后利用其估算出17次M_L3.8—5.4芦山强余震的拐角频率;最后,参考其它研究给出的地震矩震级,计算出地震应力降。结果显示,芦山强余震的拐角频率主要介于1.0—2.0 Hz之间,应力降平均值为9.98 MPa,且地震应力降表现出明显的震级相关性。      

由强震动数据分析芦山地震地面运动持时及周期特征

以芦山地震震中距100km内的20个台站60条原始强震记录为依据,研究芦山地震动卓越周期、反应谱峰值周期、地震动持续时间和震源持续时间。研究发现:EW、SN和UD3个方向卓越周期的变化范围为0.013~0.275s;5%阻尼比反应谱峰值周期变化范围为0.03~0.65s;芦山地震动90%持续时间变化范围为5.1~35.9s,震源持续时间在EW、SN和UD 3个方向的均值分别约为6.41s、6.05s和5.47s。以台站地震动持时为依据,回归地震动持续时间的空间变化曲线,并与Bommer(2009)地震动持时预测曲线进行比较,发现芦山地震动持续时间高于Bommer(2009)预测曲线,表明芦山地震中的震源作用时间大于目前的预测平均值。芦山地震动持时和震源持续时间在NE向有较明显的方向性,尤其是震源持续时间;而该方向上的余震分布弱于震中SW方向,表明发震构造在主震发生过程中NE方向能量释放较多,所以余震分布较弱。地震动持时和震源持时无上盘效应可能与隐伏发震构造相关。     

汶川地震强地面运动模拟

基于确定性震源模型的方法主要用于计算低频(11 Hz)地震动常用经验格林函数法或随机方法,对低频地震动模拟不够准确．本文在确定性震源模型方法基础上,尝试采用分解给定的震源模型的方法来模拟宽频带(0.1——10 Hz)强地面运动,即采用分级离散断层面和分解断层面破裂单元上升时间的方法,增加震源时间函数中的高频信号,从而避免了对地震记录丰富程度和准确性的依赖．文中模拟计算了汶川MS8.0地震在8个地震观测台的地震动,将模拟结果和观测记录进行了加速度时程曲线和傅里叶振幅谱的对比分析．对比结果显示,模拟估计的地震动峰值加速度和持续时间与观测记录的数据基本在plusmn;50%的精度范围内相同,傅里叶振幅谱显示模拟结果有得到10 Hz左右的高频成分. 四川盆地中的台站模拟结果高频衰减比观测记录要快,原因是模拟过程没有考虑场地效应.对强地震动模拟还是要综合考虑震源、传播路径和场地的影响．研究结果表明,此改变震源输入的确定性方法可应用于模拟近断层宽频强地面运动．      

经验格林函数方法模拟地震动研究现状

经过30多年的发展,格林已经成为被地震工作者普遍接受的模拟地震动的方法.回顾发展历史,该方法可以分为3个阶段:第1阶段是经验格林函数理论建立阶段;第二阶段是经验格林函数方法逐步得到完善阶段;第三个阶段是经验格林函数方法的广泛应用阶段,该阶段提出了用该方法预测未来地震动的一般步骤,标志着用经验格林函数方法预测未来地震动,为建筑物的抗震设防提供依据已经到了比较成熟阶段.最后,对用经验格林函数方法和随机方法模拟地震动的优势和不足进行了比较.     

大震近场地震动数值模拟的不确定性研究

地震动数值模拟是研究地震发生、发展特征和规律性的基础研究,涉及震源、传播路径、场地三方面的一系列基本问题。研究结果可以直接用于强地震动基本参数的标定和地震动的工程预测以及地震危险性分析,是防震减灾中一项很重要的基础工作。当前考虑断层模型中的几个重要震源参数模拟未来地震动成为热点。例如,不管是随机有限断层方法,还是经验格林函数方法,都特别重视震源参数的选取,但当前对于强地面运动的数值模拟大多是对已发生地震的反演,反演所需的震源参数多是根据地震发生后的观测资料、经验关系推断或者人为的假定得到的。地震动特征预测的好坏,关键也在于震源参数的选取是否可靠,针对工程应用,对未来近断层地面运动的预测应当考虑参数选取的不确定性。强震动数值模拟方法已经比较成熟,比如随机有限断层法、经验格林函数法、混合方法等。但这些模拟地震动的方法多是确定性方法,对模拟过程中各个环节的不确定性因素考虑较少,尤其是断层震源参数的选取都是精确的数值,模拟结果自然不能综合表征未来地震发生时的地震动特征。另外一种地震动特性表征形式是利用地震动预测方程来表现参数的衰减关系,但是预测方程也会受观测数据、模型选取、随机因素的影响,不确定性无法避免。在地震动预测中每个步骤都带有很大的不确定性,主要包括认知的不确定性和随机的不确定性。震源模型的建立依赖各种专家的观点和意见,这种观点和意见往往相互差别很大,有很大的人为不确定性。地震资料的完整性、可靠性,以及理论模型等的建立也会造成很大的不确定性。因此,本论文的研究内容重点在预测未来大震近场地震动时震源环节的不确定性因素,系统地研究影响大震近场地震动数值模拟不确定性的影响因素,以期建立描述影响大震近场地震动数值模拟不确定性影响因素的模型,并将考虑不确定性因素的震源模型应用在破坏性大地震的地震动特征预测中去。针对地震发生时震源参数具有的随机不确定性,通过统计学方法得到震源参数间的经验关系。同时由于地球内部的不可入性以及观测技术手段的限制,地震震源模型还存在很多的认知不确定性因素,我们采用逻辑树方法处理认知不确定性因素,引入强震生成区的概念,重点分析Asperity的认知不确定性因素并在此基础上建立Asperity震源模型,取多种方案结果的最优值作为地震动的预测结果。论文具体的研究内容有以下5点:(1)在Hiroe Miyake工作的基础上,用经验格林函数法模拟了1997年日本九州鹿儿岛县的MJMA6.5地震,验证了经验格林函数法模拟强地面运动的有效性。数据使用的是日本K-NET强震台站数据。选择此次地震以及用日本K-NET数据的原因是日本的地震记录比较丰富,台站收集到的地震记录质量相对较好。同时分析了地震动模拟值中在工程领域使用较多的参数,表明模拟结果较好地反映了真实的地震动记录。证明经验格林函数预测未来强地震动是切实可行的。(2)本论文在大量地震记录及文献调研的基础上,运用统计学方法,通过统计回归分析得到地震密集区域震源参数间的经验关系,重点统计了包含龙门山断裂带在内的中国西南地区震源参数间的经验关系。为了得到更加适应于中国大陆区域的局部震源参数的经验关系,我们将地震样本局限在一定的区域内,缩小统计区域的范围,同时得到多个Asperity等局部震源参数的经验关系。(3)研究和分析各不确定性因素对地震动结果的影响程度对提高地震动预测的可靠性和准确性意义重大。本论文以2016年发生的日本熊本7.3级地震为例,通过参数敏感性分析技术筛选了对拟合结果有决定性影响的关键模型参数。未来地震动预测时会考虑在这几个不确定性因素中加入更高的权重,以提高地震动预测的可信度和准确性。(4)基于汶川地震资料,研究分析Asperity的认知不确定性对地震动预测的影响,运用逻辑树方法建立了考虑不确定性因素的Asperity震源模型。采用Hiroe Miyake提出的强震生成区的概念,将Asperity区域等同于强震生成区域,主要从Asperity的面积、数量、应力降、上升时间等方面进行了认知不确定性的分析研究,得到了针对强震的可能性较大的Asperity的面积和数量的最佳搭配形式。本文的创新点将概率地震危险性分析中不确定性因素的处理方法应用到未来大震近场地震动数值模拟中去。随机不确定性因素用统计学方法处理,认知不确定性用逻辑树方法处理。模拟结果能更好地表征未来大震近场的地震动特征,地震动参数的模拟结果在未来抗震设防及灾害预防中会有更好、更高的参考价值。考虑不确定性因素预测未来破坏性大地震,从而定量地评价发生地震时近场地面运动的综合特征,为防震减灾提供科学的依据,为工程设计决策提供参考依据。     

随机有限断层方法的改进和芦山地震动模拟

随着随机有限断层方法被逐步推广至宽频带地震动的震前预测与震后模拟,该方法存在的2类问题凸显:(1)模拟地震动的频率非平稳特性缺失;(2)高低频地震动的模拟精度不一致。文中探索上述问题的改进措施:(1)引入区域地震波传播群速度模型以生成频率非平稳地震动相位;(2)引入王国新改进震源谱代替Brune震源谱以实现高低频地震动辐射水平的调整。将改进后的随机有限断层方法应用于2014年芦山地震(M_w6.7)的强地震动模拟,模拟结果验证了改进措施的有效性:新引入相位使得模拟记录的频率非平稳特性趋近于真实记录;新引入震源谱模型在一定程度上可实现高低频地震动的一致精度模拟。     

1995年7月20日怀来盆地ML4.1地震序列震源参数的精确测定

1995年7月20日在北京西北的怀来盆地发生了一次ML4.1地震．这次地震震中位于40.326N，115.448E，震源深度5.5ｋｍ．在此主震之后该地区微震活动变得十分活跃．中-欧合作怀来数字地震台网记录并精确定位了这一地震序列．以８次余震的记录作为经验格林函数提取了ML4.1地震的震源时间函数，并通过叠加得到了信噪比较高的平均结果．结果表明，这次ML4.1地震由一强一弱两次事件组成．各个台站的震源时间函数显示出明显的地震多普勒效应．用试错法得到了ML4.1地震的第一次事件的破裂长度为0.44ｋｍ，破裂速度为4.0km/s；第二次事件破裂长度为0.25km，破裂速度3.0km/s．两次事件的破裂传播方向与破裂面走向的夹角（逆时针为正）分别为140和90，两次事件相距0.57km，第二次破裂发生于第一次破裂开始之后0.09ｓ．用单个经验格林函数提取了怀来盆地ML4.1地震序列中另外13次ML2.1地震的震源时间函数．结果表明，这些小震均由单次事件组成，震源时间函数宽度为0.05～0.16ｓ．用远场地动位移频谱测量法得到了0.9ML4.1的２５次地震的地震矩、应力降和破裂半径．地震矩和应力降都呈现出随震级而单调递增的规律性变化．在0.9ML2.4范围内，求得的破裂半径与震级没有显示出明显的相关性，可以认为，在这样小的震级范围内用本文所用资料已无法准确分辨小震的破裂尺度．      

基于加速度包络的芦山地震高频辐射反演

基于经验格林函数和加速度包络,利用差分进化法反演了芦山地震断层面上高频地震波(1~10 Hz)辐射区域,其高频地震波的辐射强度利用震源时间函数的幅值表示。结果表明,芦山地震断层面上高频辐射分布很不均匀,其中凹凸体部分是高低频辐射最强的区域,另外,凹凸体周边以及断层破裂停止的区域也存在较强的高频地震波辐射,而这些区域的低频地震波辐射则很弱。这表明芦山地震地震波辐射过程是极其复杂的,研究结果对宽频带地震动模拟芦山地震震源模型的建立具有一定的应用价值。     

考虑凹凸体理论的经验格林函数方法

以15次地震资料为基础,回归得到断层上所有凹凸体地震矩与凹凸体面积的关系;研究经验格林函数方法中凹凸体子源的应力降,提出基于断层凹凸体理论的经验格林函数方法;以唐山地震余震（Ms=6.9）作为经验格林函数,合成1976年唐山地震（Ms=7.8）的地震动,比较合成结果与观测记录表明,该方法符合断层破裂过程,计算结果可靠．