arXiv:2610.33415v1 Announce Type: new Abstract: We study the uplink of a user-centric (UC) cell-free massive MIMO system assisted by multiple beyond-diagonal reconfigurable intelligent surfaces (BD-RISs), under imperfect channel state information (CSI) and pilot reuse. The BD-RISs are modeled as reciprocal, lossless, group-connected scattering networks. We first derive the minimum mean-square error (MMSE) estimation statistics for each user's aggregate channel to its serving access points (APs). To avoid frequent CSI acquisition and RIS reconfiguration, we optimize the BD-RIS scattering matrices using long-term channel statistics, allowing their configurations to remain fixed over multiple coherence blocks. Specifically, we minimize the sum of normalized mean-square errors (NMSEs) of the MMSE aggregate channel estimators, subject to the symmetry and unitarity constraints of the group-connected architecture. We develop a product-manifold Riemannian optimization algorithm and analyze its convergence and computational complexity. To evaluate uplink spectral efficiency (SE), we separately derive a closed-form expression for the effective SINR under maximum ratio (MR) combining at the central processing unit (CPU) and establish a lower bound on the achievable ergodic SE. The SINR derivation accounts for the statistical dependence among channels associated with different APs due to shared BD-RIS-assisted propagation paths, including cross-AP second-order statistics under pilot reuse. Numerical results show that moderate group connectivity retains most of the NMSE improvements achieved by full connectivity, with larger optimization gains as BD-RIS-assisted paths strengthen relative to direct AP-user links. They also reveal that lower estimation NMSE does not necessarily lead to higher uplink SE.