结合波形互相关的双差定位方法在三峡水库地震中的应用1

罗佳宏 马文涛 李春政



结合波形互相关的双差定位方法在三峡水库地震中的应用1

罗佳宏 马文涛 李春政

（中国地震局地质研究所活动构造与火山重点实验室，北京 100029）

本文采用基于波形互相关算法的双差定位方法对三峡水库地震进行精定位和地震活动性分析。首先使用双谱互相关方法分析了三峡库区加密台网于2009年3月至2010年12月观测到的地震波形数据，并对波形互相关分析的结果进行了评价。基于结合获得的波形互相关数据使用双差定位方法对地震事件进行精定位研究，结果表明使用双谱法验证的波形互相关数据的定位精度要高于其他数据的结果，其东西向震源位置平均误差为3.2m、南北向为3.9m、垂直向为6.2m。重定位震中结果显示巴东神龙溪两岸微震分布明显呈现出3条近东西向的线性条带状，与地表小规模断裂和碳酸盐岩地层走向一致，揭示了库水主要沿着溶洞或者地下暗河渗透进而诱发地震活动，较强地震可能是微小地震贯穿活动面的结果。

双谱验证 波形互相关 双差定位 三峡水库

引言

地震定位是地震学中最基本、最经典的问题之一，其核心内容是最大限度地提高地震定位精度。影响地震定位精度的因素很多，比如记录台站的分布、震相拾取的精度、速度模型的准确性等，其中影响最大的是速度模型误差和震相到时读取误差。基于波形互相关技术的双差定位方法能够很好地克服这2种误差（黄媛，2008），提高了数字地震波形应用范围。使用互相关技术对P波和S波走时进行校正，减少到时读数误差的影响，很大程度上降低由于直接读取震相到时所造成的误差，并解决因速度模型不准确造成的定位结果分散问题，提高了地震定位的精度。

波形互相关分析分为时域和频域2种计算方法（Poupinet等，1984）。时域互相关分析操作简单，获得了较广泛的应用（Schaff等，2004；2005；Waldhauser等，2000），针对垂直分量互相关计算的局限性，发展了时域多通道相关检测函数并用于计算波形互相关走时差（王清东等，2015）。在时域中选取互相关系数的阈值需要有一定的经验，如果阈值设定过高，可用数据就会很少，达不到高精度定位的要求，但如果阈值取太低，则会得到一些非真实震相信息。一个比较稳妥的方法便是使用双谱法（Du等，2004）对互相关系数结果进行验证，可以更有效地增强互相关的可用信息并提高互相关数据的可信度。双谱法能有效地压制高斯互相关噪声，在三阶谱域内同时计算原始和滤波后波形的时间延时，并用这2个时间延时来验证互相关系数的可靠性。因其具有精度高的特点，自提出以来便得到了一定的应用（Bannister等，2011；Bourguignon等，2015；Hansen等，2013；张广伟等，2015；赵翠萍，2006）。

本文首先利用双谱法软件包（BCSEIS）对三峡库区地震事件的波形数据进行互相关分析，并用双差定位方法对地震事件进行精定位，对不同类型数据的定位结果进行比较，获得了三峡库区高精度微震定位结果。

1 地震波形数据处理

三峡水库地震监测台网由26个地震台站组成，包括21个英国L-22型三分向短周期速度摆和5个美国Guralp CMG-3ESPC宽频地震计，数采均为Reftek130且采样率为200Hz（马文涛等，2010），在2009年3月至2010年12月的观测时间范围内一共获得了5275次地震事件观测报告及其波形数据。

在本研究中，我们首先对地震波形数据进行1.0—10Hz带通滤波，得到原始和滤波后波形数据，然后分别在P波和S波到时前后截取一段波形进行互相关分析。我们取P波时窗长128个采样点，P波到时前30个采样点，之后97个采样点；取S波时窗长192个采样点，S波到时前50个采样点，之后141个采样点。图1显示了1个事件对的波形记录，2个地震事件的编号分别为242和244，图中只绘制了发震时刻到S波时窗范围内的波形记录，此外的波形不参与计算波形互相关系数，事件242和244的详细信息见表1。波形互相关结算结果说明事件对之间具有高度相似，2次地震事件的波形互相关系数大于0.9，震相延时小于0.04s。将地震目录震相的到时差和波形互相关延时差进行比较（图2），波形互相关延时可以降低由于震相拾取误差带来的偏差，能获得更高精度的到时差数据，说明用双谱验证方法提取的波形互相关系数具有较高的精度。

在双谱计算过程中我们使用3个互相关系数阈值来检验计算结果：中央范围系数、下限和上限。一般情况下，相当于传统互相关延时计算中的阈值，在实际检测时我们使用以下原则：

表1 图1中地震事件对详细信息

在双差定位或双差层析成像方法中，地震对之间的最大距离对反演结果有很大的影响。如果该值取得越小，则建立联系的地震对之间距离越小，定位精度越高，但是太小则能够建立联系的地震就越少（黄媛等，2006）。在进行双差定位或双差层析成像之前，需要评价地震对之间的距离。一般可以计算地震对之间的距离并统计深度分布直方图，从而可以快速地确定合适的地震距，我们计算了研究范围内地震事件之间的距离（图4），表明地震对之间的距离取20km较为合理，这个值与程序ph2dt给出的“强连接”事件对的最大距离也是一致的。但是该地震距是否与互相关地震对之间的距离匹配，需要进行详细的分析。根据双谱验证法获得的波形互相关数据，统计所有台站记录到的事件对震源距和相应的波形互相关系数，得到了互相关系数与事件对震源距分布图（图5）。结果表明，单纯地设置地震对之间距离阈值为20km会降低互相关数据的使用质量，定位结果受横向不均匀性的影响较大，因此需要适当地减小互相关地震对之间距离的阈值，这里我们选取6km，以保证大多数互相关系数和距离的一致性，并确保互相关数据的使用效率和定位的精度。

2 三峡库区双差定位实现

根据2009年3月至2010年12月中国地震局地质研究所三峡水库加密台网观测数据结果，一共得到5275个地震事件，根据地震目录获得P波走时数据43538条和S波走时数据43385条，平均每个地震事件有7个清晰台站记录。为排除震相判读错误和其他信息的干扰，使用最小二乘拟合走时曲线，剔除误差较大的数据，同时使用和达曲线（即S-P曲线）对结果进行检验（蔺永等，2014）。图6表明数据拟合前存在一系列的干扰数据，对数据进行走时曲线最小二乘拟合，剔除误差较大的震相后，得到误差较小的地震事件信息。经过挑选后，得到5244个地震事件，共41423条P波和42354条S波震相数据。图7为5244个地震事件的信息分布图，分别给出了沿着经度、纬度和深度的地震事件分布图。结果显示，地震事件主要分布在8km深度以上的范围，且地震分布相对集中，长江水系外围分布有离散地震事件。

根据三峡库区已有的速度模型研究结果（李强等，2009；赵旭等，2007），选定速度模型的水平深度分别取0、2、5、8、11、14和20km，对应的P波速度为4.8、5.4、5.65、5.8、6、6.15和6.5km/s，P/S值根据图6的拟合结果取为1.73。

3 三峡水库地震精定位结果及分析

根据BCSEIS波形互相关分析的结果和地震目录数据，我们使用hypoDD双差定位软件对三峡库区水库地震事件进行3组测试：第一组，仅利用地震目录数据进行定位（CAT）；第二组，利用标准的波形互相关方法获得的波形延时和地震目录到时差数据（CC+CAT）进行定位；第三组，利用双谱验证获得的波形延时数据和地震目录到时差数据（WCC+CAT）进行定位。为了对比定位结果的差异性，对3组测试数据设置相同的参数和速度模型，最终获得不同数据类型的地震平面位置（图8）和深度误差分布图（图9）。图8表明，第二组和第三组的定位结果明显优于第一组数据的结果，具有较小的定位误差和优势深度（图9）。

重定位后的地震震中分布更加集中，呈线性化分布，显示出地震事件向内收敛的趋势。在巴东神龙溪两岸表现出3条明显的近东西向条状分布，震源深度较浅，平均深度在5km左右，震中分布与地层走向基本一致，进一步验证了水库蓄水后库水从神龙溪等地下暗河渗入而诱发地震的解释（马文涛等，2010）。中间一条地震分布线东口于2013年12月16日发生了S5.1级地震，从2009年3月—2010年12月的微震定位结果可以说明，该地震可能是在库水的作用下，贯穿整个东西向的小规模断裂引发的较大地震活动。在泄滩西区域存在的一个陡立东西向的微震活动带，其可能和长江水系的渗透有关。而在秭归香溪河口附近，地震分布在仙女山断裂带上，震源深度向下延深到10km，揭示了库水渗透的最大深度。该处于2014年3月发生2次4级以上地震，可能也是库水渗透引起断裂局部活动的结果。

3组不同类型数据的定位误差统计结果（图10）表明使用双谱验证获得波形互相关延时与地震目录到时差数据的hypoDD定位精度明显要高于其他2种数据的精度，地震震源位置的平均误差在东西向为3.2m、南北向为3.9m、垂直向为6.2m。

4 结论

本文应用双谱分析方法，对三峡库区2009年3月至2010年12月的地震观测数据进行了波形互相关分析，获得了高质量的双谱验证波形互相关数据，并使用双差定位方法对地震进行了重定位，获得了高精度的定位结果。

基于三阶谱域的双谱验证波形互相关算法能获得比传统阈值判定更高质量的波形延时数据，得到的互相关地震事件对的距离更加符合实际情况，能为双差定位或者双差层析成像提供更高质量的数据资料。

根据双谱验证方法获得的传统阈值判定和双谱验证的波形互相关数据，分别使用双差定位程序进行了定位，结果表明结合双谱验证的波形互相关数据的定位结果精度高于其他数据的结果，其东西向震源位置平均误差为3.2m、南北向为3.9m、垂直向为6.2m。

对三峡水库地震重定位结果表明，巴东库区神龙溪两岸微震分布明显呈现出3条近东西向的线性分布特征，与新构造时期产生的小规模断裂和碳酸盐岩地层走向一致，揭示了库水主要沿着溶洞或地下暗河渗透而诱发地震活动，其中中间1条地震分布线端点正是发生2013年12月16日湖北巴东5.1级地震的震中位置，说明它是在库水作用下沿着东西向小规模断裂的一次较大岩体错动事件；泄滩西区东西向垂直条带状分布预示地震活动与长江库水渗透之间的关系；而在秭归香溪河口附近，地震分布在仙女山断裂带上，震源深度向下延深到10km，揭示了库水渗透的最大深度。

致谢：评审专家对本文的修改提出了宝贵的意见，文中大部分图件使用GMT绘图软件包绘制（Wessel等，1995），在此一并表示感谢。

黄媛，2008．结合波形互相关技术的双差算法在地震定位中的应用探讨．国际地震动态，（4）：29—34．

黄媛，杨建思，张天中，2006．2003年新疆巴楚—伽师地震序列的双差法重新定位研究．地球物理学报，49（1）：162—169．

李强，赵旭，蔡晋安等，2009．三峡水库坝址及邻区中上地壳P波速度结构．中国科学 D辑：地球科学，39（4）：427—436．

蔺永，马文涛，2014．基于Matlab的三峡水库地震数据处理与分析．震灾防御技术，9（3）：447—453．

马文涛，徐长朋，李海鸥等，2010．长江三峡水库诱发地震加密观测及地震成因初步分析．地震地质，32（4）：552—563．

王清东，朱良保，苏有锦等，2015．2012年9月7日彝良地震及余震序列双差定位研究．地球物理学报，58（9）：3205—3221．

张广伟，雷建设，2015．2011年云南腾冲5.2级双震发震机理．地球物理学报，58（4）：1194—1204．

赵翠萍，2006．1997—2003年新疆伽师震源区特征的地震学方法研究．北京：中国地震局地球物理研究所．

赵旭，李强，蔡晋安，2007．三峡库首区最小一维速度模型研究．大地测量与地球动力学，27（专刊）：1—7．

Bannister S., Fry B., Reyners M., et al., 2011. Fine-scale Relocation of Aftershocks of the 22 FebruaryW6.2 Christchurch Earthquake using Double-difference Tomography. Seismological Research Letters, 82(6): 839—845.

Bourguignon S., Bannister S., Henderson C. M., et al., 2015. Structural heterogeneity of the midcrust adjacent to the central Alpine Fault, New Zealand: Inferences from seismic tomography and seismicity between Harihari and Ross. Geochemistry, Geophysics, Geosystems, 16(4): 1017—1043.

Du W. X., Thurber C. H., Eberhart-Phillips D., 2004. Earthquake relocation using cross-correlation time delay estimates verified with the bispectrum method. Bulletin of the Seismological Society of America, 94(3): 856—866.

Hansen S. E., DeShon H. R., Moore-Driskell M. M., et al., 2013. Investigating the P wave velocity structure beneath Harrat Lunayyir, northwestern Saudi Arabia, using double-difference tomography and earthquakes from the 2009 seismic swarm. Journal of Geophysical Research: Solid Earth, 118(9): 4814—4826.

Poupinet G., Ellsworth W., Frechet J., 1984. Monitoring velocity variations in the crust using earthquake doublets: An application to the Calaveras Fault, California. Journal of Geophysical Research, 89(B7): 5719—5731.

Schaff D. P., Bokelmann G. H. R., Ellsworth W. L., et al., 2004. Optimizing correlation techniques for improved earthquake location. Bulletin of the Seismological Society of America, 94(2): 705—721.

Schaff D. P., Waldhauser F., 2005. Waveform cross-correlation-based differential travel-time measurements at the Northern California Seismic Network. Bulletin of the Seismological Society of America, 95(6): 2446—2461.

Waldhauser F., Ellsworth W. L., 2000. A double-difference earthquake location algorithm: Method and application to the northern Hayward fault, California. Bulletin of the Seismological Society of America, 90(6): 1353—1368.

Wessel P., Smith W. H. F., 1995. New version of the generic mapping tools. Eos, Transactions American Geophysical Union, 76(33): 329.

Application of Double-difference Relocation Method Combined with Waveforms Cross-correlation on Earthquakes in the Three Gorges Reservoir Area

Luo Jiahong, Ma Wentao and Li ChunZheng

(Key Laboratory of Active Tectonics and Volcano, Institute of Geology, China Earthquake Administration, Beijing 100029, China)

In this paper, we applied the double difference location method based on waveform cross-correlation algorithm for earthquake precision positioning of the Three Gorges Reservoir (TGR). First we used bispectrum cross-correlation method to analyze the seismic waveform data of TGR encrypted networks in March 2009 to December 2010, and evaluate the quality of waveform cross-correlation analysis. Combining the waveform cross-correlation of data obtained, we used the double difference method to relocate the position of quakes. The results show that the location by using bispectrum verified waveform cross-correlation data is higher than other type of data, and the mean 2sig-errors in EW, NS and UD are 3.2m, 3.9m and 6.2m, respectively. The results also show that the Badong and Shenlong River quake in reservoir distribution is characterized by linear distribution of three nearly east-west, which is in accordance with the small faults and carbonate strata line of new tectonic period, revealing reservoir water main along the underground river or cave penetration induced seismic activity. A strong earthquake may be the result of a small earthquakes that broken through the active plane.

Bispectrum verification；Waveform cross-correlation；Double-difference location；The Three Gorges Reservoir

1 基金项目 国家科技支撑计划（2008BAC38B04）和中国地震局专项业务项目（16A44ZX282）共同资助。

2016-06-08

罗佳宏，男，生于1989年。硕士研究生。主要从事地震层析成像、近场地震学、诱发地震等研究。E-mail：ljh1771@163.com

马文涛，男，生于1958年。副研究员。主要从事理论地震学和水库诱发地震研究。E-mail：wentaoma_1@126.com

罗佳宏，马文涛，李春政，2017．结合波形互相关的双差定位方法在三峡水库地震中的应用．震灾防御技术，12（1）：56—67． doi：10.11899/zzfy20170106