Июн 14
Уважаемые коллеги,
В эклипсе при расчёте выдает такой вот warning:
NEGATIVE COMPRESSIBILITY FOUND IN GAS
PRESSURE TABLE 1 AND OIL PRESSURE TABLE 1 AT A SAMPLE PRESSURE VALUE 176.27586. ADJUST SATURATED FLUID PROPERTY VALUES AT THIS PRESSURE. NEGATIVE COMPRESSIBILITIES OCCUR FOR GAS SATURATIONS LESS THAN 0.12173.
Но в PROPS у меня не стоят отрц. значения. С чего он вдруг такой warning выдает? Серьезно ли это и как можно решить эту проблему, если кто сталкивался с таким.
С Уважением,
ECLIPSE NOOB.
Опубликовано
25 Июн 2014
Активность
15
ответов
5886
просмотров
6
участников
17
Рейтинг
Eclipse интерполирует и экстраполирует табличные данные (линейно). Покажите сами PVT данные из модели.
скорее всего диапазон давлений в таблице не покрывает весь диапазон изменения давления во время расчета модели, поэтому при линейной экстраполяции сжимаемость уходит за горизонт.
Все пропало - ничего сделать нельзя, смирись =))
А если серьезно, то попробуй сделать DEBUG с выводом PVT. Он покажет табличку и место где у него отрицательная сжимаемость получается. Помню на египетском форуме качал эксельку, которая помогает по твоим пивитям сделать все правильно. Там реально хватает изменить третий знак после запятой иногда и все норм. Но нужно смотреть конкретный твой случай.
Как правило, когда у меня начинал писать такие сообщения, через пару лет расчета выскакивал BUG.
Он, вроде как, выскакивает за ось когда экстраполирует значения 1/Bg или 1/Bo и у него получаются отрицательные значения.
ROCK -- Generated : Petrel
287.5 2.0E-005 /
PVTW -- Generated : Petrel
287.5 1.017 4.69E-005 0.362 0 /
PVTO -- Generated : Petrel
0.0002 1 1 4.906 /
44.93 50 1.2061 2.131 /
89.76 100 1.4122 1.39 /
112.18 125 1.5152 1.197 /
134.59 150 1.6182 1.017 /
140.68 157.9 1.6553 0.961
200 1.6355 0.967
250 1.6106 0.974
287.5 1.5920 0.98
300 1.5858 0.982
350 1.5610 0.989
400 1.5361 0.997 /
/
PVDG -- Generated : Petrel
1 1.1239 0.0118
20 0.054007 0.012
40 0.025953 0.0124
60 0.016687 0.0129
80 0.01213 0.0136
100 0.0094628 0.0144
120 0.0077445 0.0155
140 0.0065676 0.0168
160 0.0057257 0.0181
180 0.0051029 0.0196
200 0.0046289 0.0211
250 0.00385 0.025
300 0.0033601 0.029
/
DENSITY -- Generated : Petrel
829 1020 0.949 /
INCLUDE -- Generated : Petrel
'HM_SWAT_INIT.GRDECL' /
FILLEPS -- Generated : Petrel
SWOF -- Generated : Petrel
0.2 0 1 13.333
0.35 0 1 10.833
0.4000 0.0025 0.9375 10.000
0.4300 0.0040 0.9000 9.500
0.5100 0.0105 0.4000 8.167
0.5200 0.0115 0.3600 8.000
0.5500 0.0150 0.2230 7.500
0.5900 0.0200 0.1000 6.833
0.7000 0.0374 0.0130 5.000
0.777 0.15 0 4.725
1 1 0 0.000
/
SGOF -- Generated : Petrel
0 0 1 0
0.06 0 0.47484 0
0.091 0.010264 0.281 0
0.128 0.02323 0.132 0
0.167 0.041598 0.059 0
0.209 0.068069 0.025 0
0.252 0.094 0 0
0.6 0.516 0 0
/
Попробую DEBUG, спасибо)
Это получается нужно будет ставить этот код в PROPS секции и просмотреть .PRT потом?
DEBUG
0*8 1 0*51 /
/
Debug можно ставить в любом месте где хотите, насколько я помню. Создается в корневой папке с моделью файл *.DBG. Его и отрываете в текстовом редакторе =) Там все будет.
Еще можно попробовать:
EXTRAPMS
3 /
тогда он какую то часть инфы выведет в PRT файл в виде ворнингов или типо того.
A complete table of total compressibility evaluation can be written to the .DBG
file by activating switch 21 of the DEBUG keyword in the PROPS section
Проще пойти по пути, который предложил AGA
Включаешь 21 параметр в DEBUG, отвечающий за выдачу total compressibility check. DEBUG указываешь в любой секции. После расчета получаешь файл с расширением DBG, открываешь его, ищешь отрицательные значения в Tgas или в Toil. В соответсвующей строке значения Bg или Bo, которые нужно изменить.
Чтобы подобрать значение Bg или Bo, можно воспользоваться формулами, по которым считается Tgas и Toil.
Все множители будут содержаться в таблице DBG.
Затем подобранное значение переносишь к себе в соответствующую пвт таблицу
Можно конечно ничего из перечисленного не делать, а сразу пробовать изменять Bg в данном случае (сотые, тысячные доли) напротив соответствующего давления. Изменяя значение, не забывать о монотонности.
Максимальное давление по вашим таблицам 400, соответственно до этого значение симулятор строит таблицу сжимаемостей. Таблица сжимаемостей строится для насыщенного состояния, а у вас максимальное давление для насыщенных состояний по таблице 157, т.е. дальше он экстраполирует. На давлении 176 спотыкается и получает отрицательную сжимаемость. Если у вас такое давление ни в одной ячейке модели не достигается при расчете, то ничего страшного не произойдет. Иначе надо править (проще всего дополнить таблицу PVTO фиктивными значениями при больших давлениях насыщения).
Всё верно до этих слов. У него отрицательная сжимаемость при экстраполяции. Ему надо не править, а дописать новые строчки в PVTO. А газ ни при чем, он задан до давлений много выше 176, и при этом сухой.
Понял, буду использовать вышеперечисленные методы!
Всем огромное спасибо :)
Вообще-то странная модель, пластовое по-видимому там 287.5, при этом судя по всему недонасыщенная нефть с давлением насыщения 157.9, и при этом же сухой газ. Обычно в моделях с газовой шапкой нефть в насыщенном состоянии (или близко к тому).
а при чем здесь газовая шапка? начальное пластовое 287,5 - состояние однофазное жидкое (нефть). давление насыщения 157,9, при этом газ который стал свободным принимается сухим. по крайней мере я так понял
все верно