臨清坳陷西部任縣凸起基底電性結(jié)構(gòu)研究及地?zé)衢_發(fā)建議

胥博文 劉志龍 張敏

摘? 要：河北省任縣地區(qū)地?zé)豳Y源具有儲量豐富、埋深適宜、開發(fā)潛力大等優(yōu)點。本文應(yīng)用大地電磁測深法（MT法）對臨清坳陷西部任縣凸起深部地層結(jié)構(gòu)進(jìn)行探測研究，共布設(shè)了5條MT剖面，獲得了216個測深點。通過對大地電磁測深點的處理及反演，獲得了任縣周邊區(qū)域4km以淺電性結(jié)構(gòu)模型。阻抗張量分解結(jié)果表明任縣凸起淺部一維特征明顯，深部二維特征增強(qiáng)，主構(gòu)造方向以近南北向為主。任縣凸起4km以淺電性結(jié)構(gòu)差異明顯，從淺到深電阻率逐漸升高，具有縱向分層結(jié)構(gòu)，據(jù)此將研究區(qū)地層從上到下劃分為第四系、第三系、石炭系-二疊系、寒武系-奧陶系，并勾畫出各層頂、底面起伏情況。結(jié)合區(qū)域地質(zhì)資料，認(rèn)為任縣地區(qū)地?zé)衢_發(fā)有利區(qū)域位于研究區(qū)南東側(cè)，目的層位以基巖裂隙巖溶熱儲為佳。

關(guān)鍵詞：大地電磁測深;二維反演;地?zé)峥辈?河北省任縣

中圖分類號：P314.4;P631.325? ? 文獻(xiàn)標(biāo)識碼：A? ? 文章編號：1007-1903（2019）02-0036-06

Abstract：Geothermal resources in Renxian area of Hebei Province have the advantages of abundant reserves， suitable burial depth and great potential for development. In this paper， magnetotelluric sounding （MT） method is applied to detect and study the deep stratum structure of Renxian uplift in western Linqing depression. Five MT profiles are laid out and 216 sounding points are obtained. Through processing and inversion of MT points， a shallow electrical structure model of 4 km in the surrounding area of Renxian County is obtained. The results of impedance tensor decomposition show that the shallow part of Renxian uplift has obvious one-dimensional characteristics， while the deep part has enhanced two-dimensional characteristics. The main structural direction is nearly north-south. Renxian uplift has obvious difference in shallow electrical structure with increasing resistivity from shallow to deep， and has a vertical stratified structure. According to this， the strata in the study area are divided into four layers： Quaternary， Neogene， Carboniferous-Permian and Cambrian-Ordovician， and the top and bottom undulations of each layer are delineated. Combining with regional geological data， it is considered that the favorable geothermal development area in Renxian is located on the southeastern side of the study area， and the target horizon is the bedrock fissure karst thermal reservoir.

Keywords： Magnetotelluric sounding; Two-dimensional inversion; Geothermal exploration; Renxian area of Hebei Province

0 引言

地?zé)豳Y源作為一種新型能源，開發(fā)利用過程具有低成本、可持續(xù)利用以及環(huán)保等優(yōu)點，這對于改善能源結(jié)構(gòu)，解決日趨嚴(yán)重的全球環(huán)境問題具有重要意義。臨清坳陷是中-新生代渤海灣盆地的一部分，位于渤海灣盆地西南端，屬于平原沉降區(qū)，總面積18000km2，呈北北東向展布，油氣資源豐富（姜慧超等，2008;呂大煒等，2012）。該區(qū)開展了多項地質(zhì)、地球物理勘查工作（徐守余等，2005;呂希學(xué)，2006;周建勛等，2006），闡明了臨清坳陷的構(gòu)造發(fā)展史、構(gòu)造樣式以及地球物理場特征等，但研究內(nèi)容多集中在臨清坳陷東部地區(qū)，主要用來評價油氣資源，對地?zé)豳Y源涉及較少。目前任縣地區(qū)面臨節(jié)能減排的艱巨任務(wù)，尋找清潔、可持續(xù)的替代能源成為當(dāng)務(wù)之急，任縣地區(qū)地?zé)衢_發(fā)利用尚屬空白區(qū)，因此有必要對任縣所處區(qū)域地?zé)豳Y源賦存情況進(jìn)行摸底排查，以期為地?zé)衢_發(fā)利用提供有效支撐。但是現(xiàn)有的地質(zhì)、地球物理認(rèn)識，并不足以對小范圍的任縣區(qū)域地?zé)豳Y源開發(fā)利用提供有效支撐，如需開發(fā)利用任縣地區(qū)地?zé)豳Y源，需進(jìn)行大比例尺的勘探工作。

大地電磁測深法作為天然源電磁法，在查明基底電性結(jié)構(gòu)特征、劃分?jǐn)嗔褬?gòu)造、評價地?zé)豳Y源等方面具有獨特優(yōu)勢。由于方法理論和計算機(jī)技術(shù)的發(fā)展，二維大地電磁反演已經(jīng)非常成熟（李焱等，2010;陳小斌等，2014;霍光譜等，2015）。本次工作將根據(jù)在任縣周邊采集的216個MT測點為基礎(chǔ)，經(jīng)過系統(tǒng)處理及反演，獲得研究區(qū)域4km以淺導(dǎo)電性結(jié)構(gòu)模型，并進(jìn)行地質(zhì)解釋，為地?zé)豳Y源開發(fā)提供依據(jù)。

1 區(qū)域地質(zhì)背景

任縣凸起位于臨清坳陷西部，是臨清坳陷的次級構(gòu)造單元，北以隆堯斷層為界與隆堯凸起相鄰，西以邢臺-安陽斷裂為界與太行山隆起區(qū)相鄰，東以邢家灣斷裂帶為界與巨鹿凹陷相接，南以臨漳-魏縣斷裂為界與邯鄲凹陷相鄰（圖1）。任縣凸起形成于印支-燕山早期，新近紀(jì)活動強(qiáng)烈，整體下沉與華北平原連為一體。該區(qū)基底以石炭—二疊系為主，靠近山前斷裂帶見少量三疊系殘留，缺失侏羅—白堊系沉積，古近系沉積不發(fā)育，新近系及第四系遍布全區(qū)。任縣凸起地溫分布主要受地質(zhì)構(gòu)造格局的控制，地溫高低與下伏基底構(gòu)造的起伏基本一致，凸起部位地溫梯度大，凹陷部位地溫梯度小，新生界地溫與基巖埋深關(guān)系密切。

研究區(qū)自上而下主要有3個熱儲層，即新近系明化鎮(zhèn)組孔隙型熱儲層、新近系館陶組孔隙型熱儲層以及寒武—奧陶系巖溶裂隙型熱儲層（河北省環(huán)境地質(zhì)勘察院，2003）。明化鎮(zhèn)組以棕黃、淺棕黃、淺棕、棕紅色泥巖、砂質(zhì)泥巖與淺棕色砂巖、含礫砂巖交互沉積為主，巖性具有上粗下細(xì)的特征;館陶組上部為淺棕、黃綠色泥巖與黃棕、淺棕色泥質(zhì)粉砂巖、粉砂巖互層，中部以棕、黃綠色泥巖、粉砂質(zhì)泥巖為主，夾灰白、灰綠色粉砂巖、細(xì)砂巖及含礫砂巖;下部以棕色泥巖、粉砂質(zhì)泥巖與淺棕、灰綠、灰白色泥質(zhì)粉砂巖、細(xì)砂巖互層為主;寒武—奧陶系以灰色、深灰色灰?guī)r、灰質(zhì)白云巖、白云巖為主（河北省地質(zhì)礦產(chǎn)局，1996）。

2 大地電磁數(shù)據(jù)采集、處理與分析

2.1 數(shù)據(jù)采集與處理

此次MT工作共布設(shè)5條MT測線，其中RX_L1及RX_L2測線方位為北西西向，RX_L3及RX_L4測線方位為北北東向，RX_L5為近南北向。測點間距為300m，共完成測點216個，實測點位如圖1。數(shù)據(jù)采集使用加拿大鳳凰公司的MTU-5A系統(tǒng)，配備MTC-80H磁傳感器。同時對電場和磁場進(jìn)行測量，共有6套儀器同時采集，采集頻段0.01～320Hz。每個測點數(shù)據(jù)采集時間不小于3h，利用高精度RTK測量電極距，控制精度為厘米級。依據(jù)DZ/T 0173-1997以及SY/T5820-2014標(biāo)準(zhǔn)對MT數(shù)據(jù)進(jìn)行質(zhì)量評價。在所有頻點都參與評級的情況下，其中一級點176個，占測點總數(shù)81.5%，二級點40個，占測點總數(shù)18.5%，數(shù)據(jù)質(zhì)量符合規(guī)范要求，可以進(jìn)行后續(xù)處理解釋。經(jīng)過處理后典型數(shù)據(jù)如圖2所示。

2.2 維性分析

在進(jìn)行二維反演之前，需對數(shù)據(jù)進(jìn)行阻抗張量分解，獲得一些維性參數(shù)及電性主軸參數(shù)等。MT阻抗張量分解方法多種多樣，主要有Swift分解（Swift，1967）、Bahr分解（Bahr，1991）、GB分解（Groom et al，1989）和相位張量分解（Caldwell et al，2004）等。不同的阻抗張量分解手段得到的參數(shù)不同，但都可以獲得二維偏離度和主軸方位角，二維偏離度是判斷研究區(qū)構(gòu)造二維性的參數(shù)，主軸方位角結(jié)合區(qū)域地質(zhì)情況可以判斷區(qū)域構(gòu)造的走向，通常以玫瑰圖形式給出。

圖3為任縣地區(qū)5條測線0.1～320Hz Swift及Bahr法分解的二維偏離度值。通常認(rèn)為二維偏離度值小于0.3，則區(qū)域構(gòu)造可以近似為二維構(gòu)造。從圖3可見：整體上說兩種阻抗張量分解方法獲得的二維偏離度值在頻率大于1Hz時幾乎都小于0.2，在0.01～0.1Hz時多介于0.2～0.4之間。表明研究區(qū)構(gòu)造較簡單，可以近似看作二維模型從而進(jìn)行二維反演。研究區(qū)在任縣凸起內(nèi)部，淺表沉積較厚的第四系及新近系，局部可能受到內(nèi)丘-任縣斷裂及邢臺-安陽斷裂影響，因此高頻部分（>1Hz）一維特征明顯，低頻部分（<1Hz）二維特征增強(qiáng)。

3 二維反演

根據(jù)阻抗張量分解結(jié)果及區(qū)域構(gòu)造地質(zhì)特征，對5條線分別進(jìn)行了二維反演，反演采用基于非線性共軛梯度算法（Rodi et al，2001）的winglink軟件。反演過程中試算了不同模式、不同參數(shù)的多個模型，最終選取的模型設(shè)置參數(shù)為：100Ω·m的均勻半空間;TM模式;頻率范圍0.1～320Hz;正則化因子tau=5。RX_L1反演擬合差為2.0056，RX_L2反演擬合差為2.1814。RX_L3反演擬合差為2.3510，RX_L4反演擬合差為2.4313，RX_L5反演擬合差為1.7303，反演結(jié)果如圖4。

以RX_L2線為例，該剖面的視電阻率最小值約1Ω·m，最大值約956Ω·m，整個剖面呈現(xiàn)的是高—低—高3層結(jié)構(gòu)。地面以下400m左右為第一電性層，視電阻率在1～956Ω·m之間，此層不穩(wěn)定，中間夾有高阻圈閉。深度400～1700m之間為第二電性層，從整體趨勢看，此層呈較穩(wěn)定的層狀結(jié)構(gòu)，視電阻率在3～36Ω·m之間。深度1700m以下為第三電性層，此層較穩(wěn)定，整體呈層狀展布，視電阻率隨深度穩(wěn)定增加，視電阻率在20～190Ω·m之間。

評價反演結(jié)果是否可靠，除了反演擬合差不能太大外，反演后的擬斷面圖與原始觀測數(shù)據(jù)擬斷面圖的對比也是重要的參考資料，如果兩者相一致，反演模型從至少從數(shù)學(xué)上來說是可靠的。以RX_L2為例，做了二維反演前后視電阻率及頻率擬斷面圖。如圖5，視電阻率及阻抗相位擬斷面圖的色標(biāo)各保持一致。通過對比可以看出，無論是視電阻率還是阻抗相位，反演模型響應(yīng)數(shù)據(jù)均與原始觀測數(shù)據(jù)高度一致，僅個別頻點有差別，這充分說明了本次二維反演結(jié)果可靠。

4 地電結(jié)構(gòu)模型討論及地?zé)衢_發(fā)建議

從反演得到的電性結(jié)構(gòu)模型（圖4）可以看出，研究區(qū)4km以淺電性結(jié)構(gòu)具有明顯的分層性，淺部以低阻為特征，往深處電阻率值逐漸升高，電阻率變化范圍從幾Ω·m到幾百Ω·m。圖5的擬斷面圖也顯示了研究區(qū)電性結(jié)構(gòu)的成層性及上部低阻、下部高阻特征。根據(jù)其他地質(zhì)、地球物理資料并結(jié)合本次二維反演得到的電性結(jié)構(gòu)模型，把研究區(qū)4km以淺地層劃分為4層結(jié)構(gòu)（圖6），從上到下分別是第四系（Q）、第三系（N+E）、石炭系—二疊系（C-P）、寒武系—奧陶系（?+O）。各層劃分依據(jù)及電性特征分述如下：

第四系（Q）巖性主要為亞粘土、亞砂土等，在反演圖上表現(xiàn)為低阻特征，電阻率值小于10Ω·m，底面埋深400m左右，區(qū)域上厚度變化不大。第三系包括新近系（N）和古近系（E）兩部分，主要巖性為泥巖、粉砂巖及細(xì)砂巖等。電阻率值變化較大，反演圖上表現(xiàn)為低阻—高阻的過渡，底板埋深在2000m左右，研究區(qū)內(nèi)該面有一定的起伏，RX_L1線較RX_L2線第三系厚度小，底板埋深小，RX_L4線較RX_L3線第三系厚度小，底板埋深小，因此第三系底面可以看作是東南往西北方向的單傾面。晚古生界包括石炭系（C）和二疊系（P）兩套地層，巖性主要為碎屑巖、粘土巖、炭質(zhì)頁巖夾煤層等。電阻率值較高，一般大于100Ω·m，厚度800m左右，研究區(qū)內(nèi)厚度變化不大，頂面與上覆第三系不整合接觸。早古生界包括寒武系（?）和奧陶系（O）兩套地層，巖性主要為灰?guī)r、白云巖及頁巖等。在反演圖上表現(xiàn)為高阻特征電阻率值較晚古生界地層更高，頂面與上覆晚古生界平行不整合接觸。

新近系熱儲層分為明化鎮(zhèn)組和館陶組上下兩層，主要熱儲層為館陶組。根據(jù)區(qū)域地?zé)岬刭|(zhì)資料分析，該層為本區(qū)熱儲的主要開采層，出水溫度可達(dá)50℃，水量一般大于1000m3/d，水化學(xué)類型屬HCO3·Cl-Na或Cl-Na型水。大地電磁測深推斷結(jié)果顯示，館陶組熱儲層具有分布面積廣、層位穩(wěn)定、連續(xù)型好的特點，適宜進(jìn)行規(guī)模性開發(fā)。

研究區(qū)寒武—奧陶系為厚層灰?guī)r及白云質(zhì)灰?guī)r層，頂面埋深均在2700m以上，厚度一般大于1000m，上覆地層為石炭—二疊系。任縣凸起早古生代經(jīng)歷了兩次沉降和兩次抬升運動，沉降時期形成了以碳酸鹽巖為主的地層，抬升時期則使碳酸鹽巖地層發(fā)生一定剝蝕。一般來說沉積間斷持續(xù)的時間越長，遭受風(fēng)化淋蝕越強(qiáng)，碳酸鹽巖層孔、縫、洞就越發(fā)育。因此，研究區(qū)寒武—奧陶系具備一定形成良好熱儲層的潛力。根據(jù)本次大地電磁測深結(jié)果可知，研究區(qū)東南部寒武—奧陶系埋深較淺，基底高阻特征明顯，電阻率較為穩(wěn)定，因此，建議在下一步地?zé)徙@探過程中，應(yīng)根據(jù)勘探實際需求，選擇合理的熱儲目的層，將會取得較好的效果。

5 結(jié)論

（1）任縣凸起任縣地區(qū)電性結(jié)構(gòu)分層特征明顯，從上到下電阻率值逐漸升高，4km深度范圍地層可劃分4層：從上到下分別是第四系（底面埋深400m左右）、第三系（底面埋深1200～1500m之間）、石炭系—二疊系（底面埋深2700～2900m之間）和寒武系—奧陶系。

（2）任縣地區(qū)地層成層狀分部，第三系底面埋深南東淺、北西深，早古生代寒武系—奧陶系頂面埋深南東淺、北西深。

（3）任縣地區(qū)地?zé)衢_發(fā)有利區(qū)域位于研究區(qū)南東側(cè)，目的層位以基巖裂隙巖溶熱儲為首選。

參考文獻(xiàn)

陳小斌，蔡軍濤，王立鳳，等，2014.大地電磁資料精細(xì)處理和二維反演解釋技術(shù)研究（四）——阻抗張量分解的多測點-多頻點統(tǒng)計成像分析[J].地球物理學(xué)報，57（6）： 1946-1957.

河北省地質(zhì)礦產(chǎn)局，1996. 河北省巖石地層[M]. 武漢：中國地質(zhì)大學(xué)出版社.

河北省環(huán)境地質(zhì)勘察院， 2003.邢臺市地?zé)豳Y源調(diào)查評價報告[R].

霍光譜，胡祥云，黃一凡，等，2015.帶地形的大地電磁各向異性二維模擬及實例對比分析[J].地球物理學(xué)報，58（12）： 4696-4708.

姜慧超，張勇，任鳳樓，等，2008.濟(jì)陽、臨清坳陷及魯西地區(qū)中新生代構(gòu)造演化對比分析[J].中國地質(zhì)，（5）： 963-974.

李焱，胡祥云，吳桂桔，等，2010.基于MPI的二維大地電磁正演的并行計算[J].地震地質(zhì)，32（3）： 392-401.

呂大煒，趙洪剛，李增學(xué)，等，2012.渤海灣盆地臨清坳陷晚古生代古地理特征[J].古地理學(xué)報，14（4）： 437-450.

呂希學(xué)，2006.渤海灣盆地臨清坳陷東部熱演化歷史和成烴史[J].地質(zhì)科學(xué)，（4）： 676-687.

徐守余，嚴(yán)科，2005.渤海灣盆地構(gòu)造體系與油氣分布[J].地質(zhì)力學(xué)學(xué)報，（3）： 259-265.

周建勛，周建生，2006.渤海灣盆地新生代構(gòu)造變形機(jī)制：物理模擬和討論[J].中國科學(xué).D輯：地球科學(xué)，（6）： 507-519.

Bahr K，1991. Geological noise in magnetotelluric data：a classification of distortion types[J].Physics of the Earth and Planetary Interiors，66（1-2）：24-38.

Caldwell T G， Bibby H M， Brown C，2004. The magnetotelluric phase tensor[J]. Geophysical Journal International，158（2）：457-469.

Groom R W， Bailey R C，1989. Decomposition of magnetotelluric impedance tensors in the presence of local three-dimensional galvanic distortion[J]. Journal of Geophysical Research， 94（B2）：1913-1925.

Rodi W，Mackie R L，2001.Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion[J]. Geophysics，（1）：174-187.

Swift C M，1967.A magnetotelluric investigation of an electrical conductivity anomaly in the southwestern United States[D].Cambridge：Massachusetts Institute of Technology.