鄂爾多斯南緣地區(qū)重力變化場源特征

張永奇 韓美濤 鄭增記 曹建平

1）中國西安 710068 陜西省地震局2）中國北京 100029 中國地震局地質(zhì)研究所地震動力學(xué)國家重點實驗室

0 引言

地震孕育過程中產(chǎn)生的地殼變形和震源介質(zhì)變化在地表會引起重力變化，通過開展連續(xù)重力觀測或定期相對重力復(fù)測，有可能捕捉到地震的這些前兆信息。重力場具有明確的地球物理含義，屬于構(gòu)造運動的一種表現(xiàn)形式，一般分為“靜態(tài)”重力場和“動態(tài)”重力場?！办o態(tài)”重力場是指，基于地球重力模型或?qū)崪y重力數(shù)據(jù)，獲取布格重力及均衡重力異常信息，用來研究地質(zhì)構(gòu)造特征、地殼物性密度變化及地殼均衡狀態(tài)等地質(zhì)時間尺度的地學(xué)問題（陳石等，2011a，2011b，2014）；“動態(tài)”重力場是指，在年時間尺度內(nèi)，利用流動作業(yè)方式獲取的重力異常變化信息，主要用來認識強震孕育機理、捕捉地震前兆信息和判定地震危險區(qū)（祝意青等，2009，2017，2020；申重陽等2010，2011；陳石等，2015；Chen et al，2016a）。

當(dāng)前主要基于“場”模式開展地震監(jiān)測預(yù)測，而應(yīng)用地表觀測的重力場變化信息來開展異常源的研究，即“以場求源”問題的研究仍不夠深入。眾所周知，地表觀測的重力場變化是測量誤差、地表水影響、測點三維形變、地下介質(zhì)密度變化及物質(zhì)遷移等信息的綜合反映（陳運泰等，1980；陳石等，2011c）。準(zhǔn)確分離并識別這些因素對重力場變化信息的影響，是比較困難的事情，原因如下：①深部介質(zhì)變化的先驗信息及參數(shù)難以觀測或獲?。虎诂F(xiàn)有重力儀器觀測精度偏低，信噪比不高。此外，在重力場變化的相關(guān)研究中，經(jīng)常采用大地測量學(xué)方法進行數(shù)據(jù)計算和解釋，間接造成了上述困難無法克服。為此，嘗試將三維歐拉反褶積方法引入重力場變化分析。該方法是Reid 等（1990）結(jié)合Thompson（1982）在歐拉齊次方程的研究基礎(chǔ)上建立起來的，其最大優(yōu)勢是，在先驗信息缺乏情況下，仍能較好完成位場數(shù)據(jù)反演與解釋工作（范美寧，2006）。

文中對鄂爾多斯南緣地區(qū)EGM2008 重力場模型數(shù)據(jù)及2014—2017 年相對重力觀測數(shù)據(jù)進行精細化處理。在此基礎(chǔ)上，以布格重力異常、1 年尺度差分重力場和3 年尺度累積重力場為研究對象，分析上述重力場異常變化的場源特征。設(shè)計一個與實際重力觀測相似的理論模型模擬地表重力場變化，同時利用三維歐拉反褶積方法對場源參數(shù)進行反演，優(yōu)選與模型接近的構(gòu)造指數(shù)及滑動窗口大小等參數(shù)，對相對重力觀測網(wǎng)數(shù)據(jù)進行歐拉反演，并利用水平梯度濾波法對發(fā)散的反演結(jié)果進行優(yōu)化，結(jié)合地震目錄資料，對獲取的三維場源信息進行分析與解釋。

1 測區(qū)概況與數(shù)據(jù)處理

鄂爾多斯南緣地區(qū)位于我國青藏塊體與華北塊體交匯區(qū)域，地質(zhì)構(gòu)造形式多樣，動力學(xué)環(huán)境復(fù)雜。同時，該地區(qū)位于我國重要的汾渭地震帶和南北地震帶北段，歷史上曾發(fā)生多次6 級以上大震，屬于我國地震活動重點監(jiān)視區(qū)。該區(qū)發(fā)育多條活動斷裂，其中鄂爾多斯南緣斷裂、渭河斷裂及秦嶺北緣斷裂控制著該區(qū)主要構(gòu)造變形演化格局，區(qū)內(nèi)其他斷裂進一步把鄂爾多斯南緣地區(qū)分隔為不同樣式的地質(zhì)構(gòu)造單元，單條或多條斷裂邊界或交匯處成為地震孕育的主要場所。為了監(jiān)測鄂爾多斯南緣地區(qū)地震活動潛在危險性，尤其是渭河斷陷盆地內(nèi)主要斷裂的地球物理場異常變化，陜西省地震局于20 世紀(jì)80 年代建立了相對重力觀測網(wǎng)，每年觀測1 期。2008 年汶川8.0 級地震后，依托“震后恢復(fù)重建”等項目，對原有相對重力觀測網(wǎng)進行優(yōu)化改造，改造后的重力測網(wǎng)覆蓋陜西全境，加上“陸態(tài)網(wǎng)絡(luò)”一、二期建立的重力站點，重力站點總計達200 余座，點間距基本在20—30 km，站點分布見圖1。這些重力站點覆蓋渭河盆地及鄰區(qū)大部分?jǐn)嗔?，且分布相對合理，為監(jiān)測斷裂活動引起的重力變化奠定了堅實的基礎(chǔ)（曹建平等，2014；張永奇等，2014）。

圖1 鄂爾多斯南緣地區(qū)重力點及斷裂分布Fig.1 Distribution of gravity points and faults in the southern margin of the Ordos

基于重力資料開展鄂爾多斯南緣地區(qū)的研究，以往集中在數(shù)據(jù)精細化處理、重力場變化及重力網(wǎng)的映震能力分析等方面，針對深部場源特征的研究相對較少。因此，嘗試基于EGM2008 重力場模型數(shù)據(jù)及相對重力數(shù)據(jù)，利用歐拉反褶積方法，對該區(qū)進行場源特征反演。反演前將相對重力觀測數(shù)據(jù)進行平差計算，具體流程為：①采用中國地震局實用化攻關(guān)推廣軟件——LGADJ，對多期重力觀測資料進行獨立平差計算，解算重力儀一次項周期誤差等參數(shù)的區(qū)域適定解，計算各臺儀器觀測精度，合理確定其先驗方差；②以四川廣元、西安子午、甘肅天水等絕對測點為起算基準(zhǔn)，采用經(jīng)典平差方法進行平差計算，對其中因周圍環(huán)境改變而出現(xiàn)異常的個別測點進行改正或刪除；③采用克里金插值算法擬合離散重力場，獲得不同時間尺度的重力場。

2 方法原理

歐拉反褶積方法不依賴場源密度參數(shù)，僅與場源位置的幾何參數(shù)相關(guān)，使得反演問題得以簡化（Reid et al，1990）。通過求解歐拉方程，可以確定場源位置的三維空間位置，并可對不同類型地質(zhì)體給出參數(shù)估計，公式如下

構(gòu)造指數(shù)N反映了位場異常強度隨場源深度變化的衰減率，一般與場源幾何結(jié)構(gòu)有關(guān)。N的選取是場源位置確定的關(guān)鍵。根據(jù)不同構(gòu)造特征，針對重力場異常數(shù)據(jù)，N的范圍一般取0—2 較為合適，其中斷層、巖脈、臺階等構(gòu)造指數(shù)皆不同（范美寧，2006；高倩等，2015；劉芳等，2016；翟麗娜等，2020）。三維重力歐拉反褶積方法的反演流程如下：①對實測離散重力位場數(shù)據(jù)進行網(wǎng)格化；②計算重力位場異常在x、y、z三個方向的導(dǎo)數(shù)；③選擇適當(dāng)?shù)臉?gòu)造指數(shù)N及窗口大小。根據(jù)異常體性質(zhì)選擇構(gòu)造指數(shù)N，窗口大小根據(jù)異常體大小進行選擇，一般應(yīng)覆蓋異常范圍，不得低于3 倍測點間距。根據(jù)精度要求確定允許誤差；④解方程組，求出場源位置參數(shù)及背景值；⑤根據(jù)解的匯聚位置、深度調(diào)整參數(shù)。如果對反演結(jié)果滿意，結(jié)束計算，否則重新選擇參數(shù)重復(fù)步驟②、③，直至對結(jié)果滿意為止（高倩等，2015；劉芳等，2016；翟麗娜等，2020）。

3 模型反演

歐拉反褶積方法不需要已知場源先驗信息，即可對場源參數(shù)進行反演與解釋。反演計算時，應(yīng)選擇合適的構(gòu)造指數(shù)N，并求出異常在x、y、z三個方向的導(dǎo)數(shù)，即可組成線性方程組進行求解（范美寧，2006）。在實際應(yīng)用中，地質(zhì)情況復(fù)雜，場源形態(tài)也不確定，所以如何選擇構(gòu)造指數(shù)，需要對不同模型進行反復(fù)計算來確定（魯寶亮等，2009）。此外，滑動窗口大小的選取，一般需考慮測點間距、異常形態(tài)特征等，在一個窗口內(nèi)應(yīng)有足夠多的測點滿足方程組的求解，并要覆蓋一個場源異常形態(tài)。對于反演結(jié)果的評價，需要通過求得場源位置參數(shù)的水平誤差和垂直誤差大小，及反演場源位置解的三維空間聚集程度和與地質(zhì)認識是否符合等判斷，即場源位置越集中，其精度及收斂性越好，反之解的收斂性差（高倩等，2015；劉芳等，2016）。

基于對歐拉反褶積法的認識，根據(jù)鄂爾多斯地塊南緣地區(qū)現(xiàn)有重力觀測網(wǎng)形態(tài)及規(guī)模，設(shè)計一個相對符合實際觀測情況的理論模型。以3 個直面正六面體組合模型為研究模型，假設(shè)其分別位于上地殼、中地殼和下地殼。其中模型1（Model1）位于上地殼，單元尺寸為20 km × 20 km × 20 km，中心埋深為5.5 km；模型2（Model2）位于中地殼，單元尺寸為10 km×10 km×10 km，中心埋深為11.5 km；模型3（Model3）位于下地殼，單元尺寸為30 km×30 km×30 km，中心埋深為22.5 km。3 個模型密度變化統(tǒng)一取為3.0 g·cm-3，正演計算模型在地表重力測網(wǎng)的理論異常值，用于模擬觀測的重力變化信號和優(yōu)選反演參數(shù)，模型異常見圖2，矩形框表示場源模型在地表的水平投影。地表重力測網(wǎng)點間距設(shè)計為5 km，測網(wǎng)長度在X、Y方向均為200 km、深度方向為50 km，3 個不同深度的場源模型參數(shù)見表1，可見3個模型的尺度特征不同，埋深也不同，正演結(jié)果顯示理論異常范圍在(0—120)×10-8m·s-2。

圖2 不同深度場源組合模型及正演異常Fig.2 Combination models of different depth field sources and forward modeling anomalies

表1 理論場源模型參數(shù)Table 1 The theoretical field source model parameters

基于上述理論測試結(jié)果，實際測試時設(shè)計多個參數(shù)組合，其中構(gòu)造指數(shù)分別選擇0.5、1、1.5、2，滑動窗口分別選擇4—8 倍的測點距離。通過對多次試驗結(jié)果進行對比分析，認為選擇構(gòu)造指數(shù)N=1，滑動窗口大小在4 個測點距離時，反演結(jié)果的聚集程度更符合模型的幾何結(jié)構(gòu)，組合模型的反演深度誤差范圍限制在10%以內(nèi)。圖3 給出計算所得一組最佳反演結(jié)果，圖中彩色圓圈為反演的場源參數(shù)，共獲得467 組解，黑色矩形框為理論組合模型位置，可見聚集程度相對較好，場源特征明顯。其中模型2（Model2）和模型3（Model3）反演結(jié)果最好，模型1（Model1）反演結(jié)果稍差。

圖3 基于歐拉反褶積的反演結(jié)果空間分布Fig.3 Spatial distribution of inversion results based on Euler deconvolution

4 歐拉反演解優(yōu)化

歐拉反演的基本條件是模型滿足歐拉齊次方程的假設(shè)條件，此時反演結(jié)果應(yīng)該是準(zhǔn)確的。然而實際上地質(zhì)體的位置和構(gòu)成情況復(fù)雜，無法用任何一種模型準(zhǔn)確表達，這就會導(dǎo)致反演解出現(xiàn)發(fā)散現(xiàn)象（姚長利等，2004；管志寧等，2005；石磊等，2008），因此需要尋找一種減弱或消除反演解發(fā)散的方法，研究表明水平梯度濾波法是一種比較有效的方法。在利用歐拉反褶積方法開展反演研究時，當(dāng)窗口滑動到異常體邊緣時，梯度值變小甚至趨近于零，此時反演容易出現(xiàn)發(fā)散解。理論上，水平梯度決定了場源的水平位置，構(gòu)造指數(shù)決定了場源的深度，因此構(gòu)造指數(shù)越大，反演的深度越深，結(jié)果將更加發(fā)散（姚長利等，2004；鄭秋月等，2021），使用水平梯度濾波法即可克服上述問題。

設(shè)Hi=為研究區(qū)任意測點i處的重、磁異常水平梯度模，則該區(qū)域內(nèi)重磁異常的水平梯度模為n為研究區(qū)網(wǎng)格點總數(shù)），以Hi≥kHA為濾波準(zhǔn)則（k為常數(shù)）。若k值太小，將無法剔除發(fā)散解區(qū)域；若k值太大，會導(dǎo)致有效場源信息丟失。一般，k值在[0.5，1.5]范圍內(nèi)選擇（姚長利等，2004）。

水平梯度濾波法優(yōu)勢在于：①可剔除水平梯度較小的區(qū)域網(wǎng)格點，選擇水平梯度模相對較大的區(qū)域參與反演；②可去除因測網(wǎng)不均勻插值帶來的虛假網(wǎng)格點信息（姚長利等，2004；管志寧等，2005；鄭秋月等，2021）。文中在計算場源位置時，采用水平梯度濾波法，對傳統(tǒng)的歐拉反演結(jié)果予以優(yōu)化，其中k值取為1。

5 重力場變化及場源特征反演與解釋

5.1 鄂爾多斯南緣地區(qū)“靜態(tài)”重力場變化特征

地震孕育發(fā)生是區(qū)域構(gòu)造應(yīng)力積累與釋放的過程，震源位置和震源機制是巖石圈動力學(xué)過程的反映，通過地震可以了解地下深部的地質(zhì)構(gòu)造和應(yīng)力狀態(tài)，同時可為分析重力異常變化與地震分布的關(guān)系提供便利。為此，統(tǒng)計鄂爾多斯南緣地區(qū)歷史地震目錄，繪制歷史地震分布圖，見圖4。

圖4 鄂爾多斯南緣地區(qū)地震分布Fig.4 Earthquake distribution in the southern margin of the Ordos area

重力位場異常可以表征地質(zhì)構(gòu)造運動過程，常用于劃分構(gòu)造單元、識別斷裂形態(tài)和盆地邊界（馬宗晉等，2006）。一般認為，自由空氣重力異?？梢员碚鞯叵挛镔|(zhì)盈虧程度；布格重力異?？梢苑从成畈繕?gòu)造特征，被稱為構(gòu)造運動的“活化石”（陳石等，2011b）。在分析鄂爾多斯南緣地區(qū)“動態(tài)”重力場變化之前，分析“靜態(tài)”重力場變化特征?；贓GM2008 地球重力模型，計算自由空氣重力異常和布格重力異常（Sandwell et al，2009；Pavlis et al，2012），結(jié)果見圖5。

圖5 鄂爾多斯地塊南緣地區(qū)重力異常（a）自由空氣異常；（b）布格重力異常Fig.5 Gravity anomaly in the southern margin of the Ordos block

圖5(a)為自由空氣重力異常，異常范圍基本在(-140—110)×10-5m·s-2。該異常與地形具有較強的相關(guān)性，在高海拔的山地，重力異常表現(xiàn)為正值，在低海拔的平原和盆地，重力異常表現(xiàn)為負值，自由空氣異?；痉从沉嗽搮^(qū)構(gòu)造單元的大體輪廓。圖5(b)為布格重力異常，研究區(qū)整體表現(xiàn)為負值異常，范圍在(-240—-80)×10-5m·s-2。研究區(qū)西北端靠近青藏高原東北緣地區(qū)，其布格重力異常呈現(xiàn)極小負值，向SE，重力異常呈逐漸增大趨勢，異常整體形態(tài)與斷裂及構(gòu)造區(qū)分布具有一定對應(yīng)關(guān)系。結(jié)合圖4 所示的地震分布，發(fā)現(xiàn)大部分地震震中與斷裂帶發(fā)育位置具有較好的一致性，但有些地震與斷裂之間的對應(yīng)關(guān)系不明顯。此外，重力異常的變化特征與地震分布也具有一定相關(guān)性。進一步說明，重力異常變化、斷裂活動與地震孕育發(fā)生具有較好的協(xié)調(diào)性和一致性。

5.2 “靜態(tài)”重力場歐拉反演與解釋

為了全面掌握鄂爾多斯南緣地區(qū)重力異常變化的場源特征，對該區(qū)布格重力異常進行歐拉反演計算。反演前，對布格重力異常進行坐標(biāo)轉(zhuǎn)換，將大地坐標(biāo)轉(zhuǎn)換為平面直角坐標(biāo)，采用克里格插值算法進行插值，并采用高斯濾波進行平滑。在此基礎(chǔ)上，對布格重力異常數(shù)據(jù)進行歐拉反演，策略如下：①構(gòu)造指數(shù)取1.5；②X、Y方向點距均設(shè)置為10 km；③水平梯度濾波倍數(shù)k=0.5；④滑動窗口設(shè)為4—8 倍點距。由此獲得基于布格重力異常的鄂爾多斯南緣地區(qū)場源特征，見圖6。從“靜態(tài)場”角度分析鄂爾多斯南緣地區(qū)場源特征，由圖6 可知，在4—8 倍的歐拉反演結(jié)果中，場源分布的水平位置和深度基本相同，聚集性相對較好的區(qū)域集中在斷裂邊界帶附近，尤其是研究區(qū)西部的海原—六盤山地區(qū)以及南北地震帶北部地區(qū)，異常源聚集性較好且與斷裂走向基本一致，呈NS 向分布。渭河盆地內(nèi)部場源呈EW 向分布特征，主要沿斷裂走向展布。在鄂爾多斯地塊內(nèi)部出現(xiàn)平?jīng)觥獩艽ê蛻c城—富縣地區(qū)2 個場源特征較顯著區(qū)域。秦嶺造山帶場源集中在勉—略大斷裂、月河斷裂和山陽斷裂附近。從場源深度看，0—30 km 深度的場源占整個反演結(jié)果的90%以上，表明該區(qū)布格重力異常變化主要由淺部異常源引起。結(jié)合圖4 可知，地震分布的三維位置與基于布格重力異常計算的場源位置吻合程度較高。

圖6 基于布格重力異常的鄂爾多斯南緣地區(qū)場源特征Fig.6 Field source characteristics in the southern margin of Ordos based on Bouguer gravity anomaly

5.3 鄂爾多斯南緣地區(qū)“動態(tài)”重力場變化特征

選取研究區(qū)2014—2017 年相對重力數(shù)據(jù)，計算“動態(tài)”重力異常變化，結(jié)果見圖7。其中(a)—(c)圖分別表示2014年8月—2015年8月、2015年8月—2016年7月、2016年7月—2017 年7 月的1 年尺度重力差分變化結(jié)果，(d)圖表示2014 年8 月—2017 年7 月3 年尺度的重力累積變化結(jié)果。圖中灰色線表示斷裂帶，藍色等值線表示負變化，紅色等值線表示正變化，黑色等值線表示“0”值線，等值線間距10×10-8m·s-2。因部分地區(qū)重力站點分布不均，采用克里格插值算法，對平差后重力數(shù)據(jù)進行插值及高斯濾波，等值線覆蓋范圍與測網(wǎng)范圍基本一致。在測點稀疏地區(qū)，重力場變化的可信度偏低，不做深入分析。

由圖7(a)可見，2014 年8 月—2015 年8 月，該區(qū)重力變化基本以西安為界，以西大部分地區(qū)以負值變化為主，以東地區(qū)以正值變化為主。在成縣—太白斷裂與秦嶺北緣斷裂交匯處出現(xiàn)局部重力正值變化，量值約20×10-8m·s-2。另一個重力正值變化區(qū)位于海原—六盤山斷裂帶以東固原、慶城一帶，量值在30×10-8m·s-2。以渭河斷裂與長安—臨潼斷裂交界為中心出現(xiàn)不規(guī)則重力四象限分布特征，整體來看，重力場展布特征受主要斷裂控制。由圖7(b)可見，2015 年8 月—2016 年7 月，研究區(qū)重力以負值變化為主，在西秦嶺北緣斷裂、六盤山—寶雞斷裂帶附近出現(xiàn)重力變化梯度帶，差異變化達70×10-8m·s-2。其他地區(qū)重力場變化相對平穩(wěn)，變化量基本在(-40— -10)×10-8m·s-2。由圖7(c)可見，2016 年7 月—2017 年7 月重力場變化相對復(fù)雜。與圖6(a)相比，重力場出現(xiàn)反向變化，基本以西安為界，以西地區(qū)重力表現(xiàn)為大范圍正值變化，而東部則為負值變化。在六盤山—寶雞斷裂附近不僅存在重力變化梯度帶，還具有四象限分布特征。在秦嶺中部及鄂爾多斯東南部地區(qū)也出現(xiàn)了重力變化梯度帶。由圖7(d)可見，2014 年8 月—2017年7 月，研究區(qū)重力變化相對平穩(wěn)，且以負值變化為主，變化量基本在(-40— -10)×10-8m·s-2，僅在海原—六盤山斷裂以西地區(qū)出現(xiàn)局部正值變化，變化量最大為50×10-8m·s-2，可見3 年尺度重力累積變化不明顯，未出現(xiàn)大震、中強震的孕育特征。

由圖7 可見，1 年尺度的重力場變化在某些區(qū)域存在一定繼承性，但大部分地區(qū)的重力變化存在明顯差異。海原—六盤山—隴縣斷裂帶附近重力變化梯度帶位置在不同時間尺度上基本一致，但是變化量存在差異，重力變化等值線的展布和斷裂分布具有一定對應(yīng)關(guān)系。研究表明：半年尺度，甚至1 年尺度的重力變化依然存在正負波動的變化特征。這種波動可能由近地表局部影響、觀測誤差所導(dǎo)致，目前尚無法有效分離各種場源引起的重力變化。理論上，若為深部場源引起的區(qū)域重力變化，其在不同時間尺度內(nèi)應(yīng)具有相似性和一致性，此區(qū)域即為比較可靠的異常源位置，也是物質(zhì)匯集、應(yīng)力積聚區(qū)域，未來發(fā)生地震的可能性相對較大。

5.4 “動態(tài)”重力場歐拉反演結(jié)果與解釋

基于常規(guī)歐拉反褶積方法，對不同時間尺度的重力變化進行反演，尋找聚集性好、場源重復(fù)性好的區(qū)域。此類區(qū)域往往物質(zhì)積累相對集中、構(gòu)造活動頻繁。實際反演時設(shè)定網(wǎng)格化間距為10 km，構(gòu)造指數(shù)為1，窗口大小為5—8 倍測點間距。反演完成后，對反演結(jié)果進行統(tǒng)計，發(fā)現(xiàn)反演解個數(shù)存在差異，具體結(jié)果見表2。

表2 改正前后基于不同滑動窗口獲得的反演結(jié)果個數(shù)Table 2 The number of inversion results obtained based on different sliding windows before and after correction

對不同反演結(jié)果進行對比分析，將6 倍測點間距獲得的結(jié)果作為最優(yōu)解來進行分析，具體反演結(jié)果見圖8(a)。圖中圓圈表示反演的場源位置，圓圈大小與場源深度相對應(yīng)，采用紅、黃、綠、藍4 種顏色表示2014—2017 年不同時間段的反演結(jié)果。由圖7 可知，從引起重力場變化的不同時間尺度結(jié)果來看，場源位置相對分散，無明顯規(guī)律性。總體而言，研究區(qū)邊緣地區(qū)場源深度反演結(jié)果相對較深，內(nèi)部區(qū)域場源則相對較淺。在不同時間尺度下，海原—六盤山—隴縣斷裂帶、成縣—太白斷裂、華山山前斷裂以及韓城斷裂附近地區(qū)場源位置集聚性較高，而鄂爾多斯地塊南部及秦嶺造山帶南部場源位置相對分散。作為一種地震前異常信號，場源位置的分散型分布特征表明，鄂爾多斯南緣地區(qū)近年來不具備中強地震孕育與發(fā)生的條件，同時說明，2014—2017 年重力場變化屬于一種隨機的、正常的波動變化，地震緊迫性相對較弱。此外，統(tǒng)計得到不同深度的反演解個數(shù)，具體結(jié)果見表3。

表3 改正前后不同深度的反演結(jié)果個數(shù)Table 3 The number of inversion results at different depths before and after correction

由表3 可知，2014—2017 年，發(fā)生重力變化的場源深度集中在0—50 km，其中場源深度為5—30 km 的反演結(jié)果約占60%。對比發(fā)現(xiàn)，2014 年8 月—2015 年8 月，場源深度0—10 km 的反演解數(shù)量較少，深度50—100 km 的反演解數(shù)量較多，說明該期重力場變化受深部物質(zhì)影響較大。在該時間范圍之前曾于2013 年7 月22 日發(fā)生岷縣—漳縣MS6.6 地震，該地震距研究區(qū)西部僅150 km，因此重力變化可能受到地震遠程應(yīng)力調(diào)整及物質(zhì)運移的影響。其他3 期的反演結(jié)果顯示，大部分場源深度在30 km 以上，與該區(qū)震源深度基本吻合。

利用常規(guī)歐拉反演方法可獲得比較合理的場源特征，但反演結(jié)果發(fā)散性嚴(yán)重，采用水平梯度濾波法對反演結(jié)果進行優(yōu)化，發(fā)現(xiàn)歐拉反演結(jié)果更加聚集收斂，發(fā)散現(xiàn)象顯著減弱，場源特征更加明顯，具體結(jié)果見圖8(b)，與圖8(a)對比可知，研究區(qū)邊緣大部分區(qū)域反演結(jié)果數(shù)量明顯減少。這是因為，水平梯度法刪除了測點稀疏造成的不真實反演結(jié)果?？傮w來看，研究區(qū)西北部的海原—六盤山地區(qū)出現(xiàn)場源叢集現(xiàn)象，鄂爾多斯地塊內(nèi)部場源集聚性較好地區(qū)分布在慶城一帶。渭河斷陷盆地內(nèi)部反演結(jié)果集聚性較好區(qū)域基本分布在主要斷裂附近，且場源展布與斷裂走向基本一致。秦嶺造山帶的場源特征相對分散，且場源位置相對較深。同時，統(tǒng)計水平梯度濾波反演結(jié)果，與常規(guī)歐拉反演結(jié)果（表2，表3）進行對比，可知：從5—8 倍測點間距滑動窗口結(jié)果看，采用水平梯度濾波法效果較好，反演解數(shù)量較常規(guī)歐拉反演結(jié)果至少減少50%，剔除了大部分不真實的場源信息；采用水平梯度濾波法處理后，對場源深度反演結(jié)果改正效果顯著，進一步表明相對重力觀測結(jié)果主要反映深度在30 km 以內(nèi)的異常源信息。

圖8 鄂爾多斯南緣地區(qū)2014—2017 年重力變化場源特征(a)常規(guī)歐拉反演結(jié)果；(b)水平梯度濾波結(jié)果Fig.8 Field source characteristics of gravity changes in the southern margin of Ordos from 2014 to 2017

5.5 重力變化場源特征與地震的關(guān)系

鄂爾多斯南緣地區(qū)在地質(zhì)構(gòu)造劃分上主要包含鄂爾多斯地塊南部、渭河斷陷盆地及秦嶺造山帶，在地震帶劃分上主要包含汾渭地震帶的渭河地震帶及南北地震帶北部。該研究區(qū)歷史上曾發(fā)生多次6 級以上強震，具備強震發(fā)生的歷史背景和構(gòu)造條件。據(jù)歷史地震記錄和震源機制解結(jié)果，鄂爾多斯南緣地區(qū)震源深度相對較淺，多在30 km 以內(nèi)（王衛(wèi)東，2009；杜建軍，2016）。對圖4、圖6、圖8 及表2 進行細致分析，發(fā)現(xiàn)基于布格重力異常數(shù)據(jù)計算的“靜態(tài)”場源位置和基于相對重力觀測數(shù)據(jù)計算的“動態(tài)”場源位置，均與研究區(qū)歷史地震分布位置具有高度一致性，表明場源集聚性較好的地區(qū)有利于地震的孕育與發(fā)生（劉芳等，2016；鄭秋月等，2021）。本研究結(jié)果表明，海原—六盤山斷裂帶、隴縣—寶雞斷裂帶、渭河斷裂帶西段以及華山山前斷裂東段場源積聚性較好，未來發(fā)生地震的可能性較高。其他構(gòu)造區(qū)異常源呈現(xiàn)相對性和隨機性，地震危險性相對較弱。此外，據(jù)歐拉反演深度，重力場源集中在30 km 以內(nèi)，位于中上地殼，與王衛(wèi)東（2009）、杜建軍（2016）等研究所得震源深度一致。

6 結(jié)論

作為我國東、西部不同構(gòu)造單元的連接部位，鄂爾多斯南緣地區(qū)不僅是物質(zhì)流通的重要通道，也是青藏塊體和華北塊體相互作用的關(guān)鍵部位。受2 種不同構(gòu)造應(yīng)力的影響，該區(qū)山地隆升，盆地下陷，斷裂發(fā)育，地震頻發(fā)，歷史上發(fā)生多次6 級以上強震，甚至8 級巨震。因此，鄂爾多斯南緣地區(qū)的地震危險性歷來是地學(xué)研究的重點關(guān)注問題之一。近年來，隨著地震觀測手段的日益豐富、觀測技術(shù)的不斷提高，對孕震機理、發(fā)震構(gòu)造、前兆信號的研究日趨成熟，并取得長足進展（陳石等，2011c，2014，2015；祝意青等，2017，2020），如基于重力觀測資料開展地震監(jiān)測預(yù)測，取得豐碩成果（祝意青等，2009；申重陽等，2010；Chen et al，2016a，2016b）。但是，大多數(shù)工作是基于重力異?！皥觥遍_展的，針對重力異?！霸础遍_展的工作相對較少，主要原因是缺少足夠的先驗信息。為此，本研究嘗試采用在地球物理勘探工作中應(yīng)用廣泛的三維反褶積方法，進行重力變化場源特征研究，取得以下認識。

（1）利用歐拉反褶積方法進行反演，關(guān)鍵在于選擇合適的構(gòu)造指數(shù)、滑動窗口大小等參數(shù)。通過構(gòu)建模型，正演得到基本可以反映真實重力變化的模型參數(shù)。通過不同實驗方案進行測試，發(fā)現(xiàn)當(dāng)構(gòu)造指數(shù)為1，滑動窗口為4—8 倍點位間距時，反演效果最佳。

（2）反演結(jié)果表明：海原—六盤山—隴縣斷裂、隴縣—寶雞斷裂、渭河斷裂和華山山前斷裂附近，在不同時間尺度下場源特征一致性較好。結(jié)合圖4，可知以上地區(qū)也是中小地震頻發(fā)區(qū)域，即更易于應(yīng)力積累、能量積聚而發(fā)生地震的位置。由圖8、表3 可知，場源反演深度集中在5—30 km 范圍內(nèi)，屬中上地殼。據(jù)地震震源機制解，可知該深度范圍也是鄂爾多斯南緣地區(qū)的主要震源深度。

（3）歐拉反褶積方法適用于反演重力場變化異常的場源信息，不同的場源位置及深度揭示了不同的構(gòu)造信息，多期重力變化反演三維場源特征一致性較好，表明該區(qū)存在物質(zhì)積聚和應(yīng)力增強；反之，則該區(qū)地震危險性較弱。圖8(a)顯示，在不同時間尺度下，場源位置相對分散，表明大部分地區(qū)的重力變化是隨機的物質(zhì)運動現(xiàn)象，屬于較短期的地殼活動現(xiàn)象或者與觀測誤差不確定性有關(guān)。但圖8(b)顯示，經(jīng)水平梯度濾波的歐拉反演結(jié)果聚集性和收斂性更強，場源特征更加明顯。因此，在實際反演時，需要對歐拉反演結(jié)果做進一步濾波處理。

需要說明的是，選擇2014—2017 年相對重力觀測數(shù)據(jù)開展研究，主要是因為，2013 年陜西重力測網(wǎng)進行了優(yōu)化改造，增建大量站點，而2018 年進行了二次優(yōu)化改造，刪減了大量站點。因此，陜西重力測網(wǎng)在該時間范圍內(nèi)點位密度相對較高，分布較均勻。此外，本研究在建立正演模型時，僅考慮了直立六面體組合，實際地質(zhì)構(gòu)造更復(fù)雜，而文中測試模型相對簡單。盡管如此，本研究結(jié)果仍能驗證歐拉反褶積方法的可靠性和適用性。希望通過本研究，為鄂爾多斯南緣地區(qū)地震危險性研判提供一種新的研究思路和技術(shù)途徑。

文中流動重力數(shù)據(jù)由國家重力臺網(wǎng)中心提供，撰寫過程中多次與云南省地震局鄭秋月工程師進行討論，匿名審稿專家及編輯部老師對本文提出的修改意見與建議，對文章質(zhì)量的提升大有益處，在此一并表示感謝！