/DepthConversion

Seismic Depth Conversion using Quadratic Functions

Primary LanguagePython

DepthConversion

Seismic Depth Conversion using Quadratic Functions

  • Given a oneway time depth pairs, a quadratic function is fitted
    resulting in 3 coefficients, a0, a1, and a2

  • Given seismic horizon as an x, y, z file, ztqsegy.py
    will convert the horizon from time to depth or depth to time.

  • Quadratic coefficients can be just one set or supplied as three
    columns in one file or in three seperate files.

  • Probabilities in the depth is computed from the computed
    P10, P50, and P90 of the corresponding coefficients.

  • ztqsegy.py can also read in a segy in time or depth and convert
    that outputting another segy in the other domain.

To integrate this program with the output of Machine learning, an option to supply a sonic dt segy is given. This sonic dt segy would have been generated with trace2log2mb program.

If a time segy and a sonic segy is given then depth conversion will use both to generate a depth segy.

Because the sonic segy is usually noisy, you might want to smooth the resulting depth segy by using --smoothradius option. The radius would spatially smooth all the samples around the center and use the median.