2. 中国地质科学院 地质力学研究所,北京 100081
2. Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing 100081, China
2008年5月12日,中国四川省汶川县发生了8.0级特大地震,造成了巨大的人员伤亡和经济损失,震感波及全国[1]。地震导致大量房屋倒塌,并诱发了强烈的山崩、滑坡、塌方和泥石流等次生地质灾害,造成近8万人死亡,37万多人受伤,4500多万人流离失所,导致上万亿元的财产损失。我国地球科学工作者心系祖国和人民,为减少地质灾害开展了全方位的调查研究工作[2~5]。大地震后的快速钻探是研究地震机制和捕捉余震直接信息的有效方法之一,已为世界地学界所公认[6]。2008年8月,在许志琴院士的倡导下,国家科技支撑“汶川断裂带科学钻探”项目在龙门山地震断裂带上迅速展开,从而揭开了我国对地震快速反应、地震发震机理的研究序幕。国内外的一些研究表明,通过对现今小震震源位置群体特征和震源机制的研究,可得出历史大震震源断裂的位置及其构造应力场作用下的运动方式,且单次中小地震的发生具有随机性,而监测并统计它们的群体特征,可以提取共有的信息。大量地震的震源位置和震源机制解既可以描述震源断层的空间取向和应力场的状态,也能够用来探讨大范围的构造应力场的作用,因此,利用数字化地震台网的高精度资料去研究已经发生大震的震源构造是一条可行的途径[7]。
1 台网建设意义地震发生后,余震分布区沿北东向龙门山活动断裂带分布,长度近300 km[8]。有限元模拟计算和实际余震发展表明,由震中区向北东方向延伸直到青川、宁强、陇南一带为应力增强区,5级以上强余震分布较多,应加强监视; 震中区西南方向主要为应力下降区,小部分为增强区[9]。另外,余震自四川大邑向陕西宁强呈北东向发展(见图 1),是否沿青川断裂向汉中方向发展,是否还会发生强烈余震成为社会关注的问题。为了查明汶川8.0级地震后青川断层库仑应力的变化,判明余震的发展趋势[10],中国地质科学院地质力学研究所实施安装了“龙门山断裂带WFSD数字地震台网”,开展龙门山断裂带地震监测工作,为监测龙门山断裂带余震分布、迁移趋势以及判断断层活动程度提供基础资料,进一步为搞清发震机理和地震复发周期提供科学依据。
龙门山地震破裂带位于人口密集的四川省成都盆地西缘,是四川盆地的西北界; 北西与松潘—甘孜褶皱系为邻,南东与扬子地台—四川盆地接壤,北东为秦岭构造带所截,南西与康滇地轴相接,为我国西部地槽区与东部地台区接壤地带,是北起广元,南达天全,长约500 km,宽为50~80 km的狭长地带。龙门山中段是指界于北川和都江堰之间的部分,以后山带出露彭灌杂岩和前山带发育飞来峰群为显著特征[11]。龙门山推覆构造的主体由3条近于平行的北东—南西向的区域性大断裂构成,分别为汶川—茂县断裂(后山断裂)、映秀—北川断裂(中央断裂)和灌县—江油断裂(前山断裂)。
2008年5.12汶川特大地震发生以及龙门山向南东方向推覆的动力来源是印度板块与欧亚大陆碰撞及其向北的推挤的结果(见图 2)[10]。这一板块间相对运动导致了亚洲大陆内部大规模的构造变形,造成了青藏高原地体向东挤出、地貌隆升。但龙门山断裂带的地壳水平位移速率极低[2, 12],而与刚性的阿拉善地块、鄂尔多斯地块和扬子地块接壤的整个青藏高原北缘和东缘,GPS测量均显示水平位移速率为低值(小于3 mm / a),而汶川8.0级特大地震就发生在青藏高原的东缘———龙门山断裂带上。汶川地震究竟在什么样的应力环境下产生?其发震成因机制如何?这些问题需要通过诸如数字微地震台网地震观测、科学钻探、地应力监测等等多种渠道去探索。
龙门山断裂带数字地震台网由1个数据处理中心和11个数字地震台站组成(见表 1),台网东西约110 km,南北约120 km。地震台站主要分布在龙门山的东北缘和西南端(见图 3)。数字地震台站由地震计、数据采集器、通信设备和供电系统组成。地震发生后,大安台站首先于2008年6月建成并传回地震数据(见图 4)。
① 地震计是将地面运动直接转换成电压型的设备。在龙门山断裂带地震台网中除了纸厂台站(ZC)采用三分向一体安装的宽频带地震计(频带范围在120 s~60 Hz)外,其他10个台站均采用了短周期三分量地震计(频带范围在2 s~50 Hz)。
② 数据采集器是将地震计输出的模拟电压量转换为数字化电压量的装置。龙门山断裂带数字地震台网均采用北京港震公司生产的EDAS-24型低功耗地震数据采集器,它具有24位A / D转换器、分辨率高、动态范围大的特点,其动态范围大于130 dB (采样率为50Hz时)。
③ 通信设备:地震台站至台网中心的八角庙台站(BJM)采用DDN专线传输,传输速率9600 bps; 其它各台站均采用CDMA无线传输,带宽25 KC,传输速率为4800 bps。
④ 供电系统:由于大部分台站放置在野外,附近没有接交流电的条件,所以对各台站的供电采用了市电、太阳能+蓄电池相结合2种供电方式。
⑤ 台网中心系统:具有对多路地震波形数据收集、地震台站数据采集器监管、地震事件检测、地震波形数据记录及地震事件处理等功能。
综合上述,整个数字地震观测网络采用了先进的数字观测技术,具有宽频带、动态范围大、分辨率高等特点,可获得高精度的地震波形数据。从运行的实践看,对龙门山断裂带附近地区的微小地震以及<1000 km的ML>4.0的中远震,地震波记录清晰。
4 地震台站台基背景噪声分析WFSD台网是要真实地记录龙门山地区地震引起的地面运动过程,为项目研究震后余震的发生发展趋势、震后断层活动、震后地面运动提供充分和完善的基础数据。其任务是监测震区乃至整个龙门山地区的微震及小震,进行定位和震源机制分析,以确定震后区内所发生地震的空间分布和震源的力学机制。台网对区内地震的观测和监控能力对项目研究十分重要,需要对台网的监测能力进行分析。因此,台网部署完成后,对龙门山断裂带上所有台站进行了台基背景噪声分析,以此为基础计算各个台站观测动态范围和地震监控能力,并进一步结合台网实际记录的地震事件进行了检验。
影响台网地震监控能力及观测动态范围的主要因素有:地震台观测系统观测灵敏度、台基背景噪声、仪器动态范围等[13]。选择2009年4月的台站记录,在日常观测数据中随机选取10天,将每天的观测数据分为6个观测时段,从每个观测时段中抽取60 s连续记录。在选取数据时尽量避开了地震、标定和其他偶然干扰(如车辆、爆破、工业生产振动和电磁等人为干扰以及起风、下雨、气压骤变等天气干扰),共选取360000点数据(相当于1个小时的数据),组合为分析样本。
功率谱密度计算采用经典谱估计的周期图法,它把随机信号f (x)的N 点观测数据视为能量有限信号,直接取f (x)的傅立叶变换,得F (u) :
(1) |
然后再取其幅值的平方作为对f (x)真实的功率谱P (u)的估计[14]。
(2) |
台站台基地动噪声有效值RMS与功率谱有如下关系:
(3) |
式中: PSD———功率谱密度; fc———分度倍频程中心频率; RBW———相对带宽。
相对带宽RBW由下式计算:
(4) |
式中: fu———分度倍频程上限频率; fd———分度倍频程下限频率。
台站分析样本功率谱的具体计算步骤如下:
① 去直流偏移;
② 将Counts转换为速度;
③ 采用窗函数抑制频谱泄漏;
④ 傅立叶变换;
⑤ 扣除仪器响应,最后得到功率谱密度函数。
计算过程利用北京港震公司提供的Noise_ psd 011噪声功率谱密度测定程序对样本进行分析,得到各台站1 /3倍频程(1~20 Hz)频带内的台基地平均噪声功率谱密度值(见图 5)和地动噪声的RMS值(UD向) (见表 2)。龙门山活动断裂带数字地震台网所有台站的台基背景噪声均在10 -8 m / s数量级上,最优的达到了1.548 × 10 -8 m / s。
台站观测的动态范围反映了台站记录大地震事件完整波形的能力。数字地震记录受观测仪器、传输方式和场地脉动噪声的影响,而模拟地震台站信号在传输过程中易受到干扰,观测动态范围一般在45 dB左右[13]。数字台站采用数字化传输,主要是线路带宽对采样率的限制和影响,在带宽允许的条件下,考虑到数字地震台站根据本地区地动噪声不同,如使用24位数据采集器,观测动态范围可达到90~130 dB。
目前常用的地震数据采集器为24位,量程为± 2 23 (± 8388607),本身的动态范围达到138.4 dB。考虑到台站台基地动噪声的影响,台站观测的有效动态范围为:
(5) |
式中: R———EDAS-24IP地震数据采集器量程,2 23 counts; RMS———台站台基地动噪声有效值,m / s; S———系统灵敏度,counts /μm / s; P———数采的量程修正值(1.33);
本文中8个典型观测台站的动态范围全部接近或超过100 dB,最优的达到了114 dB,完全符合规范要求(见表 3)。由于各台站台基背景噪声水平中等,在取得较好观测动态范围的情况下,台网仍保持了较高的系统观测灵敏度。
台网的监测能力取决于各地震台的布置位置、密度、台址背景噪声、所用仪器的频带和仪器响应灵敏度设置的高低。选择优良的台基,尽量安装在基岩和洞中,避开地表杂波和气象因素及人类活动的干扰,合理地搭配仪器频带范围和响应灵敏度,才能有效地发挥数字地震仪的动态范围,做到既能较清晰记录微小地震,又能使完整记录较大地震的能力得到延伸[15]。通过以上各台的背景噪声值与动态范围的分析,本文抽取样本计算了台网监测地震的下限范围,其上限范围需要根据震中距的增大、台网记录不限幅地震的能力而定。
6.1 样本计算台网地震监控能力通常是指对微小地震的监测能力,主要取决于地震计的响应灵敏度、观测动态范围和台基地动噪声平均水平[13]。
由于地动噪声的存在,地震仪分辨小于噪声的信号是困难的,尽管目前已有各种从噪声中提取信号的技术,但由于噪声的随机性和不确定性,目前地震研究一般不用小于噪声的信号。因此,安装现代高灵敏地震仪的台站,台网对微小地震的监测能力主要取决与台基地动噪声水平。
本文仍然使用上述计算台基地动噪声水平时的样本,选取东西和南北向数据滤去直流分量后,计算地动噪声幅值的有效值; 然后将2个水平向的地动噪声有效值进行平均,再扩大6倍,作为能识别S震相的阀值,此时P震相也能较为清晰地分辨,这样就使得对地震记录有足够充分的分辨率。
根据S波最大振幅确定近震震级的计算公式,用量规函数反推单台控制距离(震中距Δ)。取4个以上台站监测某一震级地区的相关区域,确定台网的监控能力。具体计算方法为:
根据短周期地震仪近震震级公式:
(6) |
将
(7) |
式中: VNS———南北向S波最大速度振幅值,μm; TNS———南北向S波最大速度振幅的周期值; VEW———东西向S波最大速度振幅值,μm; TEW———东西向最大速度振幅的周期值; R (Δ)———量规函数,其物理意义是地震波随距离的衰减; S———台站改正值,一般情况下取0。
通过上述公式计算出台网中每个子台可检测某一震级地震的最大距离(见表 4)。以各台站为圆心,以各台站检测到同一震级的最大距离为半径画圆,有4个或4个以上台站的圆弧包围的区域就是台网对某一震级的地震能够控制的监测范围。将不同震级的监测范围叠加起来,就反映出整个台网的监控能力。依据上述理论,对龙门山断裂带数字WFSD地震监测台网的监控范围进行了计算和划分,结果表明在整个龙门山断裂带上,具有监测ML 1.5级地震的能力(见图 6)。
但应强调的是,此方法在计算中的设定条件和实际地动噪声水平间有一定差异,地震监测还受各种环境和人为干扰因素影响,加上地震计的动态影响,所以获得结果只能粗略地反映台网对地震的控制范围,与实际监测可能有一定差异。
6.2 监测结果检验龙门山地震断裂带数字台网在2008年6月至2009年8月计14个月的运行期间,记录了数千个ML 1.5~ML 6.0级地震,其中2009年3月25日8: 00汶川台记录到草坡ML 1.5级地震波形(见图 7)。
本文对龙门山地震断裂带数字台网地震记录波形及定位结果与国家地震局台网数据中心定位结果进行了对比,结果见表 5。从表 5中可见:两个台网地震三要素定位基本吻合,所以龙门山地震断裂带数字台网能够满足地震活动监测和地震预报的需要。
龙门山地震断裂带数字台网所记录的地震波形完整,震相容易分辩,如:四川省平武县ML4.2级地震波形(见图 8)。
WFSD数字地震台网建成后,通过后期实践检验,具有布局合理、宽频带、动态范围大、精度高等特点。该数字地震台网将会提高WFSD项目辅助观测地震技术水平,对加快特大地震机理及发震机制研究有着重要作用。
7 结论(1) 通过对龙门山断裂带地震观测样本分析,获得各台站1 /3倍频程(1~20 Hz)频带内的台基地平均噪声功率谱密度值和地动噪声的RMS值(UD向)。表明各台站的台基背景噪声均在10 -8 m / s数量级上,最优的达到了1.548 × 10 -8 m / s。
(2) 所有观测台站的动态范围全部接近或超过100 dB,最优的达到了114 dB,完全符合规范要求。由于各台站台基背景噪声水平中等,在取得较好观测动态范围的情况下,台网保持了较高的系统观测灵敏度。
(3) 在整个龙门山断裂带上,龙门山断裂带WFSD数字地震监测台网具有监测ML 1.5级地震的能力。
(4) 利用WFSD数字台网地震数据记录进行定位和定级,其结果与中国地震局台网中心发布的结果基本吻合。
致谢 WFSD数字地震台网建设,得到四川、陕西两省有关市、县、镇领导的支持及相关机构同仁的帮助,在此表示衷心感谢。
[1] |
彭华, 马秀敏, 姜景捷. 山丹地应力监测站体应变仪的地震效应[J]. 地质力学学报, 2008, 14(2): 97-108. PENG Hua, MA Xiumin, JIANG Jingjie. Analysis of the volume strain data from the Shandan in-situ stress Monitoring Station[J]. Journal of Geomechanics, 2008, 14(2): 97-108. (in Chinese) |
[2] |
张永双, 石菊松, 孙萍, 等. 汶川地震内外动力耦合及灾害实例[J]. 地质力学学报, 2009, 15(2): 131-141. ZHANG Yong-shuang, SHI Ju-shong, SUN Ping, et al. Coupling between endogenic and exogenic geological processes in the Wenchuan earthquake and example analysis of geo-hazards[J]. Journal of Geomechanics, 2009, 15(2): 131-141. (in Chinese) |
[3] |
谭成轩, 孙叶, 吴树仁, 等. 5 12汶川MS 8.0大地震后关于我国区域地壳稳定性评价的思考[J]. 地质力学学报, 2009, 15(2): 142-150. TAN Cheng-xuan, SUN Ye, WU Shu-ren, et al. A consideration on regional crustal stability assessment after MS 8.0 Wenchuan strong earthquake in China[J]. Journal of Geomechanics, 2009, 15(2): 142-150. (in Chinese) |
[4] |
王连捷, 崔军文, 周春景, 等. 汶川5.12地震发震机理的数值模拟[J]. 地质力学学报, 2009, 15(2): 105-113. WANG Lian-jie, CUI Jun-wen, ZHOU Chun-jing, et al. Numerical modeling for Wenchuan earthquake mechanism[J]. Journal of Geomechanics, 2009, 15(2): 105-113. (in Chinese) |
[5] |
姚鑫, 张永双. 基于差分干涉雷达的汶川地震同震形变特点[J]. 地质力学学报, 2009, 15(2): 151-161. YAO Xin, ZHANG Yong-shuang. Co-seismic Deformation of "5.12 " Wenchuan earthquake based on D-INSAR[J]. Journal of Geomechanics, 2009, 15(2): 151-161. (in Chinese) |
[6] |
许志琴, 李海兵, 吴忠良. 汶川地震和科学钻探[J]. 地质学报, 2008, 82(12): 1613-1622. XU Zhi-qin, LI Hai-bing, WU Zhong-liang. Wenchuan earthquake and scientific drilling[J]. Acta Geologica Sinica, 2008, 82(12): 1613-1622. DOI:10.3321/j.issn:0001-5717.2008.12.001 (in Chinese) |
[7] |
张宏志, 刁桂苓, 陈祺福. 1976年唐山7.8级地震震区现今地震震源机制分析[J]. 地震研究, 2008, 31(1): 1-6. ZHANG Hong-zhi, DIAO Gui-ling, CHEN Qi-fu. Focal mechanism analysis of the recent earthquakes in Tangshan seismic region of M 7.6 in 1976[J]. Journal of Seismological Research, 2008, 31(1): 1-6. (in Chinese) |
[8] |
张永双, 雷伟志, 石菊松, 等. 四川5.12地震次生地质灾害的基本特征初析[J]. 地质力学学报, 2008, 14(2): 110-116. ZHANG Yong-shuang, LEI Wei-zhi, SHI Ju-song, et al. General characteristics of 5.12 earthquake induced geo-hazards in Sichuan[J]. Journal of Geomechanics, 2008, 14(2): 110-116. (in Chinese) |
[9] |
王连捷, 周春景, 孙东生, 等. 汶川5.12地震引起的库仑应力变化及其对周边地震活动的影响[J]. 地质力学学报, 2008, 14(3): 193-200. WANG Lian-jie, ZHOU Chun-jing, SUN Dong-sheng. Coulomb stress changes caused by Wenchuan earthquake and its influence on seismic activity in the adjacent area[J]. Journal of Geomechanics, 2008, 14(3): 193-200. (in Chinese) |
[10] |
彭华, 马秀敏, 姜景捷. 龙门山北端青川断层附近应力测量与断层稳定性[J]. 地质力学学报, 2009, 15(2): 114-130. PENG Hua, MA Xiu-min, JIANG Jing-jie. Stability and stress measurement near the Qingchuan fault in the northern Longmen Mountains[J]. Journal of Geomechanics, 2009, 15(2): 114-130. (in Chinese) |
[11] |
刘树根, 田小彬, 李智武. 龙门山中段构造特征与汶川地震[J]. 成都理工大学学报:自然科学版, 2008, 35(4): 388-397. LIU Shu-gen, TIAN Xiao-bin, LI Zhi-wu. Structural features of the central Longmen Mountains and the Wenchuan earthquake in Sichuan, China[J]. Journal of Chengdu University of Technology: Science & Technology Edition, 2008, 35(4): 388-397. (in Chinese) |
[12] |
Zhang P Z, Shen Z, Wang M, et al. Continuous deformation of the Tibetan Plateau from Global Positioning System data[J]. Geology, 2004, 32: 809-812. DOI:10.1130/G20554.1 |
[13] |
李涌, 张艺, 戴仕贵. 阿坝州数字地震台网的动态范围与台网监控能力计算[J]. 四川地震, 2007, 4: 37-41. LI Yong, ZHANG Yi, DAI Shi-gui. Dynamic range and monitoring ability computation results of the Digital Seismic Network of Aba Prefecture[J]. Earthquake Research in Sichuan, 2007, 4: 37-41. (in Chinese) |
[14] |
邵玉平, 王翠芳, 宋澄. 地震台网勘选数据的自动化处理[J]. 四川地震, 2008, 1: 9-12. SHAO Yu-ping, WANG Cui-fang, SONG Cheng. The automation of site survey data processing for seismic network[J]. Earthquake Research in Sichuan, 2008, 1: 9-12. DOI:10.3969/j.issn.1001-8115.2008.01.003 (in Chinese) |
[15] |
唐明帅, 段天山. "十五"新疆数字地震台网建设[J]. 地震地磁观测与研究, 2007, 28(6): 58-62. TANG Ming-shuai, DUAN Tian-shan. An introduction of " the Tenth Five-year Plan " digital earthquake network construction in Xinjiang[J]. Seismological and Geomagnetic Observation and Research, 2007, 28(6): 58-62. (in Chinese) |