基于有限差分的油水两相渗流方程求解油藏数值模拟毕业设计.docx
- 文档编号:812407
- 上传时间:2022-10-13
- 格式:DOCX
- 页数:34
- 大小:490.87KB
基于有限差分的油水两相渗流方程求解油藏数值模拟毕业设计.docx
《基于有限差分的油水两相渗流方程求解油藏数值模拟毕业设计.docx》由会员分享,可在线阅读,更多相关《基于有限差分的油水两相渗流方程求解油藏数值模拟毕业设计.docx(34页珍藏版)》请在冰豆网上搜索。
基于有限差分的油水两相渗流方程求解油藏数值模拟毕业设计
本科毕业设计(论文)
基于有限差分的油水两相渗流方程求解
学生姓名:
学号:
专业班级:
指导教师:
2012年06月20日
摘要
为了保证油藏的稳定产油量以及最终采收率,以获得最大的经济效益,必须对油藏的压力以及饱和度等参数进行监控,因而需要对油藏进行模拟,以确定合适的开采时间、注水量、开采速度等问题,因而油藏数值模拟对于提高油田效益是至关紧要的。
油水两相渗流方程描述了原油开发过程中的油水两相流动过程,是一对耦合的偏微分方程,针对这组方程采用有限差分方法实现其数值求解,以对油藏的压力、饱和度进行即时监控。
本文应用隐式压力——显式饱和度解法(IMPES方法)来对油藏油水两相一维渗流进行模型建立与求解,这一方法来自于Sheldon等人,以及Stone和Gardner的著作,基本思想是合并流体方程得到一个只含有压力的方程,某一时间步的压力求解出来后,饱和度采用显式更新。
经过Matlab编写程序,实现对一维油水两相渗流方程的数值求解。
矚慫润厲钐瘗睞枥庑赖。
关键词:
油水两相;有限差分;压力;含水饱和度;油藏模拟
ABSTRACT
Displacingoilbywaterisalwaysusedintheexploitationoftheoilreservoir.Inordertoobtainthemaximumeconomicbenefits,theparametersofreservoirpressureandsaturationshouldbemonitoredtoensurethestabilityofoilproductionandtheultimaterecoveryofthereservoir.Thusthereservoirmustbesimulatedsoastodeterminetheappropriaterecoverytime,waterinjectionrateandrecoveryrate.Itisnecessarytoestablishtheoil-watertwo-phaseflowmodelfortheoil-watertwo-phaseflowequation,coupledpartialdifferentialequations,candescribetheoil-watertwo-phaseflowintheprocessofcrudeoildevelopment.Achievingthenumericalsolutionaccordingtothisgroupofequationsandusingthefinite-differencemethodcandorealtimemonitoringforreservoirpressureandsaturationsothatappropriatemeasurescouldbetakentoensuretheminingofoilproductioninhighefficiencyandsustainable.Inthisthesis,implicitpressure-explicitsaturationsolution(IMPESmethod)whichisinventedbySheldon,etc.andfromworksofStone’sandGardner’sisappliedforthebuildingandsolvingofreservoir’sone-dimensionalwater-oiltwo-phaseseepagemodel.TThebasicideaiscombiningthefluidequationstoachieveaequationwhichcontainsonlythepressure.Aftersolvingthepressureattheverytimeandstep,thesaturationwillbedisplayedupdated.WritingprogramsthroughMatlab,one-dimensionalwater-oiltwo-phaseseepageequationswouldbesolved聞創沟燴鐺險爱氇谴净。
Keywords:
water-oiltwophase;finitedifference;pressure;watersaturation;reservoirsimulation残骛楼諍锩瀨濟溆塹籟。
第1章前言
1.1油藏数值模拟在油田开发中的重要意义
油气藏是在单一圈闭中具有统一压力系统的油气聚集单元。
在原始条件下,油气藏处于平衡状态;当受到干扰(如打井、生产)时,原来的不平衡状态被打破,油气藏处于动态变化中。
油气藏从投入生产到最后废弃就是一个不断变化的动态过程。
猫虿驢绘燈鮒诛髅貺庑。
描述或实现油气藏动态变化的过程叫做模拟,油藏数值模拟是用数值的方法来求解描述油藏中流体渗流特征的数学模型,是一门将计算机、应用数学、油藏工程等结合起来的综合性工程应用学科,在油气田开发方案设计和动态分析中有十分重要的作用,为了保证油藏的稳定产油量以及最终采收率,以获得最大的经济效益,必须对油藏的压力以及饱和度等参数进行监控,因而需要对油藏进行模拟,以确定合适的开采时间、注水量、开采速度等问题,需建立油水两相渗流模型,以对油藏的压力、饱和度进行即时监控,及时采取合适的开采措施来保障油田的高效、高产、持久生产,它能够从油田实际出发对油藏进行精确地描述,以最少的投资、最合适的速度去获得最高的最终采收率,也就是要获得最大的效益。
锹籁饗迳琐筆襖鸥娅薔。
1.2国内外研究方法现状
油藏模型中有压力的、和饱和度、两组未知量,目前基本上有两类求解方法。
一类是顺序求解法,即先求压力项,后求饱和度;另一类是联立求解法,即同时求解压力项和饱和度项。
同时,由于方程中含非线性系数,他们依赖于压力和饱和度的变化,在求解数学模型时有如下几种处理方法:
構氽頑黉碩饨荠龈话骛。
(1)隐式压力显式饱和度方法(IMPES)
即隐式求解压力方程,显式求解饱和度方程。
它属于顺序求解法的一种,是数值模拟中最常用的、最简单的一种方法。
具有占内存小、计算工作量小、方法简便等优点。
但该方法存在两个问题:
第一,达西项的系数处理是显式的,因此对如锥进的问题,由于井底周围流速高,压差变化大,而存在较大误差,对于强非线性问题的适应性也差;第二,饱和度的计算是显式的,当时间步长较大时,会出现解的不稳定性。
因此,IMPES方法只适用于一般的弱非线性渗流问题,对于某些非线性问题如注气、气锥或水锥等问题,IMPES方法无能为力,即使时间步长取得很小,仍会出现解的振荡或算出的压力和饱和度为负值的情况,以至于模拟计算无法正常进行。
輒峄陽檉簖疖網儂號泶。
(2)半隐式方法
半隐式方法属于联立求解的方法的一种,也是数值模拟中常用的一种方法。
其基本思路是:
联立求解油相方程和水相方程,同时求出压力和饱和度,因此压力和饱和度都是隐式求解。
在计算过程中,半隐式方法对方程右端项的处理与IMPES方法完全相同,不同之处在于对方程左端项的处理。
它需要对方程左端的达西系数、产量项及毛管压力进行泰勒级数展开,并忽略二阶小量,一阶导数项用n时刻的值。
由于系数是处理近似的,并未真正用n+1时刻的值,所以这种方法叫做半隐式方法。
尧侧閆繭絳闕绚勵蜆贅。
(3)隐式压力隐式饱和度法(IMPIMS)
隐式压力隐式饱和度法实际上是IMPES方法和半隐式方法的混合变种,它也是属于顺序求解方法的一种。
它既有半隐式方法求解饱和度的特点,又保留了IMPES方法省内存、省工作量的特点。
识饒鎂錕缢灩筧嚌俨淒。
IMPIMS方法(TheImplicitPressureImplicitSaturationMethod)实际上是IMPIMS方法和半隐式方法的混合和变种,它也是属于顺序求解方法的一种。
它既有半隐式方法求解饱和度的特点,又保留了IMPES方法省内存、省工作量的特点。
凍鈹鋨劳臘锴痫婦胫籴。
IMPIMS方法的求解思路是:
压力和饱和度分开顺序求解,求解压力时可直接利用IMPES方法的压力求解方法,然后将求出的压力值代入半隐式方法的水相差分方程中,将该方程化为只有含水饱和度的一个变量的差分方程,再用隐式计算格式求解即可。
恥諤銪灭萦欢煬鞏鹜錦。
(4)全隐式方法
对于某些强非线性渗流问题如高速渗流等,即使使用半隐式方法也会引起计算结果的波动,或者时间步长只能取到很小。
为此提出了全隐式方法,这里简单介绍一种与半隐式方法相类似的全隐式方法,即达西系数项也用泰勒级数展开并忽略二阶小量,但一阶导数项不用时刻的值而用时刻的值,由此构成一个非线性的代数方程组,可用牛顿迭代法或其他非线性代数方程组的解法进行求解。
鯊腎鑰诎褳鉀沩懼統庫。
1.3本文章节安排
论文根据质量守恒原理建立了一维油水两相渗流的数学模型,经过逐步简化得到实验室进行单管模型一维水驱油实验时的数学模型,利用IMPES方法对此模型进行处理。
第一部分首先列述了油藏渗流力学的相关知识,对油藏的油水渗流机理进行了详细的介绍以及油水两相渗流方程的推导过程,然后对有限差分法的基础知识进行了详细的介绍;第二部分主要是油水两相渗流模型的建立过程以及应用此方法对实例油藏的模拟。
最后综述了油水两相渗流方程的建立过程中所遇到的问题、模型建立的理想条件及其准确性、针对实例模拟得出的结论。
硕癘鄴颃诌攆檸攜驤蔹。
第2章相关基础知识
2.1油藏渗流力学相关知识
(一)、渗流过程中的力学分析
油气水之所以能在岩石空隙中渗流,是各种力共同作用的结果。
1、流体的重力和重力势能
流体具有质量,在中立场中受到地球的吸引力而具有重力。
在渗流过程中,流体的重力和它的相对位置联系起来,就表现为重力势能。
这种重力势能用压力表示为:
阌擻輳嬪諫迁择楨秘騖。
(2-1)
式中---表示重力势能的压力,Pa
---流体密度,
---重力加速度,
---相对位置高差,
重力有时是动力,有时为阻力。
2、流体的质量和惯性力
惯性是物体所固有的一种物理特征,其大小取决于质量。
流体质量的大小一般用密度表示。
流体由于具有质量,因此也具有惯性。
当流体运动时,惯性使其总要维持原状,因而惯性力在渗流过程中多表现为阻力。
氬嚕躑竄贸恳彈瀘颔澩。
根据牛顿第二定律,惯性力用压力表示为:
(2-2)
式中,---惯性力
---质量
v---速度
---加速度
由于渗流时渗流速度通常很小,因此常忽略惯性力。
3、流体的粘度及粘滞力
粘滞力是流体的一种特殊属性。
在流动的流体中,如果各层的流速不同,将有一对作用力及反作用力,使原来快的流层减速,而慢的加速。
这一对等值而又反向的力阻碍着流层的相对运动。
流体的这种属性叫粘滞力。
度量粘滞性大小的参数叫粘度。
由牛顿内摩擦定律表示为:
釷鹆資贏車贖孙滅獅赘。
(2-3)
式中,A---两流层的接触面积,
---沿流层法线方向的流速梯度,
F---内摩擦力,N
---粘滞系数,粘度的单位通常用表示:
1
粘度的单位以表示时成为泊(P):
1P=100cP(厘泊)
在渗流中,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 有限 油水 两相 渗流 方程 求解 油藏 数值 模拟 毕业设计