CN112346125B - 一种地震数据处理vda双参数分析方法 - Google Patents

一种地震数据处理vda双参数分析方法 Download PDF

Info

Publication number
CN112346125B
CN112346125B CN202011230836.2A CN202011230836A CN112346125B CN 112346125 B CN112346125 B CN 112346125B CN 202011230836 A CN202011230836 A CN 202011230836A CN 112346125 B CN112346125 B CN 112346125B
Authority
CN
China
Prior art keywords
double
point
parameter analysis
parameter
seismic data
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
CN202011230836.2A
Other languages
English (en)
Other versions
CN112346125A (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 Earthquake Disaster Prevention Center
Original Assignee
China Earthquake Disaster Prevention Center
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 Earthquake Disaster Prevention Center filed Critical China Earthquake Disaster Prevention Center
Priority to CN202011230836.2A priority Critical patent/CN112346125B/zh
Publication of CN112346125A publication Critical patent/CN112346125A/zh
Application granted granted Critical
Publication of CN112346125B publication Critical patent/CN112346125B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • 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
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • 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
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time

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

本发明提出了一种地震数据处理VDA双参数分析方法，包括：布置地震数据采集观测系统，包括设置炮点、检波点的位置以及二者的排列关系，利用地震数据采集观测系统采集地震数据，并设置地震波传播相关参数；计算VDA双参数分析方法的变换方程；通过上述变换方程，计算生成速度和地层倾角双参数分析界面，包括(1)生成双参数谱，用于双参数分析估算速度和地层倾角双参数值；(2)将t转换为t0，实现地震数据成像，获得叠加能量剖面，辅助双参数分析。

Description

一种地震数据处理VDA双参数分析方法
技术领域
本发明涉及地震数据处理技术领域，特别涉及一种地震数据处理VDA双参数分析方法。
背景技术
在地震数据处理过程中，有多种零偏移距成像和速度分析方法。
在几乎所有地震数据处理系统赖以为基础的共中心点(CMP)方法中，CMP速度VCMP其实只是一个信号叠加参数，没有地质意义，这是因为所有实际数据远不是理想的水平层状均匀介质，而只有在水平层状均匀介质的情况下才能得到合理的均方根速度VRMS。反射界面的稍微倾斜或者出现速度非均质性，都会导致VCMP的急剧变化。
还有其它一些得到零偏移距成像剖面的方法，如倾角时差校正(DMO)叠加、共反射面元(CRS)叠加和多聚焦(MF)成像等，这些方法都是经过各自不同的变换进行时距曲线叠加。但是，时距曲线的反射信号属于反射界面上的不同点。
利用不同方法求得的零偏移距成像速度，它们往往都与地下反射界面的倾角和偏移距有关，这就不可避免地导致在叠加时得到的是等效速度，与真正的均方根速度vrms存在差别，有些差别很大。众所周知，根据等效速度只有在一种情况下能够得到正确的层速度，那就是理想水平层状均匀介质。对于具有倾向倾角各异的若干反射界面等复杂地质情况来说，不管用哪种方法，得到的速度的准确性更低，这个速度往往只是一个叠加参数，没有实际地质意义，这势必影响成像的准确性，无法实现真正意义上的共反射点零偏移距成像，获得真正的零偏移距时间剖面。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此，本发明的目的在于提出一种地震数据处理VDA双参数分析方法。
为了实现上述目的，本发明的实施例提供一种地震数据处理VDA双参数分析方法，包括如下步骤：
步骤S1，布置地震数据采集观测系统，包括设置炮点、检波点的位置以及二者的排列关系，利用地震数据采集观测系统采集地震数据，并设置地震波传播相关参数，其中，所述地震波传播相关参数，包括：l为炮点S与检波点R之间的距离，v为地震波传播速度，l0表示反射点O处反射界面的法线在地面出射点C与炮点的距离，β为入射波射线和反射波射线之间夹角的一半，θ为反射点O处反射界面的地层倾角；
步骤S2，通过炮点S、检波点R和反射点O构建一个圆，对于圆上每一个反射点，入射线和反射线之间的夹角2β固定不变，在每个反射点处做角平分线，该角平分线即为反射点处反射界面的法线，所有法线相交于一点，该交点为极点；
相等的角β将圆弧SR分成相等的两段弧，该圆的圆心位于炮点S-检波点R线段的中垂线上；令t表示沿入射波射线的旅行时t1和反射波射线的旅行时t2的总旅行时；其中反射点O处法线与炮检点的连线SR相交于点C，C即为反射点O对应的成像道的地表位置，t0表示中间法向射线OC段的双程旅行时，令
Figure BDA0002765147510000021
表示中间法向射线延长到P点即OP段对应的双程旅行时，v表示地震波速度，计算的VDA双参数分析方法的变换方程为：
Figure BDA0002765147510000022
步骤S3，通过上述变换方程，计算生成速度和地层倾角双参数分析界面，包括(1)生成双参数谱，用于双参数分析估算速度和地层倾角双参数值；(2)将t转换为t0，实现地震数据的展开等时线的相应振幅的绝对值叠加，获得叠加能量剖面，该剖面可以确定主要地层的层数和大致位置，辅助双参数分析。
进一步，所述θ进一步为法向射线与垂向线的夹角。
进一步，所述计算VDA双参数分析方法的变换方程，包括如下步骤：
已知如下数学关系：
Figure BDA0002765147510000031
由圆的性质，有如下关系式：
Figure BDA0002765147510000032
把(2)中的关系代入(1),得到:
Figure BDA0002765147510000033
通过计算可以求取圆的直径d：
d＝l/(2sin2β) (4)
进而可以推出如下关系式，
Figure BDA0002765147510000034
利用公式(1)-(4)，可以推出
Figure BDA0002765147510000035
结合椭圆展开变换方程，有
Figure BDA0002765147510000036
公式(6)和(7)联立,可以推出
Figure BDA0002765147510000037
地层倾角θ与射线参数的关系式为：
Figure BDA0002765147510000038
进而变换得出关于速度v和地层倾角θ的新方程
Figure BDA0002765147510000041
进一步，所述步骤S3，还包括：根据所述VDA双参数分析方法的变换方程，生成综合地层倾角和速度的双参数谱，并估算不同位置、不同成像时间处相应的速度和地层倾角双参数的值，以获得更理想的成像剖面。
根据本发明实施例的地震数据处理VDA双参数分析方法，基于VDA双参数成像原理，估算制作用于双参数分析的叠加能量剖面以及速度和地层倾角双参数谱，通过分析获得参数速度和地层倾角，最终实现双参数共反射点成像。该方法没有地下反射层为水平的假设前提条件，地下反射层可以倾斜或弯曲。该方法可以估算不同位置、不同成像时间处相应的速度和地层倾角双参数的值，考虑此双参数的叠加成像方法就消除了传统单一速度分析方法估算的叠加速度与实际地震波射线速度不同导致的成像计算误差，修正仅使用速度一个参数进行叠加时的不聚焦现象，提高了成像质量。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为根据本发明实施例的地震数据处理VDA双参数分析方法的流程图；
图2为根据本发明实施例的地震波传播路径示意图；
图3为根据本发明实施例的时VDA双参数成像及双参数分析原理图；
图4为根据本发明实施例的VDA双参数分析界面图；
图5为根据本发明实施例的VDA双参数分析获得的地层倾角剖面图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
本发明提供一种地震数据处理VDA双参数分析方法，考虑地层倾角和速度两个参数进行零偏移距成像，修正仅使用速度一个参数进行零偏移距成像时的不聚焦现象，提高成像质量。
如图1所示，本发明实施例的地震数据处理VDA(VelocityandFormationDipAngle)双参数分析方法，包括如下步骤：
步骤S1，布置地震数据采集观测系统，包括设置炮点、检波点的位置以及二者的排列关系，利用地震数据采集观测系统采集地震数据，并设置地震波传播相关参数，其中，所述地震波传播相关参数，包括：l为炮点S与检波点R之间的距离，v为地震波传播速度，l0表示反射点O处反射界面的法线在地面出射点C与炮点的距离，β为入射波射线和反射波射线之间夹角的一半，θ为反射点O处反射界面的地层倾角。θ进一步为法向射线与垂向线的夹角。
反射波速度与入射波速度相等时,地震波传播路径如图2所示。
步骤S2，通过炮点S、检波点R和反射点O构建一个圆，对于圆上每一个反射点，入射线和反射线之间的夹角2β固定不变，在每个反射点处做角平分线，该角平分线即为反射点处反射界面的法线，所有法线相交于一点，该交点为极点。
相等的角β将圆弧SR分成相等的两段弧，该圆的圆心位于炮点S-检波点R线段的中垂线上；令t表示沿入射波射线的旅行时t1和反射波射线的旅行时t2的总旅行时；其中反射点O处法线与炮检点的连线SR相交于点C，C即为反射点O对应的成像道的地表位置，t0表示中间法向射线OC段的双程旅行时，令
Figure BDA0002765147510000051
表示中间法向射线延长到P点即OP段对应的双程旅行时，v表示地震波速度，计算的VDA双参数分析方法的变换方程为：
Figure BDA0002765147510000052
图3展示了VDA双参数分析方法原理。如图3所示，设置炮点S、检波点R，其位置坐标已知。它们对应地下某任意反射地层上的一个反射点O，三者满足地震波反射定律。O处相应地层倾角为θ。下面推导用于VDA双参数分析的变换方程。炮点S、检波点R和反射点O三点确定一个圆。根据反射定律，在O点处作反射界面的法线OC，其在炮检点连线SR上的出露点为C，与圆的交点为P。OC即入射波射线SO和反射波射线OR的夹角的角平分线，将角∠SOR分为入射角和反射角，二者相等，记为β。C点即为反射点O对应的零偏移距成像道的地表位置。
该圆具有很好的几何性质：该圆可以看作是一组可能反射点的轨迹，对于该圆上的每一个反射点，入射线和反射线之间的夹角固定不变；在每个反射点处作角平分线，该角平分线即是各个反射点处反射界面的法线，所有的法线都恰巧相交于点P，P点也是圆和纵向垂直的直径的交点，称之为极点；该圆的圆心位于SR线段的中垂线上。
令t表示地震波沿入射波射线的旅行时t1和反射波射线的旅行时t2的总旅行时。t0表示法向射线OC段的双程旅行时。令
Figure BDA0002765147510000061
表示法向射线延长到P点即OP段对应的双程旅行时，v表示地震波速度。l表示炮点S与检波点R之间的距离，l0表示反射点O处反射界面的法线在SR上的出射点C与炮点S的距离。
计算VDA双参数分析方法的变换方程，包括如下步骤：
已知如下数学关系：
Figure BDA0002765147510000062
由圆的性质，有如下关系式：
Figure BDA0002765147510000063
把(2)中的关系代入(1),得到:
Figure BDA0002765147510000064
通过计算可以求取圆的直径d：
d＝l/(2sin2β) (4)
进而可以推出如下关系式，
Figure BDA0002765147510000071
利用公式(1)-(4)，可以推出
Figure BDA0002765147510000072
结合椭圆展开变换方程，有
Figure BDA0002765147510000073
公式(6)和(7)联立,可以推出
Figure BDA0002765147510000074
地层倾角θ与射线参数的关系式为：
Figure BDA0002765147510000075
进而变换得出关于速度v和地层倾角θ的新方程
Figure BDA0002765147510000076
公式(10)是均匀介质时的VDA双参数分析的变换方程，该方程中含有速度和地层倾角两个参数，对于均匀介质来说，它绝对准确，通过该方程，可以实现速度和地层倾角双参数分析，估算拾取速度和地层倾角，用于零偏移距成像。步骤S3，通过上述变换方程，可以把t转换为t0，实现地震数据的展开等时线的相应振幅的绝对值叠加，获得叠加能量剖面，该剖面可以确定主要地层的层数和大致位置，辅助双参数分析。
具体的，公式(10)是VDA双参数分析方法的变换方程，该方程中含有速度和地层倾角两个参数，绝对准确，通过该方程，可以实现速度和地层倾角双参数分析。
步骤S3，通过上述变换方程，计算生成速度和地层倾角双参数分析界面，包括(1)生成双参数谱，用于双参数分析估算速度和地层倾角双参数值；(2)将t转换为t0，实现地震数据的展开等时线的相应振幅的绝对值叠加，获得叠加能量剖面，该剖面可以确定主要地层的层数和大致位置，辅助双参数分析。
图4和图5为某数据的VDA双参数分析界面。
图4中左图为叠加能量剖面与零偏移距成像剖面的叠合显示，从浅至深，无论是断面、低倾角地层、陡倾角地层、地质异常体和背斜构造，二者吻合一致，各反射层能量集中，干涉噪音小，地层接触关系清晰。右图为CP1380处的双参数谱的时间-速度剖面，用于显示速度变化趋势和选择要人工分析的双参数谱的时间。中间小图为CP1380处t0为1328ms时刻的双参数谱的速度-地层倾角剖面，横向尺度为速度，纵向尺度为地层倾角。通过该谱可以实现双参数分析，估算获得双参数的值。
此外，步骤S3还包括：根据所述VDA双参数分析方法的变换方程，生成综合地层倾角和速度的双参数谱，并估算不同位置、不同成像时间处相应的速度和地层倾角双参数的值，以获得更理想的成像剖面。具体的，图5为本发明自动估算的地层倾角剖面与零偏移距成像剖面的叠合显示。从浅至深，无论是低倾角地层还是陡倾角地层，估算的地层倾角的变化趋势与实际零偏移距成像剖面上反射层的形态基本一致，较好的说明了该方法的有效性。右下角小图为CP1380处t0为1328ms时刻的双参数谱速度-地层倾角剖面，横向尺度为速度，纵向尺度为地层倾角。通过该谱可以实现双参数分析，估算获得不同位置、不同零偏移距成像时间处相应的速度和地层倾角双参数的值，以获得更理想的零偏移距成像剖面。
根据本发明实施例的地震数据处理VDA双参数分析方法，基于VDA双参数成像原理，估算制作用于双参数分析的叠加能量剖面以及速度和地层倾角双参数谱，通过分析获得参数速度和地层倾角，最终实现双参数共反射点成像。该方法没有地下反射层为水平的假设前提条件，地下反射层可以倾斜或弯曲。该方法可以估算不同位置、不同成像时间处相应的速度和地层倾角双参数的值，考虑此双参数的叠加成像方法就消除了传统单一速度分析方法估算的叠加速度与实际地震波射线速度不同导致的成像计算误差，修正仅使用速度一个参数进行叠加时的不聚焦现象，提高了成像质量。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。

Claims (3)

1.一种地震数据处理VDA双参数分析方法，其特征在于，包括如下步骤：
步骤S1，布置地震数据采集观测系统，包括设置炮点、检波点的位置以及二者的排列关系，利用地震数据采集观测系统采集地震数据，并设置地震波传播相关参数，其中，所述地震波传播相关参数，包括：l为炮点S与检波点R之间的距离，v为地震波传播速度，l0表示反射点O处反射界面的法线在地面出射点C与炮点的距离，β为入射波射线和反射波射线之间夹角的一半，θ为反射点O处反射界面的地层倾角；
步骤S2，通过炮点S、检波点R和反射点O构建一个圆，对于圆上每一个反射点，入射线和反射线之间的夹角2β固定不变，在每个反射点处做角平分线，该角平分线即为反射点处反射界面的法线，所有法线相交于一点，该交点为极点；
相等的角β将圆弧SR分成相等的两段弧，该圆的圆心位于炮点S-检波点R线段的中垂线上；令t表示沿入射波射线的旅行时t1和反射波射线的旅行时t2的总旅行时；其中反射点O处法线与炮检点的连线SR相交于点C，C即为反射点O对应的成像道的地表位置，t0表示中间法向射线OC段的双程旅行时，令
Figure FDA0003953245490000011
表示中间法向射线延长到P点即OP段对应的双程旅行时，v表示地震波速度，计算的VDA双参数分析方法的变换方程为：
Figure FDA0003953245490000012
所述计算VDA双参数分析方法的变换方程，包括如下步骤：
已知如下数学关系：
Figure FDA0003953245490000013
由圆的性质，有如下关系式：
Figure FDA0003953245490000021
把(2)中的关系代入(1),得到:
Figure FDA0003953245490000022
通过计算可以求取圆的直径d：
d＝l/(2 sin 2β) (4)
进而可以推出如下关系式，
Figure FDA0003953245490000023
利用公式(1)-(4)，可以推出
Figure FDA0003953245490000024
结合椭圆展开变换方程，有
Figure FDA0003953245490000025
公式(6)和(7)联立,可以推出
Figure FDA0003953245490000026
地层倾角θ与射线参数的关系式为：
Figure FDA0003953245490000027
进而变换得出关于速度v和地层倾角θ的新方程
Figure FDA0003953245490000028
公式(10)是均匀介质时的VDA双参数分析的变换方程，该方程中含有速度和地层倾角两个参数，通过该方程，实现速度和地层倾角双参数分析，估算拾取速度和地层倾角，用于零偏移距成像；
步骤S3，通过上述变换方程，计算生成速度和地层倾角双参数分析界面，包括(1)生成双参数谱，用于双参数分析估算速度和地层倾角双参数值；(2)将t转换为t0，实现地震数据的展开等时线的相应振幅的绝对值叠加，获得叠加能量剖面，该剖面可以确定主要地层的层数和大致位置，辅助双参数分析。
2.如权利要求1所述的地震数据处理VDA双参数分析方法，其特征在于，所述θ进一步为法向射线与垂向线的夹角。
3.如权利要求1所述的地震数据处理VDA双参数分析方法，其特征在于，所述步骤S3，还包括：估算不同位置、不同成像时间处相应的速度和地层倾角双参数的值，以获得更理想的成像剖面。
CN202011230836.2A 2020-11-06 2020-11-06 一种地震数据处理vda双参数分析方法 Active CN112346125B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011230836.2A CN112346125B (zh) 2020-11-06 2020-11-06 一种地震数据处理vda双参数分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011230836.2A CN112346125B (zh) 2020-11-06 2020-11-06 一种地震数据处理vda双参数分析方法

Publications (2)

Publication Number Publication Date
CN112346125A CN112346125A (zh) 2021-02-09
CN112346125B true CN112346125B (zh) 2023-03-14

Family

ID=74429468

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011230836.2A Active CN112346125B (zh) 2020-11-06 2020-11-06 一种地震数据处理vda双参数分析方法

Country Status (1)

Country Link
CN (1) CN112346125B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101285894A (zh) * 2008-05-30 2008-10-15 中国科学院地质与地球物理研究所 起伏地表下采集的地震资料的直接叠前时间偏移方法
CN101900832A (zh) * 2010-06-25 2010-12-01 恒泰艾普石油天然气技术服务股份有限公司 地震资料处理真地表条件下的椭圆展开成像方法和装置
CN104570069A (zh) * 2013-10-12 2015-04-29 中国石油集团东方地球物理勘探有限责任公司 一种从起伏地表直接成像的地震成像方法和装置
CN106526664A (zh) * 2015-09-09 2017-03-22 娄兵 一种时空域几何变换水平叠加地震成像方法
CN106950595A (zh) * 2017-02-09 2017-07-14 恒泰艾普集团股份有限公司 一种角变换成像方法与装置
CN106970417A (zh) * 2017-03-23 2017-07-21 恒泰艾普集团股份有限公司 椭圆展开转换波速度分析方法与系统
CN106990434A (zh) * 2017-03-23 2017-07-28 恒泰艾普集团股份有限公司 椭圆展开转换波成像方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2014200209A1 (en) * 2013-01-15 2014-07-31 Cgg Services Sa Seismic data processing including data-constrained surface-consistent correction

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101285894A (zh) * 2008-05-30 2008-10-15 中国科学院地质与地球物理研究所 起伏地表下采集的地震资料的直接叠前时间偏移方法
CN101900832A (zh) * 2010-06-25 2010-12-01 恒泰艾普石油天然气技术服务股份有限公司 地震资料处理真地表条件下的椭圆展开成像方法和装置
CN104570069A (zh) * 2013-10-12 2015-04-29 中国石油集团东方地球物理勘探有限责任公司 一种从起伏地表直接成像的地震成像方法和装置
CN106526664A (zh) * 2015-09-09 2017-03-22 娄兵 一种时空域几何变换水平叠加地震成像方法
CN106950595A (zh) * 2017-02-09 2017-07-14 恒泰艾普集团股份有限公司 一种角变换成像方法与装置
CN106970417A (zh) * 2017-03-23 2017-07-21 恒泰艾普集团股份有限公司 椭圆展开转换波速度分析方法与系统
CN106990434A (zh) * 2017-03-23 2017-07-28 恒泰艾普集团股份有限公司 椭圆展开转换波成像方法及系统

Also Published As

Publication number Publication date
CN112346125A (zh) 2021-02-09

Similar Documents

Publication Publication Date Title
WO2015196779A1 (zh) 利用叠加速度求取高精度地震波速度的方法
US4736347A (en) Multiple stacking and spatial mapping of seismic data
CN108196305B (zh) 一种山地静校正方法
CN109738945B (zh) 一种利用叠前深度偏移成果直接生成构造图的方法
US5083297A (en) Method of improving the seismic resolution of geologic structures
CN108710148B (zh) 三维倾角域稳相叠前深度偏移方法和装置
CN109444956B (zh) 三维起伏观测面地震斜率层析成像方法
CN104297789A (zh) 一种三维倾角域稳相叠前时间偏移方法及系统
CN102012521A (zh) 一种地震储层预测中叠前裂缝的检测方法
CN104090301B (zh) 一种求取三维高频静校正量的方法
CN104459794A (zh) 共反射点道集时变时间差值的校正方法及装置
CN104570102A (zh) 一种近地表速度模型与中深层速度模型的融合方法
JPH0216486A (ja) 地震データの処理方法
CN107656308B (zh) 一种基于时间深度扫描的共散射点叠前时间偏移成像方法
CN104977615B (zh) 一种基于模型统计拾取的深水obc资料多次波压制方法
CN106950595B (zh) 一种角变换成像方法与装置
CN112346125B (zh) 一种地震数据处理vda双参数分析方法
US5050131A (en) Quantitative method for evaluating velocity model quality
CN112346123B (zh) 一种地震数据处理via双参数分析方法
CN102798888A (zh) 一种利用非零井源距数据计算纵横波速度比的方法
CN112346124B (zh) 一种地震数据处理via双参数成像方法
CN112180452A (zh) 基于探地雷达和三维速度谱的地下管线埋深估计方法
CN111399031A (zh) 一种山地地震资料的采集及处理方法
CN112346122B (zh) 一种地震数据处理vda双参数成像方法
US6131070A (en) Method of generating a fold distribution and of evaluating a seismic survey

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