ОЦЕНКА ВОЗМОЖНОСТИ ПОВЫШЕНИЯ ЭФФЕКТИВНОСТИ ИНЖЕНЕРНОЙ СЕЙСМОРАЗВЕДКИ ПУТЁМ РЕАЛИЗАЦИИ КОДОИМПУЛЬСНОГО НАКАПЛИВАНИЯ КОЛЕБАНИЙ

Ю.П.Кострыгин, М.Н.Удалый

Аннотация

В данной статье рассмотрена возможность реализации для инженерной сейсморазведки кодоимпульсного накапливания колебаний с использованием маломощных малогабаритных переносных кодоимпульсных излучателей. Показано, что для реализации кодоимпульсного метода с целью решения задач инженерной сейсморазведки весьма перспективным является применение однополярных кодовых функций ЛИЧ с максимальной частотой следования импульсов, равной 150-200 Гц, либо двухполярных дополнительных последовательностей с раздельной отработкой импульсов с положительной и отрицательной полярностью.

Очевидно, что эффективное решение задач, стоящих перед инженерной сейсморазведкой, представляется возможным лишь при использовании высокочастотных зондирующих сигналов, позволяющих существенно повысить разрешающую способность метода. При этом многочисленные теоретические и экспериментальные исследования указывают на существование устойчивой тенденции к повышению спектрального состава формируемых колебаний по мере уменьшения энергии воздействия импульсного источника на грунт [8]. В частности, поэтому в инженерной сейсморазведке для возбуждения упругих колебаний, в основном, применяются маломощные излучатели типа ручной кувалды массой 4 кг и 8 кг, характеризующиеся энергией воздействия, равной 60 – 140 Дж [7].

Использование в качестве источника сигналов ручной кувалды позволяет существенно расширить спектр колебаний и в значительной степени облегчает проведение работ, например, в горных условиях, где применение транспортных средств становится проблематичным. Вместе с тем сейсмические импульсы, возбуждаемые кувалдой, в принципе не могут иметь стабильную динамику, что в значительной степени снижает эффективность сейсмических наблюдений по технологии многократных перекрытий. Более того, во многих случаях даже при кратностях статистического накапливания 12-15 при работе с кувалдой не удаётся обеспечить необходимую помехоустойчивость сейсмического метода.

В этой связи безусловный интерес представляют исследования возможности применения для инженерной сейсморазведки маломощных малогабаритных переносных кодоимпульсных излучателей. Такие излучатели обеспечили бы формирование более стабильных по динамике мало амплитудных широкополосных сигналов. Вместе с тем помехоустойчивость кодоимпульсного метода, как известно, на порядок превышает помехоустойчивость сейсмического метода в режиме редко следующих импульсов [2]. Указанные, весьма значимые преимущества технологии сейсмических наблюдений с использованием маломощных кодоимпульсных излучателей, позволили бы существенно повысить эффективность инженерной сейсморазведки.

При построении технологии сейсмических наблюдений с использованием кодоимпульсного метода важнейшим является вопрос по выбору кодовых функций последовательностей. На основании выполненных к настоящему времени исследований показано, что для нефтегазовой сейсморазведки наибольший интерес представляют однополярные последовательности с линейным по времени изменением частоты следования импульсов в развёртках (коды ЛИЧ), с линейным изменением периода следования импульсов (коды ЛИП), а также двухполярные дополнительные импульсные последовательности [1, 2, 6]. Вместе с тем очевидно, что для инженерной сейсморазведки параметры эффективных кодовых посылок могут отличаться от соответствующих параметров для нефтегазовой сейсморазведки, т.к. при инженерных изысканиях характеристики волновых полей существенно отличаются от аналогичных характеристик при проведении сейсмических исследований на нефть и газ.

Так, на рис. 1 приведены сейсмограммы и соответствующие им спектры, полученные в процессе инженерных изысканий в районе Сочи с использованием вертикальных (а) и горизонтальных (б) сейсмоприёмников. Для возбуждения колебаний в данном случае применялась кувалда массой 4 кг. Из рис. 1 следует, что в отличие от нефтегазовой сейсморазведки при инженерных изысканиях в связи с небольшим расстоянием от точек приёма до источника продольные полезные колебания (Р) характеризуются существенно более высокочастотным составом, что, безусловно, следует учитывать при выборе кодовых функций. В конкретном случае, преобладающая частота колебаний равнялась 93 Гц, а правая граничная частота на уровне 0,1 от максимального значения спектра, 175 Гц.

В непосредственной близости от продольных полезных колебаний на сейсмограмме, представленной на рис. 1,а можно видеть многофазный цуг низкоскоростных поверхностных волн R. Такие волны присутствуют при инженерных изысканиях повсеместно, характеризуются повышенной амплитудой и по этой причине при реализации кодоимпульсного метода могут формировать интенсивный корреляционный фон, снижающий качество выделения относительно слабых полезных колебаний. Указанное обстоятельство также следует учитывать при определении технологии кодоимпульсного накапливания колебаний для инженерной сейсморазведки.

И, наконец, важно отметить, что при проведении инженерных изысканий методами сейсмической разведки, как правило, предполагается исследование не только структурных форм геологических объектов, но и динамических модулей упругости пород. Для решения таких задач существует необходимость возбуждения и регистрации не только продольных, но и поперечных колебаний. Это обстоятельство также необходимо учитывать при определении методики инженерных изысканий с использованием кодоимпульсного метода.

В данной статье по результатам проведённых модельных исследований обсуждается возможность решения задач инженерной сейсморазведки путём использования кодоимпульсного метода с применением однополярных импульсных последовательностей ЛИЧ, однополярных последовательностей ЛИП и двухполярных дополнительных последовательностей [1, 2, 4, 5, 6]. При моделировании кодоимпульсные развёртки были составлены из импульсов Берлаге (1)

F (t) = a∙tⁿ∙ e ^-^βt sin2πf_видt t ≥ 0 (1)

где f_вид – видимая (преобладающая) частота колебаний, β – коэффициент затухания импульса.

При этом для продольных волн f_вид = 93 Гц, β = 300, n = 2,1, а для низкоскоростных поверхностных волн f_вид = 26 Гц, β =90 , n =1,9 . Выбранные значения указанных параметров обеспечивали наибольшее соответствие спектров наблюдённых и теоретических импульсов.

Временное сжатие импульсных последовательностей X(t) осуществлялось путём корреляционной свёртки X(t) * K(-t), где K(t) – функция кода последовательности.

Для кодов ЛИЧ с начальной и конечной частотами следования F_н и F_к при длительности развёртки Т, временное положение n-го импульса, т.е. функция кода, определялась исходя из известного выражения [2, 4]

_____________

t_n= (T/DF)(ÖF_н²+ 2nDF/T –F_н), (2)

где t_n –временное положение n-го импульса в развёртке (начальному импульсу присвоен номер «0»), DF – ширина частотного диапазона развёртки.

При этом результаты исследований, приведённые в работе [4], показывают, что сохранить высокочастотные гармоники, участвующие в формировании единичных импульсов, входящих в развёртку ЛИЧ, представляется возможным лишь при условии, если максимальная частота развёртки располагается в области правой граничной частоты спектра единичных импульсов. Принимая же во внимание относительно высокочастотный состав продольных колебаний при инженерных изысканиях можно предположить, что максимальная частота развёртки ЛИЧ для инженерной сейсморазведки с использованием продольных волн должна соответствовать диапазону 150-200 Гц.

Так, на рис. 2 приведены главные максимумы ФВК импульсных последовательностей ЛИЧ и соответствующие им спектры, рассчитанные при различных значениях F_max. Расчёты в данном случае были выполнены с учётом динамических характеристик полезных продольных колебаний, регистрируемых при инженерных изысканиях. Максимальная частота развёртки F_max изменялась от 60 Гц до 200 Гц, а в качестве единичного импульса последовательности использовался импульс Берлаге с параметрами: f_вид = 93 Гц, β = 300, n = 2,1.

Из рисунка 2 следует, что для высокочастотных продольных колебаний с верхней граничной частотой f_гр ≈ 175 Гц спектр ФВК кодоимпульсных последовательностей ЛИЧ при F_max = 60 Гц, т.е. при значительном несоответствии спектральных характеристик функции кода и единичного импульса, имеет весьма сложную форму и два экстремума, сформированные соответственно на частотах ≈60 Гц и ≈110 Гц. С увеличением F_max форма спектра ФВК становится более простой, однако из записи практически исключаются гармонические составляющие, расположенные правее F_max. Таким образом, очевидно, что эффективное применение кодоимпульсного накапливания колебаний с использованием развёрток ЛИЧ представляется возможным при условии, если максимальная частота следования F_max равна или превышает правую граничную частоту f_гр спектра единичных импульсов. Соблюдение указанного условия позволяет сохранить высокочастотную часть спектра продольных колебаний при инженерных изысканиях и сформировать предельно компактные корреляционные импульсы (см. рис. 2).

Более того, в связи с тем, что спектры функций кода ЛИЧ, как известно, характеризуются равномерным увеличением спектральной плотности с ростом частоты [3], корреляционные импульсы для последовательностей ЛИЧ при достаточно больших значениях F_max ≥ f_гр отличаются несколько более высокочастотным составом по сравнению с единичными импульсами, формирующими эти последовательности.

Итак, высокую эффективность инженерных изысканий с использованием кодоимпульсных последовательностей ЛИЧ, можно обеспечить при условии, если применяемые сейсмические излучатели позволяют формировать импульсные последовательности с максимальной частотой следования 150-200 Гц. Очевидно, что разработка и создание таких излучателей является достаточно сложной технической задачей. Однако, как уже отмечалось, кодоимпульсный метод по сравнению с методом статистического накапливания в режиме редко следующих импульсов, отличается значительно более высокой помехоустойчивостью, что в принципе позволяет без ущерба для качества сейсмических записей существенно снизить динамические нагрузки на грунт и тем самым существенно упростить задачу по разработке высокочастотного излучателя.

Так, если предположить, что кодоимпульсный излучатель развивает такое же динамическое усилие на грунт, как и в случае удара кувалды, то при кодоимпульсном накапливании с использованием последовательности ЛИЧ с частотным диапазоном 10-180 Гц и длительностью 15 с помехоустойчивость сейсмического метода будет в ≈ 10 раз превышать помехоустойчивость метода при 15-кратном возбуждении колебаний кувалдой. Т.е. динамическое усилие, развиваемое кодоимпульсным излучателем, можно без ущерба для качества, как минимум, уменьшить в 5 раз по сравнению с силовым воздействием кувалдой на грунт. Учитывая, что силовое воздействие на грунт при ударе кувалдой соответствует диапазону 15-25 кН [7], усилие, развиваемое кодоимпульсным излучателем можно ограничить значениями 3-5 кН.

Таким образом, кодоимпульсный излучатель, используемый для решения задач, стоящих перед инженерной сейсморазведкой, при возбуждении импульсных последовательностей ЛИЧ должен обеспечивать максимальную частоту следования воздействий, равную 150-200 Гц при динамическом усилии на грунт 3-5 кН.

Существенно уменьшить максимальную частоту следования сейсмических импульсов в однополярных кодоимпульсных развёртках при проведении инженерных изысканий можно путём использования сигналов с линейным изменением периода следования воздействий (коды ЛИП). Величина периода следования импульсов Т_n в развёртке ЛИП определяется условием

Т_n=T₀- ∆T(n-1), (3)

где T₀– начальный период следования импульсов;

∆T – приращение периода следования импульсов;

n - номер импульса в развёртке.

Используя формулу суммы членов арифметической прогрессии и полагая, что начальному импульсу присвоен номер «0» можно определить относительное запаздывание t_n для n-го импульса в развёртке ЛИП, т.е. функцию кода последовательности

t_n = n T₀- n(n-1) ∆T/2 (4)

Как было показано в работе [5] в отличие от импульсных последовательностей ЛИЧ, последовательности ЛИП могут обеспечивать формирование компактных корреляционных импульсов, характеризующихся широкополосными спектральными характеристиками, при использовании кодовых функций с существенно меньшими значениями максимальной частоты развёртки F_max. Так, на рис. 3 приведены главные максимумы ФВК импульсных последовательностей ЛИП и соответствующие им спектры, рассчитанные при различных значениях F_max. Материалы, представленные на рис. 3, подтверждают возможность весьма значительного уменьшения параметра F_max при проведении инженерных изысканий с применением кодоимпульсного метода, основанного на использовании последовательностей ЛИП. Действительно, как можно видеть из рис. 3, даже при условии F_max= 30 Гц главный максимум и амплитудный спектр ФВК импульсной последовательности ЛИП практически не отличаются от аналогичных характеристик единичных импульсов последовательности. Некоторые, весьма не существенные различия, заключаются в незначительном увеличении веса низкочастотных составляющих для корреляционных импульсов последовательностей ЛИП. Физически это связано с тем, что значения спектральной плотности функции кода ЛИП при увеличении частоты несколько уменьшаются [3].

Очевидно, что обсуждаемая возможность уменьшения параметра F_maxвследствие применения при инженерных изысканиях последовательностей ЛИП, позволяет весьма существенно упростить разработку кодоимпульсных излучателей.

Как уже отмечалось выше, на сейсмических записях, регистрируемых при инженерных изысканиях, в непосредственной близости от продольных полезных колебаний присутствуют многофазные цуги высокоамплитудных низкоскоростных поверхностных волн. Такие помехи могут формировать интенсивный корреляционный фон, снижающий качество выделения более слабых продольных колебаний. В этой связи при выборе наиболее приемлемых для инженерных изысканий кодовых функций последовательностей следует учитывать не только динамику корреляционных импульсов полезных сигналов, но и степень мешающего влияния помех корреляционного преобразования, соответствующих наиболее интенсивным помехам.

С целью оценки мешающего влияния на полезную запись корреляционного фона, связанного с доминирующими поверхностными волнами, были рассчитаны графики динамического диапазона корреляционного преобразования D(τ) для последовательностей ЛИЧ и ЛИП (см. рис. 4). Единичные импульсы в развёртках задавались в данном случае импульсами Берлаге с параметрами f_вид = 26 Гц, β = 90, n = 1,9. Динамический диапазон преобразования рассчитывался по формуле

D(τ) = 20 lg А_max/σ(τ) , (5)

где А_max – пиковая амплитуда главного максимума ФВК кодоимпульсного сигнала; σ(τ) – среднеквадратичная амплитуда корреляционного фона во временном окне 50 мс. Значения D(τ) определялись вдоль корреляционной функции с дискретностью 1 мс.

Из рис. 4, а можно видеть, что для наиболее оптимальных для инженерных изысканий импульсных последовательностей ЛИЧ значения динамического диапазона помех корреляционного преобразования, формируемых низкоскоростными помехами, соответствуют 40 – 53 дБ при τ = 0,11 – 0,6 с, что в принципе позволяет использовать такие последовательности при условии, если амплитуда помех не более, чем на ≈ 30 дБ превышает уровень полезных сигналов.

Для кодоимпульсных же сигналов ЛИП в аналогичном временном интервале значения динамического диапазона преобразования помех соответствуют 25 – 50 дБ (см. рис. 4, б), что не позволяет рассчитывать на эффективное использование кодоимпульсного метода для решения задач инженерной сейсморазведки. Так, например, для сейсмической записи, приведённой на рис. 1, соотношения между амплитудами продольных полезных колебаний, регистрируемых на временах 50 – 150 мс, и амплитудами доминирующих поверхностных волн приближались к 30 дБ. Т.е. в случае использования в таких, совершенно типичных условиях, импульсных последовательностей ЛИП относительно слабые полезные сигналы были бы соизмеримы по интенсивности с корреляционным фоном, формируемым поверхностными волнами.

Уменьшить интенсивность корреляционного фона при кодоимпульсном накапливании можно, как известно, путём применения двухполярных дополнительные последовательностей [1].

Импульсные последовательности F_п(t) и ¢F_п(t) называются дополнительными, если им соответствуют функции кода K(t) и ¢K(t), для которых выполняется условие

2 при t=0,

U₀(t) + ¢U₀(t) = (6) 0 при t= ±1, ±2,…, ±(N-1),

где U₀(t) и ¢U₀(t) – нормированные ФАК функций K(t) и ¢K(t) .

Т.е. для таких последовательностей сумма ФАК функций кода равна нулю везде, кроме t=0.

Как уже отмечалось, дополнительные последовательности относятся к классу так называемых двухполярных последовательностей. Вместе с тем используемые в сейсморазведке импульсные излучатели являются излучателями одностороннего действия. В работе [6] был предложен весьма простой способ, позволяющий формировать двухполярные последовательности излучателями одностороннего действия. В предложенном способе разнополярную последовательность импульсов реализуют односторонними воздействиями за два сеанса. Сначала возбуждается развёртка, соответствующая временному закону следования импульсов с положительной полярностью, затем развёртка, соответствующая закону следования импульсов с отрицательной полярностью. Вычитая вторую запись из первой, получают аналог виброграммы, которая бы регистрировалась при использовании кодоимпульсных излучателей двухстороннего действия силы.

Оценивая эффективность применения дополнительных последовательностей для решения задач инженерной сейсморазведки необходимо учитывать, что в отличие от «идеальных» кодоимпульсных сигналов, характеризующихся постоянной амплитудой единичных импульсов, последовательности, формируемые реальными сейсмическими источниками, обладают некоторой амплитудной неидентичностью импульсов в развёртке, зависящей от частоты следования F. Так, например, экспериментальная зависимость А_имп (F), полученная для серийного кодоимпульсного источника ИКИ-10/40, удовлетворительно аппроксимируется при F ≥ 5 Гц степенной функцией [2]

А_имп=k·F^-^p. (7)

Вместе с тем, при раздельной отработке импульсов с положительной и отрицательной полярностью, входящих в дополнительные развёртки, частота следования не остаётся постоянной, что приводит к амплитудной неидентичности, ориентировочно определяемой зависимостью (7). При такой технологии формирования дополнительных последовательностей амплитуду импульсов можно рассчитать исходя из выражения (8)

А_n~[1/(t_n-t_n_-1)] ^-^p, (8)

где n – номер импульсов для однополярных развёрток, получаемых вследствие раздельной отработки импульсов дополнительных последовательностей с положительной и отрицательной полярностью.

Очевидно, что в случае повышенных значений параметра р амплитудная неидентичность в двухполярных последовательностях, формируемых путём раздельной отработке импульсов с положительной и отрицательной полярностью, может приводить к уменьшению динамического диапазона корреляционного преобразования и к искажениям динамики главного максимума ФВК. Анализ взаимосвязи между амплитудной неидентичностью импульсов в развёртках, значениями динамического диапазона преобразования и динамикой главного максимума ФВК был выполнен с использованием дополнительных последовательностей K(t) и ‛K(t), каждая из которых включала в себя по 256 импульсов.

На рис. 5 представлены графики динамического диапазона корреляционного преобразования дополнительных последовательностей, формируемых путём раздельной отработки импульсов с положительной и отрицательной полярностью, рассчитанные для доминирующей поверхностной волны (f_вид = 26 Гц, β = 90, n = 1,9) при различных значениях параметра р. Из приведённого рисунка можно видеть, что для р ≤ 0,12-0,16, т.е. при условии если в частотном диапазоне развёртки 10 – 80 Гц амплитуда импульсов изменяется не более чем в 1,3 - 1,4 раза, или же в частотном диапазоне развёртки 10-30 Гц не более чем в 1,15 – 1,2 раза, динамический диапазон корреляционного преобразования доминирующих поверхностных волн в интервале τ = 0,3 – 1,0 с составляет не менее 65-62 дБ. Такие значения динамического диапазона преобразования в ближней зоне корреляционных импульсов поверхностных помех позволяют практически исключить мешающее влияние корреляционного фона на целевые отражённые волны при инженерных изысканиях.

При значениях же параметра р = 0,32, т.е. при двух - кратном изменении амплитуды импульсов в пределах частотного интервала развёрток 10-80 Гц, динамический диапазон преобразования дополнительных последовательностей для поверхностных помех в указанном временном отрезке составляет не менее 56 дБ, что также превышает динамический диапазон преобразования при работе с однополярными кодовыми функциями. Однако очевидно, что эффективность кодоимпульсного метода определяется не только уровнем корреляционного фона, но и динамикой главных максимумов ФВК полезных сигналов.

На рис. 6 представлены главные максимумы ФВК дополнительных последовательностей, и соответствующие им амплитудные спектры, рассчитанные для целевых отражений, регистрируемых при инженерных изысканиях (f_вид = 93 Гц, β = 300, n = 2,1). При моделировании значения параметра р изменялись в данном случае от нуля до 0,32, тактовый период следования импульсов τ₀ равнялся 30 мс. Анализ графиков, представленных на рис. 6, свидетельствует о том, что при использовании дополнительных последовательностей главные максимумы ФВК целевых отражений и соответствующие им спектры практически не искажаются если р ≤ 0,12-0,16. При более высоких же значениях параметра р наблюдается заметное усиление боковых пульсаций главных максимумов ФВК, что в свою очередь обуславливает заметные искажения формы амплитудных спектров корреляционных импульсов.

Заключение

1. Безусловный интерес представляют исследования возможности и целесообразности применения для инженерной сейсморазведки кодоимпульсного метода с использованием малогабаритного переносного излучателя, возбуждающего не только продольные, но и поперечные колебания. Технология кодоимпульсного накапливания колебаний, отличающаяся высокой помехоустойчивостью, позволила бы повысить эффективность инженерных изысканий путём формирования стабильных по динамике мало амплитудных широкополосных сигналов.

2. Для реализации кодоимпульсного метода с целью решения задач инженерной сейсморазведки весьма перспективным при возбуждении поля продольных колебаний является применение однополярных кодовых функций ЛИЧ с максимальной частотой следования импульсов, равной 150-200 Гц. При таких параметрах кодоимпульсного сигнала спектр корреляционных импульсов целевых отражений для продольных колебаний при инженерных изысканиях сохраняет высокочастотные гармонические составляющие, а динамический диапазон корреляционного преобразования для доминирующих поверхностных помех соответствуют 40 – 53 дБ при τ = 0,11 – 0,6 с, что позволяет использовать такие последовательности при условии, если амплитуда помех не более, чем на ≈ 30 дБ превышает уровень полезных сигналов.

3. Однополярные последовательности ЛИП отличаются от последовательностей ЛИЧ значительно более высоким уровнем корреляционного фона в ближней зоне ФВК. Так, в интервале времён τ = 0,11 – 0,6 с динамический диапазон преобразования низкоскоростных помех соответствуют 25 – 50 дБ, что не позволяет применять такие развёртки для решения задач инженерной сейсморазведки.

4. Уменьшить интенсивность корреляционного фона при кодоимпульсном накапливании можно, как известно, путём применения двухполярных дополнительных последовательностей с раздельной отработкой импульсов с положительной и отрицательной полярностью. Для формирования дополнительных последовательностей могут быть использованы излучатели, для которых амплитуда импульсов уменьшается не более чем в 1,3-1,4 раза с увеличением частоты следования от 10 Гц до 80 Гц, или не более чем в 1,15-1,2 раза с увеличением частоты следования от 10 Гц до 30 Гц. При таких параметрах излучателя динамика корреляционных импульсов дополнительных последовательностей, соответствующих целевым отраженным продольным волнам, регистрируемым при инженерных изысканиях, практически не искажается, а корреляционный фон, формируемый доминирующими поверхностными помехами, практически не оказывает мешающего влияния на полезную запись.

5. При использовании двухполярных дополнительных развёрток технический руководитель работ может задавать относительно низкие значения тактовой частоты, что существенно облегчает технические требования при разработке кодоимпульсных излучателей.

6. При формировании поля поперечных колебаний при инженерных изысканиях с применением кодоимпульсного метода перспективными являются двухполярные дополнительные последовательности наряду с однополярными развёртками ЛИЧ с максимальной частотой следования импульсов F_max= 60 – 70 Гц.

Литература

1. Варакин Л.Е., 1970, Теория сложных сигналов. М.: Сов. радио, 376 с.

2. Кострыгин Ю.П., 2002, Сейсморазведка на сложных сигналах. Тверь: Изд-во ГЕРС, 416 с.

3. Кострыгин Ю.П., Колесников Д.А., 2011, Оценка влияния амплитудной неидентичности импульсов в развёртках на эффективность кодоимпульсного накапливания сейсмических колебаний//Технологии сейсморазведки. №2. С. 97-108.

4. Кострыгин Ю.П., Колесников Д.А., 2010, Оценка возможности повышения эффективности кодоимпульсного накапливания сейсмических колебаний с использованием однополярных кодов ЛИЧ// Технологии сейсморазведки. №3. С. 110-114.

5. Кострыгин Ю.П., Колесников Д.А., 2010, Сравнительный анализ эффективности кодоимпульсного накапливания сейсмических колебаний при использовании однополярных кодов ЛИЧ и ЛИП// Технологии сейсморазведки. №4. С. 46-52.

6. Роман В.И., Шпортюк Г.А., Михайличенко О.Г., Королюк П.А., 1986, Перспективы применения дополнительных последовательностей импульсов в невзрывных методах сейсморазведки, ГСЗ и сейсмологии// Современные геодинамические процессы и их изучение в связи с проблемой прогноза землетрясений. Киев: Наук. думка, 22-26 Рудомаха Н.Н., Горшков.

7. Рудомаха Н.Н., 2012, Оценка эффективности сейсмоисточника "КВАНТ" при инженерной сейсморазведке // Материалы 2-й Международной научно-практической конференции "Сочи-2012", С. 4.

8. Шнеерсон М.Б., Майоров В.В., 1988, Наземная невзрывная сейсморазведка. М., Недра, 237 с.

КОРОТКО ОБ АВТОРАХ

КОСТРЫГИН Юрий Петрович – доктор технических наук, директор ООО «Новоросморгео». В 1967 г. окончил Грозненский нефтяной институт. Научные интересы - сейсморазведка с использованием сложных сигналов. Автор 90 опубликованных работ, в том числе двух монографий и 14 авторских свидетельств и патентов на изобретение.

E-mail: kostr_p@mail.ru

УДАЛЫЙ Михаил Николаевич – ведущий геофизик ООО «Ильинской геофизической экспедиции». В 2004 г. окончил Кубанский государственный университет. Область научных интересов – технология полевых сейсморазведочных работ. Автор двух опубликованных работ.

Контакты просим осуществлять через Удалого М.Н.

Тел. 8-918-451-77-01.

e-mai:, mihail-222@mail.ru