[发明专利]基于相控的序贯协同模拟地质建模方法

权利要求书

1.一种基于相控的序贯协同模拟地质建模方法，其特征在于，包括以下步骤：

A、获取原始井数据、已进行反演或模拟的数据体及层位数据，构建相模型，设定随机模拟的次数、插值层数及分别在主变量和次变量邻域找点个数；

B、对井数据和层位数据进行插值；

C、将井数据作为主变量，数据体作为次变量，分别构建主变量和次变量的变差函数；

D、采用最小二乘法计算各沉积相克里金估值对待估值点的贡献比例，具体包括以下分步骤：

D1、设定待估值点总个数为m，各相点个数的下限，邻域找点个数下限；

D2、随机选择待估值点，判断待估值点是否满足对F相找点要求；若是则对待估值点进行邻域内搜索；若否则放弃该点，重新随机选择待估值点；

D3、判断步骤D2中进行邻域内搜索得到的点是否分属三种不同的沉积相；若是则进行下一步骤；若否则返回步骤D2；

D4、判断步骤D3中分属三种不同的沉积相的点个数是否满足预设要求；若是则保存该次搜索结果；若否则返回步骤D2；

D4、分别对分属三种不同的沉积相的搜索结果进行克里金估值；

D5、判断待估值点个数是否满足预设要求；若是则采用最小二乘法计算各相克里金估值贡献比例；若否则返回步骤D2；

E、在每个模拟平面生成随机路径；

F、根据模拟平面上的随机路径进行协同模拟。

2.如权利要求1所述的基于相控的序贯协同模拟地质建模方法，其特征在于，所述步骤D2中对待估值点进行邻域内搜索具体采用空间回环式找点方式，在空间进行上下扩充，将回环路径上的每一个点都纳入搜索范围。

3.如权利要求1所述的基于相控的序贯协同模拟地质建模方法，其特征在于，所述步骤D5中采用最小二乘法计算各相克里金估值贡献比例的计算公式具体为(α,β,γ)＝(X'^TX')^-1X'^TY'

其中，α,β,γ分别为各相克里金估值贡献比例，Y'＝Y′_est＝(y'₀,y′₁,…,y'_m)^T，X'₀,X′₁,X'₂分别为分属三种不同的沉积相的搜索结果的克里金估值向量，Y′_est为待估值点。

4.如权利要求1所述的基于相控的序贯协同模拟地质建模方法，其特征在于，所述步骤F根据模拟平面上的随机路径进行协同模拟具体包括以下分步骤：

F1、选择待估值点；

F2、寻找待估值点邻域范围内的已知点；

F3、对待估值点进行基于相控的协同克里金估值；

F4、构造高斯分布模型并得到模拟数值。

5.如权利要求4所述的基于相控的序贯协同模拟地质建模方法，其特征在于，所述步骤F2寻找待估值点邻域范围内的已知点具体为采用空间回环式找点方式，分别搜索待估值点邻 域范围内的主变量和次变量。

6.如权利要求4所述的基于相控的序贯协同模拟地质建模方法，其特征在于，所述步骤F3对待估值点进行基于相控的协同克里金估值具体为利用步骤C中主变量和次变量的变差函数计算主变量、次变量及主变量和次变量之间的变异值，得到各相的克里金估值，再按照步骤D中各沉积相克里金估值对待估值点的贡献比例对各相的克里金估值进行叠加，得到协同克里金估值。

7.如权利要求4所述的基于相控的序贯协同模拟地质建模方法，其特征在于，所述步骤F4构造高斯分布模型并得到模拟数值具体为将步骤F2中搜索得到的待估值点邻域范围内的已知点纳入条件累计分布函数中，构造高斯分布模型，将高斯分布模型的均值作为协同克里金估值，将高斯分布模型的方差作为协同克里金方差，并从高斯分布模型中随机抽取一个值作为待估值点的随机模拟值。