Способ определения волнового поля, свободного от ложных отраженных сигналов, излучаемого морским сейсмоисточником (варианты)

СПОСОБ ОПРЕДЕЛЕНИЯ ВОЛНОВОГО ПОЛЯ, СВОБОДНОГО ОТ ЛОЖНЫХ ОТРАЖЕННЫХ СИГНАЛОВ, ИЗЛУЧАЕМОГО МОРСКИМ СЕЙСМОИСТОЧНИКОМ В изобретении сейсмические данные получают путем активации первого сейсмоисточника в первый момент времени и одного или нескольких сейсмоисточников в характерные для них собственные моменты времени относительно времени записи сигналов, при этом источники существенно связаны и расположены на различных глубинах. Определяют первое волновое поле,которое возникло бы, если бы первый источник был активирован в момент времени, выбранный относительно времени начала выполнения записей, и корректируют по времени с учетом глубины воды. Определяют одно или несколько дополнительных волновых полей, которые возникли бы,если бы один или каждый из нескольких дополнительных источников были активированы в указанный момент времени, выбранный по отношению к указанному времени начала выполнения записей, и корректируют по времени с учетом глубины одного или нескольких дополнительных источников в воде. Первое волновое поле объединяют с одним или несколькими дополнительными волновыми полями с целью определения волнового поля с подавленными ложными отражениями,соответствующего активации единственного сейсмоисточника. Хегна Стиан (NO), Паркес Грегори Эрнст (GB) Хмара М.В., Рыбаков В.М. (RU)(71)(73) Заявитель и патентовладелец: ПГС ГЕОФИЗИКАЛ АС (NO) Область техники, к которой относится изобретение Настоящее изобретение в целом относится к области сбора и обработки морских сейсмических данных. Конкретнее, изобретение относится к способам уменьшения эффекта отражения от водной поверхности волн, излучаемых морскими сейсмоисточниками. Предшествующий уровень техники При проведении сейсмических исследований сейсмические данные получают, направляя акустическую волну в геологическую среду вблизи поверхности Земли и регистрируя акустическую волну, которая отражается от границ между различными пластами формаций геологической среды. Акустическая волна отражается, когда существует разность акустических импедансов на границе между соседними пластами. Сигналы, представляющие зарегистрированную акустическую волну, интерпретируются, чтобы вывести заключение о строении и составе формаций геологической среды. При проведении морских сейсморазведочных исследований сейсмоисточник, например пневмопушка или группа таких пневмопушек, используется обычно для направления акустической волны в формации, расположенные ниже подошвы водного слоя. Пневмопушка или группа пневмопушек активируется в воде на выбранной глубине, как правило, в то время как пневмопушка или группа буксируется судном. Это же или другое судно буксирует в воде один или несколько сейсмоприемных кабелей, называемых "сейсмоприемными косами". В общем случае сейсмоприемная коса тянется за судном вдоль направления, в котором она буксируется. Как правило, сейсмоприемная коса включает набор гидрофонов, расположенных в пространственно разнесенных известных местах по длине кабеля. Гидрофоны, как известно специалистам, представляют собой датчики, которые генерируют оптический или электрический сигнал, соответствующий давлению воды или временному градиенту (dp/dt) давления в воде. Судно, которое буксирует одну или несколько сейсмоприемных кос, как правило, везет записывающее оборудование для выполнения индексированных по времени записей сигналов, генерируемых гидрофонами в ответ на зарегистрированную акустическую волну. Запись сигналов, как объяснялось ранее, подвергается обработке, чтобы вывести заключение о строении и составе формаций геологической среды, находящихся ниже места проведения сейсморазведки. Морские сейсмические данные включают эффект, который ограничивает точность выводов о строении и составе формаций геологической среды. Этот эффект, известный как образование волн-спутников со стороны источника, возникает вследствие того, что вода имеет плотность и скорость распространения волн давления, существенно отличающиеся от характеристик воздуха, находящегося над поверхностью воды. Образование волн-спутников источника можно понять следующим образом. При активировании пневмопушки или группы пневмопушек акустическая волна излучается в основном во внешнем направлении от пневмопушки или группы пневмопушек. Половина энергии волны распространяется в нижнем направлении и проходит через подошву водного слоя в формации геологической среды. Другая половина энергии акустической волны распространяется в верхнем направлении относительно пневмопушки или группы, при этом большая часть этой энергии отражается от поверхности воды, после чего распространяется вниз. Отраженная акустическая волна отстает по времени, а также сдвигается по фазе приблизительно на 180 относительно прямой нисходящей акустической волны. Отраженная от поверхности нисходящая акустическая волна обычно называется "волной-спутником". Волна-спутник интерферирует с полем прямой нисходящей волны, вызывая конструктивную (усиливающую) интерференцию в некоторых частях диапазона частот и деструктивную (ослабляющую) интерференцию в других частях диапазона частот. Деструктивная интерференция создает последовательность провалов в спектре, равномерно разнесенных по частоте, включая провал на нулевой частоте (0 Гц). Частоты этих провалов в зарегистрированном акустическом сигнале связаны с глубиной, на которой располагается пневмопушка или группа пневмопушек, как хорошо известно специалистам. Эффект образования волн-спутников со стороны источника, как правило, называется "волной-спутником источника". Сейсмические волны, излучаемые источником, ослабевают с расстоянием, пройденным волной, изза геометрического расхождения, потерь при распространении и поглощения. Поглощение энергии в высокочастотной части диапазона в большей степени по сравнению с низкочастотной энергией хорошо известно специалистам. Поэтому для более глубокого проникновения желательно максимально увеличить энергию, излучаемую источником на низких частотах. Поскольку волна-спутник источника создает провал в записанном сигнале на нулевой частоте, она ограничивает энергию в низкочастотной части сигнала. Количество энергии на низких частотах можно увеличить, буксируя источники на большей глубине. Однако при этом провалы, вызванные волной-спутником в спектре, возникают на более низких частотах,и тем самым ограничивают высокочастотные части спектра, необходимые для построения высокоразрешающего изображения для менее глубоко залегающих объектов исследований. Кроме того, при использовании пневмопушки(-ек) в качестве сейсмоисточника основная частота пушки(-ек) возрастает с увеличением глубины. Таким образом, увеличение энергии на низких частотах при буксировании пневмопушки(-ек) на большей глубине из-за волны-спутника сводится на нет ростом основной частоты пневмопушки(-ек). Известный способ повышения уровня сигнала, излучаемого источником в полосе частот при использовании пневмопушки(-ек), состоит в увеличении общего объема воздуха, выпускаемого пневмо-1 021738 пушкой(-ами), и (или) в увеличении рабочего давления. Однако максимальный объем воздуха, который может быть выпущен при каждом выстреле, и максимальное воздушное давление воздуха ограничены имеющимся оборудованием источника и системой подачи воздуха. Изменение максимального объема и давления может представлять собой очень дорогостоящую и трудоемкую задачу. Кроме того, повышение мощности источника может оказать воздействие на морскую флору и фауну. Таким образом, максимальное использование сигнала, излучаемого источником, может представлять собой большую ценность и уменьшить необходимость в повышении уровня энергии, излучаемой источником. Благодаря выделению полей восходящих (волна-спутник) и прямых нисходящих волн из источника, эффекты волны-спутника источника подавляются, а сигнал вблизи всех провалов, вызванных волной-спутником, усиливается,включая сигнал в области низких частот вблизи от провала на частоте 0 Гц. Эти разделенные волновые поля могут быть также сдвинуты во времени относительно поверхности моря или общей опорной глубины с использованием известной глубины источника(-ов). Затем, за счет применения фазового сдвига на 180 к волне-спутнику, разделенные волновые поля можно конструктивно (с усилением) суммировать. При этом используется почти вся энергия, излучаемая источником. Известный специалистам способ выделения волны-спутника источника описан в работе М. Egan etal., Full deghosting of OBC data with over/under source acquisition, 2007 Annual Meeting, San Antonio, TX,Society of Exploration Geophysicists. Методика, описанная в публикации Egan и др., включает буксирование первого сейсмоисточника на первой глубине в воде, а также буксирование второго сейсмоисточника на второй глубине в воде. Источники представляют собой пневмопушки или их группы. Второй источник также буксируется на выбранном расстоянии позади первого источника. Первый источник активируется, и выполняется запись сейсмических сигналов, соответствующих активации первого источника. После того как буксирующее судно переместится так, что второй источник окажется расположенным, по существу, в том же геодезическом положении, где первый источник находился в момент активации, активируется второй источник, и снова выполняется запись сейсмических сигналов. Набор сейсмических данных с подавленными волнами-спутниками (или "очищенных от ложных отражений") получают с использованием методики, более подробно описанной в публикации Egan и др. Одна из основных проблем при использовании методики с источниками, расположенными сверху/снизу, которая описана в упомянутой публикации Egan и др., состоит в том, что количество позиций для выстрела получается в два раза меньшим по сравнению с традиционными способами активации источника, что приводит к половинному перекрытию. Другая проблема при использовании этого способа,если сейсмоприемники буксируются за судном и, таким образом, перемещаются от одной позиции для выстрела до другой, состоит в том, что, пока происходит активация источников на разных глубинах,приемники проходят значительное расстояние между позициями. Чтобы сохранить количество позиций для выстрела и перекрытие таким же, как при традиционном сборе морских сейсмических данных, и свести к минимуму различие в позициях приемников во время активации источников на разных глубинах,желательно иметь способ выделения волны-спутника источника, который обеспечит активацию источников, буксируемых на разных глубинах, во время регистрации сейсмограммы каждого выстрела (ОПВ,общего пункта возбуждения). Специалистам также известно использование двух или нескольких сейсмоприемных кос, буксируемых на различных глубинах, для получения записанного сигнала, который, по существу, свободен от волн-спутников, вызванных отражениями от водной поверхности в непосредственной близости от сейсмоприемных кос; см., например, публикацию Posthumus, В. [1993] Deghosting of twin streamer configuration, Geophysical Prospecting, vol. 41, p. 267-286. Результаты непосредственной адаптации способа, описанного в публикации Posthumus, к очистке от ложных отражений источников, расположенных "сверху/снизу", когда источники на различных глубинах активируются с временными задержками в тех же сейсмических записях, оказались неточными. Соответственно сохраняется потребность в усовершенствованных технологиях для уменьшения воздействия образования волн-спутников со стороны источника на морские сейсмические данные. Сущность изобретения В одном из вариантов осуществления настоящего изобретения предлагается способ определения волнового поля с подавленными ложными отражениями, генерируемого морским сейсмическим источником, по записям сигналов сейсмических данных, полученных путем активации первого сейсмоисточника в первый момент времени и одного или более дополнительных сейсмоисточников в характерные для них собственные моменты времени относительно времени индексации записей сигналов, при этом источники, по существу, связаны и расположены на различных глубинах в водоеме. По записям сигналов определяют первое волновое поле, которое возникло бы, если бы первый источник был активирован в момент времени, выбранный относительно времени начала выполнения записей сигналов. Первое волновое поле корректируют по времени с учетом глубины первого источника в воде. По записям сигналов определяют одно или более дополнительных волновых полей, которые возникли бы, если бы один или каждый из дополнительных источников были активированы в указанный момент времени, выбранный относительно времени начала выполнения записей сигналов. Одно или более дополнительных волновых полей корректируют по времени с учетом глубины одного или более дополнительных источников в воде. Первое волновое поле объединяют с одним или более дополнительными волновыми полями с целью определения волнового поля с подавленными ложными отражениями, соответствующего активации единственного сейсмоисточника. Прочие особенности и преимущества настоящего изобретения будут ясны из нижеследующего описания и прилагаемой формулы изобретения. Перечень фигур чертежей Настоящее изобретение и его преимущества легче понять из следующего подробного описания и прилагаемых чертежей: на фиг. 1 показано получение сейсмических данных в поперечном разрезе с целью демонстрации примера расположения сейсмоисточников; на фиг. 2 показано получение сейсмических данных в горизонтальной проекции с целью демонстрации примера расположения сейсмоприемных кос; на фиг. 3 А-3 С показаны соответственно примеры форм импульсов двух источников, активируемых в различные моменты времени и на различных глубинах, разложение с помощью методики, относящейся к известному уровню техники (способ, описанный в публикации Posthumus [1993]), и разложение с помощью способа согласно настоящему изобретению; на фиг. 4 А и 4 В показаны соответственно спектры волновых полей с подавленными ложными отражениями, представленных на фиг. 3, при использовании способа Posthumus и способа настоящего изобретения; на фиг. 5 А и 5 В показаны соответственно спектральная погрешность в волновых полях с подавленными ложными отражениями при использовании способа Posthumus и способа настоящего изобретения; на фиг. 6 А и 6 В показан соответственно спектр неотфильтрованного и точного знаменателя в уравнении согласно способу настоящего изобретения и спектр полосового фильтра, применяемого к двум последним членам в таком знаменателе; на фиг. 7 А и 7 В показан соответственно спектр знаменателя в уравнениях, используемых в способеPosthumus, и спектр знаменателя в варианте с полосовым фильтром согласно способу настоящего изобретения; на фиг. 8 А и 8 В показан соответственно спектр двух последних членов знаменателя в уравнении согласно способу настоящего изобретения без полосовой фильтрации и с применением полосового фильтра; на фиг. 9 представлена блок-схема, иллюстрирующая вариант осуществления настоящего изобретения для определения волнового поля с подавленными ложными отражениями, генерируемого морским сейсмическим источником. Хотя настоящее изобретение будет описано в связи с предпочтительными вариантами осуществления, понятно, что изобретение не ограничивается только ими. Напротив, изобретение предназначено для того, чтобы охватить все альтернативные варианты, модификации и эквиваленты, которые могут быть включены в объем настоящего изобретения согласно прилагаемой формуле изобретения. Сведения, подтверждающие возможность осуществления изобретения На фиг. 1 показан пример схемы получения сейсмических данных в поперечном разрезе в соответствии с настоящим изобретением. Сейсморазведочное судно 10 движется по поверхности 11 А водоема 11, например озера или океана. Судно 10, как правило, везет оборудование, обозначенное в целом позицией 12 и для удобства называемое "записывающей системой". Записывающая система 12 может включать устройства (не показанные отдельно) для выборочной активации сейсмоисточников 14, 16 (объясняемых ниже) для записи сигналов, генерируемых датчиками или приемниками 20 (объясняемыми ниже) в ответ на сейсмическую волну, направляемую в воду 11 (источниками 14, 16) и тем самым в формации 19, 21, расположенные ниже подошвы водного слоя 13, и для определения геодезического положения судна 10, сейсмоисточников 14, 16 и каждого из набора сейсмических датчиков или приемников 20 в любой момент времени. Судно 10 показано буксирующим два сейсмоисточника 14, 16 только в целях наглядности. Изобретение предназначено для использования с двумя или несколькими сейсмоисточниками и не ограничивается только двумя сейсмоисточниками. Сейсмоисточники 14, 16 могут принадлежать к любому типу морских источников энергии, включая, без ограничений, пневмопушки и гидропушки или группы таких источников энергии. В примере, показанном на фиг. 1, источники 14, 16 буксируются за судном 10, по существу, на одинаковом расстоянии и на различных глубинах в воде 11. Такое расположение источников называется "сверху/снизу". В других примерах источники 14, 16 могут буксироваться другим судном(не показано) или могут находиться в неподвижном положении (при условии, что глубины различны, как показано на фиг. 1). Поэтому буксирование исследовательским судном 10 источников 14, 16 не ограничивает объем настоящего изобретения. Судно 10 также показано буксирующим сейсмоприемную косу 18, включающую набор установленных на ней сейсмоприемников 20, пространственно разнесенных в продольном направлении. Однако настоящее изобретение в целом относится к сейсмоисточникам, и поэтому может использоваться с любым типом расположения сейсмоприемников. В различных вариантах осуществления изобретение может использоваться с буксируемыми сейсмоприемными косами в любой конфигурации, кабелями, размещаемыми на морском дне, датчиками, размещаемыми в скважинах, и т.д., и с любым типом приемного датчика, включая, без ограничений, датчики давления, датчики временного градиента давления, датчики скорости, акселерометры и т.д. или любую их комбинацию. В выбранный момент времени при работе системы сбора данных, показанной на фиг. 1, система сбора данных 12 активирует первый из сейсмоисточников, например источник 14. Выбранный момент времени может индексироваться по моменту начала записи откликов приемников или любому другому моменту начала отсчета, связанному с временем записи сигналов, и, как правило, характеризуется первой выбранной задержкой относительно начала записи. Сейсмическая волна, генерируемая первым источником 14, распространяется от него во внешнем направлении, как показано позициями 22 и 24. Часть сейсмической волны распространяется вниз, 22, где она отражается от границ акустических импедансов, например подошвы водного слоя 13, и от границ 15, 17 между различными формациями 19, 21. Для наглядности на фиг. 1 показаны только отражения от подошвы водного слоя. Восходящие части сейсмической волны от первого источника 14, обозначенные позицией 24, отражаются от поверхности воды 11 А, как показано на фиг. 1. Записывающая система 12 выполнена с возможностью активации второго сейсмоисточника, например источника 16, в конце второй временной задержки относительно начала записи сейсмических данных, или, альтернативно, после выбранного периода времени перед или после активации первого источника 14. Сейсмическая волна, распространяющаяся во внешнем направлении от второго источника 16, движется по той же траектории, что и сейсмическая волна от первого источника 14. В настоящем изобретении описанные выше временные задержки выбирают таким образом, чтобы сейсмическая волна от обоих источников регистрировалась приемниками 20 в каждой записи данных (называемой "сейсмограммой ОПВ"). Предполагается, что время активации источников на различных глубинах на каждой из сейсмограмм ОПВ может отличаться произвольным, полупроизвольным или систематическим образом. На фиг. 2 схема расположения согласно фиг. 1 показана в горизонтальной проекции, чтобы проиллюстрировать буксирование набора сейсмоприемных кос 18, пространственно разнесенных в поперечном направлении. Сейсмоприемные косы 18 можно удерживать в их относительном поперечном и продольном положении по отношению к судну 10 с помощью буксировочного оборудования 23, принадлежащего к типам, хорошо известным специалистам. Как и в предыдущем случае, судно 10 показано буксирующим два сейсмоисточника 14, 16 только в целях наглядности. Изобретение предназначено для использования с двумя или несколькими сейсмоисточниками и не ограничивается только двумя сейсмоисточниками. Кроме того, на фиг. 2 показано, что первый источник 14 и второй источник 16 могут быть смещены в поперечном направлении и/или в продольном направлении в других вариантах осуществления, чтобы избежать воздействия диспергированного воздуха в воде 11 от первого источника 14 на восходящую сейсмическую волну от второго источника 16, и наоборот, в том случае, когда источники 14, 16 представляют собой пневмопушки или их группы. Предполагается, что поперечное и/или продольное смещение составляет только несколько метров, так что источники 14, 16 обеспечивают сейсмическую волну, эквивалентную той, которая возникла бы, если бы источники 14, 16 находились в одной вертикальной плоскости и на одинаковом продольном расстоянии позади судна или, выражая это иными словами, по существу, в одном и том же геодезическом положении. Благодаря предотвращению присутствия диспергированного воздуха над источниками при их активации, влияние водной поверхности (11 А на фиг. 1) на сейсмическую волну, излучаемую каждым источником 14, 16, будет, с учетом глубины воды,по существу, таким же, как влияние водной поверхности на волну, излучаемую другим из источников 14,16. Активация источников и запись сигналов, объяснявшиеся выше, могут многократно повторяться,пока судно 10, источники 14, 16 и сейсмоприемные косы 18 движутся по воде 11. Как объяснялось выше,каждая такая сейсмограмма ОПВ для каждого приемника 20 будет включать сигналы, соответствующие сейсмической волне, произведенной как первым источником 14, так и вторым источником 16. На фиг. 1 и 2 представлены два источника, каждый из которых работает на различной глубине. Это сделано только для наглядности. Настоящий способ не ограничивается двумя источниками, работающими на различной глубине, но может включать множество источников, работающих на разных глубинах. Сейсмические данные, полученные при проведении сейсморазведки и характерные для геологической среды, обрабатываются для получения информации, относящейся к геологическому строению и свойствам формаций геологической среды в исследуемой зоне. Обработанные сейсмические данные обрабатываются для отображения и анализа потенциального группового состава углеводородов в этих подземных формациях. Целью обработки сейсмических данных является извлечение из сейсмических данных максимального объема информации о подземных формациях для получения адекватного изображения геологической среды. При определении участков геологической среды, где существует вероятность обнаружения залежей нефти, крупные средства расходуются на сбор, обработку и интерпретацию сейсмических данных. Процесс построения поверхностей отражающих горизонтов, определяющих исследуемые подземные пласты, по записанным сейсмическим данным позволяет получить изображение среды по глубине или времени. Изображение строения геологической среды получают с целью предоставления интерпретатору возможности выбора участков, где с наибольшей вероятностью могут находиться залежи нефти. Чтобы удостовериться в наличии нефти, необходимо пробурить скважину. Бурение скважин для определения присутствия нефтяных залежей представляет собой чрезвычайно дорогостоящую и трудоемкую задачу. По этой причине сохраняется потребность в улучшении обработки и отображения сейсмических данных для получения изображения строения геологической среды, которое повысит возможности интерпретатора оценивать вероятность существования скопления нефти на определенном участке геологической среды независимо от того, выполняется ли интерпретация компьютером или человеком. Обработка и отображение собранных сейсмических данных способствует принятию более точных решений при определении необходимости и места бурения, тем самым снижая риск бурения непродуктивных скважин. Объяснив методику сбора данных, которая может использоваться с настоящим изобретением, поясним теперь способ настоящего изобретения в отношении разложения волнового поля. В следующем примере варианта осуществления два источника на двух глубинах применяются только для простоты иллюстрации. Настоящее изобретение охватывает использование любого количества источников на различных глубинах. Сначала во временной области операторы волн-спутников gn (функции, которые определяют форму импульса волны-спутника соответствующих сейсмоисточников) могут быть определены для каждого источника следующим образом. Для первого сейсмоисточника В приведенных выше выражениях d1 - это рабочая глубина первого источника; d2 - рабочая глубина второго источника;- угол испускания сейсмической волны, распространяющейся от источника, по отношению к вертикали;- акустическая скорость в воде, которая составляет приблизительно 1500 м/с;оператор -функции Дирака. Пусть w(t) представляет волновое поле, которое было бы записано на выбранном приемнике от первого источника при t1= 0 с (в момент начала записи), если бы не было ложного отражения со стороны источника и если бы форма импульса источника представляла собой идеальный -импульс. Фактически записанный сигнал r(t), возникающий в результате активации первого источника, активируемого в некоторый фактический момент времени после начала записи t=t1, и второго источника, активируемого в момент времени t=t2, может быть выражен следующим уравнением: где n(t) представляет помехи, s1(t) и s2(t) представляют свободные от ложных отражений формы импульсов первого и второго источников соответственно, апредставляет свертку во временной области. Сдвиг r(t) назад во времени на t1 секунд в уравнении (1) и, таким образом, внесение в него поправки на фактическое время активации первого источника относительно начала записи (время t=0) в сейсмограмме ОПВ позволяет получить выражение для первого скорректированного по времени сигнала o1(t) Преобразование уравнения (2 а) в частотную область позволяет получить следующее выражение: в котором N1 представляет собой преобразование n(t+t1), и В уравнении (2b) и последующих выражениях написание каждой из переменных прописной буквой представляет ее преобразование из временной области в частотную область. Сдвиг r(t) назад во времени на t2-(d2-d1)cos/ секунд в уравнении (1) и, таким образом, внесение в него поправки на время активации второго источника вместе с компенсацией его более близкого расположения относительно отражающих горизонтов по сравнению с первым источником дает следующее выражение для второго скорректированного по времени сигнала o2(t): Преобразование уравнения (3 а) в частотную область позволяет получить следующее выражение: Если создаются две копии (как показано здесь, хотя можно создать и большее количество копий) записанных волновых полей и O1 определяется как записанные данные, скорректированные на время активации первого источника, а О 2 - как записанные данные, скорректированные на время активации второго источника (как объяснялось выше), при этом в обе записи вносятся поправки на величину небольшой разности между временами вступления, вызванной различными рабочими глубинами каждого из двух источников, результат, как объяснялось выше, может быть выражен следующим уравнением: где T2 представляет собой время активации второго источника минус время активации первого источника и плюс разность времен, которая обусловлена различными рабочими глубинами источников, aT1, аналогичным образом, время активации первого источника минус время активации второго источника и плюс разность времен, которая обусловлена различными рабочими глубинами источников. По существу, уравнение (4 а) представляет данные первого источника, скорректированные по времени, и может быть описано как волновое поле, создаваемое первым источником, плюс волновое, создаваемое вторым источником, с временной задержкой, зависящей от разности времени активации двух источников и разности глубин источников, плюс помехи. Уравнение (4b) аналогичным образом представляет записанные данные второго источника, скорректированные по времени. Умножение уравнений (4 а) и (4b) на комплексно сопряженную функцию волн-спутников для соответствующих глубин источников дает следующие уравнения: в которыхпредставляет комплексно сопряженную функцию соответствующих операторов волнспутников в частотной области. Сложение уравнений (5 а) и (5b) позволяет получить следующее выражение: Из уравнения (6) можно вывести следующее выражение для вычисления волнового поля с подавленными ложными отражениями со стороны источника с поправкой на форму импульсов источника В способе, предложенном в публикации Posthumus (ссылка на которую приводится в разделе"Предшествующий уровень техники" настоящего изобретения) для очистки от ложных отражений с помощью приемников, расположенных сверху/снизу, а не источников, расположенных сверху/снизу и активируемых в пределах тех же сейсмических записей на различных глубинах, для вычисления волнового поля с подавленными ложными отражениями на сейсмоприемных косах используется следующее выражение: Сравнение уравнения (7) с уравнением (8) показывает, что эти два уравнения очень похожи. Различие между двумя уравнениями заключается в знаменателе, при этом вычисление волнового поля с подавленными ложными отражениями со стороны источника с поправкой на форму импульсов источника включает свободные от ложных отражений формы импульсов источников и два дополнительных члена. В составе двух дополнительных членов G1G2 и G2G1 эквивалентны взаимной корреляции двух функций волн-спутников во временной области, a exp(-iT1) и exp(-iT1) представляют собой операторы временной задержки для введения поправок на разность между моментами времени активации и рабочими глубинами источников (T1 и Т 2, как объяснялось выше). Установлено, что использование двух дополнительных членов, как объясняется выше, позволяет существенно уменьшить погрешность при разложении сигнала в отличие от непосредственного применения способа Posthumus для приемников, расположенных сверху/снизу и активируемых для записи на одной и той же сейсмограмме ОПВ. Кроме того, если два (или несколько) источников, работающих на различных глубинах, не рассчитаны специально на то, чтобы быть близкими к идентичным, формы импульсов источников необходимо включать в выражение, чтобы учесть разность форм импульсов, описываемую уравнением (7). Поскольку эти взаимно коррелированные функции волн-спутников сдвинуты вперед и назад по времени, в спектре знаменателя будут наблюдаться пики и провалы, связанные с разностью времен активации двух источников. Эти пики и провалы могут потенциально вызывать усиление помех в некоторых частях спектра, особенно на высоких и низких частотах. Таким образом, при вычислении волнового поля с подавленными ложными отражениями желательно ослаблять или фильтровать эти два дополнительных члена. Уравнения (7) и (8) аналогичны, за исключением этих двух дополнительных членов в знаменателе уравнения 7. Ослабление двух дополнительных членов в уравнении (7) позволяет приблизить результат применения уравнения (8) к данным метода "сверху/снизу" при возбуждении двух источников с записью на одних и тех же сейсмограммах с различными моментами времени возбуждения. Однако даже если для того, чтобы избежать усиления помех в нижней и верхней частях спектра, необходима полосовая фильтрация двух дополнительных членов в уравнении (7), вносимая погрешность будет очень невелика при условии, что фильтрация выполняется в тех частях спектра, где взаимная корреляция двух волнспутников обладает ограниченной энергией, и/или в частях спектра с ограниченной энергией сейсмических волн. Такой подвергнутый фильтрации вариант уравнения (7) можно представить в следующем виде:a F представляет собой полосовой фильтр. Уравнения (7) и (9) можно легко распространить на случаи возбуждения трех или более источников, расположенных на трех или большем количестве глубин, с записью на одних и тех же сейсмограммах с различными временными задержками. Поскольку полностью отделить помехи от сигнала, как правило, невозможно, только первый член в правой части уравнений (7) или (9) обычно используется для вычисления волнового поля с подавленными ложными отражениями со стороны источника с поправкой на форму импульсов источника. Следует отметить, что представленные выше варианты осуществления изобретения включают определение функций взаимной корреляции операторов волн-спутников первого и одного или нескольких дополнительных источников. В общем случае, способ настоящего изобретения можно распространить на любое количество источников, каждый из которых работает на различной глубине в воде. При этом учитываются формы импульсов отдельных источников. На фиг. 3 А примеры форм импульсов двух источников, активируемых в различные моменты времени и на различных глубинах, показаны на кривой 40, где свободные от ложных отражений формы импульсов двух источников идентичны и, вследствие этого, не принимаются во внимание. На кривой 42,представленной на фиг. 3 В, показана форма импульса источника с подавленными ложными отражениями, вычисляемая с помощью разложения в соответствии с методикой Posthumus. Погрешности формы импульса источника отчетливо видны на кривой 42 как в более ранние, так и в более поздние моменты времени относительно сигнала с подавленными ложными отражениями. Временная разность связана с разностью времен активации источников, расположенных на различных глубинах, и разностью времен вступления, вызванной разностью глубин источников. Разложение по способу настоящего изобретения показано на кривой 44, фиг. 3 С. Видно, что при использовании способа настоящего изобретения погрешность формы импульса источника с подавленными ложными отражениями очень мала. На фиг. 4 А и 4 В показаны соответственно спектры волновых полей с подавленными ложными отражениями, представленных на фиг. 3, при использовании способа Posthumus и способа настоящего изобретения. На фиг. 5 А и 5 В показаны соответственно спектральная погрешность в волновых полях с подавленными ложными отражениями при использовании способа Posthumus и способа настоящего изобретения. На фиг. 6 А и 6 В показан соответственно спектр неотфильтрованного и точного знаменателя в уравнении (например, уравнении (7 согласно способу настоящего изобретения и спектр полосового фильтра, применяемого к двум последним членам в таком знаменателе, как объяснялось выше. На фиг. 7 А и 7 В показан соответственно спектр знаменателя в уравнениях, используемых в способеPosthumus, и спектр знаменателя в варианте с полосовым фильтром согласно способу настоящего изобретения. На фиг. 8 А и 8 В показан соответственно спектр двух последних членов знаменателя в уравнении согласно способу настоящего изобретения без полосовой фильтрации и с применением полосового фильтра. На фиг. 9 представлена блок-схема, иллюстрирующая вариант осуществления настоящего изобретения для определения волнового поля с подавленными ложными отражениями, генерируемого морским сейсмическим источником. В блоке 90 сейсмические данные получают по записям сигналов, активируя первый сейсмоисточник в первый момент времени и один или несколько дополнительных сейсмоисточников в характерные для них собственные моменты времени относительно времени индексации записей сигналов, при этом источники существенно связаны и расположены на различных глубинах в водоеме. В блоке 91 по записям сигналов в блоке 90 определяют первое волновое поле, которое возникло бы,если бы первый источник был активирован в момент времени, выбранный по отношению к времени начала выполнения записей сигналов, при этом первое волновое поле корректируют по времени с учетом глубины первого источника в воде. В одном из вариантов осуществления настоящего изобретения выбранным моментом времени может быть время начала записи. В блоке 92 по записям сигналов в блоке 90 определяют одно или несколько дополнительных волновых полей, которые возникли бы, если бы один или каждый из нескольких дополнительных источников были активированы в указанный момент времени, выбранный по отношению к времени начала выполнения записей сигналов, при этом одно или несколько дополнительных волновых полей корректируют по времени с учетом глубины одного или нескольких дополнительных источников в воде. В блоке 93 первое волновое поле из блока 91 объединяют с одним или несколькими дополнительными волновыми полями из блока 92 с целью определения волнового поля с подавленными ложными отражениями, соответствующего активации единственного сейсмоисточника. Способы настоящего изобретения могут обеспечить повышение амплитуды сигналов в морских сейсмических данных без необходимости увеличения размера используемых сейсмоисточников. В настоящем изобретении предлагается методика разложения записанных сейсмических сигналов на сигналы, которые были бы записаны с использованием только одного источника, в то время как два или несколько источников работают на различных глубинах и активируются с записью на одних и тех же сейсмограммах ОПВ. Волновое поле морских сейсмоисточников обладает существенной чувствительностью к гидростатическому давлению, которое, в свою очередь, является функцией глубины источника. Кроме того, может оказаться желательным использовать различные конфигурации источников для различных рабочих глубин. Таким образом, может потребоваться введение дополнительной поправки на характеристики формы импульса источника, как показано в уравнениях (7) и (9). Следует отметить, что такая поправка не была бы необходимой, если бы характеристики отдельных источников, за исключением волнспутников, были специально рассчитаны на то, чтобы быть близкими к идентичным. Существуют разнообразные известные методики проектирования, измерения или вычисления волновых полей сейсмоисточников, обладающие различной степенью точности. Волновое поле или выбранные позиции в волновом поле могут быть измерены непосредственно (например, измерение в дальней зоне поля), либо вычисление волнового поля может быть выполнено на основе физических моделей источника. Существуют также различные способы контроля источника, которые определяют волновое поле группы источников от выстрела к выстрелу с помощью различных датчиков, расположенных в группе сейсмоисточников. Сюда входит так называемый "Способ значимого источника", описанный в публикации Ziolkowski, A.,Parkes, G., Hatton, L. and Haugland, Т. [1982] The signature of an airgun array: Computation from near-fieldSource Monitoring Using Modeled Source Signatures with Calibration Function, патент США 7218572, выданный Parkes и находящийся в совместном владении с настоящим изобретением. Хотя изобретение описано с использованием ограниченного числа вариантов осуществления, специалисты, воспользовавшись раскрытым изобретением, смогут вывести из описания другие варианты осуществления, не отступающие от объема раскрытого изобретения. Соответственно объем настоящего изобретения ограничивается только прилагаемой формулой. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ определения волнового поля с подавленными ложными отражениями, генерируемого морским сейсмическим источником, по записям сигналов сейсмических данных, полученных путем активации первого сейсмоисточника в первый момент времени и одного или более дополнительных сейсмоисточников в различные моменты времени относительно времени индексации записей сигналов, при этом источники, по существу, связаны и расположены на различных глубинах в водоеме, включающий следующие шаги: определяют первое волновое поле, которое возникло бы, если бы первый источник был активирован в момент времени, относящийся к времени начала выполнения записей сигналов, при этом первое волновое поле определяют по записям сигналов, а указанное время корректируют с учетом глубины пер-9 021738 вого источника в воде; определяют одно или более дополнительных волновых полей, которые возникли бы, если бы один или каждый из дополнительных источников были активированы в указанный момент времени, относящийся к времени начала выполнения записей сигналов, при этом одно или более дополнительных волновых полей определяют по записям сигналов, а указанное время корректируют с учетом глубины одного или более дополнительных источников в воде; объединяют первое волновое поле с одним или более дополнительными волновыми полями с целью определения волнового поля с подавленными ложными отражениями, соответствующего активации единственного сейсмоисточника. 2. Способ по п.1, отличающийся тем, что первый источник и один или каждый из дополнительных источников содержат по меньшей мере по одной пневмопушке. 3. Способ по п.1, отличающийся тем, что включает внесение поправки на формы импульсов источников, буксируемых на различных глубинах. 4. Способ по п.1, отличающийся тем, что включает определение функций взаимной корреляции операторов волн-спутников первого и одного или более дополнительных источников. 5. Способ по п.4, отличающийся тем, что включает полосовую фильтрацию функций взаимной корреляции. 6. Способ по п.1, отличающийся тем, что первый источник и один или более дополнительных источников расположены в различных вертикальных плоскостях, чтобы избежать воздействия воздуха при определении отдельных волновых полей. 7. Способ по п.1, отличающийся тем, что включает проектирование источников, буксируемых на различных глубинах, таким образом, чтобы они имели, по существу, идентичные свободные от ложных отражений формы импульсов. 8. Способ по п.1, отличающийся тем, что указанным моментом времени является время начала записи. 9. Способ по п.3, отличающийся тем, что объединение первого волнового поля с одним или более дополнительными волновыми полями с целью определения волнового поля с подавленными ложными отражениями включает применение следующего уравнения:S1 и S2 - свободные от ложных отражений формы импульсов для первого и второго источников соответственно;N1 и N2 - представляющие помехи члены для первого и второго сейсмоисточников соответственно, где t1 и t2 - времена активации первого и второго сейсмоисточников соответственно;d1 и d2 - рабочие глубины первого и второго сейсмоисточников соответственно;- угол испускания сейсмической волны, распространяющейся от источника, по отношению к вертикали;- акустическая скорость в воде. 10. Способ по п.5, отличающийся тем, что объединение первого волнового поля с одним или более дополнительными волновыми полями с целью определения волнового поля с подавленными ложными отражениями с использованием полосовой фильтрации включает применение следующего уравнения:S1 и S2 - свободные от ложных отражений формы импульсов для первого и второго источников соответственно;N1 и N2 - представляющие помехи члены для первого и второго сейсмоисточников соответственно;F - полосовой фильтр, где t1 и t2 - времена активации первого и второго сейсмоисточников соответственно;d1 и d2 - рабочие глубины первого и второго сейсмоисточников соответственно;- угол испускания сейсмической волны, распространяющейся от источника, по отношению к вертикали;- акустическая скорость в воде. 11. Способ определения волнового поля с подавленными ложными отражениями, генерируемого морским сейсмическим источником в водоеме, включающий следующие шаги: активируют первый сейсмоисточник на первой глубине в воде в первый момент времени, относящийся к времени начала записи сейсмических сигналов; активируют один или более дополнительных сейсмоисточников, по существу, связанных с первым сейсмоисточником, при этом все дополнительные источники расположены в воде на разных глубинах и активируются в различные моменты времени относительно времени начала записи сейсмических сигналов, при этом сейсмическая волна от дополнительных источников, активированных в различные моменты времени, присутствует в записях сейсмических сигналов; записывают сейсмические сигналы, возникающие в результате активации первого источника и одного или более дополнительных источников; определяют первое волновое поле, которое возникло бы, если бы первый источник был активирован в момент времени, относящийся к времени начала выполнения записей сигналов, при этом первое волновое поле определяют по записям сигналов, а указанное время корректируют с учетом глубины первого источника в воде; определяют одно или более дополнительных волновых полей, которые возникли бы, если бы один или каждый из дополнительных источников были активированы в указанный момент времени, относящийся к времени начала выполнения записей сигналов, при этом каждое дополнительное волновое поле определяют по записям сигналов и корректируют согласно различному времени с учетом глубины соответствующего источника в воде; объединяют первое волновое поле с одним или более дополнительными волновыми полями с целью определения волнового поля с подавленными ложными отражениями, соответствующего активации единственного сейсмоисточника. 12. Способ по п.11, отличающийся тем, что первый источник и один или каждый из дополнительных источников содержат по меньшей мере по одной пневмопушке. 13. Способ по п.11, отличающийся тем, что первый источник и один или более дополнительных источников расположены в различных вертикальных плоскостях, чтобы избежать воздействия воздуха при определении отдельных волновых полей. 14. Способ по п.11, отличающийся тем, что включает внесение поправки на формы импульсов источников, буксируемых на различных глубинах. 15. Способ по п.11, отличающийся тем, что включает определение функций взаимной корреляции операторов волн-спутников первого и одного или более дополнительных источников. 16. Способ по п.11, отличающийся тем, что включает полосовую фильтрацию функций взаимной корреляции. 17. Способ по п.11, отличающийся тем, что включает проектирование источников, буксируемых на различных глубинах, таким образом, чтобы они имели, по существу, идентичные свободные от ложных отражений формы импульсов. 18. Способ по п.11, отличающийся тем, что указанным моментом времени является время начала записи. 19. Способ по п.14, отличающийся тем, что объединение первого волнового поля с одним или более дополнительными волновыми полями с целью определения волнового поля с подавленными ложными отражениями включает применение следующего уравнения:S1 и S2 - свободные от ложных отражений формы импульсов для первого и второго источников соответственно;N1 и N2 - представляющие помехи члены для первого и второго сейсмоисточников соответственно, где t1 и t2 - времена активации первого и второго сейсмоисточников соответственно;d1 и d2 - рабочие глубины первого и второго сейсмоисточников соответственно;- угол испускания сейсмической волны, распространяющейся от источника, по отношению к вертикали;- акустическая скорость в воде. 20. Способ по п.16, отличающийся тем, что объединение первого волнового поля с одним или более дополнительными волновыми полями с целью определения волнового поля с подавленными ложными отражениями с использованием полосовой фильтрации включает применение следующего уравнения:S1 и S2 - свободные от ложных отражений формы импульсов для первого и второго источников соответственно;N1 и N2 - представляющие помехи члены для первого и второго сейсмоисточников соответственно; где t1 и t2 - времена активации первого и второго сейсмоисточников соответственно;d1 и d2 - рабочие глубины первого и второго сейсмоисточников соответственно;- угол испускания сейсмической волны, распространяющейся от источника, по отношению к вертикали;- акустическая скорость в воде.