Typical seizure prediction models aim at discriminating interictal brain activity from pre-seizure electrographic patterns. Given the lack of a preictal clinical definition, a fixed interval is widely used to develop these models. Recent studies reporting preictal interval selection among a range of fixed intervals show inter- and intra-patient preictal interval variability, possibly reflecting the heterogeneity of the seizure generation process. Obtaining accurate labels of the preictal interval can be used to train supervised prediction models and, hence, avoid setting a fixed preictal interval for all seizures within the same patient. Unsupervised learning methods hold great promise for exploring preictal alterations on a seizure-specific scale. Multivariate and univariate linear and nonlinear features were extracted from scalp electroencephalography (EEG) signals collected from 41 patients with drug-resistant epilepsy undergoing presurgical monitoring. Nonlinear dimensionality reduction was performed for each group of features and each of the 226 seizures. We applied different clustering methods in searching for preictal clusters located until 2 h before the seizure onset. We identified preictal patterns in 90% of patients and 51% of the visually inspected seizures. The preictal clusters manifested a seizure-specific profile with varying duration (22.9 ± 21.0 min) and starting time before seizure onset (47.6 ± 27.3 min). Searching for preictal patterns on the EEG trace using unsupervised methods showed that it is possible to identify seizure-specific preictal signatures for some patients and some seizures within the same patient.