本发明专利技术公开了一种砾岩水力裂缝扩展方向的预测方法，包括以下步骤：获取地质参数和压裂施工参数；通过砾石形状及方位、砾石空间位置分布、砾石粒径、砾石含量生成含砾岩的几何模型；对基质和砾石的参数分别赋值；建立压裂裂缝扩展的岩石变形模型；建立压裂裂缝扩展的流体流动模型；建立水力裂缝与砾石相交作用准则；通过改变砾石大小和砾石含量探究砾石对水力裂缝扩展的影响规律。本发明专利技术基于二维位移不连续法，建立了水力裂缝与砾石相交作用准则，能够准确的模拟水力裂缝与砾石相交后的复杂轨迹，有助于揭示水力裂缝与砾石相互作用机理。理。理。
【技术实现步骤摘要】
一种砾岩水力裂缝扩展方向的预测方法
[0001]本专利技术涉及一种砾岩水力裂缝扩展方向的预测方法，属于油气田增产改造领域。
技术介绍
[0002]水平井体积压裂已成为新疆玛湖砾岩油藏等非常规致密油藏的主要开发方式，准确模拟砾石影响下裂缝扩展规律，对满足国内油气资源和保障国家能源安全至关重要。
[0003]砾岩储层具有非均质性强、基础物性差、低孔、低渗、孔隙类型多样和孔隙结构复杂等特点，常规的水力压裂不能获得有效的工业产能，因此需要水平井体积压裂技术，该技术在页岩储层开发已经非常成熟，但砾岩储层由于砾石的存在对水力裂缝延伸轨迹具有一定影响。因此需要探究砾石对水力裂缝的影响规律，地层中砾石的表征也尤为关键，一般包括砾石形状、砾石方位、砾石空间位置分布、砾石粒径以及砾石含量等。水力裂缝与砾石相交后会产生多种复杂的相互作用结果，即出现止砾、穿砾、绕砾等多种复杂压裂裂缝特征，目前尚未建立一套较为完善的裂缝扩展判断准则。
技术实现思路
[0004]本专利技术为了克服现有技术中未形成一套较为完善的裂缝扩展判断准则，提出了一种砾岩水力裂缝扩展方向的预测方法。
[0005]本专利技术解决上述技术问题所提供的技术方案是：一种砾岩水力裂缝扩展方向的预测方法，包括以下步骤：
[0006]步骤1：获取地质参数和压裂施工参数；
[0007]步骤2：通过砾石形状及方位、砾石空间位置分布、砾石粒径、砾石含量生成含砾岩的几何模型；
[0008]步骤3：对基质和砾石的参数分别赋值；
[0009]步骤4：建立压裂裂缝扩展的岩石变形模型；
[0010]步骤5：建立压裂裂缝扩展的流体流动模型；
[0011]步骤6：建立水力裂缝与砾石相交作用准则；
[0012]步骤7：通过改变砾石大小和砾石含量探究砾石对水力裂缝扩展的影响规律。
[0013]进一步的技术方案是，所述步骤1的地质参数包括储层厚度、砾石形状及方位、砾石空间位置分布、砾石粒径、砾石含量、水平最大主应力、水平最小主应力、杨氏模量、泊松比、基质和砾石断裂韧性、胶结面断裂韧性。压裂施工参数包括注入速率、压裂液黏度和压裂液密度。
[0014]进一步的技术方案是，所述步骤2的砾石形状及方位采用不规则多边形来表征，在极坐标系统中，确定多边形的形状由顶点数m，极角θ
i
以及极半径R
i
分别表示为：
[0015]θ
i
＝η1×
2π
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0016]R
i
＝A0+(2η
i
‑
1)
×
A1ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0017]式中：η1为0～1之间的随机数；A0为极半径R
i
的平均值，mm；A1为极半径R
i
的增加值，
mm；η
i
为0～1之间的随机数。
[0018]砾石空间位置分布可表示为：
[0019][0020]式中：W
min
为研究区域横坐标最小值，m；W
max
研究区域横坐标最大值，m；H
min
为研究区域纵坐标最小值，m；H
max
研究区域纵坐标最大值，m；λ1，λ2为0～1之间的随机数。
[0021]砾石粒径分布的概率密度函数表示为：
[0022][0023]式中：L
a
为长轴大小，mm；u为砾石长轴大小的均值，mm；σ为砾石长轴大小的方差，mm2。
[0024]砾石含量可表示为：
[0025][0026]式中：η为研究区域内砾石含量，％；M为研究区域内砾石数量；S
i
为第i个砾石的面积，m2；S为研究区域的总面积，m2。
[0027]进一步的技术方案是，所步骤4的压裂裂缝扩展的岩石变形模型为：
[0028][0029]式中：σ
n,i
,σ
s,i
分别为裂缝单元i受到的法向应力和切向应力，MPa；N表示裂缝单元总个数；B表示边界应变影响系数矩阵；D
s,j
表示裂缝单元j处的切向位移不连续量；D
n,j
表示裂缝单元j处的法向位移不连续量。
[0030]进一步的技术方案是，所述步骤5的压裂裂缝扩展的流体流动模型为：
[0031][0032][0033][0034]式中：q(s,t)为当前时刻s断面流量，m3/min；q
t
(s,t)为当前时刻s断面单位缝长的压裂液滤失速率，m2/min；A(s,t)当前时刻s断面裂缝截面积，m2；t为压裂施工时长，min；p为水力裂缝中的流体摩阻；n表水流体幂律指数；k表示流体粘度指数；h为储层的厚度，m；w表示水力裂缝宽度，m；C
t
为压裂液滤失系数，m/min
0.5
；t0为裂缝开启时间，min。
[0035]进一步的技术方案是，所述步骤6水力裂缝与砾石相交作用准则。
[0036]水力裂缝止砾的判断准则可表达为：
[0037]K
e
＜min(K
IC_C
,K
IC_G
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)
[0038]式中：K
e
为裂缝尖端等效应力强度因子，MPa
·
m
0.5
；K
IC_C
为基质与砾石之间胶结面的临界断裂韧性，MPa
·
m
0.5
；K
IC_G
为砾石的临界断裂韧性，MPa
·
m
0.5
。
[0039]水力裂缝穿过砾石的判断准则可表达为：
[0040][0041]式中：K
v1
，K
v2
表示水力裂缝沿基质与砾石之间胶结面扩展相对容易与相对较难的虚拟等效应力强度因子，其计算表达式为：
[0042][0043]式中：为水力裂缝与砾石相交时的逼近角，
°
。
[0044]一般发生绕砾的情况又可分为4种，水力裂缝沿基质与砾石之间胶结面单向绕砾的判断准则可表达为：
[0045][0046]水力裂缝沿基质与砾石之间胶结面双向绕砾的判断准则可表达为：
[0047][0048]水力裂缝沿胶结面扩展时发生单向偏移或双向偏移行为的判断准则为：
[0049][0050]式中：K
v0
为水力裂缝沿当前胶结面扩展方向的虚拟等效应力强度因子，即MPa
·
m
0.5
。
[0051]如果不满足上述偏移条件，则水力裂缝可能继续沿胶结面扩展或扩展到胶结面端本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种砾岩水力裂缝扩展方向的预测方法，其特征在于，包括以下步骤：步骤1：获取地质参数和压裂施工参数；步骤2：通过砾石形状及方位、砾石空间位置分布、砾石粒径、砾石含量生成含砾岩的几何模型；步骤3：对基质和砾石的参数分别赋值；步骤4：建立压裂裂缝扩展的岩石变形模型；步骤5：建立压裂裂缝扩展的流体流动模型；步骤6：建立水力裂缝与砾石相交作用准则；步骤7：通过改变砾石大小和砾石含量探究砾石对水力裂缝扩展的影响规律。2.根据权利要求1所述的一种砾岩水力裂缝扩展方向的预测方法，其特征在于，所述步骤2的具体过程为：在极坐标系统中，确定多边形的形状由顶点数m，极角θ
i
以及极半径R
i
分别表示为：θ
i
＝η1×
2π
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)R
i
＝A0+(2η
i
‑
1)
×
A1ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)式中：η1为0～1之间的随机数；A0为极半径R
i
的平均值，mm；A1为极半径R
i
的增加值，mm；η
i
为0～1之间的随机数；砾石空间位置分布可表示为：式中：W
min
为研究区域横坐标最小值，m；W
max
研究区域横坐标最大值，m；H
min
为研究区域纵坐标最小值，m；H
max
研究区域纵坐标最大值，m；λ1，λ2为0～1之间的随机数；砾石粒径分布的概率密度函数表示为：式中：L
a
为长轴大小，mm；u为砾石长轴大小的均值，mm；σ为砾石长轴大小的方差，mm2；砾石含量可表示为：式中：η为...
【专利技术属性】
技术研发人员：李小刚，张景强，易多，杨长鑫，杨兆中，易良平，何建冈，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。