張家口—渤海地震帶天津—塘沽段的電性結(jié)構(gòu)及其意義

高武平 林逸 張文朋 徐新學(xué) 閆成國(guó) 曹井泉 劉紅艷

摘要：

華北北部的張家口—渤海地震帶是一條新生地震斷裂帶。位于其中段的薊運(yùn)河斷裂構(gòu)成了唐山菱形塊體的西邊界，一直被認(rèn)為與1979年寧河M6.9地震密切相關(guān)。然而，該區(qū)域的深部結(jié)構(gòu)與孕震動(dòng)力學(xué)環(huán)境尚未得到有效揭示。大港—玉田大地電磁測(cè)深剖面直接穿越寧河M6.9地震的震中區(qū)，大致與張家口—渤海斷裂帶相垂直。探測(cè)結(jié)果顯示，區(qū)域電性結(jié)構(gòu)具有顯著的“橫向分塊、縱向分層”特征。大地電磁測(cè)深剖面還揭示了張家口—渤海斷裂帶中段的海河斷裂、漢沽斷裂、薊運(yùn)河斷裂等區(qū)域性斷裂的深部構(gòu)造特征，且漢沽斷裂可能為切穿地殼的深斷裂。深部低阻體邊緣的高低阻梯度帶與寧河M6.9地震和漢沽兩次M6.2地震的震中具有良好的對(duì)應(yīng)關(guān)系，揭示了研究區(qū)中強(qiáng)地震的發(fā)生與深部顯著高低阻梯度帶的發(fā)育密切相關(guān)。

關(guān)鍵詞：

大地電磁測(cè)深; 電性結(jié)構(gòu); 張家口—渤海地震帶; 薊運(yùn)河斷裂

中圖分類號(hào)： P315????? 文獻(xiàn)標(biāo)志碼：A?? 文章編號(hào)： 1000-0844（2023）02-0410-11

DOI：10.20000/j.1000-0844.20211130003

Electrical structure of the Tianjin—Tanggu segment of the

Zhangjiakou—Bohai seismic belt and its importance

GAO Wuping1， LIN Yi1， ZHANG Wenpeng1， XU Xinxue2， YAN Chengguo1， CAO Jingquan1， LIU Hongyan1

（1. Tianjin Earthquake Agency， Tianjin 300201， China;

2. Tianjin North China Geological Exploration Bureau， Tianjin 300170， China）

Abstract：

The Zhangjiakou-Bohai seismic belt is a seismic fault zone in the northern part of North China. The Jiyunhe fault zone in its middle segment constitutes the western boundary of the Tangshan rhombus block， and it is considered the seismogenic fault of the Ninghe M6.9 earthquake in 1979. To effectively reveal the deep structure and seismogenic dynamic environment of this region， we conducted a magnetotelluric sounding profile in the Dagang—Yutian area， which directly crosses the epicenter of the Ninghe M6.9 earthquake and is almost perpendicular to the Zhangjiakou—Bohai fault zone. Detection results show that the electrical structure in the study area is layered vertically and partitioned horizontally. Additionally， the magnetotelluric sounding profile reveals the deep structural characteristics of some regional faults， i.e.， the Haihe fault， Hangu fault， and Jiyunhe fault in the middle section of the Zhangjiakou—Bohai fault zone， in which the Hangu fault may be a deep fault that runs through the crust. The good correlation between the high and low resistivity gradient zones and the epicenters of the Ninghe M6.9 earthquake and the other two M6.2 earthquakes in Hangu suggests that the moderate-strong earthquakes in the study area are closely related to the development of obvious high and low resistivity gradient zones in the deep area.

Keywords：

magnetotelluric sounding; electrical structure; Zhangjiakou—Bohai seismic belt; Jiyunhe fault

0 引言

張家口—渤海地震帶（后文簡(jiǎn)稱張渤帶）西起太行山北端，橫貫燕山與華北盆地，東入渤海，是中國(guó)東部重要的NW向地震活動(dòng)帶［1］。帶內(nèi)發(fā)育有20多條傾向不一的NW向斷裂，單條長(zhǎng)約40～60 km，大多斷裂同時(shí)具有左旋走滑和正斷活動(dòng)分量，是一條典型的第四紀(jì)活動(dòng)斷裂帶［2］。徐杰等［1］根據(jù)張渤帶與NE向斷裂帶的交切關(guān)系將其劃分為張北—懷來(lái)、南口—三河、天津—塘沽、渤海中部和蓬萊—煙臺(tái)5個(gè)構(gòu)造交匯段（圖1），并指出張渤帶中的大部分6級(jí)以上強(qiáng)震均發(fā)生在與NE向斷裂帶交匯部位。

薊運(yùn)河斷裂、漢沽斷裂與海河斷裂是天津—塘沽段十分重要的三條斷裂（圖1）。薊運(yùn)河斷裂一直被認(rèn)為與1979年寧河M6.9地震的發(fā)生密切相關(guān)［3］;而漢沽斷裂則被認(rèn)為與1976—1977年漢沽附近的兩次M6.2地震有關(guān)聯(lián)［4］;在海河斷裂與滄東斷裂交匯部位也曾在1815年發(fā)生葛沽M5地震［5］。這些地震與這幾條斷裂具有怎樣的聯(lián)系？孕震的深部背景和動(dòng)力學(xué)機(jī)制是什么？這些問(wèn)題并未得到有效的揭示。20世紀(jì)80年代以來(lái)，多條人工地震反射、折射剖面初步揭示了本地區(qū)的地殼基本輪廓［6］。近年來(lái)，亦有不少學(xué)者通過(guò)布格重力和地震數(shù)據(jù)等資料反演分析華北地區(qū)地殼結(jié)構(gòu)特征［7-8］，為更加深入地了解本地區(qū)地殼速度結(jié)構(gòu)、低速層分布及深度、重要斷裂位置與展布、中強(qiáng)地震孕震環(huán)境等提供了重要的數(shù)據(jù)支持。但前人研究對(duì)于本地區(qū)深部構(gòu)造細(xì)節(jié)的刻畫還不夠精細(xì)，無(wú)法為深部孕震動(dòng)力學(xué)環(huán)境分析提供更多參考。

大地電磁測(cè)深（magnetotelluric sounding，簡(jiǎn)稱MT）是研究地球深部結(jié)構(gòu)的有效手段之一，對(duì)刻畫巖石圈和深部熱結(jié)構(gòu)特征有獨(dú)到之處。電性結(jié)構(gòu)的差異歸根到底是由地質(zhì)體的物質(zhì)組成、結(jié)構(gòu)、構(gòu)造、流體含量、熱狀態(tài)等諸多因素引起的。通過(guò)對(duì)電性結(jié)構(gòu)剖面的分析，可以較為客觀地獲得相應(yīng)的地質(zhì)信息。大量研究表明，深大斷裂往往構(gòu)成構(gòu)造單元的分界線，不同構(gòu)造單元往往具有不同的電性結(jié)構(gòu)特征［10-11］。為獲取張渤帶天津段的精細(xì)殼幔電性結(jié)構(gòu)、尋找不同構(gòu)造單元間的結(jié)構(gòu)差異、分析區(qū)域深部斷裂特征及地震孕震環(huán)境，本研究擬通過(guò)大地電磁測(cè)深勘探及相應(yīng)的數(shù)據(jù)處理與反演以獲取電性結(jié)構(gòu)剖面，結(jié)合已有的地質(zhì)與地球物理等資料成果開展對(duì)比分析，以期為探索華北平原北部地區(qū)的深部構(gòu)造背景提供新的證據(jù)。

1 大地電磁測(cè)深剖面

1.1 剖面位置與數(shù)據(jù)采集

天津市地震局在“十二五”期間沿大港—玉田一線開展了長(zhǎng)約135 km的大地電磁測(cè)深勘探，剖面走向約30°，與海河斷裂、漢沽斷裂、薊運(yùn)河斷裂等重要斷裂近垂直交切。構(gòu)造上跨燕山隆起和華北斷陷盆地2個(gè)二級(jí)構(gòu)造單元，穿越薊寶褶皺帶、滄縣隆起和黃驊坳陷3個(gè)三級(jí)構(gòu)造單元，并穿過(guò)1976年寧河6.9級(jí)地震及漢沽兩次6.2級(jí)地震震源區(qū)（圖2）。剖面共布置90個(gè)測(cè)點(diǎn)，平均點(diǎn)距1.5 km，并在薊運(yùn)河斷裂兩側(cè)進(jìn)行了加密。其中長(zhǎng)周期有效觀測(cè)點(diǎn)30個(gè)，寬頻大地電磁測(cè)點(diǎn)60個(gè)。

本研究采用具有GPS同步功能的鳳凰系列MTU-5A型/V8型MT儀器系統(tǒng)進(jìn)行四分量張量觀測(cè)。通過(guò)開工極距、采集時(shí)長(zhǎng)等試驗(yàn)，確定極距范圍為80～100 m，普通MT測(cè)點(diǎn)數(shù)據(jù)采集時(shí)間為20 h，長(zhǎng)周期MT測(cè)點(diǎn)為70 h。由于數(shù)據(jù)采集使用了4臺(tái)V5-2000、V8儀器及探頭，為保證數(shù)據(jù)可靠性，在觀測(cè)前后分別進(jìn)行了儀器標(biāo)定和一致性對(duì)比試驗(yàn)。相鄰兩次標(biāo)定相對(duì)誤差均小于2%，且各套儀器觀測(cè)計(jì)算同一極化模式的視電阻率和相位的相對(duì)均方誤差也小于5%，滿足觀測(cè)需求。

1.2數(shù)據(jù)處理與反演

本研究采用遠(yuǎn)參考和互參考技術(shù)，分別使用SSMT-2000和MT-Editor軟件進(jìn)行數(shù)據(jù)處理與編輯。原始時(shí)間序列經(jīng)過(guò)快速傅里葉變換由時(shí)間域轉(zhuǎn)化為頻率域;在此基礎(chǔ)上，通過(guò)遠(yuǎn)參考Robust估計(jì)等處理技術(shù)［12-13］獲得與頻率有關(guān)的阻抗張量信息。

1.2.1 視電阻率與相位曲線

典型測(cè)點(diǎn)的視電阻率與相位曲線（圖3）顯示，不同構(gòu)造單元的曲線類型存在明顯差異。黃驊坳陷（圖3中116～169號(hào)）顯示的曲線類型為KH型，而薊寶褶皺帶（圖3中209～244）的曲線類型為K型。曲線類型的不同直接反映出不同構(gòu)造單元電性結(jié)構(gòu)的巨大差異。

1.2.2 構(gòu)造走向與維性特征

通過(guò)Groom-Bailey（簡(jiǎn)稱GB）分解［14-15］進(jìn)行構(gòu)造走向分析，繪制了7個(gè)頻段所有測(cè)點(diǎn)電性主軸方位的玫瑰圖（圖4）。由圖可知，中頻部分與總的電性統(tǒng)計(jì)類似，高頻和低頻部分則略有差異。盡管每個(gè)頻段電性主軸方位的統(tǒng)計(jì)結(jié)果不盡相同，但綜合來(lái)看電性主軸方位為10°或100°。

在一定條件下，三維地質(zhì)情況可以近似為二維或一維進(jìn)行處理。本研究中，剖面布置大致垂直于構(gòu)造走向［16］。采用Swift方法計(jì)算二維偏離度（圖5），可以判別大地電磁測(cè)深數(shù)據(jù)是否滿足二維性要求。通常認(rèn)為二維偏離度的值小于0.3時(shí)，可以近似為二維。由圖5可見，盡管在各構(gòu)造單元邊界，大地電磁測(cè)深數(shù)據(jù)在低頻部分表現(xiàn)出一定的三維性，但剖面整體表現(xiàn)基本滿足二維性假設(shè)，適于進(jìn)行二維反演。

1.2.3 數(shù)據(jù)反演

本研究采用二維非線性共軛梯度法（Nonlinear conjugate gradient，簡(jiǎn)稱NLCG）進(jìn)行反演計(jì)算，根據(jù)TE+TM聯(lián)合模式設(shè)置不同的參數(shù)進(jìn)行反演。視電阻率誤差級(jí)數(shù)10%，相位誤差級(jí)數(shù)5%。反演過(guò)程中使用了0.001～10 000 s共七個(gè)數(shù)量級(jí)的MT數(shù)據(jù)，初始模型為100 Ω·m均勻半空間。

正則化因子τ分別取1，3，5，7，10，30，100時(shí)，均方根誤差（root mean square error，簡(jiǎn)稱RMS）與模型光滑度對(duì)應(yīng)關(guān)系曲線如圖6。當(dāng)τ=1，3，5，7，10時(shí)，RMS隨著模型光滑度的增加緩慢增大;而當(dāng)τ=30，100時(shí)，RMS突然陡增，故反演選擇正則化因子τ=10。停止反演時(shí)的RMS擬合差為2.79。

2 電性結(jié)構(gòu)特征分析

2.1 剖面總體電性特征

經(jīng)二維反演計(jì)算獲得的電性結(jié)構(gòu)剖面顯示，研究區(qū)在NNE方向上表現(xiàn)出明顯的縱向分層和橫向分塊特征（圖7）。

縱向上，從上至下大致分為淺表低阻層、中間高阻層和深部低阻層。這一分層特征與地殼不同深度的物質(zhì)組成、流體含量和溫壓條件等強(qiáng)烈相關(guān)。淺表低阻層的平均厚度約為3～4 km，該層內(nèi)主要分布的是第四紀(jì)松散沉積物。由于松散沉積物孔隙度高、流體含量豐富，該層整體表現(xiàn)出明顯的低阻特征。在淺表松散沉積層之下的基巖主要是中生代火山沉積巖系和前中生代碳酸鹽巖，在更深處則是厚度巨大的基底變質(zhì)巖系。相對(duì)于淺表的松散沉積物，中深層內(nèi)的巖石具有更低的孔隙度和流體含量，因此總體表現(xiàn)為高電阻特征。隨著深度的增加，溫壓條件也發(fā)生顯著變化，尤其是溫度的升高誘使巖石發(fā)生化學(xué)反應(yīng)，使得造巖礦物中的化學(xué)束縛水析出，從而降低了巖石的電阻率。因此，在電性結(jié)構(gòu)剖面的深部總體表現(xiàn)為低阻特征。

橫向上，中間高阻層的分塊特征明顯。從區(qū)域尺度看，大致以薊運(yùn)河斷裂（F5）為界，電性結(jié)構(gòu)剖面可分為西南段的低阻區(qū)和東北段的高阻區(qū)。西南段華北斷陷盆地區(qū)域的電阻率總體較低，僅在漢沽斷裂和薊運(yùn)河斷裂之間出現(xiàn)相對(duì)的高阻體;而東北段燕山隆起區(qū)的電阻率則相對(duì)較高。這一顯著差異正是兩個(gè)不同構(gòu)造單元間巨大理化條件差異的體現(xiàn)。前人研究顯示，華北斷陷盆地區(qū)的平均熱流值明顯高于燕山隆起區(qū)［17-18］，這意味著高熱流值地區(qū)巖石可能具有個(gè)更高的溫度，從而引起流體的增加和電阻率的降低。這種地下熱環(huán)境的不同可能是不同構(gòu)造單元間電性結(jié)構(gòu)差異的重要因素之一。此外，各構(gòu)造單元內(nèi)部也同樣存在明顯的分塊特征，主要表現(xiàn)為縱向高阻體和低阻體的相間排列，且在大型低阻體區(qū)域幾乎都有斷裂的發(fā)育，可能暗示了構(gòu)造單元內(nèi)部的分塊特征與區(qū)域性斷裂存在一定的關(guān)聯(lián)。事實(shí)上，這一特征并不難理解。斷裂破碎帶的存在，直接為地下流體提供運(yùn)移通道，從而降低了該區(qū)域的電阻率。

莫霍面及軟流圈頂界深度引自文獻(xiàn)［19-20］，結(jié)合本文獲取剖面的電性結(jié)構(gòu)特征綜合推測(cè)給出其大致位置;圖中圓圈表示地震，為本次測(cè)線兩側(cè)50 km范圍的地震震中投影。圓圈大小與震級(jí)為正比，圖中最大圓圈為寧河M6.9地震震中投影，兩個(gè)次大圓圈為漢沽附近兩次M6.2地震震中投影。

2.2 各構(gòu)造單元的電性特征

2.2.1 華北斷陷區(qū)

剖面西南段穿過(guò)華北斷陷區(qū)的黃驊坳陷和滄縣隆起（圖2，圖7），各次級(jí)斷陷盆地中的新生界厚度多大于2 km，部分地區(qū)可達(dá)6～7 km［9］。電性結(jié)構(gòu)剖面（圖7）顯示，除了與燕山隆起相鄰近的寧河凸起一帶外，華北斷陷區(qū)主要表現(xiàn)為100 Ω·m以下低電阻率。

在以第四紀(jì)松散沉積為主的淺部區(qū)域，低阻層厚度由西南向東北從3～4 km逐漸減薄至1 km以下。這一特征在深地震反射剖面（圖8）上也有很好的揭示［21］，與本區(qū)的次級(jí)構(gòu)造單元地質(zhì)情況相吻合。在中深層區(qū)域，華北斷陷區(qū)主要呈現(xiàn)出低阻主導(dǎo)、高阻間隔出現(xiàn)的縱向排布特征。此外，電阻率由西南向東北方向總體呈增大趨勢(shì)。

北塘凹陷下方發(fā)育一個(gè)巨大的“Y”型低阻體，電阻率小于100 Ω·m。在地殼中，該低阻體的西南支與北塘凹陷沉積盆底中心相連，東北支與漢沽斷裂形成有效連通（圖7）。兩支低阻體大致在25 km深處合而為一，并一直向下延伸至100 km處。寧河凸起下方的高阻體內(nèi)，上部電阻率可達(dá)10 000 Ω·m，中下部電阻率逐漸下降至100 Ω·m左右。苗莊凸起下方的低阻體向下延伸到70 km深處，并沿軟流圈頂界附近向北轉(zhuǎn)折，延伸至燕山隆起下方（圖7）。

2.2.2 燕山隆起區(qū)

剖面東北段主要穿過(guò)燕山隆起區(qū)的薊寶褶皺帶（圖2，圖7），該區(qū)內(nèi)的第四系厚度較小，并向北部山區(qū)進(jìn)一步減薄。電性結(jié)構(gòu)剖面顯示，薊寶褶皺帶以高電阻率為特征。在淺部第四系覆蓋區(qū)，低阻層相比華北斷陷區(qū)更薄。在5～40 km的中深層區(qū)域，則表現(xiàn)為顯著的高阻特征與橫向分塊特征（圖7）。該處發(fā)育有3個(gè)電阻率約為100～10 000 Ωm次級(jí)高阻體，與介于其間的中高阻區(qū)縱向間隔排布，并向下延伸至60 km深處。

豐臺(tái)—野雞坨斷裂（F6）下方為一規(guī)模較小的局部低阻體，與下方80 km左右的低阻體相對(duì)應(yīng)。該處90 km以下又出現(xiàn)了明顯的橫向分布的低阻體，與北塘凹陷及薊運(yùn)河斷裂帶下方的低阻體似乎有貫通趨勢(shì)，形成了一個(gè)變形的“Y”型低阻體，交叉點(diǎn)在100 km左右，并一直延伸到120 km以下。

2.3 主要斷裂的電性特征

區(qū)域性斷裂往往構(gòu)成不同地質(zhì)體之間的邊界，其兩側(cè)的地質(zhì)體通常在物質(zhì)組成、結(jié)構(gòu)、構(gòu)造、流體含量等方面存在較大差異，因此它們?cè)陔娦越Y(jié)構(gòu)上的表現(xiàn)也不盡相同。在斷裂帶淺部，巖石受應(yīng)力作用發(fā)生破碎，引起孔隙度和流體增加，則可能導(dǎo)致斷裂帶電阻率降低;而在斷裂帶深部，巖石受韌性剪切作用發(fā)生礦物定向排列，亦可能改變其各向異性特征［22］。因此，結(jié)合已有地震勘探、深鉆孔及地質(zhì)勘查等資料開展電性結(jié)構(gòu)特征分析，可以在剖面上識(shí)別出主要斷裂帶的位置及其構(gòu)造特征。本次獲取的電性結(jié)構(gòu)剖面顯示，斷裂主要分布在華北斷陷區(qū)內(nèi)的高低阻梯度帶附近。

海河斷裂（F1）：海河斷裂是一條鏟式正斷層，向下延伸至8 km左右［23］。在電性剖面上，斷裂頂端位于測(cè)點(diǎn)124～126之間，切割深度較淺，是一條典型的殼內(nèi)斷裂。斷裂向深部延伸與高低阻梯度帶逐漸重合。受斷層影響，斷裂西南側(cè)5 km深處的電阻率等值線發(fā)生了明顯的扭曲（圖7），其下凹形態(tài)與石油資料揭示的板橋凹陷形態(tài)吻合。

茶淀斷裂（F2）：茶淀斷裂是一條傾向S的正斷層，斷裂頂端位于測(cè)點(diǎn)138～140之間，下部位于高低阻梯度帶上，并延伸至10 km深度附近（圖7）。斷裂南北兩側(cè)分別為電阻率小于20 Ω·m的低阻體和電阻率大于100 Ω·m的高阻體。此外，南側(cè)低阻體與下方的顯著低阻體存在貫通趨勢(shì)。

漢沽斷裂（F3）：測(cè)點(diǎn)158～160之間的漢沽斷裂是北塘凹陷和寧河凸起的分界斷裂，正好位于其下方低阻體邊緣的高低阻梯度帶上。剖面顯示，該斷裂可能是一條切穿莫霍面的深斷裂（圖7）。此外，淺層人工地震勘探結(jié)果表明，漢沽斷裂的淺部結(jié)構(gòu)單一，斷層面傾角較大［4］，與本次獲取的深部大傾角斷層相協(xié)調(diào)。

薊運(yùn)河斷裂（F5）：薊運(yùn)河主斷裂位于測(cè)點(diǎn)181～183之間，是黃驊坳陷和薊寶褶皺帶的分界線。斷裂受到上地殼高阻體阻擋未明顯下切，其下方發(fā)育有延伸至上地幔的低阻體。寧河M6.9地震震源區(qū)正好位于該低阻體南側(cè)的高低阻梯度帶上（圖7）?？赡苁苌系貧じ咦梵w影響，震中能量向上傳播受阻，此次地震并未造成明顯的地表破壞。淺層人工地震剖面［24］顯示，薊運(yùn)河斷裂的淺部結(jié)構(gòu)表現(xiàn)為典型的正花狀構(gòu)造，可能是深部斷裂走滑活動(dòng)的表現(xiàn)，其潛在地震最大震級(jí)也可能因此而提升，地震風(fēng)險(xiǎn)值得關(guān)注。

豐臺(tái)—野雞坨斷裂（F6）：測(cè)點(diǎn)205～206之間的豐臺(tái)—野雞坨斷裂位于高阻體內(nèi)的相對(duì)低阻帶上（圖7）。斷裂控制了前新生界基巖的埋深，其兩側(cè)第四系厚度之差可達(dá)300～400 m，并在豐潤(rùn)西南發(fā)現(xiàn)有第四系厚度高梯度帶［25］。在深地震反射剖面［21］上（圖8），該斷裂在淺部錯(cuò)斷地層明顯，而向深處的延伸并不清晰。

3 孕震環(huán)境區(qū)域動(dòng)力學(xué)背景討論

3.1 孕震區(qū)域背景

強(qiáng)震的發(fā)生往往與深部構(gòu)造背景有關(guān)。探測(cè)剖面兩側(cè)50 km范圍內(nèi)的2級(jí)以上地震震中投影顯示，剖面附近的三個(gè)6級(jí)以上地震均發(fā)生在區(qū)域性斷裂附近，暗示了區(qū)內(nèi)中強(qiáng)地震的發(fā)生與斷裂活動(dòng)相關(guān)。從電性結(jié)構(gòu)上看，這三個(gè)地震的震中均位于大型低阻體邊緣電性結(jié)構(gòu)變化較劇烈的梯度帶附近，與翁愛華等［26］給出的結(jié)論基本一致。此外，另一些震級(jí)相對(duì)較大的地震也主要發(fā)生在大型低阻體或高阻體邊緣的高低阻梯度帶附近，僅有極少數(shù)位于大型低阻體或高阻體內(nèi)部。

3.1.1 寧河M6.9地震孕震環(huán)境

寧河M6.9地震發(fā)生在華北斷陷盆地的東北部，震中位于薊運(yùn)河斷裂正下方［27-28］或南側(cè)［29］。其震源機(jī)制解的兩個(gè)節(jié)面走向分別為NNW向（329°）和NEE向（59°）。張之立等［28］認(rèn)為，NNW向節(jié)面是可能的斷層面，且斷層以走滑為主，與薊運(yùn)河斷裂的走向、斷層性質(zhì)一致。李志義等［30］認(rèn)為，唐山斷裂帶向西南延伸受到薊運(yùn)河斷裂阻擋，產(chǎn)生新的應(yīng)力集中區(qū)并發(fā)生破裂，從而引發(fā)寧河M6.9地震。

本研究給出的大地電磁測(cè)深剖面顯示，寧河M6.9地震與大型低阻體邊緣的高低阻梯度帶密切相關(guān)（圖7）。在薊運(yùn)河斷裂下方存在一個(gè)源自軟流圈的低阻體，該低阻體在燕山隆起區(qū)下方近水平展布，在大約70 km深處轉(zhuǎn)為近直立，并一直向上延伸至5 km深處。寧河M6.9地震即發(fā)生在低阻體頂端與其南側(cè)高阻體交匯部位的高低阻梯度帶上，并與薊運(yùn)河斷裂的空間位置十分接近。

此外，在薊運(yùn)河斷裂下方低阻體與第四系松散沉積之間存在厚約3～5 km的高阻體（圖7）。該高阻體的存在很可能對(duì)下方低阻體的向上發(fā)展和薊運(yùn)河斷裂的下切產(chǎn)生了顯著的影響。結(jié)合寧河M6.9地震的宏觀烈度偏低、沒有明顯的地震斷層出露、中深層人工地震資料揭示下切不明顯［21］以及鉆孔地層對(duì)比揭示上斷點(diǎn)未達(dá)到全新世［24］等情況看，薊運(yùn)河斷裂的深、淺部破裂可能尚未形成有效貫通。因此，薊運(yùn)河斷裂是否為寧河M6.9地震的發(fā)震構(gòu)造還需更多探測(cè)結(jié)果支持，但該斷裂潛在的中強(qiáng)地震風(fēng)險(xiǎn)仍然需要持續(xù)關(guān)注。

3.1.2 漢沽M6.2地震孕震環(huán)境

在唐山地震期間，漢沽斷裂南側(cè)發(fā)生了兩次M6.2地震，大地電磁測(cè)深剖面清晰顯示了這兩次地震的孕震環(huán)境。漢沽斷裂下方發(fā)育的巨大低阻體，自100 km深處一直向上延伸至地表附近，并在莫霍面附近分為南北兩支，分別與茶淀斷裂和漢沽斷裂貫通（圖7）。兩次M6.2地震均發(fā)生在電阻率梯度更大的北側(cè)分支附近。

研究顯示，唐山斷裂帶基本止于漢沽附近［31］，且唐山余震基本止于漢沽斷裂。這表明漢沽斷裂可能是唐山斷裂帶向南活動(dòng)阻擋斷裂。本次研究顯示，漢沽斷裂割深度可能達(dá)到了莫霍面附近（圖7），極有可能是兩次M6.2地震的發(fā)震構(gòu)造。

3.2 區(qū)域孕震動(dòng)力環(huán)境討論

研究區(qū)位于渤海灣盆地西岸，處在華北克拉通巖石圈減薄的中心地帶，大規(guī)模的構(gòu)造變形和大地震發(fā)生是其重要表現(xiàn)。多年的研究認(rèn)為，華北克拉通破壞是西太平洋俯沖帶遠(yuǎn)程效應(yīng)的結(jié)果。古太平洋板塊的西向俯沖，導(dǎo)致了軟流圈地幔的上涌，形成多期次的巖漿活動(dòng)并導(dǎo)致了巖石圈底部和上地殼的同時(shí)減薄［32-33］。黃驊坳陷區(qū)巖石圈厚度僅為70～80 km，而燕山褶皺帶的巖石圈厚度達(dá)到了110 km，進(jìn)一步顯示巖石圈減薄在中心地帶是十分顯著的。

本次研究發(fā)現(xiàn)了兩個(gè)十分值得關(guān)注的大型低阻體，這兩個(gè)低阻體均源自軟流圈及上地幔，幾乎貫通了整個(gè)地殼。研究顯示，低阻體與部分熔融或熱流體密切相關(guān)［34-35］。這些大型低阻體的發(fā)育暗示著較大規(guī)模部分熔融或熱流體的存在，可能與華北克拉通破壞過(guò)程中的上地幔熱物質(zhì)上涌［36］有關(guān)。同時(shí)，低阻體內(nèi)較大規(guī)模斷裂帶的存在，也為這些流體的運(yùn)移提供了重要通道，為低阻體的向上發(fā)展創(chuàng)造了便利條件。

渤海灣盆地巖石圈結(jié)構(gòu)是典型的“冷殼熱?！苯Y(jié)構(gòu)，由于上地幔熱物質(zhì)上涌，導(dǎo)致上部冷的地殼拉伸變形甚至破裂。同時(shí)，隨著底侵作用與上地幔減壓作用的發(fā)生，上地幔熱物質(zhì)通過(guò)破裂或深大斷裂帶等薄弱部位發(fā)生侵入并誘使冷的地殼發(fā)生部分熔融，在某些部位與淺部先存或新生斷層形成了有效貫通，或形成塌陷盆地，進(jìn)而形成對(duì)淺部地質(zhì)構(gòu)造的強(qiáng)烈控制［36］。這一動(dòng)力演化過(guò)程引起區(qū)域構(gòu)造應(yīng)力狀態(tài)的不斷調(diào)整［37］，即表現(xiàn)為大量地震的發(fā)生。漢沽兩次M6.2地震和寧河M6.9地震，乃至2003—2004年塘沽—北塘小震群的發(fā)生均是這種區(qū)域應(yīng)力調(diào)整的結(jié)果。因此，漢沽斷裂與薊運(yùn)河斷裂的地震活動(dòng)性需進(jìn)一步關(guān)注。

4 結(jié)論

本文利用跨越張渤地震帶天津—塘沽分段的大地電磁測(cè)深剖面，獲取了精細(xì)殼幔電性結(jié)構(gòu)。通過(guò)詳細(xì)分析主要構(gòu)造單元及主要斷裂的電性結(jié)構(gòu)特征，得到以下認(rèn)識(shí)：

（1） 華北平原東北部電性結(jié)構(gòu)呈現(xiàn)出顯著的“橫向分塊、縱向分層”特征，揭示了不同構(gòu)造單元在物質(zhì)組成和地質(zhì)構(gòu)造等方面的巨大差異。

（2） 漢沽斷裂可能是一條下切莫霍面的深大斷裂，薊運(yùn)河斷裂的發(fā)育可能與深部的走滑構(gòu)造活動(dòng)密切相關(guān)，兩條斷裂的地震風(fēng)險(xiǎn)值得關(guān)注。

（3） 寧河M6.9地震和漢沽兩次M6.2地震的發(fā)生與大型低阻體邊緣的高低阻梯度帶密切相關(guān)。

致謝：評(píng)審專家對(duì)本文的修改提出了大量寶貴意見和建議，提升了稿件質(zhì)量，在此表示衷心感謝！

參考文獻(xiàn)（References）

［1］ 徐杰，宋長(zhǎng)青，楚全芝.張家口—蓬萊斷裂帶地震構(gòu)造特征的初步探討［J］.地震地質(zhì)，1998，20（2）：146-154.

XU Jie，SONG Changqing，CHU Quanzhi.Preliminary study on the seismotectonic characters of the Zhangjiakou—Penglai fault zone［J］.Seismology and Geology，1998，20（2）：146-154.

［2］ 方穎，張晶.張家口—渤海斷裂帶分段活動(dòng)性研究［J］.地震，2009，29（3）：136-140.

FANG Ying，ZHANG Jing.Study of the segmentary activities in the Zhangjiakou—Bohai fault system［J］.Earthquake，2009，29（3）：136-140.

［3］ 國(guó)家地震局.一九七六年唐山地震［M］.北京：地震出版社，1982：113-121.

State Seismological Bureau.Tangshan earthquake in 1976［M］.Beijing：Seismological Press，1982：113-121.

［4］ 閆成國(guó)，王志勝，陳宇坤，等.天津漢沽斷裂準(zhǔn)確位置及活動(dòng)性的綜合探測(cè)研究［J］.中國(guó)地震，2014，30（4）：501-513.

YAN Chengguo，WANG Zhisheng，CHEN Yukun，et al.Comprehensive survey and investigation on the exact spatial location and activity of the Hangu fault in Tianjin［J］.Earthquake Research in China，2014，30（4）：501-513.

［5］ 國(guó)家地震局震害防御司.中國(guó)歷史強(qiáng)震目錄（公元前23世紀(jì)—公元1911年）［M］.北京：地震出版社，1995.

The Earthquake Disaster Prevention Department of the State Seismological Bureau.Catalogue of Chinese history strong earthquakes （23rd century B.C.to 1911 A.D.）［M］.Beijing：Seismological Press，1995.

［6］ 劉昌銓，嘉世旭.唐山地震區(qū)地殼上地幔結(jié)構(gòu)特征：二維非均勻介質(zhì)中理論地震圖計(jì)算和結(jié)果分析［J］.地震學(xué)報(bào)，1986，8（4）：341-353.

LIU Changquan，JIA Shixu.Structural property of the crust and the upper mantle in the Tangshan earthquake region：a method of computing body wave theoretical seismogram and the result of analysis of 2D inhomogeneous media［J］.Acta Seismologica Sinica，1986，8（4）：341-353.

［7］ 徐志萍，方盛明，李德慶，等.利用布格重力資料研究華北裂陷盆地地殼結(jié)構(gòu)特征［J］.大地測(cè)量與地球動(dòng)力學(xué)，2017，37（3）：246-250.

XU Zhiping，F(xiàn)ANG Shengming，LI Deqing，et al.Characteristics of the crustal structure in North China rift-depression basin using bouguer gravity data［J］.Journal of Geodesy and Geodynamics，2017，37（3）：246-250.

［8］ 曲中黨，張訓(xùn)華，賀日政，等.華北盆地邊緣及鄰區(qū)地殼S波速度結(jié)構(gòu)及其地震孕育機(jī)制［J］.地球物理學(xué)進(jìn)展，2018，33（3）：957-968.

QU Zhongdang，ZHANG Xunhua，HE Rizheng，et al.S wave velocity structure of the crust and the mechanisms of earthquake occurrence in the North China Basin and its adjacent areas［J］.Progress in Geophysics，2018，33（3）：957-968.

［9］ 天津市地質(zhì)調(diào)查研究院.中國(guó)區(qū)域地質(zhì)志（天津志）［M］.北京：地質(zhì)出版社，2018.

Tianjin Institute of Geological Survey.The regional geology of China：Tianjin［M］.Beijing：Geological Publishing House，2018.

［10］ 謝成良，葉高峰，魏文博，等.藏北高原主要斷裂帶電性結(jié)構(gòu)特征［J］.地球物理學(xué)報(bào)，2012，55（12）：3991-4002.

XIE Chengliang，YE Gaofeng，WEI Wenbo，et al.Electrical features of the main faults beneath Northern Tibetan Plateau［J］.Chinese Journal of Geophysics，2012，55（12）：3991-4002.

［11］ 詹艷，趙國(guó)澤，王立鳳，等.河北石家莊地區(qū)深部結(jié)構(gòu)大地電磁探測(cè)［J］.地震地質(zhì)，2011，33（4）：913-927.

ZHAN Yan，ZHAO Guoze，WANG Lifeng，et al.Deep structure in Shijiazhuang and the vicinity by magnetotellurics［J］.Seismology and Geology，2011，33（4）：913-927.

［12］ 張泉，李永飛，田靖，等.復(fù)雜條件下大地電磁資料二維反演及初始模型構(gòu)建研究［J］.地球物理學(xué)進(jìn)展，2018，33（4）：1645-1651.

ZHANG Quan，LI Yongfei，TIAN Jing，et al.Two-dimensional inversion and initial model construction of magnetotelluric data in complex conditions are studied［J］.Progress in Geophysics，2018，33（4）：1645-1651.

［13］ 鄧琰，湯吉.大地電磁測(cè)深方法數(shù)據(jù)處理進(jìn)展［J］.地球物理學(xué)進(jìn)展，2019，34（4）：1411-1422.

DENG Yan，TANG Ji.Advances in magnetotelluric data processing［J］.Progress in Geophysics，2019，34（4）：1411-1422.

［14］ GROOM R W，BAILEY R C.Decomposition of magnetotelluric impedance tensors in the presence of local three-dimensional galvanic distortion［J］.Journal of Geophysical Research，1989，94（B2）：1913.

［15］ GROOM R W，BAILEY R C.Analytic investigations of the effects of near-surface three-dimensional galvanic scatterers on MT tensor decompositions［J］.Geophysics，1991，56（4）：496-518.

［16］ 張樂(lè)天，金勝，魏文博，等.青藏高原東緣及四川盆地的殼幔導(dǎo)電性結(jié)構(gòu)研究［J］.地球物理學(xué)報(bào)，2012，55（12）：4126-4137.

ZHANG Letian，JIN Sheng，WEI Wenbo，et al.Electrical structure of crust and upper mantle beneath the eastern margin of the Tibetan Plateau and the Sichuan Basin［J］.Chinese Journal of Geophysics，2012，55（12）：4126-4137.

［17］ 汪洋，鄧晉福，汪集旸，等.中國(guó)大陸熱流分布特征及熱-構(gòu)造分區(qū)［J］.中國(guó)科學(xué)院研究生院學(xué)報(bào)，2001，18（1）：51-58.

WANG Yang，DENG Jinfu，WANG Jiyang，et al.Terrestrial heat flow pattern and thermo tectonic domains in the continental area of China［J］.Journal of the Graduate School Academia Sinica，2001，18（1）：51-58.

［18］ NABELEK J，CHEN W P.The Tangshan earthquake sequence and its implications for the evolution of the North China Basin［J］.J Geophys Res，1987，92（B12）：12615-12628.

［19］ CHEN L，TAO W，ZHAO L，et al.Distinct lateral variation of lithospheric thickness in the northeastern North China Craton［J］.Earth and Planetary Science Letters，2008，267（1-2）：56-68.

［20］ 段永紅，王夫運(yùn)，張先康，等.華北克拉通中東部地殼三維速度結(jié)構(gòu)模型（HBCrust1.0）［J］.中國(guó)科學(xué)：地球科學(xué)，2016，46（6）：845-856.

DUAN Yonghong，WANG Fuyun，ZHANG Xiankang，et al.Three-dimensional velocity structure model of the crust in the central and eastern North China Craton （HBCrust1.0）［J］.Scientia Sinica （Terrae），2016，46（6）：845-856.

［21］ 閆成國(guó)，曹井泉，陳宇坤，等.深地震反射剖面揭示的天津地區(qū)張渤帶地殼精細(xì)結(jié)構(gòu)［J］.地球物理學(xué)報(bào)，2020，63（12）：4431-4439.

YAN Chengguo，CAO Jingquan，CHEN Yukun，et al.Fine crustal structures of Zhangjiakou—Bohai tectonic zone in Tianjin area revealed by a deep seismic reflection profile［J］.Chinese Journal of Geophysics，2020，63（12）：4431-4439.

［22］ WECKMANN U，RITTER O，HAAK V.A magnetotelluric study of the Damara belt in Namibia［J］.Physics of the Earth and Planetary Interiors，2003，138（2）：91-112.

［23］ 趙國(guó)敏，趙根模，陳宇坤.天津海河斷裂活動(dòng)性研究［M］.北京：地震出版社，2006.

ZHAO Guomin，ZHAO Genmo，CHEN Yukun.Study on the activity of Haihe fault in Tianjin［M］.Beijing：Seismological Press，2006.

［24］ 張文朋，張春麗，高武平，等.用淺層地震勘探資料研究薊運(yùn)河斷裂的第四紀(jì)活動(dòng)特征［J］.地震工程學(xué)報(bào)，2022，44（1）：183-191.

ZHANG Wenpeng，ZHANG Chunli，GAO Wuping，et al.Quaternary activity characteristics of the Jiyunhe fault revealed by shallow seismic prospecting data［J］.China Earthquake Engineering Journal，2022，44（1）：183-191.

［25］ 虢順民，李志義，程紹平，等.唐山地震區(qū)域構(gòu)造背景和發(fā)震模式的討論［J］.地質(zhì)科學(xué)，1977，12（4）：305-321.

GUO Shunmin，LI Zhiyi，CHENG Shaoping，et al.Discussion：on the recional structural background and the seismogenic model of the Tangshan earthquake［J］.Chinese Journal of Geology，1977，12（4）：305-321.

［26］ 翁愛華，李建平，范小平，等.大地電磁測(cè)深揭示的1668年郯城8.5級(jí)地震震中地殼精細(xì)結(jié)構(gòu)［J］.地震地質(zhì)，2018，40（2）：396-409.

WENG Aihua，LI Jianping，F(xiàn)AN Xiaoping，et al.Fine electrical structure beneath the epicenter of 1668 Tancheng MS8.5 earthquake revealed by MT sounding［J］.Seismology and Geology，2018，40（2）：396-409.

［27］ NBLEK J，CHEN W P，YE H.The Tangshan earthquake sequence and its implications for the evolution of the North China Basin［J］.Journal of Geophysical Research，1987，92（B12）：12615.

［28］ 張之立，李欽祖，谷繼成，等.唐山地震的破裂過(guò)程及其力學(xué)分析［J］.地震學(xué)報(bào)，1980，2（2）：111-129.

ZHANG Zhili，LI Qinzu，GU Jicheng，et al.The fracture processes of the Tangshan earthquake and its mechanical analysis［J］.Acta Seismologica Sinica，1980，2（2）：111-129.

［29］ 中國(guó)地震局震害防御司.中國(guó)近代地震目錄（公元1912年—1990年，MS≥4.7）［M］.北京：中國(guó)科學(xué)技術(shù)出版社，1999.

Department of Earthquake Disaster Prevention of the China Earthquake Administration.Catalogue of Chinese earthquakes （1912 A.D. to 1990 A.D.，MS≥4.7）［M］.Beijing：Science and Technology of China Press，1999.

［30］ 李志義，虢順民.從地震構(gòu)造看寧河6.9級(jí)余震與唐山主震的關(guān)系［J］.地震地質(zhì)，1979，1（4）：27-35.

LI Zhiyi，GUO Shunmin.On the relationship between the Ninghe 6.9 aftershock and the Tangshan earthquake from the viewpoint of seismotectonics［J］.Seismology and Geology，1979，1（4）：27-35.

［31］ 郭慧，江娃利，謝新生.對(duì)1976年河北唐山MS7.8地震地表破裂帶展布及位移特征的新認(rèn)識(shí)［J］.地震地質(zhì)，2011，33（3）：506-524.

GUO Hui，JIANG Wali，XIE Xinsheng.New evidence for the distribution of surface rupture zone of the 1976 MS7.8 Tangshan earthquake［J］.Seismology and Geology，2011，33（3）：506-524.

［32］ LIU J L，SHEN L，JI M，et al.The Liaonan/Wanfu metamorphic core complexes in the Liaodong Peninsula：two stages of exhumation and constraints on the destruction of the North China Craton［J］.Tectonics，2013，32（5）：1121-1141.

［33］ ZHU R X，XU Y G，ZHU G，et al.Destruction of the North China Craton［J］.Science China Earth Sciences，2012，55（10）：1565-1587.

［34］ 李立，金國(guó)元.攀西裂谷帶及龍門山斷裂帶地殼上地幔的大地電磁測(cè)深研究［J］.物探與化探，1987，11（3）：161-169.

LI Li，JIN Guoyuan.Telluric electromagnetic sounding study of crust and upper mantle in the Panxi “rift zone” and the Longmenshan faulted zone［J］.Geophysical and Geochemical Exploration，1987，11（3）：161-169.

［35］ 魏文博，金勝，葉高峰，等.藏南巖石圈導(dǎo)電性結(jié)構(gòu)與流變性：超寬頻帶大地電磁測(cè)深研究結(jié)果［J］.中國(guó)科學(xué)（D輯：地球科學(xué)），2009，39（11）：1591-1606.

WEI Wenbo，JIN Sheng，YE Gaofeng，et al.Conductive structure and rheology of lithosphere in southern Tibet：research results of ultra-wideband magnetotelluric sounding［J］.Science in China （Series D：Earth Sciences），2009，39（11）：1591-1606.

［36］ ZHU R X，CHEN L，WU F Y，et al.Timing，scale and mechanism of the destruction of the North China Craton［J］.Science China Earth Sciences，2011，54（6）：789-797.

［37］ 漆家福.渤海灣新生代盆地的兩種構(gòu)造系統(tǒng)及其成因解釋［J］.中國(guó)地質(zhì)，2004，31（1）：15-22.

QI Jiafu.Two tectonic systems in the Cenozoic Bohai Bay Basin and their genetic interpretation［J］.Geology in China，2004，31（1）：15-22.