ARTICLE · 1270658
测井方面的纯个人思维构造的插件
01
PART
CAD 里出图：一条命令搭起规范图框
AUTOLISP · DRAWING
第一套是 AutoCAD 里的 AutoLISP 程序。用法极简：APPLOAD加载主程序，命令行敲ZHCJ，弹出对话框勾选要画的曲线，确定——图框、表头、刻度、深度标尺一次性画完。
数据只要一种格式：一个 txt 一条曲线
每个文件两列——深度、数值，文件名就是曲线名。gr.txt 是自然伽马，sp.txt 是自然电位，声波、电阻率、井温、井径、密度各一个。想换一批数据，覆盖文件即可，程序不用改。
0.81 15.20
0.91 16.70
1.01 15.94
测井成图里面包括了不同比例尺成图， 用cad进行分取值的深度值，可以加上不同测井解释，绘制测井曲线可以选择本次有的相关测井曲线，进行限定的最小和最大。我根据以往测井的也可以把后面激电测井和三分量测井进行说明。
这是得到的数据进行origin 3d绘图的中间可以调整。
读拐点，做分层取值
分层的边界深度写在 fc.txt 里，一行一个拐点深度。点对话框里的【分层取值】，程序按「顶板深度 → 底板深度」取区间内全部数据点，逐层算出每条曲线的最大值、最小值、平均值，直接写成一张对齐的表格。
02
PART
Python 流水线：五个页签走完全程
PYQT5 · WELLY · MATPLOTLIB
CAD 出图强，也没法批量做质量检查。所以又写了第二套：PyQt5 做界面、matplotlib 出图、welly 读 LAS，五个页签对应五道工序，从原始曲线一路走到成果图和 DXF。
— Tab1 数据预处理：加载 6 条 DAT，右侧按道预览，窗口越大图越长
加载 + 网格对齐
LAS 和 DAT 都能读。几份数据采样间隔往往不一致，先统一插值到 0.1 米，后面的交会图、分层统计才有共同的分母。
去尖峰 + 平滑
去尖峰用的是 k 倍 MAD 判据，窗口可选 3 到 11 点；平滑是滑动平均。两步都只改被选中的那一条曲线，改完立刻重画预览。
井斜轨迹
四种算法（最小曲率、曲率半径、正切、平衡正切）任选，按设定间隔重采样，算出垂深和水平位移，可出轨迹投影图，也能导出轨迹 DAT。
— Tab2 全井 QC 与井径分段，下半部分是井径相对钻头的比值曲线
Tab2 最实用的是井径分段。一口井常常有几段变径，若全井只拿一个钻头直径去比，浅部整段都会被误判成扩径。
现在的做法是先按井径众数自动找出变径点，吸附到标准钻头尺寸，再每一段各用各的钻头算比值，判扩径、缩径，顺带给出这段的孔隙度和饱和度是否可信。因为扩径的井段测出的值不是该段岩石的真实值。
后期交会图可以根据不同的参数进行二维交汇和三维交汇。不同圈代表不同的岩石的物理性质，如果标志层和其它岩层分的很开就可以依据这个进行很好的判别标志层。
本井实测：自动识别出 104.8 / 90 / 76 mm 三段钻头，扩径 2.8%、缩径 0.3%、正常 96.9%，井眼质量评为「优」，同时生成可信度曲线 QI（均值 0.99）。
— Tab3 分层-曲线综合图：图头、岩性柱、9 道曲线、深度刻度都在一张图上
Tab3 负责分层。导入分层 DAT（顶深、底深、岩性），一键统计每层各曲线的最大、最小、平均，再按纵向比例出图。图头信息填一次就记住——项目台头、井号、井别、坐标、海拔、日期，PNG 和 DXF 都带着。
Tab4 是一条计算链：Vsh → 孔隙度 → 渗透率 → 饱和度。泥质含量用 GR 基线法，孔隙度优先走密度法，缺密度才用声波估；算之前会做骨架自检，把骨架波速和全井岩石波速对一遍，矛盾值当场标出来。
本井结果：Vsh 均值 0.316，总孔隙度 0.463，有效孔隙度 0.427，含水饱和度 0.332（Archie，RW=2.0）。
— Tab5 综合绘图：左边排道，右边所见即所得，纵向 1:500
Tab5 是给「按自己的习惯排图」准备的：一列一条曲线，可增删、上下移动、改中文道名、单独配色、逐道改宽度（默认 50 mm），不选就不画。比例可选 1:200 到 1:2000，出图按真实毫米计算——本井 1479.8 米，1:500 出 6 道，纸高约 2958 mm。
保存图片、转 DXF 一键完成。DXF 里每条曲线一个独立图层，颜色就是界面上选的那个颜色，拿到 CAD 里能继续编辑、拼图、套图框。
