Using a 3-D finite difference method with a rotated-staggered-grid (RSG) scheme we generated synthetic seismograms for a reservoir model consisting of three horizontal layers with the middle layer containing parallel, equally spaced fractures. By separating and analyzing the backscattered signals in the FK domain, we can obtain an estimate of the fracture spacing. The fracture spacing is estimated by taking one-half of the reciprocal of the dominant wavenumber of the backscattered energy in data acquired normal to the fractures. FK analysis for fracture spacing estimation was successfully applied to these model results, with particular focus on PS converted waves. The method was then tested on data from the Emilio Field. The estimated fracture spacing from the dominant wavenumber values in time windows at and below the reservoir level is 25-40m. A second approach for fracture spacing estimation is based on the observation that interference of forward and backscattered energy from fractures introduces notches in the frequency spectra of the scattered wavefield for data acquired normal to the fracture strike. The frequency of these notches is related to the spacing of the fractures. This Spectral Notch Method was also applied to the Emilio data, with the resulting range of fracture spacing estimates being 25-50m throughout the field. The dominant spacing fracture spacing estimate is about 30-40 m, which is very similar to the estimates obtained from the FK method.