论文--采油过程中产砂趋势预测及采油最佳压力降选择的3维Abaqus数值分析.docxVIP

采油过程中产砂趋势预测及采油最佳压力降选择的3维Abaqus数值分析 摘要：本文简要介绍了石油生产过程中产砂分析及最佳采油压力降选择的数值计算原理、流程、及应用实例。针对基于塑性的产砂分析模型，分析了相关的建模过程中遇到的问题，提出了相应的解决措施。对有关的塑性力学概念和力学量进行了简要阐述和使用分析。 关键词：产砂，油层，塑性，渗流，完井 引言： 一般地说，采油过程中产砂是由于地层的地应力大、而岩石强度小。如果岩石强度足够强，或者地应力足够小，那么不会有产砂的问题。地应力大而岩石强度小的时候，地层的力学行为特征之一就是岩石的塑性变形，而且很多的室内试验及工程经验都证明了采油过程中产砂和岩石塑性变形有密切关系[1,2,3]。在过去的25年里，对基于塑性力学理论框架的采油过程中产砂预测方法的研究及其应用有很多。文献[4]介绍了使用耦合的储层岩石力学模型对产砂的起始与等效塑性应变的关系进行研究的成果。文献[5] 报道了采用一系列的渗流-材料参数的函数建立的解析模型进行产砂预测解析求解的成果。文献[6]使用霍克-布朗屈服准则而不是摩尔-库仑准则进行解析计算产砂问题的成果。文献[7]介绍的‘等效塑性应变’和‘米泽思应力’是塑性力学分析的两个经典概念，本文的产砂分析将会用到。本文的方法是基于塑性的采油过程中产砂分析方法。 数值分析中采用的理论体系属于‘孔隙弹塑性固结分析’的范畴。这个固结分析主要关注井孔周围（裸眼完井）和射孔周围（套管完井）的渗流和变形行为。近年来有若干相关研究。文献[8][9][10][11]报道了采用数值计算手段对产砂及相关的不完全固结地层侵蚀的研究及应用实例。文献[12] 报道了他们做的2维数值计算，他们分析了裸眼完井和套管完井两种情况下的产砂及侵蚀情况。文献[13]报道了采用解析和数值两种手段、对裸眼及套管两种完井方式下的产砂进行分析的结果。文献 [14] 报道了采用文献[13]介绍的3维数值方法进行产砂分析的应用实例。 本文的目标是简要介绍采油过程中产砂分析及最佳压力降选择的数值计算原理、流程。这个最佳采油、压力降是‘能够保证生产过程中不出现明显产砂现象的最大压力降’的数值。此外，数值分析结果还将给出在指定的压力降数值下生产过程中产砂的严重程度（按照交通红绿灯的形式，分为绿、黄、红三种），给之后的完井设计提供理论指导和参考依据。 计算原理： 本课题采用的方法是‘基于塑性分析的3维数值方法’，也就是在产砂趋势预测时，首先先进行孔隙弹塑性有限元分析，根据模型的塑性变形的数值结果的特性，得出结论：产砂程度轻微、或者中等、或者严重。之后，在此基础上，选择最合适的压力降、以及相应的完井设计。 这里所说的‘塑性变形的数值结果的特性’，包括两个方面：一个是数值结果中塑性变形的最大值：这个塑性变形越大，产砂趋势越严重；若没有塑性变形，则完全没有产砂风险。另一个是塑性变形沿井壁或者射孔的内表面的分布：如果塑性区在井壁周围/射孔表面形成闭合环，则认为产砂趋势严重，否则即认为不严重。 下面的分析中，采用了石油工程中关于应力正负号的惯例，即压应力为正。所谓最小水平地应力是指最小水平压应力。 当数值结果中，压力降引起的塑性变形的最大值达到或者超过一个临界值（叫做临界塑性应变值Critical Plastic Strain CPS）的时候，这个压力降就是‘没有明显产砂风险的最大需用压力降’,也叫临界压力降Critical Drawdown Pressure CDP。超过这个CDP值的压力降将会在生产过程中有明显的产砂现象。 数值计算流程： 根据上述原理，产砂趋势预测及最佳压力降选择的3维岩石力学Abaqus数值分析的数值流程如下： 进行测井数据分析，得出地层各处的岩石力学材料特性参数。需要的测井数据包括：Gama 射线，电阻率，声波包括压缩波和剪切波，密度测井，等。得出的岩石力学材料特性参数将包括弹性模量、泊松比、剪切强度、内摩擦角。这个过程包括数据整理、准备等。 初始地应力预测。将在内容1）的基础上、并结合影像测井数据及钻井报告进行。这个过程得出的结果为初始地应力三个分量的大小及方向。 建立油田区块有限元3维模型。建立的网格模型包括断层、向斜、背斜等地质结构。 建立井眼尺度的有限元3维子模型。建立的子模型将包括实际过程的细节，具体模拟钻井过程。对于套管完井，将在钻井模拟的基础上，进一步模拟射孔过程。 施加压力降，模拟生产过程中压力降产引起的模型的孔隙弹塑性变形行为。 根据塑性变形的特征（数值大小、分布范围），确定产砂风险大小、以及确定最佳压力降的数值。 下面图1给出了产砂趋势预测及最佳压力降选择的3维Abaqus分析的流程图。图中第一个阶段包括了上述的步骤1）、2）、3），共3个内容。 应用实例： 这个例子是为某国际石油公司进行的产