Results of a statistical variation of total ion density observed in the vicinity of epicenters as well as around magnetically conjugated points of earthquakes are presented in this paper. Two data sets are used: the ion density measured by DEMETER during about 6.5 years and the list of strong earthquakes (MW ≥ 4.8) occurring globally during this period (14,764 earthquakes in total). First of all, ionospheric perturbations with 23–120 s observation time corresponding to spatial scales of 160–840 km are automatically detected by a software (64,287 anomalies in total). Second, it is checked if a perturbation could be associated either with the epicenter of an earthquake or with its magnetically conjugated point (distance < 1500 km and time < 15 days before the earthquake). The index Kp < 3 is also considered in order to reduce the effect of the geomagnetic activity on the ionosphere during this period. The results show that it is possible to detect variations of the ionospheric parameters above the epicenter areas as well as above their conjugated points. About one third of the earthquakes are detected with ionospheric influence on both sides of the Earth. There is a trend showing that the perturbation length increases as the magnitude of the detected EQs but it is more obvious for large magnitude. The probability that a perturbation appears is higher on the day of the earthquake and then gradually decreases when the time before the earthquake increases. The spatial distribution of perturbations shows that the probability of perturbations appearing southeast of the epicenter before an earthquake is a little bit higher and that there is an obvious trend because perturbations appear west of the conjugated point of an earthquake.