1.一种基于相控非参数各向异性变差函数构建方法,其特征在于,包括:S1、工区沉积相划分,采用指示模拟对工区沉积相进行划分;
S2、构建不同沉积相的变差函数,具体包括以下分步骤:S21、初始化:导入已经反演好的波阻抗数据和测井数据;
S22、相控建模:利用PETREL软件针测井数据进行相控建模,得到每个点对应的标签phase_label;
S23、对角度进行等分,设置角度区间数量为M等分,记等分M=degree_num,间隔记为θ,θ=360/M,则第n个角度区间对应斜率为[arctan((n-1)·θ),arctan(n·θ));
S24、随机选取点对,将不同标签下的点进行分类,任意两点随机组合成点对,将每组点对都进行标号;
S25、计算所选取点对的斜率,根据其取值将该点对放入对应的角度区间;
S26、对各的沉积相分别重复步骤S33至步骤S35,直到各沉积相的各个角度区间都有足够的采样点;
S27、根据不同沉积相下的采样数据分别计算变差函数;
S3、拟合变差函数参数求解,包括:
S31、对候选解{x1,x2,…}的每一变量xi用字长为N的二进制码串{bNbN-1…b1b0}进行编码,根据下面的公式进行解码:
其中,b∈{0,1},j=1,2…N,bN-1为最高位,b0为最低位,变量xi的左边界为实数值ximin,右边界为实数值ximax,z表示二进制码串对应的十进制整数值的左边界为实数值;
S32、将待拟合参数转化为有向图的形式;
S33、采用指数模型进行参数拟合。
2.根据权利要求1所述的一种基于相控非参数各向异性变差函数构建方法,其特征在于,所述步骤S1之前还包括:指示变换,具体为:根据门限值,把连续确定的原始数据离散成布尔量0或1。
3.根据权利要求1所述的一种基于相控非参数各向异性变差函数构建方法,其特征在于,所述步骤S32具体为:定义有向图G=(C,L),其中顶点集C为
其中,vs为起始点,顶点 和 分别表示二进制码串中位bj取值为0和1的状态,j=1,2…N,c0(vs), 均表示顶点集C中的元素,e=1,2,…,2N。
4.根据权利要求1所述的一种基于相控非参数各向异性变差函数构建方法,其特征在于,所述步骤S33具体包括以下分步骤:S331、nc=0,各τij和Δτij的初始化,将m个蚂蚁置于n个顶点上;
其中,nc为迭代步数或搜索次数,τij为t时刻在i,j连线上残留的信息量, 为蚂蚁k在边弧(i,j)上留下的单位长度轨迹信息素数量;
S332、将各蚂蚁的初始出发点置于当前解集中,对每个蚂蚁k(k=1,…,m),按概率 移至下一顶点j;将顶点j置于当前解集;
S333、计算各蚂蚁的目标函数值Zk,记录使得目标函数值Zk最小的解;k=1,…m;
S334、设Path*(t)为第t搜索周期内的最佳路径,该最佳路径对应的目标函数值为f*(t),边弧(i,j)中顶点i对应候选解的第Κ位,则蚁群搜索的信息素按照下式更新
其中,f*(t+1)表示第t+1搜索周期内的最佳路径对应的目标函数值,τij(t,k)表示t时刻的信息素,τij(t+1,k)表示t+1时刻的信息素,L为候选解的二进制编码的编码长度,为正整数,k表示第k个蚂蚁;
S335、对各弧边(i,j),置Δτij=0;nc=nc+1;
S336、若nc<预定的迭代次数,转步骤S332。