[发明专利]基于流场表征的湖盆地貌定量预测方法

说明书

技术领域

本发明涉及油气勘探开发技术领域，更具体地说它是基于流场表征的湖盆地貌定量预测方法。

背景技术

在传统的地质研究中，对于湖盆地貌的表征往往采用逐层回剥法，将地层逐层的去压实，逐层的剥蚀从而恢复古地貌，这种方法往往采用的是对地质过程的反演得到研究区的古地貌，采用这种将今论古的方法恢复的地貌从而恢复沉积后沉积体的形态，这样的方法往往只能通过寻求与现代沉积相类似的沉积体来进行约束表征，这个过程是一种定性的表征，不仅繁琐而且得到的结果也是差强人意。

发明内容

本发明的目的是提供一种基于流场表征的湖盆地貌定量预测方法，可以对不同的沉积环境进行模拟分析，模拟的成本低、结果可以反复调试。

为了实现上述目的，本发明的技术方案为：基于流场表征的湖盆地貌定量预测方法，其特征在于：包括如下步骤：

步骤1：确定主要控砂地质因素：从沉积砂体的沉积动力学特征出发，开展沉积学调查研究，确定研究区沉积砂体的构造地质特征和与构造地质特征相关的沉积动力控制因素；

步骤2：根据研究区构造特征，运用地层埋藏史，对盆地进行恢复；在盆地恢复的基础上，建立盆地沉积砂体体系动力学模拟系统，分析物源、水动力条件参数和构造特征参数对砂体分布演化的影响；

首先依据Navier-Stocks方程的并行有限元算法模拟研究区的流场分布，将流场分布数据与地质资料进行对比分析，调整水动力条件参数，使研究区的流场分布符合地质条件；Navier-Stocks方程在控制方程的基础上对包括粘度、密度、重力加速度的水动力学影响因素进行表征；控制方程包括连续性方程、动量方程，

连续性方程如下：

其中：ζ为自由表面到基准面的距离，m；

d为基准面以下的水深，m；

U、V为平均速度，m/s；

Q为流量，m³/s；

动量方程如下：

其中：ζ为自由表面到基准面的距离，m；

d为基准面以下的水深，m；

U、V为平均速度，m/s；

ρ为流体密度，kg/m³；

P为压强，Pa；

H为总共的水深H＝d+ζ，m；

E为系统的蒸发量，mm；

v_V为垂直涡流粘度，m²/s；

fv、fu分别为垂直和水平方向的摩擦系数；

步骤3：按照颗粒搬运方式的不同，将沉积物分为粘性泥沙、非粘性泥沙，利用二种沉积物建立的模拟方法，得到沉积组分的模拟表征，

当沉积物为粘性泥沙时，采用扩散方程进行包括传输运移、沉积的沉积模拟，对于粘性泥沙，设定为一固定值(粘性的泥沙的沉积剥蚀学术界没有确定的值，使用时根据地质工作者的经验进行设定)，只对粘性泥沙的沉积和剥蚀条件加以约束，

扩散方程公式如下：

其中：c：沉积组分的浓度，kg/m³；

u,v,ω：水流的速度分量，m/s；

ε_s,x,ε_s,y,ε_s,z：沉积组分的扩散系数，m²/s；

ωs：沉积组分的沉降速度，m/s；

对于非粘性泥沙采用Van Rijn的沉积物的沉降函数进行沉积模拟，公式如下：

其中：s沉积组分的相对密度ρ_s/ρ_w；

D_s：沉积组分代表粒径，μm；

υ：水体的运动粘度，m²/s；

步骤4：利用步骤3得到沉积组分的模拟表征，在分析沉积物的沉积、剥蚀、路过不留三种沉积条件的判断，计算沉积物的沉积量，运用岩石的密度公式ρ＝m/v(ρ为密度、m为质量、v为体积)计算沉积物的体积，

根据计算沉积区的大小计算沉积物的沉积厚度，将沉积厚度叠加到初始底形上，作为下一时刻的沉积底形；

步骤5：重复步骤2至步骤4，直至模拟结束实现对沉积盆地的定量表征；

步骤6：在模拟的基础上，将沉积物的厚度叠加到底型上，最后投影到水平面上，从而实现湖盆地貌的定量预测。