三维重构算法简述

系统的控制函数（特征函数）：

1、单点概率函数，即孔隙率，重构前后2个模型具有相同的孔隙率（空隙相元素个数相等）。

2、三点概率函数，对于三维模型需要分别计算沿着x轴方向、y轴、z轴方向的三点概率函数。

3、线性路径函数，与三点概率函数相似，也需要计算三维空间的三个正交方向的线性路径函数。

计算过程：

1、计算参考模型的各个控制函数。

2、生成初始模型

随机生成一个512*512*200的“0、1”矩阵，“0”表示岩石基质， “1”表示空隙。

3、计算当前模型（或称为系统）的控制函数。

4、计算系统的能量E1——是一个关于三点概率函数和线性路径函数及其导数的表达式。

5、更新系统，即随机从中选取一个0和1对换位置，形成一个新系统，计算新系统的能量E2。

6、判断是否满足系统更新条件，满足转到5，否则转到7。

系统更新条件：记系统在第K+1次搜索时的状态的接受概率为P，且

其中，T为程序中的一个控制参数，称为“温度”，那，表示随机生成么系统状态接受准则定义为：[0,1）区间的任意数。若当前系统的P满足接受准则，则接受当前系统，需在其基础上进行第k+2次搜索，即转到5，否则转到7。

7、判断是否满足降温条件，满足转到8，否则转到5。

降温条件：记

度下更新失败的次数，

值，如果 ，其中为某一温度下失败率，为该温为该温度下系统更新的总次数。设定一个大于它，那么就认为该温度下，大部分性能优异的解都已经找到，继续搜素对结果没有影响，可以对系统进行降温处理，即转到8，否则转到5。

8、系统降温。降温方案为：取>0.9。

9、判断是否满足终止条件，满足输出最优系统，否则转到5。 终止条件：记nE为系统某一个状态连续更新失败的次数，当nE>10000次时则认为系统更新下去基本上不会找到比当前性能更好的状态，此时的状态可以被认为是满足要求的全局最优解的近似解，终止系统更新，并输出最终系统状态，否则转到5。 ，为降温速率，0<<1,通常

本文来源：https://www.bwwdw.com/article/d9x1.html