A computer-implemented method of geophysical exploration within a predetermined volume of the earth containing subsurface structural and stratigraphic features conducive to the generation, migration, accumulation, or presence of hydrocarbons, wherein is provided two seismic surveys each having a plurality of seismic traces associated therewith, comprising the steps of: a. accessing a first of said two seismic surveys; b. reading a selected number of input traces from said plurality of seismic traces associated with said first one of said two seismic surveys; c. accessing a second of said two seismic surveys; d. reading said selected number of input traces from said plurality of seismic traces associated with said second one of said two seismic surveys; the method being characterized in e. calculating a forward dual-tree complex wavelet transform, DCWT, on said read input traces from said first seismic survey, thereby forming a first DCWT dataset; f. calculating a forward DCWT on said read input traces from said second seismic survey, thereby forming a second DCWT dataset; g. matching said first DCWT dataset to said second DCWT dataset to produce a matched DCWT dataset; h. calculating an inverse DCWT on said matched DCWT dataset to produce a plurality of matched seismic traces; and, i. using said plurality of matched seismic traces in seismic exploration for subsurface hydrocarbons within said predetermined volume of the earth, wherein step (g) comprises the step of: (gl) matching said first DCWT dataset and said second DCWT dataset to produce a matched DCWT dataset using by calculating the quantity Z ′ t x y s o ri = Z 1 t x y s o ri ∗ Z 2 t x y s o ri ∗ Envp t x y s o, where, Envp t x y s o = 1 Z 1 t x y s o real 2 + Z 1 t x y s o imaginary 2, and where Z ' (t,x,y,s,o,ri) is said matched DCWT dataset, Z1 (t,x,y,s,o,ri) is said DCWT first dataset, Z 2 (t,x,y,s,o,ri) is said DCWT second dataset, t is time, x is an X coordinate vector associated with said at least sixteen first seismic survey traces, y is an Y coordinate vector associated with said at least sixteen first seismic survey DCWT traces, o is an orientation vector, s is a scale vector, and, ri is a vector that contains either a real or an imaginary component depending on the context. a. accessing a first of said two seismic surveys; b. reading a selected number of input traces from said plurality of seismic traces associated with said first one of said two seismic surveys; c. accessing a second of said two seismic surveys; d. reading said selected number of input traces from said plurality of seismic traces associated with said second one of said two seismic surveys; the method being characterized in e. calculating a forward dual-tree complex wavelet transform, DCWT, on said read input traces from said first seismic survey, thereby forming a first DCWT dataset; f. calculating a forward DCWT on said read input traces from said second seismic survey, thereby forming a second DCWT dataset; g. matching said first DCWT dataset to said second DCWT dataset to produce a matched DCWT dataset; h. calculating an inverse DCWT on said matched DCWT dataset to produce a plurality of matched seismic traces; and, i. using said plurality of matched seismic traces in seismic exploration for subsurface hydrocarbons within said predetermined volume of the earth, wherein step (g) comprises the step of: (gl) matching said first DCWT dataset and said second DCWT dataset to produce a matched DCWT dataset using by calculating the quantity Z ′ t x y s o ri = Z 1 t x y s o ri ∗ Z 2 t x y s o ri ∗ Envp t x y s o, where, Envp t x y s o = 1 Z 1 t x y s o real 2 + Z 1 t x y s o imaginary 2, and where Z ' (t,x,y,s,o,ri) is said matched DCWT dataset, Z1 (t,x,y,s,o,ri) is said DCWT first dataset, Z 2 (t,x,y,s,o,ri) is said DCWT second dataset, t is time, x is an X coordinate vector associated with said at least sixteen first seismic survey traces, y is an Y coordinate vector associated with said at least sixteen first seismic survey DCWT traces, o is an orientation vector, s is a scale vector, and, ri is a vector that contains either a real or an imaginary component depending on the context. Z ' (t,x,y,s,o,ri) is said matched DCWT dataset, Z1 (t,x,y,s,o,ri) is said DCWT first dataset, Z 2 (t,x,y,s,o,ri) is said DCWT second dataset, t is time, x is an X coordinate vector associated with said at least sixteen first seismic survey traces, y is an Y coordinate vector associated with said at least sixteen first seismic survey DCWT traces, o is an orientation vector, s is a scale vector, and, ri is a vector that contains either a real or an imaginary component depending on the context.