ScaleResfusionResidual Rectified Flow based on Residual Vector Field

Zhenning Shi1,*, Chen Xu2,*, Junhao Zhang3, Kefei Zhang1, Linjie Liu4, Zhedong Zheng2,†, Tao Li1,†
1Nankai University  ·  2University of Macau  ·  3CSIRO Data61  ·  4Beihang University
*Equal contribution    †Corresponding authors
Real-world low-quality inputs (top) and ScaleResfusion restorations (bottom)
Top: real-world low-quality inputs. Bottom: ScaleResfusion restorations with the FLUX.2-klein-4B backbone and only 4 sampling steps.

Abstract

Real-world Image Restoration (Real-IR) aims to recover high-quality (HQ) images from complex and unknown degradations. Recent diffusion-based methods have substantially improved perceptual quality, yet two obstacles remain: methods that sample from Gaussian noise require many steps and are often less faithful to the degraded input, whereas residual-based methods that start from the low-quality (LQ) image typically train task-specific models from scratch, with optimization objectives coupled to a particular noise scheduler, and therefore cannot reuse modern pre-trained generative priors.

We present ScaleResfusion, which rewrites residual restoration as a scheduler-independent adaptation interface for pre-trained text-to-image rectified-flow models. Its core, Residual Rectified Flow (RRF), inserts the residual term R into the linear transport path of Rectified Flow, so that sampling starts from noisy LQ at an exact acceleration point, where the signal-to-noise ratio of the starting state is continuously controlled by the residual ratio γ. The resulting optimization target, the residual vector field, contains no scheduler-specific coefficients and differs from the pre-trained rectified-flow target only by the residual offset γR; adapting a frozen billion-scale backbone therefore reduces to fitting this compact residual correction with LoRA-only training. A knowledge-distillation pipeline built around RRF further reduces sampling to as few as 4 steps. Experiments on real-world super-resolution across multiple benchmarks show that ScaleResfusion achieves state-of-the-art restoration quality and transfers consistently across pre-trained rectified-flow backbones from 2B to 9B parameters.

Residual Rectified FlowStart sampling from the noisy LQ image at the exact acceleration point t* = 1/(1+γ), with the SNR of the start controlled by γ.
Scheduler-independent targetThe residual vector field resv = vRF + γR differs from the pre-trained objective only by a residual offset.
LoRA-only adaptationFrozen billion-scale backbones (SD3, FLUX.2, Z-Image) are adapted with ~60M trainable parameters.
4 steps, 646 msState-of-the-art fidelity and distribution alignment on the real-world DRealSR / RealSR benchmarks, and consistently competitive on DIV2K-Val / LSDIR-Val, at sub-second latency.

Method

Residual Rectified Flow

Given an LQ–HQ pair, we define the residual R = x̂0 − x0. RRF inserts the weighted residual into the linear path of Rectified Flow while keeping the transport straight:

xt = (1 − t) x0 + tγR + tε,  resv = γR + ε − x0.

Substituting R shows that the coefficient of the unknown HQ image vanishes at t* = 1/(1+γ): the state there depends only on the LQ image and Gaussian noise, so inference starts from xt* = γ/(1+γ) x̂0 + 1/(1+γ) ε and integrates the learned residual vector field back to t = 0 over a shortened interval. Because resv = vRF + γR carries no scheduler-specific coefficients, adapting a pre-trained rectified-flow model only requires fitting the compact residual correction, which we realise with LoRA alone.

Geometric interpretation of Residual Rectified Flow
(a) RRF (purple) intersects the implicit noise-to-LQ path at the acceleration point. (b) Standard Rectified Flow transports Gaussian noise to HQ images. (c) RRF starts from the noisy-LQ state and recovers the HQ image over a shortened sampling interval.

Why start from the residual?

The residual term is compact in both the spatial and frequency domains
The residual R is a low-energy, edge-like correction in both the spatial and frequency domains, so RRF only needs to learn a compact update around the pre-trained transport — a natural fit for LoRA.

The reverse ODE is initialization-sensitive: a start closer to the restoration manifold is transported to a terminal state closer to the HQ distribution. Gaussian initialization is task-agnostic and far from the observed content; deterministic one-step LQ→HQ mappings lose the sampling capability of diffusion models; noisy-LQ DDPM sampling is biased toward the degradation without an explicit correction direction.

ScaleResfusion instead starts from a noisy LQ state and follows a residual-oriented transport, giving a closer initialization and more coherent convergence toward the target distribution.

Initialization comparison from a transport perspective
(a) Gaussian initialization starts far from the restoration manifold. (b) One-step restoration collapses into a deterministic LQ→HQ mapping. (c) Noisy-LQ DDPM sampling is biased toward the input degradation. (d) ScaleResfusion starts from a noisy LQ state and follows a residual-oriented transport.

Knowledge-distilled parameter-efficient training

The restored-image generator predicts the residual vector field under RRF. LQ features are injected as conditioning (a ReferenceNet-style feature extractor for SD3 / Z-Image, token concatenation for FLUX.2), and DAPE provides textual constraints. A DMD loss and an optional GAN loss regularise the restored distribution toward natural HQ images. Generation and regularization branches share the same frozen backbone weights; only LoRA adapters are trained.

ScaleResfusion training pipeline
Overview of the training pipeline with a generation module and a regularization module.

Results

Quantitative comparison with the default FLUX.2-klein-4B backbone and 4 sampling steps. Red = best, blue = second best. The first group lists GAN-based methods, the second group diffusion-based methods.

MethodPSNR ↑SSIM ↑LPIPS ↓DISTS ↓FID ↓NIQE ↓MUSIQ ↑MANIQA ↑
BSRGAN28.700.800.290.21155.616.5457.150.48
Real-ESRGAN28.610.810.280.21147.666.7054.270.49
LDL28.200.810.280.21155.517.1453.940.49
FeMaSR26.870.760.320.22157.725.9153.700.44
StableSR28.040.750.330.23144.156.6058.530.56
SUPIR25.090.650.420.28169.487.3958.790.55
TSD-SR27.770.760.300.21134.985.9166.620.59
AddSR26.680.740.370.26164.827.8065.360.60
CCSR28.240.780.320.23157.306.8166.280.61
DiffBIR25.900.620.470.29180.336.3366.130.62
OSEDiff27.920.780.300.22135.416.4664.690.59
PASD28.020.780.320.23174.766.7257.230.51
ResShift27.050.740.390.26159.908.6551.240.47
SeeSR28.070.770.320.23147.376.4165.090.61
FluxSR20.880.610.390.25144.667.1166.380.57
Ours (w/o GAN)29.770.820.250.20118.186.9562.420.61
Ours (w/ GAN)28.260.780.290.21124.036.2165.160.64
MethodPSNR ↑SSIM ↑LPIPS ↓DISTS ↓FID ↓NIQE ↓MUSIQ ↑MANIQA ↑
BSRGAN26.380.770.270.21141.245.6463.280.54
Real-ESRGAN26.650.760.270.21136.295.8560.450.55
LDL25.280.760.280.21142.745.9960.920.55
FeMaSR25.060.740.290.23141.015.7759.050.49
StableSR24.620.700.310.22128.545.7865.480.62
SUPIR23.650.660.350.25130.386.1162.090.58
TSD-SR24.810.720.270.21114.455.1371.190.63
AddSR22.650.650.380.27154.186.6271.410.67
CCSR25.920.750.280.21122.845.7369.180.64
DiffBIR24.830.650.360.24130.755.8469.280.65
OSEDiff25.150.730.290.21123.535.6569.080.63
PASD26.040.740.280.21135.485.7160.030.56
ResShift25.660.740.330.25128.038.0756.890.51
SeeSR25.150.720.300.22125.305.4069.810.65
FluxSR23.830.690.320.23120.636.5770.070.61
Ours (w/o GAN)27.050.780.240.20104.746.2167.250.64
Ours (w/ GAN)25.780.750.260.20106.225.2969.860.68
MethodPSNR ↑SSIM ↑LPIPS ↓DISTS ↓FID ↓NIQE ↓MUSIQ ↑MANIQA ↑
BSRGAN24.580.630.350.2349.554.7561.680.50
Real-ESRGAN24.020.640.320.2138.874.8360.380.54
LDL23.830.630.330.2242.284.8660.040.53
FeMaSR22.450.590.340.2241.974.8757.940.48
StableSR23.270.570.310.2024.954.7765.780.62
SUPIR22.130.530.390.2331.405.6863.860.59
TSD-SR23.020.580.270.1829.164.3271.690.62
AddSR22.370.560.380.2334.915.8469.150.63
CCSR24.300.630.300.2030.845.3469.530.61
DiffBIR23.140.540.370.2232.714.9969.870.64
OSEDiff23.720.610.290.2026.344.7167.960.61
PASD24.010.610.380.2237.064.9863.750.55
ResShift24.590.620.310.2130.816.9258.900.53
SeeSR23.680.600.320.2025.894.8168.660.62
Ours (w/o GAN)24.900.640.270.1822.905.0966.100.63
Ours (w/ GAN)23.630.600.260.1719.474.4469.300.68
MethodPSNR ↑SSIM ↑LPIPS ↓DISTS ↓FID ↓NIQE ↓MUSIQ ↑MANIQA ↑
BSRGAN20.820.540.250.1646.374.2168.940.63
Real-ESRGAN20.580.550.240.1541.284.1869.520.64
LDL20.310.530.250.1644.754.3668.610.64
FeMaSR19.870.510.270.1748.634.0967.850.61
StableSR20.310.550.310.1854.765.0762.960.61
SUPIR20.350.500.240.1543.814.8171.470.67
TSD-SR19.050.490.210.1445.663.8674.450.68
AddSR19.200.450.340.2079.804.9974.200.70
CCSR20.760.530.260.1656.284.2572.570.66
DiffBIR20.510.490.270.1658.454.4473.260.68
OSEDiff20.390.520.270.1659.574.0372.340.66
PASD20.930.520.310.1759.113.8069.290.62
ResShift21.230.550.230.1438.985.3268.560.61
SeeSR20.690.520.250.1552.064.1073.270.68
Ours (w/o GAN)21.350.570.230.1441.694.2170.610.67
Ours (w/ GAN)20.640.550.210.1443.563.8272.590.69

Efficiency

Number of function evaluations (NFE) and average inference time per 512×512 image on a single NVIDIA RTX A6000.

MethodStableSRSUPIRSeeSRResShiftCCSROSEDiffOurs
NFE200505015614
Inference time (ms)12,03625,2524,445848516266646

Consistent transfer across backbones

RRF adapts pre-trained rectified-flow backbones that differ in architecture and scale — SD3 (2B), FLUX.2-klein (4B / 9B) and Z-Image (6B) — with the same LoRA-only recipe. Results on DRealSR; the upper block is w/o GAN, the lower block w/ GAN.

BackbonePSNR ↑SSIM ↑LPIPS ↓DISTS ↓FID ↓
SD3 (2B)28.770.790.300.23149.67
FLUX2-Klein (4B)29.770.820.250.20118.18
Z-Image (6B)29.350.800.280.22132.31
FLUX2-Klein (9B)30.170.820.250.19109.28
SD3 (2B), w/ GAN27.860.760.320.22146.93
FLUX2-Klein (4B), w/ GAN28.260.780.290.21124.03
Z-Image (6B), w/ GAN28.340.780.280.21127.40
FLUX2-Klein (9B), w/ GAN28.770.790.280.19116.69
Consistent transfer of ScaleResfusion across backbones
Fidelity and perceptual quality improve from SD3 (2B) to FLUX2 (9B); FLUX2-4B offers the best quality–efficiency trade-off.

Few-step sampling

Separately trained 4-, 2- and 1-step FLUX.2-4B models on DRealSR. More steps improve fidelity, while even the 1-step model remains competitive thanks to the RRF acceleration point.

NFEPSNR ↑SSIM ↑LPIPS ↓FID ↓
4-step29.770.820.25118.18
2-step29.030.790.27115.43
1-step28.930.780.28117.56
Intermediate restoration results of the multi-step model
Intermediate predictions of a single 4-step model: early steps recover the global structure, later steps refine details.

Visual Comparison

Visual comparison with existing diffusion-based restoration methods
Left: ScaleResfusion result. Right: local crops of LQ, ResShift, StableSR, CCSR, AddSR, SeeSR, SUPIR, OSEDiff, ScaleResfusion and HQ.

User Study

Pairwise model arena: in each round a participant sees an LQ input with the anonymised results of two models and picks the preferred one (or tie / both bad). 48 evaluators contributed 1,756 pairwise comparisons over benchmark images and in-the-wild photos.

Pairwise arena user study results
Left: ScaleResfusion wins 71.8% of cross-model comparisons. Right: it tops the Random Arena ranking with a 75.8% overall score.

Video

A short video introduction of ScaleResfusion is available online.

Watch the video

BibTeX

@article{shi2026scaleresfusion,
  title   = {ScaleResfusion: Residual Rectified Flow based on Residual Vector Field},
  author  = {Shi, Zhenning and Xu, Chen and Zhang, Junhao and Zhang, Kefei and Liu, Linjie and Zheng, Zhedong and Li, Tao},
  journal = {arXiv preprint arXiv:2607.25275},
  year    = {2026}
}