LASIO库在油气储藏方面的应用

1. 写在前面

我一开始搜集资料的时候发现，有laspy等等的库可以用，我也打算用laspy，但是安装阶段就受了苦了，后面才发现：安装laspy必须先安装其他的几个库。好不容易安装完了，去读取我的las文件居然失败了，报错是说什么signature error。后面才发现，laspy支持的las文件版本有要求。我的是2.0版本，我去查什么可以读2.0版本的文件，就推荐了lasio，于是我就用了lasio。

但是我是一个地道的小白，上面的所谓的“经验”也是碰巧运气，瞎猫碰上死耗子的事儿，如果各位有懂的，我一定认认真真听取建议。

2. lasio库的使用经验

1. 在利用lasio库读取一个las文件时，以下代码可以在读取之后查看las文件的标签，数据尺寸，以及各个标签的数据。

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
# @Time    : 2023/5/7 17:21
# @Author  : TuoTuo
# @FileName: 读取las文件.py
# @Software: PyCharm
# @Blog    ：20205775@cqu.edu.cn

import lasio
import numpy as np
import matplotlib.pyplot as plt

path = "your las file path"#你自己的las文件地址

#读取las文件
las_file = lasio.read(path)
#查看数据尺寸
print(las_file.data.shape)
#查看数据标签
print(las_file.keys())
#获取相应标签下的内容
print(las_file['DEPTH'])
print(las_file['a'])
print(las_file['b'])
#获取索引，一般都是你的测井深度数据
print(las_file.index)

#初步查看数据
#一般测井数据都是深度作为纵坐标
fig = plt.figure()
trackNum=las_file.data.shape[1]
for i in np.arange(1,trackNum):
    ax = fig.add_subplot(1,trackNum-1,i)
    ax.plot(las_file.data[:,i], las_file.index)
    ax.set_xlabel(las_file.keys()[i])
    ax.xaxis.tick_top()
    ax.invert_yaxis()
plt.show()

2. lasio还可以将你的las文件转为pandas的DataFrame结构，但一定注意：转化之后，默认把你的第一列（一般都是测井深度）作为索引，如果你不希望这样，可以运行下面的代码。

df2 = las.df().reset_index()

3.待补充。