本发明涉及一种基于多点地质统计的地质属性预测方法，包括以下步骤：导入各个实际地质勘探点的观测数据，观测数据至少包括每个实际地质勘探点的三维坐标及相应的地质属性；输入查询点的三维坐标，根据Voronoi算法或基于欧氏距离的包络圆算法来判别所述查询点的空间邻接关系图；确认查询点的邻接点，邻接点为落入所述空间邻接关系图内的实际地质勘探点；计算出每个邻接点处的三维地质属性的概率分布；采用加权算法及权值函数对所有邻接点的概率分布加权求和，得到查询点处的三维地质属性概率分布函数；依据概率分布函数，确定查询点处三维地质属性的极大似然估计并得出输入点的地质属性。该方法原理简单、修改维护方便，便于一线从业人员使用。