一种应用mapgis和excel软件结合生成LAS数据的方法

张松平

【摘 要】为了将工程技术人员从手绘LAS曲线中解放出来，应用了一种mapgis和excel软件相结合对测井综合解释纸质剖面数字化，提取声波、密度曲线数据并生成LAS文件的方法，最终用于合成记录制作。通过实践证明，该方法确实可行，能有效解决地震解释中井震结合时深转换问题。

【关键词】Mapgis；数字化；Las曲线

0 引言

在地震解释工作中经常通过LAS曲线中声波和密度曲线制作地震合成记录，由地震合成记录与目的层地震波组特征的匹配关系进行目的层地质特征标定，但在没有LAS曲线或者只能纸质曲线的情况下，要在地震与地质上进行目的层标定，而数字化生成LAS曲线工作起到了桥接作用。以往工程技术人员绘图多采用人工手绘成图，即采用以纸为介质的方式绘制各种地质图件，后期再使用Mapgis软件用矢量化的方法描绘扫描后的纸制图件，以达到图件更统一和美观的目的。在本文中，介绍了一种应用mapgis、excel相结合的方法[1-5]，对测井综合解释纸质剖面数字化，并提取声波与密度曲线。通过最终制作合成记录的应用与时深转换情况分析，证明此方法是确实可行，能明显改善时深转换应用效果。

1 数字化生成LAS曲线的过程

Mapgis的数字化过程大体上可以分为如下几个步骤：（1）底图整理；（2）图形校正与坐标变换；（3）矢量化；（4）数字文件输出。以上四个步骤中第（2）（3）步相对重要一些，下面对这一方法作进一步介绍。

1.1 曲线底图整理

对图纸测井资料进行扫描生成曲线底图（图1所示），要求格式为TIFF栅格数据格式，并对扫描的底图初步校正，图片分辨率一般为300DPI，原始图件的变形要很小、清晰度要较高，曲线连续性好，即能满足一般的精度要求。

1.2 图形校正与坐标变换

本步骤是将系统中的控制点位内插、坐标平移旋转，对图框坐标定位并校正图片等。如图1所示，首先将图1所示的TIFF格式图片转为mapgis内部MSI格式文件，删去所有自动控制点。接着手动增加控制点（图2所示），图2中控制点一共有6个，给定第1点坐标为（0，1190），按照图中框格给定第6点坐标为（200，1210），其他坐标依次类推，依照控制点确定其他点位坐标。再采用镶嵌融合处理，对图片采用二次多项方式进行影像几可校正，最终保存为MSI格式文件。

1.3 矢量化

矢量化是指把栅格数据转换成矢量数据的过程，要求以栅格图像的中心线进行矢量化处理，矢量化后回放的图形与原图的线重合最大误差不应大于0.15mm。mapgis中设置系统参数，采用结点半径0.004，矢量化参数中抽稀因子一般为1，矢量化范围为全图范围。新建曲线DEN.wl并设置线参数（颜色、线宽）并交互矢量化，图3中蓝色曲线AB段为已经矢量化后的曲线，黑色曲线BC段为将要矢量化的曲线，为了显示方便可以将曲线加粗显示。图4中为数字化以后曲线剖面，其中蓝色曲线为DEN曲线，所示密度平均在2.6克每立方厘米左右，煤层密度约为1.4克每立方厘米。同样方式交互矢量化生成绿色AC曲线。

1.4 数字提取与LAS文件制作

设置曲线DEN.wl文件编辑状态，并弧段提取曲线，选中图4中CD 段后，文件转换中输出mapgis明码格式文件K10。将文件K10用excel表格打开，如图5所示，采用等比例计算得对应深度的密度值，对数据内插重采样，深度域采集间隔为0.05米，再用LAS编辑器输入DEPTH和DEN两列数据生成LAS文件（图6所示），同理可以生成AC曲线，其步骤与DEN曲线生成过程一样。

将生成的LAS曲线文件导入解释系统，再制作合成记录，并时深转换，对比有K10井曲线（图7所示）与无K10井曲线的时深转换结果（图8所示）发现，有K10井数据应用时深转换结果更加可靠合理，说明以上步骤生成的测井曲线是对无井资料数据的一个校正与补充，方法确实可行，而且具有明显的应用价值。

2 总结

通过应用mapgis交互矢量化功能，将纸质测井曲线数字化，应用EXCEL进行数据插值、重新采样，将重采样数据制作成地震地质解释平台通用的LAS曲线数据，并成功应用于生产中，对现有地震解释方案进行了校正与调整，使得解释成果更加可靠，说明以上方法确实可行，具有较好的应用价值。应用中需要注意如下两点：（1）底图整理后曲线连续性好，要求分辨率300DPI以上。（2）矢量化后回放的图形与原图的线重合最大误差不应大于0.15mm。

【参考文献】

[1]靳皇玉，崔亚茹，许祥彬.应用EXCEL和MAPGIS快速生成地质钻孔柱状图[J].吉林地质，2009，28（2）：102-105.

[2]罗义.应用EXCEL和MAPGIS快速生成地浸砂岩型铀矿钻孔综合柱状图[J].中国新技术新产品，2014（18）：25-27.

[责任编辑：王楠]