In addressing the challenges of high energy consumption and low efficiency in fertilization borehole drilling for clayey soils in southern orchards, this study utilizes the Discrete Element Method to establish a simulation model for clayey soils. Through this approach, we identify an optimal set of operational parameters that significantly reduces energy consumption. By utilizing simulation technology to model the drilling process, we analyzed the impact of rotation speed and feed rate on the torque and resistance of the drilling apparatus. Initially, this paper describes field measurements of particle parameters in soils from southern orchards. Subsequently, utilizing the Discrete Element Method and particle contact theory, we established a simulation model to represent the interactions between soil and soil, as well as soil and auger in the soil environment of the southern region. For the Southern orchard clay with a moisture content of 16.8% and a measured angle of repose of 35.55°, parameter calibration was performed. The contact model “Hertz-Mindlin with Johnson-Kendall-Roberts” was selected in EDEM. Using Design Expert, a regression model variance analysis was applied to the discrete element model parameters, leading to the determination of the optimal values for significant soil model parameters. The soil JKR surface energy is 5.85 J·m−2, with a soil–soil restitution coefficient of 0.65 and a soil–steel static friction coefficient of 0.5. Subsequently, discrete element simulation experiments on the drilling apparatus were conducted in EDEM, considering various rotation speeds and feed rates. The simulation analysis indicates that the torque consistently increases with higher rotation speeds, with a maximum relative error of 7%. The torque initially rises from zero to a maximum value, then gradually decreases to a low value, followed by a rapid increase to a higher value, and finally drops back down. This cycle repeats in the observed pattern. The total force experienced reaches its minimum average value of 200 N at a feed rate of 0.05 m/s. Simulation test results indicate that, among the three forces acting on the auger (vertical resistance, horizontal resistance, and lateral resistance), vertical resistance is the primary factor contributing to power consumption. As the rotation speed increases, the maximum value of vertical resistance continues to rise, while horizontal resistance and lateral resistance exhibit a declining trend. As the feed rate increases, the maximum values of resistance in all three directions also increase. When the feed rate exceeds 0.05 m/s, the maximum lateral resistance experiences a sharp increase. Through comprehensive analysis, the optimal operational parameters for borehole fertilization are determined to be a rotation speed of 100 r·min−1 and a feed rate of 0.05 m/s. The aim of this study is to reduce the energy consumption of borehole fertilization operations, minimize carbon emissions, and promote the sustainable development of orchard production.