曾方祿

(中國煤炭科工集團 西安研究院有限公司,陜西 西安 710077)

基于海洋MCSEM監(jiān)測的海底油氣后期開發(fā)的可行性

曾方祿

(中國煤炭科工集團 西安研究院有限公司,陜西 西安 710077)

為了探討海洋可控源電磁法(MCSEM)監(jiān)測海底油氣后期開發(fā)的可行性,根據(jù)油氣運移及開采過程中儲層電阻率變化特點,設計并利用MCSEM一維正演程序試算不同的海底電阻率模型,包括海底沉積層僅有儲層的簡單模型和儲層上方有高阻覆蓋層的復雜模型,分別將接收極沉放海底和油氣鉆井,分析電偶極子源在海底和井中的響應特征.結果表明:傳統(tǒng)深海頻率域接收模式MVO曲線能較好地反應儲層電阻率的變化;井中接收頻率域電場對于高阻層位的變化非常敏感;淺海下MCSEM在不同偏移距的時域電場信號既能分離空氣波干擾,又能很好地反應儲層電阻率變化.海洋可控源電磁法對于海底油氣后期監(jiān)測具有潛在優(yōu)勢.

海底油氣開發(fā);可控源電磁法 (MCSEM);電阻率模型;海洋電磁;一維正演;監(jiān)測

0 引言

海洋可控源電磁法(MCSEM)最早于20世紀70～80年代開始進行海洋試驗[1-2],但受限于當時硬件技術,試驗并未成功.隨著20世紀末MCSEM儀器研發(fā)的成功,海洋電磁進入飛速發(fā)展時期,在海底油氣前期勘探中發(fā)揮愈來愈重要的作用,并且已經(jīng)開發(fā)一維到三維的正反演程序[3-7].近幾年,由于MCSEM相比于海洋地震成本低,將其應用在海洋油氣后期開發(fā)的監(jiān)測及油氣運移已逐漸成為研究熱點.Martha Lien和Arnold Orange等探討利用MCSEM監(jiān)測海底油氣儲層電阻率變化的可行性[8-9];Alireza Shahin等[10]討論海洋電磁監(jiān)測海底油氣開發(fā)的可能性;Kong F N等[11]在傳統(tǒng)MCSEM基礎上,將電磁接收極放置在海底油氣井中,固定在油氣出現(xiàn)的層位上觀測電場,該方法有利于監(jiān)測海底油氣開發(fā);Carsten S等[12]突破傳統(tǒng)模式,將電磁發(fā)射源放置在海洋油氣井中,將電磁接收極放在海底,該方法能夠有效監(jiān)測海底油氣因開采而發(fā)生的運移;沈金松等定量分析空氣波的干擾[13-14],為深入利用MCSEM勘探及監(jiān)測提供理論指導;Kang S等研究利用MCSEM監(jiān)測海底二氧化碳隔離的可行性[15].

由于海洋可控源電磁勘探技術的應用時間較短,在監(jiān)測油氣上的研究更少.我國傳統(tǒng)疆域界限內(nèi)的海域面積約為300×104km2,海域盆地面積約為140×104km2[16].近些年,不僅在成熟油氣田周邊進行滾動勘探擴大儲量規(guī)模[17],還取得一系列新區(qū)、新領域的突破[18-19].與海洋地震比較,盡管MCSEM成本較低且實用,但是我國海洋電磁技術發(fā)展還沒有取得實質(zhì)性進展,有必要研究適合我國的海洋電磁監(jiān)測技術.筆者基于海洋可控源電磁基本原理,試算多個電阻率模型,分析頻率域及時間域MCSEM在油氣監(jiān)測應用的可行性.

1 MCSEM一維正演算法原理

在地球物理電磁勘探中,由于頻率一般小于105Hz,因此可以忽略位移電流,得到麥克斯韋方程

組[20]:

式中:ω為圓頻率;σ為電導率;μ為磁導率;Js為外加電源項[11].

利用磁矢勢A,可將電場E和磁場B[21]表示為

由式(1)和(2)變換得到磁矢勢亥姆霍茲方程,并做X和Y方向的二維傅氏變換,令,得到波數(shù)域下的磁矢勢方程:

由于海洋電磁采用水平電偶極子,假設電偶極子沿Y方向布置,則磁矢勢A=(0,Ay,Az),相應的波數(shù)域下的磁矢勢采用漢克爾變換,得到

式中:J0為零階貝塞爾函數(shù).求解磁矢勢A后,通過式(2)得到電場和磁場的頻率域響應.

式中:ai、bi、ci、di為各層衰減系數(shù),在迭代計算中可以得到每層衰減系數(shù);z為觀測點深度;zs為源所在深度;源所在的層位δij=1,無源層位δij=0.

余弦數(shù)字濾波算法[22-23]計算式為

式中:Re為電場實部;Δ為采樣間隔,在電磁余弦變換中一般取為第n個余弦變換濾波系數(shù);t為當前時刻.

計算得到電場頻率域響應后,利用式(9)將頻率域結果轉(zhuǎn)換到電場沖激信號的時間域響應.

2 模型試算結果

海洋電磁發(fā)射接收系統(tǒng)見圖1.假設海水電阻率恒為0.3Ω·m,海平面Z坐標為0 m,沉積層電阻率恒為1Ω·m,海底頂面Z坐標為2 000 m,儲層上頂面Z坐標為3 000 m,下底面Z坐標為3 100 m.水平電偶極子發(fā)射機沉放在海水中,距離海底50 m,測線分別布置在海底和井中,激發(fā)頻率為1 Hz.海底接收極沿Y軸每隔500 m測量1次,終點在與發(fā)射機水平距離20 km處,總共為41個接收點.井中接收極沿Z軸每隔20 m測量1次,終點在Z坐標為4 000 m的位置,總共為101個接收點.

圖1 海洋電磁發(fā)射接收系統(tǒng)示意Fig.1 Schematic diagram of transmitter and receiver system of MCSEM

在海底接收時,可以同時接收水平電場和磁場,在垂直井中只能接收垂直電場.在海底油氣后期開發(fā)過程中,由開采引起的油氣運移使得儲層電阻率改變.在儲層電阻率變化時,為研究接收電磁響應變化規(guī)律及監(jiān)測開發(fā)能力,設置儲層電阻率從100Ω·m依次變化為80、60、40、20Ω·m.

2.1 頻率域海底接收

目前最常用的海洋電磁勘探技術是將水平電偶極子放置在海底上方幾十米內(nèi),接收極布置在海底上探測頻率域水平電磁場.為研究海洋電磁勘探技術在油氣開發(fā)過程中的監(jiān)測能力,首先計算在海水深度為2 000 m時,海底接收到的電磁響應在儲層電阻率變化時水平電場及磁場的變化.極化方向的水平電場和垂直極化方向的磁場能夠有效探測到海底油氣,因此只計算電場Y分量振幅Ey和磁場X 分量振幅B x

[3].根據(jù)電阻率模型及儲層變化規(guī)律,計算結果(見圖2(a)、(b))顯示水平電場和磁場發(fā)生明顯變化.

由圖2(a)、(b)可以看出,在靠近發(fā)射源時,所有的曲線重合,說明接收極所得到的信號基本來自于發(fā)射源;隨著收發(fā)距的逐漸拉大,振幅曲線分開,并且儲層電阻率變化越大,相應曲線變化越大,說明將該技術應用于油氣后期開發(fā)監(jiān)測有很大潛力;隨著收發(fā)距的繼續(xù)增大,由于受到空氣波及電磁波衰減等因素的影響,所有曲線再次重合.因此,應該選擇合適的收發(fā)距,避免無用的信號采集.

為驗證該方法在復雜區(qū)域的監(jiān)測能力,在圖1的電阻率模型里,添加一個電阻率為30Ω·m、離海底為100 m、層厚為50 m的高阻層.儲層電阻率從100Ω·m依次變化為80、60、40、20Ω·m,計算結果見圖2(c)、(d).由圖2(c)、(d)可以看出,電場及磁場在復雜區(qū)域也能明顯反映儲層變化,說明海洋電磁勘探技術在監(jiān)測領域的實用性.

2.2 頻率域井中接收

除了將接收電極沉放在海底方式外,也有將接收極沉放在海底油氣井中的發(fā)射方式,以探測頻率域信號.為檢驗該方法監(jiān)測儲層電阻率變化的效果,設發(fā)射機離井距離為1 000 m,發(fā)射機放在離海底50 m處,發(fā)射接收系統(tǒng)見圖1,儲層電阻率變化從100Ω·m依次變化為80、60、40、20Ω·m.

在垂直井中,只能觀測垂直電場,顯示理論計算的垂直電場,為了突出由儲層電阻率變化造成的電場變化,只顯示井深為2 600～4 000 m的垂直電場振幅.根據(jù)電阻率模型,計算結果見圖3(a),垂直電場在目標層有非常強烈的變化.為驗證該方法的實用性,在圖1的電阻率模型里,添加一個電阻率為30Ω·m、離海底為100 m、層厚為50 m的高阻層,計算結果見3(b).

由圖3(a)、(b)可以看出,在儲層出現(xiàn)的層位,井中接收到的垂直電場出現(xiàn)明顯增強的現(xiàn)象,說明海底激發(fā)、井中接收的發(fā)射方式能夠精確地確定儲層層位位置.在地質(zhì)條件簡單或者相對復雜的區(qū)域,儲層電阻率的變化能夠引起相應深度接收到的垂直電場振幅的強烈變化,說明該種發(fā)射方式具有非常大的潛力.

圖2 圖1模型和復雜模型的電磁響應曲線Fig.2 The EM response curve of fig.1 and complex model

圖3 圖1模型和復雜模型的Ez響應曲線Fig.3 The response curve of fig.1 and complex model

2.3 淺海時間域接收

由圖2可以看出,在收發(fā)距大于一定距離后,所有曲線因空氣波干擾而發(fā)生再次重合,海水越淺,再次重合越早[2].目前,由于受到空氣波的強烈干擾,淺海海洋頻率域MCSEM應用受到極大限制,需要在淺海使用時間域MCSEM以解決淺?？諝獠ǜ蓴_的問題.利用余弦數(shù)字濾波算法,將頻率域MCSEM轉(zhuǎn)換為時間域MCSEM,將圖1的海水深度減小為200 m,儲層的頂面Z坐標為900 m,儲層厚度為100 m,海水及沉積層電阻率不變.發(fā)射源在海面以下10 m,接收極在海面以下100 m,分別計算收發(fā)距3 000 m和5 000 m時的電場沖激響應曲線(見圖4).

由圖4可以看出,接收到的時間域信號剛開始出現(xiàn)重合,原因是受到時域空氣波的影響,空氣波的干擾在前期很快結束,后期沖激響應曲線分開,隨著儲層電阻率的降低,沖擊響應曲線幅度也隨之降低,與頻率域的變化有些類似.這反映儲層電阻率的變化能明顯地引起時間域信號的變化,說明時間域淺海海洋電磁勘探技術同樣也可以應用在淺海海洋的油氣監(jiān)測.

圖4 不同收發(fā)距時淺海海洋的Ey沖激響應曲線Fig.4 The Ey impulse reponse of different offset

3 結論

(1)基于海洋電磁水平電偶極子原理,試算海底接收和井中接收方式下的簡單及復雜電阻率模型的電磁響應,研究海洋電磁在海底油氣后期開發(fā)中的應用潛力.

(2)采用頻率域深海海底接收方式時,采集極化水平電場和垂直極化方向的水平磁場需要根據(jù)不同的情況,選擇相應的收發(fā)距,以避免無效數(shù)據(jù)采集.淺海海洋時間域電磁不同收發(fā)距也能夠較好地反映儲層電阻率的變化.

(3)采用垂直井中接收方式時,接收信號雖然受限于垂直電場,但試算結果說明對于精確確定儲層電阻率的變化有優(yōu)勢,能夠為反演建立可靠的模型.

(4)兩種觀測模式能比較好地反映儲層電阻的變化,說明頻率域和時間域海洋電磁技術不僅可以用于海底油氣的前期勘探中,在后期開發(fā)的過程中也能發(fā)揮監(jiān)測油氣運移的作用,使用海洋可控源電磁法為監(jiān)測油氣開發(fā)提供一種新的手段和思路.

[1] 何展翔,余剛.海洋電磁勘探技術及新進展[J].勘探地球物理進展,2008,31(1):1-8.

He Zhanxiang,Yu Gang.Marine EM survey technology and its new advances[J].Progress in Exploration Geophysics,2008,31(1):1-8.

[2] Steven Constable.Ten years of marine CSEM for hydrocarbon exploration[J].Geophysics,2010,75(5):75A67-75A81.

[3] Kerry Key.1D inversion of multicomponent,multifrequency marine CSEM data:methodology and synthetic studies for resolving thin resistive layers[J].Geophysics,2009,74(2):F9-F20.

[4] Yutaka Sasaki,Max A M.Useful characteristics of shallow and deep marine CSEM responses inferred from 3D finite-difference modeling[J].Geophysics,2009,74(5):F67-F76.

[5] Karen Weitemeyer,Guozhong Gao,Steven Constable.The pratical application of 2D inversion to marine controlled-source electromagnetic data[J].Geophysics,2010,75(6):F199-F211.

[6] Nuno Vieira da Silva,Joanna V Morgran.A finite element multifrontal method for 3D CSEM modeling in the frequency domain[J].Geophysics,2012,77(2):E101-E115.

[7] Yutaka Sasaki.3D inversion of marine CSEM and MT data:An approach to shallow-water problem[J].Geophysics,2013,78(1):E59-E65.

[8] Martha Lien,Trond Mannseth.Sensitivity study of marine CSEM data for reservoir production monitoring[J].Geophysics,2008,73(4):F151-F163.

[9] Arnold Orange,Kerry Key,Steven Constable.The feasibility of reservoir monitoring using time-lapse marine CSEM[J].Geophysics,2009,74(2):F21-F29.

[10] Alireza Shahin,Kerry Key,Paul S,et al.Petro-electric modeling for CSEM reservoir characterization and monitoring[J].Geophysics,2012,77(1):E9-E20.

[11] Kong F N,Roth F.Casing effects in the sea to borehole electromagnetic method[J].Geophysics,2009,74(5):77-87.

[12] Carsten S,Edwards R N.Marine downhole to seafloor dipole-dipole electromagnetic methods and the resolution of resistive targets[J].Geophysics,2007,72(2):39-49.

[13] 沈金松,詹林森,王鵬飛,等.海洋電磁勘探中空氣波相互作用機制和水深影響分析[J].中國石油大學學報,2011,35(6):51-58.

Shen Jinsong,Zhan Linsen,Wang Pengfei,et al.Analysis of air wave interaction mechanism and effect of sea water depth on electromagnetic prospecting[J].Journal of China University of Petroleum,2011,35(6):51-58.

[14] 沈金松,詹林森,王鵬飛,等.可控源海洋電磁勘探中空氣波影響的理論分析和數(shù)值模擬[J].地球物理學報,2012,55(7):2473-2488.

Shen Jinsong,Zhan Linsen,Wang Pengfei,et al.Theoretic analysis and numerical simulation of effects of air wave interaction on the marine controlled source electromagnetic exploration[J].Chinese Journal of Geophysics,2012,55(7):2473-2488.

[15] Kang S,Seol S,Byun J.A feasibility study of CO2sequestration monitoring using the CSEM method at a deep brine aquifer in a shallow sea[J].Geophysics,2012,77(2):E117-E126.

[16] 孫龍德,撒利明,董世泰.中國未來油氣新領域與物探技術對策[J].石油地球物理勘探,2013,48(2):317-324.

Sun Longde,Sa Liming,Dong Shitai.New challenges for the future hydrocarbon in China and geophysical technology strategy[J].Oil Geophysical Prospecting,2013,48(2):317-324.

[17] 蔡希源.中國石化天然氣勘探實踐與勘探方向[J].大慶石油學院學報,2010,34(5):1-8.

Cai Xiyuan.Sinopec's natural gas exploration practice and priority[J].Journal of Daqing Petroleum Institute,2010,34(5):1-8.

[18] 朱偉林,米立君,高陽東,等.領域性突破展現(xiàn)中國近海油氣勘探前景——2011年中國近海油氣勘探回顧[J].中國海上油氣,2012,24(1):1-5.

Zhu Weilin,Mi Lijun,Gao Yangdong,et al.Exploratory field breakthroughs exhibit the future of hydrocarbon exploration offshore China:A review of hydrocarbon exploration offshore China in 2011[J].China Offshore Oil and Gas,2012,24(1):1-5.

[19] 朱偉林.中國近海油氣勘探的回顧與思考[J].中國工程科學,2011,13(5):4-9.

Zhu Weilin.China offshore oil and gas exploration:review and thinking[J].Engineering Sciences,2011,13(5):4-9.

[20] Yuguo,Kerry Key.Marine controlled-source electromagnetic methods 2D marine controlled-source electromagnetic modeling:Part 1-An adaptive finite-element algorithm[J].Geophysics,1972,72(2):WA51-WA62.

[21] 考夫曼A A.頻率域和時間域電磁測深[M].北京:地質(zhì)出版社,1987:153-154.

Kaufman A A.Frequency and transient sounding[M].Beijing:Geological Publishing House,1987:153-154.

[22] 李予國,Stenven Constable.淺水區(qū)的瞬變電磁法:一維數(shù)值模擬結果分析[J].地球物理學報,2010,53(3):737-742.

Li Yuguo,Stenven Constable.Transient electromagnetic in shallow water:insights from 1D modeling[J].Chinese Journal of Geophysics,2010,53(3):737-742.

[23] 王華軍.正余弦變換的數(shù)值濾波算法[J].工程地球物理學報,2004,1(4):329-335.

Wang Huajun.Digital filter algorithm of the sine and cosine transform[J].Chinese Journal of Engineering Geophysics,2004,1(4):329-335.

TE331.3;TE332

A

2095- 4107(2014)01- 0031- 06

DOI 10.3969/j.issn.2095-4107.2014.01.005

2013- 11- 12;編輯:任志平

國家自然科學基金青年基金項目(41304116)

曾方祿(1965-),男,高級工程師,主要從事電磁法勘探方面的研究.