海水侵蚀环境下礁灰岩地层三维离散元数值模型构建方法

本发明为一种海水侵蚀环境下礁灰岩地层三维离散元数值模型构建方法。包括如下步骤：建立由基础球形颗粒单元体堆积而成的初始模型；根据K‑Means聚类算法遍历所有颗粒、选取K个溶蚀种子颗粒，循环寻找所有溶蚀种子距离最小的一个颗粒，并定义为聚类颗粒；定义孔隙率并设置目标孔隙率；每聚类一次计算当前孔隙率，如果达到目标孔隙率则跳出循环，反之继续聚类；当孔隙率达到目标值，聚类结束并删除所有聚类及溶蚀种子颗粒，剩余颗粒构成岩体部分；赋予数值模型接触参数，验证接触参数标定的有效性；构建多尺度礁灰岩三维离散元数值模型。本发明构建了多级地层礁灰岩三维离散元数值模型，为实际工程中礁灰岩地层硐室稳定性等力学问题提供了依据。