JCZ-360超寬頻帶地震計實用性分析

朱厚林 陸棟梁 應允翔 陶方宇 周志 劉浩

摘 ?要：該課題旨在對金寨臺JCZ-360超寬頻帶地震計數據實用性進行詳細分析，首先，介紹了基于Matlab的JCZ-360超寬頻帶地震計LP通道固體潮提取方法，給出了VP寬頻帶傾斜儀功率譜計算方法，并將VP垂直擺數據同JCZ-360加速度輸出（LP通道）數據進行精確分析，全方位反映出該儀器的監測效能，事實上，其作為形變觀測儀具有更好的穩定性，避免了頻繁進山洞的人為干擾，JCZ-360完全可作為綜合觀測儀。

關鍵詞：超寬頻帶地震計 ?功率譜 ?JCZ-360 ?固體潮

中圖分類號：P315 ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標識碼：A ? ? ? ? ? ? ? ? ? 文章編號：1672-3791（2020）10（b）-0048-04

Abstract： This topic to JinZhai JCZ-360 ultra-broadband seismometer data practical detailed analysis. First of all， based on Matlab was introduced JCZ-360 ultra-wide band seismometer LP tidal channel extraction method， gives the VP wide-band inclinometer power spectrum calculation method， and the VP vertical pendulum acceleration output data with JCZ-360 data for precise analysis （LP）， all reflect the effectiveness of monitoring， in fact， as the deformation observation instrument has better stability， avoid the frequent into the cave of human disturbance， JCZ-360 can be used as a completely integrated viewer.

Key Words： Uwb seismometer; Power spectrum; JCZ-360; Earth tide

隨著儀器計量的標準化越來越規范，儀器種類越來越繁多，1997年，中國地震局地震研究所呂永清研究員通過研究發現JCZ-1的LP通道輸出可以觀測到地球自由震蕩。2007年，蔡亞先研究員在文章中提出JCZ-1T作為綜合觀測儀器的拓展運用取得了理想的成果。隨著產品的不斷優化，JCZ-360的性能更加穩定，目前安裝在安徽金寨地震臺試運行，目前儀器運行正常，且采用美國Q330HRS數采，目前采樣率使用的是100Hz，且安裝了2臺JCZ-360，因此要充分利用此數據，與金寨臺VP寬頻帶傾斜儀進行數據對比分析，驗證其綜合性價值。

1 ?數據提取與數據對比分析

1.1 JCZ-360超寬頻帶地震計LP通道固體潮提取

筆者利用Matlab提取了JCZ超寬頻帶地震計LP通道固體潮，提取程序如下：

Je11=ffit（（1：length（Je1）），Je1，1）;

figure; plot（Je11）

Je111=Je11*2.384e-6/20000

figure; plot（Je111）

其中Je1表示第一臺JCZ東西分量原始數據。提取后的固體潮見圖1。

提取的JCZ（LP通道）固體潮與VP固體潮對比分析見圖2。

由上圖可知對于近震的固體潮記錄，VP比JCZ記錄到的更為清晰，因為VP使用的擺是2s的擺，而JCZ超寬頻帶地震計使用的是360s的擺，試想一下 JCZ如果改為2s的擺記錄固體潮是不是更清楚，且遠比VP本身性能穩定？后續我們從功率譜圖中更能得到證明，其作為綜合觀測儀具有更好的穩定性，避免了頻繁進山洞的人為干擾，且地震計每月標定，滿足標定要求，且無量程限制。

1.2 形變儀器計算功率譜數據處理方法

現如今，對于形變儀器數據處理更多的是利用小波變換分析，現筆者介紹一種關于形變儀器功率譜計算方法，以供參考，反映儀器監測能力，數據可靠性是研究一切的基礎，通過功率譜的計算，更能夠反映出儀器性能，使我們相信儀器監測到的數據，且通過對JCZ功率譜、VP寬頻帶地震計功率譜的計算對比，找到儀器之間的通性、緊密性，事實上二者的數據是互為反映、相互驗證的。

如圖3可知：以VP寬頻帶傾斜儀為例，已知弧度為，值為vpx1，即形變儀器所測值，單位為ms（角秒），

故：

其中，vpx1即為VP寬頻帶傾斜儀的加速度值。

由此，我們可以計算出該儀器的功率譜。接下來執行：

vpx1=9.81*tan（vpx/1000/3600*pi/180）;

nfft=680280/2;

vpy1=9.81*tan（vpy/1000/3600*pi/180）;

[p11，fre1]=pwelch（vpx1，nfft，round（nfft*0.833），nfft，1）;[p12，fre2]=pwelch（vpy1，nfft，round（nfft*0.833），nfft，1）;

semilogx（fre1，10*log10（p11））;semilogx（fre2，10*log10（p12））;

hold on; fgnoise（0，1）

可得圖4VP擺NS分量、EW分量功率譜圖

1.3 金寨臺JCZ-360功率譜分析

圖5為Z臺JCZ-360LP通道固體潮功率譜圖，圖6是2臺JCZ和VP二分量功率譜對比圖。

2 ?總結分析

通過功率譜計算可知，JCZ-360加速度通道（LP）通道在1/360s之后未能記錄到清晰的固體潮，低頻記錄不到，而VP在0.5s之后不能記錄到高頻信號，這是由于擺體周期原因所致，從功率譜圖對比分析可知，形變儀器和地震計之間是相通的，JCZ-360超寬頻帶地震計記錄到的信息更全，不僅能夠清晰地記錄地球物理事件，更能記錄固體潮信息，且不會因為現如今形變儀器頻繁進洞調零、擺系故障造成數據不完整性，其作為綜合觀測儀更具有優勢，只需開發一套數據采集輸出系統即可。

參考文獻

[1] 張幸興，毛國良，常亮，等.易縣臺BBVS-120地震計記錄重力固體潮能力分析[J].華北地震科學，2019，37（2）：56-60，81.

[2] 張小飛，袁順，張嘉敏，等.庫爾勒地震臺JCZ-1T地震計LP通道固體潮數據分析[J].地震地磁觀測與研究，2016，37（4）：131-137.

[3] 陳珊樺，全建軍，蔡佩蕊，等.泉州地震臺M2166超寬頻帶地震儀記錄的地球球型自由振蕩信息研究[J].華南地震，2019，39（S1）：102-108.

[4] 曲保安，范曉勇，林秀娜，等.利用JCZ-1T地震計LP通道數據檢測斐濟M_W8.2深震激發的地球自由振蕩[J].地震地磁觀測與研究，2019，40（3）：87-92.

[5] 荊濤，孫藝，方禹心，等.JCZ-1T超寬頻帶地震計常見故障分析[J].地震地磁觀測與研究，2020，41（1）：85-90.

[6] 皮譽洋，周云耀，呂永清.Morlet變換在超寬頻帶地震儀記錄的汶川震前數據處理中的應用[J].大地測量與地球動力學，2018，38（9）：974-978.