王家進(jìn) （中國(guó)石油集團(tuán)鉆井工程技術(shù)研究院，北京102206）

隨鉆遙測(cè)成像技術(shù)[1，2]是隨鉆測(cè)井中高端技術(shù)，用于識(shí)別地層界面變化，進(jìn)行地質(zhì)導(dǎo)向作業(yè)，具有常規(guī)測(cè)井無(wú)法比擬的優(yōu)勢(shì)。然而由于圖像信息量較大，常規(guī)無(wú)線遙測(cè)系統(tǒng)[3，4]信息傳輸速率通常較低，信道容量無(wú)法滿足大信息量圖像的實(shí)時(shí)傳輸。因此，要實(shí)現(xiàn)井底地層測(cè)量參數(shù)的實(shí)時(shí)成像，就需要對(duì)原始圖像信息進(jìn)行處理，即高度壓縮后傳輸。目前國(guó)外已有成熟的產(chǎn)品可以應(yīng)用[5]，但其技術(shù)細(xì)節(jié)仍屬商業(yè)秘密。國(guó)內(nèi)關(guān)于隨鉆遙測(cè)成像技術(shù)的研究，還僅限于方位地層參數(shù)的隨鉆測(cè)量技術(shù)及方法研究[6～8]，由于傳輸信道的局限性，還遠(yuǎn)不能達(dá)到連續(xù)實(shí)時(shí)成像的程度。此外，多數(shù)關(guān)于成像測(cè)量技術(shù)的研究多側(cè)重于成像原理和成像結(jié)果的應(yīng)用[9～13]，可應(yīng)用的成熟技術(shù)多針對(duì)于常規(guī)有纜測(cè)井成像[14～16]，并不適用于無(wú)線隨鉆實(shí)時(shí)傳輸。無(wú)線傳輸信道具有低速特點(diǎn) （通常為每秒數(shù)個(gè)比特），所能實(shí)時(shí)傳輸?shù)男畔⒘糠浅Ｓ邢蓿Ｒ?guī)圖像壓縮方法[14～16]難于湊效。為此，筆者從研究隨鉆遙測(cè)成像測(cè)井的特征入手，通過(guò)提取圖像特征點(diǎn)，將圖像高度壓縮，然后對(duì)經(jīng)由無(wú)線遙測(cè)系統(tǒng)傳到地面的高度壓縮圖像數(shù)據(jù)進(jìn)行圖像重構(gòu)，來(lái)解決隨鉆遙測(cè)成像的實(shí)時(shí)傳輸問(wèn)題。

1 隨鉆遙測(cè)成像測(cè)井圖像特征提取

1.1 隨鉆遙測(cè)成像測(cè)井?dāng)?shù)據(jù)圖像描述方法

地層剖面特性決定了隨鉆井壁或井周成像普遍具有如下特點(diǎn)：連續(xù)、緩變、周向閉合。緩變的圖像可以由控制點(diǎn)描述，以實(shí)現(xiàn)圖像信息的高度壓縮?？刂泣c(diǎn)的選取以能反映緩變圖像邊界為準(zhǔn)，如果圖像突變，可以反映圖像的明顯變化。由于通常同一地層內(nèi)圖像變化幅度不大，僅在地層邊界附近圖像會(huì)有劇烈變化，因此這些控制點(diǎn)多能反映地層界面或?qū)傩缘母淖?，這對(duì)于實(shí)現(xiàn)地質(zhì)導(dǎo)向鉆井作業(yè)具有重要意義。

常規(guī)測(cè)井圖像是由若干測(cè)量參數(shù)數(shù)據(jù)點(diǎn)經(jīng)過(guò)映射著色后繪制而成，直接傳輸圖像并不合適，因?yàn)閳D像信息要比數(shù)據(jù)信息大很多。因此，該次研究在討論圖像壓縮時(shí)，并不是傳統(tǒng)意義上的色彩圖像，而是一個(gè)數(shù)據(jù)圖像，由數(shù)據(jù)構(gòu)成，而不是色點(diǎn)。測(cè)量探頭在某一時(shí)點(diǎn)或深度點(diǎn)完成一周測(cè)量后獲得的數(shù)據(jù)點(diǎn)在極坐標(biāo)系內(nèi)可以構(gòu)成一條閉合曲線s，不同時(shí)點(diǎn)或深度點(diǎn)所測(cè)得的一系列閉合曲線就形成了一個(gè)類似柱面的柱形曲面，如圖1所示。隨著時(shí)間推移或深度變化，圖像是沿時(shí)間軸 （t）或深度軸 （l）不斷延伸。

通常井底儀器只能得到以時(shí)點(diǎn)為縱軸的圖像，需要進(jìn)行時(shí)－深轉(zhuǎn)換，才能得到反映地層空間變化的圖像。在實(shí)際傳輸時(shí)，進(jìn)行圖像特征提取只能基于時(shí)軸坐標(biāo)。在進(jìn)行數(shù)據(jù)圖像特征提取時(shí)，以該時(shí)點(diǎn)的掃描采樣曲線為對(duì)象，抽取表征該曲線的特征點(diǎn)進(jìn)行傳輸，在隨后的傳輸中，選擇某時(shí)點(diǎn)要傳輸?shù)某闃狱c(diǎn)時(shí)，剔除與上次傳輸抽樣點(diǎn)位置和數(shù)值都相近的點(diǎn)。為完整描述數(shù)據(jù)圖像在周向與縱向二緯尺度上的變化，可分2個(gè)步驟進(jìn)行圖像特征的提?。阂皇菍?duì)某一時(shí)間點(diǎn) （對(duì)應(yīng)某一深度）的周向掃描采樣曲線進(jìn)行特征提??；二是對(duì)某一方向不同時(shí)點(diǎn) （對(duì)應(yīng)該方向不同深度）進(jìn)行縱向曲線特征提取。

1.2 周向采樣曲線特征提取

為避免采樣噪聲干擾，在處理數(shù)據(jù)集之前都要經(jīng)過(guò)濾波處理，由于地層屬性的測(cè)量數(shù)值通常是緩變的，即便在邊界處測(cè)量數(shù)值也存在過(guò)渡帶，因此一個(gè)適當(dāng)?shù)牡屯V波器是去除噪聲干擾的好方法。經(jīng)過(guò)濾波平滑后的周向采樣曲線在極坐標(biāo)系內(nèi)表現(xiàn)為一條閉合曲線，對(duì)于均質(zhì)地層這條閉合曲線可以近似一個(gè)圓 （圖2（a））；對(duì)于各向異性地層，則會(huì)出現(xiàn)扭曲 （圖2（b））。在地層界面處，該曲線的形狀與井眼軸線和地層界面夾角有關(guān)，非垂直相交時(shí)，大體上表現(xiàn)為類似梨形的不規(guī)則橢圓形狀 （圖2（c））。因此該類曲線特征的提取以能夠反映曲線曲率的突變?yōu)橹?。如圖2（c）所示，可取曲率變化梯度拐點(diǎn)（圖2中1～4）為特征點(diǎn)，而且如果曲線具有對(duì)稱性，可以只取一半曲線的特征點(diǎn)來(lái)描述該曲線，進(jìn)一步壓縮曲線信息。為防止曲線抖動(dòng)增加曲線描述的復(fù)雜性，可以設(shè)置一個(gè)閾值，僅將曲率變化量連續(xù)（一定區(qū)間）超過(guò)閾值的拐點(diǎn)標(biāo)記為特征點(diǎn)。事實(shí)上，可以設(shè)計(jì)一個(gè)合理的積分－微分濾波器，該濾波器可以過(guò)濾小的抖動(dòng)、保留大的起伏，濾波后超過(guò)一定數(shù)值的曲線極點(diǎn)可以作為曲線特征點(diǎn)的提取依據(jù)。

圖1 隨鉆測(cè)井?dāng)?shù)據(jù)圖像

圖2 周向掃描采樣曲線

1.3 縱向采樣曲線特征提取

與常規(guī)測(cè)井曲線一樣，縱向曲線的變化特征與地層沿井眼軸向的層序變化有關(guān)，在均質(zhì)地層內(nèi)表現(xiàn)直線，非均質(zhì)地層內(nèi)存在波動(dòng)，地層界面處有起伏 （圖3）。因此縱向曲線的特征提取以能夠識(shí)別曲線起伏為主。同樣，為防止曲線抖動(dòng)增加曲線描述的復(fù)雜性，也可以設(shè)置一個(gè)閾值，僅將數(shù)值變化量連續(xù) （一定區(qū)間）超過(guò)閾值的拐點(diǎn)標(biāo)記為特征點(diǎn)。與提取周向采樣曲線特征點(diǎn)的方法類似，也可以用濾波的方法提取特征點(diǎn)。如圖3所示，所提取的特征點(diǎn)1～4用于描述曲線起伏形態(tài)。

1.4 特征點(diǎn)合并

圖3 縱向采樣曲線

經(jīng)過(guò)上述2個(gè)步驟，提取的周向曲線與縱向曲線特征點(diǎn)相比原始測(cè)量數(shù)據(jù)已經(jīng)大幅降低，但這2個(gè)途徑所得到特征點(diǎn)中可能包含重復(fù)點(diǎn)或位置相近且數(shù)值相近的點(diǎn)，可以將這些點(diǎn)合并為一個(gè)點(diǎn)，進(jìn)一步壓縮數(shù)據(jù)。基于地層測(cè)量數(shù)據(jù)的緩變特征，經(jīng)過(guò)上述步驟處理后的數(shù)據(jù)被大幅度壓縮，可以經(jīng)由數(shù)據(jù)傳輸速率較低的隨鉆遙測(cè)系統(tǒng)傳輸通道上傳到地面上來(lái)。

2 成像數(shù)據(jù)的圖像重構(gòu)

如上所述，井下所測(cè)量的地層信息被高度壓縮，這個(gè)過(guò)程損失許多地層細(xì)微信息，屬于有損壓縮，因此在地面接收到這些壓縮數(shù)據(jù)后，需要對(duì)地層信息進(jìn)行重構(gòu)，這個(gè)過(guò)程是圖像復(fù)原過(guò)程。與常規(guī)圖像重構(gòu)過(guò)程[17～20]有所不同，所接收到的數(shù)據(jù)為描述地層的特征點(diǎn)，這些點(diǎn)表現(xiàn)為一系列稀落的不規(guī)則散點(diǎn)，在進(jìn)行成像著色之前，需要進(jìn)行網(wǎng)格化插值處理。網(wǎng)格化插值方法有很多種，如Shepard法、反距離平均法、線性插值三角網(wǎng)格法、Kriging法等[21～25]。這些插值方法各有特點(diǎn)，該次研究不討論這些插值方法的細(xì)節(jié)，只討論圖像重構(gòu)過(guò)程。在利用得到的不規(guī)則散點(diǎn)進(jìn)行圖像重構(gòu)時(shí)，采用了三角網(wǎng)/線形插值法，因?yàn)樵摬逯捣椒軌驀?yán)密遵守特征數(shù)據(jù)點(diǎn)約束，而且計(jì)算簡(jiǎn)單，容易實(shí)現(xiàn)。由于特征數(shù)據(jù)點(diǎn)較稀少，嚴(yán)格尊重?cái)?shù)據(jù)點(diǎn)的約束對(duì)于復(fù)原原始圖像是十分必要的。另外，該插值方法便于實(shí)現(xiàn)局部插值，而不必對(duì)整體數(shù)據(jù)進(jìn)行復(fù)雜計(jì)算，這對(duì)于實(shí)時(shí)圖像處理而言具有重要意義。唯一不足的是得到的圖像在插值點(diǎn)處不光滑，不過(guò)該現(xiàn)象可以通過(guò)一個(gè)高斯平滑濾波處理來(lái)消除。

需要注意的是，曲面是周向閉合的，而在進(jìn)行插值處理時(shí)是在展開曲面上進(jìn)行的，即以方向（角）為橫坐標(biāo)、時(shí)間或深度為縱坐標(biāo)的二維空間上進(jìn)行。插值處理后的圖像不能夠保證在周向邊界連續(xù)，為避免該現(xiàn)象，根據(jù)特征數(shù)據(jù)點(diǎn)在圓周上具有周期性特征，將數(shù)據(jù)點(diǎn)沿方向坐標(biāo) （周向）進(jìn)行周期延拓，比如分別向左、右延拓一個(gè)周期 （圖4中－360～0°，360～720°區(qū)域），然后對(duì)延拓后數(shù)據(jù)點(diǎn)集進(jìn)行插值處理，并限定在延拓前的周期 （圖4中0～360°區(qū)域）內(nèi)的插值結(jié)果為最終結(jié)果，則可以確保該周期內(nèi)的插值結(jié)果在橫向邊界上的連續(xù)性。為了減少計(jì)算量，可不對(duì)延拓周期進(jìn)行網(wǎng)格化計(jì)算，僅利用該區(qū)域內(nèi)的散點(diǎn)數(shù)據(jù)作為插值約束點(diǎn)，網(wǎng)格化處理僅對(duì)延拓前的區(qū)域進(jìn)行。

圖4 圖像重構(gòu)數(shù)據(jù)周期延拓

3 應(yīng)用效果

數(shù)據(jù)圖像經(jīng)過(guò)著色處理后即可成像。圖5（a）是原始數(shù)據(jù)成像結(jié)果，該圖由8×16個(gè)原始數(shù)據(jù)繪制而成，圖5（b）是根據(jù)上述方法提取到的22個(gè)特征點(diǎn) （圖5（b）中白色點(diǎn)）重構(gòu)而成；對(duì)比圖5（a）和圖5（b）不難發(fā)現(xiàn)，重構(gòu)圖像與原始圖像基本一致，其數(shù)據(jù)信息壓縮比超過(guò)1∶5，大部分變化幅度較小的數(shù)據(jù)被過(guò)濾掉。

4 結(jié)論與認(rèn)識(shí)

通過(guò)在周向和縱向上對(duì)圖像進(jìn)行特征提取，可以壓縮掉大部分?jǐn)?shù)據(jù)變化程度不明顯的數(shù)據(jù)信息。壓縮后圖像數(shù)據(jù)量大幅度降低。壓縮數(shù)據(jù)在周向上進(jìn)行周期延拓，可以確保壓縮數(shù)據(jù)網(wǎng)格化插值處理后在橫向邊界保持連續(xù)性。特征點(diǎn)提取需要反映圖像的變化，如果圖像變化復(fù)雜，則特征點(diǎn)增多，數(shù)據(jù)壓縮效果會(huì)降低。由特征點(diǎn)重構(gòu)的圖像只保留了圖像變化明顯的信息，不明顯的信息被壓縮掉，因此該圖像壓縮方法屬于有損壓縮。

圖5 原始數(shù)據(jù)成像與特征點(diǎn)重構(gòu)圖像對(duì)比

[1]吳鵬程，陳一健，楊琳，等 .成像測(cè)井技術(shù)研究現(xiàn)狀及應(yīng)用 [J].天然氣勘探與開發(fā)，2007，30（2）：36～40.

[2]陳瓊，王偉，葛輝 .成像測(cè)井技術(shù)現(xiàn)狀及進(jìn)展 [J].國(guó)外測(cè)井技術(shù)，2007，22（3）：8～10.

[3]Beattie M S，Abdallah A H.Mud pulse telemetry[P].US6421298B1，2002－07－16.

[4]MacLeod N C.Apparatus and method for down－h(huán)ole EM telemetry while drilling [P].US4739325A，1988－04－19.

[5]Orion II Data Compression Platform [EB/OL].http：//www.slb.com/services/drilling/mwd＿lwd/mwd/orion＿ii＿mwd.aspx，2013－12－11.

[6]王偉 .CGDS172NB近鉆頭地質(zhì)導(dǎo)向鉆井技術(shù)在江漢油田的應(yīng)用 [J].鉆采工藝，2011，34（5）：34～36.

[7]岳愛忠，張健民，何彪，等 .隨鉆方位密度成像測(cè)井?dāng)?shù)值模擬研究 [A].中國(guó)核科學(xué)技術(shù)進(jìn)展報(bào)告——中國(guó)核學(xué)會(huì)2009年學(xué)術(shù)年會(huì)論文集 [C].北京：原子能出版社，2010.304～308.

[8]侯爽 .隨鉆方位密度測(cè)井方法基礎(chǔ)研究 [D].青島：中國(guó)石油大學(xué) （華東），2011.8～14，56～57.

[9]范斐，龐巨豐，徐佳，等 .井周聲波成像測(cè)井儀原理與應(yīng)用 [J].計(jì)量與測(cè)試技術(shù)，2009，36（8）：74～76.

[10]楊世奪，雷霄，蔡軍，等 .隨鉆電阻率成像測(cè)井在北部灣碳酸鹽巖儲(chǔ)層中的綜合應(yīng)用 [J].測(cè)井技術(shù)，2010，34（2）：177～182.

[11]何小胡，張迎朝，張道軍，等 .成像測(cè)井技術(shù)在重力流沉積研究中的應(yīng)用 [J].測(cè)井技術(shù)，2013，37（1）：103～109，113.

[12]宋代文，曹民權(quán)，張峻清，等 .PeriScope隨鉆成像測(cè)井技術(shù)在水平井中的應(yīng)用 [J].國(guó)外油田工程，2010，26（4）：45～47.

[13]趙軍，閆爽，宋帆，等 .隨鉆成像測(cè)井資料在水平井井周構(gòu)造解釋中的應(yīng)用 [J].石油地球物理勘探，2007，42（s1）：76～78.

[14]汪鋼，王燕 .一種超聲成象測(cè)井圖象無(wú)損壓縮方法 [J].石油儀器，2010，24（1）：71～73.

[15]漢澤西，郭楓，呂飛 .測(cè)井?dāng)?shù)據(jù)的無(wú)損壓縮算法研究 [J].現(xiàn)代電子技術(shù)，2004，27（22）：94～96.

[16]楊子賢 .基于壓縮感知的井下圖像壓縮 [D].合肥：中國(guó)科學(xué)技術(shù)大學(xué)，2012.7～45.

[17]尹麗 .圖象去霧和超分辨率圖像重構(gòu)技術(shù)研究 [D].武漢：武漢理工大學(xué)，2012.27～48.

[18]Douglas P J，Zeng K，Yin Z.System and method for image reconstruction [P].US20080340787，2012－10－09.

[19]楊揚(yáng)，劉哲，張萌 .一種基于全變差模型的欠采樣圖像重構(gòu)方法 [J].紅外與毫米波學(xué)報(bào)，2012，31（2）：153～158.

[20]解成俊，張鐵山 .基于壓縮感知理論的圖像重構(gòu)算法研究 [J].計(jì)算機(jī)應(yīng)用與軟件，2012，29（4）：49～52.

[21]劉兆平，楊進(jìn)，武煒 .地球物理數(shù)據(jù)網(wǎng)格化方法的選取 [J].物探化探，2010，34（1）：93～97.

[22]Lianghui Guo，Xiaohong Meng，Lei Shi.Gridding aeromagnetic data using inverse interpolation [J].Geophysical Journal International，2012，189 （3）：1353～1360.

[23]張偉，覃慶炎，簡(jiǎn)興祥 .自然鄰點(diǎn)插值算法及其在二維不規(guī)則數(shù)據(jù)網(wǎng)格化中的應(yīng)用 [J].物探化探計(jì)算技術(shù)，2011，33（3）：291～295.

[24]孟小紅，侯健全，梁宏英，等 .離散光滑插值方法在地球物理位場(chǎng)中的快速實(shí)現(xiàn) [J].物探與化探，2002，26（4）：302～306.

[25]劉國(guó)峰，孟小紅，張穎 .基于Delaunay剖分和高斯小波函數(shù)的地球物理數(shù)據(jù)規(guī)格化 [J].工程地球物理學(xué)報(bào)，2005，2（6）：403～407.