CN114252381B - 一种裂缝性储层水平井钻井液固相污染后污染程度评价方法 - Google Patents

一种裂缝性储层水平井钻井液固相污染后污染程度评价方法 Download PDF

Info

Publication number
CN114252381B
CN114252381B CN202111571801.XA CN202111571801A CN114252381B CN 114252381 B CN114252381 B CN 114252381B CN 202111571801 A CN202111571801 A CN 202111571801A CN 114252381 B CN114252381 B CN 114252381B
Authority
CN
China
Prior art keywords
reservoir
pollution
calculating
permeability
fluid
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.)
Active
Application number
CN202111571801.XA
Other languages
English (en)
Other versions
CN114252381A (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.)
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 CN202111571801.XA priority Critical patent/CN114252381B/zh
Publication of CN114252381A publication Critical patent/CN114252381A/zh
Application granted granted Critical
Publication of CN114252381B publication Critical patent/CN114252381B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明公开了一种裂缝性储层水平井钻井液固相污染后污染程度评价方法,主要通过考虑钻井过程中天然裂缝在井底压力及泥浆固相堵塞影响下的宽度和流动能力变化以及泥浆固相对储层基质渗透率的伤害,最终评价钻井过程中泥浆固相对储层的伤害程度并获取储层渗透率分布;然后在储层模型中设置固定出口模拟生产井位置,设定出口处流体压力,计算稳态条件下储层内流体压力分布,计算过程中根据天然裂缝内流体压力值修正天然裂缝宽度及渗透率,并迭代直至收敛;然后,分别计算污染条件下产能指数和未污染条件下产能指数,最终计算出储层污染程度。本发明采用了一种相对直接、简便的方法准确评价水平井钻井后泥浆污染程度。

Description

一种裂缝性储层水平井钻井液固相污染后污染程度评价方法
技术领域
本发明涉及石油工程领域,尤其是一种裂缝性储层水平井钻井液固相污染后污染程度评价方法。
背景技术
裂缝性碳酸盐岩气藏是国内当前和未来天然气开发的主阵地,酸化是实现该类油气藏高效开发的主要手段之一。酸化的主要目标是解除由于钻井过程中钻井液固相对储层的伤害,因此准确评价钻井液固相对储层、天然裂缝的污染程度,是支撑酸化参数优化设计、油气井获得高产的基础。
目前国内外针对钻井泥浆污染的评价方法一般分为三类,分别为测井解释类、实验评价类与数值模拟类。基于测井资料可估算钻井液侵入深度,但无法建立天然裂缝参数与钻井液固相侵入特征的关系。实验可评价钻井液对岩心、裂缝的污染特征,但实验结果可表征的尺度较小,在现场设计中的应用受限。数值模拟是描述此类储层最为适用的方法,但现有的研究多针对钻井液侵入天然裂缝的深度开展研究(李松等,基于钻井液漏失侵入深度预测的裂缝性碳酸盐岩储层改造优化.钻采工艺,2018,41(2):(42-45),使得在后续的酸化设计过程中仅能基于泥浆在天然裂缝中的侵入深度开展设计,而无法有效的解除岩石内部的钻井液固相污染。
发明内容
本发明的目的是提供一种相对直接、简便的方法准确评价裂缝性储层水平井钻井后泥浆污染程度。
本发明提供的裂缝性储层水平井钻井液固相污染后污染程度评价方法,步骤如下:
S1、建立天然裂缝钻井液污染模型,赋值储层孔隙度、渗透率、天然裂缝宽度、天然裂缝密度等特征参数。
S2、时间步开始,计算基质及天然裂缝内流体压力及x、y方向的流速,具体如下;
根据公式(1)计算流体压力:
根据公式(2)计算天然裂缝网格渗透率:
根据压力通过立方定律计算流体流速:
式中,p为单元格内流体压力,Pa;kx、ky为单元格x、y方向的平均渗透率,m2;u、v为x、y方向的流速,m/s;μ为流体粘度;φ为单元格内平均孔隙度,无因次;cl为流体压缩系数,Pa-1;ρ为流体密度,kg/m3;knf为天然裂缝渗透率,m2;wnf为天然裂缝宽度,m。
其中,为模拟水平井动态钻进过程,边界条件在水平井位置处设置为定压边界,入口边界范围从模型左边入口开始以等于钻速的速度向模型右侧移动,入口处压力则设置为钻井设计井底压力。
S3、基于流体流速计算结果计算泥浆固相浓度分布,公式如下:
式中,C为由钻井液固相浓度,mol/m3;De为泥浆固相扩散系数,m2/s;t为时间,x、y分别是单元格的x、y方向。
其中,为模拟水平井动态钻进过程,边界条件在水平井位置处设置为定浓度边界,入口边界范围从模型左边入口开始以等于钻速的速度向模型右侧移动。
S4、根据天然裂缝向壁面的滤失速度及天然裂缝网格内钻井液固相浓度计算本时间步内由天然裂缝内向壁面滤失的钻井液固相质量MGX
MGX=CnfvnfΔyhΔt (6)
式中,Cnf为天然裂缝内泥浆固相的质量浓度,kg/m3,由公式(5)得出;Δy为y方向单元格长度,m;h为储层厚度,m;vnf为天然裂缝向壁面的滤失速度,由公式(1)和(4)计算得出;Δt为时间步长度,s。
S5、根据钻井液粒径图和储层孔喉直径分布图,计算进入储层钻井液固相比例fGX,公式如下:
式中,rGX(D)为不同粒径钻井液固相占比;Dm为储层基质孔喉最大直径,m;DGX为钻井液固相粒径,m。
S6、对于不能滤失进入储层孔喉的固相,附着于天然裂缝壁面形成泥饼,计算附着于天然裂缝壁面的泥饼厚度hGX
式中,h为储层厚度,m;Δy为y方向单元格长度,m;ρGX为钻井液固相密度,kg/m3;φGX为泥饼孔隙度,实验测得。
S7、结合泥饼渗透率、泥饼厚度及基质渗透率,计算污染后天然裂缝壁面网格平均渗透率kwr,c':
式中,km为储层基质原始渗透率、Δx为x方向单元格长度,m;kGX'为滤饼渗透率,实验测得。
S8、计算泥饼生成后天然裂缝宽度w'nf
式中,wnf为天然裂缝宽度;hGX,L和hGX,R分别为附着于天然裂缝左、右壁面的泥饼厚度;pt和pt-1分别是本时间步和上时间步的流体压力;Kn为裂缝刚强系数,MP/m。
S9、更新天然裂缝宽度、天然裂缝壁面网格平均渗透率,重复步骤S2~S8。具体是:将步骤S7计算的kwr,c'更新至公式(1)中,将kx、ky均替换为kwr,c',将步骤S8计算的w'nf更新至公式(2)中,重复步骤S2~S8,直至泥浆污染过程结束,获得污染后裂缝性储层渗透率分布特征。
S10、在储层模型中设置固定出口模拟生产井位置,设定出口处流体压力为0.1×106Pa,计算稳态条件下储层内流体压力分布:
其中,计算过程中根据天然裂缝内流体压力值修正天然裂缝宽度及渗透率,并迭代直至收敛,收敛判据取1~103Pa:
计算稳态条件下出口处的总流量Q:
式中,k为生产边界上每个网格的渗透率,m2;h为储层厚度,m;Δx为x方向单元格长度,m;
计算污染条件下产能指数PI:
式中,为模型平均压力;pwf为井底流压;
S11、将储层及天然裂缝渗透率恢复为步骤S1的初始值,重复步骤S10,其中将公式(11)替换为公式(15):
计算未污染条件下产能指数PIini
式中,为模型平均压力;pwf,ini为井底流压。
S12、计算储层污染程度:
与现有技术相比,本发明的有益之处在于:
本发明采用直接、简便的方法来准确评价水平井钻井后泥浆污染程度,重点考虑钻井过程中天然裂缝在井底压力及泥浆固相堵塞影响下的宽度和流动能力变化以及泥浆固相对储层基质渗透率的伤害,最终评价钻井过程中泥浆固相对储层的伤害程度并获取储层渗透率分布,为酸化优化设计提供支撑。
本发明的其它优点、目标和特征将部分通过下面的说明体现,部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。
附图说明
图1、天然裂缝钻井液污染模型示意图。
图2、天然裂缝钻井液污染模型分析示意图。
图3、储层流体压力分布图。
图4、y方向流动速度分布图。
图5、泥浆固相浓度分布图。
图6、由天然裂缝内向壁面滤失的钻井液固相质量分布图。
图7、附着于天然裂缝壁面的泥饼厚度分布图。
图8、天然裂缝壁面网格污染后y方向渗透率分别图。
图9、天然裂缝宽度变化量。
图10、污染后井筒剖面流体流速分布图。
图11、未污染井筒剖面流体流速分布图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
一种裂缝性储层水平井钻井液固相污染后污染程度评价方法,具体步骤如下:
S1、建立天然裂缝钻井液污染模型,如图1和2所示,模型尺寸300m×50m。赋值储层孔隙度0.03、渗透率0.1×10-15m2、天然裂缝宽度0.0001m、天然裂缝密度0.2条/m。
S2、时间步开始,计算基质及天然裂缝内流体压力及x、y方向的流速。得到如图3所示的储层流体压力分布图。如图4所示的y方向流动速度分布图。
S3、基于流速计算结果计算泥浆固相浓度分布,结果如图5所示。
S4、根据天然裂缝向壁面的滤失速度及天然裂缝网格内钻井液固相浓度计算由天然裂缝内向壁面滤失的钻井液固相质量MGX。计算结果如图6所示。
S5、根据钻井液粒径图和储层孔喉直径分布图,确定能够进入储层钻井液固相比例fGX=0.1。
S6、对于不能滤失进入储层孔喉的固相,附着于天然裂缝壁面形成泥饼,首先计算附着于天然裂缝壁面的泥饼厚度。计算结果见图7。
S7、结合泥饼渗透率、泥饼厚度及基质渗透率计算污染后天然裂缝壁面网格平均渗透率。如图8所示,是天然裂缝壁面网格污染后y方向渗透率计算结果图。
S8、计算泥饼生成后天然裂缝宽度,结果见图9。
S9、更新天然裂缝宽度、天然裂缝壁面网格平均渗透率,即将步骤S7计算的kwr,c'更新至公式(1)中,将kx、ky均替换为kwr,c',重复步骤S2~S8,直至泥浆污染过程结束,即可获得污染后裂缝性储层渗透率分布特征。
S10、在储层模型中设置固定出口模拟生产井位置,设定出口处流体压力为0.1×106Pa,污染后井筒剖面流速分布如图10所示。计算污染条件下的产能指数PI为1.8840×10-12
S11、将储层及天然裂缝渗透率恢复为初始值渗透率0.1×10-15m2,重复步骤S10,其中将公式(11)替换为公式(15):
未污染井筒剖面流速分布如图11所示。计算未污染条件下产能指数PIini为9.9690×10-10
S12、计算储层污染程度
计算得到储层污染程度高达99.8%。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。

Claims (4)

1.一种裂缝性储层水平井钻井液固相污染后污染程度评价方法,其特征在于,步骤如下:
S1、建立天然裂缝钻井液污染模型,赋值储层孔隙度、渗透率、天然裂缝宽度、天然裂缝密度这些特征参数;
S2、时间步开始,计算基质及天然裂缝内流体流速,具体如下;
根据公式(1)计算流体压力:
根据公式(2)计算天然裂缝网格渗透率:
根据公式(3)和(4)计算流体流速:
式中,p为单元格内流体压力,Pa;kx、ky为单元格x、y方向的平均渗透率,m2;u、v为x、y方向的流速,m/s;μ为流体粘度;φ为单元格内平均孔隙度,无因次;cl为流体压缩系数,Pa-1;ρ为流体密度,kg/m3;knf为天然裂缝渗透率,m2;wnf为天然裂缝宽度,m;
S3、基于流体流速计算结果计算泥浆固相浓度分布;
S4、计算本时间步内由天然裂缝内向壁面滤失的钻井液固相质量MGX
MGX=CnfvnfΔyhΔt (6)
式中,Cnf为天然裂缝内泥浆固相的质量浓度,kg/m3;Δy为y方向单元格长度,m;h为储层厚度,m;vnf为天然裂缝向壁面的滤失速度,由公式(1)和(4)计算得出;Δt为时间步长度,s;
S5、计算进入储层钻井液固相比例fGX
S6、计算附着于天然裂缝壁面的泥饼厚度hGX
式中,h为储层厚度,m;Δy为y方向单元格长度,m;ρGX为钻井液固相密度,kg/m3;、φGX为泥饼孔隙度,实验测得;
S7、计算污染后天然裂缝壁面网格平均渗透率kwr,c':
式中,km为储层基质原始渗透率、Δx为x方向单元格长度,m;kGX'为滤饼渗透率,实验测得;
S8、计算泥饼生成后天然裂缝宽度w'nf
式中,wnf为天然裂缝宽度;hGX,L和hGX,R分别为附着于天然裂缝左、右壁面的泥饼厚度;pt和pt-1分别是本时间步和上时间步的流体压力;Kn为裂缝刚强系数,MP/m;
S9、更新天然裂缝宽度、天然裂缝壁面网格平均渗透率,重复步骤S2~S8,直至泥浆污染过程结束,获得污染后裂缝性储层渗透率分布特征;
S10、在储层模型中设置固定出口模拟生产井位置,设定出口处流体压力为0.1×106Pa,计算稳态条件下储层内流体压力分布:
其中,计算过程中根据天然裂缝内流体压力值修正天然裂缝宽度及渗透率,并迭代直至收敛,收敛判据取1~103Pa:
计算稳态条件下出口处的总流量Q:
式中,k为生产边界上每个网格的渗透率;
计算污染条件下产能指数PI:
式中,为模型平均压力;pwf为井底流压;
S11、将储层及天然裂缝渗透率恢复为步骤S1的初始值,重复步骤S10,其中将公式(11)替换为公式(15):
计算未污染条件下产能指数PIini
式中,为模型平均压力;pwf,ini为井底流压;
S12、计算储层污染程度:
2.如权利要求1所述的裂缝性储层水平井钻井液固相污染后污染程度评价方法,其特征在于,所述步骤S3中,计算泥浆固相浓度分布公式如下:
式中,C为由钻井液固相浓度,mol/m3;De为泥浆固相扩散系数,m2/s。
3.如权利要求1所述的裂缝性储层水平井钻井液固相污染后污染程度评价方法,其特征在于,所述步骤S5中,计算进入储层钻井液固相比例fGX的公式如下:
式中,rGX(D)为不同粒径钻井液固相占比;Dm为储层基质孔喉最大直径,m;DGX为钻井液固相粒径,m。
4.如权利要求1所述的裂缝性储层水平井钻井液固相污染后污染程度评价方法,其特征在于,所述步骤S9具体是:将步骤S7计算的kwr,c'更新至公式(1)中,将kx、ky均替换为kwr,c',将步骤S8计算的w'nf更新至公式(2)中,重复步骤S2~S8,直至泥浆污染过程结束,获得污染后裂缝性储层渗透率分布特征。
CN202111571801.XA 2021-12-21 2021-12-21 一种裂缝性储层水平井钻井液固相污染后污染程度评价方法 Active CN114252381B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111571801.XA CN114252381B (zh) 2021-12-21 2021-12-21 一种裂缝性储层水平井钻井液固相污染后污染程度评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111571801.XA CN114252381B (zh) 2021-12-21 2021-12-21 一种裂缝性储层水平井钻井液固相污染后污染程度评价方法

Publications (2)

Publication Number Publication Date
CN114252381A CN114252381A (zh) 2022-03-29
CN114252381B true CN114252381B (zh) 2023-09-22

Family

ID=80796317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111571801.XA Active CN114252381B (zh) 2021-12-21 2021-12-21 一种裂缝性储层水平井钻井液固相污染后污染程度评价方法

Country Status (1)

Country Link
CN (1) CN114252381B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116402189A (zh) * 2023-03-02 2023-07-07 西南石油大学 一种深层裂缝性致密储层钻完井液漏失损害预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0733780B1 (en) * 1995-03-23 2000-07-12 Schlumberger Limited (a Netherland Antilles corp.) Method and apparatus for locally measuring flow parameters of a multiphase fluid
CN111219175A (zh) * 2020-01-08 2020-06-02 西南石油大学 考虑应力敏感的裂缝性碳酸盐岩酸压裂缝匹配性优化方法
CN113685163A (zh) * 2021-08-26 2021-11-23 杜海峰 微生物菌剂油井压裂前置液新工艺

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020023754A1 (en) * 2000-08-28 2002-02-28 Buytaert Jean P. Method for drilling multilateral wells and related device
US6712144B2 (en) * 2000-08-28 2004-03-30 Frank's International, Inc. Method for drilling multilateral wells with reduced under-reaming and related device

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0733780B1 (en) * 1995-03-23 2000-07-12 Schlumberger Limited (a Netherland Antilles corp.) Method and apparatus for locally measuring flow parameters of a multiphase fluid
CN111219175A (zh) * 2020-01-08 2020-06-02 西南石油大学 考虑应力敏感的裂缝性碳酸盐岩酸压裂缝匹配性优化方法
CN113685163A (zh) * 2021-08-26 2021-11-23 杜海峰 微生物菌剂油井压裂前置液新工艺

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
页岩气储层钻井液损害评价研究;李红梅;《探矿工程(岩土钻掘工程)》;第43卷(第7期);全文 *

Also Published As

Publication number Publication date
CN114252381A (zh) 2022-03-29

Similar Documents

Publication Publication Date Title
CN109594968B (zh) 一种页岩气多段压裂水平井压后裂缝参数评价方法及系统
CN111197476B (zh) 一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法
CN104695950B (zh) 火山岩油藏产能预测方法
CN109209333B (zh) 页岩气多井群高效开采间距优化方法
CN110826142B (zh) 一种裂缝性地层封堵承压能力的预测方法
CN105631184B (zh) 一种油气井产量递减分析方法和系统
CN105569646A (zh) 一种油气井技术可采储量预测方法
CN108661631A (zh) 一种页岩气井产量预测方法
CN114252381B (zh) 一种裂缝性储层水平井钻井液固相污染后污染程度评价方法
CN113034003B (zh) 一种页岩气井产能快速评价方法
CN107313759A (zh) 低渗稠油油藏直井热采压力分布预测方法及系统
CN108959767A (zh) 一种窄河道型气藏不同井型凝析油伤害数值模拟方法
CN104847314A (zh) 高温高压油气直井单相流射孔完井参数优化方法
CN114818533A (zh) 基于排采数据的页岩油气藏裂缝参数确定方法及装置
CN114427432A (zh) 一种气藏剩余气开发潜力确定方法
CN116335654B (zh) 一种模拟页岩气特殊机理的压裂水平井产量预测方法
CN115828785A (zh) 一种深部煤层气藏嵌套介质气水两相生产动态预测方法
CN115713163A (zh) 深层页岩气压裂水平井产量及可采储量预测方法
CN115270662A (zh) 一种欠饱和储层煤层气井生产边界定量预测的计算方法
CN111946321B (zh) 填砂暂堵压裂的支撑剂参数设计方法
CN114239365B (zh) 一种天然裂缝钻井液污染后壁面等效渗透率的计算方法
CN108664678A (zh) 一种页岩气井产量预测方法
CN109657299B (zh) 页岩气藏开采方法
CN115906677A (zh) 一种非均质性边水气藏水体活跃程度预测方法
CN109812237B (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