页岩气开采中的井底压力并行计算.pdfVIP

第31卷第3期 计 算 力 学 学 报 Vo1．31，No．3 2014年 6月 ChineseJournalofComputationalMechanics June2014 文章编号 ：1007～4708(2014)03—0396—06 页岩气开采中的井底压力并行计算 牛 聪 ， 黄 涛， 王 磊 ， 郑海军， 卢德唐 (中国科学技术大学渗流力学实验室 近代力学系，合肥 230027) 摘 要：根据页岩气流动特点建立了考虑混合气体高压物性参数、渗透率与孔隙度随压力变化的页岩气流动方 程．通过定义拟压力函数将页岩气流动的偏微分方程线性化。针对页岩气开发采用水平井多段压裂技术 ，采用 Newman乘积原理得到地层拟压力流动方程 的解析解表达式。依据解析解 的特征将解析解分解成适合并行计算 的无限求和及积分形式，提出了一套基于CUDA 的页岩气地层压力算法，将地层拟压力函数解析解划分为多个 并行度较高的步骤，利用GPU的并行计算能力，设计每个步骤的CUDA核函数，在英特尔i3540CPU(3．07GHz 主频 ，4GB内存)和 NVIDIA 的GTX550显卡上 ，计算 了页岩气的井底压力，分析 了井底压力特征。实验结果表 明，页岩吸附影响曲线变化剧烈程度，而扩散主要影响曲线发生变化的时间，在 GPU上的页岩气压力计算可达近 40倍 的加速 比。 关键词 ：页岩气 ；拟压力；GPU计算 ；流场计算 中图分类号 ：0359．i 文献标志码 ：A doi：10．7511／jslx2O14O3O19 传导水平井多段裂缝 中井底压力表达式 ，分析了裂 l 引 言 缝条数及矩形区域对井底压力的影响。在 Aissa模 页岩储层是 由天然裂缝或人工诱导裂缝和基 型中，每条裂缝半长及位置可变 ；由于 Laplace空 质构成的典型双重介质系统，基质是页岩气的主要 间上的解非常复杂，Aissa方法的计算速度较慢 。 储集空间，而裂缝是主要的流动通道 l【]。页岩气 2010年，Ozkan[6研究了裂缝沿水平井分布的 在页岩中有其特殊的赋存运移机理 ，与常规气藏最 水平井多段压裂井底压力计算 问题 。通过对井底 主要的区别在于页岩中不仅有游离气，还存在赋存 压力分析 ，Ozkan给出了水平井多段压裂可能会 出 于页岩的基质孔隙中的吸附气 ，其流人生产井筒需 现裂缝径 向流、系统线性流、径 向一线性或双线性 要经历三个过程[3“]：(1)在存在压降时，裂缝系统 流、地层线性流、拟径向流、系统线性流和系统径 向 中的页岩气流向生产井筒并且基质系统 中的页岩 流六种流动，并分析了每种流态的压力及导数形状。 气在基质表面进行解析；(2)在浓度差的作用下，页 由于 Aissa和 Ozkan方法中，每条裂缝半长都 岩气由基质系统向裂缝系统进行扩散 ；(3)在流动 可不相同，而裂缝之间的距离也可各不相 同，计算 势的作用下 ，页岩气通过裂缝系统流向生产井筒 。 量大 ，计算速度很慢。因此，目前商业软件一般采 页岩气 的开采主要采用水平井多段压力技术， 用 Horner和 Temeng模型 ]，模型中假设每条裂 对这类问题的求解可使用点源叠加得到源汇解 ，使 缝半长相同且裂缝间距也相同，这时可采用一条裂 用 Laplace变换考虑井筒存储及表皮 的影响，再使 缝及边缝解叠加得到水平井多段压裂的井底压力 ， 用数值 Laplace变换得到井底压力随时间变化 的 计算速度可大大提高。Ragharan和 Cipo