СПОСОБЫ ДЛЯ АППРОКСИМАЦИИ ОПЕРАЦИИ УМНОЖЕНИЯ ГЕССИАНА НА ВЕКТОР В ПОЛНОЙ ИНВЕРСИИ ВОЛНОВОГО ПОЛЯ

1. Способ для определения дискретной модели физических свойств области недр, которая может упоминаться как модель, посредством итеративного инвертирования измеренных геофизических данных, полученных из области недр, содержащий:
использование матрицы Гессе целевой функции, которая затем умножается на вектор, что называется «умножение гессиана на вектор», чтобы определять обновление для модели, при этом умножение гессиана на вектор аппроксимируется, с использованием компьютера, с помощью одиночного смоделированного распространения прямой волны и одиночного вычисления градиента целевой функции, в модифицированной модели недр.
2. Способ по п. 1, в котором модифицированная модель недр является линейной комбинацией текущей модели недр и упомянутого вектора.
3. Способ по п. 2, в котором линейная комбинация = текущая модель недр + {(скалярная константа) x (вектор)}, при этом 0< скалярная константа .
4. Способ по п. 3, в котором скалярная константа выбирается согласно тому, будет ли матрица Гессе вычисляться точно, т.е. с учетом многократного рассеяния, или будет ли аппроксимироваться посредством матрицы Гессе Гаусса-Ньютона.
5. Способ по п. 4, в котором, если должна использоваться аппроксимация Гаусса-Ньютона, скалярная константа выбирается так, что {(скалярная константа) x (вектор)} равняется приблизительно 1% в величине текущей модели недр; но если матрица Гессе должна вычисляться точно, то скалярная константа выбирается так, что {(скалярная константа) x (вектор)} равняется приблизительно 10% в величине текущей модели недр.
6. Способ по п. 3, в котором вектор в первой итерации способа является градиентом целевой функции по отношению к параметрам модели, при этом целевая функция измеряет несоответствие между измеренными геофизическими данными и соответствующими синтетическими геофизическими данными, смоделированными с использованием модели, при этом упомянутое несоответствие является выбранной нормой остатка данных, при этом остаток данных равняется разности между измеренными геофизическими данными и соответствующими синтетическими геофизическими данными, смоделированными с использованием модели; при условии, что для упомянутого одиночного вычисления градиента целевой функции, измеренные геофизические данные в остатке данных заменяются на синтетические данные, смоделированные с использованием модифицированной модели недр, и остаток данных делится на скалярную константу.
7. Способ по п. 1, в котором измеренные геофизические данные, которые обращаются, являются полным волновым полем сейсмических данных.
8. Способ по п. 7, в котором модель физических свойств является моделью, по меньшей мере, одной из группы, состоящей из 21 упругих констант, включающих в себя модуль объемной упругости; плотности; или любой комбинации двух или более из них.
9. Способ для определения дискретной модели физических свойств области недр посредством итеративного инвертирования измеренных геофизических данных, полученных из области недр, содержащий использование матрицы Гессе, соответствующей целевой функции, которая затем умножается на вектор, что называется «умножение гессиана на вектор», чтобы определять обновление для начальной модели, при этом умножение гессиана на вектор аппроксимируется, с использованием компьютера, с помощью одиночного смоделированного распространения прямой волны и одиночного вычисления градиента целевой функции в модифицированной модели недр, содержащий:
формирование масштабированного вектора посредством умножения вектора на число ε, где 0<ε≤1, и добавление масштабированного вектора к начальной модели;
моделирование первого синтетического набора данных с использованием модели с масштабированным вектором, добавленным к ней;
формирование целевой функции, чтобы измерять несоответствие между вторым синтетическим набором данных, смоделированным с использованием начальной модели, и упомянутым первым синтетическим набором данных;
вычисление градиента сформированной целевой функции по отношению к параметрам физических свойств, деление градиента на ε, и использование этого в качестве оценки умножения гессиана на вектор.
10. Способ по п. 9, в котором вектор в первой итерации итеративной инверсии является градиентом в пространстве модели целевой функции, при этом целевая функция измеряет несоответствие между измеренными геофизическими данными и соответствующими смоделированными по модели геофизическими данными.
11. Компьютерный программный продукт, содержащий нетранзиторный компьютерно-используемый носитель, имеющий машиночитаемый программный код, воплощенный там, при этом упомянутый машиночитаемый программный код выполнен с возможностью для исполнения, чтобы осуществлять способ для определения дискретной модели физических свойств области недр, которая может упоминаться как модель, посредством итеративного инвертирования измеренных геофизических данных, полученных из области недр, при этом упомянутый способ содержит:
использование матрицы Гессе целевой функции, которая затем умножается на вектор, что называется «умножение гессиана на вектор», чтобы определять обновление для модели, при этом умножение гессиана на вектор аппроксимируется с помощью одиночного смоделированного распространения прямой волны и одиночного вычисления градиента целевой функции, в модифицированной модели недр.
12. Компьютерный программный продукт по п. 11, в котором вектор в первой итерации итеративной инверсии является градиентом в пространстве модели целевой функции, при этом целевая функция измеряет несоответствие между измеренными геофизическими данными и соответствующими смоделированными по модели геофизическими данными.