CN112558150A - 一种高效的频率空间域粘声介质高斯束正演*** - Google Patents
一种高效的频率空间域粘声介质高斯束正演*** Download PDFInfo
- Publication number
- CN112558150A CN112558150A CN202011372486.3A CN202011372486A CN112558150A CN 112558150 A CN112558150 A CN 112558150A CN 202011372486 A CN202011372486 A CN 202011372486A CN 112558150 A CN112558150 A CN 112558150A
- Authority
- CN
- China
- Prior art keywords
- frequency
- gaussian beam
- ray
- unit
- forward modeling
- 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.)
- Pending
Links
- 238000004613 tight binding model Methods 0.000 claims abstract description 71
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 238000001514 detection method Methods 0.000 claims abstract description 23
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 16
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims description 22
- 239000000523 sample Substances 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 12
- 125000004122 cyclic group Chemical group 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000010521 absorption reaction Methods 0.000 description 7
- 239000002356 single layer Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000010410 layer Substances 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 235000021384 green leafy vegetables Nutrition 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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/282—Application of seismic models, synthetic seismograms
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
本发明提供一种高效的频率空间域粘声介质高斯束正演***，包括以下模块：粘声介质中高斯束表征的格林函数生成模块和基于高斯束格林函数的波场合成模块；所述格林函数生成模块，用于生成粘声介质中高斯束的格林函数；格林函数生成模块的输出结果预先存储至内存中；所述波场合成模块，基于所述粘声介质中高斯束的格林函数，实现高斯束地震波的正演；波场合成模块采用OpenMP共享内存，基于格林函数生成模块的结果完成地震波正演。本发明的有益效果是：通过OpenMP共享内存，将所有炮点和检波点处的格林函数进行预先计算并调入内存中，避免了在多次覆盖观测***下检波点格林函数的重复计算问题，切实提高了粘声介质中高斯束正演的实现效率。
Description
技术领域
本发明涉及地震勘探领域，尤其涉及一种高效的频率空间域粘声介质高斯束正演***。
背景技术
高斯束是横向振幅随高斯函数衰减的光束，通常用来描述单色电磁波辐射，也被广泛应用于油气勘探进行地震波传播。Cerveny等人在1982年在《Geophysical JournalInternational》上发表的论文“Computation of wave fields in inhomogeneous media-Gaussian beam approach”中推导了高斯束正演的表达式，指出高斯束是射线中心坐标系中波动方程的高频近似解。Yue等人在2019年《Chinese Journal of Geophysics》上发表的论文“声波介质一次散射波场高斯束Born正演”中先是利用高斯束合成束中心位置处的局部平面波，再利用逆倾斜叠加将其转化为检波点处的时空域散射波场，从而提出一种基于高斯束的一阶散射波场Born正演方法。
经典的声介质中地震波正演通常忽略了地下介质对地震波的吸收衰减，基于粘声介质进行地震波传播能更好的模拟实际地震记录。Zhu等人在2013年《GeophysicalProspecting》上发表“Approximating constant-Q seismic propagation in the timedomain”在时间域中用波动方程方法来描述粘声介质中地震波传播的吸收衰减。但波动方程正演计算量巨大，不如射线类正演方法灵活高效，而高斯束正演既继承了射线类正演方法的灵活高效的优势，也可以像波动方程正演得到高精度的地震数据。Yue等人在2020年《Geophysics》文章“Gaussian Beam Born Modeling for Single-Scattering Waves inVisco-Acoustic Media”中推导粘声介质中的高斯束正演的公式，准确模拟粘声介质中地震波场的传播。
与声介质中高斯束正演不同，粘声介质中衰减是与频率有关，因此，需要将时空域转化到频率域，在频率域中计算格林函数。由于在多次覆盖观测***下，检波点格林函数会出现重复计算的问题，使计算量大大增加。
发明内容
有鉴于此，本发明针对粘声介质中高斯束正演的检波点格林函数会出现重复计算，使计算量大大增加的问题，提出一种高效的频率空间域粘声介质高斯束正演***，该***将所有炮点和检波点处的格林函数进行预先计算并调入内存中，明显降低了计算量，并通过OpenMp共享内存，以多个节点为基础，利用节点多核共享解决内存需求大的问题，从而实现高效的粘声介质中的高斯束正演。
本发明提出的一种高效的频率空间域粘声介质高斯束正演***，具体包括以下模块：
粘声介质中高斯束表征的格林函数生成模块和基于高斯束格林函数的波场合成模块；
所述格林函数生成模块，用于生成粘声介质中高斯束的格林函数；
所述波场合成模块，基于所述粘声介质中高斯束的格林函数，实现高斯束地震波的正演；
所述格林函数生成模块，包括四个循环单元，分别为炮检循环单元、射线循环单元、第一频率循环单元和第一空间循环单元；所述炮检循环单元位于最外层，且每个循环单元依次内嵌，即射线循环单元内嵌于炮检循环单元、第一频率循环单元内嵌于射线循环单元、第一空间循环单元内嵌于第一频率循环单元，炮检循环单元结束，则最终生成粘声介质中高斯束的格林函数；
所述波场合成模块，包括四个循环单元，分别为炮点循环单元、检波点循环单元、第二频率循环单元和第二空间循环单元；所述炮点循环单元位于最外层，且每个循环单元依次内嵌，即检波点循环单元内嵌于所述炮点循环单元、第二频率循环单元内嵌于检波点循环单元、第二空间循环单元内嵌于第二频率循环单元；
所述炮点循环单元结束，则最终完成高斯束地震波的正演。
进一步地，所述炮检循环单元，用于决定粘声介质中高斯束的格林函数生成循环总次数，在粘声介质高斯束正演***中，共有nshot个炮检点，即发出nshot次模拟地震波。
进一步地，每次模拟地震波，产生多个频率段以及多个不同方向的射线束；频率段共有nw个，iw为其中任意一个频率段；射线束方向有np个，ip为任意一个射线束方向；射线束方向上有多个点，ix，iz为射线方向ip上的任意一个点。
进一步地，所述射线循环单元，其循环条件为计算达到每次模拟地震波的射线束方向总数np；
所述射线循环单元在一次计算过程中，对每次模拟地震波中任一方向ip的射线束进行时空域正演，若该方向射线束有正演结果，即得到该方向射线束对应的格林函数，并标记该方向，直至达到射线束方向总数np，则所述射线循环单元结束。
进一步地，所述第一频率循环单元，其循环条件为计算次数达到每个射线方向的地震频率总数nw；
所述第一频率循环单元在一次计算过程中，对一个射线束方向ip的任意一个频率iw进行频域变换，得到该射线束方向对应该频率的频域格林函数，直至达到每个射线方向的地震频率总数nw，则所述第一频率循环单元结束。
进一步地，所述第一空间循环单元，其循环条件为计算次数达到射线束方向总数np的所有点的总数nx*nz；
所述第一空间循环单元在一次计算过程中，连续对射线束方向ip的任一个点ix，iz进行该射线束方向对应频率iw的频域格林函数积分，得到射线束方向ip，频率iw对应的格林函数；
在循环计算过程中，对所有射线束方向，频率iw对应的格林函数进行叠加，得到所有射线束方向叠加后频率iw对应的格林函数Green(iz,ix,iw,ishot)；
Green(iz,ix,iw,ishot)表示一次炮检循环过程中，频率iw对应的所有射线束方向叠加的频域格林函数。
进一步地，所述炮点循环单元、检波点循环单元、第二频率循环单元和第二空间循环单元均进行OpenMP处理。
所述炮点循环单元，用于决定高斯束地震波正演过程的总循环次数，在粘声介质高斯束正演***中，共有nshot个炮点；
所述检波点循环单元，用于接收炮点发出的地震波，共有nreceiver个检波点；
所述第二频率循环单元，对炮点产生的地震波任一频率进行循环，共有nw个频率；
所述第二空间循环单元，对地震空间内的ix和iz网格点利用所述Green(iz,ix,iw,ishot)进行积分，完成高斯束地震波正演。
本发明提供的有益效果是：通过OpenMP共享内存，将所有炮点和检波点处的格林函数进行预先计算并调入内存中，避免了在多次覆盖观测***下检波点格林函数的重复计算问题，切实提高了粘声介质中高斯束正演的实现效率。
附图说明
图1是本发明高效的频率空间域粘声介质高斯束正演***的结构示意图；
图2为单层模型为例，不同介质不同品质因素的波场快照；
图3是以单层模型为例，不同介质不同品质因素下的地震道振幅对比和相位对比示意图；
图4是以saltbag模型为例，输入速度模型和Q模型的示意图；
图5是以saltbag模型为例不同介质的波长快照示意图；
图6是以saltbag为例，不同介质不同品质因素下的地震道振幅对比和相位对比示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述。
首先，对粘声介质中高斯束的格林函数表示为式(1)：
其中，x表示地下空间点，x'表示炮检点，p＝(px,py,pz)T是初始射线参数。Ac(x,x′)是粘声介质中射线束振幅，τc(x,x′)表示粘声介质中的复值旅行时，该复旅行时包含吸收衰减旅行时τQ(x,x′)，它们分别表示如下
Ac(x,x′)＝AR(x,x′)+iAI(x,x′), (2)
其中，AR和AI分别表示振幅的实部和虚部，τR和τI分别表示声介质中旅行时的实部和虚部。Q表示品质因子，是粘声介质中地震波吸收衰减的表征参数。
在公式(1)基础上，模拟地震波传播并得到地震波正演数据，高斯束Born正演公式可以表示为：
D(xr,xs,ω)＝2ω2∫Tm(x)G(x,xs,ω)G(x,xr,ω)dx (5)
其中，D(xr,xs,ω)表示地表炮检点位于(xr,xs)的地震波正演数据，m(x)表示地下模型空间，G(x,xs,ω),G(x,xr,ω)分别表示震源和检波点处的背景格林函数，x,xs,xr分别表示地下空间点,震源和检波点，T表示一系列地下散射点。
将公式(1)代入公式(5)，高斯束正演公式也可以表示为：
其中，ps＝(psx,psy,psz)T,pr＝(psx,psy,psz)T分别表示震源和检波点处的初始射线参数，Ac(x,xs),Ac(x,xr)分别表示粘声介质中震源和检波点处的射线振幅，τc(x,xs),τc(x,xr)分别表示粘声介质中震源和检波点至地下空间的复值旅行时。
与声介质中高斯束正演不同，粘声介质中衰减是与频率有关，因此，需要将时空域转化到频率域，在频率域中计算格林函数。
由于在多次覆盖观测***下，检波点格林函数会出现重复计算的问题，使计算量大大增加。
本发明针对此问题，将所有炮点和检波点处的格林函数进行预先计算并调入内存中，明显降低了计算量，与此同时，对内存的需求也增加。所以本发明实现一种通过OpenMP共享内存，以多个节点为基础，利用节点多核可以共享并且内存大的优势，解决了这一问题，从而实现高效。
请参考图1，一种高效的频率空间域粘声介质高斯束正演***，包括以下模块：
粘声介质中高斯束表征的格林函数生成模块和基于高斯束格林函数的波场合成模块；
所述格林函数生成模块，用于生成粘声介质中高斯束的格林函数；
即通过所述格林函数生成模块，将所有炮点和检波点处的格林函数预先进行计算，并存储至内存当中；
所述波场合成模块，基于所述粘声介质中高斯束的格林函数，实现高斯束地震波的正演；
由于在格林函数生成模块，将格林函数预先计算，并存储到了内存当中，使得内存占用增大，故在所述波场合成模块，采用OpenMP共享内存的方法，利用节点多核共享内存，解决内存需求问题，使得波场合成相比传统方法而言，显得高效。
所述格林函数生成模块，包括四个循环单元，分别为炮检循环单元、射线循环单元、第一频率循环单元和第一空间循环单元；所述炮检循环单元位于最外层，且每个循环单元依次内嵌，即射线循环单元内嵌于炮检循环单元、第一频率循环单元内嵌于射线循环单元、第一空间循环单元内嵌于第一频率循环单元，炮检循环单元结束，则最终生成粘声介质中高斯束的格林函数；
所述炮检循环单元，用于决定粘声介质中高斯束的格林函数生成循环总次数，在粘声介质高斯束正演***中，共有nshot个炮检点，即发出nshot次模拟地震波。
每次模拟地震波，产生多个频率段以及多个不同方向的射线束；频率段共有nw个，iw为其中任意一个频率段；射线束方向有np个，ip为任意一个射线束方向；射线束方向上有多个点，ix，iz为射线方向ip上的任意一个点。
所述射线循环单元，其循环条件为计算达到每次模拟地震波的射线束方向总数np；
所述射线循环单元在一次计算过程中，对每次模拟地震波中任一方向ip的射线束进行时空域正演，若该方向射线束有正演结果，即得到该方向射线束对应的格林函数，并标记该方向，直至达到射线束方向总数np，则所述射线循环单元结束。
所述第一频率循环单元，其循环条件为计算次数达到每个射线方向的地震频率总数nw；
所述第一频率循环单元在一次计算过程中，对一个射线束方向ip的任意一个频率iw进行频域变换，得到该射线束方向对应该频率的频域格林函数，直至达到每个射线方向的地震频率总数nw，则所述第一频率循环单元结束。
所述第一空间循环单元，其循环条件为计算次数达到射线束方向总数np的所有点的总数nx*nz；
所述第一空间循环单元在一次计算过程中，连续对射线束方向ip的任一个点ix，iz进行该射线束方向对应频率iw的频域格林函数积分，得到射线束方向ip，频率iw对应的格林函数；
在循环计算过程中，对所有射线束方向，频率iw对应的格林函数进行叠加，得到所有射线束方向叠加后频率iw对应的格林函数Green(iz,ix,iw,ishot)；
Green(iz,ix,iw,ishot)表示一次炮检循环过程中，频率iw对应的所有射线束方向叠加的频域格林函数。
优选的，关于格林函数生成模块的一个实施例具体如下：
在进行格林函数生成生成时，申请两个结构体数组，cells(nz,nx)和Green(nz,nx,nw,nshot)，分别用于存储储个方向高斯束对应的格林函数和最终生成的格林函数；
其中，z和x指空间点二维坐标，w为频率，nw是指频率的循环次数，shot指炮点，nshot是指炮点的循环次数。将cells(nz,nx)及涉及的循环相关变量和临时变量设置为私有化变量。下面是所述格林函数生成模块的四个依次内嵌的循环单元
炮检循环单元(nshot)：
即对编号范围从1到nshot多炮地震数据进行计算，共计算nshot次。
射线方向循环单元(np)：
对某一个方向ip的射线束进行时空域正演，形成该方向的cells(nz,nx),并对射线束有值位置做好标记。
考虑到粘声介质中吸收衰减的影响，根据射线追踪得到吸收衰减旅行时τQ(x,x′)，以及粘声介质中的地震波振幅。
第一频率循环单元(nw)：
对某一频率ip的有标记的射线束进行计算，得到频率域格林函数相关信息。
第一空间循环单元(nx*nz)：
对模型空间内ix和iz的网格点，进行格林函数的积分计算，叠加了，得到Green(iz,ix,iw,ishot)，同时，考虑到粘声介质中吸收衰减的影响，加入振幅衰减项和相位畸变项，即合成衰减效应。
所述波场合成模块，包括四个循环单元，分别为炮点循环单元、检波点循环单元、第二频率循环单元和第二空间循环单元；所述炮点循环单元位于最外层，且每个循环单元依次内嵌，即检波点循环单元内嵌于所述炮点循环单元、第二频率循环单元内嵌于检波点循环单元、第二空间循环单元内嵌于第二频率循环单元；
所述炮点循环单元，用于决定高斯束地震波正演过程的总循环次数，在粘声介质高斯束正演***中，共有nshot个炮点；
所述检波点循环单元，用于接收炮点发出的地震波，共有nreceiver个检波点；
所述第二频率循环单元，对炮点产生的地震波任一频率进行循环，共有nw个频率；
所述第二空间循环单元，对地震空间内的ix和iz网格点利用所述Green(iz,ix,iw,ishot)进行积分，完成高斯束地震波正演
所述炮点循环单元结束，则最终完成高斯束地震波的正演。
优选的，关于波场合成模块的一个实施例具体如下：
在波场合成模块中，对模块内的各个循环单元进行OpenMP处理，将涉及的循环相关变量和临时变量设置为私有化变量，共享格林函数数组Green(nz,nx,nw,nshot)。
炮点循环单元(nshot)：
即对编号范围从1到nshot多炮地震数据进行计算，共计算nshot次。
检波点循环单元(nreceiver)
对编号范围从1到nreceiver多道地震数据进行计算，共计算nreceiver次，即对任意检测地震波的检波器ireceiver进行nreceiver次循环。
第二频率循环单元(nw)
对地震波有效频带范围内任一频率iw进行nw次循环
第二空间循环单元(nx*nz)
对模型空间内ix和iz的网格点利用已知格林函数Green(nz,nx,nw,nshot)进行积分，至此完成了对格林函数Green(nz,nx,nw,nshot)的高斯束地震波正演。
为了更直观解释本发明的技术效果，下面以单层模型和saltbag模型为例进行地震波正演。
(1)单层模型
以单层模型为例，地震波在地下介质中速度为3000m/s时，输入单炮，每炮301道，道间距10m，采样点数为2000，时间采样间隔为0.001s。在声介质以及粘声介质Q＝10，30，70的情况下输出0.33s时刻的波场图如图2所示；图2(a)声介质中的波场快照，其中炮点位于(2256,0)；(b)粘声介质Q＝10中的波场快照，其中炮点位于(2256,0)；(c)粘声介质Q＝30中的波场快照，其中炮点位于(2256,0)；(d)粘声介质Q＝70中的波场快照，其中炮点位于(2256,0)。
为了便于更清楚观察图2(a)、(b)、(c)、(d)之间的变化，提取图(1)中震源附近第2256处地震道振幅、相位信息作对比，请参考图3；图3在声介质以及粘声介质Q＝10，30，70下第2256处地震道(a)振幅对比，(b)相位对比。
由图可知，粘声介质中地震波会发生振幅衰减和相位畸变的现象，并且随着品质因子Q越小，振幅衰减和相位畸变越严重。
(2)saltbag模型
以saltbag模型为例，输入单炮，每炮1251道，道间距8m，采样点数为2000，时间采样间隔为0.001s。输入的速度模型和Q模型如图4(a)、(b)所示，在声介质和粘声介质中输出0.4s时刻的波场如图5(a)、(b)所示。
图4中，(a)为输入的速度模型，(b)输入的Q模型；
图5(a)声介质中的波场快照，其中炮点位于(8945,0)；图5(b)粘声介质中的波场快照，其中炮点位于(8945,0)。
便于清楚观察理想情况下和吸收衰减情况下的波场之间的区别，提取震源附近第8895道地震道的振幅和相位信息作对比，如图6(a)、(b)所示。
由图可知，粘声介质中地震波相对于声介质中地震波会发生振幅衰减和相位畸变的现象。
本发明提供的有益效果是：通过OpenMP共享内存，将所有炮点和检波点处的格林函数进行预先计算并调入内存中，避免了在多次覆盖观测***下检波点格林函数的重复计算问题，切实提高了粘声介质中高斯束正演的实现效率。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种高效的频率空间域粘声介质高斯束正演***，其特征在于：具体包括以下模块：
粘声介质中高斯束表征的格林函数生成模块和基于高斯束格林函数的波场合成模块；
所述格林函数生成模块，用于生成粘声介质中高斯束的格林函数；
所述波场合成模块，基于所述粘声介质中高斯束的格林函数，实现高斯束地震波的正演；
所述格林函数生成模块，包括四个循环单元，分别为炮检循环单元、射线循环单元、第一频率循环单元和第一空间循环单元；所述炮检循环单元位于最外层，且每个循环单元依次内嵌，即射线循环单元内嵌于炮检循环单元、第一频率循环单元内嵌于射线循环单元、第一空间循环单元内嵌于第一频率循环单元，炮检循环单元结束，则最终生成粘声介质中高斯束的格林函数；
所述波场合成模块，包括四个循环单元，分别为炮点循环单元、检波点循环单元、第二频率循环单元和第二空间循环单元；所述炮点循环单元位于最外层，且每个循环单元依次内嵌，即检波点循环单元内嵌于所述炮点循环单元、第二频率循环单元内嵌于检波点循环单元、第二空间循环单元内嵌于第二频率循环单元；
所述炮点循环单元结束，则最终完成高斯束地震波的正演。
2.如权利要求1所述的高效的频率空间域粘声介质高斯束正演***，其特征在于：所述炮检循环单元，用于决定粘声介质中高斯束的格林函数生成循环总次数，在粘声介质高斯束正演***中，共有nshot个炮检点，即发出nshot次模拟地震波。
3.如权利要求2所述的高效的频率空间域粘声介质高斯束正演***，其特征在于：每次模拟地震波，产生多个频率段以及多个不同方向的射线束；频率段共有nw个，iw为其中任意一个频率段；射线束方向有np个，ip为任意一个射线束方向；射线束方向上有多个点，ix，iz为射线方向ip上的任意一个点。
4.如权利要求3所述的高效的频率空间域粘声介质高斯束正演***，其特征在于：所述射线循环单元，其循环条件为计算达到每次模拟地震波的射线束方向总数np；
所述射线循环单元在一次计算过程中，对每次模拟地震波中任一方向ip的射线束进行时空域正演，若该方向射线束有正演结果，即得到该方向射线束对应的格林函数，并标记该方向，直至达到射线束方向总数np，则所述射线循环单元结束。
5.如权利要求4所述的高效的频率空间域粘声介质高斯束正演***，其特征在于：所述第一频率循环单元，其循环条件为计算次数达到每个射线方向的地震频率总数nw；
所述第一频率循环单元在一次计算过程中，对一个射线束方向ip的任意一个频率iw进行频域变换，得到该射线束方向对应该频率的频域格林函数，直至达到每个射线方向的地震频率总数nw，则所述频率循环单元结束。
6.如权利要求5所述的高效的频率空间域粘声介质高斯束正演***，其特征在于：
所述第一空间循环单元，其循环条件为计算次数达到射线束方向总数np的所有点的总数nx*nz；
所述第一空间循环单元在一次计算过程中，连续对射线束方向ip的任一个点ix，iz进行该射线束方向对应频率iw的频域格林函数积分，得到射线束方向ip，频率iw对应的格林函数；
在循环计算过程中，对所有射线束方向，频率iw对应的格林函数进行叠加，得到所有射线束方向叠加后频率iw对应的格林函数Green(iz,ix,iw,ishot)；
Green(iz,ix,iw,ishot)表示一次炮检循环过程中，频率iw对应的所有射线束方向叠加的频域格林函数。
7.如权利要求1所述的高效的频率空间域粘声介质高斯束正演***，其特征在于：
所述炮点循环单元、检波点循环单元、第二频率循环单元和第二空间循环单元均进行OpenMP处理。
8.如权利要求6所述的高效的频率空间域粘声介质高斯束正演***，其特征在于：
所述炮点循环单元，用于决定高斯束地震波正演过程的总循环次数，在粘声介质高斯束正演***中，共有nshot个炮点；
所述检波点循环单元，用于接收炮点发出的地震波，共有nreceiver个检波点；
所述第二频率循环单元，对炮点产生的地震波任一频率进行循环，共有nw个频率；
所述第二空间循环单元，对ix、iz利用所述Green(iz,ix,iw,ishot)进行积分，完成高斯束地震波正演。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011372486.3A CN112558150A (zh) | 2020-11-30 | 2020-11-30 | 一种高效的频率空间域粘声介质高斯束正演*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011372486.3A CN112558150A (zh) | 2020-11-30 | 2020-11-30 | 一种高效的频率空间域粘声介质高斯束正演*** |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112558150A true CN112558150A (zh) | 2021-03-26 |
Family
ID=75046772
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011372486.3A Pending CN112558150A (zh) | 2020-11-30 | 2020-11-30 | 一种高效的频率空间域粘声介质高斯束正演*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112558150A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113466934A (zh) * | 2021-06-17 | 2021-10-01 | 中国地质大学（武汉） | 基于点扩散函数的Hessian矩阵的计算与存储方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103675894A (zh) * | 2013-12-24 | 2014-03-26 | 中国海洋石油总公司 | 一种基于三维高斯束射线追踪和频率域合成地震记录方法 |
| CN110927784A (zh) * | 2019-11-29 | 2020-03-27 | 中国地质大学（武汉） | 一种基于OpenMP的频率域高斯束格林函数计算方法 |
| CN111694051A (zh) * | 2020-06-17 | 2020-09-22 | 电子科技大学 | 一种基于高斯束的黏声介质地震波正演方法 |
-
2020
- 2020-11-30 CN CN202011372486.3A patent/CN112558150A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103675894A (zh) * | 2013-12-24 | 2014-03-26 | 中国海洋石油总公司 | 一种基于三维高斯束射线追踪和频率域合成地震记录方法 |
| CN110927784A (zh) * | 2019-11-29 | 2020-03-27 | 中国地质大学（武汉） | 一种基于OpenMP的频率域高斯束格林函数计算方法 |
| CN111694051A (zh) * | 2020-06-17 | 2020-09-22 | 电子科技大学 | 一种基于高斯束的黏声介质地震波正演方法 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113466934A (zh) * | 2021-06-17 | 2021-10-01 | 中国地质大学（武汉） | 基于点扩散函数的Hessian矩阵的计算与存储方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9632192B2 (en) | Method of processing seismic data by providing surface offset common image gathers | |
| CN108802813B (zh) | 一种多分量地震资料偏移成像方法及*** | |
| CN103105623B (zh) | 一种地震勘探中的数据波形处理方法 | |
| CN107272058A (zh) | 成像方法、成像装置以及计算机存储介质 | |
| CN110927784B (zh) | 一种基于OpenMP的频率域高斯束格林函数计算方法 | |
| CN114839673B (zh) | 多震源高效采集波场分离方法、分离***及计算机设备 | |
| CN109633749A (zh) | 基于散射积分法的非线性菲涅尔体地震走时层析成像方法 | |
| KR101820850B1 (ko) | 직접 파형 역산의 반복 적용을 이용한 탄성파 영상화 장치 및 방법 | |
| Rickett et al. | Passive seismic imaging applied to synthetic data | |
| CN112558150A (zh) | 一种高效的频率空间域粘声介质高斯束正演*** | |
| CN103823241A (zh) | 一种在炮检距域中求取折射剩余静校正量的方法 | |
| CN109738944B (zh) | 基于广角反射的地震采集参数确定方法及装置 | |
| CN116719086B (zh) | 基于点扩散函数的稀疏海底四分量数据高分辨率成像方法 | |
| CN109581484B (zh) | 地震纵波分量多次波总照明度指标分析方法及装置 | |
| CN114002742B (zh) | 一种欧拉高斯束偏移成像方法及装置 | |
| CN113866821B (zh) | 一种基于照明方向约束的被动源干涉偏移成像方法和*** | |
| CN115421195A (zh) | 地震探测中速度场的生成方法、装置、设备及存储介质 | |
| CN106556860B (zh) | 布设vsp观测***的方法和装置 | |
| Reshef | Prestack depth imaging of three-dimensional shot gathers | |
| CN106597536A (zh) | 一种地震成像的处理方法及装置 | |
| CN110873894A (zh) | 基于高斯束反偏移的炮记录获取方法及*** | |
| Tang et al. | Subsalt imaging by target-oriented inversion-based imaging: A 3D field-data example | |
| CN111665546A (zh) | 用于可燃冰探测的声学参数获取方法 | |
| Hoang | Random-shot full-waveform inversion of shallow seismic surface waves | |
| CN111999769B (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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210326 |