考虑介质微观非均质特征的乳状液渗流模型

油水乳化渗流是三元复合驱和热力采油过程中常见的现象,地层介质的微观孔隙结构特征对乳状液流动有着重要影响.现有描述乳状液渗流的理论模型都属于确定性方法,只能反映出介质孔隙结构的体积平均效果.当介质内部微观非均质性相比其尺寸不能被忽略时,采用确定性方法描述会与实验结果存在偏差.基于连续时间随机游走理论建立了描述乳状液渗流的随机理论模型.该模型引入反映液滴微观运动特征的跃迁时间和跃迁位移两个概率分布函数来反映多孔介质微观非均质特征.研究结果表明该模型能很好地刻画实验曲线中出现的与介质尺度相关的拖尾现象,可作为更一般的过滤模型.     

基于规则网络的碳酸盐岩多尺度网络模型构建方法研究

碳酸盐岩油藏非均质性强,孔隙大小变化可达好几个数量级,描述碳酸盐岩油藏多尺度孔隙特征具有重要意义.本文首先基于三维规则网络模型建立了不同物理尺寸的溶洞网络、大孔隙网络和微孔隙网络;然后提出一种耦合算法,以溶洞网络为基础,通过添加适当比例的大孔隙和微孔隙,构建出碳酸盐岩多尺度网络模型;最后对比分析了各网络模型的几何性质、拓扑性质和绝对渗透率.结果表明,碳酸盐岩多尺度网络模型能够同时描述不同尺度孔隙的几何和拓扑特征;且相比各单一尺度的孔隙网络模型,多尺度网络模型有着较高的绝对渗透率,这是由于各尺度孔隙之间的相互连通极大地提高了网络的整体连通性和流动能力,为碳酸盐岩油藏微观渗流模拟提供了重要的研究平台.     

多重介质渗流研究进展

从已发现的油气田来看,裂缝性油气田通常都是大油气田,其储量很大,单井产量每天数百吨到上万吨不等,经济价值极高。油、气、水在裂缝性孔隙介质中的流动特征与一般孔隙介质的特征有本质上的差异:各重介质问液量交换性,多重介质的非均质性以及各向异性等.因此,研究多重介质渗流有较大的理论和实际意义。      

孔隙网络模型在渗流力学研究中的应用

孔隙网络模型广泛地用于同多孔介质有关的微观模拟中. 本文简介了利用网络模型研究储层渗流规律的基本思想、主要步骤以及网络模型同其它微观模型相比较的优缺点, 总结了在微观渗流研究中使用的两大类网络模型(准静态网络模型和动态网络模型)的特征及其适用范围, 综述了目前网络模型在研究渗流中的应用现状以及国内外研究比较活跃的几个方面,最后分析了网络模型今后的发展趋向.      

混凝土梁三点弯曲断裂模式转变的近场动力学模拟

混凝土结构的宏观损伤开裂与其非均质微观结构紧密相关。底部带切口的混凝土梁在进行三点弯曲破坏时,随着切口的位置由梁中向梁边转移,裂纹由从切口处萌生并生长转变为从梁的中部萌生。本文采用半均质化近场动力学（IH-PD）模型和全均质化近场动力学（FH-PD）模型,分别对混凝土梁三点弯断裂问题进行模拟研究。IH-PD模型根据混凝土中骨料体积分数随机生成不同键的组合方式,将微观尺度的非均质性引入模型,无需详细描绘骨料形状和分布即可考虑混凝土非均质性。本文将IH-PD与FH-PD模型得到的断裂模式随切口位置的变化关系,与实验结果对比,分析微观结构对混凝土梁开裂的影响;基于非均质材料特征尺寸与IH-PD模型网格参数的相关性,模拟骨料大小对混凝土梁断裂模式的影响;另外,通过在IH-PD模型中设置预损伤的方式引入随机分布的孔隙,探讨孔隙率对混凝土断裂模式的影响。     

渗流气体滑脱现象与渗透率变化的关系

气体在致密多孔介质中低速渗流时存在着因气体分子碰撞岩壁而引起的滑脱现象，它由介质的孔隙结构和气体分子的平均自由程共同决定。该现象使气测渗透率大于孔隙介质的绝对渗透率。介质中气体的低速渗流为黏滞流与滑脱流组成，各自所占比例与气体分子按自由程的分布有关。理论计算得到了低速气体渗流的气测渗透率Kg与绝对渗透率K0比值的关系式，实验结果与理论分析吻合。     

页岩气在超低渗介质中的渗流行为

基于页岩气藏复杂孔隙结构和页岩气在纳米孔隙表面的滑脱和吸附-解吸附等现象,通过引入表观渗透率,修正Darcy渗流模型,得到了页岩气渗流本构方程. 将计算结果与Darcy模型计算结果进行了比照,结果表明,在产量定常情形下,基于非Darcy渗流模型得到的井底压力高于Darcy流结果; Darcy流模型得到的压力衰减速度较快,改进后的渗流模型更能准确描述与表征页岩气渗流过程;研究结果可为页岩气藏的经营与管理提供基础参数.      

海陆过渡相页岩气藏不稳定渗流数学模型

海陆过渡相页岩常与煤层和砂岩呈互层状产出, 储层连续性较差、横向变化快、非均质性强, 水力压裂技术是其获得经济产量的关键手段. 然而, 目前缺乏有效的海陆过渡相页岩气藏不稳定渗流数学模型, 对其渗流特征分析及储层参数评价不利. 针对这一问题, 首先建立海陆过渡相页岩气藏压裂直井渗流数学模型, 其次采用径向复合模型来反映强非均质性, 采用Langmuir等温吸附方程来描述气体的解吸和吸附, 分别采用双重孔隙模型和边界元模型模拟天然裂缝和水力裂缝, 建立并求解径向非均质的页岩气藏压裂直井不稳定渗流数学模型, 分析海陆过渡相页岩气藏不稳定渗流特征, 并进行数值模拟验证和模型分析应用. 分析结果表明, 海陆过渡相页岩气藏不稳定渗流特征包括流动早期阶段、双线性流、线性流、内区径向流、页岩气解吸、内外过渡段、外区径向流及边界控制阶段. 将本模型应用在海陆过渡相页岩气试井过程中, 实际资料拟合效果较好, 其研究成果可为同类页岩气藏的压裂评价提供一些理论支撑, 具有较好应用前景.      

脆性岩石破裂过程损伤与渗流耦合数值模型研究

大量的实验结果表明，脆性岩石的渗透性不是一个常量，而是应力和应力诱发损伤破裂的函数．建立了一个描述非均匀岩石渗流-应力-损伤耦合数学模型(FSD Model)，开发出岩石破裂过程渗流-应力-损伤耦合分析计算系统(F-RFPA^2D)．在该系统中，单元的力学、水力学性质根据统计分布而变化，以体现材料的随机不均质性，材料在开裂破坏过程中流体压力传递通过单元渗流-损伤耦合迭代来实现．该系统能够对岩石试件在孔隙水压力和双轴荷载作用下裂纹的萌生、扩展过程中渗透率演化规律及其渗流-应力耦合机制进行模拟分析．最后，给出两个算例：算例1模拟载荷作用下岩石应力应变-渗透率全过程．模拟结果表明，非均匀性对岩石的应力峰值强度、峰值前后其渗透性演化规律及其破裂机制影响十分明显，模拟结果和实验结果较为一致；算例2模拟孔隙水压力作用下岩石拉伸断裂过程，通过和物理实验对比验证，验证了模型计算结果的可靠性。     

不同围压作用下非均匀岩石水压致裂过程的数值模拟

从岩石细观非均匀性的特点出发，提出一个描述非均匀材料渗流和破裂相互作用的数值模型。在这个数值模型中，单元的力学、水力学性质根据统计分布而变化，以体现材料的随机不均质性，材料在开裂破坏过程中流体压力传递通过单元渗流，损伤耦合迭代来实现。算例表明，该模型能较好地模拟出岩石类材料在水力压裂作用下，微结构非均匀分布和不同围压比对破裂模式、失稳压力的影响，非均匀性导致试件的开裂压力、失稳压力明显不同，裂纹扩展路径不规则发展，模拟结果和实验结果较为一致。     

锦屏一级水电站地下厂区破坏成因分析

锦屏一级水电站地下厂区位于雅砻江右岸一套大理岩地层中,埋深350m,地下厂区洞室在开挖过程中,围岩出现各种变形破坏现象。通过分析破坏现象,总结了破坏现象的规律性：主厂房及主变室破坏现象主要集中在下游边墙顶拱,破坏方式以岩体开裂为主; 母线洞及压力管道下平段破坏现象主要集中在顶拱外侧,破坏方式以弯折破坏为主。最后利用数值模拟结果解释了径向应力大于岩体强度是引起锚索超限的原因。     

汶川大地震后水坝建设中若干问题的思考

应用参数摄动法对可压缩N-S方程进行渐近展开,并取其零阶近似对高压下微管道液体流动 特性进行了分析．对任意截面形状和面积的微管道,在等温流动假设下将其截面形状、滑移 长度等对解的贡献转化为求解该截面的格林函数,并给出等截面圆形微管道流动的零阶近似 解．以此分析可压缩性、黏性以及壁面滑移等因素对高压下液体微管道流动特性的影响,进 一步揭示了高压驱动下液体微管道流动偏离经典Hagen-Poiseuille(HP)理论的原因．     

四川泸定磨西台地第四纪冰水台地边坡地质灾害易发性研究

河流冲刷对第四纪软弱地层岸坡的地质灾害起到控制性作用。本文以磨西台地冰水堆积物岸坡为例,充分考虑河流与岸坡地质灾害之间的内在联系,基于对实测地质灾害的分析结果,采用Logistic监督分类方法,评价各环境因子（尤其是河流冲蚀特性因子）在河流岸坡稳定性评价中的作用,进而采用高相关的环境因子进行地质灾害易发性预测。研究表明,河流弯曲度、河流流量、岸坡坡度、河床宽度、岸坡宽度与地质灾害的发育相关性强,采用这些参数和Logistic模型预测的地质灾害易发性能很好的反映现状和未来地质灾害的发育特点,进而提出了磨西台地岸坡地质灾害防治的初步建议。     

瀑布沟水电站库首右岸深部裂缝成因分析

水电工程常常建在高山峡谷地带,其天然岸坡通常由坡面向内有一个强卸荷带和弱卸荷带以及相应的强风化和弱风化带,内侧则为完整新鲜的岩石。对涉及的工程岸坡在正常卸荷带以内发育的一系列张性破裂或破裂带,称之为“深拉裂缝”。瀑布沟水电站库首右岸存在两个拉裂变形体,通过对其岸坡深部拉裂缝空间发育分布、变形特征的考察,综合分析造成深部裂缝发育规律与变形特征的因素。在此基础上提出,库首右岸深部拉裂缝是岸坡快速卸荷条件下浅表生改造的产物,其形成时期相当于河谷由宽谷深切为峡谷这一转换时期。     

大瑞铁路澜沧江大桥工程边坡岩体结构特征研究

岩体中结构面的展布及其组合特征决定了岩体的工程地质性质和力学性状,影响着岩体的破坏方式。大瑞铁路澜沧江大桥所处的澜沧江构造带,具有高地震烈度、活跃的新构造运动、活跃的外动力地质条件、活跃的岸坡浅表改造过程等地质特征。工程区两岸节理裂隙等弱面发育,岸坡浅表改造过程强烈,导致岩体结构复杂。通过野外现场调查,对两岸边坡岩体结构面数据进行统计分析,采用^x2检验与K-S检验法对结构面进行概率密度拟合,分析岩体结构特征。结果表明：(1)右岸的优势结构面有三组; (2)左岸的优势结构面有四组; (3)各组优势结构面产出状态数据均服从正态分布。在统计分析结果的基础上,建立两岸岩体结构的统计模型,并对边坡破坏模式进行定性分析,认为：(1)右岸边坡变形破坏模式为滑移 压致拉裂; (2)左岸结构面贯通性不好,主要的破坏模式是局部块体滑移。
     

黄河上游戈龙布滑坡高速下滑成因机制及堵江分析

黄河上游戈龙布滑坡位于积石峡水电站库区内,经过现场调查可知,戈龙布滑坡为一大型滑坡,整个滑坡体积达到5040万m3,一共分为4个区。研究表明,戈龙布滑坡成因机制为滑移－拉裂型,滑坡在黄河右岸经过高速远程下滑过后撞击上游山体,部分滑体越过黄河堆积于左岸并堵塞黄河达数百年,堰塞湖内有厚度达30余米的纹泥沉积,初步分析滑坡坝的高度在100m左右,库水深在70m左右,滑坡堵江时代在Q4时期。滑坡坝溃决后导致黄河改道,残留的滑体分别位移黄河的左岸和右岸。     

一种高边坡岩体卸荷分带方法的探讨

岩体卸荷是水利水电工程建设中常见的工程地质问题之一。本文介绍了高边坡岩体卸荷的背景条件、作用机理、卸荷岩体特征以及常见的卸荷带划分方法。针对河谷斜坡岩体从表部到深部 ,岩体裂隙发育程度、张开程度逐渐降低的特征 ,以黄河上游某水电站高边坡印支期花岗岩体的现场调查资料为依托 ,应用岩体裂隙发育程度、张开程度量化指标对斜坡岩体卸荷进行分带 ,较好地将岩体卸荷分带与岩体结构、岩体工程特性紧密结合起来 ,为以后其它边坡岩体的工程实践提供了借鉴.     

汶川地震极重灾区地质背景及次生斜坡灾害空间发育规律

5·12汶川大地震造成大量的次生斜坡灾害,本次研究区域为汶川大地震的11个重灾区,包括汶川、北川、青川、安县、平武、茂县、江油、彭州、什邡、绵竹、理县等市县。通过对重灾区航片、卫片、雷达图像的解译研究发现,重灾区次生斜坡灾害的主要灾种表现为崩塌、滑坡以及崩塌、滑坡高速运动解体形成的碎屑流（个别地方由于水的参与表现为泥石流）以及它们堵江形成的堰塞湖。研究发现地震次生斜坡灾害的发育具有明显的丛集性规律。从区域上看,次生斜坡灾害明显呈带状,沿龙门山断裂带展布,并主要受北川—映秀断裂控制。各灾种的发育在不同地段发育的规模、频率差别较大。以灾害分布面积来排序,汶川县灾害面积最大,为131.55km2,其次为北川县,为45.57km2,其余9个县（市）灾害面积相差不大,均介于6~17km2,其中理县灾害面积最小,为6.25km2。各灾种的发育在不同地段发育的规模、频率差别较大。青川县、平武县灾种主要为滑坡,汶川县、茂县、安县、理县灾种主要表现崩塌转化的碎屑流,北川的主要灾种则为碎屑流,其次为滑坡,什邡、彭州、绵竹、江油等地主要灾种为崩塌。
灾种发育的这种地域性差别主要受控于地层岩性,除此而外,还与构造特征、地形地貌等因素紧密相关。研究表明：岩性对灾害种类的展布有决定性控制作用。统计发现,岩性越坚硬,崩塌、碎屑流发育率越高,滑坡则在软岩地区、较软岩地区和较坚硬区发育率最高,泥石流则在软岩地区最为发育。地形地貌对次生斜坡灾害的发育有重要影响,统计表明,崩塌、碎屑流以及泥石流在1200～2000m坡段范围内发育率最高,其次为800～1200m坡段;而滑坡则在800～1200m坡段范围发育率最高。对坡度而言,除11°~20°坡度范围外,崩塌和碎屑流的发育率总体具有随坡度增高而增大的特点;而滑坡和泥石流的发育率呈现典型的单峰特征,在1°~20°范围内发育率最大。坡向对地震次生斜坡灾害的发育影响不明显。
地震次生斜坡灾害的发育规律表明,地震斜坡灾害的发生主要受控于活动构造本身,并沿活动构造呈带状展布,同时受场地条件如岩性、地形地貌等因素的强烈控制。     

汶川大地震山地灾害发育的控制因素分析

汶川大地震震后地质灾害调查结果表明,龙门山山地灾害无论是空间分布还是在剖面分布上均具有明显的不均一性：平面上集中在发震断裂及非发震断裂及其附近,剖面上集中在坡顶、坡折及孤立的突出地形上;运动形式具有加速起飞-滑翔-俯冲着陆三阶段。主要受控因素有发震断裂、高烈度、地震持续时间、斜坡地形效应、岩性及结构等。     

秦岭公路隧道2号竖井地应力与岩爆分析

利用水平圆柱形激波管对激波驱动的可压缩性气固两相流进行了试验研究. 利用压电式压力传感器、电荷放大器、示波器及计算机组成的压力信号测试系统, 对激波 与颗粒作用前后的气相参数进行测量及分析. 试验中测得了激波在管中的传播速 度, 波后气流的压力, 反射激波、透射激波的压力和速度等. 分别考察颗粒、装载 比、驱动气源以及入射激波马赫数等因素的差异对气相参数的影响. 试验结果表明: 激波与颗粒群相互作用时, 会产生反射激波和透射激波, 其强度与驱动气源、颗粒大小、颗粒装载比等参数有关; 激波衰减率随着装载比、马赫数的增大而减小. 研究指出, 在颗粒群被激波加速的初始阶段, 颗粒间的弹性碰撞起着重要的作用.