Patent US4760563A successfully refreshed from Google Patents.
US4760563A
Expired / Lapsed
US
⏱️ 0.0 years remaining
Seismic exploration using exactly invertible discrete transformation into tau-p space
Assignee: Schlumberger Technology Corp
•
Filing: 1986-01-09
•
Granted: 1988-07-26
•
Expires: 2006-01-09
Abstract
Abstract Disclosed are a method and a system for discrete transformation of measurements such as seismic data into and out of tau-p space which is both exact and practical in terms of processing time. The measurements can be filtered or otherwise processed in tau-p space in ways which are not practical or possible in their original space. Since the transformations into and out of tau-p space are exact, the filtered and transformed measurements are free of errors and distortions that perturb known approximate transforms which can be performed in a reasonable time. When the transformation process is carried out in frequency space, it is done frequency-by-frequency, and when it is carried out in the spatial domain, it can utilize a transformation matrix having a block-circulant structure. In each case, the transformation process and matrix have a structure which substantially reduces storage and processing requirements as compared with known prior art.
Inventors
- 👤 Gregory Beylkin
Bibliographic Details
Application Number
—
Jurisdiction
US
Data Source
google_patents
Last Synchronized
2026-10-27 06:07
CPC Classifications
G; G01; G01V; G01V1/00; G01V1/28; G01V1/36; G01V1/364; G01V1/32; G01V2210/00; G01V2210/20; G01V2210/40; G01V2210/47
Claims Analysis
Total Claims: 53
Claim 1
Independent Claim
96 words
A method of seismic exploration of subsurface formations comprising: deriving a first seismogram of selected subsurface formations; transforming said first seismogram to tau-p space through a discrete transformation process which permits exact discrete inverse-transformation and which in processing in the frequency domain can be completed frequency-by-frequency and in processing in the spatial domain can utilize a transformation matrix having a block-circulant structure; processing the tau-p space seismogram to enhance desirable features thereof and effect an exact discrete inverse transformation from tau-p space to generate an improved seismogram; and providing a tangible representation of the improved seismogram.
Claim 2
Independent Claim
96 words
A method of seismic exploration of subsurface formations comprising: deriving a first seismogram of selected subsurface formations; transforming said first seismogram to tau-p space through a discrete transformation process which permits exact discrete inverse-transformation and which in processing in the frequency domain can be completed frequency-by-frequency and in processing in the spatial domain can utilize a transformation matrix having a block-circulant structure; processing the tau-p space seismogram to enhance desirable features thereof and effect an exact discrete inverse transformation from tau-p space to generate an improved seismogram; and providing a tangible representation of the improved seismogram.
Claim 3
Dependent Claim
47 words
A method as in claim 1 in which said processing comprises processing the tau-p space seismogram to enhance desirable features thereof and thereby generate an improved tau-p space seismogram, and inverse-transforming said improved tau-p space seismogram using an exact discrete inverse transformation to generate said improved seismogram.
Claim 4
Dependent Claim
20 words
A method as in claim 2 including the step of providing a tangible representation of said improved tau-p space seismogram.
Claim 5
Dependent Claim
18 words
A method as in claim 1 including the step of performing a spectrum shift upon said first seismogram.
Claim 6
Dependent Claim
16 words
A method as in claim 4 in which said spectrum shift is incorporated in said transformation.
Claim 7
Dependent Claim
24 words
A method as in claim 1 in which the tangible representation of the improved seismogram is in the form of a visible display thereof.
Claim 8
Dependent Claim
37 words
A method as in claim 1 in which the inverse transformation transforms the improved tau-p space seismogram from tau-p space to the space in which said first seismogram was represented prior to its transformation into tau-p space.
Claim 9
Dependent Claim
48 words
A method as in claim 1 in which the transforming step comprises a transformation in which selected combinations of traces of said first seismogram and slopes of said tau-p space seismogram are weighted to enhance or suppress the effect of the corresponding portion of the tau-p space seismogram.
Claim 10
Dependent Claim
19 words
A method as in claim 1 in which the deriving step comprises deriving said first seismogram in time-offset space.
Claim 11
Dependent Claim
19 words
A method as in claim 1 in which the deriving step comprises deriving said first seismogram in time-depth space.
Claim 12
Dependent Claim
28 words
A method as in claim 1 in which the transforming step comprises a transformation in which selected curves of said first seismogram transform to points in tau-p space.
Claim 13
Dependent Claim
17 words
A method as in claim 11 in which said selected curves of said first seismogram are nonlinear.
Claim 14
Independent Claim
119 words
A method of seismic exploration of subsurface formations comprising: deriving a seismogram of subsurface formations which is in discrete form in a selected space; transforming the seismogram through the use of a discrete form of a transform to generate a seismogram in tau-p space, in a process which comprises utilizing a transformation matrix having the property of being block circulant when represented in the spatial domain; filtering the seismogram in tau-p space to generate a filtered seismogram; inverse-transforming the filtered seismogram from tau-p space through the use of an exact discrete form of an inverse transform to thereby generate a filtered seismogram in the selected space; and providing a tangible representation of the filtered seismogram in the selected space.
Claim 15
Independent Claim
119 words
A method of seismic exploration of subsurface formations comprising: deriving a seismogram of subsurface formations which is in discrete form in a selected space; transforming the seismogram through the use of a discrete form of a transform to generate a seismogram in tau-p space, in a process which comprises utilizing a transformation matrix having the property of being block circulant when represented in the spatial domain; filtering the seismogram in tau-p space to generate a filtered seismogram; inverse-transforming the filtered seismogram from tau-p space through the use of an exact discrete form of an inverse transform to thereby generate a filtered seismogram in the selected space; and providing a tangible representation of the filtered seismogram in the selected space.
Claim 16
Dependent Claim
19 words
A method as in claim 13 in which said tangible representation is in the form of a visible display.
Claim 17
Dependent Claim
53 words
A method as in claim 13 in which the transforming step comprises a transformation in which selected combinations of traces of the seismogram which is in said selected space and slopes of the seismogram in tau-p space are weighted to enhance or suppress the effect of the corresponding portion of the tau-p seismogram.
Claim 18
Dependent Claim
33 words
A method as in claim 13 in which the transforming step comprises a transformation in which selected curves of the seismogram which is in said selected space transform to points in tau-p space.
Claim 19
Dependent Claim
13 words
A method as in claim 16 in which said selected curves are nonlinear.
Claim 20
Independent Claim
129 words
A method of seismic exploration of subsurface formations comprising: deriving an original seismogram of selected subsurface formations which is represented as a set of traces each comprising discrete samples of a selected seismic parameter over time for a selected spatial parameter; transforming the original seismogram into tau-p space by transforming the original seismogram into the frequency domain, multiplying the frequency domain version of the original seismogram by selected transformation matrices for respective frequencies to derive respective product arrays, and inverse transforming the product arrays out of the frequency domain to derive a representation of the original seismogram in discrete form in tau-p space and allowing exact discrete inverse-transformation; and providing a tangible representation of the seismogram in tau-p space in form suitable for modifying or display in tau-p space.
Claim 21
Independent Claim
129 words
A method of seismic exploration of subsurface formations comprising: deriving an original seismogram of selected subsurface formations which is represented as a set of traces each comprising discrete samples of a selected seismic parameter over time for a selected spatial parameter; transforming the original seismogram into tau-p space by transforming the original seismogram into the frequency domain, multiplying the frequency domain version of the original seismogram by selected transformation matrices for respective frequencies to derive respective product arrays, and inverse transforming the product arrays out of the frequency domain to derive a representation of the original seismogram in discrete form in tau-p space and allowing exact discrete inverse-transformation; and providing a tangible representation of the seismogram in tau-p space in form suitable for modifying or display in tau-p space.
Claim 22
Dependent Claim
22 words
A method as in claim 18 including modifying the seismogram in tau-p space to thereby derive a modified seismogram in tau-p space.
Claim 23
Dependent Claim
119 words
A method as in claim 19 including transforming the modified seismogram into the frequency domain, multiplying the frequency domain version of the modified seismogram by the adjoint matrices of said transformation matrices to derive respective adjoint product arrays, using the adjoint matrices and the transformation matrices to derive selected matrices for filtering, solving for the array a linear relationship between the array and adjoint product arrays and said matrices for filtering, and inverse-transforming the last recited array out of frequency space to thereby generate a seismogram of said selected subsurface formations which corresponds to said modified seismogram in tau-p space but is represented as a set of traces each comprising discrete samples of a selected seismic parameter over time.
Claim 24
Dependent Claim
38 words
A method as in claim 18 including weighting selected combinations of traces of the original seismogram and slopes of the seismogram in tau-p space to enhance or suppress the effect of the corresponding portion of the tau-p seismogram.
Claim 25
Dependent Claim
28 words
A method as in claim 18 in which the transforming step comprises a transformation in which selected curves of the original seismogram transform to points in tau-p space.
Claim 26
Dependent Claim
13 words
A method as in claim 22 in which said selected curves are nonlinear.
Claim 27
Independent Claim
207 words
A method comprising: deriving a two-dimensional array x l (n) of seismic traces in (l,n) space, where l is an integer index which identifies a respective one of (2L+l) traces which are rows in the array and ranges from -L to +L, L being an integer, each row of the array comprises discrete samples of a selected seismic parameter, each sample being identified by an index n=1, 2, . . . , N, and x l (n) is a periodic sequence with the period N, N being an integer; deriving transformation matrices R jl (k), where j is an integer index which identifies a respective one of (2J+l) rows of a transform matrix and ranges from -J to +J, l is an integer index which identifies a respective one of (2L+l) columns of a transformation matrix and ranges from -L to +L, and where k is a frequency index and J and L are integers; transforming the array x l (n) into the frequency domain to derive the array x l (k+1); matrix-multiplying the array x l (k+1) by the matrices R jl (k) to derive the array y j (k+1) allowing exact discrete inverse-transformation; and modifying the array y j (k+1) to derive a modified array.
Claim 28
Independent Claim
207 words
A method comprising: deriving a two-dimensional array x l (n) of seismic traces in (l,n) space, where l is an integer index which identifies a respective one of (2L+l) traces which are rows in the array and ranges from -L to +L, L being an integer, each row of the array comprises discrete samples of a selected seismic parameter, each sample being identified by an index n=1, 2, . . . , N, and x l (n) is a periodic sequence with the period N, N being an integer; deriving transformation matrices R jl (k), where j is an integer index which identifies a respective one of (2J+l) rows of a transform matrix and ranges from -J to +J, l is an integer index which identifies a respective one of (2L+l) columns of a transformation matrix and ranges from -L to +L, and where k is a frequency index and J and L are integers; transforming the array x l (n) into the frequency domain to derive the array x l (k+1); matrix-multiplying the array x l (k+1) by the matrices R jl (k) to derive the array y j (k+1) allowing exact discrete inverse-transformation; and modifying the array y j (k+1) to derive a modified array.
Claim 29
Dependent Claim
42 words
A method as in claim 24 in which the transformation matrices R jl (k) are at least in part a function of the array x l (n) and have the property of having a block-circulant structure when represented in the spatial domain.
Claim 30
Dependent Claim
55 words
A method as in claim 24 in which the modifying step comprises inverse-transforming the array y j (k+1) out of frequency space to derive the array y j (n) and to represent it as a seismogram in (j,n) space allowing exact discrete inverse-transformation, and modifying the array y j (n) to derive a modified seismogram.
Claim 31
Dependent Claim
112 words
A method as in claim 24 including: matrix-multiplying the modified array y j (k+1) with adjoints of said transformation matrices to generate an array z l (k+1); generating matrices H ll' (k) by matrix multiplying said transform and adjoint matrices, where l' is an integer index and ranges from -L to L, L being an integer; solving the systems z l (k+1)=H ll' (k)x l' (k+1) for the array x l' (k+1), wherein said systems are linear; and inverse-transforming the array x l' (k+1) out of the frequency domain and into (l',n) space to derive a modified array x l' (n) and to represent same as a modified seismogram in (l',n) space.
Claim 32
Dependent Claim
75 words
A method as in claim 21 wherein the modifying step comprises inverse-transforming the array y j (k+1) out of frequency space to derive the array y j (n) and to represent it as a seismogram in (j,n) space allowing exact discrete inverse-transformation, and modifying the array y j (n) to derive a modified seismogram, and including transforming the modified array y j (n) into the frequency domain to generate said modified array y j (k+1).
Claim 33
Dependent Claim
42 words
A method as in claim 24 including weighting selected combinations of traces of the seismogram in (l,n) space and slopes of the seismogram in (j,n) space to enhance or suppress the effect of the corresponding portion of the seismogram in (j,n) space.
Claim 34
Dependent Claim
25 words
A method as in claim 34 in which selected curves of the seismogram in (l,n) space transform to points in the seismogram in (j,n) space.
Claim 35
Dependent Claim
13 words
A method as in claim 30 in which said selected curves are nonlinear.
Claim 36
Independent Claim
108 words
A method comprising: deriving an original record of measurements of a selected parameter of energy propagating in a medium and sampled in time and in space; transforming the original record to a record in tau-p space through a discrete transformation process which permits exact discrete inverse-transformation and which in processing in the frequency domain can be completed frequency-by-frequency and in processing in the spatial domain can utilize a transformation matrix having a block-circulant structure; processing the tau-p space record to enhance desirable features thereof and effect an exact discrete inverse transformation from tau-p space to generate an improved record; and providing a tangible representation of the improved record.
Claim 37
Independent Claim
108 words
A method comprising: deriving an original record of measurements of a selected parameter of energy propagating in a medium and sampled in time and in space; transforming the original record to a record in tau-p space through a discrete transformation process which permits exact discrete inverse-transformation and which in processing in the frequency domain can be completed frequency-by-frequency and in processing in the spatial domain can utilize a transformation matrix having a block-circulant structure; processing the tau-p space record to enhance desirable features thereof and effect an exact discrete inverse transformation from tau-p space to generate an improved record; and providing a tangible representation of the improved record.
Claim 38
Dependent Claim
50 words
A method as in claim 32 in which said processing comprises processing the tau-p space record to enhance desirable features thereof and thereby generate an improved tau-p record to enhance desirable features thereof and thereby improved tau-p space record using an exact discrete inverse transformation to generate said improved record.
Claim 39
Dependent Claim
20 words
A method as in claim 33 including the step of providing a tangible representation of said improved tau-p space record.
Claim 40
Dependent Claim
18 words
A method as in claim 32 including the step of performing a spectrum shift upon said first record.
Claim 41
Dependent Claim
16 words
A method as in claim 35 in which said spectrum shift is incorporated in said transformation.
Claim 42
Dependent Claim
28 words
A method as in claim 32 in which the inverse transformation transforms the improved tau-p space record from the tau-p space to the space of the original record.
Claim 43
Dependent Claim
39 words
A method as in claim 32 including weighting selected combinations of traces of the original record and slopes of the record in tau-p space to enhance or suppress the effect of the corresponding portion of the tau-p space record.
Claim 44
Dependent Claim
31 words
A method as in claim 32 in which the transforming step comprises a transformation in which selected curves of the original record transform to points in the record in tau-p space.
Claim 45
Dependent Claim
17 words
A method as in claim 39 in which said selected curves of the original record are nonlinear.
Claim 46
Independent Claim
145 words
A method comprising: deriving an original record of data for a selected parameter sampled with respect to two independent variables defining an original domain; transforming the original record to a record in a different domain in which each point results from a sum of values of said parameter taken over a corresponding pattern in said original domain, through a discrete transformation which permits exact discrete inverse-transformation and which in processing in a frequency domain corresponding to one of said independent variables can be completed frequency-by-frequency and in processing in said original domain with respect to the other of said independent variables can utilize a transformation matrix having a block-circulant structure; processing said different domain record to enhance desirable features thereof and effect an exact discrete inverse transformation from said different domain to generate an improved record; and providing a tangible representation of the improved record.
Claim 47
Independent Claim
145 words
A method comprising: deriving an original record of data for a selected parameter sampled with respect to two independent variables defining an original domain; transforming the original record to a record in a different domain in which each point results from a sum of values of said parameter taken over a corresponding pattern in said original domain, through a discrete transformation which permits exact discrete inverse-transformation and which in processing in a frequency domain corresponding to one of said independent variables can be completed frequency-by-frequency and in processing in said original domain with respect to the other of said independent variables can utilize a transformation matrix having a block-circulant structure; processing said different domain record to enhance desirable features thereof and effect an exact discrete inverse transformation from said different domain to generate an improved record; and providing a tangible representation of the improved record.
Claim 48
Dependent Claim
47 words
A method as in claim 41 in which said processing comprises processing said different domain record to enhance desirable features thereof and thereby generate an improved different domain record, and inverse-transforming said improved different domain record using an exact discrete inverse transformation to generate said improved record.
Claim 49
Dependent Claim
20 words
A method as in claim 42 including the step of providing a tangible representation of said improved different domain record.
Claim 50
Dependent Claim
25 words
A method as in claim 41 in which the inverse transformation transforms the improved different domain record from said different domain to said original domain.
Claim 51
Dependent Claim
48 words
A method as in claim 41 in which the transforming step comprises a transformation in which selected combinations of values of said original record and slopes of said different domain record are weighted to enhance or suppress the effect of the corresponding portion of the different domain record.
Claim 52
Dependent Claim
12 words
A method as in claim 41 in which said pattern is non-linear.
Claim 53
Dependent Claim
12 words
A method as in claim 41 in which said pattern is disjoint.
References Cited (Backward Citations)
US4562540A
Cited
USPTO Annuity & Maintenance Windows
| Annuity Stage | Window Opens | Payment Deadline | Surcharge Period Until | Estimated Fee | Status |
|---|---|---|---|---|---|
| 3.5 Years (Year 4) | 1991-07-26 | 1992-01-25 | 1992-07-26 | $2000.00 | Window Closed / Lapsed |
| 7.5 Years (Year 8) | 1995-07-26 | 1996-01-25 | 1996-07-26 | $3760.00 | Window Closed / Lapsed |
| 11.5 Years (Year 12) | 1999-07-26 | 2000-01-25 | 2000-07-26 | $7700.00 | Window Closed / Lapsed |