楊錦玲

1.閩江學(xué)院測繪工程系,福建 福州 350108; 2. 中國地震局地球物理研究所,北京 100081

陸地時變重力觀測距離場源近,觀測精度高,目前廣泛應(yīng)用于監(jiān)測地殼內(nèi)部物質(zhì)的變化。陸地重力觀測通過在陸地表面固定測點布設(shè)重力儀進(jìn)行重力測量工作以獲取地球重力場數(shù)據(jù)。但在測量過程中會受到觀測儀器誤差、水文和環(huán)境干擾等多種因素的綜合影響。對相對重力儀的儀器參數(shù)誤差、測網(wǎng)的場源監(jiān)測能力等不確定性因素的量化和研究有助于獲得更可靠的地下深部場源信號。論文基于貝葉斯重力平差方法和時變重力等效源時空平滑方法,從儀器參數(shù)的不確定性估計、測網(wǎng)的場源監(jiān)測能力評估以及重力場源因素和模型解釋方面開展了陸地時變重力觀測研究。論文主要工作如下。

(1) 利用貝葉斯重力平差方法對華南地區(qū)2015—2018年的陸地多期重力觀測數(shù)據(jù)質(zhì)量進(jìn)行評估,量化分析了相對重力儀的非線性漂移和格值系數(shù)的影響。結(jié)果表明4臺相對重力儀存在非線性漂移,且儀器間差異性明顯。段差的殘差結(jié)果和統(tǒng)計參數(shù)顯示貝葉斯平差結(jié)果優(yōu)于經(jīng)典方法。初始與優(yōu)化的格值系數(shù)偏差在10-5～10-4量級。華南測網(wǎng)觀測資料由于相對重力儀非線性漂移和格值系數(shù)不準(zhǔn)確共可引起約20×10-8m/s2的不確定性,貝葉斯重力平差方法可以有效抑制儀器誤差。開展0.25°、0.5°和1° 3種不同尺度的檢測板模型試驗結(jié)果表明,現(xiàn)今的華南重力測網(wǎng)可以有效識別測網(wǎng)內(nèi)空間尺度超過0.5°×0.5°的場源體,對應(yīng)的場源空間分辨率約為55 km。

(2) 對2017年九寨溝Ms 7.0地震前的時變重力場開展研究?；?014—2017年的7期重復(fù)高精度混合重力觀測數(shù)據(jù),采用貝葉斯平差方法獲得了2017年九寨溝Ms 7.0地震前震源區(qū)及周邊空間分辨率優(yōu)于75 km尺度的微伽級時變重力信號。絕對重力交叉驗證結(jié)果表明,貝葉斯重力平差后的絕對點偏差在±20×10-8m/s2范圍,貝葉斯重力平差方法可以有效減小儀器格值和漂移引起的不確定性。采用等效源時空平滑方法對青藏高原東緣重力測網(wǎng)的檢測板測試結(jié)果顯示青藏高原東緣測網(wǎng)的場源分辨能力優(yōu)于75 km。噪聲測試結(jié)果表明:在5%的噪聲水平下(對應(yīng)地表約12×10-8m/s2的高頻噪聲),等效源時空平滑方法可以有效地恢復(fù)場源參數(shù)?；谠摲椒▽?.75°分辨率的場源體等效到15 km、深度1 km厚度水平層內(nèi),對應(yīng)場源體0.59‰～0.74‰的視密度變化。在九寨溝Ms 7.0地震前,顯著的重力場變化區(qū)域從西部向東北變化,逐漸出現(xiàn)顯著的與構(gòu)造方向一致的重力變化梯度帶。對震中周邊等效場源結(jié)果的分析表明,在震中西南區(qū)域震前兩年可能存在地殼流體的遷移。視密度變化的梯度帶與活動地殼邊界區(qū)重合,顯示出引起重力場變化的場源主要受構(gòu)造控制。地殼垂向變形和近地表水變化對重力變化的影響較小。在2017年九寨溝Ms 7.0地震前,與構(gòu)造相關(guān)的區(qū)域重力場空間變化過程和變化特征主要反映了該區(qū)域地殼深部的流體和物質(zhì)的遷移過程。

(3) 對2016年門源Ms 6.4地震前的時變重力場開展研究。平差結(jié)果顯示2016年門源Ms 6.4地震發(fā)生在重力差異變化顯著的高梯度帶,且位于零值線拐彎部位上。采用等效源時空平滑方法對青藏高原東緣重力測網(wǎng)的檢測板測試和噪聲測試表明,青藏高原東北緣測網(wǎng)的場源分辨能力約為100 km,可觀測的重力異常量級約為±80×10-8m/s2。通過等效源時空平滑反演方法,1°分辨率的場源體可以等效到10 km、深度1 km厚度等效層內(nèi)對應(yīng)0.74‰～0.81‰的場源體視密度變化。區(qū)域GNSS速度場顯示青藏高原東北緣的水平運動整體呈北東方向和東向,量級在2.6～11.8 mm/a,從西南部向東北部速度呈逐步減小的趨勢。冷龍嶺斷裂帶具有強烈的活動特征,斷裂帶兩側(cè)的GNSS運動速度量級差異較大,在門源震中附近應(yīng)變積累較顯著。區(qū)域主應(yīng)變特征的優(yōu)勢分布為主壓應(yīng)變北東向,總體上表現(xiàn)為擠壓型變形,面應(yīng)變率主要為收縮和擠壓的特征。門源震中位于高剪應(yīng)變和高面應(yīng)變的峰值附近,其最大剪應(yīng)變年變化率約為1.4×10-8mm,面收縮年變化率為0.9×10-8mm。區(qū)域近土壤水和雪水等近地表水引起的重力年變化速率約為0.05×10-8m/s2。降雨引起的重力最大變化為6×10-8m/s2,最大地表垂直形變速率引起的重力變化為1.17×10-8m/s2。水文因素和地表垂直形變引起的時變重力變化大約為8×10-8m/s2。體波成像等結(jié)果顯示在門源地震下方30～40 km處存在著低速層,推測該區(qū)地殼內(nèi)可能存在著流體及流體的遷移。電磁測深的結(jié)果顯示冷龍嶺斷裂處于青藏高原北東向擠壓推擠作用下,斷裂南北兩側(cè)存在顯著的物性差異。2016年門源Ms 6.4地震前的重力場變化可能與地殼深部分布的低速層物質(zhì)向淺部遷移有關(guān)。