RU2445651C2 - Method to build seismic depth and/or time section "kong-macro" (versions) - Google Patents
Method to build seismic depth and/or time section "kong-macro" (versions) Download PDFInfo
- Publication number
- RU2445651C2 RU2445651C2 RU2010122578/28A RU2010122578A RU2445651C2 RU 2445651 C2 RU2445651 C2 RU 2445651C2 RU 2010122578/28 A RU2010122578/28 A RU 2010122578/28A RU 2010122578 A RU2010122578 A RU 2010122578A RU 2445651 C2 RU2445651 C2 RU 2445651C2
- Authority
- RU
- Russia
- Prior art keywords
- seismic
- points
- point
- displayed
- given
- Prior art date
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Область техникиTechnical field
Изобретение относится к области геофизических методов сейсморазведки и может быть использовано при поисках нефти и газа, других полезных ископаемых, а также при инженерных исследованиях грунтов под строительство сооружений.The invention relates to the field of geophysical methods of seismic exploration and can be used in the search for oil and gas, other minerals, as well as in engineering studies of soils for the construction of structures.
Уровень техникиState of the art
Известен способ построения сейсмических временных разрезов A(x,t) по методу общей глубинной точки (МОГТ), где А - интенсивность документирования отображаемой точки на разрезе, x - координата ОТ вдоль профиля, t - время регистрации сейсмического сигнала. МОГТ является основным и повсеместно применяемым способом документирования результатов сейсморазведки, запатентован Мэйном в 1950 году в США и описан во всех учебниках и справочниках [Сейсморазведка. Справочник геофизика. М.: Недра. 1981, с.224-226]. По способу МОГТ на временном разрезе при различных углах засветки из n пунктов воздействий (ПВ) в координатах «длина расстановки - время» отображаемую точку (ОТ) документируют под серединой расстановки (расстоянием между ПВ и сейсмоприемником), а время определяют вычислением соответствующего виртуального годографа общей глубинной точки ОГТ. Причем трасса хода лучей всегда предполагается в виде опрокинутого равнобедренного треугольника, на вершине которого визуализируют ОТ, а среда предполагает зеркальные отражения от горизонтально залегающих границ между слоями горных пород. Затем используют различные процедуры миграции для перерасчета временного разреза в глубинный.There is a known method for constructing seismic time sections A (x, t) using the common depth point method (MOGT), where A is the intensity of documenting the displayed point in the section, x is the OT coordinate along the profile, t is the time of registration of the seismic signal. MOGT is the main and widely used way of documenting the results of seismic surveys, patented by Maine in 1950 in the USA and described in all textbooks and reference books [Seismic exploration. Handbook of geophysics. M .: Subsoil. 1981, p. 224-226]. According to the MOGT method, in the time section at various exposure angles from n exposure points (PV) in the coordinates “placement length - time”, the displayed point (OT) is documented under the middle of the placement (the distance between the PV and the geophone), and the time is determined by calculating the corresponding virtual travel time curve for the total the deepest point of the CDP. Moreover, the ray path is always assumed in the form of an overturned isosceles triangle, at the top of which an OT is visualized, and the medium assumes mirror reflections from horizontally lying boundaries between rock layers. Then, various migration procedures are used to recalculate the time section into the deep one.
Недостатками известного способа являются отображение наклонных и криволинейных сейсмических горизонтов со сносом относительно их истинного положения на разрезе и небольшая кратность накопления регистрируемых отраженных от общей глубинной точки сейсмических сигналов, равная кратности профилирования, то есть числу углов засветки от n источников сейсмических воздействий (ПВ). Кроме того, по способу МОГТ из всего множества сигналов на сейсмограмме для документирования ОТ на временном разрезе используют лишь одну единственную величину амплитуды сигнала (m=ρ=1, где m - число сейсмотрасс, используемых для отображения точки, ρ - число отсчетов фаз цуга колебаний на сейсмотрассе). Это приводит к получению недостоверных результатов разведки сложно построенных сред и месторождений нефти и газа, неудовлетворительному латеральному отображению сейсмических границ и использованию менее 0.001% информации зарегистрированного на сейсмограммах волнового поля. Кроме того, при миграции временного разреза в глубинный A(x,h), где h - глубина, многочисленные пересечения годографов различных типов волн существенно искажают отображение глубинного разреза.The disadvantages of this method are the display of inclined and curved seismic horizons with drift relative to their true position in the section and the small multiplicity of accumulation of recorded seismic signals reflected from the common depth point, which is equal to the profiling ratio, i.e. the number of exposure angles from n sources of seismic effects (PS). In addition, according to the MOGT method, from the entire set of signals in a seismogram, only one single value of the signal amplitude is used to document the OT in the time section (m = ρ = 1, where m is the number of seismic traces used to display the point, ρ is the number of phase samples of the oscillation train on seismic trace). This leads to unreliable results of exploration of complexly constructed media and oil and gas fields, unsatisfactory lateral display of seismic boundaries and the use of less than 0.001% of the information recorded in the seismograms of the wave field. In addition, during the migration of a temporary section into the deep A (x, h), where h is the depth, numerous intersections of the hodographs of various types of waves significantly distort the display of the deep section.
Наиболее близким к предлагаемому изобретению является метод регулируемого направленного приема (МРНП) [Сейсморазведка. Справочник геофизика. М.: Недра. 1981, с.222-224], основанный на представлении о том, что от сейсмических границ отражаются плоские интерференционные и регулярные волны различных типов, одновременно приходящие к профилю по разным направлениям. Средством разрешения и расщепления волн и регистрации их на временном разрезе A(x,t) является разновременное прямолинейное суммирование каждой отдельной (до ρ) фазы цуга колебаний (где ρ - число отсчетов фаз цуга колебаний на сейсмотрассе) не только на зарегистрированных m сейсмотрассах, но и так называемое «многомерное» дополнительное разновременное суммирование фаз колебаний на коротких приемных базах, состоящих из m суммотрасс, зарегистрированных методом многократного профилирования и предварительно просуммированных по способу метода общей глубинной точки (МОГТ).Closest to the proposed invention is a method of controlled directional reception (MRNP) [Seismic exploration. Handbook of geophysics. M .: Subsoil. 1981, p.222-224], based on the notion that plane interference and regular waves of various types are reflected from seismic boundaries at the same time coming to the profile in different directions. A means of resolving and splitting the waves and recording them in the time section A (x, t) is a straight-time summation of each separate (up to ρ) phase of the oscillation train (where ρ is the number of phase samples of the oscillation train at the seismic trace) not only on the registered m seismic traces, but and the so-called “multidimensional” additional multi-simultaneous summation of the oscillation phases at short receiving bases, consisting of m summates registered by the multiple profiling method and previously summed by the method a common deep point (MOGT).
Недостатками способа МРНП являются отображение наклонных сейсмических горизонтов на временном разрезе со сносом относительно их истинного положения на разрезе и небольшая кратность накопления регистрируемых сейсмических сигналов, обусловленная требованием прямолинейности участка годографа на протяжении всей базы суммирования. Кроме того, по способу МРНП из всего множества сигналов на сейсмограмме для документирования ОТ на временном разрезе используют лишь одну единственную величину амплитуды фазы сигнала (ρ=1, где ρ - число отсчетов фаз цуга колебаний на сейсмотрассе). Эти недостатки приводят к получению недостоверных результатов разведки сложно построенных сред и месторождений нефти и газа, неудовлетворительному латеральному отображению сейсмических границ и использованию менее 0.01% информации зарегистрированного на сейсмограммах волнового поля.The disadvantages of the MRI method are the display of inclined seismic horizons in the time section with drift relative to their true position in the section and the small accumulation rate of the recorded seismic signals, due to the requirement of straightness of the hodograph section throughout the summation base. In addition, according to the MRI method, out of the entire set of signals in the seismogram, only one single value of the signal phase amplitude is used to document the OT in the time section (ρ = 1, where ρ is the number of phase samples of the oscillation train at the seismic trace). These shortcomings lead to unreliable results of exploration of complexly constructed media and oil and gas fields, unsatisfactory lateral imaging of seismic boundaries and the use of less than 0.01% of the information recorded on the seismograms of the wave field.
К недостаткам известного способа МРНП относится также и то, что в результате сейсмического воздействия из источника (ПВ) на m пунктов предполагаемой сейсмической границы для построения сейсмического разреза при использовании из сейсмограммы m сейсмотрасс, зарегистрированных m сейсмоприемниками, не предусмотрен способ задания поиска их местоположения относительно координат отображаемой точки OT(x0,h). Это приводит к уменьшению степени достоверности построения разреза и дополнительным затратам ресурсов на оценку и упорядочение информации, появляющейся при произвольном задании положения и размеров зоны обзора отражающей наклонной сейсмической границы. Кроме того, необходимы затраты ресурсов на поиск и обнуление (мьютинг) избыточной информации, препятствующей оперативному построению глубинных и/или временных разрезов, а также то, что при миграции временного разреза в глубинный многочисленные пересечения и петли годографов различных типов волн существенно искажают отображение глубинного разреза.The disadvantages of the known method of MRI also include the fact that as a result of seismic action from the source (PV) at m points of the alleged seismic boundary for constructing a seismic section when using m seismic records recorded by m seismic receivers, there is no way to set the search for their location relative to the coordinates displayed point OT (x 0 , h). This leads to a decrease in the reliability of the construction of the section and the additional cost of resources for evaluating and organizing the information that appears when the position and size of the viewing area of the reflecting inclined seismic boundary are arbitrarily set. In addition, it requires resources to search for and zero out (excess) information that impedes the operational construction of deep and / or temporal sections, as well as the fact that when a temporal section is migrated to a deep section, numerous intersections and hodograph loops of various types of waves significantly distort the display of a deep section .
Данный способ МРНП избран в качестве прототипа.This method of MRNP is selected as a prototype.
Прототип и заявленный способ имеют следующие общие признаки:The prototype and the claimed method have the following common features:
- используют данные о глубинном строении недр из n сейсмограмм, полученных методом многократного профилирования, обеспечивающим засветку отображаемых точек под n углами из n пунктов воздействий,- use data on the deep structure of the subsoil from n seismograms obtained by the multiple profiling method, providing illumination of the displayed points at n angles from n impact points,
- строят виртуальные годографы для каждой из m сейсмотрасс из зарегистрированных n сейсмограмм,- build virtual hodographs for each of m seismic traces from registered n seismograms,
- производят суммирование отдельных отсчетов амплитуд фазы с кратностью до mn, полученных после построения виртуальных годографов для каждой из n сейсмотрас из зарегистрированных n сейсмограмм.- summarize individual samples of phase amplitudes with a multiplicity of up to mn obtained after constructing virtual hodographs for each of n seismic traces from registered n seismograms.
Раскрытие изобретенияDisclosure of invention
Задачей изобретения является разработка СПОСОБА ПОСТРОЕНИЯ СЕЙСМИЧЕСКИХ ГЛУБИННЫХ И/ИЛИ ВРЕМЕННЫХ РАЗРЕЗОВ, обеспечивающего высокую достоверность и точность отображения сейсмических объектов и сред большой сложности, при котором для обзора полученных в результате воздействия от источника (ПВ) сейсмограмм волнового поля в виде совокупности m сейсмотрасс, зарегистрированных m сейсмоприемниками (СП), оптимальным путем осуществляется однозначная (а не произвольная) привязка координат СП в зависимости от местоположения отображаемой точки (ОТ) на предполагаемой зеркальной сейсмической границе с заданным углом наклона в ОТ и особенностей геолого-геофизической обстановки.The objective of the invention is to develop a METHOD of CONSTRUCTION of SEISMIC DEPTH AND / OR TIME SECTIONS, providing high reliability and accuracy of the display of seismic objects and medium complexity, in which to review the resulting wave field seismograms in the form of a combination of m seismic records recorded m geophones (SP), in the optimal way, an unambiguous (and not arbitrary) binding of the coordinates of the SP is carried out depending on the location of the displayed point and (OT) at the supposed mirror seismic boundary with a given angle of inclination in the OT and features of the geological and geophysical situation.
Технический результат состоит в увеличении в десятки раз отношения сигнал/шум и степени подавления случайных и регулярных волн/помех, исключении влияния сноса и петель на годографах, сокращении затрат ресурсов и в увеличении степени достоверности документирования заданных отрезков наклонных сейсмических границ при построении глубинного разреза или его фрагмента.The technical result consists in increasing the signal-to-noise ratio tens of times and the degree of suppression of random and regular waves / noise, eliminating the influence of drift and loops on the hodographs, reducing the cost of resources and increasing the degree of reliability of documenting specified segments of inclined seismic boundaries when constructing a deep section or its fragment.
Поставленная задача решена в двух вариантах за счет того, что в известном способе, при котором используют данные о глубинном строении недр из n сейсмограмм, полученных методом многократного профилирования, обеспечивающим засветку отображаемых точек под n углами из n пунктов воздействий, строят виртуальные годографы для каждой из m сейсмотрасс из зарегистрированных n сейсмограмм, производят суммирование отдельных отсчетов амплитуд фазы с кратностью до mn, полученных после построения виртуальных годографов для каждой из m сейсмотрасс из зарегистрированных n сейсмограмм, согласно предлагаемому способу построение сейсмического глубинного и/или временного разреза осуществляют последовательной разверткой отображаемых точек (ОТ) с заданным шагом по столбцам или строкам с использованием априорно заданного скоростного разреза V=V0(x,h). Задают диапазон поиска наиболее вероятного значения углов наклона в ОТ отражающего сейсмического объекта (сейсмической границы) заданной конфигурации q=q0+Δq(1, 2, …, N/2, …, N), где q - текущее значение угла наклона, q0 - начальное априорное значение угла наклона отражающей границы в ОТ, Δq - шаг приращения угла наклона, N - текущий номер угла наклона. Для каждого из N ракурсов объекта рассчитывают (генерируют) виртуальный годограф заданного типа сейсмической волны с учетом местоположения (координат) ПВ, ОТ и закона (функции) расположения последовательности из m прилегающих к ОТ пунктов (точек) отражения.The problem is solved in two versions due to the fact that in the known method, which uses data on the deep structure of the subsoil from n seismograms obtained by the multiple profiling method, which illuminates the displayed points at n angles from n exposure points, build virtual hodographs for each of m seismic traces from registered n seismograms, summation of individual samples of phase amplitudes with a multiplicity of up to mn obtained after constructing virtual hodographs for each of m seismic traces due to egistrirovannyh n seismograms, according to the proposed method for the construction of seismic depth and / or time section is carried out sequential scanning of display dots (ON) at a predetermined pitch along the columns or rows with a priori predetermined velocity profile V = V 0 (x, h). The search range of the most probable value of the tilt angles in the OT of the reflecting seismic object (seismic boundary) of the given configuration q = q 0 + Δq (1, 2, ..., N / 2, ..., N) is set, where q is the current value of the tilt angle, q 0 is the initial a priori value of the angle of inclination of the reflecting boundary in OT, Δq is the increment step of the angle of inclination, N is the current number of the angle of inclination. For each of the N angles of the object, a virtual hodograph of a given type of seismic wave is calculated (generated) taking into account the location (coordinates) of the PV, OT and the law (function) of the location of a sequence of m reflection points (points) adjacent to the OT.
Находят координаты местоположения m сейсмоприемников, на которых зарегистрированы сейсмические сигналы на m сейсмотрассах после прихода виртуальных восстающих лучей из m пунктов (точек), прилегающих к ОТ вдоль искомой сейсмической границы.The coordinates of the location of m seismic receivers are found, on which seismic signals are recorded at m seismic traces after the arrival of virtual rising rays from m points (points) adjacent to the OT along the desired seismic boundary.
Последовательность из m сейсмоприемников задают на сейсмическом профиле по варианту 1 так, чтобы центр х▲ц заданной численно последовательности из m СП располагался симметрично относительно точки сноса 
По варианту 2 для обеспечения обзора прилегающего к ОТ участка сейсмической границы на заданном от точки сноса
На каждой сейсмотрассе находят с учетом скоростного разреза маркерные отметки времени прихода сигнала после каждой генерации соответствующего виртуального годографа. На зарегистрированных m сейсмотрассах каждой сейсмограммы относительно маркерной отметки времени (например, до и после нее) считывают по ρ компонентов (обычно ρ~20) в виде псевдосинфазных отсчетов величин амплитуд и производят направленное псевдосинфазное накопление по ρm отсчетов фаз цугов колебаний из каждой из m сейсмотрасс для каждого из q углов наклона. После синфазного накопления по ρm амплитуд колебаний сигналов из каждой из n сейсмограмм (с привязкой ко всем маркерным отметкам времени mnN годографов) сигналы группируют по признаку равенства номера шага позиционирования N, в каждой из N групп производят супернакопления из ρmn амплитуд сигналов после введения поправочных коэффициентов за угол засветки известными способами. Для каждого из N супернакоплений по pmn амплитуд, специфичных для заданных q позиций объекта в ОТ, известным способом вычисляют N функционалов статистической обработки, сравнивают их между собой, по максимальной величине одного из функционалов определяют наиболее вероятный номер N углов наклона объекта. Максимум величин функционалов сравнивают с заданным порогом обнаружения, при превышении порога обнаружения ОТ наносят на глубинный сейсмический разрез с интенсивностью, пропорциональной величине максимального функционала. Аналогично разверткой по столбцам и/или строкам наносят все остальные ОТ данной версии глубинного разреза с интенсивностью, соответствующей максимуму вероятности обнаружения на зарегистрированных сейсмограммах отсчетов амплитуд волн заданного типа. После построения данной версии глубинного разреза в необходимых случаях приступают к получению иных его версий с другим законом зависимости скорости от глубины, иного типа сейсмогеологической модели с измененными параметрами обзора и/или генерации годографов других типов волн.At each seismic trace, taking into account the velocity section, marker marks of the signal arrival time after each generation of the corresponding virtual hodograph are found. On the registered m seismic traces of each seismogram relative to the marker of the time (for example, before and after it), ρ components (usually ρ ~ 20) are read in the form of pseudo-phase samples of amplitudes and directional pseudo-phase accumulation of ρm samples of phases of the oscillation trains from each of m seismic surveys is produced for each of q tilt angles. After the in-phase accumulation over ρm of the amplitudes of the oscillations of the signals from each of the n seismograms (with reference to all marker time stamps mnN of the hodographs), the signals are grouped according to the equality of the positioning step number N, in each of the N groups produce super-accumulations of the ρmn signal amplitudes after introducing correction factors for flare angle by known methods. For each of the N super accumulations of pmn amplitudes specific for the given q positions of the object in the RT, N functionals of statistical processing are calculated in a known manner, they are compared with each other, the most probable number N of the object tilt angles is determined by the maximum value of one of the functionals. The maximum values of the functionals are compared with a given detection threshold; when the detection threshold is exceeded, the RTs are applied to the deep seismic section with an intensity proportional to the maximum functional. Similarly, by scanning in columns and / or rows, all the other FROM a given version of a deep section are plotted with an intensity corresponding to the maximum probability of detection of recorded amplitudes of waves of a given type on recorded seismograms. After constructing this version of the deep section, if necessary, they begin to obtain other versions of it with a different law of the dependence of speed on depth, of a different type of seismic-geological model with modified parameters for the survey and / or generation of hodographs of other types of waves.
Способ осуществляется следующим образом.The method is as follows.
Построение сейсмического глубинного и/или временного разреза включаетThe construction of a seismic depth and / or time section includes
- использование данных о глубинном строении недр из n сейсмограмм, полученных методом многократного профилирования, обеспечивающим засветку отображаемых точек под n углами из n пунктов воздействий,- the use of data on the deep structure of the subsoil from n seismograms obtained by the multiple profiling method, providing illumination of the displayed points at n angles from n impact points,
- построение виртуальных годографов для каждой из m сейсмотрасс из зарегистрированных n сейсмограмм,- construction of virtual hodographs for each of m seismic traces from registered n seismograms,
- суммирование кратностью до mn отдельных отсчетов амплитуд фазы, зарегистрированных на заданных m сейсмотрассах каждой из n сейсмограмм.- summation up to mn of individual samples of phase amplitudes recorded on given m seismic traces of each of n seismograms.
Документирование всех ОТ разреза осуществляют путем последовательной развертки отображаемых точек (ОТ) с заданным шагом по столбцам или строкам с использованием априорно заданного скоростного разреза V=V0(x,h). Задают диапазон поиска наиболее вероятного значения углов наклона в ОТ отражающего сейсмического объекта (сейсмической границы) заданной конфигурации q=q0+Δq(1, 2, …, N/2, …, N), где q текущее значение угла наклона, q0 - начальное априорное значение угла наклона отражающей границы в ОТ, Δq - шаг приращения угла наклона, N - текущий номер угла наклона. Для каждого из N ракурсов объекта рассчитывают (генерируют) виртуальный годограф заданного типа сейсмической волны с учетом местоположения (координат) ПВ, ОТ и закона (функции) расположения последовательности из m прилегающих к ОТ пунктов (точек) отражения.Documentation of all OT sections is carried out by sequentially scanning the displayed points (OT) with a given step in columns or rows using an a priori specified velocity section V = V 0 (x, h). The search range of the most probable value of the tilt angles in the OT of the reflecting seismic object (seismic boundary) of the given configuration q = q 0 + Δq (1, 2, ..., N / 2, ..., N) is set, where q is the current value of the tilt angle, q 0 is the initial a priori value of the angle of inclination of the reflecting boundary in OT, Δq is the increment of the angle of inclination, N is the current number of the angle of inclination. For each of the N angles of the object, a virtual hodograph of a given type of seismic wave is calculated (generated) taking into account the location (coordinates) of the PV, OT and the law (function) of the location of a sequence of m reflection points (points) adjacent to the OT.
Находят координаты местоположения m сейсмоприемников, на которых зарегистрированы сейсмические сигналы на m сейсмотрассах после прихода виртуальных восстающих лучей из m пунктов (точек), прилегающих к ОТ вдоль искомой сейсмической границы.The coordinates of the location of m seismic receivers are found, on which seismic signals are recorded at m seismic traces after the arrival of virtual rising rays from m points (points) adjacent to the OT along the desired seismic boundary.
Последовательность из m сейсмоприемников задают на сейсмическом профиле по варианту 1 так, чтобы центр x▲ц заданной численно последовательности из m СП располагался симметрично относительно точки сноса R↑ x (точки выхода на профиль восстающего из ОТ отраженного луча вторичной волны, вызванной ПВ).The sequence of m seismic receivers is set on the seismic profile according to option 1 so that the center x ▲ ленно of the given numerical sequence from m SP is located symmetrically with respect to the drift point R ↑ x (the exit point to the profile of the secondary wave arising from the reflected beam from the OT caused by the PV).
По варианту 2 для обеспечения обзора прилегающего к ОТ участка сейсмической границы на заданном от точки сноса R↑ x расстоянии применяют вынос x▲e последовательности СП, величину и знак которого обуславливают особенностями геолого-геофизической обстановки. Например, изменяя величину выноса, по изменению интенсивности документирования ОТ судят о шероховатости или степени соответствия заданного априорно угла наклона и радиуса кривизны его истинному значению на дистанции выноса относительно ОТ.According to option 2, to provide an overview of the seismic boundary section adjacent to the OT at a distance set from the drift point R ↑ x, the x ▲ e offset of the SP sequence is used, the magnitude and sign of which are determined by the characteristics of the geological and geophysical situation. For example, by changing the amount of drift, by changing the intensity of documentation of RTs, one judges the roughness or degree of compliance of a given a priori tilt angle and radius of curvature with its true value at the distance of the drift relative to RT.
На каждой сейсмотрассе находят с учетом скоростного разреза маркерные отметки времени прихода сигнала после каждой генерации соответствующего виртуального годографа. На зарегистрированных m сейсмотрассах каждой сейсмограммы относительно маркерной отметки времени (например, до и после нее) считывают по ρ компонентов (обычно ρ~20) в виде псевдосинфазных отсчетов величин амплитуд и производят направленное псевдосинфазное накопление по ρm отсчетов фаз цугов колебаний из каждой из m сейсмотрасс для каждого из q углов наклона. После синфазного накопления по ρm амплитуд колебаний сигналов из каждой из n сейсмограмм (с привязкой ко всем маркерным отметкам времени на mnN годографах) сигналы группируют по признаку равенства номера шага позиционирования N, в каждой из N групп производят супернакопления из ρmn амплитуд сигналов после введения поправочных коэффициентов за угол засветки известными способами. Для каждого из N супернакоплений по ρmn амплитуд, специфичных для заданных q позиций объекта в ОТ, известным способом вычисляют N функционалов статистической обработки, сравнивают их между собой, по максимальной величине одного из функционалов определяют наиболее вероятный номер N углов наклона объекта. Максимум величин функционалов сравнивают с заданным порогом обнаружения, при превышении порога обнаружения ОТ наносят на глубинный сейсмический разрез с интенсивностью, пропорциональной величине максимального функционала. Аналогично разверткой по столбцам и/или строкам наносят все остальные ОТ данной версии глубинного разреза с интенсивностью, соответствующей максимуму вероятности обнаружения на зарегистрированных сейсмограммах отсчетов амплитуд волн заданного типа. После построения данной версии глубинного разреза в необходимых случаях приступают к получению иных его версий с другим законом зависимости скорости от глубины, иного типа сейсмогеологической модели с измененными параметрами обзора и/или генерации годографов других типов волн.At each seismic trace, taking into account the velocity section, marker marks of the signal arrival time after each generation of the corresponding virtual hodograph are found. On the registered m seismic traces of each seismogram relative to the marker of the time (for example, before and after it), ρ components (usually ρ ~ 20) are read in the form of pseudo-phase samples of amplitudes and directional pseudo-phase accumulation of ρm samples of phases of the oscillation trains from each of m seismic surveys is produced for each of q tilt angles. After the in phase accumulation over ρm of the amplitudes of the oscillations of the signals from each of the n seismograms (with reference to all marker timestamps on the mnN hodographs), the signals are grouped according to the equality of the positioning step number N, in each of the N groups produce super accumulations of ρmn signal amplitudes after the introduction of correction factors around the exposure angle by known methods. For each of the N super-accumulations of ρmn amplitudes specific for the given q positions of the object in the RT, N functionals of statistical processing are calculated in a known manner, they are compared with each other, the most probable number N of the object tilt angles is determined by the maximum value of one of the functionals. The maximum values of the functionals are compared with a given detection threshold; when the detection threshold is exceeded, the RTs are applied to the deep seismic section with an intensity proportional to the maximum functional. Similarly, by scanning in columns and / or rows, all the other FROM a given version of a deep section are plotted with an intensity corresponding to the maximum probability of detection of recorded amplitudes of waves of a given type on recorded seismograms. After constructing this version of the deep section, if necessary, they begin to obtain other versions of it with a different law of the dependence of speed on depth, of a different type of seismic-geological model with modified parameters for the survey and / or generation of hodographs of other types of waves.
Пути реализации патентуемого способа приводятся в нижеследующих описаниях 2-х примеров. Пример 1 поясняется фигурой 1а), пример 2 - фигурой 1б).Ways to implement the patented method are given in the following descriptions of 2 examples. Example 1 is illustrated by figure 1a), example 2 - by figure 1b).
Пример 1 (вариант 1). Реализация способа, при котором центр x▲ц заданной последовательности из m сейсмоприемников размещается симметрично относительно точки сноса R↑ x восстающего луча вторичной отраженной волны из ОТ, расположенной на наклонной сейсмической границе.Example 1 (option 1). The implementation of the method in which the center x ▲ c of a given sequence of m geophones is placed symmetrically with respect to the drift point R ↑ x of the rising beam of the secondary reflected wave from the OT located on the inclined seismic boundary.
Описание примера поясняется рисунком 1а), где:The description of the example is illustrated in Figure 1a), where:
OT(x0,h) - отображаемая точка на наклонной сейсмической границе;OT (x 0 , h) - the displayed point on the inclined seismic boundary;
1 - линия сейсмического профиля;1 - seismic line;
2 - отражающая сейсмическая граница;2 - reflective seismic boundary;
3 - пункт воздействий (ПВ);3 - point impacts (PV);
4 - первый из заданной последовательности сейсмоприемник СП;4 - the first of a given sequence of the seismic receiver SP;
5 - последний сейсмоприемник из заданной последовательности СП;5 - the last seismic receiver from a given sequence of SP;
6 - падающий из ПВ луч на ОТ;6 - incident beam from the PV at the OT;
7 - глубина залегания ОТ по координате h;7 - depth of OT along the coordinate h;
8 - дистанция (величина) сноса R↑ x вдоль профиля восстающего из ОТ луча 9;8 - distance (magnitude) of drift R ↑ x along the profile of ray 9 rising from OT;
9 - восстающий из ОТ отраженный луч;9 - a reflected ray rising from OT;
10 - падающий луч на первую исследуемую прилегающую к ОТ точку 14 на отражающей границе;10 - incident beam to the first investigated adjacent to OT point 14 on the reflective border;
11 - падающий луч на последнюю из заданной последовательности исследуемую точку 15, прилегающую к ОТ;11 - incident beam to the last of a given sequence of the investigated point 15 adjacent to the OT;
12 - восстающий из точки 14 отраженный луч;12 - the reflected ray rising from the point 14;
13 - восстающий из последней прилегающей к ОТ точки 15 луч;13 - rising from the last ray adjacent to the OT point 15;
14 - первая из прилегающих к ОТ точек;14 - the first of the points adjacent to the OT;
15 - последняя из прилегающих к ОТ точек;15 - the last of the points adjacent to the OT;
16 - исследуемый отрезок сейсмической границы на дистанции от 14 до 15 точки;16 - the investigated segment of the seismic boundary at a distance of 14 to 15 points;
17 - половина расстояния между первым (4-м) и последним (5-м) сейсмоприемниками;17 - half the distance between the first (4th) and last (5th) geophones;
18 - центр х▲ц заданной последовательности сейсмоприемников.18 - center x ▲ q of a given sequence of geophones.
Геометрия сейсмической расстановки на фиг.1а) состоит из сейсмического профиля 1, отражающей границы 2, на которую падает сейсмическое излучение из источника воздействий 3 (ПВ). Последовательность используемых для исследований отражений состоит из сейсмоприемников «4-5». Из ПВ 3 луч 6 падает на OT(x0,h), расположенную по вертикали 7 на глубине h, и, отражаясь от границы 2, выходит на профиль 1 со сносом 8 восстающего луча 9. Для реализации супернакоплений сигналов от сейсмической границы необходимо, чтобы сигналы от точек, прилегающих к ОТ, регистрировались симметрично относительно их расположения вблизи ОТ на сейсмической границе. Для решения этой задачи в соответствии с предлагаемым способом необходимо, чтобы падающие граничные лучи 10 и 11 отражались в виде восстающих лучей 12 и 13 из крайних точек 14 и 15 исследуемого отрезка 16 сейсмической границы на последовательность сейсмоприемников «4-5» симметрично относительно восстающего луча 9. Поэтому по заданным координатам х0 и h отображаемой точки ОТ для каждого угла наклона границы в ОТ известными способами сначала рассчитывают координаты до точки выхода на профиль восстающего луча 9. Тогда находят расстояние 8 его точки выхода до х0, которое является дистанцией сноса
Далее на каждой сейсмотрассе находят с учетом скоростного разреза маркерные отметки времени прихода сигнала после каждой генерации соответствующего виртуального годографа. На зарегистрированных m сейсмотрассах каждой сейсмограммы относительно mN маркерных отметок времени считывают ρ компонентов (обычно ρ~20) в виде псевдосинфазных отсчетов величин амплитуд на mN годографах и производят направленное псевдосинфазное накопление по ρm отсчетов фаз цугов колебаний из каждой из m сейсмотрасс для каждого из q углов наклона. После синфазного накопления по ρm амплитуд колебаний сигналов из каждой из n сейсмограмм (с привязкой ко всем маркерным отметкам времени mnN годографов) сигналы группируют по признаку равенства номера шага позиционирования N, в каждой из N групп производят супернакопления из ρmn амплитуд сигналов после введения поправочных коэффициентов за угол засветки известными способами. Для каждого из N супернакоплений по ρmn амплитуд, специфичных для заданных q позиций объекта в ОТ, известным способом вычисляют N функционалов статистической обработки, сравнивают их между собой, по максимальной величине одного из функционалов определяют наиболее вероятный номер N позиции объекта. Максимум величин функционалов сравнивают с заданным порогом обнаружения, при превышении порога обнаружения ОТ наносят на глубинный сейсмический разрез с интенсивностью, пропорциональной величине максимального функционала. Аналогично разверткой по столбцам и/или строкам наносят все остальные ОТ данной версии глубинного разреза с интенсивностью, соответствующей максимуму вероятности обнаружения на зарегистрированных сейсмограммах отсчетов амплитуд волн заданного типа. После построения данной версии глубинного разреза в необходимых случаях приступают к получению иных его версий с другим законом зависимости скорости от глубины, иного типа сейсмогеологической модели с измененными параметрами обзора и/или генерации годографов других типов волн.Then, at each seismic trace, taking into account the velocity section, marker marks of the signal arrival time after each generation of the corresponding virtual hodograph are found. On the registered m seismic traces of each seismogram with respect to mN marker timestamps, ρ components (usually ρ ~ 20) are read in the form of pseudo-phase samples of amplitudes on mN hodographs and directional pseudo-phase accumulation of ρm phase samples of oscillation trains from each of m seismic tracks for each of q angles is performed tilt. After the in-phase accumulation over ρm of the amplitudes of the oscillations of the signals from each of the n seismograms (with reference to all marker time stamps mnN of the hodographs), the signals are grouped according to the equality of the positioning step number N, in each of the N groups produce super-accumulations of the ρmn signal amplitudes after introducing correction factors for flare angle by known methods. For each of the N super accumulations of ρmn amplitudes specific for given q positions of the object in the RT, N functionals of statistical processing are calculated in a known manner, they are compared with each other, the most probable number N of the position of the object is determined by the maximum value of one of the functionals. The maximum values of the functionals are compared with a given detection threshold; when the detection threshold is exceeded, the RTs are applied to the deep seismic section with an intensity proportional to the maximum functional. Similarly, by scanning in columns and / or rows, all the other FROM a given version of a deep section are plotted with an intensity corresponding to the maximum probability of detection of recorded amplitudes of waves of a given type on recorded seismograms. After constructing this version of the deep section, if necessary, they begin to obtain other versions of it with a different law of the dependence of speed on depth, of a different type of seismic-geological model with modified parameters for the survey and / or generation of hodographs of other types of waves.
Пример 2 (вариант 2). Реализация способа, при котором осуществляют построение отражающих наклонных элементов глубинного разреза, расположенных на заданной величине удаления (дистанции) от ОТ (с выносом х▲е центра последовательности m СП).Example 2 (option 2). The implementation of the method, in which the construction of reflective inclined elements of the deep section are carried out, located at a predetermined distance (distance) from the OT (with offset x ▲ e the center of the sequence m SP).
Предлагаемый способ поясняется фиг.1б), где OT(x0,h), наименование, назначение и нумерация позиций «1-18» та же, что и на фиг.1а);The proposed method is illustrated figb), where OT (x 0 , h), the name, purpose and numbering of positions "1-18" is the same as in figa);
19 - величина (дистанция) выноса заданного центра 18 х▲ц последовательности сейсмоприемников СП «4-5» относительно точки выхода восстающего из ОТ луча 9 со сносом 8.19 - the magnitude (distance) of the removal of a given center of 18 x ▲ q sequence of seismic receivers SP "4-5" relative to the exit point of the beam 9 rising from the OT with drift 8.
В процессе генерации виртуальных годографов при интерпретации сейсморазведочных разрезов возникает задача прослеживания исследуемой наклонной сейсмической границы со смещением центра симметрии прилегающих к ОТ точек «14-15» на заданное расстояние. По величине вероятности (достоверности) обнаружения на заданной дистанции от ОТ точек «14-15» определяют стабильность простирания и шероховатость данной границы, измеряют радиус кривизны, крутизну падения/восстания пластов горных пород. Кроме того, при наличии зон недоступности для размещения СП (здания, сооружения, реки и прочие препятствия) также необходимо смещение центра 18 СП относительно выхода восстающего из ОТ луча 9 со сносом 8. В таких случаях в соответствии с предлагаемым способом применяют вынос 19 центра последовательности СП на заданную величину (дистанцию) относительно сноса 8 (точки выхода из ОТ восстающего луча 9 на линию профиля). Тогда обеспечивается обзор волнового поля участка сейсмической границы, расположенного на заданном расстоянии от ОТ. Обзор расстояния между последовательностью точек «14-15» на границе обеспечивают последовательностью СП «4-5» при задании N ракурсов сейсмического объекта с определенным шагом поворота Δq с центром вращения в ОТ. Для каждого ракурса объекта q=q0+Δq(1, 2, …, N/2, …, N) рассчитывают виртуальный годограф заданного типа сейсмической волны с учетом местоположения (координат) ПВ, ОТ и заданного закона (функции) расположения последовательности из m расположенных на заданном расстоянии от ОТ последовательности точек «14-15»; производят направленное синфазное накопление амплитуд многокомпонентных колебаний на зарегистрированных m сейсмотрассах последовательности «4-5»; далее по примеру 1 на каждой сейсмотрассе находят с учетом скоростного разреза маркерные отметки времени прихода сигнала после каждой генерации соответствующего виртуального годографа…In the process of generating virtual hodographs when interpreting seismic sections, the problem arises of tracking the investigated inclined seismic boundary with a shift of the center of symmetry of the 14-15 points adjacent to the OT by a given distance. The probability (reliability) of detection at a given distance from the points “14-15” determines the stability of the strike and the roughness of this boundary, measures the radius of curvature, the steepness of the fall / uprising of rock formations. In addition, if there are zones of inaccessibility for the location of the joint venture (buildings, structures, rivers and other obstacles), it is also necessary to shift the center of the 18 joint ventures relative to the exit of ray 9 rising from the OT from the drift 8. In such cases, in accordance with the proposed method, the center of the sequence 19 SP for a given amount (distance) relative to the drift 8 (exit points from the uprising beam 9 to the profile line). Then an overview of the wave field of the seismic boundary located at a given distance from the OT is provided. A review of the distance between the sequence of points “14-15” at the boundary is provided by the sequence of the joint venture “4-5” when setting N angles of the seismic object with a certain rotation step Δq with the rotation center in the OT. For each aspect of the object q = q 0 + Δq (1, 2, ..., N / 2, ..., N), a virtual hodograph of a given type of seismic wave is calculated taking into account the location (coordinates) of the PV, OT and the given law (function) of the location of the sequence from m located at a given distance from the sequence of points “14-15”; directional in-phase accumulation of amplitudes of multicomponent oscillations is performed on recorded m seismic traces of the sequence “4-5”; further according to example 1, at each seismic trace, taking into account the velocity section, marker marks of the signal arrival time after each generation of the corresponding virtual hodograph are found ...
Способ апробирован на синтетических и производственных тестовых примерах и разрезах и показал высокую эффективность.The method was tested on synthetic and industrial test examples and sections and showed high efficiency.
Claims (2)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2010122578/28A RU2445651C2 (en) | 2010-06-02 | 2010-06-02 | Method to build seismic depth and/or time section "kong-macro" (versions) |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2010122578/28A RU2445651C2 (en) | 2010-06-02 | 2010-06-02 | Method to build seismic depth and/or time section "kong-macro" (versions) |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2010122578A RU2010122578A (en) | 2011-12-10 |
| RU2445651C2 true RU2445651C2 (en) | 2012-03-20 |
Family
ID=45405220
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2010122578/28A RU2445651C2 (en) | 2010-06-02 | 2010-06-02 | Method to build seismic depth and/or time section "kong-macro" (versions) |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2445651C2 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2559303C1 (en) * | 2014-02-06 | 2015-08-10 | Федеральное государственное унитарное предприятие "Нижне-Волжский научно-исследовательский институт геологии и геофизики" (ФГУП "НВНИИГГ") | Seismic prospecting method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4534019A (en) * | 1981-12-29 | 1985-08-06 | Mobil Oil Corporation | Common-depth-point method for determining and displaying the shear-velocity reflectivities of a geologic formation |
| RU2324205C1 (en) * | 2006-09-27 | 2008-05-10 | Федеральное государственное унитарное предприятие Сибирский научно-исследовательский институт геологии, геофизики и минерального сырья | Seismic data processing technique |
-
2010
- 2010-06-02 RU RU2010122578/28A patent/RU2445651C2/en active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4534019A (en) * | 1981-12-29 | 1985-08-06 | Mobil Oil Corporation | Common-depth-point method for determining and displaying the shear-velocity reflectivities of a geologic formation |
| RU2324205C1 (en) * | 2006-09-27 | 2008-05-10 | Федеральное государственное унитарное предприятие Сибирский научно-исследовательский институт геологии, геофизики и минерального сырья | Seismic data processing technique |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2559303C1 (en) * | 2014-02-06 | 2015-08-10 | Федеральное государственное унитарное предприятие "Нижне-Волжский научно-исследовательский институт геологии и геофизики" (ФГУП "НВНИИГГ") | Seismic prospecting method |
Also Published As
| Publication number | Publication date |
|---|---|
| RU2010122578A (en) | 2011-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102369752B1 (en) | Detection of underground structures | |
| Hobiger et al. | Site characterization of Swiss strong‐motion stations: The benefit of advanced processing algorithms | |
| EA026658B1 (en) | Extracting sv shear data from p-wave marine data | |
| RU2449322C1 (en) | Method of constructing seismic depth profile | |
| Louie et al. | Subsurface imaging of the Garlock fault, Cantil Valley, California | |
| RU2415449C1 (en) | Building method of seismic depth section | |
| CN105093314A (en) | Method for measuring and determining micro-seismic focus | |
| Hrubcová et al. | Shallow crustal discontinuities inferred from waveforms of microearthquakes: Method and application to KTB Drill Site and West Bohemia Swarm Area | |
| RU2445651C2 (en) | Method to build seismic depth and/or time section "kong-macro" (versions) | |
| RU2692001C1 (en) | Method of determining inclination angles of reflecting boundaries in a geological environment from profile seismic survey data 2d | |
| RU2454682C1 (en) | Method of constructing seismic dynamic and/or total time section | |
| Schmedes et al. | Imaging the shallow crust using teleseismic tomography | |
| RU2463628C1 (en) | Method of constructing seismic depth profile | |
| Ebel | The effect of crustal scattering on observed high-frequency earthquake seismograms | |
| US3354985A (en) | Seismic prospecting method of locating a suspected fault | |
| Bruno | High-resolution seismic imaging in complex environments: A comparison among common-reflection-surface stack, common-midpoint stack, and prestack depth migration at the Ilva-Bagnoli brownfield site, Campi Flegrei, Italy | |
| US2024921A (en) | Method of determining slope of subsurface rock beds | |
| RU2221262C1 (en) | Process of seismic prospecting for examination of sedimentary mantle in presence of heavily broken-up acoustically stiff boundaries ( variants ) | |
| RU2760889C1 (en) | Method for borehole seismic exploration | |
| Woelz et al. | Seismic prospecting in archaeology: a 3D shear‐wave study of the ancient harbour of Miletus (Turkey) | |
| RU2339978C1 (en) | Method of borehole sesmic exploration works | |
| RU2559303C1 (en) | Seismic prospecting method | |
| RU2148838C1 (en) | Seismic data processing technique | |
| Lin et al. | Three-dimensional P-wave velocity structure of the Bear Valley region of central California | |
| Louis et al. | Conventional and modern seismic investigations for rock quality determination at a dam site—a case history1 |