EA017177B1 - Эффективный способ инверсии геофизических данных - Google Patents
Эффективный способ инверсии геофизических данных Download PDFInfo
- Publication number
- EA017177B1 EA017177B1 EA201070981A EA201070981A EA017177B1 EA 017177 B1 EA017177 B1 EA 017177B1 EA 201070981 A EA201070981 A EA 201070981A EA 201070981 A EA201070981 A EA 201070981A EA 017177 B1 EA017177 B1 EA 017177B1
- Authority
- EA
- Eurasian Patent Office
- Prior art keywords
- data
- model
- sample
- encoded
- source
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 106
- 230000000704 physical effect Effects 0.000 claims abstract description 27
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 8
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 8
- 230000001360 synchronised effect Effects 0.000 claims description 58
- 238000004088 simulation Methods 0.000 claims description 15
- 238000002939 conjugate gradient method Methods 0.000 claims description 4
- 239000004215 Carbon black (E152) Substances 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 2
- 230000004913 activation Effects 0.000 claims 1
- 238000005094 computer simulation Methods 0.000 claims 1
- 238000005553 drilling Methods 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 6
- 230000035508 accumulation Effects 0.000 abstract 1
- 238000009825 accumulation Methods 0.000 abstract 1
- 208000037516 chromosome inversion disease Diseases 0.000 description 86
- 238000004364 calculation method Methods 0.000 description 24
- 238000004422 calculation algorithm Methods 0.000 description 22
- 238000012360 testing method Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 238000004880 explosion Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000005305 interferometry Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000007670 refining Methods 0.000 description 2
- 241000899793 Hypsophrys nicaraguensis Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000004613 tight binding model Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
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
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
-
- 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/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Способ для эффективной инверсии измеренных геофизических данных из области недр для разведки углеводородов. Выборки измеренных данных (40) кодируют (60) с использованием набора неэквивалентных кодирующих функций (30). Затем все записи данных в каждой кодированной выборке, соответствующей одному приемнику, суммируют (60) с повторением для каждого приемника для формирования синхронной кодированной выборки (80). В способе используют итерационную локальную оптимизацию функции стоимости для инвертирования кодированных выборок данных синхронного источника. Метод сопряженных функций применяют для расчета градиентов функции стоимости, необходимых для процесса локальной оптимизации (100). Инвертированные данные дают модель физических свойств области недр, которая после итерационного уточнения может показывать наличие месторождений углеводородов.
Description
Изобретение в целом относится к области геофизической разведки и, в частности, к обработке геофизических данных. Конкретно, изобретение представляет собой способ инверсии данных, получаемых из разных геофизических источников, таких как сейсмические источники, включающий геофизическое моделирование, в котором за одно выполнение моделирования вычисляют данные из многих одновременно активных геофизических источников.
Предпосылки изобретения
При геофизической инверсии предпринимается попытка найти модель свойств недр, которая оптимально объясняет экспериментальные данные и удовлетворяет геологическим и геофизическим ограничениям. Существует большое количество хорошо известных способов геофизической инверсии. Эти хорошо известные способы делятся на две группы: итеративные инверсии и неитеративные инверсии.
Неитеративная инверсия представляет собой инверсию, выполняемую посредством принятия некоторой простой фоновой модели и уточнения модели на основе входных данных. В данном способе не используют уточненную модель в качестве входных данных для следующего этапа инверсии. В случае применения к сейсмическим данным эти способы обычно называют визуализацией, миграцией, дифракционной томографией или инверсией Борна.
Итеративная инверсия представляет собой инверсию, включающую в себя повторяющиеся уточнения модели свойств недр, при которых найденная модель будет удовлетворительно объяснять экспериментальные данные. Если инверсия сходится, то конечная модель будет лучше объяснять экспериментальные данные и будет более точно аппроксимировать фактические свойства недр. Итеративная инверсия обычно дает более точную модель, чем неитеративная инверсия, однако, она требует значительно больших вычислительных затрат.
Итеративная инверсия, как правило, является предпочтительней, чем неитеративная инверсия, так как она приводит к более точной модели свойств недр. К сожалению, итеративная инверсия является настолько затратной по вычислениям, что она оказывается неприменимой для решения многих представляющих интерес задач. Такая затратность по вычислениям является следствием того, что все способы инверсии требуют большого количества вычислительно интенсивных моделирований. Время вычисления любого индивидуального моделирования пропорционально количеству инвертируемых источников, а геофизические данные обычно содержат большое количество источников. Задача осложняется для итеративной инверсии, так как количество моделирований, которое необходимо вычислить, пропорционально количеству итераций инверсии, при этом, как правило, требуется от нескольких сотен до нескольких тысяч итераций.
Наиболее часто применяемым в геофизике способом итеративной инверсии является оптимизация функции стоимости. Оптимизация функции стоимости предполагает в зависимости от модели М итеративную минимизацию или максимизацию значения функции стоимости 8(М), которая является мерой несоответствия между расчетными и экспериментальными данными (иногда ее называют целевой функцией), где расчетные данные моделируются на компьютере с использованием текущей модели геофизических свойств и физических свойств, обуславливающих распространение сигнала источника в среде описываемой данной моделью геофизических свойств. Вычисления моделирования могут быть выполнены любым из многих численных методов, включая, но не ограничиваясь, способ конечных разностей, метод конечных элементов или метод трассировки лучей. Вычисления моделирования могут выполняться как в частотной, так и во временной области.
Методы оптимизации функции стоимости являются либо локальными, либо глобальными [3]. Глобальные методы просто предполагают вычисление функции стоимости 8(М) для совокупности моделей {Мь М2, М3, ...} и выбор набора из одной или нескольких моделей из совокупности, которые приближенно минимизируют 8(М). Если требуется дальнейшее улучшение, то этот новый выбранный набор моделей можно использовать в качестве базиса для генерации новой совокупности моделей, которые опять можно проверить относительно функции стоимости 8(М). Для глобальных способов каждую модель в тестовой совокупности можно рассматривать как итерацию или, на более высоком уровне, каждый проверяемый набор совокупностей можно рассматривать как итерацию. Хорошо известные методы глобальной инверсии включают в себя метод Монте-Карло, метод модельной закалки, генетические и эволюционные алгоритмы.
К сожалению, глобальные методы оптимизации обычно сходятся крайне медленно и поэтому большинство геофизических инверсий основано на локальной оптимизации функции стоимости. Локальная оптимизация функции стоимости кратко сформулирована в алгоритме 1.
- 1 017177
1. выбор начальной модели,
2. вычисление градиента функции стоимости 3 (М) с учетом параметров, описывающих модель,
3. поиск уточненной модели, являющейся возмущением начальной модели, которая лучше объясняет экспериментальные данные, в направлении отрицательного градиента.
Алгоритм 1 - алгоритм выполнения локальной оптимизации функции стоимости.
Эту процедуру повторяют с использованием новой уточненной модели в качестве начальной для другого градиентного поиска. Процесс продолжается до тех пор, пока не обнаружится, что уточненная модель удовлетворительно объясняет экспериментальные данные. Часто используемые способы локальной инверсии функции стоимости включают в себя метод градиентного поиска, метод сопряженных градиентов и метод Ньютона. Далее эта общая информация будет раскрыта несколько подробнее.
Локальная оптимизация функции стоимости сейсмических данных в акустическом приближении является распространенной задачей геофизической инверсии и, как правило, характерна для других типов геофизической инверсии. При инвертировании сейсмических данных в акустическом приближении функцию стоимости можно записать в виде
М, * , , ) = £ Σ Σ IV(^е {М, г, ί, ^)- ψοίχ (г, {, )) (1) «=1 г=1 1=1 где 8 - функция стоимости,
М - вектор N параметров, (т1, т2, ..., пк). описывающих модель недр, д - индекс выборки, \\д - функция источника для выборки д, являющаяся функцией пространственных координат и времени, для точечного источника является δ-функцией пространственных координат,
Ν6 - количество выборок, гд - индекс приемника в выборке,
N - количество приемников в выборке, ΐ - индекс временной выборки в трассе, N - количество временных выборок,
- функция критерия минимизации (мы обычно выбираем функцию ^(х)=х2, которая соответствует критерию наименьших квадратов (Ь2)), ψ03ι0 - вычисленные данные сейсмических давлений, ψ,,ι„ - измеренные данные сейсмических давлений.
В качестве выборок сейсмических данных в уравнении 1 может выступать любой тип выборки, который можно смоделировать за один прогон программы прямого сейсмического моделирования. Как правило, выборки соответствуют сейсмическому взрыву, хотя взрыв может быть более общим, чем точечный источник. Для точечных источников индекс д выборки соответствует положению индивидуальных точечных источников. Для источников плоских волн, индекс д будет соответствовать различным направлениям распространения плоской волны. Эти данные обобщенного источника, ψο|8, могут или быть получены в полевых условиях, или быть синтезированными из данных, полученных с использованием точечных источников. С другой стороны, вычисленные данные ψ03ι0 обычно можно вычислить напрямую с использованием функции обобщенного источника в прямом моделировании. Для многих видов прямого моделирования, включая моделирование методом конечных разностей, время вычисления, необходимое для обобщенного источника, примерно равно времени вычисления, необходимого для точечного источника.
Уравнение 1 можно упростить до
Я,
5(М) = £)Г(4М,>Р8)) (2) й=1 где сумма по приемникам и временным выборкам теперь выражена не явно, и
Целью инверсии посредством оптимизации функции стоимости является попытка уточнить модель М таким образом, чтобы значение 8(М) являлось минимальным. Это можно выполнить посредством локальной оптимизации функции стоимости, которая уточняет данную модель М(к) следующим образом: = М® - а^м8(м) <4) где к - номер итерации, а - скалярный размер уточнения модели и νΜ8(Μ) - градиент функции несоот
- 2 017177 ветствия, взятый с учетом параметров модели. Возбуждения модели или значения, которыми уточняется модель, вычисляются умножением градиента целевой функции на длину а шага, которая должна циклически вычисляться.
Из уравнения 2 можно получить следующее уравнение для градиента функции стоимости:
<5>
Следовательно, для того чтобы вычислить градиент функции стоимости, необходимо отдельно вычислить градиент вклада в целевую функцию от каждой выборки, затем суммировать эти вклады. Следовательно, вычислительное усилие для вычисления νΜ8(Μ) представляет собой Ν6 раз вычислительное усилие, необходимое для определения вклада одной выборки в градиент. В геофизических задачах Ν6 обычно соответствует количеству геофизических источников (каждое местоположение аппаратовисточников рассматривается как отдельный источник) и является числом порядка от 10000 до 100000, существенно увеличивая стоимость расчета νΜ8(Μ).
Следует отметить, что вычисление νΜ^(δ) требует вычисления производной от ^(8) по каждому из N параметров модели ш1. Так как в геофизических задачах N обычно очень велико (количество различных параметров, умноженное на количество ячеек модели, в которых параметрам должны быть присвоены значения, обычно более миллиона), то это вычисление может быть крайне длительным, если его необходимо выполнять для каждого индивидуального параметра модели.
Существует потребность в более эффективном способе вычисления градиента функции стоимости без существенного снижения точности локальной оптимизации функции стоимости. Настоящее изобретение удовлетворяет эту потребность.
Сущность изобретения
В одном из вариантов осуществления изобретение представляет собой осуществляемый на компьютере способ инверсии измеренных геофизических данных для определения модели физических свойств для области недр, содержащий:
(a) получение группы из двух или более кодированных выборок измеренных геофизических данных, причем каждая выборка ассоциирована с одним обобщенным источником или, используя взаимность источника и приемника, с одним приемником, и причем каждая выборка кодирована разными кодирующими функциями, выбираемыми из набора неэквивалентных кодирующих функций;
(b) суммирование кодированных выборок в группе посредством суммирования всех записей данных в каждой выборке, соответствующих одному приемнику (или источнику, если используется взаимность), и повторение для каждого другого приемника, приводящее к синхронной кодированной выборке;
(c) принятие модели физических свойств области недр, предоставление указанной моделью значений по меньшей мере одного физического свойства в местах по всей области недр; и (б) инвертирование измеренных геофизических данных по одной синхронной кодированной выборке за раз, используя принятую модель физических свойств в качестве начальной модели, и итеративное уточнение указанной модели для минимизирования функции стоимости, измеряющей степень несоответствия между моделируемыми данными и измеренными геофизическими данными, для формирования уточненной модели физических свойств, где корректирование модели осуществляют с использованием градиента функции стоимости в отношении по меньшей мере одного параметра модели, при этом градиент вычисляют из интегрирования по времени произведения данных кодированных синхронных источников, моделируемых вперед по времени, и данных кодированных синхронных источников, моделируемых назад по времени.
Уточненная модель физических свойств обычно загружается с компьютера или сохраняется в запоминающее устройство компьютера и, в конечном счете, используется для заключения о наличии или отсутствии месторождений углеводородов в области недр.
Технический результат, на достижение которого направлено данное изобретение, состоит в уменьшении времени вычислений, необходимого для оценки параметров физических свойств посредством инверсии данных.
Краткое описание чертежей
Настоящее изобретение и его преимущества станут более понятными при обращении к нижеследующему подробному описанию и к прилагаемым чертежам.
На фиг. 1 показана блок-схема последовательности операций, отражающая основные этапы в одном из вариантов осуществления настоящего изобретения;
на фиг. 2 - блок-схема последовательности операций, отражающая основные этапы в альтернативном варианте осуществления настоящего изобретения;
на фиг. 3 - основная скоростная модель, используемая для формирования измеренных данных, использованная в примере, представленном в настоящем документе, с масштабной полосой шкалы уровней серого, показывающей скорость в м/с;
на фиг. 4 - выборка одного источника для данных последовательного источника из модели, показанной на фиг. 3;
- 3 017177 на фиг. 5 - данные синхронного кодированного источника, сформированные на этапе 260 фиг. 2;
на фиг. 6 - начальная модель (210 на фиг. 2) для примера инверсии;
на фиг. 7 - соответствие модели нормы Ь2 в зависимости от итерации инверсий последовательного источника и кодированного синхронного источника в примере;
на фиг. 8 - результаты двенадцатой итерации инверсии последовательного источника из примера и на фиг. 9 - результаты двадцать первой итерации инверсии кодированного синхронного источника в примере.
Изобретение будет описано по отношению к его предпочтительным вариантам осуществления. Однако в тех случаях, когда следующее ниже подробное описание является характерным для определенного варианта осуществления или определенного использования изобретения, оно предназначено только для иллюстрации, и его не следует трактовать как ограничивающее объем изобретения. Наоборот, оно предназначено для включения всех альтернатив, модификаций и эквивалентов, которые можно включить в объем изобретения, как определено прилагаемой формулой изобретения.
Подробное описание иллюстративных вариантов осуществления
Настоящее изобретение воплощает в себе понимание того, что вычислительное время, требуемое для локальной оптимизации функции стоимости, можно значительно сократить за счет использования способа, называемого методом сопряженных уравнений, для расчета градиентов функции стоимости, необходимых для процесса локальной оптимизации, и потом применяя этот подход для инверсии кодированных выборок данных синхронного источника, тем самым вычисляя градиент функции стоимости кодированного синхронного источника.
Метод сопряженных уравнений можно использовать для эффективного выполнения (см. описание, следующее за уравнением 5) вычисления производной Α(δ) относительно каждого из N параметров модели т1, для всех параметров модели сразу [1]. Метод сопряженных уравнений для целевой функции наименьших квадратов и сеточной параметризации модели можно свести к следующему алгоритму:
1. | Вычисление | прямого моделирования данных | с | использованием |
текущей модели и сигнатуры выборки как | источника для | |||
Получения 1|)са1с (М(к), ν3) , | ||||
2. | Вычитание | экспериментальных данных | из | моделируемых |
данных, что | дает б(М(к>,ид) , | |||
3. | Вычисление | обратного моделирования | (т. | е. назад во |
времени) с | использованием δ{Μ(Ιι),νσ) в качестве источника. | |||
дающий ФаЩои-Л !М1?·’,и3) , | ||||
4. | Вычисление | интеграла по времени | от | произведения |
^са!с ΪΜ* *, Ид, | и Фаазохис (М'к| ,иэ) для получения | νΜΗ (δ (М, Ид) ) . |
Алгоритм 2 - алгоритм вычисления градиента функции стоимости наименьших квадратов в сеточной модели с использованием метода сопряженных уравнений.
Хотя вычисление градиентов с использованием метода сопряженных уравнений является эффективным по сравнению с другими методами, оно все еще является очень дорогостоящим. В частности, метод сопряженных уравнений требует проведения двух моделирований, одно вперед по времени и одно назад по времени, и для геофизических задач эти моделирования обычно требуют большого объема вычислений. Также, как обсуждалось выше, эти вычисления методом сопряженных уравнений следует выполнять индивидуально для каждой выборки измеренных данных, увеличивая вычислительную стоимость в Ν6 раз.
Как указано в ссылке [4], вычислительную стоимость всех категорий инверсии можно уменьшить посредством инверсии данных от комбинаций источников вместо инверсии источников по отдельности. Это можно назвать инверсией синхронного источника. Известно несколько типов комбинаций источников, включая когерентно суммирующие близко расположенные источники для получения эффективного источника, который производит волновой фронт желаемой формы (например, плоскую волну); суммирующие далеко расположенные источники, или полное или частичное суммирование данных перед инверсией.
Уменьшение вычислительной стоимости, полученное за счет инвертирования комбинированных источников, по крайней мере, частично компенсируется тем фактом, что инверсия комбинированных данных обычно дает менее точную инвертированную модель. Эта потеря точности происходит за счет того факта, что при суммировании отдельных источников теряется информация, и, таким образом, суммированные данные не так строго ограничивают инвертированную модель, как несуммированные данные. Эту потерю информации при суммировании можно минимизировать посредством кодирования каждой точечной записи перед суммированием. Кодирование перед комбинацией сохраняет значительно больше информации в данных синхронных источников и, таким образом, лучше ограничивает инвертирование [4]. Кодирование также допускает комбинацию близко расположенных источников, таким обра
- 4 017177 зом, позволяя комбинирование большего количества источников для данной вычислительной области. В этом способе можно использовать различные схемы кодирования, включая кодирование временного сдвига и кодирование случайной фазы. Ниже представлены краткие обзоры различных опубликованных геофизических способов синхронного источника как кодированного, так и некодированного.
Уап Мапеп [5] предлагает использовать способ сейсмической интерферометрии для ускорения прямого моделирования. Способ сейсмической интерферометрии осуществляется посредством размещения источников повсюду на границе интересующей области. Эти источники моделируют индивидуально и для всех точек, в которых необходима функция Грина, записывают волновое поле. Тогда функцию Грина между любыми двумя записанными точками можно вычислить путем кросс-коррелирования трасс, полученных в двух записанных точках, и суммирования по всем граничным источникам. Если инвертируемые данные содержат большое количество источников и приемников, находящихся внутри интересующей области (в отличие от нахождения одного или другого на границе), тогда этот способ оказывается очень эффективным для вычисления необходимых функций Грина. Однако в случае сейсмических данных довольно редко источник и приемник данных, подлежащие инверсии, находятся внутри интересующей области. Таким образом, это улучшение имеет очень ограниченную применимость для задачи инверсии сейсмических данных.
Веткйои! [6] и Ζ1ιαη§ [7] предполагают, что инверсию в большинстве случаев можно улучшить посредством инвертирования некодированных синхронных источников, которые когерентно суммируются для получения некоторого желаемого волнового фронта внутри определенной области недр. Например, данные точечных источников можно суммировать с временными сдвигами, являющимися линейной функцией положения источников, для получения нисходящей плоской волны под некоторым углом относительно поверхности. Эту методику можно применять ко всем категориям инверсии. Проблема этого способа состоит в том, что когерентное суммирование выборок трасс источников неизбежно уменьшает количество информации в данных. Так, например, суммирование для получения плоской волны удаляет всю информацию в сейсмических данных, относящуюся к зависимости времени распространения от расстояния источник - приемник. Эта информация является критической для уточнения медленно меняющейся модели фоновой скорости, и таким образом способ Беркхаута не является хорошо ограниченным. Для преодоления этой проблемы можно инвертировать многие другие когерентные суммы данных (например, множество плоских волн с разными направлениями распространения), но тогда эффективность теряется, поскольку цена инверсии пропорциональна количеству различных инвертируемых сумм. Авторы изобретения называют такие когерентно суммируемые источники обобщенными источниками. Таким образом, обобщенный источник может быть как точечным источником, так и суммой точечных источников, создающих волновой фронт необходимой формы.
Уап КТе1 [8] предлагает инверсию при помощи некодированного суммирования или частичного суммирования (в зависимости от расстояния источник - приемник) входных сейсмических данных с последующим определением функции стоимости по отношению к этим суммированным данным, которую будут оптимизировать. Таким образом, эта публикация предлагает улучшение инверсии, основанной на функции стоимости, с использованием некодированных синхронных источников. Как и в способе инверсии синхронного источника Беркхаута [6], суммирование, предлагаемое в этом способе, уменьшает количество информации в инвертируемых данных, и таким образом, инверсия является менее ограниченной, чем она была бы при использовании оригинальных данных.
Мога [9] предлагает инвертирование данных, являющихся суммой далеко расположенных источников. Таким образом, в этой публикации предложено улучшение эффективности инверсии с использованием моделирования некодированных синхронных источников. Преимуществом суммирования далеко расположенных источников является сохранение значительно большего количества информации, чем при когерентном суммировании, предложенном Беркхаутом. Однако суммирование далеко расположенных источников подразумевает, что апертуру (инвертируемая область модели), которую необходимо использовать в инверсии, нужно увеличить для размещения всех далеко расположенных источников. Так как вычислительное время пропорционально площади этой апертуры, способ Мора не дает такого прироста эффективности, которой можно достичь при суммировании источников, расположенных близко друг к другу.
ОЬет [10] предлагает ускорение сейсмической миграции, частного случая неитеративной инверсии, посредством использования синхронных кодированных источников. После тестирования различных способов кодирования Обер обнаружил, что результирующие мигрированные изображения имели существенно сниженное соотношение сигнал-шум за счет того факта, что широкополосные кодирующие функции неизбежно только приблизительно ортогональны. Таким образом, при суммировании более 16 записей качество инверсии было неудовлетворительным. Так как неитеративная инверсия сначала не очень затратна, и так как необходима инверсия с высоким соотношением сигнал-шум, то этот способ не является широкоприменимым в геофизической отрасли.
1ке11е [11] предлагает способ для быстрого прямого моделирования при помощи синхронно моделируемых точечных источников, которые активируются (в моделировании) на различных временных интервалах. Также обсужден способ декодирования этих сдвинутых по времени моделированных данных
- 5 017177 синхронного источника обратно в раздельные моделирования, которые были бы получены из отдельных точечных источников. Эти декодированные данные далее можно использовать как часть любой обычной процедуры инверсии. Проблема способа 1ке11е состоит в том, что предложенный способ декодирования будет производить разрозненные данные, имеющие уровень шума, пропорциональный разности между данными от соседних источников. Этот шум будет существенным для моделей недр, которые не являются постоянными по латерали, например от моделей, содержащих горизонтальные отражающие границы. Следовательно, шум будет увеличиваться пропорционально числу синхронных источников. Вследствие этих трудностей подход синхронного источника 1ке11е может привести к неприемлемым уровням шума при использовании в инвертировании недр, не являющихся постоянными в поперечном направлении. КтеЬк е! а1. [4] показывает, что функцию стоимости кодированного синхронного источника можно вычислить более эффективно, чем традиционные функции стоимости при обеспечении точных инверсий. Функцию стоимости синхронного источника определяют здесь как (сравнить с уравнением 2 выше)
0=1 I К®о )) >
где суммирование по приемникам и временным выборкам подразумевают как в уравнении 2, и
Σ-ΣΣ определяет сумму по выборкам при помощи подгрупп выборок,
- функция стоимости для данных синхронных источников,
С - группы синхронных обобщенных источников,
Νι; - количество групп, сд - функции времени, которые сворачиваются (®) с сигнатурой источника каждой выборки для кодирования выборок, эти кодирующие функции выбираются разными, то есть неэквивалентными (то есть без ограничений, приблизительно ортогональными), для каждого индекса д выборки (то есть различные реализации функций случайной фазы).
Внешнее суммирование в уравнении 6 ведется по группам синхронных обобщенных источников, соответствующих типу выборки (то есть точечных источников для выборок общего взрыва). Внутреннее суммирование по д выполняется по выборкам, которые сгруппированы для синхронного вычисления. Для некоторых способов прямого моделирования, таких как моделирование конечных разностей, вычисление прямой модели для суммированных обобщенных источников (внутренняя сумма по де С) можно выполнять за то же время, что и вычисления для одного источника. Таким образом, как показано в КтеЬк е! а1. [4], 8(М, Тсд0\\д) можно вычислять очень эффективно, используя такой подход, как алгоритм 3
Моделирование фОа1с (М,2сд®вдд) с использованием одного прогона моделирующей программы, используя 2сэ®ид в качестве источника.
Сворачивание каждой выборки измеренных данных с сд кодирующими функциями, и последующее суммирование результирующих кодированных выборок (то есть
2сд®фоЬБ (ν3)},
Вычитание результата этапа 2 из результата этапа 1.
Алгоритм 3 - алгоритм вычисления кодированной функции стоимости синхронного источника.
Кроме того, как показано в КтеЬк е! а1. [4], по этому алгоритму можно вычислить 8мт(М) в Ν6/Ν; раз быстрее, чем 8(М) из уравнения 2.
Приведенный выше алгоритм или его варианты, которые очевидны специалистам в технической области, можно использовать для очень эффективного вычисления 8мт(М), и, таким образом, они могут значительно ускорять инверсии, которые требуют только вычисления функции стоимости, такие как способы глобальной инверсии. Фактически, посредством приведенного выше алгоритма можно также ускорить вычисление градиента функции стоимости для многих градиентных методов оценки. Например, вычисление градиента методом конечных разностей предполагает вычисление функции стоимости для модели М и близкой модели М+5(М), в которой один из параметров возмущен на малую величину.
Однако чрезвычайно большое число параметров модели, характерное для геофизической инверсии, означает, что вычисление градиента 8мт(М) будет по-прежнему очень дорогим, если будут использоваться традиционные методы вычисления градиента, такие как метод конечных разностей. В настоящем изобретении эта проблема решена за счет вычисления градиента 88ИП(М) с использованием метода сопряженных уравнений. Таким образом, цель настоящего изобретения заключается в вычислении значения
- 6 017177
νΜ8,,ΙΙΙΙ(Μ). которое определено как (л/)=2νΜψ\ <? м, £ се ® <3=1 ( I εεσ
В настоящем изобретении. величины νΜ^(δ(Μ.Σε6®^6)) можно очень эффективно вычислять с использованием способа типа сопряженных уравнений. такого как изложенный в алгоритме 4
Алгоритм 4 - алгоритм для вычисления градиента функции стоимости кодированного синхронного источника с использованием способа типа сопряженных уравнений.
В практическом осуществлении изобретение. как правило. можно использовать для выполнения инверсии локальной функции стоимости согласно блок-схеме последовательности операций одного из вариантов осуществления изобретения. представленной на фиг. 1. Как и в традиционной итеративой инверсии. процесс начинается с выбора начальной модели 10. которая является аппроксимацией фактических физических свойств недр (уточненная модель 20 заменяет начальную модель в последующих итерациях инверсии). На следующем этапе инверсии (30) строятся кодирующие функции; их используют для кодирования измеренных данных 40 на этапе 60 и сигнатуры источников 50 на этапе 70. Как описано ниже и в ссылке [4]. на этапе 30 можно получить кодирующие функции многих различных видов. На этапах 60 и 70 получают кодированные измеренные данные (суммированные по местоположениям источников для каждого приемника) 80 и кодированные сигнатуры источника 90. где и те. и другие кодированы одним и тем набором кодирующих функций. полученных на этапе 30. где способ для читателей. не являющихся специалистами данной области. подробнее описан в работе [4]. На этапе 100 с использованием уравнения 7 и алгоритма 4. описанных выше. или их тривиальных вариаций. вычисляют градиент. На этапе 110 на основании градиента. вычисленного на этапе 100. начальную модель уточняют. как правило. посредством добавления или вычитания масштабированного варианта градиента к начальной модели. Этап 120 является проверкой схождения инверсии. Его можно осуществлять посредством вычисления функции стоимости кодированного синхронного источника для уточненной модели с использованием уравнения 6 и алгоритма 3. и сравнения значения функции стоимости со значением из предыдущей итерации. Если сделан вывод. что инверсия не сошлась. то процесс повторяют с этапа 30. используя уточненную модель 20 из этапа 110. Когда проверка 120 показывает. что итерация сошлась. или выполняется другое условие остановки. то процесс завершают (130).
В уравнениях 6 и 7 можно использовать множество типов кодирующих функций сд. включая. но не ограничиваясь следующими:
линейное. случайное. кодирование линейной частотной модуляцией и модифицированное частотнозависимое фазовое кодирование линейной частотной модуляцией. как представлено в Вошего с1 а1. [12];
частотно-независимое фазовое кодирование. как представлено в Лид е1 а1. [13]; кодирование случайного сдвига по времени;
мультиплексирование с разделением частот (ΕΌΜΑ). мультиплексирование с разделением времени (ΤΌΜΑ) и мультиплексирование с кодовым разделением (СЭМА). используемые в телекоммуникациях.
- 7 017177
Одни из этих способов кодирования работают лучше других в зависимости от области применения, также некоторые из них можно комбинировать. В частности, хорошие результаты можно получать с использованием частотно-зависимого кодирования случайных фаз и также посредством комбинирования частотно-независимого кодирования близкорасположенных источников и частотно-зависимого кодирования случайных фаз для более далеко расположенных источников. Показатель сравнительной ценности различных кодирований можно получать, выполняя тестовые инверсии с каждым набором кодирующих функций для определения того набора, при котором инверсия сойдется быстрее.
Следует отметить, что способ синхронного кодированного источника можно использовать для различных типов функции стоимости инверсии. В частности, его можно использовать для функции стоимости, основанной на нормах, отличающихся от нормы Ь2, которая обсуждалась выше. Его также можно использовать с более сложными функциями стоимости, чем функция, представленная в уравнении 2, включая регуляризированные функции стоимости. Наконец, способ синхронного кодированного источника можно использовать с любым видом глобального или локального метода инверсии функции стоимости, включая метод Монте-Карло, алгоритм модельной закалки, генетический алгоритм, эволюционный алгоритм, линейный градиентный поиск, метод сопряженных градиентов и метод Ньютона.
Способ по настоящему изобретению также можно использовать в сочетании с различными видами способов обобщенного источника, такими как способ, предложенный Вегкйои! [6]. В этом случае вместо кодирования сигнатур выборки различных точечных источников следует кодировать сигнатуру различных синтезированных плоских волн.
Некоторые разновидности представленного выше варианта осуществления включают следующее.
Кодирующие функции сд могут изменяться для каждой итерации инверсии. По крайней мере, в некоторых случаях это приводит к более быстрому схождению инверсии.
В некоторых случаях (например, когда выборка источников плотнее, чем выборка приемников) выгодным может являться использование взаимности для рассматривания фактических приемников в качестве вычислительных источников и кодирование приемников вместо источников.
Настоящее изобретение не ограничивается однокомпонентными точечными приемниками. Например, приемники могут являться массивами приемников, либо они могут являться многокомпонентными приемниками.
Настоящий способ можно улучшить посредством оптимизации кодирования для достижения максимально высокого качества инверсии. Например, кодирующие функции можно оптимизировать для уменьшения количества локальных минимумов в функции стоимости. Кодирующие функции можно оптимизировать как посредством ручной проверки тестов, выполненных с использованием различных кодирующих функций, так и с использованием автоматизированных процедур оптимизации.
Сбор данных синхронных кодированных источников может привести к значительному снижению затрат по сбору геофизических данных.
Для морских сейсмических исследований может оказаться очень эффективным сбор данных кодированного источника от одновременно действующих морских вибраторов, которые работают постоянно во время движения.
Можно использовать другие определения функции стоимости, включая использование различных норм (например, норма Ь1 (абсолютное значение) вместо нормы Ь2) и дополнительных членов для упорядочивания и стабилизации инверсии (например, члены, которые будут штрафовать модели, которые не являются гладкими, или модели, не являющиеся разреженными).
Фиг. 2 представляет собой блок-схему последовательности операций альтернативного варианта осуществления настоящего способа по изобретению для инверсии с использованием локальной функции стоимости. Эта блок-схема последовательности операций, по существу, аналогична блок-схеме, представленной на фиг. 1 (условные обозначения соответствующих полей увеличены на 200), за исключением следующих изменений:
на этапе 230 кодирующие функции между итерациями могут изменяться;
добавлен этап 305. На этапе 305 градиент с этапа 300 согласовывается с использованием любого хорошо известного метода, и для лучшего схождения согласованный градиент используется для вычисления сопряженного градиента. Широко известно, что в схемах (локальной) оптимизации сходимость алгоритма можно улучшать за счет выбора в качестве хорошего направления (градиента) линейной комбинации предыдущего и текущего хороших направлений. Этот метод известен как метод сопряженных градиентов. Согласование (или предсогласование) представляет собой классически известную стратегию улучшения сходимости;
на этапе 310 используют линейный поиск в направлении сопряженного градиента для уточнения модели. Линейный поиск является частью стандартных стратегий оптимизации и состоит из старта в текущей точке, выбора направления для вычисления значений функции стоимости, движения вдоль этого выбранного направления и перерасчета функции стоимости в каждой новой точке, и выбора точки, дающей минимальное значение используемой функции стоимости. Направление (в данном случае) является градиентом или сопряженным градиентом.
Дополнительно в данном конкретном иллюстративном варианте осуществления настоящего спосо
- 8 017177 ба по изобретению сделан следующий выбор во время инверсии.
1. Входные сейсмические выборки измеренных данных 240 являются общими выборками точечных источников.
2. В качестве кодирующих сигнатур 230 выбраны сигнатуры случайных фаз из работы Вотсго с1 а1. Такие сигнатуры можно получать просто посредством создания последовательности, состоящей из временных выборок, которая является равномерной псевдослучайной последовательностью.
3. На этапах 260 и 270 выборки и сигнатуры кодированы посредством сворачивания каждой трассы в выборке с кодирующей сигнатурой этой выборки.
4. На этапах 300 и 320 моделирование выполняется посредством кода сейсмического моделирования конечными разностями в пространственно-временной области.
5. На этапах 300 и 320 функция стоимости вычисляется с использованием нормы Ь2.
Пример.
Фиг. 3-9 представляют комплексный пример инвертирования акустических сейсмических данных постоянной плотности с использованием настоящего изобретения и для сравнения с традиционным способом последовательного источника.
Для этого примера использован вариант осуществления, описанный на фиг. 2, за исключением того, что исключен этап 305. Включение этапа 305 улучшает результаты по сравнению с этим примером, так что это является консервативным примером преимуществ данного изобретения.
Фиг. 3 представляет собой базовую скоростную модель, использованную для формирования измеренных данных (фиг. 4) для данного примера. Успешная инверсия измеренных данных, моделируемая данной моделью, должна выглядеть как ограниченная по полосе версия фиг. 3.
На фиг. 4 показана выборка одного источника измеренных данных, моделируемых на основании модели, приведенной на фиг. 3. Всего смоделировано 128 источников (из которых на фиг. 4 представлен первый) с источниками, расположенными возле поверхности модели, и источниками, равномерно распределенными по модели. Также равномерно по модели распределены 1024 сейсмических приемника. Тот же набор приемников использовали для каждого расположения источника (то есть приемники не перемещались с источниками).
На фиг. 5 показаны данные синхронного кодированного источника, полученные на этапе 260 фиг. 2. Кодирующей функцией являлась трехсекундная функция случайных фаз. Все 128 источников кодированы различными функциями случайных фаз, после чего суммированы для получения данного синхронного кодированного источника. Следует отметить, что данные кодированного синхронного источника на фиг. 5 выглядят как случайный шум по сравнению с фиг. 4.
На фиг. 6 показана модель, использованная в качестве начальной модели (210 на фиг. 2) в данном примере инверсии. Данная модель получена посредством интенсивного сглаживания базовой модели, показанной на фиг. 3. Сглаживание выбрано достаточно большим для удаления всех отражений и мелкомасштабных аномалий скорости, в то же время сохраняя гладкий фон модели.
Затем данные на фиг. 4 и 5 инвертировали. Понятно, что настоящее изобретение применяли для инверсии данных синхронного кодированного источника из фиг. 5. Следует отметить, что так как в данных последовательного источника существует 128 выборок источников, то инверсия этих данных должна требовать приблизительно в 128 раз больше вычислительных затрат за итерацию, чем это необходимо для инвертирования данных кодированного синхронного источника. Это соотношение в 128 раз связано с тем фактом, что градиент необходимо вычислять отдельно для каждого последовательного источника, и потом эти градиенты суммируют для получения полного градиента последовательного источника. С другой стороны, градиент для данных кодированного синхронного источника вычисляют только из одной выборки.
На фиг. 7 показана зависимость подгонки модели от числа итерации для двух инверсий. Инверсия последовательного источника (черные квадраты) сходится примерно в два раза быстрее, чем инверсия синхронного кодированного источника (белые квадраты). (Подгонка модели для итерации 21 инверсии кодированного синхронного источника примерно равна подгонке модели для итерации 12 инверсии последовательного источника). Однако так как каждая итерация инверсии синхронного кодированного источника использует только 1/128 часть компьютерных затрат, то инверсия кодированного синхронного источника является в 64 раза более эффективной, чем инверсия последовательного источника.
На фиг. 8 и 9 показаны результаты обеих инверсий. Для отображения выбраны двенадцатая инверсия последовательного источника и двадцать первая инверсия кодированного синхронного источника, потому что у них приблизительно одинаковые подгонки модели. Следует заметить, что результаты обеих инверсий выглядят приблизительно как ограниченная по полосе версия базовой модели на фиг. 3. Также следует обратить внимание на то, что инверсии кодированного синхронного источника очень похожи на инверсии последовательного источника, за исключением добавления небольшого количества шума. Этот шум можно легко удалить или посредством обработки изображения результата инверсии, или посредством применения согласования (этап 305 на фиг. 2) при инверсии.
Описанная выше заявка относится к частным вариантам осуществления настоящего изобретения с целью его иллюстрации. Однако для специалиста в данной области очевидно, что для вариантов осуще
- 9 017177 ствления, описываемых в данном документе, возможно множество модификаций и вариаций. Полагают, что все такие модификации и вариации входят в объем настоящего изобретения, как определено в прилагаемой формуле изобретения. Специалисты в данной области легко поймут, что в предпочтительных вариантах осуществления изобретения, по меньшей мере, некоторые этапы настоящего способа проводятся на компьютере, то есть изобретение осуществляется на компьютере. В таких случаях результирующую уточненную модель физических свойств можно либо загрузить, либо сохранить на запоминающим устройстве компьютера.
Claims (19)
1. Осуществляемый на компьютере способ инверсии измеренных геофизических данных для определения модели физических свойств для области недр, содержащий:
(a) получение группы из двух или более кодированных выборок измеренных геофизических данных, причем каждая выборка ассоциирована с одним обобщенным источником или с одним приемником, и причем каждая выборка кодирована разными кодирующими функциями, выбираемыми из набора неэквивалентных кодирующих функций;
(b) суммирование кодированных выборок в группе посредством суммирования всех записей данных в каждой выборке, соответствующих одному приемнику (или источнику, если используется взаимность), и повторение для каждого другого приемника, приводящее к синхронной кодированной выборке;
(c) принятие модели физических свойств области недр, предоставление указанной моделью значений по меньшей мере одного физического свойства в местах по всей области недр и (б) инвертирование измеренных геофизических данных по одной синхронной кодированной выборке за раз, используя принятую модель физических свойств в качестве начальной модели, и итеративное уточнение указанной модели для минимизирования функции стоимости, измеряющей степень несоответствия между моделируемыми данными и измеренными геофизическими данными, для формирования уточненной модели физических свойств, где корректирование модели осуществляют с использованием градиента функции стоимости в отношении по меньшей мере одного параметра модели, при этом градиент вычисляют из интегрирования по времени произведения данных кодированных синхронных источников, моделируемых вперед по времени, и данных кодированных синхронных источников, моделируемых назад по времени;
(е) загрузку уточненной модели физических свойств или сохранение ее на запоминающем устройстве компьютера.
2. Способ по п.1, в котором инвертирование синхронной кодированной выборки измеренных геофизических данных содержит:
(ί) компьютерное моделирование синхронной кодированной выборки, соответствующей синхронной кодированной выборке измеренных данных, используя принятую модель физических свойств, где
- 10 017177 моделирование использует сигнатуры источников, кодированные теми же кодирующими функциями, которые использовались для кодирования синхронной кодированной выборки измеренных данных, где всю синхронную кодированную выборку моделируют в одной операции моделирования; и (ίί) вычисление функции стоимости, измеряющей степень несоответствия между синхронной кодированной выборкой измеренных данных и смоделированной синхронной кодированной выборкой.
3. Способ по п.1, дополнительно содержащий использование различных наборов неэквивалентных кодирующих функций по меньшей мере для одной итерации итеративного уточнения модели.
4. Способ по п.1, в котором инвертирование измеренных геофизических данных содержит:
(ί) вычисление прямого моделирования данных кодированного синхронного источника за один прогон моделирования, используя текущую модель физических свойств М и используя в качестве источника сигнатуру выборки синхронного источника, кодированную теми же кодирующими функциями сд, которые использовались для кодирования измеренных данных для нахождения ψ^;
(ίί) вычисление функции стоимости δ путем вычитания синхронной кодированной выборки измеренных геофизических данных из результата этапа (1);
(ΐϊϊ) вычисление обратного моделирования (то есть назад во времени), используя δ в качестве источника моделирования, производящий ψ^ο^ί (ίν) вычисление интеграла по времени от произведения ψ^ на ψα,ι|οιιι1 для нахождения градиента функции стоимости и (ν) использование градиента функции стоимости для настройки и уточнения модели М.
5. Способ по п.1, в котором указанные кодированные выборки измеренных данных кодируют посредством временного сворачивания всех трасс из выборки с соответствующей кодированной сигнатурой источника, причем указанная кодированная сигнатура источника является свертыванием функции источника для выборки с кодирующей функцией, отобранной для выборки.
6. Способ по п.1, в котором две или более кодированные выборки измеренных данных получают посредством получения выборок данных из геофизического исследования, в котором данные получают от множества одновременно работающих, уникально кодированных устройств источников.
7. Способ по п.1, в котором измеренные геофизические данные являются данными сейсморазведки.
8. Способ по п.7, в котором обобщенные сейсмические источники являются либо все точечными источниками, либо все источниками плоских волн.
9. Способ по п.2, в котором измеренные геофизические данные включают в себя измеренные или предполагаемые сигнатуры активации каждого источника, и кодированные сигнатуры источника, используемые в операциях моделирования, являются сигнатурами, полученными посредством временного свертывания измеренных или предполагаемых сигнатур источника с теми же кодирующими функциями, которые использовались для кодирования соответствующей измеренной выборки на этапе (а).
10. Способ по п.1, в котором кодирующие функции относятся к типу, выбранному из группы, состоящей из линейного кодирования, кодирования случайной фазы, кодирования линейной частотной модуляцией, модифицированного кодирования линейной частотной модуляцией, кодирования случайного временного сдвига и частотно-независимого фазового кодирования.
11. Способ по п.1, в котором кодирующие функции относятся к одному типу для одних источников и к другому типу для других источников, причем каждый тип выбирают из группы, состоящей из линейного кодирования, кодирования случайной фазы, кодирования линейной частотной модуляцией, модифицированного кодирования линейной частотной модуляцией, кодирования случайного временного сдвига и частотно-независимого фазового кодирования.
12. Способ по п.1, в котором кодирующие функции оптимизированы по отношению к используемой функции стоимости.
13. Способ по п.1, в котором операции прямого или обратного моделирования на этапе (б) выполняют по коду моделирования конечных разностей, конечных элементов или конечных объемов.
14. Способ по п.7, в котором модели физических свойств являются моделями скорости сейсмических волн, параметров сейсмической упругости, параметров сейсмической анизотропии или параметров сейсмической неупругости.
15. Способ по п.1, в котором для уточнения модели используют метод локальной оптимизации функции стоимости, такой как линейный градиентный поиск, метод сопряженных градиентов или метод Ньютона.
16. Способ по п.1, в котором функция стоимости является функцией стоимости с нормой Ь1 либо функцией стоимости с нормой Ь2, и функция стоимости может содержать члены регуляризации.
17. Способ по п.1, в котором неэквивалентные кодирующие функции являются практически ортогональными функциями.
18. Способ инверсии измеренных геофизических данных из области недр для разведки углеводородов, содержащий:
(а) получение группы из двух или более кодированных выборок измеренных геофизических данных, при этом каждая выборка ассоциирована с одним обобщенным источником или с одним приемни- 11 017177 ком, и при этом каждая выборка кодирована разными кодирующими функциями, выбираемыми из набора неэквивалентных кодирующих функций;
(b) суммирование кодированных выборок в группе посредством суммирования всех записей данных в каждой выборке, соответствующих одному приемнику (или источнику, если используется взаимность), и повторение для каждого приемника, приводящее к синхронной кодированной выборке;
(c) принятие модели физических свойств области недр, предоставление указанной моделью значений по меньшей мере одного физического свойства в местах по всей области недр и (й) инвертирование измеренных геофизических данных, по одной синхронной кодированной выборке за раз, используя принятую модель физических свойств в качестве начальной модели, и итеративное уточнение указанной модели для минимизирования функции стоимости, измеряющей степень несоответствия между моделируемыми данными и измеренными геофизическими данными, для формирования уточненной модели физических свойств, причем корректирование модели осуществляют с использованием градиента функции стоимости в отношении по меньшей мере одного параметра модели, при этом градиент вычисляют из интегрирования по времени произведения данных кодированных синхронных источников, моделированных вперед по времени, и данных кодированных синхронных источников, моделированных назад по времени;
(е) использование уточненной модели физических свойств для разведки углеводородов в области недр.
19. Способ получения углеводородов из области недр, содержащий:
(a) выполнение геофизического исследования области недр;
(b) получение модели физических свойств, причем указанная модель построена с использованием способа, указанного в п.1, который здесь включен в качестве ссылки;
(c) использование модели физических свойств для определения углеводородсодержащей зоны в области недр;
(й) бурение скважины в зоне и добыча углеводородов из скважины.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US7029708P | 2008-03-21 | 2008-03-21 | |
PCT/US2009/032010 WO2009117174A1 (en) | 2008-03-21 | 2009-01-26 | An efficient method for inversion of geophysical data |
Publications (2)
Publication Number | Publication Date |
---|---|
EA201070981A1 EA201070981A1 (ru) | 2011-02-28 |
EA017177B1 true EA017177B1 (ru) | 2012-10-30 |
Family
ID=41091213
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EA201070981A EA017177B1 (ru) | 2008-03-21 | 2009-01-26 | Эффективный способ инверсии геофизических данных |
Country Status (12)
Country | Link |
---|---|
US (1) | US8812282B2 (ru) |
EP (1) | EP2260331B1 (ru) |
KR (1) | KR101535407B1 (ru) |
CN (1) | CN101952743B (ru) |
AU (1) | AU2009226041B2 (ru) |
BR (1) | BRPI0908156A2 (ru) |
CA (1) | CA2711126C (ru) |
EA (1) | EA017177B1 (ru) |
ES (1) | ES2651923T3 (ru) |
MY (1) | MY157278A (ru) |
NO (1) | NO2260331T3 (ru) |
WO (1) | WO2009117174A1 (ru) |
Families Citing this family (54)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NZ575497A (en) | 2006-09-28 | 2011-08-26 | Exxonmobil Upstream Res Co | Iterative inversion of data from simultaneous geophysical sources |
CN102112894B (zh) | 2008-08-11 | 2015-03-25 | 埃克森美孚上游研究公司 | 用地震表面波的波形评估土壤性质 |
AU2010308501B2 (en) * | 2009-10-23 | 2016-03-17 | Exxonmobil Upstream Research Company | Method for optimization with gradient information |
US8830789B2 (en) * | 2010-02-10 | 2014-09-09 | Snu R&Db Foundation | Apparatus and method for imaging a subsurface using accumulated energy of wavefield |
US8223587B2 (en) * | 2010-03-29 | 2012-07-17 | Exxonmobil Upstream Research Company | Full wavefield inversion using time varying filters |
US8694299B2 (en) | 2010-05-07 | 2014-04-08 | Exxonmobil Upstream Research Company | Artifact reduction in iterative inversion of geophysical data |
EP2606452A4 (en) * | 2010-08-16 | 2017-08-16 | Exxonmobil Upstream Research Company | Reducing the dimensionality of the joint inversion problem |
RU2570827C2 (ru) * | 2010-09-27 | 2015-12-10 | Эксонмобил Апстрим Рисерч Компани | Гибридный способ для полноволновой инверсии с использованием способа одновременных и последовательных источников |
US8437998B2 (en) * | 2010-09-27 | 2013-05-07 | Exxonmobil Upstream Research Company | Hybrid method for full waveform inversion using simultaneous and sequential source method |
CN103119552B (zh) * | 2010-09-27 | 2016-06-08 | 埃克森美孚上游研究公司 | 同时源编码和源分离作为全波场反演的实际解决方案 |
BR112013008503A2 (pt) | 2010-12-01 | 2016-08-16 | Exxonmobil Upstream Res Co | inversão de fonte simultânea para dados de streamer marinho com função objetiva de correlação cruzada |
US10088586B2 (en) | 2011-02-25 | 2018-10-02 | University Of Florida Research Foundation, Inc. | Detection of sinkholes or anomalies |
WO2012134621A1 (en) * | 2011-03-30 | 2012-10-04 | Exxonmobil Upstream Research Company | Convergence rate of full wavefield inversion using spectral shaping |
CA2825469C (en) | 2011-03-31 | 2017-06-13 | Exxonmobil Upstream Research Company | Method of wavelet estimation and multiple prediction in full wavefield inversion |
US9158018B2 (en) | 2011-04-05 | 2015-10-13 | Westerngeco L.L.C. | Waveform inversion using a response of forward modeling |
WO2013016470A2 (en) * | 2011-07-27 | 2013-01-31 | Schlumberger Technology B.V. | Multi-well anisotropy inversion |
ES2640824T3 (es) | 2011-09-02 | 2017-11-06 | Exxonmobil Upstream Research Company | Utilización de la proyección sobre conjuntos convexos para limitar la inversión del campo de onda completa |
US9176930B2 (en) | 2011-11-29 | 2015-11-03 | Exxonmobil Upstream Research Company | Methods for approximating hessian times vector operation in full wavefield inversion |
US8862405B2 (en) * | 2011-12-06 | 2014-10-14 | Schlumberger Technology Corporation | System and method for producing look-ahead profile measurements in a drilling operation |
US9453928B2 (en) | 2012-03-06 | 2016-09-27 | Westerngeco L.L.C. | Methods and computing systems for processing data |
US10012745B2 (en) | 2012-03-08 | 2018-07-03 | Exxonmobil Upstream Research Company | Orthogonal source and receiver encoding |
US9091784B2 (en) * | 2012-06-15 | 2015-07-28 | Westerngeco L.L.C. | Determining an output representing a target structure based on encoded source and receiver data |
WO2014084945A1 (en) | 2012-11-28 | 2014-06-05 | Exxonmobil Upstream Resarch Company | Reflection seismic data q tomography |
MY169125A (en) | 2013-05-24 | 2019-02-18 | Exxonmobil Upstream Res Co | Multi-parameter inversion through offset dependent elastic fwi |
US10459117B2 (en) | 2013-06-03 | 2019-10-29 | Exxonmobil Upstream Research Company | Extended subspace method for cross-talk mitigation in multi-parameter inversion |
US9702998B2 (en) | 2013-07-08 | 2017-07-11 | Exxonmobil Upstream Research Company | Full-wavefield inversion of primaries and multiples in marine environment |
DK3036566T3 (en) | 2013-08-23 | 2018-07-23 | Exxonmobil Upstream Res Co | SIMILAR SOURCE APPLICATION DURING BOTH SEISMIC COLLECTION AND SEISMIC INVERSION |
US10036818B2 (en) | 2013-09-06 | 2018-07-31 | Exxonmobil Upstream Research Company | Accelerating full wavefield inversion with nonstationary point-spread functions |
US9910189B2 (en) | 2014-04-09 | 2018-03-06 | Exxonmobil Upstream Research Company | Method for fast line search in frequency domain FWI |
US10670752B2 (en) | 2014-04-14 | 2020-06-02 | Cgg Services Sas | Method for iterative inversion of data from non-encoded composite sources |
CN103926622B (zh) * | 2014-05-06 | 2015-03-11 | 赵婧文 | 一种基于l1范数多道匹配滤波压制多次波的方法 |
AU2015256626B2 (en) * | 2014-05-09 | 2017-10-12 | Exxonmobil Upstream Research Company | Efficient line search methods for multi-parameter full wavefield inversion |
US10185046B2 (en) | 2014-06-09 | 2019-01-22 | Exxonmobil Upstream Research Company | Method for temporal dispersion correction for seismic simulation, RTM and FWI |
WO2015199800A1 (en) | 2014-06-17 | 2015-12-30 | Exxonmobil Upstream Research Company | Fast viscoacoustic and viscoelastic full-wavefield inversion |
US10838092B2 (en) | 2014-07-24 | 2020-11-17 | Exxonmobil Upstream Research Company | Estimating multiple subsurface parameters by cascaded inversion of wavefield components |
US10422899B2 (en) | 2014-07-30 | 2019-09-24 | Exxonmobil Upstream Research Company | Harmonic encoding for FWI |
US10386511B2 (en) | 2014-10-03 | 2019-08-20 | Exxonmobil Upstream Research Company | Seismic survey design using full wavefield inversion |
US10108762B2 (en) * | 2014-10-03 | 2018-10-23 | International Business Machines Corporation | Tunable miniaturized physical subsurface model for simulation and inversion |
EP3210050A1 (en) | 2014-10-20 | 2017-08-30 | Exxonmobil Upstream Research Company | Velocity tomography using property scans |
WO2016065247A1 (en) * | 2014-10-24 | 2016-04-28 | Westerngeco Llc | Travel-time objective function for full waveform inversion |
CN105891883A (zh) * | 2014-12-08 | 2016-08-24 | 北京恒欣泰克科技有限公司 | 一种基于成本函数全局控制的三维地层模型构建方法 |
AU2015363241A1 (en) | 2014-12-18 | 2017-06-29 | Exxonmobil Upstream Research Company | Scalable scheduling of parallel iterative seismic jobs |
US10520618B2 (en) | 2015-02-04 | 2019-12-31 | ExxohnMobil Upstream Research Company | Poynting vector minimal reflection boundary conditions |
AU2015382333B2 (en) | 2015-02-13 | 2018-01-04 | Exxonmobil Upstream Research Company | Efficient and stable absorbing boundary condition in finite-difference calculations |
KR20190075176A (ko) | 2015-02-17 | 2019-06-28 | 엑손모빌 업스트림 리서치 캄파니 | 다중반사파 없는 데이터 세트를 생성하는 다단식 전 파동장 역산 프로세스 |
US9784867B2 (en) * | 2015-04-01 | 2017-10-10 | Schlumberger Technology Corporation | Seismic data processing |
RU2017145541A (ru) | 2015-06-04 | 2019-07-09 | Эксонмобил Апстрим Рисерч Компани | Способ построения свободных от многократных волн сейсмических изображений |
US10838093B2 (en) | 2015-07-02 | 2020-11-17 | Exxonmobil Upstream Research Company | Krylov-space-based quasi-newton preconditioner for full-wavefield inversion |
CA2998522A1 (en) | 2015-10-02 | 2017-04-06 | Exxonmobil Upstream Research Company | Q-compensated full wavefield inversion |
US20170103144A1 (en) * | 2015-10-08 | 2017-04-13 | Schlumbeger Technology Corporation | Well trajectory adjustment |
EP3362823B1 (en) | 2015-10-15 | 2019-10-09 | ExxonMobil Upstream Research Company | Fwi model domain angle stacks with amplitude preservation |
US10768324B2 (en) | 2016-05-19 | 2020-09-08 | Exxonmobil Upstream Research Company | Method to predict pore pressure and seal integrity using full wavefield inversion |
CN106199695B (zh) * | 2016-06-29 | 2018-03-09 | 中国石油化工股份有限公司 | 基于空变目标函数的叠前三参数反演方法 |
US11487036B2 (en) * | 2017-01-12 | 2022-11-01 | Cgg Services Sas | Reflection full waveform inversion methods with density and velocity models updated separately |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6021094A (en) * | 1998-12-03 | 2000-02-01 | Sandia Corporation | Method of migrating seismic records |
US6549854B1 (en) * | 1999-02-12 | 2003-04-15 | Schlumberger Technology Corporation | Uncertainty constrained subsurface modeling |
US6944546B2 (en) * | 2003-10-01 | 2005-09-13 | Halliburton Energy Services, Inc. | Method and apparatus for inversion processing of well logging data in a selected pattern space |
US20070203673A1 (en) * | 2005-11-04 | 2007-08-30 | Sherrill Francis G | 3d pre-stack full waveform inversion |
Family Cites Families (72)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3812457A (en) | 1969-11-17 | 1974-05-21 | Shell Oil Co | Seismic exploration method |
US4168485A (en) | 1974-08-12 | 1979-09-18 | Continental Oil Company | Simultaneous use of pseudo-random control signals in vibrational exploration methods |
US4545039A (en) | 1982-09-09 | 1985-10-01 | Western Geophysical Co. Of America | Methods for seismic exploration |
US4562540A (en) | 1982-11-12 | 1985-12-31 | Schlumberger Technology Corporation | Diffraction tomography system and methods |
US4594662A (en) | 1982-11-12 | 1986-06-10 | Schlumberger Technology Corporation | Diffraction tomography systems and methods with fixed detector arrays |
US4715020A (en) | 1986-10-29 | 1987-12-22 | Western Atlas International, Inc. | Simultaneous performance of multiple seismic vibratory surveys |
FR2589587B1 (fr) | 1985-10-30 | 1988-02-05 | Inst Francais Du Petrole | Procede de prospection sismique marine utilisant un signal vibratoire code et dispositif pour sa mise en oeuvre |
US4707812A (en) | 1985-12-09 | 1987-11-17 | Atlantic Richfield Company | Method of suppressing vibration seismic signal correlation noise |
US4823326A (en) | 1986-07-21 | 1989-04-18 | The Standard Oil Company | Seismic data acquisition technique having superposed signals |
US4686654A (en) | 1986-07-31 | 1987-08-11 | Western Geophysical Company Of America | Method for generating orthogonal sweep signals |
US4953657A (en) | 1987-11-30 | 1990-09-04 | Halliburton Geophysical Services, Inc. | Time delay source coding |
GB2322704B (en) | 1994-07-07 | 1998-12-09 | Geco As | Method of Processing seismic data |
US5924049A (en) * | 1995-04-18 | 1999-07-13 | Western Atlas International, Inc. | Methods for acquiring and processing seismic data |
WO1996033425A1 (en) | 1995-04-18 | 1996-10-24 | Western Atlas International, Inc. | Uniform subsurface coverage at steep dips |
US5721710A (en) | 1995-09-29 | 1998-02-24 | Atlantic Richfield Company | High fidelity vibratory source seismic method with source separation |
US5838634A (en) | 1996-04-04 | 1998-11-17 | Exxon Production Research Company | Method of generating 3-D geologic models incorporating geologic and geophysical constraints |
US5798982A (en) | 1996-04-29 | 1998-08-25 | The Trustees Of Columbia University In The City Of New York | Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models |
GB9612471D0 (en) | 1996-06-14 | 1996-08-14 | Geco As | Method and apparatus for multiple seismic vibratory surveys |
US5878372A (en) | 1997-03-04 | 1999-03-02 | Western Atlas International, Inc. | Method for simultaneous inversion processing of well log data using a plurality of earth models |
US6014342A (en) | 1997-03-21 | 2000-01-11 | Tomo Seis, Inc. | Method of evaluating a subsurface region using gather sensitive data discrimination |
FR2765692B1 (fr) | 1997-07-04 | 1999-09-10 | Inst Francais Du Petrole | Methode pour modeliser en 3d l'impedance d'un milieu heterogene |
GB2329043B (en) | 1997-09-05 | 2000-04-26 | Geco As | Method of determining the response caused by model alterations in seismic simulations |
US5999488A (en) | 1998-04-27 | 1999-12-07 | Phillips Petroleum Company | Method and apparatus for migration by finite differences |
US6388947B1 (en) | 1998-09-14 | 2002-05-14 | Tomoseis, Inc. | Multi-crosswell profile 3D imaging and method |
US6574564B2 (en) | 1998-10-01 | 2003-06-03 | Institut Francais Du Petrole | 3D prestack seismic data migration method |
FR2784195B1 (fr) | 1998-10-01 | 2000-11-17 | Inst Francais Du Petrole | Methode pour realiser en 3d avant sommation, une migration de donnees sismiques |
US6058073A (en) | 1999-03-30 | 2000-05-02 | Atlantic Richfield Company | Elastic impedance estimation for inversion of far offset seismic sections |
FR2792419B1 (fr) | 1999-04-16 | 2001-09-07 | Inst Francais Du Petrole | Methode pour obtenir un modele optimal d'une caracteristique physique dans un milieu heterogene, tel que le sous-sol |
GB9927395D0 (en) | 1999-05-19 | 2000-01-19 | Schlumberger Holdings | Improved seismic data acquisition method |
US6327537B1 (en) | 1999-07-19 | 2001-12-04 | Luc T. Ikelle | Multi-shooting approach to seismic modeling and acquisition |
ATE337566T1 (de) | 1999-10-22 | 2006-09-15 | Jason Geosystems B V | Verfahren zur bestimmung der elastischen parameter und felszusammensetzung von unterirdischen formationen mit hilfe von seismischen daten |
US6480790B1 (en) | 1999-10-29 | 2002-11-12 | Exxonmobil Upstream Research Company | Process for constructing three-dimensional geologic models having adjustable geologic interfaces |
FR2800473B1 (fr) | 1999-10-29 | 2001-11-30 | Inst Francais Du Petrole | Methode pour modeliser en 2d ou 3d un milieu heterogene tel que le sous-sol decrit par plusieurs parametres physiques |
US6836448B2 (en) | 2000-01-21 | 2004-12-28 | Schlumberger Technology Corporation | System and method for seismic wavefield separation |
CN1188710C (zh) | 2000-01-21 | 2005-02-09 | 施鲁博格控股有限公司 | 估算地震介质特性的系统和方法 |
US6826486B1 (en) | 2000-02-11 | 2004-11-30 | Schlumberger Technology Corporation | Methods and apparatus for predicting pore and fracture pressures of a subsurface formation |
FR2805051B1 (fr) | 2000-02-14 | 2002-12-06 | Geophysique Cie Gle | Methode de surveillance sismique d'une zone souterraine par utilisation simultanee de plusieurs sources vibrosismiques |
GB2359363B (en) | 2000-02-15 | 2002-04-03 | Geco Prakla | Processing simultaneous vibratory seismic data |
US6687659B1 (en) | 2000-03-24 | 2004-02-03 | Conocophillips Company | Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications |
US6317695B1 (en) | 2000-03-30 | 2001-11-13 | Nutec Sciences, Inc. | Seismic data processing method |
EP1327162A2 (en) | 2000-10-17 | 2003-07-16 | WesternGeco, L.L.C. | Method of using cascaded sweeps for source coding and harmonic cancellation |
AU2002239619A1 (en) | 2000-12-08 | 2002-06-18 | Peter J. Ortoleva | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
FR2818753B1 (fr) | 2000-12-21 | 2003-03-21 | Inst Francais Du Petrole | Methode et dispositif de prospection sismique par emission simultanee de signaux sismisques obtenus en codant un signal par des sequences pseudo aleatoires |
US6545944B2 (en) | 2001-05-30 | 2003-04-08 | Westerngeco L.L.C. | Method for acquiring and processing of data from two or more simultaneously fired sources |
GB2379013B (en) | 2001-08-07 | 2005-04-20 | Abb Offshore Systems Ltd | Microseismic signal processing |
US6842701B2 (en) | 2002-02-25 | 2005-01-11 | Westerngeco L.L.C. | Method of noise removal for cascaded sweep data |
GB2387226C (en) | 2002-04-06 | 2008-05-12 | Westerngeco Ltd | A method of seismic surveying |
US6832159B2 (en) * | 2002-07-11 | 2004-12-14 | Schlumberger Technology Corporation | Intelligent diagnosis of environmental influence on well logs with model-based inversion |
FR2843202B1 (fr) | 2002-08-05 | 2004-09-10 | Inst Francais Du Petrole | Methode pour former un modele representatif de la distribution d'une grandeur physique dans une zone souterraine, affranchi de l'effet de bruits correles entachant des donnees d'exploration |
US6859734B2 (en) | 2002-10-04 | 2005-02-22 | Paradigm Geophysical Corporation | Method and system for limited frequency seismic imaging |
US6999880B2 (en) | 2003-03-18 | 2006-02-14 | The Regents Of The University Of California | Source-independent full waveform inversion of seismic data |
US7072767B2 (en) * | 2003-04-01 | 2006-07-04 | Conocophillips Company | Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data |
US6882938B2 (en) | 2003-07-30 | 2005-04-19 | Pgs Americas, Inc. | Method for separating seismic signals from two or more distinct sources |
US6901333B2 (en) | 2003-10-27 | 2005-05-31 | Fugro N.V. | Method and device for the generation and application of anisotropic elastic parameters |
US7046581B2 (en) | 2003-12-01 | 2006-05-16 | Shell Oil Company | Well-to-well tomography |
US20050128874A1 (en) | 2003-12-15 | 2005-06-16 | Chevron U.S.A. Inc. | Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources |
US7646924B2 (en) | 2004-08-09 | 2010-01-12 | David Leigh Donoho | Method and apparatus for compressed sensing |
GB2422433B (en) | 2004-12-21 | 2008-03-19 | Sondex Wireline Ltd | Method and apparatus for determining the permeability of earth formations |
US7373251B2 (en) | 2004-12-22 | 2008-05-13 | Marathon Oil Company | Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data |
EP1859301B1 (en) | 2005-02-22 | 2013-07-17 | Paradigm Geophysical Ltd. | Multiple suppression in angle domain time and depth migration |
WO2006122146A2 (en) | 2005-05-10 | 2006-11-16 | William Marsh Rice University | Method and apparatus for distributed compressed sensing |
GB2436626B (en) | 2006-03-28 | 2008-08-06 | Westerngeco Seismic Holdings | Method of evaluating the interaction between a wavefield and a solid body |
US20070274155A1 (en) | 2006-05-25 | 2007-11-29 | Ikelle Luc T | Coding and Decoding: Seismic Data Modeling, Acquisition and Processing |
US7725266B2 (en) | 2006-05-31 | 2010-05-25 | Bp Corporation North America Inc. | System and method for 3D frequency domain waveform inversion based on 3D time-domain forward modeling |
GB2455664B (en) * | 2006-09-13 | 2011-02-16 | Exxonmobil Upstream Res Co | Rapid inversion of electromagnetic reconnaissance survey data |
NZ575497A (en) * | 2006-09-28 | 2011-08-26 | Exxonmobil Upstream Res Co | Iterative inversion of data from simultaneous geophysical sources |
EP2104869B1 (en) | 2007-01-20 | 2012-01-25 | Spectraseis AG | Time reverse reservoir localization |
US7636275B2 (en) * | 2007-02-06 | 2009-12-22 | Conocophillips Company | Direct time lapse inversion of seismic data |
US7756642B2 (en) * | 2007-06-27 | 2010-07-13 | Schlumberger Technology Corporation | Characterizing an earth subterranean structure by iteratively performing inversion based on a function |
JP2009063942A (ja) | 2007-09-10 | 2009-03-26 | Sumitomo Electric Ind Ltd | 遠赤外線カメラ用レンズ、レンズユニット及び撮像装置 |
US9638830B2 (en) * | 2007-12-14 | 2017-05-02 | Westerngeco L.L.C. | Optimizing drilling operations using petrotechnical data |
EP2105765A1 (en) | 2008-03-28 | 2009-09-30 | Schlumberger Holdings Limited | Simultaneous inversion of induction data for dielectric permittivity and electric conductivity |
-
2009
- 2009-01-26 ES ES09723513.9T patent/ES2651923T3/es active Active
- 2009-01-26 KR KR1020107020908A patent/KR101535407B1/ko active IP Right Grant
- 2009-01-26 AU AU2009226041A patent/AU2009226041B2/en not_active Ceased
- 2009-01-26 NO NO09723513A patent/NO2260331T3/no unknown
- 2009-01-26 BR BRPI0908156-9A patent/BRPI0908156A2/pt not_active Application Discontinuation
- 2009-01-26 EA EA201070981A patent/EA017177B1/ru not_active IP Right Cessation
- 2009-01-26 US US12/810,721 patent/US8812282B2/en active Active
- 2009-01-26 EP EP09723513.9A patent/EP2260331B1/en active Active
- 2009-01-26 CA CA2711126A patent/CA2711126C/en not_active Expired - Fee Related
- 2009-01-26 CN CN200980106257.2A patent/CN101952743B/zh not_active Expired - Fee Related
- 2009-01-26 WO PCT/US2009/032010 patent/WO2009117174A1/en active Application Filing
- 2009-02-23 MY MYPI20090712A patent/MY157278A/en unknown
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6021094A (en) * | 1998-12-03 | 2000-02-01 | Sandia Corporation | Method of migrating seismic records |
US6549854B1 (en) * | 1999-02-12 | 2003-04-15 | Schlumberger Technology Corporation | Uncertainty constrained subsurface modeling |
US6944546B2 (en) * | 2003-10-01 | 2005-09-13 | Halliburton Energy Services, Inc. | Method and apparatus for inversion processing of well logging data in a selected pattern space |
US20070203673A1 (en) * | 2005-11-04 | 2007-08-30 | Sherrill Francis G | 3d pre-stack full waveform inversion |
Also Published As
Publication number | Publication date |
---|---|
MY157278A (en) | 2016-05-31 |
EA201070981A1 (ru) | 2011-02-28 |
CA2711126C (en) | 2016-07-19 |
US20110000678A1 (en) | 2011-01-06 |
BRPI0908156A2 (pt) | 2015-08-11 |
EP2260331A4 (en) | 2016-09-28 |
KR101535407B1 (ko) | 2015-07-09 |
CA2711126A1 (en) | 2009-09-24 |
EP2260331B1 (en) | 2017-10-11 |
EP2260331A1 (en) | 2010-12-15 |
CN101952743B (zh) | 2013-05-01 |
AU2009226041A1 (en) | 2009-09-24 |
CN101952743A (zh) | 2011-01-19 |
NO2260331T3 (ru) | 2018-03-10 |
WO2009117174A1 (en) | 2009-09-24 |
ES2651923T3 (es) | 2018-01-30 |
KR20110000634A (ko) | 2011-01-04 |
US8812282B2 (en) | 2014-08-19 |
AU2009226041B2 (en) | 2013-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EA017177B1 (ru) | Эффективный способ инверсии геофизических данных | |
RU2587498C2 (ru) | Инверсия одновременных источников для данных сейсмоприемной косы с взаимнокорреляционной целевой функцией | |
KR101409010B1 (ko) | 표면 아래 영역에 대한 물리적 특성 모델을 결정하기 위해 측정된 지구물리학 데이터의 반전을 컴퓨터로 수행하는 방법 및 표면 아래 영역으로부터 탄화 수소를 생성하는 방법 | |
RU2582480C2 (ru) | Кодирование одновременных источников и разделение источников в качестве практического решения по инверсии полного волнового поля | |
RU2570827C2 (ru) | Гибридный способ для полноволновой инверсии с использованием способа одновременных и последовательных источников | |
Przebindowska | Acoustic full waveform inversion of marine reflection seismic data | |
Decker et al. | A variational approach for picking optimal surfaces from semblance-like panels |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s) |
Designated state(s): AM AZ BY KZ KG MD TJ TM |
|
MM4A | Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s) |
Designated state(s): RU |