估算地震数据道集能量值的中心位置的方法和装置的制造方法
【技术领域】
[0001] 本发明涉及地震勘探技术领域，特别涉及一种估算地震数据道集能量值的中心位 置的方法和装置。
【背景技术】
[0002] 通过地震勘探采集获取的数据一般有两种记录方式：共炮点道集和共检波点道 集。其中，共炮点道集，是指所有的地震道的地震信号都来源于同一个炮点，这些地震道按 照设计的观测系统分布于炮点周围，在正常情况下，炮点位于这些地震道（地震接收道）的 中心处；共检波点道集，只有一个接收点，这个接收点所记录的地震信号来源于周围的多个 炮点。
[0003] 以前，常用的采集地震数据的方式是有线仪器采集，这种方式需要的采集人力成 本很高，同时，施工难度也很大。目前，常用的是采集方式是无线节点采集，相对于有线仪器 采集，无线节点采集不但占用人员数量少，而且在复杂地表可以很自如灵活地摆放，减少了 跨越和攀爬，尤其在高风险地区体现更为明显，进一步的，无线节点采集由于野外暴露目标 小，在高风险地区不易丢失和也不易遭到破坏。
[0004] 根据设计的观测系统，炮点或者检波点都有其确定的位置，但是采集过程中的炮 点或者检波点，存在着不准确的可能，炮点或者检波点可能被人为移动，而其移动后的位置 没有准确地被记录下来，在海上或者湖上进行地震勘探的过程中，接收点被投放到水中，通 过二次定位所确定的位置可能不准确，因此，在地震勘探过程中，炮点或者检波点的位置的 质量监控是一个重要的环节。
[0005] 对于炮点的位置的质量控制，一般采用线性动校正的方法，即，给定线性动校正的 速度，对地震数据进行时差校正，如果炮点位置正确，则初至处于水平状态，如果炮点位置 不准确，则初至会出现偏差。在这种情况下，就需要接收点均匀分布于一条直线上或者近似 于这种情况，且这种方式仅对有线仪器记录系统才有用。
[0006] 针对如何确定任意方式采集得到的地震数据道集能量值的中心位置，目前尚未提 出有效的解决方案。
【发明内容】
[0007] 本发明实施例提供了一种估算地震数据道集能量值的中心位置的方法，以达到快 速准确确定地震数据道集能量值的中心位置的目的，该方法包括：
[0008] 采集得到地震数据道集；
[0009] 确定所述地震数据道集中每一地震道的振幅绝对值的最大值，将该最大值作为该 地震道的能量值；
[0010] 从所有地震道的能量值中确定出最大值，并记录该最大值的位置坐标；
[0011] 将该最大值的位置坐标周围预定范围内的所有地震道组成预分析道集，并记录所 述预分析道集中各地震道的坐标和能量值；
[0012] 根据所述预分析道集中各地震道的坐标和能量值，估算出所述地震数据道集的能 量值的中心位置。
[0013] 在一个实施方式中，当所述地震数据道集为共检波点道集时，所述能量值的中心 位置为接收点的位置，当所述地震数据道集为共炮点道集时，所述能量值的中心位置为激 发点的位置。
[0014] 在一个实施方式中，根据所述预分析道集中各地震道的坐标和能量值，估算出所 述地震数据道集的能量值的中心位置，包括：
[0015] 根据所述预分析道集中各地震道的坐标和能量值，采用二分法确定所述预分析道 集范围内的最大的能量衰减统计值的位置；
[0016] 将确定的最大能量衰减值的位置作为所述地震数据道集的能量值的中心位置。
[0017] 在一个实施方式中，计算所述预分析道集范围内一点的能量衰减统计值，包括：
[0018] 计算该点与所述预分析道集中各地震道的距离；
[0019] 按照距离从小到大对所述预分析道集中的各个地震道从0到η编号；
[0020] 按照以下公式计算该点的能量衰减统计值：
[0021]
[0022] 其中，s表示该点的能量衰减统计值，i表示按照距离从小到大对所述预分析道集 中的各个地震道进行排序后各个地震道的编号，η表示所述预分析道集中地震道的道数，N 1 是在地震道i+Ι到地震道η的范围内，能量值小于E1的地震道的道数，E i表示地震道i的 能量。
[0023] 在一个实施方式中，采用二分法确定所述预分析道集范围内的最大的能量衰减统 计值的位置，包括：
[0024] 给定初始位置；
[0025] 采用二分法确定所述初始位置所在的水平方向的最大的能量衰减统计值的第一 位置；
[0026] 固定所述第一位置，采用二分法确定所述第一位置所在的垂直方向的最大的能量 衰减统计值的第二位置；
[0027] 将所述第二位置作为所述预分析道集范围内的最大的能量衰减统计值的位置。
[0028] 本发明实施例还提供了一种估算地震数据道集能量值的中心位置的装置，以达到 快速准确确定地震数据道集能量值的中心位置的目的，该装置包括：
[0029] 采集模块，用于采集得到地震数据道集；
[0030] 振幅确定模块，用于确定所述地震数据道集中每一地震道的振幅绝对值的最大 值，将该最大值作为该地震道的能量值；
[0031] 最大值确定模块，用于从所有地震道的能量值中确定出最大值，并记录该最大值 的位置坐标；
[0032] 记录模块，用于将该最大值的位置坐标周围预定范围内的所有地震道组成预分析 道集，并记录所述预分析道集中各地震道的坐标和能量值；
[0033] 中心位置估算模块，用于根据所述预分析道集中各地震道的坐标和能量值，估算 出所述地震数据道集的能量值的中心位置。
[0034] 在一个实施方式中，当所述地震数据道集为共检波点道集时，所述能量值的中心 位置为接收点的位置，当所述地震数据道集为共炮点道集时，所述能量值的中心位置为激 发点的位置。
[0035] 在一个实施方式中，所述中心位置估算模块包括：
[0036] 能量衰减统计值计算单元，用于根据所述预分析道集中各地震道的坐标和能量 值，采用二分法确定所述预分析道集范围内的最大的能量衰减统计值的位置；
[0037] 中心位置确定单元，用于将确定的最大能量衰减值的位置作为所述地震数据道集 的能量值的中心位置。
[0038] 在一个实施方式中，所述能量衰减统计值计算单元包括：
[0039] 距离计算子单元，用于计算该点与所述预分析道集中各地震道的距离；
[0040] 编号子单元，用于按照距离从小到大对所述预分析道集中的各个地震道从0到η 编号；
[0041] 计算子单元，用于按照以下公式计算该点的能量衰减统计值：
[0042]
[0043] 其中，s表示该点的能量衰减统计值，i表示按照距离从小到大对所述预分析道集 中的各个地震道进行排序后各个地震道的编号，η表示所述预分析道集中地震道的道数，N 1 是在地震道i+Ι到地震道η的范围内，能量值小于E1的地震道的道数，E i表示地震道i的 能量。
[0044] 在一个实施方式中，所述能量衰减统计值计算单元包括：
[0045] 初始位置给定子单元，用于给定初始位置；
[0046] 水平方向计算子单元，用于采用二分法确定所述初始位置所在的水平方向的最大 的能量衰减统计值的第一位置；
[0047] 垂直方向计算子单元，用于固定所述第一位置，采用二分法确定所述第一位置所 在的垂直方向的最大的能量衰减统计值的第二位置；
[0048] 位置确定子单元，用于将所述第二位置作为所述预分析道集范围内的最大的能量 衰减统计值的位置。
[0049] 在本发明实施例中，将地震数据道集中每一地震道的振幅绝对值的最大值作为该 地震道的能量值，然后从所有地震道的能量值中确定出最大值，并根据该最大值的位置坐 标周围预定范围内的所有地震道的坐标和能量值，估算出该地震数据道集的能量值的中心 位置。通过上述方式不仅可以实现对有线仪器采集的数据的能量值中心位置的监控，也能 实现对无线接收节点采集的数据的能量值中心位置的监控，且估算结果更为准确。
【附图说明】
[0050] 此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不 构成对本发明的限定。在附图中：
[0051] 图1是根据本发明实施例的估算地震数据道集能量值的中心位置的方法流程图；
[0052] 图2是根据本发明实施例的共检波点道集提取能量值三维立体示意图；
[0053] 图3是根据本发明实施例的