对于高含蜡原油在注冷水开发时，原油会因为遇冷粘度上升，会在注水井周围迅速形成高粘流动阻力区影响注水井的注水效果。传统做法会定义一个复杂的Asphaltene组分模型来明确各种蜡组分和温度的响应关系来模拟整个析蜡、结蜡、絮凝、沉淀从而引起粘度剧烈上升的过程，但这种做法非常复杂需要大量的原油组分测量参数且运算速度很慢，这里我们可以通过Thermal的关键字在黑油模型里根据原油粘温曲线和相渗端点值随温度变化的简化等效处理方式即可达到同样的开发模拟效果，并且数模运算速度和效率相比Asphaltene组分模型都有大幅提高。以下为相关关键字定义：

RUNSPEC

BLACKOIL

THERMAL （调用变温度油藏模拟功能）

GRID

THCROCK （定义岩石热传导率）

10800*206 /

THCOIL （定义原油热传导率）

10800*30.24 /

THCWATER （定义地层水热传导率）

10800*81.818 /

HEATCR（定义岩石的热容系数1）

10800*1918.17 /

HEATCRT（定义岩石的热容系数2）

10800*1.0 /

（说明：岩石的热容值为，其中包括两个系数：CR0通过HEATCR定义，CR1通过HEATCRT定义。）

ROCKPROP

1 86 173.33 206 /

Item1：岩石类型，最大值定义在RUNSPEC的ROCKDIMS中

Item2：初始油藏温度

Item3： 岩石热传导率

Item4：岩石热容值

ROCKCON

3 1 20 7 20 27 27 'K+' /

Item1：岩石类型，最大值定义在RUNSPEC的ROCKDIMS中

Item2：与连接的网格I方向下限编号

Item3：与连接的网格I方向上限编号

Item4：与连接的网格J方向下限编号

Item5：与连接的网格J方向上限编号

Item6：与连接的网格K方向下限编号

Item7：与连接的网格K方向上限编号

Item8：连接的网格的传热方向

PROPS

ENPTVT

--Temperature Swi Swc Swmax Sgi Sgc Sgmax Sorw Sorg

24 1* 1* 1* 1* 1* 1* 0.45 1*

38 1* 1* 1* 1* 1* 1* 1* 1*

49 1* 1* 1* 1* 1* 1* 1* 1*

60 1* 1* 1* 1* 1* 1* 1* 1*

71 1* 1* 1* 1* 1* 1* 1* 1*

86 1* 1* 1* 1* 1* 1* 0.28 1*

/

（说明：上述关键字定义了相渗曲线中随温度而改变的端点值，第一列是温度，第二列是束缚水饱和度，第三列是临界含水饱和度，第四列是最大含水饱和度，第五列是束缚气饱和度，第六列是临界含气饱和度，第七列是最大含气饱和度，第八列是临界油水饱和度，第九列是临界油气饱和度。）

ENKRVT

24 1 1 1 0.15 1* 1* 1

38 1 1 1 1* 1* 1* 1

49 1 1 1 1* 1* 1* 1

60 1 1 1 1* 1* 1* 1

71 1 1 1 1* 1* 1* 1

86 1 1 1 0.28 1* 1* 1

/

（说明：上述关键字定义了另一组相渗曲线中随温度而改变的端点值，第一列是温度，第二列是最大水相相对渗透率，第三列是最大气相相对渗透率，第四列是最大油相相对渗透率，第五列临界含油饱和度下的水相相对渗透率，第六列临界含油饱和度下的气相相对渗透率，第七列临油相相对渗透率（1-最小含水饱和度-临界气饱和度），第八列临油相相对渗透率（1-临界含水饱和度-最小含气饱和度）。）

SPECROCK

25.0 1918.17

86.0 1918.17

/

（说明：第一列是温度，第二列是对应温度下的岩石热容系数）

SPECHEAT

86 2.09 4.18 2.0/

Item1：油藏温度；

Item2：油相热比值；

Item3：水相热比值；

Item4：气相热比值；

OILVISCT

--Temp OIL

25 1418.0

30 913.5

35 446.8

40 118.3

45 19.0

50 9.5

55 7.2

60 5.2

65 4.7

70 3.9

86 3.2

/

（说明：第一列是温度，第二列是对应温度下的地层原油粘度）

WATVISCT

15.5 1

26.6 0.9

38 0.8

49 0.7

60 0.6

71 0.5

86 0.4

/

（说明：第一列是温度，第二列是对应温度下的地层水粘度）

VISCREF

217 /

（说明：指定地层流体的标准压力）

RTEMP

86 /

（说明：指定地层流体的初始温度）