As a brain disorder, epilepsy is characterized with abnormal hypersynchronous neural firings. It is known that seizures initiate and propagate in different brain regions. Long-term intracranial multichannel electroencephalography (EEG) reflects broadband ictal activity under seizure occurrence. Network-based techniques are efficient in discovering brain dynamics and offering finger-print features for specific individuals. In this study, we adopt link prediction for proposing a novel workflow aiming to quantify seizure dynamics and uncover pathological mechanisms of epilepsy. A dataset of EEG signals was enrolled that recorded from 8 patients with 3 different types of pharmocoresistant focal epilepsy. Weighted networks are obtained from phase locking value (PLV) in subband EEG oscillations. Common neighbor (CN), resource allocation (RA), Adamic-Adar (AA), and Sorenson algorithms are brought in for link prediction performance comparison. Results demonstrate that RA outperforms its rivals. Similarity, matrix was produced from the RA technique performing on EEG networks later. Nodes are gathered to form sequences by selecting the ones with the highest similarity. It is demonstrated that variations are in accordance with seizure attack in node sequences of gamma band EEG oscillations. What is more, variations in node sequences monitor the total seizure journey including its initiation and termination.