含泥质夹层油藏蒸汽注采过程的自适应网格算法研究
|
摘要
根据泥质夹层的低渗特性及空间分布,本文提出了一种含泥质夹层油藏网格渗透率的粗化计算方法,并在此基础上,将自适应网格算法应用于含泥质夹层油藏的数值模拟,提升其计算效率.在计算过程中,网格的动态划分仅依据流体物理量的变化,泥质夹层区域不全部采用细网格,仅针对流动锋面处的泥质夹层采用细网格,其余泥质夹层处采用不同程度的粗网格.相较于传统算法,网格数大幅下降.数值算例表明,自适应网格算法的计算结果精度与全精细网格一致,能够准确模拟出泥质夹层对于流体的阻碍作用,同时计算效率得到大幅提升,约为全精细网格算法的3～7倍.
Considering the low permeability and spatial distribution of muddy intercalation, a coarsening algorithm for grid permeability is proposed in this article. Based on this, the adaptive mesh refinement(AMR) algorithm is applied to the numerical simulation for the reservoir with muddy intercalations. In the AMR method, refining the dynamic grids or not only depends on the local gradient magnitude of the relevant physical quantities. It means that though fine grids are used, but not for all the muddy intercalation area. Only at the muddy intercalation area where the flow front passes over, the fine grids will be used locally. At the other muddy intercalation area, different levels of coarse grids are retained. In this way, the total number of grids is greatly reduced.Numerical examples show that the calculation results of the modified AMR method are consistent with the full fine grid method.The hindering effect of muddy intercalation on the flow can be presented accurately, and at the same time, the computational efficiency is improved by approximately 3～7 times compared to the traditional full fine grid method.
Considering the low permeability and spatial distribution of muddy intercalation, a coarsening algorithm for grid permeability is proposed in this article. Based on this, the adaptive mesh refinement(AMR) algorithm is applied to the numerical simulation for the reservoir with muddy intercalations. In the AMR method, refining the dynamic grids or not only depends on the local gradient magnitude of the relevant physical quantities. It means that though fine grids are used, but not for all the muddy intercalation area. Only at the muddy intercalation area where the flow front passes over, the fine grids will be used locally. At the other muddy intercalation area, different levels of coarse grids are retained. In this way, the total number of grids is greatly reduced.Numerical examples show that the calculation results of the modified AMR method are consistent with the full fine grid method.The hindering effect of muddy intercalation on the flow can be presented accurately, and at the same time, the computational efficiency is improved by approximately 3～7 times compared to the traditional full fine grid method.
引文
[1]张吉,张烈辉,胡书勇.陆相碎屑岩储层隔夹层成因、特征及其识别[J].测井技术,2003,27(3):221-224.
[2]徐樟有,熊琦华.一种定量研究薄泥质夹层的分布特征的方法[J].石油大学学报,1992,16(3):8-12.
[3]丁世梅,季民.泥岩隔夹层类型及剩余油控制研究[J].江汉石油学院学报,2004,26(4):130-134.
[4]朱东亚,曹学伟.临南油田隔层类型划分及其分布规律研究[J].地球科学-中国地质大学学报,2004,29(2):211-218.
[5]谢俊,张金亮,梁会珍,等.濮城油田末端扇储层隔夹层成因及分布特征[J].中国海洋大学学报,2008,38(4):653-656.
[6]李明川,姚军,葛家理.渗流力学进展与前沿[J].力学季刊,2012,33(1):74-80.
[7]孙红霞,姚军,曹志伟,等.一种低数值耗散的求解油水驱替过程的数值方法[J].力学季刊,2017,38(4):629-640.
[8]施安峰,王晓宏,贾江涛,等.复杂断层油藏蒸汽注采过程的自适应网格法[J].华中科技大学学报,2012,40(3):118-122.
[9]BERGER M J,OLIGER J.Adaptive mesh refinement for hyperbolic partial differential equations[J].Journal of Computational Physics,1984,53(3):484-512.
[10]BERGER M J,COLELLA P.Local adaptive mesh refinement for shock hydro-dynamics[J].Journal of Computational Physics,1989,82(1):64-84.
[11]WANG X H,QUINTARD M,DARCHE G.Adaptive mesh refinement for one-dimensional three-phase flow with phase change in porous media[J].Numerical Heat Transfer B,2006,50(3):231-268.
[12]王晓宏.蒸汽热采稠油数值模拟的自适应网格法[C]//中国工程热物理学会多相流学术会议,2006,564-573.
[13]罗海山,王晓宏,施安峰.非均匀裂缝性油藏蒸汽驱的自适应网格计算[J].辽宁工程技术大学学报,2009,28(S1):30-33.
[14]唐明云,施安峰,王晓宏.含轻质组分稠油油藏蒸汽注采过程自适应网格法[J].力学季刊,2013,34(1):25-31.
[15]KING P R.The use of renormalization for calculating effective permeability[J].Transport in Porous Media,1989,4(1):37-58.
[2]徐樟有,熊琦华.一种定量研究薄泥质夹层的分布特征的方法[J].石油大学学报,1992,16(3):8-12.
[3]丁世梅,季民.泥岩隔夹层类型及剩余油控制研究[J].江汉石油学院学报,2004,26(4):130-134.
[4]朱东亚,曹学伟.临南油田隔层类型划分及其分布规律研究[J].地球科学-中国地质大学学报,2004,29(2):211-218.
[5]谢俊,张金亮,梁会珍,等.濮城油田末端扇储层隔夹层成因及分布特征[J].中国海洋大学学报,2008,38(4):653-656.
[6]李明川,姚军,葛家理.渗流力学进展与前沿[J].力学季刊,2012,33(1):74-80.
[7]孙红霞,姚军,曹志伟,等.一种低数值耗散的求解油水驱替过程的数值方法[J].力学季刊,2017,38(4):629-640.
[8]施安峰,王晓宏,贾江涛,等.复杂断层油藏蒸汽注采过程的自适应网格法[J].华中科技大学学报,2012,40(3):118-122.
[9]BERGER M J,OLIGER J.Adaptive mesh refinement for hyperbolic partial differential equations[J].Journal of Computational Physics,1984,53(3):484-512.
[10]BERGER M J,COLELLA P.Local adaptive mesh refinement for shock hydro-dynamics[J].Journal of Computational Physics,1989,82(1):64-84.
[11]WANG X H,QUINTARD M,DARCHE G.Adaptive mesh refinement for one-dimensional three-phase flow with phase change in porous media[J].Numerical Heat Transfer B,2006,50(3):231-268.
[12]王晓宏.蒸汽热采稠油数值模拟的自适应网格法[C]//中国工程热物理学会多相流学术会议,2006,564-573.
[13]罗海山,王晓宏,施安峰.非均匀裂缝性油藏蒸汽驱的自适应网格计算[J].辽宁工程技术大学学报,2009,28(S1):30-33.
[14]唐明云,施安峰,王晓宏.含轻质组分稠油油藏蒸汽注采过程自适应网格法[J].力学季刊,2013,34(1):25-31.
[15]KING P R.The use of renormalization for calculating effective permeability[J].Transport in Porous Media,1989,4(1):37-58.