需求如题，操作如下：
\begin{enumerate}
\item[1.] 对每一炮，完成以下操作：
\begin{enumerate}
\item[1.1 ] 进行地震波场正演模拟，计算每一时刻的正向传播的波场值u；
\item[1.2 ] 计算检波点位置实际观测的波场与正演模拟的波场残差：$ u_cal - u_obs $；
\item[1.3 ] 将上述\textbf{残差波场}反向传播，并得到反向传播数据：$ u^{2}(x,t:x_s) $；
\item[1.4 ] 求取反传波场与正传波场在时间上的二阶导数在时间上的内积，即求取$g(v)$，可得到单炮的近似梯度。
\end{enumerate}
\item[2. ] 叠加所有炮求取的梯度值，得到模型空间的全局梯度。
\end{enumerate}
效果：
图1：枚举嵌套效果图
