Smoothening Noisy FFT

I've got the following FFT with a sampling frequency of 192 kHz that has around 17.000.000 samples. The signal is a hiss of escaping gas. This looks noisy and i want to determine the characteristic "hiss frequency" that should lie between ±40 kHz (according to other studies).

Since most of the noise in industrial areas is broadband, and the hissing characteristics can be observed accross a wide frequency band including lower frequencies (e.g. from 0-10 kHz an), I can filter out anything below 20 kHz. This is done with a bandpass designed with the help of scipy's cookbook on Butterworth. I've used a 30th order Bandpass. Below is the resulting FFT.

Using a "noise-only" (leakage-free) measurement I have determined the FFT of the noise level.

Now, according to the noise level 3, the noise should exist mainly under 20 kHz. However, this can change with the environment such that the leakage can be "consumed" by noise. How can I solve this issue?

I have plotted the PSD 4 of the non-filtered data using Welch's algorithm with a Tukey window (which is apparently good for transient data). The plot can be examined below.

Am I using the correct window for this? How can I smoothen the FFT and reduce the noise, such that I get a better understanding of the overall data, since I have other measurements that i have to compare with each other?