提出了一种研究二维弹性动力学散射问题的边界元法，其中，散射体或侵入体散场以Ｈｕｙｇｅｎｓ原理积分形式表达，虚拟场源以适配边界条件来求取，该间接边界元法可以高精度计算包括多次散射在内的全波场。通过对不同空间与尺度分布裂缝所诱导的地震波散射数值计算表明：该方法具有计算多次散射波场的能力，尾波（多次散射波）频率及延续长度受控于裂缝位置及尺度；等效介质弹性参数受控于散射体不同的统计分布。