本帖最后由 张鹏83 于 2020-11-9 09:16 编辑
油藏数值模拟可以用来研究复杂的地下变化过程。储层流体受到地下流体长期的冲刷可能会造成储层物性不同程度的变化，储层初始的渗透性和相渗曲线已经不能反应当前油气水的流动特征，常规的数值模拟技术可以通过修改不同时间步的传导率来模拟物性的变化，但是复杂的情况难以考虑到，例如以过水倍数作为参考条件对相渗曲线、渗透率等进行二次或多次调整。
那么，在tNavigator中如何实现时变模拟呢？本文以通过过水倍数为例，表征流体对储层不同程度的冲刷程度，进而进行时变模拟。
一、相渗时变
1、在SOLUTION关键字部分添加下面的关键字，输出累积过水量。
RPTMAPS
FLOWW /
2、然后用tNavigator重新加载模型并运算，在Calculated Maps中可以看到Total flux water（FLOWW）的场图，即累积过水量，然后用这个参数除以孔隙体积即为过水倍数。过水量场图查看：
通过自定义map计算查看过水倍数场图，过水倍数=过水量除以孔隙体积，如下：
3、模型中准备两条及上的相渗曲线
4、在SCH文件中起始时间步输入下面的内容
--2井为工区投产的第一口井，只要2井的井底流压大于-1时，在此后的任意时间步，当过水倍数大于0.2倍时，对应网格自动调用第二条相渗曲线，井底流压的判断标准设置为负数，是因为关井后的WBHP为0，这么写是为了无限降低时变判断的触发条件。
ACTIONW
A1 2 WBHP > -1 9999 /
ARITHMETIC
SATNUM =IF (FLOWW/PORV>0.2, 2, SATNUM) /
/
ENDACTIO
另，如果还要同时考虑饱和度和压力作为约束条件，上述ARITHMETIC可以这样写:
SATNUM =IF (FLOWW/PORV>0.2 & SOIL < 0.5 & PRESSURE < 120, 2, SATNUM) /
以下为相渗时变全过程视频：
https://v.qq.com/x/page/s0941u7bwud.html
二、传导率时变
如果要对传导率进行时变，以X方向为例，上述ARITHMETIC可以这样写：
MULTX =0.00608*(FLOWW/PORV)^-1.1522 /
三、渗透率时变
如果要对渗透率进行时变，以PERMX为例，上述ARITHMETIC可以这样写：
PERMX =PERMX*0.00608*(FLOWW/PORV)^-1.1522 /
细心的用户发现，时变的渗透率是在每一个时间步都会用上一个时间步的渗透率乘以“0.00608*(FLOWW/PORV)^-1.1522 ”
会导致渗透率场被无限倍乘而放大，可以通过下面的方法解决：
1、初始PERMX场-右键-export-以关键字ARRPERMX的场图格式输出；
2、输出场图作为include文件包含到GRID部分，成功后可以在软件界面Initial场图中看到ARRPERMX；
3、然后在ARITHMETIC中这样写即可：
PERMX =ARRPERMX*0.00608*(FLOWW/PORV)^-1.1522 /
也就是每一步都是调用初始渗透率参与时变换算，而不是调用上一步的渗透率场时变。
4、还可以在此基础上设置上限防止渗透率溢出理论范围：
ARITHMETIC
PERMX = IF (PERMX >10000, 10000, PERMX) /
/
如果需要演示工区，请联系我们：
support@colchispetro.com
更多tNavigator技巧及应用，欢迎加入tNavigator技术交流QQ群：777385939。 更多石油行业前沿资讯与应用，请关注微信公众号“科吉思石油技术”！
| 
[复制链接]
发表于 2020-3-28 22:02:54
收藏
