Кодирование одновременных источников и разделение источников в качестве практического решения по инверсии полного волнового поля

Кодирование одновременных источников и разделение источников в качестве практического решения по инверсии полного волнового поля

Реферат

Перекрестные ссылки на родственные заявки

Данная заявка притязает на приоритет предварительной заявки на патент (США) 61/386831, поданной 27 сентября 2010 года, которая содержится в данном документе по ссылке во всех разрешенных юрисдикциях. Данная заявка связана с конкретной заявкой, озаглавленной "Hybrid Method for Full waveform Inversion Using Simultaneous and Sequential source Method", и притязает на приоритет предварительной заявки на патент (США) 61/836828, поданной 27 сентября 2010 года. Эта родственная заявка также содержится по ссылке в данном документе во всех разрешенных юрисдикциях.

Область техники, к которой относится изобретение

Изобретение, в общем, относится к области техники геофизических исследований, а более конкретно - к обработке геофизических данных. В частности, изобретение представляет собой способ для инверсии данных, полученных из нескольких геофизических источников, таких как сейсмические источники, заключающий в себе геофизическое моделирование, которое вычисляет данные из многих одновременно активных геофизических источников при одном выполнении моделирования.

Уровень техники

Геофизическая инверсия [1, 2] нацелена на нахождение модели свойств геологической среды, которая оптимально объясняет наблюдаемые данные и удовлетворяет геологическим и геофизическим ограничениям. Существует большое число известных способов геофизической инверсии. Эти известные способы разделяются на одну из двух категорий - итеративную инверсию и неитеративную инверсию. Ниже приведены определения того, что обычно подразумевается под каждой из двух категорий.

Неитеративная инверсия является инверсией, которая выполняется посредством допущения некоторой простой фоновой модели и обновления модели на основе входных данных. Этот способ не использует обновленную модель в качестве ввода в другой этап инверсии. Для случая сейсмических данных эти способы обычно упоминаются как построение изображений, миграция, дифракционная томография или инверсия Борна.

Итеративная инверсия является инверсией, заключающей в себе повторяющееся улучшение модели свойств геологической среды таким образом, что находится модель, которая удовлетворительно объясняет наблюдаемые данные. Если инверсия сходится, то конечная модель должна лучше объяснять наблюдаемые данные и должна более близко аппроксимировать фактические свойства геологической среды. Итеративная инверсия обычно формирует более точную модель, чем неитеративная инверсия, но является гораздо более затратной в вычислительном отношении.

Итеративная инверсия, в общем, является предпочтительной по сравнению с неитеративной инверсией, поскольку она дает в результате более точные параметрические модели геологической среды. К сожалению, итеративная инверсия является настолько вычислительно затратной, что непрактично ее применение ко многим интересующим задачам. Эти значительные вычислительные ресурсы являются результатом того факта, что все технологии инверсии требуют множества моделирований с большим объемом вычислений. Время вычисления любого отдельного моделирования является пропорциональным числу источников, которые должны быть инвертированы, и типично присутствует большое число источников в геофизических данных, при этом термин "источник" при использовании выше означает местоположение активации аппаратной системы-источника. Проблема обостряется для итеративной инверсии, поскольку число моделирований, которые должны быть вычислены, является пропорциональным числу итераций в инверсии, и число требуемых итераций типично составляет порядка нескольких сотен или тысяч.

Типично используемым способом итеративной инверсии, применяемым в геофизике, является оптимизация функции затрат. Оптимизация функции затрат заключает в себе итеративную минимизацию или максимизацию значения, относительно модели M, функции затрат S(M), которая является мерой несоответствия между расчетными и наблюдаемыми данными (она также иногда упоминается как целевая функция), при этом расчетные данные моделируются на компьютере с использованием текущей модели геофизических свойств и физики, влияющей на распространение сигнала источника в среде, представленной посредством данной модели геофизических свойств. Моделирующие (имитационные) вычисления могут выполняться посредством любого из нескольких численных методов, включающих в себя, но не только, конечную разность, конечный элемент или трассировку лучей. Моделирующие вычисления могут быть выполнены в частотной или временной области.

Способы оптимизации функции затрат являются или локальными, или глобальными [3]. Глобальные способы просто заключают в себе вычисление функции затрат S(M) для совокупности моделей {M₁, M₂, M₃,...} и выбор набора из одной или более моделей из этой совокупности, которые приблизительно минимизируют S(M). Если требуется дополнительное улучшение, этот новый выбранный набор моделей затем может быть использован в качестве основы, чтобы формировать новую совокупность моделей, которые могут быть снова протестированы относительно функции затрат S(M). Для глобальных способов каждая модель в тестовой совокупности может рассматриваться в качестве итерации, или на верхнем уровне каждый тестируемый набор совокупностей может считаться итерацией. Известные глобальные способы инверсии включают в себя метод Монте-Карло, метод имитации отжига, генетические и эволюционные алгоритмы.

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

1. Выбор начальной модели.

2. Вычисление градиента функции затрат S(M) относительно параметров, которые описывают модель.

3. Поиск обновленной модели, которая является возмущением начальной модели в направлении отрицательного градиента, которая лучше объясняет наблюдаемые данные.

Алгоритм 1. Алгоритм для выполнения локальной оптимизации функции затрат

Эта процедура итеративно выполняется посредством использования новой обновленной модели в качестве начальной модели для другого градиентного поиска. Процесс продолжается до тех пор, пока не находится обновленная модель, которая удовлетворительно объясняет наблюдаемые данные. Обычно используемые локальные способы инверсии функции затрат включают в себя поиск градиента, сопряженные градиенты и метод Ньютона.

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

S ( M ) = ∑ g = 1 N g ∑ r = 1 N r ∑ t = 1 N t W ( ψ c a l c ( M , r , t , w g ) − ψ o b s ( r , t , w g ) ) , уравнение 1

где S является функцией затрат,

M является вектором N параметров (m₁, m₂,..., m_N), описывающим модель геологической среды,

g является индексом сейсмограммы,

w_g является функцией источника для сейсмограммы g, которая является функцией от пространственных координат и времени, для точечного источника она является дельта-функцией от пространственных координат,

N_g является числом сейсмограмм,

r является индексом приемника в сейсмограмме,

N_r является числом приемников в сейсмограмме,

t является индексом временной выборки в трассе,

N_t является числом временных выборок,

W является функцией критериев минимизации (обычно выбирается W(x)=x², которая представляет собой критерии по методу наименьших квадратов (L2)),

ψ_calc является расчетными данными сейсмического давления из модели M,

ψ_obs является измеренными данными сейсмического давления.

Сейсмограммы могут быть любым типом сейсмограммы, которая может быть моделирована в одном проходе программы сейсмического прямого моделирования. Обычно сейсмограммы соответствуют сейсмическому взрыву, хотя взрывы могут быть более общими, чем точечные источники. Для точечных источников индекс g сейсмограммы соответствует местоположению отдельных точечных источников. Для плоских сейсмоволн источники g должны соответствовать различным направлениям распространения плоских сейсмоволн. Эти данные обобщенных источников, ψ_obs, могут или быть обнаружены в поле, или могут быть синтезированы из данных, обнаруженных с использованием точечных источников. Расчетные данные ψ_calc, с другой стороны, обычно могут вычисляться непосредственно посредством использования функции обобщенных источников при прямом моделировании. Для многих типов прямого моделирования, включающих в себя моделирование на основе конечных разностей, время вычисления, необходимое для обобщенного источника, примерно равно времени вычисления, необходимому для точечного источника.

Уравнение (1) может быть упрощено до:

S ( M ) = ∑ g = 1 N g W ( δ ( M , w g ) ) , уравнение 2

причем теперь подразумевается сумма по приемникам и временным выборкам, и

δ ( M , w g ) = ψ c a l c ( M , w g ) − ψ o b s ( w g ) , уравнение 3

Инверсия пытается обновлять модель M таким образом, что S(M) является минимумом. Это может быть выполнено посредством локальной оптимизации функции затрат, которая обновляет данную модель M^(k) следующим образом:

M ( k + 1 ) = M ( k ) − α ( k ) ∇ M S ( M ) , уравнение 4

где k является номером итерации, α является скалярным размером обновления модели, и ∇_MS(M) является градиентом функции несоответствия, принятой относительно параметров модели. Возмущения модели или значения, посредством которых обновляется модель, вычисляются посредством умножения градиента целевой функции на длину α шага, которая должна вычисляться многократно.

Из уравнения (2) следующее уравнение может извлекаться для градиента функции затрат:

∇ M S ( M ) = ∑ g = 1 N g ∇ M W ( δ ( M , w g ) ) , уравнение 5

Таким образом, чтобы вычислять градиент функции затрат, следует отдельно вычислять градиент доли каждой сейсмограммы в функции затрат, затем суммировать эти доли. Следовательно, трудоемкость вычислений, требуемая для вычисления ∇_MS(M), в N_g раз превышает трудоемкость вычислений, требуемую для того, чтобы определять долю одной сейсмограммы в градиенте. Для геофизических задач N_g обычно соответствует числу геофизических источников и составляет порядка 10000-100000, значительно повышая затраты на вычисление ∇_MS(M).

Следует отметить, что вычисление ∇_MW(δ) требует вычисления производной W(δ) относительно каждого из N параметров модели m_i. Поскольку для геофизических задач N обычно является очень большим (обычно больше одного миллиона), это вычисление может быть чрезвычайно длительным, если оно должно быть выполнено для каждого отдельного параметра модели. К счастью, сопряженный способ может быть использован для того, чтобы эффективно выполнять это вычисление для всех параметров модели сразу [1]. Сопряженный способ для целевой функции на основе метода наименьших квадратов и параметризации модели с координатной привязкой обобщается посредством следующего алгоритма.

1. Вычисление прямого моделирования данных с использованием текущей модели и формы импульса w_g сейсмограммы в качестве источника, чтобы получать ψ_calc(M^(k),w_g).

2. Вычитание наблюдаемых данных из моделированных данных, дающее δ(M^(k),w_g).

3. Вычисление обратного моделирования (т.е. назад во времени) с использованием δ(M^(k),w_g) в качестве источника, формирующего ψ_adjoint(M^(k),w_g).

4. Вычисление интеграла во времени произведения ψ_calc(M^(k),w_g) и ψ_adjoint(M^(k),w_g), чтобы получать ∇_MW(δ(M,w_g)).

Алгоритм 2. Алгоритм для вычисления градиента функции затрат на основе метода наименьших квадратов модели с координатной привязкой с использованием сопряженного способа

Хотя вычисление градиентов с использованием сопряженного способа является эффективным относительно других способов, оно по-прежнему является очень затратным. В частности, сопряженные способы требуют двух моделирований, одного вперед во времени и одного назад во времени, и для геофизических задач эти моделирования обычно подвержены очень большому объему вычислений. Кроме того, как пояснено выше, это вычисление на основе сопряженного способа должно быть выполнено для каждой сейсмограммы с измеренными данными по отдельности, повышая вычислительные затраты на коэффициент N_g.

Вычислительные затраты всех категорий инверсии могут быть уменьшены посредством инвертирования данных из комбинаций источников вместо инвертирования источников по отдельности. Это может называться инверсией одновременных источников. Известно несколько типов комбинаций источников, включающих в себя: когерентное суммирование близко расположенных источников, чтобы формировать эффективный источник, который формирует фронт волны некоторой требуемой формы (например, плоскую сейсмоволну), суммирование широко разнесенных источников либо полное или частичное накопление данных перед инверсией.

Уменьшение вычислительных затрат, полученное посредством инвертирования комбинированных источников, по меньшей мере, частично смещается вследствие того факта, что инверсия комбинированных данных обычно формирует менее точную инвертированную модель. Эти потери в точности обусловлены тем фактом, что информация теряется, когда суммируются отдельные источники, и, следовательно, суммированные данные не ограничивают инвертированную модель так строго, как несуммированные данные. Эти потери информации во время суммирования могут быть минимизированы посредством кодирования каждой записи по взрыву перед суммированием. Кодирование перед комбинированием сохраняет значительно больше информации в данных одновременных источников и, следовательно, лучше ограничивает инверсию [5]. Кодирование также дает возможность комбинирования близко расположенных источников, тем самым давая возможность комбинирования большего числа источников для данной вычислительной области. Для этой технологии могут быть использованы различные схемы кодирования, включающие в себя кодирование на основе сдвига по времени и случайное фазовое кодирование. Оставшаяся часть этого раздела "Уровень техники" кратко анализирует различные опубликованные технологии на основе геофизических одновременных источников, как кодированных, так и некодированных.

Van Manen [6] предлагает использование способа сейсмической интерферометрии для того, чтобы ускорять прямое моделирование. Сейсмическая интерферометрия работает посредством размещения источников по всей границе интересующей области. Эти источники моделируются по отдельности, и записывается волновое поле во всех местоположениях, для которых требуется функция Грина. Функция Грина между любыми двумя записанными местоположениями затем может быть вычислена посредством взаимной корреляции трасс, обнаруженных в двух записанных местоположениях, и суммирования по всем граничным источникам. Если данные, которые должны быть инвертированы, имеют большое число источников и приемников, которые находятся в интересующей области (в противоположность наличию одних или других на границе), это представляет собой очень эффективный метод для вычисления требуемых функций Грина. Тем не менее, для случая сейсмических данных редко когда как источник, так и приемник для данных, которые должны быть инвертированы, находятся в интересующей области. Следовательно, это улучшение имеет очень ограниченную применимость для задачи сейсмической инверсии.

Berkhout [7] и Zhang [8] предлагают то, что инверсия, в общем, может быть улучшена посредством инвертирования некодированных одновременных источников, которые суммируются когерентно, чтобы формировать некоторый требуемый фронт волны в некоторой области геологической среды. Например, данные точечных источников могут быть суммированы со сдвигами по времени, которые являются линейной функцией от местоположения источников, чтобы формировать нисходящую плоскую сейсмоволну под некоторым конкретным углом относительно поверхности. Эта технология может применяться ко всем категориям инверсии. Проблема этого способа заключается в том, что когерентное суммирование сейсмограмм источников обязательно уменьшает объем информации в данных. Таким образом, например, суммирование для того, чтобы формировать плоскую сейсмоволну, удаляет всю информацию в сейсмических данных, связанную со временем пробега в зависимости от выноса источников от приемников. Эта информация является критически важной для обновления медленно варьирующейся фоновой модели скорости, и, следовательно, способ Berkhout не подчиняется оптимальным ограничениям. Чтобы преодолевать эту проблему, могут инвертироваться множество различных когерентных сумм данных (например, множество плоских сейсмоволн с различными направлениями распространения), но в таком случае теряется эффективность, поскольку затраты на инверсию являются пропорциональными числу различных инвертированных сумм. В данном документе такие когерентно суммированные источники называются обобщенными источниками. Следовательно, обобщенный источник может быть либо точечным источником, либо суммой точечных источников, которая формирует фронт волны некоторой требуемой формы.

Van Riel [9] предлагает инверсию посредством некодированного накопления или частичного накопления (относительно выноса источников от приемников) входных сейсмических данных, а затем задания функции затрат относительно этих накопленных данных, которые должны быть оптимизированы. Таким образом, эта публикация предлагает улучшение инверсии на основе функции затрат с использованием некодированных одновременных источников. Как и для способа инверсии одновременных источников по Berkhout [7], накопление, предлагаемое посредством этого способа, уменьшает объем информации в данных, которые должны быть инвертированы, и, следовательно, инверсия подчиняется субоптимальным ограничениям по сравнению с исходными данными.

Mora [10] предлагает инвертирование данных, которые являются суммой широко разнесенных источников. Таким образом, эта публикация предлагает повышение эффективности инверсии с использованием моделирования на основе некодированных одновременных источников. Суммирование широко разнесенных источников обеспечивает преимущество сохранения гораздо большего объема информации, чем когерентная сумма, предложенная Berkhout. Тем не менее, суммирование широко разнесенных источников подразумевает, что апертура (инвертируемая область модели), которая должна быть использована в инверсии, должна быть увеличена, чтобы охватывать все широко разнесенные источники. Поскольку время вычисления является пропорциональным площади этой апертуры, способ Mora не формирует такой прирост эффективности, который может достигаться, если суммированные источники находятся рядом друг с другом.

Ober [11] предлагает ускорение сейсмической миграции, частного случая неитеративной инверсии, посредством использования одновременных кодированных источников. После тестирования различных способов кодирования Ober выяснил, что результирующие мигрированные изображения значительно уменьшают отношение "сигнал-шум" вследствие того факта, что функции широкополосного кодирования обязательно являются только приблизительно ортогональными. Таким образом, при суммировании более 16 взрывов качество инверсии не является удовлетворительным. Поскольку неитеративная инверсия не является слишком затратной для начального уровня, и поскольку требуется инверсия высокого отношения "сигнал-шум", эта технология не используется широко на практике в геофизической отрасли.

Ikelle [12] предлагает способ для быстрого прямого моделирования посредством одновременного моделирования точечных источников, которые активируются (при моделировании) с варьирующимися временными интервалами. Также поясняется способ для декодирования этих сдвинутых по времени моделированных на основе одновременных источников данных обратно в раздельные моделирования, которые получены из отдельных точечных источников. Эти декодированные данные затем могут быть использованы в качестве части любой традиционной процедуры инверсии. Проблема способа Ikelle заключается в том, что предложенный способ декодирования должен формировать разделенные данные, имеющие уровни шума, пропорциональные разности между данными из смежных источников. Этот шум становится значительным для моделей геологической среды, которые не являются горизонтально постоянными, например из моделей, содержащих наклонные отражающие границы. Кроме того, этот шум должен расти пропорционально числу одновременных источников. Вследствие этих трудностей подход на основе одновременных источников Ikelle может приводить к недопустимым уровням шума, если используется при инвертировании геологической среды, которая не является горизонтально постоянной.

Одновременная инверсия нескольких кодированных источников предлагается Krebs и др. в публикации PCT-заявки на патент номер WO 2008/042081 [5, 18], которая содержится в данном документе по ссылке во всех разрешенных юрисдикциях, представляет собой очень эффективный с точки зрения затрат способ для того, чтобы инвертировать данные полного волнового поля. (Идентичный подход одновременной инверсии кодированной сейсмограммы должен работать для приемников либо через обратимость источников-приемников, либо посредством кодирования фактических местоположений приемников в сейсмограммах общих источников данных.) Для стационарных приемников прямые и сопряженные вычисления должны выполняться только для одного эффективного источника; см. публикацию PCT-заявки на патент номер WO 2009/117174, которая содержится в данном документе по ссылке во всех разрешенных юрисдикциях. Учитывая тот факт, что сотни взрывов записываются для типичных геометрий двумерного обнаружения и тысячи в случае трехмерных съемок, экономия вычислительных ресурсов от этого способа является довольно значительной. На практике предположение о стационарных приемниках не является строго достоверным для наиболее распространенных геометрий обнаружения полевых данных. В случае данных морского сейсморазведочного кабеля источники и приемники перемещаются для каждого нового взрыва. Даже при съемках, в которых местоположения приемников являются стационарными, практика зачастую такова, что не все приемники "прослушивают" каждый взрыв, и приемники, которые прослушивают, могут варьироваться в зависимости от взрыва. Это также нарушает "предположение о стационарных приемниках". Помимо этого вследствие логистических проблем трудно записывать данные близко к источнику, и это означает то, что данные минимальных выносов типично отсутствуют. Это является истинным и для морских, и для наземных съемок. Оба из этих факторов означают, что для сейсмограммы одновременных источников каждое местоположение приемника представляет собой отсутствующие данные для некоторых исходных взрывов. В общих словах в инверсии одновременных кодированных источников, для данной одновременной кодированной сейсмограммы, данные требуются во всех местоположениях приемников для каждого взрыва, и это может упоминаться как предположение о стационарных приемниках инверсии одновременных кодированных источников. В WO 08/042081 некоторые раскрытые варианты осуществления могут работать лучше, чем другие, когда не удовлетворяется предположение о стационарных приемниках. Следовательно, должно быть преимущественным иметь приспособление или регулирование для прямого применения инверсии одновременных кодированных источников (и/или приемников), что должно повышать производительность, когда скомпрометировано предположение о стационарных приемниках. Настоящее изобретение предоставляет способы для осуществления этого. Haber и др. [25] также описывает подход для задачи перемещения приемников в инверсии одновременных кодированных источников с использованием способа стохастической оптимизации и применяет его к задаче удельного сопротивления по постоянному току.

Сущность изобретения

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

В более конкретном варианте осуществления изобретение представляет собой машинореализованный способ для инверсии полного волнового поля измеренных геофизических данных, чтобы определять модель физических свойств для области геологической среды, причем упомянутый способ использует компьютер для того, чтобы выполнять этапы, содержащие: (a) составление начальной модели физических свойств и использование ее для того, чтобы моделировать синтетические данные, соответствующие измеренным геофизическим данным, при этом источники и/или приемники при моделировании кодируются и моделируются одновременно, тем самым формируя моделированные кодированные данные; (b) разделение моделированных кодированных данных согласно источнику, что приводит к моделированным данным последовательных источников; (c) обработку моделированных данных последовательных источников таким образом, что обеспечивается их соответствие измеренным геофизическим данным для комбинаций источников и приемников, которые отсутствуют в измеренных геофизических данных; (d) вычисление разностей, называемых остатками данных, между обработанными моделированными данными последовательных источников и измеренными геофизическими данными для комбинаций источников и приемников, которые присутствуют в измеренных геофизических данных; (e) кодирование остатков данных с использованием кодирования, идентичного кодированию, используемому в (a), или другого кодирования; (f) использование кодированных остатков данных, чтобы вычислять обновление начальной модели физических свойств; и (g) комбинирование обновления с начальной моделью, чтобы формировать обновленную модель физических свойств. Вследствие линейности порядок этапов (d) и (e) может меняться местами.

В вариации этого последнего варианта осуществления следующие этапы могут быть выполнены в любое время перед этапом (c) вычисления разностей: кодирование измеренных геофизических данных с использованием кодирования источников, идентичного или различного относительно кодирования источников, используемого при формировании моделированных кодированных данных; затем разделение кодированных измеренных данных согласно источнику, предпочтительно с использованием идентичного алгоритма разделения источников, который использован при разделении моделированных кодированных данных согласно источнику; и использование разделенных измеренных данных при вычислении остатков данных.

Во втором более конкретном варианте осуществления настоящее изобретение представляет собой машинореализованный способ для инверсии полного волнового поля измеренных геофизических данных, чтобы определять модель физических свойств для области геологической среды, причем упомянутый способ использует компьютер для того, чтобы выполнять этапы, содержащие: (a) составление начальной модели физических свойств и использование ее для того, чтобы моделировать синтетические данные, соответствующие измеренным геофизическим данным, при этом источники и/или приемники при моделировании кодируются с использованием выбранных функций кодирования и моделируются одновременно, тем самым формируя моделированные кодированные данные; (b) разделение моделированных кодированных данных согласно источнику и/или приемнику с использованием алгоритма разделения источников или приемников, что приводит к моделированным данным последовательных источников и/или приемников; (c) обработку моделированных данных последовательных источников и/или приемников, чтобы подавлять эти данные, называемые отсутствующими данными, соответствующие комбинациям источников и приемников, которые отсутствуют; (d) кодирование отсутствующих данных с использованием функций кодирования, которые являются идентичными или отличающимися относительно упомянутых выбранных функций кодирования; (e) кодирование измеренных геофизических данных с использованием функций кодирования, которые являются идентичными функциям кодирования, используемым при кодировании отсутствующих данных, и вычитание их согласно позиции источников и приемников из моделированных кодированных данных, тем самым формируя кодированные остатки данных, которые являются некорректными вследствие отсутствующих комбинаций источников и приемников в измеренных геофизических данных; (f) вычитание кодированных отсутствующих данных из кодированных остатков данных, чтобы формировать скорректированные кодированные остатки данных; (g) использование скорректированных кодированных остатков данных, чтобы вычислять обновление начальной модели физических свойств; и (h) комбинирование обновления с начальной моделью, чтобы формировать обновленную модель физических свойств.

В третьем более конкретном варианте осуществления настоящее изобретение представляет собой машинореализованный способ для инверсии полного волнового поля измеренных геофизических данных, чтобы определять модель физических свойств для области геологической среды, причем упомянутый способ использует компьютер для того, чтобы выполнять этапы, содержащие: (a) составление начальной модели физических свойств и использование ее для того, чтобы моделировать синтетические данные, соответствующие измеренным геофизическим данным, при этом источники и/или приемники при моделировании кодируются с использованием выбранных функций кодирования и моделируются одновременно, тем самым формируя моделированные кодированные данные; (b) использование упомянутых выбранных функций кодирования, чтобы кодировать измеренные геофизические данные, и вычитание их согласно позиции источников и приемников из моделированных кодированных данных, тем самым формируя кодированные остатки данных, которые включают в себя некорректные значения вследствие отсутствующих комбинаций источников и приемников в измеренных геофизических данных; (c) разделение кодированных остатков данных согласно источнику и/или приемнику с использованием алгоритма разделения источников или приемников, что приводит к последовательным остаточным данным; (d) обработку последовательных остаточных данных, чтобы оценивать моделированные данные, называемые отсутствующими данными, соответствующие комбинациям источников и приемников, которые отсутствуют в измеренных геофизических данных; (e) кодирование отсутствующих данных с использованием упомянутых выбранных функций кодирования; (f) вычитание кодированных отсутствующих данных из кодированных остатков данных, чтобы формировать скорректированные кодированные остатки данных; (g) использование скорректированных кодированных остатков данных, чтобы вычислять обновление начальной модели физических свойств; и (h) комбинирование обновления с начальной моделью, чтобы формировать обновленную модель физических свойств.

Краткое описание чертежей

Вследствие ограничений на основе правил выдачи патентов на использование цвета прилагаемые чертежи являются черно-белыми репродукциями цветных оригиналов. Копии американской взаимодополняющей заявки или публикации с 20 цветными чертежами могут быть получены из Патентного ведомства (США) после запроса и выплаты необходимых платежей.

Настоящее изобретение и его преимущества должны лучше пониматься посредством рассмотрения нижеприведенного подробного описания и прилагаемых чертежей, на которых:

фиг.1 является блок-схемой последовательности операций способа, показывающей базовые этапы в настоящем изобретаемом способе для инверсии одновременных источников для задачи геометрии нестационарных приемников в варианте осуществления, в котором кодированные прогнозированные данные разделяются с использованием способа разделения источников;

фиг.2 является блок-схемой последовательности операций способа, показывающей базовые этапы в настоящем изобретаемом способе для инверсии одновременных источников для задачи геометрии нестационарных приемников в варианте осуществления, в котором кодированные прогнозированные данные, а также кодированные измеренные данные разделяются с использованием способа разделения источников, чтобы уменьшать ошибку в процессе разделения источников; и

фиг.3 и 4 являются блок-схемами последовательности операций способа, показывающими базовые этапы в настоящем изобретаемом способе для инверсии одновременных источников для задачи геометрии нестационарных приемников в двух вариантах осуществления, в которых некорректные кодированные остатки данных корректируются посредством кодированной части отсутствующих данных, полученных из способа разделения источников.

Изобретение описывается в связи с примерными вариантами осуществления. Тем не менее, в степени, в которой нижеприведенное подробное описание является характерным для конкретного варианта осуществления или конкретного варианта применения изобретения, оно должно быть только иллюстративным и не долж