陳力剛 張曉宇 邢 康 賈漯昭

（河南省地震局, 鄭州 450018）

引言

國家地震烈度速報(bào)與預(yù)警工程項(xiàng)目（簡(jiǎn)稱“預(yù)警項(xiàng)目”）是一項(xiàng)具有防震減災(zāi)實(shí)效的民生工程和社會(huì)工程，建成后將有效提高減災(zāi)和社會(huì)服務(wù)能力。河南省預(yù)警項(xiàng)目中一般站烈度計(jì)設(shè)備于2021 年上半年安裝完畢，相關(guān)記錄數(shù)據(jù)可通過網(wǎng)絡(luò)遠(yuǎn)程訪問下載，為初步計(jì)算一般站的背景噪聲創(chuàng)造了條件。預(yù)警項(xiàng)目針對(duì)臺(tái)站數(shù)據(jù)包格式做出了新定義，各廠家生產(chǎn)的硬件采集設(shè)備已按照新的格式規(guī)程生產(chǎn)數(shù)據(jù)，而現(xiàn)有計(jì)算軟件2北京港震機(jī)電有限責(zé)任公司，童汪練，2017. 數(shù)字地震儀參數(shù)測(cè)定軟件使用手冊(cè)，北京港震機(jī)電有限責(zé)任公司資料.附帶的數(shù)據(jù)分道軟件無法直接打開預(yù)警數(shù)據(jù)采集器產(chǎn)出的數(shù)據(jù)，因此，本文通過對(duì)預(yù)警項(xiàng)目臺(tái)站數(shù)據(jù)包格式規(guī)程的解讀，采用C 語言開發(fā)了預(yù)警項(xiàng)目臺(tái)站波形數(shù)據(jù)解壓軟件，產(chǎn)出相應(yīng)的分道文本數(shù)據(jù)文件，并通過具體實(shí)例對(duì)預(yù)警項(xiàng)目連續(xù)波形數(shù)據(jù)包格式和相應(yīng)的程序?qū)崿F(xiàn)給出詳細(xì)說明。

1 預(yù)警項(xiàng)目臺(tái)站波形數(shù)據(jù)包組成

1.1 波形數(shù)據(jù)包格式

預(yù)警項(xiàng)目臺(tái)站數(shù)據(jù)包格式規(guī)程中定義了多種數(shù)據(jù)包格式，其標(biāo)識(shí)有“cc”、“si”、“ti”、“wc”、“wt”和“ws”，本文著重研究“wc”數(shù)據(jù)包，即連續(xù)波形數(shù)據(jù)包?！皐c”數(shù)據(jù)包以miniSEED 數(shù)據(jù)卷為基礎(chǔ)，做了必要的修訂，主要變化為：① miniSEED 數(shù)據(jù)卷固定頭段區(qū)的第0 至第1 個(gè)字節(jié)定義為“wc”數(shù)據(jù)包標(biāo)識(shí)位置；② miniSEED 數(shù)據(jù)卷固定頭段區(qū)的第2 個(gè)字節(jié)內(nèi)高5 位bit 到第5 個(gè)字節(jié)定義為數(shù)據(jù)包順序編號(hào)位置，用于數(shù)據(jù)傳輸時(shí)斷點(diǎn)續(xù)傳，而第2 個(gè)字節(jié)內(nèi)低3 位bit 定義為數(shù)據(jù)包長(zhǎng)度索引；③ miniSEED 數(shù)據(jù)卷中“純數(shù)據(jù)子塊[ 1000]”后增加了8 個(gè)字節(jié)的數(shù)據(jù)標(biāo)識(shí)區(qū)，用于標(biāo)識(shí)通道信息順序、數(shù)據(jù)量綱和靈敏度信息；④其他少許的新定義，如第36 個(gè)字節(jié)第7 bit 定義為入網(wǎng)標(biāo)記等。由此發(fā)現(xiàn)，預(yù)警項(xiàng)目對(duì)mini-SEED 數(shù)據(jù)卷僅做了少量修訂。為提取“wc”類連續(xù)波形數(shù)據(jù)具體內(nèi)容，需對(duì)該數(shù)據(jù)卷格式信息做出解讀，在此僅對(duì)提取波形使用到的信息做出說明。miniSEED 被定義為SEED 格式標(biāo)準(zhǔn)中的數(shù)據(jù)記錄格式，一個(gè)完整的預(yù)警項(xiàng)目miniSEED 數(shù)據(jù)卷包括固定頭段區(qū)、可變頭段區(qū)和數(shù)據(jù)區(qū)（中國地震局，2003）。

固定頭段區(qū)包含使用數(shù)據(jù)所需的最低限度信息，主要信息按字節(jié)從前向后依次為波形類別標(biāo)識(shí)、數(shù)據(jù)包序號(hào)、質(zhì)量指示符、臺(tái)站標(biāo)識(shí)、位置標(biāo)識(shí)、通道標(biāo)識(shí)、臺(tái)網(wǎng)標(biāo)識(shí)、記錄起始時(shí)間、樣本數(shù)目、采樣率、子塊等。前20 個(gè)字節(jié)除第2～5 個(gè)字節(jié)外，均采用ASCII 碼表示，后28 個(gè)字節(jié)采用二進(jìn)制表示（張旸等，2007），固定頭段區(qū)共48 個(gè)字節(jié)。

可變頭段區(qū)由多個(gè)可選的子塊構(gòu)成，在預(yù)警項(xiàng)目中選擇的子塊類型僅為“純數(shù)據(jù)子塊 [ 1000]”。該子塊包含了對(duì)數(shù)據(jù)區(qū)解碼的必要信息，主要信息有編碼格式代碼，一般為0 B，代表Steim2 壓縮格式；字序代碼一般為1，代表大端字節(jié)序，高位字節(jié)在前，低位字節(jié)在后；標(biāo)識(shí)數(shù)據(jù)記錄長(zhǎng)度的指數(shù)代碼如果為8，則數(shù)據(jù)記錄長(zhǎng)度為28字節(jié)。純數(shù)據(jù)子塊 [ 1000] 為8 個(gè)字節(jié)，均采用二進(jìn)制表示。

純數(shù)據(jù)子塊 [ 1000] 區(qū)域后的8 個(gè)字節(jié)是最具有預(yù)警項(xiàng)目特色定義的區(qū)域，位于可變頭段區(qū)內(nèi)，因此在此區(qū)域增加自定義信息最適合。本塊信息主要涉及有關(guān)數(shù)據(jù)標(biāo)識(shí)的內(nèi)容，包括數(shù)據(jù)區(qū)數(shù)據(jù)的通道順序，如“ENZ”表示數(shù)據(jù)按照“**E”通道數(shù)據(jù)、“**N”通道數(shù)據(jù)、“**Z”通道數(shù)據(jù)的順序排列；數(shù)據(jù)量綱取值范圍為0～3，0 表示無量綱，1 表示位移，2 表示速度，3 表示加速度；靈敏度因子記錄儀器記錄的靈敏度因子取值范圍為0～0x3FFFFFFF。預(yù)警項(xiàng)目定義的數(shù)據(jù)標(biāo)識(shí)內(nèi)容多數(shù)為二進(jìn)制，實(shí)際上這8 個(gè)字節(jié)內(nèi)容可視為1 個(gè)變相的子塊放在可變頭段區(qū)內(nèi)，以增加對(duì)數(shù)據(jù)區(qū)的說明。

固定頭段區(qū)和可變頭段區(qū)數(shù)據(jù)包含了臺(tái)網(wǎng)、臺(tái)站、通道、采樣率、樣本數(shù)目、時(shí)間起點(diǎn)、編碼格式、字節(jié)序、數(shù)據(jù)記錄長(zhǎng)度等數(shù)據(jù)描述信息（王曉磊等，2016），占據(jù)了每條miniSEED 數(shù)據(jù)卷的前64 個(gè)字節(jié)區(qū)域，這些信息是波形數(shù)據(jù)區(qū)的輔助信息（楊周勝等，2019），也是編程提取和使用的主要參數(shù)。

數(shù)據(jù)區(qū)包含實(shí)際的時(shí)間序列數(shù)據(jù)（中國地震局，2003），數(shù)據(jù)區(qū)位置從固定頭段區(qū)第44 至第45 個(gè)字節(jié)中內(nèi)容指定的字節(jié)序數(shù)處開始，這部分?jǐn)?shù)據(jù)塊一般采用Steim2 算法壓縮（何加勇等，2009），由若干數(shù)據(jù)幀組成，每個(gè)數(shù)據(jù)幀為64 個(gè)字節(jié)。波形數(shù)據(jù)區(qū)具體大小及數(shù)據(jù)通道、采樣率、樣本數(shù)目、量綱、字節(jié)序等信息均已在上述固定頭段區(qū)和可變頭段區(qū)中定義。

1.2 Steim2 壓縮算法

預(yù)警項(xiàng)目中的波形數(shù)據(jù)區(qū)數(shù)據(jù)為Steim2 格式，獲取觀測(cè)數(shù)據(jù)后，首先需獲取打開Steim2 格式數(shù)據(jù)壓縮的鑰匙才能使用。Steim2 壓縮算法是借用語音壓縮中最簡(jiǎn)單的差分壓縮算法發(fā)展而成的（王翠芳等，2011），其特點(diǎn)是無損壓縮，目前國內(nèi)外地震數(shù)據(jù)傳輸和本地存儲(chǔ)普遍采用了壓縮算法。地震數(shù)據(jù)采用無損壓縮算法的基本原理為：在正常地震背景情況下，地震計(jì)輸出的時(shí)間序列數(shù)據(jù)通常具有較強(qiáng)的關(guān)聯(lián)性，即前一個(gè)采樣點(diǎn)相比后一個(gè)采樣點(diǎn)變化小，差值小，該差值占用計(jì)算機(jī)內(nèi)存小，可能遠(yuǎn)遠(yuǎn)小于32 bit，而保存1 個(gè)原始值卻需32 bit，因此若多個(gè)差值按規(guī)則填充滿32 bit，可極大地節(jié)省存儲(chǔ)空間，這種方式可理解為數(shù)據(jù)壓縮。后期為恢復(fù)原始值，須在特定位置保存一個(gè)前向積分常數(shù)。筆者認(rèn)為Steim2 數(shù)據(jù)格式在壓縮時(shí)具有高保真、高效的特征，符合地震動(dòng)波形記錄“長(zhǎng)期平靜”和“長(zhǎng)期無震”的特點(diǎn)（王翠芳等，2011），但此算法在地震事件發(fā)生時(shí)，壓縮效率顯著降低，這是無損壓縮算法的特點(diǎn)。

根據(jù)Steim2 格式數(shù)據(jù)壓縮基本思想，合成壓縮數(shù)據(jù)時(shí)，首先利用直接觀測(cè)數(shù)據(jù)計(jì)算差值序列，然后按照規(guī)則利用差值序列整合重建32 bit 數(shù)據(jù)序列。獲取數(shù)據(jù)后，用戶更關(guān)心數(shù)據(jù)解壓過程，其步驟與壓縮相反，具體如下：

（1）按規(guī)則將32 bit 數(shù)據(jù)序列拆分為差值序列

數(shù)據(jù)拆分的規(guī)則鑰匙隨著壓縮數(shù)據(jù)進(jìn)行傳輸，壓縮數(shù)據(jù)幀的結(jié)構(gòu)如圖1 所示。

改革開放尤其是十八大以來,黨中央提出了關(guān)于生態(tài)文明建設(shè)的一系列新理念新要求。在生態(tài)文明理念方面,明確提出了“尊重自然、順應(yīng)自然、保護(hù)自然”“綠水青山就是金山銀山”“山水林田湖草是一個(gè)生命共同體”等理念。在生態(tài)文明與經(jīng)濟(jì)社會(huì)發(fā)展的關(guān)系方面,提出了“守住生態(tài)與發(fā)展兩條底線”“良好生態(tài)環(huán)境是最公平的公共產(chǎn)品,是最普惠的民生福祉”等理念。在生態(tài)文明實(shí)現(xiàn)路徑方面,強(qiáng)調(diào)要樹立“劃定紅線,守住底線與資源上限”“像保護(hù)眼睛一樣保護(hù)生態(tài)環(huán)境,像對(duì)待生命一樣對(duì)待生態(tài)環(huán)境”等理念。

圖1 中每個(gè)字母長(zhǎng)度均為32 位，解壓的鑰匙保存在每行行首W0和其后W1～W15中最高的兩位中，拆分對(duì)象是W1～W15中后30 位二進(jìn)制數(shù)據(jù)。鑰匙（Ck）須與鎖（Wk）一一對(duì)應(yīng)，才能順利解鎖，如圖2 所示。

圖1 數(shù)據(jù)幀結(jié)構(gòu)Fig. 1 Data frame structure

圖2 Ck 與Wk 的關(guān)系及dnib 位置Fig. 2 Relationship between Ck and Wk and dnib location

Ck為2 bit，取值為00、01、10、11，若壓縮數(shù)據(jù)為Steim1 格式，這4 類鑰匙可拆分Wk數(shù)據(jù)，但預(yù)警項(xiàng)目中規(guī)定使用Steim2 格式壓縮數(shù)據(jù)，需另外4 類鑰匙才可打開第2 道門，從而拆分Wk數(shù)據(jù)，而打開第2道門的鑰匙隱藏在Wk的最高2 位中。將Ck稱為“廣義字節(jié)編碼”，Wk的高端2 位（dnib）稱為“細(xì)解碼”，其組合使用如下：

根據(jù)上述規(guī)則，按照Ck和dnib 組合定義的拆分方法，可將Wk從最低位開始按照4、5、6、8、10、15、30 位的方式進(jìn)行切割，分離出的數(shù)據(jù)即為差值。

（2）將差值序列還原為原始數(shù)據(jù)

獲取差值后，通過前向積分常數(shù)還原數(shù)據(jù)，前向積分常數(shù)保存在數(shù)據(jù)幀0 中的X0內(nèi)。依次算出原始數(shù)據(jù)，計(jì)算步驟如下：

在數(shù)據(jù)幀0 中的Xn內(nèi)還保存有這組數(shù)據(jù)的反向積分常數(shù)，利用該反向積分常數(shù)與最后計(jì)算出的樣本原值進(jìn)行比較。如果二者不同，說明數(shù)據(jù)出現(xiàn)了錯(cuò)誤。

2 數(shù)據(jù)包解析實(shí)例

以預(yù)警項(xiàng)目一般站MI3000-C 烈度計(jì)數(shù)據(jù)采集器產(chǎn)出數(shù)據(jù)為例，對(duì)連續(xù)波形包數(shù)據(jù)進(jìn)行解析，給出的實(shí)際二進(jìn)制（用十六進(jìn)制表示）數(shù)據(jù)如圖3 所示。

圖3 MI3000-C 頭段和數(shù)據(jù)幀0 說明Fig. 3 MI3000-C Header and Data frame 0 description

從固定頭段區(qū)和可變頭段區(qū)可提取諸多參數(shù)信息，如可變頭段區(qū)的“C0”字節(jié)低2 bit 中包含數(shù)字0，說明數(shù)據(jù)單位為無量綱，此處解壓后的數(shù)據(jù)單位為count。如果“C0”字節(jié)低2 bit 中包含數(shù)字為1 或2，則必須給出確定的靈敏度因子，可直接輸出物理量值，但目前獲取的數(shù)據(jù)為count 值，所以需要在后端服務(wù)器軟件進(jìn)行實(shí)際物理量轉(zhuǎn)換。在本實(shí)例中，可獲取前向積分常數(shù)X0（FFFFFFA9）為十進(jìn)制?87，反向積分常數(shù)Xn（FFFFFF7A）為十進(jìn)制?134。隨機(jī)抽取W6為樣例，W0和W6數(shù)據(jù)解析及新合成的W6數(shù)據(jù)如圖4 所示。

圖4 W0 和W6 數(shù)據(jù)解析及新合成的W6 數(shù)據(jù)表Fig. 4 W0 and W6 data analysis and newly synthesized W6 data

利用上述分析結(jié)果編寫代碼后才能實(shí)用化，對(duì)數(shù)據(jù)進(jìn)行逐個(gè)解析，首先確定各數(shù)據(jù)的定位，需設(shè)置一個(gè)參考起點(diǎn)，根據(jù)第44 至第45 個(gè)字節(jié)的內(nèi)容獲取數(shù)據(jù)區(qū)的開始偏移量并進(jìn)行定位。筆者通過第15 至第17個(gè)字節(jié)的通道標(biāo)識(shí)符進(jìn)行定位，其目的是按通道解壓數(shù)據(jù)，并形成相應(yīng)通道的文本文件，以配套CAL79_噪聲功率譜計(jì)算軟件使用。

程序流程如圖5 所示，分為數(shù)據(jù)讀取、字節(jié)序轉(zhuǎn)換、解壓、校驗(yàn)、文本文件形成等部分，此處著重對(duì)數(shù)據(jù)拆分的實(shí)現(xiàn)過程進(jìn)行說明。

圖5 程序流程Fig. 5 Program flow chart

由程序代碼可知較關(guān)鍵的步驟為數(shù)據(jù)定位、以數(shù)據(jù)幀為單位執(zhí)行循環(huán)、采用“按位與”的方法提取信息、對(duì)出現(xiàn)負(fù)值的情況給予處理。筆者通過大量二進(jìn)制的位操作提高執(zhí)行效率。

經(jīng)實(shí)際測(cè)試，該軟件可對(duì)河南省一般站MI3000-C、REMOS-SIT4 烈度計(jì)產(chǎn)出的小時(shí)波形文件進(jìn)行分道解碼，且可對(duì)基準(zhǔn)站和基本站使用的EDAS-24 GN（EEW）、TDE-324 FI/CI 數(shù)據(jù)采集器產(chǎn)出的小時(shí)波形文件進(jìn)行分道解碼，并在程序解碼過程中設(shè)置了數(shù)據(jù)驗(yàn)證功能，即前向積分常數(shù)X0加各項(xiàng)差值等于反向積分常數(shù)Xn，程序判斷為真才能輸出結(jié)果，增加了數(shù)據(jù)解碼的準(zhǔn)確性。在獲取解壓后的分道文本數(shù)據(jù)后，可順利使用噪聲功率譜計(jì)算軟件進(jìn)行臺(tái)站背景噪聲環(huán)境評(píng)估。

3 結(jié)論

本文對(duì)預(yù)警項(xiàng)目臺(tái)站波形數(shù)據(jù)包格式進(jìn)行了研究，得出以下結(jié)論：

（1）預(yù)警項(xiàng)目對(duì)數(shù)據(jù)采集器的產(chǎn)出波形數(shù)據(jù)格式做出了統(tǒng)一規(guī)定，對(duì)數(shù)據(jù)的部分字節(jié)內(nèi)容進(jìn)行了重新定義，準(zhǔn)確把握各字節(jié)位的定義和關(guān)系是數(shù)據(jù)解析與使用的前提。

（2）連續(xù)波形數(shù)據(jù)包結(jié)構(gòu)中固定頭段區(qū)、可變頭段區(qū)和數(shù)據(jù)區(qū)是有機(jī)整體，解壓算法需相關(guān)區(qū)域的參數(shù)參與計(jì)算，解壓的關(guān)鍵在于二進(jìn)制數(shù)據(jù)按位拆分。

（3）列出了實(shí)用的程序代碼，利用C 語言的位操作方法提取有效信息，并拆分?jǐn)?shù)據(jù)。筆者查閱了部分文獻(xiàn)，發(fā)現(xiàn)針對(duì)地震數(shù)據(jù)包結(jié)構(gòu)和壓縮的原理描述較多，且較抽象，可操作性不強(qiáng)，具體實(shí)現(xiàn)過程多用一個(gè)函數(shù)包一筆帶過，對(duì)實(shí)現(xiàn)細(xì)節(jié)的解讀較少，本文給出了具體實(shí)現(xiàn)方法，稍加修改即可運(yùn)行。