МИНОБРНАУКИ РОССИИ
Федеральное государственное бюджетное
образовательное учреждение высшего образования
«Ухтинский государственный технический университет»
(УГТУ)
На правах рукописи
Кожевникова Полина Валерьевна
РАЗРАБОТКА МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ
ИНФОРМАЦИОННО-ЭКСПЕРТНОЙ ОЦЕНКИ АНАЛИЗА
ДОСТОВЕРНОСТИ ПРОМЫСЛОВО-ГЕОЛОГИЧЕСКИХ МОДЕЛЕЙ
НА ОСНОВЕ НЕЧЕТКОГО МОДЕЛИРОВАНИЯ
05.13.18 – Математическое моделирование, численные методы
и комплексы программ
Диссертация на соискание ученой степени кандидата технических наук
Научный руководитель:
доктор физико-математических наук,
профессор, Кобрунов А.И.
Ухта – 2018
Оглавление
Введение ....................................................................................................................... 4
Глава 1. Петрофизическое моделирование ............................................................. 9
1.1.
К
формированию
аналитических
зависимостей
различных
петрофизических моделей .......................................................................................... 9
1.2.
Проблемы, возникающие при моделировании петрофизических
сред
..................................................................................................................... 21
1.3.
Классификация петрофизических моделей ............................................ 22
1.4.
Корреляционно-регрессионные петрофизические модели .................. 24
1.5.
Заключение ................................................................................................ 29
Глава 2. Фазификация нечетких отношений, основанная на алгоритме
сжатия информации .................................................................................................. 31
2.1.
Нечеткое моделирование ......................................................................... 31
2.2.
Формирование функций принадлежности ............................................. 35
2.2.1
Алгоритм сжатия информации ................................................................ 37
2.2.2
Выбор аналитического вида функции принадлежности ...................... 41
2.2.3
Аппроксимация многомерных связей..................................................... 46
2.2.4
Описание алгоритма конструирования функции принадлежности ..... 49
2.3.
Заключение ................................................................................................ 55
Глава 3. Анализ нечетких моделей ........................................................................ 56
3.1.
Теоретическая основа анализа функций принадлежности ................... 56
3.2.
Алгоритм анализа функций принадлежности........................................ 57
3.3.
Экспериментальная основа 𝜶-анализа .................................................... 59
3.4.
Заключение ................................................................................................ 67
Глава 4. Нечеткий логический вывод Мамдани ................................................... 68
2
4.1.
Описание алгоритма нечеткого логического вывода ............................ 68
4.2.
Свойства нечеткого логического вывода ............................................... 72
4.2.1
Коммутативность ...................................................................................... 72
4.2.2
Возрастание неопределенности. Аналог подстановки уравнений ...... 77
4.3.
Прогнозирование емкостных параметров .............................................. 80
4.4.
Альтернативные виды композиций ........................................................ 91
4.5.
Заключение ................................................................................................ 95
Глава 5. Информационная экспертиза ................................................................... 97
5.1.
Метод пространственного продолжения ................................................ 98
5.2.
Дефазификация ....................................................................................... 105
5.3.
Сечения .................................................................................................... 107
5.4.
Заключение .............................................................................................. 116
Заключение .............................................................................................................. 117
Список литературы ................................................................................................. 119
3
Введение
Во многих отраслях науки и техники, в том числе нефтегазовой отрасли
актуальной
оказывается
неопределенности.
задача
Данная
оценки
задача
параметров
традиционно
в
условиях
решается
методами
статистической обработки информации. Статистические методы заменяют
исходные данные законом в виде уравнений зависимости параметров.
Статистика предполагает возможность повторного проведения эксперимента,
что невозможно при решении задач геологии, геофизики и нефтегазового дела.
Поэтому статистическая обработка информации нередко считается основой
ошибочных
заключений
в
условиях
неопределенности.
В
условиях
неопределенности данных и когда повторное проведение эксперимента
невозможно,
необходимо
использовать
технологию
анализа
оценки
достоверности. Данная технология учитывает неоднородную структуру данных
и тем самым, позволяет снижать технико-экономические риски. Такой
технологией является технология нечеткого моделирования, которая требует:
адаптации к конкретным задачам;
приспособления к конкретным проблемам;
развития математического аппарата.
Технология нечеткого моделирования позволила бы пересчитывать
неопределенность, присущую реальным данным и связям между реальными
данными, и определять достоверность пространственного распределения
одновременно измеренных значений параметров моделей.
Особенностями измеренных параметров в нефтегазовой отрасли,
лежащих
в
основе
создания
математических
моделей,
являются
неопределенность, фрагментарность и неоднородность. Данные обладают
стохастической природой, связанной с неоднородным внутренним строением
среды [76, 81, 82, 126]. Данная неопределенность включает также
неопределенность в связях между петрофизическими параметрами и в
результате наследуется в конечных построениях.
4
Нечеткие величины подчиняются своей алгебре – алгебре нечетких
величин. Задача анализа оценки достоверности должна быть перестроена под
эту алгебру, где роль умножения матриц играет нечеткий логический вывод,
где роль сложения играет логические операции объединения, а роль
умножения – логическая операция пересечения. Вся технология, таким
образом,
превращается
в
технологию
нечеткого
моделирования,
представляющую данные в виде нечетких величин, а зависимости – в виде
нечетких отношений. Адаптация, развитие технологии в приложениях,
решающих крупные отраслевые проблемы, является важнейшей задачей
математического моделирования.
За последние года вопросу изучения нечеткого моделирования при
решении задач в условиях неопределенности посвящены труды А. Е.
Алтунина, В. П Бочарникова, В. В. Борисова, И. П. Жабрева, М. В. Семухина,
Я. И. Хургина, Yi Ke Guo, Xu Dong Jing, Jose Finol, A. Bardossy, H. Paasche,
A. Ouenes, T. Tani, M. Sakoda, K. Tanaka. В данных трудах описаны общие
теоретические основы.
В работах А. И. Кобрунова, А. В. Григорьевых, А. С. Могутова,
А. Н. Дорогобед
описаны
результаты
адаптации
аппарата
нечеткого
моделирования, апробация на тестовых данных. Тем не менее, остался круг
вопросов, требующих развития новых методов решения этих задач.
Цель диссертационной работы: разработать теорию, численные
методы и математическую модель создания информационно-экспертных
оценок достоверности параметров в условиях неопределенности на основе
технологии нечеткого моделирования.
Основные задачи диссертационной работы:
1.
Разработать математические методы представления зависимостей
между реально измеренными данными в форме нечетких отношений,
выбирая для этой цели оптимальные базовые функции для аппроксимации
нечетких отношений.
5
2.
Проанализировать влияние значения эффективных параметров на
разброс данных при формировании нечеткого отношения на основе анализа
изменения сечений функции принадлежности.
3.
Разработать метод аналитического продолжения для получения
информации о распределении достоверности при удалении от заданной
точки.
4.
Разработать математический метод и комплекс программ,
позволяющие оценивать обеспеченность трехмерных моделей объективными
данными.
Объект
исследования:
связи
между
физико-геологическими
параметрами пористых сред.
Предмет исследования: неоднородный характер связей между
физико-геологическими параметрами, методы моделирования данных связей,
учет неопределенности исходных данных при прогнозировании промысловогеофизических параметров по геофизическим параметрам.
Методы исследований. В диссертационной работе используются
методы математического моделирования, методы регуляризации, решения
интегральных уравнений, нечеткого моделирования, нечеткого логического
вывода.
Научная новизна диссертационной работы:
предложено
анализировать
достоверность
экспертируемых
моделей на основе информационной модели, представляющей собой
нечеткое отношение между парой исходных и итоговых параметров, которые
строятся на основании композиции Мамдани промежуточных этапов,
определяющих граф связей между параметрами;
разработан критерий выбора эффективных параметров рассеяния
при аппроксимации нечетких отношений, позволяющий конструировать
нечеткие петрофизические модели адекватные исходным данным;
6
разработан математический метод интерполяции локальных
прогнозных функций принадлежности, определенных на функциональном
пространстве, позволяющий получать информацию о ранжированной
достоверности значений параметров;
разработан комплекс программ информационной экспертизы
результатов прогнозирования на основе неопределенных исходных данных.
Основные положения, выносимые на защиту:
Математическая
1.
принадлежности,
модель
основанная
на
представления
нормальном
законе
функции
распределения,
адекватна реальным данным и позволяет получать дифференцированную
оценку достоверности компонент модели, образованной из прогнозного
параметра.
2.
Математический
метод
нечеткого
логического
вывода,
оснащенный правилами представления исходных данных в форме нечетких
величин и связей в форме нечетких отношений, может служить основой для
экспертизы достоверности прогнозируемых параметров.
3.
Математический метод интерполяции локальных прогнозных
функций принадлежности, основанный на экспоненциальной модели,
позволяет выполнять пространственное распределение достоверности от
заданной точки по всей площади.
4.
Разработанный комплекс программ позволяет получать оценку
достоверности распределения параметров в построенной трехмерной модели
на основе анализа их соответствия имеющимся данным.
Теоретическая и практическая значимость.
Разработанные методы, алгоритмы и программы моделирования на
основе
принципов
нечеткого
моделирования
применяются
при
дифференцированной оценке подсчетов запасов углеводородов.
Разработанные методы, алгоритмы и программы моделирования
позволяют осуществлять проверку на достоверность компонент построенной
7
геологической модели на основе анализа их соответствия имеющимся
геолого-геофизическим данным.
Реализация результатов работы.
Результаты, полученные при выполнении работы, используются в
учебном процессе при подготовке специалистов и магистров по геологии
нефти и газа.
Степень достоверности и апробация работы.
Апробация разработанных методов выполнялась на фактических
материалах, относящихся к нескольким из месторождений ТиманоПечорской нефтегазоносной провинции. Результаты апробации показали
возможность использования разработанных методов с целью оценки на
достоверность построенных физико-геологических моделей имеющимся
исходным данным.
Основные
положения
диссертационного
исследования
были
представлены на международном семинаре им. Д. Г. Успенского (Москва,
2017;
Казань,
2018),
«Рассохинские
чтения»
(Ухта,
2015-2018),
международной молодежной научной конференции «СеверГеоЭкоТех»
(Ухта, 2015), III Школе-конференции «Гординские чтения» (Москва, 2015),
VI Международной студенческой электронной научной конференции
«Студенческий научный форум 2015», научных семинарах Ухтинского
государственного технического университета.
Публикации.
Основные положения диссертации напечатаны в 21 научной работе, в
том числе 5 в журналах из перечня ведущих периодических изданий,
рекомендованных ВАК.
Структура и объем работы.
Диссертационная работа включает введение, пять глав, заключение,
список литературы, состоящий из 151 наименования. Работа изложена на 134
листах, содержит 66 рисунков, 22 таблицы.
8
Глава 1.
Петрофизическое моделирование
Под петрофизическим моделированием можно понимать:
процесс построения объемного распределения в геологическом
пространстве различных физических параметров, характеризующих главные
структурно-вещественные комплексы изучаемого месторождения полезного
ископаемого либо другого геологического объекта;
анализ всей информации (петрофизических параметров породы),
проводимый с целью увязки и оценки ее непротиворечивости [1].
В
данной
работе
под
петрофизическим
моделированием
предполагается моделирование взаимосвязей между физико-геологическими
параметрами, которые имеют характер корреляционно-регрессионный,
аналитический и др. (например, в виде решающих правил экспертной
зависимости).
Петрофизическое моделирование необходимо в первую очередь при
подсчете запасов углеводородов для определения взаимосвязи промысловогеофизических параметров по геофизическим параметрам.
В первой главе на основе работ Добрынина В. М., Вендельштейна Б.
Ю., Кожевникова Д. А., Дахнова В. Н, Кобрунова А.И. и др. рассмотрены
связи
между
коллекторов,
физико-геологическими
проблемы,
параметрами
возникающие
при
в
разных
типах
моделировании,
типы
петрофизических моделей, их преимущества и недостатки.
1.1.
К
формированию
аналитических
зависимостей
различных петрофизических моделей
Коллектор – горная порода, имеющая пустоты, которая может иметь
флюиды (нефть, газ, воду) [103, 122, 128]. В основном выделяют три типа
коллекторов нефти и газа: терригенный, карбонатный и смешанный [113].
Терригенными коллекторами являются породы, состоящие из зерен
минералов. По литологии этот тип коллектора отличается размером зерен. У
9
алевролитов размер зерен варьируется в пределах от 0.05 до 0.1 мм, у
мелкозернистых песков – 0.1-0.25 мм, у крупнозернистых песков – 0.25-1 мм.
Терригенными коллекторами являются осадочные породы: песчаники, пески,
глины [121, 111]. Около 80% всех месторождений углеводородов России и
60% месторождений нефти в мире являются терригенными коллекторами
(почти все коллектора нефти и газа Западно-Сибирского бассейна [104, 110,
112, 114]).
В карбонатных коллекторах горные породы образуются за счет двух
основных минералов – доломита и кальцита. Поровое пространство образуют
трещины, каверны. К основным процессам, формирующим поровое
пространство
в
выщелачиванием
карбонатах,
и
относятся
карстообразованием
биогенные
или
с
накопления
с
тектоническими
напряжениями, которые привели к образованию сети каверн, трещин,
микротрещин и т.д. [127, 116, 118, 120]. Около 69% месторождений
углеводородов в мире являются коллекторами карбонатного типа (ВолгоУральская область [107] и Тимано-Печорская провинция, Прикаспийская
впадина).
Довольно часто встречаются смешанные коллектора – трещиноватые
породы, пустотное пространство которых образует как система трещин, так и
поровое пространство блоков, а также каверны и карст.
Примеры горных пород терригенного типа: пески, песчаники, глины (а,
б, в); и карбонатного типа: трещинные (в, г) и трещинно-кавернозные (д) [2]
представлены на рисунке 1.
10
Рисунок 1 – Примеры порового пространства в обломочных и глинистых
породах
а – хорошо отсортированный песчаник с высокой пористостью, б – плохо
отсортированный песчаник с низкой пористостью, в – глины, г –
трещиноватая порода (𝑘т ≈ 0.15%), д – трещиноватая порода (𝑘т ≈ 0.3%), е
– трещинно-кавернозная порода (𝑘т + 𝑘кав ≈ 5%)
Свойства коллекторов [109, 115, 119, 123, 131] включать (пористость) и
пропускать
(проницаемость)
флюиды
называются
фильтрационно-
емкостные. К таким свойствам коллекторов относятся:
пористость;
проницаемость;
водонасыщенность.
Пористость – свойство горных пород, отражающее объем пустот,
которые способны вмещать в себя газ и жидкости [108]. Как правило,
пористость измеряют двумя способами: при геофизических исследованиях
скважин (ГИС) [124, 130, 117] и в лабораторных условиях на образках керна.
Пористость бывает трех типов: общая, открытая и эффективная.
Общая пористость рассчитывается как сумма объема всех пор (𝑉п).
Коэффициент общей пористости (Кп) соответственно как отношение объема
всех пор на весь объем образца (𝑉):
Кп =
𝑉п
𝑉
11
Открытая пористость – сумма объема сообщающихся пор (𝑉по ).
Коэффициент
открытой
пористости
(Кпо )
–
отношение
объема
сообщающихся пор на весь объем образца:
𝑉по
Кп =
𝑉
о
Эффективная пористость – сумма объема пустот (𝑉пэ ), участвующих в
фильтрации, из которых может быть извлечены углеводороды. Коэффициент
эффективной пористости (Кпэ ) рассчитывается по формуле:
𝑉пэ
Кп =
𝑉
Для терригенных коллекторов с низким содержанием цемента
э
существенной разницы значения выше перечисленных типов пористости не
имеют. В коллекторах с высоким количеством изолированных пор, значения
пористости существенно различаются. Значения коэффициентов пористости
всегда соответствуют неравенству:
Кп > Кпо > Кпэ
Рассмотрим изменение открытой пористости при изменении глубины в
разных типах коллекторов. Петрофизические отношения коэффициента
открытой пористости и глубины для разных типов коллекторов приведены на
рисунках 2 и 3.
12
Рисунок 2 – Изменение значения открытой пористости с глубиной песчаноглинистых пород (скв СГ-1 Аралсорская)
а – песчаники, б – алевролиты, в – аргиллиты и глинисто-алевритовые
породы, г – сводная зависимость.
Рисунок 2 отображает динамику изменения коэффициента открытой
пористости при изменении глубины залегания песчано-глинистых пород.
Данные были получены при исследовании кернового материала Аларсоркой
скважины, расположенной в Прикаспийской низменности. Уменьшение
значения открытой пористости при изменении глубины объясняется
увеличением давления при проседании осадочных пород [2, 125].
13
Рисунок 3 – Динамика изменения значения открытой пористости с глубиной
залегания известняков
а – для «скелетных», б – для «бесскелетных» известняков
Рисунок 3 отображает динамику изменения значения коэффициента
открытой пористости известняков при изменении глубины залегания пород.
Данные были получены при исследовании известняков Оренбургского свода
«скелетных» и «бесскелетных». «Скелетные» известняки характеризуются
преимуществом форменных элементов над цементирующей основой, а
«бесскелетные» – базальной перекристаллизованной карбонатной массой. В
первом случае в осадках твердые частицы создают устойчивый каркас,
который не поддается уплотнению, поэтому изменение значения пористости
при изменении глубины не существенно. Во втором случае значение
пористости при изменении глубины уменьшается, так как жесткого скелета
нет, из-за чего происходит уплотнение. На глубине ниже 4 км в данных видах
осадочных пород наблюдаются лишь единичные поры [3].
Значение
пористости
в
терригенных
породах
больше,
чем
в
карбонатных (редко превышает 15 %).
Проницаемость – свойство горных пород пропускать через себя
флюиды при наличии давления, иначе, площадь сечения пористой среды,
14
через которую идет фильтрация. Высокую проницаемость, как правило,
имеют: пески, песчаники, доломиты, алевролиты; низкую – глины,
известняки.
Коэффициент проницаемости породы (Кпр, численная характеристика)
определяется по линейному закону фильтрации Дарси: скорость фильтрации
флюидов пропорциональна градиенту давления и обратно пропорциональна
динамической вязкости. Таким образом, коэффициент проницаемости
рассчитывается по формуле:
Кпр =
𝑄𝜇𝐿
,
𝑆∆𝑝
где 𝑄 – расход жидкости или газа и их смесей с вязкостью 𝜇,
протекающих через объем породы 𝑆 под действием градиента давления
Проницаемость
бывает
трех
типов:
абсолютная,
∆𝑝
𝐿
фазовая
[4]
и
относительная.
Абсолютная
проницаемость
–
проницаемость
пористой
среды,
определенной одной фазой, химически и физически инертной к среде.
Фазовая проницаемость – проницаемость пористой породы для
определенной жидкости или газа при наличии еще одной или нескольких
фаз.
Относительная проницаемость – отношение фазовой проницаемости к
абсолютной проницаемости [129].
Рассмотрим изменение значения коэффициента проницаемости при
изменении пористости для разных типов коллектора. На рисунках 4 и 5
представлены петрофизические отношения коэффициента проницаемости и
эффективной пористости для терригенных и карбонатных типов коллекторов
соответственно.
15
Рисунок 4 – Петрофизическое отношение коэффициента проницаемости и
эффективной пористости для группы сложных терригенных коллекторов
одного из месторождений Западной Сибири [2]
На рисунке 4 видно, что проницаемость возрастает в более пористой
среде. В терригенных коллекторах коэффициент проницаемости значим даже
при малых значениях эффективной пористости.
Рисунок 5 – Зависимость коэффициента проницаемости от открытой (а) и
эффективной пористости (б) Беркутовской газоносной зоны Мраковской
депрессии Предуральского краевого прогиба Волго-Уральской
нефтегазоносной провинции
Данные: 1 – 1980 г., 2 – 1984 г.
16
Из рисунка 5 видно, также как и в терригенных коллекторах,
коэффициент проницаемости увеличивается при увеличении коэффициента
пористости. По данной зависимости видно, что в карбонатных коллекторах к
проницаемым породам можно отнести породы с пористостью, превышающей
2,5 % [5].
Водонасыщенность – свойство породы, отображающее: сколько
остаточной
воды
содержится
в
порах.
Коэффициент
остаточной
водонасыщенности рассчитывается по формуле:
Ков =
𝑉ов
,
𝑉п
где 𝑉ов – объем остаточной воды, 𝑉п – объем пор
Предполагается,
что
до
появления
углеводородов,
поровые
пространства были заполнены водой. При появлении углеводородов, вода
уходит вниз, вытесняемая нефтью и газами. Для подсчета углеводородов,
необходимо рассчитать остаточную водонасыщенность [6].
Рассмотрим изменение остаточной водонасыщенности при изменении
пористости для разных типов коллектора. На рисунках 6 и 7 представлены
зависимости коэффициента остаточной водонасыщенности от коэффициента
пористости
для
терригенных
и
карбонатных
типов
коллекторов
соответственно.
17
Рисунок 6 – Отношение между параметрами «остаточная
водонасыщенность» и «пористость» для верейских песчаников из скважины
555 Покровского месторождения [7]
На
рисунке
6
отображена
динамика
изменения
остаточной
водонасыщенности от значения пористости в терригенных породах. Чем
больше пористость пород, тем меньше остаточной воды в порах.
Рисунок 7 – Отношение между остаточной водонасыщенностью и пористостью
Беркутовской газоносной зоны Мраковской депрессии Предуральского
краевого прогиба Волго-Уральской нефтегазоносной провинции
1 – плотные непроницаемые породы, 2 – трещинные коллекторы, 3 – поровотрещинные и поровые коллекторы, 4 – усредненная кривая для поровотрещинных и поровых коллекторов, 5 – расчетная кривая для вторичных пор [5]
18
Характер
отношения
между
параметрами
«остаточная
водонасыщенность» и «пористость» для карбонатных пород (см. Рисунок 7)
такой же, как и для терригенных пород.
При подсчете запасов нефти и газа с использованием данных
геофизических исследований главным вопросом является увязка данных ГИС
с данными, полученными в лабораторных условиях (на образах керна). При
определении пористости наблюдаются отклонения в данных, полученных
разными методами. Ниже приведены данные зависимости для карбонатных
пород (см. Рисунок 8) и терригенных пород (см. Рисунок 9).
Рисунок 8 – Отношение между параметрами пористость по ГИС и
пористость по керну карбонатных пород Салюткинского месторождения
Тимано-Печорской нефтегазовой провинции
19
Рисунок 9 – Отношение между параметрами пористость по ГИС и
пористость по керну терригенных отложений Арыстановского
Месторождения [8]
В нефтеносных коллекторах нефтенасыщенность пород (Кн), как
правило, определяется через коэффициент остаточной водонасыщенности по
формуле:
Кн = 1 − Ков
Аналогично
рассчитывается
коэффициент
газонасыщения
(Кг)
газоностных коллекторов:
Кг = 1 − Ков
Выше
описанные
отношения
между
физико-геологическими
параметрами, независимо от типа коллектора, имеют неоднородный,
неопределенный характер, обусловленный неоднородностью строения среды,
что приводит к ряду существенных проблем.
20
Проблемы,
1.2.
возникающие
при
моделировании
петрофизических сред
Петрофизическая модель адекватна описанию изучаемого объекта,
если модель и объект должны иметь одинаковую размерность физических
величин. Это одно из главных требований, которое должно всегда
выполняться, и может служить критерием проверки на адекватность модели.
При
анализе
различных
размерности
параметров
обосновываются
моделируемой
гипотезы
системы.
о
взаимосвязи
Зачастую,
анализ
размерностей используется для определения формул, за счет составления
выражения из параметров, определяющих изучаемый объект, имеющего
нужную размерность. Анализ размерностей напрямую связан с понятием о
подобии.
Петрофизическое
моделирование
должно
быть
основано
на
математической теории подобия. В теории подобия предполагается, что
модель должна достаточно хорошо отображать процессы моделируемого
объекта. Теория подобия позволяет определить самые информативные
параметры,
определяющие
процессы
в
системе,
масштабы
данных
параметров, при которых характеристики объекта и моделируемой системы
подобны. Использование теории подобия в петрофизическом моделировании
позволяет определить взаимосвязи между петрофизическими параметрами и
объединить переменные величины в безразмерные параметры, которые в
большей степени связаны с физической природой моделируемого объекта
[9].
Принцип подобия заключается в том, что в изучаемом объекте и
модели процессы развиваются аналогично, если модель и объект подобны.
Модель и изучаемый объект подобны, если и их параметры подобны.
Подобие параметров определяется коэффициентом подобия.
Коэффициентом
подобия
называется
отношение
определенного
параметра объекта к такому же параметру в модели. Если коэффициенты
подобия всех параметров равны, то модель и объект подобны.
21
К видам подобия относятся:
геометрическое подобие;
кинематическое подобие;
динамическое подобие.
Геометрическим подобием обладают объекты и их модели, в которых
отношения всех линейных характеристик (длина, площадь, объем и т.д.)
равны.
Кинематическим подобием обладают объекты и их модели, имеющие
геометрическое
подобие,
и
сохраняющие
отношения
скоростей
соответствующих точек.
Динамическим подобием обладают объекты и их модели, обладающие
кинематическим подобием, и сохраняющие отношения действующих сил на
соответствующие точки [10].
Классификация петрофизических моделей
1.3.
Петрофизические модели можно разделить на:
а.
физико–аналитические;
b.
корреляционно–регрессионные;
c.
экспертные правила;
d.
нечеткие модели.
В качестве физико-аналитической петрофизической модели чаще всего
выступает интегро-дифференциальное уравнение физических параметров
коллектора
[40,
106].
Данное
уравнение
обеспечивает
нахождение
промыслово-геофизических свойств коллектора по измеренным параметрам
ГИС и включает большое количество параметров, которые невозможно
вычислить по данным экспериментов. Эти параметры подбираются из
экспериментальных данных. Таим образом зависимости, которые находятся в
основе физико-аналитической модели, являются параметрически заданным
выражением для уравнения регрессии. Подобранные параметры достаточно
22
часто не имеют отношения к физическому смыслу, который имеет параметр
регрессии. Основанием этого является возможная компенсация невязок
экспериментальных
данных
и
теоретических
моделей,
описываемых
комплексом параметров. Невязка каждого параметра определена алгоритмом
минимизации, который лежит в основе метода регрессии. Если больше двух
параметров, появляется взаимное влияние и возможен перенос невязок
одного параметра на другой. Разделить данную невязку чаще всего не
удается.
Корреляционно-регрессионными моделями являются регрессионные
зависимости, которые сопровождаются оценкой тесноты связи [51, 80]. Они
конструируются на основе имеющихся выборок экспериментальных данных.
При конструировании регрессионных уравнений принимается их
представление в виде аналитической модели с параметрами, которые
необходимо определить. Данная
модель строится на
базе физико-
аналитической петрофизической модели либо визуального анализа данных.
Сконструированная модель всегда сопровождается оценкой тесноты и
достоверности
связи
между
исходными
данными
и
рассчитанной
регрессионной зависимостью.
Прогнозирование
промыслово-геофизических
параметров
по
измеренным геофизическим параметрам осуществляется за счет построенных
регрессионных зависимостей. Промыслово-геофизические параметры могут
быть
использованы
нефтегазоностности.
для
прогнозирования
Погрешности
в
параметров
промежуточных
перспектив
зависимостях
передаются к итоговому прогнозному параметру. Для снижения техникоэкономических рисков, относящихся к принятию ошибочных решений по
разработке месторождения, необходимо рассчитывать неопределенность
итогового прогноза. Такую возможность корреляционно-регрессионные
модели не дают.
Петрофизическими моделями на основе экспертных правил являются
систематизированные заключения экспертов в виде решающих правил, о
23
физико-геологических параметрах на базе количественных оценок и
качественных сравнений измерений неоднородных данных геофизики и
гидродинамики. Экспертные заключения распространяются на тот объект
или группе объектов, на которые распространяется квалификация экспертов.
Экспертные правила разных экспертов, вне зависимости от квалификации,
могут
иметь
противоречивый
характер,
так
как
эксперты
могут
придерживаться разных геологических гипотез. На основе экспертных
правил формируется заключение о характере флюидонасыщения.
Нечеткие петрофизические модели описывают экспериментальные
данные в форме функций принадлежности для нечеткого отношения между
физико-геологическими
параметрами.
Данный
вид
моделей
является
нетрадиционным и для их построения необходимо использовать большой
объем
данных.
Но,
несмотря
на
данные
недостатки,
нечеткие
петрофизические модели имеют высокую степень наглядности, а также
позволяют получать для каждого прогнозного параметра результат,
дифференцированный по достоверности. Таким образом, данные модели
позволяют выполнить проверку на обеспеченность исходной информацией,
получить оценку достоверности построенных геологических моделей. Это
дает возможность многовариантных построений геологических моделей,
которые отображают реальную информацию, а также управления рисками,
относящимися к принятию решений по разработке месторождения.
1.4.
Корреляционно-регрессионные
петрофизические
модели
Регрессионный
анализ
является
главным
инструментом
при
построении корреляционно-регрессионных петрофизических моделей для
установления связей между физико-геологическими параметрами. Общий
вид данных моделей – уравнения регрессии, в которых необходимо
24
подобрать параметры, сформированные на базе физико-аналитических
моделей и анализа рассеяния исходных данных.
Распространенными видами уравнений регрессии являются:
линейные зависимости;
экспоненциальные;
гиперболические;
логарифмические.
Параметры
в
зависимостях
подбираются
при
помощи
метода
наименьших квадратов. В случае линейной зависимости после подбора
параметров оценивается теснота связи полученного уравнения и исходных
данных путем вычисления коэффициента корреляции (𝑟):
𝑛
𝑟𝑥𝑦
1
=
∑(𝑥𝑖 − 𝑥̅ )(𝑦𝑖 − 𝑦̅)
𝑛𝜎𝑥 𝜎𝑦
𝑖=1
где 𝑛 – количество наблюдений; 𝑥̅ , 𝑦̅ – выборочные средние; 𝜎𝑥 , 𝜎𝑦 –
это среднеквадратическая величина отклонения.
Достоверность оценки регрессионных моделей в случае линейной
регрессии выполняется за счет вычисления коэффициента детерминации
(𝑅2 ):
∑𝑛𝑖=1(𝑦𝑖 − 𝑦)2
𝑅 = 𝑛
∑𝑖=1(𝑦𝑖 − 𝑦̅)2
2
где 𝑦 – предсказанные (вычисленные на основе уравнения регрессии)
данные; 𝑦𝑖 – исходные экспериментальные значения данных.
Вычисления коэффициента детерминации в случае нелинейной
регрессии выполняется по формуле:
∑𝑛𝑖=1(𝑦𝑖 − 𝑦)2
𝑅 =1− 𝑛
∑𝑖=1(𝑦𝑖 − 𝑦̅)2
2
Чем ближе коэффициент детерминации к единице, тем лучше
регрессионная зависимость описывает исходные данные.
25
Построение
корреляционно-регрессионных
моделей
предполагает
распределения ошибок в данных в форме нормального закона, которому не
подчиняются распределения измеренных значений, возникающие в реальных
данных. Построение таких моделей по реальным экспериментальным
данным, отображающим физические свойства горных пород, может быть
причиной ошибочных результатов оценивания параметров.
Приведем некоторые примеры расчета зависимостей и рассмотрим
результаты, выявляющие их недостатки.
На рисунке 10 представлена зависимость акустического импеданса от
пористости образцов. Со степенью достоверности 𝑅2 = 0.96 можно
использовать приведенную
на рисунке
линейную
зависимость:
𝑦=
−0.2437𝑥 + 11.696, где 𝑥 – пористость, 𝑦 – акустический импеданс. Данная
зависимость также характеризуется высоким коэффициентом корреляции
(𝑟 = 0.98),
что
при
высокой
степени
достоверности,
позволяет
рекомендовать использование данной зависимости. Тем не менее, данные
неоднородны: при низкой пористости данных мало, с увеличением
пористости данных становится больше (см. Рисунок 11). Достоверность
прогноза должна быть различной на разных участках.
Рисунок 10 – Зависимость акустического импеданса от пористости образцов
26
9
Количество данных
8
7
6
5
4
3
2
1
0
Пористость, %
Рисунок 11 – Количество измеренных значений отношения «Пористость –
акустический импеданс» на разных диапазонах пористости
На рисунке 12 представлена зависимость водонасыщенности от
пористости образцов. В данном случае регрессионная зависимость имеет
логарифмический характер: lg(𝑦 − 6) = 2.75 − 0.095𝑥, где 𝑦 – коэффициент
водонасыщенности,
𝑥
–
коэффициент
пористости.
Со
степенью
достоверности 𝑅 2 = 0.92 можно использовать данную зависимость. Также
как и в предыдущем примере, данные неоднородны: при низкой пористости
данных мало, с увеличением пористости данных становится больше, на
диапазоне от 12 до 17 наблюдается количество данных (см. Рисунок 13).
Достоверность прогноза должна быть различной на разных участках.
Рисунок 12 – Зависимость остаточной водонасыщенности от пористости
27
3,5
Количество данных
3
2,5
2
1,5
1
0,5
0
Пористость, %
Рисунок 13 – Количество измеренных значений отношения «Пористость –
водонасыщенность» на разных диапазонах пористости
Рассчитав регрессионные зависимости для отношений «пористость –
акустический импеданс» и «пористость - водонасыщенность», подставляем
одно уравнение в другое и получаем зависимость между параметрами
«акустический импеданс» и «водонасыщенность» (см. Рисунок 14): 𝑧 =
10−1.86+0.39𝑥 + 6, где 𝑥 – акустический импеданс, 𝑧 – остаточная
водонасыщенность.
90
Водонасыщенность, %
80
70
60
50
40
30
20
10
0
5,5
6,5
7,5
8,5
Акустический импеданс,
9,5
10,5
106кг/м2с
Рисунок 14 – Зависимость остаточной водонасыщенности от акустического
импеданса
28
Корреляционно-регрессионная
модель
зависимости
остаточной
водонасыщенности от акустического импеданса характеризуется достаточно
низкой достоверностью, так как исходные данные характеризуются
различной плотностью на разных диапазонах переменных. Это является
причиной
различной
меры
доверия
к
компонентам корреляционно-
регрессионной модели, а также изменчивости достоверности прогноза на
разных диапазонах значения параметров.
Таким образом, необходимо развитие аппарата, дифференцирующего
результаты прогноза параметров по достоверности на основе обеспеченности
данными
[150,
151].
нефтегазоностности
При
по
прогнозировании
геофизическим
параметров
параметрам
перспектив
через
цепочку
подстановок плотность данных уменьшается. Данное обстоятельство не
учитывается в процессе построения корреляционно-регрессионных моделей,
при которых все данные заменяются уравнениями без учета плотности на
разных участках. Не однородность в данных объясняется особенностями
месторождения:
неоднородностью
залежей,
разрывами
и
разломами.
Достоверной картины прогноза не может быть получено, но следует
оперировать
данными,
которые
были
получены.
В
данном
случае
необходимо получать прогноз с оценкой достоверности на разных участках,
что возможно при построении нечетких петрофизических моделей.
1.5.
Заключение
Традиционно
статистической
задача
прогноза
обработки
параметров
информации,
решается
которые
методами
заменяют
экспериментальный материал рассчитанным законом в форме уравнений
отношения между параметрами, которые переносятся на исследуемый
объект. Данный путь достаточно часто является основой ложных, неточных
геологических решений о запасах и ресурсах углеводородного сырья (итог
сравнительного анализа возможностей вероятностно-статистических методов
29
(Монте-Карло [53]) и теории нечетких множеств с целью оценки
неопределенностей в процессе подсчета запасов нефти и газа объемным
методом представлен в работе [11]). По этой причине предлагается
исследовать данные в полном объеме, не заменяя уравнениями, а
представлять их в форме нечетких множеств.
30
Глава 2.
Фазификация нечетких отношений, основанная на
алгоритме сжатия информации
При одновременном измерении двух или более параметров на одном
образце появляется «ложная» информация, включающая также ошибки при
измерении. «Ложная» информация не редко присутствует в экспериментах.
От таких ошибок, присутствующих в наборе измеренных значений,
необходимо
избавляться.
Оставшаяся
необнаруженной
«ложная»
информация искажает результаты измерений. Для выхода из данной
ситуации предлагается рассматривать исходные данные взаимосвязи между
параметрами
как
нечеткие
отношения
и
конструировать
функции
принадлежности с помощью аппарата, основанного на алгоритме сжатия
информации. Алгоритм сжатия подразумевает, что данные об одновременно
измеренных значениях параметров исходят от источников, которые в свою
очередь позволяют существенно снизить размерность измеряемых величин.
2.1.
Нечеткое моделирование
Началом развития нечеткого моделирования можно считать 1965 год, в
котором вышла статья профессора Калифорнийского университета Лотфи
Заде «Fuzzy Sets» (Нечеткие множества). В ней были введены базовые
определения и операции над нечеткими множествами [12]. Далее теорию
аппарата нечетких множеств развивали ученые Лофти Заде [67-71], Э.
Мамдани [19, 66], Д. Дюбуа [62], Р. Беллман [60, 61], А. Кофман [27], А. В.
Алексеев [58, 59], А. В. Леоненков [64], А. Н. Кирсанов [63], А. И. Орлов
[65].
Методы, основанные на нечеткой алгебре, нечеткой логике и теории
нечеткого
логического
вывода,
становятся
распространенными
в
нефтегазопромысловом деле как среди российских ученых в работах А. Е.
Алтунина [11, 30, 31, 84-90], В. П Бочарникова [56], В. В. Борисова [75], И.
П. Жабрева [54], А. И. Кобрунова, М. В. Семухина [11, 30, 31, 72, 73], Я. И.
Хургина [55], так и зарубежных Yi Ke Guo, Xu Dong Jing, Jose Finol [57], A.
31
Bardossy [28], H. Paasche [29], Ouenes A. [38], Tani T., Sakoda M., Tanaka K.
[39]. Вероятная оценка запасов углеводородов отражена в трудах [32-34].
Как правило, нечеткое моделирование использовалось исключительно
для решения вспомогательных задач и не использовалось для определения
подсчетных параметров.
В работах А. И. Кобрунова, А. В. Григорьевых, А. С. Могутова, А. Н.
Дорогобед было предложено рассматривать отношения между измеренными
параметрами
в
геофизических
исследованиях
в
качестве
нечётких
отношений, устанавливаемых между нечёткими переменными. Построение
функций
принадлежности
экспоненциальной
модели,
предлагалось
а
выполнять
исключение
общих
на
основе
переменных
при
конструировании цепочки уравнений осуществлять с помощью нечеткого
логического вывода (композиции Мамдани). [13, 14, 15, 16, 17, 41, 42, 48].
Недостатками данных работ являлось то, что: при конструировании функции
принадлежности рассмотрения ввелись по интервалам на оси X, из-за чего на
каждом интервале выделялся максимум, и не учитывались все данные
отношения одновременно; не использовался алгоритм для выявления
«ложной»
информации;
не
были
подробно
рассмотрены
свойства
композиции Мамдани.
Рассмотрим
основные
понятия
и
определения
нечеткого
моделирования, а также операции над нечеткими множествами.
Нечеткая переменная может быть описана параметрами: именем – 𝑎;
универсумом, содержащим все значения 𝑎, – 𝑋; нечетким подмножеством
универсума 𝑋, для каждого элемента которого определена функция
принадлежности 𝜇𝐴 (𝑥) – 𝐴.
Нечеткое
отношение
𝑅
–
нечеткое
подмножество
декартова
произведения универсальных множеств 𝑈1 , 𝑈2 , … 𝑈𝑛 .
32
Нечеткое множество 𝑀 – множество, каждый элемент которого
принадлежит этому множеству с определенной достоверностью. Можно
записать как 𝑀 = {(𝑥, 𝜇𝑀 (𝑥))|𝑥𝜖𝑋}.
Функция принадлежности 𝜇𝐴 (𝑥) определяет степень принадлежности
каждого элемента множества [67, 62, 77, 78, 79].
Основные операции над нечеткими множествами:
1.
(см.
Дополнением (логическое отрицание) нечеткого множества 𝔄
Рисунок
15)
является
нечеткое
множество
̅
𝔄
с
функцией
принадлежности:
𝜇𝔄̅ (𝑥) = 1 − 𝜇𝔄 (𝑥)
Рисунок 15 – Операция дополнения нечеткого множества
2.
Объединением (логическое «или») нечетких множеств 𝔄 и 𝔐 (см.
Рисунок 16) является нечеткое множество с функцией принадлежности:
𝜇𝔄∪𝔐 (𝑥) = max{𝜇𝔄 (𝑥), 𝜇𝔐 (𝑥)}
33
Рисунок 16 – Операция объединения нечетких множеств
3.
Пересечением (логическое «и») нечетких множеств 𝔄 и 𝔐 (см.
Рисунок 17) является нечеткое множество с функцией принадлежности:
𝜇𝔄∩𝔐 (𝑥) = min{𝜇𝔄 (𝑥), 𝜇𝔐 (𝑥)}
Рисунок 17 – Операция пересечения нечетких множеств
В
методе
рассматриваются
нечеткого
как
моделирования
нечёткие
экспериментальные
величины,
используемые
данные
для
конструирования функций принадлежности, содержащей меру доверия
каждому возможному значению.
Задачей фазификации служит представление всех используемых
данных и зависимостей в форме нечетких величин и нечетких отношений.
Поскольку и нечеткие величины, и нечеткие отношения характеризуются
своими
функциями
принадлежности,
задачей
фазификации
служит
34
нахождение конкретных функций принадлежности нечетких величин
адекватных исходным данным.
Фазификация состоит в представлении данных в форме нечетких
величин. В качестве данных выступают результаты измерений, как значений
одного переменного, так и одновременного измерения системы переменных
относящихся к одному объекту. В последнем случае, данные относятся к
созданию нечеткой петрофизической модели. Задача фазификации состоит в
адекватном подборе функции принадлежности, характеризующей нечеткость
результатов измерения переменной, или нечеткость регистрируемой связи
между параметрами, определяемой как нечеткого отношения между
переменными [132, 135, 136].
Формирование функций принадлежности
2.2.
Пусть
{𝑠𝑖 , 𝑖 = 1 ÷ 𝑁}
в
N-мерном
фазовом
экспериментально
пространстве
измерены
𝑆
значения
параметров
𝐬=
𝐬𝑗 𝜖𝑆, 𝑗 = 1 ÷ 𝑀,
образующие в нем подмножество 𝔄 ⊂ 𝑆. Данные 𝔄 ⊂ 𝑆 являются основанием
для формирования экспертного правила прогноза одних параметров 𝛈 из
числа 𝑠𝑖 , 𝑖 = 1 ÷ 𝑁, по измеренным другим, обозначаемым ξ: 𝐬 = ξ × 𝛈.
Традиционно в качестве такого правила в задачах промысловой геофизики
служат петрофизические зависимости [18], представляющие данные 𝔄 в
форме регрессионных зависимостей одних параметров от других. Такой
способ достаточно конструктивен и широко применяется в настоящее время,
но он не учитывает реальные свойства неоднородностей в изучаемых
параметрах.
̅ из
Экспертное правило 𝔄[ξ̅] заключения о прогнозируемых параметрах 𝛈
информации об измеренных ξ̅, согласованное с обучающими данными 𝔄 ⊂ 𝑆,
реализуется: определением по обучающей выборке 𝔄 нечеткого отношения
𝜇𝔄 (ξ, 𝛈) = 𝜇𝔄 (𝐬) [12] между параметрами 𝛈 и ξ; представлением конкретных
35
измерений ξ̅ в форме нечеткой величины и последующим использованием
алгоритма нечеткого логического вывода [19] о нечеткой величине 𝛈.
Принципом представления измеренных значений параметра ξ̅, по
которым выполняется прогноз, в форме нечетких величин, а взаимосвязей
между параметрами 𝔄 в форме нечетких отношений является понятие поле
рассеяния.
Полем рассеяния для данных 𝐬𝑗 𝜖𝑆, 𝑗 = 1 ÷ 𝑀 назовем функцию 𝔄𝜀 (𝐬),
такую, что для каждой подобласти ∆𝑆 из разбиения 𝑆:
𝜀
max |𝔄 (𝐬)∆𝑆 − 𝔄(∆𝑆)| ≤ 𝜀,
(1)
∆𝑆𝜖𝑆
где 𝔄(∆𝑆) – число значений из 𝔄 ⊂ 𝑆, целиком лежащее в ∆𝑆.
Величина 𝔄(∆𝑆) называется картой плотности данных.
Функция
принадлежности
𝜇𝔄 (𝐬)
для
измеренных
значений
параметров 𝐬𝜖𝑆 как нечетких величин есть нормированное к единице поле
рассеяния 𝔄𝜀 (𝐬):
𝔄𝜀 (𝐬)
𝜇𝔄 (𝐬) =
.
max[𝔄𝜀 (𝐬)]
(2)
𝐬
Данные, по которым проводится большая часть экспериментов,
представленных в данной работе, приведены на рисунках 18 и 19. Первым
отношением
была
параметрами
пористости
интерпретации
выбрана
материалов
по
взаимосвязь
ГИС
между
(определяется
геофизических
петрофизическими
по
исследований
результатам
скважин)
и
пористости по керну (определяется по результатам анализа кернового
материала) карбонатных пород Салюткинского месторождения ТиманоПечорской нефтегазовой провинции (данные, взяты из [8]). Второе
отношение – взаимосвязь между параметрами «пористость по керну» и
«остаточная водонасыщенность пород» (определяются по результатам
анализа кернового материала) без учета класса коллектора Беркутовской
газоносной зоны Мраковской депрессии Предуральского краевого прогиба
Волго-Уральской нефтегазоносной провинции (данные взяты из [5]).
36
Отношения
между
данными
параметрами
являются
основными
и
Пористость по керну, %
необходимыми, характеризующие среду на нефтегазоперспективность.
18
16
14
12
10
8
6
4
2
0
0
2
4
6
8
10
12
14
16
Пористость по ГИС, %
а
б
Рисунок 18 – Исходные данные отношения «пористость по ГИС – пористость
Остаточная
водонасыщенность, %
по керну» (а) и соответствующая карта плотности данных (б)
100
90
80
70
60
50
40
30
20
10
0
0
2
4
6
8 10 12 14 16 18 20
Пористость по керну, %
а
б
Рисунок 19 – Исходные данные отношения «пористость по керну –
остаточная водонасыщенность» (а) и соответствующая карта плотности
данных (б)
2.2.1 Алгоритм сжатия информации
Для построения поля рассеяния 𝔄𝜀 (𝐬) разобьём 𝑆 сеткой ∆𝑆𝑗 , 𝑗 = 1 ÷ 𝑀
так, чтобы была покрыта вся область значений параметров из 𝔄,
∑𝑗=1÷𝑀|𝔄(∆𝑆𝑗 )| = 𝑁 и поставим задачу поиска функции рассеяния по (1):
37
𝜀
max |𝔄 (𝐬)∆𝑆𝑗 − 𝔄 (∆𝑆𝑗 )| ≤ 𝜀
(3)
∆𝑆𝜖𝑆
Выберем базисную систему функций 𝐾(𝐬 𝑘 , 𝐬) и поставим задачу
нахождения 𝔄𝜀 (𝐬) в виде:
𝐾
𝔄
𝜀 (𝐬)
= ∑ 𝜑(𝐬 𝑘 )𝐾(𝐬 𝑘 , 𝐬),
(4)
𝑘=1
где 𝜑(𝐬 𝑘 ), требующая нахождения из (4), функция.
Соотношение (4) означает, что поле рассеяния рассчитывается в форме
линейной комбинации базовых функций 𝐾(𝐬 𝑘 , 𝐬), свойства которых
определяют принцип аппроксимации поля рассеяния 𝔄𝜀 (𝐬) и функции
принадлежности 𝜇𝔄 (𝐬). Физический смысл задачи (4) заключается в
нахождении 𝜑(𝐬 𝑘 ) по заданной 𝔄𝜀 (𝑠) и основан на представлении о том, что
рассчитанное
значение
функции
принадлежности,
следующее
из
экспериментальных данных, это результат диффузии, который привел к
сглаживанию более точной зависимости 𝜑(𝐬 𝑘 ). Величина этого сглаживания
контролируется эффективным параметром рассеяния ζ. Задача (4) это
нечеткое уравнение Фредгольма первого рода с ядром:
|𝐬 − ξ|2
𝐾(𝐬 − ξ, 𝜁) =
exp (−
),
𝜁2
√𝜋𝜁
1
(5)
которая относится к классу некорректных, и требует для своего решения
приближенных устойчивых процедур [20, 21, 49]. В качестве такой
процедуры
примем
основанную
на
оптимизационном
принципе
компактификации источников.
В уравнении (4) подынтегральным выражением являются источники
𝜑(𝐬 𝑘 ), а левая часть – отклик на них. Принимаем в качестве отклика
найденные значения 𝜇𝔄𝜀 (𝐬) =
𝔄𝜀 (𝐬)
max𝐬 [𝔄𝜀 (𝐬)]
по введенной сетке ⋃𝑖 ∆𝑆𝑖 =
𝑆; ⋂𝑖≠𝑗 ∆𝑆𝑖 = ∅. Функция 𝜇𝔄𝜀 (𝐬) есть значения в ячейках ∆𝑆𝑖 нормированной
функции рассеяния 𝔄𝜀 (𝐬). Таким образом, левая часть в (4), это набор
дискретных значений по ячейкам ∆𝑆, так же, как и искомая 𝜑(𝐬 𝑘 ). Принцип
38
компактификации
источников
формулируется
как
выбор
такого
распределения источников 𝔄min (𝐬 𝑘 ), 𝐬 𝑘 ⊂ ∆𝑆𝑖 , и ей соответствующей
функции принадлежности с 𝜇𝔄min (𝐬), которые характеризуются минимальным
числом 𝐾 носителей (ненулевых значений 𝔄min (𝐬 𝑘 ) в покрытии ∆𝑆𝑖 ),
выбранной сеткой фазового пространства 𝑆, при условии приближенного
выполнения соотношения (4). Это соответствует задаче целочисленного
программирования:
‖
1
𝐾
∑ [𝜇𝔄min (𝐬𝒌 ) ∗ exp (−
√𝜋𝜁 𝑘=1
𝐾 → min.
|2
2
|𝐬 − 𝐬𝑘
)] − 𝜇𝔄𝜀 (𝐬)‖ ≤ 𝜀,
2
𝜁
(6)
Решение данной задачи выполняется численно, перебором значений 𝐾
от максимального, определенного построенными значениями 𝔄𝜀 (𝐬), в
ячейках сетки ⋃𝑖 ∆𝑆𝑖 = 𝑆; ⋂𝑖≠𝑗 ∆𝑆𝑖 = ∅ до минимального значения, при
котором
условие
неравенства
еще
выполняется.
Результат
решения
𝜇𝔄min (𝐬𝒌 ), 𝑘 = 1 ÷ 𝐾 называется картой плотности источников данных или
полем источников данных. Алгоритм построения карты источников
называется алгоритмом сжатия информации. Исходя из карты плотности,
источников поля рассеяния, могут быть получены параметризованные
эффективным параметром рассеяния ζ модели функции принадлежности для
характеристики нечеткого отношения.
Для приведенных на рисунках 18 и 19 исходных данных, и им
соответствующих карт плотности данных (поля рассеяния), были построены
карты плотности источников поля рассеяния, которые приведены на
рисунках 20 и 21 соответственно.
39
а
б
Рисунок 20 – Карта плотности данных для нечеткого отношения «пористость
по ГИС – пористость по керну» (а) и соответствующая карта источников
данных(б)
а
б
Рисунок 21 – Карта плотности данных для нечеткого отношения «пористость
по керну – остаточная водонасыщенность» (а) и соответствующая карта
источников данных (б)
Переход к плотности источников поля рассеяния при аппроксимации
функции принадлежности позволяет существенно снизить размерность
величин.
40
2.2.2 Выбор аналитического вида функции принадлежности
В настоящей работе рассматривается три модели, лежащие в основе
конструирования функции принадлежности: модель конуса (является
традиционной [24]), экспоненциальная модель и модель обратных квадратов
[133, 141].
Рассмотрим их более подробно:
Модель конуса – традиционный прием при аппроксимации функции
принадлежности. Информация в данном случае равномерно расходится от
источников. Базисная система функций задана уравнением:
0 𝐬 ≤ 𝑎, 𝑐 ≤ 𝐬
𝐬−𝑎
𝑎<𝑠≤𝑏
𝐾(𝐬, a, b, c) = 𝑏 − 𝑎
𝑐−𝐬
𝑏<𝐬<𝑐
{𝑐 − 𝑏
Рассмотрим данную модель при условии, что коэффициенты 𝑎 и 𝑐
равноудалены от коэффициента 𝑏 и обозначим Δ𝑟 = 𝑏 − 𝑎 = 𝑐 − 𝑏. Тогда для
данных, приведенных на рисунке 18, функция принадлежности примет вид,
представленный в таблице 1, а для данных приведенных на рисунке 19,
функция принадлежности примет вид, представленный в таблице 2.
Таблица 1. Функции принадлежности для нечеткого отношения «пористость
по ГИС – пористость по керну», построенные по модели конуса с
различными значениями эффективного параметра
Δ𝑟 = 1
Δ𝑟 = 1,25
Δ𝑟 = 1,5
Δ𝑟 = 2
41
Таблица 2. Функции принадлежности для нечеткого отношения «пористость
по керну – остаточная водонасыщенность», построенные по модели конуса с
различными значениями эффективного параметра
Δ𝑟 = 1
Δ𝑟 = 1,25
Результат
свертки
Δ𝑟 = 1,5
двух
нечетких
отношений
Δ𝑟 = 2
примет
вид,
представленный в таблице 3.
Таблица 3. Функции принадлежности для нечеткого отношения «пористость
по ГИС – остаточная водонасыщенность», построенные по функциям
принадлежности с одинаковыми значениями эффективных параметров
Δ𝑟 = 1
Δ𝑟 = 1,25
Экспоненциальная
принадлежности,
модель
которая
Δ𝑟 = 1,5
–
имеет
модель
принцип
для
нахождения
максимальной
Δ𝑟 = 2
функции
энтропии.
Информация, полученная от источников – это итог диффузии, которая
длилась определенное время (𝜏). Базисная система функций задана
уравнением:
𝐾(𝐬 𝑘 , 𝐬) =
где 𝜁𝑘 =
𝜁
min (𝐬 )
√𝜇𝔄
𝒌
1
√𝜋𝜁𝑘
exp [−
|𝐬 𝑘 − 𝐬|
𝜁𝑘 2
2
],
– эффективный параметр рассеяния, изменяющийся
от величины источника, 𝜁 = 2𝑎√𝜏.
42
Функция принадлежности, сконструированная по экспоненциальной
модели, имеет смысл диффузионного рассеяния в бесконечном однородном
фазовом пространстве параметров точечных источников, расположенных в
𝐬 𝑘 [44].
По данным петрофизических отношений, приведенным на рисунке 18,
сконструированы функции принадлежности, представленные в таблице 4, а
по данным петрофизических отношений, приведенным на рисунке 19,
сконструированы функции принадлежности, представленные в таблице 5.
Таблица 4. Функции принадлежности для нечеткого отношения «пористость
по ГИС – пористость по керну», построенные по экспоненциальной модели с
различными значениями эффективного параметра
ζ=1
ζ=1,25
ζ=1,5
ζ=2
Таблица 5. Функции принадлежности для нечеткого отношения «пористость
по
керну
–
экспоненциальной
остаточная
модели
с
водонасыщенность»,
различными
построенные
значениями
по
эффективного
параметра
ζ=1
Результат
ζ=1,25
свертки
двух
ζ=1,5
нечетких
отношений
ζ=2
примет
вид,
представленный в таблице 6.
43
Таблица 6. Функции принадлежности для нечеткого отношения «пористость
по ГИС – остаточная водонасыщенность», построенные по функциям
принадлежности с одинаковыми значениями эффективных параметров
ζ=1
Модель
равномерного
ζ=1,25
обратных
рассеяния
квадратов
ζ=1,5
–
информации
модель,
о
ζ=2
имеющая
достоверности
принцип
параметра
пропорционально обратному квадрату расстояния от измеренного параметра.
Базисная система функций задана уравнением:
𝑟
𝐾(𝐬 𝑘 , 𝐬) =
,
|𝐬 − 𝐬 𝑘 |2 + 𝑟 2
где 𝑟 – сглаживающий параметр.
Построенная функция принадлежности, в основе которой лежит модель
обратных квадратов, означает, что достоверность значения измеренного
параметра распределяется вдоль окружности радиуса 𝑅𝑘 = |𝐬 − 𝐬 𝑘 | + 𝑟. При
удалении от измеренного значения достоверность данного параметра
распространяется на все большие и большие окружности.
Для данных, приведенных на рисунке 18, функция принадлежности
примет вид, представленный в таблице 7, а для данных приведенных на
рисунке 19, функция принадлежности примет вид, представленный в
таблице 8.
44
Таблица 7. Функции принадлежности для нечеткого отношения «пористость
по ГИС – пористость по керну», построенные по модели обратных квадратов
с различными значениями эффективного параметра
r=1
r =1,25
r =1,5
r =2
Таблица 8. Функции принадлежности для нечеткого отношения «пористость
по керну – остаточная водонасыщенность», построенные по модели
обратных квадратов с различными значениями эффективного параметра
r=1
Результат
r =1,25
свертки
двух
r =1,5
нечетких
отношений
r =2
примет
вид,
представленный в таблице 9.
Таблица 9. Функции принадлежности для нечеткого отношения «пористость
по ГИС – остаточная водонасыщенность», построенные по функциям
принадлежности с одинаковыми значениями эффективных параметров
r=1
r =1,25
r =1,5
r =2
45
Результаты
построенных
функций
принадлежности
по
выше
рассмотренным моделям расходятся друг от друга в большей степени при
увеличении эффективных параметров рассеяния Δ𝑟, ζ и r. При увеличении
значений эффективных параметров рассеяния, наиболее адекватными
реальным
данным,
являются
результаты,
полученные
на
базе
экспоненциальной модели. На этом основании, а также учитывая, что
экспоненциальная модель имеет содержание, связанное с нормальным
законом распределения, примем ее за основу при конструировании функций
принадлежности для отношений между параметрами.
2.2.3 Аппроксимация многомерных связей
На основе описанных ранее методов рассмотрим расчет аппроксимации
функции принадлежности на примере трехмерной зависимости
sj, j 3
[137,
138].
На рисунке 22 представлено исходное отношение «абсолютная глубина
(X) – двойное время ВСП (Y) – средняя скорость (Z)», лежащее в основе
проведенного эксперимента. Для наглядности отображения исходные данные
(одновременно измеренные значения параметров) были спроецированы на
плоскость XY.
Рисунок 22 – Исходные данные отношения «абсолютная глубина – двойное
время ВСП – средняя скорость»
46
Рассчитанная карта плотности данных отношения «абсолютная глубина
– двойное время ВСП – средняя скорость», отображающая количество
данных, попавших в выделенный интервал, представлена на рисунке 23.
Рисунок 23 – Карта плотности исходных данных нечеткого отношения
«абсолютная глубина – двойное время ВСП – средняя скорость»
Результатом алгоритма сжатия является построенная карта источников
данных, представленная на рисунке 24.
Рисунок 24 – Карта источников данных для нечеткого отношения
«абсолютная глубина – двойное время ВСП – средняя скорость»
47
По
экспоненциальной
модели
была
рассчитана
функция
принадлежности. Построенная функция принадлежности, для наглядности
ранжированная по достоверности, представлена в таблице 10.
Таблица 10. Функция принадлежности отношения «абсолютная глубина –
двойное время ВСП – средняя скорость», ранжированная по достоверности
α=0.9
α=0.8
α=0.7
α=0.6
48
Продолжение таблицы 10
α=0.5
α=0.4
α=0.3
α=0.2
2.2.4 Описание
алгоритма
конструирования
функции
принадлежности
Все алгоритмы, рассмотренные в данной работе, реализованы в среде
Matlab 2014b. Блок-схемы построены в Visio.
Рассмотрим данный алгоритм на примере экспоненциальной модели,
так как она выбрана наиболее оптимальной.
49
Алгоритм расчета функции принадлежности состоит из следующих
этапов:
Этап 1. Расчет и построение карты плотности:
Шаг 1. Все фазовое пространство параметров 𝑆𝑖 = (𝑥, 𝑦) покрыть
сеткой.
Центральной проблемой служит нахождение разбиения (размера сетки)
(𝑘𝑥 * 𝑘𝑦) с минимальным числом ненулевых ячеек.
Проведя ряд экспериментов, было решено, что оптимальный размер
сетки при условии, что количество измеренных значений 100, должен быть
50*50. Таким образом, размер сетки 𝑘𝑥 * 𝑘𝑦 задается вручную в зависимости
от количества данных (N).
Шаг 1.1. Найти минимум (𝑀𝑖𝑛𝑥, 𝑀𝑖𝑛𝑦) и максимум (𝑀𝑎𝑥𝑥, 𝑀𝑎𝑥𝑦) по
каждому параметру петрофизической модели.
Шаг 1.2. Рассчитать шаг по оси X (первому параметру) и по оси Y
(второму параметру):
𝑑𝑥 =
𝑀𝑎𝑥𝑥−𝑀𝑖𝑛𝑥
𝑘𝑥
; 𝑑𝑦 =
𝑀𝑎𝑥𝑦−𝑀𝑖𝑛𝑦
𝑘𝑦
Шаг 1.3. Рассчитать сетку:
𝑥 = 𝑀𝑖𝑛𝑥: 𝑑𝑥 ∶ 𝑀𝑎𝑥𝑥; 𝑦 = 𝑀𝑖𝑛𝑦: 𝑑𝑦 ∶ 𝑀𝑎𝑥𝑦
Шаг 2. Пройти по всем измеренным значениям параметров и записать,
в какую ячейку сетки они попали (см. Рисунок 25).
Шаг 3.
Вывести графическое изображение полученной карты
плотности.
50
Начало
Исходные
данные
Sj, j=1÷N
N>100
нет
kx = 20;
ky=20
да
Задаем размер
сетки
kx = 50;
ky=50
Находим минимум
и максимум по
каждому
параметру
Minx=Min(Sj,1);
Maxx=Maxx(Sj,1);
Miny=Min(Sj,2);
Maxy=Max(Sj,2)
dx = (Maxx-Minx)/kx;
dy = (Maxy-Miny)/ky
Рассчитываем
шаг по первому
параметру и по
второму
x = Minx:dx:Maxx;
y = Miny:dy:Maxy
Рассчитываем
сетку
j=1..N
i=1..kx
k=1..ky
нет
S(j,1) > x(i) & S(j,1) <= x(i+1)
& S(j,2) > y(k) & S(j,2) <= y(k+1)
да
map_plotn(i,k) =
map_plotn(i,k)+1
Считаем, сколько
данных попало в
каждую ячейку
сетки
map_plotn
Конец
Рисунок 25 – Алгоритм расчета карты плотности данных
51
Этап 2. Расчет и построение карты источников.
Шаг 1. Задать уровень погрешности ε.
Шаг 2. Выполнить итерационный процесс:
Шаг 2.1. Количество источников (𝐾) задать равным количеству ячеек, в
которые попали данные
Шаг 2.2. Решить систему линейных уравнений, в процессе которой
найти значения источников 𝜇𝔄min (𝐬𝒌 ), = 1 ÷ 𝐾
Шаг 2.3. Проверить условие неравенства:
‖𝜇𝔄𝜀 (𝐬) −
1
𝐾
∑ [𝜇𝔄min (𝐬𝒌 ) ∗
√𝜋𝜁0 𝑘=1
exp (−
|𝐬 − 𝐬𝑘 |2
𝜁0 2
)]‖ ≤ 𝜀;
𝐾 → min, 𝜁0 → min.
Шаг 2.4. Если условие выполняется, то количество источников
уменьшить на один и вернуться к шагу 2.2, если не выполняется вернуться к
предыдущему количеству источников, при котором выполнялось условие.
Шаг 3.
Вывести графическое изображение полученной карты
источников.
52
Начало
Карта плотности,
погрешность и размер
сетки
map_plotn, ε, kx, ky
N=0
i=1..kx
j=1..ky
нет
map_plotn(i,j)>0
да
N=N+1;
map_ist(N,2)=map_plotn(i,j);
map_ist(N,3)=map_plotn(i,j)
Считаем количество
ячеек, в которых есть
данные, и начальные
координаты
источников
Задаем начальное
количество источников
K=N+1
K2=K
Запоминаем предыдущее
количество
K=K-1
Уменьшаем количество
источников
A=map_plotn/
max[map_plotn]
Рассчитываем правую
часть системы уравнений
k=1..K
i=1..kx
j=1..ky
B(k,i+3*(j-1))=Exp(-Abs((imap_ist(k,1))^2+(j-map_ist(k,2))^2)/2*ζ0^2
Рассчитываем левую часть
системы уравнений
map_ist(:,1)=B\A
Находим значения в
источниках путем
решения системы
линейных уравнений
методом обратной
матрицы
k=1..ky
i=1..kx
A`=0
j=1..K
A` = A`+map_ist(j,3)*Exp(-Abs((i-map_ist(j,1))^2+(kmap_ist(j,2))^2)/2*ζ0^2
да
Max[map_plotn - (1/
ζ0*Sqrt(2*π))*A`]<=ε
K=K2
Рассчитываем функцию
принадлежности при
минимальном значении
эффективного параметра
по найденным источникам
нет
Если условие не
выполняется возвращаем
предыдущее значение
map_ist
Конец
Рисунок 26 – Алгоритм сжатия информации
53
Этап 3. Расчет функции принадлежности.
Шаг 1. Задать значение эффективного параметра ζ.
Шаг 2. Для каждой ячейки по каждому источнику:
Шаг 2.1. Рассчитать переменный эффективный параметр:
𝜁
𝜁𝑘 =
√𝜇𝔄min (𝐬𝒌 )
Шаг 2.2. Рассчитать поле рассеяния по формуле (экспоненциальной
модели):
𝔄𝜀 (𝐬) =
1
𝐾
∑ 𝜇𝔄min (𝐬𝒌 ) ∗ exp [−
√𝜋𝜁 𝑘=1
|𝐬 𝑘 − 𝐬|
𝜁𝑘 2
2
],
Шаг 3. Рассчитать функцию принадлежности по формуле:
𝔄𝜀 (𝐬)
𝜇𝔄 (𝐬) =
max[𝔄𝜀 (𝐬)]
𝐬
Шаг 4. Вывести графическое изображение полученной функции
принадлежности.
54
Начало
Карта плотности источников,
эффективный параметр, размер сетки
и количество источников
map_ist, ζ, kx, ky, K
k=1..ky
i=1..kx
sum=0
j=1..K
Рассчитываем эффективный
параметр
Рассчитываем функцию
принадлежности в каждой ячейке
по каждому источнику
ζ0=ζ/Sqrt(map_ist(j,3))
sum = sum+map_ist(j,3)*Exp(-Abs((imap_ist(j,1))^2+(k-map_ist(j,2))^2)/2*ζ0^2
func(i,k)=(1/ζ*Sqrt(2*π))*sum
func
Конец
Рисунок 27 – Алгоритм расчета функции принадлежности
2.3.
Заключение
Разработанный метод фазификации отношений между одновременно
измеренными параметрами, основанный на алгоритме сжатия информации,
позволяет снизить размерность величин, а также избавиться от «случайной»
информации, которая может являться причиной неверной интерпретации
данных. В качестве аналитического вида функции принадлежности была
выбрана экспоненциальная модель, имеющая содержание, связанное с
нормальным законом распределения и наиболее адекватно отражающая
исходные
данные.
Построенные
функции
принадлежности
нечетких
отношений, в отличие от статистических методов обработки, отражают
исходные данные, ранжируя их по достоверности.
55
Глава 3.
Анализ нечетких моделей
Результаты
построенных
нечетких
моделей
могут
сильно
варьироваться в зависимости от эффективного параметра рассеяния,
заложенного в модель. В данной главе описывается подход к выбору
оптимального эффективного параметра, а также способ анализа пригодности
построенных функций принадлежности отношений между параметрами для
дальнейшего использования.
Теоретическая
3.1.
основа
анализа
функций
принадлежности
Анализ функции принадлежности подразумевает построение α-сечений
нечеткого отношения, подсчет площади, занимаемой каждым сечением, и
вывод графика зависимости площади от уровня 𝛼 [148].
Система α-сечений функции принадлежности 𝜇(𝐬) строится по
формуле:
𝑆𝛼 = {𝐬 ∈ 𝑆: 𝜇(𝐬) > 𝛼}; 𝑆0 ≃ 𝑆; 𝑆1 = ∅,
где 𝑆𝛼 – площадь α-сечения.
Таким образом, α-сечение – множество принадлежащее функции
принадлежности, достоверность которого выше значения 𝛼.
Мера
этого
множества:
Относительная мера 𝔐𝛼 =
𝑚𝑒𝑠{𝑆𝛼 }
𝑚𝑒𝑠{𝑆0 }
𝑚𝑒𝑠{𝑆𝛼 }
характеризует
его
размеры.
характеризует относительные размеры 𝑆𝛼 .
По мере возрастания параметра 𝛼 𝔐𝛼 убывает и в пределе есть ноль,
которым
характеризуется
пустое
множество.
График
зависимости
относительной меры от α-сечения характеризует степень выраженности,
локализации, контрастности функции принадлежности 𝜇(𝐬). Чем ближе 𝔐𝛼 к
постоянному
значению
в
интервале
значений
𝛼 = [0,1]
тем
более
локализована функция принадлежности и, если 𝔐𝛼 = 1 на всем интервале
𝛼 = [0,1],
𝜇(𝐬)
описывает
зависимость
с
нечетким
отношением
–
функциональную зависимость между параметрами. В этой связи α-анализ,
56
состоящий в изучении функции 𝔐𝛼 , позволяет оценить пригодность 𝜇(𝒔) для
прогноза и правильно выбрать 𝛼-сечения для характеристики интервалов
неопределенности.
Так как для построения функций принадлежности была выбрана
экспоненциальная модель, то должен быть выбран такой эффективный
параметр рассеяния, при котором график зависимости площади от α-сечения
будет иметь характер плавно убывающей кривой. При этом из всех
возможных значений эффективного параметра рассеяния должно быть
выбрано минимальное. Данное обстоятельство объясняется тем, что
эффективный параметр рассеяния влияет на «размытие» данных, которое
повлияет на конечный результат.
3.2.
Алгоритм анализа функций принадлежности
Алгоритм анализа функции принадлежности (уточнение эффективного
параметра рассеяния ζ) выполняется по следующим шагам и представлен на
рисунке 28:
Шаг 1. Задать уровень 𝛼 от 0.1 до 0.9 с шагом 0.1
Шаг 2. Посчитать, какую площадь занимает функция принадлежности
с достоверностью больше, чем 0 (𝑆0 , количество ячеек, в которых значение
функции принадлежности больше нуля):
𝑆0 = {𝐬 ∈ 𝑆: 𝜇(𝐬) > 0};
Этап 3. Для каждого уровня 𝛼:
Этап 3.1. Построить график части функции принадлежности,
достоверность которой больше уровня 𝛼.
Этап 3.2. Посчитать, какую площадь занимает функция
принадлежности с достоверностью больше, чем 𝛼 (количество ячеек, в
которых значение функции принадлежности больше значения 𝛼):
𝑆𝛼 = {𝐬 ∈ 𝑆: 𝜇(𝐬) > 𝛼}.
57
Этап 3. Найти относительную площадь, поделив площадь текущего
α-сечения на площадь нулевого 𝛼-сечения:
𝔐𝛼 =
𝑚𝑒𝑠{𝑆𝛼 }
.
𝑚𝑒𝑠{𝑆0 }
Этап 4. Построить график зависимости относительной площади от αсечения.
Начало
Функция
принадлежности и
размер сетки
func, kx, ky
S=0; S`=0;
func_a=0
i=1..kx
j=1..ky
нет
func(i,j)>0
да
S0=S0+1
Считаем площадь
нулевого а-сечения
a=0,1..0,9
i=1..kx
j=1..ky
нет
func(i,j)>a
да
S`a=S`a+1
func_a(i,j)=a
Sa=S`a/S0
graf(a,1)=Sa;
graf(a,2)=a
Считаем площадь всех
а-сечений
Находим часть функции
принадлежности,
достоверность которой
больше а
Находим
относительную площадь
Расчитываем зависимость
относительной площади
от значения а
graf, func_a
Конец
Рисунок 28 – Алгоритм анализа функции принадлежности
58
В данном случае, на выходе мы получаем двумерный массив graf,
первое измерение которого содержит значения относительной площади, а
второе значение 𝛼. На графике по оси абсцисс откладываем значения 𝛼, а на
оси ординат – значения относительной площади (см. Рисунок 29).
Рисунок 29 – Схема построения зависимости относительной площади от 𝛼
3.3.
Экспериментальная основа α-анализа
Для демонстрации разработанной методики α-анализа проведем
эксперимент по данным, представленным на рисунках 18 и 19.
Результат рассчитанной системы α-сечений для нечеткого отношения
«пористость по ГИС – пористость по керну», построенного по четырем
разным значениям эффективного параметра, представлен в таблице 11, где в
первом столбце перечислены значения параметров 𝛼, меняющиеся от 0 до 1 с
шагом 0.1, а в первой строке эффективные параметры рассеяния: ζ =1, 1.25,
1.5, 2, по которым построены анализируемые функции принадлежности.
Для сравнения динамики изменения площади α-сечения от значения 𝛼
все графики представлены на одном рисунке (см. Рисунок 30).
59
Таблица 11. α-сечения функции принадлежности нечеткого отношения
между параметрами «пористость по ГИС» и «пористость по керну»
ζ=1,25
ζ=1,5
ζ=2
𝛼 > 0,6
𝛼 > 0,5
𝛼 > 0,4
𝛼 > 0,3
𝛼 > 0,2
𝛼 > 0,1
ζ=1
60
𝛼 > 0,9
𝛼 > 0,8
𝛼 > 0,7
Продолжение таблицы 11
Рисунок 30 – График зависимостей площади α-сечения от значения α
функции принадлежности нечеткого отношения «пористость по ГИС –
пористость по керну»
С увеличением значения α относительная мера α-сечения уменьшается,
а при увеличении эффективного параметра рассеяния относительная мера αсечения увеличивается. При ζ=1 график зависимости площади α-сечения от
61
значения α имеет характер ломаной линии. С увеличением данного
параметра график зависимости площади от α-сечения становится более
плавным. При ζ=1,25 график начинает плавно меняться.
Результат рассчитанной системы α-сечений для нечеткого отношения
«пористость по керну – остаточная водонасыщенность», построенного по
четырем разным значениям эффективного параметра, представлен в таблице
12, где в первом столбце перечислены значения параметров 𝛼, меняющиеся
от 0 до 1 с шагом 0.1, а в первой строке эффективные параметры рассеяния: ζ
=1,
1.25,
1.5,
2,
по
которым
построены
анализируемые
функции
принадлежности.
Таблица 12. α-сечения функции принадлежности нечеткого отношения
между
параметрами
«пористость
по
керну»
и
«остаточная
водонасыщенность»
ζ=1,25
ζ=1,5
ζ=2
𝛼 > 0,4
𝛼 > 0,3
𝛼 > 0,2
𝛼 > 0,1
ζ=1
62
𝛼 > 0,9
𝛼 > 0,8
𝛼 > 0,7
𝛼 > 0,6
𝛼 > 0,5
Продолжение таблицы 12
Для сравнения динамики изменения площади α-сечения от значения 𝛼
все графики представлены на одном рисунке (см. Рисунок 31).
Также, как и в предыдущем примере, с увеличением значения α
относительная мера α-сечения уменьшается, а при увеличении эффективного
параметра рассеяния относительная мера α-сечения увеличивается. При ζ=1
график зависимости площади α-сечения от значения α имеет характер
ломаной линии. С увеличением данного параметра график зависимости
площади от α-сечения становится более плавным. При ζ=1,25 график
начинает плавно меняться.
63
Рисунок 31 – График зависимостей площади α-сечения от значения α
функции принадлежности нечеткого отношения «пористость по керну –
остаточная водонасыщенность»
Проведем α-анализ функции принадлежности полученной в результате
композиции
двух
сконструированных
функций
принадлежности
петрофизического отношения пористости по ГИС и пористости по керну и
петрофизического
отношения
пористости
по
керну
и
остаточной
водонасыщенности.
Результат системы α-сечений для рассчитанного нечеткого отношения
«пористость по ГИС – остаточная водонасыщенность», полученного при
использовании
исходных
нечетких
отношений
с
одинаковыми
эффективными параметрами рассеяния, представлен в таблице 13, где в
первом столбце перечислены значения параметров 𝛼, меняющиеся от 0 до 1 с
шагом 0.1, а в первой строке эффективные параметры рассеяния: ζ =1, 1.25,
1.5, 2, по которым построены функции принадлежности, лежащие в основе
анализируемого отношения.
Для сравнения динамики изменения площади α-сечения от значения 𝛼
все графики представлены на одном рисунке (см. Рисунок 32).
64
Таблица 13. α-сечения функции принадлежности нечеткого отношения
между параметрами «пористость по ГИС» и «остаточная водонасыщенность»
ζ=1,25
ζ=1,5
ζ=2
𝛼 > 0,6
𝛼 > 0,5
𝛼 > 0,4
𝛼 > 0,3
𝛼 > 0,2
𝛼 > 0,1
ζ=1
65
𝛼 > 0,9
𝛼 > 0,8
𝛼 > 0,7
Продолжение таблицы 13
Рисунок 32 – График зависимостей площади α-сечения от значения α
функции принадлежности нечеткого отношения «пористость по ГИС –
остаточная водонасыщенность»
При ζ=1 график зависимости площади α-сечения от значения α имеет
характер ломаной линии. С увеличением данного параметра график
зависимости площади от α-сечения становится более плавным. При ζ=1,25
график начинает плавно меняться.
66
3.4.
Заключение
При конструировании функции принадлежности нечетких отношений
центральным
вопросом
является
выбор
оптимального
значения
эффективного параметра. Данную задачу можно решить, воспользовавшись
разработанным методом α-анализа. Данный метод заключается в построении
системы α-сечений, а выбор оптимального значения эффективного параметра
– в анализе графиков зависимости относительной площади 𝛼-сечения от
самого значения 𝛼. В данной главе приведены результаты экспериментов,
отображающие динамику изменения относительной площади каждого 𝛼сечения в зависимости от значения эффективного параметра.
67
Глава 4.
Нечеткий логический вывод Мамдани
В случаях, когда необходимо найти взаимосвязь одного параметра от
другого через промежуточные параметры при том, что исходные отношения
заданы в виде набора одновременно измеренных значений, необходимо
представлять исходные данные в форме нечетких отношений, а затем
использовать алгоритм нечеткого логического вывода Мамдани [134, 140,
142, 143]. В главе приведены: сравнение композиции Мамдани с другими
видами композиций, свойства композиции, а также теоретическая основа
нечеткого логического вывода Мамдани.
Описание алгоритма нечеткого логического вывода
4.1.
Вычисление функции принадлежности геофизического параметра 𝑦 по
исходному параметру 𝑥, заданному на основе фазификации данных 𝔐, и
функции принадлежности для отношения между 𝑥 и 𝑦, определенной
фазификацией данных 𝔄, осуществляется путем построения функции
принадлежности
и
последующего
использования
правил
нечеткого
логического вывода нечетких отношений (традиционной максиминной
нечёткой свёртки или композиции Мамдани) [139]:
𝜇𝔄∗𝔐 (𝐲) = max{min[𝜇𝔐 (x), 𝜇𝔄 (x, y)]} = ∪(𝜇𝔐 (x) ∩ 𝜇𝔄 (x, y))
x
Текст
x
По рассчитанной функции принадлежности 𝜇𝔄∗𝔐 (y) и заданному
нечеткому отношению 𝜇ℜ (y, z) вычисляем:
𝜇𝔄∗𝔐∗ℜ (𝐳) = max{min[𝜇𝔄∗𝔐 (y), 𝜇ℜ (y, z)]} = ∪(𝜇𝔄∗𝔐 (y) ∩ 𝜇ℜ (y, z))
y
Подставляя
Текст
y
функцию
принадлежности
𝜇𝔄∗𝔐 (y)
в
последнее
соотношение получаем:
68
𝜇𝔄∗𝔐∗ℜ (𝒛) = ∪ (∪(𝜇𝔐 (x) ∩ 𝜇𝔄 (x, y)) ∩ 𝜇ℜ (y, z))
y
x
Текст
= ∪ (∪ ((𝜇𝔐 (x) ∩ 𝜇𝔄 (x, y)) ∩ 𝜇ℜ (y, z)))
x
y
= ∪ {[∪(𝜇𝔄 (x, y) ∩ 𝜇ℜ (y, z))] ∩ 𝜇𝔐 (x)} .
x
y
Обозначив:
𝜇𝔄∗ℜ (𝐱, 𝐳) = [∪(𝜇𝔄 (x, y) ∩ 𝜇ℜ (y, z))],
y
(7)
получаем функцию принадлежности 𝜇𝔄∗𝔐∗ℜ (𝐳).
Заменив в соотношении (7) логические операции пересечения и
объединения на эквивалентные операции матричной алгебры, получаем:
𝜇𝔄∗ℜ (𝐱, 𝐳) = max[min(𝜇𝔄 (x, y), 𝜇ℜ (y, z))].
y
(8)
Соотношение (8) известно как композиция Мамдани нечетких
отношений 𝜇𝔄 (x, y) и 𝜇ℜ (y, z), которая имеет смысл подстановки уравнений с
исключением общих повторяющихся переменных [50, 51, 52].
Процедура вычисления 𝜇𝔄∗𝔐 (𝐲) называется композицией нечеткого
отношения и нечеткой величины. Другое ее название – правило нечеткого
логического вывода. Аналогом такой задачи в «четком» случае служит
расчет значений переменной 𝐲, по известной связи между x и 𝐲, и заданной
нечеткой переменной x.
Рассмотрим пример композиции нечеткого отношения «пористость по
ГИС – нефтенасыщенность» и нечеткой величины «пористость по ГИС». На
рисунке 33 представлены исходные функции принадлежности.
69
б
а
Рисунок 33 – Нечеткая величина «пористость по ГИС», при значении
параметра равном 10 (а); нечеткое отношение «пористость по ГИС –
нефтенасыщенность»
Результатом композиции является функция принадлежности нечеткой
величины «нефтенасыщенность», представленная в двух вариантах на
рисунке 34.
а
б
Рисунок 34 – Функция принадлежности прогнозного параметра
«нефтенасыщенность»
Алгоритм нечеткого логического вывода можно разбить на следующие
шаги (см. Рисунок 35):
Шаг 1. Рассчитать функции принадлежности с общим параметром (см.
Рисунок 27):
70
𝜇𝔄 (x, y), 𝜇ℜ (y, z)
Шаг 2. Для каждой пары значений x и z попарно по y сравнить и
выбирать меньшее из двух функций принадлежности. В результате получить
для каждой пары x и z последовательность значений.
Шаг
3.
Для
каждой
пары
значений
x
и
z
выбирать
из
последовательности значений наибольшее, получить требуемое значение
функции принадлежности 𝜇𝔄∗ℜ (𝐱, 𝐳).
Начало
func, kx, ky,
func2, ky2
Две функции
принадлежности и
размеры сеток
F=0
i=1..kx
j=1..ky2
f=0
k=1..ky
f(k)=min[func(i,k)
,func2(k,j)]
F(i,j)=max[f]
Находим минимум
Находим максимум
F
Конец
Рисунок 35 – Алгоритм нечеткого логического вывода Мамдани
71
4.2.
Свойства нечеткого логического вывода
4.2.1 Коммутативность
При
прогнозировании
физико-геологических
параметров
𝒘
по
сейсмическим атрибутам 𝒙 через промежуточные параметры 𝒚, 𝒛, … 𝒉
необходимо выполнять несколько композиций [145]:
𝜇𝔄×𝔐 (𝐱, 𝐳) = max[min(𝜇𝔄 (𝐱, 𝐲), 𝜇𝔐 (𝐲, 𝐳))]
𝐲
Текст
𝜇𝔄×𝔐×ℜ (𝐱, 𝐰) = max[min(𝜇𝔄×𝔐 (𝐱, 𝐳), 𝜇ℜ (𝐳, 𝐰))]
𝐳
Последовательность композиции нечетких отношений определяется
конкретной технологией и может быть различной, но результат при этом не
меняется:
𝜇𝔐×ℜ (𝐲, 𝐰) = max[min(𝜇𝔐 (𝐲, 𝐳), 𝜇ℜ (𝐳, 𝐰))]
𝐳
Текст
𝜇𝔄×𝔐×ℜ (𝐱, 𝐰) = max[min(𝜇𝔄 (𝐱, 𝐲), 𝜇𝔐×ℜ (𝐲, 𝐰))]
𝐲
Таким образом, с помощью нечеткого логического вывода можно
формировать правила вычисления цепочек композиций промежуточных
отношений для получения функции принадлежности итоговой нечеткой
величины и итогового нечеткого отношения [37]. Данное обстоятельство
исключительно важно для формирования правил нечеткой математики при
прогнозировании физико-геологических параметров.
Для проверки гипотезы о коммутативности
композиции была
выполнена свертка нечетких отношений в разной последовательности. В
качестве исходных данных были использованы отношения «пористость по
керну – пористость по ГИС» (см. Рисунок 18), «пористость по керну –
остаточная водонасыщенность» (см. Рисунок 19) и зависимость «остаточная
водонасыщенность – нефтенасыщенность», которая задана уравнением
Кн = 1 − Ков .
Выполним композицию нечетких отношений в следующем порядке:
72
керн
𝜇𝔄×𝔐 (КГИС
[min(𝜇𝔄 (КГИС
п , Ков ) = max Ккерн
п , Кп
п
керн
), 𝜇𝔐 (Кп
, Ков ))]
(см.
Рисунок 36);
ГИС
𝜇𝔄×𝔐×ℜ (КГИС
п , Кн ) = max Ков [min(𝜇𝔄 (Кп , Ков ), 𝜇𝔐×ℜ (Ков , Кн ))]
(см.
Рисунок 37).
а
б
в
Рисунок 36 – Построение промежуточной функции принадлежности: первая
функция принадлежности «пористость по ГИС – пористость по керну» (а);
вторая функция принадлежности «пористость по керну – остаточная
водонасыщенность» (б); промежуточная функция принадлежности
«пористость по ГИС – остаточная водонасыщенность» (в)
73
а
б
в
Рисунок 37 – Построение итоговой функции принадлежности:
промежуточная функция принадлежности «пористость по ГИС – остаточная
водонасыщенность» (а); третья функция принадлежности «остаточная
водонасыщенность – нефтенасыщенность» (б); итоговая функция
принадлежности «пористость по ГИС – нефтенасыщенность» (в)
Затем
проведем
свертку
нечетких
отношений
в
иной
последовательности:
керн
𝜇𝔐×ℜ (Кп
керн
, Кн ) = maxКов [min(𝜇𝔐 (Кп
, Ков ), 𝜇ℜ (Ков , Кн ))]
(см.
Рисунок 38);
керн
𝜇𝔄×𝔐×ℜ (КГИС
[min(𝜇𝔄 (КГИС
п , Кн ) = max Ккерн
п , Кп
п
керн
), 𝜇𝔐×ℜ (Кп
, Кн ))]
(см. Рисунок 39).
74
б
а
в
Рисунок 38 – Построение промежуточной функции принадлежности: первая
функция принадлежности «пористость по керну – остаточная
водонасыщенность» (а); вторая функция принадлежности «остаточная
водонасыщенность – нефтенасыщенность» (б); промежуточная функция
принадлежности «пористость по керну – нефтенасыщенность» (в)
75
а
б
в
Рисунок 39 – Построение итоговой функции принадлежности: третья
функция принадлежности «пористость по ГИС – пористость по керну» (а);
промежуточная функция принадлежности «пористость по керну –
нефтенасыщенность» (б); итоговая функции принадлежности «пористость по
ГИС – нефтенасыщенность» (в)
Результаты эксперимента по построению функций принадлежности
нечеткого отношения «пористость по ГИС – нефтенасыщенность» в разных
последовательностях (см. рисунок 37в и рисунок 39в), подтверждают, что
композиция
Мамдани
обладает
свойством
коммутативности.
При
вычислении функции от каждой из переменных, входящих в нечеткое
отношение, порядок композиции не имеет значения, т.е. отношение является
инвариантом, что особо важно при формировании графов прогноза
параметров.
76
4.2.2 Возрастание неопределенности. Аналог подстановки
уравнений
Прогноз параметров выполняется на основе построенных функций
принадлежности по правилам нечеткого логического вывода (композиции
Мамдани). За основу построения функций принадлежности была выбрана
экспоненциальная модель, обладающая максимальной энтропией [25].
Достоверность информации от точки измерения значения параметра к точке
оценивания достоверности в такой модели распространяется подобно
процессу диффузии. Из-за этого, возникает процесс возрастания уровня
нечеткости результатов при увеличении длины цепочек из элементов
композиции [144, 147].
Продемонстрируем данный процесс на примере. Для эксперимента
использовались тестовые данные: петрофизическое отношение пористости
по ГИС и пористости по керну, петрофизическое отношение пористости по
керну
и
остаточной
керн
уравнениями: Кп
соответственно
и
водонасыщенности,
которые
были
заменены
керн
= 1.064 ∗ КГИС
− 0.275 и Ков = 88.15 ∗ 𝑒 (−0.159∗Кп
п
петрофизическая
зависимость
)
остаточной
водонасыщенности от нефтенасыщенности Кн = 1 − Ков .
Выполним композицию петрофизической зависимости пористости по
ГИС и пористости по керну и петрофизической зависимости пористости по
керну и остаточной водонасыщенности (см. Рисунок 40).
Затем, выполним композицию полученного результата и зависимости
«коэффициент
остаточной
водонасыщенности
–
коэффициент
нефтенасыщенности». Результат представлен на рисунке 41.
77
б
а
в
Рисунок 40 - Расчет функции принадлежности петрофизического отношения
пористости по ГИС и остаточной водонасыщенности (в) с помощью
композиции нечеткого петрофизического отношения пористости по ГИС и
пористости по керну (а) и нечеткого петрофизического отношения
пористости по керну и остаточной водонасыщенности (б)
а
б
в
Рисунок 41 - Расчет функции принадлежности петрофизического отношения
пористости по ГИС и остаточной водонасыщенности (в) с помощью
композиции нечеткого петрофизического отношения пористости по ГИС и
пористости по керну (а) и нечеткого петрофизического отношения
пористости по керну и остаточной водонасыщенности (б)
78
По результатам проведенных экспериментов видно, что при удлинении
цепочки выполненных композиции нечетких отношений, диффузия как мера
размытости результата композиции, а, следовательно, мера нечеткости
данных - возрастает.
Если воспользоваться подстановкой уравнений одно в другое для
исключения общих переменных, остаточную водонасыщенность можно
ГИС +0.043725
выразить уравнением Ков = 88.15 ∗ 𝑒 −0.169176Кп
. График данного
уравнения имеет вид, представленный на рисунке 42.
Рисунок 42 – График зависимости остаточной водонасыщенности от
пористости по ГИС
Уравнение зависимости нефтенасыщенности от пористости по ГИС
примет
вид
ГИС +0.043725
Кн = 1 − 88.15 ∗ 𝑒 −0.169176Кп
.
График
данной
зависимости представлен на рисунке 43.
79
Рисунок 43 – График зависимости нефтенасыщенности от пористости по
ГИС
Результаты
композиции
нечетких
отношений
зависимостей
и
подстановки уравнений дают аналогичные результаты: сохраняется динамика
изменения зависимости (см. Рисунок 40, Рисунок 42 и Рисунок 41, Рисунок
43). Данное обстоятельство позволяет использовать нечеткий логический
вывод для получения отношения итоговых параметров с исходными
параметрами.
4.3.
Проблема
Прогнозирование емкостных параметров
прогнозирования
фильтрационно-емкостных
свойств
и
параметров нефтегазоносности по сейсмическим атрибутам оказывается
крайне актуальной и научно значимой, из-за того, что она лежит в основе
подсчета запасов углеводородов [35]. Данная проблема подразделяется на
ряд этапов. Первым этапом является определение параметров, наилучшим
образом обеспечивающих их пригодность для прогноза фильтрационноемкостных свойств по данным на эталонной площади, которые оцениваются
величиной ошибки прогноза. Данная процедура может выполняться с
применением
факторного
анализа
[91].
Эффективность
метода
в
80
значительной степени заключается в правильном подборе комбинаций
сейсмических параметров, и в настоящее время данному вопросу уделяется
большое внимание [92, 93, 94, 98, 99, 100, 105]. Данные методы также
применяются для решения близких задач, которые связанны с картированием
зон трещиноватости [95]. Вторым этапом является определение и настройка
обучающих правил прогноза на эталонной площади [96,97]. Фактически,
данный этап задает метод прогнозирования, который далее используется
[146].
После выбора параметров, и настройки обучающих правил прогноза под
конкретную
обучающую
выборку,
появляется
задача
сохранения
неопределенности исходных данных в результатах прогноза параметров,
выражающуюся в разбросе измеренных значений параметра по отношению к
получаемым значениям в результате прогноза. Данное обстоятельство
наиболее рельефно выражается при использовании регрессионных моделей.
Рассеяние точек по отношению к линии регрессии является компонентой,
которую невозможно устранить, и присутствует в конечных результатах
прогноза фильтрационно-емкостных свойств даже тогда, когда оно не
проявляется отчетливо. Без сомнения в конечных результатах прогноза будет
присутствовать неопределенность, которую стоило бы учитывать. Еще одной
проблемой является то обстоятельство, что параметры расположены не
равномерно в обучающей выборке. В одних областях значений параметров
заключается большое количество данных, а в других областях плотность
данных меньше. То обстоятельство, если для разных интервалов принимать
одно и то же правило – регрессионную модель – приводит к возникновению
неконтролируемых ошибок, которые чреваты большими последствиями.
Лишенный этих недостатков является подход, основанный на
технологии нечеткого моделирования состоящий в том, что данные для
обучения представляются в форме отношений нечетких величин, которые в
последующем используются при нечетком логическом выводе Мамдани [19].
Использование данной технологии позволяет получать согласованные с
81
фактической неопределенностью результаты прогноза, наследующие в себе
различный уровень достоверности данных по всем интервалам обучающей
модели.
Современное
развернутое
описание
технологии
нечеткого
моделирования представлено в работе [27].
Согласно технологии нечеткого моделирования нечеткая величина
представляется в форме функции принадлежности 0 1 , которая имеет
значение меры достоверности – возможности того, что определение этой
величины будет равно значению . Функции принадлежности должна
соответствовать условию нормировки:
отличием
функции
max 1 ,
принадлежности
от
что является ключевым
вероятностной
меры.
Это
обстоятельство обеспечивает возможность организации вычислительных
процедур, основанных на теории нечеткого логического вывода, а также
композиции нечетких отношений с целью моделирования конечных функций
принадлежности физико-геологических моделей.
Основой алгебры для нечетких величин, составляющей вычислительную
основу конструирования, функций принадлежности служат операции:
объединение x x max x , x ;
пересечение x x min x , x ;
разность x x x , x x ;
симметрическая разность x x .
Для
получения
одномерных
двухмерной
функции
используется
принадлежности
операция
из
двух
импликации
, min , .
Двухмерная функция принадлежности характеризует отношение между
двумя одномерными нечеткими переменными. Две двухмерные функции
принадлежности x, y ; y, могут быть использованы при композиции:
x, x, y y, max min x, y , y, ,
(9)
y
82
которая имеет смысл матричного умножения, с заменой операции
алгебраического на логическое умножение – пересечения, а операции
алгебраического на логическое суммирования – объединение. Получается,
выражение (9) является аналогом подстановки одного уравнения в другое,
исключающим одну переменную (в данном случае
y ). По аналогии
выполняется композиция одномерной и двухмерной нечеткой величины:
x x, max min x, , x ,
(10)
x
аналогичная
расчету
из
уравнения
двух
переменных
x, значения
переменной по известному значению величины x . Правило (10) - это
композиция Мамдани, реализующая правило нечеткого логического вывода о
значениях нечеткой переменной (прогнозной) , по отношению между
двумя нечеткими переменными x, , и функцией принадлежности x
для нечеткой переменной x - аргумента прогноза.
Представление
данных
по
эталонному
в
форме
функций
принадлежности для отношений между нечеткими параметрами разработано
в работе [136], где приведены и правила логического вывода для прогноза
параметров.
Рассмотрим некоторые примеры.
Пример 1. Для данных, по одной из тестовых площадей определяющих
зависимости «скорость продольных волн и коэффициент пористости»,
«коэффициент
пористости
и
коэффициент
нефтенасыщенности»,
представленных на рисунке 44, были построены функции принадлежности,
определяющие соответствующие нечеткие отношения, которые изображены
на рисунке 45.
83
Рисунок 44 – исходные данные: А – зависимость «скорость продольных волн
v – коэффициент пористости kп », Б – зависимость «коэффициент
пористости – коэффициент нефтенасыщенности kн »
Обратим внимание на проявляющуюся в данных (см. Рисунок 44) явную
неоднородность по степени их группируемости. Это влечет за собой
различия в достоверности фиксируемых функцией принадлежности (см.
Рисунок 45). На этом же рисунке приведен результат композиции
(соотношение (10)) в соответствии с правилом логического вывода,
исключающий промежуточный параметр – коэффициент пористости, и
устанавливает неизвестное ранее отношение между нечеткой переменной
«коэффициент
нефтенасыщенности
и
нечеткой
переменной
скорость
продольных волн».
84
Рисунок 45 – расчет зависимости «скорость продольных волн – коэффициент
нефтенасыщенности»
А – поле рассеяния зависимости «скорость продольных волн –
коэффициент пористости»
Б – поле рассеяния зависимости «коэффициент пористости –
коэффициент нефтенасыщенности»
В – поле рассеяния зависимости «скорость продольных волн –
коэффициент нефтенасыщенности»
Как видно, из приведенного результата, данные характеризуются
существенной
неоднородностью,
что
необходимо
учитывать
для
дальнейшего прогноза.
В таблице 14 приведены сечения найденного отношения, определяющие
значения нечеткой переменной «коэффициент нефтенасыщенности» для
разных значений скорости продольных волн.
85
Таблица
14.
Прогнозирование
нечеткой
переменной
«коэффициент
нефтенасыщенности»
Скорость = 2.7
Скорость = 3
Скорость = 3.2
При прогнозировании нечеткой переменной также имеют существенное
влияние неопределенности исходных данных, и учет влияния на результат
обязателен.
Пример 2. Тестовые данные, определяющие зависимости «плотность и
пористость», «пористость и скорость продольных волн», представлены на
рисунке 46. Соответствующие им функции принадлежности, определяющие
нечеткие отношения приведены на рисунке 47.
Рисунок 46 – исходные данные: А – зависимость «плотность ( ) –
пористость kп », Б – зависимость «пористость – скорость»
86
Рисунок 47. А – поле рассеяния зависимости «плотность - пористость»
Б – поле рассеяния зависимости «пористость - скорость»
Композиция двух
нечетких
отношений
для исключения
общей
переменной в данном случае физического смысла не имеет, поскольку
прогнозным
параметром
служит
пористость.
Рассмотрим
задачу
прогнозирования пористости при заданном распределении двух параметров
(плотность, скорость).
Построение функции принадлежности трехмерной зависимости –
скорость, плотность, и им соответствующее пересечение пористости может
быть выполнено по формуле:
v, k п , min v, k п , k п ,
(11)
Последнее соотношение следует понимать как импликацию v, k п и
k п , ,
и рассчитывать следующим образом. Фиксируется тройка
параметров v, k п , , для которой значением функции принадлежности
служит наименьшее из v, k п , y, .
Для нахождения пористости необходимо найти пересечение срезов
нечетких отношений зависимостей «плотность - пористость» и «пористость скорость». В таблице 15 приведены срезы с максимальным пересечением.
87
Таблица 15. Пересечение срезов нечетких отношений зависимостей
«плотность - пористость» и «пористость - скорость»
𝑣 = 5.2, 𝜎 = 2.53
𝑣 = 4.6, 𝜎 = 2.4
𝑣 = 4, 𝜎 = 2.1
В связи с тем, что зависимость пористости имеет вид убывающей
функции относительно
скорости
и плотности, представим срез по
рассчитанной функции принадлежности плоскостью, заданной уравнением
kп
55
55
1922
v
3
7
21
(см. Рисунок 48).
Рисунок 48 – Срез функции принадлежности зависимости «плотность –
скорость – пористость»
Для более наглядного представления были выполнены также сечения
функции принадлежности плоскостями, параллельными плоскостям XY, YZ,
XZ (см. Рисунок 49).
88
Рисунок 49. А – Срез функции принадлежности плоскостями, параллельными
плоскости XY; Б – Срез функции принадлежности плоскостями,
параллельными плоскости YZ; В – Срез функции принадлежности
плоскостями, параллельными плоскости XZ.
Для рассмотрения структуры нечеткого отношения трехмерной
зависимости
выполним
M v,k п , : v,k п ,
построение
α-сечений
(см. Таблица 16).
89
Таблица 16. Последовательность α-сечений
α=0.8
α=0.5
α=0.3
α=0.1
Выводы
Представление исходных данных в форме нечетких величин, а связей
между параметрами в форме нечетких отношений, позволит, используя
правила нечеткого логического вывода и правила импликации, оснащенные
построением
сечений для использования независимых параметров,
выполнять прогноз параметров, объективно учитывающих особенности
реальных данных
Технология
нечеткого
моделирования
при
прогнозе
параметров
позволит получить дифференцированные по достоверности результаты
прогноза,
согласованные
с
реальными
данными,
содержащими
неопределенность.
90
Альтернативные виды композиций
4.4.
Логически непротиворечиво и математически корректно рассчитывать
отношения итоговых параметров с исходными
многими способами.
Рассмотрим результаты, полученные при использовании различных видов
композиций.
В качестве исходных данных используем петрофизическое отношение
пористости по ГИС и пористости по керну и петрофизическое отношение
пористости по керну и остаточной водонасыщенности. Исходные данные
петрофизических отношений и соответствующие им сконструированные
керн
функции принадлежности 𝜇𝔄 (КГИС
п , Кп
керн
), 𝜇𝔐 (Кп
, Ков ) представлены на
рисунках 50 и 51 соответственно.
б
а
Пористость по керну, %
18
16
14
12
10
8
6
4
2
0
0
2
4
6
8
10
12
14
16
Пористость по ГИС, %
Рисунок 50 – Нечеткое отношение «пористость по ГИС – пористость по
керну»: исходные данные петрофизического отношения (а), функция
принадлежности (б)
а
б
91
Остаточная
водонасыщенность, %
100
80
60
40
20
0
0
2
4
6
8 10 12 14 16 18 20
Пористость по керну, %
Рисунок 51 – Нечеткое отношение «пористость по керну – остаточная
водонасыщенность»: исходные данные петрофизического отношения (а),
функция принадлежности (б)
Тогда композиция Мамдани 𝜇𝔄×𝔐 (КГИС
п , Ков ) приводит к результату:
Рисунок 52 – Композиция Мамдани
Используя Max-prod композицию:
керн
𝜇𝔄×𝔐 (КГИС
{[𝜇𝔄 (КГИС
п , Ков ) = max
п , Кп
керн
Кп
керн
) ∗ 𝜇𝔐 (Кп
, Ков )]},
получаем результат:
92
Рисунок 53 – Max-prod композиция
Min-max композиция:
керн
𝜇𝔄×𝔐 (КГИС
[max(𝜇𝔄 (КГИС
п , Ков ) = min
п , Кп
керн
Кп
керн
), 𝜇𝔐 (Кп
, Ков ))]
приводит к результату:
Рисунок 54 – Min-max композиция
Min-average композиция:
керн
𝜇𝔄×𝔐 (КГИС
[𝜇𝔄 (КГИС
п , Ков ) = 0.5 ∙ max
п , Кп
керн
Кп
керн
) + 𝜇𝔐 (Кп
, Ков )]
приводит к результату:
93
Рисунок 55 – Min-average композиция
Min-min композиция
керн
𝜇𝔄×𝔐 (КГИС
[min(𝜇𝔄 (КГИС
п , Ков ) = min
п , Кп
керн
Кп
керн
), 𝜇𝔐 (Кп
, Ков ))]
приводит к результату:
Рисунок 56 – Min-min композиция
Max-max композиция
керн
𝜇𝔄×𝔐 (КГИС
[max(𝜇𝔄 (КГИС
п , Ков ) = max
п , Кп
керн
Кп
керн
), 𝜇𝔐 (Кп
, Ков ))]
приводит к результату:
94
Рисунок 57 – Max-max композиция
Данные эксперименты приводят к выводу о явном преимуществе
композиции Мамдани в сравнении с другими возможными композициями.
Именно эта композиция, теоретически обоснованная адекватности
алгебраическим операциям, берется за основу в конструировании методов
нечеткого прогнозирования.
Эти, и многочисленные другие экспериментальные исследования,
допустимых внутренней логикой алгебры нечетких величин, композиций
лежащих в основе нечеткого логического вывода, показали несомненное
преимущество и адекватность задачам прогноза фильтрационно – емкостных
параметров композиции Мамдани.
4.5.
Заключение
Для прогнозирования зависимости одних параметров по другим через
промежуточные параметры целесообразно представлять исходные данные в
виде нечетких моделей и использовать алгоритм нечеткого логического
вывода
(композицию
Мамдани).
Композиция
Мамдани
адекватна
подстановке уравнений одно в другое с целью избавления от общих
переменных, а также имеет существенные преимущества перед другими
видами композиций. При удлинении цепочки промежуточных параметров
неопределенность результата возрастает, поэтому при равных условиях
95
лучше выбирать более короткую цепочку, при этом порядок композиции не
имеет значения. Также алгоритм нечеткого логического вывода можно
использовать для вычисления нечеткой прогнозной величины. Данная
процедура называется композицией нечеткого отношения и нечеткой
величины.
96
Глава 5.
Информационная экспертиза
Информационно-экспертная
оценка
результатов
промыслово-
геологического моделирования в нефтегазовой геологии заключается в
контроле достоверности построенных моделей, анализируя их соответствие
имеющимся исходным петрофизическим данным [149].
Построение геолого-промысловой модели месторождения, пригодной
для последующего использования при проектировании схемы разработки
месторождения основано на структурно-фациальной и тектонической модели
месторождения
[35,
наполненной
74]
содержательными
физико-
геологическими параметрами, такими как нефтегазонасыщенность, фазовая
проницаемость, общая и динамическая пористость и так далее [46, 47, 48].
Итогом служит
куб
физико-геологических
параметров,
в
частности
подсчетных параметров для оценки запасов месторождения углеводородов
(проницаемость, нефтегазонасыщение и т.д.). Этап создания структурно
фациальной и тектонической основы, обеспечивающий геометрическую
основу формирования фациальных сегментов в кубе параметров основан на
результатах предшествующих зональных и локальных геологических
изысканий, сейсмических методов, дополненные данными гравиметрии и
других полевых методов с привлечением скважинных данных для
корреляции структурно-фациальных границ. Наполнение куба модели
содержательными физико-геологическими параметрами это второй этап. Он
основан на использовании предварительно построенных петрофизических
моделей коллектора [18, 36] относительно рассматриваемых параметров и
прогнозирования
требуемых
преимущественно
геофизических
Различная
степень
петрофизической
параметров
исследований
обеспеченности
модели,
ее
в
реальными
возможная
кубе
по
скважин
данными
неполнота
для
результатам
(ГИС)
[22].
фрагментов
отдельных
измеренных и используемых при прогнозировании значений параметров ГИС
приводит к тому, что объективно различные фрагменты прогнозируемой
модели (куба параметров) определяются с различной степенью надежности,
97
допускают многовариантное толкование. Не учет этих обстоятельств
повышает технико-экономические риски принятия решений о проекте
разработки
залежей
[23]
и
требует
объективной
оценки,
меры
информационной обеспеченности выполненных построений различных
фрагментов. При всем доверии к опыту и интуиции геолога-интерпретатора
требуется
ранжировать
компоненты
модели
по
обеспеченности
объективными данными, указывая на долю интерполяции и субъективизма.
Решение этой задачи требует нового подхода и нового аппарата анализа
данных.
На основе нечеткого моделирования создан аппарат информационного
анализа, который отображает данные, используемые для построения куба
физико-геологической модели, в форме, ранжированной по достоверности.
Для реализации данного аппарата было введено понятие нечеткой
петрофизической
модели,
которая
участвует
в
построении
поля
достоверности на основе нечеткого логического вывода преобразования
измеренных параметров в прогнозные. Разработанный аппарат включает
анализ сечений поля достоверности экспертируемой модели.
Метод пространственного продолжения
5.1.
Замеры
скважине,
физико-геологических
при
достоверности
построении
прогнозных
параметров
пространственного
параметров
осуществляются
распределения
(коэффициент
по
поля
пористости,
коэффициент нефтенасыщенности) необходимо рассчитывать значения вдали
от нее на расстоянии ∆R. То есть, для каждого пропластка функция
принадлежности определена там, где проходит скважина. Таким образом,
возникает
задача
интерполяции
функции
принадлежности
по
пространственным координатам, для решения которой была разработана
технология
информацию
аналитического
о
продолжения,
ранжированной
позволяющая
достоверности
значений
получать
физико-
геологических параметров в исходной геологической модели.
98
При удалении от скважины информация должна убывать до
критического расстояния (𝑅𝑘). Предполагается, что дальше этого расстояния
мы не можем оценивать достоверность. При этом значение достоверности в
пределах допустимой области должно меняться по формуле:
μ(z, ∆R) = max min[μ(x, z, ξ1 (∆R)), μ(x, ξ2 (∆R))]
x
где ξ(∆R) – эффективный параметр рассеяния, ∆R – расстояние от
скважины до ячейки, для которой рассчитывается значение.
В точке замера (скважине) эффективный параметр должен быть
оптимальным, по которому были построены петрофизические модели, и
убывать до
1
ξ .
3 0
Следовательно, эффективный параметр рассеяния при
удалении от скважины рассчитывается по формуле:
ξ(∆R) =
ξ0
2∆R
1+
Rk
где ξ0 – оптимальный эффективный параметр рассеяния, при котором
были построены функции принадлежности, выбранный на этапе α-анализа.
В случае, когда скважина одна, критическое расстояние выбирается
равное одному километру, если скважин несколько критическое расстояние
рассчитывается по формуле:
m
1
Rk = ∑ Ri
m
i=1
где R i – расстояние между парой скважин, m – количество скважин.
Алгоритм аналитического продолжения для получения информации о
достоверности значений физико-геологических параметров в физикогеологической модели:
Входные данные: сетка – матрица, состоящая из 0 и 1, отображающая
структуру пропластка, содержащая координаты ячеек сетки; координаты
скважин; измеренные значения геофизического и принятого прогнозного
99
параметра по каждой скважине для каждого пропластка; рассчитанные карты
источников; значения эффективных параметров для каждого отношения.
Шаг 1. Найти координаты ячеек сетки x_s, y_s.
Шаг 2. Найти количество скважин.
Шаг 3. Рассчитать критическое расстояние от скважины (Rk), после
которого достоверность становится равной 0.
В случае, когда скважина одна, критическое расстояние выбирается
равное одному километру, если скважин несколько критическое расстояние
рассчитывается по формуле:
𝑚
1
𝑅𝑘 = ∑ 𝑅𝑖
𝑚
𝑖=1
где 𝑅𝑖 – расстояние между парой скважин, 𝑚 – количество скважин.
Шаг 4. Для каждой ячейки сетки:
Шаг 4.1 Если значение исходной сетки равно 0, то достоверность
приравнять 0 и перейти к следующей ячейке. Если значение исходной сетки
равно 1, то перейти к следующему шагу.
Шаг 4.2 Для каждой скважины:
Шаг 4.2.1 Рассчитать расстояние от текущей скважины до ячейки.
Шаг 4.2.2 Если расстояние от текущей скважины до ячейки больше
критического расстояния, то достоверность приравнять 0, перейти к
следующей скважине. Если расстояние от текущей скважины до ячейки
меньше критического расстояния, то перейти к следующему шагу.
Шаг 4.2.3 Рассчитать значение эффективного параметра, который
изменяется по формуле:
𝜉(∆𝑅) =
𝜉0
2∆𝑅
1+
𝑅𝑘
где 𝜉0 – оптимальный эффективный параметр рассеяния, при котором
были построены функции принадлежности, выбранный на этапе α-анализа;
100
∆𝑅 – расстояние от ячейки, для которой рассчитывается значение, до
ближайшей скважины.
Шаг 4.2.4 Рассчитать функции принадлежности (см. Рисунок 27) с
данным эффективным параметром
Шаг 4.2.5 Выполнить композицию (цепочку композиций) всех
функций принадлежности (см. Рисунок 35).
Шаг 4.2.6 Получить функцию принадлежности для прогнозного
параметра, выполнив композицию полученного нечеткого отношения и
геофизического параметра.
Шаг 4.2.7 По значению принятого прогнозного параметра найти
распределение достоверности от каждой скважины.
Шаг 4.3 Найти распределение достоверности каждой ячейке как
максимум от всех скважин.
Шаг
5.
Вывести
графическое
отображение
распределения
достоверности по пропластку.
Начало
Сетка пропластка, координаты скважин, значения измеренного
геофизического параметра и принятого прогнозного парметра, карты
источников исходных данных, оптимальные эффективные параметры,
диапазон значений геофизического параметра и прогнозного
setka, koord_skv,
par_x_sloi, par_z_sloi,
map_ist, ζ0, x, z
x_s=setka(1,:)
y_s=setka(:,1)
setka=setka(2:y_s,2:x_s)
K=kol(koord_skv)
нет
Находим координаты ячеек сетки
пропластка, удаляем координаты ячеек
из сетки, находим количество скважин
да
K>1
i=1..K-1
Rk=1000
Рассчитываем Rk: если одна скважина,
Rk=1000, если несколько, Rk =
среднему расстоянию между
скважинами
j=i+1..K
R=R+Sqrt((koord_skv(i,1)koord_skv(j,1))^2+(koord_skv(i,2)-(koord_skv(j,2))^2)
Rk=R/K
1
Рисунок 58 – Начало алгоритма аналитического продолжения
101
1
i=1..kol(x_s)
j=1..kol(y_s)
да
нет
Setka(i,j)==1
k=1..K
R=Sqrt((koord_skv(k,1)x_s(i))^2+(koord_skv(k,2)-y_s(j))^2)
Находим расстояние между скважиной
и ячейкой
да
нет
R<Rk
Рассчитываем значение эффективного
параметра
ζ=ζ0/(1+2*R/Rk)
func=Функция принадлежности(map_ist,ζ)
Пересчитываем функции
принадлежности
F=Композиция(func)
Пересчитываем композиции
F_prognoz=Композиция(F,param_x_sloi)
Находим распределения
достоверности прогнозного
параметра
l=1..kol(z)
нет
dost(k)=0
z(l)=
=param_z_sloi
да
dost(k)=F_prognoz(l)
rasp_dost(i,j)=0
rasp_dost(i,j)=max(dost)
Определяем значение
достоверности принятого
прогнозного параметра
Определяем достоверность в
каждой ячейке сетки
rasp_dost
Конец
Рисунок 59 – Конец алгоритма аналитического продолжения
Процесс расчета достоверности по всему пропластку занимает большое
количество времени, т.к. для каждой ячейки сетки пересчитываются
исходные петрофизические модели, и выполняется композиция данных
моделей. Для оптимизации процесса выполним интерполяцию достоверности
распределения параметра насыщенности по формуле:
102
𝜇(𝑧, ∆𝑅) = 𝜇(𝑧) ∗ 𝑒
−(
∆𝑅 𝑛
) ∗𝑘
𝑅𝑘
,
где 𝑛 = 2 и 𝑘 = 1
Данный процесс состоит из следующих шагов:
Шаг 1. Найти координаты ячеек сетки 𝑥_𝑠, 𝑦_𝑠.
Шаг 2. Найти количество скважин.
Шаг 3. Рассчитать критическое расстояние от скважины (Rk), после
которого достоверность становится равной 0.
В случае, когда скважина одна, критическое расстояние выбирается
равное одному километру, если скважин несколько критическое расстояние
рассчитывается по формуле:
𝑚
1
𝑅𝑘 = ∑ 𝑅𝑖
𝑚
𝑖=1
где 𝑅𝑖 – расстояние между парой скважин, 𝑚 – количество скважин.
Шаг 4. Для каждой ячейки сетки:
Шаг 4.1 Если значение исходной сетки равно 0, то достоверность равна
0 и переходим к следующей ячейке. Если значение исходной сетки равно 1,
то переходим к следующему шагу.
Шаг 4.2 Для каждой скважины:
Шаг 4.2.1 Считаем расстояние от текущей скважины до ячейки.
Шаг 4.2.2 Если расстояние от текущей скважины до ячейки больше
критического расстояния, то достоверность приравнять 0, переходим к
следующей скважине. Если расстояние от текущей скважины до ячейки
меньше критического расстояния, то перейти к следующему шагу.
Шаг 4.2.3 По значению принятого прогнозного параметра найти
распределение достоверности по формуле:
𝜇(𝑧, ∆𝑅) = 𝜇(𝑧) ∗ 𝑒
−(
∆𝑅 2
)
𝑅𝑘
Шаг 4.3 Найти распределение достоверности каждой ячейке как
максимум от всех скважин.
103
Шаг
5.
Вывести
графическое
отображение
распределения
достоверности по пропластку.
Начало
Сетка пропластка, координаты скважин, значения измеренного
геофизического параметра и принятого прогнозного парметра, карты
источников исходных данных, оптимальные эффективные параметры,
диапазон значений геофизического параметра и прогнозного
setka, koord_skv,
par_x_sloi, par_z_sloi,
map_ist, ζ0, x, z
x_s=setka(1,:)
y_s=setka(:,1)
setka=setka(2:y_s,2:x_s)
K=kol(koord_skv)
нет
Находим координаты ячеек сетки
пропластка, удаляем координаты ячеек
из сетки, находим количество скважин
Рассчитываем Rk: если одна скважина,
Rk=1000, если несколько, Rk =
среднему расстоянию между
скважинами
да
K>1
i=1..K-1
Rk=1000
j=i+1..K
R=R+Sqrt((koord_skv(i,1)koord_skv(j,1))^2+(koord_skv(i,2)-(koord_skv(j,2))^2)
Rk=R/K
i=1..kol(x_s)
j=1..kol(y_s)
да
нет
Setka(i,j)==1
k=1..K
R=Sqrt((koord_skv(k,1)x_s(i))^2+(koord_skv(k,2)-y_s(j))^2)
Находим расстояние между скважиной
и ячейкой
да
нет
R<Rk
dost(k)=0
rasp_dost(i,j)=0
dost(k)=param_z_sloi*exp(-(R/Rk)^2)
rasp_dost(i,j)=max(dost)
Определяем распределение
достоверности по каждой скважине
Определяем достоверность в
каждой ячейке сетки
rasp_dost
Конец
Рисунок 60 – упрощенный алгоритм аналитического продолжения
104
5.2.
Дефазификация
Дефазификация
состоит
в
представлении
полученного
поля
достоверности 𝜇(𝐳(𝐑)) в форме привычной для геологических приложений –
в
качестве
распределения
параметра
геологической
модели
𝐳(𝐑).
Дефазификация - это визуализация и представление поля достоверности в
форме
числовых
и
графических
объектов,
характеризующих
поле
достоверности. Дефазификация может быть выполнена многими различными
способами, основанными на желаемых к визуализации компонент модели
геологической среды.
Для каждого значения 𝐳(𝐑) функция принадлежности 𝜇(𝐳(𝐑)) задает
нечеткую функцию, которая характеризует распределение достоверности
данного параметра, как функции координаты 𝐑.
Пусть ℵ – банахово пространство ограниченных функций 𝐳(𝐑) в
области 𝑉 евклидова пространства. Данной областью может являться
геометрический контур месторождения, в рамках которого осуществлялся
прогноз параметра 𝐳 в области его пересечения со скважинами. В таком
случае, функцию принадлежности 𝜇(𝐳(𝐑)) следует считать нелинейным
ограниченным оператором, отражающим ℵ в себя, область значения которого
определена функциями 𝜇𝐳 (𝐑), обладающими смыслом достоверности
значения параметра 𝐳(𝐑) в точке 𝐑 ∈ 𝑉. Каждой пространственной модели
𝐳(𝐑) распределения изучаемого параметра функцией 𝜇𝐳 (𝐑) рассчитывается
покоординатная оценка достоверности значения параметра 𝐳 в этой модели.
В
данном
контексте
функция
принадлежности
𝜇(𝐳(𝐑))
является
информационным полем для оценки достоверности пространственного
распределения параметра 𝐳. Она также называется полем достоверности на
пространстве ℵ. Такое поле служит основой для построения вариантов
моделей
из
соответствующих,
априори,
введенным
принципам
оптимальности. Данное конструирование является дефазификацией функции
принадлежности 𝜇(𝐳(𝐑)).
105
Определение 1. Оптимальной дефазификацией поля 𝜇(𝐳(𝐑)) называется
такая модель 𝐳 𝑑 (𝐑), что для каждого 𝐑 ∈ 𝑉:
𝜇 (𝐳 𝑑 (𝐑)) = max 𝜇(𝐳(𝐑)),
𝐳
где 𝐳 𝑑 (𝐑) – оптимальная дефазифированная модель.
Утверждение
1.
Оптимальная
дефазифицированная
модель
однозначна, если для каждого 𝐑 𝑖 𝜇(𝐳, 𝐑 𝑖 ) сильно выпукла по параметру 𝐳.
Утверждение 2. Поле
дефазификаций
𝜇(𝐳(𝐑))
𝐳 𝑘 (𝐑), 𝑘 = 1 ÷ 𝐾,
допускает
образованных
систему локальных
функциями
из
совпадающими в 𝐑 𝑖 со значениями 𝐳 𝑘 (𝐑 𝑖 ), 𝑘 = 1 ÷ 𝐾, если для одного или
нескольких 𝐑 𝑖 𝜇(𝐳, 𝐑 𝑖 ) имеет систему локальных максимумов 𝐳 𝑘 (𝐑 𝑖 ), 𝑘 =
1 ÷ 𝐾, отличных от 𝐳 𝑑 (𝐑 𝑖 ).
Определение
2.
Локальными
дефазификациями
поля
𝜇(𝐳(𝐑))
называется система 𝐳 𝑘 (𝐑), 𝑘 = 1 ÷ 𝐾.
Процесс построения локальных дефазификаций состоит в:
упорядочении и определении корреляций между 𝐳 𝑘 (𝐑 𝑖 ), 𝑘 = 1 ÷
𝐾, для разных 𝐑 𝑖 ;
интерполяции состоявших во взаимоотношении локальных
экстремумов 𝜇(𝐳, 𝐑 𝑖 );
собственно,
конструировании
локальных
дефазификаций
𝐳 𝑘 (𝐑), 𝑘 = 1 ÷ 𝐾.
Обычно, величина 𝐾 не превышает трех, и, как правило, равна двум.
При условии отсутствия максимумов в 𝜇(𝐳(𝐑)), отличающихся от 𝐳 𝑑 (𝐑 𝑖 ),
дефазификаций нет и 𝐾 = 1. Оптимальной дефазифицированной моделью
является одна из локальных дефазификаций. Таким образом, в реальных
условиях, упорядочение и корреляция локальных экстремумов происходит
для
того,
чтобы
построить
максимум
две
–
три
локальные,
дефазифицированные модели, совпадающих в большей части области 𝑉 с
оптимально дефазифицированной моделью.
106
Локальные дефазификации есть варианты построения моделей, и
составляют основу многовариантности метода оценки достоверности
промыслово-геологической
модели
на
основе
принципов
нечеткого
моделирования.
Локальные
дефазификации
𝐳 𝑘 (𝐑), 𝑘 = 1 ÷ 𝐾
отличаются
своей
функцией принадлежности соответствующей модели 𝐳 𝑘 (𝐑), и являются
реализацией нечеткой функции из ансамбля ℵ. Смысл 𝜇z𝑘 (R) = 𝜇 (𝐳 𝑘 (𝐑))
состоит в том, что данная величина равна достоверности значения параметра
𝐳 𝑘 (𝐑) в каждой точке R [101, 102].
5.3.
Обозначим
Сечения
𝑆 𝛼 (𝜇 (𝐳 𝑘 (𝐑))) = 𝑆 𝛼 (𝜇z𝑘 (R))
построенной
𝛼-сечение
модели 𝐳 𝑘 (𝐑), как система таких значений R, что функция принадлежности
𝜇 (𝐳 𝑘 (𝐑)) > 𝛼:
𝑆 𝛼 (𝜇z𝑘 (R)) = {𝐑: 𝜇 (𝐳 𝑘 (𝐑)) > 𝛼}.
𝑆 𝛼 (𝜇z𝑘 (R))
−
это
подобласть в 𝑉, внутри которой значение поля достоверности 𝜇 (𝐳 𝑘 (𝐑)) не
меньше, чем 𝛼. Сужение 𝐳 𝑘 (𝐑) с 𝑉 на 𝑆 𝛼 (𝜇z𝑘 (R)) обозначим (𝐳 𝑘 )𝛼 (𝐑).
Данное сужение выделяет фрагменты модели, которые соответствуют ее αсечению – имеющие достоверность не ниже 𝛼.
Увеличив параметр 𝛼, 𝑆 𝛼 (𝜇z𝑘 (R)) уменьшится и для 𝛼2 ≥ 𝛼1
𝑆 𝛼1 (𝜇 (𝐳 𝑘 (𝐑))) ⊆ 𝑆 𝛼2 (𝜇 (𝐳 𝑘 (𝐑))). Таким образом, с возрастанием 𝛼,
область
𝑆 𝛼 (𝜇 (𝐳 𝑘 (𝐑)))
уменьшается
и
(𝐳 𝑘 )𝛼 (𝐑)
определяется
с
достоверностью 𝛼 на более узком множестве. При значении 𝛼 близком к
максимально возможному от 𝐳 𝑘 (𝐑) в (𝐳 𝑘 )𝛼 (𝐑), остается только узкая область
определения.
Анализ функции принадлежности 𝜇(𝐳(𝐑)), определяющую нечеткую
функцию 𝐳(𝐑),
состоит
в
построении
совокупности
ее
локальных
107
дефазификаций 𝐳 𝑘 (𝐑), 𝑘 = 1 ÷ 𝐾, а затем системы α-сечений для них
𝑆 𝛼 (𝜇z𝑘 (R)). Значение уровня доверия 𝛼 определяется в диапазоне от нуля до
единицы и определяет величину достоверности не ниже которой определено
соответствующее сечение.
Для демонстрации построения α-сечений рассмотрим распределение
достоверности пористости и нефтенасыщенности тестового месторождения.
Исходными данными для прогноза являлись одновременно измеренные
значения:
пористости
интерпретации
остаточной
результатов
(определяется
геофизических
водонасыщенности
анализа
керна)
по
исследований
пород
и
материалам
скважин),
пористости
и
по
материалам
зависимость
параметров
(определяется
принятая
результатов
нефтенасыщенности и остаточной водонасыщенности.
На рисунках 61 и 62 представлен исходный трехмерный куб
распределения пористости и рассчитанное распределение достоверности
пористости.
а
б
Рисунок 61 – Трехмерный куб распределения пористости (а) и ее
достоверности (б)
Если увеличить 𝛼 величину, то область 𝑆 𝛼 (𝜇z𝑘 (R)) уменьшится и для
𝛼2 ≥ 𝛼1 𝑆 𝛼1 (𝜇 (𝐳 𝑘 (𝐑))) ⊆ 𝑆 𝛼2 (𝜇 (𝐳 𝑘 (𝐑))). При возрастании значения 𝛼,
область 𝑆 𝛼 (𝜇(𝐳(𝐑))) уменьшится в размерах, и 𝐳(𝐑) с достоверностью 𝛼
определится на более узком множестве. Из чего следует, что при значении 𝛼
108
близком к максимальному значению, область 𝑆 𝛼 (𝜇(𝐳(𝐑))) имеет малый
размер. Последовательность α-сечений приведена в таблицах 17 и 18.
Таблица 17. Последовательность α-сечений распределения достоверности для
пористости
𝛼 = 0.1
𝛼 = 0.2
𝛼 = 0.3
𝛼 = 0.4
𝛼 = 0.5
𝛼 = 0.6
𝛼 = 0.7
𝛼 = 0.8
𝛼 = 0.9
Для более наглядного отображения здесь и далее трехмерные кубы
представлены в ином графическом исполнении.
а
б
Рисунок 62 – Трехмерный куб распределения пористости (а) и ее
достоверности (б) (вид сверху)
109
Таблица 18. Последовательность α-сечений распределения достоверности для
пористости (вид сверху)
𝛼 = 0.1
𝛼 = 0.2
𝛼 = 0.3
𝛼 = 0.4
𝛼 = 0.5
𝛼 = 0.6
𝛼 = 0.7
𝛼 = 0.8
𝛼 = 0.9
На рисунках 63 и 64 представлен исходный трехмерный куб
распределения
нефтенасыщенности
достоверности
нефтенасыщенности.
и
рассчитанное
распределение
Последовательность
α-сечений
приведены в таблицах 19 и 20.
110
б
а
Рисунок 63 – Трехмерный куб распределения нефтенасыщенности (а) и ее
достоверности (б)
Таблица 19. Последовательность α-сечений распределения достоверности для
нефтенасыщенности
𝛼 = 0.1
𝛼 = 0.2
𝛼 = 0.3
𝛼 = 0.4
𝛼 = 0.5
𝛼 = 0.6
𝛼 = 0.7
𝛼 = 0.8
𝛼 = 0.9
111
а
б
Рисунок 64 – Трехмерный куб распределения нефтенасыщенности (а) и ее
достоверности (б) (вид сверху)
Таблица 20. Последовательность α-сечений распределения достоверности для
нефтенасыщенности (вид сверху)
𝛼 = 0.1
𝛼 = 0.2
𝛼 = 0.3
𝛼 = 0.4
𝛼 = 0.5
𝛼 = 0.6
𝛼 = 0.7
𝛼 = 0.8
𝛼 = 0.9
112
Расчет
достоверности
пересечение
запасов
распределений
углеводородов
достоверности
вычисляется
пористости
как
и
нефтенасыщенности: min(Кп , Кн ). Результат расчета достоверности запасов
представлен на рисунках 65 и 66, а последовательность его 𝛼-сечений в
таблицах 21 и 22.
Рисунок 65 – Трехмерный куб распределения достоверности запасов
углеводородов
Рисунок 66 – Трехмерный куб распределения достоверности
нефтенасыщенности (вид сверху)
113
Таблица 21. Последовательность α-сечений распределения достоверности для
запасов углеводородов
𝛼 = 0.1
𝛼 = 0.2
𝛼 = 0.3
𝛼 = 0.4
𝛼 = 0.5
𝛼 = 0.6
𝛼 = 0.7
𝛼 = 0.8
Анализ структуры α-сечений дает объективное представление об
информационном обеспечении для построенной модели. Это позволяет:
управлять рисками – оценивать технологические и экономические риски при
введении месторождения в эксплуатацию и принимать соответствующие
управленческие решения; выделять области сравнительной недоразведоности
месторождения и обеспечить оптимальное планирование работ по доразведке
месторождений, тем самым, с одной стороны обеспечивает управление
рисками в принятии решений о технологии разработки месторождении, а с
другой, выделяет объективно не обеспеченные информацией зоны модели,
требующие доразведки, определяя тем самым направление и характер работ
по доразведке месторождения.
114
Таблица 22. Последовательность α-сечений распределения достоверности для
запасов углеводородов (вид сверху)
𝛼 = 0.1
𝛼 = 0.2
𝛼 = 0.3
𝛼 = 0.4
𝛼 = 0.5
𝛼 = 0.6
𝛼 = 0.7
𝛼 = 0.8
Как видно из представленного эксперимента, явно прослеживается
неоднородность
компонент
модели
по
уровню
информационной
обеспеченности принятой модели. Достаточно высокая дифференциация по
уровню
достоверности
может
свидетельствовать
не
только
о
недоразведанности месторождения, но и о недостатках в петрофизическом
обеспечении, которые рельефно проявляются при рассмотрении нечеткой
петрофизической модели.
115
5.4.
Заключение
Метод пространственного продолжения основан на экспоненциальной
модели, результатом которого является распределение достоверности
принятого прогнозного параметра по пропластку. Полученное распределение
достоверности накладывается на физико-геологическую модель.
Информационно-экспертная
оценка
результатов
промыслово-
геологического моделирования в нефтегазовой геологии заключается в
контроле достоверности построенных моделей, анализируя их соответствие
имеющимся исходным петрофизическим данным. Модель реализуется на
основе представления результатов моделирования значений параметров в
кубе в форме, ранжированной по достоверности системы его значений.
Ранжирование следует, и согласовано с неопределенностью: исходной
информации; правил прогнозирования значений параметров.
116
Заключение
Основные результаты выполненных диссертационных исследований
состоят в следующем:
1. Разработан метод конструирования функции принадлежности на
основе
точечных
источников.
Данный
метод
позволяет
исключать
«случайные» данные, при этом сохраняя структуру зависимости исходных
данных.
2. Проведен сравнительный анализ эффективности базовых функций,
используемых при построении функции принадлежности. По результатам
анализа выбрана экспоненциальная модель, при которой обеспечивается
максимальная неопределенность в интерполируемых параметрах.
3. Разработан метод анализа функций принадлежности на основе
построения системы α-сечений для поиска оптимального эффективного
параметра.
4. Доказано, что для прогнозирования фильтрационно-емкостных
параметров по геофизическим параметрам следует использовать нечеткий
логический вывод.
5. Разработана
технология
пространственного
распределения
достоверности по пропластку в физико-геологической модели.
6. Разработан
метод
экспертной
оценки
информационной
обеспеченности на основе конструирования последовательности α-сечений
построенной функции принадлежности на пространстве моделей. Анализ их
структуры
дает
объективное
представление
об
информационном
обеспечении компонент построенной модели.
7. Разработан
программный
комплекс
построения
нечетких
петрофизических моделей с разными базовыми функциями, композиции
нечетких
отношений,
интерполяции
распределения
достоверности
прогнозных значений по пропластку.
117
Перспективами
дальнейшей
разработки
темы
заключаются
в
проведении научно исследовательских работ с ведущими нефтегазовыми
предприятиями
Республики
Коми,
направленных
на
определение
результативности использования разработанных рекомендаций и методов.
118
Список литературы
1.
Шевнин В. А., Рыжов А. А., Делгадо-Родригес О. Оценка
петрофизических параметров грунтов по данным метода сопротивлений.–
М.: Геофизика, 2006.– № 4.– С. 37-43
2.
Добрынин В. М., Вендельштейн Б. Ю., Кожевников Д. А.
Петрофизика (физика горных пород): Учебник для вузов. 2-е изд./Под ред. Д.
А. Кожевникова.– М.: Нефть и газ, 2004.– 368 с.
3.
Киркинская В. Н. и Смехов Е. М. Карбонатные породы-
коллекторы нефти и газа.– Л.: Недра, 1981.– 255 с.
4.
Дахнов В. Н. Геофизические методы определения коллекторских
свойств и нефтегазонасыщения горных пород.– М.: Недра, 1985.– 310 с.
5.
сложных
Рыскаль О. Е., Ахметов Р. Т., Малинин В. Ф. Типы пористости
карбонатных
коллекторов по
результатам
петрофизических
исследований // Геология нефти и газа.– 1987.– №7.– С.47-50.
6.
Головин К. Б., Головин Б. А., Калинникова М. В. Типы и виды
коллекторов:
Учебно-методическое
пособие.–
Саратов:
Саратовский
государственный университет им. Н. Г. Чернышевского, 2014.– 60 с.
7.
Борисов Б. Ф. Результаты анализа верейских песчаников
Покровского месторождения, отобранных с применением бурового раствора
на углеводородной основе/ Б. Ф. Борисов, Н. В. Калери, О. Ю. Худякова,
А. П. Моргунов // Геология нефти и газа.– 1981.– №10.– С.47-50.
8.
Ниязова А. Т., Умирова Г. К., Ахметов Е. М. Характеристика
коллекторских свойств продуктивных горизонтов по керну и гис в
терригенных
отложениях
месторождения
Арыстановское
//
труды
Международных сатпаевских чтений «Роль и место молодых ученых в
реализации новой экономической политики Казахстана».– 2015.– Том III.– С.
39-44
9.
Леонтьев Е. И. Моделирование в петрофизике.– М.: Недра, 1978.–
126 с.
119
10.
Кобрунов А. И. Математические методы моделирования в
прикладной геофизике (избранные главы). В 2-х ч. Ч. 1. Функциональноаналитические основы: учебное пособие / А. И. Кобрунов. – Ухта : УГТУ,
2014. – 224 c.
11.
Алтунин,
А.
Е.
Сравнительный
анализ
использования
вероятностных и нечётких методов оценки неопределённостей и рисков при
подсчёте запасов и ресурсов углеводородов / А. Е. Алтунин, М. В. Семухин //
Нефтяное хозяйство.– 2011. - № 9.– С. 44-49
12.
Zadeh L. A. Fuzzy sets// Information and Control.– 1965.– Vol.8, N
3.– P. 338-353.
13.
Кобрунов А. И., Григорьевых А. В. Методы нечёткого
моделирования
при
изучении
взаимосвязей
между
геофизическими
параметрами // Геофизика.– 2010.– №. 2.– С. 17-23.
14.
Кобрунов А. И., Григорьевых А. В., Художилова А. Н.
Технология фазификации метода нечетких петрофизических композиций //
Материалы
III-ей
международной
научно-практической
конференции
«Нефтегазовая геофизика-нетрадиционные ресурсы».– Ивано-Франковск,
2013. – С. 25-30.
15.
Кобрунов А. И., Григорьевых А. В., Художилова А. Н. Нечеткие
модели в задачах корреляционно-регрессионного анализа // Материалы 37-й
сессии международного семинара им. Д.Г. Успенского.– М.: ИФЗ РАН, 2010.
– С. 185-190.
16.
Кобрунов А. И. Метод нечетких петрофизических композиций
при прогнозировании петрофизических параметров / Кобрунов А. И.,
Кулешов В. Е, Могутов А. С., Художилова А. Н. // Вестник института
Геологии Коми НЦ УРО РАН.– 2011. – № 9. – С. 18–24
17.
Кобрунов А. И. Моделирование неопределенности при прогнозе
параметров геолого-промысловых моделей месторождений углеводородов /
Кобрунов А. И., Кулешов В. Е., Могутов А. С., Дорогобед А. Н. // Нефтяное
хозяйство.– 2014. – № 7. – С.78–80
120
18.
Вендельштейн Б. Ю., Резванов Р. А. Геофизические методы
определения параметров нефтегазовых коллекторов: при подсчете запасов и
проектировании разработки месторождений.– М.: Недра, 1978.– 318 с.
19.
Mamdani E. H. Application of fuzzy algorithms for control of simple
dynamic plant // Electrical Engineers, Proceedings of the IEE.– 1974.– 121(12)– P.
1585-1588.
20.
Кобрунов А. И. Математические основы теории интерпретации
геофизических данных: учебное пособие.– М. : ЦентрЛитНефтеГаз, 2008.–
288 с.
21.
Тиханов А. Н., Арсенин В. Я. Методы решения некорректных
задач.– М.: Наука, 1979.– 285 с.
22.
Добрынин В. М., Вендельштейн Б.Ю., Резванов Р.А., Африкян
А.Н. Промысловая геофизика: Учебник для вузов / Под ред. В. М.
Добрынина.– М.: Недра, 1986.– 342 с.
23.
Еремин Н. А. Моделирование месторождений углеводородов
методами нечеткой логики.– М: Наука, 1994.– 462 c.
24.
Жирабок А. Н. Нечеткие множества и их использование для
принятия решений // Соросовский образовательный журнал. – 2001. – Т. 7, №
2. – С. 109–115.
25.
Park S. Y., Bera A. K. Maximum entropy autoregressive conditional
heteroskedasticity model // Journal of Econometrics. – 2009. – Т. 150, № 2. – С.
219–230.
26.
Сабельников И. С., Потехин Д. В. Геостатистические методы
определения
достоверности
трехмерной
литологической
модели
месторождения [Электронный ресурс] // Website EAGE : Geomodel 2014 –
16th EAGE science and applied research conference on oil and gas geological
explorationand
development
:
аbs.
–
URL:
http://www.earthdoc.org/publication/publication details/?publication=77944.
27.
Кофман А. Введение в теорию нечетких множеств. – М. : Радио и
связь, 1982. – 432 с.
121
28.
Bardossy A., Duckstein L. Fuzzy Rule-Based Modeling with
Applications to Geophysical, Biological and Engineering Systems. – Boca Raton,
FL : CRC Press, 1995. – 256 p.
29.
Paasche H., Tronicke J. Cooperative inversion of 2D geophysical data
sets: A zonal approach based on fuzzy c-means cluster analysis // Geophysics. –
2007. – Vol. 72, № 3. – P. A35-A39.
30.
Алтунин А. Е., Семухин М. В. Модели и алгоритмы принятия
решений в нечетких условиях.– Тюмень : Изд-во ТГУ, 2000.– 352 с.
31.
Алтунин А. Е., Семухин М. В., Ядрышникова О. А. Оптимизация
инвестиционных планов проведения геолого-разведочных работ на основе
нечеткого математического программирования // Нефтяное хозяйство.–
2009.– № 10.– С. 30-32.
32.
Пороскун В. И., Стернин М. Ю., Шепелев Г. И. Вероятностная
оценка запасов на начальных стадиях изучения залежей нефти и газа //
Геология нефти и газа.– 1999.– № 5-6.– С. 59-63.
33.
Старобинец А. Е. Методика вероятностной оценки достоверности
запасов и ресурсов нефти и газа // Геология, геофизика и разработка
нефтяных месторождений.– М. : ВНИИОЭНГ, 2000.– № 4.– С. 14-16.
34.
Старобинец А. Е. Методика оценки неопределенностей в
значениях подсчетных параметров при подсчете ресурсов и запасов
углеводородов на основе вероятностного подхода // Процесс принятия
управленческих решений на основе экономического анализа работ по
поискам и разведке нефти и газа.– М. : ВНИИОЭНГ, 2001.– С. 128-147.
35.
Ампилов
Ю.
П.
От
сейсмической
интерпретации
к
моделированию и оценке месторождений нефти и газа.– М. : Газоил пресс,
2008. – 385 с.
36.
Добрынин В. М. Геофизические исследования скважин : учебник
для вузов / В. М. Добрынин, Б. Ю. Вендельштейн, Р. А. Резванов, А. Н.
Африкян.– М. : Нефть и газ, 2004.– 400 с.
122
37.
Кобрунов А. И. Математические методы моделирования в
прикладной геофизике (избранные главы). В 2 ч. Ч. 2: Системный анализ и
моделирование в условиях неопределенности.– Ухта : УГТУ, 2014.– 155 с.
38.
Ouenes A. Practical application of fuzzy logic and neural networks to
fractured reservoir characterization //Computers & Geosciences. – 2000. – Т. 26,
№. 8. – С. 953-962.
39.
Tani T., Sakoda M., Tanaka K. Fuzzy modeling by ID3 algorithm and
its application to prediction of heater outlet temperature //Fuzzy Systems, 1992.,
IEEE International Conference on. – IEEE, 1992. – С. 923-930.
40.
Кожевников Д. А. Проблемы интерпретации данных ГИС //
Научно-технический
вестник
Ассоциации
исследователей
скважин
«Каротажник».– 1997.– №34.– С. 7-27.
41.
Кобрунов А. И., Кулешов В. Е., Могутов А. С. Адаптация метода
нечётких петрофизических композиций для определения подсчётных
параметров Низевого месторождения // Нефтегазовое дело.– 2011.– №6.– С.
307-315.
42.
Кобрунов А. И., Кулешов В. Е., Могутов А. С. Повышение
достоверности подсчёта запасов углеводородов на основе метода нечётких
петрофизических композиций. SPE – 162038.
43.
Методические рекомендации по подсчету геологических запасов
нефти и газа объемным методом / Под редакцией В. И. Петерсилье, В. И.
Пороскуна,
Г.
Г.
Яценко.–
Москва-Тверь:
ВНИГНИ,
НПЦ
«Тверьгеофизика».– 2003.– 261 с.
44.
Кобрунов А. И., Бурмистрова О. Н. Полугрупповые свойства
распределения Гаусса // Фундаментальные исследования.– 2015.– №2, Ч. 14.–
С. 3044-3047
45.
Агаев Х. Б. Прогнозирование модели скоростей распространения
поперечных волн по данным геофизических исследований скважин и
сейсморазведки с применением нейронных сетей //Геоiнформатика.– 2012.–
№4(44).– С. 46-52
123
46.
Гришин Ф. А. Промышленная оценка месторождений нефти и
газа: Издание 2-е перераб. и дополненное.– М.: Недра, 1985.– 277 с.
47.
Латышова М. Г., Дьяконова Т. Ф., Цирульников В. П.
Достоверность геофизической и геологической информации при подсчёте
запасов нефти и газа.– М.: Недра, 1986.– 121 с.
48.
Вендельштейн
Б.
Ю.
Геофизические
методы
изучения
подсчетных параметров при определении запасов нефти и газа / Б. Ю.
Вендельштейн, Г. М. Золоева, Н. В. Царева и др.– М.: Недра, 1985.– 248 с.
49.
Некорректные задачи математической физики и анализа / М.М.
Лаврентьев и др. – М.: Наука, 1980. – 286с.
50.
Методические
параметров
залежей
исследований
рекомендации
нефти
скважин
с
и
газа
по
по
применением
определению
материалам
результатов
подсчетных
геофизических
анализа
керна,
опробований и испытаний продуктивных пластов // Под редакцией
Вендельштейна Б. Ю., Козяра В. Ф., Яценко Г. Г. – Калинин: НПО
«Союзпромгеофизика», 1990. – 261 с.
51.
Никитин А. А. Теоретические основы обработки геофизической
информации: Учебник для вузов.– М: Недра, 1986.– 342 с.
52.
Матерон Ж. Основы прикладной геостатистики.– М.: Мир, 1968.–
53.
Белонин М. Д. и др. Метод Монте-Карло в нефтяной геологии:
408 с.
Обзор. – М.: ВИЭМС, 1981. – 269 с.
54.
Жабрев И. П. Нечеткая математическая модель при подсчете
запасов / Жабрев И. П., Хургин Я. И. // Геология нефти и газа. – 1993. – №
11. – С. 45-51.
55.
Хургин Я. И. Проблемы неопределенности в задачах нефти и
газа. – Москва-Ижевск: Институт компьютерных исследований, 2004. – 320 с.
56.
Бочарников В. П., Цыганок А. В. Нечеткая технология анализа и
управления рисками //Зб. наук. праць за матеріалами І Всеукр. наук.-практ.
124
конф.«Проблеми економічного ризику: аналіз та управління»(26-28 жовтня
1998 р.).– К.: КНЕУ, 1998. – С. 8-9.
57.
Finol J., Guo Y. K., Jing X. D. A rule based fuzzy model for the
prediction of petrophysical rock parameters //Journal of Petroleum Science and
Engineering.– 2001.– Т. 29, №. 2.– С. 97-113.
58.
Алексеев А. В. Применение нечеткой математики в задачах
принятия решений // Методы и системы принятия решений. – Рига: РПИ,
1983. – C.38-42.
59.
Алексеев А.
В.
Проблемы
разработки
математического
обеспечения выполнения нечетких алгоритмов // Модели выбора альтернатив
в нечеткой среде.– Рига, 1984.– C. 79-82.
60.
Беллман Р. Принятие решений в расплывчатых условиях /
Беллман Р., Заде Л. // Вопросы анализа и процедуры принятия решений.– М.:
Мир, 1976.– С.172-215.
61.
Беллман Р.
Прикладные
задачи
динамического
программирования / Беллман Р., Дрейфус С.– М: Наука, 1965.– 460 с.
62.
Дюбуа Д. К анализу и синтезу нечетких отображений / Дюбуа Д.,
Прад А. // Нечеткие множества и теория возможностей.– 1986.– C. 229-240.
63.
Кирсанов А. Н. Подсчет запасов газа с использованием теории
нечетких множеств. Петрофизическое обеспечение подсчета запасов нефти и
газа / Кирсанов А. Н., Семухин М. В., Адлер В. В. // Труды института
ЗапСибНИГНИ. – Тюмень: ЗапСибНИГНИ, 1989. – C. 86-107.
64.
Леоненков А. В. Нечеткое моделирование в среде MATLAB и
fuzzyTECH / Леоненков А. В.– СПб.: БХВ-Петербург, 2003.– 736 с.
65.
Орлов А. И. Теория принятия решений: учебное пособие.– М.:
Март, 2004.– 656 с.
66.
Mamdani E. H. Twenty years of fuzzy control: experiences gained and
lessons learned, IEEE Internal // Conf on Fuzzy Systems.– 1993.– P. 339-344.
67.
Zadeh L. A. Fuzzy sets as a basis for a theory of possibility // Fuzzy
sets and systems.– 1978.– №1.– P. 3-28.
125
68.
Zadeh L. A theory of approximate reasoning (AR) //Machine
Intelligence.– 1979.– V. 9.– P. 149-194
69.
Zadeh L. A. Fuzzy Logic, Neural Network and Soft Computing //
Communication of the ACM.– 1994.– V. 37, № 3.– P. 77-84.
70.
Zadeh L. A. Fuzzy Logic = Computing With Words // IEEE
Transactions of Fuzzy Systems.– 1996.– V.4.– P. 103-111.
71.
Zadeh L. A. Toward a perception-based theory of probabilistic
reasoning with imprecise probabilities // Journal of Statistical Planning and
Inference.– 2002.– V.105.– P. 233-264.
72.
Семухин М. В. Нечеткие оценки запасов нефти // Сборник
докладов международной конференции по мягким вычислениям. – С-Пб:
СПбГЭТУ «ЛЭТИ», 2003.– Т.2.– С. 164-167.
73.
Семухин М. В. Теория нечетких множеств: учебно-методическое
пособие / Семухин М. В.– Тюмень: ТюмГУ, 1999.– 50 с.
74.
Ампилов Ю. П. Методы геолого-экономического моделирования
ресурсов и запасов нефти и газа с учетом неопределенности и риска /
Ампилов Ю. П.– М.: Геоинформмарк, 2002.– 201с.
75.
Борисов В. В. Нечеткие модели и сети / Борисов В. В.,
Круглов В. В., Федулов А. С..– М.: Горячая линия - Телеком, 2007.– 284 с.
76.
Иванов М. К. Петрофизические методы исследования кернового
материала (Терригенные отложения): учебное пособие в 2-х книгах. Кн. 1. /
Иванов М. К., Бурлин Ю. К., Калмыков Г. А., Карнюшина Е. Е., Коробова Н.
И.– М.: Изд-во Моск. ун-та, 2008.– 112 с.
77.
Кандель А. Нечеткие множества, нечеткая алгебра, нечеткая
статистика / Кандель А., Байатт У. Дж. // Труды американского общества
инженеров-радиоэлектроников.– 1978.– Т. 66, № 12.– С. 37-61.
78.
Заде
Л.
А.
Размытые
множества
и
их
применение
в
распознавании образов и кластер-анализе // Классификация и кластер.–
1980.– С. 208-247.
126
79.
Ягер P. P. Множества уровня для оценки принадлежности
нечетких подмножеств // Нечеткие множества и теория возможностей.
Последние достижения: пер. с англ./под ред. Р. Р. Ягера.– М.: Радио и связь,
1986.– С.71-78.
80.
Никитин А. А. Теоретические основы обработки геофизической
информации / Никитин А. А, Петров А. В.– М: РГТУ, 2008.– 112 с.
81.
Рац М. В. Неоднородность горных пород и их физических
свойств.– М.: Наука, 1968.– 110 c.
82.
Рац М. В. Трещиноватость и свойства трещиноватых горных
пород / Рац М. В., Чернышев С. Н. – М.: Недра, 1970. – 164 c.
83.
Жданов М. А. Нефтегазопромысловая геология и подсчёт запасов
нефти и газа : учеб. пособие для вузов / 2-е изд., перераб. и доп. – М. : Недра,
1981. – 453 с.
84.
Алтунин Е. А. Нечеткие методы идентификации и управления
процессами нефтегазодобычи: дис. канд. техн. наук. : 05.13.01 / Алтунин
Евгений Александрович. – Тюмень: ТГНГУ, 2002. – 203 с.
85.
Алтунин А. Е. Методы определения функций принадлежности в
теории размытых множеств / Алтунин А. Е., Востров Н. Н. // Труды
ЗапсибНИГНИ.– Тюмень: ЗапсибНИГНИ, 1980. – Вып.154. – С. 62-72.
86.
Алтунин А. Е., Семухин М. В. Вероятностный и нечеткий
подходы к оценке запасов // Сборник трудов международного научнотехнического семинара «Информационные системы и технологии в геологии
и нефтегазодобыче». – Тюмень: Вектор Бук, 2005. – C. 143-145.
87.
Алтунин А. Е. , Семухин М. В. Расчеты в условиях риска и
неопределенности в нефтегазовых технологиях: Монография.– Тюмень: Издво Тюмен. гос. ун-та, 2005. – 220 с.
88.
Алтунин А. Е. Вероятностные и нечеткие оценки запасов нефти /
Алтунин А. Е., Семухин М. В., Алтунин Е. А., Ядрышникова О. А. //
Автоматизация, телемеханизация и связь в нефтяной промышленности. –
2003. – №9. – C.27-31.
127
89.
Алтунин А. Е. Вероятностные и нечеткие оценки запасов нефти /
Алтунин А. Е., Семухин М. В., Алтунин Е. А., Ядрышникова О. А. //
Нефтепромысловое дело. – 2003. – № 10. – C. 54-58.
90.
Алтунин А. Е. Применение теории нечеткости для оценивания
технологических параметров в АСУ ВПО «Тюменгазпром» / Алтунин А. Е.,
Семухин М. В., Чуклеев С. Н. // Проблемы нефти и газа Тюмени : Труды
ЗапсибНИГНИ. – Тюмень: ЗапсибНИГНИ, 1983. – Вып. 58. – С. 57-59.
91.
Ампилов Ю. П. Сейсмическая интерпретация: опыт и проблемы.–
М.: Геоинформмарк, 2004.– 277 c.
92.
Ампилов Ю. П. и др. Сопоставление альтернативных методов
прогноза фильтрационно-емкостных свойств по данным сейсморазведки
//Geomodel 2008-10th EAGE science and applied research conference on oil and
gas geological exploration and development. – 2008
93.
Лаврик С. А., Логинов Д. В. Оценка значимости сейсмических
атрибутов на основе многослойных нейронных сетей //Геонауки – от новых
идей к новым открытиям: труды 3-ей международной конференции и
выставки.– СПб, 2008.– 106 с.
94.
для
Лаврик С. А., Логинов Д. В. Применение статистических методов
определения
информативного
набора
сейсмических
атрибутов
//Geomodel 2008-10th EAGE science and applied research conference on oil and
gas geological exploration and development. – 2008.
95.
Mendrii I., Dovbnich M. Seismic Attribute Analysis for Mapping of
Fracture Zones in Coal-bearing Formations //76th EAGE Conference and
Exhibition 2014. – 2014.
96.
Лаврик С. А. Результаты анализа эффективности и применимости
статистических
методов
для
определения
информативного
набора
сейсмических атрибутов //Seismic Technology.– 2009.– Т. 6, №. 1.– С. 36-44.
97.
Kobrunov A., Priezzhev I. Hybrid combination genetic algorithm and
controlled gradient method to train a neural network //Geophysics.– 2016.– Т. 81,
№. 4.– С. IM35-IM43.
128
98.
Ампилов Ю. П. и др. Почти все о сейсмической инверсии
//Технологии сейсморазведки. – 2009. – №. 04. – С. 3-16.
99.
Ампилов Ю. П., Лапо А. В. Анализ геолого-экономических
показателей, применяемых при оценке эффективности разведки и освоения
участков недр //Минеральные ресурсы России. Экономика и управление. –
2010. – №. 5. – С. 29-34.
100. Ампилов Ю. П. Многофакторная система оценки месторождений
углеводородов //Минеральные ресурсы России. Экономика и управление. –
2015. – №. 4. – С. 35-43.
101. Кобрунов А.И., Бурмистрова О.Н. Нечеткие геологические
модели и их дефазификации // Успехи современного естествознания. – 2016.
– № 1. – С. 116-119
102. Кобрунов А. И. Доопределение обратных задач гравиметрии
нечеткой информацией // Успехи современной науки. – 2016. – №11-6. – С.
13-18
103. Котяхов Ф. И. Физика нефтяных и газовых коллекторов.– М.:
Недра, 1977.– 287 с.
104. Макаревич В. Н., Искрицкая Н. И., Богословский С. А. Ресурсный
потенциал тяжелых нефтей Российской Федерации: перспективы освоения
//Нефтегазовая геология. Теория и практика. – 2010. – Т. 5, №. 2. – С. 1-13.
105. Головинский В. И. и др. Сейсмофации и нефтегазоносность
верхнепалеогеновых и неоген-четвертичных отложений российской части
Черного моря //Геология, геофизика и разработка нефтяных и газовых
месторождений. – 2002. – №. 2. – С. 35-42.
106. Кожевников Д. А., Коваленко К. В., Арсибеков А. А.
Инвариантность петрофизических связей в адаптивной интерпретации
данных ГИС //Каротажник. – 2009. – №. 7. – С. 69-86.
107. Алиев М. М., Батанова Г. П., Хачатрян Р. О. Девонские
отложения Волго-Уральской нефтегазоносной провинции.– М.: Недра, 1978.–
216 с.
129
108. Ромм Е. С. Структурные модели порового пространства горных
пород. – Л. : Недра. Ленингр. отд-ние, 1985. – С. 157.
109. Тиаб Д. Петрофизика: теория и практика изучения коллекторских
свойств горных пород и движения пластовых флюидов / Перевод с
английского.– М.: ООО «Премиум инжиниринг», 2009.– 868 с.
110. Нестеров
И.
И.
Нефтегазоносность
битуминозных
глин
баженовской свиты Западной Сибири //Советская геология. – 1980. – Т. 64,
№. 2. – С. 115.
111. Добрынин В. М., Мартынов В. Г. Коллектор нефти в
нефтематеринских глинистых толщах // Геология нефти и газа.– 1979.– №7.–
С. 36-43.
112. Конторович А. Э. Геология нефти и газа Западной Сибири /
Конторович А. Э., Нестеров И. И., Салмонов Ф. К. и др.– М.: Недра, 1975
113. Нестеров И. И. Новый тип коллектора нефти и газа //Геология
нефти и газа. – 1979. – №. 10. – С. 26-29.
114. Нестеров И. И. Новые открытия месторождений нефти и газа в
Западной Сибири //Нефть и газ Тюмени. – 1969. – №. 3. – С. 1.
115. Климчук А. Б. Основные особенности и проблемы гидрогеологии
карста: спелеогенетический подход. Сообщ. 1: Пористость и проницаемость
карстовых коллекторов //Пещеры. – 2001. – №. 27. – С. 28-51.
116. Чернышев С. Н. Трещины горных пород.– М.: Наука, 1983.–
125 с.
117. Дахнов
В.
Н.
Интерпретация
результатов
геофизических
исследований разрезов скважин.– М.: Недра, 1982.– 448 с.
118. Голф-Рахт
Т.
Д. Основы нефтепромысловой
геологии и
разработки трещиноватых коллекторов. – М.: Недра, 1986.– 608 с.
119. Бурханов Р. Н., Ханнанов М. Т., Фаррахов И. М. Особенности
оптических свойств нефти горизонтальных скважин //Геология, геофизика и
разработка нефтяных и газовых месторождений. – 2013. – №. 4. – С. 50-54.
130
120. Баюк И. О. и др. Математическое моделирование анизотропных
эффективных упругих свойств карбонатных коллекторов сложного строения
//Технологии сейсморазведки. – 2012. – № 03. – С. 42-55.
121. Авчян Г. М., Матвеенко А. А., Стефанкевич З. Б. Петрофизика
осадочных пород в глубинных условиях.– М.: Недра, 1979. – 224 с.
122. Бека К., Высоцкий И. Геология нефти и газа. – М.: Недра, 1976.–
592 с.
123. Инженер-геолог В., ЦНИГРИ В. Физические свойства горных
пород и полезных ископаемых (петрофизика). Справочник геофизика
//Труды. – 2017.
124. Горбачев Ю. И. Геофизическое исследование скважин. – М.:
Недра, 1990.– 398 с.
125. Абасов М. Т. и др. О влиянии пластового давления на изменение
фильтрационно-емкостных
свойств
терригенных
пород-коллекторов
в
процессе разработки месторождений нефти и газа //Геология нефти и газа. –
1997. – Т. 5. – С. 34-39.
126. Меркулов
В.
П.,
Краснощёкова
Л.
А.
Исследование
пространственной литолого-петрофизической неоднородности продуктивных
коллекторов
месторождений
нефти
и
газа
//Известия
Томского
политехнического университета. – 2002. – Т. 305, №. 6. – С. 296-303.
127. Золоева Г. М. и др. Изучение карбонатных коллекторов методами
промысловой геофизики.– М.: Недра, 1977.– 150 с.
128. Коршак А. А., Шаммазов А. М. Основы нефтегазового дела.
Учебник для ВУЗов.– Уфа.: ООО «ДизайнПолиграфСервис, 2001. – 544 с.
129. Абасов М. Т. и др. О методе оценки изменения относительной
фазовой
проницаемости
пород-коллекторов
в
процессе
разработки
//Геология, геофизика и разработка нефтяных и газовых месторождений. –
2010. – №. 4. – С. 54-57.
131
130. Дахнов В. Н., Долина Л. П. Геофизические методы изучения
нефтегазоносных коллекторов. – М.: Государственное научно-техническое
издательство нефтяной и горно-топливной литературы, 1959.– 268 с.
131. Баюк И. О. Теоретические основы определения эффективных
физических свойств коллекторов углеводородов //Акустика неоднородных
сред. Ежегодник РАО. – 2011. – №. 12. – С. 107-120
132. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В. Нечеткие
модели прогнозирования физико-геологических и продуктивных параметров
неоднородных
геологических
сред
//
Материалы
42-й
сессии
Международного семинара им. Д. Г. Успенского. – Пермь : Горный ин-т УрО
РАН, 2015. – С. 101-105
133. Кобрунов А. И.,
Кожевникова П. В. Выбор и обоснование
функции принадлежности при прогнозировании параметров геологических
сред в условиях неопределенности // Материалы международного семинара
«Рассохинские чтения». – Ухта : УГТУ, 2015. – ч. 2. – С. 153-163
134. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В.
Многовариантное
моделирование
на
основе
технологий
нечеткого
логического вывода в задачах нефтегазовой отрасли // Материалы III Школыконференции «Гординские чтения». – Москва : ИФЗ РАН, 2015. – С. 102-106
135. Дорогобед А. Н., Кожевникова П. В. Нечеткое отношение
параметров физико-геологической модели и техника их прогноза по
наблюдаемым данным // Международный Студенческий Научный вестник.–
2015. – № 3, ч. 3. – С. 395-397
136. Кобрунов А. И., Кожевникова П. В. Теоретические основы при
прогнозировании
параметров
геологических
сред
в
условиях
неопределенности // Фундаментальные исследования. – 2015. – № 5, ч. 3. –
С. 506-510
137. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В.
Аппроксимация трехмерных связей на основе технологий многовариантного
132
моделирования // Фундаментальные исследования. – 2015. – № 10, ч. 2. – С.
268-273
138. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В.
Использование технологий многовариантного моделирования на основе
теории нечётких множеств для аппроксимации многомерных связей (на
примере трёхмерных) // Ресурсы Европейского Севера. Технологии и
экономика освоения, ISSN 2412-9976. – 2015. – № 2. – С. 23-30.
139. Кобрунов А. И., Бурмистрова О. Н., Кожевникова П. В. Нечеткие
подстановки и принцип Мамдани // Успехи современного естествознания. –
2016. – №1. – С. 96-101
140. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В. Решение
задач нефтегазовой отрасли в условиях неопределенности на основе
нечеткого
моделирования
//
Материалы
Международного
семинара
«Рассохинские чтения».– Ухта: УГТУ, 2016.– №2.– С. 123-128
141. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В. О проблеме
выбора вида функций принадлежности при нечетком моделировании
петрофизических зависимостей // Материалы 43-й сессии Международного
семинара им. Д.Г. Успенского.– Воронеж : Научная книга, 2016.– С. 102-106
142. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В. Результаты
прогноза на основе нечеткого моделирования // Материалы 43-й сессии
Международного семинара им. Д.Г. Успенского.– Воронеж : Научная книга,
2016.– С. 99-101
143. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В. Метод
нечеткого
логического
вывода
и
информационная
обеспеченность
результатов моделирования в нефтегазовой геологии // Геоинформатика.–
2016.– №2.– С. 35-40
144. Кобрунов А. И., Бурмистрова О. Н., Кожевникова П. В.
Увеличение диффузии при нечеткой свертке нечетких отношений //
материалы 44-й сессии Международного семинара им. Д. Г. Успенского. –
Москва: ИФЗ РАН, 2017.– С. 163-165
133
145. Кобрунов А. И., Бурмистрова О. Н., Кожевникова П. В.
Экспериментальная
проверка
коммутативности
композиции
нечетких
отношений с использованием замены переменных // Успехи современной
науки.– 2016.– № 11-6.– С. 24-28
146. Кобрунов А. И., Кожевникова П. В., Ломинский Д. О. Метод
нечеткого логического вывода для прогноза емкостных параметров // Успехи
современной науки.– 2016.– №12-5.– С. 18-24
147. Кобрунов А. И., Бурмистрова О. Н., Кожевникова П. В.
Возрастание неопределенности при нечеткой свертке нечетких отношений //
Успехи современной науки и образования. – 2016. – № 11-7. –С. 93-95
148. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В. Анализ αсечений функций принадлежности для нечетких петрофизических моделей //
Материалы Международного семинара «Рассохинские чтения».– Ухта:
УГТУ, 2017.– №2.– С. 111-117
149. Кобрунов А. И., Дорогобед А. Н., Кожевникова П. В. Элементы
информационной экспертизы результатов геологического моделирования в
нефтегазовой геологии // Геофизика.– 2017.– № 1.– С. 16-21
150. Кожевникова П.В. Методы построения петрофизических моделей
// Геология и полезные ископаемые Западного Урала.– 2017.– № 17.– С. 152155
151. Кожевникова П. В. Моделирование связей между физикогеологическими параметрами // Успехи современной науки и образования. –
2017. – № 8-2. –С. 83-87
134
Отзывы:
Авторизуйтесь, чтобы оставить отзыв