US20110158046A1 - Seismic data acquisition and processing quality control - Google Patents
Seismic data acquisition and processing quality control Download PDFInfo
- Publication number
- US20110158046A1 US20110158046A1 US12/649,869 US64986909A US2011158046A1 US 20110158046 A1 US20110158046 A1 US 20110158046A1 US 64986909 A US64986909 A US 64986909A US 2011158046 A1 US2011158046 A1 US 2011158046A1
- Authority
- US
- United States
- Prior art keywords
- error attribute
- residua
- computing
- decomposed
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
- G01V2200/10—Miscellaneous details
- G01V2200/14—Quality control
Definitions
- the present application relates generally to the field of seismic data acquisition and processing, and particularly to using surface consistent decomposition for seismic data acquisition and processing quality control.
- Geophysical prospecting has been used extensively in the search for underground resources such as oil, gas, and minerals.
- Common techniques used for exploration include seismic, gravity, magnetic, and electrical methods.
- Seismic is historically the most widely used and can be subcategorized into seismic reflection and seismic refraction methods.
- the seismic reflection method the structure of subsurface formations is mapped by measuring the times required for a seismic wave, generated in the earth by a near-surface explosion, mechanical impact, vibration, or air gun, for example, to return to the surface after reflection from interfaces between formations having different physical properties.
- the reflections are recorded by detecting instruments responsive to ground motion or pressure waves.
- reflection methods one can locate and map, for example, such features as anticlines, faults, salt domes, and reefs.
- the recorded data generally are processed using computers prior to being interpreted.
- the basic objective of seismic processing is to convert the information recorded in the field into a form that best facilitates geological interpretation.
- the field data are transformed into corrected record sections.
- One object of the processing is to eliminate or reduce noise. Another is to present the reflections with the greatest possible resolution.
- Digital filtering is one commonly used processing method.
- One type of digital filtering is known as deconvolution. Because the reflected signal is effectively filtered or convolved by the earth as it passes from the source to the receiver, an inverse filter is produced that ideally cancels the effect of the earth's filtering. That is, the deconvolution operation converts the waveform of a reflection modified by the filtering of the earth into a simple pulse representing the reflection waveform before the filtering took place.
- Seismic reflection data have historically been deconvolved using single-channel spiking deconvolution. That technique, however, is based on the false assumption that the reflectivity series are “white”. That is, the deconvolution filter is designed to produce a reflection signal that is a spike, meaning all frequency components of the reflection signal's frequency spectrum have equal amplitudes. A spike would provide the highest resolution possible. However, absorption of higher frequencies by the earth broadens the source signal so that it is no longer white. Single-channel spiking deconvolution is also sensitive to random noise because any deconvolution program that raises the high frequency response could magnify existing noise in the record.
- Surface consistent deconvolution is a multichannel deconvolution using a geometric mean and gathers of common shots, common receivers, common midpoints, and common offsets. This process provides a deconvolution operator for every shot and receiver gather, corresponding respectively to the source function and the receiver impulse response, as well as an estimate of the reflectivity functions at every common midpoint (CMP) location.
- CMP common midpoint
- the present invention relates to producing a quality control measure for use during data acquisition and/or data processing of, preferably, seismic data.
- a surface consistent decomposition of the data is performed. From the surface consistent decomposition, one may compute a decomposed logarithmic spectra, and from the decomposed logarithmic spectra, one may compute one or more residua.
- An error attribute based on the one or more residua can be formulated, analyzed, and output. The error attribute can be used as a quality control measure or the analysis result can be used to produce a quality control measure.
- FIG. 1A is a scatter plot of an error attribute based on the source signal strength in accordance with the present disclosure.
- FIG. 1B is a superposition of the scatter plot of FIG. 1A and the geological map of FIG. 1C .
- FIG. 1C is a geological map corresponding to the location where the data used to produce the scatter plot of FIG. 1A were obtained.
- FIG. 2 is a flowchart showing the steps in one embodiment described in the present disclosure.
- a seismic trace P ij (t) created by a shot number i and recorded at receiver j can be written as:
- the offset operator is introduced to account for the nonlinearity of Q and its approximated estimation.
- the offset operator helps make the source and receiver operators in equation (4) physically reasonable, despite the nonlinearity of the CMP processing.
- Equation (4) Estimating the source, receiver, offset, and CMP terms (S i , R j , H p , Y k ) from the seismic data using equation (4) is known as “surface consistent decomposition”, and removing the source, receiver, and offset operators from the data is called “surface consistent deconvolution”.
- the goal of surface consistent deconvolution is to exploit data redundancy to obtain more robust statistical estimates of the wavelets.
- Equation (4) can be further modified to include a global response term, g, which represents a geometrically averaged amplitude spectrum for all traces (seismograms):
- This time domain function can be transformed into the frequency domain using the Fourier transform. Transforming to the frequency domain converts the convolution operators to multiplication operators, as is well known in the art. Surface consistent deconvolution assumes the power spectrum of the trace can be decomposed as the product of the spectra of component wavelets. This may be represented as:
- Equation (7) The decomposed logarithmic spectra used in the surface consistent spectral analysis can be used to create a quality control measure. Because in practice the right-hand side of Equation (7) will not in fact always equal the left-hand side, there are residua given by:
- R ( s,r,c,o ) L ( src )+ L ( rcv )+ L ( cmp )+ L (offset)+ Lg ⁇ L (trace). (8)
- An error attribute can be defined based on those residua. For example, the RMS (root-mean-square) of all residua for a common source could be used. Similarly, the RMS of all residua for the other terms could be used. Other quantities or desired statistics computed from the residua could also be used to define the error attribute.
- L 1 /L 2 decomposition method One method used to solve the linear system of Equation (7) is an L 1 /L 2 decomposition method. (See Kirchheimer, F., Ferber R., 2001: Robust Surface Consistent Deconvolution With Noise Suppression, 71st Ann. Meeting, SEG, San Antonio, Tex.) This method is robust and relatively immune from the effects of outliers produced by inconsistent noise. Because the L 1 /L 2 decomposition is largely unaffected by those outliers, they are visible in the residua and the error attribute can be a measure of the overall fit of the surface consistent model to those traces arising, say, from a common source. If the residua for other terms are used, the error attributes would similarly measure the overall fit to those traces sharing that common parameter.
- the error attributes are preferably written to the deconvolution database and thus can be retrieved, analyzed, and extracted exactly as the regular (spectral) attributes by a special processing module.
- a special processing module By analyzing and extracting features from the decomposed logarithmic amplitude spectra created by the surface consistent deconvolution process, one may obtain scatter maps of the source, detector, etc. stations, color coded by the feature in question.
- the error attributes and the application of spatial filtering one may obtain a multitude of possible representations which either indicate acquisition or data processing issues or correlate with known geological or cultural features of the near-surface.
- Some of the various data acquisition problems that can be identified by the error attribute include: (1) coupling problems, such as variation of sensor coupling with time; (2) additive noise (e.g., electrical noise) that shows up in a single sensor or otherwise locally within a spread; (3) ground roll that is not uniform within an analysis window for all traces of a common shot gather; and (4) excessive air blast noise that becomes visible only in some analysis gates.
- High values of the error attributes may indicate groups or even individual sensors with non-surface consistent noise.
- Dedicated hardware may receive the data, perform the data quality check using the present method, and output an alarm or other indication that a potential problem has been detected.
- the surface consistent spectral analysis can be updated after a certain number of seismic shots have been recorded to monitor the performance of groups of seismic sensors or even individual sensors. In the extreme case, the surface consistent spectra can be updated after each shot.
- the method can also be used for data quality assurance in later stages of data processing. For example, the method can be performed subsequent to a data pre-conditioning process to alert the geophysicist or data processor of potential processing-related problems.
- the error attributes can be analyzed using a special processing module.
- This analysis module may perform the following exemplary operations, though other operations are also possible.
- a scalar value from a frequency range within the spectra may be extracted using a variety of norm functions. Typically a simulated power criterion is used as norm, but RMS, median, and extreme values can be used as well.
- One may also compute the difference of two frequency ranges, not necessarily from the same attribute or analysis window, either sample by sample (then extracting the scalar from the difference vector) or after extraction. This lends itself to relative analyses.
- Spatial filters (highpass or lowpass) may be applied to the fields of extracted scalar values. Under the assumption that geological features normally exhibit smooth variation in space, as opposed to acquisition related effects, the spatial filtering is useful to selectively enhance either of those in the displays.
- the output of this module can be processed by a multitude of programs that create color coded scatter maps or plots.
- the regular spectra (optionally with spatial high pass filtering) show irregularly strong or weak stations, or indicate specific types of noise, if a characteristic frequency band is extracted or a difference of two such bands is displayed. There is also a strong correlation with the locations of roads and similar shallow underground structures.
- FIGS. 1A , 1 B, and 1 C a data acquisition quality control error attribute for the source signal strength (spatially low-pass filtered source power) superimposed on a geological map. This attribute shows the extent of source variation during the seismic survey.
- FIG. 1A shows a scatter plot of an error attribute based on the source signal strength at various source locations
- FIG. 1C shows a geological map corresponding to the location where the data used to produce the scatter plot of FIG. 1A were obtained.
- FIG. 1B is a superposition of those figures, and shows a correlation between the error attribute and the geology.
- FIG. 2 shows steps that are taken according to one embodiment of the present method.
- One begins by obtaining new or existing data having reflectivity information ( 200 ).
- a surface consistent decomposition on the data is performed ( 202 ).
- a decomposed logarithmic spectra from the surface consistent decomposition is computed ( 204 ), and one or more residua using the decomposed logarithmic spectra are computed ( 206 ).
- An error attribute is computed based on the one or more residua ( 208 ).
- the error attribute is output ( 210 ) and used as a quality control measure ( 212 ).
- the method can be used for both land and marine data.
Abstract
Description
- Not applicable.
- 1. Field of the Invention
- The present application relates generally to the field of seismic data acquisition and processing, and particularly to using surface consistent decomposition for seismic data acquisition and processing quality control.
- 2. Background
- Geophysical prospecting has been used extensively in the search for underground resources such as oil, gas, and minerals. Common techniques used for exploration include seismic, gravity, magnetic, and electrical methods. Seismic is historically the most widely used and can be subcategorized into seismic reflection and seismic refraction methods. With the seismic reflection method, the structure of subsurface formations is mapped by measuring the times required for a seismic wave, generated in the earth by a near-surface explosion, mechanical impact, vibration, or air gun, for example, to return to the surface after reflection from interfaces between formations having different physical properties. The reflections are recorded by detecting instruments responsive to ground motion or pressure waves. With reflection methods, one can locate and map, for example, such features as anticlines, faults, salt domes, and reefs.
- The recorded data generally are processed using computers prior to being interpreted. The basic objective of seismic processing is to convert the information recorded in the field into a form that best facilitates geological interpretation. The field data are transformed into corrected record sections. One object of the processing is to eliminate or reduce noise. Another is to present the reflections with the greatest possible resolution.
- Digital filtering is one commonly used processing method. One type of digital filtering is known as deconvolution. Because the reflected signal is effectively filtered or convolved by the earth as it passes from the source to the receiver, an inverse filter is produced that ideally cancels the effect of the earth's filtering. That is, the deconvolution operation converts the waveform of a reflection modified by the filtering of the earth into a simple pulse representing the reflection waveform before the filtering took place.
- Seismic reflection data have historically been deconvolved using single-channel spiking deconvolution. That technique, however, is based on the false assumption that the reflectivity series are “white”. That is, the deconvolution filter is designed to produce a reflection signal that is a spike, meaning all frequency components of the reflection signal's frequency spectrum have equal amplitudes. A spike would provide the highest resolution possible. However, absorption of higher frequencies by the earth broadens the source signal so that it is no longer white. Single-channel spiking deconvolution is also sensitive to random noise because any deconvolution program that raises the high frequency response could magnify existing noise in the record.
- Surface consistent deconvolution is a multichannel deconvolution using a geometric mean and gathers of common shots, common receivers, common midpoints, and common offsets. This process provides a deconvolution operator for every shot and receiver gather, corresponding respectively to the source function and the receiver impulse response, as well as an estimate of the reflectivity functions at every common midpoint (CMP) location.
- The present invention relates to producing a quality control measure for use during data acquisition and/or data processing of, preferably, seismic data. While or after obtaining the data, a surface consistent decomposition of the data is performed. From the surface consistent decomposition, one may compute a decomposed logarithmic spectra, and from the decomposed logarithmic spectra, one may compute one or more residua. An error attribute based on the one or more residua can be formulated, analyzed, and output. The error attribute can be used as a quality control measure or the analysis result can be used to produce a quality control measure.
-
FIG. 1A is a scatter plot of an error attribute based on the source signal strength in accordance with the present disclosure. -
FIG. 1B is a superposition of the scatter plot ofFIG. 1A and the geological map ofFIG. 1C . -
FIG. 1C is a geological map corresponding to the location where the data used to produce the scatter plot ofFIG. 1A were obtained. -
FIG. 2 is a flowchart showing the steps in one embodiment described in the present disclosure. - It is to be understood that the drawings are to be used for the purpose of illustration only, and not as a definition of the metes and bounds of the invention, the scope of which is to be determined only by the scope of the appended claims.
- Specific embodiments of the invention will now be described with reference to the figures. Like elements in the various figures will be referenced with like numbers for consistency. In the following description, numerous details are set forth to provide an understanding of the present invention. However, it will be understood by those skilled in the art that the present invention may be practiced without many of these details and that numerous variations or modifications from the described embodiments are possible.
- A seismic trace Pij(t) created by a shot number i and recorded at receiver j can be written as:
-
P ij(t)=S i(t)*R j(t)*Γij(t), (1a) - where the * represents convolution in the time domain. S is a source function, R is the receiver impulse response, and Γ is the Green's function. The goal of deconvolution is to remove the effects of the wavelet Wij=Si*Rj(t) from this seismic trace.
- For more convenient notation, the data can be sorted by common midpoint (CMP) and offset. If k is the CMP number (k=i+j when shot spacing is equal to receiver spacing) and p is the offset index (p=i+j), we have:
-
P ij(t)=S i(t)*R j(t)*Γkp(t). (1b) - Assuming CMP gathering is valid, meaning all traces gathered at a particular CMP location contain essentially the same subsurface information, this information represents the normal incidence reflectivity. If Q represents a processing sequence that approximately transforms each trace into its corresponding normal incidence reflectivity, then
-
Q(Γkp(t))≈Yk(t), (2) - where Yk(t) is the normal incidence reflectivity at CMP position k. We now have
-
Q(P ij(t))=S i(t)*R j(t)*Y k(t). (3) - An additional term, Hp(t), representing an offset dependent operator may be included in equation (3) to yield:
-
Q(P ij(t))=S i(t)*R j(t)*Y k(t)*H p(t) (4) - The offset operator is introduced to account for the nonlinearity of Q and its approximated estimation. The offset operator helps make the source and receiver operators in equation (4) physically reasonable, despite the nonlinearity of the CMP processing.
- Estimating the source, receiver, offset, and CMP terms (Si, Rj, Hp, Yk) from the seismic data using equation (4) is known as “surface consistent decomposition”, and removing the source, receiver, and offset operators from the data is called “surface consistent deconvolution”. The goal of surface consistent deconvolution is to exploit data redundancy to obtain more robust statistical estimates of the wavelets.
- Equation (4) can be further modified to include a global response term, g, which represents a geometrically averaged amplitude spectrum for all traces (seismograms):
-
Q(P ij(t))=S i(t)*R j(t)*Y k(t)*H p(t)*g. (5) - This time domain function can be transformed into the frequency domain using the Fourier transform. Transforming to the frequency domain converts the convolution operators to multiplication operators, as is well known in the art. Surface consistent deconvolution assumes the power spectrum of the trace can be decomposed as the product of the spectra of component wavelets. This may be represented as:
-
P(trace)=P(src)×P(rcv)×P(cmp)×P(offset)×Pg. (6) - Assuming all factors are positive and taking the logarithm of both sides of Equation (6) converts the product to a sum:
-
L(trace)=L(src)+L(rcv)+L(cmp)+L(offset)+Lg. (7) - Thus, the equation reduces to a linear problem for the logarithmic power spectra.
- The decomposed logarithmic spectra used in the surface consistent spectral analysis can be used to create a quality control measure. Because in practice the right-hand side of Equation (7) will not in fact always equal the left-hand side, there are residua given by:
-
R(s,r,c,o)=L(src)+L(rcv)+L(cmp)+L(offset)+Lg−L(trace). (8) - An error attribute can be defined based on those residua. For example, the RMS (root-mean-square) of all residua for a common source could be used. Similarly, the RMS of all residua for the other terms could be used. Other quantities or desired statistics computed from the residua could also be used to define the error attribute.
- One method used to solve the linear system of Equation (7) is an L1/L2 decomposition method. (See Kirchheimer, F., Ferber R., 2001: Robust Surface Consistent Deconvolution With Noise Suppression, 71st Ann. Meeting, SEG, San Antonio, Tex.) This method is robust and relatively immune from the effects of outliers produced by inconsistent noise. Because the L1/L2 decomposition is largely unaffected by those outliers, they are visible in the residua and the error attribute can be a measure of the overall fit of the surface consistent model to those traces arising, say, from a common source. If the residua for other terms are used, the error attributes would similarly measure the overall fit to those traces sharing that common parameter.
- The error attributes are preferably written to the deconvolution database and thus can be retrieved, analyzed, and extracted exactly as the regular (spectral) attributes by a special processing module. By analyzing and extracting features from the decomposed logarithmic amplitude spectra created by the surface consistent deconvolution process, one may obtain scatter maps of the source, detector, etc. stations, color coded by the feature in question. By including the error attributes and the application of spatial filtering, one may obtain a multitude of possible representations which either indicate acquisition or data processing issues or correlate with known geological or cultural features of the near-surface.
- Some of the various data acquisition problems that can be identified by the error attribute include: (1) coupling problems, such as variation of sensor coupling with time; (2) additive noise (e.g., electrical noise) that shows up in a single sensor or otherwise locally within a spread; (3) ground roll that is not uniform within an analysis window for all traces of a common shot gather; and (4) excessive air blast noise that becomes visible only in some analysis gates. High values of the error attributes may indicate groups or even individual sensors with non-surface consistent noise. Dedicated hardware may receive the data, perform the data quality check using the present method, and output an alarm or other indication that a potential problem has been detected.
- For data acquisition quality control, the surface consistent spectral analysis can be updated after a certain number of seismic shots have been recorded to monitor the performance of groups of seismic sensors or even individual sensors. In the extreme case, the surface consistent spectra can be updated after each shot. Performing the present method concurrently with data acquisition allows identification of acquisition-related problems that might otherwise go undiscovered for weeks, or even months. The method can also be used for data quality assurance in later stages of data processing. For example, the method can be performed subsequent to a data pre-conditioning process to alert the geophysicist or data processor of potential processing-related problems.
- As stated above, the error attributes can be analyzed using a special processing module. This analysis module may perform the following exemplary operations, though other operations are also possible. A scalar value from a frequency range within the spectra may be extracted using a variety of norm functions. Typically a simulated power criterion is used as norm, but RMS, median, and extreme values can be used as well. One may also compute the difference of two frequency ranges, not necessarily from the same attribute or analysis window, either sample by sample (then extracting the scalar from the difference vector) or after extraction. This lends itself to relative analyses. Spatial filters (highpass or lowpass) may be applied to the fields of extracted scalar values. Under the assumption that geological features normally exhibit smooth variation in space, as opposed to acquisition related effects, the spatial filtering is useful to selectively enhance either of those in the displays. The output of this module can be processed by a multitude of programs that create color coded scatter maps or plots.
- The regular spectra (optionally with spatial high pass filtering) show irregularly strong or weak stations, or indicate specific types of noise, if a characteristic frequency band is extracted or a difference of two such bands is displayed. There is also a strong correlation with the locations of roads and similar shallow underground structures.
- Another application is the integration of the frequency and error analysis with geographical information and satellite imaging. In examples studied, the low pass filtered displays of the overall power of all terms exhibit a very good correlation with the geology (including some deeper features), and the high pass versions of the detectors show strong correlations to cultural features such as roads. As an example we show in
FIGS. 1A , 1B, and 1C a data acquisition quality control error attribute for the source signal strength (spatially low-pass filtered source power) superimposed on a geological map. This attribute shows the extent of source variation during the seismic survey. Specifically,FIG. 1A shows a scatter plot of an error attribute based on the source signal strength at various source locations, andFIG. 1C shows a geological map corresponding to the location where the data used to produce the scatter plot ofFIG. 1A were obtained.FIG. 1B is a superposition of those figures, and shows a correlation between the error attribute and the geology. -
FIG. 2 shows steps that are taken according to one embodiment of the present method. One begins by obtaining new or existing data having reflectivity information (200). A surface consistent decomposition on the data is performed (202). A decomposed logarithmic spectra from the surface consistent decomposition is computed (204), and one or more residua using the decomposed logarithmic spectra are computed (206). An error attribute is computed based on the one or more residua (208). The error attribute is output (210) and used as a quality control measure (212). The method can be used for both land and marine data. - While the invention has been described with respect to a limited number of embodiments, those skilled in the art, having benefit of this disclosure, will appreciate that other embodiments can be envisioned that do not depart from the scope of the invention as disclosed herein. Accordingly, the scope of the invention shall be limited only by the attached claims.
Claims (28)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/649,869 US8441891B2 (en) | 2009-12-30 | 2009-12-30 | Seismic data acquisition and processing quality control |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/649,869 US8441891B2 (en) | 2009-12-30 | 2009-12-30 | Seismic data acquisition and processing quality control |
Publications (2)
Publication Number | Publication Date |
---|---|
US20110158046A1 true US20110158046A1 (en) | 2011-06-30 |
US8441891B2 US8441891B2 (en) | 2013-05-14 |
Family
ID=44187432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/649,869 Expired - Fee Related US8441891B2 (en) | 2009-12-30 | 2009-12-30 | Seismic data acquisition and processing quality control |
Country Status (1)
Country | Link |
---|---|
US (1) | US8441891B2 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102967878A (en) * | 2012-11-08 | 2013-03-13 | 中国石油天然气股份有限公司 | Method and device for predicting oil and gas reservoirs through passive seismic wave field |
CN103364831A (en) * | 2013-07-29 | 2013-10-23 | 中国石油集团川庆钻探工程有限公司 | Physical property parameter quantification method based on neural network algorithm |
CN104765063A (en) * | 2015-03-13 | 2015-07-08 | 中国石油天然气股份有限公司 | Oil-gas detection method and device for calculating absorption attenuation attribute based on frequency spectrum |
CN105425288A (en) * | 2015-12-10 | 2016-03-23 | 哈尔滨工业大学 | Collapse response spectra based on displacement evaluating indicators |
US10176407B1 (en) * | 2015-10-02 | 2019-01-08 | Hrl Laboratories, Llc | Nonlinear sparse representation-based classification for foveated analysis of spectral data with distortions |
CN111596350A (en) * | 2020-04-20 | 2020-08-28 | 江苏省地震局 | Seismic station network waveform data quality monitoring method and device |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103592683B (en) * | 2013-10-17 | 2017-01-04 | 中国石油天然气集团公司 | A kind of earth surface consistency vibration amplitude compensation method for converted wave seismic data and device |
CN104656141B (en) * | 2015-03-12 | 2016-01-20 | 中国矿业大学(北京) | Attitude road collection earth surface consistency vibration amplitude compensation method and device altogether |
US10481287B2 (en) | 2016-08-05 | 2019-11-19 | Saudi Arabian Oil Company | Surface consistent statics solution and amplification correction |
US10852450B2 (en) | 2017-05-03 | 2020-12-01 | Saudi Arabian Oil Company | Refraction-based surface-consistent amplitude compensation and deconvolution |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5144590A (en) * | 1991-08-08 | 1992-09-01 | B P America, Inc. | Bed continuity detection and analysis using crosswell seismic data |
US5491669A (en) * | 1994-10-06 | 1996-02-13 | Western Atlas International | Method for identifying and suppressing anomalous frequency components in seismic data |
US5671136A (en) * | 1995-12-11 | 1997-09-23 | Willhoit, Jr.; Louis E. | Process for seismic imaging measurement and evaluation of three-dimensional subterranean common-impedance objects |
US5771203A (en) * | 1994-09-23 | 1998-06-23 | Compagnie Generale De Geophysique | Seismic prospection method with application of a self-deconvoluted prediction error filter |
US5774417A (en) * | 1996-10-25 | 1998-06-30 | Atlantic Richfield Company | Amplitude and phase compensation in dual-sensor ocean bottom cable seismic data processing |
US20030200030A1 (en) * | 1998-09-11 | 2003-10-23 | Den Norske Stats Oljeselskap A.S. | Method of processing acquired seismic data |
US6772068B2 (en) * | 2000-06-30 | 2004-08-03 | Andreas Laake | Quality control of data |
US20080294393A1 (en) * | 2007-05-24 | 2008-11-27 | Laake Andreas W | Near Surface Layer Modeling |
US20100292973A1 (en) * | 2008-01-21 | 2010-11-18 | Arkex Limited | Geophysical data processing systems |
US8139440B2 (en) * | 2007-12-27 | 2012-03-20 | Westerngeco L.L.C. | Spectral conditioning for surface seismic data |
-
2009
- 2009-12-30 US US12/649,869 patent/US8441891B2/en not_active Expired - Fee Related
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5144590A (en) * | 1991-08-08 | 1992-09-01 | B P America, Inc. | Bed continuity detection and analysis using crosswell seismic data |
US5771203A (en) * | 1994-09-23 | 1998-06-23 | Compagnie Generale De Geophysique | Seismic prospection method with application of a self-deconvoluted prediction error filter |
US5491669A (en) * | 1994-10-06 | 1996-02-13 | Western Atlas International | Method for identifying and suppressing anomalous frequency components in seismic data |
US5671136A (en) * | 1995-12-11 | 1997-09-23 | Willhoit, Jr.; Louis E. | Process for seismic imaging measurement and evaluation of three-dimensional subterranean common-impedance objects |
US5774417A (en) * | 1996-10-25 | 1998-06-30 | Atlantic Richfield Company | Amplitude and phase compensation in dual-sensor ocean bottom cable seismic data processing |
US20030200030A1 (en) * | 1998-09-11 | 2003-10-23 | Den Norske Stats Oljeselskap A.S. | Method of processing acquired seismic data |
US6772068B2 (en) * | 2000-06-30 | 2004-08-03 | Andreas Laake | Quality control of data |
US20080294393A1 (en) * | 2007-05-24 | 2008-11-27 | Laake Andreas W | Near Surface Layer Modeling |
US8139440B2 (en) * | 2007-12-27 | 2012-03-20 | Westerngeco L.L.C. | Spectral conditioning for surface seismic data |
US20100292973A1 (en) * | 2008-01-21 | 2010-11-18 | Arkex Limited | Geophysical data processing systems |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102967878A (en) * | 2012-11-08 | 2013-03-13 | 中国石油天然气股份有限公司 | Method and device for predicting oil and gas reservoirs through passive seismic wave field |
CN103364831A (en) * | 2013-07-29 | 2013-10-23 | 中国石油集团川庆钻探工程有限公司 | Physical property parameter quantification method based on neural network algorithm |
CN104765063A (en) * | 2015-03-13 | 2015-07-08 | 中国石油天然气股份有限公司 | Oil-gas detection method and device for calculating absorption attenuation attribute based on frequency spectrum |
US10176407B1 (en) * | 2015-10-02 | 2019-01-08 | Hrl Laboratories, Llc | Nonlinear sparse representation-based classification for foveated analysis of spectral data with distortions |
CN105425288A (en) * | 2015-12-10 | 2016-03-23 | 哈尔滨工业大学 | Collapse response spectra based on displacement evaluating indicators |
CN111596350A (en) * | 2020-04-20 | 2020-08-28 | 江苏省地震局 | Seismic station network waveform data quality monitoring method and device |
Also Published As
Publication number | Publication date |
---|---|
US8441891B2 (en) | 2013-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8441891B2 (en) | Seismic data acquisition and processing quality control | |
US7953556B2 (en) | Geophone noise attenuation and wavefield separation using a multi-dimensional decomposition technique | |
Tselentis et al. | Strategy for automated analysis of passive microseismic data based on S-transform, Otsu’s thresholding, and higher order statistics | |
US11880011B2 (en) | Surface wave prediction and removal from seismic data | |
US8315124B2 (en) | System and method for suppression of seismic multiple reflection signals | |
GB2479347A (en) | A process for characterising the evolution of a reservoir | |
US7616524B1 (en) | Wavelet based intercept attribute for seismic exploration | |
CA2855734A1 (en) | Coherent noise attenuation | |
US9234977B2 (en) | Processing collected survey data | |
US20140078860A1 (en) | Interference noise attenuation method and apparatus | |
MX2010005019A (en) | Method for calculation of seismic attributes from seismic signals. | |
Dai et al. | Effects due to aliasing on surface-wave extraction and suppression in frequency-velocity domain | |
AU2014101578A4 (en) | System and method for seismic adaptive optics | |
CN113805237B (en) | Method and system for offset land crossrange seismic using compressed sensing models | |
Lambert et al. | Spectral analysis of ambient ground-motion—Noise reduction techniques and a methodology for mapping horizontal inhomogeneity | |
Cheng et al. | Q-estimation using seismic interferometry from vertical well data | |
Pevzner et al. | Seismic attenuation from VSP and well log data: approaches, problems and relative contribution of scattering | |
Strobbia et al. | Surface waves–Use them then lose them | |
GB2590177A (en) | Methods and devices performing adaptive subtraction of multiples | |
GB2375606A (en) | Angle dependent surface multiple attenuation for two - component marine bottom sensor data | |
Behura | Virtual real source | |
Massin et al. | A better automatic body-wave picker with broad applicability | |
Chopra et al. | Volumetric curvature-attribute applications for detection of fracture lineaments and their calibration | |
Shang et al. | A bandwidth enhancement workflow through wavelet analysis | |
Orji et al. | A Comparative Study of Thin-Bed Interpretation using Spectral and Cepstral Transform Techniques in Dense 3D Seismic Amplitude Data in Niger Delta |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: WESTERNGECO L. L. C., TEXAS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:FERBER, RALF;KIRCHHEIMER, FRANZ;LAAKE, ANDREAS;SIGNING DATES FROM 20090223 TO 20091105;REEL/FRAME:023718/0418 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
LAPS | Lapse for failure to pay maintenance fees |
Free format text: PATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20210514 |