所述两期时移地震数据中炮点与预设点之间的误差距离。
[0123]在一些实施例中，可以分别计算所述两期时移地震数据中炮点与预设点之间的误差距离。具体的，可以包括:对两期施工重叠范围进行预设网格剖分，将所述两期施工重叠范围划分为预设数量的预设网格区间；计算所述预设网格区间内炮点与所述预设网格区间中的预设点的距离。
[0124]S202:统计所述误差距离的数值分布，得到所述两期时移地震数据中炮点位置的误差分布数据，将所述炮点位置的误差分布数据作为质量监测数据。
[0125]在一些实施例中，在步骤S201获得所述炮点与预设点之间的距离之后，可以统计所述误差距离的数值分布，得到所述两期时移地震数据中炮点位置的误差分布数据，将所述炮点位置的误差分布数据作为质量监测数据。具体的，可以统计所述第一期的误差距离的数值分布，得到所述第一期时移地震数据中炮点位置的误差分布数据，将所述炮点位置的误差分布数据作为质量监测数据；可以统计所述第二期的误差距离的数值分布，得到所述第二期时移地震数据中炮点位置的误差分布数据，将所述炮点位置的误差分布数据作为质量监测数据。
[0126]如图3所示的是本申请实施例提供的一种所述两期时移地震数据中炮点位置的误差分布数据图。图中包括所述第一期时移地震数据中炮点位置的误差分布数据和所述第二期时移地震数据中炮点位置的误差分布数据。
[0127]由此可见，通过本申请第一方式获得的质量监测数据可以很好的反映炮点位置偏差情况，从而可以利用所述质量监测数据对采集的时移地震数据进行有效的质量监测。
[0128]图4是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图，结合附图4，第二方式包括:
[0129]S203:分别计算两期时移地震数据的羽角值。
[0130]在一些实施例中，可以分别计算两期时移地震数据的羽角值。图5是本申请实施例提供的分别计算两期时移地震数据的羽角值过程的流程图；具体的，结合附图5，可以包括:
[0131]S510:获取所述两期时移地震数据的测线上炮点所对应的检波点羽角值。
[0132]在一些实施例中，可以依次获取所述两期时移地震数据的每条测线上每个炮点所对应的检波点羽角值。所述检波点的羽角值可以包括航向矢量与由所述检波点坐标至其所在拖缆测线上第一个检波点坐标确定的检波矢量之间的夹角。
[0133]进一步的，所述航向矢量可以包括由采集观测系统设计时确定的当前拖缆测线上拖缆船的行进方向矢量。
[0134]进一步的，所述拖缆测线上第一个检波点可以包括所述拖缆上按距离拖船由近至远排列的第一个检波点。
[0135]进一步的，所述航向矢量与所述检波矢量间夹角范围为一 180°至+180°。具体的，当由航向矢量在180°范围内按顺时针方向可旋转至所述检波矢量时，可以将所述航向矢量与所述检波矢量间夹角记作正；当由所述航向矢量在180°范围内按逆时针方向可旋转至所述检波矢量时，可以将所述航向矢量与所述检波矢量间夹角记作负。
[0136]S520:计算所述炮点所对应的检波点羽角值的均值，将所述检波点羽角值的均值作为所述炮点的羽角值。
[0137]在实际应用中，一般一个炮点对应多个检波点，可以计算所述炮点所对应的检波点羽角值的均值，将所述检波点羽角值的均值作为所述炮点的羽角值。
[0138]S530:计算所述两期时移地震数据中测线上炮点的羽角值的均值，将所述炮点的羽角值的均值作为所述两期时移地震数据的羽角值。
[0139]在实际应用中，一般一条测线上对应多个炮点，可以计算所述两期时移地震数据中测线上炮点的羽角值的均值，将所述炮点的羽角值的均值作为所述两期时移地震数据的羽角值。
[0140]S204:统计所述测线羽角值的数值分布，得到所述两期时移地震数据的羽角值的数值分布数据，将所述羽角值的数值分布数据作为质量监测数据。
[0141]在一些实施例中，在步骤S203之后，可以统计所述测线羽角值的数值分布，得到所述两期时移地震数据的羽角值的数值分布数据，将所述羽角值的数值分布数据作为质量监测数据。具体的，可以统计所述第一期测线羽角值的数值分布，得到所述第一期时移地震数据的羽角值的数值分布数据；可以统计所述第二期测线羽角值的数值分布，得到所述第二期时移地震数据的羽角值的数值分布数据。
[0142]如图6所示的是本申请实施例提供的一种两期时移地震数据的羽角分布图；图中包括第一期时移地震数据的羽角分布和第二期时移地震数据的羽角分布。由图中可见，第一期采集过程中羽角值主要分布于±10°以内，第二期期采集过程中羽角值主要分布于±5°以内。
[0143]由此可见，通过本申请第二方式获得的质量监测数据可以很好的反映检波点位置偏差情况，从而可以利用所述质量监测数据对采集的时移地震数据进行有效的质量监测。
[0144]图7是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图，结合附图7，第三方式包括:
[0145]S205:获取所述第一期时移地震数据中第一炮点的羽角值。
[0146]在一些实施例中，可以获取所述第一期时移地震数据中第一炮点的羽角值。具体的，所述第一炮点可以包括所述第一期时移地震数据中的任一炮点。
[0147]S206:确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点，获取所述第二炮点的羽角值。
[0148]在一些实施例中，可以确定所述第二期时移地震数据中与所述第一炮点位置相对应的第二炮点，获取所述第二炮点的羽角值。具体的，可以将所述第二期时移地震数据中与所述第一炮点距离最近的炮点作为第二炮点。
[0149]S207:计算所述第一炮点的羽角值与所述第二炮点的羽角值的差值，得到所述两期时移地震数据的羽角差数据。
[0150]在一些实施例中，在步骤S206之后，可以计算所述第一炮点的羽角值与所述第二炮点的羽角值的差值，得到所述两期时移地震数据的羽角差数据。
[0151]S208:统计所述羽角差数据的数值分布，得到所述两期时移地震数据的羽角差分布数据，将所述羽角差分布数据作为质量监测数据。
[0152]在一些实施例中，在步骤S207之后，可以根据所述第一期时移地震数据的炮点羽角差确定所述两期时移地震数据的羽角差分布，所述羽角差分布作为质量监测数据。
[0153]由此可见，通过本申请第三方式获得的质量监测数据可以很好的反映两期检波点位置之间的偏差情况，从而可以利用所述质量监测数据对采集的时移地震数据进行有效的质量监测。
[0154]图8是本申请一种海上拖缆采集的时移地震数据的质量监测方法的实施例的一种流程图，结合附图8，第四方式包括:
[0155]S209:根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围。
[0156]在一些实施例中，可以根据所述第一期时移地震数据中的第一炮点的位置确定出预设处理范围。具体的，可以包括:从所述第一期时移地震数据中获取第一炮点的位置；将以所述第一炮点的位置为圆心，且以预先设置的炮点误差阈值为半径的圆作为预设处理范围。
[0157]进一步的，所述第一炮点可以包括所述第一期时移地震数据中的任一炮点。
[0158]S210:获取所述预设处理范围内的第二期时移地震数据。
[0159]在一些实施例中，在步骤S209之后，可以获取所述预设处理范围内的第二期时移地震数据。
[0160]S211:确定所述第一炮点的第一炮检对，从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对。
[0161]在一些实施例中，在步骤S210之后，可以确定所述第一炮点的第一炮检对，从所述预设处理范围内的第二期时移地震数据中获取与所述第一炮检对相对应的第二炮检对。具体的，可以包括:判断所述预设处理范围内的第二期时移地震数据中的炮检对的信息是否与所述第一炮检对的信息相匹配；当判断结果为不匹配时，从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对。当判断结果为匹配时，将所述预设处理范围内的第二期时移地震数据中与所述第一炮检对的信息相匹配的炮检对的信息作为第二炮检对。具体的，所述第一炮检对可以包括所述第一炮点中的任一炮检对。
[0162]进一步的，所述从所述预设处理范围内的第二期时移地震数据中选取符合预设条件的炮检对的信息作为第二炮检对可以包括:计算所述第一炮检对的炮点与所述预设处理范围内的第二期时移地震数据中第三炮检对的炮点间的炮点距离，计算所述第一炮检对的检波点与所述第三炮检对的检波点间的检波点距离，计算所述炮点距离与所述检波点距离之和；当所述炮点距离小于所述预先设置的炮点误差阈值，且所述检波点距离小于预先设置的检波点误差阈值，且所述炮点距离与所述检波点距离之和小于预先设置的炮检对误差阈值时，将数值最小的所述炮点距离与所述检波点距离之和所对应的炮检对信息作为第二炮检对。
[0163]进一步的，所述预先设置的炮点误差阈值可以包括两期时移地震数据中炮点的最大位置间距误差值，且当两个炮点间的距离小于所述炮点误差阈值时，可以判断两个炮点的位置为重合。具体的，所述炮点误差阈值可以根据应用场景进行预先设定。在一个具体的实施例中，所述炮点误差阈值可以设置为70m，但本申请实施例并不以此为限。
[0164]进一步的，所述预先设置的检波点误差阈值可以包括所述两期时移地震数据中检波点的最大位置间距误差值，且当两个检波点间的距离小于所述炮点误差阈值时，可以判断两个检波点的位置为重合。具体的，所述检波点误差阈值可以根据应用场景进行预先设定。在一个具体的实施例中，所述检波点误差阈值可以设置为50m，但本申请实施例并不以此为限。
[0165]进一步的，所述预先设置的炮检对误差阈值可以包括两期时移地震数据中炮