胡春橋,謝偉彪,殷秋麗,劉迪仁
(1.油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室長江大學(xué),湖北 武漢430100;2.長江大學(xué)地球物理與石油資源學(xué)院,湖北 武漢430100;3.冀東油田勘探開發(fā)研究院,河北 唐山063004)
高分辨率感應(yīng)測井(HRI)采用全對稱線圈系結(jié)構(gòu)[1],縱向分辨率為2ft*非法定計(jì)量單位,1ft=12in=0.3048m,下同,深、中感應(yīng)徑向探測深度分別為91、39in[2],在徑向探測深度和縱向分辨率上都優(yōu)于傳統(tǒng)雙感應(yīng)測井儀。它利用同相和正交電導(dǎo)率分量對圍巖和趨膚效應(yīng)同時進(jìn)行校正[3],測量得到的結(jié)果更接近地層真實(shí)情況。在滲透性地層,徑向探測深度不同的電阻率測井曲線往往會形成幅度差。如果是淡水泥漿,一般在油層處會形成正幅度差,在水層處會形成負(fù)幅度差,可以利用幅度差定性認(rèn)識油水層。1999年柯式鎮(zhèn)等[4]分析了雙側(cè)向測井曲線幅度差的主要影響因素,并在此基礎(chǔ)上討論了引起雙側(cè)向“雙軌”的主要原因。2003年劉國民等[5]分析了利用雙感應(yīng)測井曲線的正、負(fù)幅度差識別水淹層的方法。2004年鄧少貴等[6]分析了泥漿侵入地層中雙感應(yīng)測井曲線正負(fù)差異特性,深、淺電阻率曲線的正負(fù)差異不僅與儲層含油性有關(guān),并受到泥漿濾液、地層水礦化度、地層束縛水飽和度、殘余油氣飽和度等諸多因素的共同影響。對于高分辨率感應(yīng)測井,對其幅度差特性以及“雙軌”方面的研究比較少,因此有必要對高分辨率感應(yīng)測井的幅度差進(jìn)行研究。本文通過數(shù)值模擬的方法,計(jì)算并討論了圍巖、井眼、泥漿侵入等3個方面對高分辨率感應(yīng)測井曲線幅度差的影響。
為了討論高分辨率感應(yīng)測井幅度差的影響因素,建立了具有旋轉(zhuǎn)對稱性的3層介質(zhì)模型(見圖1)。圖1中,h為目的層厚度,dm為井眼半徑,di為侵入深度,Rs為圍巖電阻率,Rm為泥漿電阻率,Ri為侵入帶電阻率,Rt為原狀地層電阻率;R為主接收線圈,TD和Td為1對深感應(yīng)發(fā)射線圈,TM和Tm為1對中感應(yīng)發(fā)射線圈,RB和Rb為1對深、中感應(yīng)屏蔽線圈,所有線圈均繞在同1根玻璃鋼芯棒上。
在垂直井中,可以將儀器發(fā)射線圈看作是Z方向上的磁偶極子(見圖1)。雙線圈系的視電導(dǎo)率為
式中,NT、ST分別為發(fā)射線圈和接收線圈的匝數(shù)。接收線圈中的二次磁場HR可以表示為散射場HS和背景場Hb的迭加[7],即
由積分方程可知,散射場HS可表示為
式中,rR、rT分別為接收線圈和發(fā)射線圈的坐標(biāo)位置;Δ為空間電導(dǎo)率差;GH為磁并矢Green函數(shù)[8-9]
式中,I=diag(1,1,1)為單位并矢,k2=-iωμσ。散射場求解時按積分方程法原理將整個求解區(qū)間剖分成N個小體積單元Ωn(n=1,2,…,N),假設(shè)每個小體積單元內(nèi)介質(zhì)及其電磁場分布均勻。背景場有解析解,電場可以通過單位磁偶極子在均勻介質(zhì)中產(chǎn)生的Hertz勢計(jì)算得到[10]。由式(2)可得到二次磁場,由此依據(jù)高分辨率感應(yīng)測井儀器的復(fù)合線圈系得到視電導(dǎo)率。
在圖1所示地層模型下,不考慮儀器自身影響,計(jì)算高分辨率深感應(yīng)與高分辨率中感應(yīng)在不同地層、井眼條件下的比值(RHRd/RHRm),用兩者比值的變化說明幅度差的變化情況。
圖2是無泥漿侵入(非滲透層)時,在不同的圍巖電阻率Rs情況下,RHRd/RHRm隨地層厚度h的變化曲線。其中,井眼半徑dm為0.1m,原狀地層電阻率Rt和泥漿電阻率Rm均為100Ω·m。
從圖2可見,曲線變化過程中會出現(xiàn)極值點(diǎn),極值點(diǎn)的極性與圍巖電阻率的值有關(guān),低電阻率圍巖會出現(xiàn)極小值點(diǎn),高電阻率圍巖會出現(xiàn)極大值點(diǎn)。高分辨率深感應(yīng)受圍巖影響比中感應(yīng)大,當(dāng)Rs>Rt且地層較薄時,高電阻率圍巖會使深感應(yīng)的值變大。隨著地層厚度逐漸變厚,圍巖對深感應(yīng)的這種影響又會減弱,因此RHRd/RHRm隨h變化時是先增大后減小,繼而出現(xiàn)極大值點(diǎn)。這個過程中高分辨率深、中感應(yīng)測井曲線會出現(xiàn)正幅度差(RHRd>RHRm),并且幅度差是先增大后減小;當(dāng)?shù)貙雍穸冗_(dá)到5m時,RHRd/RHRm基本不再隨地層厚度的變化而變化,2條測井曲線基本保持重合。當(dāng)Rs<Rt且地層較薄時,低電阻率圍巖會使深感應(yīng)的值變小。隨著地層厚度逐漸變厚,圍巖對深感應(yīng)的這種影響又會減弱。因此,RHRd/RHRm隨h變化時是先減小后增大的,繼而出現(xiàn)極小值點(diǎn),這個過程中高分辨率深、中感應(yīng)測井曲線會出現(xiàn)負(fù)幅度差(RHRd<RHRm),并且幅度差先增大后減小。當(dāng)?shù)貙雍穸冗_(dá)到5m時,RHRd/RHRm基本上不再隨地層厚度的變化而變化,但是此時當(dāng)圍巖電阻率太低時,負(fù)幅度差在5~10m范圍內(nèi)依然存在,這是由于圍巖電阻率太低使深感應(yīng)電阻率比中感應(yīng)電阻率低很多造成的。
圖2 不同Rs情況下h對RHRd/RHRm的影響
圖3是無泥漿侵入(非滲透層)時在不同地層厚度h情況下RHRd/RHRm隨圍巖電阻率Rs的變化曲線。從圖3可以看出,當(dāng)Rs<Rt時,高分辨率深、中感應(yīng)測井曲線隨著圍巖電阻率變化時會出現(xiàn)一定程度的負(fù)幅度差,這是由于低電阻率圍巖使得深感應(yīng)測井值比中感應(yīng)測井值小造成的;當(dāng)Rs=Rt時,高分辨率深、中感應(yīng)測井曲線會重合;當(dāng)Rs>Rt時,高分辨率深、中感應(yīng)測井曲線會出現(xiàn)一定程度的正幅度差,這是由于高電阻率圍巖使得深感應(yīng)測井值比中感應(yīng)測井值大造成的。當(dāng)Rs/Rt<0.1或者Rs/Rt>10時,RHRd/RHRm基本隨Rs變化很平緩,此時RHRd/RHRm受Rs的影響較小。
圖3 不同h情況下Rs對RHRd/RHRm的影響
在無泥漿侵入時,同時考慮圍巖電阻率和地層厚度對幅度差的影響可得到圖4,即RHRd/RHRm隨圍巖電阻率Rs和地層厚度h變化的三維圖形。其中,參數(shù)設(shè)置與圖2相同。圖4中,RHRd/RHRm=1的部分是高分辨率深、中感應(yīng)測井曲線重合的情況,此時基本上不存在幅度差,RHRd/RHRm>1的部分表示正幅度差,RHRd/RHRm<1的部分表示負(fù)幅度差。
圖4 RHRd/RHRm隨Rs和h變化的三維圖形
由以上分析,圍巖電阻率和地層厚度對高分辨率感應(yīng)測井曲線的幅度差都會產(chǎn)生影響。一般情況下,高電阻率圍巖會產(chǎn)生正幅度差,低電阻率圍巖會產(chǎn)生負(fù)幅度差,在一定范圍內(nèi),圍巖電阻率與目的層電阻率之間差別越大,所產(chǎn)生的幅度差越大。地層越厚,圍巖的影響越小,地層越薄,圍巖的影響越大,當(dāng)?shù)貙雍穸仍?~5m范圍內(nèi)時,需要對測井資料進(jìn)行校正。
圖5是無泥漿侵入(非滲透層)時,在不同泥漿電阻率Rm情況下RHRd/RHRm隨井眼半徑dm的變化曲線。其中,地層厚度h為5m,原狀地層電阻率Rt和圍巖電阻率Rs均為100Ω·m。
圖5 不同Rm情況下dm對RHRd/RHRm的影響
從圖5可見,當(dāng)泥漿電阻率很低(Rm/Rt=0.01)時,RHRd/RHRm隨著井眼半徑dm的變化曲線上出現(xiàn)了一個極大值。由于高分辨率中感應(yīng)受井眼影響比深感應(yīng)大,當(dāng)dm較小時,低電阻率泥漿使得中感應(yīng)值比深感應(yīng)值小,隨著dm的增大,中感應(yīng)與深感應(yīng)之間的差距會變大;當(dāng)dm增加到一定的程度,低電阻率泥漿也會對深感應(yīng)造成較大的影響,致使中感應(yīng)與深感應(yīng)之間的差別又逐漸變小,因此會出現(xiàn)極大值點(diǎn)。低電阻率泥漿(Rm/Rt<1)會使RHRd/RHRm>1,高分辨率深、中感應(yīng)測井曲線會出現(xiàn)正幅度差;當(dāng)Rm/Rt≥1時,RHRd/RHRm基本上不隨dm的變化而變化,高分辨率深、中感應(yīng)測井曲線基本上是重合的。
圖6是無泥漿侵入(非滲透層)時,在不同的井眼半徑dm情況下,RHRd/RHRm隨泥漿電阻率Rm的變化曲線。其中,參數(shù)設(shè)置與圖5相同。當(dāng)井眼半徑dm=0.1m 時,RHRd/RHRm基本上不隨泥漿電阻率Rm的變化而變化,這是由于在比較小的井眼半徑范圍內(nèi),井眼對深感應(yīng)和中感應(yīng)的影響相近造成的,這時高分辨率深、中感應(yīng)基本上不存在幅度差,2條曲線基本重合。當(dāng)Rm/Rt<1且dm≥0.2m時,RHRd/RHRm隨泥漿電阻率Rm逐漸變小,這是由于中感應(yīng)受井眼影響比深感應(yīng)大的原因,低電阻率泥漿使得中感應(yīng)的值比深感應(yīng)的值小,因此在這個范圍內(nèi)使得RHRd/RHRm>1,隨著泥漿電阻率Rm的逐漸增大,中感應(yīng)的值會逐漸增大,因此在這個范圍內(nèi)RHRd/RHRm會隨泥漿電阻率Rm逐漸變小。當(dāng)Rm/Rt>1且dm≥0.2m時,RHRd/RHRm基本上不隨泥漿電阻率Rm的變化而變化,此時井眼對高分辨率深、中感應(yīng)造成的影響基本是一樣的,2條曲線基本重合。
圖6 不同dm情況下Rm對RHRd/RHRm的影響
在無侵入模型(非滲透層)下,同時考慮泥漿電阻率和井眼半徑對幅度差的影響即可得到圖7,即RHRd/RHRm隨泥漿電阻率Rm和井眼半徑dm變化的三維圖形。其中,參數(shù)設(shè)置與圖5相同。圖6中,RHRd/RHRm=1表示高分辨率深、中感應(yīng)測井曲線重合,2曲線之間基本上不存在幅度差,RHRd/RHRm>1的部分表示深、中感應(yīng)測井曲線之間有正幅度差。
由以上分析可見,泥漿電阻率和井眼半徑都會對高分辨率感應(yīng)測井曲線的幅度差產(chǎn)生影響。一般情況下,泥漿電阻率和井眼半徑在一定范圍內(nèi)變化時,會使得高分辨率感應(yīng)測井曲線產(chǎn)生正幅度差,泥漿電阻率越低,井眼半徑越大,產(chǎn)生的幅度差越大。當(dāng)井眼半徑大于0.2m,且Rm/Rt<0.5時,需要對測井資料進(jìn)行校正。
圖7 RHRd/RHRm隨Rm和dm變化的三維圖形
圖8是在有泥漿侵入(滲透層)時,在不同的侵入深度di情況下,RHRd/RHRm隨Ri的變化曲線。其中,井眼半徑dm為0.1m,原狀地層電阻率Rt、泥漿電阻率Rm和圍巖電阻率Rs均為100Ω·m。
從圖8可見,當(dāng)di=0時,也就是不存在侵入時,RHRd/RHRm基本上是一條平直線,此時不存在幅度差。當(dāng)h=1且Ri/Rt<1時,RHRd/RHRm隨著Ri逐漸變小,并且高分辨率深、中感存在正幅度差,因?yàn)橹懈袘?yīng)比深感應(yīng)受侵入帶影響要大,所以在這個范圍內(nèi)深感應(yīng)的值要比中感應(yīng)的值大,因此會出現(xiàn)正幅度差;隨著Ri的增大,使中感應(yīng)值在逐漸增大,因而RHRd/RHRm會逐漸變小,2條曲線之間的幅度差也在逐漸減小。當(dāng)h=1且Ri/Rt>1時,RHRd/RHRm幾乎不再隨Ri而變化,高分辨率深、中感應(yīng)測井曲線基本上沒有幅度差,侵入帶不再對幅度差產(chǎn)生影響。當(dāng)h=10m時,RHRd/RHRm幾乎不隨Ri的變化而變化,高分辨率深、中感應(yīng)測井曲線基本上是重合的,因?yàn)榇藭r地層厚度較厚,形成的侵入帶也就相對較大,侵入帶對深、中感應(yīng)的影響基本上是一樣的,因此2條曲線基本上不會產(chǎn)生幅度差。
圖8 不同di情況下Ri對RHRd/RHRm的影響
圖9是有泥漿侵入(滲透層)時,在不同的侵入帶電阻率Ri情況下,RHRd/RHRm隨侵入深度di的變化曲線。從圖9中可以看出,當(dāng)?shù)貙雍穸萮=10m時,RHRd/RHRm基本上不隨di的變化而變化的。當(dāng)?shù)貙雍穸萮=1m且Ri/Rt<1時,RHRd/RHRm的變化曲線上會出現(xiàn)一個極大值點(diǎn),高分辨率深、中感應(yīng)測井曲線會出現(xiàn)正幅度差。這是由于侵入帶對中感應(yīng)影響要比深感應(yīng)大,導(dǎo)致中感應(yīng)的值比深感應(yīng)小,隨著侵入深度di在一定范圍內(nèi)增大,侵入帶對中感應(yīng)的影響會越來越明顯,因此RHRd/RHRm變化曲線在一定范圍內(nèi)是上升的;但是當(dāng)侵入深度超過一定值時,侵入帶對深感應(yīng)的影響也會增大,因此RHRd/RHRm變化曲線會隨著di的增大而減??;當(dāng)di增大到一定程度時,RHRd/RHRm不再隨di變化而變化,此時侵入帶對深、中感應(yīng)造成的影響基本上是一樣的。隨著Ri/Rt的增大,RHRd/RHRm曲線上的極大值點(diǎn)會變得不明顯。當(dāng)?shù)貙雍穸萮=1m且Ri/Rt>1,RHRd/RHRm基本上不隨di的變化而變化的,此時高分辨率深、中感應(yīng)測井曲線基本上是重合的。
圖9 不同Ri情況下di對RHRd/RHRm的影響
圖10 RHRd/RHRm隨Ri和di變化的三維圖形
在有泥漿侵入時,同時考慮侵入帶電阻率和侵入半徑對幅度差的影響即可得到圖10,即RHRd/RHRm隨侵入帶電阻率Ri和侵入半徑di變化的三維圖形。其中,圖10(a)的地層厚度h為1m,圖10(b)的地層厚度h為10m,其他參數(shù)與圖8相同。圖10(a)中,RHRd/RHRm=1表示高分辨率深、中感測井曲線重合,此時2條曲線基本上不存在幅度差;RHRd/RHRm>1的部分表示正幅度差,RHRd/RHRm<1的部分表示負(fù)幅度差。在圖10(b)中雖然有向下凹陷的部分,但是RHRd/RHRm只在1.04~1.05之間變化,相比圖10(a),變化幅度非常小,所以圖10(b)顯示的基本上是平直的,在這個范圍內(nèi)高分辨率深、中感應(yīng)2條曲線基本上是重合的。
由以上分析可以看出,侵入帶電阻率和侵入深度都會對高分辨率感應(yīng)測井的幅度差產(chǎn)生影響。一般情況下,侵入帶電阻率和侵入深度在一定范圍內(nèi)變化時,會使得高分辨率感應(yīng)測井曲線產(chǎn)生正幅度差。侵入帶電阻率越低,產(chǎn)生的幅度差越大。當(dāng)?shù)貙虞^薄,侵入深度大于0.2m,且Ri/Rt<0.02時,需要對測井資料進(jìn)行校正。
(1)圍巖、井眼、侵入等因素都會對高分辨率感應(yīng)測井曲線的幅度差產(chǎn)生影響。在一定情況下,高電阻率圍巖會產(chǎn)生正幅度差,低電阻率圍巖會產(chǎn)生負(fù)幅度差,在一定范圍內(nèi),圍巖電阻率與目的層電阻率之間差別越大,所產(chǎn)生的幅度差越大,且地層越厚,圍巖的影響越小,地層越薄,圍巖的影響越大。高分辨率感應(yīng)測井隨侵入和井眼條件變化時,在一定程度上會產(chǎn)生正幅度差。侵入帶電阻率越低,產(chǎn)生的幅度差越大;泥漿電阻率越低,井眼半徑越大,產(chǎn)生的幅度差越大。
(2)高分辨率感應(yīng)測井不僅僅在滲透層會產(chǎn)生幅度差,在非滲透層(無侵入模型)也會出現(xiàn)曲線分離現(xiàn)象,通常稱為“雙軌”現(xiàn)象?!半p軌”的產(chǎn)生有可能是儀器故障產(chǎn)生的,也有可能是圍巖層厚和井眼等環(huán)境因素造成的。當(dāng)目的層較薄、井眼擴(kuò)徑嚴(yán)重、泥漿電阻率非常低、圍巖電阻率與目的層電阻率相差較大等情況下都有可能會造成“雙軌”。由環(huán)境因素影響造成的“雙軌”可以通過感應(yīng)測井反演或校正等方法加以克服。
(3)測井響應(yīng)是地層中諸多因素的綜合體現(xiàn),不同的地層條件和井眼環(huán)境會產(chǎn)生不同的幅度差。利用高分辨率感應(yīng)測井深、淺電阻率曲線的正、負(fù)幅度差來評判地層特性時,需要綜合考慮圍巖、井眼、泥漿侵入等各方面因素變化。
[1]Paul L Sinclair,Robert W Strickland.Coil Array for a High Resolution Logging Tool and Method of Logging in Earth Formation:United States Patent,5065099[P].1991-11-12.
[2]Carlos Silva,David Spooner.High Resolution Induction Logging Comparison with Conventional Induction as Used in Thin Sands in the Gulf Coast Region[J].SPWLA 32nd Annual Logging Symposium,June 16-19,1991.
[3]Paul L Sinclair,Robert W Strickland.Method and Apparatus for Determining the Conductivity of Subsurface Earth Formations by Filtering and Summing Inphase and Quadrature Conductivity Signals with Correction Values:United States Patent,5146167[P].1992-09-08.
[4]柯式鎮(zhèn),馮啟寧,孫艷茹,等.雙側(cè)向測井幅度差的影響因素與“雙軌”現(xiàn)象[J].測井技術(shù),1999,23(2):123-126.
[5]劉國民,阮民,張魯惠.利用雙感應(yīng)曲線正負(fù)幅度差識別水淹層[J].測井技術(shù),2003,27(4):309-312,348.
[6]鄧少貴,范宜仁,謝關(guān)寶,等.泥漿侵入地層雙感應(yīng)測井曲線正負(fù)差異特性分析[J].測井技術(shù),2004,28(6):496-498.
[7]陳桂波,汪宏年,姚敬金,等.用積分方程法模擬各向異性地層中三維電性異常體的電磁響應(yīng)[J].地球物理學(xué)報(bào),2009,52(8):2174-2181.
[8]魏寶君,田坤,張旭,等.用并矢Green函數(shù)的矢量本征函數(shù)展開式評價(jià)偏心對隨鉆電磁波電阻率測井響應(yīng)的影響[J].中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,2010,34(5):57-62.
[9]張庚驥.交流電測井響應(yīng)的積分方程解法[J].中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,2001,25(1):80-86.
[10]魏寶君,王甜甜,王穎.用磁流源并矢Green函數(shù)的遞推矩陣方法計(jì)算層狀各向異性地層中多分量感應(yīng)測井響應(yīng)[J].地球物理學(xué)報(bào),2009,52(11):2920-2928.