刘骞
中图分类号：TE357 文献标识码：A 文章编号：1009-914X（2017）02-0129-02
油藏数值模拟中，人工压裂缝的处理方法主要分为两种：网格加密和近井模拟（Near Wellbore Module）。但是由于近井模拟处理起来较繁琐，在模型中计算起来较慢且容易不收敛，所以目前人工裂缝的模拟主要是通过局部网格加密（Local Grid Refinemnt）的方法来完成[1]。局部网格加密的方法很多，有基于笛卡尔网格的加密方法，有基于非结构网格的加密方法，如Pebi网格加密[2-6]。
笛卡尔网格的加密方法也分为两种，一种是运用不均匀加密网格（HXFIN），另一种是利用不均匀网格（DX）。以往的研究表明，将裂缝所在的网格用“等效导流能力”法处理，在“缝宽”小于1.0m的情况下井的产量变化不大，该方法实用可靠。
本文分别用笛卡尔网格的加密两种方法模拟裂缝，比较了它们的区别。
一、不均匀加密网格
在不均匀加密网格方法中，我们设置的模型维数为110*40*10，DX=DY=DZ=5m，模型中心顺着X方向有一水平井生产（图1-图2）。我们分别在X方向第10列、20列、30列、40列、50列、60列、70列、80列、90列、100列网格中分别进行不均匀加密，这样在每一列加密的网格中加密出一条0.5m宽的网格，我们用这列0.5m宽的网格来模拟分段压裂的裂缝（图）。这样就形成一个长度为550m，等间距分段压裂了10段0.5m宽裂缝的水平井在生产（图3）。
二、不均匀网格
在不均匀网格方法中，我们设置的模型维数为119*41*10，除了X方向第10、21、32、43、54、65、76、87、98、109个网格DX=0.5外，其他的DX=5m（图4-图5）。所有的DY=DZ=5m。这样也形成一个长度为550m，等间距分段压裂了10段0.5m寬的裂缝的水平井在生产（图6-图7）。
模拟过程中，两种方法所用的基本参数都相同，参数取值如表1：
两种方式均以定液（50t）方式生产，效果如图8-图11：
三、结论
1）虽然两种方法都是比较常用的裂缝模拟方法，但是结果差别较大；
2）不均匀加密网格方法中，水平井生产140天，日产油由50t/d缓慢递减至0t/d；
3）不均匀网格方法中，水平井生产140天，前50天保持日产油50t/d不变，之后10天日产油迅速降至0t/d；
4）两种方法虽然日产油不同，但是生产140天后的累产油相同都为2590m3。
5）在计算时间上，由于不均匀网格方法的网格数少，计算速度要远远快于不均匀加密网格方法。
参考文献
[1] 韩大匡.油藏数值模拟基础[M].北京：石油工业出版社，1993.HAN Da-kuang.Reservoir Simulation Foundamental[M].Beijing：Petroleum Industry Press，1993.
[2] 刘立明，廖新维，陈钦雷.混合PEBI网格精细油藏数值模拟应用研究[J].石油学报，2003，24（3）：64-67.LIU Li-ming，LIAO Xin-wei，CHEN Qin-lei.Usage of hy-brid PEBI grid in fine reservoir simulation[J].ActaPetroiei Sinica，2003，24（3）：64-67.
[3] 蔡强，杨钦，孟宪海，等.二维PEBI网格的生成[J].工程图学学报，2005（2）：69-72.CAI Qiang，YANG Qin，MENG Xian-hai，et al.Research on2D PEBI grid generation[J].Journal of EngineeringGraphics，2005（2）：69-72.
[4] 安永生，吴晓东，韩国庆.基于混合PEBI网格的复杂井数值模拟应用研究[J].石油学报，2007，31（6）：60-63.AN Yong-sheng，WU Xiao-dong，HAN Guo-qing.Application of numerical simulation of complex well based on hy-brid PEBI grid[J].Acta Petroiei Sinica，2007，31（6）：60-63.
[5] 向祖平，张烈辉，陈中华，等.油藏任意约束平面域PEBI网格的生成算法[J].西南石油学院学报，2006，28（2）：32-35.XIANG Zu-ping，ZHANG Lie-hui，CHEN Zhong-hua，et al.Algorithm for constructing PEBI mesh of arbitrarily shapedand constrained planar domains in oil reservoir[J].Journal of Southwest Petroleum Institute，2006，28（2）：32-35.
[6] Patrick W von Pattay，Leonhard J Ganzer.Reservoir simulation model for fractured and partially fractured reservoirs based on PEBI grids[C].SPE 66384，2001.
