The presence of pharmaceuticals in drinking water, even at very low concentrations, has raised concerns among stakeholders such as drinking-water regulators, governments, water suppliers and the public, with regard to the potential risks to humans. Despite this, the occurrence and the fate of pharmaceuticals in drinking waters of many countries (e.g. in Poland) remains unknown. There is a lack of sufficiently sensitive and reliable analytical methods for such analyses and a need for more in-depth hydrogeological analysis of the possible sources of drug residues in drinking water. In this paper, a multi-residual method for the simultaneous determination of seventeen human pharmaceuticals in drinking waters has been developed. Large-volume extractions using Speedisk extraction disks, and derivatization prior to GC-MS-SIM analysis using a new silylating agent DIMETRIS were applied. The method detection limits (MDLs) ranged from 0.9 to 5.7 ng/L and the absolute recoveries of the target compounds were above 80% for most analytes. The developed method was successfully applied in the analysis of the target compounds in drinking water collected in Gdansk (Poland), and of the 17 pharmaceuticals, 6 compounds were detected at least once. During the investigation, the geomorphology of the site region was taken into account, possible sources of pharmaceuticals in the analysed drinking water samples were investigated, and the presence of the drugs in ground and surface waters, raw and treated drinking waters was determined. Concentrations were also compared with those observed in other countries. As a result, this study has not only developed a new analytical method for determining pharmaceuticals in drinking waters as well as rendering missing information for Poland (a country with one of the highest consumptions of pharmaceuticals in Europe), but it also presents a modelled in-depth hydrogeological analysis of the real sources of drugs in drinking waters.