CN112379409A - 一种陆上无线节点地震仪器桩号匹配方法及装置 - Google Patents
一种陆上无线节点地震仪器桩号匹配方法及装置 Download PDFInfo
- Publication number
- CN112379409A CN112379409A CN202011201851.4A CN202011201851A CN112379409A CN 112379409 A CN112379409 A CN 112379409A CN 202011201851 A CN202011201851 A CN 202011201851A CN 112379409 A CN112379409 A CN 112379409A
- Authority
- CN
- China
- Prior art keywords
- point
- geometric center
- center point
- distance
- wireless node
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 238000004590 computer program Methods 0.000 claims description 19
- 238000003860 storage Methods 0.000 claims description 10
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 12
- 238000009826 distribution Methods 0.000 abstract description 9
- 238000010586 diagram Methods 0.000 description 18
- 238000012545 processing Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000007621 cluster analysis Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000005433 ionosphere Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/14—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring distance or clearance between spaced objects or spaced apertures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Radar, Positioning & Navigation (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例提供的陆上无线节点地震仪器桩号匹配方法及装置，利用每个点距中心点的距离与距离平均值做比较，如果大于距离平均值则标记当前点状态为1，否则标记当前点状态为0。利用状态标记为1的点重新计算几何中心点。该方法通过标记点的状态，标记出离散点、飞点的状态为0，将这些点不参与新几何中心点的计算。完全消除了离散点、飞点对计算新几何中心点的影响，克服了反距离权重算法中离散点、飞点对桩号匹配的影响。通过循环迭代计算，几何中心点逐渐向点分布密度大的地方移动。等循环终止，以最后计算的几何中心点位置最优，以该中心点为圆心以特定距离为半径(通常为半个道距)的圆包含的点集数量最多。
Description
技术领域
本发明涉及地质勘探技术领域，具体涉及一种陆上无线节点地震仪器桩号匹配方法及装置。
背景技术
陆上石油地球物理勘探接收地震波信号仪器从有线检波器向陆上无线节点方向发展。陆上无线节点在陆上石油地球勘探中应用越来越广。陆上无线节点仪器在野外作业时，从激活状态开始就持续接收GPS信号，形成节点坐标记录点集。为了将节点采集到的坐标数据和地震测线SPS检波点进行关联，需要计算陆上无线节点坐标的几何中心位置，然后地震测线SPS桩号进行关联匹配。
陆上无线节点地震仪器在野外工作时，陆上无线节点仪器采集的坐标点因系统初始化、卫星信号的多路径效应、卫星信号的电离层干扰、周围高压电线电磁波影响等因素干扰，以及陆上无线节点一起从作业现场回收到装载到上位机上下传数据时，依然在进行坐标记录，因此造成坐标点分布呈发散状、离散状。
因上述情况的存在，节点设备记录了干扰噪声坐标数据，这些陆上无线节点发散的、离散的点位坐标影响了桩号匹配精度和成功率。
发明内容
针对现有技术中的问题，本发明实施例提供一种陆上无线节点地震仪器桩号匹配方法及装置，能够解决陆上地震勘探使用陆上无线节点地震数据采集业务中，陆上无线节点号与物理检波点号之间的匹配问题。
一方面，一种陆上无线节点地震仪器桩号匹配方法，包括：
解析陆上无线节点QC文件记录的点位坐标，并根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离；
根据所述平均距离标记点集中点的状态，得到第一状态点和第二状态点；
根据所述第二状态点重新计算几何中心点，若重新计算得到的更新几何中心点与原始的几何中心点之间的距离大于设定阈值，则用所述更新几何中心点替换原始的几何中心点，直至两者小于设定阈值后，将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。
在优选的实施例中，所述点位包括：坐标东坐标值和北坐标值，所述根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离，包括：
累加每个陆上无线节点的东坐标值，除以点集总数，计算出东坐标的平均值；累加每个陆上无线节点的北坐标值，除以点集总数，计算出北坐标的平均值；
以东坐标和北坐标平均值构造几何中心点；
计算点集中每个点和几何中心点的距离，逐个累加距离值并除以总数，求出点集中每个点距离中心点的距离平均值。
在优选的实施例中，所述第一状态点与所述几何中心点的距离大于第一设定距离，所述第二状态点与所述几何中心点的距离小于所述第一设定距离。
在优选的实施例中，在根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离之前，还包括：
将节点记录的大地坐标转换为平面直角坐标。
另一方面，一种陆上无线节点地震仪器桩号匹配装置，包括：
几何中心点计算模块，解析陆上无线节点QC文件记录的点位坐标，并根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离；
标记模块，遍历点集中所有点据中心点的距离，与上所述平均距离做比较，分离、标记点集中各个点的状态，得到第一状态点和第二状态点；
迭代匹配模块，根据所述第二状态点重新计算几何中心点，若重新计算得到的更新几何中心点与原始的几何中心点之间的距离大于设定阈值，则用所述更新几何中心点替换原始的几何中心点，直至两者小于设定阈值后，将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。
在优选的实施例中，所述点位包括：坐标东坐标值和北坐标值，所述几何中心点计算模块，包括：
累加单元，累加每个陆上无线节点的东坐标值，除以点集总数，计算出东坐标的平均值；累加每个陆上无线节点的北坐标值，除以点集总数，计算出点集北坐标的平均值；
几何中心点构造单元，以点集的东坐标和北坐标平均值构造几何中心点；
距离平均值计算单元，计算点集中每个点和几何中心点的距离，逐个累加距离值并除以总数，求出点集中每个点距离中心点的距离平均值。
在优选的实施例中，所述第一状态点与所述几何中心点的距离大于第一设定距离，所述第二状态点与所述几何中心点的距离小于所述第一设定距离。
在优选的实施例中，还包括：
坐标转换模块，将大地坐标转换为平面直角坐标。
再一方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任一实施例所述陆上无线节点地震仪器桩号匹配方法的步骤。
又一方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任一实施例所述陆上无线节点地震仪器桩号匹配方法的步骤。
本发明实施例提供的陆上无线节点地震仪器桩号匹配方法及装置，利用每个点距中心点的距离与距离平均值做比较，如果大于距离平均值则标记当前点状态为1，否则标记当前点状态为0。利用状态标记为1的点重新计算几何中心点。该方法通过标记点的状态，标记出离散点、飞点的状态为0，将这些点不参与新几何中心点的计算。完全消除了离散点、飞点对计算新几何中心点的影响，克服了反距离权重算法中离散点、飞点对桩号匹配的影响。通过循环迭代计算，几何中心点逐渐向点分布密度大的地方移动。等循环终止，以最后计算的几何中心点位置最优，以该中心点为圆心以特定距离为半径(通常为半个道距)的圆包含的点集数量最多。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1是本发明一个实施例提供的基于改进的kmeans算法的陆上无线节点桩号匹配计算流程。
图2是本发明一个实施例提供的00045站采集4393个点坐标分布图。
图3是本发明一个实施例提供的00045采集站通过聚类分析计算中心点图。
图4是本发明一个实施例提供的采集站5米范围内包含4245点数示意图。
图5是本发明一个实施例提供的00086站采集3555个点坐标分布图。
图6是本发明一个实施例提供的00086采集站通过聚类分析计算中心点图。
图7是本发明一个实施例提供的采集站5米范围内包含3471点数示意图。
图8是本发明一个实施例提供的计算的节点坐标中心点示意图。
图9是本发明一个实施例提供的地震测线SPS桩号分布图。
图10是本发明一个实施例提供的节点中心点和SPS桩号匹配效果图。
图11是本发明一个实施例提供的陆上无线节点地震仪器桩号匹配方法流程示意图。
图12是本发明一个实施例提供的实现陆上无线节点地震仪器桩号匹配方法的装置结构示意图。
图13是本发明又一个实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
图11是本发明一个实施例提供的陆上无线节点地震仪器桩号匹配方法的流程示意图，如图11所示，本发明实施例提供的陆上无线节点地震仪器桩号匹配方法，包括：
S1：解析陆上无线节点QC文件记录的点位坐标，并根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离；
S2：根据所述平均距离标记点集中点的状态，得到第一状态点和第二状态点；
S3：根据所述第二状态点重新计算几何中心点，若重新计算得到的更新几何中心点与原始的几何中心点之间的距离大于设定阈值，则用所述更新几何中心点替换原始的几何中心点，直至两者小于设定阈值后，将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。
本发明实施例提供的陆上无线节点地震仪器桩号匹配方法，利用每个点距中心点的距离与距离平均值做比较，如果大于距离平均值则标记当前点状态为1，否则标记当前点状态为0。利用状态标记为1的点重新计算几何中心点。该方法通过标记点的状态，标记出离散点、飞点的状态为0，将这些点不参与新几何中心点的计算。完全消除了离散点、飞点对计算新几何中心点的影响，克服了反距离权重算法中离散点、飞点对桩号匹配的影响。通过循环迭代计算，几何中心点逐渐向点分布密度大的地方移动。等循环终止，以最后计算的几何中心点位置最优，以该中心点为圆心以特定距离为半径(通常为半个道距)的圆包含的点集数量最多。
在一些实施例中，所述点位包括：坐标东坐标值和北坐标值，所述根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离，包括：
累加每个陆上无线节点的东坐标值，除以点集总数，计算出东坐标的平均值；累加每个陆上无线节点的北坐标值，除以点集总数，计算出北坐标的平均值；
以东坐标和北坐标平均值构造几何中心点；
计算点集中每个点和几何中心点的距离，逐个累加距离值并除以总数，求出点集中每个点距离中心点的距离平均值。
在一些实施例中，所述第一状态点与所述几何中心点的距离大于第一设定距离，所述第二状态点与所述几何中心点的距离小于所述第一设定距离。
在一些实施例中，在根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离之前，还包括：
将大地坐标转换为平面直角坐标。
具体而言，解析陆上无线节点QC文件记录的点位坐标，将大地坐标转换为平面直角坐标。
步骤1：计算点集的几何中心点
累加每个点的东坐标值，除以点集总数，计算出东坐标的平均值；累加每个点的北坐标值，除以点集总数，计算出北坐标的平均值。以东坐标、北坐标平均值构造几何中心点centerPt。
步骤2：计算每个点距离几何中心点的平均距离
计算点集中每个点和几何中心点centerPt的距离dist_i，逐个累加dist_i距离值、除以总数，求出点集中每个点距离中心点的距离平均值。
步骤3：标记点集中点的状态
将点集中每个点和几何中心点的距离dist_i与距离平均值进行比较，如果dist_i小于距离平均值，则标记该点的状态为1，否则标记该点的状态为0。
步骤4：重新计算几何中心点
遍历点集中标记为1的点，累加标记点的东坐标值，除以标记点总数，计算出东坐标的平均值；累加标记点的北坐标值，除以标记点总数，计算出北坐标的平均值。以东坐标平均值和北坐标平均值重新构造几何中心点newCenterPt。
步骤5：计算本次中心点newCenterPt距离上次中心点centerPt之间的距离，如果两个点距离差别大于阈值，将本次中心点newCenterPt的值赋予centerPt，进行循环计算。
如果两个点距离差别小于阈值，则终止循环。
经过多次循环迭代，计算的几何中心点newCenterPt作为陆上无线节点点集的几何中心点。
将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。该方法的流程图如图1所示，陆上无线节点仪器点集和计算的几何中心点如图2至图7所示。
常规的方法是利用陆上无线节点记录的大量点位坐标，通过求点集的东坐标和北坐标的算术平均值当作点集的几何中心，然后将几何中心点与地震测线SPS桩号进行比对，完成桩号匹配流程。
改进的算法是利用反距离权重算法求取点集的几何中心，然后将几何中心点与地震测线SPS桩号进行比对，完成桩号匹配流程。
利用算术平均值求取的几何中心与地震测线SPS桩号进行桩号关联匹配。但该算法抗噪声干扰能力弱，受边缘点、离散点和飞点的影响大。计算的几何中心点与预期不符，偏离地震测线SPS桩号情况严重。造成地震数据处理人工干预工作量大、效率低。
反距离权重算法采用点距离的倒数作为权重，使离中心点距离大的点用于计算时占的权重小，离中心点距离小的点用于计算时占的权重大。该算法通过引入权重降具有良好的抗噪声能力，削弱了边缘点、离散点和飞点对计算几何中心点的影响。计算的几何中心点相比利用算术平均值法得到改善，提升了地震数据处理效率。
反距离权重的算法削弱了边缘点对中心点的影响，但并没有消除边缘点对计算的几何中心点的影响。造成陆上无线节点坐标和SPS桩号坐标距离差别仍然存在，桩号匹配率不能达到100％，需要少量人工进行手动干预，制约着地震数据自动化处理效率的提升。在2020年4月28日至2020年5月8日，把该方法应用到华北物探处某探区进行现场桩号匹配，成功率为100％，具体效果如图8至图10所示。
本发明利用每个点距中心点的距离与距离平均值做比较，如果大于距离平均值则标记当前点状态为1，否则标记当前点状态为0。利用状态标记为1的点重新计算几何中心点。该方法通过标记点的状态，标记出离散点、飞点的状态为0，将这些点不参与新几何中心点的计算。完全消除了离散点、飞点对计算新几何中心点的影响，克服了反距离权重算法中离散点、飞点对桩号匹配的影响。
通过循环迭代计算，几何中心点逐渐向点分布密度大的地方移动。等循环终止，以最后计算的几何中心点位置最优，以该中心点为圆心以特定距离为半径(通常为半个道距)的圆包含的点集数量最多。
图12是本发明又一个实施例提供的实现陆上无线节点地震仪器桩号匹配的支付端的结构示意图，如图12所示，本发明实施例提供的实现陆上无线节点地震仪器桩号匹配的支付端包括：
几何中心点计算模块1，解析陆上无线节点QC文件记录的点位坐标，并根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离；
标记模块2，根据所述平均距离标记点集中点的状态，得到第一状态点和第二状态点；
迭代匹配模块3，根据所述第二状态点重新计算几何中心点，若重新计算得到的更新几何中心点与原始的几何中心点之间的距离大于设定阈值，则用所述更新几何中心点替换原始的几何中心点，直至两者小于设定阈值后，将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。
本发明实施例提供的陆上无线节点地震仪器桩号匹配装置，利用每个点距中心点的距离与距离平均值做比较，如果大于距离平均值则标记当前点状态为1，否则标记当前点状态为0。利用状态标记为1的点重新计算几何中心点。该方法通过标记点的状态，标记出离散点、飞点的状态为0，将这些点不参与新几何中心点的计算。完全消除了离散点、飞点对计算新几何中心点的影响，克服了反距离权重算法中离散点、飞点对桩号匹配的影响。通过循环迭代计算，几何中心点逐渐向点分布密度大的地方移动。等循环终止，以最后计算的几何中心点位置最优，以该中心点为圆心以特定距离为半径(通常为半个道距)的圆包含的点集数量最多。
在上述各实施例的基础上，进一步地，所述点位包括：坐标东坐标值和北坐标值，所述几何中心点计算模块，包括：
累加单元，累加每个陆上无线节点的东坐标值，除以点集总数，计算出东坐标的平均值；累加每个陆上无线节点的北坐标值，除以点集总数，计算出北坐标的平均值；
几何中心点构造单元，以东坐标和北坐标平均值构造几何中心点；
距离平均值计算单元，计算点集中每个点和几何中心点的距离，逐个累加距离值并除以总数，求出点集中每个点距离中心点的距离平均值。
在上述各实施例的基础上，进一步地，所述第一状态点与所述几何中心点的距离大于第一设定距离，所述第二状态点与所述几何中心点的距离小于所述第一设定距离。
在上述各实施例的基础上，还包括：坐标转换模块，将大地坐标转换为平面直角坐标。
图13是本发明又一个实施例提供的电子设备的实体结构示意图，如图13所示，该电子设备可以包括：处理器(processor)1201、通信接口(Communications Interface)1202、存储器(memory)1203和通信总线1204，其中，处理器1201，通信接口1202，存储器1203通过通信总线1204完成相互间的通信。处理器1201可以调用存储器1203中的逻辑指令，以执行如下方法：
S1：解析陆上无线节点QC文件记录的点位坐标，并根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离；
S2：根据所述平均距离标记点集中点的状态，得到第一状态点和第二状态点；
S3：根据所述第二状态点重新计算几何中心点，若重新计算得到的更新几何中心点与原始的几何中心点之间的距离大于设定阈值，则用所述更新几何中心点替换原始的几何中心点，直至两者小于设定阈值后，将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。
此外，上述的存储器1203中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本实施例公开一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法实施例所提供的方法，例如包括：
S1：解析陆上无线节点QC文件记录的点位坐标，并根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离；
S2：根据所述平均距离标记点集中点的状态，得到第一状态点和第二状态点；
S3：根据所述第二状态点重新计算几何中心点，若重新计算得到的更新几何中心点与原始的几何中心点之间的距离大于设定阈值，则用所述更新几何中心点替换原始的几何中心点，直至两者小于设定阈值后，将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。
本实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机程序，所述计算机程序使所述计算机执行上述各方法实施例所提供的方法，例如包括：
S1：解析陆上无线节点QC文件记录的点位坐标，并根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离；
S2：根据所述平均距离标记点集中点的状态，得到第一状态点和第二状态点；
S3：根据所述第二状态点重新计算几何中心点，若重新计算得到的更新几何中心点与原始的几何中心点之间的距离大于设定阈值，则用所述更新几何中心点替换原始的几何中心点，直至两者小于设定阈值后，将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在本说明书的描述中，参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种陆上无线节点地震仪器桩号匹配方法，包括：
解析陆上无线节点QC文件记录的点位坐标，并根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离；
根据所述平均距离标记点集中点的状态，得到第一状态点和第二状态点；
根据所述第二状态点重新计算几何中心点，若重新计算得到的更新几何中心点与原始的几何中心点之间的距离大于设定阈值，则用所述更新几何中心点替换原始的几何中心点，直至两者小于设定阈值后，将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。
2.根据权利要求1所述的陆上无线节点地震仪器桩号匹配方法，其特征在于，所述点位包括：坐标东坐标值和北坐标值，所述根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离，包括：
累加每个陆上无线节点的东坐标值，除以点集总数，计算出东坐标的平均值；累加每个陆上无线节点的北坐标值，除以点集总数，计算出北坐标的平均值；
以东坐标和北坐标平均值构造几何中心点；
计算点集中每个点和几何中心点的距离，逐个累加距离值并除以总数，求出点集中每个点距离中心点的距离平均值。
3.根据权利要求1所述的陆上无线节点地震仪器桩号匹配方法，其特征在于，所述第一状态点与所述几何中心点的距离大于第一设定距离，所述第二状态点与所述几何中心点的距离小于所述第一设定距离。
4.根据权利要求1所述的陆上无线节点地震仪器桩号匹配方法，其特征在于，在根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离之前，还包括：
将大地坐标转换为平面直角坐标。
5.一种陆上无线节点地震仪器桩号匹配装置，包括：
几何中心点计算模块，解析陆上无线节点QC文件记录的点位坐标，并根据所述点位坐标计算点集的几何中心点以及每个点距离几何中心点的平均距离；
标记模块，根据所述平均距离标记点集中点的状态，得到第一状态点和第二状态点；
迭代匹配模块，根据所述第二状态点重新计算几何中心点，若重新计算得到的更新几何中心点与原始的几何中心点之间的距离大于设定阈值，则用所述更新几何中心点替换原始的几何中心点，直至两者小于设定阈值后，将计算出的几何中心点和地震测线SPS桩号进行坐标匹配。
6.根据权利要求5所述的陆上无线节点地震仪器桩号匹配装置，其特征在于，所述点位包括：坐标东坐标值和北坐标值，所述几何中心点计算模块，包括：
累加单元，累加每个陆上无线节点的东坐标值，除以点集总数，计算出东坐标的平均值；累加每个陆上无线节点的北坐标值，除以点集总数，计算出北坐标的平均值；
几何中心点构造单元，以东坐标和北坐标平均值构造几何中心点；
距离平均值计算单元，计算点集中每个点和几何中心点的距离，逐个累加距离值并除以总数，求出点集中每个点距离中心点的距离平均值。
7.根据权利要求5所述的陆上无线节点地震仪器桩号匹配装置，其特征在于，所述第一状态点与所述几何中心点的距离大于第一设定距离，所述第二状态点与所述几何中心点的距离小于所述第一设定距离。
8.根据权利要求5所述的陆上无线节点地震仪器桩号匹配装置，其特征在于，还包括：
坐标转换模块，将大地坐标转换为平面直角坐标。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一项所述方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至4任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011201851.4A CN112379409A (zh) | 2020-11-02 | 2020-11-02 | 一种陆上无线节点地震仪器桩号匹配方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011201851.4A CN112379409A (zh) | 2020-11-02 | 2020-11-02 | 一种陆上无线节点地震仪器桩号匹配方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112379409A true CN112379409A (zh) | 2021-02-19 |
Family
ID=74576573
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011201851.4A Pending CN112379409A (zh) | 2020-11-02 | 2020-11-02 | 一种陆上无线节点地震仪器桩号匹配方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112379409A (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080125161A1 (en) * | 2006-04-19 | 2008-05-29 | Mustafa Ergen | Method and system for hybrid positioning using partial distance information |
| US20100248637A1 (en) * | 2009-03-31 | 2010-09-30 | Zafer Sahinoglu | Method for Estimating Location of Nodes in Wireless Networks |
| US20150219790A1 (en) * | 2012-08-16 | 2015-08-06 | Magseis As | Ocean bottom seismic node system |
| CN107544092A (zh) * | 2017-07-27 | 2018-01-05 | 中国石油天然气集团公司 | 垂直地震剖面采集中的数据处理方法及装置 |
| CN209606627U (zh) * | 2019-01-18 | 2019-11-08 | 潜能恒信能源技术股份有限公司 | 一种节点感应器阵列无人机巡检系统 |
| CN111400851A (zh) * | 2018-12-29 | 2020-07-10 | 中国石油天然气集团有限公司 | 节点布设方法及系统 |
| CN114779299A (zh) * | 2022-04-21 | 2022-07-22 | 中石化石油工程技术服务有限公司 | 一种节点地震仪器精准定位方法 |
-
2020
- 2020-11-02 CN CN202011201851.4A patent/CN112379409A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080125161A1 (en) * | 2006-04-19 | 2008-05-29 | Mustafa Ergen | Method and system for hybrid positioning using partial distance information |
| US20100248637A1 (en) * | 2009-03-31 | 2010-09-30 | Zafer Sahinoglu | Method for Estimating Location of Nodes in Wireless Networks |
| US20150219790A1 (en) * | 2012-08-16 | 2015-08-06 | Magseis As | Ocean bottom seismic node system |
| CN107544092A (zh) * | 2017-07-27 | 2018-01-05 | 中国石油天然气集团公司 | 垂直地震剖面采集中的数据处理方法及装置 |
| CN111400851A (zh) * | 2018-12-29 | 2020-07-10 | 中国石油天然气集团有限公司 | 节点布设方法及系统 |
| CN209606627U (zh) * | 2019-01-18 | 2019-11-08 | 潜能恒信能源技术股份有限公司 | 一种节点感应器阵列无人机巡检系统 |
| CN114779299A (zh) * | 2022-04-21 | 2022-07-22 | 中石化石油工程技术服务有限公司 | 一种节点地震仪器精准定位方法 |
Non-Patent Citations (7)
| Title |
|---|
| ZHIYUAN YIN 等: "Seismic Exploration Wireless Sensor System Based on Wi-Fi and LTE", SENSORS, vol. 20, no. 04 * |
| 严皓 等: "数字化地震队系统(DSS)在国外某混采项目的应用", 石油管材与仪器, vol. 6, no. 01, pages 81 - 84 * |
| 孙洁 等: "Z700节点系统在生产中的应用", 物探装备, vol. 22, no. 03, pages 204 - 209 * |
| 封召鹏 等: "节点地震勘探仪器eSeis 1.0性能探讨", 石油管材与仪器, vol. 4, no. 04, pages 38 - 43 * |
| 曲华东 等: "采集中无法录制或丢失SPS文件的恢复方法", 物探装备, vol. 25, no. 04, 31 August 2015 (2015-08-31), pages 262 - 265 * |
| 郑永明 等: "基于坐标排序的陆上无线节点桩号匹配方法及实现", 中国石油学会2021年物探技术研讨会论文集, pages 1511 - 1514 * |
| 郭飞 等: "GSR滑动扫描中坐标系统问题的解决", 物探装备, vol. 25, no. 01, pages 29 - 32 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102243320A (zh) | 一种地震波初至拾取方法 | |
| CN106650192B (zh) | 一种火山岩型铀矿床磁性界面反演方法 | |
| CN110286416B (zh) | 一种基于物性函数的快速二维密度反演方法 | |
| CN112084672B (zh) | 一种基于分形维数判断地下水污染的方法 | |
| Ni et al. | Lineament length and density analyses based on the segment tracing algorithm: a case study of the gaosong field in gejiu tin mine, China | |
| CN108646292A (zh) | 裂缝密度预测方法、装置及计算机存储介质 | |
| CN106054252B (zh) | 一种叠前时间偏移的方法及装置 | |
| CN112114358B (zh) | 一种基于三维地震资料表征的地下火山通道识别方法 | |
| CN112379409A (zh) | 一种陆上无线节点地震仪器桩号匹配方法及装置 | |
| Yang et al. | A novel algorithm with heuristic information for extracting drainage networks from raster DEMs | |
| Ibrahim et al. | Comparison between inverse distance weighted (IDW) and Kriging | |
| CN112463899A (zh) | 一种车辆轨迹点纠偏方法、系统、服务器和存储介质 | |
| Hou et al. | Automatic Extraction of Drainage Networks from DEMs Base on Heuristic Search. | |
| US11965997B2 (en) | Determining fault surfaces from fault attribute volumes | |
| CN116027454A (zh) | 一种砂岩厚度等值线图的成图方法及智能绘制装置 | |
| CN114779299A (zh) | 一种节点地震仪器精准定位方法 | |
| CN109387872A (zh) | 表面多次波预测方法 | |
| EP2877934A2 (en) | Systems and methods for estimating opportunity in a reservoir system | |
| CN111751886A (zh) | 一种基于微地震监测数据的页岩气藏裂缝建模方法 | |
| King | Sequential Gaussian simulation vs. simulated annealing for locating pockets of high-value commercial trees in Pennsylvania | |
| CN111781661A (zh) | 少井地区沉积微相平面展布预测方法及装置 | |
| CN117870651B (zh) | 基于rtk-slam技术的地图高精度采集方法、存储器及存储介质 | |
| CN102455272A (zh) | 确定流体密度及确定固体骨架密度的方法 | |
| Šešum | Genetic algorithms and smoothing filters in solving the geophysical inversion problem | |
| CN112014875B (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 |