CN109359420A - 不同工况下射孔对封隔器冲击压力预测方法及装置 - Google Patents

不同工况下射孔对封隔器冲击压力预测方法及装置 Download PDF

Info

Publication number
CN109359420A
CN109359420A CN201811351465.6A CN201811351465A CN109359420A CN 109359420 A CN109359420 A CN 109359420A CN 201811351465 A CN201811351465 A CN 201811351465A CN 109359420 A CN109359420 A CN 109359420A
Authority
CN
China
Prior art keywords
perforation
packer
surge
pressure
under
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
CN201811351465.6A
Other languages
English (en)
Other versions
CN109359420B (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 University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201811351465.6A priority Critical patent/CN109359420B/zh
Publication of CN109359420A publication Critical patent/CN109359420A/zh
Application granted granted Critical
Publication of CN109359420B publication Critical patent/CN109359420B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种不同工况下射孔对封隔器冲击压力预测方法及装置,该方法包括:获取不同射孔工况下的射孔数据;根据射孔数据确定不同射孔工况下的射孔爆炸载荷计算模型;获取不同射孔工况下的封隔器处峰值压力;将射孔爆炸载荷和封隔器处峰值压力进行拟合,确定射孔冲击波衰减模型;基于射孔爆炸载荷计算模型、射孔冲击波衰减模型和射孔冲击波对封隔器冲击压力变化方程,确定不同射孔工况下的封隔器所受冲击压力计算模型;获取实际射孔工况下的实际射孔数据,基于所述封隔器所受冲击压力计算模型,确定实际射孔工况下封隔器所受冲击压力。该方案可以更加准确地评估不同射孔工况下封隔器的安全性。

Description

不同工况下射孔对封隔器冲击压力预测方法及装置
技术领域
本发明涉及油气井工程射孔技术领域,特别涉及一种不同工况下射孔对封隔器冲击压力预测方法及装置。
背景技术
射孔作业的目的在于使井筒与油气层之间形成通路,是油气田开采的关键环节。近些年来,射孔测试联作广泛应用于油气井完井试油过程中,作业时射孔枪常常结合封隔器进行联合作业,其中,封隔器连接在管柱上,在封隔器完全坐封后,射孔枪起爆时产生巨大爆轰波,有部分爆轰波会向井筒狭长空间内释放形成动态冲击载荷,一方面直接作用在枪管上,并通过枪管将这部分载荷传递给与之相连的减震器、油管、筛管、封隔器等其他管柱结构,引起管柱系统强烈冲击振动;另一方面冲击载荷会造成管柱外环液体压力在短时间内发生剧烈的变化,以冲击波形式在射孔液中传播,瞬间造成井内液体大变形、高速剧烈运动,影响整个管柱系统的结构稳定性以及局部结构强度。而封隔器作为连接在管柱上的重要部件,在这种复杂的环境很容易自动解封,导致射孔工艺失败,直接影响了射孔作业的安全性。因此,有必要针对不同射孔工况下封隔器的安全性展开研究工作。早期射孔技术研究多集中在射孔参数方面,而针对射孔爆炸冲击问题研究相对较少,随着射孔强度和井深的增加,射孔工程问题越来越多,射孔爆轰压力及其对管柱及封隔器安全影响问题逐渐引起人们的重视,但针对不同射孔工况下封隔器的安全性尚无详细的研究。
发明内容
本发明实施例提供了一种不同工况下射孔对封隔器冲击压力预测方法及装置,可以评估不同射孔工况下封隔器的安全性。
本发明实施例提供了一种不同工况下射孔对封隔器冲击压力预测方法,该方法包括:
获取不同射孔工况下的射孔数据,所述射孔数据包括射孔弹数量、单发装药量、油管长度、地层压力和井筒压力;
根据所述不同射孔工况下的射孔数据,确定射孔爆炸载荷计算模型;
获取不同射孔工况下的封隔器处峰值压力;
将不同射孔工况下的射孔爆炸载荷和不同射孔工况下的封隔器处峰值压力进行拟合,确定射孔冲击波衰减模型;
基于射孔爆炸载荷计算模型、射孔冲击波衰减模型和射孔冲击波对封隔器冲击压力变化方程,确定封隔器所受冲击压力计算模型;
获取实际射孔工况下的实际射孔数据,基于所述封隔器所受冲击压力计算模型,确定实际射孔工况下封隔器所受冲击压力。
本发明实施例还提供了一种不同工况下射孔对封隔器冲击压力预测装置,该装置包括:
射孔数据获取模块,用于获取不同射孔工况下的射孔数据,所述射孔数据包括射孔弹数量、单发装药量、油管长度、地层压力和井筒压力;
射孔爆炸载荷计算模型确定模块,用于根据所述不同射孔工况下的射孔数据,确定射孔爆炸载荷计算模型;
封隔器处峰值压力获取模块,用于获取不同射孔工况下的封隔器处峰值压力;
射孔冲击波衰减模型确定模块,用于将不同射孔工况下的射孔爆炸载荷和不同射孔工况下的封隔器处峰值压力进行拟合,确定射孔冲击波衰减模型;
封隔器所受冲击压力计算模型确定模块,用于基于射孔爆炸载荷计算模型、射孔冲击波衰减模型和射孔冲击波对封隔器冲击压力变化方程,确定封隔器所受冲击压力计算模型;
封隔器所受冲击压力计算模块,用于获取实际射孔工况下的实际射孔数据,基于所述封隔器所受冲击压力计算模型,确定实际射孔工况下封隔器所受冲击压力。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述所述不同工况下射孔对封隔器冲击压力预测方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述所述不同工况下射孔对封隔器冲击压力预测方法的计算机程序。
在本发明实施例中,从压力变化的角度,根据冲击波衰减规律,从理论上提出了预测不同工况下射孔对封隔器冲击压力大小的方法,最终建立了可应用于现场评估不同射孔工况下封隔器安全性的计算模型。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种不同工况下射孔对封隔器冲击压力预测方法流程图;
图2是本发明实施例提供的一种射孔管柱系统模型示意图;
图3是本发明实施例提供的一种射孔管柱系统模型网格划分示意图;
图4是本发明实施例提供的一种射孔冲击波在封隔器不同面的反射和折射示意图;
图5是本发明实施例提供的一种不同封隔器设置距离下压力随时间变化示意图;
图6是本发明实施例提供的一种不同工况下射孔对封隔器冲击压力预测装置框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明实施例中,提供了一种不同工况下射孔对封隔器冲击压力预测方法,如图1所示,该方法包括:
步骤101:获取不同射孔工况下的射孔数据,所述射孔数据包括射孔弹数量、单发装药量、油管长度、地层压力和井筒压力;
步骤102:根据所述不同射孔工况下的射孔数据,确定射孔爆炸载荷计算模型;
步骤103:获取不同射孔工况下的封隔器处峰值压力;
步骤104:将不同射孔工况下的射孔爆炸载荷和不同射孔工况下的封隔器处峰值压力进行拟合,确定射孔冲击波衰减模型;
步骤105:基于射孔爆炸载荷计算模型、射孔冲击波衰减模型和射孔冲击波对封隔器冲击压力变化方程,确定封隔器所受冲击压力计算模型;
步骤106:获取实际射孔工况下的实际射孔数据,基于所述封隔器所受冲击压力计算模型,确定实际射孔工况下封隔器所受冲击压力。
在本发明实施例中,在执行步骤101之前,还需要基于实际射孔情况,针对射孔管柱系统进行合理的简化。
(1)简化原因:在实际射孔作业中,不同的井况条件下射孔管柱长度不同,可能数十米到上千米不等,射孔管柱系统上有着不同的组件(比如,点火头、接头、减震器、油管、筛管、封隔器等),在对射孔管柱系统进行数值分析时必须进行简化。
(2)简化依据:根据现场射孔工艺及配套工具的规范,射孔作业过程中对整个封隔器以下管柱系统动态响应过程进行分析并提取关键环节作为模型。
由于点火头、接头和减震器等壁厚较大,屈服强度较高,简化时可以忽略,而油管柱及封隔器中心杆等屈服强度要低些,最容易出现整体屈曲和断裂,需要重点考虑,不能忽略;
(3)简化模型:在不影响仿真结果的前提下,对管柱系统结构进行合理的简化(包括射孔弹、射孔枪、筛管、减震器、油管、封隔器、套管等),模型简化后主要包括射孔枪、油管和套管,射孔枪内除弹药外剩余空间充满空气,油管内、环空内则充满射孔液。井筒内射孔管柱上端受到封隔器约束,周围受到套管的限制。
另外,还要对地层条件进行简化,简化为只考虑地层压力。
在本发明实施例中,步骤102按照如下方式实现:
(1)基于上述简化的测试管柱系统,在建立多套有限元仿真模型时,对不同参数的选择采用控制变量法,即针对某个参数研究时,其他参数保持不变。以钢级N80油管进行建模,屈服强度为552MPa,杨氏模量206GPa,剪切模量79.4GPa,泊松比0.3,密度7846kg/m3;变量参数主要为:射孔枪型号、射孔弹数量、单发装药量、油管长度、地层压力和井筒压力等,其中5"(5寸)射孔枪模型中,射孔弹数量120、150、180、210、234变化,单发装药量25、30、37、45g变化,油管长度6、10、14、20米变化,井筒压力58、54、50、46、42MPa,地层压力40、50、60、70MPa,射孔枪长度、井底口袋长度、射孔相位等参数保持不变;7"(7寸)射孔枪模型中,射孔弹数量162、198、234、276变化,单发装药量37、45、53、61g变化,油管长度6、10、14、20米变化,井筒压力58、54、50、46、42MPa,地层压力40、50、60、70MPa,射孔枪长度、井底口袋长度、射孔相位等参数保持不变。有限元仿真的射孔管柱系统模型示意图如图2所示。
(2)采用HYPERMESH对有限元仿真的射孔管柱系统模型进行了网格划分,如图3所示。为了节省仿真计算时间,在保证计算精度的前提下,增大了网格尺度,使网格平均间距在4~5mm之间。同时,为保证流体物质流动及能量在各部分网格之间有效传递,则必须保证各部分材料在相接界面上的网格共节点,并且又要保证各部分材料结构自身的界面形状,整个计算模型共1162724个节点、1050895个单元。最终完整的仿真模型是以k文件的形式输入到LS-DYNA程序中进行计算的。
(3)针对上述建立的多套仿真模型,利用ANSYS/LS-DYNA模块展开大量数值模拟计算,而后利用高级有限元前后处理软件LS-PREPOST提取不同射孔弹数量、单发装药量、油管长度、地层压力及井筒压力的射孔爆炸载荷数据并建立相应的数据库,利用MATLAB软件针对建立的数据库进行多元非线性回归,拟合得到不同射孔工况下爆炸载荷输出大小经验公式(即射孔爆炸载荷计算模型):
其中,PL为射孔爆炸载荷峰值压力(未经过反射和投射的),MPa;x1为射孔弹数量,枚;x2为单发装药量,g;x3为地层压力,MPa;x4为井筒压力,MPa;x5为油管长度,m;k、a、b、c为拟合相关系数。其中,针对5in(英寸)及7in射孔枪,k、a、b、c的数值如下表1所示:
表1
k a b c
5in 1.82 0.185 0.13 0.14
7in 1.87 0.212 0.134 0.165
上述数学模型仅作为预测方法的一种举例,本方法不仅限于该表达式。因此该数学模型不应视为本发明的限制。
在本发明实施例中,步骤103按照如下方式实现:
(1)针对某一射孔工况,设置封隔器到射孔枪顶部的距离分别为10,15…30…55,60米;
(2)基于数值模拟结果,利用高级有限元前后处理软件LS-PREPOST提取该射孔工况下压力随时间变化数据,如图4所示。从图4中可以得到压力峰值点,进而得到封隔器处峰值压力;
(3)然后针对不同的射孔工况,同样设置封隔器到射孔枪顶部的距离分别为10,15…30…55,60米,然后基于数值模拟结果,利用高级有限元前后处理软件LS-PREPOST提取不同射孔工况压力随时间变化数据,然后找出峰值点,进而得到不同射孔工况下的封隔器处峰值压力(有多个数值)。
在本发明实施例中,步骤104按照如下方式实现:
(1)基于水下爆炸试验,当装药在水介质中发生爆炸时,装药本身体积内形成高温高压的爆轰产物,其压力远大于周围介质的静压力,产生水中冲击波和气泡脉动。在通常情况下,水中的冲击波和反射波的传播可以近似看作符合声学理论的规律,即冲击波传播的衰减服从于指数衰减。
(2)借鉴水下爆炸模型,将不同射孔工况下的射孔爆炸载荷和不同射孔工况下的封隔器处峰值压力进行拟合,可以拟合得到射孔冲击波衰减公式(即射孔冲击波衰减模型):
PR=PL×e-αR
其中,PR为封隔器处峰值压力,MPa;PL为射孔爆炸载荷峰值压力,MPa;α为拟合所得的衰减系数;R是封隔器安全设置距离,m。
在本发明实施例中,步骤105按照如下方式实现:
(1)由于封隔器置于液体介质中,因此存在两个界面(如图5所示,接触面Ⅰ和接触面Ⅱ),而射孔冲击波在井筒液体中传播时,遇到封隔器会发生反射和透射,封隔器所受冲击压力会发生变化。
(2)根据射孔冲击波反射和透射定律,可以确定射孔冲击波对封隔器冲击压力变化方程:
其中,P为反射及透射后射孔冲击波对封隔器冲击压力;PF为反射压力;PT为透射压力;(ρc)p为常态下橡胶制封隔器的冲击阻抗;(ρc)f为常态下水介质的冲击阻抗;PR为封隔器处峰值压力,MPa。由此可得,封隔器对压力的反射会使压力增加。
(3)基于射孔爆炸载荷计算模型、射孔冲击波衰减模型和射孔冲击波对封隔器冲击压力变化方程,确定封隔器所受冲击压力计算模型:
其中,PR为封隔器处峰值压力,MPa;x1为射孔弹数量,枚;x2为单发装药量,g;x3为地层压力,MPa;x4为井筒压力,MPa;x5为油管长度,m;R是封隔器安全设置距离,m;A、a、b、c为拟合相关系数;α为拟合所得的衰减系数。
结合现场实际射孔情况,根据上述建立的不同射孔工况下封隔器所受冲击压力计算模型,即可简单快速地计算预测不同射孔工况下封隔器所受射孔冲击力大小。
基于同一发明构思,本发明实施例中还提供了一种不同工况下射孔对封隔器冲击压力预测装置,如下面的实施例所述。由于不同工况下射孔对封隔器冲击压力预测装置解决问题的原理与不同工况下射孔对封隔器冲击压力预测方法相似,因此…装置的实施可以参见不同工况下射孔对封隔器冲击压力预测方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图6是本发明实施例的不同工况下射孔对封隔器冲击压力预测装置结构框图,如图6所示,包括:
射孔数据获取模块601,用于获取不同射孔工况下的射孔数据,所述射孔数据包括射孔弹数量、单发装药量、油管长度、地层压力和井筒压力;
射孔爆炸载荷计算模型确定模块602,用于根据所述不同射孔工况下的射孔数据,确定射孔爆炸载荷计算模型;
封隔器处峰值压力获取模块603,用于获取不同射孔工况下的封隔器处峰值压力;
射孔冲击波衰减模型确定模块604,用于将不同射孔工况下的射孔爆炸载荷和不同射孔工况下的封隔器处峰值压力进行拟合,确定射孔冲击波衰减模型;
封隔器所受冲击压力计算模型确定模块605,用于基于射孔爆炸载荷计算模型、射孔冲击波衰减模型和射孔冲击波对封隔器冲击压力变化方程,确定封隔器所受冲击压力计算模型;
封隔器所受冲击压力计算模块606,用于获取实际射孔工况下的实际射孔数据,基于所述封隔器所受冲击压力计算模型,确定实际射孔工况下封隔器所受冲击压力。
下面对该结构进行说明。
在本发明实施例中,所述射孔爆炸载荷计算模型如下:
其中,PL为射孔爆炸载荷峰值压力,MPa;x1为射孔弹数量,枚;x2为单发装药量,g;x3为地层压力,MPa;x4为井筒压力,MPa;x5为油管长度,m;k、a、b、c为拟合相关系数。
在本发明实施例中,所述射孔冲击波衰减模型如下:
PR=PL×e-αR
其中,PR为封隔器处峰值压力,MPa;PL为射孔爆炸载荷峰值压力,MPa;α为拟合所得的衰减系数;R是封隔器安全设置距离,m。
在本发明实施例中,所述射孔冲击波对封隔器冲击压力变化方程如下:
其中,P为反射及透射后射孔冲击波对封隔器冲击压力;PF为反射压力;PT为透射压力;(ρc)p为常态下橡胶制封隔器的冲击阻抗;(ρc)f为常态下水介质的冲击阻抗;PR为封隔器处峰值压力,MPa。
在本发明实施例中,所述封隔器所受冲击压力计算模型如下:
其中,PR为封隔器处峰值压力,MPa;x1为射孔弹数量,枚;x2为单发装药量,g;x3为地层压力,MPa;x4为井筒压力,MPa;x5为油管长度,m;R是封隔器安全设置距离,m;A、a、b、c为拟合相关系数;α为拟合所得的衰减系数。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述所述不同工况下射孔对封隔器冲击压力预测方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述所述不同工况下射孔对封隔器冲击压力预测方法的计算机程序。
综上所述,本发明提出的不同工况下射孔对封隔器冲击压力预测方法及装置,从压力变化的角度,根据冲击波衰减规律,从理论上提出了预测不同工况下射孔对封隔器冲击压力大小的方法,最终建立了可应用于现场评估不同射孔工况下封隔器安全性的计算模型。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种不同工况下射孔对封隔器冲击压力预测方法,其特征在于,包括:
获取不同射孔工况下的射孔数据,所述射孔数据包括射孔弹数量、单发装药量、油管长度、地层压力和井筒压力;
根据所述不同射孔工况下的射孔数据,确定射孔爆炸载荷计算模型;
获取不同射孔工况下的封隔器处峰值压力;
将不同射孔工况下的射孔爆炸载荷和不同射孔工况下的封隔器处峰值压力进行拟合,确定射孔冲击波衰减模型;
基于射孔爆炸载荷计算模型、射孔冲击波衰减模型和射孔冲击波对封隔器冲击压力变化方程,确定封隔器所受冲击压力计算模型;
获取实际射孔工况下的实际射孔数据,基于所述封隔器所受冲击压力计算模型,确定实际射孔工况下封隔器所受冲击压力。
2.如权利要求1所述的不同工况下射孔对封隔器冲击压力预测方法,其特征在于,所述射孔爆炸载荷计算模型如下:
其中,PL为射孔爆炸载荷峰值压力,MPa;x1为射孔弹数量,枚;x2为单发装药量,g;x3为地层压力,MPa;x4为井筒压力,MPa;x5为油管长度,m;k、a、b、c为拟合相关系数。
3.如权利要求1所述的不同工况下射孔对封隔器冲击压力预测方法,其特征在于,所述射孔冲击波衰减模型如下:
PR=PL×e-αR
其中,PR为封隔器处峰值压力,MPa;PL为射孔爆炸载荷峰值压力,MPa;α为拟合所得的衰减系数;R是封隔器安全设置距离,m。
4.如权利要求1所述的不同工况下射孔对封隔器冲击压力预测方法,其特征在于,所述射孔冲击波对封隔器冲击压力变化方程如下:
其中,P为反射及透射后射孔冲击波对封隔器冲击压力;PF为反射压力;PT为透射压力;(ρc)p为常态下橡胶制封隔器的冲击阻抗;(ρc)f为常态下水介质的冲击阻抗;PR为封隔器处峰值压力,MPa。
5.如权利要求1所述的不同工况下射孔对封隔器冲击压力预测方法,其特征在于,所述封隔器所受冲击压力计算模型如下:
其中,PR为封隔器处峰值压力,MPa;x1为射孔弹数量,枚;x2为单发装药量,g;x3为地层压力,MPa;x4为井筒压力,MPa;x5为油管长度,m;R是封隔器安全设置距离,m;A、a、b、c为拟合相关系数;α为拟合所得的衰减系数。
6.一种不同工况下射孔对封隔器冲击压力预测装置,其特征在于,包括:
射孔数据获取模块,用于获取不同射孔工况下的射孔数据,所述射孔数据包括射孔弹数量、单发装药量、油管长度、地层压力和井筒压力;
射孔爆炸载荷计算模型确定模块,用于根据所述不同射孔工况下的射孔数据,确定射孔爆炸载荷计算模型;
封隔器处峰值压力获取模块,用于获取不同射孔工况下的封隔器处峰值压力;
射孔冲击波衰减模型确定模块,用于将不同射孔工况下的射孔爆炸载荷和不同射孔工况下的封隔器处峰值压力进行拟合,确定射孔冲击波衰减模型;
封隔器所受冲击压力计算模型确定模块,用于基于射孔爆炸载荷计算模型、射孔冲击波衰减模型和射孔冲击波对封隔器冲击压力变化方程,确定封隔器所受冲击压力计算模型;
封隔器所受冲击压力计算模块,用于获取实际射孔工况下的实际射孔数据,基于所述封隔器所受冲击压力计算模型,确定实际射孔工况下封隔器所受冲击压力。
7.如权利要求6所述的不同工况下射孔对封隔器冲击压力预测装置,其特征在于,所述射孔爆炸载荷计算模型如下:
其中,PL为射孔爆炸载荷峰值压力,MPa;x1为射孔弹数量,枚;x2为单发装药量,g;x3为地层压力,MPa;x4为井筒压力,MPa;x5为油管长度,m;k、a、b、c为拟合相关系数。
8.如权利要求6所述的不同工况下射孔对封隔器冲击压力预测装置,其特征在于,所述射孔冲击波衰减模型如下:
PR=PL×e-αR
其中,PR为封隔器处峰值压力,MPa;PL为射孔爆炸载荷峰值压力,MPa;α为拟合所得的衰减系数;R是封隔器安全设置距离,m。
9.如权利要求6所述的不同工况下射孔对封隔器冲击压力预测装置,其特征在于,所述射孔冲击波对封隔器冲击压力变化方程如下:
其中,P为反射及透射后射孔冲击波对封隔器冲击压力;PF为反射压力;PT为透射压力;(ρc)p为常态下橡胶制封隔器的冲击阻抗;(ρc)f为常态下水介质的冲击阻抗;PR为封隔器处峰值压力,MPa。
10.如权利要求6所述的不同工况下射孔对封隔器冲击压力预测装置,其特征在于,所述封隔器所受冲击压力计算模型如下:
其中,PR为封隔器处峰值压力,MPa;x1为射孔弹数量,枚;x2为单发装药量,g;x3为地层压力,MPa;x4为井筒压力,MPa;x5为油管长度,m;R是封隔器安全设置距离,m;A、a、b、c为拟合相关系数;α为拟合所得的衰减系数。
CN201811351465.6A 2018-11-14 2018-11-14 不同工况下射孔对封隔器冲击压力预测方法及装置 Active CN109359420B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811351465.6A CN109359420B (zh) 2018-11-14 2018-11-14 不同工况下射孔对封隔器冲击压力预测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811351465.6A CN109359420B (zh) 2018-11-14 2018-11-14 不同工况下射孔对封隔器冲击压力预测方法及装置

Publications (2)

Publication Number Publication Date
CN109359420A true CN109359420A (zh) 2019-02-19
CN109359420B CN109359420B (zh) 2021-02-02

Family

ID=65345140

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811351465.6A Active CN109359420B (zh) 2018-11-14 2018-11-14 不同工况下射孔对封隔器冲击压力预测方法及装置

Country Status (1)

Country Link
CN (1) CN109359420B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112926242A (zh) * 2021-02-04 2021-06-08 长江大学 一种油气井射孔动态冲击过程的模拟方法及装置
CN113218610A (zh) * 2021-04-08 2021-08-06 哈尔滨工程大学 一种强冲击大脉宽冲击载荷模拟装置及其控制方法
CN113255245A (zh) * 2021-05-21 2021-08-13 长江大学 一种井下射孔冲击波预测方法、装置及可读存储介质
CN113324448A (zh) * 2021-05-12 2021-08-31 西安近代化学研究所 一种弹药爆炸场风动压测试方法及装置
CN113742966A (zh) * 2021-08-24 2021-12-03 长江大学 一种确定射孔冲击载荷影响因素主次关系的方法
CN113761766A (zh) * 2021-08-20 2021-12-07 长江大学 一种预测射孔爆炸剩余能量的方法、装置及存储介质
CN115324538A (zh) * 2022-08-09 2022-11-11 西南石油大学 一种用于油气勘探的射孔管柱动力学系统与分析方法
WO2023124938A1 (zh) * 2021-12-31 2023-07-06 中国石油天然气集团有限公司 测试完井管柱的射孔爆轰分析方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040029739A1 (en) * 2001-01-18 2004-02-12 Hashem Mohamed Naguib Determining the pvt properties of a hydrocarbon reservoir fluid
CN105484687A (zh) * 2015-11-27 2016-04-13 中国海洋石油总公司 一种射孔瞬时井液压力的预测方法
CN107423471A (zh) * 2017-04-25 2017-12-01 西南石油大学 基于封隔器受力分析的井下射孔测试工具串优化方法
CN108625828A (zh) * 2018-03-28 2018-10-09 中国石油大学(北京) 预测射孔爆炸载荷输出大小的方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040029739A1 (en) * 2001-01-18 2004-02-12 Hashem Mohamed Naguib Determining the pvt properties of a hydrocarbon reservoir fluid
CN105484687A (zh) * 2015-11-27 2016-04-13 中国海洋石油总公司 一种射孔瞬时井液压力的预测方法
CN107423471A (zh) * 2017-04-25 2017-12-01 西南石油大学 基于封隔器受力分析的井下射孔测试工具串优化方法
CN108625828A (zh) * 2018-03-28 2018-10-09 中国石油大学(北京) 预测射孔爆炸载荷输出大小的方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHENG-BING LI等: "Numerical Computation of Wellbore Dynamic Pressure during Perforating", 《2012 FOURTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL AND INFORMATION SCIENCES》 *
李柏杨: "射孔冲击载荷下封隔器中心管强度安全性分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112926242A (zh) * 2021-02-04 2021-06-08 长江大学 一种油气井射孔动态冲击过程的模拟方法及装置
CN113218610A (zh) * 2021-04-08 2021-08-06 哈尔滨工程大学 一种强冲击大脉宽冲击载荷模拟装置及其控制方法
CN113218610B (zh) * 2021-04-08 2022-09-06 哈尔滨工程大学 一种强冲击大脉宽冲击载荷模拟装置及其控制方法
CN113324448A (zh) * 2021-05-12 2021-08-31 西安近代化学研究所 一种弹药爆炸场风动压测试方法及装置
CN113324448B (zh) * 2021-05-12 2022-08-19 西安近代化学研究所 一种弹药爆炸场风动压测试方法及装置
CN113255245A (zh) * 2021-05-21 2021-08-13 长江大学 一种井下射孔冲击波预测方法、装置及可读存储介质
CN113761766A (zh) * 2021-08-20 2021-12-07 长江大学 一种预测射孔爆炸剩余能量的方法、装置及存储介质
CN113742966A (zh) * 2021-08-24 2021-12-03 长江大学 一种确定射孔冲击载荷影响因素主次关系的方法
WO2023124938A1 (zh) * 2021-12-31 2023-07-06 中国石油天然气集团有限公司 测试完井管柱的射孔爆轰分析方法及装置
CN115324538A (zh) * 2022-08-09 2022-11-11 西南石油大学 一种用于油气勘探的射孔管柱动力学系统与分析方法

Also Published As

Publication number Publication date
CN109359420B (zh) 2021-02-02

Similar Documents

Publication Publication Date Title
CN109359420A (zh) 不同工况下射孔对封隔器冲击压力预测方法及装置
Lu et al. An equivalent method for blasting vibration simulation
Ge et al. Experimental and numerical investigation of the centrifugal model for underwater explosion shock wave and bubble pulsation
CN101501298B (zh) 分析和控制在钻孔中水锤产生的波传播的方法和装置
Yuan et al. Experimental and numerical study on the effect of water-decoupling charge structure on the attenuation of blasting stress
CN108625828B (zh) 预测射孔爆炸载荷输出大小的方法及装置
Lu et al. Experimental and numerical investigation of the effect of blast-generated free surfaces on blasting vibration
CN103559383B (zh) 核电站防波堤抗强震性能预测与评价方法
Gannon Simulation of underwater explosions in close-proximity to a submerged cylinder and a free-surface or rigid boundary
CN105484687B (zh) 一种射孔瞬时井液压力的预测方法
CN103869363B (zh) 微地震定位方法及装置
Liu et al. Pressure field investigation into oil&gas wellbore during perforating shaped charge explosion
CN112926242A (zh) 一种油气井射孔动态冲击过程的模拟方法及装置
Li et al. Application of coupled analysis methods for prediction of blast-induced dominant vibration frequency
CN114427445A (zh) 缝洞型油藏无限大地层动态产能计算方法和系统
Deng et al. Effects of perforation fluid movement on downhole packer with shock loads
Gu et al. Experimental and theoretical study on influence of different charging structures on blasting vibration energy
Zhao et al. Effects of working angle on pneumatic down-the-hole hammer drilling
CN113255245A (zh) 一种井下射孔冲击波预测方法、装置及可读存储介质
Xie et al. Large-scale field experiments on blast-induced vibration and crater in sand medium
Fan et al. Influence of water depth on the range of crushed zones and cracked zones for underwater rock drilling and blasting
Jianhua et al. Discussion on blasting vibration monitoring for rock damage control in rock slope excavation
CN113742966A (zh) 一种确定射孔冲击载荷影响因素主次关系的方法
Sun et al. Numerical investigation of fracturing of coal with cleats caused by supercritical CO2 explosion
Zhou et al. The effect of blast-induced vibration on the stability of underground water-sealed gas storage caverns

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