CN104502967A - 快速获取地震勘探观测系统面元信息的方法及装置 - Google Patents
快速获取地震勘探观测系统面元信息的方法及装置 Download PDFInfo
- Publication number
- CN104502967A CN104502967A CN201410814786.0A CN201410814786A CN104502967A CN 104502967 A CN104502967 A CN 104502967A CN 201410814786 A CN201410814786 A CN 201410814786A CN 104502967 A CN104502967 A CN 104502967A
- Authority
- CN
- China
- Prior art keywords
- bin
- logic
- work area
- cautious
- maxhalfoffset
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种快速获取地震勘探观测系统面元信息的方法及装置，该方法包括以下步骤：构建工区内待计算面元的逻辑炮检点矩形网格；基于所述逻辑炮检点矩形网格建立所述工区内待计算面元的面元矩阵；对所述面元矩阵进行并行计算实现快速获取所述旋转后工区的面元信息。本发明大大提高了获取面元信息效率。
Description
技术领域
本发明涉及地球物理勘探技术领域，尤其是涉及一种快速获取地震勘探观测系统面元信息的方法及装置。
背景技术
地震勘探设计是地震勘探过程中的首要任务之一，在设计过程中需要对地震勘探观测系统进行充分地论证，用面元的覆盖次数、炮检对中点坐标、炮检对距离以及炮检连线方位角等数据信息进行全面的分析论证。面元信息计算就是找出地震勘探工区中每个面元的这些属性值，供技术人员分析使用，因此，面元信息是石油勘探中观测系统优化设计的重要数据信息。
目前石油勘探已经进入了精细勘探时代，三维和高分辨勘探已成为主要的勘探施工方案，观测系统的优化设计中就需要计算大量的各种类型的面元信息数据。这些观测系统一般都有上亿个炮检对，计算出的面元信息达到几十GB乃至几百GB的数据量。而采用传统的面元信息计算方法是：如计算第(i,j)面元的信息，则从炮点开始循环搜索，依据炮检关系，循环搜索检点，测试该炮检对的中点是否在该面元内，若是则计算其他信息，否则继续下个点。这种计算方法耗费大量的计算时间用于查找符合条件的炮检对，从而给后续面元分析工作带来了瓶颈问题，对工作效益造成巨大影响。
发明内容
本发明的目的在于提供一种快速获取地震勘探观测系统面元信息的方法及装置，以提高面元信息的获取效率和灵活性。
为达到上述目的，一方面，本发明提供了一种快速获取地震勘探观测系统面元信息的方法，包括以下步骤：
构建工区内待计算面元的逻辑炮检点矩形网格；
基于所述逻辑炮检点矩形网格建立所述工区内待计算面元的面元矩阵；
对所述面元矩阵进并行计算获取所述旋转后工区的面元信息。
本发明的获取地震勘探观测系统面元信息的方法，所述构建工区内待计算面元的逻辑炮检矩形网格，具体包括：
将待计算检点面元的工区按预设角度旋转，使旋转后工区内的检点的主测线方向平行于坐标轴，获得旋转后工区内的逻辑检点；
平移所述逻辑检点的分布，使所述逻辑检点的检线距及检点距均匀分布；
依据平移后的逻辑检点的分布，通过添加虚拟检点的方式使每条检线上的逻辑检点相同，获取逻辑检点数据结构；
依据上述方法获取逻辑炮点数据结构；
将逻辑检点坐标及逻辑炮点坐标对应读入所述逻辑检点数据结构及所述逻辑炮点数据结构，并依据所述逻辑检点的主测线方向重新调整所述逻辑炮点分布，使得调整后的逻辑炮点的炮线方向与所述逻辑检点的检线相同，从而获取工区内待计算面元的逻辑炮检矩形网格。
本发明的获取地震勘探观测系统面元信息的方法，所述面元矩阵为[nx,ny]，其中，nx表示X方向的面元数，ny表示Y方向的面元数，X方向的面元大小表示为dfx，Y方向的面元大小表示为dfy。
本发明的获取地震勘探观测系统面元信息的方法，所述对所述面元矩阵进并行计算获取所述旋转后工区的面元信息，具体包括：
根据公式isl0＝(i_fex-maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的首炮线isl0；
根据公式isl1＝(i_fex+maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的终炮线isl1；
根据公式isp0＝(j_fey-maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的首炮点isp0；
根据公式isp1＝(j_fey+maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的终炮点isp1；
根据公式irl0＝(i_fex-maxHalfOffset)/drl计算所述旋转后工区的第i行第j列面元相关的首检线irl0；
根据公式irl1＝(i_fex+maxHalfOffset)/drl计算所述旋转后工区的面元的终检线irl1；
根据公式irp0＝(j_fey-maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的首检点irp0；
根据公式irp1＝(j_fey+maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的终检点irp1；
其中，dsl、dsp分别表示炮线间隔、炮点间隔，drl、drp分别表示检线间隔、检点间隔，maxHalfOffset表示最大炮检距的一半，i_fex、j_fey分别表示第i行第j列面元的横坐标和纵坐标。
本发明的获取地震勘探观测系统面元信息的方法，所述在对所述面元矩阵进并行计算获取所述旋转后工区的面元信息之后，还包括：
剔除虚拟检点面元信息和虚拟炮点面元信息，获得有效炮检点面元信息。
另一方面，本发明还提供了一种获取地震勘探观测系统面元信息的装置，包括：
炮检点矩形网格构建模块，用于构建工区内待计算面元的逻辑炮检点矩形网格；
面元矩阵建立模块，用于基于所述逻辑炮检点矩形网格建立所述工区内待计算面元的面元矩阵；
面元信息计算模块，用于对所述面元矩阵进并行计算获取所述旋转后工区的面元信息。
本发明的获取地震勘探观测系统面元信息的装置，所述构建工区内待计算面元的逻辑炮检矩形网格，具体包括：
将待计算检点面元的工区按预设角度旋转，使旋转后工区内的检点的主测线方向平行于坐标轴，获得旋转后工区内的逻辑检点；
平移所述逻辑检点的分布，使所述逻辑检点的检线距及检点距均匀分布；
依据平移后的逻辑检点的分布，通过添加虚拟检点的方式使每条检线上的逻辑检点相同，获取逻辑检点数据结构；
依据上述方法获取逻辑炮点数据结构；
将逻辑检点坐标及逻辑炮点坐标对应读入所述逻辑检点数据结构及所述逻辑炮点数据结构，并依据所述逻辑检点的主测线方向重新调整所述逻辑炮点分布，使得调整后的逻辑炮点的炮线方向与所述逻辑检点的检线相同，从而获取工区内待计算面元的逻辑炮检矩形网格。
本发明的获取地震勘探观测系统面元信息的装置，所述面元矩阵为[nx,ny]，其中，nx表示X方向的面元数，ny表示Y方向的面元数，X方向的面元大小表示为dfx，Y方向的面元大小表示为dfy。
本发明的获取地震勘探观测系统面元信息的装置，所述对所述面元矩阵进并行计算获取所述旋转后工区的面元信息，具体包括：
根据公式isl0＝(i_fex-maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的首炮线isl0；
根据公式isl1＝(i_fex+maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的终炮线isl1；
根据公式isp0＝(j_fey-maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的首炮点isp0；
根据公式isp1＝(j_fey+maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的终炮点isp1；
根据公式irl0＝(i_fex-maxHalfOffset)/drl计算所述旋转后工区的第i行第j列面元相关的首检线irl0；
根据公式irl1＝(i_fex+maxHalfOffset)/drl计算所述旋转后工区的面元的终检线irl1；
根据公式irp0＝(j_fey-maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的首检点irp0；
根据公式irp1＝(j_fey+maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的终检点irp1；
其中，dsl、dsp分别表示炮线间隔、炮点间隔，drl、drp分别表示检线间隔、检点间隔，maxHalfOffset表示最大炮检距的一半，i_fex、j_fey分别表示第i行第j列面元的横坐标和纵坐标。
本发明的获取地震勘探观测系统面元信息的装置，所述在对所述面元矩阵进并行计算获取所述旋转后工区的面元信息之后，还包括：
剔除虚拟检点面元信息和虚拟炮点面元信息，获得有效炮检点面元信息。
本发明首先构建工区内待计算面元的逻辑炮检点矩形网格，然后基于逻辑炮检点矩形网格建立工区内待计算面元的面元矩阵以利于后续并行处理，最后对面元矩阵进并行计算获取旋转后工区的面元信息，从而大大提高了获取面元信息效率。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1为本发明实施例的快速获取地震勘探观测系统面元信息的方法流程图；
图2为本发明实施例中一个工区的地震勘探观测系统的布设图；
图3为图2所示地震勘探观测系统经过旋转平移修正后的逻辑地震勘探观测系统的布设图；
图4为本发明实施例的获取地震勘探观测系统面元信息的装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
下面结合附图，对本发明的具体实施方式作进一步的详细说明。
参考图1所示，本发明实施例的快速获取地震勘探观测系统面元信息的方法包括以下步骤：
步骤S1、构建工区内待计算面元的逻辑炮检点矩形网格。具体步骤如下：
将待计算检点面元的工区(如图2所示)按预设角度旋转，使旋转后工区内的检点的主测线方向平行于坐标轴(例如X轴或Y轴)，获得旋转后工区内的逻辑检点；
平移所述逻辑检点的分布，使所述逻辑检点的检线距及检点距均匀分布；
依据平移后的逻辑检点的分布，通过添加虚拟检点的方式使每条检线上的逻辑检点相同，获取逻辑检点数据结构；
同理依据上述方法可获取逻辑炮点数据结构；
将逻辑检点坐标及逻辑炮点坐标对应读入所述逻辑检点数据结构及所述逻辑炮点数据结构，并依据所述逻辑检点的主测线方向重新调整所述逻辑炮点分布，使得调整后的逻辑炮点的炮线方向与所述逻辑检点的检线相同，从而获取工区内待计算面元的逻辑炮检矩形网格(如图3所示)，其中，用网格编号表示实际面元。因此本步骤能够得到：检线距drl，检点距drp；炮线距dsl，炮点距dsp；工区内的检点范围：minrx,minry,maxrx,maxry；炮点范围：minsx,minsy,maxsx,maxsy等。
步骤S2、基于所述逻辑炮检点矩形网格建立所述工区内待计算面元的面元矩阵。假设面元矩阵为[nx,ny]，其中，nx表示X方向的面元数，ny表示Y方向的面元数，X方向的面元大小表示为dfx，Y方向的面元大小表示为dfy。
步骤S3、对所述面元矩阵进并行计算获取所述旋转后工区的面元信息。具体的依据面元矩阵中各个面元的位置(i,j)计算它的覆盖次数、炮检对的中点、炮检方位角、炮检距等等面元属性数据。其中，一些具体计算可参见如下内容：
根据公式isl0＝(i_fex-maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的首炮线isl0；
根据公式isl1＝(i_fex+maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的终炮线isl1；
根据公式isp0＝(j_fey-maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的首炮点isp0；
根据公式isp1＝(j_fey+maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的终炮点isp1；
根据公式irl0＝(i_fex-maxHalfOffset)/drl计算所述旋转后工区的第i行第j列面元相关的首检线irl0；
根据公式irl1＝(i_fex+maxHalfOffset)/drl计算所述旋转后工区的面元的终检线irl1；
根据公式irp0＝(j_fey-maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的首检点irp0；
根据公式irp1＝(j_fey+maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的终检点irp1；
其中，dsl、dsp分别表示炮线间隔、炮点间隔，drl、drp分别表示检线间隔、检点间隔，maxHalfOffset表示最大炮检距的一半，i_fex、j_fey分别表示第i行第j列面元的横坐标和纵坐标。
本发明实施例中，在对所述面元矩阵进并行计算获取所述旋转后工区的面元信息之后，还可以依据地震观测系统的定义规则以及面元的定义规则，剔除虚拟检点面元信息和虚拟炮点面元信息，获得有效炮检点面元信息。
本发明实施例首先构建工区内待计算面元的逻辑炮检点矩形网格，然后基于逻辑炮检点矩形网格建立工区内待计算面元的面元矩阵以利于后续并行处理，最后对面元矩阵进并行计算获取旋转后工区的面元信息，从而大大提高了获取面元信息效率。
结合图4所示，与本发明实施例的获取地震勘探观测系统面元信息的方法对应，本发明实施例的获取地震勘探观测系统面元信息的装置包括炮检点矩形网格构建模块41、面元矩阵建立模块42和面元信息计算模块43，其中：
炮检点矩形网格构建模块41，用于构建工区内待计算面元的逻辑炮检点矩形网格。具体的：
将待计算检点面元的工区(如图2所示)按预设角度旋转，使旋转后工区内的检点的主测线方向平行于坐标轴(例如X轴或Y轴)，获得旋转后工区内的逻辑检点；
平移所述逻辑检点的分布，使所述逻辑检点的检线距及检点距均匀分布；
依据平移后的逻辑检点的分布，通过添加虚拟检点的方式使每条检线上的逻辑检点相同，获取逻辑检点数据结构；
同理依据上述方法可获取逻辑炮点数据结构；
将逻辑检点坐标及逻辑炮点坐标对应读入所述逻辑检点数据结构及所述逻辑炮点数据结构，并依据所述逻辑检点的主测线方向重新调整所述逻辑炮点分布，使得调整后的逻辑炮点的炮线方向与所述逻辑检点的检线相同，从而获取工区内待计算面元的逻辑炮检矩形网格(如图3所示)，其中，用网格编号表示实际面元。因此本步骤能够得到：检线距drl，检点距drp；炮线距dsl，炮点距dsp；工区内的检点范围：minrx,minry,maxrx,maxry；炮点范围：minsx,minsy,maxsx,maxsy等。
面元矩阵建立模块42，用于基于所述逻辑炮检点矩形网格建立所述工区内待计算面元的面元矩阵。假设面元矩阵为[nx,ny]，其中，nx表示X方向的面元数，ny表示Y方向的面元数，X方向的面元大小表示为dfx，Y方向的面元大小表示为dfy。
面元信息计算模块43，用于对所述面元矩阵进并行计算获取所述旋转后工区的面元信息。具体的依据面元矩阵中各个面元的位置(i,j)计算它的覆盖次数、炮检对的中点、炮检方位角、炮检距等等面元属性数据。其中，一些具体计算可参见如下内容：
根据公式isl0＝(i_fex-maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的首炮线isl0；
根据公式isl1＝(i_fex+maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的终炮线isl1；
根据公式isp0＝(j_fey-maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的首炮点isp0；
根据公式isp1＝(j_fey+maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的终炮点isp1；
根据公式irl0＝(i_fex-maxHalfOffset)/drl计算所述旋转后工区的第i行第j列面元相关的首检线irl0；
根据公式irl1＝(i_fex+maxHalfOffset)/drl计算所述旋转后工区的面元的终检线irl1；
根据公式irp0＝(j_fey-maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的首检点irp0；
根据公式irp1＝(j_fey+maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的终检点irp1；
其中，dsl、dsp分别表示炮线间隔、炮点间隔，drl、drp分别表示检线间隔、检点间隔，maxHalfOffset表示最大炮检距的一半，i_fex、j_fey分别表示第i行第j列面元的横坐标和纵坐标。
本发明实施例首先构建工区内待计算面元的逻辑炮检点矩形网格，然后基于逻辑炮检点矩形网格建立工区内待计算面元的面元矩阵以利于后续并行处理，最后对面元矩阵进并行计算获取旋转后工区的面元信息，从而大大提高了获取面元信息效率。
本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块、单元和步骤可以通过硬件、软件或两者的结合来实现。至于是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本发明实施例保护的范围。
本发明实施例中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器，数字信号处理器，专用集成电路(ASIC)，现场可编程门阵列或其它可编程逻辑装置，离散门或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置来实现。
本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中，本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种快速获取地震勘探观测系统面元信息的方法，其特征在于，包括以下步骤：
构建工区内待计算面元的逻辑炮检点矩形网格；
基于所述逻辑炮检点矩形网格建立所述工区内待计算面元的面元矩阵；
对所述面元矩阵进行并行计算实现快速获取所述旋转后工区的面元信息。
2.根据权利要求1所述的快速获取地震勘探观测系统面元信息的方法，其特征在于，所述构建工区内待计算面元的逻辑炮检矩形网格，具体包括：
将待计算检点面元的工区按预设角度旋转，使旋转后工区内的检点的主测线方向平行于坐标轴，获得旋转后工区内的逻辑检点；
平移所述逻辑检点的分布，使所述逻辑检点的检线距及检点距均匀分布；
依据平移后的逻辑检点的分布，通过添加虚拟检点的方式使每条检线上的逻辑检点相同，获取逻辑检点数据结构；
依据上述方法获取逻辑炮点数据结构；
将逻辑检点坐标及逻辑炮点坐标对应读入所述逻辑检点数据结构及所述逻辑炮点数据结构，并依据所述逻辑检点的主测线方向重新调整所述逻辑炮点分布，使得调整后的逻辑炮点的炮线方向与所述逻辑检点的检线相同，从而获取工区内待计算面元的逻辑炮检矩形网格。
3.根据权利要求1所述的快速获取地震勘探观测系统面元信息的方法，其特征在于，所述面元矩阵为[nx,ny]，其中，nx表示X方向的面元数，ny表示Y方向的面元数，X方向的面元大小表示为dfx，Y方向的面元大小表示为dfy。
4.根据权利要求1所述的快速获取地震勘探观测系统面元信息的方法，其特征在于，所述对所述面元矩阵进行并行计算实现快速获取所述旋转后工区的面元信息，具体包括：
根据公式isl0＝(i_fex-maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的首炮线isl0；
根据公式isl1＝(i_fex+maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的终炮线isl1；
根据公式isp0＝(j_fey-maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的首炮点isp0；
根据公式isp1＝(j_fey+maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的终炮点isp1；
根据公式irl0＝(i_fex-maxHalfOffset)/drl计算所述旋转后工区的第i行第j列面元相关的首检线irl0；
根据公式irl1＝(i_fex+maxHalfOffset)/drl计算所述旋转后工区的面元的终检线irl1；
根据公式irp0＝(j_fey-maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的首检点irp0；
根据公式irp1＝(j_fey+maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的终检点irp1；
其中，dsl、dsp分别表示炮线间隔、炮点间隔，drl、drp分别表示检线间隔、检点间隔，maxHalfOffset表示最大炮检距的一半，i_fex、j_fey分别表示第i行第j列面元的横坐标和纵坐标。
5.根据权利要求1所述的快速获取地震勘探观测系统面元信息的方法，其特征在于，所述在对所述面元矩阵进行并行计算实现快速获取所述旋转后工区的面元信息之后，还包括：
剔除虚拟检点面元信息和虚拟炮点面元信息，获得有效炮检点面元信息。
6.一种快速获取地震勘探观测系统面元信息的装置，其特征在于，包括：
炮检点矩形网格构建模块，用于构建工区内待计算面元的逻辑炮检点矩形网格；
面元矩阵建立模块，用于基于所述逻辑炮检点矩形网格建立所述工区内待计算面元的面元矩阵；
面元信息计算模块，用于对所述面元矩阵进行并行计算实现快速获取所述旋转后工区的面元信息。
7.根据权利要求6所述的快速获取地震勘探观测系统面元信息的装置，其特征在于，所述构建工区内待计算面元的逻辑炮检矩形网格，具体包括：
将待计算检点面元的工区按预设角度旋转，使旋转后工区内的检点的主测线方向平行于坐标轴，获得旋转后工区内的逻辑检点；
平移所述逻辑检点的分布，使所述逻辑检点的检线距及检点距均匀分布；
依据平移后的逻辑检点的分布，通过添加虚拟检点的方式使每条检线上的逻辑检点相同，获取逻辑检点数据结构；
依据上述方法获取逻辑炮点数据结构；
将逻辑检点坐标及逻辑炮点坐标对应读入所述逻辑检点数据结构及所述逻辑炮点数据结构，并依据所述逻辑检点的主测线方向重新调整所述逻辑炮点分布，使得调整后的逻辑炮点的炮线方向与所述逻辑检点的检线相同，从而获取工区内待计算面元的逻辑炮检矩形网格。
8.根据权利要求6所述的快速获取地震勘探观测系统面元信息的装置，其特征在于，所述面元矩阵为[nx,ny]，其中，nx表示X方向的面元数，ny表示Y方向的面元数，X方向的面元大小表示为dfx，Y方向的面元大小表示为dfy。
9.根据权利要求6所述的快速获取地震勘探观测系统面元信息的装置，其特征在于，所述对所述面元矩阵进并行计算获取所述旋转后工区的面元信息，具体包括：
根据公式isl0＝(i_fex-maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的首炮线isl0；
根据公式isl1＝(i_fex+maxHalfOffset)/dsl计算所述旋转后工区的第i行第j列面元相关的终炮线isl1；
根据公式isp0＝(j_fey-maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的首炮点isp0；
根据公式isp1＝(j_fey+maxHalfOffset)/dsp计算所述旋转后工区的第i行第j列面元相关的终炮点isp1；
根据公式irl0＝(i_fex-maxHalfOffset)/drl计算所述旋转后工区的第i行第j列面元相关的首检线irl0；
根据公式irl1＝(i_fex+maxHalfOffset)/drl计算所述旋转后工区的面元的终检线irl1；
根据公式irp0＝(j_fey-maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的首检点irp0；
根据公式irp1＝(j_fey+maxHalfOffset)/drp计算所述旋转后工区的第i行第j列面元相关的终检点irp1；
其中，dsl、dsp分别表示炮线间隔、炮点间隔，drl、drp分别表示检线间隔、检点间隔，maxHalfOffset表示最大炮检距的一半，i_fex、j_fey分别表示第i行第j列面元的横坐标和纵坐标。
10.根据权利要求6所述的快速获取地震勘探观测系统面元信息的装置，其特征在于，所述在对所述面元矩阵进并行计算获取所述旋转后工区的面元信息之后，还包括：
剔除虚拟检点面元信息和虚拟炮点面元信息，获得有效炮检点面元信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410814786.0A CN104502967B (zh) | 2014-12-23 | 2014-12-23 | 快速获取地震勘探观测系统面元信息的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410814786.0A CN104502967B (zh) | 2014-12-23 | 2014-12-23 | 快速获取地震勘探观测系统面元信息的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104502967A true CN104502967A (zh) | 2015-04-08 |
| CN104502967B CN104502967B (zh) | 2017-03-08 |
Family
ID=52944378
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410814786.0A Active CN104502967B (zh) | 2014-12-23 | 2014-12-23 | 快速获取地震勘探观测系统面元信息的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104502967B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105738949A (zh) * | 2016-03-01 | 2016-07-06 | 中国海洋石油总公司 | 一种用于时移地震的九面元一致性并行处理方法 |
| CN110068863A (zh) * | 2019-04-04 | 2019-07-30 | 中国科学院计算机网络信息中心 | 一种共深度点cdp的确定方法及装置 |
| CN112305599A (zh) * | 2019-07-30 | 2021-02-02 | 中国石油天然气集团有限公司 | 高密度宽方位的地震数据分类方法及装置 |
| CN112748458A (zh) * | 2019-10-30 | 2021-05-04 | 中国石油天然气集团有限公司 | 用于现场精细变观的在线实时面元获取方法及装置 |
| CN112946736A (zh) * | 2019-11-26 | 2021-06-11 | 中国石油天然气集团有限公司 | 一种三维观测系统重建方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090147620A1 (en) * | 2007-12-11 | 2009-06-11 | Naide Pan | Method for processing marine towed streamer seismic data from regular multi-azimuth surveys |
| CN102455435A (zh) * | 2010-10-22 | 2012-05-16 | 中国石油集团东方地球物理勘探有限责任公司 | 块状模板观测系统满覆盖自动布设方法 |
| CN103675893A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 海上拖缆观测系统复杂模板满覆盖布设方法 |
-
2014
- 2014-12-23 CN CN201410814786.0A patent/CN104502967B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090147620A1 (en) * | 2007-12-11 | 2009-06-11 | Naide Pan | Method for processing marine towed streamer seismic data from regular multi-azimuth surveys |
| CN102455435A (zh) * | 2010-10-22 | 2012-05-16 | 中国石油集团东方地球物理勘探有限责任公司 | 块状模板观测系统满覆盖自动布设方法 |
| CN103675893A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 海上拖缆观测系统复杂模板满覆盖布设方法 |
Non-Patent Citations (1)
| Title |
|---|
| 李振勇 等: "海洋地震数据处理技术探讨", 《石油地球物理勘探》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105738949A (zh) * | 2016-03-01 | 2016-07-06 | 中国海洋石油总公司 | 一种用于时移地震的九面元一致性并行处理方法 |
| CN110068863A (zh) * | 2019-04-04 | 2019-07-30 | 中国科学院计算机网络信息中心 | 一种共深度点cdp的确定方法及装置 |
| CN112305599A (zh) * | 2019-07-30 | 2021-02-02 | 中国石油天然气集团有限公司 | 高密度宽方位的地震数据分类方法及装置 |
| CN112748458A (zh) * | 2019-10-30 | 2021-05-04 | 中国石油天然气集团有限公司 | 用于现场精细变观的在线实时面元获取方法及装置 |
| CN112748458B (zh) * | 2019-10-30 | 2024-05-28 | 中国石油天然气集团有限公司 | 用于现场精细变观的在线实时面元获取方法及装置 |
| CN112946736A (zh) * | 2019-11-26 | 2021-06-11 | 中国石油天然气集团有限公司 | 一种三维观测系统重建方法及系统 |
| CN112946736B (zh) * | 2019-11-26 | 2024-06-25 | 中国石油天然气集团有限公司 | 一种三维观测系统重建方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104502967B (zh) | 2017-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Park et al. | Parallel computation of skyline and reverse skyline queries using mapreduce | |
| JP5425541B2 (ja) | マルチプロセッサ・システム上でデータ・セットを区分化およびソートするための方法および装置 | |
| US20060139349A1 (en) | Applications of interval arithmetic for reduction of number of computations in ray tracing problems | |
| CN104502967A (zh) | 快速获取地震勘探观测系统面元信息的方法及装置 | |
| CN104346384A (zh) | 一种小文件处理方法及装置 | |
| CN102591855A (zh) | 一种数据标识方法及系统 | |
| CN105556474A (zh) | 管理数据操作的存储器和存储空间 | |
| CN103455531A (zh) | 一种支持高维数据实时有偏查询的并行索引方法 | |
| CN114283099A (zh) | 一种图处理的方法，系统以及装置 | |
| CN108646297A (zh) | 海洋气枪震源远场子波模拟方法及装置 | |
| US20240078255A1 (en) | Method and apparatus for determining spatial two-tuple, computer device, and storage medium | |
| Nikitin et al. | Cache-efficient parallel eikonal solver for multicore CPUs | |
| CN105359142B (zh) | 哈希连接方法和装置 | |
| CN109657197A (zh) | 一种叠前深度偏移计算方法及系统 | |
| CN104572862A (zh) | 一种海量数据存储访问方法及系统 | |
| CN104182334A (zh) | Numa系统的内存测试方法及系统 | |
| CN104570077A (zh) | 一种逆时偏移偏移距域共成像点道集提取方法 | |
| CN109298450A (zh) | 地震采集观测系统属性分析方法及装置 | |
| Li et al. | Losha: A general framework for scalable locality sensitive hashing | |
| CN105700998A (zh) | 一种对并行程序的性能进行监测分析的方法及装置 | |
| CN106484532B (zh) | 面向sph流体模拟的gpgpu并行计算方法 | |
| Beier et al. | GiST scan acceleration using coprocessors | |
| CN106407233A (zh) | 一种数据处理方法和设备 | |
| US11507290B2 (en) | Method and apparatus for storage device management | |
| CN112597662A (zh) | 一种建筑模型正误校验方法和系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |