注水开发储层物性参数变化数值模拟研究
|
摘要
老油田高含水期取心资料和矿场测井动态资料表明,储层渗透率、原油粘度和岩石润湿性等物性参数随注水开发发生连续变化。现有的油藏数值模拟软件没有考虑这些参数的动态变化,这与油田开发实际不符,导致预测的结果和矿场实际结果存在一定偏差。因此,系统全面地研究开发因素引起储层物性参数的变化规律,建立符合地下实际情况的动态地质模型进行精细油藏数值模拟,准确预测高含水期剩余油的分布规律,在当前油气资源日益紧张的形势下对指导油田的进一步挖潜、提高采收率具有重要意义。
本文在全面总结储层物性参数变化原因的基础上,重点分析了注入水动力作用对储层物性参数变化的影响,基于储层物性参数变化与注入水冲刷孔隙体积倍数之间关系建立了考虑渗透率、原油粘度及润湿性连续变化的时变描述模型,并将其引入常规三维三相黑油模型中,建立参数时变数学模型。利用有限差分方法将该数学模型离散化得到数值模型,应用IMPES方法进行求解,在此基础上应用FORTRAN 90和Visual C++两种计算机语言开发了参数时变数值模拟软件,弥补了目前商业油藏数值模拟器的不足,提高了油藏数值模拟的精度。
应用研制的可考虑储层参数连续变化的参数时变数值模拟软件,对参数时变概念模型和参数不变概念模型进行数值模拟研究,将模拟结果进行对比,分析储层各参数动态变化对开发效果的影响。研究结果表明:储层物性参数的变化对水驱开发效果影响显著;储层渗透率增大及原油粘度增大均导致水驱开发效果变差;润湿性的变化使得水驱开发效果变好。
Both old oilfield coring data in high water cut stage and mines dynamic logging data show that reservoir permeability, crude oil viscosity and rock wettability have changed greatly after long-term flood development compared with the early. But the change of these parameters is not considered in existing reservoir numerical simulation software. It’s innconsistent with the actual oilfield development and will cause prediction results of poor accuracy which have big deviation to actual results. Therefore, it has vital significance to clear the distribution of remaining oil in high water cut stage to guide further oilfield development and enhance oil recovery in the increasingly tense situation of oil and gas resources cullently, which bases on establishment dynamic geological model conforming to the actual situation of the underground for fine reservoir numerical simulation after systematic and comprehensive research on reservoir parameters change rules due to development factors.
The effect of injection water power function on parameters change has been focused on in this paper. After that, parameters time-varying description models are established on the basis of the relationship between reservoir parameters change and flood water pore volume multiple and introduced in the mathematical model used in conventional three-dimensional and three-phase black oil model. Using finite difference methods, the mathematical model is discretized to numerical model, The numerical model is soluted by IMPES method. At last, Parameter time-variance numerical simulation software is compiled by using FORTRAN 90 and Visual C++ program languages. This software makes up for the current commercial reservoir numerical simulator and improves the precision of the reservoir numerical simulation.
Research shows that the change of reservoir parameters has significant effects on the flood development effect: the increase of reservoir permeability and crude oil viscosity causes water flooding development to become worse; the rock wettability change generates a better effect.
本文在全面总结储层物性参数变化原因的基础上,重点分析了注入水动力作用对储层物性参数变化的影响,基于储层物性参数变化与注入水冲刷孔隙体积倍数之间关系建立了考虑渗透率、原油粘度及润湿性连续变化的时变描述模型,并将其引入常规三维三相黑油模型中,建立参数时变数学模型。利用有限差分方法将该数学模型离散化得到数值模型,应用IMPES方法进行求解,在此基础上应用FORTRAN 90和Visual C++两种计算机语言开发了参数时变数值模拟软件,弥补了目前商业油藏数值模拟器的不足,提高了油藏数值模拟的精度。
应用研制的可考虑储层参数连续变化的参数时变数值模拟软件,对参数时变概念模型和参数不变概念模型进行数值模拟研究,将模拟结果进行对比,分析储层各参数动态变化对开发效果的影响。研究结果表明:储层物性参数的变化对水驱开发效果影响显著;储层渗透率增大及原油粘度增大均导致水驱开发效果变差;润湿性的变化使得水驱开发效果变好。
Both old oilfield coring data in high water cut stage and mines dynamic logging data show that reservoir permeability, crude oil viscosity and rock wettability have changed greatly after long-term flood development compared with the early. But the change of these parameters is not considered in existing reservoir numerical simulation software. It’s innconsistent with the actual oilfield development and will cause prediction results of poor accuracy which have big deviation to actual results. Therefore, it has vital significance to clear the distribution of remaining oil in high water cut stage to guide further oilfield development and enhance oil recovery in the increasingly tense situation of oil and gas resources cullently, which bases on establishment dynamic geological model conforming to the actual situation of the underground for fine reservoir numerical simulation after systematic and comprehensive research on reservoir parameters change rules due to development factors.
The effect of injection water power function on parameters change has been focused on in this paper. After that, parameters time-varying description models are established on the basis of the relationship between reservoir parameters change and flood water pore volume multiple and introduced in the mathematical model used in conventional three-dimensional and three-phase black oil model. Using finite difference methods, the mathematical model is discretized to numerical model, The numerical model is soluted by IMPES method. At last, Parameter time-variance numerical simulation software is compiled by using FORTRAN 90 and Visual C++ program languages. This software makes up for the current commercial reservoir numerical simulator and improves the precision of the reservoir numerical simulation.
Research shows that the change of reservoir parameters has significant effects on the flood development effect: the increase of reservoir permeability and crude oil viscosity causes water flooding development to become worse; the rock wettability change generates a better effect.
引文
[1]刘定珺,谢俊,张金亮.我国剩余油技术研究现状与进展[J].西北地质,2004,37 (4):1-6
[2]穆龙新.油藏描述的阶段性及特点[J].石油学报,2000,21(5):103-108
[3]大港油田科技丛书编委会.大港油田科技大港油田开发实践[M].北京:石油工业出版社,1999:279-289
[4]方宏长,冯明生.中国东部几个主要油田高含水期提高水驱采收率的方向[J].石油勘探与开发,1999,26(1):12-16
[5]孙焕泉.油藏动态模型和剩余油分布模式[D].北京:中国科学院地质与地球物理研究所,2002
[6]孙焕泉,孙国,吴素英,等.储集层参数动态地质模型的建立——以胜坨油田二区沙二段12层为例[J].石油勘探与开发,2004,31(增):89-91
[7]邓玉珍,徐守余.三角洲储层参数动态模型研究——以胜坨油田二区下第三系沙二段83小层为例[J].石油学报,2003,24(2):61-64
[8]宋万超,孙焕泉,孙国等.油藏开发流体动力地质作用[J].石油学报,2002;(3): 52-55
[9]李健,李红南.油藏开发流体动力地质作用——以胜坨油田二区为例[J].石油学报, 2002,23(3):52-55
[10]陈永生.油藏流场[M].北京:石油工业出版社,1998:26-77
[11]杜玉洪,张继春,侯翠芬.储层宏观物性参数随注水开发动态演化模式研究[J].特种油气藏,2004,11(5):52-55
[12] Y.Cinar, S.Marquez, F.M.Orr. Effect of IFT Variation and Wettability on Three-Phase Relative Permeability[J]. Society of Petroleum Engineers Inc., 2004.9: 26-29
[13] Faru.K.Civan. Characterization and prediction of formation damage in two-phase flow systems[R]. SPE: 25429, 1993
[14] Martin A.Fern, Arne Graue. Numerical Sinulation and Sensityvity Analysis of In-Situ Fluid Flow in MRI Laboratory Waterfloods of Fracturde Carbonate Rocks at Different Wettabilities. SPE: 116145, 2008
[15] Matthew D.Jackson, Per H.Valvatne, Martin J.Blunt. Prediction of wettability variation and its impact on waterflooding using pore-to reservoir-scale simulation[J]. SPE: 77543, 2002
[16] P.Behrenbruch, H.M.Goda. Two-Phase Relative Permeability Prediction:A Comparison of the Modified Brooks-Corey Methodology with a New Carman-Kozeny Based Flow Formulation[J]. SPE: 101150, 2006
[17] Zh.Hongxing, Liu Qingnan, Li Fengqin. Variations of Petrophysical Parameters After Sandstone Reservoirs Watered Out in Daqing Oil Field[J]. SPE: 30884, 1995
[18] LIU Xinghui, FARU K Civan.Characterization and prediction of formation damage in Two-phase flow systems[R]. SPE: 25429, 1993.
[19] Helge Holden, Kenneth Hvistendahl Karlsen, Knut Andreas Lie. Operator splitting methods for degenerate convection-diffusion equations II:Numerical examples with emphasis on reservoir simulation and sedimentation. Canadian Mathematical Society, Conference proceedings
[20]鞠斌山.油藏渗流系统物性变化机理与数学模型研究[D].北京:中国地质大学, 2006
[21]徐守余.渤海湾盆地胜坨油田二区储层微观渗流场演化研究[J].石油实验地质, 2003,20(1):62-65
[22]林光荣,陈付星,邵创国,等.马岭油田长期注水对油层孔隙结构的影响[J].西安石油学院学报(自然科学版),2001,24(2):31-35
[23]单华生,姚光庆,周锋德.储层水先后结构变化规律研究[J].海洋石油,2004,24 (1):62-66
[24]李兴国,孙春红,赵云飞.喇嘛甸油田储层微观结构特征变化规律研究[J].大庆石油地质与开发,2007,26(1):79-82
[25]郭莉,王延斌,刘伟新,等.大港油田注水开发过程中油藏参数变化规律分析[J],石油实验地质,2006,l28(1):85-90
[26]邓玉珍,吴素英.注水开发过程中储层物理特征变化规律研究[J].油气采收率技术, 1996,3(4):44-52
[27]张伟峰,刘守军,李拥安,等.孤岛馆陶组注水开发储层性质动态变化特征研究[J].地球科学与环境学报,2004,26(2):51-54
[28]刘志宏.胜二区沙1-2开发期储层变化特征及特高含水期改善水驱效果研究[D].成都:西南石油学院,2002
[29]龙雨丰,张晓晖,刘鸿,等.低渗油藏注水开发中储层参数变化研究与应用[J].江汉石油学院学报,2001,23(增):86-87
[30]杨晓蓓,冯毅,郭恩常,等.砂砾岩油藏特高含水期储层参数变化规律[J].西南石油学院学报,2000,22(4):22-25
[31]吴素英.胜利油田二区沙二段83层储层润湿性变化规律及对开发效果的影响[J].油气地质与采收率,2006,13(2):71-74
[32]孙焕泉,屈岩涛,房会春,等.胜利油区砂岩储集层敏感性特征研究[J].石油勘探与开发,2000,27(5):72-75
[33]李阳.陆相断陷湖盆油藏流场宏观参数变化规律及动态模型[J].石油学报,2005, 26(2):65-68
[34]陈亮,吴胜和,刘宇红.胡状集油田胡十二块注水开发过程中储层动态变化研究[J].1999,22(4):78-81
[35]李存贵,徐守余.长期注水开发油藏的孔隙结构变化规律[J].石油勘探与开发.2003, 30(2):94-96
[36] Aziz K, Settari A. Petroleum Reservoir Simulation[M]. Applied Science Publishers, 1979
[37]程林松,李丽,要春兰.考虑水驱油藏润湿性变化的数值模拟方法[J].水动力学研究与进展,2003,18(6):86-97
[38]高博禹,彭仕宓,黄述旺.胜坨油田二区沙二段3砂层组分阶段油藏数值模拟[J].石油勘探与开发,2004,31(6):82-84
[39]盖英杰,吕德灵,郭元灵.高含水期油藏分阶段数值模拟[J].油气采收率技术,2000, 7(1):54-56
[40]姜瑞忠,陈月明,邓玉珍,等.胜坨油田二区油藏物理特征参数变化数值模拟研究[J].油气采收率技术,1996,3(2):50-56
[41]崔传智,赵晓燕.考虑油藏储层参数变化的数值模拟研究.水动力学研究与进展, 2004,19(A):912-915
[42]高博禹,彭仕宓,黄述旺,等.高含水期油藏精细数值模拟研究[J].石油大学学报(自然科学版),2005,29(2):11-15
[43]谷建伟.高含水开发期基于微观渗流机理的宏观油藏数值模拟研究[J].中国海洋大学:2004
[44]刘红卫,刘义坤,潘晓军.萨中开发区考虑油藏润湿性影响数值模拟[J].内蒙古石油化工,2008,9:221-222
[45]张洪波,王志章,戴胜群,等.水驱前后油藏参数变化机理研究[J].石油天然气学报(江汉石油学院学报).2005,27(A):665-666
[46]王志章,蔡毅,杨蕾.开发中后期油藏参数变化规律及变化机理[M] .北京:石油工业出版社,1999:50-60
[47]姜瑞忠.储层特征参数变化对油藏开发效果的影响[J].油气田地面工程,2005, 24(4):32-33
[48]秦积舜,李爱芬.油层物理学[M].山东:石油大学出版社,2001:94-178
[49]吴素英,孙国,程会明,等.长期水驱砂岩油藏储层参数变化机理研究[J].油气地质与采收率,2004,11(2):9-11
[50]张大为,闫宝珍,曾照英.注水开发油藏油水相对渗透率变化规律[J].重庆大学学报(自然科学版),2000,23(增):77-78
[51]刘中云,曾庆辉,唐周怀,等.润湿性对采收率及相对渗透率的影响[J].2000,21 (2):940-942
[52]崔浩哲,姚光庆,周锋德.低渗透砂砾岩油层相对渗透率曲线的形态及其变化特征[J].地质科学情报,2003,22(1):87-91
[53]鞠斌山,樊太亮,张金川,等.水驱开发油藏原油黏度变化规律及其对开发效果的影响[J].石油勘探与开发,2006,33(1):99-102
[54]鞠斌山,王春田,李师涛,等.变原油黏度油水两相流驱替特征研究[J].西安石油学院(自然科学版),2005,18(1):17-20
[55]邓瑞健.濮城油田注入水体积流量试验特征[J].江汉石油学院学报,2003,25(3): 111-112
[56]孙士孝,韩锦文,郭云尧,等.胜坨油田岩石水洗前后物性变化的试验研究.1996, 20:33-35
[57]刘月田,徐明旺,彭道贵,等.各向异性渗透率油藏数值模拟[J].计算物理,2007, 24(3):295-300
[58]刘慧卿.油藏数值模拟方法专题[M].山东:石油大学出版社,2001:79-84
[59]皮斯曼D.W.油藏数值模拟基础[M].北京:石油工业出版社,1982:1-20
[60]韩大匡,陈钦雷,闫存章.油藏数值模拟基础[M].北京:石油工业出版社,1999: 30-40
[61]关治,陆金甫,等.数值分析基础[M].北京:高等教育出版社,2002:310-323
[62]李庆阳,关治,等.数值计算原理[M].北京:清华大学出版社,2000:170-236
[63]封建湖.数值分析原理[M].北京:科学技术出版社,2001:173~175
[64]李庆阳,关治.数值计算原理[M].北京:清华大学出版社,2000:171~226
[65]李维国,黄炳家,同登科,等.数值计算方法[M].山东东营:石油大学出版社, 2004:173~180
[2]穆龙新.油藏描述的阶段性及特点[J].石油学报,2000,21(5):103-108
[3]大港油田科技丛书编委会.大港油田科技大港油田开发实践[M].北京:石油工业出版社,1999:279-289
[4]方宏长,冯明生.中国东部几个主要油田高含水期提高水驱采收率的方向[J].石油勘探与开发,1999,26(1):12-16
[5]孙焕泉.油藏动态模型和剩余油分布模式[D].北京:中国科学院地质与地球物理研究所,2002
[6]孙焕泉,孙国,吴素英,等.储集层参数动态地质模型的建立——以胜坨油田二区沙二段12层为例[J].石油勘探与开发,2004,31(增):89-91
[7]邓玉珍,徐守余.三角洲储层参数动态模型研究——以胜坨油田二区下第三系沙二段83小层为例[J].石油学报,2003,24(2):61-64
[8]宋万超,孙焕泉,孙国等.油藏开发流体动力地质作用[J].石油学报,2002;(3): 52-55
[9]李健,李红南.油藏开发流体动力地质作用——以胜坨油田二区为例[J].石油学报, 2002,23(3):52-55
[10]陈永生.油藏流场[M].北京:石油工业出版社,1998:26-77
[11]杜玉洪,张继春,侯翠芬.储层宏观物性参数随注水开发动态演化模式研究[J].特种油气藏,2004,11(5):52-55
[12] Y.Cinar, S.Marquez, F.M.Orr. Effect of IFT Variation and Wettability on Three-Phase Relative Permeability[J]. Society of Petroleum Engineers Inc., 2004.9: 26-29
[13] Faru.K.Civan. Characterization and prediction of formation damage in two-phase flow systems[R]. SPE: 25429, 1993
[14] Martin A.Fern, Arne Graue. Numerical Sinulation and Sensityvity Analysis of In-Situ Fluid Flow in MRI Laboratory Waterfloods of Fracturde Carbonate Rocks at Different Wettabilities. SPE: 116145, 2008
[15] Matthew D.Jackson, Per H.Valvatne, Martin J.Blunt. Prediction of wettability variation and its impact on waterflooding using pore-to reservoir-scale simulation[J]. SPE: 77543, 2002
[16] P.Behrenbruch, H.M.Goda. Two-Phase Relative Permeability Prediction:A Comparison of the Modified Brooks-Corey Methodology with a New Carman-Kozeny Based Flow Formulation[J]. SPE: 101150, 2006
[17] Zh.Hongxing, Liu Qingnan, Li Fengqin. Variations of Petrophysical Parameters After Sandstone Reservoirs Watered Out in Daqing Oil Field[J]. SPE: 30884, 1995
[18] LIU Xinghui, FARU K Civan.Characterization and prediction of formation damage in Two-phase flow systems[R]. SPE: 25429, 1993.
[19] Helge Holden, Kenneth Hvistendahl Karlsen, Knut Andreas Lie. Operator splitting methods for degenerate convection-diffusion equations II:Numerical examples with emphasis on reservoir simulation and sedimentation. Canadian Mathematical Society, Conference proceedings
[20]鞠斌山.油藏渗流系统物性变化机理与数学模型研究[D].北京:中国地质大学, 2006
[21]徐守余.渤海湾盆地胜坨油田二区储层微观渗流场演化研究[J].石油实验地质, 2003,20(1):62-65
[22]林光荣,陈付星,邵创国,等.马岭油田长期注水对油层孔隙结构的影响[J].西安石油学院学报(自然科学版),2001,24(2):31-35
[23]单华生,姚光庆,周锋德.储层水先后结构变化规律研究[J].海洋石油,2004,24 (1):62-66
[24]李兴国,孙春红,赵云飞.喇嘛甸油田储层微观结构特征变化规律研究[J].大庆石油地质与开发,2007,26(1):79-82
[25]郭莉,王延斌,刘伟新,等.大港油田注水开发过程中油藏参数变化规律分析[J],石油实验地质,2006,l28(1):85-90
[26]邓玉珍,吴素英.注水开发过程中储层物理特征变化规律研究[J].油气采收率技术, 1996,3(4):44-52
[27]张伟峰,刘守军,李拥安,等.孤岛馆陶组注水开发储层性质动态变化特征研究[J].地球科学与环境学报,2004,26(2):51-54
[28]刘志宏.胜二区沙1-2开发期储层变化特征及特高含水期改善水驱效果研究[D].成都:西南石油学院,2002
[29]龙雨丰,张晓晖,刘鸿,等.低渗油藏注水开发中储层参数变化研究与应用[J].江汉石油学院学报,2001,23(增):86-87
[30]杨晓蓓,冯毅,郭恩常,等.砂砾岩油藏特高含水期储层参数变化规律[J].西南石油学院学报,2000,22(4):22-25
[31]吴素英.胜利油田二区沙二段83层储层润湿性变化规律及对开发效果的影响[J].油气地质与采收率,2006,13(2):71-74
[32]孙焕泉,屈岩涛,房会春,等.胜利油区砂岩储集层敏感性特征研究[J].石油勘探与开发,2000,27(5):72-75
[33]李阳.陆相断陷湖盆油藏流场宏观参数变化规律及动态模型[J].石油学报,2005, 26(2):65-68
[34]陈亮,吴胜和,刘宇红.胡状集油田胡十二块注水开发过程中储层动态变化研究[J].1999,22(4):78-81
[35]李存贵,徐守余.长期注水开发油藏的孔隙结构变化规律[J].石油勘探与开发.2003, 30(2):94-96
[36] Aziz K, Settari A. Petroleum Reservoir Simulation[M]. Applied Science Publishers, 1979
[37]程林松,李丽,要春兰.考虑水驱油藏润湿性变化的数值模拟方法[J].水动力学研究与进展,2003,18(6):86-97
[38]高博禹,彭仕宓,黄述旺.胜坨油田二区沙二段3砂层组分阶段油藏数值模拟[J].石油勘探与开发,2004,31(6):82-84
[39]盖英杰,吕德灵,郭元灵.高含水期油藏分阶段数值模拟[J].油气采收率技术,2000, 7(1):54-56
[40]姜瑞忠,陈月明,邓玉珍,等.胜坨油田二区油藏物理特征参数变化数值模拟研究[J].油气采收率技术,1996,3(2):50-56
[41]崔传智,赵晓燕.考虑油藏储层参数变化的数值模拟研究.水动力学研究与进展, 2004,19(A):912-915
[42]高博禹,彭仕宓,黄述旺,等.高含水期油藏精细数值模拟研究[J].石油大学学报(自然科学版),2005,29(2):11-15
[43]谷建伟.高含水开发期基于微观渗流机理的宏观油藏数值模拟研究[J].中国海洋大学:2004
[44]刘红卫,刘义坤,潘晓军.萨中开发区考虑油藏润湿性影响数值模拟[J].内蒙古石油化工,2008,9:221-222
[45]张洪波,王志章,戴胜群,等.水驱前后油藏参数变化机理研究[J].石油天然气学报(江汉石油学院学报).2005,27(A):665-666
[46]王志章,蔡毅,杨蕾.开发中后期油藏参数变化规律及变化机理[M] .北京:石油工业出版社,1999:50-60
[47]姜瑞忠.储层特征参数变化对油藏开发效果的影响[J].油气田地面工程,2005, 24(4):32-33
[48]秦积舜,李爱芬.油层物理学[M].山东:石油大学出版社,2001:94-178
[49]吴素英,孙国,程会明,等.长期水驱砂岩油藏储层参数变化机理研究[J].油气地质与采收率,2004,11(2):9-11
[50]张大为,闫宝珍,曾照英.注水开发油藏油水相对渗透率变化规律[J].重庆大学学报(自然科学版),2000,23(增):77-78
[51]刘中云,曾庆辉,唐周怀,等.润湿性对采收率及相对渗透率的影响[J].2000,21 (2):940-942
[52]崔浩哲,姚光庆,周锋德.低渗透砂砾岩油层相对渗透率曲线的形态及其变化特征[J].地质科学情报,2003,22(1):87-91
[53]鞠斌山,樊太亮,张金川,等.水驱开发油藏原油黏度变化规律及其对开发效果的影响[J].石油勘探与开发,2006,33(1):99-102
[54]鞠斌山,王春田,李师涛,等.变原油黏度油水两相流驱替特征研究[J].西安石油学院(自然科学版),2005,18(1):17-20
[55]邓瑞健.濮城油田注入水体积流量试验特征[J].江汉石油学院学报,2003,25(3): 111-112
[56]孙士孝,韩锦文,郭云尧,等.胜坨油田岩石水洗前后物性变化的试验研究.1996, 20:33-35
[57]刘月田,徐明旺,彭道贵,等.各向异性渗透率油藏数值模拟[J].计算物理,2007, 24(3):295-300
[58]刘慧卿.油藏数值模拟方法专题[M].山东:石油大学出版社,2001:79-84
[59]皮斯曼D.W.油藏数值模拟基础[M].北京:石油工业出版社,1982:1-20
[60]韩大匡,陈钦雷,闫存章.油藏数值模拟基础[M].北京:石油工业出版社,1999: 30-40
[61]关治,陆金甫,等.数值分析基础[M].北京:高等教育出版社,2002:310-323
[62]李庆阳,关治,等.数值计算原理[M].北京:清华大学出版社,2000:170-236
[63]封建湖.数值分析原理[M].北京:科学技术出版社,2001:173~175
[64]李庆阳,关治.数值计算原理[M].北京:清华大学出版社,2000:171~226
[65]李维国,黄炳家,同登科,等.数值计算方法[M].山东东营:石油大学出版社, 2004:173~180