CN113917536A - 积分法成像实现ovg道集直接输出的方法 - Google Patents

积分法成像实现ovg道集直接输出的方法 Download PDF

Info

Publication number
CN113917536A
CN113917536A CN202010661120.1A CN202010661120A CN113917536A CN 113917536 A CN113917536 A CN 113917536A CN 202010661120 A CN202010661120 A CN 202010661120A CN 113917536 A CN113917536 A CN 113917536A
Authority
CN
China
Prior art keywords
imaging
ovg
gather
offset
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.)
Granted
Application number
CN202010661120.1A
Other languages
English (en)
Other versions
CN113917536B (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 Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
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 Petroleum and Chemical Corp, Geophysical Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN202010661120.1A priority Critical patent/CN113917536B/zh
Publication of CN113917536A publication Critical patent/CN113917536A/zh
Application granted granted Critical
Publication of CN113917536B publication Critical patent/CN113917536B/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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G01V2210/512—Pre-stack

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)
  • Image Processing (AREA)

Abstract

本发明提供一种积分法成像实现OVG道集直接输出的方法，该积分法成像实现OVG道集直接输出的方法包括：步骤1，获取偏移计算的参数；步骤2，根据成像目标线位置及偏移孔径，获取所需的预处理后CMP道集；步骤3，根据参与计算的进程数目，将获取的CMP道集数据均分到各个进程；步骤4，读取偏移所需的速度场；步骤5，计算该目标线成像所需的旅行时场；步骤6，计算当前输入地震道对应的OVG成像输出参数信息；步骤7，输出该成像目标线的OVG道集；步骤8，循环计算下一成像目标线。该积分法成像实现OVG道集直接输出的方法大大降低了成像前和成像后地震数据的读写量和磁盘占有量，对于高密度地震数据具有非常好的适应性。

Description

积分法成像实现OVG道集直接输出的方法
技术领域
本发明涉及地震勘探资料处理与解释技术领域，特别是涉及到一种积分法成像实现OVG道集直接输出的方法。
背景技术
Kirchhoff积分偏移成像方法是目前地震勘探中最主流的成像技术。Kirchhoff积分法叠前偏移通常输出偏移距域成像道集，并已成为地震资料处理中必备的常规处理流程。随着高密度地震采集的出现，Kirchhoff积分偏移生成的OVG(Offset Vector Gather)成像道集，能够更好地用于层析速度反演，提高偏移速度模型和地震成像质量；用于叠前属性反演与分析，提高油气预测精度。但是，在当前计算机条件下输出OVG道集却仍有较大压力。这是因为输出OVG道集需要输出I(nz,ncxline,ninline,nhx,nhy)的一个五维成像道集，其中nz为时间或深度方向采样点数，ncxline为横测线(cxline)方向的采样点数，ninline为纵测线(inline)方向的采样点数，nhx为纵测线(inline)方向偏移距采样点数，nhy为横测线(cxline)方向偏移距采样点数。该五维数据体I(nz,ncxline,ninline,nhx,nhy)的数据规模量级，在当前主流的高精度三维地震采集条件会高达数TB，在高密度三维地震采集条件下会达到十数TB甚至数十TB规模。
目前，地震勘探中通常采用的Kirchhoff积分法偏移OVG道集输出技术，首先需要将预处理后CMP域道集数据分选为许多个独立的OVT(Offset Vector Tile)道集，然后对分选后OVT道集逐个进行偏移成像，并得到对应每个OVT偏移的成像体，最终再对所有OVT成像数据体进行重新分选，得到最后的共成像点OVG道集。上述方法需要将偏移前的基础地震数据和偏移后的成像数据进行多次道集分选和重排序，实现过程繁琐，效率极低。特别是针对当前高密度地震采集的海量地震数据，道集重排序不仅产生极大的数据读写和计算耗时，同时也带来极大的存储资源需求压力。因此，亟需形成一种能够以常规预处理后CMP道集为成像所需基础数据，并直接偏移得到OVG道集的成像方法。
在申请号：201510924779.0的中国专利申请中，涉及到一种基于双路径积分的速度模型的建立方法及系统，所述方法包括：采集目标数据；确定处理所述目标数据所需的多个待选速度模型；获取每个待选速度模型对应的共成像点道集；基于共成像点道集，获取每个待选速度模型对应的曲率参数及偏移结果；对每个所述偏移结果利用其对应的曲率参数进行两次加权叠加，得到两个路径积分成像剖面；分别计算每个待选速度模型对应的成像点道集的相干值和路径积分成像剖面的相干值；基于相干值，确定待选速度模型中的目标速度模型，目标速度模型用以建立最终偏移速度模型。
在申请号：201510924867.0的中国专利申请中，涉及到一种基于双路径积分的速度模型的建立方法及系统，所述方法包括：采集目标数据；确定处理所述目标数据所需的多个待选速度模型；获取每个待选速度模型对应的共成像点道集，进而获取每个待选速度模型对应道集的相干谱及偏移结果；对每个偏移结果利用其对应的相干谱进行两次加权叠加，得到每个两个路径积分成像剖面；分别计算每个待选速度模型对应的共成像点道集的相干值和路径积分成像剖面的相干值；基于相干值，确定待选速度模型中的目标速度模型，目标速度模型用以建立速度模型。
上述两专利所述内容为一种速度分析与建模的方法，虽涉及到偏移成像内容，但其重点为一种速度分析与建模方法，偏移成像及成像道集输出，仅为其实现过程中数据准备的一个基本操作步骤。上述专利中，作为判断成像道集同相轴相干性的准则，相干值和曲率参数是等价的，或者说二者是对同一参量的进行描述的两种形式。因此，两专利所要申请保护的内容是一致的。利用成像道集同相轴相干谱或者曲率参数进行速度分析，是速度分析与建模中的一种常用做法。通过多个待选模型及其对应的偏移结果，并利用一定的判断标准确定最终结果的办法，是蒙特卡洛类方法的具体表现形式。现有技术中没有一种能够以常规预处理后CMP道集为成像所需基础数据，并直接偏移得到OVG道集的成像方法。
为此我们发明了一种新的积分法成像实现OVG道集直接输出的方法，解决了以上技术问题。
发明内容
本发明的目的是提供一种更加高效地实现Kirchhoff积分法偏移OVG道集直接输出的方法。
本发明的目的可通过如下技术措施来实现：积分法成像实现OVG道集直接输出的方法，该积分法成像实现OVG道集直接输出的方法包括：步骤1，获取偏移计算的参数；步骤2，根据成像目标线位置及偏移孔径，获取所需的预处理后CMP道集；步骤3，根据参与计算的进程数目，将获取的CMP道集数据均分到各个进程；步骤4，读取偏移所需的速度场；步骤5，计算该目标线成像所需的旅行时场；步骤6，计算当前输入地震道对应的OVG成像输出参数信息；步骤7，输出该成像目标线的OVG道集；步骤8，循环计算下一成像目标线。
本发明的目的还可通过如下技术措施来实现：
在步骤1，获取的偏移计算的参数包括偏移孔径、成像目标线范围、成像时间或深度采样信息、计算进程个数。
在步骤2，在获取CMP道集数据后，根据偏移孔径及成像后数据线间距计算所获取输入CMP数据的线范围：
Figure BDA0002578499790000031
其中，Nline表示所需获取的CMP道集数据的线范围，Lapt为偏移计算孔径，Dline为成像后数据线间距。
在步骤2，判断当前成像目标线与前一目标线是否存在偏移孔径内的地震数据重叠；当存在重叠部分时，从数据存储介质获取孔径内新增的非重叠部分地震数据，重叠部分的地震数据则直接从内存中相应获取，以此节省地震数据读取时间。
在步骤3中，根据参与计算的进程数目Nproc，将Nline条线的地震数据均分到Nproc个进程，每个计算进程的CMP道集数据线数量Nline_per_node为：
Nline_per_node＝Nline/Nproc。
在步骤4，根据当前成像目标线需要的CMP数据线范围，读取偏移所需的速度场；判断当前目标线与前一目标线是否存在偏移孔径内的速度场重叠；当存在重叠部分时，从数据存储介质获取孔径内新增的非重叠部分速度场，重叠部分的速度场则直接从内存中相应获取，以此节省速度场读取时间。
在步骤5，根据当前成像目标线需要的CMP数据线范围，计算该目标线成像所需的旅行时场；判断当前目标线与前一目标线是否存在偏移孔径内的旅行时场重叠；当存在重叠部分时，仅计算孔径内非重叠部分的旅行时场，重叠部分的旅行时场则直接从内存中相应获取，以此节省旅行时场的计算时间。
在步骤6，根据所获取的地震道数据的炮点和检波点相对位置关系，利用如下公式计算当前输入地震道对应的OVG成像输出参数信息：
Figure BDA0002578499790000041
Figure BDA0002578499790000042
Figure BDA0002578499790000043
Figure BDA0002578499790000044
其中，hx、hy分别为偏移距向量在x、y方向的偏移距分量；dhx、dhy为x、y方向的偏移距增量；ihx、ihy分别为x、y方向的偏移距向量片分组号；
Figure BDA0002578499790000045
为地表方位角；offset为偏移距向量的模，即绝对偏移距；
根据上述OVG道集输出参数信息，在偏移过程中将成像结果自动投影到对应位置，依次实现该成像目标线的OVG道集输出。
在步骤7，针对当前成像目标线，完成孔径内全部相关输入数据的偏移计算，并对相同成像位置点成像结果求和，实现该成像目标线的OVG道集输出。
在步骤8，循环计算下一条成像目标线的OVG成像结果，重复步骤2至步骤7，直至完成所有成像目标线的OVG道集输出。
本发明中的积分法成像实现OVG道集直接输出的方法，在高密度地震数据Kirchhoff积分法偏移过程中以共中心点(CMP，Common Mid-point)域地震道集为成像所需基础数据，高效偏移并直接生成共成像点OVG(Offset Vector Gather)道集。首先，该方法以成像输出目标线作为主循环，逐条成像线依次运算。针对某条成像输出线，根据偏移孔径确定其成像计算所需数据的纵测线(inline)范围。然后，将所有需要的纵测线(inline)范围内的地震数据平均分配到各个计算进程的内存中。然后，根据当前成像输出线所需要的地震数据坐标范围，计算地震旅行时场，并将对应该条成像目标线的旅行时场分发给所有的计算进程。最后，实现每道数据在成像目标线上的偏移计算，并根据地震数据的纵、横向偏移距，将像投影到该条线的OVG空间，从而得到OVG道集。本发明基于线偏移成像策略，均衡考虑地震数据、速度场、旅行时场、OVG成像道集等对计算机内存的要求，达到利用当今主流计算机群实现高密度地震数据高效偏移和OVG成像道集直接输出的目的。
在综合考虑内存使用量、地震数据读写量、旅行时数据读写量、成像结果读写量和偏移成像精度等影响因素基础上，设计了一个面向大规模海量地震数据体的三维Kirchhoff积分法偏移及OVG道集输出的并行实现方案。该方案基于当前主流高性能计算集群特点，以CMP道集作为成像所需基础数据，利用输入道观点实现Kirchhoff积分法叠前偏移，并输出全方位角度域OVG道集。
本发明的积分法成像实现OVG道集直接输出的方法，由于每次只对一条目标线进行成像，即使在成像阶段将成像结果分选到OVG域，成像体的规模仍旧是可接受的，并不会存在过大的内存需求而造成无法计算。在所有进程完成当前目标线的成像后，规约求和得到成像道集，并写入到指定输出文件，即完成一条地震线的OVG道集成像。在完成一条线的数据成像后，内存中去除对下一条目标线成像无贡献的地震数据，并增量读入对下一条目标线成像有贡献的CMP道集数据，重复偏移计算和OVG道集输出步骤，实现新目标线的成像。
本发明的积分法成像实现OVG道集直接输出的方法，相较于输出道思想实现的Kirchhoff积分法偏移，所采用策略增加了一次获取的数据量，减少了输出成像体的数据量。在对成像目标线逐条偏移的情况下，通常规模的三维地震数据体，一条目标线成像所需的地震数据获取量通常在500GB以内，如果50个以上节点参与运算，每个节点需要读入的地震数据在10GB以内。对输出OVG道集，每条线需要的存储也在10GB以内。总的内存需求在20GB左右。并且，随着参与的计算节点的增加，分担地震数据的节点增多，单个进程的内存需求会进一步降低，从而进一步降低存储需求的压力，提高计算效率。
本发明针对整个三维工区OVG数据体总量庞大的现象，设计一次输出一条目标线的OVG道集数据。相较整个三维工区的OVG道集数据量，内存需求降低两个数量级以上，大大提高了本发明技术实现的可行性和实用性。同时，本发明在成像过程中直接得到共成像点OVG道集。不需要事先抽取OVT数据片并逐个偏移，然后再进行重新排序得到OVG道集，避免了两次极其耗时的大规模数据重排序工作量，大大降低了成像前和成像后地震数据的读写量和磁盘占有量，对于高密度地震数据具有非常好的适应性。
附图说明
图1为本发明的积分法成像实现OVG道集直接输出的方法的一具体实施例的流程图；
图2为本发明的一具体实施例中在某三维实际资料实施例中应用得到的快慢维为方位角-偏移距排序的OVG道集的示意图；
图3为本发明在某三维实际资料实施例中应用得到的快慢维为偏移距-方位角排序的OVG道集的示意图；
图4为本发明在某三维实际资料实施例中应用得到OVG叠加剖面的示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。
本发明的积分法成像实现OVG道集直接输出的方法包括了以下步骤：
步骤1，获取偏移计算的参数，包括偏移孔径、成像目标线范围、成像时间或深度采样信息、计算进程个数等；
步骤2，根据成像目标线位置及偏移孔径，获取所需的预处理后CMP道集。然后，根据偏移孔径及成像后数据线间距计算所获取输入CMP数据的线范围：
Figure BDA0002578499790000071
其中，Nline表示所需获取的CMP道集数据的线范围，Lapt为偏移计算孔径，Dline为成像后数据线间距。
在本步骤2，判断当前成像目标线与前一目标线是否存在偏移孔径内的地震数据重叠；当存在重叠部分时，从数据存储介质获取孔径内新增的非重叠部分地震数据，重叠部分的地震数据则直接从内存中相应获取，以此节省地震数据读取时间。
本步骤所需获取的基础地震数据为常规预处理后的CMP道集，不需要将地震数据重排序为OVT(Offset Vector Tile)道集，以降低成像前基础数据准备过程的时间开销。
步骤3，根据参与计算的进程数目Nproc，将Nline条线的地震数据均分到Nproc个进程。每个计算进程的CMP道集数据线数量为：
Nline_per_node＝Nline/Nproc
步骤4，根据当前成像目标线需要的CMP数据线范围，读取偏移所需的速度场；判断当前目标线与前一目标线是否存在偏移孔径内的速度场重叠；当存在重叠部分时，从数据存储介质获取孔径内新增的非重叠部分速度场，重叠部分的速度场则直接从内存中相应获取，以此节省速度场读取时间。
步骤5，根据当前成像目标线需要的CMP数据线范围，计算该目标线成像所需的旅行时场。其中，需判断当前目标线与前一目标线是否存在偏移孔径内的旅行时场重叠。当存在重叠部分时，仅计算孔径内非重叠部分的旅行时场，重叠部分的旅行时场则直接从内存中相应获取，以此节省旅行时场的计算时间。
步骤6，根据所获取的地震道数据的炮点和检波点相对位置关系，利用如下公式计算当前输入地震道对应的OVG成像输出参数信息，
Figure BDA0002578499790000081
Figure BDA0002578499790000085
Figure BDA0002578499790000082
Figure BDA0002578499790000083
其中，hx、hy分别为偏移距向量在x、y方向的偏移距分量；dhx、dhy为x、y方向的偏移距增量；ihx、ihy分别为x、y方向的偏移距向量片分组号；
Figure BDA0002578499790000084
为地表方位角；offset为偏移距向量的模，即绝对偏移距。
根据上述OVG道集输出参数信息，在偏移过程中将成像结果自动投影到对应位置，依次实现该成像目标线的OVG道集输出；
不同于常规以OVT道集为输入的成像方法是将地震数据成像到相应的OVT域，并在所有OVT数据成像结束后将成像结果重新分选到OVG域。本发明根据输入地震道数据的纵、横向偏移距，将其成像结果直接投影到对应的OVG域，避免成像后道集数据的重排序。
步骤7，针对当前成像目标线，完成孔径内全部相关输入数据的偏移计算，并对相同成像位置点成像结果求和，实现该成像目标线的OVG道集输出；
本发明以成像目标线为主循环，一次只需要计算并存储当前成像目标线的旅行时场和OVG成像道集，降低旅行时场和OVG成像道集的数据规模。在主流并行计算机框架下，所需地震数据均分后被分发给各个计算进程分散承担，有效降低本发明对计算内存的需求。
步骤8，循环计算下一条成像目标线的OVG成像结果，重复步骤2至步骤7，直至完成所有成像目标线的OVG道集输出。
在应用本发明的一具体实施例中，如图1所示，图1为本发明的积分法成像实现OVG道集直接输出的方法的流程图。
步骤101，获取偏移计算的参数。包括偏移孔径、成像目标线范围、成像时间或深度采样信息、计算进程个数等。
步骤102，根据成像目标线位置及偏移孔径，获取所需的预处理后CMP道集。然后，根据偏移孔径及成像后数据线间距计算所获取输入CMP数据的线范围：
Figure BDA0002578499790000091
其中，Nline表示所需获取的CMP道集数据的线范围，Lapt为偏移计算孔径，Dline为成像后数据线间距。
在本步骤，需判断当前成像目标线与前一目标线是否存在偏移孔径内的地震数据重叠；当存在重叠部分时，从数据存储介质获取孔径内新增的非重叠部分地震数据，重叠部分的地震数据则直接从内存中相应获取，以此节省地震数据读取时间。
步骤103，根据参与计算的进程数目Nproc，将Nline条线的地震数据均分到Nproc个进程。每个计算进程的CMP道集数据线数量为：
Nline_per_node＝Nline/Nproc
步骤104，针对当前需要计算的成像目标线，根据步骤102中确定的成像所需CMP线范围，读取偏移所需的速度场，判断当前目标线与前一目标线是否存在偏移孔径内的速度场重叠。当存在重叠部分时，从数据存储介质获取孔径内新增的非重叠部分速度场，重叠部分的速度场则直接从内存中相应获取。将最终速度场分发给各个计算进程。
步骤105，针对当前需要计算的成像目标线，根据步骤102中确定的成像所需CMP线范围，判断当前目标线与前一目标线是否存在偏移孔径内的旅行时重叠。当存在重叠部分时，仅计算孔径内非重叠部分的旅行时场，重叠部分的旅行时场则直接从内存中相应获取，以此节省旅行时场的计算时间。然后将该最终旅行时场分发到各个计算进程。
步骤106，每个独立的计算进程，基于其自身对应获取的地震数据、速度场和旅行时场进行偏移运算，并将偏移结果相应位置叠加后得到成像结果。同时，针对每道地震数据，根据其偏移前的纵向偏移距hy和横向偏移距hx参量，依据下列公式计算的成像参数信息，将成像结果投影到对应的OVG域。
Figure BDA0002578499790000101
Figure BDA0002578499790000102
Figure BDA0002578499790000103
Figure BDA0002578499790000104
其中，hx、hy分别为偏移距向量在x、y方向的偏移距分量；dhx、dhy为x、y方向的偏移距增量；ihx、ihy分别为x、y方向的偏移距向量片分组号；
Figure BDA0002578499790000105
为地表方位角；offset为偏移距向量的模，即绝对偏移距。
每个共成像点OVG域道集存在两种道排序方式：一种道集排序方式如图2所示，其快维为绝对偏移距offset，慢维为方位角azimuth；另一种道集排序方式如图3所示，其快维为方位角azimuth，慢维为绝对偏移距offset。其中，绝对偏移距
Figure BDA0002578499790000106
方位角azimuth为当前道炮检点连线与正北方向的顺时针旋转角度。
不同于常规基于OVT道集的成像方法将地震数据成像到相应的OVT域，并在所有OVT数据成像结束后将其重新分选到OVG域。本发明根据CMP域地震道数据的纵、横向偏移距，将其成像结果直接投影到对应的OVG域，省却了成像后OVG道集数据重排序的时间。
步骤107，针对处在当前循环中的成像目标线，各个进程分别完成所负责CMP道集偏移计算后，对各进程负责的各个独立成像结果进行规约，实现该成像目标线上各个共成像点的OVG道集输出。图4为本发明在某三维实际资料实施例中应用得到某条目标线的OVG道集叠加剖面。
由于本发明所述方法以成像目标线为主循环，每次只需要计算并存储当前成像目标线的地震数据、速度场、旅行时场和OVG道集，大大降低了偏移计算过程中的数据规模。同时，由于本发明采用多节点多进程并行计算策略，在偏移过程中，设计将所需地震数据、速度场、旅行时场合理分发给各个进程分散承担，有效降低了对计算内存的需求。相较于常规Kirchhoff积分法偏移，本发明所采用策略虽一定程度上增加了单次获取的地震数据量，但很大程度上减少了输出成像体的数据量。在对成像目标线逐条偏移的情况下，通常规模的三维地震数据体，一条目标线成像所需的地震数据获取量通常在500GB以内，如果50个以上节点参与运算，每个节点需要读入的地震数据在10GB以内。对于要输出的OVG道集，每条线需要的存储也在10GB以内。总的内存需求在20GB左右。并且，随着参与的计算节点的增加，分担所获取数据的节点数目增多，单个进程的内存需求会进一步降低，从而进一步降低了对存储需求的压力，提高计算效率。
本发明针对整个三维工区的OVG道集数据体总量庞大的现象，设计一次输出一条目标线的OVG道集数据。相较整个三维工区的OVG道集数据量，内存需求降低两个数量级以上，大大提高了本发明技术实现的可行性和实用性。同时，本发明在成像过程中直接得到共成像点OVG道集。不需要事先抽取OVT数据片并逐个偏移，然后再进行重新排序得到OVG道集，避免了两次极其耗时的大规模数据重排序工作量，大大降低了成像前和成像后地震数据的读写量和磁盘占有量，对于高密度地震数据具有非常好的适应性。
步骤108，循环计算下一条目标线的成像结果，重复步骤2至步骤6，直至完成所有成像目标线的OVG道集输出。
本发明的种积分法成像实现OVG道集直接输出的方法，在整体设计过程中，综合考虑了内存使用量、地震数据读写量、旅行时数据读写量、成像结果读写量和偏移成像精度等因素影响，优化形成了一个面向大规模海量地震数据，特别是高密度地震采集数据的三维Kirchhoff积分法偏移，以及OVG道集直接输出的并行实现方案。该方案基于当前主流高性能计算集群特点，以CMP道集作为成像所需基础数据，利用输入道观点实现Kirchhoff积分法叠前偏移，并输出全方位角度域共成像点OVG道集。

Claims (10)

1.积分法成像实现OVG道集直接输出的方法，其特征在于，该积分法成像实现OVG道集直接输出的方法包括：
步骤1，获取偏移计算的参数；
步骤2，根据成像目标线位置及偏移孔径，获取所需的预处理后CMP道集；
步骤3，根据参与计算的进程数目，将获取的CMP道集数据按线数量均分到各个进程；
步骤4，读取偏移所需的速度场；
步骤5，计算该目标线成像所需的旅行时场；
步骤6，计算当前输入地震道对应的OVG成像输出参数信息；
步骤7，输出该成像目标线的OVG道集；
步骤8，循环计算下一成像目标线。
2.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法，其特征在于，在步骤1，获取的偏移计算的参数包括偏移孔径、成像目标线范围、成像时间或深度采样信息、计算进程个数。
3.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法，其特征在于，在步骤2，在获取CMP道集数据后，根据偏移孔径及成像后数据线间距计算所获取输入CMP数据的线范围：
Figure FDA0002578499780000011
其中，Nline表示所需获取的CMP道集数据的线范围，Lapt为偏移计算孔径，Dline为成像后数据线间距。
4.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法，其特征在于，在步骤2，判断当前成像目标线与前一目标线是否存在偏移孔径内的地震数据重叠；当存在重叠部分时，从数据存储介质获取孔径内新增的非重叠部分地震数据，重叠部分的地震数据则直接从内存中相应获取，以此节省地震数据读取时间。
5.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法，其特征在于，在步骤3中，根据参与计算的进程数目Nproc，将Nline条线的地震数据均分到Nproc个进程，每个计算进程的CMP道集数据线数量Nline_per_node为：
Nline_per_node＝Nline/Nproc。
6.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法，其特征在于，在步骤4，根据当前成像目标线需要的CMP数据线范围，读取偏移所需的速度场；判断当前目标线与前一目标线是否存在偏移孔径内的速度场重叠；当存在重叠部分时，从数据存储介质获取孔径内新增的非重叠部分速度场，重叠部分的速度场则直接从内存中相应获取，以此节省速度场读取时间。
7.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法，其特征在于，在步骤5，根据当前成像目标线需要的CMP数据线范围，计算该目标线成像所需的旅行时场；判断当前目标线与前一目标线是否存在偏移孔径内的旅行时场重叠；当存在重叠部分时，仅计算孔径内非重叠部分的旅行时场，重叠部分的旅行时场则直接从内存中相应获取，以此节省旅行时场的计算时间。
8.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法，其特征在于，在步骤6，根据所获取的地震道数据的炮点和检波点相对位置关系，利用如下公式计算当前输入地震道对应的OVG成像输出参数信息：
Figure FDA0002578499780000021
Figure FDA0002578499780000022
Figure FDA0002578499780000023
Figure FDA0002578499780000024
其中，hx、hy分别为偏移距向量在x、y方向的偏移距分量；dhx、dhy为x、y方向的偏移距增量；ihx、ihy分别为x、y方向的偏移距向量片分组号；
Figure FDA0002578499780000025
为地表方位角；offset为偏移距向量的模，即绝对偏移距；
根据上述OVG道集输出参数信息，在偏移过程中将成像结果自动投影到对应位置，依次实现该成像目标线的OVG道集输出。
9.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法，其特征在于，在步骤7，针对当前成像目标线，完成孔径内全部相关输入数据的偏移计算，并对相同成像位置点成像结果求和，实现该成像目标线的OVG道集输出。
10.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法，其特征在于，在步骤8，循环计算下一条成像目标线的OVG成像结果，重复步骤2至步骤7，直至完成所有成像目标线的OVG道集输出。
CN202010661120.1A 2020-07-10 2020-07-10 积分法成像实现ovg道集直接输出的方法 Active CN113917536B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010661120.1A CN113917536B (zh) 2020-07-10 2020-07-10 积分法成像实现ovg道集直接输出的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010661120.1A CN113917536B (zh) 2020-07-10 2020-07-10 积分法成像实现ovg道集直接输出的方法

Publications (2)

Publication Number Publication Date
CN113917536A true CN113917536A (zh) 2022-01-11
CN113917536B CN113917536B (zh) 2024-06-25

Family

ID=79231970

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010661120.1A Active CN113917536B (zh) 2020-07-10 2020-07-10 积分法成像实现ovg道集直接输出的方法

Country Status (1)

Country Link
CN (1) CN113917536B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050088913A1 (en) * 2003-04-09 2005-04-28 Norsar Method for simulating local prestack depth migrated seismic images
US20050143921A1 (en) * 2003-12-12 2005-06-30 Exxonmobil Upstream Research Company Method for seismic imaging in geologically complex formations
US20090310443A1 (en) * 2008-06-11 2009-12-17 Baker Hughes Incorporated Vector Migration of Virtual Source VSP Data
CN101957455A (zh) * 2010-09-20 2011-01-26 中国海洋石油总公司 三维保幅叠前时间偏移方法
CN102692644A (zh) * 2012-05-30 2012-09-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 生成深度域成像道集的方法
CN104570081A (zh) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 一种积分法叠前时间偏移地震资料处理方法及系统
CN104849751A (zh) * 2015-05-15 2015-08-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 叠前地震资料成像的方法
WO2016041185A1 (zh) * 2014-09-19 2016-03-24 杨顺伟 一种高效叠前时间偏移速度分析方法
US20160202369A1 (en) * 2015-01-13 2016-07-14 Cgg Services Sa Using an offset vector tile gather to image a subsurface
US20160334528A1 (en) * 2013-12-31 2016-11-17 Cgg Services Sa Systems and methods for characterizing subterranean formations utilizing azimuthal data
CN106526666A (zh) * 2016-09-29 2017-03-22 中国石油天然气集团公司 叠前深度偏移方法、装置及系统
CN110187383A (zh) * 2019-05-27 2019-08-30 中海石油(中国)有限公司 一种海上宽方位地震数据cov道集快速分选方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050088913A1 (en) * 2003-04-09 2005-04-28 Norsar Method for simulating local prestack depth migrated seismic images
US20050143921A1 (en) * 2003-12-12 2005-06-30 Exxonmobil Upstream Research Company Method for seismic imaging in geologically complex formations
US20090310443A1 (en) * 2008-06-11 2009-12-17 Baker Hughes Incorporated Vector Migration of Virtual Source VSP Data
CN101957455A (zh) * 2010-09-20 2011-01-26 中国海洋石油总公司 三维保幅叠前时间偏移方法
CN102692644A (zh) * 2012-05-30 2012-09-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 生成深度域成像道集的方法
CN104570081A (zh) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 一种积分法叠前时间偏移地震资料处理方法及系统
US20160334528A1 (en) * 2013-12-31 2016-11-17 Cgg Services Sa Systems and methods for characterizing subterranean formations utilizing azimuthal data
WO2016041185A1 (zh) * 2014-09-19 2016-03-24 杨顺伟 一种高效叠前时间偏移速度分析方法
US20160202369A1 (en) * 2015-01-13 2016-07-14 Cgg Services Sa Using an offset vector tile gather to image a subsurface
CN104849751A (zh) * 2015-05-15 2015-08-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 叠前地震资料成像的方法
CN106526666A (zh) * 2016-09-29 2017-03-22 中国石油天然气集团公司 叠前深度偏移方法、装置及系统
CN110187383A (zh) * 2019-05-27 2019-08-30 中海石油(中国)有限公司 一种海上宽方位地震数据cov道集快速分选方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
MIKE PERZ等: "Considerations in Azimuthal Processing and Velocity inversion for Unconventional plays", UNCONVENTIONAL RESOURCES TECHNOLOGY CONFERENCE, HOUSTON, TEXAS, 25 July 2018 (2018-07-25), pages 4182 *
王华忠;刘少勇;孔祥宁;蔡杰雄;方伍宝;: "大规模三维地震数据Kirchhoff叠前深度偏移及其并行实现", 石油地球物理勘探, no. 03, pages 404 - 413 *
王延光: "地震叠前深度偏移技术进展及应用问题与对策", 油气地质与采收率, vol. 24, no. 04, pages 1 - 8 *
王狮虎;罗国安;寇芹;田振平;张丽梅;聂红梅;: "GeoEast系统积分法叠前时间偏移技术", 石油工业计算机应用, no. 03, 20 September 2016 (2016-09-20), pages 20 - 25 *
袁刚: "宽方位数据的炮检距向量片域处理及偏移道集校平方法", 石油物探, vol. 55, no. 01, 31 January 2016 (2016-01-31), pages 84 - 90 *
邓笋根;: "克希霍夫偏移成像分块策略对文件I/O的影响", 科研信息化技术与应用, no. 05, pages 25 - 29 *

Also Published As

Publication number Publication date
CN113917536B (zh) 2024-06-25

Similar Documents

Publication Publication Date Title
AU2013296416B2 (en) Method for editing a multi-point facies simulation
CN109100795B (zh) 一种面元的炮检点布设方法、装置及系统
CN109471171B (zh) 一种混叠地震数据分离的方法、装置及系统
AU2020328458B2 (en) Correlation matrix for simultaneously correlating multiple wells
EP2880470A2 (en) Hybrid method of combining multipoint statistic and object-based methods for creating reservoir property models
CN107204040B (zh) 多点地质统计学建模方法及装置、计算机存储介质
Ghosh et al. A strategy for adaptive sampling of multi-fidelity gaussian processes to reduce predictive uncertainty
Ding et al. Performance evaluation of GPU-accelerated spatial interpolation using radial basis functions for building explicit surfaces
Hou et al. Full tensor gravity gradiometry data inversion: Performance analysis of parallel computing algorithms
CN103235879A (zh) 基于小生境主从式并行遗传算法的双相介质参数反演方法
CN113917536A (zh) 积分法成像实现ovg道集直接输出的方法
CN113360187B (zh) 一种基于CUDA与OpenMP的三维Kriging算法协同加速方法
CN111159869B (zh) 一种砂岩铀矿床与古河道间定量参数的自动计算方法
CN115128671A (zh) 一种基于深度学习技术的三维地震相自动识别方法
CN112241020A (zh) 稀疏地震数据采集中确定欠采样率的方法和装置
CN113419278A (zh) 一种基于状态空间模型与支持向量回归的井震联合多目标同时反演方法
CN107730586B (zh) 一种地层建模的方法和系统
CN108537883B (zh) 一种基于MapReduce框架的地质建模方法
CN108646288B (zh) 近地表模型的建立方法和装置
Dossary et al. Self-Organizing Maps for Regions Exploring and Identification Based on Geological Signatures Similarities and Anomalies
CN113568047B (zh) 一种叠前成像道集生成方法和装置
CN112965112B (zh) 地下储层微构造三维散点的生成方法
CN116976218A (zh) 多磁偶极子反演方法、装置和电子设备
Snyder et al. Algebraic Topological Methods for the Analysis and Modeling of Protein Data
CN114428318A (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