单个裂隙中油运移实验及特征分析

利用狭窄平行板裂隙模型来模拟和分析单个裂隙内油运移的过程。实验中采用表面起伏的玻璃板制作平行板裂隙模型,通过改变水中蔗糖的浓度控制两相流体间的密度差,来分析不同密度差和注入速率条件下油的运移过程和特征。观察发现,油在饱含水的狭窄裂隙空间内的运移特征与在孔隙介质中的运移有相似之处:运移的路径只占通道的一部分,都具有不规则的分形特征,也可形成活塞式、指进式和路径式等模式。而这些模式的形成也同样受各种因素的影响。利用孔隙介质中油运移实验结果总结出来的、以Bo数和Ca数为纵横坐标运移相图完全可以用来表征单个裂隙中油的运移模式。     

三维交叉裂隙渗流传质特性数值模拟

对交叉裂隙渗流传质特性的定量描述是研究整个裂隙网络渗透传质特性的基础。为真实模拟水流及溶质在三维交叉裂隙中的运移过程,首先通过三维轮廓仪获取天然岩石裂隙表面的形貌数据,再应用三维重构技术生成相应的三维交叉裂隙模型,随后求解Navier-Stokes方程,假定溶质运移满足Fick定律,模拟水流和溶质在三维交叉裂隙中的运移过程。通过对比粗糙裂隙模型与平行平板模型的模拟结果发现：粗糙度对流体的分布及流动状态存在显著的影响;不同进、出口工况下的流体流动及溶质运移状态亦表明：裂隙交叉的几何形貌会显著地影响溶质混合行为。这些结果表明,目前被广泛采用的平行平板模型在评估岩体内特别是交叉口的物质运移特性时将导致较大的偏差,在将来的研究中有必要针对裂隙交叉口的几何特征建立修正的模型以提高评估的准确性。     

松辽盆地增强型地热系统(EGS)地热能开发热-水动力耦合过程

增强型地热系统地热能开发涉及到热和水动力的耦合,对应的温度和压力场时空变化特征是评价地热开发效果的关键问题。基于松辽盆地徐家围子深部地质条件,采用TOUGH2进行了地热能开发过程中裂隙-孔隙介质系统中温度和压力变化的数值模拟,分析了不同埋深水平情况下地热能开发的差别,研究了孔隙基质和裂隙介质的渗透率和孔隙度、岩石导热系数、井径、注入压力、注入温度及裂隙周围基质因素对地热能开发的影响。结果表明：采用定压力开发时生产井抽出控制整个区域的压力分布,压力梯度在注入井区域较大,并随着开发的进行,注入井的注入对压力的影响逐渐增大;温度由注入井到生产井逐渐增大,并随着开发的进行温度降低范围逐渐向生产井扩大;质量和热提取速率随时间逐渐减小。不同埋深位置的模拟结果显示,埋深大的温度相对较高,水的流动性较强,质量和热提取速率较高,压力和温度变化幅度均较大。裂隙系统的渗透率、注入井/生产井压力和注入温度、井径对深部地热开采过程中的压力和温度影响较大,从而影响热的提取效率;而孔隙基质的渗透率和孔隙度、裂隙介质的裂隙度和岩石的热传导系数的影响并不明显。     

油气幕式成藏及其驱动机制和识别标志

连续稳态流动和周期性(幕式)瞬态流动是沉积盆地流体的 2种流动方式。幕式流体流动是压力和应力的作用引起地层周期性破裂或断裂、先存裂隙周期性开启的结果。幕式成藏是沉积盆地中油气与地层水组成的混相、不连续流体的多期次充注/聚集过程,超压顶界面附近、底辟和深断裂附近是幕式成藏的有利场所。与油气稳态连续聚集过程相比,幕式成藏更快,大 中型油气田可在较短的时间内形成,根据传统模式难以成藏的年轻圈闭可成为有效的勘探目标。幕式流体流动的主要特征是流体成份和流动过程的不连续性、流体流动过程中温度、压力的快速变化及流体流动的多期性和周期性。流体的时空非均质性、流体流动的瞬态温度响应和运移相态分异及其揭示的多期流体相互作用是幕式成藏的有效识别标志。     

非均质裂隙介质中重非水相流体运移

由于裂隙介质具有强烈非均质性,使得重非水相流体(DNAPLs)在裂隙介质中的运移行为较孔隙介质更为复杂.基于指示模拟算法模拟裂隙介质的二元变量(粗糙面接触和裂隙开口),综合模拟退火算法构建裂隙渗透率随机场.采用T2VOC模拟DNAPLs在裂隙介质的运移,探讨粗糙面接触的相关长度、各向异性比、倾角以及非均质性程度对DNAPLs运移分布的影响.数值分析结果表明,粗糙面接触的相关性越差,DNAPLs污染范围越大;粗糙面接触的各向异性比和倾角增大,将导致更多DNAPLs残留;而相关长度和各向异性比增大以及倾角减小,都会引起DNAPLs锋面运移速率增大;随着渗透率非均质性增强,会导致局部粗糙面接触上蓄积的DNAPLs饱和度增大,运移速率减小.     

多孔介质中非水相流体运移的数值模拟

针对多孔介质中水、气和非水相流体(NAPLs)的多相流动特点,建立了非水相流体(NAPLs)污染物迁移模型,分析了非水相流体在土壤非饱和区和地下水系统中的运移规律。通过有限元数值解对轻非水相流体和重非水相流体在土壤系统中的迁移过程进行模拟,得到了污染物的时空分布特征和污染范围。计算结果表明,数值模拟方法能够合理地描述非水相流体的运移过程和污染特征。土体渗透性和污染物残余饱和度是其重要影响因素。     

马岭油田北三区污水回注的实验研究

根据马岭油田北三区油田水性质，结合污水与地层水混合的静态实验结果分析，发现污水回注对地层可能引起不同程度的伤害。同时，模拟地层条件做岩心流动实验。将采出的污水回注到地层与地层水混合后，亦会使地层渗透率下降，且随着污水含量的增加，渗透率伤害程度增大。     

基于微米焦点CT技术的煤岩数字岩石物理分析方法——以沁水盆地伯方3号煤为例

数字岩石物理技术既可实现煤体孔隙的三维空间表征及重构,又可模拟孔裂隙系统中流体的流动和绝对渗透率的计算。以伯方矿3号煤层为研究对象,基于微米焦点CT扫描与图像处理技术相结合,建立了真实的三维数字岩心模型;应用MATLAB编程及AVIZO软件内含的多种形态学算法,进行了数字岩心孔隙结构量化和表征,建立了等价孔隙网络模型;将AVIZO与COMSOL完美对接,实现了孔隙尺度的渗流模拟、绝对渗透率计算,并探讨了流体运移过程中的压力场及速度场的分布。研究方法丰富了现有数字岩石物理研究手段,为微观尺度上的煤体孔隙结构及流体运移研究提供了一种新的途径。      

海底多相流作用下的含水合物沉积地层静水力学性质分析

海底多相流动区域沉积物孔隙内流体迁移-甲烷输运-水合物形成是一种普遍模式,形成的水合物在孔隙内沉淀并与多孔介质骨架胶结从而改变当地的地层结构和性质。水合物的不断形成将减少沉积地层孔隙度,改变孔隙内各相间界面张力,增大当地孔隙的进入压力及毛细压力,增强地层滞后效应,降低地层渗透率,同时多相流体流动前缘气液分离带变厚而使得气柱变长。建立了在这类环境里水合物-水-气-盐共同作用下的水合物成藏模型,选择合适的参数分析了水合物形成对沉积地层静水力学性质等的影响关系。最后根据资料估算了南海北部神狐海域沉积物内甲烷气柱的分布,结果表明：随着水合物在沉积物孔隙内逐渐饱和,临界甲烷气柱长度将在接近海底面处达到最大,约为09 m。     

山东淄博市大武水源地裂隙岩溶水中污染物运移的数值研究

在分析研究淄博市大武水源地裂隙岩溶含水层的水力性质和污染物运移特征的基础上 ,对裂隙岩溶水的水头和污染物运移进行数值研究。目前国内外对裂隙岩溶水进行数值计算时 ,通常用等价多孔介质模型 ,但裂隙岩溶介质和多孔介质有很大不同。裂隙岩溶介质的储水和导水空间为裂隙网络 ,导水系数大 ,地下水的实际平均流速比孔隙水大得多 ,但给水度和贮水系数小。当用等价多孔介质模型进行模拟时应考虑这些特点。对于污染物运移的模拟 ,要同时求解水头方程和对流弥散方程 ,可采用MODFLOW和MT3D软件进行模拟。研究区裂隙岩溶水水头的数值计算表明 ,等效多孔介质模型水头的拟合误差能满足国标GB/T144 97- 93的要求。各时段地下水水量均衡计算的精度也满足要求。对流弥散方程的数值计算 ,由于Peclet数高达 95 .6 7,对流占绝对优势 ,可能存在数值弥散和数值振荡 ,因而采用多种方法进行了比较。对于同一问题 ,同时采用上游有限差分法 (UFDM) ,混合的欧拉拉格朗日方法 (特征线法MOC、改进特征线法MMOC和混合特征线法HMOC) ,总变异消减法(TVD)进行计算 ,并比较其结果。结果表明 ,混合特征线法 (HMOC)和总变异消减法 (TVD)比较适合于对流占优势的运移问题计算。由于渗透系数K和有效孔隙度θ对溶质运移结果的影响很大 ,?