EA007881B1 - Оценивание времени вступления сейсмической волны - Google Patents

Оценивание времени вступления сейсмической волны Download PDF

Info

Publication number
EA007881B1
EA007881B1 EA200501270A EA200501270A EA007881B1 EA 007881 B1 EA007881 B1 EA 007881B1 EA 200501270 A EA200501270 A EA 200501270A EA 200501270 A EA200501270 A EA 200501270A EA 007881 B1 EA007881 B1 EA 007881B1
Authority
EA
Eurasian Patent Office
Prior art keywords
wave
time
estimated
entry
calculated
Prior art date
Application number
EA200501270A
Other languages
English (en)
Other versions
EA200501270A1 (ru
Inventor
Роберт Хьюз Джоунз
Original Assignee
Шлюмбергер Текнолоджи Б.В.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Шлюмбергер Текнолоджи Б.В. filed Critical Шлюмбергер Текнолоджи Б.В.
Publication of EA200501270A1 publication Critical patent/EA200501270A1/ru
Publication of EA007881B1 publication Critical patent/EA007881B1/ru

Links

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. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G01V1/305—Travel times
    • G01V1/01—
    • 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. analysis, for interpretation, for correction
    • G01V1/284—Application of the shear wave component and/or several components of the seismic signal

Abstract

Средство обработки данных предназначено для вычисления расчетного времени вступления сейсмической или микросейсмической Р- или S-волны на приемную станцию на основе отношения скоростей Р- и S-волн, вычисленного расчетного времени начала события и, когда должно быть вычислено расчетное время вступления Р-волны, выбранного времени вступления S-волны на приемную станцию, или, когда должно быть вычислено расчетное время вступления S-волны, выбранного времени вступления Р-волны на приемную станцию.

Description

Область техники
Настоящее изобретение относится к способу вычисления расчетного времени вступления сейсмической или микросейсмической Р- или 8-волны на приемную станцию к средству обработки данных для осуществления способа оценки и к носителю данных, содержащему считываемое компьютером программное средство для осуществления указанного способа, когда он выполняется на компьютере.
Предшествующий уровень техники
Типичное сейсмическое или микросейсмическое событие создает сигналы, соответствующие как Рволне, так и 8-волне. Эти отдельные сигналы проходят через среду с различными скоростями. Для широкого класса пород земной коры отношение скоростей Р-волны и 8-волны хорошо известно и не изменяется в значительной степени. Безусловно, во многих случаях априори невозможно получить точную скорость 8-волны, и поэтому распределение скорости 8-волны для района, представляющего интерес, часто получают в предположении постоянного отношения при наличии распределения скорости Р-волны. Поэтому на практике для распределений двух скоростей сейсмических волн в пределах объема, представляющего интерес, указывается приблизительно постоянное отношение друг к другу, даже когда значения скоростей обоих типов изменяются.
Обычный способ оценивания местоположения и времени (в значениях координат х, у, ζ и начала (10) отсчета времени), в которое сейсмическое или микросейсмическое событие произошло, включает в себя позиционирование нескольких приемных станций, содержащих геофоны или аналогичные элементы, на различных местах в зоне охвата, и измерение времен вступлений создаваемых Р- и 8-волн на каждую из указанных станций. После того как вступление каждой волны на каждую станцию зарегистрировано (то есть выбрано на основе сейсмических данных, зарегистрированных на рассматриваемой станции), делают предположение о местоположении и времени начала сейсмического или микросейсмического события, а известные или предполагаемые скорости Р- и 8-волн используют для вычисления моментов времени, в которые Р- и 8-волны достигают каждой станции, допуская, что это предположение было правильным. Эти вычисленные времена вступлений, основанные на предполагаемом местоположении и времени начала, затем сравнивают с зарегистрированными временами вступлений (которые, как отмечалось выше, основаны на сейсмических данных, действительно зарегистрированных), а местоположение и время начала корректируют с помощью математической обработки, известной как метод наименьших квадратов, так что конечное расчетное положение и время начала события оказываются наилучшим образом согласованными в смысле наименьших квадратов с зарегистрированными временами вступлений волн на станции.
Хотя указанный выше способ обычно считается наилучшим способом определения местоположения события из доступных в настоящее время, надежность и точность определения местоположения данного сейсмического или микросейсмического события, регистрируемого в определенной сети, сильно зависит от точности и надежности выбора времени вступления. В некоторых случаях качество сейсмических данных может вообще затруднить идентификацию отдельного времени вступления. Если выбранное время вступления очень неточное или еще хуже, тип волны неверно идентифицирован (например, когда Р-волна выбрана как 8-волна или обмен волны с Р на 8 выбран как 8-волна), итоговое расчетное местоположение и время начала события, вероятно, будут значительно отличаться от истинных местоположения и времени. Кроме того, обработка местоположения методом наименьших квадратов означает, что влияние одного неправильного измерения времени распространяется и проявляет себя во всех остальных каналах и в итоговом несоответствии положения, что делает невозможной идентификацию неверного снятого значения.
Краткое изложение существа изобретения
Согласно одному аспекту настоящего изобретения предложено средство обработки данных для вычисления расчетного времени вступления сейсмической или микросейсмической Р- или 8-волны на приемную станцию, содержащее процессор данных, выполненный с возможностью
a) вычисления расчетного времени начала сейсмического или микросейсмического события, генерирующего Р- и 8-волны, на основе отношения скоростей Р- и 8-волн и зарегистрированных времен вступлений Р- и 8-волн на другую приемную станцию, а не на ту приемную станцию, для которой должно быть вычислено расчетное время вступления Р- или 8-волны,
b) вычисления расчетного времени вступления Р- или 8-волны на основе отношения скоростей Р- и 8-волн, вычисления расчетного времени начала события и, когда расчетное время вступления Р-волны должно быть вычислено, вычисляется выбранное время вступления 8-волны на приемную станцию, для которой подсчитывается расчетное время вступления Р-волны, или, когда расчетное время вступления 8волны должно быть вычислено, вычисляется выбранное время вступления Р-волны на приемную станцию, для которой подсчитывается расчетное время вступления 8-волны.
Процессор данных может быть выполнен с возможностью вычисления расчетных времен вступлений как для Р- , так и для и 8-волн на приемную станцию.
Предпочтительно, чтобы процессор данных был выполнен с возможностью вычисления большого количества расчетных времен вступлений Р- и/или 8-волн на приемную станцию на основе большого количества расчетных времен начала для микросейсмического события, вычисленных по выбранным
- 1 007881 временам вступлений Р- и 8-волн на большое количество других приемных станций, а не на той приемной станции, для которой должны быть вычислены расчетные времена вступлений. Процессор данных может быть дополнительно выполнен с возможностью отображения выбранных времен вступлений и расчетных времен вступлений во взаимосвязи друг с другом, так что может быть проанализирована кластерная картина времен вступлений. Когда требуется именно этот анализ, процессор данных также может обеспечивать возможность отображения информации, касающейся вычисления любого конкретного расчетного времени вступления, в ответ на выбор пользователем указанного расчетного времени вступления.
Предпочтительно, чтобы процессор данных был выполнен с возможностью вычисления одного или нескольких расчетных времен вступлений для Р- и/или 8-волн на каждую приемную станцию в сети приемных станций, при этом для каждой приемной станции необходимое расчетное время или времена начала вычисляются по выбранным временам вступлений Р- и 8-волн на одну или несколько других станций в указанной сети.
В некоторых случаях для процессора данных может потребоваться, чтобы пользователь выбрал времена вступлений для Р- и 8-волн на различные приемные станции путем исследования сейсмических данных, зарегистрированных на каждой станции, и затем снабдил процессор данных выбранными временами вступлений. Альтернативно, процессор данных может быть выполнен с возможностью получения сейсмических данных от приемных станций и выполнения на основе указанных сейсмических данных выбора времен вступлений для Р- и 8-волн на каждую приемную станцию.
В тех случаях, когда может быть выбрано несколько возможных времен вступлений волны на приемную станцию (либо пользователем, либо самим процессором данных), процессор данных может быть выполнен с возможностью сравнения указанных возможных времен вступлений с любыми оценками, вычисленными для времени вступления указанной волны на указанную станцию, чтобы определить, какое из возможных выбранных времен вступлений наиболее вероятно соответствует истинному времени вступления указанной волны на указанную станцию. Процессор данных может быть дополнительно выполнен с возможностью выбора или изменения одного из указанных возможных времен вступлений, чтобы получить в итоге выбранное время вступления, которое на основе указанного определения представляется наиболее правдоподобно соответствует истинному времени вступления указанной волны. Альтернативно, процессор данных может быть выполнен с возможностью отображения, какое из возможных времен вступлений должно быть выбрано или изменено, чтобы получить в итоге выбранное время вступления, которое на основе указанного определения наиболее правдоподобно соответствует истинному времени вступления.
Согласно другому аспекту настоящего изобретения предложен носитель данных, содержащий считываемое компьютером программное средство, чтобы компьютер мог функционировать в качестве средства обработки данных, согласно изобретению.
Согласно еще одному аспекту настоящего изобретения предложен способ вычисления расчетного времени вступления сейсмической или микросейсмической Р- или 8-волны на приемную станцию, при этом способ содержит следующие этапы:
a) вычисление расчетного времени начала сейсмического или микросейсмического события, генерирующего Р- и 8-волны, на основе отношения скоростей Р- и 8-волн и выбранных времен вступлений Ри 8-волн на другую приемную станцию, а не на ту приемную станцию, для которой расчетное время вступления Р- или 8-волны должно быть вычислено,
b) вычисление расчетного времени вступления Р- или 8-волны на основе отношения скоростей Р- и 8-волн, вычисление расчетного времени начала события и, когда расчетное время вступления Р-волны должно быть вычислено, вычисляется выбранное время вступления 8-волны на приемную станцию, для которой подсчитывается расчетное время вступления Р-волны, или, если расчетное время вступления 8волны должно быть вычислено, вычисляется выбранное время вступления Р-волны на приемную станцию, для которой подсчитывается расчетное время вступления 8-волны.
Следовательно, согласно настоящему изобретению предложено средство для получения одной или нескольких оценок времен вступления Р- и/или 8-волн на одну или несколько приемных станций без необходимости первоначального определения местоположения сейсмического или микросейсмического события. Полученные таким путем расчетные времена вступлений могут быть использованы различным образом. В тех случаях, когда сейсмические данные, зарегистрированные на станции, не могут быть интерпретированы с достаточной надежностью, чтобы осуществить выбор, оценка времени вступления может быть использована вместо значения, выбор которого не выполнен. В тех случаях, когда имеется несколько приемных станций, может быть вычислен разброс расчетных времен вступлений, чтобы снабдить пользователя информацией, достаточной для идентификации возможных неудовлетворительных выбранных значений до попытки определить местоположение события. Такие разбросы оценок также могут быть использованы для анализа распределения и отношения скоростей Р- и 8-волн в части процесса автоматического определения местоположения и/или для оценки неопределенности выбора.
Ниже описан предпочтительный вариант осуществления изобретения.
Согласно описываемому варианту средство обработки данных содержит выполняемое на компью
- 2 007881 тере программное обеспечение для вычисления необходимых расчетных времен начала и вступления (хотя равным образом это может быть осуществлено с помощью подходящего аппаратного обеспечения) путем использования определенных алгоритмов, вывод и применение которых поясняются ниже.
Если предположить, что сеть приемных станций (каждая из которых может быть сейсмическим многокомпонентным приемником, описанным, например, в ПВ-А-2275337) развернута в изотропной среде произвольной сложности с постоянным или почти постоянным отношением скоростей Р- и 8-волн, расстояние (И) микросейсмического события от станции может быть выражено в виде
и/или
Р«(Тр-То)Ур где Уз - скорость 8-волны; Ур - скорость Р-волны; Тр - время вступления Р-волны на станцию; Тз время вступления 8-волны на станцию; и То - время начала микросейсмического события, создающего указанные Р- и 8-волны. Умножение двух уравнений дает
ТвУв-ТоУв = В и/или
ТрУр - ТоУр - В.
Исключив И, расстояние до источника, получим одно уравнение ТзУз - ТоУз - ТрУр - ТоУр.
Группировка всех членов, включающих в себя То, время начала, дает ТоУр - То Уе - ТрУр - ТвУз,
ТоУр/Ув - То=ТрУр/Ув - Те,
То(Ур/Уз -1) - ТрУр/Ув - Тз.
где Ур/Уз представляет собой отношение скорости Р-волны к скорости 8-волны, которое будем обозначать как К.
Подстановка К вместо Ур/У з приводит уравнение к виду
Теперь мы имеем оценку времени начала землетрясения в членах времен (Тр и Тз) вступлений Р- и 8-волн на станцию и отношения К между полями скоростей Р-волны и 8-волны.
В этом уравнении К присутствует дважды, и для его преобразования с тем, чтобы К фигурировало только один раз, можно добавить члены +Тр и -Тр в числитель уравнения
Преобразовав уравнение, получим
Разделение членов числитель приводит к виду
Теперь можно исключить член (К-1) из первой части правой стороны уравнения
После окончательного преобразования получают
То=Тр-(Тв-ТрУ(К-1).
Поэтому, когда времена вступлении Р- и 8-волн на станцию могут быть выбраны с относительно высокой степенью определенности, и при наличии модели среды с постоянным отношением скоростей, на основании данных этой станции можно получить оценку времени начала события, основанную только на выбранных временах вступлений Р- и 8-волн и на предполагаемом отношении полей скоростей Р- и 8волн. Отметим, что для этих уравнений не требуется знание действительного местоположения микро
- 3 007881 сейсмического события или полей скоростей для Р- и 8-волн.
В случае, когда станция, для которой получают выбранные времена вступлений Р- и 8-волн, обозначена как станция 1, можно добавить нижний индекс 1, чтобы показать, что эти члены применяются к станции 1, чтобы получить
ВД “ ВД ~ (Τβι - ВД)/(В. - 1), а) где То1 - время начала, оцениваемое путем использования результатов измерений времени на станции 1.
Точно так же, если времена вступлений были выбраны для другой станции, обозначенной как станция 2, можно добавить нижний индекс 2, чтобы получить
ВД»ВД-(ВД-ВД№1). (2)
Полученное в предположении изотропной среды отношение полей скоростей остается постоянным для обеих станций.
Однако времена начала, вычисленные по данным станций 1 и 2, должны быть идентичными, поскольку одно и то же событие генерирует волны, регистрируемые на обеих станциях. Поэтому в условиях, когда выбор 8-волны для станции 2 невозможен или точность выбора представляется неопределенной, так что желательно получить оценку для времени вступления 8-волны на станцию 2, то она может быть получена путем подстановки То2 вместо То1 и Тзе2 вместо Тз2 в уравнение (2)
ВД = ВД-(ВД-ВДУ(В-1),
ВД(н-1) “Тр^-1) - (ВД - ВД
-вдСк-О-вд+вд,
ВД « ВД(В.-1) - Το,ίΚ-1) +Тр2 “ ВДВ. - ВД - ВДВ - ВД + ВД “ВДВ.-ВДВ.-ВД «<ВД-ВД)К-ВД.
ВД“НСвд -ВД)-ВД, где Тзе2 обозначает оценку для времени вступления 8-волны на станцию 2.
Отметим, что для вычисления расчетного времени вступления 8-волны на станцию 2 требуются только выбранные Р- и 8-волны со станции 1 (чтобы вычислить То1), выбранная Р-волна со станции 2 и предполагаемое постоянное отношение К.
Точно так же, когда желательно получить оценку для времени вступления Р-волны на станцию 2, ее можно получить путем подстановки в уравнение (2) То2 вместо То1 и Тре2 вместо Тр2
ВД=Трй-^-ВДУ(М),
ВДВ-Ц’ВДМ-СВД-ВД) ~ВД(М)-ВД + ВД,
ВД(В-1)+Те, “ ВД(В-1) * ВД
ВДК-ВД+ВД “ВД&вд-ао^-п+т^ув (ВДВ. -ВД +Тв^В.
- ВДВ/В+(Та, - ВДУВ, вд-вд+свдвдув,
- 4 007881 где Тре2 - расчетное время вступления Р-волны на станцию 2.
Отсюда следует, что, когда выбранные Р-волна и 8-волна для станции 2 получены, процесс, описанный выше, может быть обращен, чтобы получить оценки для выбранных Р-волны и 8-волны для станции 1. Действительно, используя указанный выше процесс, по данным каждой станции, для которой получены выбранные Р-волна и 8-волна, можно вычислить расчетное время вступления для Р- и/или 8волн на другую станцию в сети (предполагая, что время вступления Р- или 8-волны, не будучи оцененным, выбрано для указанной станции).
Поэтому для сети из шести станций в случае, если выбранные Р- и 8-волны получены для каждой станции, пять расчетных времен вступлений для каждой волны на каждую станцию также можно вычислить. Если выбор не может быть сделан для времени вступления одной волны на одну из станций (вследствие того, что время вступления волны не может быть идентифицировано по сейсмическим данным), то одно из пяти расчетных времен вступления для указанной волны на указанную станцию (или среднее из них) может быть использовано для последующего определения местоположения микросейсмического события. Точно так же, если все выбранные значения получены, но одно из выбранных значений существенно отличается от истинного времени прихода (опять вследствие трудности идентификации истинного времени вступления по зарегистрированным данным), то оценки, осуществляемые на основании неудовлетворительного выбранного значения, будут существенно отличаться от других оценок, в которые не включено неудовлетворительное выбранное значение, что позволяет идентифицировать и скорректировать неудовлетворительное выбранное значение (или на самом деле, заменить оценкой) до определения местоположения. Здесь имеется отличие от попытки исследовать аномалии канала, создаваемые при выполнении процесса определения местоположения, чтобы идентифицировать неудовлетворительное выбранное значение, поскольку каждая аномалия канала находится под влиянием всех выбранных значений, используемых при определении местоположения, а также отношения распределений скоростей Р- и 8-волн и трех измерений времени.
Расчетные времена вступлений, которые получают, также могут быть использованы для повышения состоятельности выбора времен вступлений. Если предположение постоянства отношения скоростей Р- и 8-волн справедливо, то все расчетные времена вступлений для волны должны образовывать плотный кластер вокруг выбранного времени вступления для волны. Если такая кластерная картина не наблюдается, то имеются три причины, по которым это может случиться:
1) предположение постоянства отношения скоростей Р- и 8-волн не верно для рассматриваемой области,
2) используемое отношение скоростей Р- и 8-волн неверно,
3) одно или несколько выбранных времен вступлений не соответствуют действительности в значительной степени.
Причины (1) и (2) связаны с распределением скоростей. Если причиной является причина (1), то ее можно легко устранить путем изменения скорости Р- или 8-волны, пока не будет достигнута хорошая кластеризация. Если причиной является причина (2), то отношение скоростей Р- и 8-волн можно изменять до установления оптимального отношения, при котором обеспечивается наилучшая кластеризация вокруг группы многочисленных измерений времени от многих событий. Поэтому процесс оценивания времен вступлений, описываемый в настоящей заявке, также может быть применен для оптимизации распределения скорости. Обычно это относится к скорости 8-волны, которая хуже всего разрешается, и поэтому в большинстве случаев она будет скоростью, которую необходимо изменять, чтобы получить оптимальное отношение скорости и структуры. Можно провести различие между причинами (1) и (2), заключающееся в том, что когда причина (2) действует (т.е. когда отношение является постоянным, но неподходящим), то наряду с оцененными выбранными значениями наблюдается хорошая кластеризация, но не вокруг выбранных времен вступлений. Когда действует причина (1), этот эффект будет менее уловимым, а оцененные выбранные значения не будут образовывать кластер той же самой протяженности. На практике это труднее идентифицировать, но в большинстве случаев изменения отношения скоростей Р- и 8-волн в пределах конкретного района, вероятно, будут небольшими, так что в большинстве случаев они будут эффектом второго порядка. После того как местоположения многочисленных событий определены, на основе полученных данных можно исследовать взаимосвязь между местоположениями и отношением скоростей Р- и 8-волн.
В большинстве случаев действие причины (3) не вызовет образования кластерной картины. Чтобы помочь пользователю определить, какие станции дают неудовлетворительную оценку выбранных времен вступлений, математическое обеспечение позволяет пользователю выбрать любую из отображаемых оценок, используя мышь или клавиши управления курсором, и в ответ на это выделяется станция, что позволяет уточнить результат. Таким образом, при выполнении программы на компьютере пользователь определяет станцию с «худшими» выбранными значениями. Кроме того, этот процесс можно повторять до тех пор, пока кластеризация расчетных и измеряемых выбранных значений не даст удовлетворительных результатов.
Программа, выполняемая на компьютере, также может быть рассчитана на использование оцененных выбранных значений, как часть процесса автоматического определения местоположения. Как указа
- 5 007881 но выше, «удовлетворительное» выбранное значение может быть идентифицировано, исходя из того, что оно имеет вокруг себя кластер оцененных выбранных значений. Во время процесса автоматического измерения времени в компьютере обычно сначала производится отбор нескольких возможных выбранных значений для времени вступления волны на станцию. При использовании указанного выше процесса формирования оценок для каждого выбранного значения и анализа кластеризации оценок вокруг каждого возможного автоматически выбранного значения может быть разработана программа для исключения возможных автоматически выбранных значений, которые с меньшей вероятностью соответствуют истинному времени вступления. В программе также можно использовать степень кластеризации для принятия решения, если волна выделена достаточно хорошо, о целесообразности сохранения и использования в последующем процессе определения местоположения.
Когда оценивание неопределенности или точности выбранного значения является трудным и субъективным процессом, то можно также использовать расчетные времена вступлений, чтобы получить показатель неопределенности выбранного значения. Каждое расчетное время вступления содержит два других времени вступления. Например, время вступления 8-волны на станцию η включает в себя следующие выбранные значения:
Тр! Тз1 Трп
Т₽1
Те,
Тз4
т%
где Τρι обозначает выбранное время вступления для Р-волны на станцию 1; Τδι обозначает выбранное время вступления для 8-волны на станцию 1; Трп обозначает выбранное время вступления Р-волны на станцию η и т.д. Из приведенного описания должно быть понятно, что Трп является постоянной величиной во всех этих оценках. Поэтому разброс в оценках Тзп будет функцией неопределенности выбора при выборе Р- и 8- волн для каждого канала.

Claims (12)

  1. ФОРМУЛА ИЗОБРЕТЕНИЯ
    1. Средство обработки данных для вычисления расчетного времени вступления сейсмической или микросейсмической Р- или 8-волны на приемную станцию, содержащее процессор данных, выполненный с возможностью вычисления расчетного времени начала сейсмического или микросейсмического события, генерирующего Р- и 8-волны, на основе отношения скоростей Р- и 8-волн и выбранных времен вступлений Р- и 8-волн на другую приемную станцию, а не на ту приемную станцию, для которой должно быть вычислено расчетное время вступления Р- или 8-волны, вычисления расчетного времени вступления Р- или 8-волны на основе отношения скоростей Р- и 8волн, а также расчетного времени начала события и, когда расчетное время вступления Р-волны должно быть вычислено, вычисляется выбранное время вступления 8-волны на приемную станцию, для которой подсчитывается расчетное время вступления Р-волны, или, когда расчетное время вступления 8-волны должно быть вычислено, вычисляется выбранное время вступления Р-волны на приемную станцию, для которой подсчитывается расчетное время вступления 8-волны.
  2. 2. Средство обработки данных по п.1, отличающееся тем, что процессор данных выполнен с возможностью вычисления расчетных времен вступлений как для Р-, так и для 8-волн на приемную станцию.
  3. 3. Средство обработки данных по любому из пп.1 или 2, отличающееся тем, что процессор данных выполнен с возможностью вычисления большого количества расчетных времен вступлений Р- и/или 8волн на приемную станцию на основе большого количества расчетных времен начала для микросейсмического события, вычисленных по выбранным временам вступлений Р- и 8-волн на большом количестве других приемных станций, а не на той приемной станции, расчетные времена вступлений на которую должны быть вычислены.
  4. 4. Средство обработки данных по п.3, отличающееся тем, что процессор данных дополнительно выполнен с возможностью отображения выбранных времен вступлений и расчетных времен вступлений во взаимосвязи друг с другом, так что может быть проанализирована кластерная картина времен вступлений.
  5. 5. Средство обработки данных по п.4, отличающееся тем, что процессор данных выполнен с возможностью отображения информации, касающейся вычисления любого конкретного расчетного времени вступления, в ответ на выбор пользователем расчетного времени вступления.
  6. 6. Средство обработки данных по любому из предшествующих пп.1-5, отличающееся тем, что процессор данных выполнен с возможностью вычисления одного или нескольких расчетных времен вступлений для Р- и/или 8-волн на каждую приемную станцию в сети приемных станций, при этом для каждой
    - 6 007881 приемной станции необходимое расчетное время или времена начала вычисляют по выбранным временам вступлений Р- и δ-волн на одну или несколько из других станций в указанной сети.
  7. 7. Средство обработки данных по любому из предшествующих пп.1-6, отличающееся тем, что процессор данных выполнен с возможностью получения сейсмических данных от приемных станций и выбора времен вступлений для Р- и δ-волн на каждую станцию на основе указанных сейсмических данных.
  8. 8. Средство обработки данных по любому из предшествующих пп.1-7, отличающееся тем, что, когда могут быть выбраны несколько возможных времен вступлений для волны на приемную станцию, процессор данных выполнен с возможностью сравнения указанных возможных времен вступлений с любыми оценками, вычисленными для времени вступления указанной волны на указанную станцию, чтобы определить, какое из возможных выбранных времен вступлений наиболее вероятно соответствует истинному времени вступления указанной волны на указанную станцию.
  9. 9. Средство обработки данных по п.8, отличающееся тем, что процессор данных дополнительно выполнен с возможностью выбора или изменения одного из возможных времен вступлений, чтобы получить в итоге выбранное время вступления, которое на основе указанного определения наиболее соответствует истинному времени вступления указанной волны.
  10. 10. Средство для обработки данных по п.8, отличающееся тем, что процессор данных выполнен с возможностью отображения, какое из возможных времен вступлений должно быть выбрано или изменено, чтобы получить в итоге выбранное время вступления, которое на основе указанного определения наиболее соответствует истинному времени вступления.
  11. 11. Носитель данных, содержащий считываемое компьютером программное средство для адаптации компьютера к функционированию в качестве средства обработки данных по любому из пп.1-10.
  12. 12. Способ вычисления расчетного времени вступления сейсмической или микросейсмической Рили δ-волны на приемную станцию, заключающийся в том, что вычисляют расчетное время начала сейсмического или микросейсмического события, генерирующего Р- и δ-волны, на основе отношения скоростей Р- и δ-волн и выбранных времен вступлений Р- и δволн на другую приемную станцию, но не на ту приемную станцию, для которой расчетное время вступления Р- или δ-волны должно быть вычислено, вычисляют расчетное время вступления Р- или δ-волны на основе отношения скоростей Р- и δволн, расчетное время начала события и, когда расчетное время вступления Р-волны должно быть вычислено, вычисляют выбранное время вступления δ-волны на приемную станцию, для которой расчетное время вступления Р-волны вычисляется, или, когда расчетное время вступления δ-волны должно быть вычислено, вычисляют выбранное время вступления Р-волны на приемную станцию, для которой расчетное время вступления δ-волны должно быть вычислено.
EA200501270A 2003-02-08 2004-01-14 Оценивание времени вступления сейсмической волны EA007881B1 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GB0302898A GB2398124B (en) 2003-02-08 2003-02-08 Estimating the time of arrival of a seismic wave
PCT/GB2004/000126 WO2004070424A1 (en) 2003-02-08 2004-01-14 Estimating the time of arrival of a seismic wave

Publications (2)

Publication Number Publication Date
EA200501270A1 EA200501270A1 (ru) 2006-06-30
EA007881B1 true EA007881B1 (ru) 2007-02-27

Family

ID=9952669

Family Applications (1)

Application Number Title Priority Date Filing Date
EA200501270A EA007881B1 (ru) 2003-02-08 2004-01-14 Оценивание времени вступления сейсмической волны

Country Status (7)

Country Link
US (1) US7457195B2 (ru)
EP (1) EP1590689A1 (ru)
CA (1) CA2515345C (ru)
EA (1) EA007881B1 (ru)
GB (1) GB2398124B (ru)
NO (1) NO20053954L (ru)
WO (1) WO2004070424A1 (ru)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104407381A (zh) * 2014-12-17 2015-03-11 中国石油大学(华东) 从井中弹性波速的径向变化获取地层岩石脆裂性质的方法

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7391675B2 (en) 2004-09-17 2008-06-24 Schlumberger Technology Corporation Microseismic event detection and location by continuous map migration
US20070223308A1 (en) * 2006-03-21 2007-09-27 Frivik Svein A Methods of range selection for positioning marine seismic equipment
US7676326B2 (en) * 2006-06-09 2010-03-09 Spectraseis Ag VH Reservoir Mapping
GB2439571B (en) 2006-06-28 2008-11-12 Schlumberger Holdings Method for updating a model of the earth using microseismic measurements
US7590491B2 (en) 2006-06-30 2009-09-15 Spectraseis Ag Signal integration measure for seismic data
US8902707B2 (en) * 2007-04-09 2014-12-02 Baker Hughes Incorporated Analysis of uncertainty of hypocenter location using the combination of a VSP and a subsurface array
GB2463591B (en) * 2007-05-17 2012-04-11 Spectraseis Ag Seismic attributes for reservoir localization
US8908473B2 (en) 2008-12-23 2014-12-09 Schlumberger Technology Corporation Method of subsurface imaging using microseismic data
US8599644B2 (en) * 2009-02-04 2013-12-03 Schlumberger Technology Corporation Velocity models for a single well and for a set of wells
US8670288B2 (en) * 2009-02-04 2014-03-11 Schlumberger Technology Corporation Velocity model for well time-depth conversion
US8639443B2 (en) * 2009-04-09 2014-01-28 Schlumberger Technology Corporation Microseismic event monitoring technical field
US20110093203A1 (en) * 2009-10-21 2011-04-21 Chevron U.S.A. Inc. System and method for clustering arrivals of seismic energy to enhance subsurface imaging
US9410421B2 (en) 2009-12-21 2016-08-09 Schlumberger Technology Corporation System and method for microseismic analysis
US9945970B1 (en) * 2011-08-29 2018-04-17 Seismic Innovations Method and apparatus for modeling microseismic event location estimate accuracy
US11774616B2 (en) 2011-08-29 2023-10-03 Seismic Innovations Method and system for microseismic event location error analysis and display
US9835017B2 (en) 2012-09-24 2017-12-05 Schlumberger Technology Corporation Seismic monitoring system and method
US20140142854A1 (en) * 2012-11-16 2014-05-22 Conocophillips Company Method for locating a microseismic event
US20140223010A1 (en) * 2013-02-01 2014-08-07 David Alan Hayner Data Compression and Encryption in Sensor Networks
US20150195693A1 (en) * 2014-01-04 2015-07-09 Ramin Hooriani Earthquake early warning system utilizing a multitude of smart phones
CN104375180B (zh) * 2014-11-19 2017-04-05 中国石油天然气股份有限公司 一种地震数据处理方法、装置及系统
US9939541B2 (en) * 2015-01-09 2018-04-10 Chevron U.S.A. Inc. Layered linear inversion techniques for locating microseismic activity
US10338246B1 (en) 2015-08-31 2019-07-02 Seismic Innovations Method and system for microseismic event wavefront estimation
CN105487114B (zh) * 2015-12-08 2018-02-02 中南大学 一种微震信号p波初至点综合拾取方法
CN105527650B (zh) * 2016-02-17 2017-10-17 中国科学院武汉岩土力学研究所 一种工程尺度下微震信号及p波初至自动识别算法
CN109521221B (zh) * 2018-12-10 2020-09-29 东北大学 一种钻爆法施工硬岩隧道微震波波速实时获取方法
CN111596347B (zh) * 2019-02-21 2023-08-22 中国石油天然气集团有限公司 一种快速求取表层纵横波速度比的方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6067275A (en) * 1997-12-30 2000-05-23 Schlumberger Technology Corporation Method of analyzing pre-stack seismic data
WO2001031364A1 (fr) * 1999-10-28 2001-05-03 Compagnie Generale De Geophysique Procede de prospection sismique mettant en oeuvre un traitement sur les ondes converties
WO2002059647A1 (en) * 2001-01-25 2002-08-01 Westerngeco Seismic Holdings Limited A method of processing marine seismic data and a method of seismic surveying

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4596005A (en) * 1983-04-20 1986-06-17 Chevron Research Company Method of seismic collection utilizing multicomponent processing receivers and processing resultant conventional and converted P- or S-wave data
US4611311A (en) * 1983-04-20 1986-09-09 Chevron Research Company Method of seismic processing involving converted P- or S-wave data
US4686657A (en) * 1984-09-24 1987-08-11 Mobil Oil Corporation Matching P wave and shear wave sections in seismic exploration
US4881209A (en) * 1988-10-12 1989-11-14 Mobil Oil Corporation Method of processing converted mode seismic data
GB2275337B (en) 1993-02-17 1997-01-22 Csm Associates Limited Improvements in or relating to seismic detectors
US5996726A (en) * 1998-01-29 1999-12-07 Gas Research Institute System and method for determining the distribution and orientation of natural fractures
FR2800170B1 (fr) * 1999-10-22 2002-01-11 Geophysique Cie Gle Perfectionnement aux procedes de traitement sismique
GB0018480D0 (en) * 2000-07-27 2000-09-13 Geco Prakla Uk Ltd A method of processing surface seismic data
US6611764B2 (en) * 2001-06-08 2003-08-26 Pgs Americas, Inc. Method and system for determining P-wave and S-wave velocities from multi-component seismic data by joint velocity inversion processing
US6760667B1 (en) * 2001-09-14 2004-07-06 Emerald Geoscience Research Corporation Method for generating P-S and S-S seismic data and attributes from P-P seismic data

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6067275A (en) * 1997-12-30 2000-05-23 Schlumberger Technology Corporation Method of analyzing pre-stack seismic data
WO2001031364A1 (fr) * 1999-10-28 2001-05-03 Compagnie Generale De Geophysique Procede de prospection sismique mettant en oeuvre un traitement sur les ondes converties
WO2002059647A1 (en) * 2001-01-25 2002-08-01 Westerngeco Seismic Holdings Limited A method of processing marine seismic data and a method of seismic surveying

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
GAISER J.E.: "Multicomponent Vp/Vs correlation analysis". GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS. TULSA, US, vol. 61, no. 4, July 1996 (1996-07), pages 1137-1149, XP002143571, ISSN: 0016-8033, the whole document *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104407381A (zh) * 2014-12-17 2015-03-11 中国石油大学(华东) 从井中弹性波速的径向变化获取地层岩石脆裂性质的方法
CN104407381B (zh) * 2014-12-17 2017-01-18 中国石油大学(华东) 从井中弹性波速的径向变化获取地层岩石脆裂性质的方法

Also Published As

Publication number Publication date
GB0302898D0 (en) 2003-03-12
US7457195B2 (en) 2008-11-25
NO20053954D0 (no) 2005-08-25
EP1590689A1 (en) 2005-11-02
WO2004070424A1 (en) 2004-08-19
CA2515345A1 (en) 2004-08-19
GB2398124B (en) 2006-10-25
CA2515345C (en) 2014-12-09
GB2398124A (en) 2004-08-11
US20060092765A1 (en) 2006-05-04
NO20053954L (no) 2005-11-04
EA200501270A1 (ru) 2006-06-30

Similar Documents

Publication Publication Date Title
EA007881B1 (ru) Оценивание времени вступления сейсмической волны
Ochoa et al. Fast magnitude determination using a single seismological station record implementing machine learning techniques
US6791900B2 (en) Method of calculating a throw volume for quantitative fault analysis
Mahanti et al. Small lunar craters at the Apollo 16 and 17 landing sites-morphology and degradation
US5568450A (en) System and processor for real-time extraction of ocean bottom properties
US5475651A (en) Method for real-time extraction of ocean bottom properties
Kummerow Using the value of the crosscorrelation coefficient to locate microseismic events
US7720609B2 (en) Method of seismic signal processing
Küperkoch et al. Automated determination of S-phase arrival times using autoregressive prediction: application to local and regional distances
Andrioli et al. Improved analysis of all-sky meteor radar measurements of gravity wave variances and momentum fluxes
Tsuruoka et al. GRiD MT (grid-based real-time determination of moment tensors) monitoring the long-period seismic wavefield
Radi et al. A multisignal wavelet variance-based framework for inertial sensor stochastic error modeling
RU2014135771A (ru) Способ и устройство для обработки сейсмических данных
Tarpø et al. Automated reduction of statistical errors in the estimated correlation function matrix for operational modal analysis
KR101635791B1 (ko) P파 도달시간을 이용한 진원 및 진원시 결정방법
US20200271810A1 (en) Hybrid Residual Moveout Error Estimation
CN108304869A (zh) 一种基于多个传感器的综合幅值信息的融合方法及系统
Coll et al. Positioning with stationary emitters in a two-dimensional space-time
Khadhraoui et al. Real-time detection and localization of microseismic events
Zhang et al. On improving the precision of localization with gross error removal
JP3810323B2 (ja) 時間差方位検出装置
KR101035295B1 (ko) 확률 분포를 이용한 펄스열 분리 방법
Kouba et al. Ionospheric drift measurements: Skymap points selection
Mellor et al. The Modified Turning Bands (MTB) model for space-time rainfall. III. Estimation of the storm/rainband profile and a discussion of future model prospects
CN111221043A (zh) 一种被动源电场法的工作参数优化方法

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 BY KG MD TJ TM

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): KZ

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): AZ RU