ИСТИНА

Представлено описание метода инверсии сейсмических данных, основу которых составляют волновые формы приёмных функций (ПФ): результат специальной обработки сейсмограмм удалённых землетрясений. ПФ содержат обменные фазы, образованные на границах среды под сейсмостанцией при прохождении первичной или вторичной волн (P-ПФ и S-ПФ, соответственно). Определение строения среды под сейсмостанцией осуществляется в терминах плоскослоистой модели. Синтетические сейсмограммы рассчитываются методом Томсона-Хаскелла. По ряду причин невозможно определить единственную модель, наилучшим образом объясняющая данные наблюдений. Для снижения степени неопределённости решения вместе с ПФ в инверсии используются другие сейсмические данные: дисперсионные кривые поверхностных волн, невязки времен пробега телесейсмических фаз и др. Поэтому для оценки сейсмической модели под станцией используется вероятностный подхода. Центральную роль здесь играет плотность апостериорного распределения вероятности и теорема Байеса. Апостериорное распределение позволяет не только определить наиболее вероятные значения параметров модели, но также их доверительный интервал. В настоящий момент наиболее распространён подход к построению апостериорного распределения, основанный на аппарате цепей Маркова. Он имеет ряд недостатков, главные из которых связано с ограничением на количество параметров модели, сложностью настройки и пр. В нашем подходе построение апостериорного распределения осуществляется на основе комбинации методов Монте-Карло и квазиньютоновского градиентного метода. К достоинствам метода следует отнести минимальное количество подстроечных параметров и слабая зависимость решения от параметризации. Разработанный метод был опробован как на синтетических данных, так и на реальных наблюдениях