[发明专利]融合单井敏感性局域化的EnKF油藏辅助历史拟合方法

权利要求书

1.一种基于单井敏感性局域化EnKF的油藏自动历史拟合方法，其特征在于，所述基于单井敏感性局域化EnKF的油藏自动历史拟合方法包括：

将FMM引入到历史拟合中，通过利用油藏模型的静态参数场信息，计算油藏中三维单井动态敏感性区域；

以单井敏感性区域为基础，结合EnKF数据同化方法，建立消除梯度伪相关的新的油藏历史拟合方法FMM-EnKF；

构建基于单井动态敏感性区域的相关系数矩阵，利用FMM方法计算得到网格飞行时间及单井敏感性区域后，给定单井敏感性区域之外网格的参数与单井生产数据之间的敏感性相关系数为0；单井敏感性区域之内网格的参数与单井生产数据之间的敏感性相关系数采用Gaspari与Cohn提出的关系函数进行计算；将相关系数矩阵与EnKF的协方差矩阵进行Schur乘积减小或消除原EnKF方法中协方差矩阵的伪相关数据，进行局域化处理，消除敏感性区域之外的相关性，矫正EnKF数据同化方法的梯度，去除伪相关；

单井动态敏感性区域表征包括：

在油藏多孔介质中地层压力一维扩散方程为

令由上式进一步表示程函方程为

F(x)|▽τ(x)|＝1,x∈Rⁿ

式中：F为传播速度，τ为飞行时间，x为均空间位置的函数；k为渗透率，mD；μ为粘度，mPa·s；C_t为油藏综合压缩系数，MPa^-1；φ为孔隙度；

基于逆风差分近似求解程函方程，FMM法计算飞行时间的方程表示为

式中：F_I、F_J、F_z为网格x、y、z方向上的速度，ft·s^-0.5；Δx、Δy、Δz为x、y、z方向网格尺寸，ft；τ、τ₁、τ₂、τ₃分别表示

τ＝τ_i,j,k

τ₁＝min(τ_i-1,j,k,τ_i+1,j,k)

τ₂＝min(τ_i,j-1,k,τ_i,j+1,k)

τ₃＝min(τ_i,j,k-1,τ_i,j,k+1)；

其中，i、j、k表示网格体系中某一网格x、y、z方向位置索引；τ_i,j,k表示该网格飞行时间；τ₁表示网格体系中x方向上该网格左右两个网格飞行时间的最小值；τ₂表示网格体系中y方向上该网格左右两个网格飞行时间的最小值；τ₃表示网格体系中z方向上该网格上下两个网格飞行时间的最小值。

2.如权利要求1所述的基于单井敏感性局域化EnKF的油藏自动历史拟合方法，其特征在于，基于协方差局域化的EnKF方法更新公式为：

式中：维数与相同，为N_y×N_d；符号表示矩阵Schur乘积，是两个矩阵中的元素对应相乘；敏感性相关系数矩阵为如下形式

式中：ρ_ij表示状态向量中第i个油藏网格参数与第j个生产数据观测值之间的敏感性相关系数。

3.如权利要求1所述的基于单井敏感性局域化EnKF的油藏自动历史拟合方法，其特征在于，

利用FMM方法计算得到网格飞行时间及单井敏感性区域后，给定单井敏感性区域之外网格的参数与单井生产数据之间的敏感性相关系数为0，单井敏感性区域之内网格的参数与单井生产数据之间的敏感性相关系数采用关系函数进行计算，关系函数表达式如下：

式中：τ_max为单井敏感性区域之内网格的飞行时间最大。

4.一种实现权利要求1～3任意一项所述基于单井敏感性局域化EnKF的油藏自动历史拟合方法的信息数据处理终端。