-
内蒙古锡林郭勒区域科学台阵的勘选工作来源于“一带一路”地震监测台网项目—内蒙古子项目科学台阵建设项目(以下简称“内蒙古科学台阵项目”),该项目旨在建立综合探测监测网的地震监测体系,使我国地震监测能力和地震灾害安全保障水平得以大幅提升。锡林郭勒区域科学台阵的站点共62个,目前已完成该区所有站点的实地踏勘和场地仪器测试工作。
地震台阵起初用于监测核爆实验,由于其具有特定的台站布局和数据处理方式,表现出突出的监测能力,后来被人们广泛应用于地震监测领域[1-2]。该技术不仅可以大大提升地震监测能力,同时能够帮助人们探寻更精细的地球内部结构[3]。相比地震台网,地震台阵可以汇集多个相关子台的数据,通过专用的处理技术,达到抑制地噪声,加强有效地震信号的效果,从而拥有更高的地震监测能力,尤其是对微小地震信号的识别[4-5]。
地震台阵的建设涉及多个重要环节,其中,台阵的勘测选址对台阵建设尤为重要。为建设对微小地震拥有较高监测能力的台阵,在台阵踏勘选址过程中,应尽可能选取构造块体较为完整稳定、地质结构较为均一、地脉动背景噪声较小的区域[3]。由于地震信号是指波场中的有效成分,而噪声则属于会干扰地震信息需要压制的部分,因此若台址拟选区地脉动背景噪声较大,则势必会降低台阵的地震监测能力。所以,开展台阵台址勘选背景噪声研究不仅可评价台阵的勘选质量,也可为优化台阵后期建设提供参考依据与基础资料。本文以内蒙古锡林郭勒区域为研究区,利用该区域科学台阵台址勘选的地脉动噪声数据,开展62个子台背景噪声的计算分析工作。
-
内蒙古科学台阵整体分布于内蒙古东北部区域,涉及锡林郭勒盟、赤峰市、通辽市、兴安盟、呼伦贝尔市在内的5个行政区域,共包括201个子台站点(图1)。内蒙古科学台阵的建设,可大大提高区域地震监测台网的地震监控能力,有利于建立由综合探测监测网、数据传输系统、地震观测数据处理与共享平台和运行管理保障系统组成的分级部署、业务协调的地震监测体系,进而提升中国对东北亚、南亚方向陆地以及黄渤海区域的地震监测能力。
锡林郭勒区域的科学台阵主要分布于该区43°N以北,共包括62个子台站点。台阵布局呈规则网格状,各子台之间的横纵间隔为0.5°(图1)。根据地质研究,整个锡林郭勒盟地跨东北块体和蒙古弧形构造带两大构造单元。其中,42°N以北地区属天山蒙古地槽褶皱带,具有较为均一与完整的构造块体[6]。而且锡林郭勒地区人口稀少,经济相对落后,工业欠发展,主要以畜牧业为主,因此该区子台环境较安静,干扰因素相对较少。
-
内蒙古科学台阵项目的野外勘选工作于2019年4月16日从锡林郭勒区域开始,勘选工作人员严格按照GB/T 19531.1—2004“地震台站观测环境技术要求 第1部分:测震”的相关技术规范,开展各子台的勘探选址工作,具体勘选内容由4部分组成。①点位测试:在拟选点附近寻找环境相对安静并且存在基岩出露的地方进行测试, 适当留取2~3个候选点位;②技术分析:根据测震台站建设要求,分析台站的噪声水平与可能的干扰源,确保子台噪声水平与距干扰源的最小距离均满足技术规范[7];③点位确定:在保证建台技术要求的基础上,权衡基础建设、供水供电、通信、交通、日后运维等条件,最终挑选出各方面均优秀的点位;④观测分析:对选择点位连续观测48 h,分析处理观测数据, 若分析结果满足台站观测环境的技术要求,则认为该点为台站建设点,完成点位的勘选工作。
在观测地动噪声时选用工作频带包含0.05~20 Hz的宽频带数字地震仪[8],实际测试使用地震计为BBVS-60(观测频带为60 s/0.017 Hz~50 Hz)和CMG-3TB(观测频带为120 s/0.008 Hz~50 Hz),后者在低频端更低,各测试场地的连续记录波形形态可能略有差别,但两种地震计均能覆盖大部分噪声频段,数据采集器使用大动态范围的24位设备,且仪器参数设置为能尽量记录小的信号,保证数据采集的有效性。计算软件使用北京港震机电有限公司开发的“Cal79_20160301C_user_地震数据分析处理系列软件包”[9],实际使用时需精确填写系统灵敏度等参数以实现扣除仪器响应。通过以上勘选测试设备和处理软件,可对测试台基垂直、南北与东西向的地脉动速度值进行测定,能够实时查看数据波形记录、存储观测数据,并计算噪声功率谱用以评估地噪声水平。
-
为了较为客观地评估测试场地的环境地噪声水平,同时为降低人类活动和其他地表干扰,在安放宽频带地震计进行测试前,一般均须清理测试场地的地表风化层。宽频带地震计安放位置以底部接触到坚实基岩为佳,或者坚实硬土层,测试基坑深度不小于50 cm,测试基坑直径以容纳地震计防护罩为宜,注意在地震计防潮保温及布设的线缆保持松弛状态下进行填埋式观测。地震计开锁后进行调零和系统标定操作,标定计算结果正常后可进行连续观测记录。在以上专业处理的基础上,对台站进行48 h(三通道连续记录,100点/s采样率)连续记录。并使用edas-ias-view-1.2软件浏览记录波形,查验勘选期间记录到的事件、统计非天然地震事件发生频度N和非天然地震事件持续时间占记录时间的百分比R[10]。在浏览地脉动波形记录时, 需要根据实际情况,判断台基是否有受到干扰,并及时判断干扰来源,对于观测系统所带来的干扰,要通过反复检查观测设备,及时摒去偶发干扰。锡林郭勒区域的台站勘选工作共收集了62个台站的48 h连续地脉动数据。
-
分别对锡林郭勒区域科学台阵62个子台勘选的48 h连续地脉动噪声数据进行了处理与分析,挑选出未记录到地震事件且无显著性偶发干扰的噪声数据,计算各个台站台址勘选时带宽在1~20 Hz(1/3倍频程)的地动速度均方根值(RMS值)、噪声功率谱密度值[11-14],分析各台站台址勘选背景噪声情况。
-
1)地动速度均方根值
地脉动噪声的速度均方根值(RMS)是反映台址背景噪声水平的一个重要指标[7, 15-16],该值的计算公式为:
$$RMS = \sqrt {\frac{{\displaystyle\sum\limits_{i = 1}^n {{{({V_i} - \bar V)}^2}} }}{n}} ,\bar V = \frac{1}{n}\sum\limits_{i = 1}^n {{V_i}} 。$$ (1) 式中:
$\bar V$ 为实测地脉动速度均值(m/s);Vi为某点实测地脉动速度值(m/s)。根据公式(1),计算地脉动噪声的速度均方根值。数字地震计记录的原始值为Counts值,因此要将Counts值转换为地脉动的速度值,二者间的转换公式为:$$V = \frac{{N \times {K_{{\rm{ad}}}}}}{{G \times {S_{{\rm{s}}0}}}}。$$ (2) 式中,V为实测地脉动速度(m/s);N为实际记录噪声的数字值;U为模拟输入峰值电压(V);R为仪器分辨率(Counts); Kad为数采实际增益;Ss0为地震仪工作灵敏度(V·s /m)。
根据公式(1)与公式(2),可计算出台站每小时的地脉动噪声RMS值[17],并取48 h地脉动噪声RMS值的算术平均值作为台站台址勘选时的RMS均值。
2)噪声功率谱
功率谱密度(Power Spectral Density,简称PSD) 是可用于衡量地脉动噪声水平的另一个重要指标[7, 11, 18-19],一般用实函数|
$F(\omega )$ |2表示,其计算公式为:$$F(\omega ) = \int\limits_{ - \infty }^{ + \infty } {f(t){e^{ - i\omega t}}} dt。$$ (3) 在物理学中,功率谱密度的意义是指在以
$\omega $ 为中心的单位频谱宽度内,信号的频率分量对功率的贡献[7]。一般通过对f(t)作傅立叶变化得到$F(\omega )$ ,然后再求取|$F(\omega )$ |2,或先求取f(t)的自相关函数,然后计算其频谱值[20]。为更加形象地凸显噪声在较长时间内的变化特征,通常以曲线图的形式展示噪声功率谱密度,即噪声功率谱密度曲线图。该图中一般会绘出地球高噪声模型(NHNM)曲线与低噪声模型(NLNM) 曲线。
基于以上原理,利用北京港震机电技术有限公司提供的地震数据分析处理系列软件,可计算台站勘选时的噪声功率谱密度(PSD)。
3)台站环境地噪声水平等级划分
按照国标《地震台站观测环境技术要求 第1 部分:测震》中对测震台站台基环境地噪声水平的相关规定,测震台站环境地噪声水平可分为5级[21],划分依据是台基背景噪声水平,即噪声的速度RMS值[22],主要通过Bnlv值与BnldB值体现(表1)。其中,Bnlv为由速度记录数据计算得到的RMS值,而分贝值BnldB则由BnlV计算得到,二者间转化公式为:
表 1 测震台站环境地噪声水平分级标准(依据台基处背景噪声水平)
台站环境地噪声水平等级 Ⅰ级 Ⅱ级 Ⅲ级 Ⅳ级 Ⅴ级 BnldB/dB <−150 [−150, −140) [−140, −130) [−130, −120) [−120, −110) BnlV/×10−8 m/s <3.16 [3.16, 10) [10, 31.6) [31.6, 100) [100, 3 160) $$Bn{l_{{\rm{dB}}}} = 10\lg \left[ {Bn{l_{\rm{v}}}/1{{({\rm{m}} \cdot {{\rm{s}}^{ - 1}})}^2}/{\rm{H}}{{\rm{z}}^{ - 1}}} \right]$$ (4) -
通过对锡林郭勒地区科学台阵62个子台台址勘选48 h连续观测数据的处理分析,计算得到了各子台环境地噪声UD、EW与NS三分向的 RMS 平均值与功率谱密度平均值,并将所有台站的台基地噪声水平进一步分为4类。
由计算得到的所有台站噪声功率谱密度值,绘制出了各台站的噪声功率谱曲线。因涉及台站较多,只选择了部分具有代表性的台站,展示其噪声功率谱曲线(图2)。根据图2可以看出,台站的噪声功率谱大都分布于高噪声模型(NHNM)曲线与低噪声模型(NLNM) 曲线之间,无明显偏离高、低噪声模型曲线形态的现象。基于功率谱密度值,并根据各子台的RMS 计算结果以及台站环境地噪声分级标准(表1),锡林郭勒地区科学台阵的62个子台可划分为4类(图3)。Ⅰ级台站(Bnlv< 3.16×10−8 m/s)共有4个;Ⅱ级台站(3.16×10−8 m/s≤Bnlv<1.0×10−7 m/s)有25个;Ⅲ级台站(1.0×10−7 m/s≤Bnlv<3.16×10−7 m/s)有31个;Ⅳ级台站(3.16×10−7 m/s≤Bnlv<1.0×10−6 m/s)有3个。其中,Ⅱ级台站中台站代码分别为15B06、15A89的2个子台,其RMS值仅稍大于Ⅰ级台站标准值,Ⅲ级台站中台站代码分别为15A64、15A54、15A68、15B06的子台,其RMS值也只稍高于Ⅱ类台站标准值,这些台站可在后期建设过程中通过改善或优化相应防干扰设施,使其台基质量进一步提高一个级别,优于勘选时水平。
根据62个子台在连续观测48 h内台基背景噪声的RMS均值与功率谱密度均值的分布认为,锡林郭勒地区科学台阵在台址勘选阶段,绝大多数台站的背景地噪声水平较为良好。在分析观测数据时,也发现有个别时间段内的噪声水平较高,这可能是在勘选过程中某个时段存在较大干扰源所致,结合台址勘选时周围环境,认为应是台址附近草场牛羊群活动,以及牧民驾驶摩托车等交通工具所产生的偶发噪声影响了观测数据的质量。
基于以上计算结果与研究分析,内蒙古锡林郭勒地区科学台阵建成后,将有半数以上台站的台基质量均可满足甚至高于Ⅱ级台站标准,这反映出锡林郭勒地区科学台阵的勘选工作取得了预期的成果。
-
根据对锡林郭勒地区科学台阵台址勘选时的背景噪声分析,得出以下结论:
1)通过对62个子台台址勘选4 h连续观测数据的处理分析, 计算得到了各子台背景噪声UD、EW与NS三分向的RMS平均值与功率谱密度平均值,并依据台站环境地噪声水平分级标准,将所有台站环境地噪声分为四级:Ⅰ级、Ⅱ级、Ⅲ级与Ⅳ级。从噪声功率谱曲线图中可看出,绝大多数台站(约85%)的噪声功率谱密度曲线大都分布于高噪声模型(NHNM)曲线之下,认为锡林郭勒地区科学台阵在台址勘选阶段,绝大多数台站的背景噪声水平较为良好。
2)根据各子台在48 h噪声RMS均值与功率谱密度均值的变化情况,以及台址勘选时的周围环境,认为勘选过程中干扰观测环境的噪声源主要来自于台址附近草场牛羊群活动,以及牧民驾驶摩托车等交通工具所产生的噪声。
3)根据研究结果,认为内蒙古锡林郭勒地区科学台阵建成后,将有半数以上台站其台基环境地噪声水平满足甚至高于Ⅱ级标准,这说明内蒙古锡林郭勒地区科学台阵的勘选工作取得了预期的勘选成果。
4)根据地震台站观测环境技术要求,不同环境地噪声水平的台站安装地震仪的类型也有所不同。内蒙古自治区为中国大陆背景地噪声区域划分中的A类地区,因此锡林郭勒区域科学台阵中环境地噪声水平为Ⅰ级、Ⅱ级的台站可依次安放甚宽频带地震仪、宽频带地震仪;对于环境地噪声水平为Ⅲ级与Ⅳ级的台站,建台时须采取有效降低环境地噪声干扰的措施,使环境地噪声水平至少达到Ⅱ级,以利于安装宽频带地震仪或短周期数字地震仪。
5)在实际建台时,建议对环境地噪声水平为Ⅲ级与Ⅳ级的台站进行降低环境地噪声干扰的相关措施:①在勘选点附近适当调整台站位置,优选山坡中部可避免风力、牛羊群等干扰的地方;②提高部分台站安放地震计的埋藏深度,以减少气流等外界环境的干扰;③可为地震计添加防护罩,使地震计的工作环境尽可能保持恒适温、稳气流的状态。
Background Noise Analysis of the Site Exploration of Scientific Seismic Array in Xilinguole Region, Inner Mongolia
-
摘要: 为评价内蒙古科学台阵各子台的环境地噪声水平,利用锡林郭勒区域62个子台台址勘选时的环境地噪声数字记录,计算分析各台站环境地噪声水平;根据台阵各子台勘选地脉动噪声的速度均方根值、噪声功率谱密度等数据以及相关国标规定,得出锡林郭勒区域62个子台台址勘选环境地噪声水平较为理想,均在标准范围内,并将其划分为4个等级,认为各子台的选址结果较好。通过对科学台阵台址勘选背景噪声的分析,不仅有利于科学台阵建成以后各子台场地响应及震级校正研究的开展,同时可对内蒙古科学台阵地震监测能力的评估提供基础资料。Abstract: In order to evaluate the site-selecting quality of each seismic station, the background noise of each station is calculated and analyzed by using the digital record of the background noise of the 62 stations in Xilinguole area. According to the calculating results of root- mean-square (RMS) and power spectrum density of these seismic stations, and combining with the international specification for digital seismic observation techniques, it is concluded that the background noise level of each station is relatively well, they are all within the standard range and are divided into four levels. The analysis of the background noise of the site selection of the scientific seismic array is not only conducive to the study of the site response and magnitude correction of these stations after the construction of the scientific seismic array, but also to the evaluation of the seismic monitoring capability of the scientific seismic array in Inner Mongolia.
-
表 1 测震台站环境地噪声水平分级标准(依据台基处背景噪声水平)
台站环境地噪声水平等级 Ⅰ级 Ⅱ级 Ⅲ级 Ⅳ级 Ⅴ级 BnldB/dB <−150 [−150, −140) [−140, −130) [−130, −120) [−120, −110) BnlV/×10−8 m/s <3.16 [3.16, 10) [10, 31.6) [31.6, 100) [100, 3 160) -
[1] 寿海涛. 上海台阵精确定位及应用[D]. 合肥: 中国科学技术大学, 2005. [2] 陈建涛, 谢剑波, 劳谦. 广东阳江小孔径井下型地震监测台阵介绍[J]. 华南地震, 2017, 37(1): 42-48. [3] 佟玉霞. 地震台阵的设计与数据处理研究[D]. 北京: 中国地震局地球物理研究所, 2003. [4] 彼得·鲍曼. 新地震观测实践手册[M]. 中国地震局监测预报司, 译. 北京: 地震出版社, 2006. [5] Rost S, Thomas C. Array seismology: methods and applications[J]. Reviews of Geophysics, 2002, 40(3): 1008. doi: 10.1029/2000RG000100 [6] 额日德木图. 浅谈锡林郭勒盟地质构造及地震概况[J]. 内蒙古科技与经济, 2002(12): 131. [7] 吴永权, 林伟, 何寿清. 广东省"十五"项目测震台站台址勘选结果分析[J]. 地震地磁观测与研究, 2007, 28(1): 90-95. doi: 10.3969/j.issn.1003-3246.2007.01.015 [8] 陈继锋. 兰州台网记录的青藏高原影区地震的分析和研究[D]. 兰州: 中国地震局兰州地震研究所, 2009. [9] 李晓东, 李姗姗, 周远新, 等. 库尔勒国家数字地震台近震震级偏差分析[J]. 地震地磁观测与研究, 2017, 38(4): 45-48. doi: 10.3969/j.issn.1003-3246.2017.04.008 [10] 胡文灼, 林伟, 杨选, 等. 曲江台野外试记分析[J]. 防灾科技学院学报, 2007, 9(3): 77-80. doi: 10.3969/j.issn.1673-8047.2007.03.020 [11] 郑培玲, 陈力刚, 张亚琳, 等. 河南数字台网子台台基背景噪声分析[J]. 防灾减灾学报, 2014, 30(1): 56-59. doi: 10.3969/j.issn.1674-8565.2014.01.011 [12] 李媛媛, 吴东. 山西数字遥测地震台网各子台台址地脉动噪声分析[J]. 内陆地震, 2004, 18(2): 182-186. doi: 10.3969/j.issn.1001-8956.2004.02.013 [13] 陈祥熊, 陈文明, 江燕, 等. 福建数字地震台网台址背景噪声分析[J]. 地震地磁观测与研究, 1998, 19(5): 13-22. [14] 何彦, 王宝柱, 宋太成, 等. 新疆数字地震台站观测动态范围和台基噪声的分析[J]. 内陆地震, 2006, 20(2): 178-182. doi: 10.3969/j.issn.1001-8956.2006.02.013 [15] 安祥宇, 赵倩, 王承伟, 等. 辽宁新建台站勘选过程概述及结果分析[J]. 防灾减灾学报, 2016, 32(3): 74-81. [16] 杨晶琼, 杨周胜. 云南地震数字遥测台网子台地动噪声分析[J]. 地震研究, 2005, 28(1): 86-89. doi: 10.3969/j.issn.1000-0666.2005.01.017 [17] 石玉燕, 颜启, 穆娟, 等. 山东“十五”测震台网子台地动躁声初步分析[J]. 内陆地震, 2008, 22(1): 73-77. doi: 10.3969/j.issn.1001-8956.2008.01.010 [18] 谢江涛, 林丽萍, 谌亮, 等. 地震台站台基噪声功率谱概率密度函数Matlab实现[J]. 地震地磁观测与研究, 2018, 39(2): 84-89. doi: 10.3969/j.issn.1003-3246.2018.02.012 [19] 王宝柱, 冯磊. 新疆测震台站台址的勘选及影响因素分析[J]. 内陆地震, 2014, 28(3): 254-260. doi: 10.3969/j.issn.1001-8956.2014.03.009 [20] 陈长兴, 高晓光, 刘昌云. 基于傅立叶变换的积分算法研究[J]. 系统仿真学报, 2004, 16(4): 817-819. doi: 10.3969/j.issn.1004-731X.2004.04.062 [21] GB/T 19531.1-2004, 地震台站观测环境技术要求第1部分: 测震[S]. 北京: 中国标准出版社, 2004. [22] 张玲, 薛振岳, 高跃雄. 山西“十五”测震台网子台环境地噪声分析[J]. 地震地磁观测与研究, 2012, 33(5/6): 225-230. -