一種應(yīng)用mapgis和excel軟件結(jié)合生成LAS數(shù)據(jù)的方法

張松平

【摘 要】為了將工程技術(shù)人員從手繪LAS曲線中解放出來，應(yīng)用了一種mapgis和excel軟件相結(jié)合對(duì)測(cè)井綜合解釋紙質(zhì)剖面數(shù)字化，提取聲波、密度曲線數(shù)據(jù)并生成LAS文件的方法，最終用于合成記錄制作。通過實(shí)踐證明，該方法確實(shí)可行，能有效解決地震解釋中井震結(jié)合時(shí)深轉(zhuǎn)換問題。

【關(guān)鍵詞】Mapgis；數(shù)字化；Las曲線

0 引言

在地震解釋工作中經(jīng)常通過LAS曲線中聲波和密度曲線制作地震合成記錄，由地震合成記錄與目的層地震波組特征的匹配關(guān)系進(jìn)行目的層地質(zhì)特征標(biāo)定，但在沒有LAS曲線或者只能紙質(zhì)曲線的情況下，要在地震與地質(zhì)上進(jìn)行目的層標(biāo)定，而數(shù)字化生成LAS曲線工作起到了橋接作用。以往工程技術(shù)人員繪圖多采用人工手繪成圖，即采用以紙為介質(zhì)的方式繪制各種地質(zhì)圖件，后期再使用Mapgis軟件用矢量化的方法描繪掃描后的紙制圖件，以達(dá)到圖件更統(tǒng)一和美觀的目的。在本文中，介紹了一種應(yīng)用mapgis、excel相結(jié)合的方法[1-5]，對(duì)測(cè)井綜合解釋紙質(zhì)剖面數(shù)字化，并提取聲波與密度曲線。通過最終制作合成記錄的應(yīng)用與時(shí)深轉(zhuǎn)換情況分析，證明此方法是確實(shí)可行，能明顯改善時(shí)深轉(zhuǎn)換應(yīng)用效果。

1 數(shù)字化生成LAS曲線的過程

Mapgis的數(shù)字化過程大體上可以分為如下幾個(gè)步驟：（1）底圖整理；（2）圖形校正與坐標(biāo)變換；（3）矢量化；（4）數(shù)字文件輸出。以上四個(gè)步驟中第（2）（3）步相對(duì)重要一些，下面對(duì)這一方法作進(jìn)一步介紹。

1.1 曲線底圖整理

對(duì)圖紙測(cè)井資料進(jìn)行掃描生成曲線底圖（圖1所示），要求格式為TIFF柵格數(shù)據(jù)格式，并對(duì)掃描的底圖初步校正，圖片分辨率一般為300DPI，原始圖件的變形要很小、清晰度要較高，曲線連續(xù)性好，即能滿足一般的精度要求。

1.2 圖形校正與坐標(biāo)變換

本步驟是將系統(tǒng)中的控制點(diǎn)位內(nèi)插、坐標(biāo)平移旋轉(zhuǎn)，對(duì)圖框坐標(biāo)定位并校正圖片等。如圖1所示，首先將圖1所示的TIFF格式圖片轉(zhuǎn)為mapgis內(nèi)部MSI格式文件，刪去所有自動(dòng)控制點(diǎn)。接著手動(dòng)增加控制點(diǎn)（圖2所示），圖2中控制點(diǎn)一共有6個(gè)，給定第1點(diǎn)坐標(biāo)為（0，1190），按照?qǐng)D中框格給定第6點(diǎn)坐標(biāo)為（200，1210），其他坐標(biāo)依次類推，依照控制點(diǎn)確定其他點(diǎn)位坐標(biāo)。再采用鑲嵌融合處理，對(duì)圖片采用二次多項(xiàng)方式進(jìn)行影像幾可校正，最終保存為MSI格式文件。

1.3 矢量化

矢量化是指把柵格數(shù)據(jù)轉(zhuǎn)換成矢量數(shù)據(jù)的過程，要求以柵格圖像的中心線進(jìn)行矢量化處理，矢量化后回放的圖形與原圖的線重合最大誤差不應(yīng)大于0.15mm。mapgis中設(shè)置系統(tǒng)參數(shù)，采用結(jié)點(diǎn)半徑0.004，矢量化參數(shù)中抽稀因子一般為1，矢量化范圍為全圖范圍。新建曲線DEN.wl并設(shè)置線參數(shù)（顏色、線寬）并交互矢量化，圖3中藍(lán)色曲線AB段為已經(jīng)矢量化后的曲線，黑色曲線BC段為將要矢量化的曲線，為了顯示方便可以將曲線加粗顯示。圖4中為數(shù)字化以后曲線剖面，其中藍(lán)色曲線為DEN曲線，所示密度平均在2.6克每立方厘米左右，煤層密度約為1.4克每立方厘米。同樣方式交互矢量化生成綠色AC曲線。

1.4 數(shù)字提取與LAS文件制作

設(shè)置曲線DEN.wl文件編輯狀態(tài)，并弧段提取曲線，選中圖4中CD 段后，文件轉(zhuǎn)換中輸出mapgis明碼格式文件K10。將文件K10用excel表格打開，如圖5所示，采用等比例計(jì)算得對(duì)應(yīng)深度的密度值，對(duì)數(shù)據(jù)內(nèi)插重采樣，深度域采集間隔為0.05米，再用LAS編輯器輸入DEPTH和DEN兩列數(shù)據(jù)生成LAS文件（圖6所示），同理可以生成AC曲線，其步驟與DEN曲線生成過程一樣。

將生成的LAS曲線文件導(dǎo)入解釋系統(tǒng)，再制作合成記錄，并時(shí)深轉(zhuǎn)換，對(duì)比有K10井曲線（圖7所示）與無K10井曲線的時(shí)深轉(zhuǎn)換結(jié)果（圖8所示）發(fā)現(xiàn)，有K10井?dāng)?shù)據(jù)應(yīng)用時(shí)深轉(zhuǎn)換結(jié)果更加可靠合理，說明以上步驟生成的測(cè)井曲線是對(duì)無井資料數(shù)據(jù)的一個(gè)校正與補(bǔ)充，方法確實(shí)可行，而且具有明顯的應(yīng)用價(jià)值。

2 總結(jié)

通過應(yīng)用mapgis交互矢量化功能，將紙質(zhì)測(cè)井曲線數(shù)字化，應(yīng)用EXCEL進(jìn)行數(shù)據(jù)插值、重新采樣，將重采樣數(shù)據(jù)制作成地震地質(zhì)解釋平臺(tái)通用的LAS曲線數(shù)據(jù)，并成功應(yīng)用于生產(chǎn)中，對(duì)現(xiàn)有地震解釋方案進(jìn)行了校正與調(diào)整，使得解釋成果更加可靠，說明以上方法確實(shí)可行，具有較好的應(yīng)用價(jià)值。應(yīng)用中需要注意如下兩點(diǎn)：（1）底圖整理后曲線連續(xù)性好，要求分辨率300DPI以上。（2）矢量化后回放的圖形與原圖的線重合最大誤差不應(yīng)大于0.15mm。

【參考文獻(xiàn)】

[1]靳皇玉，崔亞茹，許祥彬.應(yīng)用EXCEL和MAPGIS快速生成地質(zhì)鉆孔柱狀圖[J].吉林地質(zhì)，2009，28（2）：102-105.

[2]羅義.應(yīng)用EXCEL和MAPGIS快速生成地浸砂巖型鈾礦鉆孔綜合柱狀圖[J].中國新技術(shù)新產(chǎn)品，2014（18）：25-27.

[責(zé)任編輯：王楠]