Las-py

las-py is a zero-dependency Python library for parsing .Las file (Geophysical/Canadian well log files).

Currently supports only version 2.0 of LAS Specification. For more information about this format, see the Canadian Well Logging Society web page

• What's new in 1.1.0

  • Export to csv
  • Export to csv without rows containing null values
  • Bug fixes

• To Install

      $pip insatll las-py
  

• Usage

      from las_py import Laspy
  

      my_las = Laspy('path_to_las_file.las')
  

• Read data

     data = my_las.data
     print(data)
     #[[2650.0, 177.825, -999.25, -999.25], [2650.5, 182.5, -999.25,-999.25], [2651.0,180.162, -999.25, -999.25], [2651.5, 177.825, -999.25, -999.25], [2652.0, 177.825, -999.25, -999.25] ...]
  

      # get data with rows that has null value stripped
     data = my_las.data_stripped
     print(data)
     #[[2657.5, 212.002, 0.16665, 1951.74597], [2658.0, 201.44, 0.1966, 1788.50696], [2658.5, 204.314, 0.21004, 1723.21204], [2659.0, 212.075, 0.22888, 1638.328], [2659.5, 243.536, 0.22439, 1657.91699]...]
  

• Get the log headers

      headers = my_las.header
      print(headers)
      # ['DEPTH', 'GR', 'NPHI', 'RHOB']
  

• Get the log headers descriptions

      hds_and_desc = my_las.header_and_descr
      print(hds_and_desc)
      # {DEPTH': 'DEPTH', 'GR': 'Gamma Ray', 'NPHI': 'Neutron Porosity','RHOB': 'Bulk density'}
  

• Get a particular column, say Gamma Ray log

      GR = my_las.column('GR')
      print(GR)
      # [-999.25, -999.25, -999.25, -999.25, -999.25, 122.03, 123.14, ...]
  

      # get column with null values stripped
      GR = my_las.column_stripped('GR')
      print(GR)
      # [61.61, 59.99, 54.02, 50.87, 54.68, 64.39, 77.96, ...]
  

  Note this returns the column, after all the data has been stripped off their null values, which means that valid data in a particular column would be stripped off if there is another column that has a null value at that particular row

• Get the Well Parameters

  Presents a way of accessing the details individual well parameters.

  The details include the following:

    1. descr - Description/ Full name of the well parameter
    2. units - Its unit measurements
    3. value - Value
  

    start = my_las.well.STRT.value # 1670.0
    stop = my_las.well.STOP.value #  1669.75
    null_value = my_las.well.NULL.value #  -999.25
    # Any other well parameter present in the file, canbe gotten with the same syntax above
  

• Get the Curve Parameters

  Presents a way of accessing the details individual log columns.

  The details include the following:

    1. descr - Description/ Full name of the log column
    2. units - Unit of the log column measurements
    3. value - API value of the log column
  

    NPHI = my_las.curve.NPHI.descr # 'Neutron Porosity'
    RHOB = my_las.curve.RHOB.descr # 'Bulk density'
    # This is the same for all log column present in the file
  

• Get the Parameters of the well

  The details include the following:

    1. descr - Description/ Full name of the log column
    2. units - Unit of the log column measurements
    3. value - API value of the log column
  

    BHT = my_las.param.BHT.descr # 'BOTTOM HOLE TEMPERATURE'
    BHT_valaue = my_las.param.BHT.value # 35.5
    BHT_units = my_las.param.BHT.units # 'DEGC'
    # This is the same for all well parameters present in the file
  

• Get the number of rows and columns

      rows = my_las.row_counts # 4
      columns = my_las.column_counts # 3081
  

• Get the version and wrap

      version = my_las.version # '2.0'
      wrap = my_las.wrap # 'YES'
  

• Get other information

      other = my_las.other
      print(other)
      # Note: The logging tools became stuck at 625 metres causing the data
      # between 625 metres and 615 metres to be invalid.
  

• Export to CSV

  This writes a csv file to the current working directory, with headers of the well and data section only.

      my_las.to_csv('result')
      # result.csv has been created Successfully!
  

  result.csv

  DEPT	RHOB	GR	NPHI
  0.5	-999.25	-999.25	-0.08
  1.0	-999.25	-999.25	-0.08
  1.5	-999.25	-999.25	-0.04
  ...	...	...	...
  1.3	-999.25	-999.25	-0.08

  Or get the version of csv with null values stripped

      my_las.to_csv_stripped('clean')
      # clean.csv has been created Successfully!
  

  clean.csv

  DEPT	RHOB	GR	NPHI
  80.5	2.771	18.6	-6.08
  81.0	2.761	17.4	-6.0
  81.5	2.752	16.4	-5.96
  ...	...	...	...
  80.5	2.762	16.2	-5.06

• Support

  las-py is an MIT-licensed open source project. You can help it grow by becoming a sponsor/supporter. Donate on Patreon