Процедура
MASW_ZV
Построение модели ВЧР в каждой точке ОГТ после
редактирования в пункте меню Processing->MASW
Analysis.
Краткое описание метода MASW и ссылки на публикации смотрите в разделе
Теория и практика
применения метода MASW.
Входные данные определяются автоматически по имени
входного сейсмического файла и кода сортировки сейсмограмм:
1. При сортировке по ОПВ или ОМГ
inFileName_crs_S_pk_New.txt.
2. При сортировке по ОГТ
inFileName_crs_pk_New.txt.
Результаты выводятся в текстовые файлы:
inFileName_vels_calc.dat - модель
ВЧР в каждой ОГТ.
inFileName_disp_calc.dat -
начальные и рассчитанные дисперсионные кривые в каждой точке ОГТ.
inFileName_masw_out.txt -
протокол прохождения задания.
Полученные результаты можно визуализировать утилитой MASW Plot.
Параметры:
Random seed - Случайное начальное число
(произвольное положительное нечетное целое число,
состоящее из 9 цифр.
Например 999999999 или 777777711).
Применяется при отсутствии начальных значений
мощностей H и скоростей Vs
в файле Init model Flile
(Параметр N_initial Flag=0).
Начальные значения мощностей H и скоростей Vs вычисляются как случайные числа в заданных
интервалах. Поскольку результат зависит от исходных скоростей, заданных
случайными числами, настоятельно рекомендуется применять эту программу
несколько раз с различными значениями этого параметра.
Следующие три параметра описывают сужение коридора
поиска значений H и Vs на
каждой "глобальной итерации".
Т = Т0 * ехр (-(C * K ^
A)), где K - номер итерации в пределах одной "глобальной
итерации".
Init Reduction - Начальный коэффициент сужения на каждой
"глобальной итерации".
Рекомендуется 1.0
Coefficient A - Коэффициент сужения A.
Рекомендуется 0.6
Coefficient C - Коэффициент сужения C. Рекомендуется
1.3
Reduction Nmax - Максимально допустимое значение количества
"глобальных итераций".
Параметр должен быть в пределах 5-10.
Iteration Number
- Максимально допустимое значение количества итераций на каждой "глобальной
итерации".
Параметр должен быть в пределах 10-32000.
Общее максимально допустимое количество итераций
определяется как:
MaxIter=Iteration
Number * Reduction Nmax * Nlayers,
где Nlayers -
количество слов в модели ВЧР.
Averaged Dev
- Пороговое значение стандартного отклонения решения
на К-ой итерации от решения на К-1-ой итерации.
При достижении этого значения итерации прекращаются.
N_vp Flag
- Флаг применения эмпирического соотношения скоростей Vp и Vs по Ludwig и
др (1970).
2 - применяется формула Vp = Vs
1 - применяется формула Vp = 1.11*Vs+1.29
0 = Vp прикрепляется
к начальным значениям в файле "Init Model Flile".
N_roh Flag
- Флаг применения эмпирического соотношения плотностей и скоростей Vp по Kitzunezaki и др (1990),
1 - применяется формула roh=1.2475+0.399*Vp‐0.026*Vp*Vp
0 = Vp прикрепляется
к начальным значениям в файле "Init Model Flile".
N_initial Flag
-
0 = начальная модель устанавливается с помощью
случайных чисел.
1 = начальная модель берется из файла Init Model Flile.
2 = начальная модель берется из файла Init Model Flile. и на каждой следующей точке в качестве начальной
модели используется рассчитанная на предыдущей точке.
Режим N_initial Flag=2 пока не работает.
Init Model Flile - Текстовый файл с начальной моделью ВЧР.
Формат файла определяется параметром N_initial Flag.
При N_initial Flag=1,2
Iwaki City, SS-1,22/12/2012 :Комментарий
5 :количество слоев
1.9 1.5 0.002 0.010 0.06
0.18 0.0022 0.069 :den,Vp,hmin,hmax,vmin,vmax,hini,vini
1.9 1.5 0.002
0.010 0.08 0.15 0.0036 0.092 :den,Vp,hmin,hmax,vmin,vmax,hini,vini
1.9 1.5 0.002
0.010 0.06 0.18 0.0049 0.160 :den,Vp,hmin,hmax,vmin,vmax,hini,vini
1.9 1.5 0.002
0.010 0.06 0.20 0.0035 0.175 :den,Vp,hmin,hmax,vmin,vmax,hini,vini
2.0 1.70 980.0
999.0 0.16 0.30 990.0 0.280
:den,Vp,hmin,hmax,vmin,vmax,hini,vini
den - плотность каждого слоя.
Vp - скорость продольных волн.
(hmin, hmax): Минимум и максимум интервала поиска мощностей H каждого слоя.
Для самого глубокого слоя должны быть (980.0, 999.0).
(vmin, vmax): Минимум и максимум интервала поиска скоростей Vs каждого слоя.
(hini,vini): Начальные значения мощностей H и скоростей Vs для каждого
слоя.
При N_initial Flag=0
Iwaki City, SS-1,22/12/2012 :Комментарий
5 :количество слоев
1.9 1.5 0.002 0.010 0.06 0.18 :den,Vp,hmin,hmax,vmin,vmax
1.9 1.5 0.002
0.010 0.08
0.15
:den,Vp,hmin,hmax,vmin,vmaxi
1.9 1.5 0.002
0.010 0.06
0.18
:den,Vp,hmin,hmax,vmin,vmax
1.9 1.5 0.002
0.010 0.06 0.20
:den,Vp,hmin,hmax,vmin,vmax
2.0 1.70 980.0
999.0 0.16
0.30
:den,Vp,hmin,hmax,vmin,vmax
den - плотность каждого слоя.
Vp - скорость продольных волн.
(hmin, hmax): Минимум и максимум интервала поиска мощностей H каждого слоя.
Для самого глубокого слоя должны быть (980.0, 999.0).
(vmin, vmax): Минимум и максимум интервала поиска скоростей Vs каждого слоя.
Примечание
При возникновении ошибок в процессе решения обратной
задачи процедура выводит в файл inFileName_masw_out.txt коды
этих ошибок:
Ier=1 - Медленная сходимость решения.
Ier=2 - Корень не найден.
Ier=5 - Количество итераций превысило заданный лимит.
Ier=-100*N - Ошибка во
входных данных для слоя с номером N.
Ier=-1 - Ошибка во входных данных обнаружена в функции
RAYDSP.
Ier=-2 - Ошибка во входных данных обнаружена в функции
RAYDSP.
При отрицательном значении ошибки процесс на данной
точке останавливается.