基于中国地震台网中心2013 MS7.0芦山地震余震数据我们首先确定了余震空间分布范围并根据G-R关系计算了主震后半小时内的完备震级Mc=3.5, 并且得到了ML≥3.5和ML≥3.0的地震在2001年至芦山地震前的背景场地震发生率.通过Omori-Ustu经验定律和两种Dieterich模型对芦山地震余震发生率的拟合, 我们发现阶梯型Dieterich模型只能模拟p=1的情况, 从而造成了模拟曲线与观测数据的差别;前人研究表明震后滑移同样是产生余震的原因, 如果假设余震序列由主震静态剪应力Δτ和震后滑移共同作用所产生, 我们数值模拟得到的对数型Dieterich模型能够较好地推断余震发生率R随时间t增加而衰减的趋势, 能够从物理机制上解释MS7.0芦山地震余震序列衰减指数大于1这一现象.通过对数型Dieterich模型的拟合并结合Andrews的方法, 我们还得到MS7.0芦山地震Aσ约为0.155 MPa, ta约为8.4年, 这一值与前人研究结果十分接近.