利索能及
我要发布
收藏
专利号: 201510569783X
申请人: 电子科技大学
专利类型:发明专利
专利状态:已下证
更新日期:2025-12-01
缴费截止日期: 暂无
联系人

摘要:

权利要求书:

1.一种基于非参数各向异性变差函数的随机建模方法,其特征在于,包括以下步骤:A、对确定性反演得到P波阻抗体进行抽样计算变差函数;

B、根据步骤A中得到的变差函数,采用序贯高斯模拟方法进行随机建模。

2.如权利要求1所述的基于非参数各向异性变差函数的随机建模方法,其特征在于,所述步骤A对确定性反演得到P波阻抗体进行抽样计算变差函数,具体包括以下分步骤:A1、计算P波阻抗体的网格的最大距离,并对网格的最大距离和角度进行等分;

A2、随机选取一个网格点作为圆心,以圆心为起始点在一定角度范围的扇形区域内进行随机抽样;

A3、计算各个方向的变差参数,得到基台值和变程;

A4、根据所有角度范围对应的基台值和变程,拟合得到基台值函数和变程函数;

A5、采用局部平均加权平滑方法对步骤A4中得到的基台值函数和变程函数进行平滑处理,得到变差函数。

3.如权利要求2所述的基于非参数各向异性变差函数的随机建模方法,其特征在于,所述步骤A1中对网格的最大距离和角度进行等分,具体为:将网格的最大距离进行N等分,等分间隔记为h,将网格的角度进行M等分,等分间隔记为θ。

4.如权利要求3所述的基于非参数各向异性变差函数的随机建模方法,其特征在于,所述步骤A2随机选取一个网格点作为圆心,以圆心为起始点在一定角度范围的扇形区域内进行随机抽样,具体为:随机选取一个网格点i作为圆心,以圆心为起始点在角度范围(l-1)θ~lθ的扇形区域内抽取一定数量的点集Cil,使得抽取的点集均匀地落在扇形区域所包含的区间内。

5.如权利要求4所述的基于非参数各向异性变差函数的随机建模方法,其特征在于,所述步骤A3计算各个方向的变差参数得到基台值和变程,具体包括以下分步骤:A31、计算角度范围(l-1)θ~lθ的扇形区域内所有点与圆心i之间的欧氏距离dij及半方差vij,计算公式分别具体为:2

vij=0.5*(zi-zj),j∈Cil

其中,i,j为网格点序号,x,y为网格点坐标,z为网格点(x,y)的值;

A32、统计落在距离区间(k-1)h~kh内的所有点对(i,j),对半方差取平均值,记为vk'dij∈((k-1)h,kh)且j∈CilA33、用角度均值αl代替角度范围(l-1)θ~lθ,用区间均值dk'代替区间(k-1)h~kh,得到αl方向的基台值cl和变程al,具体为:其中,αl=(2l-1)*θ/2,dk'=(2k-1)*h/2。

6.如权利要求4所述的基于非参数各向异性变差函数的随机建模方法,其特征在于,所述步骤A5中采用局部平均加权平滑方法对基台值函数和变程函数进行平滑处理,具体为:将180°~360°的函数曲线部分置于0°~360°函数曲线前,将0°~180°的函数曲线部分置于0°~360°函数曲线后,再对组合得到的函数曲线进行局部加权平均平滑处理,局部加权平均平滑公式具体为:其中,xn为被平滑点,xn1、xn2、xn3和xn4分别为xn前后的取样点。

7.如权利要求1所述的基于非参数各向异性变差函数的随机建模方法,其特征在于,所述步骤B中采用序贯高斯模拟方法进行随机建模,具体包括以下分步骤:B1、生成一条随机访问路径;

B2、根据步骤B1中随机访问路径的顺序求解待模拟点,从而建立储层随机模型。

8.如权利要求7所述的基于非参数各向异性变差函数的随机建模方法,其特征在于,所述步骤B2中根据随机访问路径的顺序求解待模拟点,具体包括以下分步骤:B21、根据随机访问路径依次选取一个待模拟点,进行克里金估计;

B22、根据步骤B21中克里金估计得到的均值和方差,构造条件高斯分布函数,并从高斯分布中抽取样点作为该点的值。