Abstract: Based on a sample of 2.7×10<sup>9 ψ(3686) events collected by the BESIII detector operating at the BEPCII collider, an analysis of the decay ψ(3686)→γχcJ,χcJ→ηηη<sup>′ is performed. The decay modes χc1 and χc2→ηηη<sup>′ are observed for the first time, and their corresponding branching fractions are determined to be B(χc1→ηηη<sup>′)</sup>=(1.40±0.13(stat.)±0.09(sys.))×10<sup>−4 and B(χc2→ηηη<sup>′)</sup>=(4.18±0.84(stat.)±0.48(sys.))×10<sup>−5. An upper limit on the branching fraction of χc0→ηηη<sup>′ is set as 2.59×10<sup>−5 at 90\% confidence level (CL). A partial wave analysis (PWA) of the decay χc1→ηηη<sup>′ is performed to search for the 1<sup>−+ exotic state η1(1855). The PWA result indicates that the structure in the ηη<sup>′ mass spectrum is mainly attributed to the f0(1500), while in the ηη mass spectrum, it is primarily the 0<sup>++ phase space. The upper limit of $\mathcal{B}(\chi_{c1}\to\eta_{1}(1855)\eta) \cdot \mathcal{B}(\eta_{1}(1855)\to\eta\eta<sup>{\prime})<</sup> 9.79 \times 10<sup>{-5}$ is set based on the PWA at 90\% CL.