This study deals with computational modeling of defect clustering effects observed in bilayer phospholipid membranes. Two defect clustering models (algorithms) are presented and compared with the random defect distribution approach. Specific defect distribution instances are evaluated using a simple methodology based on Voronoi diagrams and statistical properties of their sector areas. Computational experiments are performed by using the models to generate synthetic defect distributions with different parameter combinations. The proposed methodology is also validated against atomic force microscopy images of real membranes with defects.