一种基于逆变换采样的非均质性储存层渗透率获取方法

1.一种基于逆变换采样的非均质性储存层渗透率获取方法，其特征在于，包括以下步骤：S1、结构性网格构建步骤：获取当前研究区域的地质数据，根据地质数据，构建研究区域的三维空间直角坐标系，并确定结构性网格尺寸、数量和各方向上的相关距离长度，构建由各结构性网格之间的距离计算得到的相关矩阵；S2、随机场生成步骤：确定随机种子数，利用S1中得到的相关矩阵分解出下三角矩阵，从而计算生成随机场；S3、随机场转化步骤：获取当前研究区域的渗透率均值和变异系数，将S2中生成的随机场转化为渗透率随机场；S4、逆变换步骤：对每个结构性网格进行逆变换抽样，获得逆变换后的渗透率随机场；S5、空间分布导入步骤：将S4中得到的逆变换后的渗透率随机场导入非结构性网格获得渗透率数据；S4中具体内容包括：S401获取渗透率分布的概率密度函数f(x)；式中，x为自变量，此时为渗透率；μ'为对数正态分布均值；σ'为对数正态分布方差；S402计算概率密度函数f(x)的累积分布函数F(x)得到概率值随渗透率变化的函数；式中，erf为误差函数，此时自变量为渗透率，因变量为概率值；S403计算逆累积分布函数F -1 (x)，得到渗透率随概率值变化的函数；式中，κ'为最终所得非结构性渗透率场；erfinv为逆误差函数；此时因变量为渗透率，自变量为概率值且通过u(x,y,z)获得，u(x,y,z)为一个三维空间内0～1均匀分布的随机函数，通过COMSOL软件直接生成，在空间坐标系各节点中均匀随机抽取对应的概率值。 2.根据权利要求1所述的一种基于逆变换采样的非均质性储存层渗透率获取方法，其特征在于，S1中构建相关矩阵：L·L T ＝ρ (2)其中，ρ为m×m规格的相关矩阵，由每个网格与其他网格之间的空间距离计算得到，其中m为网格总数量；Δx、Δy和Δz分别为三个方向上不同网格之间的距离，δ x 、δ y 和δ z 分别为三个方向上的相关距离，由不同岩性地层分别设置，与相距的网格数和网格尺寸有关；L和L T 分别为相关矩阵ρ通过Cholesky分解得到的下三角矩阵及其转置矩阵。 3.根据权利要求1所述的一种基于逆变换采样的非均质性储存层渗透率获取方法，其特征在于，S2中利用一个随机标准正态序列初步生成随机场：式中，G为初步生成的随机场，Z为随机标准正态序列。 4.根据权利要求1所述的一种基于逆变换采样的非均质性储存层渗透率获取方法，其特征在于，S3中变异系数为渗透率标准差与均值的比值，利用下列式子将随机场转化为渗透率随机场：κ＝exp(μ+σ·G)     (4)式中，κ为生成的结构性网格可用的渗透率场；μ为对数正态分布均值；σ为对数正态分布方差；μ N 为整个模拟区域的渗透率均值；σ N 为地区渗透率数值的标准差。