when conducting exploration in accordance with seismic methods , reservoir zones are often characterized by strong perturbations in elastic parameters , in particular compressional wave velocities . lateral variations of compressional and shear wave velocities may be determined in the forward scattering ( transmission ) direction using 3 - c vsp data acquisition geometry . such an analysis of the transmitted seismic wavefield helps to assess how well log data relate to the surrounding reservoir zone . as will be understood to persons of ordinary skill in the art , the effects of heterogeneities on seismic wave propagation can be described in terms of different propagation regimes ( wu , 1989 ): quasi - homogeneous for heterogeneities too small to be seen by seismic waves , rayleigh scattering , mie scattering and small - angle scattering . these scattering regimes cause characteristic amplitude , phase and travel time fluctuation , which can be used to obtain estimates of scale length and thereby assess lateral continuity of lithological units and structure in the immediate vicinity of boreholes . horizontal resolution of exploration seismic data is often discussed in terms of the fresnel zone . for surface and vsp data , the fresnel zone radius increases with increasing depth of investigation . the lateral resolution scale obtained from resonance scattering in accordance with embodiments of the present invention , however , is limited only by the effective frequency content of the seismic signal . in accordance with the operation of embodiments of the invention in an exemplary embodiment depicted in fig2 , seismic sources at the surface generate seismic energy ( compressional and shear wave energy ). part of the seismic energy propagates through the earth to a 3 - c receiver in a borehole ( i . e . along a direct path ) and part of the energy encounters a small - scale heterogeneity in the proximity of the borehole . as such , scattered energy also propagates to and will be recorded at the 3 - c receiver . thus , fig2 illustrates schematically and on location at a volume of earth to be seismically explored an exemplary system 100 adapted for performance in accordance with a method aspect of embodiments of the present invention . system 100 comprises seismic sources 102 , sensors 104 , data acquisition equipment 106 and a data processing system 108 . seismic sources 102 may comprise any broadband controlled ( repeatable ) sources such as vibroseis , explosives , airguns and water guns , sparker , weight - drop , etc . seismic sources 102 may be located at or beneath the surface of the earth in the region of the volume of the earth to be explored . in some embodiments , sensors 104 comprise 3 - c geophones or accelerometers . it is noted that omni - directional sensors such as hydrophones ( detectors sensitive to pressure changes ) are not generally suitable to produce 3 - c signal data . sensors may be positioned in the borehole in accordance with known techniques for gathering 3 - c vsp data . single 3 - c receivers or multiple 3 - c receiver arrays may be used . the resonance scattering technique in accordance with embodiments of the present invention is suitable for zero - offset vsps , walk - away vsps , azimuthal vsps and walk above vsps in vertical or deviated boreholes . as such the seismic sources 102 and sensors 104 may be positioned accordingly . data acquisition equipment 106 may comprise amplifiers , filters , analog to digital signal conversion and seismic data recording equipment for measurement of elastic waves ( compressional and shear waves ) excited by broadband seismic sources 102 at or beneath the surface with single - or multilevel 3 - c sensors 104 as described above . 3 - c sensors [ 104 ] and data acquisition equipment [ 106 ] may be separated by wireline or may be collocated . the resonance scattering technique typically uses full waveform seismic data ( including direct ( down going ), reflected ( up going ) and scattered arrivals ). thus true amplitude 3 - c recordings are utilized . data acquisition equipment 106 records the signal data it acquires to a computer readable medium such as disk , tape or other device ( not shown ). data processing system 108 may comprise a microcomputer programmed in accordance with method aspects of embodiments of the invention to perform resonance scattering processing operations . data processing system 108 comprises , in this exemplary embodiment , at least one processor 109 adapted for 3 - c data signal processing in accordance with a method described further herein with reference to fig3 , such as by software programming or other configuration techniques ; one or more resonance scattering visualization devices 110 ( such as a display screen and / or printer ); and a data integration sub - system 111 including a relational data base with complimentary information . data processing system 108 typically includes one or more input devices such as a pointing device , keyboard or other device for interacting with the system 108 ( all not shown ). data system 108 may be coupled to data acquisition equipment 106 as per fig2 for directly receiving acquired signal data from this equipment 106 or may receive it through other means such as via recorded media or a network ( all not shown ). many alternate embodiments of data processing system 108 may be used . while the data processing system 108 is shown including data integration sub - system 111 , persons skilled in the art will appreciate that such a sub - system 111 may be remotely located and coupled via a network or other connection . further , various display or other visualization devices may be used or coupled remotely to system 108 . system 108 ( i . e . processor 109 ) may produce resonance scattering analysis data for visualisation on a visualisation device coupled to a different data processing system ( not shown ). signal processing ( for example , by data acquisition equipment 106 ) in some embodiments comprises adaptive noise removal and instrument response de - convolution . the resonance scattering technique of the embodiments of the present invention works well when applied to raw ( unconditioned ) 3 - c seismic data ; however , some seismic data processing techniques that are commonly applied to raw data before it is examined may adversely affect this resonance scattering method . therefore it is desirable that wiener filters ( spike and predictive convolution techniques ) and conventional vsp wavefield separation techniques are not applied to the signals for resonance scattering analysis . referring now to fig3 there is shown a flow chart of operations 150 for an exemplary embodiment of a method aspect of the present invention . such operations may be performed by a computer programmed in accordance with programming techniques known to those of ordinary skill in the art to implement the embodiments of the present invention as described herein above . for each depth interval of interest and availability in a set of true amplitude 3 - c vsp data , at step 152 , the 3 - c vsp data is rotated into the ray coordinate system . one component of the three , the ray component , is rotated to point into the direction of the incident pressure wave . windowing of the first arrival signal ( compressional wave ) is performed on the radial component ( that is , the component that points in the direction of the ray ) ( step 154 ). at step 156 , further windowing is performed for the signals on the other two ( i . e . transverse ) directional components while using the same time window applied to the first ( radial ) component . the signals thus share a time component but have respective transverse directional components . at step 158 , a tapering of signals may be performed in order to prevent discontinuities at the edges of the window . tapering may be performed in accordance with gaussian taper techniques . amplitude spectra for each of the 3 - c vsp data is computed ( step 160 ) such as by the use of a fast fourier transform . at step 162 , resonance spectra for the two transverse directions are obtained by dividing respectively the two transverse amplitude spectra components by the radial amplitude spectra component . the operations 150 are repeated in accordance with the various time / depth parameters of the signal data captured to create a multidimensional map of resonance spectra for the volume of earth that is the subject of the seismic exploration . resonance spectra obtained at step 162 characterises lateral heterogeneity in the vicinity of the sensors . resonance spectra obtained from the two transverse components detect small scale heterogeneities in the proximity of the borehole . signal processing as described in steps 152 to 162 for first arrival compressional waves may also be applied to first arrival shear waves . resonance spectra from the map of such generated data may then be integrated with well log , petrophysical and other seismic data using , for example , data integration sub - system 111 to build an enhanced three dimensional geological computerized model around the borehole in accordance with data integration techniques known to those skilled in the art . data generated in accordance with the resonance scattering method described herein may be used to enhance the accuracy and / or scope of these models and the charts , graphs , tables , and other visualisations that may be generated therefrom . in accordance with known techniques , the resonance scattering data may be displayed for further analysis and interpretation . the data may be displayed in various horizontal and vertical time and depth slices . the size of a heterogeneity detected may be determined by dividing an estimated velocity of the seismic wave used to generate the seismic data by the peak frequency of the amplitude spectra . velocity may be estimated using known techniques and knowledge of the medium being explored . further , proximity of a heterogeneity may be determined using a wavelength of the peak frequency of the amplitude spectra . the heterogeneity will typically be within a distance less than 2 to 3 times the length of wavelength . this is a reasonable assumption due to the rate at which a scattered wave decays with distance . direction is easily determined using the known location of the sensor and the seismic source . as noted previously , in contrast to prior art techniques , the lateral resolution scale obtained for resonance scattering is limited only by the effective frequency content of the seismic signal and not in terms of the fresnel zone radius . fig4 shows the lateral resolution for resonance scattering ( obtained from forward scattered wavefield data ) and radius of fresnel zone for surface sources having model parameters : v = 2400 m / s , fmin = 20 hz and fmax = 120 hz . for surface seismic sources , the radius of the fresnel zone increases with increasing depth of investigation . in contrast , the radius of investigation for resonance scattering analysis is controlled by the effective bandwidth of the seismic signal and remains constant at depth . thereby resonance scattering analysis closes the resolution gap between conventional well logs ( decimetre resolution ) and surface seismic methods ( hundred ( s ) of meters ). fig5 a , 5 b and 5 c show resonance spectra for three synthetic reservoir models . the spectra are computed for seismic frequencies from 20 to 120 hz . the vertical scale length for all models is 50 m and the models can be tied to an existing borehole log ; the horizontal scale of the reservoir model has a correlation length of 2000 m ( fig5 a ), 500 m ( fig5 b ) and 50 m ( fig5 c ). as expected , seismic waveform data for models with large horizontal scale length show no evidence of resonance scattering ( fig5 a , 5 b ). for short wavelength horizontal scales lengths , prominent resonance peaks are observed in the resonance spectra ( fig5 c ). resonance scattering analysis of 3 - c vsp data detects evidence for local anomalous physical properties in proximity of the borehole . although well suited to hydrocarbon and mineral exploration , this borehole seismic technique is readily adaptable to geotechnical and environmental problems such as detection of local cavities , faults or fractures , heterogeneities in foundations and tunnels . this is a fundamental consequence of the underlying physics of this method , which ultimately rests on local contrasts in physical material properties ( such as densities , compressional modulus or shear modulus ). although the above description relates to specific embodiments as presently contemplated by the inventors , it is understood that the inventive subject matter in its broad aspect includes mechanical and functional equivalents of the elements described herein . the embodiment ( s ) of the invention described above is ( are ) intended to be exemplary only . the scope of the invention is therefore intended to be limited solely by the scope of the appended claims .