Модификации свд: СНАЙПЕРСКИЕ ВИНТОВКИ СВД, СВДС И СВУ СИСТЕМЫ ДРАГУНОВА

Содержание

СВД - Warface

Снайперская Винтовка Драгунова

СВД
Снайперская винтовка
Общая информация
СтоимостьБесплатно
КатегорияСтартовое (Обычное)
Оружие доступно классам
ТТХ
Урон108
Дальность60
Темп стрельбы370
Точность "в прицел"104
Точность "с бедра"57
Емкость магазина10
Достижения

Некуда бежать

Убить 10000 врагов из 'СВД'

Доступные модификации
Ствол

Снайперская винтовка системы Евгения Драгунова (СВД), принятая на вооружение Советской Армии в 1963 году, представляет собой самозарядное оружие с газоотводной автоматикой. Автоматика винтовки действует за счет отвода пороховых газов через боковое отверстие в стенке ствола. Специально для СВД был создан "снайперский" патрон с пулей со стальным сердечником, однако винтовка может использовать все типы российских патронов 7.62х54R. Все варианты винтовки Драгунова оснащены несъемными открытыми прицельными приспособлениями в виде мушки в намушнике и регулируемого целика, расположенного перед крышкой ствольной коробки. Кронштейн типа "Ласточкин хвост" для оптического прицела крепится к ствольной коробке слева. Изначально цевье и приклад рамочной конструкции были выполнены из дерева, на более современных вариантах эти детали выполняются из пластика. Одна из характерных особенностей СВД - наличие на стволе прилива для установки штык-ножа.

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

Особенности

СВД - входит в состав начального вооружения снайпера. Предоставляется игроку бесплатно для ознакомления с классом Снайпер. У винтовки практически отсутствует увод ствола в сторону после выстрела. Движение происходит только в вертикальной плоскости, по "оси роста", это даёт возможность снайперу делать два последовательных выстрела с минимальной паузой. Например, произведя первый выстрел в грудь, второй выстрел пойдет в район шеи и головы противника - последствия для него будут фатальными. Несмотря на некоторые недостатки, в умелых руках СВД является достаточно опасным оружием.

Плюсы:

  • Не требует починки.
  • Выдаётся бесплатно.
  • Вместительный магазин.
  • Маленькая отдача.
  • Возможность установки глушителя.
  • Хорошая скорострельность.
  • Хорошая точность стрельбы "от бедра".

Но есть у этой винтовки и минусы:

  • 1 слот модификации.
  • Невысокий урон.
  • Малая дальность поражения цели.
Ствол

СВД - Снайперская Винтовка Драгунова в Warface

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

  • Цена — бесплатно
  • Урон — 150
  • Дальность — 60
  • Темп стрельбы — 250
  • Точность в прицеле — 104
  • Точность от бедра — 63
  • Ёмкость магазина — 12
  • Боезапас — 60
  • Отдача — Слабая
  • Модификации — ствол, цевьё, прицел
  • Принцип стрельбы — одиночными без передёргивания затвора
  • Способ перезарядки — смена магазина

[adspost]

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

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

Лучше всего применять на самом легком уровне; или же весьма опытным игрокам. Одной из немногих положительных характеристик СВД можно считать то, что она практически не нуждается в ремонте. Бонусный модуль только один — вы можете оборудовать глушитель на ствол винтовки.

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

Модификации СВД

СВД Сталь
(День ВДВ)
Урон — 170 (+20)
Дальность — 70 (+10)

СВД Джунгли
(Подарок другу)
Урон — 190 (+40)
Отдача — средняя

Видео обзор СВД

Снайперская винтовка Драгунова

Калибр7,62 мм
Начальная скорость пули830 м/с
Максимальная дальность стрельбы1300 м
Скорострельность30 в/мин

7,62-мм снайперская винтовка Драгунова

СВД, Индекс ГРАУ — 6В1

Самозарядная снайперская винтовка, созданная в 1957—1963 годах группой конструкторов под руководством Евгения Драгунова.

История

В середине 1960-х годов в технические процессы производства винтовки СВД были внесены изменения: оружейники И. А. Самойлов и В. Никитин создали новую технологию изготовления ствола.

В 1990-е годы винтовку начали комплектовать пластмассовым цевьем. Кроме того, был освоен выпуск конверсионного варианта винтовки — самозарядного карабина «Тигр» (конструктивно отличается более коротким стволом, отсутствием пламегасителя, газового регулятора и прилива для крепления штыка, измененной фурнитурой).

Боеприпасы и комплектация

Для стрельбы из СВД используются винтовочные патроны 7,62х54 мм R с обыкновенными, трассирующими и бронебойно-зажигательными пулями, а также снайперские патроны (7Н1, 7Н14), способна также стрелять патронами с экспансивными пулями JHP и JSP.

Огонь из СВД ведется одиночными выстрелами. Подача патронов при стрельбе осуществляется из коробчатого магазина ёмкостью 10 патронов. На дульной части ствола крепится пламегаситель с пятью продольными прорезями, маскирующий также выстрел в ходе ночных операций и предохраняющий от загрязнения ствола. Наличие газового регулятора для изменения скоростей отката подвижных частей обеспечивает надёжность винтовки в работе.

Винтовка комплектуется оптическим прицелом ПСО-1М2, есть возможность установки ночных прицелов НСПУМ или НСПУ-3.

Принцип действия

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

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

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

Для производства очередного выстрела необходимо отпустить спусковой крючок и нажать на него снова. После освобождения спускового крючка тяга продвигается вперед и её зацеп заскакивает за шептало, а при нажатии на спусковой крючок зацеп тяги поворачивает шептало и разъединяет его с боевым взводом курка. Курок, поворачиваясь на своей оси под действием боевой пружины, наносит удар по ударнику, а последний продвигается вперед и производит накол капсюля-воспламенителя патрона. Происходит выстрел.

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

Точность и кучность

в соответствии с «Наставлением по стрелковом делу» кучность боя винтовки считается нормальной, если при стрельбе четырьмя выстрелами из положения лёжа на дальность 100 м все четыре пробоины умещаются в круг диаметром 8 см. Стрельба ведётся взятыми из герметичной укупорки патронами из одной партии с обыкновенными пулями со стальным сердечником.

Изначально СВД производилась с шагом нарезов 320 мм, аналогичным спортивному оружию и обеспечивающим высокую кучность стрельбы, однако при принятии винтовки на вооружение обнаружилось, что с таким шагом рассеивание бронебойно-зажигательных пуль Б-32 увеличивается в два раза. В результате было решено изменить шаг нарезов на 240 мм, что привело к увеличению норматива по рассеиванию на дальности 100 м с 8 см до 10 см (что, однако, не нашло отражения в НСД). При стрельбе снайперским патроном 7Н1 рассеивание составляет (в зависимости от шага нарезов) не более 10—12 см на дистанции 300 м.

Дальность прямого выстрела составляет:

-по головной фигуре, высотой 30 см — 350 м,
-по грудной фигуре, высотой 50 см — 430 м,
-по бегущей фигуре, высотой 150 см — 640 м.
Прицел ПСО-1 рассчитан на стрельбу до 1300 метров, однако на такой дальности можно эффективно стрелять только по групповой цели, либо вести беспокоящий огонь.

Главной сложностью при стрельбе на большие дальности являются ошибки подготовки исходных данных для стрельбы (это справедливо для всех снайперских винтовок). На дальности 600 метров срединная ошибка по высоте (в определении дальности, равной 0,1 % дальности) — 63 см, срединная ошибка в боковом направлении (определение скорости бокового ветра, равной 1,5 м/с) — 43 см. Для сравнения, срединное отклонение рассеивания пуль для лучших снайперов для 600 м — по высоте 9,4 см, боковое 8,8 см.

Хорошая подготовка снайпера позволяет вести прицельный огонь даже по таким целям, как вертолёты и малоскоростные самолёты. Например, известен случай, когда бойцу партизанского отряда ФНОФМ удалось сбить выстрелом из СВД реактивный штурмовик ВВС Сальвадора. Это произошло 12-го ноября 1989 года у деревни Сан-Мигель. Заходивший в атаку самолёт Cessna A-37B удачно вписался в прицел и был поражён (позже удачливый снайпер рассказывал, что целился в кабину экипажа). Пуля попала в лётчика, после чего самолёт потерял управление и разбился. Похожим образом применяли СВД иракские боевики, заявлявшие об уничтожении огнём снайперских винтовок малых разведывательных БПЛА RQ-11 Raven.

Модификации

-СВДС

Вариант СВД со складным прикладом и укороченным стволом.

-СВДК

Крупнокалиберный вариант СВД под патрон 9,3х64 мм со складным прикладом, такому же как и у СВДС.

-СВУ

Вариант СВД с компоновкой булл-пап.

-ТСВ-1

Тренировочная винтовка под патрон .22 Long Rifle, созданная Евгением Драгуновым для первоначального обучения снайперов. По сути это самостоятельное оружие, лишь повторяющее в общих чертах внешний вид СВД.

Страны-эксплуатанты

-СССР
-Россия
-Азербайджан


-Армения
-Албания
-Афганистан
-Белоруссия
-Болгария
-Вьетнам

-Венгрия
-Венесуэла: приобретены для вооружённых сил Венесуэлы
-Грузия
-Индия: производится по лицензии.
-Ирак: Al-Qadissiya или Al-Gadissiya — вариант местного производства.


-Иран: Nakhjir — вариант местного производства
-Казахстан
-Киргизия
-КНР: Тип 79, клон производства Norinco, а также модернизированная версия Тип 85 и несколько вариантов для гражданского рынка.
-Никарагуа
-Монголия
-Польша
-Румыния: производится по лицензии.
-Словакия
-Таджикистан
-Туркмения
-Турция: эксплуатируется армейской жандармерией.
-Узбекистан
-Украина


-Финляндия: под обозначением 7.62 TKIV Dragunov.
-Чехия: на вооружении армии

ТТХ

-Масса, кг:
-4,3 (СВД, ранний выпуск, без штыка-ножа, с оптическим прицелом, неснаряжённым магазином и щекой приклада)
-4,5 (СВД, современный вариант, без штыка-ножа, с оптическим прицелом, неснаряжённым магазином и щекой приклада)
-4,68 (СВДС с оптическим прицелом и неснаряженным магазином)
-0,21 (магазин)
-0,26 (штык-нож без ножен)
-0,58 (прицел ПСО-1)


-Длина, мм:
-1225 (СВД без штыка-ножа)
-1370 (СВД со штыком-ножом)
-1135/875 (СВДС с разложенным/сложенным прикладом)
-Длина ствола, мм:
-620 (СВД, общая)
-547 (СВД, нарезная часть)
-565 (СВДС)
-Ширина, мм: 88
-Высота, мм: 230
-Патрон: 7,62х54 мм R
-Калибр, мм: 7,62
-Принципы работы: Поворотный затвор, отвод пороховых газов
-Скорострельность, выстрелов/мин: 30 (боевая)


-Начальная скорость пули, м/с:
-830 (СВД)
-810 (СВДС)
-Прицельная дальность, м:
-1200 (открытый прицел)
-1300 (оптический прицел)
-300 (ночные прицелы НСПУМ и НСПУ-3)
-Максимальная дальность, м:
-800 (эффективная)
-3800 (убойное действие пули)
-Вид боепитания: коробчатый магазин на 10 патронов
-Прицел: открытый секторный (резервный), длина прицельной линии — 587 мм, имеется крепление для установки оптических (например, ПСО-1) или ночных (например, НСПУ-3 или НСПУМ) прицелов

Калибр7,62 мм
Начальная скорость пули830 м/с
Максимальная дальность стрельбы1300 м
Скорострельность30 в/мин
Снайперская винтовка

СВД v1_5 в Fallout New Vegas

Об этом моде

Добавляет стандартную и уникальную версию снайперской винтовки СВД с 3 модами.

Требования

У этого мода нет никаких известных зависимостей, кроме базовой игры.

Разрешения и кредиты

Кредиты и разрешение на распространение

  • Активы других пользователей Все ресурсы в этом файле принадлежат автору или взяты из бесплатных ресурсов моддера.
  • Разрешение на загрузку Вы можете загрузить этот файл на другие сайты, но вы должны указать меня как создателя файла.
  • Разрешение на изменение Вы должны получить от меня разрешение, прежде чем вам будет разрешено изменять мои файлы для улучшения
  • Разрешение на преобразование Вам не разрешается конвертировать этот файл для работы в других играх ни при каких обстоятельствах.
  • Разрешение на использование активов Вы должны получить от меня разрешение, прежде чем вам будет разрешено использовать какие-либо активы в этом файле.
  • Разрешение на использование активов в проданных модах / файлах Вам не разрешается использовать активы из этого файла в каких-либо модах / файлах, которые продаются за деньги в Steam Workshop или на других платформах.
  • Разрешение на использование активов в модах / файлах, зарабатывающих баллы пожертвования Вы должны получить разрешение на зарабатывание очков пожертвования за свои моды, если они используют мои активы

Заметки автора

Если вы хотите использовать их в личных целях, прекрасно.Если вы хотите распространять их, использовать в коммерческих целях и т. Д., Пожалуйста, сначала спросите меня. Спасибо.

Авторские права на файл

Этот автор никого не указывал в этом файле

Система баллов пожертвований

Этот мод не включен для получения очков пожертвования

********************************
* Винтовка СВД для Fallout: New Vegas *
****** **************************

ПОМНИТЕ, что вы используете этот мод на свой страх и риск, я не несу ответственности за любой ущерб вашему компьютеру, игра, себя и т. д.
Также модели, включенные в этот мод, довольно высокополигональны, поэтому вы можете заметить некоторое падение производительности на более медленных машинах.

Этот мод добавляет две модели российской снайперской винтовки СВД. Вы можете купить их обоих у Gunrunners.
Включены модели:
Обычная СВД: "Винтовка Драгунова СВД"
Уникальная СВД: "Винтовка Драгунова нестандартная"
Они стреляют патронами калибра .308 и немного мощнее снайперской винтовки, но они тяжелее и дороже. Более того, обычный SVD можно модифицировать. Для него есть 3 модификации:
Синтетический приклад и крышка ствола
Высокомощный прицел 1p59
Глушитель (полностью собранный, но его не обязательно должно быть в RL, не так ли?)

****** ********
* Установка *
**************
Просто распакуйте архив в папку с игрой.

------------------------------------------------ ----------------------------
Для идей, комментариев, предложений и т. Д. Пишите мне:
[адрес электронной почты защищен]

Ковариация модификаций подпространственных баз (технический отчет)

Харрис Д. Б. Ковариантные модификации основ подпространства . США: Н. П., 2008. Интернет. DOI: 10,2172 / 945871.

Харрис, Д. Б. Ковариационные модификации основ подпространства . Соединенные Штаты. https://doi.org/10.2172/945871

Харрис, Д. Б. Ср. «Ковариантные модификации основ подпространства». Соединенные Штаты. https://doi.org/10.2172/945871. https://www.osti.gov/servlets/purl/945871.

@article {osti_945871,
title = {Ковариационные модификации основ подпространства},
author = {Харрис, Д. Б.},
abstractNote = {Адаптивные алгоритмы обработки сигналов, которые полагаются на представления подпространств сигнала и шума, часто требуют обновления этих представлений, когда становятся доступными новые данные.Представления подпространств часто оцениваются на основе имеющихся данных с помощью разложения по сингулярным значениям (SVD). Обновления подпространства требуют модификации этих декомпозиций. Обновления могут быть выполнены недорого, если они невысоки. Существует обширная литература по обновлениям SVD, часто фокусирующаяся на обновлениях ранга 1 (см., Например, [Karasalo, 1986; Comon and Golub, 1990, Badeau, 2004]). В этих методах матрицы данных изменяются путем добавления или удаления строки или столбца, или матрицы ковариации данных изменяются путем добавления внешнего произведения нового вектора.В недавней статье Brand [2006] представлен общий и эффективный метод произвольного обновления ранга SVD. Цель этой заметки - описать тесно связанный метод для приложений, где не требуются правые сингулярные векторы. В этом примечании также описываются обновления SVD для конкретного сценария обработки сигналов сейсмической группы, представляющего интерес. Конкретное приложение включает обновление представления широкополосного подпространства, используемого в сейсмических детекторах подпространства [Harris, 2006]. Эти детекторы подпространства обобщают алгоритмы корреляции формы волны для обнаружения сигналов, которые лежат в подпространстве форм волны размерности d {ge} 1.Они потенциально представляют интерес, поскольку расширяют диапазон изменения формы сигнала, в котором применяются эти чувствительные детекторы. Детекторы подпространства работают, проецируя данные формы сигнала из окна обнаружения в подпространство, заданное набором базисных векторов ортонормированной формы сигнала (называемого шаблоном). Шаблоны подпространства состоят из набора нормализованных, выровненных осциллограмм основных событий, которые могут быть получены одним датчиком, трехкомпонентным датчиком, массивом таких датчиков или сетью датчиков. Процесс проектирования шаблона влечет за собой построение матрицы данных, столбцы которой содержат данные формы сигнала основного события, а затем выполнение разложения по сингулярным значениям в матрице данных для извлечения ортонормированной основы для набора сигналов. Шаблон обычно состоит из подмножества левых сингулярных векторов, соответствующих большим сингулярным значениям. Приложение включает в себя обновление шаблона подпространства, когда становятся доступными новые данные, то есть когда обнаруживаются новые определяющие события для конкретного источника.Часто бывает так, что формы сигналов, соответствующие конкретному источнику, дрейфуют во времени [Harris, 2001]. Функции Грина, описывающие распространение, могут быть изменены из-за изменений в области источника. Например, если источником является шахта, сигналы от взрывов могут постепенно изменяться по мере расширения ямы (смещения источника) или изменения топографии рассеяния в результате раскопок. Это мотивирует корректировку отслеживания для представления подпространства. В этой заметке также комментируются обновления SVD для связанной проблемы.В реальных конвейерных операциях часто бывает так, что данные из одного или нескольких каналов массива непригодны для использования (мертвые каналы, каналы с частыми выпадениями и т. Д.). В таких случаях желательно изменить шаблон подпространства массива для работы с данными только из оставшихся используемых каналов. Кроме того, желательно изменять шаблоны напрямую, не обращаясь к исходной матрице данных. Обычно процесс разработки шаблона отделен от применения шаблона в детекторе к непрерывному потоку данных.Следовательно, исходная матрица данных может быть недоступна для модификации шаблона во время работы детектора.},
doi = {10.2172 / 945871},
url = {https://www.osti.gov/biblio/945871}, journal = {},
number =,
объем =,
place = {United States},
год = {2008},
месяц = ​​{11}
}

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

Оценщик на основе SVD Мы переформулируем этот алгоритм на основе SVD ниже. \ top, \) где \ ({\ tilde {K}} = \ max \ big \ {K + 1, \ mathop {{\ text {arg max}}} \ limits _k \ {\ sigma _k \ ge 1.01 \ sqrt {N} \} \ big \} \).

  • 4.

    Пусть \ ({\ hat {X}} = ({\ hat {x}} _ {ij}) _ {N \ times J} \) определено как

    $$ \ begin {align} {\ hat { x}} _ {ij} = {\ left \ {\ begin {array} {ll} \ epsilon _ {N, J}, \ quad \ text {if} x_ {ij} <\ epsilon _ {N, J} , \\ x_ {ij}, \ quad \ text {if} \ epsilon _ {N, J} \ le x_ {ij} \ le 1- \ epsilon _ {N, J}, \\ 1- \ epsilon _ { N, J}, \ quad \ text {if} x_ {ij}> 1 - \ epsilon _ {N, J}.\ top \), где \ ({{\ hat {\ sigma}}} _ 1 \ ge \ cdots \ ge {{\ hat {\ sigma}}} _ J \ ge 0 \) - сингулярные значения, а \ ({ {\ hat {{\ mathbf {u}}}}} _ j \) s и \ ({{\ hat {{\ mathbf {v}}}}} _ j \) s являются левым и правым сингулярными векторами соответственно.

  • 8.

    Выходные данные \ ({\ hat {A}} = \ frac {1} {\ sqrt {N}} ({{\ hat {\ sigma}}}} _ 1 {\ hat {{\ mathbf {v}}}} _ 1 , \ ldots, {\ hat {\ sigma}} _ K {\ hat {{\ mathbf {v}}}} _ K), {\ hat {\ Theta}} = \ sqrt {N} ({\ hat {{\ mathbf {u}}}} _ 1, \ ldots, {\ hat {{\ mathbf {u}}}} _ K). \)

  • Примечание 1

    SVD - это мощный инструмент факторизации прямоугольных матриц, который широко используется в многомерной статистике для уменьшения размерности данных (Wall et al.2003 г.). Благодаря математическим свойствам SVD оценка, данная алгоритмом 1, является аналитической и не имеет проблем сходимости. С другой стороны, поскольку целевые функции CJMLE и маргинальной оценки максимального правдоподобия (MMLE; Bock and Aitkin 1981) невыпуклые, нет никакой гарантии нахождения их глобального оптимума. Кроме того, этот подход SVD также намного быстрее, чем другие средства оценки, включая CJMLE и MMLE. В частности, вычисление MMLE на основе алгоритма максимизации ванильного ожидания недоступно, когда скрытый размер K имеет средний размер (например,g., \ (K \ ge 5 \)). Даже стохастические алгоритмы для MMLE (Cai 2010a; 2010b; Zhang et al. 2020) и алгоритм чередующейся минимизации для CJMLE (Chen et al. 2019b; 2019c) намного медленнее, чем алгоритм SVD, поскольку для этих алгоритмов обычно требуется большое количество итераций для схождения. Сравнение скорости приведено в исследовании моделирования между методом SVD и CJMLE.

    Замечание 2

    Алгоритм 1 можно рассматривать как обобщение PCA на двоичные данные.PCA - это алгоритм на основе SVD (например, глава 14, Friedman et al. 2001), который является быстрым и обычно используется для исследовательского линейного факторного анализа. К сожалению, PCA не может применяться к исследовательской IFA из-за нелинейной функции связи в моделях IFA. В отличие от PCA, который применяет SVD только один раз, алгоритм 1 применяет SVD дважды. Первое применение SVD и обратного преобразования (шаги 2–5) снижает шум и линеаризует данные. Затем второе приложение SVD (шаги 6–7) по существу выполняет PCA для линеаризованных данных.

    Замечание 3

    Подобно CJMLE (Chen et al. 2019b; 2019c), эта оценка на основе SVD не требует, чтобы латентное распределение F было известно или принимало параметрическую форму, как требуется в подходе MMLE. Более того, исследовательский IFA, основанный на тетрахорических / полихорических корреляциях (Muthén 1984; Lee et al. 1990; Lee et al. 1992; Jöreskog 1994) или составной вероятностной оценке (Katsikatsou et al. 2012) требует, чтобы F был многомерным. нормально, при первом подходе дополнительно требуется, чтобы обратная линия f была пробит.В этом смысле для оценки на основе SVD и CJMLE требуется меньше допущений модели, чем для других оценок. В качестве цены их согласованность требует более строгих условий, а именно режима двойной асимптотики, при котором расходятся как N , так и J .

    Замечание 4

    Шаги 2–4 алгоритма по существу следуют той же процедуре, что и у Чаттерджи. 2} \) обозначает норму Фробениуса матрицы \ (X = (x_ {ij}) \).*, {\ hat {A}}) \) сходится к нулю по вероятности. Условия регулярности и результат согласованности формально описаны в теореме 1, а два частных случая обсуждаются в дальнейшем. Подобные двойные асимптотические параметры рассматривались в психометрических исследованиях, включая анализ одномерных моделей IRT (Haberman 1977; 2004) и диагностических классификационных моделей (Chiu et al. 2016). Следующие условия регулярности необходимы для нашего основного результата в теореме 1. Как будет показано ниже, эти условия мягкие.* \) независимы и одинаково распределены (i.i.d.) в соответствии с распределением F , которое имеет среднее значение \ ({\ mathbf {0}} \) и положительно определенную матрицу ковариации \ (\ Sigma. \)

  • A3.

    Функция обратной связи f является строго монотонно возрастающей, непрерывно дифференцируемой и непрерывной по Липшицу с постоянной Липшица L . Мы также предполагаем, что

    $$ \ begin {выровнено} \ lim _ {x \ rightarrow - \ infty} f (x) = 0, ~~ \ text {и} ~~ \ lim _ {x \ rightarrow \ infty} е (х) = 1.*, {\ hat {A}}) \ overset {pr} {\ rightarrow} 0 \), как \ (N, J \ rightarrow \ infty. \)

    Замечание 6

    Отметим, что понятие согласованности для оценки матрицы нагрузки слабее, чем в традиционном смысле, поскольку функция потерь (2) представляет собой среднее значение потерь на входе при росте J . Пусть \ ({\ tilde {O}} \) минимизирует правую часть (2), и пусть \ ({\ tilde {A}}: = ({\ tilde {a}} _ {jk}) _ {J \ times K} = {\ hat {A}} {\ tilde {O}} \).* - {\ tilde {a}} _ {jk} \ vert> \ epsilon \}}}) / {JK} \) также сходится к 0. То есть доля неточно оцененных параметров нагрузки сходится к нулю по вероятности при оптимальное вращение. Из-за двойной асимптотики наш теоретический результат предполагает разумное использование алгоритма на основе SVD только тогда, когда размер выборки N и количество элементов J являются большими.

    Замечание 7

    Хорошо известно, что PCA может последовательно оценивать линейную факторную модель при аналогичных условиях двойной асимптотики (Stock and Watson 2002), что обеспечивает теоретическое обоснование использования PCA в исследовательском линейном факторном анализе.Теорему 1 можно рассматривать как аналогичный результат для исследовательского факторного анализа.

    Замечание 8

    Мы предлагаем некоторые обсуждения условий регулярности, требуемых в теореме 1. Предположение A1 требует, чтобы параметры каждого элемента, включая параметры точки пересечения и наклона, не были слишком большими. То есть наличие экстремального элемента может исказить анализ. Предположение A2 - очень стандартное предположение для исследовательской IFA.Он более гибкий, чем многие предварительные настройки IFA, поскольку не требует, чтобы распределение F было многомерным нормальным.

    Предположение A3 удовлетворяется функциями логистической и пробит-ссылки, двумя наиболее часто используемыми функциями связи в исследовательском IFA, но оно исключает, например, многомерную версию трехпараметрической логистической модели как особый случай. Предположение A4 требует, чтобы элементы были достаточно разнообразными. Такое же предположение требуется и у Chen et al.* \) являются i.i.d. выборки из распределения K -вариант, ковариационная матрица которого невырождена. Наконец, предположение A5 практически разумно, поскольку при крупномасштабных измерениях размер выборки обычно больше, чем количество элементов. Поскольку люди и предметы почти математически симметричны в модели IFA, аналогичные асимптотические результаты могут быть получены при \ (J \ ge N \).

    Замечание 9

    Мы также даем некоторые интуитивные предположения о том, почему алгоритм работает.*, {\ hat {A}}) \) с большой вероятностью мала.

    Замечание 10

    Уравнения (3) и (4) являются требованиями к параметру усечения \ (\ epsilon _ {N, J} \), который зависит как от хвоста распределения F , так и от свойств функции обратной связи. Грубо говоря, уравнение (3) говорит, что \ (\ epsilon _ {N, J} \) не может быть слишком большим. Это потому, что для F и f вероятность в (3) увеличивается в \ (\ epsilon _ {N, J} \).* \) больше, чем \ (h (\ epsilon _ {N, J}) \). Функция \ (h (\ cdot) \) преобразует усечение на \ (x_ {ij} \) в усечение на \ ({\ tilde {m}} _ {ij} \). Использование \ (h (2 \ epsilon _ {N, J}) \) вместо \ (h (\ epsilon _ {N, J}) \) по техническим причинам.

    Уравнение (4) требует, чтобы \ (\ epsilon _ {N, J} \) не был слишком маленьким, так как левая часть (4) убывает в \ (\ epsilon _ {N, J} \). Это требование тоже интуитивно понятно. Обратите внимание, что \ (| {\ tilde {m}} _ {ij} | \ le h (\ epsilon _ {N, J}) \), где \ (h (\ epsilon _ {N, J}) \) - убывает по \ (\ epsilon _ {N, J} \).*) | \) мала из-за локальной плоскости функции обратной связи.

    Замечание 11

    Мы берем стохастический план для истинных параметров человека и фиксированный план для истинных параметров предмета, следуя соглашению о факторном анализе предмета (например, Bartholomew et al. 2008). Стоит отметить, что выбор стохастического или фиксированного плана не является существенным в нашем режиме двойной асимптотики. Например, непротиворечивый результат теоремы 1 остается в силе, если мы можем заменить условие A2 соответствующим фиксированным планом, как в Chen et al.(2019b).

    После обсуждения \ (\ epsilon _ {N, J} \) в Замечании 10 мы рассматриваем две конкретные настройки, при которых требования к \ (\ epsilon _ {N, J} \) становятся более конкретными. Эти результаты приведены в предложениях 1 и 2.

    Предложение 1

    Предположим, что F имеет компактную опору. {\ dagger} \) в алгоритме 1 установлена ​​фиксированной (т.е., независимо от N и J ), но не обязательно равного истинному количеству факторов, существует константа \ (\ delta> 0 \) такая, что для истинного количества факторов K ,

    $$ \ begin {align} \ lim _ {N, J \ rightarrow \ infty} \ Pr \ left (\ frac {{\ hat {\ sigma}} _ {K}} {\ sqrt {NJ}}> \ delta \ right) = 1, \ text {и} \ frac {{\ hat {\ sigma}} _ {K + 1}} {\ sqrt {NJ}} \ overset {pr} {\ rightarrow} 0, \ end { выровнено} $$

    как N и J одновременно растут до бесконечности.{\ dagger} \) не влияет на асимптотику, пока она не растет с N и J . Однако для относительно небольших N и J , X , полученных на шаге 3 алгоритма, может не зарезервировать достаточно информации, когда входной размер меньше \ (K + 1 \), что может привести к недооценке количество факторов. Таким образом, в практических приложениях мы рекомендуем выбирать размер входных данных немного больше, чем максимальное количество факторов, которые, как предполагается, существуют в данных.

    Рис. 1

    Осыпной график для выбора количества факторов. На оси y показаны стандартизированные сингулярные значения \ (\ hat {\ sigma} _k / \ sqrt {NJ} \), где \ ({\ hat {\ sigma}} _ k \) s получены на шаге 7 из Алгоритм 1. Данные моделируются из модели IFA с \ (K = 5 \), \ (J = 200 \) и \ (N = 4000 \). Входной размер установлен равным 10 в алгоритме 1. Разрыв сингулярного значения может быть найден между 5-м и 6-м сингулярными значениями

    Статистическая эффективность Далее мы отмечаем, что вычислительное преимущество оценщика на основе SVD приходится расплачиваться.Чтобы уточнить этот момент, мы сравним его с CJMLE (Chen et al. 2019b; 2019c). CJMLE рассматривает как параметры элемента, так и скрытые факторы как фиксированные параметры и максимизирует совместную функцию правдоподобия по всем фиксированным параметрам. 2 / NJ, \), которая определяет сходимость \ ({\ hat {A}} \).{-1}), \) где \ ({\ hat {X}} _ {JML} \) обозначает CJMLE. Этот результат предполагает, что оценщик на основе SVD сходится гораздо медленнее, чем CJMLE.

    «Параллельная реализация разложения по сингулярным значениям с использованием Open», Бхушан Райрикар

    Абстрактные
    Блоки обработки графики общего назначения

    (GPGPU) обладают возможностями массовых параллельных вычислений. Низкая стоимость и простота программирования делают их популярным выбором среди других параллельных архитектур, таких как большие кластеры и ускорители, такие как программируемые вентильные матрицы (FPGA).Зрелые среды программирования для GPGPU, такие как CUDA от Nvidia и OpenCL от Khronos Group, сокращают время обучения и время разработки для программирования архитектур GPGPU. OpenCL, относительно новый отраслевой стандарт для параллельных вычислений, позволяет писать единую программу для разнородных платформ, переносимую на несколько платформ, включая GPGPU и многоядерные процессоры, с минимальными модификациями кода. Архитектура
    GPGPU успешно использовалась для ускорения многих вычислительно дорогостоящих задач, включая множество алгоритмов линейной алгебры, которые по своей природе параллельны.Разложение по сингулярным значениям (SVD) - это дорогостоящий в вычислительном отношении метод разложения матриц линейной алгебры, который имеет множество приложений, включая сжатие данных, распознавание лиц и решение системы уравнений. По мере увеличения размеров матрицы вычисление SVD становится все более трудоемким. Поскольку SVD является основной частью некоторых алгоритмов, таких как Eigenfaces (алгоритм распознавания лиц, основанный на анализе основных компонентов), общее время выполнения этих алгоритмов сильно зависит от времени выполнения SVD.Следовательно, для реализации эффективных приложений на основе SVD, например распознавания лиц в реальном времени, желательно ускорить алгоритм SVD.
    В этой работе подробно обсуждается параллельная реализация разложения по сингулярным значениям. Он использует многие базовые методы линейной алгебры, такие как умножение матрицы на вектор, векторные нормы и векторные внешние произведения. В этой работе основное внимание уделяется методам реализации, методам оптимизации (особенно для реализации GPGPU) и их влиянию на общую производительность алгоритма.Мы представляем анализ производительности этого алгоритма на графическом процессоре NVIDIA Tesla C2050 по сравнению с однопоточной последовательной реализацией, выполненной на процессоре Intel Q9450 с тактовой частотой 2,66 ГГц. Мы сообщаем об ускорении до 20 раз при параллельном вычислении SVD. Результаты, обсуждаемые в этой диссертации, демонстрируют потенциал вычислительных ресурсов, доступных с GPGPU.

    Рекомендуемое цитирование

    Райрикар, Бхушан, «Параллельная реализация разложения сингулярных значений с использованием OpenCL» (2011). Все тезисы . 1292.
    https://tigerprints.clemson.edu/all_theses/1292

    нелинейный усеченный SVD-дисперсионный анализ и анализ разрешения двумерных магнитотеллурических моделей | Международный геофизический журнал

    Сводка

    Представлен новый подход к оценке дисперсии и разрешающей способности двумерных моделей удельного электрического сопротивления, полученных на основе магнитотеллурических измерений. Основанный на схеме усеченного сингулярного разложения (TSVD) на локальном подпространстве, он частично учитывает нелинейность обратной задачи.Анализ разрешения и дисперсии TSVD выполняется для одной ячейки за раз. Выбирается порог дисперсии и определяется результирующее разрешение модели. В качестве усовершенствования по сравнению с существующими схемами вводятся нелинейные полуоси для описания нелинейной доверительной поверхности в направлениях собственных векторов модели, и они заменяют обратные сингулярные значения, входящие в стандартное выражение дисперсии модели. Дисперсия модели рассматриваемой ячейки оценивается из суммы квадратов нелинейных полуосей до заданного порога дисперсии. Это, в свою очередь, дает уровень усечения TSVD, и строка матрицы разрешения модели, принадлежащая рассматриваемой ячейке, может быть вычислена из собственных векторов модели TSVD. Информация, содержащаяся в матрице разрешения, сжимается до легко воспринимаемых величин, таких как центр разрешения и длины горизонтального и вертикального разрешения.

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

    Анализируется синтетическая модель с проводящим блоком в однородном полупространстве. Анализ TSVD для ячеек модели на верхнем крае блока и вне блока иллюстрирует, как работает процесс усечения. Обычно линейные и нелинейные полуоси почти равны до определенного числа сингулярного значения, после которого нелинейные полуоси увеличиваются намного меньше, чем линейные полуоси. Этот важный результат указывает на то, что разрешение двумерных магнитотеллурических моделей значительно лучше, чем ранее предлагалось линейными схемами для вычисления дисперсии модели и свойств разрешения.Полевой пример из района Скидига (Швеция) показывает, что распределение удельного электрического сопротивления песчано-гравийных образований, ограниченных только латерально проводящими глиняными линзами, относительно хорошо разрешено, в то время как разрешение перехода между слоем песка и гравия и слоем подвал под глиняным покровом.

    1 Введение

    Обратную задачу можно разделить на несколько частей (Parker 1977): (1) исследование существования и единственности решения, (2) построение решения, (3) анализ устойчивости проблемы и (4) оценка значимости решения.

    Здесь мы описываем новый метод оценки значимости решения с помощью разрешения модели и дисперсионного анализа в духе Бэкуса и Гилберта (1968).

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

    Было разработано несколько обратных методов для построения решения неуникальных и некорректно поставленных, то есть нестабильных задач. Юпп и Возофф (1975), а затем Лайнс и Трейтель (1984) описывают технику усечения сингулярных чисел и метод Марквардта как адекватные средства для решения некорректно поставленных обратных задач. Совсем недавно Педерсен (2004) представил одномерную обратную схему для магнитотеллурических данных, которая определяет оптимальный уровень усечения для разложения по сингулярным значениям на основе общей среднеквадратичной ошибки оцениваемой модели.Благодаря важной статистической информации, которую дает разложение по сингулярным значениям, информация о дисперсии модели и разрешении легко доступна из этих обратных схем (Jupp & Vozoff 1975; Lines & Treitel 1984; Pedersen 2004). Виггинс (1972) предлагает построить обобщенную обратную задачу для решения линейной обратной задачи так, чтобы дисперсия модели оставалась ниже заданного порога. Вследствие устрашающего времени вычисления разложения по сингулярным значениям по мере увеличения количества данных и параметров модели эти схемы не используются широко в двумерных или трехмерных обратных вычислениях.Как в первоначальном варианте двумерной схемы Оккама (de Groot-Hedlin & Constable, 1990), так и в ее формулировке подпространства данных (Siripunvaraporn & Egbert 2000) вычисляется матрица чувствительности, которая может использоваться непосредственно для оценки разрешающей способности и ковариации модели. Mackie & Madden (1993) предлагают метод релаксации сопряженного градиента, который позволяет избежать явного построения матрицы чувствительности и требует только знания о влиянии матрицы чувствительности на произведение матрица-вектор.Следовательно, использование метода сопряженных градиентов требует других методов разрешения модели и оценки ковариации (см. Ниже).

    Несмотря на то, что разрешение и ковариация модели могут быть относительно легко получены из обратных схем, которые вычисляют матрицу чувствительности и ее обобщенную обратную, в прошлом такой способ интерпретации обычно не применялся. Частично это связано с тем, что линейная аппроксимация обратной задачи не учитывает все модели, которые соответствуют данным с одним и тем же допуском.Следовательно, большинство подходов, использованных для оценки свойств разрешения и ковариации, были ограничены получением максимальной глубины исследования, относящегося к рассматриваемому набору данных, либо линейными, либо так называемыми нелинейными исследованиями, либо их комбинацией. В этом контексте нелинейные исследования включают модели, которые значительно отличаются от предпочтительной модели, полученной с помощью линеаризованного обратного процесса. Для простого случая одномерной инверсии данных LOTEM Kalscheuer et al. (2006) использовал различные ограничения гладкости схемы Оккама, чтобы получить максимальную глубину исследования из нелинейно связанных моделей.Ольденбург и Ли (1999) выступают за использование нелинейных исследований для расчета индекса глубины исследования (DOI) для наборов данных удельного сопротивления постоянному току и наведенной поляризации. После вычисления предпочтительной модели м *, которая минимизирует целевую функцию модели с учетом ограничений данных, они модифицируют целевую функцию модели для получения других моделей, которые соответствуют данным в той же степени и близки, но не линейно близки к м *. Предполагается, что эта процедура гарантирует более полное исследование модельного пространства, чем дает линейный анализ, такой как анализ Бэкуса-Гилберта (Parker 1977), который рассматривает только модели, линейно близкие к м *.В зависимости от выбранной целевой функции модели индекс DOI является либо клеточной нормализованной разницей модели, либо полученным на основе взаимной корреляции. Ячейки модели выше определенной изолинии DOI не считаются разрешенными данными.

    Schwalenberg et al. (2002) рекомендуют сочетание линейных и нелинейных исследований чувствительности. При линейном подходе вычисляются суммы по столбцам матрицы чувствительности, нормализованные соответствующими ошибками данных, то есть чувствительность всех точек данных для одной ячейки модели входит в сумму.Вдохновленный Oldenburg & Li (1999), изолиния чувствительности, которая считается подходящей, используется для исключения клеток с более низкой нормализованной чувствительностью из интерпретации. Из-за недостатка линеаризованного обратного процесса для поиска моделей, которые значительно отличаются от предполагаемой модели, Schwalenberg et al. (2002) предлагает систематические исследования нелинейной чувствительности на основе прямого моделирования, а также включение априорных моделей в прямое моделирование и обратные вычисления.Поскольку прямое моделирование не учитывает взаимозависимость между параметрами модели, использование априорной информации модели в обратных вычислениях считается более подходящим для получения информации о глубине исследования (Schwalenberg et al. 2002).

    Фридел (2003) предлагает более строгий анализ разрешения применительно к данным томографии удельного сопротивления, который не ограничивается оценкой глубины исследования. Он представляет одношаговую линеаризованную обратную процедуру, основанную на усеченном разложении по сингулярным числам.Из-за используемого приближения с низким контрастом разрешение модели, ковариация модели и оценки разрешения данных вычисляются a priori , то есть на основе SVD матрицы чувствительности для исходной модели. Friedel (2003) предлагает вычислить апостериорных оценок , если используется нелинейная обратная схема с несколькими итерационными шагами. Friedel (2003) также описывает простые меры разрешения как радиус разрешения и флаг искажения (см. Ниже). Из ковариационной матрицы модели извлекается процент шума изображения.

    Минкофф (1996) использует алгоритм сопряженного градиента для решения обратных сейсмических задач. Поскольку матрица чувствительности и, что более важно, матрица нормалей никогда явно не формируются во время инверсии, Минкофф (1996) формирует приближенную матрицу разрешения модели, оценивая собственные значения и собственные векторы матрицы нормалей с помощью метода Ланцоша во время обратной процедуры. Аналогичным образом Yao et al. (1999) выводит как приближенное разрешение модели, так и матрицы ковариации для итеративного решателя LSQR данных сейсмической томографии.

    Для электромагнитных проблем Alumbaugh & Newman (2000) вычислили разрешение модели и матрицы ковариации как для прямых, так и для итерационных обратных схем. Однако для обоих случаев разрешение и ковариация оцениваются из обобщенной обратной задачи, относящейся к регуляризованной обратной задаче Тихонова. Таким образом, подход, принятый Alumbaugh & Newman (2000), дает оценки разрешения и ковариации, которые частично отражают регуляризацию Тихонова обратной схемы.

    Ни в одном из приведенных выше анализов (Minkoff 1996; Yao et al. 1999; Alumbaugh & Newman 2000; Friedel 2003) - это компромисс между разрешением модели и ковариацией (Menke 1989) или влиянием нелинейности обратной задачи на оценку дисперсии модели и учитываемых свойств разрешения.

    Здесь мы представляем дисперсию модели и анализ разрешения моделей удельного сопротивления по двумерной обратной схеме с помощью нелинейного усеченного SVD-анализа в локальном подпространстве. Анализ проводится после последней итерации обратного процесса.Не делается никаких предположений об ограничениях гладкости или других средствах уменьшения некорректности, используемых во время обратного процесса, и анализ основан на SVD якобиана, а также на полуосях вдоль собственных векторов модели до заданного уровня нелинейная доверительная поверхность. Нелинейные полуоси ранее использовались Йохансеном (1977) и Педерсеном и Расмуссеном (1989) для оценки доверительных границ нелинейных параметров одномерных моделей удельного сопротивления, которые были получены из данных удельного сопротивления постоянному току и магнитотеллурических данных соответственно.Впервые мы применяем эту технику к двумерным моделям удельного сопротивления и, что более важно, используем дисперсию параметров нелинейной модели для определения уровня усечения SVD, для которого вычисляется ядро ​​разрешения рассматриваемой ячейки. Этот подход позволяет эффективно изучать влияние нелинейности на разрешающую способность модели, аспект геофизической обратной теории, которым в значительной степени пренебрегали. Вычисление нелинейной дисперсии до заранее определенного порога дисперсии модели дает уровень усечения SVD и в то же время фиксирует точку на кривой компромисса между разрешением модели и ковариацией.Из матрицы разрешения, соответствующей определенному таким образом уровню усечения, оцениваются характерные меры разрешения в виде центра разрешения, а также длины горизонтального и вертикального разрешения. Мы добавили наш анализ разрешения и дисперсии к известному двумерному магнитотеллурическому обратному коду REBOCC Сирипунварапорна и Эгберта (2000).

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

    2 Теория

    В обратных задачах вектор значений измерения N d = ( d 1 ,…, d N ) T с соответствующими стандартными ошибками σ = ( σ 1 ,…, σ N ) T связано с вектором M параметров модели m = ( m 1 ,…, m M ) T через обычно нелинейный прямой оператор F : 1, где обозначение a W = Wa обозначает умножение на диагональную матрицу весов W = диаг { σ −1 1 ,…, σ −1 N }, который содержит стандартные ошибки.

    Вектор ошибок e описывает качество соответствия прямого отклика F [ m ] текущей модели m к измеренным данным d . Основная стратегия улучшения подгонки - минимизировать кумулятивную нормированную квадратичную ошибку Q = e T W e W относительно параметров модели (Lines & Treitel 1984; Menke 1989), то есть требовать, чтобы градиент Q относительно м исчез, ∇ м Q = 0.

    При определении вектора несоответствия данных исходная задача минимизации нелинейной величины Q = e T W e W заменяется минимизацией линеаризованная величина Q lin = ( e lin W ) T e lin W = (Δ43 d - J W Δ м ) T d W - J W Δ44 м м м м

    С определением обобщенных параметров α p = V T p Δ m и обобщенными данными β p 3 43 U p Δ d , решение обратной задачи (уравнение 4) можно записать как 5

    2.1 Разрешение модели и ковариация

    Матрица разрешения модели R связывает модель m est , полученную в результате инверсии TSVD, с неизвестной истинной моделью m (Menke 1989): 6 Строку k -й матрицы разрешения можно рассматривать как набор коэффициентов, описывающих, с какими весами параметры истинной модели m входят в k -й параметр оценочной модели m est .Если, например, k -й элемент строки k -й является единственным ненулевым элементом, то только k -й параметр истинной модели m будет вносить вклад в k -й параметр. оценочной модели m est и k -й параметр считается полностью разрешенным. В отличие от этого, если имеется большое количество ненулевых элементов в строке k с приблизительно равными амплитудами, то параметр k -й оценочной модели m est будет средним за широкий диапазон параметров модели и, как говорят, плохо разрешается.Матрица разрешения, идентичная единичной матрице, то есть R = I , подразумевает идеальное разрешение всех параметров модели и может быть получена только для переопределенных задач наименьших квадратов (Menke 1989). с шахматной сеткой, выбранной для решения прямой задачи, как видно из отношения к ядру непрерывного разрешения или ядру усреднения r (Backus & Gilbert 1968), как определено ниже. В интегральной формулировке расчетный параметр модели м est в позиции ( y , z ) связан с непрерывной истинной моделью м на площади (Backus & Gilbert 1968; Parker 1977) через 7 При переходе к дискретному случаю, когда индекс ячейки k назначается ячейке в координате ( y , z ), а индекс ячейки i назначается ячейке в ( y ′ , z ′), ур.(7) становится 8, где Δ y i и Δ z i - ширина и высота ячейки с номером i , соответственно. Однако для ячейки с номером k из ур. (6) 9 Следовательно, элементы матрицы разрешения модели R связаны с элементами дискретизированного разрешающего ядра на 10, то есть матрица разрешения R содержит отсчеты разрешающего ядра r , масштабированные по ячейке область.

    На графиках разрешения в разделах 3 и 4 мы изображаем элементы r ki разрешающего ядра ячейки k , масштабированные по их максимальному модулю.

    Ковариационная матрица модели описывает, как ошибки данных, заданные ковариационной матрицей данных [cov d ], распространяются на ошибки оценочной модели m est (Menke 1989): 11 На практике [cov d ] ] содержит только диагональные элементы, то есть квадраты общих ошибок данных.Однако, поскольку мы нормализовали якобиан на общие ошибки данных, ковариационная матрица модели упрощается до 12, а стандартные ошибки ste параметров модели являются квадратными корнями из диагональных элементов уравнения. (12), то есть 13
    2.1.1 Итерационные вычисления
    Из-за ортогональности собственных векторов модели матрица разрешения и ковариации модели может быть вычислена итеративно с увеличением сингулярного числа p . Матрица разрешения, содержащая первые p + 1 собственных векторов модели, может быть затем выражена как матрица разрешения первых собственных векторов модели p плюс внешнее векторное произведение ( p + 1) -го собственного вектора модели: 14 Увеличение количество собственных векторов модели улучшает разрешение.Точно так же ковариацию модели можно вычислить итеративным способом: 15 Однако нас интересует только дисперсия каждой ячейки. Для номера ячейки модели k мы получаем: 16, где v k p + 1 - k -й компонент собственного вектора модели ( p + 1). Увеличение количества собственных векторов и сингулярных значений означает увеличение дисперсии.
    2.1.2 Компромисс между разрешением и ковариацией

    Согласно уравнениям (14) и (16), разрешение улучшается одновременно с увеличением дисперсии, если количество включенных собственных значений и собственных векторов увеличивается.Такое поведение указывает на компромисс между разрешением модели и ковариацией (Backus & Gilbert 1968; Parker 1977; Menke 1989). Компромисс оправдывает два разных подхода. Можно зафиксировать разрешение для рассматриваемых ячеек, а затем определить дисперсию ячеек или установить порог дисперсии для всех ячеек и вычислить свойства разрешения.

    В нашем подходе мы сохраняем отклонения параметров модели ниже заранее определенного порога. Согласно ур. (16), это соответствует вычислению уровня усечения p отсечки SVD.Для этого уровня усечения мы затем вычисляем свойства разрешения из ур. (14).

    В двумерном магнитотеллурическом инверсном коде REBOCC (Siripunvaraporn & Egbert 2000) параметры модели представляют собой логарифмы удельных сопротивлений ρ по основанию 10. Это означает, что рассматриваемая оценка дисперсии для ячейки k равна квадрат логарифмической ошибки удельного сопротивления, то есть var ( м est k ) = (Δlog 10 ρ k ) 2 .Диапазон удельных сопротивлений, который мы допускаем при анализе разрешения путем определения порога дисперсии, затем определяется как log 10 ρ ± (Δlog 10 ρ) порог . Ошибка Δlog 10 ρ может быть понята как коэффициент f на удельное сопротивление ρ, если принять Δlog 10 ρ = log 10 f . Отсюда следует log 10 ρ ± (Δlog 10 ρ) порог = log 10 ρ ± log 10 f порог = log 10 (ρ · f ± 1 обмолота ).Следовательно, значения удельного сопротивления в диапазоне ρ / f thresh ,…, f thresh · ρ, заданные с коэффициентом f thresh , допускаются в нашем анализе дисперсии и разрешения модели.

    2.2 Нелинейный анализ

    Более адекватный выбор оценок дисперсии модели и разрешения отражает нелинейность основной обратной задачи, по крайней мере, частично. Упрощенное описание нелинейной доверительной поверхности в пространстве параметров представляет собой псевдогиперэллипсоид (Johansen 1977) с полуосями, которые отражают, на каком расстоянии от оптимальной модели достигается данное несоответствие в направлении собственных векторов модели.В этом контексте мы называем модель с минимальным несоответствием для определенного уровня усечения оптимальной моделью. Псевдогиперэллипсоид может иметь разную длину полуосей в положительном и отрицательном направлении конкретного собственного вектора модели. Таким образом, он напоминает гиперэллипсоид линейной доверительной поверхности с полуосями одинаковой длины.

    Упрощенное описание нелинейной доверительной поверхности может быть получено путем определения длин фактических полуосей s ± i как в положительном, так и в отрицательном направлениях собственного вектора v i (Йохансен 1977; Педерсен 1979; Педерсен и Расмуссен 1989).Фактическая форма нелинейной доверительной поверхности в направлениях, отличных от направлений собственных векторов, не учитывается и, следовательно, аппроксимируется как псевдогиперэллипсоид. На рис. 1 изображен псевдогиперэллипсоид для простого примера двух параметров модели. Фактические полуоси s ± 1 и s ± 2 соответствуют точкам в пространстве модели, где псевдогиперэллипсоидальная поверхность (тонкая линия) пересекает истинную нелинейную доверительную поверхность ( жирная линия) для данного Δ Q .Псевдогиперэллипсоид состоит из четырех эллипсоидных секций с полуосями ( s + 1 , s - 2 ), ( s + 3 1 + 2 ), ( s - 1 , s + 2 ) и ( s - 1 , s - 2 .

    Рисунок 1.

    Иллюстрация псевдогиперэллипсоидальной поверхности для двух параметров м 1 и м 2 как приближение к истинной доверительной поверхности нелинейной задачи (изменено после Pedersen & Rasmussen 1989). Толстая красная линия показывает истинную поверхность уверенности. Псевдогиперэллипсоид описывается тонкой синей линией и построен из участков четырех эллипсов с полуосями ( s + 1 , s - 2 ), ( s + 1 , с + 2 ), ( с - 1 , с + 2 ) и ( с - 9048 с - 2 ).Оси a 1 и a 2 указывают в направлениях собственных векторов v 1 и v 2 , соответственно. Крайние точки м мин 1 и м макс 1 из м 1 на псевдогиперэллипсоиде переоценивают и недооценивают истинную изменчивость модели, соответственно.

    Рисунок 1.

    Иллюстрация псевдогиперэллипсоидальной поверхности для двух параметров м 1 и м 2 как приближение к истинной доверительной поверхности нелинейной задачи (модифицировано по Педерсену и Расмуссену. 1989).Толстая красная линия показывает истинную поверхность уверенности. Псевдогиперэллипсоид описывается тонкой синей линией и построен из участков четырех эллипсов с полуосями ( s + 1 , s - 2 ), ( s + 1 , с + 2 ), ( с - 1 , с + 2 ) и ( с - 9048 с - 2 ).Оси a 1 и a 2 указывают в направлениях собственных векторов v 1 и v 2 , соответственно. Крайние точки м мин 1 и м макс 1 из м 1 на псевдогиперэллипсоиде переоценивают и недооценивают истинную изменчивость модели, соответственно.

    Задача вычисления фактических полуосей может быть сформулирована как нахождение нулей 22 для всех собственных векторов v i и m = m 0 ± s ± i v i .Требуется указать желаемое отклонение Δ Q от минимального несоответствия Q 0 . Уравнение (22) является нелинейным уравнением с одним неизвестным (либо + s + i или - s - i ) и решается с помощью безопасного метода, использующего Метод секущих как быстро сходящийся метод, но сохраняет скобку метода деления пополам вокруг решения на тот случай, если следующее приближенное решение метода секанса выпадает за пределы интервала заключения в скобки (Heath 2002).Для собственного вектора v i с i > 1, начальные скобки метода деления пополам находятся путем вычисления уравнения. (22) при ± с ± i = ± с ± i −1 с последующим увеличением или уменьшением ± с ± i постоянный коэффициент, пока не изменит знак. Для первого собственного вектора v 1 поиск исходных скобок начинается с ± s lin 1 .

    Из-за нелинейности длины нелинейных и линейных полуосей могут быть разными, особенно для самых маленьких сингулярных значений (Johansen 1977; Pedersen & Rasmussen 1989).

    После определения линейных / нелинейных полуосей доверительные границы параметров модели могут быть вычислены путем нахождения их крайних точек на поверхности (псевдо) гиперэллипсоида, который характеризуется особым выбором Δ Q . Крайние точки параметра м k находятся путем требования, чтобы градиент Δ Q был параллелен нормали n гиперплоскости м k = const (Johansen 1977 ; Pedersen 1979), поскольку гиперплоскость, соответствующая n , касается (псевдо) гиперэллипсоидальной поверхности в крайних точках (рис.1). Компонент k -й компонент n устанавливается равным +1 для максимальной точки или -1 для минимальной точки, тогда как все остальные компоненты устанавливаются равными нулю. Кроме того, Δ м должен находиться на (псевдо) гиперэллипсоидальной поверхности, характеризуемой Δ Q . Для линейного случая это дает (Johansen 1977; Pedersen 1979; Pedersen & Rasmussen 1989) 23. Для выбора Δ Q = 1, ур. (23) упрощается до стандартной ошибки параметров модели, как в ур.(13) и соответствует уровню уверенности 68%. Вычисление нелинейных экстремальных точек требует, чтобы s ± i были заменены на s lin i в уравнении. (23) согласно следующему правилу (Johansen 1977). Учитывая, что необходимо определить максимальную точку м k на псевдогиперэллипсоиде, необходимо заменить s + i на s lin 90 if ( V T n ) i > 0 и s - i для s lin i 9047 9047 9047
    T n ) i V T n ) i = n · v
    905 905 905 v
    i вдоль n , полуоси s + i положительного направления v i должны быть заменены будет выполнено, если n · v i > 0 и наоборот.Аналогичное правило должно применяться для минимальной точки м k .

    Только в том случае, если оценочная модель m est равна оптимальной модели m 0 , можно ожидать, что нелинейные и линейные оценки полуосей будут иметь равные значения для наибольших сингулярных значений. . Для широко используемых сглаженных обратных операторов обратный процесс не попадает в оптимальную модель Q и м est смещен от м 0 .Следовательно, также s + i и s - i будут систематически сдвинуты с s lin i таким образом s + i < s lin i затем s - i > s Следовательно, нелинейные эффекты невозможно отличить от смещения от оптимальной модели м 0 . Поэтому в нашем подходе мы используем комбинацию инверсии среднеквадратичной ошибки (MSE) (Pedersen, 2004) и инверсии TSVD (Jupp & Vozoff, 1975), чтобы перейти к оптимальной модели до того, как будет проведен окончательный анализ дисперсии модели и разрешения. Во время построения обратного решения неизвестно, сколько пар нелинейных полуосей должно быть включено в следующий нелинейный анализ, чтобы достичь допустимого порога дисперсии.По этой причине мы выполняем начальную инверсию с помощью метода MSE, который обычно включает меньше собственных векторов, чем требуется для следующего начального нелинейного анализа. Включение большего количества собственных векторов в нелинейный анализ означает, что включены определенные собственные векторы, которые считались частью нулевого пространства (Menke 1989) во время предыдущей инверсии MSE. Следовательно, окончательная инверсия TSVD выполняется с уровнем усечения, равным уровню исходного нелинейного анализа, а после этого выполняется окончательный нелинейный анализ.Как правило, предполагаемый уровень усечения не сильно меняется от начального к окончательному нелинейному анализу, и в большинстве случаев достаточно двух итераций улучшения модели и последующего нелинейного анализа.

    2.3 Инверсия большинства квадратов

    В нашем нелинейном анализе мы исследуем нелинейность только в направлениях собственных векторов модели. Достоверность наших оценок дисперсии можно проверить с помощью инверсии большинства квадратов (Jackson 1976; Meju & ​​Hutton 1992), которая вычисляет экстремальные сопротивления ячейки для заданного порога несоответствия и уровня усечения.

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

    Цель инверсии наибольших квадратов - найти крайние значения функционала m T n при ограничении, что несоответствие Q не превышает некоторого порога Q t .Если n выбрано в качестве нормали гиперплоскости м k = const, k -й параметр модели будет максимизирован или минимизирован. Функционал затрат устойчивой итерационной схемы большинства квадратов может быть записан как (Meju & ​​Hutton 1992) 24, где - множитель Лагранжа, β - коэффициент затухания Марквардта, а L 2 ограничивает энергию модели. шаг коррекции Δ м . В нашем случае пороговое несоответствие выбрано как Q t = Q [ м 0 ] + Δ Q .Экстремум функционала стоимости (уравнение 24) соответствует шагу коррекции модели (Meju & ​​Hutton 1992) 2526, который состоит из коррекции модели методом наименьших квадратов Δ м l s и компенсирующего члена, связанного с к максимальному или минимизированному параметру. Множитель Лагранжа μ определяется квадратичным ограничением: линеаризованное несоответствие Q lin = (Δ d W - J W Δ ) T d W - J W Δ м ) равняется пороговому несоответствию, то есть Q 9047 9047 лин т .Это ограничение дает (Meju & ​​Hutton 1992) 27, где Q lin l s - это линеаризованная функция несоответствия, относящаяся к поправке модели наименьших квадратов Δ m l s .

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

    Коэффициенты демпфирования следует выбирать тщательно, поскольку они влияют на направления шагов коррекции модели. Хотя коэффициент демпфирования можно выбрать так, чтобы доверительная поверхность, соответствующая Q 0 + Δ Q , была достигнута за одну итерацию, нет гарантии, что удельное сопротивление рассматриваемой ячейки будет максимизировано или минимизировано таким образом. Следовательно, необходимо применять итеративную реализацию алгоритма большинства квадратов. Мы ограничиваем диапазон тестируемых коэффициентов демпфирования теми, которые не приводят к несоответствиям, превышающим Q 0 + Δ Q , и выбираем тот, который дает наибольшее отношение изменения удельного сопротивления рассматриваемой ячейки к среднему изменению удельного сопротивления. модели.

    2,4 Подпространство модели

    Использование метода подпространства мотивировано тем фактом, что SVD очень дорогостоящее с точки зрения процессорного времени. Чтобы найти сингулярные значения и собственные векторы модели, время вычислений масштабируется как M 3 (Голуб и ван Лоан 1996), где M - количество параметров модели. Следовательно, мы ограничиваемся рассмотрением подпространства вокруг интересующей модельной ячейки. Протяженность подпространства определяется формой подпространства и радиусами горизонтального и вертикального подпространства s y и s z .На практике наиболее подходящими оказались подпространства полуэллиптической формы с центром на поверхности и радиусы подпространства, выбранные в качестве глубины до рассматриваемой ячейки. Этот выбор оправдан диффузионным поведением электромагнитных полей. С увеличением глубины электромагнитные поля распространяются на больший объем, и разрешение уменьшается. Обзор различных вариантов подпространства приведен в Приложении A.

    2.5 Меры разрешения

    Чтобы получить быстрый обзор свойств разрешения, желательно сжать разрешающее ядро ​​в сокращенный набор мер разрешения.Alumbaugh и Newman (2000) предлагают использовать 50-процентную ширину разброса, которая представляет собой горизонтальную или вертикальную ширину разрешающего ядра между точками, в которых амплитуда составляет 50 процентов от ее максимума. Friedel (2003) сокращает информацию, содержащуюся в матрице разрешения модели, до нескольких репрезентативных показателей, таких как радиус разрешения, полученный из диагональных элементов матрицы разрешения, и флаг искажения, который описывает, является ли положение максимального разрешения для определенной ячейки смещение от самой ячейки.Однако оба подхода обеспечивают ограниченные меры разрешения, поскольку возможное наличие боковых лепестков в разрешающем ядре не допускается. Также не учитывается положение центра разрешения по отношению к рассматриваемой ячейке. Мы следуем определению центров разрешения и длин как интегралов по разрешающему ядру Педерсеном (1977). Этот выбор учитывает форму полного разрешающего ядра на подпространстве. Подробности можно найти в Приложении B.

    2.6 Общая схема

    Общую схему нашего анализа разрешения и дисперсии TSVD можно резюмировать следующим образом.

    • (1)

      Выберите ячейки сетки, для которых должен выполняться анализ разрешения / дисперсии TSVD.

    • (2)

      Выберите пороговое значение дисперсии модели в виде коэффициента f порог для удельных сопротивлений модели (см. Конец раздела 2.1.2).

    • (3)

      Для каждой выбранной ячейки:

    • (i)

      Определите радиусы подпространства вокруг выбранной ячейки, например, на основе глубины ячейки (см. Раздел 2.5 и Приложение А).

    • (ii)

      Для шага анализа = 1,2

    • (1)

      Если шаг анализа = 1, выполнить инверсию MSE в подпространстве, дающую уровень усечения p MSE (см. Конец Раздел 2.2). Если шаг анализа = 2, выполните инверсию усеченного SVD на подпространстве с уровнем усечения p TSVD = p 1 шага анализа 1 (см. Конец раздела 2.2).

    • (2)

      Масштабировать якобиан подпространства последней модели по суммарным ошибкам данных.

    • (3)

      Вычислите SVD якобиана масштабированного подпространства, следуя уравнению. (3).

    • (4)

      Определите нелинейные полуоси в направлении собственных векторов модели с помощью комбинации методов секущей и пополам (раздел 2.2) и замените нелинейные полуоси на линейные полуоси в экв. (23) для дисперсии параметров модели до тех пор, пока не будет достигнут заданный порог дисперсии. Это определяет уровень усечения p 1,2 шага анализа 1 или 2, соответственно, и обычно p 1 > p MSE .

    • (iii)

      Для конечного уровня усечения p 2 этапа анализа 2, вычислите матрицу разрешения подпространства (уравнение 14) и меры разрешения (раздел 2.5).

    • (iv)

      Проверьте отклонения модели с помощью инверсии наибольших квадратов (раздел 2.3).

    3 Синтетический пример

    Синтетические прямые данные простой двухмерной модели с блоком 10 Ом · м в однородном полупространстве 1000 Ом · м (рис.2) были рассчитаны для поперечной магнитной (TM) и поперечной электрической (TE) мод (см. E.грамм. Брюитт-Тейлор и Уивер, 1976; Aprea et al. 1997). Отклики были рассчитаны на 10 приемных участках в течение 12 периодов в диапазоне от 0,0031 до 6,4 с, что дало в общей сложности 480 точек данных. Гауссовский белый шум, соответствующий 2,5% модуля вычисленных импедансов, был добавлен к прямым характеристикам обеих поляризаций.

    Рисунок 2.

    Двухмерная модель с блоком 10 Ом · м, встроенным в полупространство 1000 Ом · м.

    Рис. 2.

    Двухмерная модель с блоком 10 Ом · м, встроенным в полупространство 1000 Ом · м.

    После этого была проведена инверсия синтетического набора данных по инверсной схеме REBOCC. Предполагалось, что минимальная погрешность соответствует 2,5% модуля импедансов, а исходная модель представляла собой однородное полупространство 300 Ом · м. В конце седьмой итерации была получена модель (рис. 3), которая соответствует данным со среднеквадратичным несоответствием 0,88. Дополнительные итерации с REBOCC не привели к дальнейшему снижению среднеквадратичного несоответствия.

    Рисунок 3.

    Результат инверсии 2-D REBOCC синтетических данных из 2-D блочной модели на рис.2. После семи итераций достигается среднеквадратичное несоответствие 0,88. Для ячеек, отмеченных I – III, вычисляются дисперсия модели и свойства разрешения.

    Рисунок 3.

    Результат инверсии 2-D REBOCC синтетических данных из 2-D блочной модели на рис. 2. После семи итераций достигается среднеквадратичное несоответствие 0,88. Для ячеек, отмеченных I – III, вычисляются дисперсия модели и свойства разрешения.

    На основе этой модели REBOCC анализ разрешающей способности и дисперсии TSVD был выполнен для набора из трех ячеек, где первая ячейка расположена прямо над проводящим блоком (отмечена ячейкой I на рис.3), вторая ячейка расположена сбоку от проводника (ячейка II на фиг. 3), а последняя ячейка находится непосредственно под проводящим блоком (ячейка III на фиг. 3). Форма подпространств во всех случаях была выбрана полуэллиптической (см. Уравнение A1). Радиусы горизонтального и вертикального подпространства s y и s z , соответственно, были выбраны как глубина до центра ячейки. Результаты анализа представлены в таблице 1.

    Таблица 1.

    Результаты нелинейного анализа дисперсии и разрешения в примере синтетических данных.

    Таблица 1.

    Результаты нелинейного анализа дисперсии и разрешения примера синтетических данных.

    Сначала обсуждаются свойства дисперсии и разрешения ячейки I поверх проводящего блока. Всего в подпространство входит 168 ячеек. На рис. 4 показано разрешающее ядро ​​для ячейки I. Ячейка I отмечена белым ромбом, а граница подпространства показана пунктирной линией.Значения разрешающего ядра ячеек, относящихся к подпространству, нанесены в соответствии с нормализованной цветовой шкалой справа от графика. Начальная инверсия MSE дает консервативный уровень усечения 71. Для анализа разрешения и дисперсии порог относительной дисперсии модели был установлен на f thresh = 3.0. Как в начальном, так и в конечном анализе дисперсии используется всего 168 собственных векторов и полуосей без превышения заданного порогового значения дисперсии.Линейные и нелинейные полуоси изображены на рис. 5. Линейные полуоси (красная кривая) имеют динамический диапазон от ≈2 × 10 –2 до ≈3 × 10 +3 , тогда как не -линейные полуоси (синие и зеленые кривые) по величине не превышают 1. До числа сингулярного числа 50 линейные и нелинейные полуоси почти равны, что указывает на то, что нелинейность становится важной только для более высоких значений сингулярных чисел. Поскольку нелинейные полуоси увеличиваются очень медленно для сингулярных чисел, превышающих 50, очевидно, что большее количество собственных векторов и полуосей может быть включено во время вычисления дисперсии, если нелинейность разрешена.Следовательно, больше собственных векторов модели входит в оценки нелинейного разрешения, чем в оценки линейного разрешения, а нелинейность имеет свойство фокусировки разрешения. На практике мы обнаружили, что SVD необходимо усекать, когда полуоси имеют порядок 1, если требуются отклонения, соответствующие факторам f thresh = 3. Поскольку нелинейные полуоси меньше 1, этот порог не может быть достигнут, даже если включены все 168 собственных векторов и сингулярных значений, как здесь.Результирующие коэффициенты на удельное сопротивление исследуемой ячейки равны f - = 1,50 в направлении уменьшения и f + = 1,49 в направлении увеличения сопротивления. Поскольку все собственные векторы включены, ячейка разрешается однозначно, и разрешающее ядро ​​отличается от нуля только в самой ячейке I (рис. 4). Этот результат подтверждается инверсией большинства квадратов, которая дает множители f - MSQ = 2,06 и f + MSQ = 1.47 соответственно. Хотя первый коэффициент несколько выше, чем его нелинейный эквивалент f , он явно остается ниже заданного порога. Для сравнения описаны результаты анализа дисперсии линейной модели и разрешения ячейки I. Имея только 83 из 168 собственных векторов модели и линейных полуосей, коэффициент f lin = 2,66 на удельное сопротивление ячейки I приближается к f thresh = 3,0. Линейные полуоси, вычисленные во время линейного анализа, и линейные полуоси, вычисленные во время нелинейного анализа (рис.5) довольно похожи. Длина линейной полуоси для уровня усечения p отсечка = 83 на рис. 5 составляет приблизительно 2, что соответствует приведенному выше эмпирическому правилу. Согласно линейному анализу, элементы разрешающего ядра на рис.6 не равны нулю за пределами ячейки I, а расчетное удельное сопротивление ячейки I преимущественно представляет собой среднее значение удельного сопротивления ячеек, близких к ячейке I, и ячеек, близких к поверхности. . Однако также оценка линейного разрешения показывает, что ячейка I довольно хорошо разрешена.На рис. 6 красными полосами обозначены предполагаемые длины разрешения по вертикали и горизонтали, которые вычисляются в соответствии с уравнением. (БИ 2). Точка пересечения полос - это центр разрешения, рассчитанного по формуле. (B1). Центр разрешения совпадает с ячейкой I, а длина разрешения примерно в 3 раза больше, чем размеры ячейки I. По сравнению с оценками нелинейного разрешения, оценки линейного разрешения относительно плохо разрешенных ячеек имеют тенденцию ухудшаться более резко, чем таковые из достаточно хорошо разрешенных ячеек.Это в основном связано с тем, что элементы в собственных векторах модели довольно хорошо разрешенных ячеек велики в собственных векторах первой модели, которые уже участвуют в линейном анализе, и малы в собственных векторах более поздних моделей (см. Примеры в Jupp & Vozoff 1975 ; Vozoff & Jupp 1975). В отличие от этого, элементы плохо разрешенных ячеек, по-видимому, более сконцентрированы в собственных векторах модели с большим числом сингулярных значений, которые, возможно, не входят в линейный анализ.

    Рисунок 4.

    Ядро разрешения подпространства для ячейки I (отмечено белым ромбом) на верхнем крае проводящего блока. Пунктирная линия показывает границу подпространства, в котором были вычислены свойства дисперсии и разрешения. Значения ядра разрешения вне ячейки I равны нулю, а ячейка I полностью разрешена.

    Рисунок 4.

    Ядро разрешения подпространства для ячейки I (отмечено белым ромбом) на верхнем крае проводящего блока. Пунктирная линия показывает границу подпространства, в котором были вычислены свойства дисперсии и разрешения.Значения ядра разрешения вне ячейки I равны нулю, а ячейка I полностью разрешена.

    Рисунок 5.

    Линейные (красная кривая) и нелинейные полуоси (синяя и зеленая кривые) для ячейки I на верхнем крае блока. До сингулярного числа 50 наблюдается примерно равное увеличение линейной и нелинейной полуосей. После сингулярного числа 50 нелинейные полуоси увеличиваются очень мало, тогда как линейные полуоси продолжают увеличиваться. Следовательно, большее количество полуосей и собственных векторов может быть включено в вычисление дисперсии нелинейной модели до тех пор, пока не будет достигнут заданный порог дисперсии, а нелинейность имеет тенденцию фокусировать разрешение.

    Рисунок 5.

    Линейные (красная кривая) и нелинейные полуоси (синяя и зеленая кривые) для ячейки I на верхнем крае блока. До сингулярного числа 50 наблюдается примерно равное увеличение линейной и нелинейной полуосей. После сингулярного числа 50 нелинейные полуоси увеличиваются очень мало, тогда как линейные полуоси продолжают увеличиваться. Следовательно, большее количество полуосей и собственных векторов может быть включено в вычисление дисперсии нелинейной модели до тех пор, пока не будет достигнут заданный порог дисперсии, а нелинейность имеет тенденцию фокусировать разрешение.

    Рис. 6.

    Ядро разрешения подпространства для ячейки I, полученное из анализа дисперсии и разрешения линейной модели. В отличие от результатов нелинейного анализа, удельное сопротивление ячейки I разрешено не полностью (см. Рис. 4), а среднее удельное сопротивление ячеек вблизи ячейки I и ячеек вблизи поверхности. Значения разрешающего ядра представляют собой веса, с которыми удельное сопротивление каждой ячейки входит в среднее значение.Однако также оценка линейного разрешения показывает, что ячейка I довольно хорошо разрешена. Красные полоски отображают расчетную длину разрешения по вертикали и горизонтали (см. Уравнение B2). Точка пересечения полос определяет центр разрешения согласно уравнению. (B1).

    Рисунок 6.

    Ядро разрешения подпространства для ячейки I, полученное в результате анализа дисперсии и разрешения линейной модели. В отличие от результатов нелинейного анализа, удельное сопротивление ячейки I не решается полностью (см. Рис.4), но среднее удельное сопротивление ячеек вблизи ячейки I и ячеек вблизи поверхности. Значения разрешающего ядра представляют собой веса, с которыми удельное сопротивление каждой ячейки входит в среднее значение. Однако также оценка линейного разрешения показывает, что ячейка I довольно хорошо разрешена. Красные полоски отображают расчетную длину разрешения по вертикали и горизонтали (см. Уравнение B2). Точка пересечения полос определяет центр разрешения согласно уравнению. (B1).

    Затем анализируется ячейка II в резистивном хозяине слева от проводящего блока.Соответствующее подпространство содержит 167 ячеек (рис. 7). Уровень усечения начальной инверсии MSE составляет 49. Следующий анализ начальной дисперсии предсказывает уровень усечения 133 для коэффициентов удельного сопротивления немного меньше 3 (см. Таблицу 1). Этот предполагаемый уровень усечения остается удивительно постоянным после окончательной инверсии усеченного SVD и окончательного дисперсионного анализа. Окончательный дисперсионный анализ дает коэффициенты f ± = 2,99 в направлениях как уменьшения, так и увеличения удельного сопротивления.Длина разрешения как по вертикали, так и по горизонтали примерно в 4 раза больше, чем соответствующие размеры ячейки. Центр разрешения совпадает с рассматриваемой ячейкой. Однако дисперсии занижены, что доказывает инверсия большинства квадратов, которая сообщает изменчивость как f - MSQ = 4,70 и f + MSQ = 18,4. Следовательно, наши оценки разрешения слишком оптимистичны, и разрешающее ядро, соответствующее фактору f ± ≈ 3, меньше сосредоточено на рассматриваемой ячейке, чем показано на рис.7.

    Рисунок 7.

    Ядро разрешения подпространства для ячейки II слева от проводящего блока. Расчетное удельное сопротивление ячейки II представляет собой средневзвешенное значение истинных сопротивлений всех ячеек в подпространстве. Преимущественно резистивность клеток, близких к клетке II, и клеток, близких к поверхности, вносит вклад в оценочное сопротивление клетки II.

    Рисунок 7.

    Ядро разрешения подпространства для ячейки II слева от проводящего блока.Расчетное удельное сопротивление ячейки II представляет собой средневзвешенное значение истинных сопротивлений всех ячеек в подпространстве. Преимущественно резистивность клеток, близких к клетке II, и клеток, близких к поверхности, вносит вклад в оценочное сопротивление клетки II.

    Для процессов электромагнитной индукции обычно считается, что разрешение ниже толстого проводника относительно низкое из-за сильного затухания электромагнитных полей в самом проводе. Мы исследуем разрешающие свойства ячейки III с центром под проводящим блоком.Разрешающее ядро ​​показано на рис. 8. Для максимального периода 6,4 с и удельного сопротивления полупространства 10 Ом · м глубина скин-слоя составляет около 4000 м. Поскольку проводящий блок имеет толщину около 5800 мкм, можно ожидать, что разрешающая способность исследуемой ячейки относительно мала. Всего в подпространство включено 284 ячейки. Расширение подпространства полностью включает в себя проводящий блок (см. Рис. 8). Модель финальной инверсии подпространства TSVD представлена ​​на рис.9. По сравнению с гладкой моделью REBOCC (рис. 3), модель подпространства имеет немного улучшенное среднеквадратичное соответствие данных 0,86. Верхняя и боковая границы блока кажутся одинаково хорошо определенными в обеих моделях. Также в обеих моделях нижний край размыт в большую глубину. Не существует режима или диапазона периодов, в котором одна из моделей показывает систематическое улучшение соответствия данных по сравнению с другой моделью (не показано). Количество включенных собственных векторов одинаково для начального и окончательного дисперсионного анализа (таблица 1), то есть 245.Окончательный дисперсионный анализ дает коэффициенты f - = 2,86 в направлении уменьшения и f + = 2,88 в направлении увеличения удельного сопротивления. В соответствии с нашим ожиданием разрешающее ядро ​​(рис. 8) сильно искажено. Многие клетки, расположенные близко к поверхности, способствуют разрешению ядра, указывая на то, что расчетное удельное сопротивление клетки III является средним по большой площади. Рис. 4 предполагает, что ячейки в верхней части проводящего блока очень хорошо разрешены, в то время как разрешение ухудшается для ячеек глубже в блоке.Поскольку матрица разрешения симметрична, вклад ячеек в верхней части блока в удельное сопротивление ячейки III незначителен. Это приводит к небольшим элементам таких клеток для разрешающего ядра клетки III (см. Рис. 8). Ячейки в более глубокой части проводящего блока менее разрешены, чем ячейки в верхней части, что означает, что их разрешающие ядра имеют ненулевые значения, разбросанные по ячейкам в их окрестностях, включая ячейки под блоком. Следовательно, ячейки оценочной модели под проводящим блоком будут иметь более низкое удельное сопротивление, чем соответствующие ячейки истинной модели, и толстые проводники будут размываться в глубину в оценочной модели (см.рис.8). Из-за больших значений разрешающего ядра ячеек, близких к поверхности, центр разрешения смещен от исследуемой ячейки, а длина разрешения по горизонтали и вертикали превышает размеры ячейки. Инверсия наибольших квадратов дает коэффициенты f - MSQ = 2,73 и f + MSQ = 2,75 (таблица 1), которые подтверждают наши оценки нелинейной дисперсии.

    Рисунок 8.

    Ядро разрешения подпространства для ячейки III, которое расположено под проводящим блоком.Расчетное сопротивление ячейки III сильно зависит от сопротивления ячеек в верхних 600 м модели (то есть ячеек над проводящим блоком), ячеек в нижней части блока и других ячеек, близких к ячейке III. Вклад ячеек в нижней части проводящего блока в удельное сопротивление ячейки III объясняет, почему проводник размывается на большую глубину в оценочной модели (рис. 9).

    Рисунок 8.

    Ядро разрешения подпространства для ячейки III, расположенной под проводящим блоком.Расчетное сопротивление ячейки III сильно зависит от сопротивления ячеек в верхних 600 м модели (то есть ячеек над проводящим блоком), ячеек в нижней части блока и других ячеек, близких к ячейке III. Вклад ячеек в нижней части проводящего блока в удельное сопротивление ячейки III объясняет, почему проводник размывается на большую глубину в оценочной модели (рис. 9).

    Рис. 9.

    Двухмерная модель, полученная в результате окончательной инверсии усеченного SVD в подпространстве, относящемся к ячейке III.Модель соответствует синтетическим данным со среднеквадратичным несоответствием 0,86. В результате уменьшенного несоответствия, полученного с помощью инверсии TSVD, эта модель немного грубее, чем модель REBOCC на рис. 3. Ячейка сверху с черным заполнением имеет удельное сопротивление более 1800 Ом · м, как указано в легенде.

    Рис. 9.

    Двухмерная модель, полученная в результате окончательной инверсии усеченного SVD в подпространстве, относящемся к ячейке III. Модель соответствует синтетическим данным со среднеквадратичным несоответствием 0,86. В результате уменьшенного несоответствия, полученного с помощью инверсии TSVD, эта модель немного грубее, чем модель REBOCC на рис.3. Верхняя ячейка с черной заливкой имеет удельное сопротивление более 1800 Ом · м, как указано в легенде.

    4 Пример полевых данных из района Скидига (Швеция)

    Pedersen et al. Отчет № (2005) по проблеме определения геометрии водоносной системы ледниково-речного горизонта, состоящей из песчано-гравийных отложений в Скедиге (Швеция), с помощью метода радиомагнитотеллурии с контролируемым источником. Нижняя граница системы водоносных горизонтов - кристаллический фундамент, перекрытый системой линз глины.Из-за проводящих глиняных линз сложно определить структуру водоносного горизонта и особенно переход в резистивный фундамент.

    Модель, полученная после четырех обратных шагов с использованием детерминантного режима (Pedersen & Engels 2005) в обратном коде REBOCC (Siripunvaraporn & Egbert 2000), показана на рис. 10. Число точек данных составляет 528, и среднеквадратичное соответствие составляет 1,47. Между профильными измерителями 10 и 100, а также 130 и 220 установлены две токопроводящие глиняные линзы со средним сопротивлением 10 Ом · м (желтого цвета).Pedersen et al. (2005) интерпретирует изолинию 30 Ом · м (т.е. переход между двумя зеленоватыми цветами как log 10 30 1,5) как нижнюю границу глиняных линз. Среднее удельное сопротивление нижележащего водоносного горизонта составляет 100 Ом · м. Примерно на отметке 115 профиля песчано-гравийные образования водоносного горизонта простираются до поверхности. Судя по скважинам в районе профиля, переход от водоносного горизонта к нижележащему кристаллическому фундаменту происходит на глубине около 30 м.В нескольких скважинах глубина до фундамента очень хорошо коррелирует с контурной линией 300 Ом · м на рис. 10 и предполагает удельное сопротивление фундамента более 300 Ом · м.

    Рис. 10.

    2-мерная модель удельного сопротивления района Скедига, построенная как распределение log 10 ρ ( y , z ), где y и z - координаты в направлении профиля и с глубина соответственно. Две токопроводящие глиняные линзы глубиной 3–20 м и удельным сопротивлением до 30 Ом · м (обозначены переходом между двумя зеленоватыми цветами как log 10 30 1.5) видны в левой и правой половинах модели. Примерно на отметке 115 профиля нижележащие отложения песка и гравия простираются до поверхности. По скважинным данным, глубина до фундамента составляет примерно 30 метров. Анализ дисперсии модели и разрешения выполняется для ячеек, помеченных A – D.

    Рис. 10.

    Двухмерная модель удельного сопротивления района Скидиги, построенная как распределение log 10 ρ ( y , z ), где y и z - координаты в направлении профиля, а с глубиной соответственно.Две токопроводящие глиняные линзы глубиной 3–20 м и удельным сопротивлением до 30 Ом м (обозначенные переходом между двумя зеленоватыми цветами как log 10 30 1,5) видны в левой и правой половинах модели. . Примерно на отметке 115 профиля нижележащие отложения песка и гравия простираются до поверхности. По скважинным данным, глубина до фундамента составляет примерно 30 метров. Анализ дисперсии модели и разрешения выполняется для ячеек, помеченных A – D.

    Свойства дисперсии и разрешения четырех ячеек, обозначенных A – D на рис. 10, исследуются в подпространствах, которые имеют форму полуэллипсов (см. Приложение A). Результаты дисперсионного анализа и инверсии большинства квадратов сведены в Таблицу 2. Чтобы четко различать диапазоны удельного сопротивления, разрешенные для глинистых линз и водоносного горизонта, используется пороговое значение коэффициента удельного сопротивления f thresh = 3.0. Для этого выбора диапазон удельного сопротивления глиняных линз составляет {3.3,…, 30 Ом · м} при среднем удельном сопротивлении 10 Ом · м, тогда как диапазон удельного сопротивления водоносного горизонта составляет {33,…, 300 Ом · м} при среднем удельном сопротивлении 100 Ом · м. Следовательно, допустимые диапазоны удельного сопротивления не перекрываются.

    Таблица 2.

    Результаты нелинейного анализа дисперсии и разрешения на примере месторождения Скидига.

    Таблица 2.

    Результаты нелинейного анализа дисперсии и разрешения на примере месторождения Скидига.

    На рис. 11 показаны свойства разрешения ячейки А, расположенной внутри глиняной линзы.Количество ячеек в подпространстве - 184. После включения всех нелинейных полуосей в начальный анализ, идеальное разрешение ячеек и дисперсия, соответствующая коэффициенту f + = 1,75 в направлении увеличения и f - = 1,76 в направлении уменьшения удельного сопротивления ячейки. После окончательного нелинейного анализа коэффициенты удельного сопротивления равны f - = 1,76 и f + = 1,74. Инверсия большинства квадратов подтверждает этот результат, давая коэффициент f - MSQ = 1.42 в направлении уменьшения удельного сопротивления и f + MSQ = 1,46 в направлении увеличения удельного сопротивления.

    Рисунок 11.

    Ядро разрешения подпространства ячейки A (отмечено белым ромбиком), которое находится внутри глиняной линзы, наложенное на log 10 ρ ( y , z ) -изолинии показанной модели на рис. 10. Значения разрешения вне ячейки A равны нулю, и ячейка A имеет идеальное разрешение.

    Рисунок 11.

    Ядро разрешения подпространства ячейки A (отмечено белым ромбом), которое расположено внутри глиняной линзы, наложенное на log 10 ρ ( y , z ) -изолинии модели, показанной на рис. • Значения разрешения за пределами ячейки A равны нулю, и ячейка A полностью разрешена.

    Для ячейки B, которая, как считается, находится в нескольких верхних метрах подвала и под глиняной линзой, разрешающее ядро ​​показано на рис. 12. Первоначальный дисперсионный анализ дает коэффициенты f ± = 2.93 для 214 из 260 собственных векторов. Количество собственных пар остается удивительно стабильным в окончательном дисперсионном анализе, в результате которого было получено 215 включенных собственных векторов и факторов f - = 2,81 и f + = 2,84, соответственно. Результирующее разрешающее ядро ​​(рис. 12) распространяется на пять ячеек как в вертикальном, так и в горизонтальном направлении и центрируется вокруг интересующей ячейки, что предполагает умеренное разрешение. Однако инверсия по наибольшему количеству квадратов дает коэффициент f - MSQ = 24.04 в направлении уменьшения сопротивления и f + MSQ = 46,59 в направлении увеличения сопротивления, что означает, что разрешение сильно завышено.

    Рис. 12.

    Ядро разрешения подпространства ячейки B, расположенной ниже глиняной линзы и в пределах нескольких верхних метров фундамента. Красные полоски отображают расчетную длину разрешения по вертикали и горизонтали (см. Уравнение B2). Точка пересечения полос является центром разрешения согласно уравнению.(B1). Ядро разрешения предполагает, что удельные сопротивления ячеек из водоносного горизонта и фундамента усредняются по расчетному сопротивлению ячейки B.

    Рисунок 12.

    Ядро разрешения подпространства ячейки B, расположенной под глиняной линзой и в пределах нескольких верхних метров фундамента. Красные полоски отображают расчетную длину разрешения по вертикали и горизонтали (см. Уравнение B2). Точка пересечения полос является центром разрешения согласно уравнению. (B1). Ядро резолютона предполагает, что удельное сопротивление клеток из водоносного горизонта и фундамента усредняется по отношению к расчетному сопротивлению клетки B.

    Менее серьезное демпфирование и, следовательно, лучшее разрешение перехода между водоносным горизонтом и основанием можно ожидать в области, где песчано-гравийные образования простираются до поверхности. Исследуем ячейку C в верхней части окна между глиняными линзами (рис. 13). Первоначальный дисперсионный анализ дает уровень усечения 142 из 192 собственных векторов, который увеличивается до 153 с помощью окончательного дисперсионного анализа. Факторы f - = 2,59 и f + = 2.95 получены в направлении уменьшения и увеличения удельного сопротивления соответственно. Инверсия большинства квадратов подтверждает наши нелинейные оценки, давая коэффициент f - MSQ = 3,88 в направлении уменьшения удельного сопротивления и только f + MSQ = 1,53 в направлении увеличения удельного сопротивления. . Следовательно, удельное сопротивление лучше определяется в последнем направлении, которое можно отнести к моде TM, входящей в детерминант.

    Рис. 13.

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

    Рисунок 13.

    Ядро разрешения подпространства ячейки C под окном между линзами глины и внутри песчано-гравийных образований.Значения ядра разрешения, которые значительно отличаются от нуля, обнаруживаются только в шести ячейках, относящихся к водоносному горизонту. Это означает, что глиняные линзы и фундамент имеют незначительное влияние на расчетное удельное сопротивление ячейки C.

    Аналогичное поведение наблюдается для ячейки D наверху фундамента ниже окна между глиняными линзами (рис. 14). Окончательный анализ дисперсии дает уровень усечения 222 с конечными факторами f - = 2,85 и f + = 2.80. Следующая инверсия большинства квадратов показывает, что удельное сопротивление этой ячейки также лучше ограничено в направлении увеличения удельного сопротивления с коэффициентами f - MSQ = 2,03 и f + MSQ = 1,45. Фактически, данный порог дисперсии не достигается, как показывают результаты инверсии наибольших квадратов, и разрешение ячейки D, возможно, занижено. Сравнение с результатами анализа дисперсии и разрешения линейной модели подчеркивает важность учета нелинейности.Во время построения модели подпространства и линейного анализа уровень усечения p отсечки = 95 дает линейный коэффициент f lin = 2,97 на удельное сопротивление ячейки D. Согласно разрешающему ядру ( Рис. 15) из этого линейного анализа, расчетное удельное сопротивление ячейки D больше не является средним для ячеек, которые принадлежат исключительно водоносному горизонту и верхней части фундамента, но также включает ячейки от линз глины и ячейки, расположенные близко к поверхности.Как объяснялось в разделе 3, исключение собственных векторов модели с существенными элементами для ячейки D привело к значительному ухудшению оценки разрешения ячейки D.

    Рис. 14.

    Ядро разрешения подпространства ячейки D под окном между глиняными линзами и в пределах нескольких верхних метров фундамента. Разброс ядра разрешения аналогичен разбросу для ячейки B (рис. 12). В отличие от ячейки B, оценки дисперсии по инверсии большинства квадратов ближе к нелинейным оценкам (см. Таблицу 2), что указывает на то, что ячейка D разрешена лучше, чем ячейка B.

    Рис. 14.

    Ядро разрешения подпространства ячейки D под окном между глиняными линзами и в пределах нескольких верхних метров фундамента. Разброс ядра разрешения аналогичен разбросу для ячейки B (рис. 12). В отличие от ячейки B, оценки дисперсии по инверсии большинства квадратов ближе к нелинейным оценкам (см. Таблицу 2), что указывает на то, что ячейка D разрешается лучше, чем ячейка B.

    Рисунок 15.

    Ядро разрешения подпространства ячейки D, полученное в результате линейного анализа.Ядро разрешения более распространено, чем у его нелинейного аналога (рис. 14), поскольку собственные векторы модели с существенным вкладом в ячейку D не включены.

    Рисунок 15.

    Ядро разрешения подпространства ячейки D, полученное в результате линейного анализа. Ядро разрешения более распространено, чем у его нелинейного аналога (рис. 14), поскольку собственные векторы модели с существенным вкладом в ячейку D не включены.

    5 Обсуждение и выводы

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

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

    В первом случае, когда априорных моделей и априорных модельных ковариаций являются неотъемлемой частью стохастической формулировки обратной задачи, матрицы ковариации и разрешения апостериорных моделей непосредственно задаются хорошо известными формулами ( Менке 1989).

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

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

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

    Надежность нашего дисперсионного анализа можно проверить с помощью инверсии наибольших квадратов, которая обеспечивает независимую и улучшенную оценку изменчивости модели. Согласно результатам анализа наибольших квадратов, порог дисперсии, который используется в нашем нелинейном анализе, может быть скорректирован для получения оценки наибольших квадратов, которая ближе к некоторой дисперсии, которая считается подходящей для задачи.В случае ячейки B модели Skediga, например, порог дисперсии f thresh должен быть уменьшен ниже 3,0, чтобы приблизиться к f ± MSQ ≈ 3,0. Таким образом, уровень усечения будет уменьшен, и будет получена улучшенная оценка разрешения модели.

    Для мелкомасштабной обратной задачи, такой как наш пример поля RMT с примерно 2800 модельными блоками и 500 точками данных, анализ дисперсии модели для подпространства примерно из 200 ячеек занимает около одного часа на современном 32-битном процессоре с частотой ЦП. 2 ГГц.Основные вклады во время вычислений - это вычисление матрицы Якоби и ее подпространства SVD во время инверсии подпространства и большое количество прямых вычислений во время поиска нелинейных полуосей. Для трехмерных или больших двумерных задач время для выполнения прямого вычисления и вычисления матрицы Якоби становится очень большим, а ограничения памяти ограничивают возможность хранения полной матрицы Якоби в памяти. Чтобы сделать наш метод более эффективным, было бы желательно вычислить якобиан только для подпространства.

    Благодарности

    Вирачай Сирипунварапорн из Университета Махидол, Таиланд и Гэри Эгберт из Университета штата Орегон сделали код REBOCC доступным для модификации. Мы также хотели бы поблагодарить трех рецензентов и редактора за их комментарии, которые помогли улучшить рукопись.

    Список литературы

    ,

    2000

    .

    Оценка изображений для двумерной и трехмерной электромагнитной инверсии

    ,

    Geophysics

    ,

    65

    ,

    1455

    -

    1467

    .

    ,

    1997

    .

    Прямая задача электромагнитной индукции: точные конечно-разностные аппроксимации для двумерных дискретных границ с произвольной геометрией

    ,

    Geophys. J. Int.

    ,

    129

    ,

    29

    -

    40

    .

    ,

    1968

    .

    Разрешающая способность валовых данных о Земле

    ,

    Geophys. J. R. astr. Soc.

    ,

    16

    ,

    169

    -

    205

    .

    ,

    1976

    .

    О конечно-разностном решении двумерных индукционных задач

    ,

    Geophys. J. R. astr. Soc.

    ,

    47

    ,

    375

    -

    396

    .

    ,

    1990

    .

    Инверсия Оккама для создания гладких двумерных моделей для магнитотеллурических данных

    ,

    Geophysics

    ,

    55

    ,

    1613

    -

    1624

    .

    ,

    2003

    .

    Разрешающая способность, стабильность и эффективность томографии сопротивления, оцененные с помощью обобщенного обратного подхода

    ,

    Geophys. J. Int.

    ,

    153

    ,

    305

    -

    316

    .

    ,

    1996

    .

    Матричные вычисления, Исследования Джона Хопкинса в области математических наук

    , 3-е изд.,

    Johns Hopkins University Press, Балтимор

    .

    ,

    2002

    .

    Научные вычисления - вводный обзор

    , 2-е изд.,

    МакГроу-Хилл, Нью-Йорк

    .

    ,

    1976

    .

    Most Squares Inversion

    ,

    J. geophys. Res.

    ,

    81

    ,

    1027

    -

    1030

    .

    ,

    1977

    .

    Человеко-компьютерная система интерпретации результатов зондирования удельного сопротивления над горизонтально стратифицированной землей

    ,

    Geophys. Проспект.

    ,

    25

    ,

    667

    -

    691

    .

    ,

    1975

    .

    Простые итерационные методы обращения геофизических данных

    ,

    Geophys.J. R. astr. Soc.

    ,

    42

    ,

    957

    -

    976

    .

    ,

    2006

    .

    Электромагнитные свидетельства существования обода древней лавинной кальдеры на южном склоне горы Мерапи, Индонезия

    ,

    J. Volc. Геотер. Res.

    , в печати.

    ,

    1984

    .

    Учебное пособие: обзор инверсии методом наименьших квадратов и ее применения к геофизическим задачам

    ,

    Geophys.Проспект.

    ,

    32

    ,

    159

    -

    186

    .

    ,

    1993

    .

    Трехмерная магнитотеллурическая инверсия с использованием сопряженных градиентов

    ,

    Geophys. J. Int.

    ,

    115

    ,

    215

    -

    229

    .

    ,

    1992

    .

    Итерационная инверсия большинства квадратов: приложение к магнитотеллурическим данным

    ,

    Geophys. J. Int.

    ,

    108

    ,

    758

    -

    766

    .

    ,

    1989

    .

    Анализ геофизических данных: дискретная обратная теория

    , том

    45

    , Международная серия по геофизике,

    Academic Press

    ,

    London

    .

    ,

    1996

    .

    Расчетно выполнимая матрица приближенного разрешения для обратных сейсмических задач

    ,

    Geophys. J. Int.

    ,

    126

    ,

    345

    -

    359

    .

    ,

    1999

    .

    Оценка глубины исследования в исследованиях удельного сопротивления постоянному току и IP

    ,

    Geophysics

    ,

    64

    ,

    403

    -

    416

    .

    ,

    1977

    .

    Понимание обратной теории

    ,

    Ann. Преподобный "Планета Земля". Sci.

    ,

    5

    ,

    35

    -

    64

    .

    ,

    1977

    .

    Интерпретация данных потенциального поля - обобщенный обратный подход

    ,

    Geophys. Проспект.

    ,

    25

    ,

    199

    -

    230

    .

    ,

    1979

    .

    Ограниченная инверсия данных потенциального поля

    ,

    Geophys. Проспект.

    ,

    27

    ,

    726

    -

    748

    .

    ,

    2004

    .

    Определение уровня регуляризации инверсии усеченного сингулярного разложения: случай одномерной инверсии МТ данных

    ,

    Geophys. Проспект.

    ,

    52

    ,

    261

    -

    270

    .

    ,

    2005

    .

    Обычная двухмерная инверсия магнитотеллурических данных с использованием определителя тензора импеданса

    ,

    Geophysics

    ,

    70

    ,

    G33

    -

    G41

    .

    ,

    1989

    .

    Инверсия магнитотеллурических данных: нелинейный метод наименьших квадратов

    ,

    Geophys. Проспект.

    ,

    37

    ,

    669

    -

    695

    .

    ,

    2005

    .

    Разведка подземных вод с использованием комбинированных методов контролируемого источника и радиомагнитотеллурии

    ,

    Geophysics

    ,

    70

    ,

    G8

    -

    G15

    .

    ,

    2002

    .

    Исследования чувствительности применительно к двумерной модели удельного сопротивления Центральных Анд

    ,

    Geophys. J. Int.

    ,

    150

    ,

    673

    -

    686

    .

    ,

    2000

    .

    Эффективный метод инверсии подпространства данных для двумерных магнитотеллурических данных

    ,

    Geophysics

    ,

    65

    ,

    791

    -

    803

    .

    ,

    1975

    .

    Совместная инверсия геофизических данных

    ,

    Geophys. J. R. astr. Soc.

    ,

    42

    ,

    977

    -

    991

    .

    ,

    1972

    .

    Общая линейная обратная задача: влияние поверхностных волн и свободных колебаний на структуру Земли

    ,

    Rev. Geophys. Space Phys.

    ,

    10

    ,

    251

    -

    285

    .

    ,

    1999

    .

    Расчет разрешающей способности и ковариационных матриц для сейсмической томографии методом LSQR

    ,

    Geophys. J. Int.

    ,

    138

    ,

    886

    -

    894

    .

    Приложения

    Приложение A: Построение подпространства

    На основе радиусов подпространства s y и s z , модельные ячейки со средней точкой ( y , z ), которые являются частью подпространства вокруг интересующей ячейки с средняя точка ( y k , z k ) должны соответствовать одному из следующих критериев эллипса:

    • (i) полуэллипс (A1)
    • (ii) 1-норма (A2)
    • (iii) 2-норма (A3)

    Подпространство в форме полуэллипса (ур.A1) с центром на поверхности ( z = 0) и s y = s z = z k - наш предпочтительный выбор по нескольким причинам. В однородном полупространстве он объясняет, что индукционное пространство сферически увеличивается с периодом (Schwalenberg et al. 2002). Поскольку чувствительность, как правило, наибольшая непосредственно под участками приемника и в проводниках, часто требуется включить относительно большое количество ячеек под соседними участками приемника в подпространстве.Кроме того, полуэллиптическое подпространство дает более полную картину разрешения для моделей с высокопроводящими поверхностными структурами (например, осадочные бассейны в МТ-исследованиях или глиняные линзы в приповерхностных исследованиях), поскольку в этих случаях нижележащие резисторы часто плохо разрешаются и их разрешение сильно искажается проводниками. Остальные варианты формы подпространства в основном используются для уменьшения количества ячеек в подпространстве. Однако радиусы усечения s y и s z должны быть выбраны так, чтобы все ячейки вне подпространства пренебрежимо связывались с номером ячейки k .Существенно, что радиусы подпространства выбраны достаточно большими. Радиусы подпространства, превышающие требуемые задачей, не влияют на наши оценки с опозданием. Их основной эффект - увеличение времени вычислений.

    Приложение B: Меры разрешения

    Журнал Авторы 507 Составление © 2007 РАН

    Разложение по сингулярным значениям

    Анализ факторов разложения по сингулярным значениям

    Факторный анализ - это многомерный статистический метод для извлечения информационного содержания из набора данных.В случае Olis мы имеем большую матрицу, состоящую из оптической плотности и т. Д. Как функции длины волны и времени. Факторный анализ, в случае разложения по сингулярным значениям (SVD) Olis, сокращает это до 3 меньших матриц, которые содержат всю информацию, присутствующую в данных. Кинетическая информация содержится в кинетической матрице собственных векторов. Опыт показал, что для описания данных необходимо лишь небольшое количество собственных векторов. Олис использует 6. Различные процессы скорости хранятся по небольшому количеству собственных векторов, и кинетическая модель может быть адаптирована к ним для извлечения констант скорости.Точно так же спектральная информация хранится в наборе из 6 собственных векторов. Наконец, создается вектор сингулярных значений. Обычно они начинаются с высокого значения и быстро падают до значения, близкого к нулю. Число сингулярных значений выше нуля указывает на количество существенно присутствующих видов и полезно для выбора модели.

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

    SVD - это один из нескольких методов факторного анализа, который, как правило, считается наиболее надежным с математической точки зрения. Модификации Olis к SVD никоим образом не изменяют способ его использования или сделанные на его основании выводы. Единственное отличие - большой прирост скорости.

    Линейная алгебра, лежащая в основе поисковых систем

    Модель векторного пространства обеспечивает основу для большинства алгоритмов поиска информации, используемых сегодня. Однако одной этой самой базовой модели векторного пространства недостаточно.Для ускорения работы базовой модели было изобретено множество модификаций и эвристик, что привело к появлению популярной модели под названием модель Latent Semantic Indexing (LSI) (Berry, Drmac & Jessup, 1999).

    Основным математическим объектом в модели LSI также является матрица. Мы уже обсуждали, как модель векторного пространства, описанную на предыдущих страницах, можно рассматривать как матрицу. Теперь давайте соединим 3 миллиарда векторов документов, связанных с Интернетом (каждый из которых представляет собой вектор размером 300 000 D), друг с другом, чтобы сформировать матрицу размером 300 000 на 3 000 000 000, которая называется поэлементной матрицей .Эта матрица может выглядеть как

    , где ( i, j ) -элемент матрицы содержит 1, если слово для позиции i в документе j является ключевым словом, и 0 в противном случае. Столбец 1 в этой матрице - это вектор документа 1, столбец 2 - вектор документа 2 и так далее.

    Эта огромная поэлементная матрица может содержать избыточную информацию. Берри, Дрмак и Джессап (1999) приводят следующий пример:

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

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

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

    Примечание: Ниже мы кратко представим SVD; Берри, Дрмак и Джессап (1999) и Мейер (2000) содержат подробные объяснения SVD матрицы.

    По сути, усеченный SVD постаментной матрицы сокращает матрицу m на n A до аппроксимации A k , которая хранится компактно с меньшим количеством векторов. Нижний индекс k представляет, какую часть исходной матрицы A мы хотим сохранить. Большое значение k означает, что A k очень близко к A . Однако ценность этого хорошего приближения должна быть сопоставлена ​​с желанием сжатия данных.Фактически, усеченное представление A k требует хранения ровно k скаляров, k m -мерных векторов и k n -мерных векторов. Очевидно, что с увеличением k увеличиваются и требования к хранилищу.

    Давайте сделаем паузу для формального определения SVD с последующим числовым примером. Перед определением SVD мы также рассмотрим некоторые другие термины матричной алгебры. Следующие определения и теоремы взяты из Meyer (2000).

    Определение Ранг матрицы м на n A - это количество линейно независимых строк или столбцов A .

    Определение Ортогональная матрица является реальной матрицей m на n P ранга r , столбцы (или строки) которой составляют ортонормированную основу для R r . То есть набор векторов столбцов (или строк) { p 1 , p 2 ,…, p r }, где r = ранг ( P ), имеет свойство, что p i T p j = 1 для i = j , 0 в противном случае.

    Теорема Для каждых м на n матрица A ранга r , есть ортогональные матрицы

    U м на м и V n на n

    и диагональная матрица

    D r -by- r = диаг (σ 1 , σ 2 ,…, σ r )

    такое, что

    с σ 1 ≥ σ 2 ≥… ≥ σ r > 0.

    Определение Значения σ и в предыдущей теореме называются ненулевыми сингулярными значениями из A . Когда r p = min ( m, n ), говорят, что A имеет p - r дополнительных нулевых сингулярных значений . Факторизация в теореме называется разложением по сингулярным числам (SVD) A , а столбцы U (обозначены u i ) и столбцы V (обозначены v i. ) называются левыми и правыми сингулярными векторами из A , соответственно.В другом варианте SVD выражает A как сумму матриц ранга 1:

    .


    Определение Усеченный SVD - это

    Теорема (Meyer, 2000, стр. 417) Если σ 1 ≥ σ 2 ≥… ≥ σ r ненулевые сингулярные значения A m -by- n , то для каждого k расстояние от A до ближайшего ранга - матрицы k равно


    Таким образом, σ k + 1 является мерой того, насколько хорошо A k приближается к A .Фактически, для A m -by- n ранга r , A k - лучшее приближение матрицы ранга k к A (Eckart & Young, 1936 ). Цель использования A k для приближения A состоит в том, чтобы выбрать k достаточно большим, чтобы дать хорошее приближение к A , но достаточно маленьким, чтобы k r требовало гораздо меньше памяти.

    На рис. 6 в целом показана экономия памяти, полученная с помощью усеченного SVD.

    Рисунок 6: Сжатие данных с использованием усеченного SVD


    Числовой пример может прояснить эти концепции. Предположим, что исходная поэлементная матрица A имеет размер 100 на 300. После усеченного SVD предположим, что A 20 хорошо аппроксимирует A , т. Е. Мы берем k = 20. Чтобы вычислить экономию памяти, мы подсчитываем количество элементов, которые должны быть сохранены в обоих представлениях - полная матрица 100 на 300 A (30 000 элементов) по сравнению с усеченным представлением SVD, которое требует хранения вектора 20 на 1, матрицы 100 на 20 и матрицы 20 на 300 (20 × 1 + 100 × 20 + 20 × 300 = 8,020 элементов).Количество элементов в полной матрице Представление почти в 4 раза больше, чем количество элементов, необходимых для усеченного SVD.

    Большие коллекции данных часто позволяют значительно сэкономить на хранении. В коллекции Medlars используется поэлементная матрица размером 5831 × 1033, которая может быть очень хорошо аппроксимирована усеченным SVD с k = 70, что обеспечивает экономию памяти примерно в 13 раз меньше, чем исходная матрица.

    Мы завершаем этот раздел небольшим примером, взятым из Berry & Browne (1999), показывающим, что усеченный SVD может хорошо аппроксимировать матрицу.Предположим, что A - это матрица 9 × 7

    .

    Значение усечения k = 3 дает усеченное представление SVD.

    Элементы A 3 лишь немного напоминают элементы A - некоторые даже отрицательны - но это не проблема. Нас интересует только способность A 3 моделировать отношения между векторами из A . Предположим, что вектор запроса

    Сравните ранжирование векторов документов с использованием A с ранжированием, полученным с использованием A 3 :
    Таблица 3: Косинусы углов между векторами документа и запросом с использованием A и A3


    По мере увеличения k усеченный SVD лучше отражает суть A , но экономия на хранении уменьшается.Например, это A 6 :

    Обратите внимание, что A 6 гораздо лучше справляется с захватом взаимосвязи между векторами A , чем A 3 . Теперь сравните ранжирование векторов документов с использованием A с ранжированием, полученным с использованием A 6 . Также сравните рейтинг из A 3 с рейтингом из A 6 .

    Таблица 4: Абсолютные косинусы углов между векторами документа

    и запрос с использованием A, A3 и A6

    Обратите внимание, что ранжирование A 6 намного ближе к истинному ранжированию, полученному из полной матрицы A , чем ранжирование A 3 .Но объем памяти, необходимый для A 6 в его более компактной форме ( U 6 , V 6 , σ 6 ), намного превышает объем памяти, необходимый для менее точного A . 3 ( U 3 , V 3 , σ 3 ). Выбор подходящего и разумного значения усечения k - это область активных исследований VSM (Berry, Drmac & Jessup, 1999, Berry & O'Brien, 1998, Letsche & Berry, 1997, Ding, 1999, Zha, Marques & Simon , 1998).

    Продемонстрировав потенциал усеченного SVD для экономии памяти (примечание), давайте теперь рассмотрим, как вычисляются угловые меры с использованием U k , D k , V k . Вспомните формулу для вычисления косинуса угла & theta i между вектором документа d i и вектором запроса q :

    Поскольку d i - это столбец i th для A , давайте заменим d i на A e i , где e i - столбец вектор всех нулей, кроме 1, в позиции i th .Таким образом,

    Используя усеченный SVD, заменим A на его аппроксимацию A k и получим

    Поскольку A k можно разложить на три маленькие матрицы, мы заменяем A k на U k D k V k T :

    Чтобы еще больше сократить объем вычислений во время выполнения, мы положим s i = e i T V k D k .Затем,

    с || U k s i T || 2 = || s i T || 2 для ортогональной матрицы U k , и || s i T || 2 = || s i || 2 .

    Полезно рассмотреть приведенное выше уравнение, подчеркнув преимущества усеченного SVD. Во-первых, U k , s i и || s i || 2 необходимо вычислить только один раз, и его можно повторно использовать для каждого запроса q .Во время работы только U k T q и || q || 2 необходимо вычислить. Векторы в U k и V k вычисляются один раз, и их можно определить без вычисления всего SVD для A , поскольку существует алгоритм для вычисления первых k сингулярных значений и векторов итеративно (Голуб и Ван Лоан, 1996).

    На практике k может быть намного меньше, чем r , и все же A k по-прежнему хорошо приближается к A .Одна из причин этого связана с шумом в A из-за многозначности и синонимичности слов в английском языке. Усечение SVD снижает этот шум. Еще одна причина, по которой k r хорошо работает в приложениях для поиска информации, заключается в том, что числовая точность не имеет большого значения. Например, трех знаков точности в косинусоиде может быть достаточно. То есть косинус 0,231 имеет такое же значение, как и более точный числовой результат 0,23157824. Усеченный SVD дает преимущества при хранении и при этом сохраняет адекватную точность при расчетах угловых измерений.

    Упражнение: поэкспериментируйте с поисковой системой LSI, встроенной в Matlab

    Тест: линейная алгебра

    .
  • Добавить комментарий

    Ваш адрес email не будет опубликован. Обязательные поля помечены *