CN103064110B - 一种波动方程叠前偏移中的分层延拓成像方法 - Google Patents
一种波动方程叠前偏移中的分层延拓成像方法 Download PDFInfo
- Publication number
- CN103064110B CN103064110B CN201110317615.3A CN201110317615A CN103064110B CN 103064110 B CN103064110 B CN 103064110B CN 201110317615 A CN201110317615 A CN 201110317615A CN 103064110 B CN103064110 B CN 103064110B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- mfrac
- mover
- imaging
- 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 109
- 230000005012 migration Effects 0.000 title claims abstract description 34
- 238000013508 migration Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 95
- 239000010410 layer Substances 0.000 claims description 119
- 230000006870 function Effects 0.000 claims description 20
- 239000002356 single layer Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 abstract description 7
- 238000005070 sampling Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明为一种波动方程叠前偏移中的分层延拓成像方法,先根据节点内存大小,对成像空间进行分层;后从首个单平面波数据开始,依次对每一个单平面波数据进行延拓和成像操作;直至单平面波全部各层的单频数据延拓成像完成;所有单平面波延拓成像结束输出成像结果。本发明针对常规波动方程叠前偏移方法内存需求过大及并行节点利用率不高的问题,根据计算节点可用内存空间大小自动调整偏移计算的内存需求,做到内存开支与CPU运算均衡,充分利用计算机集群的并行计算能力,实现大规模地震数据的叠前偏移成像。
Description
技术领域
本发明属于三维波动方程叠前偏移成像处理方法,针对常规方法内存需求过大及并行节点利用率不高的问题,涉及地震叠前偏移成像和高性能并行计算等,可应用于油气勘探地震资料成像处理。
背景技术
三维有限差分波动方程叠前偏移存在计算过程中内存开支过大的问题,特别是当今三维地震勘探采集数据量越来越大,内存问题成为困扰波动方程叠前偏移方法实用化的一种重要因素。假设有一个三维工区inline方向点数为Nx,crossline方向点数Ny,成像深度点数为Nτ,偏移频率空间域复数类型数据体的单频个数为Nω,那么波动方程叠前偏移主进程内存需求为4×(2Nτ+2Nω)×Nx×Ny字节,从进程内存需求为4×Nτ×Nx×Ny字节。一个满覆盖面积100km2(目前数据通常可达到1000km2以上)的地震数据波动方程偏移时,主进程所需内存空间约为9G,一个从进程所需内存空间约为4G;一个8核16G内存(主流配置)的计算节点,至多可用4核进行偏移计算,如此造成并行节点无法得到有效利用。数据量越大,上述问题越严重。当数据量为上述情况的两三倍时,目前主流配置的计算机集群已无法进行波动方程叠前偏移运算。
发明内容
本发明为了解决现有技术中存在的非分层延拓成像方法占用内容过大,无法进行操作和偏移运算的技术问题,研发了一种波动方程叠前偏移中的分层延拓成像方法。本发明针对常规波动方程叠前偏移方法内存需求过大及并行节点利用率不高的问题,根据计算节点可用内存空间大小自动调整偏移计算的内存需求,做到内存开支与CPU运算均衡,充分利用计算机集群的并行计算能力,实现大规模地震数据的叠前偏移成像。
本发明为了解决上述现有技术中存在的技术问题,采用的技术方案如下,
一种波动方程叠前偏移中的分层延拓成像方法,所述延拓成像方法先根据节点内存大小,对成像空间进行分层;后从首个单平面波数据开始,依次对每一个单平面波数据进行延拓和成像操作;逐层延拓成像操作中包括主进程和从进程,所述主进程负责读写地震数据和向从进程分发数据,从进程负责该单层单频数据的延拓成像;各单层计算结束后,从进程把延拓至本层深度后的单频数据送回到主进程;由主进程再把数据发送出去;并进行下一层的延拓成像操作,直至单平面波全部各层的单频数据延拓成像完成;所有单平面波延拓成像结束输出成像结果。
具体的,所述的一种波动方程叠前偏移中的分层延拓成像方法,如下步骤:
步骤1,读取输入输出文件名、成像范围、最大可用内存等参数
步骤2,成像空间进行分层:
设每个从进程最大可用内存为Memory_max兆,假设层厚度以100点数为最小单位,如果有整数Nt_layer满足:
4×Nt_layer×100×Nx×Ny≤Memory_max×1024×1024(1)
4×(Nt_layer+1)×100×Nx×Ny>Memory_max×1024×1024(2)那么分层延拓的层厚度为Nt_layer×100。如果Nτ/(100×Nt_layer)为整数,则该值为延拓成像的层数,记为N_layer;如果Nτ/(100×Nt_layer)不为整数,则延拓成像的层数为其值取整加1。
Nx、Ny、Nτ依次是成像空间在三个坐标x,y,τ方向上的点数。Nt_layer等于层厚度除以100,若Nt_layer等于10,表示每层的层厚度为1000。N_layer是层数,假设成像深度τ方向共3000点,每层的厚度为1000个点,则层数N_layer=3。
步骤3,单平面波P循环,逐一对各个单平面波数据进行延拓和成像;
步骤4,单个平面波数据逐层延拓成像
(1)判断:是否是第一层:
如果是第一层,主进程从磁盘读取波场数据,利用快速傅里叶变换(FFT)将数据转换到频率空间域后,再利用并行库MPI_SEND函数将单频数据(τ=0表示地表,ω1表示一个频率)发送到从进程;如果不是第一层进入(2);
(2)主进程则利用并行库MPI_RECV函数接收从进程完成层延拓后的波场数据,再利用并行库MPI_SEND函数将单频数据(τ=τ1表示某延拓后的深度)发送到从进程;
(3)从进程接收主进程发送来的单频数据,按有限差分解进行波场延拓,再计算单频的成像值Iω1(x,y,τ);
(4)单层延拓和成像结束后,利用并行库MPI_REDUCE函数规约累加所有从进程的单频数据成像结果,并计入磁盘;
步骤5,判断是否层循环结束,即各层全部延拓成像是否结束,如果否转入步骤4,如果是转入步骤6;
步骤6,进行下一个平面波数据偏移,判断是否单面波循环结束,否转入步骤3,是转入步骤7;
步骤7,结束,输出波动方程偏移成像数据体结果。
其中,MPI_SEND是并行库MPI(MessagePassingInterface)中定义的一个消息发送函数,用于将一个消息(变量、数组等)从当前进程发送至目标进程。MPI_RECV是并行库MPI(MessagePassingInterface)中定义的一个消息接收函数,用于将指定源进程的一个消息接收到当前进程。MPI_REDUCE是并行库MPI(MessagePassingInterface)中定义的一个全局归约操作函数,用于将所有进程中的同一数据按设定的算术运算符(加减等)归约到一个指定进程。所述在步骤4的(3)中,
有限差分解进行波场延拓和成像采用:
(31)三维偏移距平面波方程:
其中,x和y是平面坐标,τ是时间深度坐标,ω是圆频率,是表示平面波的射线参数,为频率空间域波场,v=v(x,y,τ)为时间域介质层速度。
(32)方程(3)可以分裂为以下两个方程:
(33)方程(4)的解析解为:
其中,表示向下延拓步长Δτ后的波场。
(34)方程(5)可以表示如下x、y两个方向的差分方程,采用追赶法求解本方程,即实现波场由到延拓。
其中: I=(0,1,0),T=(-1,2,-1)。Δx和Δy表示x、y方向的差分间隔。α、a、b为差分方程系数。
(35)单平面偏移成像值可由频率空间域波场表示为:
real表示取实部。
本发明有效降低了波动方程叠前偏移中的内存需求,若总延拓成像深度点数为Nall,分层延拓中的层厚度为Nt_layer×100,则主进程内存需求减小到非分层整体延拓成像的(Nt_layer×100+Nω)/(Nall+Nω)倍,从进程内存需求减小到非分层整体延拓成像的Nt_layer×100/Nall倍。通过减小从进程内存需求,可以利用节点中更多的核参与偏移计算,提高了并行节点的有效利用率。本发明与非分层整体延拓成像相比,在基本不增加额外运行时间的条件下(仅增加部分节点通信时间),使得波动方程叠前时间偏移具备了大规模地震数据处理能力。
附图说明
图1是本发明方法的流程图;
图2是本发明实施例1的效果输出图
图3是本发明实施例2的效果输出图
将结合具体的实施方式加以说明
具体实施方式
本发明中,
(1)平面波波动方程叠前偏移方法
由时间域双平方根算子的频散关系可推导出
其中,kτ为双程旅行时对应的波数,γ为炮检点射线的半张角,v为常速度。上式在频率-空间域可以表示为:
其中,v=v(x,y,τ),
γ与平面波矢量有如下关系式:
进而得三维偏移距平面波方程:
对上式利用有限差分算法可进行单平面波的波场延拓。
利用零时间成像条件,单平面偏移成像值可用频率空间域波场表示为:
real表示取实部。
逐一对每一个单平面波进行延拓和成像,可以实现平面波波动方程偏移成像。
(2)分层延拓层厚度计算方法
波动方程偏移中,主进程负责从硬盘读写数据和向从进程分发作业,内存数组包括一个Nω×Nx×Ny的复数类型的三维数组和两个Nτ×Nx×Ny的实数类型的三维数组,这两个实数数组其中一个存放单频成像数据体,另一个存放最终成像数据体;从进程负责对接收自主进程的单频数据进行延拓处理,内存数组是一个Nτ×Nx×Ny的实数类型的三维数组。
在计算机集群运行波动方程偏移时,一般主进程单独使用一个节点,可用内存空间较大,多个从进程共用一个节点,共享一个节点的内存空间,因此利用从进程最大可用内存进行分层控制。记每个从进程最大可用内存为Memory_max兆,假设层厚度以100点数为最小单位,如果有整数Nt_layer满足:
4×Nt_layer×100×Nx×Ny≤Memory_max×1024×1024(6)
4×(Nt_layer+1)×100×Nx×Ny>Memory_max×1024×1024(7)
那么分层延拓的层厚度为Nt_layer×100。如果Nτ/(100×Nt_layer)为整数,则该值为延拓成像的层数,记为N_layer;如果Nτ/(100×Nt_layer)不为整数,则延拓成像的层数为其值取整加1。
如图1所述的方法包括:
步骤1,读取输入输出文件名、成像范围、最大可用内存等参数
步骤2,成像空间进行分层:
设每个从进程最大可用内存为Memory_max兆,假设层厚度以100点数为最小单位,如果有整数Nt_layer满足:
4×Nt_layer×100×Nx×Ny≤Memory_max×1024×1024(1)
4×(Nt_layer+1)×100×Nx×Ny>Memory_max×1024×1024(2)
那么分层延拓的层厚度为Nt_layer×100。如果Nτ/(100×Nt_layer)为整数,则该值为延拓成像的层数,记为N_layer;如果Nτ/(100×Nt_layer)不为整数,则延拓成像的层数为其值取整加1。
Nx、Ny、Nτ依次是成像空间在三个坐标x,y,τ方向上的点数。Nt_layer等于层厚度除以100,若Nt_layer等于10,表示每层的层厚度为1000。N_layer是层数,假设成像深度τ方向共3000点,每层的厚度为1000个点,则层数N_layer=3。
步骤3,单平面波P循环,逐一对各个单平面波数据进行延拓和成像;
步骤4,单个平面波数据逐层延拓成像
(1)判断:是否是第一层:
如果是第一层,主进程从磁盘读取波场数据,利用快速傅里叶变换(FFT)将数据转换到频率空间域后,再利用并行库MPI_SEND函数将单频数据(τ=0表示地表,ω1表示一个频率)发送到从进程;如果不是第一层进入(2);
(2)主进程则利用并行库MPI_RECV函数接收从进程完成层延拓后的波场数据,再利用并行库MPI_SEND函数将单频数据(τ=τ1表示某延拓后的深度)发送到从进程;
(3)从进程接收主进程发送来的单频数据,按有限差分解进行波场延拓,再计算单频的成像值Iω1(x,y,τ);
(4)单层延拓和成像结束后,利用并行库MPI_REDUCE函数规约累加所有从进程的单频数据成像结果,并计入磁盘;
步骤5,判断是否层循环结束,即各层全部延拓成像是否结束,如果否转入步骤4,如果是转入步骤6;
步骤6,进行下一个平面波数据偏移,判断是否单面波循环结束,否转入步骤3,是转入步骤7;
步骤7,结束,输出波动方程偏移成像数据体结果。
其中,MPI_SEND是并行库MPI(MessagePassingInterface)中定义的一个消息发送函数,用于将一个消息(变量、数组等)从当前进程发送至目标进程。MPI_RECV是并行库MPI(MessagePassingInterface)中定义的一个消息接收函数,用于将指定源进程的一个消息接收到当前进程。MPI_REDUCE是并行库MPI(MessagePassingInterface)中定义的一个全局归约操作函数,用于将所有进程中的同一数据按设定的算术运算符(加减等)归约到一个指定进程。
所述在步骤4的(3)中,
有限差分解进行波场延拓和成像采用:
(31)三维偏移距平面波方程:
其中,x和y是平面坐标,τ是时间深度坐标,ω是圆频率,是表示平面波的射线参数,为频率空间域波场,v=v(x,y,τ)为时间域介质层速度。
(32)方程(3)可以分裂为以下两个方程:
(33)方程(4)的解析解为:
其中,表示向下延拓步长Δτ后的波场。
(34)方程(5)可以表示如下x、y两个方向的差分方程,采用追赶法求解本方程,即实现波场由到延拓。
其中: I=(0,1,0),T=(-1,2,-1)。Δx和Δy表示x、y方向的差分间隔。α、a、b为差分方程系数。
(35)单平面偏移成像值可由频率空间域波场表示为:
real表示取实部。
图2和图3分别表示实施例1和实施例2
其中实施例1,南方某资料满覆盖面积为54.17km2,面元大小为25m×25m,采样时间为6s,采样间隔为2ms,采样点数为3000,CMP道集数据量为70G。采用非分层整体延拓成像方法,主进程内存需求为2457M,每个从进程的内存需求为1006M。采用本发明的分层延拓成像方法,将从进程最大可用内存空间设置为350M,偏移程序将成像深度分为3层,每层厚度为1000个采样点。分层延拓成像结果正确,图2是该资料利用本发明获得的一条线的成像剖面。
实施例2,西部某高密度资料,面元大小为7.5m×7.5m,采样时间为7s,采样间隔为2ms,采样点数为3500。对面积30km2的部分地震数据做波动方程叠前偏移处理,处理范围是360(线)×1450(CDP),成像深度为3000个点数。采用非分层整体延拓成像方法,主进程内存需求为12.5G,每个从进程的内存需求为5.83G,由于内存需求很大,一个8核16G内存的计算节点只能用2核进行偏移计算,并行节点利用率很低。采用本发明的分层延拓成像方法,为使一个节点的8核全部参与计算,将从进程最大可用内存空间设置为1.6G,则偏移程序将成像深度分为4层,第1到第3层的厚度为800个点,第4层厚度为600个点。该资料分层延拓成像结果见图3,同样正确。
Claims (3)
1.一种波动方程叠前偏移中的分层延拓成像方法,其特征在于,所述延拓成像方法先根据节点内存大小,对成像空间进行分层;后从首个单平面波数据开始,依次对每一个单平面波数据进行延拓和成像操作;在逐层延拓成像操作中包括主进程和从进程,所述主进程负责读写地震数据和向从进程分发数据,从进程负责单层单频数据的延拓成像;各单层计算结束后,从进程把延拓至本层深度后的单频数据送回到主进程;由主进程再把数据发送出去;并进行下一层的延拓成像操作,直至单平面波全部各层的单频数据延拓成像完成;所有单平面波延拓成像结束输出成像结果。
2.根据权利要求1所述的一种波动方程叠前偏移中的分层延拓成像方法,其特征在于,所述方法包括如下步骤:
步骤1,读取输入输出文件名、成像范围、最大可用内存参数;
步骤2,成像空间进行分层:
设每个从进程最大可用内存为Memory_max兆,假设层厚度以100点数为最小单位,如果有整数Nt_layer满足:
4×Nt_layer×100×Nx×Ny≤Memory_max×1024×1024(1)
4×(Nt_layer+1)×100×Nx×Ny>Memory_max×1024×1024(2)
那么分层延拓的层厚度为Nt_layer×100;如果Nτ/(100×Nt_layer)为整数,则该值为延拓成像的层数,记为N_layer;如果Nτ/(100×Nt_layer)不为整数,则延拓成像的层数为其值取整加1;
其中,Nx、Ny、Nτ依次是成像空间在三个坐标x,y,τ方向上的点数;Nt_layer等于层厚度除以100,若Nt_layer等于10,表示每层的层厚度为1000;N_layer是层数,假设成像深度τ方向共3000点,每层的厚度为1000个点,则层数N_layer=3;
步骤3,单平面波P循环,逐一对各个单平面波数据进行延拓和成像;
步骤4,单个平面波数据逐层延拓成像
(1)判断:是否是第一层:
如果是第一层,主进程从磁盘读取波场数据,利用快速傅里叶变换(FFT)将数据转换到频率空间域后,再利用并行库MPI_SEND函数将单频数据τ=0表示地表,ω1表示一个频率,发送到从进程;如果不是第一层进入(2);
(2)主进程则利用并行库MPI_RECV函数接收从进程完成层延拓后的波场数据,再利用并行库MPI_SEND函数将单频数据发送到从进程,τ=τ1表示某延拓后的深度;
(3)从进程接收主进程发送来的单频数据,按有限差分解进行波场延拓,再计算单频的成像值Iω1(x,y,τ);
(4)单层延拓和成像结束后,利用并行库MPI_REDUCE函数归约累加所有从进程的单频数据成像结果,并计入磁盘;
步骤5,判断是否层循环结束,即各层全部延拓成像是否结束,如果否转入步骤4,如果是转入步骤6;
步骤6,进行下一个平面波数据偏移,判断是否单面波循环结束,否转入步骤3,是转入步骤7;
步骤7,结束,输出波动方程偏移成像数据体结果;
其中,MPI_SEND是并行库MPI(MessagePassingInterface)中定义的一个消息发送函数,用于将一个消息从当前进程发送至目标进程;MPI_RECV是并行库MPI(MessagePassingInterface)中定义的一个消息接收函数,用于将指定源进程的一个消息接收到当前进程;MPI_REDUCE是并行库MPI(MessagePassingInterface)中定义的一个全局归约操作函数,用于将所有进程中的同一数据按设定的算术运算符归约到一个指定进程。
3.根据权利要求2所述的一种波动方程叠前偏移中的分层延拓成像方法,其特征在于,在所述步骤4的(3)中,
有限差分解进行波场延拓和成像采用:
(31)三维偏移距平面波方程:
其中,x和y是平面坐标,τ是时间深度坐标,ω是圆频率,是表示平面波的射线参数,为频率空间域波场,v=v(x,y,τ)为时间域介质层速度;
(32)方程(3)分裂为以下两个方程:
(33)方程(4)的解析解为:
(34)方程(5)表示如下x、y两个方向的差分方程,采用追赶法求解本方程,即实现波场由到的向下延拓;
其中: I=(0,1,0),T=(-1,2,-1);Δx和Δy表示x、y方向的差分间隔;α、a、b为差分方程系数;
(35)单平面偏移成像值由频率空间域波场表示为:
real表示取实部,逐一对每一个单平面波进行延拓和成像,实现平面波波动方程偏移成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110317615.3A CN103064110B (zh) | 2011-10-18 | 2011-10-18 | 一种波动方程叠前偏移中的分层延拓成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110317615.3A CN103064110B (zh) | 2011-10-18 | 2011-10-18 | 一种波动方程叠前偏移中的分层延拓成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103064110A CN103064110A (zh) | 2013-04-24 |
CN103064110B true CN103064110B (zh) | 2015-11-18 |
Family
ID=48106812
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110317615.3A Active CN103064110B (zh) | 2011-10-18 | 2011-10-18 | 一种波动方程叠前偏移中的分层延拓成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103064110B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104570124B (zh) * | 2013-10-29 | 2018-04-06 | 中国石油化工股份有限公司 | 一种适合井间地震大角度反射条件的延拓成像方法 |
CN104216012A (zh) * | 2014-08-27 | 2014-12-17 | 中国科学院地质与地球物理研究所 | 三维Born-Kirchhoff变步长插值成像方法 |
CN110895350B (zh) * | 2018-09-13 | 2021-12-24 | 中国石油化工股份有限公司 | 基于gpu的双程波傅里叶有限差分波场传播方法 |
CN109471173B (zh) * | 2018-10-08 | 2021-04-30 | 中国石油天然气集团有限公司 | 一种剩余静校正方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797038A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
CN101086533A (zh) * | 2007-07-06 | 2007-12-12 | 福州华虹智能科技开发有限公司 | 基于嵌入式Linux的浅层地震勘探仪器软件体系结构的装置 |
CN101839998A (zh) * | 2009-03-18 | 2010-09-22 | 中国石油天然气集团公司 | 一种高精度的叠前深度偏移方法 |
-
2011
- 2011-10-18 CN CN201110317615.3A patent/CN103064110B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1797038A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的叠前深度偏移方法 |
CN101086533A (zh) * | 2007-07-06 | 2007-12-12 | 福州华虹智能科技开发有限公司 | 基于嵌入式Linux的浅层地震勘探仪器软件体系结构的装置 |
CN101839998A (zh) * | 2009-03-18 | 2010-09-22 | 中国石油天然气集团公司 | 一种高精度的叠前深度偏移方法 |
Non-Patent Citations (2)
Title |
---|
基于单程波真振幅分步傅里叶叠前深度偏移方法;解建建 等;《物探化探计算机技术》;20110131;第33卷(第1期);36-39 * |
高精度克希霍夫三维叠前深度偏移及并行实现;陈志德 等;《大庆石油地质与开发》;20010630;第20卷(第3期);64-66 * |
Also Published As
Publication number | Publication date |
---|---|
CN103064110A (zh) | 2013-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tromp et al. | Spectral-element and adjoint methods in seismology | |
WO2020029015A1 (zh) | 一种人工源面波勘探方法、面波勘探装置及终端设备 | |
CN107479092A (zh) | 一种基于方向导数的频率域高阶声波方程正演模拟方法 | |
CN103064110B (zh) | 一种波动方程叠前偏移中的分层延拓成像方法 | |
US9063248B2 (en) | RTM seismic imaging using combined shot data | |
WO2004034088A2 (en) | Method and system for limited frequency seismic imaging | |
CN113063490B (zh) | 一种基于声压和质点振速双面测量的声场分离方法 | |
CN107894613A (zh) | 弹性波矢量成像方法、装置、存储介质及设备 | |
CN106168942B (zh) | 一种基于奇异边界法的波动类型动态数据重构方法 | |
CN103698809B (zh) | 一种无加速比瓶颈的克希霍夫叠前时间偏移并行方法 | |
CN110109177B (zh) | 基于旋转时空双变网格有限差分法的地震波正演模拟方法 | |
US20050027455A1 (en) | Systems and methods for geophysical imaging using amorphous computational processing | |
Gesemann | From particle tracks to velocity and acceleration fields using B-splines and penalties | |
AU2017366758A1 (en) | Diving wave illumination using migration gathers | |
CN109001813A (zh) | 一种压制多次波的方法、装置及系统 | |
CN106291678A (zh) | 一种地震数据采集方法及系统 | |
Chen et al. | pSIN: A scalable, Parallel algorithm for Seismic INterferometry of large-N ambient-noise data | |
CN102830431B (zh) | 真地表射线追踪自适应插值方法 | |
CN103675895A (zh) | 一种利用gpu提高波场延拓计算效率的方法 | |
CN105974471B (zh) | 一种基于异步流的地震数据多gpu快速正演计算方法 | |
CN109709602B (zh) | 一种远探测声波偏移成像方法、装置及系统 | |
CN116520418A (zh) | 一种弹性波角度域共成像点道集高效提取方法 | |
CN105738949B (zh) | 一种用于时移地震的九面元一致性并行处理方法 | |
CN104750954A (zh) | 一种在复杂各向异性介质中模拟地震波的方法及装置 | |
CN107976714B (zh) | 一种复杂地表分级空间距离加权的超道集计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |