数值模拟培训讲义.docx
《数值模拟培训讲义.docx》由会员分享,可在线阅读,更多相关《数值模拟培训讲义.docx(26页珍藏版)》请在冰豆网上搜索。
数值模拟培训讲义
数值模拟培训讲义
----Eclipse软件应用部分
第一部分:
数值模拟用数据资料准备
第二部分Eclipse简介及建模步骤
第三部分:
地质建模及前处理模块GRID的使用
第四部分:
数值模拟计算结果分析及后处理部分
第一部分数值模拟用数据资料准备
在进行数值模拟之前,需要收集一些相关的数据,以便为后面的数值模拟作准备,这些资料总体来讲可以分为两大部分,一是静态资料,二是动态资料。
为方便数值模拟操作人员更好更全的收集这些资料,下面将这些必要的资料作一总结如下:
(一)静态资料
1.小层数据表或等值线图[包括砂层厚度、有效厚度(或净毛比)、顶部深度、孔隙度、渗透率等];
2.地质储量及地层、油藏特点的总结报告;
3.油、气、水高压物性PVT数据;
4.油水、油气相渗曲线数据和毛管压力曲线数据;
5.原始地层压力、温度、压力系数数据;
6.油、气、水分布(原始饱和度)或压力分布或油水界面和油气界面;
7.井位分布图;
8.流体和岩石化验分析报告;
(二)动态资料
1.射孔完井报告;
2.井史报告、压裂等措施;
3.系统测压资料;
4.试油、试井和试采资料(压力恢复曲线);
6.油水井别,调整井位示意图;
7.油井生产(水井注水)数据报表:
日产油、日产液、日产气、综合含水、压力
累积产油(气、水、液)
日注水、累积注水
8.区块综合生产数据统计报表:
日产油(水、气、液)、采出程度、综合含水
累积产油(气、水、液)
日注水、累积注水
第二部分Eclipse简介及建模步骤
一、Eclipse简介
Eclpise是斯伦贝谢公司开发的一套数值模拟软件,它界面好,图形输出功能强大,可输出两维和三维视图,并可以进行角度变换,能够很好处理断层,并能半自动进行敏感性分析。
Eclpise不仅为各种各样的油藏和各种复杂程度(构造、地质、流体、开发方案)的油藏提供了准确、计算快速的多项选择,而且还提供了全隐式、IMPES、AIM和IMPSAT求解方法,可以在任何工作平台上运行,包括UNIX和PC等,并能够完成在多个处理器上的大型并行计算。
二、Eclipse建模步骤及基本数据
为了熟悉用Eclipse建立地质模型的基本步骤,在这一部分中以一个一维均质等厚各向同性井组为例作一介绍。
模拟井组基本情况如下:
模拟区块面积为1040m*560m,网格结点882个(21*21*2),有效厚度和顶部深度通过文件输入,两个层的X方向渗透率和Y方向渗透率分别为1.25md和2.5md,Z方向渗透率为0.1md,孔隙度分别为0.12和0.15,井网井位如下图所示:
基本步骤如下:
(一)CaseDefinition
1.SimulatorBlackOil
2.General
(1)SimulationStart:
Day1MonthJanYear1970
(2)Title:
*******
(3)ModelDimensions:
NumberofCellsinXDirection:
21
NumberofCellsinYDirection:
21
NumberofCellsinZDirection:
2
(4)Units
3.Reservoir
AquifersAnalytical
Numerical
NumberofReservoirs1
GridTypeCartesian
GeometryTypeBlackCentred
4.PVT
GasOilWaterDissolvedGas
(二)Grid(Subsection→GRIDKeywords)
1.Geometry
(1)GridAxeswrtMapCoordinates
XcoordofYaxis:
0;YcoordofYaxis:
560;
Xcoordoforigin:
0;Ycoordoforigin:
0;
XcoordofXaxis:
1040;YcoordofXaxis:
0;
(2)XGridBlockSizes:
49.5249.52
(3)YGridBlockSizes:
26.6726.67
(4)ZGridBlockSizes:
1015
(5)DepthsofTopFaces:
15201530
(6)GridDataUnits:
METRES
注:
在此可以通过读入数据文件的格式加载顶部深度、有效厚度等数据文件。
2.Properties
(1)NetThickness:
1015
(2)XPermeability:
1.252.5
(3)YPermeability:
1.252.5
(4)ZPermeability:
0.010.01
(5)Porosity:
0.120.15
(三)PVT(Section→Keywords)
1.FluidDensitiesatSurfaceConditions
(1)OilDensity:
854kg/m3
(2)WaterDensity:
1000kg/m3
(3)GasDensity:
1.0708kg/m3
2.RockProperties
(1)ReferencePressure:
350bar
(2)RockCompressibility:
0.0001351/bar
3.LiveOilPVTProperties(DissolvedGas)–PVTOkeyword
fromdatafile:
pvto.txt
4.DryGasPVTProperties(NoVapourisedOil)–PVDGkeyword
fromdatafile:
pvtg.txt
5.WaterPVTProperties
(1)ReferencePressure(Pref):
350bar
(2)WaterFVFatPref(Bw):
1.006rm3/sm3
(3)WaterCompressibility(Cw):
6.25E-51/bar
(4)WaterViscosityatPref(μo):
0.6cp
(5)WaterViscosibility:
01/bar
(四)SCAL(Section→Keywords)
(1)Water/OilSaturationFunction(SWOF):
sw.txt
(2)Gas/OilSaturationFunctionvsLiquidSaturation(SLGOF):
sl.txt
(五)Initialisation(KeywordsType→Equilibration)
1.EquilibrationDataSpecification(EQUIL)
(1)DatumDepth:
1525m
(2)PressureatDatumDepth:
166.7bar
(3)WCODepth:
1545m
(4)OWCapPressure:
168.7bar
(5)GOCDepth:
0m
(6)Rs/PbvDepthTable:
1
2.RsvDepth(RSVD)
Depth(m)Rs(sm3/sm3)
152580.2
132580.2
3.DefineAquifer→Initialisation(KeywordsType→Aquifer)→FetkovichAquiferData(AQUFET)
(1)DatumDepth(m):
1545
(2)InitialPressure(bar):
166.7
(3)InitialVolume(m3):
20000
(4)TotalCompressibility(1/bar):
0.0002
(5)ProductivityIndex(sm3/day/bar):
0.01
(6)PVTWaterPropertyTable:
1
(7)LowerIConnection:
1
(8)UpperIConnection:
21
(9)LowerJConnection:
1
(10)UpperJConnection:
21
(11)LowerKConnection:
2
(12)UpperKConnection:
2
(13)ConnectionFace:
k+
(14)SaltConcentration(kg/m3):
0
(六)Regions(空)
(七)Schedule(Events→New)
1.TimeEntryChoice(Time→Insert)
(1)TimeStep
(2)AddTimeSteptoSubsequentSteporTimes
TimeStep1yrNum10
2.DefineWells,Groups,&Connections→WellSpecification(WELSPECS)
(1)Well:
L37-37L37-39L39-37L39-39L38-38
(2)ILocation,JLocation,PreferredPhase
(3)DrainageRadius:
0.1m
3.DefineWells,Groups,&Connections→WellConnectionData(COMPDAT)
(1)Well:
L37-37L37-39L39-37L39-39L38-38
(2)ILocation,JLocation,KUpper,KLower
(3)WellBoreID:
0.1m
4.WellControlsandLimits→ProductionWellControl(WCONPROD)
(1)Well:
L37-37L37-39L39-37L39-39
(2)Control:
LRAT
(3)LiquidRate:
1.5m3/d
5.WellControlsandLimits→InjectionWellControl(WCONINJE)
(1)Well:
L38-38
(2)InjectorType:
WATER
(3)LiquidSurfaceRate:
6m3/d
(八)Summery(选择所要输出的内容)
(九)MultipleSensitivity(空)
(十)ECLIPSEOffice→RUN
第三部分:
地质建模及前处理模块GRID的使用
静态地质模型主要包括:
油层顶深、油层x、y方向的网格划分、各网格的孔隙度、渗透率数据。
建立静态地质模型的基本数据:
等值线图(油层顶深,油层厚度、孔隙度、渗透率)或者各井点的上述参数值。
用Grid建模型的主要工作是将上述数据输入计算机,有Grid输出各网格点的上述相关数据。
EditModel
a.输入Model名称
b.选择坐标系
c.选择单位
d.输入基本信息
EditMap
输入map名
EditWells
EditFaults
EditContours
EditGrid
输入Grid名称
b.输入Grid的等值线图
c.选择网格划分方法,划分网格
d.定义顶层深度,各层厚度,添加层
e.添加属性,编辑各层的属性,输出结果
f.输出网格模型
第四部分:
数值模拟计算结果分析及后处理部分
一、数据文件调试(演示)
A、用E100直接算时的查错方法介绍
错误:
如:
PVT表中压力项非单调增等;
B、用office算时的查错方法介绍
二、历史拟合
1、SCHEDULE部分需要作如下的修改:
生产井不用WCONPROD,而改用历史拟合用的关键字:
WCONHIST,注入井仍用WCONINJ。
WCONPROD与WCONHIST的区别在于:
WCONPROD主要用于历史拟合结束后进行动态预测,其所定的产量或定的压力在以后的生产中其主要的控制作用,井的生产制度不会超过这一关键字所规定的生产制度;而WCONHIST所输入的产量为油田实际的统计数据,这一产量仅仅作为一参考值,模拟过程中的产出量仅与所建立的模型的物性参数有关,模拟的产出量可比WCONHIST中所规定的值大也可比它小,因此需要修改模型相关的参数以达到两者相近的目的,也即达到历史拟合的目的。
这两种定义方式如下:
WCONHIST
'E3-6''OPEN''LRAT'20.040.00.04*/
/
WCONPROD
‘E3-6''OPEN''LRAT'3*59.285*/
/
2、历史拟合内容:
区块:
储量拟合、压力拟合、累积产液量、累积产油量、累积产水量、累积注入量、日产油量、日产水量、含水率、汽油比等的拟合
单井:
累积产液量、累积产油量、累积产水量等的拟合
3、拟合方法简介
储量的拟合:
由储量计算的经验公式:
可知其主要影响因素有:
孔隙度、油层厚度、原始含油饱和度,对于饱和度而言主要受以下两方面的影响:
●饱和度的变化来说,水压缩系数的改变对饱和度变化有较明显的影响
●油水界面的高低可以直接影响原始含油饱和度,从而影响储量
压力的拟合:
可调参数:
孔隙度、厚度、饱和度、油层综合压缩系数、渗透率、粘度、水体、注水量的分配等,但在储量已拟合的前提下,孔隙度,厚度、饱和度的修改要慎重考虑。
在给定产量的条件下,增大孔隙度或厚度,可使计算压力值升高,这两个参数的改动都会造成地质储量的改变,所以在调整这些参数时都要慎重考虑这种调整的合理性。
压力剖面呈不合理的不光滑形状,这种情况常常可能是由于该处的渗透率值过低所引起。
把该处的渗透率值乘以一个大于1的常数,即把渗透率普遍提高一个幅度,增加了流动性,就可以使压力剖面变成比较光滑的曲线。
所以水体的大小和边外渗透率的高低常常只是一个大致的估计值,可靠性较差,所以有关边外水体的参数是拟合压力时需要考虑的一个重要因素。
切割注水时,注水井排两侧区块的注入水量的分配比例应该随着这些区块的地质条件和开发历史的差异而有所不同。
原油粘度越大,井底压力也越大。
原油粘度越大,两相渗流的范围越大,表明油的饱和度值变化得越缓慢。
无论是油或水的压缩系数值的变化都对压力差值及其变化产生明显的影响。
而对饱和度的变化来说,水压缩系数的改变对饱和度变化有较明显的影响。
单井动态的拟合难度比较大,常对表皮效应进行调整、有时也调整井周围各网格的渗透率值。
含水率或气油比的拟合
含水率和气油比的拟合都主要依靠相对渗透串曲线的修改,前者需得到水的相对渗透率曲线,后者则涉及气的相对渗透率曲线,两者原理和方法基本相同。
当计算的含水率高于实测值时,应把水的相对渗透率曲线下移,反之,则应上移对于气油比的拟合,也同样处理,
至于计算见水时间的过早或过迟,主要和水相渗透率曲线的端点位置即临界饱和度的大小有关:
当计算见水时间过早,则应把水相临界饱和度值增大,即把此端点右移,反之,则应将其左移;
当计算的气窜时间过早或过迟时,处理方法也类似,只是由于油气相对渗透率曲线图上的横坐标为So,所以若计算气窜时间过早,虽然拟合时同样是增大临界饱和度值,但却是把端点左移;反之亦然。
当调整相对渗透率曲线时,由于该相流量也随之而发生相应的改变,必然导致相应的压力变化。
例如,当把水的相对渗透率曲线下移时,由于水相的流量减少,必然导致流动时的压力值减少,这也是在拟合含水率数值时应考虑的一个问题。
影响含水率及气油比计算值的因素还有油水界面和气水界面的位置。
例如当计算时输入的油水界面高于实际值时也会造成见水过早及含水上升过快。
所以在拟合时还应检查所给的油水界面或气水界面的位置是否准确,发现问题应作适当调整。
上面仅仅对历史拟合中的常见影响因素做了一点分析,其中还有很多的影响因素需要在模拟的过程中不断总结,不断探索。
总之,历史拟合是一个非常费时,而且非常难做的一部分工作,它不仅需要认真仔细有耐心,而且需要操作人员有较好的专业知识,这样才能对参数的调整作出较为准确的判断。
三、结果显示(演示)
1.需要对数据文件中的Summary部分作相应的修改,以便输出我们需要的模拟结果。
1.输出文本文件或者数据报表,用REPORT模块能很方便输出所需要的报表文件。
3.用Summary模块输出曲线图,平面图,三维图等,而且在这一模块中可以对不同时间步下的模拟结果进行动态显示,这让数值模拟人员能够非常形象的对结果进行观察。
四、结果分析
对于结果的分析主要从以下几方面来分析:
●采出程度对比(或累积采油量)
●含水率对比
●地层平均压力对比
●剩余油饱和度、储量等分布图
●年产油量对比
具体可以参见多媒体报告中的:
“后处理部分.PPT”文件。
五、后处理模块-GRAF介绍
Themainoptionsare:
0PRIMARYMENU
1ReadECLIPSEoruserdata
2Printorindexloadedinformation
3Createapicture
4Modifyapicture
5Displayapicture
6Loadorsaveworkspace
7Copyonepicturetoanother
8Writeorexecutegraphicsrunfile
9ModifyGranitesettings
10Endsession
11Additionalfacilities
-------------------------------------------------------------------------------------------------------------------------
0主菜单
1读取ECLIPSE文件或用户数据
2打印或搜索图片或工作空间的信息
3作图
4修改图
5显示图
6读取或保存工作空间
7将一张图片复制到另一图片中
8生成或执行GRF文件
9修改线型,字体,颜色等的设置
10结束
11附加工具
-------------------------------------------------------------------------------------------------------------------------
GRAFcanread,processandstorethefollowingtypesofdata:
Vector:
Asetofnumberswhichcanbeusedforconstructinglinegraphs,e.g.fieldoilproduction(FOPR)ateachtimestep.VectordatamaybereadfromECLIPSEsummaryfiles(option1.1)oruserdatafiles(option1.2).TheECLIPSEinitialfilealsocontainssomevectordata,e.g.forrelativepermabilitiesvs.saturations.
Field:
Thisdefinesthegeometryofareservoirsimulationgrid,asasetofcornerpoints.FielddataarereadfromECLIPSEgridgeometryfilesusingoption1.3.
Solution:
Thisdefinesasetofsolutionvariables(pressures,saturationsetc.)atareportstep.SolutiondataarereadfromECLIPSErestartfiles(option1.4)oruserdatafiles(option1.2).Gridpropertiessuchasporositymayalsobeloadedassolutiondatafromaninitialfileusingoption1.5.
Table:
ThisdefinesasetofdataforaVFPtableasusedforinputtoECLIPSE(option1.8).
RFT:
Vectordatadefiningdepths,pressuresandsaturationsatwellconnections(option1.9).
-------------------------------------------------------------------------------------------------------------------------
1LOADINFORMATIONTOBEDISPLAYED
0Returntoprimarymenu
1ReadfromECLIPSEsummaryfile
2Readfromuserdatafile
3ReadfromECLIPSEgridfile
4ReadfromECLIPSErestartfile
5ReadfromECLIPSEinitialfile
6Changetoformattedinput
7Changetounifiedfiles
8ReadfromVFPfile(alwaysformatted)
9ReadfromEclipseRFT