Modeling seismic impedance with Markov chains

Abstract

Abstract Acoustic impedance is modeled as a special type of Markov chain, one which is constrained to have a purely exponential correlation function. The stochastic model is parsimoniously described by M parameters, where M is the number of states or rocks composing an impedance well log. The probability mass function of the states provides M-1 parameters, and the 'blockiness' of the log determines the remaining degree of freedom. Synthetic impedance and reflectivity logs constructed using the Markov model mimic the blockiness of the original logs. Both synthetic impedance and reflectivity are shown to be Bussgang, i.e., if the sequence is input into an instantaneous nonlinear device, then the correlation of input and output is proportional to the autocorrelation of the input. The final part of the paper uses the stochastic model in formulating an algorithm that transforms a deconvolved seismogram into acoustic impedance. The resulting function is blocky and free of random walks or sags. Low-frequency information, as provided by moveout velocities, can be easily incorporated into the algorithm.