Abstract: In this paper the computation method of geoid variation with time by repeated level-ling and gravimetry is discussed by the authors.The optimization method of positional parameters of fictitious point masses is presented.The fictitious point mass method over-comes the disadvantage of stokes integral and collocation.Simulated computation shows that the computed geoid variation is very accurate.