Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,441

О ВОЗМОЖНОСТИ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ СЕЙСМО-АКУСТИЧЕСКИХ ВОЛН МЕТОДОМ ОБРАТНЫХ ЗАДАЧ

Мкртчян А.Г. 1 Мкртчян А.Р. 1 Котанджян Х.В. 1 Асланян А.А. 1 Арутюнян Э.М. 1 Арутюнян А.С. 1 Айвазян Г.А. 1 Мурадян О.Р. 1 Агабекян В.Н. 1 Миракян С.А. 1
1 Институт прикладных проблем физики Национальная академия наук Республики Армения
Работа посвящена определению основных параметров сейсмо-акустических волн, распространяющихся в земной коре, путем решения обратной задачи. Экспериментальные работы проводились на естественных и искусственных водоемах сотрудниками ИППФ НАН РА. С целью эффективной регистрации сейсмо-акустических волн были разработаны новые сверхчувствительные системы регистрации на основе явления двойной модуляции мессбауэровского излучения. Проведенный анализ показывает, что, применяя представленный математический аппарат при обработке экспериментальных данных временной мессбауэровской спектроскопии путем решения обратной задачи, можно определить основные параметры низкочастотных сейсмо-акустических колебаний и определить местонахождения источников сейсмо-акустических колебаний. Точность определения исходных параметров зависит от реальных условий экспериментальных исследований и погрешностей экспериментальной установки, которые могут привести к соответствующим статистическим ошибкам. Предварительная оценка сравнения теоретических расчетов с полученными экспериментальными данными дает хорошее согласование.
двойная модуляция
мессбауэровское излучение
сейсмо-акустические волны
амплитудно-временная спектроскопия
обратная задача
1. Вертхейм Т.К. Эффект Мессбауэра. – М.: Мир, 1966. – С. 9–24.
2. Мкртчян А., Габриелян Р., Шагинян С., Мартиросян А. О некоторых обратных задачах в спектроскопии // Изв. АН Арм. ССР., физика. – 1983. – Т. 18, Вып. 3. – С. 3–13.
3. Мкртчян А., Габриелян Р. Изменение мессбауэровских спектров поглощения, индуцируемые низкочастотными возбуждениями разных форм // Акустический журнал. – 1980. – Т. 26, Вып. 2. – С. 200–206.
4. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. – М.: Наука, 1986. – С. 288.
5. Mkrtchyan A., Nadjarian G., Gabrielyan R., Khanamiryan R. – Inverse Problem in Mössbauer Spectroscopy in case of Low Frecuency // Phys. Stat. Sol. – 1982. – b. 109. – Р. 131–136.
6. Mkrtchyan A., Gabrielyan R., Kocharyan L., Tumanian S. – Double gamma-acoustic resonance in crystals // Phys. Stat. Sol. – 1976. – b. 73 – Р. 681–687.
7. Mkrtchyan A., Harutyunyan G., Pulatov M., Mkrtchyan A. On Determination of Parameters of Mössbauer Absorption Spectra at Ultrasonic Excitation // Proceedings of NAS RA, Physics. – 1978. – Vol. 13, Iss. 3. – Р. 217–220.

В настоящей работе проводится определение основных параметров распространяющихся в земной коре сейсмо-акустических волн путем решения обратной задачи [2, 5]. Для регистрации сейсмо-акустических волн на основе явления двойной модуляции мессбауэровского излучения [6] в ИППФ НАН РА были разработаны уникальные детекторы. Полученные результаты теоретических расчетов достаточно близки с полученными экспериментальными данными.

Определение параметров и вида сейсмо-акустических волн на основании мессбауэровских спектров поглощения, модулированных акустическими колебаниями, часто сводятся к решению обратных задач. Трудности решения обратных задач во многих случаях обусловлены их некорректной постановкой [2, 4, 5].

Для определения закона возбуждения рассмотрим форму линии поглощения при внешнем периодическом возбуждении [1, 3, 6, 7]

makrchan01.wmf (1)

где F(x) – мессбауэровский спектр поглощения с некоторой ошибкой; φ(s) – закон возбуждения,

makrchan02.wmf

Умножая обе части уравнения (1) на exp(+ipx), интегрируя по x и используя обратное Фурье преобразование, получаем

makrchan03.wmf (2)

Однако вычисление Ψ(s) из (2) представляет определенные трудности, т.к. для сходимости интеграла необходимо, чтобы экспериментальная функция F(x) была гладкой функцией. Даже маленькая ошибка в F(x) приводит к большому искажению Ψ(s) из-за члена exp|p|. Эти трудности можно преодолеть при помощи метода регуляризации [5], развитой для решения некорректно поставленных обратных задач.

Для определения значения Ψ(s) рассмотрим «регуляризационный» функционал

makrchan04.wmf (3)

Функция Ψ1 является регуляризованным решением и лишена вышеуказанных недостатков.

Таким образом, учитывая, что функционал J имеет минимум и что Ψ1 (±∞) = 0 или Ψ′1(±∞) = 0, из условия минимума J имеем

makrchan05.wmf (4)

Для ψ1(s) получаем

makrchan06.wmf (5)

Как и ожидалось, это выражение совпадает с выражением (2) при условии C = 0. Анализируя далее экспериментальные случаи с четной функцией F(x), имеем

makrchan07.wmf (6)

где makrchan08.wmf которая вычисляется метод Филона и значение C определяется из экспериментальной ошибки.

При некоторых асимметричных случаях закона возбуждения задача сводится к решению интегрального уравнения

makrchan09.wmf

который решается численными методами.

Метод математической обработки модулированных временных мессбауэровских спектров для определения неизвестных параметров внешних возбуждений заключается в согласовании экспериментальных точек Nk = N(tk) некой аналитической функцией Ñ = Ñ(tk; α1, α2,…, αn), где tk – время, α1, α2,…, αn – неизвестные параметры внешних возбуждений. Точное определение параметров α1, α2,…, αn сводится к подбору согласующей функции для получения наилучших совпадений экспериментальных точек Nk с согласующими функциями Ñ(tk; α1, α2,…, αn).

Согласующую функцию для временных мессбауэровских спектров можно описать следующим выражением

makrchan10.wmf (7)

Для подбора параметров аппроксимирующей функции минимизируется следующий неотрицательный функционал

makrchan11.wmf

где n – число экспериментальных точек, для которых выполняется условие

makrchan12.wmf (8)

Отсюда для параметров αi получаем следующую систему уравнений

makrchan13.wmf

makrchan14.wmf

makrchan15.wmf

makrchan16.wmf (9)

Получается система нелинейных алгебраических уравнений для параметров α1, α2, α3 и α4, при решении которой используется метод квазилинеаризации Ньютона – Рапсона, т.е. выбираются начальные приближенные значения makrchan17.wmf и разлагается fji) в ряд в окрестности точки makrchan18.wmf. Оставляя только линейный член, получаем

makrchan19.wmf i = 1, 2, 3, 4, (10)

где makrchan20.wmf – матрица Якоби.

Решение этой системы makrchan21.wmf является первым приближением параметров αi

makrchan22.wmf (11)

где makrchan23.wmf – обратная матрица Якоби.

Найденные таким образом решения makrchan24.wmf снова подставляются в первоначальную систему уравнений (9) для нахождения второго приближения makrchan25.wmf. Решения makrchan26.wmf снова подставляются в (9) и находится следующее приближение и т.д. Окончательно для α1, α2, α3 и α4 в n-м приближении имеем

makrchan27.wmf, i,j = 1, 2, 3, 4. (12)

pic_29_1.tifа

pic_29_2.tifб

Закон движения источника гамма-квантов (а) и временной мессбауэровский спектр поглощения (б): V – скорость движения источника гамма-квантов; t – время; N – количество гамма-квантов, зарегистрированных детектором

Точность определения исходных параметров makrchan18.wmf зависит от реальных условий эксперимента и полученных экспериментальных данных, так как в реальном физическом эксперименте величины Nk определяются с некоторой погрешностью, определяемой статистическим набором гамма-квантов, которая и может привести к соответсвующим ошибкам в определении αi.

На рисунке а приведено заданное нами движение источника гамма-квантов, а на рисунке б – временной мессбауэровский спектр поглощения гамма-квантов, регистрированных с помощью детекторов, разработанных в ИППФ НАН РА. Как видно из этих рисунков, имеется достаточно хорошее совпадение измеряемого сигнала с законом движения источников гамма-квантов.

Заключение

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

Рецензенты:

Ванцян А.А., д.ф.-м.н., профессор, ведущий научный сотрудник, Институт прикладных проблем физики Национальной академии наук Республики Армения, г. Ереван;

Саарян А.А., д.ф.-м.н., профессор, Ереванский государственный университет, г. Ереван.

Работа поступила в редакцию 15.05.2014.


Библиографическая ссылка

Мкртчян А.Г., Мкртчян А.Р., Котанджян Х.В., Асланян А.А., Арутюнян Э.М., Арутюнян А.С., Айвазян Г.А., Мурадян О.Р., Агабекян В.Н., Миракян С.А. О ВОЗМОЖНОСТИ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ СЕЙСМО-АКУСТИЧЕСКИХ ВОЛН МЕТОДОМ ОБРАТНЫХ ЗАДАЧ // Фундаментальные исследования. – 2014. – № 8-1. – С. 51-54;
URL: http://www.fundamental-research.ru/ru/article/view?id=34506 (дата обращения: 13.04.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074