CC BY
58
ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
2011 Геология Вып. 2 (11)
УДК 550. 837
О повышении однозначности интерпретации электрических зондирований
В.П. Колесников
Пермский государственный университет, 614990, Пермь, ул. Букирева, 15 E-mail: poisk@psu.ru
(Статья поступила в редакцию 4 апреля 2011 г.)
Предложен новый способ повышения однозначности автоматической интерпретации площадных электроразведочных наблюдений методом ВЭЗ, основанный на использовании совокупности регуляризирующих элементов, вовлекаемых в интерпретационный процесс на основных его этапах. Приведены способы доопределения обратной задачи, общая схема построения алгоритма и приемы его реализации.
Ключевые слова: вертикальное электрическое зондирование, компьютерная интерпретация, регуляризация, технология.
Метод вертикального электрического зондирования (ВЭЗ) является одним из базовых методов электроразведки. В силу потребностей и интенсификации производства, специфики решаемых задач он находится в постоянном совершенствовании и развитии, требуя создания технологий, обеспечивающих оперативность, достоверность и геологическую содержательность получаемых результатов.
Наиболее сложным в развитии метода является этап количественной интерпретации, затрагивающий фундаментальные вопросы, связанные с некорректностью решения обратных задач. Ввиду некорректности наблюденным значениям электрического зондирования может соответствовать достаточно большой набор эквивалентных решений, отвечающих наблюденным данным с точностью до погрешностей полевых наблюдений 5 (в рамках 5-эквивалентности), но в то же время способных существенно отличаться от истинного.
Одним из путей повышения однозначности решения некорректных задач является их доопределение за счет привлече-
ния дополнительной информации, вовлекаемой в интерпретационный процесс на основе использования принципа регуляризации [2].
Практическая реализация данного подхода находится в стадии поиска и развития. Как показывает опыт работ [1], эффективность его применения требует разработки специальных способов и приемов регуляризации, максимально учитывающих специфику конкретного метода, его информационные возможности, особенности и обоснованность используемой априорной информации, вовлекаемой на тех или иных этапах интерпретационного процесса.
В данной работе рассмотрен один из таких подходов повышения однозначности автоматической интерпретации применительно к интерпретационной системе программ «Зонд» [1].
В случае метода электрического зондирования можно выделить ряд элементов регуляризации, способствующих повышению однозначности решения обратной задачи, с использованием различного рода дополнительной качественной и количе-
© Колесников В.П., 2011
ственной априорной информации об искомом решении:
1) согласованность решения с имеющимися параметрическими данными, например, с результатами бурения скважин:
А
)
I
Ф р,Р
■I
(2)
3) при условии гладкости границ:
„ (
/ ы =Х
п
где х
\
(3)
1л
С г =
< <
(4)
(5)
(1)
где р р - наблюденные и
рассчитанные по данным интерпретации кривые зондирования; р, , рскв - искомые
и определенные по скважинам мощности и удельные электрические сопротивления слоев р = = ; а - параметр
регуляризации;
2) поиск решения в наиболее достоверных пределах (к примеру, в пределах эквивалентности решений) для каждого ьго слоя (р,н, р,в):
При этом элементы регуляризации (1)-(5) должны использоваться при условии соответствия искомого решения точности полевых наблюдений 8 , т. е. в рамках 5-эквивалентности :
( )
<
И У
(6)
глубина до i -й границы в
точке с координатой х;
4) при условии квазигоризонтальности границ:
где к - предельное значение возможного угла наклона границы; к = 0.29 соответствует углу наклоны границы, равному 16о;
5) в рамках чувствительности кривой зондирования к параметрам 1-го слоя:
Практическая реализация отмеченных выше элементов регуляризации (1)-(6) может быть самой различной. Ею, в принципе, может служить любая совокупность вычислительных операций, нацеленных на получение решения в рамках погрешности наблюдений, наиболее согласующегося с имеющимися априорными представлениями об исследуемой модели геологической среды. Это можно рассматривать как своеобразный процесс ”закачки” различного рода дополнительной информации об искомом решении на соответствующих этапах автоматизированного интерпретационного процесса.
Ниже приводится один из вариантов его реализации с помощью алгоритма стабилизации решений. Блок-схема алгоритма приведена на рис. 1.
В качестве исходных данных на начальном этапе, помимо результатов полевых наблюдений - кривых зондирования и точности полевой сьемки (5о), предполагается использование геологических данных (результатов бурения скважин) и сведений о максимально возможной для данных условий степени изменчивости (градиента G) глубины залегания границ, а также количества вариантов (М) статистически генерируемых кривых зондирования для оценки параметров слоев и пределов их эквивалентности.
1.ІІсходные данные
1) наблюденные кривые чоццЦ}ования р^{г); і=1, Ь1; 5) -погрешность наблюдений;
2) ре чулыаты интерпретации параметрических ВЭЗ: к., р ;і= 1, п+1;і=1.Ь2;
3) управляющие параметры: О- максимальный градиент изменения глубины границ; М-количество статистически гснсріфусмьіх ВЭЗ
Г II
Ь2
2. Задание опорной модели
Опорная модель - модель, формируемая по количеству вовлекаемых в процесс стабилизации параметрических ВЭЗ (N5). К=(Ь2-Ы{) - количество контрольных ВЭЗ. N1 =1.
3. Выбор рядовой точки зондирования для уточнения решения
__________________________________1пл:: = шт_______________________________
4. Опенка пределов эквивалентности
1) внесение погрешности р’(г ) = р°(г )-(1 + <5 ); ц = 1, М; 2) количественная интерпретация набора кривых РІ(г \= при фиксированном числе слоев;
3) определение статистических параметров слоев: к , р1 , АА , Ар. .
I —
5. 'Задание параметров начального приближения
А, =?™кГ+Я%К: р<=?
- СКв СКЄ
РГ +?АДА
Т
6. Минтппаиия расхождения расчетной н наблюденной кривых зондирования
с учетом пределов эквивалентности, условия ква'лігоршонтальносіи границ, пар а ме тр а чу в с т в пт е ль н ос т н С,________________________________
І І=Ь1
7. Контроль решения н переход к следующему циклу стабилизации
8. Результаты стабилизации.
1) оптимизированная модель среды; 2) контрольные значения - графики (ц); карта 5 расхождения расчетного и наблюденного полей.
Рис.1. Блок-схема алгоритма стабилизации решения
Выполняемое на данном этапе согласование решения с параметрическими данными - один из наиболее ответственных этапов, так как связан с увязкой физической модели с реальным геологическим разрезом, включая выбор оптимального числа обобщенных геоэлектрических слоев, их литологическую и стратиграфическую привязку, выявление опорного горизонта и др., т.е. приведение априорной геологической информации в соответствие с геоэлектрической моделью исследуемой среды в рамках информационных возможностей метода.
Для реализации данного процесса создана специальная программа анализа параметрических зондирований (рис.2), которая позволяет на основе численного моделирования и визуального контроля подобрать решение в рамках 5-эквивалентности, отвечающее как инфор-
мативности метода, так и его наибольшей согласованности с геологическим разрезом. При наличии данных бурения по нескольким скважинам в результате такого анализа получается соответствующая совокупность опорных параметрических решений, являющаяся основой для формирования опорной геоэлектрической модели.
При отсутствии параметрического материала в качестве опорных зондирований может быть использована совокупность зондирований, формируемых с помощью отмеченной выше программы по результатам предварительной формальной интерпретации электрических зондирований, выбираемых автоматически по совокупности оценочных параметров (наилучшей сходимости расчетной и наблюденной кривых ВЭЗ, степени согласованности решений и др.).
Й и В р'г Я ;в Й- | Прв4иль:|1№ ^ Пш*г. П^ 3 *- ■* ЧЬ 1 [ЗА_________________________________________________________________________________________________________________________________________________________I
Я млсиялс |1 1ЙЧ 3 Г ^ШЮ№Ч4 мчри Ф * — "& ■ ■ Меси™* |Гк,<чктьч 3 Сели [Опминам 31
Рис.2. Пример визуализации результатов интерактивного анализа параметрического зондирования
Последующие блоки алгоритма (2-7) отображают автоматизированный процесс согласованной интерпретации совокупности зондирований с постепенным вовлечением опорных параметрических зондирований (формированием опорной модели).
На первом этапе из заданной совокупности параметрических ВЭЗ выбирается одна, принимаемая за опорную, остальные - используются в качестве контрольных. Процесс стабилизации решений начинается с выбора ближайшей к опорной рядовой точки зондирования. Для данной точки с помощью статистического моделирования выполняется оценка вероятностных значений и пределов эквивалентности параметров каждого ьго слоя. Данная процедура включает следующие операции:
1) получение репрезентативной выборки вариантов кривой зондирования в рамках заданной погрешности д0 путем
внесения погрешности в анализируемую кривую зондирования р (г ):
Р.
+
7=
Д - случайные числа, генерируемые с
помощью датчика случайных чисел; ]=1^; N - число разносов; 50 - средняя
относительная погрешность функции р , равная погрешности полевых наблюдений; ^ = 1, М;
2) количественная интерпретация набора сгенерированных кривых р , £ =
1, М, выполняемая при фиксированном числе слоев п, соответствующем числу слоев, полученных в результате анализа параметрических зондирований;
3) статистический анализ полученных решений и определение вероятностных значений параметров слоев:
ИА =
V
П
л
V
п
Ь - количество параметров, для которых частота г/ > г] ;
а . = ...- •... ;
А = - • .
Рис.3. Пример статистического моделирования при погрешности наблюдений до = 5 % и количестве вариантов решения М=100
Графический пример статистического моделирования при заданной погрешности д0 = 5 % и числе генерируемых кривых зондирования М=100 приведен на рис. 3.
На основе полученных оценочных статистических параметров слоев, результатов анализа параметрических зондирований и степени их весомости формируются начальные приближения Ъг-, р, необходимые для поиска решения в процессе минимизации расхождения расчетной и наблюденной кривых зондирования:
И,=дскв Нскв + д*Н,
интерпретации наблюдений опорной 8 и исправляемой 8 точек зондирования:
Ч
дА =
дскв + дА
дскв + дА
8 г
..jnp.ni птах
в
-I
Г =
+ 4* Р.
р = д р где _
а) /7/“’, р ж - средние весовые значения мощности и сопротивления ьго слоя (1=1,п; п - число слоев), вычисленные по совокупности параметрических ВЭЗ:
К‘=Ъ ; р =Е ; <; =
$= <Г=
1^; N - количество параметрических ВЭЗ;
б) дскв, д - весовые коэффициенты, учитывающие расстояние и погрешность
где д^ 8
о Пр.т ~ погрешности интерпретации опорной и исправляемой точек зондирования; г- расстояние между выбранной рядовой до ^-й опорной точками зондирования; г^п, гтх - минимальное и
максимальное расстояния между точками ВЭЗ и скважинами для данной площади;
в) /?/Л, р - значения мощности и сопротивления ьго слоя по результатам статистического моделирования для выбранной точки ВЭЗ (см. п. 4).
Процесс минимизации, осуществляемый с использованием метода Ньютона (наискорейшего спуска), выполняется с
г
учетом условия квазигоризонтальности границ (4) и заданного предела чувствительности кривой рк(г) к изменению сопротивления каждого слоя (5) в рамках вероятностных пределов эквивалентности Дhmax•
Если получаемое в результате минимизации решение отвечает условиям его достоверности: (5 < 5о ; hj > /?_ ; С,
>Со ; |Л | < А ), то данная точка ВЭЗ
используется в качестве опорной при анализе последующих рядовых точек зондирования.
После завершения интерпретации всех L1 зондирований границы разреза сглаживаются в рамках 5-эквивалентности путем минимизации функционала (3) при условии (6).
Полученное на данном этапе решение при использованной для стабилизации одной опорной параметрической точки зондирования в целях оценки его точности сравнивается с данными анализа параметрических зондирований, не использованных в процессе стабилизации (контрольных параметрических ВЭЗ) (рис.4).
Параметрическая точка с наибольшим расхождением решений, указывающим на наибольшую потребность в априорных данных на данном участке, переводится в
разряд опорных и интерпретации анализируемого массива кривых ВЭЗ повторяются с использованием для стабилизации решений уже двух параметрических зондирований. Данный процесс повторяется до тех пор, пока не будут исчерпаны все L2 параметрических точек зондирования. Получаемая при этом информация, отображающая динамику изменения величины расхождения количественных параметров в параметрических точках при соответствующей величине расхождения расчетного и наблюденного полей (рис.4), косвенно свидетельствует о точности и достоверности решений в рядовых точках зондирования в пределах исследуемой территории.
В качестве отладочного инструмента по развитию автоматической интерпретации предусмотрен вариант просмотра и, при необходимости, корректировки формируемой модели в режиме «рабочая модель» (рис. 5).
Как показала практика работ, по мере отладки и совершенствования алгоритма стабилизации необходимость в интерактивном режиме («рабочая модель») постепенно снижается, достигая к настоящему времени не более 10-15 %.
Рис.4. Фрагменты промежуточных этапов процесса стабилизации решений с контролем расхождения стабилизированных и контрольных решений Zср в параметрических точках, а также точности совпадения расчетных и наблюденных кривых зондирования 3
gj , Ш Пр.нмь-1154 ЛІ П».ет [ї
ftoo*Wb. 154. 4
Рис.5. Фрагмент интерактивной интерпретации в режиме “рабочая модель”
Одной из задач дальнейшего повышения информативности получаемых решений является полная автоматизация данного процесса с выбором оптимального решения на основе перебора всех возможных вариантов задания управляющих па-
Библиографический список
1. Колесников В.П. Основы интерпретации электрических зондирований. М: Научный мир, 2007. 248 с.
раметров и опорных параметрических решений, формализующая тем самым опыт и элементы интеллекта интерпретатора, полная реализация которых в рамках обычной используемой интерактивной интерпретации практически невозможна.
2. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука, 1974. 223 с.
About Increase of Unambiguity of Electric Sounding Interpretation
V.P. Kolesnikov
Perm State University, 614990, Perm, Bukirev st., 15 E-mail: poisk@psu.ru
A new method to improve the uniqueness of the automatic interpretation of the vertical electric sounding (VES) dates is proposed. It is based on the use of combination of regularizing elements, involved in the interpretive process in its basic stages. Extending of the inverse problem, the general scheme of algorithm constructing and methods for its implementation are given.
Key words: vertical electric sounding, computer interpretation, regularization , technology-
Рецензент - доктор геолого-минералогических наук Б.К. Матвеев
