Method for the interpretation of statistically-related seismic records to yield valuable characteristics, such as gas-bearing potential and lithology of strata

Claim 1

A method for determining fluid hydrocarbon-bearing potential and/or lithology of strata in the earth using high-intensity amplitude events in seismic records generated by conventional centerpoint seismic collection and processing techniques in which the presence of refraction energy within said records has been minimized, comprising the steps of: (a) generating seismic data by said conventional centerpoint seismic collection and processing techniques including recording signals from acoustic discontinuities associated with said strata of interest by positioning and employing an array of sources and detectors such that centerpoints between selected pairs of sources and detectors form a series of centerpoints along a line of survey, said recorded signals being the output of said detectors; (b) by means of automated processing means, statically and dynamically correcting said recorded signals to form corrected traces whereby each of said corrected traces is associated with a centerpoint horizontally midway between a source-detector pair from which said each corrected trace was originally derived; (c) by means of automated processing means, indexing said corrected traces so that each of said corrected traces is identified in its relationship to neighboring traces on the basis of progressive changes in common centerpoint location; (d) by a high-amplitude event of interest, determining from among a series of analytic functions of known mathematical character, a best fit to amplitude vs. horizontal offset variations of a gather of said corrected traces, said gather of traces being ientified with a common centerpoint location and a set of progressively changing horizontal offset values; (e) predicting near and far amplitude vs. time trace projections or representations of said projections for said gather of corrected traces at new offset locations based on said best fitting analytic function, said predicted near and far offset trace projections or representations being identified with offset locations falling on opposite sides of said set of changing horizontal offset values; (f) generating a first envelope of said trace amplitude projections of step (e) associated with near offset locations and a second amplitude envelope of said trace projections of step (e) associated with far offset locations and subtracting the two envelopes one from the other to form a difference envelope; (g) weighting said difference envelope of step (f) by the normalized semblance coefficients wherein after ratios of output-to-input energy are generated by stacking, events in said traces associated with a limited number of phases are better indicated; (h) displaying said weighted difference envelope of step (g) so as to depict amplitude vs. time change as a function of centerpoint coordinate values so that progressive change in said high-intensity amplitude event in said displayed traces is identified as a function of progressive change in centerpoint values.

Claim 2

A method for determining fluid hydrocarbon-bearing potential and/or lithology of strata in the earth using high-intensity amplitude events in seismic records generated by conventional centerpoint seismic collection and processing techniques in which the presence of refraction energy within said records has been minimized, comprising the steps of: (a) generating seismic data by said conventional centerpoint seismic collection and processing techniques including recording signals from acoustic discontinuities associated with said strata of interest by positioning and employing an array of sources and detectors such that centerpoints between selected pairs of sources and detectors form a series of centerpoints along a line of survey, said recorded signals being the output of said detectors; (b) by means of automated processing means, statically and dynamically correcting said recorded signals to form corrected traces whereby each of said corrected traces is associated with a centerpoint horizontally midway between a source-detector pair from which said each corrected trace was originally derived; (c) by means of automated processing means, indexing said corrected traces so that each of said corrected traces is identified in its relationship to neighboring traces on the basis of progressive changes in common centerpoint location; (d) by a high-amplitude event of interest, determining from among a series of analytic functions of known mathematical character, a best fit to amplitude vs. horizontal offset variations of a gather of said corrected traces, said gather of traces being ientified with a common centerpoint location and a set of progressively changing horizontal offset values; (e) predicting near and far amplitude vs. time trace projections or representations of said projections for said gather of corrected traces at new offset locations based on said best fitting analytic function, said predicted near and far offset trace projections or representations being identified with offset locations falling on opposite sides of said set of changing horizontal offset values; (f) generating a first envelope of said trace amplitude projections of step (e) associated with near offset locations and a second amplitude envelope of said trace projections of step (e) associated with far offset locations and subtracting the two envelopes one from the other to form a difference envelope; (g) weighting said difference envelope of step (f) by the normalized semblance coefficients wherein after ratios of output-to-input energy are generated by stacking, events in said traces associated with a limited number of phases are better indicated; (h) displaying said weighted difference envelope of step (g) so as to depict amplitude vs. time change as a function of centerpoint coordinate values so that progressive change in said high-intensity amplitude event in said displayed traces is identified as a function of progressive change in centerpoint values.

Claim 3

The method of clam 1 in which said one of said series of analytic function of known mathematical character of step (d) is selected from a group of linear and quadratic equations of the form: A(x)=C.sub.0 +C.sub.1 x; A(x)=C.sub.0 +C.sub.2 x.sup.2 ; and A(x)=C.sub.0 +C.sub.1 x+C.sub.2 x.sup.2

Claim 4

The method of clam 1 in which said one of said series of analytic function of known mathematical character of step (d) is selected from a group of linear and quadratic equations of the form: A(x)=C.sub.0 +C.sub.1 x; A(x)=C.sub.0 +C.sub.2 x.sup.2 ; and A(x)=C.sub.0 +C.sub.1 x+C.sub.2 x.sup.2

Claim 12

A method for converting an original multitrace seismic record into an improved section having increased capability as to the fluid hydrocarbon-bearing potential and/or lithologic nature of high-intensity amplitude events related to reflections from subsurface strata, said improved section being composed of a plurality of amplitude-versus-centerpoint-and-time traces, said original record consisting of a plurality of multitrace seismic traces of amplitude-versus-horizontal coordinate-and-time generated by conventional centerpoint seismic collection and processing techniques in which the presence of refraction energy within said records has been minimized, each of said traces constituting energy derived in association with a particular source-detector pair of known horizontal offset and of known centerpoint location, and representing, in part, event reflections from said subsurface strata, said conversion comprising the steps of: (a) classifying said original traces on the basis of common but progressively changing horizontal offset values and common centerpoint locations, whereby each resulting trace is identified by a centerpoint location common to at least another trace and a known horizontal offset value; (b) for a high-intensity event, determining with regard to a series of analytic functions of known mathematical character, that one of said series is the best fit to amplitude vs. time variations of said each resulting trace and said at least another trace; (c) predicting near and far amplitude vs. time trace projections for said resulting and another traces at new offset locations based on said one, best fitting analytic function; (d) generating a first series of trace amplitude projections vs. time of step (c) associated with near offset locations, and a second series of trace amplitude projections vs. time associated with far offset locations and subtracting increments of said first and second series normalized to the same time sample, one from the other to form at least a segment of said improved section; (e) weighting said subtracted increments of step (d) by semblance coefficients wherein after ratios of output-to-input energy is generated by stacking, events in said increments associated with a limited number of phases are better indicated; (f) displaying said improved section depicting amplitude change as a function of time and centerpoint coordinate so that progressive change in said high intensity amplitude event in said displayed traces is identified as a function of progressive change in centerpoint values.

Claim 13

A method for converting an original multitrace seismic record into an improved section having increased capability as to the fluid hydrocarbon-bearing potential and/or lithologic nature of high-intensity amplitude events related to reflections from subsurface strata, said improved section being composed of a plurality of amplitude-versus-centerpoint-and-time traces, said original record consisting of a plurality of multitrace seismic traces of amplitude-versus-horizontal coordinate-and-time generated by conventional centerpoint seismic collection and processing techniques in which the presence of refraction energy within said records has been minimized, each of said traces constituting energy derived in association with a particular source-detector pair of known horizontal offset and of known centerpoint location, and representing, in part, event reflections from said subsurface strata, said conversion comprising the steps of: (a) classifying said original traces on the basis of common but progressively changing horizontal offset values and common centerpoint locations, whereby each resulting trace is identified by a centerpoint location common to at least another trace and a known horizontal offset value; (b) for a high-intensity event, determining with regard to a series of analytic functions of known mathematical character, that one of said series is the best fit to amplitude vs. time variations of said each resulting trace and said at least another trace; (c) predicting near and far amplitude vs. time trace projections for said resulting and another traces at new offset locations based on said one, best fitting analytic function; (d) generating a first series of trace amplitude projections vs. time of step (c) associated with near offset locations, and a second series of trace amplitude projections vs. time associated with far offset locations and subtracting increments of said first and second series normalized to the same time sample, one from the other to form at least a segment of said improved section; (e) weighting said subtracted increments of step (d) by semblance coefficients wherein after ratios of output-to-input energy is generated by stacking, events in said increments associated with a limited number of phases are better indicated; (f) displaying said improved section depicting amplitude change as a function of time and centerpoint coordinate so that progressive change in said high intensity amplitude event in said displayed traces is identified as a function of progressive change in centerpoint values.