[发明专利]基于非参数各向异性变差函数的随机建模方法

权利要求书

1.一种基于非参数各向异性变差函数的随机建模方法，其特征在于，包括以下步骤：

A、对确定性反演得到P波阻抗体进行抽样计算变差函数，具体包括以下分步骤：

A1、计算P波阻抗体的网格的最大距离，并对网格的最大距离和角度进行等分；

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

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

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

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

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

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

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

4.如权利要求3所述的基于非参数各向异性变差函数的随机建模方法，其特征在于，所述步骤A3计算各个方向的变差参数得到基台值和变程，具体包括以下分步骤：

A31、计算角度范围(l-1)θ～lθ的扇形区域内所有点与圆心i之间的欧氏距离d_ij及半方差v_ij，计算公式分别具体为：

dij=(xi-xj)2+(yi-yj)2]]>

v_ij＝0.5*(z_i-z_j)²，j∈C_il

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

A32、统计落在距离区间(k-1)h～kh内的所有点对(i,j)，对半方差取平均值，记为v_k'

d_ij∈((k-1)h,kh)且j∈C_il

A33、用角度均值α_l代替角度范围(l-1)θ～lθ，用区间均值d_k'代替区间(k-1)h～kh，得到α_l方向的基台值c_l和变程a_l，具体为：

cl=vk′‾,k=1,2,3,...,N]]>

al=cl/vk′/dk′‾]]>

其中，α_l＝(2l-1)*θ/2，d_k'＝(2k-1)*h/2。