一种三维地震观测系统偏移噪声获取方法及装置制造方法及图纸

本发明专利技术实施例提供一种三维地震观测系统偏移噪声获取方法及装置，所述方法包括：获取三维地震观测系统数据，并获取三维网格速度模型；利用所述三维地震观测系统数据和所述三维网格速度模型，进行地震波波场延拓，获取地震波波场延拓结果；根据所述地震波波场延拓结果，分别计算炮点偏移噪声和检波点偏移噪声；利用所述炮点偏移噪声和所述检波点偏移噪声，计算所述三维地震观测系统偏移噪声。本发明专利技术实施例适用于复杂介质条件且计算精度较高。

【技术实现步骤摘要】

本专利技术涉及三维地震观测系统分析评价技术，尤其涉及一种三维地震观测系统偏移噪声获取方法及装置。
技术介绍
理想的三维地震观测系统为检波点呈网格状分布的区域观测系统，能够实现地震波场的完整采样。但在实际地震勘探中，由于受勘探时间与预算的限制，必须采用稀疏的三维地震观测系统(比如正交型、斜交型与平行型等)。稀疏的三维地震观测系统在减小勘探成本的同时，也造成了采样的不均勻，因此必须对稀疏的三维地震观测系统采样的合理性进行评价。常规的三维地震观测系统设计基于地震勘探中水平地层共中心点叠加理论，主要通过覆盖次数，炮检距分布与方位角分布等来判断三维地震观测系统设计方案的合理性。 对于倾斜地层，这种共中心点叠加理论不再有效。因此地震勘探界发展出地震共反射点分析方法，它基于宏观地震速度模型，通过地震波射线追踪手段，得到通过目的地层上面元位置的射线条数与角度。与常规共中心点分析类似，地震共反射点分析也可以得到目标面元的覆盖次数，炮检距分布与方位角分布信息。随着地震勘探叠前偏移方法逐步成为主流的地震资料处理方法，三维地震观测系统设计从地震叠加效应分析逐步转向地震叠前偏移效应分析。目前直接采用地震叠前偏移理论计算三维地震观测系统偏移噪声的方法主要有两种(1)基于Beylkin理论，结合检波点与震源点的分析方法；(2)Berkhout共聚焦分析方法。目前这两种方法均采用地震波射线追踪或地震波KirchhofT算子实现，因此仅适用于均勻介质或层状介质情况，无法适用于复杂的地下介质条件，且计算精度较低。
技术实现思路
本专利技术实施例提供一种三维地震观测系统偏移噪声获取方法及装置，以适用于复杂介质条件且计算精度较高。一方面，本专利技术实施例提供了一种三维地震观测系统偏移噪声获取方法，所述方法包括获取三维地震观测系统数据，并获取三维网格速度模型；利用所述三维地震观测系统数据和所述三维网格速度模型，进行地震波波场延拓，获取地震波波场延拓结果；根据所述地震波波场延拓结果，分别计算炮点偏移噪声和检波点偏移噪声；利用所述炮点偏移噪声和所述检波点偏移噪声，计算所述三维地震观测系统偏移噪声。另一方面，本专利技术实施例提供了一种三维地震观测系统偏移噪声获取装置，所述装置包括信息获取单元，用于获取三维地震观测系统数据，并获取三维网格速度模型；波场延拓单元，用于利用所述三维地震观测系统数据和所述三维网格速度模型，进行地震波波场延拓，获取地震波波场延拓结果；偏移噪声计算单元，用于根据所述地震波波场延拓结果，分别计算炮点偏移噪声和检波点偏移噪声，并利用所述炮点偏移噪声和所述检波点偏移噪声，计算所述三维地震观测系统偏移噪声。上述技术方案具有如下有益效果因为采用获取三维地震观测系统数据，并获取三维网格速度模型；利用所述三维地震观测系统数据和所述三维网格速度模型，进行地震波波场延拓，获取地震波波场延拓结果；根据所述地震波波场延拓结果，分别计算炮点偏移噪声和检波点偏移噪声；利用所述炮点偏移噪声和所述检波点偏移噪声，计算所述三维地震观测系统偏移噪声的技术手段，所以适用于复杂介质条件且计算精度较高。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术实施例一种三维地震观测系统偏移噪声获取方法流程图；图2为本专利技术应用实例一种三维地震观测系统偏移噪声获取方法流程图；图3为本专利技术应用实例波场延拓示意图；图4为本专利技术应用实例波场延拓的延拓层示意图；图5为本专利技术应用实例三维地震观测系统偏移噪声图；图6为本专利技术实施例一种三维地震观测系统偏移噪声获取装置结构示意图。具体实施例方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。如图1所示，为本专利技术实施例一种三维地震观测系统偏移噪声获取方法流程图， 所述方法包括101、获取三维地震观测系统数据，并获取三维网格速度模型。可选的，所述三维地震观测系统数据可以包括炮点(即地震勘探普遍采用的人工炸药震源，一般将其称为炮点)坐标、检波点(指地震勘探中在地表布设的检波器，一般至少有上千个，用于接收震源发出经地下反射后回到地表的地震波)坐标及炮点与检波点的对应关系；获取所述三维地震观测系统数据后，可以利用所述炮点与检波点的对应关系， 得到各个炮点坐标及其分别对应的检波点坐标。102、利用所述三维地震观测系统数据和所述三维网格速度模型，进行地震波波场延拓，获取地震波波场延拓结果。可选的，所述利用所述三维地震观测系统数据和所述三维网格速度模型，进行地震波波场延拓，获取地震波波场延拓结果，可以包括利用Lippmarm-Schwinger单程波动方程的离散退化算子，将各个炮点及其分别对应的检波点的地面波场分别进行向下波场延拓，得到地下目标深度的地震波波场信息，进而提取目标点水平位置及其邻近区域的波场。103、根据所述地震波波场延拓结果，分别计算炮点偏移噪声和检波点偏移噪声。可选的，所述根据所述地震波波场延拓结果，分别计算炮点偏移噪声和检波点偏移噪声，可以包括将炮点进行向下波场延拓得到的邻近区域的波场乘以目标位置波场的逆，得到所述炮点的偏移噪声矩阵；将检波点进行向下波场延拓得到的邻近区域的波场乘以目标位置波场的逆，得到所述检波点的偏移噪声矩阵。104、利用所述炮点偏移噪声和所述检波点偏移噪声，计算所述三维地震观测系统偏移噪声。可选的，所述利用所述炮点偏移噪声和所述检波点偏移噪声，计算所述三维地震观测系统偏移噪声，可以包括利用所述各个炮点的偏移噪声矩阵及所述各个炮点分别对应的检波点的偏移噪声矩阵的乘积，计算所述三维地震观测系统偏移噪声。可选的，可以根据计算获得的所述三维地震观测系统偏移噪声，绘制所述三维地震观测系统的偏移噪声图。本专利技术实施例的目的在于克服常规的基于地震波射线追踪或地震波Kirchhoff 算子的偏移噪声评估方法无法适用于复杂介质条件及其计算精度较低的缺陷，提出一种适用于任意复杂介质条件的三维地震观测系统偏移噪声快速定量评估方法。本专利技术实施例的三维地震观测系统偏移噪声分析方法，基于Lippmarm-Schwinger单程波动方程的离散退化算子表达形式，依据地下地震速度模型与三维地震观测系统设计方案，直接计算三维地震观测系统的预期偏移噪声。由于采用离散退化方法求解波动方程，该方法即使对于地下强对比度介质条件与广角成像范围也具有很高的计算精度。由于输入地震速度模型为三维网格模型，通过合理设计网格间距，该方法可适用于任意复杂的介质模型。如图2所示，为本专利技术应用实例一种三维地震观测系统偏移噪声获取方法流程图，所述方法包括201、获取三维地震观测系统数据。获取各种标准格式的三维地震观测系统数据(包含炮点坐标、检波点坐标及炮点与检波点的对应关系)，得到各个炮点坐标及其分别对应的检波点坐标。202、获取三维网格速度模型。获取三维本文档来自技高网...

【技术保护点】

【技术特征摘要】

【专利技术属性】
技术研发人员：魏伟，符力耘，
申请(专利权)人：中国科学院地质与地球物理研究所，
类型：发明
国别省市：