The Brunswick No. 6 area, which is located in the Bathurst Mining Camp, New Brunswick, Canada, is the focus of this thesis. Almost a decade ago, in order to improve the understanding of the crustal structures and explore for new mineral deposits at depth, three 2D seismic profiles totaling about 30 km and 3D seismic data covering an area of about 38 km2 were acquired from the study area. Petrophysical properties including compressional-wave velocity and density were also measured in two deep boreholes in the area. These data were recovered and reanalyzed, and the improved seismic images interpreted as the main part of this PhD thesis. A prestack DMO and poststack migration algorithm was considered for processing both 2D and 3D data. Processing of 2D data revealed shallow and deep reflections, which correlate well with surface geology. Steeply-dipping reflections, some of which could host mineral deposits, were imaged down to a depth of 6-7 km. Processing of 3D data showed similar results to the processed 2D profiles. Nevertheless, the non-orthogonal nature of the 3D survey, combined with irregular distribution of offsets, azimuths and trace midpoints, caused a severe acquisition footprint masking reflections in the DMO-corrected unmigrated stacked cube. An FK-dip filter in the wavenumber domain was designed to reduce the effects of the acquisition footprint. To better understand wave propagation and scattering effects, calculated acoustic impedance log from the available borehole data was used to estimate vertical scale length using a von Karman autocorrelation function. 2D synthetic models representative of heterogeneity in the area were generated accounting for the estimated scale length. Numerical modeling was used to study the scattering effects on the synthetic models, where some predefined targets were superimposed in the provided 2D heterogeneous medium. The effects of variable source frequency, longer horizontal scale length and petrophysical fluctuations of heterogeneous medium were also investigated. The modeling results indicate that, in the presence of large horizontal, but small vertical scale lengths (structural anisotropy), the identification of mineral deposits is possible in the unmigrated stacked sections, but can be challenging in the migrated sections.