本文研究了多种SPH的实现方法，并从计算稳定性、边界条件的实现和计算速度等方面对各种方法进行了比较，最后发展了一种较适合用于弹性波数值模拟的方法。该法通过采用SPH的核函数直接乘以动力方程两边，经过两次分部积分来求解波动方程，我们称之为极弱形式的核估值法。该法计算量相对较小，对应力边界的处理十分方便。在边界粒子体积和表面积的计算上采用类似虚粒子的方法，建立了一种新的计算粒子体积和表面积的方法，来模拟粒子真实的体积和表面积，取得了较好的结果。在此基础上实现了无反射边界条件。通过把守恒光滑逼近方法(Conservative Smoothing Approach，简称CSA)应用到三维来消除张拉不稳定性和模型的零能模态所带来的振荡，也取得了比较好的效果。同时提出了计算结果的反卷积处理，它使得人们可以使用较大的粒子间距获得较理想的结果。最后给出了一个爆破震动波在岩体中传播的算例。计算结果表明，SPH方法在模拟岩体中波的传播方面是一种有效的和很有发展潜力的方法。 岩体中的裂隙分布和岩体块度分布与岩体的力学参数有密切的联系，岩体的尺度信息的提取作为我们岩体结构识别这个大课题中一个辅助手段，本文对其进行了一定的探索。通过对破裂岩体或含节理、裂隙岩体的图像(数码照片)信息进行处理，得到岩体块度分布的表征岩体多尺度特征的分维指数和特征块度。