доктор физико-математических наук;
doctor of physical and mathematical sciences;
Поставлена задача изучения разрешающей способности метода дирекционных магнитотеллурических зондирований (ДМТЗ). Показано, что резкие изменения импедансов электрического и магнитного типов в критической области параметров R~0, J~ωµσ приводят к существенному повышению разрешающей способности метода по сравнению с традиционной интерпретацией при помощи импеданса Тихонова–Каньяра. Рассматривается проблема устойчивости решения обратной задачи метода ДМТЗ с учетом погрешностей полевых измерений импедансов, накладывающих ограничения на разрешающую способность. Определен нижний предел этих ограничений для малых вариаций проводимости ∆σ/σ. Для изучения разрешающей способности метода ДМТЗ применительно к задаче магнитотеллурического мониторинга очага землетрясения выполнен математический эксперимент для трехслойной геоэлектрической модели. Показано, что при изменении проводимости очага ∆σ~±10 % вариации коэффициента отражения электрической моды находятся в пределах 20 %, что существенно превышает погрешность полевых измерений. Обсуждается возможность построения прогностических моделей в свете полученных результатов.
The problem is to examine the resolution of directional magnetotelluric soundings (DMTS). Abrupt changes of impedances of electrical and magnetic types in critical region of parameters R~0, J~ωµσ are shown to result in significantly higher resolution of the method as compared with the traditional interpretation using Tikhonov–Cagniard impedance. The stability of solution of DMTS method inverse problem is considered subject to field measurement errors limited the resolution. The minimum of the limits is determined for the conductivity ∆σ/σ small variations. For studying DMTS resolution as applied to MT monitoring of earthquake site the mathematical experiment for three-layer geoelectric model was carried out. When changing the earthquake site conductivity of ∆σ~±10 % variations of reflection coefficient of electric mode are in the range of 20 % that is significantly more than the field measurement error. The possibility for prediction modeling in the context of obtained results is discussed.
ВВЕДЕНИЕДирекционный анализ магнитотеллурических (МТ) наблюдений [Четаев, 1985], основанный на использовании модели неоднородной плоской волны (волны Четаева), расширил возможности традиционной интерпретации при помощи модели Тихонова–Каньяра [Ваньян, 1997]. При этом ключевым является вопрос о разрешающей способности методов МТ-зондирований (МТЗ), или степени точности, с которой можно определить параметры (σi, hi) i-го слоя геоэлектрического разреза, i=1, 2, …, N, где N — число слоев в разрезе. На практике речь идет о распознавании геоэлектрических структур, вычленении зон аномальной проводимости и т. п. Как известно [Бердичевский и др., 1986], при помощи метода МТЗ, обладающего большой глубинностью, оценивается проводимость сравнительно крупных структур земной коры. Так, в работе [Мороз и др., 2008] оценивается разрешающая способность МТ-зондирований для определения наличия проводящих зон под Камчаткой и магматических зон под вулканами. Показано, что выделение на кривых МТЗ проводящих слоев с удельным сопротивлением ρ~10 Ом•м оказывается затруднительным, если мощность слоя h≤1.5 км. При этом достаточно тонкие слои мощностью от десятков до первых сотен метров проводящих рудных залежей или залежей углеводородов, обладающих высоким сопротивлением (ρ~1011 Ом•м), этим методом неразличимы. Использование импеданса Тихонова–Каньяра Z0, например,для мониторинга проводимости очага землетрясения также остается проблематичным в силу его малой чувствительности [Жданов, 1986]. Давая оценку сложившейся в магнитотеллурике ситуации, авторы работы [Морозова и др., 1998] приходят к выводу, что «точность методов естественных полей оказывается недостаточной в свете требований, предъявляемых к изучению геоэлектрического строения земной коры и вариаций электропроводности». Разумеется, речь идет о методах интерпретации МТ-полей в рамках модели однородной плоской волны, исключающей Hz- и Еz-компоненты.Первые попытки интерпретации полного (6-компонентного) МТ-поля, выполненные для Украинского кристаллического щита способом аналитического продолжения коэффициентов отражения вниз [Четаев и др., 1984], свидетельствовали о его более высокой, по сравнению с методом МТЗ, разрешающей способности. При решении обратной задачи дирекционных МТЗ (ДМТЗ) методом невязки [Cавин, Израильский, 1985] было показано, что метод в принципе не имеет ограничений на разрешающую способность. Одновременно выяснилось, что малые погрешности при измерении импедансов, превышающие первые проценты, приводят к ухудшению устойчивости обратной задачи. Следовательно, при оценке чувствительности метода нельзя игнорировать степень точности полевых импедансных измерений.Таким образом, проблема разрешающей способности ДМТЗ требует дальнейшего изучения как в применении к различным задачам геоэлектрики, так и в отношении того или иного способа решения обратной задачи.