CN106682384B - 井眼扩径率计算方法及装置 - Google Patents

井眼扩径率计算方法及装置 Download PDF

Info

Publication number
CN106682384B
CN106682384B CN201610806508.XA CN201610806508A CN106682384B CN 106682384 B CN106682384 B CN 106682384B CN 201610806508 A CN201610806508 A CN 201610806508A CN 106682384 B CN106682384 B CN 106682384B
Authority
CN
China
Prior art keywords
rock
mpa
borehole wall
wellbore
coefficient
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
CN201610806508.XA
Other languages
English (en)
Other versions
CN106682384A (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201610806508.XA priority Critical patent/CN106682384B/zh
Publication of CN106682384A publication Critical patent/CN106682384A/zh
Application granted granted Critical
Publication of CN106682384B publication Critical patent/CN106682384B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • E21B7/00Special methods or apparatus for drilling
    • E21B7/28Enlarging drilled holes, e.g. by counterboring
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

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)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供了一种井眼扩径率计算方法及装置。该方法包括:通过室内实验得到井壁的取样岩石的线弹性参数;通过室内实验得到所述井壁的取样岩石的塑性参数;通过地层资料和室内实验获取所述井壁的基础应力参数;根据所述线弹性参数、所述塑性参数、所述基础应力参数和钻井液密度,计算井眼扩径率。本发明可以在给定钻井液密度的情况下得到钻井的井眼扩径率,因而可以准确选择并配定合理的钻井液密度,为优质快速钻井提供技术支持。

Description

井眼扩径率计算方法及装置
技术领域
本发明涉及石油钻井技术领域,特别涉及一种井眼扩径率计算方法及装置。
背景技术
在石油钻井过程中,一方面为了提高钻井速度需要较小的钻井液密,另一方面为了防止井眼垮塌又需要维持合适的钻井液密度。目前通常的做法是以“坍塌压力”作为井筒钻井液的允许的最小密度,即刚好维持井壁平衡的钻井液密度,但实际情况下在井壁垮塌率在10%~20%情况下,甚至30%,仍可以维持钻井的进行,显然以“坍塌压力”作为设计最小钻井液密度依据显得保守,不利于钻井提速的需求。
针对石油钻井井眼扩径率的计算,目前只有通过有限元分析的方法,还没有实际的计算公式。已经有文献给出了可以计算维持扩径率为0%的钻井液密度(也就是坍塌压力计算公式),但是实际上一旦需要计算井眼扩径就无法满足要求。
因此,在石油钻井技术领域,缺乏一种在不同钻井液密度下有效计算井眼扩径率的方法。
发明内容
为了解决上述的技术问题,本发明提供一种井眼扩径率计算方法及装置,有效计算不同钻井液密度下的井眼扩径率,从而可以准确配置合适的钻井液密度,为优质快速钻井提供技术支持。
具体而言,包括以下的技术方案:
一种井眼扩径率计算方法,包括:
通过室内实验得到井壁的取样岩石的线弹性参数;
通过室内实验得到所述井壁的取样岩石的塑性参数;
通过地层资料和室内实验获取所述井壁的基础应力参数;
根据所述线弹性参数、所述塑性参数、所述基础应力参数和钻井液密度,计算井眼扩径率。
可选择地,所述井壁的取样岩石的线弹性参数包括:岩石弹性阶段的抗压强度、泊松比与岩石内摩擦系数。
可选择地,所述井壁的取样岩石的塑性参数包括:岩石产生裂缝时的残余强度、岩石强度软化模量梯度和塑性软化区围岩扩容梯度。
可选择地,所述井壁的基础应力参数包括:最小地应力、有效应力系数、孔隙压力、地应力系数和井深度。
可选择地,所述井眼扩径率的计算公式为:
其中:
式中:γ,为井眼扩径率,%;ρ钻井液密度,g/cm3;H为井深度,m;Pp为孔隙压力,MPa;α为有效应力系数;σh为最小地应力,MPa;λ为地应力系数;σc为岩石弹性阶段抗压强度,MPa;μ为泊松比;E为弹性模量,MPa;K岩石摩擦系数;σcs为产生裂纹时的残余强度,MPa;Q为岩石强度软化模量梯度;MPa;η1为塑性软化区围岩扩容梯度。
一种井眼扩径率计算装置,包括:
线弹性参数获取模块,用于通过室内实验得到井壁的取样岩石的线弹性参数;
塑性参数获取模块,用于通过室内实验得到所述井壁的取样岩石的塑性参数;
基础应力参数获取模块,通过地层资料和室内实验获取所述井壁的基础应力参数;
扩径率计算模块,用于根据所述线弹性参数、所述塑性参数、所述基础应力参数和钻井液密度,计算井眼扩径率。
可选择地,所述井壁的取样岩石的线弹性参数包括:岩石弹性阶段的抗压强度、泊松比与岩石内摩擦系数。
可选择地,所述井壁的取样岩石的塑性参数包括:岩石产生裂缝时的残余强度、岩石强度软化模量梯度和塑性软化区围岩扩容梯度。
可选择地,所述井壁的基础应力参数包括:最小地应力、有效应力系数、孔隙压力、地应力系数和井深度。
可选择地,所述井眼扩径率的计算公式为:
其中:
式中:γ,为井眼扩径率,%;ρ钻井液密度,g/cm3;H为井深度,m;Pp为孔隙压力,MPa;α为有效应力系数;σh为最小地应力,MPa;λ为地应力系数;σc为岩石弹性阶段抗压强度,MPa;μ为泊松比;E为弹性模量,MPa;K岩石摩擦系数;σcs为产生裂纹时的残余强度,MPa;Q为岩石强度软化模量梯度;MPa;η1为塑性软化区围岩扩容梯度。
本发明实施例提供的技术方案的有益效果:
本发明的井眼扩径率计算方法和井眼扩径率计算装置,引入井壁岩石的线弹性参数、井壁岩石的塑性参数、井壁的基础应力参数和钻井液密度来计算井眼扩径率,从而可以在给定钻井液密度的情况下得到钻井的井眼扩径率,因而可以准确选择并配定合理的钻井液密度,为优质快速钻井提供技术支持。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是根据本发明一实施例的井眼扩径率计算方法的流程图;
图2是根据本发明一实施例的井眼扩径率计算装置的框图;
图3是本发明一实例测得的井眼扩径率实际变化情况的曲线图。
具体实施方式
为使本发明的技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明一实施例提供了一种井眼扩径率计算方法,如图1所示,包括:
步骤S1:通过室内实验得到井壁的取样岩石的线弹性参数;
步骤S2:通过室内实验得到井壁的取样岩石的塑性参数;
步骤S3:通过地层资料获取井壁的基础应力参数;
步骤S4:根据线弹性参数、塑性参数、基础应力参数和钻井液密度,计算井眼扩径率。
其中,井壁的取样岩石的线弹性参数包括:岩石弹性阶段的抗压强度、泊松比、弹性模量与岩石内摩擦系数。
泊松比是指岩石在单向受拉或受压时,横向正应变与轴向正应变的绝对值的比值,也叫横向变形系数,它是反映材料横向变形的弹性常数。
岩石在弹性变形阶段,其应力和应变成正比例关系(即符合胡克定律),该比例系数称为弹性模量。
井壁的取样岩石的塑性参数包括:岩石产生裂缝时的残余强度、岩石强度软化模量梯度和塑性软化区围岩扩容梯度。
岩石产生裂缝时的残余强度是指岩石在破坏后所残留的抵抗外荷的能力。可在应力—应变全过程曲线上求得。
井壁的取样岩石的线弹性参数和塑性参数都是岩石的力学参数。岩石力学参数的测试可采用以下方法:将岩样加工成的圆柱形型岩心(例如为Φ25×50mm)进行称重记录,通过岩石三轴压缩实验测试计算得出岩石的弹性模量、泊松比、弹性抗压强度、内摩擦系数等弹性参数和岩石软化模量、扩容梯度、残余强度等。从而完成步骤S1、S2。步骤S1、S2中的室内实验实际上是同一次实验,也可以看成是一次实验的先后两部分。
井壁的基础应力参数包括:最小地应力、有效应力系数、孔隙压力、地应力系数和井深度。
孔隙压力指的是岩石内孔隙流体承受的压力。其值为孔隙水压力和孔隙气压力之和。
其中,井深度、孔隙压力、有效应力系数等可以从测井资料获得;而最小地应力、地应力系数可通过室内实验确定,例如可通过波速各向异性法和粘滞剩磁法来确定各主地应力的方向,并通过差应变法测试和计算分析得出地应力大小等参数。从而完成步骤S3。
步骤S4中,钻井液密度根据选定或配备的钻井液测定,也可预设一个钻井液密度。井眼扩径率的计算公式为:
其中:
式中:γ,为井眼扩径率,%;ρ钻井液密度,g/cm3;H为井深度,m;Pp为孔隙压力,MPa;α为有效应力系数;σh为最小地应力,MPa;λ为地应力系数;σc为岩石弹性阶段抗压强度,MPa;μ为泊松比;E为弹性模量,MPa;K岩石摩擦系数;σcs为产生裂纹时的残余强度,MPa;Q为岩石强度软化模量梯度;MPa;η1为塑性软化区围岩扩容梯度。
本领域技术人员应当理解,上述公式只是例示性的,还可以以其他公式计算,但是只要是根据线弹性参数、塑性参数、基础应力参数和钻井液密度来计算井眼扩径率,都在本发明的保护范围之内。
本实施例的井眼扩径率计算方法,引入井壁岩石的线弹性参数、井壁岩石的塑性参数、井壁的基础应力参数和钻井液密度来计算井眼扩径率,从而可以在给定钻井液密度的情况下得到钻井的井眼扩径率,因而可以准确选择并配定合理的钻井液密度,为优质快速钻井提供技术支持。
本发明另一实施例提供了一种井眼扩径率计算装置,如图2所示,包括:
线弹性参数获取模块1,用于通过室内实验得到井壁的取样岩石的线弹性参数;
塑性参数获取模块2,用于通过室内实验得到所述井壁的取样岩石的塑性参数;
基础应力参数获取模块3,通过地层资料获取所述井壁的基础应力参数;
扩径率计算模块4,用于根据所述线弹性参数、所述塑性参数、所述基础应力参数和钻井液密度,计算井眼扩径率。
其中,井壁的取样岩石的线弹性参数包括:岩石弹性阶段的抗压强度、泊松比、弹性模量与岩石内摩擦系数。
井壁的取样岩石的塑性参数包括:岩石产生裂缝时的残余强度、岩石强度软化模量梯度和塑性软化区围岩扩容梯度。
岩石产生裂缝时的残余强度是指岩石在破坏后所残留的抵抗外荷的能力。可在应力—应变全过程曲线上求得。
井壁的取样岩石的线弹性参数和塑性参数都是岩石的力学参数。岩石力学参数的测试可采用以下方法:将岩样加工成的圆柱形型岩心(例如为Φ25×50mm)进行称重记录,通过岩石三轴压缩实验测试计算得出岩石的弹性模量、泊松比、弹性抗压强度、内摩擦系数等弹性参数和岩石软化模量、扩容梯度、残余强度等。获得上述参数后,就实现了线弹性参数获取模块1和塑性参数获取模块2,例如可在相应的实验装置上加上实验数据传送接口就构成了上述模块。
井壁的基础应力参数包括:最小地应力、有效应力系数、孔隙压力、地应力系数和井深度。
其中,井深度、孔隙压力、有效应力系数等可以从测井资料获得;而最小地应力、地应力系数可通过室内实验确定,例如可通过波速各向异性法和粘滞剩磁法来确定各主地应力的方向,并通过差应变法测试和计算分析得出地应力大小等参数。获得上述参数后,就实现了基础应力参数获取模块3。
扩径率计算模块4中,钻井液密度根据选定或配备的钻井液测定,也可预设一个钻井液密度。井眼扩径率的计算公式为:
其中:
式中:γ,为井眼扩径率,%;ρ钻井液密度,g/cm3;H为井深度,m;Pp为孔隙压力,MPa;α为有效应力系数;σh为最小地应力,MPa;λ为地应力系数;σc为岩石弹性阶段抗压强度,MPa;μ为泊松比;E为弹性模量,MPa;K岩石摩擦系数;σcs为产生裂纹时的残余强度,MPa;Q为岩石强度软化模量梯度;MPa;η1为塑性软化区围岩扩容梯度。
本实施例的井眼扩径率计算装置,引入井壁岩石的线弹性参数、井壁岩石的塑性参数、井壁的基础应力参数和钻井液密度来计算井眼扩径率,从而可以在给定钻井液密度的情况下得到钻井的井眼扩径率,因而可以准确选择并配定合理的钻井液密度,为优质快速钻井提供技术支持。
下面通过一个实例来验证本发明的效果。
某区块井深3000m处,最小地应力为43MPa,最大地应力为64.5MPa(最大地应力除以最小地应力就是地应力系数λ),地层压力系数1.10(地层压力系数乘以深度可以得到空隙压力Pp)。根据该地区取芯获得的岩石力学参数计算不同钻井液密度对应的垮塌率。
表1岩石力学参数
现场使用了两种钻井液密度1.30g/cm3和1.15g/cm3分别在X井和Y井进行了试验性钻井,把上表的岩石力学参数,以及上文提及的地应力数据、深度数据、井下液柱压力(可由Pi=ρH计算),代入扩径率公式可计算密度1.30g/cm3和1.15g/cm3两种钻井液密度下的扩径率。
表2钻井液密度与预测扩径率
现场钻井钻头直径311.2mm,如果井眼不扩径,那么井筒内径为311.2mm。在实际钻井过程中:参见图3,在X井中,使用密度1.30g/cm3的钻井液的井径变化率基本在3~8%左右(预测5%);在Y井中,使用密度1.13g/cm3钻井液的井径变化率在13~17%左右(预测15%),可见在不同钻井液密度下,通过模型计算出的井眼扩径率与实际井眼扩径率很吻合,说明建立的模型的可靠性。
以上所述仅是为了便于本领域的技术人员理解本发明的技术方案,并不用以限制本发明。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.一种井眼扩径率计算方法,其特征在于,包括:
通过室内实验得到井壁的取样岩石的线弹性参数;
通过室内实验得到所述井壁的取样岩石的塑性参数;
通过地层资料和室内实验获取所述井壁的基础应力参数;
根据所述线弹性参数、所述塑性参数、所述基础应力参数和钻井液密度,计算井眼扩径率,
所述井壁的取样岩石的线弹性参数包括:岩石弹性阶段的抗压强度、泊松比与岩石内摩擦系数;所述井壁的取样岩石的塑性参数包括:岩石产生裂缝时的残余强度、岩石强度软化模量梯度和塑性软化区围岩扩容梯度;所述井壁的基础应力参数包括:最小地应力、有效应力系数、孔隙压力、地应力系数和井深度;所述井眼扩径率的计算公式为:
其中:
式中:γ,为井眼扩径率,%;Pi=ρH,为井下液柱压力;ρ为钻井液密度,g/cm3;H为井深度,m;Pp为孔隙压力,MPa;α为有效应力系数;σh为最小地应力,MPa;λ为地应力系数;σc为岩石弹性阶段抗压强度,MPa;μ为泊松比;E为弹性模量,MPa;K岩石摩擦系数;σcs为产生裂纹时的残余强度,MPa;Q为岩石强度软化模量梯度;MPa;η1为塑性软化区围岩扩容梯度。
2.一种井眼扩径率计算装置,其特征在于,包括:
线弹性参数获取模块,用于通过室内实验得到井壁的取样岩石的线弹性参数;
塑性参数获取模块,用于通过室内实验得到所述井壁的取样岩石的塑性参数;
基础应力参数获取模块,通过地层资料和室内实验获取所述井壁的基础应力参数;
扩径率计算模块,用于根据所述线弹性参数、所述塑性参数、所述基础应力参数和钻井液密度,计算井眼扩径率;所述井壁的取样岩石的线弹性参数包括:岩石弹性阶段的抗压强度、泊松比与岩石内摩擦系数;所述井壁的取样岩石的塑性参数包括:岩石产生裂缝时的残余强度、岩石强度软化模量梯度和塑性软化区围岩扩容梯度;所述井壁的基础应力参数包括:最小地应力、有效应力系数、孔隙压力、地应力系数和井深度;所述井眼扩径率的计算公式为:
其中:
式中:γ,为井眼扩径率,%;Pi=ρH,为井下液柱压力;ρ为钻井液密度,g/cm3;H为井深度,m;Pp为孔隙压力,MPa;α为有效应力系数;σh为最小地应力,MPa;λ为地应力系数;σc为岩石弹性阶段抗压强度,MPa;μ为泊松比;E为弹性模量,MPa;K岩石摩擦系数;σcs为产生裂纹时的残余强度,MPa;Q为岩石强度软化模量梯度;MPa;η1为塑性软化区围岩扩容梯度。
CN201610806508.XA 2016-09-06 2016-09-06 井眼扩径率计算方法及装置 Active CN106682384B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610806508.XA CN106682384B (zh) 2016-09-06 2016-09-06 井眼扩径率计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610806508.XA CN106682384B (zh) 2016-09-06 2016-09-06 井眼扩径率计算方法及装置

Publications (2)

Publication Number Publication Date
CN106682384A CN106682384A (zh) 2017-05-17
CN106682384B true CN106682384B (zh) 2019-02-15

Family

ID=58840283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610806508.XA Active CN106682384B (zh) 2016-09-06 2016-09-06 井眼扩径率计算方法及装置

Country Status (1)

Country Link
CN (1) CN106682384B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108316895B (zh) * 2018-02-23 2020-10-16 中国海洋石油集团有限公司 一种获取海域天然气水合物层实时钻井扩径率的方法
CN110792418B (zh) * 2018-08-03 2022-03-01 中国石油天然气股份有限公司 井筒工作液配方优化方法及装置
CN111797490B (zh) * 2019-04-04 2022-11-04 中国石油天然气股份有限公司 钻井液安全密度设计的方法及设备
CN115238861B (zh) * 2022-07-11 2023-11-17 西南石油大学 一种基于井壁垮塌程度约束的安全钻井液密度确定方法
CN115324571B (zh) * 2022-08-02 2023-10-10 西南石油大学 基于钻录井信息的复杂地层井壁垮塌定量预测方法与装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4633449A (en) * 1982-10-15 1986-12-30 Ingram John D Method and apparatus for indirect determination of shear velocity from guided modes
US7482309B2 (en) * 2003-11-24 2009-01-27 Halliburton Energy Services, Inc. Methods of drilling wellbores using variable density fluids comprising coated elastic particles
CN102128025B (zh) * 2011-01-20 2013-03-13 中国石油天然气股份有限公司 一种盐膏层造斜钻井液密度确定方法
US9134291B2 (en) * 2012-01-26 2015-09-15 Halliburton Energy Services, Inc. Systems, methods and devices for analyzing drilling fluid
CN103390112B (zh) * 2013-07-29 2016-08-10 中国石油大学(华东) 考虑煤层地应力不确定性的钻井液密度设计方法
CN103573285A (zh) * 2013-10-25 2014-02-12 山西潞安环保能源开发股份有限公司五阳煤矿 一种用于三低煤层抽取煤层气的圆柱式开采方法
CN103775074B (zh) * 2014-01-24 2015-05-20 西安石油大学 一种煤层气储层密度测井扩径影响的校正方法
CN104778303B (zh) * 2015-02-12 2017-11-14 中国石油大学(北京) 一种裂缝型地层防塌钻井液性能参数的设计方法
CN105401939A (zh) * 2015-11-30 2016-03-16 中国石油大学(北京) 一种多因素耦合作用下的煤层井壁稳定性分析方法
CN105550448B (zh) * 2015-12-15 2018-08-14 中国石油天然气股份有限公司 基于钻井轨迹设计参数的预钻井三维井眼建模方法及装置

Also Published As

Publication number Publication date
CN106682384A (zh) 2017-05-17

Similar Documents

Publication Publication Date Title
CN106682384B (zh) 井眼扩径率计算方法及装置
Liu et al. Some key technical issues in modelling of gas transport process in shales: a review
Trimonova et al. The Study of the Unstable Fracure Propagation in the Injection Well: Numerical and Laboratory Modeling
CN103278400B (zh) 土体原位环剪实验仪
CN106917623B (zh) 钻井井壁稳定性预测方法及装置
Uchida et al. Geomechanical study of the Mallik gas hydrate production field trials
CN105738252A (zh) 一种裂缝内稠油可流动开度界限的测量方法
CN102809522A (zh) 一种双重介质岩心的加工方法与双重介质岩心的应用
Chen et al. Simulation of damage–permeability coupling for mortar under dynamic loads
Abdulhadi et al. Experimental study of wellbore instability in clays
Su et al. Experiment and failure mechanism of cement sheath integrity under development and production conditions based on a mechanical equivalent theory
CN205982211U (zh) 用于测试钻井液与岩石间压力传递的实验装置
CN203201545U (zh) 一种固井水泥环完整性模拟评价试验仪
CN110096669B (zh) 一种获取碳酸盐岩油藏中裂缝体积的方法
Lian et al. Numerical simulation of cement-to-formation interface debonding during hydraulic fracturing of shale gas wells
CN106772674B (zh) 一种裂缝性砂岩储层含油饱和度的计算方法
Hu et al. Full‐Life‐Cycle Analysis of Cement Sheath Integrity
Opedal et al. Cement Bond Strength Measurements
CN115822532B (zh) 一种确定二氧化碳-原油最小混相压力的方法、装置、设备及介质
Wang et al. Study on annular pressure buildup phenomenon in subsea wells considering the effect of cement
Huang et al. Stress sensitivity analysis of fractal porous media based on the elasto-plastic thick-walled cylinder model
Li et al. Capillary pressure at the imbibition front during water–oil counter-current spontaneous imbibition
Masalmeh et al. Oil mobility in transition zone
Ladubec et al. An extended finite element method model for carbon sequestration
Wu et al. Thermal technique for controlling hydraulic fractures

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