CN110671084A - 普通稠油水溶性降粘剂驱数值模拟的方法 - Google Patents

普通稠油水溶性降粘剂驱数值模拟的方法 Download PDF

Info

Publication number
CN110671084A
CN110671084A CN201910831256.XA CN201910831256A CN110671084A CN 110671084 A CN110671084 A CN 110671084A CN 201910831256 A CN201910831256 A CN 201910831256A CN 110671084 A CN110671084 A CN 110671084A
Authority
CN
China
Prior art keywords
water
viscosity reducer
viscosity
oil
soluble
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910831256.XA
Other languages
English (en)
Other versions
CN110671084B (zh
Inventor
刘祖鹏
吴光焕
邓宏伟
韦涛
徐永辉
李洪毅
李伟
唐亮
曹秋莉
魏超平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
China Petrochemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petrochemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petrochemical Corp, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petrochemical Corp
Priority to CN201910831256.XA priority Critical patent/CN110671084B/zh
Publication of CN110671084A publication Critical patent/CN110671084A/zh
Application granted granted Critical
Publication of CN110671084B publication Critical patent/CN110671084B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16Enhanced recovery methods for obtaining hydrocarbons
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Lubricants (AREA)

Abstract

本发明涉及油田开发技术领域,特别是涉及到一种普通稠油油藏水溶性降粘剂驱数值模拟的方法。其包括:步骤1,收集水溶性降粘剂的基本性能参数和降粘实验数据;步骤2,对水溶性降粘剂降低原油粘度的实验进行趋势线拟合,确定拟合公式;步骤3,利用拟合得到的公式计算不同水溶性降粘剂质量浓度下的油水混合液粘度;步骤4,计算非线性混合法则中不同质量浓度下水溶性降粘剂的权重因子;步骤5,建立水溶性降粘剂驱油藏数值模拟模型;步骤6,调整水溶性降粘剂的权重因子,将水溶性降粘剂驱岩心驱油实验结果与油藏数值模拟的计算结果进行拟合,直到符合工程计算精度。该方法方便快捷、参数要求少、应用简单,实现了水溶性降粘剂驱的数值模拟。

Description

普通稠油水溶性降粘剂驱数值模拟的方法
技术领域
本发明涉及油田开发技术领域,特别是涉及到一种普通稠油油藏水溶性降粘剂驱数值模拟的方法。
背景技术
稠油油藏开发方法主要分为热采法和冷采法。热采方法通过加热储层降低原油粘度,提高原油在地层中的流动性。传统的热采方法包括蒸汽吞吐、蒸汽驱、火烧油层、SAGD、热力复合开发技术等。然而,由于热力采油的生产投资成本高、蒸汽锅炉的污染大、技术要求高等问题,在低油价形势下面临的挑战逐渐增加。
近年来,提高稠油产量的化学方法越来越受到重视。稠油化学驱体系主要包括碱驱、表面活性剂驱、聚合物驱和聚合物-碱、碱-表面活性剂、聚合物-表面活性剂等复合驱体系。然而,稠油化学降粘剂驱的研究不多。
化学降粘剂可分为油溶性降粘剂和水溶性降粘剂。在稠油油藏开发过程中,加入化学降粘剂能够提高近井段的流动性能,提高单井生产效果。目前只有实验室实验研究和少量工业试验取得了良好的采收率,尚未得到广泛应用,主要是由于尚没有建立水溶性降粘剂的油藏数值模拟实现方法,无法得到该技术的油藏工程设计,因此限制了该技术的发展。
发明内容
本发明的目的是提供一种具实用性强、方便快捷、参数要求少的水溶性降粘剂驱数值模拟的方法,利用该方法可以快速将室内水溶性降粘剂降低原油粘度实验数据转化为油藏数值模拟器的数值数据,有利于矿场快速开展油藏工程方案设计。
为实现上述目的,本发明采用以下技术方案:
本发明主要提供一种稠油水溶性降粘剂驱数值模拟的方法,其包括:
步骤1,收集水溶性降粘剂的基本性能参数和降粘实验数据;
步骤2,对水溶性降粘剂降低原油粘度的实验进行趋势线拟合,确定拟合公式;
步骤3,利用拟合得到的公式计算不同水溶性降粘剂质量浓度下的油水混合液粘度;
步骤4,计算非线性混合法则中不同质量浓度下水溶性降粘剂的权重因子,采用以下公式:
Figure BDA0002190624350000011
式中,f(xi)为不同水溶性降粘剂浓度下的权重因子,xi为水溶性降粘剂质量浓度,μi为加入降粘剂后的油水混合液粘度,μo为原油初始粘度,μw为水溶性降粘剂粘度;
步骤5,建立水溶性降粘剂驱油藏数值模拟模型,油藏数值模拟器采用以下公式对原油粘度的变化进行插值计算:
μi'=eb
b=f(xi)×ln(μi)+(1-f(xi))×ln(μw)
式中,μi'为数值模拟器中计算的粘度,e为自然对数,b为系数;
步骤6,调整水溶性降粘剂的权重因子,将水溶性降粘剂驱岩心驱油实验结果与油藏数值模拟的计算结果进行拟合,直到符合工程计算精度。
以上所述方法,优选地,在步骤2中,根据水溶性降粘剂质量浓度与降粘后油水混合液粘度的变化曲线,在油水混合液粘度所在的坐标轴上,做出表征粘度变化趋势特征的乘幂公式:y=A*xn,其中,A表示该乘幂方程的系数,n表示乘幂的次数;该趋势线拟合方程与降粘后油水混合液粘度变化曲线的相关系数为R2,R2越大,表明曲线与相应的乘幂方程的相关性越好,R2越小,表明曲线与相应的乘幂方程的相关性越差。
以上所述方法,优选地,当相关系数R2大于0.95时,认为该拟合方程符合工程计算精度要求,可以用于下一步数值模拟计算。
以上所述方法,优选地,在步骤3中,设定水溶性降粘剂质量浓度变化区间为[a,b],将浓度变化区间划分为10等份,形成a1,a2,…,a11共11个数据点,利用步骤2得到的公式,计算区间内油水混合液粘度μi'。
以上所述方法,优选地,在步骤5中,油藏数值模拟模型中共有油、水两相,原油、水、水溶性降粘剂三个组分,其中水和水溶性降粘剂为水相。
以上所述方法,优选地,在步骤6中,利用油藏数值模拟软件CMG建立岩心驱替数值模型,调整非线性函数的权重因子,开展实验结果与数值计算结果的拟合,直到符合工程计算精度。
与现有技术相比,本发明具有以下优势:
本发明中的普通稠油水溶性降粘剂驱数值模拟的方法,实用性强、方便快捷、参数要求少、应用简单,具有较好的操作性,有效实现了水溶性降粘剂驱的数值模拟。
附图说明
图1为本发明的普通稠油水溶性降粘剂驱数值模拟的方法的一具体实施例的流程图;
图2是本发明的一具体实施例的水溶性降粘剂的权重因子曲线;
图3是本发明的一具体实施例的水溶性降粘剂驱岩心驱油实验结果与数值模拟计算的结果拟合曲线;
图4是本发明的一具体实施例的调整后的水溶性降粘剂的权重因子曲线。
具体实施方式
应该指出,以下详细说明都是示例性的,旨在对本发明提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本发明的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作和/或它们的组合。
为了使得本领域技术人员能够更加清楚地了解本发明的技术方案,以下将结合具体的实施例详细说明本发明的技术方案。
下述实施例中所涉及的仪器、试剂、材料等,若无特别说明,均为现有技术中已有的常规仪器、试剂、材料等,可通过正规商业途径获得。下述实施例中所涉及的实验方法,检测方法等,若无特别说明,均为现有技术中已有的常规实验方法,检测方法等。
实施例
如图1所示,图1为本实施例的普通稠油水溶性降粘剂驱数值模拟的方法的。
步骤101,收集水溶性降粘剂的基本性能参数和降粘实验数据。水溶性降粘剂在常温条件浓度为1.02mPa.s,质量密度为1023kg/m3,分别测试水溶性降粘剂质量浓度为0.1%,0.3%,0.5%,0.8%,1%时,得到的油水混合液的浓度为14.5,8.4,6.5,5.3,4.7mPa.s。
步骤102,在步骤101的基础上,对水溶性降粘剂降低原油粘度的实验进行趋势线拟合,,得到了拟合的乘幂公式:
y=4.73×x0.486
步骤103,在步骤102的基础上,利用拟合得到的公式计算不同水溶性降粘剂质量浓度下的油水混合液粘度。水溶性降粘剂的质量浓度区间为[0,1%],划分的11个等分点为0.0%,0.1%,0.2%,0.3%,0.4%,0.5%,0.6%,0.7%,0.8%,0.9%,1.0%;计算得到的油水混合液粘度为1330.0,14.5,10.3,8.5,7.4,6.6,6.1,5.6,5.3,5.0,4.7mPa.s。
步骤104,在步骤103的基础上,计算非线性混合法则中不同质量浓度下水溶性降粘剂的权重因子:
Figure BDA0002190624350000041
式中,f(xi)为不同水溶性降粘剂浓度下的权重因子,xi为水溶性降粘剂质量浓度,μi为加入降粘剂后的油水混合液粘度,μo为原油初始粘度,μw为水溶性降粘剂粘度。
xi为0.0%,0.1%,0.2%,0.3%,0.4%,0.5%,0.6%,0.7%,0.8%,0.9%,1.0%,μi为1330.0,14.5,10.3,8.5,7.4,6.6,6.1,5.6,5.3,5.0,4.7mPa.s,μo为1330mPa.s,μo为1mPa.s,计算的水溶性降粘剂的权重因子与质量浓度的关系曲线如图2所示。
步骤105,在步骤104的基础上,油藏数值模拟器采用以下公式对原油粘度的变化进行插值计算。
μi'=eb
b=f(xi)×ln(μi)+(1-f(xi))×ln(μw)
式中,μi'为数值模拟器中计算的粘度,e为自然对数,b为系数。
在步骤106中,利用油藏数值模拟软件CMG建立岩心驱替数值模型,调整非线性函数的权重因子,开展实验结果与数值计算结果的拟合,直到符合工程计算精度。室内水溶性降粘剂驱实验结果与数值模拟的含水率和采出程度结果如图3所示,调整后得到的非线性含水的权重因子曲线如图4所示。
对于胜利油区多个区块的水溶性降粘剂驱的室内实验与数值模拟计算结果对比,吻合程度大于90%,满足了开发上对水溶性降粘剂驱数值模拟的需求,对于指导矿场应用具有一定的指导意义。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (6)

1.一种稠油水溶性降粘剂驱数值模拟的方法,其特征在于,其包括:
步骤1,收集水溶性降粘剂的基本性能参数和降粘实验数据;
步骤2,对水溶性降粘剂降低原油粘度的实验进行趋势线拟合,确定拟合公式;
步骤3,利用拟合得到的公式计算不同水溶性降粘剂质量浓度下的油水混合液粘度;
步骤4,计算非线性混合法则中不同质量浓度下水溶性降粘剂的权重因子,采用以下公式:
Figure FDA0002190624340000011
式中,f(xi)为不同水溶性降粘剂浓度下的权重因子,xi为水溶性降粘剂质量浓度,μi为加入降粘剂后的油水混合液粘度,μo为原油初始粘度,μw为水溶性降粘剂粘度;
步骤5,建立水溶性降粘剂驱油藏数值模拟模型,油藏数值模拟器采用以下公式对原油粘度的变化进行插值计算:
μi'=eb
b=f(xi)×ln(μi)+(1-f(xi))×ln(μw)
式中,μi′为数值模拟器中计算的粘度,e为自然对数,b为系数;
步骤6,调整水溶性降粘剂的权重因子,将水溶性降粘剂驱岩心驱油实验结果与油藏数值模拟的计算结果进行拟合,直到符合工程计算精度。
2.根据权利要求1所述的方法,其特征在于,在步骤2中根据水溶性降粘剂质量浓度与降粘后油水混合液粘度的变化曲线,在油水混合液粘度所在的坐标轴上,做出表征粘度变化趋势特征的乘幂公式:y=A*xn,其中,A表示该乘幂方程的系数,n表示乘幂的次数;该趋势线拟合方程与降粘后油水混合液粘度变化曲线的相关系数为R2,R2越大,表明曲线与相应的乘幂方程的相关性越好,R2越小,表明曲线与相应的乘幂方程的相关性越差。
3.根据权利要求2所述的方法,其特征在于,当相关系数R2大于0.95时,认为该拟合方程符合工程计算精度要求,可以用于下一步数值模拟计算。
4.根据权利要求1所述的方法,其特征在于,在步骤3中,设定水溶性降粘剂质量浓度变化区间为[a,b],将浓度变化区间划分为10等份,形成a1,a2,…,a11共11个数据点,利用步骤2得到的公式,计算区间内油水混合液粘度μi′。
5.根据权利要求1所述的方法,其特征在于,在步骤5中,油藏数值模拟模型中共有油、水两相,原油、水、水溶性降粘剂三个组分,其中水和水溶性降粘剂为水相。
6.根据权利要求1所述的方法,其特征在于,在步骤6中,利用油藏数值模拟软件CMG建立岩心驱替数值模型,调整非线性函数的权重因子,开展实验结果与数值计算结果的拟合,直到符合工程计算精度。
CN201910831256.XA 2019-09-04 2019-09-04 普通稠油水溶性降粘剂驱数值模拟的方法 Active CN110671084B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910831256.XA CN110671084B (zh) 2019-09-04 2019-09-04 普通稠油水溶性降粘剂驱数值模拟的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910831256.XA CN110671084B (zh) 2019-09-04 2019-09-04 普通稠油水溶性降粘剂驱数值模拟的方法

Publications (2)

Publication Number Publication Date
CN110671084A true CN110671084A (zh) 2020-01-10
CN110671084B CN110671084B (zh) 2021-07-06

Family

ID=69076001

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910831256.XA Active CN110671084B (zh) 2019-09-04 2019-09-04 普通稠油水溶性降粘剂驱数值模拟的方法

Country Status (1)

Country Link
CN (1) CN110671084B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007076044A2 (en) * 2005-12-22 2007-07-05 Chevron U.S.A. Inc. Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive
CN107729587A (zh) * 2016-08-12 2018-02-23 中国石油化工股份有限公司 考虑降粘剂注入的数值模拟方法
CN107869336A (zh) * 2016-09-23 2018-04-03 中国石油化工股份有限公司 超稠油化学辅助热采提高采收率方法
CN107869338A (zh) * 2017-11-03 2018-04-03 中国石油化工股份有限公司 化学驱普通稠油动态降粘评价方法
CN109214016A (zh) * 2017-06-30 2019-01-15 中国石油化工股份有限公司 一种特稠油油藏化学冷采优化方法
CN109308358A (zh) * 2017-07-26 2019-02-05 中国石油化工股份有限公司 一种稠油油藏化学驱数值模拟方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007076044A2 (en) * 2005-12-22 2007-07-05 Chevron U.S.A. Inc. Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive
CN107729587A (zh) * 2016-08-12 2018-02-23 中国石油化工股份有限公司 考虑降粘剂注入的数值模拟方法
CN107869336A (zh) * 2016-09-23 2018-04-03 中国石油化工股份有限公司 超稠油化学辅助热采提高采收率方法
CN109214016A (zh) * 2017-06-30 2019-01-15 中国石油化工股份有限公司 一种特稠油油藏化学冷采优化方法
CN109308358A (zh) * 2017-07-26 2019-02-05 中国石油化工股份有限公司 一种稠油油藏化学驱数值模拟方法
CN107869338A (zh) * 2017-11-03 2018-04-03 中国石油化工股份有限公司 化学驱普通稠油动态降粘评价方法

Also Published As

Publication number Publication date
CN110671084B (zh) 2021-07-06

Similar Documents

Publication Publication Date Title
RU2732868C1 (ru) Способ расчета продуктивности горизонтальных скважин в залежах сланцевого газа при гидроразрыве в условиях нестационарной диффузии
Saaltink et al. A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage
CN111222271A (zh) 基于基质-裂缝非稳态窜流油藏裂缝数值模拟方法及系统
CN106337680A (zh) 一种油藏生产井可采储量的确定方法
CN111104766A (zh) 基于离散裂缝模型的油水两相非达西渗流数值模拟方法
Cao et al. Starch graft copolymer and polymer gel applied in Bohai oilfield for water plugging and profile control and their mechanisms
Lin A new approximate iteration solution of Blasius equation
CN110671084B (zh) 普通稠油水溶性降粘剂驱数值模拟的方法
Huang et al. Preparation and experimental study of a low-initial-viscosity gel plugging agent
CN114154430A (zh) 一种压裂油藏co2驱油流动模拟方法
Śliwa et al. Theoretical model of borehole heat exchanger
CN108959687A (zh) 一种页岩气水平井重复压裂的评估方法
Cao et al. Study on viscosity reducer flooding technology for deep low permeability extra heavy oil reservoirs
Ding et al. Experimental investigation of the heterogeneity limit for the application of polymer flooding in reservoirs
Chai et al. Production Characteristics, Evaluation, and Prediction of CO2 Water-Alternating-Gas Flooding in Tight Oil Reservoir
CN112012728A (zh) 一种聚合物在储层隙间流速的确定方法
Han et al. Experimental study of polyacrylamide polymer gel for profile controlling in high-temperature offshore reservoirs
CN115573689A (zh) 稠油油藏降粘压驱开发合理驱替压差的确定方法
CN113027399A (zh) 一种基于微观流动模拟获取高含水区块水驱曲线方法
Wang et al. Comprehensive Effects of Polymer Flooding on Oil–Water Relative Permeabilities
CN110488776A (zh) 一种水务供水系统的随机饱和控制方法
CN114139432A (zh) 利用神经网络技术的裂缝性油藏co2驱流动模拟方法
Alrajhi et al. The relationship between electricity consumption and economic growth in the Kingdom of Saudi Arabia
CN109902354A (zh) 一种基于界面压缩的多相流界面捕获方法
CN108335016B (zh) 一种海上特/超稠油热采采出液处理流程动态模拟方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant