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