[发明专利]一种注入流体渗漏、窜流机理及防控方法的数值分析方法

权利要求书

1.一种注入流体渗漏、窜流机理及防控方法的数值分析方法，其特征在于，所述的注入流体渗漏、窜流机理及防控方法的数值分析方法包括：

步骤一、选取区域内典型的泥岩，进行不同含水率下的室内试验，获取不同含水率下泥岩的弹性模量、抗压强度、粘聚力、内摩擦角、蠕变特征、吸水膨胀特征力学参数及渗透率、孔隙率、饱和度渗流参数；

步骤二、建立根据套管变形数据反演套变点地层应力和变形的反演模型，根据现场统计的实际套变情况，反演套变处的地层应力分布情况和地层变形情况；

步骤三、对整个区域进行综合地质调查，清楚各地层的基本地质属性、区域内断层的分布情况及各地层的层间集合情况，在此基础上建立整个区块的Petrel地质模型，对其中的断层及层间结合进行精细建模，并转化为有限元模型；

步骤四、针对典型的高压注水方案，基于流-固耦合理论，进行整个区域内流固耦合计算，得出整个区域内的压力场分布情况、油水饱和度分布情况和整个区域内的地层变形趋势；在此基础上，选取典型的井组进行精细分析，获取井组内压力场分布、油水饱和度分布和地层变形情况；

步骤五、将每种可能的渗漏、窜层及沿断层流走情况作为一种工况，进行多工况组合计算；以压力、油水饱和度、地层变形数据的计算值与监测值或反演值之间的残差最小为优化目标，以工况组合为优化变量，建立优化模型，研究高效的优化算法，对该优化模型进行求解，寻找与监测数据差异最小的工况组合；

步骤六、针对渗漏、窜层及断层流失三种形式，给出不同流失位置和不同流失方式的防控措施；

步骤七、总结发生渗漏、窜层、沿断层流失时的地层条件、生产参数的基本规律，对不同类型针油藏发生渗漏、窜层、沿断层流失的可能性进行预测，并提出预防措施。

2.如权利要求1所述的注入流体渗漏、窜流机理及防控方法的数值分析方法，其特征在于，不同含水情况下的泥岩力学实验包括含水量对泥岩力学性质的影响的实验和含水量对泥岩蠕变特性的影响的实验；

含水量对泥岩力学性质的影响的实验的具体方法为选取不同岩心开展不同含水量下泥岩单轴抗压强度试验；

含水量对泥岩蠕变特性的影响的实验的具体方法为分别对不同岩心进行相同压力情况下进行不同含水量试样的蠕变试验。

3.如权利要求1所述的注入流体渗漏、窜流机理及防控方法的数值分析方法，其特征在于，考虑泥岩蠕变及流固耦合作用的理论模型中的砂岩弹塑性模型为：

砂岩采用弹塑性初应力法进行计算，砂岩弹塑性力学模型由平衡微分方程、几何方程、弹塑性本构方程组成，包括含屈服准则、流动法则、增量形式的应力应变关系；

(1)考虑有效应力作用的油藏岩土变形的平衡微分方程为：

σ_ij,j-(αδ_ijp)_,j+f_i＝0

(2)岩土变形的几何方程为：

ϵij=12(ui,j+uj,i)]]>

式中：ε_ij为应变张量；u为位移；

(3)本构方程，变形场本构方程采用弹塑性本构方程，本构方程增量形式为：

{dσ′_ij}＝[D_ep]{dε_ij}

式中：dσ′_ij为有效应力增量；[D_ep]为弹塑性矩阵；dε_ij为应变增量；

塑性准则采用Drucker-Prager准则，Drucker-Prager准则的形式为：F=αI1+J21/2-S]]>，式中

J₂为应力第二偏张量不变量，C为粘结力，为内摩擦角；

油藏岩土变形场的边界条件有两种：

(1)岩土骨架的表面力已知，表示为：

σ_ijL_j＝s_i(x,y,z)

式中：L_j为边界的方向导数；s_i为表面力分布函数；

(2)岩土骨架的表面位移已知，表示为：

u_i＝g_i(x,y,z)

式中：g_i为表面位移分布函数；

用有限元方法求解流固耦合岩土弹塑性变形问题的步骤为：首先将油藏空间离散化，以节点位移为基本未知量，利用虚功原理建立单元结点力和结点位移之间的平衡方程，解方程求得位移矢量后，再由几何方程求得应变，由本构方程求得应力分布。