[0018] 步骤三，对工区已有声波、密度测井资料进行处理、整理得到工区相关地质模型的 纵波速度、横波速度、密度参数；
[0019] 步骤四，利用已得到的纵波速度、横波速度、密度和地质模型轮廓进行地质建模并 保存；
[0020] 步骤五，对地质模型进行叠后交错网格正演；
[0021] 步骤六，对地质模型进行叠前交错网格正演。
[0022] 进一步，所述步骤五具体包括：
[0023] 第一步，根据稳定条件、频散条件结合原始地震资料设置网格大小、采样时间；
[0024] 第二步，根据实际地震资料定义二维模型的大小，定义参数并初始化；
[0025] 第三步，读取纵波速度、密度、子波，并计算边界区域纵波速度、密度及衰减系数；
[0026] 第四步，计算模型的反射系数；
[0027] 第五步，利用爆炸反射面原理，设置初始条件；
[0028] 第六步，计算内部区域正应力分量；
[0029] 第七步，计算边界区域正应力分量；
[0030] 第八步，计算内部区域速度分量；
[0031] 第九步，计算边界区域速度分量；
[0032] 第十步，更新应力分量和速度分量，进行下一次时间循环重复第五步-第十步直 到最大时间为止，输出叠后正演记录；
[0033] 第十一步，对输出的叠后正演记录进行叠后逆时偏移，输出偏移记录；
[0034] 第十二步，偏移记录和实际的叠后地震数据对比，调节地质模型重复第二步-第 十步，直到达到期望的输出为止；结束叠后正演，并保存最后的地质模型参数和子波。
[0035] 进一步，所述步骤六具体包括：
[0036] 第一步，根据稳定条件、频散条件结合原始地震资料自动设置网格大小、采样时 间；
[0037] 第二步，根据实际地震资料定义二维模型的大小，定义参数并初始化；
[0038] 第三步，读取保存地质模型参数和子波，并计算边界区域纵波速度、横波速度、密 度，及衰减系数；
[0039] 第四步，设置炮点的起始位置、根据跑间距，设置炮点循环变量；
[0040] 第五步，设置应力变量速度变量的初始值为0 ;
[0041] 第六步，根据设置的炮点位置，对相应网格点处对应力或速度给相应时间点的子 波值；
[0042] 第七步，根据设置的炮点位置，对相应网格点处对应力或速度给相应时间点的子 波值；
[0043] 第八步，计算内部区域正应力分量；
[0044] 第九步，计算边界区域正应力分量；
[0045] 第十步，计算内部区域速度分量；
[0046] 第^^一步，计算边界区域速度分量；
[0047] 第十二步，更新应力分量和速度分量，进行下一次时间循环重复第六步_第十二 步直到最大时间为止，输出相应炮点位置的叠前正演记录；
[0048] 第十三步，进行下一炮循环重复第五步-第十二步进行下一炮叠前正演，直到最 大炮点位置为止；
[0049] 第十四步，对所有炮点正演记录进行常规处理、偏移并叠加形成最终剖面。
[0050] 本发明的另一目的在于提供一种使用所述基于实际地震资料的交错网格波动方 程正演的方法的油气地震勘探正演数值模拟系统。
[0051] 本发明的另一目的在于提供一种使用所述基于实际地震资料的交错网格波动方 程正演的方法的野外地震数据采集系统。
[0052] 本发明提供的基于实际地震资料的交错网格波动方程正演的方法，是在 MicrosoftVisualC++编译环境下用C语言实现，采用高阶交错网格有限差分法进行叠后、 叠前正演并采用基于实际地震、测井、地质等资料，应用实际地震子波而不是理论子波，提 出了一套先叠后逐步迭代地质模型再叠前正演的优化正演模型的方法对实际地质模型进 行正演模拟，其可以灵活地匹配各个地区的地震资料，具有较高的计算精度和计算效率。本 发明的核心是应用实际地震子波而不是理论子波，利用多源信息（测井、地震、地质等资 料）对正演模型参数进行约束，提高正演精度，克服频散现象，提出了先叠后正演再叠前正 演优化正演模型的方法。与现有技术相比，具有以下优势：
[0053] (1)采用实际地震数据的子波而不是理论子波进行正演模拟，克服现有技术对于 实际地震资料匹配性较差的缺陷。
[0054] (2)应用高阶差分格式以及利用多源信息（测井、地震、地质等资料）对正演模型 参数进行约束，提高正演精度，克服频散现象。
[0055] (3)紧密结合已知多源资料（地震、测井、地质等资料），提出了一套先叠后逐步迭 代地质模型再叠前正演的优化正演模型的方法，优化正演技术流程与地质建模流程，提前 设置好观测系统，在正演过程中就不需要手动去修改观测系统参数，提高计算效率与计算 精度。
【附图说明】
[0056] 图1是本发明实施例提供的基于实际地震资料的交错网格波动方程正演的方法 流程图。
【具体实施方式】
[0057] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明 进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于 限定本发明。
[0058] 本发明提供一种基于实际地震资料的交错网格波动方程有限差分正演方法，它针 对二维地质模型进行波动方程正演，可用于逆时偏移、全波形反演等方面，也可有效地指导 油气地球物理勘探的采集、处理、解释。
[0059] 下面结合附图对本发明的应用原理作详细的描述。
[0060] 如图1所示，本发明实施例的基于实际地震资料的交错网格波动方程正演的方法 包括以下步骤：
[0061] SlOl:综合研究区已有地震资料获取时间采样率、空间采样间隔、覆盖次数为后面 的程序参数设置提供参考；
[0062] S102:根据实际地震资料提取子波并保存，有井的地方可以联合地震和测井提取 子波，没井的地方可以用地震数据做相关来提取子波（杨培杰，地震子波提取方法综述，石 油地球物理勘，2008);
[0063] S103 :对工区已有声波、密度测井资料进行处理、整理得到工区相关地质模型的纵 波速度、横波速度、密度参数；
[0064] S104:利用已得到的纵波速度、横波速度、密度和地质模型轮廓进行地质建模并保 存；
[0065] S105 :对地质模型进行叠后交错网格正演（周学明，交错网格高阶差分数值模拟 及叠前逆时偏移，2010年6月，长安大学硕士论文。）；
[0066] S106 :对地质模型进行叠前交错网格正演（周学明，交错网格高阶差分数值模拟 及叠前逆时偏移，2010年6月，长安大学硕士论文。）。
[0067] 本发明的具体步骤包括：
[0068] (1)综合研究区已有地震资料获取时间采样率、空间采样间隔、覆盖次数为后面的 程序参数设置提供参考。
[0069] (2)根据实际地震资料提取子波并保存，有井的地方可以联合地震和测井提取子 波，没井的地方可以用地震数据做相关来提取子波。
[0070] (3)对工区已有井资料进行处理、整理得到工区相关地质模型的纵波速度、横波速 度、密度参数；
[0071] (4)利用已得到的纵波速度、横波速度、密度和地质模型轮廓进行地质建模并保 存。
[0072] (5)对地质模型进行叠后交错网格正演。
[0073] (5-1)根据稳定条件、频散条件结合原始地震资料设置网格大小、采样时间。
[0074] (5-2)根据实际地震资料定义二维模型的大小，定义参数并初始化；
[0075] (5-3)读取纵波速度、密度、子波，并计算边界区域纵波速度、密度及衰减系数；
[0076] (5-4)计算模型的反射系数；
[0077] (5-5)利用爆炸反射面原理，设置初始条件；
[0078] (5-6)计算内部区域正应力分量；
[0079] (5-7)计算边界区域正应力分量；
[0080] (5-8)计算内部区域速度分量；
[0081] (5-9)计算边界区域速度分量；
[0082] (5-10)更新应力分量和速度分量，进行下一次时间循环重复（5-5)-(5-10)直到 最大时间为止，输出叠后正演记录；
[0083] (5_11)对输出的叠后正演记录进行叠后逆时偏移，输出