Geologic surface approximation is profoundly affected by the presence, density and location of scattered geologic input data. Many studies have recognized the importance of utilizing varied sources of information when reconstructing a surface. This paper presents an improved geologic surface approximation method using a multiquadric function and borehole data. Additional information, i.e., inequality elevation and dip-strikes data extracted from outcrops or mining faces, is introduced in the form of physical constraints that control local changes in the estimated surface. Commonly accepted hypothesis states that geologic surfaces can be approximated to any desired degree of exactness by the summation of regular, mathematically defined, surfaces: in particular displaced quadric forms. The coefficients of the multiquadric functions are traditionally found by a least squares method. The addition of physical constraints in this work makes such an approach into a non-deterministic polynomial time problem. Hence we propose an objective function that represents the quality of the estimated surface and that includes the additional constraints by incorporation of a penalty function. Maximizing the smoothness of the estimated surface and its fitness to the additional constraints then allows the coefficients of the multiquadric function to be obtained by iterative methods. This method was implemented and demonstrated using data collected from the 81’st coal mining area of the Huaibei Coal Group. Geologic surface approximation is profoundly affected by the presence, density and location of scattered geologic input data. Many studies have recognized the importance of utilizing varying sources of information when reconstructing a surface. This paper presents an improved geologic surface approximation method using a multiquadric function and borehole data. Additional information, ie, inequality elevation and dip-strikes data extracted from outcrops or mining faces, is introduced in the form of physical constraints that control local changes in the estimated surface. Commonly accepted hypothesis states that geologic surfaces can be approximated to any desired degree of exactness by the summation of regular, mathematically defined, surfaces: in particular displaced quadric forms. The addition of physical constraints in this work makes such an approach into a non-deterministic polynomia l time problem. Therefore we propose an objective function that represents the quality of the estimated surface and that includes the additional constraints by incorporation of a penalty function. Maximizing the smoothness of the estimated surface and its fitness to the additional constraints then allows the coefficients of the multiquadric function to be obtained by iterative methods. This method was implemented and demonstrated using data collected from the 81’st coal mining area of ​​the Huaibei Coal Group.