Assume an IID sample of the form $ \left\{ y^{r}_{i},\mathbf{x}_{i} \right\}$ (notice the superscript on $y$). The observed values $y^{r}_{i}$ are bounded from above by the actual, unobserved values of $y_{i}$. The values that both $y_{i}$ and $y_{i}^{r}$ obtain are in $\mathbb{Z}_{+}$ (consider the following for instance: the former represent crimes counts in some neighborhood and the latter, reported crimes counts). I postulate the following simple linear model: $$y_{i}=\mathbf{x}_{i}\cdot \beta +\epsilon_{i1}$$ and assume the observed values come from the following: $$y^{r}_{i}=\mathbf{x}_{i}\cdot \beta -a_{i}\cdot \gamma +\epsilon_{i2}$$ while $a_{i}$ is unobserved.
Is there any way to consistently estimate $\beta$? I have read about Heckman's correction but it does not seem applicable here since the observed values are not observed after some threshold. I have failed to develop anything useful by myself. Ideas and suggestions welcome.