張松平

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

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

0 引言

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

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

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

1.1 曲線底圖整理

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

1.2 圖形校正與坐標變換

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

1.3 矢量化

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

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

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

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

2 總結

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

【參考文獻】

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

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

[責任編輯：王楠]