Using the Biot equation in scope of the frequency domain and quasi-static state, we can model low-frequency loading and obtain strain-stress relations that vary with frequency for fluid-filled poroelastic materials. To solve the linear algebraic equations resulting from the finite difference discrectisation of the Biot equation, we propose an approach based on the BCSGStab iterative solver and a preconditioner technique based on a field-split approach. In this study, we calculate the approximate amount of floating point operations needed for the suggested technique as well as a direct method. Despite the direct solver's superior OMP scalability over the iterative one, numerical experiments show that the iterative technique quickly converges and is advantageous for big Biot problems.