地电场多极距观测装置系统与文安MemSsub5.1地震前首都圈地电场异常研究

马钦忠

马钦忠. 2008: 地电场多极距观测装置系统与文安MemSsub5.1地震前首都圈地电场异常研究. 地震学报, 30(6): 615-625.
引用本文: 马钦忠. 2008: 地电场多极距观测装置系统与文安MemSsub5.1地震前首都圈地电场异常研究. 地震学报, 30(6): 615-625.

地电场多极距观测装置系统与文安MemSsub5.1地震前首都圈地电场异常研究

  • 摘要: 2006年7月4日河北文安县发生了MS5.1地震,有6个地电场观测台布设在北京周围。本文阐述了地电场多极距观测系统的噪音识别原理及其在地电场数据分析中的应用方法。通过对首都圈6个地电场台观测资料的研究分析,在所选出的众多震前地电场异常信号中逐一排除,最后选定的5组地电场异常信号基本上被认为是震前地电场前兆异常信号。结果表明:①文安MS5.1地震前存在着地电场异常前兆信号;②文安地震前地电场异常前兆信号的出现存在着ldquo;敏感点rdquo;效应;③地电场异常前兆信号在不同的台站存在着不同步出现的现象,并且异常结束至发震时间间隔也不同;④离震中区较近的宝坻台记录到的震前地电场异常幅度较小,而离震中区较远的昌黎台记录到的震前地电场异常幅度则较大。
  • 产生地震动空间变化的原因多且复杂,与之密切相关的因素主要有震源性质(包括震源深度、断层的有限尺度和破裂方式)、传播途径和场地条件等(Abrahamson,1993der Kiureghian,1996),而地震动的这种空间变化对生命线工程如大坝、大跨桥梁、高速铁路、地铁和地下管线等长、大结构的地震反应有很大影响(Park et al,2009Chopra,Wang,2010Bi et al,2013Wu et al,2016Li et al,2017Liu,Gao,2018)。目前,研究地震动空间变化最常用的方法是相干函数法,各国研究人员提出了大量的描述地震动空间变化的相干函数模型(Zerva,2009),这些模型中的相干系数主要是基于高密度强震观测台网(台阵)的加速度记录统计拟合得到(Loh,1985Harichandran,Vanmarcke,1986Hao,1989Abrahamson et al,1991屈铁军等,1996刘先明等,2004李英民等,2013)。相干函数模型中的两个主要参数为频率f和台站间距离d。统计获得的所有相干函数模型均具有相同的变化规律,即随频率或距离的增大相干函数逐渐减小。在进行相干函数模型参数的拟合时,截止频率的取值对地震动空间相干函数统计结果有很大影响(丁海平等,2018),所以,不同相干函数模型基本给出了相应的频率使用范围。而由于研究地震动空间变化的强震台阵间台站距离d的变化范围较小(台站数量较少,且台站间的距离一般集中在数十米至数千米之间),因此对距离d的范围和不同的间距数量的选取也影响着相干函数的拟合结果(饶威波等, 2018)。例如,对应不同距离d的相干函数,按照假定的相干函数数学模型可以得到各自的拟合参数;而对所有不同距离的相干函数进行拟合,也可以得到一组统一的拟合参数。但是,根据某一特定距离的相干函数得到的拟合(参数)曲线,与根据所有不同距离的相干系数统一回归得到的该距离的相干函数拟合(参数)曲线是不一致的,有时甚至偏差较大。为了减小该偏差,本文拟提出一种修正方法,即对由不同距离得到的相干函数模型中的参数进行二次回归,得到一组针对距离的参数修正表达式,以期提高相干函数模型中参数的拟合精度,并以SMART-1台阵(Abrahamson et al,1987)第5次和第45次地震的加速度作为统计数据,对该方法的可行性进行验证。

    SMART-1台阵是我国台湾设置的大型密集强震观测台网,由37个台站组成(图1),中心台站(C00)设置在三个同心环的中心,同心环半径分别为200 m,1 000 m和2 000 m,沿每个圆环布置12个间隔大致相等的台站,每个台站可同时测得水平和竖向分量的地震记录。其中,I01−I12台站位于半径为200 m的圆环上,M01−M12台站位于半径为1 000 m的圆环上,O01−O12台站位于半径为2 000 m的圆环上。

    图  1  SMART-1地震台阵布置图
    Figure  1.  The layout of SMART-1 accelerograph array

    本文选取SMART-1台站的第5次和第45次地震记录的东西向(EW)和南北向(NS)水平加速度作为统计数据。第5次地震发生于1981年01月29日,震级为ML6.3,震中距30 km,震源深度25 km,截取的S波时间窗为6—11 s,数据采集间隔dt=0.01 s;第45次地震发生于1986年11月14日,震级为ML7.0,震中距79 km,震源深度7 km,截取的S波时间窗为9—29 s,数据采集间隔同样取dt=0.01 s。SMART-1台阵与两次地震的相对位置如图2所示。

    图  2  两次地震震中位置示意图
    Figure  2.  Schematic diagram of the epicenter location of two earthquakes

    基于强震台阵观测资料对地震动空间变化的研究通常是以随机振动理论为基础,采用相干函数研究空间上任意两点地震动的统计相关性,其中最常用的是迟滞相干系数(常简称为相干系数),其表达式为

    $$ \left|{\gamma }_{ jk} ( f{\text{,}}d ) \right|{\text{=}}\frac{\left|{\bar{S}}_{ jk} ( f ) \right|}{\sqrt{{\bar{S}}_{ j} ( f ) {\bar{S}}_{ k} ( f ) }} {\text{,}} $$ (1)

    式中: f为频率;d为台站间距离;$ {\bar S_{ j}} ( f ) $$ {\bar S_{ k}} ( f ) $$ {\bar S_{ jk}} ( f ) $分别为场点 jk的两个时间过程的自功率谱$ {S_{ j}} ( f ) $$ {S_{ k}} ( f ) $和互功率谱$ {S_{ jk}} ( f ) $的平滑谱,本文采用的平滑窗为11阶汉宁(Hamming)窗。

    根据相干函数的变化趋势,研究人员提出过很多相干函数模型(Loh,1985Harichandran, Vanmanmarche,1986Hao,1989Abraham son et al,1991屈铁军等,1996刘先明等,2004李英民等,2013),本文选用Loh模型(Loh,1985)对函数系数进行回归分析,该模型的数学表达式为

    $$ \left| {\gamma ( f{\text{,}} d ) } \right| {\text{=}} \exp [ - ( a{\text{+}} b{f^2} ) d ] {\text{,}} $$ (2)

    式中ab为拟合系数。

    根据台站间的相对位置,本文组合了12种不同距离的台站对。不同的台站间距和对应的台站对数量列于表1。根据式(1)可以计算不同距离的台站对相干系数(包括EW和NS两个水平方向),如图3图4中的实线所示。根据式(2)的相干函数模型,本文采用2种方法对相干函数进行拟合。一种是分别对不同台站间距d的相干函数拟合,得到的拟合参数ab列于表2,对应的拟合相干函数曲线如图3图4中以γ表示的红色实线。另一种是对所有12种不同台站间距d的相干函数进行拟合(图34蓝色虚线),得到的拟合参数 a b表2中对应的12d

    表  1  台站间距和对应的台站对及数量
    Table  1.  Distance between stations and the corresponding number of station pairs
    台站距离d/m台站对台站对个数台站距离d/m台站对台站对个数
    200 C00−I03 4 1 200 I06−M12 2
    282 I03−I06 2 1 732 M06−O09 2
    400 I03−I09 2 1 800 I03−O03 2
    800 I03−M03 2 2 000 C00−O03 4
    980 I03−M06 2 2 200 I03−O09 2
    1 000 C00−M03 4 3 000 M03−O09 2
    下载: 导出CSV 
    | 显示表格
    图  3  不同间距的相干系数γ与根据所有不同间距的相干系数γ12的拟合曲线对比(以第5次地震为例)
    Figure  3.  Comparison between fitting curves with lagged coherency γ of different distances and fitting curves with lagged coherency γ12 of all distance (taking the 5th earthquake as an example)
    (a) d=200 m;(b) d=282 m;(c) d=400 m;(d) d=800 m;(e) d=980 m;(f) d=1 000 m;(g) d=1 200 m; (h) d=1 732 m;(i) d=1 800 m;(j) d=2 000 m;(k) d=2 200 m;(l) d=3 000 m
    图  4  不同间距的相干系数γ与根据所有不同间距的相干系数γ12的拟合曲线对比 (以第45次地震为例)
    Figure  4.  Comparison between fitting curve with lagged coherency γ of different distances and fitting curve with lagged coherency γ12 of all distance(taking the 45th earthquake as an example)
    (a) d=200 m;(b) d=282 m;(c) d=400 m;(d) d=800 m;(e) d=980 m;(f) d=1 000 m;(g) d=1 200 m; (h) d=1 732 m;(i) d=1 800 m;(j) d=2 000 m;(k) d=2 200 m;(l) d=3 000 m
    表  2  基于Loh相干函数模型的拟合系数结果
    Table  2.  Fitting coefficient results based on Loh coherency function model
    台站距离d/m第5次地震 第45次地震
    ab ab
    200 0.476 49 0.000 72 0.299 25 0.001 60
    282 0.328 58 0.000 38 0.359 96 0.001 10
    400 0.300 20 0.000 29 0.461 33 0.000 88
    800 0.305 26 0.000 21 0.408 63 0.000 46
    980 0.325 09 0.000 11 0.449 52 0.000 23
    1 000 0.260 01 0.000 13 0.496 97 0.000 38
    1 200 0.208 47 0.000 08 0.471 04 0.000 16
    1 732 0.193 98 0.000 04 0.475 44 0.000 06
    1 800 0.180 25 0.000 04 0.425 02 0.000 05
    2 000 0.169 69 0.000 06 0.449 15 0.000 07
    2 200 0.161 71 0.000 04 0.374 76 0.000 05
    3 000 0.164 18 0.000 00 0.321 45 0.000 01
    12d 0.189 77 0.000 08 0.357 08 0.000 36
    下载: 导出CSV 
    | 显示表格

    表2可知,不同台站间距d的相干函数拟合参数各不相同,并且与由12种不同台站间距d的所有相干函数拟合得到的参数也有很大差别。对比图3图4中拟合曲线γγ12可以看出,两曲线完全不重合,且在台站间距d的统计范围内,当距离较大(如2200 m,3000 m)或距离较小(如200 m,282 m)时,相干系数的差别更明显。这是因为相干系数具有随间距增大而减小的趋势,当对所有12种不同台站间距d的相干系数进行拟合时,其结果相对更接近所有相干系数的平均值,也就更接近中间距离的相干系数。

    因不同台站距离的相干函数拟合系数不一致,导致对于相同的台站距离,由不同台站间距得到的相干系数拟合曲线与由所有不同台站间距得到的相干系数拟合曲线相比有误差,且对于某些距离误差较大(图34)。因此,在通常采用根据所有相干系数得到的拟合参数进行地震动场的人工合成时,就会出现短距离的模拟地震动相关性偏大,长距离的模拟地震动相关性偏小的现象。由图5可见:拟合系数ab随台站距离d的分布不在一条直线上,且有一定的离散;而由所有不同台站间距得到的相干系数拟合结果并不随距离d发生变化,实际显示为一条直线,因此,由不同的拟合系数ab得到的同一距离的相干系数存在差别的。为了减小因拟合系数的差异而造成拟合曲线γγ12的误差,本节将对不同距离相干系数的拟合系数ab进行修正,即对ab进行二次拟合。根据ab的分布规律,比较了几组ab的拟合函数形式(只与距离d有关),最终选用参数ab的拟合函数表达式分别为

    图  5  第5次(蓝色)和第45次(红色)地震的参数a (a),b (b)值随距离的变化及修正结果
    Figure  5.  Changes of parameters a (a) and b (b) with distance and modified results of the 5th (blue) and the 45th (red) earthquakes
    $$ a {\text{=}} {a_1} {d^{{a_2}}} {\text{+}} {a_3} {\text{,}} $$ (3)
    $$ b{\text{=}}{b}_{1} ( d{{\text{+}}{b}_{2} ) }^{{b}_{3}}{\text{,}} $$ (4)

    式中a1a2a3b1b2b3分别为二次拟合系数,最终得到ab的修正结果列于表3

    表  3  ab二次拟合系数的修正结果
    Table  3.  Modified results of a and b
    地震a1a2a3b1b2b3
    第5次 0.733 53 −0.133 37 −0.482 98 0.000 10 −0.114 25 −0.805 45
    第45次 −0.000 57 −3.410 92 0.432 46 0.000 82 0.573 00 −2.470 90
    下载: 导出CSV 
    | 显示表格

    将修正后的系数ab代入式(2),得到修正后的相干函数模型为

    $$ \left|{\gamma }_{{\rm{s}}}(f{\text{,}} d)\right|{\text{=}}{{\rm{exp}}} { [ - ( {a}_{1} {d}^{{a}_{2}}+{a}_{3} ) {\text{+}}{b}_{1} ( d{{\text{+}}{b}_{2})}^{{b}_{3}}{f}^{2} ] d }{\text{,}} $$ (5)

    式中系数a1a2a3b1b2b3的值列于表3。将表3中修正后的拟合系数代入式(4),可计算得到修正后的相干系数γs,并与先前的两种拟合结果进行对比,结果如图6图7所示。图中红色实线γ为各单独台站距离下相干系数的拟合结果,蓝色虚线γ12为所有12个台站距离的拟合结果,黑色虚线γs为修正后的相干系数拟合结果。从图中的比较可以看出,修正后的相干系数拟合结果γs与各单独台站距离下相干系数的拟合结果γ的符合程度得到了明显提高,保证了相干系数的拟合精度。

    图  6  基于Loh相干函数模型,不同间距的相干系数γ、由所有不同间距的相干系数γ12和修正后的相干系数γs拟合曲线对比(以第5次地震为例)
    Figure  6.  Based on Loh coherency function model,comparison of fitting curve with lagged coherency γ of different distances,fitting curve with lagged coherency γ12 of all distance and fitting curve with lagged coherency γs of modified parameter (taking the 5th earthquake as an example)
    图  7  基于LOH相干函数模型,不同间距的相干系数γ、根据所有不同间距的相干系数γ12和修正后的相干系数γs拟合曲线对比(以第45次地震为例)
    Figure  7.  Based on Loh coherency function model,comparison of fitting curve of lagged coherency γ of different distances,fitting curve of lagged coherency of all distance γ12 and fitting curve with lagged coherency of modified parameter γs (taking the 45th earthquake as an example)

    为了便于应用,对两次相干系数统一进行拟合,可得到由单独台站距离拟合的相干系数(红色实线γ)和由所有12个台站距离拟合的相干系数(蓝色虚线γ12),结果如图8所示。采用上述相同方法,对 a b进行修正,得到的统一修正结果列于表4,而相应的修正后的相干系数拟合结果如图8中黑色虚线γs所示。可以看出,修正后的结果得到了较为明显的改进。

    图  8  不同间距的相干系数经修正前γ、后γs的拟合曲线与根据所有不同间距的相干系数γ12拟合曲线的比较(以两次地震为例)
    Figure  8.  Comparison of fitting curve of lagged coherency γ of different distances,fitting curve of lagged coherency γ12 of all distance and fitting curve with lagged coherency γs of modified parameter (take two earthquakes as examples)
    表  4  a b的统一修正结果
    Table  4.  Unified modified results of a and b
    a1 a2 a3 b1 b2 b3
    −0.022 1.759 9 0.416 3 0.011 9 1.531 3 −4.367
    下载: 导出CSV 
    | 显示表格

    本文采用了SMART-1台阵两次地震的水平向分量加速度记录,计算分析了台站距离d对地震动空间相干函数的影响。首先计算了不同间距台站对的相干函数,然后选用Loh相干函数模型对不同距离的相干函数,以及所有不同距离的相干函数分别进行了参数拟合,再对由不同距离相干函数的拟合参数进行二次回归,比较了所有的拟合结果,得到结果如下:

    1) 不同台站间距强震记录的相干函数拟合参数的离散性很大,但根据各个不同距离的拟合参数计算得到的相干函数符合随距离的增大而减小的规律。

    2) 对于某一特定距离,由该距离的相干函数值得到的拟合结果γ与由所有不同距离相干函数值得到的拟合结果γ12存在明显的差异,距离越大或越小,偏差越大。这种偏差将会造成在采用由所有相干系数得到的拟合参数进行地震动场的人工合成时,短距离的模拟地震动相关性偏大,长距离的模拟地震动相关性偏小。

    3) 对不同距离相干函数的拟合参数进行二次回归,对于相同的距离,根据修正后的参数得到的拟合相干函数γs与该距离单独得到的拟合相干函数γ符合程度得到了明显提高,这说明本文的修正结果更适合应用于实际的所有不同场点间距(200—3000 m)的空间地震动的人工模拟。

    4) 给出了基于Loh相干函数模型的SMART-1台阵第5次和第45次地震的修正拟合参数,以及两次地震统一的修正拟合参数。

  • 期刊类型引用(0)

    其他类型引用(2)

计量
  • 文章访问数:  1484
  • HTML全文浏览量:  245
  • PDF下载量:  100
  • 被引次数: 2
出版历程
  • 发布日期:  2009-02-15

目录

/

返回文章
返回