This is some data and ideas for a hackathon project called t20-litho_boundary_from_gamma taking place at the Transform 2020 conference, a fully online geoscience conference. This repo is to provide a dataset for this project.
Estimate lithology boundaries from downhole geophysics GR log, using a continuous wavelet transform algorithm (CWT). Idea was to use an existing packages, found PyWavelets. Based on a paper by Evelyn Jun Hill (here) the CWT family best suited for this kind of work is Gaussian and its derivative, which needs to be used in the algorithm.
A dataset for one hole and one curve (or three holes with four curves) can be found in the folder. Following is a basic strategy to play with.
- Try and get some code to make lithology boundaries based (CWT) algorithm that read the gamma data in and outputs the interval boundaries. Trying to use the code snippet on the home page of the PyWavelets package with 'gaus2'
- Next step would be to view the gamma log, with the associated litho data and the outputted tesselations. Essentially 3 strips.
- compare the the tesselation boundaries with the actual lithology boundaries.
- Repeat until the a reasonable match is found. I think the trick here is dialling in the right.
- When we get some results, maybe even plot using DASH or something flash like that.
-
conditioning the data
-
define the boundary using CWT
-
boundary strength calculation
-
include a lithology log in the plot
-
genererate a mosaic plot (kind of solved?)