高含硫气藏数值模拟方法研究.pptxVIP

高含硫气藏数值模拟方法研究;一;超深、高温、常压、高含硫、边水，礁滩相储层;礁、滩体空间连续性差，储层构型复杂，不同构型单元物性差异大

礁、滩体非均质性强，地震响应特征复杂，表征礁、滩体内部结构难度大；;（2）高含硫气藏硫析出对开发的影响尚不清楚;基于构型的礁滩相高含硫气藏地质建模与数值模拟研究思路;一;1、礁滩相储层建模思路;元坝205;15KM;结构单元;储层垂向发育概率;通过单井相与地震阻抗数据体相关性分析，确定不同礁带储层阻抗体截断值，利用储层-阻抗相关性，建立储层发育地震概率体。;;物性概率统计;Por;建立的三维属性模型精细表征了礁滩体孔渗饱参数的非均质性特征；储层非均质性与气水分布特征符合礁滩体内部结构和储层发育特征。

礁体内部孔隙度分布符合储层发育规律：礁盖礁核礁基，礁后礁前;;利用YB273-1H水平井进行模型验证，构造、沉积相和储层物性分布基本与实钻一致。;目录;;井口气样中单质硫含量（1.1g/m3）小于地层条件下饱和硫含量，说明初始情况下，地层中未饱和单质硫，硫以游离态或化合态溶解在气相中，对初始渗透率没有影响。

建立了考虑液硫-酸气体系平衡的硫溶解度方程，计算了高温高压条件下酸气PVT相图。确定液硫析出条件（150℃、24.8MPa）。;构建新的酸性气体气-液-固三相状态方程;实现了酸性气体中硫析出的快速闪蒸计算，为酸性气藏的组分模拟奠定理论基础。;饱和液硫时液硫达到临界流动饱和度后可以流动，气液两相同流。

气体和液硫两相共渗区很窄，液硫临界流动饱和度较大。

气体流速对两相渗流曲线特征有影响（具有Hunderson效应）。;气体流速越快，气相及液相相对渗透率越高。（Hunderson，1994）

随着驱替压差的增大，气体流速增加，液硫以液滴形式被高速气体携带。;建立了酸性气藏气-液硫两相渗流模型，应用于数值模拟;建立了考虑液硫析出的酸气开采不同阶段流动模式，为渗流数学模型建立奠定基础;建立了液硫速度敏感相渗模型;;1）酸气渗流数学模型的改进;现有模拟器;网格文件;开发过程中井筒附近液硫饱和度分布变化;元坝204-1H井液硫分布分布图;目录;1、???成了考虑液硫析出的等效数值模拟方法，提出了液硫析出后通过重复改造改善开发效果的对策;利用等效模型预测，混合气产能40亿方，稳产6年

通过重复改造，可改善液硫析出后渗流环境，提高气井开发效果。;通过大湾固硫析出过程分析，建立储层温压条件下固硫沉积三个模型并应用于数值模拟模型中。;;;;;2、高含硫数值模拟下一步发展方向;普光硫沉积逐步向井筒深部及储层扩展，为保障气井正常生产，必须攻克井筒-储层硫沉积一体化预测及防治工艺技术难题，其中井筒硫沉积规律研究及建模是十分重要的一环。因此亟需研究针对硫沉积问题的准确预测方法，并与油藏数值模拟手段相结合，为实现高含硫气藏储层-井筒硫沉积预测技术的实现提供技术支撑。;2.井筒含硫颗粒多相流问题具有多尺度、多流态、多相态、多场耦合等特点，是一个十分复杂的传热传质过程。针对这一“卡脖子”问题的研究尚存在一些困难和未解决的关键技术：

;解决方案与技术路线;解决方案与技术路线;;精细数值模拟手段的优势

经济。无需搭建高昂实验室设备，时间成本和预算大幅降低。

高效。基于超算平台，可同时运行大量不同工况算例。

准确。基于先进数值方法，精细还原多相流界面破碎、聚并过程。;三;三;三;三;三;三;全井尺度宏观建模研究;1.气液控制方程建模—B-N七方程系统回顾

七方程表示方程求解的七个变量，包括气相体积分数、气相/液相的密度、速度、压力。

;气液控制方程建模—混相方程

水近似不可压，气体强可压缩，方程组特性迥异，直接采用七方程求解存在问题和困难；

一维宏观尺度条件下，可假设气液压力、温度平衡，可依此假设考虑混相运动方程，简化系统；

引入水-气漂移速度模型，考虑水气两相间存在速度分异的机制，求得混相速度后，可分解得到两相速度。;1.气液控制方程建模—混相方程组的简化

混相方程本质上是一个压力-速度隐式强耦合的系统，基于现有解耦方法求解要求单步内多次迭代使持液率方程和混相质量方程在数值层面保证一致性，否则将产生虚假的质量源造成计算崩溃；

进一步考虑稳健求解策略，首先求解气相方程，再更新水相速度及持液率分布。;2.硫输运方程建模;3.亚尺度模型嵌入;4.混相方程组的离散与求解

提出了一套半隐式迭代算法对速度-压力-温度-输运系统进行解耦。;三;4.混相方程组的离散与求解

采用交错网格技术，避免棋盘模式振荡；在网格界面存储离散速度，在网格中心存储压力、温度等变量

界面处通量构造：（1）一阶迎风格式；（2）二阶TVD格式;三;三;三;三;4.普光气井实际应用;