I am going to try to answer the question, because of what is said in the comments:
It does not seem to me that the 7 answers properly answer the question! For that reason I reopened the question.
Before delving into the explanation I think it is important to note that Orbital Angular Momentum is a physical quantity that is first defined in Classical Mechanics. Its Quantum Mechanical analogue is obtained by a (reasonable) guess that the expression should have the same components, and simply translating the classical variables to quantum operators. But it did not arise directly from the symmetry principles of QM.
Why can't $\hat{\mathbf{L}}^2$ have half-integer quantum numbers?
We first encounter orbital angular momentum in the context of the solution to the Hydrogen atom while looking for a Complete Set of Commuting Observables to express the bound states of the system.
It is defined, in the same way as the Hamiltonian, by analogy with the classical angular momentum as:
$$\hat{\mathbf{L}}\equiv\hat{\mathbf{X}}\times\hat{\mathbf{P}}=-i\hbar\mathbf{x}\times\boldsymbol{\nabla}$$
where the last equality holds in the position and momentum representation of the state vector $\Psi$. The components of $\hat{\mathbf{L}}$ can be written as:
$$\hat L_i=-i\hbar\sum\limits_{j,k}\varepsilon_{ijk}x_j\frac{\partial}{\partial x_k}$$
where $\varepsilon_{ijk}$ is the Levi-Civita Symbol. The total orbital angular momentum, the subject of your question is then given by $\hat{\mathbf{L}}^2$,
$$\hat{\mathbf{L}}^2=\sum\limits_i\hat L_i\hat L_i=-\hbar^2\sum\limits_{i,j,k,l,m}\varepsilon_{ijk}\varepsilon_{ilm}x_j\left(\frac{\partial}{\partial x_k}\right)x_l\left(\frac{\partial}{\partial x_m}\right)$$
And skipping some steps we arrive at:
$$\hat{\mathbf{L}}^2=-\hbar^2\left(r^2\nabla^2-\frac{\partial}{\partial r}r^2\frac{\partial}{\partial r}\right)$$
Expressed in spherical coordinates $(r, \theta, \phi)$ (the total angular momentum has to be rotationally invariant, so it is expected that the operator is only a function of the $r$ coordinate).
For a reasonable$^{*}$ potential $V(\mathbf x)$ we can express $\Psi$ near the origin as a power series in the basis vectors. In particular, as $r\to 0$, only the terms with the lowest order in the series "survive":
$$\Psi(\mathbf x) \to r^{\ell}Y(\theta, \phi)$$
Here $\ell$ is the order of the lowest non-vanishing term in the power series. The $Y(\theta, \phi)$ is also a $\ell$'th order (homogeneous) polynomial in $\{\mathrm{sin}(\theta)\mathrm{cos}(\phi), \ \mathrm{sin}(\theta)\mathrm{sin}(\phi), \ \mathrm{cos}(\theta)\}$ (you may identify these as the $\{x, \ y, \ z\}$ components of a unitary vector).
Acting on this $\Psi$ with $\hat{\mathbf{L}}^2$ we get,
$$\hat{\mathbf{L}}^2\Psi \to \hbar^2\ell(\ell+1)\Psi$$
Now, if $\Psi$ is an eigenfunction of $\hat{\mathbf{L}}^2$ then the eigenvalue should be the same for any and every point, so we conclude that the eigenvalue of such a function must be $\hbar^2\ell(\ell+1)$.
This is as far as we can go with angular momentum in the position-momentum representation. As we saw above all (reasonable) wavefunctions, including the eigenfunctions of $\hat{\mathbf{L}}^2$ can be expressed as $\propto r^\ell$ for small $r$. And there are no fractional power series$^{**}$, it follows that the orbital angular momentum only has (positive) integer quantum numbers.
So, to respond your question,
Does anyone have a decisive argument on why do we exclude the half-integer values from the orbital operator spectrum?
Because the orbital angular momentum was chosen to match our classical definition of angular momenta. This leads to the above expression of $\hat{\mathbf{L}}^2$ with its corresponding limitations.
At least that's how I interpret it.
How do half-integer angular momenta appear?
As we have seen, the position-momentum representation only gives us integer angular momenta, for more general angular momenta we need to rely on symmetry principles.
Symmetries are transformations that leave the equations of motion, and the physical content of the mathematical representations invariant. In classical mechanics the symmetries are encoded in the Galilean group: space rotations, translations in space and time and galilean transformations should not change the outcome of an experiment (if I perform an experiment under the same conditions it shouldn't matter if I do it today or tomorrow, or if I do it in Hansford or Livingston).
In quantum mechanics, all the transformations acting on the state vectors should be Unitary, in order to preserve the inner product (and therefore the probabilities). That is,
$$\langle\hat{U}\Psi\vert\hat{U}\Phi\rangle = \langle\Psi\vert\Phi\rangle$$
It follows that $\hat{U}^{\dagger}\hat{U}=\hat{1}$. For infinitesimal transformations (for example, rotation by an infinitesimal angle), the Unitary operators take the following form:
$$\hat U = \hat{1}+i\epsilon \hat T+\cdots$$
where $\epsilon$ is the "infinitesimal" (arbitrarily small) number. Here, the operator $\hat T$ is hermitian. In this context $\hat T$ is called the generator of the transformation (you could repeat the infinitesimal transformation an infinite number of times to get a finite transformation and the only information you need to do that is $\hat T$). You can use this to show that the momentum operator is the generator of space translations, and the Hamiltonian is the generator of time translations.
If you know what form $\epsilon$ takes (it is related to the transformation in question) you can exploit the unitary condition $\hat{U}^{\dagger}\hat{U}=\hat{1}$ to arrive at the commutator relations for the operator $\hat T$. The particular case of the angular momentum ($\hat T =\hat{\mathbf{J}}$) is shown in chapter 4, section 1 of Weinberg's Lectures on Quantum Mechanics. I will omit the details of the derivation, but what you arrive at is:
$$\left[\hat J_i, \hat J_j\right]=i\hbar\sum\limits_k\varepsilon_{ijk}\hat J_k$$
As you can check, $\hat{\mathbf{L}}$ also satisfies the commutation relations$^{***}$. In general we can express $\hat{\mathbf{J}}$ as,
$$\hat{\mathbf{J}}=\hat{\mathbf{L}}+\hat{\mathbf{S}}$$
We have already seen that $\hat{\mathbf{L}}$ has just integer angular momenta, so the only place for the half-integers is $\hat{\mathbf{S}}$. Additionally, with this you can check that $\left[\hat S_i, \hat L_j\right]=0$, and also, $$\left[\hat S_i, \hat X_j\right]=\left[\hat S_i, \hat P_j\right]=0$$
What this last equation means is that $\hat{\mathbf{S}}$ is independent of $\hat{\mathbf{X}}$ and $\hat{\mathbf{P}}$, and this is why we can't see the spin in the position-momentum representation.
In the symmetry part I omitted a lot of things not to make the answer too extensive. The source for most of my answer is chapters 2,3 and 4 of:
Weinberg, S. (2015). Lectures on Quantum Mechanics. Cambridge: Cambridge University Press.
$^{*}$ In the time independent Schrödinger equation we have $\frac{\hbar^2}{2m}\nabla^2\Psi=(V-E)\Psi$, $E$ is finite but $V$ may diverge (like in the Hydrogen atom), but if the divergence is not very singular then the first derivatives of $\Psi$ will be at least finite and $\Psi$ itself will be continuous everywhere (and differentiable everywhere except possibly at the singular point)
$^{**}$ For instance $r^{\frac{1}{2}}$ doesn't even have a series expansion at $r=0$ whereas higher order half-integers of the form $a+\frac{1}{2}$, for $a\in\mathbb N$, have their first $a$ derivatives at $r=0$ equal to zero (so it is a trivial expansion) and their $(a+1)$'th derivative diverge (so the expansion is undefined).
$^{***}$ There is an important difference. For any vector operator $\hat{\mathbf{V}}$, the commutation relation $\left[\hat J_i, \hat V_j\right]=i\hbar\sum\limits_k\varepsilon_{ijk}\hat V_k$ holds, but this is not true if we change $\hat{\mathbf{J}}$ for $\hat{\mathbf{L}}$.