基于格子Boltzmann法（LBM）与离散单元法（DEM）的基本理论，通过对格子单松弛模型演化方程进行修正，建立LBM–DEM耦合计算模型。引入浸没移动边界方法（IMB）处理复杂的流–固边界问题，利用编程软件MATLAB开发LBM–IMB–DEM耦合求解程序，对悬浮颗粒在含水介质中的运动过程进行模拟计算。通过开展层析柱强制渗流试验，模拟DKT算例，验证了计算模型与求解方法的正确性；从孔隙尺度分析地下水动力、颗粒形貌和尺寸效应对含水层渗透性能的影响。结果表明：当进口流速为5×10–5 m/s，向含水介质注入球形硅微粉悬浮液时，含水层各区域的渗流速度均呈现下降、回升、稳定的3个连续阶段；随着渗流速度的降低，颗粒滞留率增加了47.4%，所选3个断面平均流速恢复率降低。将悬浮颗粒由球形硅微粉替换为非球形硅微粉后，所选3个断面的平均流速最大降幅分别提高了15.0%、11.4%、2.6%，恢复率依次降低了37.5%、69.6%、71.9%，含水介质的颗粒滞留率增加了47.9%。基于颗粒受力平衡分析可知，球形硅微粉转动惯量较低，易于滚动并脱离含水介质表面；同等效体积的非球形硅微粉比表面积更大，而且多絮凝成团，迁移过程中出现沉积或被孔喉捕获的机率提高，并且沉积后难以发生再释放过程。随着渗流速度的降低与含水介质计算区域的增加，悬浮颗粒形貌的变化对于含水层水动力场演化过程的影响程度增强。
