CN113323658A - 一种井筒漏失压力计算方法 - Google Patents

一种井筒漏失压力计算方法 Download PDF

Info

Publication number
CN113323658A
CN113323658A CN202110753840.5A CN202110753840A CN113323658A CN 113323658 A CN113323658 A CN 113323658A CN 202110753840 A CN202110753840 A CN 202110753840A CN 113323658 A CN113323658 A CN 113323658A
Authority
CN
China
Prior art keywords
superposition
formula
pressure
crack
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.)
Granted
Application number
CN202110753840.5A
Other languages
English (en)
Other versions
CN113323658B (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 CN202110753840.5A priority Critical patent/CN113323658B/zh
Publication of CN113323658A publication Critical patent/CN113323658A/zh
Application granted granted Critical
Publication of CN113323658B publication Critical patent/CN113323658B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing 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
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes 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)
  • Geophysics (AREA)
  • Drilling And Exploitation, And Mining Machines And Methods (AREA)

Abstract

本发明公开了一种井筒漏失压力计算方法，属于石油与天然气钻井技术领域。其特征在于：首先收集相关测井和钻井资料，确定待钻井裂缝参数、计算点位置参数和工程许容漏失速率；然后计算钻井液初始密度，从而得到钻井液密度叠加更新后的井底压力；进一步建立裂缝内钻井液流动平面坐标系，从而得到单裂缝漏失速率；进一步计算水平段裂缝总漏失速率，从而判断是否达到工程许容漏失速率，最终确定出井筒计算点的漏失压力。本发明能够快速、准确预测井筒漏失压力，为裂缝性地层钻井提供技术支撑，提高井下作业的安全性。

Description

一种井筒漏失压力计算方法
技术领域
本发明涉及石油与天然气钻井技术领域，具体是涉及一种井筒漏失压力计算方法。
背景技术
随着油气资源勘探开发力度的不断加大，钻井过程中井漏所带来的危害日益突出，尤其对于裂缝性地层而言，安全密度窗口窄，承压能力低，采用过平衡方式钻井时井漏情况严重。其中，井漏是指井筒内液柱压力大于地层流体的压力时，钻井液漏入地层孔隙或裂缝的现象。而为避免井漏事故的发生，关键就在于所计算的漏失压力是否准确，若预测结果与实际地层漏失压力误差较大，将导致钻井液大量漏入地层，造成井塌、井喷、卡钻等井下复杂事故，因此建立一种便于工程实际应用，能快速、准确预测井筒漏失压力的方法是很有必要的。
目前，针对井筒漏失压力计算方法的研究中，“漏失压力确定方法、装置及存储介质”(申请公布号：CN 113027426 A)主要通过对多个井筒钻井井段进行地层承压试验以获取钻井井段不会发生漏失的安全环空压力，确定每个采样点的破裂压力和环空压耗，并结合钻井井段内钻井液的漏失状态，从而确定待钻井井筒的漏失压力，但该方法基于大量实验确定安全环空压力，时间和经济成本较大，且裂缝发育地层实钻表现为漏而不破，因而具有一定局限性。“一种新型窄安全密度窗口地层漏失压力测量的方法”(申请公布号：CN109458171 A)主要通过停泵、关闭节流阀以提高套压值至设计套压值的方法，观察套压变化情况，从而获取稳定套压值，最终计算得到地层漏失压力为液柱压力与稳定套压值之和，但该方法停泵观察时影响正常的钻进作业，增加了经济成本，同时在套压提高过程中，易导致钻井液大量漏失而发生井下复杂事故，不利于保证井下作业的安全。
发明内容
鉴于上述问题，本发明的目的在于提供一种井筒漏失压力计算方法，以解决现有的水平井在裂缝性地层钻进过程中难以准确预测漏失压力的问题，数值计算过程简便，可有效减少井下复杂事故的发生。
本发明采用以下技术方案，一种井筒漏失压力计算方法，其特征在于，包括以下步骤：
步骤一：根据所收集的邻井岩心和成像测井资料，确定待钻井水平段的裂缝参数，包括裂缝宽度ae、裂缝高度be、裂缝线密度ke和裂缝长度le；根据钻井工程设计资料，选取漏失压力计算点并确定位置参数，包括计算点井深L、计算点垂深H和水平段长度ΔL；根据邻井水平段进行堵漏作业时的漏失情况，确定工程许容漏失速率[qt]；
步骤二：将步骤一中所述计算点对应的孔隙压力当量密度ρp和给定的抽汲压力当量密度S1、安全附加值S2及钻井液密度叠加值Δρ代入公式(1)中获取钻井液初始密度ρ0；
ρ0＝ρp+S1+S2-Δρ (1)
式中：ρ0为钻井液初始密度，g/cm3；ρp为孔隙压力当量密度，g/cm3；S1为抽汲压力当量密度，g/cm3；S2为安全附加值，g/cm3；Δρ为钻井液密度叠加值，g/cm3；
步骤三：将步骤二中所述钻井液初始密度ρ0和钻井液密度叠加值Δρ代入公式(2)中获取第i次叠加后钻井液密度ρi(i＝1,2,3,…，k)；
ρi＝ρ0+Δρ(i＝1,2,3,…,k) (2)
式中：ρi为第i次叠加后钻井液密度，g/cm3；ρ0为钻井液初始密度，g/cm3；Δρ为钻井液密度叠加值，g/cm3；
步骤四：将步骤三中所述第i次叠加后钻井液密度ρi(i＝1,2,3,…,k)代入公式(3)中获取环空循环压降Pc；
Figure BDA0003146753820000021
式中：Pc为环空循环压降，MPa；f为环空流体摩阻系数，无量纲；ρi为第i次叠加后钻井液密度，g/cm3；L为计算点井深，m；Q为排量，L/s；D为井眼直径，cm；d为管柱外径，cm；
将步骤三中所述第i次叠加后钻井液密度ρi(i＝1,2,3,…,k)和所述环空循环压降Pc代入公式(4)中获取第i次叠加后井底压力PLi(i＝1,2,3,…,k)；
PLi＝0.0098ρiH+Pc(i＝1,2，3，…,k) (4)
式中：PLi为第i次叠加后井底压力，MPa；ρi为第i次叠加后钻井液密度，g/cm3；H为计算点垂深，m；Pc为环空循环压降，MPa；
步骤五：以裂缝一侧起点为坐标原点O，以钻井液流动方向为x轴，以裂缝宽度方向为y轴建立平面直角坐标系，并沿y轴方向将裂缝内钻井液划分为区域I、区域II和区域III，将计算点对应孔隙压力Pp和步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)代入公式(5)中获取区域II对应宽度tp；
Figure BDA0003146753820000022
式中：tp为区域II对应宽度，mm；τ0为流体屈服值，Pa；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；
步骤六：将步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)和步骤五中所述区域II对应宽度tp代入公式(6)中获取区域I内流体y方向上速度剖面v1(y)；
Figure BDA0003146753820000023
式中：v1(y)为区域I内流体y方向上速度剖面，m/s；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)和步骤五中所述区域II对应宽度tp代入公式(7)中获取区域II内流体y方向上速度剖面v2(y)；
Figure BDA0003146753820000031
式中：v2(y)为区域II内流体y方向上速度剖面，m/s；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)和步骤五中所述区域II对应宽度tp代入公式(8)中获取区域III内流体y方向上速度剖面v3(y)；
Figure BDA0003146753820000032
式中：v3(y)为区域III内流体y方向上速度剖面，m/s；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
步骤七：将步骤六中所述区域I内流体y方向上速度剖面v1(y)代入公式(9)中进行积分获取区域I对应流量q1；
Figure BDA0003146753820000033
式中：q1为区域I对应流量，m3/h；be为裂缝高度，m；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤六中所述区域II内流体y方向上速度剖面v2(y)代入公式(10)中进行积分获取区域II对应流量q2；
Figure BDA0003146753820000034
式中：q2为区域II对应流量，m3/h；be为裂缝高度，m；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤六中所述区域III内流体y方向上速度剖面v3(y)代入公式(11)中进行积分获取区域III对应流量q3；
Figure BDA0003146753820000041
式中：q3为区域III对应流量，m3/h；be为裂缝高度，m；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
步骤八：将步骤七中所述区域I对应流量q1、区域II对应流量q2和区域III对应流量q3代入公式(12)中获取单裂缝漏失速率qa；
qa＝q1+q2+q3 (12)
式中：qa为单裂缝漏失速率，m3/h；q1为区域I对应流量，m3/h；q2为区域II对应流量，m3/h；q3为区域III对应流量q3，m3/h；
步骤九：将步骤一中所述裂缝线密度ke和水平段长度ΔL代入公式(13)中获取第i次叠加后水平段总漏失速率qti(i＝1,2,3,…,k)；
qti＝qakeΔL (13)
式中：qti为第i次叠加后水平段总漏失速率，m3/h；qa为单裂缝漏失速率，m3/h；ke为裂缝线密度，条/m；ΔL为水平段长度，m；
步骤十：判断第i次叠加后水平段总漏失速率qti(i＝1,2,3,…,k)是否等于工程许容漏失速率[qt]，若qti≠[qt]，则返回步骤三计算第i+1次叠加后的钻井液密度ρi+1，重新进行步骤四～步骤十；若qti＝[qt]，则确定此时第i次叠加后井底压力PLi为该计算点对应漏失压力Pf，结束计算。
进一步，所述步骤二～步骤十中环空和裂缝内钻井液均属于屈服幂律流体。
进一步，所述的钻井液密度叠加值Δρ取值范围为0.01～0.05g/cm3。
本发明由于采取以上技术方案，具有以下优点：
(1)该方法综合考虑地层裂缝参数、钻井液性能等因素的影响，分区域确定裂缝内钻井液在y方向上速度剖面及流量分布情况，保证了水平段总漏失速率计算的准确性。
(2)该方法考虑裂缝性地层允许存在井漏的情况，采用叠加更新钻井液密度的方式，使得水平段总漏失速率与工程许容漏失速率相等，从而确定计算点漏失压力，增强了裂缝性地层漏失压力计算方法的适用性。
附图说明
图1是一种井筒漏失压力计算方法的流程图；
图2是裂缝内钻井液流动平面直角坐标系的示意图。
具体实施方式
下面结合附图和具体实施方式对发明进行详细的描述。
如图1所示，本发明提供的一种井筒漏失压力计算方法包括以下步骤：
步骤一：根据所收集的邻井岩心和成像测井资料，确定待钻井水平段的裂缝参数，包括裂缝宽度ae、裂缝高度be、裂缝线密度ke和裂缝长度le；根据钻井工程设计资料，选取漏失压力计算点并确定位置参数，包括计算点井深L、计算点垂深H和水平段长度ΔL；根据邻井水平段进行堵漏作业时的漏失情况，确定工程许容漏失速率[qt]；
步骤二：将步骤一中所述计算点对应的孔隙压力当量密度ρp和给定的抽汲压力当量密度S1、安全附加值S2及钻井液密度叠加值Δρ代入公式(1)中获取钻井液初始密度ρ0；
ρ0＝ρp+S1+S2-Δρ (1)
式中：ρ0为钻井液初始密度，g/cm3；ρp为孔隙压力当量密度，g/cm3；S1为抽汲压力当量密度，g/cm3；S2为安全附加值，g/cm3；Δρ为钻井液密度叠加值，g/cm3；
步骤三：将步骤二中所述钻井液初始密度ρ0和钻井液密度叠加值Δρ代入公式(2)中获取第i次叠加后钻井液密度ρi(i＝1,2，3,…,k)；
ρi＝ρ0+Δρ(i＝1,2,3,…,k) (2)
式中：ρi为第i次叠加后钻井液密度，g/cm3；ρ0为钻井液初始密度，g/cm3；Δρ为钻井液密度叠加值，g/cm3；
步骤四：将步骤三中所述第i次叠加后钻井液密度ρi(i＝1,2,3,…,k)代入公式(3)中获取环空循环压降Pc；
Figure BDA0003146753820000051
式中：Pc为环空循环压降，MPa；f为环空流体摩阻系数，无量纲；ρi为第i次叠加后钻井液密度，g/cm3；L为计算点井深，m；Q为排量，L/s；D为井眼直径，cm；d为管柱外径，cm；
将步骤三中所述第i次叠加后钻井液密度ρi(i＝1，2，3，…，k)和所述环空循环压降Pc代入公式(4)中获取第i次叠加后井底压力PLi(i＝1,2，3，…，k)；
PLi＝0.0098ρiH+Pc(i＝1，2,3,…,k) (4)
式中：PLi为第i次叠加后井底压力，MPa；ρi为第i次叠加后钻井液密度，g/cm3；H为计算点垂深，m；Pc为环空循环压降，MPa；
步骤五：以裂缝一侧起点为坐标原点O，以钻井液流动方向为x轴，以裂缝宽度方向为y轴建立平面直角坐标系，并沿y轴方向将裂缝内钻井液划分为区域I、区域II和区域III，将计算点对应孔隙压力Pp和步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)代入公式(5)中获取区域II对应宽度tp；
Figure BDA0003146753820000052
式中：tp为区域II对应宽度，mm；τ0为流体屈服值，Pa；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；
步骤六：将步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)和步骤五中所述区域II对应宽度tp代入公式(6)中获取区域I内流体y方向上速度剖面v1(y)；
Figure BDA0003146753820000061
式中：v1(y)为区域I内流体y方向上速度剖面，m/s；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)和步骤五中所述区域II对应宽度tp代入公式(7)中获取区域II内流体y方向上速度剖面v2(y)；
Figure BDA0003146753820000062
式中：v2(y)为区域II内流体y方向上速度剖面，m/s；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)和步骤五中所述区域II对应宽度tp代入公式(8)中获取区域III内流体y方向上速度剖面v3(y)；
Figure BDA0003146753820000063
式中：v3(y)为区域III内流体y方向上速度剖面，m/s；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
步骤七：将步骤六中所述区域I内流体y方向上速度剖面v1(y)代入公式(9)中进行积分获取区域I对应流量q1；
Figure BDA0003146753820000064
式中：q1为区域I对应流量，m3/h；be为裂缝高度，m；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤六中所述区域II内流体y方向上速度剖面v2(y)代入公式(10)中进行积分获取区域II对应流量q2；
Figure BDA0003146753820000071
式中：q2为区域II对应流量，m3/h；be为裂缝高度，m；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤六中所述区域III内流体y方向上速度剖面v3(y)代入公式(11)中进行积分获取区域III对应流量q3；
Figure BDA0003146753820000072
式中：q3为区域III对应流量，m3/h；be为裂缝高度，m；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
步骤八：将步骤七中所述区域I对应流量q1、区域II对应流量q2和区域III对应流量q3代入公式(12)中获取单裂缝漏失速率qa；
qa＝q1+q2+q3 (12)
式中：qa为单裂缝漏失速率，m3/h；q1为区域I对应流量，m3/h；q2为区域II对应流量，m3/h；q3为区域III对应流量q3，m3/h；
步骤九：将步骤一中所述裂缝线密度ke和水平段长度ΔL代入公式(13)中获取第i次叠加后水平段总漏失速率qti(i＝1,2,3,…,k)；
qti＝qakeΔL (13)
式中：qti为第i次叠加后水平段总漏失速率，m3/h；qa为单裂缝漏失速率，m3/h；ke为裂缝线密度，条/m；ΔL为水平段长度，m；
步骤十：判断第i次叠加后水平段总漏失速率qti(i＝1,2,3,…,k)是否等于工程许容漏失速率[qt]，若qti≠[qt]，则返回步骤三计算第i+1次叠加后的钻井液密度ρi+1，重新进行步骤四～步骤十；若qti＝[qt]，则确定此时第i次叠加后井底压力PLi为该计算点对应漏失压力Pf，结束计算。
进一步，所述步骤二～步骤十中环空和裂缝内钻井液均属于屈服幂律流体。
进一步，所述的钻井液密度叠加值Δρ取值范围为0.01～0.05g/cm3。
下面结合实施例对本发明的具体实施方式做进一步的描述，并不因此将本发明限制在所述的实施例范围之中。
实施例：
一种井筒漏失压力计算方法，包括以下步骤：
步骤一：根据所收集的邻井岩心和成像测井资料，确定待钻井裂缝宽度ae＝0.2mm、裂缝高度be＝0.5mm、裂缝线密度ke＝0.25条/m和裂缝长度le＝2.5m；根据钻井工程设计资料选取的计算点井深L＝4000m、计算点垂深H＝3800m和水平段长度ΔL＝54m；根据邻井水平段漏失速率达到5～7m3/h便需要进行堵漏作业，因此确定工程许容漏失速率[qt]＝6m3/h；
步骤二：将步骤一中所述计算点对应的孔隙压力当量密度ρp＝1.15g/cm3和给定的抽汲压力当量密度S1＝0.05g/cm3、安全附加值S2＝0.04g/cm3及钻井液密度叠加值Δρ＝0.04g/cm3代入公式(1)中得到钻井液初始密度ρ0＝1.20g/cm3；
步骤三：将步骤二中的钻井液初始密度ρ0和钻井液密度叠加值Δρ代入公式(2)中得到第1次叠加后钻井液密度ρ1＝1.24g/cm3；
步骤四：根据相关钻井资料可知，环空流体摩阻系数f＝0.03、排量Q＝9L/s、井眼直径D＝15.4cm、管柱外径d＝8.89cm，结合步骤一中计算点井深L＝4000m和步骤三中第1次叠加后钻井液密度ρ1，代入公式(3)中得到环空循环压降Pc＝2.62MPa；
将步骤三中所述第1次叠加后钻井液密度ρ1和环空循环压降Pc代入公式(4)中得到第1次叠加后井底压力PL1＝48.80MPa；
步骤五：以裂缝一侧起点为坐标原点O，以钻井液流动方向为x轴，以裂缝宽度方向为y轴建立平面直角坐标系，并沿y轴方向将裂缝内钻井液划分为区域I、区域II和区域III，已知流体屈服值τ0＝11MPa，将计算点对应孔隙压力Pp＝42.83MPa和步骤四中计算点对应井底压力PL代入公式(5)中得到区域II对应宽度tp＝0.09mm；
步骤六：已知流性指数n＝0.5和流体稠度系数K＝1Pa·Sn，将步骤四中计算点对应井底压力PL和步骤五中区域II对应宽度tp分别代入公式(6)、公式(7)和公式(8)中得到区域I内流体y方向上速度剖面v1(y)、区域II内流体y方向上速度剖面v2(y)以及区域III内流体y方向上速度剖面v3(y)；
步骤七：将步骤六中所述区域I内流体y方向上速度剖面v1(y)代入公式(9)中进行积分得到区域I对应流量q1＝0.022m3/h；
将步骤六中所述区域II内流体y方向上速度剖面v2(y)代入公式(10)中进行积分获取得到区域II对应流量q2＝0.049m3/h；
将步骤六中所述区域III内流体y方向上速度剖面v3(y)代入公式(11)中进行积分获取得到区域III对应流量q3＝0.022m3/h；
步骤八：将步骤七中区域I对应流量q1、区域II对应流量q2和区域III对应流量q3代入公式(12)中得到单裂缝漏失速率qa＝0.093m3/h；
步骤九：将步骤一中所述裂缝线密度ke和水平段长度ΔL代入公式(13)中得到第1次叠加后水平段总漏失速率qti＝1.256m3/h；
步骤十：根据步骤一中工程许容漏失速率[qt]＝6m3/h，判断出qt1≠[qt]，返回步骤三计算出第2次叠加后的钻井液密度ρ2＝1.28g/cm3，重复进行步骤四～步骤九得到第2次叠加后水平段总漏失速率qt2＝3.127m3/h，判断出qt2≠[qt]；返回步骤三计算出第3次叠加后的钻井液密度ρ3＝1.32g/cm3，重复进行步骤四～步骤十得到第3次叠加后井底压力PL3＝51.95MPa，以及水平段总漏失速率qt3＝6.048m3/h，近似判断为qt2＝[qt]，则确定该计算点对应漏失压力Pf＝PL3＝51.95MPa，结束计算。
通过上述的计算，实施例中井深L为4000m处水平段的漏失压力Pf为51.95MPa。
该方法考虑裂缝性地层实际钻井过程中允许存在一定井漏的情况，分区域确定裂缝内钻井液速度剖面及流量分布情况，同时循环叠加更新钻井液密度使得水平段漏失速率等于工程许容漏失速率，以此提出了一种适用于裂缝性地层水平段的漏失压力计算方法。该方法基于邻井岩心、成像测井资料和已知工程参数，可对裂缝性地层水平段的漏失压力进行快速、准确地预测，为裂缝性地层安全钻井提供技术支撑，避免井下复杂事故的发生。

Claims (3)

1.一种井筒漏失压力计算方法，其特征在于，包含以下步骤：
步骤一：根据所收集的邻井岩心和成像测井资料，确定待钻井水平段的裂缝参数，包括裂缝宽度ae、裂缝高度be、裂缝线密度ke和裂缝长度le；根据钻井工程设计资料，选取漏失压力计算点并确定位置参数，包括计算点井深L、计算点垂深H和水平段长度ΔL；根据邻井水平段进行堵漏作业时的漏失情况，确定工程许容漏失速率[qt]；
步骤二：将步骤一中所述计算点对应的孔隙压力当量密度ρp和给定的抽汲压力当量密度S1、安全附加值S2及钻井液密度叠加值Δρ代入公式(1)中获取钻井液初始密度ρ0；
ρ0＝ρp+S1+S2-Δρ (1)
式中：ρ0为钻井液初始密度，g/cm3；ρp为孔隙压力当量密度，g/cm3；S1为抽汲压力当量密度，g/cm3；S2为安全附加值，g/cm3；Δρ为钻井液密度叠加值，g/cm3；
步骤三：将步骤二中所述钻井液初始密度ρ0和钻井液密度叠加值Δρ代入公式(2)中获取第i次叠加后钻井液密度ρi(i＝1,2,3,…,k)；
ρi＝ρ0+Δρ(i＝1,2,3,…,k) (2)
式中：ρi为第i次叠加后钻井液密度，g/cm3；ρ0为钻井液初始密度，g/cm3；Δρ为钻井液密度叠加值，g/cm3；
步骤四：将步骤三中所述第i次叠加后钻井液密度ρi(i＝1,2，3，…，k)代入公式(3)中获取环空循环压降Pc；
Figure FDA0003146753810000011
式中：Pc为环空循环压降，MPa；f为环空流体摩阻系数，无量纲；ρi为第i次叠加后钻井液密度，g/cm3；L为计算点井深，m；Q为排量，L/s；D为井眼直径，cm；d为管柱外径，cm；
将步骤三中所述第i次叠加后钻井液密度ρi(i＝1,2,3,…,k)和所述环空循环压降Pc代入公式(4)中获取第i次叠加后井底压力PLi(i＝1,2,3,…,k)；
PLi＝0.0098ρiH+Pc(i＝1,2,3,…,k) (4)
式中：PLi为第i次叠加后井底压力，MPa；ρi为第i次叠加后钻井液密度，g/cm3；H为计算点垂深，m；Pc为环空循环压降，MPa；
步骤五：以裂缝一侧起点为坐标原点O，以钻井液流动方向为x轴，以裂缝宽度方向为y轴建立平面直角坐标系，并沿y轴方向将裂缝内钻井液划分为区域I、区域II和区域III，将计算点对应孔隙压力Pp和步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)代入公式(5)中获取区域II对应宽度tp；
Figure FDA0003146753810000021
式中：tp为区域II对应宽度，mm；τ0为流体屈服值，Pa；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；
步骤六：将步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)和步骤五中所述区域II对应宽度tp代入公式(6)中获取区域I内流体y方向上速度剖面v1(y)；
Figure FDA0003146753810000022
式中：v1(y)为区域I内流体y方向上速度剖面，m/s；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)和步骤五中所述区域II对应宽度tp代入公式(7)中获取区域II内流体y方向上速度剖面v2(y)；
Figure FDA0003146753810000023
式中：v2(y)为区域II内流体y方向上速度剖面，m/s；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤四中所述第i次叠加后井底压力PLi(i＝1,2,3,…,k)和步骤五中所述区域II对应宽度tp代入公式(8)中获取区域III内流体y方向上速度剖面v3(y)；
Figure FDA0003146753810000024
式中：v3(y)为区域III内流体y方向上速度剖面，m/s；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
步骤七：将步骤六中所述区域I内流体y方向上速度剖面v1(y)代入公式(9)中进行积分获取区域I对应流量q1；
Figure FDA0003146753810000025
式中：q1为区域I对应流量，m3/h；be为裂缝高度，m；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤六中所述区域II内流体y方向上速度剖面v2(y)代入公式(10)中进行积分获取区域II对应流量q2；
Figure FDA0003146753810000031
式中：q2为区域II对应流量，m3/h；be为裂缝高度，m；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
将步骤六中所述区域III内流体y方向上速度剖面v3(y)代入公式(11)中进行积分获取区域III对应流量q3；
Figure FDA0003146753810000032
式中：q3为区域III对应流量，m3/h；be为裂缝高度，m；n为流性指数，无量纲；PLi为第i次叠加后井底压力，MPa；Pp为计算点对应孔隙压力，MPa；le为裂缝长度，m；K为流体稠度系数，Pa·Sn；ae为裂缝宽度，mm；tp为区域II对应宽度，mm；
步骤八：将步骤七中所述区域I对应流量q1、区域II对应流量q2和区域III对应流量q3代入公式(12)中获取单裂缝漏失速率qa；
qa＝q1+q2+q3 (12)
式中：qa为单裂缝漏失速率，m3/h；q1为区域I对应流量，m3/h；q2为区域II对应流量，m3/h；q3为区域III对应流量q3，m3/h；
步骤九：将步骤一中所述裂缝线密度ke和水平段长度ΔL代入公式(13)中获取第i次叠加后水平段总漏失速率qti(i＝1,2,3,…,k)；
qti＝qakeΔL (13)
式中：qti为第i次叠加后水平段总漏失速率，m3/h；qa为单裂缝漏失速率，m3/h；ke为裂缝线密度，条/m；ΔL为水平段长度，m；
步骤十：判断第i次叠加后水平段总漏失速率qti(i＝1,2,3，…,k)是否等于工程许容漏失速率[qt]，若qti≠[qt]，则返回步骤三计算第i+1次叠加后的钻井液密度ρi+1，重新进行步骤四～步骤十；若qti＝[qt]，则确定此时第i次叠加后井底压力PLi为该计算点对应漏失压力Pf，结束计算。
2.根据权利要求1所述的一种井筒漏失压力计算方法，其特征在于：所述步骤二～步骤十中环空和裂缝内钻井液均属于屈服幂律流体。
3.根据权利要求1所述的一种井筒漏失压力计算方法，其特征在于：所述的钻井液密度叠加值Δρ取值范围为0.01～0.05g/cm3。
CN202110753840.5A 2021-07-04 2021-07-04 一种井筒漏失压力计算方法 Active CN113323658B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110753840.5A CN113323658B (zh) 2021-07-04 2021-07-04 一种井筒漏失压力计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110753840.5A CN113323658B (zh) 2021-07-04 2021-07-04 一种井筒漏失压力计算方法

Publications (2)

Publication Number Publication Date
CN113323658A true CN113323658A (zh) 2021-08-31
CN113323658B CN113323658B (zh) 2022-03-25

Family

ID=77425709

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110753840.5A Active CN113323658B (zh) 2021-07-04 2021-07-04 一种井筒漏失压力计算方法

Country Status (1)

Country Link
CN (1) CN113323658B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115951422A (zh) * 2023-03-14 2023-04-11 北京阳光杰科科技股份有限公司 构建天然裂缝漏失压力模型的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090188718A1 (en) * 2008-01-30 2009-07-30 M-I L.L.C. Methods of detecting, preventing, and remediating lost circulation
CN108825156A (zh) * 2017-05-05 2018-11-16 中国石油化工股份有限公司 一种用于控压钻井的气侵控制方法
CN110826245A (zh) * 2019-11-15 2020-02-21 西南石油大学 一种井漏与堵漏的室内模拟评价方法
CN111236925A (zh) * 2020-01-19 2020-06-05 中国石油大学(华东) 一种裂缝性地层钻井液漏失量预测方法
CN111734399A (zh) * 2020-04-01 2020-10-02 四川维泰科创石油设备制造有限公司 一种钻井过程中智能堵漏方法及***
US20210002969A1 (en) * 2017-07-11 2021-01-07 Equinor Energy As Method and System for Monitoring Influx and Loss Events in a Wellbore

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090188718A1 (en) * 2008-01-30 2009-07-30 M-I L.L.C. Methods of detecting, preventing, and remediating lost circulation
CN108825156A (zh) * 2017-05-05 2018-11-16 中国石油化工股份有限公司 一种用于控压钻井的气侵控制方法
US20210002969A1 (en) * 2017-07-11 2021-01-07 Equinor Energy As Method and System for Monitoring Influx and Loss Events in a Wellbore
CN110826245A (zh) * 2019-11-15 2020-02-21 西南石油大学 一种井漏与堵漏的室内模拟评价方法
CN111236925A (zh) * 2020-01-19 2020-06-05 中国石油大学(华东) 一种裂缝性地层钻井液漏失量预测方法
CN111734399A (zh) * 2020-04-01 2020-10-02 四川维泰科创石油设备制造有限公司 一种钻井过程中智能堵漏方法及***

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
EVGENII A. KANIN等: "A radial hydraulic fracture with pressure-dependent leak-off", 《JOURNAL OF THE MECHANICS AND PHYSICS OF SOLIDS》 *
李双贵等: "麦盖提１区碳酸盐岩地层钻井漏失分析及优化", 《石油天然气学报》 *
李大奇等: "基于漏失机理的碳酸盐岩地层漏失压力模型", 《石油学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115951422A (zh) * 2023-03-14 2023-04-11 北京阳光杰科科技股份有限公司 构建天然裂缝漏失压力模型的方法
CN115951422B (zh) * 2023-03-14 2023-05-09 北京阳光杰科科技股份有限公司 构建天然裂缝漏失压力模型的方法

Also Published As

Publication number Publication date
CN113323658B (zh) 2022-03-25

Similar Documents

Publication Publication Date Title
CN108868748B (zh) 一种页岩气水平井重复压裂裂缝开启压力的计算方法
CN109538185B (zh) 一种耦合温度场条件下的多层套管井井筒完整性分析模型
CN112627733B (zh) 深水控压钻井水力参数实时优化方法及设备
CN105840187A (zh) 致密性油藏水平井分段压裂产能计算方法
CN106522928B (zh) 一种酸化压裂后停泵测井口压降不稳定试井方法
CN110656915B (zh) 一种页岩气多段压裂水平井多工作制度产能预测方法
CN109707336B (zh) 一种基于环空抽吸泵的控压固井方法
CN108708711A (zh) 一种准确确定漏失层位的方法
CN105672997A (zh) 钻井液地层漏失量监测方法
CN113323658B (zh) 一种井筒漏失压力计算方法
CN109667564A (zh) 一种海上稠油油田蒸汽吞吐开发定向井产能的确定方法
CN105804713B (zh) 一种快速确定注水井各小层井口注水启动压力的方法
Zhang et al. Prediction model of casing wear shape and residual strength under compound modes
CN110714755B (zh) 水驱油藏剩余油二次富集速度快速预测方法
CN108487905B (zh) 一种页岩气水平井压裂参数的优化方法
Liu et al. Casing shearing failure in shale gas wells due to fault slippage caused by hydraulic fracturing-case study
US20080230221A1 (en) Methods and systems for monitoring near-wellbore and far-field reservoir properties using formation-embedded pressure sensors
CN115293066B (zh) 考虑地层渗流传热效应的气井温度场计算方法
Marbun et al. Improvement of loads calculation of the perforated liner in a geothermal production well
RU2447272C1 (ru) Способ разработки нефтяной залежи массивного типа
CN113550727B (zh) 一种裂缝发育区体积压裂套管强度计算方法
CN114810023A (zh) 一种适应于密集井眼防碰的施工方法
CN114647913A (zh) 一种基于Mohr-Coulomb准则的最大水平主地应力获取方法
Zhou et al. Analysis of leak-off tests in shallow marine sediments
Xie et al. Analysis of the applicability of a risk quantitative evaluation method to high temperature-pressure drilling engineering

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