Бесплатный автореферат и диссертация по геологии на тему

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

ВАК РФ 04.00.12, Геофизические методы поисков и разведки месторождений полезных ископаемых

Автореферат диссертации по теме "Интерпретация электрических зондирований квазигоризонтально-слоистых сред"

. ПЕРМСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

¿9-

§

^ На правах рукописи

Ч

КОЛЕСНИКОВ Владимир Петрович

ИНТЕРПРЕТАЦИЯ ЭЛЕКТРИЧЕСКИХ ЗОНДИРОВАНИИ КВАЗИГОРИЗОНТАЛЬНО-СЛОИСТЫХ СРЕД

Специальность 04.00.12—Геофизические методы поисков и разведки месторождений полезных ископаемых

Автореферат диссертации на соискание ученой степени доктора технических наук

Пермь — 1995

Работа вьгаолнеяена в Горном институте Уральского отделения Российской академии наук и Пермской сельскохозяйственном институте имени академика Д. Н. Прянишникова.

Официальные оппоненты:

доктор физико-математических иаук, профессор М. Н. Юдин; доктор физико-математических наук П. С. Мартышко; доктор технических наук, профессор Б. И. Тульбович.

Ведущая организация: "Уралгеолком"* Комитета Российской Федерации по теологии и использованию недр (г. Екатеринбург).

Защита состоится 1В мая 1095 г. в 14-00 часов иа заседании диссертационного совета Д 063.59.02 при Пермском государственном университете по адресу: 614600, г. Пермь, ГСП, ул. Бук^-.ииа, 10.

С диссертацией можно ознакомиться ь библиотеке Пермского государственного университета.

Автореферат разослан "_/" ^^ 1995 г.

Ученый секретарь диссертационного совет' канд. геол.-мив. наук, дм

В. А. Гершаиок

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность пробхеш. Метод вертикального влектр5!ческого soh-дироваикя является одним иэ базовых кетодоз элсотрорагведга. С 30-х годов нпясго столетня он пряж'.яется при геологическом картировании, поисках ч разведке нефтяпда и рудных месторождений, поиске строительных материалов, инженерно-геологических, ггвдрсггологл-ческих, мелиоративных изысканиях, реаеяии г1ях5огических задач. Благодаря яеслогиостя приемов наблядеиия, ниэкой себестоимости и экологячности полевых работ, постоянно^ ссвериенствовада» приемов и методов интерпретации, сфера приглжмости его непрерывно ркси-ряетсл.

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

Несмотря на глубокие теоретические разработки з области репе-яия прям« вадвч (Азыган, Даев,Карш:ек"$,1985; Дми'фТ'.ев.Юбб; Hpti-ев,19бб; Кусков,1979,1935; 5Екабарня,1685; Глот,1971; Юдин, »32, 1994; Ядэноч,199£>; Таборовский, 1939; Хомэя,1981 и Др. ) и создаете фундаментальных основ решения некорректных' дадач (Тихонов, 1963; 1986; Лаврентьев, 1963; Иванов,1SS4; Вейсус, Гильберт, 1 £6'/., 1 £59), спецяфюга метода не позволяет формально испояьвовата их для истолкования электрических гоидировзлк?' и требует решения целого ряда сложных вопросов: исследования информационных возможностей метода и оценки достоверности реь-еняй, формализации объекта ксследовезкя адекватно информационным возможностям кетодк, разработки гйюхтйт-ннх приемов доопределения вадачи в условиях ограниченной епркориой кнфоригации, ревеяия методологических вопресоз, обеспечивагжщх создание высокоавтоэдагиенровашак кшпгютержм технологий.

Вместе с тем концепции количественно;! интерпретации одиночках

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

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

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

Основные задачи исследований.

1. Разработка и обоснование методологических принципов Построения автоматизированных систем интерпретации площадных съемок электрического зондирования.

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

3. Разработка и совершенствование принципов автоматизации геологического истолкования-результатов зондирования.

4. Оценка искажений полевого материала и разработка эффективных способов их минимизации.

5. Исследование и равработка способов оценки неодновначности решения.

6. Разработка методики комплексной интерпретации совокупности

зондирований в' рамках единого.информационного пространства "модель-поле'1.

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

Основные научные пололенкя, выносимые на защиту:

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

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

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

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

Научная новизна.

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

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

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

3. Научно обоснован и разработан интегральный принцип построения физика-геологической модели среды, заключающийся в согласованной интерпретации совокупности электрических зондирований при разделении процесса на два отличающихся по своей природе и особенностям проявления некорректности этапа: физического и геологического, вовлекаемых в итерационный процесс на интегральном уровне "модель- поле".

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

б. Исследовано влияние поверхностных неоднородиостей (карстовых воронок, выемок и др.), отдельного вида техногенных помех (трубопроводы, кабели, линейно вытянутые объекты относительно малого сечения) на результаты электрического зондирования.

6. Ра&работакы способы минимизации уровня погрешностей наблюдений, основанные на применении вариационного принципа метода регуляризации.

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

8. Разработана автоматиз;фовадаая система и компьютерная технология обработки и интерпретации данных электрического зондирования 'и отдельных модифккиций метода сопротивлений (электрического профилирования, срединного градиента и т.п.), обеспеченная графическим сервисом для виауалиаацда объемного отображения и сканирования поля, графической редакции и получения жестких копий с выдачей их на плоттер, либо стандартный принтер в веданном масштабе.

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

Практическая ценность и реалнзатая работы.

Ивложенная в работе концепция, способы, алгоритма, компьютерная технология, ' реа^иэоганаые в виде автоматнэ j:роваиной системы, предоставляют пользователе шнрокие вовькшюсти для коделирования и ревения м=:-окого круга эадач обработки, интерпретации, визуализации и графического представления результатов. Испольвозсгше chctî-мн повволяет повысить точность, надежность и полноту ивсаечгнкя информации, обеспечить оперативность и спивать ватратн на обработку полевых материалов, система внедрена я испогьвуотся пря ревс'шя геологоравведочных, поисковых инженерно- геологических и экологических вздач во многих производственных организациях н НХН (ТГО "Уралгеология", Челябинским геолмом, Институт Геофиеига УрО РАК я др.). Об™кй этономшесжя! эффект,, подтвержденный аитемн инедрекуй, составляет около^ЮО тыс. руб. (в ценах 1990 г.). Отделы«» разработки передали для проведения ifcc.1едсз.чтельсгеих работ и обучения студентов в Псрмский а Московский государственник университеты; Уральскую государственную гсрно-геоютическу» екадею®.

Апробация работы.

Основные результаты работы излоиекн в монография, спреЕочк;т;;ЕХ по электроразведке (1S78, 1989) и дояскенн на 5-й,6-5% и 11-й Всесоюзных пколах-семинарах (г.Кизв,1е78; г.Мосюа, МГУ,1981; г.Кша, 1987), на восьмом научно-техническом сгиккаре-совеягрязш (r.Epesaa, 1985), научно-техннгаеския конференциях я региоиальякх ссведв&гях (г.Пермь,1980, 1963,1984,1905,1986,1887,1988,1S&9), m 6-Й И б-â Уральской научной конференции (г.Свердловск, 1655; г.Челябинск, 1989), меяяуяародоом семинаре по горной гесфгоике (г<Пер8Л>,1930), всероссийском научно-техническом совещании по применен»» гоо^шн-ческих методов при решении геологических, штаяерво-геологических н экологических задач (г.Пермь,1994).

Публикации.

Пй теме диссертации опубликовало 62 работа, з том числе одна монография, отдельные параграфы в двух справочниках по етФктрорав-ведке и,одной монографии, статьи и теенсы докладов.

Структура и обьем работы.

Диссертационная работа состоит из введения, аести глав к заключения. Обьем работы 293 стр. машинописного текста, включая 8 таблиц и. 61 рисунок. Библиографический список содержит 264 наименования.

Исходные материалы и личный вклад агл'ора.

Диссертация отображает ревудьтаты исследований, проводивиихся автором в период 1Й72-1995 гг. по госбюджетной и хоздоговорной тематике, саачааа - на кафедре геофизики Пермского университета, а с 1988 года в Горном институте УрО РАН.

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

На разных этапах работы автор использовал помощь и советы коллег по работе д.т.н. М.С.Чадаева, профессора Б.А.Спасского, доцента й.М.Скумбина, поддеряку и участие сотрудников Iirp "Ураагеоло-, гия" В.Ф.Сомоаа, В.М.Соловьева и Г.Л. Сомовой, . Л.М.Хайкельссна, конструктивные замечания и доброжелательную критику академика АЕН, профессора М.Н.Вердичевского и академика АЕН, профессора В.И.Дмитриева. За что автор выражает км искреннюю привнатеденость.

Автор особо благодарен профессору Б.К.Матвееву и чл.-корр. АЕН, профессору В.М.Новоселицкому за плодотворное творческое общение и научные дискуссии, которые во многом стимулировали проведение исследований и написание данной работы.

СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе, которая называется ''Теоретические основы методов интерпретации электрических зондирований", дан критический обзор существующих способов, алгоритмов и программ интерпретации

ВЗЗ, рассмотрены особенности метода электричесшго аондиров&чия с позиции теории реиенкя некорректных задач, дано обоснование методологического подхода к построению автоматизированных сисгег!.

История формирования и развития методов автоматизированной интерпретации данных электрического зондирования бэрет сгое начало с момента появления первых электронных вычислительных мазия и охватывает более чем тридцатилетний период. За это время по данной теме опублмювакы сотни научных работ, предложены и опребэваны десятки алгоритмов и программ, реализующие различные способы решения данной задачи (Иоотова.Хорев, 196S; Шлабарня, 1955,1371,1936; Колесников,1976,1991,1994; ДиансЕа,1984; Мнцик,1682,1994; Дорохова,Шев-нин,Бахирсз,1937; Рыков,Каринская,1981; Петрухин, 1986,1993; Уо-zoff,1953; Kunez,Rocroi,1970; Johansen,1975,1977; Zohdi,1974,1989; Onodera, Fukuda,Sasaki,1980; Lee,1931; Szaraniec, 1979,1982; Kece-' Ii,1933; Basokur,1984; Koefoed,1068,1970,1976,1984; И др.) , созданы автоматизированные системы, испольг.уемые в практике электро- разведочных работ (Зинеиберг, Кирзин, 1974; Колесников, 1931, 1989,1994; Шкабарня,1986, 1989 и др.). Процесс совераекствовакия алгоритмов осуществляется в соответствии с непрерывным развитием теории я стремительным нараазгаанием вычислительных и диалоговых вотжностей SBM, отобраяая общую тенденции: от автоматизации отдельных вычислительних процедур к автоматизированным диалоговым системам и интеллектуальным терминалам.

Большой теоретическим вклад в развитие машинных методоз ин-терпретвпции внесли'отечественное ученые: Л.Л.Ваньян, В.Н.Страхов, А.Н.Тихонов, В.И.Дмитриев, М.Н.Бердичевсши, В.К.Матвеев, Н.Г.Шкабарня и зарубежные исследователи: К.\'02cff, D.Ghosh, Gi'.-'unotz, О.Koefoed и др.

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

Трудности решит о «к проблем определяются как некорректность» обратных задач, таг» к особенностями геофизических исследований, овязезЕЫх о изучением физического поля, допускающего строгие гкшгоествсишз оценки решения и геологического объекта, кото-рыа трудно представить априори в формализованном веде.

Вагьш ксиэвтои г/ ревзюш некорректных, практически недоопре-деленш;; задач является вопрос о выЗоре и использовании дополнительной юфориацди. К настоящему времени сформировалось два основами напраагаша, определяют методологические подходы к исполъзс-вдош дополяктед&шй информации. Первое из них связано с рая работ-гай равличных Ег.р;заптоз метода подбора, второе - направлено па совда;ше к развитее пряшх численных способов..

В способам гштерлротацта, ссновашых на подборе, формирование модели начального прибдкт.ешя • и дополнительной информации о ней кйккянюя на предьмтерпротацяакный отап. Это облегчает формализацией жтерпротгционного процесса, который сводится к вадачэ мкнима-ездиз ссэтвотствуыдо; функционалов, достаточно полно разработанной 2 категаяг-^. Олнзко гаг-сй подход, удобный в математическом отно-В8ЮП, часто осваивается малоэффективным в условиях прьеткческого . производства (ссоСеию при дефиците априорной информации), поскольку значительная часть работы интерпретатора по существу пере-гпсигся из стадии иотоеюзепкя в стадии подготовки исходных данных, к, что главное, не обеслечгшает достаточного физического обосновался интерпретационных моделей.

Эта трудности значительно воэрмтаст при переходе к изучений сгихагопоотрозкшн сред (!Екабарпз,193б), где конструирование начального приближения модели, . су^-зстеевно йшгаизее на конечный рэ-зультет, представляет елейную задачу, свсэеннуп с определением набора модельных элементов и соотгетствувдих га блсдов. на плездэдк кссдэдоаа»», кнааззокм изменения пгра:етров локальных объектов и слоистой среды, углзанием расположения профилей относительно яеод-кородвостей и др. Нспольэуеные для этого пркеки предварительного снаяига каблядеипых п&ргг<етров далеко не всегда позволит с уверенность» подучить подобную информацию.

Ь'айолъгше успехи в развитии метода подбора для горизонтально- едоке моделей сред получеки благодаря исследованиям Е.Б.Иво-тэаой, О.А.Хор&ва (1971); Х.Йохансена (1974); А.А.Рыжова. (1981); Н.В.кЬщика (1984); Л.Н.Пороховой (1971,1987) и др., для слагаю-

построенных сред - И.Г.Шкабарни (1S86); В.В.Кускова, А.Г.Яковлева, И.Н.Модина, 8.А.Шэвнина (1089,1992) и др.

Второе направление свяеано с разработки и развитием прямых численных способов анализа, предполагахвдх двухэтапньш процесс интерпретации [1,4,31,62]: а) определение физических параметров среды непосредственно по результатам полеия наблюдений без априорного конструирования модели среды и б) построение геоэлектрэтескик моделей в рамках физически обосновании результатов с учетом имеющейся априорной информации. Созданием алгоритмов численной интерпретации занимались Н.Г.Шкабария, Т.К.Куничкина (1965,1971); В.Н.Страхов (1966,1976); G.Kunetz, J.P.Rocroi (1970); O.Koefoed (1965,1968,1984); В.П.Романов (1969); В.И.Гудоь (1972,1976); Б.К. Матвеев (1974); В.П.Колесников (1975,1976,1989); О.А.Хачай (1978); 0. Hohman, A.Raich (1987); П.С.Мартышко (1994); T.l.ee (1981); A.Kesli (1933); Лам Куанг Тхиеп, Ле Вуег Зы Хыонг (1984) и др.

Отмеченные направления, в реализации обратной задачи электрического зондирования не исключают, а, скорее, дополняют друг друга, "высвечивая" проблему с равных точек зрения.

К особенностям метода; которые необкод'лко учитывать при истолковании наблвденкш результатов, можно отвести следукщяе:

1) 'вследствие практической недоопределенности физической моделью, отвечающей информационным возможностям метода, является обобщенная модель, объединяющая некоторую совокупность геологических образований и, в силу своей природу, облвдащая свойствам! анизотропии (благодаря • чему физические модели по геометрическим параметрам могут иногда в несколько раз отличаться от их геологических, аналогов); в случае слоистых сред резением эадачи на фдая-ческом уровне строгости является совокупность обобщенных слоев с параметрами hi-Xihi, рц- pi/Xi„ pti- >iPi, где h-- кстшшая мощность обоб!ценкого слоя; Pi, pt - соответственно его среднее продольное и поперечное сопротивления, - коэффициент внквотропни, величина которого для осадочных образований достигает нес кальки», единиц' и не; определяется непосредственно по результатам наземнкх наблюдений методом ВЭЭ (Каленов,1957; Пылаев,1968; Редсаубов, 1973; Колесников 1975 и др.);

£) формирование обобщенных слоев происходит не произвольно, а йо сложным законам, определяемым "вкладом" в форкироваяне электрического поля каждого ив слоев в отдельности ti,2,9], причем часто

данному распределении электрических свойств с глубиной может соответствовать несколько вариантов расчленения разреза на обобщенные слои; это наглядно проявляется при моделировании обратной задачи ВЗЭ по данным электрического каротада, при использовании так называемых кривых Дар-Зарука, в процессе численной интерпретации, анализе параметрических зондирований 11,2,4,34};

3) определение оптимального числа обобщенных слоев, обеспечивающих максимально возможную детальность расчленения разреза при практически приемлемых пределах неоднозначности решения, представляет славшую задачу, т.к. информативность метода снижается как при попытках чревмерной детализации, так и при чрезмерном еагрублении. модели [Колесников,1981; Бердичевский, Дмитриев,1992 и др.].

Таким образом, в общем случае обратная вадача недоопределена и неустойчива как относительно числа обобщенных слоев (моделей) и их параметров, включая анизотропию, так и возможных вариантов формирования геоэлектрических горизонтов в реальных условиях. То есть наряду с математическими (некорректностью обратной' задачи), существуют и физические причины неоднозначности решения. Решением обратной задачи является не математическая, строго детерминированная, а физико-математическая модель, определяемая физическими и методическими условиями ее формирования в электрических (влектро-магнитных) полях в виде обобщенных комплексов, имеющих аномальное проявление в наблюденных полях. Эта модель зависит от методики и электрических свойств слагающих разрез пород, и потому не может иметь строгого априорного представления.

Многсфактсрщй характер недоопределенностй обратной задачи является причиной ситуаций, когда для соседних зондирований исследуемого участка могут быть получены равные варианты расчленения разреза, отвечающие принципу условно-корректной задачи, но не обеспечивайте геологически содержательных результатов. Это объясняет бесперспективность поиска оптимального решения в рамках одиночного зондирования и приводит к выводу о том, что критерий отбора множества корректности помимо требований, свяванных с математической устойчивостью решения, должен обеспечивать получение пространственно согласованной, геологически содержательной модели среды [4,52].

Этот интегральный принцип, на над ввгляд, является из

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

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

Бессмысленно искать множество корректности не определив, что следует пенимать под модельным решением задачи. Это понятие не яз-ляется однозначным. Теоретически оно определяется теоремой единственности [Langer,1933; Тихонов,1934; Друскин, 1982]. Однако в условиях практической недоопределенности его смысл меняется. Под модельным решением в этом случае будем понимать некоторый предельно ограниченный, максимально достижимый для данного набора наблюденных, априорных и апостериорных данных класс решений. По мере доопределения обратной задачи результатами других методов: электрического каротала скважин, геологической съемки и др. этот класс сузха-ется в соответствии с информативностью дополнительных данных.

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

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

Попытки формального использования разработанного аппарата численного моделирования сложнопостроенных сред для репения обратных задач нэ приводят пока к заметному повышению результативности метода. Причинами этого, наряду с рассмотренными выао, являются и такие, как отсутствие характерных, однозначных признаков распознавания моделей сред по внешним особенностям наблюденных полей и существующая при определенных условиях неразличимость моделей, являющаяся следствием неоднозначности решений в разных модельных классах СТопфер,1972; Szaranlek,1982; Schulz,Tezkan,1888 и др.}. Пересечение множеств решений, относящихся к равным модельным классам, в пределах которого эти модели неразличимы вследствие практической

недоопределенности исходных данных, будем называть множеством модельной эквивалентности С521.

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

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

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

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

Все вышесказанное определило цель проводимых исследований и суть развиваемой концепции, которая состоит в том, что в отличяе от традиционных методов интерпретации одиночных .зондирований, предлагается комплексная интерпретация совокупности аондирований, полученных по площади или профилю, в режиме "модель-поле". При .атом процесс интерпретации разбивается на два отдельных, но тесно вэакжсвяЕашшх этапа, отличающихся по своей природе и особенностям проявления некорректности: физический и геологический. Это составляет суть первого защдасгмого положения [4,18,19,25,30,35, 52].

Практическая реализация предложенной • концепции базируется на

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

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

" Вопросы распознавания ' и подавления различного рода помех изучались в работах Л.И.Альпииа (1945); Н.Г.Шкзбарпи (1990); В.Я.Колесникова (1976,1981,1993); В.В.Кусксза.Д.Г.Яксплезс (1С8й); О.Ку-фуда (1984) и др. , ЕЗирокио бозмояности в ресении данной проблемы появились 'в связи с разработкой программ численного моделирования в двухмерно- и трехмерно-неоднородных средах. Нгаболее вяачиш». результаты здесь получены благодаря >,сследовани*;ч В.Я.Дмитриева (1973,1986); В.В.Кускова (1979,1985,1989); А.Г.Яковлева (1039); Н.Н.Серебренниковой (1987) и др.

Физическое объяснение и мздгнигм обраговапш болыгтаства вндсэ искажений связаны с галыжмческими эффектами пере<>асг.р.¿деления плотности электрического тока в лротодяг.ей среде, обусловленного теми или ¡жми ocoienHocTPv.il его строения.

На моделях вертикально-неоднородной среды, среды с поверхностным полусферическим включением и однородной среда с включением двумерного проводящего тела заложено математическое моделкрозаччо с цель» количественной оценют влияния на результаты аондировакчя кеоднородностей разреза. Анализ и обобщение полученных результатов с учетом данных друпгх исследователей, заншавгяхса .изучением влияния глубинных двумерных и трехмерных кеоднородностей ГСаков-цев,1959; Сапожников, 1973; Кус:хв, Ячозлев, Ыодня,1989,1692 и др.], пгагззалл, что наибольшие искажения связали с влиянием верхней неоднородной части разреза. Эти неоднородности, з частности, являются основной причиной образования "ворот" на кривых зондирования, возникающих при изменении длины приемной линии. Наибольшие искажения связаны с линейно вытянутыми проводниками, которае, концентрируя в себе повышенную плотность зарядов и перенося их на больше расстояния, способны более еначительно ослаблять платность тока в окружающей среде по сравнению с локальными -голами, для которых, вследствие гальванического эффекта обтекания, эти акомалъ-

ные проявления менее значительны.

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

Рассмотрели приемы нормализации - устранения вертиг-альных смэ-цевий кривых зондирования и их отдельных участков, полученных при разных длинах приемной линии (Ш), связанных с влиянием горизонтальной неоднородности среды; показано, что вопрос о выборе уровня приведения решается в.зависимости от известных сведений о строении научаемого разреза: при сравнительно однородной верхней части раз-рева (речные, щзскиэ и т.п. зондирования) им шгут служить участки кривых зондирования, полученные при минимальных Ш, а при устойчив!« электрически свойствах опорного горизонта - соответственно при максимальных Ш1. В случае отсутствия таких сведений предпочтение, следует отдавать участкам кривой, соотвзтствугвдкм наибольшим разносам питающей лиши, так как в силу интегрального характера метода вероятность "выскоков" в отом случае меньше, чем для малых разносов. Задача нормализации решается путем минимизации соответствующего (оплакивающего) функционала, учитывавшего рассмотренное особенности и использующего информацию о характере и размерах перекрытий для каядой кривой вотирования.

На модели линейного проводника проведет катематические и экс-перкмектальные исследования влияния техногенных поиех типа трубопровод, кабель и т.п. Показано, что аномальные проявления связанные с влиянием таких объектов, наблюдаются в широком диапазоне разносов и могут достигать 49 % относительно фоновых значений. Эжлюри-кеоташа« наблюдения на раэии типах трубопроводов (нефтяют, г'а-еотзих, водопроводах, магистральных к местного сна-гения!, в целом подтверждая результаты математического моделирования, выявили ряд особенностей, в частности, связанных с вертикальной 'неоднородностью среды и состоянием антикоррозийного покрытия. Установлено, что при хорошем состоянии антоирроэкйного покрытия, являющегося каолятором для постоянного электрического тока,, гальванические эффект;: перетекания тока отсутствуют и влияние трубопровода при глубине его залегания 1,6 м и более практически не сказывается на результатах наблюдений.

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

С1,15,41], основаннннх на совместном анализе совокупности зондирований для исследуемого участка, и, в частности, способа объемной оптимизации наблюденного поля, использующего в качестве анализиру- ■ емой функции 1саяущейся проводимости 5(г,х,у)-г/ру.(г,х,у). Из физических'представлений [Матвеев,1674,1000; Хмелевской, 1984] поле значений 5(г,х,у) в случае квазигсризоптально-слоистых сред можно представить в виде совокупность взаимосвязанных непересекащдхся поверхностей, соответствующих какдоуу из разносов г. Отсюда задача выделения регулярной части поля сводится к согласованному послойному сглгзкиванию значений Э(г,х,у) в пределах заданного уровня погрешностей.

. Решение задачи ищется путем минимизации следующего функгиона-

ла:

а.....■»-АА[ахг!Ч

2

•»11

й

+ | £ фгас! Sj.fi - йтай ^ сМу | ,

■ ¡1

где 1, 1 - соответственно искомые и экспериментальные значения функции кажущейся проводимости для -го ра&носа Од =3 (ХьУ!));

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

Первый член функционала мшчмиэирует среднегаадратйческое отклонение исковой зависимости от оютёрзадеатальной, второй - обеспечивает условие гладкости функции 3 (х,у), третий сглагтаавт различие в геометрии между соседними поверхностями или, »паап словаг ми, максимизирует'их подобие.

В процессе миниаация функционала подбирается такое значение й, которое обеспечивает выполнение условия

где 5 - заданная степень свободы изменения S, определяемая погрешностью наблюдении.

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

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

Визуализируемыми параметрами могут служить наблюденные параметры поля, физические характеристики среды и различные их трансформанты, подчеркивающее те или иные особенности разреза, либо усиливающие проявление скрытой или неявно выраженной информации. По традиции этап истолкования, связанный с анализом наблюденных параметров поля, будем наеывать качественной интерпретацией. Вопросам разработки приемов качественной интерпретации посвящено много работ (Каленов,1957; Киричек, 1969; Матвеев, 1974,1975; Одеков и др.,1976; Колесников,1975,1981; Хмелевской и др.1988; Р.Венца-лек, 1991; и др.1.

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

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

2. Способы, основанные на прослеживании и анализе отдельных частей геоэлектрического разреза по характерным точкам и участкам наблюденных кривых B3Q (карты абсцисс и ординат экстремальных точек pk(r), равных значений кажущихся сопротивлений для определенных разносов, углов наклона правых асимптот, суммарной продольной проводимости и т.п.).

3. Способы, основанные на анализе дифференциальных эффективных параметров (разрезы, графики, карты, блоковые отображения дифференциальных кажущихся сопротивлении (проводимостей), нормированных производных и т.п.).

Рассмотрен физический смысл и приемы анализа эффективных и ко-ииественнкх характеристик разреза, позволяющие визуализировать пространственное поведение около двадцати различных вмдов графического их отображения. При визуализации результатов количественной интерпретации та значения могут быть вычислены для любой ездите:") части разреза, что необходимо при сопоставлении электроразведочных данных с результатами других методов в определенном интервале глубин [24,32,333, а также при объемных отображениях поля либо его сканировании.

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

В обцем случае задача визуализации двух- и трехмерных изображений пространственного поведения параметров среды включает разработку способа аппроксимации дискретно заданной функции Ъ (х,у) на равномерную прямоугольную сеть и визуализацию полученной сеточной функции ¿и, 1=1 ,НХ, и различных ее проекций с помогу»)

построения изолиннй, областей либо поверхностей на экране дисплея.

При алпропсжацил исходные данные рассматриваются как потенциальные нсточгшот (Аронов, 1990), принадлежащие некоторому потенциальному полю (что достигается использованием гармонических функций) , знчислегтя сводятся к репешго краевой задачи, Дг 2(хп,уп)-5(х-хп,у-уп) » где 5(1)— дельта-функция прямоугольной области с однородны?.!!! граничными условиями типа б1/дп « 0 . Расчеты ведутся численно-итерационными методами и при большом числе исходных точек Ь имеют быстроте сходимость.

Разработанные алгоритмы и программы получения одно-, дгуя- и трехнгрккх изображений, аксонометрии и скзнкровачия поля обеспечи-вкот Еогчскность_ оперативного анализа пространственного изменения пг-раметроз к свойств среди. Они содержат специальные процедуры управления и редакции графических изображений, включая: 1) фильтрацию помех в рамках взданного "окна"; 2) выделение и локализация аяомальных зон, .нанесение текстов и условных обозначений; Э) наполнение геологически символикой (цвет, геологичесгсач легенда) выделенных пачек пород; 4) хранение копий изображении и их совместный анализ; Б) получение твердых копий графических ивобракеиий (графиков, разрезов и карт) в заданном масштабе как с помощью плоттера, так к стандартного принтера и до

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

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

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

• Впервые идея численной беспараметрической интерпретации, явив- -шаяся следствием доказанной теоремы о единственности решения обратной эадачи ВЭЗ (Langer, 1933), была высказана Л.Б.Слихтером еще в 1933 году. Графо-аналитический ее вариант для горизонтально-ело- • истых сред рассмотрен Пекерисом CPekeris,1940]. Наиболее интенсивно идеи численной интерпретации стали развиваться с конца семидесятых годов в сваей с широким внедрением ЭВМ. Большой вклад в ре- * .шение этих вопросов внесли Л.Л.Ваньян, В.Н.Страхов, Н.Г.Шкабарня, Б.К.Матвеев, О.Куфуд и др.

В настоящее время известно несколько подходов к реализации ■ этих способов: численный, основанный на послойном анализе трансформанты кажущегося сопротивления R(m) (Н.Г.Шкабарня, Т.К.Кунички-на,1965,1971; В.Н.Страхов (1966); O.Koefoed (1966,1968,1984); * В.П.Романов (1969); В.И.Гудзь (1972,1976); В.К.Матвеев (1974); В.П.Колесников (1976,1976,1989); Т,Lee (1981); A.Kesli (1983); Лам Куанг Тхиеп, Ле Вует Зы Хыонг (1984) и др.), способы анализа функций "мнимых источников" (Kunetz, Rocroi,1970), способы математических трансформаций (Страхов,1976; Курдюков, Соболев,1992,1993), способы инверсии наблюденных данных в искомые характеристики среды на основе эмпирических (Шемэурин, Уваров, 1970; Варламов, Лобанов, 1981; Дианова, 1984) и аналитико-корреляционных (Колесников.-. 1989, 1994) соотношений.

В диссертационной работе рассматриваются численный £8,10,16, 31,39}, аналигико-корреляционный [45,46,51 1 и созданный на их основе статистический способы интерпретации 14,26,34].

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

Проведенные автором исследования,- в частности, по изучению механизма послойного образования, трансформации и распределения погрешностей в функциональном диапазоне анализируемых функций, позволили получить аналитические выражения, характеризующие зависимость погрешностей от основных вычислительных параметров 18,39]. С помощью этих соотношении установлено, что результирующая погрешность состоит из нескольких компонентов, обусловливающих 'закономерные, но довольно ревкие и значительные их изменения в диапазоне заданного аргумента т. Найдены сравнительно простые критерии и способы определения диапазонов т, обеспечивающих минимальную погрешность вычисления параметров, а также весьма полезные для интерпретации диапазоны эффективного проявления каждого слоя [13. Вместе с тем неучены свойства ряда информативных функций и диапазоны их значений, соответствующие веданной области существования физического решения, получены условия, характеризующие величину соотношения "сигнал-помеха", существенно модифицированы приемы послойного анализа, выполнен расчет новых фильтров для пересчета кажущихся сопротивлений в функцию (?(п0, отличающихся от известных более высокой точностью, особенно в случае повышенной контрастности слоев по удельному электрическому сопротивлению.

Дальнейшее совершенствование способа привело и необходимости разработки специального способа коррекции [10,393, реализующего итерационный' процесс контроля и уточнения результатов каждого из этапов послойного анализа, особенностью которого является то, что в каждом цикле уточняются не^только параметры, но и число слоев, а также интервалы их эффективного проявления в функциональном диапазоне кажущихся сопротивлений. В итоге бьш соэдан по существу новый способ-численной интерпретации; отличающийся..повышенной устойчивостью результатов при полной автоматизации процесса.

Результаты математического моделирования покавалн, что при отсутствии погрешностей наблюдения и оптимальном шаге дискретизации анализируемой функции 1?(т) способ позволяет восстанавливать иссле-

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

Вместе с тем теоретические и зксперименталъше оценки способа, в условиях кедоопределенности исходных данных выявили ограниченные возможности его практического применения, ся;шашше о опредеденк*-ми требованиями к экспериментальному материалу и ограничениями, накладываемыми на параметры среды: наличия хороигс выракеншм асимптотических ветвей кривых зондирования, отсутствия реэзок искажений эгяперимектального ыатериача, умеренной дкйдреициадаи слоев по сопротивлению ( 0.06< рц-1/р1 < 15) и др. Кроме того, что немаловажно, способ крайне трудоемок в его реализации. Эта неудовлетворенность и накопленный опыт в области интерпретации ылектрк-чгсккх вондировгшиа привели к необходимости понагз более простых способов, адекватных информативности метода.

В результате разрзбот-зк ноьий получивший название аналотк-ко-корреляционного способ рэвенин обратной задачи, основанный на кспсльтоваилн выявленных зпристическнх корреляционно-аналиткчес!«« связей м<5»ду параметра;«! среды и характерны».« особенностями их проявления на кривых зондирования [45,45,47,513.

Осноеой для его создании послужит результаты математического моделирования, позволившие- установить области существсвлакя реые-нж; и получить уравнения свяеи меэду координатами характерных точек кривых зондирования и нскотк: параметра«« среды. Устойчивость и надежность результатов обеспечиваются: а) поиском решения в ограниченном множестве, определенном б результате математического коделкроаедия для заданного модельного класса; 'б) спецкааыййм приемами коррекции результатов, основсанымк на испо-.ъзоааник цо»с-флцированного способа, предстс^лякщего комбинацию способов Ньютона и покоординатного спуска, с оцгкпой аппроксуэ>'ации ь ьлаяаэсье аффективного влкянел соответствующего слоя; в) модельной фильтрацией, педволг.шеи выделять регулярную часть анализируемой фушади для заданной модели среды и др.

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

Создание прямых численных способов, осуществляют,;« детерминирован:^ связь искомых характеристик среды с наблюденными парамет-рамл поля, дает возможность объективного реиемия кас данной задачи, т?к и исследования различных факторов, влияющих на результата зондирований. Разработанный на их основе способ статистичес:юго моделирования позволяет выделять полный набор решений, отзечгзсщих заданному урогию погрешностей функции южущегося сопротивления, и ка количественном уровне анализировать особенности формирования обобщенных геоэлектрических моделей сред, пределы неоднозначности, эффективность тех или иных приемов адаптационной регуляризации, наличие и тесноту связей медцу параметрами и др. С помощью данного способа получены новые результаты относительно эквивалентности решений; экспериментально обоснована многомерная некорректность (относительно числа слоев, их параметров и вариантов формирования обобщенных слоев); установлены фильтрационные свойства моделей и их влияние на оценку аппроксимации наблюденных данных рассчеташи, похаваны возможности повышения детальности расчленения ровреэа в случае слабоконтрастных по сопротивлению слоев и др.

Результаты четвертой и второй глав определяют содерчетие второго защищаемого положения 11,3,8,9,10,11,12,14,16,16,17,25,31, 38,39,41,44,45,47,43 3.

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

В силу неоднозначного соответствия меяду геологическими и физическими характеристиками пород, модельной згаивалентнссти набл»-демных электрических полей, индивидуальных неформализуешх особенностей каждого геологического об?>екта, построение фивико-геологической модели может рассматриваться яак вада*:а выбора каибодзе вероятной гипотезы строения среды в классе вовмссяшх физических решений Е191.

В общем случае эта задача включает: а) классификацию элементов, выделенных на этапе послойного расчленения разрева, по их принадлежности к разным модельным классам (кзагигорвэоэтальио-сло-истой модели, локальным неоднородноетжи и помехам); б) стратиграфическую привязку и оптимизацию геоэлвктричэсхой модели ср-едо £4,3&,42).

Задачу первого из этих этапов можно сформулировать следующим-образом. Пусть в результате физического этапа интерпретации в каждой точке зондирования получены решения, представляющие один из эквивалентных и вэаимонезависимых вариантов расчленения разреза, отображающих проявление неизвестного числа геозлектрических слоев, возможных локальных неоднородностей и помех, которые будем расс- . матривать как классы [353. Каидый класс характеризуется вектором изучаемых признаков Хц ( 1 - 1,2,...,т), где и - число признаков, 1 - номер класса. -

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

Число классов и их статистические характеристики обычно неизвестны. Такая постановка задачи в теории вероятности приводит к задаче распознавания образов на принципах самообучения.

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

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

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

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

метров и свойств среды, представляющие апостериорную информацию о геоэлектрической модели среды.

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

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

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

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

2) анализ параметрических зондирований, выполняемый с целью геологической привязки границ физической модели и определения коэффициентов анизотропии слоев;

3) трансформацию полученных в параметрических точках коэффициентов анизотропии на всю площадь исследования и получение уточненных параметров среды, согласованных с априорными представлениями о разрезе.

На основе корреляционного анализа зависимостей между электрическими характеристиками разреза и другими физическими свойствами среды (скоростью, плотностью, минерализацией и др.) (Серов, 1972; Тульбович,1979; Колесников,1983,1993 и др.) параметры геоэлектри-ческиого разреза могут наполняться соответствующим содержанием.

В итоге геологического истолкования формируется геоэлеотри-• ческая модель, едекватно отвечающая информативности метода и имеющемуся набору априорных данных.

Содержание пятой главы раскрывает суть третьего защищаемого положения [1,4,6,6,18,19,28,35,36,48, 523. -

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

Система гСГО представляет собой автоматизированное рабочее место июязиера-геофизика, располагающее необходимым арсеналом методов и средств для обработки, интерпретации, визуализации к графического представления результатов вертикального электрического-зондирования [1,4,38,44,52).

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

Математическое обеспечение системы гСШ включает более 50 отдельных программных модулей, реализующих алгоритмы различных функциональных назначений, объединенных в.семь основных блоков обработки и интерпретации (рис.1):

1. Контроль, визуализация и корректировка исходных данных.

• 2. Первичная обработка и предштерпретационная регуляризация ревультатов наблюдений.

3. Качественная интерпретация.

4. Количественная интерпретация (физический этап).

5. Геологическое истолкование результатов физического решения.

6.' Расчет геоэлектрических параметров для заданной части геоэлектрического разреза. 4 ~

7. Исследование.корреляционных связей с заданными характеристиками среды (минерализацией, скоростью, плотностью и др.) с целью решения прикладных задач.

6. Моделирование обратной задачи электрического зондирования.

Модули блока ДАННЫЕ осуществляют различные операции над зе-зультатами наблюдений: формирование банков данных, сортировку по заданным параметрам' (числу разносов, номеру профиля, и т.п.), контроль создания файлов исходной инфозмации, визуализацию', редактирование графиков зондирования и др.

Блок ОБРАБОТКА включает анализ перекрытий участков кривых -зон-

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

В блоке КАЧЕСТВЕННАЯ ИНТЕРПРЕТАЦИЯ реализованы алгоритмы расчета и визуализации разрезов, карт, графиков и объемных отображений различных эффективных параметров, включая: кажущиеся сопротивления, кажущиеся проводимости, их дифференциальные и интегральные трансформации, нормированные производные. Предусмотрены процедуры фильтрации поля в рамках задаваемого окна.

Программы блока КОЛИЧЕСТВЕННАЯ ИНТЕРПРЕТАЦИЯ реализуют алгоритмы решения обратной задачи и статистической оценки пределов неоднозначности решений, выполняемых как в автоматическом, так и в интерактивном режимах.-

Блок ГЕОЛОГИЧЕСКОЕ ИСТОЛКОВАНИЕ предназначен для выделения регулярной части разреза (квазигоризонтально-слоистой модели среды), участков проявления локальных неоднородностёй и помех.

Технологически построение фивико-геологической модели представляет итерационный процесс, выполняемый при участии интерпретатора. Работа интерпретатора сводится к заданию управляющих параметров, представляющих ограничения на пространственное изменение физических параметров среды и на пороговые значения, определяющие принадлежность результатов решения к одному из ваданных модельных классов, а также к выбору наиболее оптимального варианта конструируемой модели. На этом этапе при необходимости используется блок моделирования решения обратной задачи и повторная количественная интерпретация отдельных зондирований, не отвечающих регулярной части модели, с целью проверки наличия в рамках эквивалентности других наиболее вероятных вариантов расчленений разреэа. Аномальные эффекты, не укладывающиеся п ргмгси квазигоркБонтально-слоистой модели среды и отображаемые в виде поля остаточных аномалий, классифицируются как лекальные неоднородности либо помехи по совокупности признаков (величине коэффициентов прослеживаемости, пространственной согласованности, достоверности результатов).

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

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

Геофизические измерения ■ с гги^рсьзй регистрацией

База даягодс

Ввод данных в ЭВМ к контроль:

Визуализация к редахгдо ва&подекных [ дыягых |

Предылтерзретафюняал

1*?гул«ризация: - ворколкэацкя; • объеюия овткьою^&м

Визуализация и люлиа; рк <Г), ПОЛуч^аЗДЕЯ ДОЛ \

разных приемных ллжЛ

- гмегогра-*№$ :

.........................

Качеством дал юггерпретац;** /расчет и анализ эффективных параметров

поля/ Рк, брц, б5к, Хк и др.

Моделирование |

т

Еиауализадия и редахция пол/:

- скиыровдие;

- шсонэыетрня; -ка рта;

- разрезы

Количественная интерпретации /физический этап/

- численный способ;

• авалмтихо - хо£едяциовяыЯ способ;

- етйгмегкчвский спост-Й

I

/упорные

__дя/ггше

>

Геслоюгчесгое мстоЛсэзашиг: - фмзкю - гтхиктэтеехдч код^ль; - кор^кдефотсаьЛ иаалн4

РЕЖИМ: |—автоматический •7.7.'.'. ин?асрйктоие«гш

Визуализация и ; редахция результата ; интерпретации

; сдааочеых ЗОЕДИр^аЗЕИЙ

; Вкз^влиаация и

• р^дашфш в ^яуси» I "модель •• псле"

• - ска^р-. илйане;

• - аихоясмефия;

- карги; ■ -

[олучеане жестких кс>Ш1й, ^едадддгл и о фор млеете

Рис. 1. Технологическая схот, автоматизированной обработки и интерпретатга оевультагов электрического зондирования

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

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

Разработка и совершенствование системы автоматизированной интерпретации электрических зондирований проводились в процессе непрерывного, тесного контакта с прои8годственнш!М организациями (Территориально геологическими объединениями "Уралгеэлогия", "Са-халингесяогуя", трестом "Пермнефтегеофивика", Челябинским геолко-мем и др.) в условиях непрерывно совервенствукщегося парка ЭЗМ, с привлечением появлкотдася новых программных средств и технических возможностей [12,13,14,20,21,27,29,31,32,43,483. С конца семкдеся-тых годов с помогаю системы обработаны и проинтерпретированы десятки тысяч зондирований при решении различных геологоразведочных задач в производственных условиях.

Содержание пестой и отчасти третьей глав раскрывает суть четвертого защищаемого положения [1,3,4,12,14,16,17,18,19,27,28,90, 31,35,38,44,52).

- , ЗАКЛЮЧЕНИЕ

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

1. Предложена и научно обоснована концепция согласованной интерпретации совокупности электрических вондирований с разделением процесс» уа два этапа: физического и геологического, вовлекаемых в

итерационный процесс на интегральном уровне "модель-поле". В связи с этим введены и даны определения основопологахндих принципов и понятий, учитывающих физике-геологические особенности ' метода, а именно: понятия ква8игоризонтально-слоистой модели среды, интегрального принципа интерпретации, модельной эквивалентности; уточнены понятия множества корректности, точного модельного решения задачи, принципа регуляризации, . что играет важную роль в вопросах методологического обоснования компьютерной технологии.

2. Разработаны способы решения,обратной задачи электрического зондирования, не требующие предварительного задания начальных приближений конструируемой модели среды: прямой численный способ, аналитико-корреляционный и статистический. Устойчивость и точность решения обеспечиваются комплексом приемов адаптационной регуляризации и численного анализа, включающим: а) поиск решения на ограниченном множестве, определенном в результате математического моделирования для заданного модельного класса; б) использование диапазонов эффективного влияния каждого сдоя и критериев оптимальности решений с учетом соотношения "сигнал-помеха"; в) повышение точности фильтров трансформации рк(г) в И(гп); г) применение специальных приемов контроля и коррекции результатов на основе взаимосвязанного решения прямой и обратной задач; д) предынтерпретационную регуляризацию измеренных величин; е) модельную фильтрацию; ж) использование интегрального принципа при выделении регулярной части наблюденного поля и изучаемой модели среды и др. С помощью численного способа .экспериментально подтверждены основные положения теоремы о единственности решения обратной задачи ВЭЗ для слоистых горизонтально-однородных сред.

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

тнвлению слоев и др.

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

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

6. Разработаны и реализованы принципы автоматизации геологического этапа истолкования результатов электр!1ческого зондирования, включающие: а) ввделение квазигоризонталько-слоистон модели среды, областей влияния локальных неоднородностей и помех по результатам интерпретации одиночных зондирований на основе интег-ралького принципа интерпретации, методов распознавания образов и принципов самообучения; б) анализ параметрических зондировании, с цель» изучения особенностей формирования физико-геологических моделей; в) оптимизацию пространственного поведения физических параметров с учетом информативности наблюденных данных и априорных представлений о строении среды. Предложен методологический подход к изучении слокнопостроенных сред, основанный на поникении модельной размернссти при решении обратной задачи, путем последовательного распознавания модельных классов по признаку их аномальных отклонений относительно моделей более низкого порядка. Развитие этого направления может служить одни! из альтернативных вариантов проблемы интерпретации в условиях горизонтально-неоднородных сред.

7. Разработана автоматизированная система и высокопроизводительна? компьютерная технология обработки и интерпретации дйнтгх электрического зондирования, а также отдельных модификаций метода сопротивлений (электрического профилирования, срединного градиента

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

Использование системы при решении структурных, гидрогеологических, ив*.ено-геологичес!;их, экологических и других задач подтверждает высокую геологическую эффективность ее практического применения в условиях гетерогенных квазигоризонтально-слоистых сред.

Основные положения диссертации опубликованы в следующих работах :

монографии и справочники

1. Обработка и интерпретация результатов вертикального электрического зондирования с помощью ЭВМ. М. "Недра", 1981, 140 с.

2. В монографии В.К.Матвеева "Интерпретация электромагнитных зондирований". - М.,"Нздрз" 1S74, с. 195-198, 204-207.

3. Обработка и интерпретация кригых ВЭЗ с помоиуью ЭВМ. "Электроразведка" - справочник геофизика, М. "Недра", 1980, с. 106-110 (совместно с В.К.Матвеевым).

4. Система програл! автоматической обработки и интерпретации данных ВЭЗ - "Зонд" . В кн."Электроразведка" . Справочник геофизики. , М.,"Недра", 1989, т.1, с. 163-168.

статьи и те&исн докладов

5. Машинная интерпретация параметрических зондирований с использованием данных электрического каротажа. //Вопросы обработки и интерпретации геофизических наблюдений. - Пермь, 1974. Вып. 12. С. 98-102.

6. Привязка опорного ялектрического горизонта на основе расчета синтетических кривых ВЭЗ. //Геофизические изыскания. - Перка, 1S75. Вып'Л. С. 93-100.

7. Качественная - интерпретация результатов ВЗЗ на ЭВМ. //Геофизические изыскания. - Пермь, 1975. Вып.1'. С. 83-92. (Совместно с Б.К.Матвеевым, Г.А.Акопяном, Л.П.Сосниной).

8. Исследование ошибок j/гдинной интерпретации результатов ВЭЗ с использованием данных электрического каротака. //Вопросы обработки я интерпретации геофизических наблюдений. - Пермь, 1975'. Вып. 13. С. 112-116.

9. О восстановлении геоэлектричесгаго paspesa по донным ВЮ. Деп. В ВИНИТИ. 1976, N 2043-76. С. 1-6.

10. Способ коррекции в методе численного решения обратной задачи ВЭЭ. //Геофизические изыскания; - Пермь, 1976. Вып.2. С. 105-103.

11. Исследование факторов образовали! разрывов экспериментальных кривых ВЗЗ. //Геофизические изыскания. - Пермь, 1976. йш.2. С. 98-104.

12. Система автоматической интерпретации результатов вертикального электрического вондирозания. Тезисы докл.У-й Уральской конф. молодых геол. и геофиз . Свердловск, 1976. С. 67-69.

13. Опыт интерпретации результатов вертикального электрического зондирования на ЭВМ. //Вопросы обработки и интерпретации геофизических аномалий. - Пермь, 1977. Межвед. сб. иаучн. трудов. С.96-103.

14. Количественная интерпретация результатов вертикального электрического зондирования на ЗВМ. Сб. "Алгоритмы и программа для ЭВМ V-220M". Деп. в БИШТИ, 1978, !) 5.303-78. 15 С. (СоЕкестш с В.К.Матвеевым, Л.П.Сосни.'.ой).

15.. Алгоритм регуляризации результатов наблюдения метода вертикального зле)стркчэского зондирования. Сб. "Алгоритма и программы для ЗВМ V-220M". Деп. в ВИНИТИ, 1978, N 1803-78. 7 С. (Совместно о Н.Н.Серебренниковой).

16: 0 построении алгоритма автоматической интерпретации данных МТЗ. //Электромагн.' зондирования. Тез. докл.У-й Всессюзн. шгсо-лы-семин.- Киев, "Наукова думка", 1978. С.41-42. (Совместно с

B.К.Матвеевым).

17. Интерпретация дачных'вертикального электрического зондирования с помощью ЭВМ. Тезисы докл.обл. научно-техн.конф. Пермь, 1980.

C.62-63. (Совместно с Н.Н.Серебренниковой).

18. Площадная интерпретация данных вертикального электрического зондирования на ЭВМ. //Электромагн. зондирования. Тез. докл. У1-й Всесоювн. школы-семин.по электромагн. эонд-ям.М. МГУ, 1981. С.108. (Совместно с Б.К.Матвеевым).'

19: Автоматизация процесса истолкования формальных решений в методе ВЭЗ. //Геофиз. м-ды поиск, и разв. местор. нефти и газа. -Пермь. 1981. С. 96-100. (Совместно с Н.Н.Серебренниковой ).

20. Опыт интерпретации данных ВЭЗ с помощью ЭВМ при решении гидрогеологических задач. //Примен.матем. методов и ЗВМ при обработке

информ. на геологоразв. работах. Тез. докл. Уральской конф. -Свердловск. 1982. (Совместно с Н.Н.Серебренниковой,, 3.И.Лузиной, B.W.Соловьевым).

21. Выявление аномальных эффектов от 8алежей углеводородов по данным электромагнитных вондирований. //Трассир.геофиз. м-ми вон, перспект-х на литолого-стратигр. залежи ■нефти и газа. Тев.докл.и.-техн. конф. - Пермь, 1983. С.14-15. (Совместно с Б.К.Матвеевым, В.А.Поносовым, Ю.И. Степановым).

22. Математическое моделирование электромагнитных полей с целью поисков месторождений нефти и газа в северной части Тимано-Печерс-кой провинции. //Минеральные ресурсы Западного Урала и их народно-хоз. значение. Тез.докл.н.-техн. совещания - Пермь, 1983. С. 37-38. (Совместно с В.Л.Третьяковым).

23. О кинематическом способе интерпретации результатов электромагнитного зондирования становлением поля. //Теофизич. м-ды поисков и разведки нефти и газа. - Пермь, 1983, с. 20-27 (в соавторстве с В.К.Матвеев, Е.И.Ларионовым).

24. Определение скоростных параметров верхней части разреза по' • - данным электрического вондирования с помощью ЭВМ. //Изучение ВЧР в '

сейсморазведке МОГТ.- Деп.в ВИНИТИ. 1983, N 1245-83. С. 99-123.

25. О путях повышения однозначности истолкования данных электрического вондирования. //Методы геологических исследований. 'Тез. докл.н.-техн. совещания - Пермь. 1984. С. 100-102.

26. Алгоритм интерпретации данных электрического зондирования на основе статистического моделирования. Сб."Автоматизация научных исследований и обучения". Деп. в ВИНИТИ, N 2344-84, 1984, с. 47-53 (в соавторстве с В,А.Кутиным, Л.П.Сосниной).

27. Опыт автоматизации истолкования данных ВЭЗ при решении 1 гидрогеологических задач. //Геофиа.мет. в'гидрогеол.,инж.геол., гидротехнике. Tee. докл. Э-го науч.-техн. семин.-совещ. - Ереван, 1985. С. 62-63.

28. Вопросы автоматизации истолкования данных электрического зондирования. //Примен. мат: методов и ЭВМ при обраб. инф.на геол. -ра8в. работах. Тез.Докп.5-й Уральской конф.- Свердловск. 1985. С. 126-127.

29. Опыт автоматической интерпретации данных электрического зондирования при решении геолого-поисковых задач. //Примен. мат. методов и. ЭВМ при обраб. инф.на геол.-равв. работах. Тез.докп.5-й

Уральской конф. - Свердловск. 1985. С. 127-129. (Совместно с Соловьевым В.М., Сосниной Л.П. ).

30. Ó методологии интерпретации электромагнитных зондирований. //Автоматизация приемов обр.-ки геоф. информ.при поисках нефти и газа. Тез. докл.науч.-техн.совещания. - Пермь. 1985. С. 5-7. (Сов-' местно с Б.К.Матвеевым).

31. Автоматизированная система обработки и интерпретации данных ВЭЭ." //Автоматизация приемов обр.-ки геоф. информ.при поисках нефти и газа. Тез. докл.науч.-техн.совещания. - Пермь. 1986. С.31-33.

32. Опыт использования автоматизированной системы интерпретации ВЭЗ при изучении ВЧР. //Автоматизация приемов обр.-ки геоф. информ.при поисках нефти и гава. Тез. докл.науч.-техн.совещания. 1986. С. 75-76. (Совместно с J0.И.Степановым).

33. О методике получения статических поправок по данным электрического зондирования. //Изучение и учет верхней части раареза при геофяв. работах на нефть и газ. Тев.докл.научно-техн.конф. -Пермь. 1985. С. 9-10. (Совместно с Б.К.Матвеевым,"Ю.И.Степановым).

34. Использование статистического моделирования при интерпретации данных B3G. Деп. з ШШТЛ, 1986, N 121-В. 15 с.

35. Приемы автоматизации геологического этапа истолкования данных вертикального электрического зондирования. Деп. в ВИНИТИ, 1986, N 120-В. 11 с. (Совместно о В.А.Пугиным).

36. Использование приемов самообучения при автоматизации • истолкования данных электрического зондирования. //Нов.методы поиск, и раав. и анализа мест. пол. ископ. в связи с коми.ивуч.недр Зап. Урала. Тез.докл. научн.-техн.совещ. - Пермь. 1987. С. 87-88.

37. 0 применении геофизических исследований на Верхнекамском месторождении калийных солей. //Нов.методы поиск, разв. и анализа мест. пол. ископ. в связи с компл.изуч.недр Зап. Урала. Тез.докл. научн.-техн.совещ. - Березники. 1987. С. 74-75. (Совместно с

B.М.Новоселицким, В.А.Поносовкм, И.Е.Королевым).

38. Автоматизированная система интерпретации электрических зондирований. //Электромагнитные зондирования. Тез. докл.11-й Всесоюзной школы-семинара, - Киев. 1987. С. 45.

39. Некоторые вопросы послойного аналива данных электрического зондирования. //Ивуч.влияния верхней части раврева на парам.фиэ. полей в услов. Пермск. Прикамья. Деп. в ВИНИТИ, 1987, M 2931-В81.

C. 43-67.

40. Приме-ненке СЗП при выявлении локальных неоднсродностей ВЧР. //Поиски и разведка геофис. методами неструктурных еалежей нефти и газа. - Тез. докл. научно-техн. конф. - Пермь. 1937. С. 15. (Совместно с D.M.Степановым).

41. Использование приемов оптимального сглаглвания экспериментальных дсгзак при автоматической интерпретации • В38. Деп. в ВИШ1-ТИ, 1987,' N 4665-В87. 13 с. (Совместно с В.А.Путиным).

42. Построение физико-геологических моделей ВЧР по данным электроразведки. //Комплексное исследование кедр Зап.Урала - путь ус-кор.разй.хоз.-ва региона . Тез. докл. - Пермь. 1988. С. 102-103. (Совместно с Б.К.Матвеевым, Ю.И.Степановым).

43. Состояние и перспективы обработки и интерпретации вертикальных электрических зондирований (В38) по система "Зонд-5". //Прямей.мат. мет. и ЭВМ при обработке- инф. на геологорззв.работах. Тез.докл.У1-Й Уральской конф. - Челябинск. 1989. С. 94-9Б. (Совместно с Л.Г.Ломовым, В.М.Соловьевым).

44. Разработка автоматизированной системы интерпретации данных электрического аондирзвакия в условиях ВКМКС. //Пробл. - компл. к?уч. водоеод. толщи на местор. кал. солей. Теь. дога. per. сов. -Пермь. 1989. С. 27-28. (Совместно с В.М. Новоселищага, И.В.Онько-вым, И.Е.Королевым, Е.И.Туляевым).

45. Зкспресс-интёрпретация данных электрического зондирования на основе исследования связей мегду эффективными и истинными параметрами среды. //Пробл. кэмпл. изуч. ьодозач- тол'ци на местор. кал. солей. Тев. докл. per. сов. - Пермь. 1&?9. С. 28-29. (Совместно с И.Е.Королевым).

46. Интерпретация электрических зондирований методом особых точек. -Деп. в ВИНИТИ, 1990, N 4959-В90. 19 с.

47. Об использовании особых точек при интерпретации данных В'ЭЗ. Ден. в ВИНИТИ, 1991, N 1487-В91. S с.

48. Ре&удьтаты электроразведочных работ на поверхности шахтного поля ВКПРУ-3 Верхнекамского месторождения калийных солей. //Комплексное1 освоение недр San'. Урала / АН СССР, УрО. - Свердловск. 1991.. С. 67-61 (в соавторство с В.А.Покосовым, И.Е.Королевым).

49. Исследование влияния техногенных помех на результаты электрического вондирования. //Горная геофизика. Tes. докл. 6 Междунар.. сеьин. - Пермь. 1993. С. 29-30. (Совместно 1 с Б.А.Кутиным, И.М.Скумбиньш).

50. Воимолвости метода ВЭЗ при решении экологических задач в условиях техногенных помех. //75 лет сельскохозяйственно^ образования на Урале. Теэ. докл. - Пермь. 1993, С. 332-113, (Совместно с И.М.Скумбиннм).

51. Лналнтико-корреляционный способ интерпретации данных ВЗВ. //Применение геофмзическкх методов при решении- геологических, инженерно-геологических V экологических задач. Тез. докл. - Пер». 1994. С. 7-8.

52. Площадная компьютерная интерпретация' электркчеаэдх зондирований. //Применение геофизических методов при решении геологических, «яженерно-геологмчесютх и экологических задач. Тез.- долл. -

Пермь. 1994. С. 9-11.