Quite some time ago I made some test with colored flats for OSC. Reason being I was concerned about the noise in the red channel. I'm using an Aurora EL flat by Gerd Neumann and that thing is really blue. Neumann says this is intended to match the color of the sky. The result is a very weak signal in red in an otherwide well exposed OSC image. So I ordered 1/8 orange Lee foil, the sort that is made for daylight to light bulb conversion. I opend the flat panel and stacked as many as needed between the actual EL foil and the diffusor until the color balance was better. It is still not perfect but the red level is reasonably higher than it was before, using longer exposures of course. I have never seen any kind of color gradient since I use the flat panel. Like you I alway to background extraction at function degree 1, subtraction, because for the small patch of the sky that my telescopes can see the sky gradient can be approximated by a linear function. I think you should clearly distinguish between vignetting correction which is done by good flats during the calibration process and gradient removal in post processing. I have processed a lot of other CNer's data who do not care about flat at all. When there is a lot of nebulosity in the image PI will 'correct' the nebulosity alongside with the vignetting when background removel is used (abused?) for correction of vignetting.
Mark, how do you measure the sky background? From images or do you have a sky quality meter? I'd like to know what bortle I really have. Clearoutside estimates 4 for my location. I'm surprise you find 3 in the UK. I doubt there is any Bortle 3 nearby for me. Sigh.