CN113514887B - 确定黄土山地区地震资料共中心点参考面的方法及装置 - Google Patents

确定黄土山地区地震资料共中心点参考面的方法及装置 Download PDF

Info

Publication number
CN113514887B
CN113514887B CN202010273291.7A CN202010273291A CN113514887B CN 113514887 B CN113514887 B CN 113514887B CN 202010273291 A CN202010273291 A CN 202010273291A CN 113514887 B CN113514887 B CN 113514887B
Authority
CN
China
Prior art keywords
common
center point
correction
point
shot
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
CN202010273291.7A
Other languages
English (en)
Other versions
CN113514887A (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 CN202010273291.7A priority Critical patent/CN113514887B/zh
Publication of CN113514887A publication Critical patent/CN113514887A/zh
Application granted granted Critical
Publication of CN113514887B publication Critical patent/CN113514887B/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. analysis, for interpretation, for correction
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy

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)

Abstract

本发明提供一种确定黄土山地区地震资料共中心点参考面的方法及装置，其中，该方法包括：对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量。本发明通过对炮点野外静校正量和检波点野外静校正量进行插值处理，又通过平滑计算，更精准地确定了共中心点参考面的校正量，以共中心点参考面的校正量的形式代表共中心点参考面，与现有技术相比，提高了黄土山地区地震资料共中心点参考面的求取精度。

Description

确定黄土山地区地震资料共中心点参考面的方法及装置
技术领域
本发明涉及石油地球物理勘探处理技术领域，特别涉及一种确定黄土山地区地震资料共中心点参考面的方法及装置。
背景技术
随着油田加大勘探开发力度，油田对黄土山地区地震资料成像要求越来越高，而黄土山地区地震资料成像很大程度取决于CMP(Common Middle Point，即共中心点)参考面。
CMP参考面在地震数据应用时是以一种CMP校正量的形式存在的，如图1所示。黄土山地地震资料，尤其是巨厚黄土山地地震资料，静校正问题严重，反射波双曲线扭曲严重，无法利用速度分析求取速度值。反射波速度分析主要是在应用完炮点高频静校正量和检波点高频静校正量后的CMP道集上进行的，原始CMP道集数据被校正到CMP面上，反射波同相轴双曲线特征明显，这个速度分析的时间面为CMP参考面。理论上讲，没有经过时间上漂移的原始CMP道集求取的速度才是准确的，但是应用完炮点和检波点到CMP参考面的校正量，如图1所示，原始CMP道集数据被校正到CMP面上，反射波双曲线的t0(时距曲线数学方程的一个参数,称为自激自收时间)时就发生了时间上的垂直移动，即时间漂移。而黄土山地地震资料，尤其是巨厚黄土山地地震资料，高程起伏变化大，炮点和检波点野外校正量变化剧烈，而目前工业方法确定的CMP参考面与对应的原始地表面差距很大。
常用的CMP参考面校正量的确定方法，是通过计算同一CMP道集内所有地震道的炮点野外静校正量和检波点野外静校正量之和的平均值，这种平均效应导致CMP参考面校正量变化趋势与炮检点校正量变化趋势不一致，从而造成CMP参考面道集t0时间漂移量巨大。因此目前技术中难以精准确定黄土山地区地震资料的CMP参考面。
发明内容
本发明实施例提供了一种确定黄土山地区地震资料共中心点参考面的方法，用以精准确定黄土山地区地震资料的CMP参考面，该方法包括：
对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；
对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；
对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量。
本发明实施例还提供了一种确定黄土山地区地震资料共中心点参考面的装置，用以精准确定黄土山地区地震资料的CMP参考面，该装置包括：
插值计算模块，用于对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；
校正量求和模块，用于对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；
平滑计算模块，用于对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述所述方法。
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述所述方法的计算机程序。
本发明实施例中，对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量，可更精准地确定共中心点参考面的校正量，以共中心点参考面的校正量的形式代表共中心点参考面，与现有技术相比，提高了黄土山地区地震资料共中心点参考面的求取精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的CMP参考面与炮点和检波点野外静校正量示意图；
图2是本发明实施例提供的一种确定黄土山地区地震资料共中心点参考面的方法的流程示意图；
图3是本发明实施例提供的反射波时距曲线示意图；
图4是采用本发明实施例提供的方法得到的共中心点参考面的校正量与采用现有技术得到的共中心点参考面的校正量的对比示意图；
图5a是本发明实施例提供的黄土山地地震资料CMP点原始道集的示意图；
图5b是采用本发明实施例提供的方法得到的炮检点共中心点参考面的高频静校正量的示意图；
图5c是本发明实施例提供的采用现有技术得到的炮检点共中心点参考面的高频静校正量的示意图；
图6a是采用本发明实施例提供的方法得到的CMP参考面的叠加剖面的示意图；
图6b是本发明实施例提供的采用现有技术得到的CMP参考面的叠加剖面的示意图；
图7是本发明实施例提供的一种确定黄土山地区地震资料共中心点参考面的装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前工业方法确定的CMP参考面与对应的原始地表面差距很大，常用的CMP参考面校正量的确定方法，是通过计算同一CMP道集内所有地震道的炮点野外静校正量和检波点野外静校正量之和的平均值，这种平均效应导致CMP参考面校正量变化趋势与炮检点校正量变化趋势不一致，从而造成CMP参考面道集t0时间漂移量巨大。因此目前技术中难以精准确定黄土山地区地震资料的CMP参考面。而本发明实施例提供了一种确定黄土山地区地震资料共中心点参考面的方法，用以精准确定黄土山地区地震资料的CMP参考面，如图2所示，该方法包括：
步骤201：对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；
步骤202：对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；
步骤203：对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量。
本发明实施例中，对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量，可更精准地确定共中心点参考面的校正量，以共中心点参考面的校正量的形式代表共中心点参考面，与现有技术相比，提高了黄土山地区地震资料共中心点参考面的求取精度。
具体实施时，首先对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量。
在一实施例中，通过炮点野外静校正量和检波点野外静校正量，利用数学插值计算，分别获得CMP点位置参考面的炮点校正量和检波点校正量。
在获得共中心点位置的炮点校正量和检波点校正量之后，对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和。
在一实施例中，按如下公式，对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和：
Tj＝TSj0+TRj0
其中，Tj为共中心点位置零炮检距校正量和，TSj0为共中心点位置的炮点校正量，TRj0为共中心点位置的检波点校正量。
在得到共中心点位置零炮检距校正量和后，对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量。
在一实施例中，平滑计算可以包括：如平均值平滑，线性拟合平滑，最小二乘平滑等。
具体实施时，本发明实施例提供的确定黄土山地区地震资料共中心点参考面的方法，还可以包括：在得到共中心点参考面的校正量后，将共中心点参考面的校正量加载到数据道头中，得到炮检点共中心点参考面的高频静校正量。
而在现有技术中，黄土山地区地震资料成像，很大程度取决于速度分析的成败。地震资料数据处理主要是在CMP(Common Middle Point，即共中心点)参考面，通过速度分析，来求取叠加速度。CMP参考面在地震数据应用时，是以一种CMP校正量的形势存在的，如图1所示。黄土山地地震资料，尤其是巨厚黄土山地地震资料，静校正问题严重，反射波双曲线扭曲严重，无法利用速度分析求取速度值。反射波速度分析主要是在应用完炮点高频静校正量和检波点高频静校正量后的CMP道集上进行的，原始CMP道集数据被校正到CMP面上，反射波同相轴双曲线特征明显，这个速度分析的时间面为CMP参考面。理论上讲，没有经过时间上漂移的原始CMP道集求取的速度才是准确的，但是应用完炮点和检波点到CMP参考面的校正量(如图1所示，以下简写为炮点CMP高频量和检波点CMP高频量)后，原始CMP道集数据被校正到CMP面上，反射波双曲线的t0(时距曲线数学方程的一个参数,称为自激自收时间)时就发生了时间上的垂直移动，即时间漂移。而黄土山地地震资料，尤其是巨厚黄土山地地震资料，高程起伏变化大，炮点和检波点野外校正量变化剧烈，而目前常用的CMP参考面校正量的确定方法，是通过计算同一CMP道集内所有地震道的炮点野外静校正量和检波点野外静校正量之和的平均值，这种平均效应导致CMP参考面校正量变化趋势与炮检点校正量变化趋势不一致，从而造成CMP参考面道集t0时间漂移量巨大，从而对速度分析求取的精度影响严重。
根据表1列举的示例，t0时间的改变对速度分析求取的精度影响严重。如何确定黄土山地地震资料的CMP参考面，尤其是巨厚黄土山地区的CMP参考面，使应用完炮检点高频量后CMP道集的t0时不漂移或者漂移量很小，意义重大。图3为反射波时距曲线，其数学表达式为其中x为偏移距或炮检距，vrms为均方根速度值，t0为自激自收时间，t为双程旅行时，图3中tx为与炮点S的距离x的检波点接收到地震数据的双程旅行时，ΔtNMO为tx与t0的差值，即动校正时差，S代表炮点。如表1中所示：假设原始地表采集的反射波，t0为0.8s，tx为1.2s，偏移距x为2000m。根据反射波时距曲线公式，速度分析求得速度值为2235.1m/s。而当加入野外静校正量，地震记录做垂向上的常时漂移：
①情况1：假设t0时向下漂移0.02s，由于波形无变化，所以tx也下移动0.02s，从而t0为0.82s，tx为1.22s，这样求得速度值为2214.6.0m/s。速度变小；②情况2：假设t0时向下漂移0.1s，由于波形不变，所以tx也下移动0.1s，从而t0为0.9s，tx为1.3，这样求取的速度值为2132m/s，与真实速度比减小了103m/s；③情况3：假设t0时向上漂移0.1s，由于波形不变，所以tx也上移动0.1s，从而t0为0.7s，tx为1.1s，这样求取的速度值为2357m/s，与真实速度比增大了122m/s。从而可得出，静校正量导致t0向下漂移，求得速度值变小；静校正量导致t0向上漂移，求得速度值偏大；t0漂移量越大，速度值偏差越大。利用反射波求取的均方根速度受CMP道集t0漂移影响。
表1 t0时变化对求取均方根速度的影响
时移 t0 tx x vrms
0 0.8s 1.2s 2000m 2235.1m/s
+0.02 0.82 1.22 2000m 2214.6m/s
+0.1s 0.9s 1.3s 2000m 2132.0m/s
-0.1s 0.7s 1.1s 2000m 2357.0m/s
综上，现有技术中，因黄土山地地震资料，尤其是巨厚黄土山地区地震资料的静校正问题严重，导致反射波双曲线扭曲严重，因此无法进行速度拾取求取速度值。工业生产上利用应用完炮点高频静校正量和检波点高频静校正量后的CMP参考面的道集进行速度分析，而工业方法确定的CMP参考面与对应的原始地表面差距很大，导致应用炮检点高频静校正量的CMP道集的t0时间偏离原始t0时间很大，从而导致速度求取的不准确，降低了反射波速度分析的速度求取的精度。
进一步地，本发明实施例提供的确定黄土山地区地震资料共中心点参考面的方法，还可以包括：在得到炮检点共中心点参考面的高频静校正量后，将炮检点共中心点参考面的高频静校正量加载至共中心点道集，获得共中心点参考面的地震数据。
在本实施例中，在得到共中心点参考面的校正量后，利用GeoEast软件headmath(道头运算)模块加载到数据道头中，求取炮点和检波点的CMP面高频静校正量进行应用。将CMP参考面的校正量加载到处理系统相应地震数据的CMP参考面校正量CMP datum static道头字中，进行相应的炮点高频量Source CMP static和检波点高频量Receiver CMPstatic道头字的计算，得到炮检点共中心点参考面的高频静校正量。之后将炮点静校正量Source CMP static和检波点静校正量Receiver CMP static应用到CMP道集上，获得CMP参考面地震数据，进行后续处理。
与以往方法比，本发明实施例得到的CMP参考面地震数据道集的高频静校正量更接近野外原始道集t0时间。以应用本发明实施例的方法确定的CMP参考面进行速度分析，得到的叠加剖面与以往方法确定的CMP参考进行速度分析得到的叠加剖面相比，用本发明实施例的方法确定的CMP参考面进行速度分析，得到的叠加剖面的成像品质改善明显。
具体实施时，所述对共中心点位置零炮检距校正量和进行平滑计算，包括：根据预设的平滑计算参数，对共中心点位置零炮检距校正量和进行平滑计算；所述平滑计算参数根据黄土山地区的地表起伏情况设置。
在一实施例中，根据预设的平滑计算参数，对共中心点位置零炮检距校正量和进行平滑计算，若黄土山地区的地表起伏严重，可设置该平滑参数进行相应的减小。
本发明实施例还提供了一具体实例，在一黄土山地区如银洞子地区，银洞子地区属于典型的巨厚黄土山地地貌，地表条件复杂，高程变化剧烈。
针对银洞子地区，采用本发明实施例提供的确定黄土山地区地震资料共中心点参考面的方法确定的CMP参考面，如图4所示，采用本发明实施例提供的方法得到的共中心点参考面的校正量(图4中黑色较粗线)，与采用现有技术得到的共中心点参考面的校正量(图4中黑色最粗线)进行对比，可知，本发明实施例提供的方法确定的CMP参考面与炮点校正量线(图4中灰色细线)检波点校正量线(图4中黑色细线)的趋势一致，而以往方法与炮检点校正量差距较大。由于CMP参考面校正量为炮检点校正量之和，故图4中显示的参考面校正量为CMP(共中心点)参考面校正量的一半，主要是为了与炮检点校正量作对比，比较异同。
采用本发明实施例提供的确定黄土山地区地震资料共中心点参考面的方法确定的CMP参考面的炮检点高频量的CMP道集，和应用以往方法CMP参考面的炮检点高频量的CMP道集对比，对比结果如图5a，图5b以及图5c所示，图5a是本发明实施例提供的黄土山地地震资料CMP点原始地震数据道集的示意图；图5b是采用本发明实施例提供的方法得到的炮检点共中心点参考面的高频静校正量应用后同一CMP点地震数据道集的示意图；图5c是本发明实施例提供的采用现有技术得到的炮检点共中心点参考面的高频静校正量应用后同一CMP点地震数据道集的示意图。图5a，图5b以及图5c的横坐标均为偏移距x(单位m)，纵坐标均为双程旅行时t(单位ms)。通过读取道头，本发明实施例提供的方法确定的CMP道集的第一道的炮点高频量检波点高频量/>可计算出第一道的总的高频校正量为/>可得出t0向上移动33.99ms；而以往方法计算的CMP道集的第一道的炮点高频量/>检波点高频量/>可计算出第一道的总的高频校正量为/>可得出t0向下移动112.949ms。因此与以往方法相比较，本发明实施例提供的方法得到的CMP参考面道集，更接近野外原始道集t0时间。
对采用本发明实施例提供的确定黄土山地区地震资料共中心点参考面的方法确定的CMP参考面进行速度分析得到的叠加剖面，与以往方法确定的CMP参考进行速度分析得到的叠加剖面相比，成像品质改善明显，如图6a和图6b所示，图6a是采用本发明实施例提供的方法得到的CMP参考面的叠加剖面的示意图；图6b是本发明实施例提供的采用现有技术得到的CMP参考面的叠加剖面的示意图。图6a和图6b的横坐标均为位置，纵坐标均为双程旅行时t(单位ms)。通过图6a与图6b对比，可以看到图6b在300ms位置的同相轴更加连续、清晰；cmp范围1900到2150范围的400ms到700ms高陡构造同相轴出现，信噪比改善；cmp范围2200到3000范围600ms到1400ms同相轴连续性改善明显；cmp范围1350到1800范围的1300ms位置同相轴更加连续，清晰；cmp范围1100到1600范围1600ms到2200ms位置的绕射波更加清晰连续，双曲线特性明显。
本发明实施例中，对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量，可更精准地确定了共中心点参考面的校正量，以共中心点参考面的校正量的形式代表共中心点参考面，与现有技术相比，提高了黄土山地区地震资料共中心点参考面的求取精度。
实施例中，基于黄土山地区地震资料，首先分析以往地震数据处理CMP参考面确定的方法，然后分析t0时间发生改变的原因，结合CMP参考面校正量计算原理，提出了一种确定黄土山地区地震资料共中心点参考面的方法，能够有效减小t0改变量，从而最大程度的保持原始CMP道集双曲线的形态，提高速度分析求取速度的准确度，改善叠加剖面成像效果。
基于同一发明构思，本发明实施例中还提供了一种确定黄土山地区地震资料共中心点参考面的装置，如下面的实施例所述。由于确定黄土山地区地震资料共中心点参考面的装置解决问题的原理与确定黄土山地区地震资料共中心点参考面的方法相似，因此确定黄土山地区地震资料共中心点参考面的装置的实施可以参见确定黄土山地区地震资料共中心点参考面的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图7是本发明实施例中确定黄土山地区地震资料共中心点参考面的装置的结构示意图，包括：
插值计算模块01，用于对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；
校正量求和模块02，用于对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量；
平滑计算模块03，用于对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量。
在一个实施例中，还包括，校正量计算模块，用于：将共中心点参考面的校正量加载到数据道头中，得到炮检点共中心点参考面的高频静校正量。
在一个实施例中，还包括，地震数据确定模块，用于：将炮检点共中心点参考面的高频静校正量加载至共中心点道集，获得共中心点参考面的地震数据。
在一个实施例中，校正量求和模块，具体用于：按如下公式，对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和：
Tj＝TSj0+TRj0
其中，Tj为共中心点位置零炮检距校正量和，TSj0为共中心点位置的炮点校正量，TRj0为共中心点位置的检波点校正量。
在一个实施例中，平滑计算模块，具体用于：根据预设的平滑计算参数，对共中心点位置零炮检距校正量和进行平滑计算；所述平滑计算参数根据黄土山地区的地表起伏情况设置。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述所述方法。
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述所述方法的计算机程序。
综上所述，本发明实施例中，对炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量，可更精准地确定共中心点参考面的校正量，以共中心点参考面的校正量的形式代表共中心点参考面，与现有技术相比，提高了黄土山地区地震资料共中心点参考面的求取精度。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种确定黄土山地区地震资料共中心点参考面的方法，其特征在于，包括：
对黄土山地区地震资料中的炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；
对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；
对共中心点位置零炮检距校正量和进行平滑计算，得到黄土山地区地震资料的共中心点参考面的校正量；
所述对共中心点位置零炮检距校正量和进行平滑计算，包括：
根据预设的平滑计算参数，对共中心点位置零炮检距校正量和进行平滑计算；所述平滑计算参数根据黄土山地区的地表起伏情况设置；
还包括，在得到共中心点参考面的校正量后，将共中心点参考面的校正量加载到数据道头中，得到炮检点共中心点参考面的高频静校正量。
2.如权利要求1所述的方法，其特征在于，还包括，将炮检点共中心点参考面的高频静校正量加载至共中心点道集，获得共中心点参考面的地震数据。
3.如权利要求1所述的方法，其特征在于，按如下公式，对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和：
Tj＝TSj0+TRj0
其中，Tj为共中心点位置零炮检距校正量和，TSj0为共中心点位置的炮点校正量，TRj0为共中心点位置的检波点校正量。
4.一种确定黄土山地区地震资料共中心点参考面的装置，其特征在于，包括：
插值计算模块，用于对黄土山地区地震资料中的炮点野外静校正量和检波点野外静校正量，进行插值计算，得到共中心点位置的炮点校正量和检波点校正量；
校正量求和模块，用于对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和；
平滑计算模块，用于对共中心点位置零炮检距校正量和进行平滑计算，得到共中心点参考面的校正量；
平滑计算模块，具体用于：根据预设的平滑计算参数，对共中心点位置零炮检距校正量和进行平滑计算；所述平滑计算参数根据黄土山地区的地表起伏情况设置；
还包括，校正量计算模块，用于：将共中心点参考面的校正量加载到数据道头中，得到炮检点共中心点参考面的高频静校正量。
5.如权利要求4所述的装置，其特征在于，还包括，地震数据确定模块，用于：将炮检点共中心点参考面的高频静校正量加载至共中心点道集，获得共中心点参考面的地震数据。
6.如权利要求4所述的装置，其特征在于，校正量求和模块，具体用于：按如下公式，对共中心点位置的炮点校正量和检波点校正量求和，得到共中心点位置零炮检距校正量和：
Tj＝TSj0+TRj0
其中，Tj为共中心点位置零炮检距校正量和，TSj0为共中心点位置的炮点校正量，TRj0为共中心点位置的检波点校正量。
7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至3任一所述方法。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至3任一所述方法的计算机程序。
CN202010273291.7A 2020-04-09 2020-04-09 确定黄土山地区地震资料共中心点参考面的方法及装置 Active CN113514887B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010273291.7A CN113514887B (zh) 2020-04-09 2020-04-09 确定黄土山地区地震资料共中心点参考面的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010273291.7A CN113514887B (zh) 2020-04-09 2020-04-09 确定黄土山地区地震资料共中心点参考面的方法及装置

Publications (2)

Publication Number Publication Date
CN113514887A CN113514887A (zh) 2021-10-19
CN113514887B true CN113514887B (zh) 2023-09-26

Family

ID=78060712

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010273291.7A Active CN113514887B (zh) 2020-04-09 2020-04-09 确定黄土山地区地震资料共中心点参考面的方法及装置

Country Status (1)

Country Link
CN (1) CN113514887B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5508914A (en) * 1994-09-28 1996-04-16 Western Atlas International, Inc. Method for calculating static corrections for seismic data
CN103412333A (zh) * 2013-08-26 2013-11-27 郑鸿明 一种静校正基准面确定方法
CN105093319A (zh) * 2014-05-14 2015-11-25 中国石油化工股份有限公司 基于三维地震数据的地面微地震静校正方法
CN107831536A (zh) * 2017-11-01 2018-03-23 中国石油化工股份有限公司 沙漠地表地震资料处理浮动基准面建立方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050256648A1 (en) * 2004-05-11 2005-11-17 West Michael P Velocity determination of the near-surface layers in the earth using exploration 2D or 3D seismic data
US10386519B2 (en) * 2015-12-18 2019-08-20 Saudi Arabian Oil Company Automated near surface analysis by surface-consistent refraction methods

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5508914A (en) * 1994-09-28 1996-04-16 Western Atlas International, Inc. Method for calculating static corrections for seismic data
CN103412333A (zh) * 2013-08-26 2013-11-27 郑鸿明 一种静校正基准面确定方法
CN105093319A (zh) * 2014-05-14 2015-11-25 中国石油化工股份有限公司 基于三维地震数据的地面微地震静校正方法
CN107831536A (zh) * 2017-11-01 2018-03-23 中国石油化工股份有限公司 沙漠地表地震资料处理浮动基准面建立方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CMP-based near surface velocity model building for mountainous area;Wen Peng et al.;《SEG International Exposition and 88th Annual Meeting》;20181231;第2707-2710页 *
夏常亮等.一种新静校正方法及其在黄土塬地震资料的应用.地球物理学进展.2017,32(5),第2116页. *
赵忠华等.针对浅层气藏的地震处理技术研究与应用.2015 年物探技术研讨会.2015,第234页. *
针对浅层气藏的地震处理技术研究与应用;赵忠华等;《2015年物探技术研讨会》;20151231;第234页 *

Also Published As

Publication number Publication date
CN113514887A (zh) 2021-10-19

Similar Documents

Publication Publication Date Title
Lehmann et al. Topographic migration of georadar data: Implications for acquisition and processing
US20130060475A1 (en) Estimating Reservoir Properties from 4D Seismic Data
CN105093318B (zh) 一种自适应波动方程波场延拓静校正方法
CN108957553B (zh) 动校正量递推修正的无拉伸畸变动校正方法及装置
CN104597485B (zh) 一种微小断层检测方法及断层检测装置
CN113514887B (zh) 确定黄土山地区地震资料共中心点参考面的方法及装置
CN102162858A (zh) 利用非对称走时进行动校正速度分析的方法
CN106932823A (zh) 一种光滑基准面选择技术
CN109581494B (zh) 叠前偏移方法及装置
Shragge et al. Wave-equation migration from topography
AU739128B2 (en) A method of seismic processing, and in particular a 3D seismic prospection method implementing seismic data migration
CN112731520A (zh) 一种基于结构张量扩散滤波的全波形反演方法及系统
CN112213784B (zh) 复杂地表地震数据一次处理快速静校正方法
OA10802A (en) Improved method for migration before summation
RU2155972C2 (ru) Способ определения скоростей миграций при обработке сейсмических данных
WO2022250801A1 (en) Seismic first breaks onset times determination by interface tracking using level-sets method
CN106324677B (zh) 一种位移场剩余动校正方法及装置
CN107179548A (zh) 一种基于真地表的叠前地震成像方法
CN114415234A (zh) 基于主动源面波频散和h/v确定浅地表横波速度的方法
CN113900141A (zh) 油气分布预测方法及装置
NO20190324A1 (en) Automated mistie analysis and correction across two-dimensional ("2d") seismic surveys
CN113267808B (zh) 振幅补偿方法及装置
CN111435172A (zh) 层析静校正的方法和装置
CN113495296B (zh) 层析静校正量的确定方法、装置、设备及可读存储介质
CN111999767B (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
GR01 Patent grant