The discovery of the Higgs boson has opened the door to the next phase of exploration in particle physics. Since the Higgs boson was found in July 2012 at CERN's Large Hadron Collider (LHC), many researches and experiments tried to measure its properties, making processes involving Higgs production interesting and important in many different ways. Monte Carlo simulations play an important role in studies of particle physics. They generate theoretical predictions and can be compared with experimental data to test our understanding of physics. With the arithmetic speed of super-computers increasing, we can use computer programs to simulate combinations of next-to-leading order (NLO) matrix elements and parton showers in Higgs production processes. In this thesis, by using a Monte Carlo event generator named Herwig 7, Higgs plus two jet production processes have been simulated with two different matrix elements, HJets++ and VBFNLO. Higgs plus two jet production was analyzed at 14 TeV, 33 TeV, and 100 TeV. This project also computed the differential cross section at the leading order and next-to-leading order with matched parton showers. This research provides a comparative study of the electroweak Higgs boson at future colliders.