SU1427312A1 - Method of seismic probing for predicting geological section - Google Patents

Method of seismic probing for predicting geological section Download PDF

Info

Publication number
SU1427312A1
SU1427312A1 SU864123270A SU4123270A SU1427312A1 SU 1427312 A1 SU1427312 A1 SU 1427312A1 SU 864123270 A SU864123270 A SU 864123270A SU 4123270 A SU4123270 A SU 4123270A SU 1427312 A1 SU1427312 A1 SU 1427312A1
Authority
SU
USSR - Soviet Union
Prior art keywords
section
transverse
determined
longitudinal
geological
Prior art date
Application number
SU864123270A
Other languages
Russian (ru)
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 Южное Отделение Всесоюзного Научно-Исследовательского Института Геофизических Методов Разведки
Priority to SU864123270A priority Critical patent/SU1427312A1/en
Application granted granted Critical
Publication of SU1427312A1 publication Critical patent/SU1427312A1/en

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Description

(21)4123270/24-25(21) 4123270 / 24-25

(22)26.06.86(22) 06.26.86

(46) 30.09.88. Бюл. Р 36 (71) Южное отделение Всесоюзного научно-исследовательского института геофизических методов разведки (72) М.М.Раджабов(46) 09/30/88. Bul R 36 (71) Southern Branch of the All-Union Scientific Research Institute of Geophysical Methods of Exploration (72) MM Radjabov

(53)550.834(088.8)(53) 550.834 (088.8)

(56)Гурвич И.Н. Сейсмическа  разведка . М. : Недра, 1980, с. 87-294.(56) Gurvich I.N. Seismic exploration. M.: Nedra, 1980, p. 87-294.

.Мембей В.И. Петодика многократных перекрытий в сейсморазведке. М.: Недра , 1985, с. 99, рис. 46а.,.Membey V.I. Petodika multiple overlaps in seismic exploration. M .: Nedra, 1985, p. 99, fig. 46a.,

(54)СПОСОБ СЕПСМОЗОНДИРО ВАНИЯ ДЛЯ ПРОГНОЗИРОВАШИ ГЕОЛОГИЧЕСКОГО РАЗРЕЗА(54) METHOD OF SEPSMOZONDRO VANIYA FOR FORECASTING GEOLOGICAL CUT

(57)Изобретение относитс  к сейсмической разведке и может быть использовано в геологоразведочном деле при проектировании бурени  поисково- разведочных скважин. Цель - повышение точности и детальности прогнозировани  геологического разреза и его объемного расчленение по вертикали, до любой глубины в заданной на дневной поверхности области. Способ осно .ван на разбивке продольного профил (57) The invention relates to seismic exploration and can be used in exploration in the design of drilling exploratory wells. The goal is to improve the accuracy and detail of the prediction of the geological section and its volumetric dismemberment along the vertical, to any depth in the area specified on the day surface. The method is based on the longitudinal profile.

.с центром в заданной области, возбуждении колебаний в пунктах, расположенных вдоль этого профил  и симметрично его центру, приеме колебаний отраженных волн в пунктах, расположенных на равных отрезках совмещенных продольных и поперечных профилей, центры которых  вл ютс  -симметричными с пунктами возбуждени  относительно центра продольного профил , определении пространственного сноса и выполнении селекции вертикального пол  лучевых скоростей по двум- сечени м , что позвол ет повысить его помехоустойчивость и более точно найти распределение усредненных значений лyчeJыx скоростей и соответствующих им вертикальных времен, по которым определ ют пластовые (интервальные ) скорости,после чего по идентичным глубинам наход т отношение пластовых скоростей, определенных порознь по системам поперечных и продольных профилей, по величине которого и изменени м с глубиной прогнозируют геологический разрез и суд т о степени его неоднородности. 1 ил.With a center in a given region, the excitation of oscillations in points located along this profile and symmetrically to its center, the reception of oscillations of reflected waves in points located on equal segments of combined longitudinal and transverse profiles, the centers of which are symmetrical with the points of excitation relative to the center of the longitudinal profile, determining the spatial drift, and performing the selection of the vertical radial velocity field over two sections, which makes it possible to increase its noise immunity and more accurately find the Determination of average values of peak speeds and the corresponding vertical times, which determine reservoir (interval) speeds, after which, at identical depths, the ratio of reservoir velocities determined separately by systems of transverse and longitudinal profiles is found, according to their magnitude and changes with depth, they predict geological section and judge the degree of its heterogeneity. 1 il.

ii

(L

j ю j you

соwith

Изобретение относитс  к сейсмической разведке и может быть использовано в геологоразведочном деле дл  прогнози ровани  геологического разреза, в г том числе разреза отвода и околосква жинного пространства проектируемых,, бур щихс  и пробуренных скважин не только на интервале глубин, вскрытых скважиной, но и ниже забо . tOThe invention relates to seismic exploration and can be used in geological prospecting to predict a geological section, including the withdrawal section and the near-well surface of the designed, drilled and drilled wells not only at the depth interval discovered by the well, but also below the bottom hole. tO

Целью изобретени   вл етс  повышение детальности и точности изучени  геологического разреза в пределах площадной системы наблюдений, ограничивающей исследуемую область. 15The aim of the invention is to increase the detail and accuracy of studying a geological section within an area observation system, limiting the area under investigation. 15

На чертеже показана схема наблюдений .The drawing shows the scheme of observations.

На схеме обозначены лини  1 наблюдени , секуща  центр исследуемой области , база 2 исследуемой области, центр 3 базы, совпадающий с центром линии наблюдени , систеь а 4 зондировани , совмещенна  система 5 приема , источники 6 возбуждени .The diagram shows the line 1 of observation, the secant of the center of the area under study, the base 2 of the area under study, the center 3 of the base, which coincides with the center of the line of observation, the system and 4 probing, the combined system 5 of reception, the sources 6 of excitation.

Способ осуществл етс  следующим образом.The method is carried out as follows.

Разбивают линию 1 наблюдени , центр которой совпадает с центром исследуемой области,определ ют базу 2 исследуемой области, равную половин длин отрезков продольных и поперечных расстановок. Длину отрезков, образующих совмещенную систему 5, устанавливают из услови  необходимой точности определени  зффективной скорости, величина которой используетс  дл  определени  лучевой скорости . На точность вычислени  эффективной скорости вли ет р д факто ров, главным из которых  вл етс  длина участка годографа отраженной волны. Со отношение, св зывающее длину участка годографа с искомой величиной эффективной скорости, можно представить в видеThe observation line 1 is split, the center of which coincides with the center of the study area, and the base 2 of the study area is determined to be equal to half the lengths of the segments of the longitudinal and transverse arrangements. The length of the segments forming the combined system 5 is established from the condition of the necessary accuracy of determining the effective speed, the value of which is used to determine the radial velocity. The accuracy of calculating the effective velocity is influenced by a series of factors, the main of which is the length of the hodograph section of the reflected wave. The ratio relating the length of the section of the hodograph to the desired effective velocity can be represented as

т - V Т |cJV/V| t - V T | cJV / V |

5050

- длина участка годографа, на котором производитс  измерение эффективной скорости- the length of the section of the hodograph, where the measurement of the effective velocity

V - V -

- предельна  абсолютна  ошибка в измерении времен отраженных волн; - предельна  относительна - the absolute error in measuring the time of the reflected waves is absolute; - extreme relative

ощибка в определении скорости.Error in determining the speed.

Эффективные скорости отраженных волн возрастают с увеличением рас сто ни  от источника возбуждени , стрем сь к наибольшей пластовой скорости пц. маис изучаемого интервала разреза. Поэтому, прин в V,m The effective velocities of the reflected waves increase with increasing distance from the excitation source, tending to the highest reservoir velocity pc. maize study interval cut. Therefore, pr in V, m

формулу (1) представимformula (1) we will present

ПА. WdXC PA Wdxc

в вцде:in all:

IdtlIdtl

Tdv7vTTdv7vT

пл. макс square Max

(2)(2)

Формулой (2) определ етс  минимально допустима  длина отрезка годографа отраженной волны дл  определени  эффективной скорости.Formula (2) determines the minimum allowable length of the hodograph segment of the reflected wave to determine the effective velocity.

Дл  поперечной расстановки L - длина отрезка, отсчитываема  в одну сторону от проекции, на нее источника .возбуждени , т.е. полна  длина поперечного годографа состоит из двух отрезков. Дл  обеспечени  необходимой точности определени  эффективных скоростей достаточно прин ть минимально допустимую длину каждого отрезка расстановок, вход щих в систему приема, равнуюFor a transverse arrangement, L is the length of a segment, measured in one direction from the projection, onto which the source of excitation, i.e. The full length of the transverse locus consists of two segments. To ensure the necessary accuracy in determining the effective velocities, it is sufficient to take the minimum allowable length of each segment of the arrangements included in the reception system, equal to

2L 2L

2V2V

пл максmax max

К,TO,

(3)(3)

где К where k

IdtlIdtl

IdV/Vl IdV / Vl

Поскольку поперечна  расстановка состоит из двух отрезков, то полна  ее длина равнаSince the transverse arrangement consists of two segments, its full length is equal to

4V4V

ПА.МОКСPA.MOKS

к.to.

(4)(four)

Следовательно, в совмещенной системе минимально допустима  длина поперечной расстановки принимаетс  в два раза больше минимально допустимой длины продольной расстановки, определ емой из соотнсЛиени Consequently, in the combined system, the minimum permissible length of the transverse arrangement is twice the length of the minimum permissible length of the longitudinal distribution, determined from the relation

2 V2 v

ПЛ МО КС,PL MO KS,

к.to.

(5)(five)

5050

На практике наибольшим ошибкам подвержено врем  прихода, волн. При мем, что Idtl измен етс  в пределах от 0,02 до 0,002 с и средн   величина относительной ошибки IdV/v в определении скорости определ етс  с точностью 5%. При этих значени х 0,04 К «i:. 0,4. Минимально допустима  длина каждого отрезка системы приема в этом случае может быть определена из услови :In practice, the arrival time of the waves is subject to the greatest errors. It is understood that Idtl varies from 0.02 to 0.002 s and the average value of the relative error IdV / v in the definition of velocity is determined with an accuracy of 5%. With these values, 0.04 K? I :. 0.4. The minimum length of each segment of the reception system in this case can be determined from the condition:

0.08 У„„ L 0,8 У,,.0.08 U „„ L 0,8 U ,,.

При VWhen V

Пл, ммцPl, mmc

0,1 км.0.1 km

1,5 1.5

4 км/с.4 km / s.

x 0,15x 0.15

6 км/с 0,5 L ё. 6 km / s 0.5 L ё.

При „,. 4,8 км.With „,. 4.8 km

При работе на каждой системе 5 приема, расположенной по одну сторону от центра исследуемой области, принимают отраженные волны, возбуж- денные из источника 6 возбуждени , расположенного по другую сторону от центра исследуемой области. В дальнейшем последовательно смеща  систему 5 приема и источник 6 возбуждени  в противоположные направлени  от центра исследуемой области и образу  при этом систему 4 зондировани , сохран ют указанный пор док возбуждени  и приема отраженных волн.When operating on each reception system 5 located on one side of the center of the area under study, receive reflected waves excited from the excitation source 6 located on the other side of the center of the area under study. Subsequently, sequentially displacing the receiving system 5 and the excitation source 6 in opposite directions from the center of the study area and forming the sounding system 4, the indicated order of excitation and reception of reflected waves is maintained.

Интервал между центрами системы 5 приема, расположенными по одну сторону от центра исследуемой обласми 6 возбуждени , расположенными в противоположном направлении от него, определ етс  из услови :The interval between the centers of the reception system 5, located on one side of the center of the investigated excitation regions 6, located in the opposite direction from it, is determined from the condition:

2uZ, 2uZ,

где Л X - интервал между центрамиwhere L X - the interval between the centers

системы приема и источниками возбуждени ; uZ - мощность пласта с наимень-шей пластовой скоростью 35 разреза.reception systems and excitation sources; uZ is the thickness of the reservoir with the smallest reservoir velocity 35 cut.

Использование пласта с наименьшей пластовой скоростью разреза обуслов-. лено необходимостью вы влени  пластаThe use of the reservoir with the lowest reservoir cutting speed caused-. it is necessary to reveal the formation

и с большей пластовой скоростью. Этим достигаетс  заданна  степень детализации изучаемого разреза.and with greater reservoir velocity. This achieves a given degree of detail of the section being studied.

Между величинами Л Z и V, „ц существует следующее соотношение:Between the values of L Z and V, „c there is the following relationship:

UZ 1,4 V,UZ 1,4 V,

ПА. М«НPA M "N

К.TO.

(8) формулу (7) в виде(8) formula (7) in the form

можно предU X 2,8 V,can be preu X 2.8 V,

К.TO.

При наименьшей пластовой скорости и заданной погрешности в определении скорости, равной 5%, Jol tj К 0,04. ТогдаWith the lowest reservoir velocity and a given error in determining the speed, equal to 5%, Jol tj K 0.04. Then

,,8-V,, 8-V

п(|. минn (|. min

0,1-У„д, мин . (10)0,1-U „d, min. (ten)

Дл  отраженных волн разной прот женности , регистрируемых на каждом из профилей системы, наход т лучевые скорости дл  соответствующих глубин по формуле:For reflected waves of different lengths recorded on each of the system profiles, radial velocities are found for the corresponding depths according to the formula:

q-V q-v

))

(11)(eleven)

где q - коэффициент пересчета, прини- маемьй равным 0,98 дл  осадочного разреза и 0,95 дл  консолидированного разреза. Затем определ ют глубину дл  точек отражений, соответствующих середине годографов отраженных воли. Дл  поперечных профилей глубина определ етс  по формулеwhere q is a conversion factor, the acceptance is 0.98 for the sedimentary section and 0.95 for the consolidated section. Then, the depth is determined for the reflection points corresponding to the mid-hodograph of the reflected will. For transverse profiles, the depth is determined by the formula

Z попер 5 (tVд, )-(Z popper 5 (tVd,) - (

дл  продольных профилей поfor longitudinal profiles

Znpoq I -(t V,Znpoq I - (t V,

луч Ray

)2-Х2 ,) 2-X2,

(12)(12)

формулеformula

(13)(13)

30thirty

3535

4040

4545

5050

5555

где t - врем  прихода отраженной волны в середине годографа;where t is the arrival time of the reflected wave in the middle of the hodograph;

у-рассто ние от проекции соответствующего источника возбуждени  на поперечньш профиль до приемника где зарегистрирован приход отраженного луча, соответствующего середине годографа;a y-distance from the projection of the corresponding excitation source on the transverse profile to the receiver where the arrival of the reflected beam corresponding to the middle of the locus is registered;

R - рассто ние по перпендикул ру от соответствующего источника возбуждени  до линии поперечного профил ;R is the perpendicular distance from the corresponding excitation source to the transverse profile line;

X - рассто ние на продольном профиле , отсчитываемое от источника возбуждени , где зарегистрирован приход отраженного луча, соответствующего середине годографа.X is the distance on the longitudinal profile, measured from the excitation source, where the arrival of the reflected beam, corresponding to the middle of the hodograph, is recorded.

Найденные дл  соответствующих глубин лучевые скорости по данным поперечных и продольных профилей системы порознь представл ют в виде вертикального пол  лучевых скоростей. Дл  дальнейшего исследовани  вертикальные пол  лучевых скоростей усредн ют . С этой целью используют скольз щее по каждому из вертикальных полей лучевых скоростей заданное окно, величина которого определ етс  изThe radial velocities found for the corresponding depths according to the transverse and longitudinal profiles of the system are separately represented as a vertical radial velocity field. For further study, the vertical radial velocity fields are averaged. For this purpose, a given window is slid along each of the vertical radial-velocity fields, the value of which is determined from

51А27351A273

соотношени  (6), которое при прин тыхrelations (6), which, when received

услови х равно Д2conditions equal to D2

йХyx

0,0750.075

,05 км. При каждом усреднении за- данное окно перемещают по глубине.на его величину. Наход т усредненные значени  лучевых скоростей Удцц„g. соответственно по системам поперечных и продольных профилей, которые представл ют в виде графиков изменени  усредненных значений лучевых скоростей с глубиной, а также наход  соответствующие им вертикальные времена по формуле:, 05 km. With each averaging, the given window is moved in depth. To its value. Find the averaged values of the radial velocities of Udc „g. respectively, along the systems of transverse and longitudinal profiles, which are represented in the form of graphs of changes in averaged values of radial velocities with depth, and also found the corresponding vertical times by the formula:

(U)(U)

.p.p

где Z-глубина 5 соответствующа  усредненному значению лучевых скоростей.where Z-depth 5 corresponds to the average radial velocity value.

Найденные вертикальные времена представл ют в виде графиков изменени  по глубине вертикальных времен распространени  отраженных :волн по которым стро т графики распределени  по глубине пластовых скоростей. Затем наход т отношениеThe found vertical times are presented in the form of graphs of changes in the depth of the vertical propagation times of the reflected: waves in which graphs of the distribution over the depth of reservoir velocities are plotted. Then find the relationship

1C 1C

пл. проА,.С«ет, square proA, .C "em,

(15)(15)

VV

пл. попер, сист.square poper syst

характеризующее коэффициент неоднородности пласта. Значение коэффициента ТС представл ют в виде графика изменени  коэффшУиента неоднородно сти пласта с глубиной. Пласты считаютс  однородными, если % 1, и неоднородными, если К 5 1.characterizing the coefficient of heterogeneity of the reservoir. The value of the TC coefficient is presented in the form of a graph of changes in the coefficients of the inhomogeneity of the formation with depth. The strata are considered homogeneous if% 1, and heterogeneous if K 5 1.

Положительной эффект достигаетс  за счет регулировани  степени деталь ности изучени  геологического разрез при многократном накапливании информации .A positive effect is achieved by adjusting the degree of detail of studying a geological section with multiple accumulations of information.

Применение способа позвол ет производить вертикальное зондирование, дающее возможность прогнозировать геологический разрез с необходимой точностью до любой заданной глубиныThe application of the method allows vertical sounding, making it possible to predict a geological section with the necessary accuracy to any given depth.

с выделением пластов с тюпилче.ипьши и повьш1енными скорост ми.with the release of layers with typilch.ipshi and higher speeds.

Claims (1)

Формула изобретени Invention Formula Способ сейсмозондировани  дл  прогнозировани  геологического разре™ за, включающий возбуждение упругих колебаний в р де пунктов,располо- женньк вдоль профил , регистрацию колебаний в пунктах приема, образующих систему приема, состо щую из отрезков совмещенных продольных и поперечных профилей, при этом систему приема и источник возбухсдени  пос- ледовательно перемещают по профилю в противоположных направлени х от центра исследуемой области с одинаковым шагом перемещени  Л Х, о т- лич ающий с  тем, что, с целью повышени  де тельности и точности изучени  геологического, разрезаA seismic sounding method for predicting a geological section, including the excitation of elastic vibrations at a number of points, located along a profile, recording vibrations at reception points that form a reception system consisting of segments of combined longitudinal and transverse profiles, the receiving system and the source airways are successively moved along the profile in opposite directions from the center of the study area with the same displacement step LH, which is so that, in order to increase the activity the accuracy of the geological, section В пределах площадной системы наблюдений , ограничивающей исследуемую область,минимально допустимую длину L каждого отрезка поперечной расстановки , отсчитываемую в одну сторону от проекции на дневную поверхность источника возбуждени , принимают равной длине отрезка продольной расстановки , определ емой из услови  0,08 ,, L 0,8 V,,. д,„,, , где Vf,. плакс наибольша  пластова  скорость, ожидаема  в изучаемом интервале разреза , км/с, а шаг перемещени  Д X определ ют из услови  а X Within an area observation system limiting the area under study, the minimum allowable length L of each transverse spread segment, counted in one direction from the projection onto the day surface of the excitation source, is taken to be equal to the length of the longitudinal spread segment, determined from the condition 0.08 ,, L 0, 8 V ,,. d, „,, where Vf ,. the weeping layer has the highest formation velocity expected in the section interval under study, km / s, and the displacement step D X is determined from the condition X п„.„н . где, Vn,,..;, - наименьша  пластова  скорость разреза, км/с. По каждой системе продольных и поперечных расстановок наход т распределение усредненных.значений лучевых скоростей и соответствующих им вертикальных времен, по которым определ ют пластовые скорости, и по величине и характеру их изменени  с глубиной суд т о степени неоднородности пласта и прогнозируют геологический разрез.  n „.„ n. where, Vn ,, ..;, is the smallest formation velocity of the cut, km / s. For each system of longitudinal and transverse arrangements, the distribution of averaged values of radial velocities and the corresponding vertical times for which reservoir velocities are determined is found, and the magnitude and nature of their changes with depth judge the degree of heterogeneity of the reservoir and predict a geological section. 2 .S2 .S
SU864123270A 1986-06-26 1986-06-26 Method of seismic probing for predicting geological section SU1427312A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
SU864123270A SU1427312A1 (en) 1986-06-26 1986-06-26 Method of seismic probing for predicting geological section

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
SU864123270A SU1427312A1 (en) 1986-06-26 1986-06-26 Method of seismic probing for predicting geological section

Publications (1)

Publication Number Publication Date
SU1427312A1 true SU1427312A1 (en) 1988-09-30

Family

ID=21258772

Family Applications (1)

Application Number Title Priority Date Filing Date
SU864123270A SU1427312A1 (en) 1986-06-26 1986-06-26 Method of seismic probing for predicting geological section

Country Status (1)

Country Link
SU (1) SU1427312A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101403797B (en) * 2008-11-14 2011-03-23 北京市市政工程研究院 Advanced geological prediction system and method for underground engineering construction

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101403797B (en) * 2008-11-14 2011-03-23 北京市市政工程研究院 Advanced geological prediction system and method for underground engineering construction

Similar Documents

Publication Publication Date Title
CA2594339C (en) Method and apparatus for estimating formation slowness
US7310580B2 (en) Method for borehole measurement of formation properties
US7639563B2 (en) Method for sonic indication of voids in casing cement
Harrison et al. Acquisition and analysis of sonic waveforms from a borehole monopole and dipole source for the determination of compressional and shear speeds and their relation to rock mechanical properties and surface seismic data
US10451475B2 (en) Gauge length optimization in distributed vibration sensing
CN102109617A (en) Method for measuring Q value of near surface strata by using twin-well microlog
US5197038A (en) Method and sonic tool apparatus for investigating properties of earth formations transversed by a borehole
CN101178439A (en) Accurate ground surface investigation method
US8995224B2 (en) Real-time velocity and pore-pressure prediction ahead of drill bit
JPS62269084A (en) Seismic exporation method having extremely high resolution in horizontal well
CN112558180B (en) Method for rapidly checking seismic horizon calibration accuracy by utilizing horizontal isochronal plane
CN109541689A (en) A kind of medium compactness evaluation method based on reflected energy feature
CN1052530A (en) Pore pressure prediction method
Huang et al. Fast-forward modeling of compressional arrival slowness logs in high-angle and horizontal wells
RU2336541C2 (en) Method of low frequency seismic probing for research and survey of hydrocarbon deposits (versions)
SU1427312A1 (en) Method of seismic probing for predicting geological section
RU2305298C1 (en) Method for controlling position of drilling bit during drilling
RU2024891C1 (en) Method of seismic prospecting
Chudy et al. Cased hole acoustic logging-a solution to a problem
RU2201606C1 (en) Method of typification and correlation of oil and gas productive rocks by borehole spectral-time parameters
SU873185A1 (en) Method of oil and gas deposit prospecting
RU2339978C1 (en) Method of borehole sesmic exploration works
SU966634A1 (en) Seismic prospecting method
Mari et al. Tying surface seismic data6
RU2305856C1 (en) Land-well seismology method