CN112241023A - 陆上三维炮点二次定位及定量偏移监控方法和装置 - Google Patents

陆上三维炮点二次定位及定量偏移监控方法和装置 Download PDF

Info

Publication number
CN112241023A
CN112241023A CN201910639630.6A CN201910639630A CN112241023A CN 112241023 A CN112241023 A CN 112241023A CN 201910639630 A CN201910639630 A CN 201910639630A CN 112241023 A CN112241023 A CN 112241023A
Authority
CN
China
Prior art keywords
detection point
line
minimum value
shot
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910639630.6A
Other languages
English (en)
Other versions
CN112241023B (zh
Inventor
冯发全
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201910639630.6A priority Critical patent/CN112241023B/zh
Publication of CN112241023A publication Critical patent/CN112241023A/zh
Application granted granted Critical
Publication of CN112241023B publication Critical patent/CN112241023B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种陆上三维炮点二次定位及定量偏移监控方法和装置，对距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间，进而得到所述检波点线对应的极小值处的检波点的横纵坐标；分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号，进而得到所述炮点对应的极小值处的检波点线的纵坐标和横坐标；根据所述纵坐标和横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量，能够自动实现炮点偏移的定量监控。

Description

陆上三维炮点二次定位及定量偏移监控方法和装置
技术领域
本发明涉及地球物理勘探技术领域，尤其涉及一种陆上三维炮点二次定位及定量偏移监控方法和装置。
背景技术
炮点偏移检查是地震资料采集及处理中的一项重要基础工作，对确保地震工程质量和地震成像精度至关重要。
目前，主要采用现场处理系统在屏幕上显示数据体，人工核对的方法进行炮点偏移检查，基本属于一种定量监控。这种监控方法的工作量巨大，而且有时因疏忽造成差错。特别是三维采集时，工作量更加大。例如一天施工炮数400炮，按人工每检查一炮平均用时1.5分钟计算，则检查炮偏至少需要400/(60/1.5)＝10小时，如果一天施工炮数超过800炮，则需要另增加一台处理系统和相关人员，提高地震队的生产成本，并且不能满足地震采集项目高效质量监控和定量监控的需求。
发明内容
有鉴于此，本发明提供一种陆上三维炮点二次定位及定量偏移监控方法和装置、电子设备以及计算机可读存储介质，能够自动实现炮点偏移的定量监控，有效降低了工作量，防止人工误差，降低生产成本，满足地震采集项目高效质量监控和定量监控的需求，另外，还能够对偏移炮点进行精确二次定位，供地震队复测，满足地震队对采集质量的监控要求。
为了实现上述目的，本发明采用如下技术方案：
第一方面，提供一种陆上三维炮点二次定位及定量偏移监控方法，包括：
对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间；
根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标；
分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号；
根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；
根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；
根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
进一步地，陆上三维炮点二次定位及定量偏移监控方法还包括：
根据SPS文件获取地震初至数据文件，所述地震初至数据文件包括多个炮点的地震初至数据，每个炮点的地震初至数据均包括：用于接收所述炮点产生的地震波的检波点桩号、地震初至时间、横纵坐标及其所属的检波点线。
进一步地，陆上三维炮点二次定位及定量偏移监控方法还包括：
对所述地震初至数据文件进行基准面静校正得到校正后的地震初至数据文件，以利用校正后的地震初至文件进行陆上三维炮点二次定位及定量偏移监控。
进一步地，陆上三维炮点二次定位及定量偏移监控方法还包括：
采用多项式拟合方式剔除所述地震初至数据文件中的异常数据。
进一步地，陆上三维炮点二次定位及定量偏移监控方法还包括：
分别对SPS文件中每条检波点的横纵坐标与桩号分别进行线性拟合，得到每条检波点线的检波点横纵坐标-桩号的线性转换关系。
进一步地，所述对地震初至数据文件中距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合之前，还包括：
根据预设正负偏移距分布范围对地震初至数据文件中距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行规则化。
进一步地，所述根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标之前，还包括：
根据所述极小值处的检波点附近地震道的拟合初至与实际地震初至的时差对所述极小值处的检波点的地震初至时间进行校正。
进一步地，陆上三维炮点二次定位及定量偏移监控方法还包括：
对所述炮点对应的每条检波点线的线号与对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系。
进一步地，陆上三维炮点二次定位及定量偏移监控方法还包括：
对所述炮点对应的每条检波点线对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系。
进一步地，陆上三维炮点二次定位及定量偏移监控方法还包括：
获取偏移量大于预设误差门槛的炮点编号，以根据所述炮点编号进行二次测量得到实际的横纵坐标，并利用实际的横纵坐标替换SPS文件中所述炮点的横纵坐标。
进一步地，陆上三维炮点二次定位及定量偏移监控方法还包括：
利用所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标替换SPS文件中所述炮点的横纵坐标。
第二方面，提供一种陆上三维炮点二次定位及定量偏移监控装置，包括：
第一处理模块，对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间；
第一转换模块，根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标；
第二处理模块，分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号；
第二转换模块，根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；
第三转换模块，根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；
偏移量计算模块，根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
数据获取模块，根据SPS文件获取地震初至数据文件，所述地震初至数据文件包括多个炮点的地震初至数据，每个炮点的地震初至数据均包括：用于接收所述炮点产生的地震波的检波点桩号、地震初至时间、横纵坐标及其所属的检波点线。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
静校正模块，对所述地震初至数据文件进行基准面静校正得到校正后的地震初至数据文件，以利用校正后的地震初至文件进行陆上三维炮点二次定位及定量偏移监控。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
异常数据剔除模块，采用多项式拟合方式剔除所述地震初至数据文件中的异常数据。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
第一线性拟合模块，分别对SPS文件中每条检波点的横纵坐标与桩号分别进行线性拟合，得到每条检波点线的检波点横纵坐标-桩号的线性转换关系。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
规则化模块，根据预设正负偏移距分布范围对地震初至数据文件中距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行规则化。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
时间校正模块，根据所述极小值处的检波点附近地震道的拟合初至与实际地震初至的时差对所述极小值处的检波点的地震初至时间进行校正。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
第二线性拟合模块，对所述炮点对应的每条检波点线的线号与对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
第三线性拟合模块，对所述炮点对应的每条检波点线对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
二次定位模块，获取偏移量大于预设误差门槛的炮点编号，以根据所述炮点编号进行二次测量得到实际的横纵坐标，并利用实际的横纵坐标替换SPS文件中所述炮点的横纵坐标。
进一步地，陆上三维炮点二次定位及定量偏移监控装置还包括：
文件修正模块，利用所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标替换SPS文件中所述炮点的横纵坐标。
第三方面，提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时上述的陆上三维炮点二次定位及定量偏移监控方法的步骤。
第四方面，提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述的陆上三维炮点二次定位及定量偏移监控方法的步骤。
本发明提供一种陆上三维炮点二次定位及定量偏移监控方法和装置、电子设备以及计算机可读存储介质，该方法包括：对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间；根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标；分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号；根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量，能够自动实现炮点偏移的定量监控，有效降低了工作量，防止人工误差，降低生产成本，满足地震采集项目高效质量监控和定量监控的需求，另外，还能够对偏移炮点进行精确二次定位，供地震队复测，满足地震队对采集质量的监控要求。
为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中的服务器S1与客户端设备B1之间的架构示意图；
图2为本发明实施例中的服务器S1、客户端设备B1及数据库服务器S2之间的架构示意图；
图3是本发明实施例中的陆上三维炮点二次定位及定量偏移监控方法的流程示意图一；
图4是本发明实施例中的陆上三维炮点二次定位及定量偏移监控方法的流程示意图二；
图5是本发明实施例中的陆上三维炮点二次定位及定量偏移监控方法的流程示意图三；
图6示出了LGB地区三维某炮点排列规则化前的拟合曲线及极值点位置；
图7示出了图6中的数据应用本发明实施例提供的陆上三维炮点二次定位及定量偏移监控方法的拟合曲线及极值点位置；
图8示出了LGB地区三维某炮点排列常规方法的拟合曲线及极值点位置；
图9示出了LGB地区三维某炮点采用本发明实施例提供的陆上三维炮点二次定位及定量偏移监控方法校正后的拟合曲线及极值点位置；
图10示出了LGB地区采用常规方法与本发明实施例提供的陆上三维炮点二次定位及定量偏移监控方法监控1000炮炮点偏移的时间对比图；
图11是本发明实施例中的陆上三维炮点二次定位及定量偏移监控装置的结构框图；
图12为本发明实施例电子设备的结构图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
目前采用人工核对的方法进行炮点偏移检查，导致工作量巨大，而且有时因疏忽造成差错。
为解决现有技术中的上述技术问题，本发明实施例提供一种陆上三维炮点二次定位及定量偏移监控方法，能够自动实现炮点偏移的定量监控，实现了对整个地震采集质量的监控，有效降低了工作量，防止人工误差，降低生产成本，满足地震采集项目高效质量监控和定量监控的需求，另外，还能够对偏移炮点进行精确二次定位，供地震队复测，满足地震队对采集质量的监控要求，并且满足地震资料处理的时效要求，同时也可以将拟合得到的坐标直接用于地震资料处理，以确保地震资料的成像质量。
有鉴于此，本申请提供了一种陆上三维炮点二次定位及定量偏移监控装置，该装置可以为一种服务器S1，参见图1，该服务器S1可以与至少一个客户端设备B1通信连接，所述客户端设备B1可以将地震初至数据文件或SPS文件发送至所述服务器S1，所述服务器S1可以在线接收所述地震初至数据文件或SPS文件。所述服务器S1可以在线或者离线对获取的地震初至数据文件或SPS文件进行预处理，对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间；根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标；分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号；根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量。而后，所述服务器S1可以将述炮点的偏移量在线发送至所述客户端设备B1。所述客户端设备B1可以在线接收所述述炮点的偏移量。
另外，参见图2，所述服务器S1还可以与至少一个数据库服务器S2通信连接，所述数据库服务器S2用于存储历史地震初至数据文件或SPS文件。所述数据库服务器S2在线将所述历史地震初至数据文件或SPS文件发送至所述服务器S1，所述服务器S1可以在线接收所述历史地震初至数据文件或SPS文件，而后根据多个历史地震初至数据文件或SPS文件进行陆上三维炮点二次定位及定量偏移监控。
可以理解的是，所述客户端设备B1可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备、智能穿戴设备等。其中，所述智能穿戴设备可以包括智能眼镜、智能手表、智能手环等。
在实际应用中，进行陆上三维炮点二次定位及定量偏移监控的部分可以在如上述内容所述的服务器S1侧执行，即，如图1所示的架构，也可以所有的操作都在所述客户端设备B1中完成，且该所述客户端设备B1可以直接与数据库服务器S2进行通信连接。具体可以根据所述客户端设备B1的处理能力，以及用户使用场景的限制等进行选择。本申请对此不作限定。若所有的操作都在所述客户端设备B1中完成，所述客户端设备B1还可以包括处理器，用于进行陆上三维炮点二次定位及定量偏移监控的具体处理。
所述服务器与所述客户端设备之间可以使用任何合适的网络协议进行通信，包括在本申请提交日尚未开发出的网络协议。所述网络协议例如可以包括TCP/IP协议、UDP/IP协议、HTTP协议、HTTPS协议等。当然，所述网络协议例如还可以包括在上述协议之上使用的RPC协议(Remote Procedure Call Protocol，远程过程调用协议)、REST协议(Representational State Transfer，表述性状态转移协议)等。
图3是本发明实施例中的陆上三维炮点二次定位及定量偏移监控方法的流程示意图一。如图3所示，该陆上三维炮点二次定位及定量偏移监控方法可以包括以下内容：
步骤S100：对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间。
具体地，炮点对应的检波点线中均包括用于采集其产生的地震波的检波点，一个炮点通常对应多条检波点线。
另外，该预定偏移距范围根据实际工区的地质环境以及陆上三维炮点二次定位及定量偏移监控精度进行设置，例如，该预定偏移距范围可为2000～8000米，例如4000米或5000米。
其中，采用二次函数对距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合。二次曲线拟合可采用如下公式：
t＝ax2+bx+c
其中：t代表拟合的地震初至时间；
a是计算的二次项拟合系数；
b是计算的一次项拟合系数；
c是计算的常数项拟合系数；
x代表检波点桩号；
上述二次曲线的极小值处的检波点桩号和初至时间分别按下式计算：
Figure BDA0002131419750000101
Figure BDA0002131419750000102
其中：
xmin是拟合曲线极小值对应的桩号；
ttmin是拟合曲线极小值对应的地震初至时间；
步骤S200：根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标。
步骤S300：分别对所述炮点对应的每条检波点线的线号与该检波点线对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号。
具体地，将所述炮点的每条检波点线的线号与该检波点线对应的极小值处的检波点的地震初至时间构成一组数据对，然后采用二次曲线对所述炮点对应的多条检波点线形成的多个数据对进行拟合，并计算该二次曲线的极小值，得到该曲线极小值处对应的检波点线号和地震波初至时间；
二次拟合采用如下公式：
t＝ax2+bx+c
其中：t代表拟合的地震初至时间；
a是计算的二次项拟合系数；
b是计算的一次项拟合系数；
c是计算的常数项拟合系数；
x代表检波点线号；
上述二次曲线的极小值处的检波点线号按下式计算：
Figure BDA0002131419750000103
其中：
xmin是拟合二次曲线极小值对应的检波点线号；
步骤S400：根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；
步骤S500：根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；
步骤S600：根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
其中，该炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标即为拟合得到的该炮点的定位坐标，计算该定位坐标与SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
计算公式为：
Figure BDA0002131419750000111
其中：D为计算的炮点偏移量；
xn为所述炮点对应的极小值处的检波点线的横坐标；
xs为SPS文件中的炮点横坐标；
yn为该炮点对应的极小值处的检波点线的纵坐标；
ys为SPS文件中的炮点纵坐标；
通过采用上述技术方案，能够自动实现炮点偏移的定量监控，实现了对整个地震采集质量的监控，而且杜甫块，有效降低了工作量，防止人工误差，降低生产成本，满足地震采集项目高效质量监控和定量监控的需求，另外，还能够对偏移炮点进行精确二次定位，供地震队复测，满足地震队对采集质量的监控要求，并且满足地震资料处理的时效要求，同时也可以将拟合得到的坐标直接用于地震资料处理，以确保地震资料的成像质量。
在一个可选的实施例中，参见图4，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
步骤S10：根据SPS文件获取地震初至数据文件，所述地震初至数据文件包括多个炮点的地震初至数据，每个炮点的地震初至数据均包括：用于接收所述炮点产生的地震波的检波点桩号、地震初至时间、横纵坐标及其所属的检波点线。
具体地，根据野外生产整理工区的三维的SPS文件，即检波点文件、炮点文件和关系文件。检波点文件按检波点桩号升序排序、炮点文件按炮点桩号升序排序；然后拾取每炮的地震初至时间，得到地震初至数据文件。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
步骤S20：对所述地震初至数据文件进行基准面静校正得到校正后的地震初至数据文件，以利用校正后的地震初至文件进行陆上三维炮点二次定位及定量偏移监控。
其中，对地震勘探数据进行基准面静校正属于本领域常用的技术手段，在此不再赘述。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
步骤S30：采用多项式拟合方式剔除所述地震初至数据文件中的异常数据。
其中，可采用人机交互和多项式拟合方法剔除异常地震初至数据。
在一个可选的实施例中，参见图5，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
步骤S40：分别对SPS文件中每条检波点的横纵坐标与桩号分别进行线性拟合，得到每条检波点线的检波点横纵坐标-桩号的线性转换关系。
其中，该检波点横纵坐标-桩号的线性转换关系包括横坐标与桩号的线性转换关系和纵坐标与桩号的线性转换关系。
线性转换公式如下：
X＝a×R+b
其中：X代表横(纵)坐标；
a是计算的一次项拟合系数；
b是计算的常数项拟合系数；
R代表检波点桩号。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
根据预设正负偏移距分布范围对地震初至数据文件中距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行规则化。
即通过规则化证正负偏移距的最大绝对值基本相等。
具体地，以正负偏移距的最大绝对值中较小的偏移距为标准，将超过该偏移距预设值的地震初至数据舍弃。其中，该预设值根据实际需要选取，例如20米～120米，如40米、60米、90米等。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
根据所述极小值处的检波点附近地震道的拟合初至与实际地震初至的时差对所述极小值处的检波点的地震初至时间进行校正。
具体地，计算每条检波点线极小值处检波点两边附近几个地震道(例如3道、5道)的拟合初至与实际地震初至的误差，计算它们的平均值(或采用数据拟合得到极值点处的误差拟合值)作为系统极值误差，对极小值处的检波点的地震初至时间进行校正。
其中，系统极值误差采用如下公式计算：
Figure BDA0002131419750000131
其中：
wt为某条检波点线的拟合系统极值误差；
n是数据点数，这里n为5；
nti是检波点编号i的拟合初至；
ti是检波点编号i的地震初至；
i是数据点编号，取值范围为1到n。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
对所述炮点对应的每条检波点线的线号与对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系。
其中，线性拟合的方法参见上述，在此不再赘述。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
对所述炮点对应的每条检波点线对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系。
其中，线性拟合的方法参见上述，在此不再赘述。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
获取偏移量大于预设误差门槛(可根据实际需要任意设置，例如20米，或者，一般默认为1/2道距)的炮点编号，以根据所述炮点编号进行二次测量得到实际的横纵坐标，并利用实际的横纵坐标替换SPS文件中所述炮点的横纵坐标。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控方法还可以包括：
利用所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标替换SPS文件中所述炮点的横纵坐标。
其中，通过重复上述步骤可以对所有炮点的偏移进行定量监控计算和二次定位，最后得到炮点的新的SPS文件，用于地震资料处理。
通过在LGB地区应用本发明实施例提供的陆上三维炮点二次定位及定量偏移监控方法，LGB地区三维某炮点排列规则化前的拟合曲线及极值点位置如图6所示，应用本发明实施例提供的陆上三维炮点二次定位及定量偏移监控方法的拟合曲线及极值点位置如图7所示；LGB地区三维某炮点排列常规方法的拟合曲线及极值点位置如图8所示；LGB地区三维某炮点采用本发明实施例提供的陆上三维炮点二次定位及定量偏移监控方法校正后的拟合曲线及极值点位置如图9所示LGB地区采用常规方法与本发明实施例提供的陆上三维炮点二次定位及定量偏移监控方法监控1000炮炮点偏移的时间对比图如图10所示，通过对比可以发现，采用本发明实施例提供的陆上三维炮点二次定位及定量偏移监控方法校正后的拟合曲线及极值点位置与采用常规方法相比，提高了数据精度，另外，采用常规方法需要14小时，而采用本方法只需要20分钟(0.33小时)，炮偏检查速度提高了42倍，并且对偏移的炮点能进行二次定位，供野外复测或直接作为炮点坐标信息提供资料处理，从而很好满足了三维采集质量监控及资料整理的要求。
基于同一发明构思，本申请实施例还提供了一种陆上三维炮点二次定位及定量偏移监控装置，可以用于实现上述实施例所描述的方法，如下面的实施例所述。由于陆上三维炮点二次定位及定量偏移监控装置解决问题的原理与上述方法相似，因此陆上三维炮点二次定位及定量偏移监控装置的实施可以参见上述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图11是本发明实施例中的陆上三维炮点二次定位及定量偏移监控装置的结构框图一。如图11所示，该的陆上三维炮点二次定位及定量偏移监控装置具体包括：第一处理模块10、第一转换模块20、第二处理模块30、第二转换模块40、第三转换模块50以及偏移量计算模块60。
第一处理模块10对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间；
具体地，炮点对应的检波点线中均包括用于采集其产生的地震波的检波点，一个炮点通常对应多条检波点线。
另外，该预定偏移距范围根据实际工区的地质环境以及陆上三维炮点二次定位及定量偏移监控精度进行设置，例如，该预定偏移距范围可为2000～8000米，例如4000米或5000米。
其中，采用二次函数对距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合。二次曲线拟合可采用如下公式：
t＝ax2+bx+c
其中：t代表拟合的地震初至时间；
a是计算的二次项拟合系数；
b是计算的一次项拟合系数；
c是计算的常数项拟合系数；
x代表检波点桩号；
上述二次曲线的极小值处的检波点桩号和初至时间分别按下式计算：
Figure BDA0002131419750000161
Figure BDA0002131419750000162
其中：
xmin是拟合曲线极小值对应的桩号；
tmin是拟合曲线极小值对应的地震初至时间；
第一转换模块20根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标；
第二处理模块30分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号；
具体地，将所述炮点的每条检波点线的线号与该检波点线对应的极小值处的检波点的地震初至时间构成一组数据对，然后采用二次曲线对所述炮点对应的多条检波点线形成的多个数据对进行拟合，并计算该二次曲线的极小值，得到该曲线极小值处对应的检波点线号和地震波初至时间；
二次拟合采用如下公式：
t＝ax2+bx+c
其中：t代表拟合的地震初至时间；
a是计算的二次项拟合系数；
b是计算的一次项拟合系数；
c是计算的常数项拟合系数；
x代表检波点线号；
上述二次曲线的极小值处的检波点线号按下式计算：
Figure BDA0002131419750000163
其中：
xmin是拟合二次曲线极小值对应的检波点线号；
第二转换模块40根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；
第三转换模块50根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；
偏移量计算模块60根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
其中，该炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标即为拟合得到的该炮点的定位坐标，计算该定位坐标与SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
计算公式为：
Figure BDA0002131419750000171
其中：D为计算的炮点偏移量；
xn为所述炮点对应的极小值处的检波点线的横坐标；
xs为SPS文件中的炮点横坐标；
yn为该炮点对应的极小值处的检波点线的纵坐标；
ys为SPS文件中的炮点纵坐标；
通过采用上述技术方案，能够自动实现炮点偏移的定量监控，实现了对整个地震采集质量的监控，而且杜甫块，有效降低了工作量，防止人工误差，降低生产成本，满足地震采集项目高效质量监控和定量监控的需求，另外，还能够对偏移炮点进行精确二次定位，供地震队复测，满足地震队对采集质量的监控要求，并且满足地震资料处理的时效要求，同时也可以将拟合得到的坐标直接用于地震资料处理，以确保地震资料的成像质量。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：数据获取模块，根据SPS文件获取地震初至数据文件，所述地震初至数据文件包括多个炮点的地震初至数据，每个炮点的地震初至数据均包括：用于接收所述炮点产生的地震波的检波点桩号、地震初至时间、横纵坐标及其所属的检波点线。
具体地，根据野外生产整理工区的三维的SPS文件，即检波点文件、炮点文件和关系文件。检波点文件按检波点桩号升序排序、炮点文件按炮点桩号升序排序；然后拾取每炮的地震初至时间，得到地震初至数据文件。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：静校正模块，对所述地震初至数据文件进行基准面静校正得到校正后的地震初至数据文件，以利用校正后的地震初至文件进行陆上三维炮点二次定位及定量偏移监控。
其中，对地震勘探数据进行基准面静校正属于本领域常用的技术手段，在此不再赘述。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：异常数据剔除模块，采用多项式拟合方式剔除所述地震初至数据文件中的异常数据。
其中，可采用人机交互和多项式拟合方法剔除异常地震初至数据。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：第一线性拟合模块，分别对SPS文件中每条检波点的横纵坐标与桩号分别进行线性拟合，得到每条检波点线的检波点横纵坐标-桩号的线性转换关系。
其中，该检波点横纵坐标-桩号的线性转换关系包括横坐标与桩号的线性转换关系和纵坐标与桩号的线性转换关系。
线性转换公式如下：
X＝a×R+b
其中：X代表横(纵)坐标；
a是计算的一次项拟合系数；
b是计算的常数项拟合系数；
R代表检波点桩号。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：规则化模块，根据预设正负偏移距分布范围对地震初至数据文件中距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行规则化。
即通过规则化证正负偏移距的最大绝对值基本相等。
具体地，以正负偏移距的最大绝对值中较小的偏移距为标准，将超过该偏移距预设值的地震初至数据舍弃。其中，该预设值根据实际需要选取，例如20米～120米，如40米、60米、90米等。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：时间校正模块，根据所述极小值处的检波点附近地震道的拟合初至与实际地震初至的时差对所述极小值处的检波点的地震初至时间进行校正。
具体地，计算每条检波点线极小值处检波点两边附近几个地震道(例如3道、5道)的拟合初至与实际地震初至的误差，计算它们的平均值(或采用数据拟合得到极值点处的误差拟合值)作为系统极值误差，对极小值处的检波点的地震初至时间进行校正。
其中，系统极值误差采用如下公式计算：
Figure BDA0002131419750000191
其中：
wt为某条检波点线的拟合系统极值误差；
n是数据点数，这里n为5；
nti是检波点编号i的拟合初至；
ti是检波点编号i的地震初至；
i是数据点编号，取值范围为1到n。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：第二线性拟合模块，对所述炮点对应的每条检波点线的线号与对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系。
其中，线性拟合的方法参见上述，在此不再赘述。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：第三线性拟合模块，对所述炮点对应的每条检波点线对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系。
其中，线性拟合的方法参见上述，在此不再赘述。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：二次定位模块，获取偏移量大于预设误差门槛(可根据实际需要任意设置，例如20米，或者，一般默认为1/2道距)的炮点编号，以根据所述炮点编号进行二次测量得到实际的横纵坐标，并利用实际的横纵坐标替换SPS文件中所述炮点的横纵坐标。
在一个可选的实施例中，该陆上三维炮点二次定位及定量偏移监控装置还可以包括：文件修正模块，利用所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标替换SPS文件中所述炮点的横纵坐标。
其中，通过上述各模块配合可以对所有炮点的偏移进行定量监控计算和二次定位，最后得到炮点的新的SPS文件，用于地震资料处理。
上述实施例阐明的装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为电子设备，具体的，电子设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中电子设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现下述步骤：
对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间；
根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标；
分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号；
根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；
根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；
根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
从上述描述可知，本发明实施例提供的电子设备，可用于监控炮点偏移，能够自动实现炮点偏移的定量监控，有效降低了工作量，防止人工误差，降低生产成本，满足地震采集项目高效质量监控和定量监控的需求，另外，还能够对偏移炮点进行精确二次定位，供地震队复测，满足地震队对采集质量的监控要求。
下面参考图12，其示出了适于用来实现本申请实施例的电子设备600的结构示意图。
如图12所示，电子设备600包括中央处理单元(CPU)601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中，还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡，调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口606。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现下述步骤：
对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间；
根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标；
分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号；
根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；
根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；
根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
从上述描述可知，本发明实施例提供的计算机可读存储介质，可用于监控炮点偏移，能够自动实现炮点偏移的定量监控，有效降低了工作量，防止人工误差，降低生产成本，满足地震采集项目高效质量监控和定量监控的需求，另外，还能够对偏移炮点进行精确二次定位，供地震队复测，满足地震队对采集质量的监控要求。
在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (24)

1.一种陆上三维炮点二次定位及定量偏移监控方法，其特征在于，包括：
对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间；
根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标；
分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号；
根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；
根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；
根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
2.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，还包括：
根据SPS文件获取地震初至数据文件，所述地震初至数据文件包括多个炮点的地震初至数据，每个炮点的地震初至数据均包括：用于接收所述炮点产生的地震波的检波点桩号、地震初至时间、横纵坐标及其所属的检波点线。
3.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，还包括：
对所述地震初至数据文件进行基准面静校正得到校正后的地震初至数据文件，以利用校正后的地震初至文件进行陆上三维炮点二次定位及定量偏移监控。
4.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，还包括：
采用多项式拟合方式剔除所述地震初至数据文件中的异常数据。
5.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，还包括：
分别对SPS文件中每条检波点的横纵坐标与桩号分别进行线性拟合，得到每条检波点线的检波点横纵坐标-桩号的线性转换关系。
6.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，所述对地震初至数据文件中距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合之前，还包括：
根据预设正负偏移距分布范围对地震初至数据文件中距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行规则化。
7.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，所述根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标之前，还包括：
根据所述极小值处的检波点附近地震道的拟合初至与实际地震初至的时差对所述极小值处的检波点的地震初至时间进行校正。
8.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，还包括：
对所述炮点对应的每条检波点线的线号与对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系。
9.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，还包括：
对所述炮点对应的每条检波点线对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系。
10.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，还包括：
获取偏移量大于预设误差门槛的炮点编号，以根据所述炮点编号进行二次测量得到实际的横纵坐标，并利用实际的横纵坐标替换SPS文件中所述炮点的横纵坐标。
11.根据权利要求1所述的陆上三维炮点二次定位及定量偏移监控方法，其特征在于，还包括：
利用所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标替换SPS文件中所述炮点的横纵坐标。
12.一种陆上三维炮点二次定位及定量偏移监控装置，其特征在于，包括：
第一处理模块，对地震初至数据文件中距离一炮点对应的一检波点线上的检波点预定偏移距范围内的检波点的桩号及地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点桩号及地震初至时间；
第一转换模块，根据所述检波点线对应的极小值处的检波点桩号、预获取的所述检波点线的检波点横纵坐标-桩号的线性转换关系，得到所述检波点线对应的极小值处的检波点的横纵坐标；
第二处理模块，分别对所述炮点对应的每条检波点线的线号与对应的极小值处的检波点的地震初至时间进行二次曲线拟合，并获取拟合得到的二次曲线的极小值处的检波点线的线号；
第二转换模块，根据所述炮点对应的极小值处的检波点线的线号、预获取的所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系，得到所述炮点对应的极小值处的检波点线的纵坐标；
第三转换模块，根据所述炮点对应的极小值处的检波点线的纵坐标、预获取的所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系，得到所述炮点对应的极小值处的检波点线的横坐标；
偏移量计算模块，根据所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标以及SPS文件中的所述炮点的坐标得到所述炮点的偏移量。
13.根据权利要求12所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
数据获取模块，根据SPS文件获取地震初至数据文件，所述地震初至数据文件包括多个炮点的地震初至数据，每个炮点的地震初至数据均包括：用于接收所述炮点产生的地震波的检波点桩号、地震初至时间、横纵坐标及其所属的检波点线。
14.根据权利要求12所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
静校正模块，对所述地震初至数据文件进行基准面静校正得到校正后的地震初至数据文件，以利用校正后的地震初至文件进行陆上三维炮点二次定位及定量偏移监控。
15.根据权利要求12所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
异常数据剔除模块，采用多项式拟合方式剔除所述地震初至数据文件中的异常数据。
16.根据权利要求12所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
第一线性拟合模块，分别对SPS文件中每条检波点的横纵坐标与桩号分别进行线性拟合，得到每条检波点线的检波点横纵坐标-桩号的线性转换关系。
17.根据权利要求12所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
规则化模块，根据预设正负偏移距分布范围对地震初至数据文件中距离一炮点对应的一检波点线预定偏移距范围内的检波点的桩号及地震初至时间进行规则化。
18.根据权利要求12所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
时间校正模块，根据所述极小值处的检波点附近地震道的拟合初至与实际地震初的时差至对所述极小值处的检波点的地震初至时间进行校正。
19.根据权利要求12所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
第二线性拟合模块，对所述炮点对应的每条检波点线的线号与对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的检波点线的线号-极小值处的检波点的纵坐标的线性转换关系。
20.根据权利要求12所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
第三线性拟合模块，对所述炮点对应的每条检波点线对应的所述极小值处的检波点的横纵坐标进行线性拟合得到所述炮点的每个检波点线对应的所述极小值处的检波点的横纵坐标线性转换关系。
21.根据权利要求12所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
二次定位模块，获取偏移量大于预设误差门槛的炮点编号，以根据所述炮点编号进行二次测量得到实际的横纵坐标，并利用实际的横纵坐标替换SPS文件中所述炮点的横纵坐标。
22.根据权利要求13所述的陆上三维炮点二次定位及定量偏移监控装置，其特征在于，还包括：
文件修正模块，利用所述炮点对应的极小值处的检波点线的纵坐标、所述炮点对应的极小值处的检波点线的横坐标替换SPS文件中所述炮点的横纵坐标。
23.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至11任一项所述的陆上三维炮点二次定位及定量偏移监控方法的步骤。
24.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至11任一项所述的陆上三维炮点二次定位及定量偏移监控方法的步骤。
CN201910639630.6A 2019-07-16 2019-07-16 陆上三维炮点二次定位及定量偏移监控方法和装置 Active CN112241023B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910639630.6A CN112241023B (zh) 2019-07-16 2019-07-16 陆上三维炮点二次定位及定量偏移监控方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910639630.6A CN112241023B (zh) 2019-07-16 2019-07-16 陆上三维炮点二次定位及定量偏移监控方法和装置

Publications (2)

Publication Number Publication Date
CN112241023A true CN112241023A (zh) 2021-01-19
CN112241023B CN112241023B (zh) 2023-04-25

Family

ID=74166630

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910639630.6A Active CN112241023B (zh) 2019-07-16 2019-07-16 陆上三维炮点二次定位及定量偏移监控方法和装置

Country Status (1)

Country Link
CN (1) CN112241023B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5757722A (en) * 1997-05-21 1998-05-26 Western Atlas International, Inc. Method for verifying the location of an array of detectors
US20020075759A1 (en) * 2001-01-21 2002-06-20 Xixiang Zhou Static correction method for exploration seismic data using first arrivals of seismic waves
CN102968404A (zh) * 2012-11-02 2013-03-13 中国石油天然气股份有限公司 地震资料的初至信息处理方法及装置
CN104375181A (zh) * 2014-11-24 2015-02-25 中国石油天然气集团公司 一种快速生成三维共检波点初至道集的方法
CN104483705A (zh) * 2014-12-31 2015-04-01 中国石油天然气集团公司 一种三维剩余静校正方法
CN106324671A (zh) * 2015-07-01 2017-01-11 中国石油天然气股份有限公司 一种检查炮点偏移的方法及装置
CN108107471A (zh) * 2017-11-17 2018-06-01 中国石油天然气集团公司 一种分方位初至数据体的获取方法及装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5757722A (en) * 1997-05-21 1998-05-26 Western Atlas International, Inc. Method for verifying the location of an array of detectors
US20020075759A1 (en) * 2001-01-21 2002-06-20 Xixiang Zhou Static correction method for exploration seismic data using first arrivals of seismic waves
CN102968404A (zh) * 2012-11-02 2013-03-13 中国石油天然气股份有限公司 地震资料的初至信息处理方法及装置
CN104375181A (zh) * 2014-11-24 2015-02-25 中国石油天然气集团公司 一种快速生成三维共检波点初至道集的方法
CN104483705A (zh) * 2014-12-31 2015-04-01 中国石油天然气集团公司 一种三维剩余静校正方法
CN106324671A (zh) * 2015-07-01 2017-01-11 中国石油天然气股份有限公司 一种检查炮点偏移的方法及装置
CN108107471A (zh) * 2017-11-17 2018-06-01 中国石油天然气集团公司 一种分方位初至数据体的获取方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
FENG FAQUAN等: "A case study of near-surface modeling and high-precision statics in huge-thick loess tableland" *
冯发全等: "陆上三维炮偏定量监控方法及应用实例" *
孙哲等: "地震采集异常道实时监控技术与应用" *
蔡加铭等: "兴隆台城区三维炮点偏移校正技术" *
马强;李东安;: "SPS文件恢复与地震资料重处理" *

Also Published As

Publication number Publication date
CN112241023B (zh) 2023-04-25

Similar Documents

Publication Publication Date Title
CN114706933A (zh) 基于无人机作业的喷洒区域确定方法及装置
CN109029368B (zh) 像方补偿的遥感影像/sar影像高精度几何定位后处理方法
CN110197615B (zh) 用于生成地图的方法及装置
CN107527382B (zh) 数据处理方法以及装置
Petersen et al. Automated quality control for large seismic networks: Implementation and application to the AlpArray seismic network
CN103499834B (zh) 恢复模拟检波器地震数据低频信号的方法和装置
CN106054252B (zh) 一种叠前时间偏移的方法及装置
CN111781635B (zh) 海底四分量弹性波高斯束深度偏移方法和装置
Xu et al. Systemic comparison of seismometer horizontal orientations based on teleseismic earthquakes and ambient‐noise data
CN114463932A (zh) 非接触式施工安全距离主动动态识别预警系统及方法
CN109782355A (zh) Obs检波点漂移的检测方法及装置
CN112241023B (zh) 陆上三维炮点二次定位及定量偏移监控方法和装置
CN112379436A (zh) 混合震源激发地震资料静校正量计算方法及装置
CN113568033B (zh) 一种三维非规则采样地震采集观测系统设计方法和装置
CN111435171B (zh) 海底节点地震数据的时钟漂移校正方法、装置及存储介质
CN111337973B (zh) 地震数据重建方法、系统
CN112946736A (zh) 一种三维观测系统重建方法及系统
Banaś A review of robust estimation methods applied in surveying
CN106094035B (zh) 一种基于单机快速处理的剩余静校正计算方法
CN105928543B (zh) 一种仿生偏振导航载体测量误差的测量分析方法
CN112241022A (zh) 基于射线密度生成层析反演模型速度界面的方法和装置
CN111965729A (zh) 可控震源组合中心实时监控方法、系统及装置
CN111045090B (zh) 一种磁异常网格快速缝合方法
CN109613458B (zh) 航天器磁异常定位方法
CN111596346B (zh) 弹性波速度反演方法和装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant