I’m gonna assume you want to run 1D convolutions just in those specific bins rather than 2D convolutions around the harmonics.
If that were the case you can preset an slice operator which retrieves those harmonics, then you can run a 2D convolution with kernel size 1 for the frequency. With a proper padding you can keep the shape and then using the preset slice to “paste” back the values.