• ISSN 1673-5722
  • CN 11-5429/P

常熟台倾斜观测资料的异常特征分析

狄樑 陆德明 丁建国 刘冬冬 钱文杰

宫杰,张朋,张扬,张敏,王佳,居海华,2022. 江苏省测震台网井下地震计方位角检测. 震灾防御技术,17(1):181−187. doi:10.11899/zzfy20220119. doi: 10.11899/zzfy20220119
引用本文: 狄樑, 陆德明, 丁建国, 刘冬冬, 钱文杰. 常熟台倾斜观测资料的异常特征分析[J]. 震灾防御技术, 2017, 12(2): 415-422. doi: 10.11899/zzfy20170219
Gong Jie, Zhang Peng, Zhang Yang, Zhang Min, Wang Jia, Ju Haihua. Borehole Seismometer Azimuth Detection in Jiangsu Seismic Network[J]. Technology for Earthquake Disaster Prevention, 2022, 17(1): 181-187. doi: 10.11899/zzfy20220119
Citation: Di Liang, Lu Deming, Ding Jianguo, Liu Dongdong, Qian Wenjie. Analysis of Anomaly Characteristics of Tilt Measuremennt Data from Changshu Seismic Station[J]. Technology for Earthquake Disaster Prevention, 2017, 12(2): 415-422. doi: 10.11899/zzfy20170219

常熟台倾斜观测资料的异常特征分析

doi: 10.11899/zzfy20170219
基金项目: 

中国地震局"三结合"课题 161004

详细信息
    作者简介:

    狄樑, 男, 生于1972年。工程师。主要从事地震前兆观测与研究。E-mail:jscswyn123456789@163.com

Analysis of Anomaly Characteristics of Tilt Measuremennt Data from Changshu Seismic Station

  • 摘要: 在研究常熟台水管仪与垂直摆多年倾斜观测资料的基础上,利用潮汐因子调和分析、Nakai拟合检验、小波分析等方法对台站周边约250km范围内发生的高邮-宝应MS 4.9及如东近海MS 3.8等地震前的异常信息进行了提取,对数据异常分析方法进行评判,旨在为以后地震前及时发现地震提供参考。
  • 随着社会城市化发展进程的不断加快,地震监测环境面临严峻挑战,井下地震观测凭借其恒温、气流稳定、信噪比高等优势逐渐得到应用(宫杰等,2019樊晓春等,2020贾媛等,2020)。相比地面地震计,井下地震计受台站外部条件、磁性环境、定位设备精度、地震计底座固定情况、摆线应力释放等因素制约,方位角检测结果往往存在较大偏差,对利用地震观测数据开展各向异性、面波频散、接收函数及震源机制解等研究产生影响(陈继锋等,2016李少睿等,2016张明等,2019),因此针对井下地震计方位角的检测与校正尤为重要。

    目前,井下地震计方位角计算分析方法主要包括基于P波质点偏振原理的多地震事件信噪比加权叠加方法及基于地面井下双地震计对比观测的相关分析法(Niu等,2011李少睿等,20122016谢剑波,2014)。Aster等(1991)利用近震初至P波质点运动计算获得精度5°以内的井下地震计方位角,李少睿等(2016)通过开展不同频带地震计、不同台站址及不同距离的大量对比试验,结果表明相关分析法可使井下地震计计算精度达4°以内。

    2016年,江苏省地震局利用数据波形相关分析法对省内14个深井台站开展井下地震计方位角普查工作。此外,利用数据波形相关分析法对2018年起江苏地区新建井下地震计进行方位角检测存档。本文梳理目前江苏地区部分正在运行井下地震计方位角结果,并展开讨论,以新建溧阳2地震台站井下地震计为例,论述计算结果正确性,最后依托数据波形相关分析法,结合实际方位角检查过程,开展滤波频带研究。

    目前,江苏地区共运行35套井下地震设备,其中宝应地震台站采用CMG-3TB型120 s~50 Hz甚宽频地震计,如皋地震台站采用TBG-60B型60 s~50 Hz一体化宽频带地震计,其余均为GL-S60B型60 s~50 Hz宽频带地震计。沛县、淮安、大丰、海安、阳光岛深井台站均由JDF系列宽频带地震计更改为GL-S60B型地震计,由于原有设备底座与更改后井下地震计密封桶相匹配,因此保留原有设备底座,改造后安装深度不变。本文共涉及22个井下地震计方位角检测,基本信息如表1所示。其中,溧阳地震台站、兴化地震台站原井下设备暂且保留。

    表 1  井下地震台基本信息
    Table 1.  The fundamental information of borehole seismic stations
    台站名代码台基岩性安装深度/m安装方式地震计型号频带范围数采型号采样率/Hz
    沛县PX石英砂280卡壁JDF-220 s~20 HzEDAS-24IP100
    280卡壁GL-S60B60 s~50 HzEDAS-24GN100
    淮安HUA松散沉积层315卡壁JDF-220 s~50 HzEDAS-24IP100
    315卡壁GL-S60B60 s~50 HzEDAS-24GN100
    海安HA松散沉积层425卡壁JDF-3120 s~50 HzEDAS-24GN100
    425卡壁GL-S60B60 s~50 HzEDAS-24GN100
    大丰DF松散沉积层409卡壁JDF-3120 s~50 HzEDAS-24GN100
    409卡壁GL-S60B60 s~50 HzEDAS-24GN100
    阳光岛YGD基岩530卡壁JDF-3120 s~50 HzEDAS-24GN100
    530卡壁GL-S60B60 s~50 HzEDAS-24GN100
    泰州TZ沉积岩500落底GL-S60B60 s~50 HzEDAS-24GN100
    盐城YC松散沉积层410卡壁GL-S60B60 s~50 HzEDAS-24GN100
    宝应BY石英砂450卡壁CMG-3TB120 s~50 HzEDAS-24GN100
    射阳SY松散沉积层460卡壁GL-S60B60 s~50 HzEDAS-24GN100
    如东RD灰岩450落底GL-S60B60 s~50 HzEDAS-24GN100
    启东QD灰岩410落底GL-S60B60 s~50 HzEDAS-24GN100
    溧阳LY安山玄武岩80落底GL-S60B60 s~50 HzEDAS-24GN100
    溧阳2LY2安山玄武岩203落底GL-S60B60 s~50 HzEDAS-24GN100
    高邮GY松散沉积层452落底GL-S60B60 s~50 HzEDAS-24GN100
    南通NT石英砂159卡壁GL-S60B60 s~50 HzEDAS-24GN100
    兴化XH混合花岗岩490落底GL-S60B60 s~50 HzEDAS-24GN100
    兴化2XH2混合花岗岩510卡壁GL-CS60B60 s~50 HzEDAS-24GN100
    涟水LAS石英岩400落底GL-S60B60 s~50 HzEDAS-24GN100
    丰县FX石英砂405落底GL-S60B60 s~50 HzEDAS-24GN100
    李堡LIB松散沉积层450落底GL-CS60B60 s~50 HzEDAS-24GN100
    坪山PIS泥岩425落底GL-S60B60 s~50 HzEDAS-24GN100
    扬中YAZ泥岩416落底GL-CS60B60 s~50 HzEDAS-24GN100
    下载: 导出CSV 
    | 显示表格

    相关分析法广泛应用于描述2个随机变量之间的相似程度。首先在待测深井台站井口附近地面安装地面地震计,并以寻北仪精确指北,此地面地震计作为参考地震计与井下待测地震计进行相关分析:

    $$ x={x}_{1}\mathrm{cos}\varphi +{y}_{1}\mathrm{sin}\varphi $$ (1)
    $$ y=-{x}_{1}\mathrm{sin}\varphi +{y}_{1}\mathrm{cos}\varphi $$ (2)
    $$ {R}_{x{x}{'}}=\frac{\displaystyle\sum (x-\overline{x})({x}{'}-\overline{{x}{'}})}{\sqrt{\displaystyle\sum {(x-\overline{x})}^{2}}\sqrt{\displaystyle\sum {({x}{'}-\overline{{x}{'}})}^{2}}} $$ (3)
    $$ {R}_{y{y}{'}}=\frac{\displaystyle\sum (y-\overline{y})({y}{'}-\overline{{y}{'}})}{\sqrt{\displaystyle\sum {(y-\overline{y})}^{2}}\sqrt{\displaystyle\sum {({y}{'}-\overline{{y}{'}})}^{2}}} $$ (4)

    式中,$ {x}_{1} $$ {y}_{1} $分别为井下待测地震计EW和NS向分量数据,$ x $$ y $分别为井下待测地震计旋转$ \varphi $角度后EW和NS向分量数据,$ {x}{'} $$ {y}{'} $分别为地面参考地震计EW和NS向分量数据,$ {R}_{x{x}{'}} $$ {R}_{y{y}{'}} $分别为$ x $$ {x}{'} $$ y $$ {y}{'} $的相关系数(李少睿等,2020)。

    井下地震计方位角检测时,应尽量选择与井下地震计同频带的地震计作为地面参考,因地面与井下同频带地震计在滤波频段内相频特性曲线具有较小的相位差,可直接使用二者观测数据进行相关分析,获取较高精度的方位角检测结果。而当地面参考地震计与井下待测地震计观测频带相差较大时,若直接利用二者观测数据进行相关分析,滤波频段内相频特性曲线相位差将导致较大的方位角检测偏差,甚至出现180°的反向结果(李少睿等,2016)。此时应利用较短周期地震计传递函数,将地面与井下地震计中较宽频带地震计观测数据进行仿真处理。本文采用频域滤波法进行数据预处理,因低频段信号传播距离长、对比测试仪器记录相关性高,因此仅考虑地震观测系统低频段地震计特性,不考虑数据采集器高频段影响,计算公式如下:

    $$ y\left(t\right)=x\left(t\right)\cdot h\left(t\right) $$ (5)
    $$ Y\left(\omega \right)=X\left(\omega \right)\cdot H\left(\omega \right) $$ (6)
    $$ {Y}_{1}\left(\omega \right)={H}_{1}\left(\omega \right)\cdot \frac{Y\left(\omega \right)}{H\left(\omega \right)} $$ (7)

    式中,$ y\left(t\right) $为较宽频带地震计记录数据;$ x\left(t\right) $为地面真实运动;$ h\left(t\right) $为较宽频带地震计脉冲响应;$ Y\left(\omega \right) $$ X\left(\omega \right) $分别为$ y\left(t\right) $$ x\left(t\right) $的傅里叶变换;$ H\left(\omega \right) $为较宽频带地震计频域传递函数,$ {H}_{1}\left(\omega \right) $为较短周期地震计频域传递函数;$ {Y}_{1}\left(\omega \right) $为仿真数据的傅里叶变换谱。

    通过对$ {Y}_{1}\left(\omega \right) $进行傅里叶逆变换,可得到预处理后的时域数据。

    需强调的是,地面GL-S120型甚宽频地震计与井下GL-S60B型宽频带地震计在滤波频带内相频特性曲线具有较高的一致性,几乎不存在相位差,所以两者组合对比观测相关分析无需进行仿真处理。宝应、海安、大丰、阳光岛甚宽频井下地震计及部分GL-S60B型宽频带井下地震计选择GL-S120型甚宽频地震计作为地面参考地震计,无需进行仿真处理;剩余部分GL-S60B型宽频带井下地震计选择同频带地面GL-S60型宽频带地震计作为地面参考地震计,无需进行仿真处理。沛县、淮安JDF-2型宽频带井下地震计选择GL-S60型宽频带地震计作为地面参考地震计,因观测频带差异,需根据二者传递函数对GL-S60型地震计观测数据进行仿真处理。

    为获取较准确的井下地震计方位角,仅需对地面井下地震计观测数据进行滤波处理,以求井下待测地震计旋转一定角度后与地面参考地震计最大相关性。滤波频带的选择直接影响方位角计算精度。陆地与海洋相互作用在1~10 s频段存在明显波峰且较稳定(Peterson,1993),Lacoss等(1969)研究表明地脉动在0.2~0.3 Hz频段主要包含高阶模式瑞利面波及体波,而高阶模式瑞利面波主要能量集中于约半个波长深度内,基本覆盖井下地震计安装深度,因此在地面与井下地震计间具有较高相关性。本文主要以0.2~0.3 Hz作为滤波频段,开展井下地震计方位角检测。

    井下地震计方位角检测过程中,地面参考地震计寻北采用NV-NF301型寻北仪,精度达0.3°。地面参考地震计采用GL-S120型甚宽频地震计和GL-S60型宽频带地震计,具体参数如表2所示,分别配备的2套EDAS-24GN型数据采集器量程为±10 V,采样率为100 Hz,最小相位滤波,转换因子为1 192 nV/count。

    表 2  地面参考地震计参数
    Table 2.  The parameters of ground referenceseismometer
    地震计型号地震计序列号频带范围电压灵敏度/ V·m−1·s−1
    UD分向EW分向NS分向
    GL-S120G14408VS120 s~50 Hz1 993.761 986.201 994.00
    GL-S60G11842VS60 s~50 Hz2 002.262 000.322 005.64
    下载: 导出CSV 
    | 显示表格

    考虑地震计内部水平向分量正交性误差等因素,单台站检测结果以24 h的2个水平向分量均值为准,井下地震计方位角检测结果如表3所示。

    表 3  井下地震计方位角检测结果
    Table 3.  The azimuth detection results of borehole seismometers
    台站名测试地震计型号参考地震计型号滤波频带/Hz参考地震计方位角/°测试地震计方位角/°校正后井下方位角/°相关系数
    沛县JDF-2GL-S600.2~0.3359.9186.2186.10.921 89
    淮安JDF-2GL-S600.2~0.3359.6166.7166.30.895 56
    GL-S60BGL-S600.2~0.30.150.250.10.992 51
    海安JDF-3GL-S1200.2~0.3359.9146.2146.10.814 48
    GL-S60BGL-S600.2~0.3359.9314.8314.70.997 47
    大丰JDF-3GL-S1200.2~0.3359.8184.1183.90.831 92
    GL-S60BGL-S600.2~0.30.1122.3122.40.990 79
    阳光岛JDF-3GL-S1200.2~0.30.2162.4162.60.800 07
    GL-S60BGL-S600.2~0.3359.6340.6340.20.960 31
    泰州GL-S60BGL-S1200.2~0.30.0−4.5−4.50.928 45
    盐城GL-S60BGL-S1200.2~0.30.0−137.3−137.30.900 50
    宝应CMG-3TBGL-S1200.2~0.3359.86.86.60.963 71
    射阳GL-S60BGL-S600.1~0.2359.6346.0345.60.968 96
    如东GL-S60BGL-S600.2~0.30.3329.9330.20.937 37
    启东GL-S60BGL-S600.2~0.30.06.26.20.823 46
    溧阳GL-S60BGL-S600.2~0.30.3−2.2−1.90.917 73
    溧阳2GL-S60BGL-S600.2~0.30.0176.5176.50.997 48
    高邮GL-S60BGL-S600.2~0.3359.8299.5299.30.975 95
    南通GL-S60BGL-S1200.2~0.3359.5347.4346.90.997 02
    兴化GL-S60BGL-S600.2~0.3359.384.383.60.988 24
    兴化2GL-CS60BGL-S600.2~0.30.0299.9299.90.990 08
    涟水GL-S60BGL-S600.2~0.30.0233.3233.30.938 00
    丰县GL-S60BGL-S600.2~0.30.0301.8301.80.956 90
    李堡GL-CS60BGL-S600.2~0.3359.6147.5147.10.998 01
    坪山GL-S60BGL-S600.2~0.30.186.686.70.979 15
    扬中GL-CS60BGL-S600.2~0.30.6211.8212.40.996 99
    下载: 导出CSV 
    | 显示表格

    表3可知,井下地震计方位角受多因素制约,普遍存在较大偏差,其中沛县、淮安、海安、大丰、阳光岛深井台站使用的JDF系列宽频带地震计方位角偏差约为180°,大致反向,这可能与其安装工艺有关。需强调的是,上述5个水平向分量极性反向的JDF系列井下宽频带地震计方位角检测结果可结合测震台网观测数据远震P波极性分析结果得到验证。此外,这5个深井台站经更改后虽使用原有设备底座,但新地震计放入密封桶方位校正后安装在原有底座上,井下地震计方位角必定发生变化。因此,在深井台站井下地震计更改或维修时,如重新下井,须重新开展井下地震计方位角检测工作,及时更新、掌握运行状态。

    为深入研究JDF-2井下地震计方位角检测情况,给出沛县及淮安台站仿真与未仿真时24 h平均检测结果,如表4所示。由表可知,未仿真情况下,沛县、淮安台站JDF-2井下地震计方位角检测结果相关系数平均值约为0.83;南北向分量相关系数较大,达0.9;东西向分量相对相关系数较小,约为0.7。依据GL-S60、JDF-2型地震计传递函数对GL-S60型地震计观测数据进行仿真处理,相关系数平均值提高至0.9左右,东西向分量相关系数提高至0.8以上。利用不同频带地震计进行井下地震计方位角检测时,对较宽频带地震计进行仿真处理尤为重要。

    表 4  沛县、淮安台站24 h平均检测结果
    Table 4.  24-hour average azimuth detection results of Peixian、Huaian
    台站状态东西向分量南北向分量平均值
    相对方位角/º相关系数相对方位角/º相关系数相对方位角/º相关系数
    沛县仿真185.30.882 14187.00.961 63186.10.921 89
    未仿真188.00.697 36186.00.962 12186.80.829 74
    淮安仿真164.70.846 90168.60.944 21166.30.895 56
    未仿真160.10.729 65169.10.955 54164.60.842 60
    下载: 导出CSV 
    | 显示表格

    井下地震计方位角检测过程中,基本采用0.2~0.3 Hz滤波频带,但射阳台站除外。当采用0.2~0.3 Hz滤波频带时,射阳台站井下地震计方位角检测相关系数约为0.69,改为0.1~0.2 Hz滤波频带后,相关系数提高至约0.97。

    表3可知,利用相关分析法计算溧阳2新建GL-S60B型井下宽频带地震计方位角为176.5°,水平向分量几乎反向。为证明相关分析法在井下地震计方位角检测中的可信性,以溧阳2新建井下地震计为例,利用初至震相相对清晰的远震所在时段进行相关性分析,并结合P波进行极性分析,以进行直观比较。本文选取2020年5月6日印尼班达海7.2级地震,远震所处时段地震计方位角为176.8°,相关系数达0.998 86,地面参考地震计及井下待测地震计记录的原始波形如图1所示。

    图 1  印尼班达海7.2级地震原始波形
    Figure 1.  Indonesia Banda Sea M7.2 earthquake

    地面参考地震计经寻北仪指北校准,由图1可知地面与井下地震计初至震相清晰,其中垂直向分量具有较高相似性,初至P波同时向上;水平向分量出现差异,初至震相就振幅而言同时表现为东西向分量小、南北向分量大,但初至方向不同,地面参考地震计东西向分量初至向下,南北向分量向上,而井下待测地震计东西向分量初至向上,南北向分量向下。地面与井下地震计远震初至P波水平向分量几乎反向,这与相关分析法计算结果一致,由此印证了相关分析法应用于井下地震计方位角检测中具有较高可靠性,也验证本文利用相关分析法获取的江苏测震台网井下地震计方位角结果可信。

    当滤波频带改为0.1~0.2 Hz后,射阳台站井下地震计方位角检测相关系数得到大幅提高,可知滤波处理可获取较高的相关性,从而得到精度较高的方位角结果。台基噪声是评定地震观测环境质量的重要指标之一(宫杰等,2020),本文针对射阳台站井下地震计观测数据进行台基噪声分析,结合台基噪声功率谱密度(PSD)计算深入分析滤波频带与相关性的内在联系,选取射阳台站GL-S60B型井下宽频带地震计方位角2016年9月7日00时的检测数据,功率谱密度曲线如图2所示。

    图 2  射阳台站井下地震计台基噪声功率谱密度曲线
    Figure 2.  Borehole seismometer power spectrum curve of Sheyang station

    图2可知,射阳台站井下地震计所记录的台基噪声功率谱处于地球高噪声模型(NHNM)和地球低噪声模型(NLNM)之间,属于正常噪声水平。密度曲线垂直向分量在0.2~0.3 Hz频段内存在明显波峰;对于水平向分量,0.1~0.6 Hz频段内具有较高一致性,0.1~0.2 Hz频段内存在明显波峰,0.15 Hz左右功率谱密度显著高于两侧,0.2~0.4 Hz频段较平坦,缓慢升至0.6 Hz时达到波峰。由此可判断采用相关分析法计算井下地震计方位角时,地面与井下地震计相关性与滤波频带内台基噪声功率谱密度曲线形态有关。选择滤波频带内存在明显波峰的频段有利于提高地面与井下地震计相关性,滤波频带内功率谱密度曲线相对平坦会导致较低的相关性,从而影响方位角检测精度。射阳台站滤波频带不同于其他台站,可能与地质构造、台站环境等因素有关。对于大多数深井台站而言,选择0.2~0.3 Hz滤波频段开展井下地震计方位角检测是可行的,如出现相关性较低的情况,可从台基噪声入手,深入分析功率谱密度曲线形态,尽可能选择具有明显波峰的频段作为滤波频段。若无法解算台基噪声功率谱密度曲线,可分别选择0.1~0.2 Hz、0.2~0.3 Hz、0.3~0.4 Hz滤波频段,选取相关性较高的检测结果。

    江苏省测震台网深井台站较多,具有频段宽、安装深度大等特点,本文分析2016年以来井下地震计方位角检测结果,得出以下结论:

    (1)受台站外部条件、磁性环境、定位设备精度、地震计底座固定情况、摆线应力释放等因素制约,井下地震计方位角检测结果普遍存在较大偏差,包含JDF系列在内的部分井下地震计水平向分量甚至出现反向现象。无论是否更换底座,涉及提井时,井下地震计方位角必定发生变化,因此每次提井时须重新开展方位角检测工作,及时更新、掌握仪器运行状态。

    (2)就相关性而言,不同频带地震计组合检测井下地震计方位角,应利用二者传递函数对较宽频地震计进行仿真处理,有利于提高二者相关性,从而获得较精确的检测结果。当地面参考地震计与井下待测地震计同频带时,相关系数普遍较高。

    (3)以溧阳2台站新建井下地震计为例,通过远震P波极性分析法直观反映出地面与井下地震计水平向分量几乎反向,与相关分析法得到的结论一致,从而验证相关分析法应用于井下地震计方位角检测是正确可行的,也证明了本文获得的江苏省测震台网井下地震计方位角检测结果是真实可信的。

    (4)计算分析射阳台站井下地震计记录台基噪声功率谱密度曲线,可知应选择功率谱密度曲线具有明显波峰频段作为滤波频段,有利于获得较高的相关性,从而得到较精确的井下地震计方位角检测结果。

  • 图  1  常熟台水管仪垂直摆降雨干扰图

    Figure  1.  Interference diagrams of vertical pendulum andwater pipe system from Changshu seismic station

    图  2  常熟台垂直摆南北分量倾斜加速异常图(20100703—2010708(a)、20120807—20120809(b))

    Figure  2.  Acceleration anomaly of vertical tilt of the north and south from Changshu seismic station

    图  3  常熟台垂直摆南北分量M2波潮汐因子变化图

    Figure  3.  Variation of tidal factor of M2 wave of north-south component in vertical pendulum from Changshu seismic station

    图  4  常熟台水管仪垂直摆非潮汐应变加速度

    Figure  4.  Accelerations of non-tidal strain of water tube meter and vertical pendulum from Changshu seismic station

    图  5  常熟台水管仪垂直摆小波分解第7阶、8阶细节

    Figure  5.  The wavelet decomposition result at 7th and 8th level of water tube meter and vertical pendulum from Changshu seismic station

    表  1  常熟地震台周边5次震例目录

    Table  1.   The records of 5 earthquake examples around Changshu seismic station

    发震时间 震级/MS 震中距/km 发震地点 经度/° 纬度/°
    20100709 3.8 134 如东近海 121.36 32.30
    20100719 3.8 133 如东近海 121.36 32.27
    20120720 4.9 207 高邮-宝应 119.34 33.01
    20120825 4.3 252 南黄海 121.55 33.11
    20140118 3.7 132 如东近海 121.36 32.27
    下载: 导出CSV
  • 狄樑, 陈洪军, 冯志生等, 2003.常熟地震台形变观测仪数字化改造.地震地磁观测与研究, 24(4):92-95. http://www.cnki.com.cn/Article/CJFDTOTAL-DZGJ200304016.htm
    丁建国, 陆德明, 狄梁等, 2011.常熟台倾斜仪降雨干扰定量分析.华南地震, 31(3):83-88. http://www.cnki.com.cn/Article/CJFDTOTAL-HNDI201103010.htm
    刘建明, 李志海, 孙甲宁等, 2016.基于小波分析提取地倾斜异常特征.地震, 36(1):38-48. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZN201601005.htm
    王梅, 季爱东, 2000.山东省定点形变倾斜观测异常指标研究.地震研究, 23(4):398-404. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYJ200004005.htm
    张立, 牛安福, 汤曙恩等, 2005.大姚地震前连续形变短期前兆特征.地震, 25(2):75-82. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZN20050200A.htm
    张燕, 吴云, 刘永启等, 2003.潮汐形变资料中地震前兆信息的识别与提取.大地测量与地球动力学, 23(4):34-39. http://www.cnki.com.cn/Article/CJFDTOTAL-DKXB200304007.htm
    张雁滨, 蒋骏, 陈绍绪等, 2001.连续形变的前兆参量、判别方法及实用化研究.内陆地震, 15(1):1-10. http://www.cnki.com.cn/Article/CJFDTOTAL-LLDZ200101000.htm
  • 期刊类型引用(2)

    1. 周仁才,徐雨,刘杰,马千军. 基于多源传感器的电力新材料智慧仓储系统设计. 自动化与仪器仪表. 2024(03): 153-157 . 百度学术
    2. 姚君,邵将,朱金佗,宗威,李轩轩. 公共安全与保障视角下的应急救援装备工业设计及应用研究综述. 包装工程. 2024(18): 1-16+513 . 百度学术

    其他类型引用(2)

  • 加载中
图(5) / 表(1)
计量
  • 文章访问数:  130
  • HTML全文浏览量:  20
  • PDF下载量:  10
  • 被引次数: 4
出版历程
  • 收稿日期:  2016-08-29
  • 刊出日期:  2017-06-01

目录

/

返回文章
返回