Skip to main content

CMG热采STARS模拟器粘度非线性混合的使用

关键字用法:

  *VSMIXCOMP comp_des
  *VSMIXENDP xlow xhigh
  *VSMIXFUNC f(1)...f(11)

  comp_des  在粘度非线性混合中关键组分的名称或序号。
  xlow     对应于第一个表格元素的横坐标。允许的范围是0~xhigh。
  xhigh    对应于最后一个表格元素的横坐标。允许的范围是xlow~1。
  f(1)...f(11)   表格元素定义的非线性混合规则函数,这个函数应该是单调递增的,并且应该是平滑的。

  缺省:
    如果*VSMIXCOMP不存在,则假设对所有组分进行线性混合。
    如果*VSMIXENDP不存在,则假设xlow=0,xhigh=1。
    如果*VSMIXFUNC不存在,则表格输入f(i)=(i-1)/10,i=1到11,对应于0到1之间的线性间隔。

算法:

  首先线性混合的规则为

              
  使用表征非线性的代替,其中a为关键组分。但是仍然需要满足和为1的条件,即
              
  由上式可以得出
              
  因此对于粘度的非线性混合规则可以表示为
              
  对于任意的就有三个范围去得到相应的的值
              
              的值从输入的表中得到
              
 可以用下图表示其含义:
 通过上述非线性混合的式子可以推导出的计算方法
              
              
  将和中间等距间隔的9个值分别代入上述式子中,得到*VSMIXFUNC关键字中需要输入的11个值。
  如果,则非线性混合规则退化为线性混合规则。

查看输出:

  通过*OUTPRN和*GRID的*VISCCMP子关键字,以及*OUTSRF *GRID的*VISWCMP和*VISOCMP子关键字和*SPECIAL检查粘度非线性混合关键组分的组成。

例子:

  以下数据和计算都是在同一温度下进行。一种沥青粘度为,一种纯溶解剂的浓度为,主要组分为溶解剂。其两者混合后的粘度表示为下图所示:

  相应的数据如下表所示:
0.000300000
0.015268143
0.030191506
0.045109375
0.06050090
0.07518580
0.0905816
0.1051808
0.120824
0.135633
0.150604
1.000600
  对于上述关系,溶解剂浓度从0.0到0.15之间可认为是非线性的,根据算法中计算的公式,就可以得到相应的值:
0.00012.610
0.01512.500.018064
0.03012.160.072227
0.04511.600.162360
0.06010.820.288023
0.0759.830.447603
0.0908.670.634493
0.1057.500.822474
0.1206.700.949045
0.1356.450.991380
0.1506.400.998959
  关键字设置为:
    *COMPNAME 'WATER' 'BITUMEN' 'SOLVENT'
    *VISCTABLE
    **TEMP
     75     0     300000   600
    *VSMIXCOMP 'SOLVENT'
    *VSMIXENDP 0.0 0.15
    *VSMIXFUNC
     **0.0  0.015    0.03   0.045   0.06    0.075
       0   0.0018064 0.072227 0.162360 0.288023  0.447603
     **0.09   0.105   0.12    0.135  0.15
      0.634493 0.822474 0.949045 0.991380 0.998959
 注:此例子为CMG帮助手册中的例子。

 一个模拟程序中只有一相可以使用非线性混合规则,即水相或者油相。如果是水相的非线性混合,需要注意的是在*MODEL关键字下的组分数和*COMPNAME关键字下的组分顺序,以及相应的各组分性质关键字的设置顺序。

说明:

  粘度非线性混合选项的计算方式与*DNMIXCOMP关键字说明的密度选项相同,唯一不同的是进行加权的量是粘度的对数。

  以上内容由本人参考帮助文档及模拟经验所得,如有错误还请您不吝赐教。

Comments

Popular posts from this blog

使用PHP Webhook方式打造Telegram Bot

一、找BotFather拿到bot token     在telegram中私聊BotFather建立自己的bot,给bot取名,名字必须要以bot结尾。建好后自己的bot就有一个唯一的token,类似下面的一串字符 164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8 二、确认bot响应文件的位置     在写好bot响应文件后,要把bot放在网络上的一个位置,并且这个位置必须要加密的,即以https开头的一串网址。比如响应文件的名称为telbot.php,把它放在下面这个网址的位置: https://my.webhost.com/ 164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8 /telbot.php 上面网址中的红色设置和bot的token一样是为了确定这个唯一的位置,当然也可以任意设置。 三、告诉Telegram响应文件的位置 Telegram用下面网址的形式来设定webhook响应方式 https://api.telegram.org/bot [myauthorization-token] /setwebhook?url= [myboturl] 按照上面的网址形式,把自己创建的bot的token以及响应文件的位置填入,然后在浏览器中运行一下即可设置成功。比如: https://api.telegram.org/bot164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8/setwebhook?url=https://my.webhost.com/164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8/telbot.php 设置成功后,页面会显示下面的内容: {"ok":true,"result":true,"description":"Webhook is already set"} 四、在Telegram中给自己的bot发消息进行验证 php响应文件例子 <?php  define('BOT_TOKEN', 'YOURBOT:TOK...

telegram中的Sci-Hub机器人,又一文献下载利器

或许你看到标题会问什么是telegram,什么是Sci-Hub?请听我一一道来。 什么是Sci-Hub Sci-Hub是一个线上 数据库 ,其上提供48,000,000篇科学学术论文和文章。网站透过“.edu”代理服务器访问相关页面,每天会上传新的论文文章。2011年,哈萨克研究生亚历珊卓·艾尔巴金(Alexandra  Elbakyan)因为研究论文成本过高,每篇论文在付费墙机制下通常需要花费30美元,而决定成立Sci-Hub。2014年,学术界开始预测网站将会发展为类似Napster的服务。不过到了2015年,学术出版社爱思唯尔向纽约地方法院提交诉讼,指控Sci-Hub已经侵犯版权。纽约地方法院在2015年10月28日仍下令Sci-Hub原本使用的网域名称“Sci-Hub.org”必须终止。爱思唯尔在法院上获得胜诉后,一群研究人员、作家和艺术家则连署一封表态支持Sci-Hub和创世纪图书馆的公开信,声称这次诉讼对于世界各地的研究人员是“重大打击”,并指出:“它同样贬低我们、作者、编辑和读者。它寄生于我们的劳动,它阻挠我们为大众服务,它阻拦我们进入。”而该计划于11月因法院命令中止后,在同一个月内便改用网域名称“.io”重新上线,并开放使用Tor浏览。2016年1月时,Sci-Hub平均每天约有200,000人访问,Sci-Hub则声称网站服务每天平均有数十万次档案请求。  Sci-Hub是目前已知第一个提供大量自动且免费的付费学术论文的网站,使用者不需要事前订阅或付款,就能够使用原本存放在付费数据库的论文文章,并提供搜寻原先出版社网站内的文件档案服务。 以上介绍来源于维基百科词条 Sci-Hub Sci-Hub网站被屡次下线,但是又通过更换域名重新上线。以下三个网址经测试可以使用:  http://www.sci-hub.bz/   http://www.sci-hub.ac/   http://www.sci-hub.cc/   广大学者将自己的文章发表至学术期刊(免费或者支付版面费),然而当需要查看其他学者的文章时还需要向出版商付费,你是不是也觉得这完全阻碍了科学文化的传播。艾尔巴金在为自己辩护时援引联合国《世界人权宣言》第二十七条所提的:“人人有权自由参加社会之文化生活,欣赏艺...

CMG操作简介

前处理 数据文件的建立 查看数据文件实例 CMG通过读取一个数据文件,利用数据文件中的关键字指示进行相应的数据计算。CMG数据文件主要包括以下几个部分: 输入、输出控制段 这部分内容包括:确定输入、输出文件名;单位标准(国际标准单位、矿场单位、实验室单位);输入、输出内容等等。 油藏描述段 这部分内容主要是把建模输出的静态参数包含进来(include),定义你所应用的坐标系。如果是多个压力系统还要用关键字sector定义每个压力系统 包括的部分。 流体组分性质部分 流体组分模型;高压物性(PVT);油藏温度;岩石压缩系数;原油压缩系数;原始饱和压力;底层原油体积系数等。 岩石—流体性质部分 相对渗透率的定义等。 初始条件数据段 主要包括压力和饱和度的定义。 数值计算方法数据控制段 设置运算限制条件,如最大时间步长、最大时间步数等等,建议用缺省值。 井数据段 井位置的定义(网格定义),井类型的定义,井生产或注入的限制条件,射孔,动态数据等等。 历史拟合 为了取得跟油藏实际动态相一致的一组油藏参数,可以把模拟计算的动态跟实际动态相比较,这种方法叫历史拟合。 历史拟合过程实际是参数校正的过程,在储量拟合的基础上进行单井和井组拟合。根据实测的动态数据,主要拟合以下参数,油藏地层压力及综合含水,单井井底压力的变化,单井见水时间及含水率变化,单井生产指数变化等。 历史拟合的步骤如下: (一)、对数模参数全面检查使模型数据体通过运算 (二)、确定参数的可调范围 首先清楚哪些参数是确定的,哪些是不确定的。然后根据情况确定可调范围。 孔隙度 :此参数由测井解释和岩心分析得出,视为确定参数,允许改动范围在3%,一般不做修改。 渗透率 :渗透率在任何油、气、田都是不定参数。这不仅是由于测井解释的渗透率值和岩心分析值误差大,而且根据渗透率的特点,井间的渗透率分布也是不确定的。因此对渗透率的修...