CN109683194A - 自适应平滑半径的基准面确定方法及装置 - Google Patents

自适应平滑半径的基准面确定方法及装置 Download PDF

Info

Publication number
CN109683194A
CN109683194A CN201811532483.4A CN201811532483A CN109683194A CN 109683194 A CN109683194 A CN 109683194A CN 201811532483 A CN201811532483 A CN 201811532483A CN 109683194 A CN109683194 A CN 109683194A
Authority
CN
China
Prior art keywords
smooth
physical point
physical
radius
determining
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
CN201811532483.4A
Other languages
English (en)
Other versions
CN109683194B (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.)
Petrochina Co Ltd
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201811532483.4A priority Critical patent/CN109683194B/zh
Publication of CN109683194A publication Critical patent/CN109683194A/zh
Application granted granted Critical
Publication of CN109683194B publication Critical patent/CN109683194B/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

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)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种自适应平滑半径的基准面确定方法及装置，该方法包括：根据每个物理点，根据每个物理点的高程和预先设置的最大炮检距，确定物理点对应的平滑半径，根据物理点对应的平滑半径，计算物理点对应的平滑高程，根据物理点对应的平滑高程进行计算，得到物理点对应的时间差，根据多个时间差，确定基准面的平滑半径，根据基准面的平滑半径确定基准面。避免了根据得到的平滑半径确定基准面的平滑半径，得到的基准面不一定满足静校准对基准面起伏幅度的要求，因此需要多次计算平滑半径进行确定基准面的平滑半径，花费时间多，导致确定基准面平滑半径的效率低的问题，提高了确定基准面平滑半径的效率。

Description

自适应平滑半径的基准面确定方法及装置
技术领域
本发明涉及石油地震勘侦技术领域，具体而言，涉及一种自适应平滑半径的基准面确定方法及装置。
背景技术
野外地震勘探的基准面是由地表的高程决定，地震勘探过程受近地表起伏变化的影响，每激发一个物理点，地震波开始从该物理点向地下传播，当遇到地层的分界面地震波会产生反射并到达底面，由固定范围内的物理点接收并记录反射的地震波。每个物理点因地面起伏高程不同，使用固定范围内的不同物理点记录地震反射波到达地面的时间不同，会产生不同程度的延迟。因此，需要对延迟时间进行校正，将固定范围内不同的物理点校正在一个特定的面上，这个面称为基准面，基准面的计算在对地震资料的处理是至关重要的。
相关技术中，在确定基准面的过程中，给定一个固定半径，根据固定半径对每个物理点的高程进行计算，得到每个物理点对应的平滑半径，每个物理点对应的平滑半径不一定满足静校准对基准面起伏幅度要求，需要多次计算选择满足静校准对基准面起伏幅度要求的平滑半径，将满足静校准对基准面起伏幅度的平滑半径作为基准面的平滑半径。
但是，在计算基准面平滑半径的过程中，根据给定的根据固定半径对每个物理点的高程进行计算，得到的平滑半径确定的基准面的平滑半径时，得到的基准面不一定满足静校准对基准面起伏幅度的要求，需要多次计算平滑半径进行确定基准面的平滑半径，花费时间多，导致确定基准面的平滑半径效率低。
发明内容
本发明的目的在于，针对上述现有技术中的不足，提供一种自适应平滑半径的基准面确定方法及装置，以解决在多次进行计算确定基准面的平滑半径时，花费时间多，导致确定基准面平滑半径的效率低的问题。
为实现上述目的，本发明实施例采用的技术方案如下：
第一方面，本发明实施例提供了一种自适应平滑半径的基准面确定方法，所述方法包括：
对于每个物理点，根据所述物理点的高程和预先设置的最大炮检距，确定所述物理点对应的平滑半径；
根据所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差，所述物理点对应的平滑高程是根据所述物理点对应的平滑半径计算得到的；
根据多个时间差，确定所述基准面的平滑半径。
进一步地，所述根据所述物理点的高程和预先设置的最大炮检距，确定所述物理点对应的平滑半径，包括：
根据所述物理点的高程和预先设置的最大炮检距进行计算，得到所述物理点的均方差；
根据所述物理点的均方差进行计算，得到所述物理点对应的平滑半径。
进一步地，所述根据所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差，包括：
根据所述物理点对应的平滑半径进行计算，得到所述物理点对应的平滑高程；
根据最大炮检距和所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差。
进一步地，所述根据多个时间差，确定所述基准面的平滑半径，包括：
根据多个时间差对应的参数值，确定最大参数值；
对所述最大参数值和预先设置的时间阈值进行比较，得到比较结果；
若所述比较结果指示所述最大参数值不大于所述时间阈值，则将所述最大参数值所属物理点对应的平滑半径作为所述基准面的平滑半径。
进一步地，所述根据多个时间差，确定所述基准面的平滑半径，还包括：
若所述比较结果指示所述最大参数值大于所述时间阈值，则对每个物理点对应的平滑半径进行调整。
第二方面，本发明实施例还提供了一种自适应平滑半径的基准面确定装置，所述装置包括：
第一确定模块，用于对于每个物理点，根据所述物理点的高程和预先设置的最大炮检距，确定所述物理点对应的平滑半径；
计算模块，用于根据所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差，所述物理点对应的平滑高程是根据所述物理点对应的平滑半径计算得到的；
第二确定模块，用于根据多个时间差，确定所述基准面的平滑半径。
进一步地，所述第一确定模块，具体用于根据所述物理点的高程和预先设置的最大炮检距进行计算，得到所述物理点的均方差；根据所述物理点的均方差进行计算，得到所述物理点对应的平滑半径。
进一步地，所述计算模块，具体用于根据所述物理点对应的平滑半径进行计算，得到所述物理点对应的平滑高程；根据最大炮检距和所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差。
进一步地，所述第二确定模块，具体用于根据多个时间差对应的参数值，确定最大参数值；对所述最大参数值和预先设置的时间阈值进行比较，得到比较结果；若所述比较结果指示所述最大参数值不大于所述时间阈值，则将所述最大参数值所属物理点对应的平滑半径作为所述基准面的平滑半径。
进一步地，所述第二确定模块，还具体用于若所述比较结果指示所述最大参数值大于所述时间阈值，则对每个物理点对应的平滑半径进行调整。
综上所述，本发明实施例提供的自适应平滑半径的基准面确定方法及装置，通过获取物理点的高程和预先设置的最大炮检距，可以得到确定物理点对应的平滑半径，根据物理点对应的平滑半径计算物理点对应的平滑高程，根据物理点对应的平滑高程，得到物理点对应的时间差，根据物理点对应的时间差，确定基准面的平滑半径，进而根据确定的基准面的平滑半径确定基准面，避免了根据得到的平滑半径确定基准面的平滑半径，得到的基准面不一定满足静校准对基准面起伏幅度的要求，花费时间多，导致确定基准面平滑半径的效率低的问题，提高了确定基准面平滑半径的效率。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本发明提供的一种自适应平滑半径的基准面确定方法所涉及的应用场景示意图；
图2为本发明一实施例提供的自适应平滑半径的基准面确定方法的流程示意图；
图3为本发明又一实施例提供的自适应平滑半径的基准面确定方法的流程示意图；
图4为本发明一实施例提供的自适应平滑半径的基准面确定方法的应用场景示意图；
图5为本发明一实施例提供的自适应平滑半径的基准面确定方法的应用场景示意图；
图6为本发明一实施例提供的基准面剖面示意图；
图7为本发明一实施例提供的自适应平滑半径的基准面确定装置的示意图；
图8为本发明一实施例提供的自适应平滑半径的基准面确定装置的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。
图1为本发明提供的一种自适应平滑半径的基准面确定方法所涉及的应用场景示意图；如图1所示，该应用场景包括：多个物理点。
其中，可以将该物理点绘制为由多个物理点组成的物理点排布示意图，根据物理点排布示意图获取每个物理点对应的位置坐标和每个物理点的高程。
具体地，根据物理点绘制物理点排布示意图，确定物理点排布示意图中，各物理点的坐标和各物理点的高程，对于每个物理点根据物理点的高程和预先设置的最大炮检距，计算每一物理点对应的平滑半径。根据每一物理点对应的平滑半径，计算每一物理点对应的平滑高程，根据物理点对应的平滑高程计算，得到物理点对应的时间差，根据每一物理点对应的时间差，确定基准面的平滑半径，进而确定基准面。
野外的地震勘测中，受地表起伏变化的影响，野外的地震勘探不可能在一个水平面进行。同一物理点发出的地震波向地下传播，当遇到地震的分层面产生反射并到达地面时，可以由固定范围内的物理点进行接收，固定范围中的物理点因地表起伏的影响，接收到的地震波的时间会产生一定的延迟。因此需要确定基准面来消除延迟时间，使得对地震数据的处理更加准确。
图2为本发明一实施例提供的自适应平滑半径的基准面确定方法的流程示意图，如图2所示，该方法包括：
步骤201、对于每个物理点，根据物理点的高程和预先设置的最大炮检距，确定物理点对应的平滑半径。
野外的地震勘测中，需要根据选定的物理点绘制物理点排布示意图，可以根据物理点排布示意图中每一物理点对应的位置，获取每一物理点对应的高程，并根据每一物理点的高程计算每一物理点对应的均方差，再根据每一物理点对应的均方差和预先设置的最大炮检距，计算得到物理点排布示意图中每一物理点对应的平滑半径。
具体地，获取物理点排布示意图，在物理点排布示意图中，根据预先设置的最大炮检距确定半径，得到以任一物理点为圆心的圆形区域，可以根据圆内每个物理点的高程计算出圆内全部物理点的平均高程，根据圆内全部物理点的平均高程，计算圆心所属物理点对应的均方差，从而得到物理点排布示意图中每一物理点对应的均方差，最后根据每一物理点对应的均方差，结合预先设置的平滑半径计算公式，计算得到每一物理点对应的平滑半径。
其中，平滑半径计算公式可以为：
其中，物理点排布示意图中任一的物理点的位置用pi(x,y)表示，x和y表示物理点排布示意图中任一物理点的坐标，i是物理点的编号，σi是物理点排布示意图中任一物理点pi的均方差，σmin是物理点排布示意图中所有物理点均方差中的均方差最小值，σmax是物理点排布示意图中所有物理点均方差中均方差的最大值，rmin可以是根据预先设置的最大炮检距确定，rmax可以是预先设置的最大炮检距进行确定。ri是任一物理点pi的平滑半径。
例如：可以获取物理点排布示意图中每一物理点对应的高程，以物理点排布示意图中任一物理点p1为圆心，可以以预先设置的最大炮检距的三分之一为半径，得到一圆形区域，并计算圆内所有物理点的平均高程，根据圆内所有物理点的平均高程计算物理点p1的均方差σ1，从而得到物理点排布示意图中每一物理点对应的均方差，确定物理点排布示意图中所有物理点对应的均方差中的最大值σmax和最小值σmin，rmax可以取预先设置的最大炮检距的三分之一，rmin可以取预先设置的最大炮检距的六分之一，应用以上公式计算物理点p1的平滑半径r1，从而得到物理点排布示意图中每一物理点对应的平滑半径。
步骤202、根据物理点对应的平滑高程进行计算，得到物理点对应的时间差。
其中，物理点对应的平滑高程是根据物理点对应的平滑半径计算得到的。
在得到物理点排布示意图中每一物理点对应的平滑半径后，可以根据得到的每一物理点对应的平滑半径，计算每一物理点对应的平滑高程，并根据每一物理点对应的平滑高程，计算每一物理点对应的时间差。
具体地，得到物理点排布示意图中任一物理点对应的平滑半径后，以任一物理点为圆心，以圆心所属物理点对应的平滑半径为半径，得到一圆形区域，计算圆内全部物理点的平均高程，将全部物理点的平均高程，作为圆心所属物理点的平滑高程，从而依次得到物理点排布示意图中每一物理点对应的平滑高程，在得到物理点排布示意图每一物理点对应的平滑高程后，以预先设置的最大炮检距来确定半径，得到一圆形区域，确定圆内全部物理点中多个平滑高程中的的平滑高程最大值和平滑高程最小值，结合预先设置的速度，根据预先设置的物理点对应时间差计算公式，计算得到每一物理点对应的时间差。
其中，时间差计算公式可以为：
其中，以物理点排布示意图中任一物理点为圆心，以预先设置的最大炮检距来确定半径，得到任一物理点为圆心的圆形区域，为圆内物理点对应平滑高程中的平滑高程最大值，为圆内物理点对应平滑高程中的平滑高程最小值。Vc为预先设置的速度。Δti为物理点排布示意图中任一物理点对应的时间差，结合预先设置的物理点对应时间差计算公式，计算得到每一物理点对应的时间差。
例如：当获取物理点排布示意图中物理点p1的平滑半径后，可以以物理点p1为圆心，以物理点p1对应的平滑半径r1为半径，得到一圆形区域，可以计算圆内全部物理点的平均高程，圆内所有物理点的平均高程作为物理点p1的平滑高程，从而得到物理点排布示意图中全部物理点中每一物理点对应的平滑高程，可以以物理点p1为圆心以预先设置的最大炮检距的二分之一为半径，得到一圆形区域，在圆内物理点对应的平滑高程中确定最大平滑高程和最小平滑高程结合预先设置的物理点时间差计算公式计算得到物理点p1对应的时间差Δt1，最后计算得到每一物理点对应的时间差。
步骤203、根据多个时间差，确定基准面的平滑半径。
在得到物理点排布示意图全部物理点中每一物理点对应的时间差后，在多个时间差中确定最大时间差，根据最大时间差对应的参数值，结合地震数据反射波波形周期确定基准面的平滑半径，通过基准面的平滑半径确定基准面。
具体地，根据得到的多个时间差，在多个时间差中确定每一物理点时间差对应的参数值，确定多个时间差中的最大参数值，预先设置的时间阈值可以由地震数据反射波波形的周期确定，根据最大参数值结合预先设置的时间阈值可以确定基准面的平滑半径。若最大参数值满足下述公式，则可以将最大参数值对应的物理点的平滑半径作为基准面的平滑半径，若最大参数值不满足下述公式，则对每个物理点对应的平滑半径进行调整，重新确定基准面的平滑半径。
其中，确定基准面的平滑半径公式可以为：
其中，上述公式中t为地震数据反射波波形周期，Δtmax为最大参数值，若多个时间差对应的最大参数值满足上述公式，则将满足上述公式中最大参数值对应的物理点的平滑半径作为基准面的平滑半径。若多个时间差中的最大参数值不满足上述公式，则改变物理点对应的平滑半径，重新计算基准面的平滑半径。
例如：在得到物理点排布示意图中每一物理点对应的时间差时，在多个时间差中确定最大时间差对应的最大参数值Δtmax，预先设置的时间阈值可以是地震数据反射波波形周期的二分之一，应用上述公式进行判断，若最大参数值Δtmax不大于地震数据反射波波形周期的二分之一时，则Δtmax对应的物理点的平滑半径作为基准面的平滑半径，通过该基准面的平滑半径确定基准面。若Δtmax大于地震数据反射波波形周期的二分之一时，则改变物理点对应的平滑半径，重新计算基准面的平滑半径。
综上所述，本发明实施例提供的自适应平滑半径的基准面确定方法，通过获取物理点的高程和预先设置的最大炮检距，可以得到确定物理点对应的平滑半径，根据物理点对应的平滑半径计算物理点对应的平滑高程，根据物理点对应的平滑高程，得到物理点对应的时间差，根据物理点对应的时间差，确定基准面的平滑半径，进而根据确定的基准面的平滑半径确定基准面，避免了根据给定的根据固定半径对每个物理点的高程进行计算，得到的平滑半径确定的基准面的平滑半径，确定的基准面不一定满足静校准对基准面起伏幅度的要求，花费时间多，导致确定基准面平滑半径效率低的问题，提高了确定基准面平滑半径的效率。
图3为本发明另一实施例提供的自适应平滑半径的基准面确定方法的流程示意图，如图3所示，该方法包括：
步骤301、根据物理点的高程和预先设置的最大炮检距进行计算，得到物理点的均方差。
根据物理点的高程和预先设置的最大炮检距进行计算物理点的均方差，可以得到物理点的均方差。
具体地，获取物理点排布示意图，在物理点排布示意图中，根据预先设置的最大炮检距确定半径，得到以任一物理点为圆心的圆形区域，计算出圆内全部物理点的平均高程，根据该圆心所在的物理点的高程，结合预先设置的均方差计算公式，对圆心所在的物理点的均方差进行计算，得到圆心所在的物理点的均方差，结合预先设置的均方差计算公式，计算得到每一物理点对应的均方差。
其中，均方差计算公式可以为：
其中：
其中，物理点排布示意图中任一物理点的位置用pi(x,y)表示，x和y表示物理点排布示意图中任一物理点的坐标，任一物理点的高程用Ei表示i是物理位置的编号，σi是物理点排布示意图中任一物理点pi的均方差，μ是以物理点pi为圆心，以预先设置的最大炮检距确定的半径，得到一圆形区域，m为圆内物理点的个数，k取常数值1，计算圆内全部物理点的平均高程，得到圆内全部物理点的平均高程后，根据上式进行计算，得到物理点排布示意图中每个物理点对应的均方差。
例如：请参考图4，在物理点排布示意图中，以圆心物理点p为圆心，可以以预先设置的炮检距的三分之一为半径，得到一圆形区域，计算圆内全部物理点的平均高程，根据圆内全部物理点的平均高程结合物理点p的高程E，应用上述供述进行计算物理点p对应的均方差σ。
步骤302、根据物理点的均方差进行计算，得到物理点对应的平滑半径。根据得到物理点排布示意图中每一物理点对应的均方差，可以在多个均方差中确定最大均方差和最小均方差，结合预先设置的最大炮检距，可以计算物理点排布示意图中全部物理点中每一物理点对应的平滑半径。
具体地，根据物理点排布示意图中每个物理点对应的均方差，确定物理点排布示意图中，全部物理点均方差中的最大均方差和最小均方差，确定物理点排布示意图中任一物理点的对应的高程，根据最大均方差、最小均方差、物理点对应的高程和预先设置的最大炮检距，结合预先设置的平滑半径的计算公式，计算得到每一物理点对应的平滑半径。
其中，平滑半径计算公式可以为：
其中，ri是物理点排布示意图中任一物理点pi对应的平滑半径，σi是物理点排布示意图中任一物理点pi对应的的均方差，σmax是多个均方差中均方差的最大值，σmin是多个均方差中均方差的最小值，rmin和rmax是由预先设置的最大炮检距确定。
例如：计算物理点排布示意图中物理点p对应的平滑半径r时，可以在物理点排布示意图中每一物理点对应的均方差中确定物理点对应的最大均方差σmax和最小均方差σmin，rmax可以取预先设置的最大炮检距的三分之一，rmin可以取预先设置的最小炮检距的六分之一，根据物理点p对应的均方差σ，结合上述公式计算得到物理点p1对应的平滑半径r。
步骤303、根据物理点对应的平滑半径进行计算，得到物理点对应的平滑高程。
具体地，获取物理点排布示意图，以物理点排布示意图中任一物理点为圆心，以圆心所在的物理点对应的平滑半径为半径，得到一圆形区域，根据圆内物理点的高程，计算圆内物理点的平均高程，将圆内物理点的平均高程作为圆心物理点的平滑高程。结合预先设置的平滑高程计算公式，计算得到每一物理点对应的平滑高程。
其中，平滑高程计算公式可以为：
其中，是物理点排布示意图中任一物理点pi对应的平滑高程，以pi为圆心以pi对应的平滑半径为半径，得到一圆形区域，圆内的物理点个数记为n，j为常数值1。
步骤304、根据最大炮检距和物理点对应的平滑高程进行计算，得到物理点对应的时间差。
具体地，以物理点排布示意图中任一物理点为圆心，以最大炮检距确定的半径，得到一圆形区域，在圆内每一物理点对应的平滑高程中，确定圆内物理点对应的最大平滑高程和圆内物理点对应的最小平滑高程。结合预先设置的时间差计算公式，计算得到每一物理点对应的时间差。
其中，时间差计算公式可以为：
其中，Δti为物理点排布示意图任一物理点pi对应的时间差，以物理点排布示意图任一物理点pi为圆心以最大炮检距确定的半径，得到一圆形区域，在圆内每一物理点对应的平滑高程中，将物理点对应的平滑高程最大值记为物理点对应平滑高程度的最小值记为Vc为预先设置的速度。
例如：计算物理点p1对应的时间差时，以物理点p1为圆心，以最大炮检距的二分之一为半径，得到一圆形区域，确定圆内全部物理点中物理点对应的最大平滑高程和物理点对应的最小平滑高程，结合预先设置的速度Vc，根据上式计算物理点p1对应的时间差。
步骤305、根据多个时间差对应的参数值，确定最大参数值。
具体地，得到物理点排布示意图中每一物理点对应的时间差后，根据每一物理点对应的时间差的参数值，确定多个参数值中的最大参数值。以便后续步骤中，对基准面平滑半径的确定，进而确定基准面。
例如：物理点排布示意图中，总共有z个物理点，z个物理点中每一物理点对应一个时间差，根据每一物理点对应的时间差的参数值，确定z个参数值中的最大参数值Δtmax。
步骤306、对最大参数值和预先设置的时间阈值进行比较，得到比较结果。
具体地，在确定最大参数值后，根据最大参数值和预先设置的时间阈值进行比较，预先设置的时间阈值根据地震数据反射波波形周期确定，根据预先设置的判断公式进行比较，判断最大参数值是否满足下述公式。
其中，判断公式可以为：
其中，t为地震数据反射波波形周期，Δtmax为最大参数值。判断最大参数值是否满足下述公式，即就是判断最大参数值所属物理点对应的平滑半径作为基准面的平滑半径确定的初次基准面，是否满足静校正基准面的准则。
需要说明的是，参考图5，预先设置的时间阈值是根据地震数据反射波波形周期进行确定的，可以是地震数据反射波波形周期的二分之一，本发明实施例对此不做限定。
步骤307、若比较结果指示最大参数值不大于时间阈值，则将最大参数值所属物理点对应的平滑半径作为基准面的平滑半径。
具体地，图6为基准面的剖面示意图，根据判断若比较结果不大于时间阈值，则将最大参数所属物理点对应的平滑半径作为基准面的平滑半径，根据基准面的平滑半径，确定基准面。
需要说明的是，若最大参数值不大于时间阈值，则表明最大参数值所属的物理点对应的平滑半径作为基准面的平滑半径，进而确定的初次基准面满足静校准基准面定义的原则，则将初次基准面作为最终确定的基准面。
步骤308、若比较结果指示最大参数值大于时间阈值，则对每个物理点对应的平滑半径进行调整。
具体地，若最大参数值大于时间阈值，则需要物理点排布示意图每一物理点对应的平滑半径，根据调整后的平滑半径，可以确定基准面的平滑半径，进而确定基准面。
需要说明的是，若最大参数值大于时间阈值，则表明计算确定的最大参数值所属物理点对应的平滑半径确定的初次基准面，不满足静校准基准面准则，则需要调整物理点排布示意图中每一物理点对应的平滑半径，平滑半径可以根据物理点排布示意图物理点之间的间距进行调整，其他参数不变。
综上所述，本发明实施例提供的自适应平滑半径的基准面确定方法，通过获取物理点的高程和预先设置的最大炮检距，得到物理点的均方差，根据均方差得到物理点对应的平滑半径，应用物理点对应的平滑半径计算得到物理点的平滑高程，在多个时间差对应的参数值中确定最大参数值，对最大参数值和预先设置的时间阈值进行比较，若最大参数值不大于时间阈值，则将最大参数所属物理点对应的平滑半径作为基准面的平滑半径，确定基准面。若最大参数值大于时间阈值，则调整每个物理点对应的平滑半径，根据调整后的物理点对应的平滑半径，进行重新确定基准面。避免了根据给定的固定半径对每个物理点的高程进行计算，得到的平滑半径确定的基准面的平滑半径，确定的基准面不一定满足静校准对基准面起伏幅度的要求，花费时间多，提高了确定基准面平滑半径的效率。
图7为本发明一实施例提供的自适应平滑半径的基准面确定装置的示意图，如图7所示，该装置具体包括：
第一确定模块701，用于对于每个物理点，根据物理点的高程和预先设置的最大炮检距，确定物理点对应的平滑半径；
计算模块702，用于根据物理点对应的平滑高程进行计算，得到物理点对应的时间差，物理点对应的平滑高程是根据物理点对应的平滑半径计算得到的；
第二确定模块703，用于根据多个时间差，确定基准面的平滑半径。
可选地，第一确定模块701，具体用于根据物理点的高程和预先设置的最大炮检距进行计算，得到物理点的均方差；根据物理点的均方差进行计算，得到物理点对应的平滑半径。
可选地，计算模块702，具体用于根据物理点对应的平滑半径进行计算，得到物理点对应的平滑高程；根据最大炮检距和物理点对应的平滑高程进行计算，得到物理点对应的时间差。
可选地，第二确定模块703，具体用于根据多个时间差对应的参数值，确定最大参数值；对最大参数值和预先设置的时间阈值进行比较，得到比较结果；若比较结果指示最大参数值不大于时间阈值，则将最大参数值所属物理点对应的平滑半径作为基准面的平滑半径。
可选地，第二确定模块703，还具体用于若比较结果指示最大参数值大于时间阈值，则对每个物理点对应的平滑半径进行调整。
综上所述，本发明实施例提供的自适应平滑半径的基准面确定装置，通过获取物理点的高程和预先设置的最大炮检距，可以得到确定物理点对应的平滑半径，根据物理点对应的平滑半径计算物理点对应的平滑高程，根据物理点对应的平滑高程，得到物理点对应的时间差，根据物理点对应的时间差，确定基准面的平滑半径，进而根据确定的基准面的平滑半径确定基准面，避免了根据给定的根据固定半径对每个物理点的高程进行计算，得到的平滑半径不一定满足静校准对基准面起伏幅度的要求，花费时间多，导致确定基准面平滑半径的效率低的问题，提高了确定基准面平滑半径的效率。
上述装置用于执行前述实施例提供的方法，其实现原理和技术效果类似，在此不再赘述。
以上这些模块可以是被配置成实施以上方法的一个或多个集成电路，例如：一个或多个特定集成电路(Application Specific Integrated Circuit，简称ASIC)，或，一个或多个微处理器(digital singnal processor，简称DSP)，或，一个或者多个现场可编程门阵列(Field Programmable Gate Array，简称FPGA)等。再如，当以上某个模块通过处理元件调度程序代码的形式实现时，该处理元件可以是通用处理器，例如中央处理器(CentralProcessing Unit，简称CPU)或其它可以调用程序代码的处理器。再如，这些模块可以集成在一起，以片上系统(system-on-a-chip，简称SOC)的形式实现。
图8为本发明一实施例提供的自适应平滑半径的基准面确定装置的示意图，该装置可以集成于终端设备或者终端设备的芯片，该终端可以是具备基准面确定功能的计算设备。
该装置包括：存储器801、处理器802。
存储器801用于存储程序，处理器802调用存储器801存储的程序，以执行上述方法实施例。具体实现方式和技术效果类似，这里不再赘述。
可选地，本发明还提供一种程序产品，例如计算机可读存储介质，包括程序，该程序在被处理器执行时用于执行上述方法实施例。
在本发明所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元，可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或处理器(英文：processor)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(英文：Read-Only Memory，简称：ROM)、随机存取存储器(英文：Random Access Memory，简称：RAM)、磁碟或者光盘等各种可以存储程序代码的介质。

Claims (10)

1.一种自适应平滑半径的基准面确定方法，其特征在于，应用于基准面，所述基准面包括多个物理点，所述方法包括：
对于每个物理点，根据所述物理点的高程和预先设置的最大炮检距，确定所述物理点对应的平滑半径；
根据所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差，所述物理点对应的平滑高程是根据所述物理点对应的平滑半径计算得到的；
根据多个时间差，确定所述基准面的平滑半径。
2.根据权利要求1所述的自适应平滑半径的基准面确定方法，其特征在于，所述根据所述物理点的高程和预先设置的最大炮检距，确定所述物理点对应的平滑半径，包括：
根据所述物理点的高程和预先设置的最大炮检距进行计算，得到所述物理点的均方差；
根据所述物理点的均方差进行计算，得到所述物理点对应的平滑半径。
3.根据权利要求1所述的自适应平滑半径的基准面确定方法，其特征在于，所述根据所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差，包括：
根据所述物理点对应的平滑半径进行计算，得到所述物理点对应的平滑高程；
根据最大炮检距和所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差。
4.根据权利要求1所述的自适应平滑半径的基准面确定方法，其特征在于，所述根据多个时间差，确定所述基准面的平滑半径，包括：
根据多个时间差对应的参数值，确定最大参数值；
对所述最大参数值和预先设置的时间阈值进行比较，得到比较结果；
若所述比较结果指示所述最大参数值不大于所述时间阈值，则将所述最大参数值所属物理点对应的平滑半径作为所述基准面的平滑半径。
5.根据权利要求4所述的自适应平滑半径的基准面确定方法，其特征在于，所述根据多个时间差，确定所述基准面的平滑半径，还包括：
若所述比较结果指示所述最大参数值大于所述时间阈值，则对每个物理点对应的平滑半径进行调整。
6.一种自适应平滑半径的基准面确定装置，其特征在于，应用于基准面，所述基准面包括多个物理点，所述装置包括：
第一确定模块，用于对于每个物理点，根据所述物理点的高程和预先设置的最大炮检距，确定所述物理点对应的平滑半径；
计算模块，用于根据所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差，所述物理点对应的平滑高程是根据所述物理点对应的平滑半径计算得到的；
第二确定模块，用于根据多个时间差，确定所述基准面的平滑半径。
7.根据权利要求6所述的自适应平滑半径的基准面确定装置，其特征在于，所述第一确定模块，具体用于根据所述物理点的高程和预先设置的最大炮检距进行计算，得到所述物理点的均方差；根据所述物理点的均方差进行计算，得到所述物理点对应的平滑半径。
8.根据权利要求6所述的自适应平滑半径的基准面确定装置，其特征在于，所述计算模块，具体用于根据所述物理点对应的平滑半径进行计算，得到所述物理点对应的平滑高程；根据最大炮检距和所述物理点对应的平滑高程进行计算，得到所述物理点对应的时间差。
9.根据权利要求6所述的自适应平滑半径的基准面确定装置，其特征在于，所述第二确定模块，具体用于根据多个时间差对应的参数值，确定最大参数值；对所述最大参数值和预先设置的时间阈值进行比较，得到比较结果；若所述比较结果指示所述最大参数值不大于所述时间阈值，则将所述最大参数值所属物理点对应的平滑半径作为所述基准面的平滑半径。
10.根据权利要求9所述的自适应平滑半径的基准面确定装置，其特征在于，所述第二确定模块，还具体用于若所述比较结果指示所述最大参数值大于所述时间阈值，则对每个物理点对应的平滑半径进行调整。
CN201811532483.4A 2018-12-13 2018-12-13 自适应平滑半径的基准面确定方法及装置 Active CN109683194B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811532483.4A CN109683194B (zh) 2018-12-13 2018-12-13 自适应平滑半径的基准面确定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811532483.4A CN109683194B (zh) 2018-12-13 2018-12-13 自适应平滑半径的基准面确定方法及装置

Publications (2)

Publication Number Publication Date
CN109683194A true CN109683194A (zh) 2019-04-26
CN109683194B CN109683194B (zh) 2023-10-03

Family

ID=66187518

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811532483.4A Active CN109683194B (zh) 2018-12-13 2018-12-13 自适应平滑半径的基准面确定方法及装置

Country Status (1)

Country Link
CN (1) CN109683194B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1123328A (en) * 1965-09-13 1968-08-14 Gulf Research Development Co Method of removing the effect of multiple and other non-primary reflections in seismic reflection shooting
US6757616B1 (en) * 2003-02-04 2004-06-29 Charles Wayne Emmons Model-free refraction statics solution
CN1580817A (zh) * 2003-08-01 2005-02-16 中国石油天然气集团公司 分步静校正计算方法
CN101980054A (zh) * 2010-10-22 2011-02-23 中国石油化工股份有限公司 一种在高密度地震静校正处理中建立近地表速度模型的方法
CN101984366A (zh) * 2010-09-29 2011-03-09 北京吉星吉达科技有限公司 基于共孔径面的起伏地表叠前时间偏移方法及设备
KR20140105275A (ko) * 2013-02-22 2014-09-01 경북대학교 산학협력단 지표면 높이를 이용한 노풍도 산정 장치 및 방법
CN104459794A (zh) * 2014-12-01 2015-03-25 郑鸿明 共反射点道集时变时间差值的校正方法及装置
CN106932823A (zh) * 2015-12-31 2017-07-07 上海青凤致远地球物理地质勘探科技有限公司 一种光滑基准面选择技术

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1123328A (en) * 1965-09-13 1968-08-14 Gulf Research Development Co Method of removing the effect of multiple and other non-primary reflections in seismic reflection shooting
US6757616B1 (en) * 2003-02-04 2004-06-29 Charles Wayne Emmons Model-free refraction statics solution
CN1580817A (zh) * 2003-08-01 2005-02-16 中国石油天然气集团公司 分步静校正计算方法
CN101984366A (zh) * 2010-09-29 2011-03-09 北京吉星吉达科技有限公司 基于共孔径面的起伏地表叠前时间偏移方法及设备
CN101980054A (zh) * 2010-10-22 2011-02-23 中国石油化工股份有限公司 一种在高密度地震静校正处理中建立近地表速度模型的方法
KR20140105275A (ko) * 2013-02-22 2014-09-01 경북대학교 산학협력단 지표면 높이를 이용한 노풍도 산정 장치 및 방법
CN104459794A (zh) * 2014-12-01 2015-03-25 郑鸿明 共反射点道集时变时间差值的校正方法及装置
CN106932823A (zh) * 2015-12-31 2017-07-07 上海青凤致远地球物理地质勘探科技有限公司 一种光滑基准面选择技术

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郑鸿明 等: "基准面校正的理论研究及误差分析", 新疆地质, vol. 23, no. 1, pages 79 - 81 *

Also Published As

Publication number Publication date
CN109683194B (zh) 2023-10-03

Similar Documents

Publication Publication Date Title
KR102170105B1 (ko) 신경 네트워크 구조의 생성 방법 및 장치, 전자 기기, 저장 매체
US8254723B2 (en) System and method for extracting boundary elements of an object
US10346279B2 (en) Method for optimized canary deployments for improved customer experience
US20150109290A1 (en) Device and method for removing noise points in point clouds
Campen et al. Practical anisotropic geodesy
CN112767551B (zh) 三维模型构建方法、装置以及电子设备、存储介质
Shragge Solving the 3D acoustic wave equation on generalized structured meshes: A finite-difference time-domain approach
EP2601622A1 (en) Predicting a user behavior number of a word
US20160171761A1 (en) Computing device and method for patching point clouds of object
US20220148135A1 (en) Denoising images rendered using monte carlo renderings
JP2017010554A (ja) 曲線検出方法と曲線検出装置
US20200097439A1 (en) Systems and methods for improving the interpretability and transparency of machine learning models
US7813907B2 (en) Hybrid method for enforcing curvature related boundary conditions in solving one-phase fluid flow over a deformable domain
US11530908B2 (en) Measurement point determination method, non-transitory storage medium, and measurement point determination apparatus
CN109683194A (zh) 自适应平滑半径的基准面确定方法及装置
CN111553986B (zh) 三角网的构建方法、构建装置和数字表面模型的生成方法
CN113947668A (zh) 一种模型简化方法及其相关设备
CN111080792B (zh) 模型简化处理方法、装置以及电子设备、存储介质
US10685479B2 (en) Reconstruction of object surfaces from a point cloud by employing regional level sets
CN110866043A (zh) 数据预处理方法及装置、存储介质、终端
US20160232919A1 (en) Estimation of noise characteristics
US20190244403A1 (en) Paint reflection around a curvilinear path
CN115511724A (zh) 三维扫描方法、装置、计算机设备和存储介质
CN114139105A (zh) 一种基于多项式拟合的快速声线追踪算法
CN110207699A (zh) 一种定位方法和装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
TA01 Transfer of patent application right

Effective date of registration: 20201030

Address after: Intercontinental building, 16 ande Road, Dongcheng District, Beijing, 100007

Applicant after: PetroChina Company Limited

Address before: 830000 Institute of Geophysics, No. 397, Beijing North Road, Urumqi, the Xinjiang Uygur Autonomous Region

Applicant before: Zheng Hongming

TA01 Transfer of patent application right
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant