Способ и система для преобразования трехмерных данных буксируемой сейсмической косы с двумя типами датчиков в данные с подавленным шумом

СПОСОБ И СИСТЕМА ДЛЯ ПРЕОБРАЗОВАНИЯ ТРЕХМЕРНЫХ ДАННЫХ БУКСИРУЕМОЙ СЕЙСМИЧЕСКОЙ КОСЫ С ДВУМЯ ТИПАМИ ДАТЧИКОВ В ДАННЫЕ С ПОДАВЛЕННЫМ ШУМОМ(71)(73) Заявитель и патентовладелец: ПГС ГЕОФИЗИКАЛ АС (NO) Записи давления и записи вертикальной скорости сейсмической массы в составе данных буксируемой сейсмической косы с двумя типами датчиков преобразуются в область продольных волновых чисел. К преобразованным записям вертикальной скорости сейсмической массы при каждом значении продольного волнового числа применяется ряд масштабирующих фильтров, причем все масштабирующие фильтры в ряде фильтров вычисляются для разных диапазонов поперечных относительно косы волновых чисел в блоках продольных сейсмических трасс, в которых все сейсмические события приблизительно линейны. Спектр давления и масштабированный спектр вертикальной скорости сейсмической массы комбинируются для разделения восходящей и падающей компонент волнового поля. Разделенные восходящая и падающая компоненты волнового поля преобразуются обратно в пространственно-временную область. 014279 Область техники, к которой относится изобретение Настоящее изобретение в целом относится к области геофизических изысканий. Более конкретно,изобретение относится к области обработки данных морской сейсмической разведки. Предшествующий уровень техники Геофизические изыскания широко используются в нефтегазодобывающей промышленности для поиска и оценки подземных формаций. Технологии геофизических изысканий обеспечивают возможность получения информации о строении недр земли, используемой для поиска и добычи ценных полезных ископаемых, в частности месторождений углеводородов, таких как нефть и природный газ. Сейсмическая разведка представляет собой хорошо известную технологию геофизических изысканий. При проведении наземной сейсмической разведки сейсмический сигнал, генерируемый на поверхности земли или вблизи нее, распространяется вниз в глубинные отложения. При проведении морской сейсмической разведки сейсмический сигнал может также распространяться вниз через толщу воды, расположенную над глубинными отложениями. Сейсмические источники используются для генерации сейсмического сигнала, который распространяется вглубь земли, после чего, по меньшей мере, частично отражается от подземных сейсмических отражающих границ. Такими сейсмическими отражающими границами обычно являются границы раздела между подземными формациями, обладающими разными упругими свойствами, в частности разными значениями скорости звуковых волн и плотности породы, что приводит к возникновению разности акустических импедансов на границе раздела. Отраженная сейсмическая энергия регистрируется при помощи сейсмических датчиков (также называемых сейсмическими приемниками),установленных на поверхности земли или вблизи нее, в расположенной над дном толще воды или на известной глубине в буровых скважинах и записывается. В качестве сейсмических источников для генерирования сейсмических сигналов в наземной сейсмической разведке можно использовать взрывчатые вещества или вибрационные сейсмические источники. В морской сейсмической разведке обычно применяются подводные сейсмические источники, буксируемые судном и периодически приводимые в действие для генерирования акустического волнового поля. Используются различные типы сейсмических источников, генерирующих волновое поле, в том числе малые взрывчатые заряды, искровые разряды или электрические дуги, морские вибрационные источники и наиболее часто применяемые сейсмические пушки. Сейсмическая пушка, используемая в качестве источника сейсмических колебаний, может представлять собой гидропушку, паровую пушку или чаще всего пневмопушку. Обычно морской сейсмический источник состоит не из одного элемента источника, но представляет собой распределенную в пространстве группу элементов источника. Такую конфигурацию,в частности, применяют в случае пневмопушек, которые в настоящее время представляют собой наиболее часто используемый тип морского сейсмического источника. Подходящие виды сейсмических датчиков включают в себя датчики скорости сейсмической массы,в особенности в наземной разведке, и датчики давления воды (в качестве которых обычно применяются датчики градиента давления воды), в особенности в морской разведке. Иногда вместо датчиков скорости сейсмической массы или в дополнение к ним также используются датчики ускорения сейсмической массы. Датчики скорости сейсмической массы и датчики давления воды общеизвестны в данной области под названием соответственно геофонов и гидрофонов. Сейсмические датчики могут быть установлены по отдельности, но чаще устанавливаются в виде групп датчиков. Кроме того, при проведении морской разведки датчики давления и датчики скорости сейсмической массы могут быть установлены совместно,причем их объединяют в пары датчиков или в пары групп датчиков. При типичной морской сейсмической разведке судно сейсмической разведки, содержащее оборудование для сбора сейсмических данных, такое как оборудование навигационного управления, оборудование управления сейсмическими источниками, оборудование управления сейсмическими датчиками и записывающее оборудование, перемещается по поверхности воды, обычно со скоростью около 5 узлов. Оборудование управления сейсмическими источниками обеспечивает активацию сейсмического источника, буксируемого в воде судном сейсмической разведки, в определенные моменты времени. Сейсмические косы, также называемые сейсморазведочными кабелями, представляют собой вытянутую конструкцию в виде кабеля, которую буксирует в воде судно сейсмической разведки, буксирующее сейсмический источник, или иное судно сейсмической разведки. Как правило, судно сейсмической разведки буксирует несколько сейсмических кос. Сейсмические косы содержат датчики для регистрации отраженных волновых полей, порожденных сейсмическим источником и отраженных от отражающих границ разделов. Обычно сейсмические косы содержат датчики давления, например гидрофоны; однако предлагались и сейсмические косы, содержащие в дополнение к гидрофонам водные датчики скорости сейсмической массы, например геофоны, или датчики ускорения сейсмической массы, например акселерометры. Датчики давления и датчики движения сейсмической массы могут быть установлены на сейсморазведочном кабеле в непосредственной близости друг от друга, объединены в пары датчиков или пары групп датчиков. Конечные сейсмические данные, полученные в результате проведения разведки, обрабатываются с целью получения информации о геологическом строении и свойствах подземных образований в разведываемой области. Обработанные сейсмические данные проходят обработку для отображения и анализа-1 014279 потенциального содержания углеводородов в таких подземных формациях. Цель обработки сейсмических данных заключается в извлечении из сейсмических данных максимального количества информации о подземных образованиях с целью адекватного отражения геологического строения недр. На сбор, обработку и интерпретацию сейсмических данных с целью выявления в недрах земли местоположений, в которых возможно обнаружение залежей нефти, затрачиваются крупные денежные средства. В результате восстановления отражающих поверхностей, задающих представляющие интерес подземные слои, по записанным сейсмическим данным делает возможным получение схемы строения земли в глубину или же по времени. Схемы строения недр земли создаются, чтобы позволить интерпретатору выбрать местоположения с наибольшей вероятностью наличия залежей нефти. Для проверки наличия нефти необходимо бурение скважины. Бурение скважин с целью установления наличия или отсутствия залежей нефти представляет собой чрезвычайно дорогостоящее предприятие, отнимающее много времени. Поэтому существует постоянная потребность в совершенствовании обработки и отображения сейсмических данных с тем, чтобы получать картину строения недр земли, которая повысит возможности интерпретатора по оценке вероятности наличия в определенном месте в недрах земли залежей нефти, при этом интерпретация может производиться как компьютером, так и человеком. Сейсмические данные отражений, полученные с помощью буксируемой сейсмической косы с двумя типами датчиков, состоят из записей поля давления и поля вертикальной скорости сейсмической массы. Центральный этап в последовательности обработки сейсмических данных состоит в разделении данных на записи, содержащие только восходящую и падающую компоненты волновых полей давления. Такое разделение может быть осуществлено после преобразования данных в область частот - волновых чисел(f-kx-ky), вычисления разности и суммы измеренного частотно-волнового спектра записи давления и масштабированного частотно-волнового спектра записи вертикальной скорости сейсмической массы и деления полученных спектров на два. Следует отметить, что использование обращения вышеописанного масштабирующего фильтра позволяет получить, в ином варианте, восходящие и падающие компоненты волновых полей вертикальной скорости сейсмической массы путем вычисления суммы и разности частотно-волнового спектра записи вертикальной скорости сейсмической массы и масштабированного частотно-волнового спектра записи давления и деления полученных спектров на два. Обратное преобразование из области частот - волновых чисел в пространственно-временную область позволяет получить требуемые восходящие и падающие компоненты волновых полей. В описанной процедуре путем масштабирования изменяется только запись вертикальной скорости сейсмической массы (или, в альтернативном варианте, только запись давления). Для неисчезающей энергии она масштабируется в области частот - волновых чисел при помощи реального фильтра, систематически увеличивающегося с увеличением волнового числа для заданной частоты. Однако в морской сейсмической разведке чрезвычайно часто присутствует пространственное наложение частот (аляйсинг) в направлении, перпендикулярном сейсмической косе. В случае поперечного аляйсинга энергия "сворачивается" в меньшее волновое число ky поперечного косе направления. Вследствие этого, если такие эффекты "свертки" для волнового числа не принимать во внимание, то масштабирующий фильтр вычисляется по неправильному, слишком малому волновому числу. Таким образом, энергия наложенных частот в записи вертикальной скорости сейсмической массы (или, в альтернативном варианте, в записи давления) масштабируется с использованием коэффициентов фильтра, которые всякий раз оказываются слишком малыми (или, соответственно, слишком большими). Суперпозиция восходящих и падающих компонент волновых полей в исходных записях приводит к возникновению характерной структуры паразитных провалов на приемнике в соответствующих частотно-волновых спектрах. В каждом случае падения записанной энергии при определенной комбинации частоты и волнового числа, например, в спектре поля давления, соответствующая записанная энергия в спектре вертикальной скорости сейсмической массы максимальна. Однако такое соответствие приводит к некорректному разделению восходящих и падающих компонент волнового поля на провалах в частотно-волновом спектре записи давления (или, в альтернативном варианте, записи вертикальной скорости сейсмической массы). В этих положениях энергия наложенных частот некорректно масштабируется для записи вертикальной скорости сейсмической массы (или, в альтернативном варианте, записи давления), в результате чего неправильно вычисляются конечные разделенные компоненты волнового поля. Таким образом, существует потребность в способе разделения восходящей и падающей компонент волнового поля в трехмерных сейсмических данных буксируемой сейсмической косы с двумя типами датчиков. В качестве аналога предложенного изобретения можно рассмотреть решение согласно патентному документу США 5621699, 15.04.1997, относящееся к устройству и способу калибровки датчиков давления и вертикальной скорости сейсмической массы. Для масштабирования и комбинирования сигналов совмещенных сейсмических приемников в виде гидрофона и геофона известный способ предусматривает получение сигналов помех окружающей среды для этих датчиков и генерирование нормализованного относительного шумового сигнала, представляющего приведенную разность между шумовыми сигналами гидрофона и геофона. На основе этого нормализованного относительного шумового сигнала регули-2 014279 руется относительная амплитуда между данными гидрофона и геофона. В частности, масштабированные сигналы этих датчиков комбинируются для получения линейной комбинации данных. Сущность изобретения Изобретение предлагает способ разделения восходящей и падающей компонент в трехмерных данных буксируемой сейсмической косы с двумя типами датчиков, которые могут содержать энергию наложенных частот в поперечном направлении сейсмической косы. Записи давления и записи вертикальной скорости в составе данных сейсмической косы преобразуются в область продольных волновых чисел. К преобразованным записям вертикальной скорости сейсмической массы для каждого продольного волнового числа применяется ряд масштабирующих фильтров, причем все фильтры из ряда масштабирующих фильтров вычисляются для различных диапазонов поперечных относительно косы волновых чисел и по блокам продольных сейсмических трасс, в которых все сейсмические события приблизительно линейны. Спектр давления и масштабированный спектр вертикальной скорости сейсмической массы комбинируются с целью разделения восходящей и падающей компонент волнового поля. Разделенные восходящая и падающая компоненты волнового поля преобразуются обратно в пространственно-временную область. Перечень фигур чертежей Изобретение и его преимущества станут более ясны из нижеследующего подробного описания и прилагаемых чертежей. На фиг. 1 представлена блок-схема, иллюстрирующая этапы обработки по первому варианту осуществления способа по изобретению для разделения восходящей и падающей компонент волнового поля в трехмерных данных буксируемой сейсмической косы с двумя типами датчиков. На фиг. 2 представлена блок-схема, иллюстрирующая начальные этапы обработки по второму варианту осуществления способа по изобретению для разделения восходящей и падающей компонент волнового поля в трехмерных данных буксируемой сейсмической косы с двумя типами датчиков. На фиг. 3 представлена блок-схема, иллюстрирующая промежуточные этапы обработки по одному из вариантов осуществления изобретения, служащие для обработки окон по фиг. 2. На фиг. 4 представлена блок-схема, иллюстрирующая промежуточные этапы обработки по одному из вариантов осуществления изобретения, служащие для обработки спектральных отношений по фиг. 3. На фиг. 5 представлена блок-схема, иллюстрирующая заключительные этапы обработки по второму варианту осуществления проиллюстрированного на фиг. 1-4 способа по изобретению для разделения восходящей и падающей компонент волнового поля в трехмерных данных буксируемой сейсмической косы с двумя типами датчиков. На фиг. 6 представлена схема, иллюстрирующая третью панель, используемую в способе по изобретению. На фиг. 7 представлена схема, иллюстрирующая маски, созданные из полос третьей панели, представленной на фиг. 6. На фиг. 8 представлен график амплитудных спектров эталонной сейсмической трассы в сравнении с результатами применения способа по изобретению с разными значениями максимального порядка L очистки от наложений. На фиг. 9 представлен график амплитудных спектров эталонной сейсмической трассы в сравнении с другими результатами применения способа по изобретению с разными значениями максимального порядка L очистки от наложений. На фиг. 10 представлен график сигнала, оцененного посредством способа по изобретению, и его отличие от эталонной сейсмической трассы при порядке очистки от наложений L=0. На фиг. 11 представлен график сигнала, оцененного посредством способа по изобретению, и его отличие от эталонной сейсмической трассы при порядке очистки от наложений L=5. Хотя приведенное ниже описание изобретения касается предпочтительных вариантов его осуществления, следует понимать, что изобретение не ограничивается данными вариантами. Напротив, предполагается, что изобретение охватывает все альтернативные, видоизмененные и эквивалентные варианты,которые могут быть заключены в объеме изобретения, определенного пунктами прилагаемой формулы. Сведения, подтверждающие возможность осуществления изобретения Перед корректным масштабированием записи вертикальной скорости сейсмической массы для разделения волнового поля ее спектр необходимо должным образом очистить от наложений (от аляйсинга). Известный способ достижения этой цели заключается в интерполяции сейсмических трасс в поперечном относительно косы направлении с целью уменьшения пространственного разнесения в косе и, соответственно, увеличения поперечного волнового числа Найквиста. Способ по изобретению предусматривает альтернативный подход к обработке энергии наложенных частот в процессе разделения волнового поля без необходимости в непосредственной интерполяции сейсмических трасс. Способ по изобретению использует циклические свойства быстрого преобразования Фурье (БПФ). Однако для разделения волнового поля не требуется вычисление интерполированных сейсмических трасс. Вместо этого к исходному спектру записи вертикальной скорости сейсмической массы для каждого значения продольного (параллельного сейсмической косе) волнового числа kx применяются несколько масштабирующих фильтров. Все такие масштабирующие фильтры вычисляются в разных диапазонах-3 014279 поперечных относительно косы волновых чисел, причем каждый из них применяется только к той части энергии в спектре, которая находится или в случае энергии наложенных частот (aliased energy) должна находиться в данном диапазоне волновых чисел. Разделение волнового поля в соответствии со способом по изобретению осуществляется по блокам, заданным в продольном направлении, так, что они содержат только те сейсмические события, которые приблизительно линейны. Только для упрощения и большей наглядности описания ниже данная процедура описана для единственного блока. На фиг. 1-5 представлены блок-схемы, иллюстрирующие варианты осуществления способа по изобретению для разделения волнового поля. На фиг. 1 и 2-5 представлены соответственно два варианта осуществления способа по изобретению. На фиг. 2 представлены начальные этапы второго варианта осуществления способа по изобретению, на фиг. 3 и 4 представлены дальнейшие промежуточные этапы способа, представленного на фиг. 2, а на фиг. 5 представлены конечные этапы способа по фиг. 2-4. Фиг. 6-11 иллюстрируют некоторые из этапов, описанных в блок-схемах, поясняемых со ссылками на фиг. 1-5. На фиг. 6 и 7 представлены соответственно спектральные отношения, третья панель и образуемые по ней маски. Фиг. 8 и 9 иллюстрируют сравнение амплитудных спектров по результатам способа по изобретению для разных значений максимального порядка очистки от наложений. Фиг. 10 и 11 иллюстрируют сигналы, оцениваемые при помощи способа по изобретению для значений порядка L очистки от наложений L=0 и L=5 соответственно. На фиг. 1 представлена блок-схема, иллюстрирующая этапы обработки, используемые в первом варианте осуществления способа по изобретению для разделения восходящей и падающей компонент волнового поля в трехмерных данных буксируемой сейсмической косы с двумя типами датчиков. В данных сейсмической косы может присутствовать энергия наложенных частот в поперечном направлении, которая будет скорректирована при помощи способа по изобретению. На шаге 11 записи вертикальной скорости сейсмической массы и записи давления в составе данных буксируемой сейсмической косы с двумя типами датчиков преобразуются из пространственновременной области (t-x-y) в область (t-kx-y) продольного волнового числа, где t - время, х и y - соответственно продольная и поперечная относительно косы (перпендикулярная сейсмической косе) координаты,а kx - продольное волновое число. На шаге 12 к преобразованным на шаге 11 записям вертикальной скорости сейсмической массы применяются несколько масштабирующих фильтров при каждом значении продольного волнового числаkx. Таким образом, масштабирующие фильтры применяются к участку с постоянным значением kx. Каждый из ряда масштабирующих фильтров вычисляется для своего отличного диапазона поперечных волновых чисел, по блокам продольных сейсмических трасс, в которых все сейсмические события приблизительно линейны. Каждый из масштабирующих фильтров воздействует только на соответствующем участке энергии в частотно-волновом спектре, который находится в соответствующем диапазоне волновых чисел в случае энергии без наложений (non-aliased energy) или должен находиться в соответствующем диапазоне волновых чисел в случае энергии наложенных частот. На шаге 13 спектр давления и масштабированный спектр вертикальной скорости сейсмической массы, полученный на шаге 12, комбинируются с целью разделения восходящей и падающей компонент волнового поля. На шаге 14 разделенные восходящая и падающая компоненты волнового поля, полученные на шаге 13, преобразуются обратно в пространственно-временную область (t-x-y). На фиг. 2 представлена блок-схема, иллюстрирующая начальные этапы обработки по второму варианту осуществления способа по изобретению для разделения восходящей и падающей компонент волнового поля в трехмерных данных буксируемой сейсмической косы с двумя типами датчиков. Данные сейсмической косы могут содержать энергию наложенных частот в поперечном относительно косы направлении. Фиг. 2 развивает обсуждение, приведенное выше со ссылками на фиг. 1. На шаге 20 записи давления и вертикальной скорости сейсмической массы разбиваются на блоки в пространственно-временной области (t-x-y). Каждый из блоков содержит все сейсмические трассы по поперечной координате (у) и все точки выборки сейсмических трасс по времени (t). Однако блок ограничен по продольной координате (х) так, что все сейсмические события в блоке приблизительно линейны. Таким образом, блоки могут перекрываться по продольной координате. В случае перекрывания блоков необходимо их корректное комбинирование после разделения волнового поля. На шаге 21 выбирается блок, полученный на шаге 20. На шаге 22 записи давления и вертикальной скорости сейсмической массы в блоке, выбранном на шаге 21, преобразуются из области t-x-y в область t-kx-y. На шаге 23 выбираются два интервала с постоянным значением kx, причем один из интервалов выбирается в записи вертикальной скорости сейсмической массы, а второй - в соответствующей записи давления, обе из которых получены на шаге 22. На шаге 24 два интервала с постоянным значением kx, выбранные на шаге 23, разбиваются на несколько окон, содержащих все сейсмические трассы по поперечной координате, но не обязательно все точки выборки сейсмических трасс по времени. Окна могут перекрываться по времени. В случае перекрывания окон необходимо их корректное комбинирование после разделения волнового поля.-4 014279 На шаге 25 выбираются два соответствующих окна, полученных на шаге 24, причем одно из окон выбирается в записи давления, а второе - в записи вертикальной скорости сейсмической массы. На шаге 26 выбирается максимальный порядок очистки от наложений, выражаемый целым числомL. Исходный диапазон поперечных волновых чисел соответствует L=0. На шаге 27 окна, выбранные на шаге 25, передаются на шаг 31 с фиг. 3 для дальнейшей обработки. На шаге 28 определяется, остались ли в интервале с постоянным значением kx, выбранном на шаге 23, окна, выбираемые на шаге 25. Если еще не выбранные окна имеются, то выполнение возвращается к шагу 25. В противном случае выполнение переходит к следующему шагу 29. На шаге 29 определяется, остались ли интервалы с постоянным значением kx, выбираемые на шаге 23. Если еще не выбранные интервалы имеются, то выполнение возвращается к шагу 23. В противном случае выполнение переходит к следующему шагу 30. На шаге 30 определяется, остались ли блоки, выбираемые на шаге 21. Если еще не выбранные блоки имеются, то выполнение возвращается к шагу 21. В противном случае выполнение в соответствии с блок-схемой по фиг. 2 завершается. После этого выполнение переходит к блок-схеме по фиг. 3. На фиг. 3 представлена блок-схема, иллюстрирующая этапы обработки окон по фиг. 2 в соответствии с одним из вариантов осуществления способа по изобретению. Каждое такое окно подвергается следующей процедуре. На шаге 31 с шага 27 по фиг. 2 поступает два окна, одно - из записи давления, а второе - из записи вертикальной скорости сейсмической массы. На шаге 32 вычисляется масштабирующий фильтр для применения в разделении волнового поля при текущем значении kx и в требуемом диапазоне поперечных волновых чисел, который превышает диапазон волновых чисел Найквиста, задаваемый окнами данных, в L раз. На шаге 33 окно записи вертикальной скорости сейсмической массы, полученное на шаге 31, дополняют нулевыми сейсмическими трассами в поперечном направлении, количество которых в L раз превышает количество сейсмических трасс в окне. Здесь L - максимальный порядок очистки от наложений, выбранный на шаге 26 с фиг. 2. На шаге 34 дополненное окно записи вертикальной скорости сейсмической массы, полученное на шаге 33, дополняют нулевыми точками выборки по времени, количество которых в L раз превышает количество точек выборки в сейсмической трассе, при этом образуется первая расширенная запись. ЗдесьL - максимальный порядок очистки от наложений, выбранный на шаге 26 с фиг. 2. На шаге 35 первая расширенная запись, полученная на шаге 34, периодически обнуляется L сейсмическими трассами между каждыми двумя активными сейсмическими трассами, при этом образуется вторая расширенная запись. На шаге 36 первое и второе расширенные окна записи вертикальной скорости сейсмической массы,полученные, соответственно, на шагах 34 и 35, преобразуются из области (t-kx-y) продольных волновых чисел в область (f-kx-ky) частот - волновых чисел. Спектр первой расширенной записи, полученной на шаге 34, в дальнейшем именуется как "первый спектр", а спектр второй расширенной записи, полученной на шаге 35, в дальнейшем именуется соответственно как "второй спектр". Диапазоны частот и волновых чисел теперь содержат в (L+1) раз больше точек выборки, чем имеется в исходных спектрах окон,выбранных на шаге 25 с фиг. 2. В дальнейшем описании число точек выборки в области частот - волновых чисел исходных спектров для частотной координаты обозначено как nf, а для координаты поперечного волнового числа - как nk соответственно. На шаге 37 частотные диапазоны как первого, так и второго спектров, полученных на шаге 36,сжимают, рассматривая для всех волновых чисел наиболее внутренние nf штук выборок. На шаге 38 ко второму спектру, сжатому на шаге 37, добавляется небольшое количество белого шума. Добавление белого шума служит для предотвращения проблем, связанных с тем, что второй спектр окажется малым при его использовании в качестве делителя на следующем шаге 39. На шаге 39 вычисляется отношение первого спектра, полученного на шаге 37, и второго спектра,полученного на шаге 38. Полученное отношение первого и второго спектра в дальнейшем именуется как"третья панель". На шаге 40 масштабирующий фильтр, полученный на шаге 32, и третью панель, полученную на шаге 39, передают на шаг 41 с фиг. 4. На этом выполнение блок-схемы на фиг. 3 завершается, после чего выполнение переходит к блок-схеме по фиг. 4. На фиг. 4 представлена блок-схема, иллюстрирующая промежуточные этапы обработки в соответствии с вариантом осуществления способа по изобретению для обработки спектральных отношений, полученных в соответствии с фиг. 3. На шаге 41 с шага 40 по фиг. 3 принимаются масштабирующий фильтр и третья панель. На шаге 42 третья панель, полученная на шаге 41, разбивается на полосы, количество которых равно удвоенному числу (L+1), причем каждая полоса содержит точек выборки в области волновых чисел. Полосы нумеруются от -L до L, т.е. от отрицательных индексов L к положительным индексам L. На шаге 43 каждые две соответствующие полосы комбинируются в маску перемещением полосы с-5 014279 положительным индексом L в исходный диапазон положительных волновых чисел, а полосы с отрицательным индексом L - в исходный диапазон отрицательных волновых чисел. Если число L нечетное, положение обеих полос меняется на противоположное. На шагах 42 и 43 создаются (L+1) масок, каждая из которых содержит nf x nk точек выборки. На шаге 44 масштабирующий фильтр, полученный на шаге 41, разбивается на полосы, количество которых равно удвоенному числу (L+1), причем каждая полоса содержит точек выборки в области волновых чисел. Полосы имеют номера от -L до L. На шаге 45 каждые две соответствующие полосы комбинируются в фильтр перемещением полосы с положительным индексом L в исходный диапазон положительных волновых чисел, а полосы с отрицательным индексом L - в исходный диапазон отрицательных волновых чисел. Если число L нечетное, положение полос меняется на противоположное. На шагах 44 и 45 создается (L+1) фильтров, соответствующих (L+1) маске, созданной на шагах 42 и 43, при этом фильтры воздействуют на различные диапазоны волновых чисел масок. На шаге 46 все векторы, образованные точками выборки всех масок в выборке с одинаковой частотой и волновым числом, нормализуются. Нормализация необходима, чтобы исключить введение лишней энергии в разделенные волновые поля. На шаге 47 соответствующие фильтры и маски перемножаются, получая ряд масштабирующих фильтров. На шаге 48 ряд масштабирующих фильтров, полученных на шаге 47, передается на шаг 51 с фиг. 5 для дальнейшей обработки. На фиг. 5 представлена блок-схема, иллюстрирующая завершающие этапы обработки по второму варианту осуществления способа по изобретению для разделения восходящей и падающей компонент волнового поля в трехмерных данных буксируемой сейсмической косы с двумя типами датчиков, проиллюстрированного на фиг. 1-4. На шаге 51 с шага 48 по фиг. 4 получают ряд масштабирующих фильтров. На шаге 52 ряд масштабирующих фильтров, полученных на шаге 51, применяется к исходному спектру окна записи вертикальной скорости сейсмической массы, полученному на шаге 25 с фиг. 2. На шаге 53 все полученные фильтрованные спектры окна записи вертикальной скорости сейсмической массы с шага 52 суммируются. Данная сумма представляет собой масштабированный спектр записи вертикальной скорости сейсмической массы в окне. На шаге 54 комбинируются спектр записи давления в окне, полученном на шаге 25 с фиг. 2, и масштабированный спектр записи вертикальной скорости сейсмической массы в окне, полученный на шаге 53. В результате получаются спектры восходящего и падающего волновых полей в окне. На шаге 55 все спектры восходящего и падающего волновых полей, полученные на шаге 54, преобразуются из области (f-kx-ky) частот - волновых чисел обратно в область (t-kx-y) продольных волновых чисел. На шаге 56 спектры, полученные в результате обратного преобразования на шаге 55, комбинируются для всех окон, полученных на шаге 24 по фиг. 2 в интервале с постоянным значением kx. В результате получаются восходящее и падающее волновые поля в интервале. На шаге 57 все интервалы с шага 56, содержащие восходящие и падающие волновые поля, вычисленные для интервалов с шага 23 по фиг. 2, комбинируют в блоке. В результате получается запись восходящего и падающего волновых полей в блоке в области t-kx-y. На шаге 58 восходящая и падающая компоненты волнового поля, полученные на шаге 57, преобразуются из области (t-kx-y) продольных волновых чисел обратно в пространственно-временную область(t-x-y). В результате получаются восходящая и падающая компоненты волнового поля в блоке в пространственно-временной области. На шаге 59 комбинируются восходящая и падающая компоненты волнового поля в блоке, полученном на шаге 58, в результате чего получаются конечные восходящая и падающая компоненты волнового поля в пространственно-временной области. Для дальнейшей стабилизации способа по изобретению может быть применено вычисление ряда масок также для каждого окна записи давления и их комбинирование с рядом масок для записи вертикальной скорости сейсмической массы перед шагом нормализации. Это позволяет усилить общие для обеих записей сигналы и уменьшить влияние случайного шума. Задача данного способа заключается в предсказании отношения энергии наложенных частот и энергии без наложений по их отношению при более низких частотах. Желательно применение способа по изобретению во временных окнах, так как данный способ предполагает наличие ограниченного числа локально линейных событий. Вышеприведенное описание касается варианта осуществления способа по изобретению, в котором запись вертикальной скорости сейсмической массы масштабируется, а запись давления не масштабируется. Однако способ по изобретению предусматривает и альтернативный вариант, в котором масштабируется запись давления, а запись вертикальной скорости сейсмической массы не масштабируется. Вышеприведенное описание может быть очевидным образом изменено с включением такого альтернативного-6 014279 варианта. Для этого достаточно поменять местами упоминания записей вертикальной скорости сейсмической массы и записей давления. Настоящее изобретение описано выше в качестве способа исключительно в целях иллюстрации,однако оно также может быть осуществлено в виде системы. Система по изобретению предпочтительно реализуется при помощи компьютеров, в частности цифровых компьютеров, в сочетании с другим традиционным оборудованием для обработки данных. Такое оборудование для обработки данных, хорошо известное в данной области, может включать в себя любые подходящие сочетания сетевого или компьютерного оборудования для обработки данных, в том числе, но не ограничиваясь указанным, аппаратное обеспечение (процессоры, системы временного и постоянного хранения информации и любое другое подходящее компьютерное оборудование), программное обеспечение (операционные системы, прикладные программы, библиотеки математических программ и любое другое подходящее программное обеспечение), средства соединения (электрического, оптического, беспроводного и т.д.) и периферийное оборудование (устройства ввода и вывода, например, клавиатуры, устройства указания и сканнеры; устройства отображения, например, мониторы и принтеры; средства хранения данных, например, диски и жесткие диски, а также любое другое пригодное оборудование). Кроме того, следует отметить, что разделение восходящих и падающих компонент волнового поля,производимое в вышеприведенном описании на шаге 54 по фиг. 5, может быть произведено в любой момент после применения ряда масштабирующих фильтров к спектру вертикальной скорости сейсмической массы в окне, производимого в вышеприведенном описании на шаге 52 по фиг. 5. Любые связанные с этим изменения в порядке шагов о не выходят за рамки способа по изобретению. Процедура создания масок и результат предлагаемого способа проиллюстрированы на нижеследующих фигурах чертежей. Используемые данные были промоделированы посредством прямого моделирования с использованием двадцати трех сейсмических кос на расстоянии в 25 м. Модель состояла из одной отражающей поверхности с наклоном в 60 в поперечном относительно сейсмической косы направлении и отсутствием наклона в продольном направлении косы. На нижеследующих фигурах чертежей приведены результаты для интервала с kx=0. Разбиение на окна во временной области не применялось, так как данные содержат всего одно событие отражения и соответствующий паразитный сигнал приемника. На фиг. 6 приведена схема, иллюстрирующая пример третьей панели, используемой в способе по изобретению. На фиг. 6 изображена третья панель после сокращения до nf точек выборки для каждого волнового числа. Полосы 61, используемые для создания масок, разделены вертикальными линиями 62. Полосы пронумерованы и помечены номерами от L=-5 до L=5. На фиг. 7 приведена схема, иллюстрирующая маски, созданные из полос примера третьей панели на фиг. 6. На фиг. 7 полосы с фиг. 6 собраны в маски 71 вплоть до значения L=4. Маски представлены после нормализации, производимой на шаге 46 по фиг. 4. Каждая из масок содержит то же число точек выборки, что и исходный спектр записи вертикальной скорости сейсмической массы. На фиг. 8-11 приведено сравнение между смоделированным полем давления, свободным от паразитных сигналов (эталон) и оценкой поля давления для сейсмической трассы номер 11. На фиг. 8 приведен график спектров 81 амплитуды эталонной сейсмической трассы в сравнении с результатами 82 способа по изобретению для разных значений максимального порядка L очистки от наложений. На фиг. 9 приведен график спектров 91 амплитуды эталонной сейсмической трассы в сравнении с другими результатами 92 способа по изобретению для разных значений максимального порядка L очистки от наложений. Чем выше выбранный порядок очистки от наложений, тем выше частота, вплоть до которой спектр оцениваемой сейсмической трассы совпадает со спектром эталонной сейсмической трассы. На фиг. 10 представлен график сигнала 101, оцененного посредством способа по изобретению, и его разность 102 с эталонной сейсмической трассой при порядке очистки от наложений L=0. На фиг. 11 представлен график сигнала 111, оцененного посредством способа по изобретению, и его разность 112 с эталонной сейсмической трассой при порядке очистки от наложений L=5. Вышеизложенное описание следует рассматривать просто как подробное описание отдельных вариантов осуществления настоящего изобретения. В соответствии с раскрытием в изобретение могут быть внесены различные изменения, модификации и альтернативные варианты, не выходящие за рамки настоящего изобретения. Таким образом, вышеприведенное описание не должно рассматриваться как ограничивающее объем охраны изобретения. Объем охраны изобретения определяется только прилагаемой формулой изобретения с учетом эквивалентов. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ преобразования трехмерных данных буксируемой сейсмической косы с двумя типами датчиков в данные с подавленным шумом, в котором преобразуют трехмерные данные буксируемой сейсмической косы с двумя типами датчиков, представляющие подземные образования и полученные группой установленных вблизи поверхности земли сейсмических датчиков, в масштабированные сейсмические данные с подавленным шумом от наложения частот; преобразуют записи давления и записи вертикальной скорости сейсмической массы в составе данных сейсмической косы в область продольных волновых чисел; к преобразованным записям вертикальной скорости сейсмической массы для каждого продольного волнового числа применяют ряд масштабирующих фильтров, причем все фильтры из ряда масштабирующих фильтров вычисляют для различных диапазонов поперечных относительно косы волновых чисел по блокам продольных сейсмических трасс, в которых все сейсмические события приблизительно линейны; комбинируют преобразованные записи давления и масштабированные записи вертикальной скорости сейсмической массы с целью разделения восходящей и падающей компонент волнового поля; и осуществляют обратное преобразование разделенных восходящей и падающей компонент волнового поля в пространственно-временную область, причем обратно преобразованная восходящая компонента волнового поля представляет данные с подавленным шумом. 2. Способ по п.1, отличающийся тем, что при преобразовании записей давления и записей вертикальной скорости сейсмической массы в составе сейсмических данных в область продольных волновых чисел выполняют следующие шаги: получают записи давления и записи вертикальной скорости сейсмической массы, соответственно,из регистраторов давления и регистраторов вертикальной скорости сейсмической массы, установленных на буксируемой сейсмической косе; в данных сейсмической косы выбирают несколько блоков, причем каждый из блоков содержит все сейсмические трассы по поперечной относительно косы координате, все точки выборки сейсмических трасс по времени и минимально достаточное количество сейсмических трасс по продольной координате с тем, чтобы все сейсмические события в каждом блоке были приблизительно линейны; преобразуют записи давления и записи вертикальной скорости сейсмической массы каждого блока в область продольных волновых чисел; в каждом из преобразованных блоков выбирают несколько интервалов с постоянным значением продольного волнового числа; каждый из интервалов разбивают на несколько окон, причем каждое окно содержит все сейсмические трассы по поперечной относительно косы координате; и для каждого из окон выполняют следующие шаги: выбирают для окна максимальный порядок L очистки от наложений; вычисляют масштабирующий фильтр для продольного волнового числа для выбранного диапазона поперечных волновых чисел; дополняют нулевыми сейсмическими трассами и нулевыми точками выборки, формируя первую расширенную запись; первую расширенную запись периодически обнуляют нулевыми точками выборки, формируя вторую расширенную запись; преобразуют первую и вторую расширенную запись в область частот - волновых чисел, формируя соответственно первый и второй спектры; сжимают первый и второй спектры и вычисляют отношение первого спектра к второму спектру, формируя третью панель. 3. Способ по п.2, отличающийся тем, что выбранный диапазон поперечных волновых чисел больше диапазона волновых чисел Найквиста для окна в L раз. 4. Способ по п.2, отличающийся тем, что при дополнении окна нулевыми сейсмическими трассами и нулевыми точками выборки выполняют следующие шаги: дополняют окно нулевыми сейсмическими трассами по поперечной относительно косы координате,число которых в L раз больше числа сейсмических трасс в окне; и дополняют окно нулевыми точками выборки по времени, число которых в L раз больше числа точек выборки в окне. 5. Способ по п.2, отличающийся тем, что при периодическом обнулении первой расширенной записи нулевыми точками выборки периодически обнуляют L сейсмическими трассами первую расширенную запись между каждыми двумя активными сейсмическими трассами. 6. Способ по п.2, отличающийся тем, что при сжатии первого и второго спектра сжимают частотные диапазоны первого и второго спектров до наиболее внутренних точек выборки с частотным диапа-8 014279 зоном, равным числу частот в преобразованных записях вертикальной скорости сейсмической массы. 7. Способ по п.2, отличающийся тем, что при вычислении отношения первого спектра к второму спектру к второму спектру добавляют небольшое количество белого шума перед делением. 8. Способ по п.2, отличающийся тем, что при применении ряда масштабирующих фильтров выполняют следующие шаги: разбивают третью панель на 2(L+1) полос с номерами от -L до L; каждые две полосы с соответствующими номерами комбинируют в маску; разбивают масштабирующий фильтр на 2(L+1) полос с номерами от -L до L; каждые две полосы с соответствующими номерами комбинируют в фильтр; нормализуют маски; перемножают соответствующие фильтры и маски, формируя ряд масштабирующих фильтров; и ряд масштабирующих фильтров применяют к преобразованным записям вертикальной скорости сейсмической массы. 9. Способ по п.8, отличающийся тем, что каждая из полос содержит в два раза меньше точек выборки по поперечной относительно косы координате по сравнению с преобразованными записями вертикальной скорости сейсмической массы. 10. Способ по п.8, отличающийся тем, что при объединении каждых двух полос с соответствующими номерами в фильтр выполняют следующие шаги: полосу с положительным номером L перемещают в диапазон положительных волновых чисел; полосу с отрицательным номером L перемещают в диапазон отрицательных волновых чисел и в случае нечетного номера L полосы меняют местами. 11. Способ по п.8, отличающийся тем, что при нормализации масок нормализуют все векторы точек выборки во всех масках в выборке с одинаковыми частотой и волновым числом. 12. Способ по п.8, отличающийся тем, что при комбинировании спектра давления и масштабированного спектра вертикальной скорости сейсмической массы выполняют следующие шаги: суммируют все фильтрованные преобразованные записи вертикальной скорости сейсмической массы, формируя масштабированный спектр вертикальной скорости сейсмической массы в окне; комбинируют спектр давления и масштабированный спектр вертикальной скорости сейсмической массы в окне, формируя спектры восходящей и падающей компонент волнового поля в окне; производят обратное преобразование спектров восходящей и падающей компонент волнового поля из области частот - волновых чисел в область продольных волновых чисел; комбинируют спектры восходящей и падающей компонент волнового поля, прошедшие обратное преобразование, во всех окнах в интервал с постоянным значением kx, формируя восходящую и падающую компоненты волнового поля в интервале; комбинируют восходящую и падающую компоненты волнового поля во всех интервалах в блок,формируя восходящую и падающую компоненты волнового поля в блоке; производят обратное преобразование восходящей и падающей компонент волнового поля в блоке в пространственно-временную область и комбинируют восходящую и падающую компоненты волнового поля во всех блоках, формируя разделенные восходящую и падающую компоненты волнового поля. 13. Система для преобразования трехмерных данных буксируемой сейсмической косы с двумя типами датчиков в данные с подавленным шумом, содержащая средства преобразования трехмерных данных буксируемой сейсмической косы с двумя типами датчиков, представляющих подземные образования, в масштабированные сейсмические данные с подавленным шумом от наложения частот; регистраторы давления и регистраторы вертикальной скорости сейсмической массы, установленные на сейсмической косе и производящие запись давления и запись вертикальной скорости сейсмической массы соответственно; средства для преобразования записи давления и записи вертикальной скорости сейсмической массы в составе данных сейсмической косы в область продольных волновых чисел; средства для применения к преобразованным записям вертикальной скорости сейсмической массы для каждого продольного волнового числа ряда масштабирующих фильтров, причем все фильтры из ряда масштабирующих фильтров вычислены для различных диапазонов поперечных относительно косы волновых чисел по блокам продольных сейсмических трасс, в которых все сейсмические события приблизительно линейны; средства для комбинирования преобразованных записей давления и масштабированных записей вертикальной скорости сейсмической массы с целью разделения восходящей и падающей компонент волнового поля и средства для обратного преобразования разделенных восходящей и падающей компонент волнового поля в пространственно-временную область, причем обратно преобразованная восходящая компонента волнового поля представляет данные с подавленным шумом. 14. Система по п.13, отличающаяся тем, что средства для преобразования записи давления и записи-9 014279 вертикальной скорости сейсмической массы в составе данных сейсмической косы в область продольных волновых чисел содержат средства для выбора в данных сейсмической косы нескольких блоков, причем каждый блок содержит все сейсмические трассы по поперечной относительно косы координате, все точки выборки сейсмических трасс по времени и минимально достаточное количество сейсмических трасс по продольной координате с тем, чтобы все сейсмические события в каждом блоке были приблизительно линейны; средства для преобразования записей давления и вертикальной скорости сейсмической массы каждого блока в область продольных волновых чисел; средства для выбора в каждом из преобразованных блоков нескольких интервалов с постоянным значением продольного волнового числа; средства для разбиения каждого из интервалов на несколько окон, причем каждое окно содержит все сейсмические трассы по поперечной относительно косы координате; и средства для проведения в каждом из окон следующих операций: выбор для окна максимального порядка L очистки от наложений; вычисление масштабирующего фильтра для продольного волнового числа для выбранного диапазона поперечных волновых чисел; дополнение окна нулевыми сейсмическими трассами и нулевыми точками выборки с формированием первой расширенной записи; периодическое обнуление первой расширенной записи нулевыми точками выборки с формированием второй расширенной записи; преобразование первой и второй расширенной записи в область частот - волновых чисел с формированием соответственно первого и второго спектров; сжатие первого и второго спектров и вычисление отношения первого спектра к второму спектру с формированием третьей панели. 15. Система по п.14, отличающаяся тем, что выбранный диапазон поперечных волновых чисел больше диапазона волновых чисел Найквиста для окна в L раз. 16. Система по п.14, отличающаяся тем, что средства для дополнения окна нулевыми сейсмическими трассами и нулевыми точками выборки содержат средства для дополнения окна нулевыми сейсмическими трассами по поперечной относительно косы координате, число которых в L раз больше числа сейсмических трасс в окне; и средства для дополнения окна нулевыми точками выборки по времени, число которых в L раз больше числа точек выборки в окне. 17. Система по п.14, отличающаяся тем, что средства для периодического обнуления первой расширенной записи нулевыми точками выборки содержат средства для периодического обнуления первой расширенной записи L сейсмическими трассами между каждыми двумя активными сейсмическими трассами. 18. Система по п.14, отличающаяся тем, что средства для сжатия первого и второго спектров содержат средства для сжатия частотных диапазонов первого и второго спектров до наиболее внутренних точек выборки с частотным диапазоном, равным числу частот в преобразованных записях вертикальной скорости сейсмической массы. 19. Система по п.14, отличающаяся тем, что средства для вычисления отношения первого спектра ко второму спектру содержат средства для добавления ко второму спектру небольшого количества белого шума перед делением. 20. Система по п.14, отличающаяся тем, что средства для применения ряда масштабирующих фильтров содержат средства для разбиения третьей панели на 2(L+1) полос с номерами от -L до L; средства для комбинирования каждых двух полос с соответствующими номерами в маску; средства для разбиения масштабирующего фильтра на 2(L+1) полос с номерами от -L до L; средства для комбинирования каждых двух полос с соответствующими номерами в фильтр; средства для нормализации масок; средства для перемножения соответствующих фильтров и масок с формированием ряда масштабирующих фильтров и средства для применения ряда масштабирующих фильтров к преобразованным записям вертикальной скорости сейсмической массы. 21. Система по п.20, отличающаяся тем, что каждая из полос содержит в два раза меньше точек выборки по поперечной относительно косы координате по сравнению с преобразованными записями вертикальной скорости сейсмической массы. 22. Система по п.20, отличающаяся тем, что средства для объединения каждых двух полос с соответствующими номерами в фильтр содержат средства для перемещения полосы с положительным номером L в диапазон положительных волновых чисел; средства для перемещения полосы с отрицательным номером L в диапазон отрицательных волно- 10014279 вых чисел и средства для перемены мест полос в случае нечетного номера L. 23. Система по п.20, отличающаяся тем, что средства для нормализации масок содержат средства для нормализации всех векторов точек выборки во всех масках в выборке с одинаковыми частотой и волновым числом. 24. Система по п.20, отличающаяся тем, что средства для комбинирования спектра давления и масштабированного спектра вертикальной скорости сейсмической массы содержат средства для суммирования всех отфильтрованных преобразованных записей вертикальной скорости сейсмической массы с формированием масштабированного спектра вертикальной скорости сейсмической массы в окне; средства для комбинирования спектра давления и масштабированного спектра вертикальной скорости сейсмической массы в окне с формированием спектров восходящей и падающей компонент волнового поля в окне; средства для обратного преобразования спектров восходящей и падающей компонент волнового поля из области частот - волновых чисел в область продольных волновых чисел; средства для комбинирования спектров восходящей и падающей компонент волнового поля, прошедших обратное преобразование, во всех окнах в интервал с постоянным значением kx с формированием восходящей и падающей компонент волнового поля в интервале; средства для комбинирования восходящей и падающей компонент волнового поля во всех интервалах в блок с формированием восходящей и падающей компонент волнового поля в блоке; средства для обратного преобразования восходящей и падающей компонент волнового поля в блоке в пространственно-временную область и средства для комбинирования восходящей и падающей компонент волнового поля во всех блоках с формированием разделенных восходящей и падающей компонент волнового поля.