Given this product:
$$a(n) = \prod_{k=1}^{n} (1+p_k)$$
where $p_k$ is the $k$-th prime number and which can be interpreted also as the sum of the divisors of the $n$-th primorial (OEIS A054640), is there a way to get an asymptotic estimate for $n \to \infty$ of the $p_k$-adic valuation of $a(n)$:
$$\nu_{p_k}(a(n))$$?
As expected, lower primes seem more probable:
\begin{array}{cccc} n&\nu_2(a(n))&\nu_3(a(n))&\nu_5(a(n))\\ 1&0&1&0\\ 2&2&1&0\\ 3&3&2&0\\ 4&6&2&0\\ 5&8&3&0\\ 6&9&3&0\\ 7&10&5&0\\ 8&12&5&1\\ 9&15&6&1\\ 10&16&7&2 \end{array}
After @Denis Shatrov first comment, and some numerical tests for $n \le 1000000$ here are some guesses:
$$\nu_2(a(n)) \approx 2n$$
$$\nu_3(a(n)) \approx \frac{3}{4}n$$
$$\nu_5(a(n)) \approx \frac{5}{16}n$$
$$\nu_7(a(n)) \approx \frac{7}{36}n$$
Using more standard notation: for any fixed prime $q$, \begin{align*} v_q\biggl( \prod_{p\le x} (p+1) \biggr) &= \sum_{p\le x} v_q(p+1) \\ &= \sum_{p\le x} \sum_{\substack{k\ge 1 \\ q^k \mid (p+1)}} 1 \\ &= \sum_{k\ge1} \sum_{\substack{p\le x \\ p\equiv-1\pmod{q^k}}} 1 \\ &= \sum_{k\ge1} \pi(x;q^k,-1). \end{align*} Since $\pi(x;m,a) \sim \frac1{\phi(m)} \frac x{\log x}$ for every $a$ with $\gcd(a,m)=1$, we deduce that $$ v_q\biggl( \prod_{p\le x} (p+1) \biggr) \sim \sum_{k\ge1} \frac1{\phi(q^k)} \frac x{\log x} = \sum_{k\ge1} \frac1{q^{k-1}(q-1)} \frac x{\log x} = \frac{q}{(q-1)^2} \frac x{\log x}. $$