Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Бескровный Н.Т. Экономика и оптимизация надежности и ремонта горношахтного оборудования

.pdf
Скачиваний:
7
Добавлен:
24.10.2023
Размер:
11.13 Mб
Скачать

Т а б л и ц а 19

 

Число шахт со среднесуточной производственной

Себестоимость і т угля, руб.

 

 

мощностью, т

 

 

 

956 -

2130-

 

свыше

и т о г о

 

ДО 955

3305-

 

2130

3304

4480

4480

шахт

 

 

Свыше 22,4 .......................

6

3

2

 

 

9

18,6—2 2 ,4 ...........................

10

33

_

_

45

14,7—1 8 ,5 ...........................

18

46

11

1

_

76

10,7—1 4 ,6 ...........................

3

37

20

8

3

71

До 10,7 ...........................

4

21

15

4

44

И т о г о шахт

37

123

54

24

7

245

имеет определенное сходство с гиперболической кривой, которая часто используется (см., например, [11], [14]) при анализе себе- -стоимостн.

Множественную регрессионную зависимость между себестоимо­ стью Т т угля (У) и влияющими на нее факторами (Xj) искали в виде неполного квадратного много­

члена

(-1 і- 1

Рис. 67.' Эмпирическая линия регрессии себестоимости 1 т угля Y по средне­ суточной добыче шахты Х 7

Рассмотрено 12 факторов, которые могли оказывать влияние на себестоимость добычи 1 т угля. Средне­ суточная добыча угля в чис­ ловом значении представлена

обратной величиной

В ре­

зультате многошагового регрессионного анализа по критерию ti = = min (формула 101) полностью были исключены факторы: средне­ динамическая мощность пласта, число действующих очистных забоев, объем проведения подготовительных выработок за счет эксплуата­ ционных затрат, число и протяженность поддерживаемых выработок,

ичлены управления, содержащие квадраты некоторых факторов.

Вокончательном виде множественное регрессионное уравнение

себестоимости 1 т угля имеет вид

Y = 32,31 + - 7,60Х2 - 0,0342Х3- 0,5786Х4 + 7,83JCB+

+ І ,4 5 Х в + 0 , 3 9 2 Х 7 - ^ | р - + 7 ,9 • 1 0 - БХ 23 + 5 ,9 • 1 0 -3Х | . (1 1 0 )

150

Исключенные факторы и члены уравнения содержащие Х\, Х\, Х |, Х%, не оказывают существенного влияния на себестоимость угля, так как не несут дополнительной (по отношению к оставшимся) информации об изменчивости изучаемого показателя. Таким образом, следует считать, что в условиях комбината Донецкуголь значитель­ ное влияние на себестоимость угля по шахте оказывают факторы X L Х 7 (см. табл. 18).

Уравнение (110) получено при высоком индексе множественной корреляции R = 0,88, что соответствует очень тесной связи между функцией У и аргументами X f . Значение индекса множественной

Рис. 68. Теоретическое нормальное распределение (1) и фак­ тическое распределение плотности (2) по уравнению себестои­ мости угля

детерминации 7?? = 0,78 указывает на то, что 78%

изменчивости

зависимой

переменной

У

обусловлено

вариацией

независимых

переменных. Дисперсионное

отношение

S|p

= 4,34 показывает,

^ —

насколько

лучше выбор

в

 

^ост

уравнения (110) по-

предположениях

сравнению с оценкой функции У по ее средней У [70].

Распределение плотности ошибки Y t Y t хорошо согласуется с теоретическим нормальным законом распределения плотности, вероятностей (рис. 68).

Распределение ошибок в процентах к Y t приведено в табл. 20,. из которой видно, что более чем в 70% случаев ошибка не превышает

12%.

Средняя ошибка составляет 9,2%, что приемлемо для использо­ вания уравнения на практике.

Анализ уравнения (110) показывает, что на всем интервала изменения переменных (см. табл. 18) увеличение добычи угля,.

•15і

Т а б л и ц а 20

 

 

 

Число случаев с отклоненном,

%

Себестоимость 1 т угля, руб.

0—і 1,9

12,0-26,0

26,1-43,0

ИТОГО

 

 

случаев

Свыше 18,4 .......................

 

37

12

6

55

12,5—18,4 ...........................

 

94

19

5

118

До 12,5 ...............................

 

44

28

72

И т о г о случаев

. . . .

175

59

И

245

Уо ................................................

 

71,5

24,0

4,5

100

объемной массы угля и снижение отношения количества выдаваемой породы к добытому углю, числа шахт в административной единице и числа разрабатываемых пластов, а также увеличение средней длины лавы п подвигапия линии очистного забоя на большей части интервала изменения этих переменных приводят к уменьшению себестоимости угля, что соответствует установившимся представле­

ниям о

этих факторов на себестоимость угля.

Факторыв л и я Хо,н и и

Х 5, Ха, X-j оказывают преимущественно линейное

влияние на себестоимость угля, что было установлено в результате проведения многошагового регрессивного анализа исключением членов уравнения, содержащих квадраты этих факторов. При этом увеличение средней объемной массы угля на 0,1 т/м3 и уменьшение отношения количества выдаваемой породы к добытому углю на 0,1 приводят к снижению себестоимости 1 т угля на 0,76 руб. и 0,783 руб, соответственно. При возможности уменьшения числа шахт в административной единице и числа разрабатываемых пластов на 1 наиболее вероятное снижение себестоимости 1 т угля составит соответственно 1,45 и 0,392 руб.

Нелинейное влияние на себестоимость угля оказывают факторы Х г, Х3, Х±. Исследования показали, что на всем интервале изменения •среднесуточной производственной мощности шахты увеличение ее приводит к замедленному уменьшению себестоимости 1 т угля, а увеличение средней длины лавы и среднемесячного подвигания линии очистного забоя целесообразно доводить до 216 д 49 м соот­ ветственно так как с дальнейшим увеличением их наблюдается рост •себестоимости угля.

При подсчете экономической эффективности технических и орга­ низационных мероприятий приходится определять долю условно постоянных и условно переменных производственных затрат в се­ бестоимости продукции. В литературе 14] приводятся примерные границы участия этих затрат. Однако эти границы имеют большой диапазон и поэтому могут служить только для приблизительных расчетов. Для более точного определения доли условно постоянных и условно переменных затрат проводится детальное исследование -затрат по процессам и категориям трудящихся.

452

/

Уравнение (110) позволяет без кропотливой и трудоемкой работы по выборке затрат из неупорядоченной документации шахт опре­ делить соотношение между условно переменной и условно постоянной частями производственных затрат в себестоимости 1 т угля. В этом заключается основное отличие разработанной экономико-матема­ тической модели себестоимости угля от предложенных ранее [42], [48].

Производственные затраты 3 на добычу Х г т угля можно пред­ ставить в виде функций

 

3 = СпХ г + В,

(ill)

где СпХ г — условно

переменные

(пропорциональные

объему до­

бычи угля) затраты,

руб.; В — условно постоянные, не зависящие-

от объема добычи угля затраты,

руб.

 

На основании формулы (111) себестоимость 1 т угля определится

по формуле

 

У = С а + ^ - .

(112)

Учитывая, что

 

где Z — количество породы, выдаваемой за сутки, т, на основаниит формулы (110) коэффициенты Сп и В уравнения (112) определятся из выражений

Сп = 32,31 - 7,60Х2- 0,0342Х3 - 0 , 5786Х, + 1,45Х6+

 

+ 0,392Х7 + 7,9

• Ю-5Х! + 5,9 • 10-3Х24,

(ИЗ)

В = 4759 +

7.83Z - 48^700 .

(114)

Анализ уравнения (114) показывает, что по абсолютной величине условно постоянная часть затрат В не остается постоянной, а повы­ шается при увеличении количества выдаваемой породы и нагрузки, на шахту. Однако при увеличении нагрузки на шахту часть этих затрат в себестоимости 1 т угля уменьшается более интенсивно,, за счет чего достигается экономия на условно постоянной части,

затрат.

Q

Распределение отношения -у-, рассчитанного на основании

уравнений (110), (ИЗ) по фактическим данным шахт комбината Донецкуголь, приведено в табл. 21. В среднем отношение пропор­ циональной части затрат С„ к общей себестоимости У составляет- 0,62. Остальная доля (0,38) приходится на условно постоянные затраты.

Из табл. 21 видно, что имеется тенденция к увеличению отноше-

Q

ния -у- по шахтам с малой фактической ^себестоимостью угля. Эта: тенденция наглядно показана на рис. 69.

155

Т а б л и ц а 21

 

 

 

с

 

 

Число шахт с отношением ——

 

Себестоимость 1 т угля, руб.

 

 

 

и т о г о

0,33 -

0,50-

0,60-

0,71-

0,49

0,59

0,70

0,87

шахт

Свыше 18,4 ...............................

И

27

14

3

55

12,5—1 8 ,4 ...................................

И

46

48

13

118

До 1 2 , 5 .......................................

5

31

36

72

И т о г о шахт . . . . . . .

22

78

93

52

245

% .........................................................

9,0

31,8

38,0

21,2

100

Поскольку упорядоченного учета условно постоянных и условно переменных производственных затрат на шахтах не ведется, исполь­

 

зование

приведенной

мето­

 

дики

значительно

облегчает

 

расчет

экономической

эф­

 

фективности проводимых ме­

 

роприятий.

Особенно эффек­

 

тивно

может быть

использо­

 

вана

методика разделения

Рис. 69. Теоретическая линия регрессии

затрат

по

себестоимости

при оценке

массовых

меро­

себестоимости 1 т утля У по среднесуточ­

ной добыче шахты X t

приятий, проводимых в мас­

 

штабе комбинатов,

например

при оптимальном распределении оборудования, при

установле­

нии эффективных ремонтных нормативов

и др.

 

 

 

Использование авторегрессионных

моделей

 

 

для прогнозирования производительности

оборудования

 

очистных

забоев

 

 

 

 

 

Прогнозировать основные технико-экономические параметры ис­ пользования машин можно с помощью метода анализа экономических рядов динамики, основанного на использовании авторегрессионной модели 120].

Рассматривая ряд динамики как случайный процесс, можно сказать, что он состоит из части, детерминированной во времени (общей тенденции развития процесса и периодических колебаний) и случайных отклонений от тенденции. Последние связаны с непра­ вильным выбором аппроксимирующей кривой и с небольшим объемом выборки.

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

154

Ряд динамики представляется в виде последовательных значений

У х, У 2, . . Y t изучаемого

показателя

через равные

промежутки

времени (например, год), где

У, соответствует значению показателя

в момент времени, начиная

с

которого

производятся

измерения,

У а — в последующий момент,

и т. д.

 

 

Известно, что на изменение большинства технико-экономических показателей существенным образом может влиять изменение различ­ ных факторов. Учитывая, что изменение этих факторов случайно, можно предположить, что величины У х, У 2, .’. ., Y t связаны между собой таким линейным соотношением, при котором значение пара­

метра Y t

в каждый момент времени зависит только от предыдущих

значений

этого же

параметра.

 

 

Yt —

аіУt~i 4- aiYt-%+ . ,

(115)

где zt — значение некоторой случайной величины.

В уравнении (115) коэффициенты at отражают общую тенденцию развития и периодические колебания изучаемого показателя, а ве­ личина е{ — случайное отклонение от тенденции.

Для построения уравнения (115) значение величины У в момент времени t аппроксимируется с помощью линейной комбинации к

функций ѵг (t), и2 (£), . . .,

vk (£),

которые

принимают значения

v1{t) = Y i _ i

при

г = 1, 2,

3 . . .

Эти функции рассматриваются как независимые переменные.

 

 

 

Независимые переменные

Т а б л и ц а

2 2

Время

Зависимая

 

 

 

ui (О

Ѵ„ (/)

V, (О

.

” п- 1

измерения

переменная

 

 

 

 

 

 

1

Ух

 

 

_

_ _

_

 

У і

 

2

Y 2

Ух

 

 

3

Y 3

 

 

 

І — І

УІ-2

У и

Y t - t

~ ъ

■ Уг

t — 1

Y t - 1

Yt-2

У * - з

Y t - i

t

Yt

Y t - 1

Уі-2

Y t - з

^ 2

У і

 

Схема представления числовых данных, подлежащих обработке, приведена в табл. 22, из которой видно, что чем больше выделено' переменных для аппроксимации, тем меньше значений Yt-t> для которых определены все функции ѵг (t), i>2 (£), . . ., vk (t). На­ пример, если для аппроксимации будет выделена одна переменная,, то количество членов зависимой переменной, для которых опреде­ лены члены независимой переменной vt (£), окажется равным t — У Если будет выделено две переменных, то количество членов зависимой

155

переменной, для которых определены все члены независимых пере­

менных Ѵі

(t),

окажется равным

t — 2. И, наконец,

если

будет

выделено

п

— 1 переменных, то

окажется, что только

для

одного

члена Y t

зависимой переменной будут определены все члены неза­

висимых переменных vL(t).

 

 

 

Известно,

что для однозначного определения коэффициентов

щ уравнения

(115) необходимо определить полностью

не менее к

членов зависимой и независимой переменных. Таким образом, порядок уравнения модели к не может превышать величины

где п — число моментов времени, в которые производились замеры показателя У.

Допустим, что для прогнозирования выбрана модель к-го порядка.

В этом

случае

находим коэффициенты а х , а 2 ,

.

. . ,

а к так,

чтобы

2

=

% { Y t - a 1 Y t - 1 - a i Y l - i - . . . - a

k Y

t -

k ) a = mia.

(116)

t=/и-i

1

 

 

 

 

Взяв частные производные по a t и приравняв их к нулю, находим, что минимум суммы квадратов отклонений линейной функции (фор­ мула 116) от имеющихся наблюдений находится решением оиотемы к нормальных уравнений с неизвестными а£:

2

Y

t Y t . ! = «! І Y

U

 

+

^

2

Y

i , 2 Y t . 1

+ . . .

i=h+1

t=k+l

 

 

 

 

 

t=hbl

 

 

 

 

2

 

. . . - j - a È 2 1Y t - k Y t - i

2

 

1 У ? - 2

“Г • • •

1

=

2

 

 

 

 

 

+

fl2

 

t=h+

 

t=k-tl

Y(.jYf_2

t=h+

 

 

 

 

 

 

"

 

 

 

 

 

 

 

 

 

(117)

 

 

■ ■ • + « *

2

 

Y

t - k Y

t - 2

 

 

 

 

 

 

 

 

t=A+1

 

 

 

 

 

 

 

 

І Y

t Y

t . k = a 1 2

Y

^ Y

t

 

- b

 

+ a , 2

 

 

Y t ^ Y

t . k + . . .

t=h+1

 

t =ft+x

 

 

 

 

 

 

i =fe+l

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

■ ■ ■ J r

a k

 

2

 

Y t - к -

 

 

 

 

(

 

 

 

 

t

- к

+

і

 

 

 

 

 

Как уже

указывалось,

количество

независимых переменных

в авторегрессионной модели

не

 

может

превышать

величины

Для определения необходимого числа переменных, при котором достигается наилучшая аппроксимация рассматриваемого процесса, используют критерий, предложенный Манном и Вальдом. Он удобен тем, что использует непосредственно вычисляемые, а не оцениваемые величины. Критерий Манна и Вальда основан на сравнении остатков,

456

полученных в результате аппроксимации наблюдаемых значений ряда динамики моделью (115) порядка к и моделью порядка д, где q > к по формуле

(пg)ln-|^-,

(118)

Zih

 

где — сумма квадратов остатков, определенных по авторег­ рессионной модели (см. 115) порядка д; 2^ — сумма квадратов остатков, определенных по авторегрессионной модели (см. формулу 115) порядка к для тех же наблюдений, что и модели порядка д.

Величина (см. формулу 118) имеет распределение %с q к сте­ пенями свободы. Следовательно, при определенном уровне значи­ мости можно проверить, отвергается ли гипотеза в предположениях модели с меньшим порядком к моделью с большим порядком д. Увеличение порядка модели следует вести до тех пор, пока окажется, что модель с большим порядком не отвергает гипотезу в предположе­ ниях модели с меньшим порядком. Это дает возможность утверждать, что исследуемый процесс представляет собой авторегрессию порядка д.

Таким образом устанавливаются общая тенденция развития

ипериодические колебания во времени исследуемого показателя. Оценка той части процесса, которая характеризуется случайными

отклонениями, производится из предположения, что отклонения Е( представляют собой независимые случайные величины путем определения доверительного интервала

Et = ±ха,

(119)

где X — квантиль нормального распределения, соответствующая за­ данной вероятности; а — стандартная ошибка.

Стандартная ошибка определяется по формуле

° = ]/~Цг>

(12°)

где N = п к.

Условие независимости остатков может быть проверено двумя способами, которые требуют соответственно вычисления коэффици­ ента автокорреляции и отношения среднего квадрата последователь­ ных разностей к дисперсии [29]. Второй их этих показателей, назы­ ваемый коэффициентом Неймана, считается общепринятым при анализе экономических рядов динамики. Коэффициент Неймана

вычисляется по формуле

 

2 (е*+і—ЕЭ2

 

К = -----^ ------.

(121)

Вычисленное значение К сравнивается с табличным для назна­ ченного уровня значимости и таким образом устанавливается су­ щественность автокорреляции остатков.

157

Существенная автокорреляция остатков указывает на то, что интервал времени между двумя соседними замерами изучаемого показателя был выбран заниженным. Это приводит к преуменьшению стандартной ошибки,, а следовательно, и доверительного интервала. В этом случае необходимо увеличить интервал между замерами или соответствующим образом скорректировать величину стандарт­ ной ошибки в большую сторону. Методика выбора необходимого интервала или корректировка стандартного' уклонения здесь не рассматриваются.

Иллюстрация вычисления коэффициента Неймана и критерия Манна и Вальда, их применение и истолкование результатов будут приведены при рассмотрении конкретного примера.

Определенные на основании решения системы (117) коэффициенты а,- могут быть также использованы в практике перспективного прогнозирования на т периодов вперед.

Модель прогнозирования в общем виде записывается следующим образом

Y t+m = a i Y t+m -l Т а 2 ^і+ т -2 + • • • + a k Y t+m -k

( 1 2 2 )

Из этой формулы видно, что прогноз показателя в период t + т основывается на прогнозах этого же показателя в предыдущие периоды. В случае если используется модель первого порядка, формула (122) может быть существенно упрощена

(123)

В связи с тем что при перспективном прогнозе вероятность ошибки увеличивается, величину доверительного интервала предлагается определять с учетом удаленности прогноза от последнего момента времени, в котором проведено измерение показателя по формуле

(124)

Формула (124) учитывает возрастание неопределенности прогноза в связи с увеличением периода прогноза.

По этой методике по комбинату Донецкуголь был составлен прогноз среднемесячной производительности комбайнов на год и дан доверительный интервал этого прогноза для некоторого зна­ чения вероятности. Для этого на основании отчетных данных со­

ставлена табл. 23.

Вначале рассмотрим построение модели авторегрессии первого

порядка.

в этом случае выражается одним уравне­

Система уравнений (117)

нием с одним неизвестным

 

 

2 ^ 7

, . ! = ^ 2 Y U -

(125)

158

д

 

д

Et

%

О

и

1

2

11958

21959

31960

41961

51962

61963

71964

81965

91966

101967

Т а б л и ц а 23

Производи^ тельиость комбайна, т/месяц

Р

Р

Р

Р

Р

V, (0

Р

Р

!

 

 

 

 

 

 

 

 

 

3

4

5

6

7

8

9

10

И

12

4567

 

 

 

 

 

'--

5159

4567

5194

5159

4567

5159

5194

5159

4567

5417

5159

5194

5159

4567

6000

5417

5159

5194

5159

4567

6523

6000

5417

5159

5194

5159

4567

7020

6523

6000

5417

5159

5194

5159

4567

7156

7020

6523

6000

5417

5159'

5194

5159

4567

7590

7156

7020

6523

6000

5417

5159

5194

5159

4567

Вычисление сумм удобно производить по схеме, приведенной в табл. 24, вторая и третья колонки которой составляются на осно­ вании данных табл. 23.

Т а б л и ц а 24

а

 

 

гН

 

 

е і

 

 

941

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

Н

 

1

1

 

СЧ

f-ч

 

1

 

 

1

1

Z

 

 

 

 

 

 

 

 

 

1

4567

 

23 561 153

20 857489

 

 

 

_

_

2

5159 4567

 

3

5194 5159

26 795 846

26 615 281

4567

23 720 998

20 857489

23 561 158

4

5159

5194

26 795 846

26 977 636

5159

26 615 281

26 615 281

26 795 846

5

5417

5159

27 946 303

26 615 281

5194

28 135 898

26 977 636

26 795 846

6

6000

5417

32502 000

29 343 889

5159

30 954000

26 615 281

27 946 303

7

6523

6000

39 138 000

36 000 000

5417

35 335 091

29 343889

32 502 000

8

7020

6523

45 791 460

42 549 529

6000

42120 000

36 000000

39 138 000

9

7156

7020

50 235 120

49 280 400

6523

46 678 588

42 549 529

45 791 460

10

7590 7156

54 314040

51 208 336

7020

53 281 800

49 280 400

50235120

л- 1 0

327 079 768 309 447 841

 

 

 

Л-10

 

303 518 615 288 590 352

— 286 841 636 258 239 505 272 765 728

2

2-3

 

 

 

 

 

 

 

 

 

Л - 0

 

 

272 765 728 258 239 505

_

_

 

_

 

2

___

___

 

___

і - а

 

 

 

 

 

 

 

 

 

Подставляя данные первой итоговой строки табл.

24 в формулу

(125), находим ах = 1,0570 и составляем уравнение

авторегрессии

первого

порядка

У, = •1,05707^.

 

7

 

(126)

 

 

 

 

 

 

159