Skip to content
/ las-py Public

Python library for parsing standard well log files (Geophysical well logs)

License

Notifications You must be signed in to change notification settings

laslibs/las-py

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

10 Commits
 
 
 
 
 
 
 
 
 
 

Repository files navigation

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 product page

  • What's new in 1.1.0

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

        $pip install 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_count # 4
        columns = my_las.column_count # 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

Releases

No releases published

Packages

No packages published

Languages