CN112946736B - 一种三维观测系统重建方法及系统 - Google Patents

一种三维观测系统重建方法及系统 Download PDF

Info

Publication number
CN112946736B
CN112946736B CN201911171968.XA CN201911171968A CN112946736B CN 112946736 B CN112946736 B CN 112946736B CN 201911171968 A CN201911171968 A CN 201911171968A CN 112946736 B CN112946736 B CN 112946736B
Authority
CN
China
Prior art keywords
work area
information
detection
point
line
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.)
Active
Application number
CN201911171968.XA
Other languages
English (en)
Other versions
CN112946736A (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 CN201911171968.XA priority Critical patent/CN112946736B/zh
Publication of CN112946736A publication Critical patent/CN112946736A/zh
Application granted granted Critical
Publication of CN112946736B publication Critical patent/CN112946736B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
    • 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
    • G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/64—Geostructures, e.g. in 3D data cubes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/70—Other details related to processing
    • G01V2210/74—Visualisation of seismic data

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)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明提供了一种三维观测系统重建方法及系统，所述方法包括：根据工区范围参数确定所述工区的工区大小及工区方位角信息；通过坐标旋转公式和所述工区方位角信息将所述工区的方位角转换为直角得到直角工区；通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息；通过坐标反旋转公式将所述直角工区炮检点信息转换为实际炮检点信息，根据所述实际炮检点信息进行地震资料分析，本发明可能快速获得在给定参数范围内的全新观测系统信息，从而灵活满足地震资料处理过程中数据规则化、多维数据插值与地震道数据恢复环节中的数据准备。

Description

一种三维观测系统重建方法及系统
技术领域
本发明涉及石油勘探技术领域，尤其涉及一种三维观测系统重建方法及系统。
背景技术
地震观测系统是通过野外采集时产生的地震数据赋予实际空间关系。具体包括两大步骤：第一是直接利用野外采集的SPS文件(炮检点坐标信息以及每一炮对应的不同检波点接收关系)通过计算各种属性库，然后将这些属性库和地震数据匹配，匹配主要通过野外文件号和通道号来唯一匹配，从而使地震数据具有真实的空间关系。
但是，随着勘探力度的加大，勘探工区向山区、戈壁等地表复杂的区域进行，在许多地表恶劣的区域不得不采用蔽障布设的手段来实现观测点的布设。而以此布设的观测点与理论设计的位置可能相差较大。且这种位置的偏差使得资料处理过程难度加大，对资料的成像质量降低。
发明内容
本发明的一个目的在于提供一种三维观测系统重建方法，能快速获得在给定参数范围内的全新观测系统信息，从而灵活满足地震资料处理过程中数据规则化、多维数据插值与地震道数据恢复环节中的数据准备。本发明的另一个目的在于提供一种三维观测系统重建系统。本发明的另一个目的在于提供一种计算机设备。本发明的另一个目的在于提供一种可读介质。
为了达到以上目的，本发明一方面公开了一种三维观测系统重建方法，包括：
根据工区范围参数确定所述工区的工区大小及工区方位角信息；
通过坐标旋转公式和所述工区方位角信息将所述工区的方位角转换为直角得到直角工区；
通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息；
通过坐标反旋转公式将所述直角工区炮检点信息转换为实际炮检点信息，根据所述实际炮检点信息进行地震资料分析。
优选地，进一步包括：
根据给定参数和所述实际炮检点信息建立炮检关系。
优选地，所述炮检点参数信息包括直角工区中炮线距、检波线距、炮点点距和检波点间距。
优选地，所述直角工区炮检点信息包括炮检点的点数、线数和坐标信息。
优选地，所述给定参数包括每个炮点对应的检波线数量和每一条检波线的接收通道数量。
优选地，所述工区的炮点位于所述工区的中央。
优选地，所述通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息具体包括以下步骤的至少之一：
根据直角工区与检波线平行的工区范围信息和检波点距得到每条检波线的检波点点数；
根据直角工区与检波线垂直的工区范围信息和检波点线距得到检波点线数。
优选地，所述通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息具体还包括以下步骤的至少之一：
根据直角工区与检波线平行的工区范围信息、检波点距和炮线距得到炮点线数；
根据检波线数量、检波点线距和炮点距得到每条炮线炮点点数；
根据炮点线数和每条炮线炮点点数得到炮点总数；
根据每条检波线的检波点点数和检波点线数得到检波点总数；
根据直角工区与检波线平行的工区原点坐标、一个炮点的线号和线号距得到所述一个炮点与检波线平行的直角工区坐标；
根据直角工区与检波线垂直的工区原点坐标、一个炮点的点号和点号距得到所述一个炮点与检波线垂直的直角工区坐标；
根据直角工区与检波线平行的工区原点坐标、一个检波点的点号和点号距得到所述一个检波点与检波线平行的直角工区坐标；
根据直角工区与检波线垂直的工区原点坐标、一个检波点的线号和线号距得到所述一个检波点与检波线垂直的直角工区坐标。
本发明还公开了一种三维观测系统重建系统，包括：
工区确定单元，用于根据工区范围参数确定所述工区的工区大小及工区方位角信息；
工区转换单元，用于通过坐标旋转公式和所述工区方位角信息将所述工区的方位角转换为直角得到直角工区；
炮检点确定单元，用于通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息；
炮检点转换单元，用于通过坐标反旋转公式将所述直角工区炮检点信息转换为实际炮检点信息，根据所述实际炮检点信息进行地震资料分析。
本发明还公开了一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，
所述处理器执行所述程序时实现如上所述方法。
本发明还公开了一种计算机可读介质，其上存储有计算机程序，
该程序被处理器执行时实现如上所述方法。
本发明提供了一种简便有效的三维观测系统重建方案，在三维观测系统重建过程中利用已有的坐标旋转公式将带方位角的工区信息旋转成水平正交的信息，免去了角度对后续计算坐标的精度与计算效率影响。通过坐标计算公式计算完所有炮检点信息后再进行角度反旋转，从而得到炮、检点真实坐标信息。在确定任意炮点的接收范围时，将任意炮与给定接收范围的检波点信息进行匹配可快速确定。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出本发明一种三维观测系统重建方法一个具体实施例的流程图之一；
图2示出本发明一种三维观测系统重建方法一个具体实施例的流程图之二；
图3示出本发明一种三维观测系统重建方法一个具体例子中的工区示意图；
图4示出本发明一种三维观测系统重建方法一个具体例子中的直角工区示意图；
图5示出本发明一种三维观测系统重建方法一个具体例子中的直角工区调整后的示意图；
图6示出本发明一种三维观测系统重建方法一个具体例子中的直角工区放炮方式布设示意图；
图7示出本发明一种三维观测系统重建方法一个具体例子中新观测系统示意图之一；
图8示出本发明一种三维观测系统重建方法一个具体例子中新观测系统示意图之二；
图9示出本发明一种三维观测系统重建方法一个具体例子中炮点纵偏距和非纵偏距的示意图；
图10示出本发明一种三维观测系统重建方法一个具体例子中检波点接收范围的示意图；
图11示出本发明一种三维观测系统重建系统一个具体实施例的结构图之一；
图12示出本发明一种三维观测系统重建系统一个具体实施例的结构图之二；
图13示出适于用来实现本发明实施例的计算机设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
根据本发明的一个方面，本实施例公开了一种三维观测系统重建方法。如图1所示，本实施例中，三维观测系统重建方法具体包括：
S100：根据工区范围参数确定所述工区的工区大小及工区方位角信息。
优选的，工区范围参数可包括工区的四个角点坐标等信息，其中坐标可为预先设置的直角坐标系，本实施例中，将与检波线平行的方向定义为x方向，与x垂直的方向为y方向。根据工区的四个角点坐标可确定所述工区的工区大小以及工区方位角信息。
S200：通过坐标旋转公式和所述工区方位角信息将所述工区的方位角转换为直角得到直角工区。
优选的，可选择旋转围绕的旋转0点(旋转不动点)为四个角点坐标中y坐标最小的点。通过坐标旋转公式和所述工区方位角信息可将所述工区的方位角转换为直角得到直角工区，通过方位角转换可使得后续计算各信息如各炮检点新的坐标不需考虑夹角因素，避免了由此产生的精度误差，同时也提升了计算效率。
在一个具体例子中，坐标旋转公式可选择如下公式：
其中，代表的是旋转任意角度后的新的坐标点；而(x0，y0)代表的是以此点为圆心进行旋转的点，即旋转0点，可称为原点；(x1，y1)代表的是旋转前的任意坐标点；θ即需要旋转的角度(这里是工区方位角与x轴大于0部分的夹角)。
S300：通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息。通过给定的炮检点参数信息如新工区范围中炮线距、检波线距，炮点的点距和检波点间距，结合工区范围信息，计算出炮检点的点数和线数，坐标信息。由于工区已经旋转，该过程计算十分简便。其中，可在工区中设置多条平行的检波线和多条平行的炮点线，检波线和炮点线呈近似于垂直设置，每条检波线上设置有多个检波点，每条炮点线上设置有多个炮点。当中间放炮的模式时，炮点线设置于工区的中央。
S400：通过坐标反旋转公式将所述直角工区炮检点信息转换为实际炮检点信息，根据所述实际炮检点信息进行地震资料分析。得到炮检点信息后，再对所有的点进行反旋转，将各点坐标旋转至真实范围，至此，坐标系统重建过程已完成。
本发明提供了一种简便有效的三维观测系统重建方案，在三维观测系统重建过程中利用已有的坐标旋转公式将带方位角的工区信息旋转成水平正交的信息，免去了角度对后续计算坐标的精度与计算效率影响。通过坐标计算公式计算完所有炮检点信息后再进行角度反旋转，从而得到炮、检点真实坐标信息。在确定任意炮点的接收范围时，将任意炮与给定接收范围的检波点信息进行匹配可快速确定。且通过本发明可根据给定重建相关参数在任意范围内进行三维观测系统信息重新快速准确地建立，重建后的数据信息能为数据规则化、多维数据插值与地震道数据恢复等重要处理技术提供必要的数据支撑，为提高地震资料的成像效果打下基础。
在优选的实施方式中，如图2所示，所述方法进一步包括：
S500：根据给定参数和所述实际炮检点信息建立炮检关系。根据给定参数(每一炮的有多少检波线数接收和每一条检波线有多少道)建立炮检关系，并对该关系进行记录，供后续处理计算需要。
在优选的实施方式中，所述炮检点参数信息包括直角工区中炮线距、检波线距、炮点点距和检波点间距。
所述直角工区炮检点信息包括炮检点的点数、线数和坐标信息。
所述给定参数包括每个炮点对应的检波线数量和每一条检波线的接收通道数量。
在优选的实施方式中，所述通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息具体包括以下步骤的至少之一：
根据直角工区与检波线平行的工区范围信息和检波点距得到每条检波线的检波点点数。具体的，可通过如下公式计算每条检波线的检波点点数：
每条检波线的检波点点数＝(直角工区的最大X坐标-直角工区的原点X坐标)/(检波点间距+1)。
根据直角工区与检波线垂直的工区范围信息和检波点线距得到检波点线数。具体的，可通过如下公式计算检波点线数：
检波点线数＝(直角工区的最大Y坐标-直角工区原点Y坐标)/(检波点线距+1)。
在优选的实施方式中，所述炮点位于所述工区的中央。炮点位于所述工区的中央即为中间放炮的模式。
在优选的实施方式中，在中间放炮模式时，所述通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息具体还包括以下步骤的至少之一：
根据直角工区与检波线平行的工区范围信息、检波点距和炮线距得到炮点线数。具体的，可通过如下公式计算炮点线数：
炮点线数＝(直角工区的最大X坐标-直角工区的原点X坐标-检波点间距)/(炮线间距+1)。
根据检波线数量、检波点线距和炮点距得到每条炮线炮点点数。具体的，可通过如下公式计算每条炮线炮点点数：
每条炮线炮点点数＝(检波线数-1)×检波点线距/炮点距。
根据炮点线数和每条炮线炮点点数得到炮点总数。具体的，可通过如下公式计算炮点总数：
炮点总数＝炮点线数×一条炮线炮点点数。
根据每条检波线的检波点点数和检波点线数得到检波点总数。具体的，可通过如下公式计算检波点总数：
检波点总数＝每条检波线的检波点点数×检波点线数。
根据直角工区与检波线平行的工区原点坐标、一个炮点的线号和线号距得到所述一个炮点与检波线平行的直角工区坐标。具体的，可通过如下公式计算一个炮点与检波线平行的直角工区坐标：
一个炮点与检波线平行的直角工区坐标X＝原点X坐标+该炮点的线号×线号间距。
根据直角工区与检波线垂直的工区原点坐标、一个炮点的点号和点号距得到所述一个炮点与检波线垂直的直角工区坐标。具体的，可通过如下公式计算一个炮点与检波线垂直的直角工区坐标：
一个炮点与检波线垂直的直角工区坐标Y＝原点Y坐标+该炮点的点号×点号间距-(点号间距×0.5)。
根据直角工区与检波线平行的工区原点坐标、一个检波点的点号和点号距得到所述一个检波点与检波线平行的直角工区坐标。具体的，可通过如下公式计算一个检波点与检波线平行的直角工区坐标：
一个检波点与检波线平行的直角工区坐标X＝原点X坐标+该检波点的点号×点号间距。
根据直角工区与检波线垂直的工区原点坐标、一个检波点的线号和线号距得到所述一个检波点与检波线垂直的直角工区坐标。具体的，可通过如下公式计算一个检波点与检波线垂直的直角工区坐标：
一个检波点与检波线平行的直角工区坐标Y＝原点Y坐标+该检波点的线号×线号间距。
其中，各点的点号和各线的线号可按照坐标系的方向按炮点、检波点、炮点线和检波线的排列顺序对炮点、检波点、炮点线和检波线进行编号，以准确表示炮点、检波点、炮点线和检波线间的相对位置。
在一个具体例子中，如图3所示，首先根据数据进行工区范围的确定，根据工区的坐标可以找出该工区的四点坐标。进一步可确定出该工区的方位角，为坐标旋转做好数据准备。而按照坐标旋转公式旋转后，得到的直角工区如下图4所示，范围方正，但中间点位分布不规则，需要重建规则的炮检点坐标。按照前述炮检点坐标公式，在旋转后的坐标系下能直接求出各点对应的新坐标如图5所示，放大局部可以检查是否符合中间放炮方式布设如图6所示。最后重建信息正确后，再按照旋转公式进行反旋转，得到与原工区范围一致的新观测系统信息如图7所示，且与原工区布设点信息进行比较，可以直观看出重建的观测系统信息位置规则，弥补了实际布设情况中位置偏离的问题，如图8所示。
进一步地，根据给定参数(每一炮的有多少检波线数接收和每一条检波线有多少道)建立炮检关系，并对该关系进行记录，供后续处理计算需要。要确定任意炮点的接收范围，只需要确定接收检波点范围的最小线号，最大线号，每条线内最小点号和最大点号即可。因此并未采用依次遍历的方式进行，而是通过进行计算对应索引的方式快速确定接收检波点的范围。对于任意旋转后的工区，任意炮点可按照参数信息与重建的观测系统定义出纵偏距和非纵偏距，如图9所示。
假设是10条线接收，则每条线距500米，炮点在中间，则单边非纵偏距＝4.5条×500米＝2250米。
每条线是100道接收，每道25米，炮点在中间，则单边纵偏距＝49.5道×25＝1237.5米。
对于任意炮点的接收范围可由以下公式计算获得：
起始线编号＝(当前炮点y坐标-非纵偏距-原点坐标y)/接收线距。
终止线编号＝(当前炮点y坐标+非纵偏距-原点坐标y)/接收线距。
起始道编号＝(当前炮点x坐标-纵偏距-原点坐标x)/接收道距。
终止道编号＝(当前炮点x坐标+纵偏距-原点坐标x)/接收道距。
综上，即可快速准确的求出任意炮点对应的检波点接收范围，计算结果如图10所示。最终，按照数据格式要求记录下新的炮检点坐标信息与炮检接收关系信息，可为后续处理计算做好数据准备。
基于相同原理，本实施例还公开了一种三维观测系统重建系统。如图11所示，三维观测系统重建系统包括工区确定单元11、工区转换单元12、炮检点确定单元13和炮检点转换单元14。
所述工区确定单元11用于根据工区范围参数确定所述工区的工区大小及工区方位角信息。
优选的，工区范围参数可包括工区的四个角点坐标等信息，其中坐标可为预先设置的直角坐标系，本实施例中，将与检波线平行的方向定义为x方向，与x垂直的方向为y方向。根据工区的四个角点坐标可确定所述工区的工区大小以及工区方位角信息。
所述工区转换单元12用于通过坐标旋转公式和所述工区方位角信息将所述工区的方位角转换为直角得到直角工区。
优选的，可选择旋转围绕的旋转0点(旋转不动点)为四个角点坐标中y坐标最小的点。通过坐标旋转公式和所述工区方位角信息可将所述工区的方位角转换为直角得到直角工区，通过方位角转换可使得后续计算各信息如各炮检点新的坐标不需考虑夹角因素，避免了由此产生的精度误差，同时也提升了计算效率。
在一个具体例子中，坐标旋转公式可选择如下公式：
其中，代表的是旋转任意角度后的新的坐标点；而(x0，y0)代表的是以此点为圆心进行旋转的点，即旋转0点，可称为原点；(x1，y1)代表的是旋转前的任意坐标点；θ即需要旋转的角度(这里是工区方位角与x轴大于0部分的夹角)。
所述炮检点确定单元13用于通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息。
通过给定的炮检点参数信息如新工区范围中炮线距、检波线距，炮点的点距和检波点间距，结合工区范围信息，计算出炮检点的点数和线数，坐标信息。由于工区已经旋转，该过程计算十分简便。其中，可在工区中设置多条平行的检波线和多条平行的炮点线，检波线和炮点线呈近似于垂直设置，每条检波线上设置有多个检波点，每条炮点线上设置有多个炮点。当中间放炮的模式时，炮点线设置于工区的中央。
所述炮检点转换单元14用于通过坐标反旋转公式将所述直角工区炮检点信息转换为实际炮检点信息，根据所述实际炮检点信息进行地震资料分析。得到炮检点信息后，再对所有的点进行反旋转，将各点坐标旋转至真实范围，至此，坐标系统重建过程已完成。
在优选的实施方式中，如图12所示，所述系统进一步包括炮检关系确定单元15。
所述炮检关系确定单元15用于根据给定参数和所述实际炮检点信息建立炮检关系。根据给定参数(每一炮的有多少检波线数接收和每一条检波线有多少道)建立炮检关系，并对该关系进行记录，供后续处理计算需要。
在优选的实施方式中，所述炮检点参数信息包括直角工区中炮线距、检波线距、炮点点距和检波点间距。所述直角工区炮检点信息包括炮检点的点数、线数和坐标信息。所述给定参数包括每个炮点对应的检波线数量和每一条检波线的接收通道数量。
在优选的实施方式中，所述炮检点确定单元13进一步用于实现根据直角工区与检波线平行的工区范围信息和检波点距得到每条检波线的检波点点数或根据直角工区与检波线垂直的工区范围信息和检波点线距得到检波点线数。
具体的，可通过如下公式计算每条检波线的检波点点数：
每条检波线的检波点点数＝(直角工区的最大X坐标-直角工区的原点X坐标)/(检波点间距+1)。
具体的，可通过如下公式计算检波点线数：
检波点线数＝(直角工区的最大Y坐标-直角工区原点Y坐标)/(检波点线距+1)。
在优选的实施方式中，所述炮点位于所述工区的中央。炮点位于所述工区的中央即为中间放炮的模式。
在优选的实施方式中，在中间放炮模式时，所述炮检点确定单元13进一步用于实现根据直角工区与检波线平行的工区范围信息、检波点距和炮线距得到炮点线数；根据检波线数量、检波点线距和炮点距得到每条炮线炮点点数；根据炮点线数和每条炮线炮点点数得到炮点总数、根据每条检波线的检波点点数和检波点线数得到检波点总数；根据直角工区与检波线平行的工区原点坐标、一个炮点的线号和线号距得到所述一个炮点与检波线平行的直角工区坐标；根据直角工区与检波线垂直的工区原点坐标、一个炮点的点号和点号距得到所述一个炮点与检波线垂直的直角工区坐标；根据直角工区与检波线平行的工区原点坐标、一个检波点的点号和点号距得到所述一个检波点与检波线平行的直角工区坐标及根据直角工区与检波线垂直的工区原点坐标、一个检波点的线号和线号距得到所述一个检波点与检波线垂直的直角工区坐标中的至少之一。
具体的，可通过如下公式计算炮点线数：
炮点线数＝(直角工区的最大X坐标-直角工区的原点X坐标-检波点间距)/(炮线间距+1)。
具体的，可通过如下公式计算每条炮线炮点点数：
每条炮线炮点点数＝(检波线数-1)×检波点线距/炮点距。
具体的，可通过如下公式计算炮点总数：
炮点总数＝炮点线数×一条炮线炮点点数。
具体的，可通过如下公式计算检波点总数：
检波点总数＝每条检波线的检波点点数×检波点线数。
具体的，可通过如下公式计算一个炮点与检波线平行的直角工区坐标：
一个炮点与检波线平行的直角工区坐标X＝原点X坐标+该炮点的线号×线号间距。
具体的，可通过如下公式计算一个炮点与检波线垂直的直角工区坐标：
一个炮点与检波线垂直的直角工区坐标Y＝原点Y坐标+该炮点的点号×点号间距-(点号间距×0.5)。
具体的，可通过如下公式计算一个检波点与检波线平行的直角工区坐标：
一个检波点与检波线平行的直角工区坐标X＝原点X坐标+该检波点的点号×点号间距。
具体的，可通过如下公式计算一个检波点与检波线垂直的直角工区坐标：
一个检波点与检波线平行的直角工区坐标Y＝原点Y坐标+该检波点的线号×线号间距。
其中，各点的点号和各线的线号可按照坐标系的方向按炮点、检波点、炮点线和检波线的排列顺序对炮点、检波点、炮点线和检波线进行编号，以准确表示炮点、检波点、炮点线和检波线间的相对位置。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机设备，具体的，计算机设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中计算机设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上所述的由客户端执行的方法，或者，所述处理器执行所述程序时实现如上所述的由服务器执行的方法。
下面参考图13，其示出了适于用来实现本申请实施例的计算机设备600的结构示意图。
如图13所示，计算机设备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。
特别地，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机程序产品，其包括有形地包含在机器可读介质上的计算机程序，所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (11)

1.一种三维观测系统重建方法，其特征在于，包括：
根据工区范围参数确定所述工区的工区大小及工区方位角信息；
通过坐标旋转公式和所述工区方位角信息将所述工区的方位角转换为直角得到直角工区；
通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息；
通过坐标反旋转公式将所述直角工区炮检点信息转换为实际炮检点信息，根据所述实际炮检点信息进行地震资料分析。
2.根据权利要求1所述的三维观测系统重建方法，其特征在于，进一步包括：
根据给定参数和所述实际炮检点信息建立炮检关系。
3.根据权利要求1所述的三维观测系统重建方法，其特征在于，所述炮检点参数信息包括直角工区中炮线距、检波点线距、炮点点距和检波点间距。
4.根据权利要求1所述的三维观测系统重建方法，其特征在于，所述直角工区炮检点信息包括炮检点的点数、线数和坐标信息。
5.根据权利要求2所述的三维观测系统重建方法，其特征在于，所述给定参数包括每个炮点对应的检波线数量和每一条检波线的接收通道数量。
6.根据权利要求1所述的三维观测系统重建方法，其特征在于，所述工区的炮点位于所述工区的中央。
7.根据权利要求1-6任一项所述的三维观测系统重建方法，其特征在于，所述通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息具体包括以下步骤的至少之一：
根据直角工区与检波线平行的工区范围信息和检波点间距得到每条检波线的检波点点数；
根据直角工区与检波线垂直的工区范围信息和检波点线距得到检波点线数。
8.根据权利要求6所述的三维观测系统重建方法，其特征在于，所述通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息具体还包括以下步骤的至少之一：
根据直角工区与检波线平行的工区范围信息、检波点间距和炮线距得到炮点线数；
根据检波线数量、检波点线距和炮点点距得到每条炮线炮点点数；
根据炮点线数和每条炮线炮点点数得到炮点总数；
根据每条检波线的检波点点数和检波点线数得到检波点总数；
根据直角工区与检波线平行的工区原点坐标、一个炮点的线号和线号距得到所述一个炮点与检波线平行的直角工区坐标；
根据直角工区与检波线垂直的工区原点坐标、一个炮点的点号和点号距得到所述一个炮点与检波线垂直的直角工区坐标；
根据直角工区与检波线平行的工区原点坐标、一个检波点的点号和点号距得到所述一个检波点与检波线平行的直角工区坐标；
根据直角工区与检波线垂直的工区原点坐标、一个检波点的线号和线号距得到所述一个检波点与检波线垂直的直角工区坐标。
9.一种三维观测系统重建系统，其特征在于，包括：
工区确定单元，用于根据工区范围参数确定所述工区的工区大小及工区方位角信息；
工区转换单元，用于通过坐标旋转公式和所述工区方位角信息将所述工区的方位角转换为直角得到直角工区；
炮检点确定单元，用于通过炮检点参数信息和所述直角工区的工区范围信息得到直角工区炮检点信息；
炮检点转换单元，用于通过坐标反旋转公式将所述直角工区炮检点信息转换为实际炮检点信息，根据所述实际炮检点信息进行地震资料分析。
10.一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，
所述处理器执行所述程序时实现如权利要求1-8任一项所述方法。
11.一种计算机可读介质，其上存储有计算机程序，其特征在于，
该程序被处理器执行时实现如权利要求1-8任一项所述方法。
CN201911171968.XA 2019-11-26 2019-11-26 一种三维观测系统重建方法及系统 Active CN112946736B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911171968.XA CN112946736B (zh) 2019-11-26 2019-11-26 一种三维观测系统重建方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911171968.XA CN112946736B (zh) 2019-11-26 2019-11-26 一种三维观测系统重建方法及系统

Publications (2)

Publication Number Publication Date
CN112946736A CN112946736A (zh) 2021-06-11
CN112946736B true CN112946736B (zh) 2024-06-25

Family

ID=76224787

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911171968.XA Active CN112946736B (zh) 2019-11-26 2019-11-26 一种三维观测系统重建方法及系统

Country Status (1)

Country Link
CN (1) CN112946736B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114488287B (zh) * 2022-01-25 2024-05-28 中国海洋大学 一种由地震数据道头炮检点坐标重建sps文件的方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502967A (zh) * 2014-12-23 2015-04-08 中国石油天然气集团公司 快速获取地震勘探观测系统面元信息的方法及装置

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2265235C1 (ru) * 2004-07-29 2005-11-27 Закрытое Акционерное Общество "Геоспектр" Способ поисков и разведки нефтегазовых месторождений
JP2007322217A (ja) * 2006-05-31 2007-12-13 Yamatake Corp 地震観測システム
CN102053277A (zh) * 2009-10-30 2011-05-11 中国石油化工股份有限公司 一种利用地震资料进行检测储层裂缝发育方向的方法
CN103217711B (zh) * 2013-03-14 2015-07-08 中国石油天然气股份有限公司 获取共偏移距共方位角域道集的方法及系统
RU2563323C1 (ru) * 2014-04-02 2015-09-20 Открытое акционерное общество "Нефтяная компания "Роснефть" Способ реконструкции тонкой структуры геологического объекта и прогноза его флюидонасыщения
CN106249283A (zh) * 2015-06-15 2016-12-21 中国石油化工股份有限公司 多元化的地震资料观测系统定义方法及系统
CN106772583B (zh) * 2017-01-10 2018-09-04 中国科学院地质与地球物理研究所 一种地震绕射波分离方法和装置
CN107561574A (zh) * 2017-07-20 2018-01-09 中国石油天然气集团公司 炮点位置、接收点位置的检测方法和装置
CN109407143B (zh) * 2017-08-16 2020-06-30 中石化石油工程技术服务有限公司 基于压缩感知的地震勘探非规则观测系统设计方法
CN107678063B (zh) * 2017-09-25 2019-11-29 中国石油集团东方地球物理勘探有限责任公司 一种基于等级相关分析的多分量转换波裂缝预测方法
CN108107471B (zh) * 2017-11-17 2019-10-11 中国石油天然气集团公司 一种分方位初至数据体的获取方法及装置
CN108363097B (zh) * 2018-02-02 2019-10-15 中国石油大学(华东) 一种地震资料偏移成像方法
CN108761520B (zh) * 2018-05-08 2019-12-10 中石化石油工程技术服务有限公司 一种地震采集观测系统参数设计的方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502967A (zh) * 2014-12-23 2015-04-08 中国石油天然气集团公司 快速获取地震勘探观测系统面元信息的方法及装置

Also Published As

Publication number Publication date
CN112946736A (zh) 2021-06-11

Similar Documents

Publication Publication Date Title
CN110285754B (zh) 基于激光扫描的工件定位方法、装置、系统和存储介质
CN107820314B (zh) 基于AP选择的Dwknn位置指纹定位方法
CN112946736B (zh) 一种三维观测系统重建方法及系统
CN117664087B (zh) 垂轨环扫式卫星影像核线生成方法、系统及设备
CN110427506A (zh) 空间数据边界处理方法、装置、计算机设备及存储介质
CN113534259A (zh) 一种可控震源高效采集实时叠前时间偏移成像方法
CN111474574A (zh) 基于压缩感知的地震采集观测系统的生成方法及装置
CN112379436A (zh) 混合震源激发地震资料静校正量计算方法及装置
CN112764102B (zh) 地震资料采集方法及装置
CN111538004B (zh) 圆周扫描地基sar的多相位中心分布设计方法及装置
CN108828655B (zh) 道集记录与叠加剖面的处理方法、装置及计算机存储介质
CN113568033A (zh) 一种三维非规则采样地震采集观测系统设计方法和装置
CN114002742B (zh) 一种欧拉高斯束偏移成像方法及装置
CN111337973B (zh) 地震数据重建方法、系统
CN113805233B (zh) 一种点扩散函数的计算方法
CN112255679A (zh) 地震资料绕射深度偏移处理方法及装置
CN112444876A (zh) 一种地震数据处理方法及装置
CN113240290A (zh) 一种观测弧段规划生成方法、系统、存储介质及电子设备
CN104880732B (zh) 一种十字子集的构建方法及装置
Guan et al. Parabolic equation modeling of propagation over terrain using digital elevation model
CN111694911B (zh) 一种在指定范围内生成随机分布点的方法、设备和系统
CN115511925B (zh) 基于管道通信机制的StaMPS并行优化处理方法及产品
CN113568047B (zh) 一种叠前成像道集生成方法和装置
CN112241023B (zh) 陆上三维炮点二次定位及定量偏移监控方法和装置
CN112379410A (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