潘举玲1)
黄尚军1)
祝 杨2)
宋志强2
)
1)胜利油田有限公司物探研究院;2)胜利油田有限公司孤岛采油厂
摘要：介绍了当前国内外油藏数值模拟的现状，简述了并行算法、网格技术、粗化技术、数值解法、动态油藏模型建立、动态跟踪模拟及三维显示等技术，指出了数值模拟的发展趋势。关键词：并行算法；网格技术；网格粗化；分阶段模拟；动态跟踪模拟；数值解法中图分类号:TE319 文献标识码:B
文章编号:1009-9603(2002)04-0069-03
越战忠魂引言
近年来，随着计算机、应用数学和油藏工程学科的不断发展，油藏数值模拟方法得到不断的改进和广泛应用。通过数值模拟可以搞清油藏中流体的流动规律、驱油机理及剩余油的空间分布；研究合理的开发方案，选择最佳的开采参数，以最少的投资，最科学的开采方式而获得最高采收率及最大经济效益[1]。经过几十年的发展，该技术不断成熟和完善并呈现出一些新的特点。
1 国内外现状
1.1 并行算法
本草纲目拾遗并行算法是一些可同时执行的诸进程的集合，这些进程互相作用和协调动作从而达到给定问题的求解[2]。并行算法首先需合理地划分模块，其次要保证对各模块的正确计算，再次为各模块间通讯安排合理的结构，最后保证各模块计算的综合效果。并行机及并行软件的开发和应用将极大地提高运算速度，以满足网格节点不断增多的油藏数值模型。在并行计算机上使用并行数值解法是提高求解偏微分方程的计算速度，缩短计算时间的一个重要途径[3,4]。在共享内存的并行机上把一个按向量处理的通用油藏模拟器改写成并行处理是容易的，但硬件扩充难；分布内存并行机编程较共享式并行机困难，但硬件扩充容易，关键是搞好超大型线形代数
方程组求解的并行化。并行部分包括输入输出、节
点物性、构造矩阵、节点流动及井筒等。1.2 网格技术
为了模拟各种复杂的油藏、砂体边界或断层，渗透率在垂向或水平方向的各向异性，以及近井地区的高速、高压力梯度的渗流状态，近年来在国外普遍发展了各种类型的局部网格加密及灵巧的网格技术。这种系统大体可以分为二类：一类称控制体积有限元网格(CVFE )，这是将油藏按一定规则剖分为若干个三角形以后，把三角形的中心和各边的中点连接起来所形成的网格。另一类则称垂直等分线排比网格(PEB I )，其剖分方法是将油藏分成若干三角形后，使三角形各边的垂直等分线相交而形成网格。
这些方法在处理复杂几何形状油藏及进行局部网格加密时简单而一致。在多相流情况下，参照某一给定的几何准则时该方法是单调的，这保证了其稳定性和收敛性。这两种方法都能以直观的控制体积的概念出发并且采用一致的上游权而推导得出。这些方法对网格的方向不敏感，在某些情况下比九点差分格式的效果好。1.3 计算机辅助历史拟合技术斯伦贝谢公司的Eclipse 数模软件最新推出计算机辅助历史拟合模块(Simopt )。运用均方差、海赛(Hessian )矩阵、相关性矩阵、协方差矩阵对结果进行分析以确定敏感参数；引入梯度带分析技术对地质模型进行优化；在进行常规历史拟合后，应用置信度限制(规定需优化的参数及参数的可调范围),
收稿日期2002-02-19；改回日期2002-04-22。
作者简介：潘举玲，女，工程师,1991年毕业于石油大学油藏工程专业，现从事油藏数值模拟工作。:(0546)8778817，通讯地址:(257022)山东省东营市泰山路344号物探研究院勘探开发计算机应用研究所。
油 气 地 质 与 采 收 率
2002年8月 PETROL EUM GEOLO GY AND RECOV ER Y EFFICIENCY 第9卷 第4期
通过线性预测分析，实现计算机辅助调整参数，减少模拟次数。
1.4 网格粗化技术
对于一些油藏参数(如孔隙度、深度、饱和度等)，采用体积加权平均法；对于与流体有关的参数(如渗透率等)就不能用简单的加权平均计算得出，而要基于流动计算再进行粗化。流动算法相对精确，首先解出沿压力降方向的总流量，然后再解相同的流动方程，从而解出等效渗透率。
在垂向分层合并计算中，把相同性质的油砂体(按相同的物性、储量类型)的网格单元合并在一起，使油藏的数值模拟的网格系统反映出地质沉积特点[5]。网格合并可以按不同井组、区块进行合并计算，为井组模型和分区模拟提供数据模型。模拟还可以按不均匀网格，考虑水平方向非均质性及储量分布程度因素等进行内插计算，提供不均匀网格模型。
1.5 动态地质建模
动态地质建模是壳牌公司的K ortekass概括了当前世界上关于油藏地质建模的经验，提出的建立动态、集成化油藏模型的新概念和技术方法。其强调把动态资料以至数值模拟技术等应用于油藏建模，从而使所建立的地质模型更加符合油藏的实际情况，并且要随着油田开发中资料的增多和新资料的获得而不断更新。
这种新方法包括一系列获得和运用各种所需资料的技术和方法，包括地质、地质统计、地震、测井、岩心和流体分析、试井、驱替特征以及网格的细分和粗化，拟函数的应用等，但关键是使所建立的地质模型更加符合油藏的实际情况，而且还可以加快建模的过程。
1.6 分阶段模拟
对开发历史较长、地下储层物性和原油物性发生较大变化的油藏，把随开发时间变化的地质静态模型划分为多个不同开发阶段的地质模型。
常规的油藏数值模拟是从一个油藏(区块)投入开发时开始模拟，一直拟合到目前状况，再进行方案预测优选。我国许多老油田已进入了高含水或特高含水期。由于开发历史长、综合调整、措施次数多，地下岩石和流体的物性发生了较大的变化，这给常规模拟工作带来了极大的困难。一方面是历史拟合
计算一次所需要的机时非常多，另一方面是常规模拟无法考虑流体和岩石随时间的变化。因此，模拟结果的可信度会大大降低。分阶段模拟就是一种解决上述问题的行之有效的方法。分阶段模拟可将一个长期开发的油藏，按照一定原则划分成几个模拟阶段。
1.7 动态跟踪模拟
油田开发是一个长期过程，储层物性和原油物性随开发期的不同以及油水井措施发生变化。根据开发期及措施类型制定数值模拟的时间步长，在油水井见效初期采用较小的时间步长，进入见效稳定期后以较大的时间步长，将分析周期由常规的以年计算提高到以月或天计算。对方案实施后的效果、生产状况等再进行跟踪模拟，并提出新的方案。如此反复研究，使人们对油藏的构造、物性、油水状况及生产动态的认识更趋合理。
1.8 数值解法
古老的差分法继续得到创造性的发展和应用，也得到深入的分析和研究。差分法的另一重大发展就是全隐式、自适应隐式方法[6]。由于全隐式方法对所有方程系数进行隐式处理，所以与IMPES、半隐式和SEQ相比，稳定性好，隐式程度高，适应范围更宽。它能解决油、气、水三相渗流、注气、水气锥进、高速气渗等强非线性渗流问题。但是全隐式模型的求解需要采用牛顿迭代法，因此，其工作量和存储量比IMPES、半隐式、SEQ大。所以,20世纪80年代初期美国的托马斯等人为了解决方程隐式程
度高低和计算量大小之间的矛盾，提出了自适应隐式方法，其特点是对不同的网格节点和不同的时间步采用不同的隐式程度来处理，以便在具有同样稳定性的前提下减小计算量、加快计算速度。1985年Berteger W I等人又提出了一种近似的自适应隐式方法，其稳定性进一步增强。
差分法的又一个重大发展就是多重网格法和预处理共轭梯度法。两者的共同点都在于加速数值解的收敛性。多重网格法实质上是外推与内插技术的创造性应用。近年来，共轭梯度法与各种不完全分解预处理相结合，并采用D4网格排序及Orthomin 加速技术，从而使预处理共轭梯度法成为20世纪80年代油藏数值模拟中最引人注目的方法之一。预处理共轭梯度法克服了油藏数值模拟中由于局部网格加密、混合网格嵌套、隐式井底压力处理、大断层大裂缝处理、死节点处理以及自适应隐式方法等产生的不规则系数矩阵的求解问题。由于Or2 thomin方法收敛性好，且运算过程简单、不需要迭代参数、不需要估计矩阵的特征值，是一个稳定、有效的迭代加速方法。因此，目前该法在油藏数值模
・
7
・油 气 地 质 与 采 收 率 2002年8月
拟中获得了相当成功的运用。
1.9 三维动态显示
数值模拟结果的可视化程度高，人机交互性能强。可三维动画显示油藏中流体的流动规律，再现油藏的开发历史及剩余油的空间分布，并可任意旋转、平移、缩放、光照，多重照相、透视和透明处理；灵活的剖面切割功能，任意参数的区间显示，用户可按自己的爱好定义或修改颜和注释。
2 发展趋势
工作站由于价格昂贵以至很难普及。而今PC 机越来越强大，速度日益加快，内存越来越大且易扩展，价格便宜，应用广泛，因而将数模软件开发在PC 机上的远景极为广阔。但目前的困难是PC并行机尚没有统一的标准且易出现网络问题。
并行算法的进一步完善、各种灵巧网格的随意化、自适应化以及历史拟合的自动化将大幅度提高数值模拟的精度，减少模拟计算的工作量。
在大规模粗化前利用流管法，预先筛选地质模型。根据物质平衡方程和达西定律，在单时间步内采用隐压显饱解法计算流线上流体的流动：全隐式有限差分压力，沿流线解饱和度方程。在考虑重力影响时，先沿流线再沿重力线求解压力和饱和度。
勘探开发一体化越来越受到关注。油气勘探开发的各个环节中计算机应用技术的突飞猛进，带动了勘
探开发工作向着多学科综合集成方向发展，改变了以往以单项专业为主的顺序化工作方式①。各学科间的信息交换、数据共享、成果可视和知识继承遵照统一的标准进行综合集成：通过岩心声波测量建立油藏参数与地震响应间的关系；利用油藏数值模拟技术预测压力和含水饱和度的变化、生成不同观测数据的合成地震记录；开发地震包括了四维地震和地震反演技术，通过井的约束以及和油藏数模技术的结合，将进一步减少多解性的困扰并大幅度提高分辨率，从而扩展开发地震的实际应用范围，特别对比较复杂的砂泥岩薄互层也能够加以应用。这些技术的发展，将大大加深对井间储层非均质特征及油藏剩余油分布状况的认识。但当前面临的困难是：种类各异的数据格式、各式各样的数据库系统、互不相通的应用软件以及不同的用户界面等严重妨碍了多学科的协同工作、数据资源及研究成果的共享。
借助多井试井数据、数值模拟和优化技术，通过对试井实测压力数据的拟合，反演油藏中某些非均质参数分布，能有效地解释测试井的径向平面渗透率变化情况。
借助于多井试井数据，结合油藏的地质信息，将突破目前建立在单层、单相、均质储层等基本假设条件下的局限，为非均质、多层、多相等复杂问题提供满意的解释。
3 结论
纵观国内外油藏数值模拟技术的研究与发展，随着计算机、应用数学和油藏工程学科的不断发展，国
外的软件技术发展非常迅速。不仅模型功能齐全，且已一体化。角点坐标描述、断层处理技术、斜井及水平井模拟、局部网格加密、并行算法等功能非常完善；并形成了商业性软件包。应用广泛的软件有斯伦贝谢公司的Eclipse系列、兰德马克公司的Desktop-Vip油藏模拟系统、SSI公司的Work2 bench油藏工作平台等。国内数模软件水平较落后，与国际商业数模软件差距太大，要形成工业化并广泛推广还需在投入力度、人才培养、软件开发、知识结构等方面做出努力。
参考文献:
1 陈月明.油藏数值模拟基础.东营：石油大学出版社,1988
2 陈国良.并行计算.北京：高等教育出版社,1999
3 舒继武，归丽忠，周维四等.区域分解法解黑油数值模拟问题的并行计算.南京大学学报,1999,1(35):51～57
4 杨耀忠，韩子辰，周维四等.多层二维二相油藏数值模拟并行技术研究.油气地质与采收率,2001,8(6)52～54
5 李福垲，由军，陈立生.地质模型网格粗化技术.石油勘探与开发,1996,23(3):76～79
6 张烈辉.裂缝性油藏水平井数值模拟的进展和展望.西南石油学院学报.1997,19(4):49～52
编辑 张朝启
・
1
吉尔吉斯族7
・
第9卷 第4期 潘举玲等：油藏数值模拟技术现状与发展趋势
①韩大匡.常规油田开发技术综述.第十五届世界石油大会论
文资料.1997
XP3ν38well as example using the techniques of pyrolytic geochemical analysis and thermal evapora2 tion hydrocarbon chromatography analysis.Oilνbearing grade of the reservoir is determined and reser2 voir nature is evaluated.The methods of geochemical log and core extraction are identical on the evaluation of the distribution of the remaining oil through the comparison between the values of the oil saturation obtained by the above two methods.
K ey w ords:geochemistry,pyrolysis analysis,ther2 mal evaporation hydrocarbon chromatography
Ma Jianhua,Liang K aifang,G ong Junfeng et al. Laboratory study and application of plugging agent of Amorphophallus rivieri powder.PGRE,2002,9 (4):64～65
Viscosity,swelling capacity,gel strength,heat resis2 tance,gel performance,scouring resistance capacity and breakthrough pressure of the plugging agent of Amorphophallus rivieri powder are measured in labo2 ratory experiments,and then the performance of the plugging agent is evaluated.Five water injection wells are operated at Shengtuo oilfield in J uly2001, response has appeared in14oil wells among corre2 sponding19oil wells with2816.5t of cumulative in2 cremental oil production.
K ey w ords:laboratory experiment,plugging agent of Amorphophallus rivieri powder,profile control by water plugging,field application
Yu Yuanzhou,Yang G uangrong,Tian Jinhai et al. Improvement and application of material balance e2 quation for condensate gas reservoir.PGRE,2002, 9(4):66～68
拉赫玛尼诺夫第二钢琴协奏曲The most effective production scheme of cyclic gas in2 jection for keeping formation pressure at pre
sent is adopted generally to improve oilνgas recovery in the condensate gas reservoir.Large deviation will appear if material balance equation of the gas reservoir for calculating OIP is used directly in this production scheme.So the equation is improved,then the mate2 rial balance equation in the gas injection development for the condensate gas reservoir is obtained.
K ey w ords:condensate gas reservoir,material bal2 ance equation,deviation factor,cyclic gas injection
P an Juling,H uang Shangjun,Zhu Yang et al.State of the art and developing tendency on reservoir nu2 merical simulation.PGRE,2002,9(4):69～71 The state of the art on the reservoir numerical simula2 tion at home and abroad is introduced.A brief de2 scription is on parallel algorithm,grid technique, rough grid,numerical solution,dynamic reservoir model building,dynamic tracing simulation,3D dis2 play and some other techniques.The developing ten2 dency of numerical simulation is pointed out.K ey w ords:parallel algorithm,grid technique,rough grid,simulation stage by stage,dynamic tracing sim2 ulation,numerical solution
Li Yingchun,Q iu G uoqing,Yuan Mingqi et al. EOR experiment of inνsitu combustion in south area of Lean oilf ield.PGRE,2002,9(4):72～74 Feasibility of development by inνsitu combustion after steam s
timulation in south area of Lean oilfield is discussed comprehensively,and reservoir engineering design of inνsitu combustion is carried out through the study on thermal recovery screening,physical simulation and numerical simulation.Initial response appeared in the experiments of the inνsitu combus2 tion in south area,and it provides reference and guid2 ing for the development of the same kind of oil reser2 voirs in Shengli petroliferous province.
K ey w ords:inνsitu combustion,heavy oil reser2 voir,later period of steam stimulation,glutenite, Lean oilfield
八点听书Zhou Jiaju,G uo Zhihua,Liu Zhiping et al.Mea2 sures for tapping potential of remaining oil in extra high w ater cut period of H e31fault block.PGRE, 2002,9(4):75～76
In order to enhancing oil recovery of He31fault block further,rules of distribution and measures for tapping potential of the remaining oil in extra high water cut period are studied started with fine reservoir descrip2 tion.Many techniques are applied to adjustment and tapping potential including directed drilling,sidetrack drilling technique and hydrodynamic methods.After field practice,the effect of waterflooding development is improved obviously,and recovery efficiency in2 creased from42.2%to45.98%.It is proved that the above measures are effective on tapping the poten2 tial in the complex fault block reservoirs by field ap2 plication.
K ey w ords:complex fault block reservoir,extra high water cut period,remaining oil distribution,side2 track drilling,hydrodynamics,recovery efficiency Chang F angrui,You Xingyou.R esearch and f ield application of highνpressure packing technique. PGRE,2002,9(4):77～79
Highνpressure packing is a technique in which screen casing annulus is sealed by highνpressure packing tool of PFS,and sand control system with good continuity,high permeability and big radius will be formed at high pressure using fracturing truck group and right sandνtransport ratio in integrated formation gravel input and wire wrapped liner outside packing.The application in Gudong oilfield proves
精密输液器Ⅴ
Vol.9No.4 Petroleum G eology and Recovery Efficiency:English Abstracts
