In this research, simultaneous modeling of geoid data, topography, basic concepts of physics and mathematics and local isostasy will be used to determine the thickness of the crust and the lithosphere model in geologically active areas, which is always one of the most interesting topics of researchers. Present method not only has good accuracy and high compatibility with geological phenomena, but also has high speed and lower cost compared to other methods. Other advantages of this method include the possibility of simultaneously studying vast geographical areas and comparing several different geological structures together. It is good to know that due to the desired initial results and high speed, the present method can be used as a primal model in complex inversions that require a preliminary model to start the inversion.