本文通过三维有限元数值模拟的方法,对伊朗kashan地区Qom组A-F段的构造应力场进行了有限元建模和元素应力应变分析。其有限元线性代数方程为:在上述七个方程的理论基础上,学习并运用ANSYS软件结合已有的伊朗kashan地区Qom组的地学模型和力学模型,采用ANSYS软件中的3维8节点六面体单元模拟kashan地区Qom组A-F段的构造应力场,采用solid45实体单元模拟,对地质模型、力学模型、边界力的作用方式及约束条件、断层的处理和数学模型的建立均进行了合理简化,从而建立了kashan地区Qom组A-F段的构造应力场三维有限元模型,并进行了有限元模态分析,得到了不同层段上的应力与应变的情况。与已有的应力场的情况基本吻合,为储层裂缝定量预测提供了参考依据本文以伊朗kashan地区Qom组A-F段断层,背斜,向斜作为计算模型的主要目的层及研究对象,采用细网格的三维立体单元,有限元模型模拟范围:设模型长1200m,宽1200m,纵深500m,考虑三个主断裂系统,垂向上划分六个岩性段(A—F)盆地内岩性横向变化较快,尤以垂直盆地走向(NW—SE)NE—SW为甚。每个单元按单元组为之单独编号,共划分了2080个单元,2172个节点。在模型底部施加Z方向约束;模型的南(下方)边界,即在局部坐标的Y方向约束;模型东(右方)边界为局部坐标系X方向约束;垂向上加载自身重力;近东南向或近北北东向加载80.34Mpa或85.76Mpa。此外,本文利用三维有限元模拟应力场得出应力—应变,根据岩石的破裂准则,通过地层综合破裂值评价指标(Fy)和岩石形变能量值建立了裂缝密度的数学模型,利用多元统计方法的反复讨论和试算,得出了裂缝的不同层段的平均密度,为储层裂缝的分布规律的预测提供了充分的准备,也为后续的研究工作打下了良好的基础。