油藏数值模拟技术课件

油藏数值模拟技术与应用油藏数值模拟技术与应用常晓平河南油田勘探开发研究院油藏数值模拟技术与应用1一、油藏数值模拟技术一、油藏数值模拟技术二、油藏数值模拟的主流软件及发展方向二、油藏数值模拟的主流软件及发展方向三、实例三、实例一、油藏数值模拟技术油藏数值模拟技术与应用2（一）基本概念及作用（一）基本概念及作用 （二）所需的数据及准备方法（二）所需的数据及准备方法 （三）模型初始化方法（三）模型初始化方法 （四）生产史拟合技术（四）生产史拟合技术 （五）动态预测技术（五）动态预测技术 一、油藏数值模拟技术一、油藏数值模拟技术（一）基本概念及作用 一、油藏数值模拟技术3 油油藏藏数数值值模模拟拟技技术术是是一一门门将将油油田田开开发发重重大大决策纳入严格科学轨道的关键技术。决策纳入严格科学轨道的关键技术。从从油油田田投投产产开开始始，无无论论是是单单井井动动态态，还还是是整整个个油油田田动动态态，都都要要进进行行监监测测与与控控制制。油油藏藏数数值值模模拟拟是是油油田田开开发发最最优优决决策策的的有有效效工工具具，成为油田开发中不可缺少的一项重要工具。成为油田开发中不可缺少的一项重要工具。油藏数值模拟技术是一门将油田开发重大决策纳入严格科4一、油藏数值模拟基本概念油藏数值模拟基本概念模拟 仿真：油藏数值模拟就是把油藏在油藏数值模拟就是把油藏在三维的空间里分为许多离散的三维的空间里分为许多离散的单元，并且模拟油藏及流体在单元，并且模拟油藏及流体在空间及一系列离散的时间步里空间及一系列离散的时间步里的发展变化。的发展变化。建立油藏模拟软件，一般包括：（1）质量守恒原理（2）能量守恒原理（3）运动方程（达西定律）（4）状态方程（5）辅助方程（如饱和度方程，毛管力方程等）一、油藏数值模拟基本概念模拟 仿真：建立油藏模拟软件，一般5油藏模拟的作用油藏模拟的作用 1)1)剩余油分布研究。剩余油分布研究。2)2)优化井网、开发层系、井数和井位。优化井网、开发层系、井数和井位。3)3)选择注水方式。选择注水方式。4)4)对油藏和流体性质的敏感性进行研究对油藏和流体性质的敏感性进行研究 5)5)实施方案的可行性评价实施方案的可行性评价油藏模拟的作用 6油藏模拟的基本工作流程油藏模拟的基本工作流程数据文件准备初始化计算生产史拟合开发预测结果整理输出每一步的目的？油藏模拟的基本工作流程数据文件准备初始化计算生产史拟合开发预7（一）基本概念及作用（一）基本概念及作用 （二）所需的数据及准备方法（二）所需的数据及准备方法 （三）模型初始化方法（三）模型初始化方法 （四）生产史拟合技术（四）生产史拟合技术 （五）动态预测技术（五）动态预测技术 一、油藏数值模拟技术一、油藏数值模拟技术（一）基本概念及作用 一、油藏数值模拟技术8所需的数据及其准备方法所需的数据及其准备方法一、所需数据基本类型一、所需数据基本类型：2 2、PVTPVT及岩石性质数据及岩石性质数据 提供流体的地层体积系数、黏度、密度、气油比、岩石及油水提供流体的地层体积系数、黏度、密度、气油比、岩石及油水 的压缩系数、相对渗透率及界面毛管力曲线。的压缩系数、相对渗透率及界面毛管力曲线。3 3、初始化数据、初始化数据提供平衡区，如流体界面、参考深度、压力及毛管力。提供平衡区，如流体界面、参考深度、压力及毛管力。1 1、静态油藏描述数据、静态油藏描述数据 用离散的网格单元构造油藏的一个几何模型。用离散的网格单元构造油藏的一个几何模型。对每一个网格单元提供尺寸、海拔高程、孔隙度及渗透率。对每一个网格单元提供尺寸、海拔高程、孔隙度及渗透率。4 4、井数据、井数据 提供井、完井位置、生产提供井、完井位置、生产/注入速率、井组及其它数据如表皮注入速率、井组及其它数据如表皮 系数、井半径及井控制等系数、井半径及井控制等5 5、其它：、其它：聚合物等。聚合物等。所需的数据及其准备方法一、所需数据基本类型：2、PVT及岩石91 1、网格数据、网格数据地质模型等地质模型等2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区数据、平衡区数据5 5、油藏工程方法、油藏工程方法6 6、数据文件实例、数据文件实例数据准备方法数据准备方法 1、网格数据地质模型等数据准备方法 10(一一)网格的概念网格的概念(二二)网格设计的发展阶段网格设计的发展阶段(三三)网格的类型及描述方法网格的类型及描述方法(四四)网格数据的读取顺序网格数据的读取顺序(五五)各类网格的优缺点各类网格的优缺点网格描述网格描述网格描述网格描述(一)网格的概念网格描述11(一一)网格的概念及设计方法网格的概念及设计方法 1、网格的概念：网格的概念：2、网格描述的目的网格描述的目的:建立能较准确描述油藏建立能较准确描述油藏特征的地质模型。特征的地质模型。(一)网格的概念及设计方法 1、网格的概念：12建立地质模型的网格设计方法：建立地质模型的网格设计方法：1 1、选择模型的几何描述、选择模型的几何描述(一一)网格的概念及设计方法网格的概念及设计方法（1 1）研究区域的大小及形态）研究区域的大小及形态（2 2）需要的研究精度）需要的研究精度（3 3）所得数据的详细程度）所得数据的详细程度（4 4）断层结构的复杂性）断层结构的复杂性（5 5）断层两边地层的接触关系及属性）断层两边地层的接触关系及属性（6 6）存在倾斜及下降断层）存在倾斜及下降断层（7 7）建模的时间限制）建模的时间限制建立地质模型的网格设计方法：1、选择模型的几何描述(一)网13(一一)网格的概念及设计方法网格的概念及设计方法3 3、设计区域网格、设计区域网格1 1、基于顶面构造。、基于顶面构造。2 2、只包含一层。、只包含一层。3 3、井尽可能位于网格中心、井尽可能位于网格中心4 4、断层尽可能位于网格边界、断层尽可能位于网格边界5 5、边界为不流动边界、边界为不流动边界6 6、网格单元的大小及形态的变化随着：、网格单元的大小及形态的变化随着：（1 1）井附近的网格一般较小。）井附近的网格一般较小。（2 2）水体位置的网格一般较大）水体位置的网格一般较大 (一)网格的概念及设计方法3、设计区域网格14(一一)网格的概念及设计方法网格的概念及设计方法在垂向上网格分层基于：在垂向上网格分层基于：1 1、可得的层面数据。、可得的层面数据。2 2、渗透性随深度的变化。、渗透性随深度的变化。3 3、影响垂向上流动的隔挡：、影响垂向上流动的隔挡：如泥岩夹层的影响程度如泥岩夹层的影响程度 及范围等。及范围等。优点：夹层描述精确；弱点：网格多，占用空间大；容易与流量相关的导致不收敛问题优点：网格少，省空间；与流量相关的不收敛问题较少。弱点：垂向上描述较粗，需花较大的功夫调整垂向传导率4 4、设计分层结构、设计分层结构(一)网格的概念及设计方法在垂向上网格分层基于：优点：夹层15(二二)网格设计的发展阶段：网格设计的发展阶段：网格描述网格描述1、块中心、块中心2、角点、角点3、PEBI(二)网格设计的发展阶段：网格描述1、块中心16FF网格技术以块中心粗网格为主网格技术以块中心粗网格为主网格技术以块中心粗网格为主网格技术以块中心粗网格为主FF网格和数学解法在特殊边界网格和数学解法在特殊边界网格和数学解法在特殊边界网格和数学解法在特殊边界(断层、断层、断层、断层、边界、井周围等边界、井周围等边界、井周围等边界、井周围等)不能满足正交不能满足正交不能满足正交不能满足正交FF不能模拟复杂地质结构油气藏不能模拟复杂地质结构油气藏不能模拟复杂地质结构油气藏不能模拟复杂地质结构油气藏FF模拟对象只限于黑油，模拟结果模拟对象只限于黑油，模拟结果模拟对象只限于黑油，模拟结果模拟对象只限于黑油，模拟结果误差较大误差较大误差较大误差较大1、块中心网格、块中心网格：1960-1980 1960-1980 1960-1980 1960-19801 15 510101 1151520205 5101015152020X X X XY Y Y Y(二二)网格设计的发展阶段网格设计的发展阶段网格技术以块中心粗网格为主1、块中心网格：1960-19817FF8080年代中期发展了角点网格年代中期发展了角点网格技术技术FF角点网格着重于遵循油藏的角点网格着重于遵循油藏的几何形态和地质结构特征几何形态和地质结构特征,但但却不满足数学上的正交性却不满足数学上的正交性.导导致油藏构造越复杂致油藏构造越复杂,其模拟误其模拟误差越大差越大2、角点网格、角点网格：1981-19991981-19991981-19991981-1999FaultFault(二二)网格设计的发展阶段网格设计的发展阶段80年代中期发展了角点网格技术2、角点网格：1981-199183 3、PEBIPEBI网格：网格：1999 1999 1999 1999(二二)网格设计的发展阶段网格设计的发展阶段可模拟任何几何形状的油藏可模拟任何几何形状的油藏可模拟任何几何形状的油藏可模拟任何几何形状的油藏用用用用k-PEBIk-PEBIk-PEBIk-PEBI网格处理油藏的各向异性网格处理油藏的各向异性网格处理油藏的各向异性网格处理油藏的各向异性加密网格与基础网格能够自动耦合加密网格与基础网格能够自动耦合加密网格与基础网格能够自动耦合加密网格与基础网格能够自动耦合收敛速度快稳定性高误差极小收敛速度快稳定性高误差极小收敛速度快稳定性高误差极小收敛速度快稳定性高误差极小正交极小加速算法正交极小加速算法正交极小加速算法正交极小加速算法(Orthomin)(Orthomin)(Orthomin)(Orthomin)3、PEBI网格：1999(二)网格设计的发展阶段可19更具灵活性3 3、PEBIPEBI网格：网格：1999 1999 1999 1999 网格设计方法的发展阶段网格设计方法的发展阶段网格设计方法的发展阶段网格设计方法的发展阶段更具3、PEBI网格：1999 网格设计方法的发展阶段20(三三)网格的类型及描述方法网格的类型及描述方法网格描述网格描述到目前主要有三种网格类型：块中心、角点及PEBI网格。(三)网格的类型及描述方法网格描述到目前主要有三种网格类型21网格几何描述块中心网格：给出DX、DY、DZ及深度角点网格：需要指定组成每一个网格的四条坐标线坐标（COORD）及八个角点的深度（ZCORN）网格几何描述块中心网格：给出DX、DY、DZ及深度角点网格：22块中心网格描述倾斜构造在块中心网格系统，虽然相邻网格的深度错开，但ECLIPSE还是认为其间是相连的，有流动发生。块中心网格描述倾斜构造在块中心网格系统，虽然相邻网格的深度错23角点网格描述倾斜构造角点网格描述倾斜构造24块中心、角点网格的例子块中心、角点网格的例子25水平井水平井水平井水平井PEBIPEBI网格网格网格网格直井直井直井直井PEBIPEBI网格网格网格网格PEBI 网格生产区生产区生产区生产区PEBIPEBI网格网格网格网格水平井PEBI网格直井PEBI网格PEBI 网格生产区PEB26PEBI 网格非结构性非结构性非结构性非结构性PEBIPEBI网格网格网格网格PEBI 网格非结构性27PEBI 网格非结构性非结构性非结构性非结构性PEBIPEBI网格网格网格网格PEBI 网格非结构性28n n一口注水井一口注水井一口注水井一口注水井n n两口采油井两口采油井两口采油井两口采油井n n均质油藏均质油藏均质油藏均质油藏n n5818 15818 15818 15818 1n n断层不封闭断层不封闭断层不封闭断层不封闭n n运行运行运行运行3 3 3 3年年年年,对比水对比水对比水对比水饱和度分布饱和度分布饱和度分布饱和度分布PEBIPEBI网格技术网格技术-PEBIPEBI网格与角点网格对比网格与角点网格对比一口注水井PEBI网格技术-PEBI网格与角点网格对比29直角网格直角网格PEBIPEBI网格技术网格技术-PEBIPEBI网格与角点网格对比网格与角点网格对比直角网格PEBI网格技术-PEBI网格与角点网格对比30模拟模拟模拟模拟3 3年后含水饱和度分布年后含水饱和度分布年后含水饱和度分布年后含水饱和度分布模拟3年后含水饱和度分布31角点网格角点网格角点网格角点网格Prod1Prod2Inj1角点网格Prod1Prod2Inj132角点网格角点网格3 3年后含水饱和度分布年后含水饱和度分布年后含水饱和度分布年后含水饱和度分布角点网格3年后含水饱和度分布33加入断层后所形成的加入断层后所形成的加入断层后所形成的加入断层后所形成的PEBIPEBI 网格网格网格网格加入断层后所形成的PEBI 网格34用用用用SURE-PEBISURE-PEBI 模型模拟模型模拟模型模拟模型模拟3 3年后含水饱和度分布年后含水饱和度分布年后含水饱和度分布年后含水饱和度分布用SURE-PEBI 模型模拟3年后含水饱和度分布35笛卡尔网格读取顺序笛卡尔网格读取顺序36(一一)网格的概念网格的概念(二二)网格设计的发展阶段网格设计的发展阶段(三三)网格的类型及描述方法网格的类型及描述方法(四四)网格数据的读取顺序网格数据的读取顺序(五五)各类网格的优缺点各类网格的优缺点网格描述网格描述网格描述网格描述(一)网格的概念网格描述37(六六)各类网格的优缺点各类网格的优缺点网格描述网格描述网格描述网格描述网格类型优点缺点块中心简单、正交性好描述断层等误差较大角点描述地层、特别是断层的能力较强较复杂，正交性不好，需借助专用软件生成PEBI描述各类地层的能力强，正交性好。较复杂，需借助专用软件生成(六)各类网格的优缺点网格描述网格类型优点缺点块中心简单、381 1、网格描述、网格描述2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区、平衡区5 5、油藏工程方法、油藏工程方法6 6、数据文件实例、数据文件实例数据准备方法数据准备方法 1、网格描述数据准备方法 39PVT数据数据 来源于来源于PVT分析，数据以表的方式分析，数据以表的方式描述有关参数与压力、温度的关系。描述有关参数与压力、温度的关系。数据项：数据项：1、各相的地层体积系数（、各相的地层体积系数（Bo、Bg、Bw）2、各相的黏度、各相的黏度 3、气油比、气油比Rs 和（或）油气比和（或）油气比Rv）PVT数据 来源于PVT分析，数据以表的方式40相态图分析区域A为黑油，该区域或许经过泡点线，并且包含溶解气，但是距离临界点较远DD线：表示接近临界点的流体。这时很难确切的知道流体是液体还是气体，因为流体沿着DD经过临界点时，变成了两相混合物并且不易确定过度区。HH线：非等温过程，如在分离器中时。相态图分析区域A为黑油，该区域或许经过泡点线，并且包含溶解气41黑油模拟器黑油模拟器l严格地说，严格地说，黑油黑油在地面条件下不包含在地面条件下不包含溶解气。溶解气。l黑油模拟器最好用在远离临界点的单黑油模拟器最好用在远离临界点的单相区域。相区域。l如果流体的组成变化比较大，就不应如果流体的组成变化比较大，就不应该使用黑油模拟器。该使用黑油模拟器。黑油模拟器严格地说，黑油在地面条件下不包含溶解气。42基本概念基本概念黑油模型：黑油模型：黑油模型是描述含有非挥发组分的黑油和挥发性组分的黑油模型是描述含有非挥发组分的黑油和挥发性组分的原油溶解气在油藏中运动规律的数值模型。黑油模型也称低挥原油溶解气在油藏中运动规律的数值模型。黑油模型也称低挥发油双组分模型。发油双组分模型。组分模型：组分模型：组分模型是用来研究高挥发性烃类的系统在油气藏中运动组分模型是用来研究高挥发性烃类的系统在油气藏中运动规律的数值模型，也成多组分模型。规律的数值模型，也成多组分模型。基本概念黑油模型：组分模型：43黑油模拟器的适用条件黑油模拟器的适用条件1、因为凝析及压力降低而释放出来的气体、因为凝析及压力降低而释放出来的气体的总量只占总含烃量的很少一部分。的总量只占总含烃量的很少一部分。2、由于凝析及压力降低而释放出来气体后，、由于凝析及压力降低而释放出来气体后，流体的烃组成变化不明显。流体的烃组成变化不明显。3、路径远离临界点。、路径远离临界点。4、整个过程是等温的。、整个过程是等温的。如果上述条件不能完全满足，则需要使如果上述条件不能完全满足，则需要使用全组分模拟器。用全组分模拟器。为了保证使用黑油模拟器可以达到足够的精度，为了保证使用黑油模拟器可以达到足够的精度，油藏流体必须具备下列条件：油藏流体必须具备下列条件：黑油模拟器的适用条件1、因为凝析及压力降低而释放出来的气体441 1、网格描述、网格描述2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区、平衡区5 5、油藏工程方法、油藏工程方法6 6、数据文件实例、数据文件实例数据准备方法数据准备方法 1、网格描述数据准备方法 45l来源于来源于SCAL（special core analysis）分析，）分析，数据主要是以表的方式描述属性与饱和度的数据主要是以表的方式描述属性与饱和度的关系。关系。l对于油藏的不同部分，可分用不同的相渗曲对于油藏的不同部分，可分用不同的相渗曲线。线。l在每个区域，需要设置最大、最小及临界饱在每个区域，需要设置最大、最小及临界饱和度值。和度值。l用于定义过渡带的饱和度。用于定义过渡带的饱和度。岩心分析岩心分析 岩心分析46岩石数据岩石数据 岩石数据是特定的岩心分析试验的结果，岩石数据是特定的岩心分析试验的结果，该数据用于：该数据用于：l设置每一流体相的最大、最小饱和度，该设置每一流体相的最大、最小饱和度，该值用于定义平衡区的相饱和度。值用于定义平衡区的相饱和度。l定义过渡带的范围及属性。定义过渡带的范围及属性。l描述各相在网格块间流动时的动态表现。描述各相在网格块间流动时的动态表现。岩石数据471 1、网格描述、网格描述2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区、平衡区5 5、油藏工程方法、油藏工程方法6 6、数据文件实例、数据文件实例数据准备方法数据准备方法 1、网格描述数据准备方法 48 平衡区平衡区-概念及作用概念及作用 在原始状况下，认为在同一个在原始状况下，认为在同一个平衡区中的流体的静力学是平衡平衡区中的流体的静力学是平衡的。的。概念：概念：作用：作用：用于定义油藏初始条件下的用于定义油藏初始条件下的每一个网格单元的压力、饱和度每一个网格单元的压力、饱和度值。值。平衡区-概念及作用 49平衡区-初始平衡计算（一）初始压力计算：平衡区-初始平衡计算（一）初始压力计算：50原始水分布l也可以在PROPS 部分中使用SWATINIT 定义每个单元的原始含水饱和度值。SW变化相对较为剧烈在。原始水分布也可以在PROPS 部分中使用SWATINIT 定511 1、网格描述、网格描述2 2、PVTPVT分析分析3 3、岩心分析、岩心分析4 4、平衡区、平衡区5 5、油藏工程方法、油藏工程方法数据准备方法数据准备方法 1、网格描述数据准备方法 52网格参数场计算公式网格参数场计算公式53单井动态计算公式l井流量方程：l单井的完井传导率系数：l井的压力等值系数：如果有效的井眼半径超过了re的一半，将会出现警告；如果有效的井眼半径超过了re，Twj变成了负值从而导致了错误的结果l井的生产指数：单井动态计算公式井流量方程：单井的完井传导率系数：井的压力等54（一）基本概念及作用（一）基本概念及作用 （二）所需的数据及准备方法（二）所需的数据及准备方法 （三）模型初始化方法（三）模型初始化方法 （四）生产史拟合技术（四）生产史拟合技术 （五）动态预测技术（五）动态预测技术 一、油藏数值模拟技术一、油藏数值模拟技术（一）基本概念及作用 一、油藏数值模拟技术55初始化阶段应完成的主要工作初始化阶段应完成的主要工作l数据查错：1、文本方式 -*.PRT 2、OFFICE 环境中。l储量拟合：1、油水界面深度检查。2、油藏每一层边界检查。3、地质模型检查：Q=A*H*POR*SO初始化阶段应完成的主要工作56数模软件在初始化时的工作1、软件依次读入及处理各部分数据。在开始读入下部软件依次读入及处理各部分数据。在开始读入下部分内容前，必须对本部分各种数据进行一致性检查。分内容前，必须对本部分各种数据进行一致性检查。只有最后一部分例外，因为该部分数据是时间相关只有最后一部分例外，因为该部分数据是时间相关的，因此，总体上初始化时并不读入与处理。的，因此，总体上初始化时并不读入与处理。2、软件开始运行的第一个任务是为输入的数据分配内存。、软件开始运行的第一个任务是为输入的数据分配内存。3、处理模拟网格的几何情况及属性，得到更利于流动计算、处理模拟网格的几何情况及属性，得到更利于流动计算的三维空隙体积、传导率、单元中部深度及可产生流动的三维空隙体积、传导率、单元中部深度及可产生流动的网格间的连接。的网格间的连接。4、在初始化（在初始化（initialization）时，计算油藏每一层的静压力）时，计算油藏每一层的静压力剃度，并给每一个网格赋每一相的饱和度值剃度，并给每一个网格赋每一相的饱和度值数模软件在初始化时的工作1、软件依次读入及处理各部分数据57给定初始化数据的方法l给定平衡区数据：平衡区数据被用来初始化平衡区数据被用来初始化饱和度分布、基准深度的压饱和度分布、基准深度的压力及压力梯度。力及压力梯度。l用枚举的方式给出每一个单元的初始饱和度及压力。给定初始化数据的方法给定平衡区数据：58（一）基本概念及作用（一）基本概念及作用 （二）所需的数据及准备方法（二）所需的数据及准备方法 （三）模型初始化方法（三）模型初始化方法 （四）生产史拟合技术（四）生产史拟合技术 （五）动态预测技术（五）动态预测技术 一、油藏数值模拟技术一、油藏数值模拟技术（一）基本概念及作用 一、油藏数值模拟技术59 生产史拟合 修改现有的模拟数据（如渗透率、相渗曲线、修改现有的模拟数据（如渗透率、相渗曲线、井的生产指数等）直到能与观察数据（如压力、井的生产指数等）直到能与观察数据（如压力、产量、含水等）进行合理的对比，这种过程叫做产量、含水等）进行合理的对比，这种过程叫做生产史拟合。生产史拟合。生产史拟合 修改现有的模拟数据（如渗透率、相渗曲60生产史拟合生产史拟合 生产史拟合过程也就是一个反复修改参数、反复试算的生产史拟合过程也就是一个反复修改参数、反复试算的过程，需要耗费大量的机时和人力。过程，需要耗费大量的机时和人力。拟合结果的好坏，不但拟合结果的好坏，不但取决于人们对油藏的认识，还取决于油藏工程师的拟合经验取决于人们对油藏的认识，还取决于油藏工程师的拟合经验和处理技巧。和处理技巧。由于油藏地质情况的复杂性及反求参数的多解性，不同由于油藏地质情况的复杂性及反求参数的多解性，不同的人拟合同一油田，其试算次数和拟合结果可能均不同。的人拟合同一油田，其试算次数和拟合结果可能均不同。为了减少拟合工作的盲目性，有必要了解各项油藏参数为了减少拟合工作的盲目性，有必要了解各项油藏参数与拟合对象的相互关系，同时还应一些拟合经验和处理技巧。与拟合对象的相互关系，同时还应一些拟合经验和处理技巧。生产史拟合 生产史拟合过程也就是一个反复61生产史拟合 拟合的对象有：拟合的对象有：压力、流量、油压力、流量、油气比、油水比等气比、油水比等。生产史拟合 拟合的对象有：压力、流量、油气比、油水比等。62生产史拟合的步骤：l在油藏的历史拟合中，一般把产量作为已知条件来拟合其它动态参数。拟合步骤为：开始油藏原始平衡状态检查（零流量模拟）拟合井底压力拟合油田平均压力拟合单井压力拟合油田综合含水率拟合单井含水率油藏压力已经拟合好了？结束生产指数控制着单井流量开始下降的时间，对预测单井的动态意义较大。生产史拟合的步骤：在油藏的历史拟合中，一般把产量作为已知条件63生产史拟合 可以单个的或集体的改变几个参数。通常按下述范可以单个的或集体的改变几个参数。通常按下述范围进行修改：围进行修改：1、岩石数据的修改：、岩石数据的修改：a.渗透率渗透率 b.孔隙度孔隙度 c.厚度厚度 d.饱和度饱和度 2、流体数据的修改：、流体数据的修改：a.压缩性压缩性 b.PVT数据数据 3、相对渗透率数据：、相对渗透率数据：a.相对渗透率曲线的移位相对渗透率曲线的移位 b.临界饱和度数据的移位临界饱和度数据的移位 4、单井完井数据：、单井完井数据：a.表皮效应表皮效应 b.井底流动压力。井底流动压力。3 3、可修改的参数、可修改的参数、可修改的参数、可修改的参数生产史拟合 可以单个的或集体的改变几64生产史拟合 在历史拟合过程中，由于模型参数数量多，在历史拟合过程中，由于模型参数数量多，可调的自由度很大。为了避免或减少修改参数的可调的自由度很大。为了避免或减少修改参数的随意性，在历史拟合开始时，必须确定模型参数随意性，在历史拟合开始时，必须确定模型参数的可调范围，使模型参数的修改在合理的、可接的可调范围，使模型参数的修改在合理的、可接受的范围内，这是历史拟合的原则。受的范围内，这是历史拟合的原则。4 4 4 4、确定模型参数的可调范围、确定模型参数的可调范围、确定模型参数的可调范围、确定模型参数的可调范围生产史拟合 在历史拟合过程中，由于模型参数数量多，65 确定参数的可调范围是一项重要而细致的工确定参数的可调范围是一项重要而细致的工作，需要油藏工程师和地质师共同努力，需要收作，需要油藏工程师和地质师共同努力，需要收集和分析一切可以利用的资料。首先分清哪些参集和分析一切可以利用的资料。首先分清哪些参数是确定的，即准确可靠的，哪些参数是不确定数是确定的，即准确可靠的，哪些参数是不确定的，即不准确、不可靠的，然后根据情况，确定的，即不准确、不可靠的，然后根据情况，确定可调范围。通常只修改不确定的参数，而确定的可调范围。通常只修改不确定的参数，而确定的参数一般不允许修改。参数一般不允许修改。确定参数的可调范围是一项重要而细致的工作，需要油66确定模型参数的可调范围确定模型参数的可调范围确定模型参数的可调范围确定模型参数的可调范围 一般情况下，认为：一般情况下，认为：（1 1）孔隙度：）孔隙度：如果油层大量岩心分析资料表明，油层部分孔隙度如果油层大量岩心分析资料表明，油层部分孔隙度在在1919到到2121之间，平均为之间，平均为2020，变化范围不大。则，变化范围不大。则把孔隙度视为确定参数，不做修改，或允许改动范围在把孔隙度视为确定参数，不做修改，或允许改动范围在33。（2 2）渗透率：渗透率：渗透度在任何油田都是不定参数。这不仅是由于测渗透度在任何油田都是不定参数。这不仅是由于测井解释的渗透率值和岩心分析值误差较大，而且根据渗井解释的渗透率值和岩心分析值误差较大，而且根据渗透率的特点，井间的渗透率分布也是不确定的。因此对透率的特点，井间的渗透率分布也是不确定的。因此对渗透率的修改，允许范围较大，可放大或缩小渗透率的修改，允许范围较大，可放大或缩小2 23 3倍或倍或更多。更多。确定模型参数的可调范围 一般情况下，认为：（2）渗透率67确定模型参数的可调范围确定模型参数的可调范围确定模型参数的可调范围确定模型参数的可调范围（3）岩石与液体的压缩系数：）岩石与液体的压缩系数：液液体体的的压压缩缩系系数数是是实实验验测测定定的的，变变化化范范围围很很小小，认认为为是是确确定定的的。而而岩岩石石的的压压缩缩系系数数虽虽然然也也是是实实验验室室测测定定的的，但但受受岩岩石石内内饱饱和和液液体体和和应应力力状状态态的的影影响响，有有一一定定变变化化范范围围，考虑这部分影响，允许岩石的压缩系数可以扩大一倍。考虑这部分影响，允许岩石的压缩系数可以扩大一倍。（4）初始流体饱和度和初始压力：）初始流体饱和度和初始压力：和和通通常常做做法法一一样样，认认为为是是确确定定参参数数。必必要要时时允允许许小小范范围内修改。围内修改。确定模型参数的可调范围（3）岩石与液体的压缩系数：（468确定模型参数的可调范围确定模型参数的可调范围确定模型参数的可调范围确定模型参数的可调范围（5）相对渗透率曲线：）相对渗透率曲线：由由于于油油藏藏模模拟拟模模型型的的网网格格粗粗，网网格格内内部部存存在在严严重重非非均均质质，其其影影响响不不可可忽忽视视，这这与与均均质质岩岩心心的的情情况况不不同同。因因此此相相对对渗渗透透率率曲曲线线应应看看作作是是不不定定参参数数。在在拟拟曲曲线线的的研研究究中中，给出了较好的初始值，但仍允许做适当修改。给出了较好的初始值，但仍允许做适当修改。（6 6）油气的）油气的PVTPVT性质：性质：视为确定参数。视为确定参数。（7 7）油水界面：）油水界面：在资料不多的情况下，允许在一定范围内修改。在资料不多的情况下，允许在一定范围内修改。确定模型参数的可调范围（5）相对渗透率曲线：（6）油气的P69生产史拟合生产史拟合-油藏原始平衡状态检查油藏原始平衡状态检查 当油藏未投产时，油藏当油藏未投产时，油藏压力系统处于静止平衡状态，压力系统处于静止平衡状态，相邻点不应有压力势差；否相邻点不应有压力势差；否则会出现流动，打破静态平则会出现流动，打破静态平衡。衡。生产史拟合-油藏原始平衡状态检查 70生产史拟合-全区及单井压力拟合之一 当油藏局部区域出现高压异常与低压异常时：（1）沿高压异常区向低压异常区方向增大渗透率值，可使两类异常区同时得到改善。但修改K值应在渗透率图上成片修改，以保持岩性的连续性。（2）若修改K值效果不佳，则可减少高压异常区储量，增加低压异常区储量。生产史拟合-全区及单井压力拟合之一 当油藏局部区域出现71生产史拟合-全区及单井压力拟合之二 当油藏采出大量液体后压力无当油藏采出大量液体后压力无明显下降，可考虑是否因岩石弹性明显下降，可考虑是否因岩石弹性能量过大所致。可减少能量过大所致。可减少岩石压缩系岩石压缩系数数；当油藏压力下降过快时，也可；当油藏压力下降过快时，也可增加岩石压缩系数加以改善。但此增加岩石压缩系数加以改善。但此时还应考率到是否是油藏储量偏小时还应考率到是否是油藏储量偏小所致。若是，则应增加所致。若是，则应增加油藏储量油藏储量。生产史拟合-全区及单井压力拟合之二 当油藏采出72生产史拟合-含水率拟合一单元含水率主要受相对渗透率曲线的影响。含水率主要受相对渗透率曲线的影响。油井的见水时间和水驱前缘的拟合比较难以实现。它油井的见水时间和水驱前缘的拟合比较难以实现。它还受到模型的时间步长和网格密度的限制还受到模型的时间步长和网格密度的限制。当见水区域或油井含水上升过快时，可以通过调整油当见水区域或油井含水上升过快时，可以通过调整油相、水相相对渗透率曲线来控制。增大水相临界饱和度值相、水相相对渗透率曲线来控制。增大水相临界饱和度值可以减缓见水，而增大可以减缓见水，而增大Kro/KrwKro/Krw的值可降低含水上升速度。的值可降低含水上升速度。生产史拟合-含水率拟合一单元含水率主要受相对渗透率曲73生产史拟合-含水率拟合二 在拟合上油藏综合含水率之后，还有必要对某些在拟合上油藏综合含水率之后，还有必要对某些单井含水率进行拟合。单井含水率进行拟合。单井含水率的拟合主要通过调整单井含水率的拟合主要通过调整平面上平面上的渗透的渗透率场及率场及纵向上纵向上的生产指数两个方面来做。的生产指数两个方面来做。生产史拟合-含水率拟合二 在拟合上油藏综合含水率之74生产史拟合生产史拟合拟合井底流压拟合井底流压拟合井底流压拟合井底流压 主要是通过修改生产井的采油指数（包括表皮系数）主要是通过修改生产井的采油指数（包括表皮系数）来实现。流动压力偏小，可增大采油指数；反之，可减小来实现。流动压力偏小，可增大采油指数；反之，可减小采油指数。采油指数。生产史拟合拟合井底流压 主要是通过修改生产井的75生产史拟合 岩石数据通常是二维或三维网格矩阵输入值。修改网格数据从不按单个单元来进行，而是在需作必要变化的面积上进行的。1 1、岩石数据的修改：、岩石数据的修改：、岩石数据的修改：、岩石数据的修改：生产史拟合 岩石数据通常是二维或三维网格矩阵输入值。76生产史拟合 多孔岩石内膨胀能量的数量的大小直接与三种流体相的多孔岩石内膨胀能量的数量的大小直接与三种流体相的饱和度有关，油层能量的饱和度有关，油层能量的影响影响表现在几方面：表现在几方面：1.油层压力水平油层压力水平 2.产量产量 3.油气比油气比 生产史拟合 多孔岩石内膨胀能量的数量的大小直接与77生产史拟合流体饱和度的修改流体饱和度的修改流体饱和度的修改流体饱和度的修改 当计算中无法保持一定地区给定的产量、进而油层压力当计算中无法保持一定地区给定的产量、进而油层压力降低较严重时，通常这是由于油层中储油量不充足和流入给降低较严重时，通常这是由于油层中储油量不充足和流入给定位置的流量不够所引起的。定位置的流量不够所引起的。这时除了改变这时除了改变孔隙体积孔隙体积可以校正这种情况以外，也可通可以校正这种情况以外，也可通过过提高油饱和度提高油饱和度来进行校正。来进行校正。生产史拟合流体饱和度的修改 当计算中无法保持78生产史拟合流体数据的修改流体数据的修改流体数据的修改流体数据的修改 在模拟研究中，流体数据一般是一知的。这些数据包括：在模拟研究中，流体数据一般是一知的。这些数据包括：1.地层体积系数地层体积系数 2.粘度粘度 3.压缩性压缩性 4.溶解气数据溶解气数据 双河双河IV下数模研究下数模研究生产史拟合流体数据的修改 在模拟研究中，流79（一）基本概念及作用（一）基本概念及作用 （二）所需的数据及准备方法（二）所需的数据及准备方法 （三）模型初始化方法（三）模型初始化方法 （四）生产史拟合技术（四）生产史拟合技术 （五）动态预测技术（五）动态预测技术 一、油藏数值模拟技术一、油藏数值模拟技术（一）基本概念及作用 一、油藏数值模拟技术80动态预测动态预测 动动态态预预测测可可用用来来研研究究预预测测包包括括枯枯竭竭开开采采、注注水开发、钻加密井等各种方案的效果。水开发、钻加密井等各种方案的效果。动态预测 动态预测可用来研究预测包括枯竭开采、注水81预测方法：预测方法：1.设置单井的生产目标2.设置井组的生产目标预测方法：设置单井的生产目标82设定单井的动态控制模式及目标：设定单井的动态控制模式及目标：井底流压井口压力产油率产水率产气率产液率井组产率设定单井的动态控制模式及目标：井底流压83 压力控制一般用于预测阶段，一般不适用于生产史拟合阶段。在压力控制下，模拟器计算相应的流量；如果计算的流量超过了指定的速率，则会切换回速率控制模式下。压力控制一般用于预测阶段，一般不适用于生产史拟合阶段84示例：一注一采两口井的剖面模型：示例：一注一采两口井的剖面模型：85生产井中的控制在预测阶段自动切换的示例生产井中的控制在预测阶段自动切换的示例BHP产油量产水量产油/水量压力控制方式控制方式生产井中的控制在预测阶段自动切换的示例BHP产油量产水量产油86定生产压差求产：定生产压差求产：87 经济限制是通过修改生产策略来满足所定的最低经济目标（如产油率、油气比等）的一种方法。当抵达经济限制条件时，根据事先的定义，一系列的措施就会自动实施。以经济极限制为约束求产：经济限制是通过修改生产策略来满足所定的最低经济目标（88采用经济限制可设定：采用经济限制可设定：89经济限制可应用于：油田/井组 产量GECON井的产量WECON层的产量（ratios only）CECON经济限制可应用于：油田/井组 产量GECON90下面的措施在计算满足经济限制条件时，可自动实施：关井。关井。关层。关层。打水泥塞。打水泥塞。下面的措施在计算满足经济限制条件时，可自动实施：关井。91采用经济限制用法说明1用于在需要时自动进行作业及关井。当抵达含水率、GOR、OGR以及流量限制条件时，自动进行作业一旦井、层关闭后，缺省情况下，要一直保持关闭状态可以用WTEST通过周期性测试将井重新打开。采用经济限制用法说明1用于在需要时自动进行作业及关井。92采用经济限制用法说明2如果指定了第二含水率限制条件，井可以被再次打开。如果指定了第二含水率限制条件，井可以被再次打开。如果井达到经济限制条件，模拟可以被合法地终止。如果井达到经济限制条件，模拟可以被合法地终止。在井名中可以使用通配符。在井名中可以使用通配符。采用经济限制用法说明2如果指定了第二含水率限制条件，井可以被93井组的经济限制设置方法作业方法与作业方法与WECON一样一样井组的经济限制设置方法作业方法与WECON一样94层段的经济条件限制方法：层段的经济条件限制方法：95自动测试开井：自动测试开井：通常，井通常，井/层关闭后就不再被打开。自动测层关闭后就不再被打开。自动测试用来重新打开原来因无力或经济原因关闭的试用来重新打开原来因无力或经济原因关闭的井井自动测试开井：通常，井/层关闭后就不再被打开。自动测96自动测试开井：自动测试开井：971、井名2、测试间隔时间3、测试当时关闭井的原因及目前的有效性检查 P：物理原因（如BHP、THP）E：经济原因 G：井组或油田目标限制 D：超过最大的THP C：层的关闭原因4、测试的最大时间步5、重新开井时间，天1、井名981.设置单井的生产目标2.设置井组的生产目标预测方法：预测方法：设置单井的生产目标预测方法：99井组控制：井组控制：100用WELSPECS将井放在井组中生产井的产量用GCONPROD关键字注入井的流量用GCONINJE关键字生产井与注入经应该分别放在不同的井组中FIELD 永远存在在井组中的井应分别在WCONPROD及WCONINJE中指定控制模式为GRUP。井组控制：用WELSPECS将井放在井组中井组控制：101井组井组/油田油田 生产控制，井的生产是以相对于指导速率生产控制，井的生产是以相对于指导速率（guide rateguide rate）的比率，该比率服从于他自己的速率及压）的比率，该比率服从于他自己的速率及压力限制，指导速率缺省为该井的生产能力。力限制，指导速率缺省为该井的生产能力。井组控制：井组控制：102亏空补偿（亏空补偿（voidage replacemnetvoidage replacemnet）:如确定注采比如确定注采比（voidage replacement ratiovoidage replacement ratio）优先级控制（优先级控制（prioritypriority）:定义一个打开井次序的优先级序定义一个打开井次序的优先级序列，在完成井组或油田目标的前提下，根据需要、按顺序列，在完成井组或油田目标的前提下，根据需要、按顺序打开不同的井。打开不同的井。井组控制：亏空补偿（voidage replacemnet）:如确定注103本井组的产量在父井组中占的百分比井组中的生产控制：该井组产量是由子井组及井完成的，其中各井组及井的产量该井组产量是由子井组及井完成的，其中各井组及井的产量服从各自定义的生产量百分数（服从各自定义的生产量百分数（GUIDE RATESGUIDE RATES）。）。本井组的产量在父井组中占井组中的生产控制：该井组产量是由子井104在井组目标控制下不在井组目标控制下，以单井控制为准指定井及井组是否处于井组目标控制下：在井组目标控制下不在井组目标控制下，以单井控制为准指定井及井105定气油比生产：通过采取定气油比的方式，可以控制日产油量。定气油比生产：通过采取定气油比的方式，可以控106定气油比生产：备注：如果base GOR 为零或负数，则自动关闭该井的气油比控制。定气油比生产：备注：如果base GOR 为零或负数，则自动107WCONPROD O1 OPEN LRAT 3*60 5*/O2 OPEN LRAT 3*75 5*/O3 OPEN LRAT 3*75 5*/O4 OPEN LRAT 3*90 5*/O5 OPEN LRAT 3*60 5*/O6 OPEN LRAT 3*75 5*/GCONINJE G1 WATER VREP 3*1.2/定注采比为1.2/WCONINJE -I1 WATER OPEN GRUP 150/-I2 WATER OPEN GRUP 150 /I1 WATER OPEN GRUP /I2 WATER OPEN GRUP /定注采比生产：井组控制示例：WCONPROD 108井组控制示例：井组控制示例：109井组注入控制：亏空补偿系数：井组注入量=补偿系数*采液量-注入的其它相。亏空补偿系数被通常用来设定为注采比值。井组注入控制：亏空补偿系数：110井打开顺序优先级设置：PRIORITY设置确定每口井优先级的计算方法，格式为：如果A=D=E=F=H=0，B=C=G=1 则：P=(Q0+Qw)/Qw,就是说含水率的倒数越高，优先级越高，也即是说，含水率越小的井优先级越高。示例井打开顺序优先级设置：PRIORITY设置确定每口井优先级的111井打开顺序优先级设置：当不用PRIORITY计算优先级时，可用WELPRI逐井定义各自的优先级。WELPRI-1 2井名 优先级号井打开顺序优先级设置：当不用PRIORITY计算优先级时，可112设置在优先级控制下（GCONPRI）的井组目标：设置在优先级控制下（GCONPRI）的井组目标：113 油田产油率为10000，以PLAT开头的井组的产液率限制为6000，井按优先级顺序被打开直至达到组产量限制值。井组产气的限制值为20000，如果达到限制值，则对最差的井进行作业措施。示例 油田产油率为10000，以PLAT开头的井组的产液率限114 在动态预测阶段，根据所处的开发阶段、计划采取的开发方式以及研究目的的不同所采取的动态预测方式也有所不同。主要有：定流量求产及定压求产两种预测方式。在动态预测阶段，根据所处的开发阶段、计划采取的开发方115动态预测工作流程：1.在sch.inc模块中使用WELSPECS、WCOMPDAT以及WCONPROD/WCONINJE 分别定义井位、生产层位及流量（如果需要定义井组的流量，则需使用GCONPROD及GCONINJE来定义）。2.检查sum.inc文件中是否定义了足够的输出项3.数模计算4.数据整理动态预测工作流程：在sch.inc模块中使用WELSPECS116 在在预预测测阶阶段段，首首先先要要确确定定预预测测控控制制方方式式，通通常常情情况况下下，都都是是采采用用定定流流量量求求产产的的方方式式-即：设定所有井的注采量预测开发效果。即：设定所有井的注采量预测开发效果。在预测阶段，首先要确定预测控制方式，通常情况下，117单井定义示例：单井定义示例：118井组定义示例：注采比1.0井组定义示例：119一、油藏数值模拟技术一、油藏数值模拟技术二、油藏数值模拟的主流软件及发展方向二、油藏数值模拟的主流软件及发展方向三、实例三、实例一、油藏数值模拟技术油藏数值模拟技术与应用120油藏数值模拟主流软件系统：油藏数值模拟主流软件系统：1、ECLIPSE2、VIP3、WorkBench4、CMG5、SURE油藏数值模拟主流软件系统：1、ECLIPSE121 50 50 50 50年代，开展数值模拟研究。年代，开展数值模拟研究。年代，开展数值模拟研究。年代，开展数值模拟研究。60 60 60 60年代，气、水两相和三相黑油油藏问题的求解。年代，气、水两相和三相黑油油藏问题的求解。年代，气、水两相和三相黑油油藏问题的求解。年代，气、水两相和三相黑油油藏问题的求解。70 70 70 70年代，模拟常规递减和保持压力以外的新方法。年代，模拟常规递减和保持压力以外的新方法。年代，模拟常规递减和保持压力以外的新方法。年代，模拟常规递减和保持压力以外的新方法。到到到到80808080年年年年代代代代，油油油油藏藏藏藏模模模模拟拟拟拟进进进进入入入入商商商商品品品品化化化化阶阶阶阶段段段段：大大大大到到到到一一一一个个个个油油油油公公公公司司司司，小小小小到到到到一一一一个个个个企企企企业业业业普普普普遍遍遍遍使使使使用用用用。在在在在模模模模型型型型上上上上，形形形形成成成成一一一一系系系系列列列列可可可可以以以以处处处处理理理理各各各各种种种种各各各各样样样样复复复复杂杂杂杂问问问问题题题题的的的的模模模模型型型型，如如如如黑黑黑黑油油油油模模模模型型型型、组组组组分分分分模模模模型型型型、热热热热采采采采模模模模型型型型以以以以及及及及各种三次采油用的化学驱模型等。各种三次采油用的化学驱模型等。各种三次采油用的化学驱模型等。各种三次采油用的化学驱模型等。油藏模拟技术的发展过程油藏模拟技术的发展过程 ：50年代，开展数值模拟研究。油藏模拟技术的发展过程 122 (1)(1)(1)(1)新的网格化技术。新的网格化技术。新的网格化技术。新的网格化技术。(2)(2)(2)(2)细网格和精细模拟。细网格和精细模拟。细网格和精细模拟。细网格和精细模拟。(3)(3)(3)(3)开窗技术。开窗技术。开窗技术。开窗技术。(4)(4)(4)(4)自动历史拟合。自动历史拟合。自动历史拟合。自动历史拟合。(5)(5)(5)(5)并行算法。并行算法。并行算法。并行算法。(6)(6)(6)(6)一体化的混合模型。一体化的混合模型。一体化的混合模型。一体化的混合模型。(7)(7)(7)(7)研究非线性油藏和非线性流体的数值模拟方法。研究非线性油藏和非线性流体的数值模拟方法。研究非线性油藏和非线性流体的数值模拟方法。研究非线性油藏和非线性流体的数值模拟方法。(8)(8)(8)(8)向着油藏管理模拟系统方向发展。向着油藏管理模拟系统方向发展。向着油藏管理模拟系统方向发展。向着油藏管理模拟系统方向发展。油藏数值模拟技术的油藏数值模拟技术的 发展趋势：发展趋势：(1)新的网格化技术。油藏数值模拟技术的发展123一、油藏数值模拟技术一、油藏数值模拟技术二、油藏数值模拟的主流软件及发展方向二、油藏数值模拟的主流软件及发展方向三、实例三、实例一、油藏数值模拟技术油藏数值模拟技术与应用124现象现象：局部压力异常。：局部压力异常。目标目标：降低高压区，提高低压区。：降低高压区，提高低压区。方法方法：要实现上述目标，可通过下述一个或几个方法来完成：要实现上述目标，可通过下述一个或几个方法来完成：1 1、改变岩石渗透率，使流体从高压带流向低压带。、改变岩石渗透率，使流体从高压带流向低压带。2 2、减少高压带得储油量：、减少高压带得储油量：3 3、通过下述的任意方法增加低压带的储油量：、通过下述的任意方法增加低压带的储油量：示例示例示例示例1 1：现象：局部压力异常。示例1：125方法：方法：最可行的方法是使流体从高压区穿过油层流到低压区。渗透率修改前后对比图压力场修改前后对比图：方法：最可行的方法是使流体从高压区穿过油层流到低压区。渗透率126示例示例示例示例2 2：通过给孔隙度减小一定的比值，比如0.96：孔隙度新=孔隙度原始*比值,由图可见，整体压力水平降低了。整体压力水平降低了。现象现象：整体压力水平太高整体压力水平太高整体压力水平太高整体压力水平太高。示例2：通过给孔隙度减小一定的比值，比如0.96：现象：整体127示例示例示例示例3 3：现象现象：甚至在大量采出流体以后，模型中压力降不明：甚至在大量采出流体以后，模型中压力降不明显。显。原因原因：岩石压缩性较大，因而使饱和度变化的影响可：岩石压缩性较大，因而使饱和度变化的影响可以忽略不记，从有效压缩性方程就可看出：以忽略不记，从有效压缩性方程就可看出：Ce=(SoCo+SwCw+SgCg+Cj)/So处理处理：调整岩石的压缩系数：调整岩石的压缩系数示例3：现象：甚至在大量采出流体以后，模型中压力降不明显。128示例示例示例示例4 4：现象现象：没有任何注入的和流入的水，模型内含水饱和：没有任何注入的和流入的水，模型内含水饱和度出现增加。度出现增加。原因原因：岩石压缩性太低，引起孔隙空间自由体积的扩：岩石压缩性太低，引起孔隙空间自由体积的扩大。在模型中，这种自由体积充填着非流动相，通常大。在模型中，这种自由体积充填着非流动相，通常为水。因此