RU2116657C1 - Способ обработки сейсмических данных для подавления многократных переотражений - Google Patents
Способ обработки сейсмических данных для подавления многократных переотражений Download PDFInfo
- Publication number
- RU2116657C1 RU2116657C1 RU95111130A RU95111130A RU2116657C1 RU 2116657 C1 RU2116657 C1 RU 2116657C1 RU 95111130 A RU95111130 A RU 95111130A RU 95111130 A RU95111130 A RU 95111130A RU 2116657 C1 RU2116657 C1 RU 2116657C1
- Authority
- RU
- Russia
- Prior art keywords
- characteristic
- sign
- speed
- specified
- pressure
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000012545 processing Methods 0.000 title claims description 10
- 230000001629 suppression Effects 0.000 claims abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 24
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 230000003044 adaptive effect Effects 0.000 claims description 9
- 238000005311 autocorrelation function Methods 0.000 claims description 8
- 239000002245 particle Substances 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 2
- 239000012736 aqueous medium Substances 0.000 claims 1
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 230000006835 compression Effects 0.000 abstract description 14
- 238000007906 compression Methods 0.000 abstract description 14
- 238000004364 calculation method Methods 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 3
- 238000011835 investigation Methods 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 230000000875 corresponding effect Effects 0.000 description 10
- 230000004044 response Effects 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- 238000004880 explosion Methods 0.000 description 4
- 239000002689 soil Substances 0.000 description 4
- 238000002955 isolation Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3808—Seismic data acquisition, e.g. survey design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/42—Devices characterised by the use of electric or magnetic means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/56—De-ghosting; Reverberation compensation
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Measuring Volume Flow (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Использование: для ослабления эффектов реверберации на сейсмические волновые поля с высоким уровнем помех, записанных изнутри водонасыщенной среды. Сущность изобретения: волновое поле появляется в виде характерного признака сжатия и характерного признака скорости, в котором содержится характерный признак фонового шума. Признак сжатия адаптивно фильтруется и субтрактивно (в режиме вычитания) комбинируется со скоростным признаком, в результате чего удается выделить почти в чистом виде характерный признак помех. Этот последний почти в чистом виде вновь складывается с первоначальным характерным признаком скорости, но с противоположным знаком, чтобы исключить влияние фонового шума и получить уточненный и очищенный характерный признак скорости. Производится итерационный пересчет уточненного скоростного признака и его суммирование с характерным признаком сжатия при постоянном увеличении коэффициента пересчета при каждой итерации и автокорреляции полученной суммы. После каждой операции суммирования выполняется расчет коэффициента сходимости, и полученная в результате величина коэффициента сходимости, максимально близкая к единице, определяет оптимальный коэффициент пересчета. 9 з.п. ф-лы, 6 ил.
Description
Изобретение относится к усовершенствованному способу обработки сейсмических данных, полученных от спаренных, но принципиально разнотипных сейсмоприемников, обеспечивающему подавление интерференции, возникающей вследствие реверберации в толще воды, при проведении морских сейсмических съемок на мелководных участках.
В ходе морских сейсмических съемок используют буксирование сквозь толщу воды на длинном кабеле приемников звуковых волн. Периодически источник звуковых волн вызывает ответное звучание находящихся под поверхностью воды слоев грунта, образуя отраженное поле сейсмических волн, которое фиксируется сейсмоприемниками и преобразуется в электрические сигналы, передаваемые по кабелю на судно. Сигналы записываются, а затем обрабатываются с целью визуального отображения топографии тех слоев грунта на морском дне, на которые направлены звуковые волны.
В морской среде отраженные волновые поля проявляются в виде колебаний гидродинамического давления или колебания скоростей частиц. Сейсмоприемники могут принадлежать к принципиально разным типам аппаратуры. Говоря иначе, они могут предназначаться для измерения колебаний давления (гидрофоны) или колебаний скорости частиц (геофоны). В зависимости от характера задач, которые решаются с помощью полевых исследований, на одном и том же кабеле могут находиться сейсмоприемники обоих типов либо какого-то одного из них. Реже в аналогичных условиях применяются акселерометры.
На относительных мелководьях (при глубинах порядка 8-63 м [25-200 футов] целесообразно применять так называемые "донные сейсмические косы". При этом методе сейсмоприемная коса размещается непосредственно на морском дне, а не буксируется постоянно вслед за судном, как принято на глубоководных участках моря. Оптимальный сейсмоприемной аппаратурой для измерения скоростей частиц на морском дне считаются геофоны, установленные на универсальном шарнирном устройстве. Речь идет о развертывании на дне косы из сотен или тысяч сейсмоприемников, электрические выходные сигналы которых подаются в мультиплексном режиме в соответствующие устройства регистрации данных.
На фиг. 1 показано судно 10 для укладки донного кабеля 12, который в данном случае находится на дне водоема 16 под толщей воды 14. Известно, что поверхность воды 18 образует в таких случаях отражающую сигналы границу раздела воздух-вода. Морское дно является также хорошим отражателем, хотя это свойство зависит от величины акустического импеданса. Датчик скорости (геофон) 20 и датчик давления (гидрофон) 22 размещены рядом на морском дне и соединены с различными электрическими каналами кабеля 12, через которые их сигналы поступают в установленную на судне 10 аппаратуру архивного запоминания и обработки данных (не показана). В целях упрощения изложения здесь показаны лишь два сейсмоприемника. Периодически источник 13 создает поле звуковых волн 15, которое распространяется в донный грунт 17, где происходит его последующее отражение от подповерхностных слоев и возвращение в виде поля отраженных сигналов, такого как поле 24. Источник 13 может быть запущен с борта судна 10 или с использованием отдельной морской сейсморазведочной станции (не показана).
Приходящее поле 24 отраженных волн сжатия достигает установленного на дне геофона 20, который фиксирует его приход. В результате, согласно общепринятой в сейсморазведке эталонной схеме, формируется нарастающей электрический импульс 26 (фиг. 2). Отраженная волна продолжает восходящее движение, пока не достигнет границы раздела воздух - вода 18, откуда после сдвига фазы на 180o снова отразится вниз и уже в виде побочной волны попадет на геофон сверху. При нормальном ходе событий импульс сжатия, воздействующий на верхнюю часть геофона, создает уменьшающийся электрический заряд, однако с учетом обращения фазы на водной поверхности первая побочная отраженная волна формирует второй нарастающий электрический импульс 30, приведенный на фиг. 2, относительно произвольной шкалы амплитуд. В зависимости от величины акустического импеданса морского дна и сглаженности поверхности воды поля звуковых волн могут многократно переотражаться в прямом и обратном направлениях (реверберировать) между поверхностью и дном, подобно многократным переотражениям между зеркалами на противоположных стенах помещения. Вторая отраженная побочная волна формирует отрицательный импульс 32. Последующие импульсы чередуются по знаку, принимая то положительную, то отрицательную полярность. Реверберационные явления наблюдаются в водонасыщенных средах, например в прибрежном песке при низком приливе или в зыбучих песках.
Гидрофон 22 воспринимает переходящее отраженное поле волн сжатия и фиксирует его приход в виде импульса сжатия. И снова, по общепринятому отраслевому стандарту, гидрофон преобразует импульс сжатия в нарастающий электрический импульс такой, как импульс 36 на фиг. 3, изображенный относительно той же самой произвольной шкалы амплитуд, что и на фиг. 2. Первая отраженная от поверхности побочная волна движется по направлению вниз в виде импульса разрежения, который улавливается гидрофоном 22 как электрический импульс 40 отрицательной полярности. Благодаря второму отражению от дна к поверхности, вторая отраженная побочная волна имеет положительный знак, а последующие волны будут поочередно иметь то положительную, то отрицательную полярности.
Интервал времени между импульсами τ определяется произведением удвоенной кажущейся глубины водной толщи на величину, обратную скорости воды (1/v). При нормальном падении, которое имеет место в точках 28 и 38 (фиг. 1), кажущаяся глубина воды является ее истинной глубиной. Во всех остальных случаях, например для волнового поля 44 и пути 45, кажущаяся глубина представляет собой произведение удвоенной глубины водной толщи и обратной величины скорости воды, деленное на косинус угла падения α .
Таким образом, геофоны фиксируют поступление первой отраженной волны и первой побочной волны в виде электрических импульсов одинаковой полярности, в то время как гидрофоны воспринимают поступление первой отраженной волны и первой побочной волны в виде разнополярных электрических импульсов. Поэтому использование двух спаренных разнотипных сейсмоприемников при сочетании их отличительных характерных признаков должно создавать условия, позволяющие усиливать энергию первой приходящей волны и подавлять последующие кратные волны реверберации. Второе преимущество основано на том, что случайные помехи, воспринимаемые разнотипными сейсмоприемниками, не обязательно должны быть одинаковыми либо коррелируемыми.
Для успешного совмещения сигнала гидрофона и геофона необходимо преобразовать сигналы одного из этих типов сейсмоприемников таким образом, чтобы они соответствовали сигналам сейсмоприемников другого типа. Их передаточные функции должны быть согласованы по фазе, амплитуде, частоте и затуханию.
Наиболее целенаправленным способом такого пересчета будет способ, включающий демультиплексирование набора зарегистрированных сейсмических данных от двух сейсмоприемников, устранение инструментальных эффектов усилений сигналов и уравновешивание сигналов по амплитуде с помощью идентичных функций восстановления амплитуды для обоих наборов данных. Соотношение амплитуд соответствующих первых приходящих сигналов определяет коэффициент пересчета. Однако, если учесть наличие помех и других возмущений, то данный метод представляется сверхупрощенным.
Известна электрическая сеть для согласования передаточной функции акселерометра или гидрофона с передаточной функцией геофона [1].
Известен детерминированный алгоритм расчета передаточной функции, позволяющий согласовывать гидрофоны и геофоны, компенсируя дефекты согласования этих сейсмоприемников с морским дном [2]. Такой способ требует наличия специального калибровочного волнового поля, генерируемого путем взрыва, произведенного прямо над выбранными расположенными рядом сейсмоприемниками. По времени калибровочный взрыв, как правило, непосредственно предшествует началу сейсмической съемки.
Известна компьютерная программа для получения коэффициента пересчета для съемки с использованием гидрофона, совмещенного с геофоном, с целью ослабления ревербераций в столбе воды [3]. Сейсмоприемники устанавливаются в какой-либо точке водной толщи выше уровня дна либо на самом дне. Коэффициент пересчета, который определяется из акустического импеданса воды или морского дна, можно рассчитать как детерминированно, так и статистически. В первом случае применяется калибровочный взрыв, который производится непосредственно над сейсмоприемниками, как в ранее приведенном источнике. Статистические способы не связаны с какими-либо конкретными требованиями к местоположению взрыва относительно сейсмоприемников, что составляет непременное требование при использовании детерминированного способа. Статистические способы основаны на итерационном расчете соотношений между различными сочетаниями функций автокорреляции и кросскорреляции волновых пакетов, зарегистрированных соответствующими сейсмоприемниками. С учетом разностей времен прихода волн с различными траекториями могут вноситься поправки на направленность волновых полей.
Адаптивная фильтрация помех - понятие, которое будет введено позже, объясняется в работе Бернарда Видроу и др. "Адаптивное подавление помех: принципы и применения", опублик. в "Трудах ИЕЕЕ", т. 63, N 12, декабрь 1975, с. 1692 - 1716.
Недостатки приведенных известных способов расчета коэффициента пересчета связаны, во-первых, с проблемой случайных помех. В частности, датчики скорости (геофоны) создают гораздо более высокий уровень помех по сравнению с гидрофонами. Что же касается статистических способов, то помехи геофонов неизбежно искажают результаты процессов коррелирования. Во-вторых, в итерационном процессе определения коэффициента пересчета на основе соотношений кросскорреляции и автокорреляций не содержится указания на коэффициент сходимости, который позволил бы точно установить оптимально пригодный коэффициент пересчета. Кроме того, простое суммирование данных, полученных от двух сейсмоприемников, как рекомендовано в этом источнике, обусловливает уменьшение отношения сигнал - шум для результата суммирования по сравнению с аналогичным отношением для гидрофона.
В дальнейшем изложении будут неоднократно встречаться ссылки на "скорость воды" и "характерный признак (сигнатура) скорости". При отсутствии каких-либо иных специально оговоренных значений термин "скорость" - это условный рабочий термин, означающий скорость частиц определенной среды (грунта или воды), возникающую при прохождении поля сейсмических волн через данную среду. Если же термин "скорость" будет употребляться в каком-либо ином смысле, то это его значение будет оговорено особо. Термины "характерный признак сжатия", "характерный признак помех" и "характерный признак скорости" относятся к количественным колебаниям уровня амплитуд соответствующих параметров как функции времени. Собирательный термин "характерный признак помех" или "характеристика помех" включает все и любые нежелательные сигналы, в то же время две другие отличительные характеристики "сжатие" и "скорость" относятся к искомым сейсмическим сигналам, полезным при сейсморазведке.
Таким образом, существует потребность в усовершенствованном способе статистического определения коэффициента взаимного пересчета ответных импульсов гидрофона и геофона, которых не будут искажать помехи. Кроме того, имеется потребность в способе позитивного определения коэффициента сходимости для определения оптимального коэффициента пересчета и снижения влияния помех в геофоне на суммарные данные.
В изобретении предлагается усовершенствованный способ обработки сейсмических данных с целью подавления ревербераций в столбе воды. В соответствии с предлагаемым способом в водоеме создается отраженное волновое поле, характеризующееся признаками сжатия и скорости, которые включают также фоновый шум. Характеристики сжатия и скорости обнаруживаются одновременно в выбранных соседних местоположениях. Характеристика помех геофона рассчитывается адаптивно на основе характеристик скорости и сжатия. Полученная характеристика помех суммируется с характеристикой скорости, но с противоположной полярностью, образуя не содержащую помех характеристику скорости. К такой уточненной характеристике скорости применяется выбранный коэффициент пересчета, после чего пересчитанная уточненная характеристика скорости суммируется с характеристикой сжатия. Затем выполняется автокорреляция полученной суммарной характеристики.
Для суммарной характеристики рассчитывается варимаксная функция. Перечисленные шаги повторяются несколько раз, причем величина коэффициента пересчета увеличивается при каждой итерации, пока варимаксная функция не приблизится максимально к единице.
В одном из аспектов этого изобретения на этапе выделения помех предусматривается включение операции адаптивного подавления сигнала, обусловленного фоновым шумом.
На фиг. 1 изображена иллюстрация расположения двух разнотипных сейсмоприемников на дне мелководного участка моря; на фиг. 2 - ответный сигнал датчика скорости на первичные отражения и отраженные побочные волны; на фиг. 3 - ответный сигнал датчика давления на первичные отражения и отраженные побочные волны; на фиг. 4 - схема, иллюстрирующая логику адаптивного выделения помех, а также адаптивный фильтр, действующий в направлении от гидрофона к геофону; на фиг. 5 - последовательность операций, характеризующая способ пересчета сигналов датчика скорости в сигналы датчика давления; на фиг. 6 - блок-схема предпочтительного способа осуществления обработки.
На первый взгляд может показаться, что суммирование выходных сигналов геофона с выровненными по амплитуде выходными сигналами гидрофона способно подавить многократные отражения. Однако эта проблема имеет несколько более сложное решение. Геофонам внутренне присущ высокий уровень помех по сравнению с гидрофонами. Именно этот уровень помех тем или иным способом необходимо устранить. Сигналы, показанные на фиг. 2 и 3, идеализированы в виде функций Дирака, тогда как на самом деле эти сигналы состоят из идеализированного колебания, свернутого с колебаниями переходного процесса, отображающими передаточную функцию сейсмоприемников.
В последующем изложении предполагается, что и данные представлены в форме дискретных выборок, индексированных в функции времени распространения отраженных волн.
Обработка начинается с фильтрации характеристики давления гидрофона для адаптивного получения оптимальной оценки характеристики скорости геофона. Затем оптимальную расчетную величину характерного признака скорости геофона вычитают из соответствующей реальной характеристики геофона, в которую включен также шумовой фон. В результате этой операции произойдет подавление соответствующей характеристики (полезного сигнала) геофона, в результате практически в чистом виде выделяется характеристика помех. Параллельно осуществляют временное запоминание первоначального реального зашумленного сигнала геофона. После выделения характеристики помех (шумового сигнала) считывают первоначальный реальный зашумленный сигнал геофона из устройства запоминания и вычитают из него соответствующую характеристику помех, чтобы получить незашумленный "чистый" сигнал геофона. К чистому сигналу геофона итерационно применяют коэффициент пересчета, после чего сигнал геофона суммируют с сигналом гидрофона, и при каждой очередной итерации выполняют автокорреляцию получаемой суммы. Величину коэффициента пересчета соответственно определяют, исходя из исследования основного и первого бокового лепестков автокорреляционной функции. Итерации проводятся до тех пор, пока не будет максимизирована соответствующая мера приемлемости получаемого результата.
На фиг. 4 показан в общем виде контур адаптивного выделения помех, который работает по каждой выборке данных как функции индекса j, времени - выборки. Фильтр 47 преобразования сигналов от гидрофона к геофону включает блоки 46 и 48, в состав которых по каждой новой выборке входит адаптивный фильтр, аналогичный широко известному фильтру Винера. Фильтр обеспечивает оценку скорости по амплитуде, фазе, частоте и затуханию. Сигнал гидрофона вводится на вход блока 46 фильтра 47, где осуществляется его свертка с оценкой глубины воды WDj и динамической амплитудной характеристикой фильтра Aj с образованием оценки сигнала геофона SUMj. Вслед за тем SUMj комбинируется путем вычитания (субтрактивно) с входным сигналом геофона, имеющим высокий уровень помех, в точке суммирования 52. На выходе получаются почти в чистом виде помехи геофона εj, которые вводятся в блок 50 фильтра 47. При обработке используется коэффициент усиления μ цепи обратной связи. Величина коэффициента усиления цепи обратной связи находится в пределах узкого диапазона значений от нуля до 0,01. При коэффициенте усиления, равном нулю, фильтр прекращает функционировать. Если коэффициент усиления чересчур велик, фильтр становится неустойчивым. Поэтому необходимо настраивать коэффициент усиления цепи обратной связи, чтобы обеспечить нужную степень адаптивности. Производится приращение кажущейся глубины водной толщи для формирования нового фильтра преобразования сигналов от гидрофона к геофону. С использованием нового фильтра процесса свертки процесс свертки в блоке 48 повторяется. Если разность Δ WD результатов первоначального и обновленного фильтров улучшает оценку сигнала геофона, то опытное значение приращения глубины водной толщи принимается постоянным в подобных расчетах.
На фиг. 5 G обозначает зашумленный входной сигнал геофона (характеристика скорости), H представляет соответствующий входной сигнал гидрофона (характеристика давления). Почти чистый шумовой сигнал εj, не загрязненный полезным сигналом, вводится с противоположным знаком в зашумленный входной сигнал геофона G ("Сигнал + Помеха") в блоке 54 для получения чистого сигнала геофона, который на фиг. 5 обозначен как "Очищенный сигнал G".
В блоке 56 к сигналу геофона применяется коэффициент пересчета, и результат прибавляется в блоке суммирования 58 к сигналу гидрофона (порядок применения может быть обратным) для получения окончательной величины на выходе. Коэффициент пересчета определяется в блоке 60.
При заданных выборках из соответствующих наборов данных гидрофона и геофона, отсортированных на базе равноудаленных приемников и выровненных по амплитуде, выбирают окно зарегистрированных сигналов, которые характеризуются довольно чистым сигналом высокого уровня. Определяют значение коэффициента пересчета для уточненного сигнала, геофона, полученного на основе выбранных записей сигналов и с использованием этого коэффициента пересчета суммируют полученные величины с соответствующими зарегистрированными сигналами гидрофона. Осуществляют автокорреляционную обработку результата суммирования. Рассчитывают коэффициент сходимости и выполняют итерационную процедуру обработки, с каждой итерацией задавая приращения коэффициенту пересчета, пока не будет пройден весь необходимый диапазон коэффициентов пересчета. Для получения коэффициента пересчета для обработки всей выборки выбирают то значение коэффициента пересчета, которому соответствует максимальная величина коэффициента сходимости, т. е. значение, которое максимально приближено к единице.
Расчет коэффициента сходимости наиболее удобно выполняется на основе варимаксной функции
где
xi; i = 1, 2, 3,..., n,
есть последовательность дискретных выборок данных в пределах выбранного временного окна по функции автокорреляции, включающей область, в которой следует ожидать появления центрального и первых боковых лепестков функции автокорреляции. Если, как показано на фиг. 2 и 3, x представляет собой дискретный всплеск, то VM = 1. Значение VM уменьшается по мере того, как отклик волнового поля все менее напоминает дискретный всплеск.
где
xi; i = 1, 2, 3,..., n,
есть последовательность дискретных выборок данных в пределах выбранного временного окна по функции автокорреляции, включающей область, в которой следует ожидать появления центрального и первых боковых лепестков функции автокорреляции. Если, как показано на фиг. 2 и 3, x представляет собой дискретный всплеск, то VM = 1. Значение VM уменьшается по мере того, как отклик волнового поля все менее напоминает дискретный всплеск.
Блок-схема предпочтительного способа обработки приведена на фиг. 6. Каждый из наборов данных гидрофона и геофона сортируется по общим принимаемым выборкам. Эти данные затем преобразуются в соответствии с инверсным значением коэффициента усиления, выравниваются по амплитуде и корректируются на геометрическое расхождение обычными методами. Частично обработанные данные проходят дальнейшую обработку, предусматривающую исключение шумовых выбросов, фильтрацию (f - k) шумов и обратное преобразование свертки. Все вышеприведенные операции выполняются в соответствии со стандартными процедурами. В блоке 54 процедура адаптивного подавления помех применяется к наборам данных по способу, объяснение которого приводилось со ссылками на фиг. 4 и 5. В блоке 61 осуществляется пересчет данных относительно друг друга способом, объяснение которого приводилось со ссылками на фиг. 5, с использованием наилучшей оценки коэффициента пересчета, который определяется по максимальному значению коэффициента корреляции VM.
Следует иметь в виду, что, хотя способ был описан в отношении сейсмоприемной косы, уложенной на морское дно, его можно применять при морских сейсмических съемках любых других типов, а также при проведении работ в так называемых переходных зонах, где реверберация может иметь место в массе насыщенного влагой берегового песка. Неоднократно упоминалась операция пересчета сигналов геофона в сигналы гидрофона. При этом следует иметь в виду, что в принципе можно выполнить и обратный процесс, т.е. пересчет сигналов гидрофона и их согласование с сигналами геофона. Характеристика конкретных аспектов способа приводится здесь в целях иллюстрации на примерах, без каких-либо ограничений сущности и объема изобретения.
Claims (10)
1. Способ обработки сейсмических данных для подавления многократных переотражений, содержащий следующие приемы: создают отраженное сейсмическое волновое поле в водной среде, причем указанное сейсмическое волновое поле содержит характерный признак давления и характерный признак скорости, который включает в себя характерный признак фонового шума, осуществляют обнаружение характерного признака давления указанного отраженного сейсмического волнового поля в выбранном местоположении и одновременно осуществляют обнаружение характерного признака скорости отраженного сейсмического волнового поля в местоположении, которое является смежным с указанным выбранным положением, и осуществляют фильтрацию характерного признака давления, отличающийся тем, что выделяют указанный характерный признак шума путем субтрактивного комбинирования указанного отфильтрованного характерного признака давления с указанным характерным признаком скорости, получают уточненный характерный признак скорости сложением выделенного характерного признака шума с указанным характерным признаком скорости с обратным знаком, применяют соответствующий коэффициент пересчета для указанного уточненного характерного признака скорости и суммируют указанный характерный признак давления с указанным уточненным и пересчитанным характерным признаком скорости и формируют функцию автокорреляции для суммарного характерного признака.
2. Способ по п.1, отличающийся тем, что рассчитывают коэффициент сходимости для функции автокорреляции, осуществляют итерацию суммирования характерного признака давления с уточненным и пересчитанным характерным признаком скорости и формирования функции автокорреляции для суммарного характерного признака, задают приращение коэффициента пересчета при каждой очередной итерации, пока коэффициент сходимости не достигнет значения, максимально приближенного к единице.
3. Способ по п. 1, отличающийся тем, что при операции выделения осуществляют операцию адаптивного подавления характерного признака скорости в характерном признаке фонового шума.
4. Способ по п.1, отличающийся тем, что осуществляют операцию фильтрации указанного характерного признака давления как функцию глубины водной толщи до операции субтрактивного комбинирования.
5. Способ по п.4, отличающийся тем, что при фильтрации вводят соответствующий коэффициент усиления цепи обратной связи для характерного признака шума.
6. Способ по п.5, отличающийся тем, что осуществляют коррекцию коэффициента усиления цепи обратной связи для обеспечения соответствующей степени адаптивности.
7. Способ по п.2, отличающийся тем, что коэффициент сходимости определяют с помощью варимаксной функции.
8. Способ по п.7, отличающийся тем, что варимаксную функцию определяют по формуле
где xi, i = 1,2,3,...,n,
есть последовательность дискретных выборок данных, взятых в пределах избранного временного окна по упомянутой автокорреляционной функции, причем окно включает в себя главный и первый боковые лепестки автокорреляционной функции.
где xi, i = 1,2,3,...,n,
есть последовательность дискретных выборок данных, взятых в пределах избранного временного окна по упомянутой автокорреляционной функции, причем окно включает в себя главный и первый боковые лепестки автокорреляционной функции.
9. Способ по п.1, отличающийся тем, что пересчет характерного признака давления осуществляют относительно уточненного признака скорости частиц.
10. Способ по п. 1, отличающийся тем, что преобразование характерного признака давления для оценки характерного признака скорости частиц осуществляют по таким параметрам, как амплитуда, фаза, частота, затухание.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US08/101,949 US5365492A (en) | 1993-08-04 | 1993-08-04 | Method for reverberation suppression |
US08/101,949 | 1993-08-04 |
Publications (2)
Publication Number | Publication Date |
---|---|
RU95111130A RU95111130A (ru) | 1997-06-20 |
RU2116657C1 true RU2116657C1 (ru) | 1998-07-27 |
Family
ID=22287331
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU95111130A RU2116657C1 (ru) | 1993-08-04 | 1994-07-13 | Способ обработки сейсмических данных для подавления многократных переотражений |
Country Status (11)
Country | Link |
---|---|
US (1) | US5365492A (ru) |
EP (1) | EP0666992B1 (ru) |
JP (1) | JP2659461B2 (ru) |
AU (1) | AU671036B2 (ru) |
BR (1) | BR9405543A (ru) |
CA (1) | CA2145281C (ru) |
CO (1) | CO4290391A1 (ru) |
DE (1) | DE69420250D1 (ru) |
NO (1) | NO311779B1 (ru) |
RU (1) | RU2116657C1 (ru) |
WO (1) | WO1995004945A1 (ru) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EA021201B1 (ru) * | 2009-12-07 | 2015-04-30 | Пгс Геофизикал Ас | Способ очистки данных морской сейсмической косы от ложных отражений со стороны источника во всей полосе частот |
Families Citing this family (55)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9321125D0 (en) * | 1993-10-13 | 1993-12-01 | Geco As | Method of processing reflection data |
US5448531A (en) * | 1994-05-05 | 1995-09-05 | Western Atlas International | Method for attenuating coherent noise in marine seismic data |
US5442591A (en) * | 1994-06-21 | 1995-08-15 | Western Atlas International | Method for adaptively suppressing noise transients in summed co-sensor seismic recordings |
US5521881A (en) * | 1994-09-02 | 1996-05-28 | Exxon Production Research Company | Method of processing seismic data having multiple reflection noise |
US5774416A (en) * | 1995-04-07 | 1998-06-30 | Pgs, Tensor, Inc. | Method and device for attenuating water column reverberations using co-located hydrophones and geophones in ocean bottom seismic processing |
US5621699A (en) * | 1995-07-07 | 1997-04-15 | Pgs Ocean Bottom Seismic, Inc. | Apparatus and method of calibrating vertical particle velocity detector and pressure detector in a sea-floor cable with in-situ passive monitoring |
US5572483A (en) * | 1995-07-25 | 1996-11-05 | Western Atlas International, Inc. | Method of reducing noise in seismic signals by adaptive filtering of a noise reference |
US5724306A (en) * | 1995-12-29 | 1998-03-03 | Western Atlas International, Inc. | Method for correcting dual sensor data for imperfect geophone coupling using production seismic data |
US5737433A (en) * | 1996-01-16 | 1998-04-07 | Gardner; William A. | Sound environment control apparatus |
FR2743896B1 (fr) * | 1996-01-23 | 1998-04-10 | Geophysique Cie Gle | Procede de traitement de calibration d'une paire de capteurs hydrophone/geophone et procede de prospection sismique mettant en oeuvre ce traitement |
FR2743897B1 (fr) * | 1996-01-23 | 1998-04-10 | Geophysique Cie Gle | Procede de prospection sismique marine au moyen d'un couple de capteurs hydrophone et geophone |
US5754492A (en) * | 1996-02-12 | 1998-05-19 | Pgs Tensor, Inc. | Method of reverberation removal from seismic data and removal of dual sensor coupling errors |
US5696734A (en) * | 1996-04-30 | 1997-12-09 | Atlantic Richfield Company | Method and system for eliminating ghost reflections from ocean bottom cable seismic survey signals |
US5621700A (en) * | 1996-05-20 | 1997-04-15 | Schlumberger Technology Corporation, Geco-Prakla Div. | Method for attenuation of reverberations using a pressure-velocity bottom cable |
US5793702A (en) * | 1996-05-21 | 1998-08-11 | Western Atlas International, Inc. | Method for measuring the water bottom reflectivity |
GB9612470D0 (en) * | 1996-06-14 | 1996-08-14 | Geco As | Multiple attenuation method |
US5774417A (en) * | 1996-10-25 | 1998-06-30 | Atlantic Richfield Company | Amplitude and phase compensation in dual-sensor ocean bottom cable seismic data processing |
US5818795A (en) * | 1996-10-30 | 1998-10-06 | Pgs Tensor, Inc. | Method of reduction of noise from seismic data traces |
GB9800741D0 (en) * | 1998-01-15 | 1998-03-11 | Geco As | Multiple attenuation of multi-component sea-bottom data |
GB2337591B (en) | 1998-05-20 | 2000-07-12 | Geco As | Adaptive seismic noise and interference attenuation method |
US6151275A (en) * | 1998-09-11 | 2000-11-21 | Pgs Tensor, Inc. | Method of dual wavefield reinforcement |
US6246637B1 (en) | 1998-09-28 | 2001-06-12 | Pgs Tensor, Inc. | Method and system for combining three component seismic data |
US6094620A (en) * | 1999-01-08 | 2000-07-25 | Exxonmobil Upstream Research Company | Method for identifying and removing multiples from seismic reflection data |
US6725164B1 (en) * | 1999-03-17 | 2004-04-20 | Input/Output, Inc. | Hydrophone assembly |
GB9906456D0 (en) | 1999-03-22 | 1999-05-12 | Geco Prakla Uk Ltd | Method and system for reducing effects of sea surface ghost contamination in seismic data |
US20020118602A1 (en) | 2001-02-27 | 2002-08-29 | Sen Mrinal K. | Angle dependent surface multiple attenuation for two-component marine bottom sensor data |
WO2002097472A2 (en) * | 2001-05-25 | 2002-12-05 | Exxonmobil Upstream Research Company | Multiple suppression for ocean bottom seismic data |
US6738715B2 (en) | 2001-09-14 | 2004-05-18 | Exxonmobil Upstream Research Company | Method for attenuating noise in seismic data |
GB2389183B (en) * | 2002-05-28 | 2006-07-26 | Westerngeco Ltd | Processing seismic data |
CA2476512C (en) * | 2002-03-05 | 2012-07-31 | Robert H. Whidden | Method of transmitting electrical power |
US6903277B2 (en) | 2002-03-05 | 2005-06-07 | Robert H Whidden | Conduit for use in the transmission of electrical power |
CN1295520C (zh) * | 2002-03-20 | 2007-01-17 | 输入输出公司 | 用于地震数据采集的自适应滤波装置和方法 |
GB2415258B8 (en) * | 2002-08-30 | 2007-06-07 | Pgs Americas Inc | Method for multicomponent marine geophysical data gathering |
US7239577B2 (en) | 2002-08-30 | 2007-07-03 | Pgs Americas, Inc. | Apparatus and methods for multicomponent marine geophysical data gathering |
GB2395305B (en) * | 2002-11-15 | 2006-03-22 | Westerngeco Seismic Holdings | Processing seismic data |
US6894948B2 (en) * | 2003-01-29 | 2005-05-17 | Pgs Exploration (Uk) Limited | Method for processing dual sensor seismic data to attenuate noise |
US7561493B2 (en) * | 2003-05-30 | 2009-07-14 | Fairfield Industries, Inc. | Method and apparatus for land based seismic data acquisition |
US7310287B2 (en) * | 2003-05-30 | 2007-12-18 | Fairfield Industries Incorporated | Method and apparatus for seismic data acquisition |
GB2405473B (en) * | 2003-08-23 | 2005-10-05 | Westerngeco Ltd | Multiple attenuation method |
GB2410551B (en) * | 2004-01-30 | 2006-06-14 | Westerngeco Ltd | Marine seismic acquisition system |
US8534959B2 (en) | 2005-01-17 | 2013-09-17 | Fairfield Industries Incorporated | Method and apparatus for deployment of ocean bottom seismometers |
US7319636B2 (en) * | 2005-03-14 | 2008-01-15 | Westerngeco, L.L.C. | Calibration of pressure gradient recordings |
US8127706B2 (en) * | 2005-05-02 | 2012-03-06 | Fairfield Industries Incorporated | Deck configuration for ocean bottom seismometer launch platforms |
US7466625B2 (en) * | 2006-06-23 | 2008-12-16 | Westerngeco L.L.C. | Noise estimation in a vector sensing streamer |
US8593907B2 (en) * | 2007-03-08 | 2013-11-26 | Westerngeco L.L.C. | Technique and system to cancel noise from measurements obtained from a multi-component streamer |
US20080270035A1 (en) * | 2007-04-26 | 2008-10-30 | Ahmet Kemal Ozdemir | System and Technique to Remove Perturbation Noise from Seismic Sensor Data |
US7957906B2 (en) * | 2008-05-07 | 2011-06-07 | Pgs Geophysical As | Method for attenuating low frequency noise in a dual-sensor seismic streamer |
US8611191B2 (en) * | 2008-05-22 | 2013-12-17 | Fairfield Industries, Inc. | Land based unit for seismic data acquisition |
US8902699B2 (en) | 2010-03-30 | 2014-12-02 | Pgs Geophysical As | Method for separating up and down propagating pressure and vertical velocity fields from pressure and three-axial motion sensors in towed streamers |
US9182512B2 (en) * | 2011-10-26 | 2015-11-10 | Westerngeco L.L.C. | Processing multi-component seismic data |
US9448315B2 (en) * | 2011-12-27 | 2016-09-20 | Cgg Services Sa | Device and method for denoising ocean bottom data |
CN103487794B (zh) * | 2012-06-13 | 2016-01-06 | 中国科学院声学研究所 | 一种基于小波包变换的水底混响抑制方法 |
WO2014179282A1 (en) * | 2013-04-29 | 2014-11-06 | Westerngeco Llc | Deghosting with adaptive operators |
CN104597500B (zh) * | 2013-10-31 | 2017-05-10 | 中国石油天然气集团公司 | 一种水陆检波器地震数据匹配方法 |
US11249214B2 (en) * | 2018-05-30 | 2022-02-15 | Saudi Arabian Oil Company | Noise suppression of seafloor geophone seismic data |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4253164A (en) * | 1978-10-30 | 1981-02-24 | Western Geophysical Co. Of America | Multi-purpose seismic transducer |
US4910716A (en) * | 1989-01-31 | 1990-03-20 | Amoco Corporation | Suppression of coherent noise in seismic data |
EG19158A (en) * | 1989-08-25 | 1996-02-29 | Halliburton Geophys Service | System for attenuation of water-column reverberation |
US5235554A (en) * | 1991-03-11 | 1993-08-10 | Halliburton Geophysical Services, Inc. | Method for correcting impulse response differences of hydrophones and geophones as well as geophone coupling to the water-bottom in dual-sensor, bottom-cable seismic operations |
US5257241A (en) * | 1991-05-08 | 1993-10-26 | Atlantic Richfield Company | Method and system for acquisition of 3-dimensional marine seismic data |
US5163028A (en) * | 1991-09-27 | 1992-11-10 | Halliburton Geophysical Services, Inc. | Method for correcting impulse response differences of hydrophones and geophones as well as geophone coupling to the water-bottom in dual-sensor, bottom-cable seismic operations |
US5251183A (en) * | 1992-07-08 | 1993-10-05 | Mcconnell Joseph R | Apparatus and method for marine seismic surveying utilizing adaptive signal processing |
-
1993
- 1993-08-04 US US08/101,949 patent/US5365492A/en not_active Expired - Lifetime
-
1994
- 1994-07-13 CA CA002145281A patent/CA2145281C/en not_active Expired - Fee Related
- 1994-07-13 JP JP7506413A patent/JP2659461B2/ja not_active Expired - Lifetime
- 1994-07-13 AU AU73280/94A patent/AU671036B2/en not_active Ceased
- 1994-07-13 BR BR9405543-2A patent/BR9405543A/pt not_active Application Discontinuation
- 1994-07-13 WO PCT/US1994/007739 patent/WO1995004945A1/en active IP Right Grant
- 1994-07-13 EP EP94923411A patent/EP0666992B1/en not_active Expired - Lifetime
- 1994-07-13 RU RU95111130A patent/RU2116657C1/ru active
- 1994-07-13 DE DE69420250T patent/DE69420250D1/de not_active Expired - Lifetime
- 1994-07-29 CO CO94033344A patent/CO4290391A1/es unknown
-
1995
- 1995-03-31 NO NO19951270A patent/NO311779B1/no not_active IP Right Cessation
Non-Patent Citations (1)
Title |
---|
2. US, патент, 5163028, кл . G 01 V 1/38, 1992. 3. * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EA021201B1 (ru) * | 2009-12-07 | 2015-04-30 | Пгс Геофизикал Ас | Способ очистки данных морской сейсмической косы от ложных отражений со стороны источника во всей полосе частот |
US9025413B2 (en) | 2009-12-07 | 2015-05-05 | Pgs Geophysical As | Method for full-bandwidth source deghosting of marine seismic streamer data |
Also Published As
Publication number | Publication date |
---|---|
NO311779B1 (no) | 2002-01-21 |
BR9405543A (pt) | 1999-09-08 |
CO4290391A1 (es) | 1996-04-17 |
CA2145281C (en) | 1998-08-11 |
AU671036B2 (en) | 1996-08-08 |
NO951270L (no) | 1995-05-18 |
DE69420250D1 (de) | 1999-09-30 |
RU95111130A (ru) | 1997-06-20 |
EP0666992B1 (en) | 1999-08-25 |
JP2659461B2 (ja) | 1997-09-30 |
NO951270D0 (no) | 1995-03-31 |
EP0666992A1 (en) | 1995-08-16 |
CA2145281A1 (en) | 1995-02-16 |
US5365492A (en) | 1994-11-15 |
JPH08501395A (ja) | 1996-02-13 |
AU7328094A (en) | 1995-02-28 |
EP0666992A4 (en) | 1998-04-29 |
WO1995004945A1 (en) | 1995-02-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2116657C1 (ru) | Способ обработки сейсмических данных для подавления многократных переотражений | |
US5621700A (en) | Method for attenuation of reverberations using a pressure-velocity bottom cable | |
US6894948B2 (en) | Method for processing dual sensor seismic data to attenuate noise | |
US6654693B2 (en) | Angle dependent surface multiple attenuation for two-component marine bottom sensor data | |
EP0681193A1 (en) | Method of attenuating coherent noise in marine seismic data | |
MX2010012105A (es) | Metodo para atenuacion de ruido de baja frecuencia en un cable marino sismico de sensor-dual. | |
AU750647B2 (en) | Weighted backus filter method of combining dual sensor traces | |
EP0851243B1 (en) | Measuring sea floor reflectivity | |
US5963507A (en) | Method for improving the accuracy of ocean bottom reflectivity estimations using the inverse backus filter | |
US3350683A (en) | Optimum detector type and depth in marine seismic exploration | |
Paffenholz et al. | An improved method for deriving water-bottom reflectivities for processing dual-sensor ocean-bottom cable data | |
EP0689064B1 (en) | Methods of adaptively reducing noise transients in co-sensor seismic signals | |
CN113514889B (zh) | 一种提升海洋深反射地震数据中低频信号能量的处理方法 | |
GB2375606A (en) | Angle dependent surface multiple attenuation for two - component marine bottom sensor data | |
WO2001020363A2 (en) | Amplitude spectra estimation | |
AU765956B2 (en) | Method of deriving water bottom reflectivities in the presence of geologic interference | |
Gelius et al. | Inversion of refracted data: an automatic procedure | |
Holland et al. | High resolution geoacoustic inversion in shallow water: A joint time and frequency domain technique | |
Demoulin et al. | Operational Feasibility of the Geoacoustic Inversion of Equivalent Media: Principle and Application to the INTIMATE96 Data | |
Robinson | REMOTE DETECTION USING DUAL SENSORS |