CN113917536B - 积分法成像实现ovg道集直接输出的方法 - Google Patents
积分法成像实现ovg道集直接输出的方法 Download PDFInfo
- Publication number
- CN113917536B CN113917536B CN202010661120.1A CN202010661120A CN113917536B CN 113917536 B CN113917536 B CN 113917536B CN 202010661120 A CN202010661120 A CN 202010661120A CN 113917536 B CN113917536 B CN 113917536B
- Authority
- CN
- China
- Prior art keywords
- imaging
- offset
- ovg
- data
- gather
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 195
- 238000000034 method Methods 0.000 title claims abstract description 96
- 238000004364 calculation method Methods 0.000 claims abstract description 39
- 238000013500 data storage Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 abstract description 12
- 238000013508 migration Methods 0.000 description 12
- 230000005012 migration Effects 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
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道集直接输出的方法。
背景技术
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数据的线范围:
其中,Nline表示所需获取的CMP道集数据的线范围,Lapt为偏移计算孔径,Dline为成像后数据线间距。
在步骤2,判断当前成像目标线与前一目标线是否存在偏移孔径内的地震数据重叠;当存在重叠部分时,从数据存储介质获取孔径内新增的非重叠部分地震数据,重叠部分的地震数据则直接从内存中相应获取,以此节省地震数据读取时间。
在步骤3中,根据参与计算的进程数目Nproc,将Nline条线的地震数据均分到Nproc个进程,每个计算进程的CMP道集数据线数量Nline_per_node为:
Nline_per_node=Nline/Nproc。
在步骤4,根据当前成像目标线需要的CMP数据线范围,读取偏移所需的速度场;判断当前目标线与前一目标线是否存在偏移孔径内的速度场重叠;当存在重叠部分时,从数据存储介质获取孔径内新增的非重叠部分速度场,重叠部分的速度场则直接从内存中相应获取,以此节省速度场读取时间。
在步骤5,根据当前成像目标线需要的CMP数据线范围,计算该目标线成像所需的旅行时场;判断当前目标线与前一目标线是否存在偏移孔径内的旅行时场重叠;当存在重叠部分时,仅计算孔径内非重叠部分的旅行时场,重叠部分的旅行时场则直接从内存中相应获取,以此节省旅行时场的计算时间。
在步骤6,根据所获取的地震道数据的炮点和检波点相对位置关系,利用如下公式计算当前输入地震道对应的OVG成像输出参数信息:
其中,hx、hy分别为偏移距向量在x、y方向的偏移距分量;dhx、dhy为x、y方向的偏移距增量;ihx、ihy分别为x、y方向的偏移距向量片分组号;为地表方位角;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数据的线范围:
其中,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成像输出参数信息,
其中,hx、hy分别为偏移距向量在x、y方向的偏移距分量;dhx、dhy为x、y方向的偏移距增量;ihx、ihy分别为x、y方向的偏移距向量片分组号;为地表方位角;offset为偏移距向量的模,即绝对偏移距。
根据上述OVG道集输出参数信息,在偏移过程中将成像结果自动投影到对应位置,依次实现该成像目标线的OVG道集输出;
不同于常规以OVT道集为输入的成像方法是将地震数据成像到相应的OVT域,并在所有OVT数据成像结束后将成像结果重新分选到OVG域。本发明根据输入地震道数据的纵、横向偏移距,将其成像结果直接投影到对应的OVG域,避免成像后道集数据的重排序。
步骤7,针对当前成像目标线,完成孔径内全部相关输入数据的偏移计算,并对相同成像位置点成像结果求和,实现该成像目标线的OVG道集输出;
本发明以成像目标线为主循环,一次只需要计算并存储当前成像目标线的旅行时场和OVG成像道集,降低旅行时场和OVG成像道集的数据规模。在主流并行计算机框架下,所需地震数据均分后被分发给各个计算进程分散承担,有效降低本发明对计算内存的需求。
步骤8,循环计算下一条成像目标线的OVG成像结果,重复步骤2至步骤7,直至完成所有成像目标线的OVG道集输出。
在应用本发明的一具体实施例中,如图1所示,图1为本发明的积分法成像实现OVG道集直接输出的方法的流程图。
步骤101,获取偏移计算的参数。包括偏移孔径、成像目标线范围、成像时间或深度采样信息、计算进程个数等。
步骤102,根据成像目标线位置及偏移孔径,获取所需的预处理后CMP道集。然后,根据偏移孔径及成像后数据线间距计算所获取输入CMP数据的线范围:
其中,Nline表示所需获取的CMP道集数据的线范围,Lapt为偏移计算孔径,Dline为成像后数据线间距。
在本步骤,需判断当前成像目标线与前一目标线是否存在偏移孔径内的地震数据重叠;当存在重叠部分时,从数据存储介质获取孔径内新增的非重叠部分地震数据,重叠部分的地震数据则直接从内存中相应获取,以此节省地震数据读取时间。
步骤103,根据参与计算的进程数目Nproc,将Nline条线的地震数据均分到Nproc个进程。每个计算进程的CMP道集数据线数量为:
Nline_per_node=Nline/Nproc
步骤104,针对当前需要计算的成像目标线,根据步骤102中确定的成像所需CMP线范围,读取偏移所需的速度场,判断当前目标线与前一目标线是否存在偏移孔径内的速度场重叠。当存在重叠部分时,从数据存储介质获取孔径内新增的非重叠部分速度场,重叠部分的速度场则直接从内存中相应获取。将最终速度场分发给各个计算进程。
步骤105,针对当前需要计算的成像目标线,根据步骤102中确定的成像所需CMP线范围,判断当前目标线与前一目标线是否存在偏移孔径内的旅行时重叠。当存在重叠部分时,仅计算孔径内非重叠部分的旅行时场,重叠部分的旅行时场则直接从内存中相应获取,以此节省旅行时场的计算时间。然后将该最终旅行时场分发到各个计算进程。
步骤106,每个独立的计算进程,基于其自身对应获取的地震数据、速度场和旅行时场进行偏移运算,并将偏移结果相应位置叠加后得到成像结果。同时,针对每道地震数据,根据其偏移前的纵向偏移距hy和横向偏移距hx参量,依据下列公式计算的成像参数信息,将成像结果投影到对应的OVG域。
其中,hx、hy分别为偏移距向量在x、y方向的偏移距分量;dhx、dhy为x、y方向的偏移距增量;ihx、ihy分别为x、y方向的偏移距向量片分组号;为地表方位角;offset为偏移距向量的模,即绝对偏移距。
每个共成像点OVG域道集存在两种道排序方式:一种道集排序方式如图2所示,其快维为绝对偏移距offset,慢维为方位角azimuth;另一种道集排序方式如图3所示,其快维为方位角azimuth,慢维为绝对偏移距offset。其中,绝对偏移距方位角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 (4)
1.积分法成像实现OVG道集直接输出的方法,其特征在于,该积分法成像实现OVG道集直接输出的方法包括:
步骤1,获取偏移计算的参数;
步骤2,根据成像目标线位置及偏移孔径,获取所需的预处理后CMP道集;
步骤3,根据参与计算的进程数目,将获取的CMP道集数据按线数量均分到各个进程;
步骤4,读取偏移所需的速度场;
步骤5,计算该目标线成像所需的旅行时场;
步骤6,计算当前输入地震道对应的OVG成像输出参数信息;
步骤7,输出该成像目标线的OVG道集;
步骤8,循环计算下一成像目标线;
在步骤1,获取的偏移计算的参数包括偏移孔径、成像目标线范围、成像时间或深度采样信息、计算进程个数;
在步骤2,在获取CMP道集数据后,根据偏移孔径及成像后数据线间距计算所获取输入CMP数据的线范围:
其中,Nline表示所需获取的CMP道集数据的线范围,Lapt为偏移计算孔径,Dline为成像后数据线间距;
判断当前成像目标线与前一目标线是否存在偏移孔径内的地震数据重叠;当存在重叠部分时,从数据存储介质获取孔径内新增的非重叠部分地震数据,重叠部分的地震数据则直接从内存中相应获取,以此节省地震数据读取时间;
在步骤3中,根据参与计算的进程数目Nproc,将Nline条线的地震数据均分到Nproc个进程,每个计算进程的CMP道集数据线数量Nline_per_node为:
Nline_per_node=Nline/Nproc;
在步骤6,根据所获取的地震道数据的炮点和检波点相对位置关系,利用如下公式计算当前输入地震道对应的OVG成像输出参数信息:
其中,hx、hy分别为偏移距向量在x、y方向的偏移距分量;dhx、dhy为x、y方向的偏移距增量;ihx、ihy分别为x、y方向的偏移距向量片分组号;为地表方位角;offset为偏移距向量的模,即绝对偏移距;
根据上述OVG道集输出参数信息,在偏移过程中将成像结果自动投影到对应位置,依次实现该成像目标线的OVG道集输出;
在步骤7,针对当前成像目标线,完成孔径内全部相关输入数据的偏移计算,并对相同成像位置点成像结果求和,实现该成像目标线的OVG道集输出。
2.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法,其特征在于,在步骤4,根据当前成像目标线需要的CMP数据线范围,读取偏移所需的速度场;判断当前目标线与前一目标线是否存在偏移孔径内的速度场重叠;当存在重叠部分时,从数据存储介质获取孔径内新增的非重叠部分速度场,重叠部分的速度场则直接从内存中相应获取,以此节省速度场读取时间。
3.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法,其特征在于,在步骤5,根据当前成像目标线需要的CMP数据线范围,计算该目标线成像所需的旅行时场;判断当前目标线与前一目标线是否存在偏移孔径内的旅行时场重叠;当存在重叠部分时,仅计算孔径内非重叠部分的旅行时场,重叠部分的旅行时场则直接从内存中相应获取,以此节省旅行时场的计算时间。
4.根据权利要求1所述的积分法成像实现OVG道集直接输出的方法,其特征在于,在步骤8,循环计算下一条成像目标线的OVG成像结果,重复步骤2至步骤7,直至完成所有成像目标线的OVG道集输出。
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 CN113917536A (zh) | 2022-01-11 |
CN113917536B true 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 (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104849751A (zh) * | 2015-05-15 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 叠前地震资料成像的方法 |
CN106526666A (zh) * | 2016-09-29 | 2017-03-22 | 中国石油天然气集团公司 | 叠前深度偏移方法、装置及*** |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO322089B1 (no) * | 2003-04-09 | 2006-08-14 | Norsar V Daglig Leder | Fremgangsmate for simulering av lokale prestakk dypmigrerte seismiske bilder |
EP1721190A4 (en) * | 2003-12-12 | 2012-06-06 | Exxonmobil Upstream Res Co | SEISMIC IMAGERY FOR COMPLEX GEOLOGICAL FORMATIONS |
US8395967B2 (en) * | 2008-06-11 | 2013-03-12 | Baker Hughes Incorporated | Vector migration of virtual source VSP data |
CN101957455B (zh) * | 2010-09-20 | 2012-07-18 | 中国海洋石油总公司 | 三维保幅叠前时间偏移方法 |
CN102692644B (zh) * | 2012-05-30 | 2014-07-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 生成深度域成像道集的方法 |
CN104570081B (zh) * | 2013-10-29 | 2017-12-26 | 中国石油化工股份有限公司 | 一种积分法叠前时间偏移地震资料处理方法及*** |
WO2015101832A2 (en) * | 2013-12-31 | 2015-07-09 | Cgg Services Sa | Systems and methods for characterizing subterranean formations utilizing azimuthal data |
CN104769459A (zh) * | 2014-09-19 | 2015-07-08 | 杨顺伟 | 一种高效叠前时间偏移速度分析方法 |
US10191165B2 (en) * | 2015-01-13 | 2019-01-29 | Cgg Services Sas | Using an offset vector tile gather to image a subsurface |
CN110187383B (zh) * | 2019-05-27 | 2020-09-29 | 中海石油(中国)有限公司 | 一种海上宽方位地震数据cov道集快速分选方法 |
-
2020
- 2020-07-10 CN CN202010661120.1A patent/CN113917536B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104849751A (zh) * | 2015-05-15 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 叠前地震资料成像的方法 |
CN106526666A (zh) * | 2016-09-29 | 2017-03-22 | 中国石油天然气集团公司 | 叠前深度偏移方法、装置及*** |
Also Published As
Publication number | Publication date |
---|---|
CN113917536A (zh) | 2022-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110031896B (zh) | 基于多点地质统计学先验信息的地震随机反演方法及装置 | |
CN103824329B (zh) | 一种地质勘探三维可视化储量估算方法 | |
DK2356493T3 (en) | PROCEDURE FOR GEOLOGICAL MODELING OF SEISMIC DATA BY TRACK CORRELATION | |
CN109100795B (zh) | 一种面元的炮检点布设方法、装置及*** | |
CN108647261B (zh) | 基于气象数据离散点网格化处理的全球等值线绘制方法 | |
CN109459787B (zh) | 基于地震槽波全波形反演的煤矿井下构造成像方法及*** | |
CN109001813A (zh) | 一种压制多次波的方法、装置及*** | |
CN102692644A (zh) | 生成深度域成像道集的方法 | |
CN105403913A (zh) | 叠前深度偏移方法和装置 | |
WO2021127382A1 (en) | Full waveform inversion in the midpoint-offset domain | |
CN113917536B (zh) | 积分法成像实现ovg道集直接输出的方法 | |
CN113534259A (zh) | 一种可控震源高效采集实时叠前时间偏移成像方法 | |
CN103235879A (zh) | 基于小生境主从式并行遗传算法的双相介质参数反演方法 | |
CN104422953B (zh) | 一种提高地震叠前时间偏移计算效率的方法 | |
CN105738949B (zh) | 一种用于时移地震的九面元一致性并行处理方法 | |
CN112462427A (zh) | 多分量地震资料保幅角度域共成像点道集提取方法及*** | |
CN109581494B (zh) | 叠前偏移方法及装置 | |
CN109387872B (zh) | 表面多次波预测方法 | |
CN108072895B (zh) | 一种基于gpu的各向异性叠前逆时偏移优化方法 | |
CN105572730B (zh) | 三维复杂结构声波正演方法 | |
CN111159869B (zh) | 一种砂岩铀矿床与古河道间定量参数的自动计算方法 | |
CN105259577B (zh) | 一种确定地层界面的角度信息的方法及装置 | |
CN112241020B (zh) | 稀疏地震数据采集中确定欠采样率的方法和装置 | |
CN113568041A (zh) | 时移地震三维拖缆采集数据的可重复性分析方法及*** | |
CN115128671A (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 |