WO2010082938A1 - Data acquisition and prestack migration based on seismic visibility analysis - Google Patents
Data acquisition and prestack migration based on seismic visibility analysis Download PDFInfo
- Publication number
- WO2010082938A1 WO2010082938A1 PCT/US2009/031400 US2009031400W WO2010082938A1 WO 2010082938 A1 WO2010082938 A1 WO 2010082938A1 US 2009031400 W US2009031400 W US 2009031400W WO 2010082938 A1 WO2010082938 A1 WO 2010082938A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- visibility
- seismic
- target event
- receiver
- source
- Prior art date
Links
- 238000013508 migration Methods 0.000 title claims abstract description 37
- 230000005012 migration Effects 0.000 title claims abstract description 32
- 238000004458 analytical method Methods 0.000 title claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 52
- 238000003384 imaging method Methods 0.000 claims abstract description 8
- 230000006870 function Effects 0.000 claims description 18
- 230000015654 memory Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 description 12
- 238000010304 firing Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000005755 formation reaction Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000001902 propagating effect Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/003—Seismic data acquisition in general, e.g. survey design
-
- 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. for interpretation or for event detection
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/67—Wave propagation modeling
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/67—Wave propagation modeling
- G01V2210/673—Finite-element; Finite-difference
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/67—Wave propagation modeling
- G01V2210/675—Wave equation; Green's functions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/67—Wave propagation modeling
- G01V2210/679—Reverse-time modeling or coalescence modelling, i.e. starting from receivers
Definitions
- Seismic surveys can provide information about underground structures, including formation boundaries, rock types, and the presence or absence of fluid reservoirs. Such information greatly aids searches for water, geothermal reservoirs, and mineral deposits such as hydrocarbons and ores. Oil companies in particular often invest in extensive seismic surveys to select sites for exploratory oil wells.
- Conventional seismic surveys employ artificial seismic energy sources such as shot charges, air guns, or vibratory sources to generate seismic waves.
- the sources when fired, create a seismic "event", i.e., a pulse of seismic energy that propagates as seismic waves from the source down into the earth. Faults and boundaries between different formations create differences in acoustic impedance that cause partial reflections of the seismic waves.
- a seismic sensor array detects and records these reflections for later analysis. Sophisticated processing techniques are then applied to the recorded signals to extract an image or other representation of the subsurface structure.
- systems and methods for performing seismic visibility analysis of selected subsurface structures identify the seismic source and receiver positions that can best reveal the details of the subsurface structure. These positions can then be used as the basis for acquiring additional seismic data and/or subjecting a selected subset of the existing data to more sophisticated data processing. Because the region of data acquisition and/or processing is greatly reduced, the associated expenses are minimized.
- Some illustrative method embodiments include a seismic survey method that includes: determining visibility of a target event as a function of seismic source and receiver positions; and acquiring seismic data in a region selected at least in part to include positions having visibility values above a threshold. The target event can then be imaged based on the newly acquired seismic data.
- the illustrative method embodiments also include a seismic migration method that includes: determining visibility of a target event at the source and receiver positions of traces in an existing seismic survey; and re-migrating traces having visibility values above a threshold to image the target event.
- the visibility determination may include using a wave equation based propagator to find, for each of multiple simulated shots, a reflection wavefield from the target event in a seismic model; and to calculate, for each of multiple receiver positions, a contribution signal from each reflection wavefield.
- the visibility determination may further include converting each contribution signal into a source-receiver visibility value.
- Fig. 1 shows an illustrative seismic survey environment
- Fig. 2 shows an illustrative seismic source and receiver geometry
- Fig. 3 shows illustrative seismic traces
- Fig. 4 shows a migrated depth image for an illustrative seismic model
- Fig. 5a illustrates seismic energy propagating from a source to a target
- Fig. 5b shows an illustrative graph of receiver visibility
- Figs. 6a-6c show illustrative graphs of source visibility under different assumptions
- Fig. 7 shows an illustrative source- receiver visibility function
- Fig. 8 shows an illustrative migrated depth image for a migrated subset of seismic data
- Fig. 9 shows a flow diagram of an illustrative visibility analysis method
- Fig. 10 shows an illustrative visibility analysis system. While the invention is susceptible to various modifications and alternative forms, specific embodiments thereof are shown by way of example in the drawings and will herein be described in detail. It should be understood, however, that the drawings and detailed description thereto are not intended to limit the invention to the particular form disclosed, but on the contrary, the intention is to cover all modifications, equivalents and alternatives falling within the scope of the appended claims.
- Fig. 1 shows an illustrative seismic survey environment, in which an array of seismic receivers 102 are positioned in a spaced-apart arrangement on the earth's surface 104 to detect seismic waves.
- the receivers 102 are coupled wirelessly or via cable to a data acquisition unit 106 that receives, processes, and stores the seismic signal data collected by the receivers.
- a seismic energy source 108 e.g., a vibrator truck
- a seismic energy source 108 is triggered at multiple positions to generate seismic energy waves that propagate through the earth 110 and reflect from acoustic impedance discontinuities to reach the receivers 102.
- Such discontinuities may be created by faults, boundaries between formation beds, and boundaries between formation fluids. The discontinuities will appear as bright spots in the subsurface structure representation that is derived from the seismic signal data.
- Fig. 1 further shows an illustrative subsurface model that will be used as an example in this disclosure.
- the earth has four relatively flat formation layers with a steeply curved boundary between the third and fourth layers.
- the speed of sound in each of the layers from top to bottom is 2000 m/s, 3500 m/s, 2800 m/s, and 4000 m/s, respectively.
- Fig. 1 is not drawn to scale.
- the horizontal extent of the receiver array is expected to be limited to a couple of kilometers or less, while the horizontal dimension of the subsurface model is sixteen kilometers and the vertical dimension is five kilometers.
- Fig. 2 shows an illustrative geometry for source positions 202 and receiver positions 204, as they might be seen from an overhead view.
- Viable seismic survey geometries are infinite in variety, and can be gathered on a sector-by-sector basis, a rolling geometry basis, a moving-array basis, and various combinations thereof.
- the main message here is that the number of receiver signals acquired in response to each firing of the seismic source is fairly large, and when the number of source firings is taken into account, the resulting number of traces can easily reach into the millions.
- Fig. 3 shows some illustrative received seismic signal traces S1-S3.
- the traces can represent displacement, velocity, acceleration, pressure, or some other measure of seismic energy intensity as a function of time.
- the signal received by each receiver is typically sampled and digitized to between 8 and 32 bits of resolution at a rate of about 500 samples per second for a duration of about 30 seconds after each shot.
- the receivers sense multi-component data, further increasing the amount of sample data for each trace.
- the trace data may be filtered and compressed before storage.
- the stored seismic survey data is transported or otherwise communicated to a data processing facility.
- a network of computers at the data processing facility processes the data to estimate the volumetric distribution of sound velocities using known techniques. See, e.g., Jon F. Claerbout, Fundamentals of Geophysical Data Processing, p. 246- 56, which is hereby incorporated herein by reference.
- the velocity distribution may be independently available from other sources, e.g., well logs.
- the data processing facility "migrates" the seismic traces, translating them from functions of time to functions of depth.
- Various migration techniques exist, including ones based on the one-way wave equation migration (one-way WEM), and full-way wave equation based reverse-time migration (RTM).
- One-way WEM is a popular, widely applied technique because it is effective in many cases and is relatively inexpensive in terms of computational complexity.
- one-way WEM simply fails to image the complex structures. This situation may be exacerbated in marine seismic surveys since the usage of narrow-azimuth receivers limits the amount of cross-line offset.
- RTM is able to address such imaging problems.
- RTM has become more attractive for prestack depth imaging processing in complex media and subsalt structures.
- RTM is computationally expensive and requires the data processing facility to have computers with large memories and large disk capacities. Moreover, RTM becomes even more challenging when migrating high-frequency components of the wavefield due to the numerical dispersion of the finite-difference scheme.
- Fig. 4 shows a migrated depth image 402 for the illustrative seismic model of Fig. 1.
- One hundred shots were simulated by finite-difference forward modeling. The shot position ranged from 6 km to 14 km with an 80 meter shot interval. The receiver aperture of each shot covered the whole model. The receiver interval is 10 meters.
- the depth image of Fig. 4 was obtained by migrating all shots using the full aperture. In general, the curved event is well imaged but the amplitude of the steep dip event 404 is weak. To get a better image of this target event, it is not necessary to obtain large amounts of new survey data or to reprocess all of the existing data using RTM. Instead, we can just focus on the seismic traces that have significant contributions to imaging the target event. The seismic visibility analysis methods and systems described below will provide a quantitative identification of which traces provide such contributions.
- Fig. 5a is a schematic ray-tracing illustration of seismic propagating seismic energy.
- Ray 502 shows energy propagating from a source to one end of the target event and back to the surface
- ray 504 shows energy propagating from the source to the other end of the target event and back to the surface.
- wave equation migration is employed to propagate the seismic energy downwards as back to the surface. More specifically, visibility analysis is preferably implemented using a wave equation based propagator rather than a high frequency asymptotic ray-based approach. The wave equation method is more accurate and provides a more reliable result.
- the visibility analysis takes place in two phases. First the wavefield of a simulated shot is propagated downward and the software measures the reflection wavefield from the target event. In the second phase, the reflection wavefield is propagated back and the software measures the target's contribution to the signals recorded by each receiver.
- the simulated shot and receiver positions can be uniformly spaced throughout the model area, or they can be customized to the contemplated survey environment (e.g., a marine streamer geometry).
- V R (r) V(S,r) .
- Fig. 5b illustrates the receiver visibility 506 of target event 404 for the source firing position shown in Fig. 5a. It can be seen that for this source position, the receiver visibility is largely limited to a well-defined neighborhood above the target event.
- V s (s) is herein defined as a summation of the source-receiver visibility V(s,r) over all receiver positions ⁇ R ⁇ :
- V 8 (S) ⁇ V(s,r) .
- Fig. 6 shows the source visibility of target event 404 under three different assumptions.
- Fig. 6a shows the source visibility with full receiver aperture (i.e., all receivers can respond to all source firings).
- Fig. 6b shows the source visibility assuming a ⁇ 2 km aperture (i.e., only the receivers within 2 kilometers of the source position can respond to the source firing).
- Fig. 6c shows the zero-offset source visibility (i.e., the receiver is co-located with the source). [1] In each case, the visibility varies with source position, and those shots in the neighborhood above and to the right of the target event contribute most to the visibility of the target event. A comparison of these source visibility functions enables the effects of receiver aperture to be readily quantified.
- Fig. 7 shows a map of the source- receiver visibility function V(s,r) for the target event 404.
- the horizontal axis denotes the receiver location, while the vertical axis denotes the source location.
- the traces for the source- receiver positions found within region 702 contribute the most to the visibility of event 404.
- the visibility function along the line 704 is the receiver visibility function V R (r) for the source position shown in Fig. 5a.
- the zero-offset source visibility function shown in Fig. 6c is found along the 45° line 708.
- the source visibility shown in Fig. 6c is found by integrating horizontally across the entire figure.
- the source visibility function shown in Fig. 6b is found by integrating horizontally between lines 706.
- the maximum receiver aperture (which corresponds to the cable length in seismic surveys) can be selected by adjusting the spacing between the lines 706 to capture the bulk of the nonzero area under the visibility function.
- the source positions can then be selected to capture the bulk of the nonzero area under the source visibility function.
- Fig. 8 shows a depth migration image using only the selected data, which demonstrates that the amplitude behavior of the visibility controlled image is more balanced.
- Fig. 9 shows an illustrative visibility analysis method that can be carried out by a computer system automatically or under the guidance of a user such as a reservoir analyst.
- the system obtains data representing a depth-migrated image. In many cases, this data will have been generated by the system itself based on a previous seismic survey.
- the system identifies one or more target events, e.g., features that have been inadequately imaged. In some implementations, the system identifies the target events by displaying the depth-migrated image to a user and soliciting input from the user about which areas appear to be adequately or inadequately defined.
- the system selects a migration method that is more sophisticated than the one used to generate the original data migration.
- the original migration could have employed one-way WEM, but the system may be capable of implementing full-wave RTM.
- the user may select the desired migration method.
- the system determines source-receiver visibility V(s,r) using the selected migration method to simulate shots in the tentative structure identified in the original depth-migrated data. As previously mentioned, the visibility is determined by calculating reflection wavefields from the target events for each of
- the zero-offset visibility calculations can be simplified using the "exploding reflector" assumption, in which the target events are treated as a distributed multiple source positions, and then measuring the signal contributions from these reflection wavefields to the signals measured at each of multiple receiver positions.
- the system identifies those existing traces whose source- receiver positions have target event visibilities above a given threshold.
- the threshold can be preset, based on a peak visibility value, or selected to capture a predetermined fraction (e.g., 90%) of the area under the multidimensional visibility surface.
- the system applies the selected migration method to the identified high-visibility traces. Because the identified traces are expected to represent a small subset of the available data, the use of the more sophisticated migration method may be eminently feasible.
- the system determines whether the target events have been adequately imaged, and if so, the method jumps ahead to block 922. In some implementations, the system makes this determination by displaying the depth- migrated data to a user and soliciting user feedback. If the target event is still inadequately imaged, it is expected that additional data acquisition will be needed. Consequently, in block 916, the system identifies a survey region and other survey parameters based at least in part on the source- receiver visibility calculations. In some implementations, the range of desirable source and receiver positions can be determined by drawing a rectangle (for land surveys) or a parallelogram (for marine surveys where the receiver position varies with source position) that encloses the substantial bulk of the high-visibility value region.
- Fig. 10 shows an illustrative visibility analysis system in the form of a computer 50 having one or more input devices 54 and one or more output devices 56 through which it can interact with a user.
- Software illustrated as portable information storage disks 52 configures the computer's processor(s) to receive user commands and responsively retrieve data from network or internal storage,
- the software can typically distribute the burden of processing the data across multiple computers interconnected by a network.
- the foregoing description relies on a 2D seismic model for explanatory purposes.
- a 3D volume is being imaged, and that each of the source and receiver positions are specified in terms of at least two spatial coordinates.
- the source- receiver visibility map (see Fig. 7) is expected to have at least four spatial dimensions. Nevertheless, the underlying principles are the same.
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Claims
Priority Applications (8)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200980154752.0A CN102282481B (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis |
CA2747146A CA2747146C (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis |
MX2011007084A MX2011007084A (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis. |
AU2009337134A AU2009337134B2 (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis |
PCT/US2009/031400 WO2010082938A1 (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis |
US13/142,478 US9329288B2 (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis |
BRPI0924056-0A BRPI0924056A2 (en) | 2009-01-19 | 2009-01-19 | Seismic prospecting, seismic migration and visibility analysis system |
EP09838520.6A EP2376947B1 (en) | 2009-01-19 | 2009-01-19 | Prestack migration based on seismic visibility analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2009/031400 WO2010082938A1 (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2010082938A1 true WO2010082938A1 (en) | 2010-07-22 |
Family
ID=42340028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2009/031400 WO2010082938A1 (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis |
Country Status (8)
Country | Link |
---|---|
US (1) | US9329288B2 (en) |
EP (1) | EP2376947B1 (en) |
CN (1) | CN102282481B (en) |
AU (1) | AU2009337134B2 (en) |
BR (1) | BRPI0924056A2 (en) |
CA (1) | CA2747146C (en) |
MX (1) | MX2011007084A (en) |
WO (1) | WO2010082938A1 (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8116168B1 (en) | 2008-06-18 | 2012-02-14 | Halliburton Energy Services, Inc. | Hybrid one-way and full-way wave equation migration |
WO2012166737A2 (en) | 2011-05-27 | 2012-12-06 | Conocophillips Company | Two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
WO2012166733A1 (en) | 2011-05-27 | 2012-12-06 | Conocophillips Company | Two-way wave equation targeted data selection for improved imaging of prospects among complex geologic structures |
WO2013048586A2 (en) | 2011-09-28 | 2013-04-04 | Conocophillips Company | Reciprocal method two way wave equation targeted data selection for seismic acquisition of complex geologic structures |
WO2013048585A1 (en) * | 2011-09-28 | 2013-04-04 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for improved imaging of complex geologic structures |
US8750072B2 (en) | 2009-04-16 | 2014-06-10 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a fast target-oriented illumination calculation |
US8830788B2 (en) | 2011-02-24 | 2014-09-09 | Landmark Graphics Corporation | Sensitivity kernal-based migration velocity analysis in 3D anisotropic media |
CN104570055A (en) * | 2015-02-04 | 2015-04-29 | 中国海洋石油总公司 | Seismic source control system of air gun |
US9329288B2 (en) | 2009-01-19 | 2016-05-03 | Landmark Graphics Corporation | Data acquisition and prestack migration based on seismic visibility analysis |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2443436B (en) * | 2006-11-04 | 2010-07-21 | Westerngeco Seismic Holdings | Method of determing properties of the earth using turning wave data |
US10061046B2 (en) * | 2012-02-06 | 2018-08-28 | Ion Geophysical Corporation | Integrated passive and active seismic surveying using multiple arrays |
GB2510873A (en) | 2013-02-15 | 2014-08-20 | Total Sa | Method of modelling a subsurface volume |
GB2510872A (en) * | 2013-02-15 | 2014-08-20 | Total Sa | Method of modelling a subsurface volume |
US20150226868A1 (en) * | 2014-02-10 | 2015-08-13 | Microseismic, Inc. | Method for microseismic event moment magnitude estimation |
EP3371749A1 (en) | 2015-11-06 | 2018-09-12 | Google LLC | Regularizing machine learning models |
FR3059700A3 (en) * | 2016-12-02 | 2018-06-08 | Al Khatib Habib Agissant Au Nom Et Pour Le Compte De La Soc Spotlight En Cours De Formation | METHOD OF IMPROVING THE ULTRALIGHT ACTIVE SEISMIC DETECTION SYSTEMS |
US11137509B2 (en) * | 2018-02-02 | 2021-10-05 | Fairfield Industries, Inc. | Seismic imaging with a temporal decomposition imaging condition |
CN112782754A (en) * | 2019-11-06 | 2021-05-11 | 中国石油天然气集团有限公司 | Multisource synchronous acquisition data simulation method and device |
US11474267B2 (en) | 2020-06-11 | 2022-10-18 | China Petroleum & Chemical Corporation | Computer-implemented method and system employing compress-sensing model for migrating seismic-over-land cross-spreads |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5784334A (en) * | 1996-03-13 | 1998-07-21 | Atlantic Richfield Company | Method and system for detecting hydrocarbon reservoirs using amplitude versus offset analysis of seismic signals |
US6311131B1 (en) * | 1996-05-08 | 2001-10-30 | Lloyd Peardon | Method of and apparatus for determining the quality of seismic data |
US20020033832A1 (en) * | 2000-09-18 | 2002-03-21 | Rafail Glatman | Method for computer modeling of visual images and wave propagation |
US20050207278A1 (en) * | 2002-11-12 | 2005-09-22 | Landmark Graphics Corporation | Seismic analysis using post-imaging seismic anisotropy corrections |
US20050270537A1 (en) * | 2004-06-08 | 2005-12-08 | Mian Zahid F | Image-based visibility measurement |
Family Cites Families (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5014230A (en) | 1986-06-05 | 1991-05-07 | Xox Corporation | Solid-modeling system using topology directed subdivision for determination of surface intersections |
US4890242A (en) | 1986-06-05 | 1989-12-26 | Xox Corporation | Solid-modeling system using topology directed subdivision for determination of surface intersections |
US4953142A (en) * | 1989-01-06 | 1990-08-28 | Marathon Oil Company | Model-based depth processing of seismic data |
US5138584A (en) | 1989-09-06 | 1992-08-11 | Chevron Research & Technology Company | Migration of seismic turning waves |
US5274605A (en) | 1992-06-26 | 1993-12-28 | Chevron Research And Technology Company | Depth migration method using Gaussian beams |
US5490120A (en) | 1993-04-30 | 1996-02-06 | Union Oil Company Of California | Method of imaging overturned waves |
US5530679A (en) | 1993-05-10 | 1996-06-25 | Western Atlas International, Inc. | Method for migrating seismic data |
US5544126A (en) | 1993-11-04 | 1996-08-06 | Berryhill; John R. | System and method of seismic shot-record migration |
US6002642A (en) | 1994-10-19 | 1999-12-14 | Exxon Production Research Company | Seismic migration using offset checkshot data |
US6446007B1 (en) | 1998-11-05 | 2002-09-03 | Exxonmobil Upstream Research Company | Method for controlled-amplitude prestack time migration of seismic data |
US6021094A (en) | 1998-12-03 | 2000-02-01 | Sandia Corporation | Method of migrating seismic records |
US6687618B2 (en) | 2000-08-07 | 2004-02-03 | 3D Geo Development, Inc. | Typing picks to horizons in migration velocity analysis |
WO2002023222A1 (en) | 2000-09-15 | 2002-03-21 | Nutec Sciences, Inc. | Illumination weighted imaging condition for migrated seismic data |
US6611761B2 (en) | 2000-12-19 | 2003-08-26 | Schlumberger Technology Corporation | Sonic well logging for radial profiling |
US6864890B2 (en) | 2002-08-27 | 2005-03-08 | Comoco Phillips Company | Method of building and updating an anisotropic velocity model for depth imaging of seismic data |
US6778909B1 (en) | 2002-10-05 | 2004-08-17 | 3Dgeo Development, Inc. | Seismic data processing systems and methods |
US6819628B2 (en) | 2003-04-07 | 2004-11-16 | Paradigm Geophysical (Luxembourg) S.A.R.L. | Wave migration by a krylov space expansion of the square root exponent operator, for use in seismic imaging |
FR2858063B1 (en) * | 2003-07-21 | 2005-09-16 | Geophysique Cie Gle | METHOD FOR ESTIMATING THE ILLUMINATION COVERAGE RATE IN THE MIGRE DOMAIN |
US7065004B2 (en) | 2004-04-22 | 2006-06-20 | Pgs Americas, Inc. | Horizon-based residual depth migration velocity analysis |
US7376517B2 (en) | 2005-05-13 | 2008-05-20 | Chevron U.S.A. Inc. | Method for estimation of interval seismic quality factor |
US20060270537A1 (en) * | 2005-05-27 | 2006-11-30 | Karol Fisher | Exercising system with trellis |
US20070162249A1 (en) | 2006-01-06 | 2007-07-12 | Min Lou | Traveltime calculation in three dimensional transversely isotropic (3D TTI) media by the fast marching method |
US7196969B1 (en) | 2006-02-09 | 2007-03-27 | Pgs Geophysical As | Three-dimensional two-way acoustic wave equation pre-stack imaging systems and methods |
US7952960B2 (en) | 2006-10-03 | 2011-05-31 | Bp Corporation North America Inc. | Seismic imaging with natural Green's functions derived from VSP data |
US8120991B2 (en) | 2006-11-03 | 2012-02-21 | Paradigm Geophysical (Luxembourg) S.A.R.L. | System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems |
US7400553B1 (en) * | 2006-11-30 | 2008-07-15 | Shengwen Jin | One-return wave equation migration |
US20080137480A1 (en) | 2006-12-07 | 2008-06-12 | Woodside Energy Limited | Method of Building a Subsurface Velocity Model |
US8194498B2 (en) | 2008-01-18 | 2012-06-05 | Westerngeco L.L.C. | Using a wave propagator for transversely isotropic media |
US7663972B2 (en) | 2008-02-22 | 2010-02-16 | Pgs Geophysical As | Method for three dimensional seismic travel time tomography in transversely isotropic media |
US8116168B1 (en) | 2008-06-18 | 2012-02-14 | Halliburton Energy Services, Inc. | Hybrid one-way and full-way wave equation migration |
US8296069B2 (en) | 2008-10-06 | 2012-10-23 | Bp Corporation North America Inc. | Pseudo-analytical method for the solution of wave equations |
CN102282481B (en) | 2009-01-19 | 2014-11-12 | 兰德马克图形公司 | Data acquisition and prestack migration based on seismic visibility analysis |
WO2010120301A1 (en) | 2009-04-16 | 2010-10-21 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a fast target-oriented illumination calculation |
US8332156B2 (en) | 2009-07-10 | 2012-12-11 | Chevron U.S.A. Inc. | Method for propagating pseudo acoustic quasi-P waves in anisotropic media |
US8406081B2 (en) | 2009-09-25 | 2013-03-26 | Landmark Graphics Corporation | Seismic imaging systems and methods employing tomographic migration-velocity analysis using common angle image gathers |
AU2009354764B2 (en) | 2009-11-02 | 2014-03-20 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a 3D reverse time migration with tilted transverse isotropy |
US8830788B2 (en) | 2011-02-24 | 2014-09-09 | Landmark Graphics Corporation | Sensitivity kernal-based migration velocity analysis in 3D anisotropic media |
-
2009
- 2009-01-19 CN CN200980154752.0A patent/CN102282481B/en not_active Expired - Fee Related
- 2009-01-19 EP EP09838520.6A patent/EP2376947B1/en active Active
- 2009-01-19 AU AU2009337134A patent/AU2009337134B2/en not_active Ceased
- 2009-01-19 MX MX2011007084A patent/MX2011007084A/en active IP Right Grant
- 2009-01-19 US US13/142,478 patent/US9329288B2/en active Active
- 2009-01-19 BR BRPI0924056-0A patent/BRPI0924056A2/en not_active IP Right Cessation
- 2009-01-19 CA CA2747146A patent/CA2747146C/en active Active
- 2009-01-19 WO PCT/US2009/031400 patent/WO2010082938A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5784334A (en) * | 1996-03-13 | 1998-07-21 | Atlantic Richfield Company | Method and system for detecting hydrocarbon reservoirs using amplitude versus offset analysis of seismic signals |
US6311131B1 (en) * | 1996-05-08 | 2001-10-30 | Lloyd Peardon | Method of and apparatus for determining the quality of seismic data |
US20020033832A1 (en) * | 2000-09-18 | 2002-03-21 | Rafail Glatman | Method for computer modeling of visual images and wave propagation |
US20050207278A1 (en) * | 2002-11-12 | 2005-09-22 | Landmark Graphics Corporation | Seismic analysis using post-imaging seismic anisotropy corrections |
US20050270537A1 (en) * | 2004-06-08 | 2005-12-08 | Mian Zahid F | Image-based visibility measurement |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8116168B1 (en) | 2008-06-18 | 2012-02-14 | Halliburton Energy Services, Inc. | Hybrid one-way and full-way wave equation migration |
US9329288B2 (en) | 2009-01-19 | 2016-05-03 | Landmark Graphics Corporation | Data acquisition and prestack migration based on seismic visibility analysis |
US8750072B2 (en) | 2009-04-16 | 2014-06-10 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a fast target-oriented illumination calculation |
US8830788B2 (en) | 2011-02-24 | 2014-09-09 | Landmark Graphics Corporation | Sensitivity kernal-based migration velocity analysis in 3D anisotropic media |
WO2012166733A1 (en) | 2011-05-27 | 2012-12-06 | Conocophillips Company | Two-way wave equation targeted data selection for improved imaging of prospects among complex geologic structures |
EP2715404A2 (en) * | 2011-05-27 | 2014-04-09 | ConocoPhillips Company | Two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
EP2715403A1 (en) * | 2011-05-27 | 2014-04-09 | ConocoPhillips Company | Two-way wave equation targeted data selection for improved imaging of prospects among complex geologic structures |
WO2012166737A3 (en) * | 2011-05-27 | 2014-05-01 | Conocophillips Company | Two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
WO2012166737A2 (en) | 2011-05-27 | 2012-12-06 | Conocophillips Company | Two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
US9164184B2 (en) | 2011-05-27 | 2015-10-20 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
EP2715404A4 (en) * | 2011-05-27 | 2015-01-28 | Conocophillips Co | Two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
EP2715403A4 (en) * | 2011-05-27 | 2015-01-28 | Conocophillips Co | Two-way wave equation targeted data selection for improved imaging of prospects among complex geologic structures |
US9279896B2 (en) | 2011-05-27 | 2016-03-08 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for improved imaging of complex geologic structures |
WO2013048585A1 (en) * | 2011-09-28 | 2013-04-04 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for improved imaging of complex geologic structures |
EP2761330A4 (en) * | 2011-09-28 | 2015-05-06 | Conocophillips Co | Reciprocal method two-way wave equation targeted data selection for improved imaging of complex geologic structures |
EP2761328A4 (en) * | 2011-09-28 | 2015-05-06 | Conocophillips Co | Reciprocal method two way wave equation targeted data selection for seismic acquisition of complex geologic structures |
WO2013048586A2 (en) | 2011-09-28 | 2013-04-04 | Conocophillips Company | Reciprocal method two way wave equation targeted data selection for seismic acquisition of complex geologic structures |
CN104570055A (en) * | 2015-02-04 | 2015-04-29 | 中国海洋石油总公司 | Seismic source control system of air gun |
Also Published As
Publication number | Publication date |
---|---|
EP2376947A1 (en) | 2011-10-19 |
CN102282481B (en) | 2014-11-12 |
EP2376947A4 (en) | 2015-09-02 |
US20110273959A1 (en) | 2011-11-10 |
BRPI0924056A2 (en) | 2015-07-14 |
AU2009337134A1 (en) | 2011-07-07 |
CA2747146A1 (en) | 2010-07-22 |
US9329288B2 (en) | 2016-05-03 |
EP2376947B1 (en) | 2021-08-25 |
CA2747146C (en) | 2016-10-04 |
AU2009337134B2 (en) | 2013-11-21 |
CN102282481A (en) | 2011-12-14 |
MX2011007084A (en) | 2011-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2009337134B2 (en) | Data acquisition and prestack migration based on seismic visibility analysis | |
US11698472B2 (en) | Method and system for separating seismic sources in marine simultaneous shooting acquisition | |
US10989825B2 (en) | Method and system for determining source signatures after source ghost removal | |
US7286939B2 (en) | Method for estimating porosity and saturation in a subsurface reservoir | |
EP1563324B1 (en) | Seismic analysis using post-imaging seismic anisotropy corrections | |
US6128580A (en) | Converted-wave processing in many-layered anisotropic media | |
US7952960B2 (en) | Seismic imaging with natural Green's functions derived from VSP data | |
EA020278B1 (en) | Method for determining seismic data quality | |
AU2009212927A1 (en) | Removing non-physical wavefields from interferometric Green's functions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 200980154752.0 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 09838520 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009337134 Country of ref document: AU |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2747146 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13142478 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: MX/A/2011/007084 Country of ref document: MX |
|
ENP | Entry into the national phase |
Ref document number: 2009337134 Country of ref document: AU Date of ref document: 20090119 Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009838520 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01A Ref document number: PI0924056 Country of ref document: BR |
|
ENP | Entry into the national phase |
Ref document number: PI0924056 Country of ref document: BR Kind code of ref document: A2 Effective date: 20110711 |