陳挺，嚴(yán)迪，楊劍，廖國忠，武斌，冉中禹

(1.四川省地質(zhì)礦產(chǎn)勘查開發(fā)局物探隊(duì)，四川 成都 610072； 2.中國地質(zhì)調(diào)查局 成都地質(zhì)調(diào)查中心，四川 成都 610082)

0 引言

CG-5相對重力儀是加拿大Scintrex公司制造生產(chǎn)的的新型數(shù)字重力儀[1]，分辨率為0.001×10-5m/s2，重復(fù)率優(yōu)于0.005×10-5m/s2。廣泛應(yīng)用于重力研究、地震預(yù)報(bào)等領(lǐng)域[2]。當(dāng)?shù)卣鸹顒?dòng)突然爆發(fā)時(shí)，地下巖石的空間位置以及密度分布形態(tài)發(fā)生改變，從而引起觀測重力值的變化。吳云成等使用拉科斯特重力儀在遼寧西部地區(qū)進(jìn)行重力復(fù)測，分析了朝陽4.9級(jí)地震前后重力異常變化情況[3]。王武星分析了昆侖山口西8.1級(jí)地震和印尼8.9級(jí)地震前后部分重力觀測的變化，并對這些采樣間隔在1 min以下的重力數(shù)據(jù)進(jìn)行了高通濾波和頻譜分析[4]。車兆宏分析了張北6.2級(jí)地震前后華北地區(qū)水準(zhǔn)、GPS、重力、定點(diǎn)形變、流動(dòng)形變、地應(yīng)變觀測資料及前兆的演化過程[5]。祝意青等通過流動(dòng)重力測量，觀測到較可靠的重力隨時(shí)間變化數(shù)據(jù)[6-10]。梁偉鋒等基于河西地區(qū)流動(dòng)重力觀測資料獲得的重力場時(shí)空變化結(jié)果，研究發(fā)現(xiàn)重力場在活動(dòng)斷裂附近變化較大。這些觀測記錄表明重力變化能較好地反映地下深部物質(zhì)運(yùn)移與地殼密度變化等構(gòu)造活動(dòng)信息[11]。

長寧6.0級(jí)地震發(fā)生前，筆者正在云南省鹽津縣普洱鎮(zhèn)進(jìn)行重力靜態(tài)觀測，成功記錄到長寧6.0級(jí)地震發(fā)生前后的重力靜態(tài)觀測數(shù)據(jù)、重力儀實(shí)時(shí)觀測標(biāo)準(zhǔn)差(S.D.)以及舍棄數(shù)量(REJ)數(shù)據(jù)。截止2019年7月23日，在鹽津縣普洱鎮(zhèn)和珙縣上羅鎮(zhèn)重力儀成功記錄到地震活動(dòng)21次，說明長寧地震具有主震突出、余震豐富的特點(diǎn)。筆者在對長寧 6.0級(jí)地震前后重力儀觀測標(biāo)準(zhǔn)差數(shù)據(jù)和舍棄數(shù)量值分析的基礎(chǔ)上，認(rèn)為標(biāo)準(zhǔn)差值包括了重力儀自身噪聲和環(huán)境震動(dòng)等重力高頻干擾信號(hào)，可以用來衡量信號(hào)噪聲的大小，同時(shí)舍棄數(shù)量(REJ)的突變也指示了高頻強(qiáng)震動(dòng)干擾的存在。通過對比地震發(fā)生時(shí)重力儀在不同時(shí)刻記錄的標(biāo)準(zhǔn)差值，認(rèn)為標(biāo)準(zhǔn)差結(jié)合舍棄數(shù)量值可用于地震活動(dòng)不確定性測量，對長寧地區(qū)發(fā)震識(shí)別與記錄等具有科學(xué)與實(shí)際意義。

1 長寧震區(qū)地質(zhì)背景

根據(jù)1∶20萬區(qū)域地質(zhì)資料，長寧地區(qū)主要出露中、古生界地層(見圖1)，將區(qū)內(nèi)地層系統(tǒng)列于表1。巖漿巖主要出露晚二疊系玄武巖。玄武巖主要由致密狀、斑狀、杏仁狀玄武巖、凝灰角礫巖、凝灰?guī)r組成，礦物成分主要為拉長石、單斜輝石，主要分布于測區(qū)西部巡場鎮(zhèn)—壩底村一帶。長寧6.0級(jí)地震震中主要出露寒武系地層向奧陶系地層過渡的白云巖和灰?guī)r，珙縣5.1級(jí)地震震中主要出露志留系羅惹坪組灰?guī)r。

長寧地區(qū)主要地質(zhì)構(gòu)造為長寧背斜，該背斜東起敘永縣，西至珙縣天星橋一帶，呈不對稱的復(fù)式大背斜，背斜軸寬闊平緩，地表上在雙河鄉(xiāng)一帶。長寧背斜是一個(gè)多斷裂構(gòu)造的背斜，斷裂多發(fā)育在由古生界地層組成的背斜核部或兩翼，以小斷裂為主，呈高角度壓性逆沖斷層，震中附近斷層傾角最高達(dá)85°(圖1)。斷裂與褶皺緊密相伴，主要斷裂往往伴隨褶皺而發(fā)育，構(gòu)造的延伸方向與褶皺線形跡一致[12-13]。

圖1 長寧地區(qū)區(qū)域地質(zhì)圖(地質(zhì)底圖來源于91衛(wèi)圖助手)Fig.1 Regional geological map of Changning area(the geological base map is derived from the assistant of 91 satellite map)

表1 長寧地區(qū)地層系統(tǒng)Table 1 Stratigraphic system table of Changning area

2 觀測S.D.值數(shù)據(jù)可靠性分析

標(biāo)準(zhǔn)差在重力儀讀數(shù)記錄里用“S.D.”表示，重力儀進(jìn)行靜態(tài)觀測時(shí)，S.D.值每秒鐘更新一次。S.D.值的計(jì)算是根據(jù)已采集到的觀測重力值求取標(biāo)準(zhǔn)差，假設(shè)重力儀單點(diǎn)觀測持續(xù)時(shí)間設(shè)置為60 s，則重力儀記錄的S.D.值為60 s內(nèi)所有瞬時(shí)重力值讀數(shù)記錄的標(biāo)準(zhǔn)差。

為了保證重力儀讀數(shù)的準(zhǔn)確性和可重復(fù)性，當(dāng)打開“Auto Rejection”選項(xiàng)時(shí)，大于4～6倍標(biāo)準(zhǔn)差的高頻噪聲將被自動(dòng)拒絕。因此當(dāng)激活震動(dòng)濾波功能時(shí)，得到的重力觀測值是濾波處理后的結(jié)果。在外業(yè)觀測時(shí)，S.D.指標(biāo)主要用來衡量單點(diǎn)觀測質(zhì)量，當(dāng)S.D.較小時(shí)，表示環(huán)境震動(dòng)噪聲處于較低水平，設(shè)置觀測時(shí)間內(nèi)的數(shù)值較接近平均值，觀測質(zhì)量合格；當(dāng)S.D.較大時(shí)，表示環(huán)境震動(dòng)噪聲處于較高水平，觀測時(shí)間內(nèi)出現(xiàn)很多高頻重力信號(hào)，數(shù)據(jù)質(zhì)量不合格，需要重新觀測。正是由于S.D.值包括了重力儀自身噪聲和環(huán)境震動(dòng)等重力高頻干擾信號(hào)，它可以用來衡量信號(hào)噪聲的大小。通過對比地震發(fā)生時(shí)重力儀在不同時(shí)刻記錄的S.D.值，認(rèn)為S.D.值可用于地震活動(dòng)不確定性測量。

為了說明觀測S.D.值數(shù)據(jù)可靠性，在珙縣使用兩臺(tái)CG-5重力儀在同一地點(diǎn)進(jìn)行兩次觀測，重力儀編號(hào)分別是40242#、41036#。兩臺(tái)CG-5重力儀在使用前均進(jìn)行了格值標(biāo)定及各項(xiàng)性能試驗(yàn)，重力儀的一致性試驗(yàn)精度為0.022×10-5m/s2，試驗(yàn)精度結(jié)果滿足規(guī)范要求。

2019年7月22日16時(shí)26分在四川宜賓市長寧縣發(fā)生4.1級(jí)地震，震源深度9 km，震中位于北緯28.35°，東經(jīng)104.91°。兩臺(tái)CG-5重力儀布置在距離震中25.5 km的珙縣進(jìn)行靜態(tài)數(shù)據(jù)采集，表2是發(fā)震時(shí)兩臺(tái)重力儀記錄的不同時(shí)刻靜態(tài)觀測重力值數(shù)據(jù)和實(shí)時(shí)觀測重力值數(shù)據(jù)的標(biāo)準(zhǔn)差值(S.D.)，以及高頻噪聲舍棄數(shù)量(REJ)統(tǒng)計(jì)。

重力儀觀測持續(xù)時(shí)間DUR值均為60 s，為了對觀測S.D.值數(shù)據(jù)可靠性進(jìn)行分析，根據(jù)Excel軟件中的CORREL函數(shù)計(jì)算兩臺(tái)重力儀觀測S.D.值數(shù)據(jù)的相關(guān)系數(shù)為0.999 949，說明兩臺(tái)重力儀的觀測S.D.值數(shù)據(jù)具有較大的相關(guān)性。另外兩臺(tái)重力儀觀測開始時(shí)間僅相差1 s，從觀測S.D.值的數(shù)值大小變化來分析，40242#重力儀相對41036#重力儀在16∶26∶53、16∶28∶13、16∶29∶33三個(gè)不同時(shí)刻的相對變化率分別為4.34%、1.51%、10.34%。說明利用S.D.值來分析地震活動(dòng)性是可靠的。

2019年7月23日22時(shí)25分在四川宜賓市長寧縣發(fā)生3.0級(jí)地震，震源深度6 km，震中位于北緯28.38°，東經(jīng)104.94°。兩臺(tái)CG-5重力儀布置在距離震中30.2公里的珙縣進(jìn)行靜態(tài)數(shù)據(jù)采集，表3是發(fā)震時(shí)兩臺(tái)重力儀記錄的不同時(shí)刻靜態(tài)觀測重力值數(shù)據(jù)和實(shí)時(shí)觀測重力值數(shù)據(jù)的標(biāo)準(zhǔn)差值以及高頻噪聲舍棄數(shù)量值統(tǒng)計(jì)表。重力儀觀測持續(xù)時(shí)間DUR值均為60 s，對兩臺(tái)重力儀的觀測S.D.值數(shù)據(jù)進(jìn)行了相關(guān)分析，相關(guān)系數(shù)為0.997 503，說明兩臺(tái)重力儀的觀測S.D.值數(shù)據(jù)具有較大的相關(guān)性。表3觀測S.D.值數(shù)據(jù)的相關(guān)系數(shù)小于表2，這主要是由于表3兩臺(tái)重力儀的觀測開始時(shí)間相差20 s，即在各自的60 s觀測時(shí)段內(nèi)僅同步觀測40 s所致。

表2 長寧4.1級(jí)地震重力儀觀測數(shù)據(jù)統(tǒng)計(jì)Table 2 Statistical table of observation data of the gravimeter of Ms 4.1 earthquake in Changning

表2的S.D.值變化幅值較表3大，最大達(dá) 1.404×10-5m/s2，兩臺(tái)重力儀發(fā)震時(shí)刻經(jīng)固體潮和零漂改正后的相對重力值較發(fā)震前的變化量分別為-0.007×10-5m/s2、-0.009×10-5m/s2，說明觀測時(shí)誤差較大的高頻重力數(shù)據(jù)已經(jīng)被重力儀自帶的濾波器過濾掉。表3的S.D.值變化幅值較小，發(fā)震時(shí)刻經(jīng)固體潮和零漂改正后的相對重力值較發(fā)震前的變化量分別為0.020×10-5m/s2、0.023×10-5m/s2，說明部分由小震動(dòng)引起的小于4倍標(biāo)準(zhǔn)差的重力觀測數(shù)據(jù)未被過濾掉。據(jù)表2和表3的舍棄數(shù)量(REJ)值統(tǒng)計(jì)結(jié)果知，發(fā)震時(shí)刻大致與REJ極大值相對應(yīng)，說明發(fā)震時(shí)的高頻震動(dòng)干擾最強(qiáng)。41036#儀器在同一時(shí)刻的REJ值均大于或等于40242#儀器，且表2和表3中41036#儀器REJ值在地震發(fā)生前有明顯變化，而40242#儀器REJ值在地震發(fā)生前無明顯變化，說明REJ值與重力儀的彈簧彈性系數(shù)相關(guān)，從REJ值的統(tǒng)計(jì)結(jié)果知，41036#儀器對高頻干擾噪聲更敏感，彈簧靈敏度明顯高于40242#儀器，因此用于地震監(jiān)測的重力儀應(yīng)選擇彈簧靈敏度較高的儀器。由于REJ值與重力儀自身的彈簧彈性系數(shù)相關(guān)，如需使用REJ值來預(yù)測地震，需針對單臺(tái)重力儀的大量試驗(yàn)數(shù)據(jù)來選擇合適的REJ值判斷指標(biāo)。

表3 長寧3.0級(jí)地震重力儀觀測數(shù)據(jù)統(tǒng)計(jì)Table 3 Statistical table of observation data of the gravimeter of Ms 3.0 earthquake in Changning

3 觀測S.D.值與地震震級(jí)相關(guān)性研究

截止2019年7月23日，41036#重力儀共記錄到21次地震，見圖2和表4。

表4是發(fā)震時(shí)重力儀記錄靜態(tài)觀測數(shù)據(jù)統(tǒng)計(jì)表。在鹽津縣觀測的15次記錄中，當(dāng)時(shí)震中距介于60.1～73.2 km時(shí)，為了研究地震震級(jí)與重力儀記錄S.D值的關(guān)系，根據(jù)Excel軟件中的CORREL函數(shù)計(jì)算了地震震級(jí)與對應(yīng)S.D值的相關(guān)系數(shù)，計(jì)算得出相關(guān)系數(shù)為0.90，由于相關(guān)系數(shù)是研究變量之間線性相關(guān)程度的量，相關(guān)系數(shù)值越大，說明地震震級(jí)與S.D值之間線性相關(guān)程度越高。

在珙縣觀測的6次記錄中，當(dāng)震中距介于25.5 ～30.8 km，相關(guān)系數(shù)為0.94；當(dāng)震中距介于6.2～7.8 km，相關(guān)系數(shù)為1。需要指出的是，相關(guān)系數(shù)有一個(gè)明顯的缺點(diǎn)，即當(dāng)震中距介于6.2～7.8 km，僅有7月15日和7月19日兩組統(tǒng)計(jì)樣本數(shù)據(jù)，相關(guān)系數(shù)的絕對值為1，因此在觀測記錄較少時(shí)，僅憑相關(guān)系數(shù)較大就判定震級(jí)與重力儀記錄S.D值之間有密切的線性關(guān)系是不妥當(dāng)?shù)摹Ｍㄟ^以上相關(guān)分析認(rèn)為，當(dāng)震中距介于60.1～73.2 km、25.5～30.8 km時(shí)，地震震級(jí)與重力儀記錄S.D值呈明顯的正相關(guān)關(guān)系，即在震中距相近的情況下，S.D值越大，地震震級(jí)越高。當(dāng)震中距介于6.2～7.8 km時(shí)，由于統(tǒng)計(jì)樣本數(shù)據(jù)不足，未能準(zhǔn)確判定地震震級(jí)與重力儀記錄S.D值的相關(guān)關(guān)系。

圖2 長寧地震震中分布位置與地形疊合Fig.2 The overlap map of topography and epicenter distribution of Changning earthquake

表4 發(fā)震時(shí)刻重力儀記錄觀測數(shù)據(jù)統(tǒng)計(jì)Table 4 Statistical table of observation data recorded by gravimeter at seismic moment

表4的舍棄數(shù)量(REJ)值與震級(jí)和重力儀記錄S.D值無明顯的相關(guān)關(guān)系，這可能是由于當(dāng)觀測持續(xù)時(shí)間選擇60 s時(shí)，REJ值為60 s內(nèi)的統(tǒng)計(jì)值，無法對瞬時(shí)高頻干擾進(jìn)行更精細(xì)的反映。REJ值的突變顯示了高頻干擾的存在，在今后工作中，可以將觀測持續(xù)時(shí)間選擇為30 s，來進(jìn)一步統(tǒng)計(jì)分析REJ值與震級(jí)和重力儀記錄S.D值的相關(guān)關(guān)系。

4 利用S.D.值進(jìn)行地震預(yù)測可行性分析

2019年6月17日22時(shí)55分在四川宜賓市長寧縣發(fā)生6.0級(jí)地震， 06月17日23時(shí)36分在四川宜賓市珙縣發(fā)生5.1級(jí)地震。從表5知，22時(shí)55分重力儀記錄S.D.值高達(dá)23.066×10-5m/s2，REJ值為10%。需注意的是，發(fā)震前十分鐘即22時(shí)45分，S.D.值僅為0.021×10-5m/s2，但REJ值卻高達(dá)19%，REJ值的突變能否作為地震發(fā)生前兆的判斷依據(jù)還需更多的試驗(yàn)記錄來支撐。23時(shí)36分珙縣5.1級(jí)地震繼長寧6.0級(jí)地震后爆發(fā)，S.D.值突變?yōu)?1.708×10-5m/s2，直到06月18日00時(shí)06分，S.D.值降為0.013×10-5m/s2，表明該次地震活動(dòng)已經(jīng)平息。另外從觀測重力值來分析，長寧縣6.0級(jí)地震和珙縣5.1級(jí)地震突然爆發(fā)引起的強(qiáng)烈震動(dòng)，導(dǎo)致重力儀自帶的地震濾波功能失效，觀測重力值與地震爆發(fā)前觀測值的差值分別為4.466×10-5m/s2、-0.055×10-5m/s2。23時(shí)46分S.D.值為0.104×10-5m/s2，近似為平靜場中觀測值的5倍，但觀測重力值與地震爆發(fā)前的觀測值相差僅為-0.001×10-5m/s2，可見重力儀自帶的地震濾波功能已經(jīng)濾掉了重力高頻干擾信號(hào)。剔除發(fā)震時(shí)刻的畸變重力數(shù)據(jù)，在云南省鹽津縣觀測到的長寧縣6.0級(jí)地震、珙縣5.1級(jí)地震并未引起觀測重力觀測值的明顯變化。直至06月18日00時(shí)06分S.D.值變?yōu)?.013×10-5m/s2，說明兩次地震震動(dòng)疊加在一起的持續(xù)活動(dòng)時(shí)間約為71分鐘。

在記錄到的21次地震中，有2次在發(fā)震前重力儀記錄S.D.值有明顯突變。第1次是2019年6月20日01時(shí)44分在四川宜賓市珙縣發(fā)生的3.4級(jí)地震(見圖3)，在地震發(fā)生前的01時(shí)38分，S.D.值變?yōu)?.043×10-5m/s2，該時(shí)刻經(jīng)固體潮和零點(diǎn)位移改正后的靜態(tài)觀測重力值為4757.394×10-5m/s2。從表5知，S.D.值與上一記錄時(shí)刻相差0.024×10-5m/s2，重力值與上一記錄時(shí)刻相差-0.002×10-5m/s2，重力值無明顯變化，REJ值也無變化。由于靜態(tài)觀測時(shí)設(shè)置的讀數(shù)時(shí)間間隔是10 min，未能記錄到01時(shí)27分到01時(shí)37分之間的S.D.值，因此S.D.值可能在發(fā)震6 min之前就已發(fā)生突變。直至02時(shí)18分S.D.值變?yōu)?.015×10-5m/s2，說明此次地震活動(dòng)從發(fā)震到停止用時(shí)約34 min。

第2次是2019年7月15日16時(shí)48分在四川宜賓市珙縣發(fā)生的2.9級(jí)地震。為了更詳細(xì)地研究重力場和S.D.值隨時(shí)間的變化規(guī)律，重力儀的讀數(shù)時(shí)間間隔在鹽津縣觀測時(shí)設(shè)置為10 min，在珙縣觀測時(shí)間設(shè)置為1 min。從圖4知，在發(fā)震時(shí)刻S.D.值變?yōu)?.935×10-5m/s2，由于發(fā)震時(shí)刻S.D.值幅值變化太大，在圖4中很難觀測到發(fā)震前后S.D.值的微小變化。

圖5中剔除了16時(shí)48分的S.D.值后，可以明顯觀測到，在地震發(fā)生前的16時(shí)30分和16時(shí)32分，S.D.值發(fā)生了明顯改變，其值分別為0.048×10-5m/s2、0.040×10-5m/s2，經(jīng)固體潮和零點(diǎn)位移改正后的靜態(tài)觀測重力值分別為4 830.031×10-5m/s2、4 830.043×10-5m/s2，重力值與上一記錄時(shí)刻相差0.012×10-5m/s2，說明在此時(shí)刻由震動(dòng)產(chǎn)生的部分高頻重力異常未被重力儀自帶的濾波器濾除掉。在16時(shí)30分，S.D.值近似為平靜場中觀測值的2.5倍，即在發(fā)震前18 min重力儀已記錄到S.D.值的突變。值得注意的是16時(shí)40分，即在發(fā)震前8 min，S.D.值僅為0.022×10-5m/s2，重力值與上一記錄時(shí)刻相差-0.033×10-5m/s2，REJ值高達(dá)42%，根據(jù)重力值和REJ值的變化可對地震進(jìn)行提前預(yù)報(bào)。在直至16時(shí)54分S.D.值變?yōu)?.014×10-5m/s2，說明此次地震活動(dòng)從發(fā)震到停止用時(shí)約6 min。

圖3 珙縣3.4級(jí)地震標(biāo)準(zhǔn)差隨時(shí)間變化Fig.3 The change of S.D. value over time in Gongxian Ms 3.4 earthquake

表5 地震前后重力儀記錄靜態(tài)觀測數(shù)據(jù)統(tǒng)計(jì)Table 5 Statistical table of static observation data recorded by gravimeter before and after earthquake

圖4 珙縣2.9級(jí)地震標(biāo)準(zhǔn)差隨時(shí)間變化Fig.4 The change of S.D. value over time in Gongxian Ms 2.9 earthquake

圖5 珙縣2.9級(jí)地震標(biāo)準(zhǔn)差隨時(shí)間變化(剔除發(fā)震時(shí)刻S.D.突變點(diǎn))Fig.5 The change of S.D. value over time in Gongxian Ms 2.9 earthquake(the mutation point of S.D. at the time of shock was removed)

通過以上實(shí)際記錄S.D.值數(shù)據(jù)分析認(rèn)為，在某些特定條件下，根據(jù)S.D.值的突變并結(jié)合重力值和REJ值的變化來提前預(yù)報(bào)地震是可行的。

5 結(jié)論與建議

通過對地震發(fā)生時(shí)重力儀在不同時(shí)刻記錄S.D.值和REJ值的分析，認(rèn)為S.D.值結(jié)合REJ值可用于地震活動(dòng)不確定性測量。首先，對兩臺(tái)重力儀兩次實(shí)測S.D.值數(shù)據(jù)進(jìn)行了相關(guān)分析，相關(guān)系數(shù)為0.99，說明兩臺(tái)重力儀的觀測S.D.值數(shù)據(jù)具有較大的相關(guān)性，利用S.D.值來分析地震活動(dòng)性是可靠的。然后，通過對地震震級(jí)與發(fā)震時(shí)刻重力儀記錄S.D值的分析，認(rèn)為當(dāng)震中距介于60.1～73.2 km、25.5～30.8 km時(shí)，地震震級(jí)與重力儀記錄S.D值呈正相關(guān)關(guān)系。REJ值與重力儀自身的彈簧彈性系數(shù)相關(guān)，如需使用REJ值來預(yù)測地震，需針對單臺(tái)重力儀的大量試驗(yàn)數(shù)據(jù)來選擇合適的REJ值判斷指標(biāo)。最后，在記錄到的21次地震中，有2次在發(fā)震前重力儀記錄S.D值有明顯突變，S.D.值近似為平靜場中觀測值的2倍；有1次在發(fā)震前8 min，觀測到REJ值的突變。因此在某些特定條件下，根據(jù)S.D.值的突變并結(jié)合重力值和REJ值的變化來預(yù)報(bào)地震和研究單次地震從發(fā)震到結(jié)束的持續(xù)時(shí)間是可行的。

本次長寧6.0級(jí)地震，在地震波到達(dá)之前未能觀測到有效的重力異常變化信號(hào)，地震前后未觀測到明顯的重力值變化。這可能是由于本次地震未引起地下地質(zhì)體空間分布形態(tài)的較大改變且測站距離震中距離較遠(yuǎn)所致。在鹽津縣觀測時(shí)，重力儀觀測時(shí)間間隔設(shè)置為10 min，數(shù)據(jù)采集密度不夠，未能在地震發(fā)生前記錄到S.D.值的突變，但在發(fā)震前10 min記錄到REJ值的突變。在以后的工作中，重力儀觀測時(shí)間間隔設(shè)置為30～60 s，可以觀測到有效的重力值和S.D.值，從而進(jìn)一步研究地震活動(dòng)不確定性規(guī)律。另外在震區(qū)進(jìn)行小震中距重力靜態(tài)觀測時(shí)，地震震級(jí)與重力儀記錄S.D.值是否也呈正相關(guān)關(guān)系，在進(jìn)一步分析中是否要加入震源深度作為約束？這需要通過更多的采樣數(shù)據(jù)來證實(shí)。