CN112580271A - 水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法 - Google Patents

水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法 Download PDF

Info

Publication number
CN112580271A
CN112580271A CN202011449599.9A CN202011449599A CN112580271A CN 112580271 A CN112580271 A CN 112580271A CN 202011449599 A CN202011449599 A CN 202011449599A CN 112580271 A CN112580271 A CN 112580271A
Authority
CN
China
Prior art keywords
proppant
mpa
fracture
width
flow
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
CN202011449599.9A
Other languages
English (en)
Other versions
CN112580271B (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 CN202011449599.9A priority Critical patent/CN112580271B/zh
Publication of CN112580271A publication Critical patent/CN112580271A/zh
Application granted granted Critical
Publication of CN112580271B publication Critical patent/CN112580271B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Algebra (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明实施例涉及一种水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法，包括步骤S10，获取地层参数、施工参数、完井参数；步骤S20：设定初始的缝压和缝宽；步骤S30：建立耦合井筒流动的水平井分段多簇压裂裂缝扩展计算模型；步骤S40：建立时间步递增模型；步骤S50：建立压裂缝内支撑剂运移计算模型；步骤S60：判断该时刻是否加入支撑剂，若是，则采用高阶加权本质不震荡有限差分方法对压裂裂缝扩展过程中支撑剂运移进行计算，得到缝内支撑剂分布；若不是，则进入下一步；步骤S70：判断时间步是否达到施工时间，若没有则返回步骤S4，否则结束计算。

Description

水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模 方法
技术领域
本发明实施例涉及油气田开发技术领域，特别涉及一种水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法。
背景技术
非常规油气藏的开发依赖于水平井分段多簇压裂的发展和应用。水平井分段多簇压裂通过“分簇射孔”技术实现一段多裂缝扩展，从而进一步提高施工效率、节约施工成本。多簇压裂后在地层会形成具有一定导流能力的填砂裂缝，裂缝中的支撑剂分布形态将直接影响裂缝导流能力和油气井产能，故对压裂缝中的支撑剂运移规律进行研究十分必要。然而，目前对于水平井分段多簇压裂裂缝扩展、多裂缝非均匀扩展过程缝内支撑剂分布规律等问题研究并不充分，特别是水平井分段多簇压裂裂缝扩展与支撑剂运移一体化计算方法研究鲜有报道。
发明内容
本发明实施方式的目的在于提供一种水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法，旨在提供水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法。
为解决上述技术问题，本发明的实施方式提供了一种水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,包括：
步骤S10，获取地层参数、施工参数、完井参数；
步骤S20：设定初始的缝压和缝宽；
步骤S30：建立耦合井筒流动的水平井分段多簇压裂裂缝扩展计算模型；
步骤S40：建立时间递增模型；
步骤S50：建立压裂缝内支撑剂运移计算模型；
步骤S60：判断该时刻是否加入支撑剂，若是，则采用高阶加权本质不震荡有限差分方法对压裂裂缝扩展过程中支撑剂运移进行计算，得到缝内支撑剂分布；若不是，则进入下一步；
步骤S70：判断时间是否达到施工时间，若没有则返回步骤S4，否则结束计算。
优选地，所述步骤S10包括：
获取地应力分布、岩石力学参数、注入程序、加砂程序、液体粘度、排量、注入时间、分簇数量、裂缝间距、射孔数、射孔直径。
优选地，所述步骤S200的具体过程如下：
假设每个单元的宽度为常数，根据叠加原理，得到初始t时刻裂缝内压力与宽度的离散方程为
Figure BDA0002826268040000021
其中，N为单元总数；
p(x,t)为缝内流体压力，MPa；
σh(x)为最小主应力，MPa；
c(x,ξi)为平面三维裂缝核函数，MPa/m；
w(ξi)为裂缝宽度，m；
ξi为裂缝面网格点的坐标；
x为计算点的坐标；
通过求解式(1)可得到初始时刻的缝压与缝宽；
其中式(1)中的矩阵形式为
p-σh＝Cw (2)
其中，p为缝内压力向量，MPa；
σh为最小主应力向量，MPa；
C为影响系数矩阵，MPa/m；
w为缝宽向量。
优选地，所述步骤S300包括：
建立井筒流动模型；
建立裂缝宽度与压力的流固耦合方程；
建立裂缝动边界。
优选地，所述井筒流动模型的建立过程如下：
井筒到各簇裂缝的分流满足质量守恒和压力连续条件，式(3)和式(4)构成井筒流动模型，式(3)和式(4)如下：
Figure BDA0002826268040000031
pw＝pp,k+pc,k+pin,k (4)
式(2)中
Figure BDA0002826268040000032
式(2)中
Figure BDA0002826268040000033
其中，Qt为一段开启Nf条裂缝的情况，注入总流量，Qi为第i条分支裂缝流量；
pw为井底压力,MPa；
k＝1，2，......，Nf；
pp,k为k裂缝的射孔摩阻，MPa；
pc,k为井口到k裂缝的井筒流动摩阻，MPa；
pin,k为k裂缝的入口压力，MPa；
Qk为k裂缝的入口流量，m3/s；
ρ为液体密度，kg/m3；
nk为k射孔簇的射孔数量；
dk为k射孔簇的射孔直径，mm；
K为射孔磨蚀修正系数，无因次；
fc为沿程摩阻系数，无因次；
Dw为压裂管柱内径，m；lk为井口到k裂缝的管柱长度，m；
ε为压裂管柱的内壁粗糙度，m；
Vw为井筒内液体流速，m/s；
Re为雷诺数，Re＝DwρVw/μ；
μ为液体粘度，mPa·s；
ρ为液体密度，kg/m3。
优选地，所述裂缝宽度与压力的流固耦合方程的建立过程如下：
岩石变形通过边界元方法计算，计算式如下：
p-σh＝Cw (7)
缝内流动满足质量守恒和层流方程，对缝内流动方程的空间项进行有限体积离散，得到流动方程的一阶微分方程形式：
Figure BDA0002826268040000041
将式(7)带入式(8)，得到所述裂缝宽度与压力的流固耦合方程：
Figure BDA0002826268040000042
井筒模型得到各簇裂缝进液流量后，求解式(9)可得到更新单元宽度和压力分布；同时式(9)得到的井底压力与井筒模型进行对比，直至收敛；
其中，p为缝内压力向量，MPa；σh为最小主应力向量，MPa；C为影响系数矩阵，MPa/m；w为缝宽向量，m；θ为系数，0≤θ≤1；w0、p0分别为上一步宽度和压力分布；A(w)为流动方程的系数矩阵；S为源汇项。
优选地，所述裂缝动边界的建立过程：
建立尖端扩展判断准则，过程如下：
Figure BDA0002826268040000043
其中，d为距尖端的距离，m；
K'＝4(2/π)0.5KIc，MPa·m0.5；
KIc为I型断裂韧性，MPa·m0.5；
E'为平面应变杨氏模量，E'＝E/(1-v2)，MPa；
v-岩石泊松比，无因次；
w为缝宽向量；
式(6)可得到单元发生扩展的临界宽度，通过比较当前时刻尖端单元的宽度是否达到临界宽度，若达到则增加单元，否则单元数量不增加。
优选地，所述建立时间递增模型包括：
时间步计算公式为：
Figure BDA0002826268040000044
时间递增为：
tm＝t+Δt (12)
其中，λ为松弛因子，λ∈(0,1]；
Figure BDA0002826268040000051
为所有待检查单元的最大扩展速度，m/s；
lf为裂缝长度，m；
ε＝0.08；
smax＝301；
ΔtE为显式算法的时间步长，s；
Δx为x方向单元大小，m；
Δy为y方向单元大小，m。
优选地，所述压裂缝内支撑剂运移计算模型的建立过程包括：
裂缝内携砂液流动和支撑剂运移的质量守恒方程为
Figure BDA0002826268040000052
Figure BDA0002826268040000053
Figure BDA0002826268040000054
Figure BDA0002826268040000055
其中，qs、qp分别为携砂液和支撑剂的缝内流量，m3/s；
q1为滤失项；
p为封内流体压力,MPa；
Figure BDA0002826268040000056
为归一化的加砂浓度，无因次；
Figure BDA0002826268040000057
为归一化的支撑剂浓度，
φ为支撑剂体积浓度，无因次；
φm为最大容许浓度，无因次；
Qk为k裂缝的入口流量，m3/s；
a为支撑剂半径，m；
g为重力加速度，m/s2；
ρp、ρf分别为支撑剂和压裂液密度，kg/m3；
ez为z轴方向单位向量；
Figure BDA0002826268040000061
为描述携砂液流动的无因次函数；
Figure BDA0002826268040000062
分别为描述支撑剂对流和沉降的无因次函数；
x为缝内任一点的坐标向量；
xin,k为第k裂缝的入口坐标向量；
δ(x-xin,k)为Delta函数，输入值为0时取1，其他值为0；
B为阻塞函数，描述支撑剂在缝内砂堵现象。
优选地，所述步骤S70包括：
判断时间递增是否达到施工时间，若没有则返回步骤S40，否则结束计算。
本发明了一种水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法，首先建立耦合井筒流动的水平井分段多簇压裂裂缝扩展模型，然后建立支撑剂运移模型，并将其与裂缝扩展模型进行耦合，接着采用加权本质不震荡有限差分方法对耦合方程进行求解，最后得到缝内支撑剂分布。
进一步地，本发明提出了裂缝扩展与支撑剂运移一体化计算方法，其中裂缝扩展是平面三维多裂缝模型，支撑剂运移通过高阶加权本质不震荡方法求解，可解决支撑剂不连续分布问题。裂缝扩展模型为平面三维模型，精度高，而且缝内支撑剂运移可以自动实现砂堵、沉降、流态转变等，并基于高精度方法求解支撑剂分布问题，方法更准确可靠。该方法可应用于水力压裂加砂设计，是提高压裂设计效果和油井增产的有效设计工具。
附图说明
一个或多个实施例通过与之对应的附图中的图片进行示例性说明，这些示例性说明并不构成对实施例的限定，附图中具有相同参考数字标号的元件表示为类似的元件，除非有特别申明，附图中的图不构成比例限制。
图1为压裂井与光纤监测井的物理模型示意图；
图2为裂缝扩展形态的示意图；
图3为支撑剂分布的示意图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明，若本发明实施例中有涉及方向性指示(诸如上、下、左、右、前、后……)，则该方向性指示仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
另外，若本发明实施例中有涉及“第一”、“第二”等的描述，则该“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
本发明提供一种水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,包括：
步骤S10，获取地层参数、施工参数、完井参数；
具体地，所述步骤S10包括：
获取地应力分布、岩石力学参数、注入程序、加砂程序、液体粘度、排量、注入时间、分簇数量、裂缝间距、射孔数、射孔直径。
在本实施例中，具体参数取值参见表1，在其他实施例中，可以根据实际情况而定。
表1基本参数表
Figure BDA0002826268040000071
Figure BDA0002826268040000081
步骤S20：设定初始的缝压和缝宽；
具体地，所述步骤S200的具体过程如下：
假设每个单元的宽度为常数，根据叠加原理，得到初始t时刻裂缝内压力与宽度的离散方程为
Figure BDA0002826268040000082
其中，N为单元总数；
p(x,t)为缝内流体压力，MPa；
σh(x)为最小主应力，MPa；
c(x,ξi)为平面三维裂缝核函数，MPa/m；
w(ξi)为裂缝宽度，m；
ξi为裂缝面网格点的坐标；
x为计算点的坐标；
通过求解式(1)可得到初始时刻的缝压与缝宽；
其中式(1)中的矩阵形式为
p-σh＝Cw (2)
其中，p为缝内压力向量，MPa；
σh为最小主应力向量，MPa；
C为影响系数矩阵，MPa/m；
w为缝宽向量，MPa。
步骤S30：建立耦合井筒流动的水平井分段多簇压裂裂缝扩展计算模型；
具体地，所述步骤S300包括：
步骤S301：建立井筒流动模型；
具体地，所述步骤S301所述井筒流动模型的建立过程如下：
井筒到各簇裂缝的分流满足质量守恒和压力连续条件，式(3)和式(4)构成井筒流动模型，一段开启Nf条裂缝的情况，注入总流量等于各分支流量之和，式(3)和式(4)如下：
Figure BDA0002826268040000091
每条裂缝构成的分支回路具有相同压力降，
pw＝pp,k+pc,k+pin,k (4)
式(2)中
Figure BDA0002826268040000092
式(2)中
Figure BDA0002826268040000093
其中，Qt为一段开启Nf条裂缝的情况，注入总流量，Qi为第i条分支裂缝流量；
pw为井底压力,MPa；
k＝1，2，......，Nf；
pp,k为k裂缝的射孔摩阻，MPa；
pc,k为井口到k裂缝的井筒流动摩阻，MPa；
pin,k为k裂缝的入口压力，MPa；
Qk为k裂缝的入口流量，m3/s；
ρ为液体密度，kg/m3；
nk为k射孔簇的射孔数量；
dk为k射孔簇的射孔直径，mm；
K为射孔磨蚀修正系数，无因次；
fc为沿程摩阻系数，无因次；
Dw为压裂管柱内径，m；lk为井口到k裂缝的管柱长度，m；
ε为压裂管柱的内壁粗糙度，m；
Vw为井筒内液体流速，m/s；
Re为雷诺数，Re＝DwρVw/μ；
μ为液体粘度，mPa·s；
ρ为液体密度，kg/m3。
步骤S302：建立裂缝宽度与压力的流固耦合方程；
具体地，所述裂缝宽度与压力的流固耦合方程的建立过程如下：
将裂缝所在平面离散为一系列矩形单元，岩石变形通过边界元方法计算，计算式如下：
p-σh＝Cw (7)
缝内流动满足质量守恒和层流方程，对缝内流动方程的空间项进行有限体积离散，得到流动方程的一阶微分方程形式：
Figure BDA0002826268040000101
将式(7)带入式(8)，得到所述裂缝宽度与压力的流固耦合方程：
Figure BDA0002826268040000102
井筒模型得到各簇裂缝进液流量后，求解式(9)可得到更新单元宽度和压力分布；同时式(9)得到的井底压力与井筒模型进行对比，直至收敛；
其中，p为缝内压力向量，MPa；σh为最小主应力向量，MPa；C为影响系数矩阵，MPa/m；w为缝宽向量，m；θ为系数，0≤θ≤1；w0、p0分别为上一步宽度和压力分布；A(w)为流动方程的系数矩阵；S为源汇项。
步骤S303：裂缝动边界
具体地，建立尖端扩展判断准则，过程如下：
Figure BDA0002826268040000103
其中，d为距尖端的距离，m；
K'＝4(2/π)0.5KIc，MPa·m0.5；
KIc为I型断裂韧性，MPa·m0.5；
E'为平面应变杨氏模量，E'＝E/(1-v2)，MPa；
v-岩石泊松比，无因次；
w为缝宽向量；
式(6)可得到单元发生扩展的临界宽度，通过比较当前时刻尖端单元的宽度是否达到临界宽度，若达到则增加单元，否则单元数量不增加。
通过步骤S30建立耦合了井筒流动模型的流固耦合方程式(9)，求解流固耦合方程可得到更新单元的宽度和压力分布以及井底压力；得到的井底压力与井筒模型进行对比，直至收敛。
步骤S40：建立时间递增模型；
具体地，所述建立时间递增模型包括：
时间步计算公式为：
Figure BDA0002826268040000111
时间递增为：
tm＝t+Δt (12)
其中，λ为松弛因子，λ∈(0,1]；
lf为裂缝长度，m，可以根据裂缝面积进行近似计算；
为了避免裂缝长度太小时时间步过大，可以选择ε＝0.08；
为了避免积分步数太大而可能出现数值震荡，smax＝301，最大积分步数为301；
Figure BDA0002826268040000112
为所有待检查单元的最大扩展速度，m/s；
ΔtE为显式算法的时间步长，s；
Δx为x方向单元大小，m；
Δy为y方向单元大小，m。
步骤S50：建立压裂缝内支撑剂运移计算模型；
具体地，所述压裂缝内支撑剂运移计算模型的建立过程包括：
考虑支撑剂后，缝内流动为携砂液流动。支撑剂视为携砂液的组分，裂缝内携砂液流动和支撑剂运移的质量守恒方程为
Figure BDA0002826268040000113
Figure BDA0002826268040000114
Figure BDA0002826268040000115
Figure BDA0002826268040000121
其中，w为缝宽向量，m；
qs、qp分别为携砂液和支撑剂的缝内流量，m3/s；
q1为滤失项；
Qk为第k裂缝的入口流量，m3/s；
p为缝内流体压力；
Figure BDA0002826268040000122
为归一化的加砂浓度，无因次；
Figure BDA0002826268040000123
为归一化的支撑剂浓度，
φ为支撑剂体积浓度，无因次；
φm为最大容许浓度，无因次；
a为支撑剂半径，m；
g为重力加速度，m/s2；
μ为压裂液粘度，mPa.s；
ρp、ρf分别为支撑剂和压裂液密度，kg/m3；
ez为z轴方向单位向量；
Figure BDA0002826268040000124
为描述携砂液流动的无因次函数；
Figure BDA0002826268040000125
分别为描述支撑剂对流和沉降的无因次函数；
x为缝内任一点的坐标向量；
xin,k为第k裂缝的入口坐标向量；
δ(x-xin,k)为Delta函数，输入值为0时取1，其他值为0；
B为阻塞函数，描述支撑剂在缝内砂堵现象。
步骤S60：判断该时刻是否加入支撑剂，若是，则采用高阶加权本质不震荡(WENO，Weighted EssentiallyNon-Oscillatory)有限差分方法对压裂裂缝扩展过程中支撑剂运移进行计算，得到缝内支撑剂分布；若不是，则进入下一步；
所述步骤S60包括：
步骤S601：判断该时刻是否加入了支撑剂。因为支撑剂运移方程为非线性双曲方程，其数值求解过程中会出现间断而导致震荡误差，为了解决这一问题并同时要满足特征线稳定性要求，加上本发明中模型为结构化网格，于是采用加权本质不震荡有限差分方法求解支撑剂运移方程。另外，本发明中支撑剂模型与裂缝扩展模型采用相同网格，而裂缝扩展模型的空间步长一般为数米级别，因此支撑剂运移的计算需高阶计算格式。
具体的计算过程：调用式(13)与固体方程进行耦合得到裂缝扩展流固耦合方程，该方程通过显式高效方法求解。求解裂缝扩展模型，得到裂缝宽度、压力等结果后，缝内携砂液的流场即可确定；根据携砂液流场，求解支撑剂运移的连续性方程(14)，即可得到支撑剂分布。
步骤S602：判断无支撑剂加入，则进入下一步。
请参阅图2和图3，通过步骤S60判断该时刻是否加入支撑剂，若是，则调用步骤S50中的携砂液流动方程式(13)与固体方程进行耦合得到裂缝扩展流固耦合方程，该方程通过显式高效方法求解。求解裂缝扩展模型，得到裂缝宽度、压力等结果后，缝内携砂液的流场即可确定；根据携砂液流场，求解支撑剂运移的连续性方程(14)，即可得到支撑剂分布如图2。
步骤S70：判断时间是否达到施工时间，若没有则返回步骤S4，否则结束计算。
具体地，所述步骤S70包括：
判断时间递增是否达到施工时间，若没有则返回步骤S40，否则结束计算。
传统方法缝内支撑剂运移通常为静态裂缝，即将裂缝扩展模型与支撑剂运移模型割裂开来，并没有一体化同步模拟。由于支撑剂运移是非连续问题，本发明采用加权本质不震荡方法，可以实现对支撑剂浓度不连续的准确捕捉，实现支撑剂运移的准确求解。
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是在本发明的发明构思下，利用本发明说明书及附图内容所作的等效结构变换，或直接/间接运用在其他相关的技术领域均包括在本发明的专利保护范围内。

Claims (10)

1.一种水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，包括：
步骤S10，获取地层参数、施工参数、完井参数；
步骤S20：设定初始的缝压和缝宽；
步骤S30：建立耦合井筒流动的水平井分段多簇压裂裂缝扩展计算模型；
步骤S40：建立时间步递增模型；
步骤S50：建立压裂缝内支撑剂运移计算模型；
步骤S60：判断该时刻是否加入支撑剂，若是，则采用高阶加权本质不震荡有限差分方法对压裂裂缝扩展过程中支撑剂运移进行计算，得到缝内支撑剂分布；若不是，则进入下一步；
步骤S70：判断时间是否达到施工时间，若没有则返回步骤S4，否则结束计算。
2.如权利要求1所述的水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，所述步骤S10包括：
获取地应力分布、岩石力学参数、注入程序、加砂程序、液体粘度、排量、注入时间、分簇数量、裂缝间距、射孔数、射孔直径。
3.如权利要求1所述的水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，所述步骤S200的具体过程如下：
假设每个单元的宽度为常数，根据叠加原理，得到初始t时刻裂缝内压力与宽度的离散方程为
Figure FDA0002826268030000011
其中，N为单元总数；
p(x,t)为缝内流体压力，MPa；
σh(x)为最小主应力，MPa；
c(x,ξi)为平面三维裂缝核函数，MPa/m；
w(ξi)为裂缝宽度，m；
ξi为裂缝面网格点的坐标；
x为计算点的坐标；
通过求解式(1)可得到初始时刻的缝压与缝宽；
其中式(1)中的矩阵形式为
p-σh＝Cw (2)
其中，p为缝内压力向量，MPa；
σh为最小主应力向量，MPa；
C为影响系数矩阵，MPa/m；
w为缝宽向量。
4.如权利要求1所述的水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，所述步骤S300包括：
建立井筒流动模型；
建立裂缝宽度与压力的流固耦合方程；
建立裂缝动边界。
5.如权利要求4所述的水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，所述井筒流动模型的建立过程如下：
井筒到各簇裂缝的分流满足质量守恒和压力连续条件，式(3)和式(4)构成井筒流动模型，式(3)和式(4)如下：
Figure FDA0002826268030000021
pw＝pp,k+pc,k+pin,k (4)
式(2)中
Figure FDA0002826268030000022
式(2)中
Figure FDA0002826268030000023
其中，Qt为一段开启Nf条裂缝的情况，注入总流量，Qi为第i条分支裂缝流量；
pw为井底压力,MPa；
k＝1，2，......，Nf；
pp,k为k裂缝的射孔摩阻，MPa；
pc,k为井口到k裂缝的井筒流动摩阻，MPa；
pin,k为k裂缝的入口压力，MPa；
Qk为k裂缝的入口流量，m3/s；
ρ为液体密度，kg/m3；
nk为k射孔簇的射孔数量；
dk为k射孔簇的射孔直径，mm；
K为射孔磨蚀修正系数，无因次；
fc为沿程摩阻系数，无因次；
Dw为压裂管柱内径，m；
lk为井口到k裂缝的管柱长度，m；
ε为压裂管柱的内壁粗糙度，m；
Vw为井筒内液体流速，m/s；
Re为雷诺数，Re＝DwρVw/μ；
μ为液体粘度，mPa·s；
ρ为液体密度，kg/m3。
6.如权利要求4所述的水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，所述裂缝宽度与压力的流固耦合方程的建立过程如下：
岩石变形通过边界元方法计算，计算式如下：
p-σh＝Cw (7)
缝内流动满足质量守恒和层流方程，对缝内流动方程的空间项进行有限体积离散，得到流动方程的一阶微分方程形式：
Figure FDA0002826268030000031
将式(7)带入式(8)，得到所述裂缝宽度与压力的流固耦合方程：
Figure FDA0002826268030000032
井筒模型得到各簇裂缝进液流量后，求解式(9)可得到更新单元宽度和压力分布；同时式(9)得到的井底压力与井筒模型进行对比，直至收敛；
其中，p为缝内压力向量，MPa；
σh为最小主应力向量，MPa；
C为影响系数矩阵，MPa/m；
w为缝宽向量，m；
θ为系数，0≤θ≤1；
w0、p0分别为上一步宽度和压力分布；
A(w)为流动方程的系数矩阵；
S为源汇项。
7.如权利要求4所述的水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，所述裂缝动边界的建立过程：
建立尖端扩展判断准则，过程如下：
Figure FDA0002826268030000041
其中，d为距尖端的距离，m；
K'＝4(2/π)0.5KIc，MPa·m0.5；
KIc为I型断裂韧性，MPa·m0.5；
E'为平面应变杨氏模量，E'＝E/(1-v2)，MPa；
v-岩石泊松比，无因次；
w为缝宽向量，m；
式(6)可得到单元发生扩展的临界宽度，通过比较当前时刻尖端单元的宽度是否达到临界宽度，若达到则增加单元，否则单元数量不增加。
8.如权利要求1所述的水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，所述建立时间递增模型包括：
时间步计算公式为：
Figure FDA0002826268030000042
时间递增为：
tm＝t+Δt (12)
其中，λ为松弛因子，λ∈(0,1]；
lf为裂缝长度，m；
Figure FDA0002826268030000043
为所有待检查单元的最大扩展速度，m/s；
ε＝0.08；
smax＝301；
ΔtE为显式算法的时间步长，s；
Δx为x方向单元大小，m；
Δy为y方向单元大小，m。
9.如权利要求1所述的水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，所述压裂缝内支撑剂运移计算模型的建立过程包括：
裂缝内携砂液流动和支撑剂运移的质量守恒方程为
Figure FDA0002826268030000051
Figure FDA0002826268030000052
Figure FDA0002826268030000053
Figure FDA0002826268030000054
其中，w为缝宽向量，m；
qs、qp分别为携砂液和支撑剂的缝内流量，m3/s；
q1为滤失项；
Qk为第k裂缝的入口流量，m3/s；
Figure FDA0002826268030000055
为归一化的加砂浓度，无因次；
Figure FDA0002826268030000056
为归一化的支撑剂浓度，
φ为支撑剂体积浓度，无因次；
φm为最大容许浓度，无因次；
a为支撑剂半径，m；
g为重力加速度，m/s2；
μ为压裂液粘度，mPa.s；
ρp、ρf分别为支撑剂和压裂液密度，kg/m3；
ez为z轴方向单位向量；
Figure FDA0002826268030000057
为描述携砂液流动的无因次函数；
Figure FDA0002826268030000061
分别为描述支撑剂对流和沉降的无因次函数；
x为缝内任一点的坐标向量；
xin,k为第k裂缝的入口坐标向量；
δ(x-xin,k)为Delta函数，输入值为0时取1，其他值为0；
B为阻塞函数，描述支撑剂在缝内砂堵现象。
10.如权利要求1所述的水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法,其特征在于，所述步骤S70包括：
判断时间递增是否达到施工时间，若没有则返回步骤S40，否则结束计算。
CN202011449599.9A 2020-12-09 2020-12-09 水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法 Active CN112580271B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011449599.9A CN112580271B (zh) 2020-12-09 2020-12-09 水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011449599.9A CN112580271B (zh) 2020-12-09 2020-12-09 水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法

Publications (2)

Publication Number Publication Date
CN112580271A true CN112580271A (zh) 2021-03-30
CN112580271B CN112580271B (zh) 2022-11-04

Family

ID=75131666

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011449599.9A Active CN112580271B (zh) 2020-12-09 2020-12-09 水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法

Country Status (1)

Country Link
CN (1) CN112580271B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113216947A (zh) * 2021-05-17 2021-08-06 中国石油大学(华东) 基于监测井分布式光纤应变监测的水平井压裂过程裂缝高度确定方法
CN113283197A (zh) * 2021-06-10 2021-08-20 西南石油大学 一种基于复杂裂缝缝宽分布的加砂参数设计方法
CN114458280A (zh) * 2022-02-16 2022-05-10 西南石油大学 一种基于cfd-dem考虑颗粒粘结的覆膜支撑剂输送模拟方法
CN115637969A (zh) * 2021-07-20 2023-01-24 中国石油天然气股份有限公司 用于储层的单段多簇均匀起裂的优化方法
CN115758851A (zh) * 2022-11-28 2023-03-07 中国海洋石油集团有限公司 一种含天然裂缝地层裂缝多尺度支撑剂的选择方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089595A (zh) * 2015-05-27 2015-11-25 中国石油天然气股份有限公司 水平压裂裂缝导流作用下的油藏数值模拟方法及装置
CN109064561A (zh) * 2018-08-21 2018-12-21 北京软能创科技有限公司 基于三维拟连续介质水力压裂模型的支撑剂运移模拟方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089595A (zh) * 2015-05-27 2015-11-25 中国石油天然气股份有限公司 水平压裂裂缝导流作用下的油藏数值模拟方法及装置
CN109064561A (zh) * 2018-08-21 2018-12-21 北京软能创科技有限公司 基于三维拟连续介质水力压裂模型的支撑剂运移模拟方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HASEEB ZIA, ETAL: "Explicit versus implicit front advancing schemes for the simulation of hydraulic fracture growth", 《INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS》 *
陈铭等: "水平井分段压裂平面三维多裂缝扩展模型求解算法", 《石油勘探与开发》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113216947A (zh) * 2021-05-17 2021-08-06 中国石油大学(华东) 基于监测井分布式光纤应变监测的水平井压裂过程裂缝高度确定方法
CN113283197A (zh) * 2021-06-10 2021-08-20 西南石油大学 一种基于复杂裂缝缝宽分布的加砂参数设计方法
CN113283197B (zh) * 2021-06-10 2022-04-05 西南石油大学 一种基于复杂裂缝缝宽分布的加砂参数设计方法
CN115637969A (zh) * 2021-07-20 2023-01-24 中国石油天然气股份有限公司 用于储层的单段多簇均匀起裂的优化方法
CN114458280A (zh) * 2022-02-16 2022-05-10 西南石油大学 一种基于cfd-dem考虑颗粒粘结的覆膜支撑剂输送模拟方法
CN115758851A (zh) * 2022-11-28 2023-03-07 中国海洋石油集团有限公司 一种含天然裂缝地层裂缝多尺度支撑剂的选择方法
CN115758851B (zh) * 2022-11-28 2024-01-05 中国海洋石油集团有限公司 一种含天然裂缝地层裂缝多尺度支撑剂的选择方法

Also Published As

Publication number Publication date
CN112580271B (zh) 2022-11-04

Similar Documents

Publication Publication Date Title
CN112580271B (zh) 水平井分段多簇压裂裂缝扩展与支撑剂运移一体化的建模方法
CN112576245B (zh) 一种基于裂缝扩展模拟的分布式光纤应变监测方法
CN110334431B (zh) 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN112576240B (zh) 基于封闭井筒压力波动进行水力压裂裂缝监测的方法
CN113389534B (zh) 一种水平井密切割压裂裂缝扩展预测与设计参数优化方法
WO2022218159A1 (zh) 非常规油气藏水平井压裂缝网扩展与生产动态耦合方法
US10352146B2 (en) Formation fracture flow monitoring
US20150066446A1 (en) Connection conditions for modeling fluid transport in a well system environment
CN114201932B (zh) 一种复杂情况下的致密油藏压裂井试井模拟方法
CN113255123A (zh) 煤层顶板水平井分段多簇压裂适用性地质条件的评价方法
US11733425B1 (en) Full wellbore pressure calculation method, device and computer readable storage medium for fracturing horizontal well
US10294765B2 (en) Formation fracture flow monitoring
CN111125905B (zh) 耦合油藏流体流动的二维裂缝网络扩展模型及其模拟方法
Wen et al. Experimental analysis and numerical simulation of variable mass flow in horizontal wellbore.
CN113360984B (zh) 一种考虑壁面捕获效应的支撑剂输送数值模拟方法
CN107133373B (zh) 一种页岩气藏、井筒及地面管网的耦合模拟方法
Wang et al. Flow simulation of a horizontal well with two types of completions in the frame of a wellbore–annulus–reservoir model
CN115217467A (zh) 多重介质底水气藏水平井底水突破时间的确定方法、装置
CN111502639B (zh) 一种救援井最小压井排量的确定方法
CN114817836A (zh) 一种深井超深井钻井液漏失类型诊断方法
CN113779901A (zh) 一种多井统一计量下的单井动态储量计算方法
US20180210435A1 (en) Method and Device for Optimizing Solid Phase Transport in Pipe Flow
CN117454699A (zh) 一种页岩气水平井多簇射孔压裂下的支撑剂分配计算方法
Mohammed et al. Studying the effect of perforation parameters on vertical well performance
Bi et al. Research on Productivity Prediction Model of Three-Dimensional Directional Wells in Different Reservoirs

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