本发明公开了一种面向三维地质表面模型的四面体网格化方法，用于地质体属性建模、油藏数值模拟计算及分析。其特征在于，在对输入三维地质表面网格模型的规范化整理基础上，基于点集的德劳内三角网格生成方法，利用拓扑学原理，通过检测约束条件与德劳内四面体网格的对偶维诺图中边、面、体等元素的相交情况及建立相应加点规则，提供了一种自动可靠的地质约束边界线段及层面网格的恢复方法，生成了在输入三维地质构造模型约束下的四面体网格；方法简单易于实现，可靠性好，对于初始表面模型的复杂程度和网格类型无特殊要求；解决了复杂地质构造下的四面体网格化问题，便于在此基础上对三维地质模型进行进一步的分析与模拟计算。