叠后地震剖面的线性化波阻抗反演

本文提出了由井旁道出发逐道进行的叠后地震剖面的线性化波阻抗反演方法。每一道反演采用逐次线性化对波阻抗进行修正，反演的纵向分辨率为一个采样时间间隔，并以合成记录道与实际地震道的相似系数作为迭代收敛准则，已处理的地震道波阻抗作为下一道反演初始值。实际资料处理表明，该方法具有较高的分辨率和运算速度快的优点。     

有限差分波动方程正演模拟震源处理

为克服有限差分波动方程正演模拟中处理震源的复杂性,有两种方法可供选择:一种是从波动方程出发,利用格林公式求得波动方程的有限差分格式,从而使震源项得到简单处理;另一种是利用δ函数性质,求取δ函数在离散情况下的近似式,表示震源项简单的离散表达式。这两种方法的结果是一样的。文中最后以数值模拟例子检验了方法的有效性。     

三角网格中地震波初至旅行时的计算

本文提出一种在三角网格中计算地震波初至旅行时的精确、快速新算法，在正三角网格中，我们给求解程函方程的有限差分格式，并采用扫描算法计算出模型中离散网格点上的精确初至时间，该方法易于实现并行计算，算法中考虑了透射波，散射波及首波。因此，对于速度任意变化的模型都能得到较好的计算结果。     

叠后地震剖面的线性化波阻抗演

本文提出了由井旁道出发逐道进行的叠后地震剖面的线性化波阻抗反演方法。每一道反演采用逐次性化对波阻抗进行修正，反演的纵向分辨率为一个采样时间间隔，并拟合成记录道与实际地震道的相似系数作为迭代收敛准则，已处理的地震道波阻抗作为下一道反演初始值。     

各向异性介质中褶积微分算子法三分量地震资料的数字仿真

本文通过利用逆傅氏变换技术，得出一种设计常规微分算子的简易方法，对其截断效应引入汉宁窗加以压制；并将这种褶积微分算子用于各向异性介质中三分量弹性地震波场的数字模拟。     

快速模拟退火算法及应用

本文以广义Ｂｏｌｔｚｍａｎｎ－Ｇｉｂｂｓ统计理论为基础，采用依赖于温度的似Ｃａｕｃｈｙ分布产生新的扰动模型，建立一种新的快速模拟退火算法。文中给出了这种条件下的接收概率和降温方式的具体计算公式。新算法可以高温情况下进行大范围的搜索，在低温时仅在当前模型附近进行搜索，而且由于似Ｃａｕｃｈｙ分布有一平坦的“尾巴”，使其易于跳出局部极值，从而加快了这种模拟退火算法的收敛速度。     

有限差分法射线追踪

以往使用的求取地震波初至时间的射线追踪法具有不同程度的缺陷，因而在实际应用中存在一定的局限性。本文提出一种将局部算法与扫描算法相结合的、用有限差分解程函方程求取地震波初至时间的方法。方法中考虑了射波、首波及透射波的存在，因而可以适用于较复杂的速度模型。在求出初至时间后，利用互易性原理，能准确地得到射线路径。理论模型的试算表明，该方法计算精度较高，适宜地震偏移，层析成像以及Ｇｒｅｅｎ函数的计算。     

地震反射波波前法射线追踪

利用有限差分法解程函方程，分别求出震源点及接收到反射界面任意点的初致旅行时，将两者相加，即可得到总旅行时与界面水平坐标关系曲线。由Ｆｅｒｍａｔ原理要，该曲线上的极值点相对应的时间为反射旅行量，相对应的水平坐标值为反射点位置。本方法不但可计算出任意检皮点能接收到的所有的反射射波的旅行时及射线路径，而且具有快速，准确，适应复杂构造模型的优点。     

利用哈特莱变换进行井间声波波场正演模拟

在二维跨井声波波动方程正演模拟中,采用哈特莱变換(Hartley Transform)求取空间导数,具有占用内存比FFT少.求奇次导数时不产生尼奎斯特误差等优点,且计算速度较FFT快约一倍。文中对跨井声波波场理论模型模拟的结果表明.震源激发的位置对波场能量的分布影响很大;震源位于低速层中激发时,低速层内有较强的能量.且波场较为复杂;震源位于高速层激发时,能量很快地向周围介质扩散。