QSGS四参数随机生长法

一般来说，如果不是不可能完全描述多孔介质的微观结构是非常困难的，因为它具有复杂和随机性。人们只能获得一些基于统计的平均信息，如平均孔隙度或更好 的孔径分布。如果需要对多孔结构的全部细节进行更为 严格的处理，则必须解决此问题。 事实上，更准确地预测多孔介质的传输特性需要更详细地描述整个多孔介质 的形态，包括几何性质（如颗粒或孔形状）以及体积和 拓扑性质（如孔迂曲度和互连性）。 已经报道了几次这 样的尝试。 重建过程是一种流行的方法 再现多孔结构[。 然而，确定相关函数非常复杂。 随机当其他微观结构细节存在时，障碍物的位置是构建人造 多孔介质最简单的位置 可以忽略。 为了调整孔隙大小和连通性， Coveney等人提出了一种孔隙增长随时间模型。 通过从进一步与集群增长理论有关，我们建议 本文是一个更全面的方法，其中四个参数被确定用于控制 内部多孔 颗粒介质结构，从而形成一个称为四重结构生成集（QSGS）的集合。 这一套使我们能够生成多孔形态学特 征， 为许多真正的多孔介质的形成进程作出贡献。

clear;
max_j = 200; max_i=200; d14= 0.1; d58= d14/4; n= 0.6;
cdd=0.01; numtotal_need= （1- n） * max_i* max_j;
numsoild=0; arrgrid=0;
for i=1: max_i %第1 步 遍历所有网格 生成固相内核
   for j = 1:max_j
     if rand（ ） < cdd
       numsoild= numsoild+ 1;arrgrid（ij） = 1;
       soild（numsoild1） = i;soild（numsoild 2） = j;
       
     end % end if
   end % end j
end % end i
Tnumsoild= numsoild;
% 第2 步 从固相内核向周围8 个方向生长
while Tnumsoild< numtotal_need
  for index_soild= 1: Tnumsoild
    index_i= soild（index_soild 1）;
    index_j= soild（index_soild 2）;
     % 固相内核向右边生长 图1 中右1 方向
     if index_j< max_j
       i= index_i; j = index_j + 1; % right 1
       if （arrgrid（ i j ）== 0 & rand（ ） < d14）
          numsoild= numsoild+ 1; arrgrid（ i j ） = 1;
          soild（ numsoild 1） = i; soild（numsoild 2） = j;
       end 
     end %index_j
% 固相内核向上边生长 图1 中右2 方向
     if index_i> 1
       i= index_i-1; j = index_j; % 2
       if （arrgrid（ i j ）== 0 & rand（ ） < d14）
          numsoild= numsoild+ 1; arrgrid（ i j ） = 1;
          soild（ numsoild 1） = i; soild（numsoild 2） = j;
       end 
     end %index_j
% 固相内核向左边生长 图1 中右3 方向
     if index_j> 1
       i= index_i; j = index_j-1; % 3
       if （arrgrid（ i j ）== 0 & rand（ ） < d14）
          numsoild= numsoild+ 1; arrgrid（ i j ） = 1;
          soild（ numsoild 1） = i; soild（numsoild 2） = j;
       end 
     end %index_j
% 固相内核向下边生长 图1 中右4 方向
     if index_i< max_i
       i= index_i+1; j = index_j; % 4