CN111965694B - 地震物理点位置确定方法、装置及地震观测系统 - Google Patents
地震物理点位置确定方法、装置及地震观测系统 Download PDFInfo
- Publication number
- CN111965694B CN111965694B CN202010765897.2A CN202010765897A CN111965694B CN 111965694 B CN111965694 B CN 111965694B CN 202010765897 A CN202010765897 A CN 202010765897A CN 111965694 B CN111965694 B CN 111965694B
- Authority
- CN
- China
- Prior art keywords
- target
- grid
- xij
- range
- random number
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000006870 function Effects 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 16
- 238000001514 detection method Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 abstract description 13
- 238000003384 imaging method Methods 0.000 abstract description 9
- 238000009826 distribution Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 18
- 230000001788 irregular Effects 0.000 description 10
- 238000005070 sampling Methods 0.000 description 8
- 230000005284 excitation Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000001413 cellular effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 229940030980 inova Drugs 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
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/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种地震物理点位置确定方法、装置及地震观测系统，涉及地球物理勘探领域，该方法应用于地震观测系统，该方法包括：获取地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围；根据范围参数、目标地震物理点数量和目标地震物理点范围确定网格参数；根据网格参数生成目标网络；在目标网络的每个网格中随机生成一个目标地震物理点的坐标值。本发明可在满足目标地震物理点布设数量要求的同时，使目标地震物理点的位置在目标网络的每个网格中均匀分布，且坐标值在网格中随机生成，该坐标值可用于地震采集技术，以提高地震采集的成像精度。
Description
技术领域
本发明涉及地球物理勘探技术领域，尤其是涉及一种地震物理点位置确定方法、装置及地震观测系统。
背景技术
基于稀疏采样的地震采集技术是一种最新发展的地震采集新技术。该技术按照有利于提高成像精度的非规则采样规则，对地震数据进行噪音压制、重构，获得超越相同道密度规则采样的成像精度。地震观测系统应用于地震采集过程中，基于地震采集提高成像精度的要求，如何进行地震观测系统中的地震物理点非规则布设，是亟待解决的问题。
发明内容
本发明提供了一种地震物理点位置确定方法、装置及地震观测系统，可以通过网格确定用于布设目标地震物理点的坐标值，以便根据该坐标值布设地震物理点，有利于提高地震采集结果的成像精度。
第一方面，本发明实施例提供了一种地震物理点位置确定方法，应用于地震观测系统，该方法包括：获取所述地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围；根据所述范围参数、所述目标地震物理点数量和所述目标地震物理点范围确定网格参数；根据所述网格参数生成目标网络；在所述目标网络的每个网格中随机生成一个目标地震物理点的坐标值。
第二方面，本发明实施例还提供一种地震物理点位置确定装置，应用于地震观测系统，该装置包括：获取模块，用于获取所述地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围；确定模块，用于根据所述范围参数、所述目标地震物理点数量和所述目标地震物理点范围确定网格参数；生成模块，用于根据所述网格参数生成目标网络；位置模块，用于在所述目标网络的每个网格中随机生成一个目标地震物理点的坐标值。
第三方面，本发明实施例还提供一种地震观测系统，该系统包括上述任一种地震物理点位置确定装置。
第四方面，本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述地震物理点位置确定方法。
第五方面，本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述地震物理点位置确定方法的计算机程序。
本发明实施例带来了以下有益效果：本发明实施例提供了一种地震物理点位置确定方法、装置及地震观测系统，该方法应用于地震观测系统，该方法通过获取地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围，确定网格参数，再根据网格参数生成目标网络，最后在目标网络的每个网格中随机生成一个目标地震物理点的坐标值。本发明实施例可在满足目标地震物理点布设数量要求的同时，使目标地震物理点的位置在目标网络的每个网格中均匀分布，且坐标值在网格中随机生成，该坐标值可用于地震采集技术，以提高地震采集的成像精度。
本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的地震物理点位置确定方法流程图；
图2为本发明实施例提供的矩形网格的示意图；
图3为本发明实施例提供的正六边形网格的示意图；
图4为本发明实施例提供的一种判断点是否在正六边形内部的方法的示意图；
图5为本发明实施例提供的非规则观测系统物理点的示意图；
图6为本发明实施例提供的一种地震物理点位置确定装置结构框图；
图7为本发明实施例提供的另一种地震物理点位置确定装置结构框图；
图8为本发明实施例提供的计算机设备结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
美国Geospace公司，在2007年推出了一种新型的盲采式无线节点地震数据采集仪器——GSR(Geospace Seismic Recorder，地球空间地震记录仪)无线节点仪器，该设备受到BP公司的高度认可。BP公司通过融合ISS(International Space Station，国际空间站)技术，在实际项目中实践应用该设备，效果极为理想。此后，美国Fairfield公司开发了ZLAND点仪器，法国SERCEL公司研制推出了UNITE节点仪器，INOVA公司开发了HAWK和昆腾节点仪器，各种节点仪器不断被各个公司研发并推向市场。无线节点仪器的使用和推广，为地震勘探新技术、新方法的发展提供了坚实基础。
基于稀疏采样的地震采集技术正是在此基础上最新发展的地震采集新技术。该技术需要克服两个难题，一是怎么设计合理的非规则采样观测系统，二是非规则采样数据如何处理。
基于上述第一个问题，本发明实施例提供的一种地震物理点位置确定方法、装置及地震观测系统，该方法使该观测系统中物理点从整体分布来说是均匀的，在每一个相同大小的子区域内都有相同数目的物理点；而具体到每一个子区域内，物理点的坐标又是随机的。该方法可用于基于稀疏变换的三维非规则采样地震采集技术，能够提高地震采集的成像精度。
为便于对本实施例进行理解，首先对本发明实施例所公开的一种地震物理点位置确定方法进行详细介绍。
本发明实施例提供了一种地震物理点位置确定方法，参见图1所示的一种地震物理点位置确定方法流程图，该方法包括以下步骤：
步骤S102，获取地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围。
在本发明实施例中，地震观测系统的范围参数是地震观测系统中指用于布设地震物理点的位置范围数据。地震观测系统的范围参数包括第一方向上的最大值参数和最小值参数，以及第二方向上的最大值参数和最小值参数，从而根据地震观测系统的范围参数可以确定一个矩形范围，用于布设地震物理点。
目标地震物理点数量是计划在地震观测系统的矩形范围中布设的地震物理点的数量。根据目标地震物理点范围可以确定多个目标地震物理点在地震观测系统的矩形范围中需求的最大的面积。
需要说明的是，在本发明实施例中地震物理点可以包括检波点或炮点。
步骤S104，根据范围参数、目标地震物理点数量和目标地震物理点范围确定网格参数。
在本发明实施例中，在得到范围参数、目标地震物理点数量和目标地震物理点范围后，通过计算，可以确定建立网格所需要的相关数据，得到网格参数。
需要说明的是，由于地震物理点可以包括检波点或炮点，因此，可以分别对检波点和炮点生成网格参数。
步骤S106，根据网格参数生成目标网络。
在本发明实施例中，每个目标网络中包括多个网格。可以根据不同的网格参数生成不同的目标网络。例如，由于地震物理点可以包括检波点或炮点，可以根据网格参数生成检波点网络和激发点网络；检波点网络可用于确定目标检波点的坐标值；激发点网络可用于确定目标炮点的坐标值。根据观测系统范围及规模(检波点和激发点的数量)，生成XOY平面上的检波点网络和激发点网络，两个网络的位置可以重合(即使用同一个目标网络)，也可以错开。
步骤S108，在目标网络的每个网格中随机生成一个目标地震物理点的坐标值。
在本发明实施例中，目标网络中包括的网格数量与目标地震物理点数量相同，因此，每个网格中包括一个坐标值，从而目标地震物理点的位置按照目标网络均匀分布在地震观测系统的范围内，并且每个位置的坐标值在各自所在的网格中随机分布。
本发明实施例提供了一种地震物理点位置确定方法，该方法应用于地震观测系统，该方法通过获取地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围，确定网格参数，再根据网格参数生成目标网络，最后在目标网络的每个网格中随机生成一个目标地震物理点的坐标值。本发明实施例可在满足目标地震物理点布设数量要求的同时，使目标地震物理点的位置在目标网络的每个网格中均匀分布，且坐标值在网格中随机生成，该坐标值可用于地震采集技术，以提高地震采集的成像精度。
为了满足地质任务要求，获取地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围之前，还可以执行如下步骤：
获取地质模型的空间范围数据；根据地质模型的空间范围数据确定地震观测系统的范围参数。
在本发明实施例中，根据地质任务要求和已建立的基于空间直角坐标系的立方体地质模型，确定地质模型的空间范围(Xmin，Ymin，Zmin)、(Xmax，Ymax，Zmax)。其中，Xmin、Xmax分别是地质模型X方向的最小值、最大值，Ymin、Ymax分别是地质模型Y方向的最小值、最大值，Zmin、Zmax分别是地质模型Z方向的最小值、最大值。
获取该地质模型的空间范围数据，确定地震观测系统的范围参数。根据地质任务要求，确定观测系统在XOY平面上的范围(xmin，ymin)、(xmax，ymax)。其中，xmin、xmax分别是观测系统X方向的最小值、最大值，ymin、ymax分别是观测系统Y方向的最小值、最大值，并且观测系统不能超出地质模型的范围，即xmin≥Xmin，xmax≤Xmax，ymin≥Ymin，ymax≤Ymax。
为了合理规划范围参数所确定的面积，根据范围参数、目标地震物理点数量和目标地震物理点范围确定网格参数，可以按照如下步骤执行：
根据范围参数确定网格原点坐标位置；根据原点坐标位置、目标地震物理点数量和目标地震物理点范围确定网格数量以及网格尺寸；将网格原点坐标位置、网格数量以及网格尺寸作为网格参数。
在本发明实施例中，根据范围参数确定了矩形区域后，基于该矩形区域确定要建立的网格的远点坐标位置(x0，y0)，再基于远点坐标位置、目标地震物理点数量和目标地震物理点范围确定在第一方向上的网格数量nx，第二方向上的网格数量ny，以及网格尺寸。
为了提升该方法的适用范围，目标网络包括矩形网格或正六边形网格；确定网格尺寸，包括：确定矩形网格的长度和宽度，或确定正六边形网格的边长。
在本发明实施例中，参见图2所示的矩形网格的示意图，若建立的目标网络包括矩形网格，则确定网格尺寸包括确定矩形网格的长度dx和宽度dy。参见图3所示的正六边形网格的示意图，若建立的目标网络包括正六边形网格，则确定网格尺寸包括确定蜂窝状正六边形网格的边长a。
需要说明的是，在本发明实施例中，参见图2，矩形网格每个格子大小形状相同，均为矩形，参见图3，正六边形网格为蜂窝状，每个格子大小形状相同，均为正六边形。
在目标网络包括蜂窝状网格时，还可以执行如下步骤：
根据目标地震物理点的坐标值和正六边形网格的边长判断目标地震物理点的坐标值是否在正六边形网格内部；如果是，记录目标地震物理点的坐标值；如果否，重新根据目标网络确定目标地震物理点的坐标值。
在蜂窝状网格中，还需要判断物理点(xij，yij)是否在正六边形Gij内。参见图4所示的一种判断点是否在正六边形内部的方法的示意图，边长为a的正六边形中点为XOY坐标系原点O，判断点P(x，y)是否在正六边形内，x≥0，y≥0。
(1)判断P是否在正六边形的外包矩形内，若x≥a或者则点P不在正六边形内部，否则继续下一步判断；
(2)若MP≥MN，则P在正六边形内部(包括在边线上)；若MP＜MN，则P在正六边形外部。其中，MP＝a-x，
根据上述方法判断，若P在正六边形外部，则重新根据目标网络确定目标地震物理点的坐标值，若P在正六边形内部，则正六边形Gij内的物理点(xij，yij)即为随机物理点，记录该目标地震物理点的坐标值。
为了提升坐标值的随机性，在目标网络的每个网格中随机生成一个目标地震物理点的坐标值，可以按照如下步骤执行：
确定目标网络中每个网格的坐标范围；根据坐标范围生成第一随机数和第二随机数；根据坐标范围、第一随机数和第二随机数生成目标地震物理点的坐标值。
在本发明实施例中，生成的目标网络中每个格子编号为Gij，其中，i＝1、2......nx，为格子的列号，j＝1、2......ny，为格子的行号，根据格子的行列号即可确定每个网格的坐标范围。根据坐标范围对应每个网格生成一组随机数，包括第一随机数和第二随机数，再根据每个网格的坐标范围及对应的两个随机数生成一个地震物理点的坐标值，得到每个网格的目标地震物理点的坐标值。
在本发明实施例中，使用产生二维随机数(uxij，uyij)的方法生成非规则观测系统物理点的坐标。在网格的格子Gij中，有且仅有一个物理点，该物理点的坐标为(xij，yij)，其中，i＝1、2......nx，j＝1、2......ny，xij和yij分别独立。根据实际施工测量的精度需求以及后续的数据处理需要，xij和yij的精度可以达到米级甚至厘米级。
可使用基于下面几种随机数产生方法之一来生成物理点坐标值。
可以按照如下公式根据坐标范围生成第一随机数和第二随机数：uxij＝rand(xij)/randMax，uyij＝rand(yij)/randMax，其中，uxij是第一随机数，uyij是第二随机数，rand(xij)是根据随机函数生成的第一随机正整数，rand(yij)是根据随机函数生成的第二随机正整数，randMax是随机函数能返回的最大正整数。
在本发明实施例中，高级编程语言自带的随机函数，记为rand，根据该随机函数和上述公式进行计算，产生第一随机数和第二随机数。得到的uxij和uyij则是[0，1]上的随机数。
可以按照如下公式根据坐标范围生成第一随机数和第二随机数：vn＝(b×vn-1+c)mod(m)，un＝vn/m，其中，b、c和m为给定的常数，v0是任意非负整数，un为[0，1]上的随机数。
在本发明实施例中，可以使用小林混合同余法，按照上述公式产生第一随机数和第二随机数。其中，n＝1，2，3......，b＝314159269，c＝453806245，m＝231，mod是求模的标识，v0可以是任意非负整数，un即为[0，1]上的随机数。此方法运用了初等数论中的同余运算，产生的随机数在[0，1]上均匀分布。
可以按照如下公式根据坐标范围、第一随机数和第二随机数生成目标地震物理点的坐标值：xij＝minxij+(maxxij-minxij)×uxij×uxij，yij＝minyij+(maxyij-minyij)×uyij×uyij，其中，xij是目标地震物理点在第一方向上的坐标值，yij是目标地震物理点在第二方向上的坐标值，uxij是第一随机数，uyij是第二随机数，minxij是目标网络在第一方向上的最小值，minyij是目标网络在第二方向上的最小值，maxxij是目标网络在第一方向上的最大值，maxyij是目标网络在第二方向上的最大值。
在本发明实施例中，xij是[minxij，maxxij]上的随机数，yij是[minyij，maxyij]上的随机数。
可以按照如下公式根据坐标范围、第一随机数和第二随机数生成目标地震物理点的坐标值：xij＝minxij+(maxxij-minxij)×uxij，yij＝minyij+(maxyij-minyij)×uyij，其中，xij是目标地震物理点在第一方向上的坐标值，yij是目标地震物理点在第二方向上的坐标值，uxij是第一随机数，uyij是第二随机数，minxij是目标网络在第一方向上的最小值，minyij是目标网络在第二方向上的最小值，maxxij是目标网络在第一方向上的最大值，maxyij是目标网络在第二方向上的最大值。
在本发明实施例中，使用小林(Kobayashi)归一平方法，基于小林(Kobayashi)混合同余法得到的第一随机数和第二随机数，按照上述公式得到的xij和yij，在均匀分布的基础上，随机性更强。
下面以一个具体实施例对该方法进行说明。
(1)建立X、Y方向范围都为5000米，Z方向范围为4000米的地质模型。
(2)确定观测系统在XOY平面上的范围，检波点范围为X方向(100.0，5000.0)、Y方向(100.0，5000.0)，激发点只设置一个点，在检波点范围的中心位置。
(3)使用正六边形网格，边长为23.1，检波点网格的原点坐标为(100.0，100.0)，X方向网格数为140，Y方向网格数为120。
(4)使用Kobayashi混合同余法，生成如图5中所示的非规则观测系统物理点，图中方形点是检波点，圆点是激发点(仅一个)。
本发明实施例提供了一种地震物理点位置确定方法、装置及地震观测系统，该方法基于地质模型和观测系统规模，建立三维地质模型，确定地质模型范围；确定基于地质模型的三维观测系统的范围及规模；根据观测系统范围及规模，确定网格参数，包括原点坐标、行列数及网格尺寸；根据网格参数，使用产生二维随机数的方法生成非规则物理点。该方法可用于基于稀疏变换的三维非规则采样地震采集技术，能够提高地震采集的成像精度。
本发明实施例中还提供了一种地震物理点位置确定装置，如下面的实施例所述。由于该装置解决问题的原理与地震物理点位置确定方法相似，因此该装置的实施可以参见地震物理点位置确定方法的实施，重复之处不再赘述。参见图6所示的一种地震物理点位置确定装置结构框图，该装置包括：
获取模块71，用于获取地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围；确定模块72，用于根据范围参数、目标地震物理点数量和目标地震物理点范围确定网格参数；生成模块73，用于根据网格参数生成目标网络；位置模块74，用于在目标网络的每个网格中随机生成一个目标地震物理点的坐标值。
在一个实施例中，参见图7所示的另一种地震物理点位置确定装置结构框图，该装置还包括模型模块75，用于：获取地质模型的空间范围数据；根据地质模型的空间范围数据确定地震观测系统的范围参数。
在一个实施例中，确定模块，具体用于：根据范围参数确定网格原点坐标位置；根据原点坐标位置、目标地震物理点数量和目标地震物理点范围确定网格数量以及网格尺寸；将网格原点坐标位置、网格数量以及所述网格尺寸作为网格参数。
在一个实施例中，目标网络包括矩形网格或正六边形网格；确定模块，具体用于：确定矩形网格的长度和宽度，或确定正六边形网格的边长。
在一个实施例中，该装置还包括判断模块76，用于：根据目标地震物理点的坐标值和正六边形网格的边长判断目标地震物理点的坐标值是否在正六边形网格内部；如果是，记录目标地震物理点的坐标值；如果否，重新根据目标网络确定目标地震物理点的坐标值。
在一个实施例中，位置模块，具体用于：确定目标网络中每个网格的坐标范围；根据坐标范围生成第一随机数和第二随机数；根据坐标范围、第一随机数和第二随机数生成目标地震物理点的坐标值。
在一个实施例中，位置模块，具体用于：按照如下公式根据坐标范围生成第一随机数和第二随机数：uxij＝rand(xij)/randMax，uyij＝rand(yij)/randMax，其中，uxij是第一随机数，uyij是第二随机数，rand(xij)是根据随机函数生成的第一随机正整数，rand(yij)是根据随机函数生成的第二随机正整数，randMax是随机函数能返回的最大正整数。
在一个实施例中，位置模块，具体用于：按照如下公式根据坐标范围生成第一随机数和第二随机数：vn＝(b×vn-1+c)mod(m)，un＝vn/m，其中，b、c和m为给定的常数，v0是任意非负整数，un为[0，1]上的随机数。
在一个实施例中，位置模块，具体用于：按照如下公式根据坐标范围、第一随机数和第二随机数生成目标地震物理点的坐标值：xij＝minxij+(maxxij-minxij)×uxij×uxij，yij＝minyij+(maxyij-minyij)×uyij×uyij，其中，xij是目标地震物理点在第一方向上的坐标值，yij是目标地震物理点在第二方向上的坐标值，uxij是第一随机数，uyij是第二随机数，minxij是目标网络在第一方向上的最小值，minyij是目标网络在第二方向上的最小值，maxxij是目标网络在第一方向上的最大值，maxyij是目标网络在第二方向上的最大值。
在一个实施例中，位置模块，具体用于：按照如下公式根据坐标范围、第一随机数和第二随机数生成目标地震物理点的坐标值：xij＝minxij+(maxxij-minxij)×uxij，yij＝minyij+(maxyij-minyij)×uyij，其中，xij是目标地震物理点在第一方向上的坐标值，yij是目标地震物理点在第二方向上的坐标值，uxij是第一随机数，uyij是第二随机数，minxij是目标网络在第一方向上的最小值，minyij是目标网络在第二方向上的最小值，maxxij是目标网络在第一方向上的最大值，maxyij是目标网络在第二方向上的最大值。
本发明实施例中还提供了一种地震观测系统，该系统包括上述任一种地震物理点位置确定装置。
本发明实施例还提供一种计算机设备，参见图8所示的计算机设备结构示意框图，该计算机设备包括存储器81、处理器82及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述任一种地震物理点位置确定方法的步骤。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的计算机设备的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
本发明实施例还提供一种计算机可读存储介质，计算机可读存储介质存储有执行上述任一种地震物理点位置确定方法的计算机程序。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (19)
1.一种地震物理点位置确定方法，其特征在于，应用于地震观测系统，所述方法包括：
获取所述地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围；
根据所述范围参数、所述目标地震物理点数量和所述目标地震物理点范围确定网格参数；
根据所述网格参数生成目标网络；
在所述目标网络的每个网格中随机生成一个目标地震物理点的坐标值；
地震物理点包括检波点或炮点；
根据所述范围参数、所述目标地震物理点数量和所述目标地震物理点范围确定网格参数，包括：
根据所述范围参数确定网格原点坐标位置；
根据所述原点坐标位置、所述目标地震物理点数量和所述目标地震物理点范围确定网格数量以及网格尺寸；
将所述网格原点坐标位置、所述网格数量以及所述网格尺寸作为网格参数；
所述目标网络包括矩形网格或正六边形网格；
确定网格尺寸，包括：
确定所述矩形网格的长度和宽度，或确定所述正六边形网格的边长；
目标网络中包括的网格数量与目标地震物理点数量相同。
2.根据权利要求1所述的方法，其特征在于，获取所述地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围之前，还包括：
获取地质模型的空间范围数据；
根据所述地质模型的空间范围数据确定地震观测系统的范围参数。
3.根据权利要求1所述的方法，其特征在于，还包括：
根据所述目标地震物理点的坐标值和所述正六边形网格的边长判断所述目标地震物理点的坐标值是否在正六边形网格内部；
如果是，记录所述目标地震物理点的坐标值；
如果否，重新根据所述目标网络确定目标地震物理点的坐标值。
4.根据权利要求1所述的方法，其特征在于，在所述目标网络的每个网格中随机生成一个目标地震物理点的坐标值，包括：
确定所述目标网络中每个网格的坐标范围；
根据所述坐标范围生成第一随机数和第二随机数；
根据所述坐标范围、所述第一随机数和所述第二随机数生成目标地震物理点的坐标值。
5.根据权利要求4所述的方法，其特征在于，包括按照如下公式根据所述坐标范围生成第一随机数和第二随机数：
uxij＝rand(xij)/randMax
uyij＝rand(yij)/randMax
其中，uxij是第一随机数，uyij是第二随机数，rand(xij)是根据随机函数生成的第一随机正整数，rand(yij)是根据随机函数生成的第二随机正整数，randMax是随机函数能返回的最大正整数。
6.根据权利要求4所述的方法，其特征在于，包括按照如下公式根据所述坐标范围生成第一随机数和第二随机数：
vn＝(b×vn-1+c)mod(m)
un＝vn/m
其中，b、c和m为给定的常数，v0是任意非负整数，un为[0,1]上的随机数。
7.根据权利要求6所述的方法，其特征在于，包括按照如下公式根据所述坐标范围、所述第一随机数和所述第二随机数生成目标地震物理点的坐标值：
xij＝minxij+(maxxij-minxij)×uxij×uxij
yij＝minyij+(maxyij-minyij)×uyij×uyij
其中，xij是目标地震物理点在第一方向上的坐标值，yij是目标地震物理点在第二方向上的坐标值，uxij是第一随机数，uyij是第二随机数，minxij是目标网络在第一方向上的最小值，minyij是目标网络在第二方向上的最小值，maxxij是目标网络在第一方向上的最大值，maxyij是目标网络在第二方向上的最大值。
8.根据权利要求5或6所述的方法，其特征在于，包括按照如下公式根据所述坐标范围、所述第一随机数和所述第二随机数生成目标地震物理点的坐标值：
xij＝minxij+(maxxij-minxij)×uxij
yij＝minyij+(maxyij-minyij)×uyij
其中，xij是目标地震物理点在第一方向上的坐标值，yij是目标地震物理点在第二方向上的坐标值，uxij是第一随机数，uyij是第二随机数，minxij是目标网络在第一方向上的最小值，minyij是目标网络在第二方向上的最小值，maxxij是目标网络在第一方向上的最大值，maxyij是目标网络在第二方向上的最大值。
9.一种地震物理点位置确定装置，其特征在于，应用于地震观测系统，所述装置包括：
获取模块，用于获取所述地震观测系统的范围参数、目标地震物理点数量和目标地震物理点范围；
确定模块，用于根据所述范围参数、所述目标地震物理点数量和所述目标地震物理点范围确定网格参数；
生成模块，用于根据所述网格参数生成目标网络；
位置模块，用于在所述目标网络的每个网格中随机生成一个目标地震物理点的坐标值；
地震物理点包括检波点或炮点；
所述确定模块，具体用于：
根据所述范围参数确定网格原点坐标位置；
根据所述原点坐标位置、所述目标地震物理点数量和所述目标地震物理点范围确定网格数量以及网格尺寸；
将所述网格原点坐标位置、所述网格数量以及所述网格尺寸作为网格参数；
所述目标网络包括矩形网格或正六边形网格；所述确定模块，具体用于：
确定所述矩形网格的长度和宽度，或确定所述正六边形网格的边长；
目标网络中包括的网格数量与目标地震物理点数量相同。
10.根据权利要求9所述的装置，其特征在于，还包括模型模块，用于：
获取地质模型的空间范围数据；
根据所述地质模型的空间范围数据确定地震观测系统的范围参数。
11.根据权利要求10所述的装置，其特征在于，还包括判断模块，用于：
根据所述目标地震物理点的坐标值和所述正六边形网格的边长判断所述目标地震物理点的坐标值是否在正六边形网格内部；
如果是，记录所述目标地震物理点的坐标值；
如果否，重新根据所述目标网络确定目标地震物理点的坐标值。
12.根据权利要求9所述的装置，其特征在于，所述位置模块，具体用于：
确定所述目标网络中每个网格的坐标范围；
根据所述坐标范围生成第一随机数和第二随机数；
根据所述坐标范围、所述第一随机数和所述第二随机数生成目标地震物理点的坐标值。
13.根据权利要求12所述的装置，其特征在于，所述位置模块，具体用于：
按照如下公式根据所述坐标范围生成第一随机数和第二随机数：
uxij＝rand(xij)/randMax
uyij＝rand(yij)/randMax
其中，uxij是第一随机数，uyij是第二随机数，rand(xij)是根据随机函数生成的第一随机正整数，rand(yij)是根据随机函数生成的第二随机正整数，randMax是随机函数能返回的最大正整数。
14.根据权利要求12所述的装置，其特征在于，所述位置模块，具体用于：
按照如下公式根据所述坐标范围生成第一随机数和第二随机数：
vn＝(b×vn-1+c)mod(m)
un＝vn/m
其中，b、c和m为给定的常数，v0是任意非负整数，un为[0,1]上的随机数。
15.根据权利要求14所述的装置，其特征在于，所述位置模块，具体用于：
按照如下公式根据所述坐标范围、所述第一随机数和所述第二随机数生成目标地震物理点的坐标值：
xij＝minxij+(maxxij-minxij)×uxij×uxij
yij＝minyij+(maxyij-minyij)×uyij×uyij
其中，xij是目标地震物理点在第一方向上的坐标值，yij是目标地震物理点在第二方向上的坐标值，uxuj是第一随机数，uyij是第二随机数，minxij是目标网络在第一方向上的最小值，minyij是目标网络在第二方向上的最小值，maxxij是目标网络在第一方向上的最大值，maxyij是目标网络在第二方向上的最大值。
16.根据权利要求13或14所述的装置，其特征在于，所述位置模块，具体用于：
按照如下公式根据所述坐标范围、所述第一随机数和所述第二随机数生成目标地震物理点的坐标值：
xij＝minxij+(maxxij-minxij)×uxij
yij＝minyij+(maxyij-minyoj)×uyij
其中，xij是目标地震物理点在第一方向上的坐标值，yij是目标地震物理点在第二方向上的坐标值，uxij是第一随机数，uyij是第二随机数，minxij是目标网络在第一方向上的最小值，minxij是目标网络在第二方向上的最小值，maxxij是目标网络在第一方向上的最大值，maxyij是目标网络在第二方向上的最大值。
17.一种地震观测系统，其特征在于，包括权利要求9-16任一所述地震物理点位置确定装置。
18.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至8任一所述方法。
19.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至8任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010765897.2A CN111965694B (zh) | 2020-08-03 | 2020-08-03 | 地震物理点位置确定方法、装置及地震观测系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010765897.2A CN111965694B (zh) | 2020-08-03 | 2020-08-03 | 地震物理点位置确定方法、装置及地震观测系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111965694A CN111965694A (zh) | 2020-11-20 |
| CN111965694B true CN111965694B (zh) | 2024-01-30 |
Family
ID=73363634
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010765897.2A Active CN111965694B (zh) | 2020-08-03 | 2020-08-03 | 地震物理点位置确定方法、装置及地震观测系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111965694B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5271014A (en) * | 1992-05-04 | 1993-12-14 | International Business Machines Corporation | Method and apparatus for a fault-tolerant mesh with spare nodes |
| CN102184298A (zh) * | 2011-05-18 | 2011-09-14 | 山东大学 | 金属体积塑性成形中有限元分析刚度矩阵存储与生成方法 |
| CA2806241A1 (en) * | 2011-07-26 | 2012-02-09 | Bp Corporation North America Inc. | Method for separating independent simultaneous sources |
| CN103279991A (zh) * | 2013-04-16 | 2013-09-04 | 西南石油大学 | 一种利用离散裂缝端点变形数模提高油藏开发效果的方法 |
| CN103559740A (zh) * | 2013-08-26 | 2014-02-05 | 空气动力学国家重点实验室 | 一种可实现交互操作的棱柱网格生成方法 |
| CN104656126A (zh) * | 2015-02-26 | 2015-05-27 | 中国石油天然气集团公司 | 一种生成地震观测系统的方法 |
| CN107561574A (zh) * | 2017-07-20 | 2018-01-09 | 中国石油天然气集团公司 | 炮点位置、接收点位置的检测方法和装置 |
| CN109143338A (zh) * | 2018-08-13 | 2019-01-04 | 中国石油天然气集团有限公司 | 一种地震数据合并的方法及装置 |
| CN111352152A (zh) * | 2018-12-21 | 2020-06-30 | 中国石油天然气集团有限公司 | 地震数据观测系统快速质控方法及装置 |
| CN111474574A (zh) * | 2019-01-23 | 2020-07-31 | 中国石油天然气集团有限公司 | 基于压缩感知的地震采集观测系统的生成方法及装置 |
| WO2021216773A1 (en) * | 2020-04-23 | 2021-10-28 | The Nielsen Company (Us), Llc | Methods and apparatus to encrypt media for identification |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8559270B2 (en) * | 2008-08-15 | 2013-10-15 | Bp Corporation North America Inc. | Method for separating independent simultaneous sources |
| US11269092B2 (en) * | 2018-06-21 | 2022-03-08 | Sercel Sas | Method and system for optimizing seismic data acquisition using compressed sensing |
-
2020
- 2020-08-03 CN CN202010765897.2A patent/CN111965694B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5271014A (en) * | 1992-05-04 | 1993-12-14 | International Business Machines Corporation | Method and apparatus for a fault-tolerant mesh with spare nodes |
| CN102184298A (zh) * | 2011-05-18 | 2011-09-14 | 山东大学 | 金属体积塑性成形中有限元分析刚度矩阵存储与生成方法 |
| CA2806241A1 (en) * | 2011-07-26 | 2012-02-09 | Bp Corporation North America Inc. | Method for separating independent simultaneous sources |
| CN103279991A (zh) * | 2013-04-16 | 2013-09-04 | 西南石油大学 | 一种利用离散裂缝端点变形数模提高油藏开发效果的方法 |
| CN103559740A (zh) * | 2013-08-26 | 2014-02-05 | 空气动力学国家重点实验室 | 一种可实现交互操作的棱柱网格生成方法 |
| CN104656126A (zh) * | 2015-02-26 | 2015-05-27 | 中国石油天然气集团公司 | 一种生成地震观测系统的方法 |
| CN107561574A (zh) * | 2017-07-20 | 2018-01-09 | 中国石油天然气集团公司 | 炮点位置、接收点位置的检测方法和装置 |
| CN109143338A (zh) * | 2018-08-13 | 2019-01-04 | 中国石油天然气集团有限公司 | 一种地震数据合并的方法及装置 |
| CN111352152A (zh) * | 2018-12-21 | 2020-06-30 | 中国石油天然气集团有限公司 | 地震数据观测系统快速质控方法及装置 |
| CN111474574A (zh) * | 2019-01-23 | 2020-07-31 | 中国石油天然气集团有限公司 | 基于压缩感知的地震采集观测系统的生成方法及装置 |
| WO2021216773A1 (en) * | 2020-04-23 | 2021-10-28 | The Nielsen Company (Us), Llc | Methods and apparatus to encrypt media for identification |
Non-Patent Citations (3)
| Title |
|---|
| 利用地方震、远震走时和面波数据联合反演日本俯冲带P波和S波层析成像;X.Liu;D.P.Zhao;张晓曼;;世界地震译丛(第01期);全文 * |
| 王俊 ; 李树强 ; 刘刚 ; .基于相似度的温室无线传感器网络定位算法.农业工程学报.2013,(第22期),全文. * |
| 程铭东 ; 占亮 ; 赵东方 ; .应用于目标检测的Ad Hoc网络节点的优化部署.华中师范大学学报(自然科学版).2008,(第02期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111965694A (zh) | 2020-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xu et al. | A new computer code for discrete fracture network modelling | |
| US8027812B2 (en) | Charge trajectory calculating method, system, and program | |
| CN102057368B (zh) | 利用最大连续场在三维体积模型中分布性质 | |
| CN106896403B (zh) | 弹性高斯束偏移成像方法和系统 | |
| CN112800657A (zh) | 基于复杂地形的重力场数值模拟方法、装置和计算机设备 | |
| CN109887073B (zh) | 岩心的三维数字模型构建方法及装置 | |
| Nan et al. | Groundwater parameter estimation using the ensemble Kalman filter with localization | |
| CA3152669A1 (en) | Reservoir-based modeling method and device for pore network model | |
| CN111830557B (zh) | 基于压裂微地震的人工裂缝复杂度指数获取方法及系统 | |
| CN105388513A (zh) | 地震采集观测系统的建立方法和装置 | |
| Yildirim et al. | Stochastic simulations of ocean waves: An uncertainty quantification study | |
| CN108898013B (zh) | 一种基于布局划分特征向量的安卓应用界面相似度比较方法 | |
| CN111965694B (zh) | 地震物理点位置确定方法、装置及地震观测系统 | |
| Manchuk et al. | Implementation aspects of sequential Gaussian simulation on irregular points | |
| Fu et al. | Spatial topology identification of three-dimensional complex block system of rock masses | |
| CN109633781B (zh) | 地质属性获取方法及装置、电子设备、存储介质 | |
| CN115586575B (zh) | 一种3d被动源拟炮集计算方法、系统、设备及存储介质 | |
| CN111474574A (zh) | 基于压缩感知的地震采集观测系统的生成方法及装置 | |
| CN110687598A (zh) | 一种加速微震数值模拟的方法及装置 | |
| Sajid et al. | Box dimension of the border of Kingdom of Saudi Arabia | |
| CN114065577A (zh) | 一种直流电阻率小波伽辽金三维正演方法 | |
| CN114491897A (zh) | 地震波数值模拟方法、装置、介质及电子设备 | |
| CN112379413A (zh) | 基于能量谱等效的非规则震源表征方法和装置 | |
| CN113740908B (zh) | 地震切片的二维变差分析方法、电子设备及介质 | |
| D'Ambrosio et al. | Applications of the opencal scientific library in the context of cfd: Applications to debris flows |
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 |