Электроразведка методом сопротивлений - Шевнин В.А.
ISBN 5-211-03303-5
Скачать (прямая ссылка):
u В результате был получен алгоритм, который объединил S себе достоинства обоих методов решения обратной задачи: 1) широкая область сходимости (в качестве начального приближения достаточно указать число слоев и примерные мощности);
1«.-. 2) заметно более высокая скорость сходимости, чем у метода наискорейшего спуска.
На рис.2.4.5. сопоставляется скорость сходимости метода йьютона и наискорейшего спуска для трехслойной модели. Такое соотношение между двумя методами является типичным.
На рис. 2.4.6 показано достаточно плохое начальное при-
Начальное приближение
По методу гкввдообращения
I
MD
10
По методу Зейделя
г*
3L-
t
¦ ¦¦ ,
s>
г
—1—1-
I I 11
>
I
о
Скорость сходимости
мТтоды.
-^^Доммоинирсмм
Рис.2.4.6 Пример подбора при плохом начальном приближении
ближение - все мощности сильно занижены. Но использование метода Зейделя позволяет подобрать модель даже из такого плохого начального приближения. Необходимо обратить внимание на то, что на четвертой итерации невязка резко возросла. Большинство авторов различных алгоритмов минимизации указывают, что это является признаком расхождения алгоритма. Однако из графика сходимости предложенного метода видно, что на самом деле он все-таки сходится. Метод Зейделя используется для автоматической минимизации невязки в алгоритме подбора программы IPI (с 5-й версии), входящей в состав пакета одномерной обработки и интерпретации данных ВЭЗ IPMD [40,84].
Пакет IPI с 1990 г. прошел большую проверку на практических материалах ВЭЗ по многим регионам. Это позволило значительно усовершенствовать процесс 1D интерпретации. Основная концепция пакета IPI-1D - профильная (групповая) интерпретация ВЭЗ. Вплоть до 6 версии программа IPI работает с двумя графическими экранами: на одном изображаются разрез кажущихся сопротивлений и геоэлектрический (по результатам интерпретации), а на другом -интерпретируемая кривая ВЭЗ. Быстрый переход от кривой. ВЭЗ к разрезам по профилю (до 100 кривых ВЭЗ на одном профиле) и обратно, позволяет увидеть, как результат интер-
претации конкретной кривой вписывается в общую модель разреза по профилю, и скорректировать модель для этой кривой ВЭЗ. Профильная интерпретация позволяет получать геоэлектрические разрезы наиболее полно отвечающие геологическим требованиям и априорным данным (по критериям гладкости границ, устойчивости электрических свойств, мощностей отдельных слоев, данным буровых скважин и т.д.); По мере развития пакета программ IPMD1 продолжал совершенствоваться аппарат для автоматической интерпретации кривых ВЭЗ, один из этапов такого развития приведен в настоящем разделе. Но чем больше и точнее объем геологической априорной геологической информации о разрезе, тем меньшую помощь мог получить интерпретатор от алгоритма автоматической минимизации, тем в большей степени интерпретация становилась ручной (неавтоматической) и зависела от геологической квалификации интерпретатора. Поэтому в новой версии программы IPI основное внимание обращено на средства работы с программой, позволяющие облегчить процесс именно геологической интерпретации. Для этого интерпретируемая кривая ВЭЗ, модель среды для этой кривой и весь разрез или его часть помещаются на одном экране. Интерпретатор может менять рамки окна просмотра разреза, помещая в него часть профиля, непосредственно примыкающую к интерпретируемой кривой, или весь разрез, видеть его с учетом рельефа дневной поверхности, менять горизонтальный и вертикальный масштабы.
Собственно автоматическая интерпретация как минимизация невязки теоретической и экспериментальной кривой остается составной частью пакета интерпретации, и на отдельных этапах работы может использоваться для некоторой оптимизации решения, но не является главным средством полечения решения. Минимум невязки лишь один из критериев решения, но не единственный. Аспирантом кафедры геофизики Н.В.Мыциком в 70-х годах для ЭВМ БЭСМ-6 была разработана программа MIDI (максимального использования дополнительной информации). В алгоритм автоматической минимизации вносилось множество условий: на пределы сопротивлений и мощностей, допустимую изменчивость параметров по профилю, в том числе не только р и п, но и SmT слоев и их групп. Опыт работы с программой показал очень высокую трудоемкость задания до начала интерпретации
всех этих условий и ограничений. Многие условия формируются собственно в процессе интерпретации. Современный персональный компьютер, обладающий высокой скоростью счета и высоким качеством графического отображения данных позволяет человеку активно проводить интерпретацию с использованием большого числа различных критериев, часть которых формируется в процессе интерпретации, а в качестве исходных данных используются лишь результаты измерений (кривые ВЭЗ). Попытка передоверить компьютеру весь процесс интерпретации, т.е. опора лишь на автоматическую минимизацию невязки крайне редко может дать хотя бы приемлемый результат. Так как даже в самых спокойных условиях на данные ВЭЗ действует множество искажающих факторов и более 70% кривых ВЭЗ искажены (см.главу 3), самое точное совмещение с искаженной экспериментальной кривой не может дать точного решения.
