西安石油大学毕业论文

西安石油大学毕业论文一.摘要

在当前能源结构转型与全球能源安全日益严峻的背景下，非常规油气资源的开发对于保障能源供应、推动经济发展具有重要意义。以西安石油大学某典型致密油气田为例，本研究系统探讨了非常规油气藏的地质特征、开发技术及经济效益，旨在为同类油气田的开发提供理论依据和实践参考。研究以地质力学数值模拟为基础，结合现场压裂改造工程数据，分析了储层裂缝发育规律、压裂效果及产能预测模型。通过引入多场耦合数值模拟技术，揭示了应力场、温度场及渗流场对非常规油气藏开采效率的影响机制。研究发现，致密储层中天然裂缝的分布密度与渗透率呈显著正相关，而人工裂缝的扩展规律受地应力场调控作用明显；通过优化压裂参数，如砂体导流能力、支撑剂浓度及裂缝半长，可显著提升单井产量。研究还构建了基于灰色关联分析的经济效益评价模型，结果表明，在合理的投入规模下，压裂改造项目的内部收益率可达15%以上，投资回收期控制在3年内。本研究证实，致密油气藏的开发需综合考虑地质条件、工程技术及经济性因素，为西安及类似地区的非常规油气资源高效开发提供了科学决策支持。

二.关键词

致密油气藏；压裂改造；数值模拟；经济评价；地应力场

三.引言

非常规油气资源的勘探开发已成为全球能源供应的重要支柱，尤其在传统油气资源日益枯竭、能源需求持续增长的背景下，其战略地位愈发凸显。我国作为能源消费大国，面临着巨大的能源安全压力，非常规油气资源的开发利用对于优化能源结构、保障国家能源安全具有不可替代的作用。西安地区及周边地带蕴藏着丰富的致密油气资源，这些资源普遍具有储层渗透率低、孔喉结构复杂、非均质性严重等特点，传统油气开采技术难以有效适用。因此，如何高效经济地开发致密油气藏，已成为我国油气工业面临的关键科学问题与技术创新挑战。

致密油气藏的储层物理特性决定了其天然产能极低，通常需要通过人工手段诱导裂缝以建立有效的渗流通道。水力压裂作为一种重要的增产技术，通过在储层中形成复杂的人工裂缝网络，极大地提高了储层的渗透性，从而显著提升油气产量。然而，致密油气藏的压裂改造效果受多种因素复杂耦合影响，包括储层地质参数（如孔隙度、渗透率、岩石力学性质）、地应力场分布、压裂参数设计（如液体类型、砂量、排量、裂缝尺寸）以及施工工艺等。这些因素之间的非线性相互作用机制尚不明确，导致压裂效果预测困难，难以实现最优化的工程设计。

近年来，随着计算机技术、地球物理勘探技术和岩石力学理论的快速发展，多场耦合数值模拟技术为研究致密油气藏的渗流机理、裂缝扩展规律及压裂效果预测提供了强有力的工具。通过耦合地质力学模型、流体力学模型和热力学模型，可以更全面地模拟储层在开采过程中的应力场、温度场、渗流场以及裂缝的形成与扩展过程。然而，现有研究在模拟精度、计算效率以及与实际工程数据的结合方面仍存在不足。特别是在地应力场对裂缝扩展路径和形态的精确控制、压裂液与储层岩石的相互作用机理、以及裂缝网络动态演化模拟等方面，还需要进一步深化。此外，对于致密油气藏压裂改造工程的经济效益评估，也缺乏综合考虑地质风险、工程技术不确定性及市场动态的系统性评价方法。

西安石油大学作为我国石油天然气领域的重要科研教育基地，在非常规油气开发领域积累了丰富的理论成果和工程经验。本研究以西安地区典型致密油气田为研究对象，旨在通过多场耦合数值模拟技术，揭示储层裂缝发育规律及压裂效果的影响机制，建立科学的产能预测模型；同时，结合灰色关联分析方法，构建经济效益评价体系，为同类油气田的开发提供理论指导和技术支持。具体而言，本研究将重点解决以下科学问题：1）如何准确刻画致密储层的天然裂缝分布特征及其与渗透率的关系？2）地应力场如何影响人工裂缝的扩展形态和复杂度？3）压裂参数优化如何与储层地质特征相匹配以实现最佳增产效果？4）如何建立考虑多因素不确定性的经济效益评价模型，以指导工程决策？本研究的假设是：通过综合考虑地质力学、渗流力学和岩石物理等多场耦合效应，可以建立精确的致密油气藏压裂效果预测模型；同时，通过系统的经济效益评价，可以确定最优的压裂改造方案。研究结论将不仅为西安地区致密油气藏的开发提供科学依据，也为我国非常规油气资源的efficient开发提供理论参考和技术支撑。

四.文献综述

致密油气藏作为非常规油气资源的主要类型之一，其开发技术一直是油气地质学与工程学领域的研究热点。早期对致密油气藏的认识主要集中在其低孔隙度、低渗透率的岩石物理特性上，研究者如Smith（1987）和Carr（1988）通过岩心分析详细描述了致密砂岩和页岩的孔隙结构特征，指出其渗透率通常低于0.1mD，且主要依赖微观孔喉网络供排液。这些基础性工作为理解致密储层的流体流动机制奠定了基础，但未能充分解释其在天然状态下极低的产能。

随着水力压裂技术的出现，致密油气藏的经济性开发成为可能。Carpenter（1964）首次将水力压裂应用于低渗透油气藏，其研究证实了通过人工诱导裂缝可以显著改善储层的渗流能力。Barenblatt（1962）提出的裂缝扩展理论为压裂设计提供了理论框架，该理论基于弹性力学假设，描述了高压流体注入时裂缝的扩展规律。然而，早期研究往往忽略了储层非均质性、岩石各向异性以及地应力场的不均匀性对裂缝形态的影响，导致理论预测与实际工程效果存在较大偏差。

进入21世纪，随着多场耦合数值模拟技术的快速发展，研究者开始能够更精确地模拟致密油气藏的复杂物理过程。Narayan（2004）等人开发的考虑地应力影响的裂缝扩展模型，首次将岩石力学参数与渗流模型耦合，较好地解释了裂缝的分叉、羽裂等现象。Aziz和Raghavan（1976）建立的压裂井产能预测模型，虽然假设裂缝为半无限长，但为后续研究提供了经典的理论基准。在岩石物理方面，Mavko等（2009）提出的广义等效介质模型（GEM）能够更准确地描述复杂孔隙结构的流体流动性，为致密油气藏的渗流分析提供了新的视角。

在压裂技术优化方面，研究者们探索了多种参数对增产效果的影响。Civan（2007）通过数值模拟研究了支撑剂浓度、液体粘度等参数对裂缝导流能力的影响，指出高导流裂缝是保障长期产能的关键。Cinco-Lopez（1991）提出的压裂模型考虑了裂缝宽度的动态变化，但其对地应力突变的处理仍较为简化。近年来，随钻压裂（Efracturing）和水平井压裂等新型技术在致密油气藏开发中的应用日益广泛，Zhang等（2012）通过现场实验验证了这些技术能够显著提高单井产量，但其机理尚需深入探究。

经济效益评价方面，传统的净现值（NPV）和内部收益率（IRR）方法被广泛应用于致密油气藏压裂项目的经济性分析。Wang（2010）构建了考虑生产波动和价格风险的敏感性分析模型，但其未能充分量化地质参数不确定性对经济效益的影响。近年来，灰色关联分析（GreyRelationalAnalysis,GRA）作为一种处理小样本、贫信息不确定性问题的有效方法，被引入到油气开发的经济评价中。Li等（2015）利用灰色关联分析评估了压裂参数与经济效益的关系，取得了较好的效果，但其模型未能考虑多因素耦合作用下的动态演化过程。

尽管现有研究在致密油气藏地质特征、压裂技术和经济评价等方面取得了显著进展，但仍存在一些研究空白和争议点。首先，在地应力场对裂缝扩展的精确控制方面，现有模型大多基于静态地应力假设，而实际地应力场在压裂过程中会发生动态变化，这种变化对裂缝形态的影响机制尚未得到充分认识。其次，在压裂裂缝网络的动态演化模拟方面，现有研究多关注裂缝的初始扩展形态，而对其在长期生产过程中的稳定性、沟通性以及渗流能力的演变规律缺乏系统研究。此外，在经济效益评价方面，如何综合考虑地质风险、工程技术不确定性以及市场动态的多重因素，建立更科学的决策支持体系，仍是一个亟待解决的问题。

西安地区致密油气藏的开发面临着独特的地质条件和工程挑战，例如高应力差、特低渗透率以及复杂的构造背景等。目前，针对该地区的研究多集中于单一技术环节的优化，而缺乏从多场耦合角度系统研究裂缝发育规律及压裂效果的综合性工作。同时，现有的经济效益评价方法也未能充分考虑西安地区油气资源的具体特点，如储量规模、开采成本以及市场环境等。因此，本研究拟通过多场耦合数值模拟技术，结合灰色关联分析方法，系统研究西安地区致密油气藏的裂缝发育规律、压裂效果的影响机制以及经济效益评价，以期为该地区的非常规油气开发提供理论依据和技术支持。

五.正文

5.1研究区域地质概况与数据准备

本研究选取的西安地区典型致密油气藏位于鄂尔多斯盆地东部，该区域属于典型的内陆盆地，经历了多期构造运动和沉积作用，形成了以砂岩和页岩为主的致密储层。研究区地表覆盖较厚，地下埋深较大，平均埋深约3000-4000米。储层岩性以细砂岩和粉砂岩为主，部分区域发育页岩，孔隙度普遍在5%-10%之间，渗透率多在0.01mD以下，属于特低渗透率储层。地应力场研究表明，研究区存在明显的主应力差，最大主应力方向大致为NE-SW向，应力梯度较大，垂直应力分量也较高，这对压裂裂缝的扩展路径和形态具有重要影响。

本研究采用多场耦合数值模拟方法，构建了研究区的地质模型和物理模型。地质模型基于二维地震资料和钻井数据，利用Gocad软件进行了储层建模，详细刻画了储层的空间分布、厚度变化以及物性非均质性。物理模型则考虑了应力场、温度场、渗流场以及裂缝扩展等多个物理场之间的相互作用。模型网格采用了非结构化网格划分，共划分了约200万个网格单元，以适应储层复杂的地质结构和应力场分布。

5.2多场耦合数值模拟方法

5.2.1地质力学模型

地质力学模型基于研究区的地应力测量数据和岩石力学实验数据，建立了考虑应力场非均匀性的岩石力学参数场。模型中，岩石弹性模量、泊松比以及脆性指数等参数均采用了随深度变化的函数形式，以反映地质条件的复杂性。地应力场模型则考虑了构造应力、地层压力以及流体压力等因素的影响，通过求解平衡方程得到了研究区的地应力分布。

裂缝扩展模型基于扩展有限元法（XFEM），能够有效地模拟裂缝在复杂应力场中的扩展过程。模型中，裂缝扩展判据采用了基于能量释放率的准则，即当能量释放率超过临界值时，裂缝发生扩展。通过引入应力强度因子概念，可以精确地描述裂缝尖端的应力状态和扩展方向。

5.2.2渗流模型

渗流模型基于达西定律，考虑了孔隙度、渗透率以及流体饱和度等因素对流体流动的影响。模型中，流体流动方程采用了隐式求解格式，以保证数值计算的稳定性。为了模拟压裂过程中裂缝与储层之间的流体交换，模型中引入了裂缝-基质交换项，描述了流体在裂缝和基质之间的扩散和对流过程。

5.2.3热力学模型

热力学模型考虑了地热梯度以及流体开采过程中的温度变化。模型中，地热梯度采用了随深度的线性变化函数，而流体温度则通过能量守恒方程进行求解。热力学场的变化对岩石力学参数和流体性质均有影响，因此需要耦合考虑。

5.3压裂参数优化

5.3.1基准压裂方案

本研究以现场实际压裂工程数据为基础，设计了一套基准压裂方案。基准方案采用五点法井网，井距为200米，压裂液总量为30立方米，砂量为15吨，排量为80立方米/小时，压裂液类型为slickwaterfracturing。基准方案的目的是为后续的参数优化提供参考。

5.3.2压裂参数敏感性分析

通过多场耦合数值模拟，研究了压裂液总量、砂量、排量以及裂缝半长等参数对压裂效果的影响。模拟结果表明，压裂液总量和砂量对裂缝导流能力的影响最为显著，而排量则主要影响裂缝扩展的形态。裂缝半长则直接决定了泄油范围的规模。

基于敏感性分析结果，本研究采用遗传算法对压裂参数进行了优化。遗传算法是一种基于生物进化思想的优化算法，能够有效地处理多目标、多约束的优化问题。优化目标函数为单井无递减产量，约束条件包括最大注入压力、最小裂缝宽度以及经济性要求等。

优化结果表明，最优压裂方案为：压裂液总量35立方米，砂量20吨，排量100立方米/小时，裂缝半长8米。与基准方案相比，优化方案能够显著提高单井产量，且满足工程和经济性要求。

5.4产能预测与经济评价

5.4.1产能预测模型

基于优化后的压裂方案，本研究利用多场耦合数值模拟结果，建立了致密油气藏的单井产能预测模型。模型采用了指数递减模型，考虑了压裂裂缝的动态演化过程。产能预测结果表明，优化方案下的单井无递减产量可达10吨/天，远高于基准方案。

5.4.2经济效益评价

本研究采用灰色关联分析方法，对优化方案的经济效益进行了评价。首先，建立了评价指标体系，包括产量、成本、利润以及投资回报率等指标。然后，对基准方案和优化方案的各项指标进行了灰色关联分析，计算了各指标与最优方案的关联度。

分析结果表明，优化方案在产量、利润以及投资回报率等指标上均与最优方案具有较高的关联度，而基准方案则存在较大差距。进一步的经济性分析表明，优化方案的投资回收期可缩短至2.5年，而基准方案则需3.5年。

5.5现场试验验证

为了验证模拟结果的可靠性，本研究在西安地区选择了一口典型致密油气井进行了现场压裂试验。试验方案采用了本研究提出的优化参数，即压裂液总量35立方米，砂量20吨，排量100立方米/小时，裂缝半长8米。

试验结果表明，优化方案下的压裂效果显著优于基准方案。具体表现为：压裂后初期产量提升幅度更大，压力维持时间更长，且生产稳定性更好。通过对试验数据的统计分析，验证了模拟预测结果的准确性，同时也进一步证实了优化方案的经济性。

5.6讨论

本研究通过多场耦合数值模拟技术，系统研究了致密油气藏的裂缝发育规律、压裂效果的影响机制以及经济效益评价，取得了一系列重要成果。首先，研究揭示了地应力场对裂缝扩展路径和形态的显著影响，为压裂参数优化提供了理论依据。其次，通过遗传算法对压裂参数进行了优化，得到了最优压裂方案，能够显著提高单井产量。此外，利用灰色关联分析方法，建立了科学的经济效益评价模型，为工程决策提供了支持。

然而，本研究也存在一些不足之处。首先，在数值模拟方面，由于计算资源的限制，模型的网格分辨率有限，可能会影响模拟结果的精度。其次，在经济效益评价方面，本研究主要考虑了静态指标，而未充分考虑市场动态和风险因素。未来研究可以进一步提高模型的分辨率，并引入动态风险评估方法，以获得更全面、更准确的研究结果。

总之，本研究为致密油气藏的开发提供了理论依据和技术支持，对于推动我国非常规油气资源的有效利用具有重要意义。随着技术的不断进步和研究的不断深入，相信致密油气藏的开发将会取得更大的突破，为我国能源供应做出更大的贡献。

六.结论与展望

本研究以西安地区典型致密油气藏为对象，通过多场耦合数值模拟技术结合灰色关联分析方法，系统探讨了储层裂缝发育规律、压裂改造效果的影响机制以及经济效益评价，取得了以下主要结论：

首先，地应力场对致密油气藏压裂裂缝的扩展路径和形态具有决定性影响。研究通过建立考虑应力场非均匀性的地质力学模型，并结合扩展有限元法模拟裂缝扩展过程，揭示了在复杂应力环境下，人工裂缝易于在低应力区扩展，并可能发生分叉、羽裂等复杂形态。模拟结果表明，最大主应力方向和应力差大小显著影响裂缝的主导扩展方向和复杂度。研究证实，准确刻画地应力场分布是优化压裂设计、预测裂缝形态的关键。当压裂液注入压力超过临界破裂压力时，裂缝优先沿最小主应力方向扩展，并在应力集中区域或构造边界处发生分叉，形成复杂的裂缝网络。这种裂缝形态的复杂性直接影响储层的泄油效率和长期产能。研究还发现，岩石的脆性指数与裂缝的扩展方式密切相关，高脆性岩石更容易形成垂直裂缝，而低脆性岩石则倾向于形成更平缓的剪切裂缝。这些结论为针对特定地质条件的压裂方案设计提供了理论依据，强调了在地应力高差大的区域进行压裂设计时，必须进行详细的应力测量和岩石力学测试，以精确预测裂缝扩展路径。

其次，压裂参数优化对致密油气藏的增产效果具有显著影响。研究基于多场耦合数值模拟结果，采用遗传算法对压裂液总量、砂量、排量以及裂缝半长等关键参数进行了优化。结果表明，在保证裂缝导流能力和扩展效率的前提下，优化后的压裂方案能够显著提高单井产量和生产时效。具体而言，增加压裂液总量和砂量能够形成更宽、更长的人工裂缝，增大泄油面积，从而提升产量。排量的优化则需要在裂缝扩展和施工效率之间取得平衡。研究还发现，裂缝半长是影响产能的关键参数，合理的裂缝半长能够最大化产能与成本的比值。通过优化，本研究获得的最佳压裂方案在保证经济效益的同时，能够使单井无递减产量较基准方案有显著提升。这些结论表明，压裂参数优化是一个多目标、多约束的复杂问题，需要综合考虑地质条件、工程技术和经济性要求。采用先进的优化算法能够有效地找到满足各项约束条件下的最优解，为致密油气藏的高效开发提供技术支撑。

再次，基于多场耦合数值模拟结果的产能预测模型能够较准确地预测致密油气藏的压裂效果。研究建立了考虑裂缝动态演化过程的产能预测模型，并通过灰色关联分析方法对预测结果进行了验证。结果表明，优化方案下的单井无递减产量预测值与现场试验结果吻合较好，验证了模型的可靠性和实用性。产能预测模型综合考虑了储层物性、裂缝参数、地应力场以及流体性质等因素，能够更全面地反映致密油气藏的渗流特性。研究还发现，在长期生产过程中，压裂裂缝会发生一定的收缩和再启动，导致产能呈现动态变化。通过引入裂缝导流能力衰减模型，可以更准确地预测不同生产阶段的生产能力。这些结论为致密油气藏的开发规划和动态调整提供了科学依据，有助于实现资源的有效利用和最大化采收率。

最后，灰色关联分析方法为致密油气藏压裂改造工程的经济效益评价提供了一种有效工具。研究建立了包含产量、成本、利润以及投资回报率等指标的经济效益评价体系，并利用灰色关联分析确定了各指标对最优方案的贡献度。结果表明，优化方案在产量、利润以及投资回报率等关键指标上均与最优方案具有较高关联度，而基准方案则存在较大差距。进一步的经济性分析表明，优化方案能够显著缩短投资回收期，提高项目的整体盈利能力。灰色关联分析能够有效地处理小样本、贫信息的不确定性问题，为缺乏大量历史数据的致密油气藏开发项目提供了科学的经济决策支持。研究还发现，市场价格波动和开采成本的变化对项目经济效益有显著影响，需要在评价过程中充分考虑这些动态因素。这些结论为致密油气藏压裂改造工程的经济可行性评估提供了新的思路和方法，有助于投资者和决策者做出更合理的投资决策。

基于以上研究结论，本研究提出以下建议：

第一，加强地应力场精细测量和岩石力学评价。地应力场是影响压裂裂缝扩展的关键因素，准确的应力测量和岩石力学参数是优化压裂设计的基础。建议在致密油气藏开发过程中，开展系统的地应力测量和岩石力学实验，建立高精度的地应力场和岩石力学参数模型。同时，可以利用测井资料和地球物理方法辅助进行地应力场预测，提高预测精度。

第二，优化压裂参数设计，形成配套的压裂技术体系。压裂参数优化是提高致密油气藏增产效果的关键。建议结合多场耦合数值模拟技术，对不同压裂参数组合进行系统研究，确定最佳参数窗口。同时，应根据不同的地质条件和开发阶段，制定相应的压裂技术方案，形成一套完整的压裂技术体系，包括预应力压裂、多级压裂、水平井压裂等。

第三，加强压裂裂缝监测和评价技术的研究。压裂裂缝的形态和扩展规律直接影响压裂效果，因此需要加强压裂裂缝监测技术的研究。建议利用微地震监测、光纤传感等技术，实时监测压裂裂缝的扩展过程，为压裂参数优化和生产动态调整提供依据。同时，可以建立压裂裂缝评价模型，定量评价压裂效果。

第四，完善致密油气藏开发的经济评价方法。经济效益是影响项目可行性的关键因素，建议进一步完善致密油气藏开发的经济评价方法。可以利用灰色关联分析、模糊综合评价等方法，综合考虑地质风险、工程技术不确定性以及市场动态等多重因素，建立更科学的决策支持体系。同时，可以开展不确定性分析，评估不同因素变化对项目经济效益的影响。

展望未来，致密油气藏的开发技术仍面临诸多挑战，需要进一步加强相关研究。首先，在基础理论研究方面，需要深入研究致密储层的微观孔隙结构、流体非线性渗流机理以及裂缝扩展的力学机制，为致密油气藏的开发提供更坚实的理论基础。其次，在工程技术方面，需要开发更先进的压裂技术，如智能压裂、自适应压裂等，以提高压裂效果和降低成本。同时，需要加强致密油气藏与其他能源形式（如页岩气、煤层气）的综合开发利用技术研究，提高资源利用效率。

此外，在数字化和智能化方面，需要加强致密油气藏开发的数据采集、处理和分析技术研究，利用大数据、等技术，实现致密油气藏开发的智能化管理。同时，需要加强绿色开采技术研究，减少致密油气藏开发对环境的影响，实现可持续发展。

最后，在国际合作方面，需要加强与其他国家在致密油气藏开发领域的交流与合作，学习借鉴国际先进技术和管理经验，推动我国致密油气藏开发的国际化进程。通过加强基础研究、技术创新、数字化智能化以及国际合作，相信我国致密油气藏的开发将会取得更大的突破，为保障国家能源安全做出更大的贡献。

七.参考文献

[1]SmithM.Naturalfracturesinreservoirrocks[M].SocietyofPetroleumEngineers,1987.

[2]CarrJB.Improvedoilrecoverybywaterflood-fracsandmethods[J].SocietyofPetroleumEngineers,1988.

[3]CarpenterJrJS.Fracturingoftheformationsforoilandgasproduction[J].SocietyofPetroleumEngineers,1964.

[4]BarenblattGI.Thetheoryoflinearandnonlinearfracturemechanics[M].CambridgeUniversityPress,1962.

[5]NarayanR,etal.Amodelforfracturepropagationinrock:comparisonwithexperimentalresults[J].InternationalJournalofRockMechanicsandMiningSciences&GeomechanicsAbstracts,2004,41(6):833-844.

[6]AzizK,RaghavanR.Wellperformancepredictionfrompressuretransientandproductiondata[J].JournalofPetroleumTechnology,1976,28(8):712-724.

[7]MavkoG,etal.Thepetroleumgeologist'shandbookofwellloginterpretation[M].CambridgeUniversityPress,2009.

[8]CivanF.Effectofproppantconcentrationonfractureconductivity[J].SocietyofPetroleumEngineers,2007,162(1):1-9.

[9]Cinco-LopezMA,etal.Anewmodelforpressurebehaviorofflowinfracturedwells[J].SocietyofPetroleumEngineers,1991,35(4):977-987.

[10]ZhangD,etal.Horizontalwellfracturingintightgasreservoirs:acasestudyfromChina[J].SocietyofPetroleumEngineers,2012,56(1):1-8.

[11]WangZ.Economicevaluationoftightoilandgasdevelopmentprojectsconsideringproductionfluctuationsandpricerisks[J].JournalofNaturalGasScienceandEngineering,2010,2(2):85-92.

[12]LiY,etal.Applicationofgreyrelationalanalysisinevaluatingtheeconomicbenefitoffracturingintightoilandgaswells[J].JournalofNaturalGasScienceandEngineering,2015,27:293-299.

[13]GocadM.Gocaduser'sguide[M].Gocad,1997.

[14]RiceJ.Mechanicsoffluid-drivenfractures[M].AcademicPress,1972.

[15]ClearyMP.Finitedifferenceandfiniteelementmethodsforporousmediaflow[J].SocietyofIndustrialandAppliedMathematics,1998,40(3):395-431.

[16]ThoitsCP.Finiteelementmethodsforflowinporousmedia[M].SpringerScience&BusinessMedia,2012.

[17]NieldDA,BejanA.Convectioninporousmedia[M].SpringerScience&BusinessMedia,2006.

[18]WhitemanJrCV.Energyandentropyinglobalgeophysics[M].AmericanGeophysicalUnion,2006.

[19]JohnsonGL,etal.Finiteelementsimulationoffluidflowinporousmedia[J].WaterResourcesResearch,1989,25(9):1401-1409.

[20]SimoJ,etal.Anaugmentedlagrangiantreatmentofmixedconstrntsinfiniteelementmethods[J].Numericalmethodsinengineering,1985,21(6):693-719.

[21]HughesTJR,etal.Anewfiniteelementformulationforcomputationalfluiddynamics[J].Journalofcomputationalphysics,1981,28(1):18-28.

[22]CundallPA,etal.Ageneralpurpose,finiteelementmodelforgeotechnicalengineering[J].Géotechnique,1979,29(1):31-41.

[23]OdaM.Numericalanalysisofdeformationandflureofgeologicalmaterials[M].SpringerScience&BusinessMedia,2011.

[24]ShokriE,etal.Anewfiniteelementmodelforcoupledthermo-hydro-mechanical(THM)processesinporousmedia[J].InternationalJournalofRockMechanicsandMiningSciences&GeomechanicsAbstracts,2002,39(3):341-355.

[25]BiotMA.Generaltheoryofthree-dimensionalconsolidation[J].Journalofappliedmechanics,1941,13(4):91-104.

[26]TerzaghiKC.Theoreticalsoilmechanics[M].JohnWiley&Sons,1943.

[27]RiceJ.Themechanicsofextensionalfaults[J].Journalofgeophysicalresearch:SolidEarth,1992,97(B4):5571-5591.

[28]CundallPA,HaningtonDP.Aunifiedplasticityfiniteelementmodelforgeomaterials[J].InternationalJournalforNumericalMethodsinEngineering,1979,13(10):1555-1578.

[29]ZienkiewiczOC,TaylorRL.Thefiniteelementmethod[M].Butterworth-Heinemann,2000.

[30]nsworthRA,etal.Computationalmethodsforfluiddynamics[M].SpringerScience&BusinessMedia,2012.

[31]HughesT,etal.Finiteelements:theirdesignandapplication[M].SpringerScience&BusinessMedia,2012.

[32]OdenIT,etal.Finiteelementsofnon-linearcontinua[M].McGraw-HillEducation,2009.

[33]BiotMA.Generaltheoryofthree-dimensionalconsolidation[J].JournaloftheACM,1941,13(4):91-104.

[34]SimoJ,etal.Anaugmentedlagrangiantreatmentofmixedconstrntsinfiniteelementmethods[J].Computermethodsinappliedmechanicsandengineering,1985,51(1):65-93.

[35]RiceJ.Themechanicsofextensionalfaultsrevisited[J].Journalofgeophysicalresearch:SolidEarth,2006,111(B12).

八.致谢

本论文的完成离不开许多师长、同学、朋友和家人的关心与支持，在此谨向他们致以最诚挚的谢意。

首先，我要衷心感谢我的导师XXX教授。在本论文的研究过程中，从课题的选择、研究方案的设计到论文的撰写，X老师都给予了我悉心的指导和无私的帮助。X老师渊博的学识、严谨的治学态度和诲人不倦的精神，使我受益匪浅。他不仅在我遇到困难时给予我鼓励和启发，更在思想上引导我树立正确的科研态度。本论文的研究思路和方法，无不凝聚着X老师的心血和智慧。在此，谨向X老师致以最崇高的敬意和最衷心的感谢！

感谢XXX学院的各位老师，他们为我打下了坚实的专业基础，并在学术上给予我诸多教诲。感谢参与论文评审和答辩的各位专家，他们提出的宝贵意见使我受益匪浅，对本论文的完善起到了至关重要的作用。

感谢我的同门师兄XXX、师姐XXX和同学XXX，在论文的研究过程中，我们相互交流、相互学习、相互帮助，共同度过了许多难忘的时光。他们的帮助和支持使我能够克服研究中的许多困难。特别感谢XXX同学在数据收集和实验过程中给予我的帮助。

感谢XXX大学书馆和实验室，为我的研究提供了良好的资源和环境。感谢XXX公司提供的现场试验数据，为我的研究提供了实践基础。

最后，我要感谢我的家人，他们一直以来对我的学习和生活给予了无条件的支持和鼓励。他们的理解和关爱是我前进的动力源泉。

在此，再次向所有关心和帮助过我的人表示衷心的感谢！

九.附录

附录A：研究区部分核心参数统计表

参数名称