Pre-existing models for thermal history modelling have shown deficiency in explicit algorithms to establish the quantitative relationship between maturity indices and thermal gradients in some sedimentary basins that experienced multi-episodic rifting evolution. In this study, a forward and inverse combination model(FICM) is proposed to estimate the vitrinite reflectance(Ro) and thermal gradients. The forward module is used to calculate Ro values. It couples the EASY%Ro model with burial history reconstruction with consideration of thermal gradient variations during basin evolution. The inverse module reconstructs histoical thermal gradients by calibrating cmputed Ro against measured Ro data. The time-temperature series is a necessary input for both forward and inverse modules. Sample density is a profound factor influencing the accuracy of modelling results. In order to obtain satisfying outputs, a sufficient sample density is required. Thermal gradients are assumed to vary linearly between two given samples. Modelling results of case studies indicate that the sensitivity of heating time to Ro evlution is differnt with thermal gradients depending on geolgoical setting. Three difffernt districts, which include the time-sensitive area, the temperature-sensitive area and the non-sensitive area, can be recognized on the the relationship map among Ro variations, heating time and geothermal gradients. This model can be applied to reconstruct the thermal history and maturation evolution in a basin that has undergone complex multi-episodic rifting.
