CN108397184B - 一种自支撑裂缝导流能力的数值计算方法 - Google Patents
一种自支撑裂缝导流能力的数值计算方法 Download PDFInfo
- Publication number
- CN108397184B CN108397184B CN201810483072.4A CN201810483072A CN108397184B CN 108397184 B CN108397184 B CN 108397184B CN 201810483072 A CN201810483072 A CN 201810483072A CN 108397184 B CN108397184 B CN 108397184B
- Authority
- CN
- China
- Prior art keywords
- fracture
- crack
- matrix
- fluid
- fluid particle
- 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
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 16
- 239000012530 fluid Substances 0.000 claims abstract description 92
- 239000011159 matrix material Substances 0.000 claims abstract description 75
- 239000002245 particle Substances 0.000 claims abstract description 60
- 239000000693 micelle Substances 0.000 claims abstract description 47
- 238000006073 displacement reaction Methods 0.000 claims abstract description 30
- 239000011435 rock Substances 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000005315 distribution function Methods 0.000 claims description 27
- 238000009826 distribution Methods 0.000 claims description 19
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 8
- 239000006185 dispersion Substances 0.000 claims description 7
- 230000035772 mutation Effects 0.000 claims description 7
- 239000010454 slate Substances 0.000 claims description 5
- 230000002123 temporal effect Effects 0.000 claims description 4
- 208000010392 Bone Fractures Diseases 0.000 description 76
- 230000006835 compression Effects 0.000 description 6
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 description 4
- 239000003795 chemical substances by application Substances 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000007356 Fracture Dislocation Diseases 0.000 description 1
- 230000005483 Hooke's law Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004215 lattice model Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
- E21B47/00—Survey of boreholes 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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
- E21B43/267—Methods for stimulating production by forming crevices or fractures reinforcing fractures by propping
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)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种自支撑裂缝导流能力的数值计算方法，依次包括以下步骤：（A）对目标储层裂缝面进行单元离散分析，建立裂缝受力变形模型，绘制裂缝应力‑位移图版；（B）利用应力‑位移图版确定目标储层裂缝在已知闭合应力作用下的裂缝变形量；（C）对目标储层岩板表面进行激光扫描，计算裂缝宽度矩阵；（D）根据格子玻尔兹曼方法，计算裂缝内流体粒子微团的流动速度、密度及压力；（E）计算裂缝内流体流量Q、裂缝入口压力P1和裂缝出口压力P2；（F）使用达西公式计算裂缝导流能力。本发明原理可靠，操作简便，能够简单、准确地计算自支撑裂缝的导流能力，对水力压裂具有重要的指导意义。
Description
技术领域
本发明涉及石油领域，尤其是页岩清水压裂过程中一种自支撑裂缝导流能力的数值计算方法。
背景技术
水力压裂技术是低渗透油气藏增产改造的重要措施。水力压裂是利用地面高压泵组，以超过地层吸收能力的排量将压裂液泵入地层来产生裂缝，然后继续注入带有支撑剂(砂粒)的压裂液，使裂缝继续延伸并在其中充填支撑剂，当压裂液返排后，在地层压力作用下，支撑剂在裂缝中起到支撑裂缝的作用，阻止裂缝闭合，从而在地层中形成具有一定长度和导流能力的填砂裂缝。
清水压裂是水力压裂的一种形式，被广泛应用于页岩油气藏的增产改造中。它具体是指在压裂过程中不加入支撑剂(砂粒)，仅通过将低粘度的压裂液泵入地层来产生裂缝。地下岩石性质差异较大，压裂形成的裂缝表面一般凹凸不平，同时还会在剪切作用下发生错位，因此即使不加入支撑剂，裂缝表面的凸点之间也可以相互支撑，形成自支撑裂缝，自支撑裂缝在地层压力作用下仍能保持一定开启程度，并具有一定的导流能力，从而达到改善油气流动条件和油气井增产的目的。因此，自支撑裂缝的导流能力是评价水力压裂(清水压裂)成功与否的重要指标之一。
发明内容
本发明的目的在于提供一种自支撑裂缝导流能力的数值计算方法，该方法将目标储层页岩加工为长方体岩样，并劈裂分为两个具有粗糙表面的岩板，将两个岩板的粗糙表面错位放置形成错位裂缝，通过对错位裂缝进行变形分析和数值建模，最终计算出裂缝导流能力，该方法原理可靠，操作简便，能够简单、准确地计算自支撑裂缝的导流能力，对水力压裂具有重要的指导意义。
为达到以上技术目的，本发明采用以下技术方案：对自支撑错位裂缝的受力形态进行建模，绘制裂缝的“应力-位移”图版，基于图版计算裂缝受力变形后的裂缝宽度数据，再基于格子玻尔兹曼方程计算裂缝中的流体流速分布及密度分布数据，进而得到裂缝中的流体流量和压力数据，最终利用达西公式计算出裂缝导流能力。
利用三维激光扫描仪按专利CN201510319382.9所述的方法对岩板裂缝面进行扫描，获取裂缝表面的三维数据，并建立岩板的三维坐标系：取岩板长度方向为横轴方向，宽度方向为纵轴方向，高度方向为竖轴方向，且横轴和纵轴方向所在的面为水平面。
一种自支撑裂缝导流能力的数值计算方法，依次包括以下步骤：
(A)对目标储层裂缝面进行单元离散分析，建立裂缝受力变形模型，绘制裂缝应力-位移图版；
(B)利用步骤(A)绘制的应力-位移图版确定目标储层裂缝在已知闭合应力σf作用下的裂缝变形量Zf；
(C)对目标储层岩板表面进行激光扫描，基于步骤(B)计算裂缝宽度矩阵Wf；
(D)根据格子玻尔兹曼方法，利用下式计算裂缝内流体粒子微团的流速、密度及压力数据：
式中：α—流体粒子微团的规定运动方向；α＝0，1，2……18；
δt—离散时间步长；
τ0—松弛时间，一般计算为τ0＝3μ+0.5，μ为流体粘度；
ρ—时间为t时，流体粒子微团密度；
p—时间为t时，流体粒子微团的压力；
(E)基于步骤(D)，利用下式计算裂缝内流体流量Q，裂缝入口压力P1和裂缝出口压力P2：
式中：Wfin—裂缝入口处的裂缝宽度矩阵；
Ufx—流体粒子微团在裂缝长度方向上的流速分布矩阵；
X—对目标储层裂缝面进行单元离散分析时，离散单元的边长；
wfinj—矩阵Wfin中的元素；
ufxj—矩阵Ufx中的元素；
pinj—裂缝入口处的压力分布矩阵Pin中的元素；
poutj—裂缝出口处的压力分布矩阵Pout中的元素；
k—对裂缝面进行扫描后，纵轴方向共有k列个扫描数点；
(F)使用以下达西公式计算裂缝导流能力：
式中：Fcd—裂缝导流能力；
μ—流体粘度；
L—裂缝(岩板)长度；
B—裂缝(岩板)宽度。
本发明中，所述步骤(A)对目标储层裂缝面进行单元离散分析，建立裂缝受力变形模型，绘制裂缝应力-位移图版，过程如下：
对扫描裂缝面三维数据进行坐标变换，使两个裂缝表面刚好接触，对上裂缝面施加向下的初始位移量Z0，在竖轴方向上选定一个离散单元，离散单元截面为正方形，边长为X，上裂缝面高度Z1，竖轴方向上所受的压力为ΔFz，压缩量为ΔZ1，在竖轴方向所受的应力值为：
式中：σm为裂缝岩体的抗压强度；
Mc为岩体的应力突变系数(本发明中，推荐页岩应力突变系数Mc取值0.5)；
υ为裂缝岩体的泊松比。
从而求得上裂缝面所受的应力值σ：
式中：m—对裂缝面进行立方单元离散后，横轴方向共有m排离散单元；
n—对裂缝面进行立方单元离散后，纵轴方向共有n列离散单元；
Δσi,j—横轴方向第i排，纵轴方向第j列的离散单元所受的应力值。
以Z表示上裂缝面的向下位移量，以σ表示上裂缝面所受的应力值，通过对上裂缝面施加不同的向下位移量，不断重复上述步骤，直到获得一定数量的(σ，Z)数据点，绘制σ-Z图版，得到Z＝f(σ)曲线。
本发明中，所述步骤(B)利用步骤(A)绘制的应力-位移图版确定目标储层裂缝在已知闭合应力σf作用下的裂缝变形量Zf，过程如下：在步骤(A)中，上裂缝面的向下位移量即为实际裂缝的变形量，所以，通过σ-Z图版，给定任意应力值σ，利用Z＝f(σ)曲线即可得到相应的裂缝变形量Z，从而可得到目标储层裂缝在已知闭合应力σf作用下的裂缝变形量Zf。
本发明中，步骤(C)对目标储层岩板表面进行激光扫描，基于步骤(B)计算裂缝宽度矩阵Wf，过程如下：对扫描裂缝面三维数据进行坐标变换，使两个裂缝表面刚好接触，下裂缝面表面每个点的竖轴数据组成下裂缝面高度矩阵Zd，上裂缝面表面每个点的竖轴数据组成上裂缝面高度矩阵Zu。通过步骤(B)已确定了在闭合应力σf作用下，裂缝变形量为Zf，即上裂缝面的向下位移量为Zf。则在σf作用下，上裂缝面高度矩阵变为Z′u：
Z′u＝(zu′ij)＝(zuij-Zf)，(i＝1,2,3,...,h；j＝1,2,3,...,k) (3)
式中：z′uij—Z′u矩阵中，横轴方向第i排，纵轴方向第j列的点的高度值；
zuij—Zu矩阵中，横轴方向第i排，纵轴方向第j列的点的高度值。
获得裂缝形态变形矩阵Zc：
Zc＝Z′u-Zd (4)
从而得到裂缝宽度矩阵Wf为：
式中：wfij为裂缝中横轴方向第i排，纵轴方向第j列的点的裂缝宽度值；
zcij为裂缝形态变形矩阵Zc矩阵中横轴方向第i排，纵轴方向第j列的点的高度值。
本发明中，步骤(D)根据格子玻尔兹曼方法，计算裂缝内流体流速、密度及压力数据，过程如下：
玻尔兹曼(Boltzmann)方程(杨大勇.格子玻尔兹曼方法[M].电子工业出版社.2015)是描述粒子微团速度分布函数(简称“分布函数”)
在时间和空间上变化的守恒方程，裂缝中的流体流动满足离散格子Boltzmann方程耦合BGK方程(简称LBGK方程)：
本发明选用常见的三维格子玻尔兹曼速度模型D3Q19单速模型进行计算，该模型具有离散速度的对称性，粒子演化机制简单，是目前广泛使用的三维格子玻尔兹曼模型。
在D3Q19单速模型条件下，式(6)中：
α—流体粒子微团的规定运动方向；α＝0，1，2……18；
δt—离散时间步长；
τ0—松弛时间，一般计算为τ0＝3μ+0.5，μ为流体粘度；
式(6)描述了流体粒子微团的微观运动，通过对式(6)进行求解可获得流体粒子微团在某一位置上的分布函数
求得分布函数之后，可以利用分布函数，根据以下三个公式求得流体粒子微团的密度、流速和压力：
式中：ρ—时间为t时，流体粒子微团密度；
p—时间为t时，流体粒子微团的压力。
要对式(6)进行求解，首先需要求得式(6)中的平衡分布函数
其表统一表达式为：
式中：ωα—权重系数；
CS—格子声速。
在D3Q19单速模型格子模型中：
式(12)中，c为声速，即340m/s。
式(8)与式(11)中，取
在横轴方向的分向量代入式(8)计算得到流体粒子微团在裂缝长度方向上的流速
同理，取
在纵轴方向的分向量代入式(8)计算得到流体粒子微团在裂缝宽度方向上的流速
对裂缝表面进行网格划分，在裂缝长度方向上划分h(i＝1,2,3,...,h)个节点，裂缝宽度方向上划分k个节点，保证网格划分数量与三维激光扫描仪扫描后得到的网格数量一致。在每个网格的节点处即是一个流体微团，利用计算机编程在离散时间步长为δt条件下迭代计算。根据格子玻尔兹曼发展经验，设定初始状态(t＝0时刻)下，裂缝内流体密度ρt＝0为1g/cm3，此时的分布函数
计算式为：
根据式(8)计算得到每个节点处流体微团流动速度
并根据式(10)计算t＝0时刻的平衡分布函数：
t＝0起，每个流体微团开始运动，采用定流量边界条件(Lian-PingWang,and ChengPeng.LatticeBoltzmannsimulationofparticle-ladenturbulentchannelflow[J].Computers and Fluids.2015:1-11.)计算t＝δt时刻的流体微团分布函数：
利用
结合式(7)、式(8)、式(9)计算得到t＝δt时刻的密度分布数据
速度分布数据
以及压力分布数据
根据式(10)计算t＝δt时刻的平衡分布函数，以此开始第二次迭代。以此反复迭代计算，直到t＝nδt时刻(n为自然数)达到式(17)的收敛条件：
式中：
—时间为t＝(n-1)δt时，流体粒子微团在裂缝长度方向上的流速；
本发明中，所述步骤(E)计算裂缝内流体流量Q，裂缝入口压力P1和裂缝出口压力P2，过程如下：
获得裂缝宽度矩阵Wf后，在矩阵Wf中取i＝1条件下的所有元素组成裂缝入口处的裂缝宽度矩阵Wfin。基于达到收敛条件的流体微团分布函数
可计算出网格数i＝1条件下，所有节点处流体粒子微团在裂缝长度方向上的流速分布数据并组成矩阵Ufx，矩阵Ufx即为裂缝入口处流体粒子微团在裂缝长度方向上的流速分布矩阵。则裂缝内流量Q为：
式中：wfinj—矩阵Wfin中的元素；
ufxj—矩阵Ufx中的元素；
X—离散单元截面边长。
基于达到收敛条件的流体微团分布函数
可计算出网格数i＝1条件下，所有节点处流体粒子微团的压力数据并组成矩阵Pin，矩阵Pin即为裂缝入口处流体粒子微团的压力分布矩阵。同理，基于达到收敛条件的流体微团分布函数
可计算出网格数i＝h条件下，所有节点处流体粒子微团的压力数据并组成矩阵Pout，矩阵Pout即为裂缝出口处流体粒子微团的压力分布矩阵。采用求平均值方法得到裂缝入口处力P1和出口处压力P2：
式中：pinj—矩阵Pin中的元素；
poutj—矩阵Pout中的元素。
附图说明
图1为裂缝空间直角坐标系的任意一个横-竖轴剖面示意图。
图2为选定离散单元受力分析示意图。
具体实施方式
下面根据附图进一步说明本发明。
本发明中，对错位裂缝受力变形进行分析并建立模型，具体过程为：
如图1所示，对扫描裂缝面三维数据进行坐标变换，将下裂缝面原始数中所有点的高度值减去下裂缝面的最低点高度值，即：
Zd＝(zdij)＝[z0ij-min(z0ij)],(i＝1,2,3,...,h；j＝1,2,3,...,k) (20)
上式中：Zd—下裂缝面高度矩阵；
zdij—下裂缝面高度矩阵中，横轴方向第i排，纵轴方向第j列的点的高度值；
z0ij—原始扫描数据度矩阵中，横轴方向第i排，纵轴方向第j列的点的高度值；
h—对裂缝面进行扫描后，横轴方向共有h列个扫描数点；
k—对裂缝面进行扫描后，纵轴方向共有k列个扫描数点；
min(z0ij)—原始扫描数据高度矩阵中最低点的高度值。
同样，将上裂缝面原始数中所有点的高度值减去上裂缝面的最低点高度值，再将处理后上裂缝面数据点经过坐标变换按图1所示放入下裂缝面的坐标系中，使两个裂缝表面刚好接触。在图1所示的坐标系中，下裂缝面表面每个点的竖轴数据即组成下裂缝面高度矩阵Zd，上裂缝面表面每个点的竖轴数据即组成上裂缝面高度矩阵Zu。
保持下裂缝面不变，对上裂缝面施加向下的初始位移量Z0后，上裂缝面高度矩阵变为Z′u：
Z′u＝(z′uij)＝(zuij-Z0)，(i＝1,2,3,...,h；j＝1,2,3,...,k) (21)
上式中：z′uij—Z′u矩阵中，横轴方向第i排，纵轴方向第j列的点的高度值；
zuij—Zu矩阵中，横轴方向第i排，纵轴方向第j列的点的高度值；
利用式(4)计算获得裂缝形态变形矩阵Zc
在矩阵Zc中，当元素zcij(Zc矩阵中，横轴方向第i排，纵轴方向第j列的点的高度值)大于零时，表示该点在位移量为Z0时，两个裂缝面没有发生接触；当zcij小于或等于零时，表示该点在位移量为Z0时两个裂缝面发生了接触变形，且压缩量为|zcij|。
定义上裂缝面与下裂缝面的总压缩量矩阵为Zt，矩阵中的元素为ztij：
在竖轴方向上对裂缝面进行立方单元离散，选定一个离散单元，如图2所示，已知：上裂缝面高度Z1，下裂缝面高度Z2，并通过矩阵Zt可得到上裂缝面与下裂缝面的总压缩量ΔZ，离散单元截面为正方形，边长为X，裂缝岩体杨氏模量为E，泊松比为υ，抗压强度为σm。
对上裂缝面离散单元进行受力分析，假设为纯弹性应变，其竖轴方向上所受的压力为ΔFz，应力为Δσ，压缩量为ΔZ1，应变为εz1，根据胡克定律有：
假设上裂缝面离散单元横向变形为ΔX1，则变形后的离散单元在竖轴方向受力面积为ΔA，则有：
将公式(24)代入公式(23)中，得到上裂缝面离散单元的受力变形方程：
同理，假设下裂缝面离散单元横向变形为ΔX2，可推出下裂缝面离散单元的受力变形方程为：
设离散单元上裂缝面在横轴方向的应变为εx1，根据泊松比定义得：
由于上、下粗糙面为同一性质材料，因此有：
离散单元上裂缝面与下裂缝面的总压缩量为ΔZ：
ΔZ1+ΔZ2＝ΔZ (29)
由公式(25)、(26)、(28)、(29)可建立如下计算方程组：
上述方程组含五个方程，可求解得到五个未知数(ΔZ1、ΔZ2、ΔX1、ΔX2、ΔFz)。
由公式(28)可得到ΔX1与ΔZ1关系以及ΔX2与ΔZ2关系
将公式(31)代入公式(25)中
上式可整理为：
同理可得：
由式(34)、式(35)可知，求解出ΔZ1或ΔZ2即可计算获得微元受力ΔFz值。
因此，联立公式(34)、(35)、(29)可得ΔZ1计算表达式：
上式可整理出公式(37)如下：
上式可使用牛顿迭代法进行数值求解得到选定离散单元中上裂缝面的压缩量ΔZ1。迭代过程中由物理对象客观条件即计算值为正且变形量不超过计算边界，设定迭代起始值及迭代精度。本发明推荐迭代起始值为1，迭代精度10-8。
至此，可利用公式(34)计算出选定离散单元中上裂缝面所受的压力ΔFz。
公式(37)和公式(34)的推导过程假设了离散单位为纯弹性应变，实际上，随着岩体所受应力的增加，当岩体达到抗压强度后，岩样受到瞬间破坏，应力不再遵循纯弹性应变规律，故需要进一步判断选定离散单元所受的应力是否达到抗压强度：
通过公式(37)和公式(34)可计算获得ΔFz、ΔX1，则选定离散单元在竖轴方向上所受的应力值Δσ为：
将公式(31)代入公式(38)中得：
裂缝岩体的抗压强度为当σm时，当Δσ＜σm时，岩体发生线弹性变形，应力值为
当σ0≥σm时，页岩发生应力损伤破坏，应力值为σmMc。此处，Mc为岩体的应力突变系数，该系数定义为：岩体受力超过岩体抗压强度σm后，岩体发生瞬间破裂失效，并保持某一残余应力值，失效后的残余应力值与抗压强度的比值为应力突变系数Mc。故选定离散单元在竖轴方向所受的应力值为：
本发明在实验室采用大量页岩样本进行应力—位移测试，得到的应力突变系数分布在0.4-0.6范围内。因此本发明推荐页岩应力突变系数Mc取值0.5。
可利用求和公式计算上裂缝面所受的应力值：
上式中：σ—上裂缝面所受的应力值；
m—对裂缝面进行立方单元离散后，横轴方向共有m排离散单元；
n—对裂缝面进行立方单元离散后，纵轴方向共有n列离散单元；
Δσi,j—横轴方向第i排，纵轴方向第j列的离散单元所受的应力值；
此时，已经计算出了在一个向下的初始位移量Z0条件下，上裂缝面所受的应力值σ，现给定一个位移量步长T，将初始位移量增加T，继续计算上裂缝面所受的应力值。
以Z表示上裂缝面的位移量，以σ表示上裂缝面所受的应力值，不断重复上述步骤，直到获得一定数量的(σ，Z)数据点，利用(σ，Z)数据点绘制σ-Z图版，即可得到Z＝f(σ)曲线，在σ-Z图版中，给定任意应力值σ，利用Z＝f(σ)曲线得到相应的裂缝变形量Z。
自支撑裂缝的宽度一般是毫米级，为了尽可能获取更多的数据点，本发明建议初始位移量取值0.1mm，位移量步长T取值0.1mm。
一种自支撑裂缝导流能力的数值计算方法，依次包括以下步骤：
(A)对目标储层裂缝面进行单元离散分析，建立裂缝受力变形模型，绘制裂缝应力-位移图版；
(B)利用步骤(A)绘制的应力-位移图版确定目标储层裂缝在已知闭合应力σf作用下的裂缝变形量Zf；
(C)对目标储层岩板表面进行激光扫描，基于步骤(B)计算裂缝宽度矩阵Wf；
(D)根据格子玻尔兹曼方法，计算裂缝内流体流速、密度及压力数据；
(E)基于步骤(D)，计算裂缝内流体流量Q，裂缝入口压力P1和裂缝出口压力P2；
(F)使用以下公式计算裂缝导流能力：
式中：Fcd—裂缝导流能力；
μ—流体粘度；
L—裂缝(岩板)长度；
B—裂缝(岩板)宽度。
Claims (4)
1.一种自支撑裂缝导流能力的数值计算方法，依次包括以下步骤：
(A)对目标储层裂缝面进行单元离散分析，建立裂缝受力变形模型，绘制裂缝应力-位移图版，过程如下：对扫描裂缝面三维数据进行坐标变换，使两个裂缝表面刚好接触，对上裂缝面施加向下的初始位移量Z0，在竖轴方向上选定一个离散单元，离散单元截面为正方形，边长为X，上裂缝面高度Z1，竖轴方向上所受的压力为ΔFz，压缩量为ΔZ1，在竖轴方向所受的应力值为：
式中：σm为裂缝岩体的抗压强度；
Mc为岩体的应力突变系数；
υ为裂缝岩体的泊松比；
利用下式求得上裂缝面所受的应力值σ：
式中：m—对裂缝面进行立方单元离散后，横轴方向共有m排离散单元；
n—对裂缝面进行立方单元离散后，纵轴方向共有n列离散单元；
Δσi,j—横轴方向第i排，纵轴方向第j列的离散单元所受的应力值；
以Z表示上裂缝面的向下位移量，以σ表示上裂缝面所受的应力值，通过对上裂缝面施加不同的向下位移量，获得一定数量的(σ，Z)数据点，绘制σ-Z图版，得到Z＝f(σ)曲线；
(B)利用步骤(A)绘制的应力-位移图版确定目标储层裂缝在已知闭合应力σf作用下的裂缝变形量Zf；
(C)对目标储层岩板表面进行激光扫描，基于步骤(B)计算裂缝宽度矩阵Wf；
(D)根据格子玻尔兹曼方法，利用下式计算裂缝内流体粒子微团的流动速度、密度及压力：
式中：α—流体粒子微团的规定运动方向；α＝0，1，2……18；
δt—离散时间步长；
τ0—松弛时间，一般计算为τ0＝3μ+0.5，μ为流体粘度；
ρ—时间为t时，流体粒子微团密度；
p—时间为t时，流体粒子微团的压力；
(E)基于步骤(D)，利用下式计算裂缝内流体流量Q，裂缝入口压力P1和裂缝出口压力P2：
式中：Wfin—裂缝入口处的裂缝宽度矩阵；
Ufx—流体粒子微团在裂缝长度方向上的流速分布矩阵；
X—对目标储层裂缝面进行单元离散分析时，离散单元的边长；
wfinj—矩阵Wfin中的元素；
ufxj—矩阵Ufx中的元素；
pinj—裂缝入口处的压力分布矩阵Pin中的元素；
poutj—裂缝出口处的压力分布矩阵Pout中的元素；
k—对裂缝面进行扫描后，纵轴方向共有k列个扫描数点；
(F)使用以下公式计算裂缝导流能力：
式中：Fcd—裂缝导流能力；
μ—流体粘度；
L—裂缝长度；
B—裂缝宽度。
2.如权利要求1所述的一种自支撑裂缝导流能力的数值计算方法，其特征在于，所述步骤(C)对目标储层岩板表面进行激光扫描，基于步骤(B)计算裂缝宽度矩阵Wf，过程如下：对扫描裂缝面三维数据进行坐标变换，使两个裂缝表面刚好接触，下裂缝面表面每个点的竖轴数据组成下裂缝面高度矩阵Zd，上裂缝面表面每个点的竖轴数据组成上裂缝面高度矩阵Zu，上裂缝面的向下位移量为Zf，上裂缝面高度矩阵变为Z′u：
Z′u＝(z′uij)＝(zuij-Zf)，(i＝1,2,3,...,h；j＝1,2,3,...,k)
式中：z′uij—Z′u矩阵中，横轴方向第i排，纵轴方向第j列的点的高度值；
zuij—Zu矩阵中，横轴方向第i排，纵轴方向第j列的点的高度值；
裂缝形态变形矩阵Zc：
Zc＝Z′u-Zd
裂缝宽度矩阵Wf为：
式中：wfij为裂缝中横轴方向第i排，纵轴方向第j列的点的裂缝宽度值；
zcij为裂缝形态变形矩阵Zc矩阵中横轴方向第i排，纵轴方向第j列的点的高度值。
3.如权利要求1所述的一种自支撑裂缝导流能力的数值计算方法，其特征在于，所述步骤(D)计算裂缝内流体粒子微团的流动速度、密度及压力，过程如下：
裂缝中的流体流动满足LBGK方程：
α—流体粒子微团的规定运动方向；α＝0，1，2……18；
δt—离散时间步长；
τ0—松弛时间，一般计算为τ0＝3μ+0.5，μ为流体粘度；
通过求解获得流体粒子微团在某一位置上的分布函数
根据以下公式求得流体粒子微团的密度ρ、流动速度
和压力p：
4.如权利要求1所述的一种自支撑裂缝导流能力的数值计算方法，其特征在于，所述步骤(E)计算裂缝内流体流量Q，裂缝入口压力P1和裂缝出口压力P2，过程如下：
在裂缝宽度矩阵Wf中取i＝1条件下的所有元素组成裂缝入口处的裂缝宽度矩阵Wfin；基于达到收敛条件的流体微团分布函数
计算出网格数i＝1条件下，所有节点处流体粒子微团在裂缝长度方向上的流速分布数据，组成裂缝入口处流体粒子微团在裂缝长度方向上的流速分布矩阵Ufx，通过下式计算裂缝内流量Q：
式中：wfinj—矩阵Wfin中的元素；
ufxj—矩阵Ufx中的元素；
X—离散单元截面边长；
基于达到收敛条件的流体微团分布函数
计算出网格数i＝1条件下，所有节点处流体粒子微团的压力数据并组成裂缝入口处流体粒子微团的压力分布矩阵Pin；计算出网格数i＝h条件下，所有节点处流体粒子微团的压力数据并组成裂缝出口处流体粒子微团的压力分布矩阵Pout，通过下式得到裂缝入口处力P1和出口处压力P2：
式中：pinj—矩阵Pin中的元素；
poutj—矩阵Pout中的元素。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810483072.4A CN108397184B (zh) | 2018-05-18 | 2018-05-18 | 一种自支撑裂缝导流能力的数值计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810483072.4A CN108397184B (zh) | 2018-05-18 | 2018-05-18 | 一种自支撑裂缝导流能力的数值计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108397184A CN108397184A (zh) | 2018-08-14 |
| CN108397184B true CN108397184B (zh) | 2021-04-20 |
Family
ID=63102325
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810483072.4A Active CN108397184B (zh) | 2018-05-18 | 2018-05-18 | 一种自支撑裂缝导流能力的数值计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108397184B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111305806B (zh) * | 2018-11-27 | 2022-06-03 | 中国石油天然气股份有限公司 | 自支撑裂缝导流能力的分析方法及装置 |
| CN110396399A (zh) * | 2019-06-18 | 2019-11-01 | 中国石油天然气股份有限公司 | 一种油水井大漏失套损段封堵材料及封堵方法 |
| CN110502825B (zh) * | 2019-08-19 | 2023-04-07 | 青岛理工大学 | 一种提取三维破裂面的方法 |
| CN115032368B (zh) * | 2022-06-07 | 2023-08-04 | 西南石油大学 | 一种压裂裂缝自支撑导流能力全过程评价方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5074359A (en) * | 1989-11-06 | 1991-12-24 | Atlantic Richfield Company | Method for hydraulic fracturing cased wellbores |
| AU2006336479B2 (en) * | 2006-01-27 | 2011-03-31 | Schlumberger Technology B.V. | Method for hydraulic fracturing of subterranean formation |
| US9677393B2 (en) * | 2013-08-28 | 2017-06-13 | Schlumberger Technology Corporation | Method for performing a stimulation operation with proppant placement at a wellsite |
| CN105298488A (zh) * | 2015-12-03 | 2016-02-03 | 中国石油集团川庆钻探工程有限公司 | 非连续充填方式下导流能力测试方法 |
| CN107085638B (zh) * | 2017-04-17 | 2019-11-12 | 西南石油大学 | 一种水力压裂支撑剂参数优化方法 |
| CN107545113B (zh) * | 2017-09-08 | 2020-02-07 | 西南石油大学 | 非常规油气藏水力压裂复杂缝网形成过程模拟方法 |
-
2018
- 2018-05-18 CN CN201810483072.4A patent/CN108397184B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108397184A (zh) | 2018-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108397184B (zh) | 一种自支撑裂缝导流能力的数值计算方法 | |
| CN106874544B (zh) | 一种页岩储层改造体积的地质表征方法 | |
| Psakhie et al. | Approach to simulation of deformation and fracture of hierarchically organized heterogeneous media, including contrast media | |
| Ma et al. | Enhancement factors for grounded ice and ice shelves inferred from an anisotropic ice-flow model | |
| US20180355701A1 (en) | Hydraulic fracturing simulation | |
| CN109505576A (zh) | 页岩水力压裂三维全耦合离散裂缝网络模拟方法及系统 | |
| Miguel et al. | Influence of size on the constitutive equations of concrete or rock dowels | |
| Dündar et al. | Three-dimensional fracture and fatigue crack propagation analysis in structures with multiple cracks | |
| CN108489809B (zh) | 利用实验手段计算应力作用下粗糙错位裂缝变形量的方法 | |
| Marder et al. | Simple models of the hydrofracture process | |
| CN113887045B (zh) | 一种暂堵裂缝动态压力与扩展轨迹的预测方法 | |
| CN107578471B (zh) | 一种自支撑裂缝初始形态构建方法 | |
| CN113033049B (zh) | 一种地层尺度下的粗糙裂缝内支撑剂输送数值模拟方法 | |
| Holtzman et al. | Mechanical properties of granular materials: A variational approach to grain‐scale simulations | |
| Nie et al. | Continuous-discontinuous element method for three-dimensional thermal cracking of rocks | |
| Zhu et al. | A fracture conductivity model for channel fracturing based on lattice-Boltzmann-method and computational-fluid-dynamics | |
| CN116702539A (zh) | 滑坡冲击作用下建筑物易损性的评估方法 | |
| Azizan et al. | Numerical prediction of stress and displacement of ageing concrete dam due to alkali-aggregate and thermal chemical reaction | |
| Weng et al. | Exploring the evolution of lateral earth pressure using the distinct element method | |
| Huang et al. | A coupled nonlocal damage model for hydraulic fracture propagation | |
| Bikulov et al. | Prediction of the permeability of proppant packs under load | |
| Tabiei et al. | Evaluation of various numerical methods in LS-DYNA® for 3D Crack Propagation | |
| Maliaris et al. | Modeling of open cell structures geometry and mechanical response applying the Voronoi tessellation algorithm | |
| Palassi et al. | Comparison of 2D and 3D distinct element analyses in stability of convex jointed rock slopes | |
| Jarosch et al. | Brief communication: Stalagmite damage by cave-ice flow quantitatively assessed by fluid-structure-interaction simulations |
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 |