[发明专利]基于旋转交错网格的双相介质地质数据获取方法和装置

权利要求书

1.基于旋转交错网格的双相介质地质数据获取方法，其特征在于，包括：

根据双相介质运动方程，确定目标区域的地质数据计算式；

使用旋转交错网格对所述地质数据计算式进行离散化处理；

将所述离散化处理后的地质数据计算式在所述目标区域所对应的模型下进行数值模 拟，以确定所述目标区域的优化地质数据，所述优化地质数据包括优化速度值和优化应力 值。

2.根据权利要求1所述的基于旋转交错网格的双相介质地质数据获取方法，其特征在 于，所述使用旋转交错网格对所述地质数据计算式进行离散化处理包括：

按照将速度和密度置于整网格点处，将弹性系数和应力置于半网格点处的方式，计算 所述地质数据计算式中，初步地质数据；所述初步地质数据包括初步速度值和初步应力值。

3.根据权利要求2所述的基于旋转交错网格的双相介质地质数据获取方法，其特征在 于，所述按照将速度和密度置于整网格点处，将弹性系数和应力置于半网格点处的方式，计 算所述地质数据计算式中，初步速度值和初步应力值包括：

将速度和密度赋值于整网格点处，将弹性系数和应力赋值于半网格点处；

使用所述整网格点沿旋转交错网格的对角线的四个应力分量，通过中心差分计算的方 式计算速度对时间的偏导数，即初步速度值；

使用所述半网格点沿旋转交错网格的对角线的四个速度分量，通过中心差分计算的方 式计算应力对时间的偏导数，即初步应力值。

4.根据权利要求1所述的基于旋转交错网格的双相介质地质数据获取方法，其特征在 于，所述将所述离散化处理后的地质数据计算式在所述目标区域所对应的模型下进行数值 模拟包括：

按照随地震波传播时间变化的顺序，依次循环迭代调整所述离散化处理后的地质数据 计算式，每次迭代后，均保存迭代计算出的速度值和应力值，以确定优化速度值和优化应力 值。

5.根据权利要求1所述的基于旋转交错网格的双相介质地质数据获取方法，其特征在 于，所述地质数据计算式为：

ρ11∂ui∂t2+ρ12∂Ui∂t2=σij′j-bij(∂Uj∂t-∂uj∂t)]]>

ρ12∂ui∂t2+ρ22∂Ui∂t2=s′i+bij(∂Uj∂t-∂uj∂t);]]>

其中，i和j分别表示水平方向x,y或垂直方向z中的一个，u_j和U_j分别是固相和流相的位 移在j方向的分量，b_ij为流体相对固体骨架运动时的耗散系数，σ_ij'j为作用在固相应力分量 在j方向上的偏导，s为作用在流体单元侧面上的应力，ρ₁₁和ρ₂₂分别表示介质单位体积内固 相和流相部分的有效质量，ρ₁₂为视质量，即流相相对固相运动时的质量耦合系数。