CN116525021A - 一种基于改进双尺度模型的碳酸盐岩酸化模拟方法 - Google Patents

Abstract

本发明涉及一种基于改进双尺度模型的碳酸盐岩酸化模拟方法，步骤如下：(1)通过旋转岩盘实验，获得实验条件下酸岩反应速率与酸液浓度参数；(2)通过岩心物理实验，获得实验条件下目标储层岩石信息；(3)利用达西尺度模型构建酸液在碳酸盐岩储层多孔介质内的流动反应方程；(4)构建多级酸岩反应方程；(5)构建酸蚀过程中碳酸盐岩储层结构变化方程；(6)对不同级数碳酸盐岩酸岩反应进行数值模拟。本发明通过函数分段拟合的方法降低了经典酸化双尺度模型由于反应级数近似处理所带来的误差。提高了酸岩反应级数不为1时的酸化理论模型计算准确度，提高了碳酸盐岩酸化现场施工的精细化程度。

Description

一种基于改进双尺度模型的碳酸盐岩酸化模拟方法

技术领域

本发明涉及一种基于改进双尺度模型的碳酸盐岩酸化模拟方法，属于基质酸化数值模拟技术领域。

背景技术

我国碳酸盐岩储层油气资源量丰富，截至2021年，我国累计探明碳酸盐岩油藏地质储量达40.66×10⁸t。碳酸盐岩储层增产改造措施中体系较为成熟且应用较为广泛的一类方法便是酸化增产，通过将酸液注入近井地带的地层，与储层胶结物质和岩石碎屑等堵塞物发生化学反应，形成不同形态的酸蚀蚓孔以提高地层渗透率，从而改善油气藏生产能力。酸岩反应动力学的计算结果可以预测储层改造效果与最优施工条件，故酸岩反应动力学是酸化设计的研究重点。

反应级数是酸岩反应动力学中计算反应速率的重要参数。在酸化数值模拟中，经典双尺度模型默认酸岩反应级数为一级，而实际酸岩反应级数常不为一级，因此经典模型反应速率理论值常偏离实际值，影响酸化设计的精度。通过改进经典酸化数值模型，围绕反应级数进行修正，可提高数值模型对酸岩表面反应速率的计算准确度。反应速率计算准确度的提高有助于降低酸液注入速率及酸液黏度、酸液类型等酸化施工参数的预测偏差。因此，考虑反应级数的模型修正是提升酸化数值模型模拟精准性的重要渠道。

数值模拟是研究酸岩反应动力学的主要方法之一，目前常用的数值模拟模型为双尺度连续介质模型(Two-Scale Continuum Model)。经典双尺度连续介质模型计算酸岩反应速率时将反应级数m近似等于1，如下式所示

R(C_s)＝k_sC_s ^m

式中，R(C_s)为酸岩反应速率，mol/(m²·s)；k_s为表面反应速率常数，mol^1-m·m^{3m -1}·s^-1；C_s为酸液相内浓度，mol/L；m为反应级数，无因次。

然而实际酸化情况中，由于碳酸盐岩矿物组分、酸液种类等因素的复杂性，酸岩反应级数m常不为1级。因此经典双尺度模型对反应速率的计算存在误差值，进而影响后续对最优酸液注入速率和酸液种类优选的判断。但多级酸岩反应模型对反应级数m进行了修正，模拟结果更贴近真实值。如中国专利文件(公开号CN109882164A)公开了一种碳酸盐岩酸化数值模拟方法，制取碳酸盐岩岩心，利用MATLAB软件建立碳酸盐岩酸化数值模型，物理实验测定岩心密度、孔隙度、渗透率等岩石物性参数初始值，应力试验测定岩石裂缝参数初始值，结合取心与地震资料获取岩石溶洞参数初始值，最终建立碳酸盐岩酸化控制方程和辅助方程。然而此过程所使用的数值模型中将反应级数近似为1，因此模型计算酸化反应速率的准确性仍存在提升空间，有必要围绕反应级数对酸化模型进行修正。

目前大多数学者所使用的碳酸盐岩酸化模型仍将反应级数近似为1，对于反应级数对反应速率的误差影响少有提及。因此，针对酸化数值模型精细化模拟难题，亟需提出一种适用于计算不同级数酸岩反应的数值模型。结合酸化物模实验参数与修正的酸化数值模型，降低碳酸盐岩酸化模拟误差。

发明内容

针对现有技术的不足，本发明提供一种基于改进双尺度模型的碳酸盐岩酸化模拟方法，基于双尺度连续介质模型建立多级酸岩反应模型，模拟了不同反应级数碳酸盐岩酸化酸蚀过程。

本发明的技术方案如下：

一种基于改进双尺度模型的碳酸盐岩酸化模拟方法，包括步骤如下：

步骤(1)：通过旋转岩盘实验，获得实验条件下酸岩反应速率与酸液浓度参数，包括酸岩反应速率J、反应速率常数k_s、酸液黏度μ、酸液浓度C、酸溶蚀能力α等参数；

步骤(2)：通过岩心物理实验，获得实验条件下目标储层岩石信息，包括岩石渗透率K、孔隙度岩石密度ρ_s、比表面积a₀、岩石矿物组成及含量分布等参数。

步骤(3)：利用达西尺度模型构建酸液在碳酸盐岩储层多孔介质内的流动反应方程；

步骤(4)：构建多级酸岩反应方程，多级酸岩反应方程修正对象为反应级数m≠1情况下，经典双尺度模型的酸岩反应动力学项R(C_s)＝k_sC_s ^m(默认m＝1)计算出的理论R(C_s)值与实际R(C_s)值之间的误差；

步骤(5)：构建酸蚀过程中碳酸盐岩储层结构变化方程，为表征储层结构物性，引入孔隙尺度模型，孔隙尺度模型使用与局部孔隙度相关的半经验关系来表征渗透率、孔隙度、界面反应面积与其初始值的内在联系；

步骤(6)：使用有限元软件COMSOL Multiphysics对不同级数碳酸盐岩酸岩反应进行数值模拟，通过设置不同的酸液注入速度U₀、H⁺扩散系数D_m和酸液浓度C₀模拟不同酸化条件下的碳酸盐岩酸化反应进程。

根据本发明优选的，步骤(1)中，旋转岩盘实验是将切好的岩心固定在旋转岩盘旋转柄上，反应釜内的酸液旋转，在设定的时间间隔里测取酸液浓度变化，根据酸液的消耗速度获得实验条件下酸岩反应速率与酸液浓度参数，反应釜内使用酸液质量分数为5％、10％、15％、20％的盐酸依次实验，反应温度为90℃，转速为500r/min，岩心为含96％白云石碳酸盐岩。旋转岩盘实验为现有常规实验，其中使用的旋转岩盘等均为现有器件。

根据本发明优选的，步骤(2)中，岩石物理实验是通过实验测量不同岩性、物性、含油性状态下岩石物理性质的方法和技术总称。通过气测孔渗实验测量岩石渗透率、孔隙度和比表面积，称量法测量岩石密度，光谱测量仪测定矿物组成及含量分布。

根据本发明优选的，步骤(3)中，流动反应方程包括酸液流动方向控制方程式(2)、酸岩化学反应对岩石孔隙度影响方程式(3)、酸岩反应动力学源项定义方程式(4)和酸液流动对岩石孔隙度影响方程式(5)；

k_c(C_f-C_s)＝R(C_s) (3)

式中，U是达西速度矢量；K是渗透率张量；P是压力，单位为Pa；ε是孔隙度；C_f是液相酸的浓度，单位为mol·L^-1；C_s是液固界面酸的浓度，单位为mol·L^-1；D_e是有效扩散张量；k_c是局部传质系数，单位为m·s^-1；a_v为每单位体积介质可参与反应的界面面积m²·m^-3；ρ_s是固相的密度，单位为kg·m^-3；α是酸的溶解能力，定义为每摩尔反应酸溶解的固体克数，单位为g·mol^-1，反应动力学函数用R(C_s)表示，对于一级反应，R(C_s)的值取决于k_sC_s，其中k_s是表面反应速率常数，m·s^-1。

根据本发明优选的，步骤(4)中，本发明通过将曲线离散处理后，分段拟合函数的方法降低这一误差，对于酸液浓度0＜C_s≤1mol/L范围内的酸岩反应，仍使用经典双尺度模型；对于C_s＞1mol/L范围内的酸岩反应，选择回归拟合方法获取修正函数表达式。

根据本发明进一步优选的，C_s＞1mol/L范围内的酸岩反应函数曲线近似于线性函数，所以，假设C_s＞1mol/L范围的反应表达式如式(6)所示(反应级数m＝n时)

R_n′(C_s)＝k_s′_nC_s+R_n (6)

式中，R_n′(C_s)为反应级数为n时的修正反应动力学函数；k_s′_n为反应级数为n时的修正酸岩反应速率常数，m·s^-1；R_n为反应速率误差修正常数，mol·m^-2·s^-1。

根据最小二乘法，可以得知式(6)中k_s′_n表达式为

式中，R_n′_i(C_s)为反应级数为n时的第i个修正反应动力学函数数据点，i＝1,2,3,…；为数据点相内酸液浓度的总体平均值，mol·m^-3；/>为反应级数为n时的修正反应动力学函数数据点总体平均值；

求得k_s′_n后，根据待定系数法，由式(8)得到R_n值，

在得到最终表达式前，还需确定分段函数分段点C_seg，考虑到函数连续性，取两段曲线交点作为函数分段点，由此得知修正后的酸岩反应动力学表达式

利用公式R(C_s)＝k_c(C_f-C_s)进行等式联立，可求得修正的多级反应C_f-C_s关系式为

由此得到修正的多级酸岩反应双尺度连续介质模型表达式为

其中，源项source表达式为

以(11)式替换经典双尺度模型中的(5)式，由此，仅需酸岩反应级数及反应速率即可得到任意反应级数(m＝n，n≥0)下的修正酸岩反应动力学模型。

根据本发明优选的，步骤(5)中，碳酸盐岩储层结构表征渗透率、孔隙度、界面反应面积与其初始值的联系公式如下：

式中，ε是孔隙度；ε_o是初始孔隙度；K是局部渗透率，D；K_o是初始渗透率，D；r_p是平均孔隙半径，m；r_o是初始平均孔隙半径，m；a_o是初始比表面积，m²；a_v是比表面积，m²；β是孔隙结构相关常数。

根据本发明优选的，步骤(6)中，根据碳酸盐岩酸化现场实际施工情况，酸液注入速度区间选取7.5×10^-5m/s-7.5×10^-3m/s；由于碳酸盐岩酸化常用以HCl为基础的酸液体系，H⁺扩散系数D_m取值区间为1.5×10^-10m²/s-9×10^-9m²/s；酸化施工常用酸液浓度范围为5％-25％质量分数，酸液浓度C₀取1.4mol/L-6.4mol/L。

本发明的有益效果在于：

1、经典双尺度模型为简化计算将酸岩反应级数近似为1，影响了酸岩反应项的计算精度，导致与酸岩反应动力学结果相关的酸液排量与酸液种类参数出现偏差。本发明通过函数分段拟合的方法构建了一种适用于计算多级酸岩反应的酸化数值模型，降低了经典酸化双尺度模型由于反应级数近似处理所带来的误差。提高了酸岩反应级数不为1时的酸化理论模型计算准确度，提高了碳酸盐岩酸化现场施工的精细化程度。

2、本发明使用分段拟合函数方法处理多级酸岩反应动力学函数曲线，相较于需要对酸液浓度的非整数次方进行非线性微分计算的方法，本方法将一段曲线近似为多段线性函数，保证计算准确性的同时降低了计算复杂度，提高了计算效率。

附图说明

图1为本发明的技术流程图。

图2为本发明实施例1中白云石含量96％碳酸盐岩酸岩反应速率曲线与经典一级模型酸岩反应速率曲线对比图。

图3为本发明实施例1中修正反应级数后的反应速率曲线与白云石含量96％碳酸盐岩酸岩反应速率曲线对比图。

图4为本发明实施例1中96％白云岩含量碳酸盐岩不同酸液流速对应酸蚀孔隙度场。

图5为本发明实施例1中96％白云岩含量碳酸盐岩在不同H⁺扩散系数下的酸蚀孔隙度场。

图6为本发明实施例1中96％白云岩含量碳酸盐岩在不同酸液浓度下的酸蚀孔隙度场。

具体实施方式

下面通过实施例并结合附图对本发明做进一步说明，但不限于此。

实施例1：

如图1所示，本实施例提供一种基于改进双尺度模型的碳酸盐岩酸化模拟方法，包括步骤如下：

步骤(1)：通过旋转岩盘实验，获得实验条件下酸岩反应速率与酸液浓度参数，包括酸岩反应速率J、反应速率常数k_s、酸液黏度μ、酸液浓度C、酸溶蚀能力α等参数；

旋转岩盘实验是将切好的岩心固定在旋转岩盘旋转柄上，反应釜内的酸液旋转，在设定的时间间隔里测取酸液浓度变化，根据酸液的消耗速度获得实验条件下酸岩反应速率与酸液浓度参数，反应釜内使用酸液质量分数为5％、10％、15％、20％的盐酸依次实验，反应温度为90℃，转速为500r/min，岩心为含96％白云石碳酸盐岩。旋转岩盘实验为现有常规实验，其中使用的旋转岩盘等均为现有器件。

步骤(2)：通过岩心物理实验，获得实验条件下目标储层岩石信息，包括岩石渗透率K、孔隙度岩石密度ρ_s、比表面积a₀、岩石矿物组成及含量分布等参数，具体如表1所示：

表1实验参数数据表

步骤(3)：利用达西尺度模型构建酸液在碳酸盐岩储层多孔介质内的流动反应方程；流动反应方程包括酸液流动方向控制方程式(2)、酸岩化学反应对岩石孔隙度影响方程式(3)、酸岩反应动力学源项定义方程式(4)和酸液流动对岩石孔隙度影响方程式(5)；

k_c(C_f-C_s)＝R(C_s) (3)

式中，U是达西速度矢量；K是渗透率张量；P是压力，单位为Pa；ε是孔隙度；C_f是液相酸的浓度，单位为mol·L^-1；C_s是液固界面酸的浓度，单位为mol·L^-1；D_e是有效扩散张量；k_c是局部传质系数，单位为m·s^-1；a_v为每单位体积介质可参与反应的界面面积m²·m^-3；ρ_s是固相的密度，单位为kg·m^-3；α是酸的溶解能力，定义为每摩尔反应酸溶解的固体克数，单位为g·mol^-1，反应动力学函数用R(C_s)表示，对于一级反应，R(C_s)的值取决于k_sC_s，其中k_s是表面反应速率常数，m·s^-1。

步骤(4)：

通过实验可测得酸岩反应速度和浓度，以浓度的对数为横坐标，以反应速度的对数为纵坐标绘制曲线，对散点拟合直线得到碳酸盐岩与盐酸反应的反应速率表达式为：

lgJ＝lgk+mlgC (16)

式中，m为斜率，lgk为截距，J为酸岩反应速度，mol/(s·cm²)；V为反应酸液体积，L；S为反应面积，cm²；k为反应速度常数，(mol/L)^-m·mol/(s·cm²)；C为酸浓度，mol/L；m为反应级数，无因次，由式(6)得出酸岩反应级数及反应速率；

将本实施例测得的白云石含量96％碳酸盐岩酸岩反应的反应速率J和酸液浓度C取对数后得到的曲线图绘制在同一坐标系中，对散点拟合直线得到式(16)：

lgJ＝0.4259lgC-6.2685 (17)

将公式(16)与公式(17)中参数对应可以得到纯白云岩(白云石含量96％)与盐酸反应的反应速率表达式(18)：

J＝5.39×10^-7C^0.4259 (18)

由此可以得知白云石含量96％碳酸盐岩与盐酸的酸岩反应级数m为0.42级，酸岩表面反应速率k_s为5.39×10^-6m/s，白云石含量96％碳酸盐岩酸岩反应速率曲线与经典一级模型酸岩反应速率曲线对比图如图2所示。将1mol/L＜C_s≤6.02mol/L的曲线数据离散化，每隔0.01mol/L浓度值取一点作为离散点，获取拟合点数据组。由于拟合目标为一次函数曲线，因此假设拟合后函数表达式为：

R₀′_.42(C_s)＝k_s′_0.42C_s+R_0.42 (19)

根据最小二乘法，可以得知

求得k_s′_0.42＝1.332×10^-6m/s，根据待定系数法，由式(21)得到R_0.42＝5.395×10^{- 6}mol/(m²·s)，

由此得知修正后的酸岩反应动力学表达式

利用公式R(C_s)＝k_c(C_f-C_s)进行等式联立，可求得修正的m＝0.42时的C_f-C_s关系式为

由此得到修正的96％白云岩-盐酸反应双尺度连续模型表达式为

其中，源项source表达式为

修正后的反应速率曲线，如图3所示。

步骤(5)：构建酸蚀过程中碳酸盐岩储层结构变化方程，为表征储层结构物性，引入孔隙尺度模型，孔隙尺度模型使用与局部孔隙度相关的半经验关系来表征渗透率、孔隙度、界面反应面积与其初始值的内在联系，具体公式如下：

式中，ε是孔隙度；ε_o是初始孔隙度；K是局部渗透率，D；K_o是初始渗透率，D；r_p是平均孔隙半径，m；r_o是初始平均孔隙半径，m；a_o是初始比表面积，m²；a_v是比表面积，m²；β是孔隙结构相关常数。

步骤(6)：使用有限元软件COMSOL Multiphysics对不同级数碳酸盐岩酸岩反应进行数值模拟，通过设置不同的酸液注入速度U₀、H⁺扩散系数D_m和酸液浓度C₀模拟不同酸化条件下的碳酸盐岩酸化反应进程。

根据碳酸盐岩酸化现场实际施工情况，酸液注入速度区间选取7.5×10^-5m/s-7.5×10^-3m/s；由于碳酸盐岩酸化常用以HCl为基础的酸液体系，H⁺扩散系数D_m取值区间为1.5×10^-10m²/s-9×10^-9m²/s；酸化施工常用酸液浓度范围为5％-25％质量分数，酸液浓度C₀取1.4mol/L-6.4mol/L。

酸蚀孔隙度场如图4-6所示，图4(a)-图4(d)分别为酸液注入速度U₀＝7.5×10^-5m/s、U₀＝3×10^-4m/s、U₀＝3×10^-3m/s、U₀＝7.5×10^-3m/s时的酸蚀孔隙度场，由图4可知随酸液流速增加，酸岩反应溶蚀形态逐渐由锥形溶蚀到蚓孔溶蚀到均匀溶蚀变化；图5(a)-图5(d)分别为H⁺扩散系数D_m＝1.5×10^-10m²/s、D_m＝9×10^-10m²/s、D_m＝3×10^-9m²/s、D_m＝9×10^-9m²/s时的酸蚀孔隙度场，由图5可知随H⁺扩散系数增加，酸岩反应溶蚀形态逐渐由均匀溶蚀到蚓孔溶蚀到锥形溶蚀变化；图6(a)-图6(d)分别为酸液浓度C₀＝1.4mol/L、C₀＝2.4mol/L、C₀＝4.4mol/L、C₀＝6.4mol/L时的酸蚀孔隙度场，由图6可知随酸液浓度变化，酸岩反应溶蚀形态无明显变化。

最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (8)

1.一种基于改进双尺度模型的碳酸盐岩酸化模拟方法，其特征在于，包括步骤如下：

步骤(1)：通过旋转岩盘实验，获得实验条件下酸岩反应速率与酸液浓度参数，包括酸岩反应速率J、反应速率常数k_s、酸液黏度μ、酸液浓度C和酸溶蚀能力α；

步骤(2)：通过岩心物理实验，获得实验条件下目标储层岩石信息，包括岩石渗透率K、孔隙度岩石密度ρ_s、比表面积a₀、岩石矿物组成及含量分布参数；

步骤(3)：利用达西尺度模型构建酸液在碳酸盐岩储层多孔介质内的流动反应方程；

步骤(4)：构建多级酸岩反应方程，多级酸岩反应方程修正对象为反应级数m≠1情况下，经典双尺度模型的酸岩反应动力学项R(C_s)＝k_sC_s ^m计算出的理论R(Cs)值与实际R(Cs)值之间的误差；

步骤(5)：构建酸蚀过程中碳酸盐岩储层结构变化方程，为表征储层结构物性，引入孔隙尺度模型，孔隙尺度模型使用与局部孔隙度相关的半经验关系来表征渗透率、孔隙度、界面反应面积与其初始值的内在联系；

步骤(6)：对不同级数碳酸盐岩酸岩反应进行数值模拟，通过设置不同的酸液注入速度U₀、H⁺扩散系数D_m和酸液浓度C₀模拟不同酸化条件下的碳酸盐岩酸化反应进程。

2.如权利要求1所述的基于改进双尺度模型的碳酸盐岩酸化模拟方法，其特征在于，步骤(1)中，旋转岩盘实验是将切好的岩心固定在旋转岩盘旋转柄上，反应釜内的酸液旋转，在设定的时间间隔里测取酸液浓度变化，根据酸液的消耗速度获得实验条件下酸岩反应速率与酸液浓度参数，反应釜内使用酸液质量分数为5％、10％、15％、20％的盐酸依次实验，反应温度为90℃，转速为500r/min，岩心为含96％白云石碳酸盐岩。

3.如权利要求2所述的基于改进双尺度模型的碳酸盐岩酸化模拟方法，其特征在于，步骤(2)中，通过气测孔渗实验测量岩石渗透率、孔隙度和比表面积，称量法测量岩石密度，光谱测量仪测定矿物组成及含量分布。

4.如权利要求3所述的基于改进双尺度模型的碳酸盐岩酸化模拟方法，其特征在于，步骤(3)中，流动反应方程包括酸液流动方向控制方程式(2)、酸岩化学反应对岩石孔隙度影响方程式(3)、酸岩反应动力学源项定义方程式(4)和酸液流动对岩石孔隙度影响方程式(5)；

k_c(C_f-C_s)＝R(C_s) (3)

式中，U是达西速度矢量；K是渗透率张量；P是压力，单位为Pa；ε是孔隙度；C_f是液相酸的浓度，单位为mol·L^-1；C_s是液固界面酸的浓度，单位为mol·L^-1；D_e是有效扩散张量；k_c是局部传质系数，单位为m·s^-1；a_v为每单位体积介质可参与反应的界面面积m²·m^-3；ρ_s是固相的密度，单位为kg·m^-3；α是酸的溶解能力，定义为每摩尔反应酸溶解的固体克数，单位为g·mol^-1，反应动力学函数用R(C_s)表示，对于一级反应，R(C_s)的值取决于k_sC_s，其中k_s是表面反应速率常数，m·s^-1。

5.如权利要求4所述的基于改进双尺度模型的碳酸盐岩酸化模拟方法，其特征在于，步骤(4)中，对于酸液浓度0＜Cs≤1mol/L范围内的酸岩反应，仍使用经典双尺度模型；对于Cs＞1mol/L范围内的酸岩反应，选择回归拟合方法获取修正函数表达式。

6.如权利要求5所述的基于改进双尺度模型的碳酸盐岩酸化模拟方法，其特征在于，Cs＞1mol/L范围内的酸岩反应函数曲线近似于线性函数，所以，假设Cs＞1mol/L范围的反应表达式如式(6)所示

R_n′(C_s)＝k_s′_nC_s+R_n (6)

式中，R_n′(C_s)为反应级数为n时的修正反应动力学函数；k_s′_n为反应级数为n时的修正酸岩反应速率常数，m·s-1；R_n为反应速率误差修正常数，mol·m-2·s-1；

根据最小二乘法，可以得知式(6)中k_s′_n表达式为

式中，R_n′_i(C_s)为反应级数为n时的第i个修正反应动力学函数数据点，i＝1,2,3,…；为数据点相内酸液浓度的总体平均值，mol·m^-3；/>为反应级数为n时的修正反应动力学函数数据点总体平均值；

求得k_s′_n后，根据待定系数法，由式(8)得到R_n值，

在得到最终表达式前，还需确定分段函数分段点Cseg，考虑到函数连续性，取两段曲线交点作为函数分段点，由此得知修正后的酸岩反应动力学表达式

利用公式R(C_s)＝k_c(C_f-C_s)进行等式联立，可求得修正的多级反应C_f-C_s关系式为

由此得到修正的多级酸岩反应双尺度连续介质模型表达式为

其中，源项source表达式为

以(11)式替换经典双尺度模型中的(5)式，由此，仅需酸岩反应级数及反应速率即可得到任意反应级数下的修正酸岩反应动力学模型。

7.如权利要求6所述的基于改进双尺度模型的碳酸盐岩酸化模拟方法，其特征在于，步骤(5)中，碳酸盐岩储层结构表征渗透率、孔隙度、界面反应面积与其初始值的联系公式如下：

式中，ε是孔隙度；ε_o是初始孔隙度；K是局部渗透率，D；K_o是初始渗透率，D；r_p是平均孔隙半径，m；r_o是初始平均孔隙半径，m；a_o是初始比表面积，m²；a_v是比表面积，m²；β是孔隙结构相关常数。

8.如权利要求7所述的基于改进双尺度模型的碳酸盐岩酸化模拟方法，其特征在于，步骤(6)中，根据碳酸盐岩酸化现场实际施工情况，酸液注入速度区间选取7.5×10^-5m/s-7.5×10^-3m/s；H⁺扩散系数D_m取值区间为1.5×10^-10m²/s-9×10^-9m²/s；酸化施工常用酸液浓度范围为5％-25％质量分数，酸液浓度C₀取1.4mol/L-6.4mol/L。