Ò à á ë è ö à 7.21
Участок кри- |
Период |
Вид модели |
Изменение добычи, |
Темп ввода сква- |
âîé |
106 ì3/ìåñ |
æèí, ñêâ/ìåñ |
||
I |
01.1976–12.1976 |
1734 + 1760t – 0,03t3 |
1760–1760 |
0,3 |
II |
12.1976–04.1979 |
2056 + 1406t – 6t2 |
1262–926 |
–0,03 |
III |
04.1979–12.1982 |
17317 + 782t |
926–782 |
0 |
IV |
12.1981–12.1982 |
–66 + 1564t – 6t2 |
1124–1051 |
–0,17 |
Во втором периоде (12.1976–04.1979 гг.) на месторождении Восточный Шатлык наблюдается уменьшение фонда добывающих скважин на 2 единицы, темп снижения составил 0,03 скв/мес. В табл. 7.22 эта величина приведена со знаком минус. Отмеченное нашло свое отражение и в уменьшении темпа роста добычи газа с 1760 106 ì3/мес в первом периоде до 1262 106 ì3/мес в начале второго. В то же время вид кривой отался таким же, как и в первом периоде. Следует отметить, что уменьшение фонда скважин привело к тому, что темп роста добычи газа снижается в течение всего второго периода и в конце составил 926 106 ì3/ìåñ.
Âтретьем периоде (04.1979–12.1982 гг.) отмечается стабилизация фонда скважин, поэтому модель, аппроксимирующая кривую добычи газа в этом периоде, есть прямая.
Âэто время стабилизация фонда скважин привела к некоторому росту темпа добычи газа, поэтому, несмотря на то, что в четвертом периоде (12.1981–
12.1982 гг.) произошло снижение фонда скважин (темп снижения равен 0,17 скв/мес), темп роста добычи газа в начале периода повысился по сравнению с третьим периодом до 1124 106 ì3/мес. Однако к концу четвертого периода он опять снизился – до 1051 106 ì3/ìåñ.
Анализ приведенных выше результатов позволил определить, что изменение фонда скважин приводит к изменению вида моделей или значений их параметров, описывающих кривую суммарной добычи газа. Увеличение количества скважин вызывает рост темпа добычи газа. Так, по месторождению Западный Шатлык, несмотря на снижение темпа ввода скважин, темп роста добычи постоянно повышается. В то же время уменьшение количества вводимых скважин вызывает также и уменьшение темпа роста добычи газа.
Для подтверждения полученных результатов была проведена серия математических экспериментов, суть которых заключалась в следующем.
Для гипотетической залежи задавался вид изменения отбора газа во времени Q (t). Далее по уравнениям материального баланса для периода нарастающей добычи газа в условиях газового режима без учета реальных свойств газа определялся закон изменения средневзвешенного по объему давления в залежи pñð. Зная депрессию ∆p, определялись значения забойного давления pç и по известным коэффициентам фильтрационных сопротивлений À è Â находились дебиты средних скважин. После этого оценивалось необходимое число
скважин. Для гипотетической залежи принимались следующие параметры: ðí = = 30 ÌÏà, ∆ð = 2 ÌÏà, αλ = 1 109 ì3, À = 0,25 ÌÏà2 ñóò/òûñ. ì3,  =
=0,0000044 (ÌÏà ñóò/òûñ. ì3)2.
Âкачестве базовых моделей были использованы зависимости (7.40–7.44).
Результаты расчетов показывают следующее:
при постоянстве фонда скважин изменение отбора газа во времени описывается линейной моделью;
при резком уменьшении числа добывающих скважин изменение отбора газа описывается гиперболической моделью;
671
|
|
|
Ò à á ë è ö à |
7.22 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Вариант рас- |
1 |
2 |
3 |
|
4 |
5 |
6 |
7 |
÷åòà |
|
|||||||
|
|
|
|
|
|
|
|
|
Вид модели |
(7.41) |
(7.44) |
(7.41)– |
|
(7.45) |
(7.42) |
(7.44), |
(7.45), |
(номер фор- |
|
|
(7.44) |
|
|
|
(7.45) |
(7.46), |
ìóëû) |
|
|
|
|
|
|
|
(7.42) |
|
|
|
|
|
|
|
|
|
равномерный ввод новых скважин соответствует параболической модели изменения отбора газа;
неравномерное увеличение или уменьшение фонда скважин соответствует экспоненциальным моделям изменения отбора газа с положительным или отрицательным показателем степени соответственно.
Результаты, полученные для модельных расчетов, полностью соответствуют результатам, полученным для месторождения Западный Шатлык.
Следует отметить, что знание моделей, описывающих кривую суммарной добычи газа, позволяет осуществлять прогноз. Так, кривая суммарной добычи газа по месторождению Западный Шатлык на последнем участке описывается уравнением
Q |
ã |
(t) = 106 |
(3 + 31å−0,025t ), |
(7.47) |
позволяющим осуществлять прогноз добычи газа.
Для проверки точности модели было проведено сравнение расчетных и фактических значений Qã по первым шести месяцам 1983 г. Погрешность прогноза не превышала при этом 1,5 %. В табл. 7.22 приведены виды исходных кривой и их сочетаний, а на рис. 7.25 – кривые изменения числа скважин в соответствии с вариантами расчета, указанными в табл. 7.22. Как видно из рис. 7.23–7.25, при параболическом росте суммарной добычи газа во времени отме- чается снижение числа скважин или незначительный рост, при экспоненциальном росте добычи газа наблюдается резкий рост числа вводимых скважин, а стабильность фонда добывающих скважин соответствует прямой зависимости роста добычи газа от времени.
Данные выводы находятся в полном соответствии с полученными выше результатами.
Рис. 7.25. Кривые, полученные по уравнению материального баланса:
1, 2 è 4 – номера вариантов из табл. 7.22
672
Следует отметить, что знание моделей, описывающих кривую суммарной добычи газа, позволяет производить прогноз. Так, для месторождения Восточ- ный Шатлык ошибка прогноза на первые 6 мес 1983 г. не превышала 0,6 %, а для месторождения Западный Шатлык, как уже отмечалось, она составляла за тот же период не более 1,5 %.
7.4. ИНФОРМАЦИОННОЕ ОБЕСПЕЧЕНИЕ АНАЛИЗА РАЗРАБОТКИ МЕСТОРОЖДЕНИЙ С БОЛЬШИМ ФОНДОМ СКВАЖИН
Эффективность разработки газовых и газоконденсатных месторождений с большим фондом скважин во многом определяется имеющимся информационным массивом промысловых данных, таких как дебиты скважин, пластовые, забойные, устьевые давления и т.д. Систематическое обследование фонда скважин, проведение необходимых замеров требуют как наличия достаточно большого набора измерительных приборов и значительного числа рабочих бригад, так и крупных материальных затрат. Естественно, что на крупных месторождениях с большим фондом добывающих скважин возникают трудности с замером показателей работы индивидуальных скважин. Так, добыча газа на Оренбургском месторождении осуществляется 650 скважинами. Число их на некоторых сборных пунктах (УКПГ) достигает 80–90. Каждый УКПГ оборудован одним контрольным сепаратором, при помощи которого замеряется дебит газа, конденсата и отбираются пробы продукции по каждой скважине. Исследования показывают, что для проведения комплекса замеров на одной скважине необходима ее эксплуатация через контрольный сепаратор в течение 5—7 сут. Таким образом, детальные исследования скважин и оперативный контроль, например, за дебитом газа индивидуальных скважин несовместимы. Поэтому необходимо применять методы, позволяющие по ограниченному количеству замеров восстанавливать информационный массив в необходимом объеме по всему фонду скважин. Использование аппарата порядковых статистик позволяет значительно сократить объем замеров показателей работы скважин. Расчеты проведены на примере Оренбургского газоконденсатного месторождения для замеров дебитов газа.
Рассмотрены две задачи: 1) определение дебитов скважин, подключенных к одному УКПГ; 2) определение дебитов по УКПГ в целом.
Определение дебитов индивидуальных скважин рассмотрено на примере УКПГ-2 (по второму эксплуатационному объекту). Предположим, что на некоторую дату известны дебиты всех скважин, подключенных к данному УКПГ. Требуется на следующую дату (например, через месяц) определить дебиты всех скважин, проводя замеры дебитов лишь в нескольких из них. Задача решается следующим образом. По известным дебитам скважин определяется закон распределения случайной величины — дебита скважины. Предполагается, что в будущем вид закона распределения не изменяется, могут измениться лишь параметры этого закона. Такое предположение, несмотря на кажущуюся справедливость при неизменных условиях эксплуатации, требует периодической проверки.
673
Ò à á ë è ö à 7.23
|
Номер |
Дебит газа |
|
Дебит газа |
|
Номер |
Дебит газа |
|
Дебит газа |
|
|||||
|
(ìàðò), |
Погреш- |
(апрель), |
Погреш- |
(январь – март), |
Погреш- |
(апрель), |
Погреш- |
|||||||
Ðàíã R |
скважи- |
скважи- |
|||||||||||||
òûñ. ì3/ñóò |
ность, % |
òûñ. ì3/ñóò |
ность, % |
òûñ. ì3/ñóò |
ность, % |
òûñ. ì3/ñóò |
ность, % |
||||||||
|
íû |
|
|
|
|
|
|
íû |
|
|
|
|
|
|
|
|
замер |
расчет |
|
замер |
расчет |
|
замер |
расчет |
|
замер |
расчет |
|
|||
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2006 |
118 |
– |
– |
99 |
– |
– |
2006 |
103 |
– |
– |
99 |
– |
– |
|
2 |
2008 |
118 |
– |
– |
99 |
– |
– |
2008 |
103 |
– |
– |
99 |
– |
– |
|
3 |
118 |
296 |
352 |
+19 |
247 |
242 |
–2 |
2010 |
299 |
357 |
+19 |
297 |
247 |
–17 |
|
4 |
803 |
395 |
383 |
–3 |
396 |
281 |
–29 |
803 |
399 |
387 |
–3 |
396 |
283 |
–29 |
|
5 |
809 |
395 |
413 |
+5 |
297 |
320 |
+8 |
809 |
399 |
417 |
+5 |
297 |
319 |
+7 |
|
6 |
2010 |
395 |
444 |
+12 |
297 |
360 |
+21 |
118 |
400 |
447 |
+12 |
247 |
354 |
+43 |
|
7 |
152 |
444 |
474 |
+7 |
396 |
400 |
+1 |
152 |
449 |
478 |
+6 |
396 |
390 |
–2 |
|
8 |
104 |
475 |
505 |
+6 |
475 |
439 |
–8 |
104 |
479 |
508 |
+6 |
475 |
426 |
–10 |
|
9 |
148 |
474 |
535 |
+13 |
445 |
479 |
+8 |
148 |
479 |
538 |
+12 |
474 |
462 |
–3 |
|
10 |
101-Ä |
592 |
566 |
–4 |
396 |
518 |
+31 |
101 |
599 |
568 |
–5 |
396 |
498 |
+26 |
|
11 |
2011 |
592 |
596 |
+1 |
495 |
558 |
+12 |
2011 |
599 |
598 |
0 |
495 |
534 |
+8 |
|
12 |
2007 |
641 |
626 |
–2 |
594 |
598 |
+1 |
2007 |
648 |
629 |
–3 |
594 |
570 |
–4 |
|
13 |
112 |
691 |
657 |
–5 |
693 |
637 |
–8 |
112 |
712 |
659 |
–7 |
693 |
606 |
–13 |
|
14 |
139 |
710 |
687 |
–3 |
673 |
677 |
+1 |
139 |
718 |
689 |
–4 |
673 |
641 |
+5 |
|
15 |
119 |
730 |
718 |
–2 |
732 |
716 |
–2 |
124 |
731 |
719 |
–2 |
594 |
676 |
+14 |
|
16 |
2012 |
740 |
748 |
+1 |
692 |
755 |
+9 |
119 |
738 |
750 |
+2 |
732 |
714 |
–2 |
|
17 |
124 |
829 |
779 |
–6 |
594 |
795 |
+34 |
2012 |
748 |
780 |
+4 |
692 |
749 |
+8 |
|
18 |
101 |
839 |
804 |
–4 |
841 |
835 |
–1 |
103 |
781 |
810 |
+4 |
742 |
785 |
+6 |
|
19 |
103 |
839 |
840 |
0 |
742 |
874 |
+18 |
101 |
848 |
840 |
–1 |
841 |
821 |
–2 |
|
20 |
106 |
849 |
870 |
+2 |
851 |
914 |
+7 |
106 |
858 |
871 |
+2 |
851 |
857 |
+1 |
|
21 |
108 |
947 |
901 |
–5 |
950 |
953 |
0 |
108 |
958 |
901 |
–6 |
950 |
893 |
–6 |
|
|
|
Qmin = 261 |
|
Qmin = 123 |
|
|
Qmin = 266 |
|
Qmin = 139 |
|
|||||
|
|
Qmax = 931 |
|
Qmax = 993 |
|
|
Qmax = 931 |
|
Qmax = 929 |
|
|||||
Суммарный де- |
12109 |
12516 |
+3 |
11004 |
11715 |
+6 |
|
12048 |
12568 |
+4 |
11004 |
11211 |
+2 |
||
áèò ïî ÓÊÏÃ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для определения дебитов скважин в следующие периоды времени необходимо уметь проводить ранжирование дебитов, т.е. расположение их в порядке возрастания или убывания, не зная их истинной величины. Процедуру ранжирования можно провести любым возможным способом, например, по устьевым давлениям, перепадам давлений, мнению работников промысла и т.п. Может быть использовано и предположение о неизменности рангов — дебиты скважин со временем изменяются, но их соотношение остается постоянным. В любом случае операция ранжирования может быть проведена с некоторой ошибкой, однако, как показывает анализ промысловых данных, получаемые результаты имеют приемлемую для практических целей точность. Далее по некоторым скважинам замеряли дебиты, и с учетом их рангов определяли неизвестные параметры распределения. После этого восстанавливаются значения дебитов всех скважин.
Дебиты скважин УКПГ-2 за март 1983 г. ранжировали по абсолютной величине (табл. 7.23). Для последующих интервалов времени дебиты скважин ранжировали по порядку для предшествующего значения времени. Для сравнения в графе 4 табл. 7.23 приведены расчетные значения дебитов в предположении, что известны ранги дебитов и значения дебитов скважин, имеющие ранги 5 и 15 (Q5 = 395 òûñ. ì3/ñóò, Q15 = 730 òûñ. ì3/сут). Функция распределения имеет в данном случае вид
|
0, |
|
Qmin >Q; |
|
|
|
Q − Qmin |
|
|
|
|
Ô(Q)= |
, |
Qmin <Q <Qmax ; |
(7.48) |
||
|
|||||
Qmax − Qmin |
Q >Qmax . |
|
|||
|
1 |
|
|
||
|
|
|
Значения Qmin è Qmax приведены в табл. 7.23. Замеры дебитов скважин за апрель, ранжированные по предыдущим (мартовским) замерам, приведены в графе 6, а результаты расчетов и погрешности — в графах 7 и 8.
Как видно, даже при ошибке в ранжировании точность восстановления дебита в большинстве случаев удовлетворительна. При этом во всех случаях точность определения суммарного дебита по УКПГ достаточно высока, хотя добыча газа в апреле снижена по отношению к предыдущим месяцам на 9 —
12%.
Âграфах 9—15 в табл. 7.23 приведены результаты аналогичных расчетов, когда за основу определения вида распределения и ранжирования дебитов взяты дебиты скважин за январь — март. Как видно, точность восстановления де-
битов возросла.
Из результатов проверки гипотезы о равномерном распределении дебитов газа скважин по десяти УКПГ (табл. 7.24) следует, что только для скважин УКПГ-2 и УКПГ-7 распределение дебитов газа носит неравномерный характер. В этом случае для нахождения дебитов можно использовать определенную последовательность действий.
Ранжируют каким-либо образом дебиты газа скважин. Проводят замеры дебита, например, каждой пятой скважины. На рис. 7.26 такие скважины обозначены звездочками. По ним строят зависимость дебита от ранга. Дебиты остальных скважин экстраполируют по полученной зависимости.
На рис. 7.26 показана зависимость дебита газа от ранга для скважин УКПГ-2. Видно, что она хорошо описывается двумя прямыми, т.е. распределение дебитов для скважин УКПГ-2 включает два равномерных закона распределения, так как в случае равномерного закона зависимость дебит — ранг пред-
675
|
|
Ò à á ë è ö à |
7.24 |
|
|
||
|
|
|
|
|
|
|
|
ÓÊÏÃ |
Число |
Распределение |
|
ÓÊÏÃ |
Число |
Распределение |
|
скважин |
выборки |
|
скважин |
выборки |
|||
|
|
|
|
||||
1 |
23 |
Равномерное |
|
|
8 |
18 |
Равномерное |
2 |
51 |
Неравномерное |
|
|
9 |
38 |
« |
2 |
21 |
Равномерное |
|
|
12 |
23 |
« |
3 |
33 |
« |
|
|
14 |
27 |
« |
6 |
27 |
« |
|
|
15 |
22 |
« |
7 |
38 |
Неравномерное |
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 7.26. Зависимость дебита газа от ранга
äëÿ ÓÊÏÃ-2
Рис. 7.27. Гистограмма распределения деби-
676
|
|
|
Ò à á ë è ö à 7.25 |
òîâ äëÿ ÓÊÏÃ-2 |
|
||
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
Ðàíã R |
Номер УКПГ |
Дебит газа (ян- |
Дебит газа (апрель), млн. м3 |
|
Погрешность, |
||
(номер объекта) |
âàðü – ìàðò), |
замер |
расчет |
|
% |
||
|
|
ìëí. ì |
3 |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
10(2) |
3,673 |
|
2,967 |
|
|
|
2 |
9(1) |
7,283 |
|
5,067 |
4,3 |
|
–15 |
3 |
15(3) |
9,665 |
|
10,923 |
8,7 |
|
–20 |
4 |
1(1) |
13,765 |
12,126 |
13,3 |
|
10 |
|
5 |
8(3) |
14,707 |
10,663 |
18,1 |
|
70 |
|
6 |
7(1) |
16,133 |
13,440 |
23,1 |
|
77 |
|
7 |
9(3) |
27,024 |
23,212 |
28,4 |
|
22 |
|
8 |
3(1) |
30,077 |
28,528 |
33,8 |
|
18 |
|
9 |
14(1) |
30,903 |
27,903 |
39,6 |
|
44 |
|
10 |
14(3) |
39,518 |
50,906 |
45,6 |
|
–10 |
|
11 |
15(1) |
50,230 |
46,325 |
52,0 |
|
13 |
|
12 |
2(1) |
72,076 |
59,655 |
58,7 |
|
–2 |
|
13 |
7(2) |
85,732 |
68,398 |
65,8 |
|
–3 |
|
14 |
10(1) |
113,089 |
80,533 |
73,4 |
|
–1 |
|
15 |
6(3) |
115,362 |
101,382 |
81,5 |
|
–20 |
|
16 |
8(2) |
138,642 |
101,161 |
90,3 |
|
–11 |
|
17 |
1(3) |
148,600 |
123,359 |
99,7 |
|
–19 |
|
18 |
15(2) |
153,232 |
147,152 |
110,0 |
|
25 |
|
19 |
6(2) |
154,625 |
130,478 |
121,0 |
|
–7 |
|
20 |
3(3) |
155,484 |
135,618 |
133,5 |
|
–1 |
|
21 |
12(2) |
171,951 |
149,312 |
147,4 |
|
–1 |
|
22 |
1(2) |
244,198 |
147,815 |
163,0 |
|
10 |
|
23 |
12(3) |
296,721 |
228,088 |
181,2 |
|
–21 |
|
24 |
9(2) |
327,919 |
260,72 |
202,7 |
|
–22 |
|
25 |
14(2) |
333,027 |
303,487 |
229,0 |
|
–24 |
|
26 |
2(2) |
430,852 |
296,772 |
262,8 |
|
–11 |
|
27 |
7(3) |
437,283 |
325,462 |
310,5 |
|
–5 |
|
28 |
3(2) |
458,283 |
429,954 |
392,0 |
|
–9 |
|
29 |
2(3) |
458,74 |
402,773 |
|
|
|
|
Суммарный дебит по всем УКПГ |
4160 |
|
3724 |
3396 |
|
–9 |
|
|
|
|
|
|
|
|
|
ставляется одной прямой. Дебиты по промежуточным рангам, обозначенные на рис. 7.26 точками, достаточно хорошо описываются полученной зависимостью.
При решении задачи об определении дебитов УКПГ учитывалось, что к УКПГ подключены скважины, эксплуатирующие разные объекты. В табл. 7.25 приведены проранжированные данные о дебитах УКПГ по разным объектам за январь — март 1983 г. Гистограмма распределения дебитов, приведенная на рис. 7.27, соответствует экспоненциальному закону распределения Ф(Q) = 1 — å–λQ. В целях повышения надежности результатов величина Q определялась по замерам дебитов за апрель, соответствующих рангам R 4, 8, 20 и 27. В итоге получено λ = 0,0898 (среднее значение). С целью оценки устойчивости λ также определялось для дебитов, соответствующих рангам R 5, 13, 19 и 26. В этом случае λ = 0,0950, т.е. значение λ устойчиво для данной выборки дебитов. В табл. 7.25 приведены расчетные дебиты объектов и их погрешности. Расчетный суммарный дебит отличается от фактического на 9 %.
677
7.5.ПРОГНОЗИРОВАНИЕ ДОБЫЧИ ГАЗА
ИКОНДЕНСАТА
Современный этап разработки газовых и газоконденсатных месторождений требует совершенствования методов оценки конечной газо- и конденсатоотдачи. Особую актуальность эти проблемы приобретают при долгосрочном планировании добычи газа и конденсата, находящейся в определенной зависимости от разведанных запасов. А точный прогноз извлекаемых запасов обеспечит эффективность разработки газоконденсатных месторождений.
Прогноз можно использовать как на этапе проектирования схемы разработки новых объектов (выбор рациональной схемы разработки), так и на этапе эксплуатации. В частности, в начальный период падающей добычи можно прогнозировать конечный коэффициент газоконденсатоотдачи. Решение данной задачи рассмотрим на примере месторождения Южный Мубарек для горизонтов XII и XIII, которые находятся на поздней стадии разработки, и, следовательно, когда имеется достаточное количество информации для прогнозирования конечных газо- и конденсатоотдачи. Здесь приводится несколько методов прогнозирования, поскольку они дают различные результаты в зависимости от количества и характера исходной информации.
Наиболее распространенный и простой путь выявления тенденции развития — сглаживание (выравнивание) динамического ряда. Данный метод основан на применении взвешенных скользящих средних. При этом каждому уровню в пределах интервала сглаживания приписывается вес, зависящий от расстояния, измеряемого от данного уровня до середины интервала сглаживания. Пусть для m последовательных уровней ряда со сдвигом во времени на один шаг подбирают полиномы вида:
ói = ai + bi i + ci i2 + … |
(7.49) |
Остановимся для примера на параболе |
|
ói = ai + bi i + ci i. |
(7.50) |
Тогда средневзвешенные скользящие будут определяться по формуле
ói = 0,1(2ói+2 + ói+1 – ói–1 – 2ói–2), |
(7.51) |
т.е. веса симметричны относительно центрального уровня.
Расчет по этой формуле приводит к тому, что сглаженная кривая в значи- тельной мере сохраняет различные изгибы кривой тренда.
Иногда при анализе динамики возникает необходимость определения не самого тренда, а средней скорости изменения среднего прироста. Такая задача возникает при выборе функции для аналитического выравнивания динамиче- ского ряда.
Формула для вычисления среднего прироста при m = 5 такова:
|
|
(t) = (−2yi −2 − yi −1 + yi+1 + 2y i+2 ) /10. |
(7.52) |
U |
Остановимся теперь на вопросе выравнивания динамических рядов с оценкой кривой роста. Для выравнивания наиболее часто применяются много- члены, экспоненты и другие кривые. Многочлены имеют следующий вид:
678
I степени yt = a0 + a1t; |
|
II степени yt = a0 + a1t + a2t2; |
(7.53) |
III степени yt = a0 + a1t + a2t2 + a3t3, |
|
n-й степени yt = a0 + a1t + a2t2 + a3t3 + … |
+ antn, |
ãäå ai (i = 1, 2, … ) — параметры многочленов; t — независимая переменная (время).
Многочлен I степени характеризуется постоянством прироста ординат и поэтому применяется для описания равномерно развивающихся во времени процессов.
Многочлен II степени (парабола II степени) описывает движение с равномерным изменением приростов, причем приросты больше нуля для одной ветки и меньше нуля для другой. Легко показать, что приросты можно охарактеризовать уравнением прямой.
U (1) |
= y |
t |
− y |
t −1 |
= a + 2a |
t − a |
3 |
= |
(a |
− a |
3 |
) + 2a |
t. |
(7.54) |
|
t |
|
|
1 |
2 |
|
|
1 |
|
2 |
|
|
||||
Приросты II порядка — постоянные числа: |
|
|
|
|
|
|
|
||||||||
|
|
|
|
Ut(2) = Ut |
−Ut −1 |
= 2a2 . |
|
|
|
|
(7.55) |
Для параболы III степени знак прироста может изменяться 1 или 2 раза:
U (1) |
= y |
t |
− y |
t −1 |
= 3a |
3 |
t2 + (2a |
2 |
− 3a |
3 |
)t − |
(a |
1 |
+ a |
2 |
+ a |
3 |
); |
(7.56) |
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
Ut(2) = Ut |
−Ut −1 = 2a3 |
+ 6a3 t. |
|
|
|
|
|
|
(7.57) |
Из всех экспонент самый простой вид имеет кривая y = abt, характеризующаяся постоянным темпом прироста.
Логарифмируя последнее выражение, получаем
lg yt = lg a + t lg b,
далее
α = lg a è β = lg b |
(7.58) |
è
lg yt = α + βt.
Расчет по этой формуле приводит к тому, что сглаженная кривая в значи- тельной мере сохраняет различные изгибы кривой тренда.
При выравнивании динамических рядов широко применяются законы роста функции. Для этого наиболее часто используются многочлены, различного рода экспоненты и логистические кривые.
Многочлены имеют вид формул (7.53).
Наиболее усложненный вариант экспоненциальной кривой — логарифми- ческая парабола
y |
t |
= abt ct2 . |
(7.59) |
Проведя несложные преобразования, получим
679
lg yt = lg a + t lg b + t2 lg c. |
(7.60) |
Далее, если найти темп прироста в виде отношения производной yt к ординате, то производная равна
y′ |
= abt ct 2 ln b + 2abt ct2 |
ln c. |
(7.61) |
t |
|
|
|
Тогда темп прироста будет |
|
|
|
|
τ0 = lg b + 2t lg c. |
|
(7.62) |
Отсюда видно, что темп прироста линейно зависит от времени.
Для процессов, характеризующихся насыщением, вид кривой — модифицированная экспонента
yt = K + abt, |
(7.63) |
имеющая горизонтальную асимптоту y = K ïðè t → ∞.
Особенность модифицированной экспоненты — постоянство отношений приростов:
Ut |
/Ut |
= Ut |
3 |
/Ut |
2 |
= Ut /Ut −1 = ... = b = const. |
|
2 |
1 |
|
|
|
|
||
Кроме того, логарифмы приростов линейно зависят от времени: |
|
||||||
Ut = yt − yt −1 = K + abt − K − abt −1 = abt −1 (b −1), |
(7.64) |
||||||
откуда |
|
|
|
|
|
|
|
lg Ut |
= lg a + t lg(b – 1) + (t – 1) lg b. |
(7.65) |
Очень большое распространение на практике получили S-образные кривые, одним из видов которых является кривая, описываемая зависимостью yt = = Kabt, имеющая характерную особенность — отношение последовательных приростов ординат в логарифмах постоянно:
|
lg yt +1 − lg yt |
= |
lg a(bt +1 − bt ) |
= b = const. |
(7.66) |
|
lg yt − lg yt−1 |
|
lg a(bt − bt −1) |
|
|
Прологарифмировав yt = Kabt, получим |
|
||||
|
lg yt = lg K + bt lg a, |
(7.67) |
ò.å. lg yt представляет собой модифицированную экспоненту, а логарифм прироста есть линейная функция от времени.
Второй тип S-образной кривой — логистическая кривая
1 |
= K + abt , |
(7.68) |
yt |
|
которую иногда называют кривой Перла – Рида. Часто ее записывают в виде
yt = |
|
K |
, |
(7.69) |
|
+ bef (t) |
|||
1 |
|
|
ãäå f(t) — некоторая функция от t, f(t) = – at. Тогда
680