Научный журнал
Успехи современного естествознания
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

СПОСОБЫ ПОСТРОЕНИЯ ФУНКЦИЙ ПРИНАДЛЕЖНОСТИ НЕЧЕТКИХ ВЕЛИЧИН И НЕЧЕТКИХ ОТНОШЕНИЙ ГЕОЛОГО-ГЕОФИЗИЧЕСКИХ, ПРОМЫСЛОВЫХ ПАРАМЕТРОВ

Кобрунов А.И. 1 Мотрюк Е.Н. 1
1 ФГБОУ ВО «Ухтинский государственный технический университет»
В статье рассматривается применение методов нечеткого моделирования для оценки неопределенности геолого-геофизических, промысловых параметров. Необходмость применения данного аппарата связана с тем, что различные фрагменты прогнозируемой модели месторождения определяются с различной степенью надежности. Неучет этих обстоятельств повышает технико-экономические риски принятия решений о проекте разработки залежей и требует объективной оценки меры информационной обеспеченности выполненных построений. В связи с этим предлагается строить нечеткие петрофизические модели с оценкой достоверности компонент прогнозируемой модели. Построение модели включает несколько этапов, один из которых – дефазификация данных и отношений, заключающаяся в том, что с имеющимися данными сопоставляется численное значение меры доверия, выраженное в функции принадлежности. Разработаны способы построения функций принадлежности нечетких величин и нечетких отношений изучаемых параметров. Построение функции принадлежности происходит на основе карт плотностей, в которых учитывается неоднородность имеющейся информации. На основе теоретических и практических исследований наилучшим образом определяющей меру доверия является экспоненциальная функция принадлежности, выведенная из уравнения диффузии в однородной бесконечной среде с заданными начальными условиями с использованием метода Грина. В статье описан также способ выбора эффективного параметра для этих функций. При его выборе следует ориентироваться на то, чтобы поле рассеяния для моделируемой функции принадлежности покрывало исходный полигон данных на 90 %. Приведен пример построения функции принадлежности для зависимости пористость по ГИС – пористость по керну и данных пористости по ГИС.
нечеткие модели
моделирование
функция принадлежности
фазификация
прогноз неопределенности
геолого-геофизические параметры
нечеткие величины
нечеткие отношения
1. 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.
2. Заде Л.А. Нечеткие множества // Нечеткие системы и мягкие вычисления. 2015. Т. 10. Вып. 1. С. 7–22.
Zade L.A. Fuzzy Sets // Fuzzy Systems and Soft Computing. 2015. V. 10. Is. 1. P. 7–22 (in Russian).
3. Штовба С.Д. Введение в теорию нечетких множеств и нечеткую логику. [Электронный ресурс]. URL: http://matlab.exponenta.ru/fuzzylogic/book1/index.php (дата обращения: 28.02.2019).
Shtovba S.D. Introduction to the theory of fuzzy sets and fuzzy logic. [Electronic resource]. URL: http://matlab.exponenta.ru/fuzzylogic/book1/index.php (date of access: 28.02.2019) (in Russian).
4. Алтунин А.Е., Семухин М.В., Ядрышникова О.А. Вероятностные и нечеткие оценки запасов нефти // Нефтепромысловое дело. 2003. № 10. C. 54–58.
Altunin A.E., Semukhin M.V., Yadryshnikova O.A. Probabilistic and fuzzy estimates of oil reserves // Oilfield business. 2003. № 10. Р. 54–58 (in Russian).
5. Чернов В.Г. Основы теории нечетких множеств: учеб. пособие. Владимир: ВГУ, 2010. 96 с.
Chernov V.G. Fundamentals of the theory of fuzzy sets: studies. manual. Vladimir: VGU, 2010. 96 p. (in Russian).
6. Кобрунов А.И., Дорогобед А.Н., Кожевникова П.В. Элементы информационной экспертизы результатов геологического моделирования в нефтегазовой геологии // Геофизика. 2017. № 1. С. 16–21.
Kobrunov A.I., Dorogobed A.N., Kozhevnikova P.V. Elements of information expertise of geological modeling results in oil and gas geology // Geofizika. 2017. № 1. P. 16–21 (in Russian).
7. Кожевникова П.В. Моделирование связей между физико-геологическими параметрами // Успехи современной науки. 2017. № 8–2. С. 83–87.
Kozhevnikova P.V. Modeling the links between physical and geological parameters // Advances in modern science. 2017. № 8–2. P. 83–87 (in Russian).

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

Цель исследования: разработка способов построения функций принадлежности нечетких величин и нечетких отношений геолого-геофизических, промысловых параметров с использованием методов нечеткого моделирования, основы которых изложены в [1–3]. Это позволит оценивать меру неопределенности параметров при проведении прогноза и учитывать полученные результаты при разработке схем проектирования месторождений.

Материалы и методы исследования

Материалом для исследования являются геолого-геофизические, промысловые параметры, а методами исследования – аппарат теории нечетких множеств [4, 5].

Результаты исследования и их обсуждение

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

Фазификация отношений состоит в выражении зависимостей между нечеткими параметрами – переменными, в форме нечетких отношений. Исходными данными формирования функции принадлежности kobr01.wmf для отношения между нечеткими величинами kobr02.wmf служат результаты измеренных пар значений kobr03.wmf для всех из рассматриваемых отношений. Для определенности выберем задачу установления отношения между нечеткими величинами x и y на основе эксперимента kobr04.wmf приведшего к результатам kobr05.wmf, где число измерений равно n. Рассмотрения автоматически переносятся на любые другие пары. Пусть, например, x – экспериментально измеренная пористость по данным ГИС, y – пористость того же образца по керну. Полигон данных kobr06.wmf, полученных в результате измерений kobr07.wmf, приведен на рис. 1, а. Картой плотности данных, или полем рассеяния для связи между параметрами называется сетка, значений плотности данных в каждой ячейке вычисляемых по правилу

kobr08.wmf (1)

где Mij – количество значений, попавших в ячейку, а n – количество всех измеренных значений (рис. 1, б).

kobr1a.tif

kobr1b.tif

kobr1d.tif

в)

Рис. 1. Полигон данных kobr09.wmf: а) зависимость «пористость по ГИС – пористость по керну» и соответствующая ей карта плотности; б) поле рассеяния для связи между параметрами; в) легенда по параметру количества данных, попавших в ячейку, г) легенда по максимальным значениям по оси Y на выбранном интервале по оси Х

а)

б)

kobr1c.tif

г)

kobr2a.tif

а)

kobr2b.tif

б)

Рис. 2. Сопоставление карты плотности (а) и карты источников плотности для поля рассеяния данных (б)

Выбор аналитического вида функции принадлежности одновременно означает формирование гипотезы о причинах рассеяния экспериментальных данных в облако, для которого, конечно, можно вычислить и центр, и разброс, но оно само по себе несет еще скрытую информацию о характере неоднородности среды, приведшей к рассеянию экспериментальных данных [6, 7]. Экспоненциальная модель для нахождения аналитического вида функции принадлежности определяется из уравнения диффузии в однородной бесконечной среде с заданными начальными условиями с использованием метода Грина. Среда состоит из совокупности частей, различающихся между собой одновременными наборами параметров kobr10.wmf. Плотность концентрации элементов в этой смеси характеризуется функцией концентрации μ(s). Искомая функция принадлежности kobr11.wmf – это итог диффузии μ(s), происходящей с некоторым коэффициентом а, и длящейся некоторое время τ kobr12.wmf. Решение принимает вид

kobr13.wmf (2)

ζ – эффективный параметр рассеяния для моделируемой функции принадлежности kobr14.wmf; kobr15.wmf – стартовое распределение источников информации, которому приписывается стартовый параметр рассеяния ζ0, относительно которого далее происходит выбор эффективного параметра рассеяния. Практический выбор kobr16.wmf, определен плотностью источников поля рассеяния.

Определим значения kobr17.wmf по введенной сетке kobr18.wmf. Функция kobr19.wmf суть значения в ячейках ΔSi нормированной функции рассеяния kobr20.wmf. Таким образом, правая часть в (2) – это набор дискретных значений по ячейкам ΔSi, так же как и искомая kobr22.wmf. Принцип компактификации источников формулируется как выбор такого распределения источников kobr23.wmf и ей соответствующей функции принадлежности с kobr24.wmf, которые характеризуются минимальным числом K носителей (ненулевых значений kobr25.wmf в покрытии ΔSk), выбранной сеткой фазового пространства S, при условии приближенного выполнения соотношения (2). Это соответствует задаче целочисленного программирования:

kobr26.wmf (3)

Ее решение выполняется численно, перебором значений K от максимального, определенного построенными значениями kobr27.wmf, в ячейках сетки kobr28.wmf, до минимального, при котором условие неравенства еще выполняется. Результат решения kobr29.wmf называется картой плотности источников данных или полем источников данных. Алгоритм построения карты источников, плотности данных называется алгоритмом сжатия информации на основе уравнения (2). Исходя из карты плотности, источников поля рассеяния, могут быть получены параметризованные эффективным параметром рассеяния ζ модели функции принадлежности для характеристики нечеткого отношения.

Для приведенного на рис. 1, а, примера полигона данных и ему соответствующей карты плотности данных (поля рассеяния), его сопоставление с соответствующей картой плотности источников поля рассеяния приведено на рис. 2.

Переход к плотности источников поля рассеяния позволяет существенно снизить размерность величин.

Основываясь на полученной карте плотности источников kobr30.wmf, по формуле

kobr31.wmf (4)

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

kobr35.wmf (5)

При выборе ζ следует ориентироваться на то, чтобы поле рассеяния для моделируемой функции принадлежности покрывает исходный полигон данных на 90 %. Мера покрытия полем рассеяния исходного полигона служит критерием для подбора эффективного параметра рассеяния. Многочисленные эксперименты показали, что оптимально в качестве функции принадлежности брать именно экспоненциальный вид [6, 7].

Фазификация входных данных состоит в построении функций принадлежности kobr36.wmf для исходного параметра – аргумента x в локальных интервалах Ri, в которых будет выполняться прогнозирование на основе нечеткой петрофизической модели. Если измерена серия значений kobr37.wmf параметров, в одной и той же пространственной точке, то функция принадлежности kobr38.wmf для x в точке Ri примет вид

kobr39.wmf (6)

где σi – оцененная погрешность измерений параметра x в точках пространства Ri. Схема фазификации соответствует объективному представлению данных ГИС в форме функций принадлежности, отражающих структуру и достоверность информации (рис. 3).

kobr3.tif

Рис. 3. Преобразование прогнозных данных по ГИС в функции принадлежности

Пространственное положение единицы прогноза – пропластка, либо залежи, для которой приняты построенные функции принадлежности отношений и функции принадлежности параметра x, служащего основанием для прогноза, контролируется пространственной переменной Ri. Она не участвует в вычислительных формулах, проявляясь лишь в оценке погрешности измеренных параметров, отнесённых именно к этой точке. Пространственное положение прогнозных параметров далее используется при интерполяции моделей.

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

Заключение

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


Библиографическая ссылка

Кобрунов А.И., Мотрюк Е.Н. СПОСОБЫ ПОСТРОЕНИЯ ФУНКЦИЙ ПРИНАДЛЕЖНОСТИ НЕЧЕТКИХ ВЕЛИЧИН И НЕЧЕТКИХ ОТНОШЕНИЙ ГЕОЛОГО-ГЕОФИЗИЧЕСКИХ, ПРОМЫСЛОВЫХ ПАРАМЕТРОВ // Успехи современного естествознания. – 2019. – № 4. – С. 57-62;
URL: https://natural-sciences.ru/ru/article/view?id=37108 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674