Seismic modeling is pervasive in exploring the subsurface structure. The propagation of elastic waves in homogenous medium has to be modeled to create synthetic seismograms. A numerical solution of partial differential equations describes the propagation phenomenon in elastic medium under the initial and boundary condition that is Clayton Engquist (CE). The subsurface discontinuities like fractures effect the seismic reflections that are used for subsurface imaging. A fractured velocity model with shales and marls sedimentary rocks is built and common depth point (CDP) seismograms with single shot are preprocessed by automatic gain control. The subsurface reflections are modeled by using the first-order autoregressive (AR(1)) model. A comparison of synthetic and real data AR model is made on the basis of average reflectivity, R2 and means square error (MSE). The real data has smaller average reflectivity, -1.80e-10, 93.966 explained variation i.e. R2 and 1.71e-07 minimum MSE. © 2016 Author(s).