EA001212B1 - Способ сейсмического мониторинга нефтегазовых месторождений с использованием разновременных съемок - Google Patents
Способ сейсмического мониторинга нефтегазовых месторождений с использованием разновременных съемок Download PDFInfo
- Publication number
- EA001212B1 EA001212B1 EA199900195A EA199900195A EA001212B1 EA 001212 B1 EA001212 B1 EA 001212B1 EA 199900195 A EA199900195 A EA 199900195A EA 199900195 A EA199900195 A EA 199900195A EA 001212 B1 EA001212 B1 EA 001212B1
- Authority
- EA
- Eurasian Patent Office
- Prior art keywords
- data
- seismic
- trace
- sets
- parameters
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 135
- 238000012544 monitoring process Methods 0.000 title description 11
- 238000012545 processing Methods 0.000 claims abstract description 94
- 238000001228 spectrum Methods 0.000 claims abstract description 46
- 230000008569 process Effects 0.000 claims abstract description 21
- 238000012935 Averaging Methods 0.000 claims abstract description 20
- 238000013508 migration Methods 0.000 claims abstract description 14
- 230000005012 migration Effects 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000012937 correction Methods 0.000 claims description 17
- 238000003491 array Methods 0.000 claims description 12
- 238000007781 pre-processing Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 5
- 238000013461 design Methods 0.000 claims description 2
- 230000001629 suppression Effects 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 claims 1
- 230000001419 dependent effect Effects 0.000 abstract description 5
- 238000012986 modification Methods 0.000 abstract 4
- 230000004048 modification Effects 0.000 abstract 4
- 230000001131 transforming effect Effects 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 238000005311 autocorrelation function Methods 0.000 description 8
- 239000012530 fluid Substances 0.000 description 8
- 239000011148 porous material Substances 0.000 description 5
- 230000003068 static effect Effects 0.000 description 5
- 239000011435 rock Substances 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000005316 response function Methods 0.000 description 3
- 230000001932 seasonal effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical group [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 241001486234 Sciota Species 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002542 deteriorative effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000012958 reprocessing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 102220036542 rs144706057 Human genes 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 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/288—Event detection in seismic signals, e.g. microseismics
-
- 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
-
- 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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/612—Previously recorded data, e.g. time-lapse or 4D
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
- G01V2210/632—Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measuring Volume Flow (AREA)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
Abstract
Предлагается способ сравнения данных многократных сейсмических съемок резервуара, когда первая и вторая сейсмические съемки были проведены в разное время для обнаружения изменений в резервуаре в промежуток времени между съемками. В общих чертах способ включает расчет корректирующей функции для данных второй съемки на основе сравнения данных от не подвергшихся изменениям участков резервуара. Предлагается также способ обработки данных каждой съемки с использованием одинакового графа обработки.
Description
Область техники
Данное изобретение относится к области эксплуатации нефтяных и газовых месторождений, и, более конкретно, к обнаружению физических изменений на географической площади, в частности в резервуаре, на основании данных, полученных при сейсмическом мониторинге с использованием разновременных съемок.
Уровень техники
Знание характеристик нефтегазовых месторождений и контроль за изменением их параметров являются важным фактором при эксплуатации месторождений и добыче углеводородов. Эффективная эксплуатация месторождений является главной целью компаний, производящих энергоресурсы, поскольку они стараются уменьшить поисковые расходы, оптимизировать расположение буровых скважин и увеличить финансовую прибыль. Одним из методов, который пытается решить эту проблему, является сейсмический мониторинг с использованием разновременных съемок (также известный как 4Л). По мере извлечения нефти или газа, выработки месторождения, или его заполнения при закачке воды или газа в процессе эксплуатации или восстановления месторождения, появляются изменения эффективных упругих свойств пород в резервуаре. Возможность контролировать изменения в резервуаре как функции времени при помощи сейсмических методов позволит улучшить расположение эксплуатационных и нагнетательных скважин, даст возможность локализовать невыработанные зоны и более эффективно обслуживать месторождение, увеличивая, таким образом, общую стоимость продуктивного участка.
В двумерном варианте возможности сейсмомониторинга исследовались при проведении межскважинного просвечивания. Однако было проведено только качественное сравнение повторных результатов, на основании которого был сделан вывод, что данные двумерных разновременных съемок пока не содержат той информации, которая необходима при современных методах разработки месторождений (смотри, например, работу РаиИкоп е! а1., 1994, Т11С Ьеайтд Ейде, включенную в описание изобретения посредством ссылки на нее). 3Л технологии также опробовались, при разновременных съемках, но они содержат сложные моделирующие процедуры и требуют больших объемов обработки, не используя при этом реальные данные, получаемые в процессе проведения съемок (смотри, например, патент США № 4969130, включенный в описание изобретения путем ссылки на него).
Два конкретных примера применения разновременного сейсмического мониторинга, включенные в описание изобретения путем ссылки на них, опубликованы в работе Огеауек Ри1р (Сеор11у81С8. 1987), авторы которой контролировали выгорание пластов в Западном Техасе, и в работе Шеп апй Лии1ор (РпЫ Вгеак, 1993), описывающей контроль газового заполнения в ОкеЬегд Р1е1й, расположенном в норвежском секторе Северного моря. Еще одна статья, включающая материалы по паровому заполнению в Лиг: Р1е1й в Индонезии, была представлена Ьит1еу е! а1. на конференции 8оае1у οί Ехр1ога!юп Сеор11у8Ю8 (8ЕО) в Хьюстоне 1995 г. (она может быть получена через 8ЕО).
Одна из проблем при обработке данных сейсмического мониторинга с использованием разновременных съемок заключается в том, что изменяется не только состояние месторождения (резервуара), но и многие условия при проведении разновременных съемок, например, расположение источника и приемника при повторных работах, несомненно, будет отличаться от первоначальных. Кроме того, прилив при повторной морской съемке может быть выше или ниже, так же как и температура воздуха и воды. Аналогично, конкретные характеристики источников и приемников, используемых при повторных работах, также будут отличаться. Другие различия, кроме изменений в состоянии резервуара, также встречаются, например, разница в способах обработки данных, полученных при разновременных съемках. Таким образом, существует необходимость в разработке такого способа обработки данных разновременных съемок, чтобы различия в способах обработки не оказывали вредного влияния на результаты сравнения.
Например, при получении сейсмических данных используется источник для возбуждения сейсмических волн, которые отражаются от геологических границ в земле (например, границ пластов) и принимаются приемниками. В некоторых случаях сигнал источника считается кратковременным импульсом, хотя в действительности форма сигнала может существенно отличаться от идеальной. В процессе распространения через горные породы и отражения от сейсмических границ форма сигнала изменяется, и поэтому отраженный сигнал, регистрируемый приемниками, больше не является кратковременным импульсом или даже похожим на него. Деконволюция является процедурой, при помощи которой форма отраженного сигнала отбеливается с целью восстановления формы кратковременного импульса в сейсмических данных.
В другом примере используется широкополосный, ограниченный частотной полосой сигнал, который является нуль-фазовым. Деконволюция в этом случае используется для удаления искажений, связанных с распространением сигнала в земле.
Еще в одном примере, при выполнении деконволюции в частотной области, все частотные составляющие умножаются на определенные коэффициенты для приведения их к одинаковому уровню, следуя предположению, что сигнал, возбуждаемый в источнике, является минимально-фазовым, т.е. мгновенно нарастает до максимального значения и затем также быстро спадает до нуля. Для выполнения такой деконволюции во временной области производится вычисление функции автокорреляции сейсмической трассы, которая имеет форму симметричного волнового импульса. Энергетический спектр этого импульса анализируется для определения коэффициентов, на которые необходимо умножить каждую частотную составляющую импульса для выравнивания его частотного спектра. Эта операция выполняется в определенных окнах как вдоль сейсмической трассы (временные окна), так и поперек (интервал профиля) сейсмограммы (здесь термин сейсмограмма используется для альтернативных совокупностей записей, например, для сборок записей с общим пунктом приема (сейсмограммы ОПП), общей средней точкой (сейсмограммы ОСТ), общим пунктом возбуждения (сейсмограммы ОПВ), временных разрезов, полученных после суммирования трасс и т.д.). Автокорреляционная функция вычисляется в различных окнах, затем определяется среднее значение полученных результатов с целью последующего вычисления спектра. На основе этого спектра определяется оператор, необходимый для его сглаживания и формирования горизонтальной характеристики. Этот оператор затем используется при обработке всех входных трасс. Обычно временное окно выбирается равным 10 длинам создаваемого оператора и измеряется числом выборок. Операция деконволюции и расчет оператора деконволюции хорошо известны в сейсморазведке и не ограничиваются рассмотренным выше примером выполнения ее в частотной области. Она также часто реализуется и во временной области. Смотри, например, Υίΐιηαζ 1пуекИдайоик ίη Оеорйуыск Уо1. 2, 8е1кт1с Эа1а Ргосеккшд, 8ос1е1у о! Ехр1ота1юи Оеорйущскк (1987) и приведенные в этом источнике ссылки.
При выполнении деконволюции очень важно рассчитать оператор деконволюции с использованием данных конкретной съемки для учета особенностей формы возбуждаемого источником сигнала и других аппаратурных искажений сейсмических импульсов. Поэтому полученные при каждой съемке данные исправляются путем применения специального, оптимального именно для этой съемки, оператора деконволюции, который может не подходить для обработки данных других съемок. В результате использования различных операторов деконволюции при обработке данных разновременных съемок на разностной записи при вычитании сейсмограмм, полученных при этих съемках, будут появляться какие-то структуры. Этот результат является нежелательным. Однако до сих пор никто не предложил практического решения этой задачи.
Сущность изобретения
Задачей, решаемой настоящим изобретением, является преодоление рассмотренной выше проблемы.
Было установлено, что вопреки прежним убеждениям можно использовать один оператор деконволюции для обработки множества сейсмических данных, не только не ухудшая результаты, но и повышая качество обработки, проводимой для сравнения данных разновременных съемок. В соответствии с одним аспектом настоящего изобретения, предлагается способ получения разностной записи для данных, по меньшей мере, двух сейсмических съемок, включающий подготовку данных, по меньшей мере, двух сейсмических съемок на одной и той же географической площади и осуществление процесса деконволюции полученного множества сейсмических данных. Процесс деконволюции согласно изобретению включает: расчет оператора деконволюции на основе использования, по меньшей мере, двух совокупностей сейсмических данных, причем не менее двух совокупностей сейсмических данных, которые были зарегистрированы в разное время; использование этого оператора деконволюции при обработке, по меньшей мере, двух совокупностей сейсмических данных; и проведение дальнейшей обработки данных разновременных съемок для получения разностной записи.
Согласно еще одному из аспектов изобретения предлагается способ обнаружения физических изменений в месторождении (резервуаре), включающий подготовку данных, по меньшей мере, двух сейсмических съемок резервуара, первая из которых проведена в первый период времени, а вторая - во второй период времени, и осуществление процесса сравнения данных первой и второй съемок. При этом проведение сравнения данных разновременных съемок предусматривает: выделение первой группы отраженных сигналов (например, волновых импульсов) в совокупности сейсмических данных, полученных при первой съемке, и имеющей соответствующую вторую группу отраженных сигналов в совокупности сейсмических данных, полученных при второй съемке, причем первая и вторая группы отраженных сигналов характеризуют неизмененную часть геологической структуры или находящихся вблизи резервуара пластов, при этом первая группа отраженных сигналов представлена первой системой параметров, а вторая группа отраженных сигналов второй системой параметров. Затем определяется значение допустимой пороговой разностной функции для сейсмических данных с первой и второй системами параметров с последующим определением корректирующей функции для обработки данных со второй системой параметров.
Далее предусматривается, что корректирующую функцию вычисляют таким образом, что после ее применения к данным со второй системой волновых параметров определяется скорректированная система параметров, и разность между данными с первой и скорректированной системами волновых параметров будет ниже пороговой разностной функции. Затем корректирующую функцию применяют при обработке третьей группы отраженных сигналов, которая выделяется на сейсмических данных, полученных при второй съемке, при этом формируют скорректированную третью группу отраженных сигналов, которой соответствует четвертая группа отраженных сигналов на совокупности сейсмических данных первой съемки, причем третья и четвертая группы отраженных сигналов представляют изменяющуюся часть резервуара.
Сравнение скорректированной третьей группы отраженных сигналов с четвертой группой отраженных сигналов производится путем вычитания скорректированной третьей группы отраженных сигналов из четвертой, которое и дает желаемую информацию.
Согласно более конкретному варианту, приводимому в качестве примера, указанное выше определение допустимой пороговой разностной функции включает:
итеративный подбор измененных параметров группы отраженных сигналов для второй системы параметров, применение измененных параметров при обработке данных со второй системой параметров, при этом формируется новая запись с измененной системой параметров, сравнение данных с измененной системой параметров с данными с первой системой параметров, причем указанный итеративный подбор продолжают до тех пор, пока не будет достигнуто схождение результатов сравнения. При этом допустимая пороговая разностная функция включает измененную систему параметров группы отраженных сигналов, при которой было получено схождение. Примерная система параметров включает любую комбинацию амплитуды, фазы, частотной полосы, времени вступления или в отдельности любые из вышеприведенных параметров.
Согласно другому варианту осуществления изобретения, приводимому в качестве примера, определение корректирующей функции включает: итеративный подбор измененных параметров группы отраженных сигналов для второй системы параметров, применение измененных параметров при обработке данных со второй системой параметров волны, при этом формируется измененная система параметров, сравнение данных с измененной системой параметров с данными с первой системой параметров и определение допустимой пороговой разности, причем указанный итеративный подбор продолжают до тех пор, пока в результате проведения указанной операции сравнения, являющейся разностью данных с первой системой параметров с данными с измененной системой параметров, она не станет ниже допустимой пороговой разности.
Согласно еще одному варианту, приведенному в качестве примера, указанное определение допустимой пороговой разностной функции включает: определение разности между записью сигналов во временном окне первой трассы из совокупности сейсмических данных, полученных при первой съемке, и записью сигналов во временном окне второй трассы, полученной при проведении второй сейсмической съемки, причем вторая трасса включает отраженные сигналы, соответствующие отраженным сигналам на первой трассе, а временное окно второй трассы соответствует временному окну первой трассы; вычисление отношения разности сигналов в выделенных окнами участках первой и второй трасс к сигналу во временном окне первой трассы; выбор допустимого значения пороговой разности, которое было бы меньше этого отношения.
Сигналы во временных окнах, соответствующих как постоянным, так и изменяющимся частям резервуара, имеют похожие спектральные характеристики. Например, если сейсмические сигналы, приходящие от пластов в изменяющейся части резервуара, имеют преобладающую частоту, равную 30 Гц, то временное окно должно быть выбрано там, где преобладающая частота сейсмических сигналов от неизмененной части объекта исследований ближе всего к 30 Гц. Аналогично, фазовые изменения сейсмических сигналов, приходящих из резервуара и неизмененной части месторождения, должны быть возможно более близкими.
Предпочтительнее, однако, ошибаться в более высокочастотную сторону. Например, если преобладающая частота сигналов из резервуара равна 30 Гц, то окно, содержащее сигналы с частотой 35 Гц, считается более подходящим, чем окно с сигналами, имеющими преобладающую частоту в 25 Гц. Установлено, что при расхождении преобладающих частот не более чем на 25%, результаты получаются адекватными. Наилучшие результаты будут получены, когда расхождение в частотных полосах отраженных сигналов не превышает 10%.
В еще одном варианте указанное определение допустимой пороговой разностной функции включает: определение разности между квадратом записи сигналов во временном окне первой трассы из данных первой съемки и квадратом записи сигналов во временном окне второй трассы из данных второй съемки, причем вторая трасса содержит отраженные сигналы, соответствующие отраженным сигналам на первой трассе, и временное окно второй трассы соответствует временному окну на первой трассе; вычисление отношения разности квадратов записей сигналов в выделенных временными окнами участках записи первой и второй трасс к квадрату записи сигналов во временном окне первой трассы; и выбор допустимого порогового значения разности меньшим, чем это отношение.
Согласно другому варианту указанное применение корректирующей функции к третьей группе отраженных сигналов во второй совокупности сейсмических трасс включает свертку корректирующей функции с третьей группой отраженных сигналов во второй совокупности трасс; в этих вариантах указанная первая совокупность сейсмических данных включает трассу от сейсмоприемника. В альтернативном варианте указанная первая и вторая совокупности сейсмических данных включают просуммированные сейсмические трассы от сейсмоприемников, т.е. суммограммы трасс (временные разрезы) или сейсмограммы ОСТ (общая средняя точка) и ОПВ (общий пункт возбуждения). В других альтернативных вариантах указанная первая и вторая совокупности сейсмических данных включают сейсмограммы после предварительной обработки, предшествующей построению временных разрезов, или мигрированные данные.
Во многих вариантах указанная первая и вторая совокупности сейсмических данных подвергаются эквивалентной предварительной обработке, например, кроме описанной выше деконволюции, в некоторых вариантах первая и вторая совокупности данных имеют одинаковую процедуру восстановления формы сигнала, одинаковый граф обработки для подавления помех и многократных волн. Кроме того, во многих вариантах используют одинаковые ΌΜΟ операторы, а миграцию в первой и второй совокупностях сейсмических данных проводят с использованием эквивалентной скоростной модели среды.
И, наконец, согласно еще одному аспекту изобретения предлагается способ формирования сигнала, характеризующего географическую площадь по результатам проведения разновременных сейсмических съемок. Данный способ включает подготовку данных по меньшей мере двух сейсмических съемок этой географической площади, первая из которых проведена в первый период времени, а вторая - во второй период времени, и обработку этих сейсмических данных. Указанная обработка включает: выполнение определенного графа обработки входных данных, полученных при первой съемке, выполнение этого же графа обработки входных данных, полученных при второй съемке, причем граф обработки зависит от системы параметров сейсмических сигналов; выбор не менее одной системы параметров путем их итеративного подбора на основе данных обеих съемок, и использование не менее одной системы параметров по меньшей мере в одном графе при обработке данных, полученных как при первой, так и при второй съемках.
Краткий перечень фигур чертежей
Для более полного понимания сущности настоящего изобретения и связанных с его применением преимуществ в следующем разделе описания будут сделаны ссылки на прилагаемые графические материалы, на которых:
Фиг. 1А является изображением первой записи сейсмической съемки, проведенной в первой раз без применения деконволюции.
Фиг. 1В является изображением первой записи, показанной на фиг. 1А, с проведенной деконволюцией, причем оператор деконволюции был рассчитан на основе использования данных первой записи.
Фиг. 1С является изображением первой записи, показанной на фиг. 1А, с выполненной деконволюцией, причем оператор деконволюции рассчитывался на основе данных первой записи и данных другой съемки.
Фиг. 2А является изображением второй записи второй съемки, проведенной на той же самой географической площади, что и первая съемка, но в другое время; деконволюция при получении этой записи не применялась.
Фиг. 2В является изображением второй записи, приведенной на фиг. 2А, с выполнением деконволюции, причем оператор деконволюции был рассчитан на основе данных второй записи.
Фиг. 2С является изображением второй записи, показанной на фиг. 2А, с применением деконволюции, причем оператор деконволюции был рассчитан на основе использования данных первой и второй записей, приведенных на фиг. 1А и фиг. 2 А.
Фиг. ЗА является изображением разности записей, приведенных на фиг. 1А и 2С.
Фиг. ЗВ является изображением разности записей, приведенных на фиг. 1В и 2В.
Фиг. ЗС является изображением разности записей, приведенных на фиг. 1С и 2С.
Следует, однако, отметить, что приложенные материалы иллюстрируют только наиболее типичные варианты осуществления данного изобретения и поэтому не могут считаться ограничивающими его объем, так как изобретение допускает другие, в равной степени эффективные варианты его осуществления.
Сведения, подтверждающие возможность осуществления изобретения
Применяемый здесь термин корректирующая обработка является универсальным термином, объединяющим согласующую фильтрацию, выравнивание амплитуд и введение статических поправок, необходимых для проведения разновременных региональных сейсмических наблюдений (сейсмического мониторинга). По существу, определяют оператор или операторы волнового импульса для формирования и согласования отраженных сигналов, полученных при разных съемках. Обычно операторы рассчитывают для сигналов от не под вергшихся изменению отражающих горизонтов, находящихся вне резервуара, в коллекторах которого могут происходить значительные изменения (в порах) состояния флюидов (т.е. текучих сред, жидкости или газа). Теоретически, разность между двумя массивами данных после корректирующей обработки должна равняться нулю везде (внутри и вне выбранного окна), за исключением тех участков в резервуаре, где происходят изменения. Все отраженные сигналы от неизмененных (статических) или находящихся вне резервуара границ должны быть полностью вычтены, останутся только изменения, соответствующие динамическим участкам резервуара (т.е. местам, где происходит замещение флюидов внутри пор).
Возможность обнаружения изменений в состоянии резервуара сейсмическим методом определяется его литологией и акустическими свойствами поровых флюидов (которые зависят от температуры и давления), типа эксплуатационного или восстановительного процесса, разрешающей способности (пространственной и временной) сейсмических данных, не считая требований к повторяемости сейсмической технологии.
В некоторых резервуарах с газовым нагнетанием введение газа может уменьшить акустический импеданс (нефть замещается газом) в достаточной степени для появления яркого пятна (смотри также статью Ри1р и Сгсаусх. 0Е0РНУ81С8, 1987. включенную в описание изобретения посредством ссылки на нее, где рассматривается возникновение яркого пятна при выгорании в пластах). В других случаях заполнение газом или водой (нагнетание воды) может привести к противоположной реакции и образовать тусклое пятно. Там, где производится заполнение паром, может наблюдаться снижение скорости упругих колебаний, которое часто используется для локализации фронта пара в породах резервуара, с целью определения, до какой степени продвинулось восстановление (смотри, например, статью Ьит1еу с1 а1. 8ЕО Ехраибеб ЛЬ51гас15. 1995, включенную в описание изобретения посредством ссылки на нее). Из вышеизложенного следует, что необходима петрофизическая и сейсмическая интерпретация для того, чтобы понять, какие характерные свойства являются существенными и какие изменения в породах резервуара можно обнаружить.
Обработка
В рассматриваемых далее примерах применялся стандартный граф обработки, выполняемый с целью получения пространственно скорректированных записей полезных сигналов, с сохранением относительных амплитуд входных данных, который включает ΝΜ0/ΌΜ0 процедуру и предварительную (с нулевым сдвигом) миграцию. Другие процедуры предварительной обработки опытными специалистами могут использоваться и другие процедуры предварительной обработки, включая уплотнение данных, а также способы выборки данных для уменьшения трудоемкости вычислительных операций и сохранения качества предварительной обработки. После такой обработки будут получены ОСТ сейсмограммы с сохранением относительных амплитуд и предварительной миграцией, для которых выполняют корректирующую обработку и разностные измерения данных, полученных при разновременных съемках. В зависимости от характеристик эксплуатируемого резервуара, согласно альтернативному варианту осуществления изобретения, ОСТ, ЛУ0 свойства, а также временные разрезы с миграцией или без нее, полученные при обработке данных каждой съемки, подвергают корректирующей обработке и затем вычитают (определяют разность) для обнаружения внутрипорового замещения флюидов.
Рекомендуется, чтобы данные каждой разновременной съемки обрабатывались возможно более близкими способами (что обычно не выполняется на практике). До настоящего изобретения использовались любые обработанные данные, находящиеся в обрабатывающих центрах (с отличающимися процедурами и алгоритмами предварительной обработки), а удаление имеющихся в результате этого различий в сейсмических данных пытались устранить путем корректирующей обработки. Лучшая корректирующая обработка и более существенные значения разности сейсмических данных будут получены при повторной обработке имеющихся сейсмических данных с той же самой методологией и программным обеспечением, что и последняя съемка в сейсмическом мониторинге. Таким образом, предпочтительнее выполнять операции восстановления формы сигнала, деконволюции, подавления помех и кратных волн с одинаковыми параметрами обработки для каждой съемки, точно также важно использовать одинаковые ΌΜ0 операторы и выполнять миграцию данных с применением эквивалентной скоростной модели для получения результатов, пригодных для проведения дальнейших измерений.
В связи с этим, согласно одному аспекту изобретения, предлагается способ формирования сигнала, характеризующего географическую площадь по результатам проведения разновременных сейсмических съемок. Способ включает подготовку данных, по меньшей мере, двух сейсмических съемок этой географической площади, первая из которых проведена в первый период времени, а вторая - во второй период времени, и обработку этих сейсмических данных, в процессе которой выполняют граф обработки входных данных, полученных при первой съемке, граф обработки входных данных, полученных при второй съемке, причем граф обработки определяется заданной систе мой параметров сейсмического сигнала. Данный способ включает: (а) выбор не менее одной системы параметров путем их итеративного подбора на основе данных обеих съемок и (б) использование не менее одной системы параметров, по меньшей мере, в одном графе при обработке данных, полученных как при первой, так и при второй съемках.
Например, было обнаружено, что улучшенная деконволюционная обработка данных обеих съемок достигается при расчете деконволюции оператора на основе использования множества совокупностей данных. Этот оператор деконволюции применяется затем при обработке данных обеих съемок. Такой оператор может не являться оптимальным для обработки данных любой отдельной съемки, тем не менее, он будет оптимальным при сравнении результатов разновременных съемок.
Предлагается также способ получения разностной записи для данных, по меньшей мере, двух сейсмических съемок, включающий подготовку данных, по меньшей мере, двух сейсмических съемок на одной и той же географической площади и осуществление процесса деконволюции полученных данных. При этом процесс деконволюции предусматривает расчет оператора деконволюции на основе использования, по меньшей мере, двух совокупностей сейсмических данных, причем не менее двух совокупностей сейсмических данных были зарегистрированы в разное время; использование этого оператора деконволюции при обработке, по меньшей мере, двух совокупностей сейсмических данных; и проведение дальнейшей обработки данных разновременных съемок для получения разностной записи. В некоторых случаях имеются не менее трех совокупностей сейсмических данных, причем каждая из них представлена записями, проведенными в разное время, при этом конструкция оператора деконволюции зависит от данных, полученных при каждой съемке. В некоторых из этих случаев выполняются дополнительные процедуры обработки, включающие расчет второго оператора деконволюции на основе использования данных, полученных при первой и третьей съемках, причем второй оператор деконволюции используется при обработке данных первой и третьей съемок. С одинаковой эффективностью в качестве первой совокупности данных может быть принята совокупность, взятая до или после второй совокупности.
Согласно более конкретному варианту, приводимому в качестве примера, расчет первого оператора деконволюции включает: осреднение спектра мощности для первой совокупности трасс данных первой съемки; осреднение спектра мощности для второй совокупности трасс данных второй съемки; осреднение средних значений спектров мощности первой и второй совокупностей трасс; и расчет оператора декон волюции для обработки данных первой и второй съемок на основе осредненных значений спектров мощности первой и второй совокупностей сейсмических трасс.
В некоторых обстоятельствах рассмотренное выше осреднение спектров мощности для каждой последующей съемки может оказаться чересчур дорогим. Поэтому в соответствии с другим вариантом осуществления изобретения, при обработке данных многократных съемок запоминаются осредненные значения данных первой съемки, так что когда производится более поздняя съемка, автокорреляционные функции данных первой съемки не должны вычисляться заново. Согласно этому варианту расчет первого оператора деконволюции включает:
осреднение спектра мощности для первой совокупности трасс данных первой съемки;
обратное преобразование усредненного значения с получением функции, отображающей усредненный спектр мощности первой совокупности трасс во временной области; запоминание этой временной функции; усреднение спектра мощности второй совокупности трасс данных второй съемки; обратное преобразование усредненного значения с получением функции, отображающей усредненный спектр мощности второй совокупности трасс во временной области;
осреднение временных функций, представляющих средние значения спектров мощности первой и второй совокупностей трасс; и расчет оператора деконволюции для обработки данных первой и второй съемок с использованием усредненного значения временных функций, отображающих средние значения спектров мощности первой и второй совокупностей.
Рассмотрим теперь фиг. 1А, иллюстрирующую конкретный пример совокупности сейсмических трасс, который будет далее обсуждаться. На фиг. 1А изображен первый предварительный временной разрез (без ΌΜΘ, миграции и различных способов подавления помех, только учет сферического расхождения и геометрическая коррекция), построенный по данным первой съемки без выполнения процедуры деконволюции. На фиг. 1В изображен временной разрез, приведенный на фиг. 1А, после проведения деконволюции, причем оператор деконволюции был рассчитан на основе данных предварительного временного разреза и использован при обработке сейсмических данных до суммирования трасс. Были использованы следующие параметры:
Дискретность 2 мс
Длина оператора 140 мс
Интервал предсказания единичный импульс 2 мс
Окно анализа зависит от расстояния источник-приемник, для ближних расстояний 3005000 мс, для дальних 3600-5000 мс
Добавленный белый шум 0,5%
Окно для проведения деконволюции 0-600 мс.
На фиг. 2А изображен второй предварительный временной разрез, построенный по данным второй съемки, проведенной на той же географической площади, что и первая съемка, но в другое время; при построении разреза не выполнялась операция деконволюции; на фиг. 2В показан временной разрез, изображенный на фиг. 2А, после проведения деконволюции, причем оператор деконволюции рассчитывался на основе данных второго предварительного временного разреза и использовался до суммирования трасс. Здесь была использована та же система параметров, которая применялась для обработки данных первой съемки.
Фиг. ЗА является разностью временных разрезов, приведенных на фиг. 1А и 2А. Фиг. ЗВ представляет разность временных разрезов, изображенных на фиг. 1В и 2В. В идеальном случае на разностном разрезе, изображенном на фиг. ЗА, не должно быть ничего, кроме шумовых помех в области, соответствующей неизмененной части геологического разреза, если в остальном обе съемки были идентичны. Однако можно видеть, что неизмененные структуры проявляются на разностном разрезе, изображенном на фиг. ЗВ.
На фиг. 1С приведен первый временной разрез, показанный на фиг. 1А, после проведения деконволюции до построения временного разреза, причем оператор деконволюции был рассчитан на основе использования данных первой и второй совокупностей трасс. На фиг. 2С приведен второй временной разрез, изображенный на фиг. 2А, после деконволюции, проведенной до построения временного разреза, причем оператор деконволюции рассчитывался на основе данных первой и второй совокупностей трасс. Применявшаяся система параметров для обработки данных на фиг. 1С и 2С была такой же, как и для данных на фиг. 1А.
На фиг. ЗС изображена разность временных разрезов, приведенных на фиг. 1С и 2С. Отметим отсутствие неизмененных структур.
Для расчета деконволюционного фильтра используется несколько вариантов;
в одном из них обрабатывается первый (т.е. базовый) массив входных данных, причем осредненные значения функций автокорреляции сейсмических трасс внутри каждой совокупности трасс (т.е. сейсмограммы ОПВ, ОПП, ОСТ и т.д.) используют путем преобразования в частотной области для оценки спектра мощности.
ί=0
Ρ(ζΛ) = ΣΡ0·^)’7 = 1’-β ί=1 где N - число сейсмических трасс по оси X, О - индекс частоты Найквиста,
Р(1,к) - оценка спектра мощности для индекса частоты ί и номера трассы к, и
_) - индекс разновременной съемки.
Затем производят обратное преобразование Фурье, и автокорреляционные функции трасс для каждой совокупности сохраняются. Такая же оценка спектра мощности производится для массива данных, полученных в другое время (т.н. контрольный массив), посредством вычисления _)=1 Р(1, к) где _)=1 является индексом первой контрольной съемки.
После этого все полученные для разновременных съемок усредненные автокорреляционные функции осредняются для вычисления новой автокорреляционной функции с использованием следующей формулы:
___ м рУ Авм
.)=1,...М, где М - номер разновременной съемки.
Полученные автокорреляционные функции используются для расчета частотных фильтров, которые применяются при обработке данных всех разновременных съемок.
Другие способы расчета деконволюционных фильтров (например, переменная во времени деконволюция, деконволюция, учитывающая поверхностные условия, деконволюция в общей области (ОПВ, ОПП, ОСТ, или равные удаления источник-приемник)), которые могут применяться опытными специалистами, также основаны на использовании данных обеих съемок, и не выходят за рамки настоящего изобретения.
Далее, аналогично рассмотренной выше операции деконволюции, параметры других процедур обработки для реализации преимуществ настоящего изобретения также должны рассчитываться на основе данных, полученных при всех съемках, а не какой-либо одной. Например, корректирующая обработка (более детально будет обсуждаться ниже), ί-к частотная фильтрация, и фильтрация в Радон-области, вычисление статических поправок, удаление кратных волн в РК, в τ-ρ и в Радон-областях.
На практике имеются различия в условиях регистрации сейсмических данных, независимо от того, насколько хорошо эти работы были спланированы или выполнены. В идеальном случае, когда регистрация данных при выполнении базовой и контрольной съемок проводится в одинаковых условиях, обработка также должна быть идентичной. В реальных условиях это не выполняется. Изменения в условиях регистрации данных при проведении разных съемок могут появиться в результате использовании различных регистрирующих систем, которые применяются при реализации настоящего изобретения (например, постоянно установленные системы, ΘΒ3 системы и морские сейсмические косы). Такие изменения могут быть вызваны отклонением косы от линии профиля, смещени ем кабеля, азимутом воздействий, изменениями в источниках, состоянием моря в момент проведения наблюдений и т.д., в результате которых могут появиться различия в условиях проведения съемок.
Корректирующая обработка
Согласно одному из возможных вариантов осуществления изобретения имеются четыре корректирующих элемента в используемой корректирующей функции. Это введение временных поправок, выравнивание среднеквадратичного значения энергии, нормализация частотной полосы и фазовое согласование. Каждый элемент является существенным при построении переходной функции или функции импульсной реакции при обработке данных двух съемок или, альтернативно, двух трасс каждой съемки. Результирующая трасса после корректирующей обработки (ΐΧΕρ) вычисляется следующим образом:
ΙχΕΟ = I 1 (^СОГГ, Г1П5соГ|. тСО|т- Рсогг) (1) где I - входная трасса, * - обозначает свертку, а 5СОГГ, тш5согг, тсогг и рсогг являются элементами корректирующей обработки (в функции импульсной реакции 1), соответствующими времени, амплитуде, частотной полосе и фазе.
В соответствии с альтернативным вариантом осуществления изобретения функции импульсной реакции вычисляют последовательно, с учетом поверхностных условий и с использованием данных основной съемки (база) и данных контрольной съемки (повторение). Согласно настоящему варианту, приведенному в качестве примера, данные контрольной съемки преобразуются (в результате корректирующей обработки) таким образом, чтобы соответствовать данным базовой съемки, и функции импульсной реакции рассчитывают для горизонтальных временных окон, расположенных в интервале записи вне зоны резервуара, в которой ожидаются изменения.
Далее влияние каждой составляющей будет пояснено на простых моделях и на примере разновременных исследований при выполнении трехмерного проекта.
Элементы корректирующей обработки.
Временные поправки (5СОГГ).
Кроме бурного моря в шторм и течения, которые могут быть причиной отклонения косы от линии профиля, также играют роль приливно-отливные условия и температурные изменения. Эти изменения могут быть несущественными при регистрации, обработке или интерпретации сейсмических данных для одиночной 3Ό съемки. Однако, когда несколько 3Ό съемок используют для сейсмического мониторинга, существуют потенциальные предпосылки для непостоянства состояния моря при разных съемках.
Сезонные и кратковременные (связанные со штормом) изменения температуры и солености воды также могут служить причиной разли чия во временах прихода отражений, зарегистрированных при разных съемках. Степень изменения является функцией, зависящей от изменений температуры и солености воды на линии профиля, и глубины этих изменений в толще воды. Например, съемки, проведенные в Восточно-Китайском море, показывают, что в течение 10-летнего периода наблюдений сезонные изменения являлись ощутимыми, и они могли изменять (до 2 мс) время пробега отраженной волны от опорной сейсмической границы в течение сезонных максимумов (т.е. если проводить съемку в январе и июне) на участках, где глубина воды равна или превышает 100 м.
Что касается приливных изменений, экстремальный пример был бы получен при проведении сейсмических съемок в заливе Фанди, между канадскими атлантическими провинциями Нью Брансуик и Новая Шотландия. Здесь изменения прилива составляют около 15 м, что эквивалентно 20 мс двойного времени пробега сейсмической волны (принимая значение скорости в воде, равной 1500 м/с). Если одна съемка проведена при низком приливе, а вторая - при высоком (в предположении, что продолжительность съемок была очень короткой), то вычисление разности сейсмических данных без соответствующей временной коррекции, не даст удовлетворительных результатов.
В общем случае ошибки увеличиваются с увеличением задержки сейсмических сигналов, при этом использование более широкополосных сейсмических импульсов приведет к увеличению ошибок при вычислении разности даже при небольших задержках.
Выравнивание среднеквадратичных значений энергии (тт5СОГГ).
В наиболее совершенном эксперименте при осуществлении разновременного сейсмического мониторинга одна и та же сейсмическая партия должна проводить каждую контрольную съемку с точно таким же оборудованием и при одинаковых условиях (состояние моря) с ранее проведенной первоначальной базовой съемкой. Однако разница в условиях регистрации появляется даже в хорошо контролируемых ситуациях. Поэтому приведение данных сейсмических съемок к эквивалентному среднеквадратичному уровню необходимо, особенно, если она (среднеквадратичная энергия) значительно изменяется в период между съемками.
В одном примере было проанализировано около 100 сейсмограмм ОСТ, полученных при проведении разновременных морских 3Ό съемок с интервалом около 18 месяцев, геометрия приемных устройств и аппаратура были почти идентичными, за исключением того, что одна съемка проводилась в летние месяцы, а другая зимой. Анализ средних значений амплитуд сейсмических данных внутри дискретных временных окон длительностью 200 мс показал, что существует амплитудное смещение между данными основной и контрольной съемок по всему профилю. После применения сглаженной функции выравнивания амплитуд, входящей в состав корректирующей обработки внутри каждого временного окна, было отмечено, что расхождение между данными базовой и контрольной съемок уменьшилось.
Не вызывает сомнений, что такие изменения необходимо скорректировать прежде, чем разностные разрезы разновременных съемок могут быть проинтерпретированы.
Нормализация частотной полосы (тсо|1)
Уравнивание частотных полос сейсмических сигналов также является частью вычисления сейсмической разности. Если два сейсмических массива должны быть вычтены, то лучше это делать с сейсмическими данными, имеющими одинаковые спектры. Корректирующая обработка (если она достаточно надежна) должна исправить спектр одного сигнала таким образом, чтобы он согласовывался со вторым, но если корректирующая обработка не оптимизирована, появляются расхождения. Чем значительнее рассогласование спектров, тем больше ошибки. Тесты показывают, что остаточная энергия отраженных сигналов (т.е. энергия, которая остается после вычитания временных разрезов разных съемок и которая в идеальном случае должна быть полностью удалена), колеблется в пределах от 7% до 40% от амплитуды входных сигналов для импульсов с центральной частотой, различающейся на 2,5 Гц и 12,5 Гц (для сейсмического импульса Риккера с центральной частотой, смещавшейся дискретно через 2,5 Гц, и вычитавшегося из базового импульса Риккера с центральной частотой 30 Гц). Отклонения центральной частоты в диапазоне 12,5 Гц вполне приемлемы для самых ранних 3Ό съемок, но иногда встречаются большие отклонения частоты, и в этих случаях остаточная энергия отраженных сигналов увеличивается. Поэтому в соответствии с предлагаемыми вариантами обработки данных таких съемок выполняют корректирующую обработку для нормализации частотных спектров сейсмических сигналов с целью уменьшения возможных ошибок при определении разности сейсмических данных.
Фазовое согласование (рсогг)
И, наконец, ошибки могут быть также связаны с фазовым рассогласованием сейсмических сигналов. Например, для сейсмического 30 Гц импульса Риккера при повороте фазы на 5, 10, 15, 30, 60 и 90 градусов относительно нулевой фазы остаточная энергия, вызванная поворотом фаз, может достигать 20% от амплитуды входного сигнала при фазовом рассогласовании всего в 15 градусов. Кроме того, человеческий глаз может с трудом уловить фазовые различия в 15 градусов и меньше, поэтому важно, чтобы оператор корректирующей обработки исправлял небольшие фазовые расхождения так же, как и значительные. На ошибки, появляющиеся в результате фазового рассогласования, не влияют расхождения в частотных полосах сейсмических сигналов. Было отмечено, что остаточные ошибки в 20% могут появиться для сейсмических импульсов, различающихся по фазе на 15 градусов.
Разность сейсмических данных
Главная цель разновременного сейсмического мониторинга заключается в обнаружении изменений, которые связаны с внутрипоровым движением флюидов или отсутствии этого движения в интервале времени между съемками. В настоящее время это выполняется путем вычитания данных контрольной и базовой съемок после исправления данных в результате корректирующей обработки. На разностных данных с небольшим количеством остаточной энергии отраженных сигналов с более высокой вероятностью возможно обнаружение движения флюидов по сравнению с данными, имеющими высокий уровень остаточной энергии. Для иллюстрации эффективности корректирующей обработки будет приведено несколько разностных разрезов с корректирующей обработкой различных компонент, полученных при обработке данных разновременных съемок в Северном море. Мигрированные временные разрезы базовой и контрольной съемок, пересекающие площадь вне пределов месторождения, представлены на фиг. 4 и 5 соответственно, и преобладающая серия сильных статических отражений видна в интервале между 2000 и 2600 мс. В этом примере использовалось окно в интервале между 1900 и 2700 мс. Был разработан оператор длиной 800 мс, кроме того, применялась полосовая фильтрация 3/8 - 35/55 Гц для обоих массивов данных до корректирующей обработки. Поскольку это были отражения от статических отражающих границ (т.е. они являлись отраженными сигналами, представляющими неизменяющуюся часть геологической структуры внутри или около резервуара), то вычитание данных базовой и контрольной съемок должно было бы в идеальном случае дать очень маленькую остаточную энергию отраженных сигналов, так как обе съемки были проведены в идентичных условиях. Однако разностный разрез, полученный без применения какой-либо корректирующей обработки, показал, что вычитания нет. Остаточная энергия отраженных сигналов значительна. С таким большим количеством остаточной энергии, оставшимся на разностном разрезе, было бы трудно обнаружить любые перемещения флюидов (в других местах) в резервуаре.
Сравнение разностного разреза, полученного без корректирующей обработки, с разрезами, в которых корректировался только амплитудный спектр или только фазовый спектр, или оба этих параметра, показало, что фазовая коррекция является наиболее существенной состав ляющей в общей корректирующей процедуре, что подтверждается также результатами математического моделирования.
Рассмотренная выше обработка производилась на платформе с множеством параллельных процессоров (например, ΙΒΜ ΞΡ2, 1п1е1 Рагадоп) и с программным обеспечением, совместимым с таким оборудованием (например, Ρ6Ξ Тепйог δ СИВЕ ΜΑΝΆΟΕΕ ΤΜ).
Вышеприведенное описание дано только как пример, и другие варианты могут реализоваться опытными геофизиками без отклонения от объема изобретения, которое определяется приведенной ниже формулой изобретения.
Claims (45)
- ФОРМУЛА ИЗОБРЕТЕНИЯ1. Способ обнаружения физических изменений в резервуаре, включающий подготовку данных, по меньшей мере, двух сейсмических съемок резервуара, первая из которых проведена в первый период времени, а вторая - во второй период времени; и осуществление процесса сравнения данных первой сейсмической съемки с данными второй сейсмической съемки, предусматривающий выделение первой группы отраженных сигналов (например, волновых импульсов) в совокупности сейсмических данных, полученных при первой съемке, и имеющей соответствующую вторую группу отраженных сигналов в совокупности сейсмических данных, полученных при второй съемке, при этом:первая и вторая группы отраженных сигналов характеризуют неизмененную часть геологической структуры в резервуаре или вблизи него, первая группа отраженных сигналов представлена первой системой параметров, а вторая группа отраженных сигналов - второй системой параметров;определение допустимой пороговой разностной функции между первой и второй системами параметров волны;определение корректирующей функции для второй системы параметров, причем корректирующая функция характеризуется тем, что при приложении ее ко второй системе волновых параметров определяется скорректированная система параметров, при этом разность между данными первой и скорректированной систем волновых параметров будет ниже пороговой разностной функции;применение корректирующей функции при обработке третьей группы отраженных сигналов, которая связана с сейсмическими данными, полученными при второй съемке, при этом:формируется скорректированная третья группа отраженных сигналов, третьей группе отраженных сигналов соответствует четвертая группа отраженных сигна лов на совокупности сейсмических данных первой съемки, третья и четвертая группы отраженных сигналов представляют изменяющуюся часть резервуара;сравнение скорректированной третьей группы отраженных сигналов с четвертой группой путем вычитания третьей группы из четвертой.
- 2. Способ по п.1, отличающийся тем, что указанное выше определение допустимой пороговой разностной функции включает:итеративный подбор измененных параметров группы отраженных сигналов для второй системы параметров, применение измененных параметров при обработке данных со второй системой параметров, при этом формируется измененная система параметров, сравнение данных с измененной системой параметров с данными с первой системой параметров, при этом указанный итеративный подбор продолжается до тех пор, пока не будет достигнуто схождение результатов сравнения, и допустимая пороговая разностная функция включает измененную систему параметров группы отраженных сигналов, при которой было получено схождение.
- 3. Способ по п.2, отличающийся тем, что один из указанных параметров волны включает амплитуду.
- 4. Способ по п.2, отличающийся тем, что один из указанных параметров волны включает фазу.
- 5. Способ по п.2, отличающийся тем, что один из указанных параметров волны включает частотную полосу.
- 6. Способ по п.2, отличающийся тем, что один из указанных параметров волны включает время.
- 7. Способ по п.2, отличающийся тем, что указанные параметры волны включают амплитуду, фазу, частотную полосу и время.
- 8. Способ по п.1, отличающийся тем, что определение корректирующей функции включает итеративный подбор измененных параметров группы отраженных сигналов для второй системы параметров, применение измененных параметров при обработке данных со второй системой параметров волны, при этом формируется измененная система параметров, сравнение данных с измененной системой параметров с данными первой системы параметров и определение допустимой пороговой разности, причем указанный итеративный подбор продолжают до тех пор, пока в результате проведения указанной операции сравнения разность данных с первой системой параметров с данны ми с измененной системой параметров не станет ниже допустимой пороговой разности.
- 9. Способ по п.1, отличающийся тем, что указанное определение допустимой пороговой разностной функции включает определение разности между записью сигналов во временном окне первой трассы из совокупности сейсмических данных, полученных при первой съемке, и записью сигналов во временном окне второй трассы, полученной при проведении второй сейсмической съемки, при этом:вторая трасса включает отраженные сигналы, соответствующие отраженным сигналам на первой трассе, временное окно второй трассы соответствует временному окну первой трассы, и вычисление отношения разности сигналов в выделенных окнами участках первой и второй трасс к сигналу во временном окне первой трассы, и выбор значения допустимой пороговой разности меньшим, чем это отношение.
- 10. Способ по п.9, отличающийся тем, что разница в частотной полосе сигналов во временном окне первой трассы и временном окне второй трассы не превышает 25%.
- 11. Способ по п.10, отличающийся тем, что разница в частотной полосе сигналов во временном окне второй трассы и временном окне первой трассы менее 10%.
- 12. Способ по п.1, отличающийся тем, что указанное вычисление допустимой пороговой разностной функции включает определение разности между квадратом записи сигналов во временном окне первой трассы из данных первой съемки и квадратом записи сигналов во временном окне второй трассы из данных второй съемки, причем вторая трасса содержит отраженные сигналы, соответствующие отраженным сигналам на первой трассе, и временное окно второй трассы соответствует временному окну на первой трассе, и вычисление отношения разности квадратов записей сигналов в выделенных временными окнами участках записи первой и второй трасс к квадрату записи сигналов во временном окне первой трассы, и выбор допустимого порогового значения разности меньшим, чем это отношение.
- 13. Способ по п.7, отличающийся тем, что длительность временного окна равна длительности не менее двух импульсов отраженных сигналов.
- 14. Способ по п.13, отличающийся тем, что длительность временного окна равна длительности не менее пяти импульсов отраженных сигналов.
- 15. Способ по п.1, отличающийся тем, что указанное применение корректирующей функции к третьей группе отраженных сигналов во второй совокупности сейсмических трасс включает свертку корректирующей функции с третьей группой отраженных сигналов на второй совокупности трасс.
- 16. Способ по п.1, отличающийся тем, что указанная первая совокупность данных включает трассу от сейсмоприемника.
- 17. Способ по п. 1, отличающийся тем, что указанная вторая совокупность данных включает просуммированные сейсмические трассы от сейсмоприемников.
- 18. Способ по п.1, отличающийся тем, что указанные первая и вторая совокупности данных включают просуммированные трассы от расстановки скважинных приемников.
- 19. Способ по п.1, отличающийся тем, что первая и вторая совокупности данных включают данные после предварительной обработки.
- 20. Способ по п.19, отличающийся тем, что предварительно обработанные данные включают сейсмограммы ОСТ.
- 21. Способ по п.19, отличающийся тем, что указанные предварительно обработанные данные включают сейсмограммы ОПВ.
- 22. Способ по п.19, отличающийся тем, что указанные предварительно обработанные данные включают мигрированные данные.
- 23. Способ по п.1, отличающийся тем, что указанные первая и вторая совокупности данных подвергаются эквивалентной предварительной обработке.
- 24. Способ по п. 1, отличающийся тем, что первая совокупность данных подвергается той же процедуре восстановления формы сигнала, что и вторая совокупность.
- 25. Способ по п.1, отличающийся тем, что первая совокупность данных подвергается той же деконволюционной обработке, что и вторая совокупность трасс.
- 26. Способ по п. 1, отличающийся тем, что первая совокупность данных подвергается обработке с использованием тех же процедур подавления помех, что и вторая совокупность данных.
- 27. Способ по п. 1, отличающийся тем, что и первая, и вторая совокупности данных подвергаются одинаковой обработке с целью подавления многократных волн.
- 28. Способ по п.1, отличающийся тем, что одинаковый ΌΜΟ оператор используется при обработке первой и второй совокупностей данных.
- 29. Способ по п.1, отличающийся тем, что миграция первой и второй совокупностей данных производится с эквивалентной скоростной моделью среды.
- 30. Способ по п. 1, отличающийся тем, что миграция первой и второй совокупностей данных производится с использованием одинакового оператора миграции.
- 31. Способ по п.1, отличающийся тем, что фильтрация первой и второй совокупностей данных производится с использованием одинаковых фильтров.
- 32. Способ получения разностной записи для данных, по меньшей мере, двух сейсмических съемок, включающий подготовку данных, по меньшей мере, двух сейсмических съемок на одной и той же географической площади; и осуществление процесса деконволюции данных, по меньшей мере, двух сейсмических съемок, предусматривающий расчет деконволюционного оператора на основе использования, по меньшей мере, двух совокупностей сейсмических данных, из которых не менее двух совокупностей сейсмических данных были зарегистрированы в разное время;использование этого оператора деконволюции при обработке, по меньшей мере, двух совокупностей сейсмических данных; и проведение дальнейшей обработки данных разновременных съемок для получения разностной записи.
- 33. Способ по п.32, отличающийся тем, что не менее двух совокупностей данных включают не менее трех массивов сейсмических данных, причем каждый из этих массивов представляет записи, полученные при съемках, проведенных в разное время, при этом расчет оператора деконволюции производится на основе использования данных каждого массива.
- 34. Способ по п.32, отличающийся тем, что не менее двух совокупностей данных включает не менее трех массивов сейсмических данных, при этом каждый из этих массивов представляет записи, полученные при съемках, проведенных в разное время, расчет включает определение первого оператора деконволюции на основе первой и второй совокупностей данных, и первый деконволюционный оператор используется для обработки данных первой и второй совокупностей, дальнейшая обработка включает расчет второго оператора деконволюции на основе первой и третьей совокупностей данных, причем второй оператор деконволюции применяется для обработки первой и третьей совокупностей данных.
- 35. Способ по п.34, отличающийся тем, что первая совокупность данных представляет записи, полученные при съемке, проведенной до съемки, представленной вторым и третьим массивами данных.
- 36. Способ по п.34, отличающийся тем, что расчет первого оператора деконволюции включает осреднение спектра мощности для первой совокупности трасс данных первой съемки, обратное преобразование усредненных значений с получением функции, отображающей усредненный спектр мощности первой совокупности трасс во временной области, запоминание этой временной функции, осреднение спектра мощности для второй совокупности трасс данных второй съемки, обратное преобразование усредненных значений с получением функции, отображающей усредненный спектр мощности второй совокупности трасс во временной области, осреднение временных функций, представляющих средние значения спектров мощности первой и второй совокупностей трасс и расчет оператора деконволюции для обработки данных первой и второй съемок с использованием усредненных значений временных функций, отображающих средний спектр мощности первой и второй совокупностей.
- 37. Способ по п.32, отличающийся тем, что расчет оператора деконволюции включает:осреднение спектра мощности для первой совокупности трасс данных первой съемки;осреднение спектра мощности для второй совокупности трасс данных второй съемки;осреднение средних значений спектров мощности первой и второй совокупностей трасс и расчет оператора деконволюции для обработки данных первой и второй съемок на основе осредненных значений спектров мощности первой и второй совокупностей сейсмических трасс.
- 38. Способ формирования сигнала, характеризующего географическую площадь по результатам при проведении разновременных сейсмических съемок, включающий:подготовку данных, по меньшей мере, двух сейсмических съемок этой географической площади, первая из которых проведена в первый период времени, а вторая - во второй период времени; и осуществление процесса обработки этих сейсмических данных, предусматривающий выполнение графа обработки входных данных, полученных при первой съемке;выполнение этого же графа обработки входных данных, полученных при второй съемке;причем граф обработки определяется системой параметров сейсмических сигналов;выбор не менее одной системы параметров путем их подбора на основе данных обеих съемок, и использование не менее одной системы параметров, по меньшей мере, в одном графе при обработке данных, полученных как при первой, так и при второй съемках.
- 39. Способ по п.38, отличающийся тем, что, по меньшей мере, одна из систем параметров включает оператор деконволюции.
- 40. Способ по п.39, отличающийся тем, что он дополнительно включает расчет оператора деконволюции на основе использования, по меньшей мере, двух совокупностей сейсмических данных, причем не менее двух совокупностей сейсмических данных были зарегистрированы в разное время;использование этого оператора деконволюции при обработке, по крайней мере, двух совокупностей сейсмических данных; и проведение дальнейшей обработки данных разновременных съемок для получения разностной записи.
- 41. Способ по п.40, отличающийся тем, что не менее двух совокупностей данных включают не менее трех массивов сейсмических данных, причем каждый из них представлен записями, проведенными в разное время, при этом конструкция оператора деконволюции зависит от данных, полученных при каждой съемке.
- 42. Способ по п.40, отличающийся тем, что не менее двух совокупностей данных включают не менее трех массивов сейсмических данных, при этом:каждый из этих трех массивов представляет записи, полученные при проведении съемок в разное время, расчет первого оператора деконволюции основан на использовании первого и второго массива данных, и первый оператор деконволюции применяют при обработке первой и второй совокупностей данных, дальнейшая обработка включает расчет второго оператора деконволюции на основе использования данных, полученных при первой и третьей съемках, второй оператор деконволюции используют при обработке данных первой и третьей съемок.
- 43. Способ по п.42, отличающийся тем, что первая совокупность данных представлена записями, полученными при съемке, проведенной до съемки, представленной вторым и третьим массивами данных.
- 44. Способ по п.42, отличающийся тем, что расчет первого оператора деконволюции включает осреднение спектра мощности для первой совокупности трасс данных первой съемки;обратное преобразование усредненных значений с получением функции, отображающей усредненный спектр мощности первой совокупности трасс во временной области;запоминание этой временной функции;усреднение спектра мощности для второй совокупности трасс данных второй съемки, обратное преобразование усредненных значений с получением функции, отображающей усредненный спектр мощности второй совокупности трасс во временной области;осреднение временных функций, представляющих средние значения спектров мощности первой и второй совокупностей трасс; и расчет оператора деконволюции для обработки данных первой и второй съемок с использованием усредненного значения временных функций, отображающих средние спектры мощности первой и второй совокупностей.
- 45. Способ по п.40, отличающийся тем, что расчет оператора деконволюции включает осреднение спектра мощности для первой совокупности сейсмических трасс данных первой съемки;осреднение спектра мощности для второй совокупности сейсмических трасс данных второй съемки, осреднение средних значений спектров мощности первой и второй совокупностей сейсмических трасс данных; и расчет оператора деконволюции для обработки данных первой и второй съемок на основе осредненного значения спектров мощности первой и второй совокупностей сейсмических трасс данных.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US71394896A | 1996-09-13 | 1996-09-13 | |
PCT/US1997/016054 WO1998011455A1 (en) | 1996-09-13 | 1997-09-11 | Method for time lapse reservoir monitoring |
Publications (2)
Publication Number | Publication Date |
---|---|
EA199900195A1 EA199900195A1 (ru) | 1999-12-29 |
EA001212B1 true EA001212B1 (ru) | 2000-12-25 |
Family
ID=24868193
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EA199900195A EA001212B1 (ru) | 1996-09-13 | 1997-09-11 | Способ сейсмического мониторинга нефтегазовых месторождений с использованием разновременных съемок |
Country Status (9)
Country | Link |
---|---|
US (1) | US6438069B1 (ru) |
EP (1) | EP0937267A4 (ru) |
CN (1) | CN1230260A (ru) |
AU (1) | AU735749B2 (ru) |
EA (1) | EA001212B1 (ru) |
NO (1) | NO991147L (ru) |
OA (1) | OA10992A (ru) |
SE (1) | SE515924C2 (ru) |
WO (1) | WO1998011455A1 (ru) |
Families Citing this family (57)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5946271A (en) * | 1997-03-21 | 1999-08-31 | Western Atlas International, Inc. | Calibration system for use in time lapse tomography |
GB9813760D0 (en) * | 1998-06-25 | 1998-08-26 | Geco Prakla Uk Ltd | Seismic data signal processing method |
WO2000019240A2 (en) * | 1998-09-28 | 2000-04-06 | Pgs Seres As | Amplitude variation as a function of offset attribute and rock property contrast analysis for seismic survey data |
GB2372567B (en) * | 2001-02-22 | 2003-04-09 | Schlumberger Holdings | Estimating subsurface subsidence and compaction |
WO2002075363A1 (en) * | 2001-03-15 | 2002-09-26 | Shell Internationale Research Maatschappij B.V. | Time-lapse seismic surveying at sea |
US6906982B2 (en) | 2002-04-16 | 2005-06-14 | Shell Oil Company | Marine time-lapse seismic surveying |
GB2390904B (en) | 2002-07-16 | 2004-12-15 | Univ Southampton | Electromagnetic surveying for hydrocarbon reservoirs |
US6739165B1 (en) | 2003-02-05 | 2004-05-25 | Kjt Enterprises, Inc. | Combined surface and wellbore electromagnetic measurement system and method for determining formation fluid properties |
US7027354B2 (en) * | 2003-09-22 | 2006-04-11 | 4Th Wave Imaging Corp. | Method of obtaining pore pressure and fluid saturation changes in subterranean reservoirs by forward modeling |
WO2005040858A1 (en) * | 2003-10-24 | 2005-05-06 | Shell Internationale Research Maatschappij B.V. | Time-lapse seismic survey of a reservoir region |
GB2409900B (en) | 2004-01-09 | 2006-05-24 | Statoil Asa | Processing seismic data representing a physical system |
GB2414299B (en) | 2004-05-21 | 2006-08-09 | Westerngeco Ltd | Interpolation and extrapolation method for seismic recordings |
WO2006003137A1 (en) * | 2004-06-30 | 2006-01-12 | Shell Internationale Research Maatschappij B.V. | Determining in-situ the relation between seismic velocity and state of stress in an underground formation |
GB2420408B (en) | 2004-11-19 | 2008-03-12 | Geophysique Cie Gle | Method for processing at least two sets of seismic data |
US7523003B2 (en) * | 2006-07-12 | 2009-04-21 | Westerngeco L.L.C. | Time lapse marine seismic surveying |
US7492665B2 (en) | 2006-07-24 | 2009-02-17 | Westerngeco L.L.C. | Seismic data acquisition and source-side derivatives generation and application |
GB2442749B (en) | 2006-10-12 | 2010-05-19 | Electromagnetic Geoservices As | Positioning system |
KR100861084B1 (ko) * | 2006-12-08 | 2008-09-30 | 한국지질자원연구원 | 물리탐사 자료의 4차원 역산 방법 및 이를 이용한지하구조의 4차원 영상화 방법. |
US7773455B2 (en) * | 2006-12-13 | 2010-08-10 | Westerngeco L.L.C. | Time-lapse seismic acquisition |
US7636275B2 (en) * | 2007-02-06 | 2009-12-22 | Conocophillips Company | Direct time lapse inversion of seismic data |
WO2008140655A1 (en) * | 2007-05-09 | 2008-11-20 | Exxonmobil Upstream Research Company | Inversion of 4d seismic data |
CA2717095C (en) * | 2008-03-07 | 2017-03-28 | Shell Internationale Research Maatschappij B.V. | Method of marine time-lapse seismic surveying |
GB2459887B (en) * | 2008-05-09 | 2012-04-25 | Total Sa | Combining seismic data sets with overlapping bandwidths |
AU2009333601B2 (en) * | 2008-12-17 | 2014-08-21 | Exxonmobil Upstream Research Company | Method for imaging of targeted reflectors |
WO2010077569A1 (en) * | 2008-12-17 | 2010-07-08 | Exxonmobil Upstream Research Company | System and method for reconstruction of time-lapse data |
US8724429B2 (en) | 2008-12-17 | 2014-05-13 | Exxonmobil Upstream Research Company | System and method for performing time-lapse monitor surverying using sparse monitor data |
AU2014253465B2 (en) * | 2008-12-17 | 2015-07-30 | Exxonmobil Upstream Research Company | Method for imaging of targeted reflectors |
US8451683B2 (en) * | 2009-04-03 | 2013-05-28 | Exxonmobil Upstream Research Company | Method for determining the fluid/pressure distribution of hydrocarbon reservoirs from 4D seismic data |
US8332154B2 (en) | 2009-06-02 | 2012-12-11 | Exxonmobil Upstream Research Company | Estimating reservoir properties from 4D seismic data |
GB2474022B (en) * | 2009-09-30 | 2012-03-07 | Statoilhydro Asa | Improved estimation of time shift based on multi-vintage seismic data |
US8693283B2 (en) * | 2009-09-30 | 2014-04-08 | Statoil Petroleum As | Estimation of time shift based on multi-vintage seismic data |
BR112012017275A2 (pt) | 2010-02-12 | 2016-04-19 | Exxonmobil Upstream Res Co | método e sistema para divisão de modelos de simulação paralelos |
US9367564B2 (en) | 2010-03-12 | 2016-06-14 | Exxonmobil Upstream Research Company | Dynamic grouping of domain objects via smart groups |
CN102096101B (zh) * | 2010-11-24 | 2014-11-12 | 中国石油天然气集团公司 | 提高地震数据处理精度的方法及装置 |
WO2012102784A1 (en) | 2011-01-26 | 2012-08-02 | Exxonmobil Upstream Research Company | Method of reservoir compartment analysis using topological structure in 3d earth model |
CN103534614B (zh) | 2011-03-23 | 2017-08-04 | 离子地球物理学公司 | 用于分析时移地球物理调查中的数据的方法和设备 |
WO2012128784A1 (en) * | 2011-03-23 | 2012-09-27 | Ion Geophysical Corporation | Method and apparatus for analyzing data in subsequent geophysical surveys |
US9772415B2 (en) | 2011-08-05 | 2017-09-26 | Saudi Arabian Oil Company | Correcting time lapse seismic data for overburden and recording effects |
US10139505B2 (en) | 2011-08-09 | 2018-11-27 | Pgs Geophysical As | Digital sensor streamers and applications thereof |
US8717845B2 (en) | 2011-08-24 | 2014-05-06 | Pgs Geophysical As | Quality-based steering methods and systems for 4D geophysical surveys |
US9081118B2 (en) * | 2011-11-21 | 2015-07-14 | Cggveritas Services Sa | Device and method for computing depth velocity variations |
US9354338B1 (en) | 2012-02-22 | 2016-05-31 | Westerngeco L.L.C. | Generating synthetic seismic traces |
US9229122B2 (en) * | 2013-03-04 | 2016-01-05 | Cggveritas Services Sa | Image-domain 4D-binning method and system |
US10310122B2 (en) * | 2013-03-04 | 2019-06-04 | Cgg Services Sas | Increasing similarity between seismic datasets |
US10338250B2 (en) | 2013-03-14 | 2019-07-02 | Pgs Geophysical As | Method of removing incoherent noise |
US10386516B2 (en) | 2013-05-15 | 2019-08-20 | Conocophillips Company | Time-lapse 4D scattering for imaging hydraulically induced fractures |
EP3008281A2 (en) | 2013-06-10 | 2016-04-20 | Exxonmobil Upstream Research Company | Interactively planning a well site |
US9651695B2 (en) | 2013-09-19 | 2017-05-16 | Pgs Geophysical As | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields |
US9864098B2 (en) | 2013-09-30 | 2018-01-09 | Exxonmobil Upstream Research Company | Method and system of interactive drill center and well planning evaluation and optimization |
CN104749620A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 一种方位角抽取道集一致性处理方法 |
CN105277984B (zh) * | 2014-07-04 | 2018-02-02 | 中国石油化工股份有限公司 | 一种时移地震互约束频率一致性处理方法 |
US9695684B2 (en) | 2014-10-23 | 2017-07-04 | Cgg Services Sas | System and method for predicting the front arrival time in reservoir seismic monitoring |
CN105738949B (zh) * | 2016-03-01 | 2017-11-17 | 中国海洋石油总公司 | 一种用于时移地震的九面元一致性并行处理方法 |
CN106054242B (zh) * | 2016-05-04 | 2018-11-13 | 中国地质大学(北京) | 三维各向异性衰减介质波场模拟方法 |
CN107678068A (zh) * | 2016-08-01 | 2018-02-09 | 中国石油天然气股份有限公司 | 一种基于重力数据向下延拓的成像方法及装置 |
US10739484B2 (en) * | 2017-03-10 | 2020-08-11 | Exxonmobil Upstream Research Company | Curvelet 4D: 4D denoise in curvelet domain |
CN114427435A (zh) * | 2020-09-22 | 2022-05-03 | 中国石油化工股份有限公司 | 一种三维油藏模型更新方法、装置、电子设备及存储介质 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3921126A (en) * | 1974-08-21 | 1975-11-18 | Continental Oil Co | Method of seismically delineating the permeable extent of oil-containing formations |
US4715021A (en) | 1984-08-08 | 1987-12-22 | Exxon Production Research Co. | Method for filtering and combining seismic data having different spectral characteristics |
US4688198A (en) * | 1984-12-24 | 1987-08-18 | Schlumberger Technology Corporation | Entropy guided deconvolution of seismic signals |
US4969130A (en) * | 1989-09-29 | 1990-11-06 | Scientific Software Intercomp, Inc. | System for monitoring the changes in fluid content of a petroleum reservoir |
US5010976A (en) * | 1989-10-04 | 1991-04-30 | Atlantic Richfield Company | Characterization of the full elastic effect of the near surface on seismic waves |
US5309360A (en) * | 1991-05-23 | 1994-05-03 | Halliburton Geophysical Services, Inc. | Method for attenuating undesirable data, such as multiples, using constrained cross-equalization |
US5173879A (en) | 1992-06-25 | 1992-12-22 | Shell Oil Company | Surface-consistent minimum-phase deconvolution |
US5400299A (en) | 1993-08-20 | 1995-03-21 | Exxon Production Research Company | Seismic vibrator signature deconvolution |
US5586082A (en) | 1995-03-02 | 1996-12-17 | The Trustees Of Columbia University In The City Of New York | Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging |
-
1997
- 1997-09-11 EP EP97941023A patent/EP0937267A4/en not_active Withdrawn
- 1997-09-11 AU AU42670/97A patent/AU735749B2/en not_active Ceased
- 1997-09-11 EA EA199900195A patent/EA001212B1/ru not_active IP Right Cessation
- 1997-09-11 CN CN97197846.8A patent/CN1230260A/zh active Pending
- 1997-09-11 WO PCT/US1997/016054 patent/WO1998011455A1/en not_active Application Discontinuation
- 1997-09-11 US US09/242,232 patent/US6438069B1/en not_active Expired - Fee Related
-
1999
- 1999-03-10 NO NO991147A patent/NO991147L/no not_active Application Discontinuation
- 1999-03-10 SE SE9900849A patent/SE515924C2/sv not_active IP Right Cessation
- 1999-03-10 OA OA9900053A patent/OA10992A/en unknown
Also Published As
Publication number | Publication date |
---|---|
WO1998011455A1 (en) | 1998-03-19 |
EP0937267A4 (en) | 2001-09-05 |
NO991147D0 (no) | 1999-03-10 |
SE9900849L (sv) | 1999-03-10 |
EP0937267A1 (en) | 1999-08-25 |
SE515924C2 (sv) | 2001-10-29 |
AU4267097A (en) | 1998-04-02 |
US6438069B1 (en) | 2002-08-20 |
OA10992A (en) | 2001-11-07 |
CN1230260A (zh) | 1999-09-29 |
AU735749B2 (en) | 2001-07-12 |
SE9900849D0 (sv) | 1999-03-10 |
EA199900195A1 (ru) | 1999-12-29 |
NO991147L (no) | 1999-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EA001212B1 (ru) | Способ сейсмического мониторинга нефтегазовых месторождений с использованием разновременных съемок | |
EP2103959B1 (en) | Method for interpolating seismic data by anti-alias, anti-leakage fourier transform | |
Velis | Stochastic sparse-spike deconvolution | |
Lavergne et al. | Inversion of seismograms and pseudo velocity logs | |
EP2108980B1 (en) | Method for deghosting marine seismic streamer data with irregular receiver positions | |
US5661697A (en) | Method and apparatus for detection of sand formations in amplitude-versus-offset seismic surveys | |
WO1998011455A9 (en) | Method for time lapse reservoir monitoring | |
KR102021276B1 (ko) | 진폭 보존을 갖는 fwi 모델 도메인 각도 스택들 | |
US9234977B2 (en) | Processing collected survey data | |
EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
EA001509B1 (ru) | Способ обработки сигналов сейсмических данных | |
US6430508B1 (en) | Transfer function method of seismic signal processing and exploration | |
MacBeth et al. | Methods of measurement for 4D seismic post‐stack time shifts | |
Adamczyk et al. | Full-waveform inversion of conventional Vibroseis data recorded along a regional profile from southeast Poland | |
Kolkman-Quinn | Time-lapse VSP monitoring of CO2 sequestration at the CaMI Field Research Station | |
Zhang et al. | Interval Q inversion based on zero-offset VSP data and applications | |
Landrø et al. | Target-oriented AVO inversion of data from Valhall and Hod fields | |
Liu et al. | Two‐dimensional SH‐wave and acoustic P‐wave full waveform inversion: A Midland Basin case study | |
US20240142648A1 (en) | Method and system for estimating converted-wave statics | |
Yu et al. | Joint OBN and 3D DAS VSP Data Acquisition and Processing in Offshore Abu Dhabi | |
US20240142650A1 (en) | Efficient impulse removal in fourier anti-leakage seismic data interpolation | |
Gadallah et al. | Seismic Interpretation | |
WO2023027593A1 (en) | Method of simulating seismic data | |
Sayed et al. | Specialized processing of a massive distributed acoustic sensing (DAS) 3D VSP in a producing well | |
Bustamante et al. | Estimation of the displacements among distant events based on parallel tracking of events in seismic traces under uncertainty |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s) |
Designated state(s): AM AZ BY KZ KG MD TJ TM RU |