注入水中悬浮微粒导致储层伤害网络模拟研究*

2014-06-07 05:59冯其红韩晓冬王守磊周文胜

西南石油大学学报(自然科学版) 2014年3期

关键词：孔喉喉道微粒

冯其红，韩晓冬，2，王守磊，张欣，周文胜

1.中国石油大学（华东）石油工程学院，山东青岛 266555 2.中海石油（中国）有限公司天津分公司，天津塘沽 300452 3.中海油研究总院，北京东城 100027

注入水中悬浮微粒导致储层伤害网络模拟研究*

冯其红1，韩晓冬1，2，王守磊3，张欣1，周文胜3

1.中国石油大学（华东）石油工程学院，山东青岛 266555 2.中海石油（中国）有限公司天津分公司，天津塘沽 300452 3.中海油研究总院，北京东城 100027

注入流体中的悬浮固体微粒随流体进入储层后会对储层造成伤害，导致储层渗透率降低，明确不同因素对储层伤害的影响规律对现场储层伤害的预防和治理有十分重要的意义。因此应用网络模拟方法对不同条件下储层伤害变化规律及孔喉变化规律进行了研究。模拟结果表明：随驱替的不断进行，孔喉半径总体逐渐减小，且距离注入端面越近，孔喉半径减小幅度越大；注入流量越小、注入流体内微粒浓度越大、流体黏度越小、微粒粒径越大，越有利于微粒的沉积，造成的储层伤害越严重。

水驱油藏；悬浮微粒；储层伤害；微粒捕集；网络模拟

冯其红，韩晓冬，王守磊，等．注入水中悬浮微粒导致储层伤害网络模拟研究[J]．西南石油大学学报：自然科学版，2014，36（3）：179–184．

Feng Qihong,Han Xiaodong,Wang Shoulei,et al．Network Simulation of Formation Damage Due to Suspended Particles in Injection Water[J]．Journal of Southwest Petroleum University：Science&Technology Edition，2014，36（3）：179–184．

引言

注水开发是目前国内外应用最广泛、最成熟的开采方式[1]。注水开发过程中，陆上油田主要利用地层采出水进行回注，而海上油田主要注海水进行开发[2]。注入水中富含固体颗粒、油滴及细菌等杂质，随流体进入储层后会堵塞部分孔喉，使储层渗透性降低，造成储层伤害[3-6]。故注入水水质对注水井的生产动态及注水开发效果有十分重要的影响。国内外很多学者通过岩芯驱替实验或者数值模拟对储层伤害进行了研究。研究发现，注入水中微粒浓度、微粒粒径及注入水流速等是影响储层伤害程度的主要因素[5-8]。

孔隙网络模拟通过对微观孔隙结构的三维定量描述反映孔隙结构特征，模拟流体在网络模型中的流动来研究多孔介质中的驱替过程[8-11]。本文综合考虑了注入水中微粒在孔喉内的捕集、运移等变化机理，建立了符合地层实际孔喉结构的三维网络模型，运用网络模拟方法对不同条件下的储层伤害程度及其变化规律进行研究。

1 微粒变化机理

注入水中的微粒随流体进入储层，产生沉积、捕集等现象导致孔喉堵塞或尺寸减小，导致储层伤害，渗透率降低[12]。本节分析总结了注入水驱替过程中孔喉内流体中微粒的不同微观变化机理，并对微粒的捕集机理进行定量描述。

1.1 捕集机理

随注入流体进入孔喉的固体微粒有可能在重力、静电力等作用下在孔喉壁面沉积，导致孔喉尺寸减小，甚至会直接堵塞孔喉，使储层连通性降低。本文模型中考虑的微粒捕集机理主要包括直接堵塞、壁面沉积和架桥堵塞。

直接堵塞：对于注入流体中悬浮的微粒流经孔喉时，若微粒半径大于孔喉半径则会直接将孔喉堵塞，使配位数降低，网络连通性变差[7]。

壁面沉积：在流体流动过程中，流体中悬浮的微粒会有一部分在重力等的作用下在孔喉壁面沉积，使孔喉尺寸减小。关于微粒的沉积，国内外学者提出过很多表征方法，本文采用的微粒沉积速率计算公式为[13-14]

式中：rc—单位孔喉长度下微粒捕集速率，kg/（m3·m·s）；

β—捕集系数，无因次；

r—孔喉半径，m；

dp—微粒半径，m；

cf—孔喉中流体的微粒浓度值，kg/m3；

L—孔喉长度，m；

u—孔喉内流体流速，m/s；

µ—流体黏度，Pa·s。

架桥堵塞：当微粒半径小于喉道半径时，单个微粒无法将孔喉直接完全堵塞，但多个微粒有可能通过架桥堵塞的方式将流经的喉道堵塞，降低孔隙网络的连通性。通常认为，微粒半径大于喉道半径的1/3时，则可通过架桥的方式将喉道堵塞[15]。

模拟孔喉内微粒捕集时，首先判断微粒是否会通过直接堵塞将流经孔喉堵塞。直接堵塞不会发生时，计算流经喉道的架桥堵塞临界浓度值，判断是否会发生架桥堵塞。最后计算未发生堵塞孔喉内微粒的沉积量，并更新该时间步结束时孔喉内微粒浓度以及孔喉尺寸大小。

1.2 运移机理

实验研究发现，当微粒与孔喉半径之比很大时，微粒在孔喉内运移的速率与孔喉内流体的速度差仍可以忽略不计，即可认为孔喉内微粒的流速和流体流速相同[7，10]。本文假设微粒在流体中均匀分布，且流速和流体流速相同。

2 储层伤害网络模拟方法建立

2.1 网络模型构建

网络模型研究中，可将多孔介质中的渗流空间看成由具有较大空间的孔隙和连通孔隙的喉道组成，孔隙和喉道用轴向均匀不变，截面形状为规则几何形状的柱体来代替[10]。假定孔隙和喉道的截面形状均为圆形，建立的模型大小为10×10×10，具有1 000个孔隙，喉道半径为1～11µm，孔喉比为1.5～3.0。实际储层的孔隙和喉道尺寸符合一定的分布规律。模型中孔隙的半径、长度、孔喉比均采用截断威布尔分布随机赋值，以喉道半径值的确定为例，其计算公式为

式中：

rmax，rmin—设定的喉道半径最大和最小值，m；

δ、γ—截断威布尔分布特征参数；

x—0和1之间的随机数。

采用形状因子表征孔喉的截面形状[16]

式中：G—形状因子，无因次；对于圆形截面，形状因子值为1/4π。

A—横截面积，µm2；

P—周长，m。

孔隙的内切圆半径为

式中：rp—孔隙内切圆半径，m；

ri—与孔隙相连不同喉道的半径值，m；

α—孔隙的孔喉比。

2.2 数学模型

根据多孔介质的渗流行为和研究需要，做出如下假设：（1）孔喉以及流体均等温、不可压缩；（2）流体渗流符合层流Poiseuille方程，管壁处无滑脱；（3）微粒在孔隙中的浓度分配按流入不同喉道的流量大小确定；（4）单个孔隙或者喉道内，流体中微粒的浓度相同。

对于单个孔喉，考虑微粒运移时的物质平衡方程，满足

式中：Qi，Qj—流入和流出孔喉的流体流量，m3/s；

V—孔喉体积，m3；

ci—流入孔喉内流体的微粒浓度，kg/m3；

Δt—时间步长，s。

微粒捕集后孔喉的有效半径为

式中：ρ—微粒密度，kg/m3；

reff—孔喉有效半径，m。

对于两相邻孔隙，通过其连接喉道的流量

式中：Qij—通过连接第i和第j个孔隙的喉道的流量，m3/s；

pi，pj—相邻孔隙处的压力值，Pa；

Lij—两相邻孔隙间距，m。

对于非边界处的所有孔隙节点，满足流量守恒

式中：Zi—孔隙i的配位数，即与之相连的喉道个数。

上述方程组中，各孔隙处的压力为未知量。应用高斯–赛德尔迭代法求解方程组，得到各孔隙处的压力值，进而求取不同孔喉内流体流速以及模型渗透率等参数。

3 不同条件下储层伤害变化规律模拟

3.1 含悬浮微粒注入水驱替前后储层孔喉尺寸变化

通过储层伤害机理分析可知，注入水中的悬浮微粒进入储层后，由于堵塞、沉积等作用使孔喉连通性降低是导致储层伤害的主要原因。为研究注入微粒对储层孔喉尺寸的影响，应用之前建立的网络模拟方法，得到如表1所示的结果。模拟过程中，注入水悬浮微粒浓度设置为1.0 kg/m3，定注入量驱替且流量为1.0×10−11m3/s。

表1 驱替前后孔喉半径分布对比Tab.1 Distribution of pore-throat radius before and after water flooding

由表1可知，驱替后孔喉半径分布范围有小幅变化，且孔喉半径总体呈减小趋势。分析可得，随着驱替的不断进行，由于注入水的不断注入，悬浮微粒就会持续进入模型中。随着流体中微粒的不断捕集，导致流经的孔喉尺寸逐渐减小，故模型孔喉尺寸总体呈减小趋势。

另外，将模型沿垂直流体流动方向均分为十个截面，统计驱替前后距注入端面不同位置处截面上孔喉尺寸的减小程度情况，结果如图1所示。

图1 距离入口端不同位置处截面的孔喉尺寸变化Fig.1 Variation of pore-throat size at different cross-sections which has various distance from interface

图1表明，受注入水中悬浮微粒的影响，驱替之后距离端面不同位置截面处的孔喉尺寸均出现不同程度的减小。且离注入端面越远，孔喉尺寸降低的程度越小，即储层伤害的程度越低。分析可得，注入流体在流经不同孔喉时，由于微粒的不断捕集，导致流体中微粒浓度逐渐降低，即离注入端面越远，孔喉内流体中的悬浮微粒浓度越小。故距离微粒入口越远，微粒沉积、捕集的概率越小，储层伤害的程度会越低。

3.2 注入微粒浓度对储层伤害的影响

储层伤害的主要表现就是随着驱替的不断进行，导致储层渗透率的不断降低。本文通过设置不同的注入流体中微粒浓度值，来研究微粒浓度对储层伤害的影响。模拟过程中注入端微粒浓度分别为1.0，2.5，4.0 kg/m3，模拟结果如图2所示。

图2 不同注入微粒浓度下渗透率变化规律Fig.2 Variation of permeability under different particle concentration in injection water

由图2可知，任一浓度条件下，随注入的不断进行，渗透率逐渐降低，且前期下降较快，后期逐渐变缓。注入微粒浓度越大，渗透率下降速度越快，储层伤害程度更大。分析可得，随驱替不断进行，储层渗透率逐渐降低，在定流量条件下，孔喉内的流体流速逐渐增大，越来越不利于微粒沉积。故随着驱替进行，渗透率下降曲线逐渐变缓。注入微粒浓度越大，微粒在孔喉内的捕集速度越大，孔喉尺寸减小幅度越大，故储层伤害程度越大。

3.3 微粒粒径对储层伤害的影响

注入水中的悬浮微粒粒径一般在微米级别，且符合一定的分布规律。对不同微粒粒径条件下的储层伤害进行模拟，粒径分别设置为1，2，3µm，模拟结果如图3所示。

图3 不同微粒粒径下渗透率变化规律Fig.3 Variation of permeability under different particle sizes in injection water

图3表明，微粒粒径越大，导致的储层伤害程度越大，即渗透率的下降幅度越大。分析可得，微粒粒径越大，其受重力等力的影响越大，越不容易在流体中悬浮，即越容易沉降。另外，微粒粒径越大，在流经孔喉时产生直接堵塞或者架桥堵塞的可能性更大，故其造成的储层伤害越严重。

3.4 注入流量对储层伤害的影响

对不同注入流量下的储层伤害规律进行模拟，且模拟过程中设置注入流量分别为5.0×10−12，1.0×10−11，1.5×10−11m3/s，结果如图4所示。

图4 不同注入流量下渗透率变化规律Fig.4 Variation of permeability at different injection rates

如图4所示，注入流体流量越大时，注入水中微粒导致的储层伤害程度越小，即储层渗透率降低越少。分析可知，注入流体流量越大时，孔喉内的流体流速越快，越不利于流体中微粒在孔喉壁面的沉积。故流体流量越大，储层伤害越小。

3.5 注入流体黏度对储层伤害的影响

通常情况下，流体黏度越大，黏滞阻力越大，其中的悬浮微粒稳定性越好，越不容易产生沉降。对不同注入流体黏度下的储层伤害情况进行模拟，流体黏度分别设置为1，3，5 mPa·s，模拟结果如图5所示。

图5 不同注入流体黏度下渗透率变化规律Fig.5 Variation of permeability under different viscosity of the injection fluid

在其他条件相同的情况下，注入流体黏度越大，储层渗透率的降低幅度越小，即储层伤害程度越小。分析可知，流体黏度越大，微粒在流体中越稳定，即越不容易发生沉降，导致储层伤害。故黏度越大，储层伤害程度越小。

4 结论

（1）注入水中悬浮微粒由于在孔喉中发生捕集，会造成储层伤害。定注入流量条件下，驱替前期渗透率下降较快，驱替后期逐渐变缓。

（2）由于注入微粒的捕集，随驱替不断进行，储层的孔喉尺寸逐渐减小；且距离注入端面越近的孔喉，其尺寸降低幅度越大。

（3）注入流体流量越小、注入流体中微粒浓度越大、流体黏度越小、微粒粒径越大，越有利于微粒的沉积，在相同的注入体积倍数下造成的储层伤害越严重。

[1]吴素英.长期注水冲刷储层参数变化规律及其对开发效果的影响[J].大庆石油地质与开发，2006，25（4）：35–37. Wu Suying.Variation rule of oil layer parameters washed by long-term injected water and its impact on development effect[J].Petroleum Geology&Oilfield Development in Daqing，2006，25（4）：35–37.

[2]王平，姜瑞忠，王公昌，等.低矿化度水驱研究进展及展望[J].岩性油气藏，2012，24（2）：106–110. Wang Ping，Jiang Ruizhong，Wang Gongchang，et al.Research advance and prospect of low salinity water flooding[J].Lithologic Reservoirs，2012，24（2）：106–110.

[3]胡雪滨，徐永高.油田注入水引起储层伤害的试验评价[J].江汉石油学院学报，2002，24（3）：53–55. Hu Xuebin，Xu Yonggao.Experiment evaluation of formation damage induced by water injection[J].Journal of Jianghan Petroleum Institute，2002，24（3）：53–55.

[4]张光明，汤子余，姚红星，等.注入水水质对储层的伤害[J].石油钻采工艺，2004，26（3）：46–48. Zhang Guangming，Tang Ziyu，Yao Hongxing，et al.Formation damage caused by injection water[J].Oil Drilling &Production Technology，2004，26（3）：46–48.

[5]孙丽丽，侯吉瑞，赵凤兰，等.苏丹油田注水伤害及对策研究[J].油田化学，2011，28（2）：137–140. Sun Lili，Hou Jirui，Zhao Fenglan，et al.Sudan Oil field injecting water damage and the counter measures[J].Oilfield Chemistry，2011，28（2）：137–140.

[6]李欢，王小琳，任志鹏，等.注入水水质对姬塬长8油藏储层伤害的影响研究[J].科学技术与工程，2012，12（24）：6003–6007. Li Huan，Wang Xiaolin，Ren Zhipeng，et al.The impact of injected water quality on reservoir damage of chang 8 member in Jiyuan[J].Science Technology and Engineering，2012，12（24）：6003–6007.

[7]Rege S D，Fogler H S.Network model for straining dominated particle entrapment in porous media[J].Chemical Engineering Science，1987，42（7）：1553–1564.

[8]Siqueira A G，Petrobras S A.A 3D network model of rock permeability impairment due to suspended particles in injection water[C].SPE 82232，2003.

[9]姚军，陶军，李爱芬.利用三维随机网络模型研究油水两相流动[J].石油学报，2007，28（2）：94–97，101. Yao Jun，Tao Jun，Li Aifen.Research on oil water twophase flow using 3D random network model[J].Acta Petrolei Sinica，2007，28（2）：94–97，101.

[10]冯其红，韩晓冬，王守磊，等.水驱砂岩油藏储层参数变化三维网络模拟[J].油气地质与采收率，2013，20（4）：79–82.

[11]冯其红，白军伟.驱替压力梯度对相渗曲线影响的网络模拟[J].大庆石油地质与开发，2011，30（2）：84–88. Feng Qihong，Bai Junwei.Network simulation on effect of displacement pressure gradient on relative permeability curve[J].Petroleum Geology and Oilfield Development in Daqing，2011，30（2）：84–88.

[12]尤启东，陆先亮，栾志安.疏松砂岩中微粒迁移问题的研究[J].石油勘探与开发，2004，31（6）：104–107，111.

[13]Jalel O，Jean F V.A two-dimensional network model to simulate permeability decrease under hydrodynamic effect of particle release and capture[J].Transport in Porous Media，1999，37（3）：303–325.

[14]Changhong Gao.Understanding capture of non-Brownian particles in porous media with network model[J].Asia-Pacific Journal of Chemical Engineering，2008，3（3）：298–306.

[15]Faruk C.Formation damage mechanisms and their phenomenological modeling-an overview[C].SPE 107857，2007.

[16]Valvatne P H.Predictive pore-scale modelling of multiphase flow[D].London，Imperial College London，2004.

编辑：牛静静

编辑部网址：http：//zk.swpuxb.com

Network Simulation of Formation Damage Due to Suspended Particles in Injection Water

Feng Qihong1,Han Xiaodong1，2,Wang Shoulei3,Zhang Xin1,Zhou Wensheng3
1.College of Petroleum Engineering，China University of Petroleum（East China），Qingdao，Shandong 266555，China 2.Tianjin Branch，CNOOC，Tanggu，Tianjin 300452，China 3.CNOOC Research Institute，Dongcheng，Beijing 100027，China

Particles suspended in injection water could damage the formation and decrease the permeability.And it’s of great importance to study the impact of various factors on the formation impairment for its prevention and treatment.In our study here，pore scale network modeling method is applied to the study of the formation damage patterns and the change of the porethroat radius under various conditions.In network models，different microcosmic particle variation mechanisms are taken into consideration.The results indicate that the pore-throat radius will decrease with the process of waterflooding and pore-throats that are closer to the inlet face will have a much higher decrease of their radius.Besides，the lower the flow rate is，the higher the particle concentration in injection water is，the lower the fluid viscosity and bigger particle size，it will be beneficial for the formation damage.

waterflooding reservoir；suspended particle；formation damage；particle capture；network modeling

http：//www.cnki.net/kcms/doi/10.11885/j.issn.1674-5086.2012.08.30.13.html