Способ идентификации и удаления многократных отражений для построения изображения с помощью пучков

СПОСОБ ИДЕНТИФИКАЦИИ И УДАЛЕНИЯ МНОГОКРАТНЫХ ОТРАЖЕНИЙ ДЛЯ ПОСТРОЕНИЯ ИЗОБРАЖЕНИЯ С ПОМОЩЬЮ ПУЧКОВ Настоящее изобретение заключает в себе использование модельноуправляемых и управляемых данными методик для подавления многократных волн в сейсмических данных, используемых в модели прогнозирования, которая включает в себя многократно отраженные и отраженные от поверхности сейсмические волны. Настоящее изобретение включает в себя методики пучков и свертку спрогнозированного многократного пучка с помощью сегмента смоделированного неполнократного пучка для получения свернутого многократного пучка. Свернутый многократный пучок может затем быть подвергнут обратной свертке для подавления многократных отражений,которые присутствуют в исходном входном пучке. 017379 Область техники Настоящее изобретение относится к сейсмической разведке и обработке, и, более конкретно, к построению изображений с помощью пучков, и способу прогнозирования многократных отражений на основании однократных отражений путем сочетания модельноуправляемой методики и методики, управляемой данными. Предшествующий уровень техники В нефтяной промышленности технологии сейсморазведки обычно используются для облегчения поиска и оценки подземных углеводородных месторождений. В сейсморазведке один или более источников сейсмической энергии излучает волны в подземный интересующий регион, такой как геологическая формация. Эти волны входят в формацию и могут быть рассеяны, например, отражены или преломлены подземными сейсмическими рефлекторами, например, границами между подземными формациями,имеющими различные упругие свойства. Отраженные сигналы могут быть отобраны или измерены одним или более приемниками, и полученные данные записаны. Записанные образцы могут называться сейсмическими данными или набором "сейсмических трасс". Сейсмические данные могут быть проанализированы для извлечения подробностей структуры и свойств исследуемого региона толщи пород. Сейсморазведка состоит из трех раздельных стадий: сбор данных, обработка данных и интерпретация данных. Успешность операции сейсмической разведки зависит от удовлетворительного завершения всех трех стадий. В основном целью сейсморазведки является картирование или изображение части подземной области толщи пород (формации) путем передачи энергии в землю и регистрации отражений или эха, которые возвращаются от нижележащих слоев породы. Передаваемая в формацию энергия является обычно звуковой энергией или энергией поперечной волны. Распространяющаяся вниз звуковая энергия может быть получена из различных источников, таких как взрывы или наземные сейсмические вибраторы или пневматические пушки в морских средах. Сейсморазведка обычно использует один или более источников энергии и обычно большое количество датчиков или детекторов. Датчики, которые могут быть использованы для обнаружения возвращающейся сейсмической энергии, являются обычно сейсмоприемниками (наземная разведка) или гидрофонами (морская разведка). Во время наземной сейсморазведки источник энергии может быть расположен в одном или более местоположениях близко к поверхности земли над интересующей геологической структурой или формацией, известных как точки взрыва. Каждый раз, когда источник приводится в действие, источник генерирует сейсмический сигнал, который распространяется вниз через толщу пород и, по меньшей мере, частично отражается от подземных неоднородностей различного типа, включающих в себя отражения от границ слоев породы. В основном частичное отражение сейсмических сигналов может происходить каждый раз, где имеется изменение упругих свойств подземных материалов. Отраженные сейсмические сигналы передаются обратно на поверхность земли, где они регистрируются в нескольких местоположениях в виде функции от времени распространения. Возвращающиеся сигналы оцифровываются и регистрируются как функция от времени (амплитуда относительно времени). Одной из распространенных проблем, связанных с сейсмической энергией, записанной приемниками во время этапа сбора данных, является то, что сейсмические трассы часто содержат и желательные сейсмические отражения ("однократные" отражения) и нежелательные многократные отражения, которые могут перекрывать или подавлять однократные сейсмические отражения. Однократное отражение является звуковой волной, которая проходит из источника в приемник с одним отражением от подземного сейсмического рефлектора. Многократное отражение является волной, которая отразилась по меньшей мере три раза (вверх, вниз и опять вверх) до того, как была обнаружена приемником. В зависимости от их времени задержки относительно основных событий, с которыми они связаны, многократные отражения обычно описываются как короткие, интерферирующие с их собственным однократным отражением, или как длинные, проявляющиеся как отдельные события. Существует также множество многократных событий, которые хорошо известны в данной области техники. Существуют сигналы, которые "попадают в ловушку" в водном слое между двумя жесткими рефлекторами, свободной поверхностью и дном водного слоя. Существуют "неполнократные события",которые являются отражениями, которые отличаются наличием дополнительного прохождения через водный слой сразу после испускания или непосредственно перед обнаружением. Существуют "остаточные" поверхностно-связанные события, где первое и последнее отражения вверх происходят ниже первого (водного) слоя, и имеется по меньшей мере одно отражение от свободной поверхности между ними. Также имеется "межслойное" многократное отражение, которое имеет отражение вниз, возникающее от подземного рефлектора. В большинстве случаев многократные отражения не содержат какой-либо полезной информации,которая была бы более полезна, чем извлеченная из однократных отражений. Более того, многократные отражения от морского дна известны как наиболее серьезные проблемы шума в обработке сейсмических данных во многих морских областях. Многократные отражения могут значительно маскировать однократные события отражения для структурного представления и засорять информацию об амплитуде относительно смещения (AVO). В связи с этим удаление многократных отражений или, по меньшей мере,-1 017379 подавление многократных отражений является необходимой частью этапа обработки сейсмических данных во многих средах, особенно в морских условиях, где многократные отражения особенно сильно влияют на однократные отражения. В случае глубоководных данных подавление многократных отражений и неполнократных отражений первого порядка и нескольких следующих порядков имеет большую важность. Эти довольно сильные многократные отражения могут иметь то же время распространения, что и однократные отражения от целевых рефлекторов. В предшествующем уровне техники существует несколько способов для подавления многократных отражений в зависимости от атрибутов используемых многократных отражений. Один класс способов ослабления многократных отражений представляет собой способы прогнозирования, в которых многократные отражения прогнозируются, исходя из соответствующих им однократных отражений. Методики прогнозируемого ослабления многократных отражений из предшествующего уровня техники могут быть в основном разделены на две категории: модельноуправляемые методики и методики, управляемые данными. Модельноуправляемые методики обычно используют модель толщи пород и зарегистрированные данные для прогнозирования или моделирования многократных отражений с использованием оцененной функции отражения от морского дна и рассчитанных функций амплитуд для моделирования многократных отражений в водном слое, при этом указанные спрогнозированные многократные отражения затем вычитаются из исходных данных. Другие модельноуправляемые методики используют модель толщи пород и модель отражающей среды для прогнозирования установившихся многократных отражений. Модельноуправляемые методики используют тот факт, что однократные отражения и многократные отражения являются физически зависимыми через сверточное отношение и прогнозируют многократные отражения путем кросс-свертки релевантных однократных отражений так, чтобы они содержали установившиеся составляющие для многократных отражений. Методики, управляемые данными, могут обычно поддерживать сложные геометрии и требуют мало или вообще не требуют использования информации о свойствах недр. Технологии, основанные на моделях, обычно являются менее затратными по сравнению с технологиями, управляемыми данными, однако последние обычно являются более гибкими. Некоторые модельноуправляемые методики требуют использования структурной информации, т.е. информации о подземной структуре, определение которой является основанием для выполнения сейсморазведки в первую очередь. Другие модельноуправляемые методики требуют использования формы импульса источника, которая не будет чистой дельта-функцией из-за реверберации и ограничения частотной полосы пропускания. Некоторые модельноуправляемые методики требуют использования информации о структуре и форме импульса, в то время как другие используют согласующий фильтр для учета искаженной формы импульса источника. Методики, управляемые данными, зависят от прогнозирования многократных отражений из компонентов однократных отражений. Фактически, эта методика использует существующие сейсмические данные для генерации многократных отражений, и эти сгенерированные отражения затем вычитаются из существующих данных. Одна такая методика предшествующего уровня техники является управляемой данными и известна как "подавление поверхностно-связанных многократных отражений" или "SRME". Если кратко, то этот способ работает, используя существующие данные для создания набора данных,который содержит только прогнозы многократных отражений, которые существуют в данных. В особенности способ стремится спрогнозировать сейсмическое выражение многократных отражений и после адаптации к существующим многократным отражениям в данных спрогнозированные многократные отражения вычитаются из исходных данных, оставляя (по меньшей мере, теоретически) только основную энергию.SRME методики, управляемые данными, являются привлекательными решениями для прогнозирования многократных отражений в сложных геологических условиях, они не требуют какого-либо априорного знания недр (отражательной способности, структур, скоростей). Однако эти способы требуют использования одного местоположения взрыва для каждого местоположения приемника, чего не требуется для большинства трехмерных ("3D") геометрий сбора данных. Применение SRME методик обычно вызвано сложными трехмерными многократными отражениями вследствие большого расстояния между взрывами, узкой зоной распространения и/или широкой кабельной сетью. Пропущенные данные могут быть интерполированы или экстраполированы из существующих данных, но интерполяция и экстраполяция имеют проблемы с наложением сейсмических данных, вызванным большим расстоянием между взрывами и/или приемниками. Усовершенствованные способы интерполяции и экстраполяции также могут иметь трудности в применении и быть дорогими. Общей причиной этих сложных трехмерных данных, которые бросают вызов трехмерным SRME способам, является наличие неровностей наверху солевого пласта. Но любой тип сложной перекрывающей толщи может вызывать сложные трехмерные сейсмические данные, которые трудно интерполировать. Другие методики, управляемые данными, используют прогнозируемое обращение свертки, которое является фильтрующим способом, который предполагает, что многократные отражения являются периодическими, в то время как однократные отражения таковыми не являются. Это предположение обычно встречается для данных из водных глубин менее 500 мс (приблизительно 1200 футов или 365,76 м) и-2 017379 приблизительно прослоенной подземной геологии. В областях с водными глубинами более 500 мс, где разность между скоростями однократных и многократных отражений является значительной, могут быть использованы способы фильтрации скоростей (в противоположность прогнозируемым способам), такие как tau-p и f-k фильтрация, где переменная f представляет частоту, k представляет волновое число, р представляет параметр луча и tau представляет интервал времени с нулевым смещением. Однако фильтрующие способы обычно требуют использования определения или, по меньшей мере,обоснованного предположения о кажущихся скоростях распространения волны в подземной среде, через которую отраженные сейсмические волны проходят на своем пути от сейсмического источника к приемнику. Эти скорости могут значительно различаться вследствие сочетания различий подземной структуры и свойств породы. В дополнение, прогнозируемое обращение свертки часто приводит к случайным повреждениям однократных отражений вследствие трудности разделения многократных отражений и однократных отражений. Более того, прогнозируемое обращение свертки часто не способно учесть фактор нелинейности при отражении, что в основном вызвано неполнократными отражениями. Один способ предшествующего уровня техники, который имеет расширенное прогнозируемое обращение свертки для глубоководных приложений, использует технику пучков. Этот способ применяет локальное наклонное суммирование (или другие ограничивающие по глубине способы) к данным для декомпозиции зарегистрированных волновых полей на компоненты пучков. Эти компоненты проходят путь приблизительно вдоль траектории луча. Простая трассировка в водном слое описывает длиннопериодические реверберации и связанные с ними однократные и многократные события, возникающие в компонентах пучков в волновом поле. На основании информации из трассировки временные ряды компонента пучка однократного отражения могут быть сдвинуты в соответствии с временами пробега трассировки и затем проанализированы с помощью многоканального фильтра прогнозирования. Спрогнозированные временные ряды считаются энергиями многократных отражений и удаляются из лучевых компонентов исходных данных после многоканальной согласованной фильтрации. Фиг. 1 иллюстрирует блок-схему для одного примера способа предшествующего уровня техники, в котором обращение свертки используется вместе с методикой пучков для подавления многократных отражений. Способ предшествующего уровня техники включает в себя использование модели 4 толщи пород, которая относится к интересующему геологическому региону, и инициирование набора 6 данных пучка, который был определен из сейсмических данных интересующего геологического региона. Способ предшествующего уровня техники включает в себя ряд циклов, в которых выбирают входной пучок 8,генерирующую многократные отражения поверхность 10 и временное окно 12. Выбирают один или более временных окон для гарантии того, что сигнал в пределах каждого селектора является установившимся. Затем проводят пробные лучи из точки 14 детектора и определяют 16 установившееся неполнократное отражение. Установившееся неполнократное отражение является неполнократным отражением,удовлетворяющим закону Снеллиуса для отражения от генерирующей многократные отражения поверхности. Получают 18 однократный пучок, соответствующий установившемуся неполнократному отражению, и однократный пучок трансформируется в спрогнозированный многократный пучок 20. Спрогнозированный многократный пучок затем подвергается обратной свертке для удаления многократных отражений, которые присутствуют во входном пучке 22. Несмотря на то что техника пучков имеет улучшенные методики подавления многократных отражений из предшествующего уровня техники, все еще имеется необходимость в улучшенном способе,который обеспечивает более точное прогнозирование многократных отражений и, следовательно, допускает более точное вычитание этих многократных отражений из данных. Предшествующие техники пучков предполагают, что существует единственная доминантная генерирующая многократные отражения поверхность 10, и спрогнозированные многократные пучки относятся только к этой генерирующей многократные отражения поверхности 10 и не содержат спрогнозированных многократных отражений от других генерирующих многократные отражения поверхностей. Настоящее изобретение улучшает техники пучков предшествующего уровня техники, включая в себя спрогнозированные многократные пучки от генерирующих многократные отражения поверхностей, которые явно не использовались для определения установившихся неполнократных отражений. Сущность изобретения Настоящее изобретение преодолевает описанные выше и другие недостатки предшествующего уровня техники путем обеспечения нового и улучшенного способа для прогнозирования многократных отражений на основе однократных отражений, которые объединяют признаки из модельноуправляемых методик и методик, управляемых данными. Это достигается путем определения модельноуправляемого установившегося прогноза на основании модели толщи пород и улучшения этого прогноза с помощью управляемого данными прогноза вокруг устойчивого прогноза. Следует понимать, что модельноуправляемый устойчивый прогноз может быть заменен априорным определением устойчивых прогнозов, таких как допущение слоистой модели. Один вариант осуществления настоящего изобретения включает в себя способ для генерации модели прогнозирования многократно отраженных поверхностно-связанных сейсмических волн, который включает в себя этапы, на которых инициализируют модель толщи пород, относящуюся к геологическо-3 017379 му объему, и выбирают набор данных пучка, полученный из сейсмических данных, относящихся к геологическому объему. Способ также включает в себя этап, на котором выбирают входной пучок из набора данных пучка, генерирующую многократные отражения поверхность из модели толщи пород и временное окно. Установившееся неполнократное отражение определяют с использованием входного пучка,генерирующей многократные отражения поверхности и временного окна. Затем получают однократный пучок, который соответствует установившемуся неполнократному отражению, и получают смоделированный неполнократный пучок, относящийся к однократному пучку. Смоделированный неполнократный пучок свертывается с однократным пучком для генерации свернутого многократного пучка. Свернутый многократный пучок сравнивается со входным пучком для удаления многократных отражений во входном пучке путем согласованной фильтрации. Свернутый многократный пучок используется для обеспечения более точного способа прогнозирования и удаления многократных отражений, по сравнению со способами предшествующего уровня техники. Этап свертки смоделированного неполнократного пучка с однократным пучком не включен в способы предшествующего уровня техники и позволяет настоящему изобретению более точно прогнозировать многократные отражения. Например, в одном варианте осуществления настоящего изобретения этап, на котором свертывают однократный пучок и смоделированный неполнократный пучок для получения свернутого многократного пучка, включает в себя этапы, на которых преобразую, однократный пучок в спрогнозированный многократный пучок путем смещения времени пробега, и сворачивают спрогнозированный многократный пучок с сегментом смоделированного неполнократного пучка для получения свернутого многократного пучка, при этом сегмент смоделированного неполнократного пучка начинается на генерирующей многократные отражения поверхности и заканчивается в месте расположения детектора."Свертка" является известной данной области техники. В целом, она является математической операцией над двумя функциями, которая представляет процесс линейной фильтрации. Свертка может быть применена к двум функциям от времени или пространства (или от других переменных) для получения третьей функции, являющейся результатом свертки. Несмотря на то что математическое определение является симметричным по отношению к двум входным функциям, в обработке сигналов является общеупотребительным говорить, что одна из функций является фильтрующей по отношению к другой функции. Ответ многих физических систем может быть представлен математически в виде свертки. Например, свертка обычно используется для моделирования фильтрации сейсмической энергии различными слоями породы в толще пород. Один вариант осуществления настоящего изобретения определяет установившееся неполнократное отражение с использованием лучей, которые проводятся из места расположения детектора, при этом место расположения детектора основывается на входном пучке, и установившееся неполнократное отражение относится к одному из лучей, проведенных из места расположения детектора. Специалист в данной области техники оценит, что фраза "временное окно" используется здесь для описания всего пучка или одного из сегментов, включенных в пучок. Следует также понимать, что использование методик пучков для обработки сейсмических данных является хорошо известным в данной области техники, и эти методики находятся в объеме настоящего изобретения. Пучки в основном определены как энергетические компоненты, которые частично локализованы в пространстве и глубине. Некоторыми примерами пучков являются гауссовы пучки и другие негауссовы пучки, такие как пучки со сложными лучами, управляемые пучки и пучки с конечноразностным решением некоторых версий волнового уравнения. В одном варианте осуществления настоящего изобретения сегмент смоделированного неполнократного пучка свертывается со спрогнозированным многократным пучком для генерации свернутого многократного пучка. Сегмент смоделированного неполнократного пучка начинается на генерирующей множественные отражения поверхности и включает в себя временное окно, которое задается явно или определяется другим горизонтом модели. В другом варианте осуществления настоящего изобретения свернутый многократный пучок может быть получен путем прямой свертки временных рядов спрогнозированных многократных отражений и спрогнозированных смоделированных неполнократных пучков. Следует также понимать, что настоящее изобретение предназначено для использования с системой,которая включает в себя в основном компьютерную конфигурацию, включающую в себя по меньшей мере один процессор, по меньшей мере одно запоминающее устройство для хранения программного кода или других данных, видеомонитор или другое устройство отображения (т.е. жидкокристаллический дисплей) и по меньшей мере одно устройство ввода. Процессор является предпочтительно платформой на основе микропроцессора или микроконтроллера, которая способна отображать изображения и выполнять сложные математические алгоритмы. Запоминающее устройство может включать в себя память с произвольным доступом (RAM) для хранения события или других данных, сгенерированных или используемых во время конкретного процесса, связанного с настоящим изобретением. Запоминающее устройство может также включать в себя память только для чтения, для хранения программного кода, для управления и процессов настоящего изобретения.-4 017379 Дополнительные признаки и преимущества настоящего изобретения описаны и станут понятны из следующего подробного описания изобретения и чертежей. Краткое описание чертежей Эти и другие предметы, признаки и преимущества настоящего изобретения станут более понятны вместе со следующим описанием, формулой изобретения и прилагающимися чертежами, на которых фиг. 1 иллюстрирует блок-схему способа предшествующего уровня техники для подавления многократных отражений в сейсмических данных; фиг. 2 - блок-схему одного варианта осуществления настоящего изобретения для подавления многократных отражений в сейсмических данных; фиг. 3 - один вариант осуществления настоящего изобретения, который включает в себя преобразование однократного пучка в спрогнозированный многократный пучок; фиг. 4 - один вариант осуществления настоящего изобретения, в котором лучи проводятся из места расположения детектора и отражаются от поверхности водного дна в морской среде; фиг. 5 - один вариант осуществления настоящего изобретения, в котором однократный пучок, соответствующий неполнократному, определяется с помощью трассировки лучей; фиг. 6 - способ, используемый одним вариантом осуществления настоящего изобретения, для определения установившегося неполнократного отражения; фиг. 7 - схематическое изображение указанного выше способа определения установившегося неполнократного отражения, используемого в одном варианте осуществления настоящего изобретения; фиг. 8 - один вариант осуществления настоящего изобретения, в котором однократный пучок сворачивается с неполнократным пучком; фиг. 9 - схематическое изображение однопучковой обратной свертки, которая используется одним вариантом осуществления настоящего изобретения; фиг. 10 - блок-схему другого варианта осуществления настоящего изобретения для подавления многократных отражений в сейсмических данных; фиг. 11 - блок-схему еще одного варианта осуществления настоящего изобретения для подавления многократных отражений в сейсмических данных; фиг. 12 - блок-схему одного варианта осуществления настоящего изобретения для подавления многократных отражений в сейсмических данных; фиг. 13 - блок-схему другого варианта осуществления настоящего изобретения для подавления многократных отражений в сейсмических данных и фиг. 14 - один вариант осуществления настоящего изобретения, в котором однократный пучок, соответствующий неполнократному, определяется без трассировки лучей. Подробное описание изобретения Несмотря на то что это изобретение допускает существование вариантов осуществления во многих различных формах, предпочтительные варианты осуществления показаны на чертежах и описаны здесь в подробностях с пониманием того, что настоящее раскрытие подразумевается как пример принципов изобретения и не предназначено для ограничения широкого аспекта изобретения проиллюстрированными вариантами осуществления. Один вариант 30 осуществления настоящего изобретения проиллюстрирован на фиг. 2. Вариант осуществления включает в себя использование модели толщи пород, которая соответствует интересующей геологической области 32, и выбор набора данных пучка, полученного из сейсмических данных,собранных в интересующей геологической области 34. Выбирают входной пучок из набора 36 данных пучка, генерирующая многократные отражения поверхность 38 из модели толщи пород и временное окно 40. Установившееся неполнократное отражение определяют с использованием входного пучка, генерирующей многократные отражения поверхности и временного окна 42. Получают 44 однократный пучок, соответствующий установившемуся неполнократному отражению. Получают 46 смоделированный неполнократный пучок, связанный с однократным пучком. Смоделированный неполнократный пучок свертывают с однократным пучком для определения свернутого многократного пучка 48. Выполняют обратную свертку многократного пучка со входным пучком для удаления многократных отражений во входном пучке 50. В другом варианте осуществления настоящего изобретения спрогнозированный однократный пучок может быть непосредственно свернут со смоделированным неполнократным пучком. В еще одном варианте осуществления настоящего изобретения лучи проводят из места расположения детектора, и установившееся неполнократное отражение определяют на основании одного из лучей. Место расположения детектора в этом варианте осуществления основано на входном пучке. В дополнительном варианте осуществления настоящего изобретения сегмент смоделированного неполнократного пучка свертывают с однократным пучком для получения свернутого многократного пучка. Сегмент смоделированного неполнократного пучка начинается на генерирующей многократные отражения поверхности и включает в себя временной интервал, который задается явно или определяется другим горизонтом модели. Свернутый многократный пучок затем подвергают обратной свертке со входным пучком для удаления многократных отражений во входном пучке.-5 017379 Как описано выше, способы предшествующего уровня техники использовали способы наклонного суммирования или другие ограничивающие по глубине способы для сейсмической трассировки для разделения записанного волнового поля на компоненты пучка, и эти способы известны в данной области техники. Настоящее изобретение использует локальное наклонное суммирование для разделения записанного волнового поля на компоненты, которые локализуются и по положению, и по глубине. Эти компоненты являются теми, которые будут зарегистрированы в центре вступления 54, 56 пучка в местоположениях А 58 и В 60, как проиллюстрировано на фиг. 3. Энергия 54 пучка, который достигает местоположения В 60, отражается в слое 64 воды, и предполагается, что он поступает в виде многократных отражений 62 в пучке 56, зарегистрированном в местоположении А 58. Сдвиг пучка В 54 на время прохождения TAB 66 трассировки от местоположения В 60 до местоположения А 58 выстраивает события в пучке В 54 с многократными отражениями 68, 70 в пучке А 56. Как только события были выровнены и многократные отражения 68, 70 были идентифицированы, многократные отражения 68, 70 могут быть удалены. Когда описанное выше локальное наклонное суммирование используют для трехмерного сбора данных, зарегистрированная энергия не может быть полностью направлена в пучки, потому что волновое поле не является плотно измеренным вдоль всех регистрируемых направлений, так что возникает проблема потери данных или сильного наложения, которая требует точного определения времени распространения TAB трассировки. Например, как могут оценить специалисты в данной области техники, локальное наклонное суммирование выполняют в области общего смещения, а не в области общей глубинной точки. В одном варианте осуществления настоящего изобретения сделано допущение о том, что суммируемые скорости описывают глубину основных событий в сейсмограммах общей глубинной точки("ОГТ"). В основном суммируемые скорости являются приемлемым описанием однократных отражений от геологических структур над подземными солевыми формациями, которые могут отражаться в морские среды, становясь сильнейшими многократными отражениями. Фиг. 4-6 иллюстрируют способ, которым этот вариант осуществления настоящего изобретения вычисляет время распространения трассировки для определения установившегося неполнократного отражения. Входной пучок выбирают с местоположением "S" 74 источника и местоположением "D" 76 детектора, как проиллюстрировано на фиг. 4. Определяют генерирующую многократные отражения поверхность, в этом варианте осуществления это водное дно (wb") 78. Из детектора D 76 проводят лучи 80,которые отражаются обратно от генерирующей многократные отражения поверхности wb 78 и достигают свободной поверхности 82. Угловой интервал для лучей 80 является заранее определенным как dpx в направлении х (вертикально вниз от свободной поверхности 82) и dpy в направлении у (перпендикулярно листу). Отдельный луч 84 выбирают с заданными параметрами луча или глубины pd траектории луча от детектора, и точку вступления Q 86 луча 84 определяют на свободной поверхности 82, как проиллюстрировано на фиг. 5. Направление вступления (параметр луча pq' 88) луча 84 определяют в точке вступленияQ 86. Определяют отражение в точке Q 86 и рассчитывают параметр исходящего луча pq 90. Точки S 74 иQ 86 и их параметры луча ps 92 и pq 90 определяют однократное отражение, соответствующее неполнократному отражению с точками D 76 и Q 86 и параметрами луча pd 84 и pq' 88. Нахождение установившихся траекторий, которые описывают отражение в точке Q 86, требует использования траекторий. Поиск выполняют для нахождения отдельных отражений, которые появляются в различных точках S 74, Q 86 и D 76. В этом конкретном варианте осуществления поиск сравнивает трассированный ph с рассчитанным ph (ph является смещением глубины). Трассированный ph рассчитывается как где pm (среднеглубинная точка) в заданной точке соответствует конкретному пучку. Рассчитанныйph получают из уравнения нормального приращения ("NMO") где t является временем,h является половиной смещения,V является скоростями NMO и Т 0 является временем распространения с нулевым смещением.V и Т 0 выводят из суммируемых скоростей, полученных из модели толщи пород. Фиг. 6 иллюстрирует определение ph, где ph является локальным наклоном 96 кривой 98 на h-t графике 94. Пучок с ближайшим совпадением между трассированным ph и рассчитанным ph выбирают как установившееся неполнократное отражение. Фиг. 7 иллюстрирует схему способа, в котором установившееся неполнократное отражение определяют этим вариантом осуществления настоящего изобретения. Выбирают 102 глубинную точку pm, и-6 017379 пробное направление Pd выбирают из точки D детектора, которая соответствует выбранной pm 104. Параметр Pd используют для выполнения трассировки из точки D детектора до точки Q, выбранной между точкой S источника и точкой D детектора 106. Результатом трассировки является трассированное смещение глубины ph, которое сравнивают с рассчитанным смещением глубины ph 108, и если трассированный ph близко совпадает с рассчитанным ph, то этот процесс завершается, и пучок, соответствующий pm,выбирают в качестве установившегося неполнократного отражения 110. Если трассированный ph не совпадает с рассчитанным ph, то выбирают другое пробное направление pd, и процесс выполняют опять 112 до тех пор, пока не будет получено удовлетворительное совпадение между трассированным Ph и рассчитанным Ph. Этот вариант осуществления настоящего изобретения также включает в себя свертывание спрогнозированного многократного пучка с сегментом смоделированного пучка для получения свернутого многократного пучка. Как проиллюстрировано на фиг. 8, спрогнозированный многократный пучок В 114 свертывают с сегментом смоделированного неполнократного пучка С 116, начинающегося на генерирующей многократные отражения поверхности 118. Результатом свертки является пучок Е 120, который является свернутым многократным пучком, который подвергается обратной свертке со входным пучком А 122 для удаления многократных отражений, которые присутствуют во входном пучке А 122. Схематическая иллюстрация этого варианта осуществления представлена на фиг. 9, где пучок В 114 используют как прогноз 124 со стороны источника, и пучок С 116 используют как прогноз 126 со стороны детектора. Этот вариант осуществления использует фильтр Винера 128 и входы из пучков А 122, В 114 и С 116 для генерации оценки существования многократных отражений в пучке А 122. Фильтр Винера известен в данной области техники. В основном он является причинным фильтром, который преобразует вход в требуемый выход как можно ближе, подвергая определенным ограничениям. Специалисты в данной области техники оценят, что имеются другие фильтры или средство, которое может выполнять эту конкретную функцию, и они предназначены быть в объеме настоящего изобретения. Как только многократные отражения в пучке А 122 были определены, многократные отражения затем удаляют 130 из пучка А 122. Один вариант осуществления настоящего изобретения проиллюстрирован на фиг. 10, где инициализируется 134 модель толщи пород, которая связана с конкретным интересующим геологическим регионом. Набор данных 136 пучка, который определяли из сейсмических данных интересующего геологического региона, также инициализируют. Этот вариант осуществления настоящего изобретения включает в себя ряд циклов, в которых выбирают входной пучок 138, генерирующую многократные отражения поверхность 140 и временное окно 142. Лучи проводят из точки детектора на основании входного пучка 144, и установившееся неполнократное отражение выбирают из одного из лучей 146. Получают 148 однократный пучок, соответствующий неполнократному отражению, и однократный пучок преобразуют в спрогнозированный многократный пучок 150 путем сдвига, соответствующего времени распространения луча, соответствующего установившемуся неполнократному отражению. Затем генерируют 152 смоделированный неполнократный пучок, относящийся к спрогнозированному многократному пучку. Сегмент смоделированного неполнократного пучка, который начинается на генерирующей многократные отражения поверхности, свертывают со спрогнозированным многократным пучком для получения свернутого многократного пучка 154. Свернутый многократный пучок затем или аккумулируют, или подвергают обратной свертке со входным пучком для удаления многократных отражений из входного пучка 162. Аккумулированные пучки могут быть использованы для реконструкции многократного прогнозирования как сейсмического следа или могут быть использованы для выполнения обратной свертки со входным пуском позже. Этот вариант осуществления настоящего изобретения позволяется для нескольких различных точек после свертки или аккумулировать свернутый многократный пучок, или подвергать обратной свертке свернутый многократный пучок для удаления многократных отражений из входного пучка 162. В зависимости от обрабатываемых данных эти этапы могут возникать до конца цикла для выбора временного окна 142-156 или непосредственно после этого цикла 142-156. Эти этапы также могут возникать после цикла для выбора генерирующей многократные отражения поверхности 140-158 или после цикла для выбора входного луча 138-160. Другой вариант осуществления настоящего изобретения проиллюстрирован на фиг. 11, и этот вариант осуществления включает в себя выбор неполнократного отражения, который находится в узком диапазоне около выбранного установившегося неполнократного отражения 176. Получают 178 однократный пучок, который соответствует выбранному неполнократному отражению, и этот однократный пучок преобразуют в спрогнозированный многократный пучок 180. Получают 182 смоделированный неполнократный пучок, который относится к спрогнозированному многократному пучку. Спрогнозированный многократный пучок затем свертывают с сегментом смоделированного неполнократного пучка для получения свернутого многократного пучка. Сегмент смоделированного неполнократного пучка, используемый этим этапом, начинается на генерирующей многократные отражения поверхности и заканчивается в точке детектора. В цикле для выбора неполнократного отражения 176-186 складывают прогнозы или свернутый многократный пучок для узкого диапазона около установившегося неполнократного отражения.-7 017379 Этапы аккумулирования свернутого многократного пучка или обратной свертки свернутого многократного пучка для удаления многократных отражений во входном пучке 194 могут возникать в этом варианте осуществления после цикла для выбора неполнократного отражения 176-186, цикла для выбора временного окна 172-188, цикла для выбора генерирующей многократные отражения поверхности 170-190 или цикла для выбора входного пучка 168-192. Другой вариант осуществления настоящего изобретения проиллюстрирован на фиг. 12, где выбранное 208 неполнократное отражение не связано с установившимся неполнократным отражением. Это неполнократное отражение используют для получения соответствующего однократного пучка 210. Однократный пучок преобразуют в спрогнозированный многократный пучок 212 и получают 214 смоделированный неполнократный пучок. Сегмент смоделированного неполнократного пучка, который начинается на генерирующей многократные отражения поверхности, сворачивают со спрогнозированным многократным пучком для получения свернутого многократного пучка 216. Свернутые многократные пучки,сгенерированные в цикле для неполнократного отражения 208-218, аккумулируют. В этом варианте осуществления этапы аккумулирования свернутого неполнократного пучка или обратной свертки свернутого многократного пучка для удаления многократных отражений во входном пучке 226 могут возникать после цикла для выбора неполнократного отражения 208-218, цикла для выбора временного окна 204220, цикла для выбора генерирующей многократные отражения поверхности 202-222 или цикла для выбора входного пучка 200-224. И установившееся неполнократное отражение, и неполнократное отражение могут быть определены для различных многократных отражений. Один вариант осуществления настоящего изобретения определяет установившееся неполнократное отражение или неполнократное отражение для многократного отражения со стороны источника. Другой вариант осуществления определяет установившееся неполнократное отражение или неполнократное отражение для многократного отражения со стороны детектора. Дополнительный вариант осуществления определяет установившееся неполнократное отражение или неполнократное отражение для многократного отражения и со стороны источника, и со стороны детектора. Специалист в данной области техники должен оценить, что могут быть ситуации, в которых модель толщи пород не может быть легко доступной, в таких случаях настоящее изобретение все еще способно прогнозировать и подавлять многократные отражения. Один вариант осуществления настоящего изобретения не включает в себя использование модели толщи пород. В этом варианте осуществления входной пучок выбирают 230 из открытого набора данных 228 пучка, как проиллюстрировано на фиг. 13. Выбирают временное окно 232, и выбор 236 неполнократного отражения облегчается априорным определением установившихся прогнозов 234, таких как допущение слоистой модели толщи пород или определение локального покрытия точек источника и детектора. Выбирают 236 неполнократное отражение и получают 238 соответствующий неполнократному отражению однократный пучок. Вычисляют 240 смоделированный неполнократный пучок, соответствующий однократному пучку, и однократный пучок свертывают со смоделированным неполнократным пучком 242. Описанные выше этапы могут быть повторены 230-248, 232-246, 336-244 с аккумулированным свернутым многократным пучком, или свернутый многократный пучок подвергают обратной свертке со входным пучком для удаления многократных отражений из входного пучка 250. Описанный выше вариант осуществления свертывает множество однократных и неполнократных пучков для диапазона точек Q 252 и диапазона значений параметров лучей Pq 254 и Ps 256, как проиллюстрировано на фиг. 14. Диапазон Q 252 и диапазоны Pq 254 и Ps 256 (или Pd 258) определяют из анализа входных данных. Диапазоны определяют так, чтобы достаточное количество пучков было включено для содержания установившейся составляющей для смоделированных многократных отражений, которое усиливается после суммирования установившейся составляющей с другими не установившимися составляющими. На фиг. 14 выбирают пучок с заданными параметром луча Pm=Ps 280+Pd 258, точкой Q 252 и параметрами луча Pq 254 и Ps 256 (или Pd 258). Определяют однократный пучок Pm1=Ps 256+Pq 254 и неполнократный пучок Pm2=Pd 258+Pq 254. В этом варианте осуществления сделано допущение, что отражение от поверхности в Q 252 находится в приделах от Pq 254 до Pq 260. Однократный пучок затем может быть свернут с неполнократным пучком. В этом варианте осуществления однократный пучок и неполнократный пучок сворачивают для диапазона значений точек Q 252, начинающегося в Q 252=(S 262+264)/2, или в каком либо другом заранее определенном месте. В дополнение, однократный пучок и неполнократный пучок свертывают для диапазона значений параметра Pq 254. Как описано выше, варианты осуществления настоящего изобретения, изображенные на фиг. 10-13,включают в себя циклы, которые иллюстрируют, что конкретные этапы этих вариантов осуществления могут быть повторены в зависимости от обрабатываемых данных. Конкретные варианты осуществления настоящего изобретения, описанные выше, включают в себя проведение лучей от точки детектора для выбора установившегося неполнократного отражения или неполнократного отражения, следует понимать, что есть другое средство для выбора установившегося неполнократного отражения или неполнократного отражения, например случайный выбор, и это средство считается находящимся в объеме настоящего изобретения.-8 017379 Несмотря на то что предшествующей спецификацией это изобретение было описано в отношении конкретных предпочтительных вариантов его осуществления и было изложено много подробностей в иллюстративных целях, специалистам в данной области техники следует понимать, что это изобретение может подвергаться изменению и что конкретные другие изложенные здесь подробности могут значительно меняться без выхода за основные принципы изобретения. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Реализуемый с помощью компьютера способ подавления многократных отражений в многократно отраженных поверхностно-связанных сейсмических волнах в подземной толще пород с использованием компьютерной системы, содержащей по меньшей мере один процессор, сконфигурированный для связи по меньшей мере с одним запоминающим устройством и сконфигурированный для исполнения компьютерного программного кода для выполнения способа, содержащего этапы, на которых:(a) инициализируют модель толщи пород с помощью по меньшей мере одного процессора компьютерной системы, относящуюся к геологическому объему в подземной толще пород;(b) задают набор данных пучка, сохраненный по меньшей мере на одном запоминающем устройстве, полученный из сейсмических данных, относящихся к геологическому объему;(c) задают входной пучок из набора данных пучка;(d) задают генерирующую многократные отражения поверхность из модели толщи пород;(f) обеспечивают выполнение компьютерного программного кода по меньшей мере на одном процессоре компьютерной системы, сконфигурированного для определения установившегося неполнократного отражения с использованием входного пучка, генерирующей многократные отражения поверхности и временного окна;(g) обеспечивают выполнение компьютерного программного кода по меньшей мере на одном процессоре компьютерной системы, сконфигурированного для получения однократного пучка, который соответствует установившемуся неполнократному отражению;(h) обеспечивают выполнение компьютерного программного кода по меньшей мере на одном процессоре компьютерной системы, сконфигурированного для получения смоделированного неполнократного пучка, который относится к однократному пучку;(i) обеспечивают выполнение компьютерного программного кода по меньшей мере на одном процессоре компьютерной системы, сконфигурированного для свертывания однократного пучка и смоделированного неполнократного пучка для получения свернутого многократного пучка; и(j) обеспечивают выполнение компьютерного программного кода по меньшей мере на одном процессоре компьютерной системы, сконфигурированного для осуществления обратной свертки свернутого многократного пучка с входным пучком для удаления многократных отражений во входном пучке;k) этапы (c)-(j) повторяют для множества пучков, этапы (d)-(j) повторяют для множества генерирующих многократные отражения поверхностей и этапы (e)-(j) повторяют для множества временных окон, для формирования изображения геологического объема. 2. Способ по п.1, в котором свертка однократного пучка и смоделированного неполнократного пучка для получения свернутого многократного пучка включает в себя этапы, на которых преобразуют однократный пучок в спрогнозированный многократный пучок и свертывают спрогнозированный многократный пучок с сегментом смоделированного неполнократного пучка для получения свернутого многократного пучка, при этом сегмент смоделированного неполнократного пучка начинается на генерирующей многократные отражения поверхности и заканчивается на более глубокой поверхности или в конце заранее определенного временного интервала. 3. Способ по п.1, в котором свернутый многократный пучок получают как прямую свертку однократного пучка со смоделированным неполнократным пучком. 4. Реализуемый с помощью компьютера способ подавления многократных отражений в многократно отраженных поверхностно-связанных сейсмических волнах в подземной толще пород с использованием компьютерной системы, содержащей по меньшей мере один процессор, сконфигурированный для связи по меньшей мере с одним запоминающим устройством и сконфигурированный для исполнения компьютерного программного кода для выполнения способа, содержащего этапы, на которых:(а) инициализируют модель толщи пород с помощью по меньшей мере одного процессора компьютерной системы, относящуюся к геологическому объему в подземной толще пород;(b) задают набор данных пучка, сохраненный по меньшей мере на одном запоминающем устройстве, полученный из сейсмических данных, относящихся к геологическому объему;(c) задают входной пучок из набора данных пучка;(d) задают генерирующую многократные отражения поверхность из модели толщи пород;(g) обеспечивают выполнение компьютерного программного кода по меньшей мере на одном про-9 017379 цессоре компьютерной системы, сконфигурированного для получения однократного пучка, который соответствует неполнократному отражению;(h) обеспечивают выполнение компьютерного программного кода на по меньшей мере одном процессоре компьютерной системы, сконфигурированного для получения смоделированного неполнократного пучка, который относится к однократному пучку;(i) обеспечивают выполнение компьютерного программного кода по меньшей мере на одном процессоре компьютерной системы, сконфигурированного для свертывания однократного пучка и неполнократного пучка для получения свернутого многократного пучка; и(j) обеспечивают выполнение компьютерного программного кода по меньшей мере на одном процессоре компьютерной системы, сконфигурированного для осуществления обратной свертки многократного пучка со входным пучком для удаления многократных отражений во входном пучке, и этапы (c)-(j) повторяют для множества пучков, этапы (d)-(j) повторяют для множества генерирующих многократные отражения поверхностей и этапы (е)-(j) повторяют для множества временных окон. 5. Способ по п.4, в котором определяют установившееся многократное отражение с использованием входного пучка генерирующей многократные отражения поверхности и временного окна, и заданное неполнократное отражение относится к установившемуся многократному отражению. 6. Способ по п.5, в котором неполнократное отражение задают путем задания диапазона значений параметра луча около значения параметра луча установившегося неполнократного отражения. 7. Способ по п.5, в котором неполнократное отражение задают путем задания диапазона местоположений около местоположения вступления установившегося неполнократного отражения. 8. Способ подавления многократных отражений в многократно отраженных поверхностносвязанных сейсмических волнах, содержащий этапы, на которых:(a) задают входной пучок из набора данных пучка;(d) получают однократный пучок, который соответствует неполнократному отражению;(f) свертывают однократный пучок и неполнократный пучок для получения свернутого многократного пучка и(g) осуществляют обратную свертку свернутого многократного пучка с входным пучком для удаления многократных отражений во входном пучке. 9. Способ по п.8, в котором этапы (с)-(g) повторяют для множества пучков, этапы (d)-(g) повторяют для множества генерирующих многократные отражения поверхностей и этапы (е)-(g) повторяют для множества временных окон. 10. Способ по любому из пп.1, 4 или 8, в котором аккумулируют свернутые многократные пучки до осуществления обратной свертки входного пучка со свернутым многократным пучком для удаления многократных отражений из входного пучка. 11. Способ по любому из пп.1, 4 или 8, в котором этап задания временного окна включает в себя этап, на котором задают весь входной пучок. 12. Способ по любому из пп.1, 4 или 8, в котором этап задания временного окна включает в себя этап, на котором задают часть входного пучка. 13. Способ по любому из пп.4 или 8, в котором неполнократное отражение задают для многократного отражения со стороны источника. 14. Способ по любому из пп.4 или 8, в котором неполнократное отражение задают для многократного отражения со стороны детектора. 15. Способ по любому из пп.4 или 8, в котором неполнократное отражение задают для многократного отражения как со стороны источника, так и со стороны детектора. 16. Способ по любому из пп.1 или 5, в котором установившееся неполнократное отражение определяют для многократного отражения со стороны источника. 17. Способ по любому из пп.1 или 5, в котором установившееся неполнократное отражение определяют для многократного отражения со стороны детектора. 18. Способ по любому из пп.1 или 5, в котором установившееся неполнократное отражение определяют для многократного отражения как со стороны источника, так и со стороны детектора. 19. Способ по любому из пп.1, 4 или 8, в котором пучки определяют как энергетические компоненты, которые частично локализованы в пространстве и глубине. 20. Способ по любому из пп.1, 4 или 8, в котором набор данных пучка включает в себя гауссов пучок. 21. Способ по любому из пп.1, 4 или 8, в котором набор данных пучка включает в себя пучки со сложными лучами. 22. Способ по любому из пп.1, 4 или 8, в котором набор данных пучка включает в себя управляемые пучки. 23. Способ по любому из пп.1, 4 или 8, в котором набор данных пучка включает в себя конечноразностные пучки.