[发明专利]三维多孔连续介质离散模型的定排量挤酸计算方法

权利要求书

1.三维多孔连续介质离散模型的定排量挤酸计算方法，其特征在于，包括以下步骤：

S1、建立初始孔隙空间分布模型，并进行网格划分；

S2、通过三维多孔连续介质离散模型计算m+1时刻的压力场，再由压力场得到流速场；所述压力场为通过迭代循环得到的真实压力分布的近似解，所述压力场的计算方法包括：

S201、假设一个初始注入压力P_guess，然后求得当前注入压力下的压力分布；

S202、通过达西公式计算入口端的酸液注入流量Q_c；

S203、迭代循环，直至计算流量Q_c与真实流量Q之差的绝对值小于给定的误差系数，所述迭代循环包括以下步骤：

S2031、如果|Q-Q_c|≤ε或者Q_c＜0，跳至步骤S2032；

否则直接修正酸液注入压力为并跳至步骤S2035；

S2032、若Q_c＜0，或|Q-Q_c|≤ε且为第一次迭代循环，则酸液注入压力P_inj修正为：

若|Q-Q_c|≤ε且不是第一次迭代循环，则酸液注入压力P_inj修正为：

S2033、将压力修正次数n加1，用Q_c替代之前的Q_old，P_old替代之前的P_old1，P_inj’替代之前的P_inj；

S2034、用P_inj’作为恒定注入压力边界，求压力分布；然后计算入口端的酸液注入流量Q_c，再跳至S2031；

S2035、以计算得到的P_inj’为目标结果，迭代循环结束；

其中，Q—给定的注入流量，m³/s；Q_c—计算出的注入流量，m³/s；ε—误差系数，无因次；n—压力修正次数，无因次；Q_old—上一次通过修正计算出的注入流量，m³/s；P_old—上一次修正得到的注入压力，MPa；P_old1—P_old之前一次修正得到的注入压力，MPa；w—权重，无因次；

S3、通过三维多孔连续介质离散模型计算m+1时刻各空间节点处的传质速度和H⁺扩散速度；

S4、通过三维多孔连续介质离散模型计算m+1时刻的酸浓度场；

S5、通过三维多孔连续介质离散模型依次计算出m+1时刻的孔隙度、渗透率、比表面、孔吼半径，并更新各空间节点的物性属性；

S6、判断：如果蚓孔突破岩心则计算完毕，否则回到步骤S2；

其中，m在三维多孔连续介质离散模型的时间长度内取任意单元时刻。