RU2008151147A - Способ вычисления точного импульсного отклика плоского акустического отражателя для точечного акустического источника при нулевом смещении - Google Patents

Способ вычисления точного импульсного отклика плоского акустического отражателя для точечного акустического источника при нулевом смещении Download PDF

Info

Publication number
RU2008151147A
RU2008151147A RU2008151147/28A RU2008151147A RU2008151147A RU 2008151147 A RU2008151147 A RU 2008151147A RU 2008151147/28 A RU2008151147/28 A RU 2008151147/28A RU 2008151147 A RU2008151147 A RU 2008151147A RU 2008151147 A RU2008151147 A RU 2008151147A
Authority
RU
Russia
Prior art keywords
equation
reflector
source
reflected signal
wave velocity
Prior art date
Application number
RU2008151147/28A
Other languages
English (en)
Inventor
Шанти Кумар ГОШ (IN)
Шанти Кумар ГОШ
Original Assignee
Каусел Оф Сайнтифик Энд Индастриал Рисерч (In)
Каусел Оф Сайнтифик Энд Индастриал Рисерч
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 Каусел Оф Сайнтифик Энд Индастриал Рисерч (In), Каусел Оф Сайнтифик Энд Индастриал Рисерч filed Critical Каусел Оф Сайнтифик Энд Индастриал Рисерч (In)
Publication of RU2008151147A publication Critical patent/RU2008151147A/ru

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Aerials With Secondary Devices (AREA)

Abstract

1. Способ вычисления импульсного отклика [uS(t)] и отраженного сигнала [r(t)] плоского акустического отражателя при нулевом смещении для точечного акустического источника, описанного временной функцией s(t) источника, содержащий этапы, на которых: ! а. обеспечивают уравнение: ! ! где ! uS(ω) - преобразование Фурье от uS(t), импульсный отклик; ! ω - угловая частота; ! i - мнимая единица; ! x - координата; ! с - эталонная волновая скорость, около которой имеет место возмущение скорости; ! b. заменяют α(x3), как задано на этапе (а), при этом рассматривают бесконечную плоскую границу раздела на x3=h, на которой волновая скорость изменяется от с до c+Δc, которая составляет отражатель, эквивалент специальному случаю, где α(x3)=α0H(x3-h), где H обозначает функцию Хевисайда; ! с. изменяют переменную согласно t=2x3/с, что вместе с этапом (b) сводит уравнение этапа (а) к ! , ! где !d. выводят из уравнения на этапе (с), что , причем «↔» обозначает, что члены слева и справа от этого символа составляют пару Фурье; ! е. используют уравнения на этапах (с) и (d) для получения уравнения ! , ! f. используют уравнение на этапе (е) для дальнейшего вывода выражения: ! , ! g. идентифицируют левую часть выражения на этапе (f) как преобразование Фурье от (t-τ)H(t-τ) и выводят уравнение ! , ! h. выводят далее из уравнения, как оно задано на этапе (g), следующее выражение !, ! i. преобразуют уравнение из этапа (h) в следующее уравнение ! , ! j. модифицируют уравнение из этапа (i) в следующее уравнение ! , ! k. упрощают уравнение, заданное на этапе (j), до ! , ! где uS(t) представляет собой рассеянный отклик импульса, ! l. устраняют ограничение аппроксимации Берна, присущее уравнению на этапе (k) с помощью аргумент�

Claims (7)

1. Способ вычисления импульсного отклика [uS(t)] и отраженного сигнала [r(t)] плоского акустического отражателя при нулевом смещении для точечного акустического источника, описанного временной функцией s(t) источника, содержащий этапы, на которых:
а. обеспечивают уравнение:
Figure 00000001
где
Figure 00000002
uS(ω) - преобразование Фурье от uS(t), импульсный отклик;
ω - угловая частота;
i - мнимая единица;
x - координата;
с - эталонная волновая скорость, около которой имеет место возмущение скорости;
b. заменяют α(x3), как задано на этапе (а), при этом рассматривают бесконечную плоскую границу раздела на x3=h, на которой волновая скорость изменяется от с до c+Δc, которая составляет отражатель, эквивалент специальному случаю, где α(x3)=α0H(x3-h), где H обозначает функцию Хевисайда;
с. изменяют переменную согласно t=2x3/с, что вместе с этапом (b) сводит уравнение этапа (а) к
Figure 00000003
,
где
Figure 00000004
d. выводят из уравнения на этапе (с), что
Figure 00000005
, причем «↔» обозначает, что члены слева и справа от этого символа составляют пару Фурье;
е. используют уравнения на этапах (с) и (d) для получения уравнения
Figure 00000006
,
f. используют уравнение на этапе (е) для дальнейшего вывода выражения:
Figure 00000007
,
g. идентифицируют левую часть выражения на этапе (f) как преобразование Фурье от (t-τ)H(t-τ) и выводят уравнение
Figure 00000008
,
h. выводят далее из уравнения, как оно задано на этапе (g), следующее выражение
Figure 00000009
,
i. преобразуют уравнение из этапа (h) в следующее уравнение
Figure 00000010
,
j. модифицируют уравнение из этапа (i) в следующее уравнение
Figure 00000011
,
k. упрощают уравнение, заданное на этапе (j), до
Figure 00000012
,
где uS(t) представляет собой рассеянный отклик импульса,
l. устраняют ограничение аппроксимации Берна, присущее уравнению на этапе (k) с помощью аргументов на основе асимптотической теории распространения волн, а также принципа причинности, и далее выводят уравнение
Figure 00000013
,
m. окончательно получают из уравнения, заданного на этапе (l), следующее уравнение
Figure 00000014
с помощью принципа свертки, где
t - время наблюдения, которое начинается при подрыве источника,
s(t) - временная функция источника,
R=Δc/(2c+Δc),
с - эталонная волновая скорость,
Δc - изменение волновой скорости,
h - кратчайшее расстояние между источником и отражателем,
τ = 2h/c, a * означает свертку,
H - функция Хевисайда, обозначаемая через H(t-τ) и определяемая как
H(t-τ)=1, t>τ
=0, t<τ,
n. вычисляют импульсный отклик плоского акустического отражателя при нулевом смещении для точечного акустического источника, как заявлено на этапе m, при этом упомянутый способ содержит этапы, на которых:
i) вводят расстояние (h) от отражателя, эталонную волновую скорость (с), изменение (Δc) волновой скорости, временную функцию s(t) источника и время (Т) усечения импульсного отклика отражателя в уравнение, заявленное на этапе (l) по п.1,
ii) получают импульсный отклик uS(t) плоского акустического отражателя при нулевом смещении для точечного акустического источника согласно
Figure 00000015
,
где R=Δc/(2c+Δc),
t - время наблюдения, которое начинается при подрыве источника,
δ(t) - дельта-функция Дирака,
с - эталонная волновая скорость,
Δc - изменение волновой скорости,
h - кратчайшее расстояние между источником и отражателем,
τ = 2h/c, a * означает свертку,
H - функция Хевисайда, обозначаемая через H(t-τ) и определяемая как
H(t-τ)=1, t>τ
= 0, t<τ,
iii) получают отраженный сигнал [r(t)] плоского акустического отражателя при нулевом смещении для точечного акустического источника, описанного временной функцией s(t) источника, из выражения, полученного на этапе (m) по п.1,
Figure 00000016
,
где R=Δc/(2c+Δc),
τ = 2h/c, а * означает свертку,
Н - функция Хевисайда,
t - время наблюдения, которое начинается при подрыве источника,
s(t) - временная функция источника,
с - эталонная волновая скорость,
Δc - изменение волновой скорости,
h - кратчайшее расстояние между источником и отражателем,
о. проверяют достоверность алгоритма сейсмического численного моделирования (известный способ), при этом упомянутый способ содержит этапы, на которых:
i) получают отраженный сигнал при нулевом смещении от плоского акустического отражателя, вычисленный по алгоритму сейсмического численного моделирования;
ii) получают отраженный сигнал при нулевом смещении от плоского акустического отражателя, вычисленный по способу, заявленному на этапе (n); и
iii) сравнивают данные, полученные на этапе (i), с данными, полученными на этапе (ii), для проверки достоверности алгоритма сейсмического численного моделирования.
2. Способ вычисления импульсного отклика [uS(t)] и отраженного сигнала [r(t)] плоского акустического отражателя по п.1, в котором отраженный сигнал при нулевом смещении кругового отражателя вычисляют с учетом наблюдений на его центральной оси и в течение ограниченного временного интервала.
3. Способ вычисления импульсного отклика [uS(t)] и отраженного сигнала [r(t)] плоского акустического отражателя по п.1, в котором точное сейсмическое сечение при нулевом смещении вычисляют для наклонного или горизонтального плоскостного отражателя для подтверждения достоверности алгоритмов сейсмической миграции и инверсии.
4. Способ вычисления импульсного отклика [uS(t)] и отраженного сигнала [r(t)] плоского акустического отражателя по п.1, характеризующийся тем, что упомянутый способ проверяет достоверность интерпретации отражателя как плоскостной структуры.
5. Способ вычисления импульсного отклика [uS(t)] и отраженного сигнала [r(t)] плоского акустического отражателя по п.1, характеризующийся тем, что упомянутый способ вычисляет временную функцию s(t) сейсмического источника путем введения в качестве входных данных отраженного сигнала плоского отражателя при нулевом смещении.
6. Способ вычисления импульсного отклика [uS(t)] и отраженного сигнала [r(t)] плоского акустического отражателя по п.1, характеризующийся этапом (b), на котором подставляют α(x3) на этапе (а) и рассматривают бесконечную плоскую границу раздела при x3=h, на которой волновая скорость изменяется от с до c+Δc и которая составляет отражатель, эквивалентный для специального случая, когда α(x3)=α0H(x3-h).
7. Способ вычисления импульсного отклика [uS(t)] и отраженного сигнала [r(t)] плоского акустического отражателя по п.1, характеризующийся этапом (l), на котором удаляют ограничение аппроксимации Берна, присущее уравнению на этапе (k), с помощью аргументов на основе асимптотической теории распространения волн, а также принципа причинности, и далее выводят уравнение
Figure 00000017
.
RU2008151147/28A 2006-12-07 2007-02-21 Способ вычисления точного импульсного отклика плоского акустического отражателя для точечного акустического источника при нулевом смещении RU2008151147A (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
IN2626DE2006 2006-12-07
IN2626/DEL/2006 2006-12-07

Publications (1)

Publication Number Publication Date
RU2008151147A true RU2008151147A (ru) 2010-06-27

Family

ID=39226987

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2008151147/28A RU2008151147A (ru) 2006-12-07 2007-02-21 Способ вычисления точного импульсного отклика плоского акустического отражателя для точечного акустического источника при нулевом смещении

Country Status (8)

Country Link
US (1) US7940601B2 (ru)
EP (1) EP2089738B1 (ru)
AU (1) AU2007330350A1 (ru)
CA (1) CA2659782C (ru)
MX (1) MX2009001681A (ru)
NO (1) NO20092590L (ru)
RU (1) RU2008151147A (ru)
WO (1) WO2008068769A1 (ru)

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8694299B2 (en) 2010-05-07 2014-04-08 Exxonmobil Upstream Research Company Artifact reduction in iterative inversion of geophysical data
US8892413B2 (en) 2011-03-30 2014-11-18 Exxonmobil Upstream Research Company Convergence rate of full wavefield inversion using spectral shaping
WO2012134609A1 (en) 2011-03-31 2012-10-04 Exxonmobil Upstream Research Company Method of wavelet estimation and multiple prediction in full wavefield inversion
US9140812B2 (en) 2011-09-02 2015-09-22 Exxonmobil Upstream Research Company Using projection onto convex sets to constrain full-wavefield inversion
US9176930B2 (en) 2011-11-29 2015-11-03 Exxonmobil Upstream Research Company Methods for approximating hessian times vector operation in full wavefield inversion
EP2823335A4 (en) 2012-03-08 2016-01-13 Exxonmobil Upstream Res Co ORTHOGONAL SOURCE AND RECEIVER CODING
US10317548B2 (en) 2012-11-28 2019-06-11 Exxonmobil Upstream Research Company Reflection seismic data Q tomography
SG11201508195PA (en) 2013-05-24 2015-12-30 Exxonmobil Upstream Res Co Multi-parameter inversion through offset dependent elastic fwi
US10459117B2 (en) 2013-06-03 2019-10-29 Exxonmobil Upstream Research Company Extended subspace method for cross-talk mitigation in multi-parameter inversion
US9702998B2 (en) 2013-07-08 2017-07-11 Exxonmobil Upstream Research Company Full-wavefield inversion of primaries and multiples in marine environment
WO2015026451A2 (en) 2013-08-23 2015-02-26 Exxonmobil Upstream Research Company Simultaneous sourcing during both seismic acquisition and seismic inversion
US10036818B2 (en) 2013-09-06 2018-07-31 Exxonmobil Upstream Research Company Accelerating full wavefield inversion with nonstationary point-spread functions
US9910189B2 (en) 2014-04-09 2018-03-06 Exxonmobil Upstream Research Company Method for fast line search in frequency domain FWI
SG11201608175SA (en) 2014-05-09 2016-11-29 Exxonmobil Upstream Res Co Efficient line search methods for multi-parameter full wavefield inversion
US10185046B2 (en) 2014-06-09 2019-01-22 Exxonmobil Upstream Research Company Method for temporal dispersion correction for seismic simulation, RTM and FWI
AU2015280633B2 (en) 2014-06-17 2018-07-19 Exxonmobil Upstream Research Company Fast viscoacoustic and viscoelastic full-wavefield inversion
US10838092B2 (en) 2014-07-24 2020-11-17 Exxonmobil Upstream Research Company Estimating multiple subsurface parameters by cascaded inversion of wavefield components
US10422899B2 (en) 2014-07-30 2019-09-24 Exxonmobil Upstream Research Company Harmonic encoding for FWI
US10386511B2 (en) 2014-10-03 2019-08-20 Exxonmobil Upstream Research Company Seismic survey design using full wavefield inversion
EP3210050A1 (en) 2014-10-20 2017-08-30 Exxonmobil Upstream Research Company Velocity tomography using property scans
US11163092B2 (en) 2014-12-18 2021-11-02 Exxonmobil Upstream Research Company Scalable scheduling of parallel iterative seismic jobs
US10520618B2 (en) 2015-02-04 2019-12-31 ExxohnMobil Upstream Research Company Poynting vector minimal reflection boundary conditions
US10317546B2 (en) 2015-02-13 2019-06-11 Exxonmobil Upstream Research Company Efficient and stable absorbing boundary condition in finite-difference calculations
MX2017007988A (es) 2015-02-17 2017-09-29 Exxonmobil Upstream Res Co Proceso de inversion de campo ondulatorio completo de multifase que genera un conjunto de datos libres de multiples.
CA2985738A1 (en) 2015-06-04 2016-12-08 Exxonmobil Upstream Research Company Method for generating multiple free seismic images
US10838093B2 (en) 2015-07-02 2020-11-17 Exxonmobil Upstream Research Company Krylov-space-based quasi-newton preconditioner for full-wavefield inversion
CN108139499B (zh) 2015-10-02 2020-02-14 埃克森美孚上游研究公司 Q-补偿的全波场反演
MX2018003495A (es) 2015-10-15 2018-06-06 Exxonmobil Upstream Res Co Apilados angulares de dominio de modelo de fwi con conservacion de amplitud.
US10768324B2 (en) 2016-05-19 2020-09-08 Exxonmobil Upstream Research Company Method to predict pore pressure and seal integrity using full wavefield inversion
RU2760508C1 (ru) * 2020-10-01 2021-11-25 Общество с ограниченной ответственностью «Научно-производственный центр неразрушающего контроля «ЭХО+» Способ определения типа отражателя по амплитуде рассеянных им ультразвуковых импульсов

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2293010B (en) * 1994-07-07 1998-12-09 Geco As Method of processing seismic data
NO320271B1 (no) * 2000-05-26 2005-11-14 Norsar V Daglig Leder Fremgangsmate og anordning for a finne refleksjonsstyrken til reflektorer i undergrunnen mens refleksjonskoeffisienten for en malhorisont holdes konstant

Also Published As

Publication number Publication date
AU2007330350A1 (en) 2008-06-12
US7940601B2 (en) 2011-05-10
US20080137477A1 (en) 2008-06-12
MX2009001681A (es) 2009-02-26
CA2659782C (en) 2015-04-21
NO20092590L (no) 2009-09-04
CA2659782A1 (en) 2008-06-12
EP2089738B1 (en) 2019-01-02
WO2008068769A1 (en) 2008-06-12
EP2089738A1 (en) 2009-08-19

Similar Documents

Publication Publication Date Title
RU2008151147A (ru) Способ вычисления точного импульсного отклика плоского акустического отражателя для точечного акустического источника при нулевом смещении
AU2018208683B2 (en) Flow meter configuration and calibration
EP2138862B1 (en) Object detection by evaluating a forward scattered sound wave (FSSW) by comparing the passive-FSSW-phase conjugate with the autocorrelation of the direct, unscattered sound wave
CN105102924A (zh) 对壁表面的变化的超声波检测
CN104062645B (zh) 一种测量参量阵差频波与同频小振幅波信号相位差的方法
CN105115552B (zh) 基于侧音测相的超声波流量测量方法及装置
RU2011130315A (ru) Гидроакустический комплекс для измерения азимутального угла на источник звука в мелком море
JP2007040734A (ja) 送信モードリコメンドシステム、ソーナーシステム、送信モードリコメンド方法、そのプログラム及び記録媒体
Aunsri et al. Stochastic description and evaluation of ocean acoustics time-series for frequency and dispersion estimation using particle filtering approach
Zhang et al. A passive acoustic positioning algorithm based on virtual long baseline matrix window
WO2019069804A1 (ja) 伝搬時間測定器、気体濃度測定装置、伝搬時間測定プログラム、および伝搬時間測定方法
CN115236592B (zh) 一种基于单阵元时频曲线匹配的冰声定位方法
CN113126029B (zh) 适用于深海可靠声路径环境的多传感器脉冲声源定位方法
Prior et al. Modeling acoustical pressure and particle acceleration close to marine seismic airguns and airgun arrays
CN109188378B (zh) 一种超声多普勒回波信号仿真方法
JP4249332B2 (ja) 周波数測定方法およびドップラソナー
JP2021038945A (ja) 距離推定装置、距離推定方法、および距離推定プログラム
CN104678380A (zh) 一种高速小目标检测中基于lms算法的直达波抑制方法
Bonito Acoustic system for ground truth underwater positioning in DEEC's test tank
RU2308054C2 (ru) Гидроакустическая синхронная дальномерная навигационная система
CN113191307B (zh) 一种基于贝叶斯估计的浅海波导简正波阶数判定方法
CN111337881B (zh) 一种利用螺旋桨噪声的水下目标探测方法
CN112114364B (zh) 偶极横波反射波补偿方法及装置
Kim et al. Identification of the impact location in a gas duct system based on acoustic wave theory and the time frequency
Yardim et al. Particle smoothers in sequential geoacoustic inversion