$\chi^2$ is a valid statistic for any data where you expect the errors to be normally distributed. It is commonly used for counting data, where you expect poisson errors of order $\sqrt{O_i}$ (which become gaussian in the large count limit), but as @a_statistician notes the more general formula is
$$\chi^2 = \sum_i\frac{(O_i-E_i)^2}{Var(O_i)}$$.
Restating your problem, it sounds like you have some model parameters, let's call them $\alpha$ and $\beta$, which are used to predict $E_i(\alpha,\beta)$, and some observed values $O_i$, along with some prediction about the size of errors, $Var(O_i)$. Therefore you can calculate $\chi^2(\alpha,\beta)$ as a function of $\alpha$ and $\beta$.
You have then (correctly) estimated $\alpha$ and $\beta$ by minimizing this function, and now want to obtain the confidence intervals for these estimates. As the paper you linked explains, one way to do this is to consider the parameters themselves to be free, and lookup to corresponding cumulative-probability function for a chi-squared distribution with that number of degrees of freedom (using one of the many handy tables available online).
For the example you gave, with two parameters, the 95% cumulative probability point is 5.991 (corresponding to (2,0.05) on the linked table). Therefore, one would call the 95% confidence region the set of all points $(\alpha, \beta)$ in parameter space such that
$$\chi^2(\alpha,\beta) - \chi^2_{min} < 5.991$$
For a good, well-constrained model, this should give you something fairly nice, like an ellipse. This ellipse corresponds to the covariance matrix for your parameters, as this article nicely explains, and therefore standard errors and other commonly used measures of uncertainty.