1、第一章 油藏数值模拟方法分析1.1 油藏数值模拟1.1.1 油藏数值模拟简述油藏数值模拟是根据油气藏地质及开发实际情况，通过建立描述油气藏中流体渗流规律的数学模型，并利用计算机求得数值解来研究其运动变化规律。其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。其基础理论是基于达西渗流定律。油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来模拟实际的油田开采的一个过程。基本原理是把生产或注人动态作为确定值，通过调整模型的不确定因素使计算的确定值(生产动态) 与实际吻合。其数学模型，是通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。充分考虑了油藏
2、构造形态、断层位置、油砂体分布、油藏孔隙度、渗透率、饱和度和流体 PVT 性质的变化等因素。这组流动方程组由运动方程、状态方程和连续方程所组成。油藏数值模拟是以应用数学模型为基础的用来再现油田实际生产动态的过程。具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学，借助大型计算机为介质条件建立三维底层模型参数场中，对数学方程求解重现油田生产历史，解决实际问题。油藏数值模拟技术从 50 年代的提出到 90 年代间历经 40 年的发展，日益成熟。现在进入另外一个发展周期。近十年油藏数值模拟为油田开发研究和解决实际决策问题提供强有力的支持。在油田开发好坏的衡量、投资预测及油田开发方案的优选、
3、评价采收指标等应用非常广泛。油藏数值模拟功能包括两大部分：复杂渗流力学研究，实际油气藏开发过程整体模拟研究，且可重复、周期短、费用低。流体的 P V T 数据 、 相渗曲线 、 岩石数据地质静态参数网格数据化油水井产量 、 井史数据建立地质模型建立网格参数场表格数据 动态模拟初始化拟合生产历史拟合方案预测运算结果输出及分析地质储量拟合区块 、 单井压力拟合含油边界拟合非井点地质静态参数拟合区块 、 单井压力拟合生产指数拟合图 1 油藏数值模拟流程图1.1.2 油藏数值模拟的类型油藏数值模拟类型的划分方法有多种，划分时最常用的标准是油藏类型、需要模拟的油藏流体类型和目标油藏中发生的开采过程，也可
4、以根据油气藏特性及开发时需要处理的各种各样的复杂问题而设定，油气藏特性和油气性质不同，选择的模型也不同，还可以根据油藏数值模拟模型所使用的坐标系、空间维数和相态数来划分。以油藏和流体类型来划分，其模型有：气体模型、黑油模型和组分模型；以开采过程来划分，其模型包括：常规油藏、化学驱、热采和混合驱模型。以油藏和流体描述为基础的油藏模型分为两类：黑油模型和组分模型。（1 ）黑油模型，是常规油田开发应用的油藏数值模型，用于开采过程中，对油藏流体组分变化不敏感的情况，是最完善、最成熟的。黑油模型假设质量转移完全取决于压力变化，适应于油质比较重的油藏类型，在这些模型中，流体性质 Bo、B g、R s 决定
5、 PVT 的变化，如普通稠油及中质油的油气藏。（2 ）组分模型，应用于开采过程中对组分变化敏感的情况。这些情况包括：挥发性油藏和凝析气藏的一次衰竭采油阶段，以及压力保持阶段。同时，多次接触混相过程通常也采用组分模型进行模拟。在组分模型中，适用于油质比较轻、气体组分比较高的油气藏，使用三次状态方程表示 PVT 变化，如轻质油或凝析气藏。（3 ）根据一些特殊开采方式的需要而形成的其他类型的数值模型，如热采模型、注聚合物驱油模型、化学驱油模型、裂缝模型等。其中：热采模型使用了物质平衡方程和能量守恒方程，也用到了组分模型的方法；化学驱油模型与黑油模型有所不同，其附加的守恒方程是用来表示驱油工程中的某种
6、化学用剂的变化。这些模型都是以黑油模型或组分模型为基础演变而来的，采用了最灵活的假设条件和计算方法，只不过在编程过程中加入了一些与特殊开采相应的方程。（4 ）按照研究的需要和地质模型的维数来区分，一般油藏数值模拟模型可分为一维的平面和垂直模型，二维的平面x,y和剖面x,z 模型、柱坐标y,z中二维 2D 的锥进模型，直角坐标系X,Y,Z 中的三维 3D 模型。一维模型可用于有关实验室岩心驱替的模型，二维的平面模型主要研究平面上油水运动规律，不考虑层间的影响；二维的剖面模型主要用来研究层间非均质的影响；二维的锥进模型主要用来研究油井不同射开井段、不同开采强度对气、水锥进的影响；三维模型主要研究一
7、个井组、一个区块甚至整个油藏的开发动态和开发指标。1.2 油藏数值模拟方法油藏数值模拟方法是利用计算机技术模拟地下油气藏开采、驱替的过程，是石油地质人员科学认识、评价油藏的重要技术手段。例如，中石油公司进行的前处理的地质建模工作、清华大学核研院研发的油藏数值模拟管理平台（PNSMP ） 、大庆油田有限责任公司勘探开发研究院研发的 VIP 和 Simbest 格式数据文件相互转换的程序等。油、气、水三相流广泛存在于石油工业中，对于三相流的测量具有重要的意义。现在，油藏数值模拟方法已可用于解决大量的复杂油藏工程问题。如砂岩油藏中考虑油层中各种非均质变化以及重力、毛管力、弹性力等各种作用力的三维三相
8、多井系统的渗流问题，考虑多相、多组分间相平衡关系和传质现象的多相、多组分三维渗流问题等。油藏数值模拟方法不仅在理论上用于探讨各种复杂渗流问题的规律和机理，而且普遍用于开发设计、动态预测、油层参数识别、工程技术问题的优化设计以及重大开发技术政策的研究等。油藏数值模拟方法已经普遍应用于各种油气藏开发过程，成为油气田开发不可或缺的方法和工具，被称作“现代油藏工程” 。跟其它油藏研究方法相比，有着不可替代的优势。油藏数值模拟方法的局限性主要在于：模拟误差，结果不唯一。误差主要来自两方面，一是模型本身有误差，二是油气藏资料不全或不准。油藏数值模拟理论是利用计算机模拟地下油气水的流动，给出某一时刻油气水分
9、布，预测油藏动态。1.3 常用求解方法数学模型建立后，线性方程组的求解是油藏数值模拟方法中最核心的步骤之一，而对模型进行数值求解的第一步是偏微分方程离散化，使之产生线性的或非线性的代数方程组，方程组的线性或非线性是由问题本身的性质以及有限差分近似的性质（系数的显式或隐式处理） ，以便进行数值计算。其方法有：有限差分法和有限单元法。求解线性代数方程组所用的方法有直接法和迭代法两大类，直接法常用的有高斯消元法、高斯-约当降阶法、Crout 分解法、主元素法、D4 方法等；迭代法常用的有交替方向隐式方法、超松弛方法、强隐式方法等。线性的偏微分方程式或者方程组可以直接求解。但油藏模拟中的多相渗流方程组
10、常常是非线性的，即使通过有限差分近似法得到的是个非线性方程组，也可以通过线性化方法来将其转化为线性形式，或者还可以用某种迭代的方法进行求解。在油藏数值模拟中，常用的线性化方法有显式方法、半隐式方法、全隐式方法等。当前最常用的两种求解方法是 IMPES 方法和 Newton-Raphson 方法。在 Newton-Raph-son 方法中，流动方程的有限差分形式中的各项展开成当前迭代级的各项之和，再加上一项在迭代过程中与初始未知变量有关的各项的变化量。为了计算这些变化量，必须计算方程中各项的导数数值解或解析解。这些导数存储在加速矩阵或 Jacobian 矩阵中。第二章 油藏数值模拟技术发展趋势2
11、.1 油藏数值模拟发展现状并行算法并行算法是一些可同时执行的诸进程的集合,这些进程互相作用和协调动作从而达到给定问题的求解。并行算法首先需合理地划分模块, 其次要保证对各模块的正确计算, 再次为各模块间通讯安排合理的结构, 最后保证各模块计算的综合效果。并行机及并行软件的开发和应用将极大地提高运算速度, 以满足网格节点不断增多的油藏数值模型。在并行计算机上使用并行数值解法是提高求解偏微分方程的计算速度, 缩短计算时间的一个重要途径。在共享内存的并行机上把一个按向量处理的通用油藏模拟器改写成并行处理是容易的, 但硬件扩充难; 分布内存并行机编程较共享式并行机困难, 但硬件扩充容易, 关键是搞好超
12、大型线形代数方程组求解的并行化。并行部分包括输入输出、节点物性、构造矩阵、节点流动及井筒等。网格技术为了模拟各种复杂的油藏、砂体边界或断层,渗透率在垂向或水平方向的各向异性, 以及近井地区的高速、高压力梯度的渗流状态, 近年来在国外普遍发展了各种类型的局部网格加密及灵巧的网格技术。这种系统大体可以分为二类: 一类称控制体积有限元网格( CVFE) , 这是将油藏按一定规则剖分为若干个三角形以后, 把三角形的中心和各边的中点连接起来所形成的网格。另一类则称垂直等分线排比网格( PEBI) , 其剖分方法是将油藏分成若干三角形后, 使三角形各边的垂直等分线相交而形成网格。这些方法在处理复杂几何形状
13、油藏及进行局部网格加密时简单而一致。在多相流情况下, 参照某一给定的几何准则时该方法是单调的, 这保证了其稳定性和收敛性。这两种方法都能以直观的控制体积的概念出发并且采用一致的上游权而推导得出。这些方法对网格的方向不敏感, 在某些情况下比九点差分格式的效果好。计算机辅助历史拟合技术斯伦贝谢公司的 Eclipse 数模软件最新推出计算机辅助历史拟合模块( Simopt) 。运用均方差、海赛( Hessian) 矩阵、相关性矩阵、协方差矩阵对结果进行分析以确定敏感参数; 引入梯度带分析技术对地质模型进行优化; 在进行常规历史拟合后, 应用置信度限制( 规定需优化的参数及参数的可调范围) , 通过线
14、性预测分析, 实现计算机辅助调整参数, 减少模拟次数。网格粗化技术对于一些油藏参数( 如孔隙度、深度、饱和度等) , 采用体积加权平均法 ; 对于与流体有关的参数( 如渗透率等) 就不能用简单的加权平均计算得出,而要基于流动计算再进行粗化。流动算法相对精确, 首先解出沿压力降方向的总流量 , 然后再解相同的流动方程, 从而解出等效渗透率。在垂向分层合并计算中, 把相同性质的油砂体( 按相同的物性、储量类型 ) 的网格单元合并在一起,使油藏的数值模拟的网格系统反映出地质沉积特点。网格合并可以按不同井组、区块进行合并计算, 为井组模型和分区模拟提供数据模型。模拟还可以按不均匀网格, 考虑水平方向非
15、均质性及储量分布程度因素等进行内插计算, 提供不均匀网格模型。动态地质建模动态地质建模是壳牌公司的 Kortekass 概括了当前世界上关于油藏地质建模的经验, 提出的建立动态、集成化油藏模型的新概念和技术方法。其强调把动态资料以至数值模拟技术等应用于油藏建模, 从而使所建立的地质模型更加符合油藏的实际情况, 并且要随着油田开发中资料的增多和新资料的获得而不断更新。这种新方法包括一系列获得和运用各种所需资料的技术和方法, 包括地质、地质统计、地震、测井、岩心和流体分析、试井、驱替特征以及网格的细分和粗化, 拟函数的应用等, 但关键是使所建立的地质模型更加符合油藏的实际情况, 而且还可以加快建模
16、的过程。分阶段模拟对开发历史较长、地下储层物性和原油物性发生较大变化的油藏, 把随开发时间变化的地质静态模型划分为多个不同开发阶段的地质模型。常规的油藏数值模拟是从一个油藏( 区块) 投入开发时开始模拟, 一直拟合到目前状况, 再进行方案预测优选。我国许多老油田已进入了高含水或特高含水期。由于开发历史长、综合调整、措施次数多,地下岩石和流体的物性发生了较大的变化, 这给常规模拟工作带来了极大的困难。一方面是历史拟合计算一次所需要的机时非常多, 另一方面是常规模拟无法考虑流体和岩石随时间的变化。因此, 模拟结果的可信度会大大降低。分阶段模拟就是一种解决上述问题的行之有效的方法。分阶段模拟可将一个
17、长期开发的油藏, 按照一定原则划分成几个模拟阶段。动态跟踪模拟油田开发是一个长期过程, 储层物性和原油物性随开发期的不同以及油水井措施发生变化。根据开发期及措施类型制定数值模拟的时间步长, 在油水井见效初期采用较小的时间步长, 进入见效稳定期后以较大的时间步长, 将分析周期由常规的以年计算提高到以月或天计算。对方案实施后的效果、生产状况等再进行跟踪模拟, 并提出新的方案。如此反复研究, 使人们对油藏的构造、物性、油水状况及生产动态的认识更趋合理。数值解法古老的差分法继续得到创造性的发展和应用,也得到深入的分析和研究。差分法的另一重大发展就是全隐式、自适应隐式方法。由于全隐式方法对所有方程系数进
18、行隐式处理, 所以与 IMPES、半隐式和 SEQ 相比, 稳定性好, 隐式程度高, 适应范围更宽。它能解决油、气、水三相渗流、注气、水气锥进、高速气渗等强非线性渗流问题。但是全隐式模型的求解需要采用牛顿迭代法, 因此, 其工作量和存储量比 IMPES、半隐式、SEQ 大。所以, 20 世纪 80 年代初期美国的托马斯等人为了解决方程隐式程度高低和计算量大小之间的矛盾, 提出了自适应隐式方法, 其特点是对不同的网格节点和不同的时间步采用不同的隐式程度来处理, 以便在具有同样稳定性的前提下减小计算量、加快计算速度。1985 年 Berteger W I 等人又提出了一种近似的自适应隐式方法, 其
19、稳定性进一步增强。差分法的又一个重大发展就是多重网格法和预处理共轭梯度法。两者的共同点都在于加速数值解的收敛性。多重网格法实质上是外推与内插技术的创造性应用。近年来, 共轭梯度法与各种不完全分解预处理相结合, 并采用 D4 网格排序及 Orthomin 加速技术, 从而使预处理共轭梯度法成为 20 世纪 80 年代油藏数值模拟中最引人注目的方法之一。预处理共轭梯度法克服了油藏数值模拟中由于局部网格加密、混合网格嵌套、隐式井底压力处理、大断层大裂缝处理、死节点处理以及自适应隐式方法等产生的不规则系数矩阵的求解问题。由于 Ort homin 方法收敛性好, 且运算过程简单、不需要迭代参数、不需要估
20、计矩阵的特征值, 是一个稳定、有效的迭代加速方法。因此, 目前该法在油藏数值模拟中获得了相当成功的运用。三维动态显示数值模拟结果的可视化程度高, 人机交互性能强。可三维动画显示油藏中流体的流动规律, 再现油藏的开发历史及剩余油的空间分布, 并可任意旋转、平移、缩放、光照, 多重照相、透视和透明处理; 灵活的剖面切割功能 , 任意参数的区间显示, 用户可按自己的爱好定义或修改颜色和注释。2.2 油藏数值模拟技术发展趋势工作站由于价格昂贵以至很难普及。而今 PC 机越来越强大, 速度日益加快, 内存越来越大且易扩展, 价格便宜, 应用广泛, 因而将数模软件开发在 PC 机上的远景极为广阔。但目前的
21、困难是 PC 并行机尚没有统一的标准且易出现网络问题。并行算法的进一步完善、各种灵巧网格的随意化、自适应化以及历史拟合的自动化将大幅度提高数值模拟的精度, 减少模拟计算的工作量。在大规模粗化前利用流管法, 预先筛选地质模型。根据物质平衡方程和达西定律, 在单时间步内采用隐压显饱解法计算流线上流体的流动: 全隐式有限差分压力, 沿流线解饱和度方程。在考虑重力影响时, 先沿流线再沿重力线求解压力和饱和度。勘探开发一体化越来越受到关注。油气勘探开发的各个环节中计算机应用技术的突飞猛进, 带动了勘探开发工作向着多学科综合集成方向发展, 改变了以往以单项专业为主的顺序化工作方式。各学科间的信息交换、数据
22、共享、成果可视和知识继承遵照统一的标准进行综合集成: 通过岩心声波测量建立油藏参数与地震响应间的关系; 利用油藏数值模拟技术预测压力和含水饱和度的变化、生成不同观测数据的合成地震记录; 开发地震包括了四维地震和地震反演技术, 通过井的约束以及和油藏数模技术的结合, 将进一步减少多解性的困扰并大幅度提高分辨率, 从而扩展开发地震的实际应用范围 ,特别对比较复杂的砂泥岩薄互层也能够加以应用。这些技术的发展, 将大大加深对井间储层非均质特征及油藏剩余油分布状况的认识。但当前面临的困难是: 种类各异的数据格式、各式各样的数据库系统、互不相通的应用软件以及不同的用户界面等严重妨碍了多学科的协同工作、数据
23、资源及研究成果的共享。借助多井试井数据、数值模拟和优化技术, 通过对试井实测压力数据的拟合, 反演油藏中某些非均质参数分布, 能有效地解释测试井的径向平面渗透率变化情况。借助于多井试井数据, 结合油藏的地质信息, 将突破目前建立在单层、单相、均质储层等基本假设条件下的局限, 为非均质、多层、多相等复杂问题提供满意的解释。第三章 断块油藏数值模拟技术3.1 断块油藏地质及开发特征3.1.1 地质特征（1 ）断块复杂，断块小而多。复杂断块油藏的形成主要来源 于众多的断层切割作用,断层多是这类油田共同的、最突出的地质特征。一个较大的复杂断块油田,往往有几十条甚至上百条的断层分布,由于油藏被众多断层切
24、割, 结果形成许多被切割的不连通小断块单元油藏。（2 ）含油层系多，常具有多套油水系统，油水关系复杂。复杂断块油藏大多含有多套含油层系。在这类油藏中,凡有断裂系统到达的层位,只要有储集岩发育,差不多都有油层发育。复杂断块油藏内含油层系虽多,但富集程度差别较大。通常只有一套或少数几套分布广泛,油层厚度相对较大,富集程度相对较大的主力含油层系,是小断块油 藏单元 的最主要的开发层系。（3 ）储层变化大，物性较差，层间非均质严重。在复杂断块油藏中,一般没有一个层位能够在全油藏规模上连片含油。就单个含油层系来说, 常常只在油藏的一部分面积上含油,甚至只在很小的一部分面积上含油。纵向上,砂体类型多、变化
25、大、连通差、非均质严重。（4 ）原油性质变化大，生产动态差异大。大多数复杂断块油田内部原油性质变化幅度大。常引起油层产能和动态的差别。此外,产能和动态特点还受到储层物性和其它地质条件的影响,复杂断块油藏中含油层系多,每套层系的岩石性质,沉积环境和成岩作用 又各不相同。即使同一套含油层系,由于不 同断块的埋深不同,物性也会有很大的差异。因此在一个复杂断块油藏 内部常常存在物性差别很大的储层, 从而引 起产能和 生 产动态出现明显差别。3.1.2 开发特征（1 ）开发单元小而多，单元间差异大。复杂断块油藏 的大多数断块 含油面积都很小,而且同一个断块里各个油层位置常常相互错开, 叠合在一起能 自然
26、地组合成一个开发单元的储量少,所以开发单元一般都很小。由于开发单元小, 所以与同样规模的油藏来比较,复杂 断块油藏的开发单元数就比大油藏多得多。每个开发单元的动态特征不一 样,针对性的有效开发措施也不一样。（2 ）初期产量高，产量递减快，很快就会进入中后期能量低产生产。由于开发单元小,油井初期的较高产量只是短暂的, 随后产量迅速降低,当能量消耗较多时,就只有维持低能量低产生产。断块油藏控制单元井间连通性差、边界封闭、地层能量有限，开采中地层能量下降很快，难以采用早起注水等方式大规模开采。用 CO2 吞吐强化开采是有效的方法。3.2 二氧化碳吞吐机理3.2.1 二氧化碳机理二氧化碳吞吐机理是通过
27、与地层原油之间的多次接触非混相过程达到驱替和携带原油的目的，其机理包括以下几个方面：二氧化碳溶解引起原油体积膨胀。大量室内和现场试验表明,原油中充分溶解 CO2 后可使原油体积膨胀 10 %40 % ,注入CO2 后原油的体积增加,其结果不仅增加了原油的内动能, 而且也大大减少了原油流动过程中的毛管阻力和流动阻力,从而提高了原油的流动能力。二氧化碳溶解降低原油粘度。当原油中的 CO2 溶解气饱和后,能够大大降低原油的粘度。在地层条件下,压力越高,CO 2在原油中的溶解度就越高,原油的粘度降低越显著。二氧化碳溶解气具有弹性驱动作用。油层中的 CO2 溶解气,在井下随着温度的升高部分游离汽化,以压
28、能的形式储存部分能量。当油层压力降低时,大量的 CO2 则从原油中游离, 将原油驱入井筒, 起到溶解气驱的作用,由于气体具有较高的运移速度,从而将油层堵塞物返吐出来。二氧化碳时原油中轻质烃萃取和汽化轻质烃与 CO2 间具有很好的互溶性,当压力超过一定的值( 此值与原油性质及温度有关) 时, CO2 能使原油中的轻质烃萃取和气化, 这种现象对轻质原油表现得尤为突出。CO 2 对原油中的轻质烃的萃取和汽化现象是注入 CO2 增油的主要机理之一。改善原油和水的流度比由于大量注入的 CO2 在原油和水中的溶解,地层水的碳酸化,使原油流度增加, 而水的流度降低, 从而使原油和水的流度趋于接近, 使水的驱
29、油能力提高,同时也进一步扩大了水驱的波及面积,大大提高了扫油效率。酸化解堵。CO2 溶解于水后略呈酸性,与地层基质相应地发生反应, 从而酸解一部分杂质, 尤其在碳酸盐岩中能将部分岩石溶解,生成易溶于水的碳酸氢盐, 从而提高碳酸盐岩层的渗透性。由于注入 CO2 气的酸化作用导致油层渗透性提高,在一定压差下,一部分游离气对油层的堵塞物具有较强的冲刷作用, 可有效地疏通因二次污染造成的地层堵塞.3.2.2 CO2 吞吐数值模拟3.3 数值模拟3.3.1 建立油藏数值模型在油藏地质模型研究的基础上，从 PETREL 建模软件中直接输出粗化模型，纵向上以单砂体为单元， ，其它高压物性数据、相渗数据等直接
30、应用该断块该层系实验数据，时间单元以月为单位，确保模拟时能满足因生产中层位和频繁作业对时间步的要求。3.3.2 历史拟合历史拟合的主要目标是完善和验证油藏模拟模型，初始模拟输入数据在一定程度上不能拟合油藏历史动态，以达到精确预测未来动态的目的。历史拟合过程总是使对油藏过程有更好的理解。含水面支撑的平面，流体迁移的路径，以及原油绕流的区域都能在历史拟合过程中得到认识。通过历史拟合，一方面可以求解油藏参数，另一方面可以对油藏动态进行预测，为开发方案的调整提供依据。（1 ）历史拟合的指标历史拟合的指标主要是那些能反应油藏及单井动态的指标。油藏数值模拟的指标包括：油藏原始地质储量、油藏累计产油量、累计
31、生产气油比、综合含水率、地层平均压力、累计注水量等随是火箭的变化。油井动态指标包括：瞬时产油量、含水率、产水量、瞬时气油比、累计产油量、累计产水量、井底压力等随时间的变化。水井动态指标包括：瞬时注水量、累计注水量、井底压力随时间的变化。（2 ）历史拟合参数调整范围历史拟合的过程中,由于模型参数数量多, 可调自由度大。而实际油藏动态数据的种类有限,不是可以唯一确定油藏模型的?数,为了避免修改参数的随意油层厚度 :经过多次的划分和认定, 一般情况下油层厚度误差不大, 可以视为确定参数,但由于软件插值带来一些误差, 也允许做一些调整。孔隙度:根据油藏的単井孔隙度统计结果,其值变化范围比较大,在拟合过
32、程中允许个别调整。渗透率:渗透率在任何油田都是不确定参数,这不仅因为解释和岩芯分析的误差较大,而且根据渗透率的特点,井间渗透率分布也是不确定的, 另外由于速敏、水敏、酸敏、盐敏等影响,开发过程中出现增渗或降渗现象, 高渗层渗透率越来越大 ,而低渗透层渗透率越来越小, 级差进一步加大。因此,对渗透率的修改可调范围相对较大。岩石和流体压缩系数:它们一般变化范围很小,可以作为确定参数处理,但在一定范围内,也可以适当做些微调。相对渗透率曲线:由于油藏模拟网格较粗,网格内部存在着非均质,其影响不容忽视, 这与均质的岩芯不同, 因此, 相对渗透率应作为不定参数。有效厚度：主要是由于部分小夹层没有被完全分离
33、出来，所以有效厚度的调整为0% 30%。油水油气界面：资料不多的情况下，允许在一定范围内修改。油、气、水的高压物性资料：是实验室测定的，为确定性参数。一般认为是准确的，在一定情况下也可以小范围的修改，但需要在各种资料对比的基础上进行修改。（3）历史拟合参数调整方法油藏是一个复杂的系统，各个参数的调整都是相互影响的，因此要拟合一个指标，首先要知道哪些参数会影响到这个指标，然后又正对性的调整。1）油藏储量的拟合数值模拟方法计算储量和油藏精细描述计算地质储量都是采用容积法，影响储量的因素有含油面积、有效厚度、孔隙度、原始含油饱和度、地面原油粘度、原油体积系数等。含有面积的差异是在建立模型时差值条件不
34、同，导致砂体厚度差值时出现油藏范围扩大或者缩小的现象，这时需要对模型边界进行控制，添加有效厚度零线，或者修改断层边界。有效厚度的修改，应该求出地质模型中每个小层的平均有效厚度，并将其与油藏精细描述的结果进行对比，以确定是否进行修改。修改时建议对有效厚度进行整体的校正，而不应该从单井厚度方面进行修改。地层孔隙度修正方法是：首先求出各个小层的平均孔隙度，然后进行对比以确定是否需要修改。原始含油饱和度容易出现差异的地方是油水过渡带，对于平衡初始化，可以通过修正毛管压力曲线来校正油水过渡带的油水分布，对于非平均初始化的油水分布，可以手工进行修改。地面原油密度和原油体积系数是通过高压物性实验获得的，一般
35、比较准确，最好不要进行修改，如果需要修改，需要提出确实充分的证据。2）油层压力的拟合在数值模拟过程中，需要拟合的压力包括：油层平均压力、压力分布、井底流压。影响整体地层压力水平的因素主要有：综合压缩系数、边界性质、产气量、劈产数据的准确性。3）油藏综合含水率（气油比）的拟合含水率和气油比的拟合主要反应了油水和油气之间的分流能力，通过调整相渗曲线来实现。含水率的拟合包括两部分：拟合初始含水阶段的含水率，主要是无水采油期，拟合综合含水率变化趋势。4）单井含水率的拟合影响单井含水率的主要因素有：初始状况下井底周围的流体分布、周围注水井注入量、层间的非均质性、井底丢东参数等。 、单井含水率的拟合需要拟
36、合初始的含水率和整个含水率的变化趋势，初始含水率的大小主要取决于单井周围初始的流体饱和度和井所在区域采用的相对渗透率曲线。5）单井井底压力的拟合模拟井底压力的大小主要取决于两个方面：意识井底周围地层的性质，二是相邻注水井的注水对应状况。井底周围的地层性质主要指地层的渗透率。生产井井底的渗透率越小，流动耗费的能量越大，相应的井底压力越低；反之，井底压力越高。根据计算的井底压力的高低趋势调整井所在网格的地层渗透率，如果计算的井底压力偏高，则可以适当降低渗透率；反之，应该增加网格的渗透率。除了地层渗透率的影响外井底的污染状况也显著影响井底压力的计算。这时可以靠修改表皮系数来达到拟合井底压力的目的。在
37、定液量寸产的条件下，增加表皮系数会降低井底的流压，减小表皮系数会增加井底的流压。针对复杂断块油藏此次油藏数值模拟拟合原则有三个方面：第一，逐井逐层分析注采关系及液流方向进行重点拟合；第二，对开采时间长、注采量大的井进行重点拟合；第三，重点对高含水开发后期进行历史拟合。3.3.3 研究剩余油分布运用数值模拟软件直接计算出的小层含油饱和度值即为开发至目前地下重新分布的剩余油。第一类为注采井网不完善造成的剩余油; 第二类为层间、层内干扰造成的差层剩余油; 第三类为封闭性断层附近造成的剩余油; 第四类为剩余油分布零散型。1.立足的尺度域大不同尺度范围内的油水分布结果对油田开发的指导意义是不相同的。数值
38、模拟可以考虑 一个区块内的流体分布特征，其尺度范围在 km 级以上，而且可以得到油藏范围内的流体分布特征，尤其体现在油藏注水开发后期寻找平面和纵向上的剩余油潜力时。在纵向上，数值模拟结果可以提供砂层组、小层，甚至沉积时间单元内的流体分布特征，有助于分析层间矛盾，给开发调整提供可靠的依据。在平面上，数值模拟可以提供井间的油水分布特征。2.考虑的地层条件多一般计算油水分布的方法都在计算过程中简化了多种因素，但是数值模拟却可以考虑油藏中的复杂地质特征和渗流环境。影响剩余油分布的因素很多简中来说可以分为地质因素和开发因素。在这两者中，先天的地压条件是主号因素，开发因素是依据先天因素改变而发生变化的。地
39、质因素上要包括油藏的非均质性、构造特征、断层的分布以及流体的性质处异等；开发因素主要包括注采系统的完善性、井网方式、采取的措施类型以及采油方式等。这两大类的影响因素控制了地层中的油水分布，因此形成的剩余油形式和数员也合较大差界。从地质特征的角度来看，剩余油主要形成在以下的区域：断层及油层边角地带的滞留区、构造高部位及正向微构造区、储集层物性较差的部位、原油粘度较高的区域等。由于地历特征反映了油藏的本质特征而是这些特征一般不会随开发的过程有所改变，受到的人为因素的影响比较小，因此成为影响剩余油分市的主要因素。在高含水开发期，如果找到这种因素控制的剩余油，往往采用加密井网的开发效果比较好。3.提供
40、烦人计算结果的信息量大通过数值模拟计算，可以提供地层中的油、气、水在不同开发时刻的饱和度数值和分布特征，而且还可以计算出地层中的剩余储量分布、剩余储量丰度、剩余可采储量丰度的分布规律，用于指导开发后期的调整。通过对计算结果的整理分析，还可以得到地层中的含水率分布特征，该特征不同于现场上采用油藏工程方法得到的含水率分布特征，它能比较真实地反映地层中的流体分布。此外，还可以计算纵向上生产井的产液剖面、含水剖面，从而有效地指导堵水工作。对于三次采油的数值模拟计算，还可以得到地层中各种化学剂的分布规律，监测、计算驱替前缘的饱和度变化。同样，也可以得到地层中流体特征随开发过程的变化。井网完善程度、单井控制储量程度都可以通过数值模拟的结果计算出来。4.结果的可验性对于数值模拟方法计算的结果，可以利用不同时期、不同形式的资料给予验证。例如生产测试的产液剖面和吸水剖面可以对计算的单井剖面进行验证。此外，还可以利用已经有的生产资料对计算出的动态结果进行对比和分析。计算结果的可验证性还表现在利用数值模拟可以计算不同情况下的油藏动态并对已经计算的结果进行对比分析，从中优选出比较合理的方案和措施。
