CN112065379A - 一种注水井结垢伤害预测方法 - Google Patents

一种注水井结垢伤害预测方法 Download PDF

Info

Publication number
CN112065379A
CN112065379A CN202011083371.2A CN202011083371A CN112065379A CN 112065379 A CN112065379 A CN 112065379A CN 202011083371 A CN202011083371 A CN 202011083371A CN 112065379 A CN112065379 A CN 112065379A
Authority
CN
China
Prior art keywords
scaling
calculating
formation
radius
scale
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.)
Pending
Application number
CN202011083371.2A
Other languages
English (en)
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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202011083371.2A priority Critical patent/CN112065379A/zh
Publication of CN112065379A publication Critical patent/CN112065379A/zh
Pending legal-status Critical Current

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
    • 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
    • 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
    • E21B43/20Displacing by water
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Marine Sciences & Fisheries (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Animal Husbandry (AREA)
  • Agronomy & Crop Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明涉及一种用于油气田注水过程中一种注水井结垢伤害预测方法。它基于开展在地层条件下的注水岩心伤害评价实验拟合的结垢量与井径关系曲线,解决了现今注水井结垢伤害预测中对结垢量模拟准确度不够且计算复杂的问题。其技术方案是:首先开展地层条件下的注入水岩心结垢伤害实验,拟合出结垢量与结垢半径关系式,再通过把地层沿径向划分成若干段来计算单段半径范围内的结垢厚度,再得到结垢后地层孔喉大小;同时,根据结垢体积算出结垢后孔隙度大小,最后根据毛细管模型公式算出地层结垢后渗透率。本方法依托实验拟合得到地层垢分布,再计算孔渗变化,能有效预测、直观了解地层结垢伤害情况。

Description

一种注水井结垢伤害预测方法
技术领域
本发明涉及一种注水井结垢伤害预测方法,属于油气田开发领域。
背景技术
我国新疆、渤海、南海等地区大量油田正处在注水开发阶段。在注水过程中,当注入水流入地层后与地层水混合后,由于其离子成分的差异,易发生沉淀反应,即结垢,从而堵塞地层中的渗流通道,使储层孔隙度、渗透率明显降低,增大表皮系数,严重影响油田生产。在众多垢物中,以碳酸钙最为常见,因此,正确碳酸钙的结垢伤害规律、预测结垢后的伤害程度,进而采取相应的防治或解堵措施变得愈发重要。
预测注水井结垢方法较多,常用的标准方法是SY/T0600-2009《油田水结垢趋势预测》,标准提供的公式中引入的与离子强度和温度有关的修正系数K需要内插法查表获取,过程复杂且易造成人为误差。专利号为CN201710299676.9的《一种用于碳酸钙结垢趋势的预测方法》对油田标准方法做了改进,通过对水样离子成分的的回归分析,得到修正过后的数学预测公式,免去了人为查图版的误差,精度上有一定提升,但需拟合较多数据,且不能判断伤害程度;除此之外,还利用综合了热力学、动力学和流体力学等方面建立的软件如Scale Chem、SDCQPC等来进行预测,但都局限于简单判断是否结垢;目前,专利号为201911170132.8的《一种评价注水井结垢趋势的方法》结合注采井网压力场、温度场、孔渗分布模型建立起的结垢预测模型是较为全面的预测方法,但须考虑注采井网分布,计算量较大且复杂。
发明内容
本发明目的是:为了解决现今注水井结垢后伤害程度预测方法较少,部分方法计算复杂,推广性不强等问题,本发明依托岩心实验数据拟合得到的结垢分布规律公式,结合地层毛细管模型,通过对地层径向单元结垢后引起的毛管半径减小量的计算得到结垢后的地层孔隙度,最后算出地层结垢后的渗透率,本方法计算简单,效果较好,可推广性强。
为实现上述目的,本发明提供了一种注水井结垢伤害预测方法,该方法包括下列步骤:
第一,开展地层条件下的注入水岩心结垢伤害实验,拟合出结垢量与结垢半径关系式;
第二,计算地层毛细管表面积;
第三,应用毛细管模型公式计算生成垢的体积,厚度;
第四,计算结垢后的地层毛管半径,孔隙度;
第五,再应用毛细管模型公式计算渗透率。
上述一种注水井结垢伤害预测方法中,所述毛细管模型公式为
Figure BDA0002719452830000021
式中,K为地层渗透率,单位为D;
Figure BDA0002719452830000022
为地层孔隙度;R为地层毛细管半径,单位为10-3μm。
上述一种注水井结垢伤害预测方法中,所述结垢体积公式为
Figure BDA0002719452830000023
式中,Vd为生成垢的体积,单位为m3;V为地层孔隙体积,单位为m3;y为垢质量浓度,单位为g/m3;ρ为垢密度,单位为g/m3
上述一种注水井结垢伤害预测方法中,所述毛细管表面积计算公式为
Figure BDA0002719452830000024
式中,S为结垢区域毛细管表面积,单位为cm2/cm3;r为结垢区域半径,单位为m;H为储层厚度。
上述一种注水井结垢伤害预测方法中,所述计算垢厚度方法为
h=Vd/S
上述一种注水井结垢地层渗透率预测方法中,从结垢量到结垢后地层渗透率的具体计算步骤为,
第一,将地层沿径向按一定步长划分为若干段,根据岩心伤害实验拟合出的方程,计算出该步长的平均结垢浓度,公式为
Figure BDA0002719452830000025
式中,y为结垢浓度,单位为mg/L;r2与r1分别为计算步长两端点径向半径,单位为m;
第二,计算该步长内孔隙体积,公式为
Figure BDA0002719452830000026
式中,H为储层厚度,单位为m;
第三,计算成垢体积,公式为
Figure BDA0002719452830000027
式中,ρ为垢密度,单位为g/cm3
第四,计算结垢后地层孔隙度,计算方法如下
Figure BDA0002719452830000031
第五,计算成垢后毛管半径;
第六,根据毛管模型公式计算出结垢后地层渗透率;
第七,根据离子浓度变化公式计算出下一步长结垢离子浓度然后重复以上计算过程。
与现有技术相比,本发明具有以下有益效果:(1)计算过程简单;(2)预测快;(3)可推广性强。
附图说明
在附图中:
图1是本方法的技术路线图。
图2是本方法的结垢量与井径的拟合图。
图3是本方法的预测结垢后孔隙度分布图。
图4是本方法的预测结垢后渗透率分布图。
具体实施方式
下面结合实施方式和附图对本发明做进一步说明。
本发明提供了一种页岩气井返排率预测方法,图1为本方法的技术路线图,该方法包括下列步骤:
第一,模拟地层温度与注水速度,开展注入水岩心伤害实验,将伤害后的岩心均匀切分为若干段,并用盐酸滴定伤害后的岩心,过滤并测定Ca2+含量,计算出CaCO3的平均含量,用Excel进行回归分析,得到结垢量与结垢半径之间的关系式;
第二,分段计算各段环形单元内总结垢质量,首先应用实验拟合公式计算的平均结垢量,再计算单元体积内的孔隙体积,最后相乘的到结垢质量;
第三,根据每段得到的结垢质量除以垢密度即得到垢体积,假设垢在孔隙体积中均匀分布,用垢体积除以毛管表面积则得垢厚度,从而得到成垢后毛管孔径,毛管结垢后孔隙度;
第四,应用毛管模型孔径、渗透率、孔隙度的关系式计算出渗透率;
第五,重复以上步骤,直到计算出计算范围等于结垢总半径。
进一步的,所述岩心伤害实验需采用环氧胶结的人造岩心以排除岩心自身Ca2+的影响。
进一步的,实验中注入流体速度的设定需结合现场日配注量大小来确定,计算公式如下,
Figure BDA0002719452830000041
式中,v为流体流速,单位为m/s;Qpeizhu为现场配注量,单位为m3/d;Spaoyan为井筒射孔总表面积,单位为m2
进一步的,将滴定得到的Ca2+质量除以注入时间通过岩心的总液量得到垢浓度分布。
进一步的,各段孔隙体积计算公式应为
Figure BDA0002719452830000042
式中,rn=n*步长+rw,rw为井筒半径单位为m,n表示目前正计算第几段。
进一步的,结垢后,地层毛孔半径大小计算式为
rl=r-h
式中,rl为结垢后地层毛孔半径,单位为m;r为地层原始孔吼半径,单位为m;h为结垢厚度,单位为m。
以某油田地层A1注水井为例,该井配注量为50m3/d,井筒半径为0.1m,注水层段厚度为18m,地下射孔总面积为0.007m2,地层渗透率50mD,孔隙度为20%,地层温度为80℃;
首先制备直径为40mm,长为200mm的环氧胶结人造岩心,测得其空气渗透率为30-60mD,用相应过滤后的地层水饱和岩心,并将岩心夹持器置于80℃的恒温水浴中,以0.008m/s的速度持续注入过滤后的注入水直到在最大压力下注不进为止,将岩心取下并均匀切成五段,粉碎后用盐酸滴定过滤液中的Ca2+含量,算出CaCO3的平均含量,如图2所示,拟合结垢量与结垢半径关系式为100y=40×136.05×r-0.664
进一步的,根据曲线分析,当半径为10m时,结垢量已接近于0,因此可以将10m作为结垢半径。以为伤害半径为10m,设定计算单元半径0.1m,则需要从井筒开始,重复上述计算步骤,共计算100段。
根据上述步骤首先应用实验拟合公式计算的平均结垢量,再计算单元体积内的孔隙体积,然后相乘得到结垢质量;根据每段得到的结垢质量除以垢密度即得到垢体积,假设垢在孔隙体积中均匀分布,用垢体积除以毛管表面积则得垢厚度,从而得到成垢后毛管孔径,最后算出毛管结垢后孔隙度,分段计算得到沿井半径方向的结垢后孔隙度分布如图3所示;接着应用毛管模型孔径、渗透率、孔隙度的关系式计算出渗透率分布,如图4所示;
进一步的,上述计算假设了水垢在单位步长内是均匀分布的,结垢量是基于实验数据拟合,其中结垢量是在注入水压力较大时注不进后测定值,所以一定程度上反映出了当注水井由于结垢导致注水井停注后的最大结垢量。
与现有技术相比,本发明具有以下有益效果:(1)结垢分布依托实验数据,较为准确和真实;(2)需要数据较少;(3)可推广性强。
最后所应说明的是:以上实施例仅用以说明而非限制本发明的技术方案,尽管参照上述实施例对本发明进行了详细说明,本领域的普通技术人员应该理解:依然可以对本发明进行修改或者等同替换,而不脱离本发明的精神和范围的任何修改或局部替换,其均应涵盖在本发明的权利要求范围当中。

Claims (4)

1.一种注水井结垢伤害预测方法,其特征在于,该方法包括下列步骤:
S100、开展地层条件下的注入水岩心结垢伤害实验,拟合出结垢量与结垢半径关系式;
S200、计算地层毛细管表面积;
S300、应用毛细管模型公式计算生成垢的体积,厚度;
S400、计算结垢后的孔吼半径,孔隙度;
S500、再应用毛细管模型公式计算渗透率。
2.根据权利要求1所述的一种注水井结垢伤害预测方法,其特征在于:实验过程采用环氧树脂胶结岩心,考虑影响结垢的两个主要因素温度和压力,根据地层温度及注水井配注量,设定岩心伤害实验的温度和驱替速度,最后根据将岩心均匀分段并用盐酸滴定,得出各段CaCO3生成量,在用Excel拟合出结垢量与井半径分布关系式。
3.根据权利要求1所述的一种注水井结垢伤害预测方法,其特征在于:计算生成垢后的孔吼半径的步骤为,
S401、用下式计算单个径向步长的平均结垢量,
Figure FDA0002719452820000011
式中,y为结垢浓度,单位为mg/L;r2与r1分别为计算步长两端点径向半径,单位为m;
S402、计算该步长孔隙体积,
Figure FDA0002719452820000012
式中,H为储层厚度,单位为m;
Figure FDA0002719452820000013
为地层孔隙度;V为单个径向步长孔隙总体积,单位为m3
S403、计算生成垢体积,
Figure FDA0002719452820000014
式中,Vd为生成垢的体积,单位为m3;V为地层孔隙体积,单位为m3;y为垢质量浓度,单位为g/m3;ρ为垢密度,单位为g/m3
S404、计算单个步长内的毛细管表面积,
Figure FDA0002719452820000021
式中,S为结垢区域毛细管表面积,单位为cm2/cm3;H为储层厚度;
S405、计算生成垢的厚度,
h=Vd/S
S406、计算生成垢后孔喉半径公式为,
Rl=R-h
式中,Rl为结垢后地层毛孔半径,单位为m;R为地层原始孔吼半径,单位为m;h为结垢厚度,单位为m。
4.根据权利要求1所述的一种注水井结垢伤害预测方法,其特征在于:计算渗透率的方法为,
S501、计算出结垢后孔隙度,公式为
Figure FDA0002719452820000022
S502、应用毛管模型公式,带入结垢后孔喉半径,结垢后孔隙度,用下式计算出结垢后地层渗透率,
Figure FDA0002719452820000023
式中,K为地层渗透率,单位为D;
Figure FDA0002719452820000024
为地层孔隙度;R为地层毛细管半径,单位为10-3μm。
CN202011083371.2A 2020-10-12 2020-10-12 一种注水井结垢伤害预测方法 Pending CN112065379A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011083371.2A CN112065379A (zh) 2020-10-12 2020-10-12 一种注水井结垢伤害预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011083371.2A CN112065379A (zh) 2020-10-12 2020-10-12 一种注水井结垢伤害预测方法

Publications (1)

Publication Number Publication Date
CN112065379A true CN112065379A (zh) 2020-12-11

Family

ID=73655103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011083371.2A Pending CN112065379A (zh) 2020-10-12 2020-10-12 一种注水井结垢伤害预测方法

Country Status (1)

Country Link
CN (1) CN112065379A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113177378A (zh) * 2021-05-28 2021-07-27 中煤科工集团杭州研究院有限公司 一种用于预测矿井水在反渗透膜组件内部结垢的分析方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1846040A (zh) * 2003-06-25 2006-10-11 罗地亚化学公司 使用不同的水垢抑制剂使油田增产的方法
CN101993685A (zh) * 2009-08-12 2011-03-30 中国石油天然气股份有限公司 一种复合防垢剂
US20120118572A1 (en) * 2006-09-18 2012-05-17 Schlumberger Technology Corporation Methods of limiting leak off and damage in hydraulic fractures
CN104677775A (zh) * 2015-02-06 2015-06-03 李年银 一种测试含有机碱的螯合剂溶液除钡锶垢性能及对岩心渗透率影响的方法
CN108756817A (zh) * 2018-05-14 2018-11-06 中国石油大学(华东) 判断产水气井井筒结垢风险及确定防垢剂注入时机的方法
US20190176101A1 (en) * 2016-06-10 2019-06-13 OPEC Remediation Technologies Pty Limited Method And Apparatus For Separation Of A Substance From Groundwater
CN111101925A (zh) * 2019-11-26 2020-05-05 西南石油大学 一种评价注水井结垢趋势的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1846040A (zh) * 2003-06-25 2006-10-11 罗地亚化学公司 使用不同的水垢抑制剂使油田增产的方法
US20120118572A1 (en) * 2006-09-18 2012-05-17 Schlumberger Technology Corporation Methods of limiting leak off and damage in hydraulic fractures
CN101993685A (zh) * 2009-08-12 2011-03-30 中国石油天然气股份有限公司 一种复合防垢剂
CN104677775A (zh) * 2015-02-06 2015-06-03 李年银 一种测试含有机碱的螯合剂溶液除钡锶垢性能及对岩心渗透率影响的方法
US20190176101A1 (en) * 2016-06-10 2019-06-13 OPEC Remediation Technologies Pty Limited Method And Apparatus For Separation Of A Substance From Groundwater
CN108756817A (zh) * 2018-05-14 2018-11-06 中国石油大学(华东) 判断产水气井井筒结垢风险及确定防垢剂注入时机的方法
CN111101925A (zh) * 2019-11-26 2020-05-05 西南石油大学 一种评价注水井结垢趋势的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李年银等: "注采单元地层结垢机理及数值模拟研究", 《油气藏评价与开发》 *
田野: "结垢对注水井注入压力影响研究", 《中国优秀硕士学位论文数据库工程科技1辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113177378A (zh) * 2021-05-28 2021-07-27 中煤科工集团杭州研究院有限公司 一种用于预测矿井水在反渗透膜组件内部结垢的分析方法
CN113177378B (zh) * 2021-05-28 2022-05-13 中煤科工集团杭州研究院有限公司 一种用于预测矿井水在反渗透膜组件内部结垢的分析方法

Similar Documents

Publication Publication Date Title
Dayalan et al. CO2 corrosion prediction in pipe flow under FeCO3 scale-forming conditions
CN107045671A (zh) 产水气井积液风险预测方法
CN107437127B (zh) 一种油井停喷地层压力预测方法
CN106644385B (zh) 一种地表水与地下水潜流交换自循环试验装置及使用方法
CN108868731B (zh) 一种裂缝性储层酸压动态综合滤失系数的计算方法
CN105160071B (zh) 一种适合气液同产水平井井下工况的判别方法
CN106869918A (zh) 海上油井产能测试实时调整方法
CN105443093B (zh) 用于注聚井的井口组合测试装置及其方法
CN111502652A (zh) 一种三孔介质气藏水平井产量递减及生产动态预测方法
CN107271328A (zh) 一种聚合物溶液抗剪切性能评价装置及方法
CN112065379A (zh) 一种注水井结垢伤害预测方法
CN108071392A (zh) 一种海上异常高压气藏动态储量计算方法
Marty et al. The effect of flow rate on the in-situ gelation of a chrome/redox/polyacrylamide system
CN109339768B (zh) 一种钻井微溢流随钻监测方法
CN110685677A (zh) 一种非均质油藏水平井控水模拟实验装置及其实验方法
RU2189443C1 (ru) Способ определения характеристик скважины, призабойной зоны и пласта
CN115060875A (zh) 一种基于达西定律的水合物储层生产压力区间的确定方法
CN111287741B (zh) 一种致密油藏体积压裂改造区渗透率的快速计算方法
CN110410044B (zh) 气驱co2、n2开发方式下区块产油量计算方法
CN113177363A (zh) 一种储层大孔道参数定量表征的油藏工程方法
CN116796655A (zh) 一种定量预测注水井筒动态结垢厚度的方法
CN112012728A (zh) 一种聚合物在储层隙间流速的确定方法
CN206300878U (zh) 一种单相稳定渗流装置
CN106404631B (zh) 致密油压裂缝网油水两相导流能力测试系统及测试方法
CN110175379B (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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20201211

WD01 Invention patent application deemed withdrawn after publication