范曄, 陳小斌, 湯吉, 崔騰發(fā), 孫翔宇, 王培杰, 劉鐘尹
1 中國(guó)地震局地質(zhì)研究所地震動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100029 2 中國(guó)地震臺(tái)網(wǎng)中心, 北京 100045 3 應(yīng)急管理部國(guó)家自然災(zāi)害防治研究院, 北京 100085 4 中國(guó)地震局地震預(yù)測(cè)研究所, 北京 100036
郯廬斷裂帶是中國(guó)東部長(zhǎng)約2400 km的巨型斷裂帶,沿北北東向穿越了多個(gè)大地構(gòu)造單元,控制著中國(guó)東部的主要構(gòu)造格局.位于蘇魯造山帶西側(cè)及以南的區(qū)段被劃分為郯廬斷裂帶中南段(晁洪太等,1997; 施煒等,2003).在中生代中晚期,揚(yáng)子板塊和華北板塊發(fā)生碰撞,古太平洋板塊向西北方向俯沖(Zhu et al., 2008; Zhu and Zheng, 2009),郯廬斷裂帶發(fā)生大規(guī)模左旋走滑運(yùn)動(dòng),將蘇魯造山帶和大別造山帶左旋錯(cuò)開(kāi)400多公里(萬(wàn)天豐,1996;徐嘉煒和馬國(guó)鋒,1992;張鵬等,2007).之后郯廬斷裂帶中南段經(jīng)歷了中生代晚期的斷陷和新生代早期的側(cè)向擠壓等多期次演化(朱光等,2001;鄭朗蓀等,1988),在新生代晚期以來(lái)表現(xiàn)為右旋走滑兼逆沖的活動(dòng)特征(方仲景等,1980;張鵬等,2010).
按照斷裂的構(gòu)造習(xí)性和地震活動(dòng)性,一般以嘉山為界將郯廬斷裂帶中南段劃分為中段和南段(鄭朗蓀等,1988;徐嘉煒和馬國(guó)鋒,1992;張鵬等,2007).嘉山以北的中段是第四紀(jì)活動(dòng)斷裂集中分布的段落,由4~5條近似平行的主要斷裂組成;嘉山以南的南段由2條主要斷裂組成,新構(gòu)造變形較弱(晁洪太等,1997; 李臘月等,2020).南北兩側(cè)斷層性質(zhì)和形態(tài)的分段可能受控于斷層的空間展布和深部物性結(jié)構(gòu)特征.歷史地震表明,郯廬斷裂帶中段發(fā)生過(guò)郯城8.5級(jí)地震、渤海7級(jí)地震等強(qiáng)震,南段少有6級(jí)以上地震發(fā)生.數(shù)值模擬結(jié)果表明郯廬斷裂帶渤海以南的中南段可能處于高閉鎖狀態(tài)(李彥川等,2016),較低的b值也顯示郯廬斷裂帶南段積累了較高的應(yīng)力(朱艾斕等,2018).從小震精定位空間分布來(lái)看(朱艾斕等,2018),嘉山以南地震活動(dòng)性較弱.造成這些現(xiàn)象的原因究竟是郯廬斷裂帶南段處于斷層閉鎖造成的地震沉默期,還是此區(qū)域缺乏強(qiáng)震孕育的深部環(huán)境或外部驅(qū)動(dòng)力?由于目前缺少該區(qū)段有針對(duì)性的三維高分辨率深部資料,無(wú)法給出確切答案.此外,郯廬斷裂帶合肥—宿遷段是連接大別造山帶與蘇魯造山帶的樞紐,其深部結(jié)構(gòu)的研究對(duì)認(rèn)識(shí)郯廬斷裂帶的形成和分段差異性亦有重要意義.
目前已開(kāi)展的地震成像結(jié)果(劉保金等,2015;秦晶晶等,2018;孟亞峰等,2019;范小平等,2017;Bem et al., 2020; 顧勤平等,2020;Luo et al., 2021)給出了包含研究區(qū)的深部速度模型,并根據(jù)速度模型對(duì)斷裂帶進(jìn)行了分段分析,但有關(guān)斷裂帶展布范圍、延伸深度、分段特征等方面還有不同認(rèn)識(shí).重力和磁法勘察結(jié)果表明,郯廬斷裂帶在磁場(chǎng)上呈現(xiàn)出一條NNE走向線性磁異常帶,重力場(chǎng)表現(xiàn)為規(guī)模較大的布格重力異常梯度帶,并且斷裂帶兩側(cè)磁場(chǎng)和重力場(chǎng)在場(chǎng)值、異常展布形態(tài)等方面有明顯區(qū)別(霍光輝等,2009;王鑫等,2016),居里面深度在兩側(cè)也有明顯差異(邵世德和嚴(yán)曉紅,1992),表現(xiàn)出西深東淺的特征.這些研究都著重于大結(jié)構(gòu)(重磁和地震層析成像),或者單條剖面所途徑的區(qū)域(地震反射),對(duì)于分析郯廬帶合肥—宿遷段這個(gè)局部區(qū)段,所能提供的信息比較有限.
在郯廬斷裂帶中南段已經(jīng)開(kāi)展了一些大地電磁研究工作.張繼紅等(2010)依據(jù)郯廬帶諸城段的深部電性結(jié)構(gòu),認(rèn)為郯廬帶中段存在著與殼內(nèi)高導(dǎo)體分布相關(guān)的強(qiáng)震孕育環(huán)境.肖騎彬等(2008)通過(guò)臨沂段的剖面結(jié)果,發(fā)現(xiàn)蘇魯超高壓變質(zhì)巖以大高阻特征為標(biāo)志,且深入到地幔,郯廬斷裂帶在臨沂附近是一條殼內(nèi)斷裂.葉高峰等(2009)在泰安—日照剖面的結(jié)果中,發(fā)現(xiàn)沂沭斷裂帶淺部的兩塹夾一壘特征顯著.宋國(guó)奇(2006,2007)依據(jù)渤海和魯蘇皖段幾條深部電性結(jié)構(gòu)剖面發(fā)現(xiàn),沂沭斷裂帶盆地與深部的殼內(nèi)高導(dǎo)體存在對(duì)應(yīng)關(guān)系.胡惟等(2014)對(duì)比分析了強(qiáng)震組、弱震組和無(wú)震組的電性結(jié)構(gòu)特征,也對(duì)郯廬帶中段的孕震環(huán)境進(jìn)行了討論.翁愛(ài)華等(2018)布置了通過(guò)郯城地震震中的E-W向剖面,分析了主干斷裂深部延伸狀態(tài)和郯城地震震中位置.這幾項(xiàng)研究集中于研究區(qū)以北,目前在研究區(qū)內(nèi)部的大地電磁研究結(jié)果很少.陳小斌和郭春玲(2017)利用阻抗張量分解統(tǒng)計(jì)成像技術(shù),發(fā)現(xiàn)在一條穿過(guò)郯廬斷裂帶南段的壽縣—溧水大地電磁剖面中,中高頻和中低頻的電性結(jié)構(gòu)主軸方位存在很大的系統(tǒng)性差異,推測(cè)研究區(qū)中上地殼和下地殼之間存在結(jié)構(gòu)不整合性,不適合頻率整體性的二維反演;在此基礎(chǔ)上,郭春玲和陳小斌(2018)發(fā)展了交錯(cuò)結(jié)構(gòu)大地電磁二維反演技術(shù),結(jié)果發(fā)現(xiàn)郯廬斷裂帶具有聯(lián)通地幔的高導(dǎo)體.然而,對(duì)于上下主軸差異很大的地區(qū),這種交錯(cuò)模型二維反演技術(shù)是否能得到合理可靠的結(jié)果,還需要驗(yàn)證.綜上,目前已開(kāi)展的工作主要集中在研究區(qū)以北和南側(cè),且都是通過(guò)大地電磁二維反演技術(shù)獲得單條剖面的結(jié)果,對(duì)比性不強(qiáng).
由于數(shù)學(xué)上的復(fù)雜性和計(jì)算量問(wèn)題,大地電磁數(shù)據(jù)長(zhǎng)期處于二維反演解釋階段,但在實(shí)際條件下,尤其在構(gòu)造復(fù)雜區(qū),大地介質(zhì)多是三維結(jié)構(gòu)的.最近幾年,隨著三維反演技術(shù)的蓬勃發(fā)展,三維帶地形反演已經(jīng)成熟,有效提高了大地電磁探測(cè)結(jié)果的可靠性和可解釋性.本文中,我們利用覆蓋郯廬斷裂帶合肥—宿遷段的較為均勻的大地電磁陣列數(shù)據(jù),通過(guò)三維反演技術(shù)獲得了郯廬斷裂帶合肥—宿遷段高精度電性結(jié)構(gòu),為深入認(rèn)識(shí)郯廬斷裂帶合肥—宿遷段的構(gòu)造演化、分段性和可能的發(fā)震構(gòu)造提供深部電性結(jié)構(gòu)依據(jù).
郯廬斷裂帶合肥—宿遷段位于郯廬斷裂帶的中南段(圖1a),其西側(cè)為華北板塊,東側(cè)的蘇魯造山帶和揚(yáng)子板塊以淮陰—響水?dāng)嗔褳榻缦噙B.蘇魯造山帶內(nèi)部又以海泗斷裂為界,分成膠南榴輝巖帶(超高壓變質(zhì)帶)和張八嶺藍(lán)片巖帶(高壓變質(zhì)帶).張八嶺隆起緊鄰郯廬斷裂帶東側(cè),其中分布的藍(lán)片巖帶北與蘇魯造山帶南部的藍(lán)片巖帶相連,南與大別造山帶東側(cè)的藍(lán)片巖帶相連(Xu et al.,2009; 徐樹(shù)桐等,1994;圖1b).郯廬斷裂帶西側(cè)的華北板塊中,由北向南分布著徐淮弧形逆沖推覆構(gòu)造、固鎮(zhèn)盆地、蚌埠隆起、肥東拗陷、合肥盆地等多個(gè)構(gòu)造單元(安徽省地質(zhì)礦產(chǎn)局, 1987).斷裂帶兩側(cè)構(gòu)造的主要走向差異明顯,東側(cè)的主要斷裂如淮陰—響水?dāng)嗔选⒑闈伞獪隙諗嗔?、海泗斷裂、滁河斷裂、江浦—六合斷裂和江南斷裂多以北東走向?yàn)橹鳎鱾?cè)的渦河斷裂、肥中斷裂等多以東西走向?yàn)橹?圖1c).
郯廬斷裂帶宿遷以北的沂沭段由4~5條主干斷裂組成,自東向西分別是昌邑—大店斷裂、安丘—莒縣斷裂、白芬子—浮來(lái)山斷裂、沂水—湯頭斷裂和鄌郚—葛溝斷裂(劉保金等,2015;張鵬等,2010).沂水—湯頭斷裂和鄌郚—葛溝斷裂被新近系地層覆蓋,地震活動(dòng)性較弱;安丘—莒縣斷裂是新生代以來(lái)產(chǎn)生的一條發(fā)震斷裂(晁洪太等, 1997),與昌邑—大店斷裂在第四紀(jì)仍有右行走滑活動(dòng),地震危險(xiǎn)性較大(王華林等,2017).江蘇段主要由4條主干斷裂組成,自東向西分別是山左口—泗洪斷裂、新沂—新店斷裂、墨河—陵城斷裂和紀(jì)集—王集斷裂,其中,山左口—泗洪斷裂和紀(jì)集—王集斷裂向北分別與昌邑—大店斷裂和鄌郚—葛溝斷裂相連(王小鳳等,2000;張鵬等,2011).嘉山以南由西側(cè)的池河—太湖斷裂和東側(cè)的嘉山—廬江斷裂組成(湯加富和許衛(wèi), 2002; 張交東等, 2010).研究區(qū)南北兩側(cè)地形存在差異,嘉山以北以平原地形為主,多被第四紀(jì)構(gòu)造覆蓋;嘉山以南地形起伏明顯,含蚌埠隆起、張八嶺隆起、肥東坳陷等構(gòu)造(圖1c).
研究區(qū)近30年(1990—2020)地震活動(dòng)性較弱,以中小震為主,沒(méi)有5級(jí)以上地震發(fā)生(朱艾斕等,2018;圖1c).這些地震主要分布在斷裂帶兩側(cè)的華北板塊和揚(yáng)子板塊上,郯廬斷裂主干上地震較少,只有在肥中斷裂與郯廬斷裂交界處地震較集中.
圖1 研究區(qū)構(gòu)造簡(jiǎn)圖及測(cè)點(diǎn)位置圖
我們?cè)陂L(zhǎng)約315 km的郯廬斷裂帶合肥—宿遷段及其兩側(cè)區(qū)域完成了180個(gè)寬頻帶大地電磁測(cè)點(diǎn)的數(shù)據(jù)采集工作.這些測(cè)點(diǎn)構(gòu)成的大地電磁臺(tái)陣基本均勻地分布在整個(gè)研究區(qū),覆蓋了斷裂帶兩側(cè)寬約120 km的范圍,一些測(cè)點(diǎn)沿NWW-NEE方向延伸超過(guò)240 km,平均測(cè)點(diǎn)間距約10 km,測(cè)點(diǎn)位置如圖1c所示.
所有測(cè)點(diǎn)的野外數(shù)據(jù)采集工作均使用加拿大Phoenix公司生產(chǎn)的V5-2000寬頻帶大地電磁測(cè)深系統(tǒng)完成,測(cè)點(diǎn)布置采用“十”字形沿正南北東西向布極,X軸指向磁北方向,Y軸指向磁東方向,Z軸指向垂向.記錄包含2個(gè)相互垂直的水平電場(chǎng)分量和3個(gè)兩兩垂直的磁場(chǎng)分量的時(shí)間序列,每個(gè)測(cè)點(diǎn)觀測(cè)超過(guò)36 h.對(duì)受干擾的測(cè)點(diǎn)采用了偏移復(fù)測(cè)、延長(zhǎng)采集時(shí)間等方式改善觀測(cè)質(zhì)量,保證測(cè)點(diǎn)的有效周期在2000 s以上.在數(shù)據(jù)處理中,對(duì)于主要受隨機(jī)噪聲干擾測(cè)點(diǎn)采用遠(yuǎn)參考道和Robust方法(Parker and Booker, 1996)進(jìn)行處理,對(duì)于主要受強(qiáng)干擾人文噪聲測(cè)點(diǎn)采用遠(yuǎn)參考道和非Robust方法(張赟昀和陳小斌,2020)進(jìn)行處理,這些技術(shù)有效地抑制了噪聲干擾,提高了數(shù)據(jù)質(zhì)量.使用大地電磁數(shù)據(jù)管理、處理和分析軟件MTPioneer(陳小斌等,2004a)建立測(cè)區(qū)大地電磁數(shù)據(jù)集,然后進(jìn)行詳細(xì)處理、分析.
觀測(cè)獲得的視電阻率和阻抗相位數(shù)據(jù)是地下深部結(jié)構(gòu)的原始響應(yīng).分析不同構(gòu)造單元上大地電磁測(cè)點(diǎn)的視電阻率和阻抗相位曲線特征,可以初步定性分析各構(gòu)造單元的電性結(jié)構(gòu)特征.圖2給出研究區(qū)不同構(gòu)造單元部分測(cè)點(diǎn)的視電阻率和相位曲線,具體測(cè)點(diǎn)位置在圖1c中用綠色三角形表示.
圖2 郯廬斷裂帶合肥—宿遷段及鄰區(qū)典型MT測(cè)點(diǎn)的視電阻率和相位曲線圖
郯廬斷裂帶西側(cè),北部(測(cè)點(diǎn)1,2)和南部(測(cè)點(diǎn)5,6)的視電阻率值較小,基本上在100 Ωm以下,曲線較為平緩.中間兩個(gè)位于蚌埠隆起區(qū)的測(cè)點(diǎn)(測(cè)點(diǎn)3,4)視電阻率值高于南北兩側(cè)測(cè)點(diǎn),低頻和高頻視電阻率和相位曲線起伏較大,與蚌埠隆起區(qū)較復(fù)雜的構(gòu)造特征相對(duì)應(yīng),表明蚌埠隆起區(qū)下方可能存在局部高阻體.郯廬斷裂帶內(nèi)部,南端(測(cè)點(diǎn)11,12)和北端(測(cè)點(diǎn)7,8)的視電阻率值相對(duì)較??;中部(測(cè)點(diǎn)9)視電阻率值較大且存在明顯的升高趨勢(shì),這表明郯廬斷裂帶合肥—宿遷段存在明顯的結(jié)構(gòu)不均勻性,中部存在局部高阻體且分布較深.在郯廬斷裂帶東側(cè),除測(cè)點(diǎn)18和測(cè)點(diǎn)13的xy模式外,各測(cè)點(diǎn)視電阻率曲線均存在從高頻到低頻逐漸升高的趨勢(shì),表明郯廬斷裂帶東側(cè)深部可能存在大范圍的高阻異常體,與郯廬斷裂帶東側(cè)的蘇魯造山帶、張八嶺隆起相對(duì)應(yīng).阻抗相位值小與大基本上與視電阻率的起與伏相對(duì)應(yīng).
四個(gè)不同頻率下的視電阻率和阻抗相位的等頻率分布情況如圖3所示.阻抗相位反映視電阻率的變化,一般高相位表示視電阻率相對(duì)較低,低相位表示視電阻率相對(duì)較高.蘇魯造山帶地區(qū)阻抗相位在3 Hz時(shí)很低,與該區(qū)域視電阻率在高頻段迅速升高對(duì)應(yīng).在低頻時(shí)相位處于中值,與深部此區(qū)域視電阻率達(dá)到平穩(wěn)狀態(tài)對(duì)應(yīng).郯廬斷裂帶西側(cè)在中低頻阻抗相位存在高值,表明此區(qū)域視電阻率在深部下降,可能存在高導(dǎo)體.
圖3 視電阻率、阻抗相位等頻率分布圖
視電阻率是對(duì)地下不均勻電性結(jié)構(gòu)的均勻等效,從分布情況可以分析地下結(jié)構(gòu)的變化.在全部四個(gè)頻段上張八嶺隆起和蚌埠隆起東側(cè)呈現(xiàn)高值,合肥盆地為低值.淮陰—響水?dāng)嗔褞б员钡奶K魯造山帶在中低頻呈現(xiàn)高值.在0.00293Hz視電阻率等頻率圖中,郯廬斷裂帶東側(cè)連續(xù)存在顯著的高異常,表明郯廬斷裂帶兩側(cè)物質(zhì)在深部存在明顯的電性差異,郯廬斷裂帶是一條明顯的電性分界.
相位張量法(Caldwell et al., 2004)沒(méi)有預(yù)設(shè)的二維性假設(shè)條件,用其進(jìn)行阻抗張量分解可以有效消除局部小異常體的影響,得到的二維偏離度(β)和相位張量不變量(φ2)參數(shù)來(lái)分析深部結(jié)構(gòu)的維性和優(yōu)勢(shì)走向等信息.用橢圓表示相位張量時(shí),橢圓的長(zhǎng)短軸表示構(gòu)造的兩個(gè)可能走向,顏色代表二維偏離度角和相位張量不變量.通過(guò)這些信息不僅可以定性分析深部結(jié)構(gòu)特征,也可以印證反演結(jié)果的可靠性.圖4給出了通過(guò)相位張量分解獲得的二維偏離度(圖4a—4d)、相位張量不變量(圖4e—4h)和實(shí)感應(yīng)矢量(圖4i—4l)的等頻率分布圖.
圖4 相位張量橢圓和感應(yīng)矢量
一般認(rèn)為二維偏離度β小于3°表示介質(zhì)接近一維或者二維,β越大表示此區(qū)域越趨于三維結(jié)構(gòu)(Caldwell et al., 2004; Booker et al., 2014).從圖4a—4d展示的四個(gè)不同頻率的相位張量結(jié)果,可以看出在3 Hz二維偏離度普遍小于3°,僅有個(gè)別測(cè)點(diǎn)二維偏離度大于3°.同時(shí)相位張量橢圓長(zhǎng)短軸接近,說(shuō)明淺部結(jié)構(gòu)表現(xiàn)近一維.在0.56 Hz時(shí)沿郯廬斷裂帶宿遷—嘉山段存在高值,表明沿該段可能存在三維結(jié)構(gòu).在0.0234 Hz時(shí),嘉山以北橢圓長(zhǎng)軸普遍指向NE向,且以淮陰—響水?dāng)嗔鸭把酉驗(yàn)榻纾辈慷S偏離度普遍小于3°,南部二維偏離度大于3°的分布.在0.00293 Hz,除了蘇魯造山帶內(nèi)部,其他區(qū)域都呈現(xiàn)高值.綜合來(lái)看,在中頻段和低頻段都存在β較大的區(qū)域,研究區(qū)的整體構(gòu)造較復(fù)雜,深部以三維結(jié)構(gòu)為主.
相位張量不變量(φ2)是相位張量最大值和最小值的幾何平均值(Heise et al., 2008),是衡量電阻率隨深度變化趨勢(shì)的參數(shù),高φ2對(duì)應(yīng)電阻率值下降,一般情況下對(duì)應(yīng)于相對(duì)低阻區(qū),低φ2對(duì)應(yīng)電阻率值升高,對(duì)應(yīng)于相對(duì)高阻區(qū).從各個(gè)測(cè)點(diǎn)的φ2值隨頻率的分布來(lái)看(圖4e—4h),郯廬斷裂帶西側(cè),高頻段φ2值較低,低頻段φ2值較高,說(shuō)明西側(cè)淺部電阻率較高,深部電阻率下降,斷裂帶兩側(cè)φ2值差異性大,說(shuō)明兩側(cè)電性結(jié)構(gòu)有明顯不同.
感應(yīng)矢量由傾子中兩個(gè)復(fù)數(shù)的實(shí)部組成,模的大小表示電阻率橫向變化的梯度,方向指向電流匯聚方向(Parkinson,1962; 陳小斌等,2004b).在二維各向同性介質(zhì)中,感應(yīng)矢量指向構(gòu)造的傾向.圖4i—4l顯示,3 Hz時(shí)感應(yīng)矢量幅值較小,方向亂,表明淺層可能受小異常體影響,電性結(jié)構(gòu)復(fù)雜,低頻時(shí)各測(cè)點(diǎn)幅值變大,方向沿不同地質(zhì)塊體表現(xiàn)不同特性,表明深部為電性結(jié)構(gòu)差異明顯的構(gòu)造塊體.
通過(guò)對(duì)視電阻率、阻抗相位、相位張量橢圓和感應(yīng)矢量的分析,初步定性認(rèn)識(shí)了郯廬斷裂帶合肥—宿遷段的電性結(jié)構(gòu)特征.個(gè)別地質(zhì)單元的走向單一,但整個(gè)研究區(qū)域以三維為主.想要對(duì)區(qū)域深部電性結(jié)構(gòu)進(jìn)行更精細(xì)認(rèn)識(shí),需要利用大地電磁三維反演技術(shù)對(duì)數(shù)據(jù)集進(jìn)行反演,獲得該區(qū)域三維電阻率結(jié)構(gòu)模型.
大地電磁野外布極是沿著南北和東西方向進(jìn)行的,譜數(shù)據(jù)的編輯工作也基于這兩個(gè)方向進(jìn)行,這種情況下旋轉(zhuǎn)網(wǎng)格和數(shù)據(jù)會(huì)導(dǎo)致數(shù)據(jù)質(zhì)量變差,影響反演擬合.因此本文使用沒(méi)有進(jìn)行旋轉(zhuǎn)的數(shù)據(jù)進(jìn)行反演.在大地電磁阻抗張量中,反對(duì)角元素對(duì)應(yīng)于大地電磁響應(yīng)主分量元素,從一維、二維反演時(shí)代開(kāi)始,主分量元素特征經(jīng)過(guò)長(zhǎng)時(shí)期的研究分析,質(zhì)量好壞易于甄別.雖然對(duì)角元素(次分量元素)也包含了部分深部結(jié)構(gòu)信息,但其在一維、二維模型中不存在,故其響應(yīng)特征目前還不是很了解,數(shù)據(jù)較難進(jìn)行編輯和選擇.相對(duì)于直接對(duì)阻抗數(shù)據(jù)進(jìn)行選擇,對(duì)視電阻率和阻抗相位單獨(dú)進(jìn)行數(shù)據(jù)選擇,可以保留更多數(shù)據(jù).同時(shí),由于本研究中所用的大地電磁臺(tái)陣分布較密且均勻,只用主分量元素(xy和yx分量)也可以為反演提供足夠的約束.
反演中使用了有效頻率周期為320~0.000067 Hz之間的46個(gè)頻點(diǎn)的數(shù)據(jù).反演前利用Roplus一維反演方法(Parker and Booker, 1996)對(duì)視電阻率和阻抗相位的匹配度進(jìn)行檢測(cè),剔除影響視電阻率和阻抗相位匹配的數(shù)據(jù).
模型網(wǎng)格劃分的方向和數(shù)據(jù)觀測(cè)時(shí)的方位一致,為正南北東西方向,初始模型網(wǎng)格的數(shù)量為Nx×Ny×Nz=81×125×128.核心區(qū)y軸網(wǎng)格間距為2.5 km,x軸網(wǎng)格間距為5.9 km,核心區(qū)以外水平網(wǎng)格以1.5的比例因子向外擴(kuò)展10個(gè)網(wǎng)格.垂直網(wǎng)格單層最大厚度為50 km.為了讓模型更接近實(shí)際,模型中還加入了13個(gè)地形層以模擬實(shí)際地形,每層厚度30 m.模型網(wǎng)格如圖5所示.
圖5 三維反演網(wǎng)格
本研究使用基于非線性共軛梯度法的ModEM算法(Egbert and Kelbert, 2012)進(jìn)行三維反演研究工作.大地電磁三維反演軟件系統(tǒng)toPeak(Liu et al.,2018)為本次研究中三維反演數(shù)據(jù)的管理,帶地形的初始模型建立,反演參數(shù)的設(shè)置等工作提供了極大便利.反演中視電阻率、阻抗相位所用的門(mén)檻誤差分別為5%、1.4324.每一輪反演中,正則化因子從1000開(kāi)始,當(dāng)均方根誤差的相對(duì)變化小于0.0025時(shí),以0.2的比例因子下降,當(dāng)正則化因子小于最小值1×10-8時(shí)反演結(jié)束.
反演中使用了印模重構(gòu)法(葉濤等,2013),以增加初始模型的合理性和反演的穩(wěn)定性.首先使用100 Ωm的均勻半空間作為初始模型進(jìn)行反演并獲得反演結(jié)果,然后通過(guò)印模重構(gòu)保留前一步反演結(jié)果宏觀輪廓,去掉模型中的一些極值,使其作為下一步反演的初始模型.經(jīng)過(guò)三次反演和兩次印模重構(gòu),分別迭代了96、64、44次,最終的均方根誤差(RMS misfit)為1.962,獲得了郯廬斷裂帶合肥—宿遷段三維電阻率模型(圖6).圖7顯示了研究區(qū)南、北兩條較長(zhǎng)剖面上的擬斷面圖,從圖中可以看出兩條剖面上觀測(cè)獲得的視電阻率和阻抗相位與基于結(jié)果模型計(jì)算獲得的視電阻率和阻抗相位十分相似,表明數(shù)據(jù)擬合的較好.
圖6 郯廬斷裂帶合肥-宿遷段立體電性結(jié)構(gòu)分布圖
圖7 兩條長(zhǎng)剖面上觀測(cè)和基于結(jié)果模型計(jì)算的視電阻率和阻抗相位擬斷面圖
我們也嘗試了利用其它不同類(lèi)型的數(shù)據(jù)組合(視電阻率和阻抗相位加傾子、全阻抗數(shù)據(jù)、全阻抗加傾子數(shù)據(jù))進(jìn)行反演,分別迭代了86、62、71次,均方根誤差分別為2.407、7.98、6.13.對(duì)比這些數(shù)據(jù)類(lèi)型的反演結(jié)果后發(fā)現(xiàn),結(jié)果模型之間十分近似,其中視電阻率和阻抗相位數(shù)據(jù)反演的深部分辨率相對(duì)較高,數(shù)據(jù)擬合的最好.(Newman et al, 2008)曾發(fā)現(xiàn)因?yàn)閤x和yy元素?cái)?shù)據(jù)基數(shù)較低,信噪比較低,使用全阻抗數(shù)據(jù)組合反演會(huì)降低反演性能.本文研究區(qū)位于我國(guó)東部發(fā)達(dá)地區(qū),環(huán)境干擾較多,數(shù)據(jù)處理時(shí)功率譜的選擇是以xy和yx分量為標(biāo)準(zhǔn),對(duì)于xx和yy分量未做任何有效選擇,導(dǎo)致噪聲影響很多隱藏在xx和yy分量中.目前關(guān)于對(duì)角元素?cái)?shù)據(jù)正確形態(tài)的認(rèn)識(shí)尚不明確,反演頻點(diǎn)的選取難以把握,可能會(huì)影響反演的可靠性.同時(shí)過(guò)多的數(shù)據(jù)可能會(huì)引起反演中的欠擬合,使數(shù)據(jù)無(wú)法有效地約束模型.而研究區(qū)的均勻面狀測(cè)點(diǎn)分布,可以彌補(bǔ)由于缺少對(duì)角元素帶來(lái)的約束不夠,我們認(rèn)為選用xy和yx反對(duì)角元素進(jìn)行三維反演是比較穩(wěn)妥的辦法.因此本文選用視電阻率和阻抗相位反演結(jié)果作為最終結(jié)果.
大地電磁探測(cè)深度受到電磁波趨膚深度和地下電性結(jié)構(gòu)的共同影響,在進(jìn)行電性特征分析和地質(zhì)解釋之前,有必要進(jìn)行可靠性測(cè)試,以確定模型受觀測(cè)數(shù)據(jù)約束的情況.嘉山附近的相對(duì)低阻區(qū)LR分隔了南北兩側(cè)的高阻體R1和R2,合肥盆地深部存在低阻異常體C2,它們分布在模型核心區(qū)域,對(duì)認(rèn)識(shí)該地區(qū)的深部結(jié)構(gòu)有非常重要的意義.如圖8所示,對(duì)LR和C2進(jìn)行可靠性測(cè)試.
圖8 三維電性結(jié)構(gòu)可靠性驗(yàn)證區(qū)域與曲線擬合情況
在LR位置將深度范圍3~30 km(圖8a)中小于1000 Ωm的電阻率值修改為與周?chē)咦梵w電阻率值相當(dāng)?shù)?000 Ωm,修改后的電阻率模型如圖8b所示.在C2位置將深度范圍20~60 km(圖8c)中小于100 Ωm的電阻率值修改為與鄰區(qū)接近的100 Ωm,修改后的電阻率模型如圖8d所示.
對(duì)修改后的模型進(jìn)行三維正演,得到預(yù)測(cè)的視電阻率和阻抗相位響應(yīng)曲線.從圖8e可以看出,模型修改后嘉山區(qū)域相對(duì)低阻LR里幾個(gè)測(cè)點(diǎn)的視電阻率在中頻段開(kāi)始出現(xiàn)抬升,阻抗相位在中頻段下降,擬合程度與模型修改前相比變差,表明我們得到的電阻率模型中相對(duì)低阻體LR是受到數(shù)據(jù)良好約束的可靠結(jié)構(gòu).對(duì)C2修改后,其范圍內(nèi)的幾個(gè)測(cè)點(diǎn)在低頻段預(yù)測(cè)數(shù)據(jù)與觀測(cè)數(shù)據(jù)的擬合相比與模型修改前也變差了,表明電阻率模型中低阻體C2也是被觀測(cè)數(shù)據(jù)約束的可靠結(jié)構(gòu).
上述三維反演結(jié)果和可靠性分析,表明我們得到了郯廬斷裂帶合肥—宿遷段較為可靠的三維電阻率模型.從三維電性立體圖可以看出沿郯廬斷裂帶走向從淺至深呈現(xiàn)高阻特征(圖6a),嘉山附近是郯廬斷裂帶內(nèi)的相對(duì)低阻區(qū)LR;同時(shí)也發(fā)現(xiàn)研究區(qū)內(nèi)郯廬斷裂帶東西兩側(cè)電性結(jié)構(gòu)有明顯差異,郯廬斷裂帶作為蘇魯高阻體R1與張八嶺高阻體R2的西邊界,向下延伸較深,傾角近直立(圖6b).此外,反演結(jié)果還揭示了斷裂帶東側(cè)高阻由南(R2)向北(R1)呈現(xiàn)加深加寬的形態(tài);西側(cè)淺部高低阻分布不均,下地殼高導(dǎo)體C1與C2較明顯.反演結(jié)果的主要電性特征與前文數(shù)據(jù)定性分析的認(rèn)識(shí)基本一致,印證了反演結(jié)果的可靠性.
由電性水平截面可看出(圖9),研究區(qū)淺部受小斷裂和構(gòu)造地塊發(fā)育影響高低阻相間,而深部逐漸形成以郯廬斷裂帶為界的西低東高的電性結(jié)構(gòu).在15 km以上深度,郯廬斷裂帶兩側(cè)均有高阻分布,而到15 km深度以下,郯廬斷裂帶東側(cè)的高阻結(jié)構(gòu)R1和R2一直顯著,斷裂帶西側(cè)的高阻逐漸消減,并在嘉山以北沿郯廬斷裂帶西側(cè)出現(xiàn)了一條明顯的高導(dǎo)異常C1,西南側(cè)出現(xiàn)明顯的高導(dǎo)體C2,且C2在30 km深度分布最廣.在40 km深度,嘉山以北斷裂帶西側(cè)的高導(dǎo)條帶C1消減,合肥盆地下方的高導(dǎo)體C2、蘇魯造山帶高阻體R1和張八嶺高阻體R2依然存在.值得注意的是在5~30 km的深度上,郯廬斷裂帶東側(cè)蘇魯造山帶內(nèi)的高阻體R1與張八嶺高阻體R2之間,嘉山周?chē)嬖谝黄鄬?duì)低阻區(qū)LR將兩個(gè)高阻分離.
圖9 郯廬斷裂帶合肥—宿遷段電性結(jié)構(gòu)平面圖
郯廬斷裂帶在深部主要作為蘇魯高阻體R1和張八嶺高阻體R2的西邊界出現(xiàn)(圖10).AA′截面中,R1邊界從淺到深靠近郯廬斷裂帶中東側(cè)位置,淮陰響水?dāng)嗔炎鳛樘K魯造山帶的南邊界,分隔了北側(cè)中上地殼的高阻體與南側(cè)的高導(dǎo)異常.BB′截面中R1邊界向西陡傾,淺部靠近斷裂帶東側(cè),深部靠近斷裂帶西側(cè);郯廬斷裂帶西側(cè)10 km深度以上有高阻出現(xiàn),低阻體C1分布深度在12~40 km之間.CC′截面中,R1邊界從淺到深靠近斷裂帶的西側(cè),截面西側(cè)20 km深度以上有低阻C1和高阻分布.三個(gè)截面中,R1寬度在40~60 km之間,最大深度超過(guò)50 km,電阻率在500Ωm以上;C1分布寬度在40 km以?xún)?nèi).DD′截面穿過(guò)高阻體R1與R2分隔區(qū)域,斷裂帶下方呈現(xiàn)相對(duì)低阻(低于R1和R2),深部的高阻分布可能是受高阻體R1和R2的影響,西側(cè)高阻體R3分布在20 km以上深度;EE′截面中,張八嶺高阻R2邊界從淺到深靠近郯廬斷裂帶的東側(cè),西側(cè)深部為高導(dǎo)體C2;FF′截面中高阻體R2西邊界在淺部靠近郯廬斷裂帶西側(cè),西側(cè)20 km以上呈現(xiàn)高阻,深部為高導(dǎo)C2,在向東接近郯廬斷裂帶時(shí)C2向上延伸至淺部.最南側(cè)兩個(gè)截面中,R2寬度在20 km左右,分布深度變化范圍較大,電阻率在30 km以上深度超過(guò)500 Ωm,更深處則逐漸變低;C2的最大分布寬度超過(guò)100 km.
圖10 郯廬斷裂帶合肥—宿遷段電性結(jié)構(gòu)垂直剖面圖
從上面電性結(jié)構(gòu)模型中,可以看出研究區(qū)存在4個(gè)大的電性分界或陡變帶,電性結(jié)構(gòu)有6個(gè)典型電性分區(qū),它們的特征分別是:
郯廬斷裂帶:是華北板塊與揚(yáng)子板塊的構(gòu)造單元邊界,具有南北分段特點(diǎn)(王小鳳等,2000).在深部電性上呈現(xiàn)高阻與低阻之間的陡立邊界,是典型的平移斷裂特征,也是切割深度超過(guò)地殼的深大斷裂.在斷裂帶內(nèi)嘉山南北有明顯不同,嘉山以北主斷裂由4~5支組成,在電性結(jié)構(gòu)上對(duì)應(yīng)北部較大的高阻R1,嘉山以南由2支組成,電性結(jié)構(gòu)上對(duì)應(yīng)南部較窄的高阻R2,兩個(gè)高阻體之間的嘉山地區(qū)是相對(duì)低阻區(qū)LR.依據(jù)獲得的電性結(jié)構(gòu),嘉山地區(qū)是郯廬斷裂帶中段和南段的分段點(diǎn),與一些學(xué)者從構(gòu)造習(xí)性、演化歷史、地震成像等方面的分段結(jié)果一致(晁洪太等,1997;熊振等,2016).
海泗斷裂:是蘇魯造山帶中超高壓變質(zhì)帶(西北側(cè))和高壓變質(zhì)帶(東南側(cè))的分界線.從水平切片圖9看,5 km深度一個(gè)條帶狀高阻體沿海泗斷裂西側(cè)分布;10 km深度,沿海泗斷裂的高阻條帶逐漸消減.
淮陰—響水?dāng)嗔眩菏翘K魯造山帶與楊子板塊相接斷裂,在嘉山周?chē)c郯廬斷裂帶交匯.從AA′截面看,是蘇魯造山帶高阻與蘇北盆地低阻的分界,產(chǎn)狀較緩,向NE緩傾.由于嘉山周邊測(cè)點(diǎn)密度不夠,蘇魯造山帶與張八嶺兩個(gè)高阻體在此區(qū)域的切斷,是否和淮陰—響水?dāng)嗔延嘘P(guān)有待驗(yàn)證.
潁上—定遠(yuǎn)斷裂:呈北西西方向展布,處于郯廬斷裂帶西側(cè)并與其垂直,為徐淮地塊與合肥盆地的分界斷裂.在重力和磁場(chǎng)勘察結(jié)果中,均表現(xiàn)為顯著的異常梯度帶(安徽省地質(zhì)礦產(chǎn)局,1987).DD′截面沿該斷裂截取,處于郯廬主斷裂由4條變2條的關(guān)鍵區(qū)域,在電性結(jié)構(gòu)垂直剖面圖中斷裂下方中上地殼呈現(xiàn)高阻特征,是否阻隔了西側(cè)2條主斷裂的繼續(xù)南延值得進(jìn)一步思考.
蘇魯造山帶:西側(cè)以郯廬斷裂為界與華北地塊相鄰,南側(cè)以淮陰—響水?dāng)嗔褳榻缗c楊子板塊相接.蘇魯造山帶西側(cè)下方為高阻結(jié)構(gòu),其深度范圍超過(guò)莫霍面,與研究區(qū)北部蘇魯造山帶高阻深度范圍超過(guò)80 km的結(jié)果(肖騎彬等,2008;翁愛(ài)華等,2018)比較一致.
張八嶺隆起:緊鄰郯廬斷裂帶東側(cè)寬約20 km的高阻條帶R2與張八嶺隆起對(duì)應(yīng),從淺至深一直比較顯著(圖9).在郯廬斷裂帶向南延伸的合肥—廣濟(jì)段的深部電性結(jié)構(gòu)研究結(jié)果中,張八嶺隆起下方的高阻條帶與大別造山帶中上地殼的高阻體相連(崔騰發(fā),2019).從巖性結(jié)果(Xu et al.,2009; 徐樹(shù)桐等,1994)也表明張八嶺隆起在變形前可能是大別—蘇魯碰撞帶的一部分.張八嶺高阻條帶R2可能是在郯廬斷裂帶強(qiáng)烈的左旋走滑運(yùn)動(dòng)拖拽下,彎曲變形并沿著斷裂帶東側(cè)分布.
嘉山盆地:夾持于郯廬斷裂帶內(nèi),南端以潁上—定遠(yuǎn)斷裂為界與合肥盆地東北隅相接(韓樹(shù)棻,1996).嘉山盆地南部在電性上呈現(xiàn)相對(duì)低阻區(qū)LR(DD′截面),阻隔了張八嶺隆起與蘇魯造山帶下方的高阻體,控制了郯廬斷裂帶結(jié)構(gòu)不均勻性,此處可能是張八嶺條帶構(gòu)造在變形過(guò)程中的一個(gè)薄弱區(qū)域.作為高阻體之間的相對(duì)低阻區(qū),其強(qiáng)度低于南北兩側(cè)的高阻體,在新構(gòu)造時(shí)期可能阻礙了斷層中應(yīng)力與應(yīng)變的連續(xù)傳遞,造成了斷層活動(dòng)在嘉山南北兩側(cè)的分段以及斷層幾何結(jié)構(gòu)的變化(徐嘉煒和馬國(guó)鋒, 1992; 張鵬等, 2007).
蚌埠隆起:位于郯廬斷裂帶西側(cè),基底為五河群,屬低角閃巖相區(qū)域變質(zhì)巖系(安徽省地質(zhì)礦產(chǎn)局, 1987).其在深部對(duì)應(yīng)郯廬斷裂帶西側(cè)的高阻體R3(DD′截面),R3深度與橫波速度結(jié)構(gòu)中蚌埠隆起下方約為20 km深度的高速異常體吻合(孟亞鋒等, 2019).
固鎮(zhèn)盆地:位于蘇魯造山帶和郯廬斷裂帶西側(cè),在電性結(jié)構(gòu)中位置與低阻體C1對(duì)應(yīng).
合肥盆地:東界為郯廬斷裂帶的嘉山—廬江斷裂,北界為潁上—定遠(yuǎn)斷裂,因不同時(shí)期盆地?cái)U(kuò)張、分解和遷徙,發(fā)育范圍變化很大(韓樹(shù)棻,1996).其中下地殼與電性結(jié)構(gòu)中的高導(dǎo)體C2對(duì)應(yīng),與此處地震結(jié)果下地殼低速異常一致(Zheng et al., 2017),可能是巖漿作用造成地殼發(fā)生熔融或破碎所致.高導(dǎo)體C2在肥中斷裂與郯廬斷裂帶之間延伸至淺部,同時(shí)是地震活動(dòng)比較密集區(qū)域.合肥盆地中部電性結(jié)構(gòu)成層性明顯(FF′截面).
從研究區(qū)三維電性結(jié)構(gòu)來(lái)看,以主要斷裂帶為界,電性結(jié)構(gòu)呈現(xiàn)明顯橫向分塊特征,并且高低阻分布特征與區(qū)域構(gòu)造基本一致,不同地塊的電性結(jié)構(gòu)存在顯著差異.郯廬斷裂帶西側(cè)華北板塊中上地殼的電性結(jié)構(gòu)較為破碎,這可能和華北克拉通東部乃至整個(gè)中國(guó)東部在早白堊世的伸展活動(dòng)和強(qiáng)烈的巖漿活動(dòng)有關(guān)(Wu et al., 2005; Zhu et al., 2012).軟流圈熱流上涌促使地殼底部發(fā)生底侵,引起地殼巖石的部分熔融,可能是高導(dǎo)體C1和C2形成的主要原因.而張八嶺高阻條帶和蘇魯造山帶高阻體分布較深,在伸展活動(dòng)和巖漿活動(dòng)中未被明顯改造,從而形成了郯廬斷裂帶兩側(cè)明顯的電性結(jié)構(gòu)差異.研究區(qū)以北郯廬斷裂帶山東段的大地電磁三維反演結(jié)果,同樣顯示了斷裂帶兩側(cè)西低東高的電阻率特征(Ye et al., 2021).
研究區(qū)域內(nèi)的地震以淺源小震為主,斷裂帶兩側(cè)的華北與揚(yáng)子板塊內(nèi)部地震多于郯廬斷裂帶主干斷裂(圖1c).小震一般發(fā)生在高阻體的邊緣或其它地區(qū),高阻體內(nèi)部很少有地震發(fā)生(圖9).一般認(rèn)為高阻巖石具有較高強(qiáng)度,不易變形,在積累較大應(yīng)力后可以產(chǎn)生脆性破裂;低阻的巖石相對(duì)破碎,強(qiáng)度較弱,積累應(yīng)力能力弱(趙國(guó)澤等,2008; Cai et al., 2017).因此區(qū)域內(nèi)的地震主要發(fā)生在高阻體與周?chē)鷺?gòu)造相互作用的邊界上.主干斷裂上的地震較少,只在肥中斷裂與郯廬斷裂帶交匯處地震較集中,但是對(duì)斷層現(xiàn)今形變的研究結(jié)果表明,郯廬斷裂帶合肥—宿遷段斷層b值低(朱艾斕等, 2018),可能閉鎖程度高,具備了較高的應(yīng)力水平(李彥川等,2016; 李臘月等,2020).
斷層的地震危險(xiǎn)性不僅和現(xiàn)今的變形特征有關(guān),還受深部結(jié)構(gòu)的制約.郯廬斷裂帶在深部作為蘇魯造山帶高阻體R1的西邊界,傾角近直立,與地震反射剖面的結(jié)果非常一致(劉保金,2015).震源機(jī)制解顯示宿遷—嘉山段主壓應(yīng)力方向(NWW向)垂直于斷裂的走向(朱艾斕等,2018),斷面摩擦力增大,易于應(yīng)力積累,造成斷層閉鎖.在宿遷—泗洪段,高阻體R1邊界位于斷裂帶東側(cè);而在泗洪—嘉山段,蘇魯造山帶高阻體R1南端向西發(fā)生偏轉(zhuǎn),郯廬斷裂帶穿過(guò)了高阻體R1.在最近一期的GPS觀測(cè)中,泗洪—嘉山段斷層的平行斷層滑動(dòng)虧損率明顯高于宿遷—泗洪段,而兩段具有相近的垂直斷層滑動(dòng)虧損率(李臘月等,2020),這可能與凸出的局部高阻結(jié)構(gòu)R1對(duì)泗洪—嘉山段斷層平行滑動(dòng)的阻礙有關(guān).當(dāng)斷層面上的剪應(yīng)力超過(guò)斷層的抗剪強(qiáng)度時(shí)會(huì)發(fā)生巨大位移,因此泗洪—嘉山段的地震危險(xiǎn)性值得關(guān)注.郯廬斷裂帶在嘉山以南的主壓應(yīng)力方向以NWW和EW向?yàn)橹?朱艾斕等,2018),局部應(yīng)力狀態(tài)為擠壓兼走滑.郯廬斷裂帶嘉山—合肥段深部為高阻體R2的陡立邊界,但有較高的平行斷層滑動(dòng)虧損率(李臘月等,2020),易于積累應(yīng)力,斷層容易閉鎖.造成此現(xiàn)象原因可能是受高阻體R2向西南方向延伸的具體結(jié)構(gòu)的限制(崔騰發(fā),2019).
本文基于由180個(gè)寬頻帶大地電磁測(cè)點(diǎn)組成的觀測(cè)臺(tái)陣,利用三維反演技術(shù)得到了郯廬斷裂帶合肥—宿遷段的三維電性結(jié)構(gòu),詳細(xì)分析了該區(qū)域電性結(jié)構(gòu)分布情況,并結(jié)合地質(zhì)、大地測(cè)量、地震學(xué)等資料,綜合分析了郯廬斷裂帶合肥—宿遷段電性結(jié)構(gòu)與構(gòu)造分布,地震活動(dòng)孕育環(huán)境的關(guān)系,具體結(jié)論如下:
(1)郯廬斷裂帶是一個(gè)直立型、切穿地殼的深大斷裂帶,切割深度達(dá)莫霍面,淺部電性結(jié)構(gòu)受小斷裂和構(gòu)造發(fā)育影響高低阻相間,而深部逐漸形成以郯廬斷裂帶為界西低東高的兩側(cè)具有明顯差異的電性特征.
(2)從深部電性結(jié)構(gòu)上論證了郯廬斷裂帶具有分段性特點(diǎn).南北兩處高阻體在嘉山附近的分離,引起郯廬斷裂帶深部電性結(jié)構(gòu)明顯的分段變化,高阻體之間的相對(duì)低阻區(qū)阻礙了斷層上應(yīng)力的傳遞,造成了嘉山北、南兩側(cè)斷裂帶構(gòu)造變形強(qiáng)、弱的變化;
(3)沿郯廬斷裂帶東側(cè)分布的蘇魯造山帶高阻體在泗洪—嘉山段向西偏轉(zhuǎn),使得斷裂帶穿過(guò)了高阻體,可能會(huì)成為阻礙斷層活動(dòng)的障礙體,引起應(yīng)力積累,需重點(diǎn)關(guān)注泗洪—嘉山段.
致謝感謝朱艾斕研究員和房立華研究員提供的地震定位結(jié)果.感謝審稿專(zhuān)家和編輯對(duì)本文提出的重要修改意見(jiàn).