This is not an answer but an extended comment
As @Peter already observed, the fluctuations are heavy. So I tried to reproduce your statistic, hoping that the $S=$ number-of-even-steps in the transformation rule "$T(x):$ $3x+1$ when $x$ is odd, $x/2$ when $x$ is even" is what you've evaluated.
I only could go up to $10 \; 000 \; 000$ values at all, so I show the scattering in bins of size=1000 (blue), size=10000 (green), size=100000 (red). (The trendline is in red, based only on the 1000-size bins).
Different from your picture I scale the x-axis logarithmically and get a near linear trend. However, the fluctuations increase with the index of the bins - independently, how large the bins are: an interesting observation on its own.

After the fluctuations increase with the index-number of the bins, I think, a variable/an increasing binsize might be appropriate. Here I show the average number of even steps, when the size of bins are doubled with the index. Moreover, I took as x-coordinate the geometric mean value of the bin-bounds. See this:

Here we get a nice nearly linear trend. The slope $7.216$ divided by $\ln(2)$ gives about $10.4104$ which agrees roughly with the slope that has been found in the OP.
Again: the increasing fluctuations over the bins of equal size seem a remarkable effect to me.