• 专利附录
  • 专利说明
  • 权利要求
  • 类似技术
  • 同族专利
  • 引用文献
  • 法律状态
  • 优先权
    • 专利摘要:
    Methods of Exploration of Hydrocarbons in the Interior of a Predetermined Volume of the Earth Containing Structural and Stratigraphic Characteristics Conducive to the Generation, Migration, Accumulation or Presence of Said Hydrocarbons. In accordance with a preferred aspect of the present invention, there is provided herein a system and method for extending the zero-shift or stacked wave equation illumination analysis into an angle-gathering domain, wherein it becomes a suitable tool for evaluating the effects of complex overloads on the AVA response. A preferred method of doing this involves first creating an angle meeting that has a perfect AVA response (that is, a constant amplitude as a function of angle). This gathering is then ideally used as a reflectivity map that is fed into a migration process that creates modeled data that by construction carry with them a completely flat reflectivity signature. Migration of this dataset then results in a pool where any variation in amplitude is more likely to be a measure of lighting effects alone. The resulting VLE signature on the meeting can then be used to assess the validity of (...)'s response.
    公开号:BR112012006931B1
    申请号:R112012006931-4
    申请日:2010-10-04
    公开日:2022-02-15
    发明作者:Uwe Albertin；Ole Joran Askim；Mariana Gherasim
    申请人:Bp Corporation North America Inc；
    IPC主号:
    专利说明:

    DESCRIPTIVE REPORT Related Case
    [001] This Application claims the benefit of US Provisional Patent Application, Serial No. 61/248,222, filed October 2, 2009, and incorporates said provisional application for reference in its exposition as if fully set forth at this point. Technical Field
    [002] This invention relates to the general subject of seismic exploration and, in particular, to methods for estimating seismic and other signals that are representative of the subsurface in areas of complex subsurface structure. Background of the Invention
    [003] A seismic survey represents an attempt to image or map the earth's subsurface by sending sound energy down into the ground and recording the "echoes" that return from the rock layers below. The source of the sound energy traveling downwards could come, for example, from explosions or seismic vibrators over land or air cannons in marine environments. During a seismic survey, the energy source is placed at various locations close to the earth's surface above a geological structure of interest. Each time the source is activated, it generates a seismic signal that travels down through the earth, is reflected, and, on its return, is recorded in a multitude of many locations on the surface. Multiple source/record combinations are then combined to create an approximately continuous subsurface profile that can span many miles. In a two-dimensional (2D) seismic survey, the recording locations are generally established along a single line, while in a three-dimensional (3D) survey the recording locations are distributed across the surface in a grid pattern. In the simplest terms, a 2D seismic line can be thought of as providing a cross-sectional (vertical slice) image of the earth's layers as they exist directly beneath the recording sites. A 3D prospect produces a “cube” or volume of data that is, at least conceptually, a 3D image of the subsurface that lies beneath the prospect area. In reality, however, both 2D and 3D surveys interrogate some volume of land lying beneath the area covered by the survey.
    [004] A seismic survey is made up of a very large number of individual seismic records or traces. In a typical 2D prospect, there will usually be several tens of thousands of traces, while in a 3D prospect the number of individual traces can extend to multiples of millions of traces. Chapter 1, pages 9 - 89, of 'Seismic Data Processing' by Ozdogan Yilmaz, Society of Exploration Geophysicists, 1987, contains general information relating to conventional 2D processing and this exposition is incorporated here for reference. General Background information pertaining to 3D data acquisition and processing can be found in Chapter 6, pages 384-427, by Yilmaz, the disclosure of which is also incorporated herein by reference.
    [005] A seismic trace is a digital record of acoustic energy that is reflected from inhomogeneities or discontinuities in the subsurface, a partial reflection that occurs each time there is a change in the elastic properties of subsurface materials. Digital samples are usually taken at 0.002 second (2 millisecond or “ms”) intervals, although 4 millisecond and 1 millisecond sampling intervals are also common. Each discrete sample in a conventional digital seismic trace is associated with a discrete sampling of the time-reflected acoustic wave field. Many variations of the conventional source-receiver arrangement are used in practice, for example VSP surveys (vertical seismic profiles), ocean floor surveys, etc. Also, the surface location of each trace in a seismic survey is carefully tracked and is usually made a part of the trace itself (as part of the trace header information). This allows the seismic information contained within the traces to be further correlated to specific surface and subsurface locations, thus providing a means for posting and contouring seismic data – and attributes extracted from them – on a map (i.e. “mapping”). .
    [006] The data in a 3D prospect is responsible for observation in a number of different ways. First, horizontal “constant time slices” can be taken extracted from a stacked or unstacked seismic volume by collecting all of the digital samples that reflect from a given subsurface location after correcting those samples for propagation effects. acoustics. This operation results in a horizontal 2D plane of seismic data. By animating a series of 2D planes it is possible for the interpreter to move through the volume, giving the impression that successive layers are being torn away so that the information that is situated underneath can be observed. Similarly, a vertical plane of seismic data can be taken at an arbitrary azimuth through the volume by collecting and displaying the seismic traces that lie along a particular line. This operation, in effect, extracts an individual 2D seismic line from within the 3D data volume. It should also be noted that a 3D dataset can be thought of as being made up of a 5D dataset that has been reduced in dimensionality by stacking it into a 3D image. Dimensions are typically time (or “z” depth), “x” (e.g. North-South), “y” (e.g. East-West), source-receiver offset in the x direction, and source-to-receiver offset in the x direction. -receiver in the y direction. While the examples given here may focus on 2D and 3D cases, extending the process to four or five dimensions is straightforward.
    [007] Seismic data that have been properly obtained and processed can provide a wealth of information for the explorer, one of the individuals within an oil company whose task is to locate potential drilling sites. For example, a seismic profile provides the explorer with a broad view of the subsurface structure of rock layers and often reveals important features associated with hydrocarbon trapping and storage, such as faults, folds, anticlines, nonconformities, and salt domes and reefs. subsurface, among many others. During computer processing of seismic data, estimates of subsurface rock velocities are routinely generated and near-surface inhomogeneities are detected and displayed. In some cases, seismic data can be used to directly estimate rock porosity, water saturation, and hydrocarbon content. Less obviously, seismic waveform attributes, such as phase, peak amplitude, peak-to-valley ratio, and a host of others, can often be empirically correlated with known hydrocarbon occurrences and this correlation applied to seismic data collected on new ones. exploration targets.
    [008] Many variations of the conventional source-receiver arrangement are used in practice, for example VSP (vertical seismic profile) prospecting, ocean floor prospecting etc.
    [009] Seismic attributes such as amplitude versus displacement (“AVO”) or amplitude versus angle of incidence (“AVA”) analysis can produce important information about the contents of subsurface rock formations. Although hydrocarbons cannot generally be observed directly in the subsurface using seismic, variations in reflectivity with angle of incidence have increasingly been used as an attribute or indicator of the presence of gas in the subsurface. See, for example, Castagna and Swan, "Principles of AVO Crossplotting," The Leading Edge, April 1997, the disclosure of which is incorporated herein by reference. However, deeper targets pose a number of problems for this technology, not the least of which is related to the distortion that can be introduced by the subsurface structure and/or the processing methods that are used to image this structure.
    [0010] One of the main aspects in the continuous development of these areas of complex geology is good planning, which often must be done in geological settings where obtaining good seismic images can be a challenge. Since AVA is often used to determine the potential for the well site, any irregularities in the AVA response due to uneven acoustic lighting resulting from complex sterile deposits introduce substantial risk to the AVA analysis and could very adversely affect well placement.
    [0011] Imaging the subsurface in regions of complex structure is problematic because the seismic wave field can be distorted significantly when it passes through such complexity. Of particular interest for the purpose of the present exhibit is imaging in the presence of subsurface salt. Seismic surveys that include subsurface salt features (eg, salt domes) can yield data that is marked by uneven illumination of reflectors below the salt (or other structure). This, in turn, can make VLE-type analyzes difficult to interpret and/or simply dubious. In the case of a salt dome, the distortion in the wave field is caused by the large velocity contrast between the salt and the surrounding rock (i.e., the salt typically has a seismic velocity that is much greater than that of the surrounding sedimentary rocks). This velocity contrast results in large amounts of radius bending and rays that are normal to the target reflector will tend to critically pass at the sediment salt interface. Conventional seismic imaging methods do not adequately compensate for this uneven lighting, which can distort observed trace amplitudes and can make AVO/AVA analysis questionable.
    [0012] Thus, what is needed is a method of compensating seismic collections for lighting irregularities caused by the structure, the effects of obtaining available space, and wave propagation effects in complex structural areas while simultaneously preserving the reflectivity signature of AVA.
    [0013] Until now, as is well known in the seismic processing and seismic interpretation arts, there has been a need for a method of obtaining better estimates of the effect of AVA in areas with a complex geological subsurface structure. Therefore, it must now be recognized, as has been recognized by the present inventor, that there is, and has been for some time, a very real need for a method of processing seismic data that would address and solve the problems described above.
    [0014] Before proceeding with the description of the present invention, however, it should be noted and remembered that the description of the invention that follows, together with the accompanying drawings, is not to be understood as limiting the invention to the examples (or preferred embodiments) shown. and described. This is so because those skilled in the art to which the invention pertains will be able to devise other forms of this invention within the scope of the appended Claims. Summary of the Invention
    [0015] In accordance with a preferred aspect of the present invention, there is provided here a system and method for extending zero displacement or stacked wave equation illumination analysis into the angle grouping domain, where it becomes an effective tool for to evaluate the effects of complex sterile deposits on AVA response. A preferred method of accomplishing this involves first creating a grouping of angles (see, for example, USPN 4,646,239, the exposure of which is incorporated for reference) that has a perfect AVA response (that is, a constant amplitude as a function of angle).
    [0016] This cluster is then preferably used as a reflectivity map that is used in a process of defacing or modeling the output from which it is modeled since, by construction, it carries with itself a completely flat reflectivity signature. “Disfigurement” is a process by which a depth-migrated dataset is used to calculate an estimate of an original common displacement section from which it could have been obtained. Note that for the purposes of the present exposition, the terms “modeling” and “unit reflector modeling” should be understood to be a process by which a synthetic seismic dataset is produced from a model of the earth bearing a plane reflectivity signature. (that is, one in which there is no variation in reflection amplitude with varying angle of incidence).
    [0017] Note that this definition should be broadly interpreted to include the unit-amplitude reflector type discussed above, as well as models that are generated where density reflectors are used in two-way acoustic time modeling.
    [0018] Those of ordinary skill in the art will recognize that the migration operator adjunct can be defined as defacement. Of course, the deface operation can be used to convert a depth-migrated section into a time-domain dataset that approximates the original obtained data.
    [0019] Thus, re-emigration of a demigrated or modeled dataset produced in accordance with the present invention then produces a cluster over which any amplitude variation is most likely to be a measure of lighting effects alone. The resulting AVA signature over the cluster can then be used to help distinguish whether amplitude variation over the angle clusters is due to lighting effects during propagation, or actual rock properties. This proposal also preferably produces an AVA confidence analysis that can help the explorer determine when the AVA signatures are relatively free from lighting effects.
    [0020] The foregoing has outlined in broad terms the most important features of the invention set forth herein so that the detailed description which follows can be more clearly understood, and so that the present inventor's contribution to the art can be better taken into account. The present invention is not to be limited, in its application, to the details of construction and arrangements of components set out in the following description or illustrated in the drawings. Rather, the invention is capable of other embodiments and of being practiced and performed in various other ways not specifically enumerated here. Finally, it is to be understood that the phraseology and terminology employed herein are for the purpose of description and should not be construed as limiting unless the specification specifically limits the invention. Brief Description of Drawings
    [0021] Other objects and advantages of the invention will become evident when reading the following detailed description and referring to the drawings, in which:
    [0022] FIG. 1 illustrates the general environment of the present invention;
    [0023] FIG. 2 illustrates a seismic processing sequence suitable for use with the present invention;
    [0024] FIG. 3 contains a schematic illustration of how the present invention could be used in an exploration configuration;
    [0025] FIG. 4 is a preferred operating logic suitable for use with the present invention;
    [0026] FIG. 5 is a continuation of the preferred operating logic from Figure 4 with additional details regarding the backprojection link;
    [0027] FIGS. 6A-6C contain a schematic representation of the migration/disfiguring/shaping process of the present invention. Detailed Description
    [0028] While this invention is capable of embodiment in many different ways, some specific embodiments of the present invention are shown in the drawings and will be described hereinafter in detail. It should be understood, however, that the present disclosure is to be considered an exemplification of the principles of the invention and is not intended to limit the invention to the specific modalities or algorithms described herein. General Environment of the Invention
    [0029] Figure 1 illustrates the general environment in which the present invention would typically be used. A seismic survey is projected at 110 by the explorer to cover an area of economic interest. Field acquisition parameters (e.g. shot spacing, line spacing, bending, etc.) are typically selected in conjunction with this step, although it is common to modify optimal design parameters slightly (or substantially) in the field to accommodate the realities of conducting prospecting.
    [0030] Seismic data are collected at field 120 over a subsurface target of potential economic importance and are typically sent onward to a processing center 150 where they will be processed for use in exploration. In some cases, some initial data processing can be done in the field and this is becoming more common and possible given the computing power that is available to field teams.
    [0031] In the processing center, a variety of preparatory processes 130 are applied to the seismic traces to make them ready for use by the methods set forth hereinafter. The processed traces would then be made available for use by the present invention and could be stored, by way of example only, on hard disk, magnetic tape, magneto-optical disk, DVD disk, or other mass storage devices.
    [0032] The methods set forth here would best be implemented in the form of a computer program 140 that has been loaded into a general-purpose programmable computer 150 where it is accessible by a seismic interpreter or processor. Note that a general purpose computer 150 would typically include, in addition to central processing units and workstations, computers that provide parallel and massively parallel computations, where the computational load is distributed among two or more processors. As is also illustrated in Figure 1, in the illustrated arrangement some type of digitized zone-of-interest model 160 can be specified by the user and provided as feed to the processing computer program. In the case of a 3D seismic section, the zone-of-interest model 160 would typically include specifics regarding the lateral extent and thickness (which could be variable and could be measured in time, depth, frequency, etc.) of a subsurface target. The exact means by which such zones are created, captured, digitized, stored and subsequently read during program execution are not important to the present invention and those skilled in the art will recognize that this could be done in any number of ways.
    [0033] A program 140 embodying the present invention could be transported to the computer that must execute it by means of, for example, a floppy disk, a magnetic disk, a magnetic tape, a magneto-optical disk, an optical disk, a CD-ROM, a DVD disk, a RAM card, flash RAM, a RAM card, a PROM chip, or loaded over a network. In a typically seismic processing environment, the methods of the present invention would be made part of a larger package of software modules that is designed to perform many of the processing steps listed in Figure 2. After processing by the present methods, the resulting traces would then typically be sorted into clusters, stacked, and displayed either on a high resolution color computer monitor 170 or in hardcopy form as a section or a printed seismic map 180. The seismic interpreter would then use the displayed images to assist in identifying of subsurface features conclusive for the generation, migration, or accumulation of hydrocarbons.
    [0034] As previously indicated, the present invention will preferably be made separate and incorporated into a conventional seismic processing sequence of the type generally described in Figure 2. Those of ordinary skill in the art will recognize that the processing steps illustrated in Figure 2 are only broadly representative of the types of processes that could be applied to such data, and the choice and order of processing steps and the particular algorithms involved can vary markedly depending on the individual seismic processor, signal source (dynamite, vibrator, etc.) , the location where the data is prospected (land, sea, etc.), the company that processes the data, etc.
    [0035] As a first step, and as is generally illustrated in Figure 2, a 2D or 3D seismic survey is conducted over a particular volume of the earth's subsurface (step 210). The data collected in the field consists of unstacked (ie, not summed) seismic traces that contain digital information representative of the volume of land that is situated beneath the prospect. Methods by which such data are obtained and processed in a form suitable for use by seismic processors and interpreters are well known to those of ordinary skill in the art.
    [0036] The purpose of a seismic survey is to obtain a collection of spatially related seismic traces over a subsurface target of some potential economic importance. Data that are suitable for analysis by the methods set forth herein could consist of, for illustrative purposes only, an unstacked 2-D seismic line, an unstacked 2-D seismic line extracted from a 3D seismic survey, or, preferably , a non-stacked 3D portion of a 3D seismic survey, or a 4D or 5D survey, etc. The invention set forth herein is most effective when applied to a group of seismic features that have an underlying spatial relationship to some subsurface geological feature. Again, for illustrative purposes only, the discussion that follows will be based in terms of features contained in a 3-D prospect (stacked or unstacked, as the discussion warrants), although any group of installed seismic features, spatially related, could conceivably be used.
    [0037] After the seismic data are obtained (step 210), they are typically taken to a processing center where some initial or preparatory processing steps are applied to them. As illustrated in Figure 2, a common pre-step 215 is designed to edit incoming seismic data in preparation for subsequent processing (e.g. demux, gain recovery, small wave shaping, faulty trace removal, etc.). This could be followed by specifying prospect geometry (step 220) and storing shot/receiver number and one and one surface location as part of each seismic trace header. Once the geometry has been specified, it is customary to perform a velocity analysis consisting of NMO analysis if the processing is time processing, or NMO followed by RMS for interval velocity conversion, followed by depth migration and tomography to obtain a initial velocity model for depth migration.
    [0038] After the initial pre-stack processing is completed, it is customary to condition the seismic signal on the unstacked seismic traces before creating stacked (or summed) data volumes (step 230). In Figure 2, step 230 contains a typical "Signal Processing / Conditioning / Imaging" processing sequence, but those skilled in the art will recognize that many alternative processes could be used in place of those listed in the figure. In any case, the ultimate goal from the explorer's point of view is the production of a seismic volume or, in the case of 2D data, a seismic line for use in exploration for hydrocarbons within the subsurface of the earth.
    [0039] In some preferred arrangements, the present invention could be better used in connection with step 230. That being said, those of ordinary skill in the art will recognize that there are many other points in a typical processing sequence that could be improved through of the use of the present invention.
    [0040] As is further suggested in Figure 2, any digital sample within a seismic volume is uniquely identified by a vector (X, Y, OFFSETX,OFFSETY,TIME), with the X and Y coordinates representing some position on the surface from the earth, the OFFSETX and OFFSEDADE coordinates specifying the distance between source and receiver, and the time coordinate measuring an arrival time recorded within the seismic trace (step 240). For specificity purposes, it will be assumed that the X direction corresponds to the “in line” direction, and the Y measurement corresponds to the “cross line” direction, as the terms “in line” and “cross line” are generally understood in the art. time is a preferred and more common vertical axis unit, those skilled in the art understand that other units are certainly possible which could include, for example, depth or frequency. Additionally, it is well known to those skilled in the art that it is possible to convert seismic traces from one axis unit (e.g. time) to another (e.g. depth) using standard mathematical conversion techniques. In addition, depending on whether the volume is imaged or unimaged, a sample in the volume can be determined by surface displacement (i.e. OFFSETX and OFFSETY) whether the volume is unimaged or imaged in the form of offset clusters, or alternatively by reflection aperture angle and azimuth, if the volume is imaged in angle grouping mode.
    [0041] After an image volume is stacked, the explorer can make an initial interpretation 250 of the resulting stacked volume, where he or she locates and identifies major reflectors and faults as they occur in the dataset. This could be followed by further data enhancement 260 of the stacked or unstacked seismic data and/or attribute generation (step 270) therefrom. In many cases, the explorer will revisit their original interpretation in light of additional information gained from improved data and attribute generation steps (step 280). As a final step, the explorer will typically use information gained from the seismic data in conjunction with other types of data (magnetic prospects, gravity prospects, LANDSAT data, regional geological studies, well logging, well cores, etc.) to locate subsurface structural or stratigraphic features conclusive for the generation, accumulation, or migration of hydrocarbons (ie, prospect generation 290). Preferred Modalities
    [0042] According to a first preferred aspect of the present invention, there is provided a system and method for extending zero displacement or stacked wave equation illumination analysis to the angle clustering domain, which becomes an appropriate tool for evaluating the effects of complex sterile deposits on the AVA response. A preferred method of accomplishing this involves first creating an angle cluster that has a perfect AVA response (i.e., a response where the reflected seismic events have a constant amplitude as a function of the angle of incidence with respect to the subsurface layers). This “perfect” cluster will then be used as a reflectivity map in a defacing or modeling process that creates modeled data, which, by construction, carries with it a completely flatAVA response. Re-migrating this dataset will then result in a cluster in which any amplitude variation is more likely to be a measure of lighting effects alone. The resulting VLE signature over the cluster can then be used to determine the validity of the VLE response on modeled or actual data, resulting in a useful VLE risk analysis.
    [0043] By way of general background, the amount of seismic energy that is reflected from a subsurface reflector at a subcritical angle varies depending (at least in part) on its angle of incidence with respect to the reflector. Furthermore, the magnitude of this effect is much more pronounced at the interface between a formation that contains gas and one that does not contain gas. This effect has made it possible to use AVA techniques to estimate subsurface elastic parameters from seismic data. Thus, it is customary to include this effect in synthetic seismic traces that are produced by conventional modeling programs. However, it should be noted that the present method specifically excludes such computations in forming its modeled seismic traces.
    [0044] The preferred embodiment of the present invention can be generally understood as follows. Assume that the observed seismic data can be represented by the familiar conceptual equation Seismic data = forward propagation * reflectivity, or D = FR. Conceptually, the F operator represents all of the effects of the current propagation through the earth. In the common practice of imaging seismic data, this operator cannot be found directly, and is instead approximated with a simpler modeling operator M, so that D ~ MR. It is therefore common to make a second approximation to obtain a seismic image. Because the M operator cannot be easily inverted, the seismic image is often obtained by applying an adjunct M* of the M operator rather than its inverse. Those of ordinary skill in the art will recognize that the "adjoint" of a square matrix is defined to be its conjugate transpose. In general, the migration process can be thought of as the process of applying the adjunct operator M* to the seismic data.
    [0045] Given the precedent, the migration process to obtain a subsurface reflectivity image can be written as: I = M*D. In order to improve this approximation, one proposal would be to treat the problem as a least squares problem, instead of using the M* operator to migrate the data. Using this proposal, an improved migration will be obtained: I = (M*M)-1M*D
    [0046] So what is needed is a means to commute the inverse of M*M. The inverse of M*M contains combined information about illumination at all slopes and aperture angles as well as migration resolution and amplitude fidelity. Additionally, if increased with the acquisition geometry, it contains information related to the available acquisition space.
    [0047] The M*M quantity contains information related to three aspects of recorded seismic data, ie illumination, imaging amplitude fidelity and available acquisition space. However, for simplicity, the term “illumination information” will be used here to refer to all three types of information. Note that even if the operator (M*M) were known in its entirety, the computation of its inverse would be costly in terms of computational resources and thus unfavorable in most cases.
    [0048] M will be referred to here as a disfigurement operator, since M* represents a migration operator. A brute force proposition to find (M*M) is a difficult proposition since the arrays involved are very large (e.g. nxm*nym*nh*nt*nx*ny*nz elements where nxm is the number of center point locations in the x direction, nym is the number of trace locations in the y direction, nh is the number of shifts in a narrow azimuth geometry, nt is the number of time samples, and (nx,ny,nz) are the dimensions of the reflectivity model).
    [0049] As a consequence, the following techniques will preferably be employed to implement the present method.
    [0050] An approximation to (M*M) can be obtained by applying this operator to appropriate subsets of unit amplitude of the reflectivity model. The space information required for lighting obtained in this way depends on the nature of the dataset used. The final product obtained from the demigrated/remigrated dataset is an aperture-angle cluster, obtained either directly from the migration process, or via inclined stacking in the case of a migration process that provides subsurface shifted clusters. An alternative product is a surface displacement grouping for in the case of a migration that provides surface displacement groupings. Those of ordinary skill in the art will recognize that a “tilted stack” (Radon transform, tau-p transform, etc.) is a seismic plane wave decomposition method. It can be calculated by applying a series of linear outputs to an unstacked seismic array and summing each output over displacement. Of course, there are computationally efficient means of calculating the slanted stack that displacement/brute force sum and those of common knowledge in the art will be familiar with the same.
    [0051] By way of explanation, suppose that a single unit amplitude sample representing a point diffractor is placed at some location (x,y,z,h) in the reflectivity model as a function of the three spatial dimensions (x,y ,z) and a subsurface displacement (h), and is then demigrated and re-emigrated (ie the (M*M) operator is applied). Through the use of this proposal, illumination information is obtained for all slopes and opening angles, since the inclined stack of a spot diffractor in all dimensions (x,y,z,h) generates the slope and opening angle of the point diffractor.
    [0052] In the preferred embodiment, an AVA confidence mapping method is taught. According to this embodiment, illumination information as a function of aperture angle is obtained by extending the previously discussed point diffractor to form a surface embedded in the reflectivity volume (xyz). Thus, if, for example, the spot diffractor is replaced by a flat surface of a datum in line and slope in transverse line, but remains a spot diffractor positioned at zero displacement in the subsurface displacement direction, the effect of extending the object to a plane in amounts of (x,y,z) for selecting lighting information for a single sloped stack component, or a single slope, in reflectivity space (x,y,z). Since the input dataset is still a point diffractor in the subsurface displacement direction, however, illumination for all aperture angles is retained. If this flat surface is now deformed to follow geological structure, it follows that the derived illumination information will be aperture angle illumination information specific to the current geological structure. In another preferred embodiment, reflectors of appropriate amplitude density are placed on the subsurface in a density model for use in two-way acoustic or elastic modeling. Forward modeling is then performed, and the resulting data is then migrated to form angle clusters, either directly or via inclined stacking, as indicated above.
    [0053] As is generally suggested in Figures 6A, 6B and 6C, the present invention preferably begins by creating a subsurface reflectivity or density/velocity model. As described in the previous section, a reflectivity or density surface that follows the interpreted structure will preferably be created in (x,y,z) space. The model will then be demigrated, or forward modeling will be applied to the model for generated modeled data ( Figure 6B) according to methods well known to those of ordinary skill in the art and then re-emigrated into the subsurface (Figure 6C) to form an angle cluster. This dataset results in lighting information as a function of aperture angle. The information used in the deface or shaping process corresponds to a “perfect” angle cluster, for example a cluster with no amplitude mark as a function of aperture angle (ie angle of incidence). After reviewing, the resulting amplitude mark on the angle cluster will thus be a function of variable illumination due to propagation through the complex structure, the amplitude handling of the migration algorithm itself and the available space for obtaining.
    [0054] The preferred modality described above is suitable for those migration processes that naturally produce subsurface displacement clusters, such as wave equation migration, or migrations that produce angle clusters directly to other imaging algorithms, such as Kirchhoff migration, the preferred output is a cluster in which each stripe of surface vector displacement is imaged independently, resulting in a cluster whose traces represent independent images from each vector displacement, or a "shift" as opposed to a “angle” grouping. For this type of imaging process, the input reflectivity will preferably be chosen to be the geological surface of the unit amplitude in the reflectivity space (xyz). The preferred embodiment of the invention is to then demigrate and remigrate this input reflectivity for each vector shift independently. Amplitude variation across the resulting cluster is then a direct indication of illumination variation. Note that in the case of forward data modeled from density reflectors, the migrated data output may take the form of a surface displacement cluster. The surface displacement cluster containing the lighting information can then be converted to an aperture angle cluster using standard surface displacement for subsurface angle techniques.
    [0055] These concepts can be extended to an AVA analysis of the lighting mark as follows: - Select a subsurface event. - Create a perfect synthetic grouping. For migration processes, for which the natural clustering output is subsurface displacement, the perfect cluster is a unit-amplitude geological surface embedded in the reflectivity space (x,y,z,h), with a peak positioned at the displacement of zero subsurface (h). The sloped stack of this cluster contains no amplitude variation as a function of the angle of incidence. In the case of forward modeling, a reflection experiment containing no variation in amplitude, as a function of aperture angle, can be created by inserting horizons of some appropriate density amplitude. Those of ordinary skill in the art will understand how to create such models. For migration processes, for which the natural output cluster is vector surface displacement, the perfect cluster is the unit-amplitude geological surface embedded in the reflectivity space (x,y,z), duplicated for all vector displacements of surface. - Demigrate the perfect cluster or perform forward modeling, and then re-migrate the resulting data back to the reflectivity volume. This perfect grouping now has all the effects of wave propagation, uneven lighting, acquisition geometry, etc., contained within it. For purposes of the present disclosure, such a "perfect" synthetic seismic dataset will be referred to herein as a calibration seismic dataset. - Perform AVO/AVA amplitude analysis on this new calibration cluster to obtain the lighting AVO/AVA amplitude mark, and use the results to “renormalize” amplitudes in the current seismic cluster.
    [0056] According to the preferred embodiment, a one-way wave equation disfigurement to model the data or, alternatively, a two-way modeling of an appropriate subsurface model will be used. Consequently, it is useful to contrast this proposal with other modeling options that are conventionally available. Table 1 contains such a comparison. Table 1: Lighting modeling options.

    [0057] The first column of the matrix indicates in a general way (from “++” / relatively fast to “-” / relatively slow) the relative computational speed of the associated algorithm. The second column generally indicates how accurate the associated algorithm is in the presence of an abnormally subsurface velocity, such as a salt body, with “-” indicating “relatively imprecise”. The “Multiple” column indicates whether the associated algorithm can accommodate multiples. The next column (ie “Model Complexity”) indicates how complex the input subsurface velocity model can realistically be.
    [0058] Finally, the last column indicates what control (if any) the user of the associated algorithm has over the angle-versus-amplitude signature that is inherent in a particular modeling method. In most two-way time modeling methods, the AVA signature is inherent in the process, and is controlled by the model parameters. In these methods, reflections are created as a direct result of applying differentiation operators in space and time. This is not the case for methods that are based on Born or Kirchhoff scattering. Here, the user has control over the AVA signature inherent in the modeling method, which should be described in more detail below. An important aspect of the preferred modeling proposal is that it is possible to avoid mixing lighting AVA effect with AVA characteristics inherent in the reflector due to the physical properties of rock or hydrocarbons. That being said, any one of the aforementioned techniques could prove to be useful in a particular situation and the previous table is intended to give a general indication of the advantages and disadvantages of each technique, and is not intended to exclude any particular proposal.
    [0059] Returning below for a detailed discussion of a preferred algorithm, let M be a modeling operator that approximates the seismic experiment ahead, for which the collected data is seismic data ^(s, r, a) as a function of source, receiver and frequency.
    [0060] Migration is typically defined as the Modeling-ahead adjunct M* in the Born or Kirchhoff approximation. These types of migration methods create a reflectivity map a(x, θ) which is a function of space and aperture angle, θ. The first order least squares correction for imaging is then given by
    For the purposes of the present exposition, an approximation for the M*M operator will be used, which provides lighting information for AVA. In what follows, the mathematics will be understood to be illustrative only, and details of derivations which are well known to those skilled in the art will be omitted for purposes of clarity. In the Born approximation, the equation describing forward scattering can be written schematically as follows:
    where s' and r' are chosen to be close to the reflection point x, a is the angle-dependent reflectivity, and S is an operator that converts the angle-dependent reflectivity to a matrix operator that scales the contribution to the data. modeled on the points (s', r') in the reflection process. The adjunct to the modeling equation is the migration equation
    The S* operator isolates energy as a function of aperture angle, and for one-way wave equation migration, this is often implemented as an inclined stack operator that converts subsurface displacement to aperture angle. For those migrations that generate angle grouping information directly, the S operator can be omitted.
    [0061] In a typical inversion process, one would normally try to solve for the angle-dependent reflectivity α, for which the modeled data matches the actual data. In the present case, however, only illumination information is sought, which suggests that the reflectivity signature should be eliminated from the model calculation, and instead the reflectivity that is independent of angle in the Born modeling equation should be used. . Since S* is, in effect, a skewed stack in the preferred embodiment, it follows that S is an adjunct skewed stack. Realizing that the angle-adjoint sloped stack of an angle-independent quantity provides a delta function on the zero subsurface displacement, the following prescription for computing angle-dependent illumination can be obtained.
    [0062] In summary and according to another preferred embodiment, the present invention will preferably operate as follows: - A seismic survey will be conducted over a subsurface region of interest. - A depth model (structural, stratigraphic etc.) of the subsurface will preferably be constructed, which includes a better reference configuration of the subsurface reflectors and their respective velocities (and densities, if available). Velocities, densities, etc., which would correspond to the presence of hydrocarbons in the subsurface will be excluded from this model, one goal being to determine the AVA/AVO response in the absence of such hydrocarbons. At this stage, the model need not be overly detailed, but rather will preferably at least reflect the main features of the subsurface. - A unit-amplitude reflector model of the subsurface will preferably be constructed by converting the structural model information into a unit-amplitude reflectivity model that follows the interpreted horizons in depth. Alternatively, density reflectors of an appropriate amplitude will be created in a density model to be used in two-way acoustic or elastic front modeling. - The unit amplitude reflector model will then preferably be demigrated using a background velocity model and the acquisition geometry. Alternatively, forward modeling using the density reflectors will be performed using two-way acoustic or elastic modeling. This will produce an estimate of the unmigrated seismic data that could have given rise to the artificially constructed reflectors. - Then the data from the previous step will preferably be migrated (ie re-migrated). In conjunction with this step, downward continued displacement (DCO) clusters (using, for example, CAWE (common azimuth wave equation) or trigger-record algorithms for NATS (narrow azimuth towed floating seismographic cable) and WATS (Wide Azimuth Towed Floating Seismographic Cable) will preferably be calculated and saved. Note that preferably DCO clusters will have unit amplitude at zero displacement and zero amplitude at non-zero displacement. This corresponds to clusters of plane angles with unit amplitude at all angles. Preferably, the output will take the form of non-stacked synthetic seismic clusters representing seismic characteristics that would be expected if the subsurface did not contain hydrocarbons. - Due to the effects of lighting, migration, and obtaining, however, in practice, the Output DCO clusters will not be perfect peaks (filtered), and consequently after stacking slopes In order to obtain groupings of angles, the data no longer has a constant amplitude as a function of the angle. However, these data now represent the best seismic data that can be obtained, assuming, of course, that the velocity model / subsurface model provided is accurate and the acquisition configuration has been correctly specified. - Next, the DCO data will preferably be converted into groupings of angles. - Alternatively, in the case of one-way or reverse-time migration, where right angle cluster output is possible, DCO cluster formation can be skipped, and angle cluster output directly. - On the other hand, in the case of forward modeling from density reflectors, a model containing horizons with appropriate density contrasts can be constructed. Forward modeling is then performed to create synthetic data, which is then migrated to form angle clusters containing lighting information.
    [0063] The angle clusters produced by the above process provide a baseline against which observed seismic data can be compared. More particularly, to the extent to which the seismic data produced according to the method discussed above differs from the observed seismic data, this difference could possibly be attributed to the presence of subsurface hydrocarbons, as the effects of this would preferably have been excluded from the model. . On the other hand, where the synthetic and current data sets are comparable, that is an indication that the amplitude variations between the model and current seismic data are not due to hydrocarbons but, in contrast, are more likely to be due to the effects of enlightenment, migration or acquisition.
    [0064] Of course, there are other theories that could account for observed discrepancies (eg noise, imperfections in the original subsurface model, absorption, mode conversion, etc.). However, a known systematic difference between current data and demigrated/migrated data will be the exclusion of hydrocarbon type reflectivities and velocities from the model data.
    [0065] Returning below to the discussion of how the present method could be implemented in practice, as indicated previously, a model that contains synthetic reflectors that follow the subsurface structure and have zero reflectivity at all displacements except the displacement of zero subsurface will preferably be created. Alternatively, a set of density reflectors will be created.
    [0066] Next, preferably, the synthetic reflectivity model will be demigrated using the Born dispersion equation, or alternatively, the density model will be modeled ahead to create synthetic data. In both cases, the data is then remigrated to obtain angle clusters. Since the modeled data is designed to have an AVA signature, which is angle-independent, any variation in amplitude over the resulting cluster must be due to uneven lighting alone. Those of ordinary skill in the art will recognize that for 3D data, "angle clustering" will be understood to include the traditional notion of the same, as well as angle clusters that are a function of subsurface azimuth, angle clusters obtained by stacking, substacking or capture over subsurface azimuth, etc.
    [0067] Based on the foregoing, the preferred algorithmic flow should be started with a filtered unit reflectivity model with slopes that generally mimic those currently found in the subsurface, or alternatively reflectors of appropriate amplitude density. Preferably, the reflectivity or density model will be one that generally honors the current subsurface structure, as indicated by the seismic data. Methods of creating such models from seismic data, well logs, gravity and magnetic data, etc., are well known to those of ordinary skill in the art.
    [0068] Through the application of the disfigurement/re-emigration or forward modeling/migration scheme taught here, the model/reference dataset will be subjected to the transmission effects of the (typically) complex sterile deposits represented by the model. Since the resulting dataset is multidimensional, it can be useful to display your information using a variety of display schemes including horizontal time slices, vertical cross-sections, angle groupings, various pseudo-color enhanced versions of them, etc.
    [0069] The present inventors have determined that abrupt changes in subsurface slope (e.g. rock units near the base of a salt dome), and particularly in cases where the base becomes steeper relative to the top, have the potential to create significant changes in lighting as a function of angle. This can potentially create blank spots in angle clusters where no reflected energy is present.
    [0070] More interesting for an AVA interpreter, however, is the spatial signature of the lighting mark as a function of aperture angle. Often, certain angle substacks are better than the full stack, and in the case where azimuth information is available, that is, certain azimuths may be preferred in different areas. In some cases, it is instructive to examine lighting displays as a function of aperture angle which can preferably be obtained by capturing amplitudes across events created by partially stacking angle groupings and coloring the events according to amplitude. Areas of strong lighting can then be seen in, for example, white, while areas that are poorly lit will be lower in amplitude and could be designated to be black in the display. From these types of illumination maps, the interpreter can often determine which partial stack has the best spatial interpretability, and which will be the most successful in well mapping. The display of variability in subsurface reflector illumination using images such as those described above can provide a basis for an objective discussion concerning which dataset could potentially be best for mapping the complex structure. In complex exploration work, the interpreter may be confronted with using several different sub-stack volumes in order to fully map the surroundings of a complex salt structure, and views such as the preceding ones could help select such volumes for incorporation into the active exploration dataset.
    [0071] After merging the structural interpretation from different volumes, the next preferred step is often the assessment of drilling targets based on seismic amplitudes. Those of ordinary skill in the art will recognize that caution should be exercised in using fused volume for amplitudes, however, because amplitudes in the fused volume will tend, at the target level, to a spatially heterogeneous angle contribution. In such a case, the interpreter would be well advised to exercise caution when analyzing bright spots. For example, in an area where very shifted angles have strong illumination compared to close displacements, amplitudes may be strong over long displacements compared to close ones, leading the interpreter to possibly classify the event as a Class III AVA anomaly, when, in fact, no such anomaly exists.
    [0072] By way of explanation, those of ordinary skill in the art will understand that AVA responses are broadly characterized in the type with a Class I AVA response being one where the top of the reservoir is represented by an increase in impedance (i.e., a seismic peak). These types of carbonated sands tend to show “obscurations” in the pile data. A Class III AVA is one that is characterized by a reflection where the top of the reservoir has a decrease in impedance (trench) compared to the rock above it. These types of reservoirs tend to produce the classic “bright spots” of the Gulf Coast. Finally, a Class II AVA response is one where the top of the reservoir is represented by a peak that decreases in amplitude, changing to a ditch at far displacements/angles. In a stacked seismic section, these types of gaseous sands can be nearly invisible due to amplitude cancellations caused by summing together the near and far displacements.
    [0073] In view of the foregoing, the present inventors believe that it is prudent to devise an AVA assessment workflow that quantitatively analyzes how lighting problems impact a detailed seismic reservoir response. In general, once the lighting dataset has been obtained, different methods can be devised that integrate the lighting response with standard reservoir modeling workflows, depending on the seismic analysis objective.
    [0074] For example, in a preferred embodiment, an initial reservoir elastic modeling can be performed to determine the AVA response in the absence of any lighting effects. Since the lighting volume contains only lighting effects, these effects can be transferred to the reservoir model response through the use of, for example, a matching filter. The filter will preferably be designed in such a way that its application to a perfect or calibration cluster (ie, one with no AVA signature) reproduces the lighting response. The filter will then be applied to the modeled response to the reservoir, thus transferring the lighting effect to the modeled response. Consequently, by having a “perfect”/calibration dataset, and one with detrimental effects from uneven lighting, matching filters can be designed, which will locally transfer lighting degradation to a systematically created dataset (or dataset). real) of choice. Note that the terms “match filter” or “match filter” must be understood to be one or more 1-D, 2D, 3-D etc filters, and should not be limited in interpretation to be a single 1-D filter. D that is applied to the entire dataset.
    [0075] From this type of information, it is now possible to indicate regions (eg polygons) within which an AVO inversion is expected to accurately reflect lithology or fluid content. Knowing this, it will become possible to better determine those areas where the reservoir filled with oil does not appear to be different than the reservoir filled with salt water. Furthermore, with these types of lighting filters in accordance with the present invention being readily available, different reservoir geometries can be tested and investigated to determine how non-exclusive the assessment and payout distribution is for a given saline geometry.
    [0076] A similar proposal is also applicable to synthetic data created directly from a reservoir model.
    [0077] For example, suppose a class III AVO behavior is observed in the original reservoir modeled data. Once lighting filters have been applied, it is possible that a Class III AVO could be degraded or eliminated. Obviously, if this is observed, the reliability of the observed seismic response could be questioned. In short, where there is good seismic illumination of a subsurface rock unit, confidence in the seismic data, and calculations made from it (e.g. AVO/AVA) will be highest compared to rock units with lower grades. of lighting. Naturally, where confidence in the calculations is highest, the explorer will be more inclined to place credit and confidence in the results.
    [0078] Returning next to Figure 3, this Figure illustrates some of the many ways the present invention (i.e., “Angle Grouping Lighting Assessment” or “AGILE” in this Figure) could be utilized in an exploration setting. . Generally speaking, and as indicated in this figure, data that are suitable for feeding into and/or use with the present invention could take one of two forms. In a first case, data such as the seismic survey acquisition geometry, a velocity model (simple or complex) of the subsurface, peaks/depths of important horizons etc. will preferably be used as input or feed for the present invention. The output could include, for example, lighting weights that could be further used in improved AVA modeling, AVO/AVA confidence surface/volume, improved image/cluster volumes from field data. All of this type of information could then be useful in, for example, creating a seismic characterization of a reservoir, looking at seismic attributes that are associated with, or predictable, of a reservoir, etc.
    [0079] Additionally, and as a second general source of data, the present invention could utilize information from well logs that was taken from wells that are close to the prospect. Additionally, reservoir models and geological models could be used to create binary models, specify 3D reservoir geometry, and/or create AVA models. AVA models created from well logs could then be compared with improved AVA models obtained via the present invention (e.g. where lighting effects were taken into account) to determine (among other factors) the confidence level in the AVA estimates that were calculated from field data. Again, all of this type of information will likely prove useful in characterizing and predicting reservoir occurrence and extent, etc.
    [0080] Returning next to Figure 4, this Figure contains a preferred operating logic, which would be suitable for use with the present invention. As a first step 400, preferably, seismic data will be obtained on a subsurface target of interest. Next, and preferably, the raw data will be edited, subjected to initial processing (eg geometry specification, filtering, etc.), which processes are intended to bring the data into condition for use in imaging.
    [0081] As a preferred next step, the explorer will conduct a velocity analysis (either automatic or manual) and build an initial velocity model (step 410) according to any number of methods well known to those of ordinary skill in the art.
    [0082] As a preferred next step 450, the data will be further processed for imaging purposes, which processing could include deconvolution, multiple removal, etc.
    [0083] Next, and preferably, the present invention will continue by updating the velocity model using, for example, seismic tomography or some other velocity estimation method (step 420).
    [0084] Next, the present invention will preferably continue by inserting additional overhead complexity into the velocity model (step 425). Since the initial model could be relatively simple (eg a flat or “layered pie” model), it is usually desirable to modify it to more accurately reflect the true subsurface layer configuration. Obviously, the more accurate the overloaded model, the more reliable the modeling results will be. In the case of the overloaded model, such additional information could come from the previously mentioned tomographic analysis (step 420) or from well logs, VSP prospects, from the seismic data itself, etc.
    [0085] As an example of the type of upgrade that is contemplated in step 425, where a salt dome (or other structure) is expected or known to be present below the seismic prospecting area, velocities that are typical of this type of structure could be added to the initial velocity model, with the idea that the additional complexity will at least roughly mimic this subsurface finding under prospect. Note that it is expected that this step could be performed where there is limited knowledge of the current configuration and extent of a salt dome or other structure. Thus, further refinements may be necessary or desirable.
    [0086] As a preferred next step 430, the present invention will obtain a final velocity model, preferably to refine the updated model from the previous step.
    [0087] Next, the present invention will preferably continue as the explorer has to provide an interpretation of the events in the migrated section. The purpose of this is to obtain information that can be used to further refine the velocity model in preparation for step 440, which is a perfect cluster/calibration debugging, or alternatively, synthetic data modeling.
    [0088] Next, and preferably, the present invention will continue to re-emigrate the demigrated pool (step 445).
    [0089] As a preferred next step, grouping weights will be calculated (step 450). Preferably, these weights will be based on the amount of seismic energy illuminating subsurface reflectors that are recorded within a trace. In other words, the amplitudes within the traces produced by the migration/modeling/disfiguration of the model subsurface will tend to have higher values where a reflector has well-formed images and lower (or zero) values where there is less reflected seismic energy. This suggests that seismic traces in the original prospect could be weighted according to the energy contained in the uprooted/modeled/migrated model data, as such energy represents the amount of illumination, to which the model trace was exposed.
    [0090] In the following, depending on the wishes of the explorer, the present invention could continue along two different lines. In some cases, an AVO analysis will be performed (step 480), as indicated in greater detail in Figure 5.
    [0091] In other preferred embodiments, improvement of the invention 460 will be performed further and an initial migrated image prepared (step 465).
    [0092] In some preferred embodiments, an enhanced image will be prepared using the final velocity model and clustering weights calculated from step 450 (step 470). Naturally, the enhanced image could then be used to search for subsurface features according to methods well known to those of ordinary skill in the art.
    [0093] Figure 5 contains additional details of a preferred embodiment of step 480 of Figure 4. As a first step, data from step 405 will preferably be migrated with a final velocity model of the type obtained from step 430. Naturally, the migrated data is useful without further processing at this point and the data from step 505 could be used in any number of ways. The steps that follow in Figure 5 are intended to be examples, rather than limitations of uses, to which this invention could be applied.
    [0094] Steps 510, 515, 520 (i.e., the left branch) illustrate how data processed by the present invention could be used in the form of a conventional AVA/AVO analysis (step 510), where weights from step 450 could be used to create an improved AVA/AVO analysis (step 515) and calculation of an AVA/AVO confidence map (step 520). Generally speaking, steps 510 - 520 would be the most frequently used in scanning.
    [0095] On the other hand, steps 525 - 545 contain applications of the present invention that would be more appropriate when a specific reservoir or other target has been selected. For example, steps 525 - 535 are concerned with using data that have been processed in accordance with the present invention in determining reservoir attributes, formulating a reservoir model, and modeling ahead of this reservoir model to create synthetic seismic data. . Such steps are, of course, ancient and well known in the art. However, step 540 uses the clustering weights from step 450 to create improved synthetic seismic data sets and step 545 could come into play when the explorer compares the improved synthetic AVA/AVO data with the current subsurface response when measured. by the seismic data.
    [0096] It would be prudent in some cases to design an AVA assessment workflow that quantitatively analyzes how lighting issues impact a detailed seismic reservoir response. In general, once the lighting dataset has been obtained, different proposals can be used to integrate the lighting response with standard reservoir modeling workflows, depending on the seismic analysis objective.
    [0097] For example, in some preferred embodiments, a detriment matching filter could be obtained from the lighting data with the intention that this will be applied to a detailed reservoir response. In this workflow, an initial reservoir model could be created to determine the AVA response in the absence of any lighting effects. Since the lighting volume contains only lighting effects, these effects can be transferred to the reservoir model response through the use of a matching filter. The filter would be designed in such a way that its application to a calibration cluster such as the one taught here (ie with no AVA signature) would reproduce the lighting response. The filter would then be applied to the modeled response to the reservoir, thus transferring the lighting effect to the modeled response. Consequently, by having a “perfect”/calibration dataset (i.e., one with no AVA signature), and one with harmful effects from uneven lighting, it would be possible to design match filters that would locally transfer lighting degradation. for the user's choice of synthetically created dataset. conclusions
    [0098] The present exhibit presented a method for evaluating lighting effects as a function of aperture angle below complex sterile deposits. Through the use of matching filters, these lighting effects can be transferred to modeled AVA responses, leading to an effective method for evaluating lighting effects on AVA and the associated risk in well mapping and placement.
    [0099] Although the term “unit amplitude” model (that is, one that has unit amplitude reflectivity at all angles), has been used repeatedly here to characterize the seismic model that is used to calculate lighting weights, those of common knowledge in the art will recognize that all of the model amplitudes need not be equal to unity, although this could conventionally be done. In general, a “unit amplitude” model should be broadly understood to be an “equal amplitude” model, where the reflection amplitudes at all angles are all at least approximately equal and could be unitary or have some other constant value depending on the angle. of the wants and needs of the processor or programmer. Furthermore, those of ordinary skill in the art will recognize that a unit-amplitude reflectivity model can be constructed in many ways.
    [00100] Additionally, it should be noted that an important aspect of the present invention is the use of a calibration seismic dataset to normalize current seismic data to correct for imperfections in lighting so that a subsequent AVA analysis is more reliable. In the preferred embodiment, the calibration dataset will be obtained by migrating/disfiguring a unit amplitude model or by inserting density reflectors into a seismic model. That being said, what is important is to create a synthetic dataset that has lighting effects but no AVO effects, however this dataset could be calculated.
    [00101] Note that when the term “AVA” is used here, this term should be broadly interpreted to also include “AVO” analyses, although typically these are often considered to be a subset of the former.
    [00102] Output from the present invention could be used to create enhanced image or cluster volumes, calculate AVO/AVA confidence surfaces or volumes, perform enhanced AVA modeling (preferably using feed from the property modeling module of rock) etc.
    [00103] Finally, any or all of the foregoing could be used to characterize the seismic expression of the target reservoir. Seismic attributes calculated from the reservoir model could then be searched for and identified in the current seismic data as a means of mapping the extent of the reservoir and/or determining the presence or absence of hydrocarbons.
    [00104] While the invention set forth herein has been discussed almost exclusively in terms of seismic traces that can be organized into trigger, receiver, or “CMP” groupings, this has been done for specificity purposes only and not with any intent to limit the present invention for operation on this type of seismic data only. Thus, within the text of this exposition, the terms seismic trace and CMP cluster are used in the broadest possible sense of those terms, and they are understood to apply to conventional 2D and 3D traces and conventional CPM clusters, as well as to others. types of clusters which could include, without limitation, CRP clusters, CCP clusters (i.e. “common conversion point” clusters), CACP clusters (“asymptotic common conversion point”) clusters, common shift clusters, clusters trigger/receiver etc., the most important aspect of a “cluster” being that it represents an organized collection of unstacked seismic traces from either a 2D or 3D prospect, all of which have at least one point of common subsurface image. Note that the term migration is used in its broadest sense consistent with the goals of the present invention, which could include time or depth Kirchhoff migration, beam migration, one-way migrations, or two-way acoustic or elastic migration. lanes, with any variation of extrapolated data, including bridge source, in-line source, or plane wave extrapolations.
    [00105] As mentioned previously, typical seismic surveys can be thought of as consisting of five dimensions, ie time (or z), x, y, shot-receiver offset in x, and shot-receiver offset in y. This method taught here can easily be extended from two dimensions to three, four, five, six, or even more dimensions by those of common knowledge in the art (e.g., a 6D prospect is a 3D space-lapse prospect). pre-stacked time). Using more dimensions typically allows for a more accurate and robust interpolation method, although the computation cost tends to increase with dimensionality.
    [00106] Yet, in the previous discussion, the language was expressed in terms of operations performed on conventional seismic data. However, it is understood by those skilled in the art that the invention described herein could be advantageously applied to other subject areas, and used to locate subsurface minerals other than hydrocarbons. By way of example only, the same proposal described here could potentially be used to process and/or analyze multi-component seismic data, shear wave data, converted mode data, cross well prospecting data, full waveform sonic profiles , controlled source or other electromagnetic data (CSEM, t-CSEM etc.) or model-based digital simulations of any of the foregoing. Additionally, the methods claimed hereinbelow can be applied to mathematically transformed versions of these same data traces including, for example: filtered data traces, migrated data traces, Fourier frequency domain transformed data traces; transformations by discrete orthogonal transforms; instantaneous phase data traces, instantaneous frequency data traces, quadrature traces, analytical traces; etc. In summary, the process exposed here can potentially be applied to a wide variety of geophysical time series types, but is preferably applied to a collection of spatially related time series. Thus, when the term “seismic trace” or “trace” is used here, those terms should be broadly interpreted to include traditional seismic traces as well as any of the foregoing.
    [00107] While the inventive device has been described and illustrated herein with reference to certain preferred embodiments in relation to the accompanying drawings, various alterations and other modifications, apart from those shown or suggested herein, may be made therein by those skilled in the art, without departing from the spirit of the invention concept, the scope of which must be determined by the following Claims.
    权利要求:
    Claims (1)
    [0001]
    1. Method of Exploration of Hydrocarbons in the Interior of a Predetermined Volume of the Earth Containing Structural and Stratigraphic Characteristics Conducive to the Generation, Migration, Accumulation or Presence of the aforementioned Hydrocarbons, characterized in that it comprises the steps of: a. access a digital representation of a seismic survey that forms the images of at least a part of said predetermined volume of the earth; B. create a subsurface model that is at least representative of the predetermined volume of land; ç. using at least said subsurface model to create a calibration seismic dataset, said calibration seismic dataset being formed from a seismic modeling program comprising a lighting modeling calculation, wherein the creation of a calibration dataset comprises: (i) demigrating (440), or applying forward modeling to the subsurface model to generate modeled data based on data comprising a plane reflectivity signature, where the plane reflectivity signature corresponds to a signature in which there is no variation in the reflection amplitude with a change incidence angle and (ii) remigrating (445) the modeled data to the subsurface model to provide the calibration seismic dataset; d. using at least a portion of said set of seismic calibration data to normalize at least a portion of said digital representation of said seismic survey; and. conducting an angle of incidence analysis (AVA) (515) using at least a portion of said normalized portion of said digital representation of said seismic survey; and f. using said AVA analysis to explore hydrocarbons within said predetermined volume of land.
    类似技术:
    公开号 | 公开日 | 专利标题
    BR112012006931B1|2022-02-15|Method of exploration for hydrocarbons within a predetermined volume of land containing structural and stratigraphic characteristics leading to the generation, migration, accumulation or presence of said hydrocarbons
    US7952960B2|2011-05-31|Seismic imaging with natural Green's functions derived from VSP data
    US6904368B2|2005-06-07|Seismic analysis using post-imaging seismic anisotropy corrections
    US6826484B2|2004-11-30|3D prestack time migration method
    US9829592B2|2017-11-28|Seismic imaging with visco-acoustic reverse-time migration using pseudo-analytical method
    EP1865340B1|2010-09-29|A process and program for characterising evolution of an oil reservoir over time
    EP1611461B1|2013-05-29|Method for simulating local prestack depth migrated seismic images
    US6128580A|2000-10-03|Converted-wave processing in many-layered anisotropic media
    BRPI0918020B1|2019-09-10|seismic exploration methods
    EP1047957A1|2000-11-02|Method of prestack 3-d migration
    EP2548052B1|2019-10-30|System and method of 3d salt flank vsp imaging with transmitted waves
    BR112012021391A2|2020-08-25|methods of geophysical exploration within predetermined volume of land
    AU2014227491B2|2019-06-06|Construction and application of angle gathers from three-dimensional imaging of multiples wavefields
    AU2015337052A1|2017-05-04|Structure tensor constrained tomographic velocity analysis
    Laurain et al.2004|Simulated migration amplitude for improving amplitude estimates in seismic illumination studies
    Chen et al.2020|Machine Learning Method to Determine Salt Structures from Gravity Data
    Buske et al.2009|Imaging and inversion—Introduction
    Rathy2020|Separation of Diffractions from Reflections by the Double Square Root Operator
    BR112014024617B1|2021-12-21|METHODS AND SYSTEMS FOR OPTIMUM STACKING OF SEISMIC DATA
    Jarillo Michel0|Elastic waveform inversion of microseismic data in anisotropic media
    Purnell et al.2002|Wave-equation migration of common-receiver gathers for AVO analysis
    WO2014198347A1|2014-12-18|Method and apparatus for determining rock properties
    同族专利:
    公开号 | 公开日
    EA201200460A1|2012-09-28|
    EP2483712B1|2021-07-07|
    CN102597809A|2012-07-18|
    BR112012006931A2|2020-12-15|
    CA2775561A1|2011-04-07|
    CN102597809B|2016-06-08|
    AU2010300304B2|2014-07-10|
    EP2483712A1|2012-08-08|
    AU2010300304A1|2012-04-26|
    EA027554B1|2017-08-31|
    IN2012DN03212A|2015-10-23|
    US8553499B2|2013-10-08|
    MX2012003959A|2012-05-08|
    WO2011041782A1|2011-04-07|
    US20110103187A1|2011-05-05|
    引用文献:
    公开号 | 申请日 | 公开日 | 申请人 | 专利标题
    US4646239A|1984-07-20|1987-02-24|Standard Oil Company|Method of seismic exploration including processing and displaying seismic data to quantitatively distinguish among seismic events|
    US5838634A|1996-04-04|1998-11-17|Exxon Production Research Company|Method of generating 3-D geologic models incorporating geologic and geophysical constraints|
    US6999879B2|2003-05-14|2006-02-14|Exxonmobil Upstream Research Company|Method for controlling seismic coverage using decision theory|
    FR2858063B1|2003-07-21|2005-09-16|Geophysique Cie Gle|METHOD FOR ESTIMATING THE ILLUMINATION COVERAGE RATE IN THE MIGRE DOMAIN|
    US7280918B2|2005-08-08|2007-10-09|Knowledge Systems, Inc.|Method and system for combining seismic data and basin modeling|
    CN101581791A|2008-05-15|2009-11-18|中国石油天然气集团公司|Method for extracting pre-stack hydrocarbon detection section of seismic data with low signal-to-noise ratio|AU2011261849B2|2010-06-02|2014-10-23|Exxonmobil Upstream Research Company|Efficient computation of wave equation migration angle gathers|
    US8301385B2|2010-06-29|2012-10-30|Chevron U.S.A. Inc.|Shot gather data beamer and debeamer|
    AU2012212520B2|2011-01-31|2014-09-11|Chevron U.S.A. Inc.|Extracting geologic information from multiple offset stacks and/or angle stacks|
    US9625593B2|2011-04-26|2017-04-18|Exxonmobil Upstream Research Company|Seismic data processing|
    US9772415B2|2011-08-05|2017-09-26|Saudi Arabian Oil Company|Correcting time lapse seismic data for overburden and recording effects|
    US9176930B2|2011-11-29|2015-11-03|Exxonmobil Upstream Research Company|Methods for approximating hessian times vector operation in full wavefield inversion|
    CA2858378A1|2011-12-20|2013-06-27|Conocophillips Company|Critical reflection illumination analysis|
    US10977396B2|2012-01-13|2021-04-13|Schlumberger Technology Corporation|Determining an elastic model for a geologic region|
    US10162070B2|2012-04-05|2018-12-25|Westerngeco L.L.C.|Converting a first acquired data subset to a second acquired data subset|
    US9448316B2|2012-04-19|2016-09-20|Cgg Services Sa|Seismic data processing including predicting multiples using a reverse time demigration|
    US9665604B2|2012-07-31|2017-05-30|Schlumberger Technology Corporation|Modeling and manipulation of seismic reference datumin a collaborative petro-technical application environment|
    US9268060B2|2013-03-14|2016-02-23|Bp Corporation North America Inc.|System and method for computational geology|
    US9562983B2|2014-04-17|2017-02-07|Saudi Arabian Oil Company|Generating subterranean imaging data based on vertical seismic profile data|
    US10267937B2|2014-04-17|2019-04-23|Saudi Arabian Oil Company|Generating subterranean imaging data based on vertical seismic profile data and ocean bottom sensor data|
    WO2016001750A2|2014-06-30|2016-01-07|Cgg Services Sa|Seismic data processing using matching filter based cost function optimization|
    US20170139067A1|2014-07-03|2017-05-18|Westerngeco Llc|Parallel processing seismic wavefield data|
    CA2961572C|2014-10-20|2019-07-02|Exxonmobil Upstream Research Company|Velocity tomography using property scans|
    MX2017007988A|2015-02-17|2017-09-29|Exxonmobil Upstream Res Co|Multistage full wavefield inversion process that generates a multiple free data set.|
    WO2017066072A1|2015-10-15|2017-04-20|Bp Corporation North America Inc.|Interactive image weighting by illumination in seismic imaging|
    US10634804B2|2015-12-21|2020-04-28|Chevron U.S.A. Inc.|System and method for dip-guided seismic image stacking|
    US10739481B2|2016-03-30|2020-08-11|Exxonmobil Upstream Research Company|2D multiline seismic reflection tomography with seismic-tie constraint|
    US20180059276A1|2016-08-25|2018-03-01|Waveseis LLC|System and method for focusing seismic images|
    CA3062544A1|2017-05-09|2018-11-15|Chevron U.S.A. Inc.|System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic data|
    US20210311223A1|2020-04-07|2021-10-07|Chevron U.S.A. Inc.|System and method for seismic inversion|
    法律状态:
    2020-12-29| B06F| Objections, documents and/or translations needed after an examination request according [chapter 6.6 patent gazette]|
    2021-07-06| B07A| Application suspended after technical examination (opinion) [chapter 7.1 patent gazette]|
    2021-08-31| B07C| Technical examination (opinion): republication [chapter 7.3 patent gazette]|Free format text: REPUBLIQUE-SE POR INCORRECAO NO QUADRO 1 |
    2021-12-07| B350| Update of information on the portal [chapter 15.35 patent gazette]|
    2022-01-11| B09A| Decision: intention to grant [chapter 9.1 patent gazette]|
    2022-02-15| B16A| Patent or certificate of addition of invention granted [chapter 16.1 patent gazette]|Free format text: PRAZO DE VALIDADE: 20 (VINTE) ANOS CONTADOS A PARTIR DE 04/10/2010, OBSERVADAS AS CONDICOES LEGAIS. PATENTE CONCEDIDA CONFORME ADI 5.529/DF, QUE DETERMINA A ALTERACAO DO PRAZO DE CONCESSAO. |
    优先权:
    申请号 | 申请日 | 专利标题
    US24822209P| true| 2009-10-02|2009-10-02|
    US61/248,222|2009-10-02|
    PCT/US2010/051321|WO2011041782A1|2009-10-02|2010-10-04|Migration-based illumination determination for ava risk assessment|
    [返回顶部]