RU2722861C1 - Static corrections calculation method - Google Patents
Static corrections calculation method Download PDFInfo
- Publication number
- RU2722861C1 RU2722861C1 RU2019121747A RU2019121747A RU2722861C1 RU 2722861 C1 RU2722861 C1 RU 2722861C1 RU 2019121747 A RU2019121747 A RU 2019121747A RU 2019121747 A RU2019121747 A RU 2019121747A RU 2722861 C1 RU2722861 C1 RU 2722861C1
- Authority
- RU
- Russia
- Prior art keywords
- data
- field
- section
- seismic
- formation
- Prior art date
Links
- 230000003068 static effect Effects 0.000 title claims abstract description 30
- 238000012937 correction Methods 0.000 title claims abstract description 27
- 238000004364 calculation method Methods 0.000 title description 3
- 238000000034 method Methods 0.000 claims abstract description 42
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 12
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 5
- 230000005284 excitation Effects 0.000 claims abstract description 3
- 239000011435 rock Substances 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 3
- 230000008859 change Effects 0.000 abstract description 2
- 239000000126 substance Substances 0.000 abstract 1
- 238000005553 drilling Methods 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
- G01V11/007—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00 using the seismo-electric effect
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/081—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Acoustics & Sound (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Изобретение относится к комплексу методов геофизической разведки, включающему сейсморазведку методом отраженных волн общей глубинной точки (MOB ОГТ) и электроразведку методом малоглубинных зондирований становлением поля в ближней зоне (мЗСБ), и может быть использовано для учета скоростных аномалий верхней части разреза (ВЧР).The invention relates to a set of methods for geophysical exploration, including seismic exploration by the method of reflected waves of the common deep point (MOB OGT) and electrical exploration by the method of shallow sounding by the formation of a field in the near zone (MSSB), and can be used to account for high-speed anomalies of the upper part of the section (VChR).
При обработке и интерпретации материалов сейсморазведочных работ на территориях Сибири и Арктики, где развиты многолетнемерзлые породы (ММП), необходимо учитывать их влияние на материалы сейсморазведки. Источником такого влияния выступают скоростные аномалии, сосредоточенные в относительно тонком, но неоднородном по толщине приповерхностном интервале залегания слоев ММП. Недоучет влияния ММП на форму отражающих горизонтов может привести к значительным ошибкам в структурных построениях и ухудшению волновой картины во всем временном диапазоне.When processing and interpreting seismic data in Siberia and the Arctic, where permafrost is developed, it is necessary to take into account their influence on seismic data. The source of such an effect is velocity anomalies concentrated in a relatively thin but nonuniform in thickness near-surface interval of the occurrence of IMF layers. Underestimation of the influence of the IMF on the shape of reflecting horizons can lead to significant errors in the structural structures and the deterioration of the wave pattern in the entire time range.
Известен способ скважинной сейсморазведки (см. патент РФ на изобретение №2292063, МПК G01V 1/40), в котором определяют скоростные характеристики верхней части разреза - зоны малых скоростей (ЗМС). Сущность способа: в процессе бурения возбуждают упругие колебания путем воздействия породоразрушающего инструмента мобильной буровой установки на исследуемую среду. Одновременно регистрируют упругие колебания датчиком опорного сигнала, состоящим из четырех преобразователей, равномерно распределенных по периметру площадки рамы буровой установки, и наземным приемным устройством. Наземное приемное устройство устанавливают на дневной поверхности на расстоянии от устья скважины не менее 5-10 м на стержнях, заглубленных в грунт на глубину, превышающую мощность почвенного слоя. Выбирают диапазон рабочих частот от 100 Гц до 350 Гц, в пределах которого выделяют полезные сигналы. Формируют взаимокорреляционные функции и определяют по ним сейсмические скорости и положение сейсмических границ. Технический результат: повышение точности и достоверности построения скоростной характеристики исследуемой среды. Для проведения данного способа необходимо наличие скважины глубокого бурения.A known method of downhole seismic exploration (see RF patent for the invention No. 2292063, IPC G01V 1/40), which determines the speed characteristics of the upper part of the section - the zone of low speeds (ZMS). The essence of the method: in the process of drilling, elastic vibrations are excited by the influence of the rock cutting tool of the mobile drilling rig on the medium under study. At the same time, elastic vibrations are recorded by a reference signal sensor consisting of four transducers uniformly distributed along the perimeter of the platform frame of the drilling rig and a ground receiving device. The ground receiving device is installed on the day surface at a distance from the wellhead of at least 5-10 m on the rods buried in the ground to a depth exceeding the thickness of the soil layer. A range of operating frequencies from 100 Hz to 350 Hz is selected, within which useful signals are extracted. Intercorrelation functions are formed and seismic velocities and the position of seismic boundaries are determined from them. Effect: increase the accuracy and reliability of the construction of the speed characteristics of the investigated environment. To carry out this method, it is necessary to have a deep hole.
Известны способы определения статических поправок методом регистрации волны, преломленной на подошве зоны малых скоростей (см. патент Великобритании №32090405, МПК G01V 1/28 и Шариф Р., Гелгарт Л., Сейсморазведка, т. 1, М.: Мир, 1987). Однако данный способ дает усредненное значение физических характеристик ВЧР до подошвы зоны малых скоростей и не учитывает наличие локальных неоднородностей в ней, а это крайне важно в случае мощной зоны ММП (400-600 м).Known methods for determining static corrections by registering a wave refracted on the sole of a low-speed zone (see UK patent No. 32090405, IPC G01V 1/28 and Sharif R., Gelgart L., Seismic exploration, t. 1, M .: Mir, 1987) . However, this method gives an average value of the physical characteristics of the VChR to the bottom of the low-velocity zone and does not take into account the presence of local inhomogeneities in it, and this is extremely important in the case of a thick IMF zone (400-600 m).
Наиболее близким к предлагаемому решению является способ определения статических поправок (см. патент РФ на изобретение №2411547, МПК G01V 1/36, G01V 11/00), который принят нами за прототип, в котором прогнозирование статических поправок включает проведение сейсморазведки, обработку и интерпретацию полученных данных, отличающийся тем, что дополнительно на совмещенных профилях проводят электроразведку для изучения строения верхней части разреза в зоне малых скоростей, регистрируют изменения электромагнитного поля и годографы электромагнитной волны, определяют продольное электрическое сопротивление, по полученным данным выделяют единые стратиграфически увязанные геоэлектрические комплексы, строят геоэлектрическую модель верхней части разреза, затем с использованием данных геофизических исследований скважин или микросейсмического каротажа устанавливают взаимосвязь между временами регистрации электромагнитного и сейсмического полей, в каждой точке электроразведочных наблюдений пересчитывают электромагнитные годографы в псевдосейсмические, по ним в пределах каждого геоэлектрического комплекса рассчитывают значения прогнозных интервальных скоростей, для выбранного интервала разреза строят схемы распределения прогнозных значений интервальных скоростей и его толщины, рассчитывают значения статических поправок. Данный способ имеет недостаток, связанный с тем, что представленный способ не дает достаточной глубинности изучения верхней части разреза, а охватывает лишь ЗМС (глубины от дневной поверхности до 100 м). Также ввиду того, что сигнал становления электромагнитного поля является гладкой экспоненциальной функцией, использование кривых мЗСБ для расчета скоростных моделей приводит к тому, что резкие изменения разреза (аномалии ВЧР) не детектируются.Closest to the proposed solution is a method for determining static corrections (see RF patent for invention No. 2411547, IPC G01V 1/36, G01V 11/00), which we adopted as a prototype, in which forecasting static corrections includes seismic surveying, processing and interpretation obtained data, characterized in that in addition to the combined profiles, electrical exploration is carried out to study the structure of the upper part of the section in the low-velocity zone, changes in the electromagnetic field and hodographs of the electromagnetic wave are recorded, longitudinal electrical resistance is determined, according to the obtained data, single stratigraphically linked geoelectric complexes are distinguished, and the geoelectric is built model of the upper part of the section, then using the data of geophysical studies of wells or microseismic logging, the relationship between the times of registration of electromagnetic and seismic fields is established, at each point of electrical exploration observations I recount t electromagnetic hodographs into pseudo-seismic ones, according to them, within each geoelectric complex, the values of the predicted interval velocities are calculated, for the selected interval of the section, the distribution patterns of the predicted values of the interval velocities and its thickness are constructed, and the values of the static corrections are calculated. This method has the disadvantage that the presented method does not provide sufficient depth of study of the upper part of the section, but covers only the ZMS (depths from the surface to 100 m). Also, due to the fact that the electromagnetic field formation signal is a smooth exponential function, the use of the MZSB curves for calculating velocity models leads to the fact that sharp changes in the section (VLF anomalies) are not detected.
Задачей заявленного авторами способа является разработка эффективного подхода к расчету скоростной модели ВЧР на основе данных зондирования становлением поля в ближней зоне в малоглубинной модификации (мЗСБ) с глубинностью от дневной поверхности до 400-600 м.The objective of the invented method is to develop an effective approach to calculating the velocity model of the VChR based on sounding data by the formation of a field in the near field in a shallow modification (mSZB) with a depth from the surface to 400-600 m.
Технический результат заявляемого решения заключается в повышении точности, информативности способа определения статических поправок за счет учета локальных неоднородностей строения ВЧР в мощной зоне ММП до 600 м.The technical result of the proposed solution is to increase the accuracy, information content of the method for determining static corrections by taking into account local heterogeneities of the structure of the VChR in the powerful zone of the permafrost zone up to 600 m.
Поставленная задача решается за счет того, что в способе определения статических поправок, включающем в себя возбуждение сейсмического сигнала, регистрацию волновой картины сейсмического поля, обработку и интерпретацию полученных данных, проведение дополнительно на совмещенных профилях электроразведки для изучения строения верхней части разреза в зоне малых скоростей, регистрируют изменение электромагнитного поля, при этом электроразведку выполняют методом малоглубинного зондирования становлением поля в ближней зоне по высокоплотным сетям наблюдений, после обработки данных малоглубинного зондирования становлением поля в ближней зоне выполняют инверсию полученных данных, по результатам инверсии определяют продольное электрическое сопротивление, затем по полученным данным выделяют единые стратиграфически увязанные геоэлектрические комплексы, проводят структурную интерпретацию верхней части разреза и строят геоэлектрическую модель верхней части разреза, затем используя эмпирическую зависимость Фауста, коэффициенты которой определяются на основе данных геофизических исследований скважин или вертикального сейсмического профилирования, геоэлектрические модели, полученные в результате инверсии данных малоглубинного зондирования становлением поля в ближней зоне, пересчитывают в акустические модели, по полученным акустическим моделям строят куб скоростей верхней части разреза, по которому рассчитывают статические поправки для последующего использования при обработке данных сейсморазведки метода отраженных волн общей глубинной точки.The problem is solved due to the fact that in the method for determining static corrections, which includes the excitation of a seismic signal, registration of the wave pattern of the seismic field, processing and interpretation of the obtained data, additionally conducted on combined profiles of electrical exploration to study the structure of the upper part of the section in the low-velocity zone, a change in the electromagnetic field is recorded, while electrical exploration is performed by the method of shallow sounding by establishing a field in the near zone using high-density observation networks, after processing the data of shallow sounding by establishing a field in the near zone, the obtained data are inverted, the longitudinal electrical resistance is determined by the results of the inversion, then unified stratigraphically linked geoelectric complexes, conduct structural interpretation of the upper part of the section and construct a geoelectric model of the upper part of the section, then using empirical dependencies Faust bridge, the coefficients of which are determined on the basis of data from geophysical surveys of wells or vertical seismic profiling, geoelectric models obtained as a result of the inversion of shallow sounding data by the formation of a field in the near zone are converted into acoustic models, the cube of velocities of the upper section is constructed from the obtained acoustic models, using which is calculated static corrections for subsequent use in the processing of seismic data from the method of reflected waves at a common depth point.
Изобретение иллюстрируется следующими чертежами:The invention is illustrated by the following drawings:
Фиг. 1. Граф расчета статических поправок.FIG. 1. Graph for calculating static corrections.
Фиг. 2. График анализа ошибки определения скорости продольных волн Vp, иллюстрирующий подбор коэффициентов уравнения Фауста с использованием данных ВСП.FIG. 2. Graph of the analysis of the error in determining the velocity of longitudinal waves Vp, illustrating the selection of the coefficients of the Faust equation using VSP data.
Фиг. 3. График сходимости скоростей продольных волн, иллюстрирующий подбор коэффициентов уравнения Фауста с использованием данных ВСП.FIG. 3. A graph of convergence of longitudinal wave velocities illustrating the selection of the coefficients of the Faust equation using VSP data.
Фиг. 4. Результат ввода статических поправок в сейсмограммы ОГТ, при этом статические поправки вводились в двух вариантах: вариант, рассчитанный по методу первых вступлений (традиционная методика), график «Статические поправки по первым вступлениям» и вариант, полученный на основе скоростной модели по данным мЗСБ (заявляемый способ), график «Статические поправки по данным мЗСБ».FIG. 4. The result of introducing static corrections into the CDP seismograms, while the static corrections were introduced in two versions: the variant calculated by the first arrivals method (traditional method), the graph “Static corrections by first arrivals” and the variant obtained on the basis of the velocity model according to the MSSS (the claimed method), the schedule "Static amendments according to the Ministry of Health and Social Security".
Фиг. 5. Временной сейсмический разрез, иллюстрирующий в сравнении вариант расчета статических поправок по методу первых вступлений (традиционная методика), график «Временной сейсмический разрез (статика по первым вступлениям)» и вариант согласно заявляемому способу, график «Временной сейсмический разрез (статика по данным мЗСБ)».FIG. 5. Temporary seismic section, illustrating in comparison the option of calculating static corrections according to the method of first arrivals (traditional methodology), the graph “Temporary seismic section (statics of the first arrivals)” and the variant according to the claimed method, graph “Temporary seismic section (statics according to the data of the Ministry of Safety and Social Protection) ) ".
На основе полученных геоэлектрических моделей строят геолого-геоэлектрическую модель верхней части разреза, проводят структурную интерпретацию.Based on the obtained geoelectric models, a geological-geoelectric model of the upper part of the section is built, and a structural interpretation is carried out.
Следующим шагом необходимо осуществить переход от удельного электрического сопротивления разреза к скорости продольных волн. Для того чтобы от геоэлектрических свойств пород перейти к акустическим, возможно использование эмпирических зависимостей. Зависимость между удельным электрическим сопротивлением и скоростью продольных волн впервые была представлена Л. Фаустом в 1951 г. (фиг. 1 Граф расчета статических поправок). Актуальность работы была обусловлена необходимостью расчета карт средних скоростей продольных сейсмических волн для решения задачи структурных построений.The next step is to make the transition from the electrical resistivity of the section to the velocity of longitudinal waves. In order to move from the geoelectric properties of rocks to acoustic, it is possible to use empirical relationships. The relationship between the electrical resistivity and the velocity of longitudinal waves was first presented by L. Faust in 1951 (Fig. 1 Graph for calculating static corrections). The relevance of the work was due to the need to calculate the average velocity maps of longitudinal seismic waves to solve the problem of structural constructions.
В результате подбора параметров модели была получена следующая эмпирическая зависимость:As a result of the selection of model parameters, the following empirical dependence was obtained:
где α - константа, Z - глубина залегания, Т - возраст пород изучаемого разреза. По данным, приведенным в работе, погрешность вычислений скорости продольных сейсмических волн по приведенной выше зависимости составила 2-5%.where α is a constant, Z is the depth, T is the age of the rocks of the studied section. According to the data presented in the work, the error in calculating the velocity of longitudinal seismic waves from the above dependence was 2-5%.
В 1953 г. Л. Фауст описывает результаты вычислений скорости уже с учетом влияния литологической составляющей [Faust, 1953]. Для учета литологии использовались данные электрического каротажа. В ходе экспериментов сопоставлялись значения скорости продольных волн, полученных по данным сейсмического каротажа, и значения удельного сопротивления горной породы. В результате анализа была получена формула:In 1953, L. Faust describes the results of velocity calculations taking into account the influence of the lithological component [Faust, 1953]. To account for lithology, electrical logging data was used. During the experiments, the values of the longitudinal wave velocity obtained from seismic logging data and the resistivity of the rock were compared. As a result of the analysis, the formula was obtained:
где R - удельное сопротивление пласта (Ом⋅м). Оценка погрешности прогноза скорости по данному уравнению составила 1-2.5%.where R is the formation resistivity (Ohm⋅m). The estimation of the error of the forecast of speed according to this equation was 1-2.5%.
В упрощенном виде зависимость Фауста можно представить, как:In a simplified form, the dependence of Faust can be represented as:
где sonic - интервальное время пробега продольной волны, TVD - глубина в футах, resistivity - удельное электрическое сопротивление, const и ехр - эмпирические коэффициенты уравнения.where sonic is the longitudinal wave travel time, TVD is the depth in feet, resistivity is the electrical resistivity, const and exp are the empirical coefficients of the equation.
Далее выполняется расчет и калибровка эмпирических коэффициентов уравнения Фауста. Д ля расчета коэффициентов необходимо наличие акустического каротажа или данных вертикального сейсмического профилирования (ВСП). Первым шагом формируется литологически-генерализованная модель, в рамках каждой литологической разности производится многократный перебор эмпирических коэффициентов уравнения (1). Результатом является уникальная пара коэффициентов для каждой литологической разности. При достижении суммарного коэффициента корреляции r=0.9 и более коэффициенты считаются подобранными (на графике на фиг. 2. показан пример подбора коэффициентов зависимости Фауста).Next, the calculation and calibration of the empirical coefficients of the Faust equation are performed. To calculate the coefficients, the presence of acoustic logging or vertical seismic profiling (VSP) is necessary. The first step is the formation of a lithologically generalized model, and within each lithological difference, multiple empirical coefficients of equation (1) are enumerated. The result is a unique pair of coefficients for each lithological difference. Upon reaching the total correlation coefficient r = 0.9 or more, the coefficients are considered matched (the graph in Fig. 2. shows an example of the selection of the Faust dependence coefficients).
Далее выполняется преобразование геоэлектрических моделей* полученных по данным мЗСБ, в акустические модели.Next, the geoelectric models * obtained from the Ministry of Health and Social Protection are transformed into acoustic models.
Благодаря тому, что метод мЗСБ является надежным инструментом для изучения структурных особенностей разреза, с его помощью возможно достоверно восстанавливать структуру ММП. Владея информацией о распределении скоростей ВЧР и достоверной структурной моделью, возможно с высокой степенью достоверности прогнозировать статические поправки.Due to the fact that the MZSB method is a reliable tool for studying the structural features of the section, it is possible to reliably restore the structure of the permafrost layer with its help. Possessing information on the distribution of HF velocity and a reliable structural model, it is possible to predict static corrections with a high degree of certainty.
Существенными отличиями заявляемого способа в сравнении с известными техническими решениями являются:Significant differences of the proposed method in comparison with the known technical solutions are:
- дополнительно с сейсморазведкой на совмещенных профилях проводят электроразведку в мЗСБ с глубинностью до 400-600 м, но по высокоплотным сетям наблюдений, комплексно интерпретируют полученные данные, строят геоэлектрическую модель ВЧР;- in addition to seismic exploration on combined profiles, electrical prospecting is carried out in the MZSB with a depth of up to 400-600 m, but using high-density observation networks, they comprehensively interpret the data obtained and construct a geoelectric model of the VChR;
- с использованием данных геофизических исследований скважин (ТИС) или вертикального сейсмического профилирования (ВСП) производят расчет эмпирических коэффициентов уравнения Фауста;- using the data of geophysical research of wells (TIS) or vertical seismic profiling (VSP), the empirical coefficients of the Faust equation are calculated;
- для расчета статических поправок используется комплексная геоэлектрическая и структурная модель ВЧР в совокупности с распределением скоростей ВЧР, полученных по данным мЗСБ.- to calculate the static corrections, a complex geoelectric and structural model of the VChR is used in conjunction with the distribution of the velocities of the VChR obtained according to the data of the MSSS.
Способ реализуется следующим образом:The method is implemented as follows:
В пределах изучаемой площади одновременно (либо последовательно) с сейсморазведкой проводят электроразведку в малоглубинной модификации зондированием становлением поля (мЗСБ) по линиям сейсмических профилей, с высокой пространственно-временной плотностью.Within the study area, simultaneously (or sequentially) with seismic exploration, electrical exploration is carried out in shallow modification by sounding the formation of the field (mSSS) along the lines of seismic profiles with a high spatial-temporal density.
Выполняют обработку и инверсию сигналов мЗСБ. Строят геоэлектрическую модель ВЧР. Восстанавливают скоростную модель ВЧР за счет применения эмпирической зависимости Фауста.Perform the processing and inversion of the signals mzsb. Build a geoelectric model of the VChR. Restore the high-speed model of the VChR through the use of the empirical dependence of Faust.
На основе полученных данных рассчитывают статические поправки.Based on the data obtained, static corrections are calculated.
Заявляемый способ был опробован в различных геологических условиях. Положительный опыт при геологоразведочных работах получен в пределах Западной и Восточной Сибири.The inventive method was tested in various geological conditions. Positive experience in exploration was obtained within Western and Eastern Siberia.
ПРИМЕР.EXAMPLE.
Заявляемый способ был опробован на одном из месторождений углеводородов юга Сибирской платформы.The inventive method was tested on one of the hydrocarbon deposits in the south of the Siberian platform.
На участке исследований была выполнена сейсмическая съемка методом MOB ОГТ 3D по широкоазимутальной сети наблюдений. Параллельно выполнена электроразведка методом мЗСБ, проведена обработка, инверсия, получена геоэлектрическая модель. На исследуемом участке выполнено ВСП в одной из скважин.A seismic survey using the MOB CDP 3D method using a wide-azimuth observation network was performed at the research site. In parallel, electrical prospecting was performed using the MZSB method, processing, inversion was performed, and a geoelectric model was obtained. VSP in one of the wells was performed at the studied site.
Следующим шагом согласно графу, изображенному на фиг. 1, выполнялась калибровка коэффициентов уравнения Фауста. С использованием данных ВСП производился подбор коэффициентов уравнения Фауста, результат подбора иллюстрируют графики на фиг. 2 и фиг. 3. В таблице 1 представлены результаты расчета эмпирических коэффициентов на основе уравнения Фауста, где УЭС - удельное электрическое сопротивление, Vp - скорость продольных волн).The next step according to the graph shown in FIG. 1, the coefficients of the Faust equation were calibrated. Using the VSP data, the coefficients of the Faust equation were selected, and the graphs in FIG. 2 and FIG. 3. Table 1 presents the results of calculating empirical coefficients based on the Faust equation, where resistivity is the electrical resistivity, Vp is the velocity of longitudinal waves).
Далее, согласно полученным коэффициентам, выполнен пересчет всего объема геоэлектрических моделей в скоростные модели. Рассчитаны статические поправки.Further, according to the obtained coefficients, the entire volume of geoelectric models was recalculated into high-speed models. Static corrections are calculated.
Следующим шагом был выполнен ввод статических поправок в сейсмограммы ОГТ (фиг. 4). Статические поправки вводились в двух вариантах:The next step was the introduction of static corrections into the CDP gathers (Fig. 4). Static corrections were introduced in two versions:
1. Вариант, рассчитанный по методу первых вступлений (традиционная методика).1. Option calculated by the method of first introductions (traditional method).
2. Вариант, полученный на основе скоростной модели по даннйм мЗСБ.2. The option obtained on the basis of the speed model according to the data of the Ministry of Health and Social Security.
Анализ временных разрезов показывает, что на этапе учета статики за рельеф и среднечастотных поправок отмечается существенное улучшение прослеживаемости отражающих горизонтов при использовании модели по данным мЗСБ. На разрезе, полученном с использованием модели по данным преломленных волн (по методу первых вступлений), отмечается наличие аномалий-теней, проходящих через весь разрез. На разрезе, полученном по данным мЗСБ, такие аномалии удается в значительной степени подавить.The analysis of time sections shows that at the stage of taking statics for relief and mid-frequency corrections into account, there is a significant improvement in the traceability of reflecting horizons when using the model according to the MSSB data. In the section obtained using the model according to the data of refracted waves (by the method of first arrivals), the presence of shadow anomalies passing through the entire section is noted. In the section obtained according to the MZSB, such anomalies can be largely suppressed.
После учета высокочастотной составляющей на разрезе, полученном с учетом данных мЗСБ, наблюдается улучшение динамики и прослеживаемости отражающих горизонтов (см. фиг. 5).After taking into account the high-frequency component in the section obtained taking into account the data of the MSSS, an improvement in the dynamics and traceability of reflecting horizons is observed (see Fig. 5).
Полученные экспериментальные данные показывают высокую эффективность предложенного способа расчета статических поправок и фактическое повышение достоверности структурных построений при проведении сейсморазведочных работ на участках со сложным строением верхней части разреза.The obtained experimental data show the high efficiency of the proposed method for calculating static corrections and the actual increase in the reliability of structural structures during seismic surveys in areas with a complex structure of the upper part of the section.
Claims (1)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2019121747A RU2722861C1 (en) | 2019-07-09 | 2019-07-09 | Static corrections calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2019121747A RU2722861C1 (en) | 2019-07-09 | 2019-07-09 | Static corrections calculation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| RU2722861C1 true RU2722861C1 (en) | 2020-06-04 |
Family
ID=71067838
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2019121747A RU2722861C1 (en) | 2019-07-09 | 2019-07-09 | Static corrections calculation method |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2722861C1 (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2740509C1 (en) * | 2020-05-25 | 2021-01-14 | Антон Тимурович Семашев | Method for determining optimum depth of charge immersion and constructing model of upper part of the section based on drilling data during seismic survey |
| CN112305624A (en) * | 2020-10-30 | 2021-02-02 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | Non-contact obstacle-crossing high-density measurement method |
| CN112394412A (en) * | 2020-10-30 | 2021-02-23 | 中国石油天然气集团有限公司 | Shear wave static correction method and device based on shallow layer superposition modeling |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU335646A1 (en) * | В. С. Карапыш Трест Краснодарнефтегеофизика | METHOD FOR DETERMINING STATIC AMENDMENTS | ||
| WO2008106828A1 (en) * | 2007-03-08 | 2008-09-12 | Junchang Liu | A high resolution electromagnetic method of earth for removing static frequency domain |
| CN101393270A (en) * | 2007-09-21 | 2009-03-25 | 中国石油集团东方地球物理勘探有限责任公司 | Continuum area surface layer high precision static correction method |
| RU2411547C1 (en) * | 2009-07-13 | 2011-02-10 | Общество с ограниченной ответственностью "Инженерно-геофизическая компания" | Method of defining static correstions |
| CN104133245A (en) * | 2014-07-22 | 2014-11-05 | 中国石油天然气集团公司 | Seismic data static correction method and system |
| RU2614346C2 (en) * | 2015-08-06 | 2017-03-24 | Федеральное государственное унитарное предприятие "Всероссийский научно-исследовательский институт геологии и минеральных ресурсов Мирового океана имени академика И.С. Грамберга" | Method of seismic survey and electrical survey data complex interpretation when searching for offshore hydrocarbon deposits |
-
2019
- 2019-07-09 RU RU2019121747A patent/RU2722861C1/en active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU335646A1 (en) * | В. С. Карапыш Трест Краснодарнефтегеофизика | METHOD FOR DETERMINING STATIC AMENDMENTS | ||
| WO2008106828A1 (en) * | 2007-03-08 | 2008-09-12 | Junchang Liu | A high resolution electromagnetic method of earth for removing static frequency domain |
| CN101393270A (en) * | 2007-09-21 | 2009-03-25 | 中国石油集团东方地球物理勘探有限责任公司 | Continuum area surface layer high precision static correction method |
| RU2411547C1 (en) * | 2009-07-13 | 2011-02-10 | Общество с ограниченной ответственностью "Инженерно-геофизическая компания" | Method of defining static correstions |
| CN104133245A (en) * | 2014-07-22 | 2014-11-05 | 中国石油天然气集团公司 | Seismic data static correction method and system |
| RU2614346C2 (en) * | 2015-08-06 | 2017-03-24 | Федеральное государственное унитарное предприятие "Всероссийский научно-исследовательский институт геологии и минеральных ресурсов Мирового океана имени академика И.С. Грамберга" | Method of seismic survey and electrical survey data complex interpretation when searching for offshore hydrocarbon deposits |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2740509C1 (en) * | 2020-05-25 | 2021-01-14 | Антон Тимурович Семашев | Method for determining optimum depth of charge immersion and constructing model of upper part of the section based on drilling data during seismic survey |
| CN112305624A (en) * | 2020-10-30 | 2021-02-02 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | Non-contact obstacle-crossing high-density measurement method |
| CN112394412A (en) * | 2020-10-30 | 2021-02-23 | 中国石油天然气集团有限公司 | Shear wave static correction method and device based on shallow layer superposition modeling |
| CN112305624B (en) * | 2020-10-30 | 2024-02-06 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | Non-contact obstacle-crossing high-density measurement method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Maraschini et al. | A new misfit function for multimodal inversion of surface waves | |
| US10061046B2 (en) | Integrated passive and active seismic surveying using multiple arrays | |
| Foti et al. | Experiments of joint acquisition of seismic refraction and surface wave data | |
| RU2722861C1 (en) | Static corrections calculation method | |
| EP2972502B1 (en) | System for seismic surveying of a subsurface volume | |
| Groves et al. | Use of geophysical methods for soil profile evaluation | |
| US20170285195A1 (en) | Integrating vertical seismic profile data for microseismic anisotropy velocity analysis | |
| RU2563323C1 (en) | Method of reconstruction of fine structure of geological object and forecast of its fluid saturation | |
| Tounkara et al. | Analyzing the seismic attributes, structural and petrophysical analyses of the Lower Goru Formation: A case study from Middle Indus Basin Pakistan | |
| Gajek et al. | Results of the downhole microseismic monitoring at a pilot hydraulic fracturing site in Poland—Part 1: Event location and stimulation performance | |
| Muhamad et al. | Analysis of borehole geophysical data from the Mora area of the Siljan Ring impact structure, central Sweden | |
| RU2482519C2 (en) | Geophysical survey method | |
| Huang et al. | Fast-forward modeling of compressional arrival slowness logs in high-angle and horizontal wells | |
| Yang et al. | Application of AVO analysis to gas hydrates identification in the northern slope of the South China Sea | |
| Jahani et al. | Limits of 3D detectability and resolution of LWD deep-sensing borehole electromagnetic measurements acquired in the Norwegian Continental Shelf | |
| RU2411547C1 (en) | Method of defining static correstions | |
| Ezersky et al. | Integrated study of the sinkhole development site on the Western shores of the Dead Sea using geophysical methods | |
| Melnikov et al. | Seismic studies of frozen ground in Arctic areas | |
| Kurlenya et al. | Procedure and evidence of seismic research into physical properties of cohesive soils | |
| Farghal | Fault and fracture identification and characterization in 3D seismic data from unconventional reservoirs | |
| RU2758148C1 (en) | Method for searching and controlling hydrocarbons by a complex of geophysical methods | |
| RU2718137C1 (en) | Method of estimating wave field attenuation parameter for determining hydrocarbon saturation of formation in borehole space in constructing geological model | |
| Larson et al. | Surface-wave-based inversions of shallow seismic structure | |
| RU2201606C1 (en) | Method of typification and correlation of oil and gas productive rocks by borehole spectral-time parameters | |
| Ormos et al. | Parallel inversion of refracted travel times of P and SH waves using a function approximation |