有限元网格生成是数值分析的首要环节,也是应用地学机理与过程模型模拟和分析等地学计算的必要前提,网格质量直接影响数值分析和地学计算的精度和效率,因此有限元网格生成方法是一个长久不衰的热点课题。在众多的网格单元类型中,由于六面体网格相对于四面体网格在网格数量、计算精度和抗畸变性上具有明显的优势,而得到了广泛的使用。然而,六面体网格生成方法的研究却不如四面体成熟,六面体网格生成依然存在许多技术难题尚待解决。传统的六面体网格生成方法较少考虑复杂特征约束,基于复杂特征约束的六面体网格生成算法在设计上有所不足。而事实上,特征约束是地学研究对象中实际存在的一类重要的地物,对复杂特征约束的支持也是一个优秀的六面体网格生成算法必备的要素之一。本文对基于复杂特征约束生成高质量六面体网格这一核心问题上展开研究。通过特征约束的分类与映射生成特征约束数据集,作为六面体网格生成的前置数据。基于该前置数据,本文研究了结构化网格和非结构化网格的生成方法,给出了特征约束处理模板,在此基础上进行网格质量优化方法的探究,以生成顾及复杂特征约束的六面体网格。本文的主要研究内容与成果如下：(1)特征约束的配置与映射机制。基于传统特征约束的概念和地学对象自身的特征,本文给出了六面体网格剖分特征约束的定义,将特征约束按维度、获取方式以及对网格的影响等特征进行分类,通过特征约束的映射与匹配生成网格剖分前置数据集,将特征约束作为网格剖分的数据源,使得特征约束的处理与网格生成紧密相连。(2)基于特征约束的六面体网格生成方法。针对地表面复杂的特征约束,综合考虑各类网格剖分方法的特点,在传统扫掠法的基础上增加特征约束处理机制,以生成基于特征约束的结构化六面体网格。通过对比栅格法和前沿法的特点,在吸收二者优势的基础上抵消各自的不足,对地学研究区域内部复杂特征约束进行处理,提出基于特征约束的非结构化六面体网格生成方法。(3)特征约束的处理方案。针对特征约束数据量较大,空间分布不均匀的特点,对六面体网格中出现的各种特征约束(点、线、面、体)进行后置处理,主要通过模板和表面拟合的方法,对前文无法处理的复杂的特征约束进行处理,实现特征约束的精确拟合,生成能够满足特征约束分布的有限元网格。针对特征约束处理过程中出现的网格节点偏移与畸变现象,对网格单元质量进行优化,使用拓扑关系优化和几何形状优化的方法提高网格单元的质量。