系统的控制函数（特征函数）：
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，且