## 1 Introduction
[13] implemented uniaxial compression experiments in the vicinity of the iron ore-marble contact zone. The composite rock body exhibited a lower strength than the rock body, and the mechanical properties were found to be correlated with the contact angle of the rock samples. The establishment of a comprehensive prediction technology system, which encompasses the principles of distance and proximity, qualitative analysis, quantitative identification, structural localization, and water-bearing tectonic identification, enabled [PERSON] et al. [14] to successfully construct curtain grouting in the water-rich zone of Sichuan karst tunnels. Using curtain grouting technology for engineering treatment, [PERSON] et al. [15] analyzed the grouting effect by monitoring the change in water influx before and after the grouting operation at the quarry's water intake site. In order to reduce the water volume to less than 85% of the initial volume, [PERSON] et al. [16] implemented curtain grouting and infill in a water-rich mine voided area. [PERSON] et al. [17] conducted an analysis of the causes of the cement surge disaster, proposed the use of curtain grouting for disaster management, and evaluated it based on the inspection hole survey and P-Q-t analysis. The study of the variation of grouting pressure, grout flow rate, and the law of grout diffusion, as well as the analysis of the factors affecting grout diffusion, was conducted by [PERSON] et al. [18] using a self-developed grout diffusion visualization platform and a three-dimensional grouting experimental system. [PERSON] et al. [19] summarized the factors that must be considered during the fissure grouting process, investigated the diffusion mechanism of the slurry under the influence of the grout, and proposed an optimization method for groutingparameters and a grout quality assessment technique. Methods for optimizing grouting quality and techniques for evaluating grouting quality were suggested. Using numerical simulation, [PERSON] et al. [20] investigated the spatial evolution characteristics of plastic damage and grout diffusion in the presence of mining disturbances. The findings indicated that the diffusion range of the grout in each spatial direction in the surrounding rock increased as the ratio of lateral pressure coefficients perpendicular and parallel to the axial direction of the roadway increased. Nevertheless, the same porous media slurry diffusion is the focus of past scholars' research on the slurry diffusion law, and there are fewer distinct media between the slurry diffusion law studies. Concurrently, the majority of the aforementioned studies pertain to rock mechanical differences in the ore-rock contact zone or the curtain grouting diffusion mechanism in a specific aspect. Relatively few scholars have conducted research on curtain grouting in the ore-rock contact zone. In this paper, we utilize the theoretical equations of [PERSON]'s law for the mutually incompatible seepage of slurry water in porous media and the Navier-Stokes equations for laminar two-phase flow, in conjunction with level set functions to conduct exhaustive research on the planar diffusion law of the slurry in the contact zone of the mine rock under the water-rich condition. The context of the actual engineering conditions of the 2# vein-piercing channel at the \(-\)380 m level of the air shaft of the Zhongjiu Iron Mine of Anhui Ma Steel is in view. Numerical simulation was employed to investigate the slurry diffusion pattern in the ore-rock contact zone under water-rich conditions, as well as the impact of single-hole and group-hole grouting curtain formation. Two optimization methods for the group-hole grouting parameters were proposed based on the simulation results. Ultimately, this research will be implemented in the field of engineering to evaluate the impact of grouting on the change in water influx in exploratory openings prior to and following curtain grouting.