聪明文档网

聪明文档网

最新最全的文档下载
当前位置: 首页> 安全验证

安全验证

时间:2024-01-11 14:14:04    下载该word文档
河北永清M4.3地震北京烈度仪台网记录分析
白立新;成云辉;张杰;马士振;何荣帅;冯刚
【摘要】2018212日河北永清发生M4.3地震,北京烈度仪台网共35个台站获得有效记录,记录最大峰值加速度的台站为丰台地震台(D0001,峰值加速度为10.76cm/s2,该台距震中约61km.本文对北京烈度仪台网获得的此次地震动记录进行了初步分析,并对3个同台址的强震动仪记录的数据进行了对比分析.结果表明,度仪记录计算得出的仪器烈度与强震动仪记录计算的仪器烈度结果具有较高的一致,可用于烈度速报.%Anearthquakewithmagnitude4.3occurredon12February2018inYongqingCounty,LangfangCityofHebeiProvince.Thereare35seismicstationsthatreceivedtheaccelerationdataintheBeijingseismicintensitymeternetwork.TheneareststationisFengtaistation,withanepicentraldistance61km,andpeakgroundaccelerationwas10.7cm/s2.ThefartheststationisSankuaishistation,withanepicentraldistance158km,andthepeakgroundaccelerationwas2.1cm/s2.Thepreliminaryanalysisshowsthatthedatarecordedbytheseismicintensitymeterswereclearandcomplete.Datarecordedbytheseismicintensitymeterandtheaccelerationmetersatthesamestationarecomparabletoeachother.
【期刊名称】《震灾防御技术》【年(,期】2019(014001【总页数】10(P210-219

【关键词】烈度仪;强震动仪;噪声;仪器烈度【作者】白立新;成云辉;张杰;马士振;何荣帅;冯刚
【作者单位】北京市地震局,北京100080;北京市地震局,北京100080;北京市地震,北京100080;北京市地震局,北京100080;北京市地震局,北京100080;北京市地震局,北京100080【正文语种】引言国内外重大地震灾害救援行动的经验和教训表明,震后政府的快速有效应急响应是破坏性地震发生后保持社会稳定的重要因素。破坏性地震发生后,应急管理部门、公用事业和其他相关管理机构最迫切的需求是了解地震造成破坏的相关信息,并依此做出快速的应急响应和科学的救灾决策,最大程度减少人员灾害损失,尽可能在较短的时间内恢复社会秩序。传统的强震动台网是获取地震破坏分布信息的手段之一,但是,大规模、高密度强震动台网存在着诸如建设成本高昂和运维难度大等困难。Evans等(2005)提出了基于低成本微机电系统(MEMS)烈度仪建设地震台网的设想,用以弥补传统强震动台网由于台站间距大、获取的地震灾害分布信息分辨率低的缺点。近年来,研究人员对烈度仪与传统强震动仪记录的数据进行了大量对比研究,探讨了烈度仪用于地震监测的可行性和可靠性。观测实验结果表明,烈度仪具备中强以上地震观测的能力(D’Alessandro等,2013;王浩等,2013Evans等,2014)。由于烈度仪价格低廉、布设方便、维护成本低,一些国家和地区近年来相继开展了高密度烈度仪观测台网建设,以期通过建设高密度烈度仪台网提升获取地震烈度分布信息的时效性和可靠性(Cochran等,2011Hsieh等,
2014Clayton等,2015;于海英等,2017Wu等,2018)。
2015年北京市地震局开展了“简易烈度仪地震预警示范项目”建设,在北京建设40个烈度仪台站。2018212日河北永清发生M4.3地震,北京市行政区范围内普遍有感。地震发生后,北京市地震局立即对烈度仪台网获得的地震动记录数据进行了处理,其中,有35个烈度仪台站获得了有效的震动记录。本文将对北京烈度仪台网获得的此次地震记录进行了初步分析。1烈度仪台站烈度仪台站均建设在土层场地上,布设地点主要为北京地震台网内已经停测或无法实时传输的强震动台站,烈度仪采用地面安装方式。利用无线网络进行数据实时传输,数据采样率为100Hz,并通过JOPENS软件系统实现数据汇集。2记录情况2018212日河北永清M4.3地震中,北京烈度仪台网35个台站共获得105条记录,另外5个台站由于电源故障或设备死机,未获取记录。在此次地震中,北京烈度仪台网记录峰值加速度最大的台站是丰台地震台(D0001),峰值加速度10.76cm/s2,方向为东西向,该台站距震中约61km。距震中最远的烈度仪台站为怀柔三块石台(M0004),距震中约158km,峰值加速度为2.13cm/s21为丰台地震台和三块石台烈度仪记录的地震波形。此次地震35个台站记录的初至震相均清晰可辨认,显示烈度仪台站具备良好的获取地震记录的能力。1烈度仪台加速度记录Fig.1Accelarationrecordfromtheseismicintensitymeterstation3仪器烈度计算根据《仪器地震烈度计算暂行规程》(中国地震局监测预报司,2015),首先对获取的记录应进行基线校正处理,采用4阶巴特沃斯(Butterworth)带通滤波器对每个分向记录进行0.110Hz带通滤波,通过积分加速度记录获得速度时程。
仪器烈度计算公式如下:

其中,PGA为合成三方向加速度峰值,PGV为三方向合成速度峰值。如果IPGAIPGV均大于等于6.0,则仪器地震烈度IIPGV值,否则IIPGAIPGV的算术平均值。
通过计算,产出了本次地震的仪器烈度分布图,结果显示烈度仪台站仪器烈度最大的为丰台地震台,仪器烈度为Ⅲ度,分布于平原区的烈度仪台站仪器烈度多为Ⅱ度,北部山区仪器烈度为Ⅰ度(图2),显示出仪器烈度随震中距增加而逐渐减小的特征。
2烈度仪台站仪器烈度分布Fig.2Thedistributionofintensitiesfromintensitymeterstations4与强震动台站记录对比分析为验证烈度仪台站获取地震动记录的能力以及记录结果的可靠性,在南口中学台11NKZ)、大华山台(11DHS)和三块石台(11SKS)安装了烈度仪,烈度仪台编号为K0001(南口中学台)、N0004(大华山台)和M0004(三块石台)。3个台站的强震动仪数据采集器型号均为Reftek130-ren,加速度传感器型号为SLJ-100;烈度仪为北京港震仪器设备有限公司生产的GL-P2B型烈度仪。烈度仪和强震动仪主要性能指标见表1。在此次地震中3个台站的强震动仪和烈度仪均获得了记录,11NKZ11DHS11SKS台站的震中距分别为108km106km158km
1烈度仪和强震动仪主要参数指标Table1Themainspecificationsofseismicintensitymetersandaccelerationmeters技术指标GL-P2BSLJ-100类型MEMS传感器力平衡式加速度计传感器频带范围DC-80HzDC-80Hz态范围大于等于95dB≥120dB满量程水平向2.0g,垂直向-13g2.0gADC
2824数据采集采样率100Hz100Hz数字滤波器最小相位FIR线性相位FIR3为强震动仪与烈度仪的原始记录波形,数据记录按同一起始时刻截取。由图3可见,强震动仪与烈度仪记录的地震波形及峰值加速度十分相近,相应通道记录的峰值最大偏差均小于0.1cm/s2。为定量评价强震动仪与烈度仪获取的地震记录在时域和频域上的相关程度,采用式(3)计算强震动仪与烈度仪记录的时域皮尔森Pearson)相关系数,用式(4)计算强震动仪与烈度仪记录的幅度平方相干系数。

其中,rxy为相关系数,反映的是两过程在时域上的线性相关程度,相关系数越接近于1,相关性越强,convxy)为xy的协方差,σx、σy分别为xy标准偏差;Cxy为幅度平方相干系数,反映的是两过程在各频率上分量间的线性相关程度,Cxy1时,表示2个时程信号完全相干,Pxyxy的互功率谱密度,PxxPyy分别为xy的自功率谱密度。
通过计算得到南口中学台(11NKZK0001EWNSUD分向的皮尔森相关系数分别为0.980.980.97;大华山台(11DHSN0004EWNSUD分向的相关系数分别为0.960.970.92;三块石台(11SKSM0004EWNSUD分向的相关系数分别为0.970.960.89。计算结果表明强震动仪与烈度仪的地震记录在时域上是强相关的。
43个台站的强震动仪与烈度仪记录在频域上的相干性结果,其中南口中学台和大华山台强震动仪与烈度仪的3个分向的波形记录在120Hz频率范围内具有强相关性,而三块石台强震动仪与烈度仪的3个分向的波形记录在220Hz率范围内具有强相关性。三块石台站(11SKSM0004)的相干频带小于其它2个台站,可能与震中距较大(158km)、地震记录信号衰减过大有关。

3同台址烈度仪与强震动仪原始记录波形Fig.3Uncorrectedrecordofseismicintensitymeterandaccelerationmeteratsamestation4同台址强震动仪与烈度仪原始记录幅度平方相干性结果Fig.4Magnitudesquaredcoherencebetweenseismicintensitymeterandaccelerationmeteratsamestation5对比了同台址强震动仪与烈度仪记录的3个分向数据的信噪比,结果显示,0.120Hz范围内,3个烈度仪台站事件(用实线表示,下同)与噪声(用虚线表示,下同)的信噪比远低于同台址的强震动仪记录事件与噪声的信噪比,尤其是低频信号记录能力明显比同台址的强震动仪低。
依据《仪器地震烈度计算暂行规程》(中国地震局监测预报司,2015)对数据记录进行0.120Hz带通滤波,计算了3个台站的仪器烈度。11NKZK0001台站仪器烈度计算结果为Ⅱ度,11DHSN000411SKSM0004台站仪器烈度为Ⅰ度,表明强震动仪与烈度仪的仪器烈度具有较高的一致性。
为了解烈度仪获取加速度记录的可积分性,通过积分计算3个台站的强震动仪和烈度仪记录的速度和位移时程,并采用4阶巴特沃斯带通滤波器对数据记录进行滤波。依据幅度平方相干性结果,11NKZK000111DHSN0004台站滤波频带取120Hz11SKSM0004台站滤波频带取220Hz。滤波后,强震动仪与烈度仪的峰值加速度均有所降低,表明地震记录存在一定程度的高频成分。采用式(3)计算了3个台站的的强震动仪与烈度仪加速度记录积分后获得的速度和位移的相关系数,其中11NKZK0001台站EWNSUD分向的速度记录相关系数为0.990.990.98;位移记录相关系数为0.990.990.9911DHSN0004台站EWNSUD分向的速度记录相关系数为0.990.960.97;位移记录相关系数为0.960.970.9711SKSM0004台站EWNSUD向的速度记录相关系数为0.980.980.96;位移记录相关系数为0.970.97
0.97
5同台址未校正烈度仪与强震仪记录的信号及噪声谱Fig.5Singnalandnoisespectraofintensitymeterandaccelerationmeteratsamestation计算结果表明,强震动仪与烈度仪的地震记录积分后获得的速度和位移在时域上是强相关的,烈度仪加速度时程具有一定的可积分性。图6给出了3个台站的垂直向校正加速度记录以及积分得到的速度和位移时程。
6烈度仪与强震动仪校正记录波形Fig.6Correctedrecordofseismicintensitymeterandaccelerationmeteratsamestation5结论本文对2018212日河北永清M4.3地震的北京烈度仪台网记录情况进行了讨论,对同台址强震动仪和烈度仪的记录波形进行了对比,结果表明:
1)距震中158km的台站清晰地记录了此次地震,北京烈度仪台网具备记录中强地震的能力。
2)同台址强震动仪和烈度仪地震记录在时域上是强相关的,南口中学台11NKZK0001)和大华山台(11DHSN0004)强震动仪与烈度仪3个分向的波形记录在120Hz频率范围内具有强相关性,而三块石台(11SKSM0004强震动仪与烈度仪3个分向的波形记录在220Hz频率范围内具有强相关性。三块石台(11SKSM0004)的相干频带小于其它2个台站,可能是因其震中距较大,导致地震信号衰减过大。
3)强震动仪的信噪比远高于同台址的烈度仪,烈度仪记录低频信号能力明显低于同台址强震动仪。
4)烈度仪产出的仪器烈度结果与强震动仪一致。5)烈度仪加速度时程具有一定的可积分性。
综上所述,在烈度仪设备稳定运行的前提下,烈度仪台可用于进行烈度速报,通过
建设价格相对低廉的烈度仪台网提高烈度速报时效性是一条可行的途径。在建设烈度仪台网时,需考虑烈度仪的动态范围指标,避免烈度仪台站功能单一化。致谢:本文采用了北京市强震动台网提供的地震数据,中国地震局工程力学研究所国家强震动台网中心提供的强震数据处理软件,并利用matplotlib开源软件进行图件绘制;审稿专家提出了宝贵的修改建议,在此一并表示感谢。参考文献


【相关文献】
王浩,丁炜,2013.MEMS加速度计与传统地震加速度计的比较研究.大地测量与地球动力学,33S2):9395.于海英,周宝峰,王家行等,2017.强震动观测仪器面临的机遇和挑战.震灾防御技术,121):6877.中国地震局监测预报司,2015.仪器地震烈度计算暂行规程.北京:地震出版社,110.ClaytonR.W.,HeatonT.,KohlerM.,etal.,2015.Communityseismicnetwork:adensearraytosenseearthquakestrongmotion.SeismologicalResearchLetters,86(5:13541363.CochranE.S.,LawrenceJ.F.,KaiserA.,etal.,2011.Comparisonbetweenlow-costandtraditionalMEMSaccelerometers:acasestudyfromtheM7.1Darfield,NewZealand,aftershockdeployment.AnnalsofGeophysics,54(6:728737.D’AlessandroA.,D’AnnaG.,2013.Suitabilityoflow-costthree-axisMEMSaccelerometersinstrong-motionseismology:testsontheLIS331DLH(iPhoneaccelerometer.BulletinoftheSeismologicalSocietyofAmerica,103(5:29062913.EvansJ.R.Hamstra,R.H.,KündigC.,etal.,2005.TREMOR:awirelessMEMSaccelerographfordensearrays.EarthquakeSpectra,21(1:91124.EvansJ.R.,AllenR.M.,ChungA.I.,etal.,2014.Performanceofseverallow-costaccelerometers.SeismologicalResearchLetters,85(1:147158.HsiehC.Y.,WuY.M.,ChinT.L.,etal.,2014.LowcostseismicnetworkpracticalapplicationsforproducingquickshakingmapsinTaiwan.TerrestrialAtmosphericandOceanicSciences,25(5:617624.WuY.M.,MittalH.,HuangT.C.,etal.,2018.PerformanceofaLow-CostEarthquakeEarly
WarningSystem(P-AlertandShakeMapProductionduringthe2018MW6.4Hualien,Taiwan,Earthquake.SeismologicalResearchLetters,90(1:1929.

免费下载 Word文档免费下载: 安全验证

  • 29.8

    ¥45 每天只需1.0元
    1个月 推荐
  • 9.9

    ¥15
    1天
  • 59.8

    ¥90
    3个月

选择支付方式

  • 微信付款
郑重提醒:支付后,系统自动为您完成注册

请使用微信扫码支付(元)

订单号:
支付后,系统自动为您完成注册
遇到问题请联系 在线客服

常用手机号:
用于找回密码
图片验证码:
看不清?点击更换
短信验证码:
新密码:
 
绑定后可用手机号登录
请不要关闭本页面,支付完成后请点击【支付完成】按钮
遇到问题请联系 在线客服