CN111781639B - 针对obs多分量数据的炮检互易弹性波全波形反演方法 - Google Patents
针对obs多分量数据的炮检互易弹性波全波形反演方法 Download PDFInfo
- Publication number
- CN111781639B CN111781639B CN202010499751.8A CN202010499751A CN111781639B CN 111781639 B CN111781639 B CN 111781639B CN 202010499751 A CN202010499751 A CN 202010499751A CN 111781639 B CN111781639 B CN 111781639B
- Authority
- CN
- China
- Prior art keywords
- shot
- component
- data
- wave
- inversion
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 238000004587 chromatography analysis Methods 0.000 claims abstract description 6
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 239000011435 rock Substances 0.000 claims abstract description 4
- 238000003384 imaging method Methods 0.000 claims abstract description 3
- 238000007689 inspection Methods 0.000 claims abstract description 3
- 238000007670 refining Methods 0.000 claims abstract description 3
- 238000006073 displacement reaction Methods 0.000 claims description 36
- 238000001514 detection method Methods 0.000 claims description 22
- 230000005284 excitation Effects 0.000 claims description 13
- 239000000523 sample Substances 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 5
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 2
- 229910052709 silver Inorganic materials 0.000 description 2
- 239000004332 silver Substances 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- 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
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，其包括以下步骤：1)对原始地震数据进行预处理；2)通过初至层析获取初始纵波速度模型，并通过声波波形反演细化纵波速度模型；3)以初至层析结果或声波波形反演结果，根据岩石物理关系建立横波以及密度的初始模型，并设定反演参数，包括反演频率范围和间隔；4)根据波形和初始弹性参数进行炮检互易弹性波全波形反演，最终获取地下纵横波速度模型并成像。与现有技术相比，本发明具有计算效率高且不影响反演精度、实施简单、可移植性强等优点。
Description
技术领域
本发明涉及海底地震仪多分量数据弹性波全波形反演领域，尤其是涉及一种针对OBS多分量数据的炮检互易弹性波全波形反演方法。
背景技术
近年来，随着海洋勘探的兴起，海底地震仪(OBS)凭借着其接收到的信噪比更高、孔径更宽、更高质量的多分量数据，针对OBS多分量数据的弹性波波形反演的实际应用逐渐引起学者的关注，在实际应用中，OBS观测系统的震源(气枪震源)的数量远大于检波器的数量，炮检数比高达几十倍(二维勘探)至几百倍(三维)，使用常规的基于伴随状态法的全波形反演，每轮迭代所需正演次数正比于两倍的震源数，计算量非常高。
现有的基于伴随状态法的声波全波形反演，可以直接交换炮点和检波点的位置以实现互易，但弹性波全波形反演无法依据此方式直接应用。频率域直接求解算子的方式通过复LU分解使得全波形反演中的正演模拟不依赖于震源参数(Marfurt，1984)，但是这种方式在应对实际数据时，由于内存的巨大消耗，无法适应大规模模拟。
Bunks等(Bunks et al.,1995)提出通过少数关键频率反演可以得到良好的反演结果(Sirgue et al.,2004；Operto et al.,2006)，一定程度上提高计算和存储效率，但计算量仍然巨大，Roemro等(Romero et al.,2000)进一步提出采用相位编码同时源进行反演，Krebs等(Krebs et al.,2009)提出快速相位编码全波形反演算法，将多炮进行组合以提高计算效率，但这样会引入随机噪声，影响反演效果。Li等(Li et al.,2012)和Warner等(Warner et al.,2013)提出随机选择数据子集来进行全波形反演，以减少反演中所需正演模拟炮数。zhu等(zhu et al.,2017)提出在压缩感知框架下基于稀疏约束的全波形反演方法，在减少计算量的同时获得较好的反演结果。此外，Liu等(Liu et al.,2015)提出改进的散射积分算法全波形反演，将计算梯度时所需要的正演次数从两倍的炮数减少到炮数加两倍的检波器数，针对OBS观测系统，计算效率得到提升，但是针对OBS观测系统，由于炮非常密集，计算量仍然非常大。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种针对OBS多分量数据的炮检互易弹性波全波形反演方法。
本发明的目的可以通过以下技术方案来实现：
一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，包括以下步骤：
1)对原始地震数据进行预处理；
2)通过初至层析获取初始纵波速度模型，并通过声波波形反演细化初始纵波速度模型；
3)根据初至层析结果或声波波形反演结果，利用岩石物理关系建立横波以及密度的初始模型，并设定反演参数，反演参数包括反演频率范围和间隔；
4)根据波形和初始弹性参数，对初始模型进行炮检互易弹性波全波形反演，获取地下纵横波速度模型并成像。
所述的步骤1)中，预处理包括重定位、多分量数据的旋转以及去噪。
所述的步骤4)具体包括以下步骤：
41)根据所需要的空间精度，建立与原始观测系统等效的炮检互易观测系统；
42)在等效的炮检互易观测系统下，分别进行以x方向单向位移震源的时间域正演以及以z方向单向位移震源的时间域正演；
43)对正演得到的波场信息分量记录进行线性叠加，判断组合得到的理论合成数据与观测数据是否匹配，并采用二范数的目标函数进行检验，若目标函数值满足要求，则将当前纵横波速度模型作为最终纵横波速度模型，若否，则执行步骤44)；
44)在检波点处分别进行x方向和z方向的单向位移震源激发，采用时间域交错网格有限差分正演，分别获取所有检波点端波场，并将x分量残差和z分量残差作为残差震源，分别在炮点端以差分格式组合反传，获取炮点端的线性组合波场；
45)通过x方向的单向位移震源产生的检波点端波场与x分量残差作为残差震源产生的炮点端的线性组合波场的互相关求取对应梯度，累加形成整个观测系统所对应的x分量梯度，并通过z方向的单向位移震源产生的检波点端波场与z分量残差作为残差震源产生的炮点端的线性组合波场的互相关求取对应梯度，累加形成整个观测系统所对应的z分量梯度；
46)计算最速下降方向、方向更新步长和模型更新量后更新当前纵横波速度模型；
47)重复步骤41-46)，直到获得匹配的最终纵横波速度模型。
所述的步骤41)中，等效的炮检互易观测系统具体为：
将震源的激发从在原始炮点位置激发改为在原始检波点位置激发:将记录的接收位置从在原始检波点位置处接收改为在原始炮点处位置接收，且在炮点及炮点邻近点处组合接收。
所述的步骤42)中，以x方向单向位移震源的时间域正演以及以z方向单向位移震源的时间域正演具体为：
针对接收的多分量数据为x和z分量记录的二维数据的情况，在等效的炮检互易观测系统下将x方向与z方向单向位移震源的正演分开单独进行，并记录各自所产生的波场信息分量记录uxx、uzx、uxz和uzz，其中，uij表示以j方向单向位移震源激发，i方向所接收到的数据分量，且j＝x或z，i＝x或z。
所述的步骤43)中，将多空间点组合接收的分量记录按照空间差分的格式进行线性叠加，即可合成与观测数据所匹配的理论合成数据，具体包括以下步骤：
431)将uxx和uzx根据以下公式进行叠加，对应于实际观测数据的x分量所匹配的数据ux(xr,zr；xs,zs)：
其中，c(k)为交错网格有限差分系数，k为空间网格上移动的计数，(xr,zr)为检波点位置；
432)将uxz和uzz根据以下公式进行叠加，对应于实际观测数据的z分量所匹配的数据uz(xr,zr；xs,zs)
433)将叠加得到的合成数据与实际数据对应分量按照二范数进行匹配，则有：
其中，E(m)表示在地下模型参数m下所计算出的目标函数值，uij(xs,zs；xr,zr)表示在(xr,zr)位置处以j方向单向位移震源激发，在(xs,zs)位置处所接收到的i方向数据分量，di(xr,zr；xs,zs)表示在(xs,zs)位置处激发应力P波源，在(xr,zr)位置处所接收到的i方向数据分量，m为地下模型参数，ns为炮点总数，nr为检波点总数。
所述的步骤44)中，检波点处激发产生的波场以及炮点端线性组合的反传波场的计算步骤为：
441)在检波点端激发x方向单向位移震源，存储波场信息分量记录uxx，uzx，将x分量观测记录的残差作为震源在(xs-1,zs)，(xs-2,zs)，…，(xs-M/2-1,zs)以及(xs,zs)，(xs+1,zs)，…，(xs+M/2,zs)处分别施加在x方向位移上，同时将x分量观测记录的残差作为震源在(xs,zs-1)，(xs,zs-2)，…，(xs,zs-M/2-1)以及(xs,zs)，(xs,zs+1)，…，(xs,zs+M/2)处施加在z方向位移上，计算波场并存储；
442)在检波点端激发z方向单向位移震源，存储波场信息分量记录uxz，uzz，将z分量观测记录的残差作为震源在(xs-1,zs)，(xs-2,zs)，…，(xs-M/2-1,zs)以及(xs,zs)，(xs+1,zs)，…，(xs+M/2,zs)处分别施加在x方向位移上，同时将z分量观测记录的残差作为震源在(xs,zs-1)，(xs,zs-2)，…，(xs,zs-M/2-1)以及(xs,zs)，(xs,zs+1)，…，(xs,zs+M/2)处施加在z方向位移上，计算波场并存储。
所述的步骤45)中，将时间域计算得到的波场通过傅里叶变换得到频率域波场。
根据以下公式计算梯度：
其中，Re为取实部操作，*为复共轭算子，T表示转置，A为频率域负阻抗矩阵，u*i表示检波点处i方向激发的单分量位移震源所产生的波场在炮点处x和z方向接收的位移分量记录，ri表示检波点处观测记录i方向的残差。
与现有技术相比，本发明具有以下优点：
一、计算效率高：本发明基于互易的观测系统完成合成数据的模拟，所需正演次数为两倍的检波器总数，突破了现有方法数据匹配时需要正演两倍于炮数的限制；使得OBS多分量数据反演的计算时间大大减少。
二、实施简单且可移植性强：本方法采样特定方式于检波点端计算正传波场，在炮点端反传残差波场，仅需对现有的基于伴随状态法的程序作简单修改就可以实现，可移植性强。
三、计算效率提高的同时不影响反演精度：本方法在提高计算效率的同时，不会影响反演结果，与传统方法的反演精度基本一致。
附图说明
图1为本发明的方法流程图。
图2为本发明的硬件结构示意图。
图3为实施例1的真实理论模型图。
图4为实施例1的初始模型图，图(4a)为纵波速度，图(4b)为横波速度。
图5为获取观测数据的观测系统。
图6为实施例1的地震记录水平分量和垂直分量图，其中，图(6a)为地表水平4km的地震记录水平分量图，图(6b)为地表水平4km的地震记录垂直分量图。
图7为实施例1的等效炮检互易观测系统示意图。
图8为实施例1的Reciprocal-FWI反演结果图，图(8a)为纵波速度反演结果，图(8b)为横波速度反演结果。
图9为实施例1的AD-FWI反演结果图，图(9a)为纵波速度反演结果，图(9b)为横波速度反演结果。
图10为实施例1的SI-FWI反演结果图，图(10a)为纵波速度反演结果，图(10b)为横波速度反演结果。
图11为实施例1的真实模型、初始模型、与Reciprocal-FWI、AD-FWI和SI-FWI在水平方向6km处的速度对比图，其中，图(11a)为纵波速度对比图，图(11b)为纵波速度对比图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
随着OBS数据在海上勘探中应用逐渐广泛，声压分量的波形反演也逐渐成熟，但OBS四分量数据的利用仍然不够充分。考虑到使用弹性波多参数全波形反演能够更精细地对地下介质参数作估计，因此将弹性波多参数全波形反演引入到OBS数据处理中，然而，相较于声波波形反演，弹性波波形反演的计算效率低是首要解决的问题。
针对OBS观测系统，由于炮数远远大于检波器数，最快速高效的处理方式是通过炮检互易来将反演所需计算量从正比于炮数降至正比于检波器数。通过使用本发明提出的针对海底地震仪多分量数据的炮检互易弹性波全波形反演方法，相比于传统的反演算法，计算效率上得到了极大地提升。
实施例1：
本实施例将以二维修改后的Overthrust模型作为真实模型(如图3所示)，该模型共有801×187个网格，网格间距为10m×10m，100m深度以上为水层，纵波速度为1500m/s，横波速度为0m/s，密度为1g/cm3(本实例中密度已知)。在该模型上进行弹性波正演模拟，模拟OBS多分量数据接收，海面以下10m放炮，共正演400炮，炮间距为20m，第一炮在20m的位置，每炮有100个检波器，于海底(深度100m)接收，道间距为80m，如图5所示。在应用本发明提出的针对海底地震仪多分量数据的炮检互易弹性波全波形反演方法(Reciprocal-FWI)的同时，基于传统的伴随状态法(AD-FWI)和基于改进的散射积分算法波形反演(SI-FWI)也同时被应用，以对比突出本发明的有效性和优越性。
具体实施方式如下：
(1)数据采集器1采集地震波信号，并对原始地震数据进行重定位、多分量数据的旋转以及去噪等预处理(如图6所示)；
(2)输入设备4根据地下先验信息建立初始速度模型(如图4所示)，并设定反演频率范围与间隔等反演参数；
(3)处理器2在当前模型上，采用等效的炮检互易观测系统(如图7)分别进行x方向和z方向单向位移源的弹性波正演，通过前述线性叠加方式获得理论合成数据；
(4)处理器2判断理论合成数据与观测数据的匹配程度，如果“吻合”则保留当前模型并退出，并通过显示器显示模型(如图8、9、10、11所示)，否则继续执行以下步骤；
(5)处理器2计算并存储所有频率的检波点端波场和炮点端线性组合的反传波场，计算梯度，并累加形成整个观测系统所对应的梯度；
(6)处理器2计算最速下降方向；
(7)处理器2计算方向更新步长；
(8)处理器2计算模型更新量，并更新模型；
(9)重复步骤(4)-(8)；
利用Reciprocal-FWI、AD-FWI，SI-FWI进行反演得到的结果分别如图8、图9和图10所示。三种方法的反演结果均较好且基本一致，水平方向6km处垂直速度剖面对比如图11表明，反演结果仅仅有空间差分精度误差带来的细微区别。表1展示了三种算法每轮循环的计算时间，实验环境为英特尔xeon Silver 4110处理器，显然Reciprocal-FWI的计算效率要高于AD-FWI和SI-FWI，且不损失精度，体现了本发明的优势。
表1：三种方法每轮循环所需计算时间对比
| 方法 | Reciprocal-FWI | AD-FWI | SI-FWI |
| 计算耗时 | 6177seconds | 12327seconds | 9252seconds |
实施例2：
本实施例将本发明提出的针对海底地震仪多分量数据的炮检互易弹性波全波形反演方法应用到我国东海秋月工区获取的OBS实际资料中。模型共有1776×800个网格，网格间距18.75m×9m。共1669炮，炮间距18.75m。OBS检波器共计84台，间距100m。具体实施流程与实施例1相似，其中步骤(1)中还需要作三维到二维观测系统的数据转换，步骤(2)中的初始模型使用声波波形反演获取纵波速度模型，然后依据岩石物理经验关系转换横波速度和密度模型，主要对比计算效率，同样使用英特尔xeon Silver 4110，如表2所示，相较于传统的AD-FWI，加速比达到了10倍，相较于SI-FWI，加速比达到了5.5倍，也进一步验证了本发明相比于传统弹性波波形反演方法计算效率更高。
表2：三种方法每轮循环所需计算时间对比
Claims (8)
1.一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，其特征在于，包括以下步骤：
1)对原始地震数据进行预处理；
2)通过初至层析获取初始纵波速度模型，并通过声波波形反演细化初始纵波速度模型；
3)根据初至层析结果或声波波形反演结果，利用岩石物理关系建立横波以及密度的初始模型，并设定反演参数，反演参数包括反演频率范围和间隔；
4)根据波形和初始弹性参数，对初始模型进行炮检互易弹性波全波形反演，获取地下纵横波速度模型并成像，具体包括以下步骤：
41)根据所需要的空间精度，建立与原始观测系统等效的炮检互易观测系统；
42)在等效的炮检互易观测系统下，分别进行以x方向单向位移震源的时间域正演以及以z方向单向位移震源的时间域正演；
43)对正演得到的波场信息分量记录进行线性叠加，判断组合得到的理论合成数据与观测数据是否匹配，并采用二范数的目标函数进行检验，若目标函数值满足要求，则将当前纵横波速度模型作为最终纵横波速度模型，若否，则执行步骤44)；
44)在检波点处分别进行x方向和z方向的单向位移震源激发，采用时间域交错网格有限差分正演，分别获取所有检波点端波场，并将x分量残差和z分量残差作为残差震源，分别在炮点端以差分格式组合反传，获取炮点端的线性组合波场；
45)通过x方向的单向位移震源产生的检波点端波场与x分量残差作为残差震源产生的炮点端的线性组合波场的互相关求取对应梯度，累加形成整个观测系统所对应的x分量梯度，并通过z方向的单向位移震源产生的检波点端波场与z分量残差作为残差震源产生的炮点端的线性组合波场的互相关求取对应梯度，累加形成整个观测系统所对应的z分量梯度；
46)计算最速下降方向、方向更新步长和模型更新量后更新当前纵横波速度模型；
47)重复步骤41-46)，直到获得匹配的最终纵横波速度模型。
2.根据权利要求1所述的一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，其特征在于，所述的步骤1)中，预处理包括重定位、多分量数据的旋转以及去噪。
3.根据权利要求1所述的一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，其特征在于，所述的步骤41)中，等效的炮检互易观测系统具体为：
将震源的激发从在原始炮点位置激发改为在原始检波点位置激发； 将记录的接收位置从在原始检波点位置处接收改为在原始炮点处位置接收，且在炮点及炮点邻近点处组合接收。
4.根据权利要求1所述的一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，其特征在于，所述的步骤42)中，以x方向单向位移震源的时间域正演以及以z方向单向位移震源的时间域正演具体为：
针对接收的多分量数据为x和z分量记录的二维数据的情况，在等效的炮检互易观测系统下将x方向与z方向单向位移震源的正演分开单独进行，并记录各自所产生的波场信息分量记录uxx、uzx、uxz和uzz，其中，uij表示以j方向单向位移震源激发，i方向所接收到的数据分量，且j＝x或z，i＝x或z。
5.根据权利要求4所述的一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，其特征在于，所述的步骤43)中，将多空间点组合接收的分量记录按照空间差分的格式进行线性叠加，即可合成与观测数据所匹配的理论合成数据，具体包括以下步骤：
431)将uxx和uzx根据以下公式进行叠加，对应于实际观测数据的x分量所匹配的数据ux(xr,zr；xs,zs)：
其中，c(k)为交错网格有限差分系数，k为空间网格上移动的计数，(xr,zr)为检波点位置；
432)将uxz和uzz根据以下公式进行叠加，对应于实际观测数据的z分量所匹配的数据uz(xr,zr；xs,zs)
433)将叠加得到的合成数据与实际数据对应分量按照二范数进行匹配，则有：
其中，E(m)表示在地下模型参数m下所计算出的目标函数值，uij(xs,zs；xr,zr)表示在(xr,zr)位置处以j方向单向位移震源激发，在(xs,zs)位置处所接收到的i方向数据分量，di(xr,zr；xs,zs)表示在(xs,zs)位置处激发应力P波源，在(xr,zr)位置处所接收到的i方向数据分量，m为地下模型参数，ns为炮点总数，nr为检波点总数。
6.根据权利要求5所述的一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，其特征在于，所述的步骤44)中，检波点处激发产生的波场以及炮点端线性组合的反传波场的计算步骤为：
441)在检波点端激发x方向单向位移震源，存储波场信息分量记录uxx，uzx，将x分量观测记录的残差作为震源在(xs-1,zs)，(xs-2,zs)，…，(xs-M/2-1,zs)以及(xs,zs)，(xs+1,zs)，…，(xs+M/2,zs)处分别施加在x方向位移上，同时将x分量观测记录的残差作为震源在(xs,zs-1)，(xs,zs-2)，…，(xs,zs-M/2-1)以及(xs,zs)，(xs,zs+1)，…，(xs,zs+M/2)处施加在z方向位移上，计算波场并存储；
442)在检波点端激发z方向单向位移震源，存储波场信息分量记录uxz，uzz，将z分量观测记录的残差作为震源在(xs-1,zs)，(xs-2,zs)，…，(xs-M/2-1,zs)以及(xs,zs)，(xs+1,zs)，…，(xs+M/2,zs)处分别施加在x方向位移上，同时将z分量观测记录的残差作为震源在(xs,zs-1)，(xs,zs-2)，…，(xs,zs-M/2-1)以及(xs,zs)，(xs,zs+1)，…，(xs,zs+M/2)处施加在z方向位移上，计算波场并存储。
7.根据权利要求6所述的一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，其特征在于，所述的步骤45)中，将时间域计算得到的波场通过傅里叶变换得到频率域波场。
8.根据权利要求7所述的一种针对OBS多分量数据的炮检互易弹性波全波形反演方法，其特征在于，根据以下公式计算梯度：
其中，Re为取实部操作，*为复共轭算子，T表示转置，A为频率域负阻抗矩阵，u*i表示检波点处i方向激发的单分量位移震源所产生的波场在炮点处x和z方向接收的位移分量记录，ri表示检波点处观测记录i方向的残差。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010499751.8A CN111781639B (zh) | 2020-06-04 | 2020-06-04 | 针对obs多分量数据的炮检互易弹性波全波形反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010499751.8A CN111781639B (zh) | 2020-06-04 | 2020-06-04 | 针对obs多分量数据的炮检互易弹性波全波形反演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111781639A CN111781639A (zh) | 2020-10-16 |
| CN111781639B true CN111781639B (zh) | 2021-06-04 |
Family
ID=72753949
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010499751.8A Active CN111781639B (zh) | 2020-06-04 | 2020-06-04 | 针对obs多分量数据的炮检互易弹性波全波形反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111781639B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113447981B (zh) * | 2021-06-18 | 2022-07-19 | 同济大学 | 一种基于共成像点道集的反射全波形反演方法 |
| CN113777654B (zh) * | 2021-08-06 | 2023-07-04 | 同济大学 | 一种基于伴随状态法初至波走时层析的海水速度建模方法 |
| CN117420596A (zh) * | 2022-07-19 | 2024-01-19 | 中国石油天然气集团有限公司 | 地下介质速度的确定方法、装置、设备 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2479200A (en) * | 2010-04-01 | 2011-10-05 | Statoil Asa | Interpolating pressure and/or vertical particle velocity data from multi-component marine seismic data including horizontal derivatives |
| MY158191A (en) * | 2010-09-27 | 2016-09-15 | Exxonmobil Upstream Res Co | Hybrid method for full waveform inversion using simultaneous and sequential source method |
| US9702998B2 (en) * | 2013-07-08 | 2017-07-11 | Exxonmobil Upstream Research Company | Full-wavefield inversion of primaries and multiples in marine environment |
| US10942286B2 (en) * | 2017-08-24 | 2021-03-09 | Chevron U.S.A. Inc. | System and method for image-domain full waveform inversion |
| CN109507726A (zh) * | 2017-09-15 | 2019-03-22 | 中国石油化工股份有限公司 | 时间域弹性波多参数全波形的反演方法及系统 |
| CN109541681B (zh) * | 2017-09-22 | 2020-07-17 | 中国海洋大学 | 一种拖缆地震数据和少量obs数据联合的波形反演方法 |
| US11353613B2 (en) * | 2017-11-17 | 2022-06-07 | Cgg Services Sas | Seismic exploration using image-based reflection full waveform inversion to update low wavenumber velocity model |
| CN108345031B (zh) * | 2018-01-11 | 2020-01-17 | 吉林大学 | 弹性介质主动源和被动源混采地震数据全波形反演方法 |
| CN109387868A (zh) * | 2018-09-28 | 2019-02-26 | 中国海洋石油集团有限公司 | 一种基于地震波同相轴斜率信息的三维层析成像方法 |
| CN110764146B (zh) * | 2019-10-24 | 2020-05-19 | 南京信息工程大学 | 一种基于声波算子的空间互相关弹性波反射波形反演方法 |
| CN111158049B (zh) * | 2019-12-27 | 2020-11-27 | 同济大学 | 一种基于散射积分法的地震逆时偏移成像方法 |
-
2020
- 2020-06-04 CN CN202010499751.8A patent/CN111781639B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111781639A (zh) | 2020-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111781639B (zh) | 针对obs多分量数据的炮检互易弹性波全波形反演方法 | |
| CN108345031B (zh) | 弹性介质主动源和被动源混采地震数据全波形反演方法 | |
| KR101797451B1 (ko) | 상호상관 목적 함수를 통한 해양 스트리머 데이터에 대한 동시 소스 반전 | |
| US8296069B2 (en) | Pseudo-analytical method for the solution of wave equations | |
| US20100054082A1 (en) | Reverse-time depth migration with reduced memory requirements | |
| CN111158049B (zh) | 一种基于散射积分法的地震逆时偏移成像方法 | |
| BRPI0918020B1 (pt) | métodos de exploração sísmica | |
| WO2012139082A1 (en) | Event selection in the image domain | |
| CN109738952B (zh) | 基于全波形反演驱动的被动源直接偏移成像方法 | |
| CN109541681B (zh) | 一种拖缆地震数据和少量obs数据联合的波形反演方法 | |
| US20110276273A1 (en) | Image domain signal to noise estimate | |
| CN111045077B (zh) | 一种陆地地震数据的全波形反演方法 | |
| Lin et al. | Guided wave tomography based on supervised descent method for quantitative corrosion imaging | |
| EP3014304A2 (en) | System and method for estimating repeatability using base data | |
| CA3024240A1 (en) | System and method for real-time passive seismic event localization | |
| CN105093278A (zh) | 基于激发主能量优化算法的全波形反演梯度算子的提取方法 | |
| US9658354B2 (en) | Seismic imaging systems and methods employing correlation-based stacking | |
| CN116068493A (zh) | 一种深海大深度垂直分布式水听器被动声源定位方法 | |
| CN113740915A (zh) | 一种球坐标系重力和接收函数联合反演地壳结构参数的方法 | |
| CN112285775B (zh) | 一种海底浅部沉积层剪切波速度探测装置及方法 | |
| CN109188527B (zh) | 滩浅海地区快速建立三维近海底速度模型的方法 | |
| CN111665563A (zh) | 基于聚焦分析的叠前偏移垂向分辨率评价方法 | |
| US9465122B2 (en) | Method for acquiring seismic data | |
| CN110873895A (zh) | 一种变网格微地震逆时干涉定位方法 | |
| CN111239828B (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 | ||
| GR01 | Patent grant |