We consider the space \(\textrm{D}^\alpha _n(\Omega )\) consisting of functions u(x), harmonic in a bounded domain \(\Omega \subset \mathbb {R}^{d+1}\) with smooth boundary \(\Gamma \) , satisfying \(|\nabla ^n u(x)|^2 \rho (x)^{\alpha }\in L_1(\Omega )\) , \(\alpha >-1\) , where \(\rho (x)\) is the distance from x to the boundary. For a Borel measure \(\mu \) on \(\Gamma \) and a weight function V, \(\mu \) -measurable, we study the operator defined by means of the quadratic form \(\pmb {\mu }_V[u]=\int V(x) |\gamma u(x)|^2 \mu (dx)\) , where \(\gamma \) is, properly defined, operator of restriction of functions \(u\in \textrm{D}^\alpha _n(\Omega )\) to \(\Gamma \) . Main interest is directed to the case of a singular measure \(\mu \) possessing some Ahlfors regularity properties. For such operators, we establish two-sided estimates for singular values, and, under some geometrical conditions, Weyl type eigenvalue asymptotics.