自然電位測井?dāng)?shù)值計(jì)算中的激勵源分布

范曉敏，李舟波

吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院，長春 130026

自然電位測井?dāng)?shù)值計(jì)算中的激勵源分布

范曉敏，李舟波

吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院，長春 130026

為實(shí)現(xiàn)自然電位測井的高效正演計(jì)算，以定量研究各種因素影響，針對自然電位的激勵源——偶電層的分布進(jìn)行了研究。通過對自然電位產(chǎn)生原因及偶電層對自然電位貢獻(xiàn)的分析，認(rèn)為只有泥巖井壁、侵入帶中砂泥巖界面和侵入帶與原狀地層分界面處的偶電層對自然電位測井值有貢獻(xiàn)，而研究砂巖層處的自然電位異常，則僅考慮由侵入帶與原狀地層界面處的偶電層即可。通過直流電場與靜電場的類比，推出了偶極源的面電流密度表達(dá)式，并通過計(jì)算研究了均勻介質(zhì)中侵入深度的影響。結(jié)果顯示，自然電位曲線形狀隨侵入深度增加而變得平緩，異常幅度減小。

自然電位測井；偶電層；電化學(xué)電動勢；侵入深度影響

自然電位測井是實(shí)際生產(chǎn)中使用最多的常規(guī)測井方法之一，目前在儲層定性識別和泥質(zhì)含量估計(jì)方面起重要作用，隨著數(shù)值模擬研究的深入，自然電位還將為侵入帶研究和儲層參數(shù)計(jì)算提供有用信息［1］。同時(shí)，儲層定量評價(jià)的深入發(fā)展和充分利用所有測井信息解決復(fù)雜儲層評價(jià)工作也要求自然電位測井的高精度數(shù)值模擬。自然電場所滿足的方程與直流電阻率法相同，但自然電位測井中場源的處

理較為復(fù)雜，這也是造成其數(shù)值模擬工作遠(yuǎn)遠(yuǎn)少于電阻率測井的一個(gè)原因。一般自然電位的場源可以歸結(jié)為偶電層來考慮，但對于鉆井穿過砂泥巖剖面后偶電層的分布，文獻(xiàn)中有不同的假設(shè)：一種認(rèn)為只是在有濃度差的位置才有偶電層；另一種認(rèn)為除了濃度差，一些巖性界面，包括泥巖表面也存在偶電層。對于后一種假設(shè)，要處理砂巖和泥巖的整個(gè)接觸面。Guyod［2］給出了3種等效模型，但沒有相應(yīng)的分析和解釋。自然電位測井的數(shù)值模擬可采用有限差分法［3－4］、有限元法［5］和模式匹配法［6］等，這些方法都需要在模型中給定偶電層的位置和電動勢數(shù)值或激勵源的電流密度值。筆者通過對影響自然電位數(shù)值的雙電層分布進(jìn)行分析，認(rèn)為只有存在濃度差的偶電層才能影響自然電位測井的數(shù)值，而其他偶電層在計(jì)算中不必考慮。

1 偶電層的產(chǎn)生和分布

關(guān)于石油測井中自然電位產(chǎn)生的原因，目前公認(rèn)的理論認(rèn)為是電化學(xué)電動勢引起的，包括砂巖中由于濃度差引起的擴(kuò)散電位和泥巖中的薄膜電位。擴(kuò)散電位是在砂巖兩側(cè)溶液存在濃度差時(shí)，離子會從濃度高的一側(cè)向濃度低的一側(cè)遷移，由于2種離子的遷移速度不同，造成正、負(fù)離子濃度差異，于是在砂巖兩側(cè)產(chǎn)生電位差，亦稱電動勢。薄膜電勢稍復(fù)雜些，黏土礦物因低價(jià)離子置換高價(jià)離子等原因而帶負(fù)電，孔隙溶液中陽離子吸附在黏土礦物表面，從而維持電荷平衡，因此這些陽離子又稱平衡離子。由于這些陽離子被吸附在黏土礦物的周圍，導(dǎo)致礦物顆粒之間的空隙中溶液的陽離子濃度大于陰離子濃度，形成所謂“擴(kuò)散層”，泥巖中孔隙尺寸小于“擴(kuò)散層”厚度。泥巖兩側(cè)溶液出現(xiàn)濃度差時(shí)，離子會通過泥巖從濃度高的一側(cè)向濃度低的一側(cè)擴(kuò)散移動，而發(fā)生移動的陽離子數(shù)多于陰離子，結(jié)果造成泥巖兩側(cè)正、負(fù)離子濃度差異，于是產(chǎn)生電位差［7］。因?yàn)樵诖诉^程中既涉及濃度差引起的擴(kuò)散，又與黏土吸附陽離子有關(guān)，也有稱其為擴(kuò)散吸附電動勢［8］。早期研究發(fā)現(xiàn)，空井（無泥漿）中泥巖層相對于砂巖層具有更高的電位［9］，因此有人以此來解釋泥巖處自然電位測井值較高的現(xiàn)象。在偶電層分布問題上，目前有2種觀點(diǎn)：一種認(rèn)為只有在濃度差存在的區(qū)域才發(fā)生偶電層（圖1a）；另一種則認(rèn)為除了濃度差，在砂泥巖界面處也有偶電層存在（圖1b）。2種觀點(diǎn)的差別就在于，自然電位異常中是否包括砂泥巖地層界面處偶電層的貢獻(xiàn)。

2 有效偶電層及其電位差

對于自然電位測井，有效偶電層是指對自然電位測井值有貢獻(xiàn)的那些偶電層。自然電位測井測量的是井中測量電極與地面參考電極之間的電位差，它反映井下測量電極的電位隨深度的變化，得到的結(jié)果是自然電位測井曲線。電位隨深度變化的原因是在某種自然電場作用下，鉆井泥漿中形成了電流，該電流的流動造成井孔中縱向電位變化；因此，需要找出能夠?qū)е戮须娏髁鲃拥碾妱觿荨?/p>

圖1 兩種典型的偶電層分布模型Fig.1 Two typical models of electric double layers

為分析砂巖和泥巖界面處的偶電層對自然電位測井值的影響，首先考慮沒有濃度差的情況，即砂巖孔隙和井孔中都是相同的溶液。這時(shí)砂泥巖層界面處的偶電層與泥巖井壁處的偶電層大小相等，方向相反，互相抵消了；由于沒有濃度差，擴(kuò)散電動勢也不存在了，于是井中沒有自然電流流動，無法測到自然電位異常。泥巖表面圍著一個(gè)閉合的偶電層面，因?yàn)樵撆茧妼觾蓚?cè)的電位差恒定，它對外并不產(chǎn)生電場，因而對自然電位異常沒有貢獻(xiàn)。

圖2 有效偶電層和電流分布示意圖Fig.2 Distribution of effective electric double layer and current

其次考慮地下溶液中存在濃度差（比如泥漿濾液與地層水的濃度不同）的情況，這時(shí)2種濃度的溶液之間會出現(xiàn)不同符號的電荷聚集而具有電位差。只要濃度差得到維持，電位差就會存在。因?yàn)殡姾删奂碗娢徊钪怀霈F(xiàn)在濃度變化的區(qū)域，并且正、負(fù)電荷所帶電量相等，所以可認(rèn)定在水中離子濃度有差異的地方發(fā)生偶電層（圖2）。在砂泥巖剖面上，滲透性砂巖中偶電層發(fā)生在沖洗帶與原狀地層之間的過渡帶（圖2中CE），不具滲透性的泥巖層則發(fā)生在井壁（圖2中AB）處。前者偶電層兩側(cè)的電位差是擴(kuò)散電動勢，后者為薄膜電動勢，合在一起就是電化學(xué)電動勢。這時(shí)泥巖與砂巖界面處相對于上面的無濃度差情形沒有發(fā)生任何變化，電位差也沒有變化；但泥巖井壁處因?yàn)槿芤海酀{）濃度變化，界面偶電層的電位差發(fā)生變化，于是整個(gè)泥巖表面的偶電層不再是均勻的了，對外產(chǎn)生電場，對于自然電位異常有貢獻(xiàn)。為研究問題方便，可以把泥巖井壁的界面電位差看成兩部分的和：一部分是由于濃度差造成的電位差，即薄膜電動勢；另一部分是沒有濃度差情況下的電位差。這樣，后一部分可以與砂泥巖地層界面處的偶電層相抵消，對外界電場起作用的就只有泥巖井壁處的薄膜電動勢：

式中：Em為薄膜電動勢，mV；Km為薄膜電動勢系數(shù)，mV，與溫度和離子成分有關(guān)；cw和cmf分別為地層水和泥漿濾液的離子質(zhì)量分?jǐn)?shù)，10－6。

砂巖和泥巖層界面在侵入帶部分與原狀地層部分不同，因?yàn)樵撉秩霂Р糠稚皫r孔隙中溶液往往不同于地層水。本文不區(qū)分沖洗帶和侵入帶，或者說，侵入帶就是沖洗帶。這對某些地層會給計(jì)算帶來一些誤差，但因過渡帶較薄，誤差不會很大。一般情況下，泥漿濾液的電阻率略低于泥漿，其礦化度因泥漿的不同而不同。這是否意味著侵入帶這段地層界面（圖2中BC）具有不確定性呢？其實(shí)泥漿是由泥漿濾液與呈懸浮狀態(tài)的固體顆粒組成的懸濁液，真正在電化學(xué)作用中起作用的溶液只是泥漿濾液，就像砂巖中起作用的是地層水一樣。所以，泥漿濾液與泥漿是相同的溶液，BC段界面的電位差與井壁泥巖是一樣的，即Em。

由于泥漿滲入地層而在井壁處產(chǎn)生過濾電動勢或稱動電電動勢一般較小，且在泥巖井壁和泥餅（砂巖井壁）處都存在，至少大部分可相互抵消，故對自然電位測井值貢獻(xiàn)很小，本文不考慮該電動勢。這樣，自然電位異常由薄膜電動勢和擴(kuò)散電動勢組成，分布在砂巖層中侵入帶和原狀地層之間的擴(kuò)散電動勢為

式中：Ed為擴(kuò)散電動勢，mV；Kd為擴(kuò)散電動勢系數(shù)，mV，與溫度和離子成分有關(guān)。

3 偶電層分布的簡化

現(xiàn)在已經(jīng)明確，對自然電位測井值有貢獻(xiàn)的偶電層分布在泥巖井壁和侵入帶與原狀地層的分界面（圖3）。圖中以灰色標(biāo)出的是有效偶電層，其深淺差異表明電位差不同。計(jì)算時(shí)需要處理整個(gè)井筒，其中砂巖層的侵入帶與原狀地層界面還要向地層內(nèi)縮進(jìn)一段，工作量較大。下面研究在等效前提下的進(jìn)一步簡化。首先分析無砂巖（無限厚泥巖）情況下的自然電位。在同一油水系統(tǒng)中（地層水成分不變），這時(shí)即使井徑隨深度變化，由于偶電層的電動勢為恒定值，井內(nèi)自然電位是定值，自然電位（VSP）曲線是直線。然后再考慮向地層剖面植入砂巖層，于是出現(xiàn)了自然電位異常，砂巖VSP曲線的異常是在原直線基礎(chǔ)上向左（泥漿濾液礦化度低于地層水情況下）偏移的幅度。也就是說，實(shí)測自然電位值減去厚層泥巖處的讀數(shù)（泥巖基線），得到的就是自然電位異常，而這個(gè)異常恰好是砂巖層侵入帶與原狀地層分界面處的偶電層產(chǎn)生的電場造成的。如果在圖3a中偶電層的位置整體減去薄膜電動勢Em，結(jié)果井壁和侵入帶中砂泥巖地層界面偶電層被減掉，只剩下砂巖層侵入帶與原狀地層分界面的偶電層（圖3b），但該偶電層的電位差變成

式中：Eec為電化學(xué)電動勢，mV；Kec是電化學(xué)電動勢系數(shù)，mV。對于NaCl溶液，Em為正值，而Ed為負(fù)值，結(jié)果Kec為負(fù)值（18℃時(shí)為－69.6mV）。必須說明，上述簡化是在地層水礦化度恒定的條件下進(jìn)行的，如果地層水發(fā)生變化，泥巖基線也會移動，不能這樣處理。

圖3 偶電層模型化簡Fig.3 Simplication of electric double layer model

4 電偶極層面的電流密度和曲線計(jì)算

前面已經(jīng)給出了偶極源的電動勢，在一些數(shù)值計(jì)算方法中，直接使用該電位差即可，而有些方法則是用電流源更方便，下面利用直流電場與靜電場的類比，得出偶極層的電流分布。在地下產(chǎn)生自然電位的源是偶電流層，可看作是由2個(gè)貼近的電流面構(gòu)成［10］。它們的面電流密度分別為i和－i（A／ m2），距離恒定為l（m），則偶電流面的單位面積偶極矩為τ＝il，單位A／m。偶極層兩側(cè)電位差為ΔV時(shí)，可得到層矩為τ＝ε·ΔV，其中ε為介電常數(shù)，F(xiàn)／m。通過類比，可知在穩(wěn)恒電流場中層矩τ＝σ· ΔV＝ΔV／ρ，其中σ是介質(zhì)電導(dǎo)率（S／m），ρ是電阻率（Ω·m）。于是得到面電流密度：

在均勻介質(zhì)中，偶電層產(chǎn)生的電位可以用偶電層對觀察點(diǎn)所張開的立體角來計(jì)算［11］，空間任一點(diǎn)電位為

式中，Ω是偶層面對觀測點(diǎn)所張開的立體角，sr。

圖4 不同侵入深度條件下的自然電位測井曲線Fig.4 SP logs for different depths of invasion

利用式（5），可計(jì)算出均勻介質(zhì)條件下不同侵入深度的自然電位曲線，反映侵入深度對VSP曲線的影響。模型為泥巖中間夾一層厚度為2m的砂巖，即111～113m，Eec＝－69.6mV，并假設(shè)泥漿、砂巖和泥巖的電阻率相同。圖4中的VSP1、VSP3和VSP3分別為侵入帶半徑等于0.1m、0.3m和0.5m條件下計(jì)算的自然電位測井曲線?？梢钥闯觯S著侵入深度增加，自然電位測井曲線異常幅度減小，曲線變得平緩，界面反映趨于模糊。

5 結(jié)論

1）真正對自然電位測井值有影響的是分布于泥巖井壁、侵入帶中砂泥巖界面和侵入帶與原狀地層分界面處的偶電層，前2個(gè)偶電層的電位差為薄膜電動勢，最后一個(gè)為擴(kuò)散電動勢。

2）在地層水礦化度不變的條件下，如果僅研究砂巖層的自然電位曲線異常，只要計(jì)算侵入帶與原狀地層分界面處的偶電層的貢獻(xiàn)即可，其電位差在數(shù)值上等于其電化學(xué)電動勢。

3）影響自然電位曲線的偶電層實(shí)際上是2個(gè)符號相反的面電流層面，它們的面電流密度為i＝Eec／（lp）。

4）僅考慮侵入帶與原狀地層分界面偶電層的計(jì)算方法可以正確計(jì)算出自然電位曲線異常；泥漿侵入越深，自然電位曲線異常幅度越小，曲線越平緩。

（References）：

［1］ 張沖，毛志強(qiáng)，肖亮，等.利用測井資料計(jì)算陽離子交換量Qv的方法對比［J］.吉林大學(xué)學(xué)報(bào)：地球科學(xué)版，2010，40（5）：1194－1198.

Zhang Chong，Mao Zhi－qiang，Xiao Liang，et al.Comparation of Calculating Methods of Cation Exchange Capacity Qv Using Logging Data［J］.Journal of Jilin University：Earth Science Edition，2010，40（5）：1194－1198.

［2］ Guyod H.An Investigation of the Factors Affecting the SP in Soft Formations［C］／／Symposium Transaction of SPWLA.Houston：SPWLA，1964.

［3］ 王斌，范曉敏，陳啟艷，等.用有限差分法模擬自然電位測井的侵入帶影響［J］.測井技術(shù)，2008，32（2）：164－168.

Wang Bin，F(xiàn)an Xiao－min，Chen Qi－yan，et al.Modeling Effect of Invaded Zone on SP Log with Finite Difference Method［J］.Well Logging Technology，2008，32（2）：164－168.

［4］ 潘克家，譚永基.復(fù)雜地層中自然電位測井的高效數(shù)值模擬［J］.石油地球物理勘探，2009，44（3）：371－376.

Pan Ke－jia，Tan Yong－ji.An Efficient Algorithm for Simulation of Spontaneous Potential Log in Complex Formations［J］.Oil Geophysical Prospecting，2009，44（3）：371－376.

［5］ 孔峰，譚永基.自然電位測井中的一個(gè)數(shù)值算法［J］.復(fù)旦學(xué)報(bào)：自然科學(xué)版，2010，49（4）：483－490.

Kong Feng，Tan Yong－ji.A Numerical Method for Spontaneous Potential Well－Logging［J］.Journal of Fudan University：Natural Science，2010，49（4）：483－490.

［6］ 汪功禮，張庚驥.數(shù)值計(jì)算的SP－3圖版［J］.測井技術(shù)，1999，23（1）：24－28.

Wang Gong－li，Zhang Geng－ji.SP－3Charts by Numerical Method［J］.Well Logging Technology，1999，23（1）：24－28.

［7］ 李舟波.鉆井地球物理勘探［M］.北京：地質(zhì)出版社，2006.

Li Zhou－bo.Borehole Geophysical Exploration［M］.Beijing：Geological Publishing House，2006.

［8］ 張庚驥.電法測井［M］.北京：石油大學(xué)出版社，1996.

Zhang Geng－ji.Electrical Logging［M］.Beijing：Petroleum University of China Press，1996.

［9］ Dickey P A.Natural Potential in Sedimentary Rocks［C］／／Transactions of AIME.［S.l.］：AIME，1944.

［10］ 薛琴坊.場論［M］.北京：地質(zhì)出版社，1978.

Xue Qin－fang.Field Theory［M］.Beijing：Geological Publishing House，1978.

［11］ Vengelishtien B Y.Study on Oil and Gas well Section by Spontuneous Potential Method［M］.Moscow：Geocenter Press，1966.

Distribution of Source in Numeric Modeling of Spontaneous Potential Log

Fan Xiao－min，Li Zhou－bo

College of GeoExploration Science and Technology，Jilin University，Changchun 130026，China

In order to facilitate efficient modeling of SP－log，and investigate the influences factors quantitatively，investigation on the stimulating source for SP，the electrical double layer（EDL）was performed.Through analyzing of the physical origin of spontaneous potential and the contribution of EDL to SP，it is clear that only the EDL at the borehole wall of shale layer，sand－shale boundary with invaded zone and that between invaded zone and virgin formation would contribute to the value of SP－log.In case only the anomalies of sand layer are concerned，to calculate the contribution of the EDL between invaded zone and virgin formation is sufficient.By analog of direct current electric field and static electric field，expression for the surface current density of EDL is deduced.SP－logs for different depth of invasion were calculated for media with uniform resistivity.The results show that the form of SP－log becomes flatter with the increasing depth of invasion，the amplitude of anomaly become smaller at the same time.

spontaneous potential logging；electric double layer；electrochemical EMF；influence of invasion depth

book=2012,ebook=565

P631.81

A

1671－5888（2012） 04－0895－05

2012－04－18

國家自然科學(xué)基金項(xiàng)目（40874057）；國家“863”計(jì)劃項(xiàng)目（2009AA09A202－01－01）

范曉敏（1962－），男，教授，主要從事地球物理測井方法與應(yīng)用研究，E－mail：fanxm＠jlu.edu.cn。