The seismic velocity at the formation varies widely with physical properties in the layers. These features on seismic shot gathers are not capable of reproducing normally by numerical modeling of homogeneous medium, so that we need that of random inhomogeneous medium instead. In this study, we conducted Gaussian autocorrelation function (ACF), exponential autocorrelation function and von Karman autocorrelation function for getting inhomogeneous velocity model and applied a simple geological model. According to the results, von Karman autocorrelation function showed short wavelength to the inhomogeneous velocity medium. For numerical modeling for a gas hydrate, we determined a geological model based on field data set gathered in the East sea. The numerical modeling results showed that the von Karman autocorrelation function could properly describe scattering phenomena in the gas hydrate velocity model which contains an inhomogeneous layer. Besides, bottom-simulating-reflectors and scattered waves which appear at seismic shot gather of the field data showed properly in the inhomogeneous numerical modeling.