A computer-implemented method for calibrating depth in a well to observed seismic data in a subsurface volume of interest, the method being implemented in a computer system that includes one or more physical computer processors and non-transient electronic storage, comprising: a) obtaining, from the non-transient electronic storage, well data corresponding to the well in the subsurface volume of interest, wherein the well data comprises logged well depths specifying depth values of the well in the subsurface volume of interest; b) converting, with the one or more physical computer processors, the well data to converted well data, wherein the converted well data specifies time values corresponding to the well data; c) obtaining, from the non-transient electronic storage, a seismic trace, corresponding to the well in the subsurface volume of interest, and a corresponding wavelet, wherein the seismic trace specifies time values as a function of energy, wherein obtaining the corresponding wavelet includes generating the corresponding wavelet using the seismic trace and reflectivity based on contrast in velocity and density across a subsurface interface, the reflectivity indicating how much subsurface energy is reflected off the subsurface interface; d) generating, with the one or more physical computer processors, a transform model based on one or more probabilities and a modification parameter having scalar values corresponding to one or more of an insertion of a slope in the converted well data, a deletion of the slope in the converted well data, a movement of the slope in the converted well data, and a time shift of the converted well data, wherein the transform model includes a slope parameter and a shift parameter; e) applying, with the one or more physical computer processors, the transform model to the converted well data to transform the converted well data into transformed well data, wherein the transformed well data specifies time values corresponding to the converted well data; f) applying, with the one or more physical computer processors, the corresponding wavelet to the transformed well data to generate a modified signal; g) generating, with the one or more physical computer processors, a correspondence between the modified signal and the seismic trace, wherein the correspondence is a difference between the modified signal and the seismic trace; and h) repeating, with the one or more physical computer processors, steps d) - g) to select a given modified signal and the corresponding transformed well data to tie the well to the seismic trace, based on the correspondence between individual ones of the modified signals and the seismic trace, wherein one of the one or more probabilities used to generate a second transform model, m', is based on a probability, α: α m → m ′ = min 1, β L m ′ L m where m represents a first transform model, β represents the modification parameter, and represents a likelihood for a given transform model. a) obtaining, from the non-transient electronic storage, well data corresponding to the well in the subsurface volume of interest, wherein the well data comprises logged well depths specifying depth values of the well in the subsurface volume of interest; b) converting, with the one or more physical computer processors, the well data to converted well data, wherein the converted well data specifies time values corresponding to the well data; c) obtaining, from the non-transient electronic storage, a seismic trace, corresponding to the well in the subsurface volume of interest, and a corresponding wavelet, wherein the seismic trace specifies time values as a function of energy, wherein obtaining the corresponding wavelet includes generating the corresponding wavelet using the seismic trace and reflectivity based on contrast in velocity and density across a subsurface interface, the reflectivity indicating how much subsurface energy is reflected off the subsurface interface; d) generating, with the one or more physical computer processors, a transform model based on one or more probabilities and a modification parameter having scalar values corresponding to one or more of an insertion of a slope in the converted well data, a deletion of the slope in the converted well data, a movement of the slope in the converted well data, and a time shift of the converted well data, wherein the transform model includes a slope parameter and a shift parameter; e) applying, with the one or more physical computer processors, the transform model to the converted well data to transform the converted well data into transformed well data, wherein the transformed well data specifies time values corresponding to the converted well data; f) applying, with the one or more physical computer processors, the corresponding wavelet to the transformed well data to generate a modified signal; g) generating, with the one or more physical computer processors, a correspondence between the modified signal and the seismic trace, wherein the correspondence is a difference between the modified signal and the seismic trace; and h) repeating, with the one or more physical computer processors, steps d) - g) to select a given modified signal and the corresponding transformed well data to tie the well to the seismic trace, based on the correspondence between individual ones of the modified signals and the seismic trace,