Abstract: To accurately analyze proppant transport in rough intersecting fractures and elucidate the interaction mechanisms among liquid, particles, and rough walls, this study reconstructed a numerical model of fractures in inhomogeneous reservoirs with varying brittleness index (BI). Various auto-correlation Gaussian rough fracture models were created using Matlab to assess roughness through the fractal dimension method. This research innovatively combined Boolean operations to establish three-dimensional rough fracture models, incorporating (Computational Fluid Dynamics) CFD-DEM (Discrete Element Method) with a bidirectional method for cosimulation. The proppant transport in fractures was categorized into three zones based on the difference in the turbulent kinetic energy. Artificially induced fracture roughness increases fluid retention and turbulence, causing plugging effects and limiting proppant flow into branch fractures. Additionally, compared with the superior deposition and significant support effects of the spherical proppant, the low-sphericity proppant traveled farther under fracturing fluid, inducing more pronounced plugging near curved fracture intersections; the variation in fracture intersection angles primarily impacted the wall shear stress within the flow field, indicating smaller angles led to higher shear energy at the intersection. Compared with the intersection angle of 30°, the height and area deposited in the 90° branch fracture increased by 52.25% and 65.33%, respectively; notably, injecting proppant from smaller to larger particles (S:M:L) and a low velocity effectively ensured fracture conductivity near the wellbore at joint roughness coefficient (JRC) ≥46 while achieving satisfactory placement in the branch fracture, making it a recommended approach.