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

1 Область техники Настоящее изобретение относится к способу обработки сейсмических данных. Предшествующий уровень техники Сейсмические данные собирают с целью анализа подповерхностной области земли, в частности, для разведки углеводородов. Сейсмические данные для анализа подповерхностных структур можно собирать на земле или на воде с использованием морских судов. Для получения данных предусмотрен сейсмический источник, который может представлять собой взрывчатые вещества (на земле) или импульс сжатого воздуха или подводные пушки (на море). Сигналы сейсмических данных, отраженные различными слоями, лежащими под поверхностью земли, известных как следы, регистрируют большим количеством, обычно, сотнями датчиков, например геофонов на земле или гидрофонов на море. Отраженные сигналы регистрируют и результаты анализируют для получения наличия слоевых образований, лежащих под поверхностью. Подобные признаки можно затем использовать для получения вероятности обнаружения залежей углеводородов. Анализ результатов с целью извлечения признаков наличия слоистых пород не является,однако, простым. В частности, там, где имеет место смена материала подповерхностных областей земли в горизонтальном направлении,может существовать несколько путей прохождения сигнала между сейсмическим источником и точкой в подповерхностной области, от которой отражается сигнал. Обычно то же касается обратного пути между точкой отражения и соответствующим сейсмическим датчиком, например, геофоном или гидрофоном. В случае наличия трех различных путей в каждом направлении, сигнал может распространяться от сейсмического источника к сейсмическому датчику через одну точку отражения по девяти различным круговым маршрутам. Поскольку анализ, производимый с учетом всех возможных путей требует больших затрат, необходимы некоторые средства упрощения обработки. В публикации Функции Грина для 3 мерной миграции по глубине перед накоплением, представленном на 57 конференции EAGE(ЕАГР) Европейская ассоциация геофизиковразведчиков и технической выставке, проходившей в Глазго, Шотландия, с 29 мая по 2 июня 1995 г., обсуждаются следующие две методики упрощения анализа, известные из уровня техники. Ясно, что сигналы, приходящие от одной точки отражения, в общем случае, достигают геофона за различное время и с различными амплитудами в зависимости от длины пройденного пути и характеристик распространения звука в подповерхностных слоях, через которые проходят звуковые волны. Следовательно, существует некоторое количество траекторий луча через подповерхностную область земли, 001707 2 которые относятся к сигналу, отраженному от одной точки отражения. Один из предложенных путей упрощения анализа в случае множественных траекторий луча состоит в выборе так называемого первого приходящего сигнала. Этот приходящий сигнал соответствует наиболее быстро распространяющемуся сейсмическому сигналу. Однако эта методика страдает тем недостатком, что первый приходящий сигнал редко является самым сильным сигналом и часто несет слишком низкую энергию для обеспечения надежного и точного анализа. Тем не менее,способы вычисления времени прохождения первого приходящего сигнала, зачастую оказываются дешевле и проще других способов. Некоторые из этих способов обычно (хотя и не вполне точно) именуют способами кратчайшей траектории, хотя в действительности они предусматривают вычисление траектории,отвечающей наименьшему времени прохождения, а не траектории луча, отвечающей наименьшей физической длине. Статьи, описывающие способ этого типа можно также обнаружить в журнале Геофизика, т. 56,1, январь 1991 г., Т.Дж. Мозер (T.J. Moser) Расчет сейсмических лучей методом кратчайшей траектории (Shortest path calculation of seismic rays),стр. 59-67; Геофизика, т. 58,7, июль 1993 г.,Роберт Фишер (Robert Fisher) и др. Определение кратчайшей траектории луча с помощью разреженных графов (Shortest path ray tracingwith sparse graphs), стр. 987-996; и Геофизика, т. 59,7, июль 1994 г., Т. Дж. Мозер, Миграция с использованием способа кратчайшей траектории (Migration using the shortest-path method),стр. 1110-1120. Ссылки на процесс расчета траектории луча, отвечающей наименьшему времени прохождения (а не траектории луча, отвечающей наименьшей длине луча), можно найти в этих статьях, соответственно, на стр. 59, реферат, строка 4; стр. 987, столбец 2, строка 20; и стр. 1111, столбец 2, строка 37. Другая методика, известная из уровня техники, предусматривает выбор приходящего сигнала, имеющего максимальную амплитуду. Однако выбор этого приходящего сигнала не всегда является простым, поскольку модель подповерхностного слоя является, в общем случае,всего лишь приближением. Способ приходящего сигнала с максимальной амплитудой обеспечивает наилучший приходящий сигнал лишь до тех пор, пока верны оценки амплитуды. Другая трудность способа приходящего сигнала с максимальной амплитудой состоит в том, что алгоритм выбора приходящего сигнала может быстро переключаться между ветвями в обоих направлениях. Однако в первой вышеупомянутой ссылке из уровня техники рассмотрены преимущества использования приходящего сигнала с максимальной амплитудой перед использованием первого приходящего сигнала. 3 Цель настоящего изобретения состоит в том, чтобы предоставить способ обработки сейсмических данных, который исправляет недостатки этих методик, отвечающих уровню техники. Согласно настоящему изобретению предусмотрен способ обработки сейсмических данных с использованием модели распространения сейсмической энергии в подповерхностной области, включающий в себя задание местоположения сейсмического источника, сейсмического приемника и точки отражения в модели распространения; идентификация альтернативных траекторий луча, которые начинаются в месте расположения сейсмического источника, отражаются в месте расположения точки отражения и оканчиваются в месте расположения сейсмического приемника; выбор траектории луча, отвечающей наименьшей длине луча; и применение выбранной траектории луча в последующей сейсмической обработке. Как было отмечено выше, применение критерия кратчайшей длины лучевого пути,практически, не дает результата при низкой амплитуде сигнала. Хотя физическая длина луча изменяется в результате изменения модели скоростей, физическая длина луча менее чувствительна к таким изменениям, чем оценка амплитуды. Причина в том, что амплитуда относится к кривизне лучей. Таким образом, малые погрешности в определении траектории луча могут обуславливать большую погрешность по амплитуде. Напротив, длина луча является интегральной величиной (интегралом длины дуги вдоль луча), так что она относительно нечувствительна к малым возмущениям траектории луча. Другими словами, алгоритм выбора отдельного луча, основанный на методике настоящего изобретения, скорее всего, будет переключаться не так быстро, как алгоритм, основанный на методике максимальной амплитуды. Это обеспечивает значительно более надежные исходные посылки для последующей обработки,чем какая-либо методика, известная из уровня техники. Дополнительные преимущественные признаки настоящего изобретения изложены в прилагаемых зависимых пунктах формулы изобретения. Краткое описание чертежей Теперь перейдем к иллюстративному описанию настоящего изобретения со ссылкой на прилагаемые чертежи: где на фиг. 1 показан сейсмический источник и совокупность сейсмических приемников, размещенных над несколькими слоями, находящимися в подповерхностной области земли; на фиг. 2 - местоположение одного источника и местоположение одного изображения, а также несколько траекторий луча между ними; на фиг. 3 - двухмерное сейсмическое изображение в координатах расстояние - глубина 4 для отдельной сейсмической модели, использующей сейсмический сигнал, приходящий первым (т.е. траектории, отвечающей наименьшему времени прохождения); на фиг. 4 - сейсмическое изображение в координатах расстояние - глубина для сейсмической модели, использующей сейсмический сигнал, приходящий с максимальной амплитудой; и на фиг. 5 - сейсмическое изображение в координатах расстояние - глубина для сейсмической модели, использующей приходящий сигнал, связанный с траекторией луча, имеющей наименьшую длину. Описание предпочтительного варианта реализации изобретения Согласно фиг. 1 сейсмический источник S,например взрывчатка, размещается на поверхности земли совместно с совокупностью геофонов R1-R5. Обычно сотни геофонов размещают по поверхности земли в виде двухмерной или трехмерной матрицы. Для простоты показано только пять геофонов R1-R5. Пути распространения сейсмических сигналов (траектории лучей) между источником S и геофонами R1-R5,проходящие через подповерхностные слои L1,L2 и L3, изображены пунктирными линиями. Согласно фиг. 1 лучи отражаются от нескольких горизонтов H1, H2 и Н 3. Для ясности показаны только некоторые траектории лучей. Поскольку сейсмические сигналы проходят в слоях с различными скоростями (обычно скорость распространения возрастает с увеличением глубины), определить местоположение горизонтов H1, H2 и Н 3 далеко не просто. Данные от каждого горизонта, в общем случае, разделяются, и анализируемая модель распространения сейсмической энергии (т.е. оценка скорости звука в подповерхностной области) в подповерхностной области применяется для получения оценки фактической точки отражения (или изображения) на горизонте, находящемся под поверхностью. Процесс миграции по глубине перед накоплением включает в себя одновременно суммирование множественных выборок с целью увеличения отношения сигнал/шум в данных (аналогичное общепринятой процедуре накопления), перемещение сейсмических событий с целью компенсации расстояния смещения между источником и приемником (аналогичное общепринятой процедуре нормальной коррекции разности времен вступления сейсмических волн) и перемещение сейсмических событий с целью компенсации наклонных сейсмических отражателей (аналогичное общепринятой процедуре миграции). Отображение сейсмических данных перед накоплением обычно реализуется с использованием интегральной формулы в пространственно-временной области. Хотя в этой заявке используется термин точка отражения, более понятным был бы 5 термин точка изображения, т.е. место в подповерхностной области, освещенное искомым лучом. Аналогично, траектория луча от источника до точки изображения и от точки изображения до приемника резко изменяет направление в точке изображения. Хотя это изменение направления именуется в этой заявке отражением, можно также считать это рассеянием в точке изображения. Модель распространения сейсмической энергии может не использовать никаких предположений, касающихся углов наклонения (наклона) отражателя вблизи точки изображения, и совершенно не обязательно,чтобы угол падения на сейсмический отражатель нисходящего луча от источника к точке изображения был равен и противоположен углу падения на сейсмический отражатель восходящего луча от точки изображения к приемнику. Изображение в любой точке вычисляется интегрированием данных, не прошедших накопление, по поверхности. В общем виде интеграл можно записать следующим образом: Изображение (хi) = W Данные (хs, хr, ts + tr)dSdR где хi, хs и хr - это положения изображения, источника и приемника, W (xi, хs, хr) - это весовая функция (возможно комплексная и частотно зависимая), которая является функцией положений источника, приемника и изображения, ats(xi, хs) и tr(xi, хr) - это времена прохождения,соответственно, от источника и приемника к точке изображения. Времена прохождения и весовые функции зависят от модели, которая представляет собой оценку свойств приповерхностной области. Интегрирование производят по всем координатам источника и приемника. Времена прохождения задают траекторию в области данных, по которой производится интегрирование. Времена интегрирования можно вычислять различными способами (например, способами конечных разностей, способами сеточного приближения времени прохождения, определения траектории лучей). Большинство этих способов позволяют решить волновое уравнение в высокочастотном приближении, разделив решение на две части, т.е. сначала решить уравнение эйконала для времен прохождения, а затем решить уравнение переноса для амплитуд. Уравнение эйконала - это уравнение, полученное асимптотическим расширением волнового уравнения. Это нелинейное дифференциальное уравнение,решениями которого являются времена прохождения. Все эти способы предусматривают либо явное, либо неявное вычисление траектории луча, которая представляет собой путь, по которому энергия проходит между источником или приемником и точкой изображения. При использовании способа, отвечающего настоящему изобретению, предпочтительно использовать способ определения траектории луча. Для многих способов определения траек 001707 6 тории луча характерно вычисление расстояния прохождения, поскольку они предполагают выбор из альтернативных траекторий луча или предоставляют выходные данные, которые позволяют легко вычислить расстояния прохождения, соответствующие различным траекториям луча. Однако в усложненной модели подповерхностной области может существовать несколько траекторий луча, соединяющих источник с приемником, и это значит, что функции времени прохождения могут быть многозначными. На фиг. 2 изображен один пример усложненной модели, которая обуславливает многозначные времена прохождения. Сейсмический источникS используется для анализа изображения I в точке, находящейся ниже поверхности 10 земли. Однако по обе стороны прямого пути 18 между источником S и изображением I-1 располагаются два соляных массива 12, 14. Соляные массивы 12, 14 преломляют сейсмические лучи и обеспечивают две дополнительные траектории луча 16, 20 между источником S и изображением I. Таким образом, помимо траектории прямого луча, который не проходит через соль, образуются траектории преломленных лучей, идущих от источника к точке изображения, обусловленные двумя соляными массивами. Чтобы правильно вычислить интеграл с учетом многозначного времени прохождения,нужно произвести суммирование по всем ветвям функций времени прохождения. Однако при наличии трех траекторий от источника S до изображения I и трех траекторий от изображения до приемника, суммирование нужно проводить по 9 ветвям. Это ведет к значительному удорожанию 3-мерной миграции по глубине перед накоплением. 9-кратное увеличение сложности вычисления значительно удлиняет анализ и,следовательно, значительно увеличивает его стоимость. По мере увеличения числа траекторий, соответственно, возрастает сложность. Если не использовать результаты по всем ветвям, придется делать выбор, какую ветвь или какие ветви использовать. Выбор должен дать траекторию интегрирования, каждой точке которой соответствует значительная энергия в рассеянном волновом поле, при этом траектория интегрирования должна давать хорошее приближение к непрерывному интегралу. Если интегрирование выполняется как взвешенное суммирование по выборочным данным, это значит,что траектория интегрирования должна быть хотя бы кусочно-непрерывной. Согласно методике, известной из уровня техники, которая кратко обсуждалась выше,используется первый приходящий сигнал, поскольку он легче поддается вычислению и наверняка является непрерывной функцией. Однако, как было упомянуто выше, первый приходящий сигнал может нести очень маленькую 7 энергию и, таким образом, не всегда является подходящим выбором. На фиг. 3 показано сейсмическое изображение после анализа, где по горизонтальной оси отложено расстояние (в футах) [в метрах] вдоль поверхности от исходной точки обследования. По вертикальной оси отложена глубина (в футах) [в метрах] под поверхностью, в данном случае, слоя воды, который простирается до глубины 6000 футов [1829 м]. При построении этого изображения использовались времена прохождения, полученные решением общеизвестного уравнения эйконала способом конечных разностей. Миграция перед накоплением осуществлялась с использованием первого приходящего сигнала на сейсмическом приемнике. Следует отметить, что многие детали на этой фигуре выглядят несколько нечетко. В частности, на изображении можно видеть два слоя,расположенные, соответственно, на глубине,примерно, 17000 футов [5182 м] и 27000 футов[8230 м]. Четкость этих слоев весьма низка, и имеются значительные части изображения этих слоев, где слои совсем неразличимы. Это, в частности, относится к верхнему слою (17000 футов [5182 м]) на расстояниях примерно от 33000 футов [10058 м] до 5000 футов [15240 м] от исходной точки обследования и к нижнему слою(на 27000 футов [8230 м]) приблизительно от 3000 [9144 м] до 47000 футов [14326 м] от исходной точки обследования. Имеется два больших соляных массива, первый из которых расположен на глубине от 6000 [1829 м] до 14000 футов [4267 м] и от 10000 [3048 м] до 37000 футов [11278 м] от исходной точки обследования. Второй соляной массив находится на глубине приблизительно от 6000 [1829 м] до 10000 футов [3048 м] и от 45000 [13716 м] до 75000 футов [22860 м] от исходной точки обследования. Согласно описанному выше со ссылкой на фиг. 2 соляные массивы могут обуславливать наличие траекторий преломленных лучей и это может быть причиной низкого уровня энергии первого приходящего сигнала. Поэтому изображения слоев в этой модели слабо обозначены. Кроме того, плохо обозначены образующие сетку детали на расстоянии примерно от 40000[12192 м] до 60000 футов [18288 м] и глубине примерно от 9000 [2743 м] до 17000 футов [5182 м]. Анализ, основанный на первом приходящем сигнале, также дает низкое разрешение этих сравнительно малых деталей. Другая методика, известная из уровня техники, состоит в выборе приходящего сигнала с максимальной амплитудой. Это наилучший для использования приходящий сигнал до тех пор,пока верны оценки амплитуды. Ранее приведенная ссылка демонстрирует, насколько улучшается отображение по глубине перед накоплением, если использовать времена прохождения сигнала с максимальной амплитудой, вычисленные путем определения траектории луча. 8 На фиг. 4 показано сейсмическое изображение, охватывающее те же расстояния и глубины, полученное на основании той же сейсмической модели, что и изображение, показанное на фиг. 3. В этом случае для миграции перед накоплением и последующего накопления выбирается приходящий сигнал с максимальной амплитудой. Можно видеть, что по сравнению с фиг. 3 почти каждая деталь изображения выглядит четче. Два горизонтальных слоя, находящиеся на глубине, соответственно, примерно 17000 футов [5182 м] и 27000 футов [8230 м] также воспроизведены значительно четче, чем на изображении, построенном на основании первого приходящего сигнала. Однако те части двух слоев, которые упоминались в описании предыдущего изображения, все еще воспроизведены несколько нечетко. В частности, нижний слой (27000 футов [8230 м]) воспроизведен несколько неясно в диапазоне расстояний примерно от 32000 [9574 м] до 44000 [13411 м] футов от исходной точки обследования. Кроме того,малые детали, расположенные в виде решетки на расстоянии примерно от 40000 [12192 м] до 60000 [18288 м] футов выглядят существенно четче. К сожалению, амплитуды, вычисленные в асимптотическом приближении, не обязательно являются хорошим приближением к амплитудам волнового поля конечной частоты. Амплитуды на высоких частотах гораздо чувствительнее к мелким деталям модели скоростей, чем волновое поле с конечной полосой частот. В результате, во-первых, алгоритм выбора приходящего сигнала может быстро переключаться от одной ветви к другой, во-вторых, изменения модели оказывают очень большое влияние на выбор ветви. Поскольку в общем случае мы можем опираться лишь на приближенную модель,сделанный выбор может очень отличаться от выбора, который был бы сделан в эталонной модели. Это с очевидностью ставит под угрозу точность анализа. Чтобы исправить эти недостатки, настоящее изобретение предусматривает новый критерий выбора времени прохождения сигнала, которое используется при создании изображения. Выбирается время прохождения, соответствующее лучу, имеющему наименьшую физическую длину. Этот луч часто, хотя и не обязательно, соответствует приходящему сигналу с наивысшей амплитудой. Важнее то, что он почти никогда не соответствует первому приходящему сигналу с низкой амплитудой, связанному с преломленной энергией. Кроме того, критерий физической длины обеспечивает значительно более высокую устойчивость по отношению к малым флуктуациям модели скоростей или к изменениям модели, поскольку физическая длина, как уже было отмечено, является интегральной величиной. 9 В отличие от времени прохождения первого приходящего сигнала, время прохождения кратчайшего луча не является непрерывной функцией. Однако, она является кусочнонепрерывной. Кроме того, ввиду большей стабильности, она, в основном, состоит из нескольких больших кусков с четко очерченными границами между ними. Напротив, использование максимальных амплитуд чаще всего приводит к значительно более разрывной траектории интегрирования, которая состоит из многочисленных малых сегментов, быстро переключаясь от одной ветви к другой и обратно. Хотя желательно иметь непрерывную функцию, кусочно-непрерывная функция, обеспечиваемая кратчайшим лучом, имеет другие преимущества (например, она соответствует более высокой энергии). На фиг. 5 показано еще одно изображение,которое, как и изображения на фиг. 3 и 4, получено с помощью той же сейсмической модели. В этом примере, однако, миграция перед накоплением производилась с использованием луча наименьшей длины, и не луча, соответствующего первому приходящему сигналу или приходящему сигналу с максимальной амплитудой. Сравнивая с фиг. 3, легко видеть, что миграция перед накоплением с использованием кратчайшего луча дает значительно более четкое изображение, чем при выборе первого приходящего сигнала. Преимущество этого изображения по сравнению с полученным на основании приходящего сигнала с максимальнойамплитудой, не столь очевидно, но, тем не менее, оно существует. В частности, улучшено воспроизведение слоев, находящихся на глубине, соответственно,около 17000 футов [5182 м] и 27000 футов [8230 м]. Воспроизведение верхнего слоя имеет более высокую четкость, особенно, на расстоянии,примерно, от 32000 футов [9574 м] до 50000 футов [15240 м] от исходной точки обследования. Улучшение изображения по сравнению с фиг. 4 еще виднее в случае нижнего слоя (27000 футов [8230 м]) . Часть этого слоя, находящаяся на расстоянии, примерно от 30000 [9144 м] до 43000 футов [13106 м] от исходной точки обследования выглядит явно отчетливее. В области, лежащей под промежутком между соляными массивами, эти слои существенно более непрерывны. Хотя преимущества, которые дает смена выбора приходящего сигнала с максимальной амплитудой на выбор приходящего сигнала кратчайшей траектории, не столь заметны, как преимущества отказа от выбора первого приходящего сигнала, однако, следует признать, что в любой точке изображение имеет либо такое же качество, либо еще лучше. Во всех сечениях миграции, показанных на фиг. 3-5, изображения искривленных событий,находящихся ниже 20000 футов [6096 м], представляют собой артефакты, обусловленные мно 001707 10 гократно отраженным приходящим сигналом, и их следует игнорировать. Выбранная траектория луча используется в миграции по глубине перед накоплением, чтобы, во-первых, установить время прихода акустического импульса и вычислить оценку амплитуды приходящего сигнала. Эти значения можно использовать, например, для определения, какую входную выборку из сейсмического следа, полученного в данном месте расположения приемника и связанного с акустическим импульсом, исходящим из источника, находящегося в данном положении, следует использовать для вычисления выходной выборки, связанной с точкой изображения, а также веса, который следует присваивать входной выборке. Хотя на фиг. 3, 4 и 5 показаны двухмерные сечения подповерхностной области, они были созданы с использованием трехмерных наборов сейсмических данных. Алгоритм миграции по глубине перед накоплением, используемый для получения этих сечений, является достоверным трехмерным способом, который не основан на предположении, что акустическая энергия, передаваемая от источника к приемнику, распространяется в вертикальной плоскости, соединяющей эти точки, или, что точка изображения должна находиться в этой вертикальной плоскости. Этот вывод о возможности распространения энергии вне плоскости не только повышает качество полученного сейсмического изображения, но и повышает сложность данных, которые можно правильно обрабатывать. Выбранную траекторию луча также можно использовать для обновления модели распространения сейсмической энергии в подповерхностной области (т.е. модели скоростей). Хотя сейсмические данные, используемые для демонстрации способа, отвечающего изобретению, представляют собой сейсмические данные общепринятого режима давлениедавление (Р-Р), способ, в соответствии с изобретением, можно использовать точно так же применительно к сейсмическим данным режима давление-сдвигающее усилие (P-S) (а также режима сдвигающее усилие-сдвигающее усилие(S-S) и других режимов переноса сейсмической энергии), просто предусматривая модель распространения сейсмической энергии в подповерхностной области, которая отвечает этим альтернативным режимам переноса сейсмической энергии. Сущность и объем настоящего изобретения не ограничивается описанными вариантами реализации, но охватывает любое явно или неявно раскрытое здесь изобретение и любое его обобщение. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ обработки сейсмических данных, использующий модель распространения 11 сейсмической энергии в подповерхностной области, включающий следующие этапы: задание положений сейсмического источника, сейсмического приемника и точки отражения в упомянутой модели распространения; идентификацию совокупности альтернативных траекторий луча, согласующихся с моделью распространения, которые начинаются в месте расположения сейсмического источника,отражаются в месте расположения точки отражения и оканчиваются в месте расположения сейсмического приемника; выбор траектории луча из упомянутой совокупности альтернативных траекторий луча,имеющей наименьшую длину; и использование выбранной траектории луча в последующей сейсмической обработке. 2. Способ по п.1, согласно которому упомянутый этап выбора траектории луча включает в себя вычисление расстояний прохождения для каждой из упомянутых альтернативных траекторий луча и сравнение упомянутых вычисленных расстояний прохождения для определения,какая из упомянутых альтернативных траекторий луча имеет минимальную длину. 3. Способ по любому из предшествующих пунктов, согласно которому совокупность альтернативных траекторий луча идентифицируется посредством процедуры определения траектории луча. 4. Способ по любому из предшествующих пунктов, согласно которому упомянутая модель распространения является трехмерной и упомянутая совокупность альтернативных траекторий луча не лежит в одной плоскости. 12 5. Способ по любому из предшествующих пунктов, согласно которому выбранная траектория луча используется для вычисления времени прихода акустического импульса. 6. Способ по любому из предшествующих пунктов, согласно которому выбранная траектория луча используется для оценки амплитуды приходящего сигнала. 7. Способ по любому из предшествующих пунктов, согласно которому выбранная траектория луча используется для выбора выборки из сейсмического следа, полученного в упомянутом месте расположения сейсмического приемника, связанного с акустическим импульсом,исходящим из места расположения сейсмического источника, причем упомянутая выборка используется для создания изображения точки отражения. 8. Способ по п.7, согласно которому выбранная траектория луча используется для вычисления веса, присваиваемого выборке. 9. Способ по любому из предшествующих пунктов, согласно которому упомянутая модель распространения моделирует перенос сейсмической энергии в режиме "давление - давление"(Р-Р) в подповерхностной геофизической области. 10. Способ по любому из предшествующих пп.1-8, согласно которому упомянутая модель распространения моделирует перенос сейсмической энергии в режиме "давление - сдвигающее усилие" (P-S) в подповерхностной геофизической области.