1.一种各向异性地震成像方法，其特征在于，包括以下步骤：步骤1：读入各向异性参数模型、P波速度模型、参数文件；

步骤2：使用Runge‑Kutta法对炮点沿不同方向进行各向异性射线追踪，并计算每一条射线对应的射线束信息；

步骤3：将单炮地震记录划分为多个以窗为单元的数据体；

步骤4：计算各个窗中数据体对时间及空间的偏导数，并对窗中地震记录进行局部平面波分解；

步骤5：对窗中心沿不同方向进行各向异性射线追踪，并计算每一条射线对应的射线束信息；

步骤6：对炮点和窗中心所有射线束对按照加了权系数的成像公式进行成像计算；

所述加了权系数的成像公式为：

其中，Is表示单炮成像值；x表示成像点的位置；ps、pr分别为炮点和窗中心点发出射线的慢度参数；A表示振幅；Ds为局部平面波分解结果；L为窗中心的位置；p′与τ′为用以局部倾斜叠加的慢度和走时参数；成像公式中权重系数Wp′的表达式为：其中，ψ表示地震记录，ψx、ψt分别为地震记录对空间和时间的偏导数，W为满足慢度和走时要求点的集合，tj表示地震波走时；

步骤7：将所有射线束对成像结果进行叠加得到最终偏移成像结果。

2.根据权利要求1所述的一种各向异性地震成像方法，其特征在于，所述各向异性参数模型包含各向异性参数σ模型和各向异性参数ε模型；所述参数文件包含网格大小、初始波束宽度、地震道数及道间距、每道采样点数、最小和最大频率。

3.根据权利要求1所述的一种各向异性地震成像方法，其特征在于，所述步骤2中射线追踪方程组为：

其中，xi表示离散点的空间位置；pi、pn、pl表示慢度分量；τ表示地震波走时；anjkl通过公式anjkl＝cijkl/ρ进行计算，cijkl为弹性模量，ρ为密度；gj、gk为特征向量分量；

知道中心射线信息后，通过波束宽度计算公式获得相应射线束的信息，其中波束宽度w计算公式为：

其中，Vs为炮点处的速度值，σ为速度对射线路径的积分，Δθ为射线的角度间隔。

4.根据权利要求1所述的一种各向异性地震成像方法，其特征在于，所述步骤3中，窗中心间隔选择范围为200m到500m，窗的持续长度为初始波束宽度的1.5倍。