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

高武平 林逸 張文朋 徐新學 閆成國 曹井泉 劉紅艷

摘要：

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

關鍵詞：

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

中圖分類號： P315????? 文獻標志碼：A?? 文章編號： 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 引言

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

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

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

1 大地電磁測深剖面

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

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

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

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

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

1.2.1 視電阻率與相位曲線

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

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

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

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

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

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

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

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

2.1 剖面總體電性特征

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

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

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

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

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

2.2.1 華北斷陷區(qū)

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

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

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

2.2.2 燕山隆起區(qū)

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

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

2.3 主要斷裂的電性特征

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

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

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

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

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

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

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

3.1 孕震區(qū)域背景

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

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

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

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

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

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

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

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

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

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

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

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

4 結(jié)論

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

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

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

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

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

參考文獻（References）

［1］ 徐杰，宋長青，楚全芝.張家口—蓬萊斷裂帶地震構(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］ 方穎，張晶.張家口—渤海斷裂帶分段活動性研究［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］ 國家地震局.一九七六年唐山地震［M］.北京：地震出版社，1982：113-121.

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

［4］ 閆成國，王志勝，陳宇坤，等.天津漢沽斷裂準確位置及活動性的綜合探測研究［J］.中國地震，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］ 國家地震局震害防御司.中國歷史強震目錄（公元前23世紀—公元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ì)中理論地震圖計算和結(jié)果分析［J］.地震學報，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］.大地測量與地球動力學，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］ 曲中黨，張訓華，賀日政，等.華北盆地邊緣及鄰區(qū)地殼S波速度結(jié)構(gòu)及其地震孕育機制［J］.地球物理學進展，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)查研究院.中國區(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］.地球物理學報，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］ 詹艷，趙國澤，王立鳳，等.河北石家莊地區(qū)深部結(jié)構(gòu)大地電磁探測［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］ 張泉，李永飛，田靖，等.復雜條件下大地電磁資料二維反演及初始模型構(gòu)建研究［J］.地球物理學進展，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］ 鄧琰，湯吉.大地電磁測深方法數(shù)據(jù)處理進展［J］.地球物理學進展，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］ 張樂天，金勝，魏文博，等.青藏高原東緣及四川盆地的殼幔導電性結(jié)構(gòu)研究［J］.地球物理學報，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］ 汪洋，鄧晉福，汪集旸，等.中國大陸熱流分布特征及熱-構(gòu)造分區(qū)［J］.中國科學院研究生院學報，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］ 段永紅，王夫運，張先康，等.華北克拉通中東部地殼三維速度結(jié)構(gòu)模型（HBCrust1.0）［J］.中國科學：地球科學，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］ 閆成國，曹井泉，陳宇坤，等.深地震反射剖面揭示的天津地區(qū)張渤帶地殼精細結(jié)構(gòu)［J］.地球物理學報，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］ 趙國敏，趙根模，陳宇坤.天津海河斷裂活動性研究［M］.北京：地震出版社，2006.

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

［24］ 張文朋，張春麗，高武平，等.用淺層地震勘探資料研究薊運河斷裂的第四紀活動特征［J］.地震工程學報，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ì)科學，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］ 翁愛華，李建平，范小平，等.大地電磁測深揭示的1668年郯城8.5級地震震中地殼精細結(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］ 張之立，李欽祖，谷繼成，等.唐山地震的破裂過程及其力學分析［J］.地震學報，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］ 中國地震局震害防御司.中國近代地震目錄（公元1912年—1990年，MS≥4.7）［M］.北京：中國科學技術(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］.地震地質(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］ 郭慧，江娃利，謝新生.對1976年河北唐山MS7.8地震地表破裂帶展布及位移特征的新認識［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］ 李立，金國元.攀西裂谷帶及龍門山斷裂帶地殼上地幔的大地電磁測深研究［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］ 魏文博，金勝，葉高峰，等.藏南巖石圈導電性結(jié)構(gòu)與流變性：超寬頻帶大地電磁測深研究結(jié)果［J］.中國科學（D輯：地球科學），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］.中國地質(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.