張偉,胡蕾,張釗搏

(1.中國(guó)地質(zhì)調(diào)查局 成都地質(zhì)調(diào)查中心，四川 成都 610082; 2.電子科技大學(xué) 資源與環(huán)境學(xué)院，四川 成都 611731; 3.中科院成都信息技術(shù)股份有限公司，四川 成都 610041)

0 引言

長(zhǎng)周期大地電磁法目前是探測(cè)地球深部(>150 km)上地幔尺度電性結(jié)構(gòu)特征的主要方法，由烏克蘭Lviv Centre of Institute of Space Research研究所研制的LEMI-417儀器是最早引入到我國(guó)用于長(zhǎng)時(shí)間、穩(wěn)定、連續(xù)觀測(cè)超低頻大地電磁場(chǎng)的一套成熟裝備。2007年起，中國(guó)科學(xué)院地質(zhì)與地球物理研究所[1]、成都理工大學(xué)[2]、中國(guó)地質(zhì)大學(xué)(北京)[3-4]以及中國(guó)地震局地質(zhì)研究所[5]等科研院所開始使用該儀器探測(cè)青藏高原及鄰區(qū)、龍門山構(gòu)造帶以及川滇黔相鄰區(qū)地震活動(dòng)帶的深部地質(zhì)結(jié)構(gòu)，用于揭示其關(guān)鍵地質(zhì)帶的深部動(dòng)力學(xué)特征。與該儀器配套的數(shù)據(jù)處理軟件是由俄羅斯Ivan M. Varantsov等編寫的PRC_MTMV商業(yè)軟件，用于實(shí)現(xiàn)將大地電磁場(chǎng)從時(shí)間域到頻率域的轉(zhuǎn)換，通過傅里葉變換計(jì)算得到功率譜，并通過功率譜估計(jì)得到大地電磁阻抗要素。

眾所周知，計(jì)算機(jī)對(duì)于二進(jìn)制格式數(shù)據(jù)的讀取效率遠(yuǎn)遠(yuǎn)大于文本格式數(shù)據(jù)，因此在進(jìn)行傅里葉變換時(shí)直接讀取二進(jìn)制數(shù)據(jù)進(jìn)行計(jì)算將大大提高程序的運(yùn)行效率。但由于烏克蘭儀器生產(chǎn)廠家沒有公開給出其原始數(shù)據(jù)的二進(jìn)制格式，因此在處理LEMI-417數(shù)據(jù)時(shí)，必須首先使用廠家提供的轉(zhuǎn)換程序?qū)x器二進(jìn)制記錄數(shù)據(jù)轉(zhuǎn)換為文本格式，再使用PRC_MTMV軟件的LEMI2PRC子模塊將ASCII文本格式數(shù)據(jù)再次轉(zhuǎn)換為二進(jìn)制數(shù)據(jù)后方能進(jìn)行后續(xù)的傅里葉變換處理[3]，這無疑明顯地增加了一步不必要的處理操作。為此，作者對(duì)照V5-2000儀器的二進(jìn)制數(shù)據(jù)格式的特點(diǎn)[6]，對(duì)LEMI-417儀器的原始二進(jìn)制數(shù)據(jù)格式進(jìn)行了深入研究，在不斷試錯(cuò)后，成功剖析出了該儀器二進(jìn)制記錄數(shù)據(jù)的數(shù)據(jù)結(jié)構(gòu)。

1 LEMI-417型系統(tǒng)的原始數(shù)據(jù)格式

LEMI-417是專門針對(duì)長(zhǎng)時(shí)間采集超低頻段天然電磁場(chǎng)而研制的一款觀測(cè)系統(tǒng)，其主要特點(diǎn)是：①采用磁通門傳感器來采集磁場(chǎng)數(shù)據(jù)，觀測(cè)范圍為±70 000 nT，可以記錄到幾十萬秒的長(zhǎng)周期信號(hào)；②電道采集盒子通過專門的電路設(shè)計(jì)能夠保證在長(zhǎng)時(shí)間范圍內(nèi)保持很低的漂移，提供了4組電道接口，可以同時(shí)采集4組電極的電場(chǎng)數(shù)據(jù)；③采用高精度GPS同步授時(shí)能夠保證在長(zhǎng)時(shí)間范圍內(nèi)具有很小的時(shí)移；④整體功耗很低，能夠保證在1塊12V45ah的鉛酸電池供電下連續(xù)工作一周以上。

該儀器以1 Hz的采樣率來連續(xù)觀測(cè)天然電磁場(chǎng)，采集的數(shù)據(jù)以二進(jìn)制格式實(shí)時(shí)儲(chǔ)存到儀器內(nèi)部的CF卡中，文件名以L417_×××.B××的規(guī)則來進(jìn)行命名，其中：L417為儀器的型號(hào)名稱；_×××為單點(diǎn)的數(shù)據(jù)存放文件序號(hào)，第一個(gè)數(shù)據(jù)文件編號(hào)為001，當(dāng)GPS時(shí)間為格林威治時(shí)間的零點(diǎn)時(shí)，儀器系統(tǒng)自動(dòng)新建一個(gè)新文件來存放數(shù)據(jù)，文件序號(hào)自動(dòng)在前一個(gè)文件序號(hào)的基礎(chǔ)上累加1；B××的首字母B是二進(jìn)制(binary)的縮寫，××是采集儀器的出廠編號(hào)。

1.1 頭文件格式

通常與測(cè)點(diǎn)數(shù)據(jù)采集相關(guān)的記錄信息(如儀器型號(hào)、采集日期、經(jīng)緯度、高程、電極距等)一般以文件頭的形式存放于二進(jìn)制文件的首端。LEMI-417I型儀器共分配了32個(gè)字節(jié)的大小空間來存放信息(表1)，各字節(jié)位置信息闡述如下。

表1 LEMI-417原始二進(jìn)制記錄數(shù)據(jù)的頭文件格式(tag)

1～4字節(jié)：4個(gè)字節(jié)長(zhǎng)度，存放采集儀器的型號(hào)信息，對(duì)于LEMI-417I型儀器來說，第1個(gè)字節(jié)為ASCII字符‘4’(16進(jìn)制為34)，第2字節(jié)為ASCII字符‘1’(16進(jìn)制為31)，第3字節(jié)為ASCII字符‘7’(16進(jìn)制為37)，第4字節(jié)為ASCII字符‘I’(16進(jìn)制為49)。

5字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放采集儀器的出廠編號(hào)。

6字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放開始采集數(shù)據(jù)時(shí)的當(dāng)前年信息。

7字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放開始采集數(shù)據(jù)時(shí)的當(dāng)前月信息。

8字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放開始采集數(shù)據(jù)時(shí)的當(dāng)前日信息。

9字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放開始采集數(shù)據(jù)時(shí)的當(dāng)前小時(shí)信息(24小時(shí)制)。

10字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放開始采集數(shù)據(jù)時(shí)的當(dāng)前分信息。

11字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放開始采集數(shù)據(jù)時(shí)的當(dāng)前秒信息。

12～15字節(jié)：4個(gè)字節(jié)長(zhǎng)度，存放當(dāng)前測(cè)點(diǎn)GPS同步的緯度。

16字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放當(dāng)前測(cè)點(diǎn)坐標(biāo)是位于北緯還是南緯，以我國(guó)位于北半球?yàn)槔?，在我?guó)測(cè)量數(shù)據(jù)時(shí)應(yīng)為ASCII字符‘N’( 16進(jìn)制為4E)。

17～21字節(jié)：5個(gè)字節(jié)長(zhǎng)度，存放當(dāng)前測(cè)點(diǎn)GPS同步記錄的經(jīng)度信息。

22字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放當(dāng)前測(cè)點(diǎn)坐標(biāo)是位于東經(jīng)還是西經(jīng)，以我國(guó)位于東半球?yàn)槔?，在我?guó)測(cè)量數(shù)據(jù)時(shí)應(yīng)為ASCII字符‘E’( 16進(jìn)制為45)。

23字節(jié)：1個(gè)字節(jié)長(zhǎng)度，為空字符(對(duì)應(yīng)的16進(jìn)制為00)。

24～25字節(jié)：2個(gè)字節(jié)長(zhǎng)度，存放當(dāng)前測(cè)點(diǎn)GPS同步的高程(單位：m)。

26字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放儀器的采樣頻率。

27字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放儀器內(nèi)部的VIN(外部供電電源電壓)信息。

28字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放儀器內(nèi)部的VBAT(內(nèi)部紐扣電池電壓)信息。

29字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放電道盒子上第1對(duì)電極的極距長(zhǎng)度(單位：m)。

30字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放電道盒子上第2對(duì)電極的極距長(zhǎng)度(單位：m)。

31字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放電道盒子上第3對(duì)電極的極距長(zhǎng)度(單位：m)。

32字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放電道盒子上第4對(duì)電極的極距長(zhǎng)度(單位：m)。

1.2 時(shí)間序列格式

如圖1所示，在頭文件(Tag)后緊接著的是時(shí)間序列數(shù)據(jù)，LEMI-417以30個(gè)字節(jié)的大小空間來存放一個(gè)采樣時(shí)間片采集到的電磁場(chǎng)數(shù)據(jù)段(segment，表2所示)，由17個(gè)時(shí)間片組成一個(gè)記錄塊(record)，一個(gè)數(shù)據(jù)段結(jié)束后再接下一個(gè)文件頭數(shù)據(jù)，依次循環(huán)、直到采集結(jié)束。其中，單個(gè)時(shí)間片數(shù)據(jù)段中各字節(jié)位置信息闡述如下。

圖1 LEMI-417原始數(shù)據(jù)結(jié)構(gòu)分解示意Fig.1 Schematic diagram of original data structure decomposition

表2 LEMI-417原始二進(jìn)制記錄數(shù)據(jù)的時(shí)間序列格式(segment)

1～3字節(jié)：3個(gè)字節(jié)長(zhǎng)度，存放磁通門傳感器X方向的測(cè)量值(單位：nT)。

4～6字節(jié)：3個(gè)字節(jié)長(zhǎng)度，存放磁通門傳感器Y方向的測(cè)量值(單位：nT)。

7～9字節(jié)：3個(gè)字節(jié)長(zhǎng)度，存放磁通門傳感器Z方向的測(cè)量值(單位：nT)。

10～13字節(jié)：4個(gè)字節(jié)長(zhǎng)度，存放第1道電極間的測(cè)量值(單位：mkV/m)。

14～17字節(jié)：4個(gè)字節(jié)長(zhǎng)度，存放第2道電極間的測(cè)量值(單位：mkV/m)。

18～21字節(jié)：4個(gè)字節(jié)長(zhǎng)度，存放第3道電極間的測(cè)量值(單位：mkV/m)。

22～25字節(jié)：4個(gè)字節(jié)長(zhǎng)度，存放第4道電極間的測(cè)量值(單位：mkV/m)。

26～27字節(jié)：2個(gè)字節(jié)長(zhǎng)度，存放磁通門傳感器溫度(單位：℃)。

28～29字節(jié)：2個(gè)字節(jié)長(zhǎng)度，存放電道采集盒子的溫度(單位：℃)。

30字節(jié)：1個(gè)字節(jié)長(zhǎng)度，存放時(shí)間片采樣間隔(單位：m)。

2 計(jì)算機(jī)數(shù)據(jù)類型與程序?qū)崿F(xiàn)

2.1 頭文件數(shù)據(jù)結(jié)構(gòu)與核心代碼

在C++語言中可利用struct類型變量來結(jié)構(gòu)化地讀取和存儲(chǔ)頭文件信息，再按表1所示各信息的字節(jié)位置和字節(jié)長(zhǎng)度，定義其內(nèi)部的子變量數(shù)據(jù)類型，其核心代碼如下：

// 頭文件數(shù)據(jù)結(jié)構(gòu)定義

structTagData

{

unsignedcharstr_1; //4

unsignedcharstr_2; //1

unsignedcharstr_3; //7

unsignedcharstr_4; //儀器型號(hào)

unsignedcharno; //儀器編號(hào)

unsignedcharyear; //年

unsignedcharmonth; //月

unsignedcharday; //日

unsignedcharhour; //時(shí)

unsignedcharminute; //分

unsignedcharsecond; //秒

unsignedcharlatitude[4]; //緯度

charN_S; //北緯還是南緯

unsignedcharlongitude[5]; //經(jīng)度

charE_W; //東經(jīng)還是西經(jīng)

charseparator; //空字符

unsignedcharaltitude[2]; //高程

unsignedcharaverage; //采樣率

unsignedcharvin; // 外部電壓

unsignedcharvbat; // 內(nèi)部電壓

unsignedcharL1; //第一對(duì)極距長(zhǎng)度

unsignedcharL2; //第二對(duì)極距長(zhǎng)度

unsignedcharL3; //第三對(duì)極距長(zhǎng)度

unsignedcharL4; //第四對(duì)極距長(zhǎng)度

}

// … …由于篇幅有限省略部分非關(guān)鍵代碼

// 讀取儀器原始數(shù)據(jù)文件

TagDatatag;

FILE*fp=fopen(strFilePath,"rb");

chardataBuf[32] = "";

fread(dataBuf,sizeof(dataBuf),1,fp);

fclose(fp);

// 結(jié)構(gòu)體各子變量賦值

tag.str_1=dataBuf[0];

tag.str_2=dataBuf[1];

// … …由于篇幅有限省略部分代碼

tag.latitude[0]=dataBuf[11];

tag.latitude[1]=dataBuf[12];

tag.latitude[2]=dataBuf[13];

tag.latitude[3]=dataBuf[14];

// … …由于篇幅有限省略部分代碼

tag.L4=dataBuf[31];

2.2 時(shí)間序列數(shù)據(jù)結(jié)構(gòu)與核心代碼

同理，利用struct類型變量來結(jié)構(gòu)化地讀取和存儲(chǔ)時(shí)間片數(shù)據(jù)段信息，再按表2所示各信息的字節(jié)位置和字節(jié)長(zhǎng)度，定義其內(nèi)部的子變量數(shù)據(jù)類型，其核心代碼如下：

// 單個(gè)時(shí)間片數(shù)據(jù)段結(jié)構(gòu)定義

structSegmentData

{

charHx[3]; //Hx

charHy[3]; //Hy

charHz[3]; //Hz

charE1[4]; //E1

charE2[4]; //E2

charE3[4]; //E3

charE4[4]; //E4

charTF[2]; //TF

charTE[2]; //TE

charNo; //

}

// … …由于篇幅有限省略部分非關(guān)鍵代碼

// 定義變量

intHx=0,Hy=0,Hz=0,E1=0,E2=0,E3=0,E4=0;shortintTf=0,Te=0;

charNo;

TagDatatag;

SegmentDatasegment;

// 讀取儀器原始數(shù)據(jù)文件

fseek(infile_origin,0L,SEEK_SET);

cur_rc= 0;

intNo= 0;

// 如果當(dāng)前文件指針讀取位置不是結(jié)尾

while(cur_rc

{

if(i% 17 == 0)

{

fread(tag,sizeof(unsignedchar),32L,infile_origin);

cur_rc=ftell(infile_origin);

i++;

continue;

}

fread(segment,sizeof(unsignedchar),30L,infile_origin);

cur_rc=ftell(infile_origin);

i++;

j++;

// 將二進(jìn)制數(shù)值轉(zhuǎn)換為整型變量，詳見下文

convertSegmentData(segment,&Hx,&Hy,&Hz,&E1,&E2,&E3,&E4,&Tf,&Te,&No);

pre_rc=cur_rc;

}

2.3 二進(jìn)制數(shù)值到整型變量的轉(zhuǎn)換

LEMI-417原始數(shù)據(jù)是在儀器采集時(shí)存儲(chǔ)的二進(jìn)制數(shù)據(jù)，將它賦值到計(jì)算程序語言的整型變量時(shí)必須進(jìn)行二進(jìn)制位操作轉(zhuǎn)換。以磁通門傳感器的x方向磁場(chǎng)數(shù)據(jù)(Hx)為例，由于整型變量(int)無論是在32位還是64位操作系統(tǒng)其字節(jié)長(zhǎng)度均為32位，幸運(yùn)地避開了不同操作系統(tǒng)程序代碼的特殊處理問題，位轉(zhuǎn)換的實(shí)質(zhì)就是將原始數(shù)據(jù)中各字節(jié)位置的數(shù)值完整地賦值到整型變量的各字節(jié)位置上，其核心代碼如下：

inttemp1=0,temp2=0,temp3=0;

temp1 =segment.Hx[2];

temp1=temp1 ? 16;

temp1 =temp1 & 0x00FF0000;

temp2 =segment.Hx[1];

temp2 =temp2 ? 8;

temp2 =temp2 & 0x0000FF00;

temp3 =segment.Hx[0];

temp3 =temp3 & 0x000000FF;

Hx=temp1|temp2|temp3;

if(Hx>=0)

{

Hx=Hx& 0x00FFFFFF;

}

else

{

Hx=Hx| 0xFF000000;

}

3 應(yīng)用開發(fā)實(shí)例

在對(duì)原始數(shù)據(jù)格式成功剖析的基礎(chǔ)上，筆者針對(duì)該LEMI-417儀器在實(shí)際工作中存在的一些問題，補(bǔ)充性地開發(fā)出一套基于MFC的、可以在任何Windows系統(tǒng)運(yùn)行的界面程序，輔助一線野外工作人員對(duì)個(gè)別偶然會(huì)出現(xiàn)的數(shù)據(jù)問題進(jìn)行人為糾正。

實(shí)際應(yīng)用中，LEMI-417型儀器在無GPS信號(hào)時(shí)會(huì)出現(xiàn)記錄紊亂的現(xiàn)象[4]。在峽谷地貌或林木茂密區(qū)開展工作時(shí)可能會(huì)出現(xiàn)個(gè)別時(shí)段GPS衛(wèi)星信號(hào)偶然解鎖的現(xiàn)象，使得儀器無法同步當(dāng)前的GPS時(shí)鐘和測(cè)點(diǎn)坐標(biāo)，在對(duì)應(yīng)的文件頭信息中記錄下多組差別迥異的信息。例如，野外工作中儀器部署在深切山谷地貌區(qū)時(shí)，在部分時(shí)段會(huì)出現(xiàn)儀器GPS無法鎖定到≥3顆衛(wèi)星的情況，此時(shí)儀器記錄的測(cè)點(diǎn)經(jīng)緯度和高程均變?yōu)?，個(gè)別情況甚至?xí)霈F(xiàn)1顆衛(wèi)星都無法鎖定的情況，此時(shí)便會(huì)出現(xiàn)衛(wèi)星失鎖問題，無法同步到正確的GPS時(shí)鐘，使采集數(shù)據(jù)時(shí)序出現(xiàn)混亂，從而導(dǎo)致在后期數(shù)據(jù)導(dǎo)出時(shí)無法順利形成連續(xù)的時(shí)間序列數(shù)據(jù)，需要進(jìn)行2次返工，造成不必要的成本支出。

針對(duì)該常見問題，開發(fā)出了Lemi-417數(shù)據(jù)修改器(圖2)，可以將當(dāng)前測(cè)點(diǎn)采集的連續(xù)多天的數(shù)據(jù)加載、合并到一起，遍歷每個(gè)文件頭中的GPS時(shí)鐘、測(cè)點(diǎn)經(jīng)度、緯度、高程信息，判斷每段數(shù)據(jù)對(duì)應(yīng)的GPS時(shí)鐘是否是以設(shè)置的采樣率進(jìn)行連續(xù)觀測(cè)，會(huì)缺失的、不連續(xù)的時(shí)間序列進(jìn)行三次樣條插值修復(fù)，此外判斷每段數(shù)據(jù)對(duì)應(yīng)的經(jīng)度、緯度和高程信息是否出現(xiàn)較大的數(shù)據(jù)偏差，通過取眾數(shù)方法，識(shí)別出正常狀態(tài)下和故障狀態(tài)下的記錄，對(duì)故障狀態(tài)下的記錄進(jìn)行自動(dòng)修復(fù)。最后通過最小二乘法對(duì)正常狀態(tài)下的GPS經(jīng)度、緯度和高程進(jìn)行線性回歸擬合得到高精度的測(cè)點(diǎn)定位信息。由此，通過計(jì)算機(jī)程序解決了以往由于GPS偶然失鎖而需要重復(fù)返工的問題。

圖2 窗口程序界面截圖Fig.2 Screenshot of the window program interface

4 結(jié)論

隨著我國(guó)廣大學(xué)者對(duì)大地電磁時(shí)頻轉(zhuǎn)換算法研究的不斷深入[7-9]，國(guó)產(chǎn)化數(shù)據(jù)處理軟件[10]、儀器設(shè)備[11]的研制工作也已取得長(zhǎng)足進(jìn)步。本文深入剖析了LEMI-417I型地球深部探測(cè)主流觀測(cè)系統(tǒng)的原始數(shù)據(jù)格式，較詳細(xì)地給出了對(duì)應(yīng)的核心程序代碼，為廣大從事天然電磁場(chǎng)現(xiàn)代數(shù)字信號(hào)處理方面研究的學(xué)者提供一定的前期基礎(chǔ)支撐，對(duì)今后同類型國(guó)產(chǎn)儀器、軟件的研發(fā)工作具有一定的借鑒意義。此外，也針對(duì)LEMI-417I型儀器存在的GPS解鎖時(shí)的記錄紊亂問題，提出了一個(gè)經(jīng)濟(jì)有效的計(jì)算機(jī)輔助糾正技術(shù)方案。