CN116029167B - 一种水平井多簇压裂非均匀射孔的定量设计方法 - Google Patents

一种水平井多簇压裂非均匀射孔的定量设计方法 Download PDF

Info

Publication number
CN116029167B
CN116029167B CN202310337706.6A CN202310337706A CN116029167B CN 116029167 B CN116029167 B CN 116029167B CN 202310337706 A CN202310337706 A CN 202310337706A CN 116029167 B CN116029167 B CN 116029167B
Authority
CN
China
Prior art keywords
perforation
fracturing
parameters
stress
fracture
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
CN202310337706.6A
Other languages
English (en)
Other versions
CN116029167A (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 East China
Original Assignee
China University of Petroleum East China
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 East China filed Critical China University of Petroleum East China
Priority to CN202310337706.6A priority Critical patent/CN116029167B/zh
Publication of CN116029167A publication Critical patent/CN116029167A/zh
Application granted granted Critical
Publication of CN116029167B publication Critical patent/CN116029167B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Abstract

本发明公开一种水平井多簇压裂非均匀射孔的定量设计方法，属于油气田开发技术领域，该水平井多簇压裂非均匀射孔的定量设计方法包括如下步骤：获取压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数；根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，确定各条裂缝附近的缝间应力干扰大小，根据所述各条裂缝附近的缝间应力干扰大小形成无量纲参数λ；在不考虑段内应力非均质性的情况下，选取第一射孔参数，使步骤S20所述的无量纲参数λ的最大值小于1；获取平均进液流量，根据所述平均进液流量，计算所述第一射孔参数对应的第一射孔摩阻以及不同射孔孔数、不同射孔直径下的第二射孔摩阻。

Description

一种水平井多簇压裂非均匀射孔的定量设计方法
技术领域
本发明属于油气田开发技术领域，具体涉及一种水平井多簇压裂非均匀射孔的定量设计方法。
背景技术
水平井分段多簇压裂是致密油气、非常规储层开发的关键技术。受缝间应力干扰、段内地应力非均质性的影响，多裂缝不能均衡扩展，导致部分射孔簇缺乏有效改造，降低了射孔簇效率，制约产量上升。
为了促进多裂缝均衡扩展，通常采用极限限流压裂、非均匀射孔技术，然而现有技术中通常采用数值模拟或经验的方法来选取极限限流压裂、非均匀射孔参数，对极限限流压裂、非均匀射孔压裂过程中的多裂缝均衡扩展、各射孔簇进液情况缺乏定量分析。即对于存在段内地应力非均质性的情况下，设计非均匀射孔参数缺乏准确、高效、可靠的设计方法。
发明内容
因此，本发明所要解决的是如何提供一种水平井多簇压裂非均匀射孔的定量设计方法，旨在解决现有技术中采用数值模拟或经验的方法来选取极限限流压裂、非均匀射孔参数，对极限限流压裂、非均匀射孔压裂过程中的多裂缝均衡扩展、各射孔簇进液情况缺乏定量分析的问题。
为了实现上述目的，本发明提供一种考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法，所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法包括如下步骤：
步骤S10：获取压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数；
步骤S20：根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，确定各条裂缝附近的缝间应力干扰大小，根据所述各条裂缝附近的缝间应力干扰大小形成无量纲参数λ；
步骤S30：在不考虑段内应力非均质性的情况下，选取第一射孔参数，使步骤S20所述的无量纲参数λ的最大值小于1；
步骤S40：获取平均进液流量，根据所述平均进液流量，计算所述第一射孔参数对应的第一射孔摩阻以及不同射孔孔数、不同射孔直径下的第二射孔摩阻；根据所述第一射孔摩阻和所述第二射孔摩阻，计算射孔摩阻降低值；
步骤S50：根据段内应力非均质性的大小和位置、以及所述射孔摩阻降低值，选取第二射孔参数使相应裂缝的射孔摩阻降低值接近或等于段内应力非均质性的大小。
优选地，在所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法中，所述步骤S10中所述压裂井井筒参数包括井斜角变化、套管内径、套管内壁粗糙度、井深、以及压裂段段长；
所述压裂完井信息包括簇间距、簇数、射孔参数；
所述储层地质参数包括段内地应力非均质性的分布、岩石杨氏模量、泊松比、断裂韧性、以及滤失系数；
所述压裂施工参数包括施工排量、液体粘度。
优选地，在所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法中，所述步骤S10包括：
步骤S11：获取压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数；
步骤S12：根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，建立裂缝扩展的物理模型，包括裂缝几何模型、油藏地质模型和井筒几何模型。
优选地，在所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法中，所述步骤S20包括：
根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，确定各条裂缝附近的缝间应力干扰大小，根据所述各条裂缝附近的缝间应力干扰大小除以对应裂缝的射孔摩阻形成无量纲参数λ。
优选地，在所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法中，所述步骤S20中确定各条裂缝附近的缝间应力干扰大小的步骤包括：
在所述储层地质参数、压裂施工参数下，采用粘度主导区PKN裂缝解析解计算缝内平均净压力，计算公式如式（1）：
Figure SMS_1
（1）
将计算得到的缝内平均净压力乘以衰减系数来描述缝间应力干扰随距离的变化情况，并对每条裂缝产生的缝间应力干扰进行叠加，得到每条裂缝受到的缝间应力干扰作用的计算式：
Figure SMS_2
（2）
其中，p av,k 为第k裂缝的缝内净压力，Pa；
μ为压裂液粘度，Pa·s；
E’为平面应变杨氏模量，Pa；
Q in,k 为第k裂缝的进液流量，m3/s；
t为时间，s；
H为储层厚度，m；
ξ为归一化裂缝长度，无量纲；
σ int,k 为第k裂缝受到的缝间应力干扰大小，Pa；
n f为裂缝条数；
d ik 为第i裂缝到第k裂缝的距离，m；
δ ik 为狄拉克函数，无量纲。
优选地，在所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法中，所述步骤S20中根据所述各条裂缝附近的缝间应力干扰大小形成无量纲参数λ，包括：
计算射孔摩阻，计算式如下：
Figure SMS_3
（3）
无量纲参数λ计算式为：
Figure SMS_4
（4）
其中，
p perf,k 为第k射孔簇的射孔摩阻，Pa；
ρ为压裂液的密度，kg/m3；
n k 为第k射孔簇的射孔数目；
d p,k 为第k射孔簇的射孔孔径，m；
K为射孔磨蚀修正系数；
Q in,k 为第k裂缝的进液流量，m3/s；
t为时间，s。
优选地，在所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法中，所述步骤S30包括：
步骤S31：在不考虑段内应力非均质性的情况下，根据平均进液流量，计算不同的射孔参数对应的无量纲参数λ的最大值；
步骤S32：根据步骤S31中计算得到的不同的射孔参数对应的无量纲参数λ的最大值，选取第一射孔参数，使步骤S20所述的无量纲参数λ的最大值小于1，实现段内地应力均质情况下的多裂缝均衡扩展。
优选地，在所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法中，所述步骤S40包括：
步骤S41：在各裂缝均衡进液的条件下，根据平均进液流量，计算所述第一射孔参数对应的第一射孔摩阻；
步骤S42：根据所述根据平均进液流量，利用计算不同射孔孔数、不同射孔直径下的第二射孔摩阻；
步骤S43：所述第一射孔摩阻和所述第二射孔摩阻的差值作为射孔摩阻降低值。
优选地，在所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法中，所述步骤S50包括：
根据段内应力非均质性的大小和位置，利用步骤S40计算得到的所述射孔摩阻降低值，在相应位置的裂缝调整得到第二射孔参数，使相应位置处裂缝的射孔摩阻降低值接近或等于段内应力非均质性的大小。
优选地，在所述考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法中，所述相应位置的裂缝为段内地应力非均质性所在的裂缝，对于存在段内应力非均质的情况下，利用该方法无需进行数值模拟即可取得合适的非均匀射孔参数设计方案。
本发明具有如下有益效果：
本发明通过获取压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数；再根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，确定各条裂缝附近的缝间应力干扰大小，根据所述各条裂缝附近的缝间应力干扰大小形成无量纲参数λ；再在不考虑段内应力非均质性的情况下，选取第一射孔参数，使步骤S20所述的无量纲参数λ的最大值小于1；再获取平均进液流量，根据所述平均进液流量，计算所述第一射孔参数对应的第一射孔摩阻以及不同射孔孔数、不同射孔直径下的第二射孔摩阻；根据所述第一射孔摩阻和所述第二射孔摩阻，计算射孔摩阻降低值；最后根据段内应力非均质性的大小和位置、以及所述射孔摩阻降低值，选取第二射孔参数使相应裂缝的射孔摩阻降低值接近或等于段内应力非均质性的大小，如此通过降低段内地应力非均质所在裂缝的射孔摩阻，促使多裂缝均衡扩展，对极限限流压裂、非均匀射孔压裂过程中的多裂缝均衡扩展、各射孔簇进液情况可以定量分析，解决了现有技术中对极限限流压裂、非均匀射孔压裂过程中的多裂缝均衡扩展、各射孔簇进液情况缺乏定量分析的问题；
对于在考虑段内应力非均质性的情况下，若采用第一射孔参数，则段内地应力非均质所在裂缝由于受到段内地应力非均质的影响，难以扩展，导致多裂缝扩展极度不均衡；而本发明在不考虑段内应力非均质性的情况下，选取第一射孔参数、并计算相应的射孔摩阻。要使多裂缝扩展均衡，则需要进一步计算不同射孔孔数、不同射孔直径下的第二射孔摩阻，选取适当的第二射孔参数，使得第二射孔摩阻的数值比第一射孔摩阻的数值低，并且第一射孔摩阻与第二射孔摩阻之差应当与段内地应力非均质的数值相等或接近；
进一步地，对于段内地应力非均质所在裂缝而言，相较于不考虑段内应力非均质性的情况，其对应的地应力升高但射孔摩阻降低，一升一降，并且其地应力升高值接近或等于射孔摩阻降低值，保证了与不考虑段内应力非均质性的情况相当，从而促使多裂缝均衡扩展。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为水平井多簇压裂物理模型示意图；
图2为不同射孔孔数、射孔孔径对应的无量纲参数λ最大值曲线图；
图3为不同射孔孔数、射孔孔径的射孔摩阻降低值；
图4为X井均匀射孔下的多裂缝扩展形态；
图5为X井非均匀射孔下的多裂缝扩展形态。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
本发明实施例中术语“和/或”，描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。
本发明实施例中术语“多个”是指两个或两个以上，其它量词与之类似。
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的各实施例进行详细的阐述。然而，本领域的普通技术人员可以理解，在本发明各实施例中，为了使读者更好地理解本发明而提出了许多技术细节。但是，即使没有这些技术细节和基于以下各实施例的种种变化和修改，也可以实现本发明所要求保护的技术方案。以下各个实施例的划分是为了描述方便，不应对本发明的具体实现方式构成任何限定，各个实施例在不矛盾的前提下可以相互结合相互引用。
本发明提供一种考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法，请参阅图1至图5，该考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法包括如下步骤：
步骤S10：获取压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数。
其中，所述压裂井井筒参数包括井斜角变化、套管内径、套管内壁粗糙度、井深、以及压裂段段长；
所述压裂完井信息包括簇间距、簇数、射孔参数（例如每簇射孔数、射孔直径等）；
所述储层地质参数包括段内地应力非均质性的分布、岩石杨氏模量、泊松比、断裂韧性、以及滤失系数；在本实施例中，储层地质参数（油藏地质参数）可以是通过测井技术获取。
所述压裂施工参数包括施工排量、液体粘度。
具体地，所述步骤S10包括：
步骤S11：获取压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数；
步骤S12：根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，建立裂缝扩展的物理模型，包括裂缝几何模型、油藏地质模型和井筒几何模型。
下面以X页岩气井为例，具体油藏地质、工程参数、井筒参数如表1所示。
表1 相关参数信息
参数 数值
杨氏模量, GPa 35
泊松比, 无量纲 0.25
断裂韧性,MPa·m0.5 0.5
储层厚度, m 40
段内最小水平地应力剖面(跟端到趾端)， MPa 60-60-60-60-62
压裂液粘度， mPa·s 5
施工排量，m3/min 14
压裂簇数 5
每簇射孔数 待设计量
簇间距，m 10
孔眼直径， mm 待设计量
射孔磨蚀修正系数， 无因次 0.8
注入时间， min 90
步骤S12可以是根据表1中提供的地质和工程参数，建立X页岩气水平井多簇压裂的5簇压裂的物理模型。
步骤S20：根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，确定各条裂缝附近的缝间应力干扰大小，根据所述各条裂缝附近的缝间应力干扰大小形成无量纲参数λ；
具体地，所述步骤S20包括：
根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，确定各条裂缝附近的缝间应力干扰大小，根据所述各条裂缝附近的缝间应力干扰大小除以对应裂缝的射孔摩阻形成无量纲参数λ。
更具体地，缝间应力干扰的估算主要包括缝内平均净压力的计算、缝间应力干扰的衰减系数、多裂缝缝间应力干扰的叠加。故，所述步骤S20中确定各条裂缝附近的缝间应力干扰大小的步骤包括：
（1）计算缝内平均净压力
在所述储层地质参数、压裂施工参数下，采用粘度主导区PKN裂缝解析解计算缝内平均净压力，计算公式如式（1）：
Figure SMS_5
（1）
以表1中的参数为例进行计算，可以采用粘度主导区PKN裂缝解析解计算缝内平均净压力。
将计算得到的缝内平均净压力乘以衰减系数来描述缝间应力干扰随距离的变化情况，并对每条裂缝产生的缝间应力干扰进行叠加，得到每条裂缝受到的缝间应力干扰作用的计算式：
Figure SMS_6
（2）
其中，p av,k 为第k裂缝的缝内净压力，Pa；
μ为压裂液粘度，Pa·s；
E’为平面应变杨氏模量，Pa；
Q in,k 为第k裂缝的进液流量，m3/s；
t为时间，s；
H为储层厚度，m；
ξ为归一化裂缝长度，无量纲；
σ int,k 为第k裂缝受到的缝间应力干扰大小，Pa；
n f为裂缝条数；
d ik 为第i裂缝到第k裂缝的距离，m；
δ ik 为狄拉克函数，无量纲。
（2）计算射孔摩阻并建立无量纲参数λ
所述步骤S20中根据所述各条裂缝附近的缝间应力干扰大小形成无量纲参数λ，包括：
压裂液在流经射孔孔眼时会产生压降，即射孔摩阻，计算射孔摩阻，计算式如下：
Figure SMS_7
（3）
无量纲参数λ计算式为：
Figure SMS_8
（4）
其中，
p perf,k 为第k射孔簇的射孔摩阻，Pa；
ρ为压裂液的密度，kg/m3；
n k 为第k射孔簇的射孔数目；
d p,k 为第k射孔簇的射孔孔径，m；
K为射孔磨蚀修正系数；
Q in,k 为第k裂缝的进液流量，m3/s；
t为时间，s。
步骤S30：在不考虑段内应力非均质性的情况下，选取第一射孔参数，使步骤S20所述的无量纲参数λ的最大值小于1；
具体地，所述步骤S30包括：
步骤S31：在不考虑段内应力非均质性的情况下，根据平均进液流量，计算不同的射孔参数对应的无量纲参数λ的最大值；
在本实施例中，在不考虑段内地应力非均质性的情况下，即各条裂缝附近的最小水平地应力均为60MPa，根据平均进液流量，在不同的射孔孔数、射孔直径下利用（4）式计算所述无量纲参数λ的最大值。
步骤S32：根据步骤S31中计算得到的不同的射孔参数对应的无量纲参数λ的最大值，选取第一射孔参数，使步骤S20所述的无量纲参数λ的最大值小于1，实现段内地应力均质情况下的多裂缝均衡扩展。
在本实施例中，选取第一射孔参数（射孔孔数、射孔直径），使得所述的无量纲参数λ的最大值小于1，从而实现段内地应力均质情况下的多裂缝均匀扩展。
以图2为例，选取的方案可以为射孔孔数为7、射孔孔径为10mm，或射孔孔数为11、射孔孔径为8mm。在段内地应力均质的情况下，两种方案均能使多裂缝均衡扩展。在本具体实施方式中，选取射孔孔数为7、射孔孔径为10mm的方案。
需要说明的是，射孔参数为单个射孔簇的参数，包括：射孔孔数、射孔直径；所述平均进液流量为施工排量除以裂缝条数，即各裂缝均衡进液条件下的进液流量。
对于不存在段内应力非均质的情况，利用该方法无需进行数值模拟即可精准、快速、高效获得合适的限流压裂射孔参数设计方案。
步骤S40：获取平均进液流量，根据所述平均进液流量，计算所述第一射孔参数对应的第一射孔摩阻以及不同射孔孔数、不同射孔直径下的第二射孔摩阻；根据所述第一射孔摩阻和所述第二射孔摩阻，计算射孔摩阻降低值；
具体地，所述步骤S40包括：
步骤S41：在各裂缝均衡进液的条件下，根据平均进液流量，计算所述第一射孔参数对应的第一射孔摩阻；
以表1中的数据为例，根据表1中的施工排量，计算各裂缝均匀进液情况下的平均进液流量，根据（3）式计算步骤S30选取射孔参数的第一射孔摩阻为5.6MPa。
步骤S42：根据所述根据平均进液流量，利用计算不同射孔孔数、不同射孔直径下的第二射孔摩阻；
以表1中的数据为例，在各裂缝均匀进液条件下计算不同射孔参数下的第二射孔摩阻。
步骤S43：所述第一射孔摩阻和所述第二射孔摩阻的差值作为射孔摩阻降低值。
以表1中的数据为例，用第一射孔摩阻5.6MPa减去该第二射孔摩阻，得到射孔摩阻降低值，计算结果如图3。
步骤S50：根据段内应力非均质性的大小和位置、以及所述射孔摩阻降低值，选取第二射孔参数使相应裂缝的射孔摩阻降低值接近或等于段内应力非均质性的大小。
需要说明的是，选取第二射孔参数使相应裂缝的射孔摩阻降低值接近或等于段内应力非均质性的大小，可以为第二射孔参数使相应裂缝的射孔摩阻降低值等于段内应力非均质性的大小，也可以是第二射孔参数使相应裂缝的射孔摩阻降低值和段内应力非均质性的大小之差的绝对值小于一预设阈值，例如预设阈值0.2MPa。
具体地，所述步骤S50包括：
根据段内应力非均质性的大小和位置，利用步骤S40计算得到的所述射孔摩阻降低值，在相应位置的裂缝调整得到第二射孔参数，使相应位置处裂缝的射孔摩阻降低值接近或等于段内应力非均质性的大小。
所述相应位置的裂缝为段内地应力非均质性所在的裂缝，对于存在段内应力非均质的情况下，无需进行数值模拟即可取得合适的非均匀射孔参数设计方案。
以表1中的数据为例，根据S42和图3的计算结果，选取与段内地应力非均质性较为接近的射孔摩阻降低值对应的第二射孔参数（包括射孔孔数、射孔孔径），将其作为相应裂缝的第二射孔参数，在本实施例中，选取射孔数为7孔、射孔孔径为11mm，或射孔数为8孔、射孔孔径为10.5mm，或射孔数为9孔、射孔孔径为10mm、或射孔数目为10孔、射孔孔径为9.5mm。
根据步骤S50选取的射孔孔数、射孔孔径，设计段内非均匀射孔方案为：射孔孔径为10mm，射孔孔数为（跟端到趾端）：7-7-7-7-9孔。图4给出了均匀射孔（所有射孔簇的射孔数均为7孔）的裂缝形态，图5为非均匀射孔（S50设计方案）的裂缝形态，两图对比发现，该非均匀射孔定量设计方法能够促进多裂缝均匀扩展。
显然，上述所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下，可以做出其它不同形式的变化或变动，都应当属于本发明保护的范围。

Claims (5)

1.一种考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法，其特征在于，包括如下步骤：
步骤S10：获取压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数；
步骤S20：根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，确定各条裂缝附近的缝间应力干扰大小，根据所述各条裂缝附近的缝间应力干扰大小形成无量纲参数λ，所述步骤S20包括：
根据所述压裂井井筒参数、压裂完井信息、储层地质参数、以及压裂施工参数，确定各条裂缝附近的缝间应力干扰大小，根据所述各条裂缝附近的缝间应力干扰大小除以对应裂缝的射孔摩阻形成无量纲参数λ，所述步骤S20中确定各条裂缝附近的缝间应力干扰大小的步骤包括：
在所述储层地质参数、压裂施工参数下，采用粘度主导区PKN裂缝解析解计算缝内平均净压力，计算公式如式(1)：
Figure FDA0004261157170000011
将计算得到的缝内平均净压力乘以衰减系数来描述缝间应力干扰随距离的变化情况，并对每条裂缝产生的缝间应力干扰进行叠加，得到每条裂缝受到的缝间应力干扰作用的计算式：
Figure FDA0004261157170000012
其中，pav,k为第k裂缝的缝内净压力，Pa；
μ为压裂液粘度，Pa·s；
E’为平面应变杨氏模量，Pa；
Qin,k为第k裂缝的进液流量，m3/s；
t为时间，s；
H为储层厚度，m；
ξ为归一化裂缝长度，无量纲；
σint,k为第k裂缝受到的缝间应力干扰大小，Pa；
nf为裂缝条数；
dik为第i裂缝到第k裂缝的距离，m；
δik为狄拉克函数，无量纲；
步骤S30：在不考虑段内应力非均质性的情况下，选取第一射孔参数，使步骤S20所述的无量纲参数λ的最大值小于1，所述步骤S30包括：
步骤S31：在不考虑段内应力非均质性的情况下，根据平均进液流量，计算不同的射孔参数对应的无量纲参数λ的最大值；
步骤S32：根据步骤S31中计算得到的不同的射孔参数对应的无量纲参数λ的最大值，选取第一射孔参数，使步骤S20所述的无量纲参数λ的最大值小于1，实现段内地应力均质情况下的多裂缝均衡扩展；
步骤S40：获取平均进液流量，根据所述平均进液流量，计算所述第一射孔参数对应的第一射孔摩阻以及不同射孔孔数、不同射孔直径下的第二射孔摩阻；根据所述第一射孔摩阻和所述第二射孔摩阻，计算射孔摩阻降低值；
步骤S50：根据段内应力非均质性的大小和位置，利用步骤S40计算得到的所述射孔摩阻降低值，在相应位置的裂缝调整得到第二射孔参数，使相应位置处裂缝的射孔摩阻降低值等于段内应力非均质性的大小。
2.如权利要求1所述的考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法，其特征在于，所述步骤S10中所述压裂井井筒参数包括井斜角变化、套管内径、套管内壁粗糙度、井深、以及压裂段段长；
所述压裂完井信息包括簇间距、簇数、射孔参数；
所述储层地质参数包括段内地应力非均质性的分布、岩石杨氏模量、泊松比、断裂韧性、以及滤失系数；
所述压裂施工参数包括施工排量、液体粘度。
3.如权利要求1所述的考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法，其特征在于，所述步骤S20中根据所述各条裂缝附近的缝间应力干扰大小形成无量纲参数λ，包括：
计算射孔摩阻，计算式如下：
Figure FDA0004261157170000031
无量纲参数λ计算式为：
Figure FDA0004261157170000032
其中，
pperf,k为第k射孔簇的射孔摩阻，Pa；
ρ为压裂液的密度，kg/m3；
nk为第k射孔簇的射孔数目；
dp,k为第k射孔簇的射孔孔径，m；
K为射孔磨蚀修正系数；
Qin,k为第k裂缝的进液流量，m3/s；
t为时间，s。
4.如权利要求1所述的考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法，其特征在于，所述步骤S40包括：
步骤S41：在各裂缝均衡进液的条件下，根据平均进液流量，计算所述第一射孔参数对应的第一射孔摩阻；
步骤S42：根据所述根据平均进液流量，利用计算不同射孔孔数、不同射孔直径下的第二射孔摩阻；
步骤S43：所述第一射孔摩阻和所述第二射孔摩阻的差值作为射孔摩阻降低值。
5.如权利要求1所述的考虑段内地应力非均质性的水平井多簇压裂非均匀射孔的定量设计方法，其特征在于，所述相应位置的裂缝为段内地应力非均质性所在的裂缝，对于存在段内应力非均质的情况下，无需进行数值模拟即可取得合适的非均匀射孔参数设计方案。
CN202310337706.6A 2023-03-31 2023-03-31 一种水平井多簇压裂非均匀射孔的定量设计方法 Active CN116029167B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310337706.6A CN116029167B (zh) 2023-03-31 2023-03-31 一种水平井多簇压裂非均匀射孔的定量设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310337706.6A CN116029167B (zh) 2023-03-31 2023-03-31 一种水平井多簇压裂非均匀射孔的定量设计方法

Publications (2)

Publication Number Publication Date
CN116029167A CN116029167A (zh) 2023-04-28
CN116029167B true CN116029167B (zh) 2023-07-04

Family

ID=86072716

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310337706.6A Active CN116029167B (zh) 2023-03-31 2023-03-31 一种水平井多簇压裂非均匀射孔的定量设计方法

Country Status (1)

Country Link
CN (1) CN116029167B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109424348A (zh) * 2017-09-05 2019-03-05 中国石油化工股份有限公司 水平井多簇射孔压裂裂缝条数的判断方法及系统
CN113850029A (zh) * 2021-08-27 2021-12-28 长江大学 一种页岩气水平井密切割压裂射孔参数优化设计方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110056336B (zh) * 2019-05-31 2020-01-07 西南石油大学 一种页岩气缝网压裂施工压力曲线自动诊断方法
CN110222477B (zh) * 2019-07-08 2020-01-21 西南石油大学 维持水平井分段压裂裂缝均衡扩展的射孔参数优化方法
CN111734380B (zh) * 2020-07-30 2022-03-11 西南石油大学 一种水平井多级压裂裂缝扩展形态的快速预测方法
CN112836446B (zh) * 2021-01-28 2022-04-26 重庆科技学院 一种致密储层水平井分簇限流布孔优化设计方法
CN112949129B (zh) * 2021-03-03 2023-03-03 西南石油大学 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法
CN113255123B (zh) * 2021-05-17 2023-06-23 中国石油大学(华东) 煤层顶板水平井分段多簇压裂适用性地质条件的评价方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109424348A (zh) * 2017-09-05 2019-03-05 中国石油化工股份有限公司 水平井多簇射孔压裂裂缝条数的判断方法及系统
CN113850029A (zh) * 2021-08-27 2021-12-28 长江大学 一种页岩气水平井密切割压裂射孔参数优化设计方法

Also Published As

Publication number Publication date
CN116029167A (zh) 2023-04-28

Similar Documents

Publication Publication Date Title
CN110222477B (zh) 维持水平井分段压裂裂缝均衡扩展的射孔参数优化方法
US10294759B2 (en) Method and system for hydraulic fracturing based on skin factor analysis
CN113389534B (zh) 一种水平井密切割压裂裂缝扩展预测与设计参数优化方法
CN107288626B (zh) 一种页岩气地质甜度与工程甜度计算方法
CN109958416B (zh) 一种变孔径变孔密均匀进液进砂的多簇射孔压裂方法
CN113850029B (zh) 一种页岩气水平井密切割压裂射孔参数优化设计方法
CN109209333B (zh) 页岩气多井群高效开采间距优化方法
CN110439519B (zh) 一种基于极限限流设计的压裂方法及系统
CN111734380B (zh) 一种水平井多级压裂裂缝扩展形态的快速预测方法
CN110807541B (zh) 一种气藏开发井型的选择方法
CN109711090B (zh) 一种环空流体综合摩阻系数确定方法及装置
CN116029167B (zh) 一种水平井多簇压裂非均匀射孔的定量设计方法
CN114239308B (zh) 一种多尺度高密度压裂参数优化方法
US10837277B2 (en) Well completion system and method
CN108487905B (zh) 一种页岩气水平井压裂参数的优化方法
CN115310379A (zh) 井间干扰条件下压裂水平井的生产动态分析方法及设备
CN113530534B (zh) 储层改造方法及装置
CN113657056A (zh) 降低井底沉砂与提高射孔孔眼进砂效率的直井压裂加砂优化设计方法
CN111523696A (zh) 多层合采致密砂岩气藏的一种快速产能预测方法
WO2020122945A1 (en) System and method to optimize pumping
Rassenfoss So Many Wells, So Few Engineers: Scaling Production Engineering for all Those Shale Wells
US20090271161A1 (en) Arrangement of cutting elements on roller cones for earth boring bits
CN117113611A (zh) 深层页岩气水平井密切割压裂非均匀射孔优化设计方法
CN117350076A (zh) 一种深层页岩水平井多簇裂缝同步起裂的射孔优化方法
CN118036484A (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