RU2126984C1 - Method determining depth and speed parameters of medium and construction of its image by seismic data- prime system - Google Patents
Method determining depth and speed parameters of medium and construction of its image by seismic data- prime system Download PDFInfo
- Publication number
- RU2126984C1 RU2126984C1 RU98113007A RU98113007A RU2126984C1 RU 2126984 C1 RU2126984 C1 RU 2126984C1 RU 98113007 A RU98113007 A RU 98113007A RU 98113007 A RU98113007 A RU 98113007A RU 2126984 C1 RU2126984 C1 RU 2126984C1
- Authority
- RU
- Russia
- Prior art keywords
- medium
- parameters
- depth
- model
- account
- Prior art date
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Изобретение относится к области сейсморазведки и может быть использовано для определения геологических параметров среды и ее глубинного изображения для поиска месторождений нефти и газа в сложных сейсмогеологических условиях. The 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 fields in difficult seismic and geological conditions.
Известен способ подавления кратных волн на основе их моделирования по зарегистрированному волновому полю заданной модели среды [1]. A known method of suppressing multiple waves based on their modeling on the registered wave field of a given model of the medium [1].
Данный способ описывает построение поля кратных волн на основе волнового уравнения в простых моделях среды без учета преломления. This method describes the construction of the field of multiple waves based on the wave equation in simple models of the medium without taking into account refraction.
Наиболее близким к предложенному способу является способ определения глубинно-скоростных параметров среды и построения ее изображения по сейсмическим данным, включающий определение кинематических параметров отраженных волн и глубинно-скоростных параметров среды, миграцию временного разреза в рамках пластовой модели среды, учет негиперболичности годографов при построении динамического глубинного разреза, а также осуществление деконволюции, согласованной фильтрации и коррекции поглощения [2]. Closest to the proposed method is a method for determining the depth-velocity parameters of the medium and constructing its image from seismic data, including determining the kinematic parameters of the reflected waves and depth-velocity parameters of the medium, the migration of the time section within the reservoir model of the medium, and the non-hyperbolicity of hodographs when constructing dynamic deep section, as well as the implementation of deconvolution, consistent filtering and absorption correction [2].
В известном способе определяют глубинно-скоростные параметры пластовой модели среды и выполняют динамическую глубинную миграцию по исходным сейсмограммам. При этом, результат оценки параметров среды существенно зависит от выбора модели, но в способе не предусмотрена проверка адекватности выбранной модели реальной среде. Кроме того, миграция сейсмограмм так же опирается на построенную в системе модель и качество динамических глубинных построений существенно зависит от того, оказалась ли она удачной или нет. In the known method, the depth-velocity parameters of the reservoir model of the medium are determined and dynamic depth migration is performed according to the initial seismograms. At the same time, the result of the assessment of the environmental parameters substantially depends on the choice of the model, but the method does not provide for verification of the adequacy of the selected model for the real environment. In addition, the migration of seismograms also relies on the model built in the system and the quality of the dynamic deep structures depends significantly on whether it turned out to be successful or not.
В известном способе предусмотрено использование различных процедур фильтрации волновых полей с оцениванием амплитудных и фазовых спектров сигналов. При оценивании их в узких пространственно-временных интервалах, дисперсии погрешности оценок могут быть весьма велики. In the known method provides for the use of various procedures for filtering wave fields with the assessment of the amplitude and phase spectra of the signals. When evaluating them in narrow space-time intervals, the variance of the error of the estimates can be very large.
Техническим результатом является повышение точности оценки параметров, получение более высокой разреженности сейсмических записей, улучшение соотношения сигнал/помеха и, в целом, повышение достоверности построения сейсмических изображений. The technical result is to increase the accuracy of parameter estimation, to obtain a higher sparsity of seismic records, to improve the signal-to-noise ratio, and, in general, to increase the reliability of the construction of seismic images.
Для достижения указанного технического результата в предложенном способе определения глубинно-скоростных параметров среды и построения ее изображения по сейсмическим данным, система PRIMЕ, включающем определение кинематических параметров отраженных волн и глубинно - скоростных параметров среды, миграцию временного разреза в рамках пластовой модели среды, учет негиперболичности годографов при построении динамического глубинного разреза, а также осуществление деконволюции, согласованной фильтрации и коррекции поглощения, при определении кинематических параметров отраженных волн, используют локальный оператор преобразования для накапливания сейсмограмм, при определении глубинных скоростных параметров среды проверяют адекватность выбранной модели среды и реальных данных путем решения обратной задачи двумя способами и сопоставления их результатов между собой, миграцию выполняют послойно с учетом преломления лучей в рамках выбранной глубинной скоростной модели с учетом ее адекватности реальной среде и с возможностью управления параметрами в каждом слое, учет негиперболичности годографов при построении динамического глубинного разреза производят путем погружения исходных сейсмограмм в рамках предварительно проверенной на адекватность реальной среде глубинно-скоростной пластовой модели с учетом преломления лучей, а поле кратных волн определяют также в рамках выбранной пластовой модели, адекватной реальной среде, с помощью миграционного преобразования сейсмограмм с учетом преломления лучей и адаптивного вычитания одновременно нескольких типов кратных волн, осуществляют амплитудную и фазовую деконволюцию, а согласованную фильтрацию и коррекцию поглощения осуществляют путем параметризации амплитудных и фазовых спектров, причем параметры амплитудного спектра и поглощения оценивают одновременно, а параметры фазового спектра оценивают отдельно по неквадратическому критерию оптимизации. To achieve the specified technical result in the proposed method for determining the depth-velocity parameters of the medium and constructing its image from seismic data, the PRIME system, including determining the kinematic parameters of the reflected waves and depth-velocity parameters of the medium, the migration of the time section in the framework of the reservoir model of the medium, taking into account the non-hyperbolicity of hodographs when constructing a dynamic deep section, as well as the implementation of deconvolution, coordinated filtering and absorption correction, when the kinematic parameters of the reflected waves, they use the local transformation operator to accumulate seismograms, when determining the deep velocity parameters of the medium, the adequacy of the selected model of the medium and real data is checked by solving the inverse problem in two ways and comparing their results with each other, migration is performed layer by layer taking into account the refraction of rays within the selected deep velocity model, taking into account its adequacy to the real environment and with the ability to control parameters in each layer, the hodographs are hyperbolic in constructing a dynamic deep section by immersing the initial seismograms within the framework of a depth-velocity reservoir model previously checked for adequacy of the real medium taking into account the refraction of the rays, and the multiple wave field is also determined within the framework of the selected reservoir model that is adequate to the real medium using the migration transformation seismograms taking into account the refraction of rays and adaptive subtraction of several types of multiple waves simultaneously, carry out amplitude and phase deconvolution, and matched filtering and correction of absorption is carried out by parameterizing the amplitude and phase spectra, the parameters of the amplitude spectrum and absorption being evaluated simultaneously, and the parameters of the phase spectrum being evaluated separately according to the non-quadratic optimization criterion.
Анализ отобранных в процессе поиска известных технических решений показал, что в науке и технике нет объекта, обладающего заявленной совокупностью признаков и которые приводили бы к достижению более высокого технического результата, что позволяет сделать вывод о соответствии изобретения критериям "новизна" и "изобретательский уровень". Способ реализуют следующим образом:
1. Определяют кинематические параметры отраженных волн tо(x) и Vо(x).An analysis of the known technical solutions selected during the search showed that there is no object in science and technology that has the claimed combination of features and which would lead to a higher technical result, which allows us to conclude that the invention meets the criteria of "novelty" and "inventive step". The method is implemented as follows:
1. Determine the kinematic parameters of the reflected waves t o (x) and V o (x).
Для этого используют горизонтальные спектры скоростей суммирования, вычисляемые но накопленным сейсмограммам. Новым здесь является способ накопления сейсмограмм с помощью локального оператора. Он базируется на том, что анализируемая синфазность предварительно кореллируется па временном разрезе. Это позволяет без предположений о локальной линейности сражающей границы и постоянстве скоростей к слое между ней и дневной поверхностью (как это делается в известном операторе ДМО) вычислить положение точки накапливания сигналов. Горизонтальные спектры накопленных таким способом сейсмограмм оказываются значительно более устойчивыми и позволяют вычислять значения Vогт(x) даже в случае невысокого качества исходных данных (пропуск большого количества каналов, низкое отношение сигнал\помеха). После нахождения скоростей суммирования локальный оператор преобразования для накапливания сейсмограмм используется для построения участка временного разреза в окрестности анализируемой синфазности величиной 120-300 мс. Локальность оператора позволяет сделать оператор накапливания стационарным по времени, что в свою очередь снимает проблемы растяжения сигналов, имеющие место в существующих способах построения временного разреза.To do this, use the horizontal spectra of the summation velocities calculated but accumulated by seismograms. What is new here is the method of gathering seismograms using a local operator. It is based on the fact that the synphase being analyzed is preliminary correlated in a time section. This allows one to calculate the position of the signal accumulation point without assumptions about the local linearity of the striking boundary and the constancy of the velocities to the layer between it and the day surface (as is done in the well-known operator of the DME). The horizontal spectra of seismograms accumulated in this way turn out to be much more stable and make it possible to calculate the values of Vgt (x) even in the case of poor quality of the initial data (skipping a large number of channels, low signal / noise ratio). After finding the summation velocities, the local transformation operator for accumulating seismograms is used to construct a section of the time section in the vicinity of the analyzed common-mode value of 120-300 ms. The locality of the operator makes it possible to make the accumulation operator stationary in time, which in turn removes the problems of signal stretching that occur in existing methods for constructing a time section.
В результате существенно уточняется положение изучаемой линии tо(x) и следующая итерация, использующая новое tо(x), позволяет уточнить Vогт(x)- Две, три подобные итерации позволяют весьма точно определить кинематические параметры анализируемой волны. Этот же оператор используется при построении горизонтальных спектров скоростей для годографов вида (t
2. Определяют глубинно-скоростные параметры среды. 2. Determine the depth-velocity parameters of the environment.
Для вычисления пластовых скоростей и положения отражающих границ необходимо прежде всего задать модель, в рамках которой эти вычисления будут производиться. Сам способ расчета после этого играет вспомогательную роль (их существует множество: инверсные, топографические, оптимизационные и т. д. ). Результат существенно определяется именно выбором модели. При этом, в различных сейсмогеологических условиях модели должны быть разными. Поэтому недостаточно сформулировать общее описание модели (как это обычно делается), например, выбрать слоистую модель с постоянными пластовыми скоростями. Необходимо еще уметь проверять в каждом случае, насколько она адекватна реально обрабатываемому разрезу. Метод такой проверки разработан и используется в предлагаемом способе обработки сейсмических данных. В случае 2D наблюдений он основан на том, что обратная кинематическая задача переопределена по краевым условиям: их два - на кровле и подошве слоя, а в однородном слое для нахождения Vпл достаточно только одного. Поэтому, построены два способа решения обратной задачи, один - использующий только условие на кровле (способ взаимных точек), а другой - только на подошве (R-способ). Они применяются к одним и тем же данным и результаты сравниваются между собой. Если результаты совпадают между собой (в статистическом смысле), то выбранный слой является локально однородным и решение, полученное в этом предположении, является верным. Если результаты решения этими двумя способами не совпадают, то слой локально неоднороден и следует изменить модель (например, добавить слой). Вычисление этого критерия от слоя к слою и в каждой точке очередного исследуемого слоя по латерали позволяет построить глубинно-скоростную модель, адекватную реальной среде. В случае 3D наблюдений аналогичный критерий строится на основе R-способа и анализа поверхности tо(x,у).To calculate reservoir velocities and the position of reflecting boundaries, it is first necessary to specify the model within which these calculations will be performed. After that, the calculation method itself plays an auxiliary role (there are many of them: inverse, topographic, optimization, etc.). The result is essentially determined by the choice of model. Moreover, in different seismic and geological conditions, the models should be different. Therefore, it is not enough to formulate a general description of the model (as is usually done), for example, to choose a layered model with constant reservoir velocities. It is also necessary to be able to check in each case how adequate it is to the actually processed section. The method of such verification is developed and used in the proposed method for processing seismic data. In the case of 2D observations, it is based on the fact that the inverse kinematic problem is redefined according to boundary conditions: there are two of them on the roof and the bottom of the layer, and in a homogeneous layer only one is enough to find V pl . Therefore, two methods of solving the inverse problem are constructed, one using only the condition on the roof (the method of reciprocal points), and the other only on the sole (R-method). They apply to the same data and the results are compared with each other. If the results coincide (in a statistical sense), then the selected layer is locally homogeneous and the solution obtained under this assumption is correct. If the results of the solution in these two ways do not match, then the layer is locally heterogeneous and the model should be changed (for example, add a layer). Calculation of this criterion from layer to layer and at each point of the next studied layer in lateral allows constructing a depth-velocity model adequate to the real medium. In the case of 3D observations, a similar criterion is constructed on the basis of the R-method and analysis of the surface t о (x, у).
3. Выполняют миграцию временного разреза в пластовой модели среды. 3. Migration of a temporary section in the reservoir model of the environment.
Миграция выполняется с учетом преломления лучей в модели, адекватность которой реальной среде установлена методом, описанном в предыдущем пункте. Это позволяет получить правильное изображение глубинного строения разреза. Migration is performed taking into account the refraction of rays in a model, the adequacy of which to a real medium is established by the method described in the previous paragraph. This allows you to get the correct image of the deep structure of the section.
Миграция выполняется послойно, и в каждом слое имеется возможность управлять основным параметром, определяющим качество преобразования - апертурой. В частности, поскольку глубинно-скоростная модель к этому моменту известна достаточно точно, вычисляются точки выхода нормальных лучей от всех отражающих границ на дневную поверхность (точки накопления сигналов миграционным оператором) и только в их окрестности выполняется преобразование сигналов. Миграционное преобразование выполняется с учетом функции запаздывания синфазности относительно годографа дифрагированной волны (рассчитанного в слоистой среде), что улучшает форму и уточняет амплитуду преобразованного сигнала. В результате в сложных геологических разрезах (разломы, выклинивания. петли на временном разрезе) получается достаточно точное и высокого качества изображение среды. Migration is performed in layers, and in each layer it is possible to control the main parameter that determines the quality of the transformation - the aperture. In particular, since the depth-velocity model is known quite accurately by this moment, the points of exit of normal rays from all reflecting boundaries to the day surface (points of accumulation of signals by the migration operator) are calculated and only in their vicinity the signals are converted. Migration transformation is performed taking into account the delay function of the in-phase with respect to the travel time curve of the diffracted wave (calculated in a layered medium), which improves the shape and refines the amplitude of the converted signal. As a result, in complex geological sections (faults, pinch-outs, loops on a temporary section), a fairly accurate and high-quality image of the medium is obtained.
4. Производят погружение сейсмограмм на заданный уровень вглубь среды. 4. Submerge seismograms to a predetermined level deep into the medium.
Эта операция имеет ту же цель, что и миграция сейсмограмм, а именно: получение правильного изображения среды в сложных разрезах, когда годографы ОГТ существенно отличаются от гипербол. В отличие от широко используемой миграции до суммирования, предлагаемая операция погружения характеризуется следующими достоинствами:
а) для ее выполнения используется только верхняя часть разреза, где годографы достаточно гиперболичны. Глубинно-скоростная модель строится довольно просто, а ее адекватность реальной среде проверяется так, как описано выше;
б) оператор погружения стационарен по времени, что приводит к существенному уменьшению искажения сигналов, чем при миграции сейсмограмм и значительно ускоряет вычисления.This operation has the same goal as the migration of seismograms, namely: obtaining the correct image of the medium in complex sections, when the hodographs of the OGT are significantly different from hyperbolas. Unlike the widely used pre-stack migration, the proposed dive operation has the following advantages:
a) for its implementation, only the upper part of the section is used, where the hodographs are sufficiently hyperbolic. The deep-speed model is built quite simply, and its adequacy to the real environment is checked as described above;
b) the immersion operator is stationary in time, which leads to a significant reduction in signal distortion than with the migration of seismograms and significantly speeds up the calculation.
Погруженные сейсмограммы обрабатываются также как исходные (наблюдаем на дневной поверхности), включая такие операции, как статика, различные фильтрации и т. п. До некоторой глубины от уровня приведения годографы сраженных волн (которые были негиперболичны на дневной поверхности) значительно лучше описываются гиперболами. Поэтому операциями, которые описаны выше, снова можно надежно определить их кинематические параметры, построить глубинно-скоростную модель среды, адекватную конкретно обрабатываемому разрезу и, при необходимости, повторить операцию погружения на следующий уровень. Immersed seismograms are also processed as the initial ones (we observe on the surface), including operations such as static, various filtering, etc. Up to a certain depth from the reduction level, the travel time curves of the beaten waves (which were non-hyperbolic on the surface) are much better described by hyperbolas. Therefore, with the operations described above, it is again possible to reliably determine their kinematic parameters, build a deep-speed model of the medium adequate to the particular section being processed, and, if necessary, repeat the immersion operation to the next level.
В результате в сложных сейсмогеологических условиях получается полная глубинно-скоростная модель среды и динамическое изображение разреза значительно более высокого качества, чем при других применяемых способах. As a result, under complex seismic and geological conditions, a complete depth-velocity model of the medium and a dynamic image of the section are obtained of much higher quality than with other methods used.
5. Определяют поле кратных волн с последующим их адаптивным вычитанием. 5. Determine the field of multiple waves with their subsequent adaptive subtraction.
Для моделирования поля кратных волн предварительно строится глубинно-скоростная модель той части среды, которая содержит кратно-образующие границы. Она строится способами, описанными в предыдущих пунктах (включая проверку ее адекватности реальной среде). Затем исходные сейсмограммы преобразуются оператором миграционного типа, отличающимся от обычно применяемых тем, что в качестве преобразующего годографа используется годограф отраженной, а не дифрагированной волны. В результате этого преобразования каждая существующая кратная волна повышает свою кратность: однократная становится трехкратной, трехкратная пятикратной и т.д. Таким образом, преобразованные сейсмограммы содержат только кратные волны. Для того, чтобы вычесть эту совокупность кратных волн из исходной записи, требуется еще привести интенсивность и форму рассчитанных кратных волн к их интенсивности и форме на исходной записи (т. к. моделирование кратных выполняется без учета коэффициентов отражения на кратно-образующих границах). Это делается посредством расчета и применения многомерных фильтров, приближающих цуг смоделированных кратных к заданной сейсмической трассе, причем для адаптации используются смоделированные кратные на нескольких соседних пространственных координатах. Число параметров описывающих флуктуации формы кратных волн по пространственным координатам и их значения определяются адаптивно но изменению среднеквадратического критерия оптимизации в заданном диапазоне частот. To model the field of multiple waves, a depth-velocity model of the part of the medium that contains the multiple-forming boundaries is preliminarily constructed. It is built by the methods described in the previous paragraphs (including checking its adequacy to the real environment). Then, the initial seismograms are transformed by a migration type operator, which differs from those usually used in that the hodograph of the reflected rather than diffracted wave is used as the transforming hodograph. As a result of this transformation, each existing multiple wave increases its frequency: a single one becomes threefold, threefold fivefold, etc. Thus, transformed seismograms contain only multiple waves. In order to subtract this set of multiple waves from the original record, it is still necessary to bring the intensity and shape of the calculated multiple waves to their intensity and shape on the original record (since the simulation of multiple waves is performed without taking into account the reflection coefficients at the multiple-forming boundaries). This is done by calculating and applying multidimensional filters that approximate the train of simulated multiples to a given seismic path, and for adaptation, simulated multiples at several adjacent spatial coordinates are used. The number of parameters describing fluctuations of the shape of multiple waves in spatial coordinates and their values are determined adaptively but to a change in the mean-square optimization criterion in a given frequency range.
Операции моделирования и адаптивного вычитания кратных волн разделены между собой. Это отличие от методов, не разделяющих этаны моделирования и вычитания, позволяет точнее построить модель поля кратных и повысить эффективность вычитания этих помех, снизить опасность подавления полезных отражений. The operations of modeling and adaptive subtraction of multiple waves are divided among themselves. This difference from methods that do not separate the ethanes of modeling and subtraction allows us to more accurately construct a multiple field model and increase the efficiency of subtracting these noise, and reduce the risk of suppressing useful reflections.
6. Осуществляют амплитудную и фазовую деконволюцию, согласованную фильтрацию и коррекцию поглощения. 6. Carry out amplitude and phase deconvolution, consistent filtering and absorption correction.
Способы амплитудной и фазовой деконволюции, согласованной фильтрации и расчета коррекции поглощения основаны на оценках параметров, описывающих амплитудные и фазовые спектры сигналов и эффект частотно-зависимого поглощения. Амплитудный и фазовый спектры описываются отрезками ряда Фурье в заданном диапазоне частот. Поглощение описывается одним параметром, показателем экспоненты на которую домножается амплитудный спектр сигнала. Параметры амплитудного спектра сигнала и параметры поглощения оцениваются одновременно но нескольким пространственно-временным интервалам записей. При этом амплитудный спектр сигнала считается неизменным для всех окон, а параметры поглощения различны. Таким образом, для оценивания сравнительно небольшого количества параметров, привлекаются значительные объемы выборок, что позволяет получать оценки амплитудных спектров и параметров поглощения с высокой точностью и использовать их для различных корректирующих фильтраций. В случае оценивания параметров амплитудного спектра и поглощения используется критерий минимума дисперсии ошибки и предсказания. При оценивании параметров фазового спектра используются неквадратичные критерии максимизируются статистические моменты, в том числе и с нецелыми степенями (больше двух), либо минимизируются статистические моменты со степенями меньше двух и больше единицы. Во всех методах для оценивания параметров используются способы условной оптимизации (в случае фазовых задач учитывается многоэкстремальность функционалов). Ограничения на значения параметров позволяют сузить область их поиска и снизить возможность искажений сигналов на этапе их коррекции. The methods of amplitude and phase deconvolution, coordinated filtering and calculation of absorption correction are based on parameter estimates describing the amplitude and phase spectra of signals and the effect of frequency-dependent absorption. The amplitude and phase spectra are described by segments of the Fourier series in a given frequency range. The absorption is described by one parameter, the exponent of which is multiplied by the amplitude spectrum of the signal. The parameters of the amplitude spectrum of the signal and the absorption parameters are evaluated simultaneously for several spatio-temporal recording intervals. In this case, the amplitude spectrum of the signal is considered unchanged for all windows, and the absorption parameters are different. Thus, to evaluate a relatively small number of parameters, significant sample sizes are involved, which allows one to obtain estimates of the amplitude spectra and absorption parameters with high accuracy and use them for various corrective filtrations. In the case of estimating the parameters of the amplitude spectrum and absorption, the criterion of the minimum variance of error and prediction is used. When estimating the parameters of the phase spectrum, non-quadratic criteria are used to maximize statistical moments, including those with non-integer degrees (more than two), or to minimize statistical moments with degrees less than two and more than one. In all methods, conditional optimization methods are used to estimate parameters (in the case of phase problems, the multi-extremality of functionals is taken into account). Restrictions on the values of the parameters make it possible to narrow the scope of their search and reduce the possibility of signal distortion at the stage of their correction.
Предлагаемые способы оценивания параметров, описывающих форму сигнала и ее изменения, позволяют повысить точность их оценивания, а значит добиться более высокой разрешенности сейсмических записей, улучшить отношение сигнал\ помеха и, в целом, повысить достоверность построения сейсмических изображений. The proposed methods for estimating the parameters describing the waveform and its changes make it possible to increase the accuracy of their estimation, which means to achieve a higher resolution of seismic records, improve the signal-to-noise ratio and, in general, increase the reliability of the construction of seismic images.
Для доказательства соответствия критерию охраноспособности "промышленная применимость" приводится конкретный пример реализации способа. To prove compliance with the eligibility criterion "industrial applicability", a specific example of the implementation of the method is given.
На чертеже изображен временной разрез (Российское Заполярье). The drawing shows a time section (Russian Arctic).
Для получения временного разреза были использованы все перечисленные отличительные признаки предложенного способа. Окончательный глубинный мигрированный разрез (показанный па рисунке) демонстрирует, что основной объект разведки в этом регионе - неокомские отложения (выше Баженовской свиты на глубине 4000 м) - проявились во всех деталях. Клиноформы, которые видны на пикетах ~ 6000 м и 20000 м часто связаны с залежами нефти или газа. To obtain a temporary section, all of the above distinguishing features of the proposed method were used. The final deep migrated section (shown in the figure) demonstrates that the main exploration object in this region - the Neocomian deposits (above the Bazhenov Formation at a depth of 4000 m) - appeared in all details. Clinoforms that are visible at pickets of ~ 6,000 m and 20,000 m are often associated with oil or gas deposits.
Совокупность всех перечисленных признаков предложенного способа позволяет достичь более точного определения параметров реальной среды и высокого качества ее динамического глубинного изображения в сложных сейсмогеологических условиях, в том числе и для сейсмических данных невысокого качества, которое часто встречается на практике. The combination of all the listed features of the proposed method allows to achieve a more accurate determination of the parameters of the real environment and the high quality of its dynamic deep image in complex seismic and geological conditions, including low-quality seismic data, which is often found in practice.
Источники информации, принятые во внимание:
1. J.R. Berryhill, Y.K. Kim Deep-water peg-leg and mulliples: Emulation and Suppression. Geophysics, 1986, v.51, N 12, pp. 2177-2184.Sources of information taken into account:
1. JR Berryhill, YK Kim Deep-water peg-leg and mulliples: Emulation and Suppression. Geophysics, 1986, v. 51, N 12, pp. 2177-2184.
2. Система Geodepth. фирма Paradigm Geophysical, журнал "Нефтегаз - 98", N 1, стр. 175 (прототип). 2. Geodepth system. Paradigm Geophysical company, journal "Neftegaz - 98", N 1, p. 175 (prototype).
Claims (1)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU98113007A RU2126984C1 (en) | 1998-07-16 | 1998-07-16 | Method determining depth and speed parameters of medium and construction of its image by seismic data- prime system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU98113007A RU2126984C1 (en) | 1998-07-16 | 1998-07-16 | Method determining depth and speed parameters of medium and construction of its image by seismic data- prime system |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2126984C1 true RU2126984C1 (en) | 1999-02-27 |
RU98113007A RU98113007A (en) | 1999-04-20 |
Family
ID=20208125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU98113007A RU2126984C1 (en) | 1998-07-16 | 1998-07-16 | Method determining depth and speed parameters of medium and construction of its image by seismic data- prime system |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2126984C1 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2526794C2 (en) * | 2012-08-29 | 2014-08-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Российский государственный университет нефти и газа имени И.М.Губкина" (РГУ нефти и газа имени И.М.Губкина) | Method of determining elastic properties of rocks based on stratal adaptive inversion of seismic data |
RU2577792C1 (en) * | 2014-09-22 | 2016-03-20 | Общество с ограниченной ответственностью "ГЕОЛАБ-ИТ", ООО "ГЕОЛАБ-ИТ" | Robust process for depth imaging in seismic survey based on operator adjustment by reference seismic record |
RU2601232C2 (en) * | 2012-08-20 | 2016-10-27 | Лэндмарк Графикс Корпорейшн | Methods and systems for connection pseudosurface selected point places into speed of seismic waves models |
RU2670358C1 (en) * | 2018-01-31 | 2018-10-22 | Керим Акмамедович Хемраев | Method of dynamic correction of reflected waves with use of refracted waves |
RU2705519C2 (en) * | 2018-02-14 | 2019-11-07 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Уральский государственный горный университет" (ФГБОУ ВО "УГГУ") | Method of producing migrated seismic images of geologic media based on 2d seismic survey data |
RU2709415C1 (en) * | 2019-04-02 | 2019-12-17 | Акционерное общество "ВНИИ Галургии" (АО "ВНИИ Галургии") | Method of shaft seismic survey |
RU2710972C1 (en) * | 2019-10-28 | 2020-01-14 | Общество с ограниченной ответственностью «Сейсмотек» | Multivariate tomography method of seismic survey data |
-
1998
- 1998-07-16 RU RU98113007A patent/RU2126984C1/en not_active IP Right Cessation
Non-Patent Citations (1)
Title |
---|
Система Geodepth, фирма Paradigm Geophysics - Нефтегаз-98, N 1, с.175. * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2601232C2 (en) * | 2012-08-20 | 2016-10-27 | Лэндмарк Графикс Корпорейшн | Methods and systems for connection pseudosurface selected point places into speed of seismic waves models |
RU2526794C2 (en) * | 2012-08-29 | 2014-08-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Российский государственный университет нефти и газа имени И.М.Губкина" (РГУ нефти и газа имени И.М.Губкина) | Method of determining elastic properties of rocks based on stratal adaptive inversion of seismic data |
RU2577792C1 (en) * | 2014-09-22 | 2016-03-20 | Общество с ограниченной ответственностью "ГЕОЛАБ-ИТ", ООО "ГЕОЛАБ-ИТ" | Robust process for depth imaging in seismic survey based on operator adjustment by reference seismic record |
WO2016048194A1 (en) * | 2014-09-22 | 2016-03-31 | Общество С Ограниченной Ответственностью "Геолаб-Ит" | Method for constructing depth images on the basis of reference seismograms |
RU2670358C1 (en) * | 2018-01-31 | 2018-10-22 | Керим Акмамедович Хемраев | Method of dynamic correction of reflected waves with use of refracted waves |
RU2705519C2 (en) * | 2018-02-14 | 2019-11-07 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Уральский государственный горный университет" (ФГБОУ ВО "УГГУ") | Method of producing migrated seismic images of geologic media based on 2d seismic survey data |
RU2709415C1 (en) * | 2019-04-02 | 2019-12-17 | Акционерное общество "ВНИИ Галургии" (АО "ВНИИ Галургии") | Method of shaft seismic survey |
RU2710972C1 (en) * | 2019-10-28 | 2020-01-14 | Общество с ограниченной ответственностью «Сейсмотек» | Multivariate tomography method of seismic survey data |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10976459B2 (en) | Final statics calculation for automated near surface analysis | |
US10429530B2 (en) | Deghosting with adaptive operators | |
US8086426B2 (en) | Processing seismic data representing a physical system | |
EP0889331B1 (en) | Method for 3D modelling of the impedance of a heterogeneous volume of the subsurface | |
US7477992B2 (en) | Method for combining seismic data sets | |
US7082368B2 (en) | Seismic event correlation and Vp-Vs estimation | |
US10670750B2 (en) | Multistage full wavefield inversion process that generates a multiple free data set | |
US20130003500A1 (en) | Seismic Data Processing | |
US7616524B1 (en) | Wavelet based intercept attribute for seismic exploration | |
US9348050B2 (en) | Near-surface noise prediction and removal for data recorded with simultaneous seismic sources | |
Wang et al. | Separation and imaging of seismic diffractions using a localized rank-reduction method with adaptively selected ranks | |
CN105474048B (en) | The interbed multiple in prediction seismic data is decomposed using wave beam | |
WO2005026776A1 (en) | Wide-offset-range pre-stack depth migration method for seismic exploration | |
WO1997011390A2 (en) | Multiple suppression in geophysical data | |
RU2126984C1 (en) | Method determining depth and speed parameters of medium and construction of its image by seismic data- prime system | |
WO2002046792A2 (en) | A method of and apparatus for processing seismic data | |
AU2002234794A1 (en) | A method of and apparatus for processing seismic data | |
CN110988991B (en) | Elastic parameter inversion method, device and system | |
Li et al. | An efficient deep learning method for VSP wavefield separation: A DAS-VSP case | |
Zhou et al. | Surface diffraction noise attenuation for marine seismic data processing with mathematical morphological filtering | |
Soleimani et al. | Merging aspects of DMO correction and CRS stack to account for conflicting dip situations | |
Zhang et al. | Generalized internal multiple prediction for low relief structures | |
Zdraveva et al. | Improving interpretability of seismic images with directional image partitions and model-based techniques | |
CN113495291A (en) | Method for intelligently and quantitatively evaluating amplitude preservation of forward gather based on deep learning | |
CN117805907A (en) | Earthquake phase driven hydrocarbon source rock prediction method, device and equipment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PC4A | Invention patent assignment |
Effective date: 20080131 |
|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20080717 |
|
NF4A | Reinstatement of patent |
Effective date: 20100227 |
|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20160717 |