Способ обработки сейсмической информации

О П И СА Н И Е ИЗОБРЕТЕНИЯ К АВТОРСКОМУ СВИДЕТЕЛЬСТВУ 33 бб 28 Союз Советских Социалистических РеспубликЗависимое от авт, свидетельствач, б 01 ч 1/28 аявлено 05,Ч 1.1970 ( 1446915/26 с присоединением заявки1 риоритетОпубликовано 21 1 Ч.1972, БюллетеньДата опубликования описания 29 Л.1972 Комитет по деламизобретений и открытипри Совете МинистроеСССР К 550.834(088.8) Авторыизобретения С. А. Кац, В. Л. Киселевич и Б. М, Шубик Институт физики Земли им. О. Ю. Шмидта Заявител СПОСОБ ОБРАБОТКИ СЕЙСМИЧЕСКОЙ ИНФОРМА Изобретение относится к сеисморазведке, а именно к методике принятия решения о наличии регулярного сигнала в сейсмограмме,Обычно решение о наличии сигнала с определенным направлением прихода принимается 5 интерпретатором в результате прослеживания сигнала или какой-либо фазы сигнала по некоторому числу каналов сейсмограммы (выделение осей синфазности).Надежность выделенных осей синфазности 10 зависит от числа каналов, соотношения энергий и спектральных характеристик сигналов и помех и падает при увеличении энергии помех. В случаях, когда отношение сигнал/помеха мало, надежное выделение оси синфаз ности при ручной обработке затруднительно, Кроме того, при ручной обработке нет возможности оценить количественно надежность принятого решения о наличии регулярного сигнала. 20Цель изобретения - автоматизация процесса принятия решения, повышение надежности выделенных осей синфазности на сейсмограмме и получение количественной характеристики надежности принятого решения. 25Цель достигается тем, что решение о наличии волны в данный момент с данным направлением прихода принимают, если на нескольких каналах, число которых превышает оитическое, амплитуда сигнала в пределах 30 выбранного временного окна превосходит подобранный порог,Порог изменяют методом последовательных приближений до тех пор, пока вычисляемая по сейсмограмме оценка вероятности принятия ошибочного решения, соответствующая этому порогу, не оказывается внутри заданного, интервала значений, Надежность выделенных осей синфазности количественно оценивается вероятностью принятия ошибочного решения, соответствующей подобранному порогу.Такое решение может быть принято с помощью автоматического многоканального амплитудного порогового обнаружителя коррелированных сигналов с плавающим порогом. Критерием для подбора порога является попадание оценки вероятности ошибки, соответствующей порогу, в заданный интервал значений,Анализируемая сейсмограмма содержит Л каналов. Если иа этой сейсмограмме в некотором временном окне с центром в моментнет регулярной волны с заданным направлением прихода, то считают, что вероятность принятия ошибочного решения о наличии волны в этот момент равна вероятности случайного превышения амплитудной записи некоторого порога д не менее, чем на и каналах из Л (л (Л) в этом Временном окне.336628 М/=1МФ/=115 Составитель В. Гринь Техред Л. Евдонов Корректор Л. Бадылама Редактор И. Орлова Заказ 1499/2 Изд. Мо 607 Тираж 448 Подписное ЦНИИПИ Комитета по делам изобретений и открытий при Совете Министров СССР Москва, Ж.85, Раушская паб., д, 4/5 Типография, пр. Сапунова, 2 Вероятность такого события определяется по формуле биномиального распределения:Р(с/, и, У) =; Сдп Ры(1 - Р )(1)1=огде Р(д, п, Л) - вероятность случайного превышения порога не менее, чем на и каналах из Ж;Ра - вероятность случайного превышения порога д на одном канале;Су - число сочетаний из У по /; У - общее число каналов;и - критическое число каналов,по которому принимают решение;Ы - порог.Таким образом, задав порог д, можно вычислить соответствующую этому порогу величину вероятности принятия ложного решения при фиксированном критическом числе каналов.Путем направленного перебора можно подобрать порог д такой величины, чтобы соответствующая этому порогу вероятность ошибки попала в заданный интервал значений. Если теперь в окне с центром в момент 1 по заданному направлению прихода зафиксировано превышение этого порога д не менее, чем на и критических каналах из /Ч, то принимают решение о наличии в этот момент регулярной волны; при этом вероятность ошибки равна вероятности случайного превышения порога с/о не менее, чем на и критических каналах, Р(йо, и, У), и вычисляется по формуле (1)Трудность заключается в том, что для данного порога неизвестна величина Р вероятности случайного превышения порога на одном канале. По предлагаемому способу обработки вычисляют величину Рт непосредственно по сейсмограмме. Для этого обрабатываемую сейсмограмму разбивают на М смежных временных окон по такому направлению прихода, по которому заведомо отсутствуют регулярные волны, и для данного порога в каждом окне вычисляют величину.(2) лгде Рг/; - отношение числа каналов, на которых есть превышение порога д к общему числу каналов в /-ом окне; К; - число каналов, на которых есть превышение порога в /-ом окне;М - число смежных окон на сейсмограмме. Эта величина является оценкой вероятности случайного премышения порога на одном канале в /-ом окне. Усреднив эту оценку по всем М окнам, получают оценку вероятности 10 превышения данного порога на одном каналепо всей сейсмограмме: Эта оценка Ра используется для вычислениявероятности ошибки вместо неизвестной величины Ре.Таким образом, процедура обработки сей смограммы состоит из двух этапов. Сначалапутем направленного перебора подбирается порог 4, обеспечивающий заданный уровень вероятности ошибки. Затем по сейсмограмме в текущем окне,по заданному паправлению 25 прихода подсчитывается число каналов, накоторых есть превышение подобранного порога. Если в момент 1 зафиксировано превышение подобранного порога не менее, чем на заданном критическом числе каналов, то при нимается решение о регистрации в этот момент регулярной волны. Надежность принятого решения оценивается величиной Н = (1 - 1 ) 100 /о, где Р - вероятность принятия ошибочного решениясоответствующая по 35 рогу.Если направление прихода неизвестно, тоосуществляют перебор направлений из некоторого диапазона. 40 Предмет изобретения Способ обработки сейсмической информации путем выделения осей синфазности на сейсмограмме, отличающийся тем, что, с 45 целью повышения надежности выделенныхосей синфазности и получения количественной характеристики надежности принятого решения, решение о наличии волны в данный момент с данной кажущейся скоростью прини мают, если на нескольких каналах, число которых превышает критическое, амплитуда сигнала в пределах выбранного интервала времени превосходит хотя бы в одной точке подобранный порог,