一种考虑聚合物沉淀及吸附的注聚井伤害计算方法技术

本发明专利技术公开了一种考虑聚合物沉淀及吸附的注聚井伤害计算方法，包括如下步骤：S1、对岩心进行室内聚合物驱伤害实验，得到伤害后孔渗参数K；S2、划分岩心网格和时间步长，构建初始孔隙度场和渗透率场；S3、求解压力场和速度场；S4、求解聚合物浓度场；S5、求解孔隙度和渗透率场；S6、若计算时间等于伤害时间，则计算结束，得出平均孔渗参数和否则，执行步骤S3～S5；并计算S7、采用NSGA

【技术实现步骤摘要】
一种考虑聚合物沉淀及吸附的注聚井伤害计算方法


[0001]本专利技术涉及油气藏提高采收率领域，特别是一种考虑聚合物沉淀及吸附的注聚井伤害计算方法。

技术介绍

[0002]聚合物驱采油作为提高油田采油率的一种重要的手段，它通过水溶性的高分子聚合物和水的混合，降低油水的流度比，扩大注入流体的波及体积和提高驱油效率，该技术广泛应用于提高采收率。在上世纪70年代，大庆油田最早探索聚合物驱油的科学机理，并在1995年后大规模应用聚合物驱油技术，在经过20多年的应用，聚合物驱已覆盖地质储量10亿吨以上，并且多年聚合物驱产量保持1000万吨规模，平均提高采收率12％。随着聚合物驱广泛使用，在采油生产过程中也出现了许多问题，如注聚压力升高、注聚速度下降这些问题，严重影响增产效果。
[0003]针对储层的注聚堵塞伤害问题，目前开展的针对注聚井伤害的数值模拟研究中，没有针对某一储层聚合物驱过程中的聚合物沉淀和吸附的参数进行优化和利用优化参数进行注聚井伤害模拟计算，无法为油田注聚前的注聚工艺和近井地带堵塞伤害后的解堵工作提供相应的科学依据。

技术实现思路

[0004]本专利技术的目的是针对上述不足，提出一种考虑聚合物沉淀及吸附的注聚井伤害计算方法，具体方案如下：
[0005]步骤S1：对某油田的岩心进行室内聚合物驱伤害评价实验，得到岩心孔渗及尺寸、聚合物性质及驱替参数。初始渗透率K0和孔隙度伤害后的渗透率K1和孔隙度室内岩心伤害时间t1、稳态下的注入流量Q1、聚合物粘度μ1、岩心长度L和半径r。
[0006]步骤S2：实验岩心网格划分，实验时间步长划分，构建初始孔隙度和渗透率场，具体如下：
[0007]r方向采用等比例划分n+2个网格节点，θ方向和z方向采用等间距划分m个和h个网格节点，室内岩心伤害时间t1划分T个时间步：
[0008]r
i+1
＝α
lg
r
i
i＝1,2,
…
,n+1(1)
[0009][0010][0011][0012]式中：r
i+1
为r方向节点位置，m；
[0013]Δθ为θ方向上的间距，m；
[0014]Δz为z方向上的间距，m；
[0015]Δt为时间步长，s。
[0016]初始孔隙度场采用高斯分布模型生成，通过Garman
‑
Kozeny半经验公式构建初始渗透率场：
[0017][0018][0019]式中：f(x)为概率密度函数，无量纲；
[0020]x为孔隙度随机变量，f；
[0021]s为孔隙度标准差，无量纲；
[0022]为平均孔隙度，f；
[0023]K为渗透率，m2；
[0024]为孔隙度，f；
[0025]β为参数，由实验室测得数据，无量纲。
[0026]步骤S3：依次求解压力场和速度场：
[0027][0028]初始条件：
[0029]P|
r,θ,z
＝0(8)
[0030]边界条件：
[0031][0032]式中：P1为岩心出口端压力，Pa；
[0033][0034]式中：u为r方向流体流速，m/s；
[0035]v为θ方向流体流速，m/s；
[0036]ξ为z方向流体流速，m/s。
[0037]步骤S4：求解聚合物浓度场：
[0038][0039]式中：为聚合物吸附量，ppm；
[0040]σ为聚合物沉积量，ppm；
[0041]C为聚合物浓度，ppm；
[0042]U为聚合物流速，m/s；
[0043]D为聚合物扩散系数，m2/s；
[0044]为梯度算子。
[0045]聚合物沉积量σ：
[0046][0047]聚合物吸附量
[0048][0049]式中：λ0为聚合物过滤系数，m
‑1；
[0050]b为模型参数，无量纲；
[0051]为最大无量纲吸附聚合物浓度，ppm。
[0052]初始条件：
[0053]C|
r,θ,z
＝0(14)
[0054]边界条件：
[0055][0056]式中：C1为岩心端面入口处的聚合物浓度，ppm；
[0057]C2为岩心端面出口处的聚合物浓度，ppm。
[0058]步骤S5：考虑聚合物沉积和吸附伤害，求解聚合物驱替过程中岩心孔隙度和渗透率场：
[0059][0060][0061]式中：α为聚合物沉淀参数，ppm
‑1；
[0062]β为聚合物吸附参数，ppm
‑1。
[0063]步骤S6：若计算时间等于室内岩心伤害时间t1，则计算结束，得出岩心的平均渗透率
和平均孔隙度否则，继续执行步骤S3～S5；计算目标函数
[0064]步骤S7：采用NSGA
‑Ⅱ
多目标遗传算法优化目标函数，得到优化后的目标函数A和B(DebK,PratapA,AgarwalS,etal.Afastandelitistmultiobjectivegeneticalgorithm:NSGA
‑
II[J].IEEEtransactionsonevolutionarycomputation,2002,6(2):182
‑
197.)。
[0065]步骤S8：选取优化后的A和B所对应的α和β，参照步骤S2～S6对注聚井伤害数值模拟计算，区别如下：
[0066]步骤S2中的网格和时间步长划分应当根据注聚井的伤害范围和注聚时间确定；
[0067]步骤S3中压力场的初始条件和边界条件变为：
[0068]初始条件：
[0069]P|
r,θ,z
＝P
e
(18)
[0070]边界条件：
[0071][0072]式中：Q
inj
为油田的注聚排量，m3/s；
[0073]P
e
为油田的储层的压力，Pa；
[0074]H为油层厚度，m。
[0075]步骤S4中聚合物浓度场的初始条件和边界条件变为：
[0076]初始条件：
[0077]C|
r,θ,z
＝0(20)
[0078]边界条件：
[0079][0080]式中：C3为注聚井的注入聚合物的浓度，ppm；
[0081]步骤S6中变为直到计算时间等于注聚井注聚时间，计算结束。
[0082]输出不同注聚时间下注聚井储层渗透率、孔隙度、聚合物浓度、储层压力等数据。
[0083]上述公式中涉及的相同符号，前后符号意义一致，标注一次后，全部通用。
[0084]一种考虑聚合物沉淀及吸附的注聚井伤害计算方法的流程如图1所示。
[0085]专利技术人发现，专利CN20本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种考虑聚合物沉淀及吸附的注聚井伤害计算方法，其特征在于，包括如下步骤：S1、取待模拟的目标储层岩心，通过对岩心进行室内聚合物驱伤害评价实验，得到岩心孔渗及尺寸、聚合物性质及驱替参数：初始渗透率K0和孔隙度伤害后的渗透率K1和孔隙度室内岩心伤害时间t1、稳态下的注入流量Q1、聚合物粘度μ1、岩心长度L和半径r；S2、实验岩心r方向等比例划分、θ和z方向等间距划分网格，实验时间等步长划分，采用高斯分布模型、Garman
‑
Kozeny半经验公式分别构建初始孔隙度和渗透率场：Kozeny半经验公式分别构建初始孔隙度和渗透率场：式中：f(x)为概率密度函数，无量纲；x为孔隙度随机变量，f；s为孔隙度标准差，无量纲；为平均孔隙度，f；K为渗透率，m2；为孔隙度，f；β为参数，由实验室测得数据，无量纲；S3、依次求解压力场和速度场；初始条件：P
r,θ,z
＝0 (4)边界条件：式中：P1为岩心出口端压力，Pa；式中：u为r方向流体流速，m/s；V为θ方向流体流速，m/s；ξ为z方向流体流速，m/s；
S4、求解聚合物浓度场；S5、考虑聚合物沉积和吸附伤害，求解聚合物驱替过程中岩心孔隙度和渗透率场；S6、若计算时间等于室内岩心伤害时间t1，则计算结束，得...

【专利技术属性】
技术研发人员：刘平礼，黄城熹，杜娟，陈祥，李骏，刘金明，王铭，王冠，李策，娄凤成，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：