RU2008151147A - Способ вычисления точного импульсного отклика плоского акустического отражателя для точечного акустического источника при нулевом смещении - Google Patents
Способ вычисления точного импульсного отклика плоского акустического отражателя для точечного акустического источника при нулевом смещении Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract 15
- 238000006073 displacement reaction Methods 0.000 title 1
- 238000004422 calculation algorithm Methods 0.000 claims 4
- 230000005012 migration Effects 0.000 claims 1
- 238000013508 migration Methods 0.000 claims 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/30—Analysis
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) источника, содержащий этапы, на которых:
а. обеспечивают уравнение:
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) с помощью аргументов на основе асимптотической теории распространения волн, а также принципа причинности, и далее выводят уравнение
m. окончательно получают из уравнения, заданного на этапе (l), следующее уравнение
с помощью принципа свертки, где
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) плоского акустического отражателя при нулевом смещении для точечного акустического источника согласно
где 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,
где 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), с помощью аргументов на основе асимптотической теории распространения волн, а также принципа причинности, и далее выводят уравнение
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)
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)
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 |
-
2007
- 2007-02-21 MX MX2009001681A patent/MX2009001681A/es active IP Right Grant
- 2007-02-21 EP EP07736532.8A patent/EP2089738B1/en not_active Expired - Fee Related
- 2007-02-21 RU RU2008151147/28A patent/RU2008151147A/ru unknown
- 2007-02-21 AU AU2007330350A patent/AU2007330350A1/en not_active Abandoned
- 2007-02-21 CA CA2659782A patent/CA2659782C/en active Active
- 2007-02-21 WO PCT/IN2007/000071 patent/WO2008068769A1/en active Application Filing
- 2007-02-26 US US11/711,469 patent/US7940601B2/en not_active Expired - Fee Related
-
2009
- 2009-07-07 NO NO20092590A patent/NO20092590L/no not_active Application Discontinuation
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 |