[发明专利]基于地应力有限元的幔源油气富集区分析预测方法

权利要求书

1.基于地应力有限元的幔源油气富集区分析预测方法，其特征在于，包括以下步骤：

S1：获取被测区块的断层、层位数据以及测井数据；

S2：根据断层、层位数据以及测井数据建立三维地质实体模型为一次建模；

S3：将一次建模建立的实体模型导入到有限元软件中，并添加岩石动态物理参数；

S4：添加边界条件、划分网格，并在位移矢量云图中添加位移矢量箭头，其中位移矢量箭头的方向就是断块运动的方向；

S5：在每层的构造图高程对比：位移矢量箭头在坐标系Z轴上正方向的投影与构造图高程高的地方相符，位移矢量箭头在坐标系Z轴上负方向的投影与构造图高程低的地方相符，则表示分析结果准确，进入步骤S6，否则回到步骤S1；

S6：压应力作为地应力，输出地应力云图，在地应力云图中，深部流体充注沿着深断裂地应力低值区脉冲式充注，将被充注的浅层低地应力区分割出来，作为建立二次模型的依据；

S7：二次建模，按S3-S5的分析流程分析；

S8：压应力作为地应力，输出地应力云图，其中浅断裂模型地应力低值位置为油气富集区。

2.根据权利要求1所述的基于地应力有限元的幔源油气富集区分析预测方法，其特征在于，一次建模范围上至地表下至莫霍面，二次建模范围上至地表下至结晶基底。

3.根据权利要求1所述的基于地应力有限元的幔源油气富集区分析预测方法，其特征在于，S1中断层、层位数据通过地震解释、重、磁、电数据解释断层及地层数据获取，测井数据为被测区块所有解释过的井头、井轨迹、测井数据。

4.根据权利要求3所述的基于地应力有限元的幔源油气富集区分析预测方法，其特征在于，S2具体为：将S1得到的被测区块的断层、层位数据导入地质建模软件，生成断层面和地层面，反复清除异常点后加入测井数据，生成地质构造图；将形成的地质构造图通过地质建模软件转换为DXF格式，通过人机交互导入三维实体建模软件形成特征小平面，再通过逆向设计建立三维地质实体模型。

5.根据权利要求3所述的基于地应力有限元的幔源油气富集区分析预测方法，其特征在于，S3中所述岩石动态物理参数为垂深、密度、泊松比、杨氏模量，岩石动态物理参数数据的来源于测井数据。

6.根据权利要求5所述的基于地应力有限元的幔源油气富集区分析预测方法，其特征在于，S3还包括：

S31：采用测井数据和纵波时差计算横波时差Δt_s；

S32：通过纵波时差及横波时差计算泊松比μ；

S33：通过测井数据中、纵波时差及横波时差计算杨氏模量E；

S34：线性回归，将被测区块中各井的垂深对应的岩石动态物理参数做两个自变量和一个因变量的线性回归，两个自变量为井坐标，一个因变量为岩石动态物理参数，将回归公式整理在相应的垂深后，该数据作为有限元分析三维地质模型的岩石动态物理参数。

7.根据权利要求6所述的基于地应力有限元的幔源油气富集区分析预测方法，其特征在于，S31中横波计算公式如下：

式中：

Δt_p为岩石的纵波时差，来源校正后的解释测井数据；

ρ为岩石密度，来源于测井数据；

e为是自然对数的底数，值为e＝2.718281828459。

8.根据权利要求6所述的基于地应力有限元的幔源油气富集区分析预测方法，其特征在于，S32中泊松比计算公式如下：

式中：

Δt_p为岩石的纵波时差，来源校正后的解释测井数据；

Δt_s为岩石横波时差。