董澤義 湯 吉* 趙國澤 陳小斌 崔騰發(fā) 韓 冰 姜 峰 王立鳳

1)中國地震局地質(zhì)研究所，地震動力學國家重點實驗室，北京 100029

2)應急管理部國家自然災害防治研究院，北京 100085

3)中國地震局地震預測研究所，北京 100036

4)中國科學院南海海洋研究所，邊緣海與大洋地質(zhì)重點實驗室，廣州 510301

0 引言

地震是對人類生命與財產(chǎn)安全造成嚴重危害的突發(fā)性自然災害之一。20世紀以來，中國共發(fā)生6級以上地震800多次，死于地震的人數(shù)達55萬之多，占全球因地震死亡人數(shù)的50%以上。顯然，開展地震預測預報研究既是科學發(fā)展的需要，也是人類社會可持續(xù)發(fā)展的需要(陳颙等，2003；陳運泰，2009)。

由于地震發(fā)生的機制十分復雜，地震預測預報問題迄今尚未得到解決，仍然是世界重要的前沿科學技術(shù)難題之一。盡管如此，世界各國特別是地震災害嚴重的國家，如中國、美國、日本、意大利、希臘、俄羅斯等針對地震預測預報投入了大量人力物力，并開展了長期的探索研究，力求減輕地震災害(趙國澤等，2012)。經(jīng)過上百年的努力探索，全球建立了大量地震監(jiān)測臺站，積累了大量觀測資料，也取得了一些重要進展。大量觀測資料表明，強震前的電磁異?，F(xiàn)象確實是存在的，在眾多地震預報預測研究手段中，電磁觀測被認為是最有可能首先取得地震預測突破性進展的方法之一(Bleieretal.，2005；Zhaoetal.，2009；Uyeda，2015)，其已成為中國地震短臨預測的最有效手段之一。

中國從1966年邢臺地震開始對地震電磁觀測開展研究，至今已建設(shè)了數(shù)百個地電場、地磁場和地電阻率觀測臺站，得到了不少地震前后電、磁場和電阻率對異常變化信息(趙玉林等，1978；黃清華等，2006；高曙德等，2010)。但這些傳統(tǒng)的電磁法觀測手段主要關(guān)注恒定場成分，采樣率較低(最高采樣頻率為1Hz)，未涉及對地震事件敏感的超低頻、極低頻(SLF、ELF，3～300Hz)及其附近頻帶的電磁場觀測(趙國澤等，2012，2015)。相比傳統(tǒng)的電法或磁法，交變電磁場法在地震觀測中具有自身優(yōu)勢與潛力，也被應用于地震電磁異常的觀測研究，在國內(nèi)外取得了一些重要的研究成果(Reddyetal.，1976；張云琳等，1994；湯吉等，1998)。但大多數(shù)觀測關(guān)注的是天然源電磁場，由于受到日益嚴重的人文干擾，其應用效果也受到了限制。為了提高對各種電磁異常現(xiàn)象的觀測識別和捕捉能力，建立中國地震電磁觀測新技術(shù)系統(tǒng)，人工源極低頻電磁技術(shù)CSELF(Control Source Extremely Low Frequency)被提出和發(fā)展起來，經(jīng)過約十年的研究和試驗證明，CSELF法是當前最有發(fā)展前途的地震短臨前兆監(jiān)測方法之一，它不僅具有抗干擾能力強、信噪比高、覆蓋范圍大、探測深度大、有利于識別和捕捉地震異常信號等許多優(yōu)點，同時其既可用于研究空間電磁場異常，也可用于研究地殼電性結(jié)構(gòu)異常，使觀測成本極大降低(趙國澤等，2003a，b，2007，2010a，2012)。

“十一五”國家重大科學技術(shù)設(shè)施工程項目“極低頻探地(WEM)工程”地震預測分系統(tǒng)在地震活動較強烈的南北地震構(gòu)造帶南段的川滇地區(qū)和地震重點監(jiān)視區(qū)首都圈地區(qū)分別建設(shè)了15個CSELF電磁臺站，組成了中國乃至世界的第1個既觀測天然源電磁場(頻帶為1000～0.001Hz)也觀測人工源電磁場(頻帶為300～0.1Hz)的大陸CSELF電磁臺網(wǎng)(圖1a)。該臺網(wǎng)既能監(jiān)測地震引起的空間電磁場異常，也能監(jiān)測震源區(qū)附近地下電性結(jié)構(gòu)的變化，實現(xiàn)了真正的四維立體電磁異常信息的動態(tài)監(jiān)測，提升了川滇與首都圈地區(qū)中強以上地震電磁異常的捕捉能力，同時也提高了地震預測預報水平，推進了地震預警立體監(jiān)測系統(tǒng)的建設(shè)。

圖1 首都圈極低頻電磁臺網(wǎng)區(qū)的電性結(jié)構(gòu)測量Fig.1 Investigation of electrical resistivity structures beneath the CSELF network in the capital circle region.a 中國大陸極低頻電磁臺網(wǎng)；b 首都圈CSELF臺網(wǎng)區(qū)的地質(zhì)構(gòu)造簡圖；c 大同臺結(jié)構(gòu)測量測點分布

了解CSELF臺網(wǎng)區(qū)和各臺站地下的電性結(jié)構(gòu)背景，是研究相應地區(qū)地震預測監(jiān)測問題的基礎(chǔ)，這不僅對于認識川滇與首都圈2個臺網(wǎng)區(qū)的深部結(jié)構(gòu)特點具有重要意義，而且為識別地震活動前、后地下電性結(jié)構(gòu)可能發(fā)生的動態(tài)變化、理解地震可能發(fā)生破裂的程度和動態(tài)過程、探索地震及地震異常現(xiàn)象可能的成因機制等提供了重要數(shù)據(jù)。該項目首次在川滇和首都圈2個臺網(wǎng)區(qū)開展了深部電性結(jié)構(gòu)探測，通過數(shù)據(jù)采集、處理與分析、反演獲得了2個臺網(wǎng)區(qū)各臺站及附近區(qū)域的地下電性結(jié)構(gòu)。本文將主要介紹如何開展首都圈CSELF臺網(wǎng)區(qū)地下電性結(jié)構(gòu)探測，包括野外數(shù)據(jù)采集、數(shù)據(jù)的維性分析、數(shù)據(jù)反演以及臺網(wǎng)區(qū)的深部電性結(jié)構(gòu)特征分析等方面。文中涉及的圖件較多，由于篇幅有限，只能以部分臺站為例給出示范圖件，特此說明。

1 首都圈CSELF臺網(wǎng)的基本情況與電性結(jié)構(gòu)測量

1.1 臺網(wǎng)的組成與地質(zhì)構(gòu)造背景

中國大陸CSELF臺網(wǎng)依托中國地震局已有30個地震臺站進行建設(shè)。從地理上看，首都圈CSELF臺網(wǎng)基本覆蓋了華北中、東部地區(qū)，由天津市寶坻臺，河北省豐寧臺、懷來臺、文安臺、涉縣臺，山東省無棣臺、萊陽臺、安丘臺、莒縣臺，山西省大同臺、代縣臺、夏縣臺，遼寧省大連臺、營口臺，內(nèi)蒙古呼和浩特臺共15個臺站組成(圖1a，b)。

中國華北是地球上最古老、形成于太古宙時期的克拉通塊體之一(Liuetal.，1992)，它可分為東部華北盆地和西部鄂爾多斯地塊2部分，中間由山西斷陷帶隔開，其南、北部分別為秦嶺-大別造山帶和陰山-燕山造山帶，2條巨型線性構(gòu)造(太行山重力梯度帶與郯廬斷裂帶)穿其而過(圖1b)。已有研究表明，華北東部地區(qū)自顯生宙以來經(jīng)歷了獨特的巖石圈改造和破壞過程，并伴隨著大規(guī)模的構(gòu)造伸展、斷陷盆地的形成和廣泛的巖漿活動(Griffinetal.，1998；Fanetal.，2000；Xu，2001；Gaoetal.，2004)。華北地區(qū)擁有獨特的板內(nèi)地震活動帶，歷史上曾發(fā)生多次7級以上破壞性大地震，包括震級最大的1679年三河-平谷8級地震和震驚世界的1976年7.8級唐山大地震，該區(qū)一直是研究大陸構(gòu)造形變和動力學的天然實驗室。首都圈CSELF電磁臺網(wǎng)的大部分臺站都建立在華北地塊的構(gòu)造邊界和重要構(gòu)造帶上，對于提升整個華北區(qū)域強震的監(jiān)測預報研究水平具有重要作用。

1.2 野外測量與數(shù)據(jù)

CSELF臺網(wǎng)的觀測裝置與大地電磁法(Magnetotelluric sounding，簡稱MT)(陳樂壽等，1990)基本相同，觀測的物理量也與MT相同，不同的是其既可觀測天然源電磁信號，也可觀測人工源電磁信號。為獲得首都圈CSELF臺網(wǎng)區(qū)15個臺站及其附近的深部電性結(jié)構(gòu)特征，我們以每個臺站為中心，根據(jù)臺站的區(qū)域構(gòu)造特征，在其兩側(cè)共布設(shè)4個寬頻帶MT測點，形成1條與區(qū)域構(gòu)造基本垂直的短MT剖面，以獲取穿過臺站的電性結(jié)構(gòu)背景。圖1c 以大同臺為例給出了臺站結(jié)構(gòu)探測中臺站與寬頻帶MT測點的空間位置分布關(guān)系。按照上述方案，在首都圈的15個臺站共布設(shè)60個寬頻帶MT測點，野外數(shù)據(jù)采集工作于2016年完成。事實上，由于首都圈經(jīng)濟與社會發(fā)展太快，臺站附近的電磁干擾較之前也有所增強，為獲得高質(zhì)量的寬頻帶MT數(shù)據(jù)，對許多測點進行了多次重測，因此工作量遠大于60個測點。

大部分臺站及寬頻帶MT測點的數(shù)據(jù)質(zhì)量較好，有效頻點的周期范圍為0.001～2000s，但個別臺站離城市中心太近，受電磁干擾影響較嚴重。作為例子，圖2 給出了大同臺、夏縣臺、大連臺和豐寧臺和附近MT測點的xy與yx向視電阻率與阻抗相位曲線，測點編號均按照從北至南(或從西至東)由小變大。分析曲線發(fā)現(xiàn)，大部分臺站兩側(cè)的曲線形態(tài)發(fā)生了較大變化，可能代表著臺站位于電性邊界帶上，這與前文所述的很多臺站處于區(qū)域構(gòu)造邊界帶上的事實相符。為了更直觀地了解臺網(wǎng)區(qū)不同頻率的視電阻率分布特征，我們對臺網(wǎng)各臺站和所有寬頻帶MT測點數(shù)據(jù)進行插值，給出了首都圈CSELF臺網(wǎng)區(qū)不同頻率的xy、yx向視電阻率及兩者幾何平均的視電阻率平面分布云圖，圖3 給出了頻率為100Hz、0.01Hz的xy、yx向以及兩者幾何平均的視電阻率平面圖，按照趨膚深度的定義，這2個頻率分別代表了淺部與深部的視電阻率分布情況?？梢园l(fā)現(xiàn)，在郯廬斷裂帶、太行山重力梯度帶兩側(cè)，無論高頻還是低頻的視電阻率變化均較為明顯，暗示其不僅是構(gòu)造邊界，同時也是電性邊界；華北北部的陰山-燕山造山帶、東部的華北平原的視電阻率從高頻到低頻均呈現(xiàn)出較好的繼承性。此外，通過長期動態(tài)監(jiān)測臺網(wǎng)區(qū)的視電阻率分布云圖也可觀測臺網(wǎng)區(qū)視電阻率是否發(fā)生異常變化，這可為區(qū)域的地震監(jiān)測預報提供重要參考依據(jù)。

2 區(qū)域構(gòu)造的維性分析

圖3 首都圈CSELF電磁臺網(wǎng)區(qū)不同頻率的視電阻率分布圖Fig.3 Apparent resistivity maps at 100Hz and 0.01Hz of the CSELF network in the capital circle region.ρxy xy向視電阻率；ρyx yx向視電阻率；ρave xy與yx向視電阻率的幾何平均

圖4 不同臺站每個測點的二維偏離度Fig.4 The 2D skewness at each site for Datong，Xiaxian，Dalian and Fengning staions.

圖5 首都圈CSELF臺網(wǎng)各臺站的電性主軸方位玫瑰圖Fig.5 Rose diagrams of geoeletrical strike at each station of the CSELF network in capital circle region.

此外，根據(jù)阻抗張量分解結(jié)果，我們也給出了首都圈CSELF臺站的電性主軸方位統(tǒng)計玫瑰圖，如圖5 所示，該圖較直觀地反映了各臺站地下可能的電性主軸方位。由于玫瑰圖的電性主軸方位角存在90°的模糊性，還需結(jié)合臺站區(qū)域的地質(zhì)構(gòu)造走向確定具體的電性結(jié)構(gòu)走向?？傮w上，各臺站的電性主軸方位與所處地質(zhì)構(gòu)造單元構(gòu)造的走向較一致，且同一構(gòu)造單元的臺站電性主軸方位具有一定的相似性。例如，同處在山西斷陷帶的大同臺、代縣臺和夏縣臺，電性主軸方位基本一致。同樣地，郯廬斷裂帶附近的安丘臺與莒縣臺的電性主軸方位也基本一致。

磁感應矢量是MT數(shù)據(jù)中的一個重要參數(shù)，其定義為垂直磁場與水平磁場的復數(shù)比值，包括實感應矢量與虛感應矢量。實感應矢量主要反映地下介質(zhì)導電性分布的橫向不均勻性，矢量大小反映橫向?qū)щ娦圆町惖拇笮?，在Parkinson規(guī)范下(Parkinson，1962)，矢量方向指向電流會聚的區(qū)域，即高導體所在的位置。此外，在理想的二維條件下，實感應矢量所指方向垂直于構(gòu)造走向，可用于判斷區(qū)域電性主軸方位。我們也分析了各臺站MT剖面數(shù)據(jù)的實感應矢量，定性地了解了各臺站的地下電性結(jié)構(gòu)特征。圖6 以大同臺和豐寧臺為例給出了沿剖面MT測點的實感應矢量隨頻率的分布。結(jié)果顯示，如果用趨膚深度公式進行簡單估算，2個臺站在中上地殼表現(xiàn)出較強的橫向不均勻性，而中下地殼的介質(zhì)橫向相對均勻，且2個臺站的南段在中上地殼呈現(xiàn)相對高導，這與后文反演獲得的電性結(jié)構(gòu)揭示的特征具有一致性。

圖6 大同臺和豐寧臺沿剖面的實感應矢量擬斷面圖Fig.6 Pseudo-sections of the real induced vector along the MT profile at Datong and Fengning stations.

3 臺站及其MT數(shù)據(jù)反演與電性結(jié)構(gòu)

3.1 一維與二維反演

為獲得每個臺站及其附近的深部結(jié)構(gòu)，我們既對臺站的每個測點進行了一維反演，也對臺站及其MT測點組成的短剖面數(shù)據(jù)進行了二維反演。

一維反演研究每個臺站正下方介質(zhì)的電性結(jié)構(gòu)隨深度的變化情況。本文采用自適應正則化反演方法(陳小斌等，2005)對各臺站觀測的xy和yx向視電阻率與相位分別進行了反演，并對2個方向數(shù)據(jù)的幾何平均也進行了反演，最后在每個臺站獲取3條電阻率隨深度的變化曲線。圖7 以大同臺為例，給出了大同臺及其MT測點的數(shù)據(jù)一維反演結(jié)果。我們發(fā)現(xiàn)，在1km以淺深度，3條反演曲線基本重合，表明該深度范圍內(nèi)基本以一維構(gòu)造為主，可能反映大同盆地的沉積層厚度；大同臺與其北側(cè)2個MT測點10km深度以淺的一維結(jié)構(gòu)比較一致，而與其南側(cè)的MT測點的反演結(jié)構(gòu)有明顯差異，可能反映大同臺兩側(cè)的深部電性結(jié)構(gòu)存在較為明顯的橫向不均勻性，這一點在后面的二維反演結(jié)果中得到了證實。

圖7 大同臺及其MT測點的一維反演結(jié)果Fig.7 The 1-D resistiviy structures for Datong station and MT sites.

圖8 大同臺MT剖面的二維反演Fig.8 Two-diminsional inversion of MT profile at Datong station.a 不同光滑因子的L曲線；b 大同臺沿MT剖面的地殼電性結(jié)構(gòu)；c 二維反演數(shù)據(jù)與模型響應擬斷面圖。ρObs觀測視電阻率；ρCalc 模型響應電阻率；φObs 觀測相位；φCalc 模型響應相位

采用目前廣泛使用的非線性共軛梯度(NLCG)法(Rodietal.，2001)進行二維反演，整個反演過程均在MT-Pionner大地電磁可視化集成軟件系統(tǒng)下完成(陳小斌等，2004)。進行二維反演前，根據(jù)前文阻抗張量的分析結(jié)果，將每個臺站的剖面MT數(shù)據(jù)全部旋轉(zhuǎn)至各自的電性主軸方位，旋轉(zhuǎn)后平行于電性主軸方向的數(shù)據(jù)為TE模式，與之垂直方向上的數(shù)據(jù)為TM模式。在對實測數(shù)據(jù)進行二維反演時，不同極化模式數(shù)據(jù)的反演結(jié)果差別較大。針對如何選擇反演數(shù)據(jù)的問題，國內(nèi)外學者已進行過大量研究，一般認為，TE極化模式對深部結(jié)構(gòu)變化更靈敏，TM極化模式對表層結(jié)構(gòu)變化更敏感，采用TE+TM聯(lián)合反演能利用各自的優(yōu)勢，增加數(shù)據(jù)約束量，提高模型的可靠度(Berdichevskyetal.，1998)。但是，當數(shù)據(jù)受三維結(jié)構(gòu)影響明顯時，單獨進行TM模式反演比利用TE模式或TE+TM模式聯(lián)合反演更合理，反演結(jié)果中的虛假構(gòu)造會明顯減少，所得結(jié)果能夠比較準確地反映地下電性結(jié)構(gòu)特征。對于實測數(shù)據(jù)的二維反演，應優(yōu)先考慮采用TM模式數(shù)據(jù)進行二維反演，其次為TM+TE模式，一般不宜采用TE模式(Ledo，2006；蔡軍濤等，2010)。根據(jù)數(shù)據(jù)的定性分析，由于部分臺站數(shù)據(jù)的低頻段受三維效應的影響，我們僅利用每個臺站剖面數(shù)據(jù)的TM單極化模式數(shù)據(jù)進行二維反演，以獲得沿剖面的地下電性結(jié)構(gòu)。

為獲得可靠的二維反演模型，我們通過改變不同的反演參數(shù)進行了大量反演計算，如改變正則化因子，在0.1～1000間取十幾個值進行多次反演，通過L曲線分析確定最優(yōu)的正則化因子(圖8a)；調(diào)整視電阻率與相位的誤差門檻值；將數(shù)據(jù)旋轉(zhuǎn)不同的角度等。所有反演采用的初始模型為100Ω·m的均勻半空間，利用測點中心網(wǎng)格的自動生成技術(shù)構(gòu)建帶地形的初始模型網(wǎng)格(陳小斌等，2009)。按照上述反演方案，對首都圈15個臺站的剖面數(shù)據(jù)進行了大量反演計算，反演工作量較大。最終獲得了每個臺站及其附近地下穩(wěn)定可靠的二維電性結(jié)構(gòu)模型。

3.2 臺站地下二維電性結(jié)構(gòu)

反演得到的15個臺站的地下電性結(jié)構(gòu)模型較穩(wěn)健，能較好地描述各臺站附近地殼的結(jié)構(gòu)特征，并與區(qū)域地質(zhì)構(gòu)造具有一定的對應關(guān)系。張繼紅等(2019)對山東省無棣臺、安丘臺和莒縣臺附近區(qū)域的地下二維電性結(jié)構(gòu)模型進行了解譯。無棣臺的區(qū)域地殼電性結(jié)構(gòu)相對簡單，具有成層性，符合穩(wěn)定地塊區(qū)的電性結(jié)構(gòu)特征(趙國澤等，2010b)；而安丘臺和莒縣臺區(qū)域的地殼電性結(jié)構(gòu)相對復雜，呈現(xiàn)出活動地塊邊界帶的結(jié)構(gòu)特點，也反映了郯廬斷裂中段的活動性和電性特點。本文將以大同臺為例，簡要分析解譯該臺站附近的地下二維電阻率模型。

圖9 首都圈CSELF電磁臺網(wǎng)區(qū)二維反演模型不同深度的電阻率切片F(xiàn)ig.9 Horizontal slices of 2-D electrical resistivity structure of the CSELF network in capital circle region at different depths.

為了更好地了解臺網(wǎng)區(qū)電阻率的空間分布特征，根據(jù)各臺站剖面數(shù)據(jù)的二維反演結(jié)果，通過插值獲得了首都圈CSELF臺網(wǎng)區(qū)不同深部的電阻率分布特征。圖9 分別給出了深度為1km、5km、10km、15km、20km和30km的電性結(jié)構(gòu)平面圖。將其與臺網(wǎng)區(qū)視電阻率頻率響應結(jié)構(gòu)(圖3)進行對比發(fā)現(xiàn)，反演獲得的電性結(jié)構(gòu)與視電阻率響應結(jié)構(gòu)的基本特征大致相符。華北平原從淺至深，地殼均表現(xiàn)為低電阻率特征，而北部的陰山-燕山造山帶、東部的遼東-膠東地塊的地殼表現(xiàn)為相對高阻。值得注意的是，西部山西斷陷帶的地殼從北部的大同盆地到南部的夏縣盆地均表現(xiàn)為相對低阻的特征。電性特征的差異與不同構(gòu)造區(qū)的地形、地貌特征及表層出露巖石的電性具有較好的對應關(guān)系。華北平原地勢較平坦，沉積層比較發(fā)育，而陰山-燕山、遼東-膠東地區(qū)多為山區(qū)或丘陵地貌，主要出露變質(zhì)結(jié)晶巖和火成巖，按照巖石導電性的差異，沉積層的電阻率往往比變質(zhì)巖和火成巖低得多。橫向上，1～15km深度的電阻率橫向差異較大；15km深度以深處，電阻率的橫向差異小于淺部。太行山重力梯度帶作為華北東部地塊與西部地塊的分隔帶，在電性結(jié)構(gòu)上表現(xiàn)為電性邊界帶，兩側(cè)的電阻率結(jié)構(gòu)差異明顯，西側(cè)的電阻率小于其東側(cè)的電阻率。同樣地，郯廬斷裂帶作為中國東部的一條巨型走滑斷裂，其中段穿過了首都圈CSELF臺網(wǎng)區(qū)，兩側(cè)的地殼電性結(jié)構(gòu)也差異明顯，西側(cè)低阻、東側(cè)高阻，且切割深度較深，表現(xiàn)為超殼大斷裂。前人在首都圈CSELF臺網(wǎng)區(qū)也開展了不少大地電磁研究，如魏文博等(2006，2008)完成了應縣—商河(HB-MT01)大地電磁探測；詹艷等(2011)在河北石家莊地區(qū)開展了大地電磁探測。這些工作揭示太行山前斷裂表現(xiàn)為較大的電性邊界帶，西側(cè)太行山隆起地殼結(jié)構(gòu)為高阻區(qū)，而東側(cè)的華北平原為低阻區(qū)。穿過郯廬斷裂中段的大地電磁剖面(葉高峰等，2009；張繼紅等，2010)也反映郯廬斷裂帶東、西兩側(cè)地殼的電性結(jié)構(gòu)特征差異明顯，東低西高。以往的認識與本研究的結(jié)果基本一致。此外，華北大地電磁測深陣列觀測實驗研究所揭示的華北地區(qū)的地殼電性結(jié)構(gòu)特征與本研究的結(jié)果也有較好的一致性(Ouro-Djobo等，2018)。由于臺站比較稀疏，臺站之間的電性結(jié)構(gòu)主要通過插值獲得，這里不做過多討論。

4 臺網(wǎng)區(qū)的三維電性結(jié)構(gòu)特征

近年來，隨著計算機硬件和MT理論的不斷發(fā)展，大地電磁三維正、反演技術(shù)已逐漸走向成熟，基于數(shù)據(jù)空間OCCAM(Siripunvarapornetal.，2005)、非線性共軛梯度(NLCG)(Newmanetal.，2000)、自適應正則因子的擬牛頓算法(AR-QN)(阮帥，2015)等三維MT正、反演程序包均已實現(xiàn)，并被廣泛應用于實測資料的反演解釋中(Xiaoetal.，2015，2017；Caietal.，2017；阮帥等，2020)。

我們利用極低頻地震預測分系統(tǒng)項目研發(fā)的“大尺度三維電性結(jié)構(gòu)快速反演軟件”(toPeak)，采用基于NLCG算法的ModEM并行程序包(Kelbertetal.，2014)對首都圈CSELF臺網(wǎng)區(qū)的結(jié)構(gòu)探測數(shù)據(jù)進行了三維反演。由于臺站間距較大，考慮到網(wǎng)格剖分的問題，對每個臺站僅選擇一個測點參與反演。反演輸入數(shù)據(jù)為15個測點在1000～0.0001Hz頻段范圍內(nèi)所有49個頻點的xy和yx向視電阻率和相位，視電阻率和相位數(shù)據(jù)的門限誤差分別為5%與1.43°，自動構(gòu)建的反演網(wǎng)格有55(NS向)×70(EW向)×85(垂直方向，包括7個空氣層)個，反演初始模型采用100Ω·m 的均勻半空間。經(jīng)過100次迭代，最終模型的數(shù)據(jù)擬合均方差RMS為2.36。圖10 給出了首都圈CSELF臺網(wǎng)區(qū)三維電性結(jié)構(gòu)不同深度的平面圖。與臺網(wǎng)區(qū)二維反演電阻率平面分布相比(圖9)，三維反演呈現(xiàn)的臺網(wǎng)區(qū)不同深度的電阻率空間分布特征與之較為相似，差異主要存在于臺站間的結(jié)構(gòu)特征，這主要是由于臺站間距較大，臺站間沒有數(shù)據(jù)約束所致，因此三維反演結(jié)果中臺站間的電阻率都趨于初始模型的值。而二維反演結(jié)果中臺站間的電阻率通過插值運算獲得，雖然信息更為豐富，但可信度不高。

由于參與三維反演的數(shù)據(jù)有限，本文給出的首都圈CSELF電磁臺網(wǎng)區(qū)三維電性結(jié)構(gòu)背景是比較粗略的，對其結(jié)構(gòu)特征不做細致討論，但該模型的獲取為首都圈CSELF電磁臺網(wǎng)在地震監(jiān)測預報中發(fā)揮作用提供了重要的基礎(chǔ)資料?？偠灾?，通過對首都圈臺網(wǎng)區(qū)觀測數(shù)據(jù)的一維、二維和三維反演研究，本文揭示了各臺站及其附近的地殼電性結(jié)構(gòu)和整個臺網(wǎng)區(qū)的電阻率空間分布特征，通過與前人的研究結(jié)果進行對比與分析，本研究的結(jié)果具有較高的可靠度。

圖10 首都圈CSELF電磁臺網(wǎng)三維反演模型不同深度的切片F(xiàn)ig.10 Horizontal slices of 3-D electrical resistivity structure of the CSELF network in capital circle region at different depths.

5 結(jié)論

本研究通過野外數(shù)據(jù)采集、處理與分析，一維、二維和三維反演，獲得了可靠的首都圈CSELF電磁臺網(wǎng)區(qū)各臺站及其附近的地殼電性結(jié)構(gòu)以及整個臺網(wǎng)區(qū)不同深度的電性結(jié)構(gòu)。結(jié)果表明，電性結(jié)構(gòu)與臺網(wǎng)區(qū)的構(gòu)造具有很好的對應性，如北部陰山-燕山造山帶、東部膠-遼地塊以及中部太行山地區(qū)的臺站地殼結(jié)構(gòu)主要表現(xiàn)為高阻特征，而華北平原、山西斷陷盆地區(qū)的臺站地殼電性結(jié)構(gòu)以高導特征為主。顯然，各臺站的電性結(jié)構(gòu)特征與所處地區(qū)的地形、地貌以及區(qū)域地質(zhì)構(gòu)造特征具有較好的關(guān)聯(lián)性。整個首都圈臺網(wǎng)區(qū)的電阻率空間分布特征展示了臺網(wǎng)區(qū)的中上地殼具有較強的橫向不均勻性，特別是在太行山重力梯度帶和郯廬斷裂帶兩側(cè)，這意味著它們不僅是重要的構(gòu)造邊界帶，也是電性邊界帶。

首都圈CSELF電磁臺網(wǎng)的大部分臺站不僅位于重要的地質(zhì)構(gòu)造帶上，同時也位于地震活動帶上，如郯廬斷裂帶、張家口—渤海地震帶以及山西裂谷地震帶等。這些地區(qū)的臺站區(qū)域地殼電性結(jié)構(gòu)相對復雜，但基本具有上地殼多表現(xiàn)為高阻而中下地殼發(fā)育高導層的特征，并且區(qū)域的歷史地震震中多位于中上地殼低阻與高阻交界處的高阻一側(cè)，地震的空間分布與地下結(jié)構(gòu)具有較大關(guān)系。由于文章篇幅有限，本文僅對大同臺站及附近的地殼電性結(jié)構(gòu)進行了解譯，并討論了大同-陽原地震的深部孕震結(jié)構(gòu)。

首都圈CSELF電磁臺網(wǎng)區(qū)的地下電性結(jié)構(gòu)探測提供了整個區(qū)域的電性結(jié)構(gòu)背景，對地震活動區(qū)的孕震環(huán)境、地震電磁異常現(xiàn)象產(chǎn)生的機理以及地震預測預報研究具有重要作用。未來，我們應長期不斷地監(jiān)測臺網(wǎng)區(qū)的地下電性結(jié)構(gòu)變化特征，真正地讓它在防震減災工作中發(fā)揮作用。

致謝野外工作人員在數(shù)據(jù)采集過程中的辛勤付出為本研究提供了基礎(chǔ)數(shù)據(jù)；中國地震局地質(zhì)研究所信息中心為大地電磁三維反演工作提供了超算平臺。在此一并表示感謝！