Motivation. I was thinking in the more simple case of (6), I am saying the case $r=s=1$, that tell us this MathWorld. Then in this case one can check that the identity holds using integration methods. When I did it I've considered from Hata's integral the following different integral in this
Question. What's about $$\int_0^1\int_0^1\frac{\log(xy)xy}{-1+\log(xy)}dxdy?$$
Remark. The context was that I was trying with my imagination different expressions for an hypothetical similar (6) but now for $\sum_{k=r+1}^\infty\frac{\mu(n)}{n^3}$, where $\mu(n)$ is the Möbius function. For example I tried also consider $\int_0^1\left(\int_1^\infty\left(\int_1^\infty\frac{dx}{xy(1-(1-xy)z)}\right)dy\right)dz$, but I understand that it is science fiction to find an expression for $\sum_{k=r+1}^\infty\frac{\mu(n)}{n^3}$ from this random method.
My attempt. I can deduce the following statement $$\int_0^1\int_0^1\frac{\log(xy)xy}{-1+\log(xy)}dxdy=\frac{1}{4}+\int_0^1\int_0^1\frac{dxdy}{-1+\log(xy)},$$ and I know that $$\int\frac{dx}{-1+\log(xy)}=e\frac{\operatorname{Ei}(\log(xy)-1)}{y}+\operatorname{ constant}.$$ Then can we finish the example to get the integral in the Question in terms of particular values of special functions? Thanks in advance.
References:
Hata, A New Irrationality Measure for $\zeta(3)$, Acta Arith. 92, 47-57 (2000).