Сглаживание исходных данных
Номер пункта |
Исходные данные |
Первое сглаживание |
Второе сглаживание |
Третье сглаживание |
|||||
x |
f(x) |
f1 |
δ1 |
f2 |
2 |
f3 |
3 |
||
1 |
8,6 |
9,5 |
10,2 |
0,7 |
11,0 |
1,5 |
11,4 |
1,9 |
|
2 |
9,8 |
11,6 |
12,5 |
0,9 |
12,3 |
0,7 |
12,8 |
1,2 |
|
3 |
11,5 |
16,3 |
14,2 |
–2,1 |
15,1 |
–1,2 |
15,0 |
–1,3 |
|
4 |
13,8 |
14,8 |
18,7 |
3,9 |
17,6 |
2,8 |
17,7 |
2,9 |
|
5 |
20,0 |
25,0 |
19,9 |
–5,1 |
20,5 |
–4,5 |
20,1 |
–4,9 |
|
6 |
24,0 |
19,8 |
22,9 |
3,1 |
22,3 |
2,5 |
22,5 |
2,7 |
|
7 |
25,5 |
24,0 |
24,1 |
0,1 |
24,6 |
0,6 |
24,5 |
0,5 |
|
8 |
29,2 |
28,5 |
26,8 |
–1,7 |
26,7 |
–1,8 |
26,8 |
–1,7 |
|
9 |
35,8 |
28,0 |
29,3 |
1,3 |
29,0 |
1,0 |
28,7 |
0,7 |
|
10 |
43,8 |
31,4 |
30,9 |
–0,5 |
30,4 |
–1,0 |
30,1 |
–1,3 |
|
11 |
50,5 |
33,4 |
30,9 |
–2,5 |
31,0 |
–2,4 |
30,5 |
–2,9 |
|
12 |
55,2 |
27,9 |
31,1 |
3,2 |
30,0 |
2,1 |
29,9 |
2,0 |
|
13 |
61,5 |
32,1 |
28,0 |
–4,1 |
28,8 |
–3,3 |
28,8 |
–3,3 |
|
14 |
62,5 |
24,0 |
27,4 |
3,4 |
27,5 |
3,5 |
27,8 |
3,8 |
|
15 |
63,9 |
26,0 |
27,0 |
1,0 |
27,0 |
1,0 |
26,9 |
0,9 |
|
16 |
69,2 |
31,0 |
26,7 |
–4,3 |
26,1 |
–4,9 |
26,1 |
–4,9 |
|
17 |
73,1 |
23,2 |
24,6 |
1,4 |
24,3 |
1,1 |
24,1 |
0,9 |
|
18 |
79,6 |
19,5 |
21,6 |
2,1 |
21,8 |
2,3 |
21,9 |
2,4 |
|
19 |
82,2 |
22,1 |
19,1 |
–3,0 |
19,5 |
–2,6 |
19,9 |
–2,2 |
|
20 |
87,7 |
15,8 |
17,9 |
2,1 |
18,3 |
2,5 |
18,7 |
2,9 |
|
|
|
|
|
|
|
|
|
||
Дисперсия |
|
45,9 |
- |
7,28 |
- |
6,10 |
- |
6,71 |
|
Для сглаживания выберем окно размером в три измерения. Возьмем первые три значения f(x) и вычислим среднее значение:
получим сглаженное значение во второй строке таблицы. Переместим окно на одно измерение и вычислим следующее значение:
Когда первое сглаживание закончено, рассчитаны отклонения сглаженных значений от исходных значений 1 и найдена дисперсия отклонений (7,28), можно выполнить второе сглаживание и снова найти отклонения сглаженных значений от исходных значений f(x), а также дисперсию отклонений (6,10). Сглаживание повторяется и оканчивается после достижения минимума дисперсии отклонений. В данном примере минимум дисперсии достигнут после второго сглаживания (табл.5.5). На рис.5.11 совмещены исходные и сглаженные значения.
5.4.2. Гармонический анализ
Как упоминалось выше, в составе закономерной изменчивости часто присутствует периодическая составляющая. Ее можно выделить и вычесть из исходных данных с помощью гармонического или периодограммного анализа. Гармонический анализ позволяет ряд исходных данных представить как сумму синусоид. Но здесь, в отличие от классического гармонического анализа, необходимо выделить наиболее существенные синусоиды, которые вносят основной вклад в изменчивость пространственной переменной.
Периодические явления широко распространены в природе, они находят свое выражение в ритмичном напластовании горных пород, в развитии систем упорядоченных трещин, в формировании рудных столбов. Чтобы выявить и охарактеризовать периодичность, необходимо применить специальные математические приемы.
Предлагаемая методика характеризуется тем, что в ряду наблюдений последовательно находят синусоиды в порядке убывания их значимости и вычитают из исходных данных. Можно извлечь несколько (иногда одну) важнейших синусоид и этим ограничиться. Но можно извлекать синусоиды до тех пор, пока дисперсия отклонений суммы синусоид от исходных данных с учетом использованных степеней свободы не будет минимальной.
Каждая синусоида имеет три характеристики: амплитуду А, длину волны L (или обратную ей величину – частоту) и начальную фазу , поэтому каждая синусоида использует три степени свободы. Кроме того, часто находят дисперсию синусоиды, которая пропорциональна квадрату амплитуды. Ряд значений пространственной переменной f(х) представляют в виде суммы N гармоник (синусоид) k(x):
f(х) = 
,
(5.13)
где
– свободный член, равный среднеарифметическому
значению пространственной переменной.
В качестве первой длины волны L чаще всего принимается длина ряда наблюдений. Каждая гармоника k(x) выражается формулой
(5.14)
или в другой форме записи:
.
(5.15)
Здесь ak, bk – коэффициенты; k – номер гармоники (синусоиды); L/k – длина волны; Ak – амплитуда; k – начальная фаза гармоники; k = 1, 2, …, N.
Амплитуда и начальная фаза связаны с коэффициентами ak, bk соотношениями
.
(5.16)
Количество гармоник N может быть большим, но на практике принимается конечным. Количество гармоник не должно превышать n/2 (где n – количество измерений), но наиболее правильно ограничивать их количество по минимальной дисперсии отклонений с учетом степеней свободы, как в подразделе 3.1.7 или 4.1.3. Каждая синусоида использует три степени свободы, поэтому k синусоид поглощают 3k степеней свободы. Впрочем, для геологических целей обычно выбирают одну или две важнейшие синусоиды, пренебрегая остальными. За счет ограничения количества синусоид их сумма в формуле (5.13) приближенно равна пространственной переменной f(х).
Объединение формул (5.13) и (5.14) дает ряд Фурье (гармонический ряд):
f(х) = 
.
(5.17)
Коэффициенты ak, bk находят по формулам
);
(5.18)
.
(5.19)
Когда коэффициенты ak и bk определены, вычисляют амплитуды гармоник по формуле (5.16). Совокупность всех гармоник образует спектр амплитуд. Чем больше амплитуда, тем бόльшую роль играет соответствующая гармоника. Для количественной оценки роли гармоник используют их дисперсии Dk, совокупность которых составляет спектр дисперсий. Сумма дисперсий всех гармоник равна дисперсии пространственной переменной, что позволяет оценить роль каждой гармоники в абсолютных или относительных единицах. Вычисление дисперсии любой гармоники лучше всего осуществлять путем сравнения остаточных дисперсий до вычитания и после вычитания гармоники из значений пространственной переменной. Уменьшение остаточной дисперсии характеризует дисперсию Dk, поглощенную данной гармоникой.
Пример 5.9. Имеются результаты 25 измерений пространственной переменной f(x) – содержаний цинка в скважине по сети с шагом 1 м (табл.5.7). Очевидно, что длина отрезка (начальная длина волны) L = 25 – 1 = 24 м. Необходимо разложить пространственную переменную в ряд Фурье и оценить роль различных гармоник.
Задачу решим путем последовательного вычисления и вычитания гармоник из значений пространственной переменной (кстати, наиболее рациональный метод при неравномерной сети наблюдений). Вначале из значений пространственной переменной вычтем свободный член ряда (5.17), что даст отклонения 0(х) = f(x) – a0/2. Далее из отклонений 0(х) вычтем первую гармонику 1(х), что даст остаток 1(х). Из него вычтем вторую гармонику и этот процесс продолжим до заданного числа гармоник N = 12. Начало вычислений, включая первые две гармоники, показано в табл.5.7. Коэффициенты ak и bk каждой очередной гармоники найдем через отклонения (х) по преобразованным формулам (5.18) и (5.19):
;
(5.20)
,
(5.21)
где
– вспомогательный аргумент, вводимый
для сокращения вычислений.
Таблица 5.7