Для решения трехмерной структурной обратной задачи гравиметрии о нахождении поверхностей раздела в многослойной среде предложены линеаризованные итерационные методы градиентного типа (наискорейшего спуска и минимальной ошибки) с переменными демпфирующими множителями. На основе методов типа Ньютона, Левенберга–Марквардта и линеаризованных методов градиентного типа решения задач гравиметрии для двухслойной и трехслойной среды разработаны эффективные параллельные алгоритмы, численно реализованные на многопроцессорных системах различного типа: многопроцессорном комплексе МВС-ИММ, многоядерном процессоре Intel и графических процессорах NVIDIA, входящими в состав суперкомпьютера «Уран». Проведено исследование эффективности и оптимизация параллельных алгоритмов. Параллельные алгоритмы встроены в разработанную систему удаленных вычислений «Специализированный Веб-портал решения задач на многопроцессорных вычислительных системах». Решена модельная задача для трехслойной среды.