基于改进QSGS-LBM法的各向异性重构土体的渗流模拟

为解决传统四参数随机生长法（QSGS）重构模型与真实土体情况不符的缺陷，提出一种改进QSGS法，在综合考虑各向异性、孔隙率、土体颗粒大小情况下构建各向异性多孔介质模型，采用两点自相关函数法验证模型有效性，并结合格子玻尔兹曼方法(LBM)自编Matlab程序模拟细观渗流过程。发现各向同性模型和土体CT切片的两点相关系数的最大差值是各向异性模型的2.5倍以上，且流体粒子在孔道中线处渗流速度最大0.03，在孔壁处渗流速度接近于0。表明各向异性模型在孔隙结构和空间分布上更接近真实土体情况；流体在渗流后期出现大孔隙优先流现象和指进现象。

Numerical Simulation of Anisotropically Reconstructed Soil Seepage Based on Improved QSGS-LBM Method

In order to solve the defect that the reconstruction model generated by the traditional Quartet Structure Generation Set (QSGS) method is inconsistent with the real soil condition, an improved QSGS method is proposed. The anisotropic porous medium model was constructed by considering anisotropy, porosity and soil particle size, and then the validity of the reconstruction model was verified by the two-point autocorrelation function method. Combined with the lattice Boltzmann method(LBM), the Matlab automatic program was used to simulate the meso-seepage process. It finds that the maximum difference of the two-point correlation coefficient between the isotropic model and the CT scan slice is more than 2.5 times that of the anisotropic model. The seepage velocity of fluid particles at the middle line of the channel is the largest, and the seepage velocity at the wall of the channel is close to 0. It shows that the anisotropic model is closer to the real soil in pore structure and spatial distribution. Pore preferential flow and finger-in phenomena exists in the later stage of seepage.