Maximum cyclostationarity blind deconvolution (CYCBD) can recover the periodic impulses from mixed fault signals comprised by noise and periodic impulses. In recent years, blind deconvolution has been widely used in fault diagnosis. However, it requires a preset of filter length, and inappropriate filter length may cause the inaccurate extraction of fault signal. Therefore, in order to determine filter length adaptively, a method to optimize CYCBD by using the seagull optimization algorithm (SOA) is proposed in this paper. In this method, the ratio of SNR to kurtosis is used as the objective function; firstly, SOA is used to search the optimal filter length in CYCBD by iteration, and then it uses the optimal filter length to perform CYCBD; finally, the frequency-domain waveform is determined through Fourier transformation. The method proposed is applied to the fault extraction of a simulated signal and a test vibration signal of the closed power flow gearbox test bed, and the fault frequency is successfully extracted, in addition, using maximum correlation kurtosis deconvolution (MCKD) and multipoint optimal minimum entropy deconvolution adjusted (MOMEDA) to compare with CYCBD-SOA, which validated availability of the proposed method.