Spatial metapopulation dynamics in point pattern approach has been stochastically simulatedand analytically described by the dynamics of singlets,  pairs and triplets.  A metapopulationis  assumed  to  includes  constant  number  of  patches.   Each  patch  is  represented  as  a  point  inthe two-dimensional unit space Ω = [0,1)×[0,1) assumed to be torus with periodic boundaryand is assigned a status either 0 (empty) or 1 (occupied).  We assume that the point patternincludesNpoints with static positions and three types of point pattern in used are completespatial randomness, clumped and over-dispersed.  In our model, “colonization rate” with whichan empty patch becomes occupied plays a key role.  Colonization rate is assumed to consist oflocal  and  global  colonization,  in  which  local  colonization  occurs  with  the  rate  that  decreaseswith distance and global colonization rate occurs at the rate proportional to the total number ofoccupied patches independent of distance.  The dynamics of the metapopulation model can bemathematically described by the dynamics of singlet, pair and triplet probabilities.  The resultsof simulations and analytical derivation of equilibrium singlet and pair probabilities have beenderived and compared.