For regional studies and practical use regional-scale models of the geomagnetic field are necessary. The Bulgarian Geomagnetic Reference Field (BulGRF) for 2015.0 epoch and prediction up to 2020.0 is produced and presented in this paper. It uses data from Panagyurishte (PAG) observatory and the Bulgarian secular network. As a result, our model gives a set of six coefficients which allow the users to obtain the main field values using only the geographic coordinates and a simple calculator.