Seismic acquisition footprint removal tool: from prototype to production

[mycarta]

Slack channel

t21-hack-footprint

People

• @mycarta
• @scibbatical
• @leocd91
• @WesleyTheGeolien
• @markogauk

Official repo

https://github.com/mycarta/t21-hack-footprint

Description and background

• PowerPoint presentation here

• For a definition of acquisition footprint read the SEG Wiki

• For some examples, read this blog post, and also go to the Resources and reading material below

• Project background: Elwyn and I have done some work putting together a tool to remove acquisition footprint from seismic data; this has been my longest-lived side project.

• If you are curious about it, please read Chapter 39 of the upcoming 52 Things You Should Know About Geocomputing and then head over to the Tutorial notebook I put together to give it a spin

Goals

The long term goal (AKA, the "dream") would be to create an open-source tool, ideally part of Awesome Open Geoscience, that can be used by geophysicists and geologists to remove footprint from seismic data, WHEN 5D interpolation or other costly post-stack processing works are not an options.
The objectives of the hack are:

• The primary objective of this projects would be to find and eliminate performance bottlenecks in the existing code. We already improved computations by switching from Astropy convolution to Scipy (see here and here)
• Secondary objective, test with noisier examples (we will need to find open data with footprint, F3 has some)
• Tertiary objective would be to create documentation for the tool, either by expanding/completing the Tutorial notebook or a separate document
• Final objective, time permitting, would be to put together a VERY minimal Panel app (for example load data from numpy file >> display a time slice and its 2D spectrum >> derive filter >> save filter to numpy file) and deploy
• I am including this for information purposes only: down the road the tool will need a way to automatically recognize and segment time slices of irregular shape into components of polygonal shape (by whatever means)

Resources and reading material

• A mathematical review of the Fourier transform - SEG Wiki
• An Intuitive Explanation of Fourier Theory
• The 1-D Fourier transform - SEG Wiki
• The 2-D Fourier transform - SEG Wiki
• Elegant SciPy - Chapter 4: Frequency and the Fast Fourier Transform

Skills needed

Elwyn and I will be working on this no-matter-what. But we welcome participants that are interested in this project. Useful skills would be:

• Programming / performance optimization /profiling (Numpy broadcasting / Dask, etc.)
• Geophysics / signal analysis
• Testing
• Documentation
• Ideas - any ideas or suggestion is welcome

[filippo82]

T21 Hackathon Idea Mashup today (or tomorrow depending on your time zone) Wednesday April 14th at 16:00 UTC. Come to Slack to find the link for the meeting!