偏移记录；
[0084] (5-12)偏移记录和实际的叠后地震数据对比，调节地质模型重复（5-2)-(5-10)， 直到达到期望的输出为止；结束叠后正演，并保存最后的地质模型参数和子波。
[0085] (6)对地质模型进行叠前交错网格正演。
[0086] (6-1)根据稳定条件、频散条件结合原始地震资料自动设置网格大小、采样时间。
[0087] (6-2)根据实际地震资料定义二维模型的大小，定义参数并初始化；
[0088] (6-3)读取（5-12)保存地质模型参数和子波，并计算边界区域纵波速度、横波速 度、密度，及衰减系数；
[0089] (6-4)设置炮点的起始位置、根据跑间距，设置炮点循环变量；
[0090] (6-5)设置应力变量速度变量的初始值为0 ;
[0091] (6-6)根据设置的炮点位置，对相应网格点处对应力或速度给相应时间点的子波 值；
[0092] (6-7)根据设置的炮点位置，对相应网格点处对应力或速度给相应时间点的子波 值；
[0093] (6-8)计算内部区域正应力分量；
[0094] (6-9)计算边界区域正应力分量；
[0095] (6-10)计算内部区域速度分量；
[0096] (6-11)计算边界区域速度分量；
[0097] (6-12)更新应力分量和速度分量，进行下一次时间循环重复（6-6)-(6-12)直到 最大时间为止，输出相应炮点位置的叠前正演记录；
[0098] (6-13)进行下一炮循环重复（6-5)-(6-12)进行下一炮叠前正演，直到最大炮点 位置为止；
[0099] (6-14)对所有炮点正演记录进行常规处理、偏移并叠加形成最终剖面。
[0100] 以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精 神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种基于实际地震资料的交错网格波动方程正演的方法，其特征在于，所述基于实 际地震资料的交错网格波动方程正演的方法采用基于实际地震、测井、地质资料，利用实际 地震资料的子波，先叠后逐步迭代地质模型再叠前正演的方法对实际地质模型进行正演模 拟。2. 如权利要求1所述的基于实际地震资料的交错网格波动方程正演的方法，其特征在 于，所述基于实际地震资料的交错网格波动方程正演的方法具体包括以下步骤： 步骤一，综合研究区已有地震资料获取时间采样率、空间采样间隔、覆盖次数为后面的 程序参数设置提供参考，程序参数通过读取地震数据获得； 步骤二，根据实际地震资料提取子波并保存，有井的地方联合地震和测井提取子波，没 井的地方用地震数据做相关来提取子波 步骤三，对工区已有声波和密度测井数据进行统计分析，剔除异常数据得到工区相关 地质模型的纵波速度、横波速度、密度参数，对于没有横波测井的地方根据测井数据进行横 波反演得到； 步骤四，利用已得到的纵波速度、横波速度、密度和地质模型轮廓进行地质建模并保 存； 步骤五，读取提取的子波，对地质模型进行叠后交错网格正演； 步骤六，对地质模型进行叠前交错网格正演。3. 如权利要求2所述的基于实际地震资料的交错网格波动方程正演的方法，其特征在 于，所述步骤五具体包括： 第一步，根据稳定条件下式、频散条件结合原始地震资料设置网格空间步长、时间步 长，频散条件最短波长达到6个网格以上：其中，At，Δχ，Δ z为时间步长和X，z方向的步长，Vp，Vs为纵、横波速度； 第二步，根据实际地震资料定义二维模型的大小，定义参数并初始化； 第三步，读取地质模型参数和给定子波，计算边界区域纵波速度、横波速度、密度，及衰 减系数；边界区域的纵波速度、横波速度、密度与离它最近的内部区域网格纵波速度、横波 速度、密度相等，内部区域的衰减系数为〇,边界区域的衰减系数根据下式计算：其中Vniax为最大纵波速度，PML为匹配层宽度，X为横向距离，z为纵向距离，R为理想 的边界反射系数，R取为〇. 〇〇〇〇〇l，d(x)、d(z)不等于零时表示衰减，当d(x)、d(z)等于零 时表示不衰减； 第四步，计算模型的反射系数，反射系数的计算公式根据下式计算： R= (r2*v2-rl*vl) / (r2*v2+rl*vl)； 其中R为反射系数，vl为第一层的速度，rl为第一层的密度，v2为第二层的速度，r2 为第二层的密度； 第五步，利用爆炸反射面原理（牟永光，2007)，设置初始条件； 第六步，根据下式计算内部区域应力分量；其中弹性参数λ、μ根据纵波速度、横波速 度、密度算出：第七步，根据下面三个公式计算边界区域应力分量；其中弹性 参数λ、μ根据纵波速度、横波速度、密度算出，d(x)、d(z)由公式算出： 第八步，根据下式计算内部区域速度分量；第九步，根据下式计算边界区域速度分量；第十步，更新应力分量和速度分量，进行下一次时间循环重复第五步-第十步直到最 大时间为止，输出叠后正演记录； 第十一步，对输出的叠后正演记录进行叠后逆时偏移，输出偏移记录； 第十二步，偏移记录和实际的叠后地震数据对比，调节地质模型重复第二步-第十步， 直到达到期望的输出为止；结束叠后正演，并保存最后的地质模型参数和子波。4.如权利要求2所述的基于实际地震资料的交错网格波动方程正演的方法，其特征在 于，所述步骤六具体包括： 第一步，根据稳定条件下式、频散条件结合原始地震资料设置网格空间步长、时间步 长，频散条件最短波长达到6个网格以上：其中，At，Δχ，Δ z为时间步长和X，z方向的步长，Vp，Vs为纵、横波速度； 第二步，根据实际地震资料定义二维模型的大小，定义参数并初始化； 第三步，读取地质模型参数和给定子波，设置的炮点位置，计算边界区域纵波速度、横 波速度、密度，及衰减系数；边界区域的纵波速度、横波速度、密度与离它最近的内部区域网 格纵波速度、横波速度、密度相等，内部区域的衰减系数为〇,边界区域的衰减系数根据下式 计算：其中Vniax为最大纵波速度，PML为匹配层宽度，X为横向距离，z为纵向距离，R为理想 的边界反射系数，R取为〇. 〇〇〇〇〇l，d(x)、d(z)不等于零时表示衰减，当d(x)、d(z)等于零 时表示不衰减； 第四步，设置应力变量、速度变量的初始值为〇 ; 第五步，根据设置的炮点位置，对相应网格点处对应力或速度给相应时间点的子波 值； 第六步，根据下式计算内部区域应力分量；其中弹性参数λ、μ根据纵波速度、横波速 度、密度算出：第七步，根据下面三个公式计算边界区域应力分量；其中弹性 参数λ、μ根据纵波速度、横波速度、密度算出，d(x)、d(z)由公式算出： CN 105158797 A权利要求j_ _4/4 页 第八步，根据下式计算内部区域速度分量；第九步，根据下式计算边界区域速度分量；第十步，更新应力分量和速度分量，进行下一次时间循环重复第五步-第十步直到最 大时间为止，输出相应炮点位置的叠前正演记录； 第十一步，进行下一炮循环重复第五步-第十步进行下一炮叠前正演，直到最大炮点 位置为止； 第十二步，对所有炮点正演记录进行常规处理、偏移并叠加形成最终剖面。5. -种使用权利要求1-4任意一项所述基于实际地震资料的交错网格波动方程正演 的方法的油气地震勘探正演数值模拟系统。6. -种使用权利要求1-4任意一项所述基于实际地震资料的交错网格波动方程正演 的方法的野外地震数据采集系统。
【专利摘要】本发明公开了一种基于实际地震资料的交错网格波动方程正演的方法，采用实际地震数据的子波进行正演模拟，克服现有技术对于地震资料匹配性较差的缺陷；应用高阶差分格式以及利用多源信息对正演模型参数进行约束，提高正演精度，克服频散现象；紧密结合已知多源资料(地震、测井、地质等资料)，优化正演技术流程与地质建模流程，提前设置好观测系统，在正演过程中就不需要手动去修改观测系统参数，提高计算效率与计算精度。本发明采用基于实际地震、测井、地质资料，利用实际地震资料的子波，先叠后逐步迭代地质模型再叠前正演的方法对实际地质模型进行正演模拟，可以灵活地匹配各个地区的地震资料，具有较高的计算精度和计算效率。
【IPC分类】G01V1/28
【公开号】CN105158797
【申请号】CN201510677102
【发明人】钟雨, 阎建国, 李远娟, 吴杰, 刘海波, 闫小伟, 徐松
【申请人】成都理工大学
【公开日】2015年12月16日
【申请日】2015年10月16日