油藏数值模拟原理

Click to edit Master title style,Click to edit Master text styles,Second level,Third level,Fourth level,Fifth level,-,*,-,油藏数值模拟原理,黑油模型,1,黑油模型简介,黑油模型是指油质较重性质的油藏类型。,黑油模型是最完善、最成熟，也是应用最为广泛的模型。,是油藏数值模拟的基础，其它模型大都是黑油模型的扩展。,目的：软件这样一个“黑箱”是如何操作的，以便指导我们更好地去应用软件。,2,黑油模型简介,黑油模型的基本假设,单相渗流黑油模型基本微分方程,三维三相黑油模型基本微分方程,初边值条件,5.,网格系统,6.,黑油模型的差分方程,7.,差分方程的线性化,8.,线性代数方程组的求解,9.,黑油模型的主要数据流,3,1,、黑油模型基本假设,（1）油藏中的渗流是等温渗流。,（2） 油藏中最多只有油、气、水三相，每一相均遵守达西定律。,（3）油藏烃类只含有油、气两个组分。在油藏状态下，油气两组分可能形成油气两相，油组分完全存在于油相内，气组分则可以以自由气的方式存在于气相中，也可以以溶解气的方式存在于油相中，所以地层内油相为油组分和气组分的某种组合。在常规油田中，一般不考虑油组分向气组分挥发的现象。,（4）油藏中气体的溶解和逸出是瞬间完成的，即认为油藏中油气两相瞬时达到相平衡状态。,（5）油水之间不互溶；天然气也假定不溶于水。,4,2.,单相渗流,基本微分方程,单相渗流黑油模型基本微分方程,黑油模型基本微分方程主要是由,达西公式,(,达西定,律,),和,连续性方程,(,物质守恒定律,),所构成,一维单相渗流的,Darcy,公式为：,三维单相渗流的,Darcy,公式为：,5,渗流速度是一个空间向量，用分量表示为：,渗流速度矢量； 压力； 势函数, 绝对渗透率； 重力加速度常数；,粘度； 垂向坐标； 密度；,2.,单相渗流,基本微分方程,6,Reservoir Simulation III,Gringarten, 2002,Data from all sources,7,2.,单相渗流,基本微分方程,根据物质守恒原理，在单位时间内,流入单元内的流体质量-流出单元的流体质量,=单元内流体质量的变化,取渗流场中一个微小六面体体积单元来研究,8,9,2.,单相渗流,基本微分方程,该六面体单元的中心点坐标为（,x，y，z），,其,每一侧面的质量流速均以其侧面的中心点的质量,来代替。,方向流入质量=,方向流出质量=,方向流入质量=,方向流出质量=,方向流入质量=,方向流出质量=,单元体内流体质量的变化=,10,根据,物质守恒定律,：,将上式两边同除以 ，,2.,单相渗流,基本微分方程,11,并令 ， ，,把源汇项（井）加入到方程中，得到,单相渗流,连续性方程,：,2.,单相渗流,基本微分方程,12,用,Hamilton,算子 表示，并将达西公式代入，得,到，,单相渗流微分方程（流动方程）：,生产井,注入井,关 井,2.,单相渗流,基本微分方程,流动项,源汇项,累积项,13,这就是地下流体所遵循的渗流规律。,总之，渗流基本微分方程包含三部分内容,:,流动项,、,源汇项,以及,累积项,渗流基本微分方程由两个定律控制：,达西定律与物质守恒定律,压力场和饱和度场的分布是流体流动的结果，流,动规律由达西定律来确定，流量由物质守恒来确,定，这两个规律是数值模拟的,关键,所在。,2.,单相渗流,基本微分方程,14,将单相渗流达西定律推广到多相渗流中，得到,三相流动达西公式：,式中，,w，o，g,水、油、气诸相下标。, 相对渗透率；,3. 三维三相渗流基本微分方程,15,对油、水、气三相，每一相运用物质守恒定律，可以建立以下,连续性方程,3. 三维三相渗流基本微分方程,16,把达西公式代入连续性方程，得到,三维三相渗流基本微分方程,3. 三维三相渗流基本微分方程,相对渗透率、毛管力是饱和度函数,而粘度、体积系数孔隙度与密度均是压力的函数,17,3. 三维三相渗流基本微分方程,令,18,上述,T,的物理意义是反映相邻网格间的流动性，又叫传导率（,transmissibility),3. 三维三相渗流基本微分方程,19,3. 三维三相渗流基本微分方程,求解变量,20,3. 三维三相渗流基本微分方程,还要补充一些方程,饱和度 ：,油水毛管力：,油气毛管力：,相对渗透率和毛管力是饱和度函数：,21,Stone,公式：,式中， ,束缚水饱和度,油水两相流动时油相相对渗透率,油气两相流动时油相相对渗透率,3. 三维三相渗流基本微分方程,22,3. 三维三相渗流基本微分方程,流体,PVT,性质,水相：,气相： （,Z),油相：随压力变化，相态会发生变化。,三相状态：,两相状态：,岩石的高压物性,相渗曲线以及流体,PVT,性质，以表格形式给出，,称为静态数据。,23,3. 三维三相渗流基本微分方程,实际求解变量,初始化输入数据,ljfinit.txt,油相流体性质与压力和饱和压力的关系,24,产量项（作源汇项处理）（动态数据）,3. 三维三相渗流基本微分方程,25,3. 三维三相渗流基本微分方程,其中，, 在井基准面深度 上的井底压力。,节点等效泄油半径。, 井筒半径。,X，Y，Z,网格节点,x,，,y,和,z,方向步长,PI,叫采油指数， 通常写成,KH,可以直接输入,井数据又叫动态数据,ljfmodel.txt,26,4,.初边值条件,油气藏外边界条件,:,第一类边界 条件,P=f(x,y,z,t),第二类边界 条件,第三类边界条件,:,27,4,.初边值条件,外边界,封闭边界,:,尖灭、断层、圈闭,流动边界：边、底水,内边界，即,井的生产条件,定产边界条件（,rate constraint):,给定产液量、产油量、产气量、注水量、注气量,定压边界条件,(pressure constraint):,给定井底流压,28,4,.初边值条件,油藏初始条件,:,29,5.,网格系统,30,在一般情况下,流动方程,求不出解析解，只能用离散化方法求数值解。,离散化方法，首先要把求解区域按一定的网格系统进行剖分。网格剖分方法有很多，像直角坐标网格、径向网格、角点网格、,PEBI,网格等。目前，商用软件中大都采用,矩形网格系统,。离散化方法主要有：有限差分,（Finite Difference)、,有限元(,Finite Element),以及流线法(,Stream Line)。,比较成熟的是,有限差分方法。,有限差分方法大多空间变量采用,网格块中心（,Block Center),差分, 时间变量按,向前差分方法,。,5.,网格系统,31,5.,网格系统,:,矩形网格,32,33,5.,网格系统,:,径向网格系统,读取网格数据时，,R,方向循环最快，然后是 方向，最后是,Z,方向,径向网格数据读取规则,34,5.,网格系统,:,角点网格,Corner Point,35,块中心网格与角点网格的比较,块中心网格表示的含有断层的倾斜构造,角点网格表示的含有断层的倾斜构造,36,5.,网格系统,:,块中心网格的数据结构,模型被分为,20510,个网格区域。,网格从（,1,，,1,，,1,）开始向下倾斜。,各网格块在,X,方向上是,300ft,，在,Y,方向上是,1000ft,。,层厚度从上至下分别是,32,，,22,，,20,，,4,，,32,，,4,，,26,，,26,，,4,，,28ft,37,5.,网格系统,:,块中心网格的数据结构,-,前,20,个,TOPS,定义网格,(1, 1, 1),到,(20, 1, 1),6855.000 6865.0006875.0006885.0006895.000,6905.0006915.0006925.0006935.0006945.000,7005.0007015.0007025.0007035.0007045.000,7055.0007065.0007075.0007085.0007095.000,-,下面,20,个,TOPS,定义网格,(1, 2, 1),到,(20, 2, 1),6930.0006940.0006950.0006960.0006970.000,6980.0006990.0007000.0007010.0007020.000,7080.0007090.0007100.0007110.0007120.000,7130.0007140.0007150.0007160.0007170.000,-,下面,20,个,TOPS,定义网格,(1, 3, 1),到,(20, 3, 1),7030.0007040.0007050.0007060.0007070.000,7080.0007090.0007100.0007110.0007120.000,7180.0007190.0007200.0007210.0007220.000,7230.0007240.0007250.0007260.0007270.000,-,下面,20,个,TOPS,定义网格,(1, 4, 1),到,(20, 4, 1),7130.0007140.0007150.0007160.0007170.000,7180.0007190.0007200.0007210.0007220.000,7280.0007290.0007300.0007310.0007320.000,7330.0007340.0007350.0007360.0007370.000,-,下面,20,个,TOPS,定义网格,(1, 5, 1),到,(20, 5, 1),7205.0007215.0007225.0007235.0007245.000,7255.0007265.0007275.0007285.0007295.000,7355.0007365.0007375.0007385.0007395.000,7405.0007415.0007425.0007435.0007445.000,/,这完成了第一层,TOPS,的定义,38,角点网格的数据结构,-COORD,-,此关键字定义坐标线,-X1Y1Z1X2Y2Z2,0.0.6825.0000.0.7023.000,300.0000 0.6835.000300.00000.7033.000,600.0000 0.6845.000600.00000.7043.000,900.0000 0.6855.000900.00000.7053.000,1200.000 0.6865.0001200.0000.7063.000,.,/,对于,20 * 10,的网格模型,需要,21*11,条坐标线，即,231,条。,-,每一条坐标线都用,6,个数来定义，所以关键字,COORD,包含有,1386,个数。,ZCORN,-,这个关键字定义网格拐点的深度，数据按照,X (,或,R),方向循环最快，然后是,Y (,或,THETA),方向，最后是,Z,方向的顺序排列。,6825.0006835.0006835.0006845.0006845.000,6855.0006855.0006865.0006865.0006875.000,6875.0006885.0006885.0006895.0006895.000,6905.0006905.0006915.0006915.0006925.000,6975.0006985.0006985.0006995.0006995.000,7005.0007005.0007015.0007015.0007025.000,7025.0007035.0007035.0007045.0007045.000,7055.0007055.0007065.0007065.0007075.000,.,/,对于,1000,个网格，需要,8000,个,ZCORN,值,-,这些只是角点网格模型的数据体中很少的一部分,39,块中心网格与角点网格的比较,BC,（块中心网格）,CP,（角点网格）,网格描述简单,网格描述复杂,前处理软件不是必须的,必须使用前处理软件,与其它许多模拟器兼容,与很少的其它模拟器兼容,很难为非常规构造建立模型,能够准确的为非常规构造建立模型,网格数据量小,网格数据量大,无法区分倾斜面和断层,区分倾斜面和断层,很难如实的建立尖灭和剥蚀面的模型,如实的建立尖灭和剥蚀面的模型,沿着断层面建立的网格连接不正确，需要手工进行修改。,能够正确的模拟沿着断层面层之间的接触关系,容易建立径向网格模型,若没有前处理软件，则很难建立径向网格模型。,40,正常连接的网格关系,粗细网格交界处的,流量计算,非正常连接技术,两点流量公式：,41,非结构网格技术,传统的结构网格,-,即规则矩形剖分的网格,具有如下不足,:,存在较严重的网格取向效应,不能保证每口井都在网格中心,不能很好地描述油藏形态,粗细网格间的过渡衔接性差,研究新的网格技术势在必行,.,结构网格的不足,x,y,好,差,井,好,开发软件,42,非结构网格技术,Gurpinar, 2001,Castellini, 2001,Wolfsteiner et al., 2002,Prevost 2003,43,非结构网格技术,(Unstructured Grids),常用的非结构网格,各网格中心之间连线构成三角网,PEBI,网格（又叫垂直平分网格,Perpendicular Bisectors):,各三角形边的垂直平分线构成网格边界,CVFE,网格,-,控制体有限元网格,(Control volume finite element):,各三角形边的中点与重心的连线构成网格边界,待解决的问题,非结构网格条件下代数方程组的快速求解,CVFE,网格,PEBI,网格,开发软件,44,5.,网格系统,: PEBI,网格,PEBI,45,unstructured PEBI grids,Locally different grids,46,New Faulted Grid,Advanced Fault Modeling,Heavily Faulted Reservoir,New Fault Traces,47,非结构网格技术,Channel,s can be modelled,48,2.,非结构网格技术,(complex) faults,Geological Model of a Heavily Faulted Reservoir,Simulation Model,49,1.,概述,局部网格加密,(Local Grid Refinement,，简称,LGR),技术主要分：,静态局部网格加密,(Static Local Grid Refinement),局部加密区域在整个计算过程中固定不变的，常用于断层、尖灭、裂缝以及井附近的加密。该技术已经成熟，在多加软件公司的软件产品中都有这一功能。,动态局部网格加密,(Dynamic Local Grid Refinement),局部加密区域在整个计算过程中随流体的运动位置而变化，主要用于追踪驱替前缘。该方法在研究之中，最早是,Heinemann,将动态局部网格用于水驱油藏数值模拟 。,局部网格加密,(Local Grid Refinement),50,Productive Area 1,Productive Area 2,Sealing Interbedding,Aquifer,Aquifer,局部网格加密技术,51,局部网格加密技术,Aquifer,Productive Area,Vertical Refinement Zone,52,随驱替前缘移动的动态局部网格加密系统,化学驱驱替前缘动态追踪数值模拟研究,53,1.3,动态局部网格加密技术,(,注化学剂时的,So,）,开发软件,54,Single Well Block,杂交网格技术,(,Hybrid,Grid Refinement),Horizontal Well Grid,Computed Results,55,开窗技术,(Windowing Technique),Independent time stepping,At any location,Several different gridding options are available,within one model in the same position,56,Window Grid,Basic Grid,Window polygon,开窗技术,57,开窗技术,Automatically and implicitly calculated,in/out-flux will be used as a von Neumann boundary condition,58,开窗技术,Good match required windowing,59,差分方程,6.,差分方程,(Difference Equations),60,6.,差分方程,61,该差分方程组无论是左端还是右端都是非线性的，是一组,非线性代数方程组,。非线性最强的参数是相对渗透率和毛管力。这两个参数在时间方面的处理不同，对计算的稳定性以及计算结果影响最大，也是形成形形色色解法的关键所在。,根据隐式程度的不同，数值模拟的算法有,IMPES (Implicit Pressure Explicit Saturation),方法、,半隐式(,Semi-Implicit Method),方法,和,全隐式方法(,Fully Implicit Method ),。,IMPES,：,隐式求解压力显式求解饱和度。对流度中与时间有关的量均取上一时间步的值。,6.,差分方程,62,IMPES,方法的优点是速度快，缺点是稳定性差,尤其是锥进问题以及过泡点问题不能很好解决。因此时间步长不能太大。目前，,IMPES,方法在化学驱模型，尤其是流线法中发挥着作用。,全隐式方法,主要采用了,Newton-Raphson,迭代技术, 对方程中各项系数，包括井系数全部隐式处理。全隐式方法较之半隐式方法，具有更好的稳定性，时间步长更长，彻底解决了过泡点问题，,是八十年代数值模拟技术走向成熟的关键技术，也是大型油藏数值模拟的主流方法。,6.,差分方程,63,全隐式方法和半隐式方法的根本差别在于：,半隐式方法,是将非线性方程直接线性化形成,线性代数方程组，来求解,n+1,时刻的值。,全隐式方法,采用,Newton-Raphson,方法，迭,代地求解非线性方程组，得到,n+1,时刻的变量。,用一维非线性方程表示： 给定初始近似,对函数 在点处 进行一阶泰勒展开, 可近,似地表示为,6.,差分方程,64,半隐式：,全隐式方法：,6.,差分方程,65,下面介绍全隐式方法。为了减小计算机舍入误,差的影响， 模型中不直接求压力和饱和度，而是求其增量。,其中：,x,代表求解变量,n,时间迭代步,牛顿迭代步,6.,差分方程,66,求解变量： ， ，,对于二相状态：,对于三相状态：,6.,差分方程,67,对方程组进行线性化处理。对所有变量进行,Taylor,展开，保留线性部分。变量进行运算后也保留线性部分，舍去高阶项。形成线性代数方程组。,7.,差分方程的线性化,68,全隐式方法形成的线性代数方程组的结构,7.,差分方程的线性化,x,b,A,矩阵结构示意图,69,油藏数值模拟的线性代数方程组是大型稀疏矩阵，其求解效率在油藏数值模拟计算中举足轻重。因为，其计算量占整个计算量约70%80%。因此，研制快速求解大型稀疏矩阵的新方法成为油藏数值模拟发展的一个重要方向。,线性代数方程组的求解方法主要有三类：,1）直接法：,Gauss,消元、,D4,排序的,Gauss,消元,2）迭代法：,SOR(,松弛迭代法）,3）共轭梯度类法、预处理共轭梯度类法（,PCG,方法）,8.,线性代数方程组的求解,70,目前最先进的方法是预处理共轭梯度类方法。,BO,软件中采用,RSVP,方法。该方法计算快，精度高。,8.,线性代数方程组的求解,71,黑油模型油藏数值模拟器流程,二、黑油模型简介,72,8.,黑油模型数据流,初始程序,INIT,文件,输出结果,动态数据,动态程序,主模型程序,输出结果,输出,重启文件,输入重启文件,拷贝文件,静态数据,73,黑油模型输入数据：,主要数据流,8.,黑油模型数据流,74,一、初始化数据主要有：,1,PVT,数据包括流体性质常数，高压物性试验数据和平衡区数据。主要是：,（1）,油、气、水的地面密度和比重等常数，以及计算地层水性质的有关参数。,（2）,平衡区数据，如油水界面、油气界面。,（3）油,PVT,表，气,PVT,表。,9.,黑油模型数据流,75,2,岩石数据包括油水系统和油气系统的相对渗透率曲线和毛管压力曲线数据，以及其它的岩石性质数据。主要是：,（1）,输入岩石压缩系数及其参考压力。,（2）,油水系统以及油气系统的相对渗透率曲线和毛管压力曲线数据。,（3）毛管压力端点值与束缚水饱和度关系。,9.,黑油模型数据流,76,3,网格属性数据是用数组表述的属性数据,主要来自地质建模,它们是：,（1）,油层深度数据,（2）,孔隙度数据,（3）渗透率数据,（4）厚度数据,（5）有效厚度,9.,黑油模型数据流,77,初始化程序输入以下几类卡片：,（1）,NIT,卡,（2） 一般描述卡片,（3）,PVT,数据卡片,（4） 岩石类型数据卡片,（5） 数组数据,（6）,END,卡,每一类型卡片有引领关键字，标志者该类数据开始输入。,INIT,卡必须放在卡片叠的第一张，,END,放在卡片叠的最后一张，数组数据放在其它输入数据的后面。,ljfinit.txt,9.,黑油模型数据流,78,二、动态模型数据主要有：,1.,模型控制数据：,（1）,离散方程求解控制参数,（2）,线性代数方程组解法控制参数,（3）时间步长控制参数,2. 输出控制参数:,（1）,控制井资料报告,（2）,控制全油田、分区及分层报告,（3）控制数据场报告：压力、饱和度,9.,黑油模型数据流,79,3.,井数据：,（1）,井位数据及射孔数据,（2）,井的分类及生产方式,ljfmodel.txt,9.,黑油模型数据流,80,谢谢大家！,81,