[发明专利]一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法

权利要求书

1.一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法，其特征在于，包括以下步骤：

S1、基于假设的基质压力在时间步内预估酸液向天然裂缝、基质的总滤失量，计算裂缝扩展网格内的平均滤失速度；

S2、根据步骤S1中的平均滤失速度，通过水力裂缝扩展模型计算水力裂缝内流体压力，进一步计算水力裂缝几何形态；

S3、根据步骤S2中获得的水力裂缝几何形态以及水力裂缝内流体压力，计算基质内的压力分布；

S4、将步骤S3中计算所得的基质内压力分布与步骤S1中假设的基质内压力分布进行对比，如不符合误差预期，则将步骤S3中计算所得的基质内压力分布作为假设压力，带入步骤S1中，重复步骤步骤S1-S4，直至压力差值小于误差阈值；

S5、利用步骤S1-S4获得的水力裂缝几何形态和水力裂缝、天然裂缝内的流体压力，计算天然裂缝及水力裂缝内的流体流动速度；

S6、根据步骤S5中的流体流动速度，计算天然裂缝及水力裂缝内的酸液浓度分布；

S7、根据步骤S6中的酸液浓度分布，计算酸液对水力裂缝、天然裂缝的酸刻蚀宽度，并更新天然裂缝和水力裂缝的几何形态；

S8、根据计算获得的酸液浓度分布及平均滤失速度，计算酸蚀蚓孔扩展速率，进一步计算酸蚀蚓孔长度；基于计算的酸蚀蚓孔长度，更新基质岩体渗透率分布；

S9、将计算获得的天然裂缝、水力裂缝几何参数、水力裂缝、天然裂缝内流体压力、基质渗透率参数作为初始参数带入步骤S1，重复步骤S1-S8，直至计算时间步结束。

2.根据权利要求1所述的一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法，其特征在于，步骤S1中，酸液向天然裂缝、基质的总滤失量和裂缝扩展网格内的平均滤失速度方程分别为：

1)、向天然裂缝的滤失量：

所述天然裂缝内流体压力p_nf为：

2)、向基质的滤失量：

3)、平均滤失速度

式中，q_inf,i为时间步内单元格内滤失进入第i条天然裂缝的酸液量，m³；t为计算时间步，s；h_pay为储层厚度，m；L_nf为天然裂缝长度，m；u_lnf为水力裂缝向天然裂缝内的流体滤失速度，m/s；x为油藏长度方向坐标，m；y为油藏宽度方向位置坐标，m；z为油藏高度方向位置坐标，m；t_p为水力裂缝扩展时间步，s；w_nf为天然裂缝宽度，m；v_lm为时间步内单元格向基质的酸液滤失速度，m/s；k_mx，k_my分别为基质在油藏x、y方向油藏渗透率分布，m²；μ为流体粘度，Pa·s；p_hf为水力裂缝内流体压力，Pa；p_m为基质内流体压力，Pa；q_lm为时间步内单元格内滤失基质的酸液量，m³；h_hf为水力裂缝缝高，m；u_L为时间步内单元格内酸液的平均滤失速度，m/s；A_m为单元格水力裂缝壁面面积，m²。

3.根据权利要求1所述的一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法，其特征在于，步骤S2中，所述水力裂缝内流体压力及几何形态由以下方程迭代计算获得：

式中，A_hf为水力裂缝截面积，m²；h_hf为水力裂缝高度，m；p_hf水力裂缝内的流体压力,Pa；w_hf为水力裂缝宽度，m；h_i第i层内的水力裂缝高度，m；G为储层岩石剪切模量，Pa；σ_i第i层内的地应力，Pa；μ流体粘度，Pa·s；u_L为时间步内单元格内酸液的平均滤失速度，m/s；t_p为水力裂缝扩展时间步，s。