沈建國,王春啟,沈永進
(1.天津大學微電子學院,天津300072;2.北京華暉探測科技股份有限公司,北京101300)
瞬變電磁勘探通常在地面發(fā)射和接收(發(fā)射和接收線圈的垂直源距為0),或者在地面發(fā)射井中接收[1]。通常情況下發(fā)射線圈越大,激發(fā)的功率越強,探測的區(qū)域和深度越深[2-3],應用效果越好[4-6]。
隨著人們對深部地層勘探的需求不斷增加,在現(xiàn)有勘探原理的基礎(chǔ)上提出了油井二次開發(fā)、剩余油分布研究、幾何因子淺層勘探等方法[7-11]。因激發(fā)和接收點距離所測量地層較遠,信號弱,故制約勘探效果的因素較多。
本文利用全空間瞬變電磁響應進行深層勘探,具體闡述如下:根據(jù)套管井地層電導率測量原理,將瞬變電磁發(fā)射源(線圈或電極)置于井底,將等間距陣列接收線圈或電極置于發(fā)射井內(nèi)和相鄰的接收井內(nèi),記錄瞬變電磁場在地層中的響應過程,利用井中不同位置接收到的響應波形得到井底以下深部地層的電導率分布情況。
如果能夠通過本井和鄰井已經(jīng)獲得的測井曲線,將陣列接收波形中已知地層的特征去掉,便可利用線圈或電極激發(fā)的瞬變電磁場在全空間各個方向向外傳播、擴散的動態(tài)過程,提取不同時刻波形所攜帶的深部地層電導率及其界面信息,實現(xiàn)電阻率的精細探測。
為此,將通過套管井節(jié)箍處的瞬變電磁響應得到套管井的幾何因子[12],在裸眼井的響應中采用相減的方法去除直接耦合響應[13]。根據(jù)Doll電流環(huán)模型得到井間瞬變電磁渦流激發(fā)響應的全空間幾何因子[14],進行模型實驗[15],測量得到了不同空間位置的瞬變響應波形。井中空間小,在其中進行瞬變電磁激發(fā),主要能量均可進入地層,因此激發(fā)的深部地層有用信號較強,受影響的因素較少。
本文利用瞬變電磁滿足的基本方程及其解析解計算瞬變電磁響應波形,從相位變化反映瞬變電磁能量空間傳播情況的角度分析瞬變電磁響應波形的特點,對比線圈激發(fā)和電極激發(fā)的響應波形,分析波形中界面信息的分布方式,并且利用一維模型展現(xiàn)電導率界面的瞬變電磁響應波形。根據(jù)電磁感應原理研究電極激發(fā)的響應波形以及與地層電導率直接相關(guān)的渦流再次激發(fā)響應。響應波形每一時刻的幅度都包含距離激發(fā)源一定深度地層的渦流再次激發(fā)響應,并與該深度地層的電導率成正比,直接指示該地層電導率。不同時刻套管井測井波形的幅度與裸眼井電導率測井曲線的一致性驗證了上述結(jié)論。
任何瞬變電磁響應均可以視為正弦激勵e-iωt響應的疊加。無論是電極還是線圈,其在全空間(無限大均勻地層)中激發(fā)的瞬變電磁場均按照推遲勢reiγr/r向外擴散、傳播。其中,(地震波)電磁波(電導率為0或者高頻)的響應反映了其傳播過程,速度v為常數(shù),不同位置的響應波形形狀一樣。瞬變電磁勢函數(shù)A在球坐標系下可表示為:
(1)
(2)
(2)式所描述的瞬變電磁A既有相位移動又有幅值衰減,幅值衰減和相位移動均隨頻率改變,且衰減劇烈,因此導致不同源距的響應波形完全不同。如圖1a 所示,該響應為點源激發(fā)的響應(格林函數(shù)),用于描述導電介質(zhì)中瞬變電磁場的傳播特征,響應特征與激發(fā)源類型(線圈、電極)無關(guān),為瞬變電磁場在無限大均勻介質(zhì)中傳播后的響應特征。波形的上升沿包含了豐富的地層信息,但是對于這樣的波形,人們很少從瞬變電磁能量傳播(擴散)的角度對其相位信息及其波形上的延遲進行分析。在套管井地層電導率測井和井底瞬變電磁精細勘探中,瞬變電磁能量的傳播(擴散)展現(xiàn)了響應時間與地層深度之間的對應關(guān)系,有助于我們從不同時刻的響應中獲得不同深度的地層電導率,因此有必要對波形形狀、延遲以及對應的空間分布進行分析。
將圖1a中不同源距的響應波形放大不同倍數(shù)后繪制在一起得到圖1b,從圖1b中可以看出,不同源距的瞬變電磁響應均具有非常明顯的延遲(傳播特征),具體表現(xiàn)為:不同源距的響應起始點、上升沿和峰值等均有明顯的時間差。波形中起始點、上升沿和峰值所在的3個時刻及其對應響應形狀反映了瞬變電磁響應的主要特征。圖1a中響應波形的時間差所反映的傳播過程與電磁波和地震波的傳播過程完全不同:不同源距的響應波形差異大,并且這些差異存在于波形的各個時刻。這個差異使得我們無法像處理地震波和電磁波那樣,根據(jù)相位一致的原則從測量的波形上求取時間,也無法用響應波形的(單個)延遲時間來描述傳播過程。這是瞬變電磁場與電磁波場傳播的根本區(qū)別。
圖1 點源在無限大均勻介質(zhì)中激發(fā)得到的不同源距的響應波形(a)以及放大以后的響應波形(b)
由(2)式可知,瞬變電磁能量經(jīng)過一定的距離傳播后,因不同頻率的響應延遲時間(不是一個常數(shù))[16]和幅度衰減均不同,故所有頻率的響應疊加后形成的響應波形形狀不同。
圖1b中的時間差異由(2)式中指數(shù)的虛部(第2項)決定,相位移動作為影響瞬變電磁響應的第2個因素,隨頻率和電導率變化,并可折算為延遲時間t0:
(3)
圖2 電導率為1S/m的地層集膚深度隨頻率變化的規(guī)律
在裸眼井中激發(fā)瞬變電磁場,通過井內(nèi)的陣列接收線圈或電極接收不同源距的瞬變電磁響應。頻率越低對應的集膚深度越大,瞬變電磁能量能夠到達的地層深度越深,波形變化越慢,這些深層的信息在瞬變電磁響應的后續(xù)波形中有所反應。
圖3為線圈在無限大均勻介質(zhì)中4個不同源距下激發(fā)的瞬變電磁響應(實線、短虛線)及響應差(長虛線)。圖3a是電導率為1S/m時,沖擊函數(shù)激發(fā),源距2m處的響應波形,其延遲時間小于0.5μs,波形上升快(高頻成分豐富)、下降慢且拖尾(高頻成分部分衰減,幅度減小導致),與沖擊函數(shù)的形狀(上升、下降均極快)不一樣。因為延遲時間太短,在近源距通過測量延遲時間獲得地層電導率需要進行難度較大的高速模數(shù)轉(zhuǎn)換器采集。圖3b、圖3c和圖3d分別為:傳播10m時延遲時間接近4μs的響應波形,傳播20m時延遲時間約20μs的響應波形,傳播40m時延遲時間達到60μs的響應波形,即源距越長,響應起始時刻的延遲時間越長。這是因為源距越長,頻率成分衰減越多;頻率越高、衰減越大,高頻的成分衰減較多;剩余的低頻成分傳播速度慢,傳播延遲時間長。波形幅度起始時間(響應為零的時間長度)隨著源距的增加呈非線性增加。無論延遲多長,響應波形均呈現(xiàn)快速上升達到峰值,而后緩慢下降的形態(tài)。因不同源距的波形差異很大,故不能像地震波那樣進行相關(guān)對比和分析,只能根據(jù)波形的特點研究其它信號的處理方式。
圖3 線圈在無限大均勻介質(zhì)中4個不同源距下激發(fā)的瞬變電磁響應及響應差a r=2m; b r=10m; c r=20m; d r=40m
根據(jù)電磁感應原理,線圈激發(fā)的瞬變電磁響應隨時間變化,在地層中激發(fā)的感應電場是閉合的,會產(chǎn)生隨時間變化的渦流并再次激發(fā)響應即二次場響應。因為渦流與地層電導率成正比,故其激發(fā)的響應也與地層電導率成正比,并混疊在瞬變電磁響應中。提取實際測量波形中渦流再次激發(fā)的響應,可以用于研究地層的電導率及其空間分布情況。因此,這種再次激發(fā)的響應為有用信號。
DOLL采用假設(shè)的地層電流環(huán)模型研究了發(fā)射線圈激發(fā)的渦流及其在接收線圈產(chǎn)生的二次場響應,得到了發(fā)射、接收線圈之間的直接耦合響應和渦流再次激發(fā)的二次場響應表達式。對于線圈系結(jié)構(gòu)而言,其直接耦合響應也是確定的,因地層電導率不同導致渦流激發(fā)二次場響應存在差異。將同一源距下不同地層電導率對應的響應相減,可消除直接耦合響應,得到渦流再次激發(fā)響應的差,即不同地層渦流引起的響應差[18-23]。在響應變化最快的位置,響應差出現(xiàn)峰值(圖3)。上述結(jié)論符合電磁感應定律:渦流激發(fā)的二次場響應與渦流在接收線圈引起的磁通量隨與時間的變化率成正比(渦流引起的磁通量與瞬變電磁響應成正比)。因此,渦流激發(fā)的二次場響應與瞬變電磁響應對時間的導數(shù)成正比,此處的瞬變電磁響應對應地層(內(nèi)部)的響應,渦流激發(fā)的二次場響應與地層電導率成正比。
瞬變電磁響應中渦流激發(fā)的二次場響應對于理解響應波形中的有用信號具有重要意義。響應波形是隨時間變化的,在不同時刻有不同的響應幅度。這些響應幅度包含怎樣的有用信號?與哪個深度的地層電導率有對應關(guān)系?
圖3中各個源距的響應差峰值所在時刻是激發(fā)源在該深度位置地層電導率產(chǎn)生最大渦流激發(fā)信號的時刻,渦流再次激發(fā)的二次場響應(渦流到接收線圈的傳播時間很短)直接疊加在此刻的原始響應波形中(與瞬變電磁響應信號相比很小),或者說此時測量的波形值包含直接反映該深度位置的渦流激發(fā)的二次場響應。此時該深度地層的電導率產(chǎn)生的渦流激發(fā)響應最大,所有源距的響應中該時刻的響應值集中反映了該深度地層的電導率。其它深度的地層產(chǎn)生的渦流因為隨時間變化慢,因此激發(fā)的響應較小。當該深度的響應達到峰值,對時間的導數(shù)為0時,其渦流激發(fā)的二次場響應為零。
隨著時間推移,該深度位置的響應(達到極值后)開始下降,對時間的導數(shù)為負值,其渦流值激發(fā)的二次場響應信號為負值,疊加在后續(xù)的響應波形中。因響應下降慢,對時間的導數(shù)小,故此時該深度位置渦流激發(fā)的二次場響應小。此時瞬變電磁能量傳播到較深位置的地層,有些位置的響應剛剛開始上升,渦流開始增加,有些位置的響應早已上升,上述響應對時間的導數(shù)均為正,渦流激發(fā)的正二次場響應疊加在該時刻的響應波形中。在上升最快的深度(地層)位置,地層電導率對渦流激發(fā)的二次場響應的貢獻最大,該時刻的波形值主要用于測量該地層的電導率。
綜上所述,在所有源距的響應波形中,每個時刻的響應值均對應于距離激發(fā)源某個深度的地層區(qū)域的響應,該區(qū)域瞬變電磁響應上升最快時,產(chǎn)生的渦流對時間的導數(shù)最大,其電導率產(chǎn)生的渦流激發(fā)出最大的正二次場響應。對于每個地層激發(fā)深度而言,較其淺的地層產(chǎn)生的渦流對時間的導數(shù)較小,當導數(shù)為負時,渦流激發(fā)出負的二次場響應,較其深的各個地層響應開始增加,但因響應幅度較小,故產(chǎn)生的渦流較小,變化較慢,渦流激發(fā)的正向二次場響應也較小,這些與地層電導率成正比的二次場響應均疊加在瞬變電磁響應波形中。
所有不同源距的響應波形同一時刻的幅度值均對應距離激發(fā)源同一深度位置的地層。該地層為此刻的波形提供最大的渦流激發(fā)響應,比其更深的地層提供的響應幅度較小,比其更淺的地層則提供負的渦流激發(fā)響應。當淺部地層渦流二次激發(fā)的負響應幅度超過較深地層渦流二次激發(fā)的正響應幅度時,總的渦流二次激發(fā)的響應為負,并且隨著時間的增加逐漸減小。
以上結(jié)論適用于地層和井中任意位置的瞬變電磁響應。瞬變電磁響應變化時刻即為接收線圈接收到的渦流激發(fā)二次場響應的時刻。
根據(jù)電磁感應原理,點電極或棒狀電極激發(fā)的瞬變電磁響應在地層中也產(chǎn)生渦流,這種渦流與線圈激發(fā)的渦流存在差異。電極激發(fā)的瞬變電磁場與線圈激發(fā)的瞬變電磁場空間分布相互垂直,因此,在柱坐標系下,線圈激發(fā)的磁場沿半徑和軸線z方向分布,電場沿圓周方向形成閉合的圓環(huán);電極激發(fā)的電場沿半徑和軸線z方向分布,磁場沿圓周方向形成閉合的圓環(huán)。電極激發(fā)的瞬變電磁響應中,電場呈現(xiàn)放射狀,激發(fā)的渦流沿半徑和井軸z方向呈放射狀輻射至無窮遠,其產(chǎn)生的磁場沿圓周方向分布。因為該磁場隨時間變化,故其在接收電極位置同樣會再次激發(fā)電響應,該響應與地層(圓球形狀)電導率有關(guān)。圖4為不同源距下點電極激發(fā)的瞬變電磁響應及響應差,其響應形狀與圖3相似:先快速上升達到峰值以后再慢速下降。具體的差別表現(xiàn)為圖3中兩個不同地層電導率的響應差在響應變化最快的位置出現(xiàn)極大值之后快速減小,最后為很小的負值。在電極激發(fā)的響應中,距離激發(fā)電極某個深度地層的渦流激發(fā)的二次場響應在某個確定的時刻(其響應變化最快、即上升最快)最強,周圍其它地層響應變化比較慢,渦流激發(fā)的二次場響應幅度較小。此時距離激發(fā)電極較近的地層響應開始減小,其地層渦流在響應減小過程中二次激發(fā)的響應為負。因為響應減小時變化比較慢(拖尾),故其地層渦流二次激發(fā)的響應幅度小。距離激發(fā)電極較遠的地層此時開始有響應,且響應幅度逐漸增加,其渦流激發(fā)的響應為正。隨著時間的增加,響應變化最快的地層其渦流激發(fā)的二次場響應逐漸減弱為零并轉(zhuǎn)換極性,而距激發(fā)源較遠處的地層渦流激發(fā)的二次場響應開始達到極大值。在每個時刻,所有地層渦流激發(fā)的二次場響應均疊加在實際測量的響應波形上。每個時刻的響應波形中均包含該時刻瞬變電磁場擴散(涉及)到的全部空間的渦流激發(fā)響應。由于激發(fā)電極周圍介質(zhì)的響應始終存在,每個時刻的響應均在變化,因此空間每一點的響應波形中均包含其渦流激發(fā)響應。利用后續(xù)響應中的渦流再次激發(fā)響應獲取較遠或者深部地層的電導率時,需要將電極近處地層的渦流激發(fā)響應從波形中去除。
圖4 不同源距下點電極激發(fā)瞬變電磁響應及響應差a r=2m; b r=10m; c r=20m; d r=40m
當線圈和電極分別激發(fā)出瞬變電磁場后,由于電場強度和磁場強度均隨時間變化,根據(jù)Maxwell方程的兩個旋度公式可知,瞬變電磁場在空間的分布也是變化的,不會在空間中停留,也不會像直流電磁場那樣在空間各點形成一個穩(wěn)定的分布(不隨時間變化),而是在不同時刻分布在不同的空間中。
圖5a和圖5b是線圈激發(fā)的瞬變電磁場在0.10ms和0.12ms的空間分布情況,隨著時間的增加,瞬變電磁場分布區(qū)域逐漸擴大,向外傳播特征明顯。但當瞬變電磁能量達到一定區(qū)域后,由于地層電導率引起的衰減和相位移動,導致瞬變電磁能量在向外傳播的過程中較高的頻率成分幅度衰減得很快,隨著時間的推移,瞬變電磁場的主要能量仍然集中在發(fā)射源附近的一定區(qū)域內(nèi)(地震波或者電磁可以集中能量一起向外傳播)。由于不同頻率的瞬變電磁能量經(jīng)過了不同程度的衰減(集膚深度),到達一定的區(qū)域時,幅度均已減小,經(jīng)過相移再疊加,幅度均較小。從圖5c和圖5d可以看出,隨著時間的增加,中垂線區(qū)域內(nèi)的瞬變電磁幅度開始逐漸減小,最終消失在近源距的某個區(qū)域內(nèi)。
圖5 線圈激發(fā)的瞬變磁場在0.10ms(a)、0.12ms(b)、0.68ms(c)、0.82ms(d)時刻的分布情況
在圖5的瞬變電磁響應空間分布的等勢面(球形)上,每一點均具有相同的響應波形和渦流。該分布區(qū)域(各個輪廓線)中的瞬變電磁能量隨時間變化存在差異,某個時刻,瞬變電磁能量變化最快的區(qū)域(輪廓線)即為波形中該時刻所測量的區(qū)域。這些區(qū)域內(nèi)隨時間變化的渦流均會在接收線圈上再次激發(fā)響應,區(qū)域內(nèi)每一點渦流再次激發(fā)響應的貢獻與接收線圈所在的空間位置有關(guān),可以用空間幾何因子對其進行描述[14]。
由于瞬變電磁場在界面的響應對地下勘探具有重要意義,因此從實際測量的波形中分離出界面的響應即可確定界面的存在,也可進一步確定另外一個地層的電導率。
我們采用一維介質(zhì)模型正弦激勵eiωt對電導率界面的影響進行研究。忽略介電常數(shù)項,一維瞬變電磁場滿足的微分方程:
(4)
(5)
一維兩層模型如圖6所示,假設(shè)兩層介質(zhì)的電導率分別為σ1、σ2,磁導率μ相同。在左邊的介質(zhì)有激發(fā)源,激發(fā)瞬變電磁(R)場入射到邊界,遇到邊界后發(fā)生反射(F)和透射(T),反射場包含地層界面信息和第二層的地層電導率信息。
圖6 介質(zhì)電導率不同的一維兩層模型
為具體展示一維兩層模型中反射波的形狀和特征,將坐標原點取在圖6所示的兩層介質(zhì)交界位置,并假設(shè)左邊第一層電導率σ1為100S/m,相對磁導率為1,第二層的電導率σ2為90S/m,相對磁導率為1,取頻率步長df=0.01Hz;激發(fā)源到界面的距離為140m,第一層的接收點到界面的距離為10m,第二層的接收點到界面的距離為20m。圖7中紅線是瞬變電磁場經(jīng)過第一層的接收點時的響應波形,黑線是從界面反射后反射場到達接收點時的響應波形(放大10倍),綠線是接收點實際接收到的全波波形。將反射場繼續(xù)放大得到圖7b,可以看出:反射場較入射場有明顯的延遲,該延遲為入射場經(jīng)過接收點到達界面后,從界面反射再次傳播到接收點所經(jīng)過的時間。在此過程中,反射場經(jīng)過了20m的衰減和一次反射,因此其波形不但延遲而且形狀與入射場相比還存在較大的變化(變緩)。此波形與入射場疊加后,對接收波形會產(chǎn)生一定的影響(如圖7中紅線與綠線之間的差別)。因影響較小,并未像地震波那樣構(gòu)成獨立的反射波(地震的反射波是獨立存在的,在時間上與入射波分離),而是界面引起的反射場與幅度較大的入場波同時存在于同一時間段內(nèi),其影響貫穿于整個響應波形中。
圖7 兩層介質(zhì)模型反射波放大10倍(a)以及反射波放大50倍(b)后的入射波和反射波以及全波波形
利用井中激發(fā)瞬變電磁場,并在井內(nèi)一定的源距(區(qū)域)以外沿井軸方向進行陣列接收,瞬變電磁場經(jīng)過接收陣列時會出現(xiàn)響應波形先快速增加達到極值,之后再緩慢減小的變化過程。不論是線圈還是電極激發(fā),不同位置的響應波形均具有這樣的形狀。利用該波形形狀進行深部地層及界面勘探,必須先清楚其中的地層電導率信息分布方式,地層電導率以3種形式存在于陣列接收的響應波形中,具體闡述如下。
1) 響應波形的幅度。沿一定的方向陣列接收到的每個頻率的波形幅度均是按照1/r乘以指數(shù)規(guī)律衰減的,其指數(shù)衰減系數(shù)與地層電導率的開方呈線性關(guān)系。通過求取不同源距接收到的各個頻率的波形幅度的指數(shù)衰減系數(shù)可以獲得地層電導率。
2) 響應波形的相位。沿一定方向陣列接收到的每個頻率的波形相位隨源距呈線性變化,其斜率與地層電導率的開平方成正比。獲得各個頻率的相位隨源距變化的斜率即可獲得地層的電導率。
3) 根據(jù)電磁感應原理,瞬變電磁響應在地層中會產(chǎn)生渦流,渦流隨時間改變會再次激發(fā)新的響應,該響應與地層的電導率成正比。地層響應變化最快的時刻渦流激發(fā)的二次場響應最大。瞬變電磁響應波形的每個時刻均對應某個深度的地層,響應值集中反映該深度地層的電導率。
利用響應波形的幅度和相位獲得的地層電導率隨頻率改變,頻率越低,探測的深度越深,涉及的地層區(qū)域越大。根據(jù)不同頻率的電導率構(gòu)建一個深度剖面,在已知頻率與深度的對應關(guān)系(可以用集膚深度)的情況下,可以獲得地層電導率由淺入深的變化情況。兩種方法均是對整個波形進行處理,是整個波形中不同時刻所包含的地層電導率的綜合提取方法。
利用瞬變電磁場的擴散過程和瞬態(tài)特征,將不同時刻的響應與距離激發(fā)源的深度對應,即可得到響應波形中不同時刻包含的地層電導率信息,該方法的有效性在套管井中得到了驗證。提取套管井中實際測量的瞬變電磁響應波形在10個時刻的幅度(對應于10個徑向深度地層的電導率),并按照測井深度與裸眼井所測量的電導率曲線(RT)綜合繪制,結(jié)果如圖8 所示。圖中黑線是裸眼井測量的電導率曲線,彩色的10條線分別是10個時刻井內(nèi)瞬變電磁測量波形的幅度,兩者的深度存在1m的偏差(黑線向右移動了1m)。右邊的曲線形狀非常相似(左邊的差異是地層擴徑所致)。原始測量波形的幅度與地層電導率相關(guān),反映了地層電導率的變化。
圖8 套管井瞬變電磁測井波形10個時刻的幅度隨深度的變化情況及感應測井電導率曲線
地層電導率界面產(chǎn)生的反射波,直接反映了界面的存在。反射場與入射場疊加在一起貫穿于整個響應波形中。當電導率差異較小時,反射波的幅度也較小。
利用井中激發(fā)、陣列接收瞬變電磁響應波形頻譜的幅度、相位以及渦流再次激發(fā)的響應均可獲得不同深度地層的電導率,根據(jù)響應波形的形狀還可以獲得不同深度電導率界面。波形中每個時刻的幅度與距離激發(fā)源的某個深度地層的電導率對應,即在波形上取某個時刻,此刻某個與其對應的區(qū)域(距離激發(fā)源某個深度)地層響應變化最快,其渦流在每個位置激發(fā)的二次場響應最大。電極和線圈激發(fā)均能產(chǎn)生與地層電導率成比例的二次場響應。地層渦流激發(fā)的二次場響應在其響應(渦流)增加時為正、減小時為負,波形中的二次場響應是各個地層所激發(fā)的不同極性的二次場響應的疊加。界面反射產(chǎn)生的響應與瞬變電磁響應疊加在一起,貫穿于后續(xù)的響應波形中。線圈激發(fā)的瞬變電磁場以偶極子波陣面的形狀擴散,其上各點的渦流相同,均在接收線圈再次激發(fā)二次場(受空間位置影響),該二次場是波形中的有用信號,也是精細勘探的理論基礎(chǔ)。由于瞬變電磁響應的二次場幅度比較小,響應中與地層電導率無關(guān)的直接耦合響應幅度較大,故需要進一步研究去除直接耦合響應的方法。