We present the results of lattice Boltzmann (LB) simulations for the planar-flow of viscoplastic fluids through complex flow channels. In this study, the Bingham and Casson model fluids are covered as viscoplastic fluid. The Papanastasiou (modified Bingham) model and the modified Casson model are employed in our LB simulations. The Bingham number is an essential physical parameter when considering viscoplastic fluid flows and the modified Bingham number is proposed for modified viscoplastic models. When the value of the modified Bingham number agrees with that of the "normal" Bingham number, viscoplastic fluid flows formulated by modified viscoplastic models strictly reproduce the flow behavior of the ideal viscoplastic fluids. LB simulations are extensively performed for viscoplastic fluid flows through complex flow channels with rectangular and circular obstacles. It is shown that the LB method (LBM) allows us to successfully compute the flow behavior of viscoplastic fluids in various complicated-flow channels with rectangular and circular obstacles. For even low Re and high Bn numbers corresponding to plastic-property dominant condition, it is clearly manifested that the viscosity for both the viscoplastic fluids is largely decreased around solid obstacles. Also, it is shown that the viscosity profile is quite different between both the viscoplastic fluids due to the inherent nature of the models. The viscosity of the Bingham fluid sharply drops down close to the plastic viscosity, whereas the viscosity of the Casson fluid does not rapidly fall. From this study, it is demonstrated that the LBM can be also an effective methodology for computing viscoplastic fluid flows through complex channels including circular obstacles. (C) 2011 Elsevier B.V. All rights reserved.