缝洞型油藏三维离散缝洞数值试井模型

缝洞型碳酸盐岩油藏发育着大尺度的溶洞和裂缝,非均质性极强,缝洞型碳酸盐岩油藏问题的研究成为了世界级难题之一.由于大尺度溶洞和裂缝对储层的流体流动起主导作用,因此,基于连续介质理论的双重介质或三重介质模型已不适合其中流体流动的描述.根据大型缝洞分布地质特征,探索性地提出了一种板块组合的复合架构离散缝洞模型描述该类油藏中的流体流动,将三维空间大裂缝用板块描述,溶洞用高渗透率和高孔隙度不规则多面体团块描述.将裂缝面用二维三角形单元离散, 溶洞和基质用三维四面体离散, 利用三维混合单元有限元法对建立的不定常渗流模型进行求解,得到了三维渗流条件下的试井理论曲线及压力场分布.通过对试井理论曲线特征的分析, 获得了各敏感参数对试井曲线的影响规律.通过对1口井的实际测试资料解释结果的分析,并与实际地震反射资料及生产实际资料的对比,发现本文所建立的模型可以较好地反映裂缝和溶洞的地质动态状况,并与实际生产状况具有较好的一致性.这一结果说明了所建模型的正确性以及测试资料分析结果的可靠性.      

缝洞型多孔介质中多相流的有限体积法数值模拟

本文研究的碳酸盐岩油藏储集体属于缝洞型多孔介质.这类缝洞型多孔介质由裂缝、溶蚀孔洞和低孔隙度低渗透率的基岩组成.裂缝是空隙流体流动的主要通道;溶蚀孔洞大小从几厘米到数米不等,渗透率和孔隙度都很高,是流体主要的储集空间.由于缝洞型多孔介质空隙空间的复杂性和强非均质性,数值计算中基本控制方程的空间离散应采用非结构化网格的计算模型.本文采用有限体积法模拟缝洞型多孔介质中多相流体的流动,并给出了相应的单元中心格式有限体积法的计算公式.裂缝介质和溶洞介质中单元间多相流体的流动考虑为高速非达西流,其质量通量采用Forchheimer定律计算.非线性方程的离散选取全隐式格式,并采用Newton-Raphson迭代进行求解.通过两个二维模型注水驱油的数值模拟,验证了本文方法的有效性.     

大埋深粘土三轴高压卸载变形与强度特征

以取自山东巨野矿区埋深近 4 30～ 6 0 0m的深部粘土样为例 ,分析测试了土的常规物理力学性质参数。考虑该埋深土的天然应力水平及其在煤矿建井开挖中的卸载应力变化 ,在三轴伺服仪上 ,进行三向等压分别达 11MPa、12MPa、14MPa下的卸除围压、保持轴压的三轴高压卸载试验 ,分析了深部粘土三轴高压卸载变形和强度特点 ,得到了相关力学参数。成果丰富了土的高压卸载力学性研究 ,并为巨野矿区煤矿立井建设提供了工程地质基础参数     

基于离散缝洞网络模型的缝洞型油藏混合多尺度有限元数值模拟

针对缝洞型油藏具有多尺度特征,基于离散缝洞网络模型,建立了Darcy/Stokes-Brinkman多尺度耦合数学模型,采用多尺度混合有限元方法,对缝洞型油藏的流体流动问题进行了研究。阐述了多尺度混合有限元方法的基本原理,推导得到了Darcy/Stokes-Brinkman多尺度模型的多尺度混合有限元计算格式。数值计算结果表明,在大尺度模型上进行计算能够捕捉到小尺度上的流动特征;相对于传统有限元,多尺度混合有限元能够捕捉小尺度上的非均质特征而具有更高的计算精度,在保证计算精度的同时能够减少计算量。     

地下核爆炸引起不稳定性岩体的探测

本文通过我国某次地下核试验场地爆前及爆后钻孔中声波测试结果,计算钻孔中岩体的弹性波参数及岩体力学参数,确定因爆炸造成围岩的破裂松散区范围及地表剥离的深度,为试验安全埋深提供了地质依据。     

分形在油气田开发中的应用

分形用于测井资料分析，为油藏描述提供了新的工具，使得利用油藏数值模拟来模拟复杂结构的油气藏变得更加便利；分形油藏上的试井分析解决了一些以前难以解释的问题；用分形来描述裂缝和孔洞结构，提高了酸化压裂的设计水平；用分形理论来描述两相流，揭示了部分粘性指进的触发机理．总之，分形理论在油气田开发中的应用，促进了渗流力学基础理论的发展．     

石灰岩中不同埋深爆炸自由场直接地冲击参数的预计方法

通过现场实验,研究了石灰岩中常规装药不同埋深爆炸自由场直接地冲击参数的传播规律.采用了等效当量埋深系数的工程处理方法来描述不同爆炸埋深地冲击效应,并给出了石灰岩中不同埋深爆炸应力等效当量埋深系数和质点加速度等效当量埋深系数的计算公式,从而得到了石灰岩介质中不同埋深爆炸地冲击参数随比例距离的预计公式.     

裂缝性低渗透油藏流固耦合渗流分析

在低渗透油田的开发过程中，油藏流体渗流和储层岩土之间存在明显的耦合作用。本文首先研究给出了低渗裂缝性储层孔渗参数的等效方法，然后将渗流力学和岩土力学相结合，给出了低渗透裂缝性储层流固耦合渗流的数学模型，该模型不仅可以反映基质孔渗参数在开发中的变化，而且更能反映裂缝开度变化所引起的渗透率变化，而这对于低渗透裂缝性油田而言十分重要。最后对一实际井网进行了流固耦合油藏数值模拟，给出了开发过程中孔渗参数的变化及其耦合效应对油田开发的影响．     

金沙江中游凝灰岩强度参数综合识别研究

凝灰岩是一种典型的分布不规律、易崩解、难取样火山碎屑岩,其强度参数难以准确评估。本文依托华丽高速公路金沙江特大桥工程,判定丽江岸为反倾坡,为陡倾结构面或卸荷裂隙带剪入、缓倾结构面剪出组合滑移模式,破坏规模有限;华坪岸为顺层坡,为陡倾结构面或卸荷裂隙带剪入、凝灰岩层主滑、缓倾结构面剪出组合滑移模式,破坏规模较大。T3和T2凝灰岩软弱层为关键层,其强度参数对岸坡稳定性影响很大。针对此种情况提出了经验、试验、数值反分析相结合综合识别凝灰岩强度参数的方法。利用这种方法,在工程岩体分级标准GBT50218-2014、水利水电岩石力学参数手册、室内外试验、数值反分析的基础上,综合确定了凝灰岩带泥化参数为内摩擦角24°、粘聚力50k Pa、泥化率60%;凝灰岩层内摩擦角29°、粘聚力90k Pa,在工程经验取值的范围之内;玄武岩内摩擦角36°、粘聚力600k Pa,在水利水电工程岩石力学手册推荐值范围之内;IV陡倾结构面强度参数内摩擦角35°、粘聚力150k Pa,符合一般结构面试验经验值。经计算复核,埋深浅、局部分布的T3凝灰岩主滑模式下安全系数为1.65~2.15;埋深大、分布连续的T2凝灰岩主滑的安全系数为1.32~1.67,T2主控下的安全系数要比局部分布T3情况小约0.33~0.48,平均降低约21%。虽然T2主控下的最不利工况安全系数也大于1.30,但是在实际施工和运行过程中仍需关注T2凝灰岩层可能形成的深大滑移问题。类似和临近工程可参考该凝灰岩强度参数或强度参数综合识别方法。     

分子模拟在非常规油气开发中的应用

非常规油气资源储量丰富,具有广阔的开发前景.我国低渗/特低渗油藏以陆相沉积为主,储层特征差异性大,孔喉细小,孔隙度低,可动原油气大量储集于亚微米孔隙中,开采难度大,采收率低.岩油界面微观力学作用和限域传质力学机理是其中的关键科学问题.近年来,分子模拟技术在非常规油气开采的研究领域已经成为一种重要的研究手段.本文介绍了分子模拟的基本原理,然后给出了非常规油气藏储层流体和基质的组成及分子模型,并探讨了分子模拟在非常规油气开发中的应用及相关研究进展.     

Acid Leakoff Mechanism in Acid Fracturing of Naturally Fractured Carbonate Oil Reservoirs

In acid fracturing, excessive acid leakoff is thought to be the main reason that limits fracture propagation and live acid penetration distance. Since most carbonates are naturally fractured, we developed a new model in this paper to simulate acid leakoff into a naturally fractured carbonate oil reservoir during acid fracturing. Our model incorporates the acid-rock reaction, fracture width variation due to rock dissolution on the fractured surfaces, and fluid flow in naturally fractured carbonate oil reservoirs. Given the information of the reservoir, injected acid, and pressure in the hydraulic facture and the reservoir, the model predicts acid leakoff with time. In this study, we found that acid leakoff mechanism in naturally fractured carbonates is much different from that in reservoirs without natural fractures. Widened natural fractures by acid-rock reaction act as high-conductivity conduits allowing leakoff acid to penetrate deeper into the formation, resulting in serious leakoff. Wide natural fractures have a dominant effect on acid leakoff compared to micro-fractures and matrix.     

Flow Modeling of Well Test Analysis for Porous–Vuggy Carbonate Reservoirs

Porous–vuggy carbonate reservoirs consist of both matrix and vug systems. This paper represents the first study of flow issues within a porous–vuggy carbonate reservoir that does not introduce a fracture system. The physical properties of matrix and vug systems are quite different in that vugs are dispersed throughout a reservoir. Assuming spherical vugs, symmetrically distributed pressure, centrifugal flow of fluids and considering media that is directly connected with wellbore as the matrix system, we established and solved a model of well testing and rate decline analysis for porous–vuggy carbonate reservoirs, which consists of a dual porosity flow behavior. Standard log–log type curves are drawn up by numerical simulation and the characteristics of type curves are analyzed thoroughly. Numerical simulations showed that concave type curves are dominated by interporosity flow factor, external boundary conditions, and are the typical response of porous–vuggy carbonate reservoirs. Field data interpretation from Tahe oilfield of China were successfully made and some useful reservoir parameters (e.g., permeability and interporosity flow factor) are obtained from well test interpretation.     

Analytical Modeling of Flow Behavior for Wormholes in Naturally Fractured–Vuggy Porous Media

Acidizing technology has been widely applied when developing naturally fractured–vuggy reservoirs. So testing and evaluating acidizing wells’ pressure behavior become necessary for further improving the wells’ performance. Analyzing transient pressure data can estimate some key reservoir parameters. Generally speaking, carbonate minerals are usually composed of dolomite and calcite which are easy to be dissolved by hydrochloric acid which is often used to react with the rock to create a high conductivity channel, namely wormhole. Pressure transient behavior in fractured–vuggy reservoirs has been studied for many years; however, the models of acidizing wells with wormholes were not reported in previous studies. This article presented an analytical model for wormholes in naturally fractured–vuggy carbonate reservoirs, and wormholes solutions were obtained through point sink integral method. The results were validated accurately by comparing with previous results and numerical simulation. Then in this paper, type curves were established to recognize the flow characteristics, and flow was divided into six flow regimes comprehensively. The calculative results showed that the characteristics of type curves were influenced by inter-porosity flow factor, wormhole number, fluids capacitance coefficient. We also showed that the pressure behavior was affected by the angles between wormholes, and the pressure depletion increased as the angle decreased, because the wormholes were closer, their interaction became stronger. At the end, a reservoir example was showed to demonstrate the methodology of new type curve analysis.     

基于XFEM-MBEM的嵌入式离散裂缝模型流固耦合数值模拟方法

离散缝网的表征与模拟是目前国内外研究的热点. 在非常规油气开发过程中, 由于地应力场的存在会对裂缝的流动属性产生显著影响, 若将裂缝视为静态对象, 与矿场数据会出现极大偏差, 因此要基于动态裂缝做更深入的研究. 本文针对致密油藏应力场?渗流场耦合力学问题, 提出了一种高效的混合数值离散化方法, 其中采用扩展有限元法 (XFEM) 求解岩石的弹性形变, 采用了混合边界元法 (MBEM) 精确计算基岩与裂缝间的非稳态窜流, 这两种数值格式是完全耦合的, 并对整体计算格式的时间项进行了全隐式求解, 可准确表征致密油藏开采过程中的裂缝变形及流体流动机理. 此外, 本文采用了嵌入式离散裂缝前处理算法显式表征大尺度水力压裂缝, 并考虑了支撑剂的作用; 采用了双孔有效应力原理和双重介质隐式裂缝表征方法, 可捕捉基质与小尺度天然裂缝的动态信息; 由此, 本文所提出的混合模型综合表征了基质?天然裂缝?水力压裂缝共同组成的致密油藏复杂渗流环境, 并通过几个实例论证了模型的准确性, 研究表明: 对致密油藏压裂水平井进行产能评价时, 应力场所引起渗流参数的改变及裂缝开度降低的影响不可忽略. 本文研究可为非常规油气资源的开发提供理论指导.      

基于数字摄影测量的结构面信息表征及应用研究

结构面的存在直接影响岩体的强度、变形特征及破坏机制。为探究夏甸金矿矿区结构面的存在对岩体力学参数和岩体质量情况的影响，利用数字摄影测量技术快速获取研究区域结构面几何信息，并通过三维建模、激光扫描技术以及VR-Platform虚拟现实平台实现了结构面二维几何信息的真三维表征。基于此，利用广义Hoek-Brown强度准则分析计算了各研究区域的岩体力学参数，并根据岩体基本质量指标BQ值进行了岩体质量分级评价。结果表明，夏甸金矿岩体力学参数指标和岩体质量指标均较好；矿体下盘距离矿体方向越远，岩体完整性系数和岩体力学参数越低；利用广义Hoek-Brown强度准则计算所得岩体力学参数除内摩擦角外数值均与前期试验所得数据基本一致，计算所得岩体内摩擦角度偏大的原因，可能是其未考虑结构面自身的颗粒直径和形态等对岩体内摩擦角的影响。     

基于岩石蠕变及D-P准则的深部巷道围岩弹塑性分析

为深入分析岩石蠕变及中间主应力对深部高应力软岩巷道围岩稳定性及分区变形特性的影响机理, 通过剖析长期荷载作用下深部巷道围岩的变形特征, 揭示了基于围岩蠕变特性的深部巷道围岩四分区变形机理, 并依据Druck-Prager准则及关联性流动法则考虑了中间主应力、岩体扩容及应变软化的影响, 推导出巷道围岩各分区应力、变形和半径解析解. 结合工程算例, 通过对现场监测数据以及不同力学模型的计算结果进行对比分析, 论证了本文理论的科学性, 并揭示了不同强度准则和围岩参数对深部巷道围岩各分区形态的影响规律. 结果表明, 忽视深部巷道围岩的蠕变特性将导致初始黏聚力的取值大于实际值, 从而造成围岩的理论承载能力大于实际承载能力; 在[0, 0.7]区间内提高中间主应力系数可有效控制围岩塑性区及破碎区范围的扩展; 当初始黏聚力及内摩擦角减小时, 中间主应力系数对围岩塑性区、破碎区范围及巷道变形的影响力显著提高. 研究成果可为地下工程支护设计及围岩稳定性评估提供理论参考.      

二氧化碳在砂岩透镜体中充注封存的盖层岩石抗断裂性能分析

把CO2这一主要的温室气体注入到地下深处具有适当封闭条件的地层中进行封存和隔离,已被公认为是有效减少CO2排放量的一种比较安全的技术途径。砂岩透镜体油气藏具有良好的圈闭构造和储层物性,油气濒临枯竭的砂岩透镜体是较理想的CO2地质封存箱。基于币形裂纹模型和水力致裂原理,将纵向厚度和横向展布长度均远小于盖层岩石尺度的水平产状砂岩透镜体简化为盖层岩石中的I型币形裂纹,从岩石断裂力学角度分析封存箱盖层岩石的抗断裂性能。采用叠加原理给出了盖层岩石币形裂纹尖端（砂岩透镜体尖灭部位）应力强度因子的计算公式,在此基础上提出了断裂力学判别准则（K=KIC）和临界有效压应力判别准则（P=PC）,从岩石断裂力学角度为砂岩透镜体封存箱盖层岩石抗断裂性能分析和评价提供了一种新的研究思路。     

319国道重庆段裸洞隧道灰岩力学特性试验研究

裸洞隧道在长期使用过程中,围岩遭受风化、地下水及区域地应力等不利因素的作用,导致岩体性质恶化,影响隧道的正常使用。在调研统计的基础上,选取319国道重庆至武隆段大溪河裸洞隧道为研究对象,对其表层及深层灰岩进行了室内单轴压缩、三轴压缩力学特性试验并对单轴压缩作用下的声发射特性进行了研究。结果表明:单轴压缩作用下两种岩样均出现脆性破坏特征,三轴压缩条件下,两种岩样呈压缩为主的延性扩容破坏,声发射可较好地表征岩石单轴压缩破坏过程;获得了两种岩样的基本力学参数,深层岩样强度为表层岩样强度的2.13倍,弹性模量、泊松比、内聚力、内摩擦角分别是表层岩样的1.04倍,0.63倍,1.45倍,1.37倍。研究结果可为该地区裸洞隧道的长期安全性与稳定性提供有益的理论支撑。     

人工冻结深部软岩单轴力学性能试验及蠕变模型1）

为了了解深部软岩在冻结条件下的单轴力学性能,以东北地区的原状泥砂岩为试验对象,利用自行研制的WDT-100型人工冻土试验仪器,对其进行不同温度下的人工冻土单轴抗压强度试验和单轴蠕变试验,得到泥砂岩单轴压缩应力-应变关系曲线,各温度下试样的单轴抗压强度以及蠕变曲线.单轴压缩试验结果表明：试样在给定温度和加载速率条件下,单轴压缩应力-应变关系曲线都有较为明显的屈服点,并且都在屈服点后,强度有所提高,出现硬化现象.单轴蠕变试验结果表明：单轴压缩蠕变曲线有非线性特征,单轴压缩蠕变的等时应力-应变曲线向应变轴靠拢;单轴压缩时蠕变模量随时间的增长而降低.最后采用遗传算法优化模型参数,得出泥砂岩蠕变经验方程.与试验结果对比,发现拟合情况较好.     

开挖卸荷条件下深部圆形洞室各向同性围岩分区破裂化机理

论文提出了一种静水压力条件下动态开挖卸荷对深埋圆形洞室各向同性围岩分区破裂化影响的力学模型,获得了动态开挖卸荷条件下深埋圆形洞室围岩的弹性应力场,确定了卸荷速率和岩体动态力学参数对深埋洞室各向同性围岩分区破裂化现象的影响.利用动态Hoek-Brown准则,获得了围岩破裂区和非破裂区的数量和宽度.数值计算结果表明卸荷速率和岩体动态力学参数对破裂区与非破裂区的数量和宽度有较大影响.