張春賀,尹成明,甘貴元,王財(cái)富,孫衛(wèi)斌,劉云祥,胡祖志

(1.中國(guó)地質(zhì)調(diào)查局油氣資源調(diào)查中心,北京100029;2.中國(guó)地質(zhì)調(diào)查局油氣資源調(diào)查中心非常規(guī)油氣地質(zhì)重點(diǎn)實(shí)驗(yàn)室,北京100029;3.中國(guó)石油青海油田公司勘探開(kāi)發(fā)研究院,甘肅敦煌736202;4.中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司,河北涿州072751)

重磁電一體化勘探技術(shù)在柴達(dá)木盆地西部復(fù)雜地區(qū)的應(yīng)用

張春賀1,2,尹成明1,2,甘貴元3,王財(cái)富4,孫衛(wèi)斌4,劉云祥4,胡祖志4

(1.中國(guó)地質(zhì)調(diào)查局油氣資源調(diào)查中心,北京100029;2.中國(guó)地質(zhì)調(diào)查局油氣資源調(diào)查中心非常規(guī)油氣地質(zhì)重點(diǎn)實(shí)驗(yàn)室,北京100029;3.中國(guó)石油青海油田公司勘探開(kāi)發(fā)研究院,甘肅敦煌736202;4.中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司,河北涿州072751)

以柴達(dá)木盆地西部為代表的復(fù)雜山地地區(qū),是我國(guó)重要的油氣資源勘探遠(yuǎn)景區(qū)。由于受復(fù)雜地表?xiàng)l件、復(fù)雜地質(zhì)構(gòu)造雙重因素的影響,油氣勘探難度大。首次采用三維重磁電綜合地球物理勘探方法在柴達(dá)木盆地西部復(fù)雜地區(qū)開(kāi)展油氣勘探攻關(guān)。三維大地電磁資料采集以面元為單位代替?zhèn)鹘y(tǒng)排列式采集,重力資料采集利用正交高密度基點(diǎn)網(wǎng)法提高采集精度。針對(duì)研究區(qū)的三維重磁電資料,采用三維預(yù)處理技術(shù)壓制地形及局部不均勻體的影響,利用鉆井資料作為約束,進(jìn)行三維反演處理,獲得了研究區(qū)地層電阻率、密度的空間分布。結(jié)合鉆井、地震等資料對(duì)研究區(qū)開(kāi)展了綜合石油地質(zhì)解釋評(píng)價(jià),提出了有利的勘探目標(biāo)。

油氣勘探;復(fù)雜山地;柴達(dá)木盆地;重力勘探;磁法勘探;大地電磁勘探;反演

地震勘探一直是油氣地球物理勘探中的主導(dǎo)技術(shù)。隨著油氣勘探程度的不斷深入,勘探難度越來(lái)越大,對(duì)地震勘探的效果提出了更加嚴(yán)格的要求[1-2]。重、磁、電地球物理方法雖然在精度和分辨率兩方面都不及地震勘探,但各種方法都有其自身的特點(diǎn),重磁方法橫向分辨率較高;電磁方法,特別是大地電磁方法勘探深度大,同時(shí)具有不受高阻屏蔽以及對(duì)低阻反映靈敏等優(yōu)勢(shì)[3-6]。

近20年來(lái),綜合地球物理技術(shù)在陸域山前構(gòu)造、潛山構(gòu)造、深層火山巖、黃土塬砂體巖性油藏檢測(cè)以及海洋油氣勘探等方面都取得了顯著勘探效果[7-10]。隨著重磁電勘探技術(shù)的進(jìn)步,綜合地球物理技術(shù)在油氣勘探中正逐步由區(qū)域勘探走向目標(biāo)勘探,不斷提高識(shí)別不同地質(zhì)體的能力,降低地球物理方法自身的多解性,為進(jìn)一步深入認(rèn)識(shí)不同勘探目標(biāo)拓展了更加廣闊的空間。

位于柴達(dá)木盆地西部英雄嶺褶皺沖斷帶西南緣的花土溝—獅子溝地區(qū),地表地形起伏劇烈,溝壑縱橫,在400km2的范圍內(nèi),海拔在2800～3700m,相對(duì)高差達(dá)900m。由于受地表施工條件困難、主體構(gòu)造復(fù)雜、地震資料信噪比極低等多重因素的制約,該地區(qū)曾一度被視為地球物理勘探禁區(qū)[11-12]。自20世紀(jì)50年代以來(lái)經(jīng)歷過(guò)多輪油氣勘探和物探技術(shù)攻關(guān),一直未獲重大突破,主要原因是難以對(duì)深層構(gòu)造開(kāi)展有效評(píng)價(jià)。

2010年在該地區(qū)鉆探的砂37井獲得高產(chǎn)工業(yè)油氣流,打開(kāi)了勘探新局面。2011—2012年在該地區(qū)實(shí)施的復(fù)雜山地三維勘探攻關(guān),地震資料品質(zhì)得到了更大的提升,準(zhǔn)確發(fā)現(xiàn)和落實(shí)了多個(gè)勘探目標(biāo),據(jù)此部署的預(yù)探井和評(píng)價(jià)井均獲成功,極大地推進(jìn)了本地區(qū)的勘探步伐[13]。

2006年在柴達(dá)木盆地西部花土溝—獅子溝地區(qū)實(shí)施了405km2高精度重磁、三維大地電磁勘探工作,這是將該地區(qū)作為目標(biāo)靶區(qū)首次實(shí)施的三維重磁電采集、處理、解釋一體化油氣勘探攻關(guān)工作。經(jīng)過(guò)本項(xiàng)工作,深化了對(duì)該地區(qū)斷裂構(gòu)造、基底發(fā)育特征的認(rèn)識(shí),研究分析了構(gòu)造分布對(duì)油氣圈閉的控制作用,進(jìn)一步落實(shí)了有利圈閉。

針對(duì)柴達(dá)木盆地西部山地地區(qū)所具有的典型地表、構(gòu)造雙重復(fù)雜特點(diǎn),研究總結(jié)重磁電一體化綜合地球物理勘探技術(shù)在本地區(qū)油氣勘探的應(yīng)用效果,具有重要的現(xiàn)實(shí)意義。

1 研究區(qū)巖石物性分析

對(duì)花土溝—獅子溝地區(qū)野外地面地層露頭樣品的采測(cè)、鉆井巖心密度和磁化率進(jìn)行實(shí)測(cè),對(duì)電阻率測(cè)井等資料進(jìn)行綜合分析研究,得到本地區(qū)地層巖性、密度、磁化率、電阻率等參數(shù)值,如表1所示。

1) 密度。上油砂山組及更新地層平均密度不大于2.26g/cm3,下油砂山組地層平均密度為2.36g/cm3,上干柴溝組(N1)平均密度為2.44g/cm3,下干柴溝組(E3)以下地層平均密度不小于2.57g/cm3。上、下干柴溝組之間的密度界面構(gòu)成本區(qū)的主要密度分界面,它們之間有0.13g/cm3的密度差。古近系及中生界與古生界基底之間存在較明顯的密度界面,密度差達(dá)0.11g/cm3,是本區(qū)的另一個(gè)重要密度界面,但是它埋藏深,主要產(chǎn)生區(qū)域異常,其局部凸起產(chǎn)生的局部異常幅度可能較上、下干柴溝組之間的密度界面要弱。另外,下油砂山組與上干柴溝組之間有0.08g/cm3的密度差,為次要的密度分界面。

表1 研究區(qū)綜合物性特征

2) 磁化率。本區(qū)新生界和中生界沉積巖地層的平均磁化率均在30×10-5SI以下,屬無(wú)磁或弱磁性地層。前中生界結(jié)晶基底具有一定磁性,為本區(qū)的磁性層。其中,下元古界(Pt1dk)中的片麻巖平均磁化率為177×10-5SI;斜長(zhǎng)角閃巖的平均磁化率為352×10-5SI;奧陶系安山巖火成巖段的平均磁化率達(dá)到了534×10-5SI。另外,本區(qū)發(fā)育的花崗巖磁化率一般小于2×10-5SI,為無(wú)磁性巖石。

圖1 研究區(qū)電阻率測(cè)井曲線(xiàn)

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

在研究區(qū)完成405km2高精度重磁及三維電法勘探工作。重力、磁力測(cè)線(xiàn)67條,剖面總長(zhǎng)807.4km,物理點(diǎn)4104個(gè);電法測(cè)線(xiàn)70條,累計(jì)長(zhǎng)度843.6km,物理點(diǎn)1711個(gè)。采用Trimble 4600lsGPS定位儀和快速靜態(tài)法進(jìn)行測(cè)點(diǎn)定位。

2.1 數(shù)據(jù)采集

本研究工作中,重力和磁法采用面元采集方式,施測(cè)測(cè)網(wǎng)為0.5km×0.2km。重力使用4臺(tái)高精度LCR-G/D型重力儀,磁法使用6臺(tái)G-856/858型銫光泵磁力儀。儀器投入使用前,按規(guī)范要求進(jìn)行了靜態(tài)、動(dòng)態(tài)和一致性試驗(yàn)。重力普線(xiàn)施工前,在研究區(qū)中西部和東南部布設(shè)了重力基線(xiàn)網(wǎng),形成3個(gè)閉合圈,北部施工時(shí),由于地形條件異常復(fù)雜,基點(diǎn)隨施工逐步延伸,形成第4個(gè)閉合圈。基點(diǎn)重力值引自柴達(dá)木盆地I級(jí)重力基點(diǎn)網(wǎng)的CDM-G45,CDM-G46兩個(gè)重力基點(diǎn)。普點(diǎn)觀測(cè)面元采用單日單次觀測(cè)當(dāng)日閉合基點(diǎn)的方法進(jìn)行施測(cè)。為克服復(fù)雜地形的影響,對(duì)所有測(cè)點(diǎn)均使用森林羅盤(pán)測(cè)量了測(cè)點(diǎn)8個(gè)方位的高度角,經(jīng)轉(zhuǎn)換求取測(cè)點(diǎn)8個(gè)方位的高差,據(jù)此求取測(cè)點(diǎn)的近區(qū)地形改正值。磁法在測(cè)區(qū)共設(shè)立1個(gè)磁力日變站,有效控制半徑40km。磁力測(cè)點(diǎn)采用單日單次觀測(cè)法施測(cè),探頭高度1.5m。儀器的動(dòng)態(tài)觀測(cè)精度均小于設(shè)計(jì)要求的0.030×10-5m/s2,動(dòng)態(tài)性能良好。

大地電磁資料采集,采用小面元三維采集技術(shù)[14](圖2),顯示出了良好的勘探效果。本方法的布極方式為:電場(chǎng)的布極采用2×2或3×3電極首尾相接的陣列展開(kāi),采集站之間以GPS實(shí)現(xiàn)同步,小面元大小根據(jù)采集站數(shù)量和地形條件而定,采集站個(gè)數(shù)少及地形差時(shí),可以只布設(shè)2×2面元,即4個(gè)電場(chǎng)采集站和一個(gè)磁場(chǎng)采集站。由于電場(chǎng)同步連續(xù)采集,因此,可以在時(shí)間域?qū)崿F(xiàn)面元內(nèi)數(shù)據(jù)噪聲分析和靜態(tài)效應(yīng)的三維校正處理。同樣,磁場(chǎng)也可以做類(lèi)似的分析和處理。這種方法較電磁法具有更強(qiáng)的抗干擾能力,能有效提高資料采集質(zhì)量。

圖2 三維規(guī)則大地電磁(MT)面元多道網(wǎng)格式采集示意

2.2 數(shù)據(jù)處理

針對(duì)研究區(qū)復(fù)雜地表和地質(zhì)條件,首先采用了變密度中間層地形改正技術(shù)。由于近地表地層密度不同以及地形起伏,若使用全區(qū)統(tǒng)一的常密度中間層及地形改正時(shí),會(huì)引起山形異?；蜱R像異常等虛假重力異常。變密度重力外部改正技術(shù)可消除近地表密度與地形變化引起的虛假重力異常,應(yīng)用該項(xiàng)技術(shù)較好地消除了本區(qū)近地表地形變化與密度變化引起的重力異?；儭?/p>

其次,針對(duì)勘探目標(biāo)層采用了異常剝離技術(shù)。由于布格重力異常、剩余重力異常(常規(guī)處理)、重力垂直二次導(dǎo)數(shù)異常均不易將深、淺層構(gòu)造的影響較好地分離開(kāi),使得在研究中分析古近系目標(biāo)構(gòu)造特征存在一定困難。異常剝離技術(shù)的應(yīng)用能將與目標(biāo)層構(gòu)造因素?zé)o關(guān)的重力異常剝離開(kāi)去,使剝離出的重力異常盡可能地降低受其它層位地質(zhì)因素影響的程度,突出目標(biāo)層構(gòu)造對(duì)應(yīng)的重力異常。

最后,針對(duì)三維地質(zhì)目標(biāo)采取三維數(shù)據(jù)反演技術(shù)。采用相對(duì)密度差反演方法和多次迭代逐次逼近的方式以及極值范圍、反演空間范圍約束等方法對(duì)三維地質(zhì)目標(biāo)體進(jìn)行反演。圖3a為反演模型的重力異常,可以看到,反演結(jié)果較好地?cái)M合了目標(biāo)重力異常(圖3b),擬合精度達(dá)到0.10×10-5m·s-2。

圖3 三維反演模型(a)與實(shí)測(cè)(b)重力異常對(duì)比

在研究區(qū)三維電法勘探現(xiàn)場(chǎng)進(jìn)行的處理工作與常規(guī)預(yù)處理和二維電法勘探類(lèi)似。針對(duì)所實(shí)施的三維電法勘探,所采用的特殊處理技術(shù)主要包括兩方面:①三維規(guī)則化插值,在復(fù)雜山區(qū)施工過(guò)程中,有時(shí)需要采用不規(guī)則面元,但在數(shù)據(jù)后期處理中需要使用規(guī)則化后的數(shù)據(jù),因而需要對(duì)非節(jié)點(diǎn)上的不規(guī)則面元采集數(shù)據(jù)進(jìn)行規(guī)則化處理,三維插值在充分考慮了近區(qū)規(guī)則點(diǎn)和不規(guī)則測(cè)點(diǎn)實(shí)際值的基礎(chǔ)上,在三維空間進(jìn)行由高頻至低頻的插值計(jì)算;②三維靜校正處理,采用空間域低系數(shù)三維統(tǒng)計(jì)濾波靜校正方法,考慮了平面上局部不均勻體的影響和分布。對(duì)比靜校正前、后的視電阻率異常斷面(圖4)可以看到,靜位移引起的斷面上局部直立和鋸齒狀異常(圖4a),在靜校正后的視電阻率異常斷面上(圖4b)得到了有效壓制,表明了三維空間域統(tǒng)計(jì)濾波靜校正方法的有效性。

圖4 靜校正前(a)、后(b)視電阻率異常斷面對(duì)比

三維反演問(wèn)題是電法勘探中最具挑戰(zhàn)性的問(wèn)題之一。20世紀(jì)70年代三維問(wèn)題就已經(jīng)被提出,進(jìn)入21世紀(jì)以來(lái),隨著計(jì)算機(jī)技術(shù)的發(fā)展特別是計(jì)算機(jī)集群的出現(xiàn),三維電法勘探技術(shù)取得了實(shí)質(zhì)性的進(jìn)展。本項(xiàng)研究工作與二維MT/CEMP勘探的主要區(qū)別是實(shí)現(xiàn)了三維反演處理。核心算法是基于正演模擬算子的擬線(xiàn)性近似[15],形成一個(gè)修正的電導(dǎo)率張量線(xiàn)性方程,修正的電導(dǎo)率張量、反射率張量和復(fù)雜的異常體電導(dǎo)率成正比,利用正則化共軛梯度法解線(xiàn)性方程。具體實(shí)現(xiàn)采用非線(xiàn)性共軛梯度法大地電磁擬三維反演[16-20],該方法選取共軛梯度反演算法為擬三維反演的核心,在計(jì)算靈敏度矩陣時(shí),借鑒近似靈敏度矩陣思想,采用一維靈敏度矩陣來(lái)代替三維靈敏度矩陣,并對(duì)非測(cè)點(diǎn)的靈敏度元素提出一種近似方法。在第一次反演之后,采用擬牛頓法更新靈敏度矩陣。

在三維反演實(shí)現(xiàn)過(guò)程中,全區(qū)均一網(wǎng)格的MT測(cè)點(diǎn)結(jié)果都參與反演,反演結(jié)果無(wú)測(cè)線(xiàn)之間閉合差問(wèn)題,避免了二維反演時(shí)主測(cè)線(xiàn)與聯(lián)絡(luò)測(cè)線(xiàn)反演結(jié)果不閉合、忽略旁側(cè)效應(yīng)影響等局限性(圖5)。

圖5 二維(a)與三維(b)電法反演主測(cè)線(xiàn)與聯(lián)絡(luò)測(cè)線(xiàn)閉合差對(duì)比

對(duì)比研究區(qū)MT22線(xiàn)二維反演(圖6a)與三維反演(圖6b)切片,可以看到,兩者都存在電性層構(gòu)造異常分布,但在三維反演結(jié)果中,電性層構(gòu)造異常形態(tài)及特征更為突出,能夠?yàn)闃?gòu)造解釋工作提供重要依據(jù)。

在三維反演中考慮了各測(cè)點(diǎn)的體積效應(yīng),對(duì)局部構(gòu)造異常反應(yīng)更為靈敏,對(duì)局部構(gòu)造刻畫(huà)更精細(xì),反演的基底電性結(jié)構(gòu)更為合理。圖7為過(guò)獅23井等多口井的三維反演斷面,沿?cái)嗝娓骺诰沂镜牡貙悠鸱c反演電性層起伏相關(guān)性良好,三維反演結(jié)果與鉆井揭示的地層起伏趨勢(shì)一致。

研究區(qū)所獲得的三維反演結(jié)果表明,三維反演優(yōu)勢(shì)較之二維明顯,充分體現(xiàn)了三維反演的先進(jìn)性。三維反演結(jié)果為進(jìn)一步的資料綜合解釋奠定了基礎(chǔ)。

圖6 MT22測(cè)線(xiàn)二維(a)與三維(b)反演切片對(duì)比

圖7 過(guò)獅23井等多口井的MT22測(cè)線(xiàn)三維反演切片

3 應(yīng)用效果及地質(zhì)意義

圖8a為三維重力反演得到的密度平面分布,反映了古近系頂面至基底頂面之間密度界面的起伏變化,與常規(guī)處理得到的剩余重力異常相比,消除了其它地質(zhì)因素的影響。研究區(qū)中部目標(biāo)層密度切片中存在一個(gè)南北向高密度異常帶,與電法顯示的構(gòu)造帶對(duì)應(yīng),而西北角的局部高密度也代表了古近系頂面局部構(gòu)造的存在,油砂山構(gòu)造一帶存在局部高密度帶,顯示了油砂山深層存在構(gòu)造。

圖8b為三維磁力反演得到的等效磁化率分布,分析巖石磁性得出,該磁化率切片主要反映前中生界基底磁性差異變化,這種變化與基底結(jié)構(gòu)有關(guān)。

圖8 三維重力反演得到的密度平面分布(a)和磁力反演得到的磁化率(b)平面分布

圖9a為三維反演電阻率為16Ω·m的海拔埋深平面圖,反映了本區(qū)低阻沉積地層底界的分布特征和厚度變化。主要特征為,在工作區(qū)中部沿躍36井區(qū)至獅20井區(qū)近南北向?yàn)榈妥鑼映练e較薄的區(qū)域,底界埋深達(dá)1300～2600m,低阻沉積地層較厚的區(qū)域位于研究區(qū)東部、西部,底界埋深超過(guò)4500m。

圖9b為反演電阻率為60Ω·m的海拔埋深平面圖,基本反映了本區(qū)深層高阻基底頂面的起伏形態(tài)。主要特征與低阻層底界埋深類(lèi)似,研究區(qū)中部沿躍36井區(qū)至獅20井區(qū)為近南北向展布的高阻基底隆起區(qū),在該隆起區(qū)上,分別以獅25井、游深3井北和建參1井北區(qū)域?yàn)橹行?形成3個(gè)明顯的高點(diǎn)圈閉范圍,其高點(diǎn)埋深分別約為4400,4400,5000m,反映了該隆起帶上深層構(gòu)造較發(fā)育的特點(diǎn)。

圖10展示了三維重磁電和地震資料綜合解釋的結(jié)果。圖10a為利用地震剖面進(jìn)行建模得到的綜合解釋結(jié)果,其保留了地震反映清楚的中淺層信息,對(duì)不同的地層填充密度值,然后在綜合解釋平臺(tái)上進(jìn)行重力反演,通過(guò)反復(fù)擬合,逐步修訂深層模型,當(dāng)計(jì)算重力異常曲線(xiàn)和實(shí)測(cè)重力曲線(xiàn)擬合很好時(shí),獲得了深層地質(zhì)模型,這一步驟綜合了地震和重力的信息,為解決深層目標(biāo)結(jié)構(gòu)提供了有效的途徑;圖10b展示了利用三維重磁電資料綜合解釋的N1底構(gòu)造,從圖中可知,獅子溝構(gòu)造帶受滑脫斷裂影響,其上、下構(gòu)造層有一定的差異,上構(gòu)造層形成4個(gè)斷鼻構(gòu)造而下盤(pán)形成4個(gè)完整的局部構(gòu)造;圖10c的三維重力反演高密度體非常明顯地體現(xiàn)了獅子溝構(gòu)造帶上具有4個(gè)單獨(dú)的局部構(gòu)造,向南延伸至油砂山;圖10d是根據(jù)重磁電綜合建模信息開(kāi)展三維地震攻關(guān)后獲取的T2層構(gòu)造圖,從圖中可以看出,三維地震揭示的地質(zhì)結(jié)構(gòu)和前期開(kāi)展的三維重磁電成果(圖10b,圖10c)具有很好的可比性。

圖9 三維反演低阻電性層底界(a)與深層高阻電性層頂面(b)構(gòu)造

圖10 三維重磁電綜合解釋對(duì)比a 過(guò)三維電法39線(xiàn)的綜合解釋剖面; b N1底構(gòu)造立體顯示; c 三維重力反演高密度體分布; d 柴西南T2反射層構(gòu)造

綜合研究表明,研究區(qū)基底巖性主要由花崗巖和古生界淺變質(zhì)巖組成,花崗巖分布于研究區(qū)中部,周?chē)鸀楣派鐪\變質(zhì)巖,這種基底結(jié)構(gòu)特征對(duì)研究區(qū)的沉積巖分布和構(gòu)造發(fā)育有重要控制作用。研究區(qū)以油砂山—獅子溝大斷裂為界,分布深、淺兩套構(gòu)造層,深層構(gòu)造主要有兩個(gè)構(gòu)造帶,一個(gè)為沿花土溝—游園溝—建設(shè)溝近南北向展布的構(gòu)造帶,深淺層構(gòu)造有較好的繼承性;另一個(gè)沿獅35—獅15井北展布的北西向構(gòu)造帶,基底表現(xiàn)為凹陷特征,上部構(gòu)造具有繼承性。淺構(gòu)造層受油砂山—獅子溝大斷裂控制,形成4個(gè)斷鼻。

綜合評(píng)價(jià)表明,干柴溝—游園溝—建設(shè)溝近南北向構(gòu)造帶是進(jìn)一步勘探的有利構(gòu)造帶,游園溝高點(diǎn)為該構(gòu)造帶進(jìn)一步勘探的有利目標(biāo);獅北構(gòu)造帶是研究區(qū)較有利的勘探構(gòu)造帶。

4 結(jié)束語(yǔ)

1) 在柴達(dá)木盆地西部復(fù)雜地區(qū)開(kāi)展三維重磁電綜合勘探攻關(guān),采用小面元電法采集、變密度地形校正、重力異常剝離、三維MT插值、三維空間域統(tǒng)計(jì)濾波靜校正等處理技術(shù),配合三維重磁電綜合反演解釋,有效地克服了復(fù)雜地形、靜態(tài)效應(yīng)等因素的影響,取得了較好的應(yīng)用效果。

2) 三維MT勘探相對(duì)于二維MT而言,測(cè)點(diǎn)分布均勻,能夠更有效地通過(guò)反演獲得研究區(qū)深、淺層的構(gòu)造形態(tài)。三維MT反演考慮體積效應(yīng),其結(jié)果比二維MT結(jié)果對(duì)構(gòu)造的刻畫(huà)更精細(xì),在反映深層構(gòu)造形態(tài)和基底電性合理分布方面具有突出的優(yōu)勢(shì)。將重磁電資料與地震資料在統(tǒng)一平臺(tái)上進(jìn)行三維顯示與解釋,增強(qiáng)了三維MT資料解決中深層構(gòu)造的能力,進(jìn)一步提高了研究成果的客觀性、可靠性。

3) 在類(lèi)似復(fù)雜地區(qū)開(kāi)展油氣資源調(diào)查評(píng)價(jià)、勘探工作中,三維重磁電一體化勘探技術(shù)可作為地震勘探技術(shù)的重要補(bǔ)充,配合地震開(kāi)展構(gòu)造研究,提供有利勘探目標(biāo)。

致謝:感謝中國(guó)石油天然氣集團(tuán)公司高瑞祺、查全衡兩位專(zhuān)家對(duì)本文工作的支持、指導(dǎo)和幫助。

[1] 凌云,郭向宇,高軍,等.油藏地球物理面臨的技術(shù)挑戰(zhàn)與發(fā)展方向[J].石油物探,2010,49(4):319-335 LING Y,GUO X Y,GAO J,et al.The technical challenges on the development trend of reservoir geophysics[J].Geophysical Prospecting for Petroleum,2010,49(4):319-335

[2] 馬永生,張建寧,趙培榮,等.物探技術(shù)需求分析及攻關(guān)方向思考[J].石油物探,2016,55(1):1-9 MA Y S,ZHANG J N,ZHAO P R,et al.Requirement analysis and research direction for the geophysical prospecting technology of SINOPEC[J].Geophysical Prospecting for Petroleum,2016,55(1):1-9

[3] 傅紹乾,蘭秀麗.地震地質(zhì)條件復(fù)雜地區(qū)高精度重力勘探的應(yīng)用效果[J].石油物探,1997,36(1):116-123 FU S Q,LAN X L.The high accuracy gravity survey in complex seismic geology area[J].Geophysical Prospecting for Petroleum,1997,36(1):116-123

[4] 張瑋,詹仕凡,張少華,等.石油地球物理勘探技術(shù)進(jìn)展與發(fā)展方向[J].中國(guó)工程科學(xué),2010,12(5):97-101,112 ZHANG W,ZHAN S F,ZHANG S H,et al.Technical progress and development trend of petroleum geophysical prospecting in China[J].Engineering Science,2010,12(5):97-101,112

[5] 李琴,陳立卿.基于重磁電資料的TC盆地地層結(jié)構(gòu)特征研究[J].石油物探,2013,52(1):105-112 LI Q,CHEN L Q.Study on stratum structure characteristics in TC Basin by using gravity,magnetic and electric data[J].Geophysical Prospecting for Petroleum,2013,52(1):105-112

[6] 陳思宇,雷宛,趙倩倩,等.綜合物探方法聯(lián)合正、反演在海相碳酸鹽巖油氣遠(yuǎn)景區(qū)評(píng)價(jià)中的應(yīng)用[J].石油物探,2013,52(4):438-444 CHEN S Y,LEI W,ZHAO Q Q,et al.Application of integrated geophysical methods jointed forwarding and inversion technology for the evaluation in marine carbonate hydrocarbon potential area[J].Geophysical Prospecting for Petroleum,2013,52(4):438-444

[7] 李慶忠.石油物探的歷史功績(jī)[J].石油地球物理勘探,2013,48(增刊1):191-201 LI Q Z.The role of geophysical prospecting in Chinese petroleum exploration[J].Oil Geophysical Prospecting,2013,48(S1):191-201

[8] 王喜雙,文百紅,王曉帆,等.中國(guó)石油近年來(lái)非地震勘探技術(shù)應(yīng)用實(shí)例及展望[J].中國(guó)石油勘探,2005,10(5):34-40 WANG X S,WEN B H,WANG X F,et al.Petro-China’s latest application examples of non-seismic exploration and prospects[J].China Petroleum Exploration,2005,10(5):34-40

[9] 曲壽利,戴明剛.從第78屆SEG年會(huì)看非地震物探技術(shù)的發(fā)展[J].石油物探,2009,48(5):436-444 QU S L,DAI M G.View on the development of gravity,magnetic and electric exploration techniques from 78thSEG annual meeting[J].Geophysical Prospecting for Petroleum,2009,48(5):436-444

[10] 楊輝,文百紅,戴曉峰,等.火山巖油氣藏重磁電震綜合預(yù)測(cè)方法及應(yīng)用[J].地球物理學(xué)報(bào),2011,54(2):286-293 YANG H,WEN B H,DAI X F,et al.Comprehensive predication of hydrocarbon deposits in volcanic rock by gravity,magnetic,electrical and seismic data and its application[J].Chinese Journal of Geophysics,2011,54(2):286-293

[11] 李起程,彭章池.柴西復(fù)雜山地地震勘探效果及存在問(wèn)題分析[J].青海石油,2005,23(4):12-16 LI Q C,PENG Z C.Seismic exploration effect and the problem analysis in the complex mountainous area in Western Qaidam Basin[J].Qinghai Oil,2005,23(4):12-16

[12] 胡杰,張金崗,張立軍,等.柴達(dá)木盆地復(fù)雜山地采集技術(shù)及效果[J].中國(guó)石油勘探,2006,11(6):51-58 HU J,ZHANG J G,ZHANG L J,et al.Acquisition technique and its effect of complex mountanious region in Qaidam Basin[J].China Petroleum Exploration,2006,11(6):51-58

[13] 寧宏曉,胡杰,章多榮,等.柴達(dá)木英雄嶺復(fù)雜山地三維地震勘探技術(shù)[J].石油科技論壇,2012,31(12):26-39 NING H X,HU J,ZHANG D R,et al.3D seismic prospecting technology for complex mountainous areas in Yingxiongling,Qaidam Basin[J].Oil Forum,2012,31(12):26-39

[14] 何展翔,劉云祥,劉雪軍,等.三維綜合物化探一體化配套技術(shù)及應(yīng)用效果[J].石油科技論壇,2008,27(5):49-54 HE Z X,LIU Y X,LIU X J,et al.3D integrated geophysical-geochemical exploration technique and application effects[J].Oil Forum,2008,27(5):49-54

[15] ZHDANOV M S,FANG S,HURSAN G.Electromagnetic inversion using quasi-linear approximation[J].Geophysics,2000,65(5):1501-1513

[16] 戴世坤,徐世浙.MT二維和三維連續(xù)介質(zhì)快速反演[J].石油地球物理勘探,1997,32(3):305-317 DAI S K,XU S Z.Rapid inversion of magnetotelluric data for 2D and 3D continuous media[J].Oil Geophysical Prospecting,1997,32(3):305-317

[17] 譚捍東,余欽范,BOOKER J,等.大地電磁三維快速松弛反演[J].地球物理學(xué)報(bào),2003,46(6):850-855 TAN H D,YU Q F,BOOKER J,et al.Three dimensional rapid relaxation inversion for the magnetotelluric method[J].Chinese Journal of Geophysics,2003,46(6):850-855

[18] 何展翔,王永濤,劉云祥,等.綜合物探技術(shù)新進(jìn)展及應(yīng)用[J].石油地球物理勘探,2005,40(1):108-112 HE Z X,WANG Y T,LIU Y X.New progress and application of integrative geophysical survey techniques[J].Oil Geophysical Prospecting,2005,40(1):108-112

[19] 胡祖志,胡祥云.三維大地電磁反演綜述[J].地球物理學(xué)進(jìn)展,2005,20(1):214-220 HU Z Z,HU X Y.Review of three dimensional magnetotelluric inversion methods[J].Progress in Geophysics,2005,20(1):214-220

[20] 胡祖志,胡祥云,何展翔.大地電磁非線(xiàn)性共軛梯度擬三維反演[J].地球物理學(xué)報(bào),2006,49(4):1226-1234 HU Z Z,HU X Y,HE Z X.Pseudo three dimensional magnetotelluric inversion using nonlinear conjugate gradients[J].Chinese Journal of Geophysics,2006,49(4):1226-1234

(編輯：陳 杰)

The application of the integrated exploration with gravity,magnetic, and magnetotelluric in complex area,western Qaidam Basin

ZHANG Chunhe1,2,YIN Chengming1,2,GAN Guiyuan3,WANG Caifu4,SUN Weibin4,LIU Yunxiang4, HU Zuzhi4

(1.Oil&GasSurveyCGS,Beijing,Beijing100029,China;2.KeyLaboratoryofUnconventionalPetroleumGeology,CGS,Beijing100029,China;3.ExplorationandDevelopmentResearchInstituteofPetroChinaQinghaiOilfieldCompany,Dunhuang736202,China;4.BGPInc.,CNPC,Zhuozhou072751,China)

Some of the rugged mountains such as the western Qaidam Basin are important prospective exploration targets for hydrocarbon resources in China.Owing to the complicated topography and geo-structure in those areas,hydrocarbon exploration is difficult.To solve the problem,the integrated exploration with 3D gravity,magnetic and electromagnetic method was carried out firstly in the western Qaidam Basin.The routine area measurement is adopted in data acquisition for gravity and magnetic surveys.Data acquisition technology with independent bins was carried out in 3D MT.3D statistical low-pass filtering which based on spatial domain has been used to eliminate the influence of terrain and surface lithology.We carried out 3D inversion on gravity-magnetic-electric data and performed 3D comprehensive interpretation to 3D inversion results under the constraints of drilling data,which reveals the distribution of resistivity and density of target area.

hydrocarbon exploration,rugged mountains,Qaidam Basin,gravity prospecting,magnetic prospecting,magnetotelluric method,inversion

2016-04-01;改回日期：2017-01-19。

張春賀(1966—),男,教授級(jí)高級(jí)工程師,主要從事油氣資源調(diào)查方法技術(shù)研究和相關(guān)技術(shù)標(biāo)準(zhǔn)研制工作。

國(guó)家油氣專(zhuān)項(xiàng)“全國(guó)油氣資源戰(zhàn)略選區(qū)調(diào)查與評(píng)價(jià)(2004XQ)”和中國(guó)地質(zhì)調(diào)查局地質(zhì)礦產(chǎn)調(diào)查評(píng)價(jià)專(zhuān)項(xiàng)項(xiàng)目(DD20160094-1)聯(lián)合資助。

P631

A

1000-1441(2017)04-0607-10

10.3969/j.issn.1000-1441.2017.04.018

This research is financially supported by the National Project of Strategic Area Survey and Evaluation for Oil and gas Resources(Grant No.XQ2004)and CGS Mineral Prospecting and Assessment Program (Grant No.DD20160094-1).