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