舒心，柯式鎮(zhèn)，許巍，林小穩(wěn)，張琪

（1.中國石油大學（北京）油氣資源與探測國家重點實驗室，北京 102249；2.中國石油大學（北京）地球物理與信息工程學院，北京 102249）

高分辨率陣列側向測井井眼影響自動校正研究

舒心1，2，柯式鎮(zhèn)1，2，許巍1，2，林小穩(wěn)1，2，張琪1，2

（1.中國石油大學（北京）油氣資源與探測國家重點實驗室，北京 102249；2.中國石油大學（北京）地球物理與信息工程學院，北京 102249）

陣列側向測井具有較高的縱向分辨率，在徑向上可以提供多種探測深度的測量曲線，為地層侵入剖面的描述提供更加詳細的信息，但其測井結果不可避免地會受井眼的影響，特別是淺探測深度的測井結果受井眼影響尤為嚴重；因此，需要對其進行井眼校正。文中以斯倫貝謝公司推出的高分辨率陣列側向測井儀HRLA為例，在分析陣列側向測井井眼影響的基礎上，對HRLA的井眼校正圖版進行數(shù)字化，然后優(yōu)選Akima插值算法實現(xiàn)了圖版曲線高精度插值擬合，利用網(wǎng)絡加權分析實現(xiàn)了井眼影響的自動校正，并編制程序對實際陣列側向測井資料進行自動連續(xù)校正處理。應用結果驗證了對井眼影響規(guī)律的分析。校正結果表明，該校正方法和處理程序提高了井眼校正的精度和效率，具有良好的實用性。

陣列側向測井；井眼影響；Akima插值；自動校正

0　引言

傳統(tǒng)的雙側向測井技術一次只能提供2條徑向探測深度固定的視電阻率曲線，受薄層和高阻圍巖的影響大，不能滿足復雜油氣評價需求。陣列側向測井技術是繼雙側向測井技術后發(fā)展而來的一種新的測井技術，一次下井能提供多條不同探測深度的高分辨率測井曲線，進而可以得到地層電阻率的徑向分布、侵入剖面的二維成像等信息，使得陣列側向測井技術在薄層評價、油水識別、巖性劃分、儲量估算等方面具有獨特優(yōu)勢，尤其適用于復雜地層特別是薄層和深侵入地層的測井評價［1-2］。目前陣列側向的儀器主要有斯倫貝謝公司推出的高分辨率陣列側向測井儀High Resolution Laterolog Array Tool（HRLA）［3］，西方阿特拉斯和殼牌公司聯(lián)合推出的高分辨率側向測井儀High Definition Lateral Log（HDLL）［4］，以及中國石油集團測井公司研制的陣列側向儀器 High Resolution Array LaterologLogging Tool（HAL）［5］。無論使用哪款儀器，陣列側向的測井結果，尤其是淺探測深度，都不可避免地受到井眼（包括鉆井液電阻率、井徑和儀器偏心）的影響，特別是在儀器偏心、大井眼、地層和鉆井液電阻率對比度高的情況下，井眼影響尤其嚴重；因此，對其進行井眼校正是必不可少的環(huán)節(jié)，并且如果能在現(xiàn)場就對測井曲線進行實時快速的井眼校正，將會大大提高工作效率。

國內(nèi)目前對于陣列側向井眼校正大多是根據(jù)所用儀器的類型及具體條件，選擇由儀器制造公司所提供的相應解釋圖版或者配套服務，耗時較長，少見實時自動化校正軟件。國內(nèi)外對于測井資料中井眼的自動校正主要有2種途徑：第1種為自適應井眼校正方法，是利用有限元、有限差分或者模式匹配等算法建立井眼校正模型，然后對測井響應進行快速計算，再通過反演的方法得到井眼參數(shù)，從而實現(xiàn)井眼校正。該方法的校正效果主要取決于測井參數(shù)的選取和正、反演算法的選擇，優(yōu)點是可以根據(jù)不同的環(huán)境條件及時修改模型參數(shù)，靈活度大，缺點是計算量大、校正速度慢。該方法比較適用于定性分析和參考，不太適合對測井曲線進行實時的井眼校正［6-8］。第2種途徑是通過理論研究或實驗結果得到井眼校正的解釋圖版，對圖版進行數(shù)字化離散采樣，選取合適的算法對采樣值進行擬合，然后根據(jù)分段的或者統(tǒng)一的擬合公式編制自動化校正程序實現(xiàn)井眼校正。該方法不涉及測井響應的數(shù)值模擬，算法較為簡單，所以比較適用于實時的、快速的井眼校正。該方法采用何種算法的主要依據(jù)為解釋圖版的類型和精度，利用數(shù)學方法和公式模型進行解釋圖版的高精度擬合是自動化校正的關鍵所在［9-11］。本文采用的是第2種方法，針對斯倫貝謝公司推出的高分辨率陣列側向測井儀，闡述了陣列側向的工作原理，分析了井眼影響，對自動化校正算法進行了研究，實現(xiàn)了陣列側向井眼影響的快速自動校正。

1　儀器電極系結構及工作原理

陣列側向測井儀器是在傳統(tǒng)雙側向測井儀器結構的基礎上發(fā)展起來的，由1個發(fā)射電流電極和多個電位測量電極組成，通過多個聚焦線圈系組合，由不同的電極距組成的電極陣列來進行電位電極或梯度電極測量，從而實現(xiàn)不同徑向深度地層電阻率的探測［12-15］。斯倫貝謝公司生產(chǎn)的高分辨率陣列側向測井儀器（HRLA）可以依次得到6種不同探測深度的測井響應（RHLA0，RHLA1，RHLA2，RHLA3，RHLA4，RHLA5），其中 RHLA0探測深度最淺，主要探測鉆井液電阻率。RHLA1至RHLA5均對地層敏感，其探測深度逐漸變深，是求取地層真實電阻率主要考慮的測井響應。

2　井眼校正圖版

井眼環(huán)境對陣列側向儀器響應的影響程度主要與鉆井液電阻率、井眼大小以及儀器的偏心程度3個關鍵因素有關。斯倫貝謝公司HRLA的井眼校正圖版一共有15張，圖版中包含井徑127.0~558.8 mm、偏心塊分別為0，12.7，38.1 mm時井眼校正情況。其中偏心塊為38.1 mm時，儀器在井徑127.0，152.4 mm井眼幾何空間中無法正常作業(yè)，因此井眼井徑范圍為203.2~ 558.8 mm［16］。

圖1是井眼居中情況下模式1測井曲線RHLA1的井眼校正圖版。其中，橫坐標為視電阻率（RHLA1）與鉆井液電阻率（Rm）的比值，縱坐標為地層電阻率（Rt）與視電阻率（RHLA1）的比值，不同的曲線類型代表不同井徑的校正曲線。

由圖1可以看出，井徑越大，電阻率曲線受井眼影響越大，并且當鉆井液電阻率大于視電阻率時，校正系數(shù)小于1。

本文用圖形數(shù)字化的方法對圖版進行離散采樣，并且用Akima算法擬合得到單一井徑下不同探測深度的校正曲線。圖2是儀器居中時，井徑為203.2 mm 和 406.4 mm的井眼對測井曲線影響，圖中 j為1，2，3，4，5。從圖2a可以看出，淺探測深度的RHLA1受井眼影響最大，RHLA2次之，而其他3條探測深度較深的視電阻率RHLA3—RHLA5受井眼的影響很小，幾乎可以忽略。但對比圖2b井徑為406.4 mm時可以看出，每個探測深度均受井眼嚴重影響，并且探測深度越淺，視電阻率受井眼影響越嚴重。

在測井過程中，由于多種原因會引起儀器偏心，這時儀器不再位于井眼中央，而是有了一定的偏離，因此儀器偏心而造成的井眼影響同樣不能忽略。圖3為在203.2 mm井眼中，偏心距為12.7 mm和38.1 mm時的井眼影響。

從圖3可以看出，偏心距為12.7 mm時，淺探測RHLA1和RHLA2受井眼影響較明顯，RHLA4和RHLA5則基本不受影響；偏心距為38.1 mm時和偏心距為12.7 mm時趨勢基本一致。對比圖2a和圖3可以看出，隨著視電阻率與鉆井液電阻率比值的增大，視電阻率逐漸開始偏離地層真電阻率，但偏心距越大，視電阻率偏離地層真電阻率的趨勢變緩，并且逐漸右移。由此可以得出，儀器偏心會導致儀器測量到的視電阻率讀數(shù)降低，偏心越大視電阻率降低越嚴重。

3　自動化校正算法

為了實現(xiàn)計算機自動校正，需要對井眼校正圖版數(shù)字化，圖版數(shù)字化是將圖版中的連續(xù)模擬量轉換成離散的數(shù)字量的過程?？紤]到陣列側向的校正圖版是二維曲線，圖像分辨率高，本文采用了digXY軟件對其進行數(shù)字化，在確定坐標值范圍后，該軟件可以通過鼠標點擊將圖版中的坐標直接采集出來并且按順序保存至文本文檔中。它的不足在于，通過鼠標直接采集的坐標點是人為確定的非規(guī)律性的數(shù)據(jù)點，這樣的數(shù)據(jù)在后期進行自動校正時會產(chǎn)生較大的誤差，因此還需要對采集的坐標點進行數(shù)據(jù)重采樣，將不等間距數(shù)據(jù)轉化為等間距數(shù)據(jù)。通過重采樣能夠確定合適的數(shù)據(jù)步長來保證數(shù)據(jù)的有效性，去除冗余的數(shù)據(jù)。例如：當曲線變化平緩時只需要幾個采樣點就可達到精度要求，此時步長可以適量大一點；而在曲線轉彎處，步長應該比較小。重采樣的效果決定了后期自動校正的準確度。

對曲線進行數(shù)值逼近，即對曲線的高精度擬合，然后建立數(shù)據(jù)的重采樣模型。本文中的數(shù)字化結果要求滿足數(shù)值逼近結果，保持原有的曲線形態(tài)且曲線形態(tài)要光滑，因此適合使用插值法。

數(shù)據(jù)的內(nèi)插有各種方法，如線性內(nèi)插、多項式內(nèi)插、樣條函數(shù)插值等。線性內(nèi)插只顧及附近兩點的影響；多項式內(nèi)插時，低階多項式由于參數(shù)較少，內(nèi)插精度很低，會出現(xiàn)“龍格”現(xiàn)象：因此，又在此基礎上發(fā)展了三次樣條插值法。它避免了多項式的階數(shù)較高時數(shù)值不穩(wěn)定的缺點，因而得到了廣泛的應用。三次樣條插值方法是一種全局的插值方法，計算過程和程序設計較為復雜。另外，常用的克里金、分形等方法因為自動化程度較低也不適用于本文。Akima插值法作為一種局部的分段三次插值方法，和三次樣條函數(shù)一樣考慮了要素導數(shù)值的效應；因而得到的整個插值曲線是光滑的，并且Akima插值法具有最小模、最佳最優(yōu)逼近和收斂的特性，插值曲線不會因為數(shù)據(jù)點數(shù)較少而出現(xiàn)搖擺的現(xiàn)象［15-16］，是一種高效優(yōu)質的插值重采樣方法。鑒于此，本文提出了一種基于Akima插值法的井眼影響的自動化校正算法。

3.1Akima插值法

Akima插值法是在每2個數(shù)據(jù)點之間建立1個三次多項式曲線方程，其規(guī)定除了需要用到這2個實測值外，還需要用2點之外相鄰的2個實測點上的觀測值，以保證整條曲線上具有連續(xù)的一階導數(shù)和曲線的光滑性。設已知數(shù)據(jù)點為（xi，yi）（i=1，2，3，…，n），若在子區(qū)間 ［xk，xk+1］（k=0，1，…，n-2）上的2個斷點處有以下4個條件：

式中：ti為時間。

現(xiàn)唯一地確定一個三次多項式曲線為

由此可以計算該子區(qū)間中的插值點x（x∈［xk，xk+1］）處的函數(shù)近似值y（x）。

式（3）中uk，uk+1分別為兩點導數(shù)，一般情況下為

但在端點處，不能采用式（6）計算，需要在端點以外補足2個點，以保證曲線的原有趨勢。假設端點（xn，yn）和2個相鄰的數(shù)據(jù)點（xn-1，yn-1），（xn-2，yn-2），以及需要補上的2點（xn+1，yn+1），（xn+2，yn+2）都在拋物線y=p0+p1· （x-xn）+p2（x-xn）2（p0，p1，p2為常數(shù)）上，并設xn-2-xn= xn+1-xn=xn-xn-2，則可以得到：

按照式（3）—式（9），就可以編寫程序逐一求得三次多項式式（2）的各項系數(shù)，這樣就可以編制程序對校正點進行插值。

為檢驗Akima插值的質量，根據(jù)算法編制程序，選用三次樣條算法程序作比較，選取了測井模式為RHLA1（偏心距大小12.7 mm）以及RHLA3（儀器居中）的井眼校正圖版中所有井徑共22條曲線的數(shù)字化數(shù)據(jù)。每條曲線隨機抽取1/4的樣本點作為控制點，然后再利用這些控制點用Akima和三次樣條不同插值方法進行插值重采樣得到新的曲線，最后將得到的曲線分別與原始曲線進行相關計算，得到相關系數(shù)（見表1）。

從表1數(shù)據(jù)中可以看出，總體上，Akima插值方法比三次樣條插值法得到的相關系數(shù)大，在數(shù)據(jù)點較少的時候比三次樣條插值方法更接近原始曲線形態(tài)，且不易出現(xiàn)曲線搖擺現(xiàn)象。

3.2自動化校正算法

在進行陣列側向井眼影響的自動校正時，主要分為2個方面：1）根據(jù)15張校正圖版的數(shù)字化結果，利用Akima插值法建立數(shù)據(jù)重采樣模型。首先，設定時間間隔Δt，對Δt內(nèi)的數(shù)據(jù)進行網(wǎng)格化處理，分為m個均勻的網(wǎng)格區(qū)間；然后，將該網(wǎng)格區(qū)間內(nèi)的數(shù)據(jù)點作為第i個網(wǎng)格的數(shù)據(jù)集合；最后，利用Akima插值法計算該數(shù)據(jù)集合內(nèi)的采樣點值。2）利用網(wǎng)絡加權分析法計算井眼影響的校正值。需要根據(jù)實際地層參數(shù)（RHLAj/ Rm、儀器偏心距大小、井徑值），選擇相應校正圖版，從圖版對應的重采樣數(shù)據(jù)模型中直接計算出校正值Rt/ RHLAj，由此得到地層電阻率Rt。如果實際的儀器偏心距或者井徑不滿足圖版現(xiàn)有條件，即RHLAj/Rm落在圖版間或者曲線之間時，則需要利用網(wǎng)絡加權分析計算校正系數(shù)。具體的自動化校正流程和網(wǎng)絡加權分析流程見圖4。

當RHLAj/Rm沒有落在圖版中對應井徑的校正曲線上時，需要通過加權算法得到最后的校正系數(shù)。具體的做法是：假設實際井徑值為a，圖版中與之相鄰的2條井徑值分別為y，y，則權值，然后在y，y所1212對應的2條校正曲線上分別執(zhí)行圖4中的判據(jù)，得到2個校正系數(shù)x1與x2，那么最后的校正系數(shù)x=wx1+ （1-w）x2；當然，當實際的偏心距不是12.7 mm和38.1 mm時，則需要通過2次類似的加權流程得到相應校正系數(shù)。

4　應用實例與效果分析

圖5為同一深度不同探測深度所探測到的5個陣列側向測井響應的電阻率曲線和對應的自動校正曲線。其中，藍色曲線為實際測得的電阻率曲線，紅色曲線為井眼自動校正之后的電阻率曲線，最后一個道為由5個探測深度校正后的電阻率曲線合成得到的地層電阻率曲線。

由圖可以發(fā)現(xiàn)，RHLA1（即探測深度最淺的測井曲線）校正量大于原始數(shù)據(jù)且校正比較明顯，并且隨著探測深度的加深，校正量變小且逐漸與原曲線重合。

應用結果驗證了井眼影響規(guī)律，并且校正后的曲線形態(tài)合理，說明校正結果正確。

5　結論

1）為了保證自動化校正算法的精度，需要對通過數(shù)字化采樣得到的非等間距數(shù)據(jù)點進行重采樣，得到等間距采樣點。Akima分段三次多項式插值作為一種局部插值方法，相比作為全局算法的三次樣條插值法而言，計算過程簡單，計算量小，插值效果更好。

2）陣列側向淺探測測井曲線受井眼影響大，校正量大；隨著探測深度的增加，視電阻率受井眼影響逐漸減小，校正量也隨之減小。校正后各探測深度的測井曲線形態(tài)一致性較好，有效減小了井眼對于不同探測深度測井結果的影響。

3）本文的圖版數(shù)字化方法和自動化校正算法雖然是針對陣列側向井眼校正圖版而言，但該自動化校正方法對于其他常規(guī)測井資料的環(huán)境校正圖版也有一定的適用性。

［1］SMITS J W，BENIMELI D，DUBOURG I，et al.High resolution from a new laterolog with azimuthal imaging［R］.SPE 30584，1995.

［2］范宜仁，蔣建亮，鄧少貴，等.高分辨率陣列側向測井響應數(shù)值模擬［J］.測井技術，2009，33（4）：333-336.

［3］SMITS J W，DUBOURG I，LULING M G，et al.Improved resistivity interpretation utilizing a new array laterolog tool and associated inversion processing［R］.SPE 49328，1998.

［4］ITSKOVICH G B，MEZZATESTA A G，STRACK K M，et al.Highdefinition lateral log-resisitivity device：basic physics and resolution ［C］//Society of Petrophysicists and Well-Log Analysts.SPWLA 39thAnnual Logging Symposium，SPWLA-1998-V.Keystone，Colorado：Society of Petrophysicists and Well-Log Analysts，1998.

［5］賀飛，馬驍，馮琳偉，等.HAL6505陣列側向測井儀［J］.石油科技論壇，2013，32（2）：59-61.

［6］仵杰，謝尉尉，解茜草，等.陣列側向測井儀器的正演響應分析［J］.西安石油大學學報（自然科學版），2008，23（1）：73-76.

［7］譚茂金，高杰，鄒友龍，等.鹽水泥漿條件下定向井雙側向測井環(huán)境校正方法研究［J］.地球物理學報，2012，55（4）：1422-1432.

［8］ZHOU Z，CORLEY B，KHOKHAR R，et al.A new multi laterolog tool with adaptive borehole correction［R］.SPE 114704，2008.

［9］徐波，汪中浩，黎澤剛，等.隨鉆自然伽馬測井資料環(huán)境影響自動校正［J］.石油天然氣學報，2009，31（5）：82-85.

［10］馮琳偉，成志剛.EILog高分辨率雙側向井眼影響及校正效果分析［J］.測井技術，2009，33（6）：530-534.

［11］宋子齊，譚成仟，夏克文，等.視電阻率測井資料環(huán)境影響校正［J］.測井技術，1995，19（4）：250-261.

［12］劉紅岐，夏宏泉，郭璞，等.隨鉆電阻率測井介電效應校正圖版分析及方法［J］.西南石油大學學報（自然科學版），2010，32（2）：58-62.

［13］潘克家，王文娟，湯井田，等.高分辨率陣列側向測井的數(shù)學模型及有限元快速正演［J］.地球物理學報，2013，56（9）：3197-3211.

［14］劉振華，胡啟.陣列側向測井響應的計算及其特征［J］.西安石油學院學報，2002，17（1）：53-57.

［15］AKIM H.A new method of interpolation and smooth curve fitting based onlocalprocedures［J］.JournaloftheACM，1970，17（4）：589-602.

［16］李東升，蘇永濤，滕麗娜.直線度測量數(shù)據(jù)的阿克瑪插值法平滑處理［J］.計量技術，1999（5）：17-18.

（編輯楊會朋）

Auto-correction of borehole effect for high-resolution array laterolog

SHU Xin1，2，KE Shizhen1，2，XU Wei1，2，LIN Xiaowen1，2，ZHANG Qi1，2
（1.State Key Laboratory of Petroleum Resources and Prospecting，China University of Petroleum，Beijing 102249，China；2.College of Geophysics and Information Engineering，China University of Petroleum，Beijing 102249，China）

The array laterolog not only processes superior vertical resolution but also provides various measured curves of different depth on the radial，offering more detailed information for the description of invasion profile within the formation.However，the logging result is inevitably affected by the borehole，especially for the shallow lateral logging.Hence，borehole correction is demanded.This paper focuses on the high-resolution influence caused by the HRLA，an array lateral logging instrument provided by Schlumberger. On the basis of analyzing the result of array laterolog borehole，we digitalized the HRLA borehole correction chart and optimized Akima interpolation algorithm for automatical correction.The auto-correction of borehole effect was achieved by network weight analysis，and the automatic successive correction of practical array lateral logging data was programmed.The analysis of influence rule wasverifiedby the reasonable applicationresults.The correctionmethodandprocesserimprove the accuracy andefficiency ofborehole correction with high practicability.

array laterolog；borehole effect；Akima interpolation；auto-correction

國家科技重大專項課題“復雜儲層測井新方法與探測器研究”（2011ZX05020-009）

TE132.1+4；P631

A

10.6056/dkyqt201604014

2015-10-10；改回日期：2016-05-07。

舒心，女，1990年生，在讀碩士研究生，2013畢業(yè)于西南石油大學測控技術與儀器專業(yè)，現(xiàn)主要從事電測井儀器研究。E-mail：shuxin_cup@163.com。

引用格式：舒心，柯式鎮(zhèn)，許巍，等.高分辨率陣列側向測井井眼影響自動校正研究［J］.斷塊油氣田，2016，23（4）：470-475.

SHU Xin，KE Shizhen，XU Wei，et al.Auto-correction of borehole effect for high-resolution array laterolog［J］.Fault-Block Oil&Gas Field，2016，23（4）：470-475.