С использованием уравнений Био в частотной области в квазистатической постановке моделируется низкочастотное нагружение образца трещиновато-пористой горной породы. Уравнения Био аппроксимируются конечно-разностной схемой на сдвинутых сетках. Для решения полученной системы линейных алгебраических уравнений (СЛАУ) с числом неизвестных более 106 используется итерационный метод (стабилизированный алгоритм бисопряженных градиентов) с предобуславливателем на основе расщепления полей для разделения уравнений и переменных на две группы: описывающих деформацию твердого тела и переноса жидкости. Численные эксперименты демонстрируют быструю сходимость итерационного процесса и его преимущество по сравнению с прямым методом решения СЛАУ на больших задачах.