碳酸盐岩酸蚀蚓孔影响因素模拟分析研究.pdf

收稿日期：；修订日期：。作者简介：华青（），女，工程师，现从事天然气地质研究。：。文章编号：（）碳酸盐岩酸蚀蚓孔影响因素模拟分析研究华青，王昱珩，张娜，郑清，康志宏，杨荣相（中国石油西南油气田公司重庆气矿，重庆 ；重庆凯源石油天然气有限责任公司，重庆 ；西南石油大学石油与天然气工程学院，四川成都 ）摘要：酸蚀蚓孔生长及其形态的确定，是碳酸盐岩增产领域研究的基础。碳酸盐岩油气藏酸化工艺主要有基质酸化和酸压两大类。两类工艺对酸蚀蚓孔生长形态的调控方向不同。基质酸化的关键是促进酸蚀蚓孔的生长，酸压过程中则需要抑制蚓孔生产。基于川东某碳酸岩低渗气藏地质特征，建立了双尺度酸蚀蚓孔生长模型，通过径向模型和轴向模型两个角度，分析了孔隙度、渗透率分布等地质因素和酸液浓度、注酸排量、酸液黏度等工程因素对蚓孔生长形态的影响。研究发现，渗透率和孔隙度对蚓孔生长有明显影响，各工程因素对蚓孔生长的影响程度为：注入速度 酸液浓度 酸液黏度。关键词：碳酸盐岩；酸蚀蚓孔；数值模拟中图分类号：文献标识码：，（，；，；，）：，：；碳酸盐岩酸化过程中，在低于地层破裂压力条件下注入酸液，受到地层非均质影响，酸液沿着最小阻力方向滤失进入地层，形成蚓孔 。酸蚀蚓孔生长及其形态的确定，是碳酸盐岩增产领域研究的基础。基质酸化酸蚀蚓孔的数值模拟常用于预测酸岩反应后酸蚀蚓孔几何形态、酸液最佳注入量、蚓孔密度，通过注入量变化情况模拟岩石中蚓孔生长速度等与施工设计有关的参数 。关于蚓孔扩展规律的数值模拟按照时间发展可分为毛细管模型、网络模型、基于 法的模型、双重尺度连续模型等 。其中，双重尺度模型同时结合了达西尺度和孔隙尺度模型的优势，充分考虑 年 月石 油 地 质 与 工 程 第 卷第 期了渗流规律和微观孔隙结构变化规律，适用性较强。酸蚀蚓孔扩展数学模型采用笛卡尔坐标系下的基于双重尺度的蚓孔扩展模型进行计算。结合质量守恒方程、运动方程、酸岩反应方程及一些附加方程，通过求解压力场，得到速度场及酸浓度场的变化，从而通过孔隙度场表述蚓孔的扩展过程。使用常规酸体系进行二维的酸化模拟，设定不同酸岩反应速率常数表征不同方解石含量碳酸盐岩进行模拟实验。酸蚀蚓孔模型如图 所示。图 酸蚀蚓孔模拟模型示意图达西尺度模型：（）（）（）（）式中：为时间，；为岩心孔隙率，无量纲；为酸液密度，；为达西流速矢量，包含 、两个方向上的速度；为某方向的达西流速，；为岩心渗透率，；为酸液动力黏度，；为压力，；为梯度系数。化学场方程：（）（）（）（）（）（）（）式中：为孔隙度；为某时间孔隙内部酸液质量浓度，；为扩散系数，；为传质系数，；为孔隙比表面积，；为孔隙壁面酸液质量浓度，；为反应速度常数，；为溶蚀能力常数，；为岩石密度，。式（）式（）共同构成了双重尺度模型中的达西尺度模型。孔隙度尺度模型：（）（）（）槡（）（）式中：为初始渗透率，；为初始孔隙度；为原始孔喉半径，；为经验常数；为原始孔隙半径，；为原始比表面积，。式（）式（）共同构成了双重尺度模型中的孔隙尺度模型。注入端面和流出端面都采用定压边界，模型的初始条件为：，（）当 时，岩心两侧采用封闭边界，规定以下边界条件：，（），（）（），（）式中：为注入压力，；为岩心长度，；为标准大气压，；为岩心宽度，；为孔隙内部反应液相质量浓度，。双重尺度的蚓孔扩展模型，包含了关于时间变量和空间变量的连续偏导数。本文使用有限元方法进行编程计算，利用恒定（牛顿）法来逼近控制方程的偏导数，时间步长采用向后差分。酸蚀蚓孔数值模拟参数与表征参数?基于双重尺度蚓孔扩展模型，采用同一个基质孔隙度分布物理场进行计算分析。轴向模型设置岩心宽度为 ，长度为 ，注入端为短端，注入方石 油 地 质 与 工 程 年第 期向为垂直岩心表面，与注入方向平行的岩心端面封闭。径向模型设置岩心直径为 ，注入方向由中心沿径向流向四周。基质孔隙度场采用符合正态分布的随机数生成。表 为模拟中所使用的参数，基于川东某低渗碳酸盐岩气藏参数选取范围。基于同一基质孔隙度场，研究孔隙度状况、渗透率分布等地质因素和酸液浓度、注入速度、酸液黏度等工程因素条件下酸液突破岩心时蚓孔的扩展情况。?以宏观统计为基础，同时引入分形学刻画蚓孔的不规则特性，综合描述酸蚀蚓孔的形态特征。本文选择以突破孔隙体积倍数、分形维数、蚓孔径向扩展半径、次生分支蚓孔数四个参数作为酸蚀蚓孔数值模拟中蚓孔形态的分析指标，各指标的定义、符号和物理意义如表 所示。表 模拟参数名称单位数值孔隙度、渗透率 、温度 岩石密度 平均孔喉半径 平均比表面积 扩散系数 酸液浓度 、酸液黏度 、酸液密度 酸液溶解能力 酸液注入速率 、表 酸蚀蚓孔扩展规律表征参数信息参数名称符号单位定义表征酸蚀蚓孔的物理意义突破孔隙体积倍数 无因次突破岩心所用的酸液体积与岩心孔隙体积之比酸液消耗的程度分形维数 无因次标度不变量，描述分形粗糙、破碎、不规则、不光滑程度及复杂性的定量指标度量蚓孔分布于基质空间的程度和整体的复杂性、无序性次生分支蚓孔数条在主蚓孔四周形成的次生蚓孔条数发生滤失的潜在程度蚓孔径向扩展半径 蚓孔在径向方向上扩展的半径径向蚓孔发育情况在同一个温压系统、相同的用酸强度条件下，形成酸蚀蚓孔的深穿透，应使排量和酸浓度处于最优值 。而在本次实验设计中的最优组合应该是在最少的用酸强度（突破孔隙体积倍数）下实现对岩心的穿透，根据正交实验直观分析的原理，在宏观统计类表征参数中，以突破孔隙体积倍数为目标函数，其值越小代表实现穿透的酸液消耗量越小。在分形理论类表征参数中，以分形维数作为目标函数，分形维数越大，则表示蚓孔的理论形态越优 。酸液流动形成的蚓孔会减少鲜酸向深处流动，同时减少与岩石溶蚀的距离，导致酸液的大量滤失而无法波及深远地层，使得酸化改造仅局限于近井地带，对于酸压来讲蚓孔出现意味着酸液更多地滤失。因此，用次生分支蚓孔数表征发生滤失的潜在程度。不同因素对酸蚀蚓孔扩展规律的影响?渗透率的影响通过控制单一变量法，基于同一基质孔隙度场，研究渗透率分别为 、情况下，酸液突破岩心时蚓孔的扩展情况。图 为不同渗透率条件下蚓孔扩展形态，渗透率分布对蚓孔扩展影响明显。随着渗透率的增大，酸液在岩心中的扩展长度也越大，当渗透率为 时，酸液不能深入岩心内部，无法形成明显的蚓孔形态，只能在岩心表面反应，形成均匀溶蚀；当渗透率为 时，能够形成明显的蚓孔形态，但不能突破岩心；当渗透率为 时，酸液溶蚀岩心形成明显的酸蚀蚓孔，并且能够突破岩心。由此可知，储层渗透率太低，酸液的流动受阻，主要反应方式受到表面反应控制，产生均匀溶蚀，基质酸化形成蚓孔的可能性较小 ；随着渗透率的增大，酸液的滤失速度变快，导致形成的蚓孔半径和长度较大。孔隙度的影响通过控制单一变量法，基于同一基质渗透率分布，研究孔隙度分别为 、情况下，酸液深入岩心时蚓孔的扩展情况。图 为不同孔隙度条件下蚓孔扩展形态，孔隙度分布对蚓孔影响明显。随着孔隙度的增大，酸液在岩心中的扩展程度也越大，当孔隙度为 时，酸液无法在岩心内部扩展，溶蚀面呈锥面，未能形成明显的蚓孔形态；当孔隙度为 时，能够形成明显的蚓孔形态，但蚓孔深度未穿透岩心；当孔隙度为 时，能够形成明显的蚓孔形态，且蚓孔形态复杂，形华青等碳酸盐岩酸蚀蚓孔影响因素模拟分析研究 不同渗透率条件下蚓孔扩展形态 不同孔隙度条件下蚓孔扩展形态图 不同地质因素条件下蚓孔扩展形态成较为明显的“指进”现象。由此可知，储层孔隙度低，酸液的流动方向单一，主要流动方向为注入方向，形成的蚓孔形态单一；随着地层初始平均孔隙度的增大，蚓孔突破岩心需要溶解的岩石越少，酸液滤失量和滤失速度增加 ，蚓孔向四周扩展的难度降低，更容易产生形态复杂的酸蚀蚓孔。?酸蚀蚓孔数值模拟结果采用正交实验方法安排模拟次序，研究不同注入速度、酸液浓度及酸液黏度的情况下，酸液突破岩心时蚓孔的扩展情况。根据正交实验设计原理，设计 因素 水平实验（表 、表 ）。从川东地区碳酸盐岩储层的基本物性参数和实际施工参数出发，根据相似比尺原则，设计了相应的模拟指标，模拟中所使用的其他模拟参数见表 。模拟得到轴向角度酸蚀蚓孔拓展形态和径向角度酸蚀蚓孔拓展形态见图，正交实验见表 。表 正交实验考察的因素及水平序号 酸液质量分数 注入速度（）黏度（）水平 水平 水平 表 正交实验序号 酸液质量分数 注入速度（）黏度（）不同渗透率条件下轴向蚓孔扩展形态 不同孔隙度条件下径向蚓孔扩展形态图 不同工程因素条件下蚓孔扩展形态 各因素对蚓孔扩展规律影响）突破孔隙体积倍数的影响因素分析。模拟结果中突破所需酸液的孔隙体积的变化趋势为突破孔隙体积倍数随酸浓度的增加略有减少，随注入速度的增加而大幅度减少，随黏度的增加略有增加（图 ）。从各因素对突破孔隙体积倍数级差分析结果可以看到，对突破孔隙体积倍数影响程度最大的是注入速度，其次为酸液浓度，再次为黏度（图 ）。石 油 地 质 与 工 程 年第 期书书书表 正交实验结果序号酸液质量分数 注入速度（）黏度（）轴向模拟突破孔隙体积倍数分形维数 次生分支蚓孔数径向模拟蚓孔径向扩展半径 分形维数 图 突破孔隙体积倍数分析结果）轴向分形维数的影响因素分析。模拟结果中分形维数的变化趋势为分形维数随酸浓度的增加而增加的趋势，随注入速度的增加而大幅度增加，随黏度的升高小幅度增加（图 ）。从各因素对轴向模拟分形维数级差分析结果可以看到，对分形维数影响程度最大的是注入速度，其次为酸液浓度，再次为黏度（图 ）。）径向分形维数的影响因素分析。模拟结果中径向分形维数的变化趋势为分形维数随酸浓度的增加有先增加后减小的趋势，随注入速度的的增加而大幅度增加，随黏度的升高有小幅度增加后减小的趋势（图 ）。从各因素对径向模拟分形维数级差分析结果可以看到，对分形维数影响程度最大的是注入速度，其次为黏度，再次为酸液浓度（图 ）。图 轴向模拟分形维数分析结果华青等碳酸盐岩酸蚀蚓孔影响因素模拟分析研究图 径向模拟分形维数分析结果对模拟结果中各因素对不同酸蚀蚓孔形态表征参数的影响规律和影响程度的分析显示，注入排量对表征参数的影响最为显著，得出工程因素中影响蚓孔形成及扩展的主控因素为酸液的注入排量，其次为酸液的浓度，最后为黏度。其中，注入速度及酸液浓度呈相互制约的关系 ，酸液的黏度增大可以减小酸液突破时的酸量；较小的注入速度下不能形成有效地深穿透蚓孔，而较大的注入速度又会形成过多的次生蚓孔分支造成酸液的浪费，酸浓度以及黏度的提高又会使最优的注入排量值降低。结论）渗透率对酸蚀蚓孔生长情况有明显影响，储层渗透率太低，酸液的流动受阻，无法深入到岩心内部，酸液在岩心入口端产生均匀溶蚀，不能形成明显的酸蚀蚓孔；随着渗透率的增大，酸液的滤失速度变快，导致形成的蚓孔半径和长度较大。）孔隙度对蚓孔的扩展有明显的影响，储层孔隙度低，酸液的流动方向单一，主要流动方向为注入方向，形成的蚓孔形态单一；随着地层初始平均孔隙度的增大，蚓孔突破岩心需要溶解的岩石越少，酸液滤失量和滤失速度增加，蚓孔向四周扩展的难度降低，更容易产生形态复杂的酸蚀蚓孔。）基质酸化中各工程因素对酸蚀蚓孔生长的影响程度为：注入速度 酸液浓度 酸液黏度，过大的注入速度会形成过多的次生蚓孔分支增强酸液滤失的可能性。参考文献 ，（）：李小刚，陈雨松，邓庄，等 碳酸盐岩基质酸化中酸蚀蚓孔模拟回顾与展望 新疆石油地质，（）：赵立强，王润宇，刘平礼，等 考虑天然微裂缝的酸蚀蚓孔扩展数值模拟研究 油气藏评价与开发，（）：邓庄 碳酸盐岩基质酸化酸蚀蚓孔生长模拟 成都：西南石油大学，何春明，郭建春，王文耀，等 致密碳酸盐岩酸蚀蚓孔研究 断块油气田，（）：姚奕明，李双明，蔡鑫，等孔隙空间分布规律对酸蚀蚓孔扩展规律的影响 科学技术与工程，（）：杨兆中，邓庄，于梦男，等 碳酸盐岩基质酸化蚓孔生长机理实验研究 油气藏评价与开发，（）：（编辑赵川喜）石 油 地 质 与 工 程 年第 期