一种基于泰勒伽辽金法的缝内颗粒运移模拟方法技术

本发明专利技术公开了一种基于泰勒伽辽金法的缝内颗粒运移模拟方法，包括以下步骤：构建水力裂缝物理模型，划分裂缝网格，赋予缝宽参数；对水力裂缝形态和缝内参数初始化或更新；建立水力裂缝扩展模型，按一定时间步长模拟裂缝扩展，计算水力裂缝形态及缝内液体流速场；基于泰勒伽辽金方法建立缝内颗粒运移模型，计算缝内颗粒浓度分布；计算一个时步，判断其是否收敛，不收敛时重复第2步至第4步，收敛时计算下一个时步直至模拟完成。本方法能够准确地模拟水力裂缝内的颗粒对流输送过程，确定裂缝内颗粒浓度分布，具有广阔的市场前景。具有广阔的市场前景。具有广阔的市场前景。

【技术实现步骤摘要】
一种基于泰勒伽辽金法的缝内颗粒运移模拟方法


[0001]本专利技术涉及水力压裂
，具体涉及一种基于泰勒伽辽金法的缝内颗粒运移模拟方法。

技术介绍

[0002]水力压裂是目前使用最为广泛的油气井增产措施，可以通过实施固相颗粒在水力裂缝内对流运输，如支撑剂在裂缝中的铺置、暂堵剂的暂堵转向、降滤失剂的缝内运移沉积以及钻孔岩屑的清除等，现场取得了良好的增产效果。因此很有必要对固相颗粒在水力裂缝内的对流运输(运移)过程进行数值模拟研究，确定出缝内颗粒的浓度分布，对压裂工艺的有效实施具有一定的指导作用。
[0003]颗粒运移模拟的方法有很多，其中使用最为广泛的是计算流动力学
‑
离散元耦合方法(CFD
‑
DEM)，这种方法通过CFD技术计算流场，使用DEM方法模拟颗粒运动，考虑了颗粒的形状、材料属性和粒径分布等因素，更加准确地描述颗粒的运动情况及与流场的相互影响。虽然基于精细网格划分的CFD
‑
DEM法的计算精度高，模拟结果较为准确，但会使整个问题更复杂，需要求解规模更大，计算的成本更大。因此相比于CFD
‑
DEM法，本文提出的泰勒伽辽金法能够避免传统伽辽金有限元处理带来的明显的非物理数值震荡，满足计算精度较高的同时计算成本更低，是在有限元体系下能较好实现颗粒运移模拟的方法。

技术实现思路

[0004]针对上述问题，本专利技术提供一种基于泰勒伽辽金法的缝内颗粒运移模拟方法,综合泰勒伽辽金方法和扩展有限元来模拟颗粒在水力裂缝内的运移输送过程，能够准确地模拟颗粒运移所涉及的对流输送过程。
[0005]本专利技术采用下述的技术方案：一种基于泰勒伽辽金法的缝内颗粒运移模拟方法，包括以下步骤：
[0006]S1：构建水力裂缝物理模型，划分裂缝网格，赋予缝宽参数。
[0007]S2：对所述步骤S1中的参数进行初始化或更新。
[0008]S3：基于步骤S1建立水力裂缝扩展模型，设定x方向为水平方向，y方向为垂直方向，按一定时间步长模拟裂缝扩展,通过裂缝扩展模型计算，获得裂缝形状、缝宽分布w
f
(x,y)和压力分布p(x,y)。
[0009]根据压力分布p(x,y)计算压裂液的x方向流动速度v
fx
和y方向流动速度v
fy
：
[0010][0011][0012]式中：μ表示压裂液粘度，Pa
·
s；表示x方向上的压力梯度、表示y方向上的
压力梯度，Pa/m。
[0013]根据颗粒上浮或沉降速度计算公式计算颗粒在裂缝内的x方向运移速度v
px
和y方向运移速度ν
py
。颗粒x方向运移速度v
px
可不进行修正，其与压裂液x方向流速v
fx
相等。颗粒y方向运移速度ν
py
等于y方向上压裂液流速与颗粒上浮或沉降速度之和：
[0014]v
py
＝v
fy
+v
up
[0015]颗粒上浮或沉降速度v
up
受颗粒密度、颗粒直径、携带液粘度的影响：
[0016][0017]式中：ρ
p
表示颗粒密度，kg/m3；d
p
表示颗粒直径，m；μ
w
表示携带液粘度，Pa
·
s。
[0018]S4：基于泰勒伽辽金方法建立缝内颗粒运移模型，计算缝内颗粒浓度分布。
[0019]S401：缝内颗粒浓度分布受水力裂缝缝宽分布、颗粒运移速度(包含水平方向运移速度v
px
和垂向运移速度ν
py
)等因素影响，其颗粒运移微分方程：
[0020][0021]式中：C表示水力裂缝内颗粒的浓度，％；t表示时间，s；表示散度算子；u为压裂液流速，m/s；w
f
为水力裂缝缝宽，m。
[0022]S402：确定泰勒伽辽金法的颗粒运移微分方程。运用泰勒伽辽金法将缝宽与颗粒浓度乘积w
f
C的泰勒级数展开与所述颗粒运移微分方程相结合，得到针对颗粒运移问题的泰勒伽辽金微分方程：
[0023][0024]式中：n+1和n分别表示第n+1和第n时刻；Δt为时间步长，s；C
n+1
和C
n
分别表示第n+1和第n时刻水力裂缝内颗粒的浓度，％。
[0025]S403：确定颗粒浓度分布方程。运用伽辽金有限元对所述泰勒伽辽金微分方程进行积分，得到积分方程，再采用传统有限元形函数来构造出颗粒浓度分布的近似表达式：
[0026][0027][0028]式中：和分别表示第n+1和第n时刻水力裂缝内网格节点处的颗粒浓度，％；和均为散度算子；u
i
和u
m
为分别对应算子和的压裂液流速，m/s；表示对应颗粒浓度的传统有限元自由度；i和I分别表示任意裂缝单元对应的第i个传统自由度编
号和总共I个传统自由度；x表示计算域内任意位置坐标；N
i
表示对应传统有限元近似项的传统形函数，无量纲。
[0029]S5：计算一个时步，判断其是否收敛，不收敛时重复S2
‑
S4，收敛时进行下一个步骤，判定公式为：
[0030][0031]S6：重复步骤S2
‑
S5，直至模拟完成，确定最终水力裂缝形态和缝内颗粒浓度分布。
[0032]步骤S5和S6通过编程进行计算，得到每一个时步的颗粒浓度分布云图。
[0033]本专利技术的有益效果是：
[0034]1、传统伽辽金法的积分方程构造过程是直接对颗粒运移输送微分方程，进行传统伽辽金积分处理，计算的颗粒浓度分布存在明显的非物理数值震荡。本专利技术采用泰勒伽辽金方法替代传统伽辽金方法可以有效避免非物理数值震荡，提高计算准确度。
[0035]2、本专利技术中的泰勒伽辽金法具有迎风性质，在有限元体系下利用该方法，能够准确地模拟水力裂缝内的颗粒对流输送过程，确定裂缝内颗粒浓度分布，该方法可用于模拟分析压裂降滤失剂、支撑剂、暂堵剂等在水力裂缝内的运移输送问题，具有广阔的市场前景。
附图说明
[0036]为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例的附图作简单地介绍，显而易见地，下面描述中的附图仅仅涉及本专利技术的一些实施例，而非对本专利技术的限制。
[0037]图1为本专利技术流程框图；
[0038]图2为第11时步颗粒浓度分布云图；
[0039]图3为第20时步颗粒浓度分布云图；
[0040]图4为第40时步颗粒浓度分布云图；
[0041]图5为第60时步颗粒浓度分布云图。
具体实施方式
[0042]为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例的附图，对本专利技术实施例的技术方案进行清楚、完整地描述。显然，所描本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于泰勒伽辽金法的缝内颗粒运移模拟方法，其特征在于，包括以下步骤：S1：构建水力裂缝物理模型，划分裂缝网格，赋予缝宽参数；S2：对水力裂缝形态和缝内参数初始化或更新；S3：基于步骤S1建立水力裂缝扩展模型，按一定时间步长模拟裂缝扩展，计算水力裂缝形态及缝内液体流速场；S4：基于泰勒伽辽金方法建立缝内颗粒运移模型，计算缝内颗粒浓度分布；S5：计算一个时步，判断其是否收敛，不收敛时重复步骤S2
‑
S4，收敛时计算下一个时步；S6：重复步骤S2
‑
S5，直至模拟完成，确定最终水力裂缝形态和缝内颗粒浓度分布。2.根据权利要求1所述的一种基于泰勒伽辽金法的缝内颗粒运移模拟方法，其特征在于，所述步骤S3计算水力裂缝形态及缝内液体流速场包括：裂缝形状、缝宽分布、压力分布、压裂液在不同方向的流动速度、颗粒的上浮或沉降速度。3.根据权利要求1所述的一种基于泰勒伽辽金法的缝内颗粒运移模拟方法，其特征在于，所述步骤S4具体缝内颗粒浓度分布计算：S401：缝内颗粒浓度分布受水力裂缝缝宽分布、压裂液流速场、颗粒上浮或沉降等因素，其颗粒运移微分方程：式中：C表示水力裂缝内颗粒的浓度，％；t表示时间，s；表示散度算子；u为压裂液流速，m/s；w
f
为水力裂缝缝宽，m。S402...

【专利技术属性】
技术研发人员：江有适，唐旺，李勇明，寇焱，贾靖，陈曦宇，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：