高伊航， 沈軍輝， 蘇永軍， 匡海陽

(1.成都理工大學(xué) 環(huán)境與土木工程學(xué)院，成都 610059；2.中國地質(zhì)調(diào)查局 天津地質(zhì)調(diào)查中心，天津 300170)

0 引言

大地電磁測深法(Magnetotelluric，MT)，是以天然電磁場為場源來研究地球內(nèi)部電性結(jié)構(gòu)的一種重要的地球物理手段[1]。隨著國際合作和交流的不斷加深，大地電磁法在深部構(gòu)造的探測中已成為一種必不可少的手段[2]。我國地球科學(xué)研究者對此方法的研究與應(yīng)用付出了很大的努力[3-8]。陳龐龍等[9]將大地電磁測深方法應(yīng)用于深圳市活斷層的探測中，并采用遠(yuǎn)參考道數(shù)據(jù)觀測與處理技術(shù)消除測點附近的電磁噪音影響，運用二維反演獲得了地下電性結(jié)構(gòu)；符超等[10]運用大地電磁測深對卡拉先格爾斷裂進行測量，獲得了關(guān)于沿剖面的地殼二維電性結(jié)構(gòu)，確定了該斷裂的類型及其成礦意義。大地電磁是研究地球深部電性結(jié)構(gòu)與構(gòu)造的主要地球物理方法，被廣泛應(yīng)用于油氣勘探、資源勘查以及深部地球物理調(diào)查等領(lǐng)域，作者通過對研究區(qū)大地電磁數(shù)據(jù)的處理和解釋，查明了研究區(qū)地層結(jié)構(gòu)和隱伏斷裂構(gòu)造的空間展布特征以及研究區(qū)的電場特征，并對該區(qū)構(gòu)造單元進行重新厘定。

圖1 研究區(qū)大地構(gòu)造位置圖Fig.1 The tectonic structure of the study area

1 區(qū)域地質(zhì)概況

研究區(qū)在大地構(gòu)造位置上，處于泰沂山地丘陵北麓山前傾斜平原東段的萊州灣南岸平原，地勢自南向北傾斜，區(qū)內(nèi)地貌類型變化層次清晰，自南向北依次為山前沖洪積平原、沿海低平原和現(xiàn)代潮灘。現(xiàn)代潮灘由海岸線向海區(qū)方向延伸，為粉砂質(zhì)、泥質(zhì)海岸帶，發(fā)育有潮溝系統(tǒng)。在沖洪積平原、河海積平原上分布著黃土及次生黃土。研究區(qū)是華北板塊內(nèi)三大II級構(gòu)造單元(華北裂陷帶、膠遼地塊和魯淮地塊)的交匯地帶。整體上東臨膠遼地塊的魯東隆起，西南臨魯淮地塊的魯西隆起，跨越濟陽拗陷和沂沭斷裂帶，區(qū)域地質(zhì)背景十分復(fù)雜(圖1)。該區(qū)古近紀(jì)時期裂陷作用強烈發(fā)展，差異運動的幅度大，有厚達(dá)千米以上的內(nèi)陸河湖相沉積。新近紀(jì)和第四紀(jì)時期區(qū)域范圍內(nèi)差異構(gòu)造活動變?nèi)酰蓮娏伊严葑優(yōu)榱丝傮w緩慢沉降的坳陷區(qū)，新近系和第四系呈超覆關(guān)系蓋于古近系及早期的隆起之上。區(qū)內(nèi)主壓應(yīng)力方向為NEE-EW向?，F(xiàn)今構(gòu)造應(yīng)力場的主導(dǎo)方位為北東東到近東西向的壓應(yīng)力，地震活動不強烈，近期地震活動強度不大。

2 工作原理及測線布置

大地電磁測深法(MT)是以天然電磁場信號為場源，通過觀測同一測點不同頻率的時間域電磁場信號，計算出不同深度的視電阻率值，從而來達(dá)到探測地下不同深度地質(zhì)體電性特征、研究地球深部構(gòu)造的目的。探測的深度隨著頻率的變化而變化，通常高頻反映淺部信息，低頻反映深部信息。因此可以利用MT方法來研究深部構(gòu)造，探測地殼和上地幔的電性結(jié)構(gòu)；探測盆地高阻基底起伏和埋深，劃定盆地范圍及其次級構(gòu)造單元；研究地下空間斷裂構(gòu)造展布。

圖2 研究區(qū)斷裂、構(gòu)造單元和L0890測線位置圖Fig.2 The fracture, tectonic unit and the location of the L0890 test line in the study area

野外探測使用了加拿大鳳凰公司出產(chǎn)的電磁勘探設(shè)備V8電法工作站，該設(shè)備的發(fā)射和接收始終采用GPS同步。每道采用24位A/D轉(zhuǎn)換器，并對數(shù)字信號進行處理(DSP)，保持了高的動態(tài)范圍和分辨率，布置了23條測線，測線間距為2 km，測點間距為2 km，在局部進行適當(dāng)加密測點(圖2)，測線方位角為90°。

3 剖面反演解釋

對野外數(shù)據(jù)進行數(shù)據(jù)編輯、曲線圓滑、極化模式識別及靜態(tài)效應(yīng)等處理，并對數(shù)據(jù)進行了反演處理和圖件編制，通過對圖件分析和研究，研究區(qū)內(nèi)共推斷隱伏斷裂9條(F1：鄌郚-葛溝斷裂、F2：沂水-湯頭斷裂、F3：古城-濰河口斷裂、F4：侯鎮(zhèn)北斷裂、F5：挑溝子-沈家營斷裂、F6：柳瞳斷裂、F7：固堤斷裂、F8：神樹坡-東興斷裂、F9：老河口斷裂)，其中北東向斷裂4條，分別為F1、F2、F3、F8和F9，北西向斷裂2條為F5和F7，近東西向斷裂2條為F4、F6(圖2)。在取得野外數(shù)據(jù)后，對每個測點不同頻率的時間域電磁場信號進行處理，計算出不同深度的視電阻率值，利用視電阻率值進行網(wǎng)格化獲得視電阻率斷面圖。

圖3為本次MT勘查L0890剖面綜合剖面圖，圖3(a)為MT測深視電阻率斷面圖，圖3(b)為推斷地質(zhì)成果圖，其橫向為點號，縱向為深度。該剖面西起張家村西，向東經(jīng)三甲王村、大官莊村、高莊村至濰北農(nóng)場四分場東，呈90°方向展布，剖面長為22 km。

圖3 L0890剖面綜合剖面圖Fig.3 The comprehensive profile map of the L0890 test line(a)視電阻率斷面圖；(b)推斷地質(zhì)剖面圖

該剖面主要控制區(qū)內(nèi)鄌郚-葛溝斷裂和沂水-湯頭斷裂的空間展布特征，在69130號測點附近存在一鉆孔，該鉆孔揭露第四系地層厚度為322 m，新近系地層厚為360 m，古近系地層厚為911 m。剖面的68400號測點附近存在一明顯的等值線梯級帶，其西端為低阻異常，東端為高阻異常，為典型的斷裂異常反映，推斷為鄌郚-葛溝斷裂(F1)。該剖面的68850號測點附近標(biāo)高-600 m至標(biāo)高-2 300 m 深度范圍內(nèi)中阻電性層橫向不連續(xù)，存在一相對低阻"U"字型反映，為典型的斷裂異常反映，推斷為沂水-湯頭斷裂(F2)，此次斷裂構(gòu)造解釋位置與前人劃的斷裂位置一致。

從該剖面的垂向電性特征分析亦可劃分兩大電性層，①在67800～68350號測點之間低阻電性層視電阻率值小于10 Ω·m，等值線稀疏寬緩，推斷為第四系和新近系地層的綜合反映，沉積厚度在1 600 m左右，其下伏低阻電性層視電阻率值在10 Ω·m～100 Ω·m之間，結(jié)合區(qū)內(nèi)地質(zhì)資料，推斷該低阻電性層為古近系漸新統(tǒng)和始新統(tǒng)地層，巖性主要為泥巖、砂巖等；②在68350～70000號測點之間低阻電性層受下伏高阻電性層影響，反映不明顯，根據(jù)鉆孔資料，推斷該低阻電性層第四系和新近系地層的綜合反映，沉積厚度700 m左右，其下伏中阻電性層視電阻率值在 300 Ω·m～1 000 Ω·m之間，推斷為古近系砂巖、玄武巖地層，底界面埋深在2 200 m左右，再往下的高阻電性層視電阻率值大于1 000 Ω·m，與上覆中阻電性層之間存在明顯的等值線梯級帶，結(jié)合區(qū)內(nèi)地質(zhì)資料，推斷該高阻電性層為新元古界和新太古界地層的綜合反映。

4 三維反演結(jié)果分析

本次對大地電磁數(shù)據(jù)進行三維反演采用非線性共軛梯度法大地電磁三維反演程序，獲得了較好的效果(圖4和圖5)。圖4為大地電磁三維反演結(jié)果的垂直切片圖，圖4(a)沿南北方向垂直切片，剖面間距為3.8 km，圖4(b)為沿東西方向垂直切片，剖面間距為2 km，為了方便解釋，圖4(b)東西方向比例尺擴大為2倍。由圖4可以看出，各條斷裂對應(yīng)電阻率梯級帶明顯，在三維電阻率切片圖上異常特征反映清晰，大致可以推斷斷裂延伸深度及水平距離，這些主要斷裂基本控制了研究區(qū)內(nèi)沉積地層的分布。

根據(jù)以往該區(qū)內(nèi)研究成果，該區(qū)跨越東營潛斷陷、濰北潛斷陷兩個Ⅳ級構(gòu)造單元。根據(jù)本次MT測深推斷解譯的研究區(qū)內(nèi)地層結(jié)構(gòu)、隱伏斷裂構(gòu)造的空間展布特征以及電場特征的分析，對研究區(qū)內(nèi)構(gòu)造單元進行重新厘定，劃分出研究區(qū)內(nèi)五個Ⅴ級構(gòu)造單元(圖5)。其中A為濰縣凸起、B為候鎮(zhèn)凹陷、C為濰北凹陷、D為濰北凸起、E青東凹陷。

圖4 三維反演結(jié)果垂直切片圖Fig.4 The 3D inversion results of vertical slice map(a)沿南北方向垂直切片圖；(b)沿東西方向垂直切片圖(東西方向比例擴大2倍)

圖5 研究區(qū)構(gòu)造格架及構(gòu)造單元劃分Fig.5 The construction grid and tectonic unit division of the study area

圖6 MT反演水平切片三維顯示圖Fig.6 The MT inversion horizontal slice 3D display diagram(a)從南西往北東；(b)從南東往北西

候鎮(zhèn)凹陷(B)位于研究區(qū)西南部，其北以侯鎮(zhèn)北斷裂為界，東以F1斷裂為界。該凹陷為昌樂縣斷陷的次級單元。該凹陷整體電場特征表現(xiàn)為低阻異常特征，平面形態(tài)呈橢圓狀，向深處延伸大，對應(yīng)在布格重力異常圖是低異常區(qū)。濰北凹陷(C)位于研究區(qū)的中東部，主要受F1斷裂和F6斷裂控制。根據(jù)三維反演結(jié)果水平切片(圖6)可知，濰北凹陷向深處延伸大并且向北傾斜，淺部形態(tài)不規(guī)則，深部形態(tài)近葫蘆狀。

圖7(a)為電阻率值小于20 Ω·m的低阻體三維空間分布隨深度的變化圖，從圖7(a)中可以看出，濰北凹陷呈近乎平行的北北東向條帶狀分布，新近系地層上寬下窄向下延伸，濰北凹陷低阻地層呈由南向北逐漸加深的趨勢。結(jié)合剖面解釋和鉆孔資料，推斷濰北凹陷新近系沉積層覆蓋厚度最深在-400 m至-1 700 m左右。

圖7(b)為電阻率值小于80 Ω·m的低阻體三維空間分布隨深度的變化圖，根據(jù)地質(zhì)資料及該區(qū)物性，推斷三維地質(zhì)體為斷陷盆地，根據(jù)以往資料推斷為中生代白堊系地層，稱之為"潛凹陷"。凹陷盆地形態(tài)與第四系沉積基本一致，向下處延伸最深處達(dá)5 km，代表了古近系基底地層底界面埋深及形態(tài)。

5 結(jié)論

1)通過該研究區(qū)大地電磁測深數(shù)據(jù)的處理和研究分析，獲得了研究區(qū)隱伏斷裂的位置、性質(zhì)及空間展布。研究區(qū)內(nèi)斷裂構(gòu)造發(fā)育，結(jié)合該區(qū)地質(zhì)、鉆探等資料，對該區(qū)隱伏斷裂構(gòu)造位置進行了劃分，劃分出9條隱伏斷裂。

2)采用非線性共軛梯度法對研究區(qū)大地電磁測深數(shù)據(jù)進行三維反演，獲得較好應(yīng)用效果。推斷解譯了研究區(qū)內(nèi)地層結(jié)構(gòu)、隱伏斷裂以及研究區(qū)地層電場特征，對區(qū)內(nèi)兩個Ⅳ級構(gòu)造單元進行劃分，劃分出濰縣凸起、候鎮(zhèn)凹陷、濰北凹陷、濰北凸起和青東凹陷五個Ⅴ級構(gòu)造單元。

圖7 低阻體在不同深度上的空間分布Fig.7 The spatial distribution of low resistance in different depths(a)電阻率值小于20 Ω·m的低阻體；(b) 電阻率值小于80 Ω·m的低阻體

3)通過對三維反演不同方向，不同深度的切片圖的分析，切片圖能直觀地提供斷裂延伸方向及空間展布情況、構(gòu)造單元形態(tài)、不同電性特征的地質(zhì)體等信息，為研究區(qū)地質(zhì)解釋提供地質(zhì)信息。