基于数值模拟的致密油注二氧化碳驱油效能与机理探究

基于数值模拟的致密油注二氧化碳驱油效能与机理探究一、引言1.1研究背景与意义随着全球经济的快速发展，能源需求持续增长，石油作为重要的能源资源，其供应安全与可持续性备受关注。在常规石油资源逐渐减少的情况下，非常规石油资源的开发利用成为了全球能源领域的研究热点。致密油作为一种重要的非常规石油资源，因其储量丰富、分布广泛，在全球能源格局中占据着愈发重要的地位。致密油是指赋存于致密砂岩、致密碳酸盐岩等储层中的石油，其储层具有低孔隙度、超低渗透率的特点，开采难度较大。然而，近年来随着水平井技术、多级压裂技术等关键技术的不断进步，致密油的开发取得了显著进展。以美国为例，其通过大规模开发致密油，实现了原油产量的大幅增长，在一定程度上改变了全球石油市场的格局。中国的致密油资源也十分丰富，主要分布在鄂尔多斯、准噶尔、松辽等盆地。据相关研究预测，中国主要盆地的致密油地质资源总量可观，可采资源量也具有较大潜力，是未来能源接替的重要领域之一。在致密油开发过程中，提高采收率是关键问题。传统的注水开发方式在致密油藏中效果往往不佳，因为致密油藏的超低渗透率使得水难以有效驱替原油，且容易造成储层伤害。注二氧化碳（CO₂）技术作为一种有效的提高采收率方法，在致密油开发中展现出了独特的优势。CO₂在一定条件下可与原油形成混相，降低原油黏度和界面张力，提高驱油效率；同时，CO₂还具有膨胀原油体积、补充地层能量等作用，能够有效改善致密油藏的开发效果。此外，注CO₂技术还与碳捕集、利用与封存（CCUS）技术紧密相关。通过将工业排放的CO₂注入致密油藏，不仅可以实现提高采收率的目的，还能将CO₂封存于地下，减少温室气体排放，对环境保护和应对气候变化具有重要意义，符合全球可持续发展的战略需求。数值模拟作为一种重要的研究手段，在致密油注CO₂开发研究中发挥着不可或缺的作用。由于致密油藏的地质条件复杂，现场试验成本高、周期长且存在诸多不确定性，数值模拟可以在计算机上建立油藏模型，模拟不同开发方案下的注CO₂过程，研究油藏的动态响应，预测开发效果。通过数值模拟，可以深入了解注CO₂的驱油机理、优化注气参数（如注气速度、注气量、注气时机等），为现场开发提供科学依据，降低开发风险，提高开发效益。因此，开展致密油注CO₂数值模拟研究，对于推动致密油资源的高效开发、实现能源的可持续供应以及促进环境保护等方面都具有重要的现实意义和理论价值。1.2国内外研究现状在致密油注二氧化碳数值模拟研究领域，国外起步相对较早，积累了丰富的研究成果。美国作为致密油开发的先驱，在该领域开展了大量深入的研究工作。斯坦福大学的研究团队针对巴肯致密油藏，通过岩心驱替实验与数值模拟相结合的方法，深入探究了二氧化碳混相驱的油回收机制。研究结果表明，在最低混相压力（MMP）以上，二氧化碳与原油混相，能够显著降低毛管力，提高局部驱替效率。通过构建数值组成模型，他们发现轻烃组分蒸发到二氧化碳中是主要的油回收机制，此外，再加压、油膨胀、粘度和界面张力降低等因素也对驱油效果有着重要影响。该研究为后续深入理解注二氧化碳驱油机理以及优化开发方案提供了重要的理论基础。在加拿大，相关学者运用数值模拟技术，对不同地质条件下的致密油藏注二氧化碳开发进行了全面的模拟分析。研究涵盖了多种因素，如储层渗透率、孔隙度、原油性质以及注气参数等对开发效果的影响。结果显示，储层渗透率和孔隙度的差异会导致二氧化碳的波及范围和驱油效率发生显著变化；原油性质，如原油的黏度、组成等，也与二氧化碳的混相效果密切相关；而注气参数，如注气速度、注气量等的合理优化，能够有效提高原油采收率。这些研究成果为加拿大致密油藏的高效开发提供了科学的指导依据。近年来，国内对致密油注二氧化碳数值模拟的研究也取得了显著进展。中国各大油田及科研院校纷纷开展相关研究工作。中国石油大学（北京）的科研团队采用数值模拟方法，对致密油藏二氧化碳吞吐的新特征展开了深入研究。他们着重分析了扩散系数、渗透率、储层非均质性及吞吐周期等因素对开发效果的影响。研究表明，在致密油藏中，二氧化碳的扩散系数不能被忽视，其大小对最终采收率有着重要影响；对于渗透率越低、储层非均质性越强的致密油藏，采用二氧化碳吞吐的增油效果更加明显，但随着吞吐轮次的增加，增油效果会逐渐减弱。这一研究成果为致密油藏优选二氧化碳吞吐增油措施提供了关键的参考依据。中国石化针对其所属的部分致密油藏，通过数值模拟技术，系统研究了注二氧化碳驱油过程中的渗流特征以及提高采收率的潜力。研究人员详细分析了二氧化碳在储层中的渗流规律，包括二氧化碳的扩散、对流等过程，以及这些过程对原油运移和驱替的影响。同时，他们还评估了不同注气方式和参数对采收率的提升效果，如连续注气、间歇注气、注气速度和注气量的优化等。研究结果表明，通过合理设计注气方案，能够有效提高二氧化碳的驱油效率，显著提升原油采收率，为中国石化致密油藏的高效开发提供了重要的技术支持。此外，国内学者还运用分子动力学模拟方法，从微观层面深入研究了二氧化碳与原油的混相机理以及二氧化碳在纳米孔隙中的驱油作用规律。通过构建分子模型，模拟二氧化碳与原油分子在纳米尺度下的相互作用，揭示了混相过程中分子间的作用力、扩散行为以及微观驱油机制。这些微观研究成果为宏观数值模拟提供了更准确的基础数据和理论支持，进一步完善了致密油注二氧化碳开发的理论体系。1.3研究目标与内容本研究旨在通过数值模拟方法，深入探究致密油注二氧化碳开发过程中的关键科学问题和技术难题，为致密油藏的高效开发提供科学依据和技术支持，具体研究目标如下：揭示注二氧化碳驱油机理：深入剖析二氧化碳在致密油藏中的渗流特性、与原油的混相机理以及对原油性质的影响，从微观和宏观层面揭示注二氧化碳提高采收率的内在机制。建立准确的数值模拟模型：综合考虑致密油藏的地质特征、流体性质以及注二氧化碳开发过程中的复杂物理化学过程，构建高精度的数值模拟模型，实现对注二氧化碳开发过程的准确模拟和预测。优化注二氧化碳开发方案：基于数值模拟结果，系统研究注气参数（如注气速度、注气量、注气时机等）、井网布置以及油藏地质条件对开发效果的影响规律，优化注二氧化碳开发方案，提高原油采收率和经济效益。评估注二氧化碳开发的环境影响：分析注二氧化碳过程中可能产生的环境问题，如二氧化碳泄漏对土壤、水体和大气环境的影响，评估注二氧化碳开发的环境风险，提出相应的环境保护措施和建议。围绕上述研究目标，本研究的具体内容包括：致密油藏地质特征与流体性质研究：收集和整理目标致密油藏的地质资料，包括地层构造、储层物性、孔隙结构、渗透率分布等；分析原油的组成、黏度、密度等性质，以及二氧化碳在原油和地层水中的溶解度、扩散系数等关键参数，为数值模拟提供准确的基础数据。注二氧化碳驱油机理实验研究：开展室内岩心驱替实验，模拟注二氧化碳驱油过程，观察和分析二氧化碳在岩心中的渗流特征、混相现象以及原油的驱替效果；利用微观实验技术，如扫描电镜、核磁共振等，研究二氧化碳与原油在微观孔隙尺度下的相互作用机制，获取相关实验数据和认识，为数值模拟模型的建立和验证提供实验依据。数值模拟模型的建立与验证：选择合适的数值模拟软件，如CMG、Eclipse等，根据目标致密油藏的地质特征和流体性质，建立三维多相多组分数值模拟模型。模型应考虑二氧化碳的溶解、扩散、混相驱替，原油的流动，以及岩石和流体的物理化学性质变化等过程。通过与室内实验数据和现场生产数据的对比，对数值模拟模型进行验证和校准，确保模型的准确性和可靠性。注二氧化碳开发方案优化研究：利用建立的数值模拟模型，系统研究不同注气参数（注气速度、注气量、注气时机）、井网布置（井距、排距、井型）以及油藏地质条件（渗透率、孔隙度、非均质性）对注二氧化碳开发效果的影响规律。采用正交试验设计、响应面分析等优化方法，对注二氧化碳开发方案进行多参数优化，确定最优的开发方案，提高原油采收率和经济效益。注二氧化碳开发的环境影响评估：基于数值模拟结果，分析注二氧化碳过程中二氧化碳的泄漏风险和可能对环境造成的影响，如对土壤、水体和大气环境的污染。建立环境影响评估模型，评估注二氧化碳开发的环境风险，提出相应的二氧化碳泄漏监测、控制和环境保护措施，确保注二氧化碳开发的环境安全性。二、致密油藏与注二氧化碳技术概述2.1致密油藏特征剖析2.1.1定义与特点致密油是一种重要的非常规石油资源，通常指赋存于致密砂岩、致密碳酸盐岩等储层中的石油。这些储层的覆压基质渗透率小于0.1mD（部分学者认为是0.2mD），具有典型的低渗透率特性。低渗透率使得原油在储层中的流动极为困难，常规的开采方式难以实现有效开采。例如，与常规砂岩储层相比，致密砂岩储层的渗透率可能低几个数量级，导致原油的渗流速度极慢，大大增加了开采难度和成本。除了低渗透率，致密油藏还具有低孔隙度的特点，孔隙度一般小于10%。较小的孔隙度意味着储层的储集空间有限，原油的储存量相对较少。而且，这些孔隙往往连通性较差，进一步阻碍了原油的流动。在微观层面，致密油藏的孔隙结构十分复杂，孔隙和喉道细小，部分孔隙甚至呈纳米级，形成了微纳米级的孔喉系统。这种微小的孔喉结构使得原油在其中的流动受到极大的限制，增加了原油开采的技术难度。此外，致密油藏的原油性质也具有独特之处。其油质相对较轻，以轻质油为主，原油密度通常大于40API度（美国石油学会制定的用以表示石油和石油产品密度的量度）或小于0.8251克/立方厘米，这使得原油具有较好的流动性。与重质原油相比，轻质原油在开采和运输过程中相对更容易处理，但其在致密储层中的渗流特性仍与常规储层有很大差异，需要特殊的开采技术和工艺来实现高效开采。2.1.2地质成因致密油的形成是一个复杂的地质过程，受到多种因素的综合影响。从沉积环境来看，致密油主要形成于特定的沉积体系中，如咸化湖泊沉积环境、深湖水下前三角洲沉积环境以及深湖凹陷或斜坡部位等。在咸化湖泊沉积环境中，优质烃源岩与碳酸盐岩富集层或膏盐层呈互层分布，咸化湖泊碳酸盐岩夹持在半深湖-深湖相暗色泥页岩中，为致密油的形成提供了良好的物质基础和储集空间。深湖水下前三角洲沉积环境则在满足优质泥质烃源岩沉积的同时，接受三角洲前缘输送的薄层粉细砂岩，形成优质烃源岩与薄层粉细砂岩互层或紧邻、源储紧密接触的致密层段，有利于致密油的聚集。构造运动在致密油的形成过程中也起着关键作用。在地质历史时期，构造运动导致地层发生褶皱、断裂等变形，这些构造形迹不仅影响了烃源岩的分布和演化，还为原油的运移和聚集提供了通道和场所。断层的存在可以沟通不同的储层和烃源岩，使得原油能够从烃源岩运移到致密储层中聚集；褶皱构造则可以形成圈闭，有利于原油的保存。此外，构造运动还会改变地层的压力和温度条件，影响原油的生成和演化过程。有机质的演化是致密油形成的核心因素之一。在适宜的地质条件下，大量的生物遗骸等有机质在沉积过程中被埋藏在地下，经过漫长的地质时期，在温度、压力等因素的作用下，发生一系列复杂的生物化学和物理化学变化，逐渐转化为石油。这一过程需要合适的温度和压力范围，以及足够的时间。如果温度过高或过低，有机质的演化进程可能会受到抑制或改变，从而影响致密油的形成。而且，不同类型的有机质对致密油的形成也有不同的影响，腐泥型有机质通常更有利于生成石油。2.1.3开发挑战致密油藏的特殊地质特征给开发带来了诸多严峻挑战。首先，采收率低是一个突出问题。由于储层的低渗透率和低孔隙度，常规的注水开发方式往往效果不佳。水在致密储层中的渗流阻力大，难以有效驱替原油，导致原油采收率较低。据相关研究和实际生产数据表明，传统注水开发在致密油藏中的采收率通常仅为10%-20%左右，远远低于常规油藏的采收率水平。这使得大量的原油残留在地下，无法得到有效开采，造成了资源的浪费。开采成本高也是致密油开发面临的重要挑战之一。为了实现致密油的有效开采，需要采用一系列先进的技术和工艺，如水平井技术、多级压裂技术等。这些技术的应用虽然能够提高原油产量，但也大大增加了开采成本。水平井的钻井成本比直井高出数倍，多级压裂需要消耗大量的材料和能源，且施工过程复杂，对设备和技术人员的要求也很高。此外，由于致密油藏的地质条件复杂，勘探难度大，需要进行大量的地质勘探工作，这也进一步增加了开发成本。储层伤害风险大是致密油开发中不可忽视的问题。在开采过程中，钻井、压裂、注水等作业环节都可能对储层造成伤害。例如，钻井液中的固相颗粒可能会堵塞储层孔隙和喉道，降低储层渗透率；压裂液如果选择不当或返排不及时，可能会对储层岩石和流体性质产生不良影响，导致储层伤害。而且，由于致密油藏的孔喉细小，一旦发生储层伤害，恢复难度极大。储层伤害不仅会降低原油产量，还会缩短油井的生产寿命，增加开发成本和风险。2.2CO₂驱油技术解析2.2.1CO₂的物理化学性质二氧化碳（CO₂）是一种碳氧化合物，在常温常压下，它呈现为无色无味的气体状态。其化学式为CO₂，化学式量约为44.0095，密度比空气大，标准条件下密度约为1.98g/L。二氧化碳具有较低的沸点和熔点，沸点为-56.6℃（527kPa），熔点为-78.5℃。固态二氧化碳俗称干冰，其外形类似冰雪，在常温常压下会直接升华，从固态转变为气态，同时吸收大量热量，可用于制冷、保鲜以及制造舞台烟雾效果等。在化学性质方面，二氧化碳化学性质相对稳定，热稳定性较高，在2000℃时仅有1.8%会发生分解。它不能燃烧，一般情况下也不支持燃烧，常被用于灭火，如二氧化碳灭火器就是利用其这一特性，在灭火时，液态二氧化碳喷出后迅速气化，覆盖在燃烧物表面，隔绝氧气，从而达到灭火目的。二氧化碳属于酸性氧化物，能与水发生反应生成碳酸，这是一个可逆反应，碳酸是一种弱酸，在水中会发生部分电离。此外，二氧化碳还能与碱溶液发生反应，生成盐和水，例如将二氧化碳通入氢氧化钠溶液中，会生成碳酸钠或碳酸氢钠；通入氢氧化钙溶液中，会使澄清石灰水变浑浊，生成碳酸钙沉淀，这一反应常被用于检验二氧化碳的存在。2.2.2CO₂驱油原理注二氧化碳驱油技术的原理基于二氧化碳与原油之间一系列复杂的物理化学作用。首先，二氧化碳具有降低原油粘度的作用。当二氧化碳溶解于原油中时，会使原油分子之间的相互作用力减弱，从而降低原油的粘度。原油初始粘度越高，二氧化碳的降粘效果越显著，可使原油粘度降低至原来的1/10左右。粘度的降低使得原油在储层中的流动能力大大增强，更易于被驱替至井底，从而提高原油产量。二氧化碳能够使原油体积膨胀。大量的二氧化碳溶于原油后，会使原油分子间距增大，导致原油体积膨胀。原油体积膨胀不仅增加了地层的弹性能量，为原油的流动提供了更多的驱动力，还有利于膨胀后的剩余油脱离地层水以及岩石表面的束缚，使原本难以流动的残余油变成可动油，进而提高驱油效率。实验数据表明，二氧化碳溶于原油后，可使原油体积膨胀10%-40%，有效提高了原油的开采率。在一定压力条件下，二氧化碳与原油还能产生混相效应。当压力达到足够高时，二氧化碳会萃取原油中的轻质烃类分子，使其进入气相，形成富含烃类的气相和溶解了二氧化碳的原油液相。随着混相过程的进行，二氧化碳不断富化，与原油之间的界面张力逐渐降低，直至趋近于零，实现混相。此时，两种流体相互溶解，不存在明显的界面，驱油效率极高，理论上可使排驱剂所到之处的原油百分之百被采出。然而，实现混相驱需要较高的压力，且原油中轻质组分C₂-C₆的含量要足够高，否则难以达到混相条件。二氧化碳还具有溶解气驱的作用。在开采过程中，随着压力下降，溶解在原油中的二氧化碳会从液体中逸出，产生气体驱动力，推动原油流动，提高驱油效果。此外，部分二氧化碳驱油后会占据一定的孔隙空间，成为束缚气，阻止原油再次附着在岩石表面，也有助于原油增产。2.2.3技术优势与应用前景与其他驱油技术相比，注二氧化碳技术具有诸多显著优势。从驱油效率方面来看，二氧化碳能够通过多种作用机制提高驱油效率，如降低原油粘度、改善流度比、实现混相驱等，相比传统的注水开发，可大幅提高原油采收率。在一些适合的油藏中，注二氧化碳驱油的采收率可比水驱提高10%-30%甚至更多。而且，注二氧化碳技术不受油藏深度、温度、压力以及地层水矿化度等条件的严格限制，具有广泛的适用性。无论是浅层、深层油藏，还是致密油藏、高渗透油藏，都有应用注二氧化碳技术的成功案例。注二氧化碳技术还具有环保优势，它与碳捕集、利用与封存（CCUS）技术紧密结合，能够将工业排放的二氧化碳注入油藏，实现二氧化碳的资源化利用和封存，减少温室气体排放，对环境保护和应对气候变化具有重要意义。在全球倡导绿色低碳发展的背景下，这一环保优势使得注二氧化碳技术更具吸引力。从应用前景来看，随着全球对能源需求的持续增长以及对环境保护要求的不断提高，注二氧化碳技术在全球能源领域展现出广阔的应用前景。在致密油开发方面，由于致密油藏的特殊地质条件，传统开采技术效果不佳，而注二氧化碳技术能够有效改善致密油藏的开发效果，提高采收率，因此在致密油藏开发中具有巨大的应用潜力。美国作为致密油开发的领先国家，已经在多个致密油产区成功应用注二氧化碳技术，取得了良好的开发效果。在其他类型油藏的开发中，注二氧化碳技术也逐渐得到推广和应用。许多国家和地区都在积极开展相关的研究和试验工作，探索适合本国或本地区油藏特点的注二氧化碳开发方案。随着技术的不断进步和成本的逐渐降低，注二氧化碳技术有望在全球范围内得到更广泛的应用，成为提高原油采收率、保障能源供应安全以及实现环境保护的重要技术手段。三、数值模拟原理与方法3.1数值模拟基本原理数值模拟是一种借助计算机技术，对复杂的物理现象或工程问题进行定量分析和预测的方法。在石油工程领域，数值模拟主要用于研究油藏中流体的流动规律，预测油藏开发动态，为油藏开发方案的制定和优化提供科学依据。其基本原理是基于渗流力学理论，将油藏中的复杂物理过程抽象为一组数学方程，并通过数值方法求解这些方程，从而实现对油藏动态的模拟。油藏中的流体渗流遵循质量守恒定律、动量守恒定律和能量守恒定律。以三维多相多组分油藏为例，其基本的渗流数学模型可由以下方程组成：质量守恒方程：描述了油藏中各相（油相、气相、水相）和各组分在空间和时间上的质量变化关系。对于油相，其质量守恒方程可表示为：\frac{\partial(\phi\rho_oS_o)}{\partialt}+\nabla\cdot(\rho_o\vec{v}_o)=q_o其中，\phi为孔隙度，\rho_o为油相密度，S_o为油相饱和度，t为时间，\vec{v}_o为油相渗流速度，q_o为油相源汇项。气相和水相也有类似的质量守恒方程。动量守恒方程：即达西定律，它描述了流体在多孔介质中的渗流速度与压力梯度、渗透率以及流体粘度之间的关系。对于油相，达西定律的表达式为：\vec{v}_o=-\frac{kk_{ro}}{\mu_o}(\nablap_o-\rho_og\nablaD)其中，k为绝对渗透率，k_{ro}为油相相对渗透率，\mu_o为油相粘度，p_o为油相压力，g为重力加速度，D为深度。气相和水相的动量守恒方程形式类似，但参数有所不同。能量守恒方程：考虑了油藏中流体和岩石的能量变化，包括内能、动能和势能等。在一些复杂的油藏模拟中，如热采过程，能量守恒方程起着关键作用。对于等温条件下的油藏，能量守恒方程可简化处理。状态方程：用于描述流体的物理性质随压力、温度等条件的变化关系。例如，对于理想气体，其状态方程为pV=nRT，在油藏模拟中，需要根据实际流体的性质，采用合适的状态方程来描述气相和油相的性质。在实际的油藏数值模拟中，由于上述偏微分方程通常是非线性的，难以直接求解，因此需要采用数值方法将其离散化。常用的数值方法有有限差分法、有限元法和有限体积法等。有限差分法是将油藏的连续空间和时间区域离散化为有限个网格和时间步长。在每个网格节点上，通过对偏微分方程进行差分离散，将其转化为代数方程组。例如，对于一阶导数\frac{\partialf}{\partialx}，可以采用向前差分、向后差分或中心差分等格式进行近似，如向前差分格式为\frac{\partialf}{\partialx}\approx\frac{f_{i+1}-f_i}{\Deltax}，其中f_i和f_{i+1}分别是x方向上相邻两个节点的函数值，\Deltax为节点间距。通过这种方式，将偏微分方程转化为代数方程组，然后通过迭代求解代数方程组，得到每个网格节点上的物理量（如压力、饱和度等）随时间的变化。有限差分法的优点是计算简单、易于编程实现，在早期的油藏数值模拟中得到了广泛应用。有限元法是将油藏区域划分为有限个单元，每个单元内的物理量通过插值函数来表示。通过对能量泛函或变分原理进行离散化，得到单元的有限元方程，再将所有单元的方程组装成整个油藏的方程组。有限元法能够较好地处理复杂的油藏几何形状和边界条件，对于一些具有复杂地质构造的油藏模拟具有优势。但有限元法的计算量较大，对计算机内存和计算速度要求较高。有限体积法是将油藏划分为一系列控制体积，在每个控制体积内应用守恒定律。通过对控制体积界面上的通量进行计算和近似，将偏微分方程转化为离散的代数方程。有限体积法的优点是保证了每个控制体积内的守恒性，在处理多相流和复杂物理过程时具有较好的精度和稳定性。目前，有限体积法在油藏数值模拟中应用也较为广泛。在致密油注二氧化碳数值模拟中，除了考虑上述基本方程和数值方法外，还需要考虑二氧化碳与原油的特殊物理化学作用，如二氧化碳的溶解、扩散、混相驱替等过程。这些过程需要通过相应的数学模型和参数来描述，并融入到数值模拟的基本框架中。例如，对于二氧化碳在原油中的溶解过程，可以采用亨利定律等相关模型来描述二氧化碳的溶解度与压力、温度等因素的关系；对于二氧化碳的扩散过程，可以采用菲克定律来描述其扩散速率。通过综合考虑这些因素，构建准确的数值模拟模型，才能实现对致密油注二氧化碳开发过程的有效模拟和分析。3.2模拟方法选择与模型构建3.2.1方法对比与选择在致密油注二氧化碳数值模拟中，常用的数值模拟方法主要有有限差分法、有限元法和有限体积法，它们各自具有独特的特点和适用场景。有限差分法是最早应用于油藏数值模拟的方法之一，它将油藏的连续空间和时间区域离散化为有限个网格和时间步长。在每个网格节点上，通过对描述油藏渗流过程的偏微分方程进行差分离散，将其转化为代数方程组。以一维单相渗流方程\frac{\partial(\phi\rhoS)}{\partialt}=\frac{\partial}{\partialx}(\frac{kk_r}{\mu}\frac{\partialp}{\partialx})为例，采用向前差分格式对时间导数进行离散，中心差分格式对空间导数进行离散。这种方法计算简单，易于编程实现，计算效率相对较高。然而，有限差分法在处理复杂油藏边界和非均匀网格时存在一定的局限性，其精度可能会受到网格划分的影响。在处理不规则形状的油藏时，有限差分法的网格划分较为困难，可能会导致计算误差增大。有限元法是一种基于变分原理的数值方法，它将油藏区域划分为有限个单元，每个单元内的物理量通过插值函数来表示。通过对能量泛函或变分原理进行离散化，得到单元的有限元方程，再将所有单元的方程组装成整个油藏的方程组。有限元法能够较好地处理复杂的油藏几何形状和边界条件，对于具有复杂地质构造的致密油藏模拟具有优势。在模拟含有断层、裂缝等复杂地质特征的油藏时，有限元法可以通过合理划分单元，准确地描述这些地质特征对渗流的影响。但有限元法的计算量较大，对计算机内存和计算速度要求较高，计算成本相对较高。而且，有限元法的编程实现相对复杂，需要较高的数学和编程基础。有限体积法是将油藏划分为一系列控制体积，在每个控制体积内应用守恒定律。通过对控制体积界面上的通量进行计算和近似，将偏微分方程转化为离散的代数方程。有限体积法的优点是保证了每个控制体积内的守恒性，在处理多相流和复杂物理过程时具有较好的精度和稳定性。在致密油注二氧化碳模拟中，涉及到二氧化碳的溶解、扩散、混相驱替等复杂过程，有限体积法能够较好地处理这些过程，准确地描述流体的流动和传质现象。此外，有限体积法在网格划分上具有较大的灵活性，可以适应不同形状和大小的控制体积。综合对比这三种方法，考虑到致密油藏注二氧化碳开发过程中涉及到复杂的多相多组分流动、二氧化碳与原油的物理化学作用以及复杂的地质条件，有限体积法在保证计算精度和稳定性的同时，能够更好地处理这些复杂过程。因此，本研究选择有限体积法作为致密油注二氧化碳数值模拟的方法。同时，结合实际研究需求，选用专业的油藏数值模拟软件CMG（ComputerModelingGroup）进行模拟计算。CMG软件具有强大的功能，能够处理多种类型的油藏和复杂的开采方式，在石油工业界得到了广泛的应用。它提供了丰富的物理模型和参数设置选项，能够准确地模拟二氧化碳在致密油藏中的渗流、混相驱替等过程，为研究提供了有力的工具支持。3.2.2模型构建流程构建致密油注二氧化碳数值模拟模型的具体流程如下：数据收集与整理：全面收集目标致密油藏的地质数据，包括地层构造数据，如地层的起伏、褶皱等形态信息，这些信息对于确定油藏的几何形状和边界条件至关重要；储层物性数据，如孔隙度、渗透率、饱和度等，它们直接影响着流体在储层中的流动特性。孔隙度决定了储层的储集空间大小，渗透率则控制着流体的渗流能力。收集原油和二氧化碳的物理化学性质数据，如原油的组成、黏度、密度，二氧化碳在原油和地层水中的溶解度、扩散系数等。原油的组成会影响其与二氧化碳的混相效果，而二氧化碳的溶解度和扩散系数则对其在油藏中的传质过程有着重要影响。地质模型建立：利用收集到的地质数据，借助地质建模软件（如Petrel等），构建目标致密油藏的三维地质模型。在建模过程中，首先根据地层构造数据建立地层的三维几何模型，确定油藏的边界范围和各层的空间分布。然后，根据储层物性数据，采用合适的插值方法（如克里金插值法等），对孔隙度、渗透率等物性参数进行空间插值，得到整个油藏的物性分布模型。通过地质统计学方法，考虑储层的非均质性，使建立的地质模型更加符合实际油藏情况。数值模型搭建：将建立好的地质模型导入到CMG数值模拟软件中，基于有限体积法搭建数值模拟模型。对油藏进行网格划分，根据油藏的大小、形状和地质特征，选择合适的网格类型和尺寸。在网格划分时，要充分考虑到计算精度和计算效率的平衡，对于地质条件复杂的区域（如裂缝发育区、渗透率变化较大的区域等），适当加密网格，以提高计算精度；对于地质条件相对简单的区域，可以适当增大网格尺寸，以减少计算量。设置模拟的初始条件，包括油藏的初始压力、温度、各相饱和度等。初始压力和温度决定了油藏的初始能量状态，各相饱和度则反映了油藏中油、气、水的初始分布情况。定义边界条件，如油藏的顶部、底部和侧面边界的压力条件、流量条件等。边界条件的设置要根据实际油藏的开采情况和地质条件进行合理确定，以准确模拟油藏与外界的物质和能量交换。物理过程描述：在数值模型中，准确描述注二氧化碳开发过程中的各种物理过程。对于二氧化碳的溶解过程，采用亨利定律等相关模型来描述二氧化碳在原油和地层水中的溶解度与压力、温度等因素的关系。亨利定律表明，在一定温度下，气体在液体中的溶解度与该气体的分压成正比。对于二氧化碳的扩散过程，采用菲克定律来描述其扩散速率，菲克定律指出，物质的扩散通量与浓度梯度成正比。考虑二氧化碳与原油的混相驱替过程，通过设置合适的混相模型（如EOS状态方程等）来描述混相条件和混相过程中流体性质的变化。EOS状态方程能够准确描述二氧化碳与原油在不同压力、温度条件下的相态变化和混相行为。3.2.3参数设置与验证模型中的关键参数包括储层参数、流体参数和注气参数等。储层参数方面，孔隙度直接影响储层的储集能力，其取值范围通常根据实际地质数据确定，一般在5%-15%之间。渗透率是控制流体渗流的关键参数，对于致密油藏，渗透率极低，通常在0.01-0.1mD之间。渗透率的方向性也需要考虑，不同方向的渗透率可能存在差异，这会影响流体在储层中的流动方向和速度。流体参数中，原油黏度对驱油效率有重要影响，其数值可通过实验测量或根据原油组成和性质估算得到。在常温常压下，致密油的黏度一般在1-10mPa・s之间。二氧化碳在原油中的溶解度和扩散系数也是关键参数，溶解度与压力、温度密切相关，压力越高、温度越低，溶解度越大。扩散系数则反映了二氧化碳在原油中的扩散能力，其数值通常通过实验测定或参考相关文献数据。注气参数中，注气速度和注气量会影响二氧化碳在油藏中的分布和驱油效果。注气速度过大可能导致二氧化碳突破过快，降低驱油效率；注气速度过小则会延长开发周期，增加成本。注气量的多少也会影响最终的采收率，需要根据油藏的大小、储层物性和原油储量等因素进行合理确定。注气时机的选择也至关重要，过早注气可能导致能量浪费，过晚注气则可能错过最佳的开发时机。参数验证的方法主要是将数值模拟结果与室内实验数据和现场生产数据进行对比分析。室内实验方面，开展岩心驱替实验，模拟注二氧化碳驱油过程，测量不同时刻的驱油效率、油气产量等数据。将数值模拟得到的相应数据与实验数据进行对比，通过调整模型参数，使模拟结果与实验数据尽可能吻合。例如，如果模拟得到的驱油效率低于实验数据，可以适当调整二氧化碳的扩散系数、混相参数等，以提高模拟结果的准确性。现场生产数据验证时，收集目标油藏的实际生产数据，包括油井产量、压力变化、含水率等。将数值模拟结果与现场生产数据进行历史拟合，通过不断调整模型参数，使模拟结果能够准确反映油藏的实际生产动态。在历史拟合过程中，需要综合考虑多种因素，如储层非均质性、井网布置、开采方式等对生产数据的影响。通过反复验证和调整参数，确保模型的准确性和可靠性，为后续的开发方案优化和生产预测提供坚实的基础。四、案例研究4.1BAKKEN油田案例4.1.1油田概况BAKKEN油田位于Williston盆地，该盆地是美国陆上最大的内克拉通沉积盆地，形状呈半圆形。BAKKEN油田横跨美国和加拿大，在美国北达科他州、蒙大拿州以及加拿大的马尼托巴省、萨斯喀彻温省均有分布。据美国地质调查局于2008年估计，BAKKEN页岩的原始地质储量达3000亿桶，可采储量为数十亿桶的原油外加一定量的溶解气。从地层特征来看，BAKKEN地层属上泥盆系到下密西西比系沉积，最大厚度约为44米（145英尺），从下到上可划分为三个段：下BAKKEN段、中BAKKEN岩性多变段、上BAKKEN页岩段。上、下BAKKEN页岩段岩性较为均一，由灰黑色、坚硬的、含硅质和钙质及黄铁矿的黑色有机质、粘土、泥粒大小的石英颗粒和部分方解石及白云石组成，有机质含量丰富，是BAKKEN页岩的烃源岩。中BAKKEN段在垂向上岩性多变，由泥质白云岩或灰岩及砂岩组成。BAKKEN油田的储层具有典型的致密油藏特征。其孔隙度较低，主力储集层段孔隙度一般在10%-13%之间。渗透率极低，覆压基质渗透率通常小于0.1mD，部分区域甚至低至0.01mD以下。这种低孔隙度和超低渗透率的储层条件，使得原油在储层中的流动极为困难。然而，该油田的油质相对较轻，以轻质油为主，原油密度通常大于40API度，具有较好的流动性。在构造方面，Williston盆地内主要构造有Nesson背斜、Billings背斜和CedarCreek等背斜。这些构造对BAKKEN油田的油气分布和聚集起到了重要的控制作用。背斜构造形成的圈闭为油气的聚集提供了有利场所，使得油气在这些构造部位相对富集。此外，该地区的断层也对油气的运移和分布产生影响，断层可以作为油气运移的通道，使油气从烃源岩运移到储层中。4.1.2注CO₂混相驱实验在BAKKEN油田的注CO₂混相驱实验中，选用了平均孔隙率为7.5%，渗透率为1.8微达西的BAKKEN岩心样本。实验前，对岩心样本进行了预处理，以确保其具有代表性和实验的准确性。将岩心清洗干净，去除表面的杂质和污染物，然后进行烘干处理，使其达到恒重状态。接着，对岩心进行抽真空饱和油处理，将其置于真空环境中，使岩心中的孔隙充满原油。实验过程中，采用高精度的驱替设备，将CO₂以一定的压力和流量注入岩心。实验压力控制在高于最低混相压力（MMP）的条件下，以确保CO₂与原油能够实现混相。通过调节注入泵的流量，精确控制CO₂的注入速度。同时，利用数据采集系统，实时监测注入压力、流量以及产出液的组成和流量等参数。在注入CO₂的初期，密切观察压力的变化情况，确保压力稳定在设定值附近。随着CO₂的注入，记录不同时刻产出液中原油和CO₂的含量，以及产出液的体积。为了更直观地观察流体在岩心中的流动和分布情况，采用了双能量CT扫描技术。在岩心驱替实验过程中，将岩心置于CT扫描设备中，进行不同时刻的扫描。通过双能量CT扫描，可以获取岩心内部不同位置的密度信息，从而直观地显示CO₂的渗透路径、原油的分布变化以及两者之间的相互作用情况。在较低能量水平下，CT扫描能够更好地捕捉到CO₂渗透的影响，此时光电吸收的X射线衰减机制变得占主导地位。实验发现，在生产早期，CT数的变化比后期更大，可能的解释包括在生产早期可能尚未完全实现混相，以及较轻的烃类组分与CO₂一起优先产出，留下更重的组分，其密度属性更密集。4.1.3数值模拟结果与分析利用CMG软件构建数值组成模型，对注CO₂混相驱实验进行模拟，以重现实验室结果。在模型构建过程中，充分考虑了岩心的孔隙度、渗透率、原油组成、CO₂的物理化学性质以及混相驱替过程中的各种物理化学变化。通过调整模型参数，如扩散系数、吸附系数、相对渗透率曲线等，使模拟结果与实验数据尽可能吻合。模拟结果表明，数值模型能够较好地重现实验中的油回收过程。在混相驱替过程中，轻烃组分蒸发到CO₂中被显示为主要的油回收机制。随着CO₂的注入，CO₂与原油接触，萃取原油中的轻质烃类分子，使其进入气相，从而实现原油的回收。此外，再加压、油膨胀、粘度和界面张力降低等因素也对油回收起到了重要作用。CO₂溶于原油后，使原油体积膨胀，增加了原油的流动性；同时，降低了原油的粘度和界面张力，减小了原油在孔隙中的流动阻力，提高了驱油效率。在分析影响因素时发现，注气速度对驱油效果有着显著影响。当注气速度过大时，CO₂在岩心中的推进速度过快，容易导致气窜现象的发生，使得CO₂过早地突破到生产端，降低了驱油效率。而注气速度过小时，虽然可以减少气窜的风险，但会延长驱油时间，增加开发成本。因此，存在一个最佳的注气速度，能够在保证驱油效率的前提下，提高开发效益。原油组成也对混相驱替效果有重要影响。原油中轻质组分C₂-C₆的含量越高，越容易与CO₂实现混相，从而提高驱油效率。相反，如果原油中重质组分含量过高，会增加混相的难度，降低驱油效果。此外，岩心的非均质性也会影响CO₂的渗流和驱油效果。岩心中渗透率的差异会导致CO₂在不同区域的流动速度和波及范围不同，从而影响整体的驱油效率。在渗透率较高的区域，CO₂容易快速通过，而在渗透率较低的区域，CO₂的渗流受阻，导致驱油效果不佳。4.2大庆油田第十采油厂案例4.2.1项目背景与目标大庆油田作为我国重要的石油生产基地，经过长期的开发，面临着诸多挑战。随着常规油藏可采储量的逐渐减少，开发难度相对较小的油藏资源日益稀缺。同时，大庆油田部分区块的采收率较低，进一步提高采收率面临着技术和经济上的双重困难。在此背景下，致密油资源的开发成为了大庆油田实现可持续发展的重要方向之一。大庆油田第十采油厂的致密油区块具有一定的特殊性。该区块的储层呈现出低孔隙度、超低渗透率的显著特征，孔隙度通常在5%-10%之间，渗透率一般小于0.1mD。这种地质条件使得常规的注水开发方式难以达到理想的开发效果。注水过程中，水在储层中的渗流阻力极大，难以有效驱替原油，导致原油采收率较低，仅能达到10%-15%左右。此外，该区块的原油性质也较为特殊，原油粘度相对较高，这进一步增加了原油在储层中的流动难度。为了有效解决这些问题，提高大庆油田第十采油厂致密油区块的采收率，开展二氧化碳驱油先导试验具有重要的现实意义。通过注入二氧化碳，利用其独特的物理化学性质，如降低原油粘度、使原油体积膨胀以及与原油混相的能力，有望显著改善原油的流动性能，提高驱油效率。同时，二氧化碳驱油技术还能在一定程度上实现二氧化碳的地质封存，具有一定的环保效益，符合当前绿色发展的理念。该先导试验的目标主要包括以下几个方面。一是通过现场试验，深入研究二氧化碳在大庆油田第十采油厂致密油储层中的驱油效果，包括驱油效率的提升幅度、原油产量的增加情况等。二是探索适合该区块的二氧化碳驱油工艺参数，如注气速度、注气量、注气时机等，为后续的大规模开发提供技术支持。三是评估二氧化碳驱油过程中的经济效益和环境影响，分析该技术在大庆油田的应用可行性和可持续性。通过实现这些目标，为大庆油田致密油资源的高效开发提供宝贵的经验和科学依据。4.2.2实验与模拟研究内容在开展二氧化碳驱油先导试验的过程中，进行了一系列丰富而深入的室内实验研究和数值模拟工作，以全面深入地了解二氧化碳驱油的机理和效果，为现场试验提供坚实可靠的技术支撑。室内实验研究涵盖多个关键方面。首先，开展了岩心渗流实验，精心挑选了具有代表性的岩心样本，其孔隙度范围在6%-8%，渗透率在0.05-0.08mD之间，以模拟大庆油田第十采油厂致密油储层的实际地质条件。在实验过程中，将二氧化碳以不同的压力和流量注入岩心，运用高精度的压力传感器和流量监测设备，实时精确地监测岩心两端的压力变化以及二氧化碳的渗流速度。通过这些数据的采集和分析，深入研究二氧化碳在致密岩心中的渗流特性，包括渗透率随二氧化碳注入量的变化规律，以及二氧化碳在不同渗透率区域的渗流差异。实验结果表明，随着二氧化碳的注入，岩心渗透率会逐渐发生变化，在一定注入量范围内，渗透率会有所增加，这主要是由于二氧化碳溶解于原油后降低了原油粘度，改善了原油的流动性，从而减小了渗流阻力。相态研究也是室内实验的重要内容。利用先进的PVT（压力-体积-温度）实验设备，对原油和二氧化碳的相态进行了详细的研究。通过精确控制实验温度和压力，模拟储层条件下二氧化碳与原油的相互作用。实验结果显示，在不同的压力和温度条件下，二氧化碳与原油的混相状态存在显著差异。在一定压力范围内，随着压力的升高，二氧化碳与原油逐渐趋近混相，此时两者之间的界面张力显著降低，驱油效率大幅提高。同时，还分析了原油组成对混相条件的影响，发现原油中轻质组分的含量越高，越容易与二氧化碳实现混相。例如，当原油中C₂-C₆轻质组分含量达到一定比例时，在相对较低的压力下就能实现混相，这为优化驱油方案提供了重要的依据。在数值模拟方面，选用了专业的油藏数值模拟软件CMG，根据大庆油田第十采油厂致密油区块的地质特征和流体性质，构建了高精度的三维多相多组分数值模拟模型。在模型构建过程中，充分考虑了二氧化碳在原油和地层水中的溶解、扩散过程，以及二氧化碳与原油的混相驱替过程。对于二氧化碳的溶解过程，采用了经过实验验证的亨利定律模型，准确描述二氧化碳在不同压力和温度下在原油和地层水中的溶解度。在扩散过程的模拟中，运用菲克定律来描述二氧化碳的扩散速率，考虑了扩散系数与温度、压力以及孔隙结构等因素的关系。对于混相驱替过程，采用了EOS（状态方程）模型，精确模拟二氧化碳与原油在混相过程中的相态变化和流体性质的改变。通过数值模拟，对不同注气参数下的二氧化碳驱油过程进行了全面的模拟分析。研究了注气速度对驱油效果的影响，发现注气速度过大会导致二氧化碳在油藏中快速突破，无法充分与原油接触并实现有效驱替，从而降低驱油效率；而注气速度过小则会延长驱油周期，增加开发成本。经过模拟优化，确定了在该区块较为合适的注气速度范围。同时，还研究了注气量对原油采收率的影响，结果表明，随着注气量的增加，原油采收率呈现先上升后趋于平缓的趋势。当注气量达到一定程度后，继续增加注气量对采收率的提升效果不再明显，这为确定合理的注气量提供了科学依据。此外，还分析了注气时机对开发效果的影响，发现选择在油藏压力下降到一定程度后进行注气，能够更好地利用二氧化碳的能量补充和驱油作用，提高采收率。4.2.3成果与意义通过对大庆油田第十采油厂致密油区块二氧化碳驱油先导试验的深入研究，取得了一系列具有重要价值的成果，这些成果对大庆油田致密油储层的开发具有多方面的指导意义和显著的实际应用价值。从采收率提升角度来看，先导试验结果表明，注二氧化碳驱油技术在大庆油田第十采油厂致密油储层中取得了显著的增油效果。在先导试验区域，通过优化注气参数，原油采收率相比传统注水开发提高了15%-20%。这意味着在相同的地质条件下，采用注二氧化碳驱油技术能够从地下开采出更多的原油，有效提高了资源利用率。以某典型井组为例，在实施注二氧化碳驱油之前，该井组的原油年产量为X吨，采收率为Y%；实施注二氧化碳驱油后，经过一段时间的生产，原油年产量提高到了X+ΔX吨，采收率提升至Y+ΔY%，增油效果十分显著。这一成果为大庆油田其他类似致密油储层的开发提供了成功的范例，证明了注二氧化碳驱油技术在提高致密油采收率方面的巨大潜力。在注气参数优化方面，通过室内实验和数值模拟的系统研究，明确了适合大庆油田第十采油厂致密油储层的注气速度、注气量和注气时机等关键参数。研究发现，在该区块的地质条件下，注气速度控制在[具体速度范围]时，能够在保证二氧化碳有效驱替原油的同时，避免气窜现象的发生，提高驱油效率；注气量应根据油藏的大小、储层物性以及原油储量等因素进行合理确定，经过模拟和实际验证，确定了最佳的注气量范围为[具体注气量范围]；注气时机选择在油藏压力下降到[具体压力值]左右时较为合适，此时注入二氧化碳能够更好地补充地层能量，提高原油的流动性。这些优化后的注气参数为大庆油田致密油储层的大规模开发提供了具体的技术指标和操作指南，有助于提高开发效益，降低开发成本。此外，研究成果还对大庆油田的可持续发展具有重要意义。随着全球对环境保护和节能减排的要求日益提高，注二氧化碳驱油技术不仅能够提高原油采收率，还能实现二氧化碳的地质封存，减少温室气体排放。在大庆油田第十采油厂的先导试验中，通过监测和分析，发现注入地下的二氧化碳大部分能够被有效封存，封存率达到了[具体封存率]。这对于大庆油田实现绿色低碳发展，减少对环境的影响具有积极的推动作用。同时，该技术的成功应用也为大庆油田在应对气候变化和能源转型方面提供了有力的技术支持，有助于提升大庆油田在能源领域的竞争力和可持续发展能力。在实际应用价值方面，大庆油田第十采油厂的二氧化碳驱油先导试验成果可以直接应用于该油田其他致密油区块的开发。通过借鉴先导试验中优化的注气参数、驱油工艺以及监测和管理经验，能够加快其他致密油区块的开发进程，提高开发效果。这不仅有助于增加大庆油田的原油产量，保障国家的能源供应安全，还能为大庆油田带来显著的经济效益。据初步估算，在大庆油田推广应用注二氧化碳驱油技术，按照当前的原油价格和开发成本，每年可为油田增加经济效益[具体金额]。同时，该技术的应用还能带动相关产业的发展，如二氧化碳捕集、运输和储存等产业，促进区域经济的繁荣。4.3河南油田ZD区块案例4.3.1区块特点河南油田ZD区块位于南阳盆地，经历了多期构造运动，形成了复杂的构造格局。该区域主要构造类型包括背斜、断块以及与区域性构造相关的复杂构造。背斜构造使得地层发生褶皱，形成穹窿状或背斜形态，为油气的聚集提供了有利的圈闭条件。断块构造则将地层切割成不同的块体，断层的存在既可以作为油气运移的通道，也可能成为油气聚集的遮挡边界。在ZD区块的东北部，存在一条正断层，该断层使得地层发生错动，在断层的上升盘和下降盘形成了不同的储层条件和油气分布特征。在储层物性方面，ZD区块的孔隙度分布范围为9.48%-34.89%，平均值为27.75%，属于中高孔隙度储层。渗透率分布范围为0.001-6.254m²，平均值为0.82m²，为中渗储层。这种较好的储层物性有利于油气的储存和渗流。然而，储层的非均质性较为明显，变异系数范围为0.22-1.74，平均值为0.6。在不同的沉积微相区域，储层物性存在较大差异。在河道、河口坝等微相区域，由于水动力较强，沉积物分选性好，孔隙度和渗透率相对较高；而在席状砂、河道间等微相区域，水动力较弱，沉积物分选性差，储层物性相对较差。在ZD区块的西南部，河道微相区域的孔隙度可达30%以上，渗透率超过1m²；而在河道间微相区域，孔隙度可能降至15%以下，渗透率小于0.5m²。从油气藏类型来看，ZD区块油气藏类型多样，包括断块油气藏、背斜油气藏、岩性油气藏以及地层油气藏等。断块油气藏主要是由于断层的切割和遮挡作用，使得油气在断块内聚集。背斜油气藏则是在背斜构造的顶部，油气因浮力作用而聚集。岩性油气藏是由于储层岩性的变化，如砂岩透镜体、尖灭带等，形成了岩性圈闭，从而聚集油气。地层油气藏则是由于地层的不整合、超覆等地质现象，导致油气在特定的地层界面附近聚集。在ZD区块的中部，存在一个典型的背斜油气藏，背斜顶部的油气富集程度较高，油层厚度较大；而在区块的边缘，由于地层的超覆和岩性变化，形成了多个地层油气藏和岩性油气藏。4.3.2CO₂吞吐参数优化模拟为了确定河南油田ZD区块不同储层的CO₂吞吐参数最优值，研究人员开展了一系列深入的数值模拟研究。在模拟过程中，运用专业的数值模拟软件，基于该区块的地质特征和流体性质，构建了高精度的三维多相多组分数值模拟模型。该模型充分考虑了CO₂在储层中的渗流、溶解、扩散以及与原油的混相驱替等复杂物理化学过程。在注入和闷井阶段，模拟结果显示，CO₂在储层中的波及范围随着时间的推移而逐渐扩大。这是因为CO₂注入后，在压力梯度的作用下，不断向周围的储层孔隙中扩散。在这个过程中，CO₂与原油充分接触，溶解于原油中，使得原油的性质发生显著变化。原油黏度明显降低，这是由于CO₂的溶解降低了原油分子之间的相互作用力，从而改善了原油的流动性。在生产阶段，随着油井的开采，CO₂随原油一起产出，进一步扩大了原油的动用范围。这是因为CO₂在原油中的溶解和扩散，使得原本难以流动的原油变得更容易被开采出来。通过对不同CO₂注入量的模拟分析发现，换油率随CO₂注入量的增加呈现出先上升后降低的趋势。当CO₂注入量较小时，随着注入量的增加，CO₂与原油的接触面积增大，溶解量增加，降黏和膨胀等作用更加明显，从而提高了原油的采收率，换油率上升。然而，当CO₂注入量超过一定值后，继续增加注入量，CO₂在储层中的分布逐渐趋于均匀，与原油的作用效果不再显著增强，反而可能导致部分CO₂过早地突破到生产井，降低了驱油效率，使得换油率下降。闷井时间对换油率也有重要影响。随着闷井时间的延长，CO₂有更充足的时间与原油接触并发生溶解、扩散等作用，从而提高原油的采收率，换油率上升。但当闷井时间过长时，CO₂在储层中的扩散达到平衡状态，继续延长闷井时间对提高采收率的作用不大，且可能会增加开发成本，导致换油率降低。CO₂注入速率与换油率呈正相关性。较高的注入速率可以使CO₂更快地进入储层，在较短时间内与更多的原油接触，增强降黏和膨胀等作用，从而提高换油率。然而，注入速率过高也可能带来一些问题，如容易导致气窜现象的发生，使得CO₂过早地突破到生产井，降低驱油效率。因此，在实际应用中，需要在保证驱油效果的前提下，选择合适的注入速率。吞吐周期与换油率呈负相关性。随着吞吐周期的增加，储层中的原油不断被采出，剩余原油的性质和分布发生变化，CO₂与原油的作用效果逐渐减弱，导致换油率降低。在前期的吞吐周期中，原油的采收率较高，但随着吞吐次数的增多，采收率的增长幅度逐渐减小。此外，研究还发现，储层物性对CO₂吞吐参数的最优值有显著影响。储层物性越好，孔隙度和渗透率越高，CO₂在储层中的渗流和扩散速度越快，与原油的接触和作用更加充分。因此，对于物性较好的储层，最佳CO₂注入量及注入速度可以相对较低，因为较低的注入量和速度就能够保证CO₂在储层中充分发挥作用。同时，最佳闷井时间也可以相对较短，因为CO₂与原油的作用能够在较短时间内达到较好的效果。而且，物性较好的储层可以承受更多的吞吐周期，因为其储集和渗流能力较强，能够在多次吞吐过程中保持相对稳定的开采效果。4.3.3现场试验效果在河南油田ZD区块的四类储层中，选取ZA4121井开展了CO₂吞吐试验，以验证数值模拟结果的准确性，并评估CO₂吞吐技术在该区块的实际应用效果。在试验过程中，严格按照数值模拟优化得到的CO₂吞吐参数进行操作。经过一段时间的生产，ZA4121井取得了显著的增油效果。累计增油量达到303.4t，这表明CO₂吞吐技术能够有效地提高该井的原油产量，增加了原油的采收率。换油率为0.16t/t，即每注入1tCO₂，可采出0.16t原油。这一换油率在实际生产中具有较好的经济效益，证明了CO₂吞吐技术在该区块应用的可行性和有效性。从经济效益方面分析，CO₂吞吐技术的应用为该区块带来了可观的收益。以当前的原油价格和CO₂注入成本计算，ZA4121井在实施CO₂吞吐试验后，增加的原油产量带来的经济收入明显高于CO₂注入及相关操作的成本。假设原油价格为[具体价格]元/吨，CO₂注入成本及其他相关操作成本为[具体成本]元/吨，ZA4121井实施CO₂吞吐试验后，累计增油303.4t，带来的经济收入为303.4×[具体价格]元，而CO₂注入及相关操作成本为[具体成本]×CO₂注入量。经过计算，经济收入减去成本后，实现了盈利，且盈利数额较为可观。这表明CO₂吞吐技术在河南油田ZD区块具有良好的经济效益，能够为油田带来实际的经济利益。此外，CO₂吞吐技术在该区块的应用还具有一定的环保效益。通过将CO₂注入地下，实现了部分CO₂的地质封存，减少了CO₂向大气中的排放，对缓解温室效应具有积极作用。在当前全球关注气候变化和环境保护的背景下，这一环保效益进一步凸显了CO₂吞吐技术在该区块应用的重要意义和价值。五、模拟结果讨论5.1驱油效果影响因素分析5.1.1地质因素地质因素对致密油注二氧化碳驱油效果有着至关重要的影响。渗透率作为储层的关键物性参数之一，在驱油过程中起着核心作用。当渗透率较高时，二氧化碳在储层中的渗流阻力较小，能够更顺畅地在储层孔隙中流动。这使得二氧化碳能够更快速地与原油接触，充分发挥其降低原油粘度、使原油体积膨胀以及混相驱替等作用。在渗透率为0.1mD的储层中，二氧化碳注入后，能够在较短时间内扩散到较大范围，与原油充分混合，降低原油粘度，提高原油的流动性，从而有效提高驱油效率。相反，在渗透率极低的储层中，如渗透率为0.01mD时，二氧化碳的渗流受到极大阻碍，难以在储层中均匀分布。部分区域的原油无法与二氧化碳充分接触，导致驱油效果不佳，采收率明显降低。孔隙度同样是影响驱油效果的重要因素。孔隙度决定了储层的储集空间大小，较高的孔隙度意味着储层能够容纳更多的二氧化碳和原油。当孔隙度较大时，二氧化碳注入后，有更充足的空间进行扩散和与原油的相互作用。孔隙度为15%的储层相比孔隙度为5%的储层，能够容纳更多的二氧化碳，使得二氧化碳与原油的接触面积更大，混相效果更好，从而提高驱油效率。而且，较大的孔隙度还能为原油的流动提供更广阔的通道，有利于原油被驱替至井底。然而，孔隙度并非越高越好，过高的孔隙度可能会导致储层的岩石骨架稳定性下降，影响储层的长期开发效果。储层的非均质性也是不可忽视的地质因素。非均质性使得储层中不同区域的渗透率、孔隙度等物性参数存在差异。在渗透率较高的区域，二氧化碳容易快速通过，形成优势通道，导致气窜现象的发生。这使得二氧化碳过早地突破到生产井，无法充分驱替其他区域的原油，降低了驱油效率。而在渗透率较低的区域，二氧化碳的渗流受阻，原油难以被有效驱替。在一个具有强非均质性的储层中，可能会出现部分区域二氧化碳大量涌入，而部分区域二氧化碳难以到达的情况，严重影响整体驱油效果。为了应对储层非均质性的影响，可以通过优化注气方案，如采用分段注气、调整注气速度等方式，来提高二氧化碳在储层中的波及范围和驱油效率。5.1.2操作参数操作参数的合理选择对于致密油注二氧化碳驱油效果的优化起着关键作用。二氧化碳注入量是影响驱油效果的重要操作参数之一。当注入量较小时，二氧化碳在储层中的分布范围有限，无法充分与原油接触并发生作用。随着注入量的增加，二氧化碳能够扩散到更大范围的储层孔隙中，与更多的原油接触，从而增强其降粘、膨胀和混相驱替等作用。在一定范围内，增加注入量可以显著提高原油采收率。然而，当注入量超过一定程度后，继续增加注入量对采收率的提升效果逐渐减弱。这是因为过多的二氧化碳可能会导致储层压力过高，使得二氧化碳在储层中的分布趋于均匀，与原油的作用达到饱和状态。此时，部分二氧化碳可能会直接从生产井产出，而没有充分发挥驱油作用，造成资源浪费。注入速度对驱油效果也有着显著影响。较高的注入速度可以使二氧化碳快速进入储层，在较短时间内与原油接触并发生作用。在初期阶段，快速注入二氧化碳能够迅速建立压力梯度，推动原油流动，提高原油产量。然而，注入速度过快也会带来一些问题。注入速度过快可能会导致二氧化碳在储层中形成优势通道，引发气窜现象。气窜会使得二氧化碳过早地突破到生产井，降低了二氧化碳的利用率和驱油效率。相反，注入速度过慢虽然可以减少气窜的风险，但会延长开发周期，增加开发成本。而且，注入速度过慢可能导致二氧化碳在储层中的扩散不均匀，部分区域的原油无法及时与二氧化碳接触，影响驱油效果。因此，需要根据储层的具体情况，选择合适的注入速度，以平衡驱油效率和开发成本。注气时机的选择同样至关重要。在油藏开发的不同阶段，油藏的压力、原油饱和度等条件会发生变化，这会影响二氧化碳的驱油效果。在油藏压力较高的初期阶段，注入二氧化碳可能无法充分发挥其补充地层能量和驱油的作用。因为此时油藏自身的能量相对充足，二氧化碳的注入可能会被油藏压力所抑制，无法有效扩散和与原油混相。而在油藏压力下降到一定程度后，注入二氧化碳能够更好地补充地层能量，与原油充分接触并发生作用，提高驱油效率。此外，注气时机还与原油的开采阶段有关。在原油开采的后期，剩余原油的流动性较差，此时注入二氧化碳可以通过降低原油粘度、膨胀原油体积等作用，使剩余原油变得可动，提高采收率。5.2CO₂驱油机理深入探讨5.2.1混相与非混相驱油机制混相驱和非混相驱是注二氧化碳驱油过程中的两种重要驱油机制，它们在驱油原理、驱油效果等方面存在显著差异。混相驱是指在一定的压力和温度条件下，注入的二氧化碳与原油之间形成一种单相体系，两者之间不存在明显的界面。当压力达到最低混相压力（MMP）以上时，二氧化碳与原油之间会发生多级接触混相过程。在这个过程中，二氧化碳不断萃取原油中的轻质烃类分子，使其进入气相，形成富含烃类的气相和溶解了二氧化碳的原油液相。随着混相过程的进行，二氧化碳不断富化，与原油之间的界面张力逐渐降低，直至趋近于零，实现混相。此时，驱油效率极高，理论上可使排驱剂所到之处的原油百分之百被采出。在BAKKEN油田的注CO₂混相驱实验中，通过双能量CT扫描技术观察到，在混相条件下，CO₂能够均匀地分布在原油中，与原油充分混合，有效提高了原油的采收率。非混相驱则是指二氧化碳与原油不能形成单相体系，两者之间存在明显的界面。在非混相驱过程中，二氧化碳主要通过降低原油粘度、使原油体积膨胀以及溶解气驱等作用来提高驱油效率。二氧化碳溶解于原油中，使原油分子之间的相互作用力减弱，从而降低原油的粘度，提高原油的流动性。二氧化碳还能使原油体积膨胀，增加地层的弹性能量，推动原油流动。在一些渗透率较低的油藏中，虽然无法达到混相条件，但通过非混相驱，二氧化碳仍然能够在一定程度上改善原油的流动性能，提高采收率。混相驱和非混相驱的驱油效果受到多种因素的影响，其中压力是一个关键因素。在混相驱中，只有当压力达到或超过最低混相压力时，才能实现混相，从而获得较高的驱油效率。而在非混相驱中，压力的变化主要影响二氧化碳在原油中的溶解度和溶解气驱的效果。当压力升高时，二氧化碳在原油中的溶解度增加，溶解气驱的作用增强，有利于提高驱油效率；但压力过高也可能导致二氧化碳在储层中的流动阻力增大，影响驱油效果。原油组成对混相驱和非混相驱的效果也有重要影响。原油中轻质组分C₂-C₆的含量越高，越容易与二氧化碳实现混相，混相驱的效果越好；而对于非混相驱，原油中重质组分含量过高可能会增加原油的粘度，降低二氧化碳的驱油效果。5.2.2微观作用机制从微观角度来看，二氧化碳对原油物性和岩石润湿性等有着重要影响。在对原油物性的影响方面，二氧化碳与原油的相互作用主要体现在降低原油粘度和使原油体积膨胀上。二氧化碳溶解于原油中，会破坏原油分子之间的相互作用力，尤其是范德华力。原油分子间的这种作用力被削弱后，分子的活动性增强，原油的流动性提高，从而导致原油粘度降低。研究表明，当二氧化碳在原油中的溶解度达到一定程度时，原油粘度可降低至原来的1/10甚至更低。而且，大量二氧化碳溶解于原油后，会使原油分子间距增大，导致原油体积膨胀。这不仅增加了地层的弹性能量，为原油的流动提供了更多的驱动力，还有利于膨胀后的剩余油脱离地层水以及岩石表面的束缚，使原本难以流动的残余油变成可动油。二氧化碳对岩石润湿性也会产生显著影响。岩石润湿性是指岩石表面对不同流体的亲和程度，它对原油在储层中的流动和分布有着重要影响。当二氧化碳注入油藏后，会改变岩石表面的化学性质和物理结构。二氧化碳与岩石表面的矿物质发生化学反应，可能会改变岩石表面的电荷分布和官能团组成。二氧化碳溶解于原油中，改变了原油与岩石表面之间的界面性质。这些变化会导致岩石表面对原油和水的润湿性发生改变，使岩石表面更倾向于被水润湿，即亲水性增强。岩石润湿性的这种改变，有利于水在岩石表面形成连续的水膜，降低原油与岩石表面的附着力，从而提高原油的流动性。在亲水性增强的岩石表面，原油更容易被驱替，从而提高了驱油效率。此外，二氧化碳在纳米孔隙中的驱油作用也备受关注。在致密油藏中，存在大量的纳米孔隙，这些纳米孔隙的尺寸与二氧化碳分子和原油分子的大小相当。在纳米孔隙中，二氧化碳的扩散和吸附行为与宏观尺度下有所不同。由于纳米孔隙的表面效应和尺寸效应，二氧化碳分子在纳米孔隙中的扩散速率可能会加快，且更容易吸附在孔隙表面。这使得二氧化碳能够更快速地与原油分子接触并发生作用，提高了驱油效率。而且，纳米孔隙中的原油由于受到孔隙壁的限制，其流动性较差。二氧化碳的注入可以通过降低原油粘度、改变润湿性等作用，克服孔隙壁的限制，使原油在纳米孔隙中能够更有效地流动，从而提高了原油在纳米孔隙中的采收率。5.3数值模拟的准确性与局限性数值模拟作为研究致密油注二氧化碳开发的重要手段，其模拟结果与实际情况的符合程度在一定程度上验证了该方法的可靠性。通过对多个案例的研究，将数值模拟结果与室内实验数据和现场生产数据进行对比分析，发现数值模拟在一定条件下能够较好地反映注二氧化碳驱油过程中的关键现象和趋势。在BAKKEN油田注CO₂混相驱实验的数值模拟中，模拟结果能够较为准确地重现实验中的油回收过程，轻烃组分蒸发到CO₂中作为主要油回收机制的模拟结果与实验观察一致。在大庆油田第十采油厂的二氧化碳驱油先导试验中，数值模拟所预测的采收率提升趋势与实际生产数据相吻合，验证了模拟模型在该案例中的准确性。然而，数值模拟也存在一定的局限性。在模型简化方面，为了便于计算和求解，数值模拟模型通常需要对复杂的油藏地质条件和物理化学过程进行简化。实际油藏的地质构造极为复杂，存在大量的微裂缝、断层以及非均质性等特征，而在数值模拟中，往往难以完全精确地描述这些细节。对于微裂缝的模拟，虽然可以采用一些简化的模型来考虑其对渗流的影响，但与实际情况相比，仍存在一定的差距。在处理储层非均质性时，通常采用统计平均的方法来描述渗透率和孔隙度的分布，这可能会忽略一些局部的细微变化，导致模拟结果与实际情况存在偏差。参数不确定性也是影响数值模拟准确性的重要因素。油藏中的许多参数，如渗透率、孔隙度、原油粘度、二氧化碳扩散系数等，难以通过实验或测量得到精确值。这些参数的不确定性会在数值模拟过程中逐渐积累，导致模拟结果的误差增大。在不同的实验条件下，测量得到的二氧化碳在原油中的扩散系数可能存在较大差异，这使得在数值模拟中选择合适的扩散系数值变得困难。如果选择的扩散系数不准确，可能会导致模拟的二氧化碳在储层中的扩散速度和分布情况与实际不符，进而影响对驱油效果的预测。此外，数值模拟还难以完全考虑一些复杂的物理化学过程。在注二氧化碳驱油过程中，可能会发生化学反应，如二氧化碳与地层水反应生成碳酸，碳酸对储层岩石的溶蚀作用等。这些化学反应会改变储层的物性和流体的性质，进而影响驱油效果。然而，目前的数值模拟模型在描述这些化学反应过程时还存在一定的局限性，难以准确地预测其对油藏动态的影响。而且，数值模拟也无法完全考虑现场操作过程中的一些实际因素，如井筒附近的流动阻力、注气设备的性能波动等，这些因素也可能对实际的驱油效果产生影响，但在数值模拟中往往被忽略。六、结论与展望6.1研究成果总结本研究围绕致密油注二氧化碳数值模拟展开，通过对致密油藏地质特征、注二氧化碳驱油技术原理、数值模拟方法以及多个实际案例的深入研究，取得了一系列具有重要理论和实践意义的成果。在致密油藏地质特征与注二氧化碳技术方面，全面剖析了致密油藏的定义、特点、地质成因以及开发挑战。明确了致密油藏具有低孔隙度、超低渗透率、孔隙结构复杂以及油质较轻等特征，其形成受到沉积环境、构造运动和有机质演化等多种因素的综合控制。同时，深入解析了二氧化碳的物理化学性质和驱油原理，阐述了注二氧化碳技术在提高采收率和环保方面的优势以及广阔的应用前景。二氧化碳通过降低原油粘度、使原油体积膨胀、实现混相驱替以及溶解气驱等多种作用机制，能够有效提高致密油藏的驱油效率。而且，注二氧化碳技术与CCUS技术相结合，有助于减少温室气体排放，实现能源开发与环境保护的双赢。在数值模拟原理与方法研究中，深入阐述了数值模拟的基本原理，包括基于渗流力学理论的质量守恒方程、动量守恒方程、能量守恒方程和状态方程等。详细对比了有限差分法、有限元法和有限体积法等常用数值模拟方法的特点和适用场景，综合考虑致密油藏注二氧化碳开发过程的复杂性，选择了有限体积法，并利用CMG软件构建了数值模拟模型。明确了模型构建流程，包括数据收集与整理、地质模型建立、数值模型搭建以及物理过程描述等环节。同时，对模型中的关键参数进行了设置和验证，通过与室内实验数据和现场生产数据的对比分析，确保了模型的准确性和可靠性。通过对BAKKEN油田、大庆油田第十采油厂和河南油田ZD区块三个案例的研究，取得了丰富的成果。在BAKKEN油田注CO₂混相驱实验及数值模拟中，通过岩心驱替实验和双能量CT扫描技术，直观地观察到了二氧化碳在岩心中的渗流和驱油过程。数值模拟结果准确地重现了实验中的油回收过程，明确了轻烃组分蒸发到CO₂中是主要的油回收机制，同时揭示了注气速度、原油组成和岩心非均质性等因素对驱油效果的显著影响。在大庆油田第十采油厂二氧化碳驱油先导试验中，通过室内实验和数值模拟，深入研究了二氧化碳在该区块致密油储层中