將測震無量綱整數(shù)轉(zhuǎn)換為地動速度的快速方法

董立杰 丁 勇 趙乃千

（1、四川省地震局阿壩地震監(jiān)測中心站，四川 阿壩623300 2、四川省地震局成都地震監(jiān)測中心站，四川 成都 611730）

本文基于現(xiàn)有的轉(zhuǎn)換公式，以松潘地震臺測震記錄作為研究樣本，旨在研究一種快速的數(shù)據(jù)轉(zhuǎn)換方法，直接從數(shù)據(jù)本身獲得前放增益和地震計電壓靈敏度。通過這種方式獲得的參數(shù)更加準(zhǔn)確，并且精簡了前放增益和地震計電壓靈敏度查詢的步驟，提高了數(shù)據(jù)轉(zhuǎn)換的效率。

1 轉(zhuǎn)換公式

無量綱整數(shù)轉(zhuǎn)換為地動速度的公式為[1]：

式中，XV0為地動速度，AV0為振幅，單位為counts，Kad是A/D轉(zhuǎn)換因子，單位為μV/counts，G為前放增益，SS0為地震計電壓靈敏度，單位V/m·s-1。地震觀測系統(tǒng)不同，轉(zhuǎn)換系數(shù)也是有所不同的。以松潘地震臺為例，該地震臺目前采用的是CTS-1 型地震計和EDAS-24GN 型數(shù)據(jù)采集器的組合，A/D 轉(zhuǎn)換因子Kad等于1.192μV/counts。前放增益G和地震計電壓靈敏度SS0均為為常數(shù)，可將二者的乘積表示為一個計算因子S，則公式變?yōu)椋?/p>

2 S 值的計算

根據(jù)對S的定義可知，若要計算S 值，則需要查詢前放增益G 和地震計電壓靈敏度SS0的值。前放增益可以通過查詢相關(guān)文獻(xiàn)獲得，地震計電壓靈敏度存在默認(rèn)值。但是在實際運用中，臺站會根據(jù)標(biāo)定結(jié)果對地震計電壓靈敏度進(jìn)行相應(yīng)的修正[2]，即地震計電壓靈敏度并非一定是等于默認(rèn)值的。因此，若無法通過直接獲取準(zhǔn)確的前放增益和地震計電壓靈敏度計算S 值時，該怎么換算地動速度呢？文中提出了一個快速的方法，不需要繁瑣的查詢步驟便可以獲得準(zhǔn)確的S 值。

以松潘地震臺CTS-1 型地震計記錄的2019 年1 月15 日08 時0～120 秒的垂直分向數(shù)據(jù)為例。首先用MSDP 軟件打開數(shù)據(jù)，在數(shù)據(jù)中隨機選擇一個時間點（文中選擇的是零點）并記錄下該時間點的速度值（零點處速度值為-2.969μm/s）、振幅值（零點處振幅值為-3737.0 counts）。

依據(jù)得到的速度值XV與振幅值A(chǔ)V0，再將已知的A/D 轉(zhuǎn)換因子Kad代入公式，通過計算得到零點處S 值為1500.33816。然而這只是零點處的S 值，并不代表真實值。因此，作者通過對300 個時間點進(jìn)行抽樣采樣與計算，分析得到S 值等于1500+X，其中X∈(0,0.6)。將S 值取整為1500,代入公式計算地動速度，再將計算結(jié)果與原速度值進(jìn)行對比。其中，計算得到的速度值相較原速度值的最大誤差率為0.04%。對比結(jié)果表明：S 值取整所造成的誤差是可以忽略的，S 值取整后即為1500。這與其它分向中得到的結(jié)論是一致的。確定下S的值，那么就可應(yīng)用于整段波形的數(shù)據(jù)轉(zhuǎn)換，于是將整段0～120 秒波形共計12000 個時間點代入公式計算，求出速度值的序列并繪出波形。如圖1 所示，數(shù)據(jù)轉(zhuǎn)換后的波形形態(tài)與原波形是保持一致的。

圖1 松潘臺數(shù)據(jù)轉(zhuǎn)換后的波形

3 影響S 值的因素

在同一觀測系統(tǒng)下，采用文中方法計算并取整后的S 值都是穩(wěn)定的，如表1 所示。S 值的形式為：S=A+X,其中A 為正整數(shù)，X 為小數(shù)。通常情況下，A的值在個位數(shù)一般存在±1 左右的浮動，X 值的浮動略大，但經(jīng)過四舍五入取整后得到的S 值是一定的。值得注意的是，若計算時所選擇時間點的振幅值A(chǔ)v0過小會造成S 值出現(xiàn)明顯差異。

表1 松潘臺2018 年1-3 月不同時間點S 值計算結(jié)果

以延邊地震臺2016 年5 月13 日00 時垂直分向數(shù)據(jù)作為示例，已知延邊地震臺A/D 轉(zhuǎn)換因子Kad等于1.589μV/counts。選擇5 個時間點分別計算S 值并轉(zhuǎn)換為地動速度，結(jié)果如表2 所示。當(dāng)選擇時間點的振幅值A(chǔ)V0過小，則計算得到的S 值會不準(zhǔn)確，若以該S 值進(jìn)行數(shù)據(jù)轉(zhuǎn)換會造成速度的誤差變大。這種速度的誤差放在振幅值A(chǔ)v0較小的時間點并不明顯，但放在振幅值A(chǔ)V0較大的時間點則非常突出。當(dāng)以較大振幅值A(chǔ)V0計算得到的S 值進(jìn)行速度轉(zhuǎn)換時，對于振幅值A(chǔ)V0較小和較大的時間點均不會造成明顯的速度誤差。因此實際的應(yīng)用中，應(yīng)選擇振幅值A(chǔ)V0較大（振幅值應(yīng)在三位數(shù)以上）的時間點計算S的值，并根據(jù)計算結(jié)果合理選擇取整的整數(shù)。

表2 延邊臺不同S 值的速度轉(zhuǎn)換結(jié)果

4 實例應(yīng)用

為了證明方法的適用性，文中選擇武漢地震臺2016 年5 月13 日16 時31 分15 秒至34 分15 秒垂直分向的一段波形數(shù)據(jù)，共計時間長度為180 秒，數(shù)據(jù)序列長度為18000。已知武漢臺的A/D 轉(zhuǎn)換因子Kad等于1.589μV/counts。首先打開MSDP 軟件獲取任意5 個時間點的速度值Xv與振幅值A(chǔ)v并計算出這5個時間點的S的值，得到S的值為1969+y,其中y∈(-0.5,0.3)。將S 值進(jìn)行取整為1969，代入公式計算地動速度，再將計算結(jié)果與原速度值進(jìn)行對比。對比的結(jié)果如表3 所示，計算得到的速度值相較原速度值的最大誤差率僅為0.02%，說明S的值取整所造成的誤差是可以忽略的。將整段數(shù)據(jù)代入公式，計算速度序列并繪圖，如圖2 所示，數(shù)據(jù)轉(zhuǎn)換后的波形形態(tài)與原波形同樣是保持一致的。

表3 武漢臺地動速度及S 值計算結(jié)果

圖2 武漢臺數(shù)據(jù)轉(zhuǎn)換后的波形

5 結(jié)論

通過對測震無量綱整數(shù)轉(zhuǎn)換為地動速度的快速方法的研究，從中得到了下面幾點認(rèn)識：

5.1 若無法直接獲取前放增益和地震計電壓靈敏度的參數(shù)時，采用本文中的方法可以快速、準(zhǔn)確地實現(xiàn)數(shù)據(jù)的轉(zhuǎn)換。

5.2 在實際應(yīng)用中，應(yīng)選擇至少5 個時間點計算S 值，并根據(jù)計算結(jié)果合理選擇取整的整數(shù)。S 值的表現(xiàn)形式為：S=A+X,其中A 為正整數(shù)，X 為小數(shù)。通常X的值浮動較大，A的值通常在個位數(shù)存在±1 左右的浮動，但四舍五入取整后得到的S 值是一定的。

5.3 當(dāng)選擇時間點的振幅值A(chǔ)V0過小（振幅值為兩位數(shù)），則計算得到的S 值不準(zhǔn)確，進(jìn)而造成明顯的速度誤差，因此，應(yīng)選擇振幅值A(chǔ)V0較大（振幅值應(yīng)在三位數(shù)以上）的時間點進(jìn)行S 值的計算。