基于数字岩心技术的礁灰岩三维精细离散元数值建模方法及其应用

Three-dimensional refined discrete element numerical modeling method and its application for reef limestone based on digital core technology

Xin Wei^a, Dengxing Qu^a,b, Zhengrong Zhou^c, Xinping Li^a,b, Yingwei Zhu^a,b, Shaohua Hu^d, Wenhao Li^a,b

Abstract: The increasing scale of reef engineering construction has heightened the importance of understanding reef limestone’s mechanical properties in load-bearing strata. High-precision CT scanning provided digital core data, enabling precise characterization of pore-matrix structures through multi-level filtering and adaptive threshold segmentation. To enhance spatial mapping accuracy and computational efficiency in high - resolution modeling, we developed a three-dimensional refined discrete element method incorporating BallTree algorithm and two-level Euclidean distance logical filtering. The model was validated through pore characteristic comparison and mechanical testing before conducting systematic uniaxial compression simulations. Analysis revealed that shallow weakly-cemented compact reef limestone exhibits vesicular and banded pores, demonstrating microscopic heterogeneity alongside macroscopic homogeneity. Under uniaxial loading, the limestone undergoes five distinct stages: pore compaction, linear elasticity, crack propagation, peak failure, and residual strength. Stress transmission occurs preferentially along the 45° direction, governing crack development. The failure process initiates with tensile cracking and evolves into a combined tensile-shear failure mode. Energy analysis indicates that elastic strain energy dominates storage, while sliding friction heat represents the primary dissipation mechanism during failure. This study integrates digital core with numerical simulation to elucidate the relationship between pore structure and failure evolution in reef limestone, offering new perspectives on deformation and failure mechanisms in porous rocks.

Keywords: Reef limestone; Digital core; Discrete element method; Pore structure; MatDEM

Fig.11 Linear elastic model in discrete element method.

Fig.12 Initial dense-packed cubic particle model.

Fig.20 Crack development and evolution of reef limestone under uniaxial compression.

Fig.22  The porosity evolution of reef limestone under uniaxial compression conditions.

Xin Wei, Dengxing Qu, Zhengrong Zhou, Xinping Li, Yingwei Zhu, Shaohua Hu, Wenhao Li,Three-dimensional refined discrete element numerical modeling method and its application for reef limestone based on digital core technology, Computers and Geotechnics, Volume 185,2025,107362,ISSN 0266-352X,https://doi.org/10.1016/j.compgeo.2025.107362.