本发明提供了一种煤矿井下水力造穴工艺模拟优化方法与系统，获取目标煤层参数，建立三维煤层结构，进行模型颗粒填充，形成煤层数值模型；设置造穴参数，建立造穴系统的模型三维结构，按配比生成磨料颗粒和液体颗粒，设定边界条件；对煤、水和磨料的颗粒间进行参数标定并平衡模型；进行离散单元法和计算流体动力学耦合的水力造穴过程数值模拟计算，模拟水力造穴作业过程，记录模拟过程中目标造穴区域内局部应力、密度和孔隙率的分布情况；根据不同造穴参数下的多次数值模拟计算结果，分析卸压增透效果对各水力造穴参数的响应关系，基于分析成果，确定水力造穴参数的优化方案，得到最终水力造穴工艺方案。本发明提高了模拟的准确性和适用性。