CMG数值模拟软件培训教程(67张)课件

1、CMG数值模拟软件培训教程CMGNumericalSimulationTrainingTutorial北京华美杰润石油新科技Cmgco.,ltd2004.8*第1页，共68页。内容 Outline 前言（Background and objectives） 前处理(3D Modeling& Data Select) 数据文件的建立(File Operating) 历史拟合(History Matching) 指标预测(Index Forecast) 软件演示讲解(Ex)第2页，共68页。一、前言 Foreword利用计算机模拟地下油水的流动，给出某一时刻油水分布，预测油藏动态。油藏数值模拟定义

2、（defining）油藏数值模拟是大规模的动态仿真技术，是高性能计算在油田开发中的重要应用。通过该油藏模拟软件研究剩余油的分布，以便科学发展采油工艺和技术，确保各大油田的稳产。 意义（significance）第3页，共68页。一、前言 Foreword软件简介（introduce）加拿大CMG数值模拟软件包括黑油模型、组份模型和热采模型，稠油模型等的大型模拟软件。第4页，共68页。二、前处理(Data Select)静态资料数模文件（选取参数）动态资料（经过整理）数模文件第5页，共68页。三、CMG主文件 Main File1.1数值模拟方法 数值模拟是应用数学模型把实际的油、气藏动态重现一

3、遍，也就是通过流体力学方程借用计算机，计算数学的求解，结合油藏地质学、油藏工程学重现油、气田开发的实际过程，用来解决油、气田实际问题。第6页，共68页。1.2 CMG模型类型（1） 从油藏特性来分：黑油模型和组分模型 黑油模型是指油质较重的油藏类型模型，可以模拟油、气、水三相存在的油气藏。 组分模型是指油质较轻、气体较富的油气藏，如挥发性油藏和凝析气藏。这种类型油气藏开采过程中碳氢组分是变化的，而相间的变化是受热动力影响的，这是与黑油模型的典型不同之处。第7页，共68页。 （2） 从模型维数来分：两维平面、两维剖面、两维径向以及三维模型。 两维平面模型：主要用于较薄的油层，研究平面上油、水运动

4、规律，不考虑层间的影响。使用这种模型一般研究改变井网、井距、注水方式等平面上变化对采收率的影响，研究平面上非均质变化及断层对开采的影响。 两维剖面模型：主要用于研究层间非均质的影响，它考虑了重力、毛管力及流速对驱油效率的影响，研究纵向上油水运动规律。 两维径向模型：主要用于研究锥进动态，包括气锥和水锥。 三维模型：主要用于全区块的模拟，油层较厚、平面和纵向上非均质比较复杂的油藏。第8页，共68页。二维剖面模型三维空间模型顶面网格图单井径向模型顶面网格图第9页，共68页。1.3 CMG数模数据体的内容输入、输出控制段油藏描述段流体组分性质部分岩石流体性质部分初始条件数据段数值计算方法数据控制段井

5、数据段三、CMG主文件 Main File第10页，共68页。总体结构数模培训教程数据文件讲解RESERVOIR DESCRIPTIONCOMPONENT PROPERTIESROCK-FLUID DATA网格划分及岩石特性等参数流体模型、油气PVT表相对渗透率等参数第11页，共68页。数模培训教程NUMERICALWELL AND RUN DATA数值方法控制段（步长、最大迭代次数等）井数据段：开始时间、关井时间、生产限制、完井层位和井的几何特征等总体结构数据文件讲解INITIALIZATION泡点压力、含油含水饱和度第12页，共68页。1输入、输出控制段 这部分内容包括：确定输入、输出文件

6、名；单位标准（国际标准单位、矿场单位、实验室单位）；输入、输出内容等等。1.3 CMG数模数据体的内容第13页，共68页。2油藏描述段 这部分内容主要是把建模输出的静态参数包含进来（include），定义你所应用的坐标系。如果是多个压力系统还要用关键字sector定义每个压力系统 包括的部分。 注意的一点是建模输出的网格节点数要和数模定义的网格节点一致。1.3 CMG数模数据体的内容第14页，共68页。3流体组分性质段 这部分内容包括：流体组分模型；高压物性(PVT)；其他参数如下表。油层平均有效厚度，m地面原油密度，g/cm3平均有效渗透率，10-3m2地层水密度，g/cm3平均有效孔隙度，

7、%地层水粘度，mPa.s油藏温度，OC天然气相对密度, f岩石压缩系数，10-4/MPa原油压缩系数， 10-4/MPa原始饱和压力，MPa地层水压缩系数，10-4/MPa原始含油饱和度，%地层原油体积系数, f第15页，共68页。4岩石流体性质数据段 这部分内容主要是相对渗透率的定义。第16页，共68页。5初始条件数据段 主要是压力和饱和度的定义，方法有两种： 一、确定每个网格的压力和饱和度，这种方法适合无边、底水气藏。 二、确定油水、油 气界面深度，参考深度的参考压力，饱和压力第17页，共68页。6数值计算方法数据控制段 一些运算限制条件，如最大时间步长、最大时间步数等等，建议用缺省值。第

8、18页，共68页。7井数据段 内容包括：井位置的定义（网格定义），井类型的定义，井生产或注入的限制条件，射孔，动态数据等等。第19页，共68页。数据资料收集静态资料小层数据表或等值线图 砂层厚度 有效厚度（或净毛比） 顶部深度 孔隙度 渗透率等数模培训教程第20页，共68页。静态资料2. 地质储量及地层、油藏特点的总结报告3. 油、气、水高压物性PVT数据4. 油水、油气相渗曲线、数据5. 毛管压力曲线、数据数模培训教程数据资料收集第21页，共68页。静态资料6. 油、气、水分布（原始饱和度），压力分布或油水界面和油气界面7. 井位分布图8. 流体和岩石化验分析报告数模培训教程数据资料收集第2

9、2页，共68页。动态资料1. 射孔完井报告2. 井史报告、压裂等措施3. 系统测压资料4. 试油、试井和试采资料（压力恢复曲线） 数模培训教程数据资料收集第23页，共68页。动态资料5. 油水井别，调整井位示意图6. 油井生产数据报表： 日产油、日产液、日产气、日产水 综合含水 压力 累积产油、累积产气、累积产水、累积产液数模培训教程数据资料收集第24页，共68页。动态资料7. 注水井生产数据报表： 日注水量 累积注水量8. 区块综合生产数据统计报表： 日产油、日产液、日产气、日产水 综合含水 数模培训教程数据资料收集第25页，共68页。动态资料8. 区块综合生产数据统计报表： 采出程度 压力

10、 累积产油、累积产气、累积产水、累积产液 日注水量 累积注水量 数模培训教程数据资料收集第26页，共68页。数据资料处理静态资料数模培训教程数模文件（选取参数）动态资料（经过整理）数模文件第27页，共68页。数据资料处理数模培训教程生产数据表第28页，共68页。数据资料处理数模培训教程生产数据表拟合预测处理第29页，共68页。数据资料处理生产数据表转化程序编写Sub fore ()按格式重新排列Dim rg As String, rg1 As StringDim pre As IntegerDim times As Long, years As LongDim lettters As Stri

11、ngDim tmp As LongDim temp As Longletters = ABCDEFGHIJKLMNOPQRSTUVWXYZtimes = 1i = 1pre = 0For years = 198105 To 198301 Step 1Do If years = 198113 Then years = 198201 If years = 198213 Then years = 198301第30页，共68页。数据资料处理生产数据表转化程序编写 If Sheet1.Range(E & i).Value years Then Exit Do If (i + 1 - pre) 27 T

12、hen rg = Left(letters, (i + 1 - pre) rg = Right(rg, 1) ElseIf (i + 1 - pre) 53 Then rg = Left(letters, (i + 1 - pre) - 26) rg = A & Right(rg, 1) ElseIf (i + 1 - pre) 79 Then rg = Left(letters, (i + 1 - pre) - 52) rg = B & Right(rg, 1) Else MsgBox Error! End If rg1 = rg & times Sheet2.Range(A & times

13、).NumberFormatLocal = Sheet2.Range(B & times).NumberFormatLocal = 第31页，共68页。数据资料处理生产数据表转化程序编写Function data(year As Long, mon As Long) As LongIf mon = 2 Then If Not (year Mod 4) Then data = 29 Else data = 28 End IfEnd If If mon = 1 Or mon = 3 Or mon = 5 Or mon = 7 Or mon = 8 Or mon = 10 Or mon = 12 T

14、hen data = 31 ElseIf mon 2 Then data = 30 End If End Function第32页，共68页。 为了取得跟油藏实际动态相一致的一组油藏参数，可以把模拟计算的动态跟实际动态相比较，这种方法叫历史拟合。 历史拟合过程实际是参数校正的过程，在储量拟合的基础上进行单井和井组拟合。根据实测的动态数据，主要拟合以下参数，油藏地层压力及综合含水，单井井底压力的变化，单井见水时间及含水率变化，单井生产指数变化等。 (1) 历史拟合的定义四、历史拟合 History Matching第33页，共68页。 （一）对数模参数全面检查使模型数据体通过运算； （二）确定参

15、数的可调范围； （三）全区的储量拟合； （四）全区的地层压力拟合； （五）单井生产指数拟合。(2) 步骤第34页，共68页。 在数模数据体完成后，要经过运算生成用于再处理的其它文件。CMG软件：(一)、对数模参数全面检查使模型数据体通过运算数据体文件(文本）IMEX.OUT.IRF.MRF.RRF输出的文本文件索引文件主文件可回绕、可选择(2) 步骤第35页，共68页。 由于模型参数数量较多，可调自由度大，而实际油藏动态数据种类和数量有限，不能唯一确定油藏模型的参数。为了避免或减少修改参数的随意性，必须确定参数的可调范围。 首先分清哪些参数是确定的，哪些是不确定的。然后根据情况确定可调范围。

16、(二)、确定参数的可调范围(2) 步骤第36页，共68页。 孔隙度：此参数由测井解释和岩心分析得出，视为确定参数，允许改动范围在3%，一般不做修改。 渗透率：渗透率在任何油、气、田都是不定参数。这不仅是由于测井解释的渗透率值和岩心分析值误差大，而且根据渗透率的特点，井间的渗透率分布也是不确定的。因此对渗透率的修改允许范围较大，可以增大或缩小几倍乃至几十倍。 有效厚度：油层测井解释的有效厚度与取心井资料对比，一般偏高，主要是钙质层和泥质夹层没有完全挑出来。在调整有效厚度（范围30%）拟合储量后，拟合其他参数时，一般不要再改动。(二)、确定参数的可调范围(2) 步骤第37页，共68页。初始流体饱和

17、度和初始压力：认为是确定参数。 相对渗透率曲线：认为是不确定参数。油气PVT：认为是确定参数。油水、气水界面：允许在一定范围内修改，在储量 拟合上后，不再做改动。(2) 步骤(二)、确定参数的可调范围第38页，共68页。 影响储量的参数有油水界面、气水界面、有效厚度、孔隙度。一般通过校正油水或气水界面以及有效厚度来拟合储量，孔隙度最好不要做改动。在储量拟合好后，这几相参数基本上是确定的参数，必要时做些局部微调。 (三)、全区储量拟合(2) 步骤第39页，共68页。 影响地层压力的参数很多，几乎所有的物性参数都能影响其变化，如果实际的测压资料由于关井时间不够等因素而不够准确，那么建议拟合压力趋势

18、，不要求追求每个压力值都拟合上。影响地层压力的参数如下： (1).孔隙度和有效厚度和压缩系数：这几个参数的提高会减缓地层压力的下降，同时也会增大地质储量。 (2).此外还有其他一些参数影响地层压力，如原油体积系数、原油密度等。(四)、全区地层压力拟合(2) 步骤第40页，共68页。 这部分工作是最耗时间的工作，主要是拟合含水问题，调节的参数主要是渗透率和相渗，其它参数作为辅助。渗透率的调节参照试井分析以及压裂酸化情况等资料，下面把相渗的一些调节技巧讲一下： (五)、单井的生产指数拟合(2) 步骤第41页，共68页。油井见水过早油井见水过晚油藏在初含水时期拟合不好可以做以下处理第42页，共68页

19、。油井后期含水太低油藏在后期拟合不好可以做以下两种处理第43页，共68页。油井后期含水太高第44页，共68页。见水过早第45页，共68页。井底流压拟合综合含水拟合拟合后的结果第46页，共68页。五、 指标预测Index Forecast方案内容 在历史拟合基础上进行油田开发方案分析研究，预测阶段将提供一个灵活的、功能强的模型，首先是油田开发决策、部署论证分析，对开发原则进行调整，根据剩余饱和度提供新井井位，或对油、水井进行调层，调整液流方向等措施研究，提出不同的提高采收率的开发方案。第47页，共68页。预测前要做如下生产界限规定： (1) 单井产量（注水量）； (3)含水关井条件； (4)气油

20、比关井条件； (5)最大注入水量界限； (6)井底压力保持数据。2 布置方案注意事项指标预测第48页，共68页。 老油区经常遇到的是如何控制地层压力问题，研究剩余油分布如何挖潜油层潜力；如何钻加密井、补打新井井位；注采系统如何调整，改变液流方向，调整层间矛盾等问题。 3 方案调整方法指标预测第49页，共68页。(1)维持地层压力 油藏工程根据动态分析，指定产量预测方案。在模拟预测阶段，通过注水井来补充地层压力损耗，应用注采比可以调整压力恢复程度。(2)加密井或局部新井投产 新井的井位根据剩余油分布来确定，新井工作制度按地层条件，给定产液量再限制一个合理流压，这样新井按定产液量生产，如果达不到产

21、液量，则会按流压限制，产液量自动调整。 3 方案调整方法指标预测第50页，共68页。 ( 3)单井调整 如何调整油井见水、高含水，一般分为两类： 一类为不改变现有井开发方式，不钻新井的开发调整： 改善油、水井的水动力学特征（开采阶段补孔、压裂、酸化、表面活性剂处理）在模型上主要是调整开采层的渗透率(k) 封堵油井高含水层，在模型上关闭高含水层。 改变工作制度，强采或周期性采油，在模型上也比较好解决放大压差或改间歇开采等。 油井转注，在模型上通过定义井的类型来实现。3 方案调整方法第51页，共68页。第二类是通过局部或全面的改变注水方式钻新井调整注水动态。 新区油田开发的数模与老区区别在于新区偏

22、重于开发方案的可行性评价、井网层系模式研究、开采方式、注水 方式的选择，通过对比不同速度、不同井网模式、不同方式对开发方案的影响，提供最佳方案，而老区偏重于历史拟合和机理研究，主要目的为油田调整提供可靠的依据。指标预测第52页，共68页。指标预测第53页，共68页。 数值模拟的关键在于模型的准确性，而模型的准确性除了受地质建模提供数据体的影响，历史拟合也是主要因素。所以历史拟合在数值模拟过程是关键。 由于我们建立的油藏模型，未知参数（孔隙度、渗透率、饱和度）远比已知参数（生产数据等）多得多，这是历史拟合的困难之处，换句话说历史拟合有不唯一解。要想把模型做得准确，历史拟合过程中参数调整的合理，需

23、要有丰富的经验以及相关的地质油藏知识。指标预测第54页，共68页。六、软件演示 建立数据文件 历史拟合过程 指标预测结果分析第55页，共68页。六、软件演示一、建立数据文件 *FILENAMES *OUTPUT *REWIND-OUT *INDEX-OUT*MAIN-RESULTS-OUT *INDEX-IN cpu53-1.irf *MAIN-RESULTS-IN cpu53-1.mrf *REWIND-IN cpu53-2.rrf *TITLE1 2-PHASE (OIL-WATER) DEPLETION, CANKAO1 *TITLE2 84 WELLS DEPLETION OPERAT

24、ION *TITLE3 COMPARE WITH WTRFLD & PSEUDO-MIS.FLOOD *CASEID WAT-GAS-OIL标题段第56页，共68页。六、软件演示输入、输出控制段 *-INPUT/OUTPUT CONTROL SECTION- *INUNIT *SI *RESTART *WRST *REWIND *WRST *TIME *WPRN *WELL *TIME *WPRN *GRID *TIME *OUTPRN *WELL *LAYER *OUTPRN *GRID *IMEXMAP *SG *SO *SW *PRES *WSRF *WELL *TIME *WSRF *

25、GRID *TIME *OUTSRF *GRID *IMEXMAP *SG *SO *SW *PRES *OUTDIARY *BRIEF *PRESAQ *HEADER 20 一、建立数据文件第57页，共68页。六、软件演示油藏描述段*-RESERVOIR DESCRIPTION KEYWORDGROUP- *GRID *VARI 52 48 14 *KDIR *DOWN *DI *CON 50.0 *DJ *CON 50.0 *DEPTH-TOP *ALL *INCLUDE F:pushumomoxing3dtop.inc *DK *ALL *INCLUDE F:pushumomoxingd

26、z.inc *NETPAY *ALL *INCLUDE F:pushumomoxingdznet.inc *NULL *IJK *POR *ALL *INCLUDE F:pushumomoxingpor.inc *PERMI *ALL *INCLUDE F:pushumomoxingavk.inc *PERMJ *EQUALSI *PERMK *EQUALSI * 0.1 *PRPOR 23400 *CPOR 4.85E-7一、建立数据文件第58页，共68页。六、软件演示流体组分性质部分 *-COMPONENT PROPERTIES- *MODEL *BLACKOIL *TRES 86 *PV

27、T 1 * P RS Bo EG VISO VISG *DENSITY *OIL 847.7 *DENSITY *GAS 1.0192 *DENSITY *WATER 1000.0 *CO 1.1360E-6 *CVO 3.0120E-5 *BWI 1.030 *CW 7.3822E-7 *REFPW 23400.0 *VWI 1.00 *CVW 0.0一、建立数据文件第59页，共68页。六、软件演示岩石流体性质部分*-ROCK-FLUID DATA KEYWORD GROUP- *ROCKFLUID *RPT 1 *SWT * SW KRW KROW PCOW *SLT * SL KRG K

28、ROG PCOG 一、建立数据文件第60页，共68页。六、软件演示初始条件数据段*-INITIALIZATION- *INITIAL *VERTICAL *BLOCK_CENTER *WATER_OIL_GAS *PB *CON 20470.0 *REFDEPTH 2200.0 *REFPRES 23400.0 *DGOC 2285.0 *DWOC 2440.0 一、建立数据文件第61页，共68页。六、软件演示数值计算方法数据控制段 *-NUMERICAL- *NUMERICAL *DTMAX 30.0 *DTMIN 1E-10 *NCUTS 8一、建立数据文件第62页，共68页。六、软件演示

29、井数据段*-WELL AND RUN DATA- *RUN *DATE 1981 05 01 *DTWELL 20 *WELL 1 2-144 VERT37 10 *PRODUCER 1 *OPERATE *MAX *STO 1000.0 *OPERATE *MIN *BHP 3000.0 *SHUTIN *MONITOR *MIN *STO 0.0 *SHUTIN *MONITOR *MAX *GOR 4000.0 *SHUTIN * rad geofac wfrac skin *GEOMETRY *K 0.089 0.37 1.0 0.0 *PERFV *GEO 1 * KF FF 3:8

30、 1 10:11 1 *SHUTIN 1一、建立数据文件第63页，共68页。六、软件演示二、历史拟合区块压力、含水拟合曲线第64页，共68页。六、软件演示二、历史拟合单井产油、含水拟合曲线第65页，共68页。六、软件演示三、预测结果分析各汇总方案的日产油量曲线 各汇总方案的含水率曲线 各汇总方案的含水上升率曲线 各汇总方案的采出程度曲线 方案10为优选方案第66页，共68页。课程结束,谢谢各位！第67页，共68页。1、不是井里没有水，而是你挖的不够深。不是成功来得慢，而是你努力的不够多。2、孤单一人的时间使自己变得优秀，给来的人一个惊喜，也给自己一个好的交代。3、命运给你一个比别人低的起点是想告诉你，让你用你的一生去奋斗出一个绝地反击的故事，所以有什么理由不努力!4、心中没有过分的贪求，自然苦就少。口里不说多余的话，自然祸就少。腹内的食物能减少，自然病就少。思绪中没有过分欲，自然忧就少。大悲是无泪的，同样大悟无言。缘来尽量要惜，缘尽就放。人生本来就空，对人家笑笑，对自己笑笑，笑着看天下，看日出日落，花谢花开，岂不自在，哪里