福建地区爆破与地震的波谱差异-应用昆虫学报.DOC

福建地区爆破与地震的波谱差异 蔡杏辉1) 曹轶1)许振栋2) 陈惠芳1) 1）中国福州 350003 福建省地震局 2）中国平潭 350400 福建省地震局平潭地震台 摘要：应用福建地震台网记录的爆破与地震的数字资料，采用波谱分析的方法分析发生在同一区域的爆破与地震波谱特征。结果表明，爆破与地震的拐角频率存在差异，据此得出本区域爆破识别的定量识别判据，为今后爆破的识别提供依据。 关键词：爆破；波谱分析；拐角频率；评定阈值； 0 引言 随着社会经济的发展，爆破作业频繁，人工爆破和地震的识别贯穿于地震台网日常工作，有的爆破的爆破特征明显，较容易识别，有的爆破波形复杂，不具有常见的一些爆破特征，并且在同一区域时有微震发生，分析不当而产生的地震目录资料会为使用者带来一定的伪信息，对地震活动性分析研究有可能产生影响。目前福建台网在爆破与地震的识别上主要依靠分析人员的经验，均属于定性分析阶段，缺少定量判据。随着数字台网建立，产出了丰富的数字地震波形数据为定量识别爆破与天然地震创造良好条件。地震波频谱分析方法被认为是识别天然地震和人工地震最有希望的方法（吴忠良，1994），已有作者（张萍等，2005）利用波谱分析方法进行爆破的识别研究。本文对福建台网记录的爆破与地震资料进行波谱分析，了解本地区爆破与地震波谱差异情况，为爆破的识别提供定量的识别判据。 1 地震波谱计算方法 地震波记录图的波谱特征受震源和射线路径传播介质的影响，各频段内能量的衰减也不相同。把一个时间域的函数，经过一定的变换，转换到频率上所得到图形称为谱。把一个记录放到频率域中进行分析和处理就是所谓的谱分析。 频率域方法把记录的地震波f(t)进行傅里叶变换 根据振幅谱（巴特，1978） 上式中为实数，为频率。对于类型不同的震源取不同的时间常数1/。 采用布龙圆盘位错模式给出地震事件的高低频谱值、拐角频率等震源参数，得出远场地震P、S波位移振幅谱表达式（郑治真，1983） Ω 在低频端 Ω(0) 作者简介：蔡杏辉(1976—)，男，福建永安人，本科学历，工程师.主要从事地震监测工作.E-mail. 课题项目：福建省地震局科研基金项目（SF201105） 拐角频率值为 其中Ω(表示振幅谱，Ω(0)为低频谱值，为剪切模量，R为震源距，为应力降， r为震源的等效圆位错半径，为地震矩，V为P波或S波速度，为波谱的拐角 频率。按照上式处理地震和爆破事件的位移振幅谱方法，计算拐角频率值。 2 资料选取与处理 2.1 资料选取 地震波的频谱特性与震源到接收台站间的射线路径有关，因此研究选取地震与爆破事件地点要相近。从福建数字台网运行期间记录的爆破事件和地震事件中，选取64个爆破、64个地震作为研究样本进行分区研究，其中安溪地区20个，清流地区20个，永泰地区20个，古田地区24个，仙游地区24个，龙岩地区20个。爆破事件的震级范围为ML1.0-2.0, 地震事件的震级范围为ML1.0-2.5。 2.2波谱计算过程 采用PITSA软件进行波谱分析，该软件是用于数字地震波形分析的实用软件，并且功能较强。PITSA软件可以处理ISAM,SUDS,ASCII格式的数据文件，而目前台网波形事件都是以seed格式存储，因此在波谱分析前要先将选取事件进行格式转换，转换成ASCII格式文件。S波数据资料取爆破和地震都共有且记录较好台的南北向和东西向，窗长取S波初始震相到最大振幅间一定长度的波形数据，拐角频率取两水平向平均结果。波形资料处理步骤：对原始速度记录进行去倾斜及水平基线校正后，做积分处理；积分运算后，采用4阶巴特沃斯带通滤波器进行数字滤波处理，最后对数据取余弦窗。得出位移谱图后，采用人工拟合的方法得到S波拐角频率值。S波拐角频率的确定原则：波谱高频段与低频段的趋势渐近线交点处的频率值（陈运泰等，2000），也就是说振幅有明显衰减变化时所对应的频率。图1为波谱主要参数示意图。 图1 波谱主要参数 Fig.1 Spectrum of option legend 3 计算结果及分析 3.1计算结果 运用上述方法，对各区域爆破和地震事件的S波波谱特征进行分析，得到各区域爆破和地震平均S波拐角频率值，分析结果见表1。表1的S波拐角频率值取爆破与地震都有记录的2个近台水平向的平均结果。 表1 各区域爆破与地震S波拐角频率计算结果 Table 1 Vary area explosion and earthquake of S waveform corner frequency computed result 地区 事件类型 事件数 平均fc/Hz Δfc/Hz 安溪 爆破 10 8.7