欒曉東，底青云，安志國(guó)，許誠(chéng)，王顯祥，張文偉

?

基于CSAMT結(jié)果的川東大池干井構(gòu)造淺部電性結(jié)構(gòu)特征分析

欒曉東1, 2, 3，底青云1, 2, 3，安志國(guó)1, 2, 3，許誠(chéng)1, 3，王顯祥1, 3，張文偉1, 2, 3

(1. 中國(guó)科學(xué)院 地質(zhì)與地球物理研究所，北京，100029；2. 中國(guó)科學(xué)院 地球科學(xué)研究院，北京，100029；3. 中國(guó)科學(xué)院大學(xué)，北京，100049)

為揭示川東大池干井構(gòu)造的淺部地質(zhì)構(gòu)造及巖石的電性特征和電磁響應(yīng)規(guī)律，探索電磁法在油氣探測(cè)中的應(yīng)用，將可控源音頻大地電磁法(CSAMT)應(yīng)用于大池干井構(gòu)造油氣區(qū)，采用多種新技術(shù)進(jìn)行CSAMT數(shù)據(jù)處理，獲得地下電性體真實(shí)的電阻率特征，揭示該區(qū)的構(gòu)造具有橫向分段、縱向分層的電性特征。研究結(jié)果表明：東南部高陡背斜核部表現(xiàn)出相對(duì)高阻，此部分侏羅系地層較薄，主要是三疊系和二疊系地層呈層性分布；控制背斜的斷層在電阻率上表現(xiàn)出次低阻的特性；西北翼表現(xiàn)出大范圍低阻特性，主要是相對(duì)較厚的侏羅系上統(tǒng)和中統(tǒng)地層；在中部明顯反映出單斜存在，主要表現(xiàn)為次高阻的侏羅系下統(tǒng)和三疊系上統(tǒng)地層。通過(guò)CSAMT法研究電性差異并結(jié)合地質(zhì)構(gòu)造可達(dá)到探測(cè)油氣藏的目的。本工作可為大池干井構(gòu)造地區(qū)從電性資料識(shí)別油氣構(gòu)造圈閉提供 參考。

川東油氣區(qū)；大池干井構(gòu)造；CSAMT；淺部電性結(jié)構(gòu)

川東地區(qū)是四川盆地內(nèi)大、中型氣田最多的地 區(qū)[1]，由一系列NE?NNE向格擋式褶皺組成[2?3]。大池干井構(gòu)造處于川東褶皺帶的東部，發(fā)育許多不同類(lèi)型的良好圈閉[4?6]。位于大池干井構(gòu)造東北部的石寶寨地區(qū)逆沖斷層發(fā)育，構(gòu)造復(fù)雜，經(jīng)20多年的勘探證實(shí)，本區(qū)及鄰區(qū)高陡背斜大部為含氣構(gòu)造，主要儲(chǔ)氣地層是二疊系長(zhǎng)興組和石炭系黃龍組，其地層構(gòu)造特征屬于典型的大池干井儲(chǔ)氣構(gòu)造。為研究大池干井的淺部地質(zhì)構(gòu)造及巖石的電性特征和電磁響應(yīng)規(guī)律，探索電磁法在油氣探測(cè)中的應(yīng)用，本文作者在石寶寨地區(qū)開(kāi)展可控源音頻大地電磁法(CSAMT)的探測(cè)研究工作，結(jié)合2D反演結(jié)果分析該地區(qū)由淺至深的電阻率變化規(guī)律，以期了解工作區(qū)內(nèi)石炭系儲(chǔ)層上覆巖層的巖性分布、產(chǎn)狀及斷裂構(gòu)造情況的電性特征。

1 研究區(qū)地質(zhì)和地球物理特征

大池干井構(gòu)造帶為長(zhǎng)條形高陡背斜，呈北北東—南南西向展布，延伸達(dá)110 km左右，寬約10余km，東南翼陡，西北翼緩。在淺層，陡緩兩翼通過(guò)陡傾帶過(guò)渡到向斜帶；在中深層(中三疊—中上寒武統(tǒng))，陡緩兩翼對(duì)應(yīng)淺部構(gòu)造的倒轉(zhuǎn)帶和陡傾帶，均形成了傾軸逆斷層。從地面到地下的縱向剖面上，淺層構(gòu)造以寬大隆起為主，地層彎曲倒轉(zhuǎn)；中層(嘉四—石炭系)褶皺斷層發(fā)育，石炭系儲(chǔ)層即位于該層，分布穩(wěn)定，它與下伏直流暗色泥質(zhì)烴源巖及上覆下二疊統(tǒng)底部泥質(zhì)巖構(gòu)成完整的生、儲(chǔ)、蓋組合；深層(志留系以下)構(gòu)造起伏較緩。在大池干井構(gòu)造的主體部位淺層的背斜軸恰好對(duì)應(yīng)石炭—二疊系的主斷凹，尤其在構(gòu)造傾末端(即測(cè)區(qū)所在位置)，垂向變異基本屬于同心等厚褶皺，中層構(gòu)造與淺層構(gòu)造對(duì)應(yīng)良好，構(gòu)造完整，閉合范圍較大[6]，為該地區(qū)通過(guò)研究淺層構(gòu)造揭示石炭系油氣藏提供了可能。

石寶寨測(cè)區(qū)位于大池干井構(gòu)造帶東北部(見(jiàn)圖1[7])，其西北翼緩，地層傾角小于20°，頂部平緩，東南翼陡，地層傾角為35°，向東北端逐漸變得平緩，圈閉線最低海拔為?2 600 m，高點(diǎn)在西南端，海拔為?2 280 m[8]。橫向主要巖性(見(jiàn)表1)：北部為侏羅系蓬萊鎮(zhèn)組(J3p)、沙溪廟組(J2s)和自流井組(J1z)砂巖、砂泥巖，電阻率較低，約為100 Ω·m；中部為涼高山組(J1l)石英砂巖，含少量大安寨(J1dn)灰?guī)r，電阻率相對(duì)較高，為300~500 Ω·m，地面和淺層(下三疊統(tǒng)以上)為向斜或單斜，二疊系及以下存在構(gòu)造高帶；南部為構(gòu)造的核心，主要出露侏羅系、須家河組(T3x)石英砂巖、砂泥巖和大安寨(J1dn)段灰?guī)r、三疊系嘉陵江組(T1j)、雷口坡組(T2l)灰?guī)r，為相對(duì)高阻，可達(dá)1 000 Ω·m左右。南部方斗山構(gòu)造帶地層破碎，斷層發(fā)育。該地區(qū)已有研究者對(duì)一些物性參數(shù)進(jìn)行了研究[8?14]，統(tǒng)計(jì)結(jié)果如表1所示。表1揭示出含油氣的石炭系黃龍組為數(shù)kΩ·m的高阻，其上覆二疊系蓋層為相對(duì)次高阻，整個(gè)地區(qū)在電阻率上表現(xiàn)出呈層性的差異分布，為CSAMT法的應(yīng)用提供了物理基礎(chǔ)。

2 CSAMT探測(cè)有效性分析

1—侏羅系上統(tǒng)；2—侏羅系中統(tǒng)、下統(tǒng)；3—三疊系上統(tǒng)；4—三疊系中統(tǒng)、下統(tǒng)；5—河流；6—鄉(xiāng)鎮(zhèn)；7—測(cè)區(qū)位置。

表1 物性參數(shù)統(tǒng)計(jì)表

(a) 3D正演模型示意圖；(b) 反演結(jié)果

3 CSAMT法工作參數(shù)

測(cè)線布置如圖3所示，沿北西—南東方向共設(shè)計(jì)了6條CSAMT測(cè)線，測(cè)線L1~L6的長(zhǎng)度分別為 15 320，15 440，14 960，15 080，15 360和15 280 m，點(diǎn)距為40 m，物理測(cè)點(diǎn)數(shù)為2 286個(gè)。發(fā)射偶極位于測(cè)區(qū)西南部，偶極距長(zhǎng)為2.2 km，與最近的L1線的收發(fā)距為11 070 m，與最遠(yuǎn)的L6線的收發(fā)距為 17 140 m。發(fā)射頻率為0.125~9 600 Hz，最大發(fā)射電流為30 A。采集設(shè)備為加拿大鳳凰公司V8多功能電法儀器，在多條測(cè)線的多個(gè)測(cè)點(diǎn)用自主研制設(shè)備CLEM電磁法系統(tǒng)進(jìn)行了復(fù)測(cè)對(duì)比。在測(cè)量過(guò)程中，先用V8儀器進(jìn)行第1輪采集，采集完成后緊接著用CLEM儀器進(jìn)行第2輪采集，為得到更為準(zhǔn)確的對(duì)比結(jié)果，CLEM儀器與V8儀器共用1套接收線、不極化電極及磁傳感器，具體采集方式如圖4所示。

圖3 測(cè)線示意圖

圖4 V8儀器與CLEM儀器采集方式示意圖

4 CSAMT數(shù)據(jù)處理及結(jié)果分析

4.1 CSAMT的數(shù)據(jù)處理

對(duì)CSAMT數(shù)據(jù)的處理，包括野外數(shù)據(jù)預(yù)處理和資料后期處理2個(gè)部分。預(yù)處理包括干擾校正、近場(chǎng)校正、靜態(tài)校正；后期處理主要是多次的反演擬合。

4.1.1 干擾校正

由于測(cè)線穿過(guò)村莊、道路及民用線，外界干擾是非穩(wěn)態(tài)的。為壓制這種非穩(wěn)態(tài)干擾，對(duì)干擾點(diǎn)采用希爾伯特?黃變換(HHT)[15?17]法進(jìn)行干擾校正。HHT是一種針對(duì)非線性、非穩(wěn)態(tài)信號(hào)自適應(yīng)處理的方法，廣泛應(yīng)用于語(yǔ)音、地震、水聲等信號(hào)處理[18?20]。

在這一理論中，引入固有模函數(shù)(IMF)的概念，每個(gè)模函數(shù)都有其控制頻率。通過(guò)經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)，可將原始信號(hào)自適應(yīng)分解為有限個(gè)IMF和1個(gè)表征信號(hào)趨勢(shì)變化的殘余信號(hào)。對(duì)每個(gè)IMF進(jìn)行HHT變換，得到瞬時(shí)頻率，從而將信號(hào)表示成頻率?時(shí)間?能量的分布，稱(chēng)為Hilbert譜。為壓制干擾噪聲，一方面，剔除受噪聲干擾較強(qiáng)的IMF；另一方面，在Hilbert譜上選擇干擾噪聲最小、有用信號(hào)最強(qiáng)的時(shí)間窗進(jìn)行信號(hào)分析，最大限度地降低噪聲的干擾。

以L4線7 020 m點(diǎn)的電場(chǎng)信號(hào)為例，此段信號(hào)的發(fā)射頻率為32 Hz。該點(diǎn)位于110 kV高壓線附近，旁邊有幾戶(hù)農(nóng)家的村莊(圖3中L4線過(guò)高壓線位置)，原始信號(hào)受50 Hz干擾嚴(yán)重，能量主要集中在50 Hz附近，通過(guò)HHT變換，得到Hilbert譜(見(jiàn)圖5)，此處用了前12個(gè)IMF。從圖5可見(jiàn)：第1個(gè)IMF主要受50 Hz干擾控制，第2個(gè)IMF則顯示出所發(fā)射的32 Hz的有用信號(hào)，其他IMF則受更低頻率干擾源的影響。

以上五個(gè)方面的問(wèn)題在一定程度上影響了學(xué)生學(xué)習(xí)數(shù)據(jù)庫(kù)課程的積極性，為改善這一現(xiàn)象，筆者經(jīng)過(guò)長(zhǎng)時(shí)間的課程探索和實(shí)踐，探索出了以下解決方案。

圖5 Hilbert譜

為剔除干擾，僅選用第2個(gè)IMF作為有用信號(hào)，重建信號(hào)，如圖6所示。由圖6可見(jiàn)：經(jīng)過(guò)HHT變換后，在頻譜上50 Hz的信號(hào)被完全壓制，32 Hz信號(hào)的能量比其他頻率信號(hào)的能量高。對(duì)磁場(chǎng)亦進(jìn)行類(lèi)似處理，最后重新計(jì)算視電阻率。經(jīng)過(guò)干擾校正后，視電阻率曲線的連續(xù)性比原始曲線的好(見(jiàn)圖7)，基本能反映地下地電結(jié)構(gòu)的變化，說(shuō)明干擾噪聲得到有效壓制。

(a)時(shí)域；(b)頻域

1—校正前；2—校正后。

4.1.2 近場(chǎng)校正

CSAMT定義的視電阻率為卡尼亞視電阻率，在遠(yuǎn)區(qū)，能夠客觀地反映電斷面的垂向變化。而在過(guò)渡區(qū)和近區(qū)，由于不滿(mǎn)足場(chǎng)垂直入射的原理，卡尼亞視電阻率將發(fā)生嚴(yán)重的畸變，不能正確地反映電斷面特征，這種現(xiàn)象稱(chēng)為近場(chǎng)效應(yīng)[21?22]。因此，為得到合理的反演結(jié)果，必須對(duì)近場(chǎng)效應(yīng)進(jìn)行校正。傳統(tǒng)的“倒三角形”法[23]與場(chǎng)值校正法[24]仍有一定誤差，且場(chǎng)值定義的視電阻率與卡尼亞視電阻率有差別；湯井田 等[25?26]提出采用場(chǎng)值定義新的全區(qū)視電阻率以解決此問(wèn)題，由于場(chǎng)值受干擾影響大，會(huì)使定義的視電阻率受到較大干擾?；诶铤Q等[27]的工作，本文作者提出基于牛頓迭代法的視電阻率定義法來(lái)校正CSAMT近場(chǎng)效應(yīng)。

在均勻半空間，水平電偶源激勵(lì)產(chǎn)生的電磁場(chǎng)分量為：

迭代收斂條件為前、后2次的校正電阻率誤差小于預(yù)定數(shù)值或達(dá)到一定迭代次數(shù)。相比傳統(tǒng)方法，此方法能夠準(zhǔn)確收斂到真實(shí)值。圖8所示為L(zhǎng)1線1 000 m點(diǎn)實(shí)測(cè)數(shù)據(jù)近場(chǎng)校正前、后對(duì)比。由圖8可見(jiàn)：在校正前，視電阻率從2 Hz開(kāi)始明顯進(jìn)入近場(chǎng)，呈45°直線上升，當(dāng)視電阻率為0.125 Hz時(shí)，達(dá)1 000 Ω·m，產(chǎn)生假高阻異常；經(jīng)牛頓迭代法改正后，高頻數(shù)據(jù)與原始數(shù)據(jù)一致，低頻數(shù)據(jù)平緩收斂到100 Ω·m。

4.1.3 靜態(tài)校正

1—校正前；2—校正后。

(a) 校正前；(b) 校正后

在校正前(見(jiàn)圖9(a))，數(shù)據(jù)受靜態(tài)影響比較大，擬斷面圖上有多條縱向“直線”；經(jīng)校正后(見(jiàn)圖9(b))，靜態(tài)效應(yīng)被剔除，而且很好地保留了高阻異常。最后對(duì)校正后的數(shù)據(jù)根據(jù)式(5)還原幅值。

4.1.4 反演擬合分析

對(duì)數(shù)據(jù)進(jìn)行各項(xiàng)校正后，再進(jìn)行2D圓滑反演，用計(jì)算得到的反演結(jié)果和原始曲線進(jìn)行擬合分析，根據(jù)擬合差判斷前期處理結(jié)果的合理性。若擬合誤差較大，則說(shuō)明數(shù)據(jù)處理存在問(wèn)題，要重新進(jìn)行校正與處理，如此反復(fù)多次，直到反演結(jié)果與原始曲線擬合差達(dá)到反演結(jié)果合理要求為止。圖10所示為反演處理結(jié)果的對(duì)比。由圖10可見(jiàn)：在0 m附近、4 500~7 000 m、8 000~12 000 m以及地表出現(xiàn)不正常的假高阻異常；經(jīng)過(guò)1次處理后，0 m附近高阻異常消失，4 500~ 7 000 m的異?；貧w正常，但仍殘留部分地表及9 000~ 12 000 m的假高阻異常；再進(jìn)一步進(jìn)行反復(fù)多次處理后，地表上kΩ·m的假高阻異常消失，相比于原始數(shù)據(jù)，8 000~12 000 m的高阻異?；貧w正常形態(tài)。

4.2 淺層地電結(jié)構(gòu)分析

采用上述2D圓滑反演方法對(duì)6條測(cè)線進(jìn)行了2D圓滑反演，在反演剖面上進(jìn)行地質(zhì)分層及成果推斷解釋。圖11所示為L(zhǎng)2線的2D反演及推斷解釋斷面圖。L2測(cè)線剖面上明顯分為5個(gè)部分。

距離0~2 000 m為第1部分，此區(qū)域呈現(xiàn)出大面積的低阻，電阻率低于150 Ω·m，分為上、下層：第1層厚度約為400 m，電阻率為100 Ω·m左右；第2層一直延伸到底部，電阻率在50 Ω·m以下。從地質(zhì)資料可以推斷：此區(qū)域?yàn)橘_系砂巖，第1層可推斷為侏羅系上統(tǒng)，主要是遂寧組和地表的蓬萊鎮(zhèn)組；第2層推斷為侏羅系中統(tǒng)沙溪廟組，在北端厚度可達(dá)1 000 m以上，向東南方向逐漸變薄，與該地區(qū)西北翼平緩的地質(zhì)結(jié)構(gòu)相吻合。

(a) 原始數(shù)據(jù)反演結(jié)果；(b) 一次處理后的反演結(jié)果；(c)多次處理后的反演結(jié)果

圖11 L2線的2D反演及推斷地質(zhì)解釋斷面圖

距離2 000~5 600 m為第2部分，此部分為明顯的單斜。電阻率剖面上顯示出2段次高阻：第1段電阻率為400~500 Ω·m，總厚度約為700 m，根據(jù)地質(zhì)資料，推斷為侏羅系下統(tǒng)，主要是涼高山組的砂巖和自流井組大安寨段的部分灰?guī)r。由于夾雜灰?guī)r，電阻率相對(duì)下段稍高；緊接著第2段次高阻，電阻率約為300 Ω·m，厚度約為400 m，可推斷為三疊系上統(tǒng)須家河組。

距離5 600~8 800 m為第3部分，根據(jù)已知的地質(zhì)資料，此區(qū)域的低阻反映的是控制背斜的大斷層。

距離8 800~13 600 m為第4部分，該區(qū)域?yàn)楦叨笜?gòu)造核部，有明顯的成層性分布。為對(duì)該部分進(jìn)行準(zhǔn)確分層，根據(jù)本文反演出的電性結(jié)構(gòu)差異并結(jié)合表1的物性參數(shù)、當(dāng)?shù)氐牡刭|(zhì)結(jié)構(gòu)，可分析出該部分的地質(zhì)結(jié)構(gòu)如下：從地表到海拔200 m段，電阻率約為 100 Ω·m，顯示為低阻，推斷為侏羅系砂巖地層；海拔高程為?300~200 m，電阻率為300~500 Ω·m，顯示為次高阻，推斷為三疊系上統(tǒng)須家河組石英砂巖。下伏地層僅顯示出明顯高阻，電阻率在900 Ω·m以上，單從電阻率斷面上難以劃分。但由表1可知，嘉陵江組的電阻率稍高于上覆的雷口坡組和下伏的飛仙關(guān)組的，可進(jìn)行如下解釋?zhuān)汉０胃叱虨?500~?300 m，電阻率為1 000~1 400 Ω·m，推斷為三疊系雷口坡組；海拔高程為?1 300~?500 m，電阻率為1 400~1 800 Ω·m，推斷為三疊系嘉陵江組；海拔高程為?1 700~?1 300 m，電阻率稍降低，約為900 Ω·m，推斷為三疊系飛仙關(guān)組；海拔高程為?2 200~?1 700 m，推斷為二疊系；在最底部?2 300~?2 200 m，電阻率升高，為1 500 Ω·m以上，推斷為石炭系黃龍組的頂層。

距離13 600 m至末端為第5部分。此部分呈現(xiàn)從高阻逐漸向低阻過(guò)渡的特征，根據(jù)圖3中測(cè)線位置可看出，該處測(cè)線的地質(zhì)結(jié)構(gòu)從侏羅紀(jì)過(guò)渡到三疊紀(jì)，結(jié)合實(shí)際地質(zhì)資料可知在背斜南部斷層發(fā)育，故推斷此處為控制背斜核部南側(cè)的大斷層，該斷層與第3部分的斷層相呼應(yīng)，控制了整個(gè)高陡核部。由于斷層斷距較大，測(cè)線未能跨過(guò)此斷層，因此，只隱約表現(xiàn)出漸變過(guò)渡的特征。

對(duì)其他測(cè)線都進(jìn)行了相似的解釋工作，形成了1組垂直切片圖，如圖12所示。從圖12可以看出：在橫向剖面上，6條測(cè)線反演剖面表現(xiàn)出一致的特性。在北端表現(xiàn)為侏羅系的低阻地層；中段為次高阻單斜的特征；南段為構(gòu)造核部，顯示出高阻特性。在西南部的L1與L2線隱約可見(jiàn)石炭系黃龍組引起的高阻異常，但到東北部已無(wú)法看到，說(shuō)明石炭系地層由西南向東北逐漸加深。另外南部的高陡背斜區(qū)，電阻率由西南向東北逐漸趨于平緩，變化幅度較小，說(shuō)明此部分地層向東北端趨于平緩。以上所揭示出的地電結(jié)構(gòu)特征均與前述地質(zhì)結(jié)構(gòu)一致。此次CSAMT法僅揭示了淺部的電性構(gòu)造，對(duì)于更深層的構(gòu)造，可開(kāi)展大地電磁法(MT)等勘探深度較大的電磁方法。

圖12 CSAMT反演垂直切片圖

5 結(jié)論

1) 對(duì)干擾噪聲，采用HHT變換的方法，通過(guò)經(jīng)驗(yàn)?zāi)B(tài)分解，壓制了非穩(wěn)態(tài)的干擾信號(hào)，尤其對(duì)50 Hz工頻干擾的壓制效果明顯。

2) 提出了基于牛頓迭代法的近場(chǎng)校正方法，避免了傳統(tǒng)方法的校正不足或過(guò)校正，使近場(chǎng)校正到合理范圍。

3) 提出了基于Fourier變換的靜態(tài)校正方法，能夠有效壓制靜態(tài)效應(yīng)，凸顯深部異常特征。

4) CSAMT反演結(jié)果與已知的地質(zhì)特征相吻合，揭示試驗(yàn)區(qū)的大池干井構(gòu)造具有橫向分段、縱向分層的電性特征。南部高陡背斜核部表現(xiàn)出相對(duì)高阻，侏羅系地層較薄，主要是三疊系和二疊系地層呈層性分布；控制背斜的斷層在電阻率上表現(xiàn)出次低阻的特性；西北翼表現(xiàn)出低阻特性，主要是相對(duì)較厚的侏羅系上統(tǒng)和中統(tǒng)地層；在中部明顯反映出單斜的存在，主要表現(xiàn)為次高阻的侏羅系下統(tǒng)和三疊系上統(tǒng)的地層。

5) 對(duì)于CSAMT等電磁類(lèi)方法，直接識(shí)別儲(chǔ)層可能比較困難，可采用地質(zhì)結(jié)構(gòu)結(jié)合電性差異的方法來(lái)間接達(dá)到推斷地層，尋找儲(chǔ)層的目的。

6) 本次CSAMT結(jié)果獲得了大池干井石寶寨地區(qū)淺部電性結(jié)構(gòu)以及地質(zhì)結(jié)構(gòu)，為此區(qū)域從CSAMT資料識(shí)別與油氣有關(guān)的油氣構(gòu)造圈閉提供了可借鑒的結(jié)果。

致謝：

感謝國(guó)家重大科研裝備研制項(xiàng)目(ZDYZ2012-1-05)和國(guó)家重大科技基礎(chǔ)設(shè)施建設(shè)項(xiàng)目(WEM)工程地下資源探測(cè)分系統(tǒng)的資助，以及有關(guān)人員對(duì)本文做出的貢獻(xiàn)。

[1] 朱光有, 張水昌, 梁英波, 等. 四川盆地天然氣特征及氣源[J]. 地學(xué)前緣(中國(guó)地質(zhì)大學(xué)(北京); 北京大學(xué)), 2006, 13(3): 234?248. ZHU Guangyou, ZHANG Shuichang, LIANG Yingbo, et al. The characteristics of natural gas in Sichuan basin and its sources[J]. Earth Science Frontiers(China University of Geosciences, Beijing; Peking University), 2006, 13(3): 234?248.

[2] 胡召奇, 朱光, 劉國(guó)生. 等. 川東“侏羅山”式褶皺帶形成時(shí)代: 不整合面的證據(jù)[J]. 地質(zhì)評(píng)論, 2009, 55(1): 32?42. HU Shaoqi, ZHU Guang, LIU Guosheng, et al. The folding time of the eastern Sichuan Jura-type fold belt: evidence from unconformity[J]. Geological Review, 2009, 55(1): 32?42.

[3] 解國(guó)愛(ài), 賈東, 張慶龍, 等. 川東朱羅山式褶皺構(gòu)造帶的物理模擬研究[J]. 地質(zhì)學(xué)報(bào), 2013, 87(6): 773?788. XIE Guoai, JIA Dong, ZHANG Qinglong, et al. Physical modeling of the Jura-type folds in eastern Sichuan[J]. Acta Geologica Scienca, 2013, 87(6): 773?788.

[4] 蒲家奇, 孟德勤, 周文. 川東大池干井構(gòu)造帶形成機(jī)理分析及挖潛勘探評(píng)價(jià)[J]. 天然氣工業(yè), 1995, 15(2): 15?19. PU Jiaqi, MENG Deqin, ZHOU Wen. The analysis of formation mechanism of Dachiganjing structural belt and the evaluation of tapping exploration potentialities in east Sichuan[J]. Natural Gas Industry, 1995, 15(2): 15?19.

[5] 陳宗清. 川東大池干井構(gòu)造帶石炭系天然氣藏與再探目標(biāo)評(píng)價(jià)[J]. 中國(guó)海上油氣(地質(zhì)), 1999, 13(4): 276?281. CHEN Zongqing. Re-evaluation of carboniferous gas pools on Dachiganjing structural zone of east Sichuan[J]. China Offshore Oil and Gas (Geology), 1999, 13(4): 276?281.

[6] 胡光燦, 謝姚祥. 中國(guó)四川東部高陡構(gòu)造石炭系氣田[M]. 北京: 石油工業(yè)出版社, 1997: 30?40. HU Guangcan, XIE Yaoxiang. Carboniferous gas fields in high steep structures of eastern Sichuan[M]. Beijing: Oil Industrial Press, 1997: 30?40.

[7] 中國(guó)地質(zhì)科學(xué)院地質(zhì)研究所.中國(guó)西部及鄰區(qū)地質(zhì)圖(1:250萬(wàn))[M]. 北京: 地質(zhì)出版社, 2006: 1. The Institute of Geology, Chinese Academy of Geological Sciences. Western China and its adjacent area geological map (1:2 500 000)[M]. Beijing: Geological Publishing House, 2006: 1.

[8] 尹雙江. 萬(wàn)順場(chǎng)石炭系氣藏改建地下儲(chǔ)氣庫(kù)可行性研究[D]. 成都: 西南石油大學(xué)石油與天然氣工程學(xué)院, 2012: 13?19. YIN Shuangjiang. The feasibility study of Wanshun gas field in carboniferous gas reservoir converted to underground gas storage[D]. Chengdu:Southwest Petroleum University. College of Petroleum Engineering, 2012: 13?19.

[9] 江為為, 劉伊克, 郝天珧, 等. 四川盆地綜合地質(zhì)、地球物理研究[J]. 地球物理學(xué)進(jìn)展, 2001, 16(1): 11?23. JIANG Weiwei, LIU Yike, HAO Tianyao, et al. Comprehensive study of geology and geophysics of Sichuan Basin[J]. Progress in Geophysics, 2001, 16(1): 11?23.

[10] 馬麗芳. 中國(guó)地質(zhì)圖集[M]. 北京: 地質(zhì)出版社, 2002: 277?284. MA Lifang. Geological atlas of China[M]. Beijing: Geological Publishing House, 2002: 277?284.

[11] 李本亮, 孫巖, 陳偉. 川東層滑系統(tǒng)及其油氣地質(zhì)意義[J]. 石油與天然氣地質(zhì), 1998, 19(3): 244?247.LI Benliang, SUN Yan, CHEN Wei. Layer-slip systems in eastern Sichuan and the significance for petroleum geology[J]. Oil & Gas Geology, 1998, 19(3): 244?247.

[12] 司馬立強(qiáng), 張樹(shù)東, 劉海洲, 等. 川東高陡構(gòu)造陡翼主要構(gòu)造特征及測(cè)井解釋[J]. 天然氣工業(yè), 1996, 16(4): 25?28. SIMA Liqiang, ZHANG Shudong, LIU Haizhou, et al. The main structural characteristics of high and steep structures in East Sichuan and logging interpretation of steep wing[J]. Natural Gas Industry, 1996, 16(4): 25?28.

[13] 吉人, 吳勝, 姜艷玲, 等. 川東地區(qū)高陡復(fù)雜構(gòu)造LWD跟蹤有效儲(chǔ)層的方法[J]. 天然氣工業(yè), 2014, 32(2): 1?6. JI Ren, WU Sheng, JIANG Yanling, et al. LWD technology of effectively tracing effective reservoirs in the high-steep complicated structrues in eastern Sichuan basin[J]. Natural Gas Industry, 2014, 32(2): 1?6.

[14] 雍杰. 川東大池干井構(gòu)造老井復(fù)查及勘探潛力研究[D]. 成都: 西南石油大學(xué)地球科學(xué)與技術(shù)學(xué)院, 2014: 5?41. YONG Jie. Old well review and study of exploration potential of Dachiganjing structure in east Sichuan[D]. Chengdu: Southwest Petroleum University. School of Geoscience and Technology, 2014: 5?41.

[15] HUANG N E, SHEN Z, LONG S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[C]// Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. London, Britain: The Royal Society, 1998: 903?995.

[16] HUANG N E, WU Zhaohua.A review on Hilbert?Huang transform: method and its applications to geophysical studies[J]. Reviews of Geophysics, 2008, 46(2): RG2006-1?23.

[17] HUANG N E, SHEN S S P.Interdisciplinary mathematical sciences: Hilbert–Huang transform and its applications[M]. Hackensack: World Scientific, 2005: 107?127.

[18] ZHAO Zhidong, ZHAO Zhijin, CHEN Yuquan. Time-frequency analysis of heart sound based on HHT (Hilbert–Huang Transform)[C]// Proceedings of International Conference on Communications, Circuits and Systems. Hong Kong, China: IEEE, 2005: 926?929.

[19] WANG Yuqing, PENG Zhenming, He Yanmin. Instantaneous attributes analysis of seismic signals using improved HHT[J]. Journal of Earth Science, 2015, 26(4): 515?521.

[20] ZENG Xiangyang, WANG Shuguang. Underwater sound classification based on Gammatone filter bank and Hilbert- Huang transform[C]// Signal Processing, Communications and Computing (ICSPCC) IEEE International Conference. Beijing, China: IEEE, 2014: 707?710.

[21] 李金銘. 地電場(chǎng)與電法勘探[M]. 北京: 地質(zhì)出版社, 2005: 429?432. LI Jinming. Geoelectric field and electrical exploration[M]. Beijing: Geological Publishing House, 2005: 429?432.

[22] 湯井田, 何繼善. 可控源音頻大地電磁法及其應(yīng)用[M].長(zhǎng)沙: 中南大學(xué)出版社, 2005: 119?125. TANG Jingtian, HE Jishan. Controlled source audio-frequency magnetotellurics method its application[M]. Changsha: Central South University Press, 2005: 119?125.

[23] 周茂軍, 周玉冰. 可控源音頻大地電磁法(CSAMT)的近場(chǎng)效應(yīng)和近場(chǎng)校正[J]. 遼寧地質(zhì), 1993(3): 272?281. ZHOU Maojun, ZHOU Yubing. Near field effect and correction of control source audio-frequency magnetotelluric method[J]. Liaoning Geology, 1993(3): 272?281.

[24] 詹少全, 錢(qián)美平, 馮戔戔. CSAMT 全區(qū)視電阻率電場(chǎng)正演迭代擬合近場(chǎng)校正方法[J]. 物探與化探, 2011, 35(5): 663?665. ZHAN Shaoquan, QIAN Meiping, FENG Jianjian. The CSAMT near-field correction method for electric field apparent resistivity forward iterative fitting[J]. Geophysical & Geochemical Exploration, 2011, 35(5): 663?665.

[25] 湯井田, 何繼善. 水平電偶源頻率測(cè)深中全區(qū)視電阻率定義的新方法[J]. 地球物理學(xué)報(bào), 1994, 37(7): 543?552. TANG Jingtian, HE Jishan. A new method to define the full-zone resistivity in horizontal electric dipole frequency sounding on a layered earth[J]. Chinese Journal of Geophysics, 1994, 37(7): 543?552.

[26] 馮兵, 王珺璐, 周祥文, 等. CSAMT 探測(cè)中電場(chǎng)E全區(qū)視電阻率定義及應(yīng)用[J]. 煤田地質(zhì)與勘探, 2013, 41(6): 78?82. FENG bing, WANG Junlu, ZHOU Xiangwen, et al. Application of full-region apparent resistivity of CSAMTEin exploration[J]. Coal Geology & Exploration, 2013, 41(6): 78?82.

[27] 李鶴, 李桐林, 伍亮. CSAMT 全區(qū)視電阻率轉(zhuǎn)換及其效果分析[J]. 地球物理學(xué)進(jìn)展, 2015, 30(2): 889?893. LI He, LI Tonglin, WU liang. Transformation of all-time apparent resistivity of CSAMT and analysis of its effect[J]. Progress in Geophysics, 2015, 30(2): 889?893.

[28] TORRES?VERDIN C, BOSTICK J F X. Principles of spatial surface electric field filtering in magnetotellurics: electromagnetic array profiling (EMAP)[J]. Geophysics, 1992, 57(4): 603?622.

[29] 黃兆輝, 底青云, 侯勝利. CSAMT的靜態(tài)效應(yīng)校正及應(yīng)用[J]. 地球物理學(xué)進(jìn)展, 2006, 21(4): 1290?1295. HUANG Zhaohui, DI Qing yun, HOU Shengli. CSAMT static correction and its application[J]. Progress in Geophysics, 2006, 21(4): 1290?1295.

[30] HAMDI H, QAUSAR A M, SRIGUTOMO W. CSAMT data processing with source effect and static corrections, application of Occam’s inversion, and its application in geothermal system[C]//6th Asian Physics Symposium, Journal of Physics: Conference Series. Bandung, Indonesia: IOP Publishing, 2016: 012057.

(編輯 劉錦偉)

Shallow electric structure character analysis of Dachiganjing structure in east Sichuan based on CSAMT results

LUAN Xiaodong1, 2, 3, DI Qingyun1, 2, 3, AN Zhiguo1, 2, 3, XU Cheng1, 3, WANG Xianxiang1, 3, ZHANG Wenwei1, 2, 3

(1. Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;2. Institutions of Earth Science, Chines Academy of Sciences, Beijing 100029, China;3. University of Chinese Academy of Sciences, Beijing 100049, China)

In order to study the shallow geoelectric structure of Dachiganjing structure in east Sichuan and the electrical characteristics and electromagnetic response laws of rocks, as well as explore the application of electromagnetism in oil and gas detection, the controlled source audio-frequency magneto-telluric (CSAMT) survey over the Dachiganjing area was carried out and resistivity characteristics of underground electrical bodies were obtained by processing CSAMT data using variety of new techniques. The transverse section of study area showed bed-delineating electric characteristics and piecewise distribution laterally were revealed. The results show that the resistivity of steeply dipping anticline in southeastern part of the study area is relatively high, which means that the Jurassic stratum in this place is relatively thin and the main stratums are Triassic and Dyas. The fault controlled anticline is characterized as second-rate low-resistivity layer. The upper-middle Jurassic series, which mainly locate in the northwestern part of the study area, present as a large low-resistivity zone. There exists an obvious homocline in the middle part of the study area, showing a second-rated high resistivity and it can be inferred from the lower Jurassic and upper Triassic stratum. The potential region of oil and gas accumulation can be detected from subsurface electrical characteristics by combining a priori information of geological structure with CSAMT survey results. This work provides a way to recognize the structural traps, which is related to gas and oil accumulation through the resistivity structure of Dachiganjing area.

oil and gas area in east Sichuan; Dachiganjing structure; CSAMT; shallow electric structure

10.11817/j.issn.1672-7207.2018.05.019

P139.1+2

A

1672?7207(2018)05?1169?10

2017?05?13；

2017?07?03

國(guó)家重大科研裝備研制項(xiàng)目(ZDYZ2012-1-05) (Project(ZDYZ2012-1-05) supported by the National Major Scientific Equipment Research)

底青云，博士，研究員，從事地球電磁數(shù)據(jù)正反演、電磁儀器及隨鉆測(cè)井儀器裝備技術(shù)研發(fā)；E-mail: qydi@mail.iggcas.ac.cn