Фев 17
Коллеги,
подскажите, пожалуйста, какие алгоритмы используются при построении карт изобар с учетом разломов?
и если знаете, пожскажите, как они работают?
может есть какая либо литература, источники?
С уважнием,
Андрей Федоров
Опубликовано
08 Фев 2017
Активность
18
ответов
5859
просмотров
7
участников
3
Рейтинг
Карту по точечным скважинным данным планируете строить или по значениям давлений в гриде?
По точечным
При построении карты с разломами нужно не забыть убрать галочку - использовать Z-значения полигонов, представляющих разломы
мне тоже этот вопрос интересен. В Петреле карты изобар можно строить с учетом разломов или непрницаемых границ, обозначенных по ГДИ, интересует симбиоз точечных замеров и грида.
Не совсем понял о чем речь про симбиоз...в карту средних давлений из грида внести корректировки по точечным замерам?
Спасибо, vovlad, за ответ.
я больше имел в виду какие алгоритмы используются, то есть методы интерполяции.
я так понимаю Кригинг не учитывает разломы при построении карт, а вот чтобы разломы каким то образом учитывались, в БАСПРО и в Атласе кажется есть метод ПМФ, который как то учитывает разломы, как на Ваших скринах. вот интересно, как эти или подобные ему методы работают? + интересно в Petrel какой метод?
Насколько я понял в Петреле только Convergent Interpolation может карты с разломами строить. В других алгоритмах эта функция не реализована, почему так - другой вопрос.
да, с локальным перестроением по точечному замеру например
Посмотри эту статью - может поможет. Вроде если зарегистрироваться, то можно прочесть for free.
http://www.jstor.org/stable/44001549
Обсуждали в Телеграмме, что в мануале к Петрелу ссылаются на эту статью: M. A. Haecker. Convergent gridding: A new approach to surface reconstruction
В данном случае при построении карты в Петрел не все так очевидно... точечные данные единичны, а данных при построении с грида очень много. И даже есть указывать разные веса, точечный замер будет вносить незначительное влияние. Если только попытаться разрядить данные полученные с грида...
Подумать надо, так сразу с ходу не скажу. Может, кто еще что подскажет
В петреле нельзя строить карты давлений, т.к. давление изменятся по логарфмическому закону, а все методы больше тяготяют к линейнм зависимостям (кригинг, мин кривизны, полиномы, сплайны и пр.)
Справедливый комментарий. Но насколько помнится, гриды, выгруженные с модели, вроде все равно больше были похожи на результат минимума кривизны(?) или подобное. Может у кого по-другому было. Если подскажешь конкретный алгоритм (или статью), то можно было бы попробовать реализовать.
У меня нет для вас готового решения. Какие пути решения возникают в голове:
-пакет Изолайн, есть алгоритм "метод обратных расстояний", написано что подходит для построения карт давлений, выглядит правдоподно, но само ПО доступно немногим;
-в Петреле есть алгоритм Moving Average, опцией Point weighting можно попробовать добиться нужной картины, для учета разломов построение нескольких карт с данными по каждому блоку отдельно, далее сшивка калькуляторм;
- логарифмирование данных, построение "обычным" алгоритом, далее обратное потенцировние карты, математически вроде всё верно.
ЗЫ по этой же причине нельзя распространять заложенными алгоритмами проницаемость, можно только пересчитывать её через зависоимость Кпр=f(Кп,..)
to FullChaos: С проницаемостью как раз проще, взял логарифм, построил, а затем снова экспоненту. Хотя я абсолютно согласен с тем, что ее надо пересчитвать из пористости.
Хорошо, попробуем скоро первый и третий способы.
По моему для применения различных методов интерполяции, включая и интерполяцию любой поверхности с разломом, нет ничего лучше, чем использование программы Surfer Golden Software. Естественно, потребуется выгрузить грид из Петрел, выполнить соответствующие процедуры в Surfer, а затем загрузить полученный грид обратно. Если это конечно возможно (я не работал с Петрел).
чота вы господа намешали тут всего в кучу.
Во-первых, начнем с того что геометрия разломов может быть разной, а значит и алгоритмы для картопостроения также будут разными.
Во-вторых, при извлеченит карты из грида просто переписываются значения из центра ячеек в узлы двухмерной сетки. Экстраполяция может быть только если карта задается большего размера чем грид.
Так можно же тоже задавать логарифмичсескую трансформацию?! А если зависимость не линейная ?