Введение

Волатильность является одним из важнейших параметров в оценке опционов, управлении рисками и построении торговых стратегий. Классическая модель Блэка-Шоулза-Мертона, предполагающая постоянную волатильность, не способна отразить динамику рынка, где наблюдаются эффекты «улыбки волатильности» и кластеризации. Для более точного описания рыночных процессов разработаны модели стохастической волатильности, среди которых наиболее известными являются модель Хестона и модель SABR. Эти подходы учитывают случайный характер изменений волатильности и позволяют более адекватно оценивать деривативы.

"Улыбка волатильности" и её влияние на модели

С момента финансового краха рынка в 1987 г. было замечено, что волатильность опционов с более низкими ценами базового актива выше, чем у опционов с более высокими ценами. Это свидетельствует о стохастической природе волатильности, которая варьируется во времени и в зависимости от уровня цен базового актива. Экономисты разработали модели стохастической волатильности для учёта этого эффекта.

Улыбка волатильности (volatility smile) — это графическое представление зависимости подразумеваемой волатильности от страйк-цены для опционов с одинаковым временем до истечения. Обычно подразумеваемая волатильность имеет тенденцию к повышению как для глубоких в деньгах (deep in-the-money), так и для глубоких вне денег (deep out-of-the-money) опционов, что создает форму, напоминающую улыбку.

График улыбки волантильности
График улыбки волантильности

Небольшое отступление

В своём телеграм канале делюсь ещё больше полезным контентом по сфере деривативов и децентрализованных финансов: https://t.me/kirrya_achieves

Математически, улыбка волатильности определяется как функция подразумеваемой волатильности \sigma от страйк-цены K для фиксированного времени до истечения T:

 \sigma = \sigma(K) \quad \text{при фиксированном } T

Поверхность волатильности (volatility surface) — это трехмерное представление волатильности, отображающее изменение подразумеваемой волатильности в зависимости от страйк-цены и времени до истечения опциона. Поверхность волатильности помогает лучше понять и моделировать рыночные ожидания относительно будущих колебаний цены базового актива.
Математически, поверхность волатильности \sigma(K, T) определяется как функция страйк-цены K и времени до истечения T:

 \sigma = \sigma(K, T)

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

• Ось ? — страйк-цена ?,
• Ось ? — время до истечения ? ,
• Ось ? — подразумеваемая волатильность ?.

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

Пример поверхности волатильности
Пример поверхности волатильности

"Улыбка волатильности" отражает тот факт, что подразумеваемая волатильность опционов с разными страйк-ценами различается, хотя срок их истечения и базовый актив совпадают. Графически это проявляется в виде параболической формы с загнутыми вверх краями.

Это явление противоречит модели Блэка-Шоулза, которая предполагает, что волатильность должна быть одинаковой для всех страйков. Реальность показала, что после краха 1987 г. трейдеры осознали возможность экстремальных событий и ценовых искажений, что потребовало пересмотра моделей волатильности.

Альтернативное объяснение улыбки волатильности: маркет-мейкер, продавая опционы около at-the-money, хеджируется более out-of-the-money опционами, другие игроки, продавшие эти OTM опционы, хеджируются еще более далекими, и так далее. Таким образом, ступенчато распределяются по поверхности волатильности открытые позиции.

Проблема изменяющейся волатильности

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

  • Возврат к среднему (mean reversion): После резких всплесков волатильность имеет тенденцию возвращаться к долгосрочному среднему значению.

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

Игнорирование этих особенностей может приводить к неверной оценке опционов, ошибкам в хеджировании и арбитражных стратегиях.

"Возврат к среднему" волатильности, симуляция из модели Хестона
"Возврат к среднему" волатильности, симуляция из модели Хестона

Модель Хестона

В 1993 году Стивен Хестон в своей работе "A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options" \cite{heston1993closed} предложил свою модель для оценки стоимости деривативов. Основная идея состоит в моделировании стохастической волатильности с возвратом к среднему. Примерно такое поведение волатильности и наблюдается на реальных рынках: долгое время волатильность колебнется около среднего значения, в редких случаях "взрывается" и через некоторое время возвращается к среднему значению. Модель описывается следующей системой стохастических дифференциальных уравнений:

\begin{aligned} dS_t &= \mu S_t \, dt + \sqrt{v_t} \, S_t \, dW_t^S, \\ dv_t &= \kappa (\theta - v_t) \, dt + \sigma \sqrt{v_t} \, dW_t^v, \end{aligned}

где:

  • S_t — цена базового актива в момент времени t;

  • v_t — мгновенная вариация (квадрат волатильности);

  • μ — средняя доходность актива;

  • κ — скорость возврата волатильности к среднему значению θ;

  • σ — коэффициент волатильности самой волатильности («vol-of-vol»);

Корреляция между броуновскими процессами dW_t^S и dW_t^v описывается как:\begin{equation} E[dW_t^S \cdot dW_t^v] = \rho , dt \end{equation}

Какие преимущества у модели Хестона?

В отличие от модели Блэка-Шоулза, где волатильность фиксирована, модель Хестона учитывает её случайные колебания во времени, что делает её более реалистичной. Волатильность не остаётся высокой или низкой бесконечно, а стремится к некоторому долгосрочному среднему значению, что соответствует реальному поведению рынков.

Важным преимуществом модели Хестона является тот факт, что для европейских call, put опционов в рамках этой модели существует аналитическое решение. Этот факт позволяет проводить быструю, вычислительно не затратную калибровку поверхности волатильности европейских опционов.

Есть ли недостатки?

Одной из главных проблем модели является высокая вычислительная сложность. Несмотря на то, что для некоторых вариантов (например, для европейских опционов) можно получить аналитические решения, в ряде случаев приходится прибегать к численным методам, таким как метод Фурье-интегрирования или Монте-Карло симуляция. Это не только замедляет процесс оценки деривативов, но и усложняет процедуру калибровки параметров модели. В реальных условиях калибровка часто проводится с использованием оптимизационных алгоритмов, таких как Левенберга–Марквардта, что помогает заметно ускорить процесс калибровки

Для того чтобы вариация v_t оставалась строго положительной, должен выполняться так называемый Feller condition: [ 2\kappa\theta \geq \sigma^2. ]
На практике же часто получается, что этот баланс нарушается, особенно при калибровке модели на данных с высоким уровнем волатильности. Нарушение условия Феллера может приводить к возникновению незначительных отрицательных значений волатильности в численных симуляциях, что создает дополнительные сложности при интерпретации результатов и требует специальных корректирующих методов.

Модель SABR

Модель SABR (Stochastic Alpha, Beta, Rho) впервые введена Патриком Хаганом, Андреем Лесниевским, Дариной Вудвард в 2002 году в статье "Managing Smile Risk". Динамика модели описывается следующим образом:

Эта модель описывается следующими стохастическими дифференциальными уравнениями:

\begin{aligned}   &dF_t = \alpha_t F_t^{\beta} dW_t\\    &d\alpha_t = v \alpha_t dZ_t\\\end{aligned}

F_0 = F, \alpha_0 = \alpha

Ограничения на параметры: 0 \leq \beta \leq 1, \alpha \geq 0, v \geq 0, -1 \leq \rho \leq 1
где:
F_t — цена форварда в момент времени t,
\alpha_t — волатильность актива в момент времени t,
\nu — параметр "волатильности волатильности",
\beta — параметр, определяющий степень зависимости волатильности от цены актива.

Аналогично модели Хестона, корреляция между броуновскими процессами W_t и Z_t определяется:

\begin{equation} E[dW_t \cdot dZ_t] = \rho  dt \end{equation}

Калибровка модели SABR

Для удобства, введем переобозначения:

\begin{aligned} &F_m = \sqrt{F K}\\ &\zeta = \dfrac{v}{\alpha}F_m^{1-\beta}\log \dfrac{F}{K}\\ &X(\zeta, \rho) = \log \dfrac{\sqrt{1 - 2 \zeta \rho + \zeta^2} + \zeta - \rho}{1-\rho}\\ &q_1 = \dfrac{(\beta-1)^2 \alpha^2 F_m^{2\beta - 2}}{24}\\ &q_2 = \dfrac{\rho \beta \alpha v F_m^{\beta - 1}}{4}\\ &q_3 = \dfrac{2-3\rho^2}{24}v^2\\ &S = 1 + T(q_1 + q_2 + q_3)\\ &D = F_m^{1-\beta}\left[ 1 + \dfrac{(\beta-1)^2}{24} \log^2 \dfrac{F}{K} + \dfrac{(\beta-1)^4}{1920}\log^4 \dfrac{F}{K}\right]\\ \end{aligned}

Тогда формула для подразумеваемой волатильности примет вид:

\sigma = \dfrac{\alpha S}{D} \times \dfrac{\zeta}{X(\zeta, \rho)}

Наша задача - найти частные производные к \sigma по параметрам \alpha, \beta, \rho, v

\log \sigma = \log \alpha + \log S - \log D + \log \zeta - \log X(\zeta, \rho)\begin{aligned} &\dfrac{\partial S}{\partial \alpha} = T\left(2 \dfrac{q_1}{\alpha} + \dfrac{q_2}{\alpha} \right) = \dfrac{T}{\alpha} \left(2q_1+q_2\right)\\ &\dfrac{\partial \zeta}{\partial \alpha} = -\dfrac{\zeta}{\alpha}\\ &\dfrac{\partial X}{\partial \alpha} = \dfrac{\partial X}{\partial \zeta} \dfrac{\partial \zeta}{\partial\alpha}\\ \end{aligned}\begin{aligned} &\dfrac{\partial S}{\partial \rho} = \dfrac{T \beta \alpha v F_m^{\beta - 1}}{4} - \dfrac{\rho v^2}{4}\\ &\dfrac{\partial X}{\partial \rho} = \dfrac{1}{1-\rho} - \dfrac{1}{\sqrt{\ldots}}\dfrac{\sqrt{\ldots} + \zeta}{\sqrt{\ldots} + \zeta - \rho}\\ \end{aligned}\begin{aligned} &\dfrac{\partial S}{\partial v} = T\left( \dfrac{q_2}{v} + 2 \dfrac{q_3}{v}\right) = \dfrac{T}{v}\left(q_2 + 2q_3\right)\\ &\dfrac{\partial \zeta}{\partial v} = \dfrac{\zeta}{v}\\ &\dfrac{\partial X}{\partial v} = \dfrac{\partial X}{\partial \zeta} \dfrac{\partial \zeta}{\partial v}\\ \end{aligned}\begin{aligned} &\dfrac{\partial \zeta}{\partial \beta} = - \log F_m \zeta\\ &\dfrac{\partial X}{\partial \beta} = \dfrac{\partial X}{\partial \zeta}\dfrac{\partial \zeta}{\partial \beta}\\ &\dfrac{\partial S}{\partial \beta} = T \left( \dfrac{2 q_1}{\beta - 1} + 2 \log F_m q_1 + \dfrac{q_2}{\beta} + \log F_m q_2 \right) = T \left( 2q_1 \left[ \dfrac{1}{\beta - 1} + \log F_m \right] + q_2 \left[ \dfrac{1}{\beta} + \log F_m \right] \right)  \\ &\dfrac{\partial D}{\partial \beta} = -\log F_m D + F_m^{1-\beta} \left[ \dfrac{\beta - 1}{12} \log^2 \dfrac{F}{K} + \dfrac{(\beta - 1)^3}{480} \log^4 \dfrac{F}{K} \right]\\ &\\ &\\ \end{aligned}

Таким образом:

\begin{aligned} &\dfrac{\partial \log \sigma}{\partial \alpha} = \dfrac{1}{\alpha} + \dfrac{1}{S}\dfrac{\partial S}{\partial \alpha} + \dfrac{1}{\zeta} \dfrac{\partial \zeta }{\partial \alpha} - \dfrac{1}{X} \dfrac{\partial X}{\partial \alpha}\\ &\dfrac{\partial \log \sigma}{\partial \rho} = \dfrac{1}{S} \dfrac{\partial S}{\partial\rho} - \dfrac{1}{X} \dfrac{\partial X}{\partial \rho}\\ &\dfrac{\partial \log \sigma}{\partial v} = \dfrac{1}{S}\dfrac{\partial S}{\partial v} + \dfrac{1}{\zeta}\dfrac{\partial \zeta}{\partial v} - \dfrac{1}{X}\dfrac{\partial X}{\partial v}\\ &\dfrac{\partial \log \sigma}{\partial \beta} = \dfrac{1}{S} \dfrac{\partial S}{\partial \beta} - \dfrac{1}{D} \dfrac{\partial D}{\partial \beta} + \dfrac{1}{\zeta}\dfrac{\partial \zeta}{\partial \beta} - \dfrac{1}{X}\dfrac{\partial X}{\partial \beta}\\ \end{aligned}

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

Результаты калибровки модели SABR применительно к опционам на Ethereum на бирже Deribit
Результаты калибровки модели SABR применительно к опционам на Ethereum на бирже Deribit

Исследования показывают, что модель SABR на рынке крипто-опционов часто демонстрирует более высокую точность , особенно для опционов с короткими сроками экспирации, в то время как модель Хестона может иметь большую ошибку при неидеальной калибровке.

Преимущества модели SABR:

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

К тому же, модель SABR способна достаточно точно воспроизводить наблюдаемую на рынке улыбку волатильности. За счёт параметра\beta модель адаптируется к различным типам динамики базового актива: при (\beta = 1) она приближается к лог-нормальному распределению, а при (\beta = 0) — к нормальному. Такая гибкость позволяет использовать модель для широкого спектра финансовых инструментов, от процентных ставок до валютных деривативов.

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

Практическое применение

Стохастические модели волатильности находят применение в следующих областях:

  • Маркет-мейкинг: Точные модели помогают устанавливать оптимальные котировки, снижая риск убытков из-за резких рыночных движений.

  • Риск-менеджмент: Модели используются для прогнозирования экстремальных событий и корректировки портфелей в условиях высокой волатильности.

  • Алготрейдинг: Быстрая оценка краткосрочных изменений волатильности позволяет принимать оперативные решения в высокочастотных стратегиях.

Применение данных моделей способствует более точной оценке деривативов и снижению операционных рисков.

Заключение

Стохастические модели волатильности, такие как модели Хестона и SABR, значительно расширяют возможности классических подходов к оценке опционов. Модель Хестона предоставляет глубокое понимание динамики волатильности с учётом эффекта mean reversion, однако требует значительных вычислительных ресурсов. Модель SABR, благодаря аналитическим приближениям и более простой калибровке, демонстрирует высокую точность, особенно при краткосрочных оценках. Правильный выбор модели и её настройка позволяют обеспечить адекватное отражение рыночных процессов, что критически важно для риск-менеджмента и построения эффективных торговых стратегий.

О том, как калибровать всю поверхность волатильности, поговорим в следующих статьях.

Спасибо за ваше время!

Примечания

Комментарии (2)


  1. Mox
    03.02.2025 06:31

    Я не совсем понял для чего эти модели (

    Это преодоление проблемы что в модели Блэка-Шоулза IV зависит от страйка, а не одинаковая для всех опционов в конкретный момент времени?



    1. kirrya Автор
      03.02.2025 06:31

      Привет. Да, верно, эти стохастические модели нужны для того чтобы учесть следующее:

      • Распределение цены базового актива не является логнормальным (что предполагает модель Блека-Шоулза)

      • Подразумеваемая волатильность (implied volatility), то есть та, которую нужно подставить в формулу Блека-Шоулза, чтобы получить правильную цену опциона, формирует «улыбку» - так что нам нужны такие модели, у которых в осях «страйк - подразумеваемая волатильность» будет эта улыбка. У sabr и Heston моделей такое как раз есть.