Skip to main content

CMG操作简介

前处理

数据文件的建立

查看数据文件实例

CMG通过读取一个数据文件，利用数据文件中的关键字指示进行相应的数据计算。CMG数据文件主要包括以下几个部分：

  • 输入、输出控制段

    • 这部分内容包括：确定输入、输出文件名；单位标准（国际标准单位、矿场单位、实验室单位）；输入、输出内容等等。

  • 油藏描述段

    • 这部分内容主要是把建模输出的静态参数包含进来（include），定义你所应用的坐标系。如果是多个压力系统还要用关键字sector定义每个压力系统 包括的部分。

  • 流体组分性质部分

    • 流体组分模型；高压物性(PVT)；油藏温度；岩石压缩系数；原油压缩系数；原始饱和压力；底层原油体积系数等。

  • 岩石—流体性质部分

    • 相对渗透率的定义等。

  • 初始条件数据段

    • 主要包括压力和饱和度的定义。

  • 数值计算方法数据控制段

    • 设置运算限制条件，如最大时间步长、最大时间步数等等，建议用缺省值。

  • 井数据段

    • 井位置的定义（网格定义），井类型的定义，井生产或注入的限制条件，射孔，动态数据等等。

历史拟合

为了取得跟油藏实际动态相一致的一组油藏参数，可以把模拟计算的动态跟实际动态相比较，这种方法叫历史拟合。

历史拟合过程实际是参数校正的过程，在储量拟合的基础上进行单井和井组拟合。根据实测的动态数据，主要拟合以下参数，油藏地层压力及综合含水，单井井底压力的变化，单井见水时间及含水率变化，单井生产指数变化等。

历史拟合的步骤如下：

    (一)、对数模参数全面检查使模型数据体通过运算

    (二)、确定参数的可调范围

    首先清楚哪些参数是确定的，哪些是不确定的。然后根据情况确定可调范围。

    孔隙度：此参数由测井解释和岩心分析得出，视为确定参数，允许改动范围在3%，一般不做修改。

    渗透率：渗透率在任何油、气、田都是不定参数。这不仅是由于测井解释的渗透率值和岩心分析值误差大，而且根据渗透率的特点，井间的渗透率分布也是不确定的。因此对渗透率的修改允许范围较大，可以增大或缩小几倍乃至几十倍。

    有效厚度：油层测井解释的有效厚度与取心井资料对比，一般偏高，主要是钙质层和泥质夹层没有完全挑出来。在调整有效厚度（范围30%）拟合储量后，拟合其他参数时，一般不要再改动。

    初始流体饱和度和初始压力：认为是确定参数

    相对渗透率曲线：认为是不确定参数。

    油气PVT：认为是确定参数。

    油水、气水界面：允许在一定范围内修改，在储量 拟合上后，不再做改动。

    (三)、全区储量拟合

    影响储量的参数有油水界面、气水界面、有效厚度、孔隙度。一般通过校正油水或气水界面以及有效厚度来拟合储量，孔隙度最好不要做改动。在储量拟合好后，这几相参数基本上是确定的参数，必要时做些局部微调。

    (四)、全区地层压力拟合

    影响地层压力的参数很多，几乎所有的物性参数都能影响其变化，如果实际的测压资料由于关井时间不够等因素而不够准确，那么建议拟合压力趋势，不要求追求每个压力值都拟合上。影响地层压力的参数有：孔隙度、有效厚度、压缩系数、原油体积系数、原油密度。

    (五)、单井的生产指数拟合

    这部分工作是最耗时间的工作，主要是拟合含水问题，调节的参数主要是渗透率和相渗，其它参数作为辅助。渗透率的调节参照试井分析以及压裂酸化情况等资料。

    相渗调节的一些规律为：

指标预测

    方案内容

    在历史拟合基础上进行油田开发方案分析研究，预测阶段将提供一个灵活的、功能强的模型，首先是油田开发决策、部署论证分析，对开发原则进行调整，根据剩余饱和度提供新井井位，或对油、水井进行调层，调整液流方向等措施研究，提出不同的提高采收率的开发方案。

    布置方案注意事项

    预测前要做如下生产界限规定：

    • 单井产量（注水量）
    • 含水关井条件
    • 气油比关井条件
    • 最大注入水量界限
    • 井底压力保持数据

    方案调整方法

    老油区经常遇到的是如何控制地层压力问题，研究剩余油分布如何挖潜油层潜力；如何钻加密井、补打新井井位；注采系统如何调整，改变液流方向，调整层间矛盾等问题。

    维持地层压力：油藏工程根据动态分析，指定产量预测方案。在模拟预测阶段，通过注水井来补充地层压力损耗，应用注采比可以调整压力恢复程度。

    加密井或局部新井投产：新井的井位根据剩余油分布来确定，新井工作制度按地层条件，给定产液量再限制一个合理流压，这样新井按定产液量生产，如果达不到产液量，则会按流压限制，产液量自动调整。

    单井调整  如何调整油井见水、高含水，一般分为两类：

      第一类为不改变现有井开发方式，不钻新井的开发调整

    • 改善油、水井的水动力学特征（开采阶段补孔、压裂、酸化、表面活性剂处理）在模型上主要是调整开采层的渗透率(k)
    • 封堵油井高含水层，在模型上关闭高含水层
    • 改变工作制度，强采或周期性采油，在模型上也比较好解决放大压差或改间歇开采等
    • 油井转注，在模型上通过定义井的类型来实现

      • 第二类是通过局部或全面的改变注水方式钻新井调整注水动态

    • 新区油田开发的数模与老区区别在于新区偏重于开发方案的可行性评价、井网层系模式研究、开采方式、注水 方式的选择，通过对比不同速度、不同井网模式、不同方式对开发方案的影响，提供最佳方案，而老区偏重于历史拟合和机理研究，主要目的为油田调整提供可靠的依据

数据文件实例

RESULTS SIMULATOR STARS 200900                  **运行版本

INUNIT SI                           **单位制，SI国际单位
WSRF WELL 1                         **控制由 *OUTSRF 标识的信息写入模拟结果输出文件的频率
WSRF GRID TIME                      **将有TIME控制的时刻的网格信息输出
WSRF SECTOR TIME
OUTSRF GRID PRES SG SO SW TEMP      **输出网格信息的压力，油气水饱和度和温度信息
OUTSRF WELL LAYER NONE              
WPRN GRID 0                         **控制将信息写入输出打印文件的频率
OUTPRN GRID NONE
OUTPRN RES NONE
*****************************************************************************
**Definition of fundamental corner point grid

GRID CORNER 16 20 8
CORNERS 
 21430970.6823 2*21431000.5422 2*21431030.4020 2*21431060.2618 2*21431090.1216
 2*21431119.9814 2*21431149.8412 2*21431179.7010 2*21431209.5608 2*21431239.4206
 .............                              **此处省略N行

REFINE 8 11 1:8 INTO 5 5 1                     **局部网格加密，I方向5份、J方向5份、K方向1份，
REFINE 4 17 1:8 INTO 5 5 1     **必须在NULL和POR关键字之前使用REFINE

*INCLUDE 'NPEDC_PERM'                              **将渗透率数据文件包含在主数据体中
*INCLUDE 'NPEDC_PORO'
*INCLUDE 'NPEDC_NTG'
END-GRID


ROCKTYPE 1                                         **岩石类型，此处只定义了一种
PRPOR 15440                                        **数值参考压力
CPOR 2.0155e-5                                     **岩石压缩系数
ROCKCP 2.89e6 0
*THCONR 189604.8                                   **传热系数，R为rock
*THCONW 53481.6                                              **W为water
*THCONO 21479.04                                             **O为oil
*THCONG 1.496E5                                              **G为gas


HLOSSPROP OVERBUR 2.89e6 1.9e5
         UNDERBUR 2.89e6 1.9e5
**********************************************************************************************************
**$ Model and number of components
MODEL 3 3 3 1                               **表示用于流体火速据输入准备的每个组分类型和数目
COMPNAME 'Water' 'Dead_Oil' 'Soln_Gas'      **组分名
CMM         0       0.33719 0.0191202       **分子量 kg/mol
PCRIT       0       0       4598.86         **临界压力  kPa
TCRIT       0       0       -81.9746        **临界温度  C
KV1         0       0       529660          **气/液 K 值关系式的系数
KV2         0       0       0.00465442 
KV3         0       0       -29.3983 
KV4         0       0       -879.991 
KV5         0       0       -265.99 
PRSR 101            **参考压力  kPa
TEMR 48.9           **参考温度  C
PSURF 101           **地面条件压力    kPa
TSURF 16.85         **地面条件温度    kPa
MASSDEN 995.895 908.894 274.934         **密度
CP  4.28929e-007 2.0155e-006 2.0155e-006        **液体压缩系数    1/kPa
CT1 0.000333833 0.000719715 0.000719715         **热膨胀关系式的第一个系数  1/C
VISCTABLE                                       **粘度与温度关系表
**$      temp                                   
            5    1.69332  4.04055e+006   1099.91  
           15    1.26532        442812   341.199  
           25   0.993501       68527.7   126.993  
           35   0.802826       14054.5   54.8475  
         48.9   0.625695       2262.12   20.8222  
           50   0.614052       1984.67   19.4109  
           70   0.452454       552.113   13.8701  
           85   0.374744       159.148   6.93522  
          100   0.314717        61.787   4.22147  
          120   0.259936       27.3199   3.01388  
          150   0.205599        6.6879   1.25585  
          190   0.163056       1.49889   0.51087  
          230   0.134737      0.623812  0.301449  
          270   0.114383      0.362186  0.217298  
          310  0.0982562      0.254479  0.175683  
          350  0.0861685      0.200602   0.15222  
          390  0.0766296      0.170003  0.137766  
          430  0.0687802       0.15104  0.128282  
          470  0.0622182      0.138523  0.121759  
          510  0.0566509      0.129863  0.117107  
          550  0.0518681      0.123648  0.113691  
          590  0.0477147      0.119056  0.111123  
          630  0.0440743      0.115582  0.109153  
          670  0.0408574      0.112902  0.107617  
          710  0.0379941      0.110798  0.106401  
          750  0.0354291      0.109124  0.105426  
          790  0.0331181      0.107774  0.104636  
          830  0.0310252      0.106673  0.103989  
**确定对液相粘度的非线性混合规则
VSMIXCOMP 'Soln_Gas'
VSMIXENDP 0.0007 0.2        
VSMIXFUNC 0.0007 0.0141506 0.0292906 0.0459119 0.0637941 0.0828532 0.102978 0.124055 0.146149 0.169306 0.2 

**********************************************************************************************************
ROCKFLUID           **岩石流体数据的开始
RPT 1 WATWET        **对岩石流体数据定义岩石类型号，并且用 KRTYPE 对网格设置岩石类型号
**$        Sw       krw      krow
SWT                                                             **油水相对渗透率表
     0.371429           0          1
     0.402857      0.0005   0.658046
     0.437143       0.001   0.454023
     0.471429  0.00287356   0.341954
     0.505714  0.00574713   0.227011
     0.542857   0.0114943   0.158046
     0.577143   0.0201149  0.0890805
     0.608571   0.0316092  0.0316092
     0.682857   0.0517241          0
**$        Sl         krg        krog
SLT                                                             **油气相对渗透率表
     0.587342    0.222222           0
     0.670886   0.0740741           0
     0.711375    0.037037  0.00673401
     0.751839    0.020202    0.023569
     0.792252  0.00673401   0.0606061
     0.832596           0    0.124579
     0.872838           0    0.228956
     0.915441           0    0.400673
     0.955317           0    0.649832
            1           0           1
KRTEMTAB SWR SORW                                           **给定临界S与端点对T的依赖
**$      TEMP       SWR      SORW                           **SWR束缚水饱和度，SORW残余油饱和度
           50       0.4      0.32
          200       0.5      0.18

**********************************************************************************************************
INITIAL                                                      **初始化
VERTICAL DEPTH_AVE                                           **利用DWOC数据完成按深度平均的PC-重力垂向平衡计算

INITREGION 1                                                 **分区，此处只有一个区
REFPRES 15520                                                **参考深度下的参考压力
REFDEPTH 1170                                                **参考深度
DWOC 1238                                                    **油水界面
**********************************************************************************************************
NUMERICAL                        **用来指示后面全是数值方法控制数据
RUN

DATE 2010 8 9
DTWELL 1e-006           **最小步长
NORM PRESS 500.0                                **控制压力最大变化值为500

WELL  'I1'
   *INJECTOR *MOBWEIGHT  'I1'          **注入井定义
   *INCOMP WATER  1.0  0.0  0.0                 **注入流体定义
   *TINJW 300                                **蒸汽温度
   *QUAL .5                                  **蒸汽干度
   *OPERATE *MAX *STW   200                     **最大注水200            
   *GEOMETRY  K  0.086  0.249  1.  0.           **井参数         
 
PERF  WI  'I1'                                  **定义井位置
**$ UBA     wi       Status  Connection  
    4 17 1 / 3 3 1  100000.  OPEN    
    4 17 2 / 3 3 1  100000.  OPEN    
    4 17 3 / 3 3 1  100000.  OPEN    
    4 17 4 / 3 3 1  100000.  OPEN    
    4 17 5 / 3 3 1  100000.  OPEN    
    4 17 7 / 3 3 1  100000.  OPEN    
    4 17 8 / 3 3 1  100000.  OPEN    
** 
**
WELL  'P1'                                           
PRODUCER 'P1'                                     **生产井定义       
OPERATE  MAX  STL  40 CONT                        **最大产液为40

GEOMETRY  K  0.086  0.249  1.  0.
PERF  GEOA  'P1'
**$ UBA     ff  Status  Connection  
    4 17 1 / 3 3 1  1.  OPEN   
    4 17 2 / 3 3 1  1.  OPEN   
    4 17 3 / 3 3 1  1.  OPEN   
    4 17 4 / 3 3 1  1.  OPEN   
    4 17 5 / 3 3 1  1.  OPEN   
    4 17 7 / 3 3 1  1.  OPEN   
    4 17 8 / 3 3 1  1.  OPEN   
 SHUTIN 'P1'                                       **关井


DATE 2010 8 21

SHUTIN 'I1'
DATE 2010 8 24
OPEN 'P1'
DATE 2010 9 1  
DATE 2010 10 1
DATE 2010 11 1
DATE 2010 12 1
DATE 2011 1 1
DATE 2011 2 1 
DATE 2011 3 1
STOP

Comments

Post a Comment

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
2 comments
Read more

MatLab中patch函数的基本用法

patch是用来构建多边形的一个基本函数。 用法一 patch(X,Y,C) patch(X,Y,Z,C) patch( 'XData' ,X, 'YData' ,Y) patch( 'XData' ,X, 'YData' ,Y, 'ZData' ,Z) 1.1 说明 patch(X,Y,C)用来构建一个或者多个可填充的多边形，其使用X和Y作为每个点的坐标值，patch将会按顺序连接每个点。如果要得到一个多边形，将X和Y设置为向量；如果要得到多个多边形，将X和Y设置为矩阵，没一列对应一个多边形。C决定多边形的颜色，可以是系统认定的字符，也可以是一个数值，也可以是RGB向量。 patch(X,Y,Z,C)用来构建三维坐标下的多边形。 patch(‘XData’,X,’YData’,Y)和patch(‘XData’,X,’YData’,Y,’ZData’,Z)的用法与patch(X,Y,C)和patch(X,Y,Z,C)的用法类似，只是不设定颜色。 1.2 例子 1.2.1 x = [ 0 1 1 0 ] ; y = [ 0 0 1 1 ] ; patch(x,y, 'red' ) x和y都是1*4的向量，表示将四个点(0,0)、(1,0)、(1,1)和(0,1)依次连接，最后闭合形成一个四边形，设定颜色为红色。 1.2.2 x2 = [ 2 5 ; 2 5 ; 8 8 ] ; y2 = [ 4 0 ; 8 2 ; 4 0 ] ; patch(x2,y2, 'green' ) x2和y2都是3*2的向量，两列表示画两个多边形。第一个多边形连接的点依次是(2,4)、(2,8)和(8,4)，第二个多边形连接的点依次是(5,0)、(5,2)和(8,0)，颜色设定为绿色。 1.2.3 如果上例的三角形第一个是红色，第二个是绿色，那么patch代码修改为 x2 = [ 2 5 ; 2 5 ; 8 8 ] ; y2 = [ 4 0 ; 8 2 ; 4 0 ] ; patch(x2(:, 1 ),y2(:, 1 ), 'red' ) pat
1 comment
Read more

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/   广大学者将自己的文章发表至学术期刊（免费或者支付版面费），然而当需要查看其他学者的文章时还需要向出版商付费，你是不是也觉得这完全阻碍了科学文化的传播。艾尔巴金在为自己辩护时援引联合国《世界人权宣言》第二十七条所提的：“人人有权自由参加社会之文化生活，欣赏艺术，并共同襄享科学进步及其利益。”关于这个问题，另一位被称
3 comments
Read more