化学驱油藏数模并行化中的关键技术
油藏数值模拟方法
第一章油藏数值模拟方法分析油藏数值模拟油藏数值模拟简述油藏数值模拟是根据油气藏地质及开发实际情况，通过建立描述油气藏中流体渗流规律的数学模型，并利用计算机求得数值解来研究其运动变化规律。
其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。
其基础理论是基于达西渗流定律。
油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来模拟实际的油田开采的一个过程。
基本原理是把生产或注人动态作为确定值，通过调整模型的不确定因素使计算的确定值（生产动态）与实际吻合。
其数学模型，是通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。
充分考虑了油藏构造形态、断层位置、油砂体分布、油藏孔隙度、渗透率、饱和度和流体PVT性质的变化等因素。
这组流动方程组由运动方程、状态方程和连续方程所组成。
油藏数值模拟是以应用数学模型为基础的用来再现油田实际生产动态的过程。
具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学，借助大型计算机为介质条件建立三维底层模型参数场中，对数学方程求解重现油田生产历史，解决实际问题。
油藏数值模拟技术从50年代的提出到90年代间历经40年的发展，日益成熟。
现在进入另外一个发展周期。
近十年油藏数值模拟为油田开发研究和解决实际决策问题提供强有力的支持。
在油田开发好坏的衡量、投资预测及油田开发方案的优选、评价采收指标等应用非常广泛。
油藏数值模拟功能包括两大部分：①复杂渗流力学研究，②实际油气藏开发过程整体模拟研究，且可重复、周期短、费用低。
图1油藏数值模拟流程图油藏数值模拟的类型油藏数值模拟类型的划分方法有多种， 划分时最常用的标准是油藏类型、需要模拟的油 藏流体类型和目标油藏中发生的开采过程， 也可以根据油气藏特性及开发时需要处理的各种各样的复杂问题而设定， 油气藏特性和油气性质不同， 选择的模型也不同， 还可以根据油藏 数值模拟模型所使用的坐标系、空间维数和相态数来划分。
二氧化碳驱技术在低渗透油藏开发中提高驱油效率的研究与应用
二氧化碳驱技术在低渗透油藏开发中提高驱油效率的研究与应用摘要：在中石化总公司支持下，组建了CO2驱技术研究团队，形成了高温高盐油藏CO2驱油三次采油关键技术，解决水驱废弃油藏和低渗难动用储量的开发难题。
在国内率先开展了特高含水油藏CO2/水交替驱；深层低渗油藏CO2驱。
油田层次开展了四种油藏类型五种矿场试验。
验证该类油藏二氧化碳驱可行性，探索合理举升方式，进一步优化二氧化碳驱井网井距，验证大井距可行性，探索深层低渗稠油油藏有效开发方式，扩展二氧化碳驱应用范围以及特高含水废弃油藏二氧化碳驱提高采收率技术。
探索储层粘土含量高、水敏性强油藏二氧化碳驱提高采收率技术。
关键词：二氧化碳驱低渗油藏提高采收率换油率1、研究目的1.1 某厂低渗难动用储量涉及开发单元11个，地质储量1601.85×104t，标定采收率7.56%，目前采出程度5.54%。
涉及单元多为低孔隙、低渗透的地质特点。
2010年开始二氧化碳驱在胡1块深层低渗油藏实施先导试验，胡1井组气驱取得成效后，相继在其他五个低渗类型油藏实施气驱开发。
目前总覆盖地质储量309.5×104t。
累注气17.9×104t，累增油3.05×104t。
1.2低渗油藏水驱效率低，注采井组呈现两极分化现象，一是注水压力高油井难以见效，二是油井见效快、含水上升快、见效稳产周期短，通过二氧化碳驱提高驱油效率。
2、研究内容及成果2.1 二氧化碳驱机理上优于水驱一是超临界二氧化碳注入能力强，增大有效井距；二是CO2驱补充地层能量，可膨胀地层原油，提高驱油效率再者CO2能进入的孔喉半径比水小一个数量级（0.01μm），低渗油藏，增加驱油体积25%以上，随CO2溶解，原油体积膨胀。
毛管半径分布曲线不同驱替方式驱替压力变化曲线2.2二氧化碳驱解决注入压力过高的问题根据深层低渗油藏开发情况调查，注水压力高，注气难度不大。
从地质条件类似的胡某区块二氧化碳注入能力看，二氧化碳驱可以解决注入压力过高的问题。
中国石油勘探开发研究院各专业研究方向
中国石油勘探开发研究院各专业研究方向一、地质资源与地质工程（代码：0818）（一）矿产普查与勘探（代码：081801）1．油气成藏与含油气系统研究以油气藏形成条件与富集主控因素研究为基础，通过油气成藏静态地质要素和动态作用过程的综合分析，揭示油气成藏过程与富集规律。
研究内容包括有效烃源岩、储集层、输导层和盖层等地质要素的分布和静态评价，油气生成、运移、聚集成藏和圈闭的形成等作用的动态演化过程和时空匹配关系，以及关键时刻地质要素和动态作用组合关系，开展油气成藏综合研究与评价，明确油气富集规律，预测油气资源规模和资源空间分布。
2、非常规油气地质学非常规油气地质学是以非常规油气资源类型、细粒沉积体系形成与分布、微纳米级致密储层特征、连续型油气聚集与产出机理、“甜点区”评价方法与技术等为重点的新兴学科。
研究核心是非常规油气成藏体系的“生油气能力、储油气能力、产油气能力”；研究内容包括成藏体系的烃源性、岩性、物性、脆性、含油气性与应力各向异性“6特性”及匹配关系，研究重点是非常规油气成藏体系的分布范围与“甜点区”的分布预测与评价，确定经济有效开采的方法技术与经济发展模式。
3．盆地和构造分析以区域构造背景研究为基础，以地质、地球物理综合研究方法技术为手段，确定含油气盆地成盆演化与沉积充填历史，研究盆地性质、构造样式、类型和分布，明确盆地构造演化对油气成藏的影响。
研究内容包括成盆区域构造背景与构造动力学机制、构造运动学和几何学特征、区域构造演化、盆地构造解析、构造运动对成盆、成烃、成储、成藏的影响。
4．沉积与储层地质沉积研究是充分运用现代沉积学的理论和层序地层学、地震沉积学等研究思路与方法，明确沉积体系类型，研究沉积物的形成、搬运、沉积演化过程，确定沉积环境和沉积相、成岩作用和沉积演化特征，明确有利沉积相带。
储层学研究是以沉积研究为基础，研究储集体的岩性、物性、电性和含油气性特征，揭示与储集空间（孔、洞、缝）形成有关的成岩作用，阐明成岩历史、孔隙演化历史，构建储层地质模型，开展储层分布预测与评价，确定有利储层分布范围。
建模数模一体化技术提升油田开发研究水平-地质院
一、建模数模一体化技术发展现状
1. 建立了“集中配置、统一管理、共享应用”的油田地质综合研 究中心，支撑了油田开发综合研究信息化规模应用
128个CPU， 峰值运算速度 40万亿次/秒
DELL CR620集群
2套
地质建模
万兆交换 机
测井解释
30台
图形工作站
数值模拟 HP CP600集群
支持开发S6U大N服主务要器 业服 务务器及工地作质站制图
拓展培训
➢ 地质人员学物探 ➢ 开发人员学测井 ➢ 研究室与实验室交叉培养
精英培训
➢ 软件高级培训 ➢ 一人一块实践应用
熟练应用人员情况
建模 软件
87名 比例90%
地质人员：97
数模 软件
113名 比例87%
油藏人员：130
基础培训
➢ 技术原理及软件操作培训 ➢ 岗位锻炼
地质院建成一支深化技术研究的专家型队伍。
临时管理节点，提高运行效率。
一、建模数模一体化技术发展现状
1. 建立了“集中配置、统一管理、共享应用”的油田地质综合研 究中心，支撑了油田开发综合研究信息化规模应用
管理节点浮动技术 ：
1800
1600
1400
时 1200 间 1000
800
s 600
400
200
0
0 1
20.5
87
2
3
作业数
1683
100
98
96
94
92 90
89.5
88
86
84 2011年
96.9
2012年
时间
98.8
2013年
方案研究精细化、科学化水平大幅提升，新区产能建设达标率达到98.8%
气液两相流数值模拟方法的研究与应用
气液两相流数值模拟方法的研究与应用气液两相流是指同时存在气体和液体的复杂流动现象，广泛存在于自然界和工业生产中，如瀑布、波浪、化工反应器、石油开采等。
气液两相流的研究对于理解和控制这些现象、提高生产效率和安全性具有重要意义。
数值模拟是研究气液两相流的有效方法。
相比于实验方法，数值模拟的优势在于能够获得更多的细节信息和精确数据，同时也可以极大地降低成本并避免实验过程中的危险性和不确定性。
本文将介绍气液两相流数值模拟的方法，及其应用领域和未来挑战。
一、数值模拟方法1. 传统方法传统方法通常采用两相流模型，基于欧拉方程求解。
由于气液两相流的复杂性，这种方法常常涉及到多个物理场的耦合和相互作用，如热传递、质量传递、化学反应、多相流动力学等。
因此，该方法具有计算量大、计算时间长、计算结果不精确等缺点。
2. 基于LBM的方法LBM（lattice boltzmann method）是一种介观尺度（宏观与微观之间的中间尺度）数值模拟方法，可以直接模拟流体内部微观运动方式，适用于模拟多相流动现象。
这种方法是根据Boltzmann方程建立的，通过碰撞模型模拟流体分子的运动，以此获得整个流场在不同时间的状态。
该方法具有计算速度快、模拟精度高、易于建模及可扩展性等优点。
3. 基于CFD的方法CFD（computational fluid dynamics）是指应用计算机数值方法对流体流动进行模拟和分析的工程技术。
CFD方法通过建立流动场的数学模型并采用数值求解方法进行计算，从而得到流场的物理或数学解。
这种方法在气液两相流领域中也得到了广泛应用。
4. 其他方法此外，还有一些其他的数值模拟方法，例如基于粒子方法的SPH（smoothed particle hydrodynamics）和DEM（discrete element method）等。
这些方法基于不同的假设和算法，都有各自的优缺点，在不同的气液两相流应用场景中发挥着重要的作用。
油藏数值模拟实用技能培训
- 考虑一个典型的注采单元，二维或三维 - 跟踪井组单元内的驱替前缘、驱替效率研究
5.4 全油田模拟
- 考虑参与渗流的整个生产区域，包括油藏、井筒 - 模拟全油藏的整体开发规律
单井模型网格示意图
利用井组概念模型，进 行鱼骨刺分支水平井注 采井网开发效果研究
一能做些简单的科学运算 70年代：计算机快速升级，带动了油藏数值模拟的迅猛发展，计算速度达到
每秒100-500万次，内存增至16兆字节。黑油模型日趋成熟，已能解决中小 型油藏的模拟问题 80年代：油藏数值模似技术飞跃进步，计算方法及软件迅猛发展，超级向量 机诞生，计算机速度达到亿-几十亿次，内存高达10-20亿字节 90年代：油藏模似软件的各模块功能迅速发展，体现为向一体化方面发展， 出现了集地震、测井、钻井工程、完井工程、油藏工程、采油工程、压裂酸 化措施工艺及地面集输、经济评价等为一体的大型软件
[ Hg KKrg x g
( Pg x
g
g
D x
)]
y
[
Hg KKrg g
( Pg y
g
g
D y
)]
qg
H
t
(g
Sg
)
[ HwKKrw x w
( Pw x
wg
D)] x
[ HwKKrw y w
( Pw y
d 1
dr1 dr2 dr3
k=1 k=2 k=3 k=4
径向坐标网格系统
局部加密网格系统
混合网格系统
一、什么是油藏数值模拟 ？
8 模型的类型
8.1 根据流体分类 （1）黑油模型
- 简单的相态（油-气-水），温度的变化可以忽略 - 用于一次采油、注水、干气注入、聚合物注入
化学驱新技术及华北油田化学驱现状与展望
化学驱新技术及华北油田化学驱现状与展望
杜慧丽;游靖;李凤群;张田田;高珊珊;王玥;路淼森;赵紫印;唐雅娟
【期刊名称】《精细石油化工进展》
【年(卷),期】2024(25)2
【摘要】化学驱是大幅提升油田采收率的重要技术。
化学驱理论和技术持续创新,应用规模不断扩大,其中聚合物驱矿场提高采收率约14%,二元复合驱和三元复合驱矿场提高采收率均在18%以上,纳米驱油技术预期采收率可达80%以上,中相微乳液驱表面活性剂体系性能不断突破,矿场取得良好增油效果。
针对华北油田复杂断块砂岩油藏的实际情况,采用以可动凝胶为主的深部调驱技术,进一步提高了油田采收率。
为进一步提高华北油田采收率,可以在深化油藏剩余油分布规律认识的基础上,加强基础理论研究,推动化学驱技术的持续创新和有序接替,配套开发化学驱地面工程,以期实现油田开发的可持续发展。
【总页数】7页(P5-10)
【作者】杜慧丽;游靖;李凤群;张田田;高珊珊;王玥;路淼森;赵紫印;唐雅娟
【作者单位】华北油田公司工程技术研究院;华北油田公司第三采油厂;华北油田公司勘探开发研究院
【正文语种】中文
【中图分类】TE3
【相关文献】
1.红岗油田化学驱提高采收率技术展望
2.油田CO_(2)驱提高采收率技术及现场实践分析——评《化学驱提高石油采收率》
3.海上油田化学驱注入浓度优化调整新技术及矿场实践
4.纳米复合化学驱在低渗透油田中的驱油特性研究
5.抢占世界化学驱油技术“制高点”——记“国家卓越工程师团队”、中国石油大庆油田化学驱油技术研发团队
因版权原因，仅展示原文概要，查看原文内容请购买。
大规模油藏数值模拟软件并行计算技术及在Beowulf系统上的应用进展
２００６年６月数值计算与计算机应用第２期大规模油藏数值模拟软件并行计算技术及在ＢｅＯＷｕｌｆ系统上的应用进展¨’曹建文刘洋孙家昶姚继锋潘峰（中国科学院软件研究所并行计算实验室北京１０００８０）摘要本文主要介绍了大规模油藏数值模拟并行计算技术在国内的研究进展，提供了精细油藏模拟在国产Ｂｅｏｗｕｌｆ系统上的计算实例和应用效果，给出了百万网格点规模的油藏应用算例在不同处理器规模下的数值模拟计算结果与性能分析，并实现了一个针对海量数据可视化的三维图、二维图、表格显示的后处理显示系统．关键词：高性能计算，油藏数值模拟，并行计算，非线性方城求解器，预处理，Ｂｅ０１ｗｕｌｆ系统，可视化ＴＨＥＰＡＲＡＬＬＥＬＣＯＭＰＵＴＩＮＧＴＥＣＨＮＩＱＵＥ０ＦＬＡＲＧＥＳＣＡＬＥＰＥＴＲＯＬＥＵＭＲＥＳＥＲＶＯＩＲＳＩＭＵＬＡＴｌ０ＮＳＯＦＴＷＡＲＥＡＮＤＩＴ’ＳＡＤＶＡＮＣＥＩＮＩＮＤＵＳＴＲＩＡＬＴＥＳＴＣＡＳＥＳ０ＮＢＥＯＷＵＬＦＳＹＳＴＥＭＳＣａｕＤＪｉａｎＷｅｎＬｉｕＹａＬｌｌｇＳｕｎＪｉａｃｈａｎｇＹａｏＪｉｆｂＩｌｇＰａｎＦｅｎｇ（Ｐｏｍｌｆｅｌ∞＂妒珏抚咖三０６Ｄｍ幻阿，血ｓｔ姑钍￡ｅＤ，Ｓ妒伽。
他，Ｃ仡ｉ佗ｅｓｅＡｃｎｄｅ仃ｌ∥。
，Ｓｃ沈佗ｃｅｓ，Ｂｅ巧打姆１０００８０）ＡｂｓｔｒａｃｔＴｈｉｓｐａｐｅｒｒｅｐｏｒｔｓｓｏｍｅａｄｖａｎｃｅ８ｏｆａｐｅｔｒｏｌｅｕｍｒｅｓｅｒｖｏｉｒｓｉＩｎｕｌａ七ｉｏｎ８０ｆｔ—ｗａｒｅｉｎ１ａｒｇｅｓｃａｌｅｐａｒａＬｌｌｅｌｃｏｍｐｕｔｉｎｇ，ａｎｅｎ．ｅｃｔｉ、他ａｎｄｅｍｃｉｅｎｔｓ０１Ｖｅｒｏｆｎｏｎｌｉｎｅａｒｅ（１ｕａｔｉｏｎ８，ｐｒｅｃｏｎｄｉｔｉｏｎｉｎｇ，ａｎｄｉｎｄｕ８ｔｒｉａｌ印ｐｌｉｃａｔｉｏｎｓｏｎＢｅｏｗｕｌｆｓｙｓｔｅｍｓ．ＴｈｒｅｅｉｎｄｕｓｔｒｙｃａｓｅｓａｒｅｕｓｅｄｔｏｔｅｓｔｐｅｒｆｏｒｍａｎｃｅｏｆｔｈｉｓｐａｒａＵｅｌｒｅ８ｅｒｖｏｉｒｓｉＩｎｕｌａ七ｏｒｏｎ８ｅｖｅｒａｌｐａｒａＵｅｌｈ盯ｄｗａｒｅｐｌａ七ｆｂｒｍｓ，ａｎｄｔｈｅａ｜ｎａｌｙｓｉ８ａｎｄｅｖ砒ｕａｔｉｏｎｏｆｔｈｅｒｅ—ｐｏｒｔｅｄｒｕｎｎｉｎｇｒｅｓｕｌｔ８ａｒｅｇｉｖｅｎ．Ｔｈｅｐａｐｅｒａｌｓｏｉｎｔｒｏｄｕｃｅｓｔｈｅｒｅｑｕｉｒｅｍｅｎｔｓｏｆｌｅａｓｔｃｏｍｐｕｔａ七ｉｏｎｃｏｓｔｓｏｆ矗ｎｅｒｙｒｅｓｅｒｖｏｉｒｓｉｍｕｌａ七ｉｏｎ８，ｔｈｅｒｅｑｕｉｒｅｍｅｎｔ８ｏｆｐａｒａｌｌｅｌｃｏｍｐｕｔｅｒｈ村ｄ／ｓｏＲｗ盯ｅｓｏｆａｒｅｓｅｒ、ｒｏｉｒ８ｉＩｎｕｌａｔｉｏｎｓｙｓｔｅｍ，ｃｏｍｐｕｔａｔｉｏｎｍｅｔｈｏｄｓａｎｄｐｒｅｃｏｎｄｉｔｉｏｎｉｎｇｏｆｔｈｅｎｏｎｌｉｎｅａｒｐｒｏｂｌｅｍｓｏｌｖｅｒｏｆｒｅｓｅｒｖｏｉｒｓｉｍｕｌａ七ｉｏｎ，ｔｈｅｉｎ—ｆｌｕｅｎｃｅｓｏｆｐａｒａｌｌｅｌｃｏｍｐｕｔａｔｉｏｎｐｅｒｆｏｒｍａｎｃｅｅ砌ｕａｔｉｏｎｏｆｔｈｅｃｏｍｐｕｔｅｒｈ盯ｄｗ壮ｅｓｈ帆ｒｓｔｈａｔｔｈｉ８ｐａｒａｌｌｅｌｃｏｄｅｂｅｈａｉｖｅｓｓａ土ｉｓｆａｃｔｏｒｙ８ｃａｌａｂｉｌｉｔｙｉｎｆｒａｓｔｒｕｃｔｕｒｅ，ａｎｄ木２００４年１２月３１日收到．１）本课题得到国家重点基础研究发展规划课题“油藏模拟与波动问题及其反问题计算（Ｇ１９９９０３２８０３）”和中国２期曹建文等：大规模油藏数值模拟软件并行计算技术及在Ｂｅｏｗｕｌｆ系统上的应用进展８７ａｎｄｃｏｍｐｕｔａｔｉｏｎａｌｓｐｅｅｄｕｐｗｉｔｈｏｕｔｓａｃｒｉｆｉｃｉｎｇｔｈｅｒｅｓｏｌｕｔｉｏｎｏｆｓｏｌｖｉｎｇｐｒｏｂｌｅｍｓ．ＮｕｍｅｒｉｃａｌｒｅｓｕｌｔｓｄｉｓｐｌａｙｔｈａｔｔｈｅｔｏｔａｌｃｏｍｐｕｔｉｎｇｃａｐａｂｉｌｉｔｙｏｆｓｏｌｖｅｒＰＲＩＳＯｉｌａＢｅｏｗｕｌｆｓｙｓｔｅｍｓｐｅｅｄｓｕｐｔｏ１６００ｃｏｍｐａｒｉｎｇｗｉｔｈｔｈｅｒｅｓｕｌｔｓｏｆ１９９８，ｔｈｅｔｏｔａｌｃｏｍｐｕｔｉｎｇｃａｐａｂｉｌｉｔｙｓｐｅｅｄｓｕｐａｔｌｅａｓｔ１０ｃｏｍｐａｒｉｎｇｗｉｔｈＤａＱｉｎｇｉｍｐｏｒｔｅｄｓｉｍ—ｕｌａｔｉｏｎｓｙｓｔｅｍｉｎ２０００．ａｎｄｔｈｅｐａｒａｌｌｅｌｓｃａｌａｂｉｌｉｔｙｍａｙｒｅａｃｈｔｏ１２８－２５６ＣＰＵｓ．ＴｈｅｒｅｌａｔｉｖｅｖｉｓｕａｌｉｚａｔｉｏｎｓｙｓｔｅｍｏｆＰＲ＿ＩＳｉｓｇｉｖｅｎｏｕｔａｔｔｈｅｅｎｄｏｆｔｈｅｐａｐｅｒ．Ｋｅｙｗｏｒｄｓ：ＨｉｇｈＰｅｒｆｏｒｍａｎｃｅＣｏｍｐｕｔｉｎｇ，ＰｅｔｒｏｌｅｕｍＲｅｓｅｒｖｏｉｒＳｉｍｕｌａｔｉｏｎ，ＳｏｌｖｅｒｏｆＮｏｎｌｉｎｅａｒＥｑｕａｔｉｏｎｓ，Ｐｒｅｃｏｎｄｉｔｉｏｎｉｎｇ，ＢｅｏｗｕｌｆＳｙｓｔｅｍ，Ｖｉｓｕａｌｉｚａｔｉｏｎ§１．引言油藏数值模拟是通过求解大型非线性偏微分方程组来模拟油田开发史，进而描述地下流体的分布，以了解和控制油气田生产动态，为控制和选择优化的开发方案提供科学依据．它属于典型的多孔介质流体流动问题，是大规模非线性问题求解算法的极好的应用平台之一．正如美国ＡＳＣＩ计划的ＳＣＩＤＡＣ（ＳｃｉｅｎｔｉｆｉｃＤｉｓｃｏｖｅｒｙＴｈｏｕｇｈＡｄｖａｎｃｅｄＣｏｍｐｕｔｉｎｇ）一ＴＯＰＳ（ＴｅｒａｓｃａｌｅＯｐｔｉｍａｌＰＤＥＳｉｍｕｌａｔｉｏｎｓ）项目指出的：多物理量耦合的非线性偏微分方程组针对多个应用领域（包括加速器设计、冶炼、燃烧、地下水流动、油藏模拟、气候等）提供了共同的数学表达形式，问题的规模目前普遍呈百万量级；随着大规模计算机的出现，对求解规模的需求增长乃至数千倍；传统的求解技术远不能作到计算复杂性的优化，从而在高性能计算的应用过程中不可避免地出现了瓶颈，限制了应用软件的可扩展性．美国ＬｏｓＡｌａｍｏｓ国家实验室、ＡｍｏｃｏＰｒｏｄｕｃｔｉｏｎ公司、ＣｒａｙＲｅｓｅａｒｃｈ公司联合研制高性能的油藏数值模拟并行软件ＦａｌｃｏｎｔｌＪ时指出，大规模线性方程组的并行求解效率成为影响油藏模拟并行软件整体性能的瓶颈．充分发挥大规模线性方程组的求解性能，可以极大地提升并行软件的求解能力．如何从大规模并行计算的角度改进线性方程组求解和预处理技术成为并行油藏数值模拟的关键．以１９９７年推出的Ｆａｌｃｏｎ为标志，美国几大石油公司先后形成了商业化的油藏数值模拟并行软件．深受用户欢迎并在市场上取得成功的软件如Ｈａｌｌｉｂｕｒｔｏｎ—Ｌａｎｄｍａｒｋ公司的ＰａｒａｌｌｅｌＶＩＰ［引，Ｓｃｈｌｕｍｂｅｒｇｅｒ—ＧｅｏＱｕｅｓｔ公司的ＥＣＬＩＰＳＥｔ刮等．９０年代，许多石油公司和研究机构也形成了自己的并行软件［４－ｓＪ．中科院软件所并行计算实验室自１９９２年以来，致力于并行算法与高性能软件的研究．在为国产高性能计算机进行正确性和性能测试过程中，与中国三大石油公司ＣＮＰＣ，ＳｉｎｏＰＥＣ，ＣＮＯＯＣ在不同层面上展开合作，以８０年代末引进的黑油模型串行软件为蓝本，进行了油藏模拟软件并行化、线性问题并行求解器、非线性问题高效求解技术、针对耦合ＰＤＥ的自适应并行预处理技术等研究与开发．１９９９年针对并行油藏模拟器软件推出了基于分布式高性能计算机的并行求解器ＰＲＩＳＶ１．１，并在一系列国产高性能计算机上成功地模拟了大庆油田提供的百万网格点实际数据吵２００１年以来，针对中国科学院陆续安装的万亿次科学计算平台，科研人员对大规模油藏模拟并行软件中的高效求解及预处理技术进行了深入研究，在此基础上，实现了基于６４—２５６处理器规模的非线性问题并行求解器，并改进了并行模拟器的实现策略，形成了新版本的高效并行求解器ＰＲＩＳＶ２．１．在利用Ｂｅｏｗｕｌｆ系统进行数数值计算与计算机应用２００６焦能方面取得了突破性的进展．文章首先介绍精细油藏模拟对时效性的需求，讨论油藏模拟系统对计算机软硬件的要求，并给出本文使用的并行计算环境和工业应用实例．随后章节给出了非线性问题求解器的算法和预处理实现技术．第四节对计算机硬件配置给并行计算性能评估带来的波动进行分析．第五节给出Ｂｅｏｗｕｌｆ机群上的数值实验结果和性能统计，用实际生产数据展示ＰＲＩＳＶ２．１对大规模油藏模拟计算的求解效果和计算效率．文章的最后介绍并行模拟器的配套可视化显示系统．§２．计算环境与应用实例精细油藏模拟对于研究剩余油分布，促进我国相当部分已进入高含水后期老油田的高产、稳产具有重要意义．它的计算规模实际需要达封百万一千万量级，历史拟合时间一般超过２０年．为了解决这样一个具有重大现实意义的问题，国际上～致认为，油藏数值模拟系统应当拥有高性能计算机和高效的数值模拟应用软件．具有工业应用价值的油藏数值模拟系统对时效性有着苛刻的要求．为了精确描述地下剩余油的分布，油藏工程师需要通过对生产过程的历史拟合，使用回归分析手段来寻找油藏特性参数（渗透率、空隙度、压缩率、边界条件等）的最佳逼近．一次成功的回归分析平均需要３０次相同规模的油藏数值模拟计算吼石油界的多名技术专家指出，一次数值模拟计算耗时三天以上是不能容忍的，夕发朝至（１６小时以内）的数值模拟基本可以满足需求，朝发夕至（８小时以内）的模拟计算才真正具有工业应用价值．明显地，我们需要从高性能计算机、高效率并行软件两个方面着手去提高油藏数值模拟计算的效率．在计算机系统组成部件出现标准化、商品化趋势的情况下，拥有高性能计算机的门槛大大降低了，这集中表现为Ｂｅｏｗｕｌｆ系统的出现．它将大量商品化硬件和开放源码软件组装成“个人超级计算机”，最大特点是廉价和ＦＬＯＰＳ密集．高性价比Ｂｅｏｗｕｌｆ系统一Ｌｉｎｕｘ机群的推出，为解决我国石油开采中的大规模数值模拟问题提供了条件．文章中使用的Ｂｅｏｗｕｌｆ系统，包括了安装在中科院数学与系统科学研究院的深腾１８００机群，安装在中科院计算机网络信息中心的深腾６８００机群，以及软件所自行组装的ＲＤＣＰＳ微机机群．深腾１８００万亿次机群系统包括２５６个计算结点和两个ｉ／ｏ结点，每个计算结点拥有两个Ｉｎｔｅｌ２ＧＨｚＸｅｏｎ处理器和１ＧＢ内存，双精度浮点运算峰值性能为４Ｇｆｌｏｐｓ，同时配备了快速以太网和Ｍｙｒｉｎｅｔ２０００，其中Ｍｙｒｉｎｅｔ２０００的单向传输带宽达２Ｇｂｐｓ，短消息延迟低于１０微秒．编译器是ＧＮＵＣ和ＩｎｔｅｌＦｏｒｔｒａｎ；并行编程环境为ＭＰＩＣＨ－ＧＭ１．２．５．１０．深腾６８００为５万亿次超级计算机系统，包括用于计算的２５６个四路结点机，１０２４个Ｉｔａｎｉｕｍ２１．３Ｇｈｚ处理芯片，内存总容量２．６ＴＢ．高速连接网络为ＱｓＮｅｔ，点对点通信带宽大于每秒３００ＭＢ，延迟时间小于７微秒．编译器是ＧＮＵＣ和ＩｎｔｅｌＦｏｒｔｒａｎ，并行编程环境为ＭＰＩＣＨ１．２．４一ＱｓＮＥＴ．ＲＤＣＰＳ机群采用ＰｅｎｔｕｍＩＩＩ５００ＭＨｚ，局部内存２５６ＭＢ，采用１００Ｍｂ／ｓ的Ｅｔｈｅｒｎｅｔ网卡和ＩＢＭ８２７５—３２４１００Ｍ交换机．编译器为ＧＮＵＣ和９７７．并行编程环境为ＭＰＩＣＨ】．２．】．２期曹建文等：大规模油藏数值模拟软件并行计算技术及在Ｂｅａｗｕｌｆ系统上的应用进展８９文章主要使用了三套实际应用数据．实例ｌ源于国外某油田【１５Ｊ，属于混相驱双孔介质模型，网格剖分为１３０ｘ１２３Ｘ１０＝１５９９００，包含３７７口井，５种岩石类型，１０个数据区块，历史拟合３１．５年．实例２为国内大庆油田某块油区的１／６，属于黑油模型，网格剖分为１９９Ｘ８７Ｘ６７＝１１５９９７１，涉及２９１口井，６种岩石类型，历史拟合３１．５年【ｌｔｉＪ．该实例的特点是，生产制度变化平稳，模拟过程的平均时间步长可以达到９０天．实例３为国内胜利油田孤东七区西区块，属于黑油模型，网格点规模为３２０ｘ６４０Ｘ２７＝５５２９６００，分布有３２６口井，生产历史１４年．该实例生产制度变化剧烈，非线性问题求解困难，对预处理的要求高，数值模拟的平均时间步长为３ｌ天．§３．非线性问题并行求解器油藏数值模拟需要求解的三方程耦合在一起的偏微分方程组，展现出椭圆、双曲、抛物型交织在一起的混合性质…．进一步的研究表明，压力物理量主要表现出椭圆型的性质，油饱和度主要呈现输运方程的一阶双曲特性，气饱和度则呈现明显的对流扩散性质．偏微分方程组离散化之后，得到一个基于有限差分的非线性方程组．在非线性问题求解中，带有阻尼和回溯功能的不精确牛顿方法仍是日常使用的主流方法．拟牛顿方法可以提高求解性能，但是必须有牛顿方法的协助，需要以双倍的程序编写复杂度来换取求解效率的提高【１１Ｊ．本文的求解器采用了不精确牛顿方法，其初值选取，在离精确解偏离较大时采用了ＢＦＧＳ的拟牛顿方法来解决．牛顿法的重点在于如何高效地求解线性代数方程组．考虑到油藏模拟偏微分方程的特性，ＬＵ分解或其它的单一预处理无法起到明显的、稳定的、正的加速作用．特殊预条件子的构造成为必然．通过对Ｋｒｙｌｏｖ子空间算法和预条件子组合进行比较Ｍ，我们选取了ＦＧＭＲＥＳ（ｍ）和ｉｎｎｅｒ—ｏｕｔｅｒ双重迭代模式作为线性方程组求解算法．数值实验表明，尽管迭代预条件ＦＧＭＲＥＳ不属于Ｋｒｙｌｏｖ子空间迭代法，从计算的稳定性、计算效率的角度看，它实际上优于子空间迭代意义上的ＧＭＲＥＳ．在非线性问题求解中，与油藏模拟领域关联最紧密的是预处理．预条件子加速效果的好坏与偏微分方程的特性密不可分．我们针对偏微分方程的耦合特点，采用了解耦合的预处理选项．针对不同物理量的行为差异，采用了约束剩余预处理选项．考虑到大规模并行处理的特点，利用区域分解的思想，将求解区域剖分为若干子区域，从而采用了区域分解类的预处理选项一加法Ｓｃｌｌｗａｒｚ．每个子区域内部充分利用局部性质，采用纯代数结构的预处理不完全ＬＵ分解，子区域之间采用粗网格校正．与油藏模拟物理背景相关的预处理和与代数结构、并行处理相关的预处理之间，采用多步法的处理手段（类似于区域分解中的乘法Ｓｃｈ嗡ｒｚ，迭代算法中的Ｇａｕ８８一Ｓｅｉｄｅｌ迭代），将它们有机联系起来，使得它们之间相互配合，优势互补，对整体的预处理系统而言，呈现出叠加效果，从而形成一个整体的预条件子“系统”【１３Ｊ．§４．硬件配置对计算性能的影响以深腾１８００为代表的Ｂｅｏｗｕｌｆ并行计算系统，一个典型的特点是，每个计算节点上配置了双路处理器．这种硬件配置模式往往会给并行程序的性能评估带来波动，从而造成测试人员感觉上的困惑．在这种计算平台上考虑并行计算性能，需要区分“单节点单进程”和数值计算与计算机应用２００６年“单节点双进程”（在每个计算节点上同时运行两个进程）两种模式．事实上，这两种并行配置模式会给并行计算绝对性能的评估带来显著的影响．我们利用ＰＲＩＳＶ２．１对实例１进行数值模拟，它需要内存量仅为５１２ＭＢ，读写文件的频率较低，单进程的读文件总量小于４０ＭＢ，写文件总量小于６０ＭＢ．明显地，它不会出现内存ｓｗａｐ，Ｉ／Ｏ操作瓶颈，Ｉ／０比重大等干扰因素．计算结果见表１．处理器节点数目１进程／处理器节点２进程／处理器节点１１７．３０１８．２７２１１．９９８。
油藏工程百科知识
油藏工程是一门以油层物理、油气层渗流力学为基础，从事油田开发设计和工程分析方法的综合性石油技术科学。
它的任务是：研究油藏(包括气藏)开发过程中油、气、水的运动规律和驱替机理，拟定相应的工程措施，以求合理地提高开采速度和采收率。
20世纪30年代以前,油田开发工作处于自发阶段,缺乏理论指导，发现油田后密集钻井，浪费很大，采收率不高。
后来随着大型高产油的发现，出现了深井压力计、高压取样器等研究油、气、水在地下状态的仪器和设备，通过对油藏岩心的研究，了解油藏和油、气、水的物理性质及其随压力、温度的变化状况和流动机理，40年代形成了油、气、水在油层中的渗流理论，出现科学开发油田的概念，逐渐应用人工补给油藏能量合理驱替油气等开发方法。
油藏工程开始成为一门独立的学科。
现代大型高速电子计算机的出现，研究油田开发的数值模拟方法的应用，以及石油开发地质和海上油、气田的勘探、开发工作的发展，进一步丰富了油藏工程的内容。
油气藏开发设计油藏工程的主要工作内容。
对于油田开发方案要分析是否采用了适合油藏特点的最有效的开采机理，最合理的井网，最有效的控制开采过程中水油比、气油比的方法；比较逐年原油采出最及所能达到的采收率和投资、油田建设工作量和所需材料，原油成本和利润。
从众多的方案中选出符合油田开发方针、能获得最高的原油采收率和最大经济效益的方案。
油藏开发动态分析油田投入生产后，地下油、气、水的分布便不断发生变化。
通过生产记录和测试资料，综合分析油井压力、产量和油藏中剩余油的分布状况等预测未来动态，提供日常生产和调整开发设计的主要依据。
具体内容有：①通过油田生产实况，不断地加深对油藏的认识，核对、补充同开发地质和油藏工程有关的各项基础资料，进一步核算地质储量；②查明分区分层油、气、水饱和度和地层压力变化，研究油、气、水在储层内部的运动状况;③分析影响采收率的各项因素,预测油藏的可采储量；④根据已有的开采历史，预测未来生产状况和开发效果。
油藏数值模拟技术现状与发展趋势
油藏数值模拟技术现状与发展趋势摘要：介绍了当前国内外油藏数值模拟的现状，简述了并行算法、网格技术、粗化技术、数值解法、动态油藏模型建立、动态跟踪模拟及三维显示等技术，指出了数值模拟的发展趋势。
关键词：并行算法；网格技术；网格粗化；分阶段模拟；动态跟踪模拟；数值解法引言近年来，随着计算机、应用数学和油藏工程学科的不断发展，油藏数值模拟方法得到不断的改进和广泛应用。
通过数值模拟可以搞清油藏中流体的流动规律、驱油机理及剩余油的空间分布；研究合理的开发方案，选择最佳的开采参数,以最少的投资,最科学的开采方式而获得最高采收率及最大经济效益[1]。
经过几十年的发展,该技术不断成熟和完善并呈现出一些新的特点。
1 国内外现状1.1 并行算法并行算法是一些可同时执行的诸进程的集合这些进程互相作用和协调动作从而达到给定问题的求解[2]。
并行算法首先需合理地划分模块,其次要保证对各模块的正确计算,再次为各模块间通讯安排合理的结构,最后保证各模块计算的综合效果并行机及并行软件的开发和应用将极大地提高运算速度,以满足网格节点不断增多的油藏数值模型。
在并行计算机上使用并行数值解法是提高求解偏微分方程的计算速度,缩短计算时间的一个重要途径[3,4]。
在共享内存的并行机上把一个按向量处理的通用油藏模拟器改写成并行处理是容易的,但硬件扩充难；分布内存并行机编程较共享式并行机困难,但硬件扩充容易,关键是搞好超大型线形代数方程组求解的并行化。
并行部分包括输入输出、节点物性、构造矩阵、节点流动及井筒等。
1.2 网格技术为了模拟各种复杂的油藏、砂体边界或断层渗透率在垂向或水平方向的各向异性,以及近井地区的高速、高压力梯度的渗流状态,近年来在国外普遍发展了各种类型的局部网格加密及灵巧的网格技术。
这种系统大体可以分为二类：一类称控制体积有限元网格(CVFE),这是将油藏按一定规则剖分为若干个三角形以后,把三角形的中心和各边的中点连接起来所形成的网格。
【数值模拟】各种方法汇总
油藏数值模拟基本过程一、数值模拟发展概况30年代人们开始研究地下流体渗流规律并将理论用于石油开发；50年代在模似计算的方法方面，取得较大进展；60年代起步，人们开始用计算机解决油田开发上的一些较为简单间题，由于当时计算机的速度只有每秒几万到几十万次，实际上只能做些简单的科学运算；70 年后主要体现于计算机的快速升级带动了油藏数模的迅猛发展，大型标量机计算速度达到100--500万次，内存也高增主约16兆字节。
在理论上黑油模型计算方法更趋成熟，D. W. Peaceman的<油藏数值模似基础>以及K. Aziz和A. Settari的<油藏模似>等主要著作都是在这个阶段出版的，但仍受到计算机速度和内存的限制，使用的方法一般仅限于IMPES及半隐式等，只能解决中小型油藏的模拟应用问题；80年代则是油藏数值模似技术飞跃发展的年代，解决不同类型油藏的数模计算方法及软件相应问世，同时超级向量机的诞生，使计算机速度达到亿次，甚至几十亿次，内存高达10—20亿字节。
90年代特别是后期，油藏模似软件各模块功能也有了惊人的发展，主要体现为向一体化方面发展；即集地震、测井、油藏工程（数模）、工艺及地面集输、经济评价等为一体的大型软件方面发展。
目前油藏数值模似软件基本上形成了一套能处理各种类型油气藏和各种不同开采方式的软件系列。
?黑油模型已被广泛用于各种常规油气藏的模拟；?裂缝模型可用来解决除砂岩以外的灰岩、花岗岩、凝灰岩和变质岩的裂缝性油气藏开发问题；?组分模型用于凝析气藏、轻质油、挥发油藏的开发设计和混相驱的研究；?热采模型用于稠(重)油油藏蒸气吞吐、蒸汽驱和就地燃烧的设计；?化学驱模型用于在注入水中添加聚合物、表面活性剂、碱等各种化学剂进行三次采油提高采收率的计算和设计。
油藏数值模拟方法的新突破随着计算机运算速度的提高,向量算法的出现和应用是软件设计上一个划时代的发展。
预处理共轭梯度法更快速、有效地解各种更为复杂和困难的大型稀疏线性方程组。
油藏数值模拟方法.pptx
油藏数值模拟功能包括两大部分：①复杂渗流力学研究，②实际油气藏开发过程整体模 拟研究，且可重复、周期短、费用低。
合并在一起,使油藏的数值模拟的网格系统反映出地质沉积特点。网格合并可以按不同井组、 区块进行合并计算, 为井组模型和分区模拟提供数据模型。模拟还可以按不均匀网格, 考虑 水平方向非均质性及储量分布程度因素等进行内插计算, 提供不均匀网格模型。
⑤动态地质建模 动态地质建模是壳牌公司的 Kortekass 概括了当前世界上关于油藏地质建模的经验, 提出的建立动态、集成化油藏模型的新概念和技术方法。其强调把动态资料以至数值模拟技 术等应用于油藏建模, 从而使所建立的地质模型更加符合油藏的实际情况, 并且要随着油田 开发中资料的增多和新资料的获得而不断更新。
1.2 油藏数值模拟方法
油藏数值模拟方法是利用计算机技术模拟地下油气藏开采、驱替的过程，是石油地质人 员科学认识、评价油藏的重要技术手段。例如，中石油公司进行的前处理的地质建模工作、 清华大学核研院研发的油藏数值模拟管理平台（PNSMP ）、大庆油田有限责任公司勘探开发 研究院研发的 VIP 和 Simbest 格式数据文件相互转换的程序等。油、气、水三相流广泛存 在于石油工业中，对于三相流的测量具有重要的意义。
油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来 模 拟实际的油田开采的一个过程。基本原理是把生产或注人动态作为确定值，通过调整模 型的 不确定因素使计算的确定值(生产动态)与实际吻合。其数学模型，是通过一组方程组 ，在一 定假设条件下，描述油藏真实的物理过程。充分考虑了油藏构造形态、断层位置、油砂体分 布、油藏孔隙度、渗透率、饱和度和流体 PVT 性质的变化等因素。这组流动方程组由运动 方程、状态方程和连续方程所组成。油藏数值模拟是以应用数学模型为基础的用来再现油田 实际生产动态的过程。具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学， 借助大型计算机为介质条件建立三维底层模型参数场中，对数学方程求解重现油分析
油藏数值模拟
油藏数值模拟油藏数值模拟是随着电子计算机的出现和发展而成长的一门新学科，在国内外都取得了迅速的发展和广泛的应用。
1953年美国G..H.BUCE等人发表了《孔隙介质不稳定气体渗流的计算》后，为用数值方法计算油气藏渗流问题开辟了道路。
三十多年来，由于大型快速电子计算机的迅速发展，大大地促进了数值模拟方法的广泛应用。
20世纪60年代初期研究了多维多相的黑油模型；20世纪70年代初期研究了组分模型、混相模型和热力采油模型；20世纪70年代末期研究各种化学驱油模型。
目录编辑本段油藏数值模拟方法是迄今为止定量地描述在非均质地层中多相流体流动规律的惟一方法。
例如许多常规方法要假定油层为圆形的均匀介质，如油藏几何形状稍复杂一些，且为非均质介质，则求解非常困难，甚至无法求解。
而对油气藏数值模拟而言，计算形态复杂的非均质油藏和计算简单形态的均质油藏工作量几乎是一样的。
因此油藏数值模拟可解决其它方法不能解决的问题。
对于其它方法能解决的问题，用数值模拟方法可以更快、更省、更方便、更可靠地解决，并增加其它分析方法的可信度。
一个油气藏，在现实中只能开发一次。
但应用油藏数值模拟，可以很容易地重复计算不同开发方式的开发过程，因此人们可以从中选出最好的开发方法。
因此，对油藏工程师而言，数值模拟给动态分析提供了一种快速、精确的综合性方法；对管理者而言，数值模拟提供了不同开采计划的比较结果；对尚无经验的工程师而言，数值模拟则是有效的培训工具。
编辑本段数模研究数值模拟研究的主要工作程序对一个油气藏进行综合的数模研究，往往需要花较大的精力和较长时间（有时会达一年甚至更长的时间），同时还对计算机硬件和技术人员有很高的要求，然而尽管在不同的项目中，面对的问题会千差万别，但大多数油藏数值模拟的基本研究过程是一样的。
为了使读者一开始就对数模研究工作有一个明确的整体概念，下面简要地介绍一下油藏数值模拟的主要工作程序。
问题的定义：开展油藏数模工作的第一步，是确定研究的目标和范围。
油藏数值模拟方法.pdf
第一章油藏数值模拟方法分析1.1油藏数值模拟1.1.1油藏数值模拟简述油藏数值模拟是根据油气藏地质及开发实际情况，通过建立描述油气藏中流体渗流规律的数学模型，并利用计算机求得数值解来研究其运动变化规律。
其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。
其基础理论是基于达西渗流定律。
油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来模拟实际的油田开采的一个过程。
基本原理是把生产或注人动态作为确定值，通过调整模型的不确定因素使计算的确定值(生产动态)与实际吻合。
其数学模型，是通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。
充分考虑了油藏构造形态、断层位置、油砂体分布、油藏孔隙度、渗透率、饱和度和流体PVT性质的变化等因素。
这组流动方程组由运动方程、状态方程和连续方程所组成。
油藏数值模拟是以应用数学模型为基础的用来再现油田实际生产动态的过程。
具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学，借助大型计算机为介质条件建立三维底层模型参数场中，对数学方程求解重现油田生产历史，解决实际问题。
油藏数值模拟技术从50 年代的提出到90 年代间历经40 年的发展，日益成熟。
现在进入另外一个发展周期。
近十年油藏数值模拟为油田开发研究和解决实际决策问题提供强有力的支持。
在油田开发好坏的衡量、投资预测及油田开发方案的优选、评价采收指标等应用非常广泛。
油藏数值模拟功能包括两大部分：①复杂渗流力学研究，②实际油气藏开发过程整体模拟研究，且可重复、周期短、费用低。
图1 油藏数值模拟流程图1.1.2油藏数值模拟的类型油藏数值模拟类型的划分方法有多种，划分时最常用的标准是油藏类型、需要模拟的油藏流体类型和目标油藏中发生的开采过程，也可以根据油气藏特性及开发时需要处理的各种各样的复杂问题而设定，油气藏特性和油气性质不同，选择的模型也不同，还可以根据油藏数值模拟模型所使用的坐标系、空间维数和相态数来划分。
油藏数值模拟技术
油藏数值模拟技术2006年11月目录一、关于“油藏数值模拟技术”（一）基本概念及作用（二）数据准备（三）模型初始化（四）生产史拟合（五）生产动态预测二、油藏数值模拟的主流软件系统简介三、油藏数值模拟技术的进展及发展方向（一）进展（二）发展方向四、使用ECLIPSE软件进行油藏数值模拟的过程简介一、关于“油藏数值模拟技术”油藏数值模拟技术是一门将油田开发重大决策纳入严格科学轨道的关键技术。
从油田投产开始，无论是单井动态，还是整个油田动态，都要进行监测与控制。
油藏数值模拟是油田开发最优决策的有效工具。
油藏数值模拟技术从20世纪50年代开始研究至今，已发展成为一项较为成熟的技术，在油气藏特征研究、油气田开发方案的编制和确定、油气田开采中生产措施的调整和优化以及提高油气藏采收率方面，已逐渐成为一种不可欠缺的主要研究手段。
油藏数值模拟技术经过几十年的研究有了大的改进，越来越接近油气田开发和生产的实际情况。
油藏数值模拟技术随着在油气田开发和生产中的不断应用，并根据油藏工程研究和油藏工程师的需求，不断向高层次和多学科结合发展，将得到不断的发展和完善。
（一）基本概念及作用（1）基本概念油藏数值模拟：从地下流体渗流过程中的本质特征出发，建立描述渗流过程基本物理现象、并能描述油藏边界条件和原始状况的数学模型，借助计算机计算求解渗流数学模型，结合油藏地质学、油藏工程学重现油田开发的实际过程，用来解决实际问题。
油藏数学模型的分类，一般有四种方法：1)按流体中相的数目，划分为：单相流模型、两相流模型、三相流模型。
2)按空间维数，划分为：零维模型、一维模型、二维模型、三维模型。
3)按油藏特性类型，划分为：气藏模型、黑油模型、组分模型。
气藏模型按其组分的贫富，可以用黑油数值模型模拟，也可以用组分类型的数值模拟模型模拟。
所以，气藏模型也可以划进黑油或组分模型。
故数学模型一般分为黑油型和组分型两类模型。
4)按油藏结构特点、开采过程特征，分类为：裂缝模型、热采模型、化学驱模型、混相驱模型、聚合物驱模型等。
5 国内外油藏数值模拟软件
中国石油大学（北京） 檀朝东 tantcd@
目 录
1.油藏数值模拟技术整体概要介绍 2.国外油藏数值模拟软件介绍 3.国内油藏数值模拟软件介绍 4.各种主流油藏数值模拟软件的比较
1.油藏数值模拟技术整体概要介绍
1.1油藏数值模拟介绍、意义
模拟：基于现实的各种假设条件，通过一个模型模拟真实系统 模拟：
拉斯收购. 1996年兰德马克公司收购了西方阿特拉斯数模业务.
11
1986年INTERA收购ECL公司，1995年斯仑贝谢GeoQuest收购了
INTERA公司软件部.
12
公司。其开发的数值模拟软件为VIP。
兰德马克于1996年7月被哈里伯顿公司宣布收购，成为其下子
1.油藏数值模拟技术整体概要介绍
2.国外油藏数值模拟软件介绍
2.1 ECLIPSE油藏数值模拟软件
2.1.3
ECLIPSE出图效果
2.国外油藏数值模拟软件介绍
2.1 ECLIPSE油藏数值模拟软件
2.1.4
ECLIPSE主要特点
1.能够模拟各种类型的油气藏.包括砂岩油气藏、裂缝性双孔双渗油气藏、凝 析油气藏和低渗透油气藏。Eclipse100提供一个GI拟组份模型选项，用于模拟 凝析油气藏。 2.能够模拟多种开采方式(包括一、二、三次采油)。诸如衰竭开采、注水开采 、注气开采和循环注气开采。三次采油中有多种选项,其中包括聚合物驱、混合 相驱(包括气水交替混驱)、表面活性剂驱和泡沫驱等。 3.Eclipse串行计算可在单处理器机器上运行： Sun,SGI,IBM,PCs(Win98/NT/2000,Linux)。 4.Eclipse并行计算可在多处理器机器上运行：Sun,SGI,IBM,SP,PC。 5.具有全隐式解法，AIM和IMPES解法；
5 国内外油藏数值模拟软件
13
SimTech成立于1982年,1996年5月16日被Intera收购。Intera为
Duke工程服务公司的子公司,主要做环境及水资源评价。
14
数值模拟研究公司(RSRC)成立于1982年，1996年被Smedvig
收购。 1999年Smedvig技术公司与多流体ASA合并组成新的Roxar公司. 模拟软件为MORE, MORE为Modular Oil Reservoir Evaluation的缩写 . Hot (Heinemann Oil Technology and Engineering)公司1986年成立 于奥地利Leoben，其模拟软件为SURE,其主要特色是采用PEBI网 SURE 格.
VIP系列
GrandTM
3DSL
1.油藏数值模拟技术整体概要介绍
1.3油藏数值模拟的分类、现状和未来发展
1.3.3高级油藏数值模拟技术的发展方向：
多功能集成
在一个模拟器中整合黑油,组分， 热采 模型；整合全隐式，压力隐 式和自适应隐式等不同格式；整合 结构化和非结构化网格统；整合传 统井模型和智能井模型。
油藏岩石流体模型：
① 初始 条件
通过测试手段了解油 水气界面信息； 通过流体采样，实验 室分析了解组分信息等； 通过岩心取样，实验 室分析了解岩石信息等。
1.油藏数值模拟技术整体概要介绍
1.2油藏数值模拟的初始条件,边界条件和运动规律
1.2.1运用模型模拟现实的三大要素 ：
地质模型边界： 井：
封闭边界：
发展方向
1.油藏数值模拟技术整体概要介绍
1.3油藏数值模拟的分类、现状和未来发展
1.3.4参考文献：
1 2 3
数值模拟技术汇总
孔隙介质中通常有三种组分：油气水； 通常假定油水为非混相状态，油水之间无质量传递； 通常假定气水之间有质量传递，以随压力而变化的 溶解度的形式； 通常假定气肯定溶解在油相中； 油藏温度恒定，流体一致处在热动力平衡之中； 流体的特性只是压力的函数，不考虑组分的影响。
二、数值模拟技术的方法原理
一、数值模拟技术简介
5 软件应用分类
黑油模型模拟软件：用于水驱油藏数值模拟研究。
ECLIPSE 软件 VIP软件 CMG软件IMAX模块 SURE软件
（斯伦贝谢公司） （Landmark公司） （加拿大CMG公司） （GCC 公司）
组份模型模拟软件：用于气藏、凝析气藏数值模拟研究。
CMG-Gem 模型 Eclipse 300模块
（加拿大CMG公司） （斯伦贝谢公司）
一、数值模拟技术简介
5 软件应用分类
稠油热采模拟软件：用于稠油吞吐，蒸汽驱、火烧等热采模拟研究。
CMG软件stars模块
（加拿大CMG公司）
聚合物及化学驱模拟软件：用于聚合物驱油及复合驱油藏模拟研究。
ECLIPSE软件聚合物驱模块 VIP软件聚合物驱模块 UTCHEM化学驱软件 FACS稀体系化学驱软件
物质守恒方程
一、数值模拟技术简介
1 油藏数值模拟的基本概念
• 数学描述 ——用数学方程描述实际油藏及油藏 开采动态的过程。
━━ 物质平衡定律
原来 + 进来 ─ 出去
━━ 达西定律 ━━ 状态方程
= 目前
分析基础 时刻牢记
一、数值模拟技术简介
1 油藏数值模拟的基本概念
• 网格模型
定
一、数值模拟技术简介
2 技术发展历程
化学驱数值模拟
式中K31由实验确定
a1C r1 r 1 a1C
0 1
水相中CO2引起的碱耗r4
K 42 CO 2 OH HCO
C — OH 的浓度； r — 该现象引起的最大碱耗 ，由实验资料确定； a1 — 系数，由实验资料确定 。
0 1
可根据化学反应折算出CO2引起的碱耗 Ca++、Mg++离子引起的碱耗r5
界面张力计算 为常数 毛管数的变化 不考虑 粘度 离子交换 生成表活剂 求解方法 水相粘度为常数，油与气的 粘度随压力而改变 不考虑 不考虑 全隐
- 10 -
主要物化现象及物化参数描述 界面张力随化学剂浓度的变化：
为了真正体现多种化学剂的复合协同效应， 采用实测的界面张力等值图进行描述，对于给 定的原油和配置水矿化度，界面张力主要随碱 、表面活性剂、聚合物浓度变化：
f（C po , C p )
- 13 -
物化现象及参数描述
残余饱和度和相对渗透率 残余饱和度
各相残余饱和度与毛管数有 关，毛管数定义:
各相相对渗透率
K rj K S
0 rj ej nj
Nc
j jV j
S nj
1 j S rj
S j S rj
S rj S rj ( N c )
[ yij j
j
K krj
j
(Pj j g) S j j Dij yij ]
j 1
(i Zi j S j qi )
j
t
Bi
-8-
INPUT:参数输入 INIT: 初始化 PARAMETERS: 物化参数计算 RESOLVEP: 求解压力化 GAUSS: 高斯解法 MATRICE 系数矩阵 LSOR: 线松弛 PCG: 预处理共轭梯度
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
万方数据
软件天地中文核心期刊‘徽计算机信息＇（测控自动化｝２００７年第２３卷第１０－１期
对整个油｛葭进行网格划分后，将所有网格块的质量平衡方点，这样保证了一定的空何局部性和静态的负载平衡．每个节程鬟妻至：耋毒黑髦油。
！：兰黧；！置：善篇薯鼍亍点拥有Ⅳ‰，＝笋盟＋‰甜同格块，其中Ｎａ—ＬＩ，。
，是
一个不定常非线性偏微分组．其中未知量数目为Ｎ・（Ｎ—Ｎ一２），…’。
”出”Ⅳ７Ⅳ。
ｍⅢ
…一…’一。
’，“。
Ｎ为网格敦．其初始条件是零时刻的物理场分布，其边界条件是分配刊子区域的网格数．Ｎｎｏ【，Ｙ，ｚ）是ｘ、Ｙ、Ｚ方向上的总网格在边界上投有流体的流人和流出，利用上述给定的初始条件和数，Ｎ倒＾，ｙ，ｚ）是问题分解在ｘ、Ｙ、Ｚ方向的处理器数，Ｎ。
柚是边界条件就可以对整个方程组进行求解。
子空间在给定方向上的邻近子区域数。
一个１５ｘ２０ｘ６的油藏化学复合驱果用的是Ｉｍｐｅｓ法。
其基本思路是在一个时间划分，假定可用的处理嚣数为１２，则处理器的排列方式可能是步长内只把压力作为变量，而在方程中出现的饱和度，组分浓３ｘ４×１．那么各于区域数组的维数是７ｘ７ｘ６。
图２显示了网格度等都取前一时刻的值，即只对压力变量形成系数矩阵。
求解的划分和各子区域阿格块数的兰系。
ＵＴＣＨＥＭ的并行实现应
出压力后再利用各种化学关系显式地求解其它变量。
对于压力方程，空间差分选用一鼢向前差分和二阶中心差分，采用雅克比共轭梯度法进行求解。
压力方程求解之后．得到了全场各十网格点的压力变量，然后显式求解各点舶饱和度和备组分浓度。
图１是化学复合驱程序运行流程图。
：；：；ｈ叫＊一ｍ＊卜≈黼Ｈ＊ｍ－一ｅＨｊ：嚣：４Ｈｍ—ｍｔ！竺！ｒ弋＝旦
．．．．．．．．．．．Ｌ
Ｉｈｍ∞：ｔ＊＊嘛：：。
雾、————／ｌ％自目ｔ≮：ｉ卜｝．Ｊ—Ｔ———————７—赢亍、ｒ一‘—－————一氟；孙圈１化学复合驱程序运行过程
图２油藏网格划分简单示例
３并行模型及关键技术
化学复合驱并行模型采用单程序多数据（ＳＰＭＤ）程序模型．利用区域分解方法使得多个计苒节点同时求解一个单一的模拟问题．各计算节点通过消息传递对必要的数据进行通信。
３．Ｉ初始化
执行开始时．主节点（ＮＯＤＥ０）将整个模拟空间的网格块在各处理节点之间划分．每十节点接收必要的同格数据．来保证节点内部的计算。
该划分是有主节点控制的，这个阶段它的主要职责是保证整个网络和计算的负载均衡。
３２网格划分
并行宴现的主要目的是降低执行时间和解决大规模计算嚣要的内存不足问题。
对于第二个问题，并行程序改变了原始串行程序中的数组大小，只分配子区域计算所需的内存。
即子区域的空间加一个厨绕在其周围的信封空间佃２中的ｌ！！；影网格）．信封空间用来接收来之相邻区域的数据。
网格块划分机制将连续编号的舟格映分配给虚拟计算节用区域划分技术，模拟期间各处理节点和其相邻节点需要一定的共享数据（图中阴髟部分）。
不同子区域上的模榴可以作为一个独立的问题并行的执行，也就是在不同的数据集上同步执行相同的运算（ＳＰＭＤ）。
因为每一个子区域的计苒需要槲邻区域上的一些数据，相邻节点之间在指定的同步点上进行显式的消息传输。
３．３Ｉ／ｏ控制
在这个实现方案中，ＹＯ是通过一个指定节点（ＮＯＤＥＯ）串行控制的。
程序初始阶段。
Ｉ／０节点负责读人输人数据并将数据划分到各子区域。
在程序输出阶段。
Ｉ／Ｏ节点综合所有节点发送给它的部分运算结果，然后写入到磁盘文件。
这样的控制方式带来两个缺点：一是因为指定了Ｉ／Ｏ控制节点．在它读＾输人数据时候其它节点都在等待：二是在输出阶段、其它节点都向Ｉ／０节点发送数据，１１０节点成为消息传递的瓶颈。
３＿４基于区域分解法的线性解法器
串行ＵＴＣＨＥＭ中雅克比共轭梯度法线性解法器的并行救宰不是很好。
本文采用一种基于加法Ｓｃｈｗａｒｚ两条件子的共轭梯度区域分解算法。
对于线性系统Ａｐ＝ｂ，Ａ是压力矩阵，Ｐ是未知压力向量，ｂ是边界向量值。
在每一次求解迭代ｎ中．预条件残余向量足通过求解Ｍｚ“＝，得到，其中，是当前残余量．也就是ｒｋｂ—Ａｐ。
，一是预条件残余量，Ｍ是预条件矩阵。
在加法Ｓｃｈｗａｒｚ方法中，是不需要Ｍ矩阵的．因为通过下进步骤可以计算出矩阵Ｍ“和＾
ｆ１）将整个问题的残余和系数矩阵分配到扩展子区域上（包括藿叠边界同格块）。
在每个扩展子区域上用统一的狄利克雷边界条件求解Ａ一＝一。
在扩展子区域边界和全局区域边界一致的地方，使用原始边界条件。
｛２）计算完所有的子区域后，将局部结果相加构成全局顶条件子残余。
４测试结果及结论
试验环境为基于ＳＭＰ的千亿次计算机机群系统．其中１６个计算结点各拥有２个ＩｎｔｅｌＸｅｏｎ２８ＧＨＺＣＰＵ、２Ｇ内存．节点之间由干兆以太网同连接，该集群的计算理论峰值为１７９２亿次，ＩＮＰＡＣＫ实测峰值为９９５亿次。
经过多组测试．改造后的并行程序与原串行程序的运算结果完全一致。
被并行模块（压力方程组求解）在１６个节点取得６．４的加速比。
图３中的曲线可以看出随着网格规模的增加，所取得的加速比和效率有所提高，说蹦数据规模越大越能体现并行效果。
测试结果证明区域分解算磕应用与大规模偏微分方程组求解能够取得很好的加速效果。
但实际的应用问题往往并不是耶么简单．需要很好的策划各计算节点的负载平衡问题和节点间数据通信的代价。
尤其是节点之间的通信问题。
一２５０—３６０无１年邮局订阅号：８２．９４６万方数据
万方数据。