This will answer your questions on $\delta$:
The case $\delta=1$, sometimes referred to as the ideal or optimal LLL-Algorithm, is the most natural way to define the swap condition between vectors after the length reduction step, as it forces to swap $b_k$ and $b_{k+1}$ if and only if its orthogonal 'emergence', i.e.
orthogonal distance of $b_{k+1}$ with respect to the linear span of the $b_1,\ldots,b_{k-1}$ is smaller than the orthogonal emergence of $b_k$.
Some people refer this type of reduction to Hermite, who applies this algorithm in order to prove an upper bound for the Hermite constant in any dimension $d$.
Lovász introduced the relaxation factor $\left(\frac{1}{2}\right)^2< \delta<1$ for the swap condition
$d(b_{k+1},Span(b_1,\ldots,b_{k-1}))<\sqrt{\delta}\cdot d(b_{k},Span(b_1,\ldots,b_{k-1}))$,
in order to be able to prove an upperbound for the time complexity of the algorithm which is polynomial in the magnitude of the initial basis AND the dimension $d$.
Also in the ideal case $\delta=1$ the algorithm is still polynomial with respect to the magnitude of the basis (A. Akhavi, The optimal LLL algorithm is still polynomial in fixed dimension,Theoretical Computer Science 297 (2003),1-3, pp.3–23), but it is still an open problem whether it is also polynomial in the dimension $d$.
The factor $\delta$ can be regarded as a strictness factor: the closer it is to $1$, the more orthogonal the reduced basis will be in general, and therefore both its Hermite factor as well as its approximation factor for the minimal lattice norm $\lambda_1$ are better.
I don't know of any strict theoretical proof of this connection, but extensive experiments as made by Gama, Nguyen and Stehlé (cf. P. Nguyen, D. Stehlé, LLL on average, Algorithmic Number Theory 4076 (2006), pp.1-17 or N. Gama, P. Nguyen, Predicting lattice reduction, Proc. of the EUROCRYPT 2008, Vol. 4965, pp.31-51) show that for LLL with $\delta=0.999$ the average Hermite factor, taking samples from an equidistributed (according to a natural measure on lattices described in D.Goldstein, A.Mayer, On the equidistribution of Hecke points, Forum Mathematicum (2003) 15, 2, pp. 165-189) sequence of integer lattices, increases exponentially by
$ 0.801\cdot 1.0219^d$,
which shows much better performance than predicted by the theoretic upper bound $(4/3)^\frac{d}{4}\approx 1.074^d$.
Using what is called Gaussian heuristics, one can take the Hermite factors to estimate the average approximation factor by
$average \frac{\|b_1\|}{\lambda_1}\approx\sqrt{\frac{\pi\cdot e}{2\cdot d}}\cdot 0.801\cdot 1.0219^d$,
slightly less exponential in $d$.
According to a frequently cited draft of M.Ajtai, Generating random lattices according to the invariant distribution (2006), this heuristics is quite accurate for high dimension lattices.
Unfortunately, I have never found the draft itself or any later on published paper.
However, if you don't want to rely on Gaussian heuristics, the bounds I derived from my own experiments running LLL with $\delta=1$ on integer lattices from dimension $30$ to $100$, choosing random knapsack type bases of magnitude $2^{20d}$, are
$0.3444\cdot 1.0136^d\leq average \frac{\|b_1\|}{\lambda_1} \leq 0.6688\cdot 1.0442^d,$
the lower bound representing the Gaussian heuristics, again much better than the proven upper bound $(4/3)^{d/2}\approx 1.1547^d$.
For less ideal $\delta$, I don't know any published estimate of the average behaviour.
The experiments I've done by myself seem to verify the thumb rule: the smaller $\delta$, the faster reduction is done, but the more volatile the quality of the results, and moreover the stronger is the exponential growth of the average Hermite and approximation factors.