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.