WO2015053659A1 - Способ получения априорного годографа для выполнения литолого стратиграфической привязки - Google Patents

Способ получения априорного годографа для выполнения литолого стратиграфической привязки Download PDF

Info

Publication number
WO2015053659A1
WO2015053659A1 PCT/RU2014/000482 RU2014000482W WO2015053659A1 WO 2015053659 A1 WO2015053659 A1 WO 2015053659A1 RU 2014000482 W RU2014000482 W RU 2014000482W WO 2015053659 A1 WO2015053659 A1 WO 2015053659A1
Authority
WO
WIPO (PCT)
Prior art keywords
hodograph
values
logging
point
interval
Prior art date
Application number
PCT/RU2014/000482
Other languages
English (en)
French (fr)
Inventor
Светлана Викторовна ШОЛОХОВА
Original Assignee
Светлана Викторовна ШОЛОХОВА
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Светлана Викторовна ШОЛОХОВА filed Critical Светлана Викторовна ШОЛОХОВА
Publication of WO2015053659A1 publication Critical patent/WO2015053659A1/ru

Links

Classifications

    • 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

Definitions

  • a method for producing an a priori hodograph for performing lithological-stratigraphic binding relates to the field of geophysics and can be used to obtain seismic sections of images of the geological environment through a prior hodograph and synthetic seismogram for geological exploration purposes.
  • Chtah the magnitude and sign of the adjusted values for V 0 and a 0 are determined by the differences:
  • indices i and (I) denote the subsequent and previous a priori and current initial values of the range of changes in speed and angle of incidence at the displayed point; sequentially in block F max Cmax, Chtah) replaced automatically previous range values and on the corrected and monotonically increasing velocity with depth values are replaced only when the coincidence scanning characters displayed point on columns with V velocity ⁇ ⁇ sign increments, thus placing F (q (j-1), ⁇ ( ⁇ )) in the center of a new grid-adjusted values of the matrix windows reliability survey for subsequent qx u incidence angles and speeds are corrected values F (qj, u j), among them within the next window VERVIEW credibility finding new current value of F max, another document showing the points with an intensity due to the found parameters F max, ° max and q max, F max is placed back into the center of a regular review of the adjusted mesh window reliability, and so on until the complete construction of the deep seismic section in accordance with the specified scan mode of the displayed points in rows and / or columns (See RF patent ⁇
  • V V 0 ⁇ w (1, 2, ⁇ ⁇ , ⁇ / 2, ⁇ ), where V 0 is the a priori initial value of velocity from a given array of velocities V (x, h), x, h are the coordinates of the deep section, ⁇ ⁇ is the step of the increment of speed, and is the number of step increments of speed; document all q u values of the energies of the oscillation trains
  • indices i and (i-1) denote the subsequent and previous a priori and current initial values of the range of the speed and the angle of inclination in the displayed point; in automatic mode, successively replace the previous values of the range ( ⁇ , V) with the corrected ones, thereby, W max (q (i-1) , ⁇ (i-1) ) is placed in the center of the new adjusted grid of values of the matrix of the reliability review window, again, for all q °, the corrected values of W (q i ( ⁇ ,) are found, among them within the confidence review window, find the current value W max , document the next displayed point with an intensity corresponding to the value of W ma , place W max again in the center of the next adjusted grid of the confidence review window and repeat the process of plotting the displayed points until the deep seismic section is completely constructed and in accordance with the specified scanning mode of the displayed points in rows and
  • a disadvantage of the known methods is the high error, the complexity and lack of reliability, high cost.
  • a seismic deep section See RF patent ⁇ ° 2463628, IPC G01V 1/28, published October 10, 2012, including a scan of the displayed points (OT) of section A (x, h) with a given column step (h ) and / or lines (x), from each OT (x 0 , h 0 ), a wave field survey is performed to search for the amplitudes of the signals A of seismic boundaries, possibly belonging to the reference and / or target horizons, in the n seismograms recorded by the multiple profiling method, which ensures illumination at various angles of the displayed points from n points Procedure, suggesting the presence in the displayed seismic boundary points; the review is implemented by a series of solutions of the quasidirect beam tracing problem, they find the most probable eikonal using a series of virtual hodographs to obtain marker timestamps on seismic trails; hodographs are generated in accordance with a predefined set of speeds
  • V 0 is the a priori initial velocity value from the given velocity array
  • V (x, h) x, h are the coordinates of the deep section, ⁇ ⁇ step of the increment of speed, ⁇ - the number of incremental increments of speed)
  • a predefined specific configuration of the seismic object which is defined by q angles at the displayed point with a certain step
  • V OI V 0 (H) ⁇ 5
  • V ⁇ / 0 (I) ⁇ L V (° / 2- ° max );
  • the indices i and (i-1) denote the subsequent and previous a priori and current initial values of the range of changes in speed and angle of incidence at the displayed point; sequentially in the block F max (° max , q max ), the previous values of the range "are automatically replaced by the corrected ones, and the values of speed monotonically increasing with depth are replaced only if the signs of the scan of the displayed point in the columns with the increment sign ⁇ ⁇ of the speed V coincide, place F (q (i-) , and ⁇ I) ) in the center of the new adjusted grid of values of the matrix of the confidence review window, for the subsequent qx u incidence angles and velocities find the corrected values F (qj, u j), among them, within the following windows about Zora credibility finding new current value of F max, another document showing the points with an intensity due to the found parameters F max, u max and q max, F max is placed back into the center once adjusted grid box reliability survey, and so on to complete the construction of the deep seismic section in accordance with
  • the disadvantage of this method is the low accuracy and reliability of the construction of a deep section, high labor costs, high cost.
  • the disadvantage of this method is the high complexity and low productivity (speed of obtaining the result), accuracy and reliability.
  • a method of obtaining an a priori hodograph for performing lithological-stratigraphic binding is to reduce the complexity and cost, increase the speed (speed of obtaining the result), as well as increase the reliability and accuracy of the correspondence of the horizons of the time section and geological marks of the well, as well as efficiency ratio of costs to results.
  • the technical result is achieved by the fact that in the known method of obtaining an a priori hodograph for performing lithologic-stratigraphic binding by obtaining seismic logging data, acoustic logging, density gamma-ray gamma logging, calibration of acoustic logging values, time-scale translation, calculation of reflection coefficients and their convolution with source signal, i.e.
  • the invention includes a series of sequential actions in which independently obtain and prepare data from common depth point methods, seismic logging, vertical seismic profiling, acoustic logging, density gamma-ray gamma-ray logging, geological marks and check the quality of these data, and also obtain the reference values of the interval velocities, after which the initial travel time curve is obtained and the synthetic seismogram is calculated, then quality control is carried out, then a permanent time correction is introduced for landing on the upper reference horizon lithological-stratigraphic complex, then the synthetic seismogram is calculated again, quality control is carried out again, after which the “zero point” is entered into the obtained year count and additionally one, or more points below * downhole calculation performed synthetic seismogram and quality control, after this correction is calculated and administered to seat against the lower support lithologic horizon stratigraphic complex of the formula
  • dTi is the correction introduced at each i-th point of the last hodograph following the point of the upper reference horizon
  • tH ogt is the time value taken from the time section of the “ MOS common depth point ” method (at the well point) for the extremum of the reflected wave corresponding to the lower reference horizon;
  • H n is the value of the geological elevation of the lower reference horizon according to the diagrams of geophysical research of wells; after that, the synthetic seismogram is re-calculated with subsequent quality control, and then the points of the obtained hodograph are transferred to the nearest acoustically weak boundaries, then the synthetic seismogram with subsequent quality control is re-calculated and the a priori hodograph is obtained, it is saved, the synthetic seismogram and the acoustic stiffness curve, while as a quality control of binding, use the cross-correlation function, which is a dynamic criterion for evaluating vibrational processes, being the kinematic criterion, the reference values of the interval velocities of this interval, obtained earlier according to seismic logging data, vertical seismic profiling, in 2D or 3D coordinates using the depth values of the reference and target horizons corresponding to the recording times of the first arrivals of the transmitted wave, also use the function to control the quality of the binding cross-correlation, which is
  • Improving the accuracy and reliability of the results in the claimed method namely improving the accuracy of the correspondence of the values of the times of the extrema of the reflections from the reference and target horizons of the time section to the values of the times of the geological marks of the well on a time scale, also determines the time correction dTi, obtains the reference values of the interval speeds and quality control and preparation MOS methods ", CK, GSP, AC, GGC-P, geology tap values and the introduction of" zero "in the polar plot, landing on the upper support horizon lit.
  • the correspondence, the degree of similarity of the extremes of reflections from the reference and target horizons of the synthetic seismogram to the extremes of reflections from the reference and target the horizons of the MOS section " is achieved by using a temporary MOS section " while maintaining relative amplitudes, leading it to a nulfase pulse; using, in addition to the dynamic criterion for evaluating vibrational processes, the cross-correlation function (CVF), the kinematic criterion, the values of the interval velocity.
  • CVF cross-correlation function
  • a method for obtaining an a priori hodograph for performing lithological-stratigraphic binding is that the final hodograph can be used to recalculate the time section into a deep section, and a synthetic seismogram, an acoustic impedance curve can be used to perform inversions, or when using several such hodographs and curves from different wells, the area of the same results (deep section and inversion) in the coordinates ZD.
  • the inventive "Method of obtaining an a priori hodograph for performing lithological-stratigraphic binding" also allows you to justify geological elevations in areas of complex structure (for example, in areas with anomalous structure of the Bazhenov formation), because allows you to control the reference values of the interval velocities at the point of the well, obtained in 2D or ZD coordinates and characteristic of the corresponding interval and a certain geological area, and control can also be carried out, for example, by comparing with the range of variation of the values of the interval velocities obtained in the laboratory by measurement in the corresponding interval on rock samples in the wells of a given geological area, interval intervals are compared with averaged values x speeds previously obtained from acoustic logging data for a given interval.
  • the set of essential features of the claimed invention “A method for obtaining an a priori hodograph for performing lithologic-stratigraphic binding” is not known from the prior art, which allows us to conclude that the claimed invention “A method for obtaining an a priori hodograph for performing lithologic-stratigraphic binding” the criterion of "novelty.”
  • Technological methods embodying the claimed invention “A method for obtaining an a priori hodograph for performing lithological-stratigraphic binding” during its implementation, can ensure the achievement of the technical result perceived by the applicant, namely, reducing the complexity, improving the accuracy and reliability of operational and technical qualities, therefore, the claimed invention “Method of obtaining an a priori hodograph for performing lithological-stratigraphic binding ”meets the criterion of“ industrial enimost ".
  • a method for obtaining an a priori hodograph for performing lithological and stratigraphic bindings can be repeatedly used in the process of obtaining a priori hodographs for performing lithological and stratigraphic bindings with obtaining a technical result consisting in reducing labor input, improving accuracy, reliability, operational and technical qualities, which allows a cost-effective and fast obtaining a priori hodograph, constructing a synthetic seismogram and acoustic stiffness curve to perform lithologic-stratigraphic binding, linking seismic and geological boundaries, obtain elements of an a priori velocity model for inverting seismic data into quantitative characteristics of the reservoir, i.e., increase the reliability of integrated geological analysis in order to create geological and technological models three-dimensional digital geological hydrodynamic models of oil and gas fields.
  • FIG. 1 shows a scheme for obtaining an a priori hodograph for performing lithological-stratigraphic binding.
  • a method of obtaining an a priori hodograph for performing lithological-stratigraphic binding was performed as follows.
  • acoustic logging was carried out in the well, in which lithological-stratigraphic binding was performed.
  • dt A K f (H) - the interval travel time of the longitudinal wave.
  • Density gamma-gamma-ray logging (GGK-P) was performed, and a diagram of changes in scattered gamma radiation as a function of depth was recorded. Checked quality of this data. Prepared a chart of the interval time of acoustic logging (AK) with:
  • V MHT _CK_ ( BA) 2dH / dT, where dH is the interval power (in meters) between the values of the geological elevations corresponding to the reference exhaust gas “B” and “A” according to the well logs; dT is the temporary power of the interval (in seconds) between the reference exhaust gas “B” and “A”, which served as criteria for the accuracy of the binding.
  • dTi ⁇
  • t B orr is the time value taken from the time section MOS " (at the point of the well) for the extremum of the reflected wave corresponding to the reference exhaust gas" B ";
  • ecK is the time value taken from the last hodograph of the SK for the geological mark of the reference exhaust gas “B”;
  • N B the value of the geological elevation of the reference OG "B" according to the GIS diagrams
  • the final synthetic seismogram was calculated using the corrected hodograph and quality control of the hodograph and synthetic seismogram was carried out using FVC in the corresponding interval and compared with the reference values of the interval velocities. Then saved a priori hodograph, synthetic seismogram, acoustic stiffness curve.
  • dH is the distance in m from the value of the geological elevation of the reference exhaust gas “B” to the value of the depth of the bottomhole point;
  • T B time value taken from the time section of the MOGT and corresponding to the time of the extremum of the reflected wave from the reference exhaust gas “B”;
  • VMHT (B-A) is the reference interval velocity according to SC and VSP obtained in clause 5 for the lowest lithologic-stratigraphic complex, for example, the interval between GO "B” and GO "A"
  • the proposed "Method for obtaining an a priori hodograph for performing lithological-stratigraphic binding" reduces the complexity and cost, increases the speed of obtaining the result, and also increases the accuracy and reliability of the correspondence of the time values of the extremes of reflections from the reference and target horizons of the synthetic seismogram to the time values of the extremes of reflections from the reference and target horizons MOS " time cut and efficiency as a ratio of cost to outcome.
  • dTi is the correction introduced at each i-th point of the last hodograph following the point of the upper reference horizon;
  • t B orr is the time value taken from the MOS time section " (at the well point) for the extremum of the reflected wave corresponding to horizon" B ";
  • N B the value of the geological elevation of the reflecting horizon "B" according to the GIS diagrams
  • a hodograph point corresponding to OG “B” (N B - the value of the geological elevation of OG “B”, T B - the time value of the extremum of the reflected wave of OG “B”, taken from the time section MOS " ).
  • dH is the distance in m from the value of the geological elevation of OG "B" to the value of the depth of the bottomhole point;
  • T B the value of time taken from the time section MOS " and corresponding to the extremum from the OG" B ";
  • the time values were also calculated for the two last hodograph points entered using the reference interval velocity of this interval.
  • the proposed "Method for obtaining an a priori hodograph for performing lithological-stratigraphic binding" reduces the complexity and cost, increases the speed of obtaining the result, and also increases the accuracy and reliability of the correspondence of the time values of the geological marks of the reference and target exhaust gases in the time scale to the time values extrema respective reference and target horizons time section PMOS "and cost efficiency as a ratio to the result.

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

Способ получения априорного годографа для выполнения литолого-стратиграфической привязки относится геофизике для получения сейсмических разрезов изображений геологической среды. Для снижения трудоемкости и повышения достоверности соответствия горизонтов временного разреза и ГО скважины он включает действия, дающих данные методов общей глубинной точки, СК, ВСП, АК, ГГК-П, ГО и проверяющие их качество, дающие эталонные значения интервальных скоростей, получают исходный годограф и рассчитывают синтетическую сейсмограмму, проводят контроль качества, и вводят постоянную временную поправку для посадки на верхний ОГ литолого-стратиграфического комплекса и вновь рассчитывают синтетическую сейсмограмму и проводят контроль качества, вводят «нулевую точку» в полученный годограф, дополнительно, на одну и более точек ниже забоя, проводят расчет синтетической сейсмограммы и контроль качества, вновь рассчитывают и вводят поправку для посадки на нижний ОГ литолого-стратиграфического комплекса по формуле dTi=Hi*dTн/Hн, где dTi - поправка, вводимая в каждую i-тую точку последнего годографа, следующую за точкой верхнего ОГ; dTн=tнОГТ-tнск - поправка, вносимая в последний годограф для посадки на нижний ОГ, tногт - значение времени снятое с временного разреза МОГТ (в точке скважины) для экстремума отраженной волны соответствующей нижнему ОГ; tнск - значение времени, снятое с последнего годографа СК, ВСП для геологической отметки нижнего ОГ; Hi - значение геологической отметки i-ой точки последнего годографа СК, ВСП; Нн - значение геологической отметки нижнего ОГ по диаграммам ГИС; после этого вновь рассчитывают синтетическую сейсмограмму с последующим контролем качества, и следом переносят точки полученного годографа на ближайшие акустически слабые границы и повторно рассчитывают синтетическую сейсмограмму, получая априорный годограф.

Description

СПОСОБ ПОЛУЧЕНИЯ АПРИОРНОГО ГОДОГРАФА ДЛЯ
ВЫПОЛНЕНИЯ ЛИТОЛОГО СТРАТИГРАФИЧЕСКОЙ ПРИВЯЗКИ Заявляемое изобретение «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» относится к области геофизики и может быть использовано для получения сейсмических разрезов изображений геологической среды через априорный годограф и синтетическую сейсмограмму в геологоразведочных целях.
Известен способ построения сейсмического глубинного разреза, включающий развертку отображаемых точек (ОТ) разреза A(x,h) с заданным шагом по столбцам (h) и/или строкам (х), из каждой ОТ (х0 ,h0) производят обзор волнового поля для поиска амплитуд сигналов А сейсмических границ, возможно принадлежащих реперным и/или целевым горизонтам, на п сейсмограммах, зарегистрированных способом многократного профилирования, обеспечивающего засветку под различными углами отображаемых точек из п пунктов воздействий, предполагая наличие сейсмической границы в отображаемых точках; обзор реализуют серией решений квазипрямой задачи лучевого трассирования, находят наиболее вероятный эйконал, используя серию виртуальных годографов для получения маркерных отметок времени на сейсмотрассах; генерацию годографов производят в соответствии с предварительно заданным множеством скоростей ν=ν0±Δ V(1 , 2," ", υ/2, '", υ), (где V0 - априорно заданное первоначальное значение скорости из заданного массива скоростей V(x,h), х, h - координаты глубинного разреза, Δν шаг приращения скорости, и- число шаговых приращений скорости), а также с предварительно заданной определенной конфигурацией сейсмического объекта, который задан q ракурсами в отображаемой точке с
ЛГУ /у
определенным шагом поворота угла падения с центром вращения в отображаемой точке, для каждого ракурса объекта задают множество углов падения а=±а 0±Δ α(1 , 2, , q/2, , q),
(где и 0 - априорно заданное первоначальное значение угла падения границы в отображаемой точке, Δθί- шаг приращения угла падения, q - число пошаговых приращений угла падения), задают координаты отображаемой точки, п пунктов воздействий, m сейсмоприемников и закон (функцию) расположения последовательности из m точек, прилегающих к отображаемой точке; для каждого заданного типа волны строят qu mn виртуальных годографов, производят направленное псевдосинфазное накопление амплитуд многокомпонентных цугов колебаний в суммоцугах на зарегистрированных m сейсмотрассах вдоль маркерных отметок времени на годографах, причем для каждого многокомпонентного цуга колебаний считывают ркомпонентов (обычно р~20) в виде функционалов FP ; после накопления в каждом из суммоцугов в виде функционалов FP m по рт сигналов их группируют по признаку равенства номера ракурса q; после введения известными способами поправок на углы засветки и мьютинга, причем в соответствии с интерактивным режимом интерпретации суммоцугов отдельной сейсмограммы и коррекции заданных диапазонов V, ии и, q для каждой из υ скоростей и каждого из q суммоцугов известным способом после вычисления qu функционалов статистической обработки FP m строят для каждой ОТ двумерную матрицу из элементов их энергии Fw Рт * ,V), после соответствующих вычислений в зависимости от априорных углов падения аи априорно заданного множества скоростей V по столбцам матрицы документируют значения функционалов в зависимости от угла падения в отображаемой точке, а по строкам - значения функционала в зависимости от скорости (или наоборот); строят отображение графика из матрицы Fw Рт(* ,V), на котором визуализируют значения функционалов в виде изолиний, гистограмм или иных известных способов; график представляют в виде окна обзора достоверности и обеспечивают контроль интерпретатора за процессом опознания сейсмических границ и построением сейсмического глубинного разреза; далее уровень Fw Рт(° ,V), в пределах окна обзора достоверности сравнивают с заданным порогом обнаружения, при превышении порога обнаружения каждый последующий уровень отображаемого функционала сравнивают с предыдущим, при превышении последующего уровня над предыдущим его документируют в блоке максимума функционала статистической обработки
Fmax(°max, Qmax), в котором также хранят соответствующее максимальному уровню аргументы и тах, и qmax, изображают отображаемую точку на глубинном разрезе с интенсивностью, обусловленной найденными параметрами Fmax, u max и qmax в пределах окна обзора достоверности, переходят к документированию последующей отображаемой точки на глубинном разрезе, присваивают Fmax положение F(u /2,q/2) в центре матрицы окна обзора достоверности и определяют скорректированные значения углов падения и скоростей по смещению Fmax относительно центра окна обзора достоверности, для чего сравнивают аргументы υ/2 и q/2 с и тах и
Чтах! величину и знак скорректированных значений для V0 и а 0 определяют по разностям:
Voi =ν ν=ν0(Μ)±Δ V(° /2-и max);
аш = °4- ± δα = аои) ± Aa(q/ 2~ O
где индексами i и (И ) обозначены последующие и предыдущие априорные и текущие начальные значения диапазона изменения скорости и угла падения в отображаемой точке; последовательно в блоке FmaxCmax, Чтах) заменяют в автоматическом режиме предыдущие значения диапазона ана скорректированные, а монотонно возрастающие с глубиной значения скорости заменяют лишь при условии совпадения знаков развертки отображаемой точки по столбцам со знаком приращения δν скорости V, тем самым размещают F(q(j-1), υ (μΐ)) в центр новой скорректированной сетки значений матрицы окна обзора достоверности, для последующих qxu углов падения и скоростей находят скорректированные значения F(qj,u j), среди них в пределах следующего окна обзора достоверности находят новое текущее значение Fmax, документируют очередную отображаемую точку с интенсивностью, обусловленной найденными параметрами Fmax, °max и qmax, размещают Fmax вновь в центр очередной скорректированной сетки окна обзора достоверности, и так далее до полного построения глубинного сейсмического разреза в соответствии с заданным режимом развертки отображаемых точек по строкам и/или столбцам (См. патент РФ Ν°2463626, ΜΠΚ ΘΟΐν 1/28, опубл. 10.10.2012).
Известен также (См. патент Ν°2449322, МПК G01V 1/28, опубл.27.04.12) способ построения сейсмического глубинного разреза, включающий развертку отображаемых точек разреза по столбцам (h) и/или строкам (х), из каждой отображаемой точки ОТ (х0, h0 ) производят обзор волнового поля на сейсмограммах для поиска амплитуд сигналов А сейсмических границ, возможно принадлежащих реперным и/или целевым горизонтам; обзор реализуют серией решений квазипрямой задачи лучевого трассирования с использованием виртуальных годографов, для чего задают определенную конфигурацию сейсмического объекта и при каждом облучении отображаемой точки из пункта воздействия задают q ракурсов сейсмического объекта с
ΛίΥ ίΥ
определенным шагом поворота угла наклона с центром вращения в отображаемой точке, для каждого ракурса объекта а=±а 0±Δθί О , 2, q/2, , q), где α 0 - априорно заданное первоначальное значение угла наклона границы в отображаемой точке, рассчитывают виртуальный годограф заданного типа сейсмической волны с учетом априорно заданной скорости, координат отображаемой точки и заданного закона (функции) расположения последовательности из m точек, прилегающих к отображаемой точке, производят направленное синфазное накопление амплитуд многокомпонентных колебаний на зарегистрированных m сейсмотрассах вдоль маркерных отметок времени на годографах, причем для каждого колебания считывают P компонентов (P<20) в виде псевдосинфазных отсчетов величин амплитуд; после синфазного накопления вдоль q годографов по т амплитуд колебаний сигналов из каждой из п сейсмограмм сигналы группируют по признаку равенства номера шага позиционирования (ракурса) q, в каждой из q групп производят супернакопления из Pmn амплитуд сигналов после введения поправочных коэффициентов на угол засветки, для каждого из q супернакоплений амплитуд по Рфаз каждого из mn сигналов, специфичных для заданных позиций объекта в отображаемой точке, вычисляют q функционалов статистической обработки, сравнивают их между собой, по максимальной величине одного из функционалов определяют наиболее вероятный номер q позиции объекта, причём после определения наиболее вероятного номера q позиции объекта для каждой отображаемой точки строят двумерную сетку значений функционалов статистической обработки типа энергий супернакоплений W(w , V) в зависимости от априорных углов наклона α и априорно заданного множества скоростей V при
V=V 0±w (1 , 2, · ·, υ/2, υ), где V0 - априорно заданное первоначальное значение скорости из заданного массива скоростей V(x, h), х, h - координаты глубинного разреза, Δν - шаг приращения скорости, и- число шаговых приращений скорости; документируют все q u значений энергий цугов колебаний
, V) в узлах прямоугольной сетки в виде двумерной матрицы, по столбцам которой документируют значения энергии в зависимости от угла наклона в отображаемой точке, а по строкам - значения энергии в зависимости от скорости (или наоборот); строят двумерное отображение графика W(a , V), график представляют в виде окна обзора достоверности, что обеспечивает визуальный контроль за процессом опознания сейсмических границ и построением сейсмического глубинного разреза; в пределах графика окна обзора достоверности определяют максимум Wmax (° max, qmax), его сравнивают с заданным порогом обнаружения, при превышении порога обнаружения отображаемую точку наносят на глубинный сейсмический разрез с интенсивностью, пропорциональной величине максимального функционала в пределах окна обзора достоверности; а затем для построения последующей отображаемой точки на глубинном разрезе определяют скорректированные значения углов наклона и скоростей по смещению Wmax относительно центра окна обзора достоверности, где находят W(u 12, q/2) и сравнивают его аргументы υ/2 и q/2 с °max и qmax; величину и знак скорректированных значений для V0 и к 0 определяют по разностям:
Figure imgf000009_0001
0i= 0(i-1 )± - 0(i-1)± \ / -итах), где индексами i и (i-1 ) обозначены последующие и предыдущие априорные и текущие начальные значения диапазона изменения скорости и угла наклона в отображаемой точке; в автоматическом режиме последовательно заменяют предыдущие значения диапазона (α , V) на скорректированные, тем самым размещают Wmax (q(i-1), υ (i-1)) в центр новой скорректированной сетки значений матрицы окна обзора достоверности, снова для всех q ° находят скорректированные значения W(qi( υ,), среди них в пределах окна обзора достоверности находят текущее значение Wmax, документируют очередную отображаемую точку с интенсивностью, соответствующей величине Wma , размещают Wmax вновь в центр очередной скорректированной сетки окна обзора достоверности и повторяют процесс нанесения отображаемых точек до полного построения глубинного сейсмического разреза в соответствии с заданным режимом развертки отображаемых точек по строкам и/или столбцам
Недостатком известных способов является - высокая погрешность, трудоемкость и недостаточная достоверность, высокая себестоимость.
Поскольку достоверность проявляется в улучшении динамических характеристик разреза, т.к. интенсивность изображения точки обусловлена найденными параметрами, поэтому данные недостатки обусловлены тем, что в известных технических решениях увеличение степени достоверности построения разреза за счет учета влияния сейсмического сноса полностью относится к положению объекта в плане (координатах «X» и «У»), чего нельзя сказать о вертикальной составляющей глубины «Н», поскольку используемые виртуальные годографы, используют модельные данные скоростей, а не реальные скорости. Кроме того, контроль параметров построения осуществляется интерпретатором и интерактивно, т.е. с повторением циклов, что требует дополнительных затрат ресурсов, причем сама методика работ МОП", которая предусматривает получение и обработку данных на поверхности земли, без использования данных скважины. При реализации известных способов возможно получение нескольких вариантов разреза, что делает решение задачи неоднозначным.
Известен также способ построения сейсмического глубинного разреза (См. патент РФ Ν°2463628, МПК G01V 1/28, опубл.10.10.2012), включающий развертку отображаемых точек (ОТ) разреза A(x,h) с заданным шагом по столбцам (h) и/или строкам (х), из каждой ОТ (х0 ,h0) производят обзор волнового поля для поиска амплитуд сигналов А сейсмических границ, возможно принадлежащих реперным и/или целевым горизонтам, на п сейсмограммах, зарегистрированных способом многократного профилирования, обеспечивающего засветку под различными углами отображаемых точек из п пунктов воздействий, предполагая наличие сейсмической границы в отображаемых точках; обзор реализуют серией решений квазипрямой задачи лучевого трассирования, находят наиболее вероятный эйконал, используя серию виртуальных годографов для получения маркерных отметок времени на сейсмотрассах; генерацию годографов производят в соответствии с предварительно заданным множеством скоростей
ν=ν0±Δ V(1 , 2, , υ/2, -, υ), (где V0 - априорно заданное первоначальное значение скорости из заданного массива скоростей; V(x,h), х, h - координаты глубинного разреза, Δν шаг приращения скорости, υ- число шаговых приращений скорости), а также с предварительно заданной определенной конфигурацией сейсмического объекта, который задан q ракурсами в отображаемой точке с определенным шагом
ΛίΥ of
поворота угла падения с центром вращения в отображаемой точке, для каждого ракурса объекта задают множество углов падения 01 = ±а 0±А "(1 , 2, , q/2,— , q), (где α 0 - априорно заданное первоначальное значение угла падения границы в отображаемой точке, - шаг приращения угла падения, q - число пошаговых приращений угла падения), задают координаты отображаемой точки, п пунктов воздействий, m сейсмоприемников и закон (функцию) расположения последовательности из m точек, прилегающих к отображаемой точке; для каждого заданного типа волны строят q° mn виртуальных годографов, производят направленное псевдосинфазное накопление амплитуд многокомпонентных цугов колебаний в суммоцугах на зарегистрированных m сейсмотрассах вдоль маркерных отметок времени на годографах, причем для каждого многокомпонентного цуга колебаний считывают компонентов (обычно р~20) в виде функционалов FP ; после накопления в каждом из суммоцугов в виде функционалов FP m по рт сигналов их группируют по признаку равенства номера ракурса q; после введения известными способами поправок на углы засветки и мьютинга, причем в соответствии с интерактивным режимом интерпретации суммоцугов отдельной сейсмограммы и коррекции заданных диапазонов V, ии и, q для каждой из υ скоростей и каждого из q суммоцугов известным способом после вычисления qu функционалов статистической обработки FP m строят для каждой ОТ двумерную матрицу из элементов их энергии Fw Рт(* ,V), после соответствующих вычислений в зависимости от априорных углов падения а и априорно заданного множества скоростей V по столбцам матрицы документируют значения функционалов в зависимости от угла падения в отображаемой точке, а по строкам - значения функционала в зависимости от скорости (или наоборот); строят отображение графика из матрицы Fw Рт(* ,V), на котором визуализируют значения функционалов в виде изолиний, гистограмм или иных известных способов; график представляют в виде окна обзора достоверности и обеспечивают контроль интерпретатора за процессом опознания сейсмических границ и построением сейсмического глубинного разреза; далее уровень Fw Рт(* ,V), в пределах окна обзора достоверности сравнивают с заданным порогом обнаружения, при превышении порога обнаружения каждый последующий уровень отображаемого функционала сравнивают с предыдущим, при превышении последующего уровня над предыдущим его документируют в блоке максимума функционала статистической обработки
FmaX(° max.qmax). В KOTOpOM ТЭЮКв ХрЭНЯТ СООТВвТСТВуЮЩее максимальному уровню аргументы и тах, и qmax, изображают отображаемую точку на глубинном разрезе с интенсивностью, обусловленной найденными параметрами Fmax, и тах и qmax в пределах окна обзора достоверности, переходят к документированию последующей отображаемой точки на глубинном разрезе, присваивают Fmax положение F(° /2, q/2) в центре матрицы окна обзора достоверности и определяют скорректированные значения углов падения и скоростей по смещению Fmax относительно центра окна обзора достоверности, для чего сравнивают аргументы υ/2 и q/2 с и тах и тах! величину и знак скорректированных значений для V0 и а 0 определяют по разностям: VOI =V0(H )±5 V= \/0(И)±Л V(° /2-° max); αο» = ao(i-i) ± δα = α ι) ± Aa(q/ 2- Чшах );
где индексами i и (i-1 ) обозначены последующие и предыдущие априорные и текущие начальные значения диапазона изменения скорости и угла падения в отображаемой точке; последовательно в блоке Fmax(° max, qmax) заменяют в автоматическом режиме предыдущие значения диапазона "на скорректированные, а монотонно возрастающие с глубиной значения скорости заменяют лишь при условии совпадения знаков развертки отображаемой точки по столбцам со знаком приращения δν скорости V, тем самым размещают F(q(i- ), и {И)) в центр новой скорректированной сетки значений матрицы окна обзора достоверности, для последующих qxu углов падения и скоростей находят скорректированные значения F(qj,u j), среди них в пределах следующего окна обзора достоверности находят новое текущее значение Fmax, документируют очередную отображаемую точку с интенсивностью, обусловленной найденными параметрами Fmax, u max и qmax, размещают Fmax вновь в центр очередной скорректированной сетки окна обзора достоверности, и так далее до полного построения глубинного сейсмического разреза в соответствии с заданным режимом развертки отображаемых точек по строкам и/или столбцам.
Недостатком известного способа является низкая точность и достоверность построения глубинного разреза, большие трудозатраты, высокая себестоимость.
Данные недостатки обусловлены тем, что достоверность построения глубинного разреза за счет учета влияния сейсмического сноса полностью относится к положению объекта в плане (координатах «X» и «У»), чего нельзя сказать о вертикальной составляющей глубины «Н», т.к. используются виртуальные годографы, модельные данные скоростей, а не реальные скорости среды. Критерии оценки используют только динамическую характеристику колебательных процессов - матрицы энергии функционалов Fw Рт(* ,V), а кинематическая характеристика в качестве критерия не используется. Кроме всего этого, массив «т» и «п» данных требует немало места в оперативной памяти ЭВМ, поэтому обработка «п» сейсмограмм по «т» сейсмотрасс требует значительных затрат машинного времени, причем контроль интерпретатора за процессом опознания сейсмических границ и построением сейсмического глубинного разреза требует выезд интерпретатора в поле и больших ресурсных затрат. Следует заметить, что метод, аналогичный методу регулированного направленного приема (МРНП) связан с хорошим приемом высокочастотных сейсмических волн не дает хорошего результата, поскольку теряет низкочастотную часть спектра, что, в итоге, сказывается на разрешающей способности и обусловливает ограничение известного способа. Данный недостаток, обусловлен методикой производства работ метода регулированного направленного приема (МРНП), который является аналогом описанного способа. В МРНП "волна, у которой годограф не совпадает с линией суммирования, попадает в область подавления характеристики направленности. Это проявляется тем сильнее, чем выше ее частотный спектр. Поэтому МРНП по отношению к высокочастотным сейсмическим волнам обладает более высокой разрешающей способностью ". Бондарев В. И., Крылатков СМ. Анализ данных сейсморазведки: Учебное пособие для студентов вузов. Екатеринбург. Издательство УГГГА, 2002. - 57с.
Известен принятый за прототип способ построения синтетической сейсмограммы для одномерного моделирования (см. Гогоненков Г.Н. «Расчет и применение синтетических сейсмограмм». М., Недра, 1972, стр.46-49, 77-85), при котором «калибруют» значения акустической диаграммы метода акустического каротажа (ΤΑκ=ΐ(ΗΑκ)) по исходному годографу сейсмокаротажа (СК) или вертикального сейсмического профилирования (ВСП) (Тск_исх= ^(НСк_исх))- Затем рассчитывают последовательность коэффициентов отражения (К0ТР ) по данным метода акустического каротажа (АК). После этого производят расчет импульсной синтетической сейсмограммы, полученной посредством свертки последовательности коэффициентов отражения с исходным сигналом. Далее последовательно каждой точке глубины «Н» исходного годографа сейсмокаротажа или вертикального сейсмического профилирования находят соответствующее значение «Т» временного разреза МОП" (ищут новую зависимость ТСК_ПРИВ = f(HcK_nPHB))- Контроль качества одномерного моделирования (привязки) осуществляют с помощью функции взаимной корреляции ФВК ближайшей к скважине трассы ОГТ временного разреза и рассчитанной синтетической сейсмограммы, добиваясь максимального подобия экстремумов отражений от опорных и целевых горизонтов синтетической сейсмограммы экстремумам отражений от соответствующих опорных и целевых горизонтов временного разреза МОП".
Недостатком известного способа является высокая трудоемкость и низкие производительность (скорость получения результата), точность и достоверность.
Данный недостаток обусловлен недостаточной достоверностью, точностью, степенью подобия, соответствия экстремумов от реперных и целевых геологических горизонтов синтетической сейсмограммы экстремумам соответствующих конкретных отраженных волн временного разреза, причем в качестве критерия оценки качества используется только динамический критерий колебательных процессов - функция взаимной корреляции (ФВК). Не используются кинематические критерии - скорости, что приводит к неточности выполнения привязки, необоснованной «растяжке» или «сжатии» годографа, привязке к не той фазе временного разреза, образованию неоправданных «ступеней» на графике зависимости интервальной скорости акустического каротажа от глубины, чего не бывает на качественных исходных диаграммах акустического каротажа, особенно, если привязка выполняется только в целевых интервалах, поскольку интервальная скорость с глубиной меняется по определенному закону. Кроме того, при выполнении привязки в целевом интервале, при отсутствии контроля диапазона целевого интервала опорными горизонтами, не всегда есть возможность правильно соотнести значения скорости привязочного годографа VMHT-IIPHB И тренда скорости VCK по данным СК и ВСП в точке скважины. Следует также заметить, что в процессе получения априорного годографа для выполнения литолого-стратиграфической привязки идет последовательный подбор значений точек годографа (глубина- время), причем в, годографах ВСП количество точек годографа может достигать 200 и более, последовательный подбор значений в таком случае занимает много времени, т.к. очень трудоемок.
Техническим результатом заявляемого изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» является снижение трудоемкости и себестоимости, увеличение скорости (скорости получения результата), а также повышение достоверности и точности соответствия горизонтов временного разреза и геологических отметок скважины, а также эффективности, как отношение затрат к результату.
Поставленный технический результат достигается тем, что в известном способе получения априорного годографа для выполнения литолого-стратиграфической привязки путем получения данных сейсмокаротажа, акустического каротажа, плотностного гамма-гамма каротажа, калибровки значений акустического каротажа, перевода во временной масштаб, расчета коээфициентов отражения и свертки их с исходным сигналом, т.е. построения синтетической сейсмограммы и кривой акустической жесткости с повторением циклов, согласно изобретения, он включает в себя ряд последовательных действий, при которых независимо получают и подготавливают данные методов общей глубинной точки, сейсмокаротажа, вертикального сейсмического профилирования, акустического каротажа, плотностного гамма-гамма каротажа, геологических отметок и проверяют качество этих данных, а также получают эталонные значения интервальных скоростей, после чего получают исходный годограф и рассчитывают синтетическую сейсмограмму, затем проводят контроль качества, далее вводят постоянную временную поправку для посадки на верхний опорный горизонт литолого-стратиграфического комплекса, затем вновь рассчитывают синтетическую сейсмограмму, вновь проводят контроль качества, после чего вводят «нулевую точку» в полученный годограф и дополнительно одну, и более точки* ниже забоя скважины, проводят расчет синтетической сейсмограммы и контроль качества, вслед за этим рассчитывают и вводят поправку для посадки на нижний опорный горизонт литолого-стратиграфического комплекса, по формуле
dTi= Н|*сП"н/Нн где
dTi - поправка, вводимая в каждую i-тую точку последнего годографа, следующую за точкой верхнего опорного горизонта; dTH - = tH - tH - поправка, которую нужно ввести в последний годограф для посадки на нижний опорный горизонт, т.е. совмещения во временном масштабе значения геологической отметки нижнего опорного горизонта и значения времени экстремума соответствующей отраженной волны;
tH огт - значение времени снятое с временного разреза метода общей глубинной точки МОП" (в точке скважины) для экстремума отраженной волны, соответствующей нижнему опорному горизонту;
tH ск- значение времени, снятое с последнего годографа сейсмокаротажа, вертикального сейсмического профилирования, для геологической отметки нижнего опорного горизонта; сейсмокаротажа, вертикального сейсмического профилирования;
Нн -значение геологической отметки нижнего опорного горизонта по диаграммам геофизического исследования скважин; после этого вновь рассчитывают синтетическую сейсмограмму с последующим контролем качества, и следом переносят точки полученного годографа на ближайшие акустически слабые границы, далее повторно рассчитывают синтетическую сейсмограмму с последующим контролем качества и получают априорный годограф, сохраняют его, синтетическую сейсмограмму и кривую акустической жесткости, при этом в качестве контроля качества привязки используют функцию взаимной корреляции, которая является динамическим критерием оценки колебательных процессов, и являющиеся кинематическим критерием эталонные значения интервальных скоростей данного интервала, полученные ранее по данным сейсмокаротажа, вертикального сейсмического профилирования, в координатах 2Д или ЗД с использованием значений глубин опорных и целевых горизонтов, соответствующих времен регистрации первых вступлений проходящей волны, в качестве контроля качества привязки используют также функцию взаимной корреляции, которая является динамическим критерием оценки колебательных процессов и являющиеся кинематическим критерием эталонные значения интервальных скоростей, которые получают лабораторным путем посредством измерения в соответствующем интервале на образцах горных пород в скважинах данного геологического района, в качестве контроля качества привязки используют также функцию взаимной корреляции, которая является динамическим критерием оценки колебательных процессов и являющиеся кинематическим критерием эталонные значения интервальных скоростей, которые получают осреднением значений интервальных скоростей, полученных ранее по данным акустического каротажа на данном интервале, в качестве исходного годографа берут также осредненный годограф, как результат аппроксимации годографов сейсмокаротажа, вертикального сейсмического профилирования данной площади, в качестве исходного годографа берут также годограф, значения которого сняты со значений глубин опорных и целевых горизонтов, соответствующих времен регистрации первых вступлений проходящей волны, используемых для получения эталонных значений интервальной скорости в координатах 2Д или ЗД, в качестве исходного годографа берут также годограф, значения которого сняты с тренда изменения диаграммы интервального времени акустического каротажа, как функции глубины, в качестве исходного сигнала берут также импульс выделенный в целевом интервале вблизи скважины с суммарного временного разреза метода общей глубинной точки, полученного с сохранением относительных амплитуд и приведенного к нульфазовому импульсу, в качестве исходного сигнала берут также нульфазовый модельный импульс подобранный по частоте к частотам суммарного временного разреза метода общей глубинной точки, полученного с сохранением относительных амплитуд и приведенного к нульфазовому импульсу, при этом в скважинах, с глубинами забоя в самом нижнем литолого-стратиграфическом комплексе, фиксируют точку годографа соответствующую значению экстремума опорного отражающего горизонта, являющегося верхней границей этого литолого-стратиграфического комплекса, после чего вводят точки соответствующих глинистых горизонтов данного комплекса и рассчитывают значение времени для точки глубины забоя скважины и добавленных ниже точек годографа, после чего повторяют расчет синтетической сейсмограммы и следом осуществляют контроль качества и получают априорный годограф, сохраняют его, синтетическую сейсмограмму и кривую акустической жесткости, а в качестве функции для перевода во временной масштаб используют годографы сейсмокаротажа, вертикального сейсмического профилирования без участия данных акустического каротажа и построения синтетической сейсмограммы.
Между отличительными признаками и достигаемым техническим результатом существует следующая причинно- следственная связь.
В отличие от аналога и прототипа использование в предлагаемом изобретении «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» совокупности признаков в виде того, что он включает в себя следующий ряд последовательных действий, как то, получение данных сейсмокаротажа, акустического каротажа, плотностного гамма-гамма каротажа, калибровки значений акустического каротажа, перевода во временной масштаб, расчета последовательности коэффициентов отражения, свертки их с исходным сигналом, т.е построения синтетической сейсмограммы и кривой акустической жесткости с повторением циклов, при которых независимо получают и проводят подготовку данных методов общей глубинной точки, сейсмокаротажа, вертикального сейсмического профилирования, акустического каротажа и плотностного гамма- гамма каротажа, значений геологических отметок и проводят проверку качества полученных выше данных, при этом получают глубины опорных и целевых сейсмических горизонтов, которые соответствуют времени регистрации первых вступлений проходящей волны от опорных и целевых горизонтов, и эталонные интервальные скорости по данным сейсмокаротажа и вертикального сейсмического профилирования между этими опорными и целевыми горизонтами в координатах 2Д или ЗД, затем определяют исходный годограф после чего рассчитывают синтетическую сейсмограмму, затем проводят контроль качества, вводят постоянную временную поправку для посадки на верхний опорный горизонт литолого-стратиграфического комплекса, например ОГ «Г» и рассчитывают синтетическую сейсмограмму и приводят контроль качества, затем вводят «нулевую точку» в полученный годограф и дополнительно одну, две точки ниже забоя скважины, проводят построение синтетической сейсмограммы и контроль качества вслед за этим вводят поправку для посадки на нижний опорный горизонт литолого-стратиграфического комплекса, например, ОГ «Б» и вновь рассчитывают синтетическую сейсмограмму, причем переносят точки полученного годографа на ближайшие акустически слабые границы и повторно рассчитывают синтетическую сейсмограмму. В совокупности признаков все это значительно повышает эффективность заявленного способа, поскольку дает возможность получать окончательный годограф намного быстрее, чем в объекте-прототипе, где использовали повторение циклов (подбором Н и V) для каждой точки годографа, поскольку в заявляемом способе обе предложенные поправки вводят в точки преобразуемого годографа одновременно, рассчитывается значение времени для точки забоя и последних точек годографа, что увеличивает скорость получения априорного годографа с использованием значений эталонных интервальных скоростей и повышает эффективность способа за счет использования готовых данных методов общей глубинной точки (МОП"), сейсмокаротажа (СК), вертикального сейсмического профилирования (ВСП), акустического каротажа (АК), гамма-гамма каротажа плотностного (ГГК-П), значений геологических отметок с исключением полевого этапа, кроме того, выполняется проверка качества и подготовка методов МОП", СК, ВСП, АК, ГГК-П, значений геологических отметок, а также возможность использования в качестве исходного годографа данных, используемых для получения эталонных значений интервальных скоростей, т.е. годографа, полученного по данным СК и ВСП в координатах 2Д или ЗД; осредненного годографа по площади; годографа по данным АК, что значительно снижает себестоимость заявленного способа с учетом временной поправки dTi и подготовки данных - редакции годографов СК или ВСП с учетом проверки их качества; получения суммарного временного разреза МОП" с сохранением относительных амплитуд, приведением его к нульфазовому импульсу; моделирование кривых АК и ГГК-П и возможности использования в качестве исходного годографа данных, используемых для получения эталонных значений интервальных скоростей, т.е. глубин опорных и целевых горизонтов, соответствующих времен регистрации проходящих волн, полученных по данным СК и ВСП в координатах 2Д или ЗД; осредненного годографа по площади; годографа по данным АК. Повышение точности и достоверности результатов в заявленном способе, а именно повышение точности соответствия значений времен экстремумов отражений от опорных и целевых горизонтов временного разреза значениям времен геологических отметок скважины во временном масштабе, также определяет временная поправка dTi, получение эталонных значений интервальных скоростей и проверка качества и подготовка методов МОП", СК, ВСП, АК, ГГК-П, значений геологических отметок и введение «нулевой точки» в годограф, посадка на верхний опорный горизонт литолого-стратиграфического комплекса и перенос точек годографа на ближайшие глинистые пачки, а также построение синтетической сейсмограммы в процессе получения годографа, после выполнения определенных процедур. Кроме того, соответствие, степени подобия экстремумов отражений от опорных и целевых горизонтов синтетической сейсмограммы экстремумам отражений от опорных и целевых горизонтов разреза МОП" достигается путем использования временного разреза МОП" с сохранением относительных амплитуд, приведением его к нульфазовому импульсу; применением помимо динамического критерия оценки колебательных процессов - функции взаимной корреляции (ФВК), критерия кинематического - значения интервальной скорости. Это достигается путем: сравнения значений интервальной скорости полученного годографа после привязки VMHT TIPHB С определенными ранее по данным СК и ВСП эталонными значениями интервальных скоростей, характерными для определенного интервала и геологического района, как дополнительный технический результат это определяет возможность контроля качества интервальной скорости годографа. Кроме того, предложенный «Способ получения априорного годографа для выполнения литолого- стратиграфической привязки», в совокупности заявленных признаков, обеспечивает надежный, достоверный результат, поскольку проводится проверка качества и подготовка временного разреза МОП" с сохранением относительных амплитуд и приведением его к нульфазовому импульсу, помимо ФВК используются другие критерии оценки, а именно, данные о распределении эталонных интервальных скоростей между основными реперными и целевыми горизонтами в данном районе и данном интервале, полученные ранее по данным проведенных сейсмокаротажа и вертикального сейсмического профилирования в точке скважины в координатах 2Д или ЗД. Следует заметить, что точность и достоверность результатов повышает также то, что используют данные о диапазонах значений интервальных скоростей, полученных лабораторным путем, т.е. измерением в соответствующем интервале на образцах горных пород в скважинах нужного геологического района; использованием осредненных значений АК в определенном интервале; построение синтетической сейсмограммы в процессе получения годографа, после выполнения определенных процедур. Точность результатов достигается также за счет равномерного внесения в каждую i- тую точку годографа, следующую за точкой верхнего опорного горизонта поправки dTi, что позволяет получить значения интервальной скорости по данным АК без «ступеней», соответственно, тренду их постепенного увеличения с глубиной. Кроме всего прочего, дополнительным техническим результатом заявляемого изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» является то, что окончательный годограф может быть использован для пересчета временного разреза в глубинный разрез, а синтетическая сейсмограмма, кривая акустического импеданса могут быть использованы для выполнения инверсии, либо при использовании нескольких таких годографов и кривых от различных скважин площади тех же самых результатов (глубинного разреза и инверсии) в координатах ЗД. Заявляемый «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» позволяет также обоснованно уточнять геологические отметки в районах сложного строения (например, в районах с аномальным строением баженовской свиты), т.к. позволяет осуществить контроль по эталонным значениям интервальных скоростей в точке скважины, полученных в координатах 2Д или ЗД и характерных для соответствующего интервала и определенного геологического района, причем контроль также можно осуществлять например, путем сравнения с диапазоном изменения значений интервальных скоростей, полученных лабораторным путем, посредством измерений в соответствующем интервале на образцах горных пород в скважинах данного геологического района, сравнением с осредненными значениями интервальных скоростей, полученных ранее по данным акустического каротажа на данном интервале. Проведенный заявителем анализ уровня техники, включающий поиск по патентным и научно-техническим источникам информации и выявление источников, содержащих сведения об аналогах заявленного изобретения «Способ получения априорного годографа для выполнения литолого- страти графической привязки», позволил установить, что заявитель не обнаружил источник, характеризующийся признаками, тождественными всем существенным признакам заявленного технического решения. По имеющимся у заявителя сведениям, совокупность существенных признаков заявляемого изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки», не известно из уровня техники, что позволяет сделать вывод о соответствии заявляемого изобретения «Способ получения априорного годографа для выполнения литолого- стратиграфической привязки» критерию "новизна". Определение из перечня выявленных аналогов прототипа, как наиболее близкого по совокупности признаков аналога, позволил выявить совокупность существенных, по отношению к усматриваемому заявителем техническому результату, отличительных признаков в заявляемом изобретении «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки», изложенных в формуле изобретения. Следовательно, заявляемое изобретение «Способ получения априорного годографа для выполнения литолого- стратиграфической привязки» соответствует критерию "новизна".
Для проверки соответствия заявляемого изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» критерию
"изобретательский уровень" заявитель провел дополнительный поиск известных решений, чтобы выявить совокупность признаков, совпадающих с отличительными от прототипа признаками заявляемого изобретения «Способ получения априорного годографа для выполнения литолого- стратиграфической привязки». Результаты поиска показали, что заявляемое изобретение «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» не вытекает для специалиста явным образом из известного уровня техники, поскольку из уровня техники, определенного заявителем, не выявлено влияние предусматриваемых существенными признаками заявленного изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» преобразований для достижения технического результата. Следовательно, заявленное изобретение «Способ получения априорного годографа для выполнения литолого- стратиграфической привязки» соответствует критерию "изобретательский уровень".
Таким образом, изложенные сведения свидетельствуют о выполнении при использовании в заявленном изобретении «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» совокупности условий в том виде, как заявляемое изобретение «Способ получения априорного годографа для выполнения литолого- стратиграфической привязки» охарактеризовано в формуле изобретения, т.е. подтверждена возможность ее осуществления с помощью описанного в заявке примера конкретного выполнения. Технологические приемы, воплощающие заявленное изобретение «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» при его осуществлении, способны обеспечить достижение усматриваемого заявителем технического результата, а именно снижение трудоемкости, повышение точности и надежности эксплуатационных и технических качеств, следовательно, заявленное изобретение «Способ получения априорного годографа для выполнения литолого- стратиграфической привязки» соответствует критерию "промышленная применимость".
Совокупность существенных признаков, характеризующих сущность изобретения «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» может быть многократно использована в процессе получения априорных годографов для выполнения литолого- стратиграфических привязок с получением технического результата, заключающегося в снижении трудоемкости, повышении точности, надежности, эксплуатационных и технических качеств, что позволяет экономически выгодно и быстро получать априорный годограф, проводить построение синтетической сейсмограммы и кривой акустической жесткости для выполнения литолого-стратиграфической привязки, увязки сейсмических и геологических границ, получать элементы априорной скоростной модели для выполнения инверсии сейсмических данных в количественные характеристики резервуара, т.е повышать достоверность интегрированного геологического анализа с целью создания геолого- технологических моделей, трехмерных цифровых геолого- гидродинамических моделей нефтяных и газовых месторождений.
Сущность заявляемого изобретения «Способ получения априорного годографа для выполнения литолого- стратиграфической привязки» и поясняется схемами отражающими приемы выполнения способа:
- на фиг. 1 изображена схема получения априорного годографа для выполнения литолого-стратиграфической привязки.
- фиг.2 - схема получения априорного годографа для выполнения экспресс-посадки.
Способ получения априорного годографа для выполнения литолого-стратиграфической привязки выполняли следующим образом.
ПРИМЕР 1.
1 . Первоначально проводили сейсмокаротаж (СК) в ближайшей скважине. Получали продольный вертикальный годограф проходящей волны Tc =f (Н), где Тек - время регистрации первых вступлений проходящей волны. Проверяли качество годографа СК, получив зависимость интервальной скорости, как функцию глубины VHHT_CK = f (Н). Т.к. качество годографа удовлетворяло, брали его как исходный.
2. Затем проводили акустический каротаж (АК) в скважине, в которой выполняли литолого-стратиграфическую привязку. Записывали dtAK = f (Н)- интервальное время пробега продольной волны. Проводили плотностной гамма-гамма каротаж (ГГК-П), записывали диаграмму изменений рассеянного гамма излучения как функцию глубины. Проверяли качество этих данных. Подготавливали диаграмму интервального времени акустического каротажа (АК) при этом:
- последовательно сшивали интервалы по глубине;
- редактировали значения интервального времени в интервалах образования каверн.
- оцифровывали диаграмму интервального времени пробега АК с шагом дискретизации 1 м.
Если АК в скважине не проводился моделировали из других методов ГИС, например, из диаграммы метода кажущегося сопротивления (КС). Далее пересчитывали интервальное время в интервальную скорость как Уинт_АК = f(H) и переходили от показаний ГГК-П к кривой плотности пород δπ=ί(Η).
3. После этого получали значения отметок глубин геологических границ, к которым необходимо найти и привязать соответствующие сейсмические горизонты. Проверяли геологические отбивки по диаграммам геофизических исследований скважин (ГИС) и данным сейсморазведки.
4. Получали окончательный суммарный временной разрез метода общей глубинной точки (МОП") с сохранением относительных амплитуд, приводили его к нульфазовому импульсу и проводили корреляцию опорных (реперных) отражающих горизонтов (ОГ) «С», «Г», «М1 », «М», «Б», «А».
5. После этого через значения глубин опорных и целевых сейсмических горизонтов, соответствующие времена регистрации первых вступлений проходящей волны от опорных и целевых горизонтов получали эталонные значения интервальных скоростей по данным СК и ВСП между этими горизонтами в координатах 2Д, рассчитав VMHT_CK_(B-A) =2dH/dT, где dH - мощность интервала (в метрах), между значениями геологических отметок, соответствующих опорным ОГ «Б» и «А» по диаграммам ГИС; dT- временная мощность интервала (в секундах) между опорными ОГ «Б» и «А», которые служили критериями точности привязки.
6. Далее в качестве исходного сигнала брали импульс U(Tj) с временного разреза МОП" из выбранного целевого интервала временного разреза в районе скважины и оценивали частотный диапазон и форму импульса.
7. Вслед за этим проводили коррекцию (калибровку) значений VAK = f (Н) по сейсмокаротажу. Эта процедура необходима для исключения постоянной ошибки в данных АК и приведения данных АК и СК, в одни координаты. Полученную кривую интервального времени сохраняли. После этого получали диаграмму акустической жесткости γπ =VMHT_AK *δπ, как результат произведения диаграммы интервальной скорости VMHT_AK = f(H) и диаграммы плотности пород δπ= f(H)
8. После этого осуществляли перевод исходной информации из масштаба глубин во временной масштаб V1 (Н) в V (Т).
9. Затем рассчитывали и формировали массив коэффициентов отражения при нормальном падении волны на
Figure imgf000033_0001
границу раздела упругих сред Р^Т) К^^+^ +^+у^убт^ъ) где индекс i соответствует параметрам слоя, из которого падает волна, индекс i+1 - параметрам слоя, находящегося под границей. 10. Вслед за этим рассчитывали синтетическую сейсмограммы как результат свертки массива коэффициентов отражения Ρι(Τ) с исходным сигналом U(Tj).
Pu(Tj) = u(Ti)Pi{Tj - 1), где
q=T70.002 - число отсчетов исходного сигнала
и осуществляли контроль качества годографа и синтетической сейсмограммы, оценивали акустические свойства разреза, уточняли значение геологической отметки опорных отражающих горизонтов (ΟΓ)·
1 1. Затем осуществляли посадку на верхний опорный горизонт литолого-стратиграфического комплекса, например ОГ «Г» введением постоянной поправки с!Т_верх_опор, поскольку чаще всего данные АК и СК были получены в разных скважинах, то по этой причине, а так же из-за разных подходов к определению интервальной скорости комплексируемых методов АК, СК, и МОГТ в полученной синтетической сейсмограмме не достигается соответствия значений экстремумов отражений от первого опорного горизонта синтетической сейсмограммы значениям экстремумов соответствующего опорного горизонта временного разреза МОГТ, эту разницу dT измеряли на разрезе МОГТ и вводили в исходный годограф.
12. Вслед за этим вновь рассчитывали синтетическую сейсмограмму с исправленным годографом СК. Проверяли соответствие значений экстремумов отражений от верхнего опорного горизонта литолого-стратиграфического комплекса синтетической сейсмограммы и соответствующего опорного горизонта временного разреза МОГТ, наблюдали соответствие значений экстремумов. Далее осуществляли контроль качества годографа и синтетической сейсмограммы с помощью ФВК в ближнем интервале и сравнивали с эталонными значениями интервальных скоростей.
13. После этого в годограф СК вводили "нулевую точку" (глубина Н=0 и время 2Т =0), а также две точки ниже забоя скважины
14. Снова рассчитывали синтетическую сейсмограмму с исправленным годографом СК. Осуществляли контроль качества годографа и синтетической сейсмограммы с помощью ФВК в ближнем интервале и сравнением с эталонными значениями интервальных скоростей.
15. Следом осуществляли посадку на нижний опорный горизонт литолого-стратиграфического комплекса, например ОГ «Б». Такое же несоответствие значений экстремумов о котором упоминалось выше в п.1 1 , а именно: несоответствие значений экстремумов отражений от нижнего опорного горизонта литолого-стратиграфического комплекса синтетической сейсмограммы значениям экстремумов соответствующего опорного горизонта временного разреза МОП" наблюдали и в интервале ОГ «Б». Измеряли разницу dTB = tB огт - tB ск на временном разрезе МОП".
Рассчитывали поправку dTi = Н|*^, где dTi - поправка, вводимая в каждую i-тую точку последнего годографа, следующую за точкой верхнего опорного горизонта; dTB - = ts огт - tB ск - поправка, которую нужно ввести в последний годограф для посадки на опорный ОГ «Б», т.е. совмещения во временном масштабе значения геологической отметки ОГ «Б» и значения времени экстремума соответствующей отраженной волны;
tB orr - значение времени снятое с временного разреза МОП" (в точке скважины) для экстремума отраженной волны, соответствующей опорному ОГ «Б»;
ecK- значение времени, снятое с последнего годографа СК для геологической отметки опорного ОГ «Б»;
Hi - значение геологической отметки i- ой точки последнего годографа СК;
НБ - значение геологической отметки опорного ОГ «Б» по диаграммам ГИС;
Поправку Тиспр=Тпосл+сЛП вводили в последний полученный годограф п.13.
16. Далее вновь рассчитывали синтетическую сейсмограмму с использованием исправленного годографа. Наблюдали соответствие значений экстремумов в интервале ОГ «Б». Осуществляли контроль качества годографа и синтетической сейсмограммы с помощью ФВК в соответствующем интервале и сравнением с эталонными значениями интервальных скоростей.
17. Вслед за этим переносили точки годографа на акустически слабые границы - глины.
18. После этого рассчитывали окончательную синтетическую сейсмограмму с использованием исправленного годографа и осуществляли контроль качества годографа и синтетической сейсмограммы с помощью ФВК в соответствующем интервале и сравнением с эталонными значениями интервальных скоростей. Затем сохраняли априорный годограф, синтетическую сейсмограмму, кривую акустической жесткости.
19. Следует заметить, что для скважин с глубинами забоя в самом нижнем литолого-стратиграфическом комплексе, например, ниже ОГ «Б», сначала фиксировали точку годографа, соответствующую значению экстремума ОГ «Б». (НБ - значение геологической отметки, ТБ - значение времени экстремума отраженной волны ОГ «Б», снятое с временного разреза МОГТ).
20. Затем вводили точки годографа, соответствующие глинистым пачкам.
21. После этого рассчитывали значение времени для точки забоя скважины годографа как
Тзаб = ТБ + dT, где dT=2dH/ViiHT (Б-А), где
dH - расстояние в м от значения геологической отметки опорного ОГ «Б» до значения глубины точки забоя;
ТБ - значение времени, снятое с временного разреза МОГТ и соответствующее времени экстремума отраженной волны от опорного ОГ «Б»;
VMHT (Б-А) - эталонная интервальная скорость по данным СК и ВСП, полученная в п.5 для самого нижнего литолого- стратиграфического комплекса, например интервала между ОГ «Б» и ОГ «А»
Так же осуществляли расчет значений времени для двух последних введенных точек годографа с использованием эталонной интервальной скорости данного интервала.
22. Вслед за этим рассчитывали синтетическую сейсмограмму по окончательному годографу СК и осуществляли контроль качества годографа и синтетической сейсмограммы с помощью ФВК в соответствующем интервале и сравнением с эталонными значениями интервальных скоростей
23. После этого сохраняли априорный годограф, синтетическую сейсмограмму и кривую акустической жесткости.
Предложенный «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» снижает трудоемкость и себестоимость, увеличивает скорость получения результата, а также повышает точность и достоверность соответствия значений времени экстремумов отражений от опорных и целевых горизонтов синтетической сейсмограммы значениям времени экстремумов отражений от опорных и целевых горизонтов временного разреза МОП" и эффективность как отношение затрат к результату.
ПРИМЕР 2.
1. Первоначально проводили сейсмокаротаж (СК) в ближайшей скважине и получали исходный продольный вертикальный годограф проходящей волны TcK=f (Н), где
Тек - время регистрации первых вступлений проходящей волны. Проверяли качество годографа СК, получив зависимость интервальной скорости, как функцию глубины Уинт СК = f (Н). Если качество годографа удовлетворяло, брали его как исходный.
2. Затем получали значения отметок глубин геологических границ, которым необходимо найти соответствующие сейсмические горизонты и проверяли геологические отбивки по диаграммам геофизических исследований скважин (ГИС) и данным сейсморазведки. 3. После этого получали окончательный суммарный временной разрез метода общей глубинной точки (МОП") с сохранением относительных амплитуд, приводили его к нульфазовому импульсу и проводили корреляцию опорных (реперных) отражающих горизонтов (ОГ) «С», «Г», «М1 », «М», «Б», «А».
4. В последующем через значения глубин опорных и целевых сейсмических горизонтов, соответствующие времена регистрации первых вступлений проходящей волны от опорных и целевых горизонтов получали эталонные значения интервальных скоростей по данным СК и ВСП между этими горизонтами в координатах ЗД, рассчитав VMHT_CK_(5-A) =2dH/dT, где dH - мощность интервала (в метрах), между значениями геологических отметок, соответствующих опорным ОГ «Б» и «А» по диаграммам ГИС; dT- временная мощность интервала (в секундах) между опорными ОГ «Б» и «А», которые служили критериями точности посадки.
5. Далее осуществляли перевод исходной информации из масштаба глубин во временной масштаб V1 (Н) в V (Т).
6. Вслед за этим осуществляли посадку на верхний опорный горизонт литолого-стратиграфического комплекса, например ОГ «Г» введением постоянной поправки dT_eepx_onop. Из-за разных подходов к определению интервальной скорости комплексируемых методов СК, ВСП и МОП" при пересчете во временной масштаб не достигается соответствия значений времени геологической отметки верхнего опорного горизонта значениям времени экстремума соответствующего опорного горизонта временного разреза МОП". Эту разницу dT измеряли на разрезе МОГТ и вводили в исходный годограф.
7. После чего проверяли соответствие значений времени геологической отметки верхнего опорного горизонта полученного годографа значениям времени экстремума соответствующего опорного горизонта временного разреза МОГТ и наблюдали соответствие значений времен экстремума и геологической отметки, а также осуществляли контроль качества годографа сравнением с эталонными значениями интервальных скоростей.
8. Затем в годограф СК вводили "нулевую точку" (глубина
Н=0 и время 2Т=0), а также 2 точки ниже забоя скважины
9. Далее осуществляли контроль качества годографа сравнением с эталонными значениями интервальных скоростей.
10. После чего осуществляли посадку на нижний опорный горизонт литолого-стратиграфического комплекса, например ОГ
«Б». Такое же несоответствие времен геологической отметки во временном масштабе и значения времени экстремума нижнего опорного горизонта на временном разрезе МОГТ, о котором шла речь в п.6, а именно: несоответствие значений времени геологической отметки нижнего опорного горизонта полученного годографа значениям времени экстремума соответствующего опорного горизонта временного разреза МОГТ наблюдалось и в интервале ОГ «Б». Измеряли разницу с1ТБ = tB огт - tB ск на временном разрезе МОГТ и рассчитывали поправку dTi = Hi*^ , где
dTi - поправка, вводимая в каждую i-тую точки последнего годографа, следующие за точкой верхнего опорного горизонта; ОГТ ск
dTB - = tB - tB - поправка, которую нужно ввести в последний годограф для посадки на горизонт «Б», т.е. совмещения в масштабе времен значения геологической отметки ОГ «Б» и значения времени экстремума соответствующей отраженной волны;
tB orr - значение времени снятое с временного разреза МОП" (в точке скважины) для экстремума отраженной волны, соответствующей горизонту «Б»;
tB CK- значение времени, снятое с последнего годографа СК для геологической отметки горизонта «Б»;
Hi - значение геологической отметки i- ой точки последнего годографа СК;
НБ - значение геологической отметки отражающего горизонта «Б» по диаграммам ГИС;
Одновременно вводили поправку Тиспр=Тпосл+с1"П в последний годограф п.8.
1 1. Затем наблюдали соответствие значения времени геологической отметки во временном масштабе и значения времени экстремума отраженной волны в интервале ОГ «Б» и вновь осуществляли контроль качества годографа сравнением с эталонными значениями интервальных скоростей.
12. Далее переносили точки годографа на акустически слабые границы - глины.
13. Вновь осуществляли контроль качества годографа сравнением с эталонными значениями интервальных скоростей и сохраняли априорный годограф.
14. Причем для скважин с глубинами забоя в самом нижнем литолого-стратиграфическом комплексе, например, ниже ОГ «Б» фиксировали точку годографа, соответствующую ОГ «Б» (НБ - значение геологической отметки ОГ «Б», ТБ - значение времени экстремума отраженной волны ОГ «Б», снятое с временного разреза МОП").
15. После этого вводили точки годографа, соответствующие глинистым пачкам.
16. Затем рассчитывали значение времени для точки забоя скважины годографа как Тзаб = ТБ + dT, где dT=2dH//ViiHT (Б-А), где,
dH - расстояние в м от значения геологической отметки ОГ «Б» до значения глубины точки забоя;
ТБ - значение времени, снятое с временного разреза МОП" и соответствующее экстремуму от ОГ «Б»;
VMHT (Б-А) - интервальная скорость по данным СК и ВСП, полученная в п.4. для самого нижнего литолого- стратиграфического комплекса, например для интервала между ОГ «Б» и ОГ «А» и так же осуществляли расчет значений времени для двух введенных последних точек годографа с использованием эталонной интервальной скорости данного интервала.
17. Вслед за этим осуществляли контроль качества годографа сравнением с эталонными значениями интервальных скоростей
18. Далее сохраняли априорный годограф.
Предложенный «Способ получения априорного годографа для выполнения литолого-стратиграфической привязки» снижает трудоемкость и себестоимость, увеличивает скорость получения результата, а также повышает точность и достоверность соответствия значений времени геологических отметок опорных и целевых ОГ во временном масштабе значениям времени экстремумов соответствующих опорных и целевых горизонтов временного разреза МОП" и эффективность как отношение затрат к результату.

Claims

Ф О Р М У Л А И З О Б Р Е Т Е Н И Я
1. Способ получения априорного годографа для выполнения литолого-стратиграфической привязки, путем получения и подготовки данных акустического каротажа, 5 сейсмокаротажа, плотностного гамма-гамма каротажа, калибровки данных акустического каротажа, перевода во временной масштаб, расчета последовательности коэффициентов отражения, свертки их с исходным сигналом, т.е построения синтетической сейсмограммы и кривой акустической ю жесткости с повторением циклов, о т л и ч а ю щ и й с я тем,
что он включает в себя ряд последовательных действий, при которых независимо получают и подготавливают данные методов общей глубинной точки, сейсмокаротажа, вертикального сейсмического профилирования, акустического
15 каротажа, плотностного гамма-гамма каротажа, значения геологических отметок и проверяют качество этих данных, а также получают эталонные значения интервальных скоростей, после чего получают исходный годограф и рассчитывают синтетическую сейсмограмму, затем проводят контроль
20 качества, далее вводят постоянную временную поправку для посадки на верхний опорный горизонт литолого- стратиграфического комплекса, затем вновь рассчитывают синтетическую сейсмограмму, вновь проводят контроль качества, после чего вводят «нулевую точку» в полученный
25 годограф и дополнительно одну, и более точек ниже забоя скважины, проводят расчет синтетической сейсмограммы и контроль качества, вслед за этим рассчитывают и вводят поправку для посадки на нижний опорный горизонт литолого- стратиграфического комплекса, по формуле
Figure imgf000045_0001
dTi - поправка, вводимая в каждую i-тую точку последнего годографа, следующую за точкой верхнего опорного горизонта;
dTH - = tH огт - tH ск - поправка, которую нужно ввести в последний годограф для посадки на нижний опорный горизонт, т.е. совмещения во временном масштабе значения геологической отметки нижнего опорного горизонта и значения времени экстремума соответствующей отраженной волны;
tH огт - значение времени снятое с временного разреза метода общей глубинной точки (в точке скважины) для экстремума отраженной волны соответствующей нижнему опорному горизонту;
tH - значение времени, снятое с последнего годографа сейсмокаротажа, вертикального сейсмического профилирования для геологической отметки нижнего опорного горизонта;
Hi - значение геологической отметки i- ой точки последнего годографа сейсмокаротажа, вертикального сейсмического профилирования ;
Нн -значение геологической отметки нижнего опорного горизонта по диаграммам геофизических исследований скважины;
после этого вновь рассчитывают синтетическую сейсмограмму с последующим контролем качества, и следом переносят точки полученного годографа на ближайшие акустически слабые границы, далее повторно рассчитывают синтетическую сейсмограмму с последующим контролем качества и получают априорный годограф, сохраняют его, синтетическую сейсмограмму и кривую акустической жесткости.
2. Способ по п.1, о т л и ч а ю щ и й с я тем, что в качестве контроля качества привязки используют функцию взаимной корреляции, которая является динамическим критерием оценки колебательных процессов, и являющиеся кинематическим критерием эталонные значения интервальных скоростей данного интервала, полученные ранее по данным сейсмокаротажа и вертикального сейсмического профилирования в координатах 2Д, ЗД с использованием значений глубин опорных и целевых горизонтов, соответствующих времен регистрации первых вступлений проходящей волны.
3. Способ по п.1, о тл и ч а ю щ и й ся тем, что в качестве контроля качества привязки используют функцию взаимной корреляции , которая является динамическим критерием оценки колебательных процессов, и являющиеся кинематическим критерием эталонные значения интервальных скоростей, которые получают лабораторным путем посредством измерения в соответствующем интервале на образцах горных пород в скважинах данного геологического района.
4. Способ по п.1, о т л и ч а ю щ и й с я тем, что в качестве контроля качества привязки используют функцию взаимной корреляции, которая является динамическим критерием оценки колебательных процессов, и являющиеся кинематическим критерием эталонные значения интервальных скоростей, которые получают осреднением значений интервальных скоростей, полученных ранее по данным акустического каротажа на данном интервале.
5. Способ по п.1, о т л и ч а ю щ и й с я тем, что в качестве исходного годографа берут осредненныи годограф, как результат аппроксимации годографов сейсмокаротажа и вертикального сейсмического профилирования данной площади.
6. Способ по п.1, о т л и ч а ю щ и й с я тем, что в качестве исходного годографа берут годограф, значения которого сняты со значений глубин^ опорных и целевых горизонтов, соответствующих времен регистрации первых вступлений проходящей волны, используемых для получения эталонных значений интервальной скорости в координатах 2Д, зд.
7. Способ по п.1, о т л и ч а ю щ и й с я тем, что в качестве исходного годографа берут годограф, значения которого сняты с тренда изменения диаграммы интервального времени акустического каротажа как функции глубины.
8. Способ по п.1, о т л и ч а ю щ и й ся тем, что в качестве исходного сигнала берут импульс выделенный в целевом интервале вблизи скважины с суммарного временного разреза метода общей глубинной точки полученного с сохранением относительных амплитуд и приведенного к нульфазовому импульсу.
9 Способ по п .1, о т л и ч а ю щ и й с я тем, что в качестве исходного сигнала берут нульфазовый модельный импульс подобранный по частоте к частотам суммарного временного разреза метода общей глубинной точки, полученного с сохранением относительных амплитуд и приведенного к нульфазовому импульсу.
10. Способ по п.1, о тл и ч а ю щ и й ся тем, что в скважинах, с глубинами забоя в самом нижнем литолого- стратиграфическом комплексе, фиксируют точку годографа соответствующую значению экстремума опорного отражающего горизонта, являющегося верхней границей этого литолого- стратиграфического комплекса, после чего вводят точки соответствующих глинистых горизонтов данного комплекса и рассчитывают значение времени для точки глубины забоя скважины и добавленных ниже точек годографас использованием значений эталонных интервальных скоростей данного интервала после чего повторяют расчет синтетической сейсмограммы и следом осуществляют контроль качества и получают априорный годограф, сохраняют его, синтетическую сейсмограмму и кривую акустической жесткости.
11. Способ по п.1 , о т л и ч а ю щ и й с я тем, в качестве функции для перевода во временной масштаб используют годографы сейсмокаротажа, вертикального сейсмического профилирования без участия данных акустического каротажа и построения синтетической сейсмограммы.
PCT/RU2014/000482 2013-10-10 2014-07-02 Способ получения априорного годографа для выполнения литолого стратиграфической привязки WO2015053659A1 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
RU2013145637 2013-10-10
RU2013145637/28A RU2541091C1 (ru) 2013-10-10 2013-10-10 Способ получения априорного годографа для выполнения литолого-стратиграфической привязки

Publications (1)

Publication Number Publication Date
WO2015053659A1 true WO2015053659A1 (ru) 2015-04-16

Family

ID=52813389

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/RU2014/000482 WO2015053659A1 (ru) 2013-10-10 2014-07-02 Способ получения априорного годографа для выполнения литолого стратиграфической привязки

Country Status (2)

Country Link
RU (1) RU2541091C1 (ru)
WO (1) WO2015053659A1 (ru)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109975876A (zh) * 2019-03-20 2019-07-05 中国石油化工股份有限公司 一种基于构造层位的井震融合速度模型的建模方法
CN111458751A (zh) * 2019-01-22 2020-07-28 中国石油天然气股份有限公司 一种综合时深标定方法及装置
CN114034266A (zh) * 2021-09-30 2022-02-11 湖南天功测控科技有限公司 一种优化阻尼换能器的声波法基桩嵌岩深度检测方法
CN116201535A (zh) * 2023-02-06 2023-06-02 北京月新时代科技股份有限公司 油气藏目标井标志地层自动划分方法、装置及设备

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115407399A (zh) * 2021-05-27 2022-11-29 中国石油天然气股份有限公司 适于地震数据极性判断方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0716752B1 (en) * 1993-09-03 1998-11-11 Spiral Services Limited Evaluation method and apparatus in seismics
US20100057367A1 (en) * 2008-09-02 2010-03-04 Paal Kristiansen Processing seismic data in common group-center gathers
RU2490677C2 (ru) * 2011-11-28 2013-08-20 Александр Алексеевич Архипов Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2449322C1 (ru) * 2010-11-02 2012-04-27 Юрий Андреевич Романенко Способ построения сейсмического глубинного разреза
RU2463628C1 (ru) * 2011-04-08 2012-10-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Новосибирский национальный исследовательский государственный университет" (Новосибирский государственный университет, НГУ) Способ построения сейсмического глубинного разреза

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0716752B1 (en) * 1993-09-03 1998-11-11 Spiral Services Limited Evaluation method and apparatus in seismics
US20100057367A1 (en) * 2008-09-02 2010-03-04 Paal Kristiansen Processing seismic data in common group-center gathers
RU2490677C2 (ru) * 2011-11-28 2013-08-20 Александр Алексеевич Архипов Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
GOGENENKOV G. N.: "Raschet i primenenie sinteticheskikh seismogramm.", MOSCOW, IZDATELSTVO ''NEDRA, 1972, pages 46 - 49 , 77-85 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111458751A (zh) * 2019-01-22 2020-07-28 中国石油天然气股份有限公司 一种综合时深标定方法及装置
CN109975876A (zh) * 2019-03-20 2019-07-05 中国石油化工股份有限公司 一种基于构造层位的井震融合速度模型的建模方法
CN114034266A (zh) * 2021-09-30 2022-02-11 湖南天功测控科技有限公司 一种优化阻尼换能器的声波法基桩嵌岩深度检测方法
CN114034266B (zh) * 2021-09-30 2024-04-26 湖南天功测控科技有限公司 一种优化阻尼换能器的声波法基桩嵌岩深度检测方法
CN116201535A (zh) * 2023-02-06 2023-06-02 北京月新时代科技股份有限公司 油气藏目标井标志地层自动划分方法、装置及设备
CN116201535B (zh) * 2023-02-06 2024-02-09 北京月新时代科技股份有限公司 油气藏目标井标志地层自动划分方法、装置及设备

Also Published As

Publication number Publication date
RU2541091C1 (ru) 2015-02-10

Similar Documents

Publication Publication Date Title
US10067253B2 (en) Method for determining sedimentary facies using 3D seismic data
EP2920615B1 (en) Process for creating image gathers
US20130223187A1 (en) Geological Structure Contour Modeling and Imaging
US20090213692A1 (en) Method for three dimensional seismic travel time tomography in transversely isotropic media
Maxwell et al. Anisotropic velocity modeling for microseismic processing: Part 1—Impact of velocity model uncertainty
US7460437B2 (en) Seismic data processing method and system for migration of seismic signals incorporating azimuthal variations in the velocity
WO2015053659A1 (ru) Способ получения априорного годографа для выполнения литолого стратиграфической привязки
US10739481B2 (en) 2D multiline seismic reflection tomography with seismic-tie constraint
Ivanov et al. The inverse problem of refraction travel times, part II: Quantifying Refraction Nonuniqueness using a three-layer model
CN115877449A (zh) 用于在勘测区域内获得地下堆叠图像的计算机实现方法
Yan et al. Building subsurface models with horizon-guided interpolation and deep learning: Application to the Volve field data
Grechka et al. Velocity analysis for tilted transversely isotropic media: A physical modeling example
Gei et al. Feasibility of estimating vertical transverse isotropy from microseismic data recorded by surface monitoring arrays
RU2449322C1 (ru) Способ построения сейсмического глубинного разреза
Khoshnavaz et al. Oriented surface passive seismic location using local slopes
US5587942A (en) 3D wave equation migration of a 2D grid of seismic data
Velásquez et al. Depth-conversion techniques and challenges in complex sub-Andean provinces
CN112305595B (zh) 基于折射波分析地质体结构的方法及存储介质
RU2659753C1 (ru) Способ геофизической разведки
Xu et al. A new fracture detection method based on full-azimuth anisotropic migration and shaping regularization
RU2221262C1 (ru) Способ сейсмической разведки для изучения осадочного чехла при наличии сильно изрезанных акустически жестких границ (варианты)
US10401515B2 (en) Estimation of water properties from seismic data
Al-Heety An Evaluation between Time-term, Reciprocal Time and Refraction Tomography Analysis Methods for obtaining 2-D shallow Seismic Velocity Models over Synthetic Traveltimes
Yan et al. Building subsurface models with horizon-guided interpolation and deep learning: Applied to the Volve field
Valenta et al. 3D seismic travel time surveying a comparison of the time–term method and tomography (an example from an archaeological site)

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14852858

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 14852858

Country of ref document: EP

Kind code of ref document: A1