Full history match models in subsurface systems are challenging due to the large number of reservoir simulations required, and the need to preserve geological realism in matched models. This drawback increases significantly in big real fields due to the high heterogeneity of the geological models, the reservoir simulation computational time (which increases superlinearly). In this work, we propose a novel framework based on artificial intelligence to address these shortcomings. Our workflow is based on two main components: The first is the new combination of model order reduction techniques (e.g., principle component analysis (PCA), kernel-PCA (k-PCA)) and artificial intelligence for parameterizing complex three-dimensional (3D) geomodels, called "Geo-Net". Our new approach is able to create complex high dimensional heterogeneous reservoirs in seconds, providing better correspondence with the underlying geomodels, hard-data constraints and geological plausibility. The second component is a derivative-free optimization framework to complete the automatic history matching (AHM). This new approach allows us to perform local changes in the reservoir at the same time as we conserve geological plausibility. We have examined our methodology in a real field in Colombia. The Rubiales Oil Field is located in the Llanos Basin with original oil in place of around 6 billion barrels. The key finding here is that the Geo-Net is able to recreate the full geological workflow obtaining the same high order of statistics as traditional geostatistical techniques. Nonetheless, our Geo-Net allows us to control the full process with a low-dimensional vector and reproduces the full geological workflow 10, 000 times faster than commercial geo-statistical packages. Finally, the full optimization workflow has been applied to AHM. Results show an improvement with respect to best practice of traditional history match workflows.