Автоматическое выделение поверхностей для построения геолого-гидродинамической модели нефтегазового месторождения по сейсмическим данным

Авторы патента:

Ческис Ольга Андреевна (RU)

Трегуб Семен Леонидович (RU)

Казаров Андрей Сергеевич (RU)

G06F17/15 - вычисление корреляционных функций

Владельцы патента RU 2515081:

Общество с ограниченной ответственностью "Рок Флоу Динамикс" (RU)

Изобретение относится к области моделирования нефтегазовых месторождений. Сущность: определяют решетку модели с заданным шагом и предельно допустимой величиной на координатной сетке (1). Задают управляющие параметры поверхностей относительно шага решетки (2). Определяют набор поверхностей решетки относительно ее шага (3). Определяют область, в которой определена поверхность, и добавляют узлы решетки, в которых поверхность не определена (4). Задают параметры кубов сейсмических атрибутов - p, q, задают координаты вновь добавленного узла - X, Y, задают координаты соседнего узла - x, y, в котором поверхность определена, и задают параметр глубины поверхности в этом узле - z (x, y) (5). Упорядочивают глубины поверхностей в указанном узле решетки по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине (6). Определяют поверхности для обработки (7). Определяют функционалы и градиенты поверхностей для обработки (8). Определяют величину текущего шага смещения вдоль градиента (9). Определяют для каждой обрабатываемой поверхности с помощью величины шага смещения градиента, для какого по счету шага достигается минимальное значение функционала (10). Когда для каждой обрабатываемой поверхности минимальное значение функционала меньше шага смещения градиента, то производят смещение вдоль градиента на найденное количество шагов (11). Сохраняют в памяти величину сделанного смещения и определяют максимум сделанных смещений (12). Исключают временно из обработки поверхности, для которых достигнуто уменьшение функционала (13). Определяют поверхности, для которых не удалось уменьшить функционал, и, когда такие поверхности определены, то повторяют этапы 7-13 (14). Повторно упорядочивают глубины определенных на указанном узле решетки поверхностей по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине (15). Определяют значения текущих шагов решетки по обновленным данным поверхностей (16). Итеративно повторяют все этапы способа до тех пор, пока не будет получена модель с заданными на этапах 1 и 2 параметрами (17). Технический результат: повышение точности построения геолого-гидродинамической модели нефтегазового месторождения. 3 н.п.ф-лы, 8 з.п.ф-лы, 4 ил.

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

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

Известен способ построения трехмерных геологических структур на основе данных поверхностей (G. Caumon et al., "Surface-Based 3D Modeling of Geological Structures ", 23.09.2009), который описывает тесную интеграцию данных поверхностей при геологическом моделировании. В способе выявляют узлы поверхности границ решетки, которые пересекаются с другими поверхностями. Используют ограничения толщины и диапазон при выполнении интерполирования, чтобы смоделировать поверхности, лежащие на заданном расстоянии друг от друга. Это расстояние вычисляется на векторном поле.

Известен способ геофизической разведки для выявления нефтегазовых объектов (RU 2289829 C1, G01V 11/00, 20.12.2006), который раскрывает оптимизацию размещения глубоких скважин на нефтегазовых объектах по комплексу данных наземной сейсмической, электро-, магнито-, гравиразведки; электрического, радиоактивного, акустического, сейсмического, магнитного, гравитационного каротажа; изучения керна и испытания скважин. По данным акустического, сейсмического, электрического, радиоактивного, магнитного и гравитационного каротажа, лабораторных исследований керна формируют жесткостные, электрические, магнитные и гравиметрические модели целевого интервала геологического разреза в скважинах, рассчитывают геофизические синтетические трассы, по которым проводят СВАН, определяют модельные СВО и их спектрально-временные атрибуты (СВА). Спектрально-временные атрибуты (СВА) представляют собой отношение энергии высоких частот и больших времен к энергии низких частот и малых времен, а также произведения удельной спектральной плотности на средневзвешенные частоту и время либо на максимальные частоту и время энергетических спектров СВАН-колонки по оси частот и времен.

Известны также различные методы построения цифровых геологических моделей нефтегазовых месторождений и проектирование разработки с использованием цифровых геологических моделей, предлагаемых компанией НПЦ «Геонефтегаз» (http://www.geoneftegaz.ru/models/mod.htm). Компания предлагает решение по моделированию цифровой геологической модели, построенной с использованием технологии, позволяющей изучать коллекторские свойства и нефтеносность карбонатных и терригенных пластов в сложнопостроенных резервуарах с тектоническим и литологическим экранированием, что позволяет вести подсчет запасов нефти по промышленным категориям. Данное решение используется для расчета геолого-технологических гидродинамических моделей и оптимизации проектов разработки. Применение данной технологии прогноза коллекторских свойств и нефтеносности позволяет обоснованно строить трехмерные геологические модели при минимальном объеме глубокого бурения. Их вовлечение в процесс расчета и анализа технологических показателей разработки позволит сформировать оптимальные техсхемы и ускорить ввод новых объектов в эксплуатацию. Анализ подтверждаемости прогнозных параметров пласта - эффективной толщины, пористости, типа

флюида с данными бурения новых скважин - свидетельствует об удовлетворительной сходимости результатов.

Известен способ контроля геометрических и гидродинамических параметров гидроразрыва пласта (RU 2390805 C1, G01V 5/12, 27.05.2012), который заключается в том, что перед проведением гидроразрыва пласта на поверхности в квадрате 400×400 м с центром, совпадающим с устьем скважины, проводят определение концентрации радона в шурфах глубиной 50-70 см по сетке с шагом 50 м. По данным поверхностной съемки проводят линии геодинамических зон, которые могут быть связаны со вновь образующимися трещинными системами в пласте в результате гидроразрыва (так как развиваются большие давления). В скважине перед гидроразрывом проводят радоновые индикаторные исследования и определяют ее техническое состояние: заколонные перетоки, если таковые имеются, их направление, долевое распределение по ним закачиваемой жидкости, профиль приемистости коллектора по данным трех временных замеров и проницаемость.

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

Раскрытие изобретения

Задачей настоящего изобретения является обеспечение способа, устройства и машиночитаемого носителя данных, предназначенных для выполнения этапов, составляющих существенную часть всей задачи построения геолого-гидродинамической модели, а именно автоматического (т.е. требующего участия пользователя только на этапе ввода начальных данных) построения набора поверхностей, основанного на входных сейсмических данных. Полученные на выходе поверхности можно без дополнительной обработки использовать в построении геолого-гидродинамической решетки. Входными данными для способа являются кубы сейсмических атрибутов Р, Q; сейсмический куб (опционально).

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

В одном варианте осуществления изобретения способ включает в себя этапы на которых: 1) определяют решетку модели с заданным шагом и предельно допустимой величиной на координатной сетке; 2) задают управляющие параметры поверхностей относительно шага решетки; 3) определяют набор поверхностей решетки относительно ее шага; 4) определяют область, в которой определена поверхность, и добавляют узлы решетки, в которых поверхность не определена; 5) задают параметры кубов сейсмических атрибутов - p, q, задают координаты вновь добавленного узла - X, Y, задают координаты соседнего узла - х, у, в котором поверхность определена, и задают параметр глубины поверхности в этом узле - z( x, y); 6) упорядочивают глубины поверхностей в указанном узле решетки по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине; 7) определяют поверхности для обработки; 8) определяют функционалы и градиенты поверхностей для обработки; 9) определяют величину текущего шага смещения вдоль градиента; 10) определяют для каждой обрабатываемой поверхности с помощью величины шага смещения градиента, для какого по счету шага достигается минимальное значение функционала; причем 11) когда для каждой обрабатываемой поверхности минимальное значение функционала меньше шага смещения градиента, то производят смещение вдоль градиента на найденное количество шагов; 12) сохраняют в памяти величину сделанного смещения и определяют максимум сделанных смещений; 13) исключают временно из обработки поверхности, для которых достигнуто уменьшение функционала; 14) определяют поверхности, для которых не удалось уменьшить функционал, и, когда такие поверхности определены, то повторяют этапы 7) - 13); 15), повторно упорядочивают глубины определенных на указанном узле решетки поверхностей по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине; 16) определяют значения текущих шагов решетки по обновленным данным поверхностей; 17) итеративно повторяют этапы 3-17 способа до тех пор, пока не будет получена модель с заданными на этапах 1 и 2 параметрами.

Во втором варианте настоящее изобретение обеспечивает устройство для выполнения способа по способу, описанному в первом варианте осуществления изобретения, и который, соответственно, дополнительно не описывается. Устройство может представлять собой, но не ограничиваться: суперкомпьютер, персональный компьютер, портативный компьютер, планшетный компьютер, карманный компьютер, смартфон и тому подобное. Устройство обязательно содержит один или более процессоров, которые предназначены для выполнения компьютерных команд или кодов, хранящихся в памяти устройства с целью обеспечения выполнения способов по первому или второму вариантам осуществления настоящего изобретения, машиночитаемый носитель данных (память) и модули ввода/вывода (I/O). Модули I/O представляют собой, но не ограничиваются, типичные и известные из уровня техники средства управления устройством: манипулятор типа «мышь», клавиатура, джойстик, тачпад, трекбол, электронное перо, стилус, сенсорный дисплей и тому подобное. Так же модули I/O представляют собой, но не ограничиваются типичные и известные из уровня техники средства демонстрирования информации: монитор, проектор, принтер, графопостроитель и тому подобное. В качестве примера, а не ограничения, машиночитаемый носитель данных может включать в себя оперативную память (RAM); постоянное запоминающее устройство (ROM); электрически стираемое программируемое постоянное запоминающее устройство (EEPROM); флэш-память или другие технологии памяти; CDROM, цифровой универсальный диск (DVD) или другие оптические или голографические носители данных; магнитные кассеты, магнитную пленку, запоминающее устройство на магнитных дисках или другие магнитные запоминающие устройства, несущие волны или другой носитель данных, который может быть использован для кодирования требуемой информации и к которому может быть осуществлен доступ посредством устройства, описанного выше.

В третьем варианте настоящее изобретение обеспечивает машиночитаемый носитель данных, содержащий код программы, который побуждает процессор и/или процессоры выполнять действия по способу, описанному в первом варианте осуществления изобретения, и который, соответственно, дополнительно не описывается. В качестве примера, а не ограничения, машиночитаемый носитель данных может включать в себя оперативную память (RAM); постоянное запоминающее устройство (ROM); электрически стираемое программируемое постоянное запоминающее устройство (EEPROM); флэш-память или другие технологии памяти; CDROM, цифровой универсальный диск (DVD) или другие оптические или голографические носители данных; магнитные кассеты, магнитную пленку, запоминающее устройство на магнитных дисках или другие магнитные запоминающие устройства, несущие волны или другой носитель данных, который может быть использован для кодирования требуемой информации, и к нему может быть осуществлен доступ посредством устройства, описанного во втором варианте осуществления изобретения и, которое, соответственно дополнительно не описывается.

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

Фиг.1 иллюстрирует блок-схему выполнения этапов заявленного способа.

Фиг.2-4 иллюстрируют результаты выполнения итераций способа при построении геолого-гидродинамической модели.

Детальное описание чертежей

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

На фиг.1 изображена блок-схема выполнения этапов заявленного способа, а именно этапы: этап 1, на котором определяют решетку модели с заданным шагом и предельно допустимой величиной на координатной сетке; этап 2, на котором задают управляющие параметры поверхностей относительно шага решетки; этап 3, на котором определяют набор поверхностей решетки относительно ее шага; этап 4, на котором определяют область, в которой определена поверхность и добавляют узлы решетки, в которых поверхность не определена; этап 5, на котором задают параметры кубов сейсмических атрибутов - p, q, задают координаты вновь добавленного узла - X, Y, задают координаты соседнего узла - х, у, в котором поверхность определена, и задают параметр глубины поверхности в этом узле - z(x,y); этап 6, на котором упорядочивают глубины поверхностей в указанном узле решетки по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине; этап 7, на котором определяют поверхности для обработки; этап 8, на котором определяют функционалы и градиенты поверхностей для обработки; этап 9, на котором определяют величину текущего шага смещения вдоль градиента; этап 10, на котором определяют для каждой обрабатываемой поверхности с помощью величины шага смещения градиента, для какого по счету шага достигается минимальное значение функционала; причем на этапе 11, если для каждой обрабатываемой поверхности минимальное значение функционала меньше шага смещения градиента, то производят смещение вдоль градиента на найденное количество шагов; этап 12, на котором сохраняют в памяти величину сделанного смещения и определяют максимум сделанных смещений; этап 13, на котором исключают временно из обработки поверхности, для которых достигнуто уменьшение функционала; этап 14, на котором определяют поверхности, для которых не удалось уменьшить функционал, и, когда такие поверхности определены, то повторяют этапы 7-13; этап 15, повторно упорядочивают глубины определенных на указанном узле решетки поверхностей по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине; этап 16, на котором определяют значения текущих шагов решетки по обновленным данным поверхностей; и этап 17, на котором итеративно повторяют этапы 3-17 способа до тех пор, пока не будет получена модель с заданными на этапах 1 и 2 параметрами.

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

Успех реализации способа обеспечивается следующими критериями:

- одновременное (параллельное) выделение и последующее уточнение достаточно большого количества поверхностей;

- семейство функционалов, минимизация которых составляет основу успешного выполнения способа;

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

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

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

Настоящий способ предлагает использовать не собственно сейсмический куб, а производные данные, так называемые сейсмические атрибуты Р и Q. Эти атрибуты для каждой точки каждой трассы показывают, под каким углом находится максимально похожая на нее точка соседней трассы. Тем самым можно считать, что эти атрибуты задают поле касательных плоскостей к отражающим поверхностям. А сама задача сводится к построению поверхностей, для которых это поле плоскостей является касательным. Это делается с помощью минимизации функционала. Задается функционал, который измеряет суммарное отличие поля плоскостей, задаваемого Р, Q атрибутами и от касательных плоскостей к поверхностям набора. Метод градиентного спуска позволяет перейти к новому набору поверхностей, для которого функционал имеет меньшее значение, тем самым лучше согласованному с входными данными.

Способ в общем включает в себя этапы, использующиеся в каждой итерации выполнения способа:

1. «Растекание». Решает задачу продолжения поверхности в те области, где он не определен. Используется, в частности, при построении начальных приближений.

2. «Уточнение». Построение интегральной поверхности по полю (касательных) плоскостей, задаваемому Р, Q кубами, на основе минимизации функционала.

3. «Сортировка». Перестройка поверхностей по глубине. Решает задачу исключения пересечений поверхностей.

4. «Подтягивание» к экстремумам (максимумам или минимума). Далее приведено детальное раскрытие выполнения этапов.

Этап «Растекание»

Вход: набор поверхностей. Поверхности определены не на всей решетке, р, q - кубы сейсмических атрибутов.

Выход: набор поверхностей с большей областью определения. Выполнять в цикле:

1. Для каждой поверхности добавить узлы решетки, в которых поверхность не определена, но определена в одном из соседних узлов, к области определения этой поверхности.

2. Если добавить узлы ни для одной поверхности не удалось, выйти из цикла.

3. Для каждого из вновь добавленных узлов: А) Рассчитать значение по формуле:

$Z = \frac{\sum (z (x, y) + (X - x)) \times p (x, y, z) + (Y - y) \times q (x, y, z) / \sqrt{{(X - x)}^{2} + {(Y - y)}^{2}}}{\sum 1 / \sqrt{{(X - x)}^{2} + {(Y - y)}^{2}}}$

Здесь p, q - кубы сейсмических атрибутов, X, Y- координаты вновь добавленного узла, x, y - координаты соседнего узла, в котором поверхность определена, a z (x, y) - глубина поверхности в этом узле. Сумма распространяется по всем соседним узлам, в которых поверхность определена.

Б) Проверить, что таким образом расширенная поверхность попадает в требуемую область пространства.

Смысл указанной выше формулы следующий. Для каждого соседнего узла выражение z(x,y)+(X-x)×р(x,y,z)+(Y-y)×q(x,y,z) - это линейное приближение по формуле Тейлора значения в новом узле по значению в соседнем. Эти линейные приближения для всех соседних узлов далее усредняются с весами. Веса при этом обратно пропорциональны расстояниям между добавленным и соседним узлами.

4. Если рассчитать значение ни в одном добавленном узле не удалось, значит поверхность определена в каждом узле решетки.

Этап «Уточнение»

Вход: набор поверхностей - начальное или стартовое приближение, р, q - кубы сейсмических атрибутов.

Выполняется минимизация функционала, который представляет собой сумму квадратичных отклонений частных производных по x, y уравнений, задающих поверхности, от соответственно значений P и Q атрибутов:

$F = {\sum [\frac{\partial z}{\partial x} (x_{i}, y_{j}) - p (x_{i}, y_{j}, z_{i, j})]}^{2} + {\sum [\frac{\partial z}{\partial y} (x_{i}, y_{j}) - q (x_{i}, y_{j}, z_{i, j})]}^{2}$

Здесь суммирование происходит по всем узлам и всем поверхностям набора.

Минимизация функционала осуществляется методом наискорейшего градиентного спуска, соответствующим образом оптимизированного и адаптированного к ситуации.

Схема этапа:

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

2. Вычислить функционал.

3. Вычислить градиент.

4. Далее в цикле:

5. Вычислить величину текущего шага смещения вдоль градиента.

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

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

8. Запомнить величину сделанного смещения и вычислить максимум сделанных смещений.

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

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

11. Если не для всех, но величина текущего шага смещения достаточно мала, то также происходит выход из цикла.

12. В противном случае переход к началу цикла. (На следующем шаге цикла будут обрабатываться только те поверхности, функционал для которых не удалось уменьшить).

Этап «Сортировка»

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

Схема:

Цикл по всем узлам:

1. Упорядочить глубины поверхностей в этом узле по возрастанию.

2. Назначить глубину i-й поверхности в этом узле равной i-Pi по порядку глубине.

Этап «Подтягивание»

Простая процедура, в которой каждая точка каждой поверхности в определенной окрестности по глубине перемещается к локальному минимуму или максимуму сейсмического куба. Если в заданной окрестности экстремум найти не удалось, точка выбрасывается. Таким образом, получаем, более согласованный с сейсмическими данными набор поверхностей. Данная процедура полезна и для нейтрализации типичных проблем применения метода градиентного спуска.

Схема:

Цикл по всем узлам:

Цикл по набору поверхностей

1. Взять соответствующую узлу точку поверхности.

2. Взять отрезок соответствующей трассы сейсмического куба вокруг этой точки и найти в нем ближайший к точке нужный (максимум или минимум) локальный экстремум.

3. Если найти локальный экстремум удалось, то изменить глубину поверхности в этом узле на глубину локального экстремума.

4. Если найти локальный экстремум не удалось - в поверхности появляется дырка - в этом узле поверхность считается неопределенной.

Описание основных действий при выполнении способа автоматического выделения поверхностей для построения геолого-гидродинамической модели нефтегазового месторождения по сейсмическим данным:

Вход: набор трендовых поверхностей - начальное или стартовое приближение. При их отсутствии берутся константные по глубине поверхности.

Р, Q - кубы сейсмических атрибутов

опционально: сейсмический куб

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

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

Схема выполнения способа:

Выполняется в цикле:

Этап 1. Установить управляющие параметры алгоритма согласно текущим значениям шагов решетки, (поверхности).

Этап 2. Пересчитать имеющийся текущий набор поверхностей на решетку с текущими шагами.

Этап 3. Если задан сейсмический куб, то выполнить процедуру «подтягивания» к экстремумам.

Этап 4. Воспользоваться процедурой "растекания", чтобы распространить набор поверхностей на те точки решетки, где они не заданы.

Этап 5. Устранить взаимные пересечения поверхностей из набора процедурой «сортировки».

Этап 6. Применить процедуру «уточнения» набора поверхностей путем минимизации функционала методом градиентного спуска.

Этап 7. Устранить взаимные пересечения поверхностей из набора процедурой «сортировки».

Этап 8. Если решетка достигла максимального размера, то выход из цикла.

Этап 9. Пересчитать значения текущих шагов решетки для нового круга цикла.

Этап 10. Перейти к началу цикла.

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

1. Способ автоматического выделения поверхностей для построения геолого-гидродинамической модели нефтегазового месторождения по сейсмическим данным, заключающийся в выполнении этапов, на которых:
1) определяют решетку модели с заданным шагом и предельно допустимой величиной на координатной сетке;
2) задают управляющие параметры поверхностей относительно шага решетки;
3) определяют набор поверхностей решетки относительно ее шага;
4) определяют область, в которой определена поверхность, и добавляют узлы решетки, в которых поверхность не определена;
5) задают параметры кубов сейсмических атрибутов - p, q, задают координаты вновь добавленного узла - X, Y, задают координаты соседнего узла - x, y, в котором поверхность определена, и задают параметр глубины поверхности в этом узле - z(x, y);
6) упорядочивают глубины поверхностей в указанном узле решетки по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине;
7) определяют поверхности для обработки;
8) определяют функционалы и градиенты поверхностей для обработки;
9) определяют величину текущего шага смещения вдоль градиента;
10) определяют для каждой обрабатываемой поверхности с помощью величины шага смещения градиента, для какого по счету шага достигается минимальное значение функционала; причем
11) когда для каждой обрабатываемой поверхности минимальное значение функционала меньше шага смещения градиента, то производят смещение вдоль градиента на найденное количество шагов;
12) сохраняют в памяти величину сделанного смещения и определяют максимум сделанных смещений;
13) исключают временно из обработки поверхности, для которых достигнуто уменьшение функционала;
14) определяют поверхности, для которых не удалось уменьшить функционал, и, когда такие поверхности определены, то повторяют этапы 7-13;
15) повторно упорядочивают глубины определенных на указанном узле решетки поверхностей по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине;
16) определяют значения текущих шагов решетки по обновленным данным поверхностей;
17) итеративно повторяют все этапы способа до тех пор, пока не будет получена модель с заданными на этапах 1 и 2 параметрами.

2. Способ по п.1, отличающийся тем, что управляющие параметры представляют собой сейсмический куб.

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

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

5. Способ по п.3, отличающийся тем, что локальный экстремум является максимальным значением глубины.

6. Устройство для автоматического выделения поверхностей для построения геолого-гидродинамической модели нефтегазового месторождения по сейсмическим данным, содержащее:
один или множество процессоров;
модули ввода/вывода (I/O);
машиночитаемый носитель данных (память), который содержит код программы, который при выполнении побуждает процессор или процессоры выполнять этапы, на которых:
1) определяется решетка модели с заданным шагом и предельно допустимой величиной на координатной сетке;
2) задаются управляющие параметры поверхностей относительно шага решетки;
3) определяется набор поверхностей решетки относительно ее шага;
4) определяется область, в которой определена поверхность, и добавляют узлы решетки, в которых поверхность не определена;
5) задаются параметры кубов сейсмических атрибутов - p,q, задают координаты вновь добавленного узла - X, Y, задают координаты соседнего узла - x, y, в котором поверхность определена, и задают параметр глубины поверхности в этом узле - z(x,y);
6) упорядочиваются глубины поверхностей в указанном узле решетки по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине;
7) определяются поверхности для обработки;
8) определяются функционалы и градиенты поверхностей для обработки;
9) определяются величины текущего шага смещения вдоль градиента;
10) определяются для каждой обрабатываемой поверхности с помощью величины шага смещения градиента, для какого по счету шага достигается минимальное значение функционала; причем
11) когда для каждой обрабатываемой поверхности минимальное значение функционала меньше шага смещения градиента, то производится смещение вдоль градиента на найденное количество шагов;
12) сохраняется в памяти величина сделанного смещения и определяется максимум сделанных смещений;
13) исключаются временно из обработки поверхности, для которых достигнуто уменьшение функционала;
14) определяются поверхности, для которых не удалось уменьшить функционал, и, когда такие поверхности определены, то повторяют этапы 7-13;
15) повторно упорядочиваются глубины определенных на указанном узле решетки поверхностей по возрастанию и назначают глубину i-й поверхности в этом узле равной i-й по порядку глубине;
16) определяются значения текущих шагов решетки по обновленным данным поверхностей;
17) итеративно повторяются все перечисленные этапы до тех пор, пока не будет получена модель с заданными на этапах 1 и 2 параметрами.

7. Устройство по п.6, отличающееся тем, что управляющие параметры представляют собой сейсмический куб.

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

9. Устройство по п.8, отличающееся тем, что локальный экстремум является минимальным значением глубины.

10. Устройство по п.8, отличающееся тем, что локальный экстремум является максимальным значением глубины.

11. Машиночитаемый носитель данных, пригодный для использования в устройстве по любому из п.п.6-10, который содержит код программы, который при выполнении побуждает процессор или процессоры выполнять действия способов по любому из п.п.1-5.

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

Корреляционный измеритель временных сдвигов случайных сигналов // 2500025

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

Способ и устройство для обработки сигнала и машиночитаемый носитель информации // 2495483

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

Способ корреляционной обработки смеси гармонического сигнала с шумом // 2490701

Изобретение относится к области обработки данных и может быть использовано для выделения гармонического сигнала на фоне помех и измерения его частоты. .

Способ передачи информационных сигналов и устройство для его осуществления // 2460211

Изобретение относится к способам кодирования, декодирования и преобразования кода для обнаружения и исправления ошибок. .

Корреляционный измеритель временных сдвигов // 2445690

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

Способ оценки частотных искажений и устройство для его осуществления // 2420798

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

Способ обработки сложных квазиоптимальных сигналов и устройство для осуществления способа // 2377644

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

Телевизионное следящее устройство корреляционного типа // 2364927

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

Способ контроля уровня получаемых знаний // 2338264

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

Способ определения положения координат глубин маркера при построении геологической модели месторождения // 2530324

Изобретение относится к способу, устройству и машиночитаемому носителю данных, предназначенным для построения геологической модели нефтяного или иного месторождения, в частности, для определения коэффициентов корреляции для комплекса кривых ГИС и нахождения положений глубин маркера, для которых значение коэффициента корреляции является максимальным. Техническим результатом является повышение точности вычислений параметров, используемых при построении геологической модели расположения нефтяных или иных месторождений. Метод позволяет для маркера, уже имеющего отметки на некоторой, называемой опорной, группе скважин, вычислить их для скважин из другой группы. Для каждой скважины W, на которой ищется значение глубины маркера, выбираются скважины опорной группы, отстоящие от скважины W на заданном расстоянии, и среди них выбирается скважина с наибольшим значением коэффициента корреляции, при этом точка, в которой этот максимум достигается, назначается искомой отметки маркера. С помощью проверяющих тестов осуществляют поиск скважин, в которых функция корреляции меньше, чем максимальное значение коэффициента корреляции, а коэффициент качества корреляции больше, чем максимальное значение коэффициента корреляции. После чего добавляют найденную скважину к опорной группе скважин. 3 н. и 5 з.п. ф-лы, 9 ил.

Способ, устройство и машиночитаемый носитель данных для определения положения координат глубин маркера при построении геологической модели месторождения // 2534964

Изобретение относится к способу, устройству и машиночитаемому носителю данных, используемых при построении геологической модели нефтяного или иного месторождения. Технический результат - повышение точности вычислений параметров, используемых при построении геологической модели расположения нефтяных или иных месторождений. Изобретение позволяет для маркеров, выбранных в качестве начального решения, вычислить такие глубины маркера на каждой скважине, которые обеспечивают наилучшую суммарную корреляцию. Для каждого маркера, входящего в набор, определяется функционал, представляющий собой сумму коэффициентов корреляции комплекса методов ГИС для пар скважин, расположенных не далее заданного расстояния друг от друга. Для этого функционала вычисляются частные производные ,и полученный таким образом вектор сглаживается и используется для нахождения большего значения функционала на некотором отрезке вдоль этого вектора. Если большего значения не найдено, то последнее положение отметок маркера считается решением задачи, а если найдено, то производится сглаживание точки решения и процесс повторяется снова. На каждой итерации алгоритма производится сортировка глубин маркеров. 3 н.п. ф-лы, 9 ил.

Корреляционное устройство // 2536393

Изобретение относится к корреляционному устройству. Технический результат заключается в повышении быстродействия определения времени задержки. Корреляционное устройство содержит два регистра сдвига, блок сравнения, два триггера, счетчик, элемент ИСКЛЮЧАЮЩЕЕ ИЛИ, элемент И, генератор тактовых импульсов и двухканальный трехразрядный мультиплексор. 1 ил.

Мультиплексирующий цифровой коррелятор // 2540833

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

Способ определения границ природных очагов биогельминтозов // 2545707

Изобретение относится к области медицины, в частности эпидемиологии, и предназначено для определения границ природных очагов биогельминтозов с использованием генетических маркеров. На обследуемых территориях отлавливают рыбу из различных районов речного бассейна. Из выловленной рыбы производят забор мышечной ткани в области спинного плавника. Мышечную ткань гомогенезируют до получения однородной массы, отбирают супернатант, проводят разделение белков методом вертикального электрофореза и выявляют полиморфные белковые фракции. По каждой выявленной полиморфной белковой фракции проводят расчет популяционно-генетических параметров и определяют генетический маркер и показатель индекса генетического подобия (или сходства Нея), одновременно в исследуемых особях выявляют наличие гельминтов, а также определяют уровень инвазии. По сходству и различию популяционно-генетических параметров в пределах обследуемых территорий отслеживают возможные границы природных очагов биогильментозов. Изобретение позволяет эффективно определять границы природных очагов биогельминтозов и осуществлять противоэпидемические мероприятия среди населения. 8 з.п. ф-лы, 1 ил., 3 пр.

Способ обработки функции автокорреляции для измерения основного тона речевого сигнала // 2559710

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

Вычислитель доплеровской фазы пассивных помех // 2559750

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

Вычислитель для адаптивного режектирования помех // 2582871

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

Вычислитель радиальной скорости движущегося объекта // 2613037

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

Коррелятор // 2618941

Изобретение относится к радиоэлектронной технике и может быть использовано в системах передачи информации для определения идентичности кода принимаемого от передатчика сигнала коду, хранящемуся в памяти приемника. Техническим результатом является построение коррелятора для квазигармонических импульсов. Устройство содержит входную информационную шину входного сигнала, содержащего смесь полезного сигнала, помех и шума с квазигармонической формой импульса, АЦП, процессор сравнения, счетчик совпадений, регистр сдвига, ПЗУ, синхронизатор, схему разрешения. 2 ил.