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

СПОСОБ ОЧИСТКИ ОТ ЛОЖНЫХ ОТРАЖЕНИЙ ДАННЫХ МОРСКОЙ СЕЙСМИЧЕСКОЙ КОСЫ ВО ВСЕЙ ПОЛОСЕ ЧАСТОТ В предлагаемом способе получают сейсмические данные, записанные с помощью морской сейсмической косы, которые отражают характеристики формаций геологической среды,при этом сейсмические данные собирают посредством размещения ряда сейсмоприемников сверху оцениваемой области формаций геологической среды, причем сейсмоприемники формируют электрический и/или оптический сигнал в ответ на сейсмическую волну. Для преобразования сейсмических данных из пространственно-временной области в спектральную область используется комплексный частотный параметр по Лапласу. Чтобы получить решение нормальной системы уравнений методом наименьших квадратов для системы уравнений очистки от ложных отражений, к преобразованным сейсмическим данным применяется итеративная схема метода сопряженных градиентов с использованием физически обоснованного оператора предобусловливания. Решение подвергается обратному преобразованию в пространственновременную область для получения сейсмических данных с подавленными ложными отражениями,обеспечивающих построение изображений геологической среды. Рийанти Христина Д., Ван Борселен Роальд Дж., Фоккема Якоб Т., Ван ден Берг Петер М. (NL) Хмара М.В. (RU)(71)(73) Заявитель и патентовладелец: ПГС ГЕОФИЗИКАЛ АС (NO) Область техники, к которой относится изобретение Настоящее изобретение в целом относится к области геофизической разведки. Более конкретно,изобретение относится к области очистки от ложных отражений данных морской сейсмической косы. Предшествующий уровень техники В нефтегазовой промышленности геофизическая разведка обычно используется для содействия в поисках и оценке формаций геологической среды. Способы геофизической разведки позволяют получать сведения о структуре геологической среды, полезные для обнаружения и извлечения ценных полезных ископаемых, в особенности месторождений углеводородов, таких как нефть и природный газ. Хорошо известным способом геофизической разведки является сейсморазведка. При наземной сейсморазведке сейсмический сигнал генерируется на поверхности земли или вблизи от нее, а затем распространяется вниз в геологическую среду. При морской сейсморазведке сигнал также может распространяться вниз в водоеме, лежащем выше геологической среды. Источники сейсмических колебаний используются для генерации сейсмического сигнала, который, после распространения в земле, по меньшей мере, частично отражается сейсмоотражающими горизонтами в геологической среде. Такие сейсмоотражающие горизонты, как правило, представляют собой границы между подземными формациями с различными упругими свойствами, конкретно, скоростью упругих волн и плотностью породы, что приводит к разности акустических импедансов на границах. Отраженные сейсмические волны регистрируются сейсмическими датчиками (также именуемыми "сейсмоприемниками") на поверхности земли или вблизи от нее, в вышележащем водоеме, или на известных глубинах в скважинах. Сейсмические датчики генерируют сигналы, как правило, электрические или оптические, от зарегистрированной сейсмической волны, которые записываются для дальнейшей обработки. Результативные сейсмические данные, полученные при проведении сейсморазведки и характерные для геологической среды, обрабатываются для получения информации, относящейся к геологическому строению и свойствам формаций геологической среды в исследуемой зоне. Обработанные сейсмические данные обрабатываются для отображения и анализа потенциального группового состава углеводородов в этих подземных формациях. Целью обработки сейсмических данных является извлечение из сейсмических данных максимального объема информации о подземных формациях для получения адекватного изображения геологической среды. При определении участков геологической среды, где существует вероятность обнаружения залежей нефти, крупные средства расходуются на сбор, обработку и интерпретацию сейсмических данных. Процесс построения поверхностей отражающих горизонтов, определяющих исследуемые подземные пласты по записанным сейсмическим данным, позволяет получить изображение среды по глубине или времени. Изображение строения геологической среды получают с целью предоставить интерпретатору возможность выбора участков, где с наибольшей вероятностью могут находиться залежи нефти. Чтобы удостовериться в наличии нефти, необходимо пробурить скважину. Бурение скважин для определения присутствия нефтяных залежей представляет собой чрезвычайно дорогостоящую и трудоемкую задачу. По этой причине сохраняется потребность в улучшении обработки и отображения сейсмических данных для получения изображения строения геологической среды, которое повысит возможности интерпретатора,независимо от того, выполняется ли интерпретация компьютером или человеком, оценивать вероятность того, что скопление нефти существует на определенном участке геологической среды. Подходящие сейсмические источники для генерации сейсмического сигнала при наземной сейсморазведке могут включать взрывчатые вещества или вибраторы. При морской сейсморазведке, как правило, применяется погружной сейсмический источник, буксируемый судном и периодически активируемый для генерации поля акустических волн. Волновое поле может генерироваться сейсмическими источниками нескольких типов, включая небольшой заряд ВВ (взрывчатого вещества), электрическую искру или дугу, морской вибратор и, как правило, пушку. Пушка сейсмического источника может представлять собой гидропушку, паровую пушку и, чаще всего, пневмопушку. Как правило, морской сейсмоисточник состоит не из одного элемента-источника, а из пространственно-распределенной группы источников. Такое расположение относится в особенности к пневмопушкам, которые в настоящее время являются наиболее распространенным видом морского сейсмоисточника. Соответствующие типы сейсмических датчиков, как правило, включают датчики скорости частиц, в особенности при наземной сейсморазведке, и датчики давления воды, в особенности при морской сейсморазведке. Иногда датчики ускорения частиц или датчики градиента давления используются вместо датчиков скорости частиц или в дополнение к ним. Датчики скорости частиц и датчики давления воды обычно известны специалистам под названием геофонов и гидрофонов соответственно. Сейсмические датчики могут размещаться по отдельности, но чаще размещаются в виде групп датчиков. Кроме того, в ходе морской сейсморазведки датчики давления и датчики скорости частиц могут размещаться совместно, объединенные попарно или в пары групп. При обычной морской сейсморазведке сейсморазведочное судно движется по поверхности воды,как правило, со скоростью около 5 морских узлов и везет оборудование сбора морских сейсмических данных, такое как оборудование навигационного управления, оборудование управления сейсмическими источниками, оборудование управления сейсмическими датчиками и записывающее оборудование. Обо-1 018531 рудование управления сейсмическими источниками активирует сейсмический источник, буксируемый в водоеме сейсморазведочным судном, в выбранные моменты времени. Сейсмические косы, называемые также "сейсмоприемными кабелями", представляют собой удлиненные кабельные конструкции, буксируемые в водоеме сейсморазведочным судном, которое буксирует сейсмический источник, или другим сейсморазведочным судном. Как правило, ряд сейсмических кос буксируются за сейсморазведочным судном. Сейсмические косы содержат датчики для регистрации отраженных волновых полей, возбужденных сейсмическим источником и отраженных от отражающих границ. При традиционной методике сейсмические косы содержат датчики давления, такие как гидрофоны, однако применяются и сейсмические косы, которые в дополнение к гидрофонам содержат датчики скорости частиц воды или датчики ускорения частиц, например акселерометры. Датчики давления и датчики движения частиц обычно размещаются близко друг к другу, объединенные попарно или в пары групп вдоль сейсмического кабеля. После того как отраженная волна достигает кабеля сейсмической косы, она продолжает распространяться к границе вода-воздух на поверхности воды, от которой волна отражается в нижнем направлении и снова регистрируется гидрофонами в сейсмической косе. Поверхность воды является хорошей отражающей границей, при этом коэффициент отражения почти равен единице по величине и отрицателен по знаку для сейсмических сигналов. Волны, отраженные от поверхности, будут, таким образом,сдвинуты по фазе на 180 по отношению к волнам, распространяющимся в верхнем направлении (восходящим). Волны, распространяющиеся в нижнем направлении (нисходящие) и записываемые приемниками, обычно называются "отражением от поверхности" или "волной-спутником". В связи с отражением от поверхности водная поверхность действует как фильтр, который создает вырезы в спектре записанного сигнала, затрудняя запись данных за пределами выбранной полосы частот. Из-за влияния отражения от поверхности некоторые частоты в записанном сигнале усиливаются, тогда как другие ослабляются. В настоящее время при проведении геофизических исследований методом морской сейсморазведки судно буксирует очень длинные сейсмические косы, к которым присоединено большое количество сейсмоприемников, как правило, гидрофонов (хотя могут использоваться и другие типы приемников). Эти приемники регистрируют часть рассеянного поля акустических волн, возбужденных зондирующим сигналом сейсмоисточника. Поле акустических волн, генерируемое сейсмоисточником, рассеивается за счет отражений и преломлений в геологической среде. При традиционном сборе морских сейсмических данных приемники сейсмической косы расположены группами на конкретной глубине ниже поверхности моря. Вследствие такого расположения так называемая однократная отраженная волна, представляющая собой сигнал прямого отклика, проходящий от источника в геологическую среду, а затем к приемникам,затемняется ложным отражением из-за волны, которая проходит от источника в геологическую среду, а затем через поверхность моря к приемникам. Удаление ложного отражения из морских сейсмических данных является первой стадией предварительной обработки данных, увеличивающей разрешающую способность. Эта процедура называется "очисткой от ложных отражений" или "подавлением волнспутников". Таким образом, существует потребность в способе очистки морских сейсмических данных от ложных отражений, сохраняющем стабильность на любых частотах, даже при наличии вырезов в спектре сейсмических данных. Кроме того, способ должен работать без какого-либо предварительного знания о геологической среде. Сущность изобретения Настоящее изобретение предлагает способ преобразования сейсмических данных, записанных с помощью морских сейсмических кос, в сейсмические данные с подавленными ложными отражениями. Получают сейсмические данные, записанные с помощью морских сейсмических кос, при этом сейсмические данные отражают характеристики формаций геологической среды. Сейсмические данные собирают посредством размещения ряда сейсмоприемников сверху оцениваемой области формаций геологической среды, причем сейсмоприемники формируют электрический и/или оптический сигнал в ответ на сейсмическую волну. Для преобразования сейсмических данных из пространственно-временной области в спектральную область используется комплексный частотный параметр по Лапласу. Чтобы получить решение нормальной системы уравнений методом наименьших квадратов для системы уравнений очистки от ложных отражений, к преобразованным сейсмическим данным применяется итеративная схема метода сопряженных градиентов с использованием физически обоснованного предобусловливателя. Решение подвергается обратному преобразованию в пространственно-временную область для получения сейсмических данных с подавленными ложными отражениями, обеспечивающих построение изображений геологической среды. Краткий перечень чертежей Настоящее изобретение и его преимущества легче понять из следующего подробного описания и прилагаемых чертежей, на которых представлены: на фиг. 1 - блок-схема, иллюстрирующая первый вариант осуществления изобретения для очистки от ложных отражений данных морской сейсмической косы; на фиг. 2 - блок-схема, иллюстрирующая начальную часть второго варианта осуществления изобретения для очистки от ложных отражений данных морской сейсмической косы; на фиг. 3 - блок-схема, иллюстрирующая конечную часть второго варианта осуществления изобретения для очистки от ложных отражений данных морской сейсмической косы, начатую на фиг. 2; на фиг. 4 - блок-схема, иллюстрирующая вариант осуществления изобретения для итеративной схемы метода сопряженных градиентов с предобусловливателем, используемой для очистки от ложных отражений данных морской сейсмической косы. Хотя настоящее изобретение будет описано в связи с предпочтительными вариантами осуществления, понятно, что изобретение не ограничивается только ими. Напротив, изобретение предназначено для того, чтобы охватить все альтернативные варианты, модификации и эквиваленты, которые могут быть включены в объем настоящего изобретения согласно прилагаемой формуле изобретения. Сведения, подтверждающие возможность осуществления изобретения Одна из процедур очистки от ложных отражений данных морской сейсмической косы заключается в разложении сейсмических волновых полей для получения восходящей компоненты волнового поля. Один из способов разложения сейсмических волновых полей на восходящие и нисходящие компоненты волнового поля состоит в решении системы уравнений типа Фурье, при котором решение получают путем обращения (инверсии) системы с однозначно определенным ядром Фурье. В одном примере настоящего способа систему уравнений получают в области Лапласа, где прямое преобразование Лапласа из пространственно-временной области (x, t) в область Лапласа (x, s) в общем случае определяется как а соответствующее обратное преобразование Лапласа определяется какs - частотный параметр по Лапласу, который более подробно обсуждается ниже. В дальнейшем обсуждении настоящее изобретение будет иллюстрироваться как действующее на основе волнового поля давления (продольных волн) p(x, t), как правило, регистрируемого датчиками давления, такими как гидрофоны, хотя тип волнового поля и тип датчиков не ограничивают объем настоящего изобретения. Настоящее изобретение может применяться с любым типом соответствующего сейсмического датчика и регистрируемых им волновых полей. В одном из примеров могут применяться волновые поля движения частиц и датчики движения частиц. Следующие уравнения легко изменить для приведения в соответствие с другими волновыми полями. В преобразовании Лапласа, представленном выражениями (1) и (2), лапласовский параметр s, представляющий собой частотный параметр, при традиционной методике является чисто мнимым числом и определяется как гдепредставляет собой круговую частоту, а f - вещественную частоту. Проблема возникает при использовании традиционного частотного параметра по Лапласу s, данного в уравнении (3). Вырезы в спектре сейсмических данных появляются, когда между сейсмическими волнами возникает ослабляющая интерференция при конкретных комбинациях частот и волновых чисел. Когда система уравнений очистки от ложных отражений типа Фурье решена относительно таких частотf, при которых возникают вырезы в спектре, то измеренный сигнал для этих комбинаций частот и волновых чисел равен нулю, поскольку записывается только фоновый шум. При использовании частотного параметра по Лапласу, который является чисто мнимым, как определено уравнением (3), процесс обращения (инверсии) становится неустойчивым. Следствием этой неустойчивости являются неисправимые артефакты, вводимые в результаты разложения сейсмических волновых полей, используемых для очистки от ложных отражений. Чтобы справиться с неустойчивостью, возникающей из-за этих вырезов в спектре, в способе настоящего изобретения используется частотный параметр по Лапласу s, который является комплексной переменной с мнимой и вещественной частями. Таким образом, частотный параметр по Лапласу s задан теперь не выражением (3), а выражением гдеявляется дополнительной вещественной частью комплексного частотного параметра по Лапласу s. В дальнейшем обсуждении термин "комплексный" будет использоваться для указания на то, что частотный параметр по Лапласу имеет как вещественную, так и мнимую части, как в выражении (4). Этот комплексный частотный параметр по Лапласу s будет использоваться для преобразования области Лапласа и спектральной области, а также в уравнениях очистки от ложных отражений, рассматриваемых ниже. В одном из вариантов осуществления вещественная частькомплексного частотного параметра по Лапласу s может быть константой. В других вариантах осуществления вещественная частьможет изменяться в качестве функции таких параметров, как время, пространственное положение или частота. Таким образом, вещественная частьможет быть выражена, но не ограничивается, следующими выражениями: где x1 и x2 представляют собой горизонтальные пространственные координаты, такие как продольные и поперечные направления, а x3 - вертикальную пространственную координату, такую как глубина. Настоящее изобретение представляет собой способ очистки от ложных отражений данных морской сейсмической косы. Изобретение может применяться без каких-либо априорных знаний о геологической среде. Операция очистки от ложных отражений выполняется на основе записанных данных отдельно по каждому выстрелу. Затем для каждой частоты изобретение включает способ решения системы уравнений очистки от ложных отражений. Решение системы уравнений получают с использованием итеративного метода сопряженных градиентов с предобусловливателем. Предобусловливатель основан не на математической (как было бы при традиционном подходе), а на физической стороне рассматриваемой проблемы. В качестве эффективного предобусловливателя применяется инверсия горизонтальной сейсмической косы, поскольку при более общем подходе сейсмическую косу можно рассматривать в качестве возмущения простой горизонтальной сейсмической косы с регулярной структурой. Однако традиционный способ очистки от ложных отражений с использованием мнимого частотного параметра по Лапласу из уравнения (3) может все же становиться неустойчивым в вырезах частотной области. Поэтому в способе настоящего изобретения применяется комплексный частотный параметр по Лапласу из выражения (4), который обеспечивает устойчивость процедуры очистки от ложных отражений. Таким образом, разложение волнового поля может применяться во всем частотном диапазоне сейсмических сигналов даже при наличии вырезов в спектре сейсмических данных. Ниже в ходе дальнейшего обсуждения использование этого комплексного частотного параметра по Лапласу подробно описано со ссылкой на блок-схему, показанную на фиг. 2. На фиг. 1 представлена блок-схема, иллюстрирующая первый вариант осуществления изобретения для очистки от ложных отражений данных морской сейсмической косы. В блоке 11 получают сейсмические данные, записанные с помощью морских сейсмических кос, которые отражают характеристики формаций геологической среды, при этом сейсмические данные собирают с помощью размещения ряда сейсмоприемников, расположенных сверху оцениваемой области формаций геологической среды, причем сейсмоприемники формируют электрический и/или оптический сигнал в ответ на сейсмическую волну. В блоке 12 комплексный частотный параметр по Лапласу используется для преобразования сейсмических данных из блока 11 из пространственно-временной области в спектральную область. В блоке 13 итеративная схема метода сопряженных градиентов с использованием физически обоснованного предобусловливателя применяется к преобразованным сейсмическим данным, чтобы получить решение нормальной системы уравнений методом наименьших квадратов для системы уравнений очистки от ложных отражений. В блоке 14 сейсмические данные из блока 13 подвергаются обратному преобразованию в пространственно-временную область для получения сейсмических данных с подавленными ложными отражениями. На фиг. 2 представлена блок-схема, иллюстрирующая начальную часть второго варианта осуществления изобретения для очистки от ложных отражений данных морской сейсмической косы. Конечная часть обсуждается ниже со ссылкой на блок-схему, показанную на фиг. 3. В блоке 21 получают сейсмические данные, записанные с помощью морской сейсмической косы. Данные приемника по каждой позиции для выстрела (активации сейсмического источника) и по всем записанным временам t принимаются в качестве входных сейсмических данных. Пространственное положение приемников в декартовой системе координат задано уравнением xR=(xR1, xR2, xR3). В настоящем изобретении глубина расположения приемника x3R=x3R(xR1, xR2) представляет собой однозначную функцию горизонтальных координат xR1 и xR2. Таким образом, сейсмические косы не являются вертикальными, но и необязательно расположены горизонтально. Иначе говоря, положение приемников не ограничивает объем настоящего изобретения. В блоке 22 определяется вещественная часть s комплексного частотного параметра по Лапласу s в выражении (4). В различных вариантах осуществления вещественная часть s может быть выражена как в примерах, представленных выше выражением (5), не ограничивая при этом объем настоящего изобретения. Используя различные значения вещественной частикомплексного частотного параметра по Лапласу s, можно варьировать стабилизацию инверсии в способе настоящего изобретения, обеспечивая дополнительную стабилизацию решения приведенной ниже системы уравнений очистки от ложных отражений. В блоке 23 записанные для выстрела сейсмические данные получают по сейсмическим данным, полученным в блоке 21. В блоке 24 сейсмические данные, записанные для выстрела из блока 23, временно преобразуют из пространственно-временной области в пространственно-частотную область Лапласа. В одном из вариантов осуществления настоящего изобретения рассеянное волновое поле psct в местоположениях приемников преобразуется из пространственно-временной области (xR1, xR2, xR3, t) в пространственно-частотную область Лапласа (xR1, xR2, xR3, s) с помощью прямого преобразования Лапласа, заданного выражением (1),тем самым используя в этом преобразовании комплексный частотный параметр по Лапласу s настоящего изобретения, заданный не традиционным выражением (3), а приведенным выше выражением (4),где psct представляет собой рассеянное поле акустических волн в пространственно-частотной области;xR3,q,r - глубина расположения приемника как функция xR1,q и xR2,r, при использовании , определенного в блоке 22, в качестве вещественной части комплексного частотного параметра по Лапласу s. Такое преобразование рассеянного волнового поля psct в местоположениях приемников выполняется для каждого положения выстрела. В блоке 25 преобразованные сейсмические данные из блока 24 получают для частоты f. В блоке 26 преобразованные сейсмические данные для некоторой частоты из блока 25 преобразуются из пространственно-частотной области Лапласа в спектральную область. В одном из вариантов осуществления настоящего изобретения рассеянное волновое поле psct преобразуется из пространственно-частотной области Лапласа (xR1, xR2, xR3, s) в спектральную область (sn, sm, xR3, s) с помощью выраженияPsctn,m представляет собой рассеянное поле акустических волн в спектральной области;xR1 - интервал выборки для приемника в продольном направлении;xR2 - интервал выборки для приемника в поперечном направлении;jsn - спектральный параметр Фурье в продольном направлении;(s) - интервал спектральной выборки в продольном направлении;jsm - спектральный параметр Фурье в поперечном направлении;(s) - интервал спектральной выборки в поперечном направлении;N - общее количество приемников в продольном направлении; М - общее количество приемников в поперечном направлении. Комплексный частотный параметр по Лапласу s настоящего изобретения также используется при данном преобразовании в спектральную область. В блоке 27 определяется система уравнений для очистки от ложных отражений преобразованных сейсмических данных для частоты из блока 26. Система уравнений будет определена с использованием частотного параметра по Лапласу s настоящего изобретения для обеспечения устойчивости при решении уравнений. В общем случае, когда местоположения приемников задаются произвольно, волновое поле в спектральной области не может быть определено в явном виде. Таким образом, процедура очистки от ложных отражений настоящего изобретения содержит решение системы уравнений. Для общего случая с произвольно заданными положениями приемников xR=(xR1,q, xR2,r, xR3,q,r) в спектральной области глубина расположения приемников задается однозначной функцией горизонтальных координат приемников. Конкретно, глубина расположения приемников записывается в виде xR3,q,r=xR3(xR1,q, xR2,r). В современных системах сбора морских сейсмических данных эта информация о положении приемников регистрируется. Операция очистки от ложных отражений может затем быть представлена в виде решения системы уравнений в которой матрица системы Kq,r,n,m задается выражением Здесь неизвестный вектор Pdghn,m представляет собой поле акустических волн с подавленными ложными отражениями на нулевой глубине, который задается в спектральной области в виде где известный вектор psctq,r представляет собой записанные сейсмические данные в пространственно-частотной области, заданный выражением а вертикальный коэффициент распространения sn,m представлен как где с представляет собой скорость распространения акустической волны в воде. Матрица системы Kq,r,n,m является функцией комплексного частотного параметра по Лапласу s. В блоке 28 процесс переходит к блоку 31 на фиг. 3 для решения системы уравнений, определенной в блоке 27. Использование полного комплексного частотного параметра по Лапласу s с ненулевой вещественной частью , взятой из выражения (4) вместо выражения (3), обеспечивает устойчивость при решении уравнений. На фиг. 3 представлена блок-схема, иллюстрирующая конечную часть второго варианта осуществления изобретения для очистки от ложных отражений данных морской сейсмической косы. В рамках конечной части продолжается начальный процесс, обсуждавшийся выше со ссылкой на блок-схему на фиг. 2. В блоке 31 нормальная система уравнений определяется по системе уравнений очистки от ложных отражений в блоке 27 на фиг. 2. Поскольку известный вектор psctq,r определен в пространственночастотной области Лапласа, а неизвестный вектор Pdghn,m определен в спектральной области, в данном варианте осуществления изобретения применяется решение методом наименьших квадратов в пространственно-частотной области. Обе части уравнения (10) умножаются на комплексно-сопряженную матрицу системы Kq,r;n',m', а результаты суммируются по всем значениям номеров приемников в продольном и поперечном направлениях q и r соответственно. После перестановки результатов последнего суммирования по q и r и суммирования по n и m получаем следующую систему нормальных уравнений: Здесь матрица An',m';n,m представляет собой известную величину, заданную выражением где n и n' являются спектральными числами в продольном направлении;m и m' - спектральные числа в поперечном направлении; черта над матрицей системы Kq,r;n',m', обозначает комплексное сопряжение. В уравнении (15) настоящего изобретения предлагается способ получения устойчивого решения данной линейной системы нормальных уравнений даже при наличии вырезов в спектре сейсмических данных путем использования комплексного частотного параметра по Лапласу. Геометрическое распределение расположений приемников не ограничивает объем настоящего изобретения. Настоящее изобретение применяется к простому случаю горизонтальных сейсмических кос с регулярной структурой, а также к общему случаю негоризонтальных сейсмических кос с нерегулярной структурой. Далее изобретение иллюстрируется сначала для простого случая, а затем для общего случая. Простой случай морской сейсмической косы при традиционном способе, рассматриваемый в качестве строго горизонтальной конфигурации сейсмических кос, является частным случаем общего случая,охватываемого данным изобретением. В этом простом случае горизонтальной косы расположение приемника по глубине xR3 является постоянным, а обращение (инверсия) системы уравнений сводится к стандартной процедуре очистки от ложных отражений. Система уравнений может быть решена обратным (инверсным) преобразованием Фурье. Конкретно, система уравнений, заданных в общем виде выражениями (10) и (11), становится более простой системой уравнений: Здесь очистка от ложных отражений для константы xR3 может быть выполнена в явном виде в спектральной области. Поле с подавленными ложными отражениями в спектральной области на нулевом уровне глубины pdghn,m, определенное выражением (12), может быть получено непосредственно как Для решения общей системы линейных уравнений в нормальных уравнениях, заданных выражением (15), существует много эффективных решателей, таких как прямые решатели и итеративные решатели. Очень эффективным итеративным решателем является итеративный метод сопряженных градиентов с предобусловливателем, основанный на минимизации расхождений между измеренными данными и моделируемым откликом. В еще одном варианте осуществления настоящего изобретения рассматривается выражение (10), а решение спектрального волнового поля с подавленными ложными отражениямиPdghn,m получают с использованием итеративной схемы сопряженных градиентов, такой что погрешность минимизируется. Как показано в блоке 32 на фиг. 3, нормальные уравнения из блока 31 предобусловлены с помощью физически обоснованного предобусловливателя. Предобусловливатель часто применяется в итеративных методах сопряженных градиентов, чтобы сделать уравнения хорошо обусловленными. Стандартные предобусловливатели, такие как многоуровневые ILU (Incomplete Lower Upper decomposition, неполное разложение матрицы на нижнюю или верхнюю матрицы), например, обычно основаны на математической структуре матрицы системы. К сожалению, для решения настоящей проблемы такие предобусловливатели, основанные на математической стороне проблемы, неэффективны. В одном из вариантов осуществления настоящего изобретения используется предобусловливатель, который основан на физической стороне решаемой проблемы, в качестве которой в общем случае выступает описанная выше сейсмическая коса. Негоризонтальная сейсмическая коса с нерегулярной структурой может рассматриваться в качестве возмущения горизонтальной сейсмической косы с регулярной структурой, поэтому обращение (инверсия) матрицы системы негоризонтальной сейсмической косы с нерегулярной структурой аппроксимируется обращением матрицы системы горизонтальной сейсмической косы с регулярной структурой. Здесь обращение матрицы системы эквивалентной горизонтальной сейсмической косы принимается в качестве предобусловливателя. Следует отметить, что если сейсмическая коса является горизонтальной, то итеративная схема сопряженных градиентов завершает свою работу за одну итерацию при использовании данного предобусловливателя. Для горизонтальной сейсмической косы система нормальных уравнений выражения (15) является диагональной, при этом решение получают непосредственно в виде В одном из вариантов осуществления настоящего изобретения используется предобусловливатель при этом выражение (16) принимает вид Для удобства выражение (10) можно записать в компактной форме где n теперь обозначает (n, m), a q обозначает (q, r). Предобусловливатель является функцией комплексного частотного параметра по Лапласу s. В блоке 33 итеративная схема метода сопряженных градиентов с предобусловливателем применяется для решения системы уравнений в выражении (23), предобусловленной нормальной системы уравнений из блока 32. На каждой итерации, обозначаемой числом k, получают приближенное значение спектрального волнового поля с подавленными ложными отражениями Pdgh(k)n с остаточной погрешностьюr(k)q. Таким образом, на k-й итерации k-я остаточная погрешность задается выражением Конкретный вариант осуществления настоящего изобретения для итеративной схемы метода сопряженных градиентов описан ниже со ссылкой на обсуждение блок-схемы на фиг. 4. В блоке 34 определяется, остались ли какие-либо частоты, выбранные в блоке 25 на фиг. 2. Если частоты остаются, процесс возвращается к блоку 25 на фиг. 2. Если частот не остается, процесс переходит к шагу 35. Таким образом, решение методом итераций сопряженных градиентов повторяется для каждой частоты f. В блоке 35 определяется, остались ли какие-либо выстрелы, выбранные в блоке 23 на фиг. 2. Если выстрелы остаются, процесс возвращается к блоку 23 на фиг. 2. Если выстрелов не остается, процесс переходит к шагу 36. Таким образом, решение методом итераций сопряженных градиентов повторяется для каждого выстрела f. В блоке 36 решения с подавленными ложными отражениями из блока 33 преобразуются обратно из спектральной области в пространственно-временную область. После получения волнового поля с подавленными ложными отражениями в спектральной области волновое поле с подавленными ложными отражениями pdgh=pdgh(xR1,q, xR2,q, 0; s) при нулевой глубине в пространственно-временной области получают с помощью выражения Волновое поле с подавленными ложными отражениями pdgh является функцией комплексного частотного параметра по Лапласу s. На фиг. 4 представлена блок-схема, иллюстрирующая вариант осуществления изобретения для итеративной схемы метода сопряженных градиентов с предобусловливателем, используемой для очистки от ложных отражений данных морской сейсмической косы. Это схема, на которую приведена ссылка в блоке 33 на фиг. 3. В блоке 41 выбирается нулевое возможное начальное значение для счетчика итераций k=0. В одном из вариантов осуществления это возможное значение выражается как В других вариантах осуществления могут использоваться различные варианты возможного начального значения для схемы метода сопряженных градиентов. Например, различные варианты могут включать, помимо прочего, решение, полученное для предыдущей частоты, или решение для другого положения выстрела. В блоке 42 рассчитывается возможное значение для счетчика итераций k=1. В одном из вариантов осуществления эти возможные значения выражаются как В блоке 43 счетчик итераций k увеличивается на 1. Таким образом, выражение k = k + 1 последовательно дает k=2, 3, , . В блоке 44 возможное значение для счетчика итераций k рассчитывается с использованием возможных значений для счетчика итераций k-1. В одном из вариантов осуществления это значение выражается как: В блоке 45 рассчитывается критерий погрешности для возможного текущего итерационного значения kth остаточной погрешности r(k)q из блока 44. В одном из вариантов осуществления критерий погрешности выражается как В блоке 46 определяется, удовлетворяется ли критерий погрешности в блоке 45. Если критерий погрешности не удовлетворяется, то итеративный процесс возвращается к блоку 43 для определения нового набора возможных значений. Если критерий погрешности удовлетворяется, то итеративный процесс завершается. При использовании метода обращения (инверсии) и применении комплексного частотного параметра по Лапласу s способ разложения для очистки от ложных отражений согласно настоящему изобретению можно переформулировать в различные области, включая, помимо прочего, в частотномедленностную область, или tau-медленностную область, или пространственно-временную область. Настоящее изобретение обсуждалось выше как способ только в целях иллюстрации, но оно может также быть реализовано в виде системы. Система в соответствии с настоящим изобретением предпочтительно реализуется посредством компьютеров, конкретно, цифровых компьютеров, наряду с использованием традиционного оборудования для обработки данных. Такое оборудование для обработки данных,хорошо известное специалистам, будет состоять из любой подходящей комбинации или сети оборудования для компьютерной обработки, включая, без ограничений, аппаратное обеспечение (процессоры, устройства временного и постоянного хранения и любое иное оборудование для компьютерной обработки),программное обеспечение (операционные системы, прикладные программы, библиотеки программ математических расчетов и любое иное подходящее программное обеспечение), схемы соединений (электрические, оптические, беспроводные и др.) и периферийное оборудование (входные и выходные устройства, такие как клавиатуры, позиционирующие устройства и сканеры; устройства отображения, такие как мониторы и принтеры; машиночитаемые накопители, такие как ленты, диски, накопители на жестких дисках и любое иное подходящее оборудование). В другом варианте осуществления настоящее изобретение может быть реализовано в виде описанного выше способа, реализуемого с помощью программируемого компьютера. В другом варианте осуществления настоящее изобретение может быть реализовано в виде компьютерной программы, хранящейся на машиночитаемом носителе, при этом логика программы обеспечивает выполнение программируемым компьютером описанного выше способа. В другом варианте осуществления настоящее изобретение может быть реализовано в виде машиночитаемого носителя с записанной на нем компьютерной программой, при этом логика программы обеспечивает выполнение программируемым компьютером описанного выше способа. Следует понимать, что изложенное выше представляет собой просто подробное описание конкретных вариантов осуществления настоящего изобретения и что многочисленные изменения, модификации и варианты осуществления, альтернативные изложенным, могут быть выполнены в соответствии с раскрытым здесь изобретением без выхода за его рамки. Поэтому предшествующее описание не ограничивает объем настоящего изобретения. Объем изобретения определяется только прилагаемой формулой изобретения и ее эквивалентами. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ морской сейсморазведки, включающий следующие шаги: получают и записывают сейсмические данные с помощью морских сейсмических кос, которые отражают характеристики формаций геологической среды, при этом сейсмические данные собирают посредством размещения ряда сейсмоприемников сверху оцениваемой области формаций геологической среды, причем сейсмоприемники формируют электрический и/или оптический сигнал в ответ на сейсмическую волну; преобразуют сейсмические данные из пространственно-временной области в спектральную область,используя комплексный частотный параметр по Лапласу и применяя итеративную схему метода сопряженных градиентов с использованием физически обоснованного предобусловливателя к преобразованным сейсмическим данным, чтобы получить решение нормальной системы уравнений методом наименьших квадратов для системы уравнений очистки от ложных отражений; и подвергают решение обратному преобразованию в пространственно-временную область для получения сейсмических данных с подавленными ложными отражениями, обеспечивающих построение изображений геологической среды. 2. Способ по п.1, отличающийся тем, что использование комплексного частотного параметра по Лапласу включает следующие шаги: получают сейсмические данные, записанные с помощью морской сейсмической косы; определяют ненулевую вещественную часть комплексного частотного параметра по Лапласу; получают записанные сейсмические данные для активации сейсмического источника; используют комплексный частотный параметр по Лапласу, чтобы преобразовать записанные сейсмические данные для активации сейсмического источника из пространственно-временной области в пространственно-частотную область; получают преобразованные сейсмические данные для некоторой частоты и используют комплексный частотный параметр по Лапласу, чтобы преобразовать преобразованные сейсмические данные для некоторой частоты из пространственно-временной области в спектральную область. 3. Способ по п.2, отличающийся тем, что комплексный частотный параметр по Лапласу s задают выражением где j представляет собой мнимую единицу -1;- круговая частота;f - частота;- определенная ненулевая вещественная часть комплексного частотного параметра по Лапласу s. 4. Способ по п.3, отличающийся тем, что использование комплексного частотного параметра по Лапласу для преобразования записанных сейсмических данных включает преобразование рассеянного поля акустических волн из пространственно-временной области (xR1, xR2, xR3, t) в пространственночастотную область (xR1, xR2, xR3, s) с использованием преобразования Лапласа, при этом где psct представляет собой рассеянное поле акустических волн в пространственно-частотной области;xR3,q,r - глубина расположения приемника как функция xR1,q и xR2,r. 5. Способ по п.4, отличающийся тем, что использование комплексного частотного параметра по Лапласу для преобразования преобразованных сейсмических данных включает преобразование рассеянного волнового поля psct из пространственно-частотной области (xR1, xR2, xR3, s) в спектральную областьPsctn,m - представляет собой рассеянное поле акустических волн в спектральной области;xR1 - интервал выборки для приемника в продольном направлении;xR2 - интервал выборки для приемника в поперечном направлении;jsn - спектральный параметр Фурье в продольном направлении;(s) - интервал спектральной выборки в продольном направлении;jsm - спектральный параметр Фурье в поперечном направлении;(s) - интервал спектральной выборки в поперечном направлении;N - общее количество приемников в продольном направлении; М - общее количество приемников в поперечном направлении. 6. Способ по п.5, отличающийся тем, что определяют систему уравнений для очистки от ложных отражений преобразованных сейсмических данных для некоторой частоты; определяют нормальную систему уравнений для системы уравнений для очистки от ложных отражений и применяют физически обоснованный предобусловливатель к нормальной системе уравнений. 7. Способ по п.6, отличающийся тем, что определение системы уравнений для очистки от ложных отражений включает применение следующего выражения: в котором матрица системы Kq,r,n,m задается выражением 8. Способ по п.7, отличающийся тем, что определение нормальной системы включает применение следующего выражения:m и m' - спектральные числа в поперечном направлении; черта над матрицей системы n' обозначает комплексное сопряжение. 9. Способ по п.7, отличающийся тем, что предобусловливатель содержит обращение матрицы системы для горизонтальной сейсмической косы. 10. Способ по п.9, отличающийся тем, что применение предобусловливателя включает применение следующего выражения: 11. Способ по п.1, отличающийся тем, что применение итеративной схемы метода сопряженных градиентов включает следующие шаги: выбирают возможное начальное значение для счетчика итераций k=0 в соответствии со следующими выражениями: рассчитывают возможные значения для счетчика итераций k=1 в соответствии со следующими выражениями: выполняют следующие действия, пока критерий погрешности не будет удовлетворен: увеличивают значение счетчика итераций k=k+1 и рассчитывают возможные значения для счетчика итераций k, используя предыдущие возможные значения, рассчитанные для k-1, в соответствии со следующими выражениями: 12. Способ по п.11, отличающийся тем, что критерий погрешности задают выражением 13. Способ по п.5, отличающийся тем, что обратное преобразование решения включает применение следующего выражения: 14. Машиночитаемый носитель с записанной на нем компьютерной программой, обеспечивающей выполнение способа по п.1. 15. Носитель по п.14, отличающийся тем, что использование комплексного частотного параметра по Лапласу включает следующие шаги: получение сейсмических данных, записанных с помощью морской сейсмической косы; определение ненулевой вещественной части комплексного частотного параметра по Лапласу; получение записанных сейсмических данных для активации сейсмического источника; использование комплексного частотного параметра по Лапласу, чтобы преобразовать записанные сейсмические данные для активации сейсмического источника из пространственно-временной области в пространственно-частотную область; получение преобразованных сейсмических данных для некоторой частоты и использование комплексного частотного параметра по Лапласу, чтобы преобразовать преобразованные сейсмические данные для некоторой частоты из пространственно-временной области в спектральную область. 16. Носитель по п.15, отличающийся тем, что комплексный частотный параметр по Лапласу s задан выражением где j представляет собой мнимую единицу -1;- круговая частота;f - частота;- определенная ненулевая вещественная часть комплексного частотного параметра по Лапласу s. 17. Носитель по п.16, отличающийся тем, что использование комплексного частотного параметра по Лапласу для преобразования записанных сейсмических данных включает следующие шаги: преобразование рассеянного поля акустических волн из пространственно-временной области (xR1,R где psct представляет собой рассеянное поле акустических волн в пространственно-частотной области;xR3,q,r - глубина расположения приемника как функция xR1,q и xR2,r. 18. Носитель по п.17, отличающийся тем, что использование комплексного частотного параметра по Лапласу для преобразования преобразованных сейсмических данных включает следующие шаги: преобразование рассеянного волнового поля psct из пространственно-частотной области (xR1, xR2,RPsctn,m представляет собой рассеянное поле акустических волн в спектральной области;xR1 - интервал выборки для приемника в продольном направлении;xR2 - интервал выборки для приемника в поперечном направлении;jsn - спектральный параметр Фурье в продольном направлении;(s) - интервал спектральной выборки в продольном направлении;jsm - спектральный параметр Фурье в поперечном направлении;(s) - интервал спектральной выборки в поперечном направлении;N - общее количество приемников в продольном направлении; М - общее количество приемников в поперечном направлении. 19. Носитель по п.18, отличающийся тем, что логика программы обеспечивает выполнение программируемым компьютером следующих шагов: определение системы уравнений для очистки от ложных отражений преобразованных сейсмических данных для некоторой частоты; определение нормальной системы уравнений для системы уравнений для очистки от ложных отражений и применение физически обоснованного предобусловливателя к нормальной системе уравнений. 20. Носитель по п.19, отличающийся тем, что определение системы уравнений для очистки от ложных отражений включает применение следующего выражения: в котором матрица системы Kq,r;n,m задается выражением 21. Носитель по п.20, отличающийся тем, что определение нормальной системы включает применение следующего выражения:m и m' - спектральные числа в поперечном направлении; черта над матрицей системы n' обозначает комплексное сопряжение. 22. Носитель по п.20, отличающийся тем, что предобусловливатель содержит обращение матрицы системы для горизонтальной сейсмической косы. 23. Носитель по п.22, отличающийся тем, что применение предобусловливателя включает применение следующего выражения: 24. Носитель по п.14, отличающийся тем, что применение итеративной схемы метода сопряженных градиентов включает следующие шаги: выбор возможного начального значения для счетчика итераций k=0 в соответствии со следующим выражением: расчет возможного значения для счетчика итераций k=1 в соответствии со следующим выражением: выполнение следующих действий, пока критерий погрешности не будет удовлетворен: увеличение значения счетчика итераций k=k+1; расчет возможного значения для счетчика итераций k с использованием предыдущих возможных значений, рассчитанных для k-1, в соответствии со следующими выражениями: 25. Носитель по п.24, отличающийся тем, что критерий погрешности задан выражением 26. Носитель по п.18, отличающийся тем, что обратное преобразование решения включает применение следующего выражения: