Переход к трехмерной псевдоматрице поля Ex

Измеренное поле Ex^f,гц= , где волновое число к= = +ί в высокомных разрезах используем в стационарном смысле. A(ω), не что иное, как измеренный амплитудный спектр.

По определению, за затухание отвечает мнимая часть волнового числа (фазовые углы), а действительная часть отвечает за длину волны и фазовую (групповую) скорость [5].

Из теории функций комплексных переменных известно.

Если функцию можно разложить в ряд Тейлора, то она является аналитической. Ряд Тейлора позволяет найти функцию f(ω) в регулярных точках, если ее значение известно в точке ω=ω₀. Все другие точки, в которых ряд Тейлора не существует, называются особенностями: полюса, существенно особые точки и точки ветвления – в них функция

разрывается.

На особые точки амплитудный спектр проверяется разложением в ряд Лорана.

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

Функция e^z –универсальная функция. В плоскости комплексной частоты по формуле Коши разложим спектр F(ω)=R(ω)+iIm(ω)=А(ω)e^iφ(ω). Амплитудный спектр измерен в трех узловых частотных точках на плоскости XY. Ищем фазовый спектр из преобразований Гильберта.

Из-за многозначности фазовых спектров, которые отвечают за форму амплитудных спектров, следует учесть, что при преобразовании Гильберта фазовые углы испытывают разрывы при переходе через значения, кратные π. Поэтому нужно дополнить значения углов в точках разрыва значениями ±2π (функция unwrap).

Основная формула спектрального анализа: приводит к обратным преобразованиям Фурье измеренного спектра или кепстральным преобразованиям на комплексной плоскости.

Подготовительные работы к реализации вычислительного алгоритма.

Подсчитана необходимая статистика по трем частотным срезам ρ_ω, статистика сохранена и визуализирована для понимания типа распределения ρ_ω. Геометрия и координаты профилей и точек измерений построены на плане и визуализированы в виде карты изоом и изоаномал, импеданса Ex/Hz.

Алгоритмические шаги.

1 шаг. Составили матрицу А(i,j,k), в координатах XYF. Задаются узловые точки или плоскости F= 4,88; 625 и 1256 Гц.

Назначаем координату X(0:20:i)-вектор строка и координату Y(0:20:j)-вектор –строка, в зависимости от площади съемки и шага между профилями и точками измерений на профиле. В данном случае в основании куба задается квадратная сеть.

Назначаем третью координату F=logspace(0;3,400019635;n=2512)-вектор-строка логарифмической линейки частот от декады 10⁰ до декады 10^3,400019635 Гц. Удвоенная частота f =2·1256 Гц=2512 Гц. Для визуализации внутри этого диапазона декад применяется удвоенная линейка частот [1: 100: 2512] и достаточен логарифмический шаг 10 или 100 Гц. Затем понадобится сменить дискреты частоты, связанные с теоремой Котельникова (F_Н≤1/2t₀ Найквиста) при обратном преобразовании Фурье. В узлах объемной сети находятся три значения ρ_ω

2 шаг. Визуализация куба матрицы A с управлением прозрачностью (OpenGL) и ее вращение для рассмотрения особенностей в поведении трех частотных плоскостей ρ_ω. Принятие предварительного решения о положении будущих линий разреза.

Для приближенного представления значений ρ_ω в межузловых частотных точках применим полиномы Лежандра. Физико-математическое обоснование следующее. Пространство гармонических сферических функций ρ_ω, удовлетворяют уравнению Лапласа в сферической системы координат (кроме нуля). Следовательно, одна из производящих функций (частный случай решения телеграфного уравнения) есть функция Лежандра. Полезным является ее свойство ортогональности. Попутно заметим: есть и другие производящие функции. Здесь важно, что бы физическое поле было описано аналитическим представлением без лишних осцилляций и с минимальными остатками (аналог ряда Фурье).

3 шаг. Проводим узловую интерполяцию в матрице - кубе A(i,j,k) изображения ρ_ω полиномами Лежандра (например, функцией legandre), тем самым переводим амплитудно-частотно-дискретное пространство трех узловых плоскостей в функциональное (предпочтение - нечетным функциям).

Задача 5-го шага перевести функциональное пространство ρ_ω в пространство Гильберта.

5 шаг. Через преобразование Гильберта выделили амплитудный и фазовый спектры. Амплитудный спектр в общем-то равен измеренному изображению ρ_ω. Фазовый спектр (в радианах) вращаем функцией unwrap сначала на ±π, тем самым корректируем разрывы, свойственные преобразованию. Сохранили изображения спектров в амплитудной и фазовой матрицах A(i,j,k) и Ф(i,j,k). Визуалировали куб фазовой матрицы Ф с возможностью корректировки разрывов (возможен ввод углов ±π/2,±π/4 и других). Находим регулярные разрывы. Сохраняем матрицу Ф (ω) с фиксированными фазовыми углами- разрывами. Это путь вслепую с ручными попытками увязать разрез с геологическими отбивками.

Задача 6-го шага - создать фазовый портрет электромагнитного поля конкретной геологической среды с контролируемым источником поля. Математически обоснованный и не ручной путь ниже.

6 шаг и неоднократный возврат к 5- му шагу.

На основании свойств ортогональности проводим разложение функционального пространства ρ_ω (матрица А) полиномами Тейлора через опорные частотные точки с шагом 10 или 100, 200 Гц (предстоит выбирать, а полиноминальных членов может быть хватит 5).

Для этого в качестве смотровой площадки вводим комплексную плоскость ω= real ω+iIm ω в единичном круге. Проецируем матрицу Ф на плоскость, проверяем на верхней полуплоскости наличие нулей, полюсов и точек ветвления- ищем регулярные фазовые разрывы.

Если их нет (нет границ), то F(ω) – регулярная функция и идем к пункту 7.

Признак нуля. Если есть одиночные нули- F(ω)=0, при Im ω>0, то представляем функцию - спектр как Ф(ω)=lnF(ω)=lnA(ω)+iφ(ω). Перевычисляем фазовый спектр φ(ω)=-1/π∫ dω, и F(ω)= F₀(ω) (может быть по кепстральной функции).

В случае регулярной цепочки нулей или их концентрированного облака - F(ω)= F₀(ω)∏ . F(ω) и F₀(ω) имеют одинаковые модули, но разные фазы: перевычисляем фазовый спектр φ(ω)=φ₀(ω₀)+φ(ω₀),где φ(ω₀)=-2Σtg^{-1 (теорема вычетов).}

Признак полюса. Если есть полюс, то ряд Лорана имеет конечное число членов.

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

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

7 шаг. После проверки нулей и полюсов, сохранения регулярных разрывов в фазовой матрице проводим частотный срез в нужном сечении координат XYF с проверенными на регулярность и не регулярность амплитудным и фазовым спектром. Математически это означает: из трехмерных матриц A и Ф координатную выборку двухмерных спектральных значений: амплитудный- a(ω), фазовый-φ(ω) и их соединение в двухмерную матрицу Mf по правилу: амплитудно-фазовый спектр Mf= sqrt(a(ω)²+φ(ω)²). Линейка частот остается прежней. Частотный срез может проводиться неоднократно, исходя из геологической целесообразности или назначения нескольких линий проектных геосечений.

Задача 8-го шага вычислить фазовые задержки в радианах для получения скоростного закона

8 шаг. В двухмерной матрице φ(ω) вычислили разницу регулярных фазовых задержек на измеренных частотах 1256-625 Гц, 625-4,88 Гц и 1256-4,88 (в радианах) и сохранили в двумерной матрице Дельта размерностью 3х3.

Задача 9 шага - создание 2-х мерного временного разреза. Он может быть фазовый, амплитудный или амплитудно-фазовый. Назначить три кнопки- какой разрез выбираете?

9 шаг. Применили обратное двумерное преобразование Фурье к матрице Мf(f,j) и перешли от частотной к временной оси. Время назначается в микросекундах вектором размерности sizе, функция ifft2(Mf, sizе). Теорема отчетов в частотной области. Смена дискретов частоты, связанная с теоремой Котельникова (∆ω≤2π/T_c, T_c=1/2πf_c, (∆f≤2πf_c). Наибольшая частота, на которой измерено поле- 1256 Гц. Удвоенная (зеркальная) частота = 2512 Гц. Из-за свойств преобразований Фурье, для времени это соответствует 1/5024 =0,199044 млсек или 199,044 мксек. Шаг по оси времени или дискреты времени - из ряда 10,50 мксек. Сохранили матрицу Mt. Визуализация временного разреза-матрицы Mt с масштабированием амплитуд ρ_ω .

10 шаг. По действительной части волнового числа и фазовым задержкам находим скорость в среде. V=30000/ - где V – вектор- строка для трех фазовых сдвигов, 30 000 м/мкsec– скорость электромагнитной волны в воздухе, ε- назначается относительная диэлектрическая проницаемость, ρ_ω- из амплитудного спектра A_ω по среднеквадратичному в вертикальной плоскости направлению.

Рассчитать среднеквадратичные значения ρ_ω в указанных мышью прямоугольных окнах. Например, по первой паре фазовых задержек на высоких частотах находим скорость в среде. И по второй паре фазовых задержек на низких частотах находим так же скорость в среде. При необходимости смотрим общий фазовый сдвиг на частотах 1256-4,88 Гц. Сравниваем их, если общий сдвиг примерно равен сумме двух сдвигов, значит среда является в этом диапазоне частот линейным фазовращателем.

В случае явного скоростного закона V=V(t₀) аппроксимируем скоростные точки простейшим полиномом и его применяем при переходе к глубинному разрезу.

Принимаем решение об использовании скорости в виде одиночного значения (среда однородна) или в виде скоростного закона (среда слоистая). Время t₀ – из 9-го шага. Визуализация графика V=V(t₀) для вставки в визуализацию шага 12.

Задача 11 шага - создание 2-х мерного глубинного разреза. Он может быть только амплитудно-фазовый.

11 шаг. Скорость (скоростной закон) умножаем на временную линейку в микросекундах - переводим в матрице Mt координатную ось времен в относительную координатную ось глубин в метрах (от нуля до скин-слоя, где амплитудные значения уменьшаются в exp раз (в 524 раза). Проверяем линейку глубин по формуле эффективного проникновения поля h= на трех частотах и трех среднеквадратичных эффективных сопротивлениях в разрезе. В случае их приемлемого совпадения и близких оценках, сохраняем матрицу Mh.

Задача 12 шага - облагораживание глубинного разреза – матрицы Mh. Назначение и рисовка на разрезе скважин, точек зондирований, текстов, пикетов, масштабов съемки. Вставка других уменьшенных изображений. Масштабирование амплитуд способом переменной плотности и цветом, усиление, фильтрации и т.д. с учетом бесконечного наращивания графических возможностей. Вывод на * .dxf

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

Уф!?? Спасибо за внимание. Коллеги! Кто вникнет, оценит идею и возьмется за написание программного продукта? Бесплатно! Как элемент своей диссертации.