衰减、场地响应等地震波传播信息综合研究

朱新运

1）中国地震局地球物理研究所，北京100081

2）浙江省地震局，杭州310013

衰减、场地响应等地震波传播信息综合研究

朱新运1，2）

1）中国地震局地球物理研究所，北京100081

2）浙江省地震局，杭州310013

地震激发的波，经过传播路径、仪器、场地等因素的调制作用成为台站地震记录，即台站记录是震源、路径、仪器、场地等多种因素的耦合，而基于地震记录研究震源信息、路径衰减、场地响应等则是逆向的解耦过程。基于地震学理论、通过各种计算技术，从地震波形资料中解析出关于地震震源力学参数、地震波传播介质及台站场地响应信息是地震研究的重要内容，也是认识地震孕育、发生、发展过程，了解地震孕育背景构造特征等的有效途径。本文在地震学前沿理论方法的基础上使用最新计算机技术开展了以下3个方面的工作。

本文第一部分，在逆向双台谱比法的基础上，提出了基于谱比法的地震波衰减与场地响应联合反演方法。与逆向双台谱比法相比，该方法对数据条件限制降低，增加了可用数据量，除了能获得地震波衰减参数外，还同时获得了场地响应参数。使用该方法，选取首都圈台站记录的同区域小震资料，评估首都圈台站场地响应，计算相应研究区地震波衰减参数。研究表明，在1～6Hz频率范围内，基岩台场地响应小、平坦且稳定，土层台对低频率波放大作用高于对高频率波的放大作用，地震波衰减参数与区域地震活动性相关。同样使用该方法，选取浙江省区域地震台网台站记录的台湾地震波形资料，评估浙江省丘陵地带基岩台站场地响应，计算了穿越台湾海峡较单一路径Lg波的衰减参数，通过区域地震记录震级偏离情况对该研究场地响应结果的可靠性进行了检验。

第二部分基于弱散射的尾波理论开展地震波尾波衰减研究。文章详细阐述了基于Aki模型及Sato模型的尾波衰减参数计算的原理方法，定量分析了不同滤波器及相关性能、地震波噪声、数据采样、尾波流逝时间等对计算结果的影响。根据计算特点，针对Aki模型在数据采样方面提出了动态采样方法，针对Sato模型提出了最佳线形段自动判断算法。研究表明，利用动态采样代替固定窗长的静态采样方法，在拟合数据段较短的情况下，结果更稳定。

利用浙江、宁夏及首都圈台网近年来记录的数字地震波资料，研究相应地区地震波尾波衰减特征，并结合地震活动情况，研究了地震波衰减参数与区域地震活动的关系，对研究结果的可靠性进行了论证。研究发现，Q0的大小除了与射线穿过区域的地震活动性、地震地质条件有关外，还与台站位置所处的局部地形及地震地质构造有关。地震活动程度越高、地震波衰减越快，衰减参数Q0值越小。研究中，我们发现地震波衰减参数Q0与η具有负相关关系。小震序列对相应区域介质衰减的影响表现为，衰减参数η趋势性降低，而Q0不具有趋势性特征。

根据数据条件，使用基于Sato模型计算的衰减参数Q0构建区域衰减参数图像，分析华北区域历史大震活动与衰减参数图像的关系，结果表明，研究区MS7以上地震均发生在低Q0区或Q0高梯度递变区。

第三部分为综合研究部分。在这一部分，使用浙江珊溪水库地震序列波形资料，开展了地震定位、小震震源参数、地震波衰减、综合断层面解等研究。在地震定位方面提出了在虚波走时方程基础上的遗传算法定位，由此得到的定位结果与浙江省台网用hypo81方法定位结果相比，走时残差大幅降低，定位质量明显提高。

基于虚波走时方程，将非线性方程线性化，求解获得震源位置及虚波速度初值，之后计算初始走时残差，由走时残差的标准偏差的2倍计算在经度、纬度、深度及虚波速度产生的最大误差，构造遗传种群，利用遗传算法进行二次定位，获得地震定位的最终结果。利用该方法，对浙江珊溪水库震群进行重新定位，结果表明：地震基本上集中发生在一个近乎垂直的北西向的断层面内，断层面走向136°，倾角87°。小震综合面解获得的结果与定位结果基本一致；震源深度平均4.7km，最大深度9.5km，最浅1.7km，地震震中呈北西向狭窄条带分布，沿北西方向地震震源深度趋势性增大，在条带西北端有大约3.5km的地震不连续空段；利用Borun模型计算的小震震源参数结果表明，该震群地震震源参数没有表现出明显有别于构造地震的特征，地震应力降多在0.33MPa左右，平均应力降为0.88MPa，根据应力降空间分布的情况判断，该地震条带西北端地震活动不连续段处于应力降低值区。

根据数据处理要求，本研究自主开发了界面友好、功能完善，数据输入、输出及参数配置操作方便实用的数字地震波处理软件，软件运行节省内存、速度快、流程可控。

衰减参数；场地响应；联合反演；遗传算法

（作者电子信箱，朱新运：zhuxinyun＠126.com）

P315.3＋1；

A；

10.3969／j.issn.0235-4975.2012.01.012