Building/assets/code

Расчётные программы

Программы написаны на Python 3, команды запуска выполняются из корня репозитория. Полные исходные тексты вместе с результатами запуска приведены в отчёте в приложениях «Расчёт коэффициента формы R» и «Расчёт напряжений в осях разрезов».

stress_inversion.py — обратный расчёт девиатора напряжений

Линейная схема Майкла (Michael, 1984) по механизмам очагов землетрясений. Каталог событий и контрольные значения заданы в тексте программы, входные файлы не используются.

• запуск: python assets/code/stress_inversion.py
• зависимость: NumPy
• результат: печать в консоль обратных расчётов и их сверка с контрольными значениями Barth & Wenzel (2010)

profile_stress_ratio.py — напряжения в осях разрезов

Пересчитывает природное поле напряжений в оси вертикальных разрезов и вычисляет отношения sigma_x/sigma_y для окна Initial field stress ПК Prorock.

• запуск: python assets/code/profile_stress_ratio.py
• зависимость: SymPy
• вход: assets/tables/profile_azimuths.csv
• выход: assets/tables/profile_stresses.csv и assets/tables/section_stresses.csv (те же значения печатаются в консоль)

hu_from_slide.py — средний коэффициент Hu из моделей Slide2

Извлекает из архивов моделей Slide2 (.slim) поровое давление, посчитанное фильтрационным анализом (узловые пьезометрические напоры из файла результатов .w01, координаты узлов из .slw, рельеф и удельный вес воды из .sli), и подбирает единый коэффициент Hu формулы Prorock P = gamma_w h H_u: основной оценкой служит МНК через начало координат, справочно печатаются поточечные среднее и медиана, а также значение у изобары 1 МПа (прежний ручной способ).

• запуск: python assets/code/hu_from_slide.py "путь к .slim или каталогу"
• зависимости: нет (стандартная библиотека)
• вход: файлы .slim с посчитанными результатами (Compute)
• выход: assets/tables/hu_from_slide.csv (те же значения печатаются в консоль)