地震波形反演因其具有广泛的应用背景和潜在的巨大经济效益,近年来一直是理论和应用研究的热点和焦点。开展实用、可靠的波形反演方法研究具有十分重要的理论意义和现实意义。本文首先介绍了地震波形反演的背景及研究现状,针对现有地震模型不能很好地近似实际地层物性参数的问题,讨论粘弹性波动方程在模拟实际情况时的优势。在介绍了粘弹性介质的基本理论后,给出Kelvin粘弹性介质波动方程及其离散形式,并进行了正演模拟。在模型参数为均匀及非均匀两种不同的情况下,给出不同时刻的波场快照,从而验证了正演模拟的准确性。而后,在粘弹性波动方程正演模拟的基础上,将反演问题转化为最小二乘泛函极值问题。针对反演问题固有的不适定性和方程的非线性性,引入全变分(Total Variation或TV)正则化,使方法能更好地反演出物性参数中阶梯型间断的部分。针对粘弹性波动方程的特殊性,构造了相应的最小二乘目标泛函,进一步在离散的情形下给出了TV罚项的矩阵形式。为了求解TV正则化泛函的极值问题,将目标泛函逐步线性化,进而构造出相应的Gauss-Newton法,并给出了具体的迭代格式以及离散形式下的迭代格式。由于多重网格方法能有效地避免求解过程陷入局部极值,所以将多重网格法应用于上述非线性问题。本文中选取了二重网格迭代格式、九点延拓算子和限制算子,构造了求解粘弹性波形反演的多重网格全变分正则化方法。在数值模拟中,分别针对分层模型和孔隙模型进行了反演,给出了多重网格求解过程中参数数值的变化过程以及最终的反演结果。从得到的结果中可以看到方法很好地反演出了参数中的特殊部分,进而使得方法在实际应用中更能反映出地下介质的构造。