Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

STOCHASTIC SIMULATION MODEL ESTIMATION OF INCURRED, BUT NOT REPORTED INSURANCE CLAIMS IN TERMS QUEUEING SYSTEMS

Butov A.A. 1 Galimov L.A. 1
1 International Relations Department Ulyanovsk State University
The article describes a model estimation of incurred but not reported claims in terms queueing systems. Estimates of the number insurance events and distribution function of the insurance claims for aircraft insurance and car insurance are built. They are presented by the example of statistics insurance company. Stochastic simulation modeling was carried out ib high-level programming language C#. Estimates of distribution function claims were built in Excel and Statistica. Graphics of IBNR simulation model are presented. The benefit of this model is no need to accumulate and segment statistics of insurance claims as required by the other actuarial methods.
simulation
queueing systems
IBNR
distribution function of insurance claims
1. Borovikov V. STATISTICA. Iskusstvo analiza dannyx na komp’yutere: dlya professionalov. 2-e izd. (+CD). SPb.: Piter, 2003. 688 р.: il.
2. Butov A.A., Ravodin K.O. Teoriya sluchajnyx processov: uchebno-metodicheskoe posobie. Ul’yanovsk: UlGU, 2009 g. 62 р.
3. Ivanov S.S., Golubev S.D., Chyornaya L.A., Sharafutdinova N.E. Teoriya i praktika riskovogo straxovaniya. M.: ROSNO: Ankil, 2007. 480 р.
4. Klark S.M., Khardi M.R. i dr.; uchebnye materialy. M.:2008, 438 р.
5. Mak T., Matematika riskovogo strakhovaniya / Per. s nem. M.: ZAO «Olimp-Biznes», 2005. 432 р.:il.
6. Prikaz Ministerstva Finansov Rossijskoj Federacii ot 11 iyunya 2002 g. no. 51n «Ob utverzhdenii pravil formirovaniya straxovyx rezervov po straxovaniyu inomu, chem straxovanie zhizni». URL: http://base.garant.ru/12127460/ (data obrashheniya: 13.06.2016).
7. Federalnyj zakon no. 4015-1 ot 27 noyabrya 1992 g. no. 4015-1 (red. ot 23.05.2016) «Ob organizacii strakhovogo dela v Rossijskoj Federacii» URL: http://base.garant.ru/10100758/1/ (data obrashheniya: 13.06.2016).
8. Friedland J. Estimation unpaid claims using basic techniques, March 2009; URL: https://www.casact.org/library/studynotes/Friedland_estimating.pdf (data obrashheniya: 13.06.2016).
9. StatSoft Russia. URL: http://www.statsoft.ru (data obrashheniya: 13.06.2016).

Одним из важных условий при заключении договора страхования является своевременное извещение клиента страховой компании о страховом случае. Однако промежуток времени от наступления страхового случая до сообщения о нем страховщику в период, обозначенный в договоре, может выходить за пределы отчетного периода. В связи с этим для исполнения обязательств по таким претензиям страховщик наряду с резервами незаработанной премии, заявленных, но неурегулированных убытков и стабилизационным резервом формирует резерв произошедших, но незаявленных страховых убытков (далее – РПНУ) [6].

Согласно [9] РПНУ является оценкой обязательств страховщика по осуществлению страховых выплат, возникших в связи со страховыми случаями, произошедшими в отчетном или предшествующих ему периодах, о факте наступления которых в установленном законом или договорном порядке не заявлено страховщику в отчетном или предшествующих ему периодах. Обязанность формирования и методика расчета РПНУ предусмотрены страховым законодательством [6, 9].

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

Для расчета РПНУ наиболее широко применяются методы, основанные на использовании треугольников развития убытков (оплаченных или состоявшихся) – методы цепной лестницы, Борнхюттера – Фергюсона, Берквиста – Шермана, а также методы, основанные на ожидаемой убыточности [3–5, 7, 8].

В настоящей статье рассматривается стохастическая имитационная модель оценки резерва произошедших, но не заявленных убытков, в виде одноканальной системы массового обслуживания (на примере авиационного страхования и страхования средств наземного транспорта), описанная в семимартингальных терминах. В качестве статистических данных количества и размера произошедших убытков использовались данные «Страховой компании НИК» (далее – Компания) за период 2010–2015 г.

Стохастическое моделирование системы массового обслуживания проводилось в программной среде высокого уровня Microsoft Visual Studio 2010 на языке C#, оценки функций распределения по эмпирической выборке страховых возмещений – в среде Microsoft Excel 2010 и Statistica 10.0.

Постановка задачи

Рассмотрим одноканальную систему массового обслуживания в терминах точечных процессов [2]. Пусть процесс (At)t≥0 – точечный считающий процесс числа произошедших убытков в момент времени t ≥ 0, (Dt)t≥0 – точечный считающий процесс числа заявленных, но неурегулированных убытков (убытки, о которых известно компании на момент времени t), (Qt)t≥0 – число произошедших, но незаявленных убытков. Взаимосвязь этих процессов представлена в виде системы массового обслуживания (далее – СМО), в которой роль заявок играет наступление страховых случаев, длину очереди – процесс (Qt)t≥0, а количество обслуженных заявок – процесс (Dt)t≥0 (рис. 1). Тогда балансовое уравнение запишется в виде

Qt = Q0 + At – Dt, (1)

где A0 = 0; B0 = 0; D0 = 0.

pic_60.wmf

Рис. 1. Общая схема работы СМО

Процессы At, Dt, Qt рассматриваются на стохастическом базисе B(Ω, F, ? = (Ft)t≥0, P), где At, Dt – независимые пуассоновские процессы с интенсивностью λ > 0 и δ > 0 соответственно. В соответствии с разложением Дуба-Мейера для субмартингалов [2] процессы At, Dt могут быть представлены как

butov01.wmf butov02.wmf (2)

где butov03.wmf и butov04.wmf – квадратично-интегрируемые мартингалы на B, а butov05.wmf, butov06.wmf – компенсаторы процессов
At, Dt вида

butov07.wmf butov08.wmf (3)

Стохастическая имитационная модель расчета РПНУ может быть построена следующим образом.

Пусть ηt – размер убытка, заявленного клиентом в момент времени t, есть случайная величина с функцией распределения F(η ≤ x), тогда It – величина РПНУ в момент времени t – может быть определена как

butov09.wmf (4)

Численный расчет количества произошедших, но незаявленных убытков

Для оценки параметров butov10.wmf butov11.wmf введем следующие параметры: butov12.wmf – дата i-го страхового случая, butov13.wmf – дата уведомления об i-м страховом событии страховщику, butov14.wmf – количество дней, исчисляемое от даты страхового случая до даты уведомления, butov15.wmf, i = 1...N, где N – количество страховых случаев за рассматриваемый период.

Так как время между наступлениями страховых случаев butov16.wmf имеет экспоненциальное распределение, то значение butov17.wmf может быть вычислено методом максимального правдоподобия [5], как

butov18.wmf (5)

Аналогично вычисляется значение интенсивности обслуживания как среднее количество дней между датами страхового случая и даты заявления об убытке:

butov19.wmf (6)

На рис. 2–3 приведены иллюстративные примеры моделирования процесса At – количества произошедших убытков за период (0, T) по статистическим данным Компании для страхования средств воздушного транспорта (КАСКО ВС, N1 = 40) с оценочными параметрами butov20.wmf butov21.wmf и для страхования средств наземного транспорта (КАСКО АМ, N1 = 316) с параметрами butov22.wmf butov23.wmf.

pic_61.wmf

Рис. 2. График процессов At butov24.wmf (КАСКО ВС), дискретность t = 1 месяц, T = 63 мес.

pic_62.wmf

Рис. 3. График процессов At butov25.wmf (КАСКО АМ), дискретность t = 1 месяц, T = 40 мес.

Математическая модель размера страхового убытка

Определим по эмпирической выборке η1, η2, ..., ηN страховых возмещений, функцию распределения случайной величины η в предположении, что для каждого вида страхования выборка страховых возмещений случайна и однородна.

Используя классический аппарат подгонки функции распределения по эмпирическим наблюдениям страховых возмещений, представленный в [4, 5, 7, 8], установим, что наиболее подходящей функцией распределения страхового возмещения для КАСКО ВС и КАСКО АМ (с точки зрения критерия Колмогорова-Смирнова) является логнормальное распределение с параметрами μ > 0, σ2 > 0 и плотностью вида

butov26.wmf x > 0. (7)

Параметры логнормального распределения по эмпирической выборке оценены методом максимального правдоподобия:

butov27.wmf butov28.wmf (8)

Основные параметры оценки функции распределения, критерий Колмогорова − Смирнова

Вид страхования

Объем
выборки

butov29.wmf

butov30.wmf

butov31.wmf

Критерий Колмогорова ? Смирнова

КАСКО АМ

316

5,6415?104

1,17?105

3,5016

0,93

0,068

КАСКО ВС

40

5,8756?105

1,078?106

0,6635

1,5842

0,084

pic_63.tif

Рис. 4. График EIt, (в USD) для КАСКО ВС, t = 1 месяц, T = 5000

pic_64.tif

Рис. 5. График EIt, (в тыс. рублей) для КАСКО АМ, t = 1 месяц, T = 1000

Пример подгонки эмпирической функции распределения к теоретической представлен на рис. 4.

В предположении стационарности параметров λ, δ, μ, σ определим среднюю величину EIt при T → ∞ (время моделирования) по формуле

butov32.wmf (9)

Заключение

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

Преимуществом использования настоящей имитационной модели расчета РПНУ является высокая степень адекватности реальным данным, что подтверждается графиками 2–6. Кроме того, применение настоящей модели не требует предварительного длительного процесса накопления данных об убытках и дискретизации страховых выплат (на квартал, полугодие, год) как этого требуют цепочно-лестничные методы.

Стохастическое моделирование системы массового обслуживания расчета РПНУ реализовано в среде Microsoft Visual Studio 2010 на языке C#. Подгонка функции распределения по эмпирическим данным страховых возмещений реализована в пакете Statistica 10.0.

Настоящая имитационная модель расчета РПНУ может быть усовершенствована в случае представления процессов At, Dt в виде мультивариантных процессов [2].