解華明 ，胡光道 ，梅紅波 ，秦 鑫，李程俊，3

（1.安徽建筑工業(yè)學(xué)院 環(huán)境與能源工程學(xué)院，合肥 230601；2.中國地質(zhì)大學(xué) 數(shù)學(xué)地質(zhì)與遙感地質(zhì)研究所，武漢 430074；3.中國地質(zhì)大學(xué) 計算機(jī)學(xué)院，武漢 430074）

0 前言

測井技術(shù)作為地下地質(zhì)研究的一種手段，是從地質(zhì)體的物理性質(zhì)來認(rèn)識、了解地質(zhì)體的一種有效方法。通過測井可以詳細(xì)地研究鉆孔巖性，探測有用礦產(chǎn)，詳細(xì)提供計算儲量所必需的參數(shù)（如油層的有效厚度、孔隙度、含油氣飽和度和滲透率等）［1］。測井方法在金屬礦產(chǎn)、石油、煤，以及水文地質(zhì)、工程地質(zhì)的鉆孔工作中都得到了廣泛的應(yīng)用，并已成為地質(zhì)勘探工作中不可或缺的勘探方法之一。應(yīng)用測井方法可以減少鉆井取心工作量，提高勘探速度，降低勘探成本。

對測井?dāng)?shù)據(jù)進(jìn)行必要的預(yù)處理和成圖分析，有利于更好地解譯出鉆孔地層性質(zhì)。原始測井?dāng)?shù)據(jù)為明碼文本或二進(jìn)制數(shù)據(jù)，格式多樣，數(shù)據(jù)量大，包含有較多的無效數(shù)據(jù)，不利于直接分析、解譯，必須有專門的測井?dāng)?shù)據(jù)預(yù)處理和分析成圖軟件輔助進(jìn)行解譯。

目前，已經(jīng)有較多的測井?dāng)?shù)據(jù)預(yù)處理和成圖軟件實現(xiàn)了這些功能。但是，往往功能越強大價格越昂貴，很多功能在生產(chǎn)中不常用，有些軟件繪制生成的圖件只能打印或者保存為柵格圖，后期圖件編輯、整飾不方便。

作者在本文中，討論了基于MapGIS的地理信息系統(tǒng)平臺，進(jìn)行了石油測井?dāng)?shù)據(jù)預(yù)處理與成圖軟件設(shè)計與實現(xiàn)，并實現(xiàn)了方便的測井?dāng)?shù)據(jù)預(yù)處理，快速自動生成多種用戶自定義格式的測井柱狀圖。

1 基本思路

作者在Windows平臺上，基于MapGIS 6.7地理信息系統(tǒng)二次開發(fā)平臺，利用VC 6.0采用面向?qū)ο蠹夹g(shù)進(jìn)行開發(fā)。利用國產(chǎn)成熟的MapGIS平臺開發(fā)成圖系統(tǒng)，開發(fā)工作量小，成本低，生成圖件即為MapGIS文件格式，可以利用MapGIS對生成圖件進(jìn)行后期編輯、整飾和精美打印。

1.1 基于MapGIS開發(fā)測井繪圖系統(tǒng)可行性分析

MapGIS軟件圖形編繪功能強大，是圖表表達(dá)的有力工具，但其是為描述空間實體和對象而設(shè)計。如何利用MapGIS平臺提供的空間對象描述方法來表達(dá)繪圖元素，是開發(fā)測井成圖系統(tǒng)的一個重要問題，其直接影響著方案的可行性。

測井柱狀圖可以看成是由文字符號標(biāo)注、文字段落、復(fù)合表格、曲線、比例線段、花紋圖案所組成的一個復(fù)雜圖形。MapGIS 空間實體分為點狀實體、線狀實體、面狀實體，以及由這些實體聚合而成的聚合實體（如網(wǎng)絡(luò)實體等）?？梢越⑷鐖D1所示的利用MapGIS空間實體描述方法，來表達(dá)測井柱狀圖繪圖元素的對應(yīng)關(guān)系。

1.2 系統(tǒng)功能需求

系統(tǒng)需要實現(xiàn)的主要功能包括以下三個方面。

1.2.1 測井?dāng)?shù)據(jù)預(yù)處理功能

原始測井?dāng)?shù)據(jù)存在存儲格式多樣，不利于統(tǒng)一管理與處理的問題，如有明碼的LAS格式，二進(jìn)制的LA716格式等，在實際工作中需要能夠進(jìn)行相互轉(zhuǎn)換［2、3］。

在實際生產(chǎn)過程中，同一口井往往多次重復(fù)測井，每次測井深度不盡相同，需要根據(jù)實際情況，對測井曲線進(jìn)行分割、合并、裁剪等操作，生成一個完整的符合要求的測井?dāng)?shù)據(jù)文件。

此外，在測井過程中不可避免地存在少量的空值（Null Value）和異常值，需要檢測出這些值［4］，并利用曲線上測點值、下測點值對空值和異常值進(jìn)行插值。

1.2.2 基于數(shù)據(jù)庫的測井曲線管理

測井?dāng)?shù)據(jù)量大，獲取成本高，數(shù)據(jù)價值大，數(shù)據(jù)安全性要求高，在測井?dāng)?shù)據(jù)管理過程中需要進(jìn)行權(quán)限管理、備份。在數(shù)據(jù)分析與處理過程中，需要數(shù)據(jù)共享。因此，需要利用數(shù)據(jù)庫對測井?dāng)?shù)據(jù)進(jìn)行統(tǒng)一管理［5］。

1.2.3 測井柱狀圖自動生成

測井柱狀圖可以圖形化表達(dá)測井曲線及解譯成果，是測井分析與成果表達(dá)的有效手段，繪制測井柱狀圖對測井地質(zhì)工作具有重要意義。生成的測井柱狀圖需要能夠進(jìn)行后期的圖形編輯，這樣可以方便地進(jìn)行圖形整飾。

1.3 系統(tǒng)結(jié)構(gòu)

系統(tǒng)采用三層結(jié)構(gòu)，將數(shù)據(jù)訪問層、業(yè)務(wù)邏輯層、用戶界面層分開實現(xiàn)，如下頁圖2所示。

在數(shù)據(jù)訪問層，由各種測井?dāng)?shù)據(jù)文件存取函數(shù)、測井?dāng)?shù)據(jù)庫存取類、模板文件讀取解析類、巖性代碼存取類等構(gòu)成。實現(xiàn)與數(shù)據(jù)存取相關(guān)的功能方法，包括各種格式測井?dāng)?shù)據(jù)文件的存取，基于關(guān)系型測井?dāng)?shù)據(jù)庫的存取，柱狀圖模板文件讀取與解釋，以文件數(shù)據(jù)庫形式存在的巖性代碼庫存取。

圖1 利用MapGIS空間實體描述方法表達(dá)測井柱狀圖繪圖元素Fig.1 Element description of logging histogram by MapGIS entity

業(yè)務(wù)邏輯層是由測井?dāng)?shù)據(jù)預(yù)處理類、柱狀圖數(shù)據(jù)處理類、繪圖類等構(gòu)成。

（1）測井?dāng)?shù)據(jù)預(yù)處理類實現(xiàn)測井?dāng)?shù)據(jù)的組織、格式轉(zhuǎn)換、分割、合并、上提井下放井轉(zhuǎn)換、空值和異常值檢測、插值、各曲線統(tǒng)計分析、數(shù)據(jù)狀態(tài)檢查等功能。

（2）柱狀圖數(shù)據(jù)處理類實現(xiàn)繪制測井柱狀圖所需數(shù)據(jù)的組織、正確性驗證、巖性代號或巖性名稱到MapGIS區(qū)填度代碼轉(zhuǎn)換等工作。

（3）繪圖類依據(jù)模板類的解析結(jié)果得到柱狀圖布局等參數(shù)，根據(jù)柱狀圖數(shù)據(jù)處理類的數(shù)據(jù)，進(jìn)行測井柱狀圖繪圖，繪圖結(jié)果為MapGIS的點、線、面文件。

用戶界面層是將系統(tǒng)功能按用戶使用邏輯進(jìn)行組織，封裝為小模塊。系統(tǒng)共分為六個小模塊：①測井?dāng)?shù)據(jù)庫管理；②測井?dāng)?shù)據(jù)格式轉(zhuǎn)換；③測井?dāng)?shù)據(jù)分割、合并；④測井?dāng)?shù)據(jù)編輯；⑤測井柱狀圖繪圖；⑥巖性圖案管理。

2 關(guān)鍵技術(shù)

2.1 測井?dāng)?shù)據(jù)格式轉(zhuǎn)換

測井?dāng)?shù)據(jù)格式多樣，實現(xiàn)任意兩種格式之間的直接轉(zhuǎn)換工作量巨大，擴(kuò)展支持的數(shù)據(jù)格式也非常困難［2］。作者采用中間格式［5、6］（轉(zhuǎn)換媒介）進(jìn)行過度，實現(xiàn)多種格式測井?dāng)?shù)據(jù)格式轉(zhuǎn)換的方法，可以很好地解決這個問題（如下頁圖3所示）。

因為各種格式的頭記錄塊不盡相同，有的頭記錄塊信息豐富，有的較少，所以中間格式需盡可能地存儲頭記錄塊信息，以保證測井?dāng)?shù)據(jù)頭記錄塊信息在轉(zhuǎn)換過程中不丟失或少丟失。

2.2 測井曲線異常值檢測與插值

由于儀器等原因會造成測井曲線中存在異常值［7］，這個值不是巖性特征的物理性質(zhì)反映，必需通過統(tǒng)計方法檢測出來，再通過其相鄰的上下測值對其進(jìn)行插值處理，消除異常。

在數(shù)據(jù)分析中，異常值檢測的方法有多種，常用 的 有Dixon 檢 驗 法、Grubbs 檢 驗 法 方 法［4、8］。Dixon檢驗法和Grubbs檢驗法都可以檢驗一組測量值和剔除一組測量值中的異常值。當(dāng)一組測量值中僅有一個可疑值時，效果良好；但是當(dāng)有多個可疑值時，特別當(dāng)一組觀測值最大值和最小值同時為可疑值，或在最大值（或最小值）一側(cè)同時出現(xiàn)兩個可疑值時，這兩種方法都不是很理想。作者在本系統(tǒng)中，通過提供多種異常值檢測方法共用戶混合使用來減小該問題。

圖2 系統(tǒng)總體結(jié)構(gòu)圖Fig.2 The overall structure of the system

圖3 通過中間格式進(jìn)行數(shù)據(jù)格式轉(zhuǎn)換方法示意圖Fig.3 Data format conversion by the intermediate format

測井曲線異常值檢測與普通觀測數(shù)據(jù)異常值檢測不完全相同，因為其測點之間距離非常?。ㄒ话銥?.1m 或0.15m），具有在深度方向上觀測值局部平滑的特點。但由于整條測井曲線的觀測值變化較大，故在進(jìn)行異常值檢測時不宜采用大樣本，可以選擇以待測可疑點為中心的局部區(qū)域作為檢測樣本進(jìn)行異常值檢測。當(dāng)該局部區(qū)域內(nèi)有少量（一個或二個）空值時，則跳過空值點，根據(jù)實驗，選擇7個～11個測點作為樣本進(jìn)行檢測，效果較好。

因測井曲線具有局部平滑的特點，所以測井曲線空值與異常值可利用該值為中心的局部區(qū)域上的下測值進(jìn)行插值。插值方法可以選擇最小二乘線性插值、最小二乘二次多項式、三次多項式進(jìn)行插值。

2.3 測井曲線數(shù)據(jù)庫設(shè)計

不同的測井?dāng)?shù)據(jù)格式存儲的測井信息不完全相同，測井曲線數(shù)據(jù)庫需要盡可能地存儲所有頭記錄與曲線信息。此外，數(shù)據(jù)庫中應(yīng)能按類別來組織測井?dāng)?shù)據(jù)。

圖4（見下頁）為系統(tǒng)中測井曲線數(shù)據(jù)庫概念數(shù)據(jù)模型。在該模型中，“測井曲線”實體中“曲線列序號”字段存儲的是該曲線數(shù)據(jù)在“測井記錄”實體中的字段列序數(shù)?！皽y井記錄”實體有52 個字段，其中從“C1”至“C50”50 個字段用于存儲測井曲線值，一個測井文件一般不會超過50條曲線，50個字段足夠存儲。用這種方法將這兩個實體聯(lián)系起來而不用聯(lián)系實體［9］，是因為這種數(shù)據(jù)庫結(jié)構(gòu)具有較好的數(shù)據(jù)庫插入、查詢效率。

2.4 基于模板的測井柱狀圖自動生成

地質(zhì)柱狀圖在地質(zhì)、石油勘探等行業(yè)中是非常重要的圖件，開發(fā)地質(zhì)柱狀圖自動生成軟件時存在兩個難點：①不同行業(yè)部門的柱狀圖的格式和記錄的數(shù)據(jù)一般都不相同，如何開發(fā)一個用戶可定制格式的柱狀圖生成程序，成為軟件開發(fā)過程中的一個普遍存在的難點；②地質(zhì)柱狀圖中需要設(shè)置的參數(shù)太多，各個位置文字的字體、大小、顏色、布局要求都不相同，復(fù)合表格、曲線的線型也不盡相同，這些都需要提供給用戶設(shè)置，如果這些參數(shù)都通過對話框與用戶交互，則會使人機(jī)界面變得相當(dāng)復(fù)雜。

本系統(tǒng)測井柱狀圖是作者利用前人提出的模板和圖道思想［10～12］并作改進(jìn)，很好地解決了上述兩個問題，實現(xiàn)了用戶可定制的測井柱狀圖自動生成。可以將測井柱狀圖分解為如表1 所示（見下頁）的圖面元素類型。各個類型按其更新條件可以分為兩類：①格式型元素，只有當(dāng)柱狀圖格式修改了才需要更新；②數(shù)據(jù)型元素，只要是不同的測井?dāng)?shù)據(jù)繪制柱狀圖就需要更新。

利用MapGIS建立模板文件，對于格式型元素直接在模塊文件繪出，對于數(shù)據(jù)型元素在模板文件中繪制出元素在圖面布局中的位置，并在圖元屬性中存儲對應(yīng)的數(shù)據(jù)名稱、繪圖參數(shù)。

在自動成圖時，解析模板文件。對于格式型元素照抄到生成的圖件中；對于數(shù)據(jù)型元素，讀取其繪圖位置、繪圖參數(shù)和對應(yīng)數(shù)據(jù)項名稱，從柱狀圖繪圖數(shù)據(jù)中找到對應(yīng)的數(shù)據(jù)，利用繪圖參數(shù)進(jìn)行圖形繪制。

3 系統(tǒng)實現(xiàn)與應(yīng)用

基于上述研究內(nèi)容，作者開發(fā)并完成了測井?dāng)?shù)據(jù)預(yù)處理和測井柱狀圖生成軟件，實現(xiàn)了現(xiàn)場測井?dāng)?shù)據(jù)的預(yù)處理和快速自動生成多種格式的測井柱狀圖。該軟件在我國東海西湖凹陷石油勘探開發(fā)的初步應(yīng)用中效果良好（見后面圖5）。

圖4 測井曲線數(shù)據(jù)庫概念數(shù)據(jù)模型Fig.4 Conceptual data model of logging database

表1 測井柱狀圖元素分類Tab.1 Element classification of logging histogram

4 結(jié)論

（1）基于MapGIS開發(fā)測井?dāng)?shù)據(jù)預(yù)處理和成圖軟件，可以實現(xiàn)快速、方便地繪制多種格式測井柱狀圖，繪制出的圖件可以滿足實際生產(chǎn)需要，效果良好。

（2）采用中間格式進(jìn)行過度，可以在實際應(yīng)用中實現(xiàn)多種格式的測井?dāng)?shù)據(jù)格式相互轉(zhuǎn)換的方法，相比較傳統(tǒng)的兩兩測井格式轉(zhuǎn)換的方法，節(jié)約了編程工作量。

（3）傳統(tǒng)異常值檢測的Dixon檢驗法、Grubbs檢驗法，可以應(yīng)用于測井觀測異常值的檢測。但當(dāng)有多個可疑值時，檢測方法可能不理想，需要多種檢測方法綜合應(yīng)用。

（4）在測井?dāng)?shù)據(jù)庫設(shè)計中，對實體不能劃分過細(xì)，需要綜合考慮測井?dāng)?shù)據(jù)插入、查詢的效率。

（5）利用模板方法進(jìn)行多格式測井柱狀圖自動繪制技術(shù)是可行、實用的，可以滿足實際生產(chǎn)的需要。

圖5 軟件自動生成的兩種格式的測井柱狀圖界面Fig.5 Two types logging histogram automatically generated by the system

［1］李舟波.地球物理測井?dāng)?shù)據(jù)處理與綜合解釋［M］.長春：吉林大學(xué)出版社，2003.

［2］胡佳慶.測井?dāng)?shù)據(jù)常用格式及格式轉(zhuǎn)換［J］.海洋石油，2003，23（1）：75.

［3］胡振平，王昌德，王本奇，等.測井?dāng)?shù)據(jù)格式轉(zhuǎn)換系統(tǒng)［J］.測井技術(shù)，2005，29（4）：368.

［4］萬應(yīng)明，高峻，董建平，等.多測井曲線合成應(yīng)用方法初探［J］.石油物探，2005，44（1）：71.

［5］劉擁軍，楊永發(fā)，朱春紅，等.長慶油田測井信息數(shù)據(jù)庫的系統(tǒng)設(shè)計［J］.測井技術(shù)，2005，29（3）：233.

［6］熊曉東，郭麥成，袁洪清，等.多功能測井?dāng)?shù)據(jù)處理軟件的開發(fā)［J］.江漢石油學(xué)院學(xué)報，2002，24（3）：76.

［7］高效曾.測井?dāng)?shù)據(jù)及其誤差估計［J］.測井技術(shù)，1993，17（5）：317.

［8］GB／T 4883－2008，數(shù)據(jù)的統(tǒng)計處理和解釋－正態(tài)樣本離群值的判斷和處理［S］.

［9］王珊，薩師煊.數(shù)據(jù)庫系統(tǒng)概論第四版［M］.北京：高等教育出版社，2006.

［10］劉剛，田宜平，吳沖龍.可分幅式柱狀圖的參數(shù)化計算機(jī)輔助設(shè)計［J］.計算機(jī) 工程與設(shè)計，2005，26（9）：2277.

［11］遲文學(xué)，陳建強，許哲，等.鉆孔柱狀圖中緩沖線繪制技術(shù)［J］.地質(zhì)科技情報，2006，25（5）：87.

［12］花衛(wèi)華，鄧偉萍.一種通用鉆孔柱狀圖的自動生成技術(shù)［J］.工程圖學(xué)報，2006（4）：14.