一种循缝找洞的酸压数值模拟结果准确性评价方法与流程

技术特征：
1.一种循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，包括以下步骤：s1，在有限元处理软件中，建立目标缝洞型油气藏储层的三维地质模型；s2，根据收集的目标缝洞型油气藏储层的地应力场数据，在有限元处理软件中建立目标缝洞型油气藏储层基于所述三维地质模型的三维应力场分布模型；s3，在所述三维应力场分布模型中，标记目标缝洞型油气藏储层的溶洞位置和油气井位置，形成含天然裂缝及溶洞的三维应力场分布模型，并将包括达西尺度模型和孔隙尺度模型的多尺度数学模型赋予所述含天然裂缝及溶洞的三维应力场分布模型，形成能够进行酸压数值模拟的缝洞三维应力场地质模型；s4，计算或获取指定产量下的实际稳定供液时间；s5，在所述能够进行酸压数值模拟的缝洞三维应力场地质模型中，进行在s4中指定产量下的实际稳定供液时间的生产施工条件相同的条件进行循缝找洞酸压数值模拟计算，并确定指定产量下的理论稳定供液时间；s6,判断循缝找洞的酸压数值模拟结果准确性；若理论稳定供液时间和实际稳定供液时间的相对误差绝对值在5％之内，则数值模拟结果准确；若该相对误差绝对值大于5％，则数值模拟结果不准确。2.根据权利要求1所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述s5中，进行不同生产施工条件下的循缝找洞酸压数值模拟计算时，输入实际生产施工条件下的对应数值后，所述能够进行酸压数值模拟的缝洞三维应力场地质模型形成主裂缝并开启天然裂缝，且主裂缝和开启的天然裂缝沟通溶洞，循缝找洞酸压数值模拟计算结束后，最终确定能够沟通的溶洞。3.根据权利要求2所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述理论稳定供液时间为酸压数值模拟计算中沟通的所有溶洞的总体积与指定产量的比值。4.根据权利要求3所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述指定产量下的实际稳定供液时间，是指在与所述酸压数值模拟计算条件相同的实际生产施工条件下，按照指定产量进行实际生产时，保持该指定产量恒定的时间。5.根据权利要求4所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述理论稳定供液时间和实际稳定供液时间的相对误差为，所述理论稳定供液时间和所述实际稳定供液时间的差值的绝对值与所述实际稳定供液时间的比值。6.根据权利要求1所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述s3中，将包括达西尺度模型和孔隙尺度模型的多尺度数学模型赋予所述三维应力场分布模型，形成能够进行酸压数值模拟的缝洞三维应力场地质模型，包括如下步骤：s3.1建立达西尺度模型：达西尺度模型用于描述厘米级至微米级多孔介质的模型，酸液在厘米级至微米级多孔介质内为达西渗流；(1)基质中的酸液流动酸液以一定速度注入地层，酸液中的氢离子在流体中的宏观运动(对流作用)和浓度梯度(扩散作用)下，由孔隙介质流体中传质到碳酸盐岩表面发生反应，进而改变地层孔隙度和渗透率，酸液在基质地层中的流动由达西定律控制：
流体压力分布由不可压缩流体的连续性方程控制：氢离子在流体中的浓度分布由对流扩散方程控制，其中对流扩散方程分两种情况，情况一为酸液未将岩石完全溶蚀(ε<1)，情况二为酸液将岩石完全溶蚀(ε＝1)；酸液未将岩石完全溶蚀(ε<1)情况下的对流扩散方程，需考虑氢离子在岩石表面的消耗和孔隙度的改变：酸液将岩石完全溶蚀(ε＝1)情况下的对流扩散方程，不需考虑氢离子在岩石表面的消耗和孔隙度的改变：式中，为达西速度矢量，m/s；k为地层渗透率，m2；μ为酸液粘度，pa
·
s；p为酸液压力，pa；ε为地层孔隙度；t为反应时间，s；c
f
为岩石孔隙中的酸液浓度，mol/m3；de为酸液扩散张量，m2/s；k
c
为酸液局部传质系数，m/s；a
v
为单位体积岩石具有的孔隙面积，m2/m3；c
s
为岩石表面的酸液浓度，mol/m3；α为单位摩尔的酸液所能溶蚀的岩石质量，kg/mol；ρ
s
为岩石密度，kg/m3；(2)裂缝中的酸液流动酸液在天然裂缝中的流动与在基质中的流动对应着不同的机理，天然裂缝中流动为自由流，基质中的流动为达西定律控制的多孔介质渗流；酸液在天然裂缝中的流动根据等效渗透率的概念将天然裂缝看作渗透率较大的区域，利用上述数学模型，研究酸液在裂缝性地层中的溶蚀现象，并将对压力场影响较大的天然裂缝进行网格加密，加快计算速度并保证计算的收敛性；s3.2建立孔隙尺度模型，用于为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持；s3.3，将构建的含天然裂缝及溶洞的三维应力场分布模型进行网格化，将所述多尺度数学模型计算得到的参数赋予所述含天然裂缝及溶洞的三维应力场分布模型的网格单元、节点上；得到网格化后及赋值属性参数后的能够进行酸化数值模拟的缝洞三维应力场地质
模型。7.根据权利要求6所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，在s3.2中孔隙尺度模型为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持中，渗透率、孔隙半径、比表面积与孔隙度的关系包括：地层渗透率、孔隙半径、比表面积大小与孔隙度直接相关，采用经验公式描述地层渗透率、孔隙半径、比表面积与孔隙度之间的关系为：式中，ε(0<ε<1)、k、r
p
、a
v
分别是孔隙度、渗透率、孔隙半径、比表面积；ε0、k0、r0、a0分别是初始的孔隙度、渗透率、孔隙半径、比表面积；β是与孔隙结构有关的常数，取β＝1；一般情况下，酸液在地层中的流动速度很小，可以看作是层流运动，可将裂缝看作较细的圆管，利用圆管层流公式来考察流量；式中，q为流量，cm3/s；δp为驱动压差，0.1mpa；d为裂缝直径，cm；a为截面积，cm2；μ为流体粘度，pa
·
s；l为裂缝长度，cm；利用达西公式计算得到的流量如下：取缝宽为0.2cm时，基于公式(6)、公式(7)计算得到天然裂缝对应的等效渗透率为k＝d2/32＝125
×
103μm2；令天然裂缝孔隙度ε
max
＝0.999(ε0＝0.05、k0＝0.32)，代入公式(5)得渗透率k：k与k是同一个数量级，即从流动阻力的角度考虑，将宽度为0.2cm的天然裂缝等效为孔隙度为0.999的基质是合理的，对于缝宽超过0.2cm的天然裂缝，由于天然裂缝相对于基质其渗透率已相当高，此时缝宽不是限制裂缝导流能力(渗透率与缝宽的乘积)的因素，所以将天然裂缝等效为缝宽0.2cm的裂缝处理即可。8.根据权利要求7所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，在s3.2中孔隙尺度模型为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持中，由于地层的沉积压实作用使孔隙结构在横、纵方向上不
同，所述扩散张量包括水平扩散张量d
ex
和垂向扩散张量d
et
；所述扩散张量与分子扩散系数dm、孔隙结构、流动速度等有关，其中佩克莱特数是表示对流与扩散比值的无量纲数；扩散张量分为水平扩散张量d
ex
和垂向扩散张量d
et
的公式为：d
ex
＝(α
os
λx
·
pe
p
)dm
ꢀꢀꢀꢀ
(10)d
et
＝(α
os
λ
t
·
pe
p
)dm
ꢀꢀꢀꢀꢀꢀ
(11)式中，α
os
是与孔隙结构相关的常数，取α
os
＝0.5；由岩心扩散实验得到λx＝0.5，λ
t
＝0.1；dm是分子扩散系数；pe
p
是佩克莱特数。9.根据权利要求4所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，在s3.2中孔隙尺度模型为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持中，传质系数k
c
的计算方法为：式中，sh
∞
为渐进sherwood数；re
p
为孔隙尺度雷诺数，sc为施密特数，

技术总结
本发明所循缝找洞的酸压数值模拟结果准确性评价方法，通过建立能够进行酸压数值模拟的缝洞型三维应力场地质模型，并在该模型中进行循缝找洞酸压数值模拟计算，明确沟通溶洞的数量，根据沟通溶洞的储集体体积计算指定产量下的理论稳定供液时间；并将理论稳定供液时间与实际生产过程中稳定供液时间进行对比，将循缝找洞酸压数值模拟的准确性从定性上升到了定量，进而提高循缝找洞酸压数值模拟的准确性。性。性。


技术研发人员：赵海洋 张俊江 方裕燕 侯帆 马馨悦 罗攀登 彭明旺 应海玲 张建军 纪成 何晓波 陈定斌
受保护的技术使用者：中国石油化工股份有限公司西北油田分公司
技术研发日：2020.11.17
技术公布日：2022/5/17