缝洞型油藏溶洞内水驱流态研究.docx
( 4.5 )《缝洞型油藏溶洞内水驱流态研究.docx》由会员分享，可在线阅读，更多相关《缝洞型油藏溶洞内水驱流态研究.docx（14页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、缝洞型油藏溶洞内水驱流态研究 工业科技第35卷 缝洞型油藏由于储集体类型多样，微孔隙、微裂缝、裂缝和中型溶洞以及大型溶洞、大型裂缝的分布，在尺度和空间结构上都具有很大的差异，流体流淌类型困难，渗流、裂缝流和洞窟流都有存在12-15。赵欢等16利用油藏数值模拟方法，探讨了不同聚合物参数对“二三结合”开发模式原射孔层和补孔层的压力、含水饱和度等渗流特征的影响。刘金玉17利用COMSOL Multiphysics多物理场软件对单缝概念模型进行了水驱油数值模拟，分析了油水流态、波及面积、水驱前缘及剩余油分布的基本特征。目前对缝洞型碳酸盐岩油藏水驱的探讨多是从渗流角度动身，但是由于缝洞型油藏内部结构困难
2、，裂缝和溶洞的存在状况无法精确驾驭，特殊是针对一些大裂缝和大溶洞，内部属于未填充状态，多孔介质没有在内部分布，所以内部流体流淌不属于渗流范畴，它的流淌规律符合Navier-Stokes方程，因此需采纳流体力学对其进行探讨。 本文从流体力学角度动身，针对未填充缝洞型碳酸盐岩油藏进行了油水两相流淌模拟，探讨单裂缝-溶洞中流体流淌时油藏内原油黏度、注水速度对油水界面、波及面积的影响。 1数学模型 1.1连续方程 连续性方程是流体运动学的基本方程，是质量守恒原理的流体力学表达式。 连续方程微分形式见式。 t+=0， 式中：表示流体密度，kg/s；t表示时间，s；u表示流速，m/s。 模型设定的油水两相
3、流淌为稳态流淌，即： =0 。 油水两相都是不行压缩流体，即： t=0 。 将以上条件代入连续方程，得到油水两相流的连续方程表达式： u=0 。 式即用来描述缝洞型介质两相流淌机理的连续方程，可以通过连续方程得到流体的流速，代入纳维-斯托克斯方程。 1.2纳维-斯托克斯方程 纳维-斯托克斯方程简称N-S方程，用来描述黏性不行压缩流体动量守恒的运动方程18。 溶洞型油藏中油水两相流淌属于黏性流体的不行压缩流淌，由于须要探讨油水界面在不同黏度原油和不同入口速度下的推动状况，所以采纳张量形式的不行压缩黏性N-S方程，即： duidt=Fipijxj+2uix2j， 式中：duidt表示单位体积上的惯
4、性力；Fipijxj表示单位体积上的质量力；2uix2j表示单位体积上的黏滞力。 由于油水两项在大溶洞和大裂缝中的运动属于自由流淌，遵循N-S定律，同时为了结合VOF方法追踪溶洞中油水界面问题，所以采纳张量形式的不行压缩黏性N-S方程，用于探讨未填充溶洞中的油水两相流淌问题。 1.3VOF方法及原理 VOF方法基本原理是通过对网格单元中流体之间的体积比函数F来确定自有界面，追踪界面的推动状况。由Hirt和Nichols首次提出后经NOH等改进，在保证质量守恒的前提下，PUCKETT等采纳方向分裂算法对输送方程进行优化，利用流体体积函数方程张量表示方式如下： Ft+ujFxj=0。 由于油的密度
5、小于水的密度，两相之间存在明显界面，所以采纳VOF模型对油水界面进行跟踪。依据单位网格内油水体积比F得到油水界面形态和位置。当F=0时，说明单元内全部是油相；当F=1时，说明单元内全部是水相，已经将油驱出；当0F1时，说明单元内存在油水两相，有油水界面存在。 依据流体体积函数方程张量表示方式，通过流体流淌时间、速度以及流淌方向上的位置可以计算出单位网格内油水体积比，利用瞬态计算方法，得到不同时刻的油水界面分布状况。 2数值计算 2.1初始条件 单裂缝是缝洞型油藏的基本流淌单元，对流体在油藏中的流淌起确定性作用。假设油和水都是不行压缩流体，密度为常数。由于溶洞三度空间长度相近，同时为了减小计算量
6、，将溶洞简化为正方形二维模型，出入口位于左右兩侧中心，左端为出口，右端为入口。入口定义为速度入口，出口定义为压力出口。 2.2模型参数及网格划分 参数确定：直径大于101 mm属于巨洞；直径在10101 mm属于大洞；直径在510 mm之间属于中洞；直径在25 mm之间属于小洞。由于溶洞空隙三度空间相近，且都大于2103 m，所以选用三度尺寸均为50 mm的溶洞，缝宽为0.6 mm。实际缝洞型碳酸盐岩油藏地层原油黏度一般为230 mPas，所以选用原油黏度为5，15和25 mPas。原油密度选取860 kg/m3，水的密度为1016 kg/m3。 溶洞模型为边长50 mm的正方形二维模型，出入
7、口缝宽为0.6 mm，利用gambit对模型进行六面体网格划分，见图1。依据雷诺数计算，采纳Fluent中的层流模型，且为了便利追踪油水界面，采纳VOF方法进行瞬态计算，收敛精度为0.000 01。 首先，计算入口速度为0.08 m/s、黏度为5 mPas的溶洞内油藏的流态，得到水驱油水分界面分布，波及面积及剩余油分布，依次变更黏度为15和25 mPas，得到不同黏度下的油水界面发展状况、波及面积及剩余油分布状况。变更原油黏度，得到相同入口速度下不同黏度的油水分布状况，分析原油黏度和流态之间关系；再变更入口速度，依次设定入口水速为0.08，0.16，0.24 m/s，得到原油黏度为15 mPa
8、s时不同入口速度下的油水分布状况，分析入口水速和流态之间的关系。 3计算结果和分析 3.1肯定入口速度下黏度改变对波及面积的影响 通过图2可以看到在入口速度为0.08 m/s时完整的水驱过程波及面积改变，黏度对单缝洞型油藏中水驱波及面积的影响很大：在黏度为5 mPas时，水相趋向中心线，向洞口流淌现象明显，油藏内形成漩涡，波及面积最小，水相占据比例为27%；当黏度增大到15 mPas时，起先向活塞驱替转变，没有形成漩涡，波及面积增大，水相占据比例为52%；当黏度增大到25 mPas时，呈现水滴状，波及面积达到最大，水相占据比例为61%。 出现这种现象的缘由：在水驱入口速度为0.08 m/s、原
9、油黏度从5 mPas增大到25 mPas时，流淌性变差，对注入水的扩散阻力增大，相间磨损增加，水相推动速度变小，油水界面到达出口的时间变长，扩散作用变大；同时，由于原油黏度增大，油水界面对前推动的速度减慢，水相沿中心线纵向扩散趋势渐渐增加。因此，类射流现象随着原油黏度增加渐渐转换成活塞驱替，漩涡渐渐消逝，波及面积渐渐增大。 3.2肯定入口速度下黏度改变对油水界面的影响 为探讨水驱油过程中油水界面的发展形态，从上述3种黏度下分别取6幅不同时刻的流态图，见图3。 由图3可知，未到达出口前，3种黏度下的油水界面都是根据匀称对称的方式向前推动，到达出口后，由于压力改变，油水界面会产生“尖嘴状”。当黏度
10、为5 mPas时，水流受原油阻力较小，水相优先于中线位置通过，由于水相扩散作用产生漩涡，部分原油受到漩涡波及，在漩涡作用下向右移动，水相到达出口时间最短，耗时14 s；黏度增大到15 mPas时，水相受到的阻力增大，水相沿垂直中线方向扩散趋势增加，产生漩涡趋势减弱，油水界面匀称向前推动，水相到达出口的时间变长，耗时17 s；黏度增大到25 mPas时，在到达出口前，油水界面前缘以圆形匀称前进，靠近出口时，前缘呈现“尖嘴状”，水相到达出口的时间最长，耗时19 s。 3.3肯定黏度下入口速度改变对波及面积的影响 由图4可知，入口速度对波及面积影响显著，油水界面到达右侧缝洞出口时，波及面积随着入口速
11、度增加而减小。通过图4中完整的水驱过程波及面积图可知，入口速度对单缝洞型油藏中水驱波及面积的影响很大：在速度为0.08 m/s时，波及面积呈现水滴状，油藏内无漩涡产生，波及面积最大，水相占据比例为59%；当入口速度为016 m/s时，水相流淌有向中心线靠近的趋势，漩涡起先形成，波及面积变小，水相占据比例为43%；当入口速度为024 m/s时，射流现象明显，产生漩涡，波及面积最小，油水界面推动速度明显下降。 出现这种现象的缘由：在原油黏度为15 mPas时，水驱入口速度从0.08 m/s增加到0.24 m/s，由于注入速度与压力梯度成正比，随着入口速度增加，压力梯度增大，油水界面对前推动的速度随
12、之增加，水相推动速度变大，油水界面到达出口时间变短，从18 s削减到8 s再削减到4 s，相间扩散作用变小；同时，由于入口速度增大，水相更简单沿孔洞中间区域到达右侧出口端。因此，波及面积随着注入速度增加活塞驱替向射流发展，漩涡渐渐形成，波及面积渐渐下降。 3.4肯定黏度下入口速度改变对油水界面的影响 为了探讨水驱油过程中油水界面的发展形态，从上述3种速度下分别取6幅不同时刻的流态图，见图5。 由图5可知，未到达出口前，3种速度下的油水界面都是根据匀称对称的方式向前推动，到达出口旁边时，由于压力改变，油水界面会产生“尖嘴状”。当速度为0.08 m/s时，水相优先于中线位置通过现象并不明显，油水界
13、面到达出口时间为17 s；当速度为0.16 m/s时，水相沿中线向右侧出口推动的趋势增加，漩涡起先产生，油水界面匀称向前推动，到达出口时间为8.6 s；当速度为0.24 m/s时，油水界面以圆形匀称前进，到达出口时间为4 s，明显产生漩涡，部分原油被卷入漩涡中。 3.5精确性论证 模型采纳两套网格对比，对网格无关性进行分析，第1套网格数量为16 283、节点数量为16 024，采纳四面体结构网格。第2套网格数量为52 737、节点数量为54 238，采纳三角形网格。将两套网格条件下模型的垂直中心轴線的铅垂线上的油相体积分布进行对比和模拟验证，见图6。由图6可知，平均相符程度为101.33%，可
14、以认为两套网格的模拟结果相同，能做到网格无关性。基于运算量考虑，第1套网格更小，运算更快，因此后续运用第1套结构网格进行数值模拟分析。 4结论 在单缝洞型油藏水驱开采过程中，原油黏度和注水速度的变更会在很大程度上影响波及面积和油水界面，从而影响水驱效率。 1）相同注水速度下，黏度越大，油水界面推动速率越低，原油黏度较低时会在出口产生回流，奢侈能量达不到水驱预期效果，波及面积随着黏度增大渐渐增大，推动速度降低。 2）原油黏度肯定时，波及面积随着入口速度的增加而减小，油水界面呈对称趋势匀称推动，水速较大时会发生类射流现象并有漩涡产生，随着入口水速增大，类射流现象和漩涡现象渐渐明显。 3）由于水驱速
- 配套讲稿：
如PPT文件的首页显示word图标，表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制：
部分文档作品中含有的国旗、国徽等图片，仅作为作品整体效果示例展示，禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词：
- 缝洞型 油藏 溶洞 内水驱流态 研究
淘文阁 - 分享文档赚钱的网站所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。
限制150内