RU2300123C2 - Высокоразрешающее преобразование радона для обработки сейсмических данных - Google Patents
Высокоразрешающее преобразование радона для обработки сейсмических данных Download PDFInfo
- Publication number
- RU2300123C2 RU2300123C2 RU2004137818/28A RU2004137818A RU2300123C2 RU 2300123 C2 RU2300123 C2 RU 2300123C2 RU 2004137818/28 A RU2004137818/28 A RU 2004137818/28A RU 2004137818 A RU2004137818 A RU 2004137818A RU 2300123 C2 RU2300123 C2 RU 2300123C2
- Authority
- RU
- Russia
- Prior art keywords
- seismic data
- radon transform
- data
- coherence
- measurement
- Prior art date
Links
- 229910052704 radon Inorganic materials 0.000 title claims abstract description 118
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 title claims abstract description 118
- 238000012545 processing Methods 0.000 title claims abstract description 30
- 230000009466 transformation Effects 0.000 title abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 38
- 238000005259 measurement Methods 0.000 claims description 35
- 238000000926 separation method Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000009499 grossing Methods 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims description 4
- 230000010363 phase shift Effects 0.000 claims description 3
- 230000005713 exacerbation Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 abstract 1
- 239000013598 vector Substances 0.000 description 11
- 238000012217 deletion Methods 0.000 description 8
- 230000037430 deletion Effects 0.000 description 7
- 238000005070 sampling Methods 0.000 description 7
- 238000000844 transformation Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 3
- 230000017105 transposition Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000007717 exclusion Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Complex Calculations (AREA)
Abstract
Предложенные способ и компьютерная система относятся к геофизическим методам для обработки сейсмических данных. Сущность способа заключается в том, что находят высокоразрешающее преобразование Радона, предназначенное для применения к сейсмическим данным. Высокоразрешающее преобразование Радона регуляризируют, используя меру когерентности сейсмических данных. Предложено поддающееся оценке представление обработанных сейсмических данных. Мера когерентности сейсмических данных может включать в себя меру когерентности вдоль измерения сейсмических данных. Технический результат способа заключается в том, что сейсмические данные обрабатывают, используя высокоразрешающее преобразование Радона для улучшения требуемых признаков сейсмических данных. 4 н. и 15 з.п. ф-лы, 16 ил.
Description
Область техники
Это изобретение относится к вычислению высокоразрешающих преобразований Радона, используемых при обработке сейсмических данных.
Описание уровня техники
Преобразования Радона, наряду с преобразованием Фурье и некоторыми другими преобразованиями, являются частью способов, доступных геофизикам для моделирования и анализа сейсмических сигналов. В геофизике и в других прикладных областях были получены усовершенствованные преобразования Радона благодаря преодолению ограничений, обусловленных дискретизацией и шумовым содержанием данных.
Самые первые преобразования Радона, использованные в геофизике, были просто дискретными версиями результатов, полученных для непрерывных функций, дискретизированных на протяжении большого интервала. Алгоритмы дискретного преобразования Радона, полученные позднее (Beylkin, 1987, Hampson, 1986), дают возможность осуществлять точно инвертируемые преобразования в случае дискретных выборочных данных с ограниченной апертурой. Эти преобразования, в настоящее время называемые обычными, являются эффективными при разделении сигнала и шума при условии, что входные данные и их линейные интегралы (то есть данные в области преобразования Радона) дискретизированы без наложения зеркальных частот. На практике обработка входных данных с наложенными зеркальными частотами является важным вопросом, и усовершенствования существующей в настоящее время технологии добиваются двумя способами.
Первый способ заключается в выработке требований к дискретизации применительно к обычному преобразованию Радона (Schonewille and Duijndam, 2001, Hugonnet and Canadas, 1995, Marfurt, 1996) и соблюдение этих требований во время регистрации данных, интерполяции данных и предварительной обработки (Manin and Spitz, 1995).
Во втором способе априорную информацию относительно частей данных, которые искажены вследствие наложения зеркальных частот, включают в вычисление преобразования Радона. В самом простейшем виде априорную информацию вводят в виде диагонального члена регуляризации в решение методом наименьших квадратов для преобразования Радона. Когда все значения вдоль диагонали равны, результирующие решения в общем случае не позволяют разрешать неоднозначности, обусловленные входными данными с наложенными зеркальными частотами. Nichols (1992a, 1992b), Herrmann et al. (2000) и Hugonnet (2001) показали, что при соответствующей диагональной регуляризации получаются точные и эффективные преобразования Радона даже в случае, когда входные данные искажены вследствие наложения зеркальных частот.
В рамках оценки Байеса диагональная регуляризация может быть интерпретирована как постановка априорных распределений Гаусса в соответствие параметрам преобразования, в отношении которых также предполагается, что они должны быть статистически независимыми (Tarantola, 1987, Ulrych et al., 2001). Для других алгоритмов требуется меньше априорной информации благодаря исключению предположения относительно распределений Гаусса, но выражена тенденция к большей вычислительной сложности (Sacchi and Ulrych, 1995, Harlan et al., 1984, Thorson and Claerbout, 1985).
В предшествующем уровне техники относительно диагональной регуляризации при вычислении высокоразрешающих преобразований Радона имеются способы, которые описали Nichols (1992a), Herrmann et al. (2000) и Hugonnet et al. (2001). В этих случаях, исходя из соображений вычислительной эффективности, преобразование обычно применяют в пространственно-частотной области с зависимой от частоты диагональной регуляризацией. Однако член регуляризации не корректируют в зависимости от вычисленного преобразования Радона.
Nichols (1992a) получил веса регуляризации из меры когерентности вдоль временной оси частот данных. Мера когерентности для не искаженной зеркальными частотами сейсмической энергии обычно изменяется медленно с частотой. С другой стороны, сильные мелкомасштабные изменения меры когерентности часто используют для идентификации энергии, искаженной зеркальной частотами. Для получения весов Nichols (1992a) предложил сглаживание меры когерентности на интервале частот, содержащем частоты, которые ниже и выше частоты, подлежащей обработке.
Herrmann (2000) и Hugonnet (2001) получили веса рекурсивно на частоте, начиная с обычного преобразования Радона на самых низких частотах. Веса для частоты, на которой выполняют обработку, образуют по результатам высокоразрешающего преобразования Радона на более низких частотах.
Примерами на фиг.2, 3 и 4 иллюстрируются преобразования Радона, вычисленные при использовании различных схем регуляризации из предшествующего уровня техники. Фиг.2a-2d соответствуют обычному преобразованию Радона, регуляризированному с помощью скаляра, умноженного на единичную матрицу. На фиг.3a-3d веса регуляризации являются переменными членами вдоль диагонали, а соответствующее решение имеет более высокое разрешение. Однако на фиг.4 видно, что та же самая схема регуляризации, примененная к загрязненным шумом данным, дает субоптимальный результат.
Веса на фиг.4е не позволяют обнаруживать параболические оси синфазности волны на низких и высоких частотах. Веса надежно отображают приращения времени в данных только выше 15 кГц, тогда как в случае отсутствия шума, показанном на фиг.3е, приращения времени обнаруживаются от примерно 5 кГц. Кроме того, веса на фиг.4е имеют кажущиеся высокие значения на краях области преобразования, возможно шумы, обусловленные началом преобразования с обычного преобразования Радона (то есть с равными весами). Этими высокими значениями вызываются сильные шумы в области данных.
Описан новый способ для получения весов регуляризации, предназначенных для вычисления высокоразрешающего преобразования Радона.
В одном варианте осуществления изобретения представлен способ обработки сейсмических данных. Для применения к сейсмическим данным находят высокоразрешающее преобразование Радона. Высокоразрешающее преобразование Радона регуляризируют, используя меру когерентности сейсмических данных. Сейсмические данные обрабатывают, используя высокоразрешающее преобразование Радона для улучшения желательных признаков сейсмических данных. Предлагают поддающееся оценке представление обработанных сейсмических данных.
Мера когерентности сейсмических данных может включать в себя меру когерентности вдоль измерения сейсмических данных. Сейсмические данные могут быть подразделены в двумерный массив, который включает в себя одно измерение времени или глубины, и второе измерение из положения пространственной поверхности или угла. Кроме того, сейсмические данные могут быть подразделены в двумерный массив, включающий в себя одно измерение, включающееся в себя время или глубину, и другое измерение, выбранное из положений пространственной поверхности и углов.
Обработка сейсмических данных может включать в себя выполнение высокоразрешающего преобразования Радона относительно сейсмических данных при использовании меры когерентности сейсмических данных. Могут быть выделены первая область и вторая область в сейсмических данных, подвергнутых преобразованию Радона. Может быть выполнена инверсия высокоразрешающего преобразования Радона относительно выделенных, подвергнутых преобразованию Радона сейсмических данных.
Первая область может быть сигнальной областью, а вторая область может быть шумовой областью. Измерения сейсмических данных могут включать в себя многочисленные измерения сейсмических данных. Измерение сейсмических данных может также включать в себя частотную область.
В другом варианте осуществления изобретения обработка сейсмических данных при использовании высокоразрешающего преобразования Радона для улучшения требуемых признаков сейсмических данных может дополнительно включать в себя аппроксимацию комплексной траектории приращения времени отрезками. Сейсмические данные могут быть подразделены в локальные окна данных, согласованных с отрезками. Высокоразрешающее преобразование Радона может быть выполнено в каждом из локальных окон данных для улучшения требуемых признаков сейсмических данных, согласованных с отрезками. Осуществление высокоразрешающего преобразования Радона в каждом из локальных окон данных может включать в себя этапы вычисления преобразования с периодическими граничными условиями и применения разделения сигнала и шума, согласованного с нулевыми граничными условиями при использовании известных временных сдвигов между компонентами модели.
Регуляризация высокоразрешающего преобразования Радона при использовании меры когерентности сейсмических данных может включать в себя вычисление матрицы регуляризации, которое включает в себя применение фазового сдвига, определенного с помощью траектории приращения времени, к сейсмическим данным. Кроме того, она может включать в себя нормализацию мощности разреза вдоль траектории приращения времени для получения меры когерентности. Она также может включать в себя обострение меры когерентности или сглаживание меры когерентности по второму измерению сейсмических данных. Второе измерение сейсмических данных может включать в себя измерение, выбранное из пространственного положения и угла.
В другом варианте осуществления изобретение может быть реализовано посредством компьютерной системы, при этом компьютер включает в себя запоминающее устройство и процессор и выполняемое программное обеспечение, находящееся в запоминающем устройстве компьютера. Программное обеспечение работает совместно с процессором для нахождения высокоразрешающего преобразования Радона, предназначенного для применения к сейсмическим данным. Высокоразрешающее преобразование Радона регуляризируется при использовании меры когерентности сейсмических данных. Сейсмические данные обрабатываются при использовании высокоразрешающего преобразования Радона для улучшения требуемых признаков сейсмических данных.
Кроме того, программное обеспечение предусмотрено для осуществления высокоразрешающего преобразования Радона относительно сейсмических данных при использовании меры когерентности сейсмических данных. В сейсмических данных, подвергнутых преобразованию Радона, могут быть выделены первая область и вторая область. Относительно выделенных, подвергнутых преобразованию Радона сейсмических данных может быть использована инверсия высокоразрешающего преобразования Радона.
Программное обеспечение может также обеспечивать выполнение аппроксимации сложной траектории приращения времени отрезками. Сейсмические данные могут быть подразделены в локальные окна данных, согласованных с отрезками. Высокоразрешающее преобразование Радона может быть выполнено в каждом из локальных окон данных для улучшения требуемых признаков сейсмических данных, согласованных с большим числом отрезков.
В предложенном способе могут сохраняться характеристики разрешающей способности и исключения наложения зеркальных частот высокоразрешающих преобразований Радона и наибольшей эффективности способов, которые предложили Nichols и Herrmann, но способ может быть более устойчивым по сравнению с тем, который предложил Herrmann, и может обеспечивать лучшее разрешение по сравнению со способом, который предложит Nichols.
Описание чертежей
Другие признаки и преимущества изобретения станут более понятными при принятии во внимание нижеследующего описания с сопроводительными чертежами.
На чертежах:
Фиг.1 - (a) - синтетические тестовые данные, включающие в себя параболические оси синфазности волны с приращениями времени 0 и 400 мс при максимальном удалении 2000 м; (b) - тот же самый набор данных с добавленным случайным шумом;
Фиг.2 - иллюстрация подавления многократных волн обычным преобразованием Радона, примененным к данным из фиг.1a, свободным от шума; (А) - оцененные однократные волны; (В) - оцененные многократные волны; (С) - остаток: разность между входными данными и результатом применения к входным данным прямого и обратного преобразований Радона; и (D) - область Радона;
Фиг.3 - иллюстрация применения к данным из фиг.1a, свободным от шума, подавления многократных волн высокоразрешающим преобразованием Радона при использовании рекурсивного метода для получения весов; фиг. с (a) по (d) соответствуют фиг.2; на последней фиг.(e) показаны веса, использованные для высокоразрешающего преобразования Радона; веса отображены в зависимости от частоты (вертикальная ось) и параметра приращения времени (горизонтальная ось);
Фиг.4 - изображения, соответствующие фиг.3, но для зашумленных синтетических данных, показанных на фиг.1b;
Фиг.5 - изображения, соответствующие фиг.3, но при использовании весов, полученных на основе меры когерентности;
Фиг.6 - увеличенные изображения из фиг.4b (слева рекурсивная корректировка весов) и 5b (справа веса на основе меры когерентности);
Фиг.7 - изображения с плохой фокусировкой и шумами преобразования, связанными с обработкой гиперболической оси синфазности при использовании одного глобального окна данных и параболического преобразования Радона; (а) - остаток; (b) - данные в области параболического преобразования Радона; (с) - веса, использованные для высокоразрешающего преобразования Радона;
Фиг.8 - графики, иллюстрирующие отображение локального окна время-удаление в области параболического преобразования Радона;
Фиг.9 - результаты обработки тех же самых входных данных, что и на фиг.7, в пределах трех отдельных окон удалений; (а) - разность между входными данными и результатом прямого и обратного преобразований Радона, примененных к входным данным; (b, c, d) - данные в области параболического преобразования Радона для каждого из трех окон входных данных; (e, f, g) - веса, использованные для высокоразрешающего преобразования Радона для случаев каждого из трех окон входных данных;
Фиг.10 - пример полевых данных, обработка путем обычного преобразования Радона; (а) - входные данные; (b) - оцененные однократные волны; (с) - оцененные многократные волны; (d) - входные данные минус оцененные однократные волны; (е) - остаток от прямого и обратного преобразований: входные данные минус сумма оцененных однократных волн и многократных волн; (f) - данные в области параболического преобразования Радона;
Фиг.11 - пример полевых данных, как на фиг.10, при этом веса вычислены рекурсивно;
Фиг.12 - пример полевых данных, как на фиг.10, при этом веса определены на основе меры когерентности;
Фиг.13 - иллюстрация сравнения весов: (13а, слева) способ на основе меры когерентности и (13b, справа) рекурсивный способ;
Фиг.14 - рабочая диаграмма способа преобразования Радона;
Фиг.15 - (а) - сечения данных со схематичным отображением локальных окон данных в пределах каждого сечения, подготовленного для ввода в двумерное преобразование Радона; (b) - схематичное представление оси синфазности сейсмических данных волны (криволинейной формы) и траекторий приращения времени (прямые линии) для преобразования Радона; (с) схематичное представление области преобразования Радона. После суммирования вдоль прямых линий вертикальная линия в области преобразования будет закрашиваться, а большая часть энергии будет вноситься из траектории приращения времени, которые являются касательными по отношению к оси синфазности; и
Фиг.16 - рабочая диаграмма способа обработки сейсмических данных.
В настоящем изобретении представлена новая схема для вычисления регуляризованного высокоразрешающего преобразования Радона.
Что касается фиг.16, то сейсмические данные регистрируют 1601, используя стандартные промышленные способы и группируя их по методу общей поверхности (например, в виде сейсмограммы общей глубинной точки) или приповерхностной точки (например, в виде угловых сейсмограмм, образованных во время миграции). На фиг.15а показан двумерный массив 1502, образованный путем группирования по одному измерению в виде времени или глубины и по другому измерению в виде пространственного положения (удаления) или угла. Измерение времени или глубины обычно дискретизировано с большой плотностью и без наложения зеркальных частот, тогда как наложение зеркальных частот может наблюдаться вследствие недостаточной дискретизации вдоль пространственного положения (удаления) или углового измерения данных. В других вариантах осуществления с помощью группировок различных видов можно создавать для обработки трехмерные или более высокого порядка массивы.
Для преобразования данных траектории приращения времени для преобразования Радона выбирают 1602, используя приращение времени волн в данных и выводы, касающиеся вычислительной эффективности. Как показано на фиг.15b и 15с, типовыми вариантами траекторий приращения времени являются линейные и параболические траектории. Другие траектории приращений времени, такие, как гиперболические или более сложные траектории, также могут быть выбраны.
Во время преобразования Радона данные преобразуют 1605 в частотную область в пределах хорошо дискретизированного измерения (времени или глубины). Вектор входных данных для преобразования Радона на частоте ω обозначают d(x, ω), кратко d. Данные индексируют положениями в пространстве, x, часто представляемыми в виде расстояния (удаления) между источником и приемником. Вектор d данных может представлять данные из всей сейсмограммы или из локального (время-удаление) окна данных 1501, показанных на фиг.15а. Конкретные вопросы, относящиеся к обработке в локальных окнах данных рассмотрены ниже в разделе под названием "Критерии для разделения сигнала и шума в локальных окнах данных".
Результатом преобразования Радона на частоте ω является вектор m(q, ω) или кратко m, также называемый модельным вектором. Элементы m индексируют параметрами q, определяющими выбранные траектории приращения времени преобразования Радона. В случае параболического преобразования Радона параметры q интерпретируют как кривизны парабол при нулевом удалении. Обычно данные, вводимые в преобразование, представляют собой двумерные сечения. Однако описание методов преобразования Радона, использованных в некоторых вариантах осуществления настоящего изобретения, применимо к данным с большим числом измерений, и в этом случае параметры кривизны больше не являются скалярами, а векторами.
Преобразование Радона является линейным преобразованием относительно данных и модельных векторов. Связь между векторами модели и данных может быть выражена в матричной форме с помощью уравнения 1 (Beylkin, 1987):
где А - матрица размера (M, N);
М - число элементов (комплексных скаляров) в векторе d данных; и
N - число элементов в модельном векторе m.
Преобразование Радона можно найти как решение системы линейных уравнений, определенной в уравнении (1). Чтобы гарантировать единственное и устойчивое решение системы линейного уравнения (1), обычно задают целевую функцию, J(m), связанную с уравнением (1), которая является суммой члена регуляризации и члена подбора данных:
где верхний индекс Н означает сопряженную транспозицию вектора или матрицы; и
W-1 - матрица регуляризации.
Наиболее общим выбором для W-1 является диагональная матрица с постоянными элементами по диагонали, например W-1=αI, где I является единичной матрицей N×N. Этот выбор регуляризации связан с обычным преобразованием Радона и называется методом "затухающих наименьших квадратов". Он обеспечивает получение устойчивого решения, но требует не искаженных дискретизацией входных данных с большой апертурой.
В одном варианте осуществления настоящего изобретения веса регуляризации получают вычислением 1603 меры когерентности вдоль оси удаления:
где АН - матрица сопряженной транспозиции относительно А;
d(ω) - вектор данных может содержать данные из всей сейсмограммы или из локального окна время-удаление сейсмограммы;
α - коэффициент, который используют для повышения обострения меры когерентности; и
ε - множитель придания устойчивости и/или нормировочный.
Матрицу сопряженной транспозиции относительно А задают в членах траекторий приращения времени, а в данные вводят соответствующие фазовые сдвиги.
Затем меру когерентности можно сгладить в полосе частот, используя стандартную схему взвешенного усреднения (Marple, 1987), при этом элементы диагональной весовой матрицы W вычисляют как:
где коэффициенты αk зависят от числа частот, используемых при сглаживании, а не от действительных значений частот ω0, ω1 и т.д.
Типовыми вариантами для коэффициентов αk являются постоянные, треугольные и гауссовы веса. Предпочтительно, чтобы диапазон частот, в котором сглаживают меру когерентности, предпочтительно был широкой полосой частот и включал в себя высокие частоты. Для получения диагональной весовой матрицы, основанной на мере когерентности, могут быть использованы сглаживающие функции других видов, известные специалистам в области техники, к которой относится изобретение.
В дополнение к этому, аналогично обычной практике обработки сейсмических данных, меру когерентности можно подвергнуть фильтрации для повышения устойчивости путем пространственного усреднения по небольшому числу соседних положений общей глубинной точки или путем отбрасывания значений меры когерентности, меньших порогового значения.
После преобразования 1605 сейсмических данных при использовании высокоразрешающего преобразования Радона могут быть выделены сигнальные и шумовые области преобразованных данных (это также известно как мьютинг). Простейший способ заключается в выборе вручную геометрии разделения (функции мьютинга) сигнальных и шумовых областей и удаления данных на одной стороне геометрии. Если сейсмические данные выделяют в локальных окнах, мьютинг можно осуществлять способом, описанным в разделе под названием "Эффективное осуществление выделения (мьютинга) сигнала и шума для преобразований из локальных окон данных".
После удаления шума сейсмические данные повторно преобразуют 1607, используя обратное высокоразрешающее преобразование Радона. Затем сейсмические данные могут быть представлены 1608 в виде поддающегося оценке представления, например на компьютерной распечатке или на экране дисплея, или обработаны дополнительно до отображения или представления пользователю.
В различных вариантах осуществления этого изобретения программное обеспечение в компьютере или компьютерной системе можно использовать для осуществления высокоразрешающего преобразования Радона, представленного в этом изобретении и известного в области техники, к которой относится изобретение.
Параболическое преобразование Радона основано на волнах в данных, имеющих параболическое приращение времени и постоянную амплитуду после введения поправки на нормальное приращение времени. Однако на практике траектории приращения времени в данных являются сложными и необязательно параболическими. В другом варианте осуществления этого изобретения сложная траектория приращения времени может быть аппроксимирована путем объединения локальных отрезков парабол. Рассмотрение, приведенное ниже, относится к параболам, но применимо к любым отрезкам, обычно используемым при преобразованиях Радона.
Локальные окна данных определяют путем деления каждого измерения данных на отрезки, которые перекрывают весь диапазон этого измерения. Параболическое преобразование Радона может быть применено к данным в каждом из локальных окон. В случае, когда параболическое преобразование Радона применяют в каждом локальном окне, предположения относительно параболического приращения времени и постоянной амплитуды удовлетворяются лучше, чем в пределах всего диапазона этого измерения. В отличие от обычного преобразования Радона использование высокоразрешающего параболического преобразования Радона обеспечивает возможность отделения в локальном окне сигнала от шумовых данных. Используя отрезки парабол, которыми сложная траектория хорошо аппроксимируется в пределах локального окна, можно эффективно выразить разделение сигнала и шума в зависимости от сложной траектории приращения времени. Благодаря моделированию сложных траекторий приращения времени критерий для выбора требуемых признаков в данных может быть более эффективным.
На фиг.7 показана сейсмическая волна с остаточным гиперболическим приращением времени, которая плохо сфокусирована в случае параболического преобразования Радона, вычисленного по всему диапазону удалений. После разделения диапазона удалений на три части и вычисления параболического преобразования Радона отдельно в пределах каждой части каждое из преобразований оказывается хорошо сфокусированным, что показано на фиг.9.
Однако отметим, что при каждом из зависящих от диапазона удалений преобразований на фиг.9 гипербола фокусируется несколько в ином месте в пространстве модели. Наилучшие подобранные параболы в каждом окне удаления имеют несколько иные параметры (кривизну и время пробега по нормальному лучу). Это заключение подразумевает, что функции мьютинга, используемые для разделения волн в области преобразования, должны также изменяться как функция диапазона удаления данных, используемых для вычисления преобразования.
Для управления изменениями параметров мьютинга с удалением применим следующий подход. Во-первых, начнем с функции мьютинга, q(τ,x=0), для окна, включающего в себя нулевое удаление. При заданном времени τ пробега по нормальному лучу обозначим через q0 соответствующее положение мьютинга, q0=q(τ,x=0). Параметрами (τ,q0) задается парабола t=τ+q0x2/xmax 2, а также и гипербола с тем же временем пробега по нормальному лучу и кривизной при нулевом удалении, что и парабола:
В таком случае соотношение для мьютингов между окнами удалений:
где qref - требуемое положение мьютинга для окна данных, связанных с удалением xref, и для времени τ пробега по нормальному лучу.
Параметр qref приращения времени получают путем согласования в точке xref градиентов параболы (τ,qref) и гиперболы (τ,q0).
В более общем случае с помощью процедуры, описанной выше, находятся локальные критерии разделения сигнала и шума, согласованные с глобальным критерием разделения сигнала и шума и выраженные в членах сложной (например, гиперболической) функции мьютинга. Новым признаком по сравнению с предшествующим уровнем техники является использование траектории приращения времени (например, гиперболической) для глобального критерия разделения сигнала и шума, которая не является частью траекторий приращения времени, используемых для вычисления преобразования Радона (например, параболических). Использование этого нового признака увеличивает гибкость и повышает точность разделения сигнала и шума при сохранении вычислительной эффективности. Фактически при улучшенном разделении сигнала и шума может улучшаться разделение требуемых признаков в сейсмических данных.
Согласно другому варианту осуществления этого изобретения может быть вычислен оптимальный временной сдвиг между окнами в области (x-t) данных и в области (q, τ) преобразования Радона, как показано на фиг.8.
Рассмотрим окно время-удаление, задаваемое временным интервалом [t1, t2] и интервалом [x1, x2] удалений, и предположим, что xref является соответствующим эталонным удалением для окна (обычно это центральное удаление). Для заданного значения q времена пробега по нормальному лучу параболы, которая пересекает эталонное удаление, в пределах окна лежат в интервале [τ1, τ2], где τ1=t1-Δτ(q) и τ2=t2-Δτ(q), при этом Δτ(q)=q·xref 2/xmax 2. Можно видеть, что оптимальное окно в τ является предшествующим (для положительного q) в сопоставлении с оптимальным окном в t на величину временного сдвига Δτ(q). Отметим, что временной сдвиг является линейным от q, и поэтому прямоугольное окно в t-x описывается параллелограммом в τ-q. Когда q и xref большие, временной сдвиг Δτ может быть значительно больше, чем длина окна во временной области.
Поскольку при обработке в частотной области подразумевается, что данные и модель являются периодическими с периодом t2-t1, модель является эффективно определенной во все моменты времени, и мы просто выбираем соответствующее (в виде параллелограмма) окно из этой модели, зная сдвиг Δτ(q). Фактически процедура, описанная здесь, позволяет эффективно вычислять преобразование, используя периодические граничные условия, наряду с применением разделения сигнала и шума, как если бы граничные условия имели нулевые значения (предпочтительное описание) за пределами интервала [t1, t2] входных данных.
Для иллюстрации вариантов осуществления этого изобретения на фиг.1b показаны синтетические (с добавлением случайного шума). Для случая данных с добавленным шумом результаты отображены на фиг.5а-5е. Веса, основанные на мере когерентности, получены в частотном диапазоне от 70 до 90 Гц, и они предполагаются постоянными во всей полосе частот. В данном случае разделение однократных волн и многократных волн улучшено благодаря более точным весам. При этом на краях области преобразования высокие значения отсутствуют.
Результаты из фиг.4b и 5b в большем масштабе отражены на фиг.6. Ясно видно улучшенное разделение однократных волн и многократных волн при использовании весов на основе меры когерентности (фиг.6b), а не рекурсивно корректируемых весов (фиг.6а).
В качестве полевых данных для этого примера использована сейсмограмма общей глубинной точки (фиг.10а) с шагом между трассами 100 м и максимальным удалением около 6 км. Результаты обычного преобразования Радона (фиг.10b-10f) показывают, что встречаются обычные проблемы: плохая фокусировка в области преобразования, значительная энергия на краях области преобразования, шумы дискретизации в областях преобразования и данных.
Результаты, полученные на основании настоящего изобретения, отображены на фиг.12b-12f. Соответствующие веса на основе меры когерентности показаны на фиг.13а. Как и ожидалось, шумы дискретизации и преобразования снижены, а оцененные однократные волны (фиг.12b) улучшены по сравнению с результатами, полученными на основании обычного преобразования Радона (фиг.10b).
Результаты для высокоразрешающего преобразования Радона при использовании рекурсивно корректируемых весов отображены на фиг.11а-11f, при этом соответствующие веса показаны на фиг.13b. Эти результаты не столь хороши, как результаты, полученные с весами на основе меры когерентности. Это особенно видно на фиг.11f (область преобразования) и 11b (оцененные однократные волны), где часть энергии помех дискретизации и шумов на краях области преобразования теперь снова видна и может быть отнесена на счет субоптимальных весов, оцененных с помощью рекурсивной схемы (фиг.13b).
Описаны несколько вариантов осуществления настоящего изобретения. Тем не менее должно быть понятно, что различные модификации могут быть сделаны без отступления от сущности и объема изобретения. Поэтому другие варианты осуществления находятся в рамках объема нижеследующей формулы изобретения.
Claims (19)
1. Способ обработки сейсмических данных, включающий в себя
измерение сейсмических данных, при этом указанные данные подразделяют в двумерный массив, содержащий одно измерение, выбранное из группы, состоящей из времени и глубины, и второе измерение, выбранное из группы, состоящей из пространственного положения поверхности и угла,
нахождение высокоразрешающего преобразования Радона для применения к сейсмическим данным, включающее в себя
регуляризацию высокоразрешающего преобразования Радона при использовании меры когерентности сейсмических данных и
обработку сейсмических данных при использовании высокоразрешающего преобразования Радона для улучшения требуемых признаков сейсмических данных; и
обеспечение поддающегося оценке представления обработанных сейсмических данных.
2. Способ по п.1, в котором мера когерентности сейсмических данных содержит меру когерентности вдоль измерения сейсмических данных.
3. Способ по п.1, в котором сейсмические данные подразделяют в массив, содержащий одно измерение, выбранное из группы, состоящей из времени и глубины, и второе измерение, содержащее большее количество измерений, выбранных из группы, состоящей из положения пространственной поверхности и углов.
4. Способ по п.1, в котором обработка сейсмических данных дополнительно включает в себя
выполнение высокоразрешающего преобразования Радона относительно сейсмических данных при использовании меры когерентности сейсмических данных;
выделение первой области и второй области в сейсмических данных, подвергнутых преобразованию Радона; и
выполнение инверсии высокоразрешающего преобразования Радона относительно выделенных подвергнутых преобразованию Радона сейсмических данных.
5. Способ по п.3, в котором первая область представляет собой сигнальную область, а вторая область представляет собой шумовую область.
6. Способ по п.2, в котором измерение сейсмических данных содержит большое количество измерений сейсмических данных.
7. Способ по п.2, в котором измерение сейсмических данных содержит частотную область.
8. Способ по п.1, в котором обработка сейсмических данных при использовании высокоразрешающего преобразования Радона для улучшения требуемых признаков сейсмических данных дополнительно включает в себя
аппроксимацию сложной траектории приращения времени большим количеством отрезков;
разделение сейсмических данных в большое количество локальных окон данных, согласованных с большим количеством отрезков; и
выполнение высокоразрешающего преобразования Радона в каждом из локальных окон данных для улучшения требуемых признаков сейсмических данных, согласованных с большим количеством отрезков.
9. Способ по п.8, в котором выполнение высокоразрешающего преобразования Радона в каждом из локальных окон данных дополнительно включает в себя этапы вычисления преобразования с периодическими граничными условиями и применения разделения сигнала и шума, согласованного с нулевыми граничными условиями при использовании известных временных сдвигов между компонентами модели.
10. Способ по п.1, в котором регуляризация высокоразрешающего преобразования Радона при использовании меры когерентности сейсмических данных дополнительно включает в себя вычисление матрицы регуляризации, включающее в себя применение фазового сдвига, определенного с помощью траектории приращения времени, к сейсмическим данным.
11. Способ по п.1, в котором регуляризация высокоразрешающего преобразования Радона при использовании меры когерентности сейсмических данных дополнительно включает в себя нормализацию мощности разреза вдоль траектории приращения времени для получения меры когерентности.
12. Способ по п.1, в котором регуляризация высокоразрешающего преобразования Радона при использовании меры когерентности сейсмических данных дополнительно включает в себя обострение меры когерентности.
13. Способ по п.1, в котором регуляризация высокоразрешающего преобразования Радона при использовании меры когерентности сейсмических данных дополнительно включает в себя сглаживание меры когерентности по второму измерению сейсмических данных.
14. Способ по п.13, в котором второе измерение сейсмических данных содержит измерение, выбранное из пространственного положения и угла.
15. Способ для обработки сейсмических данных, включающий в себя
измерение сейсмических данных, при этом указанные данные подразделяют в двумерный массив, содержащий одно измерение, выбранное из группы, состоящей из времени и глубины, и второе измерение, выбранное из группы, состоящей из пространственного положения поверхности и угла,
нахождение высокоразрешающего преобразования Радона для применения к сейсмическим данным, включающее в себя
регуляризацию высокоразрешающего преобразования Радона при использовании меры когерентности вдоль измерения указанных сейсмических данных, при этом измерение сейсмических данных содержит частотную область;
обработку сейсмических данных при использовании высокоразрешающего преобразования Радона для улучшения требуемых признаков сейсмических данных, при этом обработка сейсмических данных включает в себя
выполнение высокоразрешающего преобразования Радона относительно сейсмических данных при использовании меры когерентности;
выделение сигнальной области и шумовой области в сейсмических данных, подвергнутых преобразованию Радона;
выполнение инверсии высокоразрешающего преобразования Радона относительно выделенных подвергнутых преобразованию Радона сейсмических данных и
обеспечение поддающегося оценке представления обработанных сейсмических данных.
16. Компьютерная система для обработки сейсмических данных, содержащая
компьютер, при этом компьютер содержит запоминающее устройство и процессор; и
реализующее программное обеспечение, находящееся в запоминающем устройстве компьютера, при этом программное обеспечение работает совместно с процессором для
нахождения высокоразрешающего преобразования Радона для применения к сейсмическим данным, полученным при измерении, при этом указанные данные подразделяют в двумерный массив, содержащий одно измерение, выбранное из группы, состоящей из времени и глубины, и второе измерение, выбранное из группы, состоящей из пространственного положения поверхности и угла,
при этом указанное преобразование Радона включает в себя
регуляризацию высокоразрешающего преобразования Радона при использовании меры когерентности сейсмических данных и
обработки сейсмических данных при использовании преобразования Радона для улучшения требуемых признаков сейсмических данных.
17. Система по п.16, в которой обработка сейсмических данных при использовании высокоразрешающего преобразования Радона дополнительно включает в себя
выполнение высокоразрешающего преобразования Радона относительно сейсмических данных при использовании меры когерентности сейсмических данных;
выделение первой области и второй области в сейсмических данных, подвергнутых преобразованию Радона; и
выполнение инверсии высокоразрешающего преобразования Радона относительно выделенных подвергнутых преобразованию Радона сейсмических данных.
18. Компьютерная система по п.16, в которой обработка сейсмических данных при использовании высокоразрешающего преобразования Радона дополнительно включает в себя
аппроксимацию сложной траектории приращения времени большим количеством отрезков;
разделение сейсмических данных в большое количество локальных окон данных, согласованных с большим количеством отрезков; и
выполнение высокоразрешающего преобразования Радона в каждом из локальных окон данных для улучшения требуемых признаков сейсмических данных, согласованных с большим количеством отрезков.
19. Компьютерная система для обработки сейсмических данных, содержащая
компьютер, при этом компьютер содержит запоминающее устройство и процессор; и
средство для обработки сейсмических данных на компьютере, полученных при измерении, при этом указанные данные подразделяют в двумерный массив, содержащий одно измерение, выбранное из группы, состоящей из времени и глубины, и второе измерение, выбранное из группы, состоящей из пространственного положения поверхности и угла, при использовании высокоразрешающего преобразования Радона и меры когерентности сейсмических данных для улучшения требуемых признаков сейсмических данных.
Applications Claiming Priority (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US38320902P | 2002-05-24 | 2002-05-24 | |
| US60/383,209 | 2002-05-24 | ||
| US10/165,439 | 2002-06-07 | ||
| US10/165,439 US6636810B1 (en) | 2002-05-24 | 2002-06-07 | High-resolution Radon transform for processing seismic data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2004137818A RU2004137818A (ru) | 2005-05-27 |
| RU2300123C2 true RU2300123C2 (ru) | 2007-05-27 |
Family
ID=28794021
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2004137818/28A RU2300123C2 (ru) | 2002-05-24 | 2003-05-19 | Высокоразрешающее преобразование радона для обработки сейсмических данных |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US6636810B1 (ru) |
| EP (1) | EP1508058B1 (ru) |
| CN (1) | CN100385253C (ru) |
| AT (1) | ATE514106T1 (ru) |
| AU (1) | AU2003241526B2 (ru) |
| NO (1) | NO333954B1 (ru) |
| RU (1) | RU2300123C2 (ru) |
| WO (1) | WO2003100464A1 (ru) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2549213C2 (ru) * | 2007-11-14 | 2015-04-20 | ТЕРРАСПАРК ДЖИОСАЙЕНСИЗ, ЭлЭлСи | Обработка сейсмических данных |
Families Citing this family (37)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6691039B1 (en) * | 2002-08-30 | 2004-02-10 | John M. Robinson | Removal of noise from seismic data using improved radon transformations |
| US7366054B1 (en) * | 2002-08-30 | 2008-04-29 | John M. Robinson | Tau-P filters for removal of noise from seismic data |
| US6832161B1 (en) * | 2003-08-15 | 2004-12-14 | Westerngeco, L.L.C. | Method for attenuating water layer multiples |
| RU2333515C1 (ru) | 2004-04-07 | 2008-09-10 | Вестернджеко Сайзмик Холдингз Лимитед | Быстрое трехмерное прогнозирование кратных волн от поверхности земли |
| US7561491B2 (en) * | 2005-03-04 | 2009-07-14 | Robinson John M | Radon transformations for removal of noise from seismic data |
| WO2006129204A2 (en) | 2005-05-31 | 2006-12-07 | Devgen Nv | Rnai for control of insects and arachnids |
| US7433265B2 (en) * | 2005-10-04 | 2008-10-07 | Fairfield Industries, Inc. | Converted wave energy removal from seismic data |
| US7881501B2 (en) * | 2005-12-13 | 2011-02-01 | Calgary Scientific Inc. | Local dominant wave-vector analysis of seismic data |
| GB2460013B (en) * | 2008-03-31 | 2010-10-13 | Statoilhydro Asa | A method for reflection time shift matching a first and a second set of seismic reflection data |
| CA2756231C (en) * | 2009-03-31 | 2017-11-21 | Shell Internationale Research Maatschappij B.V. | Method for stabilization of low frequencies in tau-p domain filtering and deconvolution |
| US8341984B2 (en) * | 2009-05-27 | 2013-01-01 | Westerngeco L.L.C. | Estimating velocities with uncertainty |
| US8315124B2 (en) * | 2009-07-20 | 2012-11-20 | Saudi Arabian Oil Company | System and method for suppression of seismic multiple reflection signals |
| US8774468B2 (en) * | 2009-09-08 | 2014-07-08 | Schlumberger Technology Corporation | Dynamic shape approximation |
| US9014985B2 (en) * | 2010-04-30 | 2015-04-21 | Chevron U.S.A. Inc. | System and method for compensating time and offset varying near-surface effects in seismic data background |
| CN102401908B (zh) * | 2010-09-07 | 2013-09-04 | 中国石油天然气集团公司 | 一种利用不同模加权稀疏的抛物拉东变换压制多次波的方法 |
| US20120215453A1 (en) * | 2011-02-22 | 2012-08-23 | Cggveritas Services Sa | Device and method for multi-dimensional coherency driven denoising data |
| US9405027B2 (en) | 2012-01-12 | 2016-08-02 | Westerngeco L.L.C. | Attentuating noise acquired in an energy measurement |
| US10317545B2 (en) * | 2012-03-12 | 2019-06-11 | Schlumberger Technology Corporation | Methods and apparatus for waveform processing |
| WO2013163131A1 (en) | 2012-04-23 | 2013-10-31 | Westerngeco Llc | Attenuating noise acquired in an energy measurement |
| CN102879824B (zh) * | 2012-09-07 | 2015-03-04 | 清华大学 | 一种基于迭代收缩的快速稀疏Radon变换方法 |
| EP2914981A4 (en) | 2012-11-01 | 2015-11-25 | Halliburton Energy Services Inc | DEVICE, SYSTEMS AND METHODS FOR DIFFERENTIAL PHASE SIMILARNESS |
| NO337238B1 (no) * | 2013-01-15 | 2016-02-22 | Fugro Geoteam As | Fremgangsmåte for fjerning av seismisk interferens |
| JP6542228B2 (ja) * | 2013-12-19 | 2019-07-10 | メルク、パテント、ゲゼルシャフト、ミット、ベシュレンクテル、ハフツングＭｅｒｃｋ Ｐａｔｅｎｔ ＧｍｂＨ | 複素環式スピロ化合物 |
| US20160320508A1 (en) * | 2014-01-13 | 2016-11-03 | Cgg Services Sa | Devices and methods for attenuation of turn noise in seismic data acquisition |
| US20150198729A1 (en) * | 2014-01-13 | 2015-07-16 | Cgg Services Sa | Regularization of spatially aliased seismic data |
| US10670757B2 (en) * | 2014-02-26 | 2020-06-02 | Pgs Geophysical As | Methods and systems for quantifying coherency and constraining coherency-based separation in simultaneous shooting acquisition |
| CN104849760B (zh) * | 2015-06-02 | 2017-07-14 | 中国矿业大学(北京) | 基于Tau‑p变换的矢量中值滤波方法 |
| CN104932010B (zh) * | 2015-06-09 | 2018-08-07 | 中国海洋石油集团有限公司 | 一种基于近道镶边稀疏Radon变换的绕射波分离方法 |
| CN104914466B (zh) * | 2015-06-26 | 2017-12-12 | 中国石油大学(华东) | 一种提高地震资料分辨率的方法 |
| US10032077B1 (en) | 2015-10-29 | 2018-07-24 | National Technology & Engineering Solutions Of Sandia, Llc | Vehicle track identification in synthetic aperture radar images |
| US11016212B2 (en) | 2017-04-11 | 2021-05-25 | Saudi Arabian Oil Company | Compressing seismic wavefields in three-dimensional reverse time migration |
| US10564304B2 (en) * | 2017-05-11 | 2020-02-18 | Saudi Arabian Oil Company | Processing methodology for full-waveform sonic wavefield separation |
| US10983236B2 (en) | 2017-06-20 | 2021-04-20 | Saudi Arabian Oil Company | Super-resolution radon transform based on thresholding |
| US10684382B2 (en) | 2018-01-23 | 2020-06-16 | Saudi Arabian Oil Company | Generating target-oriented acquisition-imprint-free prestack angle gathers using common focus point operators |
| US11353612B2 (en) | 2019-03-11 | 2022-06-07 | Saudi Arabian Oil Company | Nonstationary maximum likelihood method to estimate dispersion spectra for full waveform sonic logging |
| US11353609B2 (en) | 2019-12-20 | 2022-06-07 | Saudi Arabian Oil Company | Identifying geologic features in a subterranean formation using angle domain gathers sampled in a spiral coordinate space |
| US11656378B2 (en) | 2020-06-08 | 2023-05-23 | Saudi Arabian Oil Company | Seismic imaging by visco-acoustic reverse time migration |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4760563A (en) | 1986-01-09 | 1988-07-26 | Schlumberger Technology Corporation | Seismic exploration using exactly invertible discrete transformation into tau-p space |
| US5521881A (en) * | 1994-09-02 | 1996-05-28 | Exxon Production Research Company | Method of processing seismic data having multiple reflection noise |
| GB9703529D0 (en) * | 1997-02-20 | 1997-04-09 | Geco As | A method of processing seismic data signals |
| US6061298A (en) * | 1998-06-08 | 2000-05-09 | Pgs Tensor, Inc. | Method of and system for processing multicomponent seismic data |
| US6094620A (en) * | 1999-01-08 | 2000-07-25 | Exxonmobil Upstream Research Company | Method for identifying and removing multiples from seismic reflection data |
| US5970023A (en) * | 1999-02-19 | 1999-10-19 | Baker Hughes Incorporated | Reducing aliasing artifacts in seismic data processing using sharpened (TAU-P-Q) transforms |
| US6018500A (en) * | 1999-02-19 | 2000-01-25 | Baker Hughes Incorporated | Radon moveout |
-
2002
- 2002-06-07 US US10/165,439 patent/US6636810B1/en not_active Expired - Fee Related
-
2003
- 2003-05-19 AU AU2003241526A patent/AU2003241526B2/en not_active Ceased
- 2003-05-19 RU RU2004137818/28A patent/RU2300123C2/ru not_active IP Right Cessation
- 2003-05-19 WO PCT/US2003/015882 patent/WO2003100464A1/en not_active Application Discontinuation
- 2003-05-19 EP EP03731266A patent/EP1508058B1/en not_active Expired - Lifetime
- 2003-05-19 AT AT03731266T patent/ATE514106T1/de not_active IP Right Cessation
- 2003-05-19 CN CNB038119188A patent/CN100385253C/zh not_active Expired - Fee Related
-
2004
- 2004-12-23 NO NO20045613A patent/NO333954B1/no not_active IP Right Cessation
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2549213C2 (ru) * | 2007-11-14 | 2015-04-20 | ТЕРРАСПАРК ДЖИОСАЙЕНСИЗ, ЭлЭлСи | Обработка сейсмических данных |
Also Published As
| Publication number | Publication date |
|---|---|
| NO20045613L (no) | 2005-02-24 |
| AU2003241526B2 (en) | 2006-03-30 |
| EP1508058A1 (en) | 2005-02-23 |
| NO333954B1 (no) | 2013-10-28 |
| CN1656389A (zh) | 2005-08-17 |
| WO2003100464A1 (en) | 2003-12-04 |
| RU2004137818A (ru) | 2005-05-27 |
| US6636810B1 (en) | 2003-10-21 |
| ATE514106T1 (de) | 2011-07-15 |
| EP1508058B1 (en) | 2011-06-22 |
| AU2003241526A1 (en) | 2003-12-12 |
| CN100385253C (zh) | 2008-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2300123C2 (ru) | Высокоразрешающее преобразование радона для обработки сейсмических данных | |
| US9784868B2 (en) | Method and apparatus for deghosting seismic data | |
| US10281604B2 (en) | Method of removing noise in seismic reverse-time migration | |
| US6763304B2 (en) | Method for processing seismic data to attenuate multiples | |
| US8209126B2 (en) | Wavefront-defined Radon transform | |
| EP3163328A1 (en) | Seismic data least-square migration method and device | |
| US9105075B1 (en) | Enhancing seismic features using an optical filter array | |
| US11175421B2 (en) | Device and method for mitigating cycle-skipping in full waveform inversion | |
| CN107132579B (zh) | 一种保地层结构的地震波衰减补偿方法 | |
| US10393899B2 (en) | Automatic tracking of faults by slope decomposition | |
| MX2013005640A (es) | Metodo y aparato para la eliminacion de fantasma de pre-pila de datos sismicos. | |
| EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
| EP3273274A1 (en) | Device and method for estimating pre-stack wavelet model from seismic gathers | |
| US5684754A (en) | Method and system for correcting seismic traces for normal move-out stretch effects | |
| Ferguson et al. | Planned seismic imaging using explicit one-way operators | |
| WO2005017564A1 (en) | Dip value in seismic images | |
| Geng et al. | Relative time seislet transform | |
| Ventosa et al. | Window length selection for optimum slowness resolution of the local-slant-stack transform | |
| Ursin et al. | Signal and traveltime parameter estimation using singular value decomposition | |
| Sun et al. | Multiple attenuation using λ-f domain high-order and high-resolution Radon transform based on SL0 norm | |
| CN115598704A (zh) | 一种基于最小二乘逆时偏移生成保幅角道集的方法、设备及可读存储介质 | |
| Bianchin et al. | Acoustic impedance estimation from combined harmonic reconstruction and interval velocity | |
| Moore et al. | Stable, efficient, high-resolution Radon transforms | |
| WO2024016572A1 (zh) | 噪声识别方法、装置、设备及存储介质 | |
| Wang et al. | Structure-oriented edge-preserving smoothing based on accurate estimation of orientation and edges |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20070520 |
|
| NF4A | Reinstatement of patent |
Effective date: 20100727 |
|
| MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20130520 |