RU2577792C1 - Robust process for depth imaging in seismic survey based on operator adjustment by reference seismic record - Google Patents

Abstract

FIELD: physics.

SUBSTANCE: claimed process comprises seismic image plotting by reflected waves parameters and those of the medium. Note here that the amplitude of signals is balanced by automatic signal level control for processing of the initial seismic record. Additional processing of the conversion results in areas of locus functions intersection with the help of pie slice is executed with nonlinear adaptation of weight factors, hence, introducing the reference seismic record. Weight factors of migration conversion error-controlled operator are calculated for its application to reference seismic record. Migration operator resulted from application of the reference seismic record is applied to initial data or migration conversion is used to isolate the interference for it to be adaptively subtracted from the result of traditional migration.

EFFECT: higher accuracy of data.

2 dwg

Description

FIELD OF THE INVENTION

The invention relates to technologies for constructing deep seismic images of fragments of the earth's crust, and in particular to technology that allows high-precision assessment of deep-speed models, as well as to obtain high-quality and accurate dynamic deep images of the environment from seismic data in the process of prospecting and exploration of mineral deposits.

State of the art

The prior art algorithm for noise-free Radon transform ([1] Denisov, M., and Finikov, D., 2002, An alias protection scheme for Radon transform: 64th Ann. Internat. Mtg. EAGE). A method for suppressing the aliasing effect during the accumulation of seismic traces during the Radon transform is proposed in the paper (Alaizing is a hindrance arising from non-phase spatial summation of discrete signals). The algorithm is implemented in the frequency domain and uses both alternating and alternating summation. The disadvantages of the algorithm are manifested in the conditions of variability of the wave pattern along the time coordinate. The use of the Fourier transform eliminates the possibility of tuning the algorithm to local conditions of non-phase accumulation of signals.

The prior art algorithm for suppressing artifacts when summing seismic traces ([2] Denisov MS, Finikov DB, 2005. A method for suppressing sampling noise when summing seismic traces (using the example of multiple wave modeling), Geophysics, 1). In the work, a method for suppressing the aliasing effect of accumulation of traces is proposed, the problem of predicting the repeatedly reflected interference waves is considered. The algorithm is implemented in the space-time domain and uses both sign-constant and sign-alternating summation. The disadvantages of the algorithm are manifested in the conditions of signal interference, especially when crossing intense and weak waves.

Known seismic migration by shifting seismic logging data ([3] US 6002642 A, publ. 14.12.1999). The present invention relates to a geophysical exploration method that improves the accuracy of seismic migration. In the document under consideration, the inverse kinematic problem is solved by travel times, including tomographic methods. A method for migrating seismic data using bias of the surveyed measurements of seismic logging includes a wave path of similar migration for seismic data, used to determine the direct wave travel time for receivers in the well. Embodiments of the present invention provide for the direct use of wave travel time in migration or the reverse use of wave travel time in migration through the construction of a migration velocity model. Variants of the velocity model for the wave travel time provide either error correction through the use of interpolated error functions or the construction of migration error tables. The invention can be used for the wave travel time, Kirchhoff depth and migration, in two or three quantities, and in any summation. The invention can be used to carry any type of seismic data, including longitudinal wave, shear wave, and wave seismic data conversion.

A known method for determining the optimal rate of temporary migration before summation ([4] CN 101839999, publ. 09/22/2010). The method includes determining a seismic wave near the ground. A detector is located on the ground to collect reflected seismic waves, i.e. the detector collects seismic data. In this case, a change in the time velocity of the reflected seismic wave is used as the optimal speed. Using the speed interval that varies with the reflection time of the seismic wave, the step length is changed to form a group of velocities of functional objects, and temporary migration is performed until the accumulated seismic data is summed. Describe the result of changes in speed in the speed range, in accordance with the intensity of the group waves and the distribution of the spectrum of occurrence of speeds. The optimal speed is taken before summing the temporary migration. The intersecting temporal migration velocities prior to summation, determined at a common reflective point, are used to obtain a new velocity profile. This invention calculates the information that should be output at the time of the speed analysis in the effective speed range, obviously, reduces the number of calculations and improves practicality.

A known method of seismic exploration for the study of sedimentary cover in the presence of strongly indented acoustically hard boundaries ([5] RU 2221262, publ. 10.01.2004). The invention relates to geophysics, to seismic methods of mineral exploration and is intended to obtain seismic sections of high resolution and reliability in difficult geological conditions. The seismic exploration method includes conducting seismic exploration, processing the obtained initial seismic data, and constructing a seismic time section. To account for the non-hyperbolicity of the hodographs, additional processing of the initial seismograms is carried out by paleotransforming them in the time domain to the conditions of the actual existing paleo-boundary for the period of the paleogeographic situation, when the border lies subhorizontal. The paleotransformation is carried out by introducing pseudo-paleostatic corrections, their correction, carried out by selecting effective paleoscale speeds (Veff paleo). According to paleotransformed seismograms, regular waves - interference are suppressed and, taking into account data on interval velocities, a paleotemporal migrated section is constructed.

The closest analogue of the proposed technical solution is the patent ([6] RU 2126984 "Method for determining the depth-velocity parameters of the medium and constructing its image from seismic data - the prime system", 02/27/1999). This invention relates to the field of seismic exploration and can be used to determine the geological parameters of the environment and its deep image for the search for oil and gas deposits in complex seismic and geological conditions. To improve the resolution of seismic records, increase the signal-to-noise ratio and increase the reliability of seismic imaging, the kinematic parameters of the reflected waves are determined by using the local conversion operator to accumulate seismograms, and the depth-velocity parameters of the medium are determined by checking the adequacy of the selected medium model and real data by solving the inverse problem in two ways, one of which uses boundary conditions on the roof, and the other on the bottom seam, and comparing these results with each other.

SUMMARY OF THE INVENTION

The problem solved by the claimed invention is to make a high-precision assessment of deep-speed models, as well as to obtain high-quality and accurate dynamic deep images of the medium from seismic data in the process of prospecting and exploration of mineral deposits, under the conditions of interference of reflected waves.

The technical result of the invention consists in increasing the accuracy of estimation of deep-speed models by 7-18% compared with standard methods for simple seismic and geological conditions and by 45-80% for complex, as well as increasing the accuracy of dynamic deep-seated images by 12-23% for simple seismic and geological conditions and 60-75% for complex ones.

To achieve the specified technical result in the proposed method for constructing deep images in seismic exploration based on the operator’s settings using reference seismograms, including the construction of a seismic image from seismic data, namely, the parameters of the reflected waves and environmental parameters, in order to process the initial seismogram, the signal-wave amplitude is equalized by the procedure automatic level control (AGC) of the signal; carry out additional processing of the conversion results in the regions of hodographs intersection using fan filtering with non-linear adaptation of weight coefficients, thereby introducing a reference seismogram; calculate the weighting coefficients of the noise-immigration migration transformation operator, applying it to the reference seismogram; apply the migration operator obtained from the reference seismogram to the source data or apply the migration transformation in order to isolate the interference, followed by its adaptive subtraction from the result of traditional migration.

Brief Description of the Drawings

Fig. 1. A flowchart of processing based on the construction of a reference seismogram.

Fig. 2. The block diagram of the noise-free migration conversion based on the allocation of interference and its adaptive subtraction.

Disclosure of invention

Under the conditions of spatial discreteness of seismic data during the implementation of multichannel procedures, which include migration transformation of seismograms, artifacts appear, which are commonly called aliasing interference.

Seismic exploration (seismic exploration) allows you to obtain a structural map of the soil by emitting downward acoustic or elastic waves into the soil and registering "echo signals" reflected from the underlying layers of the rock. For the emission of downward acoustic or elastic waves into the ground, for example, explosions or seismic vibrators on the ground and air guns at sea can be used. In the process of seismic exploration, the wave emitter is moved along the soil surface above the investigated geological structure. Each time the emitter is excited, a downward seismic signal is generated, which propagates through the ground, is reflected and / or diffracted, and after reflection is recorded at many points on the surface. At the same time, numerous combinations of excitation and registration of sources are combined to create an almost continuous section profile, which can extend several hundred kilometers. In two-dimensional seismic exploration, the positions of the radiation and registration sources are usually located on one straight line, and in three-dimensional seismic exploration, the positions of the radiation and registration sources are usually distributed over the surface at nodes of the coordinate grid. Simply put, one can imagine that a 2-dimensional seismic profile provides an image of the cross section of soil layers with reflecting horizons located approximately in the middle between the positions of the radiation sources and geophones. Three-dimensional seismic exploration provides a “cube” of data, that is, in principle, a three-dimensional image of the geological environment under the survey area, with reflecting horizons located approximately in the middle between the positions of the radiation sources and geophones located in the nodes of the data acquisition grid.

The migrated trace (seismic trace) M (q, y, t) of the common point of excitation (OPV) seismogram is obtained as follows, where y is the lateral coordinate, q is the coordinate of the oscillation source. The point of excitation (Ndp. Point of excitation) is the point on, above or below the surface of observations, in which the excitation of seismic oscillations is performed.

Calculations are made according to the formula:

where u (q, x, t) is a seismogram, parameters A and B determine the spatial summation aperture, i.e. interval x∈ [aA, a + B] on the observation profile. The asterisk denotes the convolution procedure, f (t) is the compensating filtering operator, w (y, x) is the weighting factor, t (y, x) is the summation trajectory. Expression (1) is converted to the form

отличается от и применением фильтрации, умножением на вес и подвижкой. Введем понятие локальной суммы, полученной на базе (2К+1) каналов,Where $$\tilde{u}$$

differs from the application of filtration, multiplication by weight and movement. We introduce the concept of a local sum obtained on the basis of (2K + 1) channels,

the global sum can be obtained both by direct summation within the aperture and by adding local sums

When summing (4), there are local areas that provide signal accumulation (constructive summation), and areas in which only aliasing is formed. Therefore, the task of constructing an aliasing noise suppression algorithm is formulated as a selection of a criterion for separating these areas.

, т.е. отношение энергий локальной суммы (3), которую в дальнейшем для удобства мы будем обозначать как b⁺, и аналогичной суммы, b^-, полученной со знакопеременной весовой функцией:As such a criterion, we choose

, i.e. the ratio of the energies of the local sum (3), which for convenience will be denoted below as b ⁺ , and a similar sum, b ^- obtained with an alternating weight function:

where the variable ζ means the channel number, i.e. integer. We construct the function w (A): w (0) = 1 and w (A) → 0 as A → ∞, and moreover, w (A) is close to 1 as A∈ [0, p], then we present the aliasing noise suppression algorithm in form

So, we get the following algorithm:

1. Local summation in order to obtain b ⁺ and b ^- ;

2. Estimation of energies in a window sliding over t;

по координате x и получение искомой мигрированной трассы М.3. Summation $$\tilde{b}$$

along the x coordinate and obtaining the desired migrated route M.

In the interference of waves whose amplitudes are significantly different, such an algorithm reveals a completely natural tendency to tune in to processing an intense wave, while ignoring a weak signal. In order to overcome this situation, it seems advisable to introduce a “reference” seismogram by which you can configure the algorithm, i.e. assignment of weight coefficients, and then apply the obtained coefficients to the seismogram that needs to be processed. A flowchart of a processing procedure based on the construction of a reference seismogram is shown in Fig. one.

Consider the method of calculating the reference seismogram. The traditional way to equalize signal amplitudes is the AGC procedure (automatic level control). It is known that such a conversion successfully copes with the task in areas where there is no interference of signals, while in the vicinity of the intersection points of the hodographs, the signal amplitudes undergo local distortions. It is clear that such a property of the AGC algorithm makes it unsuitable for our task. Therefore, we propose to perform additional processing of the conversion result in order to equalize signal amplitudes in the regions of hodographs intersection. Since the distortion of the signal energy on the hodograph is local in nature, it can be eliminated by fan filtering with non-linear adaptation of weight coefficients. The expression for the transformation carried out by such a kinematic filter can be written as follows

where u (x, t) is the initial wave field, d ^α (x, t) is the directional sum, y∈ [-L, L] is the spatial base of the filter, α is the summation direction sorted through the fan, U (x, t ) is the result of filtering.

Since hodographs of waves at this stage of processing are considered unknown, difficulties will arise in estimating the kinematics of the signal, which we will eliminate by using an adaptive filter, which at each point (x, t) independently analyzes the wave field for the presence or absence of a signal or interference. Suppose that at some point there is a fan centered with respect to the tangent to the a priori hodograph. In order to decide whether a given point really belongs to reflection, the coherence of the wave field in all directions α inside the fan is calculated. According to the meaning of the problem under consideration, the fan is chosen deliberately wide. As a measure of coherence, we take the value of the semblance, that is, a parameter calculated by the formula

and accepting values from zero to one. To obtain smoother values of the samples as a function of x and t, averaging is possible within a moving window. We introduce a smooth weight function taking values close to unity for close to unity p ^α (x, t) and close to zero for close to zero p ^α (x, t). This may be, for example

where n is an integer, Q is the threshold value for deciding on the absence or presence of a coherent signal. Then we write the adaptive kinematic filtering algorithm as

where d ^α (x, t) is calculated by formula (4). In the absence of a priori information about the kinematics of the signal, this method will eliminate local disturbances in the signal dynamics caused by the application of the AGC procedure and prepare a reference seismogram for calculating the weight coefficients of the noise-protected migration conversion operator.

(1 вместо 0 и наоборот: $$\tilde{w}=1-w$$

), тогда вместо выделения зон конструктивного накапливания сигнала и подавления зон несинфазного суммирования будет наблюдаться обратный эффект. На выходе преобразования получим помеху. Одновременно с этим применим традиционную реализацию (не помехозащищенную) миграции, в результате чего имеем сигнал на фоне помехи. Теперь геофизик сможет воспользоваться всем арсеналом имеющихся у него средств адаптивного вычитания регулярных и нерегулярных помех, включая неквадратичные нормы, латеральное усреднение критерия, многооконность, многоканальность, нестационарность и т.д.It is also proposed to supplement the described method for obtaining deep images with the possibility of outputting not an anti-interference image, but the result of evaluating aliasing interference. For this purpose, in the conversion, we need to change the weight coefficients w by their inverse values $$\tilde{w}$$

(1 instead of 0 and vice versa: $$\tilde{w}=\mathrm{one}-w$$

), then instead of isolating the zones of constructive accumulation of the signal and suppressing the zones of non-phase summation, the opposite effect will be observed. At the output of the conversion, we get interference. At the same time, we apply the traditional implementation of (non-interference-protected) migration, as a result of which we have a signal against a background of interference. Now the geophysicist will be able to use the whole arsenal of the means of adaptive subtraction of regular and irregular noise available to him, including non-quadratic norms, lateral averaging of a criterion, multi-window, multi-channel, non-stationary, etc.

A block diagram of the indicated processing sequence, including the step of isolating interference with its subsequent adaptive subtraction, is shown in Fig. 2.

Information sources

1. Denisov, M., and Finikov, D., 2002. An alias protection scheme for Radon transform: 64th Ann. Internal Mtg. EAGE.

2. Denisov MS, Finikov DB, 2005. A method of suppressing sampling noise when summing seismic traces (using the example of modeling multiple waves), Geophysics, 1.

3.US 6002642 A, Seismic migration using offset checkshot data, 12/14/1999.

4. CN 101839999 A, Method for determining optimum velocity section for pre-stack time migration, 09/22/2010.

5. RU 2221262 C1. Seismic exploration method for studying sedimentary cover in the presence of severely indented acoustically hard boundaries (options), 01/10/2004.

6. RU 2126984 C1. “The method for determining the deep-speed parameters of the medium and constructing its image from seismic data is the prime system,” 02.27.1999, (Prototype).

Claims (1)

A sustainable method for constructing deep images in seismic exploration based on the operator’s settings using reference seismograms, including the construction of a seismic image from seismic data, namely, reflected wave parameters and environmental parameters, characterized in that, in order to process the initial seismogram, the signal amplitude is equalized by the automatic signal level control procedure, carry out additional processing of the conversion results in the areas of intersection of hodographs using a fan filter non-linear adaptation of the weight coefficients, thereby introducing the reference seismogram, calculate the weight coefficients of the noise-protected migration conversion operator, applying it to the reference seismogram, apply the migration operator obtained from the reference seismogram to the source data or apply the migration transformation to isolate the interference with its subsequent adaptive subtraction from the result of traditional migration.

Images