SU1689900A1 - Method of seismic prospecting - Google Patents
Method of seismic prospecting Download PDFInfo
- Publication number
- SU1689900A1 SU1689900A1 SU894686020A SU4686020A SU1689900A1 SU 1689900 A1 SU1689900 A1 SU 1689900A1 SU 894686020 A SU894686020 A SU 894686020A SU 4686020 A SU4686020 A SU 4686020A SU 1689900 A1 SU1689900 A1 SU 1689900A1
- Authority
- SU
- USSR - Soviet Union
- Prior art keywords
- wave
- reflected
- seismic
- svp
- waves
- Prior art date
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Изобретение относитс к сейсморазведке и может быть использовано при проведении нефтегазопоисковых работ методами многоволновых исследований. Цель изобретени - повышение эффективности многоволновой сейсморазведки на отраженных обменных волнах SvP. Дл этого по способу сейсмической разведки, по которому возбуждают продольные Р-и поперечные Sv-волны, регистрируют отраженные сейсмические колебани вертикальными сейсмоприемниками на линейных профил х и полученные сейсмограммы с общей точкой возбуждени (ОТВ) подвергают последующей обработке, регистрацию сейсмических колебаний провод т с выносом В, удовлетвор ющим условию 0,,8Н, и на базе Б, длина которой лежит в пределах 0,, где Н - рассто ние по нормали от пункта возбуждени до исследуемой границы . При обработке суммируют с использованием оператора подоби в качестве доа суммировани трассы указанной базы по пр мой линии с наклоном, совпадающим с наклоном годографа ОТВ-РР При этом пр молинейное суммирование начинают с времени вступлени отраженной РР-гюлны от исследуемой границы и периодически повтор ют с временными сдвигами Наход т абсолютный максимум оператора подоби , по которому идентифицируют SvP-волну от исследуемой границы, и по параметрам выделенной SvP-волны, комплексиру их с, данными от других волн, суд т о строении и вещественном составе геологического разреза . 2 ил. Г L О 00 юThe invention relates to seismic exploration and can be used in oil and gas exploration using multi-wave research methods. The purpose of the invention is to increase the efficiency of multi-wave seismic survey on reflected SvP exchange waves. To do this, the seismic survey method, which excites longitudinal P and transverse Sv waves, records the reflected seismic vibrations by vertical seismic receivers on linear profiles and the resulting seismograms with a common excitation point (TSS) are subjected to subsequent processing, the seismic vibrations are recorded with the removal B, satisfying the condition 0,, 8Н, and on the base B, the length of which lies within 0, where Н is the normal distance from the excitation point to the boundary under study. During processing, they summed using the similar operator as the additional summation of the path of the indicated base along a straight line with a slope that coincides with the slope of the OTV-PP hodograph. In this case, the linear linear summation starts from the arrival time of the reflected PP-radius from the studied boundary and periodically repeated with temporary shifts Find the absolute maximum of the similar operator, which is used to identify the SvP wave from the boundary under study, and the parameters of the selected SvP wave, complex them with data from other waves, are judged by and the material composition of the geological section. 2 Il. G L O 00 th
Description
Изобретение относитс к способам сейсмической разведки и может быть использовано при проведении многоволновых исследований с применением отраженных обменных волн SvP.The invention relates to seismic survey methods and can be used when conducting multi-wave surveys using reflected SvP exchange waves.
Целью изобретени вл етс повышение эффективности многоволновой сейсморазведки на отраженных обменных волнахThe aim of the invention is to increase the efficiency of multi-wave seismic survey on reflected exchange waves.
SvP.SvP.
На фиг.1 приведена схема осуществлени способа; на фиг.2 - результаты моделировани с оператором подоби в качестве дра суммировани .Figure 1 shows the scheme for implementing the method; Fig. 2 shows simulation results with a similarity operator as the sum kernel.
На фиг. 1 и 2 показаны дневна поверхность 1, пункт2 возбуждени , группы 3.1-З.п сейсмоприемников, соединенные сейсмической косой с бортовой станцией 4 котора св зана с ЭВМ 5, исследуема отражающа граница 6, луч 7 падающей волны, луч8отраженной волны, вынос 9 пункта 2 возбуждени база 10 размещени пунктов наблюдени , рассто ние 11 по нормали от пункта возбуждени до исследуемой границы 6, зависимость 12 коэффициента подоби от времени.FIG. Figures 1 and 2 show the daytime surface 1, excitation point 2, groups 3.1-З.п of seismic receivers connected by a seismic oblique to an onboard station 4 which is connected to the computer 5, the reflecting boundary 6 under study, the beam 7 of the incident wave, the reflected reflected wave, the carry-out of 9 point 2 excitation base 10 of location of observation points, distance 11 along the normal from the point of excitation to the boundary under study 6, dependence of the similarity factor on time.
Способ осуществл етс следующим об разом.The method is carried out as follows.
оabout
т-&t- &
Провод тс полевые работы по типовой дл исследуемого геологического района технологии метода многократных перекрытий на отраженных продольных волнах, при этом система наблюдений должна удовлетвор ть следующим услови м: вынос равен 0,,8Н, длина базы регистрации сейсмограммы с общей точкой возбуждени (ОТВ)-0,,0 Н, где Н - рассто ние по нормали от пункта возбуждени до исследуемой границы.The field work is carried out according to the method of multiple overlaps on reflected longitudinal waves, which is typical for the geological region under study, and the observation system must satisfy the following conditions: the offset is 0, ..., 8H, the base of the seismogram with a common excitation point (TSS) -0 ,, 0 Н, where Н is the normal distance from the point of excitation to the boundary under study.
Сущность способа основана на физических свойствах отраженных обменных волн SvP. Коэффициент отражени со сменой направлени пол ризации волн при малых углах падени на границу практически равен нулю, поэтому интерференционна система приема обменных волн должна быть фланговой с выносом, В этом случае прослежи- ваемость обменных волн SvP увеличитс за счет снижени интерференции обменных и монотипных волн в плоскости x-t(x- рассто ние между источником и приег ником, t - врем ),The essence of the method is based on the physical properties of the reflected exchange waves SvP. The reflection coefficient with a change in the direction of polarization of the waves at small angles of incidence on the boundary is almost zero, therefore the interference system of receiving exchange waves should be flank with the removal. In this case, the traceability of SvP exchange waves will increase due to the reduction of interference of exchange and monotypic waves in the plane xt (x is the distance between the source and the receiver, t is the time),
Кроме того, годограф ОТВ - SvP раньше выходит на асимптоту, чем годограф ОТВ-РР, т.е. первый при меньшем выносе и, следовательно , на большей базе аппроксимируетс пр мой линией. Таким образом, выделение SvP-волны целесообразно осуществл ть посредством суммировани по пр мой линии трасс ОТВ. Дл этого необходимо определить оптимальные величины выноса и базы приема, при которых волны SvP можно считать плоскими, а пр молинейное суммирование - эффективным способом их выделени . Данную задачу можно решать с помощью численного моделировани , выполненного на синтетических сейсмограммах, полученных дл различных эффективных моделей сейсмического разреза. На сейсмограммах при этом присутствуют как монотипные отраженные волны обоих типов, так и обменные отраженные обоих типов. Результаты моделировани показывают, что на выносе, лежащем в пределах 0,7-0,8-Н, и на базе в пределах 0,8:1,0 Н, волну SvP, отраженную от субгоризонтальной границы, можно считать плоской, распростран ющейс вдоль линейного профил с кажущейс скоростью VSP , равной по величине кажущейс скорости Vpp отраженной продольной волны РР, порожденной на исследуемой границе и регистрируемой на той же базе, что находитс в пр мом соответствии с теорией о свойствах годографов SvP и РР в дальней зоне.In addition, the locus OTV - SvP earlier comes to the asymptote than the locus OTV-PP, i.e. the first with a smaller offset and, therefore, on a larger base is approximated by a straight line. Thus, it is advisable to extract the SvP-wave by summing along the direct paths of the GSP. To do this, it is necessary to determine the optimal values of the offset and the base of reception, at which the SvP waves can be considered flat, and the linear summation - an effective way to isolate them. This problem can be solved using numerical simulation performed on synthetic seismograms obtained for various effective seismic section models. At the same time, on the seismograms there are both monotypic reflected waves of both types and exchange reflected of both types. The simulation results show that on the offset, lying in the range of 0.7-0.8-H, and on the base in the range of 0.8: 1.0 N, the SvP wave reflected from the subhorizontal boundary can be considered flat, propagating along a linear profile with an apparent velocity VSP equal to the apparent velocity Vpp of the reflected longitudinal wave PP generated at the boundary under study and recorded on the same base, which is in direct accordance with the theory of the SvP and PP properties in the far zone.
Совпадение наклонов годографов SvP позвол ет избежать при выделении SvPволны перебора по наклону линии суммировани , если известен наклон, иначе говор Vpp годографа волны РР. Отраженна SvP- волна приходит позже, чем отраженна РРThe coincidence of the slopes of the SvP hodographs makes it possible to avoid the selection of the SvPwave brute force by the slope of the summation line, if the slope is known, otherwise speaking Vpp of the hodograph of the PP wave. The reflected SvP wave arrives later than the reflected PP
от того же горизонта, поэтому поиск и выделение первой начинаютс пос ле времени вступлени РР-волны.from the same horizon, therefore, the search and selection of the first one begin after the time of entry of the PP wave.
Другим существенным физическим свойством SvP-волн, отличающим их от монотипных и обменных волн PSv, вл етс значительно меньший размер отражающей площадки при одинаковых размерах приемных апертур. Следовательно, Шероховатость отражающей границы приводит к меньшимAnother significant physical property of SvP waves, which distinguishes them from the monotypic and exchangeable PSv waves, is the significantly smaller size of the reflecting area with the same sizes of receiving apertures. Consequently, the Roughness of the reflecting boundary leads to less
по величине фазовым и амплитудным разбросам между наблюдаемым и аппроксимирующим годографами, чем и увеличиваетс отношение сигнал/помеха при выделении SvP-волн посредством пр молинейногоlargest phase and amplitude spreads between the observed and approximating hodographs, which increases the signal-to-noise ratio when SvP waves are emitted by means of a straight line
суммировани , а также повышаетс точность прив зки волны к отражающему горизонту .the summation, as well as the accuracy of anchoring the wave to the reflecting horizon.
Конкретные значени параметров системы приема задаютс исход из требованийSpecific values of the parameters of the reception system are based on the requirements
геологической задачи, поверхностных условий , канальности сейсмостанции и т.д.geological problem, surface conditions, channel seismic station, etc.
По волнам РР определ етс посредством стандартной цифровой обработки, включающей расчет и коррекцию кинематических поправок, врем tpp прихода волны РР, отраженной от исследуемой границы в начале базы регистрации, и ее кажуща с скорость распространени Vpp .Over the waves, PP is determined by standard digital processing, including the calculation and correction of kinematic corrections, the time tpp of the arrival of the PP wave reflected from the boundary under study at the beginning of the registration base, and its apparent speed Vpp.
В последующем работа по способу заключаетс в обработке сейсмограмм ОТВ- SvP на ЭВМ по специальной программе, в которой реализуема алгоритм поиска отраженной обменной волны SvP.Subsequently, the work according to the method consists in processing the OTV-SvP seismograms on a computer using a special program that implements an algorithm for searching for the reflected exchange wave SvP.
Сущность алгоритма следующа . Задают входные данные: врем прихода tpp и кажущуюс скорость РР волны Vpp , определенные на этапе обработки сейсмограмм ОТВ-РР, интервал суммировани ДТ по координате t, который рекомендуетс братьThe essence of the algorithm is as follows. Input data is entered: arrival time tpp and apparent velocity PP of the Vpp wave, determined at the processing stage of the HTS-PP seismograms, DT summation interval along the coordinate t, which is recommended to take
равным по величине временному интервалу , содержащему основную энергию волнового импульса и составл ющему примерно 100 мс в среднечастотной сейсморазведке. Выбирают интервал поиска волн SvPAequal in magnitude to the time interval containing the main energy of the wave pulse and approximately 100 ms in mid-frequency seismic exploration. Select the SvPA wave search interval.
по координате t, Рекомендуемое значение Аравноon t, Recommended Value
Atpp(1+y)/2y , Atpp (1 + y) / 2y,
где Atpp-средн разность времен прослеживани целевых границ на временном разрезе , построенном на отраженных продольных РР-волнах;where Atpp is the average difference in the tracking time of the target boundaries in a time section built on reflected longitudinal PP waves;
Y Vs/Vp - априорное значение отношени скоростей поперечных и продольныхY Vs / Vp - a priori value of the ratio of the transverse and longitudinal velocities
волн дл исследуемого объекта, если оно неизвестно, то у 0,5.waves for the object under study, if it is unknown, then about 0.5.
Вычисл ют в интервале поиска Дс шагом At оператор подоби П( г,т.) по формулеCalculate in the search interval Ds by step At the operator of the similarity P (r, t.) Using the formula
n(,t) n (, t)
Г:((х;)}G: ((x;)}
.)2.) 2
к - (/ 1 L- to - (/ 1 L-
где АХ - количество трасс, участвующих в суммировании;where AH is the number of runs involved in the summation;
АХ - рассто ние между каналами;AH is the distance between the channels;
r(Xi)Xi/Vpp - приращение времени по линии анализа в точке ,7V+i AX;r (Xi) Xi / Vpp is the increment of time along the line of analysis at a point, 7V + i AX;
а AT/(5t - параметр суммировани , задающий число отсчетов трассы по оси t;and AT / (5t is the summation parameter specifying the number of samples of the trace along the t axis;
AT - окно суммировани ;AT is a summation window;
5t - шаг дискретизации по t;5t is the sampling step in t;
A(tK) - амплитуда трассы в момент времени IK;A (tK) is the amplitude of the path at the moment of time IK;
to tpp+At, tpp - врем вступлени волны РР от исследуемого горизонта;to tpp + At, tpp is the time of entry of the PP wave from the studied horizon;
At - шаг перебора по оси t.At - step search in the t axis.
Вычисл ют оператор подоби Цт.о ) в точках (l-1) A/At- 5t, 1,2,..чпри заданных в пункте значени tpp, Vpp , AT и выбранных из этикетки сейсмограммы ОТВ - SvP AX, 5t. Шаг перебора At рекомендует задавать не больше 0.5 ЛТ.Calculate the operator of the similarity of Ts.o) at the points (l-1) A / At-5t, 1,2, .. when the values of tpp, Vpp, AT specified at the point and the values of OTV selected from the label of the TOV - SvP AX, 5t. The iteration step At recommends setting no more than 0.5 LT.
Определ етс врем , на которое дости- гартс абсолютный максимум коэффициента подоби The time for which the absolute maximum of the coefficient of similarity is reached is determined.
Ц1м3кс)тах(П(т.о )), (i-1) A/ At.C1m3x) max (P (to)), (i-1) A / At.
Во временном интервале, содержащем 1макс, например1макс-А1 То 1макс+Дтс шагом At (2-3)(5t, вновь вычисл етс Цт.0), to -Хмакс- j AT/, j,0,1,2,.,. и определ етс врем , на котором коэффициент подоби достигает абсолютного максимумаIn the time interval containing 1 max, for example, 1 max-A1, then 1 max + Dts, increments of At (2–3) (5t, CT0 is again calculated), to -Xmax-j AT /, j, 0,1,2,., . and determine the time at which the similarity coefficient reaches the absolute maximum
IXtsv р)тах(П(т0 )),IXtsv p) max (P (t0)),
На выходе программы получают численное значение tsv P. фиксирующее в начальной точке приема врем регистрации волны SvP, отраженной от исследуемой границы.At the output of the program, a tsv P numerical value is obtained. The time of registration of the SvP wave reflected from the boundary under study is fixed at the initial reception point.
На фиг.2 приведены результаты выделени волны SvP по указанному и программно реализованному алгоритму. Полученные значени tsv P) используютс дл типовых расчетов эффективных интегральных параметров Vs, J Vs/Vp, которые используютс в качестве начальных данных дл последующей глубокой обработки, дающей возможность оценить литологию геологического разреза с позиции наличи углеводородов с большей веро тностью, чем метод ркого п тна.Figure 2 shows the results of the SvP wave extraction by the indicated and software-implemented algorithm. The tsv P values obtained are used for typical calculations of the effective integral parameters Vs, J Vs / Vp, which are used as initial data for subsequent deep processing, making it possible to estimate the lithology of the geological section from the position of the presence of hydrocarbons more likely than the bright spot. .
Использование предлагаемого способа сейсмической разведки по сравнению с известным обеспечивает повышение эффективности исследований на отраженных обменных 5 волнах SvP благодар следующим примуще- ствам Повышаетс помехоустойчивость и, следовательно, информативность исследований ча SvP-волнзх, что обусловлено опти0 мальными (дл этих волн) параметрами системы наблюдени и способом выделени данных волн. Кроме того, в отличие от известного осуществл етс прив зка SvP- БОЛН к исследуемому объекту и, с едова5 тельно, повышаетс точность построени параметров (например, Vs, у), характеризующих литологию и структуру геологического разреза Дл реализации предлагаемого „пособа не требуетс специальный источС ник ипа кругова сила, а достаточно источника поперечных волн Sv.The use of the proposed seismic survey method as compared to the known method provides an increase in the efficiency of research on reflected SvP 5 waves due to the following advantages: Improves the noise immunity and, therefore, the informativity of the SvP waves research, which is due to the optimal (for these waves) parameters of the observation system and a way to isolate these waves. In addition, in contrast to the well-known, SvP-BOLN is attached to the object under study and, more precisely, the accuracy of the construction of parameters (for example, Vs, y) that characterize the lithology and structure of the geological section is improved. Nickname ip is a circular force, and a source of transverse waves Sv is sufficient.
Claims (1)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SU894686020A SU1689900A1 (en) | 1989-04-27 | 1989-04-27 | Method of seismic prospecting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SU894686020A SU1689900A1 (en) | 1989-04-27 | 1989-04-27 | Method of seismic prospecting |
Publications (1)
Publication Number | Publication Date |
---|---|
SU1689900A1 true SU1689900A1 (en) | 1991-11-07 |
Family
ID=21445099
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
SU894686020A SU1689900A1 (en) | 1989-04-27 | 1989-04-27 | Method of seismic prospecting |
Country Status (1)
Country | Link |
---|---|
SU (1) | SU1689900A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2735417C1 (en) * | 2020-01-13 | 2020-11-02 | Марк Аронович Завесин | Multi-wave seismic survey method with recording of elastic vibrations by vertical seismic receivers |
-
1989
- 1989-04-27 SU SU894686020A patent/SU1689900A1/en active
Non-Patent Citations (1)
Title |
---|
Пузырев Н.Н. О проблеме многоволновых сейсмических исследований./В кн.: Многоволновые сейсмические исследовани . - Новосибирск. Наука, 1987, с. 3-6, Авторское свидетельство СССР № 1373192. кл. G 01 V 1/00, 1986. * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2735417C1 (en) * | 2020-01-13 | 2020-11-02 | Марк Аронович Завесин | Multi-wave seismic survey method with recording of elastic vibrations by vertical seismic receivers |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2002226539B2 (en) | A method of processing marine seismic data and a method of seismic surveying | |
Coppens | First arrival picking on common‐offset trace collections for automatic estimation of static corrections | |
US7408836B2 (en) | Method of and apparatus for processing seismic data in the presence of anisotropy | |
EP0990174B1 (en) | High resolution determination of seismic polar anisotropy | |
US4635238A (en) | Data processing method for correcting P and S wave seismic traces | |
US20100223012A1 (en) | Method for Determination of Diffractor Locations at Sea Bottom for the Purpose of Attenuating Such Energy | |
RU2344444C2 (en) | Seismic data processing | |
US20030048696A1 (en) | Trapped water bottom multiple and peg-leg multiple suppression for ocean bottom seismic data | |
AU2002226539A1 (en) | A method of processing marine seismic data and a method of seismic surveying | |
US3696331A (en) | Automated process for determining subsurface velocities from seismograms | |
Stoffa et al. | Velocity analysis for wide aperture seismic data | |
US6430508B1 (en) | Transfer function method of seismic signal processing and exploration | |
EP0866343A2 (en) | Method for determining seismic velocities | |
US4397004A (en) | Method for seismic exploration by vertical seismic profiling and installation for its implementation | |
EA005232B1 (en) | A method and apparatus for processing seismic data | |
US4101867A (en) | Method of determining weathering corrections in seismic operations | |
Duckworth et al. | Crustal structure measurements near FRAM II in the Pole Abyssal Plain | |
US3651451A (en) | Continuous velocity estimation | |
US8144543B2 (en) | Seismic data processing method for RMO picking | |
SU1689900A1 (en) | Method of seismic prospecting | |
AU2002216298B2 (en) | A seismic method and apparatus for generating a semblance panel and computing the reflector dip | |
Basham et al. | An analysis of the “Benham” aftershock, sequence using Canadian recordings | |
Zimmerman et al. | Comparison of vertical seismic profiling techniques | |
Flinn et al. | Detection and analysis of multiple seismic events | |
Gaiser et al. | VSP fundamentals that improve CDP data interpretation |