2013年云南洱源MS5.5地震序列不同地震观测台网地震定位结果对比分析

来源 :地震研究 | 被引量 : 0次 | 上传用户:allviolet
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  摘要:基于云南地震台网和“喜玛拉雅计划”台阵探测项目流动台阵所记录到的波形资料,分3种情况对2013年云南洱源MS5.5地震进行重新定位:第一种仅利用云南地震台网数据进行定位,第二种仅利用流动台阵数据进行定位,第三种是综合云南地震台网和流动台阵数据,并结合波形互相关双差定位法进行定位。研究结果显示:在水平方向上,这3种情况得到的地震震中优势分布方向大致相同,与乔后一维西断裂NNW走向一致,但后2种情况所得到震中分布更集中,呈现的条带状更明显;在深度分布方面,后2种情况得到的震源深度比仅用云南地震台网资料进行双差定位得到的深度浅,且更加集中;在残差方面,结合波形互相关的双差定位法得到的震源位置估算误差在EW向、NS向、UD向最小。这表明台站密度更大,布局更合理,得到的定位结果误差明显降低,无论在水平方向还是深度方向分布都更为集中。
  关键词:双差定位;波形互相关;云南地震台网;流动台阵;洱源地震序列
  中图分类号:P315.6 文献标识码:A 文章编号:1000-0666(2018)01-0055-090引言
  高精度的地震定位可以更好地揭示地震序列的空間分布特征,探讨地震与断裂构造之间的关系和地震成因。Waldhauser和Ellsworth(2000)提出了一种定位精度相对较高的相对定位方法——双差地震定位法,该方法对一定空间范围内的地震进行组对,再使用地震对的走时差进行定位,这在很大程度上消除了介质横向不均匀造成的路径效应,从而获得较高精度的相对空间位置分布,因此已得到广泛应用(Waldhauser et al,2001,2002,2004;杨智娴等,2003,2004a,b;李志海等,2004;Oka-da et al,2005;王小平等,2005;黄媛等,2006;冯建刚,2008;朱艾斓等,2008;陈翰林等,2009;刘文邦等,2011;房立华等,2011)。
  2013年3月3日13时41分15秒,云南省大理州洱源县发生了MS5.5地震,之后于4月17日再次发生MS5.0地震。对于此次地震,云南地震台网记录到大量观测数据,为深入研究该区域的构造环境提供了重要条件。此外,由中国地震局地球物理研究所承担的“喜马拉雅计划”台阵探测一期项目在云南省地区布设的流动台站(以下简称流动台阵)也记录到了此次地震序列。本文基于2套台网记录的地震波形资料,首先对2013年洱源地震序列分别进行双差定位,在此基础上,综合云南地震台网和流动台阵的资料,并结合波形互相关进行重新定位,讨论3种情况下定位结果之间的差异。
  1 地震序列概况
  2013年3月3日。时52分16秒,云南省大理州洱源县(26.16°N,99.97°E)发生了ML3.6地震,随后形成序列活动,并于同日13时41分15秒以及4月17日9时45分56秒分别发生MS5.5和MS5.0地震。该地震序列共造成了云南洱源、云龙、漾濞等地60000余人受灾,20多人受伤,2500多间房屋不同程度的受损。
  截至2013年5月31日,云南地震台网共记录到该序列地震共1300多次(包括单台276次),其中ML0.0~0.9地震837次,ML1.0~1.9地震423次,ML2.0~2.9地震78次,ML3.0~3.9地震16次,ML4.0~4.9地震1次,M,,5.0~5.9地震2次(图1)。
  2 资料与速度模型选取
  2.1 资料选取
  由于受台阵观测资料的限制,本文仅选取了2013年3月3日至4月30日发生的298次ML≥1.0交切地震进行重新定位,其中,ML1.0~ML1.9地震227次,ML2.0~ML2.9地震58次,ML3.0~ML3.9地震10次,ML4.0~ML4.9地震1次,ML5.0~ML5.9地震2次。重新定位所用台站的震中距均在300km以内(图2),其中,云南地震台网固定台站32个,震中距最小的是EYA台,约为30km,流动台阵100个,震中距最小的是双漳台,约为15km。
  双差定位时,需给定P波和S波的权重值。P波震相为初至波,读数精度较高,而S波震相不易识别,加之人为因素的影响,使得其到时的读取精度比P波震相要低,所以在实际操作过程中对P波和S波赋予不同的权重值。本文选取的震相均为300km以内的近震震相,数据相对可靠,因此设定P波和S波的权重值分别为1.0和0.8。在挑选地震时,保证每个地震至少有4个震相记录,地震对之间的最大距离为10km,地震对使用的最大震相对的数目为20,事件对与台站间的最大距离参数(MAXDIST)为200km。
  2.2 速度模型
  黄媛等(2006)认为双差定位虽然可以在一定程度上消除路径效应,但仍需要可靠的速度结构作为参数。本研究所采用的分层速度模型参考了胡鸿翔等(1986)的爆破地震研究结果。该模型共分为5层,各层速度如表1所示,波速比VP/VS的值根据李永华等(2009)利用洱源台接收函数H-k方法得到的结果,设为1.75。
  3 洱源序列定位结果
  3.1 云南地震台网双差定位结果
  计算时采用共轭梯度法求解方程,并进行2轮迭代,第一轮的8次迭代由地震的初始位置和先验权重开始,第二轮的7次迭代采用标准差的5倍作为截断值,通过反复迭代,舍去残差大于截断值的震相数据,并用上一次迭代的结果更新震源位置、残差和偏导数矩阵,每次迭代得到的残差作为下一次迭代的加权函数,直到得到稳定的解。
  基于以上计算,使用双差定位法对298次地震进行精确定位后,得到了263次事件的定位结果,为原始选取地震的88%。图3为重新定位前后的地震震中分布图,从图中可以看到,该地震序列发生于维西一乔后断裂的左侧,余震的分布整体沿NNW向展布,与维西一乔后断裂的走向一致。重新定位后,外围地震向内收敛,震中分布更加集中,优势分布更为明显。   为分析地震深度分布情况,将地震震源深度沿纬度和经度方向进行投影,所得结果见图4。其中图4a、b为云南地震台网原始定位震源深度分布,图4c、d为双差定位后的震源深度分布。通过对比可以看出,双差定位后地震震源深度明显收敛。为了更加清楚,做了重新定位前后的地震震源深度分布直方图(图5)。从图5a中可以看出,重新定位前,震源深度分布较为零散,分布在2~14km的范围内,以5~12km最多,而双差定位后,深度集中分布在6~10km范围内,形成明显的优势分布层(图5b)。这与赵小艳和付虹(2014)的研究结果相同。
  重新定位后,走时残差均方根由原来的0.21s降为0.05s,震源位置的平均估算误差(2倍标准差)在IEW向为0.262km,NS向为0.294km,UD向为0.375kmo图6为EW向、NS向、UD向误差统计结果,90%以上地震误差在100~400m。
  3.2 流动台阵双差定位结果
  三维速度模型下的双差定位方法数值实验显示,台网布局、震相识别误差以及所使用的一维速度模型等因素均对双差定位结果有一定的影响;另外,双差定位方法所使用的初始地震目录对最终定位结果也有很大的影响(闫俊岗等,2013),因此用密集臺站得到的观测报告所给出的精定位结果应该比云南地震台网给出的结果可信度要高。“喜马拉雅计划”台阵探测项目在云南省地区布设的密集流动台站,为2013年3月洱源地震序列的精确定位提供了重要条件,也为深入研究该区域地震地质环境提供了可能。由于数据量大,为了方便与台网定位结果进行比较,所使用地震事件与台网精确定位所用事件相同。
  计算时采用共轭梯度法求解方程,进行2轮迭代,每轮进行8次,直到得到稳定的解。通过计算最终得到了280次事件的定位结果,为原始选取地震的93.3%。图7为重新定位前后的地震震中分布图,对比云南地震台网定位结果(图3)可以看出,两者得到的震中分布格局大致相同,这说明云南地震台网初始定位结果是较为可靠的。
  将地震震源深度沿纬度和经度方向进行投影,得到结果见图8。作重新定位前后的地震震源深度分布直方图(图9),从图8、9中可以看出,重新定位前,震源深度分散在4~18km,主要分布在5~11km,双差定位后,深度明显收敛,主要分布在5~8km,其中以6~8km最具优势。
  重新定位后,走时残差均方根由0.197。降为0.05s,震源位置的平均估算误差(2倍标准差)在EW向为0.226km,NS向为0.236km,UD向为0.24km。图10为EW向、NS向、UD向误差统计结果,可以看出90%以上地震误差分布在150~300m。
  3.3 结合波形互相关的双差定位结果
  鉴于人工拾取到时震相较多依赖于工作人员的经验,难免会出现偏差甚至误标的情况,通过波形互相关技术对P波震相到时进行校正,可以获取更加准确的到时数据,进而很大程度上解决传统方法中速度结构不够准确造成地震定位的矢量分散问题,使定位精度高达几十米甚至几米,因此在国内外都得到了广泛的运用。如在对北加州Hayward断层区域的地震定位(Waldhauser,Ellsworth,2000),对2伪8年汶川地震(黄媛,2006)等研究中,波形互相关方法都显示出其可行性和优越性。本文综合了前文中所使用的云南地震台网和流动台阵的波形数据,并利用结合波形互相关的双差定位法对序列进行重新定位,得到了270次事件的定位结果,为原始选取地震的90.6%(图11)。重新定位后的序列NNW向优势分布明显,与维西一乔后断裂走向一致,分布特征与使用流动台阵数据得到的结果(图7b)基本相同。
  将地震震源深度沿纬度和经度方向进行投影,得到震源深度分布结果如图12所示,另作震源深度分布直方图(图13)。从图13中可以看出重新定位后,地震震源深度主要分布在6~10km,其中又以7~9km最为明显。
  重新定位后,走时残差均方根由原来的0.185。降为0.043”震源位置的平均估算误差(2倍标准差)在EW向为0,118km,NS向为0.141km,UD向为0.214km。图10为EW向,NS向,UD向误差统计结果,可以看出EW向,NS向误差85%以上集中在50~150m,UD向误差主要集中在100~300m。
  4 结论和讨论
  基于云南地震台网和流动台阵的资料,分3种情况对2013年洱源MS5.5地震序列进行了双差定位,结果显示:
  (1)重定位前后,在水平方向上地震均分布在维西一乔后断裂的左侧,且优势分布方向与维西一乔后断裂的走向一致,为NNW向。但用流动台阵资料进行双差定位并综合云南地震台网和流动台阵资料,结合波形互相关方法得到的震中分布相比仅用云南地震台网得到的震中分布更为集中,优势分布方向更明显。
  (2)对于相同的地震目录,用云南地震台网资料进行双差定位得到的定位事件是263个,为原始选取地震的88%;用流动台阵资料得到的定位事件是280次,为原始选取地震的93.3%;综合云南地震台网和流动台阵资料,并结合波形互相关双差定位法得到的定位事件是270次,为原始选取地震的90.6%。
  (3)震源深度投影显示,云南地震台网定位震源深度主要分布在6~10km,流动台阵定位结果85%以上集中在6~8km,结合波形互相关的双差定位法得到的震源深度主要是7~9km,后两者的震源深度相对前者较浅,且更加集中。这说明近台记录参与联合反演更有利于震源深度的确定。
  (4)用云南地震台网资料进行双差定位,走时残差均方根由原来的0.21s降为0.05s,震源位置平均估算误差在EW向、NS向、UD向分别为0.262km、0.294km、0.375km,90%以上分布在100~400m;用台阵资料进行双差定位,走时残差均方根由原来的0.197s降为0.05s,震源位置平均估算误差在EW向、NS向、UD向分别为0.226km、0.236km、0.24km,90%以上分布在150~300m;结合波形互相关双差定位法进行定位,走时残差均方根由原来的0.185s降为0.043s,震源位置平均估算误差在EW向、NS向、UD向分别为0.118km、0.141km、0.214km,85%以上集中在50~150m,UD向主要集中在100~300m。第3种情况的定位误差相比前两种情况明显降低,说明定位效果最好。   综合以上分析认为:近台记录的增加,密集的台站分布,并结合波形互相关技术进行双差定位,得到的地震震中分布更为集中,优势方向更为明显,且测定误差无论在水平方向还是垂直方向都显著降低,这为提高该区震情的判断及掌握断层活动情况提供了基础。
  感谢中国地震局地球物理研究所“中国地震科学探测台阵数据中心”为本研究提供地震波形数据。
  参考文献:
  陈翰林,赵翠萍,修济刚,等.2009.龙滩水库地震精定位及活动特征研究[J].地球物理学报,52(8):2035-2043.
  房立华,吴建平,张天中,等.2011.2011年云南盈江MS5.8地震及其余震序列重定位[J].地震学报,33(2):262-267,
  冯建刚.2008.利用双差地震定位法对民乐—山丹地震序列重新定位[J].西北地震学报,30(1):62-65.
  胡鸿翔,陆涵行,王椿镛,等.1986:滇西地区地壳结构的爆破地震研究[J].地球物理学报,29(2):133-144.
  黄媛,杨建思,张天中.2006.2003年新疆巴楚—伽师地震序列的双差法重新定位研究[J].地球物理学报,49(1):162-169.
  李永华,吴庆举,田小波,等.2009.用接收函数方法研究云南及其邻区地壳上地幔结构[J].地球物理学报,52(1):67-80.
  李志海,赵翠萍,王海涛,等.2004.双差地震定位法在北夫山地区地震精确定位中的初步应用[J].内陆地震,18(2):146-153.
  刘文邦,王培玲,马玉虎,等.2011.2009年青海大柴旦6.4级地震序列的双差法重新定位研究[J].高原地震,23(1):24-26,
  王小平,王燕纹,李慧民.2005:结合双差地震定位法及台阵技术对江苏张家港地震序列进行精确定位[J].内陆地震,19(3):258-264,
  闫俊岗,王利兵,谭青.2013.初始震源深度对双差地震定位深度的影响分析[J].大地测量与地球动力学,33(增刊1):41-44.
  杨智娴,陈运泰,郑月军,等.2003.双差地震定位法在我国中西部地区地震精确定位中的应用[J].中国科学:地球科学,33(B4):129-134.
  杨智娴,陈运泰.2004a.用双差定位法再次精确测定1998年张北—尚义地震序列的震源参数[J].地震学报,26(2):115-119.
  杨智娴,于湘伟,郑月军,等.2004b.中国中西部地区地震的重新定位和三维地壳速度结构[J].地震学报,26(1):19-29.
  赵小艳,付虹.2014.2013年洱源MS5.5和MS5.0地震发震构造识别[J].地震学报,36(4);640-650.
  朱艾斓,徐锡伟,刁桂苓,等.2008.汶川MS8.0地震部分余震重新定位及地震构造初步分析[J].地震地質,30(3);759-768.
  Okada T,Hasegawa A,Suganomata J,et al.2005.Aftershock distributionan 3D seismic velocity structrue in and around the focal area of the2004 mid Niigata prefecture earthquake obtained by applying double-difference tomography to dense temporary seismic network data[J].Earth Planet Space,57(5):435-440.
  Waldhauser F,Ellsworth W L.2000.A double-difference earthquake lo-cation algorithm:Method and application to the Northern HaywardFault,California[J].Bull Seism Soc Am,90(6):1353-1368.
  Waldhauser F,Ellsworth W L.2002.Fault structure and mechanics of theHayward Fault,California,from double-difference earthquake loca-tion[J].J Geophys Re,3:l-15.
  Waldhauser F,Schaff D,Richards P G,et al.2004.Lop Nor revisited:Nu-cear explosion locations,1976~1996 form double-difference analy-sis of regional and teleseismic data[J].Bull Seism Soc Am,94:1879-1889a.
  Waldhauser F.2001.HypeDD-A program to compute double-differencehypocenter locations[R].US Geological Survey Open-file Report,113.
其他文献
2009年7月13日下午14-16时左右,北京最大降水量达74毫米,西三环主路被淹,丰益桥下积水水深达1米左右,造成西三环,西四环,南三环,南四环等地区发生大面积拥堵,截至晚22点左右(除丰益桥)方恢复正常。  2011年7月23日下午,北京突降暴雨,最大降水量达到173.2毫米,造成地铁1号线,4号线部分线路雨水倒灌,暂停运营,同时,西二环,西三环,西四环,丰益桥,南四环,南五环发生大面积拥堵,
期刊
摘要:为提高震后受灾人口估算精度和灾区人口空间分布可视化表达效果,以云南省大理州为研究区,结合FROM-GLC10地表覆盖数据和NPP/VIIRS夜间灯光数据,基于乡镇人口统计数据,展开人口空间化方法研究,制作了大理州50 m×50 m人口格网数据。使用相关系数、相对误差和平均相对误差对人口模拟结果进行了检验,以漾濞5.1级地震和云龙5.0级地震为例对模拟结果进行了应用分析。结果表明:①结合FRO
期刊
地震的孕育和发生与地壳形变密切相关,而地壳形变是地震过程中最直观的现象,因此精确测定区域地壳形变状态,能够为地震监测与预报提供重要的依据。我国构造动力环境特殊、孕震构造背景复杂,因此现今地壳形变观测已由传统的大地测量(如水准测量、三角测量等)发展为以GNSS、重力和InSAR技术为代表的现代大地测量和典型构造部位的区域形变测量(如跨断层观测和应变/倾斜连续观测),逐渐形成了地震大地测量学科。  近
期刊
摘要:基于2017年九寨沟MS70地震经过重定位的4 035个余震数据,利用小震确定大震断层面方法,得到该地震的断层面参数,以余震稀疏区为界将余震分为南北两段进行分段拟合,最后判定了此次地震的发震构造。结果表明:九寨沟地震的发震断层面走向1513°、倾角840°,该结果可视为发震断层参数的一级近似;南、北两段的走向和倾角分别为1474°和883°,1629°和768°,南、北两段较大的差异可能是余
期刊
摘要:采用基于人口公里格网的地震压埋人员分布评估方法,结合甘肃省人口公里格网数据,以ArcGis软件为平台,对2013年甘肃岷县漳县6.6级地震进行了压埋人员的分布评估,并结合甘肃地区的人口分布区域特点及岷县漳县6.6级地震的实际特点,对该方法中的参数进行本地化改进,并将评估结果与实际分布进行比较分析。结果表明:改进后的方法使评估结果的适用性有了一定程度的提高,可在震后快速给出震区人员压埋分布评估
期刊
摘要:2000—2001年川滇地区有大震震情趋势,应用GPS地形变测量(G),Kaiser效应法地应力测量(K)技术,采用GK联合并网方法,同步监测云南昆明—玉溪地区及邻区地形变场、地应力场;依托本构定律及岩石学组合判据,进行发震应力条件、发震几率、震级大小的分析研究;得出云南地区3~5年内不具备发生MS≥7.0大震的发震应力条件;同时划出滇东北形变高值区,10年左右有可能发生MS6.0~7.0地
期刊
摘要:基于约束条件下最小二乘原理及最速下降法,采用最小二乘配置拟合后的震间期GPS速度场作为约束,对安宁河—则木河断裂断层面滑动分布进行了反演。研究结果表明:由GPS资料计算所获取的断层面滑动分布,与地震地质研究中得到的断层面凹凸体分布具有较好的一致性,即安宁河断裂(西昌至冕宁以北)闭锁深度较深,强震危险性较大;而则木河断裂闭锁深度较浅,近期发生强震的可能性较小。  关键词:安宁河—则木河断裂;滑
期刊
摘要:选取新疆阿图什市下辖的琼哈拉峻村为研究区,将小型旋翼无人机拍摄的图像作为数据源,分别采用面向对象以及面向像元2种影像分析方法对研究区的房屋建筑进行提取,并对2种算法的提取结果进行对比,分析了各自的优势。结果表明:面向对象方法可以有效地去除椒盐噪声对分类带来的影响,保证房屋形态的完整性,但影响内部相似的光谱、纹理信息若对应多种物体则会导致影像对象的错分。在面向像元的提取方法中加入了改进的数学形
期刊
摘要:为了研究地震作用下高层框架结构的损伤程度信息,以结构刚度折减率为损伤程度指标,以结构的频率变化率为损伤程度识别参数,采用Matlab模拟结构在不同的损伤程度指标下的加速度响应数据,利用短时傅里叶变换方法对响应数据分析得到结构的模态参数,从而建立损伤程度指标与结构模态参数的函数关系。将损伤结构的模态参数代入函数关系式计算结构的损伤程度指标。采用同济大学振动台试验数据,利用此方法识别结构的损伤程
期刊
摘要:基于西秦岭构造区区域水准、长水准剖面观测资料,结合GPS速度场约束和跨断层短水准观测数据,借助三维负位错联合反演,分析2013年岷县漳县6.6级地震前相关断裂的中长期闭锁背景,讨论断裂活动的中—短临动态演变以及震后变化。结果表明:1999—2006年,岷县漳县地震前的震区所在块体即西秦岭北缘断裂与临潭—宕昌断裂之间水平面最大剪应变率、挤压性质的面应变率、垂向梯度率都较大,与构造和动力学背景一
期刊