RU2016115039A - Решающее устройство эйконала для квази продольных волн в анизотропной среде - Google Patents
Решающее устройство эйконала для квази продольных волн в анизотропной среде Download PDFInfo
- Publication number
- RU2016115039A RU2016115039A RU2016115039A RU2016115039A RU2016115039A RU 2016115039 A RU2016115039 A RU 2016115039A RU 2016115039 A RU2016115039 A RU 2016115039A RU 2016115039 A RU2016115039 A RU 2016115039A RU 2016115039 A RU2016115039 A RU 2016115039A
- Authority
- RU
- Russia
- Prior art keywords
- equation
- value
- seismic wave
- propagation time
- intermediate term
- Prior art date
Links
- 238000000034 method Methods 0.000 claims 14
- 238000011156 evaluation Methods 0.000 claims 6
- 238000000926 separation method Methods 0.000 claims 4
- QNRATNLHPGXHMA-XZHTYLCXSA-N (r)-(6-ethoxyquinolin-4-yl)-[(2s,4s,5r)-5-ethyl-1-azabicyclo[2.2.2]octan-2-yl]methanol;hydrochloride Chemical compound Cl.C([C@H]([C@H](C1)CC)C2)CN1[C@@H]2[C@H](O)C1=CC=NC2=CC=C(OCC)C=C21 QNRATNLHPGXHMA-XZHTYLCXSA-N 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
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/58—Media-related
- G01V2210/586—Anisotropic media
-
- 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/67—Wave propagation modeling
- G01V2210/675—Wave equation; Green's functions
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Acoustics & Sound (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Claims (41)
1. Способ определения времени распространения сейсмической волны в модели, включающий:
получение модели подземной области, содержащей анизотропную среду, при этом модель содержит сетку, имеющую узловые точки, представляющие местоположения в подземной области и местоположение источника;
определение уравнения эйконала для вычисления времени распространения сейсмической волны от местоположения источника через анизотропную среду по меньшей мере до одной из узловых точек;
разделение уравнения эйконала на первое уравнение и второе уравнение; при этом итерационное решение первого уравнения с использованием процессора таким образом, что определzют первое значение времени распространения сейсмической волны; и
числовую оценку второго уравнения, основанную на первом значении времени распространения сейсмической волны.
2. Способ по п.1, отличающийся тем, что разделение уравнения эйконала на первое уравнение и второе уравнение включает приравнивание части первого уравнения к промежуточному члену и приравнивание части второго уравнения к этому промежуточному члену.
3. Способ по п.2, дополнительно включающий назначение начального значения промежуточному члену прежде итерационного решения первого уравнения.
4. Способ по п.3, отличающийся тем, что числовая оценка второго уравнения включает определение первого вычисленного значения для промежуточного члена.
5. Способ по п.4, дополнительно включающий итерационное решение первого уравнения на основе первого вычисленного значения промежуточного члена таким образом, что определяют второе значение времени распространения сейсмической волны, при этом второе значение времени распространения сейсмической волны отличается от первого значения времени распространения сейсмической волны.
6. Способ по п.5, дополнительно включающий числовую оценку второго уравнения на основе расчета второго значения времени распространения сейсмической волны таким образом, что определяют второе вычисленное значение промежуточного члена.
7. Способ по п.2, отличающийся тем, что итерационное решение первого уравнения и числовая оценка второго уравнения выполняются в первой итерации, при этом способ дополнительно включает:
выполнение одной или более дополнительных итераций, включающих:
итерационное решение первого уравнения в дополнительное время на основе вычисленного значения промежуточного члена, вычисленного в предыдущей итерации, таким образом, что определяется новое значение времени распространения сейсмической волны; и
числовую оценку второго уравнения в дополнительное время выполняют таким образом, что определяют новое вычисленное значение промежуточного члена на основе нового значения времени распространения сейсмической волны.
8. Способ по п.7, дополнительно включающий экстраполяцию значения сходимости для значения времени распространения сейсмической волны на основе расчета времени распространения сейсмической волны, определенного в первой итерации, и при этом новое значение времени распространения сейсмической волны, определено по меньшей мере в одной из одной или более дополнительных итераций или их комбинаций.
9. Способ по п.1, отличающийся тем, что итерационное решение первого уравнения включает:
выбор узловой точки сети;
определение числа направлений из указанной узловой точки, на которых одна или более соседних узловых точек имеют вычисленное значение времени распространения сейсмической волны; и
если число направлений составляет по меньшей мере одно, то вычисление значения времени распространения сейсмической волны для указанной узловой точки основанное на первом уравнении и значении времени распространения сейсмической волны для одной или более соседних узловых точек.
10. Способ по п.9, отличающийся тем, что итерационное решение первого уравнения дополнительно включает определение обусловленности значения времени распространения сейсмической волны для узловой точки, по меньшей мере когда число направлений составляет по меньшей мере два.
11. Способ по п.9, отличающийся тем, что итерационное решение первого уравнения включает назначение начального значения времени распространения сейсмической волны местоположению источника.
12. Способ по п.1, отличающийся тем, что уравнение эйконала представляет собой значения времени распространения сейсмической волны в среде, имеющей анизотропию, выбранную из группы, состоящей из: наклонной и ромбической, наклонной и поперечной, моноклинной и триклинной анизотропии.
13. Способ по п.1, дополнительно включающий:
обновление модели с использованием значения времени распространения сейсмической волны и
отображение местоположения волны в модели один или более раз после того, как источник излучает или отражает волну.
14. Вычислительное устройство, содержащее:
один или более процессоров и
запоминающее устройство, содержащее один или более энергонезависимых машиночитаемых носителей, содержащих команды, которые при выполнении по меньшей мере одним из одного или более процессорами вызывают выполнение вычислительным устройством операций, включающих:
получение модели подземной области, содержащей анизотропную среду, при этом модель содержит сетку, имеющую узловые точки, представляющие местоположения в подземной области и местоположение источника;
определение уравнения эйконала для вычисления времени распространения сейсмической волны от местоположения источника через анизотропную среду по меньшей мере до одной из узловых точек;
разделение уравнения эйконала на первое уравнение и второе уравнение; при этом итерационно решается первое уравнение таким образом, что определяется первое значение времени распространения сейсмической волны; и
числовую оценку второго уравнения, основанную на первом значении времени распространения сейсмической волны.
15. Устройство по п.14, отличающееся тем, что разделение уравнения эйконала на первое уравнение и второе уравнение включает приравнивание части первого уравнения к промежуточному члену и приравнивание части второго уравнения к этому промежуточному члену.
16. Устройство по п.15, дополнительно включающее назначение начального значения промежуточному члену прежде итерационного решения первого уравнения.
17. Устройство по п.16, отличающееся тем, что числовая оценка второго уравнения включает определение первого вычисленного значения для промежуточного члена.
18. Устройство по п.17, дополнительно включающее итерационное решение первого уравнения на основе первого вычисленного значения промежуточного члена таким образом, что определяется второе значение времени распространения сейсмической волны, при этом второе значение времени распространения сейсмической волны отличается от первого значения времени распространения сейсмической волны.
19. Устройство по п.18, отличающееся тем, что операции дополнительно включают числовую оценку второго уравнения на основе второго значения времени распространения сейсмической волны таким образом, что определяется второе вычисленное значение промежуточного члена.
20. Устройство по п.15, отличающееся тем, что итерационное решение первого уравнения и числовая оценка второго уравнения выполняются в первую итерацию, при этом операции дополнительно включают:
выполнение одной или более дополнительных итераций, включающих:
итерационное решение первого уравнения в дополнительное время на основе вычисленного значения промежуточного члена, вычисленного в предыдущей итерации таким образом, что определяется новое значение времени распространения сейсмической волны; и
числовую оценку второго уравнения в дополнительное время таким образом, что определяется новое вычисленное значение промежуточного члена на основе нового значения времени распространения сейсмической волны.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361880265P | 2013-09-20 | 2013-09-20 | |
US61/880,265 | 2013-09-20 | ||
PCT/US2014/056539 WO2015042386A1 (en) | 2013-09-20 | 2014-09-19 | Eikonal solver for quasi p-waves in anisotropic media |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2016115039A true RU2016115039A (ru) | 2017-10-25 |
Family
ID=52689438
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2016115039A RU2016115039A (ru) | 2013-09-20 | 2014-09-19 | Решающее устройство эйконала для квази продольных волн в анизотропной среде |
Country Status (5)
Country | Link |
---|---|
US (1) | US20160202375A1 (ru) |
EP (1) | EP3047310A4 (ru) |
MX (1) | MX2016003639A (ru) |
RU (1) | RU2016115039A (ru) |
WO (1) | WO2015042386A1 (ru) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675905B (zh) * | 2012-09-14 | 2016-10-05 | 中国科学院地质与地球物理研究所 | 一种基于优化系数的地震波场模拟方法及装置 |
US11199644B2 (en) | 2016-09-08 | 2021-12-14 | King Abdullah University Of Science And Technology | Imaging of geological subsurface features |
CN107894613B (zh) * | 2017-10-26 | 2019-07-26 | 中国石油天然气集团公司 | 弹性波矢量成像方法、装置、存储介质及设备 |
CN109581500B (zh) * | 2018-12-18 | 2020-06-30 | 东华理工大学 | 一种反射地震记录频变速度分析方法 |
CN111948706B (zh) * | 2019-05-16 | 2023-06-30 | 中国石油天然气集团有限公司 | 正交各向异性介质地震成像方法及装置 |
CN112764040B (zh) * | 2019-11-01 | 2022-06-14 | 复旦大学 | 一种基于射线理论相位修正的合成孔径波束形成方法 |
CN111859253B (zh) * | 2020-07-08 | 2023-09-22 | 上海雪湖科技有限公司 | 一种基于fpga的高阶波动方程求解方法 |
CN112505765B (zh) * | 2020-11-18 | 2023-05-09 | 东华理工大学 | Lax Friedrichs扫描地震波旅行时间的方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5394325A (en) * | 1993-04-07 | 1995-02-28 | Exxon Production Research Company | Robust, efficient three-dimensional finite-difference traveltime calculations |
US6035256A (en) * | 1997-08-22 | 2000-03-07 | Western Atlas International, Inc. | Method for extrapolating traveltimes across shadow zones |
FR2769719B1 (fr) * | 1997-10-10 | 2000-01-07 | Geophysique Cie Gle | Procede de traitement sismique et notamment procede de prospection sismique 3d mettant en oeuvre une migration des donnees sismiques |
US6018499A (en) * | 1997-11-04 | 2000-01-25 | 3Dgeo Development, Inc. | Three-dimensional seismic imaging of complex velocity structures |
US6324478B1 (en) * | 1999-05-10 | 2001-11-27 | 3D Geo Development, Inc. | Second-and higher-order traveltimes for seismic imaging |
US6643590B2 (en) * | 2002-01-04 | 2003-11-04 | Westerngeco, L.L.C. | Method for computing finite-frequency seismic migration traveltimes from monochromatic wavefields |
CA2606686C (en) * | 2005-05-26 | 2015-02-03 | Exxonmobil Upstream Research Company | A rapid method for reservoir connectivity analysis using a fast marching method |
US20070162249A1 (en) * | 2006-01-06 | 2007-07-12 | Min Lou | Traveltime calculation in three dimensional transversely isotropic (3D TTI) media by the fast marching method |
US7460437B2 (en) * | 2007-01-03 | 2008-12-02 | Weinman Geoscience | Seismic data processing method and system for migration of seismic signals incorporating azimuthal variations in the velocity |
US8274859B2 (en) * | 2010-02-22 | 2012-09-25 | Landmark Graphics Corporation | Systems and methods for modeling 3D geological structures |
US9874648B2 (en) * | 2011-02-21 | 2018-01-23 | Exxonmobil Upstream Research Company | Reservoir connectivity analysis in a 3D earth model |
-
2014
- 2014-09-19 US US15/023,628 patent/US20160202375A1/en not_active Abandoned
- 2014-09-19 RU RU2016115039A patent/RU2016115039A/ru unknown
- 2014-09-19 MX MX2016003639A patent/MX2016003639A/es unknown
- 2014-09-19 WO PCT/US2014/056539 patent/WO2015042386A1/en active Application Filing
- 2014-09-19 EP EP14845719.5A patent/EP3047310A4/en not_active Withdrawn
Also Published As
Publication number | Publication date |
---|---|
WO2015042386A1 (en) | 2015-03-26 |
EP3047310A1 (en) | 2016-07-27 |
EP3047310A4 (en) | 2017-03-08 |
MX2016003639A (es) | 2016-10-28 |
US20160202375A1 (en) | 2016-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2016115039A (ru) | Решающее устройство эйконала для квази продольных волн в анизотропной среде | |
Xia et al. | A full-scale fluvial flood modelling framework based on a high-performance integrated hydrodynamic modelling system (HiPIMS) | |
CA2955920C (en) | Optimizing multistage hydraulic fracturing design based on three-dimensional (3d) continuum damage mechanics | |
JP2014535044A5 (ru) | ||
GB2534990A (en) | A simulation-to-seismic workflow construed from core based rock typing and enhanced by rock replacement modeling | |
RU2015143556A (ru) | Способы и системы сопоставления истории месторождений для улучшенной оценки продуктивности месторождений | |
RU2012152638A (ru) | Снижение артефактов при итерационной инверсии геофизических данных | |
RU2011149638A (ru) | Системы, компьютерно-реализуемые способы и компьютерно-считываемые программные продукты для расчета приближенного давления дренирования скважины для имитатора коллектора | |
Zhang et al. | Simultaneous estimation of relative permeability and capillary pressure for tight formations using ensemble-based history matching method | |
RU2016107165A (ru) | Способ и система калибровки статической модели недр | |
KR102476935B1 (ko) | 가상그리드를 이용한 응력확대계수 측정시스템 및 측정방법 | |
RU2014109943A (ru) | Система и способ для определения характеристик геологической среды, включая оценивание неопределенности | |
Capozzoli et al. | A comparison of fast marching, fast sweeping and fast iterative methods for the solution of the eikonal equation | |
RU2015142104A (ru) | Бессеточное моделирование условий речно-дельтовой среды | |
Yu et al. | Consistent initial conditions for the Saint-Venant equations in river network modeling | |
US11300706B2 (en) | Designing a geological simulation grid | |
RU2015134392A (ru) | Способ моделирования подземного объема | |
MX2019003730A (es) | Modelado estructural mejorado. | |
US10705235B2 (en) | Method of characterising a subsurface volume | |
MX350240B (es) | Determinación de vía de migración de fluido. | |
Benoit et al. | Fully nonlinear and dispersive modeling of surf zone waves: Non-breaking tests | |
RU2611892C1 (ru) | Способ трехмерного моделирования заданного гидрогеологического объекта, реализуемый в вычислительной системе | |
CN112231878B (zh) | 夹层影响下垂向渗透率粗化等效计算方法及装置 | |
US20190339415A1 (en) | Methods of determining front propagation within a subsurface volume | |
RU2012151408A (ru) | Способ определения положения координат глубин маркера при построении геологической модели месторождения |