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

0 引言

高精度的地震定位可以更好地揭示地震序列的空间分布特征，探讨地震与断裂构造之间的关系和地震成因。Waldhauser和Ellsworth(2000)提出了一种定位精度相对较高的相对定位方法——双差地震定位法，该方法对一定空间范围内的地震进行组对，再使用地震对的走时差进行定位，这在很大程度上消除了介质横向不均匀造成的路径效应，从而获得较高精度的相对空间位置分布，因此已得到广泛应用(Waldhauser et al，2001，2002，2004；杨智娴等，2003，2004a，b；李志海等，2004；Okada 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日0时52分16秒，云南省大理州洱源县(26.16°N，99.97°E)发生了ML3.6地震，随后形成序列活动，并于同日13时41分15秒以及4月17日9时45分56秒分别发生MS5.5和MS5.0地震。该地震序列共造成了云南洱源、云龙、漾濞等地60 000余人受灾，20多人受伤，2 500多间房屋不同程度的受损。

截至2013年5月31日，云南地震台网共记录到该序列地震共1 300多次(包括单台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次，ML5.0～5.9地震2次(图1)。

图1 2013年洱源地震序列ML≥0地震

M-t图(a)，N-t图(b)

Fig.1 M-t diagram(a)and N-t diagram(b)

for Eryuan ML≥0 earthquake sequence in 2013

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次。重新定位所用台站的震中距均在300 km以内(图2)，其中，云南地震台网固定台站32个，震中距最小的是EYA台，约为30 km，流动台阵100个，震中距最小的是双涧台，约为15 km。

(1)重定位前后，在水平方向上地震均分布在维西―乔后断裂的左侧，且优势分布方向与维西―乔后断裂的走向一致，为NNW向。但用流动台阵资料进行双差定位并综合云南地震台网和流动台阵资料，结合波形互相关方法得到的震中分布相比仅用云南地震台网得到的震中分布更为集中，优势分布方向更明显。

如今，智慧宫与国内多家出版社合作翻译阿语图书，一年翻译数百本书，内容也由文学作品扩展到社会科学等方面。随着业务越做越大，主动上门来谈合作的境外文化机构也越来越多。白鑫和他的中国合伙人专门从阿拉伯国家邀请了50多名顶尖汉学家翻译中国图书，尽可能地保持著作的原汁原味。

图2 用于重新定位的震中距300 km以内地震台站分布

Fig.2 Distribution of seismic stations with epicenter distance of ≤300 km for relocation

2.2 速度模型

重新定位后，走时残差均方根由原来的0.21 s降为0.05 s，震源位置的平均估算误差(2倍标准差)在EW向为0.262 km，NS向为0.294 km，UD向为0.375 km。图6为EW向、NS向、UD向误差统计结果，90%以上地震误差在100～400 m。

表1 双差定位采用的速度模型 Tab.1 The velocity model used in double-difference location

层序号深度/kmP波速度/(km·s-1)10.284.30213.615.93333.506.30442.906.60580.007.75

3 洱源序列定位结果

3.1 云南地震台网双差定位结果

计算时采用共轭梯度法求解方程，并进行2轮迭代，第一轮的8次迭代由地震的初始位置和先验权重开始，第二轮的7次迭代采用标准差的5倍作为截断值，通过反复迭代，舍去残差大于截断值的震相数据，并用上一次迭代的结果更新震源位置、残差和偏导数矩阵，每次迭代得到的残差作为下一次迭代的加权函数，直到得到稳定的解。

基于以上计算，使用双差定位法对298次地震进行精确定位后，得到了263次事件的定位结果，为原始选取地震的88%。图3为重新定位前后的地震震中分布图，从图中可以看到，该地震序列发生于维西—乔后断裂的左侧，余震的分布整体沿NNW向展布，与维西—乔后断裂的走向一致。重新定位后，外围地震向内收敛，震中分布更加集中，优势分布更为明显。

计算时采用共轭梯度法求解方程，进行2轮迭代，每轮进行8次，直到得到稳定的解。通过计算最终得到了280次事件的定位结果，为原始选取地震的93.3%。图7为重新定位前后的地震震中分布图，对比云南地震台网定位结果(图3)可以看出，两者得到的震中分布格局大致相同，这说明云南地震台网初始定位结果是较为可靠的。

图3 云南地震台网资料双差定位前(a)、后(b)震中分布 Fig.3 Epicenter distribution before(a)and after(b)double difference relocation by using data from Yunnan Seismic Network

图4 云南地震台网资料双差定位前(a、b)、后(c、d)震源深度分别沿经纬度方向的投影 Fig.4 The focal depth distribution along the latitude and longitude direction before(a，b)and after(c，d)double difference relocation by using data from Yunnan Seismic Network

为分析地震深度分布情况，将地震震源深度沿纬度和经度方向进行投影，所得结果见图4。其中图4a、b为云南地震台网原始定位震源深度分布，图4c、d为双差定位后的震源深度分布。通过对比可以看出，双差定位后地震震源深度明显收敛。为了更加清楚，做了重新定位前后的地震震源深度分布直方图(图5)。从图5a中可以看出，重新定位前，震源深度分布较为零散，分布在2～14 km的范围内，以5～12 km最多，而双差定位后，深度集中分布在6～10 km范围内，形成明显的优势分布层(图5b)。这与赵小艳和付虹(2014)的研究结果相同。

黄媛等(2006)认为双差定位虽然可以在一定程度上消除路径效应，但仍需要可靠的速度结构作为参数。本研究所采用的分层速度模型参考了胡鸿翔等(1986)的爆破地震研究结果。该模型共分为5层，各层速度如表1所示，波速比VP/VS的值根据李永华等(2009)利用洱源台接收函数H-k方法得到的结果，设为1.75。

基于云南地震台网和流动台阵的资料，分3种情况对2013年洱源MS5.5地震序列进行了双差定位，结果显示：

3.2 流动台阵双差定位结果

三维速度模型下的双差定位方法数值实验显示，台网布局、震相识别误差以及所使用的一维速度模型等因素均对双差定位结果有一定的影响；另外，双差定位方法所使用的初始地震目录对最终定位结果也有很大的影响(闫俊岗等，2013)，因此用密集台站得到的观测报告所给出的精定位结果应该比云南地震台网给出的结果可信度要高。“喜马拉雅计划”台阵探测项目在云南省地区布设的密集流动台站，为2013年3月洱源地震序列的精确定位提供了重要条件，也为深入研究该区域地震地质环境提供了可能。由于数据量大，为了方便与台网定位结果进行比较，所使用地震事件与台网精确定位所用事件相同。

图5 云南地震台网资料双差定位前(a)、后(b)震源深度分布柱状图 Fig.5 The histogram of focal depth before(a)and after(b)double-difference relocation by using data from Yunnan Seismic Network

图6 云南地震台网资料的双差定位结果误差统计 Fig.6 Relocation error statistics of double-difference location result by using the data from Yunnan Seismic Network

我这次来山里，是想劝父亲跟我回去。两年前，我因车祸，右胳膊落下残疾，这对以绘画为生的我来说，简直是致命的打击。幸好，在最痛苦的时候，是父亲一直鼓励着我，安慰着我，我才渐渐走出伤痛，生活又恢复平静。

图7 流动台阵资料双差定位前(a)、后(b)震中分布 Fig.7 Epicenter distribution before(a)and after(b)double-difference relocation by using data from flow temporary array

图8 流动台阵资料双差定位前(a、b)、后(c、d)深度分别沿经纬度方向的投影 Fig.8 The focal depth distribution along the latitude and longitude directions before(a，b) and after(c，d)double-difference relocation by using data from flow temporary array

将地震震源深度沿纬度和经度方向进行投影，得到结果见图8。作重新定位前后的地震震源深度分布直方图(图9)，从图8、9中可以看出，重新定位前，震源深度分散在4～18 km，主要分布在5～11 km，双差定位后，深度明显收敛，主要分布在5～8 km，其中以6～8 km最具优势。

重新定位后，走时残差均方根由0.197 s降为0.05 s，震源位置的平均估算误差(2倍标准差)在EW向为0.226 km，NS向为0.236 km，UD向为0.24 km。图10为EW向、NS向、UD向误差统计结果，可以看出90%以上地震误差分布在150～300 m。

式中：ηt为汽轮机理想循环热效率；ηi为汽轮机的相对内效率；ηm为汽轮机的机械效率，取值0.99；ηg为发电机效率，取值0.99；ηb为锅炉热效率；ηp为管道效率，取值0.99。

图9 流动台阵资料双差定位前(a)、后(b)震源深度分布柱状图 Fig.9 The histogram of focal depth before(a)and after(b)double-difference relocation by using data of flow temporary array

图10 流动台阵资料的双差定位结果误差统计 Fig.10 Relocation error statistics of double-difference location result by using the data of flow temporary array

3.3 结合波形互相关的双差定位结果

鉴于人工拾取到时震相较多依赖于工作人员的经验，难免会出现偏差甚至误标的情况，通过波形互相关技术对P波震相到时进行校正，可以获取更加准确的到时数据，进而很大程度上解决传统方法中速度结构不够准确造成地震定位的矢量分散问题，使定位精度高达几十米甚至几米，因此在国内外都得到了广泛的运用。如在对北加州Hayward断层区域的地震定位(Waldhauser，Ellsworth，2000)，对2008年汶川地震(黄媛，2006)等研究中，波形互相关方法都显示出其可行性和优越性。本文综合了前文中所使用的云南地震台网和流动台阵的波形数据，并利用结合波形互相关的双差定位法对序列进行重新定位，得到了270次事件的定位结果，为原始选取地震的90.6%(图11)。重新定位后的序列NNW向优势分布明显，与维西—乔后断裂走向一致，分布特征与使用流动台阵数据得到的结果(图7b)基本相同。

如图8所示，此综合绝缘抱杆装置需两个稳固抱箍，分别由连接块、锁紧带构成，用于将抱杆固定在电杆上。连接块由金属铸造，机械强度足够；锁紧带由布质材料制成，布带式锁紧带较原有金属链条抱箍相比，既保证机械强度同时又具有良好的绝缘效果，也减轻了配件重量，方便安装。

将地震震源深度沿纬度和经度方向进行投影，得到震源深度分布结果如图12所示，另作震源深度分布直方图(图13)。从图13中可以看出重新定位后，地震震源深度主要分布在6～10 km，其中又以7～9 km最为明显。

重新定位后，走时残差均方根由原来的0.185 s降为0.043 s，震源位置的平均估算误差(2倍标准差)在EW向为0.118 km，NS向为0.141 km，UD向为0.214 km。图10为EW向，NS向，UD向误差统计结果，可以看出EW向，NS向误差85%以上集中在50～150 m，UD向误差主要集中在100～300 m。

图12 结合波形互相关双差定位的震源深度沿纬度(a)和经度(b)方向的投影 Fig.12 The focal depth projection along altitude(a)and longtitude(b) by doubles difference method combined with waveform cross-correlation technique

图11 结合波形互相关双差定位的震中分布图 Fig.11 Epicenter distribution by using doubles difference method combined with waveform cross-correlation technique

图13 结合波形互相关双差定位 震源深度分布图 Fig.13 The focal depth distribution by using doubles difference method combined with waveform cross-correlation technique

图14 结合波形互相关的双差定位结果误差统计

Fig.14 Relocation error statistics by using doubles difference method combined with waveform cross-correlation technique

4 结论和讨论

第二，马克思主义是革命的理论。列宁指出：“马克思认为他的理论的全部价值在于这个理论‘按其本质来说，它是批判的和革命的’。后一性质的确完全地和无条件地是马克思主义所固有的，因为这个理论公开认为自己的任务就是揭露现代社会的一切对抗和剥削形式，考察它们的演变，证明它们的暂时性和转变为另一种形式的必然性，因而也就帮助无产阶级尽可能迅速地、尽可能容易地消灭任何剥削。”[1]82-83以往的一切所谓“科学”都在为人剥削人的制度进行辩护，只有马克思主义对这种人剥削人的制度进行了无情的批判，毫不掩饰自己的阶级性，公开宣布它就是要代表无产阶级利益，为无产阶级的革命运动服务。

双差定位时，需给定P波和S波的权重值。P波震相为初至波，读数精度较高，而S波震相不易识别，加之人为因素的影响，使得其到时的读取精度比P波震相要低，所以在实际操作过程中对P波和S波赋予不同的权重值。本文选取的震相均为300 km以内的近震震相，数据相对可靠，因此设定P波和S波的权重值分别为1.0和0.8。在挑选地震时，保证每个地震至少有4个震相记录，地震对之间的最大距离为10 km，地震对使用的最大震相对的数目为20，事件对与台站间的最大距离参数(MAXDIST)为200 km。

(2)对于相同的地震目录，用云南地震台网资料进行双差定位得到的定位事件是263个，为原始选取地震的88%；用流动台阵资料得到的定位事件是280次，为原始选取地震的93.3%；综合云南地震台网和流动台阵资料，并结合波形互相关双差定位法得到的定位事件是270次，为原始选取地震的90.6%。

但令人遗憾的是，时至今日福建很多地方对武术精神文化或武术人文价值的宣传仅依靠单一的传统媒体。譬如：调研发现在福鼎市和南安市的乡镇以下武术宣传仍以人际传播为主，传播方式依然是千百年来民间自然延续的口口相传，基本没有专业人士介入引导和书面整理；市(县)一级则采用报刊、广播、电视等传统媒体进行宣传，缺少“互联网+”的新媒体参与。长期以来，在一度“唯GDP论英雄”的时代背景下，当地政府对武术的社会经济价值及其对经济发展的直接贡献已经有了明确认知；但更深层次的认知，即对武术的文化内涵与潜在价值的理解仍显轻率，其外在表现是重视力度不够，不能创新和扩大宣传平台(大型国际性赛事除外)。

(3)震源深度投影显示，云南地震台网定位震源深度主要分布在6～10 km，流动台阵定位结果85%以上集中在6～8 km，结合波形互相关的双差定位法得到的震源深度主要是7～9 km，后两者的震源深度相对前者较浅，且更加集中。这说明近台记录参与联合反演更有利于震源深度的确定。

(4)用云南地震台网资料进行双差定位，走时残差均方根由原来的0.21 s降为0.05 s，震源位置平均估算误差在EW向、NS向、UD向分别为0.262 km、0.294 km、0.375 km，90%以上分布在100～400 m；用台阵资料进行双差定位，走时残差均方根由原来的0.197 s降为0.05 s，震源位置平均估算误差在EW向、NS向、UD向分别为0.226 km、0.236 km、0.24 km，90%以上分布在150～300 m；结合波形互相关双差定位法进行定位，走时残差均方根由原来的0.185 s降为0.043 s，震源位置平均估算误差在EW向、NS向、UD向分别为0.118 km、0.141 km、0.214 km，85%以上集中在50～150 m，UD向主要集中在100～300 m。第3种情况的定位误差相比前两种情况明显降低，说明定位效果最好。

综合以上分析认为：近台记录的增加，密集的台站分布，并结合波形互相关技术进行双差定位，得到的地震震中分布更为集中，优势方向更为明显，且测定误差无论在水平方向还是垂直方向都显著降低，这为提高该区震情的判断及掌握断层活动情况提供了基础。

GaAs晶体折射率由传统模型确定.将式(13-14)式代入式(12)中，即可得到不同调制因素变化条件下晶体的最佳取向，即折射率变化最大取向.图2-4分别给出的调频、调幅、调相条件下，GaAs晶体的折射率变化最大取向曲线.

感谢中国地震局地球物理研究所“中国地震科学探测台阵数据中心”为本研究提供地震波形数据。

闫俊岗，王利兵，谭青.2013.初始震源深度对双差地震定位深度的影响分析[J].大地测量与地球动力学，33(增刊1)：41-44.

陈翰林，赵翠萍，修济刚，等.2009.龙滩水库地震精定位及活动特征研究[J].地球物理学报，52(8)：2035-2043.

3) 相机采用CMOS工业相机,该相机支持UVC协议，所以在编译嵌入式linux内核时需添加对UVC相机的支持。具体步骤：make menuconfig ARCH=arm CROSS_COMPILE=arm-linux- / Device Drivers / Multimedia support / Video capture adapters / V4L USB devices / USB Video Class(UVC)、UVC input events devices support。

房立华，吴建平，张天中，等.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.

参考文献：

杨智娴，陈运泰，郑月军，等.2003.双差地震定位法在我国中西部地区地震精确定位中的应用[J].中国科学：地球科学，33(B4)：129-134.

杨智娴，陈运泰.2004a.用双差定位法再次精确测定1998年张北—尚义地震序列的震源参数[J].地震学报，26(2)：115-119.

杨智娴，于湘伟，郑月军，等.2004b.中国中西部地区地震的重新定位和三维地壳速度结构[J].地震学报，26(1)：19-29.

两组一般临床资料比较显示，两组在性别上比较，差异无统计学意义（P＞0.05）。冠心病组中吸烟、高血压、高血脂、糖尿病患者数量明显高于对照组，组间比较，差异有统计学意义（P＜0.05）。表明以上因子均为冠心病的危险因素。见表1。

赵小艳，付虹.2014.2013年洱源MS5.5和MS5.0地震发震构造识别[J].地震学报，36(4)：640-650.

招平断裂是胶西北三大控矿断裂之一，前人对其进行了大量勘查和研究工作，取得了一系列重大找矿成果。2010—2013年，随着深部找矿战略的进一步实施，胶西北地区的勘查深度从-1000m推进至-2000m，三大断裂中焦家断裂纱岭矿区主矿体以及三山岛断裂西岭村矿区主矿体相继突破-2000m，这2条断裂提交的金资源总量相继超过1200t，成为千吨级控矿断裂[1-2]，而黄金产量最高的招平断裂在深部找矿成果上却相对滞后(表1)。

朱艾斓，徐锡伟，刁桂苓，等.2008.汶川MS8.0地震部分余震重新定位及地震构造初步分析[J].地震地质，30(3)：759-768.

Okada T，Hasegawa A，Suganomata J，et al.2005.Aftershock distribution an 3D seismic velocity structrue in and around the focal area of the 2004 mid Niigata prefecture earthquake obtained by applying double-difference tomography to dense temporary seismic network data[J].Earth Planet Space，57(5)：435-440.

定向医学生学习动机的得分为(4.373±0.043)，接近中等临界值。在各个维度上的由高到低依次为物质追求(4.693±0.057)、求知进取(4.672±0.056)、社会取向(4.354±0.068)、小群体取向(4.337±0.055)、个人成就(4.058±0.067)和害怕失败(4.058±0.067)。在各维度得分高低的排列顺序上与刘洋的相关研究一致[4]。(见表1)

Waldhauser F，Ellsworth W L.2000.A double-difference earthquake location algorithm ：Method and application to the Northern Hayward Fault，California[J].Bull Seism Soc Am，90(6)：1353-1368.

Waldhauser F，Ellsworth W L.2002.Fault structure and mechanics of the Hayward Fault，California，from double-difference earthquake location[J].J Geophys Re，3：1-15.

走进新时代 面临新形势 矿产资源储量管理工作要有新举措（吕晓澜 ） ........................................................1-25

Waldhauser F，Schaff D，Richards P G，et al.2004.Lop Nor revisited：Nucear explosion locations，1976～1996 form double-difference analysis of regional and teleseismic data[J].Bull Seism Soc Am，94：1879-1889a.

Waldhauser F.2001.HypeDD-A program to compute double-difference hypocenter locations[R].US Geological Survey Open-file Report，113.

在以天然气为燃料的锅炉中，烟气中水蒸气的含量直接与燃料的原始组成以及锅炉的排烟温度密切相关[8-9]。当燃料中H元素的含量大时，产生的水蒸气量也会较多，此时，高、低位发热量差值越大。其中天然气等有机气态燃料为最大，液态燃料依次减小，其汽化潜热也逐次减少。