Июл 13
Добрый день всем!
На форуме прочитал, что можно сделать анализ неопределенности (Монте-карло) в программном продукте Петрел. Вопрос: как это сделать , если можно алгоритм действии?
Есть ли другие софты по которым можно сделать это? (софты которые можно на просторе инета найти)
Опубликовано
03 Июл 2013
Активность
20
ответов
9252
просмотра
8
участников
2
Рейтинг
Контекст
Petrel x248
В процессах Utilities -> Uncertainty and optimization
далее по мануалу
хех, по подробнее можете
подробно кнопки в хелпе хорошо расписаны, тут главное концепцию понять
а концепция состоит в следующем:
1. нужно создать воркфлоу, которое умеет считать модель от начала до конца (подсчета запасов или расчета гидродинамики). самый "тяжелый" случай - построение скоростной модели, конвертация данных из времени в глубину, построение геологической модели, построение и расчет гидродинамической модели.
2. В процессе uncertainty and optimization в настройках процессов нужные параметры заменить на переменные. Например, в процессе моделирования фаций песчанистость заменить на переменную, в процессе создания контактов - контакт заменить на переменную. в процессе создания фазовых проницаемостей вместо значений концевых точек поставить переменные, в процессе создания аквифера параметры заменить на переменные, даже в редакторе ключевых слов еклипс параметры можно заменить на переменные. Дополнительно это дает возможность посмотреть - как от варьирования переменных меняется расчет настраиваемой на историю модели.
3. названия переменных должно начинаться с $ например $NTG
4. На закладке с переменными (процессе uncertainty and optimization) нужно будет указать как варьировать каждую переменную.
5. в результате расчета петрел насчитает кучу реализаций модели (подсчет запасов или гидродинамику)
Варианты доступных задач в процессе uncertainty and optimization - оценка чувствительности, неопределенности, оптимизация (максимизация добычи, продление полки добычи, можно свою целевую функцию задать), оптимизация с учетом неопределенностей.
Можно и каменным топором (Ёкселем) все сделать... Только формулы преобразования распределения от равномерного к нужному вам задайте - и усе
вот такой ответ ожидал! Спасибо за инфу, нашел книгу от шлюмов " workflow editor and uncertainty analysis courses". Тока на английском и сздесь много чего описано. по иду по этой книге.
Вопрос: объязательно нужно сделать подсчет запасов чтобы сделать анализ или же можно обойтись тока подсчетными параметрами (poro, sw)???
просто у нас геологи старой закалки и они как сказать не доверяют программе, но есть еще иностранцы для которых этот анализ неопределнности ВАЖЕН.
на счет экселя слыхал, вот тока никогда не видел. Есть пример, или сам может когда то сделал такой анализ???
У нас нескоко месторождении и для мальенких структур думается сделать ручками чем в пети морочится, хотя если написать один раз workflow все по летит.
Для начала вам нужно понять в голове что у вас в модели обладает большей неопределенностью и какие параметры действительно окажут вляиняие. Нет смысла задавать распределение всем параметрам без исключения, вы же не хотите же кашу получить. Выберите для начала 2-3 параметра и посмотрите как меняются объемы и какой параметр более чувстителен к измениям.
Насчет запасов - да, обязательно.
про эксель видимо имелось ввиду crystal ball.
Это штука платная. Делает практически тоже самое что и петрель, тока быстрее тк оперирует не кубами а средними значениями.
Самый топорный способ это просто задать три значения для одного-двух подсчетных параметров и посмотреть как меняются запасы просто перемножив каждый вариант.
Ок, ясно!
еще можно в бесплатном екселе без кристал бола. желательно только в спец ПО (которое уважает карты) получить функцию связи GRV и положения контакта (или площадь VS положение контакта).
создать генератор случайных чисел нужного распределения с помощью генератора для равномерного очень просто - если CDF - это нужная накопленная функция распределения, то CDF^(-1)(генератор для равномерного) - это генератор для этого распределения.
Crystall Ball или @RISK в данном случае не имелся в виду. Только стандартные средства. Обратная функция CDF^(-1) это достаточно общий способ, она не всегда выражается просто, но можно заменить табулированным вариантом.
Однако если посмотреть формулы преобразования Бокса-Мюллера (нормальное распределение, ну а там и до лог-нормального рукой подать) и еще для перехода от равномерного к треугольному распределению, то этого вполне достаточно, чтобы обойтись стандартным генератором случайных чисел в Ёкселе. И даже усеченное распределение и условное Гауссово можно смоделировать.
Для совсем больших выборок конечно надо что-нибудь получше, например, на www.r-project.org (там разумеется все свободно распространяемое) можете поискать что-то вроде Gibbs Sampling.
а вот отсюда можно по подробнее... я впервые все это делаю и как сказать интерес - ОГРОМНЫЙ!!!
GRV- что за функция?
если использывать хакнутый хрустальный шарик будет легче чем в самом петреле.
Интересует любой способ. то есть самый лучший метод, в то же время не так замороченный как в петрел. Думаю все таки все концы сводятся к экселю..???
GRV - gross rock volume - имелась ввиду функция объема породы в нефтегазонасыщенной части в зависимости от положения контакта.
Самого лучшего метода не бывает. А ексель - это самы простой. Если в гугле наберешь monte carlo filetype:xls то найдешь кучу екселевских файлов с примерами использования монте карло. если xls поменять на ppt или pdf, то поиск будет по соответствующим файлам - можно лекции на эту тему найти.
тебе нужно сделать тоже самое, только для подсчета запасов (если я правильно тебя понял) - сгенерировать в екселе кучу вариантов запасов, варьируя параметры в формуле подсчета.
В конечном итоге я должен получить кучу запасов и разделить из на р10, р50, р90.
Спасибо за советы!
вы случана не Сергей Тятушкин - schlumberger?
Коллеги,
подскажите по анализу неопределенности структурной модели.
Как получить набор поверхностей, чтоб сделать картинку как в мануале?
ну т.е. я задал ст.отклонение, задал case, 5 реализаций. Во время расчета вижу, что поверхности меняются. Но после завершения у меня один грид и 5 cases, а как из этого получить горизонты как на разрезе, т.е. 5 поверхностей с расхождением в межскважинном пространстве?
Если есть более подробный мануал (материалы с одноименного курса или хотя бы что-то, но именно по софту, а не по сути) то буду премного благодарен.
В workflow выкидывайте в Input поверхность с куба на каждом шаге.
Convert horizon to Surface
Спасибо большое.
Но все-таки они могли бы и написать об этом в мануале. А то я все галочку искал какую-нибудь...
Всем добрый день!
Может кто-то знает.. Есть ли в РФ обучающие курсы по анализу неопределенностей в Petrel, которые проводятся сторонними организациями (не SLB)???
P.S. в SLB группа не набирается.. поэтому приходится искать альтернативные варианты.